r gis：找到多边形之间的边界

Question

有了一个多边形shapefile，我需要生成一个折线shapefile，它只包含多边形之间的公共边框（见图）。

我的问题类似于1和2 ，只有我需要在R执行此操作。 后一个类似的问题是指使用Shapely包for python的解决方案。 Shapely for R的类似物是rgeos 。 虽然，我自己找不到rgeos的解决方案。

---

注意：带有用于插图的边框的shapefile是使用类似问题1的解决方案在ArcGIS中生成的。现在我需要在R执行相同的操作。

Answer 1

你想要的是从溶解区域的线组和区域本身的线之间的差异线。 在rgeos包中， gUnaryUnion将解散多边形， gDifference将减去。

对于我的小欧盟子集eusub ，我可以这样做：

library(rgeos); library(sp)
borders = gDifference(
   as(eusub,"SpatialLines"),
   as(gUnaryUnion(eusub),"SpatialLines"),
   byid=TRUE)

注意需要将多边形转换为线，因为输出将是线。

然后看到这个：

plot(eusub)
plot(borders, col="red",lwd=2,add=TRUE)