Estimating plot volume using lidar height and intensity distributional parameters

This study explored the feasibility of height distributional metrics and intensity values extracted from low-density airborne light detection and ranging (lidar) data to estimate plot volumes in dense Korean pine (Pinus koraiensis) plots. Multiple linear regression analyses were performed using lidar height and intensity distributional metrics. The candidate variables for predicting plot volume were evaluated using three data sets: total, canopy, and integrated lidar height and intensity metrics. All intensities of lidar returns used were corrected by the reference distance. Regression models were developed using each data set, and the first criterion used to select the best models was the corrected Akaike Information Criterion (AIC_c). The use of three data sets was statistically significant at R² = 0.75 (RMSE = 52.17 m³ ha^?1), R² = 0.84 (RMSE = 45.24 m³ ha^?1), and R² = 0.91 (RMSE = 31.48 m³ ha^?1) for total, canopy, and integrated lidar distributional metrics, respectively. Among the three data sets, the integrated lidar metrics-derived model showed the best performance for estimating plot volumes, improving errors up to 42% when compared to the other two data sets. This is attributed to supplementing variables weighted and biased to upper limits in dense plots with more statistical variables that explain the lower limits. In all data sets, intensity metrics such as skewness, kurtosis, standard deviation, minimum, and standard error were employed as explanatory variables. The use of intensity variables improved the accuracy of volume estimation in dense forests compared to prior research. Correction of the intensity values contributed up to a maximum of 58% improvement in volume estimation when compared to the use of uncorrected intensity values (R² = 0.78, R² = 0.53, and R² = 0.63 for total, canopy, and integrated lidar distributional metrics, respectively). It is clear that the correction of intensity values is an essential step for the estimation of forest volume.     

Species stratification and upscaling of forest carbon estimates to landscape scale using GeoEye-1 image and lidar data in sub-tropical forests of Nepal

Aboveground forest biomass and carbon estimation at landscape scale is crucial for implementation of REDD+ programmes. This study aims to upscale the forest carbon estimates using GeoEye-1 image and small footprint lidar data from small areas to a landscape level using RapidEye image. Species stratification was carried out based on the spectral separability curve of GeoEye-1 image, and comparison of mean intensity and mean plot height of the trees from lidar data. GeoEye-1 image and lidar data were segmented using region growing approach to delineate individual tree crowns; and the segmented crowns (CPA) of tree were further used to establish a relationship with field measured carbon and total trees’ height. Carbon stock measured from field, individual tree crown (ITC) segmentation approach and area-based approach (ABA) was compared at plot level using one-way ANOVA and post hoc Tukey comparison test. ITC-based carbon estimates was used to establish a relationship with spectral reflectance of RapidEye image variables (NDVI, RedEdge NDVI, PC1, single band of RedEdge, and NIR) to upscale the carbon at landscape level. One-way ANOVA resulted in a highly significant difference (p-value < 0.005) between the mean plot height and lidar intensity to stratify Shorea robusta and Other species successfully. ITC carbon stock estimation models of two major tree species explained about 88% and 79% of the variances, respectively, at 95% confidence level. The ABA estimated carbon was highly correlated (R² = 0.83, RMSE = 20.04) to field measured carbon with higher accuracy than the ITC estimated carbon. A weak relationship was observed between the carbon stock and the RapidEye image variables. However, upscaling of carbon estimates from ABA is likely to improve the relationship of the RapidEye variables rather than upscaling the carbon estimates from ITC approach.     

A pseudo-waveform technique to assess forest structure using discrete lidar data

