擬合圓並求圓心(matlab)

clc
 clear all
 I=imread('a003.bmp');
 [m,n,p]=size(I);
 I=im2bw(I);
 se=strel('square',3);
 Ia=imerode(I,se);
 Iout=I-Ia;
 figure(1);
 imshow(Iout);
 Iout1=Iout;
 N=1;
 for i=1:m
     for j=1:n
         if Iout(i,j)==1
             p(N,1)=i;
             p(N,2)=j;
             N=N+1;
             Iout(i,j+5:end)=0;
         end
        
     end
 end
 figure(2);
 imshow(Iout);
 for i=1:m
     for j=n:(-1):1
         if Iout1(i,j)==1
             p(N,1)=i;
             p(N,2)=j;
             N=N+1;
             Iout1(i,1:j-5)=0;
         end
     end
 end
 figure(3);
 imshow(Iout1);
 x=p(:,1);
 y=p(:,2);
 [R,x0,y0]=circ(x,y,N-1);
 disp(['這個