Speeding up interp2

27 次查看（过去 30 天）
Philip Freeman
Philip Freeman 2022-1-20
评论： Philip Freeman 2022-1-20
I'm performing 2-D interpolation on a large gridded dataset. The method of call I'm using is
vals = interp2(V,Xq,Yq)
Where V is a large image, but Xq and Yq are fairly sparse in comparison to the image size. Xq and Yq may be invalid, so I do need to preserve the NaN returns when they are.
The code is painfully slow and running the profiler shows that line 105 in interp2 is the culprit
V = V.';
is taking 82% of the execution time! (Not sure why the author chose to transpose V rather than just swap the index arrays, but whatever, that's a Matlab "bug" I guess.)
Is there an alternate way to do this without interp2 so that I can avoid the unnecessary and wasteful transpose? Should I just call makegriddedinterp() directly the way it looks like interp2 does anyway (and swap X and Y in the call)?
  3 个评论
显示 1更早的评论
Sean de Wolski
Sean de Wolski 2022-1-20
In general using griddedInterpolant will also be faster for subsequent interpolations because you build it once, it does the underlying computations once, and then are reused. This only helps if you need to call multiple times though.
Philip Freeman
Philip Freeman 2022-1-20
Thanks, Sean. That is what I ended up doing and the performance is greatly improved.

请先登录，再进行评论。

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

回答（0 个）

类别

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

标签

产品


版本

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by