Jdi na obsah Jdi na menu
 


Převod vrstvy shapefile z WGS84 na JTSK

12. 2. 2006
Předem se omlouvám všem pokročilým a zkušeným uživatelům GIS aplikací za svoji krajně neodbornou terminologii. Hledal jsem nejlehčí ( funkční ) cestu pro převod vrstvy shapefile z WGS na JTSK. Našel jsem celou řadu aplikací, ale žádná nefungovala na 100ˇ%. Takže jsem si dovolil zvolit takovýto postup. Třeba bude někomu k užitku.

Pro převod potřebujeme dva freeware progarmy. Jsou free pouze pro nekomerční použití !

1 ) MicroDem ( http://www.usna.edu/Users/oceano/pguth/website/microdem.htm )
2) Janitor ( http://janitor.ten.cz/j2/www/j2_intro.php?lang=cze&idmn=31 )

Vlastní postup tedy je :

1) otevřeme MicroDem

Obrazek

2) stiskneme buton IN - OUT a objeví se utilita DATA MANIPULATION

Obrazek


3) zvolíme volbu RESAMPLE - SGAPEFILE - UTM to LAT/LONG

Obrazek


4) zvolíme jaké SHP budeme převádět

Obrazek


5) a jak se bude výsledná SHP jmenovat

Obrazek


6) je třeba nastavit souřadnicový systém - pro většinu ČR je UTM = 33

Obrazek


7) a zobrazení ( po stisknutí MERCATOR ) zvolit jako WGS84

Obrazek


8) následuje vlastní převod. Nevím proč, ale následně se objeví tato obrazovka - ale převod je stejně proveden a je možno MicroDem zavřít. Nová SHP je ve zvoleném adresáři.

Obrazek

V případě, že převod proběhne zcela bez závad, objeví se místo červeného pruhu pole : Done process : jméno SHP.

Obrazek


Takto jsme převedli SHP ve WGS84 z rovinného zápisu do zeměpisného. Nyní je třeba ještě ji převést na JTSK. Pro to použijeme Janitor.

1) otevřeme Janitor

Obrazek

2) vybereme SHP pro převod

Obrazek



3) kliknutím pravým tlačítkem myši na vrstvu nastavíme souřadnicový systém vrstvy na WGS84. Dále je nutné nastavit souřadnicový systém celého zobracení ( SETTINGS - VIEW ) na JTSK.

Obrazek



4) přepneme vlevo dole v TOC na LAYERS

Obrazek



5) pravým tlačítkem myši klikneme na vrstvu a nastvavíme EXPORT ( zaškrtneme pole - převést podle zobrazení - VIEW ). To je z důvodu, že naše vrstva je ve WGS84 a celé zobrazení v JTSK. Převod i výsledem tedy bude v JTSK. Nová SHP je ve stejném adresář jako původní převáděná vrstva.

Obrazek






 

Komentáře

Přidat komentář
Kontrolní kód Přehrát audio kód

Přehled komentářů

Díky

(Martin, 20. 12. 2007 21:57)

Díky za tip, podívám se časem na to. Nejsem v GISu vzdělán a hledal jsem něco relativně snadného, funkčního a free. Cokoliv jednoduššího jen potěší.

proc tak slozite? jde to i bez klikani

(zemze, 20. 12. 2007 21:47)

Mame programy jako cs2cs a ogr2ogr

JTSK->WGS84
-----------
echo "-734542.9904476 -1045586.42733531" |cs2cs +proj=krovak +ellps=bessel +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +to +init=epsg:4326

ogr2ogr -s_srs krovak.prf -t_srs '+proj=latlong +datum=WGS84' -f "ESRI Shapefile" vektor-out vektor.shp


WGS84->JTSK
-----------
echo "12d48'25.15992 49d27'8.14571" | cs2cs +proj=latlong +datum=WGS84 +to +init=esri:102065 +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56

ogr2ogr -s_srs krovak.prf -t_srs '+proj=latlong +datum=WGS84' -f "ESRI Shapefile" vektor-out sil_cz12x.shp


krovak.prf:

