Utiliser la boîte à outils des traitements pour automatiser la classification de la neige¶
Snowy Hydro Limited (SHL) est un générateur de pic d’électricité basé dans les Montagnes Neigeuses de la Nouvelle-Galles du Sud, en Australie. Bien que possédant plusieurs ouvrages de génération d’électricité à travers le pays, Snowy Hydro Limited est mieux connu pour son infrastructure emblématique en Australie, le schéma hydroelectrique des Montagnes Neigeuses.
Le schéma hydroélectrique des Montagnes Neigeuses, est rendu possible par des chutes de neige saisonnières dans la région alpine de Nouvelle-Galles du Sud. L’eau issue de la fonte de la neige de printemps est utilisée pour produire de l’électricité pour gérer les pics de demande sur le marché électrique national, tout en détournant également l’eau vers les zones plus sèches de l’intérieur de l’Australie du Sud-Est, pour l’irrigation et l’agriculture.
La classification de la couverture de neige au cours des mois d’hiver permet à SHL de contrôler le cycle de vie de la zone recouverte par la neige, d’archiver ces zones de couverture pour conserver des enregistrements historiques, et de fournir les statistiques de couverture à l’équipe des services hydrographiques pour la production d’eau et les prévisions de stockage.
Map 1.1 Carte de bassin versant de Snowy¶
Contexte¶
Chaque année pendant les mois d’hiver (du 1 juin au 31 octobre), SHL acquiert autant d’images sans nuages que possible à partir du satellite MODIS Terra (EOS AM). Ces images sont alors intégrées dans un algorithme connu sous le nom l’Index Normalisé de Différence d’Enneigement (NDSI), qui est utilisé pour indiquer les zones où la neige est présente. Alors que cet algorithme n’est pas réservé à un capteur particulier, l’instrument MODIS a été choisi en raison de sa résolution spatiale adéquate et de sa couverture quotidienne de la zone du projet (Carte 1.1 du bassin versant Snowy).
Objectif¶
Le but de ce projet est d’automatiser le processus NDSI sous la forme d’un modèle dans la boîte à outils des traitements de QGIS. La sortie attendue est un polygone vecteur montrant la zone enneigée à l’intérieur du bassin versant de Snowy ainsi que des attributs décrivant la surface totale recouverte, en kilomètres carrés.
Méthode¶
L’algorithme NDSI reste assez simple (comme montré dans la formule 1.1. Index Normalisé de Différence d’Enneigement (NDS)) et les étapes successives pour extraire la couverture neigeuse sous forme de vecteur se prêtent à être construites dans un modèle de traitement.
Formule 1.1 - Index normalisé de différence d’enneigement (NDSI)¶
L’algorithme NDSI renverra des résultats compris entre -1 et 1 avec un seuil d’enneigement qui occupera les valeurs supérieures à 0.4. Un exemple de sortie de l’algorithme est présenté ci-dessous (Carte 1.3 Sortie NDSI) ainsi que l’image en couleur naturelle (Carte 1.2 - Image MODIS en couleurs vraies) pour la date du 13 juillet 2014. Comme vous pouvez l’observer sur les images, la sortie NDSI permet de particulièrement bien discriminer les zones enneigées des autres types de couverture ou des nuages. Néanmoins, il peut toutefois mal classer des portions de surface en eau de grande taille comme zones enneigées car l’eau a des caractéristiques d’absorption dans l’infra-rouge. Cette caractéristique de NDSI indique un besoin d’incorporer un masque sur les zones en eau dans une étape du modèle.
Carte 1.2 - Image MODIS en couleurs vraies¶
Carte 1.3 Sortie NDSI¶
Carte 1.4 - Reclassification « Neige » vs « Sans Neige »¶
La sortie de NDSI est créé à l’aide d’une couche à bande unique en utilisant le module GRASS « r.mapcalculator ». Une fois la grille NDSI disponible, le seuil de 0.4 doit être appliqué et il permettra de générer une grille booléenne permettant de déterminer les zones « Neige » des zones « Sans Neige ». Pour obtenir cette classification, un autre module de GRASS « r.decode » a été utilisé et a permis de produire le résultat présenté dans la Carte 1.4 - Reclassification « Neige » vs « Sans Neige ».
La grille reclassifiée a été ensuite transformée en jeu de données polygone vecteur et les zones enneigées ont été extraites (comme présenté ci-dessous dans la Carte 1.5 - Extraction de la neige). Ces zones ont été découpées selon le bassin versant Snowy (Carte 1.6 - Neige à l’intérieur du bassin versant Snowy). Ensuite, on a appliqué le masque sur les grandes étendues d’eau (Carte 1.7 - Zones d’enneigement finales) et finalement les statistiques de zones d’enneigement sont générées.
Carte 1.5 - Extraction de la neige¶
Carte 1.6 - Neige à l’intérieur du bassin versant Snowy¶
Carte 1.7 - Zones d’enneigement finales¶
La sortie du modèle est un polygone vecteur indiquant les zones enneigées situées à l’intérieur du bassin versant Snowy. Les attributs du polygone incluent la surface couverte par la neige et le pourcentage du bassin qui est recouvert. Le modèle peut être aperçu dans son intégralité dans l’image ci-dessous: « Image 1.1 - Modeleur de traitements de QGIS pour l’Index Normalisé de Différence d’Enneigement (NDSI) ».
Image 1.1 - Modeleur de traitements de QGIS pour l’Index Normalisé de Différence d’Enneigement (NDSI)¶
Résultats/Conclusion¶
La boîte à outils de traitements de QGIS est un dispositif puissant pour automatiser un processus tel que la classification d’enneigement. Le modèle présenté ci-dessus automatise correctement le processus de classification d’enneigement et fournit un jeu de données SIG qui peut être visualisé et être retravaillé dans d’autres outils statistiques. Il permet de supprimer les erreurs subjectives des méthodes antérieurement utilisées par SHL tout en fournissant un outil cohérent et utilisable de manière répétitive. Voici quelques résultats selon les saisons.