Teledetección.

Jesús Anaya.

Corrección atmosférica y radiométrica según el modelo mejorado de

“objeto negro”, Chavez (1996).

Objetivo:

 Calcular la radiancia (L) en imágenes Landsat a partir de los niveles digitales que
almacenan las imágenes Landsat cuando se descargan de la Internet.
 Buscar valores de irradiancia (E) correspondientes a cada banda de Landsat.
 Calcular la reflectividad en la parte superior de la atmósfera (TOA) y en superficie
(SR/BOA).
 Obtener las Firmas Espectrales a partir de la reflectividad.

Introducción:

El cálculo de la radiancia es un paso fundamental para llevar las imágenes a una escala
radiométrica común, esto es crucial cuando se utilizan imágenes de distintos sensores o
distintas fechas.

Cada día la cantidad de energía que llega a la superficie terrestre varía. Cambia, entre otros,
la distancia del sol a la tierra, su ángulo de incidencia y las condiciones atmosféricas. Esto
implica que la energía reflejada por los objetos (lagos, nevados, bosques, ciudades, etc.)
aunque permanezcan constantes también varía. Lo cual resulta inconveniente dado que
podemos concluir que un objeto está cambiando, cuando en realidad lo que cambia es la
energía que llega al objeto.

Se utilizará parte del modelo COST y el método propuesto por Chander y Markham (2003)
para llevar las imágenes a una escala radiométrica común, la lectura más reciente y que se
requiere para este ejercicio es el artículo de Chander, et al 2009.

Las unidades que se utilizan en este procedimiento son:

-Irradiancia (E): watts/m2

-Radiancia (L): watts/m2 sr

-Reflectividad -rho- (ρ): sin unidades

-Transmisividad (T): sin unidades

1
Teledetección.
Jesús Anaya.

-Nivel Digital (ND): sin unidades

Las firmas espectrales son fundamentales en teledetección, pues nos permiten identificar
un objeto a partir del cambio en reflectividad en diferentes longitudes de onda.

La Reflectividad -rho- (ρ) (en inglés reflectance) se puede despejar de la ecuación de la

siguiente figura. Ver “measured radiance”. Si se despeja rho, la irradiancia pasa a dividir la
radiancia, lo cual corresponde a la ecuación 2 que se muestra en la sección de “Reflectividad
en TOA”.

2
Teledetección.
Jesús Anaya.

Radiancia
1. Calibración interna, aquí se definen los valores máximos y mínimos de radiancia por
sensor (LMAXλ, LMINλ) y se derivan, Gain (Slope) y Bias (offset). Esto relaciona la
radiancia medida por los sensores (L) con el nivel digital DN almacenado en los
píxeles de las imágenes cuando se descargan de internet, pueden ser de 1 byte (0-
255) o de dos bytes (0-65535).
La radiancia (energía) que llega al sensor es convertida en niveles digitales y se
conoce como “cuantización” (Q), el dato almacenado se conoce como ND or DN
(DN: digital numbers, ND: niveles digitales). Hay una relación directa entre la
radiancia y el ND que se puede establecer con dos parámetros: gain y bias. Estos
valores (gain, bias) existen en los metadatos de cada imagen o en la literatura.
2. Los ND deben ser llevados por el usuario nuevamente a radiancia utilizando la
ecuación que relaciona los ND con la radiancia, esto se hace con una regresión lineal
con su pendiente y su constante, es decir, gain y bias siendo ND la variable
independiente.

Entonces el proceso para estimar la radiancia que mide el sensor a partir del ND es
el siguiente: primero se debe convertir los valores de los píxeles (ND denominados

3
Teledetección.
Jesús Anaya.

QCAL, Q: quantization, cal: calibrated) a radiancia (L), para cada rango de longitud
de onda lambda (λ):

Lλ: radiancia espectral medida en el sensor W/m2*sr (W para vatio y sr para estéreo
radián)
LMINλ: radiancia espectral escalada al valor mínimo del píxel (metadatos)
LMAXλ: radiancia espectral escalada al valor máximo del píxel (metadatos)
QCAL: es el nivel digital ND, es decir cada uno de los pixeles de cada banda de la
imagen.
QCALMAX = 65535 (según el sensor, pues depende de la resolución radiométrica)

Lλ = LMINλ+ QCAL * (LMAXλ- LMINλ) / QCALMAX (eq. 1)

Tomado de: USGS (2018)Landsat7 Data Users Handbook

4
Teledetección.
Jesús Anaya.

3. Luego se debe utilizar la ecuación de reflectividad, donde se tiene en cuenta la

energía recibida por el sensor (radiancia del paso anterior) y la energía proveniente
del sol, conocida como energía incidente, es decir, la irradiancia para ese día.

Este modelo requiere muchos valores, pero en su mayoría son constantes y disponibles en
la literatura o en el metadato de la imagen, solo tres datos de entrada se deben obtener
para este cálculo: la distancia del sol a la tierra, el ángulo de elevación del sol y las bandas
Landsat con sus niveles digitales (ND o QCAL). Nótese que el proceso hay que hacerlo para
cada una de las bandas desde el azul hasta el SWIR. Las bandas termales tienen un proceso
distinto.

