Question

我想创建一个由空间数据点组成的插值网格，以便以后用于R中的克里金法。我一直在使用sp和gstat包。我正在使用的代码有效，但我的网格不包含所有的点，我想扩展它，但我似乎无法正确使用它。

我正在使用的代码在这里（我的data.frame是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中空间数据的新手。

Answer 1

由as.integer整数产生的截断可能会切断范围高端的值。请尝试改为：

x.range <- as.integer(range(kr.data@coords[,1])) + c(0,1)
y.range <- as.integer(range(kr.data@coords[,2])) + c(0,1)

在R中扩展我的插值网格以进行克里金法

1 个答案: