发明名称 基于空时显著性融合的红外运动目标检测方法
摘要 本发明公开了一种基于空时显著性融合的红外运动目标检测方法,用于解决现有基于空时显著性的运动目标检测方法获取的运动目标轮廓清晰度差的技术问题。技术方案是首先利用空域显著性方法得到红外图像的静态显著图,利用时域检测方法得到动态显著图;然后采用加权融合的方式计算空时融合显著图;最后对融合显著图做适当的阈值分割得到检测的运动目标前景结果。本发明方法利用目标的运动特征,突出运动目标的显著性。在OTCBVS公开红外数据库上测试结果表明,本发明方法保留了空域显著性特征中静止目标的显著性,抑制了时域显著性特征中运动目标的晕轮效应。
申请公布号 CN105631898A 申请公布日期 2016.06.01
申请号 CN201510998019.4 申请日期 2015.12.28
申请人 西北工业大学 发明人 郗润平;张福俊;周鑫;张艳宁
分类号 G06T7/20(2006.01)I 主分类号 G06T7/20(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种基于空时显著性融合的红外运动目标检测方法,其特征在于包括以下步骤:步骤一、将红外检测图像的颜色空间转换到Lab颜色空间,使用空域显著性检测方法-频率调谐方法计算检测帧的区域静态显著性特征,得到空域显著图Sal<sup>S</sup>;步骤二、在计算检测帧空域显著性的同时计算该帧图像的时域显著性特征;将检测帧分成子块,标记所述子块为检测帧的一个事件B(r,s,t);相应的在前N‑1帧中对应位置的事件构成集合V(r,s)={B(r,s,t‑1),B(r,s,t‑2),...,B(r,s,t‑N+1)},则该事件的时域显著性通过公式(1)计算得到;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msup><mi>Sal</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>s</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><mi>p</mi><mo>(</mo><mi>B</mi><mo>(</mo><mi>r</mi><mo>,</mo><mi>s</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mi>V</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>s</mi><mo>)</mo></mrow><mo>)</mo><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mo>-</mo><mi>log</mi><mrow><mo>(</mo><mrow><mi>p</mi><mrow><mo>(</mo><mrow><msub><mi>y</mi><mn>0</mn></msub><mo>|</mo><mi>Y</mi></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892451940000011.GIF" wi="1494" he="151" /></maths>构建一个N维向量Y={y<sub>0</sub>,y<sub>1</sub>,...,y<sub>N‑1</sub>},表示每一个子块,在此引入非参核密度估计,估计每一个事件发生的概率分布;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>f</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>K</mi><mi>H</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892451940000012.GIF" wi="1325" he="135" /></maths>式中,H表示带宽矩阵,核函数表示为K<sub>H</sub>(y)=||H||<sup>‑1/2</sup>K(H<sup>‑1/2</sup>y),式(2)具体展开表示为公式(3);<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mover><mi>f</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>K</mi><mi>H</mi></msub><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>|</mo><mo>|</mo><mi>H</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>|</mo><msup><mo>|</mo><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mi>K</mi><mrow><mo>(</mo><mrow><mi>H</mi><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892451940000013.GIF" wi="1468" he="287" /></maths>式中,H(y<sub>i</sub>)反应样本点y<sub>i</sub>的变化导致估计样本点的变化情况,带宽表示为H(y<sub>i</sub>)=h(y<sub>i</sub>)I,h(y<sub>i</sub>)表示为y<sub>i</sub>到k‑th近邻点的欧氏距离;采用对称多变量高斯密度函数作为核函数;<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><mi>f</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mrow><mo>(</mo><mrow><mn>2</mn><mi>&pi;</mi></mrow><mo>)</mo></mrow><mrow><mi>q</mi><mo>/</mo><mn>2</mn></mrow></msup><mi>N</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mrow><mo>&lsqb;</mo><mrow><msup><mrow><mo>(</mo><mrow><mi>h</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mrow><mo>-</mo><mi>q</mi><mo>/</mo><mn>2</mn></mrow></msup><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><mrow><mi>h</mi><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>I</mi></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mo>&rsqb;</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892451940000014.GIF" wi="1725" he="143" /></maths>最终检测帧中一个事件块的时域显著性表示为公式(5),<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msup><mi>Sal</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>s</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><mi>p</mi><mo>(</mo><mrow><msub><mi>y</mi><mn>0</mn></msub><mo>|</mo><mi>Y</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mo>-</mo><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><mover><mi>f</mi><mo>^</mo></mover><mo>(</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000892451940000015.GIF" wi="1485" he="178" /></maths>式中,Sal<sup>T</sup>表示检测帧的时域显著性;步骤三、将步骤一计算的静态空域显著性特征和步骤二计算的动态时域显著性特征融合,得到空时域目标显著性特征。
地址 710072 陕西省西安市友谊西路127号