Pénteken töltöttünk le az OKTVF szerveréről Európai Uniós Szlovák vízügyi állományokat, amelyeket HD72/EOV rendszerbe kellett volna transzformálnom, hogy össze lehessen kapcsolni a mi vízrajzi állományainkkal. Két scriptet is írtam pythonban a shape fájlok reprojektálásra.
arcgisscripting 9.3 ill gdal/ogr , azonban egyik sem oldotta meg a feladatot. Arcgis egyedül a GCS_ETRF_1989 dátummal ellátott shapeket alakította át helyesen, a S-JTSK_KROVAK_EAST_NORTH a GCS_ETRS_1989 és a WGS84 (!!!) is hibás volt. Az open source szoftvereknél az EPSG:23700 és az EPSG:2065 is hibás volt, ahhoz, hogy ugyanúgy váltson át, mint az Esri szoftver, ezt a transzformációt kellett használnom:
(Az ötletet innen vettem)
ogr2ogr -s_srs '+proj=krovak +lat_0=49.5 +lon_0=42.5 +alpha=30.28813972222222 +k=0.9999 +x_0=-0 +y_0=-0 +ellps=bessel +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +pm=ferro +to_meter=-1 +no_defs' -t_srs '+proj=somerc +lat_0=47.14439372222 +lon_0=19.048571778 +k_0=0.99993 +x_0=650000 +y_0=200000 +ellps=GRS67 +units=m +towgs84=52.684,-71.194,-13.975,-0.312,-0.1063,-0.3729,1.0191 +no_defs' /home/szabobal/proba/rivers.shp /home/szabobal/Asztal/Szomszedos_orszagok/SK/rivers.shp
Az alábbi térképen egyértelműen látható, hogy a világoskék (Arcgis transzformáció) és a fekete (fenti ogr2ogr) ugyanúgy vált át, ami azonban kb 50 méterrel el van csúszva a helyestől... Nem tudom, hogy mi lehet a megoldás, én arra is tippeltem, hogy direkt rosszak az adatok. Aki tudja a megoldást ne habozzon:)
Nincsenek megjegyzések:
Megjegyzés küldése