Un objet plus complexe
Le package raster fournit l’objet RasterLayer, mais aussi deux objets plus complexes : RasterStack et RasterBrick. Ces deux objets sont conçus pour stocker de nombreux rasters, tous avec les mêmes emprises et dimensions (aussi appelés rasters multibandes ou multicouches).
Vous pouvez voir RasterLayer comme une matrice, tandis que les objets RasterStack et RasterBrick ressemblent davantage à des tableaux à trois dimensions. Une chose supplémentaire à savoir pour les manipuler est la façon de spécifier une couche particulière.
Vous pouvez utiliser la sélection avec $ ou [[ sur un RasterStack ou un RasterBrick pour récupérer une couche et obtenir en retour un nouvel objet RasterLayer. Par exemple, si x est un RasterStack, x$layer_name ou x[["layer_name"]] renverra un RasterLayer ne contenant que la couche appelée layer_name.
Examinons un objet RasterStack appelé pop_by_age qui couvre la même zone que pop, mais contient désormais des couches pour la population, réparties en plusieurs groupes d’âge.
Cet exercice fait partie du cours
Visualiser des données géospatiales avec R
Instructions
- Affichez
pop_by_age. Voyez-vous les noms de toutes les couches ? - Extrayez la couche
under_1en utilisant la sélection[[. - Tracez la couche
under_1en passant votre code de l’instruction précédente àplot().
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# Print pop_by_age
# Subset out the under_1 layer using [[
# Plot the under_1 layer