Use Case 3 - Hinweiskarten zur Vitalität/de - Versionsgeschichte

FuzzyBot: Übernehme Bearbeitung einer neuen Version der Quellseite

2021-09-06T07:17:41Z

Übernehme Bearbeitung einer neuen Version der Quellseite

FuzzyBot: Übernehme Bearbeitung einer neuen Version der Quellseite

2021-09-06T07:14:10Z

Übernehme Bearbeitung einer neuen Version der Quellseite

FuzzyBot: Übernehme Bearbeitung einer neuen Version der Quellseite

2021-09-06T07:03:21Z

Übernehme Bearbeitung einer neuen Version der Quellseite

FuzzyBot: Übernehme Bearbeitung einer neuen Version der Quellseite

2021-07-08T12:32:49Z

Übernehme Bearbeitung einer neuen Version der Quellseite

FuzzyBot: Übernehme Bearbeitung einer neuen Version der Quellseite

2021-07-01T09:37:40Z

Übernehme Bearbeitung einer neuen Version der Quellseite

FuzzyBot: Übernehme Bearbeitung einer neuen Version der Quellseite

2021-07-01T09:15:14Z

Übernehme Bearbeitung einer neuen Version der Quellseite

FuzzyBot: Übernehme Bearbeitung einer neuen Version der Quellseite

2021-06-15T14:09:01Z

Übernehme Bearbeitung einer neuen Version der Quellseite

FuzzyBot: Übernehme Bearbeitung einer neuen Version der Quellseite

2021-04-06T14:04:33Z

Übernehme Bearbeitung einer neuen Version der Quellseite

FuzzyBot: Übernehme Bearbeitung einer neuen Version der Quellseite

2021-03-25T09:19:08Z

Übernehme Bearbeitung einer neuen Version der Quellseite

Neue Seite

