Nächste Seite: Punktdaten durch Interpolation in
Aufwärts: Lösung mit GRASS
Vorherige Seite: Anzeigen der Punktdaten
  Inhalt
  Index
Um bei der anschließenden Interpolation zu verhindern, daß
über das gesammte Gebiet des Mapsets eine Interpolation
gerechnet wird, maskiert man das relevante Gebiet.
Das hierzu benutzte Kommando s.hull wurde für diese Aufgabe
nach einer Frage meinerseits an die GRASS-Mailingliste innerhalb von
3 Tagen aus dem vorhandenen Kommando s.delaunay abgeleitet.
Herzlichen Dank an Markus Neteler und Andrea Aime!
Das Kommando s.hull erzeugt einen konvexen Hüllenvektor
der alle Punkte umschließt.
GRASS:/spare/gisdata/grass/bin > s.hull
OPTION: name of a sites file to be input
key: sites
required: YES
Enter the name of an existing sites file
Enter 'list' for a list of existing sites files
Hit RETURN to cancel request
> tiefen
<tiefen>
OPTION: name of a vector file to be output
key: vect
required: YES
Enter a new binary file file name
Enter 'list' for a list of existing binary file files
Hit RETURN to cancel request
> mask
<mask>
FLAG: Set the following flag?
Use all sites (do not limit to current region)?(y/n) [n]
|
|
Diese Vektorlinie kann mit d.vect map=mask im
GRASS-Monitor angezeigt werden. Dieser von der Vektorlinie
umschlossenen Fläche muß ein Wert (Label) zugewiesen werden, um
aus ihr eine Rasterfläche zu erzeugen. Hierzu wird das Modul
v.digit aufgerufen, nachdem mit d.mon ein GRASS-Monitor
gestartet wurde. Da eine Vektordatei (die eben erzeugte Datei
mask) ja schon besteht und nicht neu erzeugt werden soll,
wird der Menüepunkt 3 none aufgerufen und anschließend
der Name der Vektordatei übergeben.
GRASS:/spare/gisdata/grass/bin > v.digit
Available Digitizers
Name Description
---- -----------
[1] Calcomp Calcomp digitizer, format 23 (binary)
[2] Altek Altek digitizer, model AC30, format 8 (binary)
[3] none Run digit without the digitizer.
Hit return to use digitizer in brackets
or type in number or name of other digitizer.
Select digitizer [none] : 3
Selected digitizer is: none
Enter the name of a map to work with.
If name is entered that does not already exist, it
will be created at this time.
DIGIT FILENAME
Enter 'list' for a list of existing digit files
Hit RETURN to cancel request
>mask
|
|
Im folgenden Menü können Angaben zu der Vektordatei gemacht
werden, wichtig ist hier die Zeile Maps's scale bei der
in userem Fall 1:2400 von GRASS eingetragen wurde, bei einer
Kartendigitalisierung muß natürlich der Maßstab der zugrundeliegenden
Karte übergeben werden. Die Ausdehnung der area die ja schon mit
s.info sites=tiefen ermittelt wurde, ist in die edge of
area Zeilen einzutragen.
Die Annahme der eingetragenen Informationen geschieht mit ESC
ENTER. Das anschließende shall we continue wird mit yes
bestätigt.
Provide the following information:
Your organization _____________________________
Todays date (mon,yr)Thu Dec 13 16:02:59
Your name ___________________
Map's name Convex hull map_________________________
Map's date __________
Map's scale 1:2400_____
Other info ________________________________________
Zone 0____
West edge of area 3464438.5_____
South edge of area 5901040_______
East edge of area 3464847_______
North edge of area 5901697_______
AFTER COMPLETING ALL ANSWERS, HIT <ESC><ENTER> TO CONTINUE
|
|
Der GRASS-Monitor zeigt jetzt die Vektorhülle mit zugehörigem Maßstab
und Nordpfeil. An der Kommandozeile wird das Bedienungsmenue für
GRASS-DIGIT ausgegeben. Mit dem Kommando L für Label
wechselt man in das entsprechende Untermenü.
?------------------------------------------------------------------------?
| GRASS-DIGIT Modified 4.10 Main menu |
|------------------------------------------------------------------------|
| MAP INFORMATION | AMOUNT DIGITIZED |
| Name: Convex hull map | # Lines: 0 |
| Scale: 2400 | # Area edges: 1 |
| Person: | # Sites: 0 |
| Dig. Thresh.: 0.0300 in. | - - - - - - - - - - - - - -|
| Map Thresh.: 1.8288 meters | Total points: 24 |
| | |
|------------------------------------------------------------------------|
| OPTIONS: | |
| | |
| Digitizer: Disabled | |
| | |
| | |
| | |
| |
?------------------------------------------------------------------------?
| Digitize Edit Label Customize Toolbox Window Help Zoom Quit * ! ^|
| |
| |
|GLOBAL MENU: Press first letter of desired command.[Upper Case Only] |
?------------------------------------------------------------------------?
|
|
In diesen Menue wird mit a die Option Label Areas
ausgewählt. Ein Name kann eingegeben werden, eine Category
muß angegeben werden. Die Category kann dazu dienen die Areas
zu unterscheiden, so könnten Beispielsweise Straßenflächen einheitlich
einer Category zugewiesen werden.
l------------------------------------------------------------k
| GRASS-DIGIT Modified 4.10 Label Menu |
|------------------------------------------------------------|
| Label options: |
| a - Label Areas m - Label Multiple Lines |
| l - Label Lines M - Un-Label Multiple Lines |
| s - Label Sites |
| c - Label Contours |
| A - Un-Label Areas i - Contour interval: < 5> |
| L - Un-Label Lines |
| S - Un-Label Sites |
| |
| B - Bulk Label Remaining Lines |
| |
| h - Highlight Lines of category # |
| d - Display Areas of category # |
| |
| q - Return to main menu |
m------------------------------------------------------------j
| Digitize Edit Customize Toolbox Window Help Zoom * ! ^ |
| |
| |
|GLOBAL MENU: Press first letter of desired command. |
|[Upper Case Only] |
m------------------------------------------------------------j
?------------------------------------------------------------?
| |
| |
| |
| Enter Category Number (0 to quit):[0] 1 |
?------------------------------------------------------------|
|
|
Anschließend wird im GRASS-Monitor erst ein Punkt innerhalb der
Fläche mit einem linken Mausklick markiert und
mit der mittleren Maustaste bestätigt, danach die umschließende
Vektorlinie. Die erzeugte area wird mit einem linken Mausklick
dem GRASS-Monitor zugewiesen (Abb. 3.2).
Nachdem man wieder im Label Menu angekommen ist, wird
es durch Eingabe von q und Q verlassen.
Abbildung:
GRASS-Monitor mit Digitalisierbildschirm v.digit
|
Select point within area:
Buttons:
Left: Choose this position
Middle: Accept chosen point position
Right: Abort/Quit
Point location:
EAST: 3464633.56
NORTH: 5901339.93
|
|
Um die Vektorlinie in eine Rasterdatei umzuwandeln,
verwendet man v.to.rast, nachdem
man mit v.support mask die von GRASS benötigten
weiteren Dateien erzeugt. Die Frage nach den im Speicher zu
haltenen Zeilen kann mit 512 bestätigt werden.
GRASS:~GRASS:/spare/gisdata/grass/bin > v.support mask
V.SUPPORT:
Selected information from vector header
Organization:
Map Name:
Source Date:
Orig. Scale: 0
No snapping will be done
Reading Vector file.
100%
Building areas
Building islands
Attaching labels
Number of lines: 1
Number of nodes: 1
Number of areas: 1
Number of isles: 1
Number of atts : 0
Number of unattached atts : 0
Snapped lines : 0
|
|
GRASS:/spare/gisdata/grass/bin > v.to.rast
OPTION: vector input file
key: input
required: YES
Enter the name of an existing vector file
Enter 'list' for a list of existing vector files
Hit RETURN to cancel request
> mask
Enter a new raster file name
Enter 'list' for a list of existing raster files
Hit RETURN to cancel request
> mask
Parameters:
input vector input file
output raster output file
rows number of rows to hold in memory
default: 512
required: NO
enter option >
Loading vector information ... 0 mins 00 secs
Sorting areas by size ... 1 areas 0 mins 00 secs
Pass #1 (of 2)
Processing areas ... 1 areas 0 mins 00 secs
Processing lines ... 0 lines 0 mins 00 secs
Processing sites ... 0 sites 0 mins 00 secs
Writing raster map ... 0 mins 00 secs
Pass #2 (of 2)
Processing areas ... 1 areas 0 mins 00 secs
Writing raster map ... 0 mins 00 secs
Creating support files for raster map ... 0 mins 01 sec
Raster map <mask> done.
Total processing time: 0 mins 01 sec
|
|
Die Maske selbst wird in dem Befehl r.mask gesetzt,
mit dem alternativ auch Masken entfernt werden können.
Die erfolgte Maskierung erzeugt eine Rasterdatei MASK
die mit d.rast -o map=MASK angezeigt werden kann.
GRASS:/spare/gisdata/grass/bin > r.mask mask
MASK: Program for managing current GIS mask
current mask: none
Options:
1 Remove the current mask
2 Identify a new mask
RETURN Exit program
>2
Enter name of data layer to be used for mask
Enter 'list' for a list of existing raster files
Enter 'list -f' for a list with titles
Hit RETURN to cancel request
>mask
IDENTIFY THOSE CATEGORIES TO BE INCLUDED IN THE MASK
OLD CATEGORY NAME CAT
NUM
. . . . . . . . . . . . . . . . . . . . . . . . . . 0 0_
. . . . . . . . . . . . . . . . . . . . . . . . . . 1 1_
Next category: end__ (0 thru 1)
AFTER COMPLETING ALL ANSWERS, HIT <ESC><ENTER> TO CONTINUE
(OR <Ctrl-C> TO CANCEL)
|
|
Nächste Seite: Punktdaten durch Interpolation in
Aufwärts: Lösung mit GRASS
Vorherige Seite: Anzeigen der Punktdaten
  Inhalt
  Index
Heiko Kehlenbrink
2002-01-15