El informe “Cálculo de variables forestales para Asturias a partir de la primera y segunda cobertura LiDAR del Plan Nacional de Ortofotografía Aérea” (noviembre de 2025) presenta un inventario regional continuo de estructura forestal a 10 m de resolución, diseñado para apoyar la planificación y la gestión en un contexto de cambio climático. El documento, en cuya realización ha participado miembros de la CuCC y del Instituto de Recursos Naturales y Ordenación del Territorio (INDUROT) de la Universidad de Oviedo, parte de una premisa clara: conocer la altura, la estratificación y la continuidad de la vegetación es clave para entender productividad y resiliencia del monte, estimar biomasa y carbono, y anticipar comportamiento del fuego (continuidad vertical y horizontal del combustible), todo ello con una cartografía homogénea a escala de todo el Principado.
El objetivo técnico es estimar y cartografiar variables dasométricas y estructurales para toda Asturias en dos momentos: 2012 (1.ª cobertura PNOA-LiDAR en Asturias) y 2020 (2.ª cobertura). Sin embargo, el informe subraya que el núcleo del análisis y la interpretación se apoya en 2012, mientras que 2020 se incorpora principalmente como evaluación crítica de calidad: se generan productos, pero se evita concluir sobre cambios temporales porque se detectan problemas de consistencia en la segunda cobertura que comprometen la comparabilidad espacial y temporal.
Como base de datos, el trabajo se enmarca en el programa PNOA-LiDAR, un recurso público, abierto e interoperable. El informe resume su evolución técnica: la primera cobertura se concibe con una densidad nominal homogénea de ~0,5 puntos/m² y, en Asturias, la adquisición corresponde a 2012; la segunda cobertura (2015–2021) es más heterogénea y para Asturias se indica una densidad nominal de ~1 punto/m², con mejoras en precisión. Además, se apunta a la tercera cobertura con densidad homogénea superior y mejores especificaciones, como oportunidad para reforzar futuras comparaciones temporales.
Metodológicamente, el flujo se organiza en dos etapas. La primera es la construcción de una máscara forestal robusta. Este paso es crucial porque, si solo se aplican umbrales de altura sobre LiDAR, pueden aparecer falsos positivos (edificaciones, infraestructuras, roquedos o zonas extractivas) que “imitan” vegetación elevada. Para evitarlo, la máscara combina información de altura LiDAR con series temporales Sentinel-2 procesadas en Google Earth Engine, aplicando control radiométrico/atmosférico y calculando índices espectrales (por ejemplo NDVI, NDWI y NBR2) junto con un criterio explícito de persistencia temporal, de modo que la delimitación forestal sea más estable y fiable.
La segunda etapa calcula variables estructurales mediante métodos de masa en una malla de 10×10 m. Para garantizar consistencia altimétrica, las alturas LiDAR se normalizan con respecto al Modelo Digital del Terreno (MDT) oficial del IGN/CNIG. El procesamiento se automatiza mediante scripts (Python) y se apoya en herramientas del software FUSION (v4.40), destacando GridMetrics (métricas de distribución de alturas) y Cover (fracción de cabida cubierta). El análisis se realiza por estratos verticales definidos por altura: matorral (0,5–2 m), arbolado (2–50 m) y conjunto forestal (0,5–50 m); y se derivan métricas como percentiles/alturas dominantes (p. ej., P90 o P95), altura mínima y máxima, FCC y coeficiente de variación (CV) para capturar densidad y heterogeneidad estructural.
En la cobertura 2012, los resultados muestran patrones diferenciados por estrato. En matorral, la altura dominante alta (P95) presenta una distribución concentrada mayoritariamente entre 0,8 y 1,6 m (intervalo P05–P95), con máximos de frecuencia en torno a 1,7–1,8 m, lo que se interpreta como formaciones arbustivas relativamente desarrolladas y con dosel continuo, relevantes por su papel como combustible denso a baja altura.
En arbolado, la altura dominante (P90) refleja una clara heterogeneidad regional y la presencia de rodales altos puntuales (dominantes >25 m, pero con baja frecuencia). Para la continuidad vertical —un aspecto directamente conectado con propagación del fuego hacia copas— destaca el comportamiento de la altura mínima del arbolado, cuya distribución en 2012 se concentra en valores bajos (máximo alrededor de 2–3 m, con P05–P95 entre 2 y 7 m), sugiriendo copas con bases relativamente bajas en una fracción importante del territorio. En continuidad horizontal, la FCC del arbolado aparece fuertemente sesgada hacia valores bajos, con una mediana ~20%, lo que sugiere que, a escala de píxel, gran parte del paisaje presenta arbolado disperso, aunque coexistiendo con manchas de mayor densidad.
La evaluación de 2020 identifica un problema clave: se observan infraestimaciones y discontinuidades espaciales (cambios bruscos entre áreas adyacentes) atribuibles a la heterogeneidad de sensores y condiciones de adquisición (fechas, calibración, geometría/patrón de escaneo, densidad efectiva y/o preprocesado). A pesar de los controles de calidad y ensayos de armonización, el informe concluye que las incongruencias persisten, por lo que recomienda cautela y evita usar 2020 para análisis de cambio frente a 2012.
Finalmente, el documento explica la publicación y acceso a los ráster en el portal de datos de la Cátedra Universitaria de Cambio Climático (CuCC) (cucc-uodata.es), con organización por estrato/variable/año y un visor que permite explorar, filtrar y consultar valores por píxel. En conjunto, el informe concluye que, con control de calidad y un flujo reproducible, el LiDAR público del PNOA permite generar inventarios forestales regionales de alta utilidad para gestión y análisis territorial.
Puntos clave
Inventario regional a 10 m con variables estructurales forestales derivadas de PNOA-LiDAR.
Separación por estratos (matorral, arbolado y total) para caracterizar verticalidad y continuidad.
Máscara forestal reforzada con Sentinel-2 + Google Earth Engine para reducir falsos positivos.
Resultados 2012: matorral con P95 mayoritario entre 0,8–1,6 m y máximos en 1,7–1,8 m.
Resultados 2012: altura mínima del arbolado concentrada en valores bajos (≈2–3 m), relevante para continuidad vertical del combustible.
FCC del arbolado sesgada a valores bajos, con mediana ~20% (arbolado disperso a escala de píxel).
Cobertura 2020: discontinuidades e infraestimación; no recomendable para comparar cambios con 2012.
Referencia del informe:
De Riba Serrano, P., Menéndez Rodríguez, M., Salgado Fernández, L., Iglesias Fernández, J. C., Cabo Gómez, C., & Colina Vuelta, A. (2025). Cálculo de variables forestales para Asturias a partir de la primera y segunda cobertura LiDAR del Plan Nacional de Ortofotografía Aérea. Universidad de Oviedo; Principado de Asturias.
Visor: cucc-uodata.es.
Fecha Edición
Edita:
