基于张量扩展的半脊方向滤波器红外小目标识别方法

摘要

基于张量扩展的半脊方向滤波器红外小目标识别方法。本发明首先读取小目标图像；利用余弦函数来构造图像模板；按二维卷积规则将从图像模板上得到的张量矩阵大小扩展为9×9，选取中间5×5个值为滤波矩阵，得到的滤波矩阵构成滤波模板；确定4个滤波方向；设计每个方向上的加权滤波公式；待每个方向上的加权滤波公式遍历完图像后，保存子图；将保存的子图按照矩阵对应点相乘的方法融合成一幅图像，归一化处理后得到背景抑制后的小目标图像。本发明解决了天空、海天背景下的红外小目标识别问题。本发明方法简单、便于硬件并行处理，且漏检率低、信噪比高。适用于空天及海天背景下单个或多个小目标的红外预警、红外制导导弹智能寻的。

说明书

技术领域

本发明属于红外小目标检测技术领域，涉及基于半脊方向滤波器的选取、图 像模板的建立、张量的扩展变形及滤波模板的建立和滤波方向的选取，适用于空 天及海天背景下单个或多个小目标的红外预警、红外制导导弹智能寻的。

背景技术

随着军事高科技的飞速发展，以精确制导导弹为代表的各种精确制导武器在 现代战场中发挥了越来越强大的威力，精确制导武器也对军事要地和飞机、军舰 等高价值作战平台的安全构成严重的威胁。随着工业制造水平和隐身材料的高速 发展及逐步的成熟，隐身材料在军事上的广泛应用，使得雷达越来越不容易检测 到目标。并且由于雷达本身固有的对外辐射特性，使雷达探测系统存在容易暴露 己方作战平台从而容易受到反辐射导弹攻击的缺点。红外探测器自身不辐射能量， 仅以被动的方式接收目标辐射的热能探测目标，且标准大气环境下，地对空的探 测距离能达到100公里，外太空环境下探测距离达到数百公里。因此较雷达探测系 统而言，红外探测系统具有隐蔽性好、抗电磁干扰能力强的优点。因而红外小目 标识别成为现代战争防御系统中一项重要的技术。当前，在红外小目标检测方面， 国内外学者已经提出了多种方法，其中文献“Wang G D，Chen Ch Y，Shen X B. Facet-based in frared small target detection method[J]，Electron ics Letters.Oct.2005， 41(22):218-219.”采用的面模型的匹配滤波法是基于不同时间采集的多帧图像， 图像序列必须配准，不适用于单帧图像目标检测；文献“LùYan，SHI Lin，SU  Xin-zhu.Method for dim small target detection based on wavelet and higher-order  cumulant[J].Infrared Technology，2006，28(12):713-716.”提出的基于小波和高阶 累积量的检测方法虽能大大提高图像信噪比并且保留了目标信息，但小波算法本 身比较复杂；文献“SU Xiu-qin，LIANG Jin-feng，LU Tao，et al.IR target detection &tracking algorithm based on sea-sky background[J].Acta Photonica Sinica，2009， 38(5):1309-1312.”采用一种结合形态学和区域自适应滤波的检测算法，该方法实 时性好，但当目标点相对弱小时，容易在目标周围产生虚警点；专利“基于线性 PCA的红外点目标检测方法”将目标检测问题转化为模式分类问题，但在对目标模 型的学习过程中没有考虑噪声的干扰，使得该方法对复杂背景的适应能力很弱， 且检测过程十分复杂，实用性不强。

综上所述，许多计算量小的红外小目标识别方法，虽实时性比较好，但都是 在特定背景条件下做了一些理想假定，如差分法等，不能有效识别出背景复杂度 高的小目标图像，具有一定的使用局限性。而复杂度较高的方法，如小波变换法 等，虽具有较好的识别效果，但往往实时性比较差，不能满足实际应用的实时性 要求。

发明内容

