发明名称 一种基于纹理显著性的微光图像目标检测方法
摘要 本发明涉及一种基于纹理显著性的微光图像目标检测方法。本方法首先提取图像粗糙度特征图,进行二值化特征处理,合并属于同一目标区域,分别计算全局显著性、局部显著性和位置显著性,以此计算纹理显著性获得纹理显著图,最后利用纹理显著图进行图像目标检测。其核心内容是利用微光图像目标纹理与背景差异较大的特点,通过纹理显著性进行目标检测。利用本发明方法探测微光图像目标,目标轮廓性好、击中率高。
申请公布号 CN104346800A 申请公布日期 2015.02.11
申请号 CN201310334161.X 申请日期 2013.08.02
申请人 南京理工大学 发明人 柏连发;张毅;金左轮;韩静;岳江;陈钱;顾国华
分类号 G06T7/00(2006.01)I;G06K9/46(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 朱显国
主权项 一种基于纹理显著性的微光图像目标检测方法,其特征在于,包括以下步骤:步骤一:按如下方法提取图像粗糙度特征图:1.1计算图像中大小为4k×4k活动窗口中像素的平均灰度值,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>A</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>x</mi><mo>-</mo><mn>2</mn><mi>k</mi></mrow><mrow><mi>x</mi><mo>+</mo><mn>2</mn><mi>k</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>y</mi><mo>-</mo><mn>2</mn><mi>k</mi></mrow><mrow><mi>y</mi><mo>+</mo><mn>2</mn><mi>k</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>/</mo><msup><mrow><mo>(</mo><mn>4</mn><mi>k</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00003610142500011.GIF" wi="951" he="154" /></maths>式中,k=1,2,…,L<sub>max</sub>;L<sub>max</sub>为k的最大值,;(i,j)为像数点坐标;f(i,j)是位于(i,j)的像素灰度值;当k=0时,取窗口尺寸大小为3×3;(x,y)为各平均灰度图A<sub>k</sub>的坐标;1.2对于每个像素,分别计算它在水平和垂直方向上窗口之间的平均灰度差,水平方向平均灰度差计算方法如式(2a)所示,垂直方向平均灰度差计算方法如式(2b)所示,E<sub>k,h</sub>(x,y)=|A<sub>k'</sub>(x+ρ,y)‑A<sub>k</sub>(x,y)|     (2a)E<sub>k,v</sub>(x,y)=|A<sub>k'</sub>(x,y+ρ)‑A<sub>k</sub>(x,y)|      (2b)式中,k'=max(k‑L<sub>b</sub>,0),L<sub>b</sub>为两个窗口偏差尺度,L<sub>b</sub>=L<sub>max</sub>‑α且L<sub>b</sub>≥1;参数α取值如下:α=3,L<sub>max</sub>≥5;α=min(2,L<sub>max</sub>‑1),L<sub>max</sub><5;ρ为两个窗口偏心距,ρ=2k'+1;1.3计算每一像素点的最佳尺寸,计算公式如式(3a)所示,S<sub>best</sub>=4k<sub>max</sub>    (3a)E<sub>k</sub>=max(E<sub>k,h</sub>,E<sub>k,v</sub>)    (3b)E<sub>max</sub>=max(E<sub>k</sub>),E<sub>min</sub>=min(E<sub>k</sub>)   (3c)公式(3b)、(3c)为中间变量,式中,E<sub>k</sub>表示每个尺寸k在水平和垂直方向平均灰度差值的最大值,E<sub>max</sub>为像素点最大平均灰度差值,E<sub>min</sub>为像素点最小平均灰度差值,k<sub>max</sub>值按如下三种情形确定:(i)若k=0,E<sub>0</sub>>t<sub>M</sub>,取k<sub>max</sub>=0,其中t<sub>M</sub>为E<sub>0</sub>所有像素点局部非零极大值的均值,否则转入(ii);(ii)若Numel(DE<sub>k</sub><τ<sub>0</sub>)=L<sub>max</sub>‑1且E<sub>max</sub><t<sub>m</sub>,则取k<sub>max</sub>=L<sub>max</sub>,其中Numel表示对满足条件的k计数,DE<sub>k</sub>=|E<sub>k</sub>‑E<sub>k‑1</sub>|,参数<img file="FDA00003610142500021.GIF" wi="258" he="92" />参数<img file="FDA00003610142500022.GIF" wi="338" he="94" />为E<sub>min</sub>的平均值,否则转入(iii);(iii)k<sub>max</sub>=argmax(E<sub>k</sub>);1.4根据图像中每一像素点的最佳尺寸,计算该像素点的局部粗糙度,计算方法如式(4)所示,F<sub>crs</sub>(x,y)=S<sub>best</sub>(x,y)   (4);步骤二:对特征图F<sub>crs</sub>进行二值化处理,并根据格式塔知觉组织规则‑接近律,对属于同一目标的两个区域进行合并;步骤三:依次使用公式(5a)、(5b)计算全局显著性GS<sub>i</sub>A<sub>i</sub>=a<sub>j</sub>,A<sub>min</sub><a<sub>j</sub><A<sub>max</sub>     (5a)<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>GS</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>A</mi><mi>i</mi></msub><mrow><mi>max</mi><mrow><mo>(</mo><msub><mi>A</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mi>b</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA00003610142500023.GIF" wi="964" he="142" /></maths>式中,a<sub>j</sub>为第j特征区域总像素数,A<sub>i</sub>为重新标记后的区域,A<sub>min</sub>是对a<sub>j</sub>限定的最小值,A<sub>max</sub>是对a<sub>j</sub>限定的最大值;步骤四:计算局部显著性LS<sub>i</sub>,计算公式如式(6)所示,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>LS</mi><mi>i</mi></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>|</mo><msub><mover><mi>I</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>-</mo><mover><msub><mi>I</mi><mrow><mi>i</mi><mo>_</mo><mi>surround</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>|</mo><mo>/</mo><mi>M</mi><mo>)</mo></mrow><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00003610142500024.GIF" wi="964" he="114" /></maths><img file="FDA00003610142500025.GIF" wi="44" he="69" />为第i特征区域在粗糙度特征图中的灰度均值,M为特征图最大灰度值;为得到区域局部环境信息,将特征区域沿区域边界外延r宽度像素,<img file="FDA00003610142500026.GIF" wi="162" he="98" />是外延区域在粗糙度特征图中的灰度均值,r∈[5,7];步骤五:计算位置显著性PS<sub>i</sub>,计算公式如式(7)所示,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>PS</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mrow><mo>(</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>y</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>Y</mi><mn>0</mn></msub></mrow><msub><mi>Y</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>x</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>X</mi><mn>0</mn></msub></mrow><msub><mi>X</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mrow><mi>&eta;</mi><mo>/</mo><mn>2</mn></mrow></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00003610142500031.GIF" wi="1052" he="200" /></maths>式中,(X<sub>0,</sub>Y<sub>0</sub>)为图像中心坐标,(x<sub>i0,</sub>y<sub>i0</sub>)第i特征区域中心坐标,η为调节参数;步骤六:计算纹理显著性大小,以获得显著图,计算纹理显著性大小的方法如式(8)所示,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>TS</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mrow><mo>(</mo><msub><mi>GS</mi><mi>i</mi></msub><mo>+</mo><msub><mi>LS</mi><mi>i</mi></msub><mo>+</mo><msub><mi>PS</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00003610142500032.GIF" wi="809" he="142" /></maths>式中,GS<sub>i</sub>,LS<sub>i</sub>,PS<sub>i</sub>均已归一化;步骤七:显著图中的各显著区域即为提取原目标的感兴趣区域,对各感兴趣区域进行注意焦点转移,注意焦点为各感兴趣区域的中心,注意焦点转移的先后次序为区域显著性大小TS<sub>i</sub>的递减顺序;根据抑制返回机制,已经注意过的区域立即被抑制,TS<sub>i</sub>值归零。
地址 210094 江苏省南京市孝陵卫200号