Rapport hebdomadaire 2

Author

Alix Biton

Published

January 16, 2026

Semaine 2 (19-23/01)

Objectifs

  • Récupérer les données d’élévation sur Copernicus DEM
  • Comprendre la fonction hillShade() du package raster
  • Compute l’exposition topographique au vent sur la zone de Mayotte

Copernicus DEM

Copernicus Browser

Pour récupérer les données d’altitudes d’une zone géographique, il faut aller dans le Copernicus Browser, onglet Visualise. Il faut ensuite cadrer avec le polygone/carré la zone d’intérêt, en l’occurence ici, la zone de Mayotte. Laisser le mode Configuration sur Default et choisir dans Data Collections Copernicus DEM et choisir Copernicus 30. Pour télécharger le fichier .tiff, appuyer sur l’icône Download Image sur la barre d’icônes à droite de la page. Aller dans l’onglet Analytical et ici :

  • Image format : TIFF (32-bit float), équivalent du format DGED
  • Image resolution : CUSTOM -> 30x30, équivalent du GLO30
  • Coordinate System : WGS 84 (EPSG:4326)
  • Layers : Raw -> DEM

Cliquer sur Download et un fichier .tiff est téléchargé.

Visualisation de l’élévation de Mayotte (avec le package terra)

Compute Exposition Topographique

Sur le cas du cyclone Chido, j’ai récupéré l’exposition topographique au vent pour chaque mesure prise dans mon loi (location of interest) qui est la zone de Mayotte (~300km). Cette exposition topographique est décrite comme une variable quantitative “Exposure” et peut être lue comme un indice d’exposition au vent pour chaque cellule d’une grille (un raster). Les valeurs négatives représente les aires protégées du vent tandis que les valeurs positives représente les aires exposées au vent.

Pour calculer cette métrique, il faut récupérer les données de Profiles présentes dans la fonction spatialBehaviour() ainsi que la pente (slope) en radians et l’exposition/orientation des versants (aspect) en radians.

pf <- spatialBehaviour(sts, product = "Profiles")
slope <- terrain(dem, v ="slope", unit = "radians")
exposition <- terrain(dem, v = "aspect", unit = "radians")

Dans Ibanez et al.2024, la fonction utilisée pour calculer l’exposition topographique au vent est raster::hillShade(). Cependant, ce package n’est pas dans les dépendances du package StormR contrairement à terra. De plus le package terra est un package, plus jeune, de remplacement de raster. On utilisera donc la fonction terra::shade() pour calculer l’exposition topographique au vent.

Pour utiliser la fonction terra::shade(), il faut un angle et une direction. L’angle utilisé sera l’angle d’inflexion de 6° (Boose et al.1994). Pour la direction, on veut récupérer la direction que prend le vent lorsque sa vitesse est maximale dans un layer c’est-à-dire la direction de la cellule ayant le Vmax.

On obtient donc un raster appelé ici topo comprenant l’exposition topographique au vent à chaque mesure à partir d’une direction. On peut ainsi récupérer le layer ayant la plus haute exposition en moyenne.

Les flèches blanches ici représentent l’étendue des directions de vents tout au long du cyclone tandis que la flèche rouge représente la direction pour le layer considéré.

Requête HTTPS avec OData

Pour effectuer une requête des données Copernicus_DEM_GLO30_DGED sans utiliser Copernicus Browser, on utilise OData. Dans le terminal il faut d’abord s’identifier et obtenir un token avec la commande :

ACCESS_TOKEN=$(curl -s \
  -d "grant_type=password" \
  -d "client_id=cdse-public" \  -d "scope=openid profile email" \
  -d "username=biton.alix@gmail.com" \
  -d "password=ryWvox-fishum-qoxhu8" \
  https://identity.dataspace.copernicus.eu/auth/realms/CDSE/protocol/openid-connect/token \
  | jq -r .access_token)