Fusion of airborne LiDAR and multispectral sensors reveals synergic capabilities in forest structure characterization

Forest stand structure is an important concept for ecology and planning in sustainable forest management. In this article, we consider that the incorporation of complementary multispectral information from optical sensors to Light Detection and Ranging (LiDAR) may be advantageous, especially through data fusion by back-projecting the LiDAR points onto the multispectral image. A multivariate data set of both LiDAR and multispectral metrics was related with a multivariate data set of stand structural variables measured in a Scots pine forest through canonical correlation analysis (CCA). Four statistically significant pairs of canonical variables were found, which explained 83.0% accumulated variance. The first pair of canonical variables related indicators of stand development, i.e. height and volume, with LiDAR height metrics. CCA also found attributes describing stand density to be related to LiDAR and spectral variables determining canopy coverage. Other canonical variables pertained to Lorenz curve-derived attributes, which are measures of within-stand tree size variability and heterogeneity, able to discriminate even-sized from uneven-sized stands. The most relevant result was to find that metrics derived from the multispectral sensor showed significant explanatory potential for the prediction of these structural attributes. Therefore, we concluded that metrics derived from the optical sensor have potential for complementing the information from the LiDAR sensor in describing structural properties of forest stands. We recommend the use of back-projecting for jointly exploiting the synergies of both sensors using similar types of metrics as they are customary in forestry applications of LiDAR.     

Tree species identification in mixed coniferous forest using airborne laser scanning

This study tests the capacity of relatively low density (<1 return/m²) airborne laser scanner data for discriminating between Douglas-fir, western larch, ponderosa pine, and lodgepole pine in a western North American montane forest and it evaluates the relative importance of intensity, height, and return type metrics for classifying tree species. Collectively, Exploratory Data Analysis, Pearson Correlation, ANOVA, and Linear Discriminant Analysis show that structural and intensity characteristics generated from LIDAR data are useful for classifying species at dominant and individual tree levels in multi-aged, mixed conifer forests. Proportions of return types and mean intensities are significantly different between species ($p$-value < 0.001) for plot-level dominant species and individual trees. Classification accuracies based on single variables range from 49%–61% at the dominant species level and 37%–52% for individual trees. The accuracy can be improved to 95% and 68% respectively by using multiple variables. The inclusion of proportion of return type greatly improves the classification accuracy at the dominant species level, but not for individual trees, while canopy height improves the accuracy at both levels. Overall differences in intensity and return type between species largely reflect variations in the physical structure of trees and stands. These results are consistent with the findings of others and point to airborne laser scanning as a useful source of data for species classification. However, there are still many knowledge gaps that prevent accurate mapping of species using ALS data alone, particularly with relatively sparse datasets like the one used in this study. Further investigations using other datasets in different forest types will likely result in improvements to species identification and mapping for some time to come.     

Tree species classification using UAS-based digital aerial photogrammetry point clouds and multispectral imageries in subtropical natural forests

Tree species composition of forest stand is an important indicator of forest inventory attributes for assessing ecosystem health, understanding successional processes, and digitally displaying forest biodiversity. In this study, we acquired high spatial resolution multispectral and RGB imagery over a subtropical natural forest in southwest China using a fixed-wing UAV system. Digital aerial photogrammetric (DAP) technique was used to generate multi-spectral and RGB derived point clouds, upon which individual tree crown (ITC) delineation algorithms and a machine learning classifier were used to identify dominant tree species. To do so, the structure-from-motion method was used to generate RGB imagery-based DAP point clouds. Then, three ITC delineation algorithms (i.e., point cloud segmentation (PCS), image-based multiresolution segmentation (IMRS), and advanced multiresolution segmentation (AMRS)) were used and assessed for ITC detection. Finally, tree-level metrics (i.e., multispectral, texture and point cloud metrics) were used as metrics in the random forest classifier used to classify eight dominant tree species. Results indicated that the accuracy of the AMRS ITC segmentation was highest (F₁-score = 82.5 %), followed by the segmentation using PCS (F₁-score = 79.6 %), the IMRS exhibited the lowest accuracy (F₁-score = 78.6 %); forest types classification (coniferous and deciduous) had a higher accuracy than the classification of all eight tree species, and the combination of spectral, texture and structural metrics had the highest classification accuracy (overall accuracy = 80.20 %). In the classification of both eight tree species and two forest types, the classification accuracies were lowest when only using spectral metrics, indicated that the texture metrics and point cloud structural metrics had a positive impact on the classification (the overall accuracy and kappa accuracy increased by 1.49–4.46 % and 2.86–6.84 %, respectively).     

