Um objeto mais complicado
O pacote raster fornece o objeto RasterLayer, mas também alguns objetos mais complexos: RasterStack e RasterBrick. Esses dois objetos são projetados para armazenar muitos rasters, todos com a mesma extensão e dimensão (também chamados de rasters de múltiplas bandas ou múltiplas camadas).
Você pode pensar em RasterLayer como uma matriz, mas os objetos RasterStack e RasterBrick são mais parecidos com arrays tridimensionais. Uma coisa adicional que você precisa saber para trabalhar com eles é como especificar uma camada específica.
Você pode usar a indexação com $ ou [[ em um RasterStack ou RasterBrick para pegar uma camada e retornar um novo objeto RasterLayer. Por exemplo, se x é um RasterStack, x$layer_name ou x[["layer_name"]] vai retornar um RasterLayer contendo apenas a camada chamada layer_name.
Vamos analisar um objeto RasterStack chamado pop_by_age que cobre a mesma área que pop, mas agora contém camadas de população divididas em alguns grupos etários.
Este exercício faz parte do curso
Visualizando dados geoespaciais em R
Instruções do exercício
- Imprima
pop_by_age. Você consegue ver os nomes de todas as camadas? - Selecione a camada
under_1usando a indexação[[. - Plote a camada
under_1passando seu código da instrução anterior paraplot().
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# Print pop_by_age
# Subset out the under_1 layer using [[
# Plot the under_1 layer