Circle least squares fit for 3D data

20 次查看（过去 30 天）
Ciara Gibbs
Ciara Gibbs 2019-8-7
评论： Martin Fuchs 2020-9-24
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 个评论
显示 6更早的评论
Ciara Gibbs
Ciara Gibbs 2019-8-7
@Bruno Luong , please find the data attached in the order of x, y ,x . Thanks again
Ciara Gibbs
Ciara Gibbs 2019-8-7
@Torsten, thank you for the links I will give them a read!

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Bruno Luong
Bruno Luong 2019-8-7
编辑：Bruno Luong 2019-8-7
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 个评论
显示 2更早的评论
Ciara Gibbs
Ciara Gibbs 2019-8-7
Amazing! Thank you very much it is greatly appreciated and you have saved me a lot of time!

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Surface and Mesh Plots 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by