Tree species classification in subtropical forests using small-footprint full-waveform LiDAR data

The accurate classification of tree species is critical for the management of forest ecosystems, particularly subtropical forests, which are highly diverse and complex ecosystems. While airborne Light Detection and Ranging (LiDAR) technology offers significant potential to estimate forest structural attributes, the capacity of this new tool to classify species is less well known. In this research, full-waveform metrics were extracted by a voxel-based composite waveform approach and examined with a Random Forests classifier to discriminate six subtropical tree species (i.e., Masson pine (Pinus massoniana Lamb.)), Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.), Slash pines (Pinus elliottii Engelm.), Sawtooth oak (Quercus acutissima Carruth.) and Chinese holly (Ilex chinensis Sims.) at three levels of discrimination. As part of the analysis, the optimal voxel size for modelling the composite waveforms was investigated, the most important predictor metrics for species classification assessed and the effect of scan angle on species discrimination examined. Results demonstrate that all tree species were classified with relatively high accuracy (68.6% for six classes, 75.8% for four main species and 86.2% for conifers and broadleaved trees). Full-waveform metrics (based on height of median energy, waveform distance and number of waveform peaks) demonstrated high classification importance and were stable among various voxel sizes. The results also suggest that the voxel based approach can alleviate some of the issues associated with large scan angles. In summary, the results indicate that full-waveform LIDAR data have significant potential for tree species classification in the subtropical forests.     

亚热带森林参数的机载激光雷达估测

通过应用机载激光雷达数据,在分析云南省中部的78块样地的基础上提出2个预测森林不同生物特性的统计模型(加权平均高度的预测模型和生物量的预测模型),并讨论了预测结果及其精确性。从激光雷达数据中提取了2组变量(树冠高度变量组和植被密度变量组)作为自变量,采用逐步回归方法进行自变量选择。结果表明,激光雷达数据与森林的平均树高和地上各部分生物量有很强的相关性。对于3种不同森林类型(针叶林,阔叶林和混交林),平均树高估测均能达到比较高的精度;生物量的估测结果是针叶林优于阔叶林,混交林的生物量与激光雷达数据则没有明显相关性。最后,对回归分析的结果和影响预测精度的因素进行讨论,认为预测结果的精度可能与森林类型、激光雷达采样时间和采样密度以及坐标误差等因素有关。     

The effect of leaf-on and leaf-off forest canopy conditions on LiDAR derived estimations of forest structural diversity

Forest structural diversity metrics describing diversity in tree size and crown shape within forest stands can be used as indicators of biodiversity. These diversity metrics can be generated using airborne laser scanning (LiDAR) data to provide a rapid and cost effective alternative to ground-based inspection. Measures of tree height derived from LiDAR can be significantly affected by the canopy conditions at the time of data collection, in particular whether the canopy is under leaf-on or leaf-off conditions, but there have been no studies of the effects on structural diversity metrics. The aim of this research is to assess whether leaf-on/leaf-off changes in canopy conditions during LiDAR data collection affect the accuracy of calculated forest structural diversity metrics. We undertook a quantitative analysis of LiDAR ground detection and return height, and return height diversity from two airborne laser scanning surveys collected under leaf-on and leaf-off conditions to assess initial dataset differences. LiDAR data were then regressed against field-derived tree size diversity measurements using diversity metrics from each LiDAR dataset in isolation and, where appropriate, a mixture of the two. Models utilising leaf-off LiDAR diversity variables described DBH diversity, crown length diversity and crown width diversity more successfully than leaf-on (leaf-on models resulted in R² values of 0.66, 0.38 and 0.16, respectively, and leaf-off models 0.67, 0.37 and 0.23, respectively). When LiDAR datasets were combined into one model to describe tree height diversity and DBH diversity the models described 75% and 69% of the variance (R² of 0.75 for tree height diversity and 0.69 for DBH diversity). The results suggest that tree height diversity models derived from airborne LiDAR, collected (and where appropriate combined) under any seasonal conditions, can be used to differentiate between simple single and diverse multiple storey forest structure with confidence.     

