Alaska Satellite Facility - Distributed Active Archive Center

Mapear la Inundación Regional con SAR de banda L por el Espacio – NEEDS TESTING

Animation of inundation in Amazon rain forest. Credit: Chapman, Guritz 2016; RTC: ASF 2015; Includes Material © JAXA/METI 2007.

Source: Bruce Chapman, JPL, and Rick Guritz, ASF. Traducción a español y adaptación cruda para Mac OS X, QGIS y GIMP por Eric Anderson

Advanced (OSX): Use QGIS and GIMP to classify and map regional inundation

Introducción

Son muy conocidas las habilidades del Radar de Apertura Sintética (SAR) de adquirir imágenes de la superficie de la Tierra en todo tipo de clima, de penetrar la cobertura nubosa, y en las condiciones de luz baja. Como resultado de la disponibilidad reciente de modelos de terreno de alta resolución a los 30 metros de la Misión Topográfica de Radar volada en el Transbordador Espacial (SRTM), y un interés en hacer los datos de SAR más fácil para usar, la Facilidad de Satélite de Alaska (ASF) ahora ofrece productos de imágenes corregidos de terreno radiométrico (RTC) de la banda L de los datos de ALOS PALSAR. En esta receta de datos, describimos los métodos para (a) crear una animación de 18 productos PALSAR RTC y (b) producir un mapa de inundación que cuantifica el número de eventos de inundación con el conjunto de datos.

Prerrequisitos

Materiales para generar una animación de inundación:

  • Productos de ALOS-1 PALSAR RTC
    • Para este ejercicio, descargue de la lista de gránulos a continuación
  • QGIS con Orfeo Toolbox (Los usuarios de Windows ven la Sección Windows)
  • GIMP

Materiales para generar un mapa con color de la inundación:

  • QGIS 2.18
  • imágenes filtradas que se generó previamente (R_nnnnn_sfx)

Links a los productos de muestra de RTC

Haga clic en los nombres de los productos para descargar cada granulado individualmente.

Parte 1: Pasos de procesamiento para generar una animación de la inundación

Descargar y organizar los datos

    1. Descargar los datos a través de los URLs de los datos de muestra en la Prerrequisitos.
    2. Cree un nuevo directorio para guardar sus archivos (ej HH).
    3. Extraer los datos de los archivos .zip, y copiar la imagen HH de cada gránulo a otra carpeta separada de HH.

Convertir las imágenes a decibelios

Convertir las imágenes RTC de unidades lineales a db

    1. En QGIS, cargar todas las imágenes GeoTIFF de RTC HH de la carpeta.
    2. Abrir la herramienta Raster Calculator (Menu de Raster > Raster Calculator).
    3. Crear un directorio HH_db para guardar los resultados, y nombrar cada resultado como Rnnnnn_db.tif donde nnnnn es el número de órbita.
    4. Calcular 10 * log10 (“RTC_HH_granule”).
    5. Repetir el proceso para todos los gránulos de RTC.
Raster Calculator para un gránulo

Procesar y clasificar las imágenes

Filtrar las imágenes para reducir el moteado

    1. En el Processing Toolbox, busca “despeckle.”
    2. Hacer doble clic sobre Despeckle (lee).
    3. Elegir la imagen RTC de insumo (ej R04877_db).
    4. Crear un nuevo directorio HH_db_sf para guardar los resultados.
Herramienta Despeckle (lee) en la caja de herramientas de procesamiento
    1. Definir la ubicación y nombre del archivo de salida en Output Image como Rnnnnn_db_sf.tif, donde nnnnn es el número de órbita.
    2. Dejar las otras opciones como están por defecto, y correr la herramienta.
    3. Repetir los pasos 2–6 para todas las imágenes.

Excluir los píxeles erróneos a las orillas de las imágenes.

    1. El proceso de Despeckle deja una franja de píxeles extras con valores muy negativos a las orillas de las imágenes, que hay que quitar antes de seguir el análisis.
      1. Abrir Raster Calculator.
      2. Calcular (“Rnnnnn_db_sf@1” > -9999 ) * “Rnnnnn_db_sf@1”.
      3. Guardar la imagen resultante como Rnnnnn_db_sfx.
      4. Repetir los pasos 2-3 para todas las capas.
    2. Quitar todas las capas del proyecto para que solo queden las 18 filtradas y limpiadas (Rnnnnn_db_sfx).

