基于地基激光雷达的单木森林参数提取分析

摘要 [目的]利用激光雷达手段大范围估算森林参数。[方法]以安徽省黄山市休宁县岭南林场森林资源为研究对象，通过地面调查获取3种森林类型单木胸径和树高数据，并同步使用地基激光雷达获取点云数据，利用Hough变换检测算法提取每木胸径和树高，在Matlab2020a软件下构建基于地面调查数据和激光雷达提取值的估算模型。[结果]针叶林、阔叶林和毛竹林的胸径估测精度（R2c）分别为0.921、0.921和0.996，RMSEc分别为2.080、2.460和0.291，模型表现显著;树高估测精度分别为0.908、0.846和0.500，RMSEc分别为1.120、1.620和1.210。[结论]该研究可为运用激光雷达手段大范围估算森林参数提供依据，也可为森林类型点云识别与分割提供参考。

关键词 岭南林场;地基激光雷达;森林参数;回归分析;精度评价

中图分类号 S 758文献标识码 A文章编号 0517-6611（2022）04-0128-04

doi：10.3969/j.issn.0517-6611.2022.04.033

开放科学（资源服务）标识码（OSID）：

Extraction and Analysis of Forest Parameters for Individual Tree Based on Terrestrial Laser Scanning

FU Gen-shen1， CHEN Jian2， HUANG Qing-feng1 et al

（1. School of Forestry and Landscape Architecture， Anhui Agricultural University， Hefei， Anhui 230061; 2. KQ GEO Technologies Co.， Ltd， Beijing 100176）

Abstract [Objective]To estimate forest parameters in a wide range by lidar. [Method]The study studied by forest resources in Lingnan forestry site， Xiuning County， Anhui Province， single wood diameter at breast height and tree height data of three forest types were obtained by ground survey， and their point cloud data were obtained simultaneously by terrestrial laser scanning （TLS）. Each wood diameter at breast height and tree height were extracted by Hough transform detection algorithm， and an estimation model based on ground survey data and TLS extracted values was constructed under Matlab2020a software. [Result]The results showed that the estimation accuracy （R2c） of breast diameter was 0.921， 0.921 and 0.996 for coniferous， broadleaf and moso bamboo forests， respectively， with RMSEc of 2.080， 2.460 and 0.291， and the model performance was significant; the estimation accuracy of tree height was 0.908， 0.846 and 0.500， respectively， with RMSEc of 1.120， 1.620 and 1.210， respectively. [Conclusion]This study can provide a basis for estimating forest parameters in a large scale by TLS， and also provide a reference for forest type point cloud identification and segmentation.

Key words Lingnan forest farm;Terrestrial laser scanning;Forest parameters;Regression analysis;Precision evaluation

基金项目 安徽省自然科学基金项目（1808085QC74）。

作者简介 傅根深（1997—），男，安徽合肥人，硕士研究生，研究方向：森林经理。*通信作者，讲师，博士，从事森林资源经营管理和林业3S技术研究。

收稿日期 2021-06-15

树木胸径和树高是森林调查中常用的参数，传统方法是通过胸径卷尺和测高器进行测量获取，但是这种方法存在着工作量大、人为误差等不足。除此之外，对于地形复杂、坡度大、调查人员难以进入的区域，这些林木因子获取难度较大。遥感技术一定程度上能够弥补这些不足，光学遥感技术在二维平面能够获得较好的水平森林参数，但是在获取树高、胸径、密度等三维参数时显得不足。作为一种新兴的主动遥感技术，激光雷达发射的激光脉冲则可轻易地从树木顶端到达地表，甚至穿透地表到达土壤内部，这种获取森林垂直结构参数的技术目前大部分光学遥感传感器难以做到[1]。地基激光雷达（terrestrial laser scanning，TLS）是一种基于激光的新型遥感仪器，通过光束偏转机制获取地表物体的三维点云数据，可以同时量测较高精度的地物范围和地物之间的角度，TLS获取的点云数据具有高密度、高精度等特点，可轻易捕捉林分尤其是下层林分的结构信息，这对于精确掌握森林生态系统的变化并研究这些变化对整个陆地生态系统的影响起着举足轻重的作用[2]。但是，由于大部分森林冠层之间相互遮挡的影响，机载激光雷达（airborne LiDAR scanning，ALS）则难以准确分离乔木层下的灌木和草本，ALS的优势体现在大区域森林垂直结构的快速获取，TLS则要耗费庞大的人力物力才能做到[3]。

TLS技术在林业中的应用可以追溯到2002年，其以快速、无损、自动的方式获取胸径、树高、树冠宽度、郁闭度等森林参数[4-7]。研究表明，从TLS点云中提取胸径的方法很多，如二维曲线拟合、圆检测、圆柱拟合等，并都得到了较高的拟合精度[8-9]。Ma等[10]使用协方差特征分类得到叶片点云并实现了LAI的反演。赵方博等[11]通过对数据预处理提取出树冠点云，将其模拟为半球图像后运用球极平面投影和Lambert方位角等面积投影，通过统计面积的方法分别计算不同投影方法和图像划分方法下的孔隙率，进而计算出真实叶面积指数。骆钰波等[12]以主干点云为基础，采用迭代最小二乘拟合圆柱的方法自动提取树木位置、胸径;构建点云的八叉树结构，利用体素的空间邻接性实现点云分割，自动提取树高，估测R2超过0.9。马振宇等[13]基于自适应临近搜索法和递归特征排除法（RFE）筛选重要变量进行随机森林分类，有效提取倒木数据。张天安等[14]通过对原始点云进行枝叶分离并根据Dijkstra距离法和加权匹配法获取活立木单株点云，分别构建了含笑和樱花的枝干三维模型。杨水荣等[15]采用切片形式的Alpha Shape 3D法估算树木冠层体积，并定量反映树木冠层分布及体积变化情况。

