Založení nové GRASS Location v systému S-JTSK

Založení Location ,,ručně``

nebo

$ grass64 -gui

Vyplňte formulář:

.
                        GRASS 6.4.cvs
DATABASE: A directory (folder) on disk to contain all GRASS maps and data.

LOCATION: This is the name of a geographic location. It is defined by a
        co-ordinate system and a rectangular boundary.

MAPSET:   Each GRASS session runs under a particular MAPSET. This consists of
        a rectangular REGION and a set of maps. Every LOCATION contains at
        least a MAPSET called PERMANENT, which is readable by all sessions.

        The REGION defaults to the entire area of the chosen LOCATION.
        You may change it later with the command: g.region
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

LOCATION:   cr   (enter list for a list of locations)
MAPSET:     jachym  (or mapsets within a location)

DATABASE: /home/jachym/grassdata



        AFTER COMPLETING ALL ANSWERS, HIT <ESC><ENTER> TO CONTINUE
                            (OR <Ctrl-C> TO CANCEL)

Na další obrazovce potvrďte, že opravdu chcete založit novou Location:

LOCATION <cr> - doesn't exist

Available locations:
----------------------
LausanneDemoData                        LausanneDemoData.tar.bz2
cr-grass-jtsk0.3.0.tar.gz              cr-jtsk
cr-wgs84                                firedata
hazemburk                               imagery60
krkonose                                krtiny
latlong                                 ncspm02
ncspm0220070212.tar.gz               slovakia3d
spearfish60                             xy
----------------------

Would you like to create location <cr> ? (y/n) [y]

Potvrďte, že máte všechny potřebné informace pro založení nové Location:

To create a new LOCATION, you will need the following information:

1. The coordinate system for the database
        x,y (for imagery and other unreferenced data)
        Zeměpisná šířka a délka
        UTM
        Ostatní mapové projekce
2. The zone for the UTM database
and all the necessary parameters for projections other than
Zeměpisná šířka a délka, x,y, and UTM
3. The coordinates of the area to become the default region
and the grid resolution of this region
4. A short, one-line description or title for the location

Do you have all this information? (y/n) [y]

Zadejte, že chcete vytvořit Jinou mapovou projekci:

Please specify the coordinate system for location <cr>

A   x,y
B   Zeměpisná šířka a délka
C   UTM
D   Ostatní mapové projekce
RETURN to cancel

> D

Ostatní mapové projekce coordinate system? (y/n) [y]

Zadejte, krátký popis nové Location:

Please enter a one line description for location <cr>

> Location pro celou ceskou republiku
=====================================================
Location pro celou ceskou republiku
=====================================================
ok? (y/n) [y]

Zadejte mapovou projekci krovak. Pro Výpis všech dostupných projekcí zadejte list:

Specifikujte prosím mapovou projekci
Zadejte 'list' pro výpis dostupných mapových projekcí
Pro ukončení požadavku stiskněte RETURN
>
list

ll -- Lat/Lon
utm -- Universe Transverse Mercator
stp -- State Plane
aea -- Albers Equal Area
lcc -- Lambert Conformal Conic
merc -- Mercator
tmerc -- Transverse Mercator
leac -- Lambert Equal Area Conic
laea -- Lambert Azimuthal Equal Area
aeqd -- Azimuthal Equidistant
airy -- Airy
...

krovak

Ano, chcete zadat geodetické datum pro vaši Location:

Do you wish to specify a geodetic datum for this location?(y/n)[y] y

A nyní zadejte datum:

Specifikujte prosím výškový systém (datum)
Zadejte 'list' pro výpis dostupných výškových systémů (datumů)
or 'custom' if you wish to enter custom parameters
Pro ukončení požadavku stiskněte RETURN
>
hermannskogel

Vyberte transformační parametry:

Now select Datum Transformation Parameters
Please think carefully about the area covered by your data
and the accuracy you require before making your selection.

Enter 'list' to see the list of available Parameter sets
Enter the corresponding number, or <RETURN> to cancel request
>
list

---
3       Used in Czech Republic
        (PROJ.4 Params towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56)

---
4       Used in Slovakia
        (PROJ.4 Params towgs84=485.021,169.465,483.839,7.786342,4.397554,4.10265
5,0)

---

3

Množné číslo jednotek jsou metry, můžete ponechat:

