Stichproben aus Rasterdaten mit Hilfe von Punkten und Polygonen (QGIS3)

Viele wissenschaftliche Daten oder Umweltdaten liegen als Rasterdaten vor. Digitale Geländemodelle (DGM oder DEM) werden ebenso als Rasterdaten angeboten. In diesen Rastendaten ist der thematische Parameter im Pixelwert einer Rasterzelle hinterlegt. Oftmals benötigt man diesen Pixelwert für einen bestimmten Ort oder auch eine Aggreation der Pixelwerte für ein bestimmtes Gebiet. Diese Funktionalität ist in den Verarbeitungswerkzeugen von QGIS enthalten. Die Algorithmen heißen Rasterwerte abtasten für Punktlayer und Zonenstatistik für Polygonlayer.

Überblick über die Aufgabe

Mit einem gegebenen Rasterdatensatz der maximalen Temperatur der Kontinetalfläche der USA wollen wir die Temperatur für Punkte der Ballungsgebiete extrahieren. Außerdem wollen wir für die Fläche der Verwaltungsbezirke die mittlere Temperatur bestimmen.

Weitere Fähigkeiten, die wir erlernen

• Selektion und Entfernen mehrerer Layer aus dem QGIS Inhaltsverzeichnis.

Beschaffung der Daten

Das Climate Prediction Center der NOAA stellt GIS-Daten zu Temperatur und Niederschlag in den USA bereit: GIS data. Die aktuellen Daten können unter folgendem link heruntergeladen werden: latest grid file for maximum temperatures. Die Datei wird den folgenden Namen erhalten: us.tmax_nohads_ll_{YYYYMMDD}_float.tif

Wir werden folgende CSV Datei mit den Ballungsräumen der USA verwenden: 2018 US Gazetteer. Das Ortsverzeichnis laden wir unter:Urban Areas Gazetteer File herunter.

US Census Bureau stellt folgende Daten bereit: TIGER/Line Shapefiles. Auf der folgenden Seite kann man die Flächeneinheiten der Volkszählung Kaliforniens herunterladen: FTP site. Wir wählen die Datei Census Tracts for California.

Der Einfachheit halber können wir eine Kopie der Datensätze unter folgenden links herunterladen:

us.tmax_nohads_ll_20190501_float.tif

2018_Gaz_ua_national.zip

tl_2018_us_county.zip

Datenquellen: [NOAACPC], [USGAZETTEER] [TIGER]

Arbeitsablauf

1. Wir entpacken die Dateien 2018_Gaz_ua_national.zip und tl_2018_us_county.zip in ein lokales Verzeichnis. Öffne QGIS, gehe im QGIS Browser zur Datei us.tmax_nohads_ll_20190501_float.tif und ziehe sie in den Arbeitsbereich.

2. Wir sehen den neu geladenen Rasterlayer us.tmax_nohads_ll_20190501_float im Layer Bereich. Die Rasterdatei enthält die maximal erfasste Temperatur für jedes Pixel. Als nächstes laden wir die Punktdatei der Ballungsräume hinzu. Das Format der Datei ist Text mit dem Trennzeichen Tabulator (TSV). Klicke auf den Knopf Datenquellenverwaltung öffnen in der Datenquellenverwaltungsleiste.

3. Wir gehen zum Abschnitt Getrennte Texte. Hier klicken wir auf den Knopf … neben Dateiname und geben den Pfad zu der heruntergeladenen Textdatei an. Im Bereich Dateiformat wählen wir Benutzerdefiniert und Tabulator. Als X-Feld wählen wir INTPTLONG und als Y-Feld INTPTLAT. Klicke auf Hinzufügen und danach auf Schließen.

4. Der neue Punktlayer 2018_Gaz_ua_national wird nun im Layer Bereich angezeigt. Wir können jetzt die Werte unter den Punkten aus dem Rasterlayer extrahieren. Gehe zu Verarbeitung ‣ Werkzeugkiste.

