我正在尝试使用geoR软件包创建不规则网格以用于克里金法。有问题的地区是马拉维的国家。数据目前在UTM投影中：
data@bbox
min max
x 464479.7 814025.4
y -1894837.5 -1035058.0
我正在创建矩形网格，如下所示：
x<-seq(460000, 820000, l=51)
y<-seq(-1900000, -1000000, l=51)
grid<-expand.grid(x=x, y=y)
尝试1创建一个不规则网格，如下所示：
outline <- data@polygons[[2]]@Polygons[[1]]@coords
library(splancs)
new_grd <- grid[inout(grid,outline), ]
我没有收到任何错误消息或警告，但这是输出：
head(new_grd)
[1] x y
<0 rows> (or 0-length row.names)
任何关于可能发生的事情的想法都将受到欢迎！
谢谢， 朱莉安（Julianne）