欢迎访问林业科学,今天是

林业科学 ›› 2023, Vol. 59 ›› Issue (9): 23-33.doi: 10.11707/j.1001-7488.LYKX20210831

• • 上一篇    下一篇

点云密度对机载激光雷达大区域亚热带森林参数估测精度的影响

周相贝1,2, 李春干1, 代华兵3, 余铸1,3, 李振3, 苏凯1   

  1. 1. 广西大学林学院 南宁530004;
    2. 广西自然资源职业技术学院 扶绥 532199;
    3. 广西壮族自治区林业勘测设计院 南宁 530011
  • 收稿日期:2021-11-16 修回日期:2022-01-10 发布日期:2023-10-28
  • 通讯作者: 李春干
  • 基金资助:
    广西林业科技推广示范项目(GL2020KT02);广西壮族自治区林业勘测设计院科研业务费项目(GXLKYKJ201601)

Effects of Point Cloud Density on the Estimation Accuracy of Large-Area Subtropical Forest Inventory Attributes Using Airborne LiDAR Data

Zhou Xiangbei1,2, Li Chungan1, Dai Huabing3, Yu Zhu1,3, Li Zhen3, Su Kai1   

  1. 1. Forestry College of Guangxi University Nanning 530004;
    2. Guangxi Natural Resources Vocational and Technical College Fusui 532199;
    3. Guangxi Forest Inventory and Planning Institute Nanning 530011
  • Received:2021-11-16 Revised:2022-01-10 Published:2023-10-28

摘要: 目的 点云密度是影响机载激光雷达数据获取和预处理成本的关键因素,探明点云密度对森林参数估测精度的影响,为机载激光雷达大区域森林调查监测应用技术方案的优化提供参考依据。方法 基于我国广西一个亚热带山地丘陵区域获取的机载激光雷达和样地数据,通过系统稀疏方法,将全密度点云(4.35点·m-2)分别稀疏至4.0、3.5、3.0、2.5、2.0、1.5、1.0、0.5、0.2和0.1点·m-2,得到11个样地尺度的点云数据集,包括1个全密度和10个稀疏密度点云数据集;应用配对样本t检验方法,分析4种森林类型(杉木林、松树林、桉树林和阔叶林)中稀疏密度点云和全密度点云之间12个激光雷达变量的差异;通过变量和结构固定的多元乘幂模型式,分别采用不同密度点云数据集对林分蓄积量(VOL)和断面积(BA)进行估测,比较模型优度统计指标决定系数(R2)、相对均方根误差(rRMSE)和平均预估误差(MPE)的差异,并应用t检验方法分析稀疏密度点云VOL和BA估测值均值和全密度点云相应估测值均值的差异。结果 1) 点云密度较低时,稀疏密度点云分位数高度(ph25、ph50和ph75)的均值与全密度点云相应变量的均值存在显著性差异,但不同森林类型、不同变量出现显著性差异时的点云密度不同,各森林类型中稀疏密度点云平均高(Hmean)和点云高变动系数(Hcv)的均值与全密度点云相应变量的均值基本不存在显著性差异,但点云最大高(Hmax)的均值存在显著性差异;2) 各森林类型中,稀疏密度点云冠层覆盖度(CC)和中下层分位数密度(dh25)的均值与全密度点云相应变量的均值差异不显著(阔叶林dh25除外),但中上层分位数密度(dh50和dh75)存在显著性差异;3) 各森林类型中,稀疏密度点云平均叶面积密度(LADmean)的均值与全密度点云LADmean的均值存在显著性差异,当点云密度较低时,稀疏密度点云叶面积密度变动系数(LADcv)的均值与全密度点云LADcv的均值存在显著性差异;4) 各森林类型中,不同密度点云VOL和BA估测值差异很小,且均不存在显著性差异,但随点云密度降低,杉木林、松树林和桉树林VOL和BA估测模型的R2缓慢逐渐减小,rRMSE和MPE缓慢逐渐增大,森林参数估测精度逐渐降低,阔叶林VOL和BA估测模型的R2、rRMSE和MPE受点云密度变化影响不大。结论 点云密度降低导致激光雷达变量标准差增大是造成森林参数估测模型精度降低的主要原因,在实际机载激光雷达森林资源调查监测应用中,点云密度以大于0.5点·m−2为宜。