Estimation of aboveground biomass of Robinia pseudoacacia forest in the Yellow River Delta based on UAV and Backpack LiDAR point clouds

Forest plantations are an important source of terrestrial carbon sequestration. The forest of Robinia pseudoacacia in the Yellow River Delta (YRD) is the largest artificial ecological protection forest in China. However, more than half of the forest has appeared different degrees of dieback and even death since the 1990s. Timely and accurate estimation of the forest aboveground biomass (AGB) is a basis for studying the carbon cycle of forests. Light Detecting and Ranging (LiDAR) has been proved to be one of the most powerful methods for forest biomass estimation. However, because of an irregular and overlapping shape of the broadleaved forest canopy in a growing season, it is difficult to segment individual trees and estimate the tree biomass from airborne LiDAR data. In this study, a new method was proposed to solve this problem of individual tree detection in the Robinia pseudoacacia forest based on a combination of the Unmanned Aerial Vehicle-Light Detecting and Ranging (UAV-LiDAR) with the Backpack-LiDAR. The proposed method mainly consists of following steps: (i) at a plot level, trees in the UAV-LiDAR data were detected by seed points obtained by an individual tree segmentation (ITS) method from the Backpack-LiDAR data; (ii) height and diameter at breast height (DBH) of an individual tree would be extracted from UAV and Backpack LiDAR data, respectively; (iii) the individual tree AGB would be calculated through an allometric equation and the forest AGB at the plot level was accumulated; and (iv) the plot-level forest AGB was taken as a dependent variable, and various metrics extracted from UAV-LiDAR point cloud data as independent variables to estimate forest AGB distribution in the study area by using both multiple linear regression (MLR) and random forest (RF) models. The results demonstrate that: (1) the seed points extracted from Backpack-LiDAR could significantly improve the overall accuracy of individual tree detection (F = 0.99), and thus increase the forest AGB estimation accuracy; (2) compared with MLR model, the RF model led to a higher estimation accuracy (p < 0.05); and (3) LiDAR intensity information selected by both MLR and RF models and laser penetration rate (LP) played an important role in estimating healthy forest AGB.     

基于Landsat长时间序列数据估算树高和生物量

以Landsat长时间序列数据为研究对象,旨在以光谱序列信息反演森林参数为视角,应用Landtrendr算法从时间序列数据中提取森林扰动变量,使用随机森林计算方法建立扰动变量、反射率和GLAS激光点森林参数之间的关系模型,获取树高和生物量的空间分布信息。为多源遥感数据反演森林参数提供参考,研究证明基于Landsat长时间序列数据获得的森林扰动变量能够增强反射率和森林参数之间的相关性,可提高预测精度。     

Vegetation discrimination using irs-p3 wifs temporal dataset - a case study from rampa forests, eastern ghats, a.p.

Multitemporal data sets from coarse resolution sensors of Indian Remote Sensing Satellites provides an opportunity to classify various forest types using their phenological attributes reflected in temporal NDVI profiles. The present study attempts to classify various vegetation classes using time integrated NDVI (T-NDVI) values derived from IRS-P3 WiFS data. The algorithm explores the differential characteristics in T-NDVI values of different features and the results suggest the possible use of the methodology for forest type classification.     

Mapping forest biomass from space – Fusion of hyperspectral EO1-hyperion data and Tandem-X and WorldView-2 canopy height models

Spaceborne sensors allow for wide-scale assessments of forest ecosystems. Combining the products of multiple sensors is hypothesized to improve the estimation of forest biomass. We applied interferometric (Tandem-X) and photogrammetric (WorldView-2) based predictors, e.g. canopy height models, in combination with hyperspectral predictors (EO1-Hyperion) by using 4 different machine learning algorithms for biomass estimation in temperate forest stands near Karlsruhe, Germany. An iterative model selection procedure was used to identify the optimal combination of predictors. The most accurate model (Random Forest) reached a r² of 0.73 with a RMSE of 14.9% (29.4 t/ha). Further results revealed that the predictive accuracy depended highly on the statistical model and the area size of the field samples. We conclude that a fusion of canopy height and spectral information allows for accurate estimations of forest biomass from space.     

