Převod vrstvy shapefile z WGS84 na JTSK
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
3) zvolíme volbu RESAMPLE - SGAPEFILE - UTM to LAT/LONG
4) zvolíme jaké SHP budeme převádět
5) a jak se bude výsledná SHP jmenovat
6) je třeba nastavit souřadnicový systém - pro většinu ČR je UTM = 33
7) a zobrazení ( po stisknutí MERCATOR ) zvolit jako WGS84
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.
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
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.
4) přepneme vlevo dole v TOC na LAYERS
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.
Komentáře
Přehled komentářů
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.
Díky
(Martin, 20. 12. 2007 21:57)