我有来自美国 census.gov 的人口普查数据的 Shapefile。
我已经使用 rdgal 库中的 readOGR 将文件读入 R。
sp1 : 投票区 - SpatialPolygonsDataFrame
sp2：比 sp1 更大区域和更高分辨率的块级别 - SpatialPolygonsDataFrame
我曾尝试使用 rgeos 库中的 gIntersection，但它没有返回我需要的内容。或者我不知道如何正确使用它
如何使用 sp1 作为空间限制从 sp2 中提取较低分辨率的多边形和相关数据？换句话说，我想知道 sp1 中 sp1 和 sp2 相交的所有块。
当前工作解决方案
下面的代码返回 sp2 的索引值向量，其中多边形在 sp1 内至少有一个点，但我想知道是否有更简单的方法。
library(rgeos)
result <- vector()
for (i in 1:nrow(sp2)) {
if(gContains(sp1, sp2[i,])) {
result <- c(result, i)
}
}