发明名称 基于压缩采样矩阵分解的红外小目标检测方法
摘要 一种图像处理技术领域的基于压缩采样矩阵分解的红外小目标检测方法,包括以下步骤:步骤一:输入一幅红外小目标图像,求图像信息熵,根据图像信息熵,进行秩估计;步骤二:将图像列向量化,采用沃尔什哈达玛矩阵,压缩采样获得测量值;步骤三:进行图像重构,同时恢复出背景和目标;步骤四:对红外小目标图像进行二值化处理,检测红外小目标。本发明实现了压缩域内检测红外小目标,通过压缩采样,减少了数据处理量和存储量,重构后实现了目标和背景的分离,既得到了目标又得到了背景,所得到的目标部分具有很高的信噪比增益,有效地抑制了背景和噪声,后续处理只需简单的二值化处理就能检测出红外小目标。
申请公布号 CN103871058A 申请公布日期 2014.06.18
申请号 CN201410090156.3 申请日期 2014.03.12
申请人 北京航空航天大学 发明人 李丽;李辉;简伟健;甄红欣
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 代理人
主权项 1.一种基于压缩采样矩阵分解的红外小目标检测方法,其特征在于,包括以下步骤:步骤一:输入一幅红外小目标图像,求图像信息熵,根据图像信息熵,进行秩估计;输入一幅红外小目标图像可以表示为:f(x,y)=B(x,y)+T(x,y)+N(x,y)(1≤x≤m,1≤y≤n)其中,f(x,y)表示一幅n×m的红外图像,B(x,y)表示背景,T(x,y)表示表示目标,N(x,y)表示噪声,(x,y)表示图像中的一个像素点;背景部分缓慢变化,具有很强的相关性,可以看做是一个低秩矩阵,即rank(B)≤r,r是一个常数,r与图像的复杂度密切相关,通常背景越复杂,r值越大。信息熵可以有效地表示图像的复杂度,定义图像信息熵为:<img file="FDA0000476044270000011.GIF" wi="439" he="136" />当p<sub>s</sub>=0,定义p<sub>s</sub>logp<sub>s</sub>=0s表示图像的像素灰度值,p<sub>s</sub>表示该灰度值出现的概率(<img file="FDA0000476044270000012.GIF" wi="215" he="133" />本发明利用线性插值的方法,预先通过一些红外小目标图像,建立信息熵H和r之间的函数关系,进而通过H可以求出对应的r;红外小目标图像中目标很小,占整幅图像的像素比列0.15%左右,具有稀疏性,可以看做是一个稀疏矩阵,即||T||<sub>0</sub>≤K,||||<sub>0</sub>表示l<sub>0</sub>范数,K是一个常数,由小目标的数量和尺寸决定,K<<m×n,矩阵T中大部分元素为零。步骤二:将输入的红外图像列向量化,产生列向量F<sup>N×1</sup>,(N=n×m),然后对F进行压缩采样;y=ΦF=Φ(B+T)+e其中,y表示采样值,Φ∈R<sup>M×N</sup>(M&lt;N)为采样矩阵,e为采样噪声。本发明中采样矩阵为沃尔什哈达玛矩阵。步骤三:步骤二中得到的低维采样值包含了恢复B和T的充足信息,目标检测的任务就转化为从低维的采样值y中求解出低秩矩阵B和稀疏矩阵T,对每一幅红外图像,进行多次测量并恢复,然后进行加权平均;<maths num="0001"><![CDATA[<math><mrow><mfrac><mn>1</mn><msub><mi>N</mi><mi>u</mi></msub></mfrac><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>u</mi></msub></munderover><mi>min</mi><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&Phi;</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>B</mi><mi>i</mi></msub><mo>+</mo><msub><mi>T</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mi>F</mi></msub><mo>,</mo><mi>s</mi><mo>.</mo><mi>t</mi><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>B</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo></mrow><mo>*</mo></msub><mo>&le;</mo><mi>r</mi><mo>,</mo><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>T</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>0</mn></msub><mo>&le;</mo><mi>K</mi></mrow></math>]]></maths>其中,N<sub>u</sub>为测量次数,||||<sub>F</sub>为Frobenius范数,即<img file="FDA0000476044270000022.GIF" wi="302" he="87" />||||<sub>*</sub>为核范数即奇异值之和,在这里,用||B||<sub>*</sub>≤r代替rank(B)≤r。步骤四:对重构恢复出的目标部分T进行二值化,检测红外小目标。<img file="FDA0000476044270000023.GIF" wi="515" he="152" />M为T中的最大像素灰度值,β∈(0,1)为阈值系数。
地址 100091 北京市海淀区学院路37号