留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于三维激光点云的城市绿化树种材积及树干碳储量无损精确测算

徐诗宇 施拥军 冯晟斐

徐诗宇, 施拥军, 冯晟斐. 基于三维激光点云的城市绿化树种材积及树干碳储量无损精确测算[J]. 浙江农林大学学报, 2018, 35(6): 1062-1069. doi: 10.11833/j.issn.2095-0756.2018.06.009
引用本文: 徐诗宇, 施拥军, 冯晟斐. 基于三维激光点云的城市绿化树种材积及树干碳储量无损精确测算[J]. 浙江农林大学学报, 2018, 35(6): 1062-1069. doi: 10.11833/j.issn.2095-0756.2018.06.009
XU Shiyu, SHI Yongjun, FENG Shengfei. Nondestructive and accurate measurement of volume and stem carbon storage for urban greening tree species based on terrestrial laser scanning point cloud[J]. Journal of Zhejiang A&F University, 2018, 35(6): 1062-1069. doi: 10.11833/j.issn.2095-0756.2018.06.009
Citation: XU Shiyu, SHI Yongjun, FENG Shengfei. Nondestructive and accurate measurement of volume and stem carbon storage for urban greening tree species based on terrestrial laser scanning point cloud[J]. Journal of Zhejiang A&F University, 2018, 35(6): 1062-1069. doi: 10.11833/j.issn.2095-0756.2018.06.009

基于三维激光点云的城市绿化树种材积及树干碳储量无损精确测算

doi: 10.11833/j.issn.2095-0756.2018.06.009
基金项目: 

国家自然科学基金资助项目 31370637

国家自然科学基金资助项目 61190114

国家自然科学基金资助项目 31500520

中央财政林业推广项目 2017TS07

中国绿色碳汇基金项目 H20170049

浙江省重点科学技术项目 2015C03008

详细信息
    作者简介: 徐诗宇, 从事森林碳汇计量与监测研究。E-mail: 994653319@qq.com
    通信作者: 施拥军, 副教授, 从事林业3S技术应用、森林碳汇等研究。E-mail:syjwwh@163.com
  • 中图分类号: S758;S771

