【发布时间】:2015-03-27 22:07:28
【问题描述】:
我有一些想要合并的相邻宗地集群,但我希望保留每个集群内的宗地边界。如果我使用 Reduce(geom_list, gUnion) 之类的东西,那么相邻的边界就会消失。
这是我想要的一个简单且可重复的示例,或多或少在纯 R 中实现:https://gist.github.com/holstius/37c0ff01109e28ff817a。如果GEOS 还没有本机实现它,我会感到惊讶,但我无法确定可以满足我需求的rgeos 绑定。
rgeos::gUnion 的变体是否不会分解相邻的多边形?
(编辑:删除了%>% 以支持https://gist.github.com/holstius/37c0ff01109e28ff817a 的粗略代码示例)
【问题讨论】:
-
我不清楚“合并”是什么意思,但又不想解散。因为我不熟悉
.和%>%语法,所以您的示例草图对我没有帮助,并且缺少可重现的示例。 -
感谢您的参与。请给我几分钟时间来制作一个仅使用基本 R 和
sp和/或maptools的可重现示例。希望这将阐明意图。 -
@edzer-pebesma 请参阅gist.github.com/holstius/37c0ff01109e28ff817a 以获得简单且可重复的示例。我想使所有
Polygons唯一的IDs 的逻辑可以在sp::rbind.SpatialPolygons(..., makeUniqueIDs = TRUE)中实现。但这个想法只是为了表明几何图形是在不溶解相邻边界的情况下合并的。 -
我还是看不出你的
merge比sp 的rbind有什么改进;我也不明白你为什么称它为merge,因为base::merge和它在sp 中的对应物。