GRASS GIS Integration¶

GRASS Plugin starten¶

Um GRASS-Funktionen in QGIS zu benutzen oder GRASS Raster- und Vektorlayer zu visualisieren, muss das GRASS-Plugin im Plugin-Manager geladen werden. Wechseln Sie dazu in das Menü Erweiterungen ‣ Erweiterungen verwalten, wählen Sie GRASS aus und drücken dann auf [OK]. Nach dem Laden erscheint eine neue Werkzeugleiste mit den oben dargestellten Icons.

Sie können wie in Abschnitt GRASS Layer visualisieren beschrieben, direkt Raster- und Vektorlayer aus einer existierenden GRASS-Datenbank Location laden, oder Sie können eine neue GRASS Location erstellten (siehe Abschnitt Eine neue GRASS LOCATION erstellen), Daten in diese neue Location importieren (siehe Abschnitt import_loc_data) und über die GRASS-Werkzeugkiste mit den mehr als 300 GRASS Modulen analysieren (siehe Abschnitt Die GRASS-Werkzeugkiste).

GRASS Layer visualisieren¶

Wenn das GRASS-Plugin geladen ist, können Sie GRASS Vektor- und Rasterlayer mit den entsprechenden Knöpfen in der Werkzeugleiste laden. Als Beispiel benutzen wir den Alaska Beispieldatensatz (siehe Kapitel Beispieldaten). Dieser enthält eine kleine GRASS LOCATION mit 3 Vektor- und einem Rasterlayer mit Höheninformationen.

1. Erstellen Sie einen neuen Ordner grassdata und laden Sie den QGIS Alaska Beispieldatensatz qgis_sample_data.zip unter http://download.osgeo.org/qgis/data/ herunter und entpacken Sie diesen in den gerade erstellen Ordner.

2. Starten Sie QGIS.

3. Falls noch nicht geschehen, laden Sie das GRASS Plugin. Dazu wechseln Sie in das Menü Erweiterungen ‣ Erweiterungen verwalten und wählen GRASS aus. Die GRASS-Werkzeuge erscheinen nun in der Werkzeugleiste.

4. Klicken Sie nun auf das Icon ^{Mapset öffnen}, um eine Session in einer GRASS Location zu starten.

5. Als Gisdbase suchen und wählen Sie bitte den Pfad zum Ordner grassdata.

6. Sie sollten nun als LOCATION alaska und als MAPSET demo auswählen können.

7. Klicken Sie auf [OK]. Einige weitere, zuvor grau hinterlegte GRASS-Werkzeuge sind nun aktiv.

8. Klicken Sie auf ^{GRASS-Rasterlayer hinzufügen}, wählen Sie den Layer gtopo30 und drücken Sie auf [OK]. Die Höhendaten werden nun dargestellt.

9. Klicken Sie auf toolbtntwo ^{GRASS-Vektorlayer hinzufügen}, wählen Sie den Layer alaska und drücken Sie auf OK. Der Layer mit der Landesgrenze von Alaska wird nun über dem Rasterlayer dargestellt. Sie können nun die Layereigenschaften einstellen, wie in Kapitel vector_properies_dialog beschrieben. Z.B.: die Transparenz verändern oder andere Farben auswählen.

10. Laden Sie auch noch die weiteren Vektorlayer rivers und airports hinzu und passen Sie deren Layereigenschaften an.

Wie man sehen kann, ist es sehr einfach sich Raster- und Vektordaten aus einer existierenden GRASS-Datenbank in QGIS anzeigen zu lassen. Lesen Sie die folgenden Kapitel, um zu lernen, wie man GRASS Daten editiert und eine neue Location erstellt. Weitere Beispiel GRASS-Datenbanken finden Sie auf der GRASS Homepage unter: http://grass.osgeo.org/download/data.php.

Information zur GRASS-Datenbank¶

GRASS Daten werden in einem Ordner gespeichert, der als GISDBASE bezeichnet wird. Standardmäßig wird der Ordner grassdata genannt und er muss erstellt worden sein, bevor man beginnt, mit dem GRASS Plugin in QGIS zu arbeiten. Innerhalb dieses Ordners sind die GRASS Daten als Projekte (sog. LOCATION) in Unterordnern organisiert. Jede LOCATION ist durch ein Koordinatenbezugssystem und eine räumliche Grenze (sog. region) definiert und kann darüberhinaus weitere Unterordner MAPSETs besitzen, um die Layer der LOCATION weiter z.B. thematisch oder räumlich zu unterteilen (Neteler & Mitasova 2008 literature_and _web). Um Raster- und Vektorlayer mit den GRASS Modulen zu analysieren, müssen diese zuerst in eine passende GRASS LOCATION importiert werden. (Dies ist nicht ganz korrekt. Mit den GRASS Modulen r.external und v.external können Sie eine ‘read-only’ Verknüpfung zu externen durch GDAL/OGR-unterstützte Layer erstellen, ohne die Daten importieren zu müssen. Da dies aber nicht der normale Weg für GRASS Anfänger ist, wird auf diese Möglichkeit nicht näher eingegangen.).

Figure GRASS location 1:

GRASS data in the alaska LOCATION

Eine neue GRASS LOCATION erstellen¶

