calcul des deux points les plus éloignés d'une image

Version imprimable

19/02/2015, 13h31
hulia

calcul des deux points les plus éloignés d'une image

Bonjour,
je veux calculer l'angle entre deux objets dans une image , j'aimerai avoir des pistes , j'ai commencé par discrétiser les deux objets et calculer leur contour , mais je n'arrive pas à trouver une idée pour calculer l'angle????
19/02/2015, 14h53
Jean Dumoncel

Bonjour,

aurais-tu une image pour illustrer l'angle que tu souhaites mesurer?
19/02/2015, 14h55
hulia

1 pièce(s) jointe(s)

Pièce jointe 169496
voici l'image en gros je veux mesurer l'angle entre l'os et la prothèse
19/02/2015, 14h57
hulia

je pars de l'idée qu'il faut définir l'axe verticale de l'os puis calculer le barycentre de la prothèse pour ensuite calculer l'angle entre ce barycentre et l'axe verticale de l'os.
Mais je n'ai aucune idée de comment faire ?

au début j'ai commencé par calculer le contour de la prothèse , mais je ne crois pas que cela peut m'aider !!!

Code:

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

19/02/2015, 17h59
Jean Dumoncel

Citation:

Envoyé par hulia

je pars de l'idée qu'il faut définir l'axe verticale de l'os puis calculer le barycentre de la prothèse pour ensuite calculer l'angle entre ce barycentre et l'axe verticale de l'os.

Un angle entre une droite et un point, cela n'a pas de sens d'un point de vue mathématiques (ou alors je n'ai pas compris).
Est-ce que tu ne chercherais pas plutôt l'angle entre l'axe de l'os et l'axe de la prothèse?
19/02/2015, 20h25
hulia

oui je cherche cette angle là , mais je ne sais pas comment calculer les axes !!!!
19/02/2015, 20h42
Jean Dumoncel

Avec la segmentation de la prothèse d'un côté et de l'os de l'autre, tu peux tester le calcul de l'orientation avec la fonction regionprops (et le paramètre 'orientation').
20/02/2015, 08h00
hulia

merci pour ta réponse .
J'arrive à segmenter ma prothèse par contre je n'arrive pas à segmenter l'os!!!!!