¿Mejorando los resultados de georreferenciación? [cerrado]

13

Antecedentes Esta es mi segunda pregunta relacionada con la georreferenciación de mapas raster desnudos para volver a visualizarlos en diferentes sistemas de coordenadas y en conjunto con otras capas de datos. La pregunta anterior es Convertir una imagen de mapa arbitraria libre de metadatos en un proyecto QGIS

Problema Mi objetivo es georeferenciar este mapa:

Estepa euroasiática, Encyc.  Brit.?

Esto no parece ser Plate-Carrée. Entonces, en QGIS, creé varios puntos de control razonables, que para completar he adjuntado en la parte inferior [ref: 1]. Proporciono a QGIS Georeferencer el mismo SRS de destino que mi archivo de proyecto, EPSG: 4326. Obtengo resultados excepcionalmente pobres con Helmert y las transformaciones polinómicas, pero obtengo una imagen razonable con una ranura delgada de placa (lo que hace que el geoestimado resultante pase por mis puntos de control). Sin embargo, incluso este resultado es pobre, por ejemplo, en latitudes más altas (ver la costa rusa al norte de Japón). Esta es una captura de pantalla de mi pantalla QGIS usando un fondo Natural Earth.

Resultado de georreferencia QGIS, estría de placa delgada

Camino alternativo Intenté un ejercicio similar con la herramienta mucho más fácil de usar en MapWarper: vea el resultado y los puntos de control en http://mapwarper.net/maps/758#Preview_Map_tab donde obtengo resultados más pobres (probablemente debido al hecho que agregué menos puntos de control).

Preguntas en pocas palabras

  1. ¿Hay algún truco que me falte para obtener una buena georreferencia?
  2. ¿Es esta proyección reconocible al instante?
  3. En el Sistema de coordenadas desconocido en el dibujo anterior , gdaltransformse sugiere transformar varios puntos de coordenadas en un SRS de varios objetivos, con el objetivo de descubrir los parámetros de proyección utilizados para generar el mapa original. Intenté algo como esto: después de guardar mi lista de puntos QGIS, realicé un procesamiento de cadenas para obtener una lista de long / lats separados por espacios cat eurasian-steppe-gcp.points | tail -n+2 | cut -d, -f1-2 | sed 's/,/ /'> tmp.txte ingresando el archivo resultante en gdaltransform: gdaltransform -s_srs EPSG:3785 -t_srs EPSG:4326 < tmp.txty cambiando el s_srsyt_srsbanderas (el proyecto usa EPSG: 4326). Sé que estoy disparando en la oscuridad, con la esperanza de tener suerte, así que no me sorprendió cuando no pude entender las salidas. ¿Alguien puede ampliar cómo usaría este método para encontrar la mejor estimación de los parámetros de proyección y proyección del mapa fuente? Mi opinión detrás de esto es que, en lugar de perder el tiempo colocando una miríada de puntos de control para una buena georreferencia, ¿podría ser más fácil obtener una georreferencia casi perfecta con menos puntos de control, simplemente recorriendo todos los sistemas de coordenadas comunes? ¿Implica la validación cruzada de cada punto contra todos los demás, para cada CRS bajo prueba?

Me gustaría conocer este algoritmo o la georreferenciación para poder automatizar el proceso. Me encuentro con este problema todo el tiempo y hasta que los creadores de contenido dejen de tratar sus mapas como creaciones únicas que nunca se integrarán. con otro contenido, no espero parar.

Referencias

[ref: 1] archivo QGIS GCP:

mapX,mapY,pixelX,pixelY,enable
142.632649100000009,54.453595900000003,505.941176470588232,-95.220588235293974,1
154.934252200000003,59.559921699999997,536.411764705882206,-52.779411764705742,1
80.080158100000006,9.657192300000000,291.558823529411711,-322.661764705882206,1
10.448442600000000,57.819128900000003,21.676470588235190,-103.926470588235134,1
34.007173000000002,27.761438299999998,101.117647058823422,-244.852941176470466,1
50.950890399999999,11.862196600000001,171.852941176470495,-313.955882352941046,1
29.713217199999999,60.024133200000001,90.779411764705799,-92.499999999999829,1
60.000000000000000,0.000000000000000,208.308823529411683,-362.382352941176350,1
69.867506500000005,66.639146199999999,224.088235294117567,-33.191176470588061,1
27.276107100000001,71.049154799999997,89.147058823529306,-21.764705882352814,1
140.000000000000000,0.000000000000000,536.955882352941217,-362.926470588235190,1
20.000000000000000,0.000000000000000,43.441176470588132,-362.926470588235190,1
20.196882700000000,31.243024100000000,47.249999999999901,-231.794117647058698,1
9.171861099999999,42.848309999999998,8.073529411764603,-175.205882352941046,1
131.955786100000012,43.196468600000003,481.999999999999943,-162.691176470588090,1
73.813303700000006,45.169367200000003,256.735294117646959,-161.602941176470438,1
50.602731800000001,44.589102900000000,168.044117647058727,-167.588235294117510,1
121.394975900000006,18.941421099999999,455.882352941176407,-284.029411764705742,1
103.987047000000004,1.417439300000000,389.499999999999943,-357.485294117646959,1
109.325478599999997,55.962283100000001,380.249999999999943,-98.485294117646902,1
31.454010100000001,46.562001500000001,95.132352941176379,-158.882352941176322,1
43.639560299999999,68.844150499999998,137.573529411764611,-40.264705882352814,1

