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

林业科学 ›› 2021, Vol. 57 ›› Issue (6): 14-23.doi: 10.11707/j.1001-7488.20210602

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

基于针刺仪平稳卡尔曼滤波器的树木年龄估计

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

  1. 1. 中国林业科学研究院资源信息研究所 北京 100091
    2. 信阳师范学院数学与统计学院 信阳 464000
  • 收稿日期:2019-09-19 出版日期:2021-06-25 发布日期:2021-08-06
  • 通讯作者: 唐守正
  • 基金资助:
    国家重点研发计划项目"长白山次生林抚育更新技术研究与示范"(2017YFC0504101);中央级公益性科研院所基本科研业务费专项资金重点项目(CAFYBB2018SZ007);"十三五"国家重点研发计划项目(2017YFD060040302);青年科学基金项目(31100475);国家自然科学基金面上项目(31872704);国家自然科学基金面上项目(31470641)

Tree Age Estimation Based on Resistograph Stationary Kalman Filter

Hong Pan1,2,Jun Lu1,Xiangdong Lei1,Xuzhan Guo1,Jianfeng Yao1,Shouzheng Tang1,*   

  1. 1. Institute of Forest Resource Information Techniques, CAF Beijing 100091
    2. College of Mathematics and Statistics, Xinyang Normal University Xinyang 464000
  • Received:2019-09-19 Online:2021-06-25 Published:2021-08-06
  • Contact: Shouzheng Tang

摘要:

目的: 研究卡尔曼滤波在针刺仪抗钻阻力值序列中的应用,推导出针刺仪平稳卡尔曼滤波器估计树木年龄,为活立木年龄估计提供方法和依据。方法: 以山西省羊圈沟林场华北落叶松为研究对象,使用针刺仪在同一水平位置不同方向钻入华北落叶松获取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软件自动判定结果,具有微损、操作快、测量精度高等特点,可为华北落叶松年龄估计提供一种有效途径。

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

Abstract:

Objective: The application of Kalman filter in the resistograph drilling resistance sequence was studied. The resistograph stationary Kalman filter was derived,which provided a method and basis for estimating the age of standing trees. Method: Taking Larix principis-rupprechtii in Yangjuangou forest farm of Shanxi Province as the research object,323 sets of drilling resistance series obtained by resistograph drilling L. principis-rupprechtii in the same horizontal position and different directions were used as the research samples,and 104 stem discs were taken within 5 cm of the drilling position as reference samples. The general Kalman filter theory was applied in resistograph. By analyzing the characteristics of the filter,the parameters were simplified,and the stable Kalman filter was derived. The drilling resistance data obtained through the resistograph drilling into larch was used as input of stationary Kalman filter,and the output was the optimal estimation sequence based on the minimum mean square error as the best criterion,and the tree age was determined by the number of peaks and troughs of the optimal estimation sequence. The results estimated by the stationary Kalman filter algorithm,or auto-recognized by DECOM software were analyzed with the corresponding number of disc rings using paired data t-test. Result: The resistograph stationary Kalman filter reduced the number of initial parameters,the optimal estimation sequence of drilling resistance value only depended on parameter Rat. When the appropriate parameter Rat was selected according to their DBH size,the resistograph stationary Kalman filter applied in 323 groups of drilling resistance standing had better denoising effects,so that the algorithm estimating tree age was closed to the actual age. The relative error distribution of the algorithm was mostly between -10% and 10%,the minimum relative error was 0%,the maximum relative error was 25.69%,and the average relative error was 0.75%. Through paired data t-test,it was found that t-value was -0.468 16,indicating that there was no significant difference between the mean age estimated by the algorithm and the true mean age of trees. The relative error of the automatic analysis by DECOM software was large,mostly concentrated between -20% and -60%. The minimum relative error was -7.69%,the maximum relative error was -84.78%,and the average relative error was -40.49%. Through paired data t-test,it was found that t-value was 20.254,indicating that there was significant difference between the mean age estimated by DECOM software and the mean tree age of trees. Conclusion: The method of applying the resistograph stationary Kalman filter to the drilling resistance data sequence could accurately estimate the age of L. principis-rupprechtii,which might be better than the automatic analysis by DECOM software. This method is small invasive,fast and high accuracy,which provides an effective way for estimating the age of standing L. principis-rupprechtii.

Key words: Larix principis-rupprechtii, tree age, resistograph, drilling resistance data, the stationary Kalman filter, posterior estimation

中图分类号: