Image de données codé en 16 bit

**kevinl** · 16/04/2007, 11h57

Bonjour,
Je possède une matrice de données provenant d’un fichier codé en 16 bit (nombre compris entre 0 et 65535). Je voudrais créer une image avec. J’ai essayé avec imagesc mais les résultats que j’obtiens sont aberrants. J’aimerais donc savoir s’il y a une manière particulière sous Matlab de représenté des données qui ont été codés en 16bit.
Merci

Matlab 7.2.0.232 (R2006a)

**amarion** · 16/04/2007, 12h02

Envoyé par kevinl

Bonjour,
Je possède une matrice de données provenant d’un fichier codé en 16 bit (nombre compris entre 0 et 65535). Je voudrais créer une image avec. J’ai essayé avec imagesc mais les résultats que j’obtiens sont aberrants. J’aimerais donc savoir s’il y a une manière particulière sous Matlab de représenté des données qui ont été codés en 16bit.
Merci

Matlab 7.2.0.232 (R2006a)

Bonjour,
essaie imshow(Mat,[])
les crochets permettent l'adaptation de dynamique

**Jerome Briot** · 16/04/2007, 12h21

Envoyé par kevinl

les résultats que j’obtiens sont aberrants.

Qu'est-ce que tu entends par "aberrants" ?

**rostomus** · 16/04/2007, 13h43

Bonjour,

Essaie de faire un cast avec uint16.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

imshow(uint16(mat))

**kevinl** · 16/04/2007, 15h33

Bonjour à tous et merci pour vos pistes,

Malheureusement la fonction imshow n’est pas une commande/fonction reconnue chez moi.

Je parle de valeur aberrante parce que l’image que je devrais obtenir est une carte du monde et j’obtient une image ne représentant pas du tout une carte.

**Jerome Briot** · 16/04/2007, 15h37

Envoyé par kevinl

Je parle de valeur aberrante parce que l’image que je devrais obtenir est une carte du monde et j’obtient une image ne représentant pas du tout une carte.

OK donc le problème ne vient pas de toi alors

Quel est le format de l'image... on parle toujours du même : Fichier binaire représentant un tableau à trois dimensions

**kevinl** · 17/04/2007, 10h17

Bonjour,

Oui on parle toujours de cette image qui provient d'un fichier binaire.
La dernière explication était intéressante mais malheureusement si je considère que chaque pixel est codé avec le triplet de couleur (24bits) mon fichier devrait faire 700Mo or il en fait 445Mo ce qui correspond exactement à un codage de 16 bit par pixel.

**Jerome Briot** · 17/04/2007, 10h24

Peux-tu me donner le lien exacte vers l'image sur le site de la NASA ?

**kevinl** · 17/04/2007, 10h35

Il s'agit de la carte :
gebco_bathy.21601x10801.bin.gz.torren

au lien suivant:
http://visibleearth.nasa.gov/view_detail.php?id=8392

**Jerome Briot** · 17/04/2007, 11h08

Bittorrent a un peu de mal ce matin avec le serveur de la NASA...

Sinon ton code doit ressembler à ceci ?:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
fid=fopen('gebco_bathy.21601x10801.bin','r');
X=fread(fid,[21601,10801],'*uint16');
fclose(fid);
 
figure
colormap(gray)
imagesc(X);

Par contre, vu la taille de l'image... tu dois frôler le OUT OF MEMORY non ?

**kevinl** · 17/04/2007, 12h05

Oui, j'avais essayé plus ou moins comme ça et je suis bien OUT OF MEMORY.
Donc ce que j'ai essayé de faire c'est prendre seulement une partie de l'image, puisqu'il y a seulement la région proche de la France qui m'intéresse.
Voila le code que j’utilise pour prendre seulement une partie de l’image (je ne suis pas sure que ce soit correcte).

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
fid=fopen('D:\Mes documents\bathymetrie\gebco_bathy.21601x10801.bin','r'); 
IMImage=fread(fid,100000000,'uint16=>uint16');
fclose(fid);
 
iIM=1;
xloc=0;
yloc=0;
colorImage = zeros(300,21601);
 
