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

3

我希望创建一个由空间数据点构成的插值网格,以便在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处理空间数据方面还很陌生。

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)

成功了!只需要稍微调整一下数值就可以了。谢谢。 - Tricia

网页内容由stack overflow 提供, 点击上面的
可以查看英文原文,
原文链接