PROJCS["Krovak",
GEOGCS["Bessel 1841",
DATUM["D_unknown",
SPHEROID["bessel",6377397.155,299.1528128]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Krovak"],
PARAMETER["latitude_of_center",49.5],
PARAMETER["longitude_of_center",24.83333333333333],
PARAMETER["scale_factor",0.9999],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["Meter",-1]]

Arcdata

(MArtin, 14. 3. 2006 7:33)

Jasný, koukal jsem. Ty SHP by se převést daly, ale Arcdatům raději do jejich práv nepolezu. Pro free data ve wgs se podívejte na ( http://laker.wz.cz/dp/ggi_dp.html ) a ještě raději přímo ( http://grass.fsv.cvut.cz/wiki/index.php/Geografick%C3%A1_data_%C4%8CR_pro_GRASS ). Tam jsou data ve WGS84 pro celou Čr.

Arcdata

(Pepa, 14. 3. 2006 6:44)

Arcdata mají na serveru ukázková a výuková data dvou krajů v shp. http://www.arcdata.cz/support/download#data

Arcdata

(to Ladislav, 13. 3. 2006 21:26)

Zdravím, leč nepotěším. Nevím přesně o jaká data je zájem. Předpokláádm o administrativní členění vycházející z ARC500. To by byl problém. Za 1) patří autorská práva Arcdatům a za druhé upravení a převedení dat není v fGIS až tak jednoduché, neboť původní formát je určen pro Arc GIS a bez něj se nikam nedostanete. To by Arcdata nusela data prezentovat jako SHP. Což moudře nečiní. Takže demodata ( původní exe soubor ) není možné bez Arc GIS převést. Pro podobná data doporučuji server ESRI. Při troše trpělivosti lze získat okresy i původní ( ne současné !! ) kraje. MArtin

to Ladislav

(Pepa, 13. 3. 2006 14:54)

Chtěl jsem něco zkusit s demodaty co mají arcdata na stránkách a potřeboval bych je ve wgs84. Kdybyste měl čas a chuť a mohl mi to převést a poslat, byl bych zauzlován.
zxc@post.cz

to Ladislav

(MArtin, 9. 3. 2006 22:18)

Naše SHP vznikly na základě souřadnic z komerční GPS ( viz. jak vytváříme SHP v fGIS ). Už to samo o sobě vytváří chybu v řádu metrů. Navíc, převod na JTSK jsem zařadil jen pro možnost využití dat pro další aplikace. Zde tedy pro server Cenia. Přesnější zaměřování památek je hudba budoucnosti ( cena stanic je mimo naše nekomerční možnosti ). A i vzhledem k tomuto nepoužíváme komerční sw. Hlavní tíhu pro zpracování dat nese právě fGIS. Velké díky za nabídku, ale nevím jak bychom ji mohli využít. I když do budoucna neodmítáme nic. MArtin

Převod

(Ladislav Fofonka, 9. 3. 2006 21:51)

Jsem autorem programu pro obousměrný převod WGS84 <-> S-JTSK pro geodetické účely (s centimetrovou přesností). Bohužel neznám vámi používane shapes, pracuji s AutoCADem. Mohu nějak pomoci?
fofonka@tiscali.cz

JTSK -> WGS84

(Pepa, 2. 3. 2006 19:22)

zkusím se s tím poprat. Bohužel pracuji s hotovými daty, která jsou v jtsk a která bych potřeboval ve wgs84 a navíc mám v téhle problematice pořádný guláš. tak uvidim, jak to dopadne..

pro Pepa

(MArtin, 2. 3. 2006 14:06)

To jsem nezkoušel, neb ještě nebyla taková potřeba. SHP dělám rovnou ve WGS pomocí FGIS. Ale Jan Map umí převádět mezi WGS84 - JTSK a S-42. A Microdem ve zmíněných utilitách má jak funkci UTM to LAT/LONG, tak i LAT/LONG to UTM. Takže čistě teoreticky je postup stejný. Ale jak říkám, ještě jsem to nezkoušel.

JTSK -> WGS84

(Pepa, 2. 3. 2006 11:34)

Zajímal by mě postup pro opačný převod