Eine pragmatische Lösung für diese Aufgabe ist es, alle Punkte einer
Rasterdatei, deren (Höhen) - Wert unter einer
bestimmten Vorgabe liegen auszublenden, also durch eine Maske zu
verdecken. Hierzu bedient man sich des schon in
Kapitel 4.4 ab der Seite
vorgestellten
Programmes r.mapcalc, das direkte arithmetische und logische
Operationen auf ein Raster zulässt. In diesem Fall genügt ein r.mapcalc Ergebnissdatei=``if(Ausgangsdatei <= 1900, null(),
Ausgangsdatei)``. Hierdurch wird eine neue Ergebnissdatei erzeugt,
in der alle Werte der Ausgangsdatei, die kleiner oder gleich 1900
sind, auf Null(), also ,,nicht existent`` gesetzt werden. Alle anderen
Werte (der,, Else-Fall``) werden beibehalten. Eine Anzeige mit d.rast Ergebnissdatei zeigt das Ergebniss im aktiven
Monitor an (Abb. 6.1)
|
Soll die Ergebnissdatei als Maske gebraucht werden so kann man mit
dem Aufruf r.mapcalc Ergebnissdatei=``if(Ausgangsdatei
1900, null(), 200)``
dem nicht gefluteten Bereichen ein beliebiger Wert (hier 200) zugewiesen oder
durch r.mapcalc 'Ergebnissdatei=if(Ausgangsdatei
1900, 200, null())'
das inverse Ergebniss erzeugt werden
(alle Bereiche über 1900 werden ausgeblendet). Der Zahlenwert der Maske ist
nicht entscheident, da durch GRASS alle Gebiete maskiert werden,
die nicht null() sind.
|