contour.m

来自「对电容层析成象的ECT断层图象 进行3维显示」· M 代码 · 共 57 行

function bw=Contour(I,t);
%求轮廓,I为矩阵,t为层编号,适用三维数据处理,其中z,x为决定性二维平面

%I=[ 0     0     0     0     0     0     0     0     0     0;
    % 0     0     0     0     1     0     0     0     0     0;
    % 0     0     0     1     1     1     0     0     0     0;
     %0     0     0     0     1     0     0     0     0     0;
     %0     0     0     0     0     0     0     0     0     0;
     %0     1     1     1     0     0     1     1     1     0;
     %1     1     1     0     0     0     1     1     1     0;
     %0     0     0     0     0     0     1     1     1     0;
     %0     0     0     0     0     0     0     0     0     0;
     %0     0     0     0     0     0     0     0     0     0];
 %t=1;


bw1=zeros(10,10);
X=zeros(10,10);
Y=zeros(10,10);
Z=zeros(10,10); 
for i=1:10
    for j=1:10          
            X(i,j)=-0.9+0.2*(j-1);
            Y(i,j)=t*0.5;   
            Z(i,j)=0.9-0.2*(i-1);         
       end 
end

for i=1:10
    for j=1:10
        if(X(i,j)^2+Z(i,j)^2<=1)
            bw1(i,j)=1;
        end
    end
end 

I=I&bw1;
bw=bwperim(I);

%bw
for i=1:10
    for j=1:10
           if (bw(i,j)==1)
           z=0.9-0.2*(i-1);
            x=-0.9+0.2*(j-1);
            y=Y(i,j);
            plot3(x,y,z,'k');
            hold on
        end
    end
end
hold on
axis auto
axis equal
%X
%Y
%Z