基于压缩采样矩阵分解的红外小目标检测方法技术

一种图像处理技术领域的基于压缩采样矩阵分解的红外小目标检测方法，包括以下步骤：步骤一：输入一幅红外小目标图像，求图像信息熵，根据图像信息熵，进行秩估计；步骤二：将图像列向量化，采用沃尔什哈达玛矩阵，压缩采样获得测量值；步骤三：进行图像重构，同时恢复出背景和目标；步骤四：对红外小目标图像进行二值化处理，检测红外小目标。本发明专利技术实现了压缩域内检测红外小目标，通过压缩采样，减少了数据处理量和存储量，重构后实现了目标和背景的分离，既得到了目标又得到了背景，所得到的目标部分具有很高的信噪比增益，有效地抑制了背景和噪声，后续处理只需简单的二值化处理就能检测出红外小目标。

Infrared small target detection method based on compressed sampling matrix decomposition

An image of infrared small target detection matrix decomposition method based on the compressed sampling processing technical field, comprises the following steps: input an infrared small target image, the image according to the image information entropy, information entropy, rank estimation; step two: the image column to quantify, using the Walsh Hadamard matrix, compressed sampling the measurement values are obtained; step three: the image reconstruction, and restore the background and objectives; step four: the infrared small target image binarization processing, infrared small target detection. The invention realizes the compression of infrared small target detection domain by compressed sampling, reduce the amount of data processing and storage, the reconstruction is achieved after the separation of the target and the background, not only have the goal has been part of the background, the target with SNR gain is very high, can effectively suppress the background and noise, subsequent processing simply binarization processing can detect infrared small target.

