Nothing Special   »   [go: up one dir, main page]
Next Article in Journal
A Target Detection Algorithm Based on Fusing Radar with a Camera in the Presence of a Fluctuating Signal Intensity
Previous Article in Journal
Detecting Hailstorms in China from FY-4A Satellite with an Ensemble Machine Learning Model
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 16
Issue 18
10.3390/rs16183355
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cluster-Based Wood–Leaf Separation Method for Forest Plots Using Terrestrial Laser Scanning Data

by
Hao Tang
1,
Shihua Li
1,2,*,
Zhonghua Su
3 and
Ze He
1
1
School of Resources and Environment, University of Electronic Science and Technology of China, Chengdu 611731, China
2
Technology Innovation Center for Southwest Land Space Ecological Restoration and Comprehensive Renovation, Ministry of Natural Resources, Chengdu 610045, China
3
School of Computer and Software Engineering, Xihua University, Chengdu 610039, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(18), 3355; https://doi.org/10.3390/rs16183355
Submission received: 1 August 2024 / Revised: 31 August 2024 / Accepted: 3 September 2024 / Published: 10 September 2024
(This article belongs to the Section Forest Remote Sensing)
Browse Figures
Graphical abstract
"> Figure 1
<p>Schematics of the three locations in this study. The exact locations of the plots are marked in yellow. a, b, and c show the three areas in China, Australia, and Germany, respectively. (<b>a_1</b>–<b>c_1</b>) are remote sensing satellite maps. (<b>a_2</b>–<b>c_2</b>) are original data from representative point clouds at the locations.</p> "> Figure 2
<p>The workflow of this study. In this figure, the yellow boxes indicate data or the results of processing, and the blue boxes indicate processes or algorithms used during processing.</p> "> Figure 3
<p>Partial results after clustering by the region-growing clustering algorithm. Each color in the figure represents a cluster.</p> "> Figure 4
<p>Schematic representation of the cluster types remaining after clustering and initial wood separation and their Chaos Distance distributions. (<b>a</b>) contains four clusters, each of which contains wood. This type of cluster mainly consists of branches. (<b>b</b>) contains four clusters and denotes the type of cluster that contains one leaf or a few leaves. (<b>c</b>) is composed of wood, specifically showing the type of cluster that contains a trunk. The rough situation of the distribution of Chaos Distances for each type of cluster is indicated in (<b>d</b>), corresponding to the colors in the upper half. The horizontal axis represents the size of the <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">d</mi> <mi mathvariant="normal">c</mi> </msub> </mrow> </semantics></math> in meters, and the vertical axis represents the approximate probability density.</p> "> Figure 5
<p>Schematic representation of the distribution of the canopy structure in the forest. The point cloud colors in the figure are shown using height information.</p> "> Figure 6
<p>The results of separating wood and leaves using this method, with the overall separation on the left, the separation of the clipped sample plot on the right in the upper half, and the corresponding wood points in the lower half, where the sample plot was clipped for better visibility. The bottom half of (<b>a</b>) and the right halves of (<b>b</b>,<b>c</b>) contain the results of single trees and partial canopies in each plot after separation. These detailed figures show the separation situations and the wood results in the corresponding areas. (<b>a</b>) Chinese sample plot (Plot 1). (<b>b</b>) Australian sample plot (Plot 2). (<b>c</b>) German sample plot (Plot 3).</p> "> Figure 7
<p>The contribution of each part of the separation results.</p> "> Figure 8
<p>Separation of individual trees of different species in the plots. (<b>a</b>) is a tree in Plot 1, and (<b>b</b>) is in Plot 2. (<b>c</b>,<b>d</b>) are different tree species in Plot 3.</p> "> Figure 9
<p>The quantitative accuracy of our method, LeWoS, and the RF model was evaluated. (<b>a</b>–<b>c</b>) correspond to Plot 1, Plot 2, and Plot 3, respectively.</p> "> Figure 10
<p>Localized magnified views of the separation results. (<b>a</b>–<b>c</b>) correspond to Plots 1, 2, and 3, respectively. In each figure, the figures of overall views in (1) show partial separation results, (2) the red box figures show the locations with better separation, and (3) the blue box figures show the locations with errors.</p> ">
Versions Notes

Abstract

:
Successfully separating wood and leaves in forest plots is a prerequisite for measuring structural parameters and reconstructing 3D forest models. Terrestrial laser scanning (TLS) can distinguish between the leaves and wood of trees through precise and dense point clouds. However, most existing wood–leaf separation methods face significant accuracy issues, especially in dense forests, due to the complications introduced by canopy shading. In this study, we propose a method to separate the wood and leaves in forest plots using the clustering features of TLS data. The method first filters a point cloud to remove the ground points, and then clusters the point cloud using a region-growing algorithm. Next, the clusters are processed based on their sizes and numbers of points for preliminary separation. Chaos Distance is introduced to characterize the observation that wood points are more orderly while leaf points are more chaotic and disorganized. Lastly, the clusters’ Chaos Distance is used for the final separation. Three representative plots were used to validate this method, achieving an average accuracy of 0.938, a precision of 0.927, a recall of 0.892, and an F1 score of 0.907. The three sample plots were processed in 5.18, 3.75, and 14.52 min, demonstrating high efficiency. Comparing the results with the LeWoS and RF models showed that our method better addresses the accuracy issues of complex canopy structures.

Graphical Abstract">

Graphical Abstract

1. Introduction

