python - numpy interp 减少 xp

Question

我想插入一个 numpy 数组，np.interp 几乎完全符合我的要求:

interp(x, xp, fp, left=None, right=None)
One-dimensional linear interpolation.

除了这个位:

Does not check that the x-coordinate sequence xp is increasing. If xp is not increasing, the results are nonsense.

我的 xp 正在减少，所以哪个更好:反转 xp 和 fp 的方向:

np.interp(x, xp[::-1], fp[::-1])

或反转 x 和 xp:

np.interp(-x, -xp, fp)

或者有更好的方法吗？

Answer 1

感谢所有提供意见的人，尤其是@Jaime。

我做了一些实验，得出了这个结论:

1)除了舍入误差外，我提到的两种方法结果相同。

2) 它们花费的时间几乎相同

3) 我尝试了 scipy 版本，但它会拒绝 assume_sorted 标志。也许我的 scipy 版本已经过时了。我怀疑如果升起该标志，scipy 会在内部对数组进行排序。但是值是按相反方向排序的，因此不需要这样做。

无论如何，我将使用反向方法:

np.interp(x, xp[::-1], fp[::-1])

请记住，在这种情况下，如果需要，您还必须反转 left 和 right。