出于某种原因,我无法让@RichPauloo 提供的解决方案正常工作,并感谢一些帮助。
我有一个名为“spdf”的 SpatialPolygonsDataFrame(在下面的保管箱链接中) https://www.dropbox.com/s/ibhp5mqbgfmmntz/spdf.Rda?dl=0
我使用下面帖子中的代码来获取边界内的网格数据。
library(sp)
grd <- makegrid(spdf, n = 10000)
colnames(grd) <- c('x','y');
outline <- spdf@polygons[[1]]@Polygons[[1]]@coords
library(splancs)
new_grd <- grd[inout(grd,outline), ]
这是我得到的:
黑点是来自 makegrid 的“grd”
蓝点以“轮廓”为界
红点是“new-grd”作为边界内的网格
如您所见,它没有捕获边界内的所有数据?我究竟做错了什么?