使用自定义的Variogram Model来替代autofitVariogram函数中的cutoff选项。
例如,使用指数模型:
library(automap)
library(gstat)
# 创建一个用于演示的数据集
data(meuse)
coordinates(meuse) <- ~x+y
data(meuse.grid)
coordinates(meuse.grid) <- ~x+y
# 定义Variogram模型(这里使用指数模型)
v_model <- vgm(psill = 0.4, model = "Exp", range = 800, nugget = 0.1)
# 计算半变异函数
mv <- autofitVariogram(log(zinc)~1, meuse, model = v_model)
# 进行插值(使用Ordinary Kriging)
ok <- autoKrige(log(zinc)~1, meuse, meuse.grid, model = mv$var_model)
通过自定义Variogram模型,可以避免cutoff选项无效的问题。