r - 寻找允许选择/识别与源多边形共享线段的多边形的功能

Question

为了选择/识别 shapefile 的边界多边形，我想使用一个能够选择/识别与源多边形共享线段的多边形的函数。

有数字:

我有这种形状文件:

使用 rgeos 包中的 gUnionCascaded，我有第二个带有“轮廓多边形”的 shapefile

现在我正在寻找一个可以选择/识别边界多边形(图中阴影部分)的函数，即第一个 shapefile 的多边形与第二个 shapefile 的多边形共享一条线段。 :

Answer 1

根据 Josh O'Brien 的建议，我使用了 rgeos::gRelate() 函数。我得到 3 个 DE-9IM 箱子:

x <- gRelate(shapefile.1, shapefile.2, byid = TRUE)
table(x)
2FF10F212 2FF11F212 2FF1FF212 
       63      2470    174495        

得到的 DE-9IM 字符串代码可以解释如下:

1) 2FF1FF212:表示第一个 shapefile 中的多边形，这些多边形不与第二个 shapefile 的多边形边界相交

2) 2FF11F212:表示第一个 shapefile 中的多边形与第二个 shapefile 的多边形边界相交

3) 2FF10F212:表示第一个 shapefile 中的多边形与第二个 shapefile 的多边形边界与一个点相交

最后两个案例是我一直在寻找的边界多边形。我有他们的 ID:

poly.border <- which(x %in% c("2FF10F212","2FF11F212"))