Založení nové GRASS Location v systému S-JTSK¶
Založení Location ,,ručně``¶
http://grass.fsv.cvut.cz/wiki/S-JTSK:
$ grass64 -text
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.
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
- Jméno Location a mapsetu
- Vyberte projekci
- Další parametry
- Elipsoid
- Výsledek
- Hraniční souřadnice
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
