Raster reproyectado incorrectamente a OSGB(27700)

Question

Ubuntu 12.04, QGIS 1.8.0, GDAL/OGR 1.7.3, GEOS 3.2.2

Abra un db de Spatialite con la poli capa OS Mastermap (EPSG:27700);
Añadir OpenLayers Plugin (1.0) desde el repositorio oficial de QGIS;
Añadir una capa de Google Satellite (EPSG:3857) y hacer OS poly fills transparente;
Haz zoom y arrastra un poco el mapa hasta que se "asiente"; entonces ambos se alinearán correctamente;
Haga que la capa poly no sea visible y guarde el mapa como imagen (jpg + jpgw);
Volver a añadir el jpg guardado como capa (EPSG:3857) - sigue alineado;
Reproyectar (deformar) la capa jpg de EPSG:3857 a un Gtiff EPSG:27700;
Añade esta capa Gtiff - ya no alineada pero desplazada unos 6m al Norte.
???!!!???

(La transformación CRS sobre la marcha está activada)

He tenido este problema antes con el sabor inestable de UbuntuGIS QGIS, así que volví a la estable y todo estaba bien. Esta vez estoy recibiendo el 'cambio' en el UbuntuGIS QGIS estable y ninguna cantidad de borrar las bibliotecas (con Synaptic) y volver a instalar tiene hasta ahora deshacerse de este 'cambio'.

Cualquier ayuda o idea será bienvenida...

---

No hay problema para reproducir el error (gracias Underdark), pero si investigas un poco más, puede que acotes un poco más el problema.

La cadena de proyección proj4 CRS para OSGB/27700 en GDAL1.9.1 es:

+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +towgs84=375,-111,431,0,0,0,0 +units=m +no_defs

Lo mismo ocurre en GDAL1.7.3:

+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs

Aunque QGIS informa del CRS como 27700 con la misma cadena proj4 que el anterior GDAL1.7.3 probablemente está utilizando los parámetros posteriores de GDAL1.9.1 para la reproyección entre bastidores. Eso estaría bien si la imagen se alineara después de la reproyección, pero no es así.

Para ir un paso más allá, he creado una CRS definida por el usuario de QGIS con los mismos parámetros de proyección que la versión anterior de GDAL y la he reproyectado.

¿Son los EPSG:27700 parámetros de proyección el sospechoso probable en la última GDAL (que QGIS está utilizando) - o estoy ladrando al árbol equivocado?

Answer 1

Últimamente se ha debatido mucho sobre ello:

Error de reproyección de OSGB36 a WGS84

http://osgeo-org.1560.n6.nabble.com/OSGB-coordinates-to-WGS84-lat-lon-problem-td4965339.html

https://stackoverflow.com/questions/1426941/proj-4-library-and-osgb36

Si lo quieres muy exacto, usa el +nadgrids y el archivo .gsb mencionados. Aquí di algunos ejemplos de línea de comandos: ¿Cómo reproyectar un archivo raster en QGIS con transformación de datum?

Si desea que sus datos se ajusten entre sí, utilice siempre el mismo +towgs84 en lugar de +datum. En tu caso, actualiza el GDAL a 1.9.2.

A modo de aclaración: los elipsoides airy y WGS84 no coinciden exactamente. Por eso se inventaron los parámetros +towgs84. Pero estos cambian incluso a través del país. El +nadgrids tiene esto en cuenta.

---

EDITAR

Sólo para explicarlo un poco más en profundidad:

EPSG:27700 es originalmente sólo una definición de la proyección y el elipsoide utilizados. Las transformaciones de cambio de datum de un elipsoide a otro se enumeran por separado en la EPSG, con sus propios códigos EPSG. Hay varias transformaciones disponibles de OSB1936 a WGS84, dependiendo de la región o de la precisión deseada. Éstas son:

EPSG:1195 Accuracy: 21m
+towgs84=375,-111,431
used by GDAL 1.9.1,Qgis Lisboa/Master +towgs84:

EPSG:1196 Accuracy 10m
England
+towgs84=371,-112,434

EPSG:1197 Accuracy 21m
England, Wales, Isle of Man
+towgs84=371,-111,434

EPSG:1198 Accuracy 18m
Scotland
+towgs84=384,-111,425

EPSG:1199 Accuracy 35m
Wales
+towgs84=370,-108,434

EPSG 1314: Accuracy 2m 
+towgs84=446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894
used in proj +datum=OSGB1936
used in GDAL up to Version 1.7.3

EPSG:5339: Accuracy 0,03m
+nadgrids=OSTN02_NTv2.gsb

el datum de cuadrícula oficial puede obtenerse aquí: http://www.ordnancesurvey.co.uk/business-and-government/help-and-support/navigation-technology/os-net/ostn02-ntv2-format.html

Las transformaciones de tres parámetros fueron publicadas por el NIMA en el siglo pasado, calculadas a partir de las mediciones de un puñado de estaciones de satélite, mientras que los siete parámetros han sido publicados por Ordnance Survey más recientemente.

Por tanto, dependiendo de la versión de GDAL con la que esté georreferenciado un ráster, puede que obtenga un desplazamiento o puede que no.

---

EDIT2

el alivio está en camino:

http://trac.osgeo.org/gdal/changeset/25589

Así que buena oportunidad de tener de nuevo la transformación de siete parámetros en la próxima versión de Qgis ;-)