recognize the center and diameter of a circle

I have a 2D scatter, it is a circle but the center and diameter is not given.

I need its center and diameter but as you can see in the figure, there are some noise points around. is there anyway to do this?

Answer 1

[Akira Agata]

If you have Optimization Toolbox, you can solve this nonlinear least-square problem by simply using lsqnonlin function, like:

load('circle_coordinate.mat');

x = circle_coordinate(:,1);
y = circle_coordinate(:,2);

% Solve as nonlinear least-squares problem
f = @(a) (x-a(1)).^2 + (y-a(2)).^2 - a(3).^2;
a0 = [mean(x),mean(y),max(x)-mean(x)]; 
af = lsqnonlin(f,a0);

% Fittig circle
theta = linspace(0,2*pi);
xFit = af(1)+af(3)*cos(theta);
yFit = af(2)+af(3)*sin(theta);

% Visualize the result
figure
scatter(x,y,'.')
hold on
plot(af(1),af(2),'rx')
plot(xFit,yFit,'r-')
legend({'Data','Center','Fitting circle'})

[Yu Li]

fantastic! I have never know Matlab can do this before.

Thank you so much!

Yu