Un paquete que usa objetos Raster
Ya viste que el paquete tmap facilita la visualización de clases espaciales en sp. La buena noticia es que también funciona con las clases de raster. Simplemente pasa tu objeto Raster___ como argumento shp a la función tm_shape() y añade una capa tm_raster() así:
tm_shape(raster_object) +
tm_raster()
Cuando trabajes con un objeto RasterStack o RasterBrick, como el objeto pop_by_age que creaste en el ejercicio anterior, puedes mostrar una de sus capas usando el argumento col (abreviatura de "color") en tm_raster(), poniendo el nombre de la capa entre comillas.
Usarás tmap a lo largo del curso, pero también queremos enseñarte otro paquete, rasterVis, diseñado específicamente para visualizar objetos de raster. Hay varias funciones en rasterVis que puedes usar para crear gráficos, pero por ahora probemos solo una de ellas: levelplot().
Este ejercicio forma parte del curso
Visualización de datos geoespaciales en R
Instrucciones del ejercicio
- Usa
tmappara representar el objetopop, especificandopopcomo argumentoshpentm_shape()y añadiendo una capatm_raster(). - Usa
tmappara representar la capaunder_1depop_by_age, un objetoRasterStack. - Llama a la función
levelplot()derasterVissobrepop.
ejercicio interactivo práctico
Prueba este ejercicio completando este código de ejemplo.
library(tmap)
# Specify pop as the shp and add a tm_raster() layer
tm_shape(___) +
___
# Plot the under_1 layer in pop_by_age
library(rasterVis)
# Call levelplot() on pop