Enter plural form of units [meters]:

Zadejte hraniční souřadnice všech nových mapsetů:

            DEFINE THE DEFAULT REGION

           ====== DEFAULT REGION =======
           | NORTH EDGE:-934410.20     |
           |                           |
WEST EDGE  |                           |EAST EDGE
-905013.19 |                           |-430661.31
           | SOUTH EDGE:-1227585.7     |
           =============================

PROJECTION: 99 (Ostatní mapové projekce) ZONE: 0

                    GRID RESOLUTION
                      East-West:     30
                    North-South:     30

AFTER COMPLETING ALL ANSWERS, HIT <ESC><ENTER> TO CONTINUE
                (OR <Ctrl-C> TO CANCEL)

Potvrďte souhrnné informace:

projection: 99 (Ostatní mapové projekce)
zone: 0
north:       -934410.2
south:       -1227585.7
east:        -430661.31
west:        -905013.19

e-w res:     29.99948647  (Changed to conform to grid)
n-s res:     29.99851632  (Changed to conform to grid)

total rows: 9773
total cols: 15812
total cells: 154,530,676


Do you accept this region? (y/n) [n] > y

Location vytvořena, pokračujte opět v úvodní obrazovce, založením Mapsetu:

.
                        GRASS 6.4.cvs
DATABASE: A directory (folder) on disk to contain all GRASS maps and data.

LOCATION: This is the name of a geographic location. It is defined by a
        co-ordinate system and a rectangular boundary.

MAPSET:   Each GRASS session runs under a particular MAPSET. This consists of
        a rectangular REGION and a set of maps. Every LOCATION contains at
        least a MAPSET called PERMANENT, which is readable by all sessions.

        The REGION defaults to the entire area of the chosen LOCATION.
        You may change it later with the command: g.region
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

LOCATION:   cr   (enter list for a list of locations)
MAPSET:     jachym  (or mapsets within a location)

DATABASE: /home/jachym/grassdata



        AFTER COMPLETING ALL ANSWERS, HIT <ESC><ENTER> TO CONTINUE
                            (OR <Ctrl-C> TO CANCEL)

Víjtete v GRASSu! Ověřte nastavení projekce:

GRASS> g.region -p

projection: 99 (Krovak)
zone:       0
datum:      hermannskogel
ellipsoid:  bessel
north:      -934410.2
south:      -1227585.7
west:       -905013.19
east:       -430661.39
nsres:      29.99851632
ewres:      29.99948141
rows:       9773
cols:       15812
cells:      154530676

GRASS> g.proj -p

-PROJINFO-------------------------------------------------
name       : Krovak
datum      : hermannskogel
towgs84    : 570.8,85.7,462.8,4.998,1.587,5.261,3.56
proj       : krovak
ellps      : bessel
a          : 6377397.1550000003
es         : 0.0066743722
f          : 299.1528128000
-PROJUNITS------------------------------------------------
unit       : meter
units      : meters
meters     : 1.0

Založení Location pomocí kódu EPSG

Pokud zakládáte Location v sysétmu S-JTSK, vložte do souboru s EPSG kódy vlastní záznam proj S-JTSK

Pokud tak neučiníte, zvolte 2065 a přidejte transformační parametry.

../_images/krovak-epsg-gui.png

Ověřte projekci nové Location:

GRASS> g.proj -p

-PROJINFO-------------------------------------------------
name       : Krovak
proj       : krovak
datum      : hermannskogel
ellps      : bessel
lat0       : 49.5
lon0       : 42.5
alpha      : 30.28813972222222
k          : 0.9999
x0         : 0
y0         : 0
pm         : ferro
nodefs     : defined
towgs84    : 570.8,85.7,462.8,4.998,1.587,5.261,3.56
-PROJUNITS------------------------------------------------
unit       : metre
units      : metres
meters     : 1

GRASS> g.region -p

projection: 99 (Krovak)
zone:       0
datum:      hermannskogel
ellipsoid:  bessel
north:      1
south:      0
west:       0
east:       1
nsres:      1
ewres:      1
rows:       1
cols:       1
cells:      1

Založení Location pomocí grafického průvodce

Poznámka

Platí pro GRASS 6.4.x a novější