Effective management and conservation of forests, which are crucial for climate regulation, carbon cycling, and biodiversity preservation, require a precise understanding of their structural dynamics and compositions [1,2]. Accurate separation of wood and leaves is a critical step in forest ecology studies, particularly for calculating important parameters like the Leaf Area Index (LAI), Leaf Area Density (LAD), Aboveground Biomass (AGB), and so on. These parameters are foundational for understanding forest structure and health. Additionally, accurate separation is necessary for advanced tasks such as tree skeleton extraction and 3D reconstruction, which require high-quality input data [3]. Therefore, the precision of these subsequent tasks is closely linked to the accuracy of wood–leaf separation. Traditional methods primarily rely on manual leaf extraction, which is both costly and time-consuming. With the development of passive remote sensing technology (e.g., photogrammetric techniques), the efficiency of wood–leaf separation can be improved by applying Computer Vision (CV) theory combined with digital image processing [4,5,6]. However, these techniques have poor penetration abilities and cannot detect the interior of the canopy [7,8].
Light Detection and Ranging (LiDAR) is an emerging technology that offers high-precision active remote sensing detection. It is a powerful optical sensing technology that uses a laser as an emitting light source to measure the distance and other characteristics of a target object [9,10]. It works by emitting a laser signal in the green or near-infrared band towards a target. Then, it compares the received signal reflected from the target (target echo) with the transmitted signal. Relevant information about the target can be obtained after appropriate processing [11,12]. Compared to data collected using other remote sensing methods, LiDAR data provide more realistic depictions of object surfaces, morphologies, and other relevant details [9,13]. Terrestrial Laser Scanning (TLS) is a type of LiDAR that generates denser and more detailed point clouds compared to Backpack Laser Scanning (BLS), Uncrewed Aerial Vehicle Laser Scanning (ULS), and Airborne Laser Scanning (ALS). This high level of detail makes TLS data suitable for achieving wood–leaf separation [14].
TLS point cloud data contain 3D coordinates and intensity information for each point. At present, wood–leaf separation technology is mature for individual trees. However, the separation of the wood and leaves of individual trees cannot meet the needs of practical applications most of the time. In addition, when it comes to forest plots, the processing speed is slow due to the large amounts of data. Therefore, it is important to achieve efficient wood–leaf separation in forest plots. Methods for wood–leaf separation can be categorized into two types: those using a combination of radiation and geometric features and those using only geometric features [15]. Effective wood–leaf separation can be achieved by simultaneously analyzing the intensity and 3D coordinate differences between wood and leaves [16,17,18]. Despite its theoretical accuracy, applying this method to actual forest plots faces challenges, primarily due to the sensitivity of the intensity differences to various conditions [19]. To generalize these kinds of methods to actual forest plots, they still need to address the issue of intensity variability.
Geometric features are commonly used for wood–leaf separation and can be categorized based on their foci as graph-based features, local neighborhood features, and cluster-based features [15,20]. Graph-based features include the shortest paths and the mode point of the data. Multiple studies have shown that these parameters can be used to achieve high accuracy in wood–leaf separation [20,21,22]. Hui et al. [20] proposed a separation method based on the evolution of mode points with an average Kappa of 0.771. Tian et al. [21] proposed a graph-based wood–leaf separation (GBS) method that synthesizes the advantages of several existing wood–leaf separation methods by making full use of shortest path information. The average accuracy of GBS has reached 0.94, which is significantly better than those of other mainstream methods. However, existing graph-based methods have primarily been applied to single trees, and their application to forest plots has not yet been validated. Additionally, these methods often require empirical thresholds and are complex in principle [20,21].
The method based on local neighborhood features evaluates the geometric features of each point’s neighbors to achieve wood–leaf separation. The main parameters are the normal vector difference and the local density [23,24,25]. Separation using these features combined with machine learning has also been performed [24,25,26]. However, machine learning requires a large amount of data for manual training and a significant initial investment in terms of effort and time for program design [27,28]. The method based on local neighborhood features focuses on and fully utilizes the differences in the morphological characteristics of wood and leaves. Zhu et al. [24] proposed an adaptive radius-based local neighborhood search approach for forest location clouds combined with a Random Forest (RF) model. Moorthy et al. [26] proposed a method for wood–leaf separation that combines the geometric features defined by the radial boundary of the nearest neighbors in a machine learning model. The above studies have indicated that machine learning techniques based on local neighborhood features are feasible for wood–leaf separation and achieve high accuracy. However, these methods usually lack good transferability between different tree species and are significantly inefficient. While local neighborhood-based methods are generally easier to understand than graph-based methods, they require features to be extracted for many points, making them less efficient overall [24,26].
Cluster-based feature methods categorize a point cloud into clusters with similar attributes, then designate these clusters as either wood or leaves [29]. The main clustering methods for point clouds are Euclidean Clustering (EC), density clustering (e.g., Density-Based Spatial Clustering of Applications with Noise (DBSCAN)), hypervoxel clustering, K-Means, etc. [30,31]. Tan et al. [32] used geometric quantities that were derived by setting a threshold to separate leaves from wood. Wan et al. [33] proposed a segment-wise classification method where a point cloud was segmented into multiple small segments, and each segment was classified as wood or leaf based on its geometric features. LeWoS, an open-source method for wood–leaf separation in forest plots, has the advantage that wood–leaf separation can be completed with comparatively high accuracy [3,29]. The rationale for these methods is easier to understand compared to graph-based and local neighborhood feature-based methods. However, these methods have two obvious drawbacks: (1) they ignore differences in leaf and wood morphology, focusing only on geometric features, and (2) they are sensitive to canopy shading when identifying clusters, making separation difficult in dense forests.
To address the issues of inaccuracy in canopy shading situations and inefficiency, this study proposes a new method for separating wood and leaves in forest plots based on cluster features. This method efficiently achieves wood–leaf separation while maintaining high accuracy, particularly in dense forests. This study uses curvature-based regional growth clustering to demonstrate the effectiveness of normal vectors in wood and leaf segmentation. A new cluster parameter called Chaos Distance is proposed. Combined with cluster size and the number of points, Chaos Distance facilitates wood–leaf separation in forest plots. Additionally, we compare the proposed method with the classical LeWoS algorithm and the RF model to highlight its advantages.

2. Selection of Sample Plots

In order to verify the robustness of the proposed method, we selected data from three locations for wood–leaf separation. Overall, the three locations varied significantly, primarily in terms of geographic location, canopy density, tree species, and point cloud density. In this study, a typical plot was selected for accuracy assessment in each of the three locations, representing a broad range of broadleaf forests. The locations and topographies of the three plots are shown in Figure 1. We calculated the LAI and crown cover values in the three plots using the point cloud data processing and analysis software LiDAR360 V7.0 [34]. A summary of the sample plots is shown in Table 1.
The sample from Plot 1 was collected by the authors at the University of Electronic Science and Technology of China (UESTC) in Chengdu City, Sichuan Province, China (Figure 1a). The plot was scanned multiple times using a Leica ScanStation C10 scanner made in Heerbrugg, Switzerland. It operates with a green laser (532 nm) and a minimum spot spacing of <1 mm. The predominant tree species in this plot was Phoebe Zhennan S. K. Lee. The trees had large numbers of branchlets, which were generally thin. The plot size was 18.71 m × 18.43 m. It comprised 5,042,403 points, with an average point density of approximately 24,413 points/m2, and there were 44 clearly identifiable larger trees. This plot had an LAI of approximately 4.57 and a crown cover of approximately 0.767. The average height of the trees in this plot was approximately 8 m. The average DBH of the trees in this plot was approximately 20 cm. Because it was a planted forest, its trees were evenly and uniformly distributed, with moderate spacing between them.
The sample from Plot 2 in Victoria, Australia, corresponded to a publicly available dataset that accurately represented a native eucalypt open forest ecosystem (Figure 1b). The predominant tree species in this plot was Eucalyptus leucoxylon. This species is an evergreen woody tree. A single-station scan of Plot 2 was performed using an RIEGL VZ-400 made in Beijing, China. It works in the near-infrared band at 1550 nm. Specifically, the forest area in this dataset was classified as a dry sclerophyll box–ironbark forest [35]. A portion of the dataset was selected for experimental purposes. In the plot selected for the experiment, the average tree height was approximately 20 m and the average DBH was approximately 25 cm. The size of the plot was 52.23 m × 54.44 m. It comprised 6,388,675 points, with an average point density of approximately 4934 points/m2, and there were 42 clearly identifiable larger trees. This plot had an LAI of approximately 1.438 and a crown cover of approximately 0.272. The tree distribution was relatively sparse, with greater space between the trees than in Plot 1. Each tree had fewer leaves, resulting in a lower overall crown cover.
The sample from Plot 3 was from a publicly available dataset from the Bretten municipal forest in the federal state of Baden-Württemberg, Germany (Figure 1c) [36]. The plot was scanned by multiple stations as well as an RIEGL VZ-400. For the experiment, a portion of the dataset from BR03 was selected. In this plot, the average tree height was approximately 24 m. The thickest trees had DBH values of approximately 72 cm, and the thinner ones had values of only 15 cm. The size of the plot was 24.73 m × 19.45 m. It comprised 8,833,993 points, with an average point density of approximately 18,366 points/m2, and there were 35 clearly identifiable larger trees. This plot had an LAI of approximately 9.49 and a crown cover of approximately 0.964. Since it was a natural forest, it starkly contrasted with the neat, orderly, and unified forest of Plot 1. In this plot, the tree distribution was relatively dense, with various tree species present, and the overall canopy cover in this region was relatively high.

3. Methods

The method proposed in this study uses cluster size, the number of points, and Chaos Distance as key criteria for point cloud classification and separation. This method has four main steps (Figure 2). First, the point cloud is filtered using the Cloth Simulation Filtering (CSF) method, which mainly removes ground points to minimize their interference in the subsequent process [37]. Second, the point cloud is clustered using a region-growing algorithm. It can put wood and leaves in different clusters for better identification. Third, preliminary separation is performed based on size and the number of points, which mainly extracts obvious wood (i.e., the main trunk and the thicker branches of a tree). Finally, the refined final wood separation is accomplished based on the Chaos Distance of the clusters. Chaos Distance is employed to depict the distinct characteristics of wood and leaves in TLS data. In other words, wood tends to exhibit a more uniform morphology, whereas leaves typically display a more erratic distribution. This parameter is used to complete the classification in this step.

3.1. Extraction of Nonground Points

This step involves preprocessing the TLS data, primarily removing the ground points in the plots to reduce their interference in the subsequent process. Most traditional filtering algorithms use differences in slope, elevation, and other parameters to distinguish ground points from nonground points. However, the CSF algorithm proposed by Zhang et al. [37] introduced a new idea for filtering ground points. Currently, the CSF algorithm is widely used to remove ground points from point clouds because it has fewer parameter settings and higher accuracy [37,38].
The algorithm first inverts a point cloud and then assumes that a piece of cloth falls from above due to gravity, so the final position of the cloth represents the current sample terrain. It considers that the fabric particles are affected by both internal and external forces. In this study, the ground points in the three sample plots were effectively removed using this algorithm.

