SUMMARYIn this study, the stand volume (V, m 3 ha -1 ) of Pinus radiata plantations of the province of Lugo (NW of Spain) was estimated using LiDAR-derived metrics from the low-density countrywide LiDAR dataset provided by the Spanish National Aerial Photography Program. The estimate was based on a oriented field inventory design to improve the fit of the model and on a inferential method based on the variability of the model parameter estimation and the variability of the LiDAR variable in the target population, which allowed to estimate the population mean and the standard error of V. A total of 25 permanent field inventory plots −that represent adequately the existing range of ages, stand densities and sites of the species in the province− were measured in order to fit the empirical relationship between field-measured V and LiDAR-derived metrics, using linear and (multiplicative) power function models. Regression analysis indicated a strong relationship between V and the 70 percentile of the LiDAR normalized height distribution (h 70 ) using a power function model (R 2 = 0.899; RMSE = 63.4 m 3 ha -1 ). The selected model was applied to all timber and high-pole development stands of P. radiata with coverage of this species superior to the 60 % according to the Spanish National Forest Map. Finally, a spatial-explicit map that describes V was produced. In 2009, the estimated mean population for V was 166.3 m 3 ha -1 with a relative sampling error of 13.9 %. This methodology will allow generating thematic maps every six years (temporal resolution scheduled for the PNOA LiDAR flights).Key words: Airborne Laser Scanning (ALS), Forest Inventory (FI), Remote Sensing (RS), estimators, regression models.
RESUMENSe estimó el volumen de madera de rodal (V, m 3 ha -1 ) de Pinus radiata en la provincia de Lugo (NO España) a partir de métricas de datos LiDAR de baja resolución del PNOA (Plan Nacional de Ortofotografía de España). La estimación se basó en un diseño de muestreo orientado para mejorar el ajuste del modelo y en un método inferencial basado en la variabilidad de la estimación de los parámetros del modelo y en la variabilidad de la variable LiDAR en la población objetivo, lo que permitió estimar la media poblacional y el error en la estimación de V. Se usaron 25 parcelas permanentes, que cubren el rango existente de calidad, edad y densidad en los rodales de la especie en la provincia y se ajustaron modelos lineales y potenciales. El mejor resultado se obtuvo con un modelo potencial entre V y el percentil 70 de la altura normalizada (h 70 ) (R 2 = 0,899; REMC = 63,4 m 3 ha -1 ). Dicho modelo se aplicó a todas las teselas de P. radiata de la provincia de Lugo con estados de desarrollo de latizal y fustal y con ocupación de la especie superior al 60 % según el Mapa Forestal de España. Luego, se elaboró un mapa espacialmente explícito que describe V a escala provincial. En 2009 la media poblacional estimada para V fue de 166,3 m 3 ha -1 , con un error relativo de muestreo del 13,9 %. Esta metodología permitirá gener...