Georeferenciación de imágenes aéreas

Advertencia

Este tutorial está ahora obsoleto. Una versión nueva y actualizada está disponible en Georeferenciar Imágenes Aéreas (QGIS3)

En el tutorial: doc: georeferencing_basics cubrimos el proceso básico de georreferenciación en QGIS. Ese método implica leer las coordenadas del mapa escaneado e ingresarlas manualmente. Sin embargo, en ocasiones, es posible que las coordenadas no estén impresas en el mapa que se está intentando georeferenciar. En este caso, se puede usar otra fuente de datos georreferenciados. En este tutorial aprenderá a utilizar fuentes de datos abiertos existentes en el proceso de georreferenciación.

Información general del projecto

Vamos a georeferenciar una imágen de alta resolución de globo usando coordenadas de referencia de OpenStreetMap

Otras habilidades que aprenderás

  • ¿Cómo bajar imágenes de súper alta resolución de dominio público?

  • Por medio del complemento OpenLayers de QGIS.

  • Cómo hacer conversiones entre proyecciones diferentes usando la herramienta cs2cs en la lína de comandos.

  • Por medio de una capa ya georeferenciada para ingresar puntos GCP en la herramienta Georeferenciador.

  • Establecer un valor personalizado sin-datos para una capa.

Obtener los datos

In this tutorial, we will be using some gorgeous kite and balloon imagery collected by The Public Laboratory. They make the georeferenced versions of the images also available, but we will download a non-georeferenced JPG image and go through the process of georeferencing it in QGIS. If you like the imagery they provide, you can explore it in Google Earth as well.

Baje la imágen JPG Washington Square Park, New York. Puede hacer click derecho en el botón Save link as….

Para su comodidad, puede descargar directamente una copia del conjunto de datos del siguiente enlace:

newyorkcity-washingtonsquarepark.jpg

Procedmiento

  1. Para este tutorial, estaremos usando la capa OpenStreetMap como nuestra capa de referencia. Instalaremos el complemento OpenLayers desde Complementos ‣ Administrar e Instalar Complementos ‣ Obtener más. Vea Uso de Complementos para más información sobre el uso de complementos en QGIS.

../_images/1160.png
  1. Una vez instalado, vaya a Complementos ‣ complemento OpenLayers ‣ Agregar capa OpenStreetMap. Esto agregará una capa de teselas pre-representadas creado a partir de datos OpenStreetMap.

../_images/2133.png
  1. Ahora tiene la capa OpenStreetMap cargada en QGIS. Note el Sistema de Referencia de Coordenada (SRC) para esta capa. Está definido como EPSG 3857 Pseudo Mercator. Es importante notar esto, ya que las coordenadas que captemos de esta capa estarán en este SRC.

../_images/375.png
  1. Ahora la tarea es localizar la vecindad general del área que intentamos georeferenciar. Puedes usar herramientas de Desplazar mapa y Acercar zoom para ubicar esa área en la capa OpenStreetMap. Pero podemos aprovechar esta oportunidad para demonstrar otra herramienta que podría ayudarte en el futuro. Sabemos que la imagen que descargamos es para el parque Washington Square en Nueva York. Si buscar ese lugar, podrás ubicar la página wikipedia del mismo. Las coordenadas para el parque están listadas allí.

../_images/443.png
  1. Notarás que las coordenadas están en Grados/Minutos/Secundos y son Latitud y Longitud. Pero debido a que nuestra capa está en proyección Mercator, necesitaremos coordenadas Mercator para ubicar el parque. Aquí es donde la herramienta de línea de comando llamada cs2cs viene bien. Si has instalado QGIS desde el instalador OSGeo4W, ya la tendrás instalada en tu sistema. En Linux y Mac también, viene preinstaladas con QGIS. Abre una ventana de termina y escribe cs2cs para verificar si está disponible. Los usuarios de Windows pueden encontrar una terminar en Start ‣ OSGeo4W ‣ MSYS.

Nota

los usuarios MacOS pueden ubicar la herramienta cs2cs en /Library/Frameworks/PROJ.Framework/Programs/cs2cs

