stopr1.m

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

function  stopr1(A1,B1,C1,A2,B2,C2,m,n) 
p1x=0;
p1z=0;
p2x=0;
p2z=0;
for j=1:m
    p1x=p1x+A1(j);
    p1z=p1z+C1(j);
end
    p1x=p1x/m;
    p1z=p1z/m;
for i=1:n
    p2x=p2x+A2(i);
    p2z=p2z+C2(i);
end
    p2x=p2x/n;
    p2z=p2z/n;
%求中心点
dx=p1x-p2x;
dz=p1z-p2z;
min=sqrt((A1(1)-A2(1)-dx)^2+(B1(1)-B2(1))^2+(C1(1)-C2(1)-dz)^2);
s=1;
h=1;
a1=1;
b1=1;
a2=2;
b2=2;
for j=1:n
    d=sqrt((A1(1)-A2(j)-dx)^2+(B1(1)-B2(j))^2+(C1(1)-C2(j)-dz)^2);
  if d<min
      min=d;
       b1=j;
  end
end
h=b1;
if b1==n
   b2=1;
else 
   b2=b1+1;
end
k=0;
while  ~((a2==s&b1==h)|(a1==s&b2==h))&k<m+n
    d1=sqrt((A1(a1)-A2(b2)-dx)^2+(B1(a1)-B2(b2))^2+(C1(a1)-C2(b2)-dz)^2);
    d2=sqrt((A1(a2)-A2(b1)-dx)^2+(B1(a2)-B2(b1))^2+(C1(a2)-C2(b1)-dz)^2);
    if d1<d2
          x=[A1(a1);A2(b1);A2(b2)];
          y=[B1(a1);B2(b1);B2(b2)];
          z=[C1(a1);C2(b1);C2(b2)];
          b1=b2;
          if b1==n
             b2=1;
          else 
             b2=b1+1;
          end  
          k=k+1;   
   else 
         x=[A1(a1);A2(b1);A1(a2)];
         y=[B1(a1);B2(b1);B1(a2)];
         z=[C1(a1);C2(b1);C1(a2)];
         a1=a2;
         if a1==m
            a2=1;
         else 
            a2=a1+1;
         end
        k=k+1;
  end       
  c=[1;1;1];
  fill3(x,y,z,c);
  hold on
end   
if a1==s&b2==h
   x=[A1(a1);A2(b1);A2(h)];
   y=[B1(a1);B2(b1);B2(h)];
   z=[C1(a1);C2(b1);C2(h)]; 
   c=[1;1;1];
   fill3(x,y,z,c);
   hold on
   else
   x=[A1(a1);A2(b1);A1(s)];
   y=[B1(a1);B2(b1);B1(s)];
   z=[C1(a1);C2(b1);C1(s)]; 
   c=[1;1;1];
   fill3(x,y,z,c);
   hold on
   
        
   end