Reflectividad en TOA
Una forma de normalizar la energía que llega del sol es calcular la reflectividad (que elimina
el efecto del ángulo solar y la distancia entre el sol y la tierra) en el momento en que se
tomó la imagen. La reflectividad rho (ρ) es la proporción de la energía que llega al sensor
con respecto a la energía incidente.

Reflectividad en TOA (Top Of Atmosphere) sin corrección por efecto de la atmósfera

(eq.2):
𝐿 ∗𝜋∗𝑑
𝜌 =
𝐸𝑆𝑈𝑁 ∗ cos ( ∗ (90 − 𝜃 ))

Donde, Lλ: radiancia medida en el sensor W/m2 (eq. 1)

d2 = distancia del sol a la tierra en unidades astronómicas (Chander et al 2009), que se basa
en el día juliano (ver tabla al final) y en metadatos.

Ꝋs : Elevación del sol (sun-angle elevation, metadatos).

5
Teledetección.
Jesús Anaya.

Figura tomada de google.

ESUNλ : irradiancia W/m2 (energía que proviene del sol, tabla Chander et al, (2009))

El Objeto negro
Una dificultad en este proceso es cuantificar la energía que llega al sensor proveniente de
la atmósfera, y no del objeto de interés. Es decir, parte de la energía que llega al sensor se
debe a la energía reflejada por la atmósfera (haze ó path radiance).

 No toda la radiancia que llega al sensor Lλ es de la superficie terrestre, realmente es

la suma de la radiancia reflejada por la atmósfera Lλ,haze (La notación en la figura es
Lp –path radiance-) más la radiancia reflejada por la superficie terrestre.
 Lλ,haze se calcula como la radiancia que llega al sensor desde un objeto oscuro Lλ,min
menos la radiancia de un objeto oscuro Lλ,1%

Para cuantificar esta magnitud se propone ubicar manualmente en cada escena un “objeto
negro”, por ejemplo un lago, una sombra de relieve, sombra de nube, en la cual podamos
inferir que la energía que recibe el sensor no fue reflejada desde la superficie. Por tanto,

6
Teledetección.
Jesús Anaya.

se concluye que la energía que recibe el sensor asociada a ese objeto negro procede
mayormente de la atmósfera.

Cálculo de la energía (path radiance Lp) Lλ,haze

Se asume que la cantidad de energía que llega al sensor de un objeto negro es solo del 1%
de la radiación incidente, el resto de energía que llega al sensor viene de la atmósfera y por
tanto debe eliminarse. De esta forma la cantidad de energía que proviene de la atmósfera
para cada pixel en cada banda es definida por la siguiente ecuación:

Lλ,min = Lλ,haze + Lλ,1% por tanto, La radiancia que mide el sensor en el objeto oscuro es la
suma de la energía proveniente de la atmósfera más la energía que llega del objeto oscuro.

Despejando

Lλ,haze=Lλ,min - Lλ,1%

Donde,

Lλ,min: Radiancia W/m2*sr medida en el sensor proveniente de un objeto oscuro, basada en

la eq.1.

Lλ,min = LMINλ+ QCAL_OO * (LMAXλ- LMINλ) / QCALMAX

QCAL_OO: nivel digital del pixel del objeto oscuro en la banda λ.

Lλ,1%: Radiancia teórica de un objeto negro en función del ángulo de toma de la imagen θ,
la distancia al sol d y la energía incidente ESUNλ.

Lλ,1%= 0.01 * ESUNλ * cos((π /180)*(90-θ)) / π* d2

Esta última ecuación viene de la ecuación de reflectividad, vista arriba (eq. 2), solo que se
despeja el término de radiancia y la reflectividad ρ se reemplaza por uno por ciento (0.01).

Reflectividad en superficie SR/BOA

Esta es la última ecuación que debe utilizarse para obtener la reflectividad 𝜌 para cada
banda para cada píxel. Se conoce como Reflectividad en superficie (Surface Reflectance),
también reflectividad en BOA (Bottom of Atmosphere). Después de aplicar esta ecuación
cada píxel tendrá valores de reflectividad y por tanto se pueden generar las Firmas
Espectrales.

7
Teledetección.
Jesús Anaya.

(𝐿 − 𝐿 , )∗𝜋∗𝑑
𝜌 =
𝐸𝑆𝑈𝑁 ∗ cos ( ∗ (90 − 𝜃 ))

Para cada píxel de la imagen se puede generar una firma espectral. El resultado debe ser
algo similar al siguiente ejemplo de firma espectral de bosque.

Firma espectral de un bosque seco tropical (parque Los Colorados) basado en la

reflectividad calculada utilizando el método del objeto oscuro (naranja) y el modelo LEDAPS
utilizado por el USGS (azul). Nótese que en esta figura el eje X no está a escala. Compare
las firmas espectrales de esta figura con la primera figura que se presenta en estas notas.

Para generar reflectividad en superficie de forma operativa se recomienda el modelo

LEDAPS:

https://daac.ornl.gov/MODELS/guides/LEDAPS_V2.html

https://polipapers.upv.es/index.php/raet/article/view/3582/6061

Lectura requerida:

Chander, Marckham, Dennis, 2009. Summary of current radiometric calibration coefficients

for Landsat MSS; TM, ETM+, and EO-1 ALI sensors. Remote sensing of environment. 113,
893-903.

8
Teledetección.
Jesús Anaya.

Tabla de día Juliano. Tomado de google, https://docs.oracle.com/.