R中溶解多边形特征的边界

Question

我一直在寻找现有的 R 函数，用于聚合同一图层内共享公共边界的多边形要素（即，生成类似于 ArcGIS 中的“溶解边界”工具的输出）。

我使用 gdal_polygonizeR （ https://johnbaumgartner.wordpress.com/2012/07/26/getting-rasters-into-shape-from-r/ ）从光栅文件创建了一个多边形图层。一些多边形形状由单个栅格单元分隔，因此在形状文件中存储为不同的特征。我想将此类多边形特征组合成单个多边形特征，并创建一个新的形状文件（减少多边形元素的总数），最好具有溶解的阈值距离。

有谁知道 R 中现有的方法可以做到这一点？

更新：我认为解决方案可能涉及aggregate，然后是disaggregate。我目前正在探索这一点，特别注意确保带有孔的多边形要素与父多边形保持关联（请参阅：分割单个 SpatialPolygons 对象的多边形部分）。如果/当我找到解决方案时将再次更新。

Answer 1

library(raster)
y <- aggregate(x)