Balíček pro práci s rastrovými objekty
Viděl/a jsi, že balíček tmap usnadňuje vizualizaci prostorových tříd z sp. Dobrá zpráva je, že funguje i s třídami z balíčku raster! Stačí předat objekt Raster___ jako argument shp funkci tm_shape() a přidat vrstvu tm_raster() takto:
tm_shape(raster_object) +
tm_raster()
Při práci s objektem RasterStack nebo RasterBrick, jako je objekt pop_by_age z předchozího cvičení, můžeš zobrazit jednu z jeho vrstev pomocí argumentu col (zkratka z anglického „color") ve funkci tm_raster() — název vrstvy uzavři do uvozovek.
S balíčkem tmap budeš pracovat po celý kurz, ale chceme ti ukázat také balíček rasterVis, který je navržený přímo pro vizualizaci rastrových objektů. V rasterVis najdeš několik různých funkcí pro tvorbu grafů — teď si vyzkoušíme jednu z nich: levelplot().
Toto cvičení je součástí kurzu
Vizualizace geoprostorových dat v R
Pokyny k cvičení
- Pomocí
tmapzobraz objektpop: zadejpopjako argumentshpfunkcetm_shape()a přidej vrstvutm_raster(). - Pomocí
tmapzobraz vrstvuunder_1objektupop_by_age, což je objekt typuRasterStack. - Zavolej funkci
levelplot()z balíčkurasterVisna objektpop.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
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