Carte de l’incendie de la forêt de Houara (107 ha brûlés) – Analyse avec Sentinel-2 et QGIS

Contexte

Dans ce tutoriel, et suite à la demande de collègues, nous mettons à votre disposition la carte de l'incendie récent déclenché le 27/05/2025 dans la forêt de Houara, ainsi qu’un petit tutoriel expliquant sa réalisation. L’objectif est de montrer, étape par étape, comment cartographier un incendie — en l’occurrence celui survenu dans la forêt de Houara, près de Tanger, dans la commune de Gueznaia.

Nous allons utiliser des images satellite gratuites issues de Sentinel-2 pour identifier et quantifier la zone brûlée. Ce genre d’analyse peut s’appliquer à d’autres incendies forestiers au Maroc ou ailleurs.

Télécharger la carte

Données utilisées

• Images satellite : Sentinel-2 (ou Landsat si vous préférez)
• Source : Copernicus Open Access Hub
• Lien utile pour vérifier les passages satellites : Spectator Earth - Acquisition Plan

Exemple de sélection d'images dans Sentinel Browser

🗓️ Dates des images utilisées

Pour cet exemple, les dates sélectionnées sont :

• Avant l'incendie : 25/05/2025
• Après l'incendie : 30/05/2025

Il faut toujours attendre le prochain passage du satellite après l’extinction de l’incendie, pour obtenir une image post-feu exploitable.

Téléchargement des images

Vous pouvez suivre ce tutoriel détaillé sur le téléchargement d’images Sentinel-2 :

Comment télécharger les images Sentinel-2 – sig-maroc.com

Bandes nécessaires pour l’indice de brûlure (NBR)

Pour calculer l’indice NBR (Normalized Burn Ratio), il vous faut :

• B8 : Infrarouge proche (NIR) – résolution 10 m
• B12 : Infrarouge à ondes courtes SWIR2 – résolution 20 m

Facultatif : Vous pouvez aussi télécharger les bandes B2, B3, B4 pour créer une vue en couleurs naturelles (true color) avant/après l’incendie et visualiser les changements.

Calcul des indices dans QGIS

1. Importer vos images dans QGIS (avant et après).

Affichage en couleurs naturelles (B4, B3, B2)

Interface de la calculatrice raster dans QGIS

2. Utiliser la Calculatrice Raster pour calculer les deux indices NBR :
   • NBR = (B8 - B12) / (B8 + B12)
   (à faire une fois pour l’image avant, une fois pour l’image après)

Résultat du calcul du NBR

3. Ensuite, calculer la différence entre les deux :
   • dNBR = NBR_avant – NBR_après

Carte de différence NBR (dNBR)

Classification de la brûlure

Voici les classes standards recommandées par le United States Geological Survey (USGS) pour interpréter la dNBR :

Source : Key and Benson (2006), USGS FireMon / MTBS Program

Extraire la zone brûlée

Si vous voulez uniquement extraire la zone brûlée, vous pouvez appliquer la formule suivante dans la calculatrice raster :

"DNBR" > 0.25

Une fois que c’est fait, et pour convertir le raster en vecteur (ce qui permet de calculer la superficie par la suite), allez dans Raster → Conversion → Polygoniser. Cela vous donnera un fichier shapefile avec deux classes : 0 et 1. Éliminez ceux de la classe 0 et gardez ceux de la classe 1. N’oubliez pas de nettoyer s’il y a du bruit.

La couche vecteur sera généralement un ensemble d’entités à fusionner en une seule pour calculer la superficie facilement : allez dans Vecteur → Outils de géométrie → Collecter les géométries. Sélectionnez la couche concernée et laissez le reste par défaut.

Dans la table d’attributs de la couche résultante, ajoutez un champ superficie, puis calculez la superficie.

Résultat final avec la zone brûlée vectorisée