本发明的目的在于针对上述已有技术的不足，提出一种基于张量扩展的半脊 方向滤波器红外小目标识别方法，在以不增加检测成本的前提下，适用于更多背 景下的红外小目标图像，并提高了红外小目标检测的准确性、实时性和可靠性。

本发明是一种基于张量扩展的半脊方向滤波器红外小目标识别方法，其特征 在于，包括如下步骤：

(1)读取小目标原图，大小修正为256×256；

(2)由于小目标图像点旋转方向不变的特性，根据半脊滤波器构成的图像所具 有的方向滤波特性，利用余弦函数来构成一个左下角灰度值达到最大，右上角灰 度值最小的半脊方向滤波器的图像模板，其图像模板大小为n×n，取n＝2；

(3)利用张量的定义从图像模板上获得与之相关的4个张量矩阵，并将每个张 量矩阵扩展成大小为9×9的矩阵，选取每个矩阵中间的5×5个值为滤波矩阵，把 得到的4个5×5大小的滤波矩阵按照从图像模板中提取张量的相应位置排列构成 滤波模板；

(4)在滤波模板每个5×5大小的滤波矩阵上，分别确定4个滤波方向的角度为 atan(1/2)，π/2-a tan(1/2)，π/2+atan(1/2)，π-atan(1/2)；

(5)在步骤(3)中得到滤波模板的每个5×5大小的滤波矩阵上，分别以个 像素距离为单位在该滤波模板上按照步骤(4)中所述4个方向分别进行取值，每 个方向上取5个点，当单位距离上的点没有值时，用该点相邻最近两点和的均值 代替该点的取值；

(6)设计所述4个方向上的加权滤波公式，在每个方向上分别处理小目标原图 并保存，得到四副不同方向上的小目标子图，具体步骤有：

6.1)选取每个方向上5个点的值，分别设计4个方向上的加权滤波公式；

6.2)然后在4个滤波矩阵上只选用一个方向上的加权滤波公式同时处理小目 标原图并保存，处理后的小目标原图称为该方向上的小目标子图；

6.3)按照6.2)所述方法，在其他三个方向上对小目标图像做同样处理，并分 别保存，共得到四副小目标子图；

(7)按照矩阵对应点相乘的方法对步骤(6)中得到的四副小目标子图进行遍 历相乘，融合成一幅新的图像，并对该副图像进行归一化处理，最终得到背景抑 制后的小目标图像，识别结束。

本发明的技术方案是通过建立图像模板并提取图像模板的张量，将提取出的 张量进行扩展、变形来重新构建滤波模板，然后确定加权滤波的方向，最终得到 基于张量扩展的半脊方向滤波器来实现红外小目标识别的目的。

本发明的实现还在于：步骤(2)中，由余弦函数生成的图像模板Z(n,n)各点 灰度值的表达式为：

      $y\left(i\right)=255\times \cos (\frac{\pi }{2}·\frac{i}{2n-1}),(i=\mathrm{1,2,3},...,2n-1)$      

当i<n时：      Z(n-i+m,m)＝y(i)    (m＝1,2,3,…,i)

当i＝n时：      Z(l,k)＝y(i)     (l＝k＝1,2,3,…,n)

当i>n时：      Z(2n-i-m,n-m)＝y(i)    (m＝0,1,2,…,2n-i-1)

其中n为图像模板的大小，y(i)为灰度值函数。

本发明的该种图像模板的表达式设计，为构造n×n大小的半脊方向滤波器的 图像模板作出了解释，使其图像模板的大小不局限于实时性最好的2×2大小，为 后续当图像模板的大小大于2×2时对识别效果影响的研究进行铺垫。

本发明的实现还在于：其中步骤(3)中，依据张量的定义，对步骤(2)中 形成的大小为2×2的图像模板提取张量矩阵，得到4个张量矩阵，并将每个张量 矩阵按照二维卷积规则扩展至3×3大小，并继续扩展该3×3的矩阵直至该矩阵大 小为9×9结束，得到4个大小为9×9的扩展矩阵，选定每个扩展矩阵的中间5×5个 值为滤波矩阵并按照从图像模板中提取张量的相应位置排列构成滤波模板。

