Discussion :

[cedji]

Bonjour,
Je dispose d'un fichier raw contenant des images et j'aimerai savoir s'il y avait une méthode simple pour:
- déterminer le nombre d'images du fichier
- pour pouvoir se déplacer librement dedans

Voila le code que j'ai fait mais il n'est pas tres pratique :

Code :

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

1
2
3
4
5
6
7
8
9
10
fichier1 = 'HighEnergyCTHeadCoil.raw'
f1 = fopen(fichier1,'r');
im1 = fread(f1,[512 512],'int16','ieee-le');
 
c=0;
 
while (feof(f1)==0)
    im1 = fread(f1,[512 512],'int16','ieee-le');
    c=c+1;
end

Comme ceci jai trouve que c =305

et pour pouvoir me déplacer dans le dossier j'ai fait une boucle avec un if (par exemple pour l'image 5 jai fait : if(c==5)... end

merci de votre aide

[Jerome Briot]

- déterminer le nombre d'images du fichier

Si pour toutes les images, les valeurs des pixels sont codées sur 2 octets (16bits), prends la taille de ton fichier en octets et divise la par 524288 (512x512x2)

- pour pouvoir se déplacer librement dedans

Utilise la fonction FSEEK

[cedji]

Merci pour votre aide
Jai utilisé la fonction fseek comme vous me l'aviez suggéré :

Code :

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

1
2
3
4
f1 = fopen(fichier1,'r');
fseek(f1,268,'bof');
im1 = fread(f1,[512 512],'int16','ieee-le');
imtool(im1)

mais à l'affichage je rencontre un problème que je n'avais pas avant.
Sur l'image qui apparait on voit une partie d'une image et un bout de la suivante...

d ou vient le probleme? Est ce que le pointeur est mal initialisé??

Merci

[Jerome Briot]

A quoi correspondent ces 268 octets ?

[cedji]

ah ba c sans doute là que je me suis trompée... javais mi 268 pour afficher la 268ème image mais cela ne doit pas être ca....

[Jerome Briot]

Le second argument de FSEEK est exprimé en nombre d'octets

Donc si comme je l'ai dit, une image a pour taille 512x512x2 octets, à combien d'octet du début du fichier se trouve l'image n°5 ?

[cedji]

ah ba oui je suis bete!
pour l'image5, on doit faire :

Code :

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

fseek(f1,512*512*2*5,'bof')

Sinon pour savoir le nombres d'images dans le fichier j'ai bien fait ce que vous mavez proposé et je retombe bien sur la valeur trouvée avec la boucle while. Mais j ai recopié "à la main" la taille du fichier en regardant ses propriétés, mais y a t il une fonction sur matlab qui la donne directement? j'ai essayé avec size mais ca na pas l air detre la bonne solution...

Merci

[Jerome Briot]

Utilise la fonction DIR et regarde le champs bytes de la structure de sortie.

Sinon, tu peux toujours faire ceci :

Code :

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

1
2
3
4
5
6
fid = fopen(...,'r');
 
fseek(fid,0,1);
taille = ftell(fid);
 
fclose(fid);

[cedji]

ok ok! ca marche! Merci pour votre aide

[Jerome Briot]

ah ba oui je suis bete!
pour l'image5, on doit faire :

Code :

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

fseek(f1,512*512*2*5,'bof')

Y'avait un piège ... avec ce code, tu liras la 6ème image... puisque tu en as sautées 5

[cedji]

oups ok! C'est noté!
Merci