¿Cómo puedo utilizar los datos proyectados de GK + ETRS89 en QGIS o Proj?

Question

Tengo un conjunto de datos de Hamburgo, Alemania, en el llamado "Lagestatus 320" (LS320). Los metadatos lo describen como "Gauss-Krueger mit ETRS89 [Lagestatus 320]". Las coordenadas son como 3565575.3,5936627.0, son métricas.

La base de datos epsg no dispone de este sistema de referencia de coordenadas. La opción más cercana es EPSG:31467, también conocida como "DHDN / zona 3 de Gauss-Kruger". Por supuesto, esto utiliza el elipsoide equivocado.

Si asigno eso a mis datos, las cosas terminan a unos 600 metros al norte de su ubicación correcta.

¿Cómo puedo utilizar correctamente estos datos en QGIS o en cualquier otro software basado en proyectos?

Answer 1

Tienes que crear tu propio CRS personalizado utilizando el elipsoide correcto.

Las proyecciones de DHDN utilizan +ellps=bessel +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 .

Las proyecciones de ETRS89 utilizan +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 .

La cadena de proyecto completa para EPSG:31467 es:

+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0
+ellps=bessel +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7
+units=m +no_defs

Así que cambia lo del elipsoide por los parámetros de ETRS89 así:

+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0
+ellps=GRS80 +towgs84=0,0,0,0,0,0,0
+units=m +no_defs