TriScatteredInterp doesn't extrapolate outside scattered data points

2 views (last 30 days)
hello everybody
I wonder why Mathworks considers TriScatteredInterp(x,y,z) as an extrapolation function while it only interpolates inside the area defined by "x,y" scattered points and it doesn't calculate outside that area.
Am I doing something wrong, here is the code:
pos = [ 401.5 387.5 659.5 523.5 349.0 509.5;...
304.0 305.0 405.5 205.0 483.5 229.5];
z1 = rand(1,6);
F1 = TriScatteredInterp(pos(1,:)',pos(2,:)',z1');
gx = 300:5:700;
gy = 200:5:500;
[qx,qy] = meshgrid(gx,gy);
z_1 = F1(qx,qy);

Accepted Answer

John D'Errico
John D'Errico on 16 Nov 2015
Edited: John D'Errico on 16 Nov 2015
These tools use a triangulation of your data. They only are able to interpolate if your data is inside one of the triangles. And since a delaunay triangulation only extends out to the convex hull of your data, they do not extrapolate. Anyway, extrapolation is dangerous business.
Note that scatteredInterpolant can apparently extrapolate linearly based on the gradients at the edges, but that extrapolation will be poor in general. Do not trust it for far outside of your data. If you do not have the current release, then TriScatteredInterp does not support extrapolation at all.
Here is the help for scatteredInterpolant for extrapolation. Note that TriScatteredInterp never claims to extrapolate that I see anywhere in the help.
F = scatteredInterpolant (..., METHOD, EXTRAPOLATIONMETHOD) supports the
selection of an extrapolation method to be used outside the convex hull.
Where EXTRAPOLATIONMETHOD is one of the following:
'nearest' - Evaluates to the value of the nearest neighbor on the
boundary (default for METHOD ='nearest')
'linear' - Performs linear extrapolation based on boundary gradients
(default for METHOD = 'linear' and METHOD='natural')
'none' - Queries outside the convex hull will return NaN
  1 Comment
Abdelmoumen Bacetti
Abdelmoumen Bacetti on 16 Nov 2015
Ohh, actually i won't trust any extrapolation unless I can cross validate it.
All I want, is to have values everywhere over my grid, and only the part inside the hull interests me.
Thanks anyway, and I guess i should use the trial version of R2015 to see it working.

Sign in to comment.

More Answers (2)

Walter Roberson
Walter Roberson on 16 Nov 2015
The link you posted refers only to scatteredInterpolant(), not to TriScatteredInterp() . TriScatteredInterp() is still present in MATLAB but the recommendation is to replace it with calls to scatteredInterpolant()
  1 Comment
Abdelmoumen Bacetti
Abdelmoumen Bacetti on 16 Nov 2015
Edited: Abdelmoumen Bacetti on 16 Nov 2015
scatteredInterpolant() seems doesn't exist in R2012.
I can't test if it does extrapolation in a grid larger than the data points.

Sign in to comment.


Abdelmoumen Bacetti
Abdelmoumen Bacetti on 18 Nov 2015
Hello
I found a toolbox that does an amazing job
Thanks for your contributions, and I hope you find this toolbox useful as I did.

Categories

Find more on Bounding Regions in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!