Pakiet obsługujący obiekty Raster
Pakiet tmap świetnie radzi sobie z wizualizacją klas przestrzennych z pakietu sp. Dobra wiadomość jest taka, że działa równie dobrze z klasami raster! Wystarczy przekazać obiekt Raster___ jako argument shp do funkcji tm_shape(), a następnie dodać warstwę tm_raster() w taki sposób:
tm_shape(raster_object) +
tm_raster()
Pracując z obiektem RasterStack lub RasterBrick, takim jak pop_by_age utworzony w poprzednim ćwiczeniu, możesz wyświetlić jedną z jego warstw, używając argumentu col (skrót od „color") w tm_raster() – nazwę warstwy umieść w cudzysłowie.
Z pakietem tmap będziesz pracować przez cały kurs, ale chcemy ci też pokazać inny pakiet – rasterVis – zaprojektowany specjalnie do wizualizacji obiektów raster. Dostępnych jest kilka różnych funkcji, ale na razie wypróbujemy jedną z nich: levelplot().
To ćwiczenie jest częścią kursu
Wizualizacja danych geoprzestrzennych w R
Instrukcje do ćwiczenia
- Za pomocą
tmapwyświetl obiektpop, podającpopjako argumentshpfunkcjitm_shape()i dodając warstwętm_raster(). - Za pomocą
tmapwyświetl warstwęunder_1obiektupop_by_age– obiektu klasyRasterStack. - Wywołaj funkcję
levelplot()z pakieturasterVisna obiekciepop.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
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