Als ein Beispiel möchten wir Ihnen zeigen, wie die GRASS Location des Alaska Beispieldatensatzes erstellt wurde. Das Koordinatenbezugssystem ist Albers Equal Area mit der Einheit ‘feet’. Diese GRASS Location alaska wird für alle GRASS GIS Beispiele verwendet. Es ist also sinnvoll, sich diesen Datensatz zu installieren (siehe Abschnitt Beispieldaten).

1. Starten Sie QGIS und laden Sie das GRASS Plugin, falls dies noch nicht geschehen ist.

2. Visualisieren Sie das Shape alaska.shp (vgl. Kapitel Shape Layer laden) aus dem QGIS Beispieldatensatz.

3. Klicken Sie auf das Icon ^{Neues Mapset} in der GRASS-Werkzeugleiste, um den MAPSET Dialog zu starten.

4. Wählen Sie einen existierenden Ordner mit bereits vorhandenen Locations oder erstellen Sie einen neuen Ordner grassdata für die zu erstellende Location. Klicken Sie auf [Weiter].

5. Mit dem Dialog kann eine bereits vorhandene LOCATION durch eine weitere MAPSET ergänzt werden, oder eine neue LOCATION mit neuen MAPSETS erstellt werden (siehe Abschnitt Eine neue GRASS MAPSET erstellen). Klicken Sie nun auf den Knopf Erstelle neue Location wie in Abbildung figure_grass_location_2 zu sehen.

6. Geben Sie einen Namen für die neue LOCATION ein - wir haben den Namen alaska gewählt. Klicken Sie auf [Weiter].

7. Definieren Sie das Koordinatenbezugssystem, in dem Sie auf den Knopf Projektion klicken, um aus einer Liste auswählbare Koordinatensysteme angezeigt zu bekommen.

8. Wir verwenden die Albers Equal Area Alaska (feet) Projektion. Da wir wissen, dass es sich hierbei um den EPSG Code 2964 handelt, geben wir die Zahl in das Suchfenster ein. (Wenn Sie dies für einen anderen Layer wiederholen wollen, aber nicht den EPSG Code kennen, klicken Sie einfach auf das Icon ^KBS-Status unten rechts in der Statusleiste (vgl. Kapitel label_projstart)).

9. Geben Sie unter Filter 2964 ein um die Projektion auszuwählen.

10. Klicken Sie auf [Weiter].

11. Um die Standardausdehnung der Location zu definieren, müssen wir die Eckkoordinaten der Umrandung angeben. Da wir einen Layer geladen haben, klicken wir einfach auf den Knopf [Setze aktuelle QGIS-Ausdehung], um die Ausdehnung der Karte als Standardausdehnung für die LOCATION zu verwenden.

12. Klicken Sie auf [Weiter].

13. Wir müssen nun noch eine MAPSET für die LOCATION definieren. Wir haben uns für demo entschieden. Beim Erstellen einer neuen LOCATION, wird automatisch eine spezielle MAPSET mit dem Namen PERMANENT eingerichtet, um Kerninformationen wie etwa die räumliche Ausdehung und das Koordinatenbezugssystem der Location zu sichern. Diese wird zusätzlich zu der von Ihnen gewählten Mapset erstellt (Neteler & Mitasova 2008 literature_and _web)

14. Prüfen Sie alles, damit alles wie gewünscht ist. Klicken Sie auf [Abschliessen].

15. Die neue LOCATION alaska mit den beiden MAPSETs PERMANENT und demo sind nun erstellt worden. Ihre aktuelle Arbeitsumgebung ist die von Ihnen definierte MAPSET demo.

16. Beachten Sie, dass einige Werkzeuge des GRASS Plugins grau hinterlegt waren und nun auch zur Verfügung stehen.

Figure GRASS location 2:

Creating a new GRASS LOCATION or a new MAPSET in QGIS

Auch wenn es recht kompliziert aussieht, ist dies doch ein sehr einfacher Weg, eine neue GRASS Location zu erstellen. Unsere Location alaska ist nun bereit für den Import von Daten (siehe Abschnitt Daten in eine GRASS LOCATION importieren). Sie können an dieser Stelle aber auch die vorhandene LOCATION alaska des QGIS Beispieldatensatzes verwenden und direkt zu Kapitel Das GRASS Vektormodell weitergehen.

Daten in eine GRASS LOCATION importieren¶

In diesem Abschnitt soll an einem Beispiel gezeigt werden, wie man Raster- und Vektorlayer aus dem QGIS Beispieldatensatz in die LOCATION alaska importieren kann. Zu diesem Zweck verwenden wir den Rasterlayer landcover.img und die Vektor GML-Datei lakes.shp aus dem Beispieldatensatz Beispieldaten.

1. Starten Sie QGIS und laden Sie das GRASS Plugin, falls dies noch nicht geschehen ist.

2. Klicken Sie auf das Icon ^{Mapset öffnen} in der GRASS Werkzeugleiste, um den MAPSET Dialog zu starten.

3. Wählen Sie als Gisdbase den Ordner grassdata aus dem QGIS Alaskadatensatz, als LOCATION alaska, als Kartenset demo und klicken Sie dann auf [OK].

