- Visión General
- Definición de Regiones
- Mosaico Anuales
- Recolección de Muestras
- Clasificación
- Filtros Espaciales
- Máscara con clases de uso y cobertura del suelo
- Evaluación de la Clasificación
Las rutinas de procesamiento de imágenes y clasificación utilizadas para cartografiar las áreas quemadas en el territorio paraguayo siguieron siete etapas:
- Definición de las regiones de clasificación: las regiones se definieron según las ecorregiones del país, comparando la información de múltiples fuentes de clasificación geoespacial y la bibliografía.
- Construcción de mosaicos anuales de calidad Landsat: se generaron mosaicos anuales de alta calidad a partir de imágenes Landsat con los valores mínimos del NBR.
- Recolección de muestras de entrenamiento: se muestrearon las zonas quemadas y no quemadas para generar una biblioteca espectral.
- Mapeo de áreas quemadas (modelo de predicción DNN): el modelo DNN se entrenó utilizando las muestras recolectadas y los mosaicos anuales.
- Filtros espaciales: se aplicaron filtros espaciales para mejorar la precisión y reducir el ruido en los resultados de la clasificación.
- Máscara: se aplicaron máscaras de cobertura del suelo y uso del suelo para mejorar los resultados de la clasificación.
- Evaluación de resultados: los resultados de la clasificación se validaron con datos de referencia y se inspeccionaron visualmente para garantizar una delimitación precisa de las cicatrices de incendios.

Se definieron ecorregiones para la recolección de muestras de entrenamiento y la clasificación de áreas quemadas, con el objetivo de obtener una clasificación más precisa, basada en factores edafoclimáticos y en la vegetación regional. Se definieron de la siguiente manera:

La clasificación se realizó utilizando mosaicos de reflectancia de superficie (SR) de la Colección Landsat 2 (Nivel 1) del USGS (30 m × 30 m), construidos para cada año desde 1999 hasta 2024. Se evaluaron todas las escenas disponibles de los satélites Landsat 5 (de 1985 a 2011), Landsat 7 (de 1999 a 2021), Landsat 8 (de 2013 a 2024) y Landsat 9 (de 2022 a 2024), con un intervalo de retorno de 16 días.
Creamos una biblioteca espectral basada en la delimitación manual de polígonos de áreas quemadas y no quemadas para utilizarla como muestras de entrenamiento. Estas muestras se estratificaron por sensores Landsat (recogidas en diferentes años) y por cada ecorregión. La recolección de muestras de entrenamiento se realizó en las siete ecorregiones de clasificación, garantizando la representación de las distintas características espectrales presentes en cada una de ellas.
De esta forma, se obtuvo un conjunto de muestras por sensor y para las siete ecorregiones, para ser utilizadas en el entrenamiento del modelo de clasificación.
El modelo de clasificación utilizado fue la Red Neural Profunda (Deep Neural Network, DNN), que consiste en modelos computacionales basados en cálculos matemáticos capaces de realizar aprendizaje automático y reconocimiento de patrones visuales.
El algoritmo de mapeo de las áreas quemadas constó de dos fases: entrenamiento y predicción. Basándose en las muestras de entrenamiento de áreas quemadas y no quemadas, se utilizaron las siguientes bandas espectrales como entrada para el modelo de clasificación de áreas quemadas: rojo (ROJO – 0,65 µm), infrarrojo cercano (NIR – 0,86 µm) e infrarrojo de onda corta (SWIR 1 – 1,6 µm y SWIR 2 – 2,2 µm). Estas bandas espectrales de Landsat se eligieron por su sensibilidad a los eventos de incendio.
Los datos de entrenamiento se dividieron en dos conjuntos: el 70 % de las muestras se utilizaron para el entrenamiento y el 30 % para la prueba.
Se aplicó un filtro espacial para eliminar los píxeles aislados calculando el número de píxeles conectados y eliminando aquellos con cuatro o menos conexiones (~0,36 hectáreas). Este proceso también rellenó los huecos en o alrededor de las cicatrices cuando había más de dieciséis píxeles (~1,44 hectáreas) clasificados conjuntamente como quemados.
Además, se llevó a cabo un procesamiento posterior a la clasificación para recuperar la información de la fecha del píxel quemado a partir del mosaico anual construido a partir del NBR mínimo, identificando el mes en el que se cartografió la cicatriz del incendio.
Se aplicaron máscaras posteriores a la clasificación para reducir los errores de comisión. Para las ecorregiones del Chaco Seco, Chaco Húmedo, BAAPA y Cerrado, se empleó la clase de agua de la Colección 2 de Cobertura y Uso del Suelo de MapBiomas. Sin embargo, en el caso de la ecorregión del Pantanal, se seleccionó la clase de agua de la Colección 1, ya que sus límites espaciales mostraban una mejor idoneidad. En todas las ecorregiones se aplicó la regla de exclusión (26) correspondiente a la clase de agua permanente.
Las evaluaciones de la clasificación de cicatrices de incendios se realizaron con imágenes Landsat, mediante inspección visual, análisis estadísticos y relación con datos de cobertura y uso del suelo de MapBiomas. Además, se compararon con mapas de referencia, incluyendo FIRMS (1 km), GABAM (30 m), MODIS MCD64A1 (500 m), FIRE CCI (250 m) y focos de calor del INPE (1 km).
