1. Ermitteln der Wertebereiche aus den Datensätzen:

r.describe -r map=A

Ausgabe: 0 - 273083.302800 (Beispiel)

r.describe -r map=A2

Ausgabe: 0 - 286798.596480 (Beispiel)

Beide Karten sollen auf der Grundlage des grösseren Wertebereiches in gleiche Kategorien eingeteilt werden: Wenn jede Kategorie 10000 umfassen soll, brauchen wir 28 Zeilen Code für den Befehl r.reclass:

10000 thru 20000 = 1
20000 thru 30000 = 2
30000 thru 40000 = 3
...
260000 thru 270000 = 26
270000 thru 280000 = 27
280000 thru 290000 = 28

Die Anweisungen für r.reclass bereiten eine Menge Tipparbeit, die man sich ersparen kann.

Sehen Sie sich hierzu das awk-Skript an. Speichern Sie das Skript in ihrem Kursverzeichnis. Ändern Sie nach dem Speichern die Berechtigung mit dem Kommando

chmod +x zerschneiden

#!/usr/bin/awk -f # Bedeutung der ersten Zeile: # Diese Datei wird von awk als Programm interpretiert BEGIN { zmin = 0 zmax = 290000 inc = 10000 steps = (zmax - zmin) / inc for (i = 0; i < steps; i++) { printf "%d thru %d = %d\n", zmin+i*inc, zmin+(i+1)*inc, i+1 } }

A User's Guide for GNU Awk

Mit dem Skript zerschneiden kann das Reklassifizieren dann so gestartet werden:

./zerschneiden | r.reclass i=A o=A_cat

./zerschneiden | r.reclass i=A2 o=A2_cat

Aufgabe 1:
Passen Sie das Skript auf die Wertebereiche ihrer im letzten Kurs erstellten Ergebniskarten an.
Führen Sie das Skript dann für alle Karten (mindestens 3 verschiedene) aus. Interpretieren Sie die Ergebnisse.
d.rast .....

Aufgabe 2:
Berechnen Sie die Veränderung des Bodenabtrags bei einer hypothetischen Bodenschutzmassnahme. Die Massnahme soll mit dem Faktor 0.3 in die Berechnung eingehen - allerdings nur in den nicht bewaldeten Gebieten (ausgehend von der Ursprungskarte vegcover).
Speichern Sie dazu ihr Skript unter einem neuen Namen und passen Sie es entsprechend an. Beschreiben Sie ihr Vorgehen und die Ergebnisse.