EmpezarEmpieza gratis

Combinar datos con distintos CRS/proyecciones

Cada objeto espacial tiene asociado un sistema de referencia de coordenadas (CRS). Por lo general, se define al importar los datos y se lee directamente de los archivos espaciales. Así es como neighborhoods y nyc_tracts obtuvieron la información de su sistema de coordenadas.

Tanto los paquetes sp como raster tienen una función proj4string() que devuelve el CRS del objeto sobre el que se llama.

Intentar trabajar con datos espaciales que usan CRS diferentes es un poco como trabajar con un conjunto de datos en millas y otro en kilómetros. Miden lo mismo, pero los números no son directamente comparables.

Echemos un vistazo a nuestros dos objetos de polígonos.

Este ejercicio forma parte del curso

Visualización de datos geoespaciales en R

Ver curso

Instrucciones del ejercicio

  • Llama a proj4string() sobre neighborhoods y luego sobre nyc_tracts. Verifica que las dos cadenas son diferentes.
  • Mira el head() de las coordinates() de neighborhoods y repite para nyc_tracts. ¿Ves el problema? nyc_tracts tiene coordenadas x alrededor de -70, pero neighborhoods está alrededor de 1.000.000.
  • Representa neighborhoods y luego representa nyc_tracts con col = "red" y add = TRUE para añadirlos encima.

ejercicio interactivo práctico

Prueba este ejercicio completando este código de ejemplo.

library(sp)

# proj4string() on nyc_tracts and neighborhoods



# coordinates() on nyc_tracts and neighborhoods



# plot() neighborhoods and nyc_tracts
Editar y ejecutar código