Mi primer cuaderno: Colombia y el departamento de Bolívar

El siguiente proyecto tiene como objetivo realizar un mapa de Colombia identificando las coordenadas de sus puntos extremos, sus 32 departamentos y sus municipios. Sin embargo este proyecto va a estar enfocado en el departamento de Bolívar, conociendo así sus coordenadas, departamentos vecinos y sus municipios. Por lo tanto antes de empezar con el mapa de Colombia se van a tener unos conocimientos previos según un ejemplo en Rstudio.

Ejemplo de RStudio:

plot(cars)

Antes de comenzar con la realización de mapas, es necesario instalar previamente ciertas librerías, como lo son tidyverse, sf, GADMTools y mapview.

Llamado de librerías

library(tidyverse)
library(sf)
library(GADMTools)
library(mapview)

Un ejemplo de un mapa de Carolina del Norte, donde se muestran sus 100 condados y el número de nacimientos en cada uno en el año 1974.

system.file("gpkg/nc.gpkg", package = "sf") %>%
  read_sf() %>%
  st_transform(32119) %>%
  select(BIR74) %>%
  plot(graticule = TRUE, axes = TRUE)

A continuación una tabla donde se hace una interpretación del mapa anterior, los datos que están ahí.

(read_sf(system.file("gpkg/nc.gpkg", package = "sf")))

Ahora si, entrando en materia se van a obtener los datos administrativos del país Colombia, donde se ven las coordenadas de los estremos del país.

Como ya se instaló la librería GADMTools, se va a cargar la libreía.

library(GADMTools)
COL <- gadm_sf_loadCountries(c("COL"), level=0, basefile="./")
COL

Con los datos obtenidos:

Ahora se va a plotear a Colombia a nivel 0

gadm_plot(COL)

Posteriormente se van a extraer departamentos (i.e. level 1)

DEPTOS <- gadm_sf_loadCountries(c("COL"), level = 1, basefile = "./")
gadm_plot(DEPTOS)

Y por último se van a plotear sus municipios (nivel 2)

MUNICIPIOS <- gadm_sf_loadCountries(c("COL"), level = 2, basefile = "./")
gadm_plot(MUNICIPIOS)

Como se comentó anteriormente, ahora se va a hacer énfasis en el departamento de Bolívar.

Ahora se va a plotear al departamento Bolívar

Bolívar <- gadm_subset(DEPTOS, level=1, regions="Bolívar")
gadm_plot(Bolívar) %>% gadm_showNorth("tl") %>% gadm_showScale("bl")

Antes de continuar con el departamento, se va a obtener una lista de los departamentos y municipios de Colombia.

Lista de los departamentos

listNames(DEPTOS, 1)

Lista de los municipios de Colombia

listNames(MUNICIPIOS, 2)

Continuando con el departamento de Bolívar:

Recortar una área (departamento de Bolívar)

Se va a encerrar con una caja el departamento, para que se muestran sus límites.

(sbox <- gadm_getBbox(Bolívar))

Estructura de la caja que encierra al departamento de Bolívar

str(sbox)

Ahora se va a seguir estudiando el departamento Bolívar recortándolo en unos lítes específicos (más grandes a los de la caja)

(STUDY_AREA <- gadm_crop(DEPTOS, xmin=-75.8, ymin=6.98, xmax=-73.77, ymax=10.8))
                        

Departamentos vecinos a Bolívar

gadm_plot(STUDY_AREA)

Obtención de un simple future

DEPTOS$sf

Como ya se instaló la librería mapview, se va a llamar la librería

Llamado de librería

install.packages("mapview")

Este código es para suprimir mensajes de Warning al momemto de llamar la librería mapview y posteriormente se va a realizar un mapa

suppressPackageStartupMessages(library(mapview))
DEPTOS$sf %>% mapview(zcol = "NAME_1", legend = TRUE, col.regions = sf.colors)

Por último se va a plotear los municipios del departamento Bolívar

muncB <- gadm_subset(MUNICIPIOS, 1, "Bolívar", usevar = NULL)
gadm_plot(muncB)

muncB$sf %>% mapview(zcol = "NAME_2", legend = TRUE, col.regions = sf.colors)

Ahora se van a plotear los municipios de Bolívar

Se toma la variable de los MUNICIPIOS indicando que me voy a centrar en los municipios, luego se coloca el numero 1 y el departamento escogido. se van a extraer los municipios de Bolívar

BOL_mun <- gadm_subset(MUNICIPIOS, 1, "Bolívar", usevar = NULL)

BOL_mun$sf %>% mapview(zcol = "NAME_2", legend = TRUE, col.regions = sf.colors)

class(BOL_mun$sf)

[1] "sf"         "data.frame"