Úvodní obrazovka
../_images/location-wizard-01.png
Jméno Location a mapsetu
../_images/location-wizard-02.png
Vyberte projekci
../_images/location-wizard-03.png
Další parametry
../_images/location-wizard-04.png
Elipsoid
../_images/location-wizard-05.png
Výsledek
../_images/location-wizard-06.png
Hraniční souřadnice
../_images/location-wizard-07.png

Založení Location z existujících dat

  • Pro import rastrových dat používáme některý z modulů začínajících na r.in.*

  • Modul :command`r.in.gdal` je rozhraním ke knihovně GDAL – http://gdal.org

  • Knihovna GDAL umí číst a zapisovat z a do celé řady rastrových formátů.

  • Seznam podporovaných formátů ve vašem systému dostanete příkazem:

    $ gdalinfo --formats
  • Součástí distribuce GDAL je i knihovna OGR, která umí totéž s vektorovými daty:

    $ ogrinfo --formats
  • Programy gdalwarp a ogr2ogr slouží nejen k převodům mezi různými formáty a především i k převodům mezi souřadnými systémy.

Zdroj dat: SRTM http://srtm.csi.cgiar.org

$ gdalinfo data/Z392.TIF

Driver: GTiff/GeoTIFF
Size is 6000, 6000
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS1984",
        SPHEROID["WGS 84",6378137,298.2572235629972,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]

Origin = (10.000000000000000,55.000000000000000)
Pixel Size = (0.000833333333333,-0.000833333333333)
Metadata:
AREAORPOINT=Area
Corner Coordinates:
Upper Left  (  10.0000000,  55.0000000) ( 10d 0'0.00"E, 55d 0'0.00"N)
Lower Left  (  10.0000000,  50.0000000) ( 10d 0'0.00"E, 50d 0'0.00"N)
Upper Right (  15.0000000,  55.0000000) ( 15d 0'0.00"E, 55d 0'0.00"N)
Lower Right (  15.0000000,  50.0000000) ( 15d 0'0.00"E, 50d 0'0.00"N)
Center      (  12.5000000,  52.5000000) ( 12d30'0.00"E, 52d30'0.00"N)
Band 1 Block=6000x1 Type=Int16, ColorInterp=Gray NoData Value=-32768

Převod do S-JTSK (epsg:102067):

$ gdalwarp -t_srs +init=epsg:102067  Z392.TIF Z392-JTSK1.tif

Vypsání výsledku:

$ gdalinfo  Z392-JTSK.tif

Driver: GTiff/GeoTIFF
Size is 5533, 7862
Coordinate System is:
PROJCS["S-JTSK (Ferro) / Krovak",
    GEOGCS["S-JTSK (Ferro)",
        DATUM["SJTSKFerro",
            SPHEROID["Bessel 1841",6377397.155,299.1528128000033,
                AUTHORITY["EPSG","7004"]],
            AUTHORITY["EPSG","6818"]],
        PRIMEM["Ferro",-17.66666666666667],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4818"]],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","2065"]]}

Origin = (-1056755.290744186844677,-454594.250338208337780)
Pixel Size = (76.787259777712435,-76.787259777712435)
Metadata:
AREAORPOINT=Area
Corner Coordinates:
Upper Left  (-1056755.291, -454594.250)
Lower Left  (-1056755.291,-1058295.687)
Upper Right ( -631891.382, -454594.250)
Lower Right ( -631891.382,-1058295.687)
Center      ( -844323.337, -756444.969)
Band 1 Block=5533x1 Type=Int16, ColorInterp=Gray

Převod pomocí parametrů z aktuální GRASS Location:

GRASS> g.proj -wef

GRASS> gdalwarp -t_srs "`g.proj -wef`" Z392.TIF Z392-JTSK2.tif
GRASS> gdalinfo  Z392-JTSK2.tif

Imporr do GRASSu:

GRASS> r.in.gdal in=Z392-JTSK3.tif out=srtm-z-39-2
GRASS> r.in.gdal in=Z392-JTSK3.tif out=srtm-z-39-2 -o
GRASS> d.mon x0
GRASS> d.rast srtm-z-39-2

Informace o importovaném rastru:

GRASS> r.info srtm-z-39-2
...
Range of data:    min = -32768  max = 1244
...
GRASS> r.null srtm-z-39-2 setnull=0,-32768
GRASS> r.colors srtm-z-39-2 rules=elevation
GRASS> d.redraw
../_images/srtm-import.png