Nondestructive and accurate measurement of volume and stem carbon storage for urban greening tree species based on terrestrial laser scanning point cloud

  • 摘要: 城市森林是森林碳库的重要组成部分。为探索三维激光扫描技术(TLS)在无损精确测算城市绿化树种树干材积及碳储量方面的潜力,了解不同区分段长度对基于TLS测算树干材积的影响并对现有模型测算的树干材积及碳储量的精度进行了评估,采集30株无患子Sapindus mukorossi点云,并按径阶采集树干钻芯样本,研究在不同区分段长度下测算的树干材积及碳储量,并与现有模型的测算结果相比较。结果表明:①基于TLS能够提取高精度的立木胸径和树高(平均相对误差分别为1.97%和4.85%),可实现城市绿化树种树干材积及碳储量的无损精确测算。②基于TLS在区分段长度小于0.50 m时能测算得到高精度的树干材积。③城市绿化树种的材积和树干碳储量与野外林分存在显著差异(P < 0.05)。现有模型并不适用于城市绿化树种。
  • 图  1  无患子样木的径阶分布

    Figure  1  Diameter grade distribution of Sapindus mukorossi samples

    图  2  TLS提取的胸径、树高与实测值比较

    Figure  2  Comparison of diameter and height extracted base on TLS with actual measurement

    图  3  不同区分段长度下材积测算误差的分布及比较

    Figure  3  Distribution and comparison of volume measurement error in different section lengths

    图  4  现有模型测算的城市绿化树种材积及树干碳储量与TLS测算结果的比较

    Figure  4  Comparison of the urban greening tree species volume and stem carbon storage measured based on existing model

    表  1  不同区分段长度下测算的材积比较

    Table  1.   Comparison of volume measured based on different section lengths

    区分段长度/m TLS测算材积/dm3 标准误差(ERMSE)/dm3 乖离率(sbias)/dm3
    0.10 37.93 ± 10.96 a
    0.20 37.92 ± 10.84 a 0.44 0.02
    0.40 38.57 ± 10.86 a 1.08 0.49
    0.50 38.44 ± 11.29 a 1.15 0.44
    0.60 39.15 ± 11.02 a 1.69 0.92
    0.80 39.50 ± 11.24 a 2.08 1.38
    1.00 40.10 ± 11.57 a 2.79 2.17
    说明:*每个数据表示样本平均值±标准差;不同小写字母表示不同区分段测算的材积之间差异显著(P<0.05)
    下载: 导出CSV

    表  2  不同方法测算的材积、树干碳储量比较

    Table  2.   Comparison of volume and stem carbon storage measured by different methods

    测算方法 测算材积/dm3 测算碳储量/kg
    TLS_0.10 m区分求积 37.93 ± 10.97 c 12.38 ± 5.64 d
    TLS_1.00 m区分求积 40.10 ± 11.57 c 13.09 ± 5.77 cd
    一元材积模型 68.77 ± 36.58 ab 22.26 ± 11.80 a
    二元材积模型 54.54 ± 31.91 b 17.64 ± 10.25 bc
    碳储量模型 17.92 ± 9.59 bc
    说明:样本平均值±标准差;不同小写字母表示不同区分段测算的材积之间差异显著(P<0.05)
    下载: 导出CSV
  • [1] 罗云建.中国森林生态系统生物量及其分配研究[M].北京:中国林业出版社, 2013.
    [2] GRIMM N B, FAETH S H, GOLUBIEWSKI N E, et al. Global change and the ecology of cities[J]. Science, 2008, 319(5864):756-760.
    [3] KAYE J P, McCULLEY R I, BURKE I C. Carbon fluxes, nitrogen cycling, and soil microbial communities in adjacent urban, native and agricultural ecosystems[J]. Global Change Biol, 2005, 11(4):575-587.
    [4] McPHERSON E G. Structure and sustainability of Sacramento's urban forest[J]. J Arboric, 1998, 24(4):174-189.
    [5] SIMPSON J R. Urban forest impacts on regional cooling and heating energy use:Sacramento County case study[J]. J Arboric, 1998, 24(4):201-209.
    [6] CROWL T A, CRIST T O, PARMENTER R R, et al. The spread of invasive species and infectious disease as drivers of ecosystem change[J]. Front Ecol Environ, 2008, 6(5):238-246.
    [7] WATT P J. Measuring forest structure with terrestrial laser scanning[J]. Int J Remote Sens, 2005, 26(7):1437-1446.
    [8] 罗旭, 冯仲科, 邓向瑞, 等.三维激光扫描成像系统在森林计测中的应用[J].北京林业大学学报, 2007, 29(增刊2):82-87.

    LUO Xu, FENG Zhongke, DENG Xiangrui, et al. Application of 3D laser scanning imaging system in forest measuring[J]. J Beijing For Univ, 2007, 29(suppl 2):82-87.
    [9] 邓向瑞, 冯仲科, 罗旭.三维激光扫描系统在林业中的应用研究[J].北京林业大学学报, 2005, 27(增刊2):43-47.

    DENG Xiangrui, FENG Zhongke, LUO Xu. Application of three-dimensional laser scanning system in forestry[J]. J Beijing For Univ, 2005, 27(suppl 2):43-47
    [10] 赵阳, 余新晓, 信忠保, 等.地面三维激光扫描技术在林业中的应用与展望[J].世界林业研究, 2010, 23(4):41-45.

    ZHAO Yang, YU Xinxiao, XIN Zhongbao, et al. Application and outlook of terrestrial 3D laser scanning technology in forestry[J]. World For Res, 2010, 23(4):41-45
    [11] LIANG Xinlian, HYYPPÄJ, KANKARE V, et al. Stem curve measurement using terrestrial laser scanning[J]. IEEE Trans Geosci Remote Sen, 2013, 52(3):1739-1748.
    [12] 梁子瑜.基于TLS点云数据的野外林分调查因子测定及收获估计[D].南京: 南京林业大学, 2015.

    LIANG Ziyu. The Extraction of Forest Stand Parameters and the Estimation with Yield based on Point Cloud Data from TLS[D]. Nanjing: Nanjing Forestry University, 2015.
    [13] 尤磊.基于点云数据的树干干形测量[D].北京: 中国林业科学研究院, 2016.

    YOU Lei. Stem form Measurement based on Point Cloud Data[D]. Beijing: Chinese Academy of Forestry, 2016.
    [14] 王剑, 周国民.基于激光扫描仪的树干三维重建方法研究[J].微计算机信息, 2009, 25(24):228-230.

    WANG Jian, ZHOU Guomin. 3D reconstruction of tree stems based on laser scanning of standing[J]. Control Autom, 2009, 25(24):228-230.
    [15] 石银涛, 程效军, 贾东峰.三维激光扫描树木模型在林业中的应用[J].测绘通报, 2012(3):40-42, 52.

    SHI Yintao, CHEN Xiaojun, JIA Dongfeng, et al. Applications of 3D laser scanning based tree model toforestry[J]. Bull Surv Map, 2012(3):40-42, 52.
    [16] 韦雪花, 王永国, 郑君, 等.基于三维激光扫描点云的树冠体积计算方法[J].农业机械学报, 2013, 44(7):235-240.

    WEI Xuehua, WANG Yongguo, ZHENG Jun, et al. Tree crown volume calculation based on 3-D laser scanning point clouds data[J]. Trans Chin Soc Agric Mach, 2013, 44(7):235-240.
    [17] 史琰.中国城市建成区植被结构特征和碳吸收[D].杭州: 浙江大学, 2013.

    SHI Yan. Vegetation Structure Characteristics and Carbon Uptake of Urban Built-up Area in China[D]. Hangzhou: Zhejiang University, 2013.
    [18] 胡海清, 罗碧珍, 魏书精, 等.小兴安岭7种典型林型林分生物量碳密度与固碳能力[J].植物生态学报, 2015, 39(2):140-158.

    HU Haiqing, LUO Bizhen, WEI Shujing, et al. Biomass carbon density and carbon sequestration capacity in seven typical forest types of the Xiaoxing'an Mountains, China[J]. Chin J Plant Ecol, 2015, 39(2):140-158.
    [19] TANHUANPÄÄT, KANKARE V, SETÄLÄH, et al. Assessing above-ground biomass of open-grown urban trees:a comparison between existing models and a volume-based approach[J]. Urban For Urban Green, 2017, 21:239-246.
    [20] RAHMAN M A, BAKAR M A, RAZAK K, et al. Non-destructive, laser-based individual tree aboveground biomass estimation in a tropical rainforest[J]. Forests, 2017, 8(3):86. doi:10.3390/f8030086.
    [21] BROLLY G, KIRÁLY G. Algorithms for stem mapping by means of terrestrial laser scanning[J]. Acta Silv Lign Hung, 2009, 5(2):119-130.
    [22] TANSEY K, SELMES N, ANSTEE A, et al. Estimating tree and stand variables in a Corsican pine woodland from terrestrial laser scanner data[J]. Int J Remote Sens, 2009, 30(19):5195-5209.
    [23] WATT P J, DONOGHUE D N M. Measuring forest structure with terrestrial laser scanning[J]. Int J Remote Sens, 2005, 26(7):1437-1446.
    [24] 刘安兴, 龚赛凤, 孙源和.浙江省林木根径材积表及其数学模型[J].华东森林经理, 1992, 6(3):12-18.

    LIU Anxing, GONG Saifeng, SUN Yuanhe. Study on tree root volume table and its mathematical model in Zhejiang Province[J]. East China For Manage, 1992, 6(3):12-18.
    [25] 孟宪宇.测树学[M].北京:中国林业出版社, 2006.
    [26] SAARINEN N, KANKARE V, VASTARANTA M, et al. Feasibility of terrestrial laser scanning for collecting stem volume information from single trees[J]. ISPRS J Photogramm Remote Sens, 2017, 123:140-158. doi:10.1016/j.isprsjprs.2016.11.012.
    [27] 嵇浩翔, 史琰, 朱轶梅, 等.杭州市不同土地利用类型的树木生长和碳固存[J].生态学杂志, 2011, 30(11):2405-2412.

    JI Haoxiang, SHI Yan, ZHU Yimei, et al. Tree growth and carbon sequestration in different land-use types in Hangzhou City[J]. Chin J Ecol, 2011, 30(11):2405-2412.
  • [1] 赵凯, 唐丽华, 张姝婧.  基于OpenGL的交互式三维树木建模与可视化研究 . 浙江农林大学学报, 2019, 36(1): 138-147. doi: 10.11833/j.issn.2095-0756.2019.01.017
    [2] 周静平, 李存军, 胡海棠, 陶欢, 彭代亮, 谢春春, 葛艳.  森林植被遥感监测影像最佳分辨率选择 . 浙江农林大学学报, 2018, 35(4): 716-723. doi: 10.11833/j.issn.2095-0756.2018.04.017
    [3] 陶江玥, 刘丽娟, 庞勇, 李登秋, 冯云云, 王雪, 丁友丽, 彭琼, 肖文惠.  基于机载激光雷达和高光谱数据的树种识别方法 . 浙江农林大学学报, 2018, 35(2): 314-323. doi: 10.11833/j.issn.2095-0756.2018.02.016
    [4] 董心玉, 范文义, 田甜.  基于面向对象的资源3号遥感影像森林分类研究 . 浙江农林大学学报, 2016, 33(5): 816-825. doi: 10.11833/j.issn.2095-0756.2016.05.013
    [5] 陈金星, 岳德鹏, 冯仲科, 丁家巍, 姚炳全, 叶添雄.  手持式树径自动识别测树仪的研制与应用 . 浙江农林大学学报, 2016, 33(4): 589-598. doi: 10.11833/j.issn.2095-0756.2016.04.006
    [6] 吕常笑, 邓华锋, 王少杰, 陈振雄, 王雪军.  马尾松不同区域相容性立木材积和地上生物量模型 . 浙江农林大学学报, 2016, 33(5): 790-797. doi: 10.11833/j.issn.2095-0756.2016.05.010
    [7] 姚智, 张晓丽.  基于WebGIS平台的森林郁闭度遥感反演信息系统研建 . 浙江农林大学学报, 2015, 32(3): 392-398. doi: 10.11833/j.issn.2095-0756.2015.03.009
    [8] 王月婷, 张晓丽, 杨慧乔, 王书涵, 白金婷.  基于Landsat 8卫星光谱与纹理信息的森林蓄积量估算 . 浙江农林大学学报, 2015, 32(3): 384-391. doi: 10.11833/j.issn.2095-0756.2015.03.008
    [9] 陈婷, 施拥军, 周国模, 郑泽睿, 李翀.  毛竹碳汇林营造初期林分非空间结构年际变化特征 . 浙江农林大学学报, 2015, 32(2): 181-187. doi: 10.11833/j.issn.2095-0756.2015.02.003
    [10] 李梦, 施拥军, 周国模, 周大勇, 张娇, 张英海.  浙江省嘉兴市高速公路造林碳汇计量 . 浙江农林大学学报, 2014, 31(3): 329-335. doi: 10.11833/j.issn.2095-0756.2014.03.001
    [11] 高士增, 张怀清, 刘闽, 何清平, 罗立平.  树木枝干点云数据的等值线提取 . 浙江农林大学学报, 2013, 30(5): 648-654. doi: 10.11833/j.issn.2095-0756.2013.05.003
    [12] 王晓宁, 徐天蜀, 李毅.  利用ALOS PALSAR双极化数据估测山区森林蓄积量模型 . 浙江农林大学学报, 2012, 29(5): 667-670. doi: 10.11833/j.issn.2095-0756.2012.05.005
    [13] 吴兆艳, 汤孟平.  基于SketchUp与GIS的森林景观可视化实现 . 浙江农林大学学报, 2012, 29(3): 352-358.
    [14] 季碧勇, 陶吉兴, 张国江, 杜群, 姚鸿文, 徐军.  高精度保证下的浙江省森林植被生物量评估 . 浙江农林大学学报, 2012, 29(3): 328-334. doi: 10.11833/j.issn.2095-0756.2012.03.002
    [15] 张国江, 季碧勇, 王文武, 方炎杰, 李佐晖.  设区市森林资源市县联动监测体系研究 . 浙江农林大学学报, 2011, 28(1): 46-51. doi: 10.11833/j.issn.2095-0756.2011.01.008
    [16] 王祖华, 蔡良良, 关庆伟, 蔡霞.  淳安县森林生态系统服务价值评估 . 浙江农林大学学报, 2010, 27(5): 757-761. doi: 10.11833/j.issn.2095-0756.2010.05.019
    [17] 徐文兵, 汤孟平.  全站仪双边交会法测定树木三维坐标 . 浙江农林大学学报, 2010, 27(6): 815-820. doi: 10.11833/j.issn.2095-0756.2010.06.002
    [18] 徐文兵, 高飞, 杜华强.  几种测量方法在森林资源调查中的应用与精度分析 . 浙江农林大学学报, 2009, 26(1): 132-136.
    [19] 李亦秋, 仲科.  山东省森林资源动态变化的非等间距灰色预测 . 浙江农林大学学报, 2009, 26(1): 7-12.
    [20] 夏伟伟, 韩海荣, 伊力塔, 程小琴.  庞泉沟国家级自然保护区森林景观格局动态 . 浙江农林大学学报, 2008, 25(6): 723-727.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2018.06.009

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2018/6/1062

