CN104616299B - 一种基于空时偏微分方程的弱小目标检测方法 - Google Patents
一种基于空时偏微分方程的弱小目标检测方法 Download PDFInfo
- Publication number
- CN104616299B CN104616299B CN201510051773.7A CN201510051773A CN104616299B CN 104616299 B CN104616299 B CN 104616299B CN 201510051773 A CN201510051773 A CN 201510051773A CN 104616299 B CN104616299 B CN 104616299B
- Authority
- CN
- China
- Prior art keywords
- image
- background
- sky
- indicates
- partial differential
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
- G06T7/001—Industrial image inspection using an image reference approach
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10048—Infrared image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
Abstract
本发明公开了一种基于空时偏微分方程的弱小目标检测方法,该方法包括:(1)读入待处理图像序列的相邻三帧图像;(2)初始化参数w并计算图像的空时梯度(3)利用空时偏微分方程对图像进行背景预测;(4)对当前帧的图像和背景预测的结果进行差分,得到弱小目标的检测结果。本发明通过利用空时偏微分方程模型对远距离红外图像进行背景预测,使得算法能够对弱小目标进行检测,该方法为成像制导、安全监控、空间监测等方面提供了技术支撑,具有较强的实用价值。
Description
技术领域
本发明涉及一种基于空时偏微分方程的弱小目标检测方法,属于计算机视觉与航空航天技术领域。
背景技术
目前,弱小目标检测是成像制导、安全监控、空间监测技术中应值得重点研究的一个关键问题。随着现代科技的发展,人们希望在目标与成像传感器的距离很远的时候,即所谓的弱小目标阶段就能高可靠性的检测到目标。然而,由于自然气象条件(天气、风速、风向、气温、太阳辐射等)、目标背景环境(天空、地面、海面)以及目标结构等因素的影响,由于目标距离较远时,通常它的尺寸比较小,强度比较低。目标几何尺寸小,缺乏形状纹理信息,且目标强度弱,对比度小,也比较容易被噪声所淹没,因此,红外图象中的弱小目标检测成为成像目标检测与识别的一个难点。但是,鲁棒的弱小目标自动检测技术尚未完全突破,是当今世界高技术领域的热门研究课题,也是近年来信息处理技术中的研究热点之一。
目前,对弱小目标进行检测的方法主要可分为三大类:基于空域的检测方法、基于时域的检测。空域检测的方法包括:基于Butterworth高通滤波(BHPF)的检测方法(见文献L.Yang,J.Yang,K.Yang,“Adaptive detection for infrared small target under sea-sky complex background”,Electronics Letters,vol.40,no.17,pp.1083-1085,2004)和去最大-均值(Max-Mean)及最大-中值滤波(Max-Median)方法(见文献S.D.Deshpande,M.H.Er,V.Ronda,P.Chan,“Max-Mean and Max-Median filters for detection ofsmall-targets”,Proc.SPIE,vol.3809,pp.74-83,1999)。最大-中值滤波检测方法在处理单调背景的图像中表现良好,而Butterworth高通滤波在结构性杂波背景的弱小目标检测方面具有优势。除了空域处理之外,也有一些方法是基于时域的,如Caefer等人提出的基于小波变换和三项时域滤波(见文献C.E.Caefer,J.Silverman,J.M.Mooney,S.DiSalvo,R.W.Taylor,“Temporal filtering for point target detection in staring IRimagery:I.Damped sinusoid filters”,Proceedings of SPIE,vol.3373,pp.111-122,1998),该方法基于时域的滤波方法在目标亚像素级运动中,具有很好的表现,但是上述现有技术在复杂背景下弱小目标检测率低,而本发明能够很好地解决上面的问题。
发明内容
本发明目的在于提供了一种基于空时偏微分方程的弱小目标检测方法,该方法能够使弱小目标增强,解决了由于弱小目标几何尺寸小、强度弱、容易被噪声所淹没导致很难被检测的问题。
本发明解决其技术问题所采取的技术方案是:本发明提供了一种基于空时偏微分方程进行复杂背景预测,该方法能在有效对复杂背景中的结构性纹理进行有效的预测,能够为成像制导、安全监控、空间监测技术等方面的应用提供技术支撑。
方法流程:
步骤1:读取图像序列In,n∈[1,N],其中,n为图像序列的帧号,N为图像序列的总帧数。为了更好地进行预测,本发明将每连续的三帧图像组合在一起视为一个三维图像。每个像素点的坐标为(x,y,z)(x,y表示图像像素空间位置,z表示时间维上的位置)。因此,本发明对红外弱小目标场景图像模型进行建模为
In(x,y,z)=Bn(x,y,z)+Tn(x,y,z)+Nn(x,y,z) (1)
In表示红外图像序列中的第n帧图像;Bn和Tn分别表示相应背景图像和目标信号;Nn表示图像的噪声。背景预测就是根据原图对背景图像进行估计,估计值用表示,
其中H表示预测算法,η(x,y,z)表示像素(x,y,z)的一个特定邻域,为背景预测结果。
步骤2:构建用于背景预测的空时偏微分方程模型。假设背景预测用表示,则
其中,uk(x,y,z)表示由原图像I(x,y,z)进行迭代预测的最终结果,k为迭代次数。根据式(1),构建偏微分方程的能量函数为
其中,w为一个正参数。ux,uy,uz分别是在x、y和z方向上的一阶导数。
步骤3:根据变分原理,最小化式(4)的能量函数可通过求解Euler-Lagrange公式
其中,F=F(x,y,z,u,ux,uy,uz)。可求得模型的解为:
即认为uk为估计的背景
步骤4:通过差分将估计的背景从原始图像序列I中消除,得到残差图像,即:
残差图像f中通常包含残差背景、目标及噪声。然而由于大部分的背景能量被减去,残差图像中的目标与背景的相对差值和信杂比将增强,因此可以有效的检测出目标。
本发明具体的实现过程包括:
1.根据读取图像序列In,n∈[1,N],其中,n为图像序列的帧号,N为图像序列的总帧数。为了更好地进行预测,本发明将每连续的三帧图像组合在一起视为一个三维图像。每个像素点的坐标为(x,y,z)(x,y表示图像像素空间位置,z表示时间维上的位置)。因此,本发明对红外弱小目标场景图像模型进行建模为
In(x,y,z)=Bn(x,y,z)+Tn(x,y,z)+Nn(x,y,z) (8)
In表示红外图像序列中的第n帧图像;Bn和Tn分别表示相应背景图像和目标信号;Nn表示图像的噪声。背景预测就是根据原图对背景图像进行估计,估计值用表示,
其中H表示预测算法,η(x,y,z)表示像素(x,y,z)的一个特定邻域,B为背景预测结果。
2.假设背景预测用表示:
其中,uk(x,y,z)表示由原图像I(x,y,z)进行迭代预测的最终结果,k为迭代次数。根据式(8),构建用于背景预测的空时偏微分方程模型为
其中,w为一个正参数。ux,uy,uz分别是在x、y和z方向上的一阶导数。
3.根据变分原理,最小化式(11)的能量函数可通过求解Euler-Lagrange公式
其中,F=F(x,y,z,u,ux,uy,uz),uxx,uyy,uzz分别为二阶偏导数,uxy,uxz,uyz分别为混合二阶偏导。可求得模型的解为:
即认为uk为估计的背景
有益效果:
1、本发明应用于红外弱小目标检测,能够检测出几何尺寸小、强度弱的远距离红外弱小目标。
2、本发明能够为成像制导、安全监控、空间监测技术等方面的应用提供技术支撑,具有很强的使用价值。
附图说明
图1为本发明的方法流程图。
图2为当前帧像素的空时邻域像素示意图。
图3为原始红外图像。
图4为原图三维图。
图5为空时偏微分方程背景预测图像。
图6为弱小目标检测结果示意图。
图7为弱小目标检测三维图。
具体实施方式
下面结合附图和实例对本发明作进一步详细的说明。
如图1所示,本发明一种基于空时偏微分方程的弱小目标检测方法,该方法主要通过构建一种空时偏微分方程模型对背景杂波进行预测及抑制,进而进行红外弱小目标检测。具体包括如下步骤:
步骤1:读取图像序列I=(I1,I2,I3...In-1,In,In+1...IN-2,IN-1,IN),其中,n为图像序列的帧号,N为图像序列的总帧数。为了更好地进行预测,本发明将每连续的三帧图像组合在一起视为一个三维图像,如图2所示。其像素点的坐标为(x,y,z)(为与迭代时间t区分,用时间维上的变量用z表示)。因此,对红外弱小目标场景图像模型进行建模为
In(x,y,z)=Bn(x,y,z)+Tn(x,y,z)+Nn(x,y,z) (14)
其中(x,y,z)表示图像序列中像素位置;In表示红外图像序列中的第n帧图像;Bn和Tn分别表示相应背景图像和目标信号;Nn表示噪声。
背景预测就是根据原图对背景图像进行估计,估计值用表示,
其中H表示预测算法,η(x,y,z)表示像素(x,y,z)的一个特定邻域,为背景预测结果。
步骤2:构建用于背景预测的偏微分方程模型:背景预测就是根据原图像对背景图像进行估计,假设估计值用表示,则
其中,uk(x,y,z)表示由原图像I(x,y,z)进行迭代预测的最终结果,k为迭代次数。根据式(14),构建偏微分方程的能量函数为
其中,w为一个正参数。ux,uy,uz分别是在x、y和z向上的一阶导数。在求解过程中,其离散化形式分别为
ux=(u(x+1,y,z)+u(x-1,y,z))/2
uy=(u(x,y+1,z)+u(x,y-1,z))/2 (18)
uz=(u(x,y,z+1)+u(x,y,z-1))/2
步骤3:模型求解:根据变分原理,最小化式(17)的能量函数可通过求解Euler-Lagrange公式
其中,F=F(x,y,z,u,ux,uy,uz),且
uxx,uyy,uzz分别为二阶偏导数,uxy,uxz,uyz分别为混合二阶偏导,其离散形式分别为
uxx=u(x+1,y,z)+u(x-1,y,z)-2u(x,y,z)
uyy=u(x,y+1,z)+u(x,y-1,z)-2u(x,y,z)
uzz=u(x,y,z+1)+u(x,y,z-1)-2u(x,y,z)
最后,求得模型的解为:
即认为uk为估计的背景
步骤4:通过差分将估计的背景从原始图像序列I中消除,得到残差图像,即
可以看出,残差图像f中通常包含残差背景、目标及噪声。然而由于大部分的背景能量被减去,残差图像中的目标与背景的相对差值和信杂比将增强,因此可以有效的检测出目标。
本发明应用于红外弱小目标检测,其技术效果体现在:能够检测出几何尺寸小、强度弱的远距离红外弱小目标。本发明能够为成像制导、安全监控、空间监测技术等方面的应用提供技术支撑,具有很强的使用价值。
为了验证本发明红外小目标增强的有效性,对实拍红外小目标图像序列进行试验。图3给出了红外小目标图像序列中的第10帧原始红外图像,图4为原始红外图像的三维图。利用上述处理步骤,对图3进行空时偏微分方程进行背景预测得到预测结果图4,弱小目标检测结果和检测结果三维示意图如图5和图6所示。为更直观地对比各算法的弱小目标检测性能,使用目标信杂比(SCR)及其增益(SCR Gain)作为性能指标,可以综合衡量算法的目标检测性能和目标增强效果。
信杂比SCR定义为
其中,μs是目标亮度,σcn是邻域亮度的标准差。信杂比增益(单位dB)定义为
其中SCRin为图像序列的原始信杂比,SCRout为处理结果序列的信杂比。选择最大-均值(Max-mean)、最大-中值滤波(Max-median)及Butterworth高通滤波(BHPF)与本发明的方法进行对比。
表1:为与其他方法检测结果比较表
Max-mean | Max-median | BHPF | 本发明的方法 | |
SCR Gain | 3.8596 | 4.5199 | 4.6389 | 8.4130 |
BSF | 8.1219 | 8.6468 | 8.8528 | 21.7679 |
表1为用各种算法进行红外弱小目标检测试验得到的的信杂比和信杂比增益的结果,可以明显地看出本发明的算法要优于其他几种方法。
以上所述仅为本发明的一个实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于空时偏微分方程的弱小目标检测方法,其特征在于:所述方法是通过利用空时偏微分方程进行背景预测,包括如下步骤:
步骤1:读取图像序列In,n∈[1,N],其中,n为图像序列的帧号,N为图像序列的总帧数,将每连续的三帧图像组合在一起视为一个三维图像,每个像素点的坐标为(x,y,z),对红外弱小目标场景图像模型进行建模为:
In(x,y,z)=Bn(x,y,z)+Tn(x,y,z)+Nn(x,y,z) (1)
In表示红外图像序列中的第n帧图像;Bn和Tn分别表示相应背景图像和目标信号;Nn表示图像的噪声;背景预测就是根据原图对背景图像进行估计,估计值用表示,
其中H表示预测算法,η(x,y,z)表示像素(x,y,z)的一个特定邻域,为背景预测结果;
步骤2:构建用于背景预测的空时偏微分方程模型;假设背景预测用表示,则
其中,uk(x,y,z)表示由原图像I(x,y,z)进行迭代预测的最终结果,k为迭代次数,根据式(1),构建偏微分方程的能量函数为
其中,w为一个正参数,ux,uy,uz分别是在x、y和z方向上的一阶导数;
步骤3:根据变分原理,最小化式(4)的能量函数通过求解Euler-Lagrange公式:
其中,F=F(x,y,z,u,ux,uy,uz),能求得模型的解为:
即认为uk为估计的背景
步骤4:通过差分将估计的背景从原始图像序列I中消除,得到残差图像,即
残差图像f中包含残差背景、目标及噪声。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510051773.7A CN104616299B (zh) | 2015-01-30 | 2015-01-30 | 一种基于空时偏微分方程的弱小目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510051773.7A CN104616299B (zh) | 2015-01-30 | 2015-01-30 | 一种基于空时偏微分方程的弱小目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104616299A CN104616299A (zh) | 2015-05-13 |
CN104616299B true CN104616299B (zh) | 2019-02-19 |
Family
ID=53150732
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510051773.7A Active CN104616299B (zh) | 2015-01-30 | 2015-01-30 | 一种基于空时偏微分方程的弱小目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104616299B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104766334B (zh) * | 2015-04-21 | 2017-12-29 | 西安电子科技大学 | 红外弱小目标检测跟踪方法及其装置 |
CN108765460B (zh) * | 2018-05-22 | 2021-03-30 | 北京航空航天大学 | 基于高光谱图像的空时联合异常检测方法和电子设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1606033A (zh) * | 2004-11-18 | 2005-04-13 | 上海交通大学 | 红外图像序列中弱小目标的检测和跟踪方法 |
JP2011180717A (ja) * | 2010-02-26 | 2011-09-15 | Nagoya Institute Of Technology | 画像処理用のTotalVariationフィルタおよび画像処理プログラム |
CN102542540A (zh) * | 2011-12-31 | 2012-07-04 | 华中科技大学 | 一种基于pde的红外图像背景抑制方法 |
CN103700113A (zh) * | 2012-09-27 | 2014-04-02 | 中国航天科工集团第二研究院二O七所 | 一种下视复杂背景弱小运动目标检测方法 |
JP2014134442A (ja) * | 2013-01-10 | 2014-07-24 | Mitsubishi Electric Corp | 赤外線目標検出装置 |
CN104299229A (zh) * | 2014-09-23 | 2015-01-21 | 西安电子科技大学 | 一种基于时空域背景抑制的红外弱小目标检测方法 |
-
2015
- 2015-01-30 CN CN201510051773.7A patent/CN104616299B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1606033A (zh) * | 2004-11-18 | 2005-04-13 | 上海交通大学 | 红外图像序列中弱小目标的检测和跟踪方法 |
JP2011180717A (ja) * | 2010-02-26 | 2011-09-15 | Nagoya Institute Of Technology | 画像処理用のTotalVariationフィルタおよび画像処理プログラム |
CN102542540A (zh) * | 2011-12-31 | 2012-07-04 | 华中科技大学 | 一种基于pde的红外图像背景抑制方法 |
CN103700113A (zh) * | 2012-09-27 | 2014-04-02 | 中国航天科工集团第二研究院二O七所 | 一种下视复杂背景弱小运动目标检测方法 |
JP2014134442A (ja) * | 2013-01-10 | 2014-07-24 | Mitsubishi Electric Corp | 赤外線目標検出装置 |
CN104299229A (zh) * | 2014-09-23 | 2015-01-21 | 西安电子科技大学 | 一种基于时空域背景抑制的红外弱小目标检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104616299A (zh) | 2015-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108009515B (zh) | 一种基于fcn的无人机航拍图像的输电线定位识别方法 | |
CN105574855B (zh) | 云背景下基于模板滤波和虚警抑制的红外小目标检测方法 | |
CN106469313B (zh) | 一种管径自适应时空域滤波的弱小目标检测方法 | |
CN109003277A (zh) | 一种复杂背景下红外弱小目标检测方法及装置 | |
CN111027496B (zh) | 基于空时联合局部对比度的红外弱小目标检测方法 | |
CN103761731A (zh) | 一种基于非下采样轮廓波变换的红外空中小目标检测方法 | |
CN106530271B (zh) | 一种红外图像显著性检测方法 | |
CN104657945A (zh) | 复杂背景下多尺度时空联合滤波的红外小目标检测方法 | |
CN102819740B (zh) | 一种单帧红外图像弱小目标检测和定位方法 | |
CN107392095A (zh) | 一种基于掩码图像的红外弱小目标检测算法 | |
CN107507209B (zh) | 极化sar图像的素描图提取方法 | |
CN104899866A (zh) | 一种智能化的红外小目标检测方法 | |
CN110147778B (zh) | 稀土矿开采识别方法、装置、设备及存储介质 | |
CN105319537A (zh) | 基于空间相关性的航海雷达同频干扰抑制方法 | |
CN104268877A (zh) | 一种红外图像海天线自适应检测方法 | |
CN109711256B (zh) | 一种低空复杂背景无人机目标检测方法 | |
CN106156758B (zh) | 一种sar海岸图像中海岸线提取方法 | |
CN108648184A (zh) | 一种遥感图像高空卷云的检测方法 | |
CN110263646A (zh) | 一种基于卷积神经网络的海面弱目标探测方法和系统 | |
CN110889327A (zh) | 一种基于热红外图像的水域周边排污口智能检测方法 | |
CN105469428B (zh) | 一种基于形态学滤波和svd的弱小目标检测方法 | |
CN104616299B (zh) | 一种基于空时偏微分方程的弱小目标检测方法 | |
CN112669332B (zh) | 一种基于双向局部极大值和峰值局部奇异性判断海天条件和检测红外目标的方法 | |
CN111353496B (zh) | 一种红外弱小目标实时检测方法 | |
CN106570889A (zh) | 一种检测红外视频中弱小目标的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |