Algebra de mapas y las operaciones con imágenes de satélites

El objetivo principal de este artículo es la utilidad de la Calculadora Raster en la realización de operaciones con imágenes provenientes de sensores remotos, sin embargo es necesario ver someramente (con la salvedad que no soy un experto) algunos conceptos básicos para abordar este tema. La persona que desee profundizar, al final le dejo alguna bibliografía donde podrá obtener mas información. La información provenientes de los sensores remotos es aplicables a disciplinas como Arqueología, geografía, geología, ecología, forestal, hidrología, agricultura, áreas urbanas, inventario de recursos naturales, etc.

 

1. Teledetección:

Es el arte o ciencia de obtener información sobre un objeto, zona o fenómeno, a través del análisis de información obtenida por un dispositivo o sensor que no se encuentra en contacto con el objeto, zona o fenómeno de investigación.

 

2. Sensor remoto:

Es el dispositivo que permite capturar información de los objetos sin tener un contacto directo con ellos. Puede captar información para diferentes regiones del espectro y cada una de estas regiones se denomina canal o banda. Los diferentes canales se pueden caracterizar en función de variables como:

• Amplitud espectral: Región del espectro para la cual capta datos. Se refiere a la anchura espectral de las bandas. 
• Resolución radiométrica: Número de intervalos de intensidad que puede captar.
• Resolución espacial: Tamaño de pixel.
• Resolución temporal: Se refiere al tiempo que tarda el satélite en pasar dos veces por el mismo sitito.

En la tabla siguiente se muestra la resolución espacial y temporal así como la superficie del terreno cubierta por una imagen de algunos satélites:

 

Satélite	Operador	Pan/multi	Resolución espacial (metros)	Ancho escena (km)	Resolución temporal (días)
Spot1/2/3	CNES/SPOT	Pancromático Multiespectral	10 20	60 60	1-4 1-4
Spot 4	CNES/ SPOT	Pancromático Multiespectral	10 20 1000	60 60 2200	1-4 1-4 1
Landsat 5	Space Imaging	Multiespectral	30 80	185	16
Landsat 7	U:S Govemment	Pancromático y Multiespectral	15 30	185	16
IRS IC/D	ISRO-India	Pancromático y Multiespectral	5.8 23 188	70 150 810	5 24 3-5
Radarsat	Canadá	N/A	8- 100	50-500	3- 35
ERS ½	ESA	N/A	30-50	100-500	3- 35
Jers	Japón	N/A	18 15	75-30	44 4-5
Ikonos	Space Imaging	Pancromático y Multiespectral	1 4	11 11	3.5-5 3.5-5
Quick Biird	EarthWatch	Pancromático y Multiespetral	0.82 3.28	22	1.5-4
SPIN-2	Rusia	Pancromático	10 2	200-180	8 8
OrbView	Orbimage	Pancromático y Multiespectral	1 4	8 8	3 3

 

3. Imágenes de satélite

Una imagen raster puede tener una o más bandas de información (entiéndase Fotografía Aérea con tres bandas una por cada color visible RGB;  imágenes de Satélite multiespectrales o hiperespectrales, etc.). Las bandas espectrales, que son más frecuentemente empleadas con la tecnología actual  son las siguientes:

• Espectro Visible (0,4 a 0,7 um). Se denomina así por tratarse de la única radiación electromagnética que puede percibir el ojo humano. Hay tres bandas elementales que son: Azul (0.4um-0.5um); verde (0.5um-0.6um) y rojo (0.6um-0.7um).
• Infrarrojo próximo (0,7 a 1,3 um). Se le denomina también infrarrojo reflejado o fotográfico, puesto que parte de él puede detectarse a partir de filmes dotados de emulsiones especiales. Es importante por su capacidad para discriminar masas vegetales y concentraciones de humedad.
• Infrarrojo medio (1,3 a 8 um). En donde se entremezclan los procesos de reflexión de la luz solar y de emisión de la superficie terrestre. Resulta idóneo para estimar contenido de humedad en la vegetación y detección de focos de alta temperatura. 
• Infrarrojo lejano o térmico (8 a 14 um). Que incluye la porción emisiva del espectro, se detecta el calor proveniente de la mayor parte de las cubiertas terrestres. 
• Microondas (a partir de 1 um). Con gran interés por ser un tipo de energía bastante transparente a la cubierta nubosa.

Las características espectrales de los sensores LANDSAT-TM y MSS se muestran en la siguiente tabla.

MSS		Aplicación	TM		Aplicación
Bandas	l  (m)		Bandas	l  (m)
4	0,5-0,6	Penetración en el agua, turbidez, nieve, sedimentación en la misma, infraestructura urbana y cuerpos de agua.	1	0,45-052	Mapeo de aguas. Mapeo de tipos de bosque, diferenciación entre vegetación y suelos, identificación de rasgos culturales (caminos, urbanización, etc.)
5	0,6- 07	Estudios urbanos, infraestructura caminera y cuerpos de agua.	2	052-0,60	Diferenciación entre tipos de plantas, determinación del estado de la vegetación, identificación de rasgos culturales.
6	0,7- 08	Vegetación, redes de drenaje.	3	0,63-0,69	Diferenciación de plantas a través de mapeo de      clorofila y rasgos culturales.
7	0,8-1,1	Estudios de vegetación, suelos, humedad, contacto entre tierra y agua.	4	0,76-0,90	Determinación de tipos y estado sanitario de vegetación, delimitación de cuerpos de agua.
			5	1,55-1,75	Distinción entre nubes y nieve, determinación del contenido de humedad del suelo y la vegetación.
			6	10,40-12,50	Mapeo térmico.
			7	2,08-2,35	Mapeo hidrotérmico. Humedad retenida en las plantas.

Nota: El sensor LANDSAT de última generación es el ETM+

 

4. Análisis visual de imágenes de satélite (composición color)

La técnica más utilizada  para formar una composición en falso color, consiste en mezclar tres imágenes del mismo sitio y fecha correspondientes a distintas bandas (tabla anterior), representando a cada una de ellas con alguno de los colores rojo, verde o azul.

 

Los 256 niveles digitales de cada banda, representables como una imagen monocromática, se combinan para formar otra imagen en colores. Estos colores y tonos se emplean para la interpretación visual de las imágenes, combinando las bandas de manera que se de un primer acercamiento al contenido de la imagen. La combinación de colores permite discriminar aspectos geológicos, de vegetación, uso del suelo y morfología de zonas urbanas.

 

5. Cálculo de NDVI

 

El cálculo de índices de vegetación es una técnica de uso habitual en teledetección y es comúnmente utilizada para mejorar la discriminación entre dos cubiertas que presenten un comportamiento reflectivo muy distinto en dos o más bandas, por ejemplo para realzar suelos y vegetación en el visible e infrarrojo cercano, y para reducir el efecto del relieve (pendiente y orientación) en la caracterización espectral de distintas cubiertas (Chuvieco, 1996).

 

Los índices de vegetación son útiles porque permiten entre otras cosas.

• Identificación de áreas forestadas-deforestadas.
• Evaluación del estado de la vegetación y su grado de estrés.
• Separación entre distintos tipos de masas vegetales.
• Monitoreo de plagas.
• Evaluación de riesgos de incendio.

De igual forma se pueden derivar variables como contenido de agua en las hojas, productividad neta de la vegetación, contenido de clorofila en la hoja, dinámica fenológica, evapotranspiración potencial, etc. Un aspecto interesante del NDVI es que varía dentro de márgenes conocidos (-1 a +1), cuanto mayor sea el resultado obtenido, tanto mayor será el vigor vegetal presente en la zona observada.

 

El índice se vegetación NDVI se calcula con las Bandas 3 y 4 de una imagen y para ello se utiliza la siguiente ecuación:

 

 

Donde

B4 = Banda 4 de la imagen

B3 = Banda 3 de la imagen

 

6. Transformación Taselled Cap

Supone pasar de las 6 bandas de lansat a 3 bandas con claro sentido físico que representan el brillo de la imagen (albedo), el verdor y la humedad. Estas transformaciones permiten llevar a cabo análisis de la densidad y estado de salud de la vegetación.

 

El brillo refleja los cambios en la reflectividad total de la escena, se trata de una suma ponderara de la reflectividad de todas las bandas excepto el térmico.

 

El verdor. Permite por tanto llevar a cabo análisis de la densidad y estado de salud de la vegetación.

 

La humedad se relaciona con el contenido de agua de la vegetación y el suelo. También se ha encontrado una relación entre este componente y la madurez de la vegetación. El componente de humedad ha sido también correlacionado con la madurez de las formaciones vegetales.

 

Las ecuaciones para calcular estas transformaciones son las siguientes:

 

Brillo = C_kb1*B1 + C_kb2*B2 + C_kb3*B3 + C_kb4*B4 + C_kb5*B5 + C_kb7*B7

 

Verdor = C_kv1*B1 + C_kv2*B2 + C_kv3*B3 + C_kv4*B4 + C_kv5*B5 + C_kv7*B7

 

Humedad = C_kh1*B1 + C_kh2*B2 + C_kh3*B3 + C_kh4*B4 + C_kh5*B5 + C_kh7*B7

 

Donde:

B1… B7 corresponden a las bandas utilizadas

Los coeficientes C_kb, C_kv y C_kh son los que aparecen en la siguiente tabla

Componente B1 B2 B3 B4 B5 B7 Brillo 0.3037 0.2793 0.4743 0.5585 0.5082 0.1863 Verdor -0.2848 -0.2435 -0.5436 0.7243 0.0840 -0.1800 Humedad 0.1509 0.1973 0.3279 0.3406 -0.712 -0.4572

En ArcGIS para realizar composiciones de color se utiliza la función Composite Bands, a la cual de accede de la siguiente manera: Del ArcTool Box, damos clic en  Data Manamegent Tools, seguido de Raster y luego en Raster Processing; donde finalmente seleccionamos la opción Composite Bands y agregamos todas bandas correspondientes.  

 

Para calcular NDVI y la transformación Taselled Cap, simplemente se aplica el algebra de mapas. Dado que el objeto de este artículo es la utilización de la calculadora raster a continuación se verán dos ejemplos al respecto.

 

En la siguiente figura se ilustran 6 bandas de la imagen utilizada, la región de Urabá siempre está bastante nublada pero a manera del ejemplo que se quiere realizar será de mucha utilidad.

 

 

 

7. Calculo del NDVI

Damos clic en Spatial Analyst seguido de Raster Calculator…

 

En la ventana que nos aparece agregamos las capas correspondientes de acuerdo a la siguiente ecuación:

 

NDVI = Float([B4] - [B3]) / Float([B4] + [B3]), el operador Float, en la ecuación, se utiliza para que el resultado de la operación sea punto flotante. Si se omite Float el resultado final va ser entero (cero y uno), lo cual es incorrecto para calcular este índice de vegetación.

 

 

 

Damos clic en el botón Evaluate y el resultado obtenido es el siguiente.

 

 Después de cambiar el color, el resultado es el siguiente:

 

En esta figura se muestran en color rojo las fuentes de agua y el verde más oscuro corresponde la vegetación más saludable.

 

8. Caculo de transformación Taselled Cap

 

En la calculadora raster escribimos las siguientes ecuaciones

 

Brillo = 0.3037 * [B1] + 0.2793 * [B2] + 0.4743 * [B3] + 0.5585 * [B4] + 0.5082 * [B5] + 0.1863 * [B7]

 

El resultado es el siguiente:

Verdor = -0.2848 * [B1] -0.2435 * [B2] -0.5436* [B3] + 0.7243 * [B4] + 0.0840 * [B5] -0.1800 * [B7]

El resultado es el siguiente:

Humedad = 0.1509 * [B1] + 0.1973 * [B2] + 0.3279 * [B3] + 0.3406 * [B4] -0.712 * [B5] -0.4572 * [B7]

 

El resultado es el siguiente:

 

Bibliografía

 

http://www.educ.ar/educar/site/educar/como-leer-e-interpretar-imagenes-satelitarias.html

 

http://www.um.es/geograf/sig/teledet/

 

http://glcfapp.glcf.umd.edu:8080/esdi/index.jsp.

 

http://www.caece.edu.ar/tea/TEA.htm