【技术实现步骤摘要】
基于压缩采样矩阵分解的红外小目标检测方法
本专利技术涉及一种数字图像处理
的方法，具体是一种基于压缩采样矩阵分解的红外小目标检测方法。
技术介绍
红外成像技术具有被动，昼夜工作，高穿透力等优点，应用范围极广。红外小目标检测是红外成像技术的一个重要应用，在军事、工业、农业、医学、交通等领域具有重要意义。红外小目标的自动检测和跟踪，可不受夜间条件等限制,因此也受到了国内外的普遍关注和研究。红外小目标占有的像素少，缺少形状，尺寸，结构等纹理信息，图像信噪比低，红外小目标的准确检测已经成为了一个难题。近些年来，出现了多种红外小目标检测的方法，如最大中值滤波，TopHat滤波，以及最新的基于稀疏表达的红外小目标识别方法，这些传统的红外小目标检测方法都是是基于整幅图像处理的。随着压缩感知理论的兴起，如何在压缩域实现红外小目标的检测，即通过处理测量向量来识别红外目标，显得尤为重要。在压缩域内，通过处理压缩采样的测量值，在图像恢复重构的同时，实现目标的探测。压缩感知理论突破奈奎斯特采样定律的限制，其核心思想是在信号采样的同时实现信息的压缩，按照这一思想，信号的采样率不取决于信号的带宽，而是取决于信息在信号中的结构与内容，如果信号足够稀疏，那么其采样频率完全可以低于奈奎斯特采样频率。这一信号处理框架的突出优点是实现了信号采样与压缩的统一，节省了采样、数据传输、存储与处理的开销。同时，压缩传感过程得到的数据包含了原始信号的完整信息，因此研究压缩传感域的信号处理方法是非常有必要的，可以避免不必要的重建过程，从而降低计算和传输开销。压缩感知理论在图像处理领域得到了广泛的应用，如图像压缩，图像融合，人脸识别等，在目标识别领域压缩感知也得到了初步的探索，目标识别的最终目的找到感兴趣的目标并进行识别跟踪等后续处理，如何处理压缩采样得到的低维采样值，直接进行目标检测，即测量值在目标检测中的应用研究，是一个很值得深入研究的问题。
技术实现思路
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值越大。信息熵可以有效地表示图像的复杂度，表示为：当ps＝0,定义pslogps＝0s表示图像的灰度值，ps表示该灰度值出现的概率本专利技术利用线性插值的方法，预先通过一些红外小目标图像，建立信息熵H和r之间的函数关系，进而当知道图像熵H时，便可估计出相应的r。红外小目标很小，占整幅图像的像素比列0.15％左右，具有稀疏性，可以看做是一个稀疏矩阵，即||T||0≤K，||||0表示l0范数，K是一个常数，由小目标的数量和尺寸决定。K＜＜m×n，矩阵T中大部分元素为零。步骤二：将输入的红外图像列向量化，产生列向量FN×1，(N＝n×m)，然会对F进行压缩采样。y＝ΦF＝Φ(B+T)+e其中，y表示采样值，Φ∈RM×N(M＜N)为采样矩阵，e为采样噪声。在本专利技术中，我们用沃尔什哈达玛矩阵作为采样矩阵，因为沃尔什哈达玛变换只有加减运算，而没有乘除运算，可以大大提高运算速度，这在图像图像处理中是很重要的。步骤三：步骤二中得到的低维测量向量包含了恢复B和T的充足信息，当Φ满足约束等容条件(RIP，restrictedisometryproperty)或者秩约束等容条件(RRIP，therank-restrictedisometryproperty)，可以高概率地恢复出B和T，即：然后，目标检测的任务就转化为从低维的采样值中求解出低秩矩阵和稀疏矩阵，表达式如下：min||y-Φ(B+T)||Fs.t||B||*≤r,||T||0≤K其中，||||F为Frobenius范数，即||||*为核范数即奇异值之和。在这里，用||B||*≤r代替了rank(B)≤r，为了提高红外小目标探测的稳定性并抑制采样过程中引入的随机噪声，对每一幅红外图像，进行多次测量并恢复，然后进行加权平均。其中，Nu为测量次数。步骤四：对重构恢复出的目标部分进行二值化，检测红外小目标。M为T中的最大像素灰度值，β∈(0,1)为阈值系数。本方案的优点在于，实现了压缩域内检测红外小目标，通过压缩采样，减少了数据处理量和存储量，重构后实现了目标和背景的分离，既得到了目标又得到了背景，所得到的目标部分具有很高的信噪比增益，有效地抑制了背景和噪声，后续处理只需简单的二值化处理就能检测出红外小目标。附图说明图1为本专利技术的原理流程图图2为重构后的目标和背景A为一幅云-天空背景的红外原始图片B为重构恢复的背景部分C为重构恢复的目标部分D为目标图像的二值化结果图3为采用现有技术和本实施例方法得到的检测结果的三维显示图(A1),(B1),(C1),(D1)：不同背景的红外图像(A2),(B2),(C2),(D2)：最大中值滤波结果(A3),(B3),(C3),(D3)：TopHat滤波结果(A4),(B4),(C4),(D4)：HSMD_ALADM滤波结果(A5),(B5),(C5),(D5)：本专利技术滤波结果具体实施方式为了使本专利技术的目的、技术方案、和优点更加清晰。下面结合附图对本专利技术实施例做详细的阐述。图1为本专利技术的原理流程图。主要包括：图像输入，求图像熵，秩估计，压缩采样，重构恢复，二值化。步骤一：输入一幅红外图像，求图像熵，根据预先得到的熵和秩的函数关系，进行秩估计当ps＝0,定义pslogps＝0步骤二：对列向量化的图像，进行压缩采样y＝ΦF＝Φ(B+T)+e在这里采样率为35％，测量矩阵采用沃尔什哈达玛变换。步骤三：通过测量向量，求解下列表达式，分离出目标和背景，实施结果如图2B，2C所示。对上述表达式的求解我们采用SpaRCS方法求解，具体过程参见文献(Waters,A.E.,Sankaranarayanan,A.C.,&Baran本文档来自技高网...

【技术保护点】
一种基于压缩采样矩阵分解的红外小目标检测方法，其特征在于，包括以下步骤：步骤一：输入一幅红外小目标图像，求图像信息熵，根据图像信息熵，进行秩估计；输入一幅红外小目标图像可以表示为：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值越大；信息熵可以有效地表示图像的复杂度，定义图像信息熵为：

【技术特征摘要】
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值越大；信息熵可以有效地表示图像的复杂度，定义图像信息熵为：当ps＝0，定义pslogps＝0s表示图像的像素灰度值，ps表示该灰度值出现的概率利用线性插值的方法，预先通过一些红外小目标图像，建立信息熵H和r之间的函数关系，进而通过H可以求出对应的r；红外小目标图像中目标很小，占整幅图像的像素比例0.15％左右，具有稀疏性，可以看做是一个稀疏矩阵，即||T||0≤K，...

【专利技术属性】
技术研发人员：李丽，李辉，简伟健，甄红欣，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京,11