如何将shapefile中的多边形要素与多个多边形合并？
rbind和union只是组合了shapefile特征的行，它们实际上并没有合并多边形。
以下示例中的所需结果：
如何将以下带有重复ID_2的shapefile合并到sptemp中的单个多边形？
以下埃塞俄比亚GADM级别2的示例具有重复的shapefile ID_2列的前两行（值= 1）。我喜欢sptemp有79个功能，前两行是重复ID_2的行。 sptemp [1，]的图将显示当前sptemp [1，]和sptemp2 [2，]没有重复之间的边界，即多边形也合并。
示例代码：
下载，解压缩并加载到埃塞俄比亚2级的R GADM文件（899kb到工作目录）：
library(curl)
library(rgdal)
curl_download("http://biogeo.ucdavis.edu/data/gadm2.8/shp/ETH_adm_shp.zip",
destfile=paste0("gadmETH.zip"),
quiet=FALSE)
unzip(paste0("gadmETH.zip"), exdir="gadmETH", overwrite=FALSE)
###Load shapefile
sptemp <- readOGR(dsn="gadmETH", layer="ETH_adm2")
前两个多边形的ID_2列重复
###You'll see in the first two rows ID_2 is duplicated
df.sptemp <- as.data.frame(sptemp)
View(sptemp)
###I can't just delete one because they are separate polygons
plot(sptemp[1,], col="blue")
plot(sptemp[2,], add=TRUE, col="red" )
答案 0 :(得分：3)
注意此方法使用st_union，它结合了所有多边形＆＃39;成单个多边形。这可能不是您真正想要的结果。
如果您使用library(sf)而不是sp（它是sp的后继者），则可以使用st_union来连接几何。
您也可以在dplyr管道序列中执行此操作。
library(sf)
sptemp <- sf::st_read(dsn = "~/Desktop/gadmETH/", layer = "ETH_adm2")
library(dplyr)
sptemp %>%
group_by(ID_2) %>%
summarise(geometry = sf::st_union(geometry)) %>%
ungroup()
# Simple feature collection with 79 features and 1 field
# geometry type: GEOMETRY
# dimension: XY
# bbox: xmin: 33.00154 ymin: 3.398823 xmax: 47.95823 ymax: 14.84548
# epsg (SRID): 4326
# proj4string: +proj=longlat +datum=WGS84 +no_defs
# # A tibble: 79 x 2
# ID_2 geometry
# <dbl> <sf_geometry [degree]>
# 1 1. POLYGON ((38.85556 8.938293...
# 2 2. POLYGON ((42.15579 12.72123...
# 3 3. POLYGON ((40.17299 14.49028...
# 4 4. POLYGON ((41.11739 10.93207...
# 5 5. POLYGON ((40.61546 12.7958,...
# 6 6. POLYGON ((40.25209 11.24655...
# 7 7. POLYGON ((36.35452 12.04507...
# 8 8. POLYGON ((40.11263 10.87277...
# 9 9. POLYGON ((37.39736 11.60206...
# 10 10. POLYGON ((38.48427 12.32812...
# # ... with 69 more rows