在 R (CRAN) 中创建和裁剪网格

我需要在 gstat 库中创建一个用于空间插值的掩模网格。具体来说，我有随机分布的不同采样点，我需要创建包围这些点的最小凸多边形。然后，我必须创建一个空间网格，该网格应由计算的外壳裁剪，以将插值限制在该多边形的范围内。如果有人能向我解释详细的过程并提供一些示例，我将不胜感激。先感谢您。

我自己找到了解决方案。

library(spatstat)
library(sp)
library(plotKML)
library(maptools)

创建随机点

x<-rnorm(100,3)
y<-rnorm(100,3)
plot(x,y)
xy<-cbind(x,y)
xy<-as.data.frame(xy)

将点转换为空间点数据帧，然后转换为栅格。

coordinates(xy)=c("x","y")
pnts<-vect2rast(xy)
summary(pnts)

检查汇总像元大小值并记住它

从点创建凸包。 然后，将“owin”对象（凸包的类）转换为空间多边形（栅格创建的基本步骤）

conv<-convexhull.xy(x,y)
SpP<-as(conv,  "SpatialPolygons")
plot(SpP)
points(x,y)
attr  =  data.frame(a=1,  b=1)
SrDf  =  SpatialPolygonsDataFrame(SpP,  attr)

将“cell.size”设置为与“summary(pnts)”相同（在本例中设置为 0.085）。

rast <- vect2rast(SrDf,cell.size=0.085)

玩得开心！

plot(rast)
image(rast)
points(x,y)

笔记： 使用 vect2rast，如果没有为“rast”对象设置 cell.size，该函数会自动计算基于 un 点密度分布的最合适的像元大小。在这种情况下，多边形仅由其顶点定义，因此我们使用为我们想象的多边形包围的点计算的像元大小。