3.2. Point Cloud Clustering

This step involves clustering a preprocessed forest point cloud and grouping morphologically similar points into the same cluster. Due to the significant morphological differences between wood and leaves, clustering can be used to group them into different clusters. This step aims to separate wood and leaf points into different clusters in preparation for the subsequent separation process.
The region-growing clustering algorithm can be used to separate wood and leaves into different clusters. The fundamental principle of regional growth is to aggregate points with similar characteristics to form cohesive regions [30,39]. The resulting clusters are organized into one set or multiple sets through curvature-based region growth clustering, each representing a smooth surface. First, the normal vectors are calculated by selecting neighborhood points, computing the covariance matrix, and performing eigenvalue decomposition to identify the normal vectors based on the local point cloud geometry. The directions of the normal vectors are then determined based on the angles between the lines [30]. A region needs to grow from the point with the smallest curvature value because this point is in a flat region (growing from the flattest area can reduce the total number of line segments). Therefore, a point cloud sorts its points by curvature value [40,41].
Subsequently, for a sorted point cloud ( P sorted ), we select the point with the smallest curvature value ( P min ) as the starting seed point and initialize the current region ( R c ) and current seed list ( S c ). The algorithm starts from the seed point and finds the set of nearest neighbors through a neighborhood lookup. For each nearest neighbor point ( P j ), we check that it satisfies the following two conditions:
  • Whether the P j is still in the untreated P sorted ;
  • Whether the angle θ ij satisfies θ ij < θ th , where θ ij is the angle between the P min and the P j in the normal vectors’ direction and θ th is the preset angle threshold.
If the conditions are met, the P j is added to the R c and removed from the P sorted . Otherwise, if the curvature value of the P j ( cP j ) is less than the curvature threshold ( cP th ), the P j is added to the S c . The algorithm continuously expands the R c , aggregating points with less curvature and similar normal vector directions into the same region. The iterative process of the algorithm is expressed through the following formulas and conditions:
θ ij = cos − 1 N → P min ⋅ N → P j
cP j < cP th
Eventually, the algorithm outputs a global list of segments (R) containing localized regions with similar curvatures and normal vectors [30]. The partial results of this algorithm are visualized in Figure 3, where each color represents a local cluster.

3.3. Preliminary Wood Separation

The purpose of this step is to preliminarily isolate clusters of prominent wood parts, such as the main trunk of a tree, as well as the more obvious thicker branches.
To optimize computational efficiency, preliminary separation prioritizes the identification of distinct wood clusters using relatively simple and computationally fast parameters, such as size and the number of points. During region-growing clustering, it is essentially impossible to cluster a large number of leaf points (i.e., many leaves) into a single cluster due to the utilization of normal vector features.
Therefore, if a cluster contains a large number of points, it is considered a wood cluster. Additionally, leaf clusters are unlikely to have significant height differences along the z-axis, whereas this feature is evident in the trunk. If a cluster has a large height difference, i.e., a large difference in the z-axis direction, it is considered a wood cluster. This can be expressed as follows:
W p = C | Z c > Z t ,   or   N c > N t
where W p denotes the preliminarily separated wood, C denotes the clusters, Z t denotes the difference between the class clusters in the Z-axis direction, N c denotes the number of points in the clusters, Z t denotes the threshold value set for the class clusters in the Z-axis direction, and N t denotes the threshold value set for the number of points. In general, setting Z t to 2.5 m and N c to 5000 means that misclassification of leaves can be avoided in this step with better separation of obvious wood.

3.4. Final Wood Separation

After preliminary wood separation, a more refined extraction and separation process is conducted. The purpose of this step is to classify all clusters that were not classified in the previous step as wood or leaf clusters. After preliminary separation, the remaining cluster types can be grouped into three categories, as shown in Figure 4: clusters containing mainly branches, clusters containing single leaves or a few leaves, and clusters containing mainly trunks. The characteristic that emerges here is a gradual increase in disorganization. We find that we can express this feature with the help of root-mean-square error ( RMSE ) calculation. In statistics, RMSE is a good indicator of the precision of a measurement and is calculated as follows:
RMSE = ∑ i = 1 n X obs , i −   X model , i 2 n
In this equation, X obs , i and X model , i refer to the observed and predicted values of the ith data, respectively, and n is the number of observation points.
This equation cannot be used directly in a point cloud. There are some methods that can fit the model in a point cloud, such as least squares, random sample consensus (RANSAC), etc., but these methods require large computational resources and can directly affect efficiency. Therefore, we chose to process cluster point clouds directly. If a point cloud is considered a “model”, the error between the observed and modeled points can be expressed as follows:
E i , j = ( X i − X j ) 2 + ( Y i − Y j ) 2 + ( Z i − Z j ) 2 3
The above equation represents the error of the ith point with respect to the jth point, where X i ,   Y i ,   Z i ,   X j ,   Y j , and   Z j are the X ,   Y ,   and   Z coordinates of the ith and jth points.
Iteration continues through all points in a cluster using the above equation, and results accumulate. Metrics are calculated from these accumulated results to assess the point cloud disorder within a cluster. After bringing the dimensions (units) of all variables into the equation, the final result is expressed in meters (m). Therefore, we named it the Chaos Distance ( d c ). After removing constants that may affect computational efficiency and simplifying the formula, the following equation can be obtained:
d c = 2 ∗ ∑ i = 1 n − 1 ∑ j = i + 1 n X i − X j 2 + Y i − Y j 2 + Z i − Z j 2 n ∗ n + 1
where the unit of d c is meters ( m ); n is the number of points in the class cluster; X i ,   Y i ,   Z i ,   X j ,   Y j , and   Z j are the X ,   Y , and   Z coordinates of the ith and jth points; and i   and   j are counted in n . In practice, the calculated d c is generally multiplied by a scaling factor to amplify the differences between the data, leading to better separation. This scaling factor is generally 10.
It is intuitive to think that the d c becomes progressively larger from (a) to (c) in the types of clusters represented in Figure 4. Specifically, among similarly sized clusters, wood clusters appear denser and more organized, while leaf clusters are looser and more disorganized. Therefore, for similar-sized clusters, the d c of wood is smaller than the d c of leaves. As mentioned in Section 3.2, it is rare for a large number of leaves to be clustered together. At the same time, the class composed of trunks is significantly larger than case (b) in Figure 4 when E i , j is calculated at the more distant end; i.e., the d c of case (c) in Figure 4 is larger than that of case (b). Therefore, the d c of leaf clusters is in a threshold range. Values beyond this threshold indicate wood, as shown by the distribution in Figure 4d. These concepts can be expressed using the following equation:
W f = C | d c d St ,   or   d c d Lt ,   or   C ∈ W p
where W f refers to the final set of separated wood and d St and d Lt refer to the larger and smaller thresholds, respectively. The threshold range is determined after a specific analysis based on different situations. Generally, the threshold range is [0.5, 1.5], which should be larger if the sample plot has high canopy closure or, conversely, smaller if the canopy closure is low. See Section 5.1 for details on setting and evaluating the threshold range.

3.5. Accuracy Assessment

