Discussion :

[kamelcompte]

Bonjour tout le momde
je travail sur la segmentation des images en 3D, pour visualiser des objet que se trouvent dans le volume j'utilise la fonction isosurface du matlab.
ma question est ce que il ya une fonction ou un code matlab qui permet de calculer l'aire de la surface de l'objet en se basant sur l'isosurface.
lou précisamement

Code :

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

[faces, sommets]=isosurfaces(volume,isovaleur)

comment je peux calculer l'aire de la surface en fonction des sommets et des faces.

en merce d'avance

[Jerome Briot]

A priori les Faces sont des triangles.
Il suffit donc de calculer l'aire de chacun d'eux et de faire la somme, non ?

=> http://fr.wikipedia.org/wiki/Triangl....27un_triangle

[kamelcompte]

merci
mais comment je peux calculer l'aire des ces triangles à partir des facettes et des sommets
car lorsque j'applique la fonction isosurface j'obtient les faces qui sont des valeurs entiers et les sommets sont des valeurs décimale et que je sais pas la disposition des triangles

[Jerome Briot]

Regarde bien l'aide et les exemples de la fonctions PATCH

faces contient les indices dans sommets des points qui forment chaque triangle (un triangle pour chaque ligne de faces)

[kamelcompte]

merci j'arrive pas à calculer la surface d'un patch en plus les valeurs ses sommets ne sont pas des valeurs entiers naturels mais ils sont des valeurs réelles

[Jerome Briot]

merci j'arrive pas à calculer la surface d'un patch

Quel est les problème ?
Tu n'as toujours pas compris le rapport entre faces et sommets, ou bien tu ne sais pas calculer l'aire d'un triangle à partir de des coordonnées de ses sommets ?

en plus les valeurs ses sommets ne sont pas des valeurs entiers naturels mais ils sont des valeurs réelles

Je ne vois pas pourquoi le problème serait plus compliqué avec des nombres réels qu'avec des entiers ?

[kamelcompte]

j'ai pas accepté que les coordonnées des sommets soient des réelles parceque ces sommet correspondes à des voxels, et dans les images numériques chaque voxel est réperer par le triplet (i,j,k) dont i,j et k sont des entiers.
j'ai fait avec, et j'ai procedé comme suit:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
[nf,nv]=isosurface(volume,isovaleur);
aire=0;
for i=1:size(nf,1)  % pour chaque facette
    v1 = nv(nf(i,1),:); %sommet n°1
    v2 = nv(nf(i,2),:);% 2
    v3 = nv(nf(i,3),:);%3
    a=norm(v1-v2);  % ongueur du segment n°1 du triangle
    b=norm(v1-v3);%2
    c=norm(v2-v3);%3
    p = (a + b + c)/2;  % superperimeter
    A(i) = sqrt(p*(p - a)*(p - b)*(p - c));   % Heron's formula for triangle area 
    aire=(aire+A(i));
end

Cordialement