python - scipy interpolate 给出无界值

Question

我有一个数据集data[xi,yi,zi]我想绘制(带有插值)。使用 scipy.interpolate，一切看起来几乎完美，但是插值生成的一些值超出了输入数据的范围。例如假设 zi受 0 < zi < 1 约束， rbf 插值似乎返回超出范围的插值(例如 >1 )；这是我的简化尝试:

N=100
data=[xi yi zi]
xis = np.linspace(xi.min(), xi.max(), N)
yis = np.linspace(yi.min(), yi.max(), N)
XI, YI =  np.meshgrid(xis,yis)
rbf = scipy.interpolate.Rbf(xi, yi, zi, function='linear')
ZI=rbf(XI,YI)
print ZI.max()
->1.01357328514

有没有办法将限制传递给 Rbf 并让它知道不要超过 zi.max() 和 zi.min() ？

Answer 1

使用径向基函数进行插值可能会导致值高于给定数据值的最大值并低于给定数据值的最小值。 (Illustration)。这是该方法的一个数学特征，无法传入禁用它的选项。两种可能的解决方案:

• 使用np.clip绘制数据时，在数据的最小-最大值之间剪裁插值。
• 改用分段线性插值 ( scipy.interpolate.LinearNDInterpolator )，这样可以保证遵守数据值的最小值和最大值。