00_cubiertas 38 v1 - Asociación Española de Teledetección

40 C. Mattar, et al. / Revista de Teledetección (2012) 38, 36-50
Tabla 1. Datos utilizados: Landsat, Modis y cobertura global
Datos
Fecha
Path-
Row/Swath
Resolución
Espacial
Resolución
Temporal
Descripción
Banda µm
Reflectancia de superficie
1
2
3
4
5
6
7
0.450-0.515
0.525-0.605
0.630-0.690
0.750-0.900
1.550-1.750
10.400-12.500
2.090-2.350
Landsat-7 ETM+
SLC-Off
03-ene-05
24-feb-12
230-96 30 m 16 días
Modis Daily
Reflectance
MOD09GQ
03-ene-05
24-feb-12
250 m Diaria
Reflectancia de superficie
1 0.620-0.670
2 0.841-0.876
GLOBCOVER 2009 Dic-2010 Global 300 m Mapa Cobertura de la superficie
Metodología
La metodología llevada a cabo para estimar
el cálculo del área quemada del PNTP se desarrolló
en tres etapas. La primera etapa corresponde
a la corrección de las imágenes
ETM+ SLC-off con el objetivo de estimar el
área quemada evitando los efectos de bandeado.
La segunda etapa tiene relación con el cálculo
de índices de área quemada y la estimación
de las coberturas que fueron más
afectadas por el incendio. Finalmente, la tercera
etapa tiene relación con la comparación
en la estimación del área quemada a partir de
imágenes ETM+ y los productos de reflectancia
provenientes del sensor MODIS. En las
secciones siguientes se detallan cada una de
estas etapas.
Calibración Radiométrica ETM+
y corrección de bandas
Las imágenes ETM+ correspondientes al
3 de enero del 2005 y al 24 de febrero del
2012 fueron calibradas según el procedimiento
establecido en Chander et al. (2009).
Sin embargo, uno de los problemas más significativos
es el bandeado presente en las
imágenes efecto del SLC-off, que genera líneas
sin datos en la imagen. En Chuvieco et
al. (2012), se presenta una metodología para
determinar el área afectada por un incendio,
contemplando los vacíos que presentan estas
imágenes. Sin embargo, varios autores señalan
la posibilidad de abordar la corrección de
este efecto a través de métodos de interpolación
(Zhang et al., 2007; Pringle et al., 2009).
En este trabajo, cada escena ETM+ fue interpolada
a partir de triangulaciones de Delaunay
relacionadas con estructuras generadas
a partir del diagrama de Voronoi. Este
método aparece bien descrito en De Berg et
al. (2008) y permite la interpolación a una
superficie que no presenta datos a partir de
una triangulación basada en los datos vecinos.
Una muestra del efecto que presenta la
interpolación sobre la imagen se aprecia en
la Figura 2. Sin embargo, los valores resultantes
de esta interpolación no necesariamente
corresponderán a valores con significancia
física. Este método de interpolación
se aplicó a todas las bandas espectrales de la
escena ETM+ adquirida para las dos fechas
analizadas en este trabajo. Una vez corregidos
los efectos del SLC-off en todas las bandas,
se corrigió el efecto atmosférico con el
fin de obtener la reflectancia a nivel de la superficie
utilizando para ello el método propuesto
por Chavez (1996).