Automated Structural-level Alignment of Multi-view TLS and ALS Point Clouds in Forestry

Access to highly detailed models of heterogeneous forests from the near surface to above the tree canopy at varying scales is of increasing demand as it enables more advanced computational tools for analysis, planning, and ecosystem management. LiDAR sensors available through different scanning platforms including terrestrial, mobile and aerial have become established as one of the primary technologies for forest mapping due to their inherited capability to collect direct, precise and rapid 3D information of a scene. However, their scalability to large forest areas is highly dependent upon use of effective and efficient methods of co-registration of multiple scan sources. Surprisingly, work in forestry in GPS denied areas has mostly resorted to methods of co-registration that use reference based targets (e.g., reflective, marked trees), a process far from scalable in practice. In this work, we propose an effective, targetless and fully automatic method based on an incremental co-registration strategy matching and grouping points according to levels of structural complexity. Empirical evidence shows the method's effectiveness in aligning both TLS-to-TLS and TLS-to-ALS scans under a variety of ecosystem conditions including pre/post fire treatment effects, of interest to forest inventory surveyors

Developing Allometric Equations for Teak Plantations Located in the Coastal Region of Ecuador from Terrestrial Laser Scanning Data

Traditional studies aimed at developing allometric models to estimate dry above-ground biomass (AGB) and other tree-level variables, such as tree stem commercial volume (TSCV) or tree stem volume (TSV), usually involves cutting down the trees. Although this method has low uncertainty, it is quite costly and inefficient since it requires a very time-consuming field work. In order to assist in data collection and processing, remote sensing is allowing the application of non-destructive sampling methods such as that based on terrestrial laser scanning (TLS). In this work, TLS-derived point clouds were used to digitally reconstruct the tree stem of a set of teak trees (Tectona grandis Linn. F.) from 58 circular reference plots of 18 m radius belonging to three different plantations located in the Coastal Region of Ecuador. After manually selecting the appropriate trees from the entire sample, semi-automatic data processing was performed to provide measurements of TSCV and TSV, together with estimates of AGB values at tree level. These observed values were used to develop allometric models, based on diameter at breast height (DBH), total tree height (h), or the metric DBH2 × h, by applying a robust regression method to remove likely outliers. Results showed that the developed allometric models performed reasonably well, especially those based on the metric DBH2 × h, providing low bias estimates and relative RMSE values of 21.60% and 16.41% for TSCV and TSV, respectively. Allometric models only based on tree height were derived from replacing DBH by h in the expression DBH2 x h, according to adjusted expressions depending on DBH classes (ranges of DBH). This finding can facilitate the obtaining of variables such as AGB (carbon stock) and commercial volume of wood over teak plantations in the Coastal Region of Ecuador from only knowing the tree height, constituting a promising method to address large-scale teak plantations monitoring from the canopy height models derived from digital aerial stereophotogrammetry

Tree extraction and estimation of walnut structure parameters using airborne LiDAR data

[EN] The development of new tools based on remote sensing data in agriculture contributes to cost reduction, increased production, and greater profitability. Airborne LiDAR (Light Detection and Ranging) data show a significant potential for geometrically characterizing tree plantations. This study aims to develop a methodology to extract walnut (Juglans regia L.) crowns under leafless conditions using airborne LiDAR data. An original approach based on the alpha-shape algorithm, identification of local maxima, and k-means algorithms is developed to extract the crowns of walnut trees in a plot located in Viver (Eastern Spain) with 192 trees. In addition, stem diameter and volume, crown diameter, total height, and crown height were estimated from cloud metrics and other 2D parameters such as crown area, and diameter derived from LiDAR data. A correct identification was made of 178 trees (92.7%). For structure parameters, the most accurate results were obtained for crown diameter, stem diameter, and stem volume with coefficient of determination values (R-2) equal to 0.95, 0.87 and 0.83; and RMSE values of 0.43 m (5.70%), 0.02 m (9.35%) and 0.016 m(3) (21.55%), respectively. The models that gave the lowest R-2 values were 0.69 for total height and 0.70 for crown height, with RMSE values of 0.84 m (12.4%) and 0.83 m (14.5%), respectively. A suitable definition of the central and lower parts of tree canopies was observed. Results of this study generate valuable information, which can be applied for improving the management of walnut plantations.Estornell Cremades, J.; Hadas, E.; Marti-Gavila, J.; López- Cortés, I. (2021). Tree extraction and estimation of walnut structure parameters using airborne LiDAR data. International Journal of Applied Earth Observation and Geoinformation. 96:1-9. https://doi.org/10.1016/j.jag.2020.102273S199

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

