Incendios forestales x especies

Este mapa cruza donde hubo incendios con que arboles habia ahi para saber cuantas hectareas de cada especie se quemaron en Espana en el verano de 2025.

Que es una hectarea (ha)?
Una hectarea = 10.000 m2. Equivale aproximadamente a un campo de futbol y medio, o una manzana de ciudad. Es la unidad estandar para medir superficies de terreno forestal y agricola.

Las 5 capas del mapa

1. EFFIS — Perimetros de incendios de la UE
Poligonos que dibujan el contorno exacto de cada incendio detectado por satelite. Fuente: European Forest Fire Information System (Comision Europea / JRC). Cada poligono tiene: fecha inicio/fin, superficie en hectareas, tipo de vegetacion afectada, si esta en zona Natura 2000.
maps/EFFIS/a93770120b43423fa517fb5264afd20a.shp | CRS: EPSG:3035 | 593 poligonos

2. NASA FIRMS — Puntos de calor MODIS
Detecciones de fuego activo del satelite MODIS (Collection 6.1). Cada punto es una deteccion termica. Se convierte en un circulo de 500m de radio y se disuelven los solapados. Atributos: brillo (Kelvin), confianza (%), potencia radiativa del fuego (MW).
maps/NASA/fire_archive_M-C61_717179.shp | CRS: EPSG:4326 | 8.319 puntos → 1.391 poligonos

3. Quemado combinado (EFFIS + NASA)
Union geometrica de ambas fuentes. Se reproyectan a EPSG:25830 (ETRS89 / UTM zona 30N), se fusionan con unary_union y se separan en poligonos individuales. Esta es la vision mas completa de toda la superficie quemada.
data/output/layer_burned_combined.geojson | 1.104 poligonos | 525.532 ha

4. MFE — Mapa Forestal de Espana (MITECO)
El Mapa Forestal de Maxima Actualidad del Ministerio para la Transicion Ecologica. Cubre todo el territorio nacional parcela a parcela. Cada poligono dice: Especie1/2/3 (arbol dominante, secundario, terciario), O1/O2/O3 (% de ocupacion de cada especie), tipo de bosque (coniferas/frondosas/mixto), formacion arborea, cobertura arborea, modelo de combustible. Descargado de las 17 CCAA por separado en formato Shapefile (EPSG:25830).
maps/MITECO/mfe_andalucia/ a maps/MITECO/mfe_regiondemurcia/
911.360 poligonos forestales en total (~11 GB de datos brutos).
© Ministerio para la Transicion Ecologica y el Reto Demografico — Banco de Datos de la Naturaleza

5. Interseccion — EL RESULTADO
El cruce espacial final. Donde la zona quemada combinada se solapa con el mapa forestal = hectareas quemadas por especie.
data/output/layer_intersection.geojson | 11.417 fragmentos | 166.838 ha con especie identificada

Metodologia completa

Paso 1: Carga EFFIS. Se lee el shapefile EFFIS (EPSG:3035). Se filtra por pais=Espana y fecha entre 15-Jun y 21-Sep 2025. Se reproyecta a EPSG:25830 (ETRS89/UTM30N).

Paso 2: Carga NASA FIRMS. Se lee el shapefile de puntos MODIS (EPSG:4326). Se filtra por confianza >= 50% y mismas fechas. Se reproyecta a EPSG:25830. Cada punto se convierte en un circulo de 500m (buffer(500)). Los circulos solapados se disuelven con unary_union en poligonos unicos.

Paso 3: Union de quemado. Se concatenan EFFIS + NASA buffereado. Se disuelve la geometria completa con unary_union y se explota en poligonos individuales. Resultado: 1.104 poligonos de zona quemada, 525.532 ha.

Paso 4: Carga MFE. Se leen los 17 shapefiles del MFE (uno por CCAA). Se usa bbox espacial para cargar SOLO los poligonos que caen dentro del bounding box de las zonas quemadas (con margen de 10km). Se filtran solo parcelas con Especie1 no vacia. 911.360 poligonos forestales.

Paso 5: Interseccion espacial. Se ejecuta gpd.overlay(quemado, mfe, how='intersection'). Esto corta geometricamente la zona quemada con el mapa forestal. Cada fragmento resultante tiene la especie del MFE y la superficie real quemada de esa especie. Se recalcula area_ha con la geometria cortada.

Paso 6: Calculo por especie. Para cada fragmento de interseccion se detectan las columnas Especie1/O1, Especie2/O2, Especie3/O3. La superficie de cada especie se pondera por su % de ocupacion: ha_especie1 = area_ha × (O1/100). Se agrupa por nombre de especie y se suma. Resultado: tabla con total_ha y num_fragments por especie.

De donde salen las columnas Especie1/O1, Especie2/O2, Especie3/O3?
Son atributos nativos de los shapefiles del MFE (Mapa Forestal de Espana) publicados por MITECO. No los genera esta aplicacion; vienen dentro de los archivos originales descargados del Banco de Datos de la Naturaleza.

