snitr: un paquete de R para descargar datos geoespaciales del SNIT (Sistema Nacional de Información Territorial), Costa Rica

Descripción del paquete

snitr es un paquete de R que permite bajar datos geospaciales del Sistema Nacional de Información Territorial SNIT que tiene como objetivos generales promover y articular la generación de productos, servicios e información geográfica georreferenciada de cubrimiento nacional, regional y local de Costa Rica.

Instalación de snitr

Primero se debe instalar el paquete devtools

install.packages("devtools")

Luego proceder a instalar snitr

devtools::install_github("ManuelSpinola/snitr")

Ejemplo de uso de snitr

Cargar los paquetes o librerías snitr, tidyverse y sf

library(snitr)
library(tidyverse)
library(sf)

Listar las capas disponible de escala 1:5000

capas <- list_IGN_5k()

 [1] "IGN_5:forestal2017_5k"      "IGN_5:cultivos2017_5k"     
 [3] "IGN_5:curvas_5000"          "IGN_5:delimitacion2017_5k" 
 [5] "IGN_5:edificaciones2017_5k" "IGN_5:hidrografia_5000"    
 [7] "IGN_5:linea_costa_5000"     "IGN_5:pastos2017_5k"       
 [9] "IGN_5:urbano_5000"          "IGN_5:vias_5000"           

Escoger la capa que se quiere descargar

layer_name <- "IGN_5:delimitacion2017_5k"

Asignarle un nombre al archivo, pude ser en formato shapefile o geopackage

output_file <- "IGN_5_delimitacion2017_5k.gpkg"

Ejecutar la función get_IGN_5K para descargar la capa

get_IGN_5k(layer_name, output_file = output_file)

Reading layer `file10a8768080076' from data source 
  `/private/var/folders/96/m4trt68j16v3n97qk3zkxrk40000gn/T/RtmpzineWI/file10a8768080076.geojson' 
  using driver `GeoJSON'
Simple feature collection with 1 feature and 6 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: 156147.7 ymin: 608897.3 xmax: 658879.6 ymax: 1241118
Projected CRS: CR-SIRGAS / CRTM05
Deleting source `IGN_5_delimitacion2017_5k.gpkg' failed
Writing layer `IGN_5_delimitacion2017_5k' to data source 
  `IGN_5_delimitacion2017_5k.gpkg' using driver `GPKG'
Writing 1 features with 6 fields and geometry type Multi Polygon.

[1] "IGN_5_delimitacion2017_5k.gpkg"

Leer la capa de donde se haya guardado

cr <- read_sf("IGN_5_delimitacion2017_5k.gpkg")

Realizar un mapa con la capa descargada

ggplot() +
  theme_minimal() +
  geom_sf(data = cr, fill = "dodgerblue3")