../_images/537.png
  1. Una vez que has verificado que la herramienta cs2cs existe en tu sistema, es tiempo de convertir las coordenadas Latitud y Longitud a Mercator. La manera que esta herramienta funciona es que necesitas especificar un SRC fuente y destino. La definición de SRC podría ser una cadena PROJ4 o un código EPSG. Debido a que ya sabemos el código EPSG para nuestro SRC de entrada y salida, los usaremos. La manera más simple de usar la herramienta es proveer las coordenadas de entrada en la misma línea de comando. Note que la herramienta acepta coordenadas en el orden X Y, por lo que tenemos que ingresar Longitud Latitud. Ingrese el siguiente comando en la terminar y presione Intro (Enter). Note que necesitamos asociar las comillas dobles (») con una barra hacia atrás (\). Una vez presione intro, verá la herramienta procesar las coordenadas e imprimirá la salida de coordenadas X Y en SRC EPSG 3857.

echo "-73d59'51\" 40d43'51\"" | cs2cs +init=EPSG:4326 +to +init=EPSG:3857

-8237364.02 4972720.34 0.00

Si el comando de arriba no le funciona, por favor intente la siguiente variante

echo -73d59'51 40d43'51 | cs2cs +init=EPSG:4326 +to +init=EPSG:3857
../_images/636.png
  1. Copie estas coordenadas y pase a QGIS. Abajo en la ventana QGIS, verá una caja de texto llamada Coordenadas. Ingrese las coordenadas allí en forma X,Y. Presione Intro. Verá el mapa desplazarse un poco, pero no acercarse. Para acercarse al área, seleccione la escala 1:2500 de la lista desplegable Escala junto a la caja Coordenadas y presione Intro.

../_images/736.png
  1. Aquí está! ahora puedes ver el área del Parque Washington Square en tu pantalla. Ahora es tiempo para comenzar la georeferenciación. Abre el Georeferenciador de Raster ‣ Georeferenciador ‣ Georeferenciador. Si no los ves el elemento del menú, necesitará habilitarlo el complemento Georeferenciador GDAL de Complementos ‣ Administrar e Instalar Complementos ‣ Instalado.

../_images/835.png
  1. En la ventana Georeferenciador, vaya a Archivo ‣ Abrir raster. Navegue al archivo JPG descargado y clic Abrir.

../_images/936.png
  1. En el Selector de Sistema de Referencia de Coordenada, elige EPSG:3857 Pseudo Mercator

../_images/1043.png
  1. Ahora clic sobre el botón Agregar Punto en la barra de herramientas y selecciona una ubicación fácil de identificar en la imagen. Esquinas, intersecciones, postes, ect. resultan en buenos puntos de control.

../_images/1161.png
  1. Una vez que haces clic sobre la imagen en una ubicación de punto de control, verás una ventana emergente que pedirá ingresar coordenadas de mapa. Clic el botón Desde la pantalla del mapa.

../_images/1246.png
  1. Encuentre la misma ubicación en su capa de referencia, i.e. la capa OpenStreetMap y clic allí. Las coordenadas son auto-pobladas de su clic en la pantalla del mapa. Clic aceptar. De forma similar, escoja por lo menos 4 puntos en la imagen y agregue coordenadas de la capa de referencia.

../_images/1344.png
  1. Ahora ve a Configuración ‣ Configuración de Transformación.

../_images/1441.png
  1. Elija las configuraciones como se muestra abajo. Asegúrese que el botón Cargar en QGIS cuando esté listo esté marcado. Clic Aceptar. De vuelta en la venttana Georeferenciador, vaya a Archivo ‣ Iniciar georeferenciación. Esto iniciará el proceso de ajustado de la imagen usando los PCTs y creando el raster destino.

../_images/1537.png
  1. Una vez que termina el proceso, verás la capa georeferenciada cargada en QGIS. Si todo fue bien, la verás sobrepuesta sobre la capa OpenStreetMap.

../_images/1635.png
  1. Para hacer que nuestra salida se vea mejor, eliminemos los valores sin datos negro y blanco. Clic-derecho sobre la capa imagen y elige Propiedades.

../_images/1734.png
  1. Pase a la pestaña Transparencia. Queremos indicar que cualquier píxel negro o blanco en la imagen es un valor sin datos y debería hacerse transparente. Ingrese 0 como el Valor sin datos. También, en las Opciones personalizadas de transparencia, clic el botón + y agregue 255 como los píxeles transparentes para cada banda e ingrese 100 como el :Por ciento transparente. Clic Aceptar.

../_images/1830.png
  1. Ahora verás tu imagen georeferenciada sobrepuesta sobre la capa base.

../_images/1923.png

If you want to give feedback or share your experience with this tutorial, please comment below. (requires GitHub account)