On the use of rapid-scan, low point density terrestrial laser scanning (TLS) for structural assessment of complex forest environments

Forests fulfill an important role in natural ecosystems, e.g., they provide food, fiber, habitat, and biodiversity, all of which contribute to stable ecosystems. Assessing and modeling the structure and characteristics in forests can lead to a better understanding and management of these resources. Traditional methods for collecting forest traits, known as “forest inventory”, is achieved using rough proxies, such as stem diameter, tree height, and foliar coverage; such parameters are limited in their ability to capture fine-scale structural variation in forest environments. It is in this context that terrestrial laser scanning (TLS) has come to the fore as a tool for addressing the limitations of traditional forest structure evaluation methods. However, there is a need for improving TLS data processing methods. In this work, we developed algorithms to assess the structure of complex forest environments – defined by their stem density, intricate root and stem structures, uneven-aged nature, and variable understory - using data collected by a low-cost, portable TLS system, the Compact Biomass Lidar (CBL). The objectives of this work are listed as follow: 1. Assess the utility of terrestrial lidar scanning (TLS) to accurately map elevation changes (sediment accretion rates) in mangrove forest; 2. Evaluate forest structural attributes, e.g., stems and roots, in complex forest environments toward biophysical characterization of such forests; and 3. Assess canopy-level structural traits (leaf area index; leaf area density) in complex forest environments to estimate biomass in rapidly changing environments. The low-cost system used in this research provides lower-resolution data, in terms of scan angular resolution and resulting point density, when compared to higher-cost commercial systems. As a result, the algorithms developed for evaluating the data collected by such systems should be robust to issues caused by low-resolution 3D point cloud data. The data used in various parts of this work were collected from three mangrove forests on the western Pacific island of Pohnpei in the Federated States of Micronesia, as well as tropical forests in Hawai’i, USA. Mangrove forests underscore the economy of this region, where more than half of the annual household income is derived from these forests. However, these mangrove forests are endangered by sea level rise, which necessitates an evaluation of the resilience of mangrove forests to climate change in order to better protect and manage these ecosystems. This includes the preservation of positive sediment accretion rates, and stimulating the process of root growth, sedimentation, and peat development, all of which are influenced by the forest floor elevation, relative to sea level. Currently, accretion rates are measured using surface elevation tables (SETs), which are posts permanently placed in mangrove sediments. The forest floor is measured annually with respect to the height of the SETs to evaluate changes in elevation (Cahoon et al. 2002). In this work, we evaluated the ability of the CBL system for measuring such elevation changes, to address objective #1. Digital Elevation Models (DEMs) were produced for plots, based on the point cloud resulted from co-registering eight scans, spaced 45 degree, per plot. DEMs are refined and produced using Cloth Simulation Filtering (CSF) and kriging interpolation. CSF was used because it minimizes the user input parameters, and kriging was chosen for this study due its consideration of the overall spatial arrangement of the points using semivariogram analysis, which results in a more robust model. The average consistency of the TLS-derived elevation change was 72%, with and RMSE value of 1.36 mm. However, what truly makes the TLS method more tenable, is the lower standard error (SE) values when compared to manual methods (10-70x lower). In order to achieve our second objective, we assessed structural characteristics of the above-mentioned mangrove forest and also for tropical forests in Hawaii, collected with the same CBL scanner. The same eight scans per plot (20 plots) were co-registered using pairwise registration and the Iterative Closest Point (ICP). We then removed the higher canopy using a normal change rate assessment algorithm. We used a combination of geometric classification techniques, based on the angular orientation of the planes fitted to points (facets), and machine learning 3D segmentation algorithms to detect tree stems and above-ground roots. Mangrove forests are complex forest environments, containing above-ground root mass, which can create confusion for both ground detection and structural assessment algorithms. As a result, we needed to train a supporting classifier on the roots to detect which root lidar returns were classified as stems. The accuracy and precision values for this classifier were assessed via manual investigation of the classification results in all 20 plots. The accuracy and precision for stem classification were found to be 82% and 77%, respectively. The same values for root detection were 76% and 68%, respectively. We simulated the stems using alpha shapes in order to assess their volume in the final step. The consistency of the volume evaluation was found to be 85%. This was obtained by comparing the mean stem volume (m3/ha) from field data and the TLS data in each plot. The reported accuracy is the average value for all 20 plots. Additionally, we compared the diameter-at-breast-height (DBH), recorded in the field, with the TLS-derived DBH to obtain a direct measure of the precision of our stem models. DBH evaluation resulted in an accuracy of 74% and RMSE equaled 7.52 cm. This approach can be used for automatic stem detection and structural assessment in a complex forest environment, and could contribute to biomass assessment in these rapidly changing environments. These stem and root structural assessment efforts were complemented by efforts to estimate canopy-level structural attributes of the tropical Hawai’i forest environment; we specifically estimated the leaf area index (LAI), by implementing a density-based approach. 242 scans were collected using the portable low-cost TLS (CBL), in a Hawaii Volcano National Park (HAVO) flux tower site. LAI was measured for all the plots in the site, using an AccuPAR LP-80 Instrument. The first step in this work involved detection of the higher canopy, using normal change rate assessment. After segmenting the higher canopy from the lidar point clouds, we needed to measure Leaf Area Density (LAD), using a voxel-based approach. We divided the canopy point cloud into five layers in the Z direction, after which each of these five layers were divided into voxels in the X direction. The sizes of these voxels were constrained based on interquartile analysis and the number of points in each voxel. We hypothesized that the power returned to the lidar system from woody materials, like branches, exceeds that from leaves, due to the liquid water absorption of the leaves and higher reflectivity for woody material at the 905 nm lidar wavelength. We evaluated leafy and woody materials using images from projected point clouds and determined the density of these regions to support our hypothesis. The density of points in a 3D grid size of 0.1 m, which was determined by investigating the size of the branches in the lower portion of the higher canopy, was calculated in each of the voxels. Note that “density” in this work is defined as the total number of points per grid cell, divided by the volume of that cell. Subsequently, we fitted a kernel density estimator to these values. The threshold was set based on half of the area under the curve in each of the distributions. The grid cells with a density below the threshold were labeled as leaves, while those cells with a density above the threshold were set as non-leaves. We then modeled the LAI using the point densities derived from TLS point clouds, achieving a R2 value of 0.88. We also estimated the LAI directly from lidar data by using the point densities and calculating leaf area density (LAD), which is defined as the total one-sided leaf area per unit volume. LAI can be obtained as the sum of the LAD values in all the voxels. The accuracy of LAI estimation was found to be 90%. Since the LAI values cannot be considered spatially independent throughout all the plots in this site, we performed a semivariogram analysis on the field-measured LAI data. This analysis showed that the LAI values can be assumed to be independent in plots that are at least 30 m apart. As a result, we divided the data into six subsets, where each of the plots were 30 meter spaced for each subset. LAI model R2 values for these subsets ranged between 0.84 - 0.96. The results bode well for using this method for automatic estimation of LAI values in complex forest environments, using a low-cost, low point density, rapid-scan TLS

