高级检索

留言板

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

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

基于速度误差校正的林木应力波无损检测断层成像算法

焦治 李光辉 武夕

焦治, 李光辉, 武夕. 基于速度误差校正的林木应力波无损检测断层成像算法[J]. 北京林业大学学报, 2018, 40(1): 108-119. doi: 10.13332/j.1000-1522.20170384
引用本文: 焦治, 李光辉, 武夕. 基于速度误差校正的林木应力波无损检测断层成像算法[J]. 北京林业大学学报, 2018, 40(1): 108-119. doi: 10.13332/j.1000-1522.20170384
Jiao Zhi, Li Guang-hui, Wu Xi. Tomography imaging algorithm based on velocity error correction for stress wave nondestructive evaluation of wood[J]. Journal of Beijing Forestry University, 2018, 40(1): 108-119. doi: 10.13332/j.1000-1522.20170384
Citation: Jiao Zhi, Li Guang-hui, Wu Xi. Tomography imaging algorithm based on velocity error correction for stress wave nondestructive evaluation of wood[J]. Journal of Beijing Forestry University, 2018, 40(1): 108-119. doi: 10.13332/j.1000-1522.20170384

基于速度误差校正的林木应力波无损检测断层成像算法

doi: 10.13332/j.1000-1522.20170384
基金项目: 

江苏省重点研发计划项目 BE2016627

中央高校基本科研业务费专项 RP51635B

无锡市国际科技研发合作项目 CZE02H1706

江苏省普通高校专业学位研究生实践创新计划项目 SJLX16_0499

国家自然科学基金项目 61472368

详细信息
    作者简介:

    焦治。主要研究方向:智能无损检测技术。Email:novrule@163.com 地址: 214122 江苏省无锡市江南大学物联网工程学院

    通讯作者:

    李光辉,博士,教授。主要研究方向:传感器网络、容错计算和智能无损检测技术。Email: ghli@jiangnan.edu.cn 地址:同上

  • 中图分类号: S781;TS67