The use of airborne laser scanning systems (lidar) to describe forest structure has increased dramatically since height profiling experiments nearly 30 years ago. The analyses in most studies employ a suite of frequency-based metrics calculated from the lidar height data, which are systematically eliminated from a full model using stepwise multiple linear regression. The resulting models often include highly correlated predictors with little physical justification for model formulation. We propose a method to aggregate discrete lidar height and intensity measurements into larger footprints to create “pseudo-waves”. Specifically, the returns are first sorted into height bins, sliced into narrow discrete elements, and finally smoothed using a spline function. The resulting “pseudo-waves” have many of the same characteristics of traditional waveform lidar data. We compared our method to a traditional frequency-based method to estimate tree height, canopy structure, stem density, and stand biomass in coniferous and deciduous stands in northern Wisconsin (USA). We found that the pseudo-wave approach had strong correlations for nearly all tree measurements including height (cross validated adjusted R² (R²_cv) = 0.82, RMSE_cv = 2.09 m), mean stem diameter (R²_cv = 0.64, RMSE_cv = 6.15 cm), total aboveground biomass (R²_cv = 0.74, RMSE_cv = 74.03 kg ha^− 1), and canopy coverage (R²_cv = 0.79, RMSE_cv = 5%). Moreover, the type of wave (derived from height and intensity or from height alone) had little effect on model formulation and fit. When wave-based and frequency-based models were compared, fit and mean square error were comparable, leading us to conclude that the pseudo-wave approach is a viable alternative because it has 1) an increased breadth of available metrics; 2) the potential to establish new meaningful metrics that capture unique patterns within the waves; 3) the ability to explain metric selection based on the physical structure of forests; and 4) lower correlation among independent variables.     

A mixed pixel- and region-based approach for using airborne laser scanning data for individual tree crown delineation in Pinus radiata D. Don plantations

The aim of this study was to evaluate the use of high-resolution airborne laser scanner (ALS) data to detect and measure individual trees. We developed and tested a new mixed pixel- and region-based algorithm (using Definiens Developer 7.0) for locating individual tree positions and estimating their total heights. We computed a canopy height model (CHM) of pixel size 0.25 m from dense first-pulse point data (8 pulses m^?2) acquired with a small-footprint discrete-return lidar sensor. We validated the results of individual tree segmentation with accurate field measurements made in 37 plots of Monterey pine (Pinus radiata D. Don) distributed over an area of 36 km². Fieldwork consisted of labelling all of the trees in each plot and measuring their height and position, for posterior integration of the data from both sources (field and lidar). The proposed algorithm correctly detected and linked 59.8% of the trees in the 37 sample plots. We also manually located the trees by using FUSION software to visualize the raw lidar data cloud. However, because the latter method is extremely time-consuming, we only considered 10 randomly selected plots. Manual location correctly detected and linked 71.9% of the trees (in this subsample the algorithm correctly detected and measured 63.5% of the trees). The R² values for the linear model relating field- and lidar-measured heights of the linked trees located manually and with the automatic location algorithm were 0.90 and 0.88, respectively.     

Use of Large-Footprint Scanning Airborne Lidar To Estimate Forest Stand Characteristics in the Western Cascades of Oregon

A scanning lidar, a relatively new type of sensor which explicitly measures canopy height, was used to measure structure of conifer forests in the Pacific Northwest. SLICER (Scanning Lidar Imager of Canopies by Echo Recovery), an airborne pulsed laser developed by NASA which scans a swath of five 10-m diameter footprints along the aircraft’s flightpath, captures the power of the reflected laser pulse as a function of height from the top of the canopy to the ground. Ground measurements of forest stand structure were collected on 26 plots with coincident SLICER data. Height, basal area, total biomass, and leaf biomass as estimated from field data could be predicted from SLICER-derived metrics with r² values of 0.95, 0.96, 0.96, and 0.84, respectively. These relationships were strong up to a height of 52 m, basal area of 132 m²/ha and total biomass of 1300 Mg/ha. In light of these strong relationships, large-footprint, airborne scanning lidar shows promise for characterizing stand structure for management and research purposes.     

An airborne lidar sampling strategy to model forest canopy height from Quickbird imagery and GEOBIA