4. Nun klicken Sie auf das Icon ^{GRASS-Werkzeugkiste öffnen}, damit die GRASS Werkzeuge (siehe Abschnitt Die GRASS-Werkzeugkiste) zur Verfügung stehen.

5. Um den Rasterlayer landcover.img zu importieren, drücken Sie auf das Modul r.in.gdal im Modulbaum Reiter. Diese GRASS Module ermöglicht es, GDAL-unterstützte Rasterlayer in eine GRASS LOCATION zu importieren.

6. Wählen Sie als zu importierende Rasterdatei den Layer landcover.img aus dem Beispieldatensatz.

7. Als Name der Ausgabe wählen Sie landcover_grass und drücken Sie dann auf [Starten]. Im Reiter Ergebnis sehen Sie das gerade gestartete GRASS Kommando r.in.gdal -o input=/path/to/landcover.img output=landcover_grass.

8. Wenn Sie den Text Erfolgreich beendet sehen, klicken Sie auf den Knopf [Ergebnis visualisieren]. Der Layer landcover_grass ist nun in die aktuelle GRASS Location importiert und wird im Kartenfenster angezeigt.

9. Um den Vektor GML-Layer lakes.gml zu importieren, klicken Sie in der GRASS Werkzeugkiste auf das Modul v.in.ogr im Modulbaum Reiter. Dieses GRASS Modul ermöglicht es, OGR-unterstützte Vektorlayer in eine GRASS LOCATION zu importieren. Der Moduldialog für v.in.ogr erscheint als neuer Reiter.

10. Wählen Sie als zu importierende Vektordatei den Layer lakes.gml aus dem Beispieldatensatz.

11. Als Name der Ausgabe wählen Sie lakes_grass und drücken Sie dann auf [Starten]. Im Reiter Ergebnis sehen Sie das gerade gestartete GRASS Kommando v.in.ogr -o dsn=/path/to/lakes.shp output=lakes_grass.

12. Wenn Sie den Text Erfolgreich beendet sehen, klicken Sie auf den Knopf [Ergebnis visualisieren]. Der Layer lakes_grass ist nun in die aktuelle GRASS Location importiert und wird im Kartenfenster angezeigt.

Das GRASS Vektormodell¶

Für das Arbeiten mit GRASS-Vektordaten ist es wichtig, das Datenmodell von GRASS-Vektoren zu kennen.

GRASS nutzt ein topologisches Datenmodell.

Das bedeutet, dass Flächen nicht als geschlossene Polygone vorhanden sind, sondern als ein oder mehrere Umrandungen (Boundaries). Eine Umrandung (Boundary) zwischen zwei aneinander grenzenden Flächen ist nur einmal digitalisiert worden; beide Flächen teilen sich diese Umrandung. Umrandungen dürfen keine Lücken haben. Eine Fläche besteht also aus einer Umrandung und einem Zentroid, der diese Fläche als ein sog. Labelpunkt mit einer Attributtabelle verknüpft.

Neben den Umrandungen und Zentroiden kann eine Vektorkarte selbstverständlich auch Punkte und Linien enthalten. Alle diese Geometrieelemente können innerhalb ein und dem selben Datensatz enthalten sein. Sie werden in unterschiedlichen ‘Ebenen’ innerhalb von QGIS dargestellt. Auch wenn es möglich ist, Geometrieelemente zu mischen, so ist es eigentlich unüblich und wird normalerweise auch in GRASS GIS nur selten verwendet. Etwa bei Netzwerkanalysen. Im Normalfall sollten Sie versuchen, unterschiedliche Geometrietypen in unterschiedlichen Datensätzen (Layern) zu speichern.

Es ist auch möglich, unterschiedliche Inhalte des gleichen Geometrietyps in verschiendenen Ebenen eines Vektorlayers zu speichern. Beispielsweise können Felder, Wälder und Seen in einem Vektordatensatz gespeichert werden. Angrenzende Seen, Felder und Wälder teilen sich dann die gleiche Umrandung, jedoch haben sie separate Attributtabellen, die über ihre Ebene angesprochen wird. Darüber hinaus können Sie auch Attribute für die Umrandungen vergeben, falls eine Umrandung gleichzeitg einen Weg darstellt. In diesem Fall könnte auch die Umrandung eine separate Attributtabelle haben.

Die ‘Ebene’ eines jeden Objektes wird in GRASS intern als ‘layer’ bezeichnet. Verwechseln Sie dies bitte nicht mit dem Begriff Layer im Sinne einer Karte. Der GRASS Layer (Ebene) wird durch Zahl repräsentiert. Wenn mehr als ein Layer im Datensatz enthalten ist; beim obigen Beispiel wären es 3 Layer (Felder, Wälder und Seen), gäbe es drei Layer (1,2 und 3). Derzeit werden nur Nummern als Layername unterstützt, in kommenden GRASS-Versionen werden auch Namen möglich sein.

Attribute können in externen Datenbanktabellen abgelegt werden, beispielsweise DBase, PostgreSQL, MySQL, SQLITE3, etc.

Die Attribute in den Tabellen werden über ein sog. ‘Kategoriefeld’ an die Geometrien des Datensatzes gehängt.

Die ‘Kategorie’ (oder key, ID, etc) ist eine Ganzzahl, über die eine Verknüpfung zwischen den Geometrien und den Spalten in der Datenbanktabelle hergestellt wird.