本发明的对图像模板提取出的张量矩阵进行扩展，直至每个张量矩阵扩展到 9×9大小，为本发明所选取的4个滤波方向提供可能，同时选取每个9×9的矩阵 中间5×5个值为滤波矩阵后，降低了算法复杂度。

本发明的实现还在于：其中步骤(6)中，设计4个方向的加权滤波公式，其 中设计的atan(1/2)方向上的加权滤波公式为：

p_new(i,j)＝[T(-2,1)×p(i-2,j+1)+((T(-1,0)+T(-1,1))/2)×(p(i-1,j)+ p(i-1,j+1))/2+T(0,0)×p(i,j)+((T(1,0)+T(1,-1))/2) (p(i+1,j)+p(i+1,j-1))/2+T(2,-1)×p(i+2,j-1)]/K

K＝T(-2,1)+(T(-1,0)+T(-1,1))/2+T(0,0)+(T(1,0)+T(1,-1))/2+T(2,-1)

其中，p(i,j)为小目标图像中的像素点，则i,j为整数且i,j∈[1,256]；p_new(i,j) 为小目标图像该点处的新灰度值；T(i,j)为步骤(3)中所提到的5×5滤波矩阵， 且-2≤i,j≤2；atan(1/2)方向的加权滤波公式设计完毕。

其他三个方向的加权处理公式仅需将按照相应方向变化的p(i,j)和T(i,j)点 坐标分别带入atan(1/2)方向上的加权滤波公式即可。

本发明的该种加权滤波公式的设计具有明显的方向选择性，并且充分考虑做 加权处理的目标点周围的像素对该点的影响，大大降低了漏检率。

与现有技术相比，本发明具有如下优点：

1)本发明中因为采用的图像模板大小仅有2×2，且提出的对从图像模板中提 取张量矩阵的扩展变形算法实现方法简单，所以处理过程占用计算时间少，具有 很好的实时性；

2)本发明中因为构造的半脊方向滤波器是基于对图像模板的提取的张量矩阵 的扩展变形，所以不需要额外增加检测成本，且该方法操作简单便于硬件实现流 水线方式处理；

3)本发明中因为采用大小为atan(1/2)，π/2-atan(1/2)，π/2+atan(1/2)， π-atan(1/2)的角度作为4个滤波方向，不同于普通方向滤波器的滤波方向，而在 每个滤波方向上选定5个以个像素距离为单位的加权滤波点后，使每次运算 包含更多的原图像信息，降低漏检率；

4)实验表明，本发明可以达到良好的红外小目标识别效果。

附图说明

图1为本发明的整体流程框图；

图2为实现本发明红外小目标识别的仿真流程图；

图3为本发明中对张量矩阵扩展过程示意图；

图4为本发明在滤波矩阵上选定的4个滤波方向示意图；

图5(a)为本发明一个实施例原始红外小目标图像；

图5(b)为本发明方法应用于图5(a)的背景抑制后的图像；

图5(c)为已有的剪切波变换算法应用于图5(a)的背景抑制后的图像；

图6(a)为本发明一个实施例原始红外小目标图像；

图6(b)为本发明方法应用于图6(a)的背景抑制后的图像；

图6(c)为已有的剪切波变换算法应用于图6(a)的背景抑制后的图像；

图7(a)为本发明一个实施例原始红外小目标图像；

图7(b)为本发明方法应用于图7(a)的背景抑制后的图像；

图7(c)为已有的剪切波变换算法应用于图7(a)的背景抑制后的图像；

图8(a)为本发明一个实施例原始红外小目标图像；

图8(b)为本发明方法应用于图8(a)的背景抑制后的图像；

图8(c)为已有的剪切波变换算法应用于图8(a)的背景抑制后的图像。

具体实施方式

