NDVI en ArcGis parte 1

Ya he hablado anteriormente de como calcular el NDVI normalizado en ArcGis, sin embargo lo hice en artículo en el que se hablaba de las operaciones las calculadora raster y creo que muchos lo han pasado por alto. En este artículo lo volvemos a retomar agregando algunas posibilidades que tenemos cuando se calcula este índice.

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.

Además nos permite delimitar zonas ocupadas por agua, permitiendo una aproximación a la digitalización de ríos y lagunas. El NDVI varía de (-1 a +1), los valores de NDVI son mucho mayores para la vegetación que para el suelo descubierto, las nubes y el agua. Cuanto mayor sea el resultado obtenido, tanto mayor será el vigor vegetal presente en la zona observada; aquellas áreas ocupadas por agua dado que no tienen vegetación su valor de NDVI es negativo.

El índice se vegetación NDVI en ArcGis 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

En este caso utilizaremos una imagen de la región de los departamentos de Sucre y Bolívar (Colombia) donde se muestra el sistemas de lagunas del río Magdalena y también quiero destacar el golfo Morrosquillo donde se localiza un importante acuífero aluvial de la zona.

La composición de las bandas se ve como se muestra en la siguiente figura.

Paso 1. 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.

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

Aca pondemos ver que todo los que es agua es de color rojo y es el NDVI mas bajo, lo cual nos permite diferenciar claramente las fuentes de agua.

Haciendo un zoom en el área indicada vemos lo siguiente.

Observamos el cauce principal del río y aquellas otras zonas susceptibles de ser inundadas y deposito de sedimentos del río.

Paso 3. Si queremos cartografiar estas zonas, podemos reclasificar el NDVI calculado, extraer aquel dato relacionado con los ríos y luego convertirlo a shape para el proceso de digitalización.

Para ello nos vamos a Spatial analyst seguido de Reclassify.

Nos aparece una ventana donde incluimos la siguiente información.

• Input raster: seleccionamos NDVI.
• Reclass fiel: Seleccionamos Value.
• Ouput raster: Damos un nombre al archivo de salida.

Damos clic en el botón Classify. Nos aparece una ventana donde incluimos la siguiente información.

• En Method: seleccionamos Estándar Desviation.
• En Interval Size: seleccionamos 1/3 Std DEV.

Damos clic en Ok y luego en OK. Obtenemos el siguiente resultado.

La figura anterior nos muestras varias cosas interesantes:

Los valores del NDVI_R (NDVI Reclasificado) menores a 5 corresponden a diferentes cuerpos de agua, la variación se debe posiblemente a sedimentación, profundidad, vegetación, etc.

Los valores mayores, corresponden a vegetación en diferentes estados de salubridad, densidad, tamaño, etc.

Miremos en detalle lo que nos muestra el recuadro No 1.

Al hacer zoom sobre los resultados vemos con mayor detalle los elementos del paisaje, como ven en la figura original no se veía lo etiquetado como el canal del Dique (que a propósito por efectos del fenómeno del niño en 2010 se desbordó e inundó muchos pueblos del sur del departamento del Atlántico).

Si hacemos clic sobre cada uno de estos elementos característicos veremos el valor del NDVI_R correspondiente, por ejemplo para el río Magdalena (ríos en general) le corresponde un valor de 5. Vemos que el NDVI (reclasificado) puede tomar valores de hasta 17, cuanto mayor es este valor mas saludable es la vegetación.

En la siguiente parte de este artículo veremos como podemos utilizar los resultados de esta reclasificación para digitalizar los cuerpos de agua. Los espero entonces en NDVI en ArcGis parte 2.