Use Case 2 - Test Sommersturmschäden 2017/fr: Unterschied zwischen den Versionen

Aus Waldmonitoring Wiki
(Die Seite wurde neu angelegt: „[https://waldmonitoring.ch/stoerungen?x=829300&y=5933555&zoom=13.0000 -> vers l'application web]“)
(Übernehme Bearbeitung einer neuen Version der Quellseite)
 
(19 dazwischenliegende Versionen von einem anderen Benutzer werden nicht angezeigt)
Zeile 3: Zeile 3:
 
[https://waldmonitoring.ch/stoerungen?x=829300&y=5933555&zoom=13.0000 -> vers l'application web]
 
[https://waldmonitoring.ch/stoerungen?x=829300&y=5933555&zoom=13.0000 -> vers l'application web]
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
[[Datei:UC2_Screenshot.png|thumb|500px|Fig. 1 : Épreuve écran issue de l’application web. La légende permet de différencier les surfaces de dégât potentielles (jaune à rouge) des nuages (hachuré gris) et des surfaces sans données (grises). Le symbole entouré en bleu contient un lien vers Sentinel Playground permettant de vérifier visuellement l’image satellite de base.]]Une méthode pour fournir de manière automatique des '''cartes indicatrices sur les dégâts de tempête estivales pour toute la Suisse''' a été développée d’après les images satellite Sentinel-2. L’année 2017 a été choisie comme exemple et pour un test pratique de l’évaluation des tempêtes estivales (fin [http://www.sturmarchiv.ch/index.php?title=Extremereignisse_2017#Juli juillet] et début [https://www.sturmarchiv.ch/index.php?title=Extremereignisse_2017#August août]). L’utilisateur ou l’utilisatrice peut choisir toutes les images satellite Sentinel-2 entre le 05.07.2017 et le 19.08.2017 qui couvrent au moins en partie la Suisse. Il est alors possible, après une tempête, de vérifier rapidement si des images utilisables, c’est-à-dire avec le moins de nuages possible, sont disponibles. Les surfaces de changement potentielles disponibles s’affichent ensuite dans le visualiseur.
[[Datei:UC2_Screenshot.png|thumb|500px|Abb. 1: Screenshot aus der Webanwendung. Mithilfe der Legende können potenzielle Schadflächen (gelb bis rot), Wolken (grau schraffiert) und Flächen mit keinen Daten (grau ausgefüllt) unterschieden werden. Das blau umrandete Symbol stellt einen Link auf Sentinel Playground bereit, sodass das zugrundeliegende Satellitenbild einfach visuell überprüft werden kann.]]Im Rahmen des zweiten Use-Cases wurde eine Methode für die automatische Bereitstellung von '''schweizweiten Hinweiskarten für Sommersturmschäden''' auf Basis von Sentinel-2-Satellitenbildern entwickelt. Als Beispiel und für den Praxistest wurden die Sommerstürme (Ende [http://www.sturmarchiv.ch/index.php?title=Extremereignisse_2017#Juli Juli] und Anfang [https://www.sturmarchiv.ch/index.php?title=Extremereignisse_2017#August August]) im Jahr 2017 ausgewählt. Zwischen dem 5.7.2017 und dem 19.8.2017 können alle verfügbaren Sentinel-2-Aufnahmen, welche mindestens einen Teil der Schweiz abdecken, vom Benutzer/Benutzerin ausgewählt werden. So kann nach einem Sturmereignis rasch geprüft werden, ob eine brauchbare, sprich möglichst wolkenfreie, Aufnahme zur Verfügung steht. Daraufhin werden im Kartenviewer potenzielle Veränderungsflächen angezeigt, falls welche vorhanden sind.
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
Les surfaces de changement proposées sont déterminées d’après la diminution de l’indice de végétation Normalized Burn Ratio (NBR). Un projet précédent ([https://www.planfor.ch/de/content/schlussbericht-waldmonitoring-mit-sentinel-2-satellitenbildern Weber & Rosset, 2019]) a démontré l’utilité de cet indice pour détecter les surfaces touchées par des tempêtes.
Diese ausgewiesenen Veränderungsflächen basieren auf der Abnahme des Normalized Burn Ratio Vegetationsindexes (NBR). In einem Vorgängerprojekt ([https://www.planfor.ch/de/content/schlussbericht-waldmonitoring-mit-sentinel-2-satellitenbildern Weber & Rosset, 2019]) hat sich dieser Index als geeignet für die Detektion von Sturmschadflächen erwiesen.
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
[[Datei:orbits.png|thumb|500px|Fig. 2 : Représentation des deux orbites Sentinel couvrant la Suisse. À une date donnée, l’orbite 108 ou l’orbite 65 est relevée. La disponibilité des images est la meilleure pour la zone où elles se superposent (tous les 2 ou 3 jours).]] Les valeurs dans la légende du visualiseur représentent en moyenne pour chaque surface de changement (polygone) la '''baisse du NBR''' multipliée par 100. La multiplication par 100 ne sert ici qu’à économiser de l’espace de stockage (grâce aux chiffres entiers plutôt que décimaux). La différence résulte de la comparaison entre '''l’image à la date choisie''' et un état de référence avant l’événement présumé / la date choisie. On utilise comme référence un '''composite si possible exempt de nuage''' constitué de toutes les images disponibles datant de '''45 jours avant la date choisie'''. Des valeurs proches de -100 indiquent d’importants dégâts. Les polygones de changement sont délimités grâce à une valeur seuil (-15) et à partir d’une superficie minimum de 5 ares. Les nuages sont aussi représentés (cf. Fig. 1) afin d’éviter de conclure erronément à l’absence de dégâts. Cela vaut aussi pour les surfaces « NoData », c.-à-d. celles pour lesquelles aucune donnée Sentinel-2 n’est disponible pour la date choisie (Fig. 1).  
[[Datei:orbits.png|thumb|500px|Abb. 2: Darstellung der beiden Sentinel-Orbits, welche die Fläche der Schweiz abdecken. An einem Aufnahmedatum wird entweder Orbit 108 oder 65 erfasst. Für den überlappenden Bereich ist die Bildverfügbarkeit am höchsten (alle 2-3 Tage).]]Die Werte in der Legende des Kartenviewers stellen die '''Abnahme des NBR''' multipliziert mit 100 dar, gemittelt pro Veränderungsfläche (Polygon). Die Multiplikation mit 100 erfolgt dabei lediglich aus Speicherplatzgründen (Ganzzahl statt Dezimalzahl). Die Differenzbildung erfolgt aus dem '''Bild des ausgewählten Datums''' und einem Referenzzustand vor dem vermuteten Ereignis / ausgewählten Datum. Als Referenz wird ein möglichst '''wolkenfreies Komposit''' aller verfügbaren Bilder der '''45 Tage vor dem ausgewählten Datum''' verwendet. Werte näher bei -100 weisen auf stärkere Schäden hin. Die Veränderungspolygone wurden dabei mittels Schwellenwert (-15) und ab einer Mindestgrösse von 5 Aren ausgeschieden. Auch Wolken werden angezeigt (siehe Abb. 1), um fehlerhafte Rückschlüsse auf das Nichtvorhandensein von Schäden zu vermeiden. Das gleiche gilt für sogenannte NoData- Flächen, das heisst Flächen, für die zum gewählten Zeitpunkt keine Sentinel-2-Daten zur Verfügung stehen (siehe Abb. 1).
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
La Suisse étant couverte par deux orbites, une date de prise de vue n’en concerne toujours qu’une partie (Fig. 2). L’autre partie est indiquée comme une surface « NoData » dans le visualiseur.
Da die Schweiz durch zwei Orbits abgedeckt wird, wird an einem Aufnahmedatum immer nur ein Teil der Fläche der Schweiz erfasst (siehe Abb. 2). Der andere Teil wird im Kartenviewer jeweils als NoData-Fläche gekennzeichnet.
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
== Limitations ==
== Einschränkungen ==
+
[[Datei:Beispiel_UC2.png|thumb|500px|Fig. 3 : Délimitation correcte de dégâts de surface survenus le 02.08.2017 dans des noyers, TG. La première image sans nuage disponible est celle du 15.08.2017. Les surfaces de dégâts trouvées ont été comparées avec les données de référence du WSL (en jaune).]]
[[Datei:Beispiel_UC2.png|thumb|500px|Abb. 3: Gute räumliche Abgrenzung von flächigen Schäden für ein Sturmereignis in Nussbaumen TG am 2.8.2017. Die erste wolkenfreie Aufnahme nach dem Sturm war am 15.8.2017 verfügbar. Die ausgewiesenen Schadflächen wurden mit Referenzdaten der WSL (in Gelb) abgeglichen.]]
+
Les résultats ont été en partie validés grâce aux données de référence du WSL et les dégâts de surface ont pu être bien reconnus et délimités (Fig. 3). Certaines surfaces ont pourtant jusqu’à présent été erronément délimitées comme des surfaces de dégât. Il s’agit notamment de zones en bordure de nuages, couvertes de neige ou encore de zones d’ombre sur les versants nord abrupts. Dans ces cas, il faut être particulièrement prudent lors de l’interprétation des résultats.
Die Resultate wurden teilweise mit Referenzdaten der WSL validiert und erwiesene flächige Schäden konnten gut erkannt und abgegrenzt werden (siehe Abb. 3). Jedoch werden bis dato auch diverse Flächen fälschlicherweise als Schadflächen ausgeschieden. Dies geschieht insbesondere an Wolkenrändern, sowie im Zusammenhang mit Schnee oder Schattenwurf an steilen Nordhängen. In diesen Fällen ist bei der Interpretation der Ergebnisse besondere Vorsicht geboten.
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
Il faut si possible toujours choisir une '''date de prise de vue sans nuage'''. C’est pourquoi chaque date de prise de vue dispose d’un lien direct vers Sentinel Playground (cf. symbole entouré de bleu Fig. 1) permettant de rapidement vérifier chaque image visuellement. La date de la prise de vue et l’extrait de l’image sont sauvegardés dans le lien et directement repris. Pour les noyers de TG, il s’est par exemple écoulé 13 jours entre la tempête du 02.08.2017 et la première prise de vue sans nuage (Fig. 3).  
Nach Möglichkeit sollte immer ein '''wolkenfreier Aufnahmezeitpunkt''' ausgewählt werden. Zu diesem Zweck wird für jedes Aufnahmedatum ein Direkt-Link auf Sentinel Playground bereitgestellt (siehe blau umrandetes Symbol in Abb. 1), sodass jedes Bild einfach visuell überprüft werden kann. Aufnahmedatum und Bildausschnitt werden dabei im Link gespeichert und direkt übernommen. In Nussbaumen TG betrug die Wartezeit zwischen dem Sturm vom 2.8.2017 und der nächsten verfügbaren wolkenfreien Aufnahme zum Beispiel 13 Tage (siehe Abb. 3).
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
Le procédé de sélection des images utilisables est expliqué et illustré dans la vidéo ci-dessous.
Das genaue Vorgehen zur Auswahl einer brauchbaren Aufnahme wird im untenstehenden Video noch erklärt und veranschaulicht.
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
L’hiver complique encore la situation : l’état de la végétation, l’intensité lumineuse, les nuages, les ombres et la neige représentent autant de défis qui expliquent pourquoi les résultats des tempêtes hivernales n’ont pas encore pu être publiés. Une solution consisterait à utiliser en combinaison des données Sentinel-1, qui ne sont pas influencées par l’intensité de la couverture nuageuse ou lumineuse (cf. à ce sujet le [https://www.wsl.ch/de/projekte/sturmhinweiskarte.html projet actuel du WSL]).
Eine anspruchsvolle Situation stellt sich zudem im Winter: Vegetationszustand, Beleuchtungsintensität, Wolken, Schatten und Schnee stellen zusätzliche Herausforderungen dar, weshalb Ergebnisse zu Winterstürmen bisher nicht veröffentlicht wurden. Ein möglicher Lösungsansatz wäre eine Kombination mit Sentinel-1 Daten, welche von Bewölkungs- und Beleuchtungsintensität nicht beeinflusst werden (siehe dazu ein [https://www.wsl.ch/de/projekte/sturmhinweiskarte.html laufendes Projekt an der WSL]).
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
== Tutoriel vidéo : Fondamentaux pour une bonne utilisation ==
== Video-Tutorial : Basics für den korrekten Gebrauch ==
+
Cette vidéo montre comment utiliser l’application web « Test des dégâts des tempêtes estivales de 2017 »:  
Das folgende Video erläutert die Benutzung der Webanwendung «Test Sommersturmschäden 2017»:
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
 
 
[[File:UC2.mp4|1000px]]
 
[[File:UC2.mp4|1000px]]
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
== Données ==
== Hintergrundinformationen zur Methode ==
+
Les surfaces de changement peuvent être récupérées [https://forestmonitoring.lab.karten-werk.ch/services ici] sous la forme d’un unique layer (couche) de polygones via un service WMS ou WFS (vous trouverez [[Tutorials#Einbindung.2FDownload_der_Daten_aus_den_Use_Cases| ici]] des conseils et des informations pour intégrer et exploiter des données). En plus du moment (time), la superficie (area), la moyenne (mean), la valeur maximale (max) et le quantile 90% (90p) des valeurs différentielles du NBR sont calculés et sauvegardés comme attributs dans les données vecteur de chaque surface de changement. Comme évoqué précédemment, les différences de NBR sont multipliées par 100 pour des raisons d’espace de stockage et les surfaces de changement sont délimitées grâce à une valeur seuil (-15) et à partir d’une superficie minimum de 5 ares. L’Attribute class différencie les surfaces de changement (class = 1) des nuages (class = -1) et des NoData (class = -2).  
Für jede Sentinel-2-Aufnahme innerhalb der definierten Zeitspanne (5.7. bis 19.8.2017) wurde zuerst der NBR berechnet (NBR<sub>Aktuell</sub>). Wolken wurden dabei mittels der [https://sentinel.esa.int/web/sentinel/technical-guides/sentinel-2-msi/level-1c/cloud-masks ESA Wolkenmaske] ausgeschieden. Daraufhin wurde der NBR für die '''Referenzperiode''' berechnet, welche als das '''45-Tages-Fenster''' vor dem ausgewählten Aufnahmedatum definiert wurde.
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
== Informations contextuelles sur la méthode ==
Zu diesem Zweck wurde zuerst mittels der unter Use Case 1 beschriebenen NDVI Maximum Komposit Methode aus allen verfügbaren Aufnahmen innerhalb der 45-tägigen Zeitspanne ein möglichst wolkenfreies Komposit erstellt. Für dieses wurde dann wiederum der NBR berechnet (NBR<sub>Referenz</sub>). Aus diesen beiden Bildern wurde daraufhin die Differenz gebildet (∆NBR), und die Veränderungsflächen wurden mittels Schwellenwert ausgeschieden. Die nachfolgende Abbildung veranschaulicht das Vorgehen.
+
Au cours de la période définie (05.07 au 19.08.2017), le NBR a d’abord été calculé (NBR<sub>Aktuell</sub>). Les nuages ont été délimités grâce au [https://sentinel.esa.int/web/sentinel/technical-guides/sentinel-2-msi/level-1c/cloud-masks masque de nuage ESA]. Le NBR a ensuite été calculé pour la '''période de référence''', celle-ci étant définie comme une '''fenêtre de 45 jours''' avant la date du relevé.
</div>
 
  
 +
Pour cela, un composite avec le moins de nuages possible a été établi avec toutes les images de la période de 45 jours selon la méthode du NDVI composite maximum décrite dans le cas d’usage 1. Le NBR (NBR<sub>Referenz</sub>) a ensuite été à nouveau calculé pour celui-ci. À partir de ces deux images, on peut établir la différence (∆NBR) et les surfaces de changement sont délimitées grâce à la valeur seuil. La figure suivante illustre le procédé.
  
<div lang="de" dir="ltr" class="mw-content-ltr">
 
[[Datei:Workflow-UC2.png|thumb|750px|left|Abb. 4: Darstellung des Arbeitsflusses zu Use-Case 2. Für jede Sentinel-2-Aufnahme wird der NBR berechnet (NBR Aktuell). Daraufhin wird der NBR für die Referenzperiode berechnet, welche als das 45-Tages-Fenster vor dem ausgewählten Aufnahmedatum definiert wurde. Zu diesem Zweck wird aus allen verfügbaren Aufnahmen innerhalb der 45-tägigen Zeitspanne ein möglichst wolkenfreies Komposit erstellt, für welches wiederum der NBR berechnet wird (NBR Referenz). Aus diesen beiden Bildern wird die Differenz gebildet (∆NBR), und die Veränderungsflächen werden mittels Schwellenwert ausgeschieden.]]
 
</div>
 
  
<div lang="de" dir="ltr" class="mw-content-ltr">
+
[[Datei:Workflow-UC2.png|thumb|750px|left|Fig. 4 : Représentation du flux de travail pour le cas d’usage 2. Le NBR est calculé pour pour chaque prise de vue Sentinel-2 (NBR Aktuell). Le NBR est ensuite calculé pour la période de référence, celle-ci étant définie comme une fenêtre de 45 jours avant la date du relevé. Pour cela, un composite avec le moins de nuages possible a été établi avec toutes les images de la période de 45 jours. Le NBR est calculé pour ce composite (NBR Referenz). À partir de ces deux images, on peut établir la différence (∆NBR) et délimiter les surfaces de changement grâce à la valeur seuil.]]
 +
 
 
<br clear=all>
 
<br clear=all>
Dieser Prozess erfolgte wie beschrieben zunächst testweise für den Sommer 2017, jedoch ist die automatisierte Bereitstellung der Veränderungsflächen für die ganze Schweiz möglich. Dabei würden die Informationen für alle verfügbaren Aufnahmen innerhalb der letzten 45 Tage ab dem jeweils aktuellen Datum mittels sogenannten «rollenden Archiven» zur Verfügung gestellt werden. Die vollautomatische Analyse und Bereitstellung der Resultate sollte innerhalb von 2-5 Tagen nach Bildaufnahme möglich sein.
+
Bien que cette procédure ait d’abord été employée de manière expérimentale pour l’été 2017, comme décrit ci-avant, il est possible d’obtenir automatiquement les surfaces de changement pour toute la Suisse. Il faut pour cela disposer des informations pour toutes les images disponibles des 45 derniers jours avant la date actuelle, grâce à des « archives roulantes ». L’analyse et la mise à disposition automatisée des résultats se fait normalement en 2 à 5 jours après la prise vue.
</div>
 
 
 
<div lang="de" dir="ltr" class="mw-content-ltr">
 
== Datenbezug ==
 
Die Veränderungsflächen können [https://forestmonitoring.lab.karten-werk.ch/services hier] als ein einziger Polygon-Layer via WMS oder WFS Dienst bezogen werden (Tipps und Informationen zu Einbindung und Verarbeitung der Daten gibt es [[Tutorials#Einbindung.2FDownload_der_Daten_aus_den_Use_Cases|hier]]). Für jede ausgeschiedene Veränderungsfläche wurden neben dem Zeitpunkt (time) zudem Flächengrösse (area), sowie Mittelwert (mean), Maximalwert (max) und 90%-Quantil (p90) der NBR-Differenzwerte berechnet und als Attribute in den Vektordaten gespeichert. Wie schon erwähnt, wurden die NBR Differenzen aus Speicherplatzgründen mit 100 mulipliziert, und Veränderungslfächen mittels Schwellenwert (-15) und ab einer Mindestgrösse von 5 Aren ausgeschieden. Das Attribut class differenziert zwischen Veränderungsflächen (class = 1), Wolken (class = -1), und NoData (class = -2).
 
</div>
 

Aktuelle Version vom 6. September 2021, 07:24 Uhr

-> vers la vue d'ensemble des cas d'usage

-> vers l'application web

Fig. 1 : Épreuve écran issue de l’application web. La légende permet de différencier les surfaces de dégât potentielles (jaune à rouge) des nuages (hachuré gris) et des surfaces sans données (grises). Le symbole entouré en bleu contient un lien vers Sentinel Playground permettant de vérifier visuellement l’image satellite de base.

Une méthode pour fournir de manière automatique des cartes indicatrices sur les dégâts de tempête estivales pour toute la Suisse a été développée d’après les images satellite Sentinel-2. L’année 2017 a été choisie comme exemple et pour un test pratique de l’évaluation des tempêtes estivales (fin juillet et début août). L’utilisateur ou l’utilisatrice peut choisir toutes les images satellite Sentinel-2 entre le 05.07.2017 et le 19.08.2017 qui couvrent au moins en partie la Suisse. Il est alors possible, après une tempête, de vérifier rapidement si des images utilisables, c’est-à-dire avec le moins de nuages possible, sont disponibles. Les surfaces de changement potentielles disponibles s’affichent ensuite dans le visualiseur.

Les surfaces de changement proposées sont déterminées d’après la diminution de l’indice de végétation Normalized Burn Ratio (NBR). Un projet précédent (Weber & Rosset, 2019) a démontré l’utilité de cet indice pour détecter les surfaces touchées par des tempêtes.

Fig. 2 : Représentation des deux orbites Sentinel couvrant la Suisse. À une date donnée, l’orbite 108 ou l’orbite 65 est relevée. La disponibilité des images est la meilleure pour la zone où elles se superposent (tous les 2 ou 3 jours).

Les valeurs dans la légende du visualiseur représentent en moyenne pour chaque surface de changement (polygone) la baisse du NBR multipliée par 100. La multiplication par 100 ne sert ici qu’à économiser de l’espace de stockage (grâce aux chiffres entiers plutôt que décimaux). La différence résulte de la comparaison entre l’image à la date choisie et un état de référence avant l’événement présumé / la date choisie. On utilise comme référence un composite si possible exempt de nuage constitué de toutes les images disponibles datant de 45 jours avant la date choisie. Des valeurs proches de -100 indiquent d’importants dégâts. Les polygones de changement sont délimités grâce à une valeur seuil (-15) et à partir d’une superficie minimum de 5 ares. Les nuages sont aussi représentés (cf. Fig. 1) afin d’éviter de conclure erronément à l’absence de dégâts. Cela vaut aussi pour les surfaces « NoData », c.-à-d. celles pour lesquelles aucune donnée Sentinel-2 n’est disponible pour la date choisie (Fig. 1).

La Suisse étant couverte par deux orbites, une date de prise de vue n’en concerne toujours qu’une partie (Fig. 2). L’autre partie est indiquée comme une surface « NoData » dans le visualiseur.

Limitations

Fig. 3 : Délimitation correcte de dégâts de surface survenus le 02.08.2017 dans des noyers, TG. La première image sans nuage disponible est celle du 15.08.2017. Les surfaces de dégâts trouvées ont été comparées avec les données de référence du WSL (en jaune).

Les résultats ont été en partie validés grâce aux données de référence du WSL et les dégâts de surface ont pu être bien reconnus et délimités (Fig. 3). Certaines surfaces ont pourtant jusqu’à présent été erronément délimitées comme des surfaces de dégât. Il s’agit notamment de zones en bordure de nuages, couvertes de neige ou encore de zones d’ombre sur les versants nord abrupts. Dans ces cas, il faut être particulièrement prudent lors de l’interprétation des résultats.

Il faut si possible toujours choisir une date de prise de vue sans nuage. C’est pourquoi chaque date de prise de vue dispose d’un lien direct vers Sentinel Playground (cf. symbole entouré de bleu Fig. 1) permettant de rapidement vérifier chaque image visuellement. La date de la prise de vue et l’extrait de l’image sont sauvegardés dans le lien et directement repris. Pour les noyers de TG, il s’est par exemple écoulé 13 jours entre la tempête du 02.08.2017 et la première prise de vue sans nuage (Fig. 3).

Le procédé de sélection des images utilisables est expliqué et illustré dans la vidéo ci-dessous.

L’hiver complique encore la situation : l’état de la végétation, l’intensité lumineuse, les nuages, les ombres et la neige représentent autant de défis qui expliquent pourquoi les résultats des tempêtes hivernales n’ont pas encore pu être publiés. Une solution consisterait à utiliser en combinaison des données Sentinel-1, qui ne sont pas influencées par l’intensité de la couverture nuageuse ou lumineuse (cf. à ce sujet le projet actuel du WSL).

Tutoriel vidéo : Fondamentaux pour une bonne utilisation

Cette vidéo montre comment utiliser l’application web « Test des dégâts des tempêtes estivales de 2017 »:

Données

Les surfaces de changement peuvent être récupérées ici sous la forme d’un unique layer (couche) de polygones via un service WMS ou WFS (vous trouverez ici des conseils et des informations pour intégrer et exploiter des données). En plus du moment (time), la superficie (area), la moyenne (mean), la valeur maximale (max) et le quantile 90% (90p) des valeurs différentielles du NBR sont calculés et sauvegardés comme attributs dans les données vecteur de chaque surface de changement. Comme évoqué précédemment, les différences de NBR sont multipliées par 100 pour des raisons d’espace de stockage et les surfaces de changement sont délimitées grâce à une valeur seuil (-15) et à partir d’une superficie minimum de 5 ares. L’Attribute class différencie les surfaces de changement (class = 1) des nuages (class = -1) et des NoData (class = -2).

Informations contextuelles sur la méthode

Au cours de la période définie (05.07 au 19.08.2017), le NBR a d’abord été calculé (NBRAktuell). Les nuages ont été délimités grâce au masque de nuage ESA. Le NBR a ensuite été calculé pour la période de référence, celle-ci étant définie comme une fenêtre de 45 jours avant la date du relevé.

Pour cela, un composite avec le moins de nuages possible a été établi avec toutes les images de la période de 45 jours selon la méthode du NDVI composite maximum décrite dans le cas d’usage 1. Le NBR (NBRReferenz) a ensuite été à nouveau calculé pour celui-ci. À partir de ces deux images, on peut établir la différence (∆NBR) et les surfaces de changement sont délimitées grâce à la valeur seuil. La figure suivante illustre le procédé.


Fig. 4 : Représentation du flux de travail pour le cas d’usage 2. Le NBR est calculé pour pour chaque prise de vue Sentinel-2 (NBR Aktuell). Le NBR est ensuite calculé pour la période de référence, celle-ci étant définie comme une fenêtre de 45 jours avant la date du relevé. Pour cela, un composite avec le moins de nuages possible a été établi avec toutes les images de la période de 45 jours. Le NBR est calculé pour ce composite (NBR Referenz). À partir de ces deux images, on peut établir la différence (∆NBR) et délimiter les surfaces de changement grâce à la valeur seuil.


Bien que cette procédure ait d’abord été employée de manière expérimentale pour l’été 2017, comme décrit ci-avant, il est possible d’obtenir automatiquement les surfaces de changement pour toute la Suisse. Il faut pour cela disposer des informations pour toutes les images disponibles des 45 derniers jours avant la date actuelle, grâce à des « archives roulantes ». L’analyse et la mise à disposition automatisée des résultats se fait normalement en 2 à 5 jours après la prise vue.