在numpy中使用导数的三次样条插值

Question

5

最近我开始学习Python和使用Python进行科学计算。这是一个问题，为了避免可能已经存在的代码重复。

我有一个采样为x和y函数的字段，它们在一个规则的网格上。我想要插值数据，不仅可以得到网格上任意点的场值，还能得到一阶和二阶导数。如果我使用interp2d来进行双三次插值，我可以得到场的值。

是否有人建议使用现有的numpy或scipy函数来获取场的一阶和二阶导数？

谢谢！

- John Doe

2

你是否在寻找 numpy.diff？ - mmgp

2

http://docs.scipy.org/doc/numpy/reference/generated/numpy.gradient.html - YXD

2

虽然 UnivariateSpine 有一个返回导数的方法，但对于二维样条插值，我猜你目前最好的选择是通过 get_coeffs 获取系数，然后从那里开始。 - ev-br

1

@MrE numpy.gradient 也可能完全适用。 - John Doe

@MrE，您可以将您的评论发布为答案... - Saullo G. P. Castro

显示剩余2条评论

1个回答

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- Yannick Copin · Accepted Answer

scipy.interpolate.interp2d.__call__ 方法提供了 dx 和 dy 选项，用于在点上求取更高阶导数（至少从版本 0.14.0 开始可用）。