留言板

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

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

木材径切面内部缺陷的应力波成像算法

郑泽宇 冯海林 杜晓晨 方益明

陈诗燕, 程鸿浩, 吴筱萌, 等. 基于Fuzzy分析的茶园游猎蛛取食小贯小绿叶蝉的种间竞争作用[J]. 浙江农林大学学报, 2023, 40(5): 1008-1017. DOI: 10.11833/j.issn.2095-0756.20220753
引用本文: 郑泽宇, 冯海林, 杜晓晨, 等. 木材径切面内部缺陷的应力波成像算法[J]. 浙江农林大学学报, 2019, 36(2): 211-218. DOI: 10.11833/j.issn.2095-0756.2019.02.001
CHEN Shiyan, CHENG Honghao, WU Xiaomeng, et al. Interspecific competition of wandering spiders feeding on Empoasca onukii in tea plantations based on Fuzzy analysis[J]. Journal of Zhejiang A&F University, 2023, 40(5): 1008-1017. DOI: 10.11833/j.issn.2095-0756.20220753
Citation: ZHENG Zeyu, FENG Hailin, DU Xiaochen, et al. A stress wave tomography algorithm for internal defects in radial and longitudinal (RL) planes of wood[J]. Journal of Zhejiang A&F University, 2019, 36(2): 211-218. DOI: 10.11833/j.issn.2095-0756.2019.02.001

木材径切面内部缺陷的应力波成像算法

DOI: 10.11833/j.issn.2095-0756.2019.02.001
基金项目: 

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

浙江省重大科技专项资助项目 2015C03008

浙江省重大科技专项资助项目 2018C02013

详细信息
    作者简介: 郑泽宇, 从事木材无损检测研究。E-mail:404625184@qq.com
    通信作者: 冯海林, 教授, 博士, 从事智能信息处理、大数据、物联网等研究。E-mail:hlfeng@zafu.edu.cn
  • 中图分类号: S781.3