图(4) / 表(2)
计量
  • 文章访问数:  2976
  • HTML全文浏览量:  594
  • PDF下载量:  567
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-12-19
  • 修回日期:  2018-03-22
  • 刊出日期:  2018-12-20

基于三维激光点云的城市绿化树种材积及树干碳储量无损精确测算

doi: 10.11833/j.issn.2095-0756.2018.06.009
    基金项目:

    国家自然科学基金资助项目 31370637

    国家自然科学基金资助项目 61190114

    国家自然科学基金资助项目 31500520

    中央财政林业推广项目 2017TS07

    中国绿色碳汇基金项目 H20170049

    浙江省重点科学技术项目 2015C03008

    作者简介:

    徐诗宇, 从事森林碳汇计量与监测研究。E-mail: 994653319@qq.com

    通信作者: 施拥军, 副教授, 从事林业3S技术应用、森林碳汇等研究。E-mail:syjwwh@163.com
  • 中图分类号: S758;S771

摘要: 城市森林是森林碳库的重要组成部分。为探索三维激光扫描技术(TLS)在无损精确测算城市绿化树种树干材积及碳储量方面的潜力,了解不同区分段长度对基于TLS测算树干材积的影响并对现有模型测算的树干材积及碳储量的精度进行了评估,采集30株无患子Sapindus mukorossi点云,并按径阶采集树干钻芯样本,研究在不同区分段长度下测算的树干材积及碳储量,并与现有模型的测算结果相比较。结果表明:①基于TLS能够提取高精度的立木胸径和树高(平均相对误差分别为1.97%和4.85%),可实现城市绿化树种树干材积及碳储量的无损精确测算。②基于TLS在区分段长度小于0.50 m时能测算得到高精度的树干材积。③城市绿化树种的材积和树干碳储量与野外林分存在显著差异(P < 0.05)。现有模型并不适用于城市绿化树种。

