Capas en QGIS que deberían alinearse no lo hacen (pero sí en FME)

Question

He visto numerosas preguntas de esta naturaleza y todas sugieren la asignación incorrecta de SRI en el camino. Creo que he probado todas las combinaciones de valores CRS relevantes, pero sigo sin poder cargar dos capas vectoriales concretas de forma que se alineen. Debo estar perdiendo un truco en alguna parte.

Las dos capas provienen de OS: AddressBase Plus y Mastermap Topography. En concreto, estoy intentando alinear el conjunto de datos de direcciones con el conjunto de datos de áreas topográficas.

Puedo cargarlos en el inspector de datos de FME y los puntos de dirección se alinean con los edificios, pero al cargarlos en QGIS aparecen unos metros fuera.

Cuando compruebo el CRS en FME se reporta como EPSG:27700 para ambos conjuntos de datos, aunque extrañamente el Sistema de Coordenadas para la dirección está etiquetado como EPSG:27700 mientras que las Áreas Topográficas están etiquetadas como BritishNatGrid (comparando el WKT muestra que son idénticas, sin embargo).

Cuando cargo estos conjuntos de datos en QGIS, sin especificar el CRS, lo establece en 27700 para las Áreas Topográficas, pero 4326 para las Direcciones. No entiendo por qué.

He probado todas las combinaciones (creo que lo he hecho) de ajustar el CRS del proyecto a diferentes valores, asignando diferentes valores a cada capa, y activando y desactivando el OTF, pero sin resultado.

La mayoría de las cosas que prueban ponen las direcciones frente a la costa de África Occidental (supongo que en lat/lon 0, 0), y las Áreas Topográficas no aparecen por ningún lado. La única combinación que las pone, aunque sea vagamente, en la misma vecindad es con las Direcciones 4326 y las Áreas Topográficas 27700, pero siguen desalineadas por unos pocos metros.

¿Puede alguien sugerir un proceso sencillo para cargar estos conjuntos de datos en QGIS de manera que se alineen?

Aquí está el WKT de FME (que se alinea correctamente). Es el mismo para ambos archivos:

