Hemos visto anteriormente cómo corregir las imágenes Landsat para considerar la reflectanciaTOA (tope de la atmósfera). Aquí hay un diagrama que nos ayudará a entender este fenómeno:
La reflectancia que hemos calculado es por lo tanto, el % de luz reflejada con respecto a la totalidad de la luz visible incidente . Pero como se ha visto en el diagrama , el sensor del satélite mide dos cosas al mismo tiempo: la luz reflejada por nuestros objetivos en la superficie de la tierra ,y la luz dispersada por las partículas suspendidas en el aire.
Se puede realizar la corrección atmosférica de las imágenes satelitales para eliminar la porción de luz debida a la difusión.
Como ya hemos dicho, ésto sólo tiene interés si usted trabaja en imágenes separadas en el tiempo El porcentaje de difusión es el mismo para todo el conjunto de una imagen, realizar esta corrección es perder el tiempo.
Las instrucciones para realizar esa corrección en imágenes Landsat está
disponible en
Pero si usted es un principiante seguramente se le dificultará entender el procedimiento. Veremos aquí el procedimiento usando las herramientas de QGis.
Método DOS (Dark Object Subtraction)
El método que vamos a emplear es el más comúnmente empleado.
La resta de objetos negros es una método de corrección atmosférica empirico
simple para la imagen satelital , que asume que la reflectancia de los objetos
oscuros incluye una componente importante de la difusión atmosférica.
El concepto es simple: si se tiene un objeto completamente negro y, por lo tanto , no refleja la luz, este objeto que debería tener un valor 0 de reflectancia tendrá un valor no nulo (diferente de 0) en el sensor del satélite. Este valor será la reflectancia debida a las partículas presentes en la atmósfera.
La resta de objetos oscuros investiga en cada banda el valor de píxel más oscuro. La difusión es suprimida restando este valor a cada píxel de la banda .
Lo primero: la corrección que vamos a hacer solo aplica a la luz visible y el infrarrojo cercano. Vamos a trabajar por lo tanto , para imágenes Landsat8, en las bandas 2 (Azul), 3 (verde), 4 (Rojo) y 5 (infrarrojo cercano).
Lo segundo que hay que saber: el principio en uso es calcular la distribución para la banda 4 (rojo) y usar una fórmula de correspondencia a fin de determinar la corrección para las otras tres bandas.
Lo tercero: vamos a trabajar con las bandas 2 , 3, 4 y 5 en reflectancia TOA y también en radiancia bruta (DN) para la banda 4 (rojo). En el manual del usuario citado previamente no es evidente a primera vista . Veamos el procedimiento paso a paso:
1- se calcula la reflectancia TOA para las cuatro bandas con la fórmula
TOA banda = ( (DN * 0.00002)-0.1) / sen ( elevación solar )
2- analizando los valores brutos de radiancia (DN) de la banda roja, vamos a determinar una radiancia debida a la difusión atmosférica ( DNred )
3- tranformaremos esta difusión en reflectancia TOA
CORrojo = ( ( DNrojo* 0.00002) -0.1) / sen ( elevación solar )
4- con un ábaco se ingresa esta difusión de la banda roja y se obtiene la difusión para las bandas 2,3 y 5 ( CORazul , CORverde , CORir).
5- se restan los valores de reflectancias obtenidos de TOA para obtener los valores de reflectancia de las cuatro bandas corregido por el efecto atmosférico.
Rrojo = TOA ( rojo) – CORRrojo
Razul = TOA ( azul) – CORazul
Rverde = TOA ( verde) – CORverde
Rir = TOA ( ir ) – CORir
Veamos en detalle los pasos 2 y 4.
Cálculo de la difusión de la banda roja.
Hay varios métodos para determinar la corrección a realizar en la banda roja. El método recomendado es el conocido como Bin 5 – 0.008. Usted encontrará aquí una descripción detallada del método.
En la práctica, con QGis, usaremos el histograma de la banda para determinar el valor umbral de la difusión atmosférica. 1- Primero debe ejecutar la herramienta Estadísticas de capa raster en los datos en la banda 4 (rojo).
Una vez ejecutado, tendrá una línea de resultados en el panel del Visualizar resultados
Haga doble clic en esta línea para abrir los resultados en su navegador:
Tenga en cuenta el valor mínimo aquí 3916.
Buscaremos el valor que se asignará a la corrección de la difusión atmosférica en el histograma de la imagen. Haga clic en la capa de la imagen -> Propiedades -> Histogram
Ingrese el valor mínimo que se encuentra en las estadísticas (aquí 3916) en el cuadro Min, luego haga clic en el cuadro Max para mover el valor mínimo de la línea de puntos en el histograma.
Puede usar el formato de línea para encontrar la base del histograma, pero en este ejemplo, se usará (y se recomienda) el formato de barra. Si el histograma está en formato de línea, haga clic en Preferencias / Acciones y desactive la casilla de verificación Marcar como líneas .
Aparecerá un histograma en formato de barra. El histograma de QGIS no representa/ no muestra todos los valores como el histograma de ArcGIS, pero muestra muy bien la forma que puede ser apropiada para establecer un valor de dispersión que represente la base del histograma (en caso de un aumento relativamente abrupto de las frecuencias). Los siguientes pasos muestran cómo se puede usar el histograma de QGIS para establecer la dispersión. Pase el cursor del mouse sobre el histograma para ver una lupa. Arrastre la lupa hasta el extremo inferior del histograma para ampliar la base del mismo. Si hace demasiado zoom o si desea volver a empezar por cualquier razón, haga clic en Preferencias / Acciones y Recalcular histograma . El histograma se amplía a continuación.
Observará que la barra con el valor más bajo sigue siendo bastante superior al valor de imagen bajo que ingresó anteriormente en el área de Min (3619); el histograma de QGIS no muestra muchos valores en una cola estadística (Landsat 8 generalmente tiene largas colas estadísticas). La base del histograma es el área utilizada para establecer la dispersión. Haciendo zoom en las tres primeras barras del histograma:
En este ejemplo, los valores comprendidos, aproximadamente, entre 5570 y 5595 pueden acercarse a la base del histograma; éstos son valores muy bajos y no están demasiado en la cola estadística (este proceso puede ser ligeramente subjetivo). Para las necesidades de este documento, el DN 5569 será utilizado como base para el histograma. El histograma se ampliará una vez más para este ejemplo.
La regla a utilizar es tomar como valor de dispersión el valor mínimo de la primera barra que exceda una frecuencia de 5 (de ahí el nombre Bin 5) y que no tenga frecuencias inferiores a 5 a su derecha. En el siguiente ejemplo:
Las dos primeras barras exceden el valor 5, pero como la tercera es más baja, se retendrá el valor mínimo de la cuarta barra.
En nuestro ejemplo, el valor de radiancia que se retendrá para la corrección atmosférica de la banda roja será 5569.
Podemos pasar a la tercera etapa,
3- transformamos este resplandor en reflectancia TOA.
COR(rojo) = ((DNrojo * 0.00002) -0.1) / sen (elevación solar)
COR(rojo) = ((5569 * 0,00002) -0,1) /0,42631886 = 0.026694
Este valor corresponde al Bin 5. Dado que el método es el Bin5 – 0.008, eliminamos 0.008 de este valor: 0.018694
Es este valor el que debe restarse de los valores de TOA del rojo previamente calculado para tener los valores corregidos de la difusión atmosférica.
Corrección de otras bandas.
A partir del valor de corrección para el rojo obtenemos los valores que se aplicarán a cada una de las otras bandas:
Ir a la siguiente página Encontrará una ventana donde puede ingresar el valor de corrección para el rojo
Ingrese su valor (aquí 0.018647 ), automáticamente tendrá los valores de las otras correcciones
Así obtenemos las correcciones de
- 0.06309 para la banda azul (2)
- 0.03442 para la banda verde (3)
- 0.00626 para la banda de infrarrojo cercano (5)
Para tener las bandas corregidas usamos la calculadora raster restando a cada banda corregida TOA la constante de corrección que hemos encontrado.
