lunkuo.m

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

function bw=lunkuo(I,t);
bw1=zeros(10,10);
X=zeros(10,10);
Y=zeros(10,10);
Z=zeros(10,10);
for i=1:10
  for j=1:10
     Z(i,j)=0.9-0.2*(i-1);
     X(i,j)=-0.9+0.2*(j-1);
     Y(i,j)=t*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);
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
t = linspace(-2,2,40)*pi;
x = sin(t);
z = cos(t);
y = Y(i,j);
y=y*ones(size(t));
fill3(x,y,z,'k');

axis equal