CommencezCommencez gratuitement

Un paquet qui utilise des objets Raster

Vous avez vu que le paquet tmap facilite la visualisation des classes spatiales de sp. Bonne nouvelle : il fonctionne aussi avec les classes de raster ! Il 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 le dernier exercice, vous pouvez afficher l'une de ses couches à l'aide de l'argument col (abréviation de « color ») dans tm_raster(), en mettant le nom de la couche entre guillemets.

Vous utiliserez tmap tout au long du cours, mais nous souhaitons aussi vous présenter un autre paquet, rasterVis, conçu spécialement pour la visualisation des objets raster. Le paquet rasterVis propose plusieurs fonctions pour générer des graphiques, mais essayons-en une seule pour l'instant : levelplot().

Cette activité fait partie du cours

Visualiser des données géospatiales avec R

Voir le cours

Instructions de l’exercice

  • Utilisez tmap pour tracer l'objet pop, en indiquant pop comme argument shp de tm_shape() et en ajoutant une couche tm_raster().
  • Utilisez tmap pour tracer 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 ce code d’exemple.

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