numpy中的导数的三次插值

Question

最近移植到Python并使用Python进行科学计算。这是一个避免任何可能已经存在的代码重复的问题。

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

是否有人建议如何使用现有的numpy或scipy函数获得该字段的一阶和二阶导数？

谢谢！

Answer 1

scipy.interpolate.interp2d.__call__方法有选项dx和dy来评估高阶导数（至少从版本0.14.0开始）。