关键词: 机载激光雷达, LiDAR变量, 林分蓄积量, 断面积, 模型

Abstract: Objective Point cloud density is a critical factor affecting the cost of airborne LiDAR data acquisition and pre-processing. Therefore, exploring the influence of point cloud density on the estimation accuracy of forest inventory attributes can provide a reference for optimizing technical schemes for airborne LiDAR-based large-area forest inventory and monitoring. Method In this study, we used airborne LiDAR data and field plot data collected in a subtropical mountainous and hilly region in Guangxi, China. Firstly, the original point clouds with a density of 4.35 points·m−2 were reduced to 4.0, 3.5, 3.0, 2.5, 2.0, 1.5, 1.0, 0.5, 0.2, and 0.1 points·m−2 using a systematic thinning method, respectively, resulting in 11 plot-level point cloud datasets, including one full-density point cloud dataset and ten reduced-density point cloud datasets. Secondly, a paired sample t-test was used to analyze the differences in 12 LiDAR-derived metrics between reduced-density point clouds and full-density point clouds in four forest types (Chinese fir, pine, eucalyptus, and broad-leaved). Thirdly, using a multiplicative power model formulation with fixed variables and stable structure, the stand volume (VOL) and basal area (BA) were estimated using various density datasets of point clouds, respectively, and their goodness-of-fit statistics, including coefficient of determination (R2), relative root square error (rRMSE), and mean prediction error (MPE), were compared. Finally, a t-test was used to analyze the differences in the means of the estimates between the reduced-density point clouds and full-density point clouds. Result 1) When the point cloud density was low, the means of the 25th, 50th, and 75th height percentiles (ph25, ph50, and ph75) of the reduced-density point clouds showed statistically significant differences from those of the corresponding variables of the full-density point clouds. However, when statistically significant differences were found for different variables in various forest types, the point cloud densities differed. There were no statistically significant differences in the means of mean point cloud height (Hmean) and coefficient of variation of point cloud height distribution (Hcv) between the reduced-density point clouds and full-density point clouds in all forest types, but there were statistically significant differences in the means of maximum height (Hmax) of point clouds between the reduced-density point clouds and full-density point clouds for all forest types. 2) The means of canopy cover (CC) and 25th density percentile (dh25) of the reduced-density point clouds were not statistically significantly different from those of the corresponding variables of the full-density point clouds for all forest types (except dh25 for broadleaf forests), but statistically significant differences existed for the 50th and 75th density percentiles (dh50 and dh75). 3) The means of the mean leaf area density (LADmean) of reduced-density point clouds were statistically significantly different from those of the LADmean of full-density point clouds in all forest types, and while the means of the coefficient of variation of leaf area density (LADcv) of reduced-density point clouds were significantly different from those of the LADcv of full-density point clouds when point cloud density was low. 4) The differences in the estimates of VOL and BA for different density point clouds were small among the forest types, and none of the estimates were statistically significantly different from each other. However, as the density of point clouds decreased, the R2 of the estimation models for VOL and BA for fir, pine, and eucalyptus forests slowly decreased, and the rRMSE and MPE slowly increased, indicating that the estimation accuracy of forest inventory attributes gradually decreased. The R2, rRMSE, and MPE of the estimation models for VOL and BA for the broad-leaved forests were not obviously affected by the change in point cloud density. Conclusion The decrease in the density of point clouds leads to an increase in the standard deviation of the LiDAR-derived metrics, which is the main reason for the decrease in the estimation accuracy of forest inventory attributes. In the operational forest resources investigation and monitoring, the airborne LiDAR point cloud density should be greater than 0.5 points·m−2.

Key words: airborne LiDAR, LiDAR-derived metrics, stand volume, basal area, model

中图分类号: