[FreeGeoCZ] S-JTSK proj4 verzus OracleSpatial 10g

[Jano Kianicka <kianicka123jano AT yahoo.com>]

Ahojte,

Viem, ze hodne diskutovalo o chybnej definicii Krovakovej projekcie v kniznici proj4. Totiz testujem presnost noveho modulu v Oracle Spatial na transformaciu suradnic pomocou definicie cez EPSG kody. Ale, bohuzial Oracle je schopny transformovat len kod 2065 s nazvom "S-JTSK (Ferro) / Krovak". Ked vlozim testovacie suradnice z

http://gis.zcu.cz/studium/ugi/referaty/05/GeometrickeTransformace/index.html (uplne na spodku)

Tak mi z WGS84 long/lat (EPSG  8307):15.0092, 50.16535 vyrobí S-JTSK suradnice:

-559482,31; 1056474.27 .. pricom spravne suradnice S-JTSK su s presnstou 3m:
700000, 1040000 ako je uvedene v clanku.

Skusal som vyrobit podobnu transformaciu pre proj4:

echo "15d00'33.11  50d09'55.27" | cs2cs +proj=latlong +ellps=bessel +pm=ferro +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +to +init=epsg:2065
Ale vysledky su uplne ine. Pricom som skusal aj jednoduchsiu transformaciu WGS84 do UTM a tam sa vysledky Oracle, clanok sa zhoduju, teda plus minus 5m (EPSG 32633 zona 33).

Neviete mi poradit, co su vlastne 2065 suradnice v porovnani s klasickym S-JTSK (EPSG 4156)? Aky je ten posun, pripadne aka je definicia transformacie WGS84 -> do S-JTSK Ferro, ak hore uvedena je nespravna?

A tiez, ako presne definovat transformaciu WGS84 -> UTM zona33 pre proj4?

(Nie som s proj4 az taky skuseny ..)

Pripadne ak viete o nejakej studii co porovnava presnost transformacii klasickych GIS software (proj4, ESRI, Integraph, pripadne ten Oracle Spatial) ci sa neda niekde zohnat.

Potrebujeme totiz ist na vyzsiu presnost transformacii .. radovo az na cm.

Dakujem pekne za akukolvek radu.

Jano

---

Yahoo! oneSearch: Finally,  mobile search 
that gives answers, not web links.