如何使用gstools的外部漂移克里金插值二维地理空间数据。 Gstools只有1-D外部驱动克里格插值示例,但我想使用Gstools的外部驱动克里格插值来插值2-D地理空间数据
这是一个官方的例子:
import numpy as np
from gstools import SRF, Gaussian, krige
# synthetic condtions with a drift
drift_model = Gaussian(dim=1, len_scale=4)
drift = SRF(drift_model, seed=1010)
cond_pos = [0.3, 1.9, 1.1, 3.3, 4.7]
ext_drift = drift(cond_pos)
cond_val = ext_drift * 2 + 1
# resulting grid
gridx = np.linspace(0.0, 15.0, 151)
grid_drift = drift(gridx)
# kriging
model = Gaussian(dim=1, var=2, len_scale=4)
krig = krige.ExtDrift(model, cond_pos, x, ext_drift)
krig(gridx, ext_drift=grid_drift)
ax = krig.plot()
ax.scatter(cond_pos, cond_val, color="k", zorder=10, label="Conditions")
ax.plot(gridx, grid_drift, label="drift")
ax.legend()
我是GSTools的开发人员,当然也可以在2D和3D中执行外部漂移克里格法
在下面的示例中,我创建了x方向的人工正弦外部漂移
相关问题 更多 >
编程相关推荐