Discussion :

[goku59]

Bonjour,

je voudrais savoir comment calculer le volume défini par deux contours définis sur deux coupes scanner par les coordonnées :

54.3665 2.0225 -12.0000
54.3665 42.0225 -12.0000
4.3665 42.0225 -12.0000
4.3665 2.0225 -12.0000
54.3665 2.0225 -14.0000
54.3665 42.0225 -14.0000
4.3665 42.0225 -14.0000
4.3665 2.0225 -14.0000

je pensais déterminer le nombre de voxels présent dans chaque contour grâce à la fonction poly2mask et connaissant la taille des voxels en déduire le volume mais ça ne fonctionne pas.
J'ai ensuite cherché à déterminer le nombre de voxels présent en reconstruisant le volume (un peu comme poly2mask mais en 3D) mais je n'y suis pas arrivé
Quelqu'un aurait il une idée???

D'avance merci!

[Jean Dumoncel]

Bonjour,

je ne suis pas sûr de comprendre : tes coordonnées ont l'air de définir un cube, c'est bien ça ? et tu cherches le volume du cube ?

[goku59]

En fait j'ai d'abord essayé avec un contour simple comme le cube pour voir si je trouvais la bonne valeur.
Mon objectif est de pouvoir déterminer un volume de forme quelconque.
Je travaille dans le domaine de l'imagerie médicale et ça serait pour connaître le volume de différents organes comme le foie, le coeur. ...
Merci pour votre aide.

[Kangourou]

Bonjour,

tu as plusieurs possibilités ;
1) une première est de calculer l'aire de la section visible sur chaque coupe. Si tu as des coupes sériées, j'imagine que c'est le type de données que tu obitns en sortie. Tu peux utiliser une formule pour calculer l'aire d'un polygone. Ensuite, il suffit de multiplier par l'écart entre les coupes.
2) une autre solution est de raisonner en 3D. Dans ce cas, si les données de départ sont les contours identifiés sur chaque coupe, il faut passer par une étape de reconstruction 3D, pour calculer le maillage polygonal (trimesh ou quadmesh) correspondant à l'organe. On peut alors appliquer des operations de lissage et / ou corerction, mais il y a assez peu d'outils sous Matlab pour le faire

A++