The points in the processing results can be divided into four types: true separated wood ( W T ), true separated leaf ( L T ), false separated wood ( W F ), and false separated leaf ( L F ) [21]. To evaluate the method’s accuracy, manually separated data from the same plots were compared with the results in this study. A series of metrics were used to evaluate the accuracy of the plots’ separation results, including accuracy ( Ac ), wood precision ( Pr - Wood ), leaf precision ( Pr - Leaf ), the wood recall rate ( Re - Wood ), the leaf recall rate ( Re - Leaf ), the wood F1 score ( F 1 - Wood ), the leaf F1 score ( F 1 - Leaf ), the wood missed alarm rate ( Ma - Wood ), the leaf missed alarm rate ( Ma - leaf ), the wood false alarm rate ( Fa - Wood ), and the leaf false alarm rate ( Fa - Leaf ).
Ac refers to the proportion of measured values that meet limited conditions among multiple measured values under certain experimental conditions. It is calculated as follows [21]:
Ac = L T + W T L T + W T + L F + W F
Pr represents the proportion of samples with positive predictions and correct predictions among all samples with positive predictions [42]. High precision means that a sample recognized as positive is definitely positive. It is calculated as follows:
Pr = TP TP + FP
where TP denotes the number of positive categories predicted to be positive, FN denotes the number of positive categories predicted to be negative, FP denotes the number of negative categories predicted to be positive, and TN denotes the number of negative categories predicted to be negative [15,29]. In this experiment (including several subsequent indicators), the Pr of calculating wood is the result of wood being recorded as a positive class and leaves being recorded as a negative class. The opposite is true for the Pr of calculating leaves.
Re refers to the proportion of all positive samples that are correctly predicted to be positive. It is calculated for the original sample and means the probability of predicting a positive sample in a sample that is actually positive. A high recall rate means that there may be more false positives, but every effort will be made to find every object that should be found. Re is calculated as follows [31,43,44]:
Re = TP TP + FN
Pr and Re affect each other. Ideally, both are high, but the reality is that the two constrain each other: if a high Pr is pursued, the Re is lowered, and vice versa. The F 1 is a balance between Pr and Re , as it is the reconciled average of Pr and Re . The nature of reconciled averages is that they are only high if both values are very high. If one of them is low, the reconciled average will be pulled close to that very low number. F 1 is calculated as follows [31,43,44]:
F 1 = 2 ∗ Pr ∗ Re Pr + Re
The Fa refers to the probability of a negative example being misjudged as a positive example, and it is the ratio of false-positive examples. The Ma refers to the probability of misjudgment of positive examples and is the proportion of missed positive examples. They are a pair of tradeoff metrics, corresponding to Pr and Re , respectively. They are calculated as follows:
Fa = FP TP + FP
Ma = FN TP + FN

4. Results

To better characterize the locations and significant errors in the separation results of the three plots, we referred to other methods of vertical stratification for forest plots [45]. In this study, the canopy structures of the forests in the point cloud data were categorized into canopy, subcanopy, shrub, and ground layers, as shown in Figure 5. These layers did not have clear demarcations and only represented approximate areas. The canopy mainly contained leaves, while there were fewer branches. In the subcanopy, the proportions of branches and leaves were comparable. As for the shrub layer, it had more branches.
Figure 6a–c qualitatively represent the overall results and cropped sections of the final separation using the proposed method in the three sample plots. Red indicates separated wood, and the green indicates leaves in the figure. It can be seen that wood and leaves were successfully separated in the forest plots after processing with our method.
Figure 6a shows that leaves were mainly clustered in the canopy and subcanopy layers of the trees and were relatively neat in their growth in Plot 1. In this plot, the leaf distribution was relatively dense, and the method used in this study was more effective in separating the subcanopy layer. The main factor causing errors was the misclassification of leaves as wood in the canopy. As shown in Figure 6b, in Plot 2, leaves were primarily clustered in the canopy layer, while the bare branches in the subcanopy layer were distinct and could be extracted effectively. Due to the simple structures of the trees, it was easier to identify where errors occurred. Extracting wood from clusters with large numbers of overlapping leaves was challenging. In Plot 3, the structure was complex due to the fact that it was a natural forest with large numbers of leaves distributed in all layers. As can be seen in Figure 6c, the proposed method was accurate in identifying obvious wood, but many leaves were still misclassified as wood.
We used the experimental results and treated the manual segmentation as the reference for true separation. The results of the proposed method were compared with these manual results. The datasets were aligned, and the agreement at each point was checked. The numbers of points that were correctly or incorrectly categorized were then tallied accordingly. Figure 7 illustrates the contribution of each part, where a false leaf indicates a point where wood was incorrectly separated as a leaf. False wood is the opposite. It can be seen that the main cause of error was the misclassification of wood as leaves, in line with the response in Figure 6.
As can be seen in Figure 7, the point clouds of the three plots were dominated by leaves, ranging from 76% to 83%, and they contained relatively little wood. Therefore, when using the metrics for result assessment, clusters with larger leaf bases tended to show relatively high and less sensitive metrics for leaves, while the metrics for wood were more sensitive to errors.
For the separation results of the different tree species, Figure 8 shows the separation of the individual trees within the plots. Figure 8a,b are from Plot 1 and Plot 2, respectively, while Figure 8c,d represent different tree species within Plot 3. Among these four main tree types, the canopy structures of (a) and (c) were more complex compared to (b). There were more types of tree species in Plot 3, with (c) and (d) being two representative types. The trees belonging to type (c) had many branches, similar to (a). Due to the complex canopy, canopy shading was relatively severe. In contrast, type (d) had more small branches, and the leaves were mainly distributed at the ends of the branches, similar to (b). This type of tree had a simpler canopy structure. Despite the significant differences in tree types and structures in this experiment, good separation results were still achieved. Additionally, this method effectively addressed, to some extent, the errors commonly caused by canopy shading. In the illustrated trees, the main errors were caused by the numerous small branches present in (b) and (d).

5. Discussion

5.1. Sensitivity Analysis

In order to investigate the suitability of the selection of d c thresholds using the method proposed in this study, we selected different threshold ranges in each plot for comparison. The results of an accuracy assessment using different threshold ranges are shown in Table 2.
We chose three threshold ranges for comparison ([0, 1.5], [0.5, 1.5], and [0.5, 5]). As can be seen in the results for the three plots in Table 2, the Ac of [0.5, 1.5] was the highest. This shows that for the overall wood–leaf separation results, a d c range of [0.5, 1.5] was selected for the most favorable results. Using [0.5, 1.5] as a base, when comparing the results with the other ranges, both the Pr - wood and Re - leaf were found to be lower. In contrast, the Pr - leaf and Re - wood were the highest in all three plots.
We can explain the above phenomenon with the help of Figure 4. We can approximate it in this way: The dc of case (a) in Figure 4 was mainly distributed in the range [0, 0.5], and similarly, cases (b) and (c) were distributed in [0.5, 1.5] and [1.5, +∞], respectively, as shown in Figure 4d. Thus, it can be seen that in choosing the range [0, 1.5], we misclassified a number of (a) categories as leaves. Similarly, in choosing the range [0.5, 5], we misclassified a number of (c) categories as leaves. That is, the proportion of leaves that were correctly identified in our results decreased in this case. The proportion of all real wood that was correctly identified using our method also decreased. This is why Pr - leaf and Re - wood decreased in both ranges.
It is worth noting that this dividing line was not absolutely correct, and there was some crossover between the three categories. Crossover was the main cause of error in the [0.5, 1.5] range. The crossover led to another phenomenon: Changing the threshold correctly categorized all leaves in these crossed ranges in (b). That is, the proportion of wood that was correctly identified increased in this case. The proportion of all real leaves that were correctly identified also increased. This led to increases in Pr - wood and Re - leaf in the ranges [0, 1.5] and [0.5, 5].
As described above, different accuracy metrics increased or decreased when different threshold ranges were selected. Due to the crossover of the d c ranges in several cases, it was difficult to find a threshold range where all accuracy metrics were superior. As demonstrated by F 1 - wood and F 1 - Leaf in Table 2, after balancing the effects of both Pr and Re , the threshold range [0.5, 1.5] was clearly advantageous. Therefore, the d c range of [0.5, 1.5] was chosen for the evaluation of the results and comparison with other methods in this study.

5.2. Accuracy Validation and Method Comparison