5. Suche den Algorithmus Rasteranalyse ‣ Rasterwerte abtasten und klicke doppelt darauf, um ihn zu starten.

6. Wir wählen 2018_Gaz_ua_national als Eingabepunktlayer und us.tmax_nohads_ll_20190501_float als Raster Layer to sample. Wir erweitern den Bereich Fortgeschrittene Parameter und geben tmax als Ausgabespaltenpräfix an. Klicke auf Starte und nach Beendigung der Prozessierung auf Schließen.

7. Ein neuer Layer Abgetastete Punkte wird nun im Layer Bereich angezeigt. Wähle das Tool Objekte abfragen aus der Attributwerkzeugleiste und klicke auf einen beliebigen Punkt. Die Attribute des Punktes werden im Fenster Abfrageergebnisse angezeigt. Wir sehen ein neues Attribut tmax_1, das zu jedem Feature hinzugefügt wurde. Dabei handelt es sich um den Wert des Pixels des Rasterlayers unter dem Punkt. Die 1 steht für die Nummer des Bandes des Rasterlayers. Wenn der Rasterlayer mehrere Bänder enthalten würde, würden entsprechend viele neue Attribute im Ausgabelayer erzeugt werden.

8. Der erste Teil unserer Analyse ist abgeschlossen. Wir können die nicht mehr benötigten Layer entfernen. Halte die Taste Shift gedrückt und wähle die Layer Abgetastete Punkte und 2018_Gaz_ua_national. Klicke mit der rechten Maustaste auf die Layer und wähle Layer löschen, um die Layer aus QGIS zu entfernen. Bei der Nachfrage 2 Layer entfernen? wählen wir OK.

9. Wir werden jetzt den Layer der Verwaltungsbezirke zur Abstastung der Rasterdaten verwenden und für jeden Bezirk einen Mittelwert der Temperatur bestimmen. Wir gehen im QGIS Browser zur Datei tl_2018_us_county.shp und ziehen sie in den Arbeitsbereich.

Bemerkung

Die meisten Verarbeitungsalgorithmen lesen den Eingabelayer ein und erzeugen einen neuen Ausgabelayer. Der Algorithmus Zonenstatistik arbeitet anders. Er verändert den Eingabelayer und fügt neue Attribute zu ihm hinzu. Daher ist es wichtig, den Eingabelayer vor der Verarbeitung zu entpacken. QGIS kann Layer zwar dirket innerhalb einer zip-Datei einlesen, der Layer kann aber dann nicht verändert werden. Der Verarbeitungsalgorithmus schlägt fehl, wenn der Eingabelayer nicht verändert werden kann.

10. Der Layer tl_2018_us_county wird nun im Layer Bereich angezeigt. Wir gehen zu Verarbeitung ‣ Werkzeugkiste.

11. Wir suchen den Algorithmus Rasteranalyse ‣ Zonenstatistik und klicken doppelt darauf, um ihn zu starten.

12. Wähle us.tmax_nohads_ll_20190501_float als Rasterlayer und tl_2018_us_county als Zonenvektorlayer. Gebe tmax_ als Ausgabespaltenpräfix vor. Klicke auf … neben Zu berechnende Statistik.

13. Wir wählen nur Mittel und klicken auf OK.

14. Klicke zum Starten der Berechnung auf Starte. Die Verarbeitung des Algorithmus kann einige Minuten dauern. Klicke anschließend auf Schließen.

15. Wie vorher beschrieben erstellt der Algorithmus Zonenstatistik keinen neuen Layer. Er modifiziert statt dessen den Zonenvektorlayer. Klicke mit der rechten Maustaste auf den Layer tl_2018_us_county und wähle Attributtabelle öffnen.

16. Wir sehen die neu hinzugefügte Attributspalte tmax_mean in der Tabelle. Sie enthält die mittlere Temperatur für jedes Polygon über dem Rasterlayer. Es gibt auch einige leere Werte. Sie entstehen für Polygone die außerhalb des Rasterlayers liegen (Alaska, Hawaii und Puerto Rico).

comments powered by Disqus