High-resolution digital canopy models derived from airborne lidar data have the ability to provide detailed information on the vertical structure of forests. However, compared to satellite data of similar spatial resolution and extent, the small footprint airborne lidar data required to produce such models remain expensive. In an effort to reduce these costs, the primary objective of this paper is to develop an airborne lidar sampling strategy to model full-scene forest canopy height from optical imagery, lidar transects and Geographic Object-Based Image Analysis (GEOBIA). To achieve this goal, this research focuses on (i) determining appropriate lidar transect features (i.e., location, direction and extent) from an optical scene, (ii) developing a mechanism to model forest canopy height for the full-scene based on a minimum number of lidar transects, and (iii) defining an optimal mean object size (MOS) to accurately model the canopy composition and height distribution. Results show that (i) the transect locations derived from our optimal lidar transect selection algorithm accurately capture the canopy height variability of the entire study area; (ii) our canopy height estimation models have similar performance in two lidar transect directions (i.e., north-south and west-east); (iii) a small lidar extent (17.6% of total size) can achieve similar canopy height estimation accuracies as those modeled from the full lidar scene; and (iv) different MOS can lead to distinctly different canopy height results. By comparing the best canopy height estimate with the full lidar canopy height data, we obtained average estimation errors of 6.0 m and 6.8 m for conifer and deciduous forests at the individual tree crown/small tree cluster level, and an area weighted combined error of 6.2 m, which is lower than the provincial forest inventory height class interval (i.e., ≈ 9.0 m).     

Measuring forest structure and biomass in New England forest stands using Echidna ground-based lidar

A ground-based, upward-scanning, near-infrared lidar, the Echidna® validation instrument (EVI), built by CSIRO Australia, retrieves forest stand structural parameters, including mean diameter at breast height (DBH), stem count density (stems/area), basal area, and above-ground woody biomass with very good accuracy in six New England hardwood and conifer forest stands. Comparing forest structural parameters retrieved using EVI data with extensive ground measurements, we found excellent agreement at the site level using five EVI scans (plots) per site (R² = 0.94-0.99); very good agreement at the plot level for stem count density and biomass (R² = 0.90-0.85); and good agreement at the plot level for mean DBH and basal area (R² = 0.48-0.66). The observed variance at site and plot levels suggest that a sample area of at least 1 ha (10⁴ m²) is required to estimate these parameters accurately at the stand level using either lidar-based or conventional methods. The algorithms and procedures used to retrieve these structural parameters are dependent on the unique ability of the Echidna® lidar to digitize the full waveform of the scattered lidar pulse as it returns to the instrument, which allows consistent separation of scattering by trunks and large branches from scattering by leaves. This successful application of ground-based lidar technology opens the door to rapid and accurate measurement of biomass and timber volume in areal sampling scenarios and as a calibration and validation tool for mapping biomass using airborne or spaceborne remotely sensed data.     

Satellite lidar vs. small footprint airborne lidar: Comparing the accuracy of aboveground biomass estimates and forest structure metrics at footprint level

The use of lidar remote sensing for mapping the spatial distribution of canopy characteristics has the potential to allow an accurate and efficient estimation of tree dimensions and canopy structural properties from local to regional and continental scales. The overall goal of this paper was to compare biomass estimates and height metrics obtained by processing GLAS waveform data and spatially coincident discrete-return airborne lidar data over forest conditions in east Texas. Since biomass estimates are derived from waveform height metrics, we also compared ground elevation measurements and canopy parameters. More specific objectives were to compare the following parameters derived from GLAS and airborne lidar: (1) ground elevations; (2) maximum canopy height; (3) average canopy height; (4) percentiles of canopy height; and (5) above ground biomass. We used the elliptical shape of GLAS footprints to extract canopy height metrics and biomass estimates derived from airborne lidar. Results indicated a very strong correlation for terrain elevations between GLAS and airborne lidar, with an r value of 0.98 and a root mean square error of 0.78 m. GLAS height variables were able to explain 80% of the variance associated with the reference biomass derived from airborne lidar, with an RMSE of 37.7 Mg/ha. Most of the models comparing GLAS and airborne lidar height metrics had R-square values above 0.9.     

