library(sf)
library(geobr)
library(tidyverse)
Estou tendo dificuldades em usar o st_intersection()
do pacote sf
para determinados dados carregados com o pacote geobr
.
Por exemplo, se tento fazer uma interseção entre o sf de um estado e um sf das Unidades de Conservação do país, não consigo:
ucs <- read_conservation_units()
paraiba <- read_state(code_state = "PB")
Se tento fazer a interseção, surge a seguinte mensagem:
st_intersection(ucs, paraiba)
Error in s2_geography_from_wkb(x, oriented = oriented, check = check) :
Evaluation error: Found 5 features with invalid spherical geometry.
[402] Loop 0 is not valid: Edge 3 is degenerate (duplicate vertex)
[820] Loop 0 is not valid: Edge 3 is degenerate (duplicate vertex)
[1394] Loop 0 is not valid: Edge 3 is degenerate (duplicate vertex)
[1429] Loop 0 is not valid: Edge 3 is degenerate (duplicate vertex)
[1555] Loop 0 is not valid: Edge 3 is degenerate (duplicate vertex).
Por outro lado, se utilizo a mesma lógica para fazer a interseção entre os biomas, a função st_intersection()
funciona:
biomas <- read_biomes()
biomas_paraiba <- st_intersection(biomas, paraiba)
Verificando o biomas_paraiba
:
biomas_paraiba
Simple feature collection with 3 features and 8 fields
Geometry type: MULTIPOLYGON
Dimension: XY
Bounding box: xmin: -38.7656 ymin: -8.302955 xmax: -34.79288 ymax: -6.026569
Geodetic CRS: SIRGAS 2000
name_biome code_biome year code_state abbrev_state name_state code_region name_region geom
2 Caatinga 2 2019 25 PB Paraíba 2 Nordeste MULTIPOLYGON (((-35.22321 -...
4 Mata Atlântica 4 2019 25 PB Paraíba 2 Nordeste MULTIPOLYGON (((-35.23148 -...
7 Sistema Costeiro NA 2019 25 PB Paraíba 2 Nordeste MULTIPOLYGON (((-35.00158 -...
E o plot:
ggplot()+
geom_sf(data = biomas_paraiba, aes(fill = name_biome))
Em suma: porque consigo fazer a interseção com os biomas, mas não consigo com as unidades de conservação?
sf::st_make_valid
antes da intersecção.sf::st_make_valid
. Obrigado @CarlosEduardoLagosta