最近移民到 Python 和使用 Python 进行科学计算。这是一个避免任何可能已经存在的代码重复的问题。
我有一个字段在常规网格中作为 x 和 y 的函数进行采样。我想对数据进行插值,不仅可以获得网格上任意点的字段值,还可以获得一阶和二阶导数。如果我使用 interp2d 对它进行双三次插值,我可以获得该字段的值。
有没有人建议如何使用现有的 numpy 或 scipy 函数获得该字段的一阶和二阶导数?
谢谢!
最近移民到 Python 和使用 Python 进行科学计算。这是一个避免任何可能已经存在的代码重复的问题。
我有一个字段在常规网格中作为 x 和 y 的函数进行采样。我想对数据进行插值,不仅可以获得网格上任意点的字段值,还可以获得一阶和二阶导数。如果我使用 interp2d 对它进行双三次插值,我可以获得该字段的值。
有没有人建议如何使用现有的 numpy 或 scipy 函数获得该字段的一阶和二阶导数?
谢谢!
该scipy.interpolate.interp2d.__call__
方法具有选项dx
并dy
可以在点评估更高的导数(至少从版本 0.14.0 开始)。