Predicting individual tree attributes from airborne laser point clouds based on the random forests technique

This paper depicts an approach for predicting individual tree attributes, i.e., tree height, diameter at breast height (DBH) and stem volume, based on both physical and statistical features derived from airborne laser-scanning data utilizing a new detection method for finding individual trees together with random forests as an estimation method. The random forests (also called regression forests) technique is a nonparametric regression method consisting of a set of individual regression trees. Tests of the method were performed, using 1476 trees in a boreal forest area in southern Finland and laser data with a density of 2.6 points per m². Correlation coefficients (R) between the observed and predicted values of 0.93, 0.79 and 0.87 for individual tree height, DBH and stem volume, respectively, were achieved, based on 26 laser-derived features. The corresponding relative root-mean-squared errors (RMSEs) were 10.03%, 21.35% and 45.77% (38% in best cases), which are similar to those obtained with the linear regression method, with maximum laser heights, laser-estimated DBH or crown diameters as predictors. With random forests, however, the forest models currently used for deriving the tree attributes are not needed. Based on the results, we conclude that the method is capable of providing a stable and consistent solution for determining individual tree attributes using small-footprint laser data.     

Characterization of aboveground biomass in an unmanaged boreal forest using Landsat temporal segmentation metrics

Canada is dominated by forested ecosystems which are subject to various inventory and management practices, with more northern boreal forests subject to neither. Our objectives were to measure the capacity of temporal trajectory metrics for estimating selected forest attributes in a northern Canadian boreal forest context using Landsat imagery and investigate the importance of different types of temporal trajectory metrics. Results indicated that Wetness was the best Tasseled Cap (TC) component for aboveground biomass estimation (R² = 50%, RMSE% = 56%), and the combination of simple and complex metrics from all TC components produced the highest R² (62%) and lowest RMSE% (49%). Using a similar combination of variables, other forest attributes were estimated equally reliably with lower RMSE% values. The most important temporal trajectory metrics were simple and described TC component values at each point of change in the temporal trajectory, however the most important variables overall were environmental variables.     

基于MODIS影像多特征的CART决策树分类

以山东省为研究区域,利用2009年9月MODIS的8 d合成波段反射率产品MOD09,选择特征变量植被指数(NDVI、EVI)、NDWI、NDMI、NDSI及辅助信息DEM,通过选取其中的影像特征组合来确定分类方案,构建各波段组合的CART决策树,对MODIS影像进行分类,得到CART决策树的最优波段组合。结果表明,特征变量DEM、NDVI、EVI对分类结果贡献较大;将CART决策树的分类结果与其相对应的最大似然分类结果进行比较可知,基于影像多特征的CART决策树分类方法能明显提高分类精度。     

Classification of forest land attributes using multi-source remotely sensed data

The aim of the study was to (1) examine the classification of forest land using airborne laser scanning (ALS) data, satellite images and sample plots of the Finnish National Forest Inventory (NFI) as training data and to (2) identify best performing metrics for classifying forest land attributes. Six different schemes of forest land classification were studied: land use/land cover (LU/LC) classification using both national classes and FAO (Food and Agricultural Organization of the United Nations) classes, main type, site type, peat land type and drainage status. Special interest was to test different ALS-based surface metrics in classification of forest land attributes. Field data consisted of 828 NFI plots collected in 2008–2012 in southern Finland and remotely sensed data was from summer 2010. Multinomial logistic regression was used as the classification method. Classification of LU/LC classes were highly accurate (kappa-values 0.90 and 0.91) but also the classification of site type, peat land type and drainage status succeeded moderately well (kappa-values 0.51, 0.69 and 0.52). ALS-based surface metrics were found to be the most important predictor variables in classification of LU/LC class, main type and drainage status. In best classification models of forest site types both spectral metrics from satellite data and point cloud metrics from ALS were used. In turn, in the classification of peat land types ALS point cloud metrics played the most important role. Results indicated that the prediction of site type and forest land category could be incorporated into stand level forest management inventory system in Finland.     

机器学习算法在森林地上生物量估算中的应用

