Discussion :

[Christinita]

Bonjour ,
je voudrais convertir un ensemble d'images dicom se trouvant dans mon répertoire courant en jpg. Voici mon code

Code :

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

1
2
3
4
5
6
7
8
9
10
p=dir('GLIOMET1Gado');
 
for i=3:size(p,1)
% lire les images output de lIRM
[X, map] = dicomread(['GLIOMET1Gado\' p(i).name]);  
 
% les covertir en Dicom et les sauvergarder dans le chemin spécifié 
imwrite(X, ['T1Gado\' [p(i).name '.jpg']]);
 
end

Les images résultantes sont toutes blanches mouchetés en noir . Or je devrais trouver des images de cerveau.

Avez vous une idée de quoi il s'agit?

Merci

[Jerome Briot]

Si la variable map n'est pas vide (vérifier avec ISEMPTY), il faut essayer de l'utiliser dans IMWRITE :

Code :

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

imwrite(X, map, ...)

[Christinita]

map est vide

[Jerome Briot]

Alors il faut faire comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
X = dicomread('test.dcm');
 
figure
imshow(X,[],'InitialMagnification',100);
 
f=getframe(gca);
img = frame2im(f);
 
imwrite(img,'test.jpg','jpg')

Il peut y avoir une légère modification de taille de l'image pendant le GETFRAME qui peut être corrigée comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
X = dicomread('temp.dcm');
 
figure
imshow(X,[],'InitialMagnification',100);
 
xl=get(gca,'xlim');
yl=get(gca,'ylim');
 
f=getframe(gca,[xl(1) yl(1) xl(2)-xl(1) yl(2)-yl(1)]);
img = frame2im(f);
 
imwrite(img,'test.jpg','jpg')

Encode une fois, c'est l'idée qui est présentée ici.
Il faudra certainement la peaufiner un peu et l'adapter au contexte (boucle FOR-END)

[Christinita]

Merci Dut. C'est résolu