En la entrada de hoy volvemos a hablaros del índice NBR (Normalized Burn Ratio) sobre el que escribimos en entradas anteriores al comentar el incendio de Moguer-Doñana.

Recordando un poco sobre este índice, el NBR es el resultado de normalizar la diferencia de reflectancias entre una banda del infrarrojo cercano (IRC) y una del infrarrojo medio (SWIR), lo que nos permite identificar los daños sufridos por la vegetación. Sin embargo, lo que realmente queremos tratar en esta entrada es cómo obtener el índice NBR y su correspondiente dNBR con SNAP y qué dificultades podemos encontrar por el camino.

El índice NBR se obtiene de la ecuación: NBR= (ρ_IRC - ρ_SWIR) / (ρ_IRC + ρ_SWIR), donde ρ es la reflectancia. Para nuestras imágenes de Sentinel-2 la banda 8 corresponde al infrarrojo cercano y tiene una resolución espacial de 10 metros, y la banda 12 corresponde al infrarrojo medio, con 20 metros de resolución.

En primer lugar se debe tener descargada la imagen de la zona de interés y abrirse con SNAP. Una vez se ha hecho esto pueden ocurrir dos cosas: que SNAP pregunte si queremos hacer un remuestreo de todas las bandas para tenerlas a la misma resolución espacial, o que directamente abra el producto para empezar a trabajar con él. En el primer caso, dado que tendremos todas las bandas con la misma resolución (en concreto elegiremos 20m), será posible hacer directamente la operación para obtener el NBR con la opción Band Maths. En el segundo caso, antes de calcular el NBR deberemos hacer un remuestreo espacial, esta herramienta la encontramos en Raster → Geometric operations → Resampling. De esta manera, generaremos un nuevo producto con las resoluciones que nosotros elijamos. Es muy importante trabajar con las bandas de misma resolución espacial ya que el Band maths  no admite combinar bandas de distintos tamaños.

Una vez solucionado el problema de la resolución espacial, ya podemos obtener nuestro NBR empleando de nuevo el Band maths, este lo podemos encontrar de dos maneras: podemos pulsar el botón derecho sobre nuestra imagen remuestreada o podemos buscarlo en Raster → Band maths. Lo que vamos a generar ahora es una nueva banda, que llamaremos NBR para diferenciarla del resto. Deberemos insertar la ecuación del índice, que en SNAP quedará escrito de la siguiente manera: (B8-B12)/(B8+B12). Esto lo podemos hacer escribiendo directamente esta expresión sin identificar cada banda, porque el propio programa lo reconoce, o podemos seleccionar cada banda en la mitad izquierda de la ventana, bajo Data Source. Tras esto, ya tendremos generado nuestro nuevo índice, que aparecerá junto con el resto de bandas y que por tanto podremos visualizar simplemente haciendo doble clic.

Por otro lado tenemos el dNBR, la diferencia entre el índice NBR previo y posterior al incendio. Con este índice podemos delimitar mejor el área quemada, ya que las superficies que no se han visto afectadas por el incendio tendrán los mismos valores de NBR antes y después de este. Por tanto, al calcular la diferencia de NBRs todas estas áreas se cancelarán. Cabe mencionar que dada la variabilidad de los valores obtenidos de reflectancia en el agua, la imagen del dNBR resaltará la zona incendiada, así como todos los cuerpos de agua.

Para obtener este índice debemos partir de los índices NBR de las imágenes anterior y posterior al incendio. Después, deberemos incorporar en el Band maths la ecuación: dNBR = NBR_pre-NBR_pos, seleccionando cada banda en cada producto (en este caso trabajamos con dos productos diferentes: el anterior al incendio y el posterior). Como sucedía anteriormente, ya habremos generado nuestro nuevo índice que aparecerá junto con las demás bandas en uno de los dos productos, el que elijamos.

A continuación os mostramos la comparación entre el NBR y el dNBR obtenidos con SNAP de la zona de Moguer-Doñana: a la izquierda el NBR del 1/7/17 y a la derecha el dNBR de la diferencia entre el 21/6/17 y el 1/7/17.

Con esto concluiría este apartado sobre el uso del índice NBR en incendios. Para cualquier duda, no dudéis en contactarnos!