Einfache Visualisierung und Analyse von Rasterdaten (QGIS3)

Viele wissenschaftliche Beobachtungen und Forschungen erzeugen Rastendaten. Rasterdaten sind Gitter von Pixeln mit bestimmten Werten. Mit Hilfe von mathematischen Operationen auf den Pixelwerten kann man interessante Analysen durchführen. QGIS hat einige einfache Analysemöglichkeiten im Rasterrechner integriert. In dieser Anleitung werden wir die Möglichkeiten zur Darstellung von Rasterdaten und die Funktionalität des Rasterrechners erkunden.

Überblick über die Aufgabe

Wir werden einen Gitterdatensatz der weltweiten Bevölkerung benutzen, um eine thematische Karte der Bevölkerungsentwicklung zwischen den Jahren 2000 und 2010 zu erstellen.

Weitere Fähigkeiten die wir erlernen

  • Wie man einen Layerstil zwischen Layern kopiert

Beschaffung der Daten

Wir werden folgenden Datensatz der Columbia University benutzen: Gridded Population of the World (GPW) v4. Im Detail brauchen wir Population Count für die gesamte Erde mit einer Auflösung von 2,5 geografischen Minuten im GeoTIFF Format für die Jahre 2000 und 2010. Man benötigt einen kostenloses Nutzerkonto, um die Daten herunterzuladen: Earth Data account.

../../_images/data117.png

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

gpw-v4-population-count-rev11_2000_2pt5_min_tif.zip

gpw-v4-population-count-rev11_2010_2pt5_min_tif.zip

Datenquelle [GPW4]

Arbeitsablauf

  1. Wir öffenen QGIS und gehen im Browser zu den heruntergeladenen Dateien. Erweitere die gpw-v4-population-count-rev11_2000_2pt5_min_tif.zip Datei und ziehe die Datei gpw-v4-population-count-rev11_2000_2pt5_min.tif in den Arbeitsbereich.

../../_images/1140.png

  1. Der neue Layer gpw-v4-population-count-rev11_2000_2pt5_min wird nun im Bereich Layer angezeigt. Analog verfahren wir mit der Datei gpw-v4-population-count-rev11_2010_2pt5_min_tif.zip. Wir erweitern sie und ziehen die Datei gpw-v4-population-count-rev11_2010_2pt5_min.tif in den Arbeitsbereich.

../../_images/2108.png

  1. Wir schauen uns die Layer näher an. Klicke auf den Knopf Objekte abfragen in der Attributwerkzeugleiste. Nach Auswahl des Tools können wir in einen beliebigen Bereich der Arbeitsfläche klicken.

../../_images/359.png

  1. Der zu dem jeweiligen Pixel gehörende Wert wird in einem neuen Bereich Identifikationsergebnis angezeigt. Wir ändern im Identifikationsergebnis den Modus zu Von oben nach unten. Damit werden die Pixelwerte aller Rasterlayer angezeigt und nicht nur die des obersten Layers. Wir können die Werte der beiden Layer vergleichen. Die Auflösung der Rasterzellen beträgt in etwa 5 km x 5 km. Die Pixelwerte stellen die Bevölkerungssumme bezogen auf die Fläche des Pixels dar (25 km²).

../../_images/429.png

  1. Wir schließen das Identifikationsergebnis. Wir werden jetzt eine bessere Darstellung der Layer erzeugen. Klicke auf den Knopf Layergestaltungsfenster öffnen im Bereich Layer.

../../_images/527.png

  1. Im Fenster Layergestaltung wählen wir im Dropdown-Menü neben Symbolisierung die Darstellung Einkanalpseudofarbe.

../../_images/626.png

  1. Diese Darstellung stellt den Layer mit Hilfe eines Farbverlaufes dar. Der voreingestellte Farbverlauf ist weiß zu blau, wobei dem minimalen Wert des Layers weiß und dem maximalen Wert blau zugeordnet ist. Die dazwischen liegenden Werte werden als Farbverlauf verschiedener Blautöne dargestellt. Erweitere die Min/Max-Werteinstellungen und wähle die Option Kumulativer Zählschritt. Die Kartendarstellung sieht nun viel besser aus. Der vorgegebene Wertebereich umfasst nur die Werte zwischen 2 % und 98 % der Datenwerte. Die extremen Werte werden nicht angezeigt. Die Darstellung ist damit viel ansprechender.

../../_images/726.png

  1. Wir schließen die Layergestaltung. Eine ähnliche Darstellung können wir für den anderen Layer einstellen. Es gibt aber einen einfacheren Weg, um die Gestaltung von einem Layer auf einen anderen zu übertragen. Klicke mit der rechten Maustaste auf den Layer gpw-v4-population-count-rev11_2010_2pt5_min und wähle Stile ‣ Stil kopieren.

