地理信息系统 : find borders between polygons

Question

有了一个多边形 shapefile，我需要生成一个仅包含多边形之间公共(public)边界的折线 shapefile(见图)。

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

---

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

Answer 1

您想要的是溶解区域的线集与区域本身的线之间的差异线。在 rgeos 包中，gUnaryUnion 将分解多边形，而 gDifference 将相减。

对于我的小型 EU 子集 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)