The automatic tree detection and delineation from Airborne LiDAR

Individual tree detection from airborne laser scanning data based on supervoxels and local convexity

The rapid development of LiDAR technology has promoted great changes in forest resource surveys. The airborne LiDAR point cloud can provide the precise height and detailed structure of trees, and can estimate key forest resource indicators such as forest stock volume, diameter at breast height, and forest biomass at a large scale. By establishing relationship models between the forest parameters of sample plots and the calculated parameters of LiDAR, these developments may eventually expand the models to large-scale forest resource surveys of entire areas. In this study, eight sample plots in northeast China are used to verify and update the information using point cloud obtained by the LiDAR scanner riegl-vq-1560i. First, the tree crowns are segmented using the profile-rotating algorithm, and tree positions are registered based on dominant tree heights. Second, considering the correlation between crown shape and tree species, we use DBN classifier to identify species using features of crowns, which are extracted to 1D array. Third, when the tree species is known, parameters such as height, crown width, diameter at breast height, biomass, and stock volume can be extracted from trees, enabling accurate large-scale forest surveys based on LiDAR data. Finally, experiment results demonstrate that the F-score of the eight plots in the tree segmentation exceed 0.95, the accuracy of tree species correction exceeds 90%, and the R2 of tree height, east–west crown width, north–south crown width, diameter at breast height, aboveground biomass, and stock volume are 0.893, 0.757, 0.694, 0.840, 0.896 and 0.891, respectively. The above results indicate that the LiDAR-based survey is practical and can be widely applied in forest resource monitoring.

This paper provides a practical procedure for fusing LiDAR and photogrammetric point clouds for the extraction of tree metrics. Aerial photogrammetric point clouds are first generated using the structure-from-motion and dense-matching methods. Registration of the LiDAR and photogrammetric point clouds is then performed using an onboard global positioning system and inertial measurement unit. However, due to systematic deviations, it is impossible to directly merge the two types of point cloud. Therefore, an urban street network obtained from the OpenStreetMap digital mapping system is utilized for point cloud segmentation. After segmentation, each chunk is finely registered and merged based on the iterative closest point algorithm, allowing the two types of data to be accurately co-registered and a fused point cloud obtained. Finally, we conducted experiments to extract stand and individual tree metrics from fused point clouds created for two study plots. The height distributions of the fused point clouds were highly consistent with LiDAR data, with the 5%, 10%, 25%, 50%, 75%, 90%, and 95% height percentiles showing acceptable similarities. The height distribution of individual trees was also consistent with that of field measurements. Furthermore, the fused point clouds contain a high point density and RGB color information, which allow shape delineation and estimation of tree health status. This comprehensive analysis demonstrates that this procedure provides a practical way to inventory tree stands and individuals in urban areas.

Trees are an integral component of the urban environment and important for human well-being, adaption measures to climate change and sustainable urban transformation. Understanding the small-scale impacts of urban trees and strategically managing the ecosystem services they provide requires high-resolution information on urban forest structure, which is still scarce. In contrast, there is an abundance of data portraying urban areas and an associated trend towards smart cities and digital twins as analysis platforms. A GIS workflow is presented in this paper that may close this data gap by classifying the urban forest from LiDAR point clouds, detecting and reconstructing individual crowns, and enabling a tree representation within semantic 3D city models. The workflow is designed to provide robust results for point clouds with a density of at least 4 pts/m2 that are widely available. Evaluation was conducted by mapping the urban forest of Dresden (Germany) using a point cloud with 4 pts/m². An object-based data fusion approach is implemented for the classification of the urban forest. A classification accuracy of 95 % for different urban settings is achieved by combining LiDAR with multispectral imagery and a 3D building model. Individual trees are detected by local maxima filtering and crowns are segmented using marker-controlled watershed segmentation. Evaluation highlights the influences of both urban and forest structure on individual tree detection. Substantial differences in detection accuracies are evident between trees along streets (72 %) and structurally more complex tree stands in green areas (31 %), as well as dependencies on tree height and crown diameter. Furthermore, an approach for parameterized reconstruction of tree crowns is presented, which enables efficient and realistic city-wide modeling. The suitability of LiDAR to measure individual tree metrics is illustrated as well as a framework for modeling individual tree crowns via geometric primitives.

