r - 提取接触r中指定多边形的多边形

Question

我正在使用包含 r 中的多边形的 shapefile。我希望在 shapefile 以及周围的多边形中提取特定的多边形。即直接接触我指定多边形的多边形。

这在r中可能吗？我自己在另一个 GIS 软件中检查了数据，所以我可以看到这些周围多边形的 id 并手动提取它们

                   shape_old<-readOGR(dsn=shapefile.shp)
          shape_new<-subset(shape_old, town=="Pleasantville"

- 但我希望有一种自动化的方式来做到这一点。

我已经阅读了有关 gtouches() 的内容，但据我了解，这是一个测试多边形是否接触的测试，而不是提取它们的方法。

标签： rshapefile