next up previous contents index
Nächste Seite: Punktdaten in Rasterdaten umwandeln Aufwärts: Lösung mit GRASS Vorherige Seite: Anzeigen der Punktdaten   Inhalt   Index

Maskierung des befahrenen Hafengebietes

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.delauney 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 1:1 eingetragen wird und die Ausdehnung der area die mit s.info sites=tiefen ermittelt werden kann. 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)___________________
Your name           ___________________
Map's name          ___________________
Map's date          __________
Map's scale         1:1________
Other info          ___________________
Zone                0____
West edge of area   3464438_______
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ü.

l----------------------------------------------------------------k
|GRASS-DIGIT Modified  4.10                           Main menu  |
|----------------------------------------------------------------|
|MAP INFORMATION                    | AMOUNT DIGITIZED           |
|Name:                              |   # Lines:         0       |
|Scale:        1                    |   # Area edges:    1       |
|Person:                            |   # Sites:         0       |
|Dig. Thresh.: 0.0300 in.           | - - - - - - - - - - - - - -|
|Map Thresh.:  0.0008 meters        |   Total points:    24      |
|                                   |                            |
|----------------------------------------------------------------|
|OPTIONS:                           |                            |
|                                   |                            |
|     Digitizer:     Disabled       |                            |
|                                   |                            |
|                                   |                            |
|                                   |                            |
|                                                                |
m----------------------------------------------------------------j
|Digitize Edit Label Customize Toolbox Window Help Zoom Quit * !^|
|                                                                |
|                                                                |
|GLOBAL MENU:Press first letter of desired command.              |
|[Upper Case Only]                                               |
m----------------------------------------------------------------j

In diesen Menue wird mit a die Option Label Areas ausgewählt. Ein Name kann eingegeben werden, eine Category muß angegeben 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 (s. 4.2). Nachdem man wieder im Label Menu angekommen ist wird v.digit durch Eingabe eines kleinen und eines großen q verlassen.

Abbildung: GRASS-Monitor mit Digitalisierbildschirm v.digit
\includegraphics[width=6cm]{dem-from-sitedata/echolot-hull-label.eps}

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 gleichzeitig auch Masken entfernt werden können.

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)


next up previous contents index
Nächste Seite: Punktdaten in Rasterdaten umwandeln Aufwärts: Lösung mit GRASS Vorherige Seite: Anzeigen der Punktdaten   Inhalt   Index
Heiko Kehlenbrink 2001-11-27