Een pakket dat Raster-objecten gebruikt
Je hebt gezien dat het tmap-pakket het eenvoudig maakt om ruimtelijke klassen in sp te visualiseren. Het goede nieuws is dat het ook werkt met de raster-klassen! Je geeft simpelweg je Raster___-object door als het shp-argument aan de functie tm_shape(), en voegt daarna een tm_raster()-laag toe, zoals dit:
tm_shape(raster_object) +
tm_raster()
Als je werkt met een RasterStack- of RasterBrick-object, zoals het pop_by_age-object dat je in de vorige oefening hebt gemaakt, kun je een van de lagen weergeven met het col-argument (kort voor "color") in tm_raster(), waarbij je de laagnaam tussen aanhalingstekens zet.
Je gaat in de hele cursus met tmap werken, maar we willen je ook een ander pakket laten zien, rasterVis, dat speciaal is ontworpen voor het visualiseren van raster-objecten. Er zijn een paar verschillende functies in rasterVis om grafieken te maken, maar laten we er nu één proberen: levelplot().
Deze oefening maakt deel uit van de cursus
Georuimtelijke data visualiseren in R
Oefeninstructies
- Gebruik
tmapom hetpop-object te plotten doorpopalsshp-argument aantm_shape()te geven en eentm_raster()-laag toe te voegen. - Gebruik
tmapom deunder_1-laag vanpop_by_age, eenRasterStack-object, te plotten. - Roep de
rasterVis-functielevelplot()aan oppop.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
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