Multisensorial Close-Range Sensing Generates Benefits for Characterization of Managed Scots Pine (Pinus sylvestris L.) Stands

Terrestrial laser scanning (TLS) provides a detailed three-dimensional representation of surrounding forest structures. However, due to close-range hemispherical scanning geometry, the ability of TLS technique to comprehensively characterize all trees, and especially upper parts of forest canopy, is often limited. In this study, we investigated how much forest characterization capacity can be improved in managed Scots pine (Pinus sylvestris L.) stands if TLS point clouds are complemented with photogrammetric point clouds acquired from above the canopy using unmanned aerial vehicle (UAV). In this multisensorial (TLS+UAV) close-range sensing approach, the used UAV point cloud data were considered especially suitable for characterizing the vertical forest structure and improvements were obtained in estimation accuracy of tree height as well as plot-level basal-area weighted mean height (Hg) and mean stem volume (Vmean). Most notably, the root-mean-square-error (RMSE) in Hg improved from 0.8 to 0.58 m and the bias improved from −0.75 to −0.45 m with the multisensorial close-range sensing approach. However, in managed Scots pine stands, the mere TLS also captured the upper parts of the forest canopy rather well. Both approaches were capable of deriving stem number, basal area, Vmean, Hg, and basal area-weighted mean diameter with the relative RMSE less than 5.5% for all the sample plots. Although the multisensorial close-range sensing approach mainly enhanced the characterization of the forest vertical structure in single-species, single-layer forest conditions, representation of more complex forest structures may benefit more from point clouds collected with sensors of different measurement geometries

Multisensorial Close-Range Sensing Generates Benefits for Characterization of Managed Scots Pine (Pinus sylvestris L.) Stands

Terrestrial laser scanning (TLS) provides a detailed three-dimensional representation of surrounding forest structures. However, due to close-range hemispherical scanning geometry, the ability of TLS technique to comprehensively characterize all trees, and especially upper parts of forest canopy, is often limited. In this study, we investigated how much forest characterization capacity can be improved in managed Scots pine (Pinus sylvestris L.) stands if TLS point clouds are complemented with photogrammetric point clouds acquired from above the canopy using unmanned aerial vehicle (UAV). In this multisensorial (TLS+UAV) close-range sensing approach, the used UAV point cloud data were considered especially suitable for characterizing the vertical forest structure and improvements were obtained in estimation accuracy of tree height as well as plot-level basal-area weighted mean height (Hg) and mean stem volume (Vmean). Most notably, the root-mean-square-error (RMSE) in Hg improved from 0.8 to 0.58 m and the bias improved from −0.75 to −0.45 m with the multisensorial close-range sensing approach. However, in managed Scots pine stands, the mere TLS also captured the upper parts of the forest canopy rather well. Both approaches were capable of deriving stem number, basal area, Vmean, Hg, and basal area-weighted mean diameter with the relative RMSE less than 5.5% for all the sample plots. Although the multisensorial close-range sensing approach mainly enhanced the characterization of the forest vertical structure in single-species, single-layer forest conditions, representation of more complex forest structures may benefit more from point clouds collected with sensors of different measurement geometries

Forest structure from terrestrial laser scanning – in support of remote sensing calibration/validation and operational inventory

