matlab中用数据拟合圆心,拟合圆并求圆心(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([‘这个圆的半径是: ‘ num2str(R),‘ mm‘ ]);