发明名称 一种激光视觉图像的阈值分割方法
摘要 本发明公开了一种激光视觉图像的阈值分割方法,包括对激光视觉图像进行灰度直方图的计算、设置阈值T(k)并设定其初始值T(0)、进行第一次迭代计算、设定第二次迭代过程的阈值初始值T(N)、进行第二次迭代计算、引入权重系数a、确定背景图像和目标图像。本发明步骤简单、设计逻辑合理、可实现,采用了二次迭代算法对激光视觉图像进行分割,其抗干扰能力强、分割效果好,能够满足跟踪系统的实时性要求,能够应用于机器人的焊缝跟踪系统中,使得激光视觉传感技术可以顺利的在工业生产中广泛应用。
申请公布号 CN104835176A 申请公布日期 2015.08.12
申请号 CN201510275837.1 申请日期 2015.05.26
申请人 广西机电职业技术学院;广西柏熙科技有限公司 发明人 莫胜撼;喻宁娜;李宁;梁广瑞;宁显斌;梁耀文
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 广西南宁汇博专利代理有限公司 45114 代理人 邓晓安
主权项 一种激光视觉图像的阈值分割方法,其特征在于,包括以下步骤:步骤一,对激光视觉图像进行灰度直方图的计算,获取256个灰度等级的像素个数P(i),i∈[0,255]的整数;步骤二,设置阈值T(k),其中T(k)∈[0,255],k为≥0的整数;并根据阈值T(k)将激光视觉图像分为目标图像和背景图像:灰度等级大于T(k)的像素点组成为目标图像,灰度等级小于等于T(k)的像素点组成为背景图像;同时,设定阈值T(k)的初始值T(0),其中T(0)∈[100,200]或者T(0)为激光视觉图像中的最大灰度值和最小灰度值的算术平均值;步骤三,按照式(1)和式(2)分别计算目标图像和背景图像的区域内像素的灰度平均值μ<sub>O</sub>(k)和μ<sub>B</sub>(k);<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>T</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mn>1</mn></mrow><mn>255</mn></munderover><mi>P</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>i</mi></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>T</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mn>1</mn></mrow><mn>255</mn></munderover><mi>P</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000724239360000011.GIF" wi="628" he="293" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>T</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></munderover><mi>P</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>i</mi></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>T</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></munderover><mi>P</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000724239360000012.GIF" wi="599" he="285" /></maths>步骤四,按照式(3)计算新阈值T(k+1);<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>T</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>&mu;</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&mu;</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000724239360000013.GIF" wi="690" he="154" /></maths>步骤五,当|T(k+1)‑T(k)|≤δ时,则结束搜索,并将T(k)作为所求的阈值;当|T(k+1)‑T(k)|>δ时,则将新阈值T(k+1)代入到步骤三中,继续迭代计算,直至符合结束搜索的条件;其中,δ为自定义的精度值并且δ≥0;当搜索结束时经历的循环迭代次数设为N,N为>0的整数,此时的阈值为T(N‑1),目标图像区域内的像素的灰度平均值μ<sub>O</sub>(N‑1),此次搜索过程称之为第一次迭代过程;步骤六,令小于等于μ<sub>O</sub>(N‑1)的灰度值的像素点的个数P(i)的值为零;步骤七,进行第二次迭代,并且第二次迭代过程的阈值初始值为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>T</mi><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>[</mo><msub><mi>&mu;</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mn>255</mn><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000724239360000014.GIF" wi="791" he="128" /></maths>将式(4)计算出的阈值T(N)代入步骤三至步骤五的迭代算法中,直至符合结束搜索的条件;当搜索结束时,第二次迭代过程经历的循环迭代次数设为M,M为>0的整数,此时的阈值为T(N+M‑1),目标的灰度平均值为μ<sub>O</sub>(N+M‑1);步骤八,引入权重系数a,a∈(0,1),按式(5)对阈值T(N+M‑1)进行修正,最终得到分割阈值T′:T′=a·T(N+M‑1)+(1‑a)·μ<sub>O</sub>(N+M‑1)  (5);步骤九,将激光视觉图像中灰度等级≤T′的像素点的灰度值设置为0或255,即为背景图像;而激光视觉图像中灰度等级>T′的像素点的灰度值设置为255或0,即为目标图像。
地址 530007 广西壮族自治区南宁市大学东路101号