Cómo transformar Números Digitales (DN) en Radianzas (L) para bandas Landsat

La radianza espectral (Lλ) para cada banda Landsat 5 o Landsat 7 puede determinarse mediante la siguiente expresión:

radianza1_1

donde DN es el número digital de cada píxel, LMAX y LMIN son constantes de calibración, QCALMAX y QCALMIN son, respectivamente, el valor más alto y el más bajo del intervalo de radianzas rescaladas en DN. Las unidades para Lλ son W/m2/sr/μm.

Del archivo de la metadata, LT50030692005197COA00_MTL.txt, seleccioné los parámetros para determinar las radianzas de una banda 6 de Landsat 5 a partir de los DN; tal como se destaca en la imagen siguiente:

metadata

Con dichos valores, usados con el comando gdal_calc.py a continuación:

.
gdal_calc.py -A b6_clip.tif --calc "((15.303-1.238)/254)*(A-1)+1.238" --type Float32 --outfile b6.rad2.tif
.

se obtuvo el ráster de radianzas (b6.rad2.tif). Mediante el raster calculator de QGIS tendríamos lo siguiente:

radianza2

Esta entrada fue publicada en GDAL, QGIS, SIG, Software Libre. Guarda el enlace permanente.

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s