Maps

Importar librerias

library(ggplot2)
library(rio)
library(sf)

## Linking to GEOS 3.11.2, GDAL 3.8.2, PROJ 9.3.1; sf_use_s2() is TRUE

Graficando los departamentos según su grado de acumulación de obras paralizadas

linkdepartamento="https://github.com/Pepix2345/introgeodf/raw/main/maps/Perumap.gpkg"

sf:: st_layers(linkdepartamento)

## Driver: GPKG 
## Available layers:
##    layer_name geometry_type features fields crs_name
## 1 Peruborders       Polygon        1      1   WGS 84
## 2  Perupoints         Point        8      3   WGS 84
## 3   Perulines                      5      2   WGS 84

peru=read_sf(linkdepartamento, layer="Peruborders")

linkprovincias="https://github.com/bocchienjoyer/MaGODllanes/raw/main/ProyectoMagaDeidad/datita/peru_prov.gpkg"

sf:: st_layers(linkprovincias)

## Driver: GPKG 
## Available layers:
##   layer_name geometry_type features fields crs_name
## 1     riesgo Multi Polygon      196      3   WGS 84

provincias=read_sf(linkprovincias, layer="riesgo")

baseprovincia= ggplot(data= peru) + geom_sf(fill= "grey90") + theme_light()
basefinalrepro= baseprovincia + geom_sf(data= provincias, aes(color=escala))
basefinalrepro

linkdepa="https://github.com/bocchienjoyer/MaGODllanes/raw/main/ProyectoMagaDeidad/datita/dep_map.gpkg"

sf:: st_layers(linkdepa)

## Driver: GPKG 
## Available layers:
##     layer_name geometry_type features fields crs_name
## 1 modalidad_xd Multi Polygon       26      4   WGS 84

departamento=read_sf(linkdepa, layer="modalidad_xd")

basedepartamento= ggplot(data= peru) + geom_sf(fill= "grey90") + theme_light()
basefinaldepa= basedepartamento + geom_sf(data= departamento, aes(color=modalidad))
basefinaldepa