发明名称 用于矿山井筒环境下的模糊图像角点提取方法
摘要 本发明涉及一种用于矿井环境下的模糊图像角点提取方法,包括提出了使用内模板对图像进行区域判别,快速识别出平坦域、角点域,避免对平坦域进一步判别;使用最大熵算法结合区域特性对角点域内的图像进行分割为背景区域和前景区域,以便后继根据不同区域选择相应的角点阈值进行角点提取;根据检测点的角点域状态,使用外模板提取出候选角点;对候选角点进行去伪操作,去除边缘上的毛刺点和窄带上满足角点响应函数的伪角点,获得最终的真角点。实验结果表明,本发明在保证实时性的前提下,相比较现有的角点提取算法应用于光照不均的模糊的井筒图像具有更高的提取精度和更好的鲁棒性。
申请公布号 CN103198319B 申请公布日期 2016.03.30
申请号 CN201310124386.2 申请日期 2013.04.11
申请人 武汉大学 发明人 邢远秀;章登义;赵俭辉
分类号 G06K9/46(2006.01)I;G06T7/00(2006.01)I 主分类号 G06K9/46(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 严彦
主权项 一种用于矿井环境下的模糊图像角点提取方法,其特征在于,包括以下步骤:步骤1,使用内模板对图像进行区域判别,识别出平坦域和角点域;步骤1的实现方式为,预定义内模板,所述内模板包含均匀分布在一个圆环上的4个像素,基于内模板以图像中每个像素点为待检测点分别进行以下局部运算,设内模板的中心点位于待检测点X,内模板覆盖的图像上像素点分别为P、P'、Q和Q'点,根据预设的像素灰度差异度阈值T<sub>d</sub>计算内模板所覆盖的四个像素点P、P'、Q和Q'和待检测点X之间的差异度n<sub>d</sub>(P)、n<sub>d</sub>(P')、n<sub>d</sub>(Q)、n<sub>d</sub>(Q'),根据以下函数判断待检测点X在平坦域或角点域,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>F</mi><mrow><mi>C</mi><mi>o</mi><mi>r</mi><mi>n</mi><mi>e</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>n</mi><mi>b</mi></msub><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>n</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>P</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>+</mo><msub><mi>n</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>Q</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>n</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>Q</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>&le;</mo><mn>2</mn></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>n</mi><mi>d</mi></msub><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>n</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>P</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>+</mo><msub><mi>n</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>Q</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>n</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>Q</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>&gt;</mo><mn>2</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000865705630000011.GIF" wi="1172" he="165" /></maths>其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>n</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>|</mo><msub><mi>f</mi><mi>Y</mi></msub><mo>-</mo><msub><mi>f</mi><mi>X</mi></msub><mo>|</mo><mo>&gt;</mo><msub><mi>T</mi><mi>d</mi></msub></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>|</mo><msub><mi>f</mi><mi>Y</mi></msub><mo>-</mo><msub><mi>f</mi><mi>X</mi></msub><mo>|</mo><mo>&le;</mo><msub><mi>T</mi><mi>d</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000865705630000012.GIF" wi="572" he="150" /></maths>Y=P、P'、Q、Q',f<sub>Y</sub>和f<sub>X</sub>为像素点Y和待检测点X的像素灰度值;当F<sub>Corner</sub>(X)为1时,X位于平坦域内,当F<sub>Corner</sub>(X)为0时,X位于角点域内;步骤2,对图像的角点域进行分割得到背景区域和前景区域;步骤2的实现方式为,使用最大熵算法结合区域特性对图像的角点域进行分割,所得背景区域记为A,前景区域记为B;步骤3,针对角点域内的所有检测点分别进行角点域状态判别,判别时采用的两点间状态差异度阈值根据步骤2的分割结果设定;步骤3的实现方式为,对位于角点域的待检测点X,继续确定内模板所覆盖的四个像素点P、P'、Q和Q'与角点域的待检测点X相比的明暗状态S<sub>X→P</sub>、S<sub>X→P'</sub>、S<sub>X→Q</sub>和S<sub>X→Q'</sub>,根据以下角点域状态判别函数判断待检测点X在邻域内位于亮区域或者暗区域,S<sub>X</sub>=Σ(S<sub>X→P</sub>+S<sub>X→P'</sub>+S<sub>X→Q</sub>+S<sub>X→Q'</sub>)其中,状态比较函数为<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mi>X</mi><mo>&RightArrow;</mo><mi>Y</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mrow><msub><mi>f</mi><mi>X</mi></msub><mo>-</mo><msub><mi>f</mi><mi>Y</mi></msub><mo>&lt;</mo><mo>-</mo><msub><mi>T</mi><mi>C</mi></msub></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>|</mo><msub><mi>f</mi><mi>X</mi></msub><mo>-</mo><msub><mi>f</mi><mi>Y</mi></msub><mo>|</mo><mo>&le;</mo><msub><mi>T</mi><mi>C</mi></msub></mrow></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mrow><msub><mi>f</mi><mi>X</mi></msub><mo>-</mo><msub><mi>f</mi><mi>Y</mi></msub><mo>&gt;</mo><msub><mi>T</mi><mi>C</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000865705630000013.GIF" wi="657" he="230" /></maths>Y=P、P'、Q、Q',f<sub>Y</sub>和f<sub>X</sub>为像素点Y和待检测点X的像素灰度值,两点间状态差异度阈值<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>T</mi><mi>C</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>T</mi><mo>/</mo><mn>10</mn></mrow></mtd><mtd><mrow><mi>X</mi><mo>&Element;</mo><mi>A</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>T</mi><mo>/</mo><mn>2</mn></mrow></mtd><mtd><mrow><mi>X</mi><mo>&Element;</mo><mi>B</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000865705630000014.GIF" wi="438" he="150" /></maths>T为区域分割阈值;当S<sub>X</sub>为正数时,位于角点域的待检测点X在邻域内位于亮区域,当S<sub>X</sub>为负数时,位于角点域的待检测点X在邻域内位于暗区域,当S<sub>X</sub>为0时,将内模板顺时针旋转45度,根据新覆盖的四个像素点P、P'、Q和Q'重新计算S<sub>X</sub>;步骤4,根据步骤3所得检测点的角点域状态,使用外模板计算角点判别函数,提取出候选角点;步骤5,对步骤4所得候选角点进行去伪操作,获得最终的真角点。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学