Einen neuen GRASS Vektorlayer erstellen¶

Einen neuen GRASS Vektorlayer können Sie erstellen, indem Sie auf das Icon ^{Neuen GRASS-Vektorlayer anlegen} in der GRASS Werkzeugleiste klicken. Geben Sie einen Namen in der Textbox für den neuen GRASS Vektorlayer an und Sie können in dem Layer neue Geometrien erzeugen, wie in Abschnitt Digitalisieren und Editieren eines GRASS Vektorlayers beschrieben.

GRASS GIS erlaubt es aufgrund des topologischen Datenmodells, die verschiedenen Geometrietypen (Punkt, Linie und Fläche) in einem Vektorlayer abzuspeichern. Aus diesem Grund ist es nicht notwendig im Vorfeld einen Geometrietyp festzulegen. Dies unterscheidet sich von der Erstellung eines neuen Shapefile in QGIS, denn Shapefiles verwenden das Simple Feature Vektormodell (siehe Abschnitt Neuen Vektorlayer erstellen).

Digitalisieren und Editieren eines GRASS Vektorlayers¶

Die Digitalisierwerkzeuge für GRASS-Vektorlayer werden über den Knopf ^{GRASS-Vektorlayer bearbeiten} in der GRASS-Werkzeugleiste gestartet. Dazu müssen Sie den zu bearbeitenden Layer in der Legende auswählen, bevor Sie auf den Editier-Knopf drücken. Abbildung figure_grass_digitizing_2 zeigt den Dialog für die GRASS Digitalisierung. Die einzelnen Werkzeuge werden im folgenden Kapitel beschrieben.

Werkzeugleiste

Abbildung figure_grass_digitizing_1 zeigt die GRASS Digitalisierwerkzeuge. Tabelle table_grass_digitizing_1 beschreibt die damit zur Verfügung stehenden Funktionalitäten.

Figure GRASS digitizing 1:

GRASS Digitizing Toolbar

Icon	Werkzeug	Zweck
	Neuer Punkt	Digitalisiert neuen Punkt
	Neue Linie	Digitalisiert neue Linie
	Neue Grenze	Digitalisiert neue Grenze (zum Beenden ein neues Werkzeug wählen)
	Neuer Zentroid	Digitalisiert neuen Zentroiden (Labelpunkt für eine existierende Fläche)
	Verschiebe Vertex	Wählt einen Stützpunkt einer existierenden Linie oder Umrandung und setzt diesen an eine neue Position
	Vertex hinzufügen	Fügt einen Stützpunkt zu einer existierenden Linie hinzu
	Lösche Vertex	Löscht einen Stützpunkt von einer existierenden Linie (Bestätigung der Auswahl durch einen weiteren Klick nötig)
	Verschiebe Element	Wählt existierende Grenze, Linie oder Zentroiden und verschiebt sie an eine neue Position
	Unterteile Linie	Teilt eine existierende Linie in zwei Teile
	Element löschen	Löscht eine existierende Geometrie inklusive des Eintrags in der Attributtabelle (Bestätigung der Auswahl durch einen weiteren Klick nötig).
	Editiere Attribute	Editiert Attribute eines existierenden Objekts (Beachten Sie, dass ein Objekt mehrere Feature repräsentieren kann, siehe oben)
	Schließen	Beendet die Bearbeitung (und aktualisiert die Topolgie anschließend)

Tabelle GRASS Digitizing 1: GRASS Digitalisierwerkzeuge

Reiter Kategorie

In dem Reiter Kategorie können Sie einstellen, in welcher Weise Kategoriewerte neuen Objekten oder neue Kategoriewerte vorhandenen Objekten zugewiesen werden sollen.

Figure GRASS digitizing 2:

GRASS Digitizing Category Tab

• Modus: Welcher Kategoriewert soll der Geometrie zugewiesen werden.

  • Nächst folgender Kategoriewert - der nächste Kategoriewert, der noch nicht im Vektordatensatz benutzt wird.

  • Manueller Eintrag - Manuelles Eintragen der Kategoriewerte im ‘Kategorie’-Eingabgefeld vornehmen.

  • Keine Kategorie - digitalisiere Geometrie ohne Angabe eines Kategoriewertes. Dies wird z.B. beim Digitalisieren einer Grenzlinie eines Polygons verwendet.

• Kategorie - eine Nummer (ID), die dem digitalisierten Objekt zugewiesen wird und auf einen Eintrag in der Attributtabelle verweist.

• Layer - Objektidentifikation über die GRASS-Ebene (layer). Der Standardlayer in GRASS ist 1.

Reiter Einstellungen

Der Reiter Einstellungen erlaubt das Setzen der Fangtoleranz in Bildschirmpixeln. Dies ist der Schwellenwert in Pixeln, innerhalb dessen neu digitalisierte Knotenpunkte an vorhandene Knoten gesnappt werden. Dies hilft, Lücken oder Überlagerungen zwischen Objekten zu vermeiden. Der Standardwert ist auf 10 Pixel eingestellt.

Figure GRASS digitizing 3:

GRASS Digitizing Settings Tab

Darstellung Reiter

Der Reiter Darstellung erlaubt die Farbeinstellungen für die verschiedenen Geometrietypen und ihren Topologiestatus (z.B. offene/geschlossene Fläche).

