Assessing the Performance of Handheld Laser Scanning for Individual Tree Mapping in an Urban Area

: Precise individual tree or sample-based inventories derived from 3D point cloud data of mobile laser scanning can improve our comprehensive understanding of the structure, function, resilience, biodiversity, and ecosystem services of urban forests. This study assessed the performance of a handheld laser scanning system (HLS) for the extraction of tree position, diameter at breast height (DBH), and tree height (H) in an urban area. A total of 2083 trees of 13 species from 34 plots were analyzed. The results showed that the registration of tree positions using ground control points (GCPs) demonstrated high accuracy, with errors consistently below 0.4 m, except for a few instances. The extraction accuracy of DBH for all trees and individual species remained consistently high, with a total root mean square error (RMSE) of 2.06 cm (6.89%) and a bias of 0.62 cm (2.07%). Notably, broad-leaved trees outperformed coniferous trees, with RMSE and bias values of 1.86 cm (6%) and 0.76 cm (2.46%), respectively, compared to 2.54 cm (9.46%) and 0.23 cm (0.84%), respectively. The accuracy of H extraction varied significantly among different species, with R 2 values ranging from 0.65 to 0.92. Generally, both DBH and H were underestimated compared to ground measurements. Linear mixed-effects models (LMEs) were applied to evaluate factors affecting the performance of HLS with the plot as a random factor. LME analysis revealed that plant type and terrain significantly influenced the accuracy of DBH and H derived from HLS data, while other fixed factors such as plot area, tree density, and trajectory length showed no significance. With a large sample size, we concluded that the HLS demonstrated sufficient accuracy in extracting individual tree parameters in urban forests.


