发明名称 基于机器视觉的轴承防尘盖表面缺陷检测方法
摘要 本发明涉及一种基于机器视觉的轴承防尘盖表面缺陷检测方法,其包括如下步骤:a、获取待检测轴承的灰度图像,b、分离得到轴承防尘盖圆环图像;c、进行灰度变换、自适应中值滤波的预处理;d、利用最大类间方差法进行阈值分割,并利用Roberts边缘检测算子提取阈值分割后轴承防尘盖圆环图像的边缘;e、对上述轴承防尘盖圆环图像每隔2°时计算图像区域内数值为1的数目;f、提供模板轴承图像,并计算得到偏转角度σ;g、分离得到轴承防尘盖圆环图像的字符区域与非字符区域;h、对轴承防尘盖圆环图像中对应字符区域、非字符区域根据连通域数量及缺陷面积进行表面缺陷判定。本发明能实现自动检测,降低人工检测目测工作量,提高检测效率,安全可靠。
申请公布号 CN102636490B 申请公布日期 2014.01.15
申请号 CN201210106471.1 申请日期 2012.04.12
申请人 江南大学;无锡信捷电气股份有限公司 发明人 白瑞林;陈文达;吉峰;李新
分类号 G01N21/88(2006.01)I 主分类号 G01N21/88(2006.01)I
代理机构 北京中恒高博知识产权代理有限公司 11249 代理人 夏晏平
主权项 1.一种基于机器视觉的轴承防尘盖表面缺陷检测方法,其特征是,所述轴承防尘盖表面缺陷检测方法包括如下步骤:(a)、将待检测轴承(130)放置于黑色背景上,并获取所述待检测轴承(130)的灰度图像,且以所述灰度图像的左上角为坐标原点建立坐标系;(b)、对上述灰度图像利用边界跟踪、最小二乘法和轴承防尘盖内、外径占待检测轴承(130)外圆半径的比例,分离得到以待检测轴承(130)的圆心为中心的轴承防尘盖圆环图像;(c)、对上述分离得到的轴承防尘盖圆环图像进行灰度变换、自适应中值滤波的预处理;(d)、对上述预处理后的轴承防尘盖圆环图像利用最大类间方差法进行阈值分割,并利用Roberts边缘检测算子提取阈值分割后轴承防尘盖圆环图像的边缘;(e)、对上述轴承防尘盖圆环图像每隔2°时计算图像区域内数值为1的数目,得到1×180的一维向量数据;(f)、提供模板轴承图像,并计算得到上述轴承防尘盖圆环图像与模板轴承图像旋转重合时的偏转角度σ;(g)、根据模板轴承图像的字符区域、非字符区域与上述偏转角度σ,分离得到上述轴承防尘盖圆环图像的字符区域与非字符区域;(h)、对上述轴承防尘盖圆环图像中对应字符区域、非字符区域根据连通域数量及缺陷面积进行表面缺陷判定,并输出轴承防尘盖表面缺陷判断信息;所述步骤(a)中,以黑色绒布(140)为黑色背景,待检测轴承(130)位于检测支架(150)上;待检测轴承(130)的正上方设有安装于检测支架(150)上的CCD相机(100)及位于所述CCD相机(100)正下方的蓝色同轴光源(120);所述CCD相机(100)邻近蓝色同轴光源(120)的端部安装有定焦镜头(110);所述步骤(b)中,包括如下步骤:(b1)、在获取待检测轴承(130)的灰度图像上寻找灰度值大于25的点,得到位于待检测轴承(130)最左上方的边缘点A;(b2)、在边缘点A的右侧、右下侧、下侧及左下侧中寻找邻接边缘点B,并以边缘点B为起点按照与待检测轴承(130)边缘相对应的方向寻找邻接灰度值大于25的边缘点C;(b3)、当边缘点C与边缘点A为同一点时,则上述跟踪边界为待检测轴承(130)的边缘;当边缘点C与边缘点A不是同一点时,则以边缘点C为起点,以与边缘点B至边缘点C相反的方向寻找邻接灰度值大于25的点,直至寻找回到边缘点A为止,以搜索确定待检测轴承(130)的边缘轮廓;(b4)、上述搜索确定待检测轴承(130)边缘轮廓的点集为(x<sub>i</sub>,y<sub>i</sub>),利用最小二乘法得到:F=(N∑x<sub>i</sub><sup>2</sup>-∑x<sub>i</sub>∑x<sub>i</sub>);D=(N∑x<sub>i</sub>y<sub>i</sub>-∑x<sub>i</sub>∑y<sub>i</sub>);E=N∑x<sub>i</sub><sup>3</sup>+N∑x<sub>i</sub>Y<sub>i</sub><sup>2</sup>-∑(x<sub>i</sub><sup>2</sup>+y<sub>i</sub><sup>2</sup>)∑x<sub>i</sub> G=(N∑y<sub>i</sub><sup>2</sup>-∑y<sub>i</sub>∑y<sub>i</sub>)H=N∑x<sub>i</sub><sup>2</sup>y<sub>i</sub>+N∑Y<sub>i</sub><sup>3</sup>-∑(x<sub>i</sub><sup>2</sup>+y<sub>i</sub><sup>2</sup>)∑y<sub>i</sub>从而得到<maths num="0001"><![CDATA[<math><mrow><mi>a</mi><mo>=</mo><mfrac><mrow><mi>HD</mi><mo>-</mo><mi>EG</mi></mrow><mrow><mi>FG</mi><mo>-</mo><msup><mi>D</mi><mn>2</mn></msup></mrow></mfrac><mo>;</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>b</mi><mo>=</mo><mfrac><mrow><mi>HF</mi><mo>-</mo><mi>ED</mi></mrow><mrow><msup><mi>D</mi><mn>2</mn></msup><mo>-</mo><mi>GF</mi></mrow></mfrac><mo>;</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mi>c</mi><mo>=</mo><mo>-</mo><mfrac><mrow><mi>&Sigma;</mi><mrow><mo>(</mo><msup><msub><mi>x</mi><mi>i</mi></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>y</mi><mi>i</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><mo>+</mo><mi>a&Sigma;</mi><msub><mi>x</mi><mi>i</mi></msub><mo>+</mo><mi>b&Sigma;</mi><msub><mi>y</mi><mi>i</mi></msub></mrow><mi>N</mi></mfrac><mo>,</mo></mrow></math>]]></maths>则得到拟合的圆曲线为R<sup>2</sup>=(x-P)<sup>2</sup>+(y-Q)<sup>2</sup>,圆心坐标(P,Q)和外圆半径R:<maths num="0004"><![CDATA[<math><mrow><mi>P</mi><mo>=</mo><mfrac><mi>a</mi><mrow><mo>-</mo><mn>2</mn></mrow></mfrac></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mi>Q</mi><mo>=</mo><mfrac><mi>b</mi><mrow><mo>-</mo><mn>2</mn></mrow></mfrac></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mi>a</mi><mo>=</mo><mfrac><mrow><mi>HD</mi><mo>-</mo><mi>EG</mi></mrow><mrow><mi>FG</mi><mo>-</mo><msup><mi>D</mi><mn>2</mn></msup></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>(b5)、根据待检测轴承(130)的先验知识得知轴承防尘盖的内、外径占待检测轴承(130)外圆半径的比例,以得到待检测轴承(130)的灰度图像中轴承防尘盖的内径与外径,<img file="FDA0000363737020000027.GIF" wi="1332" he="170" /><img file="FDA0000363737020000028.GIF" wi="1309" he="163" />(b6)、根据上述得到待检测轴承(130)的灰度图像中轴承防尘盖内径与外径参数,截取上述灰度图像,以得到轴承防尘盖圆环图像;所述步骤(c)中,包括如下步骤:(c1)、对得到的轴承防尘盖圆环图像采用线性拉伸的灰度变换,以提高轴承防尘盖圆环图像的亮度;所述线性拉伸函数为s=T(w);其中,w为轴承防尘盖圆环图像灰度变换前点(x,y)的亮度值,s为轴承防尘盖圆环图像灰度变换后点(x,y)的亮度值,T为拉伸函数;(c2)、对上述灰度变换后的轴承防尘盖圆环图像进行自适应中值滤波;所述步骤(d)中,包括如下步骤:(d1)、对上述预处理后的轴承防尘盖圆环图像,利用最大类间方差法进行阈值分割,所述分割阈值为T,以获得相应的前景像素点与背景像素点;(d2)、所述前景像素点占轴承防尘盖圆环图像的比例为ω<sub>0</sub>,其平均灰度为μ<sub>0</sub>;所述背景像素点占轴承防尘盖圆环图像的比例为ω<sub>1</sub>,其平均灰度为μ<sub>1</sub>;轴承防尘盖圆环图像的总平均灰度记为<img file="FDA0000363737020000029.GIF" wi="78" he="66" />类间方差记为g;则有<maths num="0007"><![CDATA[<math><mrow><mi>g</mi><mo>=</mo><msub><mi>&omega;</mi><mn>0</mn></msub><msup><mrow><mo>(</mo><msub><mi>&mu;</mi><mn>0</mn></msub><mo>-</mo><mover><mi>&mu;</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>&omega;</mi><mn>1</mn></msub><msup><mrow><mo>(</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>-</mo><mover><mi>&mu;</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>;</mo></mrow></math>]]></maths>(d3)、根据上述公式求取类间方差g的最大值,以得到最优的分割阈值T,并根据获得最优的分割阈值T对轴承防尘盖圆环图像进行分割;(d4)、利用Roberts边缘检测算子检测上述阈值分割后的轴承防尘盖圆环图像,Roberts边缘检测算子的水平、垂直方向检测算子为:<maths num="0008"><![CDATA[<math><mrow><msub><mi>g</mi><mi>x</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><msub><mi>g</mi><mi>y</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>其中,g<sub>x</sub>为水平方向检测算子,g<sub>y</sub>为垂直方向检测算子;所述步骤(e)中,包括如下步骤:(e1)、以轴承防尘盖圆环图像的圆心为中心,水平方向0°开始,每隔2°计算相应轴承防尘盖圆环图像中数值为1的点数目;(e2)、上述进行计算一周后,得到1×180维的数据Ψ1;所述步骤(f)中,包括如下步骤:(f1)、提供模板轴承图像,并以模板轴承图像的圆心为中心,水平方向为0°开始,每个2°计算模板轴承图像中轴承防尘盖内数值为1的点数目,并在进行计算一周后,得到1×180维的数据Ψ2;(f2)、将数据Ψ1内的元素与数据Ψ2内对应的元素做作差后进行平方,并将180个平方值求和,得到误差a<sub>0</sub>;(f3)、将数据Ψ1内的元素每次移动步长p,以形成新数据Ψ1<sup>p</sup>,数据Ψ1<sup>p</sup>与数据Ψ2内对应的元素做作差后进行平方,并将180个平方值求和,得到误差a<sub>p</sub>;所述移动步长p∈[1,2,…,179];(f4)、当移动步长p取到1~179中的任意一个数值后,得到180个误差,所述误差为a<sub>0</sub>、a<sub>1</sub>、a<sub>2</sub>、……、a<sub>179</sub>;(f5)、比较并取上述误差a<sub>0</sub>~a<sub>179</sub>中的最小值,所述最小误差对应的移动步长为p,根据移动步长p得到偏转角度σ=p×2<sup>0</sup>;所述步骤(g)中,根据模板轴承图像中字符区域覆盖的角度(α,β)及偏转角度σ,得到待检测轴承(130)的轴承防尘盖圆环图像中字符区域覆盖的角度(α+σ,β+σ);所述步骤(h)中,包括如下步骤:(h1)、对上述分离出的轴承防尘盖字符区域与轴承防尘盖非字符区域进行区域标记,并根据标记区域内的像素值计算对应区域的面积,以得到轴承防尘盖字符区域内连通域是Φ个,轴承防尘盖非字符区域内连通域是Γ个;(h2)、当轴承防尘盖字符区域内超过面积S的连通域是Θ个,当Θ与Φ不相等时,则轴承防尘盖字符区域表面存在缺陷;(h3)、当轴承防尘盖非字符区域内超过面积S的连通域是Κ个,且Κ不等于0时,则轴承防尘盖非字符区域表面存在缺陷;所述CCD相机(110)采用30万像素的CCD相机;定焦焦头(120)的焦距为25mm。
地址 214122 江苏省无锡市蠡湖大道1800号