Circle least squares fit for 3D data

38 visualizzazioni (ultimi 30 giorni)
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 Commenti
Ciara Gibbs
Ciara Gibbs il 7 Ago 2019
@Bruno Luong , please find the data attached in the order of x, y ,x . Thanks again
Ciara Gibbs
Ciara Gibbs il 7 Ago 2019
@Torsten, thank you for the links I will give them a read!

Accedi per commentare.

Risposta accettata

Bruno Luong
Bruno Luong il 7 Ago 2019
Modificato: Bruno Luong il 7 Ago 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 Commenti
Ciara Gibbs
Ciara Gibbs il 7 Ago 2019
Amazing! Thank you very much it is greatly appreciated and you have saved me a lot of time!
Martin Fuchs
Martin Fuchs il 24 Set 2020
This script is just brilliant!

Accedi per commentare.

Più risposte (0)

Categorie

Scopri di più su MATLAB in Help Center e File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by