********************************************* 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 TO CONTINUE (OR TO CANCEL) Na další obrazovce potvrďte, že opravdu chcete založit novou Location:: LOCATION - 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 ? (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 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 > 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 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 TO CONTINUE (OR 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 TO CONTINUE (OR 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. .. image:: ../../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 ============================================ .. note:: Platí pro GRASS 6.4.x a novější Úvodní obrazovka .. image:: ../../images/location-wizard-01.png Jméno Location a mapsetu .. image:: ../../images/location-wizard-02.png Vyberte projekci .. image:: ../../images/location-wizard-03.png Další parametry .. image:: ../../images/location-wizard-04.png Elipsoid .. image:: ../../images/location-wizard-05.png Výsledek .. image:: ../../images/location-wizard-06.png Hraniční souřadnice .. image:: ../../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 :command:`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 :command:`gdalwarp` a :command:`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 .. image:: ../../images/srtm-import.png