Identifying bare-earth or ground returns within point cloud data is a crucially important process for archaeologists who use airborne LiDAR data, yet there has thus far been very little comparative assessment of the available archaeology-specific methods and their usefulness for archaeological applications. This article aims to provide an archaeology-specific comparison of filters for ground extraction from airborne LiDAR point clouds. The qualitative and quantitative comparison of the data from four archaeological sites from Austria, Slovenia, and Spain should also be relevant to other disciplines that use visualized airborne LiDAR data. We have compared nine filters implemented in free or low-cost off-the-shelf software, six of which are evaluated in this way for the first time. The results of the qualitative and quantitative comparison are not directly analogous, and no filter is outstanding compared to the others. However, the results are directly transferable to real-world problem-solving: Which filter works best for a given combination of data density, landscape type, and type of archaeological features? In general, progressive TIN (software: lasground_new) and a hybrid (software: Global Mapper) commercial filter are consistently among the best, followed by an open source slope-based one (software: Whitebox GAT). The ability of the free multiscale curvature classification filter (software: MCC-LIDAR) to remove vegetation is also commendable. Notably, our findings show that filters based on an older generation of algorithms consistently outperform newer filtering techniques. This is a reminder of the indirect path from publishing an algorithm to filter implementation in software.

Vertical stratification of forest canopy for segmentation of understory trees within small-footprint airborne LiDAR point clouds

Detailed forest inventory and mensuration of individual trees have drawn attention of research society mainly to support sustainable forest management. This study aims at estimating individual tree attributes from high density point cloud obtained by terrestrial laser scanner (TLS). The point clouds were obtained over single reference tree and group of trees in forest area. The reference tree is treated as benchmark since detailed measurements of branch diameter were made on selected branches with different sizes and locations. Diameter at breast height (DBH) was measured for trees in forest. Furthermore tree height, height to crown base, crown volume and tree branch volume were also estimated for each tree. Branch diameter is estimated directly from the point clouds based on semi-automatic approach of model fitting i.e. sphere, ellipse and cylinder. Tree branch volume is estimated based on the volume of the fitted models. Tree height and height to crown base are computed using histogram analysis of the point clouds elevation. Tree crown volume is estimated by fitting a convex-hull on the tree crown. The results show that the Root Mean Squared Error (RMSE) of the estimated tree branch diameter does not have a specific trend with branch sizes and number of points used for fitting process. This explains complicated distribution of point clouds over the branches. Overall cylinder model produces good results with most branch sizes and number of point clouds for fitting. The cylinder fitting approach shows significantly better estimation results compared to sphere and ellipse fitting models.

This study examines the distribution of laser pulse returns obtained from coincident airborne and terrestrial lidar surveys of a closed-canopy red pine (Pinus resinosa) plantation. The purpose of this study is to improve our understanding of laser pulse sampling within closed canopies so that estimates of forest structural variables (e.g., biomass, needle-leaf area, and base-of-live-crown) can be improved at the individual tree and stand levels using lidar. The results of this study indicate the following: (1) There is a statistically significant difference between field measurements of tree height and estimates derived from the maximum laser pulse return from airborne and terrestrial lidar. In both cases, maximum laser pulse returns underestimate tree height by 1 m, on average. (2) Both terrestrial and airborne lidar are unable to discern the base of the measured live crown. Laser pulse returns from airborne lidar are biased towards the top of the tree crown, i.e., lowest laser pulse returns occur 1.4 m on average higher in the canopy than the measured base-of-live-crown. On the other hand, terrestrial lidar captures dieback at the base of the live crown, thereby lowering the base-of-live-crown estimate by 6.6 m, on average. (3) Median airborne laser pulse returns within the canopy (20.4 m), believed to be associated with needle leaf area, occur below the maximum frequency of laser pulse returns (20.8 m) but higher in the canopy than the height of maximum crown diameter obtained from terrestrial lidar (18.0 m). The bias of airborne laser pulse reflections towards the top of the canopy with less penetration to a depth where the maximum crown diameter occurs may result in an underestimation of the needle leaf area. The results of this research suggest that future research should focus on improving our understanding of how laser pulse returns are "triggered" within vegetated environments and how canopy properties or data acquisition parameters may influence the location of this "trigger" event.

