将R内的GIS区域统一起来（例如，将县划分为基于核心的统计区域）

我有一个美国各州的形状文件，我用readOGR读入r。我需要将3000个美国县合并成更大的区域，例如，929个基于核心的统计区域提供了从县FIP到it的人行横道。在平面设计中，这将是形状的结合。我是否可以在R中以编程方式执行此操作而不使用地理信息系统？

使用函数ms\u dissole from。例如，要将县的形状合并为州：

library(rmapshaper)

# Download the shapefile for US counties from here and save in some/dir
# www2.census.gov/geo/tiger/GENZ2010/gz_2010_us_050_00_5m.zip
counties <- readOGR(dsn="some/dir",layer="gz_2010_us_050_00_5m")

# Drop Alaska, hawaii, Puerto Rico
counties <- counties[!(counties$STATE %in% c("02","15","72")), ]

# Plot counties
plot(counties)

# Unite shapes
library(rmapshaper)
states <- ms_dissolve(counties, field = "STATE")

plot(states)

县地图：

县联盟的州地图：

---

sf软件包中的函数可能就是您想要的。sf“简单要素”使处理空间数据比使用QGIS或ArcGis要容易得多。

您的链接说它可以用于一个多边形作为输出：将多个要素几何图形合并为一个，而无需合并或解析内部边界。我想要多个多边形作为输出，所以我没有尝试。这很有效，是的！我可以编辑您的答案以提供完整的代码示例吗？太好了！是的，没问题！