Investigating the impact of discrete-return lidar point density on estimations of mean and dominant plot-level tree height in Eucalyptus grandis plantations

The accuracy of lidar remote sensing in characterizing three-dimensional forest structural attributes has encouraged foresters to integrate lidar approaches in routine inventories. However, lidar point density is an important consideration when assessing forest biophysical parameters, given the direct relationship between higher spatial resolution and lidar acquisition and processing costs. The aim of this study was to investigate the effect of point density on mean and dominant tree height estimates at plot level. The study was conducted in an intensively managed Eucalyptus grandis plantation. High point density (eight points/m²) discrete-return, small-footprint lidar data were used to generate point density simulations averaging 0.25, one, two, three, four, five, and six points/m². Field surveyed plot-level mean and dominant heights were regressed against metrics derived from lidar data at each simulated point density. Stepwise regression was used to identify which lidar metrics produced the best models. Mean height was estimated at accuracy of R² ranging between 0.93 and 0.94 while dominant height was estimated with an R² of 0.95. Root mean square error (RMSE) was also similar at all densities for mean height (～1.0 m) and dominant height (～1.2 m); the relative RMSE compared to field-measured mean was constant at approximately 5%. Analysis of bias showed that the estimation of both variables did not vary with density. The results indicated that all lidar point densities resulted in reliable models. It was concluded that plot-level height can be estimated with reliable accuracy using relatively low density lidar point spacing. Additional research is required to investigate the effect of low point density on estimation of other forest biophysical attributes.     

Estimation of above-ground forest biomass using metrics based on Gaussian decomposition of waveform lidar data

