¿Cómo revertir la georeferencia con QGIS?

Tengo un mapa con una proyección "extraña" (básicamente está pintada a mano, no es exacta). Puedo georreferenciarlo fácilmente en QGIS usando el complemento Georeferencer y la transformación Thin Plate Spline con un conjunto de algunos puntos de control 2000-3000 en una simple proyección WLS 84 (EPSG: 4326) "latlong" y procesarla / mejorarla aún más con otros datos, como de Natural Earth o OpenStreetMap.

Ahora me gustaría realizar la operación inversa : transformar una de esas capas ráster (o capas vectoriales rasterizadas) nuevamente en la cuasiproyección pintada a mano, idealmente usando el mismo conjunto de puntos de control, para que las imágenes coincidan cuando se superponen Un editor gráfico. ¿Cómo lo hago?

qgis georeferencing Martin Sojka
fuente

Situación interesante. Tal vez no georeferencia el mapa pintado a mano en absoluto. Use el mapa pintado a mano como 'verdad del terreno' y obtenga sus puntos de control según el sistema de coordenadas local (por ejemplo, la esquina superior izquierda del mapa pintado a mano es 0,0 y cada píxel es una unidad). Luego 'georreferencia' sus rásteres a este sistema de coordenadas local.

pensamientos espaciales

@spatialthoughts: esto funcionaría bastante bien para los rásteres, sin embargo, también perdería la precisión de las capas vectoriales y la posibilidad de hacer cálculos con ellas (área de una entidad, longitud de un segmento, etc.). Ese no es un problema insuperable, por supuesto. Ser capaz de revertir la georreferenciación simplemente no tendría esas limitaciones.

Martin Sojka

@ Martin ¿Quiere los cálculos vectoriales en WGS84 o en las coordenadas de imagen pintadas a mano? Si lo primero, parece que su necesidad de transformar vectores en las coordenadas pintadas a mano está totalmente separada de su necesidad de hacer cálculos en el vector, y yo recomendaría la sugerencia de los pensamientos espaciales. Solo significa que debe realizar estas tareas en proyectos QGIS separados.

Jeff G

En general, no se garantiza que exista un inverso para todas las transformaciones. Con las transformaciones polinómicas, si pudiera hacer que el complemento registre los coeficientes que utilizó, podría calcular el inverso matemáticamente. Con la ranura de placa delgada, quizás también sea posible una inversa: la función GDALCreateTPSTransformer toma un argumento "bReversed": gdal.org/gdal__alg_8h.html#a245802b88a8126c138d24febe6c9822a

Jeff G

Respuestas:

Si estaba usando ArcGIS, usa la barra de herramientas de ajuste espacial. Sin embargo, QGIS ofrece una funcionalidad similar a través de complementos.

Vaya al administrador de complementos e instale 'Transformaciones afines'

Instale también 'Vector Bender'

Página de inicio: https://github.com/trenneman/qgsAffine ;

https://plugins.qgis.org/plugins/VectorBender/

Como ya tiene los puntos de control, debería poder transformar sus datos a la "proyección" original. Cargue sus puntos de control en el complemento y realice la transformación. Todavía no lo he intentado, así que quizás tengas que jugar con los puntos de control. El complemento de transformaciones afines usa constantes para realizar la transformación. Como ya tiene puntos de control, querrá usar el doblador de vectores para realizar una 'hoja de goma'.
Vector bender no está en el repositorio de complementos, por lo que deberá instalarlo manualmente agregándolo a C: \ Users.qgis2 \ python \ plugins \ VectorBender.

Segunda opción (puede que solo funcione para rásteres, no se ha probado con vectores): en el georreferenciador, puede guardar sus puntos de control en un archivo. También puede cargar puntos de control desde un archivo.

Aquí está lo que tú necesitas hacer:

Guarde el archivo cpg. Cambie la extensión a .csv. Abra el archivo en excel. Intercambie los valores en las columnas hacia / desde. Salvar. Vuelva a cambiar la extensión a .points.

Una vez que haya configurado los puntos de control, puede transformar de un lado a otro utilizando los puntos de control y su método de transformación preferido.

jbalk
fuente

¿No es la barra de herramientas de ajuste espacial parte de Arcgis Desktop? No estoy seguro de que haya algo similar en QGIS.

AndreJ

@AndreJ Por favor, vea mis ediciones.

jbalk

Tuve que hacer esto una vez para una capa de trama y me llevó un tiempo ver que no necesitaba la inversa. Para que las imágenes coincidan en la extraña proyección original, escriba un guión para crear una cuadrícula en esa proyección, correspondiente a los píxeles en su editor de gráficos. Transforme hacia adelante para cada punto en la cuadrícula para recoger el valor de la capa que tiene en WGS 84. (similar al comentario de pensamientos espaciales )

Nick Mulgan
fuente