Simulating realistic three‐dimensional (3D) forest scenes is useful in understanding the links between forest structure and ecosystem functions (e.g. radiative transfer). Light detection and ranging (LiDAR) technology provides useful 3D data for forest reconstructions since it can characterise 3D structures of individual trees and canopies. High‐density terrestrial LiDAR (terrestrial laser scanning, TLS) is suitable for fine‐scale reconstructions but is limited to smaller forest plots; low‐density airborne LiDAR (airborne laser scanning, ALS) can cover larger areas but is only suitable for coarse‐scale reconstructions. How to take advantage of TLS and ALS to enable fine‐scale forest simulations in large areas needs to be studied. We propose a new voxel‐based method for forest simulations using the integration of TLS and ALS data. TLS data of representative reference trees are used to approximate the detailed architectures of the whole forest scene, with structural information on each individual tree extracted from ALS data. The high‐density point cloud data derived from TLS and ALS data are voxelised using high resolution solid voxels for scene representation. We tested the proposed method using two virtual forests (108 m × 108 m) and a real forest (300 m × 300 m) with conifer and broadleaf species. The physically based ray tracer (PBRT) was used to visualise the true virtual forest scenes, whereas voxel‐based radiative transfer (VBRT) was used to visualise the modelled forest scenes from LiDAR data. For the real forest scene, simulated and real ALS data were compared. Our results demonstrate that the images simulated by VBRT and PBRT are similar in the virtual forest scenes, with average radiance values of 1.02 and 1.72, respectively. In the real forest scene, the distributions of points and individual tree attributes (tree height, crown radius, and tree volume) derived from real and simulated ALS also match well, with Kullback–Leibler divergence ranging from 0.006 to 0.06. We conclude that the new method is capable of modelling fine‐scale 3D forests in large areas (over 1 ha) when TLS and ALS data are available, and it has good potential in studying the process of radiative transfer in conifer and broadleaf forests.

More than 50% of the national lands in Japan have been surveyed by airborne laser scanning (ALS) data with different point densities; and developing an effective approach to take full advantage of these ALS data for forest management has thus become an urgent topic of study. This study attempted to assess the utility of ALS data for individual tree detection and species classification in a mixed forest with a high canopy density. For comparison, two types of tree tops and tree crowns in the study area were delineated by the individual tree crown (ITC) approach using the green band of the orthophoto imagery and the digital canopy height model (DCHM) derived from the ALS data, respectively. Then, the two types of tree crowns were classified into four classes—Pinus densiflora (Pd), Chamaecyparis obtusa (Co), Larix kaempferi (Lk), and broadleaved trees (Bl)—by a crown-based classification approach using different combinations of the three orthophoto bands with intensity and slope maps as follows: RGB (red, green and blue); RGB and intensity (RGBI); RGB and slope (RGBS); and RGB, intensity and slope (RGBIS). Finally, the tree tops were annotated with species attributes from the two best-classified tree crown maps, and the number of different tree species in each compartment was counted for comparison with the field data. The results of our study suggest that the combination of RGBIS yielded greater classification accuracy than the other combinations. In the tree crown classifications delineated by the green band and DCHM data, the improvements in the overall accuracy compared to the RGB ranged from 5.7% for the RGBS to 9.0% for the RGBIS and from 8.3% for the RGBS to 11.8% for the RGBIS. The laser intensity and slope derived from the ALS data may be valuable sources of information for tree species classification, and in terms of distinguishing species for the detection of individual trees, the findings of this study demonstrate the advantages of using DCHM instead of optical data to delineate tree crowns. In conclusion, the synthesis of individual tree delineation using DCHM data and species classification using the RGBIS combination is recommended for interpreting forest resources in the study area. However, the usefulness of this approach must be verified in future studies through its application to other forests.

