发明名称 一种基于视觉的连接器表面特征自动识别方法
摘要 本发明涉及了一种基于视觉的连接器表面特征自动识别方法。本方法采集基准图像和实测图像,然后使用伪Zernike矩对基准图像和实测图像的特征计算不变矩特征向量距离,然后采用鲁棒的SIFT特征匹配算法提供严密的相似特征,最后将伪Zernkie距离系数同SIFT匹配系数一起构成线性可分的二维特征向量,最后对特征向量进行分类,得出相似/不相似的结果。本发明基于视觉的手段,对连接器的表面特征进行非接触的自动识别,适应了连接器质量管控中自动化的发展要求,克服了由打标设备的误差造成的特征差异,在连接器自动化视觉检测方面具有良好的应用前景。
申请公布号 CN103729631B 申请公布日期 2017.01.18
申请号 CN201410017402.2 申请日期 2014.01.15
申请人 中国人民解放军国防科学技术大学 发明人 周朗明;张小虎;关棒磊
分类号 G06K9/00(2006.01)I;G06K9/20(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 湖南省国防科技工业局专利中心 43102 代理人 冯青
主权项 一种基于视觉的连接器表面特征自动识别方法,其特征在于,过程为:采集基准图像和实测图像,然后使用伪Zernike矩对基准图像和实测图像的特征计算不变矩特征向量距离,然后采用鲁棒的SIFT特征匹配算法提供严密的相似特征,最后将伪Zernkie距离系数同SIFT匹配系数一起构成线性可分的二维特征向量,最后对特征向量进行分类,得出相似或者不相似的结果,具体如下:第一步:采集基准图像和实测图像利用工业CCD设备采集具有指定特征的样品图像作为基准图像,实测图像指实际待检产品的图像,设采集到的基准图像和实测图像分别为G(x,y)和G′(x,y);第二步:计算基准图像和实测图像的不变矩特征向量不变矩是图像特征识别中的常用评价指标,即通过提取具有仿射不变的图像特征,采用Zernike矩、伪Zernike矩进行图像相似度的评价,Zernike矩及伪Zernike矩的定义如下:Zernike矩是在一组正交矩,重复率为m的n阶Zernike矩定义为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>A</mi><mrow><mi>n</mi><mi>m</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mi>&pi;</mi></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></msubsup><msubsup><mo>&Integral;</mo><mn>0</mn><mn>1</mn></msubsup><msub><mi>R</mi><mrow><mi>n</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mi>m</mi><mi>&theta;</mi></mrow></msup><mi>f</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>r</mi><mi>d</mi><mi>r</mi><mi>d</mi><mi>&theta;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001096607580000011.GIF" wi="1204" he="126" /></maths>其中,f(r,θ)是灰度图像的极坐标表达函数,r为映射到单位圆内的灰度图像坐标(x,y)在极坐标系下的距离参数,θ为映射到单位圆内的灰度图像坐标(x,y)在极坐标系下的角度参数,伪Zernike不变矩中的基由极坐标(r,θ),|r|<1表示,在二维离散空间中,需要将图像坐标(i,j),i,j=0,1,2...N‑1映射到单位圆内,将其映射到单位圆内的公式为:<maths num="0002"><math><![CDATA[<mrow><msub><mi>A</mi><mrow><mi>n</mi><mi>m</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mi>&pi;</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>x</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><mi>y</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>R</mi><mrow><mi>n</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>jl&theta;</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub></mrow></msup><mi>g</mi><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001096607580000012.GIF" wi="1253" he="157" /></maths>用于特征识别的不变矩为式(8)所表示的伪Zernike矩的模|A<sub>nm</sub>|,由于高阶矩中的图像信息更多,因此本方法中采用6,7,8,9,10阶不变矩[|A<sub>6,0</sub>| |A<sub>7,1</sub>| |A<sub>8,0</sub>| |A<sub>9,1</sub>| |A<sub>10,0</sub>|]作为粗筛选的特征向量,在本步中对基准图像和实测图像进行二值分割,得到特征表达区域的二值影像,记基准图像和实测图像的二值化影像分别为GB(x,y)和GB′(x,y),按照式(8)计算GB(x,y)和GB′(x,y)的高阶伪Zernike不变矩组成的特征向量A=(|A<sub>6,0</sub>| |A<sub>7,0</sub>| |A<sub>8,0</sub>| |A<sub>9,0</sub>| |A<sub>10,0</sub>|),A′=[|A<sub>6,0</sub>′| |A<sub>7,0</sub>′| |A<sub>8,0</sub>′| |A<sub>9,0</sub>′| |A<sub>10,0</sub>′|];第三步:计算不变矩匹配系数δ以不变矩特征向量的欧式距离作为相似度距离评判标准,公式如下式:δ=||A<sub>1</sub>‑A<sub>2</sub>||           (9)其中,A<sub>1</sub>=[|A<sub>6,0</sub>| |A<sub>7,1</sub>| |A<sub>8,0</sub>| |A<sub>9,1</sub>| |A<sub>10,0</sub>|]表示基准图像的伪Zernike不变矩特征向量,A<sub>2</sub>=[|A′<sub>6,0</sub>| |A′<sub>7,1</sub>| |A′<sub>8,0</sub>| |A′<sub>9,1</sub>| |A′<sub>10,0</sub>|]表示实测图像的伪Zernike不变矩特征向量,将该欧式距离作为评判标准同设定的初筛选阈值K<sub>1</sub>进行比较,若小于阈值,则认为特征相似,否则认为特征不相似,计算特征向量A同A′的欧氏距离δ,将其作为不变矩匹配系数,并同设定的初筛选阈值K<sub>1</sub>比较,对大于或等于初筛选阈值K<sub>1</sub>的样本直接判为不相似,结束本次特征识别过程,对于小于初筛选阈值K<sub>1</sub>的样本进入下一步,第四步:计算SIFT特征匹配系数φ对基准图像和实测图像进行SIFT特征提取后,以128维特征向量之间的欧式距离作为同名点匹配的相似度准则对特征点进行匹配,获取基准图像和实测图像的同名点对,仿射变换如下式:<maths num="0003"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msup><mi>x</mi><mo>&prime;</mo></msup></mtd><mtd><msup><mi>y</mi><mo>&prime;</mo></msup></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>x</mi></mtd><mtd><mi>y</mi></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>H</mi><mn>11</mn></msub></mtd><mtd><msub><mi>H</mi><mn>21</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>H</mi><mn>12</mn></msub></mtd><mtd><msub><mi>H</mi><mn>22</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>H</mi><mn>13</mn></msub></mtd><mtd><msub><mi>H</mi><mn>23</mn></msub></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001096607580000021.GIF" wi="1269" he="238" /></maths>在得到n对同名点后,可解算出上式的仿射变换系数,将匹配点转换到统一的坐标系后,计算出同名点的平均距离σ,为了增强不同特征的可分性,将匹配点数考虑进来,构建如下的匹配系数:<maths num="0004"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>&gamma;</mi><mo>=</mo><mn>10</mn><mfrac><mi>&sigma;</mi><mi>n</mi></mfrac></mrow></mtd><mtd><mrow><mi>n</mi><mo>&GreaterEqual;</mo><mn>3</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>&gamma;</mi><mo>=</mo><mn>1000</mn></mrow></mtd><mtd><mrow><mi>n</mi><mo>&lt;</mo><mn>3</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001096607580000022.GIF" wi="981" he="206" /></maths>上式表明,匹配系数与匹配点数成反比,与同名点距离中误差成正比,匹配点数小于3对的时候,将匹配系数赋为一个大值,结合实际影像的特点将该大值设为1000;第五步:分类器判定相似或者不相似按照式(9)和式(11)取不变矩匹配系数δ和SIFT匹配系数γ构成二维的特征向量(δ,φ),利用线性判别函数进行分类得出相似和不相似的结果,完成基准图像与实测图像的特征识别过程。
地址 410073 湖南省长沙市砚瓦池正街47号