林业科学

作者：潘虹1,2,卢军1,雷相东1,郭旭展1,姚建峰1,唐守正1,*

单位：1. 中国林业科学研究院资源信息研究所 北京 1000912. 信阳师范学院数学与统计学院 信阳 464000

关键词：华北落叶松,树木年龄,针刺仪,抗钻阻力值,平稳卡尔曼滤波,后验估计值

基金：国家重点研发计划项目"长白山次生林抚育更新技术研究与示范"(2017YFC0504101);中央级公益性科研院所基本科研业务费专项资金重点项目(CAFYBB2018SZ007);"十三五"国家重点研发计划项目(2017YFD060040302);青年科学基金项目(31100475);国家自然科学基金面上项目(31872704);国家自然科学基金面上项目(31470641)

目的: 研究卡尔曼滤波在针刺仪抗钻阻力值序列中的应用，推导出针刺仪平稳卡尔曼滤波器估计树木年龄，为活立木年龄估计提供方法和依据。方法: 以山西省羊圈沟林场华北落叶松为研究对象，使用针刺仪在同一水平位置不同方向钻入华北落叶松获取323组抗钻阻力值序列作为研究样本，在针刺位置5 cm内截取104个圆盘作为参考样本。将卡尔曼滤波器理论应用于针刺仪，通过分析滤波器的性质，简化参数，推导出针刺仪平稳卡尔曼滤波器。以针刺仪钻入华北落叶松获得的抗钻阻力值作为平稳卡尔曼滤波器的输入，以最小均方误差为最佳准则的最优估计值序列作为输出，根据最优估计值序列的波峰波谷数估计树木年龄，并将平稳卡尔曼滤波算法估计的树木年龄和针刺仪自带DECOM软件自动判定的树木年龄分别与相应的圆盘年轮数进行成对数据t-检验。结果: 针刺仪平稳卡尔曼滤波器缩减初始参数数量，使得抗钻阻力的最优估计值序列只依赖于参数Rat值的选取。将323组活立木抗钻阻力值序列经过针刺仪平稳卡尔曼滤波器，依据相应活立木胸径大小选择恰当的参数Rat值可以得到较好去噪效果，算法估计树木年龄与实际树木年龄比较接近。算法估计树木年龄相对误差分布大多集中在-10%~10%之间，最小相对误差为0%，最大相对误差为25.69%，平均相对误差为0.75%，经成对数据t-检验后t为-0.468 16，算法估计树木年龄均值与实际树木年龄均值之间无显著差异。而针刺仪自带DECOM软件自动判定的树木年龄结果相对误差较大，大多集中在-20%~-60%之间，最小相对误差为-7.69%，最大相对误差为-84.78%，平均相对误差为-40.49%，与实际树木年龄进行成对数据t-检验后t为20.245，差异显著。结论: 针刺仪平稳卡尔曼滤波器应用于针刺仪抗钻阻力值序列能够较准确估计华北落叶松年龄，优于针刺仪自带DECOM软件自动判定结果，具有微损、操作快、测量精度高等特点，可为华北落叶松年龄估计提供一种有效途径。

来源：2021年第6期

《林业科学》期刊编辑部

查看林业科学杂志2021年第6期