Dotaz: Převod GPS souřadnic pomocí cs2cs popř. jiného sw

[mirek19]

Zdravím všechny. Marně hledám jednoduchý sw na převod GPS souřadnic WGS-84 <-> S-42. Vygooglil jsem, že je na to vhodný cs2cs (ve Fedoře balík proj), jenže z toho kvanta blíže nepopsaných parametrů (patrně umožňujících zadávat libovolný obecný souřadnicový systém) nechápu, jak provést zrovna tyhle převody mezi WGS-84 a S-42. Nevíte někdo, jak parametry nastavit? Anebo nevíte o jiným vhodným sw? Díky moc.

Odpovědi

[Xóža Hňabaj]

Na převody souřadných systémů je nejlepší knihovna proj4.

http://proj.maptools.org/gen_parms.html

[Tomáš Zigo]

Zdravím,

Transformácia súradníc zo súradnicového systému WGS-84 -> S-42 pomocou cs2cs na základe EPSG kódov definovaných v EPSG databáze:

echo "12d48'25.15992 49d27'8.14571 559.417" | cs2cs +init=epsg:4326 +to +init=epsg:28403 \
+towgs84=26,-121,-79,0.025,0.035,-0.025,1

3341121.07	5482379.62 504.48

EPSG kód 4326 definuje súradnicový systém WGS-84 z ktorého súradnice transformujeme, pričom EPSG kód 28403 definuje súradnicový systém S-42 do ktorého súradnice transformujeme. Parameter +towgs84=26,-121,-79,0.025,0.035,-0.025,1 je spresňujúci (o málo viac presnejší ako prednastavený transformačný kľúč definovaný v EPSG kóde 28403) transformačný klúč pre transformáciu do súradnicového systému S-42.

Ekvivalentný zápis s celou definíciou oboch súradnicových systémov WGS-84 (+to) -> S-42:

echo "12d48'25.15992 49d27'8.14571 559.417" | cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 \ 
+no_defs +to +proj=tmerc +lat_0=0 +lon_0=15 +k=1 +x_0=3500000 +y_0=0 +ellps=krass \
+units=m +no_defs +towgs84=26,-121,-79,0.025,0.035,-0.025,1

3341121.07	5482379.62 504.48

Opačná transformácia zo súradnicového systému S-42 -> WGS-84:

echo "3341121.07 5482379.62 504.48" | cs2cs +init=epsg:28403 +towgs84=26,-121,-79,0.025,0.035,-0.025,1 \
+to +init=epsg:4326

12d48'25.16"E	49d27'8.146"N 559.417

Ekvivalentný zápis:

echo "3341121.07 5482379.62 504.48" | cs2cs +proj=tmerc +lat_0=0 +lon_0=15 +k=1 +x_0=3500000 +y_0=0 +ellps=krass \
+units=m +no_defs +towgs84=26,-121,-79,0.025,0.035,-0.025,1 +to +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs 

12d48'25.16"E	49d27'8.146"N 559.417

Na tej to úrovni som tomu porozumel ja, snáď ti to teda pomôže. Výchadzal som z GRASS GIS Wiki CZ.