huiduxianshi.m

对电容层析成象的ECT断层图象 进行3维显示

for i=1:10
    for j=1:10
           if (I(i,j)==0)
           y=0.9-0.2*(i-1);
            x=-0.9+0.2*(j-1);
            fill([x-0.1 x-0.1 x+0.1 x+0.1],[y-0.1 y+0.1 y+0.1 y-0.1],'b');
            hold on
       end 
   end   
end
hold on
grid on
t = linspace(-2,2,40)*pi;
x = sin(t);
y = cos(t);
fill(x,y,'k');
axis auto
hold on
t=linspace(0,1,40)*pi;
x=sin(t);
for i=1:20
 y(i)=-sqrt(1-x(i).^2);
end;
for i=21:40;
    y(i)=-1;
end;
fill(x,y,'w');
hold on
t=linspace(-1,0,40)*pi;
x=sin(t);
for i=1:20
 y(i)=-sqrt(1-x(i).^2);
end;
for i=21:40;
    y(i)=-1;
end;
fill(x,y,'w');
hold on
t=linspace(-1,0,40)*pi;
x=sin(t);
for i=1:20
 y(i)=sqrt(1-x(i).^2);
end;
for i=21:40;
    y(i)=1;
end;
fill(x,y,'w');
hold on
t=linspace(0,1,40)*pi;
x=sin(t);
for i=1:20
 y(i)=sqrt(1-x(i).^2);
end;
for i=21:40
    y(i)=1;
end;
fill(x,y,'w');
axis off