将多个多边形合并到一个有边界的多边形中答案

【问题标题】：merge multiple polygon in one bordered polygon将多个多边形合并到一个有边界的多边形中
【发布时间】：2021-06-18 17:30:59
【问题描述】：

我正在尝试将哥伦比亚的两个多边形（部门）合并为一个多边形。我的目标是只保留两个多边形的边缘，但这不会发生。下面我附上代码。感谢您的帮助。

download.file("https://github.com/nebulae-co/colmaps/raw/master/data/departamentos.rda","departamentos")
load("departamentos")

library(sf)
library(dplyr)
library(rgdal)
library(mapview)

df1 = subset(x = departamentos, depto == "Bogotá, D. C.")
df2 = subset(x = departamentos, depto == "Cundinamarca")

df1@data$id <- df2@data$id <- "a"

df3 = rbind(df1, df2, makeUniqueIDs = TRUE)
plot(df3)

绘制两个多边形

## option 1
df_12 <- sf::st_as_sf(df3) %>%
  group_by(id) %>%
  summarise(geometry = sf::st_union(geometry)) %>%
  ungroup()

plot(as_Spatial(df_12))

## option 2
plot(raster::aggregate(df3))

## option 3
plot(rgeos::gUnaryUnion(df3))

## optiom 4
plot(st_union(st_make_valid(st_as_sf(df3))))

## option 5
plot(st_as_sf(rgeos::gBuffer(as_Spatial(st_make_valid(st_as_sf(df3))), byid=F, width=0)))

我的预期输出

【问题讨论】：

标签： r polygon spatial sp

【解决方案1】：

显然问题与我使用的第一个形状有关，寻找另一个从 Arcgis 下载的形状，它可以正常工作。

Arcgis 网页

https://hub.arcgis.com/datasets/juanespo::departamentos-1?geometry=-107.669%2C-3.045%2C-40.917%2C12.241

Shapefile .zip 直接下载

https://opendata.arcgis.com/datasets/2059399a41da4dc0b7683d1822194ae0_0.zip

library(sf)
library(mapview)
library(dplyr)

departamentos <- st_read("E:/Users/Rafae/Downloads/departamentos/b3c7aaf2-9439-4bb6-92e1-aa7b3c401b20202046-1-1gl30yi.xiu1.shp")
#mapview(departamentos)

Centro <- subset(x = departamentos, NOMBRE_DPT %in% c("BOGOTÁ","CUNDINAMARCA")) %>%
  st_as_sf() %>% st_make_valid() %>% st_union()

plot(Centro)

【讨论】：