森林地上生物量是森林生产力的重要评价指标,对其进行高效监测对维持全球碳平衡和保护生态系统具有重要意义。本文首先基于冠层高度模型数据,通过分水岭分割算法得到单木冠幅边界;然后在单木冠幅范围内提取23个LiDAR变量,结合佩诺布斯科特试验森林的87组实测数据,利用随机森林和支持向量机建立森林地上生物量估算模型;最后对样地模型估算的结果进行了比较,讨论了预测结果及其精度。结果表明:本文选用的随机森林模型和支持向量机模型在估算森林地上生物量的应用中获得了较高的精度;并且,随机森林模型在基于机载雷达数据估测森林地上生物量中的估算精度更高,模型泛化能力更强,制图精度也更好,具有更好的适用性。     

Mapping beech (Fagus sylvatica L.) forest structure with airborne hyperspectral imagery

Estimating forest structural attributes using multispectral remote sensing is challenging because of the saturation of multispectral indices at high canopy cover. The objective of this study was to assess the utility of hyperspectral data in estimating and mapping forest structural parameters including mean diameter-at-breast height (DBH), mean tree height and tree density of a closed canopy beech forest (Fagus sylvatica L.). Airborne HyMap images and data on forest structural attributes were collected from the Majella National Park, Italy in July 2004. The predictive performances of vegetation indices (VI) derived from all possible two-band combinations (VI_(i,j) = (R_i − R_j)/(R_i + R_j), where R_i and R_j = reflectance in any two bands) were evaluated using calibration (n = 33) and test (n = 20) data sets. The potential of partial least squares (PLS) regression, a multivariate technique involving several bands was also assessed. New VIs based on the contrast between reflectance in the red-edge shoulder (756–820 nm) and the water absorption feature centred at 1200 nm (1172–1320 nm) were found to show higher correlations with the forest structural parameters than standard VIs derived from NIR and visible reflectance (i.e. the normalised difference vegetation index, NDVI). PLS regression showed a slight improvement in estimating the beech forest structural attributes (prediction errors of 27.6%, 32.6% and 46.4% for mean DBH, height and tree density, respectively) compared to VIs using linear regression models (prediction errors of 27.8%, 35.8% and 48.3% for mean DBH, height and tree density, respectively). Mean DBH was the best predicted variable among the stand parameters (calibration R² = 0.62 for an exponential model fit and standard error of prediction = 5.12 cm, i.e. 25% of the mean). The predicted map of mean DBH revealed high heterogeneity in the beech forest structure in the study area. The spatial variability of mean DBH occurs at less than 450 m. The DBH map could be useful to forest management in many ways, e.g. thinning of coppice to promote diameter growth, to assess the effects of management on forest structure or to detect changes in the forest structure caused by anthropogenic and natural factors.     

Individual tree crown approach for predicting site index in boreal forests using airborne laser scanning and hyperspectral data

Site productivity is essential information for sustainable forest management and site index (SI) is the most common quantitative measure of it. The SI is usually determined for individual tree species based on tree height and the age of the 100 largest trees per hectare according to stem diameter. The present study aimed to demonstrate and validate a methodology for the determination of SI using remotely sensed data, in particular fused airborne laser scanning (ALS) and airborne hyperspectral data in a forest site in Norway. The applied approach was based on individual tree crown (ITC) delineation: tree species, tree height, diameter at breast height (DBH), and age were modelled and predicted at ITC level using 10-fold cross validation. Four dominant ITCs per 400 m² plot were selected as input to predict SI at plot level for Norway spruce (Picea abies (L.) Karst.) and Scots pine (Pinus sylvestris L.). We applied an experimental setup with different subsets of dominant ITCs with different combinations of attributes (predicted or field-derived) for SI predictions. The results revealed that the selection of the dominant ITCs based on the largest DBH independent of tree species, predicted the SI with similar accuracy as ITCs matched with field-derived dominant trees (RMSE: 27.6% vs 23.3%). The SI accuracies were at the same level when dominant species were determined from the remotely sensed or field data (RMSE: 27.6% vs 27.8%). However, when the predicted tree age was used the SI accuracy decreased compared to field-derived age (RMSE: 27.6% vs 7.6%). In general, SI was overpredicted for both tree species in the mature forest, while there was an underprediction in the young forest. In conclusion, the proposed approach for SI determination based on ITC delineation and a combination of ALS and hyperspectral data is an efficient and stable procedure, which has the potential to predict SI in forest areas at various spatial scales and additionally to improve existing SI maps in Norway.     

