发明名称 一种钢板缺陷识别方法
摘要 本发明涉及一种钢板缺陷识别方法,通过对超声回波信号处理得到一波形曲线,对该波形曲线截取处理获取缺陷波形曲线,进而对波形曲线进行计算提取峰值数据,得到一缺陷类数组,对缺陷类数组进行进一步处理,最后利用缺陷阈值的判断方法及形态学分析的方法获取钢板缺陷类型。该钢板缺陷识别方法仅根据超声探头装置获取的回波信号即能自动完成钢板缺陷的识别,计算速度快,识别快速准确,且无需大量的样本数据进行训练,减小了数据计算量,加快了识别速度,从而降低了缺陷识别成本。
申请公布号 CN104502451A 申请公布日期 2015.04.08
申请号 CN201410778077.1 申请日期 2014.12.15
申请人 中国兵器科学研究院宁波分院 发明人 齐子诚;徐向群;唐盛明;乔日东;郭智敏;王晓艳;李红伟;刘子瑜;孙远东;谢宝奎
分类号 G01N29/04(2006.01)I;G01N29/06(2006.01)I;G01N29/44(2006.01)I 主分类号 G01N29/04(2006.01)I
代理机构 宁波诚源专利事务所有限公司 33102 代理人 袁忠卫
主权项 一种钢板缺陷识别方法,其特征在于:包括以下步骤:步骤一、超声仪启动并初始化,上位机向超声仪发送控制参数;步骤二、超声探头装置(1)启动并根据超声仪中的控制参数进行工作,超声探头装置(1)在被测钢板上移动,上位机获取超声探头装置(1)对钢板检测点的实时位置坐标数据,形成检测点的坐标数据组W=[(x<sub>1</sub>,y<sub>1</sub>),(x<sub>2</sub>,y<sub>2</sub>),...,(x<sub>j</sub>,y<sub>j</sub>),...,(x<sub>b</sub>,y<sub>b</sub>)],其中j和b均为正整数,1≤j≤b,b为检测点总数,(x<sub>j</sub>,y<sub>j</sub>)为第j个检测点的位置坐标,x<sub>j</sub>表示第j个检测点在钢板长度方向上的值,y<sub>j</sub>表示第j个检测点在钢板宽度方向上的值;超声探头装置(1)向被测钢板发射超声波,同时接收来自钢板的超声波回波信号,其中,超声波回波信号包括自钢板表面反射的起始信号、自缺陷位置反射的缺陷信号、自钢板底面反射的底波信号;步骤三、超声仪采集并存储超声探头装置(1)返回的各个检测点的回波信号;步骤四、超声仪根据钢板中声速和增益对回波信号进行处理,进而针对每个检测点位置,对应的获取被测钢板范围内回波信号形成的深度‑幅值波形曲线,从而形成深度‑幅值波形曲线数据组A=[A<sub>1</sub>,A<sub>2</sub>,...,A<sub>j</sub>,...,A<sub>b</sub>],其中j和b均为正整数,1≤j≤b,b为检测点总数,A<sub>j</sub>表示第j个检测点的深度‑幅值波形曲线;超声仪将深度‑幅值波形曲线数据组A中的曲线数据上传到上位机中;步骤五、上位机对深度‑幅值波形曲线数据组A进行后续处理;首先,自深度‑幅值波形曲线数据组A中获取每个检测点对应的起始信号幅值,从而形成初始信号幅值数据组I=[I<sub>1</sub>,I<sub>2</sub>,...,I<sub>j</sub>,...,I<sub>b</sub>],自深度‑幅值波形曲线数据组A中获取每个检测点对应的底波幅值,从而形成底波幅值数据组D=[D<sub>1</sub>,D<sub>2</sub>,...,D<sub>j</sub>,...,D<sub>b</sub>];计算初始信号幅值平均值<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><msub><mi>I</mi><mn>1</mn></msub><mo>+</mo><msub><mi>I</mi><mn>2</mn></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>I</mi><mi>j</mi></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>I</mi><mi>b</mi></msub></mrow><mi>b</mi></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000633747740000011.GIF" wi="599" he="140" /></maths>计算底波幅值平均值<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>D</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><msub><mi>D</mi><mn>1</mn></msub><mo>+</mo><msub><mi>D</mi><mn>2</mn></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>D</mi><mi>j</mi></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>D</mi><mi>b</mi></msub></mrow><mi>b</mi></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000633747740000012.GIF" wi="671" he="134" /></maths>其中j和b均为自然数,1≤j≤b,b为检测点总数,D<sub>j</sub>表示第j个检测点在其相应的深度‑幅值波形曲线A<sub>j</sub>上对应的底波信号幅值;对深度‑幅值波形曲线数据组A进行预处理,即对深度‑幅值波形曲线数据组A中每个检测点对应的深度‑幅值波形曲线进行截取,保留被测钢板的表面至底面深度范围内回波信号对应的波形,从而形成缺陷波形曲线数据组B=[B<sub>1</sub>,B<sub>2</sub>,...,B<sub>j</sub>,...,B<sub>b</sub>],其中j和b均为正整数,1≤j≤b,B<sub>j</sub>表示第j个检测点的缺陷波形曲线,b为检测点总数;对缺陷波形曲线数据组B中每个检测点对应的缺陷波形曲线进行求导计算,从而获取相应检测点缺陷波形曲线中的所有波峰点,从而构建波峰点信息数据组C=[C<sub>1</sub>,C<sub>2</sub>,...,C<sub>j</sub>,...,C<sub>b</sub>];C<sub>j</sub>=[C<sub>[j][1]</sub>,C<sub>[j][2]</sub>],C<sub>[j][1]</sub>=(S<sub>j0</sub>,S<sub>j1</sub>,...,S<sub>ji</sub>,...,S<sub>ja</sub>),C<sub>[j][2]</sub>=(F<sub>j0</sub>,F<sub>j1</sub>,...,F<sub>ji</sub>,...,F<sub>ja</sub>);其中j和b均为自然数,1≤j≤b,b为检测点总数,C<sub>j</sub>表示第j个检测点的缺陷波形曲线B<sub>j</sub>中包含的所有波峰点信息数据组;i和a均为自然数,0≤i≤a,a为波峰点总数,C<sub>[j][1]</sub>表示第j个检测点的缺陷深度数组,S<sub>ji</sub>为第j个检测点的缺陷波形曲线B<sub>j</sub>中第i个波峰点对应的缺陷深度值,C<sub>[j][2]</sub>为第j个检测点的缺陷幅值数组,F<sub>ji</sub>为第j个检测点的缺陷波形曲线B<sub>j</sub>中第i个波峰点对应的缺陷信号幅值;整理各个检测点的信息数据,构建一个信息数据库M=[M<sub>1</sub>,M<sub>2</sub>,...,M<sub>j</sub>,...,M<sub>b</sub>],M<sub>j</sub>=[(x<sub>j</sub>,y<sub>j</sub>),C<sub>j</sub>,D<sub>j</sub>],其中j和b均为自然数,1≤j≤b,b为检测点总数,M<sub>j</sub>表示第j个检测点对应的信息数据集;步骤六、待超声探头装置(1)对被测钢板扫查完毕后,针对任一坐标为(x<sub>n</sub>,y<sub>n</sub>)的检测点,其中1≤n≤b,n和b均为正整数,b为检测点总数;根据其相应的缺陷深度数组C<sub>[n][1]</sub>中各波峰点对应的缺陷深度值,对检测点(x<sub>n</sub>,y<sub>n</sub>)相对应的信息数据集M<sub>n</sub>进行整合;即检测点(x<sub>n</sub>,y<sub>n</sub>)对应的缺陷波形曲线B<sub>n</sub>中,依次相邻的两个波峰点中,当后一个波峰点对应的缺陷深度值减去前一个波峰点对应的缺陷深度值的差的绝对值小于深度相关性阈值q时,将后一个波峰点对应的信息数据与前一个波峰点对应的信息数据归为一个子信息数据集;否则,新建一个子信息数据集;以此类推,从而形成检测点(x<sub>n</sub>,y<sub>n</sub>)相对应的子信息数据集群:P<sub>n</sub>=[P<sub>n1</sub>,P<sub>n2</sub>,...,P<sub>nm</sub>,...,P<sub>nk</sub>],其中m和k均为正整数,1≤m≤k,k为检测点(x<sub>n</sub>,y<sub>n</sub>)对应的子信息数据集总数,P<sub>nm</sub>表示检测点(x<sub>n</sub>,y<sub>n</sub>)相对应的子信息数据集群P<sub>n</sub>中第m个子信息数据集;针对子信息数据集群P<sub>n</sub>=[P<sub>n1</sub>,P<sub>n2</sub>,...,P<sub>nm</sub>,...,P<sub>nk</sub>]中的每一个子信息数据集,获取其缺陷深度极值数据组L<sub>n</sub>=[L<sub>n1</sub>,L<sub>n2</sub>,...,L<sub>nm</sub>,...,L<sub>nk</sub>],L<sub>nm</sub>=(MaxC<sub>[nm][1]</sub>,MinC<sub>[nm][1]</sub>),其中m和k均为正整数,1≤m≤k,k为检测点(x<sub>n</sub>,y<sub>n</sub>)对应的子信息数据集总数,L<sub>nm</sub>表示检测点(x<sub>n</sub>,y<sub>n</sub>)相对应的子信息数据集群P<sub>n</sub>中第m个子信息数据集P<sub>nm</sub>中的缺陷深度极值,MaxC<sub>[nm][1]</sub>表示其缺陷深度最大值,MinC<sub>[nm][1]</sub>表示其缺陷深度最小值;将所有检测点对应的波峰点信息数据进行整合后形成新的信息数据库M'=[P<sub>1</sub>,P<sub>2</sub>,...,P<sub>n</sub>,...,P<sub>b</sub>],其中n和b均为自然数,1≤n≤b,b为检测点总数,P<sub>n</sub>表示第n个检测点(x<sub>n</sub>,y<sub>n</sub>)相对应的子信息数据集群;同时获取所有检测点对应的子信息数据集群中各子信息数据集相对应的缺陷深度极值数据集L=[L<sub>1</sub>,L<sub>2</sub>,...,L<sub>n</sub>,...,L<sub>b</sub>],其中n和b均为自然数,1≤n≤b,b为检测点总数,L<sub>n</sub>表示第n个检测点(x<sub>n</sub>,y<sub>n</sub>)相对应的缺陷深度极值数据组;针对任一坐标为(x<sub>n</sub>,y<sub>n</sub>)的检测点相对应的缺陷深度极值数据组L<sub>n</sub>=[L<sub>n1</sub>,L<sub>n2</sub>,...,L<sub>nm</sub>,...,L<sub>nk</sub>]中的每一个缺陷深度极值数据,其中m和k均为正整数,1≤m≤k,k为检测点(x<sub>n</sub>,y<sub>n</sub>)对应的子信息数据集总数;搜索检测点(x<sub>n</sub>,y<sub>n</sub>)所在八领域范围{(x<sub>n</sub>,y<sub>n‑1</sub>);(x<sub>n</sub>,y<sub>n+1</sub>);(x<sub>n‑1</sub>,y<sub>n‑1</sub>);(x<sub>n‑1</sub>,y<sub>n</sub>);(x<sub>n‑1</sub>,y<sub>n+1</sub>);(x<sub>n+1</sub>,y<sub>n‑1</sub>);(x<sub>n+1</sub>,y<sub>n</sub>);(x<sub>n+1</sub>,y<sub>n+1</sub>)}内各检测点相对应的缺陷深度极值数据组中的每一个缺陷深度极值数据,如果缺陷深度极值数据存在交集,则将其对应的检测点相对应的子信息数据集进行合并;如此,将所有检测点的子信息数据集根据缺陷深度极值数据经过搜索整合处理后,构建新的缺陷信息数据库Q=[Q<sub>1</sub>,Q<sub>2</sub>,...,Q<sub>u</sub>,...,Q<sub>z</sub>],其中u和z为正整数,1≤u≤z,z为子缺陷信息数据集总数,Q<sub>u</sub>表示第u个缺陷的子缺陷信息数据集;步骤七、计算缺陷信息数据库Q=[Q<sub>1</sub>,Q<sub>2</sub>,...,Q<sub>u</sub>,...,Q<sub>z</sub>]中各子缺陷信息数据集的相对应的缺陷参数,从而获取缺陷参数数组V=[V<sub>1</sub>,V<sub>2</sub>,...,V<sub>u</sub>,...,V<sub>z</sub>],其中u和z为正整数,1≤u≤z,z为子缺陷信息数据集总数,V<sub>u</sub>表示第u个子缺陷信息数据集的缺陷参数;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>V</mi><mi>u</mi></msub><mo>=</mo><mfrac><mrow><mi>v</mi><mo>&times;</mo><mi>x</mi><mo>&times;</mo><mi>y</mi></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>e</mi><mo>=</mo><mn>1</mn></mrow><mi>v</mi></munderover><mrow><mo>(</mo><mi>Max</mi><msub><mi>C</mi><mrow><mo>[</mo><mi>e</mi><mo>]</mo><mo>[</mo><mn>1</mn><mo>]</mo></mrow></msub><mo>-</mo><mi>Min</mi><msub><mi>C</mi><mrow><mo>[</mo><mi>e</mi><mo>]</mo><mo>[</mo><mn>1</mn><mo>]</mo></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000633747740000031.GIF" wi="675" he="198" /></maths>其中v代表字缺陷数据集Q<sub>u</sub>中包含的检测点总数,x表示钢板长度方向上的两个相邻的检测点中心之间的距离,y表示钢板宽度方向上两个相邻的检测点中心之间的距离,MaxC<sub>[e][1]</sub>表示缺陷数据集Q<sub>u</sub>中第e个坐标点对应的缺陷深度最大值,MinC<sub>[e][1]</sub>表示缺陷数据集Q<sub>u</sub>中第e个坐标点对应的缺陷深度最小值;同时计算信息数据库<img file="FDA0000633747740000032.GIF" wi="504" he="76" />中各子缺陷信息数据集的相对应的平均幅值,从而获取缺陷平均幅值数组<img file="FDA0000633747740000033.GIF" wi="571" he="92" />其中u和z为正整数,1≤u≤z,z为子缺陷信息数据集总数,<img file="FDA0000633747740000034.GIF" wi="65" he="86" />表示第u个子缺陷信息数据集的平均幅值;将缺陷参数数组V=[V<sub>1</sub>,V<sub>2</sub>,...,V<sub>u</sub>,...,V<sub>z</sub>]中的每个缺陷参数与疏松缺陷阈值R<sub>1</sub>进行比较,如果V<sub>u</sub>&lt;R<sub>1</sub>,则判断缺陷参数V<sub>u</sub>对应的子缺陷数据集Q<sub>u</sub>对应的钢板缺陷类型为疏松缺陷;.剔除掉信息数据库Q=[Q<sub>1</sub>,Q<sub>2</sub>,...,Q<sub>u</sub>,...,Q<sub>z</sub>]中疏松缺陷对应的信息数据集,将剩余的各信息数据集相对应的缺陷参数及平均幅值分别与分层缺陷阈值R<sub>2</sub>及幅值阈值f进行比较,如果R<sub>2</sub>&lt;V<sub>u</sub>&lt;1且<img file="FDA0000633747740000041.GIF" wi="184" he="90" />则判断缺陷参数V<sub>u</sub>对应的子缺陷数据集Q<sub>u</sub>对应的钢板缺陷类型为分层缺陷;步骤八、剔除掉信息数据库Q=[Q<sub>1</sub>,Q<sub>2</sub>,...,Q<sub>u</sub>,...,Q<sub>z</sub>]中疏松缺陷和分层缺陷相对应的子缺陷信息数据集,将剩余的各子缺陷信息数据集在二维平面上进行伪彩色图像重构,形成相应的子集图像,然后针对子集图像进行形态学分析,从而获取钢板缺陷类型为裂纹缺陷或者气孔缺陷;计算各子集图像的近圆度T以判断钢板的气孔缺陷,计算各子集图像的细长度G以判断钢板的裂纹缺陷;其中T=S/(π×d<sup>2</sup>/4),其中T表示近圆度,S表示子集图像包含的面积,d表示子集图像中最远两点的距离,当近圆度T≤t时,则认为该子集图像相对应的子数据集所对应的钢板缺陷为气孔型缺陷;G=S/(h×l),其中I表示细长度,S表示子集图像所占据的面积,h表示子集图像中最远两点的距离,l表示子集图像中与h方向相垂直方向上的平均宽度,当细长度G≤g时,则判断该子集图像相对应的子数据集所对应的钢板缺陷为裂纹型缺陷。
地址 315103 浙江省宁波市高新区凌云路199号