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

    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