Proyecto Peridotita 2: Tutorial: Extraer la serie temporal de NBR con Google Earth Engine
Introducción
El análisis de la regeneración vegetal en áreas afectadas por grandes incendios forestales requiere datos satelitales precisos, continuos y cuantitativos. Dentro del marco del proyecto escolar Climate Detectives de la Agencia Espacial Europea (ESA), nuestro equipo en el IES Monterroso asumió el reto de estudiar la evolución de la masa forestal en Sierra Bermeja desde el devastador incendio de 2021 hasta la actualidad.
Para evaluar la severidad del fuego y el rebrote posterior, el Índice de Quemado Normalizado ($NBR$, por sus siglas en inglés) es el indicador biofísico por excelencia. Sin embargo, el procesamiento práctico de estos datos nos enfrentó a un obstáculo técnico inesperado: las limitaciones de exportación dinámica de la interfaz web oficial de Copernicus Browser. En esta entrada explicamos cómo superamos esta barrera metodológica migrando nuestro flujo de trabajo a la computación en la nube con Google Earth Engine, permitiéndonos extraer una serie temporal limpia de cinco años con la que medir el impacto del cambio climático y la resiliencia de nuestro entorno local.
Resumen de los problemas de Copernicus Browser (para el blog)
- Bloqueo del motor de estadísticas: Copernicus Browser inhabilita automáticamente el botón Statistical Info chart para cualquier capa generada en la pestaña «Personalizar», impidiendo la extracción de series temporales de scripts propios como el NBR.
- Fallos de sincronización y caché de cuentas: El visor principal no actualiza en tiempo real los perfiles creados en el Configuration Utility, bloqueando la selección de capas nativas personalizadas a pesar de haber iniciado sesión correctamente.
- Obsolescencia de alternativas web: Tras la absorción definitiva de EO Browser en febrero de 2025, la plataforma web oficial carece de un método directo para exportar gráficas dinámicas de índices personalizados a CSV, obligando a migrar a entornos analíticos como Google Earth Engine.
Requisitos previos
- Una cuenta gratuita en Google Earth Engine Code Editor.
- El área de estudio identificada (Sierra Bermeja).
1.1. Delimitar el Área de Interés (AOI):Herramienta de mapa en la interfaz.
Entra en el Code Editor de Earth Engine, busca Sierra Bermeja en el mapa inferior y utiliza la herramienta de dibujar polígono (icono de figura geométrica en la esquina superior izquierda del mapa) para trazar la zona afectada por el incendio.
Para verificar este paso, comprueba que en la parte superior del editor de código aparece automáticamente la línea var geometry: Polygon....
2.2. Pegar y configurar el script de NBR:Editor de código (panel central).
Borra cualquier texto del panel central de código y pega el siguiente script ajustado para procesar imágenes Sentinel-2 L2A entre 2021 y 2026 con menos del 20% de nubosidad:
JavaScript
// 1. Filtrar la colección de Sentinel-2 L2A por área, fechas y nubes
var s2 = ee.ImageCollection('COPERNICUS/S2_SR_HARMONIZED')
.filterBounds(geometry)
.filterDate('2021-09-01', '2026-09-01')
.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20));
// 2. Calcular el índice NBR: (B8 - B12) / (B8 + B12)
var nbrCollection = s2.map(function(img) {
var nbr = img.normalizedDifference(['B8', 'B12']).rename('NBR');
return img.addBands(nbr);
});
// 3. Renderizar la gráfica temporal del valor medio de NBR
print(
ui.Chart.image.series(nbrCollection.select('NBR'), geometry, ee.Reducer.mean(), 10)
.setOptions({
title: 'Evolución del NBR en Sierra Bermeja (2021-2026)',
hAxis: {title: 'Fecha'},
vAxis: {title: 'Índice NBR'}
})
);
Para verificar este paso, comprueba que el editor no resalta ninguna línea con subrayado o icono de error de sintaxis.
3.3. Procesar la serie temporal:Panel derecho ‘Console’.
Haz clic en el botón Run (Ejecutar) situado en la barra superior del editor de código.
Para verificar este paso, comprueba que en la pestaña “Console” (a la derecha) se genera una gráfica interactiva con la curva del NBR desde septiembre de 2021.
4.4. Exportar los datos a CSV:Ventana emergente del gráfico.
Pasa el cursor sobre la gráfica de la consola, haz clic en el icono de la ventana emergente (esquina superior derecha del gráfico) y selecciona la opción de descargar como CSV.
Para verificar este paso, abre el archivo descargado en tu ordenador y confirma que contiene las filas con la serie temporal completa de fechas y valores promedio.
CSV de datos de NBR
ODS de datos de NBR
Conclusión
La transición desde los visores web convencionales hacia entornos de programación espacial como Google Earth Engine marca un hito fundamental en nuestra investigación dentro de las iniciativas MonteSteam y AIDARAC. Lo que comenzó como un bloqueo en la descarga de datos se transformó en una oportunidad para elevar la rigurosidad científica de nuestro proyecto: ahora disponemos de un conjunto de datos temporal homogenizado (2021-2026) que documenta con precisión métrica la respuesta del ecosistema de Sierra Bermeja.
Este proceso demuestra que la teledetección no consiste solo en observar imágenes bonitas desde el espacio, sino en saber extraer, procesar y validar los datos numéricos que hay detrás de cada píxel. Con la serie temporal de NBR finalmente en nuestras manos, el alumnado del IES Monterroso pasa a la siguiente fase del proyecto: correlacionar la regeneración de la cubierta vegetal con las variables meteorológicas locales para proponer acciones fundamentadas en favor de la conservación del medio ambiente.