早期的小目标检测主要是通过雷达检测系统来实施，而因为雷达系统本身的对 外辐射特性，容易受到反辐射导弹攻击。而红外检测设备凭借其隐蔽性强、抗电 磁干扰能力强及被动的工作方式的特点受到越来越多的重视。从上世纪70年代红 外小目标检测技术方案的提出开始，到现在红外小目标检测领域已经成为现代战 争防御系统中一个重要领域。虽然红外小目标检测技术已经得到了长足的发展， 但是现有的许多检测技术在使用过程中仍然存在着实时性差、可靠性差等问题。 本发明针对这类问题提出一种基于张量扩展的半脊方向滤波器红外小目标识别方 法。

实施例1

本发明基于张量扩展的半脊方向滤波器红外小目标识别方法，具体实现参见 图1，包括有如下步骤如下：

(1)读取小目标原图，大小修正为256×256。红外小目标图像往往具有噪声高、 信噪比低、背景复杂、目标所占像素个数少等特点，这些都对从红外图像小目标 图像中识别出小目标带了困难。

(2)由于小目标图像点旋转方向不变的特性，根据ridge滤波器构成的图像所 具有的方向滤波特性，利用余弦函数来构成一个左下角灰度值达到最大，右上角 灰度值最小的半脊方向滤波器的图像模板，其图像模板大小为n×n，经过大量实 验验证，单从检测效果上来看n的取值在大于1小于被检测小目标图像大小的任意 值时检测效果几乎没有差别，但综合考虑实际应用中实时性的要求，取n＝2。即 当采用的图像模板大小为2×2时，处理过程占用计算时间少，具有很好的实时性；

(3)利用张量的定义从图像模板上获得与图像相关的4个张量矩阵，并将每个 张量矩阵扩展成大小为9×9的矩阵，矩阵扩展方式参见图3。选取每个矩阵中间的 5×5个值为滤波矩阵，得到的4个5×5大小的滤波矩阵，将这4个滤波矩阵按照 从图像模板中提取张量的相应位置排列构成滤波模板。

(4)在滤波模板每个5×5大小的滤波矩阵上，分别确定4个滤波方向的角度， 该4个滤波方向的角度分别为atan(1/2)，π/2-atan(1/2)，π/2+atan(1/2)， π-atan(1/2)，4个滤波选取方向简称4个方向，参见图4。

(5)在步骤(3)中得到滤波模板的每个5×5大小的滤波矩阵上，分别以个像素距离为单位在该滤波模板上按照步骤(4)中所述4个方向分别进行取值， 每个方向上取5个点，当单位距离上的点没有值时，用该点相邻最近两点和的均 值代替该点的取值。参见图4，atan(1/2)方向上应选取的5个点的值，这5个点 分别为T(-2,1)、T(0.5,-1)、T(0,0)、T(0.5,-1)、T(2,-1)，其中T为滤波矩阵，T(-2,1)、 T(0,0)、T(2,-1)为该方向上单位距离上已经有值的点，而T(0.5,-1)、T(1,-0.5)两点 的值分别用(T(-1,0)+T(-1,1))/2、(T(1,0)+T(1,-1))/2代替。

(6)在滤波矩阵上设计4个方向上的加权滤波公式，在每个方向上分别处理小 目标原图并保存，得到四副不同方向上的小目标子图，具体步骤有：

6.1)选取滤波矩阵上每个方向上5个点的值，分别设计4个方向上的加权滤 波公式；

6.2)然后在4个滤波矩阵上只选用一个方向上的加权滤波公式同时处理小目 标原图并保存，处理后的小目标原图称为该方向上的小目标子图；

6.3)按照6.2)所述方法，在滤波矩阵上其他三个方向上对小目标图像做同样 处理，并分别保存，加之6.2)中所得的小目标子图共得到四副小目标子图，也就 是说在滤波模板上共得到并保存了4个不同方向的小目标子图，参见图1。