We compared the method proposed in this study with the widely used LeWoS proposed by Wang et al. [46]. LeWoS used the default parameters in this experiment because it had fewer adjustable parameters and satisfactory robustness. The proposed method was also compared with a supervised learning method based on a Random Forest model (RF model) proposed by S. M. Krishna Moorthy et al. [26]. The model was trained by the authors using five plots located at two different locations and two single trees. They could be separated using the authors’ trained model directly. The model was compared to a number of other state-of-the-art methods described in the literature and could represent most machine learning models.
The results for the quantitative evaluation metrics for our method and the methods used for comparison (LeWoS and the RF model) are shown in Figure 9. For our method, the Ac of the three plots was higher than 0.9. Pr - wood was approximately 0.92 ± 0.03. Pr - leaf was approximately 0.94 ± 0.03. Re - wood was approximately 0.79 ± 0.05. Re - leaf was 0.97. F 1 - wood and F 1 - Leaf were 0.84 ± 0.02 and 0.96 ± 0.01, respectively. The lowest F 1 - wood was in Plot 2, which also reached 0.83, indicating the high accuracy of our method.
In comparison to LeWoS, our method had more metrics that were superior. Figure 9a,c show that our method was significantly superior to LeWoS in the processing of Plots 1 and 3. While Figure 9b shows that our method had some metrics that were superior to LeWoS in Plot 2, several metrics were lower. Overall, it was about on par with LeWoS. This was also true when compared to the RF model. Our method was significantly superior to the RF model for both Pr - wood and Re - leaf in the three plots. This indicates that our method is significantly more accurate in detecting wood than the RF model when confronted with different forest types. In the three plots, the RF model showed a greater degree of uncertainty, with drastic changes in accuracy.
In Plots 1 and 3, the accuracy of our method was significantly higher than those of LeWoS and the RF model due to the complex canopy structures of these two plots. This phenomenon indicates that our method addressed the problem of difficult segregation in complex forests to a certain extent. Meanwhile, for Plot 2, which had a relatively simple canopy structure, our method was better than LeWoS for several metrics, and the overall accuracy was basically the same. Although the RF model’s Re - leaf was significantly better than those of the other two methods in Plot 2, its Pr-wood was undeniably lower, which led to its F 1 - wood and F 1 - leaf values being lower than those of the other methods. Overall, the Re - wood of our method was relatively low in all three plots. This was consistent with the contributions. It also further suggested that F 1 - wood was less accurate in the three plots compared to the other metrics.
Among the three plots, the Pr - wood and Fa - wood of our method were significantly better than those of LeWoS and the RF model. This indicates that the correct proportion of the predicted wood was relatively high with our method. That is, the problem of classifying a large number of leaves as wood was addressed to a certain extent. The Re - wood of our method was higher in Plots 1 and 3 but lower in Plot 2. This shows that for a simple canopy structure, the LeWoS treatment was superior, but our method found wood easier in a complex structure. On the other hand, while the RF model performed better in simpler canopy structures, it still had lower performance compared to LeWoS.
The better accuracy of our method was mainly due to the point cloud clustering described in Section 3.2. Our method divided the points into clusters with only leaves or wood, which made them easier to detect (Figure 3). The combination of the above caused the F 1 - wood values of Plots 1 and 3 to be significantly higher than those of LeWoS and the RF model, while these values were the same in Plot 2. For leaves, the accuracy values of all methods were high for all indexes, which was due to the large base of leaves; that is, the proportion of leaves in the point cloud was relatively high (Figure 7).

5.3. Uncertainty Analysis

The cluster-based method proposed in this study showed effectiveness in three different plots. Through the detailed qualitative and quantitative results, it was shown that this method can achieve efficient and highly accurate wood–leaf separation in actual scenarios. In these plots, the proposed method effectively handled differences in geographical locations, tree species, and point cloud densities, indicating its versatility and robustness in diverse environments.
Specifically, detailed zoomed-in views are shown in Figure 6 and Figure 10. In Figure 10, the (2) in each figure represents a region with a relatively good separation effect. As can be seen in the figure, the proposed method also extracted relatively coarse wood well and was less affected by sheltering. Quantitatively, Table 3 shows a comparison between the plot parameters and the accuracy of the separation results. Table 3 illustrates that the proposed method showed strong robustness in the face of different parameters, such as the point density, the LAI, and so on.
However, some points were still misclassified. The occurrence of such situations can be attributed to two main factors. On the one hand, errors might have arisen during the initial clustering step, where some leaf points were mistakenly clustered with wood points. Additionally, for very thin wood, which sometimes comprised only a single row of points, accurately clustering it with other wood was challenging, leading to separation difficulties during the subsequent extraction. Figure 10 shows the local areas of the point cloud data that exhibited the above characteristics. On the other hand, setting the Chaos Distance threshold could have introduced errors because it had to be adjusted based on the specific characteristics of the plots. As mentioned in Section 5.1 and shown in the lower part of Figure 4, the dividing line between leaf and wood was not clear when choosing the Chaos Distance threshold range. Inaccurate threshold settings might have led to deviations or errors in the separation process.
Other factors might have contributed to the observed situations. One such factor is the inherent complexity and variability of the natural environment. The data might have contained inherent noise, occlusions, or irregularities caused by vegetation density, lighting conditions, and sensor limitations. These factors might have introduced uncertainties and made it challenging to separate wood and leaves accurately.
Therefore, the combination of these factors, including the inherent challenges posed by clustering and the sensitivity of the parameter settings, contributed to inaccuracies in the wood–leaf separation process.

5.4. Efficiency Comparison and Analysis

In terms of the efficiency of the algorithm, the method had a time complexity of O n   log   n , where n is the total number of points in the point cloud. The actual execution duration varied depending on the size of the point cloud data and the specific parameters used for clustering and classification. C++ was the programming language used to implement the wood–leaf separation algorithm proposed in this study. The primary third-party library was the Point Cloud Library (PCL) [31]. The specific implementation algorithm was compiled using the Microsoft Visual Studio Community 2022 development environment. The open-source algorithm LeWoS, which was used for comparison, was run on the same computer using MATLAB Runtime [3]. The RF model algorithm used for the comparison was run using the Python 3.7 environment with version 0.21.3 of the scikit-learn library. The operating system was Windows 11 Professional Edition 22H2. The processor was an eleventh-generation Intel Core i5-11400H @ 2.70 GHz with 16.0 GB of RAM. Due to the large amount of data, which could easily exceed the computer’s memory, we opened an additional 30 GB of virtual memory (this was to “pretend” that a part of the hard disk space was for memory use, and the disadvantage was that this part of the space was slower to read and write). The recorded plot processing times of our method, LeWoS, and the RF model are shown in Table 4, which indicates the processing durations for the three study plots. The average processing duration was 7.82 min. It should be noted that the processing times for all methods are the times required to process the preprocessed data. All three methods used the same preprocessing method, as described in Section 3.1.
Compared with LeWoS and the RF model, the proposed method demonstrated significantly better performance when processing the same data. The main time consumption of the proposed method occurred during the estimation of the point cloud normal vectors and the regional growth clustering process, which was evidently related to the amount of data. The running time of the method in this study was directly related to the number of points in the point cloud and the number of clusters after clustering. The number of clusters was related to the characteristics of the data and the parameter settings used for clustering (such as the smoothing threshold, curvature threshold, the number of search points, etc.). Specifically, the focus was on the need to obtain point cloud neighbor information during normal vector estimation for point cloud clustering. This made it necessary to perform a traversal process on the entire point cloud. When clustering, different parameters, such as smoothing thresholds, curvature thresholds, the number of search points, etc., can lead to large differences in the cluster size, the degree of smoothing, and other characteristics. These differences may lead to significant time costs, such as variations in the frequency of updating seed points due to different parameter settings. In the subsequent separation process, in order to increase the efficiency of the algorithm, we first performed a preliminary separation of the wood in the clustering results. This effectively extracted the wood clusters that contained more points. In the final separation, we were only faced with the remaining clusters, which contained fewer points. The traversal of these clusters could be carried out efficiently to achieve the final separation. In addition, another reason why our approach was more effective than LeWoS and the RF model could be the choice of programming language. It has been shown that C/C++ is significantly more efficient than other languages in terms of runtime [47].
For the three plots in the experiment, the main reason why the time for Plot 2 was shorter than that of Plot 1 was that the proportion of ground points in Plot 2 was relatively low. After removing the ground points, there were 2,820,945 points left in Plot 2 and 3,270,721 points left in Plot 1. In addition, as shown in Figure 6, in Plot 2, the trees had a simpler structure compared to Plot 1, where the trees had more branchlets and leaves. This made the algorithm take more time to perform clustering in Plot 1, and it produced more clusters. The processing time for Plot 3 was significantly longer than for the other two plots. This was partly because the proportion of ground points in the data was smaller and because the number of points was significantly higher. Additionally, the structure of the canopy in Plot 3 was complex, with more intersections between trees and a high proportion of leaves. This resulted in more clusters, which further increased the run time.

6. Conclusions

