我希望创建一个由空间数据点构成的插值网格,以便在R中进行Kriging。我一直在使用sp和gstat包进行工作。我正在使用的代码可以正常运行,但我的网格不包括所有点,我想扩展它,但我似乎无法正确操作。
我正在使用的代码如下所示(我的数据框是kr.data):
x.range <- as.integer(range(kr.data@coords[,1]))
y.range <- as.integer(range(kr.data@coords[,2]))
data.grd <- expand.grid(x=seq(from=x.range[1], to=x.range[2], by=0.5),
y=seq(from=y.range[1], to=y.range[2], by=0.5))
coordinates(data.grd) <- ~x+y
gridded(data.grd) <- TRUE
plot(data.grd, cex=0.5)
points(kr.data, pch=1, col='red', cex=0.7)
我希望添加输出,以显示聚集在网格之外的右下角的点,但我是新手,不确定如何添加。是否有人可以提供建议,告诉我如何使用空间数据制作网格或更改我的代码?当我修改我的代码以扩展网格时,我的数据点不再可见。我在使用R处理空间数据方面还很陌生。