Advances in Waveform and Photon Counting Lidar Processing for Forest Vegetation Applications

Full waveform (FW) and photon counting LiDAR (PCL) data have garnered greater attention due to increasing data availability, a wealth of information they contain and promising prospects for large scale vegetation mapping. However, many factors such as complex processing steps and scarce non-proprietary tools preclude extensive and practical uses of these data for vegetation characterization. Therefore, the overall goal of this study is to develop algorithms to process FW and PCL data and to explore their potential in real-world applications. Study I explored classical waveform decomposition methods such as the Gaussian decomposition, Richardson–Lucy (RL) deconvolution and a newly introduced optimized Gold deconvolution to process FW LiDAR data. Results demonstrated the advantages of the deconvolution and decomposition method, and the three approaches generated satisfactory results, while the best performances varied when different criteria were used. Built upon Study I, Study II applied the Bayesian non-linear modeling concepts for waveform decomposition and quantified the propagation of error and uncertainty along the processing steps. The performance evaluation and uncertainty analysis at the parameter, derived point cloud and surface model levels showed that the Bayesian decomposition could enhance the credibility of decomposition results in a probabilistic sense to capture the true error of estimates and trace the uncertainty propagation along the processing steps. In study III, we exploited FW LiDAR data to classify tree species through integrating machine learning methods (the Random forests (RF) and Conditional inference forests (CF)) and Bayesian inference method. Results of classification accuracy highlighted that the Bayesian method was a superior alternative to machine learning methods, and rendered users with more confidence for interpreting and applying classification results to real-world tasks such as forest inventory. Study IV focused on developing a framework to derive terrain elevation and vegetation canopy height from test-bed sensor data and to pre-validate the capacity of the upcoming Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) mission. The methodology developed in this study illustrates plausible ways of processing the data that are structurally similar to expected ICESat-2 data and holds the potential to be a benchmark for further method adjustment once genuine ICESat-2 are available

Full waveform analysis for long-range 3D imaging laser radar

The new generation of 3D imaging systems based on laser radar (ladar) offers significant advantages in defense and security applications. In particular, it is possible to retrieve 3D shape information directly from the scene and separate a target from background or foreground clutter by extracting a narrow depth range from the field of view by range gating, either in the sensor or by postprocessing. We discuss and demonstrate the applicability of full-waveform ladar to produce multilayer 3D imagery, in which each pixel produces a complex temporal response that describes the scene structure. Such complexity caused by multiple and distributed reflection arises in many relevant scenarios, for example in viewing partially occluded targets, through semitransparent materials (e.g., windows) and through distributed reflective media such as foliage. We demonstrate our methodology on 3D image data acquired by a scanning time-of-flight system, developed in our own laboratories, which uses the time-correlated single-photon counting technique

Sub-canopy terrain modelling for archaeological prospecting in forested areas through multiple-echo discrete-pulse laser ranging: a case study from Chopwell Wood, Tyne & Wear

