Discussion :

[c2.pierre]

Bonjour;

Je souhaite tracer une ROI polygonale (main libre) sur une image en utilisant la fonction (roipoly)
j'ai fait une tentative en pensant à l'addition des deux images, mais sans succès.

Merci d'avance de vos conseils et aides.

voici mon code

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
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
 
repertoire = input('Répertoire: ','s');
 
rep = repertoire;
ext = '*.dcm';
 
chemin = fullfile(rep ,ext);
 
d = dir(chemin);
%%  Tracer la ROI de référence plan 20 
IM1 = I(:,:, 20);
figure (1);
imshow(IM1, []);
[x, y] = getpts;
 
x = round(x); y = round(y);
 
for n=1:numel(d)
 
   file_TEP =  fullfile(rep,d(n).name);
   info(n) = dicominfo(file_TEP);
   I(:,:,n) = dicomread(info(n));
   figure (1);
   imshow(I(:,:,n), []);
   title(['Image   ' num2str(n) ' de la séquence  ']);
 
   %% Recopie la ROI sur la séquence 
    Im = I(:,:,n) ;
    J(:,:,n) = roipoly(Im,x,y);
    Im_Roi  = J + I ; %== Addition des deux images 
   figure (2);
   imshow(Im_Roi, []);
   title(['ROI de lImage   ' num2str(n) ' de la séquence  ']);
 
   end

[Jean Dumoncel]

Bonjour,

La fonction roipoly renvoie un masque binaire, donc, par exemple, pour afficher uniquement la roi, tu peux faire :

Code :

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

1
2
J(:,:,n) = roipoly(Im,x,y);
Im_Roi  = J(:,:,n).*Im;