Calcul de volume d'un objet 3D

**mari_mimi** · 03/01/2015, 21h40

salut

,je voudrai savoir comment pourrai-je calculer le volume d'un objet 3D ,sachant que mon objet a été obtenu à partir du programme suivant:

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
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
a=dicomread('C:\Users\client\Documents\MATLAB\IM1.dcm');
% figure,imshow(a,[]);
 
b='E:\tp3\SE6\IM';
% for i=1:24
%     base=[b,num2str(i),'.dcm']
%     a=dicomread(base);
%     pause(0.001);
%    % imshow(a,[]); 
% end
 
 
   base=[b,num2str(10),'.dcm'];
    [a, map]=dicomread(base);
    %imshow(a,[]);
 % figure,imhist(a);
   a=a>650;
   a=bwareaopen(a,600);
   a1=imfill(a,'holes');
   %figure,imshow(a1,[]);
 
base=[b,num2str(8),'.dcm'];
    [a, map]=dicomread(base);
    %imshow(a,[]);
 % figure,imhist(a);
   a=a>600;
   a11=bwareaopen(a,320);
   a11=imfill(a11,'holes');
    a22=bwareaopen(a,400);
    a22=imfill(a22,'holes');
    a2=a11-a22;
 % figure,imshow(a2,[]);
 
 
   base=[b,num2str(9),'.dcm'];
    [a, map]=dicomread(base);
   % imshow(a,[]);
  %figure,imhist(a);
   a=a>600;
   a=bwareaopen(a,500);
   a3=imfill(a,'holes');
   %figure,imshow(a3,[]);
 
base=[b,num2str(11),'.dcm'];
    [a, map]=dicomread(base);
    %imshow(a,[]);
  %figure,imhist(a);
   a=a>600;
   a=bwareaopen(a,520);
   a4=imfill(a,'holes');
   %figure,imshow(a4,[]);
 
 
base=[b,num2str(12),'.dcm'];
    [a, map]=dicomread(base);
    %imshow(a,[]);
  %figure,imhist(a);
   a=a>600;
   a=bwareaopen(a,520);
   a5=imfill(a,'holes');
   %figure,imshow(a5,[]);
 
base=[b,num2str(13),'.dcm'];
    [a, map]=dicomread(base);
    %imshow(a,[]);
  %figure,imhist(a);
   a=a>600;
   a=bwareaopen(a,520);
   a6=imfill(a,'holes');
   %figure,imshow(a6,[]);
 
base=[b,num2str(14),'.dcm'];
    [a, map]=dicomread(base);
    %imshow(a,[]);
  %figure,imhist(a);
   a=a>600;
   a=bwareaopen(a,520);
   a7=imfill(a,'holes');
   %figure,imshow(a7,[]);
 
base=[b,num2str(15),'.dcm'];
    [a, map]=dicomread(base);
    %imshow(a,[]);
  %figure,imhist(a);
   a=a>600;
   a=bwareaopen(a,320);
   a8=imfill(a,'holes');
  % figure,imshow(a8,[]);
 
base=[b,num2str(16),'.dcm'];
    [a, map]=dicomread(base);
    %imshow(a,[]);
  %figure,imhist(a);
   a=a>520;
   a=bwareaopen(a,120);
   a9=imfill(a,'holes');
  %figure,imshow(a9,[]);
 
  nc=9;
  [n m]=size(a);
%        base=[b,num2str(i),'.dcm'];
   D=zeros(n,m,nc);
   D(:,:,1)=127*a1;
   D(:,:,2)=127*a2;
   D(:,:,3)=127*a3;
   D(:,:,4)=127*a4;
   D(:,:,5)=127*a5;
   D(:,:,6)=127*a6;
   D(:,:,7)=127*a7;
   D(:,:,8)=127*a8;
   D(:,:,9)=127*a9;
   Ds=smooth3(D);
   hiso = patch(isosurface(Ds,5),...
'FaceColor',[1,.75,.65],...
'EdgeColor','none');
hcap = patch(isocaps(D,5),...
'FaceColor','interp',...
'EdgeColor','none');
colormap(map)
lightangle(305,30);
set(gcf,'Renderer','zbuffer');
view(215,30);

merci d'avance

**Jean Dumoncel** · 05/01/2015, 09h57

Bonjour,

je suppose que ton objet 3d est la surface triangulée que tu affiches à la fin de ton programme. Tu devrais pouvoir trouver une fonction permettant de calculer le volume sur le file exchange :
http://www.mathworks.com/matlabcentr...rm=volume+mesh

Tu peux récupérer les sommets et facettes de ta surface en utilisant la fonction isosurface :

Code :

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