发明名称 基于概率预测的焊缝X射线自动检测伪缺陷剔除方法
摘要 基于概率预测的焊缝X射线自动检测伪缺陷剔除方法,属于焊缝X射线缺陷检测技术领域。本发明针对相邻两帧图像中检出的潜在缺陷区域,根据两帧图像拍摄范围沿焊缝方向位移与垂直于焊缝的侧向位移,对前帧图像潜在缺陷区域出现在后帧图像中不同位置上的概率进行预测,之后通过分析后帧图像潜在缺陷区域与概率预测结果,对伪缺陷进行剔除,保留实际缺陷位置。本发明能够利用前后帧图像信息实现对伪缺陷的剔除,有利于减少误检,并能够适应不同的检测装置相对焊缝移动速度的要求,并在发生侧向窜动的情况下保证准确检出,避免漏检。
申请公布号 CN103983651B 申请公布日期 2016.02.24
申请号 CN201410222780.4 申请日期 2014.05.23
申请人 清华大学 发明人 邹怡蓉;都东;潘际銮;石涵;邵家鑫;常保华;韩赞东
分类号 G01N23/00(2006.01)I;G01N23/04(2006.01)I 主分类号 G01N23/00(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 一种基于概率预测的焊缝X射线自动检测伪缺陷剔除方法,其特征在于该方法包括以下步骤:步骤1).连续拍摄焊缝X射线图像序列,记同一焊缝X射线图像序列中的相邻两帧图像分别为前帧图像I<sub>k</sub>和后帧图像I<sub>k+1</sub>,对前后两帧图像分别进行焊缝区域提取,分别计算所述后帧图像I<sub>k+1</sub>与所述前帧图像I<sub>k</sub>中的焊缝重心位置,并计算焊缝重心位置沿垂直焊缝方向的位移量,记为<img file="FDA0000510165060000011.GIF" wi="97" he="89" />步骤2).记焊缝在沿焊缝方向相对检测装置的稳定速度为V,拍摄帧频为f,X射线成像视场范围沿焊缝方向长度为L,拍摄得到X射线数字图像沿焊缝方向总像素数为l,按下式计算焊缝沿焊缝方向的位移量,记为<img file="FDA0000510165060000012.GIF" wi="99" he="89" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><msubsup><mi>d</mi><mi>k</mi><mrow><mo>/</mo><mo>/</mo></mrow></msubsup><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mi>Vl</mi><mi>fL</mi></mfrac></mrow>]]></math><img file="FDA0000510165060000013.GIF" wi="192" he="139" /></maths>其中,<img file="FDA0000510165060000014.GIF" wi="81" he="89" />表示沿焊缝方向位移量,<img file="FDA0000510165060000015.GIF" wi="82" he="89" />的矢量方向与焊缝沿焊缝方向相对检测装置运动速度V在图像坐标系中的方向相同;步骤3).在前帧图像I<sub>k</sub>中提取潜在缺陷区域,记为<img file="FDA0000510165060000016.GIF" wi="369" he="95" />其中n为潜在缺陷区域个数;步骤4).在后帧图像I<sub>k+1</sub>中提取潜在缺陷区域,记为<img file="FDA0000510165060000017.GIF" wi="404" he="95" />其中m为潜在缺陷区域个数;步骤5).按下式计算焊缝在图像坐标系中的合成位移量<img file="FDA0000510165060000018.GIF" wi="88" he="89" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><msub><mi>d</mi><mi>k</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mover><msubsup><mi>d</mi><mi>k</mi><mo>&perp;</mo></msubsup><mo>&OverBar;</mo></mover><mo>+</mo><mover><msubsup><mi>d</mi><mi>k</mi><mrow><mo>/</mo><mo>/</mo></mrow></msubsup><mo>&OverBar;</mo></mover></mrow>]]></math><img file="FDA0000510165060000019.GIF" wi="283" he="89" /></maths>将所述前帧图像I<sub>k</sub>中的潜在缺陷区域<img file="FDA00005101650600000110.GIF" wi="342" he="95" />在图像坐标系中平移<img file="FDA00005101650600000111.GIF" wi="94" he="89" />得到预测缺陷区域<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mo>{</mo><msubsup><mi>B</mi><mi>k</mi><mi>i</mi></msubsup><mo>|</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>}</mo><mo>;</mo></mrow>]]></math><img file="FDA00005101650600000112.GIF" wi="365" he="91" /></maths>步骤6).对于后帧图像I<sub>k+1</sub>上图像坐标系中任意(x,y)位置,按照下式计算后帧图像I<sub>k+1</sub>中缺陷区域的预测概率P<sub>k+1|k</sub>(x,y):<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>|</mo><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&Element;</mo><msubsup><mi>B</mi><mi>k</mi><mi>i</mi></msubsup><mo>,</mo><mi>i</mi><mo>&Element;</mo><mo>{</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>}</mo></mtd></mtr><mtr><mtd><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mi>d</mi><mn>2</mn></msup><mrow><mn>2</mn><msup><mi>&sigma;</mi><mn>2</mn></msup></mrow></mfrac></mrow></msup></mtd><mtd><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&NotElement;</mo><msubsup><mi>B</mi><mi>k</mi><mi>i</mi></msubsup><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA00005101650600000113.GIF" wi="925" he="213" /></maths>其中,d为点(x,y)到预测缺陷区域<img file="FDA00005101650600000114.GIF" wi="334" he="94" />的最短距离,σ满足σ∈[5,max(a,b)],其中a为后帧图像I<sub>k+1</sub>长度方向像素数,b为后帧图像I<sub>k+1</sub>宽度方向像素数;步骤7).构造后帧图像I<sub>k+1</sub>的潜在缺陷图像BW<sub>k+1</sub>:将后帧图像I<sub>k+1</sub>中属于潜在缺陷区域<img file="FDA00005101650600000115.GIF" wi="381" he="95" />内的像素赋值为1,将其余像素赋值为0,得到潜在缺陷图像BW<sub>k+1</sub>;步骤8).按下式计算后帧图像I<sub>k+1</sub>中任意(x,y)位置存在缺陷的后验概率R<sub>k+1</sub>(x,y):R<sub>k+1</sub>(x,y)=P<sub>k+1|k</sub>(x,y)·BW<sub>k+1</sub>步骤9).按下式对后帧图像I<sub>k+1</sub>对应的真实缺陷图像Defect<sub>k+1</sub>进行计算:Defect<sub>k+1</sub>(x,y)=R<sub>k+1</sub>(x,y)≥Th其中阈值Th满足Th∈(0,0.5],真实缺陷图像Defect<sub>k+1</sub>中,值为1的位置判断为缺陷区域,值为0的位置判断为不存在缺陷的区域。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室