(7)按照矩阵对应点相乘的方法对步骤(6)中得到的四副小目标子图进行遍 历相乘，融合成一幅新的图像，并对该副图像进行归一化处理，最终得到背景抑 制后的小目标图像，识别结束。

在本发明对小目标进行检测和识别的过程中，具体是在利用滤波器的方向选 择上采用了大小为atan(1/2)，π/2-atan(1/2)，π/2+atan(1/2)，π-atan(1/2)的 角度作为4个滤波方向，不同于普通方向滤波器的滤波方向。同时本发明在每个 滤波方向上选定5个以个像素距离为单位的加权滤波点，其中当单位距离上 的点没有值时，用该点相邻最近两点和的均值代替该点的取值，使每次加权滤波 运算过程包含更多小目标图像信息，降低漏检率。

实施例2

基于张量扩展的半脊方向滤波器红外小目标识别方法同实施例1，其中步骤 (2)中，由余弦函数生成的图像模板Z(n,n)各点灰度值的表达式为：

      $y\left(i\right)=255\times \cos (\frac{\pi }{2}·\frac{i}{2n-1}),(i=\mathrm{1,2,3},...,2n-1)$      

当i<n时：     Z(n-i+m,m)＝y(i)      (m＝1,2,3,…,i)

当i＝n时：     Z(l,k)＝y(i)     (l＝k＝1,2,3,…,n)

当i>n时：      Z(2n-i-m,n-m)＝y(i)    (m＝0,1,2,…,2n-i-1)

其中n为图像模板的大小，y(i)为灰度值函数。

实施例3

基于张量扩展的半脊方向滤波器红外小目标识别方法同实施例1-2，其中步骤 (3)中，依据张量的定义，对步骤(2)中形成的大小为2×2的图像模板提取张 量矩阵，得到4个张量矩阵，并将每个张量矩阵按照二维卷积规则扩展至3×3大 小，并继续扩展该3×3的矩阵直至该矩阵大小为9×9结束，得到4大小为9×9的 扩展矩阵，选定每个扩展矩阵的中间5×5个值为滤波矩阵，把得到的4个5×5大 小的滤波矩阵按照从图像模板中提取张量的相应位置排列构成滤波模板。

其中2×2大小的张量矩阵扩展到3×3大小的矩阵公式为：

R₁₁＝r₁₁×r₂₂，

R₁₂＝r₁₁×r₂₁+r₁₂×r₂₂，

R₁₃＝r₁₂×r₂₁，

R₂₁＝r₁₂×r₁₁+r₂₁×r₂₂，

R₂₂＝r₁₁×r₁₁+r₁₂×r₁₂+r₂₁×r₂₁+r₂₂×r₂₂，

R₂₃＝r₁₂×r₁₁+r₂₁×r₂₂，

R₃₁＝r₁₂×r₂₁，

R₃₂＝r₂₁×r₁₁+r₂₂×r₁₂，

R₃₃＝r₁₁×r₂₂，

其中r₁₁，r₁₂，r₂₁，r₂₂为上述大小为2×2张量矩阵中的4个元素，R₁₁，R₁₂， R₁₃，R₂₁，R₂₂，R₂₃，R₃₁，R₃₂，R₃₃为扩展后大小为3×3的矩阵中的新元 素。按照上述张量扩展算法，继续将扩展后的3×3矩阵扩展到5×5大小，并继续 扩展，最终将从图像模板中提取的每个2×2大小的张量矩阵都扩展成大小为9×9 的矩阵结束，得到4个9×9的矩阵，选取每个矩阵中间的5×5个值得到4个滤波 矩阵。扩展方法示意图参见图3，A框内为二维卷积方法示意图；B框内为2×2矩 阵卷积步骤示意图，按B框内所示卷积方法，将大小为2×2矩阵扩展到大小为3×3 的矩阵；C框内大小为3×3的矩阵卷积示意图；继续扩展，最终得到大小为9×9的 矩阵。

