União de feature

dúvidas
1

Olá

Tenho um arquivo do tipo shapefile de Portugal. Há vários níveis de divisão geográfica, sendo a menor as FREGUESIAS. Preciso “dissolver”, em regiões maiores, os DISTRITOS. Qual a sugestão para resolver isso? O código base é:

library(tidyverse)
library(sf)
library(ggplot2)
library(httr)

# Lendo a base de dados de dados cartográficas de Portugal

httr("http://mapas.dgterritorio.pt/wms-inspire/caop/continente?service=WMS&REQUEST=GetCapabilities&VERSION=1.3.0")

portugal <- st_read("Cont_AAD_CAOP202/Cont_AAD_CAOP2021.shp",
                      quiet = T)

glimpse(portugal)

Rows: 3,223
Columns: 9
$ Dicofre    <chr> "081504", "081504", ~
$ Freguesia  <chr> "Sagres", "Sagres", ~
$ Concelho   <chr> "Vila do Bispo", "Vi~
$ Distrito   <chr> "Faro", "Faro", "Far~
$ TAA        <chr> "ÁREA SECUNDÁRIA", "~
$ Area_T_ha  <dbl> 3436.93, 3436.93, 34~
$ Area_EA_ha <dbl> 0.04, 0.17, 0.20, 0.~
$ Des_Simpli <chr> "Sagres", "Sagres", ~
$ geometry   <POLYGON [m]> POLYGON ((-6~

ggplot(portugal)+
  geom_sf(aes())
```
2

Hidelbrando,

Como seu código não é reprodutível (a função httr() não existe e o arquivo linkado parece ser um XML), eu não tenho como testar a minha ideia. De qualquer forma, eu acho que você precisa usar a st_union() ou fazer um group_by() |> summarise(). Veja os links para exemplos do que eu estou falando.