Discussion :

[hulia]

Bonjour,
Je veux calculer le contour d'un objet dans une image pour ensuite calculer son aire , j'utilise le code suivant pour l'image ci-jointe , mais le résultats est tout autre (je veux calculer le contour de l'os et de la prothèse)

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
 
clear all
close all
clc
x=imread('face.png');
level=graythresh(x);
bw=im2bw(x,level);
h=fspecial('gaussian',5,5);
bw=imfilter(bw,h);
bw=imfill(~bw,'hole');
bw=~bw;
b=bwboundaries(bw);
imshow(x);
hold on
for i=1:length(b)
    contour=b{i};
    plot(contour(:,2),contour(:,1),'g','LineWidth',3);
end

[Gooby]

Bonjour,

Qu'est-ce qui ne va pas concrètement?

Je pense qu'au vue de ce que tu dis, l'erreur que tu effectues est dans le seuillage de ton image. Si tu fais un graythresh de ta figure, tu auras soit le fond de ton image confondu avec l'os, soit l'os confondu avec la protèse puisque ton image sera binaire après seuillage et que tu n'as pas une double inclusion mais des frontières communes pour chacun des 3 éléments de ton image.

[hulia]

Je cherche juste le moyen de calculer la surface de la prothèse et celles de l'os, donc je suis partie de l'idée de detecter le contour des deux objets puis de calculer leur surfaces.

[Gooby]

Si tu regardes le 2ème argument que renvoie la fonction bwboundaries, c'est directement la surface.

[B,L] = bwboundaries(bw);

[Gooby]

Cela dit, cela ne changera pas le fond de ton problème qui sera de discriminer la prothèse de l'os pendant ton seuillage. Je suppose qu'il y a plusieurs solutions facile à mettre en place pour cela.

[hulia]

oui , en cherchant j'ai tenté le seuillage , mais étant débutante en matlab je n'arrive pas à imposer un bon seuil pour différencier entre le gris et le noir.
en gros j'ai essayé ça

Code :

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

1
2
3
4
5
6
7
8
9
10
 
clear all
close all
clc
img=imread('nb.tif');
figure;imshow(img);
img=im2double(img);
figure;subplot(1,2,1);imshow(img);
result=(img2<1).*img;
subplot(1,2,2);imshow(result);

mais cela ne change rien !!!!

[Gooby]

Ton image de base est en niveau de gris ou en RGB?

[hulia]

en niveau de gris !!!

[Gooby]

Alors personnellement, je sais pas si c'est une méthode super rigoureuse, mais j'essayerai par tâtonnement de trouver le seuil qui te permettrait dans un premier temps de discriminer l'os (et a forciori le prothèse) du fond de ton image.
Normalises ton image entre 0 et 1 si ce n'est pas déjà fait, puis tente peu à peu de binariser l'image avec un seuil de plus en plus haut et enfin, analyse le résultat visuel pour voir si tu as bien isoler l'os et la prothèse.

Code :

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

1
2
test = im2bw(img_normalisee,monSeuil); % augmentation du seuil pas à pas
imshow(test) %regarde si tu as bien récupérer ton os avec la prothèse

Ensuite tu pourras augmenter significativement ton seuil pour n'obtenir que ta prothèse dans l'image binaire.

[hulia]

et comment normaliser ?

[hulia]

est-ce correct de normaliser ainsi ?

Code :

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

1
2
3
4
5
6
clear all
close all
clc
img=imread('nb.tif');
 
Im=double(img)/255; %ici j'ai normalisé l'image pour ne garder que les 0 et 1

[Gooby]

pour normaliser tu peux faire:

Code :

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

1
2
img = double(img);
img = (img-min(min(img)))/(max(max(img))-min(min(img)));

[hulia]

merci bcp, j'ai bien discriminé ma prothèse , maintenant je dois chercher une méthode pour sélectionner juste les points de la surface supérieur dans e but de calculer l'air ???

[Gooby]

Bah pour l'aire tu peux regarder le 2eme paramètre de sorties de bwboundaries qui te renvoie le masque de ton objet. En fait je comprends pas très bien ce que tu veux dire par "aire" dans ce cas ci. Tu veux faire un calculer de dimension? as tu la correspondance entre un pixel et une dimension réel?

[hulia]

j'ai trouvé une fonction qui le fait "bwboundaries"
voici mon code ainsi que la figure retrouvée

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
 
clear all
close all
clc
img=imread('nb.tif');
 
img = double(img);
img = (img-min(min(img)))/(max(max(img))-min(min(img)));
test = im2bw(img,0.3); % augmentation du seuil pas à pas
imshow(test) %regarde si tu as bien récupérer ton os avec la prothèse
h=fspecial('gaussian',5,5);
bw=imfilter(test,h);
bw=imfill(~bw,'hole');
bw=~bw;
b=bwboundaries(bw);
hold on
for i=1:length(b)
    contour=b{i};
    plot(contour(:,2),contour(:,1),'g','LineWidth',3);
end

Cependant je veux sélectionner uniquement certains points de surface de la prothèse ?!!!!

[hulia]

Voilà j'ai réglé mon problème concernant le contour , maintenant je veux savoir comment sélectionner juste les points de surface pour calculer la surface supérieure .

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
 
clear all
close all
clc
img=imread('nb.tif');
img = double(img);
imshow(img)
img = (img-min(min(img)))/(max(max(img))-min(min(img)));
test = im2bw(img,0.3); % augmentation du seuil pas à pas
imshow(test) %regarde si tu as bien récupérer ton os avec la prothèse
h=fspecial('gaussian',5,5);
bw=imfilter(test,h);
bw=imfill(~bw,'hole');
bw=~bw;
b=bwboundaries(bw);
hold on
for i=3:length(b)-2
    contour=b{i};
    plot(contour(:,2),contour(:,1),'g','LineWidth',3);
end

[betaplus]

Si tu pars du principe que la surface de ton objet est égale à la somme des surfaces des points qui le constituent ça peut aider, ou alors j'ai mal compris ce que tu veux faire.

En claire, mon idée est de dire que :

• tu connais le nombre de points au dessus de ton seuil
• tu connais la surface d'un point
• Donc tu multiplie le nombre de points par la surface d'un seul points

Betaplus