Data fitting using thin-plate spline/interpolation

33 次查看（过去 30 天）
Pelajar UM
Pelajar UM 2023-2-1
评论： William Rose 2023-2-3
I have several points in 3D space. I have used the surface fitting function to create a pretty accurate representation.
Now I need to generate an equation with the corresponding coefficients for this surface, so that I can use it elsewhere without having the original data points.
1) What is the equation for the thin-plate interpolation? Like when you use polynomial, sine, fourier, etc, there is a clear equation that describes the data....
2) How do I extract the coefficients? Coeffx only returns 'thinplateinterp'.
x=[0,0,0,0.5,0.5,0.5,1,1,1];
y=[0.001,0.1,10,0.001,0.1,10,0.001,0.1,10];
z=[4.35,1.37,2.31,7.04,3.58,4.14,2.08,2.35,2.43];
ylog=log10(y);
[xData, yData, zData] = prepareSurfaceData( x, ylog, z );
ft = 'thinplateinterp';
[fitresult, gof] = fit( [xData, yData], zData, ft, 'Normalize', 'on' );
Coeffx=coeffvalues(fitresult);

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

采纳的回答

William Rose
William Rose 2023-2-1
@Pelajar UM,
I too am unable to find the coefficients for the thin plate psline model.
I can evaluate the fitted model in two ways:
>> fitresult(.5,-3)
ans =
7.0400
>> evaluate(Coeffx,.5,-3)
ans =
7.0400
The reult with "evaluate(Coeffx,.5,-3)" gives the answer above only if I compute the fit with 'Normalize','off'.
If I compute the fit with 'Normlaize','On', then the output from fitresult(.5,-3) is correct, and is the value shown above, but the output from evaluate(Coeffx,.5,-3) is incorrect.
The wikipedia article for thin plate spline says that there are 2(K+3) parameters in the thin plate spline model, where K=number of fitted points. In your example, K=9, so there should be 24 parameters.
I did the fit in the CurveFitter app also, with the same results. I was also unable to find the coefficients in the curve fitter app, and when I selected the "export to workspace" option, it just exported variables like fitresult, so this did not help identify the parameters.
One can also fit a thin plate spline to the same data by:
p=1; st=tpaps([xData',yData'],zData',p)
The resulting structure, st, includes st.coefs, which is a vector of 12 coefficients. By setting p=1, tpaps() returns an exact-fitting thin plate spline. As p decreases toward 0, the fit becomes increasingly approximate. See the help for tpaps().
Good luck with your fitting.
  10 个评论
显示 8更早的评论
Pelajar UM
Pelajar UM 2023-2-3
Last update:
Long story short, it is better that you use fitrsult to generate new data points in the desired range and then just feed these into a high-degree polynomial. It's not as accurate (also depends on how many points you are willing to generate) but in my case, I got R^2 = 0.935, by almost doubling the number of my data points. However one must note that the resulting polynomial gives absolutely nonsensical predictions outside the data range .
There is a way to implement the spline (and I think it involves denoramlizing the data after fitting), but it is just too much headache for a small increase in accuracy, especially if your end goal is to replicate this outside MATLAB.

请先登录，再进行评论。

更多回答（1 个）

John D'Errico
John D'Errico 2023-2-1
Ran in:
This is a common problem people have with splines in any form. They hope to see some nice simple function they can write down. Maybe put the function into a paper, or use in some way they would use other functional forms.
Sorry. That is pretty much never the case. Splines are great things for some purposes. But if you want to write down that nice looking function, it won't happen.
You can evaluate the function at any point. There are methods provided for that purpose. But if you really need the pretty function you can write down, you will need to use other tools.
x=[0,0,0,0.5,0.5,0.5,1,1,1];
y=[0.001,0.1,10,0.001,0.1,10,0.001,0.1,10];
z=[4.35,1.37,2.31,7.04,3.58,4.14,2.08,2.35,2.43];
ylog=log10(y);
surf(reshape(x,3,3),reshape(ylog,3,3),reshape(z,3,3))
Anyway, using a thin plate spline for such minimal data is probably wild overkill. At best, a simple polynomial model may be sufficient, and about as good as you can do with such limited data.
mdl = fit([x',ylog'],z','poly22')
Linear model Poly22: mdl(x,y) = p00 + p10*x + p01*y + p20*x^2 + p11*x*y + p02*y^2 Coefficients (with 95% confidence bounds): p00 = 1.457 (-0.9711, 3.886) p10 = 9.961 (-0.05327, 19.97) p01 = -0.03542 (-1.641, 1.57) p20 = -9.753 (-19.24, -0.2679) p11 = 0.5975 (-1.079, 2.274) p02 = 0.3229 (-0.2699, 0.9158)
This gives you a nice model you can write down. Not a very good model, but your data is highly limited.
[mdl(x,ylog);z]
ans = 2×9
4.4697 1.8156 1.7447 6.1156 4.0589 4.5856 2.8847 1.4256 2.5497 4.3500 1.3700 2.3100 7.0400 3.5800 4.1400 2.0800 2.3500 2.4300
As you can see, the model does predict the data reasonably well.
[xint,yint] = meshgrid(linspace(0,1),linspace(-3,1));
surf(xint,yint,mdl(xint,yint))
hold on
plot3(x,ylog,z,'ro')
You can even plot it as a nice smooth, well behaved surface.
  2 个评论
Pelajar UM
Pelajar UM 2023-2-1
Thanks for your help. I understand that my data is limited.
But can I generate more data from the spline fit and then use that data to come up with a more accurate polynomial fit?
Also the original question remains. Where are these coeffients stored and how would one use them? I am not necessarily looking for a simple equation, but the idea is to be able to generate this surface and extract data from it without having the original data points.
William Rose
William Rose 2023-2-1
@Pelajar UM, your question is a reasonable one, and I don't know the answer to where the coefficients can be found. I demonstrated with tpaps() how to find 12 coefficients. This is only half the number of coefficients expected, according to the wikipedia article on thin plate splines. I also demonstrated three ways you can use a spline fit to predict values at other points:
fitresult(.5,-3)
evaluate(Coeffx,.5,-3) %works if the fit is done with normalize off
fnval(st,[.5,-3]) %st is obtained with tpaps()

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Smoothing 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by