The role of ground reference data collection in the prediction of stem volume with LiDAR data in mountain areas

Ground reference data collection represents an important element in the prediction of stem volume with LiDAR-derived variables, and at present it is the most expensive part of such analyses. In this paper two aspects of ground reference data collection were analyzed: (1) the positioning error of the ground plots; and (2) the optimal number of training plots. A system for the prediction of stem volume at area-based level was adopted. LiDAR data were preprocessed and 13 variables describing both height and coverage were extracted. Models were defined using a stepwise ordinary least square (OLS) regression. Three experiments were conducted: (i) the role of the plots positioning error on prediction accuracy; (ii) the influence of random downsampling of plot numbers on prediction accuracy; and (iii) the influence of a stratified downsampling of plot numbers on prediction accuracy based on LiDAR-derived variables. A dataset comprising 799 ground plots was used. They were distributed throughout a mountainous area in the Southern Alps, where the presence of a complex landscape increases the uncertainty of the Global Positioning System (GPS) accuracy, and where a large variety of tree forest species and climatic environments make it necessary to have a large number of sample plots for accurate characterization of the study area. All the experiments provided important indications for LiDAR based forest inventories: the GPS error did not significantly influence the prediction accuracy and it was possible to reduce the number of training samples without compromising the generalization ability of the prediction model. Leading on from these findings, a new ground sampling protocol based on genetic algorithms was proposed. The new protocol allowed us to obtain promising results for the considered dataset: using only 53 training plots, instead of 534 in the original dataset, we obtained the same results for the validation set. These results, obtained in a complex mountainous area, are representative of Alpine environments and allow us to infer that similar (or better) results could also be obtained within non mountainous areas.     

Can regional aerial images from orthophoto surveys produce high quality photogrammetric Canopy Height Model? A single tree approach in Western Europe

Forest monitoring tools are needed to promote effective and data driven forest management and forest policies. Remote sensing techniques can increase the speed and the cost-efficiency of the forest monitoring as well as large scale mapping of forest attribute (wall-to-wall approach). Digital Aerial Photogrammetry (DAP) is a common cost-effective alternative to airborne laser scanning (ALS) which can be based on aerial photos routinely acquired for general base maps. DAP based on such pre-existing dataset can be a cost effective source of large scale 3D data. In the context of forest characterization, when a quality Digital Terrain Model (DTM) is available, DAP can produce photogrammetric Canopy Height Model (pCHM) which describes the tree canopy height. While this potential seems pretty obvious, few studies have investigated the quality of regional pCHM based on aerial stereo images acquired by standard official aerial surveys. Our study proposes to evaluate the quality of pCHM individual tree height estimates based on raw images acquired following such protocol using a reference filed-measured tree height database. To further ensure the replicability of the approach, the pCHM tree height estimates benchmarking only relied on public forest inventory (FI) information and the photogrammetric protocol was based on low-cost and widely used photogrammetric software. Moreover, our study investigates the relationship between the pCHM tree height estimates based on the neighboring forest parameter provided by the FI program.Our results highlight the good agreement of tree height estimates provided by pCHM using DAP with both field measured and ALS tree height data. In terms of tree height modeling, our pCHM approach reached similar results than the same modeling strategy applied to ALS tree height estimates. Our study also identified some of the drivers of the pCHM tree height estimate error and found forest parameters like tree size (diameter at breast height) and tree type (evergreenness/deciduousness) as well as the terrain topography (slope) to be of higher importance than image survey parameters like the variation of the overlap or the sunlight condition in our dataset. In combination with the pCHM tree height estimate, the terrain slope, the Diameter at Breast Height (DBH) and the evergreenness factor were used to fit a multivariate model predicting the field measured tree height. This model presented better performance than the model linking the pCHM estimates to the field tree height estimates in terms of r² (0.90 VS 0.87) and root mean square error (RMSE, 1.78 VS 2.01 m). Such aspects are poorly addressed in literature and further research should focus on how pCHM approaches could integrate them to improve forest characterization using DAP and pCHM. Our promising results can be used to encourage the use of regional aerial orthophoto surveys archive to produce large scale quality tree height data at very low additional costs, notably in the context of updating national forest inventory programs.