La georreferenciación de una imagen con GDAL, basado en QGIS Georeferencer Plugin

Question

Me gustaría georreferenciar una imagen con GDAL. La imagen ya tiene una referencia espacial conjunto (EPSG:31276). Para la georreferenciación, tengo una lista de 4 puntos de control terrestre (también en EPSG:31276).

El uso de la QGIS Georeferencer Plugin, puedo correctamente georeferenciar la imagen. Yo, a continuación, seleccione "Generar GDAL Script", que se traduce en las siguientes GDAL comandos:

gdal_translate -of GTiff -gcp 6.55253e+06 -5.07402e+06 6.55273e+06 5.07395e+06 -gcp 6.55246e+06 -5.0711e+06 6.55263e+06 5.07106e+06 -gcp 6.55477e+06 -5.07382e+06 6.55511e+06 5.07368e+06 -gcp 6.55483e+06 -5.07197e+06 6.55525e+06 5.07197e+06 "C:/gis_sample_data/sample_image/image1.tif" "C:/Users/XYZ/AppData/Local/Temp/image1.tif"

gdalwarp -r near -order 1 -co COMPRESS=NONE  "C:/Users/XYZ/AppData/Local/Temp/image1.tif" "C:/gis_sample_data/sample_image/image1_georef.tif"

Los comandos se ejecutan sin un error, pero la imagen resultante no es correctamente proyectado. Aquí está una pequeña parte de la gdalinfo "C:/gis_sample_data/sample_image/image1.tif" (el original de la imagen):

Corner Coordinates:
Upper Left  ( 6551897.000, 5074277.000) ( 18d40' 4.21"E, 45d48'39.73"N)

mientras que gdalinfo "C:/gis_sample_data/sample_image/image1_georef.tif" (resultado de la GDAL de llamada) muestra

Corner Coordinates:
Upper Left  ( -534350.794,  156938.268) ( 34d51'24.74"W,  0d52'47.63"N)

¿Cómo tendría que modificar el GDAL los comandos correctamente georeferenciar la imagen?

La imagen de muestra es parte de la Baranja Hill conjunto de datos, que puede ser descargado desde aquí. Los puntos de control son estas:

mapX,mapY,pixelX,pixelY,enable
6552733.64613180514425039,5073946.04477077350020409,6552533.84482758585363626,5074021.48275862075388432,1
6552633.32449856773018837,5071059.01110315229743719,6552461.44827586226165295,5071104.32758620660752058,1
6555111.64040114637464285,5073678.52041547279804945,6554765.36206896603107452,5073821.32758620660752058,1
6555249.11819484271109104,5071965.62141833826899529,6554833.5,5071973.08620689623057842,1

Cuando se copian en un editor de texto y guardarlo con el .puntos de extensión, pueden ser importados en el QGIS Georreferencia del Plugin.

EDIT: Cambiado negativos a valores positivos en gdal_translate. Resultados de la gdalinfo:

Corner Coordinates:
Upper Left  ( -591887.923,  223739.838) ( 35d 8'28.93"W,  1d14'58.74"N)

Por cierto, también probé agregando -t_srs EPSG:31276, -s_srs EPSG:31276 y -a_srs EPSG:31276 , y diversas combinaciones de los mismos para la gdal comandos, pero no soluciona el problema.

Answer 1

Las coordenadas en el gdal_translate comando debería ser todos positivos.

El GCP lista, así que creo que esto podría ser un bug de la GDAL script creator.

---

EDITAR

Me podría replicar su error, y creo que el caso de uso no es la intención de uso de BPC.

En "normal" en la georeferenciación, el GCP puntos de píxel del mapa de coordenadas de una imagen sin referencia a la proyectada o coordenadas geográficas. El píxel de coordenadas son relativas a la esquina superior izquierda.

En su caso, tomar coordenadas proyectadas desde un georeferenciados, tif, y compararlos para "corregir" las coordenadas. Las coordenadas en el lado izquierdo no son relativas a la esquina superior izquierda, pero el origen de la proyección. Si usted comienza el proceso de proceso, esta información será distroyed y reemplazado por el GCP información. Esto no parece funcionar.

En el GDAL archivo de proceso por lotes, usted puede agregar -of VRT para obtener vrt en lugar de archivos tif. Puede abrir estos vrt en un editor de texto y ver que el juego va mal, muy mal.