Forests are an important part of the natural ecosystem, providing resources such as timber and fuel, performing services such as energy exchange and carbon storage, and presenting risks, such as fire damage and invasive species impacts. Improved characterization of forest structural attributes is desirable, as it could improve our understanding and management of these natural resources. However, the traditional, systematic collection of forest information – dubbed “forest inventory” – is time-consuming, expensive, and coarse when compared to novel 3-D measurement technologies. Remote sensing estimates, on the other hand, provide synoptic coverage, but often fail to capture the fine- scale structural variation of the forest environment. Terrestrial laser scanning (TLS) has demonstrated a potential to address these limitations, but its operational use has remained limited due to unsatisfactory performance characteristics vs. budgetary constraints of many end-users. To address this gap, my dissertation advanced affordable mobile laser scanning capabilities for operational forest structure assessment. We developed geometric reconstruction of forest structure from rapid-scan, low-resolution point cloud data, providing for automatic extraction of standard forest inventory metrics. To augment these results over larger areas, we designed a view-invariant feature descriptor to enable marker-free registration of TLS data pairs, without knowledge of the initial sensor pose. Finally, a graph-theory framework was integrated to perform multi-view registration between a network of disconnected scans, which provided improved assessment of forest inventory variables. This work addresses a major limitation related to the inability of TLS to assess forest structure at an operational scale, and may facilitate improved understanding of the phenomenology of airborne sensing systems, by providing fine-scale reference data with which to interpret the active or passive electromagnetic radiation interactions with forest structure. Outputs are being utilized to provide antecedent science data for NASA’s HyspIRI mission and to support the National Ecological Observatory Network’s (NEON) long-term environmental monitoring initiatives

Assessment of handheld mobile terrestrial laser scanning for estimating tree parameters

Sustainable forest management heavily relies on the accurate estimation of tree parameters. Among others, the diameter at breast height (DBH) is important for extracting the volume and mass of an individual tree. For systematically estimating the volume of entire plots, airborne laser scanning (ALS) data are used. The estimation model is frequently calibrated using manual DBH measurements or static terrestrial laser scans (STLS) of sample plots. Although reliable, this method is time-consuming, which greatly hampers its use. Here, a handheld mobile terrestrial laser scanning (HMTLS) was demonstrated to be a useful alternative technique to precisely and efficiently calculate DBH. Different data acquisition techniques were applied at a sample plot, then the resulting parameters were comparatively analysed. The calculated DBH values were comparable to the manual measurements for HMTLS, STLS, and ALS data sets. Given the comparability of the extracted parameters, with a reduced point density of HTMLS compared to STLS data, and the reasonable increase of performance, with a reduction of acquisition time with a factor of 5 compared to conventional STLS techniques and a factor of 3 compared to manual measurements, HMTLS is considered a useful alternative technique

Automatic Registration of Optical Aerial Imagery to a LiDAR Point Cloud for Generation of City Models

This paper presents a framework for automatic registration of both the optical and 3D structural information extracted from oblique aerial imagery to a Light Detection and Ranging (LiDAR) point cloud without prior knowledge of an initial alignment. The framework employs a coarse to fine strategy in the estimation of the registration parameters. First, a dense 3D point cloud and the associated relative camera parameters are extracted from the optical aerial imagery using a state-of-the-art 3D reconstruction algorithm. Next, a digital surface model (DSM) is generated from both the LiDAR and the optical imagery-derived point clouds. Coarse registration parameters are then computed from salient features extracted from the LiDAR and optical imagery-derived DSMs. The registration parameters are further refined using the iterative closest point (ICP) algorithm to minimize global error between the registered point clouds. The novelty of the proposed approach is in the computation of salient features from the DSMs, and the selection of matching salient features using geometric invariants coupled with Normalized Cross Correlation (NCC) match validation. The feature extraction and matching process enables the automatic estimation of the coarse registration parameters required for initializing the fine registration process. The registration framework is tested on a simulated scene and aerial datasets acquired in real urban environments. Results demonstrates the robustness of the framework for registering optical and 3D structural information extracted from aerial imagery to a LiDAR point cloud, when co-existing initial registration parameters are unavailable