../../_images/825.png

  1. Klicke jetzt auf den noch unveränderten Layer gpw-v4-population-count-rev11_2000_2pt5_min und wähle Stile ‣ Stil einfügen.

../../_images/926.png

  1. Dieselben Darstellungsparameter werden nun auf den Layer angewendet. Diese Funktionalität ist besonders nützlich, wenn man verschiedene Layer anhand derselben Kategorisierung vergleichen möchte. Durch Umschalten der Sichtbarkeit des obersten Layers, sieht man die Veränderung der Bevölkerungszahlen.

../../_images/101.gif

  1. Unsere Aufgabe ist die Erstellung einer thematischen Karte zur Änderung der Bevölkerung. Wir werden nun den Unterschied zwischen den 2 Layern berechnen und daraus einen weiteren Rasterlayer erstellen. Der neue Rasterlayer enthält für jeden Pixel die Bevölkerungsdifferenz. Gehe zu Raster ‣ Rasterrechner.

../../_images/101.gif

  1. Im Bereich Rasterkanäle kann man die Layer durch Doppelklick auswählen. Das Band wird hinter dem Layernamen und einem @ angezeigt. Da unsere beiden Rasterlayer nur ein Band enthalten, wird bei beiden Layern @1 angehängt. Der Rasterrechner führt mathematische Operationen auf den Rasterpixeln aus. In diesem Fall wollen wir eine einfache Formel verwenden, um die Bevölkerungszahl 2000 von der Bevölkerungszahl 2010 abzuziehen. Wir geben den folgenden Ausdruck ein. Als nächstes drücken wir auf den Button … neben Ausgabelayer.

"gpw-v4-population-count-rev11_2010_2pt5_min@1" - "gpw-v4-population-count-rev11_2000_2pt5_min@1"
../../_images/1235.png

  1. Gib population_change_2010_2000.tif als Ausgabelayer ein und klicke auf OK, um die Berechnung zu starten.

../../_images/1333.png

  1. Nach Beendigung der Berechnung wird der neue Layer population_change_2010_2000 im Layer-Panel hinzugefügt. Wir werden jetzt die Darstellung anpassen, so dass man die Änderung der Bevölkerung besser erkennen kann. Klicke auf den Button Layergestaltungsfenster öffnen im Bereich Layer.

../../_images/1430.png

  1. Eine Möglichkeit besteht darin, die gleiche Gestaltungstechnik wie zuvor zu verwenden und einen divergierenden Farbverlauf zu wählen. Klicke auf das Dropdown-Menü Farbverlauf und wähle Spectral. Klicke erneut auf das Dropdown-Menü und wähle Farbverlauf invertieren, um niedrigen Werten Blautöne und hohen Werten Rottöne zuzuordnen.

../../_images/1528.png

  1. Das ist eine gute Darstellung, die aber nicht leicht zu interpretieren ist. Wir erzeugen jetzt eine bessere Darstellung mit 4 Kategorien: Abnahme, Unverändert, Wachstum und Starkes Wachstum. Scrolle nach unten zur Tabelle mit der Klasseneinteilung. Halte die Shift-Taste gedrückt und selektiere alle Zeilen. Klicke auf den Button Gewählte Zeile(n) entfernen.

../../_images/1626.png

  1. Ändere den Modus bei Interpolation zu Diskret. Wir werden jetzt die Farben manuell zuweisen. Klicke auf den Button Werte manuell hinzufügen. Gib -100 als Wert und Abnahme als Beschriftung ein. Weise der Klasse blau als Farbe zu. Alle Werte kleiner als der verwendete Wert werden in der vorgegebenen Farbe der Klasse angezeigt. Wir bemerken, dass nur noch die Gebiete mit negativer Bevölkerungsentwicklung in der Karte angezeigt werden.

../../_images/1727.png

  1. Vervollständige die Farbkarte mit geeigneten Werten. Ich wählte 100, 1000 und 100000 als Obergrenzen für die Kategorien Neutral, Wachstum und Hohes Wachstum. Weise jeder der erstellten Kategorien Farben zu, zum Beispiel beige, orange und rot.

../../_images/1824.png

  1. Nach Fertigstellung der gewünschten Darstellung, können wir das Fenster Layergestaltung schließen. Wir haben jetzt unsere globale thematische Karte der Bevölkerungsentwicklung fertig gestellt.

../../_images/1919.png

If you want to give feedback or share your experience with this tutorial, please comment below. (requires GitHub account)