CommencerCommencer gratuitement

Un package qui utilise des objets Raster

Vous avez vu que le package tmap facilite la visualisation des classes spatiales de sp. Bonne nouvelle : il fonctionne aussi avec les classes de raster ! Il vous suffit de passer votre objet Raster___ comme argument shp à la fonction tm_shape(), puis d’ajouter une couche tm_raster() comme ceci :

tm_shape(raster_object) +
    tm_raster()

Lorsque vous travaillez avec un objet RasterStack ou RasterBrick, comme l’objet pop_by_age que vous avez créé dans l’exercice précédent, vous pouvez afficher l’une de ses couches en utilisant l’argument col (abréviation de « color ») de tm_raster(), en mettant le nom de la couche entre guillemets.

Vous utiliserez tmap tout au long du cours, mais nous souhaitons aussi vous montrer un autre package, rasterVis, conçu spécifiquement pour visualiser des objets raster. Il existe plusieurs fonctions dans rasterVis pour créer des graphiques, mais testons-en une pour l’instant : levelplot().

Cet exercice fait partie du cours

Visualiser des données géospatiales avec R

Afficher le cours

Instructions

  • Utilisez tmap pour afficher l’objet pop, en indiquant pop comme argument shp de tm_shape() et en ajoutant une couche tm_raster().
  • Utilisez tmap pour afficher la couche under_1 de pop_by_age, un objet RasterStack.
  • Appelez la fonction levelplot() de rasterVis sur pop.

Exercice interactif pratique

Essayez cet exercice en complétant cet exemple de code.

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
Modifier et exécuter le code