笔者以安徽省黄山市休宁县岭南林场为研究区域，通过样地设置及样木调查，获取针叶林、阔叶林、毛竹林3种森林类型的单木胸径及树高数据。综合利用GPS、全站仪、TLS获取样木的激光点云数据，通过坐标匹配实现单木与激光数据的吻合，在此基础上提取样木的胸径和树高，并进一步构建3种森林类型单木的胸径和树高估测模型，利用留一法验证模型的精度。通过对森林参数的精确估测，可为下一步应用机载激光雷达实现大区域森林参数估算提供科学依据，对开展森林蓄积量估算、森林生物量定量遥感模型构建、碳达峰和碳中和等生态监测具有重要意义。

1 研究区概况

研究数据均采集于安徽省黄山市休宁县岭南林场，林场位于安徽省最南端，与江西省交界（图1）。岭南林场地处中纬度地带，属于北亚热带湿润季风气候，年均气温16.2 ℃，日照时数1 800 h左右，无霜期220 d。林场地形表现为丘陵—低山—中山3种类型，其中丘陵面积最大。土壤类型主要为黄壤、黄棕壤。林场森林覆盖率较高，植被分布具有典型的过渡地带特征，森林类型以针叶林、针阔混交林为主。常见的主要树种包括杉木[Cumninghamia lanceolata（Lamb.）Hook.]、马尾松（Pinus massoniana Lamb.）、枫香（Liquidambar formosana Hance）、毛竹[Phyllostachys heterocycla（Carr.）Mitford cv.Pubescens]、木荷（Schima superba Gardn.et Champ.）、栓皮栎（Quercus variabilis Bl.）、苦槠（Castanopsis sclerophylla）等。

2 研究方法

2.1 样木调查

该试验涉及3种森林类型，分别为针叶林、阔叶林和毛竹林。样地设置为20 m×20 m的正方形。为了精细提取单木森林参数，每木检尺以2 cm为起测胸径，采用传统的胸径卷尺进行测量，使用激光测高器获取单木树高。利用全站仪、罗盘仪、皮尺获取样地4个角点位置，通过静态GPS获取其坐标。

2.2 激光点云数据采集

利用Trimble TX8设备进行地面样木点云数据采集。Trimble TX8具有轻便、携带方便的特点，最大扫描距离达到340 m，激光点云高达555百万点，垂直扫描角度为317°，激光波长达到1 500 nm。每个样地在样地的4个角点和中心位置共设置5个扫描站点。每站扫描时间为10 min。

2.3 数据坐标匹配

TLS、GPS、全站仪数据采集完毕后，利用TLS仪器自带Trimble RealWorks系统软件对采集到的原始TLS数据进行反射靶标识别、多站拼接、坐标转换和生成点云等预处理操作，利用后处理软件对GPS数据进行差分处理，对全站仪坐标数据进行点校正处理，以实现单木坐标的匹配[16]。

2.4 激光点云数据处理

由于受到目标物体的表面形状、方向和纹理的影响，激光脉冲被反射后形状会发生改变，尤其是扫描到目标物边缘或发射脉冲通过植被时都将较大地改变反射脉冲形状，从而产生较大偏差。原始的点云数据会接收到多个地物的回波信息，并产生密度非常高的点云数据，在进行森林参数（如胸径和树高）提取时，只需要保留树干最外部点云即可。利用TerraSolid软件的TerraScan模块对初始点云数据进行预处理，主要包括数据过滤、非树干点去除、数据抽稀[16]。

2.5 森林参数提取

利用Hough变换检测算法来估测样地内每个单株立木的胸径，然后根据树木的形态学特征，通过设置一系列阈值去除非树干处的拟合圆，可提高单木的位置与胸径识别精度[1]。最后，通过地基激光雷达扫描数据获取的单木树干生长方向与该方向上的单木垂直结构提取单木高[16]。

2.6 森林参数估测模型构建及精度验证

在Matlab2020a环境下利用回归分析方法分别构建基于实际测量值和TLS提取值的单木胸径、树高估测模型。最后采用留一法交叉验证对构建的胸径、树高估测模型进行精度验证。留一法交叉验证是指每个森林类型中采集到n株样本，每次选取1个样本作为测试集，其余样本（n-1）作为训练集。循环迭代需要运算n次训练，最后取n次训练的平均值进行检验。利用决定系数R2、均方根误差RMSE、P值等作为建模过程的评价指标。

3 结果与分析

3.1 样木数据处理结果

针叶林样地主要树种为杉木，根据现场调查及TLS提取数据，共得到有效样木70株。阔叶林样地主要树种为木荷、短柄枹、檵木、甜槠等，优势树种为木荷，有效样木119株。毛竹林有效样木320株。