Discussion :

[twix24]

Bonjours,
je veux réaliser une opération d'extraction de contours à partir d'images IRM. aprés le passage par le fitre de canny, je veux utiliser le résultat comme un masque que je vais appliquer sur mon images original. mais je rencontre des problèmes au niveau de la multiplication de deux images.

Code :

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

1
2
Error in ==> canny at 22
Ic1= IM.* BW;

voila 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
IM = dicomread('brain_010.dcm');
[m,n]= size(IM);
BW=uint8(zeros(m,n));
BW = edge(IM,'canny');
for(i=1:m)
    for(j=1:n)
        if(BW(i,j)== 0)
            BW(i,j)=1;
        else
            BW(i,j)=0;
        end
    end
end
figure
subplot(1,2,1)
imshow(IM,'DisplayRange',[])
axis image
subplot(1,2,2)
imshow(BW)
axis image
zeros(m,n);
Ic1= IM.* BW;

merci d'avance

[Jerome Briot]

Pourrait-on avoir le message d'erreur en entier ?

Note, à la place de :

Code :

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

1
2
3
4
5
6
7
8
9
for(i=1:m)
    for(j=1:n)
        if(BW(i,j)== 0)
            BW(i,j)=1;
        else
            BW(i,j)=0;
        end
    end
end

Il est plus rapide de faire :

Code :

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

BW=~BW;

[twix24]

le message d'erreurs

Code :

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

1
2
3
4
5
??? Error using ==> times
Integers can only be combined with integers of the same class, or scalar doubles.
 
Error in ==> canny at 22
Ic1= IM.* BW;

[twix24]

merci pour la réponse.

le problème de la multiplication de l'erreur au niveau de la multiplication est résolu par le code suivant

Code :

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

1
2
3
4
5
6
7
IM = dicomread('brain_010.dcm');%lire l'image IRM
[m,n]= size(IM); %taille de l,image IRM
IM=double(IM);
maximum = max(IM(:));
minimum = min(IM(:));
B=zeros(m,n)
B=uint8(255*(IM-minimum)/(maximum-minimum));

mais j'arrive pas encore a obtenir le résultat volu
mon fichier original et mon masque :

le résultat obtenu


le but de mon programme est de parcourir les deux image (L'originale est le masque binaire) et de garder les pixel qui correspont(en terme de position) a au pixel ou le contour n'est pas nul.
de cette façon je peut obtenir une reconstitution de contours avec des diffirent niveau de gris.

merci.

[Jerome Briot]

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
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
IM=rand(6);
BW=rand(6)>.5;
 
figure
colormap(gray)
 
subplot(2,2,1)
imagesc(IM)
caxis([0 1])
title('IM')
axis image
 
 
subplot(2,2,2)
imagesc(BW)
caxis([0 1])
title('BW')
axis image
 
subplot(2,2,3)
imagesc(~BW)
caxis([0 1])
title('~BW')
axis image
 
subplot(2,2,4)
IM(~BW)=0;
imagesc(IM)
caxis([0 1])
title('IM(~BW)=0')
axis image