Prueba de funcionalidades espaciales

Introducción

Este cuaderno es una prueba de algunas funcionalidades de la libreria sf (simple features) de R.

Configuración

Antes de todo, hay que instalar las librerias que vamos a usar. Esto se debe hacer desde la consola, no desde un cuaderno por que si no se instala n veces.

#install.packages("sf")
#install.packages("dplyr")
#install.packages("readr")

Vamos a cargar las librerias que necesitamos:

library(sf)
library(readr)
library(dplyr)

Vamos a definir el directorio de trabajo:

setwd("/Users/felipe/Documents/GB2/Proyecto1")

getwd()

list.files("datos")

Lectura de datos

municipios <- st_read("/Users/felipe/Documents/GB2/Proyecto1/Datos/MGN2018_MPIO")

vias <- st_read("/Users/felipe/Documents/GB2/Proyecto1/Datos/COL_rds")

(nuevo_crs <- st_crs(municipios))

vias_4686 <- st_transform(vias, 4686)

Seleccionar municipios de un departamento en esta caso Risaralda

mun_risaralda <- municipios %>% filter(DPTO_CCDGO == "66")

Recortar las vias de Risaralda

vias_risaralda <- st_intersection(vias_4686, mun_risaralda)

mun_risaralda2 <- mun_risaralda %>% select(DPTO_CCDGO, MPIO_CCDGO,MPIO_CNMBR, AREA)

head(mun_risaralda2)