CN102819740B - 一种单帧红外图像弱小目标检测和定位方法 - Google Patents
一种单帧红外图像弱小目标检测和定位方法 Download PDFInfo
- Publication number
- CN102819740B CN102819740B CN201210249215.8A CN201210249215A CN102819740B CN 102819740 B CN102819740 B CN 102819740B CN 201210249215 A CN201210249215 A CN 201210249215A CN 102819740 B CN102819740 B CN 102819740B
- Authority
- CN
- China
- Prior art keywords
- candidate target
- image
- dog
- target
- pixel
- 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
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种基于形态学带通滤波和尺度空间理论的单帧红外图像弱小目标检测和定位方法。基本思想是:首先对预处理后的红外图像进行形态学带通滤波,检测出红外弱小目标的潜在区域;然后利用DoG算子对预处理后的红外图像进行多尺度空间变换,并在目标的潜在区域通过尺度空间的极大值检测获得候选目标的位置和尺度;最后通过对候选目标的信杂比进行阈值化实现红外弱小目标的精确检测和定位。利用本发明方法,可以从复杂背景下的单帧红外图像中检测并定位出弱小目标。本发明具有很高的检测精度和较低的虚警率,同时具有较小的运算量,利于硬件***实现。
Description
技术领域
本发明涉及一种单帧红外图像弱小目标检测和定位方法,可以应用于复杂背景下的单帧红外图像弱小目标检测和定位。
背景技术
复杂背景下单帧红外图像的弱小目标检测和定位是红外搜索与跟踪***等领域的一项关键技术,也一直是该领域的研究热点和难点。受红外传感器技术、大气辐射、作用距离和探测器噪声等因素的影响,红外图像的对比度和信噪比都很低,加之远距离目标的成像面积较小,目标的成像面积仅为几个像素到几十个像素,在视场中表现为孤点或斑点,基本上无形状和纹理特征可以利用。因而,传统的基于目标大小、形状及特征的检测方法无法得到良好的检测和定位效果。
目前已有的红外图像弱小目标检测方法,概括起来可以分为基于时间域的检测方法、基于空间域的检测方法和基于变换域的检测方法三大类。基于时间域的检测方法主要包括图像熵差法、基于序列图像检测法和基于面模型的匹配滤波法等。但是,基于时间域的检测方法需要用到不同时间采集的多帧图像,且图像序列必须配准,不适用于单帧图像的目标检测。基于空间域的检测方法主要利用小目标亮度高于其邻域背景,且与背景不相关的特性,直接在空间域对图像进行处理,其代表方法主要包括中值滤波法、二维最小均方误差法和偏微分方程法。然而,在低信噪比且具有复杂背景的红外弱小目标图像中,目标和背景的空间灰度分布非常接近,导致传统的基于空间域的检测方法难以准确地区分目标和背景,从而无法取得较好的检测效果。基于变换域的检测方法主要包括频域高通滤波法和基于小波变换的检测方法,该类方法对图像数据具有更为精细的分析能力,但存在运算复杂的缺点,不利于硬件***实现。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种单帧红外图像弱小目标检测和定位方法,基于形态学带通滤波和尺度空间理论的单帧红外图像弱小目标检测和定位方法,可以自动地从低信噪比且具有复杂背景的单帧红外图像中检测并定位出弱小目标,具有很高的检测精度和较低的虚警率。
技术方案
一种单帧红外图像弱小目标检测和定位方法,其特征在于步骤如下:
步骤1预处理:采用加权平均法对输入红外图像的RGB三个分量进行加权平均得到灰度图像,然后利用灰度线性变换函数将灰度图像的灰度范围映射到[Imin,Imax]的灰度区间,得到预处理后的图像I;
其中:加权平均法计算公式为f(x,y)=0.3R(x,y)+0.59G(x,y)+0.11B(x,y),f(x,y)为加权平均法得到的灰度图像在像素点(x,y)的灰度值,R(x,y)、G(x,y)和B(x,y)分别为输入红外图像在像素点(x,y)的RGB三个分量值;灰度线性变换函数为 I(x,y)为f(x,y)经过灰度线性变换后的灰度值,fmin、fmax分别为灰度线性变换前的灰度图像的最小灰度值和最大灰度值;
步骤2形态学带通滤波检测目标潜在区域:采用大小为S1×S1的正方形结构元素对图像I进行开运算得到图像I1,采用大小为S2×S2的正方形结构元素对图像I进行开运算得到图像I2,计算I2和I1的差分图像ID=I2-I1,然后对差分图像ID采用阈值Th1进行二值化,得到二值化图像 IBW中像素值为1的像素点的集合为目标的潜在区域其中,S1是取值范围为6~12之间的整数;S2是取值范围为1~2之间的整数;阈值Th1的取值范围为0.15~0.30;
步骤3生成DoG多尺度空间:计算图像I与不同尺度因子的DoG算子的卷积得到DoG多尺度空间{Dσ|σ=σ1,...,σN};其中:σ为尺度因子,σn=σ0·kn-1,σ0是一个取值范围为0.5~0.65之间的常数,k是一个取值范围为1.15~1.25之间的常数;下标和上标n表示尺度空间的标号,1≤n≤N,N为DoG多尺度空间的总层数,N的取值范围为10~12之间的整数;Dσ表示尺度因子为σ的DoG算子对应的DoG空间;DoG(x,y,σ)表示尺度因子为σ的DoG算子,DoG(x,y,σ)=G(x,y,σ)-G(x,y,kσ);G(x,y,σ)为高斯核函数,G(x,y,σ)=(1/2πσ2)exp(-(x2+y2)/2σ2);*表示卷积运算;
步骤4候选目标检测和定位:将中像素点对应的像素值分别和它的8个相邻点、上一尺度中的9个相邻点以及下一尺度中的9个相邻点共26个像素点的值进行比较,如果为极大值,则像素点为候选目标坐标点,σm为候选目标尺度,且候选目标大小d和候选目标尺度σm的关系为由此可以得到一些列的候选目标其中:下标m表示尺度空间的标号,2≤m≤N-1;表示第j个候选目标,它的坐标为尺度为下标j表示候选目标的标号,1≤j≤J,J表示候选目标的总个数;
步骤5目标验证:计算每个候选目标的局部信杂比并通过对局部信杂比进行阈值化得到最终的目标检测结果;其中:表示尺度因子为的DoG算子对应的DoG空间中像素点的像素值;σj为第j个候选目标的背景区域的灰度标准差, μj为第j个候选目标的背景区域的灰度均值,Cj表示第j个候选目标的背景区域,Mj表示Cj中的元素个数,Aj表示第j个候选目标区域, Bj表示第j个候选目标和背景的总区域, 阈值Th2的取值范围为3~5。
有益效果
本发明提出的一种单帧红外图像弱小目标检测和定位方法,是一种基于形态学带通滤波和尺度空间理论的单帧红外图像弱小目标检测和定位方法,首先对预处理后的红外图像进行形态学带通滤波,检测出红外弱小目标的潜在区域;然后利用DoG算子对预处理后的红外图像进行多尺度空间变换,并在目标的潜在区域通过尺度空间的极大值检测获得候选目标的位置和尺度;最后通过对候选目标的信杂比进行阈值化实现红外弱小目标的精确检测和定位。本发明方法可以自动地从低信噪比且具有复杂背景的单帧红外图像中检测并定位出弱小目标。该方法具有很高的检测精度和较低的虚警率,同时具有较小的运算量,利于硬件***实现。
附图说明
图1:本发明方法的基本流程图
图2:候选目标检测和定位示意图
(a)输入图像
(b)DoG多尺度空间
(c)DoG空间极大值检测
(d)检测与定位结果
图3:本发明方法的部分检测结果
(a)原始红外图像(白色方框代表实际的弱小目标)
(b)原始红外图像的三维显示(箭头指示弱小目标)
(c)形态学带通滤波检测目标潜在区域
(d)DoG尺度空间的其中一层(σ=2.5799)
(e)候选目标检测和定位(蓝色“+”代表候选目标)
(f)最终的目标检测结果和定位(红色“+”代表最终的检测目标)
具体实施方式
现结合实施例、附图对本发明作进一步描述:
用于实施的硬件环境是:IntelPentium2.13GHzCPU计算机、2.0GB内存,运行的软件环境是:MatlabR2010b和WindowsXP。选取了100幅复杂背景下的红外图像进行红外弱小目标检测实验,其中,天空背景27幅(包含飞机类弱小目标),海空背景34幅(包含船只类弱小目标),地面背景39幅(包含车辆类弱小目标),每幅图像至少有1个弱小目标,100幅图像共包含268个弱小目标。
本发明具体实施如下:
1、预处理:采用加权平均法对每一幅输入红外图像的RGB三个分量进行加权平均得到灰度图像,然后利用灰度线性变换函数进行图像的灰度调整,使红外图像的灰度范围均处于灰度区间[0,1]内,得到得到预处理后的图像I;
所采用的加权平均法计算公式为:f(x,y)=0.3R(x,y)+0.59G(x,y)+0.11B(x,y);
所采用的灰度线性变换函数为:
其中,f(x,y)为加权平均法得到的灰度图像在像素点(x,y)的灰度值,R(x,y)、G(x,y)和B(x,y)分别为输入红外图像在像素点(x,y)的RGB三个分量值,I(x,y)为f(x,y)经过灰度线性变换后的灰度值,fmin、fmax分别为灰度线性变换前的灰度图像的最小灰度值和最大灰度值。
2、形态学带通滤波检测目标潜在区域:采用大小为10×10的正方形结构元素对图像I进行开运算得到图像I1,采用大小为2×2的正方形结构元素对图像I进行开运算得到图像I2,计算I2和I1的差分图像ID=I2-I1,然后对差分图像ID采用阈值0.2进行二值化,得到二值化图像 IBW中像素值为1的像素点的集合为目标的潜在区域
3、生成DoG多尺度空间:计算图像I与不同尺度因子的DoG算子的卷积得到DoG多尺度空间{Dσ|σ=σ1,...,σN};其中:σ为尺度因子,σn=σ0·kn-1,σ0=0.6,k=1.2;下标和上标n表示尺度空间的标号,1≤n≤N,N为DoG多尺度空间的总层数,N=11;Dσ表示尺度因子为σ的DoG算子对应的DoG空间;DoG(x,y,σ)表示尺度因子为σ的DoG算子,DoG(x,y,σ)=G(x,y,σ)-G(x,y,kσ);G(x,y,σ)为高斯核函数,G(x,y,σ)=(1/2πσ2)exp(-(x2+y2)/2σ2);*表示卷积运算。
4、候选目标检测和定位:将中像素点对应的像素值分别和它的8个相邻点、上一尺度中的9个相邻点以及下一尺度中的9个相邻点共26个像素点的值进行比较,如果为极大值,则像素点为候选目标坐标点,σm为候选目标尺度,且候选目标大小d和候选目标尺度σm的关系为由此可以得到一些列的候选目标其中:下标m表示尺度空间的标号,2≤m≤10;表示第j个候选目标,它的坐标为尺度为下标j表示候选目标的标号,1≤j≤J,J表示候选目标的总个数。
5、目标验证:设定阈值Th2=4,计算每个候选目标的局部信杂比并通过对局部信杂比进行阈值化得到最终的目标检测结果;其中:表示尺度因子为的DoG算子对应的DoG空间中像素点的像素值;σj为第j个候选目标的背景区域的灰度标准差, μj为第j个候选目标的背景区域的灰度均值,Cj表示第j个候选目标的背景区域,Mj表示Cj中的元素个数,Aj表示第j个候选目标区域, Bj表示第j个候选目标和背景的总区域,
选用正确检测率、虚警率和每幅图像的平均运行时间对本发明的有效性进行评估。其中:正确检测率定义为正确检测的弱小目标个数与总的弱小目标个数之比;虚警率定义为虚警个数与正确检测的弱小目标个数和虚警个数之和的比值。同时,将本发明所得的检测结果与基于LoG尺度空间的小目标检测算法、基于图像稀疏表示的小目标检测算法进行了对比,对比结果如表1所示。正确检测率、虚警率和每幅图像的平均运行时间三个评价指标均表明了本发明方法的有效性。
表1检测结果评价
Claims (1)
1.一种单帧红外图像弱小目标检测和定位方法,其特征在于步骤如下:
步骤1预处理:采用加权平均法对输入红外图像的RGB三个分量进行加权平均得到灰度图像,然后利用灰度线性变换函数将灰度图像的灰度范围映射到[Imin,Imax]的灰度区间,得到预处理后的图像I;
其中:加权平均法计算公式为f(x,y)=0.3R(x,y)+0.59G(x,y)+0.11B(x,y),f(x,y)为加权平均法得到的灰度图像在像素点(x,y)的灰度值,R(x,y)、G(x,y)和B(x,y)分别为输入红外图像在像素点(x,y)的RGB三个分量值;灰度线性变换函数为I(x,y)为f(x,y)经过灰度线性变换后的灰度值,fmin、fmax分别为灰度线性变换前的灰度图像的最小灰度值和最大灰度值;
步骤2形态学带通滤波检测目标潜在区域:采用大小为S1×S1的正方形结构元素对图像I进行开运算得到图像I1,采用大小为S2×S2的正方形结构元素对图像I进行开运算得到图像I2,计算I2和I1的差分图像ID=I2-I1,然后对差分图像ID采用阈值Th1进行二值化,得到二值化图像IBW中像素值为1的像素点的集合为目标的潜在区域其中,S1是取值范围为6~12之间的整数;S2是取值范围为1~2之间的整数;阈值Th1的取值范围为0.15~0.30;
步骤3生成DoG多尺度空间:计算图像I与不同尺度因子的DoG算子的卷积Dσn=DoG(x,y,σn)*I,得到DoG多尺度空间{Dσ|σ=σ1,...,σN};其中:σ为尺度因子,σn=σ0·kn-1,σ0是一个取值范围为0.5~0.65之间的常数,k是一个取值范围为1.15~1.25之间的常数;下标和上标n表示尺度空间的标号,1≤n≤N,N为DoG多尺度空间的总层数,N的取值范围为10~12之间的整数;Dσ表示尺度因子为σ的DoG算子对应的DoG空间;DoG(x,y,σ)表示尺度因子为σ的DoG算子,DoG(x,y,σ)=G(x,y,σ)-G(x,y,kσ);G(x,y,σ)为高斯核函数,G(x,y,σ)=(1/2πσ2)exp(-(x2+y2)/2σ2);*表示卷积运算;
步骤4候选目标检测和定位:将中像素点对应的像素值分别和它的8个相邻点、上一尺度中的9个相邻点以及下一尺度中的9个相邻点共26个像素点的值进行比较,如果为极大值,则像素点为候选目标坐标点,σm为候选目标尺度,且候选目标大小d和候选目标尺度σm的关系为由此可以得到一些列的候选目标其中:下标m表示尺度空间的标号,2≤m≤N-1;表示第j个候选目标,它的坐标为尺度为下标j表示候选目标的标号,1≤j≤J,J表示候选目标的总个数;
步骤5目标验证:计算每个候选目标的局部信杂比并通过对局部信杂比进行阈值化得到最终的目标检测结果;其中:表示尺度因子为的DoG算子对应的DoG空间中像素点的像素值;σj为第j个候选目标的背景区域的灰度标准差,μj为第j个候选目标的背景区域的灰度均值,Cj表示第j个候选目标的背景区域,Mj表示Cj中的元素个数,Aj表示第j个候选目标区域,Bj表示第j个候选目标和背景的总区域,
阈值Th2的取值范围为3~5;
所述的灰度区间[Imin,Imax]为[0,1]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210249215.8A CN102819740B (zh) | 2012-07-18 | 2012-07-18 | 一种单帧红外图像弱小目标检测和定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210249215.8A CN102819740B (zh) | 2012-07-18 | 2012-07-18 | 一种单帧红外图像弱小目标检测和定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102819740A CN102819740A (zh) | 2012-12-12 |
CN102819740B true CN102819740B (zh) | 2016-04-20 |
Family
ID=47303849
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210249215.8A Active CN102819740B (zh) | 2012-07-18 | 2012-07-18 | 一种单帧红外图像弱小目标检测和定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102819740B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103530629B (zh) * | 2013-11-01 | 2016-08-31 | 四川九洲电器集团有限责任公司 | 红外小信号目标增程检测方法及*** |
CN104268844B (zh) * | 2014-10-17 | 2017-01-25 | 中国科学院武汉物理与数学研究所 | 一种基于加权局部图像熵的小目标红外图像处理方法 |
CN104599273B (zh) * | 2015-01-22 | 2017-07-28 | 南京理工大学 | 基于小波多尺度交叉运算的海空背景红外小目标检测方法 |
CN104766079B (zh) * | 2015-05-05 | 2018-12-07 | 四川九洲电器集团有限责任公司 | 一种远距离红外弱小目标检测方法 |
CN104834915B (zh) * | 2015-05-15 | 2017-12-19 | 中国科学院武汉物理与数学研究所 | 一种复杂云天背景下小红外目标检测方法 |
CN108470342A (zh) * | 2017-02-23 | 2018-08-31 | 南宁市富久信息技术有限公司 | 一种红外图像边缘检测方法 |
CN109272489B (zh) * | 2018-08-21 | 2022-03-29 | 西安电子科技大学 | 基于背景抑制与多尺度局部熵的红外弱小目标检测方法 |
CN109523575A (zh) * | 2018-11-12 | 2019-03-26 | 南通理工学院 | 红外弱小目标检测方法 |
CN109741396A (zh) * | 2018-12-12 | 2019-05-10 | 天津津航技术物理研究所 | 一种极弱小红外目标检测方法 |
CN111368585B (zh) * | 2018-12-25 | 2023-04-21 | 中国科学院长春光学精密机械与物理研究所 | 一种弱小目标检测方法、检测***、存储装置及终端设备 |
CN110223344B (zh) * | 2019-06-03 | 2023-09-29 | 哈尔滨工程大学 | 一种基于形态学和视觉注意机制的红外小目标检测方法 |
CN111783656B (zh) * | 2020-06-30 | 2024-03-08 | 哈尔滨工程大学 | 一种基于统计特性预分割的自适应红外小目标检测方法 |
CN113034533B (zh) * | 2021-04-06 | 2022-05-20 | 电子科技大学 | 一种基于空时平稳性的红外小目标检测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1251144C (zh) * | 2002-09-12 | 2006-04-12 | 上海交通大学 | 复杂背景条件下红外弱小目标的检测与跟踪方法 |
CN1207687C (zh) * | 2003-07-03 | 2005-06-22 | 上海交通大学 | 低信噪比红外图像序列中点目标的多帧模糊检测方法 |
CN100385251C (zh) * | 2005-08-18 | 2008-04-30 | 上海交通大学 | 红外大视场环境下的多小目标实时检测方法 |
-
2012
- 2012-07-18 CN CN201210249215.8A patent/CN102819740B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102819740A (zh) | 2012-12-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102819740B (zh) | 一种单帧红外图像弱小目标检测和定位方法 | |
Huang et al. | A new building extraction postprocessing framework for high-spatial-resolution remote-sensing imagery | |
CN104376330B (zh) | 基于超像素散射机制的极化sar图像舰船目标检测方法 | |
CN101727662B (zh) | Sar图像非局部均值去斑方法 | |
CN101980293B (zh) | 一种基于刃边图像的高光谱遥感***mtf检测方法 | |
CN105957054B (zh) | 一种图像变化检测方法 | |
CN102096824B (zh) | 基于选择性视觉注意机制的多光谱图像舰船检测方法 | |
CN103048329B (zh) | 一种基于主动轮廓模型的路面裂缝检测方法 | |
CN103177458B (zh) | 一种基于频域分析的可见光遥感图像感兴趣区域检测方法 | |
CN106886760B (zh) | 一种基于空谱信息结合的高光谱舰船检测方法 | |
CN107392885A (zh) | 一种基于视觉对比机制的红外弱小目标检测方法 | |
CN106530271B (zh) | 一种红外图像显著性检测方法 | |
CN102867195B (zh) | 一种遥感图像多类目标检测和识别方法 | |
CN103824302B (zh) | 基于方向波域图像融合的sar图像变化检测方法 | |
CN103971364A (zh) | 基于加权Gabor小波特征和两级聚类的遥感图像变化检测方法 | |
CN103578110A (zh) | 基于灰度共生矩阵的多波段高分辨率遥感影像分割方法 | |
CN109726649B (zh) | 遥感影像云检测方法、***及电子设备 | |
CN106156758B (zh) | 一种sar海岸图像中海岸线提取方法 | |
Shaoqing et al. | The comparative study of three methods of remote sensing image change detection | |
CN106682678A (zh) | 一种基于支撑域的图像角点检测和分类方法 | |
CN105160648A (zh) | 基于小波和恒虚警率的雷达目标及阴影分割方法 | |
CN102750705A (zh) | 基于图像融合的光学遥感图像变化检测 | |
Liang et al. | Maximum likelihood classification of soil remote sensing image based on deep learning | |
CN102565062B (zh) | 基于图像灰度检测的液体浑浊度测试方法 | |
CN105512622A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |