Como Calcular NDVI en ArcGis

Un índice de vegetación puede ser definido como un parámetro calculado a partir de los valores de la reflectancia a distintas longitudes de onda, y es particularmente sensible a la cubierta vegetal, de acuerdo a Gilabert et al (1997).

El utilizar estos índices tiene su fundamento en el particular comportamiento radiométrico de la vegetación. Una cubierta vegetal en buen estado de salud, tiene una firma espectral que se caracteriza por el contraste entre la banda del rojo (entre 0,6 y 0,7 μm.), la cual es absorbida en gran parte por las hojas, y el infrarrojo cercano (entre 0,7 y 1,1 μm.), que es reflectada en su mayoría. Esta cualidad de la vegetación permite la realización de su valoración cualitativa.

El Índice de Vegetación Diferencial Normalizado se calcula mediante la siguiente expresión:

NDVI= ((IRC-R))/(IRC+R)

dónde: IRC es la reflectividad en el infrarrojo cercano y R es la reflectividad en el rojo.

El rango de valores de las reflexiones espectrales se encuentra entre el 0 y el 1; ya que, tanto la reflectividad del infrarrojo cercano como la del rojo, son cocientes de la radiación reflejada sobre la radiación entrante en cada banda espectral. Por consecuencia de estos rangos de valores, el NDVI varía su valor entre -1 y 1. El NDVI posee un gran valor en términos ecológicos, ya que es un buen estimador de la fracción de la radiación fotosintéticamente activa interceptada por la vegetación (fPAR) según Monteith (1981).

En las imágenes Landsat 8, la banda 4 (0.630 – 0.680 µm) corresponde al rojo (R) y la banda 5 (0.845 – 0.885 µm) al infrarrojo (IRC), por lo tanto para el cálculo NDVI se requiere contar con ambas bandas (Descargar imágenes Landsat 8).

Calcular NDVI en ArcGIS

Para calcular el NDVI simplemente se requiere aplicar su respectiva ecuación, cabe mencionar que dicho proceso se lo puede realizar en diferentes aplicaciones que cuenten con una calculadora de imágenes ráster (ArcGIS, Erdas, QGIS, gvSIG, Surfer, Idrisi, ENVI, etc).

El primer paso luego de abrir la aplicación ArcMap 10 es cargar las bandas 4 y 5 con la herramienta Add Data desde el directorio donde se encuentren almacenadas, seguidamente se muestra una ventana solicitando la creación de pirámides, en este caso dejar los valores por defecto y aceptar, cabe señalar que se requiere una licencia activa de la extensión Spatial Analyst, para ejecutar el Raster Calculator.

Una vez cargadas las bandas, desde la caja ArcToolbox, abrir la siguiente herramienta:

ArcToolbox > Spatial Analyst Tools > Map Algebra > Raster Calculator

Para obtener los valores NDVI en imágenes Landsat 8 usando ArcGIS se aplica la siguiente ecuación:

NDVI = Float (banda 5 – banda 4) / Float (banda 5 + banda 4)

El resultado es una imagen ráster que contiene valores que van desde -1 a 1 (siendo los valores más cercanos a 1 la vegetación más vigorosa). Generalmente la imagen del NDVI se muestra en una escala de grises, para dar un aspecto más agradable y de fácil interpretación, dirigirse a las propiedades del ráster y seleccionar una paleta de colores en la pestaña de simbología (clic derecho > Properties > Symbology)