因为从图像模板的提取的张量矩阵的扩展变形所构建的该半脊方向滤波器， 不需要额外增加检测成本，所以该方法便于硬件实现流水线方式处理

实施例4

基于张量扩展的半脊方向滤波器红外小目标识别方法同实施例1-3，其中步骤 (6)中，设计4个方向的加权滤波公式，其中设计的atan(1/2)方向上的加权滤波 公式为：

p_new(i,j)＝[T(-2,1)×p(i-2,j+1)+((T(-1,0)+T(-1,1))/2)×(p(i-1,j)+ p(i-1,j+1))/2+T(0,0)×p(i,j)+((T(1,0)+T(1,-1))/2) (p(i+1,j)+p(i+1,j-1))/2+T(2,-1)×p(i+2,j-1)]/K

K＝T(-2,1)+(T(-1,0)+T(-1,1))/2+T(0,0)+(T(1,0)+T(1,-1))/2+T(2,-1)

其中，p(i,j)为小目标图像中的像素点，则i,j为整数且i,j∈[1,256]；p_new(i,j) 为小目标图像该点处的新灰度值；T(i,j)为步骤(3)中所提到的5×5滤波矩阵， 且-2≤i,j≤2；atan(1/2)方向的加权滤波公式设计完毕。其他三个方向的加权处 理公式仅需将按照相应方向变化的p(i,j)和T(i,j)点坐标分别带入atan(1/2)方向 上的加权滤波公式即可获得。将得到的4个方向的加权滤波公式分别运用在滤波 模板上，处理小目标原图，得到每个方向的小目标子图并保存。然后将这四副小 目标子图按照矩阵对应点相乘的办法融合成一幅图像，归一化处理后，得到背景 抑制后的小目标图像。

本发明的具体实施还以通过下边的例子进行说明，该例子中的步骤也是实现 与仿真的过程：

实施例5

基于张量扩展的半脊方向滤波器红外小目标识别方法同实施例1-4，参见图 2,具体实现与执行包括如下步骤：

(1)读取小目标图像，并将其大小修正为256×256；

(2)根据本发明生成图像模板的公式，生成图像模板，以下为本发明生成的图 像模板公式：

      $y\left(i\right)=255\times \cos (\frac{\pi }{2}·\frac{i}{2n-1}),(i=\mathrm{1,2,3}...,2n-1)$      

当i<n时：          Z(n-i+m,m)＝y(i)      (m＝1,2,3,…,i)

当i＝n时：      Z(l,k)＝y(i)       (l＝k＝1,2,3,…,n)

当i>n时：       Z(2n-i-m,n-m)＝y(i)     (m＝0,1,2,…,2n-i-1)

其中n为图像模板的大小，y(i)为灰度值函数，Z(n,n)为图像模板点灰度值的 表达式；

(3)提取图像模板的张量矩阵，并按照二维卷积规则扩展张量矩阵，使每个张 量矩阵大小扩展至9×9大小，并选取每个9×9大小的扩展矩阵的中间5×5个值为 滤波矩阵，把得到的每个5×5大小的滤波矩阵按照从图像模板中提取张量的相应 位置排列构成滤波模板，矩阵扩展方式参见图3；

(4)在每个滤波模板的滤波矩阵上确定4个滤波方向的角度为atan(1/2)， π/2-atan(1/2)，π/2+atan(1/2)，π-atan(1/2)，并设计加权滤波公式，方向选 择方式参见图4。

4.1)将得到滤波模板的每个5×5大小的滤波矩阵上，分别以个像素距离 为单位在该滤波模板上按照所述4个方向分别进行取值，每个方向上取5个点， 当单位距离上没有点时，用该点相邻最近两点和的均值代替该点的取值，每个方 向的相邻最近两点的选择参见图4；

4.2)设计4个方向的加权滤波公式，其中设计的atan(1/2)方向上的加权滤波 公式为：

p_new(i,j)＝[T(-2,1)×p(i-2,j+1)+((T(-1,0)+T(-1,1))/2)×(p(i-1,j)+ p(i-1,j+1))/2+T(0,0)×p(i,j)+((T(1,0)+T(1,-1))/2) (p(i+1,j)+p(i+1,j-1))/2+T(2,-1)×p(i+2,j-1)]/K

K＝T(-2,1)+(T(-1,0)+T(-1,1))/2+T(0,0)+(T(1,0)+T(1,-1))/2+T(2,-1)

其中，p(i,j)为小目标图像中的像素点，则i,j为整数且i,j∈[1,256]；p_new(i,j) 为小目标图像该点处的新灰度值；T(i,j)为步骤(3)中所提到的5×5滤波矩阵， 且-2≤i,j≤2，atan(1/2)方向的加权滤波公式设计完毕。其他三个方向的加权处 理公式仅需将按照相应方向变化的p(i,j)和T(i,j)点坐标分别带入atan(1/2)方向 上的加权滤波公式即可。每个方向滤波点的选择方式参见图4。

(5)将所述4个方向的加权滤波公式分别处理小目标图像，并保存每个方向上 遍历小目标图像后的子图，参见图1。

5.1)将atan(1/2)方向的加权滤波公式作用于小目标图像；

5.2)判断atan(1/2)方向的加权滤波公式是否遍历完小目标图像，若没有，继 续进行步骤5.1)，若已遍历完，进行下一个步骤；

5.3)保存atan(1/2)方向上滤波完成的子图；

5.4)将π/2-atan(1/2)方向的加权滤波公式作用于小目标图像；

5.5)判断π/2-atan(1/2)方向的加权滤波公式是否遍历完小目标图像，若没有， 继续进行步骤5.4)，若已遍历完，进行下一个步骤；

5.6)保存π/2-atan(1/2)方向上滤波完成的子图；

5.7)将π/2+atan(1/2)方向的加权滤波公式作用于小目标图像；

5.8)判断π/2+atan(1/2)方向的加权滤波公式是否遍历完小目标图像，若没有， 继续进行步骤5.7)，若已遍历完，进行下一个步骤；

5.9)保存π/2+atan(1/2)方向上滤波完成的子图；

5.10)将π-atan(1/2)方向的加权滤波公式作用于小目标图像；

5.11)判断π-atan(1/2)方向的加权滤波公式是否遍历完小目标图像，若没有， 继续进行步骤5.10)，若已遍历完，进行下一个步骤；

5.12)保存π-atan(1/2)方向上滤波完成的子图。

(6)将上述得到的四副子图按照矩阵对应点相乘的方法融合成一副图像。

(7)归一化处理融合得到的图像，得到背景抑制后的小目标图像，其效果可参 见图5(b)、图6(b)、图7(b)、图8(b)，红外小目标识别结束。

本发明的技术效果还以通过仿真进一步说明：

实施例6

基于张量扩展的半脊方向滤波器红外小目标识别方法同实施例1-5，参见图 5,其中图5(a)是一个海天背景条件下的红外小目标图像，图5(b)是本发明方 法应用于图5(a)的处理结果图，图5(c)为剪切波变换算法处理图5(a)的处 理结果图。本发明的方法应用于图5(a)的时间消耗为3.7494秒，信噪比估值为 4.7034。而在相同仿真环境下，剪切波变换算法的时间消耗为4.2515秒，信噪比 估值为0.0162。对比两种方法的处理结果图可以明显看出，剪切波变换算法的处 理结果显示两目标点的位置时模糊不清楚，其中右边的点目标与周围的对比度过 低，几乎无法识别出。而本发明的处理结果清楚的显示出两个目标所在位置，成 功的进行了红外小目标的识别。所以在对图5(a)的处理效果上，本发明的方法 无论是在实时性、信噪比还是视觉表现上都明显好于剪切波变换算法。

实施例7

