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
| % liste des positions des pixels du volume
[lig,col] = find(tmp);
% nb de pixel du volume
nbPix = length(lig);
% déclaration matrice pr affichage 3D
mat = zeros(nbLig,nbCol);
% remplissage matrice pr affichage 3D
for j = 1:nbPix
% récup' coordonnées
l = lig(j);
c = col(j);
% affectation du numéro de ligne
mat(l,c) = nbLig - l;
end
% ajout matrice image au bloc matrice 3D
BW(:,:,i) = mat;
end
% affichage dans une sous-fenetre :
figure(25);
% lissage
% BW = smooth3(BW);
% affichage 3D :
isosurface(BW,1); |
Partager