In this study, we proposed a method based on cluster features for separating leaves and wood in forest plots using terrestrial LiDAR point cloud data. The main idea behind this method is to leverage the morphological differences between wood and leaves for separation. A point cloud with nonground points is clustered into groups using a region-growing clustering algorithm. Preliminary wood separation is then performed based on the cluster size and the number of points. Subsequently, the final wood extraction is completed by combining the Chaos Distance, and finally, the separation accuracy of the wood–leaf separation is verified. The proposed method is simple to implement. The Chaos Distances proposed in this study demonstrated high accuracy and efficiency, with an average recall rate of 0.89 and an average precision of 0.92. This validated the effectiveness of separating wood from leaves in forest plots. The study plots had an average of 44 trees, and all plots were treated with high efficiency. In addition, this method had good effects on different tree species and forest areas. The method comparison results show that the proposed method has higher accuracy in complex canopy structures and dense forests. Future work could focus on integrating this method with other remote sensing technologies, which could enhance its application in large-scale forest monitoring and management.

Author Contributions

H.T. and S.L. designed and performed the experiments. H.T., S.L., Z.S. and Z.H. contributed to the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China under Grants 42371382 and 41871247.

Data Availability Statement

The Plot 1 data used to support the results of this study can be obtained by contacting the authors because these data were collected by the authors and are relevant to their privacy. The Plot 2 data were obtained from http://dx.doi.org/10.4227/05/542B766D5D00D (accessed in 2024). The Plot 3 data were obtained from https://doi.org/10.1594/PANGAEA.933426 (accessed in 2024).

Acknowledgments

