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(['这个圆的半径是: ' num2s