拟合圆并求圆心(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;
           <