En las áreas de clima mediterráneo como la nuestra, las elevadas temperaturas durante el verano, unido a la sequedad de la vegetación y del ambiente, suponen un grave riesgo de incendio, con el consiguiente riesgo de pérdida de hábitats naturales.

Para identificar las áreas de vegetación natural afectadas por un incendio podemos recurrir a la información espectral proporcionada por las imágenes de satélite, ya que ésta nos permitirá identificar las áreas que hayan sufrido un cambio considerable en su respuesta espectral.

incendio_qgis

En las siguientes imágenes del satélite Landsat 8 podemos ver una misma zona afectada por un incendio en color natural (combinando las bandas del visible – rojo, verde y azul) y en falso color infrarrojo (combinando bandas del infrarrojo):

incendio_qgis_1

incendio_qgis_2

En el siguiente enlace podemos hacer un barrido entre las dos imágenes para observar con más detalle las diferencias:

http://earthobservatory.nasa.gov/NaturalHazards/view.php?id=88445

En el gráfico, a continuación, podemos ver como la vegetación sometida a algún tipo de estrés, frente a la vegetación sana, presenta un aumento de la reflectancia en la longitud de onda del rojo, debido a la menor presencia de pigmentos que absorban esa radiación, y una disminución de la reflectancia en el rango del infrarrojo cercano, debido a la menor presencia de estructuras celulares que produzcan dispersión de la radiación en ese rango.

incendio_qgis_3

Para estudiar la relación entre los valores de reflectancia de la vegetación en el rojo y el infrarrojo cercano, el índice más ampliamente utilizado es el NDVI, que tiene la forma:

incendio_qgis_4

Para realizar este ejercicio de detección de zonas afectadas por incendio en QGIS comenzamos cargando dos imágenes de satélite, una anterior y otra posterior a un incendio. En las siguientes composiciones en color natural, se observa que, aparte de la rotación de cultivos, en la zona de bosque se ha producido un cambio (círculo rojo):

incendio_qgis_5

A continuación calculamos el índice de vegetación (NDVI) de ambas imágenes introduciendo la expresión correspondiente en la “Calculadora ráster” de QGIS.

incendio_qgis_6

incendio_qgis_7

Si analizamos los valores que toman los índices de ambas fechas en la zona de bosque, podemos usar esos valores del NDVI para generar dos máscaras de superficie forestal. Para ello introducimos una expresión condicional en la “Calculadora ráster” que contenga ese rango de valores de NDVI que permiten separar el bosque sano del resto de superficies.

incendio_qgis_8

Superponiendo ambas máscaras podremos identificar las zonas que han sufrido cambios drásticos:

incendio_qgis_9

Finalmente, podríamos realizar una operación de resta con las máscaras para extraer la superficie afectada por incendio y calcular su superficie.

 

Formación de calidad impartida por profesionales

1 estrella2 estrellas3 estrellas4 estrellas5 estrellas (2 votos, promedio: 5,00 de 5)
Cargando...