Tomography imaging algorithm based on velocity error correction for stress wave nondestructive evaluation of wood

  • 摘要: 目的应力波断层成像技术已经在林木无损检测领域得到了广泛的应用, 然而木材的各向异性特征以及应力波速度反演计算的误差对断层成像的精度影响较大。因此,进一步提高断层图像的准确性非常关键。方法根据传统的应力波速度反演原理, 提出了一种基于速度误差校正的断层成像算法(ECIA)。该算法利用最小二乘QR分解法(LSQR)计算应力波在林木横截面网格单元内的速度分布, 并使用误差校正机制(ECM)优化断层图像。为了评估算法的性能, 分别选取了若干实验室内的原木试样及扬州市区古树进行无损检测实验, 利用德国PICUS应力波断层成像仪获取的应力波传播数据实现了ECIA算法, 生成了各实验样本的断层图像。结果ECIA算法较为精确地检测出了原木及活树中缺陷区域的位置及大小, 尤其在活树的健康检测中, ECIA算法对于缺陷区域尤其是轻微腐朽区域的检测精度高于PICUS检测仪。结论ECIA算法能够产生较为准确的林木断层图像,适用于林木应力波无损检测。
  • 图  1  树木横截面内应力波传播示意图

    S代表应力波发射端传感器的位置,R1R2代表应力波接收端传感器的位置;VRVT分别表示径向传播速度和弦向传播速度;β表示弦向角。

    Figure  1.  Stress wave propagation in the cross-section of tree trunk

    S is a source sensor, and R1 and R2 are receiver sensors. VT is the tangential velocity and VR is the radial velocity. β is the angle between radial direction and tangential direction.

    图  2  树木横截面网格单元划分示意图

    SR分别代表应力波发射端传感器和应力波接收端传感器的位置;1, 2, …, M代表网格单元编号;P1, P2, P3, P4代表传播路径与网格单元的交点;dij代表第i条路径被第j个网格单元分割的距离。

    Figure  2.  Grid cells of a cross-section of a tree trunk

    S is a source sensor, and R is a receiver sensor. 1, 2, …, M represent the No. of grid cells. P1, P2, P3, P4 mean the intersections of stress wave paths and grid cells. dij is segment intercepted by the ith wave path on the jth cell.

    图  3  树干横截面内网格单元的分类

    网格单元内的变量或数字代表该单元的编号。

    Figure  3.  Grid cell classification of a cross-section of a tree trunk

    The variable or number represent the No. of grid cells.

    图  4  活树应力波无损检测现场

    Figure  4.  Stress wave testing of a live standing tree

    图  5  试样1的断层图像比较

    Figure  5.  Comparison of tomography images of sample 1

    图  6  试样2的断层图像比较

    Figure  6.  Comparison of tomography images of sample 2

    图  7  试样3的断层图像比较

    Figure  7.  Comparison of tomography images of sample 3

    图  8  试样4的断层图像比较

    Figure  8.  Comparison of tomography images of sample 4

    图  9  试样4的阻力曲线图

    Figure  9.  Resistance curves of sample 4

    图  10  试样5的断层图像比较

    Figure  10.  Comparison of tomography images of sample 5

    图  11  试样5的阻力曲线图

    Figure  11.  Resistance curves of sample 5

    图  12  试样6的断层图像比较

    Figure  12.  Comparison of tomography images of sample 6

    图  13  试样6的阻力曲线图

    Figure  13.  Resistance curves of sample 6

    图  14  试样7的断层图像比较

    Figure  14.  Comparison of tomography images of sample 7

    图  15  试样7的阻力曲线图

    Figure  15.  Resistance curves of sample 7

    图  16  试样8的断层图像比较

    Figure  16.  Comparison of tomography images of sample 8

    图  17  试样8的阻力曲线图

    Figure  17.  Resistance curves of sample 8

    图  18  试样9的断层图像比较

    Figure  18.  Comparison of tomography images of sample 9

    图  19  试样9的阻力曲线图

    Figure  19.  Resistance curves of sample 9

    表  1  原木断层成像结果的定量评价

    Table  1.   Quantitative evaluation to tomography images of logs

    试样编号
    No. of
    sample
    试样横截面积
    Cross-sectional area
    of the sample/cm2
    真实S
    Actual S/cm2
    真实RS
    Actual RS/%
    ECIAPICUS
    S/cm2RS/%ΔRS/%S/cm2RS/%ΔRS/%
    1812.593 5297.284 136.584 6324.558 839.941 1-9.174 6232.514 728.613 921.787 6
    2794.226 078.019 29.823 3102.108 912.856 4-30.876 6285.921 436.000 0-266.475 6
    3804.247 70000088.467 211.000 0>0
    注:S代表试样或断层图像中的缺陷面积;RS代表试样或断层图像中缺陷面积与横截面积的比值;ΔRS代表断层图像的RS与真实RS的相对误差。Notes:S is the defective area, RS is the ratio of the defective area to the cross-section area, and ΔRS is the relative error between the RS of ECIA or PICUS and the actual RS.
    下载: 导出CSV

    表  2  活树断层成像结果的定量评价

    Table  2.   Quantitative evaluation to tomography images of standing trees

    试样编号
    Sample No.
    路径
    Path
    真实RL
    Actual RL
    ECIAPICUS
    RLΔRLRLΔRL
    42-1017.910 417.254 13.664 412.537 329.999 9
    2-913.432 812.455 27.277 70100.000 0
    5-15.970 16.275 8-5.120 518.209 0-205.003 2
    53-888.549 687.042 31.702 267.605 623.652 3
    4-783.333 377.500 07.000 053.928 635.285 7
    5-1074.626 876.417 9-2.400 168.656 77.999 9
    69-175.555 593.043 5-23.145 943.478 342.455 1
    2-1184.444 488.888 9-5.263 288.888 9-5.263 2
    71-789.937 193.548 4-4.015 495.913 9-6.645 5
    2-975.496 781.073 5-7.386 856.989 224.514 3
    4-1091.463 493.566 7-2.299 680.444 412.047 4
    85-93.125 013.437 5-330.000 014.375 0-360.000 0
    3-93.703 74.521 1-22.069 812.592 6-240.000 0
    910-632.558 131.395 33.571 527.665 615.027 0
    10-527.907 027.702 30.733 536.511 6-30.833 1
    7-1246.774 248.387 1-3.448 343.037 77.988 4
    注:RL代表试样或断层图像中路径上缺陷区域的长度与整个路径长度的比值;ΔRL代表断层图像的RL与真实RL的相对误差。Notes: RL is the ratio of the defective length to the whole length of the path under test, and ΔRL is the relative error between the RL of ECIA or PICUS and the actual RL.
    下载: 导出CSV
  • [1] Wang X P, Divos F, Pilon C, et al. Assessment of decay in standing timber using stress wave timing nondestructive evaluation tools: a guide for use and interpretation[R]. General Technical Report FPL-GTR-147. Madison, WI: United States Department of Agriculture, 2004.
    [2] Dikrallah A, Hakam A, Brancheriau L, et al. Experimental analysis of acoustic anisotropy of wood by using guided waves[C]//Proceedings of International Conference on Integrated Approach to Wood Structure, Behaviour and Application, Joint Meeting of ESWM and COST Action E35. Florence, Italy: Aalborg Universitet, 2006: 149-154.
    [3] 梁善庆, 王喜平, 蔡智勇, 等.弹性波层析成像技术检测活立木腐朽[J].林业科学, 2008, 44(5): 109-114. http://d.old.wanfangdata.com.cn/Periodical/lykx200805021

    Liang S Q, Wang X P, Cai Z Y, et al. Elastic wave tomography in standing tree decay detection[J]. Scientia Silvae Sinicae, 2008, 44(5): 109-114. http://d.old.wanfangdata.com.cn/Periodical/lykx200805021
    [4] 杨学春, 王立海.应力波在原木中传播理论的研究[J].林业科学, 2005, 41(5): 132-138. http://d.old.wanfangdata.com.cn/Periodical/lykx200505024

    Yang X C, Wang L H. Study on the propagation theories of stress wave in log[J]. Scientia Silvae Sinicae, 2005, 41(5): 132-138. http://d.old.wanfangdata.com.cn/Periodical/lykx200505024
    [5] Li G H, Wang X P, Feng H L, et al. Analysis of wave velocity patterns in black cherry trees and its effect on internal decay detection[J]. Computers and Electronics in Agriculture, 2014, 104: 32-39. doi:  10.1016/j.compag.2014.03.008
    [6] Socco L V, Sambuelli L, Martinis R, et al. Feasibility of ultrasonic tomography for nondestructive testing of decay on living trees[J]. Research in Nondestructive Evaluation, 2004, 15(1): 31-54. doi:  10.1080/09349840490432678
    [7] Du X C, Li S Z, Li G H, et al. Stress wave tomography of wood internal defects using ellipse-based spatial interpolation and velocity compensation[J]. BioResources, 2015, 10(3): 3948-3962.
    [8] Guntekin E, Emiroglu Z G, Yilmaz T. Prediction of bending properties for Turkish red pine (Pinus brutia Ten.) lumber using stress wave method[J]. BioResources, 2013, 8(1): 231-237. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=Doaj000003296196
    [9] Dackermann U, Crews K, Kasal B, et al. In situ assessment of structural timber using stress-wave measurements[J]. Materials and Structures, 2014, 47(5): 787-803. doi:  10.1617/s11527-013-0095-4
    [10] Divos F, Divos P. Resolution of stress wave based acoustic tomography[C]//Proceedings of the 14th International Symposium on Nondestructive Testing of Wood. Eberswalde: Shaker Verlag, 2005: 309-314.
    [11] Maurer H R, Schubert S, Baechle F. Application of nonlinear acoustic tomography for non-destructive testing of trees[C]//Proceedings of 14th International Symposium on Nondestructive Testing of Wood. Eberswalde: Shaker Verlag, 2005: 339-349.
    [12] Brancheriau L, Lasaygues P, Debieu E, et al. Ultrasonic tomography of green wood using a non-parametric imaging algorithm with reflected waves[J]. Annals of Forest Science, 2008, 65(7): 712-717. doi:  10.1051/forest:200851
    [13] Jackson M J, Tweeton D R. MIGRATOM: geophysical tomography using wavefront migration and fuzzy constraints[R]//United States, report of investigation: RI 9497. Minneapolis, MN: U.S. Department of the Interior, Bureau of Mines, 1994.
    [14] Paige C C, Saunders M A. LSQR: an algorithm for sparse linear equations and sparse least squares[J]. ACM Transactions on Mathematical Software (TOMS), 1982, 8(1): 43-71. doi:  10.1145/355984.355989
    [15] 杨国铮, 禹晶, 肖创柏, 等.基于形态成分分析的复杂背景SAR图像舰船尾迹检测[J].计算机辅助设计与图形学学报, 2016, 28(10): 1662-1671. http://d.old.wanfangdata.com.cn/Periodical/jsjfzsjytxxxb201610007

    Yang G Z, Yu J, Xiao C B, et al. Ship wake detection in SAR images with complex backgrounds using morphological component analysis[J]. Journal of Computer-Aided Design & Computer Graphics, 2016, 28(10): 1662-1671. http://d.old.wanfangdata.com.cn/Periodical/jsjfzsjytxxxb201610007
    [16] 谭茂金, 邹友龙. (T2, D)二维核磁共振测井混合反演方法与参数影响分析[J].地球物理学报, 2012, 55(2): 683-692. http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201202032

    Tan M J, Zou Y L. A hybrid inversion method of (T2, D) 2D NMR logging and observation parameters effects[J]. Chinese Journal of Geophysics, 2012, 55(2): 683-692. http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201202032
    [17] 程绪琦, 王建华, 刘志峰, 等. AutoCAD 2014中文版标准教程[M].北京:电子工业出版社, 2014: 182-184.

    Cheng X Q, Wang J H, Liu Z F, et al. The standard tutorial of AutoCAD 2014 in Chinese version[M]. Beijing: Publishing House of Electronics Industry, 2014: 182-184.
  • [1] 魏喜雯, 孙丽萍, 许述正, 杨扬, 杜春晓.  基于应力波传播速度模型的原木缺陷定量检测 . 北京林业大学学报, 2020, 42(5): 143-154. doi: 10.12171/j.1000-1522.20190420
    [2] 管成, 刘晋浩, 张厚江, 周卢婧.  足尺人造板力学性能无损检测研究进展 . 北京林业大学学报, 2019, 41(9): 164-172. doi: 10.13332/j.1000-1522.20180379
    [3] 陈勇平, 高甜, 李德山, 郭文静.  马尾松木材空洞开裂残损的雷达检测成像初探 . 北京林业大学学报, 2017, 39(3): 112-118. doi: 10.13332/j.1000-1522.20160313
    [4] 孟兆新, 肖定福, 尹航, 李尚, 张聪.  并联式曲线送料平台机构误差补偿方法的研究 . 北京林业大学学报, 2016, 38(9): 95-101. doi: 10.13332/j.1000-1522.20160026
    [5] 伏静, 戴思兰.  基于高光谱成像技术的菊花花色表型和色素成分分析 . 北京林业大学学报, 2016, 38(8): 88-98. doi: 10.13332/j.1000-1522.20150483
    [6] 张怡卓, 马琳, 许雷, 于慧伶.  基于小波与曲波遗传融合的木材纹理分类 . 北京林业大学学报, 2014, 36(2): 119-124.
    [7] 廖春晖, 张厚江, 王喜平, 黎冬青, 祖汉松, 张文博.  明十三陵裕陵修缮工程中的木构件无损检测 . 北京林业大学学报, 2014, 36(1): 132-137.
    [8] 刘昊, 高建民.  含水率和密度对木材应力波传播速度的影响 . 北京林业大学学报, 2014, 36(6): 154-158. doi: 10.13332/j.cnki.jbfu.2014.06.002
    [9] 宋文龙, 杨鑫, 李克新, 贾鹤鸣.  基于梯度幅度和Hough 变换的林木根系GPR 目标检测 . 北京林业大学学报, 2013, 35(6): 108-112.
    [10] 张超, 彭道黎, 涂云燕, 党永峰, 智长贵.  利用TM 影像和偏最小二乘回归方法估测三峡库区森林蓄积量 . 北京林业大学学报, 2013, 35(3): 11-17.
    [11] 胡雄武, 张平松, 于仲.  直流电法CT技术检测林木孔洞结构特征 . 北京林业大学学报, 2012, 34(3): 103-106.
    [12] 黄妍, 张厚江, 刘妍, J, ohnF.Hunt, 闫海成.  基于悬臂梁弯曲的薄板木质材料应力松弛检测 . 北京林业大学学报, 2012, 34(4): 149-153.
    [13] 张厚江, 朱磊, 孙燕良, 王喜平.  古建筑木构件材料主要力学性能检测方法研究 . 北京林业大学学报, 2011, 33(5): 126-129.
    [14] 王天龙, 陈永平, 刘秀英, 姜恩来, 余如龙.  古建筑木构件缺陷及评价残余弹性模量的初步研究 . 北京林业大学学报, 2010, 32(3): 141-145.
    [15] 张厚江, 王喜平, 苏娟, RobertJRoss.  应力波在美国红松立木中传播机理的试验研究 . 北京林业大学学报, 2010, 32(2): 145-148.
    [16] 倪成才, 刘春梅, 丁俊峰, 潘晓茹.  差分生长模型预测误差的分析 . 北京林业大学学报, 2009, 31(4): 1-6.
    [17] 王立海, 王洋, 高珊, 徐华东, 杨学春.  冻结状态下应力波在长白落叶松立木中传播速度的研究 . 北京林业大学学报, 2009, 31(3): 96-99.
    [18] 符韵林, 赵广杰, .  二氧化硅/木材复合材料的应力松弛 . 北京林业大学学报, 2008, 30(1): 119-123.
    [19] 罗彬, 殷亚方, 姜笑梅, 骆秀琴, 刘波, 郭起荣, .  3种无损检测方法评估巨尾桉木材抗弯和抗压强度性质 . 北京林业大学学报, 2008, 30(6): 137-140.
    [20] 王兰珍, 吴丽娟, 李雪华, 赵铁珍, 周传艳, 李义良, 李生宇, 耿玉清, 方升佐, 王立海, HUALi_zhong, 张冬梅, 党文杰, 刘鹏举, 朱小龙, 薛康, 李国雷, 刘剑锋, 韦艳葵, 王旭, 任强, 段文霞, 汪杭军1, 李振基, 雷加强, 杨慧敏, 李建章, 何茜, 宋永明, 余新晓, 刘勇, 方陆明, 黎明, 尹光彩, 周亮, 杨娅, 刘勇, 朱波, 周国逸, 高岚, 韩士杰, 崔同林, 阎秀峰, 周宇飞, JIANGXi_dian, 苏晓华, 玲, 王清文, 孙向阳, 徐扬, 王春林, 虞木奎, 唐小明, 张冰玉, 沈熙环, 周国逸, 徐新文, 李吉跃, 鹿振友, 赖志华, 刘锐, 柯水发, 喻理飞, 宗文君, 程云清, HEXiu_bin, 王新杰, 国庆, 3, 周晓梅, 宋爱琴, 李丙文, 王伟宏, 李俊清, 李晓兰, , 郭蓓, 茹广欣, 齐涛, 孙阁, 陈峻崎, 陈培金, 陈实, 温亚利, 张志毅, 王建林, 蒋德明, 王旭, 张可栋, 王晓静, 周玉平, 姚永刚, 唐旭利, 长山, 刘志明, 宋湛谦, 王春林, 关少华, 陈放, 赵双荣, 杨伟伟, 闫俊华, 郑凌峰.  小波和神经网络在色木孔洞缺陷超声定量检测中的应用 . 北京林业大学学报, 2007, 29(2): 128-132.
  • 加载中
