Un pacchetto che usa oggetti Raster
Hai visto che il pacchetto tmap rende facile visualizzare le classi spaziali di sp. La buona notizia è che funziona anche con le classi di raster! Ti basta passare il tuo oggetto Raster___ come argomento shp alla funzione tm_shape(), e poi aggiungere un layer tm_raster() così:
tm_shape(raster_object) +
tm_raster()
Quando lavori con un oggetto RasterStack o RasterBrick, come pop_by_age che hai creato nell'esercizio precedente, puoi mostrare uno dei suoi layer usando l'argomento col (abbreviazione di "color") in tm_raster(), racchiudendo il nome del layer tra virgolette.
Userai tmap per tutto il corso, ma vogliamo anche mostrarti un altro pacchetto, rasterVis, progettato specificamente per visualizzare oggetti raster. In rasterVis ci sono diverse funzioni che puoi usare per creare grafici, ma per ora proviamone una: levelplot().
Questo esercizio fa parte del corso
Visualizzare dati geospaziali in R
Istruzioni dell'esercizio
- Usa
tmapper tracciare l'oggettopop, specificandopopcome argomentoshpditm_shape()e aggiungendo un layertm_raster(). - Usa
tmapper tracciare il layerunder_1dipop_by_age, un oggettoRasterStack. - Chiama la funzione
levelplot()dirasterVissupop.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
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