No actualización

Análisis de van der Grinten Escribí una herramienta de Python para ajustar los GCP a cualquier proyección que Proj4 admita (a través de Pyproj) y la apliqué a las proyecciones sugeridas en las respuestas. El código fuente (algo descuidado, pido disculpas por adelantado), así como los GCP actualizados están disponibles en https://github.com/fasiha/steppe-map

Van der Grinten tiene solo 1 parámetro para sintonizar, y aquí está la imagen resultante (usando la última imagen de Britannica, muchas gracias a ellos por dar un mapa tan actualizado y de alta resolución (aunque todavía carece de datos de proyección)).

Van der Grinten fit

Van der Grinten tiene un error relativo de 0.035 entre los GCP y los puntos de mejor ajuste, que es el peor del grupo que probé, y la superposición de la costa lo confirma cualitativamente.

(Puede ayudar si abre esta imagen en su propia pestaña, es bastante alta resolución. También verá flechas verdes que indican los puntos georreferenciados (deben coincidir con puntos de referencia significativos en la imagen), así como flechas rojas que indican dónde están esos puntos están ajustados a (deben coincidir con los mismos puntos de referencia en la superposición de la costa) --- la desviación entre los dos puede ayudar al ojo a ver las diferencias entre la imagen y el ajuste).

Análisis de Albers de área igual Intentando lo mismo con la proyección de Albers de área igual (que es lo mismo que "Albers conformal Conic"? Perdón por mi ignorancia). Este ajuste, que implica un ajuste de parámetro de 4 dimensiones, es mejor, con un error relativo de 0.025, pero no obstante se ve bastante pobre.

Ajuste de área igual de Albers

Análisis de las proyecciones de Robinson y Eckert V Encajé en varias proyecciones pseudocilíndricas respaldadas por Pyproj (todo lo que pude encontrar que tenía un parámetro libre) y descubrí que las proyecciones de Robinson y Eckert V hicieron el "mejor" en términos de error relativo entre los PCG y los puntos ajustados, ambos con errores relativos de 0.015.

Aquí está el Robinson:

Robinson en forma

Y aquí está el Eckert V.

Eckert V fit

Tenga en cuenta las desviaciones de la costa ajustada de la costa de la imagen. Creo que con esto puedo concluir que el mapa no es ninguno de los anteriores.

Tripel de Winkel: ¿posible solución?

Después de probar secuencialmente todas las proyecciones en este manual de Proj de 1990 (actualizado en 2003) ftp://ftp.remotesensing.org/proj/OF90-284.pdf Finalmente llegué a la proyección de tripel de Winkel. Esto produce los errores cuantitativos más bajos (0.011) y el litoral es uniformemente bastante bueno (o equivalente, uniformemente un poco malo). Leí que esta es la proyección de la National Geographic Society, lo que significa que es famosa, y esto agrega peso a la candidatura de esta proyección para el mapa de Britannica. El SRS instalar: +units=m +lon_0=47.0257707403 +proj=wintri.

Winkel tripel fit

(Disculpas por cambiar el color de la costa a gris. Si esto ofende a alguien, puedo producir una versión azul).

Intentaré modificar mis GCP para intentar reducir el error.

Ahmed Fasih
fuente
Esa debería ser una proyección de Mercator (los paralelos no son equidistantes).
Bhell
1
Me equivoco: Mercator, por supuesto, debería tener meridianos paralelos. Tal vez van der Grinten ?
Bhell
¿Has intentado agregar algunos puntos de control en el extremo norte? Todos los demás lugares se ven bien, pero parece que falta más información.
oscuro
1
@bhell Van der Grinten es una suposición excelente, pero no encaja del todo. La proyección parece ser la misma que en este mapa (que tiene una retícula), también realizada recientemente por la Enciclopedia Británica. Es una especie de proyección pseudo policónica en un aspecto ecuatorial centrado en el primer meridiano, pero no lo reconozco. Investigar en la Enciclopedia Británica es una buena idea.
whuber
2
Ah, entonces la imagen no se puede convertir en una proyección de Van der Grinten (con la que estoy familiarizado porque se usa en un mapa 'al revés' en mi pared :) a través de cualquier transformación afín. ¡Bueno saber!
Ahmed Fasih

Respuestas:

3

Desde un primer análisis, podemos ver que solo el ecuador y el meridiano central (que pasa aproximadamente por Moscú, aproximadamente 37.36 E de longitud) son líneas rectas, mientras que todos los otros paralelos y meridianos son líneas curvas. Entre las proyecciones globales, estas características sugieren probar el Van der Grinten (como bhell hizo en un comentario de la pregunta), es decir, la proyección del mapa mundial estándar de la National Geographic Society.

