La teledetección es una herramienta imprescindible para trabajar en áreas remotas y de difícil acceso como es el caso de la mayoría de los volcanes tropicales, en especial si son activos. Dentro de la metodología del grupo de investigación, una de las mayores dificultades es diferenciar el hielo de la nieve en las imágenes de satélite, con el fin de analizar la evolución reciente de los glaciares y obtener una cartografía fiable, precisa y objetiva. Con esta finalidad, se aplicará un análisis multitemporal a las imágenes Landsat y Aster, con la ayuda de los softwares Erdas 9.2 y Arcgis 9.2 y de los criterios propuestos en la bibliografía (Dozier, 1981; Dozier, 1989; Klein, et al., 1998).

En primer lugar se preparará la imagen en el software Erdas 9.2, especializado en tratamiento de imágenes de satélite, donde se establecerá un criterio de clasificación uniforme en función de la reflectancia de cada pixel. Después, en Arcgis 9.2, se agruparán por un lado los píxeles que tienen un valor numérico que equivalga a la nieve, y por otro, los que se corresponden con el hielo glaciar. Se utilizarán imágenes correspondientes con la estación seca (noviembre-abril), momento en que el manto nival se encuentra en su mínimo anual. En concreto, el tratamiento será es el siguiente:

a) Diferenciación entre la banda 2 y la banda 5 (TM2 – TM5) de la imagen.

b) Suma de la banda 2 y la banda 5 (TM2 + TM5) de la imagen.

C) Cociente entre la diferencia de las bandas 2 y 5 y la suma de las bandas 2 y 5:

(TM2 – TM5) / (TM2 +TM5).

El resultado es una nueva capa raster donde se diferenciará la superficie glaciar y nival con un alto grado de precisión. Finalmente, la superficie cubierta por el hielo se tratará en Arcgis 10, calculándose el área de cada celda con el mismo valor numérico (igual reflectancia), y multiplicado por el número de píxeles que tengan ese valor. A partir de la superficie obtenida, será posible estimar el volumen medio de los glaciares de cada área de estudio mediante la siguiente fórmula (Klein, et al., 1998):

V = C.A^b 

donde,

V: volumen,

A: superficie glaciar

C = 0,048 (valor empírico)

b = 1,36 (valor empírico)

 

Referencias.-

Dozier, J. (1981) A method for satellite identification of surface temperature fields of subpixel resolution. Remote sensing of environment, 11: 221-229.

Dozier, J. (1989) Spectral signature of alpine snow cover from the Landsat Thematic Mapper. Remote sensing of environment, 28: 9-22.

Klein, A. G y Isacks, B. (1998) Alpine glacial geomorphological studies in the central Andes using Landsat thematic mapper images. Glacial Geology and Geomorphology; rp01/1998.