interpn griddata difference?

Question

Tahariet Sharon 2020-4-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/517853-interpn-griddata-difference

评论： Tahariet Sharon 2020-4-17

采纳的回答： Bjorn Gustavsson

Hi,

What are the differences between interpn and griddata?

Thank you.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Bjorn Gustavsson 2020-4-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/517853-interpn-griddata-difference#answer_426034

The interpn function requires input data on a plaid grid, i.e. as you would produce with meshgrid, or perhaps ndgrid, while griddata (and scatteredInterpolant and TriScatteredInterp) expects input data on irregular grids. GRIDDATA expects the inputs as 1-D vectors. See the help and documentation for the respective functions, there you should find "minimal" illustrating examples to look at.

HTH

5 个评论
显示 3更早的评论隐藏 3更早的评论

Bjorn Gustavsson 2020-4-15

在 MATLAB Online 中打开

No. interpn can only handle grids as produced by ndgrid (I recalled that it could also handle grids as produced by meshgrid, but that was not the case for the matlab-version on my laptop, this might have been the case for "very" old matlab-versions.). That means that the grids for interpn might have step-sizes that varies along each separate dimension:

[x,y,z] = ndgrid([0 1 2],[-3 -1 0 1 3],[0 1 2 4 8]);
[xi,yi,zi] = meshgrid(0:.1:2,-3:.1:3,0:.1:8);       
I = exp(-(x-1).^2-(y-1).^2-(z-3).^2);
Ii = interpn(x,y,z,I,xi,yi,zi);      
slice(xi,yi,zi,Ii,1,1,3),shading flat

That grid is still "plaid" - Each "cell" [dx,dy,dz] would still be cube-ish (dont remember name of that geometric shape in English). griddata and the other triangulation-based interpolation-functions handle interpolation of properly irregularly sampled data:

xR = 1 + randn(75,1);
yR = 2*randn(75,1);    
zR = 4 + 2.5*randn(75,1);
Ir = exp(-(xR-1).^2-(yR-1).^2-(zR-3).^2);
Igd = xi;
Igd(:) = griddatan([xR,yR,zR],Ir,[xi(:),yi(:),zi(:)]);
slice(xi,yi,zi,Igd,1,1,3),shading flat
hold on
scatter3(xR,yR,zR,32,log(Ir),'filled')

That type of data interpn cannot handle.

HTH

Bjorn Gustavsson 2020-4-16

Yeah, that might very well be a slight performance-booster (check such things with the profiler-functionality if you need speed-ups...). I was too lazy to do that since my example was rather small and I only wanted to plot the data with the slice function, which used to expect the data to be in the meshgrid format.

The one advice when it comes to this meshgrid/ndgrid option is DOCUMENT YOUR CHOICE CLEARLY (and I'm shouting mainly for me to remember this...), the speed-up should be rather minor - typically, I think I've lost more time hunting down where I've done which choice than I've gained or lost due to this performance difference.

HTH

Tahariet Sharon 2020-4-17