Introduction
Urban forests, with their minimal spatial footprint, provide significant ecosystem services for urban residents and wildlife.These services include promoting health and social well-being, enhancing children's cognitive development and educational success rate, fostering a strong economy and providing numerous resources, mitigating the urban heat island effect, storing and sequestering carbon, offering habitat and food for animals, and managing stormwater through green infrastructure [1,2].Precise inventories of individual trees or sample-based parameters are fundamental for achieving a comprehensive understanding of the structure, function, resilience, biodiversity, and ecosystem services of urban forests [3,4].In traditional plot-based forest inventories, the most frequently measured individual attributes are tree species, diameter at breast height (DBH), and tree height [5,6].These tree attributes, either individually or in combination using species-specific equations (allometric models), are employed to compute and estimate a range of characteristics at the Forests 2024, 15, 575 2 of 21 tree, plot, and stand levels, including basal area, volume, biomass, carbon stock, and others [7][8][9][10].However, the diversity of urban forests, the nuanced impacts of urban trees at a small scale, and the pronounced spatial variations in urban areas create unique site-specific differences that must be taken into account to obtain highly accurate measurements [11].To address these site-specific variations, comprehensive and high-resolution data at local scales are essential [12].
In recent decades, close-range remote sensing has undergone rapid development, fundamentally altering the landscape of in situ forest inventories [8].Three-dimensional point cloud data obtained through light detection and ranging (LiDAR) technology are currently being utilized for the precise extraction of forest characteristics, such as individual tree position, DBH, tree height, and forest biomass [7,[13][14][15][16][17][18].
Among all LiDAR sensors and platforms, the static terrestrial laser scanning (TLS) system has the highest geometric data quality at the plot level [8].It can capture highquality point clouds that provide detailed information at the millimeter level [19].Over the past decade, TLS has gained global recognition as an alternative to traditional methods for forest resource surveys [20][21][22].However, TLS lacks mobility and is unable to capture the complete profile of a tree from a single viewpoint.Additionally, it is well known that artificial reference targets are needed to combine multiple scans into a single point cloud in order to cover the entire area of interest [19,23].This inconvenience in data acquisition and processing has become a major obstacle to the widespread adoption of TLS.Recent studies have explored marker-free automated registration [20,24] to enhance the applicability of TLS in forest conditions.However, the results indicated that these methods do not achieve the same level of reliability as traditional approaches using artificial reference targets.
Mobile and portable systems greatly enhance the efficiency of data collection and have the potential to be embraced as next-generation operational tools in studying forest spatial characteristics when they achieve geometric accuracy comparable to TLS systems [8].Its platform can be a vehicle [6,25] or a person, the latter also known as personal laser scanning (PLS) [26,27] or wearable laser scanning [28].Depending on the platform they are mounted on, PLS can be further categorized into backpack laser scanning (BPLS) [29,30] or handheld laser scanning (HLS) [4,27,[31][32][33].A typical PLS system comprises a LiDAR sensor, a global navigation satellite system (GNSS) receiver, and an inertial measurement unit (IMU) [18].However, due to weak or absent GNSS signals beneath the forest canopy, there are limitations for GNSS-based PLS in positioning.The introduction of simultaneous localization and mapping (SLAM) has empowered the latest PLS systems, particularly HLS devices, with the capability to capture detailed 3D scenarios while in motion without relying on GNSS signals.By providing spatial data in a local coordinate system, SLAM software enables raw data to be quickly pre-processed and exported to various point cloud formats even on a modest laptop [4].Thus, PLS exhibits superior mobility and higher efficiency compared to TLS.This flexibility is particularly advantageous in forest inventory in challenging terrain and structurally complex forest environments [8].For example, ref. [27] reported that the data acquisition time with handheld PLS (10.96 min per plot) is 4.7 times faster than with multi-scan TLS (49.9 min per plot).Additionally, Bauwens et al. [32] conducted a comparison of survey time efficiency among HLS, TLS, and the conventional field method, reporting the survey coverage times per investigator of 50 m 2 /min, 0.85 m 2 /min, and 0.43 m 2 /min, respectively.Another advantage of HLS is the reduction in occluded areas [30].Unlike TLS, PLS uses the movement of the operator as a platform, which can minimize occlusion effects since the trajectory through the plot is equivalent to a theoretically unlimited number of scan positions [32].However, PLS systems usually have lower angular resolutions and ranging accuracy and larger beam divergences due to the integration of less expensive laser sensors in order to reduce weight and size.Consequently, more noise points are introduced, and PLS data quality is lower compared to TLS [31,34,35].
The state of the art in HLS has exhibited time efficiency and good performance in tree mapping [26][27][28][30][31][32][33]36,37].Individual tree attributes such as tree position, DBH, Forests 2024, 15, 575 3 of 21 and tree height are the primary focus of researchers.For example, ref. [27] investigated up to 2466 trees from 20 sample plots with different sizes distributed across various forest types, stand structures, and terrain characteristics in Austria.The root mean square error (RMSE) of the best DBH was 2.32 cm (12%), and the highest precision of relative bias was approximately 1%, which can be considered satisfactory for operational forest inventories (FI).However, to our best knowledge, ref. [27] performed the most comprehensive study with sufficient data.Most HLS studies often collect limited data at the tree or plot levels, which may lack strong statistical robustness.Recently, ref. [31] reviewed all studies involving HLS in FI.They reported promising results but stated that the limited number of conducted studies prohibits a definitive conclusion on the current suitability of HLS systems for operational FI.After an overview of the current status and advancements in close-range remote sensing for forest observations, ref. [8] also indicated that although promising results have been demonstrated, some challenges are still ahead.For practical applications in forests, additional studies with sufficient test data are needed to provide statistically reliable conclusions, including insights into the quantity, variability, and complexity of the test data.Furthermore, most of the studies [26,28,30,32,33,35] only consider a few tree species and do not account for the influence of tree species on the accuracy of LiDAR scanning.
On the other hand, previous studies related to FI have utilized several SLAM-based HLS devices: the ZEB1 [33], ZEB-REVO [28], ZEB-REVO-RT [26], and ZEB-HORIZON [27], all developed by GeoSLAM Ltd., Nottingham, UK.The launch dates for these devices were 2013, 2015, 2017, and 2019, respectively.Unlike the other devices, which incorporate the UTM-30LX sensor, the ZEB-HORIZON is equipped with the Velodyne VLP-16 laser scanner (Velodyne LiDAR, San Jose, CA, USA).The adoption of VLP-16 has significantly enhanced both the scanning range (up to 100 m) and acquisition rate (300,000 pts/s).Nonetheless, the accuracy of the range is moderate within a range of 3 cm, and the beam divergence is relatively large (3-5 mrad), both of which impact the localization uncertainty of individual points.As a low-cost sensor that may promote the popularization of laser scanning technology in FI, additional studies are needed to assess its performance in extracting forest parameters.
In this study, we introduced a new device called LiGrip H120, developed by Beijing Green Valley Technology Co., Ltd., Beijing, China, which integrates a Velodyne VLP-16 puck LiDAR (Velodyne LiDAR, San Jose, CA, USA) sensor and a HD camera.A relatively simple forest environment within a university campus, characterized by a clear understory, was selected as the study area.This choice served two purposes: first, it facilitated measurements and laser scanning; second, the campus harbors a variety of tree species, allowing for analysis of the impact of tree species on the device performance.
The main objectives of this study are to (1) obtain sufficient single-tree attributes in an urban area using an HLS system and a standardized workflow in data acquisition and parameter extraction; (2) investigate the accuracy of tree position, DBH, and tree height estimation with reference data; and (3) assess the impact of terrain, tree species, and other factors on the accuracy of tree parameter extraction.

Study Area and Tree Mensuration
The area selected for PLS performance assessment in this study is located on the campus of Qingdao Agricultural University (QAU), Chengyang, Qingdao, Shandong, China (Figure 1).As a coastal and mid-latitude city, Qingdao has a typical maritime climate with four distinctive seasons.The average annual temperature is 12.6 • C, with a cold winter of 0.9 • C mean temperature from December to February and a mild summer of 23.3 • C mean temperature from June to August.The average annual precipitation is about 662 mm and is mainly distributed in summer [38,39].The Chengyang campus of QAU was established in July 2001 on agricultural lands.Therefore, nearly all of the trees on the Plus, Qianxun SI, Shanghai, China) utilizing the projected China Geodetic Coordinate System 2000 (CGCS 2000).With real-time network calculations, the system can provide a positional accuracy of 1-3 cm.
The tree species were also recorded during the measurement work.Finally, a total of 2227 trees from 34 tree species were surveyed.However, many species had insufficient sample sizes, with some having only a single tree.To ensure robust statistical analysis, species with fewer than 20 samples were excluded.Consequently, 2083 trees from 13 species were retained.The summary information of the measured plots and the parameters of the tree species are shown in Tables 1 and 2.  Different from natural forests, urban trees are dispersed across various forms of vegetated areas within built-up regions, including streets, open fields in parks, affiliated areas, or residential zones [9].The characteristics of small forest patches and heterogeneous tree distribution determine that sample size, shape, and location greatly influence estimated parameters such as the number of trees per hectare, mean DBH, or tree height.Therefore, we selected certain green areas within the campus of QAU for our study, without specifically designing plot size.Instead, all trees with a DBH greater than 5 cm in these areas were surveyed, except for economic tree species that typically have multiple branches below 1.3 m.
The workflow of the methods executed in this study is shown in Figure 2. The tree measurements were conducted in December 2023.DBH was measured using a diameter tape with a precision of 1 mm.Tree height was primarily measured using a Blume Leiss (CGQ-1, Harbin Optical Instrument Factory Ltd., Harbin, China), which was developed to measure slope and tree height using trigonometric principles.Prior to measurement, the distance between the surveyor and the tree was ranged at 15 m or 20 m using a measuring tape, based on the estimated tree height.Then, the operator observed the treetop through the eyepiece of the device and triggered the measurement.Finally, the tree height was obtained by adding the operator's eye height to the value directly read from the device.For small trees below 10 m, height was measured directly by a tower ruler.Tree position was measured using a real-time kinematic (RTK) system (QianXun Xingyao X Plus, Qianxun SI, Shanghai, China) utilizing the projected China Geodetic Coordinate System 2000 (CGCS 2000).With real-time network calculations, the system can provide a positional accuracy of 1-3 cm.
The tree species were also recorded during the measurement work.Finally, a total of 2227 trees from 34 tree species were surveyed.However, many species had insufficient sample sizes, with some having only a single tree.To ensure robust statistical analysis, species with fewer than 20 samples were excluded.Consequently, 2083 trees from 13 species were retained.The summary information of the measured plots and the parameters of the tree species are shown in Tables 1 and 2.