Especie1, Especie2, Especie3 = nombre cientifico del arbol dominante, secundario y terciario de cada parcela forestal. Son identificados por los tecnicos de MITECO mediante foto-interpretacion de imagenes aereas + trabajo de campo.
O1, O2, O3 = porcentaje de ocupacion de cada especie en esa parcela (0-100%). Indican que proporcion de la parcela cubre cada arbol.

Ejemplo concreto: Si una parcela quemada mide 10 ha y tiene Especie1=Pinus halepensis con O1=60 y Especie2=Quercus ilex con O2=30, se calcula:
ha Pinus halepensis = 10 × 60/100 = 6 ha
ha Quercus ilex = 10 × 30/100 = 3 ha
De esta forma se asigna la superficie quemada proporcionalmente a cada especie segun su presencia real en el terreno.

Paso 7: Export. Se exportan 5 capas GeoJSON (simplificadas a 50m de tolerancia para el visor web) + CSV + JSON de resumen. El MFE se recorta espacialmente (sjoin) a solo poligonos que intersecten con zonas quemadas: de 911k a 11k poligonos (de 11GB a 13MB).

EFFIS (EPSG:3035, 593 poligonos) ──→ reproyectar a 25830 ──┐ ├─ UNION ──→ 1.104 poligonos NASA FIRMS (EPSG:4326, 8319 pts) ──→ buffer 500m ─────────┘ │ │ MFE (EPSG:25830, 911k polig.) ──→ filtro bbox ──→ filtro Especie1 ──────────┤ │ gpd.overlay(how='intersection') ←─────────┘ │ area_ha × O1/100 por especie │ 166.838 ha en 40+ especies

Datos del Excel del diccionario (MaximaActualidadMFE25_50)

Cada carpeta MFE contiene un Excel con 20 pestanas de metadatos. Lo mas relevante:

ESPECIES — 193 codigos de especies forestales con nombre cientifico y comun.
ModeloCombustibles — 13 modelos de combustible Rothermel. Predicen como se propaga el fuego segun el tipo de vegetacion. Van del Modelo 1 (hierba fina seca, propagacion rapida) al 13 (restos gruesos de tala).
FORARB — 73 formaciones arboreas (encinares, hayedos, pinares de pino albar, etc.).
Publicacion — Fechas de publicacion del MFE50 y MFE25 por CCAA. Aragon fue el ultimo MFE25 completado (2023).

Por que los mapas de MITECO pesan tanto?

El MFE cubre todo el territorio nacional parcela a parcela — ~500.000 poligonos con mucho detalle por CCAA (escala 1:25.000). Cada poligono almacena 50+ atributos. Los archivos EFFIS y NASA solo cubren las zonas donde hubo fuego, que son una fraccion minuscula del territorio. Por eso EFFIS pesa ~3MB y cada MFE pesa ~500MB-1.9GB.

En el visor web mostramos solo la parte del MFE que coincide con incendios: de ~11GB a ~13MB.

Controles

Flechas del teclado — Mover el mapa suavemente
🎥 Volar — Añade waypoints (draggable) y la camara vuela entre ellos. Para grabar movimiento de camara sobre el mapa.
🔥 — Particulas de fuego animadas sobre las zonas quemadas (toggle on/off)
⛶ — Pantalla completa
⚙ — Cambiar colores de las capas
◐ — Alternar tema oscuro / claro
☰ — Mostrar/ocultar paneles laterales
Especies — Click = seleccionar/deseleccionar (multi-seleccion), hover = ver foto del arbol (Wikipedia)

Tipos de mapa

IGN Espana — Mapa topografico del Instituto Geografico Nacional, nombres en castellano
OpenStreetMap — Mapa colaborativo, nombres mayormente en castellano
Satelite (Esri) — Imagen de satelite real
Terreno — Mapa topografico con relieve (OpenTopoMap)
CartoDB oscuro/claro — Mapas minimalistas de CARTO (con nombres)
Claro/Oscuro sin nombres — Mapas limpios sin etiquetas, ideal para visualizar datos sin distracciones

Fuentes de datos

EFFIS: effis.jrc.ec.europa.eu — European Forest Fire Information System, Joint Research Centre, Comision Europea
NASA FIRMS: firms.modaps.eosdis.nasa.gov — Fire Information for Resource Management System, NASA Goddard Space Flight Center
MFE: miteco.gob.es — Mapa Forestal de Espana, Banco de Datos de la Naturaleza, MITECO
© Ministerio para la Transicion Ecologica y el Reto Demografico

Tecnologia

Backend: Python + FastAPI + GeoPandas + Shapely + Fiona + PyProj
Frontend: Leaflet.js + Canvas API (particulas de fuego)
Analisis espacial: gpd.overlay (interseccion), unary_union (disolucion), sjoin (join espacial), buffer (NASA pts)
CRS: ETRS89 / UTM zona 30N (EPSG:25830) como sistema de referencia comun
Deploy: Railway + GitHub

🔥 Incendios forestales x especies
Elias Leon Siminiani, Carlos Marlaska | Artefacto

Ajustes

Colores de capas