Um pacote que usa objetos Raster
Você viu que o pacote tmap facilita a visualização de classes espaciais em sp. A boa notícia é que ele também funciona com as classes do raster! Basta passar seu objeto Raster___ como o argumento shp para a função tm_shape() e, em seguida, adicionar uma camada tm_raster() assim:
tm_shape(raster_object) +
tm_raster()
Ao trabalhar com um objeto RasterStack ou RasterBrick, como o objeto pop_by_age que você criou no último exercício, você pode exibir uma de suas camadas usando o argumento col (abreviação de "color") em tm_raster(), colocando o nome da camada entre aspas.
Você vai usar o tmap ao longo do curso, mas também queremos apresentar outro pacote, o rasterVis, projetado especificamente para visualizar objetos raster. Há algumas funções diferentes que você pode usar no rasterVis para criar gráficos, mas por enquanto vamos experimentar apenas uma delas: levelplot().
Este exercício faz parte do curso
Visualizando dados geoespaciais em R
Instruções do exercício
- Use
tmappara plotar o objetopop, especificandopopcomo o argumentoshpdetm_shape()e adicionando uma camadatm_raster(). - Use
tmappara plotar a camadaunder_1depop_by_age, um objetoRasterStack. - Chame a função
levelplot()dorasterVisempop.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
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