从包自动化中的autoKrige预测

Question

我正在尝试预测由autoKrige计算到new_data的模型。

就我而言，new_data是栅格的单元格，以R包raster计算。

autoKrige无法预测我的栅格：

library(automap)
afvg1 <- autoKrige(Sand~x+y, sand2utmOGR, 
                   model = c("Sph", "Exp", "Gau", "Ste"), new_data=rasta2)

Error in (function (classes, fdef, mtable)  :
unable to find an inherited method for function ‘is.projected’ for signature "RasterLayer"’

如何将栅格转换为autoKrige可用于预测的格式？

Answer 1

猜测...对于某些栅格r，请设为SpatialPoints，将其用于new_data，将预测反馈回栅格。

nd = as(r,"SpatialPoints")
fvg1<-autoKrige(blah blah, new_data=nd)
r[]=(predictions from fvg1) 

当你的问题出现时，我会把我的答案剪切成糊状。