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

林业科学 ›› 2021, Vol. 57 ›› Issue (8): 68-81.doi: 10.11707/j.1001-7488.20210807

• 论文与研究报告 • 上一篇    下一篇

基于机载激光雷达点云和随机森林算法的森林蓄积量估测

孙忠秋1,高金萍1,吴发云1,高显连1,胡杨2,*,高剑新1   

  1. 1. 国家林业和草原局调查规划设计院 北京 100714
    2. 宁夏大学生态环境学院 西北土地退化与生态恢复国家重点实验室培育基地 西北退化生态系统恢复与重建教育部重点实验室 银川 750021
  • 收稿日期:2020-03-03 出版日期:2021-08-25 发布日期:2021-09-30
  • 通讯作者: 胡杨
  • 基金资助:
    宁夏回族自治区自然科学基金项目(2021AAC03017);林草行业研究专项"陆地碳卫星"(2020-21-93*)

Estimating Forest Stock Volume via Small-Footprint LiDAR Point Cloud Data and Random Forest Algorithm

Zhongqiu Sun1,Jinping Gao1,Fayun Wu1,Xianlian Gao1,Yang Hu2,*,Jianxin Gao1   

  1. 1. Academy of Inventory and Planning, National Forestry and Grassland Administration Beijing 100714
    2. School of Ecology and Environment, Ningxia University Breeding Base for State Key Laboratory of Land Degradation and Ecological Restoration in Northwest China Key Laboratory for Restoration and Reconstruction of Degraded Ecosystem in Northwest China of Ministry of Education Yinchuan 750021
  • Received:2020-03-03 Online:2021-08-25 Published:2021-09-30
  • Contact: Yang Hu

摘要:

目的: 基于机载激光雷达点云数据提取的森林高度参数和郁闭度,结合分层地面样地调查数据,采用随机森林算法构建森林蓄积量估测模型,分析机载激光雷达点云数据在森林蓄积量反演方面的潜力,为森林蓄积量高效准确估测提供方法依据。方法: 以直径30 m的地面样圆离散点云数据为数据源,经数据校准等预处理后,利用LiDAR360软件提取森林高度参数(最大高、平均高等)和郁闭度,并将数据随机分成训练数据(70%)和验证数据(30%)。采用随机森林算法构建森林蓄积量估测模型,对仅用高度参数建模以及联合高度参数和郁闭度建模结果进行比较;同时运用R软件VSURF工具包筛选建模变量,对筛选后变量的建模结果进行分析。结果: 仅用高度参数建模的估测精度为R2=0.75、RMSE=40.07 m3·hm-2、MAE=29.21 m3·hm-2、MRE=49.40%,联合高度参数和郁闭度建模的估测精度为R2 =0.79、RMSE=36.23 m3·hm-2、MAE=26.16 m3·hm-2、MRE=38.35%。通过变量筛选,建模参数从24个减少至7个,可极大提高运算效率,同时R2未变化,RMSE从36.23 m3·hm-2升至36.50 m3·hm-2,rRMSE从31.92%升至32.97%,MAE从26.16 m3·hm-2降至26.08 m3·hm-2,MRE从38.35%降至38.05%。结论: 机载激光雷达点云数据可以提取森林的垂直结构信息(高度参数)和水平结构信息(郁闭度),具备三维结构参数提取能力。采用随机森林算法,增加林分郁闭度信息可显著提高森林蓄积量估测精度。通过变量筛选,虽然能够降低参数数量,但对模型精度具有一定影响,在建模精度要求较高的情况下,建议使用全变量进行蓄积量估测;而在数据量较大的情况下,建议使用筛选变量进行蓄积量估测。基于机载激光雷达点云数据估测森林蓄积量显著优于光学遥感数据,可为森林蓄积量高效准确估测提供方法依据,能够满足大范围森林蓄积量快速反演需求。

关键词: 森林蓄积量, 小光斑激光雷达, 随机森林算法, 变量筛选

Abstract:

Objective: The forest volume estimation model was constructed by using random forest algorithm based on the forest height parameters and crown density extracted from airborne LiDAR point cloud data, combined with stratified ground sample plot survey data. The potential of airborne LiDAR data in forest volume inversion was analyzed, so as to provide method basis for efficient and accurate estimation of forest volume. Method: The ground sample circle with a diameter of 30 m was set up, and the point cloud data of discrete sample circle was used as the data source. After data calibration and other preprocessing, the forest height parameters(maximum height, mean height) and crown density were extracted by LiDAR360 software. The research data were categorized as training and test datasets. The random forest algorithm was used for modeling using only height parameters as well as a combination of height parameters and crown density, the results thus obtained were analyzed. The VSURF package in R software was used to select modeling variables, and the selected variables were evaluated. Result: In the test phase, when using height parameters alone, modeling accuracy of R2=0.75, RMSE=40.07 m3·hm-2, MAE=29.21 m3·hm-2, and MRE=49.40% were obtained, respectively. In contrast, when combining height parameters with crown density, modeling accuracy of R2=0.79, RMSE=36.23 m3·hm-2, MAE=26.16 m3·hm-2, and MRE=38.35% were obtained, respectively. This indicated that crown density was an important parameter and that using height parameters alone was insufficient for estimating forest volumes via airborne LiDAR point cloud data. Using variable selection, 24 effective parameters were reduced to 7. Although R2 remained relatively unchanged, the RMSE and rRMSE increased from 36.23 m3·hm-2 to 36.50 m3·hm-2, and from 31.92% to 32.97%, respectively. Whereas MAE and MRE decreased from 26.16 m3·hm-2 to 26.08 m3·hm-2, and from 38.35% to 38.05%, respectively. Conclusion: Airborne LiDAR point cloud data could extract not only the vertical structure information of forest(height variable), but also the horizontal information of forest(crown density), which has the ability to extract forest three-dimensional structure parameters. Moreover, crown density variables could be employed to significantly improve the accuracy of estimating forest volumes using the random forest method. As variable selection reduces modeling accuracy to a certain extent, it was recommended that all the variables could be used for modeling estimation when high accuracy is required. However, for large data volumes, variable selection is still preferable. This approach of estimating forest volume using airborne LiDAR data is significantly superior to the approach utilizing optical data, thus, the proposed method serves as a basis for high-precision estimation and meets the requirements of rapid inversion of forest volume in large areas.

Key words: forest stock volume, small footprint LiDAR, random forest algorithm, variable selection

中图分类号: