Curve (Line) Fitting in 3D

24 次查看（过去 30 天）
Alec Day
Alec Day 2017-4-27
编辑： Sarvesh Singh 2018-2-15
Hi, I've noticed John D'Errico answering a few of these questions but the polyfitn doesn't seem at first trial to do what I want. But I have some data in x,y,z coordinates. See attached for coordinates. I have manually selected x,y from images based on z = [0:2:100]. However I'm not sure how to fit a curve to this in 3D in order to remove error due to manual selection.
  2 个评论
Star Strider
Star Strider 2017-4-27
What function do you want to fit to your data?
Alec Day
Alec Day 2017-5-1
I've ended up using polyfitn to find z = f(x) and z = f(y) to define a line in 3D space. seems to fit pretty well between z = 0-100.

请先登录，再进行评论。

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

采纳的回答

John D'Errico
John D'Errico 2017-4-27
编辑：John D'Errico 2017-4-27
Here is your data, plotted:
PLOT EVERYTHING.
You have what appears to be a relationship where the dominant thing happening is a curve in (x,y,z). Think of it as a path through space, with a cloud of noise around that path.
NO. You CANNOT use polyfitn here. There is insufficient information to build a viable model, of the form z(x,y). Note that z(x,y) would be a surface, defined everywhere. To try to infer the coefficients of such a model would be essentially impossible to do well. If you try to use polyfitn, the first thing you would then do is send me a plaintive e-mail, asking why polyfitn does not do well on your data! Similarly, do NOT use a tool like gridfit. It too will fail to work well here, as your data is clustered fairly tightly around this path in the (x,y,z) space.
Hey, I'm just trying to forestall the inevitable plaintive e-mails!
Again, you have what appears to be a 3-d space curve , with noise that varies along the path, sometimes called heteroscedastic noise. Yes, it hurts my eyes just to try to type the word. But it just means that the noise distribution varies spatially, so the variance changes with location.
Next, it appears as if you have a problem with errors in all three variables. Nonlinear errors in variables problems can be nasty.
The real question is, what will you do with the resulting model? Are you looking for a nice smooth curve that goes through the middle of that banana shaped blob, just for a plot? Are you looking for predictive ability of some sort? Or are you really hoping for some sort of surface? (If the latter, abandon all hope, ye who enter here. You will need to get better data for that.)
  3 个评论
显示 1更早的评论
Alec Day
Alec Day 2017-4-30
Hi John, thanks for the very detailed and prompt response. What i need is as you wrote 'a nice smooth curve that goes through the middle of that banana shaped blob' I don't need to be able to predict f(x,y) outside of the defined values for z (0-100 as plotted). But just a way of having a smooth curve somewhere through the middle, which should be able to be repeatable given the same data-points.
John D'Errico
John D'Errico 2017-5-1
As I said, this is tough to do. I once looked at some ideas for code to post for the general problem, but I could not make something work. What seems relatively easy by eye is not really so easy. That is often the case. It is sort of a 3-d errors in variables, with points projected onto a curve, which must then be a spline of some kind. It can be done, but it won't be robust/stable at all.

请先登录，再进行评论。

更多回答（1 个）

Steffen B. Petersen
Steffen B. Petersen 2017-12-16
I read your question such that you have a large number of points in 3D (x,y,z), and you want to fit a line through these points.
You should principal component analysis : First you 'whiten' your data by subtracting the mean of you data from each data point. Let us call your mean 'myaverage' - this is point in 3D space at the center of your data. Let us call the whitened data (X,Y,Z)
(X,Y,Z)=(x,y,z)-myaverage
Second you use pca to extract the coefficients - in this case a 3x3 matrix where the first column is the vector that has the largest varians over your dataset.
[coefficient]=pca(X,Y,Z);
You can now draw a line through your datapoints using
mypoint=myaverage+t*coefficient(:,1)
by varying the parameter t you find the points for your fitted line. From this point off you may want to compute the RMS deviation - but that probably is not your problem.
Regards
Steffen
  1 个评论
Sarvesh Singh
Sarvesh Singh 2018-2-15
编辑：Sarvesh Singh 2018-2-15
Hi Steffen,
I'm trying to implement this. Can you please elaborate what this parameter 't' is and what should be its typical value for the best fit. Also is it possible to fit any other curve instead of line using this method?
Regards,
Sarvesh

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Get Started with Curve Fitting Toolbox 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by