Circle least squares fit for 3D data

42 次查看（过去 30 天）

Ciara Gibbs 2019-8-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/475212-circle-least-squares-fit-for-3d-data

评论： 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更早的评论隐藏 6更早的评论

Ciara Gibbs 2019-8-7

data.csv

@Bruno Luong , please find the data attached in the order of x, y ,x . Thanks again

Ciara Gibbs 2019-8-7

@Torsten, thank you for the links I will give them a read!

请先登录，再进行评论。

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

采纳的回答

Bruno Luong 2019-8-7

4
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/475212-circle-least-squares-fit-for-3d-data#answer_386587

编辑：Bruno Luong 2019-8-7

在 MATLAB Online 中打开

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