Así que configuré el proyecto SR en QGIS como SRS personalizado World_Van_der_Grinten_I (Moscow):

+proj=vandg +lon_0=37.36 +x_0=0 +y_0=0 +R_A +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

deriva directamente de EPSG:54029( World_Van_der_Grinten_I), cambiando solo la longitud del meridiano central. Luego intenté georreferenciar la imagen usando la transformación proyectiva con los siguientes GCP elegidos en una retícula imaginaria que pasa por las marcas disponibles:

mapX,mapY,pixelX,pixelY,enable
-1930600.82502013398334384,-91.97837987072125543,43.52016129032256231,-362.54032258064518146,1
11412132.24015056900680065,-91.97837987233651802,537.54173867127497033,-362.73421498975926625,1
11804723.15208163484930992,8693410.00194966979324818,551.40707210100163138,-43.85864138853239069,1
-1651013.74878141912631691,7637895.72632359713315964,54.16126727366525984,-80.84230997717837397,1

Este es el resultado en la costumbre World_Van_der_Grinten_I (Moscow)con GCP (límite de administración del país mundial en rojo):

ingrese la descripción de la imagen aquí

Esta georreferencia podría mejorarse aún más eligiendo PCG más precisos (lo más difícil es la ausencia de marcas de verificación o una retícula), sin embargo, parece que el mapa contiene algunos artefactos, especialmente en Europa Central y del Norte. Por ejemplo, considere Alemania, que parece colapsar en referencia a Italia o la posición relativa de la isla Kolguyev (que se ajusta bastante bien) y la costa rusa (que no encaja en absoluto).

El mismo resultado reproyectado en WGS84( EPSG:4326):

ingrese la descripción de la imagen aquí

Antonio Falciano
fuente
Muchas gracias por tu consejo e iniciativa. No estoy familiarizado con las transformaciones proyectivas, así que descubriré más sobre ellas. Actualicé la publicación principal con un enfoque de ajuste de parámetros que intenta encontrar la proyección de Van der Grinten que mejor se ajusta a los GCP (estaba bastante cerca de 37.36) y estaba bastante insatisfecho con ese ajuste. Por favor, eche un vistazo y hágamelo saber.
Ahmed Fasih
Ah, veo que la transformación proyectiva en QGIS es una rotación + traducción, que debería ser fácil de probar con mis GCP.
Ahmed Fasih
En términos generales (y también en QGIS), la transformación proyectiva no es una rototraducción. La guía del usuario de QHOI en mi humilde opinión es incorrecta con respecto a este punto. Además, considere la posibilidad de que pueda haber algunos artefactos hechos a mano (?) En el mapa de baja resolución, como parece después de mi georreferencia. ASAP trataré de georreferenciar el de alta resolución.
Antonio Falciano
Gracias por la información, actualizaré mi QGIS y probaré la georreferencia de transformación proyectiva. (También intentaré mirar el código fuente de QGIS para ver qué podría estar haciendo). Si desea ahorrar algo de tiempo y usar mis GCP, están en github.com/fasiha/steppe-map/blob /master/gcp.txt (guardado por QGIS).
Ahmed Fasih
Consulte la publicación actualizada para una posible solución: Winkel tripel.
Ahmed Fasih
1

hay una fuerte distorsión en el este ("subiendo"): supongo que es una cónica conforme de Albers. Entonces el meridiano vertical parece estar a 40 ° (-> su meridiano central) y el paralelo horizontal probablemente alrededor de 40 ° N. Luego, debe medir su coordenada XY desde el eje en la intersección del meridiano 40 y el paralelo 40, luego intentarlo ...

Nota: Los paralelos no son paralelos, por lo que puede excluir las proyecciones cilíndricas (Mercator) y pseudo-cilíndricas (sinusoidales, Eckert)

radouxju
fuente
La falta de marcas en ángulo en los lados E y W del mapa sugiere que no es cónica.
Martin F
Esta es una buena observación, pero hay algo extraño si observa las marcas de verificación: no están en ángulo, pero cuando las une con una línea recta, no son líneas paralelas ni horizontales. Por lo tanto, deben ser curvos y las garrapatas deben estar en ángulo. De todos modos, si es conforme, probablemente no sea Lambert, ya que la convergencia al polo crearía distorsiones muy fuertes. Quizás Albers.
radouxju
Muchas gracias por el consejo, y lo siento, me tomó tanto tiempo intentar implementarlo. Actualicé la publicación principal con mi intento de ajustar los datos a una proyección de área igual de Albers de 4 dimensiones (¿es esto a lo que te refieres?) A los GCP que extraje. No estoy satisfecho con el ajuste, y encontré un par de otras proyecciones pseudocilíndricas (Eckert V y Robinson) que lo hacen (cuantitativamente) bastante mejor, pero aún no lo suficientemente cerca como para analizar el tema del mapa. ¿Tus pensamientos?
Ahmed Fasih