HLS Data Acquisition
The HLS data were collected from 6 to 9 April 2023, using the LiGrip H120 from jing Green Valley Technology Co., Ltd., Beijing, China (http://www.lidar360.comchives/portfolio/ligrip-h120, accessed on 20 January 2024).This period is springtime in Q dao, and the trees have just sprouted new leaves.The LiGrip H120 system integrates the dyne VLP-16 puck LiDAR sensor and a video camera with a 360° panoramic lens.The d can capture 320,000 pts/s in a single return, and the scan range is 120 m with an accuracy cm in a field view of 280° × 360°.Based on the Simultaneous Localization and Map (SLAM) and Inertial Measurement Unit (IMU) techniques, the HLS system can automat match and calibrate data reiterations.However, the performance is influenced by weather ditions (mainly wind) and obstacles from tree canopy and low shrubs.
An app installed on a smartphone or tablet is recommended for data collection advantage is that the operator can view current scanning data, coverage, the trajecto the measured area, and other details in real time.The device was first placed on a sta ary surface or platform without high pedestrian or vehicular traffic to complete th tialization.This process takes approximately 75 s.Subsequently, the operator held th ser device to walk through the sample plot with a stable speed of 1 m/s along the designed route to collect point data and record video.Effective path planning can co all information about the trees while reducing data redundancy.In this study, the size shape of the measured forest area were not specified.We conducted data collection lowing the recommended path outlined in the equipment manual (Figure 3).If the are densely packed, the path would follow the Figure 3a; otherwise, Figure 3b wa lowed when the trees are sparse.The size of the sampled area was determined base the requirement that the collection time does not exceed 30 min.Ensuring that the ning path forms a closed loop significantly enhances the reliability and precision o

HLS Data Acquisition
The HLS data were collected from 6 to 9 April 2023, using the LiGrip H120 from Beijing Green Valley Technology Co., Ltd., Beijing, China (http://www.lidar360.com/archives/portfolio/ligrip-h120, accessed on 20 January 2024).This period is springtime in Qingdao, and the trees have just sprouted new leaves.The LiGrip H120 system integrates the Velodyne VLP-16 puck LiDAR sensor and a video camera with a 360 • panoramic lens.The device can capture 320,000 pts/s in a single return, and the scan range is 120 m with an accuracy of ±3 cm in a field view of 280 • × 360 • .Based on the Simultaneous Localization and Mapping (SLAM) and Inertial Measurement Unit (IMU) techniques, the HLS system can automatically match and calibrate data reiterations.However, the performance is influenced by weather conditions (mainly wind) and obstacles from tree canopy and low shrubs.
An app installed on a smartphone or tablet is recommended for data collection.The advantage is that the operator can view current scanning data, coverage, the trajectory of the measured area, and other details in real time.The device was first placed on a stationary surface or platform without high pedestrian or vehicular traffic to complete the initialization.This process takes approximately 75 s.Subsequently, the operator held the laser device to walk through the sample plot with a stable speed of 1 m/s along the pre-designed route to collect point data and record video.Effective path planning can collect all information about the trees while reducing data redundancy.In this study, the size and shape of the measured forest area were not specified.We conducted data collection following the recommended path outlined in the equipment manual (Figure 3).If the trees are densely packed, the path would follow the Figure 3a; otherwise, Figure 3b was followed when the trees are sparse.The size of the sampled area was determined based on the requirement that the collection time does not exceed 30 min.Ensuring that the scanning path forms a closed loop significantly enhances the reliability and precision of the data.Therefore, an additional distance of 5-10 m was traveled to ensure the accurate loop closure recognized by the program.
As the device we used lacked a GNSS module, at least three ground control points (GCP) were set up along each trajectory to compare estimated tree parameters with measured results by matching absolute tree positions.During walking, the operator aligned the 'cross' hole of the marking base with each ground control point and stabilized the device to record the GCP's local coordinates in the point cloud.Once handheld marking was completed, the positions were measured using GNSS-RTK.This method is suitable for areas with poor GNSS signals, such as dense natural forests.After all, as long as a few control points with GNSS signals are identified, it is possible to correct tree positions from the point cloud.
closure recognized by the program.
As the device we used lacked a GNSS module, at least three ground control points (GCP) were set up along each trajectory to compare estimated tree parameters with measured results by matching absolute tree positions.During walking, the operator aligned the 'cross' hole of the marking base with each ground control point and stabilized the device to record the GCP's local coordinates in the point cloud.Once handheld marking was completed, the positions were measured using GNSS-RTK.This method is suitable for areas with poor GNSS signals, such as dense natural forests.After all, as long as a few control points with GNSS signals are identified, it is possible to correct tree positions from the point cloud.

Point Cloud Processing
The raw collected data were initially fused using LiFuser-BP V1.4.2 software developed by Beijing Green Valley Technology Co., Ltd., Beijing, China.Point clouds and their corresponding trajectories were generated separately for each data acquisition from the raw data.The processed output consisted of point clouds exported in LiData files, along with the corresponding calculated trajectories and marked GCPs exported as text files.After data fusion, the LiData files and trajectories were imported into the LiDAR360 V6.0 software for further processing to extract tree parameters.Initially, the point clouds within the plots were clipped according to the trajectory lines.Subsequently, a standardized procedure with a fixed parameters setting was executed.The procedure consisted of de-noising, ground point classification, normalization based on ground points, and seedbased individual tree segmentation (Figure 2): 1. De-noising.Outliers resulting from multi-path errors or inaccuracies in the laser rangefinder during the measurement process were removed to enhance data quality.
The default parameters of a neighborhood size 10 and a standard deviation multiplier 5 were adopted.2. Ground point classification.The Improved Progressive TIN Densification (IPTD) algorithm [11] was employed for ground point classification.This algorithm initially generates a sparse triangular network using seed points and iteratively densifies it layer by layer until all ground points are classified.

Point Cloud Processing
The raw collected data were initially fused using LiFuser-BP V1.4.2 software developed by Beijing Green Valley Technology Co., Ltd., Beijing, China.Point clouds and their corresponding trajectories were generated separately for each data acquisition from the raw data.The processed output consisted of point clouds exported in LiData files, along with the corresponding calculated trajectories and marked GCPs exported as text files.Figure 4d,e provide examples of point cloud profiles from two sample plots.
After data fusion, the LiData files and trajectories were imported into the LiDAR360 V6.0 software for further processing to extract tree parameters.Initially, the point clouds within the plots were clipped according to the trajectory lines.Subsequently, a standardized procedure with a fixed parameters setting was executed.The procedure consisted of denoising, ground point classification, normalization based on ground points, and seed-based individual tree segmentation (Figure 2): 1.
De-noising.Outliers resulting from multi-path errors or inaccuracies in the laser rangefinder during the measurement process were removed to enhance data quality.
The default parameters of a neighborhood size 10 and a standard deviation multiplier 5 were adopted.

2.
Ground point classification.The Improved Progressive TIN Densification (IPTD) algorithm [11] was employed for ground point classification.This algorithm initially generates a sparse triangular network using seed points and iteratively densifies it layer by layer until all ground points are classified.

3.
Normalization based on ground points.The elevation value Z of the nearest ground point was subtracted from each point in the procedure of normalization to eliminate the influence of terrain undulations on the elevation values of point cloud data.

4.
Seed-based individual tree segmentation.Parameters such as DBH and tree height were batch-extracted using the 'TLS Seed Point Editor' tool under 'TLS Forest', based on normalized ground points.DBH was determined by cylinder fitting within the slice at a height of 1.2-1.4m.Subsequently, each fitted DBH and tree height was inspected for anomalies, as well as potential commission or omission errors.In case of fitting errors, individual tree fitting would be carried out again.Finally, the fitted results were used as seed points for individual tree segmentation using the 'Point Cloud Segmentation from Seed Points' tool.
The tree positions derived from the point clouds refer to the center point coordinate of the trunk at a height of 1.3 m.They were registered to the CGCS 2000 coordinate system using the affine transformation method based on GCPs to match the positions of measured trees.This process was operated through a spatial adjustment tool in the ArcGIS 10.4 software (Environmental Systems Research Institute, Redlands, CA, USA). Figure 5 gives an example of registering the tree positions and trajectories using the affine transformation method.The scanning time for each plot was obtained from the trajectory's log file.Additionally, the trajectory files were used to calculate the trajectory length and area of each plot using the ArcGIS 10.4 software.
of fitting errors, individual tree fitting would be carried out again.Finally, the fitted results were used as seed points for individual tree segmentation using the 'Point Cloud Segmentation from Seed Points' tool.
The tree positions derived from the point clouds refer to the center point coordinate of the trunk at a height of 1.3 m.They were registered to the CGCS 2000 coordinate system using the affine transformation method based on GCPs to match the positions of measured trees.This process was operated through a spatial adjustment tool in the ArcGIS 10.4 software (Environmental Systems Research Institute, Redlands, CA, USA). Figure 5 gives an example of registering the tree positions and trajectories using the affine transformation method.The scanning time for each plot was obtained from the trajectory's log file.Additionally, the trajectory files were used to calculate the trajectory length and area of each plot using the ArcGIS 10.4 software.

Evaluation
In order to assess the accuracy of DBH and tree height derived from HLS point clouds for different species and plots, the total explained variance (R 2 ), root mean squared error (RMSE), relative RMSE (RMSE%), absolute bias (Bias), and relative bias (Bias%) were calculated using Equations ( 1)-( 5):

Evaluation
In order to assess the accuracy of DBH and tree height derived from HLS point clouds for different species and plots, the total explained variance (R 2 ), root mean squared error (RMSE), relative RMSE (RMSE%), absolute bias (Bias), and relative bias (Bias%) were calculated using Equations ( 1)-( 5): Forests 2024, 15, 575 where n is the number of individuals of a tree species, plot, or the whole number of detected trees; y i is the field-measured DBH or tree height of tree i; ŷi indicate the estimated DBH or tree height value of tree i from HLS point clouds; and y denotes the mean DBH or tree height value of n.

Statistical Analysis
In order to assess the factors affecting the performance of HLS in tree mapping, linear mixed-effects model (LME) analysis was applied (Equation ( 6)).
where Y i is the i-th response variable of RMSE, RMSE%, Bias, or Bias% for DBH and H; area is the scanning area of each plot (m 2 ); tree density is the number of trees per hectare; trajectory length is the total length of trajectory for each plot (m); plant type indicates whether the tree species is coniferous or broad-leaved; terrain indicates the topography of the plot; and ε I is the error.We considered area, tree density, trajectory length, plant type, and terrain as fixed factors and the plot as a random factor to address the repeated measurements within the plot.
The statistical analysis was conducted with R 4.2.3 [40].The lmer() function of the lme4 package [41] was used to execute the linear mixed-effects models.We tested all linear mixed-effects models with a random intercept to assess the significance of candidate factors on HLS performance.If the interaction effects were deemed non-significant, they were excluded to simplify the statistical model.The model with the lowest Akaike's Information Criteria (AIC) was selected, and the estimated p-values were computed.The r.squaredGLMM() function of the MuMIn package [42] was used to compare the R 2 in models.The significance threshold for all tests was set at 5%.

Tree Position
The tree positions obtained through the registration of HLS using ground control points showed high consistency with those obtained through RTK measurements, with a maximum error of 0.39 m shown in the box plot in Figure 6.Apart from some outliers, the positional accuracy of trees from different plots varies relatively little, mostly falling within the range of 0.08 to 0.2 m.

Estimation of DBH
The comparison between estimated DBH using HLS and measured DBH is shown in Figure 7.The results showed high precision, with an overall interpretive capability (R 2 ) of 0.98 and an RMSE of only 2.06 cm, which was about 6.89% of the quadratic mean diameter.Additionally, a slight underestimation was observed for the total of 2083 trees, with a bias of 0.62 cm (2.07%).Further analysis was conducted for coniferous and broad-leaved trees.In general, the HLS-derived DBH accuracy was better for broad-leaved trees than for coniferous trees, with the RMSE being 1.86 cm (6.01%) and 2.54 cm (9.46%), respectively.However, the bias analysis revealed that the broad-leaved trees were subject to greater underestimation compared to the coniferous trees.
For visual interpretation, specifically, 13 tree species were plotted in Figure 8. Generally, high accuracy was achieved for all tree species regarding different ranges of DBH.Except for Pinus thunbergii, R 2 ranged from 0.90 to 0.98.Additionally, RMSE varied between 1.26 and 2.79 cm, which was not consistent with its percent due to differences in the DBH range.Tree species with relative RMSE greater than 10% are typically those with smaller DBH.In terms of bias and relative bias, the HLS technique tended to overestimate DBH for Koelreuteria paniculata, Prunus cerasifera, Robinia pseudoacacia, and Pinus thunbergii, while it slightly underestimated other tree species.Notably, Prunus cerasifera exhibits a remarkably high relative bias of 7.35%.In contrast, Koelreuteria paniculata was one of the most accurately estimated species, with a relative bias of only 0.4%.For the majority of tree species, the relative bias falls between 1% and 5%.

Estimation of DBH
The comparison between estimated DBH using HLS and measured DBH is shown in Figure 7.The results showed high precision, with an overall interpretive capability (R 2 ) of 0.98 and an RMSE of only 2.06 cm, which was about 6.89% of the quadratic mean diameter.Additionally, a slight underestimation was observed for the total of 2083 trees, with a bias of 0.62 cm (2.07%).Further analysis was conducted for coniferous and broad-leaved trees.In general, the HLS-derived DBH accuracy was better for broad-leaved trees than for coniferous trees, with the RMSE being 1.86 cm (6.01%) and 2.54 cm (9.46%), respectively.However, the bias analysis revealed that the broad-leaved trees were subject to greater underestimation compared to the coniferous trees.
For visual interpretation, specifically, 13 tree species were plotted in Figure 8. Generally, high accuracy was achieved for all tree species regarding different ranges of DBH.Except for Pinus thunbergii, R 2 ranged from 0.90 to 0.98.Additionally, RMSE varied between 1.26 and 2.79 cm, which was not consistent with its percent due to differences in the DBH range.Tree species with relative RMSE greater than 10% are typically those with smaller DBH.In terms of bias and relative bias, the HLS technique tended to overestimate DBH for Koelreuteria paniculata, Prunus cerasifera, Robinia pseudoacacia, and Pinus thunbergii, while it slightly underestimated other tree species.Notably, Prunus cerasifera exhibits a remarkably high relative bias of 7.35%.In contrast, Koelreuteria paniculata was one of the most accurately estimated species, with a relative bias of only 0.4%.For the majority of tree species, the relative bias falls between 1% and 5%.

Estimation of Tree Height
Similar to DBH, the relationship between HLS-derived and measured tree height also resulted in a high R 2 value of 0.95 for all species of trees.As depicted in Figure 7, the RMSE reached 1.16 m, which was less than 9% of the mean height.Meanwhile, a slight underestimation in H was also observed, with a bias of 0.27 m (4.67%).Similar results were exhibited in broad-leaved and coniferous trees, with RMSE values of 1.13 m (8.37%) and 1.24 m (10.05%), respectively, and biases of 0.3 m (2.19%) and 0.21 (1.67%), respectively.

Estimation of Tree Height
Similar to DBH, the relationship between HLS-derived and measured tree height also resulted in a high R 2 value of 0.95 for all species of trees.As depicted in Figure 7, the RMSE reached 1.16 m, which was less than 9% of the mean height.Meanwhile, a slight underestimation in H was also observed, with a bias of 0.27 m (4.67%).Similar results were exhibited in broad-leaved and coniferous trees, with RMSE values of 1.13 m (8.37%) and 1.24 m (10.05%), respectively, and biases of 0.3 m (2.19%) and 0.21 (1.67%), respectively.
Regarding each specific tree species, the differences among them were significant (Figure 9).Firstly, the height prediction accuracy for individual tree species was much lower, ranging from a maximum of 0.92 to a minimum of 0.69, representing a 25% lower accuracy than for all trees (R 2 = 0.95).Secondly, the RMSE values ranged from 0.61 to 1.58 m, ranging from 4% to 21.4%, indicating a significant difference between species.Finally, the heights of three species (Robinia pseudoacacia, Yulia deudata, and Metasequoia glyphostroboides) tended to overestimate, while for other species, the trend was the opposite.These findings also suggest that the precision of DBH values was higher than that of tree height values.
Forests 2023, 14, x FOR PEER REVIEW 14 of 21 Regarding each specific tree species, the differences among them were significant (Figure 9).Firstly, the height prediction accuracy for individual tree species was much lower, ranging from a maximum of 0.92 to a minimum of 0.69, representing a 25% lower accuracy than for all trees (R 2 = 0.95).Secondly, the RMSE values ranged from 0.61 to 1.58 m, ranging from 4% to 21.4%, indicating a significant difference between species.Finally, the heights of three species (Robinia pseudoacacia, Yulia deudata, and Metasequoia glyphostroboides) tended to overestimate, while for other species, the trend was the opposite.These findings also suggest that the precision of DBH values was higher than that of tree height values.

9.
Comparison of the performance of HLS in estimating tree height (H) between 13 species.The red lines represent the linear regressions between measured H and HLS-derived H, while the dashed lines are the 1:1.

LME Results
In general, handheld laser scanning tends to exhibit higher accuracy in extracting DBH and tree height (H) for broad-leaved trees and trees situated in flat terrain (Figure 10).The linear mixed-effects model (LME) analysis revealed that the factors of tree density, plant type, and terrain significantly influenced the accuracy of HLS-derived DBH, whereas only the plant type significantly affected the accuracy of H (Table 3).However, in the final model, the significant influencing factors for DBH accuracy were reduced to plant type and terrain, while for accuracy of H, in addition to plant type, terrain also had a significant impact on RMSE% of H (Table 4).

LME Results
In general, handheld laser scanning tends to exhibit higher accuracy in extracting DBH and tree height (H) for broad-leaved trees and trees situated in flat terrain (Figure 10).The linear mixed-effects model (LME) analysis revealed that the factors of tree density, plant type, and terrain significantly influenced the accuracy of HLS-derived DBH, whereas only the plant type significantly affected the accuracy of H (Table 3).However, in the final model, the significant influencing factors for DBH accuracy were reduced to plant type and terrain, while for accuracy of H, in addition to plant type, terrain also had a significant impact on RMSE% of H (Table 4).

Discussion
The registration of tree positions using HLS and GCPs achieved satisfactory accuracy in this study, apart from a few tilted trees, with stem positions at 1.3 m from point clouds compared to the measured base position.There are typically two approaches to obtaining absolute positions in FI by HLS.One approach is to use tightly coupled GNSS receivers and IMUs to create positioning subsystems.These subsystems capture platform movements and sensor orientation data, providing the system's position and sensor orientation at discrete time intervals.This information enables direct georeferencing of collected data [8,43].With a tactical-grade GNSS IMU, absolute accuracy can be reached at 0.2-0.7 m after postprocessed positioning in boreal forests [25].However, this accuracy may decrease with low-cost IMUs due to high positional drift in conditions with low GNSS visibility and increased angular uncertainty.The other approach involves registering point cloud data using GCPs, with total stations commonly used to provide high-precision coordinates [44].These GCPs are used to georeference HLS data to the coordinate system defined by the GCPs during the registration of point clouds.Tupinambá-Simões et al. [18] compared the rigid and non-rigid modes for generating the 3D point clouds in the ZEB HLS system.Under the rigid option, the point cloud is rotated and adjusted, maintaining the scale factor, and GCP locations are matched without altering the relative positions of points.In contrast, the non-rigid mode may involve moving, rotating, or changing the scale between points, potentially altering their relative positions.They reported better performance of the non-rigid mode in the estimation of DBH.This is attributed to the use of total stations to obtain high-precision GCPs.However, setting up and using the total station is timeconsuming, and it requires a known control point to initiate the measurement, which further complicates the data collection process.In our study, GCPs were not used for point cloud generating but rather for the correction of tree positions derived from point cloud post-processing.GCPs do not affect the relative positions of points, so their accuracy does not affect the extraction of DBH or tree height.Therefore, it is unnecessary to measure GCPs using a total station.We utilized real-time kinematics (RTKs) for measuring the GCPs and achieved a precise match between the stem positions in the point clouds and the measured tree positions.
In terms of DBH estimation, our study indicated high accuracy using HLS in an urban area, with RMSE ranging from 1.26 to 2.79 cm.Studies conducted in urbanized areas, including city parks, urban streets, plantation sites, and other similar environments, also demonstrated good performance in estimating DBH with HLS [26,[45][46][47].This was largely due to the clean and unobstructed understory within urban forests [28].While it may seem more convenient for an operator to move through flat plots using a laser scanning device, various obstacles such as pedestrians, moving or parked cars, and landscape shrubs can still hinder the accuracy of LiDAR scanning.For instance, in a study by [48], the RMSE value of DBH for park trees was 8.95 cm, attributed to irregular trunk shapes and incomplete scanning data.In our study, factors affecting DBH accuracy mainly included inclined tree trunks, cross-section irregularities (e.g., forked, scarred, non-circular), and occlusion, which were more common in coniferous trees.The image in Figure 4c depicts an inclined tree of Pinus thunbergii, which decreased the fitted accuracy of DBH.As a result, broad-leaved trees exhibited higher accuracy values compared to coniferous trees.This observation aligns with findings from previous studies conducted in challenging forest environments [4,34].
Meanwhile, most previous studies were conducted with limited sample sizes, lacking sufficient robustness in their results, especially in urban areas.With a total of 2466 trees from 20 sample plots of various sizes distributed across diverse forest types, stand structures, and terrain characteristics in Austria, ref. [27] reported an RMSE of 2.32 cm (12%) in DBH, with a relative bias of approximately 1%.Our study included 2083 trees from 34 plots in an urban area, representing a sufficient sample size to evaluate the capability of HLS in extracting tree parameters.Hyyppä et al. [34] demonstrated that the RMSE for DBH derived from LiDAR data exhibited significant variability depending on the complexity of the forest plots.Consequently, more studies should be carried out in different circumstances to provide statistically reliable conclusions [8].
In comparison to DBH, tree heights were relatively inaccurate with lower R 2 and higher RMSE values, especially when considering individual tree species.Similarly, in a study by [28], tree heights were compared using a TLS and a portable laser scanner (PLS) in two plots with different tree densities in Spain.The results showed that the RMSE difference between the two devices was 1.34 m in the plot with lower density and 9.44 m in the plot with higher density.These variations were attributed to differences between the devices and the environmental conditions of the plots.Challenging environments not only affect the performance of HLS devices but also influence the accuracy of TLS.For example, ref. [49] classified three forests into categories of Easy, Moderate, and Hard, based on tree density and vegetation cover to evaluate the accuracy of measuring tree parameters using a TLS multi-scan approach.They reported lower average accuracy compared to our study, with the RMSE and RMSE% values for height measurements ranging from 2.4 to 4.5 m and 12 to 23% for Easy conditions, and 4.0 to 7.7 m and 28 to 57% for Hard conditions.
In general, tree height was underestimated in this study, which is consistent with previous studies [34,[50][51][52].For instance, ref. [35] reported an underestimation of tree heights with a bias of −4.61 m and an RMSE of 2.15 m compared to field reference data using the HLS method.On the contrary, HLS tree heights were slightly overestimated in the study of [53].However, it is essential to consider that indirect field measurements of tree height are complex [54].Their accuracy can be influenced by various potential sources of error, such as forest structure and complexity, tree species and crown shape, leaning trees, measuring distance, tree height, instruments, and human errors.Thus, errors in tree height accuracy from HLS encompass errors associated with both HLS itself and those related to the reference field-measured tree heights.
Most of the previous studies using LiDAR for extracting tree parameters had limited sample sizes, allowing only for qualitative analysis of the impacts of forest structure, topography, etc., on the accuracy of extraction.In our study, 34 sample plots and linear mixed-effects models were employed to quantitatively analyze the effects of forest structure and other factors on the accuracy of DBH and tree height extraction.The result revealed that the plant type (broad-leaved or coniferous) and terrain were two significant factors leading to decreased accuracy of DBH and H.As stated in a study by [55], detecting spruce trees via terrestrial-based remote sensing posed considerable challenges.The dense branching and limited visibility of the lower part, especially at breast height, make it difficult for sensors to capture their entire stems.Similarly, ref. [4] analyzed the influence of stand characteristics (such as the number of trees, stand basal area, dominant height, understorey cover, slope, etc.) on HLS data accuracy using datasets from 39 sample plots.Concerning site conditions, absolute bias percentage values of arithmetic mean diameter, stand basal area, and stand volume showed a positive correlation with both ground slope and understorey cover.The errors in estimating arithmetic mean diameter and stand basal area significantly increased with understorey coverage, while errors in stand volume estimations significantly increased with slope gradient.
Finally, handheld LiDAR devices are undergoing rapid updates.During our study, a new-generation product, the LiGrip H300, part of the LiGrip handheld series from Beijing Green Valley Technology Co., Ltd., was released.Compared to its predecessor, the LiGrip H120, it boasted an increased maximum measurement range from 100 m to 300 m, a scanning frequency from 320,000 to 640,000 pts/s, and enhanced LiDAR accuracy to 1 cm.Recently, consumer-grade laser sensors have been integrated into smartphones, offering more affordable scanning technology.These affordable sensors are anticipated to proliferate in the future, offering more handheld laser scanning technology that is user-friendly [56].Therefore, we believe that mobile and portable laser scanning systems have the potential to be embraced as next-generation operational tools.

Conclusions
This study evaluated the applicability of an HLS system for the examination of urban forest resource inventories and its performance in mapping individual trees.Tree position, DBH, and height were estimated by HLS from 34 different plots in QAU, and the results were compared with field measurements.
The registration of tree position using GCPs showed high accuracy with errors under 0.4 m, except for some outliers.The extraction accuracy of DBH for all trees and individual tree species was consistently high, with a total RMSE of 2.06 cm (6.89%) and a bias of 0.62 cm (2.07%).Broad-leaved trees exhibited better performance than coniferous trees,

Forests
planted or transplanted during the last twenty years, resulting in a heterogeneous mix in terms of species, age, and height.

Figure 1 .
Figure 1.The location of the study area in Qingdao City, Shandong, China (a), and the trajectories of handheld mobile laser scanning along with the positions of trees for each plot (b).The numbers indicate the ID of the plots.

Figure 1 .
Figure 1.The location of the study area in Qingdao City, Shandong, China (a), and the trajectories of handheld mobile laser scanning along with the positions of trees for each plot (b).The numbers indicate the ID of the plots.

Figure 2 .
Figure 2. The workflow executed in this study.

Figure 2 .
Figure 2. The workflow executed in this study.

Figure 3 .
Figure 3. Path planning for dense trees (a) and sparse trees (b).
Figure 4d,e provide examples of point cloud profiles from two sample plots.

Figure 3 .
Figure 3. Path planning for dense trees (a) and sparse trees (b).

Figure 4 .
Figure 4. Point clouds of (a) a Platanus acerifolia tree, (b) a straight-trunked Pinus thunbergii tree, (c) a curved-trunked Pinus thunbergii tree, (d) a profile from flat ground, and (e) a profile from a hilly area.

Figure 4 .
Figure 4. Point clouds of (a) a Platanus acerifolia tree, (b) a straight-trunked Pinus thunbergii tree, (c) a curved-trunked Pinus thunbergii tree, (d) a profile from flat ground, and (e) a profile from a hilly area.

ForestsForests 21 Figure 5 .
Figure 5.An example of registering tree positions and trajectory using the affine transformation method.GCP indicates the ground control point, while ligrip refers to its corresponding local coordinates marked with the HLS device in the point cloud.

Figure 5 .
Figure 5.An example of registering tree positions and trajectory using the affine transformation method.GCP indicates the ground control point, while ligrip refers to its corresponding local coordinates marked with the HLS device in the point cloud.

Figure 6 .
Figure 6.The accuracy of registration for tree position.

Figure 6 .
Figure 6.The accuracy of registration for tree position.

Figure 7 .
Figure 7. HLS-derived versus observed DBH (diameter at breast height) and H (tree height) for all species, broad-leaved, and coniferous trees.The red lines represent the linear regressions, while the dashed lines are the 1:1.

Figure 7 .
Figure 7. HLS-derived versus observed DBH (diameter at breast height) and H (tree height) for all species, broad-leaved, and coniferous trees.The red lines represent the linear regressions, while the dashed lines are the 1:1.

Forests 2024 , 21 Figure 8 .
Figure 8.Comparison of the performance of HLS in estimating DBH (diameter at breast height) among 13 tree species.The red lines represent the linear regressions, while the dashed lines are the 1:1.

Figure 8 .
Figure 8.Comparison of the performance of HLS in estimating DBH (diameter at breast height) among 13 tree species.The red lines represent the linear regressions, while the dashed lines are the 1:1.

Figure 9 .
Figure 9.Comparison of the performance of HLS in estimating tree height (H) between 13 species.The red lines represent the linear regressions between measured H and HLS-derived H, while the dashed lines are the 1:1.

Figure 10 .
Figure 10.Boxplot illustrating the accuracy of DBH (diameter at breast height) and H (tree height) for different plant types (broad-leaved or coniferous) and terrains (flat or hilly) of the plot.Flat.B and Hilly.C denote broadleaved tree species on flat ground and conifer tree species on hilly terrain, respectively.The circular dots in the figure represent outliers.

Figure 10 .
Figure 10.Boxplot illustrating the accuracy of DBH (diameter at breast height) and H (tree height) for different plant types (broad-leaved or coniferous) and terrains (flat or hilly) of the plot.Flat.B and Hilly.C denote broadleaved tree species on flat ground and conifer tree species on hilly terrain, respectively.The circular dots in the figure represent outliers.
is the number of trees, DBH indicates the diameter at breast height, mean DBH refers to the quadratic mean of DBH values, and SD indicates the standard deviation. n

Table 1 .
The summary information of the 34 measured plots on the campus of QAU (Qingdao Agricultural University).
n is the number of trees, DBH indicates diameter at breast height, H is tree height, and mean DBH refers to the quadratic mean of DBH values.

Table 2 .
Summary statistics of tree species collected from 34 plots.
n is the number of trees, DBH indicates the diameter at breast height, mean DBH refers to the quadratic mean of DBH values, and SD indicates the standard deviation.

Table 3 .
The coefficients (Coef.),standard errors (Std.Error), and p-values of linear mixed-effects models between the accuracy (RMSE, RMSE%, Bias, Bias% for DBH and H) and its candidate influence factors (tree density, trajectory length, scanning time, terrain of plot, and plant type).

Table 4 .
The coefficients (Coef.),standard errors (Std.Error), and p-values of final models between the accuracy (RMSE, RMSE%, Bias, Bias% for DBH and H) and its candidate influence factors (tree density, trajectory length, scanning time, terrain of plot, and plant type).