numpy中的导数的三次插值

Question

最近移植到Python和Python的科学计算.这是一个避免任何可能已存在的代码重复的问题.

我有一个字段,它在常规网格中作为x和y的函数进行采样.我想插入数据,不仅要获得网格上任何一点的字段值,还要获得一阶和二阶导数.如果我使用interp2d使用双三次插值对其进行插值,我可以获得该字段的值.

有没有人建议如何使用现有的numpy或scipy函数获得该领域的一阶和二阶导数？

谢谢!

Answer 1

该scipy.interpolate.interp2d.__call__方法具有选项dx并dy可以在点处评估更高的导数（至少自版本 0.14.0 以来）。