CN112508923A - 弱小目标检测方法 - Google Patents
弱小目标检测方法 Download PDFInfo
- Publication number
- CN112508923A CN112508923A CN202011478837.9A CN202011478837A CN112508923A CN 112508923 A CN112508923 A CN 112508923A CN 202011478837 A CN202011478837 A CN 202011478837A CN 112508923 A CN112508923 A CN 112508923A
- Authority
- CN
- China
- Prior art keywords
- image
- background
- gradient
- frame
- images
- 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.)
- Granted
Links
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/207—Analysis of motion for motion estimation over a hierarchy of resolutions
-
- 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/10016—Video; Image sequence
-
- 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/20024—Filtering details
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Quality & Reliability (AREA)
- Image Analysis (AREA)
Abstract
本发明旨在提供一种弱小目标检测方法,包括以下步骤:A、对多帧数的待检测图像进行读取;B、构建多个方向的梯度扩散函数,根据梯度扩散函数构建各向梯度差分函数;C、构建背景预测程序,预设其中的停止参数、步长参数及各向梯度差分函数的收敛函参数;D、运用背景预测程序结合各项异性对各帧的待检测图像进行背景预测得到背景预测结果图;E、对各帧的背景预测结果图进行背景估计建模,并将各帧的背景预测结果图与背景估计建模图分别作差得到各帧的图像差分图;F、对多帧数的图像差分图进行拉伸,得到对应的多帧数拉伸图,并进行叠加,即得到小目标的运动轨迹图。该方法克服现有技术缺陷,具有小目标能量增强明显、检测准确度高的特点。
Description
技术领域
本发明涉及图像处理领域,具体涉及一种弱小目标检测方法。
背景技术
弱小目标的检测一直是图像处理的难题,在各种复杂的条件下提取肉眼不可见的目标,对相关的算法的要求是非常高的,有些目标埋没在各种噪声当中,如果没有精良的算法是难以将其检测出来的。处理小目标不仅需要过滤各种复杂背景,还需要考虑到小目标的信号能量,这就需要设计各种算法时考虑最大程度的削弱背景噪声,且最小的影响到目标能量,甚至是在预处理背景的同时对小目标点能量进行加强。
发明内容
本发明旨在提供一种弱小目标检测方法,该方法克服现有技术缺陷,具有小目标能量增强明显、检测准确度高的特点。
本发明的技术方案如下:
一种弱小目标检测算法,包括以下步骤:
A、对多帧数的待检测图像进行读取;
B、构建多个方向的梯度扩散函数,根据梯度扩散函数构建各向梯度差分函数;
C、构建背景预测程序,预设其中的停止参数、步长参数及各向梯度差分函数的收敛函参数;
D、运用背景预测程序结合各项异性对各帧的待检测图像进行背景预测,得到各帧的背景预测结果图;
E、对各帧的背景预测结果图进行背景估计建模,得到各帧的背景估计建模图,将各帧的背景预测结果图及其对应的背景估计建模图分别作差得到各帧的图像差分图;
F、对多帧数的图像差分图进行拉伸,得到对应的多帧数拉伸图,将多帧数拉伸图进行叠加,即得到小目标的运动轨迹图。
优选地,步骤B中所述的多个方向的梯度差分函数包括四个方向的各向梯度扩散函数,分别如下:
CL=1/(1+1/(M×|(f(i,j)-f(i,j-t))|/K));
CR=1/(1+1/(M×|(f(i,j)-f(i,j+t))|/K));
CP=1/(1+1/(M×|(f(i,j)-f(i-t,j))|/K));
CD=1/(1+1/(M×|(f(i,j)-f(i+t,j))|/K));
其中CL、CR、CP、CD分别为四个方向的各向梯度扩散函数,M为限制常数,K为停止参数,f(i,j-t))、f(i,j+t)、f(i-t,j)、f(i+t,j)分别表示四个子区域距离中心像元f(i,j))的步长距离,t表示中心像元到四个方向的步长;
所述的各向梯度差分函数的表达式如下:
优选地,所述的步骤C中的背景预测程序包括以下步骤:
c1、对读取的图像进行边缘不处理;
c2、求取各个像素点四个方向上的梯度值的绝对值;
c3、把各元素四个方向的梯度绝对值带入相对应方向的各向梯度扩散函数中求得各个像素点四个方向的方向系数;
c4、对步骤c3中各个像素点四个方向的方向系数进行求和后作均值。
将各元素各方向的梯度扩散系数CL、CR、CP、CD求平均得到的平均值有利于背景均匀性的预测,保持预测背景的平稳性,进而减少各个像素间预测背景的差异,得到平滑预测背景。
优选地,所述的步骤D中的背景预测包括以下步骤:
d1、读取各元素对应各方向的均值cusm;
d2、把原图和读取到的各元素的四个方向均值做差,达到把各个噪点,目标点以及非平
稳背景抑制的效果;
d3、将抑制后的图像进行拉升输出,即得预测的背景图像。
优选地,所述的步骤E中的背景估计建模包括以下步骤:
e1、读取步骤C中完成背景预测的图像;
e2、读取多帧数的待检测图像,并把待检测图像中各个像素点的像素值与步骤C中完成背景预测的图像的对应像素点的像素值进行做差,并设定阈值,将低于阈值的差值则过滤,高于阈值得差值则保留,从而把图像中已被抑制的目标点、噪声及各类非平稳背景有效过滤,得到平滑的背景估计图,完成背景估计建模。
本发明方案采用特有的梯度扩散函数构建各向梯度差分函数,能够有效减少梯度比较低的小目标的漏检的可能性;本发明方案的各向梯度扩散函数为单调递增函数,能够在背景建模后明显增强目标点的信号,更好的保留目标点,从而使得即使在梯度很低的情况下也能够很好地对图像进行检测,同时在背景建模的时候更好的保留小目标的信号能量,因而对小目标检测有较好的效果;并且,在对背景处理时能够有效去除大部分的非平稳背景,使得目标点以较大的能量值保留下来,便于后续处理,从而提高检测成功率。
附图说明
图1为本发明实施例1提供的弱小目标检测方法的背景预测图;
图2为本发明实施例1提供的弱小目标检测方法的背景估计建模图;
图3为本发明实施例1提供的弱小目标检测方法的图像差分图(第一帧);
图4为本发明实施例1提供的弱小目标检测方法的拉伸图;
图5为本发明实施例1提供的弱小目标检测方法的目标运动轨迹图;
具体实施方式
如图1-5所示,本实施例提供的弱小目标检测算法,包括以下步骤:
A、对多帧数的待检测图像进行读取;
B、构建多个方向的梯度扩散函数,根据梯度扩散函数构建各向梯度差分函数;
C、构建背景预测程序,预设其中的停止参数、步长参数及各向梯度差分函数的收敛函参数;
步骤B中所述的多个方向的梯度差分函数包括四个方向的各向梯度扩散函数,分别如下:
CL=1/(1+1/(M×|(f(i,j)-f(i,j-t))|/K));
CR=1/(1+1/(M×|(f(i,j)-f(i,j+t))|/K));
CP=1/(1+1/(M×|(f(i,j)-f(i-t,j))|/K));
CD=1/(1+1/(M×|(f(i,j)-f(i+t,j))|/K));
其中CL、CR、CP、CD分别为四个方向的各向梯度扩散函数,M为限制常数,K为停止参数,f(i,j-t))、f(i,j+t)、f(i-t,j)、f(i+t,j)分别表示四个子区域距离中心像元f(i,j))的步长距离,t表示中心像元到四个方向的步长;
所述的各向梯度差分函数的表达式如下:
所述的步骤C中的背景预测程序包括以下步骤:
c1、对读取的图像进行边缘不处理;
c2、求取各个像素点四个方向上的梯度值的绝对值;
c3、把各元素四个方向的梯度绝对值带入相对应方向的各向梯度扩散函数中求得各个像素点四个方向的方向系数;
c4、对步骤c3中各个像素点四个方向的方向系数进行求和后作均值。
D、运用背景预测程序结合各项异性对各帧的待检测图像进行背景预测,得到各帧的背景预测结果图;
所述的步骤D中的背景预测包括以下步骤:
d1、读取各元素对应各方向的均值cusm;
d2、把原图和读取到的各元素的四个方向均值做差,达到把各个噪点,目标点以及非平
稳背景抑制的效果;
d3、将抑制后的图像进行拉升输出,即得预测的背景图像。在本实施例中得到如图2所示的平滑背景图预测图。
E、对各帧的背景预测结果图进行背景估计建模,得到各帧的背景估计建模图,将各帧的背景预测结果图及其对应的背景估计建模图分别作差得到各帧的图像差分图;所述的步骤E中的背景估计建模包括以下步骤:
所述的步骤E中的背景估计建模包括以下步骤:
e1、读取步骤C中完成背景预测的图像;
e2、读取多帧数的待检测图像,并把待检测图像中各个像素点的像素值与步骤C中完成背景预测的图像的对应像素点的像素值进行做差,并设定阈值,将低于阈值的差值则过滤,高于阈值得差值则保留,从而把图像中已被抑制的目标点、噪声及各类非平稳背景有效过滤,使得图像中只含有目标点和少许噪点,得到平滑的背景估计图,完成背景估计建模。
F、对多帧数的图像差分图进行拉伸,得到对应的多帧数拉伸图,将多帧数拉伸图进行叠加,即得到小目标的运动轨迹图;
本实施例预设停止参数为100、步长参数为4、各向梯度差分函数的收敛函参数设置为20;如图2所示,在进行背景建模后,背景建模图对目标的抑制效果很好,方便后续的处理;如图4所示,图中可清晰的观察到目标与其他少量的噪声点,如图5所示,图中清晰显示出各个目标点叠加后的运动轨迹。
Claims (5)
1.一种弱小目标检测算法,其特征在于包括以下步骤:
A、对多帧数的待检测图像进行读取;
B、构建多个方向的梯度扩散函数,根据梯度扩散函数构建各向梯度差分函数;
C、构建背景预测程序,预设其中的停止参数、步长参数及各向梯度差分函数的收敛函参数;
D、运用背景预测程序结合各项异性对各帧的待检测图像进行背景预测,得到各帧的背景预测结果图;
E、对各帧的背景预测结果图进行背景估计建模,得到各帧的背景估计建模图,将各帧的背景预测结果图及其对应的背景估计建模图分别作差得到各帧的图像差分图;
F、对多帧数的图像差分图进行拉伸,得到对应的多帧数拉伸图,将多帧数拉伸图进行叠加,即得到小目标的运动轨迹图。
2.如权利要求1所述的弱小目标检测算法,其特征在于:
步骤B中所述的多个方向的梯度差分函数包括四个方向的各向梯度扩散函数,分别如下:
CL=1/(1+1/(M×|(f(i,j)-f(i,j-t))|/K));
CR=1/(1+1/(M×|(f(i,j)-f(i,j+t))|/K));
CP=1/(1+1/(M×|(f(i,j)-f(i-t,j))|/K));
CD=1/(1+1/(M×|(f(i,j)-f(i+t,j))|/K));
其中CL、CR、CP、CD分别为四个方向的各向梯度扩散函数,M为限制常数,K为停止参数,f(i,j-t))、f(i,j+t)、f(i-t,j)、f(i+t,j)分别表示四个子区域距离中心像元f(i,j))的步长距离,t表示中心像元到四个方向的步长;
所述的各向梯度差分函数的表达式如下:
3.如权利要求2所述的弱小目标检测算法,其特征在于:
所述的步骤C中的背景预测程序包括以下步骤:
c1、对读取的图像进行边缘不处理;
c2、求取各个像素点四个方向上的梯度值的绝对值;
c3、把各元素四个方向的梯度绝对值带入相对应方向的各向梯度扩散函数中求得各个像素点四个方向的方向系数;
c4、对步骤c3中各个像素点四个方向的方向系数进行求和后作均值。
4.如权利要求3所述的弱小目标检测算法,其特征在于:
所述的步骤D中的背景预测包括以下步骤:
d1、读取各元素对应各方向的均值cusm;
d2、把原图和读取到的各元素的四个方向均值做差,达到把各个噪点,目标点以及非平稳背景抑制的效果;
d3、将抑制后的图像进行拉升输出,即得预测的背景图像。
5.如权利要求4所述的弱小目标检测算法,其特征在于:
所述的步骤E中的背景估计建模包括以下步骤:
e1、读取步骤C中完成背景预测的图像;
e2、读取多帧数的待检测图像,并把待检测图像中各个像素点的像素值与步骤C中完成背景预测的图像的对应像素点的像素值进行做差,并设定阈值,将低于阈值的差值则过滤,高于阈值得差值则保留,从而把图像中已被抑制的目标点、噪声及各类非平稳背景有效过滤,得到平滑的背景估计图,完成背景估计建模。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011478837.9A CN112508923B (zh) | 2020-12-15 | 2020-12-15 | 弱小目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011478837.9A CN112508923B (zh) | 2020-12-15 | 2020-12-15 | 弱小目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112508923A true CN112508923A (zh) | 2021-03-16 |
CN112508923B CN112508923B (zh) | 2023-06-23 |
Family
ID=74973744
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011478837.9A Active CN112508923B (zh) | 2020-12-15 | 2020-12-15 | 弱小目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112508923B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113870308A (zh) * | 2021-09-07 | 2021-12-31 | 广西科技大学 | 一种螺旋梯度优化估计的弱小目标检测方法 |
CN114155426A (zh) * | 2021-12-13 | 2022-03-08 | 中国科学院光电技术研究所 | 一种基于局部多向梯度信息融合的弱小目标检测方法 |
CN115144830A (zh) * | 2022-04-08 | 2022-10-04 | 中国人民解放军63983部队 | 一种强地物杂波干扰下的小目标雷达检测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000061261A2 (en) * | 1999-04-06 | 2000-10-19 | Massachusetts Institute Of Technology | Thermophoretic pump and concentrator |
US6850646B1 (en) * | 1997-12-31 | 2005-02-01 | Cognex Corporation | Fast high-accuracy multi-dimensional pattern inspection |
CN103279930A (zh) * | 2013-05-27 | 2013-09-04 | 辽宁工程技术大学 | 一种图像的同步去噪增强方法 |
CN105005975A (zh) * | 2015-07-08 | 2015-10-28 | 南京信息工程大学 | 基于pcnn和图像熵的各向异性扩散的图像去噪方法 |
CN106251344A (zh) * | 2016-07-26 | 2016-12-21 | 北京理工大学 | 一种基于视觉感受野的多尺度红外目标自适应检测方法 |
CN106469313A (zh) * | 2016-09-30 | 2017-03-01 | 中国科学院光电技术研究所 | 一种管径自适应时空域滤波的弱小目标检测方法 |
CN106780545A (zh) * | 2016-11-25 | 2017-05-31 | 中国科学院云南天文台 | 一种结合时空域的弱小目标能量累积增强方法 |
CN107392885A (zh) * | 2017-06-08 | 2017-11-24 | 江苏科技大学 | 一种基于视觉对比机制的红外弱小目标检测方法 |
CN110490904A (zh) * | 2019-08-12 | 2019-11-22 | 中国科学院光电技术研究所 | 一种弱小目标检测与跟踪方法 |
CN210725109U (zh) * | 2019-12-06 | 2020-06-09 | 广西科技大学 | 智能教室管理系统 |
-
2020
- 2020-12-15 CN CN202011478837.9A patent/CN112508923B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6850646B1 (en) * | 1997-12-31 | 2005-02-01 | Cognex Corporation | Fast high-accuracy multi-dimensional pattern inspection |
WO2000061261A2 (en) * | 1999-04-06 | 2000-10-19 | Massachusetts Institute Of Technology | Thermophoretic pump and concentrator |
CN103279930A (zh) * | 2013-05-27 | 2013-09-04 | 辽宁工程技术大学 | 一种图像的同步去噪增强方法 |
CN105005975A (zh) * | 2015-07-08 | 2015-10-28 | 南京信息工程大学 | 基于pcnn和图像熵的各向异性扩散的图像去噪方法 |
CN106251344A (zh) * | 2016-07-26 | 2016-12-21 | 北京理工大学 | 一种基于视觉感受野的多尺度红外目标自适应检测方法 |
CN106469313A (zh) * | 2016-09-30 | 2017-03-01 | 中国科学院光电技术研究所 | 一种管径自适应时空域滤波的弱小目标检测方法 |
CN106780545A (zh) * | 2016-11-25 | 2017-05-31 | 中国科学院云南天文台 | 一种结合时空域的弱小目标能量累积增强方法 |
CN107392885A (zh) * | 2017-06-08 | 2017-11-24 | 江苏科技大学 | 一种基于视觉对比机制的红外弱小目标检测方法 |
CN110490904A (zh) * | 2019-08-12 | 2019-11-22 | 中国科学院光电技术研究所 | 一种弱小目标检测与跟踪方法 |
CN210725109U (zh) * | 2019-12-06 | 2020-06-09 | 广西科技大学 | 智能教室管理系统 |
Non-Patent Citations (1)
Title |
---|
李甜田;张琪;张勇;张熙;: "基于动静结合的单一背景目标检测方法", 河南科技, no. 13, pages 34 - 38 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113870308A (zh) * | 2021-09-07 | 2021-12-31 | 广西科技大学 | 一种螺旋梯度优化估计的弱小目标检测方法 |
CN113870308B (zh) * | 2021-09-07 | 2024-03-22 | 广西科技大学 | 一种螺旋梯度优化估计的弱小目标检测方法 |
CN114155426A (zh) * | 2021-12-13 | 2022-03-08 | 中国科学院光电技术研究所 | 一种基于局部多向梯度信息融合的弱小目标检测方法 |
CN114155426B (zh) * | 2021-12-13 | 2023-08-15 | 中国科学院光电技术研究所 | 一种基于局部多向梯度信息融合的弱小目标检测方法 |
CN115144830A (zh) * | 2022-04-08 | 2022-10-04 | 中国人民解放军63983部队 | 一种强地物杂波干扰下的小目标雷达检测方法 |
CN115144830B (zh) * | 2022-04-08 | 2023-10-31 | 中国人民解放军63983部队 | 一种强地物杂波干扰下的小目标雷达检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112508923B (zh) | 2023-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xuan et al. | An improved canny edge detection algorithm | |
CN112508923A (zh) | 弱小目标检测方法 | |
WO2012074361A1 (en) | Method of image segmentation using intensity and depth information | |
CN110580709A (zh) | 一种基于ViBe和三帧差分融合的目标检测方法 | |
Niu et al. | Research and analysis of threshold segmentation algorithms in image processing | |
CN107085836B (zh) | 一种在运动物体检测中通用的鬼影消除方法 | |
JP7449507B2 (ja) | カメラストリームのためのマスクを生成する方法、コンピュータープログラム製品およびコンピューター可読媒体 | |
Lian et al. | A novel method on moving-objects detection based on background subtraction and three frames differencing | |
CN109255799B (zh) | 一种基于空间自适应相关滤波器的目标跟踪方法及系统 | |
CN110660048B (zh) | 一种基于形状特征的皮革表面缺陷检测方法 | |
CN116563262A (zh) | 基于多模态的建筑裂缝检测算法 | |
Liu et al. | Moving target detection based on improved ghost suppression and adaptive visual background extraction | |
CN111161308A (zh) | 一种基于关键点匹配的双波段融合目标提取方法 | |
CN111259866B (zh) | 一种基于改进背景差分法的海上舰船目标检测方法 | |
Zhang et al. | Dehazing with improved heterogeneous atmosphere light estimation and a nonlinear color attenuation prior model | |
CN114067122B (zh) | 一种两级式二值化图像处理方法 | |
CN116363064A (zh) | 融合目标检测模型和图像分割模型的缺陷识别方法及装置 | |
CN114998186A (zh) | 基于图像处理的铜始极片表面结疤缺陷检测方法及系统 | |
CN114554188A (zh) | 基于图像传感器像素阵列的手机摄像头检测方法及装置 | |
Zhao et al. | A novel method for moving object detection in intelligent video surveillance systems | |
CN111798506A (zh) | 一种图像处理方法、控制方法、终端及计算机可读存储介质 | |
CN117876379B (zh) | 基于图像特征的锚杆缺陷智能化检测方法 | |
Srikrishna et al. | Realization of Human Eye Pupil Detection System using Canny Edge Detector and Circular Hough Transform Technique | |
CN113643315B (zh) | 一种基于自适应峰值梯度下降滤波器的红外小目标检测方法 | |
Wu et al. | Moving-object Detection Based on Shadow Removal and Prospect Reconstruction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |