发明名称 基于联合显著性的红外目标分割方法
摘要 本发明公开了一种基于联合显著性的红外目标分割方法,用于解决现有基于显著性模型的可见光图像目标分割方法适应性差的技术问题。技术方案是首先将原始图像分成均匀同质的小像素块;然后根据红外图像色彩缺乏、且目标亮于背景的特点,利用像素块之间的灰度对比度和空间对比度建立像素块的全局区域对比度;进而,根据红外目标的边缘对噪声具有鲁棒性的特点,利用梯度信息建立像素块的边缘对比度显著性;接着利用这两种显著性的线性组合得到整幅图像的联合显著性映射;最后通过设定阈值进行二值化得到候选目标区域并且通过相似性滤波进一步去除误分割区域得到最终结果。测试显示,目标分割的准确率超过了90%,召回率接近80%,F-测量超过90%。
申请公布号 CN103413303A 申请公布日期 2013.11.27
申请号 CN201310323539.6 申请日期 2013.07.29
申请人 西北工业大学 发明人 魏巍;张艳宁;张磊;孟庆洁
分类号 G06T7/00(2006.01)I;G06T5/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种基于联合显著性的红外目标分割方法,其特征在于包括以下步骤:步骤一、输入一幅红外图像,使用Mean-Shift聚类方法将原始图像分割成均匀同质的N个像素块;步骤二、对每一个像素块r<sub>i</sub>(i=1,2,...,N),图像中任一个像素块r<sub>k</sub>的区域对比度显著性计算如下:统计得到整幅图像的平均灰度值<img file="FDA00003585797000011.GIF" wi="46" he="75" />和每一个像素块r<sub>i</sub>的灰度直方图,在灰度直方图结果中,<img file="FDA00003585797000012.GIF" wi="92" he="86" />用来表示像素块r<sub>i</sub>的第p种灰度的出现频率;统计每一个像素块r<sub>i</sub>的面积大小w<sub>s</sub>(r<sub>i</sub>),平均灰度值g<sub>i</sub>和质心坐标;计算像素块r<sub>k</sub>和图像中任一像素块r<sub>i</sub>(i=1,2,...,k-1,k+1,...,N)的色彩对比度C(r<sub>k</sub>,r<sub>i</sub>):<maths num="0001"><![CDATA[<math><mrow><mi>C</mi><mrow><mo>(</mo><msub><mi>r</mi><mi>k</mi></msub><mo>,</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mi>k</mi></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>q</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mi>i</mi></msub></munderover><msubsup><mi>f</mi><mi>k</mi><msub><mi>c</mi><mi>p</mi></msub></msubsup><msubsup><mi>f</mi><mi>i</mi><msub><mi>c</mi><mi>q</mi></msub></msubsup><msub><mi>D</mi><mi>c</mi></msub><mrow><mo>(</mo><msub><mi>c</mi><mi>p</mi></msub><mo>,</mo><msub><mi>c</mi><mi>q</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,n<sub>k</sub>,n<sub>i</sub>分别表示像素块r<sub>k</sub>和r<sub>i</sub>内部的灰度种类数,c<sub>p</sub>,c<sub>q</sub>分别表示像素块r<sub>k</sub>和r<sub>i</sub>内部的第p和第q种灰度,D<sub>c</sub>(c<sub>p</sub>,c<sub>q</sub>)表示两种灰度c<sub>p</sub>,c<sub>q</sub>之间差的绝对值;计算像素块r<sub>k</sub>和图像中任一像素块r<sub>i</sub>(i=1,2...,k-1,k+1,...,N)的空间对比度S(r<sub>k</sub>,r<sub>i</sub>):S(<sub>r</sub>k,r<sub>i</sub>)=exp(-D<sub>s</sub>(r<sub>k</sub>,r<sub>i</sub>)/σ<sub>s</sub>) (2)其中,D<sub>s</sub>(r<sub>k</sub>,r<sub>i</sub>)表示两个像素块r<sub>k</sub>和r<sub>i</sub>质心之间的欧式距离,σ<sub>s</sub>控制质心距离对空间对比度的影响;计算像素块r<sub>k</sub>和图像中任一像素块r<sub>i</sub>(i=1,2,...,k-1,k+1,...,N)的区域对比度R<sub>c</sub>(r<sub>k</sub>,r<sub>i</sub>):R<sub>c</sub>(r<sub>k</sub>,r<sub>i</sub>)=w<sub>s</sub>(r<sub>i</sub>)S(r<sub>k</sub>,r<sub>i</sub>)C(r<sub>k</sub>,r<sub>i</sub>)(3)计算像素块r<sub>k</sub>的全局区域对比度,即就是r<sub>k</sub>的区域对比度显著性R<sub>k</sub>:<maths num="0002"><![CDATA[<math><mrow><msub><mi>R</mi><mi>k</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><munder><mi>&Sigma;</mi><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>&NotEqual;</mo><msub><mi>r</mi><mi>i</mi></msub></mrow></munder><msub><mi>R</mi><mi>c</mi></msub><mrow><mo>(</mo><msub><mi>r</mi><mi>k</mi></msub><mo>,</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>,</mo></mtd><mtd><msub><mi>g</mi><mi>k</mi></msub><mo>&GreaterEqual;</mo><mover><mi>g</mi><mo>&OverBar;</mo></mover></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo></mtd><mtd><msub><mi>g</mi><mi>k</mi></msub><mo>&lt;</mo><mover><mi>g</mi><mo>&OverBar;</mo></mover></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤三、对于每一个像素块r<sub>k</sub>,使用Sobel边缘检测算子,得到像素块r<sub>k</sub>边界上每一个像素P<sub>i</sub>(i=1,2,...,M)的梯度绝对值|G<sub>i</sub>|,M为r<sub>k</sub>边界上的像素总个数,边界上的边缘像素点集合为<img file="FDA00003585797000015.GIF" wi="74" he="74" />像素块r<sub>k</sub>的总梯度Sum初始化为Sum=0,设与像素块r<sub>k</sub>在边界像素点P<sub>i</sub>(i=1,2,...,M)处相邻的像素块为r<sub>l</sub>,R<sub>k</sub>和R<sub>l</sub>分别为像素块r<sub>k</sub>和r<sub>l</sub>的区域对比度显著性;依次对每一个边界像素点P<sub>i</sub>作如下处理:如果<img file="FDA00003585797000021.GIF" wi="157" he="83" />如果R<sub>k</sub>>R<sub>l</sub>,则Sum=Sum+|G<sub>i</sub>|否则Sum=Sum-|G<sub>i</sub>|否则,Sum=Sum+|G<sub>i</sub>|计算像素块r<sub>k</sub>的平均梯度<img file="FDA00003585797000022.GIF" wi="96" he="97" /><maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>G</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><mo>&CenterDot;</mo><mi>Sum</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>获得像素块r<sub>k</sub>的边缘对比度显著性E<sub>k</sub>:<maths num="0004"><![CDATA[<math><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><msub><mover><mi>G</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub><mo>/</mo><msub><mi>&sigma;</mi><mi>e</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,σ<sub>e</sub>是控制平均梯度对边缘对比度显著性的影响值;步骤四、计算每一个像素块r<sub>k</sub>的联合显著性U<sub>k</sub>:U<sub>k</sub>=αR<sub>k</sub>+(1-α)E<sub>k</sub> (7)其中,R<sub>k</sub>和E<sub>k</sub>分别为像素块r<sub>k</sub>的区域对比度显著性和边缘对比度显著性;将每一个像素块r<sub>k</sub>的联合显著性U<sub>k</sub>赋值给r<sub>k</sub>内的每一个像素;步骤五、选择固定阈值T,对联合显著性映射图I(x,y)进行二值化处理:<maths num="0005"><![CDATA[<math><mrow><mi>f</mi><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><mo>,</mo></mtd><mtd><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>T</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo></mtd><mtd><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>T</mi><mn>1</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>得到候选目标区域;f(x,y)中,像素为1的区域为目标区域,0为背景区域;步骤六、对于候选目标区域进行相似滤波处理,去除候选目标区域中的非目标部分;对二值化结果中为1的区域进行8邻域的连通域标记;对于每一个连通域C<sub>k</sub>,以其质心为中心重新建立一个面积为原来区域外接矩形4倍的矩形区域A<sub>k</sub>,g<sub>1</sub>为C<sub>k</sub>在原图中的平均灰度,g<sub>2</sub>为A<sub>k</sub>在原图中的平均灰度,设定阈值T<sub>2</sub>=58:如果g<sub>1</sub>-g<sub>2</sub>&gt;T<sub>2</sub>,则C<sub>k</sub>为目标区域,保留;否则C<sub>k</sub>背景区域,去除。
地址 710072 陕西省西安市友谊西路127号