The author would like to thank the Scientific research and education program of the School of Resources and Environment (SRE) of the University of Electronic Science and Technology of China (UESTC) for providing experimental equipment and support. The authors would also like to thank the anonymous reviewers and Assigned Editor for their detailed and constructive suggestions related to this study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Soliman, S.A.M.; Hussein, K.F.A.; Ammar, A.-E.-H.A. Electromagnetic Resonances of Natural Grasslands and Their Effects on Radar Vegetation Index. Prog. Electromagn. Res. B 2020, 86, 19–38. [Google Scholar] [CrossRef]
  2. Hui, G.; Zhang, G.; Zhao, Z.; Yang, A. Methods of Forest Structure Research: A Review. Curr For. Rep. 2019, 5, 142–154. [Google Scholar] [CrossRef]
  3. Wang, D.; Momo Takoudjou, S.; Casella, E. LeWoS: A Universal Leaf-wood Classification Method to Facilitate the 3D Modelling of Large Tropical Trees Using Terrestrial LiDAR. Methods Ecol. Evol. 2020, 11, 376–389. [Google Scholar] [CrossRef]
  4. Kalyoncu, C.; Toygar, Ö. Geometric Leaf Classification. Methods Ecol. Evol. 2015, 133, 102–109. [Google Scholar] [CrossRef]
  5. Wang, X.-F.; Huang, D.-S.; Du, J.-X.; Xu, H.; Heutte, L. Classification of Plant Leaf Images with Complicated Background. Appl. Math. Comput. 2008, 205, 916–926. [Google Scholar] [CrossRef]
  6. Turkoglu, M.; Hanbay, D. Leaf-Based Plant Species Recognition Based on Improved Local Binary Pattern and Extreme Learning Machine. Phys. A Stat. Mech. Its Appl. 2019, 527, 121297. [Google Scholar] [CrossRef]
  7. Lechner, A.M.; Foody, G.M.; Boyd, D.S. Applications in Remote Sensing to Forest Ecology and Management. One Earth 2020, 2, 405–412. [Google Scholar] [CrossRef]
  8. Beland, M.; Parker, G.; Sparrow, B.; Harding, D.; Chasmer, L.; Phinn, S.; Antonarakis, A.; Strahler, A. On Promoting the Use of Lidar Systems in Forest Ecosystem Research. For. Ecol. Manag. 2019, 450, 117484. [Google Scholar] [CrossRef]
  9. Lin, Y. LiDAR: An Important Tool for next-Generation Phenotyping Technology of High Potential for Plant Phenomics? Comput. Electron. Agric. 2015, 119, 61–73. [Google Scholar] [CrossRef]
  10. Kim, I.; Martins, R.J.; Jang, J.; Badloe, T.; Khadir, S.; Jung, H.-Y.; Kim, H.; Kim, J.; Genevet, P.; Rho, J. Nanophotonics for Light Detection and Ranging Technology. Nat. Nanotechnol. 2021, 16, 508–524. [Google Scholar] [CrossRef]
  11. Guo, Q.; Su, Y.; Hu, T.; Guan, H.; Jin, S.; Zhang, J.; Zhao, X.; Xu, K.; Wei, D.; Kelly, M.; et al. Lidar Boosts 3D Ecological Observations and Modelings: A Review and Perspective. IEEE Geosci. Remote Sens. Mag. 2021, 9, 232–257. [Google Scholar] [CrossRef]
  12. Reutebuch, S.E.; Andersen, H.-E.; McGaughey, R.J. Light Detection and Ranging (LIDAR): An Emerging Tool for Multiple Resource Inventory. J. For. 2005, 103, 286–292. [Google Scholar] [CrossRef]
  13. Akay, A.E.; Oğuz, H.; Karas, I.R.; Aruga, K. Using LiDAR Technology in Forestry Activities. Environ. Monit. Assess. 2009, 151, 117–125. [Google Scholar] [CrossRef] [PubMed]
  14. Liang, X.; Hyyppä, J.; Kaartinen, H.; Lehtomäki, M.; Pyörälä, J.; Pfeifer, N.; Holopainen, M.; Brolly, G.; Francesco, P.; Hackenberg, J.; et al. International Benchmarking of Terrestrial Laser Scanning Approaches for Forest Inventories. ISPRS-J. Photogramm. Remote Sens. 2018, 144, 137–179. [Google Scholar] [CrossRef]
  15. Su, Z.; Li, S.; Liu, H.; Liu, Y. Extracting Wood Point Cloud of Individual Trees Based on Geometric Features. IEEE Geosci. Remote Sens. Lett. 2019, 16, 1294–1298. [Google Scholar] [CrossRef]
  16. Hu, C.; Pan, Z.; Zhong, T. Leaf and Wood Separation of Poplar Seedlings Combining Locally Convex Connected Patches and K-Means++ Clustering from Terrestrial Laser Scanning Data. J. Appl. Rem. Sens. 2020, 14, 018502. [Google Scholar] [CrossRef]
  17. Sun, J.; Wang, P.; Gao, Z.; Liu, Z.; Li, Y.; Gan, X.; Liu, Z. Wood–Leaf Classification of Tree Point Cloud Based on Intensity and Geometric Information. Remote Sens. 2021, 13, 4050. [Google Scholar] [CrossRef]
  18. Tan, K.; Zhang, W.; Dong, Z.; Cheng, X.; Cheng, X. Leaf and Wood Separation for Individual Trees Using the Intensity and Density Data of Terrestrial Laser Scanners. IEEE Trans. Geosci. Remote Sens. 2021, 59, 7038–7050. [Google Scholar] [CrossRef]
  19. Côté, J.-F.; Fournier, R.A.; Egli, R. An Architectural Model of Trees to Estimate Forest Structural Attributes Using Terrestrial LiDAR. Environ. Modell. Softw. 2011, 26, 761–777. [Google Scholar] [CrossRef]
  20. Hui, Z.; Jin, S.; Xia, Y.; Wang, L.; Yevenyo Ziggah, Y.; Cheng, P. Wood and Leaf Separation from Terrestrial LiDAR Point Clouds Based on Mode Points Evolution. ISPRS-J. Photogramm. Remote Sens. 2021, 178, 219–239. [Google Scholar] [CrossRef]
  21. Tian, Z.; Li, S. Graph-Based Leaf–Wood Separation Method for Individual Trees Using Terrestrial Lidar Point Clouds. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5705111. [Google Scholar] [CrossRef]
  22. Vicari, M.B.; Disney, M.; Wilkes, P.; Burt, A.; Calders, K.; Woodgate, W. Leaf and Wood Classification Framework for Terrestrial LiDAR Point Clouds. Methods Ecol. Evol. 2019, 10, 680–694. [Google Scholar] [CrossRef]
  23. Zhou, J.; Wei, H.; Zhou, G.; Song, L. Separating Leaf and Wood Points in Terrestrial Laser Scanning Data Using Multiple Optimal Scales. Sensors 2019, 19, 1852. [Google Scholar] [CrossRef] [PubMed]
  24. Zhu, X.; Skidmore, A.K.; Darvishzadeh, R.; Niemann, K.O.; Liu, J.; Shi, Y.; Wang, T. Foliar and Woody Materials Discriminated Using Terrestrial LiDAR in a Mixed Natural Forest. Int. J. Appl. Earth Obs. Geoinf. 2018, 64, 43–50. [Google Scholar] [CrossRef]
  25. Li, S.; Dai, L.; Wang, H.; Wang, Y.; He, Z.; Lin, S. Estimating Leaf Area Density of Individual Trees Using the Point Cloud Segmentation of Terrestrial LiDAR Data and a Voxel-Based Model. Remote Sens. 2017, 9, 1202. [Google Scholar] [CrossRef]
  26. Krishna Moorthy, S.M.; Calders, K.; Vicari, M.B.; Verbeeck, H. Improved Supervised Learning-Based Approach for Leaf and Wood Classification from LiDAR Point Clouds of Forests. IEEE Trans. Geosci. Remote Sens. 2020, 58, 3057–3070. [Google Scholar] [CrossRef]
  27. Wang, D.; Brunner, J.; Ma, Z.; Lu, H.; Hollaus, M.; Pang, Y.; Pfeifer, N. Separating Tree Photosynthetic and Non-Photosynthetic Components from Point Cloud Data Using Dynamic Segment Merging. Forests 2018, 9, 252. [Google Scholar] [CrossRef]
  28. Sun, S.; Cao, Z.; Zhu, H.; Zhao, J. A Survey of Optimization Methods from a Machine Learning Perspective. IEEE Trans. Cybern. 2020, 50, 3668–3681. [Google Scholar] [CrossRef]
  29. Ma, L.; Zheng, G.; Eitel, J.U.H.; Moskal, L.M.; He, W.; Huang, H. Improved Salient Feature-Based Approach for Automatically Separating Photosynthetic and Nonphotosynthetic Components within Terrestrial Lidar Point Cloud Data of Forest Canopies. IEEE Trans. Geosci. Remote Sens. 2016, 54, 679–696. [Google Scholar] [CrossRef]
  30. PCL Point Cloud Library (PCL). Available online: https://github.com/PointCloudLibrary/pcl/blob/master/doc (accessed on 9 January 2024).
  31. Rusu, R.B.; Cousins, S. 3D Is Here: Point Cloud Library (PCL). In Proceedings of the 2011 IEEE International Conference on Robotics and Automation, Shanghai, China, 9–13 May 2011; IEEE: Piscataway, NJ, USA, 2011; pp. 1–4. [Google Scholar]
  32. Tan, K.; Ke, T.; Tao, P.; Liu, K.; Duan, Y.; Zhang, W.; Wu, S. Discriminating Forest Leaf and Wood Components in TLS Point Clouds at Single-Scan Level Using Derived Geometric Quantities. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5701517. [Google Scholar] [CrossRef]
  33. Wan, P.; Shao, J.; Jin, S.; Wang, T.; Yang, S.; Yan, G.; Zhang, W. A Novel and Efficient Method for Wood–Leaf Separation from Terrestrial Laser Scanning Point Clouds at the Forest Plot Level. Methods Ecol. Evol. 2021, 12, 2473–2486. [Google Scholar] [CrossRef]
  34. LiDAR360. Available online: https://www.lidar360.com/ (accessed on 9 January 2024).
  35. Calders, K.; Newnham, G.; Burt, A.; Murphy, S.; Raumonen, P.; Herold, M.; Culvenor, D.; Avitabile, V.; Disney, M.; Armston, J.; et al. Nondestructive Estimates of Above-ground Biomass Using Terrestrial Laser Scanning. Methods Ecol. Evol. 2015, 6, 198–208. [Google Scholar] [CrossRef]
  36. Weiser, H.; Schäfer, J.; Winiwarter, L.; Krašovec, N.; Seitz, C.; Schimka, M.; Anders, K.; Baete, D.; Braz, A.S.; Brand, J.; et al. Terrestrial, UAV-Borne, and Airborne Laser Scanning Point Clouds of Central European Forest Plots, Germany, with Extracted Individual Trees and Manual Forest Inventory Measurements; PANGAEA: Bremen, Germany, 2021. [Google Scholar] [CrossRef]
  37. Zhang, W.; Qi, J.; Wan, P.; Wang, H.; Xie, D.; Wang, X.; Yan, G. An Easy-to-Use Airborne LiDAR Data Filtering Method Based on Cloth Simulation. Remote Sens. 2016, 8, 501. [Google Scholar] [CrossRef]
  38. Wang, Z.; Yu, B.; Chen, J.; Liu, C.; Zhan, K.; Sui, X.; Xue, Y.; Li, J. Research on Lidar Point Cloud Segmentation and Collision Detection Algorithm. In Proceedings of the 2019 6th International Conference on Information Science and Control Engineering (ICISCE), Shanghai, China, 20–22 December 2019; IEEE: Piscataway, NJ, USA, 2019; pp. 475–479. [Google Scholar]
  39. Vo, A.-V.; Truong-Hong, L.; Laefer, D.F.; Bertolotto, M. Octree-Based Region Growing for Point Cloud Segmentation. ISPRS-J. Photogramm. Remote Sens. 2015, 104, 88–100. [Google Scholar] [CrossRef]
  40. Li, Y.; Liu, J.; Zhang, B.; Wang, Y.; Yao, J.; Zhang, X.; Fan, B.; Li, X.; Hai, Y.; Fan, X. Three-Dimensional Reconstruction and Phenotype Measurement of Maize Seedlings Based on Multi-View Image Sequences. Front. Plant Sci. 2022, 13, 974339. [Google Scholar] [CrossRef]
  41. Yu, D.; He, L.; Ye, F.; Jiang, L.; Zhang, C.; Fang, Z.; Liang, Z. Unsupervised Ground Filtering of Airborne-Based 3D Meshes Using a Robust Cloth Simulation. Int. J. Appl. Earth Obs. Geoinf. 2022, 111, 102830. [Google Scholar] [CrossRef]
  42. Yu, D.; Li, A.; Li, J.; Xu, Y.; Long, Y. Mean Inflection Point Distance: Artificial Intelligence Mapping Accuracy Evaluation Index—An Experimental Case Study of Building Extraction. Remote Sens. 2023, 15, 1848. [Google Scholar] [CrossRef]
  43. Yang, P.; Fu, H.; Zhu, J.; Li, Y.; Wang, C. An Elliptical Distance Based Photon Point Cloud Filtering Method in Forest Area. IEEE Geosci. Remote Sens. Lett. 2022, 19, 6504705. [Google Scholar] [CrossRef]
  44. Ji, F.; Ming, D.; Zeng, B.; Yu, J.; Qing, Y.; Du, T.; Zhang, X. Aircraft Detection in High Spatial Resolution Remote Sensing Images Combining Multi-Angle Features Driven and Majority Voting CNN. Remote Sens. 2021, 13, 2207. [Google Scholar] [CrossRef]
  45. Nine Layers of the Edible Forest Garden. Available online: https://tcpermaculture.com/site/plant-index/ (accessed on 2 September 2024).
  46. LeWoS. Available online: https://github.com/dwang520/LeWoS (accessed on 5 December 2020).
  47. Pereira, R.; Couto, M.; Ribeiro, F.; Rua, R.; Cunha, J.; Fernandes, J.P.; Saraiva, J. Energy Efficiency across Programming Languages: How Do Energy, Time, and Memory Relate? In Proceedings of the 10th ACM SIGPLAN International Conference on Software Language Engineering, Vancouver, BC, Canada, 23 October 2017; ACM: New York, NY, USA, 2017; pp. 256–267. [Google Scholar]