Detection and delineation of individual trees is important for the estimation of forest information using Light Detection and Ranging (LiDAR) point clouds. However, over-segmentation of canopy height models (CHMs) is frequently caused by subcrowns and branches, whereas under-segmentation is caused by adjacent crowns, especially in deciduous forests. To eliminate such problems, an approach for delineating individual trees based on curvature attributes derived from LiDAR data is proposed. To detect accurate treetops for reducing over-segmentation, the dominant sizes of tree crowns are first determined and used to fit the dome-shaped tree crowns in CHMs, and then the minimum curvature values of fitted CHMs at different scales are calculated for treetop detection. For delineation of individual trees, tree features such as height variation and curvature are chosen to form a feature vector as a pixel-level logistic regression classifier. Then an ensemble classifier combining multiple logistic regression classifiers is trained to analyze the similarity between neighboring points to reduce under-segmentation. Finally, segment maps at different scales are integrated to generate the final segment map. Our experimental results on 2 test areas demonstrate that the proposed method has promising applications for eliminating over-segmentation and under-segmentation in deciduous forests.

Segmentation of vegetation LiDAR point clouds is an important method for obtaining individual tree structure parameters. The current individual tree segmentation methods are mainly for airborne LiDAR point clouds, which use elevation information to form a grid map for segmentation, or use canopy vertices as seed points for clustering. Side-view LiDAR (vehicle LiDAR and hand-held LiDAR) can acquire more information about the lower layer of trees, but it is a challenge to perform the individual tree segmentation because the structure of side-view LiDAR point clouds is more complex. This paper proposes an individual tree segmentation method called Shadow-cut to extract the contours of the street tree point cloud. Firstly, we separated the region of the trees using the binary classifier (e.g., support vector machine) based on point cloud geometric features. Then, the optimal projection of the 3D point clouds to the 2D image is calculated and the optimal projection is the case where the pixels of the individual tree image overlap the least. Finally, after using the image segmentation algorithm to extract the tree edges in the 2D image, the corresponding 3D individual tree point cloud contours are matched with the pixels of individual tree edges in the 2D image. We conducted experiments with the proposed method on LiDAR data of urban street trees, and the correctness, completeness, and quality of the proposed individual tree segmentation method reached 91.67%, 85.33%, and 79.19%, which were superior to the CHM-based method by 2.70%, 6.19%, and 7.12%, respectively. The results show that this method is a practical and effective solution for individual tree segmentation in the LiDAR point clouds of street trees.

Abstract. With recent government initiatives for smart cities, 3D virtual city models are in demand for managing and monitoring the urban infrastructure. 3D building models forms an important component of 3D virtual city model. LiDAR remote sensing has revolutionized the way the third dimension can be precisely mapped and proved to be an important source of data for 3D models. The model thus generated should be in an open data format to be used across various applications. CityGML is an open data model framework that enables storage and exchange of 3D models which can be used for diversified applications. The main objective of this research is to develop a methodological workflow to create 3D building models in CityGML standard from airborne LiDAR point cloud. Initially building points were isolated from the airborne LiDAR data using point cloud processing algorithms. 3D building models with levels of detail (LoD1 and LoD2) were generated from the building points in CityGML standard using commercial (ArcGISPro) well as open source packages (3dfier, Citygml4j). Results prove that the models developed using open source packages are comparable to that provided by commercial packages.

Forest species diversity mapping using airborne LiDAR and hyperspectral data in a subtropical forest in China

[Figure: see text]