图(19) / 表 (2)
计量
  • 文章访问数:  1286
  • HTML全文浏览量:  118
  • PDF下载量:  28
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-10-23
  • 修回日期:  2017-11-23
  • 刊出日期:  2018-01-01

基于速度误差校正的林木应力波无损检测断层成像算法

doi: 10.13332/j.1000-1522.20170384
    基金项目:

    江苏省重点研发计划项目 BE2016627

    中央高校基本科研业务费专项 RP51635B

    无锡市国际科技研发合作项目 CZE02H1706

    江苏省普通高校专业学位研究生实践创新计划项目 SJLX16_0499

    国家自然科学基金项目 61472368

    作者简介:

    焦治。主要研究方向:智能无损检测技术。Email:novrule@163.com 地址: 214122 江苏省无锡市江南大学物联网工程学院

    通讯作者: 李光辉,博士,教授。主要研究方向:传感器网络、容错计算和智能无损检测技术。Email: ghli@jiangnan.edu.cn 地址:同上
  • 中图分类号: S781;TS67

摘要: 目的应力波断层成像技术已经在林木无损检测领域得到了广泛的应用, 然而木材的各向异性特征以及应力波速度反演计算的误差对断层成像的精度影响较大。因此,进一步提高断层图像的准确性非常关键。方法根据传统的应力波速度反演原理, 提出了一种基于速度误差校正的断层成像算法(ECIA)。该算法利用最小二乘QR分解法(LSQR)计算应力波在林木横截面网格单元内的速度分布, 并使用误差校正机制(ECM)优化断层图像。为了评估算法的性能, 分别选取了若干实验室内的原木试样及扬州市区古树进行无损检测实验, 利用德国PICUS应力波断层成像仪获取的应力波传播数据实现了ECIA算法, 生成了各实验样本的断层图像。结果ECIA算法较为精确地检测出了原木及活树中缺陷区域的位置及大小, 尤其在活树的健康检测中, ECIA算法对于缺陷区域尤其是轻微腐朽区域的检测精度高于PICUS检测仪。结论ECIA算法能够产生较为准确的林木断层图像,适用于林木应力波无损检测。