Figure GRASS digitizing 4:

GRASS Digitizing Symbolog Tab

Reiter Tabelle

Der Reiter Tabelle gibt Informationen über die Attributtabelle des aktuellen Layers. In diesem Reiter können Sie eine Attributtabelle hinzufügen, modifizieren oder eine neue Tabelle erstellen (siehe Abschnitt Einen neuen GRASS Vektorlayer erstellen).

Figure GRASS digitizing 5:

GRASS Digitizing Table Tab

Einstellung der GRASS Region¶

Die aktuelle Region (Projektausdehnung und Pixelauflösung) ist in GRASS sehr wichtig für alle Rastermodule. Alle neu erzeugten Rasterdaten haben nämlich die Ausdehnung und Auflösung der aktuellen Region, egal wie groß die Standard-Region definiert ist. Die Region ist in der Datei $LOCATION_NAME/$MAPSET/WIND gespeichert, und definiert Nord-, Süd-, Ost- und Westausdehnung sowie Anzahl der Spalten, Anzahl der Reihen und die horizontale und vertikale räumliche Pixelauflösung.

Es ist auch möglich, die Region mit dem Knopf ^{Aktuelle GRASS-Region darstellen} aus- bzw. einzusschalten.

Um die aktuelle GRASS-Region zu ändern, benutzen Sie den Knopf ^{Aktuelle GRASS-Region editieren}. Dort haben Sie die Möglichkeit, die Region zu verändern sowie das Aussehen der Regionsdarstellung. Wenn der Dialog sichtbar ist, können Sie die Region neben der einfachen Koordinateneingabe auch interaktiv mit der Maus im Kartenbild verändern. Die Koordinaten werden dann automatisch angepasst, sobald Sie den Dialog mit dem Knopf [OK] schließen.

Das GRASS Modul g.region bietet viele weitere Optionen zur Einstellung der passenden Ausdehnung und Auflösung der Region für die Rasteranalyse. Das Modul g.region können Sie über die GRASS Werkzeugkiste nutzen, wie in Kapitel Die GRASS-Werkzeugkiste beschrieben.

Die GRASS-Werkzeugkiste¶

Die ^{GRASS-Werkzeugkiste} ermöglicht es, GRASS Module auf Daten innerhalb einer ausgewählten GRASS LOCATION und MAPSET anzuwenden. Dazu muss im Vorfeld eine GRASS LOCATION und MAPSET geöffnet werden, in der Sie Schreibrechte besitzen. Dies ist normalerweise garantiert, wenn Sie die MAPSET selbst erstellt haben und notwendig, damit die Ergebniskarten der Raster- und Vektoranalysen in der ausgewählten MAPSET gespeichert werden können.

Figure GRASS toolbox 1:

GRASS Toolbox and Module Tree

Die GRASS Shell der Werkzeugkiste bietet Zugriff auf fast alle (mehr als 330) GRASS Module über die Kommandozeile. Um eine benutzfreundliche Umgebung zu bieten, sind davon etwa 200 Module graphisch auswählbar und bieten einen Dialog in Form eines zusätzlichen Reiters in der Werkzeugkiste.

Arbeiten mit GRASS Modulen¶

Das GRASS Shell Modul in der Werkzeugkiste bietet Zugriff auf fast alle in GRASS GIS vorhandenen Module (mehr als 330). Um die Verwendung benutzerfreundlicher zu gestalten, sind etwa 200 Module über grafische Dialoge in die Werkzeugkiste integriert. Diese Module sind thematisch untergliedert, können aber auch durchsucht werden.

Eine vollständige Liste vorhandener GRASS Module, die man über das QGIS GRASS Plugin ansprechen kann finden Sie im GRASS Wiki unter der URL : http://grass.osgeo.org/wiki/GRASS-QGIS_relevant_module_list.

Es ist außerdem möglich, die GRASS Werkzeugkiste anzupassen und weitere Module zu integrieren. Die Herangehensweise ist in Abschnitt Anpassen der Module beschrieben.

Wie in Abbildung figure_grass_toolbox_1 zu sehen, können Sie nach dem passenden GRASS Modul in dem Reiter Modulbaum nachschauen oder im Reiter Modulliste suchen.

Wenn Sie auf das grafische Icon eines Modules klicken, öffnet sich ein neuer Moduldialog mit drei Reitern Optionen, Ergebnis und Handbuch.

Optionen

Der Reiter Optionen stellt Ihnen in vereinfachter Form die unbedingt notwendigen Eingabeparameter zur Verfügung, die das Modul zum Laufen benötigt.

Figure GRASS module 1:

GRASS Toolbox Module Options

Bitte beachten Sie, dass diese Eingabefelder extra so einfach wie möglich gehalten wurden, um die Struktur klar zu halten. Falls Sie alle vom Modul unterstützten Parameter benötigen, nutzen Sie die GRASS-Eingabeshell, in der Sie das Modul ebenfalls aufrufen können.

Eine neue Funktion seit der QGIS Version 1.8.0 ist der Knopf Fortgeschrittene Optionen einblenden unterhalb der vereinfachten Moduldialoge im Reiter Optionen. Momentan ist diese Funktion nur für das Modul v.in.ascii umgesetzt. Es ist aber zukünftig für weitere Module vorgesehen und wird so die kompletten Funktionen der GRASS Module grafisch bereitstellen, ohne die GRASS Shell und somit die Kommandozeile benutzen zu müssen.

