我正在使用krige()和interpolate()函数进行通用克里金法。对于predict()函数,我在“ parallel”包中使用clusterR()函数来减少进度。
此外,它还提到它也可以使用插值。但是,我无法对Kriging进行正确的编码。
eg)'预测'
begincluster(18)
clusterCub <- clusterR(muese.grid, predict, args=list(muese.Cub), filename="muese_pred.tif", format="Gtiff", progress="text")
endCluster()
but) 'interpolate'
formula <- OM ~ DEM + temp + NDVI
variogram
gUK (gstat)
UKdata <- krige(OM ~ DEM + temp + NDVI, location=meuse, model=variogram, newdata=meuse.grid)
UKintdata <- interpolate(rasters, gUK, xyOnly=False, index=1, filename="muese_pred.tif", format="GTiff")
clusterUK <- clusterR(rasters, interpolate, args=list(?????), filename="muese_pred.tif", format="GTiff", progress="text")
endCluster()
如何设置以上代码?如果我仅使用“插值”,则要花一个月的时间才能得出结果。