Circle least squares fit for 3D data

38 views (last 30 days)
Hi everyone,
I have 6000 x coordinates, y coordinates and z cooridinates that form a circle that does not perfectly occupy one plane. I am looking to fit the data with a circle but I can only find functions online that do it for x and y coordinates and not including z. I would appreciate any help in creating some code for this as I am not sure where to start (I am a beginner in MATLAB!)
Thank you in advance
Ciara
  8 Comments
Ciara Gibbs
Ciara Gibbs on 7 Aug 2019
@Bruno Luong , please find the data attached in the order of x, y ,x . Thanks again
Ciara Gibbs
Ciara Gibbs on 7 Aug 2019
@Torsten, thank you for the links I will give them a read!

Sign in to comment.

Accepted Answer

Bruno Luong
Bruno Luong on 7 Aug 2019
Edited: Bruno Luong on 7 Aug 2019
circle.png
X=csvread('data.csv');
XC = mean(X,1);
Y=X-XC;
[~,~,V]=svd(Y,0);
Q = V(:,[1 2]); % basis of the plane
Y=Y*Q;
xc=Y(:,1);
yc=Y(:,2);
M=[xc.^2+yc.^2,-2*xc,-2*yc];
% Fit ellipse through (xc,yc)
P = M\ones(size(xc));
a=P(1);
P = P/a;
r=sqrt(P(2)^2+P(3)^2+1/a); % radius
xyzc = XC' + Q*P(2:3); % center
theta = linspace(0,2*pi);
c = xyzc + r*Q*[cos(theta); sin(theta)]; % fit circle
close all
plot3(X(:,1),X(:,2),X(:,3),'.');
hold on
plot3(c(1,:),c(2,:),c(3,:),'r','LineWidth', 2);
axis equal
  4 Comments
Ciara Gibbs
Ciara Gibbs on 7 Aug 2019
Amazing! Thank you very much it is greatly appreciated and you have saved me a lot of time!
Martin Fuchs
Martin Fuchs on 24 Sep 2020
This script is just brilliant!

Sign in to comment.

More Answers (0)

Categories

Find more on MATLAB in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!