一种红外弱小目标检测方法及检测系统与流程

本发明涉及图像处理技术领域，特别涉及一种红外弱小目标检测方法及检测系统。

背景技术：

随着红外探测技术的发展，复杂背景下的红外弱小目标检测技术在军事和民用领域中都越来越重要，如何能快速而准确的发现目标已经成为预警与跟踪系统的关键。但是由于红外弱小目标常常因为无形状或边缘信息、低信杂比成像以及复杂背景等因素的影响，大大增加了检测难度。

目前红外弱小目标的检测技术主要有基于图像序列的检测和基于单帧图像的检测两大类。基于图像序列的检测理论是通过序列图像中的相关性来估计目标，该类算法的一个重要前提是必须保证连续帧之间目标与背景的一致性，以及一些先验的目标信息，这些前提与先验信息限制了这类算法的应用，另外基于序列图像的检测方法通过多帧目标轨迹确认来实现目标检测，对于低信杂比成像的场景有较好的效果，但是相比单帧检测类算法而言其实时性较差，在实际工程应用中较少。

技术实现要素：

有鉴如此，有必要针对现有技术存在的缺陷，提供一种能够提高红外弱小目标检测准确性及实时性的红外弱小目标检测方法。

为实现上述目的，本发明采用下述技术方案：

一种红外弱小目标检测方法，包括下述步骤：

计算原始图像的facet方向导数特征；

在所述facet方向导数特征图的局部内，沿着当前方向计算相对极差对比度显著图；

将各个方向上的所述相对极差对比度显著图进行融合得到显著性图像；

提取所述显著性图像的目标。

在一些较佳的实施例中，在计算原图像的facet方向导数特征的步骤中，包括下述步骤：

基于facet模型，在原始图像的5x5邻域内灰度强度曲面拟合一个二元三次多项式f(r,c)的表达式，所述表达式如下，其中，r,c为5x5邻域内的行列坐标，ki是拟合系数：

0度方向、90度方向以及任意α度方向导数特征计算公式为：

其中ki可以由原图像i与卷积核wi的卷积快速计算，wi如下：

通过所述表达式计算原图像的facet方向导数特征。

在一些较佳的实施例中，在所述facet方向导数特征图的局部内，沿着当前方向计算相对极差对比度显著图的步骤中，具体包括：

设facet_α为沿α方向的一阶导数图像，facet_α(r,c)为位置中心点p(r,c)上的灰度值，dr为距离p(r,c)为r的邻域图像；

记沿着α方向且过p(r,c)的直线上的灰度值集合为line(dr,α)，并设frontline(dr,α)为沿α方向线且在点p(r,c)之前截取到的dr内灰度值集合，backline(dr,α)为沿α方向线且在点p(r,c)之后截取到的dr内灰度值集合；

分别计算集合frontline(dr,α)内的最大值maxvar(dr,α)和backline(dr,α)内的最小值minvar(dr,α)以及集合line(dr,α)内所有灰度的均值meanvar(dr,α)，则在中心点p(r,c)位置上的局部沿α方向的相对极差为：

定义显著性度量公式如下：

c(p,dr,α)＝exp(rr(p,dr,α))；

对局部沿α方向的相对极差公式进行近似为：

其中mean[·]表示求集合·的灰度均值；

设m(range,α)为由mean[frontline(dr,α)]-mean[backline(dr,α)]值构成的图像，则m(range,α)可以由下面的卷积公式求得：

m(range,α)＝facet_αek(range,α)

其中，k(range,α)只在沿α方向的前半部分为1，后半部分为-1，其余元素全为0，其余的以此类推；

m(mean,α)可以由下面的卷积公式求得：

m(mean,α)＝facet_αek(mean,α)，

k(range,α)只在沿α方向全为1，其余元素全为0，其余的以此类推；

将m(range,α)与m(mean,α)逐像素相除并做自然指数拉伸得到对应像素的显著度值。

在一些较佳的实施例中，在将各个方向上的所述相对极差对比度显著图进行融合得到显著性图像的步骤中，具体为：

采用下述公式将各个方向上的相对极差对比度显著图进行融合得到显著性图像，

在一些较佳的实施例中，在提取所述显著性图像的目标的步骤中，具体包括：

对所述显著性图像用高斯平滑滤波以及阈值分割提取目标。

在一些较佳的实施例中，所述阈值分割通过下述公式实现，

κ是限定因子，thr是分割阈值，t是分割图像，1值代表小目标所在像素的区域，0代表背景区域。

另外，本发明还提供了一种红外弱小目标检测系统，包括：特征提取模块，计算原始图像的facet方向导数特征；

显著度值获取模块，在所述facet方向导数特征图的局部内，沿着当前方向计算相对极差对比度显著图；

图像融合模块，将各个方向上的所述相对极差对比度显著图进行融合得到显著性图像；及

目标提取模块，提取所述显著性图像的目标。

本发明采用上述技术方案的优点是：

本发明提供的红外弱小目标检测方法及检测系统，通过计算原始图像的facet方向导数特征，在所述facet方向导数特征图的局部内，沿着当前方向计算相对极差对比度显著图，将各个方向上的所述相对极差对比度显著图进行融合得到显著性图像，再提取所述显著性图像的目标，本发明提供的红外弱小目标检测方法及检测系统，模型结构简单，计算复杂度低，并且可以基于卷积核近似计算，十分方便，具有很高的信杂比增益和背景抑制能力，弱小目标正确检测概率达到98％以上，实时性算法良好。

附图说明

图1为本发明实施例1提供的红外弱小目标检测方法的步骤流程图。

图2为本发明实施例1提供的天空复杂背景红外弱小目标实际检测图。

图3为本发明实施例1提供的海面复杂背景红外弱小目标实际检测图。

图4为本发明提供的一种红外弱小目标检测系统的结构示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。

实施例1

请参阅图1，为本发明实施例1提供的红外弱小目标检测方法，包括下述步骤：

步骤s110：计算原始图像的facet方向导数特征。

在一些较佳的实施例中，在计算原图像的facet方向导数特征的步骤中，包括下述步骤：

基于facet模型，在原始图像的5x5邻域内灰度强度曲面拟合一个二元三次多项式f(r,c)的表达式，所述表达式如下，其中，r,c为5x5邻域内的行列坐标，ki是拟合系数：

0度方向、90度方向以及任意α度方向导数特征计算公式为：

其中ki可以由原图像i与卷积核wi的卷积快速计算，wi如下：

通过所述表达式计算原图像的facet方向导数特征。

步骤s120：在所述facet方向导数特征图的局部内，沿着当前方向计算相对极差对比度显著图；

在一些较佳的实施例中，在所述facet方向导数特征图的局部内，沿着当前方向计算相对极差对比度显著图的步骤中，具体包括：

设facet_α为沿α方向的一阶导数图像，facet_α(r,c)为位置中心点p(r,c)上的灰度值，dr为距离p(r,c)为r的邻域图像；

记沿着α方向且过p(r,c)的直线上的灰度值集合为line(dr,α)，并设frontline(dr,α)为沿α方向线且在点p(r,c)之前截取到的dr内灰度值集合，backline(dr,α)为沿α方向线且在点p(r,c)之后截取到的dr内灰度值集合；

分别计算集合frontline(dr,α)内的最大值maxvar(dr,α)和backline(dr,α)内的最小值minvar(dr,α)以及集合line(dr,α)内所有灰度的均值meanvar(dr,α)，则在中心点p(r,c)位置上的局部沿α方向的相对极差为：

定义显著性度量公式如下：

c(p,dr,α)＝exp(rr(p,dr,α))；

对局部沿α方向的相对极差公式进行近似为：

其中mean[·]表示求集合·的灰度均值；

设m(range,α)为由mean[frontline(dr,α)]-mean[backline(dr,α)]值构成的图像，则m(range,α)可以由下面的卷积公式求得：

m(range,α)＝facet_αek(range,α)

其中，k(range,α)只在沿α方向的前半部分为1，后半部分为-1，其余元素全为0，其余的以此类推；

例如：7×7的α＝45°方向的，k(range,α)如下：

m(mean,α)可以由下面的卷积公式求得：

m(mean,α)＝facet_αek(mean,α)，

k(range,α)只在沿α方向全为1，其余元素全为0，其余的以此类推；

例如7×7的α＝45°方向的k(mean,α)为：

将m(range,α)与m(mean,α)逐像素相除并做自然指数拉伸得到对应像素的显著度值。

步骤s130：将各个方向上的所述相对极差对比度显著图进行融合得到显著性图像；

在将各个方向上的所述相对极差对比度显著图进行融合得到显著性图像的步骤中，具体为：

采用下述公式将各个方向上的相对极差对比度显著图进行融合得到显著性图像，

步骤s140：提取所述显著性图像的目标。

在提取所述显著性图像的目标的步骤中，具体包括：对所述显著性图像用高斯平滑滤波以及阈值分割提取目标。

具体地，所述阈值分割通过下述公式实现，

κ是限定因子，thr是分割阈值，t是分割图像，1值代表小目标所在像素的区域，0代表背景区域。

请参阅图2及图3，分别为本发明实施例1提供的天空复杂背景红外弱小目标实际检测图及海面复杂背景红外弱小目标实际检测图，通过上述红外弱小目标检测方法，其采用的模型简单，计算复杂度低，并且可以基于卷积核近似计算，十分方便，经过相关数据集测试，该方法可以具有很高的信杂比增益和背景抑制能力，弱小目标正确检测概率以及算法实时性优良。

本发明提供的红外弱小目标检测方法，通过计算原始图像的facet方向导数特征，在所述facet方向导数特征图的局部内，沿着当前方向计算相对极差对比度显著图，将各个方向上的所述相对极差对比度显著图进行融合得到显著性图像，再提取所述显著性图像的目标，本发明提供的红外弱小目标检测方法，模型结构简单，计算复杂度低，并且可以基于卷积核近似计算，十分方便，具有很高的信杂比增益和背景抑制能力，弱小目标正确检测概率达到98％以上，实时性算法良好。

实施例二

请参阅图4，为本发明提供的一种红外弱小目标检测系统，包括：特征提取模块110，计算原始图像的facet方向导数特征；显著度值获取模块120，在所述facet方向导数特征图的局部内，沿着当前方向计算相对极差对比度显著图；图像融合模块130，将各个方向上的所述相对极差对比度显著图进行融合得到显著性图像；及目标提取模块140，提取所述显著性图像的目标。

以下详细描述各个模块的工作方式。

在一些较佳的实施例中，特征提取模块110包括：基于facet模型，在原始图像的5x5邻域内灰度强度曲面拟合一个二元三次多项式f(r,c)的表达式，所述表达式如下，其中，r,c为5x5邻域内的行列坐标，ki是拟合系数：

0度方向、90度方向以及任意α度方向导数特征计算公式为：

其中ki可以由原图像i与卷积核wi的卷积快速计算，wi如下：

通过所述表达式计算原图像的facet方向导数特征。

显著度值获取模块120在所述facet方向导数特征图的局部内，沿着当前方向计算相对极差对比度显著图的步骤中，具体包括：

设facet_α为沿α方向的一阶导数图像，facet_α(r,c)为位置中心点p(r,c)上的灰度值，dr为距离p(r,c)为r的邻域图像；

记沿着α方向且过p(r,c)的直线上的灰度值集合为line(dr,α)，并设frontline(dr,α)为沿α方向线且在点p(r,c)之前截取到的dr内灰度值集合，backline(dr,α)为沿α方向线且在点p(r,c)之后截取到的dr内灰度值集合；

分别计算集合frontline(dr,α)内的最大值maxvar(dr,α)和backline(dr,α)内的最小值minvar(dr,α)以及集合line(dr,α)内所有灰度的均值meanvar(dr,α)，则在中心点p(r,c)位置上的局部沿α方向的相对极差为：

定义显著性度量公式如下：

c(p,dr,α)＝exp(rr(p,dr,α))；

对局部沿α方向的相对极差公式进行近似为：

其中mean[·]表示求集合·的灰度均值；

设m(range,α)为由mean[frontline(dr,α)]-mean[backline(dr,α)]值构成的图像，则m(range,α)可以由下面的卷积公式求得：

m(range,α)＝facet_αek(range,α)

其中，k(range,α)只在沿α方向的前半部分为1，后半部分为-1，其余元素全为0，其余的以此类推；

例如：7×7的α＝45°方向的，k(range,α)如下：

m(mean,α)可以由下面的卷积公式求得：

m(mean,α)＝facet_αek(mean,α)，

k(range,α)只在沿α方向全为1，其余元素全为0，其余的以此类推；

例如7×7的α＝45°方向的k(mean,α)为：

将m(range,α)与m(mean,α)逐像素相除并做自然指数拉伸得到对应像素的显著度值。

图像融合模块130将各个方向上的所述相对极差对比度显著图进行融合得到显著性图像，具体为：采用下述公式将各个方向上的相对极差对比度显著图进行融合得到显著性图像，

目标提取模块140对所述显著性图像用高斯平滑滤波以及阈值分割提取目标。

具体地，所述阈值分割通过下述公式实现，

κ是限定因子，thr是分割阈值，t是分割图像，1值代表小目标所在像素的区域，0代表背景区域。

请参阅图2及图3，分别为本发明实施例2提供的天空复杂背景红外弱小目标实际检测图及海面复杂背景红外弱小目标实际检测图，通过上述红外弱小目标检测方法，其采用的模型简单，计算复杂度低，并且可以基于卷积核近似计算，十分方便，经过相关数据集测试，该方法可以具有很高的信杂比增益和背景抑制能力，弱小目标正确检测概率以及算法实时性优良。

本发明提供的红外弱小目标检测系统，通过计算原始图像的facet方向导数特征，在所述facet方向导数特征图的局部内，沿着当前方向计算相对极差对比度显著图，将各个方向上的所述相对极差对比度显著图进行融合得到显著性图像，再提取所述显著性图像的目标，本发明提供的红外弱小目标检测系统，模型结构简单，计算复杂度低，并且可以基于卷积核近似计算，十分方便，具有很高的信杂比增益和背景抑制能力，弱小目标正确检测概率达到98％以上，实时性算法良好。

当然本发明的红外弱小目标检测方法还可具有多种变换及改型，并不局限于上述实施方式的具体结构。总之，本发明的保护范围应包括那些对于本领域普通技术人员来说显而易见的变换或替代以及改型。