Airborne Light Detection and Ranging (LiDAR) technology is assessed for its effectiveness as a tool for measuring terrain under forest canopy. To evaluate the capability of multiple-return discrete-pulse airborne laser ranging for detecting and resolving sub-canopy archaeological features, LiDAR data were collected from a helicopter over a forest near Gateshead in July 2009. Coal mining and timber felling have characterised Chopwell Wood, a mixed coniferous and deciduous woodland of 360 hectares, since the Industrial Revolution. The state-of-the-art Optech ALTM 3100EA LiDAR system operated at 70,000 pulses per second and raw data were acquired over the study area at a point density of over 30 points per square metre. Reference terrain elevation data were acquired on-site to ‘train’ the progressive densification filtering algorithm of Axelsson (1999; 2000) to identify laser reflections from the terrain surface. A number of sites, offering a variety of tree species, variable terrain roughness & gradient and understorey vegetation cover of varying density, were identified in the wood to assess the accuracy of filtered LiDAR terrain data. Results showed that the laser scanner over-estimated the elevation of reference terrain data by 13±17 cm under deciduous canopy and 23±18 cm under coniferous canopy. Terrain point density was calculated as 4.1 and 2.4 points per square metre under deciduous and coniferous forest, respectively. Classified terrain points were modelled with the kriging interpolation technique and topographic archaeological features, such as coal tubways (transportation routes) and areas of subsidence over relic mine shafts, were identified in digital terrain models (DTMs) using advanced exaggeration and artificial illumination techniques. Airborne LiDAR is capable of recording high quality terrain data even under the most dense forest canopy, but the accuracy and density of terrain data are controlled by a combination of tree species, forest management practices and understorey vegetation

Individual tree detection and modelling aboveground biomass and forest parameters using discrete return airborne LiDAR data

Individual tree detection and modelling forest parameters using Airborne Laser Scanner data (Light Detection and Ranging (LiDAR) is becoming increasingly important for the monitoring and sustainable management of forests. Remote sensing has been a useful tool for individual tree analysis in the past decade, although inadequate spatial resolution from satellites means that only airborne systems have sufficient spatial resolution to conduct individual tree analysis. Moreover, recent advances in airborne LiDAR now provide high horizontal resolution as well as information in the vertical dimension. However, it is challenging to fully exploit and utilize small-footprint LiDAR data for detailed tree analysis. Procedures for forest biomass quantification and forest attributes measurement using LiDAR data have improved at a rapid pace as more robust and sophisticated modelling used to improve the studies. This thesis contains an evaluation of three approaches of utilizing LiDAR data for individual tree forest measurement. The first explores the relationship between LiDAR metrics and field reference to assess the correlation between LiDAR and field data at the individual-tree level. The intention was not to detect trees automatically, but to develop a LiDAR-AGB model based on trees that were mapped in the field so as to evaluate the relationships between LiDAR-type metrics under controlled conditions for the study sites, and field-derived AGB. A non-linear AGB model based on field data and LiDAR data was developed and LiDAR height percentile h80 and crown width measurement (CW) was found to best fit the data as evidenced by and Adj-R2 value of 0.63, the root mean squared error of the model of 14.8% and analysis of the residuals. This paper provides the foundation for a predictive LiDAR-AGB model at tree level over two study sites, Pasoh Forest Reserve and FRIM Forest Reserve. The second part of the thesis then takes this AGB-LiDAR relationship and combines it with individual tree crown delineation. This chapter shows the contribution of performing an automatic individual tree crown delineation over the wider forest areas. The individual tree crown delineation is composed of a five-step framework, which is unique in its automated determination of dominant crown sizes in a forest area and its adaption of the LiDAR-AGB model developed for the purpose of validation the method. This framework correctly delineated 84% and 88% of the tree crowns in the two forest study areas which is mostly dominated with lowland dipterocarp trees. Thirdly, parametric and non-parametric modelling approaches are proposed for modelling forest structural attributes. Selected modelling methods are compared for predicting 4 forest attributes, volume (V), basal area (BA), height (Ht) and aboveground biomass (AGB) at the species level. The AGB modelling in this paper is extracted using the LiDAR derived variables from the automated individual tree crown delineation, in contrast to the earlier AGB modelling where it is derived based on the trees that were mapped in the field. The selected non-parametric method included, k-nearest neighbour (k-NN) imputation methods: Most Similar Neighbour (MSN) and Gradient Nearest Neighbour (GNN), Random Forest (RF) and parametric approach: Ordinary Least Square (OLS) regression. To compare and evaluate these approaches a scaled root mean squared error (RMSE) between observed and predicted forest attribute sampled from both forest site was computed. The best method varied according to response variable and performance measure. OLS regression was to found to be the best performance method overall evidenced by RMSE after cross validation for BA (1.40 m2), V (1.03 m3), Ht (2.22 m) and AGB (96 Kg/tree) respectively, showed its applicability to wider conditions, while RF produced best overall results among the non-parametric methods tested. This thesis concludes with a discussion of the potential of LiDAR data as an independent source of important forest inventory data source when combined with appropriate designed sample plots in the field, and with appropriate modelling tools