English Abstract

焦治, 李光辉, 武夕. 基于速度误差校正的林木应力波无损检测断层成像算法[J]. 北京林业大学学报, 2018, 40(1): 108-119. doi: 10.13332/j.1000-1522.20170384
引用本文: 焦治, 李光辉, 武夕. 基于速度误差校正的林木应力波无损检测断层成像算法[J]. 北京林业大学学报, 2018, 40(1): 108-119. doi: 10.13332/j.1000-1522.20170384
Jiao Zhi, Li Guang-hui, Wu Xi. Tomography imaging algorithm based on velocity error correction for stress wave nondestructive evaluation of wood[J]. Journal of Beijing Forestry University, 2018, 40(1): 108-119. doi: 10.13332/j.1000-1522.20170384
Citation: Jiao Zhi, Li Guang-hui, Wu Xi. Tomography imaging algorithm based on velocity error correction for stress wave nondestructive evaluation of wood[J]. Journal of Beijing Forestry University, 2018, 40(1): 108-119. doi: 10.13332/j.1000-1522.20170384
  • 无损检测对于发现林木内部缺陷,防止腐朽扩散,提高木材利用率具有重要意义。应力波无损检测技术具有成本低,操作方便,无辐射,不受被测木材的大小形状影响等优点,已经在林木无损检测中得到了广泛的应用。应力波无损检测技术通过测量应力波在树木横截面中多个传播方向的传播时间,然后反演计算应力波在内部单元的传播速度,并生成木材的断层图像,从而准确、直观地检测到木材内部缺陷大小与位置[1]。国内外关于木材应力波无损检测技术的研究已有大量成果。Dikrallah等[2]通过实验分析了湿材的声学各向异性,研究了应力波速度与方向角之间的数学关系。梁善庆等[3]与美国农业部林产品实验室合作,将弹性波断层成像技术运用到活树的检测,结果表明该技术能够直观显示活树内部腐朽位置、大小、形状和严重程度。杨学春等[4]对林木内部的应力波传播规律进行了研究,建立了相关的波动方程及其数值求解模型。Li等[5]研究了健康树木横截面内部的应力波传播速度变化规律,建立了数学模型,利用美国黑樱桃树(Prunus serotina)的无损检测实验结果验证了该模型的有效性。Socco等[6]提出了一种较为复杂的数据处理方法用于对活树的健康监测,并在地震波图像算法的基础上研究了一种反演算法用来生成树干的应力波传播时间图像,该算法对一些小的缺陷(如小空洞)不敏感,而且没有很好地考虑木材的各向异性问题。Du等[7]在考虑木材各向异性的基础上,提出了一种改进型的椭圆空间插值方法,并利用速度补偿获得更精确的空间插值数据,该方法能够有效抵抗由于缺陷区域密度变化引起的信号干扰。目前,国际市场上已有不少林木应力波断层成像检测设备(例如,匈牙利产的FAKOPP,德国产的PiCUS Sonic Tomograph和Arbotom),均能够生成木材横截面的二维或三维断层图像,但成像效果和检测精度仍然有待改善[8]

    为了进一步提高林木应力波断层图像的准确性,提出了一种应力波速度反演的误差校正机制(ECM)。在此基础上,利用最小二乘QR分解法(LSQR)设计了一种新的林木应力波无损检测断层成像算法(ECIA),并且利用大量的真实样本实验验证了该算法的有效性。

    • 介绍应力波在健康树木中的传播速度模型及无损检测基本思想,并描述林木无损检测断层成像的速度反演过程。

    • 为简单起见,我们将树干视为一个圆柱体,其横截面为正圆形,如图 1所示,其中S代表应力波发射端传感器的位置,R1R2代表应力波接收端传感器的位置。在检测过程中,通过脉冲锤敲击发射端,产生应力波,并从发射端传播到接收端。检测仪测量发射端与接收端传感器之间的应力波的传播时间。SR2代表应力波沿直径的传播路径,其传播速度称为径向传播速度,用VR表示。SR1代表应力波沿与直径成β角的弦的传播路径,其传播速度称为弦向传播速度,用VT表示[9]。通常在健康树木中,应力波的径向传播速度最快,并且随着弦向角β的增大,传播速度会变慢。Li等[5]通过研究建立了应力波在健康树木中的速度模型:

      $$ {V_{\rm{T}}}/{V_{\rm{R}}} \approx 1 - 0.2{\beta ^2} $$ (1)

      图  1  树木横截面内应力波传播示意图

      Figure 1.  Stress wave propagation in the cross-section of tree trunk

      式中:VT/VR表示弦向传播速度和径向传播速度的比值,β表示弦向角。

      树木产生腐朽或空洞往往从内部开始,树干外层材质相对而言更健康,所以实际检测过程中,经常使用各相邻传感器之间的应力波传播速度平均值作为健康参考值,并根据公式(1)计算横截面内任意方向的参考速度Vref[10]。对任意的传播路径,将实际测量得到的传播速度Vtest和参考传播速度Vref进行比较,可以初步判断该路径是否存在缺陷。根据国内外学者的研究经验[1],若VtestVref的比值小于0.9,可以认为该路径存在缺陷。

    • 速度反演是指根据检测仪测量得到的不同传播路径上的应力波平均速度,反演计算应力波在横截面内不同单元的速度分布的过程[11]。目前大多数的木材断层成像反演算法都是基于地震波传播的模拟算法[12]

      图 2所示,反演算法需将树干横截面划分为若干网格单元,且每个网格单元内需要至少有一条路径穿过。从SR为一条应力波传播路径,该路径的传播时间t可以用慢度p(速度的倒数)从起点S到终点R的线积分表示:

      $$ t = \int_S^R p {\rm{d}}l $$ (2)

      图  2  树木横截面网格单元划分示意图

      Figure 2.  Grid cells of a cross-section of a tree trunk

      式中:dl表示路径长度的增量,p是应力波的传播慢度(速度的倒数)。

      为了计算慢度分布,把图 2所示的横截面划分为M个网格单元,从左上角第一个网格单元开始,自左上至右下依次编号,将连续函数p(x, y)近似成每个网格单元内离散的慢度分布。假定每个网格单元内的慢度不变,将式(2)离散表示为[13]

      $$ {t_i} = \sum\limits_{j = 1}^M {{p_j}} {d_{ij}}\quad (i = 1, 2, \cdots , N) $$ (3)

      式中:dij表示第i条路径被第j个网格单元分割的长度,N表示总的应力波传播路径数。式(3)可用矩阵的形式表示:

      $$ \mathit{\boldsymbol{T}}{\rm{ = }}\mathit{\boldsymbol{DP}} $$ (4)

      式中:TP分别代表N×1的传播时间列向量和M×1的传播慢度列向量,D代表传播路径被网格单元分割的长度距离矩阵,其大小为N×M

      对于给定的TD,只需要计算出慢度列向量P后就可以生成初步的断层图像。由于一条传播路径只会被少数几个网格单元分割,所以对于任意路径i在大多数网格单元中,被分割的路径长度dij的值是0,因此式(4)中的系数矩阵D通常是一个大型稀疏矩阵。LSQR方法是一种利用Lanczos迭代法求解最小二乘问题的方法,由于其计算量小,运用矩阵的稀疏性简化计算,因此十分适合求解大型稀疏矩阵问题[14]

      由于应力波检测仪的传感器数量有限,应力波传播路径总数N通常小于网格单元总数M,故矩阵D的行数小于列数,式(4)为欠定方程组[15]。欠定方程组有多个解,为了求欠定方程组的特解,LSQR需要先将系数矩阵D对角化[14, 16]

      $$ \left\{ \begin{array}{l} \begin{array}{*{20}{l}} {\mathit{\boldsymbol{D}} = \mathit{\boldsymbol{UB}}{\mathit{\boldsymbol{V}}^{\rm{T}}}}\\ {{\mathit{\boldsymbol{U}}_k} = \left[ {{u_1}, {u_2}, \cdots , {u_k}} \right]}\\ {{\mathit{\boldsymbol{V}}_k} = \left[ {{v_1}, {v_2}, \cdots , {v_k}} \right]} \end{array}\\ {\mathit{\boldsymbol{B}}_k} = \left[ {\begin{array}{*{20}{c}} {{\alpha _1}}&{}&{}&{}\\ {{\beta _2}}&{{\alpha _2}}&{}&{}\\ {}&{{\beta _3}}& \ddots &{}\\ {}&{}& \ddots &{{\alpha _k}}\\ {}&{}&{}&{{\beta _{k + 1}}} \end{array}} \right] \end{array} \right. $$ (5)

      通过逐列对比,得到迭代关系式,利用此迭代方法实现矩阵D的下双对角分解。

      $$ \left\{ {\begin{array}{*{20}{l}} {{\mathit{\boldsymbol{U}}_{k + 1}}\left( {{\beta _1}{e_1}} \right) = b}\\ {\mathit{\boldsymbol{D}}{\mathit{\boldsymbol{V}}_k} = {\mathit{\boldsymbol{U}}_{k + 1}}{\mathit{\boldsymbol{B}}_k}}\\ {{\mathit{\boldsymbol{D}}^{\rm{T}}}{\mathit{\boldsymbol{U}}_{k + 1}} = {\mathit{\boldsymbol{V}}_k}\mathit{\boldsymbol{B}}_k^{\rm{T}} + {\alpha _{k + 1}}{v_{k + 1}}\mathit{\boldsymbol{e}}_{k + 1}^{\rm{T}}} \end{array}} \right. $$ (6)

      式中:ek+1T表示n阶单位矩阵的第k+1行。同理,矩阵D可以分解为上双对角矩阵。系数矩阵对角化完成后,设:

      $$ \left\{ {\begin{array}{*{20}{l}} {{r_k} = b - \mathit{\boldsymbol{D}}{x_k}}\\ {{x_k} = {\mathit{\boldsymbol{V}}_k}{y_k}}\\ {{t_{k + 1}} = {\beta _1}{e_1} - {\mathit{\boldsymbol{B}}_k}{y_k}} \end{array}} \right. $$ (7)

      联立式(6)和式(7)可得:

      $$ \begin{array}{*{20}{c}} {{r_k} = b - \mathit{\boldsymbol{D}}{x_k} = {\mathit{\boldsymbol{U}}_{k + 1}}\left( {{\beta _1}{e_1}} \right) - \mathit{\boldsymbol{D}}{\mathit{\boldsymbol{V}}_k}{y_k} = }\\ {{\mathit{\boldsymbol{U}}_{k + 1}}\left( {{\beta _1}{e_1}} \right) - {\mathit{\boldsymbol{U}}_{k + 1}}{\mathit{\boldsymbol{B}}_k}{y_k} = {\mathit{\boldsymbol{U}}_{k + 1}}{t_{k + 1}}} \end{array} $$ (8)

      由于Uk+1是正交矩阵,所以通过QR方法解yk使‖tk+1‖取极小值,求解最小二乘问题min‖β1e1- Bkyk‖,通过反复迭代,最终会收敛于极小值,从而得到式(4)的最终解。

    • 经过速度反演得到每个网格单元的速度分布后,树木内部的断层成像就可以初步生成了,但是VtestVref的比值的阈值0.9仅是一个经验值,可能会造成实际情况和断层图像的偏差。提出一种速度误差校正机制,用来提高成像算法的精度。

      以往的断层成像算法中[13],计算每个网格单元的传播速度Vtest和参考速度Vref后,通过计算两者的比值Vtest/Vref直接决定该网格单元是否存在缺陷区域:若比值小于0.9,则该网格单元标记为缺陷区域,否则标记为健康区域。这种评判标准有时会因为速度反演的不准确导致成像误差。事实上,树木内部的缺陷通常会从一个区域向其邻近的区域扩散,所以邻居网格单元之间的健康状况是互相关联的。提出了一种基于邻居网格单元的误差校正机制ECM,对反演得到的δ分布进行校正,用来减少上述的潜在误差。ECM分析由LSQR方法求出的每个网格单元的速度分布及其相应的VtestVref比值δ分布(式9),对当前网格单元的δ值进行调整,从而减少误差。

      $$ {\delta _i} = {V_{{\rm{ test }}}}/{V_{{\rm{ ref }}}} $$ (9)

      式中:δi代表第i个网格单元的VtestVref比值,全部网格单元的δ值构成δ分布(用向量R表示)。

      由于ECM是基于邻居网格单元的δ值进行校正的,在横截面边缘的网格单元会产生位于横截面外部的邻居网格单元(如图 3中编号为1、2、M-1、和M的网格单元),称这种网格单元为虚拟单元。因此,不同于速度反演过程,在误差校正过程中,需要将包括虚拟单元在内的所有网格单元重新编号。由于虚拟单元内没有路径穿过,故将虚拟单元内的传播速度及参考速度均设置为0,其δ值为无穷大,但为了防止计算时产生数据异常,算法将其δ值校正为0。

      图  3  树干横截面内网格单元的分类

      Figure 3.  Grid cell classification of a cross-section of a tree trunk

      完全处于横截面内部的网格单元以及与横截面边缘相交的网格单元统称为横截面单元,ECM对两种横截面单元应用了不同的校正策略。如图 3所示,对于完全处于横截面内部的第i个网格单元,只需要考虑其4个紧邻的网格单元(mni1i2)的δ值大小。对于第j个网格单元,其4个紧邻的网格单元中有虚拟单元l,此时需要考察包括其对角邻居单元在内的所有8个邻居单元(k1kk2j1j2l1,l,l2)。这8个邻居单元中,第l1ll2个网格单元为虚拟单元,第k1kk2j1j2个网格单元为横截面单元,故将上述5个网格单元作为第j个网格单元的新邻居网格单元。

      ECM是根据当前网格单元及其邻居网格单元的δ值,对其δ值进行校正,校正方案如式(10):

      $$ {\delta _i} = \left\{ \begin{array}{l} 1, \;\;\;\;s > l/2\\ 1, \;\;\;\;s = l/2且{\delta _i} > 0.9\\ 0, \;\;\;s = l/2且{\delta _i} \le 0.9\\ 0, \;\;\;\;s < l/2 \end{array} \right. $$ (10)

      式中:δi代表第i个网格单元的δ值,s代表其邻居网格单元中δ>0.9的个数,l代表邻居网格单元的数量。该校正机制从第一个网格单元开始,逐个扫描所有网络单元,并根据公式(10)进行校正。这个过程反复进行,直到所有网格单元的δ值保持稳定为止。ECM算法的误差校正过程如下所示:

      算法名:ECM误差校正机制

      输入:反演生成的δ分布(用向量R表示);

      输出:误差校正后的δ分布(用向量R′表示);

      开始:

      1:do

      2:for i=1,2,…,M do

      3:  if第i个网格单元的四邻域均为横截面单元

      4:   then l=4;

      5:  else

      6:  l=第i个网格单元的八邻域中横截面单元的个数;

      7:  end if

      8:  根据式(10)对R进行误差校正,得到R′;

      9:end for

      10:if R′R

      11:   then flag=1;

      12:else

      13:   flag=0;

      14: end if

      15: R=R′;

      16:while(flag==1)

      结束

    • 基于速度反演和误差校正机制设计了断层成像算法(ECIA)。算法将树干横截面模拟成一个正圆,通过检测仪器测量应力波在树木内部的传播速度,利用LSQR方法反演计算应力波在树干横截面各个网格单元内的速度分布,并使用ECM校正速度误差,最终生成树木横截面的断层图像。

      ECIA算法步骤如下:

      输入.传感器个数、树木横截面周长、应力波传播时间、传感器两两间距离。

      输出.树木横截面断层图像

      Step1.根据输入数据建立坐标系,模拟树木横截面,划分网格单元,计算每个网格单元内的参考速度。

      Step2.计算应力波传播路径被网格单元分割的距离,得到距离矩阵D

      Step3.利用LSQR反演求解速度分布,并求相应的传播速度与参考速度的比值δ分布向量R

      Step4.调用ECM算法对δ分布向量R进行误差校正。

      Step5.根据各网格单元的δ值进行颜色填充。

      Step6.利用区域生长算法分割缺陷区域与健康区域,对缺陷区域边缘进行曲线拟合以优化图像。

    • 分别采用实验室原木和活树完成了两组无损检测成像实验。原木试样包括含空洞的柳杉(Cryptomeria fortunei,含水率11.5%)、含空洞的乌桕(Sapium sebiferum,含水率12.0%)和健康的柳杉(含水率12.5%)作为原木组实验对象,并分别记为试样1、2、3。此外,选取了江苏扬州市区6棵有内部缺陷的古树(香樟(Cinnamomum camphora)、枫香(Liquidambar formosana)、桑树(Morus alba)、圆柏(Sabina chinensis)、垂柳(Salix babylonica)、圆柏)作为实验对象,分别编号为试样4, 5, …, 9。实地检测的场景如图 4所示。实验过程中使用德国生产的PICUS应力波断层成像仪采集应力波传播数据,使用国产泰克曼TM410木材水分仪测量原木试样的含水率。除此之外,实验中还使用了Resistograph微钻阻力测试仪测量树木内部的阻力值,作为评估活树内部缺陷大小及位置的参考标准,从而验证成像算法的准确性。

      图  4  活树应力波无损检测现场

      Figure 4.  Stress wave testing of a live standing tree

    • 利用PICUS对3个原木试样进行无损检测,采集应力波在试样内部的传播速度,并生成断层图像,3个原木试样所测截面周长分别为101、99和100 cm。基于相同的应力波传播速度数据,利用前文设计的ECIA算法产生断层图像,并与原木试样实际缺陷以及PICUS的断层成像进行对比。本组实验中,对于两种断层图像和实际的原木横截面照片,利用AutoCAD中不规则图形面积的计算方法得到原木横截面面积及缺陷面积[17],从而计算缺陷区域的面积与整个横截面面积之比(RS),比较断层图像与试样实际情况的RS,从而定量评估算法的准确性。

      图 5~7分别展示了试样1、2、3的实际情况及两种断层成像算法所生成的图像。其中,ECIA生成的断层图像中,仅用棕色区域表示缺陷区域,PICUS生成的断层图像中,绿色、紫色和蓝色都代表不同程度的缺陷区域。因此,计算PICUS的缺陷面积时,需要计算这几部分的面积之和。

      图  5  试样1的断层图像比较

      Figure 5.  Comparison of tomography images of sample 1

      图  6  试样2的断层图像比较

      Figure 6.  Comparison of tomography images of sample 2

      图  7  试样3的断层图像比较

      Figure 7.  Comparison of tomography images of sample 3

      定量比较结果如表 1所示,其中,ΔRS表示断层图像的RS与真实RS的相对误差,计算方法为用它俩的差值除以真实的RS

      表 1  原木断层成像结果的定量评价

      Table 1.  Quantitative evaluation to tomography images of logs

      试样编号
      No. of
      sample
      试样横截面积
      Cross-sectional area
      of the sample/cm2
      真实S
      Actual S/cm2
      真实RS
      Actual RS/%
      ECIAPICUS
      S/cm2RS/%ΔRS/%S/cm2RS/%ΔRS/%
      1812.593 5297.284 136.584 6324.558 839.941 1-9.174 6232.514 728.613 921.787 6
      2794.226 078.019 29.823 3102.108 912.856 4-30.876 6285.921 436.000 0-266.475 6
      3804.247 70000088.467 211.000 0>0
      注:S代表试样或断层图像中的缺陷面积;RS代表试样或断层图像中缺陷面积与横截面积的比值;ΔRS代表断层图像的RS与真实RS的相对误差。Notes:S is the defective area, RS is the ratio of the defective area to the cross-section area, and ΔRS is the relative error between the RS of ECIA or PICUS and the actual RS.

      表 1中,S代表试样或断层图像中的缺陷面积,ΔRS≠0代表断层图像与试样的真实情况存在误差。从表 1可知,ECIA的断层图像比PICUS更接近于真实情况,但是其断层图像仍存在误检区域(如图 6)。较大的误检区域是由原木的低含水量造成的。国内外大量的研究与实验表明:应力波技术更适用于活树的健康检测[6]

    • 与原木实验类似,使用PICUS对6棵古树进行无损检测,采集应力波传播速度数据。由于无法砍伐树木获得真实的横截面图像,我们在同一高度的横截面上使用Resistograph进行微钻阻力测试。对于选定的应力波传播路径,Resistograph检测仪能够准确地测量该路径上缺陷区域的长度,经过多路径测量,可以确定缺陷区域的大体位置。图 8~19分别给出了6棵古树的实物图、ECIA和PICUS的断层图像,以及在同一横截面中若干路径的阻力曲线。断层图像和阻力曲线中SijEij的部分代表试样在第i条阻力曲线中第j段缺陷区域,SE分别代表该段缺陷区域的起点和终点。阻力曲线图中横坐标代表路径长度,蓝色部分为去掉树皮部分的真实长度,纵坐标代表相对阻力。在每条阻力曲线中,路径m-n表示与断层图像中第m个传感器到第n个传感器的路径相对应。

      图  8  试样4的断层图像比较

      Figure 8.  Comparison of tomography images of sample 4

      图  9  试样4的阻力曲线图

      Figure 9.  Resistance curves of sample 4

      图  10  试样5的断层图像比较

      Figure 10.  Comparison of tomography images of sample 5

      图  11  试样5的阻力曲线图

      Figure 11.  Resistance curves of sample 5

      图  12  试样6的断层图像比较

      Figure 12.  Comparison of tomography images of sample 6

      图  13  试样6的阻力曲线图

      Figure 13.  Resistance curves of sample 6

      图  14  试样7的断层图像比较

      Figure 14.  Comparison of tomography images of sample 7

      图  15  试样7的阻力曲线图

      Figure 15.  Resistance curves of sample 7

      图  16  试样8的断层图像比较

      Figure 16.  Comparison of tomography images of sample 8

      图  17  试样8的阻力曲线图

      Figure 17.  Resistance curves of sample 8

      图  18  试样9的断层图像比较

      Figure 18.  Comparison of tomography images of sample 9

      图  19  试样9的阻力曲线图

      Figure 19.  Resistance curves of sample 9

      从试样4的外观可以观察出在1号和2号传感器之间有明显的缺陷区域,PICUS的检测结果也是如此。但通过Resistograph的阻力检测,发现1号和2号传感器之间仅有树皮缺陷,缺陷区域向树木内部扩散不深,反而在树木内部其他位置(如图 8S12-E12S22-E22S31-E31),有若干较小的缺陷区域,ECIA通过速度反演算法检测出了这些缺陷区域,并通过ECM保证了对轻微缺陷区域的检测精度。

      图 10a所示,试样5有明显的空洞,由于这个空洞过大,使得应力波只能沿接近树皮的边缘部分传播,造成了应力波数据异常,导致PICUS生成的断层图像中缺陷区域比实际的空洞小,如路径3-8中S11-E11部分在阻力曲线图中为缺陷区域,而在图 10c中不全为缺陷。ECIA使用ECM校正了由于数据异常造成的误差,使得断层图像更精确。

      由于试样6的树干胸径较大,Resistograph无法穿透树木,所以只能取得路径上的部分阻力值数据,图 12b12c中从路径起点到L点为取得阻力值数据的部分。如图 12a,试样6外部有一处明显的缺陷区域,并且被养护人员填充了特殊的材料,这些填料造成了Resistograph测量的阻力值较高。另外,如图 13aS11-E11部分阻力值相对较低,可视为缺陷区域。ECIA生成的断层图像准确地检测出了S11-E11段的缺陷区域,但PICUS生成的断层图像中S11-E11部分仅存在长度较小的轻微腐朽,与阻力曲线图相差较大。

      图 14图 15分别给出了试样7的断层图像比较及3条应力波路径的阻力曲线。从图 14可以看出,ECIA和PICUS产生的断层图像对缺陷区域的位置及形状判断大致相同,但ECIA生成的断层图像中缺陷区域较大。路径1-7中两种算法都有不同程度增大了缺陷区域的长度,但ECIA生成的断层图像中路径1-7中的误检长度相对较小。路径2-9中,ECIA成功的检测到了S23-E23部分的缺陷而PICUS在该部分只检测到了较小长度的缺陷。同样在路径4-10中,PICUS生成的断层图像中S31-E31部分的缺陷长度比阻力曲线图 15中给出的该部分缺陷长度小,而ECIA生成的断层图像在该部分的缺陷长度与阻力曲线图中的缺陷长度较为符合。

      对于试样8,两种算法生成的断层图像(如图 16bc所示)十分相近,相比阻力曲线反映的实际情况,断层图像的路径3-9中,ECIA和PICUS对缺陷区域的位置长度判断都较为准确,在路径5-9中由于树木节子的存在(如图 17a,60~70 mm处阻力值突然增大的部分),导致应力波传播数据异常,使得两种算法生的断层图像在该条路径上的缺陷区域的长度都增大了。

      图 18a,试样9的外部存在一个空洞,虽然空洞较小但一直扩散至树干内部,两种成像算法都成功地检测出了空洞。但在路径10-6中,S12-E12部分存在较小的缺陷区域,两种算法生成的断层图像中该部分都标记为健康区域。ECIA生成的断层图像在路径10-5中的缺陷区域与阻力曲线图 19基本一致,达到了较高的精确度,但PICUS无论在路径10-5中还是路径7-12中都存在一定程度的漏检。

      活树无法像原木一样测量缺陷区域的面积,为了评估成像算法的精确度,计算了部分路径上缺陷区域的长度与路径长度的比(RL),利用Resistograph测量真实的缺陷长度,断层图像中对应路径上的缺陷长度可以通过几何计算得到。如同原木实验中的ΔRS,定义了断层图像的RL与真实RL的相对误差ΔRL。对于6棵活树试样,ECIA和PICUS生成的断层图像的实验结果如表 2所示。实验结果证明ECIA能够更准确地检测出轻微腐朽的区域。对于试样4、5、7、9,ECIA成像结果的相对误差绝对值小于8%。由于试样6的树干存在较大空洞,试样8树干中存在节子,这些特殊情况造成检测到的应力波数据异常,使得两种算法的成像结果与实际情况相比都存在较大的误差,但ECIA的相对误差值仍小于PICUS。总体来说,ECIA的精确度比PICUS高,能够较为准确地反映实际情况。但是ECIA产生的断层图像中缺陷区域的边缘真实感略显不足,需要更完善的校正机制和算法来改进断层图像的真实感。

      表 2  活树断层成像结果的定量评价

      Table 2.  Quantitative evaluation to tomography images of standing trees

      试样编号
      Sample No.
      路径
      Path
      真实RL
      Actual RL
      ECIAPICUS
      RLΔRLRLΔRL
      42-1017.910 417.254 13.664 412.537 329.999 9
      2-913.432 812.455 27.277 70100.000 0
      5-15.970 16.275 8-5.120 518.209 0-205.003 2
      53-888.549 687.042 31.702 267.605 623.652 3
      4-783.333 377.500 07.000 053.928 635.285 7
      5-1074.626 876.417 9-2.400 168.656 77.999 9
      69-175.555 593.043 5-23.145 943.478 342.455 1
      2-1184.444 488.888 9-5.263 288.888 9-5.263 2
      71-789.937 193.548 4-4.015 495.913 9-6.645 5
      2-975.496 781.073 5-7.386 856.989 224.514 3
      4-1091.463 493.566 7-2.299 680.444 412.047 4
      85-93.125 013.437 5-330.000 014.375 0-360.000 0
      3-93.703 74.521 1-22.069 812.592 6-240.000 0
      910-632.558 131.395 33.571 527.665 615.027 0
      10-527.907 027.702 30.733 536.511 6-30.833 1
      7-1246.774 248.387 1-3.448 343.037 77.988 4
      注:RL代表试样或断层图像中路径上缺陷区域的长度与整个路径长度的比值;ΔRL代表断层图像的RL与真实RL的相对误差。Notes: RL is the ratio of the defective length to the whole length of the path under test, and ΔRL is the relative error between the RL of ECIA or PICUS and the actual RL.
    • 提出了一种适用于林木应力波无损检测的断层成像算法ECIA,算法运用LSQR方法求解速度反演方程,反演计算出横截面内每个网格单元的应力波传播速度。为了减少速度反演产生的误差,提出了一种基于邻居网格单元的误差校正机制ECM,使用ECM减少了速度反演产生的误差,从而提高了断层图像的准确度。通过对9个样本的两组实验,证实了本算法的准确度高于PICUS的准确度。本算法已在扬州市区古树名木健康监测项目中得到了实际应用。今后,我们将会完善ECM及成像算法,进一步提高断层图像的真实感。

参考文献 (17)

目录

    /

    返回文章
    返回