EmpezarEmpieza gratis

Convertir de un CRS/proyección a otro

El proceso de convertir de un CRS o proyección a otro lo gestionan los métodos spTransform() del paquete rgdal. spTransform() tiene métodos para todos los objetos de sp, incluido SpatialPolygonsDataFrame, pero no funciona con objetos de raster. Esto se debe a que transformar un ráster es un poco más complicado; la cuadrícula rectangular transformada dejará de ser rectangular. Si te da curiosidad cómo transformar rásteres, puedes consultar ?raster::projectRaster.

La transformación es sencilla. El primer argumento de spTransform(), x, es el objeto espacial que quieres transformar y el segundo, CRS, es la especificación del CRS deseado. El CRS puede indicarse con una cadena PROJ4, que podrías construir a mano, pero es mucho más fácil tomarla de un objeto existente (por ejemplo, con la función proj4string()).

Es hora de poner tus dos conjuntos de polígonos en el mismo CRS.

Este ejercicio forma parte del curso

Visualización de datos geoespaciales en R

Ver curso

Instrucciones del ejercicio

  • Transforma neighborhoods para que tenga el mismo CRS que nyc_tracts usando spTransform() con el argumento CRS establecido a proj4string(nyc_tracts).
  • Verifica la transformación consultando el head() de coordinates(neighborhoods).
  • Comprueba que ahora los conjuntos de datos encajan dibujando neighborhoods, luego nyc_tracts con add = TRUE y col = "red", y por último water con add = TRUE y col = "blue".

ejercicio interactivo práctico

Prueba este ejercicio completando este código de ejemplo.

library(sp)
library(raster)

# Use spTransform on neighborhoods: neighborhoods



# head() on coordinates() of neighborhoods


# Plot neighborhoods, nyc_tracts and water
Editar y ejecutar código