Desplegar las imágenes

    1. Hacer clic derecho sobre una capa y elegir Properties / Style.
    2. Cambiar el tipo de representación (Render) a Singleband pseudocolor.
    3. Ver al histograma (Histogram) para interpretar la distribución de los valores de los píxeles en la imagen. Se nota dos máximos que indican una distribución bimodal.
    4. Regresar a la pestaña de Style. Cambiar la Interpolación a Discrete.
    5. A la par de Classify, hacer clic dos veces sobre el ‘+‘ para tener dos clases. Cambiar el valor de la primera clase a -14, y la segunda clase a inf. Hacer clic sobre Apply.

Notar: Los valores que se determina para distinguir entre el agua la y tierra pueden variarse mucho entre diferentes imágenes, y por lo general, se requiere varias iteraciones para determinar los umbrales adecuados.

Histograma que muestra la distribución bimodal

Refinar el umbral hasta que los canales de los ríos se contienen el agua, y hay muy pocos puntos azules sobre la tierra afuera de los canales.

Notar: Determinamos que el valor “-11” minimiza mejor el moteado de agua sobre la tierra.

Refinar más los umbrales para determinar una segunda clase con áreas más brillantes como resultado del doble-rebote de vegetación inundada. En esta etapa, vamos a determinar otro umbral al otro lado del máximo global del histograma.

    1. Hacer clic en el ‘+‘ a la par de Classify, para añadir otra clase. Definir un valor más alto como 0, y cambiar el color que representa los valores más alto que 0 a amarillo y Aplicar los colores.
    2. Revisar la imagen para el moteado amarillo.
    3. Disminuir este segundo umbral uno por uno hasta que el moteado amarillo es similar a lo que se observó para la clase azul.

Notas: Determinamos que el valor “-4” es un buen valor para la segunda clase. 

Para el proceso de determinar los umbrales adecuados, es una buena práctica considerar una imagen que representa la inundación más severa y otra imagen que representa la inundación menos severa. El conocimiento de la cobertura geográfica de las inundaciones estacionales ayuda confirmar que los resultados son correctos.

Clasificación de la imagen con umbrales determinados

Aplicar las mismas propiedades de los colores y las clases a todas las imágenes.

    1. En el Layers Panel, hacer clic derecho sobre la imagen con la clasificación deseada. Seleccionar Styles / Copy Style.
    2. Seleccionar las otras 17 imágenes en el Layers Panel (usando Ctrl o Shift). Hacer clic derecho, y Seleccionar Paste Style.
Todas las imágenes con la misma simbología.

Extraer una imagen procesada para cada imagen original.

    1. Después de aplicar los umbrales y los colores correctos a todas las imágenes, apagar todas las capas con la excepción de la primera (R04887_db_sfx). Acomodar bien el zoom y la ubicación de la imagen dentro de la ventana para preparar para la exportación a formato png.
    2. Selecionnar Project / Save As Image…
    3. Guardar la imagen como class_nnnnn_db.png.
    4. Hacer clic en OK.
    5. Apagar la primera imagen y seleccionar la siguiente imagen.
    6. Repetir los pasos 2-5 para todas las imágenes.

Generar y visualizar un GIF animado

Generar un GIF animado para mostrar la secuencia de imágenes clasificadas que indican la cobertura de la inundación.

    1. Abrir Gimp.
    2. Seleccionar File / Open As Layers…
    3. Navegar al directorio dónde se guardó las imágenes clasificadas (class_nnnnn_db.tif). Seleccionar todos los archivos (18).
    4. Seleccionar File / Export As…
    5. Guardar la animación como anim.gif.
    6. Marcar As animation, y “Delay between frames where unspecified 500 milliseconds.”

Navegar al archivo anim.gif en su computadora, y visualizar la animación en un browser.

Ejemplo de la Animación

Hacer clic en el siguiente imagen para ver el gif resultante:

Lamina de la animación de la inundación en el bosque de las Amazonas. Crédito: Chapman, Guritz 2016; RTC: ASF 2015; Includes Material © JAXA/METI 2007.

Parte 2: Pasos para generar un mapa con colores de la inundación

Clasificar las imágenes

Clasificar cada imagen como agua (1) o tierra (0).

    1. En QGIS, cargar cada imagen filtrada (R_nnnnn_sfx). Abrir la herramienta Reclassify values (simple). Seleccionar la primera imagen. Para, Replace condicion elegir [1] Low value < grid value < high value. Abrir el Lookup Table, y asignar los siguientes valores. Guardar el resultado como reclass_nnnnn.tif.
    2. Repetir para todas las imágenes (o usar “Run as batch process”).
    3. Quitar las imágenes originales, y cargar las imágenes clasificadas si ya no están en el proyecto.
