xueguan.rar_内切圆_立体图_血管 三维_血管三维重建_血管重建

clear all; for u=1:4 imname=sprintf('%d.png',u); % I=imread('2.png'); I=imread(imname); I1=im2bw(I,0.5); A=I1; figure; imshow(A); hold on [m,n]=size(A); t=1; X=[]; Y=[]; for i=1:m for j=1:n if(A(i,j)==0) if(A(i-1,j)==1||A(i+1,j)==1||A(i,j-1)==1||A(i,j+1)==1) X(t)=i; Y(t)=j; t=t+1; end end end end xmax=max(X); xmin=min(X); ymax=max(Y); ymin=min(Y); for i=ymin:ymax plot(i,xmin,'r'); plot(i,xmax,'r'); end for i=xmin:xmax plot(ymin,i,'r'); plot(ymax,i,'r'); end for i=1:t-1 plot(Y(i),X(i),'r'); end Rmax=0; r=[]; for i=xmin:xmax for j=ymin:ymax if((A(i,j)==0)&&(A(i-1,j)==0)&&(A(i+1,j)==0)&&(A(i,j-1)==0)&&(A(i,j+1)==0)) for k=1:t-1 r(k)=sqrt((i-X(k))^2+(j-Y(k))^2); end; Rt=min(r); if(Rmax<Rt) Rmax=Rt; Rx=i; Ry=j; Rc(1,u)=Rmax; Pxy(1,u)=Rx; Pxy(2,u)=Ry; end end end end R=Rmax; m=Rx-R:0.1:Rx+R; n=[]; n2=[]; p=length(m); for i=1:p n(i)=sqrt(R^2-(m(i)-Rx)^2)+Ry; n2(i)=-sqrt(R^2-(m(i)-Rx)^2)+Ry; end plot(n,m,'r',n2,m,'r'); plot(Ry,Rx,'r'); % for i=Rx-R:Rx+R % for j=Ry-R:Ry+R % if((Rx-i)^2+(Ry-j)^2==R^2) % plot(j,i); % hold on % end % end % end end