专利名称：基于超完备稀疏表示的红外小目标检测方法
技术领域：
本发明涉及的是一种图像处理技术领域的方法，具体是一种基于超完备稀疏表示 的红外小目标检测方法。
背景技术：
随着红外探测技术和红外传感器件的快速发展，该技术在红外天文观测、红外遥 感、红外安全检测、以及红外医疗成像等领域的应用迅速扩展。从上个世纪九十年代起，红 外成像技术就被用来辅助各种环境下的搜救工作，红外成像利用的是目标的辐射能量，因 其具有穿透烟雾、云层、雨雪等的能力，抗干扰能力强，可夜间工作，所以目前在交通和安全 领域，红外弱小目标的检测和跟踪技术受到了越来越多的关注。然而，由于红外传感器受到 大气，海面辐射和作用距离以及探测器噪声等因素影响，使得远距离的目标在红外图像上 尺寸较小，甚至呈现点状，没有明显的纹理和结构特征。此外，图像的信噪比较低，加上背景 通常情况下比较复杂，目标很容易被噪声和背景杂波所淹没，使得红外小目标的检测变得 更加困难。目前基于单帧的红外小目标检测方法主要分为两类基于图像滤波的检测方法和 基于机器学习的检测方法。基于图像滤波的检测方法的基本思想是首先对红外图像的背 景起伏分量进行估计，即背景估计，然后将原始图像与背景起伏分量相减，以得到包含目标 成分和噪声成分的图像，接着通过阈值处理或其他方法得到目标的位置。具体方法包括二 维最小均方滤波、空域高通滤波、空域匹配滤波、中值滤波、形态学滤波、非线性滤波、以及 基于小波变换的滤波等。这些方法相互间的区别仅在于背景估计部分所采用方法的不同， 因此它们又被称为基于背景估计的检测方法。基于机器学习的目标检测方法将目标检测问 题转化为模式分类问题，该类方法不需要对红外图像进行预处理，其基本思想是首先根据 红外目标特性构造目标样本库和背景样本库，接着通过具体的学习方法对目标模型和背景 模型进行训练，最后利用得到的目标模型和背景模型对输入的测试图像进行分类判别，即 分块提取输入图像的子图像，然后根据判别规则直接判别该子图像块为背景图像还是含有 目标的图像。随着机器学习和模式识别技术的快速发展，该类方法越来越受到人们的关注 和重视。经对现有文献检索发现，刘瑞明在2009年用PCA (主成分分析)和KPCA (核主成分 分析)来检测红夕卜小目标(Ruiming Liu. "Eigentargets Versus Kernel Eigentargets Detection ofInfrared Point Targets Using Linear and Nonlinear Subspace Algorithms (特征目标和核特征目标基于线性和非线性子空间的红外点目标检测方 法)”，Infrared Milli Terahz Waves30 :278_293.)，该方法通过对大量的样本进行训练， 利用最终得到的目标模型，将目标检测问题转化为模式分类问题，但对目标模型的学习过 程使得该方法变得十分复杂，实用性不强。
发明内容
本发明的目的在于克服现有技术的上述不足，提供一种基于超完备稀疏表示的红 外小目标检测方法。本发明通过二维高斯模型来生成红外小目标样本，继而构造红外目标 的超完备字典，采用基追踪方法来求解子图像块在超完备字典中的表示系数，最终根据目 标和背景的表示系数的稀疏程度显著不同来完成目标检测的任务，该方法不需要对目标模 型进行训练，不需要任何背景先验知识，具有更强的通用性和实用性。本发明是通过以下技术方案实现的，本发明包括以下步骤第一步，采用二维高斯模型生成W个红外目标样本图像，进而构造红外目标超完 备字典。所述的二维高斯模型是以mXm图像左上角的像素点为原点建立坐标系，得到坐 标为(i，j)的像素点的像素值I (i，j)为 其中Imax是(X。，y0)处的像素值，O χ为水平散布参数，σ y为垂直散布参数，(X。， Y0)为目标的坐标，所述的构造红外超完备字典，包括以下步骤1)设定目标的坐标个数为k，通过调节Imax、Ox和％三个参数生成的样本个数是 q，则得到的红外目标样本图像的个数为W = kXq ；2)分别将每个mXm图像展开为m2 X 1的一维列向量，W个m2 X 1的一维列向量组 成的矩阵就是红外目标的超完备字典。第二步，将测试图像划分为若干mXm的子图像，分别对每个子图像进行提取表示 系数处理，得到每个子图像在红外目标样本图像的超完备字典下的表示系数。所述的提取表示系数处理，是min || α || j s. t. || D α -χ || 2 ^ ε ,其中^eiT2是子图像的向量表示,DeiT2xw^超完备字典，α e Rw是子图像χ在 超完备字典中的表示系数，s. t.表示约束关系，ε是误差，min Il1是I1范数最小化。第三步，对表示系数进行指标化处理，得到每个子图像的稀疏系数。所述的指标化处理，是 其中SI(x)是向量为χ的子图像的稀疏系数，k为设定的目标坐标的个数，Cii(X) 表示在设定的第i个目标坐标下得到的向量为X的子图像的表示系数，1 < i ^ k。第四步，当子图像的稀疏系数大于阈值τ时，则该子图像内存在目标，从而得到 目标的位置。与现有技术相比，本发明的有益效果是不需要对目标模型和背景模型进行训练，只 需解相关的线性规划问题，实现起来更加简单；利用二维高斯模型来生成目标样本，继而构造 超完备目标字典，能够更加有效地抓住目标的内在几何特征，即图像子块在该字典下的表示系 数具有更加显著的差异性，因此，本方法能够更好的抑制背景，突出目标，得到更高的检测率。
图1是实施例的测试图像及其三维显示结果；其中(a)是实施例的测试图像；(b)是实施例的测试图像的三维显示结果。图2是实施例中子图像在超完备字典中的表示系数；其中(a)是不含噪声的目标样本在超完备字典中的表示系数；(b)是图1(a)中 目标4所在的子图像在超完备字典中的表示系数；(c)是背景子图像在超完备字典中的表 示系数。图3是采用现有技术方法和本实施例方法得到的检测结果图及其三维显示图；其中(al)是采用最大均值滤波器得到的检测结果图；(a2)是图(al)的三维显 示图；(bl)是采用最大中值滤波器得到的检测结果图；(b2)是图(bl)的三维显示图；(cl) 是采用Top-Hat滤波器得到的检测结果图；(c2)是图(cl)的三维显示图；(dl)是采用二维 最小均方滤波器得到的检测结果图；(d2)是图(dl)的三维显示图；(el)是采用主成分分 析方法得到的检测结果图；(e2)是图(el)的三维显示图；(fl)是采用本实施例方法得到 的检测结果图；(f2)是图(fl)的三维显示图。图4是实施例中不同的虚警概率下采用不同的检测方法得到的ROC (Receiver operatingcharacteristic，观察者操作特性)曲线示意图。
具体实施例方式下面结合附图对本发明的实施例作详细说明本实施例在以本发明技术方案为前 提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下 述的实施例。实施例本实施例对图1(a)所示的含有10个小目标的测试图像进行检测，图1(b)为图 1(a)所示的测试图像的三维显示结果，包括以下步骤第一步，采用二维高斯模型生成1476个红外目标样本图像，进而构造红外目标超 完备字典。所述的二维高斯模型是以16X16图像左上角的像素点为原点建立坐标系，得到 坐标为(i，j)的像素点的像素值I (i，j)为 其中Imax是(Χ。，y0)处的像素值，σ χ为水平散布参数，σ y为垂直散布参数，(X。， Y0)为目标的坐标，所述的构造红外超完备字典，包括以下步骤1)设定目标的坐标个数为36，通过调节Imax、0；￡和％三个参数生成的样本个数 是41，则红外目标样本图像的个数为W= 36X41 = 1476 ；2)分别将每个16 X 16图像展开为256 X 1的一维列向量，1476个256 X 1的一维 列向量组成的256X1476维矩阵就是红外目标的超完备字典。第二步，采用一个16X16的正方形滑动窗口，依照从上至下、从左至右的顺序逐 像素地提取测试图像的子图像块，并将每个子图像块展开成256维的列向量，分别对每个子图像进行提取表示系数处理，得到每个子图像在红外目标样本图像的超完备字典下的表 示系数。图像块在超完备字典中的表示系数可以通过以下优化问题来求解min Il α || 0 s. t. χ = D α(1)其中Χ是输入的子图像块的向量表示，D是红外小目标超完备字典，α为图像块 χ在超完备字典中的表示系数。由于字典是超完备的，所以α非唯一，压缩传感理论表明
当W, 时，χ在字典D中的稀疏表示α GRn就是唯一的，σ为D列向量线性无关的最
小值。公式⑴是一个NP问题，可用I1范数来替代Itl范数，这样就将公式⑴从一个非凸 问题转变成了凸优化问题min Il α || 丄 s. t. χ = D α(2)在红外小目标检测问题中，输入的测试图像通常会含有噪声和杂波，因此，为了消 除噪声的干扰，得到更加稳定的解，本实施例所述的提取表示系数处理，是min Il α || j s. t. || D α -χ || 2 ^ ε ,其中X e R256是子图像的向量表示，D e R256xi476是超完备字典，α e R1476是子 图像X在超完备字典中的表示系数，S. t.表示约束关系，ε是误差，min Il Il工是I1范数最 小化。本实施例中的ε是子图像块的标准差。第三步，对表示系数进行指标化处理，得到每个子图像的稀疏系数。所述的指标化处理，是 其中SI(x)是向量为χ的子图像的稀疏系数，k为设定的目标坐标的个数(k =36)，Qi(X)表示在设定的第i个目标坐标下得到的向量为χ的子图像的表示系数，
1 ^ i ^ 36 ο第四步，当子图像的稀疏系数大于阈值τ时，则该子图像内存在目标，从而得到 目标的位置。所述的阈值τ e (0，1)。图2是不同子图像在超完备字典中的表示系数，其中图2(a)是不含噪声的目标 样本在超完备字典中的表示系数，图2(b)是图1(a)中目标4所在的子图像在超完备字典 中的表示系数，图2(c)是背景子图像在超完备字典中的表示系数。图2说明子图像是否包 含目标，包含目标和不包含目标在超完备字典中的表示系数有着显著的差异，据此可以轻 易的将目标和背景区分开来。采用Max-Mean (最大均值滤波器)得到的检测结果图如图3 (al)所示，三维显示 图如图3(a2)所示；采用Max-Median (最大中值滤波器)得到的检测结果图如图3 (bl)所 示，三维显示图如图3 (b2)所示；采用Top-Hat滤波器得到的检测结果图如图3 (cl)所示， 三维显示图如图3(c2)所示；采用TDLMS(二维最小均方滤波器)得到的检测结果图如图 3(dl)所示，三维显示图如图3(d2)所示；采用PCA(主成分分析方法)得到的检测结果图 如图3(el)所示，三维显示图如图3(e2)所示；采用本实施例方法得到的检测结果图如图3(fl)所示，三维显示图如图3(f2)所示。比较上述图像可知本实施例方法能更好的抑制 背景，突出目标，得到更高的检测率。图4给出了在不同的虚警概率下采用不同的检测方法得到的检测概率，即 ROC (Receiveroperating characteristic，观察者操作特性)曲线，当曲线以下包含的面积 越大，则证明性能越好，通过观察图4可知本实施例方法得到的曲线最接近左上角，曲线 以下面积最大。采用不同检测方法得到的信噪比如表1所示，由该表可知本实施例方法得到的 数值最大，故再次证明其检测效果最好。表 权利要求
一种基于超完备稀疏表示的红外小目标检测方法，其特征在于，包括以下步骤第一步，采用二维高斯模型生成W个红外目标样本图像，进而构造红外目标超完备字典；第二步，将测试图像划分为若干m×m的子图像，分别对每个子图像进行提取表示系数处理，得到每个子图像在红外目标样本图像的超完备字典下的表示系数；第三步，对表示系数进行指标化处理，得到每个子图像的稀疏系数；第四步，当子图像的稀疏系数大于阈值τ时，则该子图像内存在目标，从而得到目标的位置。
2.根据权利要求1所述的基于超完备稀疏表示的红外小目标检测方法，其特征是，第 一步中所述的二维高斯模型是以mXm图像左上角的像素点为原点建立坐标系，得到坐标 为(i，j)的像素点的像素值I (i，j)为 其中是(x。，y。)处的像素值，ο x为水平散布参数，O y为垂直散布参数，(x。，y。)为 目标的坐标。
3.根据权利要求1所述的基于超完备稀疏表示的红外小目标检测方法，其特征是，第 一步中所述的构造红外超完备字典，包括以下步骤1)设定目标的坐标个数为k，通过调节1_、0；￡和Oy三个参数生成的样本个数是q， 则得到的红外目标样本图像的个数为W = kXq ；2)分别将每个mXm图像展开为m2Xl的一维列向量，W个m2Xl的一维列向量组成的 矩阵就是红外目标的超完备字典。
4.根据权利要求1所述的基于超完备稀疏表示的红外小目标检测方法，其特征是，第 二步中所述的提取表示系数处理，是min Il α Il j s. t. || D α -χ || 2 ^ ε ,其中^eiT2是子图像的向量表示,De兄"2>^是超完备字典，α e Rw是子图像χ在超完 备字典中的表示系数，s. t.表示约束关系，ε是误差，min Il Il工是I1范数最小化。
5.根据权利要求1所述的基于超完备稀疏表示的红外小目标检测方法，其特征是，第三步所述的指标化处理，是 其中SI(X)是向量为X的子图像的稀疏系数，k为设定的目标坐标的个数，CIi(X)表 示在设定的第i个目标坐标下得到的向量为χ的子图像的表示系数，1 ≤i ≤k。
全文摘要
一种图像处理技术领域的基于超完备稀疏表示的红外小目标检测方法，包括以下步骤采用二维高斯模型生成若干红外目标样本图像，进而构造红外目标超完备字典；将测试图像划分为若干子图像，分别对每个子图像进行提取表示系数处理，得到每个子图像在红外目标样本图像的超完备字典下的表示系数；对表示系数进行指标化处理，得到每个子图像的稀疏系数；当子图像的稀疏系数大于阈值τ时，则该子图像内存在目标，从而得到目标的位置。本发明不需要进行训练，实现起来更加简单；能够更加有效地抓住目标的内在几何特征，即图像子块在该字典下的表示系数具有更加显著的差异性，能够更好的抑制背景，突出目标，得到更高的检测率。
文档编号G06T7/00GK101882314SQ20101023087
公开日2010年11月10日 申请日期2010年7月20日 优先权日2010年7月20日
发明者唐峥远, 杨杰, 赵佳佳 申请人:上海交通大学