PROJCS["British National Grid (ORD SURV GB)", GEOGCS["OSGB 1936", DATUM["OSGB_1936", SPHEROID["Airy, 1830",6377563.396,299.3249612664953, AUTHORITY["EPSG","7001"]], AUTHORITY["EPSG","6277"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4277"]], PROJECTION["Transverse_Mercator"], PARÁMETRO["latitud_de_origen",49], PARÁMETRO["meridiano_central",-2], PARAMETER["scale_factor",0.9996012717], PARAMETER["false_easting",400000], PARAMETER["false_northing",-100000], UNIDAD["METRO",1], AUTHORITY["EPSG","27700"]]

Los resultados de ogrinfo son:

Metadatos: DESCRIPCIÓN=Ordnance Survey (c) Crown Copyright. Todos los derechos reservados, 2016 y producido por GeoPlace

Nombre de la capa: Dirección Geometría: Punto Recuento de rasgos: 41206 Extensión: (-0.135610, 51.323120) - (-0.062560, 51.369150) Capa SRS WKT: GEOGCS["ETRS89", DATUM["Sistema_de_referencia_terrestre_europeo_1989", SPHEROID["GRS 1980",6378137,298.257222101, AUTHORITY["EPSG","7019"]], TOWGS84[0,0,0,0,0,0,0], AUTHORITY["EPSG","6258"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4258"]]

Metadatos: DESCRIPCIÓN=Ordnance Survey, (c) Crown Copyright. Todos los derechos reservados, 2016-04-05

Nombre de la capa: TopographicArea Geometría: Polígono Recuento de rasgos: 99404 Extensión: (529039.110000, 159261.400000) - (535629.000000, 165876.100000) Capa SRS WKT: PROJCS["OSGB 1936 / British National Grid", GEOGCS["OSGB 1936", DATUM["OSGB_1936", SPHEROID["Airy 1830",6377563.396,299.3249646, AUTHORITY["EPSG","7001"]], TOWGS84[446.448,-125.157,542.06,0.15,0.247,0.842,-20.489], AUTHORITY["EPSG","6277"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4277"]], PROJECTION["Transverse_Mercator"], PARÁMETRO["latitud_de_origen",49], PARÁMETRO["meridiano_central",-2], PARAMETER["scale_factor",0.9996012717], PARAMETER["false_easting",400000], PARAMETER["false_northing",-100000], UNIT["metro",1, AUTHORITY["EPSG","9001"]], EJE["Orientación",ESTE], EJE["Northing",NORTH], AUTHORITY["EPSG","27700"]]

Answer 1

El alcance de la Adress de la base de datos es `(-0,135610, 51,323120) - (-0,062560, 51,369150). Esto es seguramente NO EPSG:27700 British National Grid (usando metros), más bien algo en grados.

Ogrinfo informa de EPSG:4258, que es casi idéntico a EPSG:4326, y parece razonable para la extensión dada. Como alternativa, podría ser EPSG:4277 OSGB36 en grados, con un cambio de datum a WGS84. Pero esto está aún más desplazado a los datos EPSG:27700.

El Topographic Area extensión bien podría ser EPSG:27700 en metros.

---

ACTUALIZACIÓN

La página de descarga del sistema operativo ofrece un conjunto de datos de muestra. Mirando las columnas CSV sin nombre, contienen tanto OSGB 27700 (en los campos 8 y 9) como grados (en los campos 10 y 11). Así que al importar como texto delimitado, depende de usted qué columnas toma como coordenadas. Ambas no coinciden exactamente (grados en verde y metros en rojo):

De la fuente de datos GML (que también tiene ambos conjuntos de coordenadas), QGIS está tomando los valores de grados, mientras que FME podría tomar los valores de metros.

Ordnance Survey utiliza un cambio de datos diferente al que utiliza QGIS para convertir entre 27700 y WGS84. Utilizan un archivo de cambio de cuadrícula, que puede descargar y guardarlo en algún lugar de su ordenador. A continuación, puede construir un CRS personalizado para EPSG:27700 basado en esa transformación:

 +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +nadgrids=F:\Karten\gdal\ntv2\OSTN02_NTv2.gsb +units=m +no_defs

Aplicando eso a los datos de 27700 en metros (con la ruta correcta al archivo), los puntos rojos se alinean perfectamente con las coordenadas de grados verdes.

Answer 2

Transformación de la red nacional OSTN02 (ETRS89-OSGB36) Para hacer frente a las distorsiones en el OSGB36 TRF, se realizan diferentes transformaciones necesario en diferentes partes del país. Por esta razón, la transformación de datos transformación de datos estándar entre OSGB36 y ETRS89 no es una simple transformación de datos Helmert.

En su lugar, Ordnance Survey ha desarrollado una transformación de estilo "hoja de goma" que funciona con una rejilla de transformación expresada en coordenadas de coordenadas de noreste y de longitud. Las cuadrículas de desplazamientos de noreste y east entre ETRS89 y OSGB36 cubren Gran Bretaña con una resolución de un kilómetro. A partir de estas cuadrículas, se ha calculado un desplazamiento de noreste y east para cada punto que se va a transformar se se obtiene mediante una interpolación bilineal. Esto se denomina "National Grid Transformation OSTN02, y está disponible gratuitamente en paquetes de software en el sitio web de Ordnance Survey GPS,

https://www.ordnancesurvey.co.uk/business-and-government/products/os-net/index.html

La Transformación de la Red Nacional no sólo hace frente al cambio de datum entre los dos sistemas de coordenadas, sino también a las distorsiones del TRF en la red de triangulación OSGB36

el tipo Helmert limitado a aplicaciones a cinco metros y más niveles de precisión. Esta transformación elimina la necesidad de estimar transformaciones Helmert locales entre ETRS89 y OSGB36 para ubicaciones particulares. Dado que la Transformación de la Red Nacional funciona con coordenadas de este y norte, otros tipos de coordenadas ETRS89 (cartesianas 3D o de latitud y longitud) deben ser convertidas primero a a coordenadas de este y norte. Esto se hace utilizando la misma proyección cartográfica que se utiliza para la National Grid (véase la sección 7), con la salvedad de que que se utiliza el elipsoide GRS80 en lugar del elipsoide Airy. Los parámetros y fórmulas de parámetros y fórmulas necesarios para obtener estas orientaciones ETRS89 se encuentran en los anexos A-C. Después de la transformación, las las orientaciones y nortes resultantes de la National Grid se pueden volver a convertir a latitud y longitud (esta vez utilizando el elipsoide de Airy) si es necesario. si es necesario.

fuente: http://www.bnhs.co.uk/focuson/grabagridref/html/OSGB.pdf

Answer 3

Cuando QGis abre el archivo AddressBase no puede determinar el SRS y por lo tanto adivina que es EPSG:4326 (este es el comportamiento por defecto) lo cual es por supuesto incorrecto. La solución sencilla es hacer clic con el botón derecho del ratón en esa capa y seleccionar set Layer CRS y elija EPSG:27700 y haga clic en Aceptar.

Sus dos capas se alinearán ahora.

Para facilitarte la vida en el futuro (suponiendo que trabajas principalmente con datos del Reino Unido) puedes establecer que la proyección por defecto para las capas desconocidas sea OSGB yendo a ajustes->opciones->crs y eligiendo allí epsg:27700.