next up previous contents index
Nächste Seite: Punktdaten durch Interpolation in 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.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
\includegraphics[width=0.8\textwidth]{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 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)


next up previous contents index
Nächste Seite: Punktdaten durch Interpolation in Aufwärts: Lösung mit GRASS Vorherige Seite: Anzeigen der Punktdaten   Inhalt   Index
Heiko Kehlenbrink 2002-01-30