Remotesensing 16 03355 g001
Figure 1. Schematics of the three locations in this study. The exact locations of the plots are marked in yellow. a, b, and c show the three areas in China, Australia, and Germany, respectively. (a_1–c_1) are remote sensing satellite maps. (a_2–c_2) are original data from representative point clouds at the locations.
Figure 1. Schematics of the three locations in this study. The exact locations of the plots are marked in yellow. a, b, and c show the three areas in China, Australia, and Germany, respectively. (a_1–c_1) are remote sensing satellite maps. (a_2–c_2) are original data from representative point clouds at the locations.
Remotesensing 16 03355 g001
Remotesensing 16 03355 g002
Figure 2. The workflow of this study. In this figure, the yellow boxes indicate data or the results of processing, and the blue boxes indicate processes or algorithms used during processing.
Figure 2. The workflow of this study. In this figure, the yellow boxes indicate data or the results of processing, and the blue boxes indicate processes or algorithms used during processing.
Remotesensing 16 03355 g002
Remotesensing 16 03355 g003
Figure 3. Partial results after clustering by the region-growing clustering algorithm. Each color in the figure represents a cluster.
Figure 3. Partial results after clustering by the region-growing clustering algorithm. Each color in the figure represents a cluster.
Remotesensing 16 03355 g003
Remotesensing 16 03355 g004
Figure 4. Schematic representation of the cluster types remaining after clustering and initial wood separation and their Chaos Distance distributions. (a) contains four clusters, each of which contains wood. This type of cluster mainly consists of branches. (b) contains four clusters and denotes the type of cluster that contains one leaf or a few leaves. (c) is composed of wood, specifically showing the type of cluster that contains a trunk. The rough situation of the distribution of Chaos Distances for each type of cluster is indicated in (d), corresponding to the colors in the upper half. The horizontal axis represents the size of the d c in meters, and the vertical axis represents the approximate probability density.
Figure 4. Schematic representation of the cluster types remaining after clustering and initial wood separation and their Chaos Distance distributions. (a) contains four clusters, each of which contains wood. This type of cluster mainly consists of branches. (b) contains four clusters and denotes the type of cluster that contains one leaf or a few leaves. (c) is composed of wood, specifically showing the type of cluster that contains a trunk. The rough situation of the distribution of Chaos Distances for each type of cluster is indicated in (d), corresponding to the colors in the upper half. The horizontal axis represents the size of the d c in meters, and the vertical axis represents the approximate probability density.
Remotesensing 16 03355 g004
Remotesensing 16 03355 g005
Figure 5. Schematic representation of the distribution of the canopy structure in the forest. The point cloud colors in the figure are shown using height information.
Figure 5. Schematic representation of the distribution of the canopy structure in the forest. The point cloud colors in the figure are shown using height information.
Remotesensing 16 03355 g005
Remotesensing 16 03355 g006
Figure 6. The results of separating wood and leaves using this method, with the overall separation on the left, the separation of the clipped sample plot on the right in the upper half, and the corresponding wood points in the lower half, where the sample plot was clipped for better visibility. The bottom half of (a) and the right halves of (b,c) contain the results of single trees and partial canopies in each plot after separation. These detailed figures show the separation situations and the wood results in the corresponding areas. (a) Chinese sample plot (Plot 1). (b) Australian sample plot (Plot 2). (c) German sample plot (Plot 3).
Figure 6. The results of separating wood and leaves using this method, with the overall separation on the left, the separation of the clipped sample plot on the right in the upper half, and the corresponding wood points in the lower half, where the sample plot was clipped for better visibility. The bottom half of (a) and the right halves of (b,c) contain the results of single trees and partial canopies in each plot after separation. These detailed figures show the separation situations and the wood results in the corresponding areas. (a) Chinese sample plot (Plot 1). (b) Australian sample plot (Plot 2). (c) German sample plot (Plot 3).
Remotesensing 16 03355 g006
Remotesensing 16 03355 g007
Figure 7. The contribution of each part of the separation results.
Figure 7. The contribution of each part of the separation results.
Remotesensing 16 03355 g007
Remotesensing 16 03355 g008
Figure 8. Separation of individual trees of different species in the plots. (a) is a tree in Plot 1, and (b) is in Plot 2. (c,d) are different tree species in Plot 3.
Figure 8. Separation of individual trees of different species in the plots. (a) is a tree in Plot 1, and (b) is in Plot 2. (c,d) are different tree species in Plot 3.
Remotesensing 16 03355 g008
Remotesensing 16 03355 g009
Figure 9. The quantitative accuracy of our method, LeWoS, and the RF model was evaluated. (a–c) correspond to Plot 1, Plot 2, and Plot 3, respectively.
Figure 9. The quantitative accuracy of our method, LeWoS, and the RF model was evaluated. (a–c) correspond to Plot 1, Plot 2, and Plot 3, respectively.
Remotesensing 16 03355 g009
Remotesensing 16 03355 g010
Figure 10. Localized magnified views of the separation results. (a–c) correspond to Plots 1, 2, and 3, respectively. In each figure, the figures of overall views in (1) show partial separation results, (2) the red box figures show the locations with better separation, and (3) the blue box figures show the locations with errors.
Figure 10. Localized magnified views of the separation results. (a–c) correspond to Plots 1, 2, and 3, respectively. In each figure, the figures of overall views in (1) show partial separation results, (2) the red box figures show the locations with better separation, and (3) the blue box figures show the locations with errors.
Remotesensing 16 03355 g010
Table 1. Summary of sample plot information.
Table 1. Summary of sample plot information.
LocationSize of Plot (m)Number of PointsPoint
Density (Points/m2)		Number of TreesLAICanopy Cover
Plot 1China18.71 × 18.435,042,40324,413444.570.767
Plot 2Australia52.23 × 54.446,388,6754934421.4380.272
Plot 3Germany24.73 × 19.458,833,99318,366359.490.964
Table 2. Accuracy evaluation results for different threshold ranges.
Table 2. Accuracy evaluation results for different threshold ranges.
Threshold RangeAccuracyPr-WoodPr-LeafRe-WoodRe-LeafF1-WoodF1-Leaf
Plot 1[0, 1.5]0.8560.9950.8410.4140.9990.5950.913
[0.5, 1.5]0.9370.8940.9500.8430.9680.8680.959
[0.5, 5]0.8630.9310.8530.4730.9890.6280.916
Plot 2[0, 1.5]0.9140.9690.9010.6910.9920.8060.944
[0.5, 1.5]0.9220.9520.9140.7380.9870.8320.949
[0.5, 5]0.9170.9710.9040.7020.9930.8150.946
Plot 3[0, 1.5]0.9340.8910.9410.7110.9820.7910.961
[0.5, 1.5]0.9530.8840.9670.8410.9770.8640.972
[0.5, 5]0.9330.9720.9280.6330.9960.7670.961
Table 3. Comparison between plot parameters and accuracy of separation results.
Table 3. Comparison between plot parameters and accuracy of separation results.
Point
Density (Points/m2)		LAICanopy CoverAccuracyPr-WoodPr-LeafRe-WoodRe-LeafF1-WoodF1-Leaf
Plot 124,4134.570.7670.9370.8940.9500.8420.9670.8670.958
Plot 249341.4380.2720.9220.9520.9140.7370.9860.8310.949
Plot 318,3669.490.9640.9530.8830.9670.8400.9760.8630.971
Table 4. Comparison of processing times for sample plots.
Table 4. Comparison of processing times for sample plots.
Number of PointsTime Cost (min)
Our MethodLeWoSRF Model
Plot 15,042,4035.188.87313.37
Plot 26,388,6753.7510.73223.37
Plot 38,833,99314.5220.41523.3
Average6,755,0247.8213.33686.68
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tang, H.; Li, S.; Su, Z.; He, Z. Cluster-Based Wood–Leaf Separation Method for Forest Plots Using Terrestrial Laser Scanning Data. Remote Sens. 2024, 16, 3355. https://doi.org/10.3390/rs16183355

AMA Style

Tang H, Li S, Su Z, He Z. Cluster-Based Wood–Leaf Separation Method for Forest Plots Using Terrestrial Laser Scanning Data. Remote Sensing. 2024; 16(18):3355. https://doi.org/10.3390/rs16183355

Chicago/Turabian Style

Tang, Hao, Shihua Li, Zhonghua Su, and Ze He. 2024. "Cluster-Based Wood–Leaf Separation Method for Forest Plots Using Terrestrial Laser Scanning Data" Remote Sensing 16, no. 18: 3355. https://doi.org/10.3390/rs16183355

APA Style

Tang, H., Li, S., Su, Z., & He, Z. (2024). Cluster-Based Wood–Leaf Separation Method for Forest Plots Using Terrestrial Laser Scanning Data. Remote Sensing, 16(18), 3355. https://doi.org/10.3390/rs16183355

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop