Lectura de datos espaciales

Visualización de datos geoespaciales en R

Charlotte Wickham

Assistant Professor at Oregon State University

Rentas medianas en el condado de Nueva York

Los distritos censales son áreas con un número similar de personas
Objetos espaciales:
- Polígonos de distritos censales
- Polígonos de barrios más grandes
- Polígonos de zonas de agua

Procedimiento

Lee archivos shape que describen barrios y vías de agua
Alinea dos sistemas de coordenadas distintos
Fusiona datos de un data frame en un SpatialPolygonsDataFrame
Pule un mapa hasta dejarlo listo para publicar

Lectura de un archivo shape

Datos vectoriales: descritos por puntos, líneas, polígonos
El archivo shape es el formato más común

Lectura de un archivo shape

# rgdal::readOGR() lee formatos vectoriales
library(rgdal) 
library(sp)

dir()

"water"

dir("water")

"water-areas.dbf" "water-areas.prj" 
"water-areas.shp" "water-areas.shx"

water <- readOGR("water", "water-areas")

Fuente de datos OGR con driver: ESRI Shapefile 
Origen: "water", capa: "water-areas"
con 20 entidades
Tiene 5 campos

Comprobación del resultado

summary(water)

Objeto de clase SpatialPolygonsDataFrame
Coordenadas:
        min       max
x -74.04731 -73.90866
y  40.68419  40.88207
¿Proyectado?: FALSE 
...

plot(water)

Lectura de un archivo ráster

library(rgdal)  # rgdal::readGDAL() lee formatos ráster a objetos sp
library(raster) # raster::raster() lee formatos ráster a objetos raster

dir()

"usgrid_data_2000"   "usgrid_data_2000_1"

dir("usgrid_data_2000")

"metadata"     "usarea00.tif"  "usba00.tif" "usfb00.tif" "usgrid-2000-variables.xls" "usp2500.tif" 
"uspop300.tif" "uspov00.tif"  "uspvp00.tif"

total_pop <- raster("usgrid_data_2000/uspop300.tif")

Comprobación del resultado

total_pop

class       : RasterLayer 
dimensions  : 3120, 7080, 22089600  (nrow, ncol, ncell)
resolution  : 0.008333333, 0.008333333  (x, y)
extent      : -125, -66, 24, 50  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 
data source : /Users/wickhamc/Documents/Projects/courses-visualizing-geospatial-
              data-in-r/data/census_grids/usgrid_data_2000/uspop300.tif 
names       : uspop300 
values      : 0, 65535  (min, max)

¡Vamos a practicar!

Visualización de datos geoespaciales en R