Ligera traducción con transformación de coordenadas de Gauss-Kruger a WGS 84. ¿Posibles razones?

Question

Tengo una serie de puntos en Gauss-Kruger, Zona 3, que quiero convertir para visualizarlos dentro de Google-Earth. Asumiendo WGS84 para Google-Earth las coordenadas transformadas difieren de las originales unos 140m al noreste.

Estoy haciendo la transformación de coordenadas utilizando la "Biblioteca de proyecciones cartográficas PROJ.4" ( http://trac.osgeo.org/proj/ ) con las siguientes cadenas de proj4:

Gauß Krüger, Zona 3:

+proj=tmerc +lat\_0=0 +lon\_0=9 +k=1.000000 +x\_0=3500000 +y\_0=0 +ellps=bessel +datum=potsdam +units=m +no\_defs **WGS84:**     +proj=longlat +ellps=WGS84 +datum=WGS84 +no\_defs - "

¿Podría alguien comprobar si mis cadenas proj4 son correctas (especialmente WGS84 para Google Earth) o sugerir otras causas del error?

Answer 1

Su proyección de Gauss-Krueger utiliza +datum=potsdam . Hasta 2012, esto estaba codificado en proj4 a un valor muy poco preciso utilizando una transformación de 3 parámetros.

Encontrará valores más exactos para las transformaciones de 7 parámetros en este tema:

http://forum.openstreetmap.org/viewtopic.php?id=12723

Existe una transformación ntv2-grid aún mejor aquí (toma el binario) que tiene que estar en la misma carpeta que la aplicación y los datos, a menos que se especifiquen los nombres de ruta completos.

Para comparar los diferentes valores posibles, hice este archivo por lotes de prueba:

echo epsg31467-epsg4326 >out.txt
cs2cs +init=epsg:31467 +to +init=epsg:4326 31467.txt >>out.txt

echo proj-Definition epsg >>out.txt
cs2cs +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs  +to +init=epsg:4326 31467.txt >>out.txt

echo proj-definition Qgis >>out.txt
cs2cs +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +towgs84=582,105,414,1.04,0.35,-3.08,8.3 +units=m +no_defs +to +init=epsg:4326 31467.txt >>out.txt

echo proj-definition nadgrid >>out.txt
cs2cs +proj=tmerc +lat_0=0 +lon_0=9 +x_0=3500000 +y_0=0 +k=1.000000 +ellps=bessel +units=m +nadgrids=./BETA2007.gsb +wktext +to +init=epsg:4326 31467.txt >>out.txt

echo epsg31467-epsg3785 >>out.txt
cs2cs +init=epsg:31467 +to +init=epsg:3785 31467.txt >>out.txt
echo proj-definition Qgis >>out.txt
cs2cs +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +towgs84=582,105,414,1.04,0.35,-3.08,8.3 +units=m +no_defs +to +init=epsg:3785 31467.txt >>out.txt

echo proj-definition nadgrid >>out.txt
cs2cs +proj=tmerc +lat_0=0 +lon_0=9 +x_0=3500000 +y_0=0 +k=1.000000 +ellps=bessel +units=m +nadgrids=./BETA2007.gsb +wktext +to +init=epsg:3785 31467.txt >>out.txt

echo epsg31467-epsg3857 >>out.txt
cs2cs +init=epsg:31467 +to +init=epsg:3857 31467.txt >>out.txt
echo proj-definition Qgis >>out.txt
cs2cs +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +towgs84=582,105,414,1.04,0.35,-3.08,8.3 +units=m +no_defs +to +init=epsg:3857 31467.txt >>out.txt

echo epsg31467-epsg900913 >>out.txt
cs2cs +init=epsg:31467 +to +init=epsg:900913 31467.txt >>out.txt
echo epsg31467-proj900913 >>out.txt
cs2cs +init=epsg:31467 +to +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +over +no_defs 31467.txt >>out.txt
echo proj31467-proj900913 >>out.txt
cs2cs +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +towgs84=582,105,414,1.04,0.35,-3.08,8.3 +units=m +no_defs +to +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +over +no_defs 31467.txt >>out.txt

con cualquier par de coordenadas Gauss-Krüger lon-lat de muestra en 31467.txt