Ergebnis

Figure GRASS module 2:

GRASS Toolbox Module Output

Der Reiter Ergebnis stellt die Ausgabe des Moduls zur Laufzeit dar. Nachdem Sie den Knopf [Los] gedrückt haben, wird auf diesen Reiter gewechselt und Sie sehen die Statusausgaben des Moduls. Wenn alles funktioniert hat, sehen Sie den Ausgabetext Erfolgreich beendet.

Handbuch

Figure GRASS module 3:

GRASS Toolbox Module Manual

Der Reiter Handbuch enthält die Hilfeseite zu jedem GRASS-Modul. Falls Sie weitergehende Informationen zum gewählten Modul benötigen, schauen Sie in diese Hilfeseite. Dabei werden Sie feststellen, dass viele Module mehr Parameter haben, als im Reiter Optionen zu sehen sind. Das ist gewollt. Um die GUI möglichst einfach zu halten, sind nur die nötigsten Optionen übernommen worden. Natürlich stehen Ihnen alle Optionen des Moduls direkt in der GRASS-Eingabeshell zur Verfügung. Am Ende jeder Hilfeseite finden Sie weitere Links zum Main Help index, dem Thematic index und dem Full index. Diese enthalten alle Informationen, die Sie ansonsten über das Modul g.manual finden.

Tipp

Ergebnisse direkt anzeigen

Wollen Sie Ihre Ergenisse direkt in der Kartenansicht ansehen, nutzen Sie den Knopf ‘Ergebnis visualisieren’ im unteren Bereich des jeweiligen Modulreiters

GRASS Beispielanwendung¶

Die folgenden Beispiele sollen die Anwendung verschiedener GRASS Module demonstrieren.

Höhenlinien aus einem DGM erstellen¶

Im ersten Beispiel sollen automatisiert Höhenlinien auf Basis eines Höhenmodells (DGM) erstellt werden. Dabei gehen wir davon aus, dass Sie die Alaska LOCATION auf dem Rechner installiert ist, wie in Abschnitt Daten in eine GRASS LOCATION importieren beschrieben.

• Als erstes öffnen Sie die Location, indem Sie auf das ^{Mapset öffnen} Icon klicken und dann die Alaska Location auswählen.

• Nun laden Sie den gtopo30 Rasterlayer, indem Sie auf das Icon ^{GRASS-Rasterlayer hinzufügen} klicken und den Layer gtopo30 aus der Alaska Location auswählen und laden.

• Als nächstes drücken Sie auf das Icon ^{GRASS-Werkzeugkiste öffnen}.

• Im Reiter Modulbaum wechseln Sie mit der Maus in den Bereich Raster ‣ Oberflächenverwaltung ‣ Vektorkonturlinien erzeugen.

• Mit einem einfachen Klick auf r.contour öffnet sich das Modul in einem neuen Reiter, wie unter Arbeiten mit GRASS Modulen beschrieben. Der gtopo30 Rasterlayer sollte nun automatisch im Feld Name der Eingaberasterkarte erscheinen.

• Geben Sie im Feld Abstand zwischen den Kontourintervallen den Wert 100 an. (Dadurch werden Höhenlinien in einem Abstand von 100m erstellt.)

• In das Feld Name der Vektorausgabekarte geben Sie den Namen hoehen_100 an.

• Drücken Sie nun auf den Knopf [Starten]. Nach einiger Zeit sollte im Ergebnisreiter die Meldung Erfolgreich beendet erscheinen. Drücken Sie nun auf den Knopf [Ergebnis visualisieren] und dann auf [Schließen].

Da die aktuelle GRASS Region ziemlich groß ist, kann es eine Weile dauern, bis der Layer vollständig dargestellt wird. Danach können Sie noch die Layereigenschaften verändern und eine Farbe für die Linien auswählen, die sich deutlich vom Höhenmodell unterscheidet, siehe Vektorlayereigenschaften.

Als nächstes zoomen Sie in einen bergigen Bereich im Zentrum Alaskas. Wenn Sie weit genug in die Karte hineingezoomt sind, werden Sie erkennen, dass die Höhenlinien teilweise sehr eckig erscheinen. Um das Erscheinungsbild zu optimieren, bietet GRASS ein Modul mit dem Namen v.generalize. Dabei wird mit Hilfe des Douglas Peuker Algorithmus und der einer Reduktion der Stützpunkte eine Glättung der Linien erreicht, ohne die Geometrien zu zerstören. Da der Ergebnislayer weniger Stützpunkte hat, ist er auch kleiner und kann schneller geladen werden. Die Analyse wird z.B. angewendet, wenn man sehr detailierte Daten nur in einem kleinen Maßstab anzeigen möchte.

Tipp

Geometrien in QGIS vereinfachen

QGIS stellt mit den fTools auch ein Werkzeug Geometrien vereinfachen zur Verfügung. Es funktioniert genau wie der Douglas-Peuker Algorithmus im GRASS Modul v.generalize.