A stress wave tomography algorithm for internal defects in radial and longitudinal (RL) planes of wood

  • 摘要: 为了获得木材径切面上的缺陷形状、大小和位置,提出一种木材径切面内部缺陷成像的方法。首先,基于应力波在木材径切面上的传播规律,提出一种木材径切面上的应力波速度修正方法。将应力波速度转换为径切面上的若干个预估点的值,结合反距离加权插值(IDW)法提出一种速度修正插值(VCI)方法。最后,使用VCI方法在不同树木样本上进行了木材径切面缺陷的二维成像实验。结果表明:①VCI方法可以重建木材内部缺陷大小以及缺陷位置,缺陷成像结果与真实的缺陷情况相吻合。②对比IDW方法的成像结果,VCI方法对缺陷的大小以及缺陷形状、位置的成像结果有较大提高。③利用混淆矩阵方法对VCI与IDW方法进行定量分析表明,VCI方法的平均准确率、平均精确度和平均查全率均高于IDW方法,说明VCI方法成像效果的可行性和有效性。
  • 小贯小绿叶蝉 Empoasca onukii在中国茶园广泛分布,是对茶树Camellia sinensis危害最大的害虫之一[1]。小贯小绿叶蝉对茶树嫩茎、嫩叶的刺吸会导致茶叶焦枯,抑制茶树正常生长,严重影响茶叶的产量和品质[23]。游猎性蜘蛛(游猎蛛)作为一种捕食性天敌,捕食量大、行动活跃、捕食范围广,对防治害虫、调节生态系统平衡都有重要作用[46]。利用天敌与害虫之间的捕食关系对茶园害虫进行生物防治的方法已有诸多报道[716],这些研究大都直接讨论天敌与害虫之间的关联,对于不同种天敌之间的竞争作用却极少研究。但单一地增加茶园害虫的优势种天敌等势必会导致整个茶园生态系统的不稳定[1718]。FAUST等[19]将2个物种之间的作用关系分为中性、竞争、偏害、寄生、捕食、偏利和互利共生等几种类型。天敌和害虫之间的关系分为捕食和寄生2种类型,天敌之间为了争夺资源多数是竞争关系[20]。除了天敌与害虫之间的捕食和被捕食关系外,各天敌之间的竞争关系也需要得到重视。因此,从生物多样性等多维度对茶园整个生态环境进行调控更有利于长期高效防治茶园害虫,以达到绿色环保的平衡调节作用[21]

    在昆虫种群数量生态学的研究中,对害虫及天敌数据的统计分析通常为抽样分析,采用的分级方式为连续等分分级,这会使得不同精度级宽下的结果具有主观性和偶然性。张书平等[22]将Fuzzy分级法与灰色关联度法相结合研究茶园与假眼小绿叶蝉 Empoasca vitis数量上关系密切的天敌种类,是一种较好的种间关系研究方法。多数学者对林木种间竞争的研究较多[23],未见对茶园各天敌种间竞争关系的报道,将数据按不同方式分级后再研究种间竞争关系的研究更少。为此,本研究运用灰色关联度法,并将Fuzzy分级法与竞争系数法相结合,研究‘农抗早’‘Nongkangzao’和‘平阳特早’‘Pingyangtezao’茶园游猎蛛取食小贯小绿叶蝉的种间竞争关系,并通过竞争强度指数验证所得结论,为茶园合理保护和利用种间竞争关系强的天敌,有效防治小贯小绿叶蝉提供科学依据。

    在安徽农业大学科技示范园(31°56′N,117°12′E)中共调查了2种茶园,分别为‘农抗早’和‘平阳特早’茶园,面积均为0.2 hm2。调查时间为2021年3月25日至11月19日,隔15 d调查1次,共14次。茶园互不相连,按常规措施管理。在春茶采摘结束和秋末进行茶树修剪,并且在秋末进行茶园耕翻,加强秋冬季管理,及时除草、修剪茶树和采摘茶叶,诱杀和人工捕杀害虫,不施用化学农药。

    采用平行跳跃法在茶园随机选取3行,每行茶树相隔1 m,在每行间隔1 m处设置1个2 m×1 m的样方,每行取10个样方,每个茶园取30个样方。用目测调查的方式在每个样方随机选取10片叶,调查一些不易振落的害虫及天敌的种类和数量,再用盘拍法对样方中所有枝条进行盘拍(盘拍所用搪瓷盘口长为40 cm,宽30 cm,搪瓷盘上喷洒稀释1 000倍的洗衣粉水溶液)。调查记载盘中害虫及天敌物种数和个体数,不能准确鉴定的物种样本编号保存,装瓶带回室内由专家鉴定。

    1.3.1   灰色关联度法

    将害虫数量(Yi)及其主要天敌数量(Xj)分别看作1个本征系统,害虫数量作为该系统的参照序列,天敌作为比较序列,并把不同样方的害虫及其主要天敌的数量作为该序列在第k个样方的效果白化值,进行双序列关系分析。数据经均值化后得:

    $$ Y_i=\left\{Y_i(1),\; Y_i(2),\; \cdots ,\; Y_i(n)\right\},\; i=1; $$ (1)
    $$ X_j=\left\{X_j(1),\; X_j(2),\; \cdots,\; X_j(n)\right\}, \;j=1,\;2,\; \cdots,\; M。 $$ (2)

    式(1)和式(2)中:n表示样方数,M表示天敌种类数。YiXj在第k个样方上的关联系数rij为:

    $$ {r}_{ij}=\frac{{\mathrm{min}}_{i}{\mathrm{min}}_{j}\left|{Y}_{i}\left(k\right)-{X}_{j}\left(k\right)\right|+\rho {\mathrm{max}}_{i}{\mathrm{max}}_{j}\left|{Y}_{i}\left(k\right)-{X}_{j}\left(k\right)\right|}{\left|{Y}_{i}\left(k\right)-{X}_{j}\left(k\right)\right|+\rho {\mathrm{max}}_{i}{\mathrm{max}}_{j}\left|{Y}_{i}\left(k\right)-{X}_{j}\left(k\right)\right|} 。 $$ (3)

    式(3)中:ρ 为分辨系数,取值区间为[0, 1],一般取 ρ = 0.5,为扩大各物种之间关联度的差异,本研究取 ρ = 0.8。$|Y_{i}\left(k\right)-{X}_{j}\left(k\right) |$ 为序列YiXj在第k点上差的绝对值;$ \mathrm{m}\mathrm{i}\mathrm{n}\left|{Y}_{i}\left(k\right)-{X}_{j}\left(k\right)\right| $为1级最小差,$ {\mathrm{min}}_{i}{\mathrm{min}}_{j}\left|{Y}_{i}\left(k\right)-{X}_{j}\left(k\right)\right| $ 为2级最小差。$ \mathrm{m}\mathrm{a}\mathrm{x}\left|{Y}_{i}\left(k\right)-{X}_{j}\left(k\right)\right| $与$ {\mathrm{max}}_{i}{\mathrm{max}}_{j}\left|{Y}_{i}\left(k\right)-{X}_{j}\left(k\right)\right| $ 分别为1级和2级最大差。利用该公式可求出第j种天敌(Xj)与害虫(Yi)数量间的关联度为$R({Y}_{i},{X}_{j})=1/n\displaystyle \sum _{k=1}^{M}{r}_{ij}\left(k\right)$,其大小反映害虫与天敌相互联系的紧密程度。天敌与害虫数量间关联度越大,表明天敌与害虫数量关系越密切[24]。本研究由数据处理系统软件(DPS系统)进行灰色关联度数据的运算。

    1.3.2   模糊分级法

    在论域[A1, A2, A3$, \;\cdots , $ An$, \;\cdots , $ Am]上按所求解问题的性质和要求规定的一个隶属函数μi,叫作Fuzzy分级隶属函数。将通常连续等分分级的频数作为原始数据,设原始数据为[a1, a2, a3$,\; \cdots , $ an$, \;\cdots , $ am],称${\hat{{a}}}_{{i}}=\displaystyle \sum _{{n}={i}-5}^{{i}+5}{{\mu }}_{{i}}{{a}}_{{n}}$为第i个Fuzzy等级的Fuzzy频数。它在论域上的分布曲线叫Fuzzy频数曲线[2526]。本研究规定Fuzzy分级隶属函数μi为:

    $$ {\mu }_{i}=\left\{\begin{array}{l}1.0\;n=i\\ 0.8\;n=i+1,\;i-1\\ 0.6\;n=i+2,\;i-2\\ 0.4\;n=i+3,\;i-3\\ 0.2\;n=i+5,\;i-5\\ 0.1\;n=i+5,\;i-5\\ 0\;n > i+5,\;n < i-5\end{array}\right. 。 $$ (4)

    式(4)中:n为原始数据项数,i为Fuzzy频数项数。则Fuzzy频数为$ {\hat{a}}_{i}=\displaystyle \sum _{n=i-5}^{i+5}{\mu }_{i}{a}_{n}={a}_{i}+ 0.8\left({a}_{i+1}+{a}_{i-1}\right)+0.6\left({a}_{i+2}+{a}_{i-2}\right)+0.4\left({a}_{i+3}+{a}_{i-3}\right)+0.2\left({a}_{i+4}+{a}_{i-4}\right)+0.1\left({a}_{i+5}+{a}_{i-5}\right) $。为了解茶园游猎蛛的自然种群动态,减少在对游猎蛛进行竞争关系分析时由种群数据导致的误差,本研究根据游猎蛛种群数量变化幅度,对游猎蛛的种群数据按照30个样方中5只游猎蛛的级宽进行等分分级统计(级宽不宜太宽也不宜太窄),算出各级出现的频数,再以30个样方中5只游猎蛛为级宽的频数作为原始数据即[a1, a2, a3$,\; \cdots , $ an$, \;\cdots , $ am]进行Fuzzy分级统计。由Excel 2019计算Fuzzy频数。

    1.3.3   竞争系数法

    一般用Levins的生态重叠公式计算竞争系数[27]

    $$ {\alpha }_{ij}=\sum _{k=1}^{n}\left({P}_{ik}{P}_{jk}\right)/\sum _{k=1}^{n}{P}_{ik}^{2} 。 $$ (5)

    式(5)中:$ {P}_{ik} $ 和$ {P}_{jk} $分别为种i和种j在第k个样方中的相对优势度,n为样地数。

    使用重叠的方法计算竞争系数的公式[28]

    $$ {\alpha }_{ij}=\sum _{k=1}^{n}\left({P}_{ik}{P}_{jk}\right)/\sqrt{\sum _{k=1}^{n}{P}_{ik}^{2}\sum _{k=1}^{n}{P}_{jk}^{2}} 。 $$ (6)

    式(6)中:αij为种j对种i的竞争系数,PikPjk分别是第i个物种和第j个物种使用的第k个资源的比例。αij = αji,PIANKA[28]评价该式是对称的,并称其为重叠值。MAY[2930]将该式与其他表达式进行比较,沿用了PIANKA的表达式,并将其作为按函数比例进行计算的竞争系数。本研究中所求得的竞争系数均按上述(6)式计算。

    1.3.4   竞争强度法

    本研究将关联度与竞争系数相结合,引入竞争强度指数的概念。不同天敌与同一害虫之间关联度的比值称为相对密切度,该相对密切度与2种天敌之间竞争系数的乘积即为竞争强度指数,种i对种j的竞争强度指数(C)为:

    $$ C={\alpha }_{ij}\left(\frac{{R}_{iy}}{{R}_{jy}}\right) 。 $$ (7)

    式(7)引入了害虫与天敌关系因素,RiyRjy分别为天敌XiXj与害虫Y数量上的关联度,αij为种j对种i的竞争系数。本研究中竞争系数和竞争强度指数均由 Excel 2019 对Fuzzy频数计算所得,再使用DPS数据处理系统软件用多重比较方法中的 Duncan 新复极差法进行竞争系数和竞争强度指数平均值和显著水平的分析[31]

    选取调查日期中2个茶园游猎蛛数量最多的7种蜘蛛作为主要天敌,分别为鞍型花蟹蛛Xysticus ephippiatus、三突花蟹蛛Ebrechtella tricuspidata、粽管巢蛛Clubiona japonicola、斑管巢蛛C. reichlini、斜纹猫蛛Oxyopes sertatus、黑色跳蛛Plexippus paykulli和条纹蝇虎P. setipes。它们的数量动态见表1

    表 1  2个茶园小贯小绿叶蝉与游猎蛛数量动态
    Table 1  Population dynamics of E. onukii and wandering spiders in two tea plantations
    茶树品种日期
    (月-日)
    害虫数量/头游猎蛛数量/头
    X1X2X3X4X5X6X7
    ‘农抗早’ 03-2531016323277
    04-161127107019
    05-07422600220
    05-233015005121
    06-042227150150
    06-2017520555050
    07-08594146551886380
    08-1020927231560202624
    08-224292319295933256
    09-0692199192357186
    09-173017171310311319
    10-13181204545
    11-0152429939419
    11-194711361535
    合计 1 937310250155177253155262
    ‘平阳特早’03-259230013174
    04-16110500162
    05-07311600060
    05-23573410050
    06-0425115030110
    06-201391581820101
    07-08674264624342270
    08-10318173174233130
    08-225348182148311029
    09-0616029121418391610
    09-17201312121937129
    10-134720020112
    11-01235283223
    11-1910812352624
    合计 2 339245230128128170118174
      说明:X1~X7分别指鞍型花蟹蛛、三突花蟹蛛、粽管巢蛛、斑管巢蛛、斜纹猫蛛、黑色跳蛛和条纹蝇虎数量(头)。
    下载: 导出CSV 
    | 显示表格

    用灰色关联度法求得小贯小绿叶蝉与游猎蛛数量之间的灰色关联度(表2)。由表2可知:与小贯小绿叶蝉在数量上关联度最大的前3位天敌,‘农抗早’茶园为斑管巢蛛、粽管巢蛛和黑色跳蛛;‘平阳特早’茶园为斑管巢蛛、粽管巢蛛和三突花蟹蛛。2个茶园均有斑管巢蛛和粽管巢蛛。

    表 2  2个茶园小贯小绿叶蝉与游猎蛛间的灰色关联度
    Table 2  Grey correlation between E. onukii and wandering spiders in two tea plantations
    游猎蛛‘农抗早’‘平阳特早’游猎蛛‘农抗早’‘平阳特早’
    灰色关联度排位灰色关联度排位灰色关联度排位灰色关联度排位
    鞍型花蟹蛛0.833 850.826 14斜纹猫蛛0.754 860.794 76
    三突花蟹蛛0.856 940.856 23黑色跳蛛0.857 030.802 65
    粽管巢蛛 0.890 220.857 12条纹蝇虎0.744 270.763 47
    斑管巢蛛 0.896 710.860 31
    下载: 导出CSV 
    | 显示表格

    图1图2可看出:2个茶园Fuzzy分级频数集中性明显大于等分分级频数。图1A中在第4、7和12个级宽处出现3个明显的峰值,而图1B中 7条曲线均由高到低一致趋于平缓。图2A中从第11个级宽处多数曲线已贴近横轴,而图2B中的曲线从第15个级宽开始才趋于横轴。通过Fuzzy分级法得到的新数据明显增加了原始数据的有效区间,并且Fuzzy频数由原始数据通过Fuzzy隶属函数求得,每个原始数据都与几个等级发生联系,而等分分级里每个原始数据只与1个等级发生联系,因此Fuzzy频数无明显的分级界限。对2种频数进行t检验,‘农抗早’茶园中鞍型花蟹蛛、三突花蟹蛛、粽管巢蛛、斑管巢蛛、斜纹猫蛛、黑色跳蛛和条纹蝇虎这7种蜘蛛的t值依次为2.333、2.379、2.281、2.646、3.937、3.059、1.975,‘平阳特早’茶园这7种蜘蛛的$ t $值依次为 2.176、2.048、2.153、2.477、3.151、2.090、1.591。自由度为34时,t0.01=2.728,t0.05=2.032,t0.10=1.091,除条纹蝇虎外,2个茶园各游猎蛛2种频数间的差异都显著,并且斜纹猫蛛的2种频数间差异极显著。Fuzzy频数的直接图示是曲线,可保持原始数据的精度,故用此法分级的数据进行后续竞争关系的分析,其结果更接近实际。

    图 1  ‘农抗早’茶园游猎蛛种群数量等分分级频数曲线(A)和Fuzzy频数曲线(B)
    Figure 1  Equally graded frequency curve (A) and Fuzzy frequency curve (B) of the number of wandering spiders in ‘Nongkangzao’ tea plantation        
    图 2  ‘平阳特早’茶园游猎蛛种群数量等分分级频数曲线(A)和Fuzzy频数曲线(B)
    Figure 2  Equally graded frequency curve (A) and Fuzzy frequency curve (B) of the number of wandering spiders in ‘Pingyangtezao’ tea plantation        

    以Fuzzy频数作为原始数据,计算2个茶园各游猎蛛之间的竞争系数并将结果列于表3。对各游猎蛛及其竞争对手之间的竞争系数进行方差分析,再用新复极差法分析天敌之间竞争系数的差异性(表4)。所得结果中,当竞争对手分别为鞍型花蟹蛛、三突花蟹蛛、粽管巢蛛和条纹蝇虎时,游猎蛛间差异显著,此时不同竞争对手下竞争力最强的游猎蛛不同,故暂时无法得出茶园竞争力最强的蜘蛛类别,但除斑管巢蛛和黑色跳蛛外,无论竞争对手为哪种游猎蛛,所有游猎蛛都与斜纹猫蛛差异显著,且斜纹猫蛛的均值都最小。因此,可得出结论:除斑管巢蛛和黑色跳蛛外,5种游猎蛛里斜纹猫蛛竞争力最弱。

    表 3  2个茶园各游猎蛛之间的竞争系数
    Table 3  Competition coefficients among wandering spiders in two tea plantations
    茶树品种竞争对手游猎蛛间竞争系数
    鞍型花蟹蛛三突花蟹蛛粽管巢蛛斑管巢蛛斜纹猫蛛黑色跳蛛条纹蝇虎
    ‘农抗早’ 鞍型花蟹蛛1.000 00.991 50.983 70.939 80.906 30.959 20.991 0
    三突花蟹蛛0.991 51.000 00.995 20.973 30.893 80.928 20.995 3
    粽管巢蛛 0.983 70.995 21.000 00.980 30.901 50.916 70.993 6
    斑管巢蛛 0.939 80.973 30.980 31.000 00.872 90.856 80.961 4
    斜纹猫蛛 0.906 30.893 80.901 50.872 91.000 00.959 80.871 7
    黑色跳蛛 0.959 20.928 20.916 70.856 80.959 81.000 00.915 9
    条纹蝇虎 0.991 00.995 30.993 60.961 40.871 70.915 91.000 0
    ‘平阳特早’鞍型花蟹蛛1.000 00.989 30.998 50.976 50.913 70.993 50.975 1
    三突花蟹蛛0.989 31.000 00.984 20.977 40.869 70.984 40.992 5
    粽管巢蛛 0.998 50.984 21.000 00.967 40.909 30.996 30.967 4
    斑管巢蛛 0.976 50.977 40.967 41.000 00.939 70.954 00.957 1
    斜纹猫蛛 0.913 70.869 70.909 30.939 71.000 00.880 60.826 5
    黑色跳蛛 0.993 50.984 40.996 30.954 00.880 61.000 00.967 7
    条纹蝇虎 0.975 10.992 50.967 40.957 10.826 50.967 71.000 0
    下载: 导出CSV 
    | 显示表格
    表 4  2个茶园各游猎蛛之间竞争系数新复极差法分析结果
    Table 4  Results of the new multiple range test analysis of competition coefficients among various wandering spiders in two tea plantations
    竞争对手F物种2个茶园竞争
    系数均值
    5%显著
    水平
    1%极显著
    水平
    竞争对手F物种2个茶园竞争
    系数均值
    5%显著
    水平
    1%极显著
    水平
    鞍型花蟹蛛7.530粽管巢蛛 0.991 1aA斜纹猫蛛1.203黑色跳蛛 0.920 2aA
    三突花蟹蛛0.990 4aA鞍型花蟹蛛0.910 0aA
    条纹蝇虎 0.983 0aA斑管巢蛛 0.906 3aA
    黑色跳蛛 0.976 4aA粽管巢蛛 0.905 4aA
    斑管巢蛛 0.958 2aAB三突花蟹蛛0.881 7aA
    斜纹猫蛛 0.910 0bB条纹蝇虎 0.849 1aA
    三突花蟹蛛11.356条纹蝇虎 0.993 9aA黑色跳蛛0.562鞍型花蟹蛛0.976 4aA
    鞍型花蟹蛛0.990 4aA粽管巢蛛 0.956 5aA
    粽管巢蛛 0.989 7aA三突花蟹蛛0.956 3aA
    斑管巢蛛 0.975 4aA条纹蝇虎 0.941 8aA
    黑色跳蛛 0.956 3aA斜纹猫蛛 0.920 2aA
    斜纹猫蛛 0.881 7bB斑管巢蛛 0.905 4aA
    粽管巢蛛 3.305鞍型花蟹蛛0.991 1aA条纹蝇虎12.037三突花蟹蛛0.993 9aA
    三突花蟹蛛0.989 7aA鞍型花蟹蛛0.983 0aA
    条纹蝇虎 0.980 5aA粽管巢蛛 0.980 5aA
    斑管巢蛛 0.973 9aA斑管巢蛛 0.959 2aA
    黑色跳蛛 0.956 5abA黑色跳蛛 0.941 8aA
    斜纹猫蛛 0.905 4bA斜纹猫蛛 0.849 1bB
    斑管巢蛛 1.610三突花蟹蛛0.975 3aA
    粽管巢蛛 0.973 9aA
    条纹蝇虎 0.959 2aA
    鞍型花蟹蛛0.958 2aA
    斜纹猫蛛 0.906 3aA
    黑色跳蛛 0.905 4aA
      说明:计算不同竞争对手下各游猎蛛在2个茶园竞争系数的平均值,5%水平上均数最大的标记为a,1%水平上均数最大的标记为A。向下比较,与之差异性不显著的标记相同字母,差异性显著的标记不同字母。
    下载: 导出CSV 
    | 显示表格

    结合关联度和竞争系数将所求的竞争强度指数列于表5。对竞争强度指数进行方差分析,用新复极差法进行比较(表6)。由表6可知:无论竞争对手为哪种游猎蛛,粽管巢蛛都与其他蜘蛛差异显著,其次为斑管巢蛛,且游猎蛛中斜纹猫蛛与两者差异极显著。因此,可得出结论:粽管巢蛛竞争力最强,斑管巢蛛次之,斜纹猫蛛竞争力最弱。

    表 5  2个茶园各游猎蛛之间的竞争强度指数
    Table 5  Competition intensity indices among wandering spiders in two tea plantations
    茶树品种竞争对手游猎蛛间竞争强度指数
    鞍型花蟹蛛三突花蟹蛛粽管巢蛛斑管巢蛛斜纹猫蛛黑色跳蛛条纹蝇虎
    ‘农抗早’ 鞍型花蟹蛛1.000 01.019 01.050 31.010 70.820 40.985 90.884 5
    三突花蟹蛛0.964 71.000 01.033 81.018 50.787 30.928 30.864 3
    粽管巢蛛 0.921 40.958 01.000 00.987 50.764 40.882 50.830 6
    斑管巢蛛 0.873 90.930 10.973 21.000 00.734 80.818 90.797 8
    斜纹猫蛛 1.001 21.014 71.063 21.037 01.000 01.089 80.859 4
    黑色跳蛛 0.933 20.928 10.952 20.896 50.845 31.000 00.795 3
    条纹蝇虎 1.110 41.146 11.188 61.158 50.884 21.054 81.000 0
    ‘平阳特早’鞍型花蟹蛛1.000 01.025 41.036 11.017 00.879 00.965 30.9012
    三突花蟹蛛0.954 41.000 00.985 20.982 00.807 20.922 70.884 9
    粽管巢蛛 0.962 30.983 21.000 00.971 00.843 10.932 90.861 7
    斑管巢蛛 0.937 60.972 80.963 81.000 00.868 10.890 00.849 3
    斜纹猫蛛 0.949 70.937 00.980 71.017 21.000 00.889 30.794 0
    黑色跳蛛 1.022 61.050 21.064 01.022 60.872 01.000 00.920 5
    条纹蝇虎 1.055 11.113 11.086 11.078 50.860 41.017 31.000 0
    下载: 导出CSV 
    | 显示表格
    表 6  2个茶园各游猎蛛之间的竞争强度指数新复极差法分析结果
    Table 6  Results of the new multiple range test analysis of competitive intensity indices among various wandering spiders in two tea plantations
    竞争对手F物种2个茶园竞争
    强度指数均值
    5%显著
    水平
    1%极显著
    水平
    竞争对手F物种2个茶园竞争
    强度指数均值
    5%显著
    水平
    1%极显著
    水平
    鞍型花蟹蛛32.992粽管巢蛛 1.043 2aA斜纹猫蛛2.142斑管巢蛛 1.027 1aA
    三突花蟹蛛1.022 2abA粽管巢蛛 1.022 0aA
    斑管巢蛛 1.013 8abA黑色跳蛛 0.989 5abA
    黑色跳蛛 0.975 6bA三突花蟹蛛0.975 8abA
    条纹蝇虎 0.892 9cB鞍型花蟹蛛0.975 4abA
    斜纹猫蛛 0.849 7cB条纹蝇虎 0.826 7bA
    三突花蟹蛛33.492粽管巢蛛 1.009 5aA黑色跳蛛1.575粽管巢蛛 1.008 1aA
    斑管巢蛛 1.000 2aAB三突花蟹蛛0.989 1aA
    鞍型花蟹蛛0.960 9abAB鞍型花蟹蛛0.977 9aA
    黑色跳蛛 0.925 5bBC斑管巢蛛 0.959 6aA
    条纹蝇虎 0.874 6cC斜纹猫蛛 0.858 6aA
    斜纹猫蛛 0.797 3dD条纹蝇虎 0.857 9aB
    粽管巢蛛 9.706斑管巢蛛 0.979 3aA条纹蝇虎10.533粽管巢蛛 1.137 3aA
    三突花蟹蛛0.970 6aAB三突花蟹蛛1.129 6aA
    鞍型花蟹蛛0.941 9aAB斑管巢蛛 1.118 5aA
    黑色跳蛛 0.907 7abABC鞍型花蟹蛛1.082 7aA
    条纹蝇虎 0.846 2bcBC黑色跳蛛 1.036 1aA
    斜纹猫蛛 0.803 8cC斜纹猫蛛 0.872 3bB
    斑管巢蛛 3.588粽管巢蛛 0.968 5aA
    三突花蟹蛛0.951 4abA
    鞍型花蟹蛛0.905 8abcA
    黑色跳蛛 0.854 4abcA
    条纹蝇虎 0.823 6bcA
    斜纹猫蛛 0.801 5cA
      说明:计算不同竞争对手下各游猎蛛在2个茶园竞争强度指数的平均值,5%水平上均数最大的标记为a,1%水平上均数最大的标记为A。向下比较,与之差异性不显著的标记相同字母,差异性显著的标记不同字母。
    下载: 导出CSV 
    | 显示表格

    本研究首先通过灰色关联度分析初步得出2个茶园与小贯小绿叶蝉数量相关性最大的游猎蛛均为斑管巢蛛和粽管巢蛛,再对7种游猎蛛的种群数量进行Fuzzy分级统计,将得出的Fuzzy频数作为原始数据进行竞争系数分析,结果显示:除斑管巢蛛和黑色跳蛛外,5种游猎蛛中斜纹猫蛛竞争力最弱。为验证结果准确性,综合灰色关联度和竞争系数结果引入了竞争强度指数概念,得出7种游猎蛛中斜纹猫蛛竞争力最弱,并且在任何竞争对手下粽管巢蛛都与其他蜘蛛差异显著,其次为斑管巢蛛,即在取食茶园小贯小绿叶蝉时粽管巢蛛和斑管巢蛛竞争力最强,斜纹猫蛛竞争力最弱。

    在进行Fuzzy分级统计时,并未对小贯小绿叶蝉数量进行分级,因为在进行竞争关系分析时,只需要7种游猎蛛的数量数据,故无需对小贯小绿叶蝉数量做同样的处理。在进行灰色关联度分析时,小贯小绿叶蝉和7种游猎蛛数量均未进行处理,原因是小贯小绿叶蝉数量与7种游猎蛛数量数据大小相差较大,选择过大的级宽会导致游猎蛛数量均处于第1级宽内,选择较小级宽会出现较多级层且在靠后多个级层里只有小贯小绿叶蝉数据,而游猎蛛数据均为0。对7种游猎蛛的种群数量进行Fuzzy分级,然后进行竞争关系的统计计算,使它们的数据集中性更加突出,弥补了抽样时造成的误差,是一种简洁有效的计算方法。

    2个茶园竞争力最强和最弱的蜘蛛相同。本研究的2个茶园均按常规措施管理,不使用化学农药,且于冬季除草修剪,修剪会影响天敌的虫口基数,而‘农抗早’和‘平阳特早’抗逆性和抗寒性强[32],因此在相近的受害程度后,‘农抗早’和‘平阳特早’恢复时间均较短,恢复效果均较好,害虫所处环境变化的速度一致性可能是研究结果相同的原因之一。此外,茶园竞争力最强的是斑管巢蛛,这可能与斑管巢蛛的生活习性有关。斑管巢蛛定居且游猎于树冠上被害卷叶或枯叶等阴暗干燥处,白天基本不出行,黄昏时刻,蜘蛛开始活动,主动巡游猎取食,沿着枝、叶逐一搜索前进,几乎无遗漏之处。合理保护和利用斑管巢蛛这类竞争力强的蜘蛛可达到有效防治小贯小绿叶蝉的目的。

    对茶园游猎蛛之间竞争作用的研究可以更好地理解游猎蛛之间的竞争如何影响害虫的数量和变化趋势以及如何影响生物防治的有效性[33]。至今利用害虫与天敌的种间关系对茶园害虫进行生物防治更多地还是停留在增加优势种天敌数量的方向上,分析天敌之间的竞争作用有利于选出最高效的天敌组合,在不破坏茶园原有生态环境的基础上高效防治害虫。

    本研究中的7种游猎蛛都是广食性天敌,但为了研究方便,把它们作为只取食小贯小绿叶蝉一种食饵的单食性天敌,若进行深入研究就需要考虑多种猎物(害虫)共存时天敌对食物的嗜食性,在此基础上应用竞争关系分析方法就能更加真实地反映天敌之间的关系。另外,本研究是分析天敌两两之间的竞争关系,实际上,在食饵不足时,7种天敌之间也存在竞争关系,这种情况有待进一步研究。

  • 图  1  木材径切面模型

    Figure  1  Radial and longitudinal plane model of wood

    图  2  应力波在不同缺陷程度的木材中传播

    Figure  2  Stress waves propagate in wood of different defect degree

    图  3  应力波测量装置与传感器的布置

    Figure  3  Stress wave measuring device and sensor arrangement

    图  4  雪松样本实验结果

    Figure  4  Sample experiment results of Cedrus deodara sample

    图  5  泡桐样本实验结果

    Figure  5  Sample experiment results of Paulownia fortunei sample

    图  6  樟树样本实验结果

    Figure  6  Sample experiment results of Cinnamomum camphora sample

    表  1  实验样本基本信息

    Table  1.   Experimental sample information

    样本编号样本种类样本周长/cm检测区域面积/cm2缺陷面积/cm2
    1雪松Cedrus deodara62.831 000111.20
    2泡桐Paulownia fortunei73.831 175427.50
    3樟树Cinnamomum camphora60.95940314.15
    下载: 导出CSV

    表  2  缺陷与预测的4种组合

    Table  2.   Four combination of real defects and predict defects

    类别预测缺陷预测完好
    真实缺陷被正确地预测为缺陷的面积C被错误地预测为完好的面积E
    真实完好被错误地预测为缺陷的面积F被正确地预测为完好的面积W
    下载: 导出CSV

    表  3  不同样本的准确率

    Table  3.   Accuracy of different samples

    样本编号准确率A/%
    IDW方法VCI方法
    189.2095.85
    286.3693.30
    372.3290.86
    下载: 导出CSV

    表  4  不同样本的精确度

    Table  4.   Precision of different samples

    样本编号精确度P/%
    IDW方法VCI方法
    167.6076.92
    277.1888.35
    365.7981.51
    下载: 导出CSV

    表  5  不同样本的查全率

    Table  5.   Recall rate of different samples

    样本编号查全率R/%
    IDW方法VCI方法
    117.9889.92
    256.1494.03
    335.8094.00
    下载: 导出CSV
  • [1] WANG Xiping. Acoustic measurements on trees and logs:a review and analysis[J]. Wood Sci Technol, 2013, 47(5):965-975.
    [2] 徐华东, 王立海, 游祥飞, 等.应力波在旱柳立木内的传播规律分析及其安全评价[J].林业科学, 2010, 46(8):145-150.

    XU Huadong, WANG Lihai, YOU Xiangfei, et al. Analysis of stress wave propagation in Hankow Willow standing trees and stability assessment[J]. Sci Silv Sin, 2010, 46(8):145-150.
    [3] 杨学春, 王立海.红松木材结构缺陷对应力波传播参数的影响[J].东北林业大学学报, 2005, 33(1):30-31.

    YANG Xuechun, WANG Lihai. Impacts of structure defaults on propagation parameters of stress wave in wood of Korean pine[J]. J Northeast For Univ, 2005, 33(1):30-31.
    [4] LIN C J, CHANG T T, JUAN M Y, et al. Stress wave tomography for the quantification of artificial hole detection in camphor trees (Cinnamomum camphora)[J]. Taiwan J For Sci, 2011, 26(1):17-32.
    [5] FENG Hailin, LI Guanghui, FU Sheng, et al. Tomographic image reconstruction using an interpolation method for tree decay detection[J]. Bioresources, 2014, 9(2):3248-3263.
    [6] DU Xiaochen, LI Saozhi, LI Guanghui, et al. Stress wave tomography of wood internal defects using ellipse-based spatial interpolation and velocity compensation[J]. Bioresources, 2015, 10(3):3948-3962.
    [7] 翁翔, 李光辉, 冯海林, 等.应力波在树木径切面内的传播速度模型[J].林业科学, 2016, 52(7):104-112.

    WENG Xiang, LI Guanghui, FENG Hailin, et al. Stress wave propagation velocity model in rl plane of standing trees[J]. Sci Silv Sin, 2016, 52(7):104-112.
    [8] 徐华东, 王立海.杨木径切面准纵波传播路径追踪及速度变化试验研究[J].声学学报, 2012(1):62-67.

    XU Huadong, WANG Lihai. Experiment study on the propagation path tracing and velocity change of quasi-P wave in the radial section of poplar wood[J]. Acta Acust, 2012(1):62-67.
    [9] LI Guanghui, WANG Xiping, FENG Hailin, et al. Analysis of wave velocity patterns in black cherry trees and its effect on internal decay detection[J]. Comput Electron Agric, 2014, 104(2):32-39.
    [10] 冯海林, 李光辉, 方益明, 等.应力波传播模型及其在木材检测中的应用[J].系统仿真学报, 2010, 22(6):1490-1493.

    FENG Hailin, LI Guanghui, FANG Yiming, et al. Stress wave propagation modeling and application in wood testing[J]. J Syst Simulation, 2010, 22(6):1490-1493.
    [11] LI Li, WANG Xiping, WANG Lihai, et al. Acoustic tomography in relation to 2D ultrasonic velocity and hardness mappings[J]. Wood Sci Technol, 2012, 46(1/3):551-561.
    [12] HU Yingcheng, NAKAO T, NAKAI T, et al. Dynamic properties of three types of wood-based composites[J]. J Wood Sci, 2005, 51(1):7-12.
    [13] BRANCHERIAU L, ASHKAN G, GALLET P, et al. Application of ultrasonic tomography to characterize the mechanical state of standing trees (Picea abies)[J]. J Phys Conf Ser, 2012, 253(1):1-13.
    [14] 贾悦, 崔宁博, 魏新平, 等.基于反距离权重法的长江流域参考作物蒸散量算法适用性评价[J].农业工程学报, 2016, 32(6):130-138.

    JIA Yue, CUI Ningbo, WEI Xinping, et al. Applicability evaluation of different algorithms for reference crop evapotranspiration in Yangtze River Basin based on inverse distance weighted method[J]. Trans Chin Soc Agric Eng, 2016, 32(6):130-138.
    [15] SAMMUT C, WEBB G I. Encyclopedia of Machine Learning[M]. New York:Springer, 2010.
    [16] FAWCETT T. An introduction to ROC analysis[J]. Pattern Recogn Lett, 2006, 27(8):861-874.
    [17] POWERS D M W. Evaluation:from precision, recall and F-Factor to ROC, informedness, markedness & correlation[J]. J Mach Learn Technol, 2011, 2(1):37-63.
  • [1] 杨凡, 杨博凯, 李荣荣.  基于图像分割和深度学习的人造板表面缺陷检测 . 浙江农林大学学报, 2024, 41(1): 176-182. doi: 10.11833/j.issn.2095-0756.20230280
    [2] 邢东, 胡建鹏, 姚利宏.  热处理木材的材性预测与质量控制研究现状与发展 . 浙江农林大学学报, 2020, 37(4): 793-800. doi: 10.11833/j.issn.2095-0756.20190449
    [3] 周凡, 周永东, 高鑫, 付宗营, 侯俊峰, 翁翔, 林韶辉.  黑木相思木材干燥特性及干燥工艺制定 . 浙江农林大学学报, 2020, 37(3): 571-577. doi: 10.11833/j.issn.2095-0756.20190394
    [4] 陈勇平, 郭文静, 唐启恒.  基于雷达检测的木材内部孔洞面积测算与修正 . 浙江农林大学学报, 2020, 37(6): 1193-1199. doi: 10.11833/j.issn.2095-0756.20190709
    [5] 汪紫阳, 尹世逵, 李颖, 李耀翔.  基于可见/近红外光谱识别东北地区常见木材 . 浙江农林大学学报, 2019, 36(1): 162-169. doi: 10.11833/j.issn.2095-0756.2019.01.020
    [6] 杨建飞, 宁莉萍, 杨了, 王杰, 钱钰滢, 陈颐萱.  黑壳楠木材构造特征及挥发性有机物成分 . 浙江农林大学学报, 2018, 35(5): 927-934. doi: 10.11833/j.issn.2095-0756.2018.05.018
    [7] 李任, 黄荣凤, 常建民, 高志强, 伍艳梅.  预热温度对层状压缩木材力学性能的影响 . 浙江农林大学学报, 2018, 35(5): 935-941. doi: 10.11833/j.issn.2095-0756.2018.05.019
    [8] 张春晓, 冯海林, 李光辉, 王燕凤, 杜晓晨.  应力波在树木不同角度纵截面的传播速度模型 . 浙江农林大学学报, 2017, 34(5): 926-933. doi: 10.11833/j.issn.2095-0756.2017.05.021
    [9] 刘光林, 李光辉, 孙晔, 方益明.  树木内部应力波传播速度模型 . 浙江农林大学学报, 2015, 32(1): 18-24. doi: 10.11833/j.issn.2095-0756.2015.01.003
    [10] 蒋明衔, 陈奶荣, 林巧佳, 曾钦志, 饶久平.  响应面法优化杨木防腐胶合板的制作工艺 . 浙江农林大学学报, 2014, 31(1): 122-128. doi: 10.11833/j.issn.2095-0756.2014.01.019
    [11] 张高品, 李光辉, 李剑, 冯海林.  木材无损检测技术中的应力波传播时延估计算法 . 浙江农林大学学报, 2014, 31(3): 394-398. doi: 10.11833/j.issn.2095-0756.2014.03.010
    [12] 陈利芳, 马红霞, 张燕君, 王剑菁, 曹永建.  浸注药剂热处理工艺对木材物理性能的影响 . 浙江农林大学学报, 2013, 30(6): 955-959. doi: 10.11833/j.issn.2095-0756.2013.06.023
    [13] 伏建国, 刘金良, 杨晓军, 安榆林, 骆嘉言.  分子生物学技术应用于木材识别的研究进展 . 浙江农林大学学报, 2013, 30(3): 438-443. doi: 10.11833/j.issn.2095-0756.2013.03.022
    [14] 单海斌, 刘杏娥, 杨淑敏, 孙建平, 杨喜.  计算机断层扫描技术(CT)在木材无损检测中的应用与发展 . 浙江农林大学学报, 2013, 30(1): 123-128. doi: 10.11833/j.issn.2095-0756.2013.01.018
    [15] 杨小军, 孙友富.  碳纤维木材复合材结合层剪切蠕变性能研究 . 浙江农林大学学报, 2012, 29(2): 192-196. doi: 10.11833/j.issn.2095-0756.2012.02.007
    [16] 尹建新, 祁亨年, 冯海林, 杜晓晨.  一种基于混合纹理特征的木板材表面缺陷检测方法 . 浙江农林大学学报, 2011, 28(6): 937-942. doi: 10.11833/j.issn.2095-0756.2011.06.017
    [17] 钱俊, 金永明, 沈哲红, 俞友明, 楼永生.  电场强度与作用时间对马尾松木材含水率梯度的影响 . 浙江农林大学学报, 2005, 22(2): 193-197.
    [18] 金春德, 张美淑, 文桂峰, 刘继生, 张鹏.  不同坡位人工林赤松木材材性的径向变异 . 浙江农林大学学报, 2004, 21(2): 119-124.
    [19] 邵千钧, 徐群芳, 王伟龙.  木材干燥过程控制策略与方法的研究 . 浙江农林大学学报, 2003, 20(3): 307-310.
    [20] 刘志坤, 李延军, 杜春贵, 文桂峰, 林勇.  刨切薄竹生产工艺研究 . 浙江农林大学学报, 2003, 20(3): 227-231.
  • 期刊类型引用(0)

    其他类型引用(1)

  • 加载中
  • 链接本文:

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

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2019/2/211

图(6) / 表(5)
计量
  • 文章访问数:  3538
  • HTML全文浏览量:  795
  • PDF下载量:  181
  • 被引次数: 1
出版历程
  • 收稿日期:  2018-03-26
  • 修回日期:  2018-07-13
  • 刊出日期:  2019-04-20

木材径切面内部缺陷的应力波成像算法

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

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

    浙江省重大科技专项资助项目 2015C03008

    浙江省重大科技专项资助项目 2018C02013

    作者简介:

    郑泽宇, 从事木材无损检测研究。E-mail:404625184@qq.com

    通信作者: 冯海林, 教授, 博士, 从事智能信息处理、大数据、物联网等研究。E-mail:hlfeng@zafu.edu.cn
  • 中图分类号: S781.3

摘要: 为了获得木材径切面上的缺陷形状、大小和位置,提出一种木材径切面内部缺陷成像的方法。首先,基于应力波在木材径切面上的传播规律,提出一种木材径切面上的应力波速度修正方法。将应力波速度转换为径切面上的若干个预估点的值,结合反距离加权插值(IDW)法提出一种速度修正插值(VCI)方法。最后,使用VCI方法在不同树木样本上进行了木材径切面缺陷的二维成像实验。结果表明:①VCI方法可以重建木材内部缺陷大小以及缺陷位置,缺陷成像结果与真实的缺陷情况相吻合。②对比IDW方法的成像结果,VCI方法对缺陷的大小以及缺陷形状、位置的成像结果有较大提高。③利用混淆矩阵方法对VCI与IDW方法进行定量分析表明,VCI方法的平均准确率、平均精确度和平均查全率均高于IDW方法,说明VCI方法成像效果的可行性和有效性。

English Abstract

陈诗燕, 程鸿浩, 吴筱萌, 等. 基于Fuzzy分析的茶园游猎蛛取食小贯小绿叶蝉的种间竞争作用[J]. 浙江农林大学学报, 2023, 40(5): 1008-1017. DOI: 10.11833/j.issn.2095-0756.20220753
引用本文: 郑泽宇, 冯海林, 杜晓晨, 等. 木材径切面内部缺陷的应力波成像算法[J]. 浙江农林大学学报, 2019, 36(2): 211-218. DOI: 10.11833/j.issn.2095-0756.2019.02.001
CHEN Shiyan, CHENG Honghao, WU Xiaomeng, et al. Interspecific competition of wandering spiders feeding on Empoasca onukii in tea plantations based on Fuzzy analysis[J]. Journal of Zhejiang A&F University, 2023, 40(5): 1008-1017. DOI: 10.11833/j.issn.2095-0756.20220753
Citation: ZHENG Zeyu, FENG Hailin, DU Xiaochen, et al. A stress wave tomography algorithm for internal defects in radial and longitudinal (RL) planes of wood[J]. Journal of Zhejiang A&F University, 2019, 36(2): 211-218. DOI: 10.11833/j.issn.2095-0756.2019.02.001
  • 木材的内部缺陷会严重影响木材的质量。由于应力波的木材无损检测成本低、携带方便及对木材无损的特性,使得其越来越受到人们的青睐。国内外已开展了应力波速度与木材内部缺陷之间的相关研究。如WANG[1]通过应力波在缺陷木材中的传播规律,表明应力波在腐烂或劣化的木材中比在完好的木材中传播速度慢。徐华东等[2]以校园内的活立木作为样本,使用应力波检测仪器对校园树木进行二维成像,表明在活立木的断层横截面上,应力波横向的波速要小于径向的波速,且应力波在木材内部的波速受制于木材内部的缺陷情况。国内学者也对应力波无损检测仪器的可行性进行了研究。杨学春等[3]利用应力波检测仪器对原木内部腐朽进行了研究,结果表明:应力波测试仪能准确判断不同树种原木内部的腐朽程度,并且能够得到原木内部腐朽基本形状的二维图像。LIN等[4]运用应力波断层成像技术,对不同大小的樟树Cinnamomum camphora人造空洞进行了研究,表明应力波断层成像技术可以反映木材的空洞大小及位置,验证了应力波无损成像技术的可行性。以上研究主要针对应力波传播规律与成像技术的应用领域,但对应力波成像算法的研究不多。FENG等[5]提出了一种基于插值的图像重建算法,利用周围点的值估计未知点的速度情况,进而重建树木内部缺陷图像。DU等[6]提出了基于椭圆的空间插值以及速度补偿的方法对木材的横截面进行成像,能够准确地判定缺陷位置,并得到了良好的成像效果。国内外对于应力波在树木内部的传播规律做了大量的研究[8-10],但是对于木材径切面传播规律和成像方法研究还比较少。翁翔等[7]利用无损检测仪器对木材的径切面进行测量,对应力波横向速度和径向速度的比值进行了回归,得到了应力波在木材径向上的传播速度模型,其结果的拟合度较高。但在木材径切面的缺陷成像方法研究较少。鉴于此,本研究基于应力波在木材径切面中沿直线上传播的假设,提出一种木材径切面内部缺陷的速度修正成像方法。

    • 木材可以用圆柱体来近似的表示,那么木材的径切面可以被抽象成一个近似的矩形模型,将n个(n一般为12)传感器放置于被测树木两侧(传感器1至n/2位于同一侧)。每个传感器距离地面的高度(或者某一水平高度的距离)分别为(H1H2H3,…,Hn);S1~Sn表示应力波传感器,传感器按照从右向左,呈U字形进行放置;D表示树木直径;θjk为每个传感器两两之间的相对角度。如图 1所示:定义一个坐标系表示树木径切面。

      图  1  木材径切面模型

      Figure 1.  Radial and longitudinal plane model of wood

      Sn/2+1所在的位置为原点坐标,建立传感器的坐标系。

      $$ S\left( {{X_i}, {Y_i}} \right) = \left\{ \begin{array}{l} {X_i} = D/2, {Y_i} = {H_i} - {H_{n/2 + 1}}, \;\;1 \le i \le n/2\\ {Y_i} = 0, {Y_i} = {H_i} - {H_{n/2 + 1}}, \;\;n/2 + 1 \le i \le n \end{array} \right.。 $$ (1)

      式(1)中:有1≤inXi表示传感器位置横坐标,Yi表示传感器位置纵坐标。根据每个传感器的位置数据,可得到每个传感器两两之间的相对角度θjk为:

      $$ {\theta _{jk}} = \left\{ {\begin{array}{*{20}{l}} {\arctan (\frac{{{H_j} - {H_k}}}{D})}&{{H_j} > {H_k}, {H_k} < {H_{n/2 + 1}}}\\ {\arctan (\frac{{{H_j} - {H_{\left( {n/2 + 1} \right)}} - \left( {{H_k} - {H_{\left( {n/2 + 1} \right)}}} \right)}}{D}), }&{{H_j} > {H_k}, {H_k} < {H_{n/2 + 1}}}\\ {\arctan (\frac{{{H_k} - {H_j}}}{D})}&{{H_j} < {H_k}} \end{array}} \right.。 $$ (2)

      式(2)中:1≤jn/2,n/2+1≤kn

      假定应力波在木材中是沿着直线传播的,通过传感器间的距离和应力波传播的时间计算出传感器间应力波的传播速度[11]

    • 为了预测木材离轴的单轴压缩强度,Hankinson公式[12]给出一种数学模型。当应力波在木材径切面上传播时,如果传播方向与垂直于木材纹理角方向所成的角度为θ,那么根据Hankinson公式可以推出如下公式:

      $$ {V_\theta } = \frac{{{V_1}{V_{\rm{r}}}}}{{{V_1}{{\sin }^2}\theta + {V_{\rm{r}}}{{\cos }^2}\theta }}。 $$ (3)

      式(3)中:V1表示平行于纹理方向上的应力波传播速度;Vr表示垂直于纹理方向上的应力波传播速度;θ表示传播方向与木材纹理方向所形成的角度;Vθ表示在角度为θ时应力波在木材径切面上的传播速度。

      翁翔等[7]通过计算得出:

      $$ \frac{{{V_\theta }}}{{{V_{\rm{r}}}}} \approx 1 + \frac{{{V_1} - {V_{\rm{r}}}}}{{{V_1}}} \times {\theta ^2}。 $$ (4)

      式(4)中:$ \frac{{{V_1} - {V_{\rm{r}}}}}{{{V_1}}}$是常数。那么如图 2所示:假设在空洞木材中有2个对应角度αβ的传播路径LαLβ,其中有α>βLα经过缺陷区域,其速度表示为VαLβ是经过正常区域的,其速度表示为Vβ。同时假设在未产生空洞的同一木材中,有L'α经过的是正常的区域,其速度为V'α。那么有:

      图  2  应力波在不同缺陷程度的木材中传播

      Figure 2.  Stress waves propagate in wood of different defect degree

      $$ V{'_\alpha } > {V_\beta }, V{'_\alpha } > {V_\alpha } $$ (5)

      根据式(5),可以得出:

      $$ \frac{{V{'_\alpha }}}{{1 + z \times {\alpha ^2}}} \approx \frac{{V{'_\beta }}}{{1 + z \times {\beta ^2}}} \approx {V_{\rm{r}}}。 $$ (6)

      式(6)中:$ z = \frac{{{V_1} - {V_{\rm{r}}}}}{{{V_1}}}$。由式(5)和式(6)可以得出:

      $$ \frac{{{V_\alpha }}}{{1 + z \times {\alpha ^2}}} < \frac{{{V_\beta }}}{{1 + z \times {\beta ^2}}}。 $$

      就可以得出速度修正的公式:

      $$ {V_{{\rm{c}}\left( {{\theta _{jk}}} \right)}} = \frac{{{V_{{\rm{d}}\left( {{\theta _{jk}}} \right)}}}}{{1 + z \times {\theta _{jk}}^2}}。 $$ (7)

      式(7)中:Vc(θjk)表示修正后的速度;Vd(θjk)表示为测量得到的应力波传播速度;k是常数;θjk表示2个传感器间与垂直于木材方向纹理所成的角度。

      因此经过缺陷区域的Vα经过速度修正后,就可以反映出应力波经过缺陷区域速度呈现出降低的现象,有效地减少了由木材的各向异性[13]导致的应力波在木材不同方向上传播速度不一致现象的影响,为木材内部缺陷成像提供更精确的数据。

      将通过传感器之间的线速度VL转换成木材径切面上交点的速度VP

      $$ {V_{\rm{P}}} = \max \left( {{V_{{\rm{L1}}}}, {V_{{\rm{L2}}}}, \cdots , {V_{{\rm{L}}n}}} \right)。 $$ (8)

      式(8)中:VL1VLn分别表示每条经过该交点的应力波传播路径上修正后的线速度;(VL1, VL2, …, VLn)∈Vc(θjk)

      利用速度修正插值法对木材的径切面进行成像。每个已知点都对插值点具有一定的影响,即权重[14]。权重随着已知点和插值点之间距离的增加而减弱,距离插值点越近的已知点的权重越大。而且当已知点在距离插值点一定范围以外时,权重可以忽略不计。在任一待插值点的值是邻域内已知点权重之和。可以表示为:

      $$ {V_{{\rm{q}}i}} = \left( {\sum\limits_{i = 1}^n {\frac{{\max \left( {\frac{{\left( {V{'_{{\rm{L1}}}}} \right)}}{{1 + z \times {\theta _{{\rm{L1}}}}}}} \right), \left( {\frac{{\left( {V{'_{{\rm{L2}}}}} \right)}}{{1 + z \times {\theta _{{\rm{L2}}}}}}} \right), \cdots , \left( {\frac{{\left( {V{'_{{\rm{L}}n}}} \right)}}{{1 + z \times {\theta _{{\rm{L}}n}}}}} \right)}}{{{d^m}_i}}} } \right)/\sum\limits_{i = 1}^n {\frac{1}{{{d^m}_i}}} 。 $$ (9)

      式(9)中:Vqi为待插值点的预估值;V'Ln为每条经过该交点的应力波传播路径上未修正的线速度,(V'L1, V'L2, …, V'Ln)∈Vd(θjk)θLn表示每条经过该交点的应力波传播路径与水平面的所成的夹角,(θL1, θL2, …, θLn)∈θjkdi为第i个已知点与待插值点的距离;n表示邻域内参与插值计算的已知点的个数;m为常数。

      成像算法步骤:

      步骤1:初始化(建立木材径切面的模型,输入传感器的坐标(x, y),计算传感器间的线速度Vd(θjk)

      步骤2:利用式(7)求出修正后的线速度矩阵Vc(θjk)

      步骤3:传感器两两连线的交点坐标(xc, yc),将邻域内已知属性值的交点的坐标以及每个交点对应的线速度V'Ln代入式(8)计算Vp的属性值;

      步骤4:通过将已知点的属性值代入式(9),计算未知点的属性值;

      步骤5:重复步骤4,直到计算出所有径切面上的每个点的属性值。

      步骤6:将预估点根据属性值的不同进行不同颜色的赋值,并生成二维图像。

    • 实验用于采集木材应力波数据的仪器是Fakopp。如图 3所示:将传感器分别放置在原木径切面两边,每边6个传感器,总共12个。实验过程中使用检测仪器自带的重锤以相同的力度,垂直于传感器方向,对各个传感器分别敲击3~5次,采集总共36~60组数据。将数据用于下文的成像中。

      图  3  应力波测量装置与传感器的布置

      Figure 3.  Stress wave measuring device and sensor arrangement

      本实验样本编号、样本种类、样本周长、样本测量高度等信息如表 1所示。

      表 1  实验样本基本信息

      Table 1.  Experimental sample information

      样本编号样本种类样本周长/cm检测区域面积/cm2缺陷面积/cm2
      1雪松Cedrus deodara62.831 000111.20
      2泡桐Paulownia fortunei73.831 175427.50
      3樟树Cinnamomum camphora60.95940314.15

      对样本的空洞面积测量方法如下:利用卷尺测量树木空洞的长和宽,假设样本的长为L,宽为W,则对于长方形的空洞缺陷面积SR计算公式为:

      $$ S_{\rm R}=LW。 $$ (10)

      假设梯形的空洞,上底长度为A,下底长度为B,缺陷的高度为H,梯形的空洞的面积ST则用公式可以表示为:

      $$ {S_{\rm{T}}} = \frac{{\left( {A + B} \right)H}}{2}。 $$ (11)

      假设圆形的缺陷孔洞的直径为R,那么缺陷的面积SY可以表示为:

      $$ S_{\rm Y}=\mathit{π} R^2。 $$ (12)
    • 图 4A可以看出:缺陷为人工挖掘并且形状为长方形。如图 4B所示:未修正的线图将经过缺陷区域的线段表示为黄色,未正确区分缺陷区域。图 4C表示的是修正后的线图,将经过缺陷区域的线段(如1号和7号、2号和8号等传感器之间的连线)错误地表示为缺陷颜色红色。从图 4D可以看出:使用IDW的成像结果,只能粗略反映缺陷位置,对于缺陷大小与形状都不能进行准确成像。对比VCI的成像结果图 4E,VCI可以较好还原缺陷位置,可以看出缺陷的位置是处于第3号、5号与8号、10号传感器之间,与图 4A所示的真实缺陷位置基本吻合。但是对于缺陷形状的成像不够准确。

      图  4  雪松样本实验结果

      Figure 4.  Sample experiment results of Cedrus deodara sample

      图 5A可以看出:样本缺陷位置处于右边2号至5号传感器的下面3 cm,以及左边11号传感器至7号传感器的下面3 cm的位置,形状为矩形。从未修正线图(图 5B)可以看到:仅有少部分经过缺陷区域的线段被表示为缺陷的红色,大部分经过缺陷区域的线段被表示为非缺陷的绿色。修正后的线图(图 5C)较好地区分缺陷区域,并以不同的颜色进行区分。可以看到IDW的成像结果虽然能够还原缺陷的位置(图 5D),但对缺陷还原的准确度较低,对缺陷位置和大小的还原参考价值较低。如图 5E所示:VCI的成像结果可以较好地反映缺陷位置,但存在对缺陷区域成像结果比真实结果偏小的情况,以及对缺陷边缘的重建不够平整,对于样本的缺陷轮廓可以较好还原,可以反映样本的缺陷特征。

      图  5  泡桐样本实验结果

      Figure 5.  Sample experiment results of Paulownia fortunei sample

      图 6A可以看出:树木的缺陷位置位于右边的3号和4号传感器之间,以及左边的9号和10号传感器之间。对比图 6B图 6C可以看到:未修正的线图将经过缺陷区域的线段(如7号、1号传感器线段以及6号、12号传感器线段)错误地表示为经过正常区域的线段,呈现绿色,与真实情况不符。而修正后的线图,能够较好区分缺陷部位。从IDW成像结果(图 6D)反映缺陷所处的低速度区域,但对于缺陷的形状无法正确成像。缺陷成像结果呈现不规则的块状分布,无法反映样本的缺陷特征。VCI的成像结果(图 6E)显示:缺陷形状呈现一个类似圆形的不规则形状,与真实的缺陷形状较为接近。成像结果表明缺陷位于3号、4号传感器与9号、10号传感器之间,与图 6A的真实缺陷位置也基本相同。

      图  6  樟树样本实验结果

      Figure 6.  Sample experiment results of Cinnamomum camphora sample

    • 使用混淆矩阵法对成像效果进行定量的误差分析,混淆矩阵包含由分类系统完成的实际和预测的分类信息[15]。FAWCETT[16]描述了真、假分类系统中4种可能组合(表 2)。

      表 2  缺陷与预测的4种组合

      Table 2.  Four combination of real defects and predict defects

      类别预测缺陷预测完好
      真实缺陷被正确地预测为缺陷的面积C被错误地预测为完好的面积E
      真实完好被错误地预测为缺陷的面积F被正确地预测为完好的面积W

      表 2的矩阵可以推出的常用度量值[17],如准确率或者正确率(A),精确度(P)和查全率(R)。

      $$ A = \frac{{C + W}}{{C + E + F + W}}; $$ (13)
      $$ P = \frac{C}{{C + F}}; $$ (14)
      $$ R = \frac{C}{{C + E}}。 $$ (15)

      式(13)~(15)中:准确率A表示对木材状态预测正确的百分比;精确度P表示对木材缺陷预测的可靠性水平高低;查全率R表示对木材缺陷的预测能力高低;CWEF所代表含义见表 2

      表 3可见:3个样本的准确率A为72.32%~95.85%。使用VCI方法的3个样本的平均准确率为93.34%,使用IDW方法的3个样本平均准确率为82.63%。IDW方法在1号、2号样本中的准确率均大于85.00%,其原因应该是IDW方法在成像结果中虽然不能很好地反映木材缺陷的特征,但是对木材完好的部位成像结果较好,所以得到较高的准确率。而3号样本只有72.32%的准确率,对比真实缺陷可以看出:IDW方法将部分正常木材错误地成像为缺陷木材,所以导致准确率下降。IDW和VCI的准确度标准差分别为7.38%和2.04%。可见,VCI在准确率上的标准差较小,说明VCI方法在不同样本中的准确率结果更稳定,准确率相较于IDW方法均有所提高。

      表 3  不同样本的准确率

      Table 3.  Accuracy of different samples

      样本编号准确率A/%
      IDW方法VCI方法
      189.2095.85
      286.3693.30
      372.3290.86

      表 4可见:3个样本的精确度P为65.79%~88.35%。使用VCI方法的3个样本的平均精确度为82.26%,使用IDW方法的3个样本平均精确度为70.19%。VCI方法中3号样本的结果要好于其他样本,其原因是3号样本中的预测缺陷大部分与真实缺陷重合,且预测缺陷大部分落在真实缺陷上,所以得到较高的精确度。IDW和VCI方法的精确度标准差分别为5.00%和4.70%。VCI方法在准确度上的标准差较小,VCI方法在精确度的稳定性略高于IDW方法。同时VCI方法在不同样本上的精确度均好于IDW方法。

      表 4  不同样本的精确度

      Table 4.  Precision of different samples

      样本编号精确度P/%
      IDW方法VCI方法
      167.6076.92
      277.1888.35
      365.7981.51

      表 5可见:3个样本的查全率R为17.98%~94.03%。使用VCI方法的3个样本的平均查全率为92.65%,使用IDW方法的3个样本平均查全率为36.64%。VCI方法中的查全率高于IDW方法。可见,VCI方法可以较好地反映木材的真实缺陷结果,且对木材的缺陷位置成像结果更好。在真实缺陷大小不变的情况下,由于VCI方法可以更好地反映木材的缺陷情况,因此得到了更高的查全率,与真实的成像结果相吻合。IDW和VCI方法的查全率标准差分别为:15.59%和1.93%。说明VCI方法的查全率标准差要远小于IDW方法,在查全率的稳定性上要好于IDW方法。

      表 5  不同样本的查全率

      Table 5.  Recall rate of different samples

      样本编号查全率R/%
      IDW方法VCI方法
      117.9889.92
      256.1494.03
      335.8094.00
    • 本研究提出了一种应用于木材无损检测领域的,基于应力波的木材径切面成像方法。通过与IDW的成像效果比较发现:不论是在定量还是定性方面,VCI方法对木材径切面的缺陷成像结果都好于IDW方法。与真实的缺陷进行对比得到较小的误差,可以准确反映木材内部的真实缺陷情况,具有较高的成像精度。利用传感器测量得到的不同应力波数据,通过对不同缺陷的样本进行实验,验证了该算法对于木材的径切面成像效果的可行性。利用混淆矩阵对成像结果进行了定量分析发现:VCI方法的平均准确度为93.34%,平均精确度为82.26%,平均查全率为92.65%,均高于IDW方法,但对于缺陷具体形状的成像结果,其准确度有待进一步提高。

参考文献 (17)

目录

/

返回文章
返回