Large-footprint waveform light detection and ranging (lidar) data have been widely used in above-ground forest biomass estimation. Waveform metrics derived from basic statistics (e.g. percentile of energy) of the lidar waveform, such as canopy height and height of median energy, have been applied to biomass estimation in numerous studies. In this study, a set of metrics based on Gaussian decomposition (GD) results were developed and evaluated for forest above-ground biomass estimation using NASA’s laser vegetation imaging sensor (LVIS) data. The GD metrics were designed to explicitly incorporate lidar intensity and vertical structures of canopy layers for biomass estimation. The proposed GD metrics used information related to the above-ground height of each Gaussian centroid and the Gaussian area index (GAI), where GAI is the area covered by a Gaussian function. Two types of novel GD metrics were developed: (1) percentile-height GAI metrics expressing the GAI summation of a subset of Gaussian centroids located within a certain percentile height range; and (2) height-weighted GAI metrics, a summation of GAIs of a waveform weighted by the corresponding heights of their Gaussian centroids. A biomass regression model was built by eight newly developed GD metrics using GAI information and five pre-existing GD-derived metrics that have not previously been used for biomass estimation. The performance of the regression model was then compared to another regression model using 12 previously published metrics (non-GD metrics). The Random Forests (RF) regression algorithm was employed for predicting biomass. The RF out-of-bag results indicated that above-ground biomass estimations using GD metrics achieved significantly better results than those derived from non-GD metrics for deciduous plots (19% lower root mean square error (RMSE), 25% higher coefficient of determination (R²), and marginally better results in coniferous plots (4% lower RSME, 6% higher R²). The combination of GD and non-GD metrics achieved comparable biomass estimation results to the model using exclusively GD metrics. GD metrics also showed strong correlation with forest attributes such as mean diameter at breast height (DBH) and stem density. This study contributes to the usage of GD results for accurate estimation of forest above-ground biomass in large-footprint lidar waveform data in temperate deciduous forests, because temperate deciduous forests have been proved challenging in regard to lidar-derived biomass estimations.     

A voxel-based lidar method for estimating crown base height for deciduous and pine trees

The overall goal of this study was to develop methods for assessing crown base height for individual trees using airborne lidar data in forest settings typical for the southeastern United States. More specific objectives are to: (1) develop new lidar-derived features as multiband height bins and processing techniques for characterizing the vertical structure of individual tree crowns; (2) investigate several techniques for filtering and analyzing vertical profiles of individual trees to derive crown base height, such as Fourier and wavelet filtering, polynomial fit, and percentile analysis; (3) assess the accuracy of estimating crown base height for individual trees, and (4) investigate which type of lidar data, point frequency or intensity, provides the most accurate estimate of crown base height. A lidar software application, TreeVaW, was used to locate individual trees and to obtain per tree measurements of height and crown width. Tree locations were used with lidar height bins to derive the vertical structure of tree crowns and measurements of crown base height. Lidar-derived crown base heights of individual trees were compared to field observations for 117 trees, including 94 pines and 23 deciduous trees. Linear regression models were able to explain up to 80% of the variability associated with crown base height for individual trees. Fourier filtering used for smoothing the vertical crown profile consistently provided the best results when estimating crown base height.     

A lidar-based hierarchical approach for assessing MODIS fPAR

The purpose of this study was to estimate the fraction of photosynthetically active radiation absorbed by the canopy (fPAR) from point measurements to airborne lidar for hierarchical scaling up and assessment of the Moderate Resolution Imaging Spectroradiometer (MODIS) fPAR product within a “medium-sized” (7 km × 18 km) watershed. Nine sites across Canada, containing one or more (of 11) distinct species types and age classes at varying stages of regeneration and seasonal phenology were examined using a combination of discrete pulse airborne scanning Light Detection And Ranging (lidar) and coincident analog and digital hemispherical photography (HP). Estimates of fPAR were first compared using three methods: PAR radiation sensors, HP, and airborne lidar. HP provided reasonable estimates of fPAR when compared with radiation sensors. A simplified fractional canopy cover ratio from lidar based on the number of within canopy returns to the total number of returns was then compared with fPAR estimated from HP at 486 geographically registered measurement locations. The return ratio fractional cover method from lidar compared well with HP-derived fPAR (coefficient of determination = 0.72, RMSE = 0.11), despite varying the lidar survey configurations, canopy structural characteristics, seasonal phenologies, and possible slight inaccuracies in location using handheld GPS at some sites. Lidar-derived fractional cover estimates of fPAR were ～ 10% larger than those obtained using HP (after removing wood components), indicating that lidar likely provides a more realistic estimate of fPAR than HP when compared with radiation sensors. Finally, fPAR derived from lidar fractional cover was modelled at 1 m resolution and averaged over 99 1 km areas for comparison with MODIS fPAR. The following study is one of the first to scale between plot measurements and MODIS pixels using airborne lidar.     

Mapping canopy height using a combination of digital stereo‐photogrammetry and lidar

Ranging techniques such as lidar (LIght Detection And Ranging) and digital stereo‐photogrammetry show great promise for mapping forest canopy height. In this study, we combine these techniques to create hybrid photo‐lidar canopy height models (CHMs). First, photogrammetric digital surface models (DSMs) created using automated stereo‐matching were registered to corresponding lidar digital terrain models (DTMs). Photo‐lidar CHMs were then produced by subtracting the lidar DTM from the photogrammetric DSM. This approach opens up the possibility of retrospective mapping of forest structure using archived aerial photographs. The main objective of the study was to evaluate the accuracy of photo‐lidar CHMs by comparing them to reference lidar CHMs. The assessment revealed that stereo‐matching parameters and left–right image dissimilarities caused by sunlight and viewing geometry have a significant influence on the quality of the photo DSMs. Our study showed that photo‐lidar CHMs are well correlated to their lidar counterparts on a pixel‐wise basis (r up to 0.89 in the best stereo‐matching conditions), but have a lower resolution and accuracy. It also demonstrated that plot metrics extracted from the lidar and photo‐lidar CHMs, such as height at the 95th percentile of 20 m×20 m windows, are highly correlated (r up to 0.95 in general matching conditions).     

Development of a robust canopy height model derived from ALS point clouds for predicting individual crown attributes at the species level

ABSTRACT

Tree crown attributes are important parameters during the assessment and monitoring of forest ecosystems. Canopy height models (CHMs) derived from airborne laser scanning (ALS) data have proved to be a reliable source for extracting different biophysical characteristics of single trees and at stand level. However, ALS-derived tree measurements (e.g., mean crown diameter) can be negatively affected by pits that appear in the CHMs. Thus, we propose a novel method for generating pit-free CHMs from ALS point clouds for estimating crown attributes (i.e., area and mean diameter) at the species level. The method automatically calculates a threshold for a pixel based on the range of height values within neighbouring pixels; if the pixel falls below the threshold then it is recognized as a pitted pixel. The pit is then filled with the median of the values of the neighbouring pixels. Manually delineated individual tree crowns (ITC) of four deciduous and two coniferous species on Colour Infrared (CIR) stereo images were used as a reference in the analysis. In addition, a variety of different algorithms for constructing CHMs were compared to investigate the performance of different CHMs in similar forest conditions. Comparisons between the estimated and observed crown area (R² = 0.95, RMSE% = 19.12% for all individuals) and mean diameter (R² = 0.92, RMSE% = 12.16% for all individuals) revealed that ITC attributes were correctly estimated by segmentation of the pit-free CHM proposed in this study. The goodness of matching and geometry revealed that the delineated crowns correctly matched up to the reference data and had identical geometry in approximately 70% of cases. The results showed that the proposed method produced a CHM that estimates crown attributes more accurately than the other investigated CHMs. Furthermore, the findings suggest that the proposed algorithm used to fill pits with the median of height observed in surrounding pixels significantly improve the accuracy of the results the species level due to a higher correlation between the estimated and observed crown attributes. Based on these results, we concluded that the proposed pit filling method is capable of providing an automatic and objective solution for constructing pit-free CHMs for assessing individual crown attributes of mixed forest stands.     

Characterizing forest canopy structure with lidar composite metrics and machine learning

A lack of reliable observations for canopy science research is being partly overcome by the gradual use of lidar remote sensing. This study aims to improve lidar-based canopy characterization with airborne laser scanners through the combined use of lidar composite metrics and machine learning models. Our so-called composite metrics comprise a relatively large number of lidar predictors that tend to retain as much information as possible when reducing raw lidar point clouds into a format suitable as inputs to predictive models of canopy structural variables. The information-rich property of such composite metrics is further complemented by machine learning, which offers an array of supervised learning models capable of relating canopy characteristics to high-dimensional lidar metrics via complex, potentially nonlinear functional relationships. Using coincident lidar and field data over an Eastern Texas forest in USA, we conducted a case study to demonstrate the ubiquitous power of the lidar composite metrics in predicting multiple forest attributes and also illustrated the use of two kernel machines, namely, support vector machine and Gaussian processes (GP). Results show that the two machine learning models in conjunction with the lidar composite metrics outperformed traditional approaches such as the maximum likelihood classifier and linear regression models. For example, the five-fold cross validation for GP regression models (vs. linear/log-linear models) yielded a root mean squared error of 1.06 (2.36) m for Lorey's height, 0.95 (3.43) m for dominant height, 5.34 (8.51) m²/ha for basal area, 21.4 (40.5) Mg/ha for aboveground biomass, 6.54 (9.88) Mg/ha for belowground biomass, 0.75 (2.76) m for canopy base height, 2.2 (2.76) m for canopy ceiling height, 0.015 (0.02) kg/m³ for canopy bulk density, 0.068 (0.133) kg/m² for available canopy fuel, and 0.33 (0.39) m²/m² for leaf area index. Moreover, uncertainty estimates from the GP regression were more indicative of the true errors in the predicted canopy variables than those from their linear counterparts. With the ever-increasing accessibility of multisource remote sensing data, we envision a concomitant expansion in the use of advanced statistical methods, such as machine learning, to explore the potentially complex relationships between canopy characteristics and remotely-sensed predictors, accompanied by a desideratum for improved error analysis.     

Geographic variability in lidar predictions of forest stand structure in the Pacific Northwest

Estimation of the amount of carbon stored in forests is a key challenge for understanding the global carbon cycle, one which remote sensing is expected to help address. However, carbon storage in moderate to high biomass forests is difficult to estimate with conventional optical or radar sensors. Lidar (light detection and ranging) instruments measure the vertical structure of forests and thus hold great promise for remotely sensing the quantity and spatial organization of forest biomass. In this study, we compare the relationships between lidar-measured canopy structure and coincident field measurements of forest stand structure at five locations in the Pacific Northwest of the U.S.A. with contrasting composition. Coefficient of determination values (r²) ranged between 41% and 96%. Correlations for two important variables, LAI (81%) and aboveground biomass (92%), were noteworthy, as was the fact that neither variable showed an asymptotic response.Of the 17 stand structure variables considered in this study, we were able to develop eight equations that were valid for all sites, including equations for two variables generally considered to be highly important (aboveground biomass and leaf area index). The other six equations that were valid for all sites were either related to height (which is most directly measured by lidar) or diameter at breast height (which should be closely related to height). Four additional equations (a total of 12) were applicable to all sites where either Douglas-fir (Pseudotsuga menziesii), western hemlock (Tsuga heterophylla) or Sitka spruce (Picea sitchensi) were dominant. Stand structure variables in sites dominated by true firs (Abies sp.) or ponderosa pine (Pinus ponderosa) had biases when predicted by these four additional equations. Productivity-related variables describing the edaphic, climatic and topographic environment of the sites where available for every regression, but only two of the 17 equations (maximum diameter at breast height, stem density) incorporated them. Given the wide range of these environmental conditions sampled, we conclude that the prediction of stand structure is largely independent of environmental conditions in this study area.Most studies of lidar remote sensing for predicting stand structure have depended on intensive data collections within a relatively small study area. This study indicates that the relationships between many stand structure indices and lidar measured canopy structure have generality at the regional scale. This finding, if replicated in other regions, would suggest that mapping of stand structure using lidar may be accomplished by distributing field sites extensively over a region, thus reducing the overall inventory effort required.     

Development of a national model of Pinus radiata stand volume from lidar metrics for New Zealand

Light detection and ranging (lidar) has been successfully used to describe a wide range of forest metrics at local scales. However, little research has tested the general applicability of this technology to describe commercially important stand dimensions, such as total stem volume (V), at national levels across broad environmental gradients.

Using an extensive national data set covering the spatial extent of Pinus radiata plantation forests in New Zealand, the key objectives of this study were to (1) develop regression models to best describe V for P. radiata from lidar metrics and (2) investigate whether these relationships could be improved using coincident environmental and stand-level information. Development of relationships between lidar metrics and forest volume are of particular importance for P. radiata, as this species constitutes approximately 90% of the 1.8 Mha plantation resource.

Using lidar mean height and the percentage of lidar ground returns, the initial model (model 1) accounted for 85% of the variance in V. Addition of stand stocking (number of stems ha^?1), measured within the plots, to the model (model 2) significantly (p < 0.001) improved predictions, with R ² increasing to 0.86 and the root mean square error declining from 80.1 m³ ha^?1 to 71.6 m³ ha^?1. For both models, partial responses show V to be most sensitive to lidar mean height, which was included in the model as a second-order polynomial.

Although environmental variables are established determinants for V, their inclusion did not significantly improve either model 1 or 2. Residual values for both models showed little apparent bias when plotted against stand-level information or a wide array of environmental variables, supporting the general applicability of these relationships.     

Predicting Eucalyptus spp. stand volume in Zululand,South Africa: an analysis using a stochastic gradient boosting regression ensemble with multi-source data sets

Accurate, reliable, and up-to-date forest stand volume information is a prerequisite for a detailed evaluation of commercial forest resources and their sustainable management. Commercial forest responses to global climate change remain uncertain, and hence the mapping of stand volume as carbon sinks is fundamentally important in understanding the role of forests in stabilizing climate change effects. The aim of this study was to examine the utility of stochastic gradient boosting (SGB) and multi-source data to predict stand volume of a Eucalyptus plantation in South Africa. The SGB ensemble, random forest (RF), and stepwise multiple-linear regression (SMLR) were used to predict Eucalyptus stand volume and other related tree-structural attributes such as mean tree height and mean diameter at breast height (DBH). Multi-source data consisted of SPOT-5 raw spectral features (four bands), 14 spectral vegetation indices, rainfall data, and stand age. When all variables were used, the SGB algorithm showed that stand volume can be accurately estimated (R² = 0.78 and RMSE = 33.16 m³ ha^?1 (23.01% of the mean)). The competing RF ensemble produced an R² value of 0.76 and a RMSE value of 37.28 m³ ha^?1 (38.28% of the mean). SMLR on the other hand, produced an R² value of 0.65 and an RMSE value of 42.50 m³ ha^?1 (42.50% of the mean). Our study further showed that Eucalyptus mean tree height (R² = 0.83 and RMSE = 1.63 m (9.08% of the mean)) and mean diameter at breast height (R² = 0.74 and RMSE = 1.06 (7.89% of the mean)) can also be reasonably predicted using SGB and multi-source data. Furthermore, when the most important SGB model-selected variables were used for prediction, the predictive accuracies improved significantly for mean DBH (R² = 0.81 and RMSE = 1.21 cm (6.12% of the mean)), mean tree height (R² = 0.86 and RMSE = 1.39 m (7.02% of the mean)), and stand volume (R² = 0.83 and RMSE = 29.58 m³ ha^?1 (17.63% of the mean)). These results underscore the importance of integrating multi-source data with remotely sensed data for predicting Eucalyptus stand volume and related tree-structural attributes.     

Patterns of covariance between forest stand and canopy structure in the Pacific Northwest

In the past decade, lidar (light detection and ranging) has emerged as a powerful tool for remotely sensing forest canopy and stand structure, including the estimation of aboveground biomass and carbon storage. Numerous papers have documented the use of lidar measurements to predict important aspects of forest stand structure, including aboveground biomass. Other papers have documented the ability to transform lidar measurements to approximate common field measures, such as cover, stand height, and vertical distributions of foliage density and light transmittance. However, only a small number of existing works have thoroughly examined relationships between comprehensive assemblages of forest canopy and forest stand structure indices. In this work, canonical correlation analysis of coincident lidar and field datasets in western Oregon and Washington is used to define seven statistically significant pairs of canonical variables, each defining an axis of variation that stand and canopy structure have in common. The first major axis relates mean stand height, and related variables, to aboveground biomass. The second relates canopy cover and volume to leaf area index and stem density. The third relates canopy height variability to mean stem diameter and the basal area of deciduous species. Of the four remaining axes, three are related to contrasts between mature and old-growth stands. Canonical correlation analysis provides a method for ranking the importance of these effects, and for placing both canopy and stand structure indices within the overall covariance structure of the two datasets. In this sense, and for the study area involved, the first three factors (mean height, cover or leaf index area, height variability) represent the same kind of enhancement of lidar data that the tasseled cap indices [Crist, C.P., R.C. Cicone, 1984. A physically-based transformation of thematic mapper data—the TM tasseled cap. IEEE Transactions on Geoscience and Remote Sensing 22, 256-263.] represent for optical remote sensing.