In diesem Beispiel wollen wir nun aber etwas anderes erreichen. Die Höhenlinien, die wir mit r.contour erstellt haben, zeigen teilweise sehr scharfe Winkel, die wir glätten möchten. Unter den Algorithmen des Moduls v.generalize befindet sich auch der Chaikens-Algorithmus (), der exakt das macht, was wir möchten. Achten Sie aber darauf, dass es passieren kann, dass Stützpunkte nicht nur eliminiert sondern auch hinzugefügt werden können. Dadurch kann der Layer wieder langsamer geladen werden.

• Öffnen Sie die GRASS-Werkzeugkiste, wechseln Sie in den Bereich Vektor ‣ Karte entwickeln ‣ Generalisierung, und klicken dann auf das v.generalize Modul.

• Stellen Sie sicher, dass ‘hoehen_100’ als Vektorlayer in Feld Name der Vektoreingabekarte erscheint.

• Aus der Liste möglicher Algorithmen wählen Sie nun Chaiken’s. Belassen Sie alle weiteren Optionen wie sie sind und gehen Sie zum Ende des Dialogs, wo Sie als Name der Vektorausgabekarte ‘hoehen_100_smooth’ angeben. Drücken Sie nun auf [Starten].

• Der Prozess dauert eine Weile. Sobald Erfolgreich beendet im Reiter Ergebnis erscheint, drücken Sie wieder auf [Ergebnis visualisieren] und dann auf [Schließen].

• Ändern Sie nun auch die Farbe des neuen Layers, damit er sich deutlich von dem Höhenmodell und den zuvor berechneten Höhenlinien abhebt. Sie werden erkennen, dass die Kanten der neuen Höhenlinien wesentlich weicher gezeichnet sind.

Figure GRASS module 4:

GRASS module v.generalize to smooth a vector map

Tipp

Vektorlinien glätten mit dem GRASS Modul v.generalize

Die oben beschriebene Anwendung kann auch in anderen Situationen verwendet werden. Wenn Sie z.B. eine Rasterkarte mit Niederschlagswerten haben, können Sie mit r.contour einen Isohyetallayer erstellen.

Erstellen eines 3D Schummerungseffekts¶

Es gibt verschiedene Methoden, um Höhenlayer anzuzeigen und ihnen einen 3D Schummerungseffekt zu verleihen. Der Gebrauch von Höhenlinien ist eine populäre Methode, die häufig angewendet wird, um topographische Karten zu erstellen. Eine andere Möglichkeit, um einen 3D Effekt zu erzeugen ist, bietet das Hillshading. Der Effekt basiert dabei auf einem Höhenmodell. Dabei wird zuerst die Hangneigung- und richtung der Zellen bestimmt und dann durch die Simulation des Sonnenstandes eine Reflexionswert erzeugt. Dadurch werden der Sonne zugewandte Bereiche aufgehellt und der Sonne abgewandte Bereiche (im Schatten) dunkler dargestellt.

• Laden Sie zuerst den Rasterlayer gtopo30. Öffnen Sie die GRASS-Werkzeugkiste, wechseln Sie in den Bereich Räumliche Analysen ‣Geländeanalyse.

• Nun klicken Sie auf r.shaded.relief, um den Modulreiter zu öffnen.

• Ändern Sie den Wert im Feld Winkel der Sonne in Grad östlich von der Nordrichtung 270 auf 315.

• Geben Sie als Name der Schummerungskarte den Namen gtopo30_shade an und klicken Sie auf [Starten].

• Nachdem die Karte berechnet wurde, visualisieren Sie sie und setzen Sie die Farbe des Rasterlayers auf Graustufen.

• Um die Schummerung und das Höhenmodell gtopo30 zu sehen, ziehen Sie die Schummerungskarte in der Legende unter das Höhenmodell. Öffnen Sie dann den Dialog Layereigenschaften der Karte gtopo30, und ändern Sie die Transparenz im Reiter Transparenz auf den Wert 25%.

Die Höhenkarte gtopo30 wird nun als Farbkarte leicht transparent über der Schummerungskarte angezeigt. Dadurch entsteht ein visueller 3D Effekt. Um den Unterschied besser zu erkennen, wechseln Sie über das Kontrollkästchen den Anzeigemodus der Schummerungskarte in der Legende und wieder zurück.

Die GRASS Kommandozeile verwenden

Das GRASS Plugin in QGIS stellt die GRASS Module oftmals in vereinfachter Form und auch nicht vollständig bereit. Es ist also grundsätzlich für Anwender gestaltet, die sich nicht so gut mit GRASS und all seinen Fähigkeiten auskennen oder nur einfacher Analysen durchführen wollen. Daher werden in den grafischen Moduldialogen oftmals nicht alle Optionen und Parameter, die das GRASS Modul bieten bereitgestellt, um den Umgang einfacher und intuitiver zu gestalten. Wer tiefer in GRASS einsteigen möchte, er hat die Möglichkeit, sämtliche Funktionalitäten und Module über die GRASS Kommandozeile (GRASS Shell) anzusprechen. In dem folgenden Beispiel soll eine zusätzliche Option des Moduls r.shaded.relief angesprochen werden, die nur über die Kommandozeile genutzt werden kann.

Figure GRASS module 5:

The GRASS shell, r.shaded.relief module