for a = 83595870:89428140
    if xloc == 21601 || xloc == 0
        xloc=0;
        yloc= yloc+1;
    end
 
    xloc= xloc+1;
    colorImage(yloc,xloc) = IMImage(a);
 
end

**Jerome Briot** · 17/04/2007, 12h11

De quelle quantité de mémoire disposes-tu ?
As-tu lu ceci : Gestion de la mémoire

**sango85** · 17/04/2007, 12h32

Est-ce que par hasard le problème ne viendrait-pas du fait que la lecture de l'image se ferai sur 8 bit, mais sur les 8 bits de points faibles (très peu représentatifs pour une image),
alors que les 8 bits de points forts afficheraient une représentation (certes moins précise que du 16 bit) mais cohérente...
(Ce n'est qu'une piste de réflexion... Peut-être totalement erronnée..)

**kevinl** · 17/04/2007, 13h40

Pour la mémoire j'ai 512Mo de RAM ce qui ne vas pas être suffisant, c'est pour cela que je voulais prendre seulement une partie de la carte.

**Jerome Briot** · 17/04/2007, 13h41

Essaie ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
W=21601
H=10801
 
r=4000;
c=8000;
 
w=6000;
h=6000;
 
fid=fopen('gebco_bathy.21601x10801.bin','r');
    fseek(fid,(W*r+c)*2,'bof');
    X=fread(fid,[h,w],sprintf('%d*uint16=>uint16',w),2*(W-w));
fclose(fid);
 
figure
imagesc(X.');axis equal

r et c correspondent à la position du pixel en haut à gauche de la partie de l'image que tu veux lire. w et h sont respectivement, le nombre de colonnes (largeur) et le nombre de ligne (hauteur) de cette même partie.

Tout ça bien sur... si je ne me suis pas trompé

**kevinl** · 17/04/2007, 14h07

Je viens d'essayer le bout de code que tu viens de me donner et ça ne donne toujours rien (voir image en pièce jointe).

D'un autre côté je commence à me demander si je fichier que je possède est bien complet. Puisque d'après les explications du site le la NASA il devrait être codé en 3*8bits et le fichier complet devrait donc être plus important.

En me baladant sur internet j'ai trouvé une autre carte qui ressemble à celle qu'on est en train d'essayer de visualiser mais elle est en format ecw et je n'ai rien trouvé sur ce type de format dans l'aide de Matlab. Est ce que ce format est supporté par Matlab

**Jerome Briot** · 17/04/2007, 14h20

Revoici le code un peu modifié :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
W=21601
H=10801
 
r=1000;
c=10000;
 
w=4000;
h=3000;
 
fid=fopen('gebco_bathy.21601x10801.bin','r');
    fseek(fid,(W*r+c)*2,'bof');
    X=fread(fid,[w,h],sprintf('%d*uint16=>uint16',w),2*(W-w));
fclose(fid);
 
figure
colormap(gray)
imagesc(X.');axis equal

**kevinl** · 17/04/2007, 15h47

Je viens d'essayé ce code. La carte est bien en noire et blanc maintenant mais par contre j'ai beau me déplacé sur la carte il n'y a rien qui ressemble a des continents.

**Jerome Briot** · 18/04/2007, 11h15

sango85 avait presque raison... il faut lire le fichier en big-endian :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 W=21601
H=10801

r=1000;
c=10000;

w=4000;
h=3000;

fid=fopen('gebco_bathy.21601x10801.bin','r','b');
    fseek(fid,(W*r+c)*2,'bof');
    X=fread(fid,[w,h],sprintf('%d*uint16=>uint16',w),2*(W-w));
fclose(fid);

figure
colormap(gray)
imagesc(X.');axis image

Je l'ai testé avec l'image, ça fonctionne parfaitement pour moi

**kevinl** · 18/04/2007, 12h31

Merci beaucoup pour ton aide.

Ca marche aussi très bien chez moi. Je vais enfin pouvoir commencer à utiliser les données.

Image de données codé en 16 bit

Images

Discussions similaires

Partager

Partager