[[Datei:Uc3_screenshot.png|thumb|right|500px|Abb. 1: Screenshot aus der [https://waldmonitoring.ch/vitalitaet?x=829300&y=5933555&zoom=13.0000 Webanwendung]. Rot deutet auf eine Abnahme der Vitalität hin, blau auf eine Zunahme. Die Farbintensität spiegelt die Stärke der Veränderung wieder.]]

[[Use-Cases Waldmonitoring | ''-> Zurück zur Übersicht über die Use Cases'']]

[https://waldmonitoring.ch/vitalitaet?x=829300&y=5933555&zoom=13.0000 -> Zur Webanwendung]

Auch die Vitalität von Waldflächen kann mit Sentinel-2-Satellitenbildern grob erfasst werden. Durch den '''Vergleich des aktuellen Vegetationszustandes mit den Vorjahren''' lassen sich daraus '''Hinweiskarten zur Veränderung der Vitalität''' erstellen. Hier handelt es sich allerdings um einen ersten Test. Eine Validierung und Verfeinerung der Methode steht noch aus.

Als Indikator für die Vitalität wurde wiederum der NDVI-Vegetationsindex verwendet. Auf der Karte dargestellt werden sogenannte '''NDVI-Anomalien'''. Dabei werden jeweils die NDVI-Medianwerte des zweimonatigen Zeitfensters (z.B. Juni - Juli) mit den Medianwerten aller Vorjahre (bis 2015) innerhalb derselben Zeitperiode verglichen (siehe Abb. 2). Negative Werte deuten auf eine Abnahme der Vitalität hin, positive Werte deuten auf eine Zunahme der Vitalität hin. Je weiter die Werte von null abweichen (Erwartungswert), desto wahrscheinlicher ist es, dass eine effektive Veränderung stattfand.

==Einschränkungen==
Ob es sich bei einer Veränderung um Borkenkäferbefall, Trockenstress oder einen Holzschlag handelt, wird nicht unterschieden. So können z.B. negative Werte sowohl einen Holzschlag wie auch eine vorzeitige Herbstverfärbung beschreiben. Wir sprechen daher von Hinweiskarten und für die Interpretation ist immer auch Expertenwissen über die Wälder und gegebenenfalls eine Feldbegehung notwendig.

Durch die Kombination von Satellitenaufnahmen innerhalb des zweimonatigen Fensters können viele Probleme mit Wolken und andere Fehlerquellen reduziert, jedoch nicht vollständig ausgeschlossen werden. Derzeit treten noch fehlerhafte Indikationen im Zusammenhang mit der Abdeckung durch unterschiedliche Orbits, sowie durch nicht ausmaskierte Wolken oder eine ungenügende Zahl wolkenfreier Aufnahmen auf. Eine verbesserte Wolkenmaske sowie eine integrierte Überprüfung der Anzahl wolkenfreier Pixel pro Komposit befinden sich in Entwicklung (Stand Januar 2021).

Weiter ist zu berücksichtigen, dass die Sentinel-2-Daten erst seit 2015 verfügbar sind, was für die Detektion von Anomalien bis jetzt ein relativ kurzer Zeitraum ist. Mit jedem zusätzlichen Jahr wird eine höhere Aussagekraft erwartet.

==Video-Tutorial : Basics für den korrekten Gebrauch==
Das folgende Video erläutert die Benutzung der Webanwendung «Hinweiskarten zur Vitalität»:

[[File:UC3.mp4|1000px]]

==Hintergrundinformationen zur Methode==
Die Berechnung der NDVI-Anomalien basiert auf dem sogenannten [https://de.wikipedia.org/wiki/Standardisierung_(Statistik) Z-Wert] (siehe z.B. [https://www.sciencedirect.com/science/article/pii/S0034425718305509 Meroni et al., 2019]):

Zi = (NDVImedian,i - NDVImedian,ref) / NDVIsd,ref

Im Gegensatz zu Meroni et al. haben wir den Median anstatt dem arithmetischen Mittel verwendet, da dieser robuster gegenüber Ausreissern ist. NDVImedian,i bezeichnet somit den Median der NDVI-Werte pro zweimonatigem Zeitfenster i. NDVImedian,ref bezeichnet den Median der NDVI-Werte innerhalb der gleichen zweimonatigen Zeitfenster (z.B. Juni - Juli) der Referenzperiode. NDVIsd,ref ist die Standardabweichung aller NDVI-Werte innerhalb der gleichen zweimonatigen Zeitfenster (z.B. Juni- Juli) der Referenzperiode.

Die Standardabweichung gibt Auskunft über die Variabilität der NDVI-Werte innerhalb der Referenzperiode. Je grösser die Variabilität ist, desto kleiner wird die Wahrscheinlichkeit, dass Veränderungen, insbesondere schwache, ausgeschieden werden, was sich in Z-Werten nahe null widerspiegelt. Die Berechnung erfolgt pixelbasiert für jeden 10 x 10 m Pixel, und die Pixelwerte in den Hinweiskarten entsprechen dem jeweiligen Z-Wert.

[[Datei:NDVI-Beispiel_UC3.png|750px|left|thumb|Abb. 2: NDVI-Zeitreihen-Beispiel zur Veranschaulichung der Z-Wert-Methode. Die blauen Punkte entsprechen den tatsächlichen NDVI-Werten an einem bestimmten Ort (pixelbasiert). Die blauen Linien stellen Interpolationen dar und haben keine weitere Bedeutung. Die Referenzperiode bezieht sich in diesem Fall auf 2016 bis 2018, dargestellt in Grün. 2019 stellt in diesem Beispiel den aktuellen Zustand dar und ist in Rot gekennzeichnet. Die gestrichelten Linien stellen jeweils die beiden Mediane dar, die miteinander verglichen werden, das heisst der Median über Juni/Juli aller Jahre der Referenzperiode (in Grün) mit dem Median über Juni/Juli des aktuellen Jahres (in Rot).]]

@@ Zeile 41: / Zeile 41: @@
 [[Datei:NDVI-Beispiel_UC3.png|750px|left|thumb|Abb. 2: NDVI-Zeitreihen-Beispiel zur Veranschaulichung der Z-Wert-Methode. Die blauen Punkte entsprechen den tatsächlichen NDVI-Werten an einem bestimmten Ort (pixelbasiert). Die blauen Linien stellen Interpolationen dar und haben keine weitere Bedeutung. Die Referenzperiode bezieht sich in diesem Fall auf 2016 bis 2018, dargestellt in Grün. 2019 stellt in diesem Beispiel den aktuellen Zustand dar und ist in Rot gekennzeichnet. Die gestrichelten Linien stellen jeweils die beiden Mediane dar, die miteinander verglichen werden, das heisst der Median über Juni/Juli aller Jahre der Referenzperiode (in Grün) mit dem Median über Juni/Juli des aktuellen Jahres (in Rot).]]
-<br clear=all>
+<br clear=all/>
 == Einzelnachweise==

@@ Zeile 22: / Zeile 22: @@
 [[File:UC3.mp4|1000px]]
 ==Hintergrundinformationen zur Methode==
@@ Zeile 36: / Zeile 41: @@
 [[Datei:NDVI-Beispiel_UC3.png|750px|left|thumb|Abb. 2: NDVI-Zeitreihen-Beispiel zur Veranschaulichung der Z-Wert-Methode. Die blauen Punkte entsprechen den tatsächlichen NDVI-Werten an einem bestimmten Ort (pixelbasiert). Die blauen Linien stellen Interpolationen dar und haben keine weitere Bedeutung. Die Referenzperiode bezieht sich in diesem Fall auf 2016 bis 2018, dargestellt in Grün. 2019 stellt in diesem Beispiel den aktuellen Zustand dar und ist in Rot gekennzeichnet. Die gestrichelten Linien stellen jeweils die beiden Mediane dar, die miteinander verglichen werden, das heisst der Median über Juni/Juli aller Jahre der Referenzperiode (in Grün) mit dem Median über Juni/Juli des aktuellen Jahres (in Rot).]]
-==Datenbezug==
+<br clear=all>
 == Einzelnachweise==

@@ Zeile 1: / Zeile 1: @@
 [[Datei:Uc3_screenshot.png|thumb|right|500px|Abb. 1: Screenshot aus der [https://waldmonitoring.ch/vitalitaet?x=829300&y=5933555&zoom=13.0000 Webanwendung]. Rot deutet auf eine Abnahme der Vitalität hin, blau auf eine Zunahme. Die Farbintensität spiegelt die Stärke der Veränderung wieder.]]
@@ Zeile 13: / Zeile 14: @@
 Ob es sich bei einer Veränderung um Borkenkäferbefall, Trockenstress oder einen Holzschlag handelt, wird nicht unterschieden. So können z.B. negative Werte sowohl einen Holzschlag wie auch eine vorzeitige Herbstverfärbung beschreiben. Wir sprechen daher von Hinweiskarten und für die Interpretation ist immer auch Expertenwissen über die Wälder und gegebenenfalls eine Feldbegehung notwendig.
-Durch die Kombination von Satellitenaufnahmen innerhalb des zweimonatigen Fensters können viele Probleme mit Wolken und andere Fehlerquellen reduziert, jedoch nicht vollständig ausgeschlossen werden. Es können fehlerhafte Indikationen im Zusammenhang mit der Abdeckung durch unterschiedliche Orbits, sowie durch nicht ausmaskierte Wolken oder eine ungenügende Zahl wolkenfreier Aufnahmen auftreten. Deshalb wurden in der finalen Hinweiskarte diejenigen Pixel ausmaskiert, für die innerhalb der Monitoring-Periode für die Berechnung des Medians weniger als 6 wolkenfreie Aufnahmen (Pixel) vorhanden sind.
+Durch die Kombination von Satellitenaufnahmen innerhalb des zweimonatigen Fensters können viele Probleme mit Wolken und andere Fehlerquellen reduziert, jedoch nicht vollständig ausgeschlossen werden. Im Zusammenhang mit der Abdeckung durch unterschiedliche Orbits treten teilweise Artefakte auf, und die Grenzen der Orbits sind auf den Hinweiskarten erkennbar. Auch können fehlerhafte Indikationen durch nicht ausmaskierte Wolken oder eine ungenügende Zahl wolkenfreier Aufnahmen auftreten. Deshalb wurden in der finalen Hinweiskarte diejenigen Pixel ausmaskiert, für die innerhalb der Monitoring-Periode für die Berechnung des Medians weniger als 6 wolkenfreie Aufnahmen (Pixel) vorhanden sind.
 Weiter ist zu berücksichtigen, dass die Sentinel-2-Daten erst seit 2015 verfügbar sind, was für die Detektion von Anomalien bis jetzt ein relativ kurzer Zeitraum ist. Mit jedem zusätzlichen Jahr wird eine höhere Aussagekraft erwartet.

@@ Zeile 6: / Zeile 6: @@
 [https://waldmonitoring.ch/vitalitaet?x=829300&y=5933555&zoom=13.0000 -> Zur Webanwendung]
-Auch die Vitalität von Waldflächen kann mit Sentinel-2-Satellitenbildern grob erfasst werden. Durch den '''Vergleich des aktuellen Vegetationszustandes mit den Vorjahren''' lassen sich daraus '''Hinweiskarten zur Veränderung der Vitalität''' erstellen. Hier handelt es sich allerdings um einen ersten Test. Eine Validierung und Verfeinerung der Methode steht noch aus.
+Auch die Vitalität von Waldflächen kann mit Sentinel-2-Satellitenbildern grob erfasst werden. Durch den '''Vergleich des aktuellen Vegetationszustandes mit den Vorjahren''' lassen sich daraus '''Hinweiskarten zur Veränderung der Vitalität''' erstellen. Hier handelt es sich allerdings um einen ersten Test. Eine ausführliche Validierung der Methode steht noch aus.
 Als Indikator für die Vitalität wurde wiederum der NDVI-Vegetationsindex verwendet. Auf der Karte dargestellt werden sogenannte '''NDVI-Anomalien'''. Dabei werden jeweils die NDVI-Medianwerte des zweimonatigen Zeitfensters (z.B. Juni - Juli) mit den Medianwerten aller Vorjahre (bis 2015) innerhalb derselben Zeitperiode verglichen (siehe Abb. 2). Negative Werte deuten auf eine Abnahme der Vitalität hin, positive Werte deuten auf eine Zunahme der Vitalität hin. Je weiter die Werte von null abweichen (Erwartungswert), desto wahrscheinlicher ist es, dass eine effektive Veränderung stattfand. Dabei kann ein '''absoluter Wert von mindestens 3,5''' als '''Indikator für eine signifikante Abweichung''' betrachtet werden <ref>Boris Iglewicz and David Hoaglin (1993), "Volume 16: How to Detect and Handle Outliers", The ASQC Basic References in Quality Control: Statistical Techniques, Edward F. Mykytka, Ph.D., Editor.</ref>.

@@ Zeile 8: / Zeile 8: @@
 Auch die Vitalität von Waldflächen kann mit Sentinel-2-Satellitenbildern grob erfasst werden. Durch den '''Vergleich des aktuellen Vegetationszustandes mit den Vorjahren''' lassen sich daraus '''Hinweiskarten zur Veränderung der Vitalität''' erstellen. Hier handelt es sich allerdings um einen ersten Test. Eine Validierung und Verfeinerung der Methode steht noch aus.
-Als Indikator für die Vitalität wurde wiederum der NDVI-Vegetationsindex verwendet. Auf der Karte dargestellt werden sogenannte '''NDVI-Anomalien'''. Dabei werden jeweils die NDVI-Medianwerte des zweimonatigen Zeitfensters (z.B. Juni - Juli) mit den Medianwerten aller Vorjahre (bis 2015) innerhalb derselben Zeitperiode verglichen (siehe Abb. 2). Negative Werte deuten auf eine Abnahme der Vitalität hin, positive Werte deuten auf eine Zunahme der Vitalität hin. Je weiter die Werte von null abweichen (Erwartungswert), desto wahrscheinlicher ist es, dass eine effektive Veränderung stattfand. Dabei kann ein absoluter Z-Wert von mindestens 3,5 als Indikator für eine signifikante Abweichung betrachtet werden <ref>Boris Iglewicz and David Hoaglin (1993), "Volume 16: How to Detect and Handle Outliers", The ASQC Basic References in Quality Control: Statistical Techniques, Edward F. Mykytka, Ph.D., Editor.</ref>.
+Als Indikator für die Vitalität wurde wiederum der NDVI-Vegetationsindex verwendet. Auf der Karte dargestellt werden sogenannte '''NDVI-Anomalien'''. Dabei werden jeweils die NDVI-Medianwerte des zweimonatigen Zeitfensters (z.B. Juni - Juli) mit den Medianwerten aller Vorjahre (bis 2015) innerhalb derselben Zeitperiode verglichen (siehe Abb. 2). Negative Werte deuten auf eine Abnahme der Vitalität hin, positive Werte deuten auf eine Zunahme der Vitalität hin. Je weiter die Werte von null abweichen (Erwartungswert), desto wahrscheinlicher ist es, dass eine effektive Veränderung stattfand. Dabei kann ein '''absoluter Wert von mindestens 3,5''' als '''Indikator für eine signifikante Abweichung''' betrachtet werden <ref>Boris Iglewicz and David Hoaglin (1993), "Volume 16: How to Detect and Handle Outliers", The ASQC Basic References in Quality Control: Statistical Techniques, Edward F. Mykytka, Ph.D., Editor.</ref>.
 ==Einschränkungen==

@@ Zeile 13: / Zeile 13: @@
 Ob es sich bei einer Veränderung um Borkenkäferbefall, Trockenstress oder einen Holzschlag handelt, wird nicht unterschieden. So können z.B. negative Werte sowohl einen Holzschlag wie auch eine vorzeitige Herbstverfärbung beschreiben. Wir sprechen daher von Hinweiskarten und für die Interpretation ist immer auch Expertenwissen über die Wälder und gegebenenfalls eine Feldbegehung notwendig.
-Durch die Kombination von Satellitenaufnahmen innerhalb des zweimonatigen Fensters können viele Probleme mit Wolken und andere Fehlerquellen reduziert, jedoch nicht vollständig ausgeschlossen werden. Derzeit treten noch fehlerhafte Indikationen im Zusammenhang mit der Abdeckung durch unterschiedliche Orbits, sowie durch nicht ausmaskierte Wolken oder eine ungenügende Zahl wolkenfreier Aufnahmen auf. Eine verbesserte Wolkenmaske sowie eine integrierte Überprüfung der Anzahl wolkenfreier Pixel pro Komposit befinden sich in Entwicklung (Stand Januar 2021).
+Durch die Kombination von Satellitenaufnahmen innerhalb des zweimonatigen Fensters können viele Probleme mit Wolken und andere Fehlerquellen reduziert, jedoch nicht vollständig ausgeschlossen werden. Es können fehlerhafte Indikationen im Zusammenhang mit der Abdeckung durch unterschiedliche Orbits, sowie durch nicht ausmaskierte Wolken oder eine ungenügende Zahl wolkenfreier Aufnahmen auftreten. Deshalb wurden in der finalen Hinweiskarte diejenigen Pixel ausmaskiert, für die innerhalb der Monitoring-Periode für die Berechnung des Medians weniger als 5 wolkenfreie Aufnahmen (Pixel) vorhanden sind.
 Weiter ist zu berücksichtigen, dass die Sentinel-2-Daten erst seit 2015 verfügbar sind, was für die Detektion von Anomalien bis jetzt ein relativ kurzer Zeitraum ist. Mit jedem zusätzlichen Jahr wird eine höhere Aussagekraft erwartet.
@@ Zeile 25: / Zeile 25: @@
 Die Berechnung der NDVI-Anomalien basiert auf dem sogenannten [https://de.wikipedia.org/wiki/Standardisierung_(Statistik) Z-Wert] (siehe z.B. [https://www.sciencedirect.com/science/article/pii/S0034425718305509 Meroni et al., 2019]):
-Z<sub>i</sub> = (NDVI<sub>median,i</sub> - NDVI<sub>median,ref</sub>) / NDVI<sub>sd,ref</sub>
+Z<sub>i</sub> = (NDVI<sub>median,i</sub> - NDVI<sub>median,ref</sub>) / NDVI<sub>MAD,ref</sub>
-Im Gegensatz zu Meroni et al. haben wir den Median anstatt dem arithmetischen Mittel verwendet, da dieser robuster gegenüber Ausreissern ist. NDVI<sub>median,i</sub> bezeichnet somit den Median der NDVI-Werte pro zweimonatigem Zeitfenster i. NDVI<sub>median,ref</sub> bezeichnet den Median der NDVI-Werte innerhalb der gleichen zweimonatigen Zeitfenster (z.B. Juni - Juli) der Referenzperiode. NDVI<sub>sd,ref</sub> ist die Standardabweichung aller NDVI-Werte innerhalb der gleichen zweimonatigen Zeitfenster (z.B. Juni- Juli) der Referenzperiode.
+Im Gegensatz zu Meroni et al. haben wir den sogenannten [https://www.ibm.com/docs/en/cognos-analytics/11.1.0?topic=terms-modified-z-score modifizierten Z-Score] verwendet, da dieser robuster gegenüber Ausreissern ist. Dabei wird das arithmetische Mittel durch den Median und die Standardabweichung durch die  [https://de.wikipedia.org/wiki/Mittlere_absolute_Abweichung_vom_Median mittlere absolute Abweichung vom Median (MAD)] ersetzt. NDVI<sub>median,i</sub> bezeichnet somit den Median der NDVI-Werte pro zweimonatigem Monitoring-Zeitfenster i. NDVI<sub>median,ref</sub> bezeichnet den Median der NDVI-Werte innerhalb der gleichen zweimonatigen Zeitfenster (z.B. Juni - Juli) der Referenzperiode. NDVI<sub>MAD,ref</sub> ist die mittlere absolute Abweichung der NDVI-Werte vom Median innerhalb der gleichen zweimonatigen Zeitfenster (z.B. Juni- Juli) der Referenzperiode.
-Die Standardabweichung gibt Auskunft über die Variabilität der NDVI-Werte innerhalb der Referenzperiode. Je grösser die Variabilität ist, desto kleiner wird die Wahrscheinlichkeit, dass Veränderungen, insbesondere schwache, ausgeschieden werden, was sich in Z-Werten nahe null widerspiegelt. Die Berechnung erfolgt pixelbasiert für jeden 10 x 10 m Pixel, und die Pixelwerte in den Hinweiskarten entsprechen dem jeweiligen Z-Wert.
+Die mittlere absolute Abweichung vom Median gibt Auskunft über die Variabilität der NDVI-Werte innerhalb der Referenzperiode. Je grösser die Variabilität ist, desto kleiner wird die Wahrscheinlichkeit, dass Veränderungen, insbesondere schwache, ausgeschieden werden, was sich in Z-Werten nahe null widerspiegelt. Die Berechnung erfolgt pixelbasiert für jeden 10 x 10 m Pixel, und die Pixelwerte in den Hinweiskarten entsprechen dem jeweiligen Z-Wert.
 [[Datei:NDVI-Beispiel_UC3.png|750px|left|thumb|Abb. 2: NDVI-Zeitreihen-Beispiel zur Veranschaulichung der Z-Wert-Methode. Die blauen Punkte entsprechen den tatsächlichen NDVI-Werten an einem bestimmten Ort (pixelbasiert). Die blauen Linien stellen Interpolationen dar und haben keine weitere Bedeutung. Die Referenzperiode bezieht sich in diesem Fall auf 2016 bis 2018, dargestellt in Grün. 2019 stellt in diesem Beispiel den aktuellen Zustand dar und ist in Rot gekennzeichnet. Die gestrichelten Linien stellen jeweils die beiden Mediane dar, die miteinander verglichen werden, das heisst der Median über Juni/Juli aller Jahre der Referenzperiode (in Grün) mit dem Median über Juni/Juli des aktuellen Jahres (in Rot).]]

← Nächstältere Version		Version vom 6. April 2021, 14:04 Uhr
Zeile 2:		Zeile 2:


−	[[Use-Cases Waldmonitoring \| ''-> Zurück zur Übersicht über die Use Cases'']]	+	[[Special:MyLanguage/Use-Cases Waldmonitoring \| ''-> Zurück zur Übersicht über die Use Cases'']]

	[https://waldmonitoring.ch/vitalitaet?x=829300&y=5933555&zoom=13.0000 -> Zur Webanwendung]		[https://waldmonitoring.ch/vitalitaet?x=829300&y=5933555&zoom=13.0000 -> Zur Webanwendung]