Leer un shapefile
Los shapefiles son una de las formas más comunes de compartir datos espaciales y se leen fácilmente en R usando readOGR() del paquete rgdal. readOGR() tiene dos argumentos importantes: dsn y layer. Lo que pases exactamente a estos argumentos depende del tipo de datos que estés leyendo. En el vídeo aprendiste que, para shapefiles, dsn debe ser la ruta al directorio que contiene los archivos que componen el shapefile y layer es el nombre del archivo concreto del shapefile (sin extensión).
Para tu mapa, quieres los límites de los barrios. Hemos descargado las Neighborhood Tabulation Areas, tal como las define la Ciudad de Nueva York, desde la Open Data Platform of the Department of City Planning. La descarga venía en un archivo zip y hemos colocado el resultado de descomprimirlo en tu directorio de trabajo.
Usarás la función dir() de R base para examinar el contenido de tu directorio de trabajo y luego leer el shapefile en R.
Este ejercicio forma parte del curso
Visualización de datos geoespaciales en R
Instrucciones del ejercicio
- Usa
dir()sin argumentos para averiguar el nombre del directorio del shapefile. - Usa
dir(), pasando la ruta al directorio del shapefile, para ver los archivos que contiene. - Ahora ya conoces el directorio y el nombre del archivo. Usa
readOGR()para leer el shapefile de barrios en un objeto llamadoneighborhoods. - Comprueba el contenido llamando a
summary()sobreneighborhoods. - Comprueba el contenido representando
neighborhoodscon un gráfico.
ejercicio interactivo práctico
Prueba este ejercicio completando este código de ejemplo.
library(sp)
library(rgdal)
# Use dir() to find directory name
# Call dir() with directory name
# Read in shapefile with readOGR(): neighborhoods
# summary() of neighborhoods
# Plot neighborhoods