Pykrige：具有指定漂移的克里金法

Question

我不知道如何在插值中包含指定的漂移。我的代码如下：

x = np.array([0, 1, 2, 3, 4])
y = np.array([0, 1, 2, 3, 4])
values = np.array([1, 2, 3, 5, 4])
drift = np.array([-14, 2, 6, 1, 3])

gridx = np.linspace(0, 4, 10)
gridy = np.linspace(0, 4, 10)

UK = UniversalKriging(
    x,
    y,
    values,
    variogram_model="linear",
    drift_terms=["specified"],
    specified_drift = [drift]
)

z, ss = UK.execute("grid", gridx, gridy)
plt.imshow(z)
plt.show()

我收到一个值错误

ValueError: Must provide drift values for kriging points when using 'specified' drift capability.

Answer 1

如果您想使用指定的漂移，您还需要在目标网格处使用

specified_drift_arrays

中的

UK.execute

参数提供它（请参阅文档）

grid_drift = ... # external drift at target locations
z, ss = UK.execute("grid", gridx, gridy, specified_drift_arrays=[grid_drift])

请注意，由于可能存在多个漂移，因此需要在列表中给出漂移。

希望这有帮助。