Tabla para reclasificar los valores.

Crear el producto final

Añadir todas las imágenes.

    1. Abrir el Raster Calculator.
    2. Hacer doble clic sobre la primera imagen, hacer clic sobre +, y.
    3. Seguir añadiendo cada imagen y el símbolo + hasta que todas las imágenes estén en la ecuación.
    4. Especificar el nombre (count.tif) y ubicación del resultado.
    5. Hacer clic sobre OK para ejecutar el cálculo.

Notar: La imagen combinada saldrá con un rango de grises entre blanco y negro, o con colores asignados al azar. Es bueno asignar los valores 0 y 1 al color negro, porque los valores de 1 normalmente son un resultado del moteado que se quedó después del filtro (y 0 indica que nunca se clasificó el píxel como inundado). Se puede cambiar los colores con las propiedades de la imagen. Elegir los colores que son distinguibles y que tengan lógico en términos de la imagen combinada (mayores valores indican inundación más frecuente).

Ejemplos de un mapa con colores de las inundaciones

Mapas de inundación de la selva amazónica. La leyenda indica la frecuencia de inundación durante los 18 gránulos. Los colores con el número 18 representa el río constante. Crédito: Anderson 2017; RTC: ASF 2015; Includes Material © JAXA/METI 2007.

Cómo acceder a Orfeo Toolbox (Windows)

Para utilizar la herramienta despeckle debe tener acceso a Orfeo Toolbox, que no se incluye con la versión independiente de Windows de QGIS 2.18.12. Para agregar la caja de herramientas a usted debe descargar la caja de herramientas usando el instalador de OSGeo4W.

Notar: Para que estos pasos funcionen correctamente, debe tener QGIS preinstalado.

  1. Navegue hasta QGIS.org.
  2. Seleccione “Download Now”.
  3. Bajo “Advanced Users” descargar OSGeo4W (64 bit).
  4. Abra la OSGeo4W setup.
  5. En la ventana de diálogo, seleccione “Advanced Install” y haga clic en “Next>”.
  6. Seleccione “Install from Internet” y haga clic en Next>.
  7. Seleccione Root Install Directory (por ejemplo, C: \ OSGeo4W) y haga clic en “Next>”.
Configuración de OSGeo4W
  1. Seleccione Local Package Directory (por defecto funciona correctamente) y haga clic en “Next”.
  2. Seleccione “Direct Connection” y haga clic en “Next>”.
  3. Haga clic en “http://download.osgeo.org” para seleccionar y haga clic en “Next>”.
  4. En el Select Packages: menú, escriba “otb” en el cuadro de búsqueda.
  5. Haga clic en “Libs” para expandir.
  6. En “Libs” haga clic en cada “Skip” una vez y haga clic en “Next>”.
  7. Deje todos los valores predeterminados y haga clic en “Next>”.
  8. Haga clic en la casilla de verificación y haga clic en “Next>”.
  9. Su descarga comenzará.
  10. Haga clic en “Finish” una vez finalizada la instalación.
  11. Abrir QGIS 2.18.12.
  12. Seleccione “Processing” en la barra de menú principal.
  13. Haga clic en “Options …”
  1. Haga clic en “Providers” para expandir.
  2. Haga clic en ‘GRASS commands’ para expandir.
  3. Haga doble clic en el espacio en blanco a la derecha de “Msys Folder”.
  4. Suprima la ruta existente para “Msys Folder”.
  5. Haga clic en “Orfeo Toolbox” para expandir.
  6. Asegúrese de activar “Activate”.
  7. Haga doble clic en el espacio en blanco a la derecha de “OTB applications folder”.
  8. Introduzca la ruta de la carpeta de las aplicaciones OTB en el cuadro (por ejemplo, “C: \ OSGeo4W \ apps \ orfeotoolbox \ applications”).
  9. Haga doble clic en el espacio en blanco a la derecha de “OTB command line tools folder”.
  10. Introduzca la ruta de la carpeta bin OSGeo4W64 en el cuadro (por ejemplo, “C: \ OSGeo4W64 \ bin”).
  11. Salga de QGIS y reinicie la aplicación.
  12. En el menú Toolbox de la derecha aparecerá ahora el Orfeo Toolbox.
  13. Haga clic aquí para regresar a la Prerrequisitos.
Processing Options: Providers