如何在R中使用自克里格将通用克里格与自定义预测空间网格相结合

我想使用R中的autokrige函数对数据集应用通用克里格法。我想为预测点创建自己的自定义空间网格（用于autokrige的新_数据参数）。我使用的是R版本3.2.2（64位）和RStudio版本0.99.486。以下是我迄今为止所做的工作：

library(automap)
library(sp)
library(gstat)
library(raster)
library(rgdal)

data(meuse)
coordinates(meuse) <- ~x + y
proj4string(meuse) <- CRS("+init=epsg:28992")

然后接收到以下错误：

model.frame.default中出错（terms.f，newdata，na.action=na.action，： 对象不是矩阵。此外：警告消息： “newdata”有3102行，但找到的变量有1行

---

我在创建网格时出错了吗（r_meuse）？有没有“更好”的方法来为预测数据创建网格？到目前为止，我找到的所有示例都使用meuse.grid数据，但我想将通用克里格法应用于其他还没有自己的网格数据的数据。

我相信这里的问题是，您在执行UK时没有预测器，

dist

t in

r\u-meuse

。这是一个问题，因为线性预测需要信息。因此，

r\u-meuse

需要是一个定义了

dist

的

空间点数据帧

ext_meuse <- as(extent(meuse), "SpatialPolygons")
r_meuse <- rasterToPoints(raster(ext_meuse, resolution = 59), spatial = TRUE)
proj4string(r_meuse) <- proj4string(meuse)

kriging_result = autoKrige(zinc~dist, meuse, r_meuse)