基于背包激光雷达构建天然林蒙古栎胸径-树高模型.pdf
《基于背包激光雷达构建天然林蒙古栎胸径-树高模型.pdf》由会员分享,可在线阅读,更多相关《基于背包激光雷达构建天然林蒙古栎胸径-树高模型.pdf(8页珍藏版)》请在咨信网上搜索。
1、第 39 卷 第 5 期2023 年 9 月森 林 工 程FOREST ENGINEERINGVol.39 No.5Sep.,2023doi:10.3969/j.issn.1006-8023.2023.05.007基于背包激光雷达构建天然林蒙古栎胸径-树高模型杨军1,王帆2,宋仲禹1(1.黑龙江省林业和草原调查规划设计院,哈尔滨 150040;2.东北林业大学 林学院,哈尔滨 150040)摘 要:为探究基于背包激光雷达(Backpack Laser Scanning,BLS)点云数据提取天然林蒙古栎单木因子的精度,确定基于BLS 点云数据构建胸径-树高模型的可行性,以 3 块天然林蒙古栎样地
2、(30 m 30 m)的 BLS 点云数据为数据源,对点云数据进行预处理、归一化和单木分割等操作,并从分割后的单木中提取胸径和树高因子。通过对比提取数据与实测数据的相关性来分析单木因子的提取精度。最后利用提取出的单木因子构建天然林蒙古栎胸径-树高模型。结果表明,蒙古栎胸径提取结果的决定系数(R2)为 0.9300.957,均方根误差(RMSE)为 0.6970.897 cm;树高提取结果的 R2为 0.9250.951,RMSE 为1.4791.683 m。基于提取数据构建的天然林蒙古栎最优胸径-树高模型为模型 1(抛物线模型),模型检验指标 R2为 0.627,RMSE 为 2.044 m。
3、由此得出,BLS 有一定能力在天然林中进行数据采集,提取的蒙古栎胸径和树高与实测数据具有很高的相关性,结果可靠。基于 BLS 建立的蒙古栎胸径-树高模型精度较高,为 BLS 应用于天然林中的单木生长建模提供了进一步的理论基础。关键词:背包激光雷达;天然林;蒙古栎;单木因子;胸径-树高模型中图分类号:S758 文献标识码:A 文章编号:1006-8023(2023)05-0057-08Development of Natural Forest Quercus mongolica Diameter-Height Curve Model Based on BLSYANG Jun1,WANG Fan2
4、,SONG Zhongyu1(1.Heilongjiang Forestry and Grassland Investigation and Planning Institute,Harbin 150040,China;2.School of Forest,Northeast Forestry University,Harbin 150040,China)Abstract:To explore the accuracy of extracting individual tree factors of Quercus mongolica in natural forests based on B
5、ackpack Laser Scanning(BLS)point cloud data,and to determine the feasibility of constructing a diameter-height curve model based on BLS point cloud data,in this study,BLS point cloud data from three natural forest Quercus mongolica sample plots(30 m 30 m)were used as the data source.Then,the point c
6、loud data were preprocessed,normalized and individual tree segmented,and the DBH and tree height factors were extracted from the segmented individual trees.The extraction accuracy of individual tree factors was analyzed by comparing the correlation between the extracted data and the measured data.Fi
7、nally,natural forest Quercus mongolica diameter-height curve model was constructed by the extracted individual tree factors.The results showed that the coefficient of determination(R2)of Quercus mongolica DBH extraction results ranged 0.930-0.957,and the root mean square error(RMSE)ranged 0.697-0.89
8、7 cm,and R2 of tree height extraction results ranging 0.925-0.951 and RMSE ranging 1.479-1.683m.The optimal diameter-height curve model for Quercus mongolica in natural forests constructed based on extracted data was model I(parabolic model)with a model test index R2 of 0.627 and RMSE of 2.044 m.It
9、was concluded that BLS had the capability to conduct data collection in natural forests,and the ex-tracted Quercus mongolica DBH and tree height had a high correlation with the measured data,and the results were reliable.BLS-based Quercus mongolica diameter-height curve models were developed with hi
10、gh accuracy,which provided a further theoretical basis for the application of BLS to individual tree growth modeling in natural forests.Keywords:BLS;natural forests;Quercus mongolica;individual tree factor;diameter-height model收稿日期:2023-03-29基金项目:国家重点研发计划项目(2022YFD2201003-02);黑龙江省嘉荫县生态系统各碳库碳储量本底调查及动
11、态分析项目(HFW220100057)。第一作者简介:杨军,高级工程师。研究方向为森林经理学。E-mail:396285554 通信作者:宋仲禹,高级工程师。研究方向为森林经理学。E-mail:81830991 引文格式:杨军,王帆,宋仲禹.基于背包激光雷达构建天然林蒙古栎胸径-树高模型J.森林工程,2023,39(5):57-64.YANG J,WANG F,SONG Z Y.Development of natural forest Quercus mongolica diameter-height curve model based on BLSJ.Forest Engineering,
12、2023,39(5):57-64.森 林 工 程第 39 卷0 引言森林资源是人类赖以生存的重要自然资源之一,惯有“地球之肺”的美称,是全球气候变化的晴雨表,且和全球生物多样化存在着密切的关系,也是生物多样性的基础1。森林空间结构信息在森林生态研究、资源监测等方面具有重要作用,是陆地生态系统研究的重要参数,尤其在提高叶面积指数、森林生物量及进行森林生态系统碳循环、森林生长演替等动态过程的研究方面有重要意义2-3。而且,森林资源是全球资源的重要组成部分,具有缓解全球气候变化及促进社会经济发展的重要作用,对于构建良好生态环境有着独特意义4。为了有效地进行森林资源的管理,调查掌握森林资源信息是非常必
13、要的。传统的森林资源调查通过人工每木检尺来实现,劳动强度大、效率低,很难获得准确并且具有实际应用意义的森林结构参数5,有明显的局限性。在过去 20 年中,激光扫描技术作为一种林地信息收集、提取参数6与模型构建7的重要手段,在森林资源调查中得到了广泛的应用8,激光雷达扫描可以高效地获取大量的森林资源信息,这样节约很多的资源和成本,且对应的工作效率也明显地提高,有利于提高研究和应用水平。激光扫描技术可以在一定程度上取代人工作业,对森林进行高精度、快速和无损数据收集的特点,极大地提高了森林信息采集的效率和精度9。并且激光扫描技术作为一种高科技应用科学,对植被的空间结构与地形具有较强的检测能力,已成功
14、应用于森林参数的测量和反演,成为森林调查的重要手段10。激光雷达遥感作为主动遥感技术可以有效穿透森林,在获取森林高度和森林垂直结构参数方面有着其他光学遥感无法比拟的优势11-12,很多研究基于机载激光雷达(Airborne laser scanning,ALS)与地基激光雷达(Terrestrial laser scanning,TLS)成功提取了森林垂直结构及水平分布参数和单木结构参数13-17。例如,骆钰波等26基于 TLS数据提取了亚热带山区森林的胸径和树高,结果显示提取胸径的决定系数(R2)为 0.996,均方根误差(RMSE)为0.66 cm;提取树高的 R2为0.972,RMSE为
15、 0.92 m。黄旭等18基于 BLS 提取了人工落叶松的胸径和树高,研究显示胸径提取结果明显大于树高的提取结果,胸径平均提取精度为 94.51%,树高平均提取精度为 80.27%。蒙古栎(Quercus mongolica),主要分布于东北地区,是我国的主要用材树种和东北次生阔叶林中的主要组成树种,有着极高的经济价值与生态效益,但是长期以来一直未受到重视。人们对蒙古栎的生长特性尚未全面的研究,且大部分蒙古栎在现实林分中多生长不良,长期以来一直被作为改造的对象,导致蒙古栎未能充分发挥其多种效益。因此,为了准确了解蒙古栎的生长状况,促进蒙古栎林的合理经营,必须构建完善的蒙古栎生长模型。目前机载激
16、光雷达和地基激光雷达在林业方面的应用和研究已经趋于成熟,但是这 2 种激光雷达在进行外业活动采集数据的过程中成本十分高昂,且数据处理较为繁琐,因此在小范围样地调查时并不是十分适用。背包激光雷达作为一种新型便携的激光扫描仪十分适用于样地尺度下的调查,但在林业方面的应用还比较少,因此,准确评估利用背包激光雷达提取单木因子的精度是亟待解决的问题。综上所述,本研究使用背包激光雷达点云数据提取蒙古栎的单木因子,在提取精度满足生产需求的情况下进一步利用提取的因子构建蒙古栎胸径-树高模型。以期为背包激光雷达在样地尺度下的应用提供进一步的实践基础,同时为蒙古栎生长模型系统的构建做出一定贡献。1 研究区概况研究
17、区位于黑龙江省佳木斯市孟家岗林场(1303213052 E,46204630 N),地处于完达山西麓余脉,地貌主要为低山丘陵;坡度平缓,介于 1020,平均海拔约 250 m。孟家岗林场属于东亚大陆性季风气候,春季干旱少雨,秋季易霜冻,夏季短暂且温暖湿润,冬季持久且寒冷干燥。最高气温可达 35.6,最低气温为-34.7,年均气温约为 2.7。无霜期大约 120 d,年平均降水量550670 mm。土壤类型主要为典型暗棕壤,还有少量草甸暗棕壤、潜育暗棕壤和原始暗棕壤;除此之外,还存在草甸土、白浆土、泥炭土以及沼泽土。该林场以人工林为主,占整个施业区面积的 2/3,天然次生林约占 1/3,森林覆盖
18、率为 81.7%。2 数据收集2.1 野外样地数据调查2021 年 4 月,在孟家岗林场的蒙古栎天然林中设置了 3 块固定标准样地,样地面积均为 0.09 hm285第 5 期杨军,等:基于背包激光雷达构建天然林蒙古栎胸径-树高模型(30 m 30 m)。在设置的标准样地的中心点用GPS 进行定位,定位精度为 3 m。每木检尺工作包括测量每株树木的树高、胸径、冠幅和第一活枝高,并且在样地中记录每棵树的相对位置。采用超声波测高器来测量树高,用胸径尺测量样地内所有胸径大于 5 cm 的树木,用测距仪来测量冠幅,用 50 m皮尺确定样地内树木的相对位置,共计测量 244 株活立木。样地信息见表 1。
19、表 1 样地信息统计表Tab.1 Statistical table of sample plot information样地号Plot 树种Species株数Number of trees胸径/cmDiameter at breast height(DBH)最大值Maximum最小值Minimum平均值Mean标准差SD树高/mTree height最大值Maximum最小值Minimum平均值Mean标准差SD1蒙古栎8932.65.117.744.3921.06.015.212.612蒙古栎8528.36.418.614.9622.06.016.012.923蒙古栎7032.78.021
20、.205.2222.48.217.423.192.2 背包激光雷达数据采集每木检尺工作结束后开始利用背包激光雷达(北京数字绿土有限公司)采集样地的点云数据。由于样地面积较大,为了保证能较为全面地获取样地内的点云数据,采用了黄旭等18的路线规划,如图 1 所示。最终采集到的点云数据,如图 2 所示。图 1 样地扫描路径规划图Fig.1 Map of scanning route planning图 2 样地点云示意图Fig.2 Schematic diagram of a sample plot point cloud3 研究方法3.1 点云数据预处理在 Lidar360 软件中加载扫描获得的.
21、ply 格式的样地点云数据以及.xyz 格式的样地采集路线行走轨迹,根据行走轨迹获取样地在点云数据中的大致位置,从而对点云数据进行裁剪,精确获取样地内的点云数据。采用高度阈值法19对样地点云进行去噪处理,随后根据渐进加密三角网滤波算法20对地面点进行分类。利用分离出的地面点生成数字高程模型从而获得归一化的点云数据。3.2 提取单木因子利用归一化的点云数据在 Lidar360 软件中进行单木自动分割。但是这种自动分割方法得到的单木点云数据易受低矮灌木干扰,并且连株树木易被分割为同一单木,对后续数据处理造成困难。因此利用软件的单木点云编辑功能对分割后的点云数据进行二次处理。在 2D 视窗中提取一定
22、范围的树木剖面,对样地数据逐片排查进行编辑。多株树木分割为一株的用创建单木功能进行再分割;一株树木被分为多株而过分割的,用合并单木功能重新划分为一株单木。对单木与周围灌木因为距离过近而分割为一株单木,影响胸径提取精度的,用创建单木功能将多余的灌木部分分割出来并使用删除单木功能将其删除;对低矮灌木成片生长被分割为多株单木的,用删除单木功能进行修正。处理后单木分割结果如图 3 所示。最后从修正后的单木点云数据中截取胸径处上下 5 cm 厚的点云进行胸径的提取。单木点云中 Z 轴上的最大值与最小值之95森 林 工 程第 39 卷差定义为树高。图 3 单木分割结果Fig.3 Segmentation
23、results for individual trees3.3 构建胸径-树高模型3.3.1 基础模型根据马武等21对天然蒙古栎单木生长模型的研究结果,从中选取了 5 种胸径-树高模型作为备选模型。模型的具体形式如下。H1=a+bD+cD2。(1)H2=1.3+a(1-e-bD)c。(2)H3=1.3+a(1-e-bDc)。(3)H4=1.3+aDbDc。(4)H5=1.3+a1+1bDc。(5)式中:H 为树高,m;D 为胸径,cm;a、b、c 分别为模型的参数。3.3.2 模型评价与检验对本研究所用的数据进行随机抽样,抽出 75%的样本用作建模数据,剩余 25%样本用来进行模型检验。利用
24、R2和 RMSE 评价模型的拟合能力。R2越大,RMSE 越小,则证明模型的拟合效果越好。利用独立样本对模型进行检验,选择平均相对偏差绝对值(RMAE)、平均绝对偏差(MAE)、预估精度(forecast precision,FP,式中为 Fp)作为评价指标。Fp越大,MAE、RMAE 越小,模型的预测效果会越好22。模型的拟合与检验都是在 R 软件中实现的。预估精度(Fp)计算方式如下。Fp=1-t0.05y-(yi-yi)2n n-p()()100%。(6)式中:yi为观测值;yi为预测值;n 为样本个数;p 为模型的参数个数;y-i为所有观测值的平均值;t0.05为95%下显著性水平,本
25、研究取值为 2.306。4 结果与分析4.1 单木因子提取精度表 2 展示了胸径与树高提取结果的统计信息。为方便后续进行建模,将提取结果按照 3 1的比例随机划分为建模数据和检验数据。图 4 展示的是3 块样地中蒙古栎点云数据提取出的胸径和树高与实测数据之间的回归关系。从总体上看,所有的散点都集中在 1 1 参考线的附近,证明提取出的胸径和树高与实测数据之间有着极强的相关性,R2最低为 0.925,最高为 0.957。胸径的提取精度整体上优于树高的提取精度,这是由于背包激光雷达获取冠下点云信息的能力大于获取冠上点云信息的能力造成的。样地 1 的胸径和树高的提取结果都最差,这可能是由于样地 1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 背包 激光雷达 构建 天然林 蒙古 胸径 模型
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。