English Abstract

徐诗宇, 施拥军, 冯晟斐. 基于三维激光点云的城市绿化树种材积及树干碳储量无损精确测算[J]. 浙江农林大学学报, 2018, 35(6): 1062-1069. doi: 10.11833/j.issn.2095-0756.2018.06.009
引用本文: 徐诗宇, 施拥军, 冯晟斐. 基于三维激光点云的城市绿化树种材积及树干碳储量无损精确测算[J]. 浙江农林大学学报, 2018, 35(6): 1062-1069. doi: 10.11833/j.issn.2095-0756.2018.06.009
XU Shiyu, SHI Yongjun, FENG Shengfei. Nondestructive and accurate measurement of volume and stem carbon storage for urban greening tree species based on terrestrial laser scanning point cloud[J]. Journal of Zhejiang A&F University, 2018, 35(6): 1062-1069. doi: 10.11833/j.issn.2095-0756.2018.06.009
Citation: XU Shiyu, SHI Yongjun, FENG Shengfei. Nondestructive and accurate measurement of volume and stem carbon storage for urban greening tree species based on terrestrial laser scanning point cloud[J]. Journal of Zhejiang A&F University, 2018, 35(6): 1062-1069. doi: 10.11833/j.issn.2095-0756.2018.06.009
  • 城市森林是森林碳储量的重要组成部分之一,对调节局部乃至全球水平的气候变化起着重要作用。绿化树种在城市森林的碳吸收、碳固存中占主导地位,其中树干约占乔木总碳储量的58.2%[1]。城市绿化树种的树干干形在复杂的环境条件和频繁的人为干扰下相比野外林分存在一定差异,传统的基于野外林分调查拟合的材积模型及树干碳储量模型并不能很好地与之相匹配[2-3];并且城市绿化树种因具有生态景观、保护生物多样性等生态功能[4-6],观赏经济价值较高,不宜通过伐倒收获进行研究。因此,无损精确地测算立木材积及树干碳储量对于评估城市森林碳汇功能以及城市森林对于缓减气候变化的贡献具有重要意义。三维激光扫描技术又称地面激光扫描技术(terrestrial laser scanning,TLS),是一种能够采集目标物体表面的三维坐标和反射强度等数据,并最终以点云形式表达在计算机上的技术[7-10]。在现有的TLS在林业调查领域的研究中,LIANG等[11]、梁子瑜[12]、尤磊[13]虽然利用建模法[13-15]和干曲线法[11-12]测算了干形通直的杉木Cunninghamia lanceotala材积,并通过与实际量测或通过材积方程测算得到的材积进行简单的比较分析,揭示了TLS在树干材积测算方面的巨大潜力。采用不同分段间隔对立木材积测算精度的影响,以及TLS在测算树干碳储量方面的潜力等研究并未有后续深入地展开。无患子Sapindus mukorossi为落叶阔叶乔木,树干通直、枝叶广展,是中国长江流域以南地区极为常见的城市绿化树种之一。本研究基于TLS结合局部采样的方法对无患子立木材积以及树干碳储量进行无损测算研究,旨在探索TLS在测算城市绿化树种材积及树干碳储量的方面潜力,分析不同区分段长度对立木材积及碳储量计算精度的影响,找出适用于材积测算的区分段长度并分析城市绿化树种树干与野外林分的异同,为中国城市绿化树种固碳能力的测算及生物量扩展因子(BEF)的研究提供依据。

    • 研究区域位于浙江省杭州市临安区的浙江农林大学校园(29°56'~30°23'N,118°51'~119°52' E)。研究区为典型的亚热带季风气候,四季分明;年平均气温为16.0 ℃,年降水量为1 614.0 mm,无霜期为237.0 d,全年光照充足。研究区有无患子,银杏Ginkgo biloba,玉兰Magnolia denudata等城市常用绿化树种;植被覆盖率高,种类丰富,人为干扰强,是典型的城市绿地类型。

    • 使用的Leica ScanStation C5脉冲式三维激光扫描仪主要通过激光器不断发射脉冲信号,通过接收接触物体表面后返回的信号计算飞行时间,结合相应信号发射时的水平角和竖直角,得到并记录目标的三维坐标点云数据。TLS特点:①精准高效。TLS在低精度下完成单测站360°×270°的全景扫描仅需2~3 min。TLS在1.50~50.00 m测距范围内采集点云数据的点位精度为±6.0 mm,测距精度为±1.0 mm。②无损。TLS通过非接触性激光对目标进行测量,不会对立木以及周边的环境产生任何影响及破坏,可以实现对立木的动态监测。③信息丰富。TLS是对目标表面的致密精细化测量,得到的点云数据具有极高的还原度,可以提取树高、胸径等丰富的测树信息。

    • 于2017年7月利用Leica ScanStation C5脉冲式三维激光扫描仪进行了单立木三维点云的采集工作。在进行点云数据采集前,以5.00 cm为起测胸径分别用胸径尺和全站仪悬高测量法对实验区内所有无患子的胸径、枝下高、树高进行测量,按《原木材积表》中所规定的2.00 cm为一径阶对无患子进行径阶的划分与统计,采用径阶等比法按比例从各径阶中选取无患子共计30株进行实验(径阶分布如图 1),由于径阶10.00 cm以下(18.00 cm以上)的样本过少,因此将10.00 cm以下(18.00 cm以上)的立木划分到18.00 cm(18.00 cm)径阶中。在选取无患子样本后,进行点云数据的采集。采集方法:本研究中为获取完整的单立木点云数据,无患子都经过3~4站·株-1的扫描,相应测站的布设近似呈等边三角形或矩形分布。并且根据地形及布设测站的通视情况,在目标立木周围放置3个公共参考标靶,确保在不同的测站都能被完整扫描到。在完成测站及标靶地布设后,进行仪器地架设并整平,并将扫描参数设置如下:扫描范围采用单立木窗扫描,即用户自定义扫描范围;采用低分辨率,并在设置相应的保存文件夹后开始扫描。扫描过程中平均扫描时间约为10 min·株-1·站-1[16]

      图  1  无患子样木的径阶分布

      Figure 1.  Diameter grade distribution of Sapindus mukorossi samples

      待样本点云数据采集完成后,在各径阶内选取3株标准木(若数量不满3株则取该径阶内全部立木),利用直径5.15 mm的瑞典手持式取芯器CO250分别在树干基部处、树高1.00 m处及1.30 m对样木钻芯进行取样,并用游标卡尺量测记录钻孔孔径[17]。将样本分径阶装袋带回实验室中,置于100 ℃下烘干至恒量(间隔2.0 h 2次称量所得的样本质量差不超过0.5%即视为恒量),用0.000 1 g精密电子天平称取各径阶的钻芯总干质量[18];利用取芯器标定的直径与记录孔径的均值作为圆柱模型的截面直径计算钻芯体积,将计算得到的均值作为该径阶的树干生物量密度。

    • 点云数据预处理及部分测树因子的提取在Cyclone软件中完成。点云数据的预处理主要包括以下几个方面:①坐标配准。利用每个测站所采集的3个公共标靶的信息可将不同测站的单立木(枝叶样本)点云数据匹配到同一坐标系下。②点云拼接。通过软件中点云拼接(unify cloud)功能将原本各测站独立的点云数据拼接成一个完整点云,并舍去重复的点和距离小于设定阈值的点(默认1.0 mm),从而减小数据冗余。点云拼接误差为1.0~2.0 mm。③点云去噪。依靠人工判别剔除在单立木点云数据采集过程中产生的噪声点、孤立点以及无用的地面点。

      提取的单立木基本测树因子包括树高、枝下高以及胸径。其中:树高采用单立木三维点云坐标中的高程直方图提取;枝下高通过在点云图像中依靠人工判别获取树冠起始位置的高程值并结合高程直方图提取[19-20];干径、胸径则采用最小二乘圆拟合算法对立木胸高处(距水平地面1.30 m)厚度为2.00 cm的切片点云进行拟合后提取[21-23]

    • 基于三维激光点云的材积测算方法采用精细区分段求积法:将树干点云以固定间距进行分段,通过最小二乘圆拟合算法提取各分段截面直径,采用圆台模型测算各分段体积,汇总后得到立木材积VS。碳储量CS则用材积VS与树干单位体积生物质量ρS及含碳率的乘积进行测算。

      $$ \begin{array}{l} {V_{\rm{S}}} = \sum\limits_{i = 1}^{n - 1} {\frac{1}{3}} {\rm{ \mathsf{ π} }}{h_i}\left( {R_i^2 + r_{i + 1}^2 + R \times r} \right);\\ \;\;\;\;\;\;\;\;\;\;{C_{\rm{S}}} = {V_{\rm{S}}} \times {\rho _{\rm{S}}} \times 0.48。 \end{array} $$

      其中:VS表示立木材积;hi表示分段间隔;Ri表示下截面半径;ri表示上截面半径。

    • 采用浙江省全省通用的阔叶树二元材积模型计算立木材积模型值VS2;而一元材积模型因受地域影响较大,因此,采用浙江省浙西北丘陵山区专用的一元材积模型计算立木材积模型值VS1[24],计算公式如下所示。碳储量测算方法则与上述2.4节中基于三维点云的树干碳储量测算方法相同。

      $$ \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{V_{{\rm{S2}}}} = 6.7453 \times {10^{ - 5}}{D^{1.96}}{H^{0.8144}};\\ {V_{{\rm{S1}}}} = 5.0479 \times {10^{ - 5}}{D^{1.9085}}{\left( {37.3295 + \frac{{ - 2.053.1848}}{{D + 60}}} \right)^{0.9908}}。 \end{array} $$

      其中:D为胸径;H为树高。

      树干碳储量参考值CS采用基于第7次全国森林资源清查资料的《中国森林植被生物量和碳储量评估》中的其他硬阔叶树类单立木碳储量模型。计算公式如下:

      $$ \begin{array}{l} {C_{{\rm{MS}}}} = 0.044{\left( {{D^2}H} \right)^{0.9169}} \times 0.48;\\ {C_{{\rm{MP}}}} = 0.023{\left( {{D^2}H} \right)^{0.7115}} \times 0.48;\\ \;\;\;\;\;\;\;\;\;\;\;\;{C_{\rm{S}}} = {C_{{\rm{MS}}}} + {C_{{\rm{MP}}}}。 \end{array} $$

      其中:D为胸径;H为树高;CMS表示树干碳储量(去皮);CMP表示树皮碳储量。

    • 使用Excel 2010和SPSS 20.0进行单因素方差分析(one-way ANOVA)检验不同处理方法间的差异显著性,显著性水平设定为α=0.05。并采用标准误差(ERMSE)和乖离率(sbias)检测不同处理方法间的误差。计算公式如下:

      $$ \begin{array}{l} {E_{{\rm{RMSE}}}}{\rm{ = }}\sqrt {\frac{{\sum\limits_{i = 1}^n {{{\left( {{y_i} - {{\hat y}_i}} \right)}^2}} }}{n}} ;\\ \;\;\;{s_{{\rm{bias}}}} = \frac{{\sum\limits_{i = 1}^n {\left( {{y_i} - {{\hat y}_i}} \right)} }}{n}。 \end{array} $$

      其中:y表示TLS测算值;$ {{\hat y}} $表示实测值。

    • 对基于TLS测算的不同径阶共计30株无患子的树高与胸径数据与实测值进行比较(图 2),发现两者均具有极高的精度。其中,胸径标准误差(ERMSE)为±0.34 cm,平均相对误差为1.97%;树高标准误差(ERMSE)为±0.36 m,平均相对误差为4.58%。

      图  2  TLS提取的胸径、树高与实测值比较

      Figure 2.  Comparison of diameter and height extracted base on TLS with actual measurement

    • 在以0.10 m区分段长度计算的材积值作为标准值的前提下(本研究中在对不同方法测算的材积、碳储量精度的比较中,均以此为前提),由表 1可见:30株无患子利用不同区分段长度计算得到的材积均值变化范围为37.93~40.10 dm3,其中各区分段长度计算材积值的标准误差(ERMSE)变化范围从最低的0.44 dm3到最高的2.79 dm3,变化幅度并不明显,各区分段长度测算得到的材积值之间均无显著差异(P>0.05)。乖离率(sbias)变化范围则从最低的0.02 dm3到最高的2.17 dm3,表明随区分段长度增加,测算的材积值逐渐偏离真实值。

      表 1  不同区分段长度下测算的材积比较

      Table 1.  Comparison of volume measured based on different section lengths

      区分段长度/m TLS测算材积/dm3 标准误差(ERMSE)/dm3 乖离率(sbias)/dm3
      0.10 37.93 ± 10.96 a
      0.20 37.92 ± 10.84 a 0.44 0.02
      0.40 38.57 ± 10.86 a 1.08 0.49
      0.50 38.44 ± 11.29 a 1.15 0.44
      0.60 39.15 ± 11.02 a 1.69 0.92
      0.80 39.50 ± 11.24 a 2.08 1.38
      1.00 40.10 ± 11.57 a 2.79 2.17
      说明:*每个数据表示样本平均值±标准差;不同小写字母表示不同区分段测算的材积之间差异显著(P<0.05)

      图 3可见:当区分段长度≤0.40 m时,材积测算误差集中、近似呈正态分布,材积测算精度高且稳定;而当区分段长度由0.40 m变为0.50 m时,平均误差及误差离散度迅速变大;当区分段长度大于0.50 m时,随着区分段长度的增加,平均误差及误差离散度持续变大,在区分段长度为1.00 m时达到最大。说明随着区分段长度的增加,测算的立木材积更易偏离真值,出现异常情况如:区分段长度为0.60 m下的普通异常点whz10和whz17,极端异常点whz16以及区分段长度为0.80 m下的普通异常点whz19和whz26。

      图  3  不同区分段长度下材积测算误差的分布及比较

      Figure 3.  Distribution and comparison of volume measurement error in different section lengths

    • 表 2对基于4种不同方法测算的30株无患子立木材积进行方差分析,发现利用现有一元材积模型和二元材积模型测算的立木材积之间并无显著差异(P>0.05)。但基于一元材积模型和二元材积模型测算的立木材积却均要显著高于基于TLS区分求积得到的立木材积(P<0.05)。其中,利用一元材积模型测算的材积值最高(68.77 dm3 ± 36.58 dm3),分别是以0.10 m为区分段长度测算的材积值(37.93 dm3 ± 10.97 dm3)的1.67倍和以1.00 m为区分段长度测算的材积值(40.10 dm3 ± 11.57 dm3)的1.58倍,可见利用一元材积模型测算城市绿化树种材积误差之大。相比较而言,利用二元材积模型测算的材积值(54.54 dm3 ± 31.91 dm3)与以区分求积得到的材积值之间的差距要小的多。

      表 2  不同方法测算的材积、树干碳储量比较

      Table 2.  Comparison of volume and stem carbon storage measured by different methods

      测算方法 测算材积/dm3 测算碳储量/kg
      TLS_0.10 m区分求积 37.93 ± 10.97 c 12.38 ± 5.64 d
      TLS_1.00 m区分求积 40.10 ± 11.57 c 13.09 ± 5.77 cd
      一元材积模型 68.77 ± 36.58 ab 22.26 ± 11.80 a
      二元材积模型 54.54 ± 31.91 b 17.64 ± 10.25 bc
      碳储量模型 17.92 ± 9.59 bc
      说明:样本平均值±标准差;不同小写字母表示不同区分段测算的材积之间差异显著(P<0.05)

      利用5种不同方法测算的30株无患子CS之间也存在不同程度的差异:其中利用一元材积模型测算的CS(22.26 kg ± 11.80 kg)最高,并与其他方式测算的CS存在显著差异(P<0.05)。而二元材积模型以及树干碳储量模型测算的CS(17.64 kg ± 10.25 kg,17.92 kg ± 9.59 kg)与以1.00 m为区分段长度测算的CS(13.09 kg ± 5.77 kg)之间无显著差异(P>0.05),但均要显著高于以0.10 m为区分段长度测算的CS(12.38 kg ± 5.64 kg)(P<0.05),表明利用二元材积模型和碳储量模型仅能实现CS的粗略测算。本研究中利用材积模型测算CS是以材积与树干生物量密度以及含碳率的乘积实现的,但研究发现利用二元材积模型测算的CS与一元材积模型测算的CS之间存在显著差异,而与以1.00 m为区分段长度测算的CS之间并无显著差异,这与上述研究中材积显著性差异的表现并不一致。

      图 4A将现有模型测算的材积和作为标准值的基于TLS测算的材积进行比较发现:利用现有模型测算的材积在立木自身材积较小时具有较高精度,但随着立木自身材积的不断增大,模型测算的立木材积逐渐开始偏离参考线y=x,即精度不断下降,并且偏离速度随材积的增大逐渐变快。由图 4B可见:利用现有模型测算的立木CS精度在随立木自身CS变化的趋势上与利用模型测算的立木材积精度具有一致性——即利用现有模型测算的CS在立木自身CS较小时具有较高精度,但随着立木自身CS的不断增大,模型测算的CS逐渐开始偏离参考线y=x,即精度不断下降,并且偏离速度随CS的增大逐渐变快。表明利用现有模型测算的材积和CS精度随立木的生长在总体上呈下降趋势。而一元材积模型测算的材积和CS精度均为现有模型中最差。

      图  4  现有模型测算的城市绿化树种材积及树干碳储量与TLS测算结果的比较

      Figure 4.  Comparison of the urban greening tree species volume and stem carbon storage measured based on existing model

    • 基于不同区分段长度测算的立木材积值之间并没有显著差异,说明在对精度没有较高要求的研究或调查中,≤1.00 m的区分段长度均可满足立木材积测算的需要,这与传统林业伐倒木的材积测算过程中通常将树干按长短以0.50或1.00 m进行分段相一致[25]

      LIANG等[11]和SAARINEN等[26]在对高大通直的杉木材积研究中的干图像反映:立木干曲线虽然随着树高的增加总体呈减小趋势,但仍存在大量比上一级分段截面干径要大的情况,即不同的干径提取间隔会对立木材积测算造成一定影响。本研究通过以0.10 m区分段长度测算的立木材积为真值,对比分析各区分段下材积测算误差曲线图,发现随区分段长度的增加材积测算的平均误差及误差的离散度也会随之增加,主要是因为随着区分段长度的增加就会忽略越多的树干形变现象,而该现象在区分段长度由0.40 m增加0.50 m时尤其明显。这表明当区分段长度大于0.40 m时,测算的立木材积与真值之间的差距会明显增大,这也为TLS在如何精细高效分段测算材积提供了重要补充。

      对比分析4种不同方法测算的立木材积,发现利用现有二元或一元材积模型测算的材积在总体上都要显著高于TLS区分求积值,并且测算的材积精度随立木胸径的增大的生长降低。这说明城市绿化树种干形与野外林分相比存在明显差异,而且随林龄的增大,其干形差异逐步趋大,现有的材积模型并不能很好地与之相匹配。史琰[17]、嵇浩翔等[27]也发现城市中复杂的环境条件、频繁的人为干扰,在不同程度上促进或抑制着立木的生长,使其生长态势相比野外林分存在极大的不确定性,如树干起始分叉位置的复杂多变,并且这种现象会随立木的不断生长而累积,与本研究结果具有一致性。

      此外,本研究对5种不同CS测算方法的分析研究表明:基于一元材积模型测算CS会产生严重高估,主要是由测算立木材积过程中仅用胸径单一因子无法得到精确的材积数据导致的,因此,不建议利用一元材积模型测算城市绿化树种的立木材积及碳储量;基于二元材积模型以及树干碳储量模型测算的CS与以1.00 m为区分段长度测算的CS无显著差异,但相比以0.10 m为区分段长度测算的CS均显著偏高,说明两者作为常用的树干碳储量测算模型仅能对城市绿化树种树干碳储量进行粗略的测算,但与真实值存在一定的差距。该问题可能主要是由两大原因所引起,其一是由于存在材积测算误差,其二则是由于不同树种的树干单位体积碳储量存在差异,甚至相同树种不同径阶之间的树干单位体积碳储量也存在微小的差异,而现有的碳储量模型为保证广泛通用性大多采用多种阔叶树的均值来拟合构建。此外,不同碳储量测算方法间的显著性差异与材积间的显著性差异并非一致也说明CS是由材积及树干单位体积碳储量共同决定的。

    • 基于TLS提取的胸径、树高能达到极高的精度。通过对比分析30株无患子基于TLS胸径和树高与外业实测值,其平均相对误差分别仅为1.97%和4.85%。表明基于TLS材积的单立木点云数据具有较高的还原度和客观真实等特点,也表明基于TLS提取的基本测树因子可在无需伐倒立木的情况下精确测算立木的材积,并结合树干钻芯取样,在对立木的正常生长不造成影响的前提下精确地测算树干碳储量。

      不同的区分段长度对基于TLS测算的树干材积结果虽无显著影响,但随着区分段长度的增加,材积的测算精度仍会不断下降。因此,在对树干材积精度有较高要求的研究中,建议利用≤0.40 m的区分段长度进行树干材积的测算。该结论同样适用于伐倒木的材积精确测算。

      城市绿化树种的主干材积及碳储量受环境因子和人为干扰影响与野外林分相比较存在显著差异。现有材积模型及碳储量模型仅能用于粗略测算城市绿化树种的材积及树干碳储量,但并不能到达高度匹配,尤其是一元材积模型。

参考文献 (27)

目录

    /

    返回文章
    返回