Georeferenciar imagen con 4 puntos

[hamipa]

Buenas tardes a todos/as. Estoy intentando georeferenciar una imagen raster de una ortofoto antigua en la que aparecen en sus cuatro esquinas las coordenadas X e Y. ¿habría alguna forma de hace una georeferenciación elástica? Lo he intentado con el comando alinear, pero solo me lo ajusta en dos puntos, y se me va en los otros dos. Gracias por su colaboración.

[LaUrA-28]

Para georeferenciar una imagen con 4 puntos de control necesitas una transformación de segundo orden o una transformación elástica (como spline), no una afín simple que solo usa 2 puntos. La herramienta Alinear de QGIS, por ejemplo, por defecto puede usar una transformación Helmert (lineal) que con solo 2 puntos no ajustará correctamente los 4.

La mejor solución actual es usar el complemento de georreferenciación integrado en QGIS. Debes abrir el panel Georreferenciador desde el menú Ráster. Añade tus puntos de control introduciendo manualmente las coordenadas X e Y para cada una de las cuatro esquinas. Es crucial que selecciones un tipo de transformación adecuado.

Para 4 puntos, debes elegir una transformación Polinómica de segundo orden. Esta transformación requiere un mínimo de 6 puntos, pero con 4 puntos puedes forzar su uso aunque la precisión no será óptima en toda la imagen. Si tu imagen tiene distorsión y solo tienes 4 puntos, la transformación de segundo orden puede funcionar. Añade los puntos con la mayor precisión posible.

Configura la transformación a Polinómica 2 y define los parámetros de resolución de salida. Luego haz clic en Comenzar georreferenciación. Si la imagen se deforma mucho, es porque los 4 puntos pueden no ser suficientes para un polinomio de segundo grado; en ese caso, si la imagen es plana y la deformación es mínima, podrías probar con una transformación Afín (lineal), que ajusta traslación, escala y rotación, pero esto asume que la imagen no tiene distorsión interna.

Una alternativa más robusta es el software GDAL. Puedes usar el comando gdal_translate con una opción -gcp para añadir cada punto de control. Necesitarás al menos 3 puntos para una transformación afín y más para polinómica. Luego usas gdalwarp con el método de transformación adecuado, como -tps (placas delgadas o spline) que es una transformación elástica perfecta para 4 puntos, ya que ajusta exactamente todos los puntos dados y deforma el resto de la imagen de manera suave. El comando sería primero gdal_translate -gcp pixelX pixelY coordX coordY imagen_entrada.tif imagen_con_gcp.tif para cada punto. Después ejecutas gdalwarp -tps imagen_con_gcp.tif imagen_salida_georef.tif.

En resumen, la forma más directa es desde QGIS usando el Georreferenciador con transformación Polinómica 2 o, si necesitas un ajuste exacto a los 4 puntos, usar GDAL desde la línea de comandos con el método de placas delgadas (TPS).