r - 合并 R 中的相邻区域(聚合空间数据)？

Question

我想我需要改写我之前措辞糟糕的问题(删除它)。这是另一个尝试。我想加入相邻的区域，以使它们的共同边界消失并且只能看到它们的外线。

这是一个可重现的示例:

require(shapefiles)
require(sp)

xx <- readShapeSpatial(system.file("shapes/sids.shp", package="maptools")[1],
                   IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))

# show all the subregions
plot(xx)

现在让我们只考虑区域 3 和 5

plot(xx[c(3,5),])

我怎样才能聚合这些区域。在实践中，我想做的就像拥有一张显示所有国家的整个大陆的 map ，并制作一张显示北美和南美的 map 。

对我来说，这看起来像是一项非常常见的任务，但到目前为止我找不到合适的功能来完成它。我只是错过了一个功能，还是可以简单地手动完成？

Answer 1

rgeos软件包提供了许多出色的处理工具Spatial*数据，可以在这种情况下使用。

例如:

library(rgeos)
regionOfInterest <- gUnion(xx[3,], xx[5,])

这也有相同的结果，并且可能对多个多边形更有用:

regionOfInterest <- gUnionCascaded(xx[c(3,5), ])

来自 plot(regionOfInterest) 的结果: