3

我想创建一个由空间数据点组成的插值网格,以便以后在 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 中处理空间数据的新手。

4

1 回答 1

3

整数导致的截断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)
于 2012-10-20T18:51:44.103 回答