Das Modul r.shaded.relief stellt einen zusätzliche Parameter zmult bereit, über den der Höhenwert relativ zu den X-Y Werten multipliziert werden kann. Dadurch wird der Schummerungseffekt noch prägnanter.

• Laden Sie das Höhenmodell gtopo30 wie im vorherigen Beispiel. Öffnen Sie die GRASS-Werkzeugkiste und klicken auf die GRASS Shell. In das Kommandozeilenfenster tippen Sie folgenden Befehl r.shaded.relief map=gtopo30 shade=gtopo30_shade2 azimuth=315 zmult=3 und drücken dann die Taste Enter.

• Wenn die Berechnung abgeschlossen ist, wechseln Sie den Reiter Browser und Doppelklicken Sie auf die neu erstellte Karte gtopo30_shade2, um Sie in QGIS anzuzeigen.

• Visualisieren Sie die Karten wie in dem obigen Beispiel beschrieben. Sie sollten dabei erkennen, dass der Schummerungseffekt stärker ist, wenn der zmult Parameter verwendet wird.

Figure GRASS module 6:

Displaying shaded relief created with the GRASS module r.shaded.relief

Rasterstatistik auf Basis eines Vektorlayer berechnen¶

Das folgende Beispiel zeigt, wie man univariate Statistik für Rasterwerte innerhalb von Vektorpolygonen berechnen kann und die Werte in neue Attributspalten des Vektorlayers hinzufügt.

• Verwenden Sie wieder die Alaska Location (siehe Abschnitt Daten in eine GRASS LOCATION importieren), und importieren Sie das Shapefile trees aus dem Ordner shapefiles nach GRASS.

• Als erstes müssen dem Layer trees Zentroide (Labelpunkte) hinzugefügt werden. Dadurch wird es ein korrekter topologischer Datensatz, bei dem an die Labelpunkte die Attribute aus der statistischen Analyse angehängt werden können.

• Wählen Sie aus der GRASS-Werkzeugkiste Vektor ‣ Karte entwickeln ‣Objekte verwalten, und öffnen Sie das Modul v.centroids.

• Geben Sie als Name für die Ausgabe-Vektorkarte ‘forest_areas’ an und starten Sie das Modul.

• Nun laden Sie den neuen Layer forest_areas und visualisieren Sie die verschiedenen Waldtypen in verschiedenen Farben - deciduous, evergreen und mixed. Dazu öffnen Sie den Dialog Eigenschaften des Layers, wechseln zum Reiter Darstellung, wählen ‘Eindeutiger Wert’ und setzen das Klassifikationsfeld auf ‘VEGDESC’ (siehe auch Abschnitt sec_symbology).

• Als nächstes öffnen Sie wieder die GRASS-Werkzeugkiste und wechseln nach Vektor ‣ Vektor mit anderen Karten aktualisieren.

• Klicken Sie auf das v.rast.stats Modul. Geben Sie als Name der Rasterkarte für die Statistiken berechnet werden sollen gtopo30 und als Name der Vektorkarte (Polygone) forest_areas an.

• Nun fehlt nur noch ein Parameter: Geben Sie als Spaltenpräfix elev an, und klicken Sie dann auf [Starten]. Die Berechnung wird ziemlich lange dauern (bis zu mehreren Stunden).

• Wenn es erfolgreich beendet wurden. öffnen Sie den Layer forest_areas erneut und lassen Sie sich die Attributtabelle anzeigen. Dort gibt es nun weitere Spalten, in denen die Ergebnisse der univariate Statistik für Rasterwerte innerhalb der Vektorpolygone angezeigt werden, z.B.: elev_min, elev_max, elev_mean.

Der GRASS Datei-Browser¶

Eine weitere interessante Funktion in den GRASS-Werkzeugen ist der Browser. In figure_grass_module_7 können Sie die aktuelle LOCATION mit den dazu gehörigen MAPSETs sehen.

In der linken Hälfte des Datei-Browsers können Sie sich alle in der aktuellen Location enthaltenen Mapsets anschauen. In der rechten Hälfte des Fensters sehen Sie weitere Metainformationen zum ausgewählten Datensatz wie beispielsweise Auflösung, BoundingBox, Datenquelle, bei Vektordaten die Attributtabelle und die COMMAND Geschichte.

Figure GRASS module 7:

GRASS LOCATION browser

Die Werkzeugleiste im GRASS Datei-Browser stellt folgende Werkzeuge für die ausgewählte LOCATION zur Verfügung:

• Ausgewählten Layer dem Kartenfenster hinzufügen

• Gewählte Karte kopieren

• Gewählte Karte umbenennen

• Gewählte Karte löschen

• Setze die Region auf die gewählte Karte

• Neu zeichnen

Die Knöpfe Gewählte Karte umbenennen und Gewählte Karte löschen stehen nur für die aktuelle MAPSET zur Verfügung. Die anderen Werkzeuge funktionieren auch in den anderen MAPSETs.

<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE qgisgrassmodule SYSTEM "http://mrcc.com/qgisgrassmodule.dtd">

<qgisgrassmodule label="Vector buffer" module="v.buffer">
        <option key="input" typeoption="type" layeroption="layer" />
        <option key="buffer"/>
        <option key="output" />
</qgisgrassmodule>