基于张量扩展的半脊方向滤波器红外小目标识别方法同实施例1-5，参见图 6,其中图6(a)是一个天空背景条件下的红外小目标图像，图6(b)是本发明方 法应用于图6(a)的处理结果图，图6(c)为剪切波变换算法处理图6(a)的处 理结果图。本发明的方法应用于图6(a)的时间消耗为3.7295秒，信噪比估值为 3.0308。而在相同仿真环境下，剪切波变换算法的时间消耗为6.0839秒，信噪比 估值为0.0082。对比两种方法的处理结果图可以明显看出，本发明的方法在处理 效果上虽然背景抑制的不全面，但本发明的方法很好的突出了小目标所在的位置， 视觉上要好于剪切波变换算法的处理效果。所以在对图6(a)的处理效果上，本 发明的方法无论是在实时性、信噪比还是视觉表现上都明显好于剪切波变换算法。

实施例8

基于张量扩展的半脊方向滤波器红外小目标识别方法同实施例1-5，参见图 7,其中图7(a)是一个天空背景条件下多目标的红外图像，图7(b)是本发明方 法应用于图7(a)的处理结果图，图7(c)为剪切波变换算法处理图7(a)的处 理结果图。本发明的方法应用于图7(a)的时间消耗为3.7355秒，信噪比估值为 0.0084。而在相同仿真环境下，剪切波变换算法的时间消耗为5.3031秒，信噪比 估值为0.0091。本发明设计的方向滤波公式，在兼顾处理效果的同时更多的采用 原图像内容，这样虽然会造成背景抑制的不全面，但同时也会大大降低漏检概率， 而在处理结果上也增强了小目标所在位置的灰度。虽然在对图7(a)的处理效果 上，本发明的方法在信噪比估值上不如剪切波变换算法，但是对比两种方法的处 理结果图可以明显看出，图7(C)的左上角的点目标与周围环境对比度低，而本 发明的方法处理效果视觉上好于剪切波变换算法的处理效果。所以在对图7(a) 的处理效果上，本发明的方法在实时性和处理结果的视觉表现上都明显好于剪切 波变换算法。

实施例9

基于张量扩展的半脊方向滤波器红外小目标识别方法同实施例1-5，参见图 8,其中图8(a)是一个天空背景条件下单目标的红外图像，图8(b)是本发明方 法应用于图8(a)的处理结果图，图8(c)为剪切波变换算法处理图8(a)的处 理结果图。观察图8(a)可以看出小目标所占像素个数少且与周围环境的灰度对 比度低。对比图8(b)和图8(c)可以看出，虽然图8(b)的背景抑制不全面， 但是由于本发明中采用大小为atan(1/2)，π/2-atan(1/2)，π/2+atan(1/2)， π-atan(1/2)的角度作为4个滤波方向，不同于普通方向滤波器的滤波方向，而在 每个滤波方向上选定5个以个像素距离为单位的加权滤波点后，使每次运算 包含更多的原图像信息，降低漏检率，所以可以有效识别出图8(a)中的小目标， 而剪切波变换算法不能有效识别出该图中的小目标。由此可以得出在处理图8(a) 中的弱小目标，本发明的方法明显优于剪切波变换算法。

综上所述，基于张量扩展的半脊方向滤波器红外小目标识别方法。具体包括以 下步骤：读取小目标图像；根据ridge滤波器构建的图像具有方向选择的特性，利 用余弦函数来构造本发明的图像模板；按照二维卷积规则将从图像模板上得到的 张量矩阵大小扩展为9×9，并选取该矩阵中间5×5个值为滤波矩阵；确定4个滤 波方向；设计每个方向上的加权滤波公式；待每个方向上的加权滤波公式遍历完 图像后，保存该方向上的子图；将保存的每个方向的子图按照矩阵对应点相乘的 方法融合成一幅图像，归一化处理后得到背景抑制后的小目标图像。本发明解决 了天空、海天背景下的红外小目标识别问题。本发明方法简单、便于硬件并行处 理，且漏检率低、信噪比高。适用于空天及海天背景下单个或多个小目标的红外 预警、红外制导导弹智能寻的。