CN104268877A - 一种红外图像海天线自适应检测方法 - Google Patents
一种红外图像海天线自适应检测方法 Download PDFInfo
- Publication number
- CN104268877A CN104268877A CN201410505799.XA CN201410505799A CN104268877A CN 104268877 A CN104268877 A CN 104268877A CN 201410505799 A CN201410505799 A CN 201410505799A CN 104268877 A CN104268877 A CN 104268877A
- Authority
- CN
- China
- Prior art keywords
- infrared image
- sea horizon
- sea
- line
- len
- 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
Classifications
-
- 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/13—Edge detection
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种红外图像的海天线自适应检测方法,包括以下步骤:(1)获取原始红外图像Forg;(2)利用Canny算子对原始红外图像进行边缘检测,得到边缘图像Fedge;(3)建立划窗模型MWin和海天线长度和模型LENsum(j);(4)使用划窗模型MWin对边缘图像Fedge进行全图遍历,求取不同行的海天线长度和;(5)将海天线长度和LENsum最大的行定义为当前图像海天线位置;(6)对当前图像海天线进行滤波,输出第p+1帧海天线滤波值Lineseasky(p+1)。
Description
技术领域
本发明涉及红外图像处理领域,特别涉及一种适合红外硬件实时实现的红外图像海天线自适应检测方法。
背景技术
海天线检测是海面舰船目标检测的关键环节,特别是红外探测器在成像距离较远时,受海面潮湿大气影响,舰船目标红外辐射衰减剧烈,在图像中呈现为弱小目标,大了舰船检测和定位的难度。而反舰导弹发射后飞行高度比较低,加上地球曲率的影响,如果海面舰船目标出现在探测器视场内,就一定位于海天线的附近区域。基于这一先验条件,如果能准确定位海天线就可以有效排除海面杂波、人工干扰的影响,缩小目标搜索范围,提高目标检测效率。
目前,海天线检测技术受到了国内外学者的广泛关注,主要包括直线拟合法、Hough变换方法、OTSU方法、基于多尺度小波变换的方法等。现有的海天线检测技术存在以下缺点:(1)直线拟合法在某些复杂背景条件下,图像受到海杂波、太阳光、云层等自然背景的干扰,难以得到满意的提取效果;(2)Hough变换检测直线,容易受天空中云团干扰,形成伪轮廓;(3)Otsu阈值分割方法抗噪声能力差,且不易确定最佳阈值;(4)基于小波变换的方法,即利用不同尺度下小波系数抑制背景和提取边缘,但检测效率不高。
发明内容
发明目的:本发明所要解决的技术问题是针对现有技术的不足,提供一种红外图像的海天线自适应检测方法。
为了解决上述技术问题,本发明公开了一种红外图像的海天线自适应检测方法,包括以下步骤:
(1)获取原始红外图像Forg;
(2)利用Canny算子对原始红外图像进行边缘检测,得到边缘图像Fedge;
(3)建立划窗模型MWin和海天线长度和模型LENsum(j);
(4)使用划窗模型MWin对边缘图像Fedge进行全图遍历,求取不同行的海天线长度和;
(5)将海天线长度和LENsum最大的行定义为当前图像海天线位置;
(6)对当前图像海天线进行滤波,输出第p+1帧海天线滤波值Lineseasky(p+1)。
本发明一种红外图像海天线自适应检测方法中,使用Canny算子计算原图图像的边缘图像Fedge。
本发明一种红外图像海天线自适应检测方法中,使用5行划窗模型对图像进行遍历。
划窗模型以当前行为中心,向上下行分别扩展,扩展行随着与中心距离的增加而减少划窗模型中的像素数量,则5行划窗模型定义如下:
MWin={{0,2},{-1,1},{0,1},{1,1},{-2,0},{-1,0},{0,0},{1,0},{2,0},{-1,-1},{0,-1},{1,-1},{0,-2}}
海天线长度和模型根据连通线的长度确定加权系数,连通线越长反映是海天线的可能性越大,则加权系数也越大。将最大加权系数定义为COEmax,图像行数定义为M,列数定义为N,连通线的长度分别为{Len(k)|k=1,2,...,A},其中,k为自变量,A为连通线数量,则第j行的海天线长度和模型LENsum(j)定义如下:
本发明一种红外图像海天线自适应检测方法中,连通线长度{Len(k)|k=1,2,...,A}定义如下:
Len(k)=Pixelright(j)-Pixelleft(j)
其中,Pixelleft(j)为第j行划窗MWin从左向右遍历时连通区域最左边的像素点位置,Pixelright(j)为第j行划窗MWin从左向右遍历时连通区域最右边的像素点位置,k表示自变量,A表示当前行的连通线数量。
本发明一种红外图像海天线自适应检测方法中,使用5行划窗模型对图像进行遍历时,判断连通线是结束的准则是:当划窗模型中5行的最右边像素同时为0时,认为当前连通线结束。
本发明一种红外图像海天线自适应检测方法中,第p+1帧海天线滤波值Lineseasky(p+1)的计算公式如下:
其中,Num表示滤波系数,Lineseasky(p)表示当前红外图像第p帧的海天线滤波值,j为步骤(5)中计算得到的当前帧红外图像海天线位置,Num=100。
本发明与现有技术相比,具有以下显著优点:(1)使用划窗模型和海天线长度和模型计算海天线,海天线检测精度高;(2)利用统计原理对海天线进行滤波处理,避免噪声干扰,能有效提高海天线检测精确;(3)抗海杂波、天空云团干扰能力强;(4)不存在高阶运算和复杂结构,算法运算量小,易于硬件实时实现。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/或其他方面的优点将会变得更加清楚。
图1为5行划窗模型图。
图2为本发明流程图。
图3a为实施例在简单背景下海天线检测结果。
图3b为实施例在复杂背景下海天线检测结果。
具体实施方式
本发明公开了一种红外图像海天线自适应检测方法,包括以下步骤:
(1)获取原始红外图像Forg;
(2)利用Canny算子对原始红外图像进行边缘检测,得到边缘图像Fedge;
(3)建立划窗模型MWin和海天线长度和模型LENsum(j);
(4)使用划窗模型MWin对边缘图像Fedge进行全图遍历,求取不同行的海天线长度和;
(5)将海天线长度和LENsum(j)最大的行定义为当前图像海天线位置;
(6)对当前图像海天线进行滤波,输出第p+1帧海天线滤波值Lineseasky(p+1)。
步骤(2)中,使用Canny算子对原始红外图像进行边缘检测,得到边缘图像。Canny算子具体包括:用高斯滤波器来对图像滤波;用高斯算子一阶微分对图像进行滤波,得到每个像素梯度的大小和方向;对梯度进行“非极大抑制”;对梯度取高低两个阈值,构造目标边缘。
步骤(3)中,选用5行划窗模型MWin,即
MWin={{0,2},{-1,1},{0,1},{1,1},{-2,0},{-1,0},{0,0},{1,0},{2,0},{-1,-1},{0,-1},{1,-1},{0,-2}}
连通线长度{Len(k)|k=1,2,...,A}定义如下:
Len(k)=Pixelright(j)-Pixelleft(j)
其中,Pixelleft(j)为第j行划窗MWin从左向右遍历时连通区域最左边的像素点位置,Pixelright(j)为第j行划窗MWin从左向右遍历时连通区域最右边的像素点位置,k表示自变量,A表示当前行的连通线数量。
第j行海天线长度和模型LENsum(j)定义如下:
其中,最大加权系数定义为COEmax,图像行数定义为M,列数定义为N,连通
线的长度分别为{Len(k)|k=1,2,...,A},其中,k为自变量,A为连通线数量。
步骤(6)中,第p+1帧海天线滤波值Lineseasky(p+1)的计算公式如下:
其中,Num表示滤波系数,Lineseasky(p)表示第p帧的海天线滤波值,j表示第p+1帧图像海天线位置。
实施例1
结合图1,下面以实例来说明本发明一种红外图像的海天线自适应检测方法。红外图像的像素个数320×256,帧频50HZ。红外热像仪的数字信号通过光纤传给DSP+FPGA架构的专用图像处理板,红外图像的海天线自适应检测方法在DSP处理器中实现,处理时间小于20ms,满足实时处理的需求,具体实施步骤如下:
(1)获取原始红外图像Forg;
原始红外图像Forg是14-bit数字图像。红外热像仪的数字信号通过光纤打包方式传给专用图像处理板的FPGA,FPGA完成光纤解包将数字信号经DSP的EMIF接口传给DSP,得到14-bit原始红外图像Forg。
(2)利用Canny算子对图像Forg进行边缘检测,得到边缘图像Fedge;
(3)建立5行划窗模型MWin,即
MWin={{0,2},{-1,1},{0,1},{1,1},{-2,0},{-1,0},{0,0},{1,0},{2,0},{-1,-1},{0,-1},{1,-1},{0,-2}}
5行划窗模型MWin如图1所示:
计算第j行连通线长度{Len(k)|k=1,2,...,A}长度:
Len(k)=Pixelright(j)-Pixelleft(j)
其中,Pixelleft(j)为第j行划窗MWin从左向右遍历时连通区域最左边的像素点位置,Pixelright(j)为第j行划窗MWin从左向右遍历时连通区域最右边的像素点位置,k表示自变量,A表示当前行的连通线数量。
计算第j行海天线长度和LENsum(j):
其中,最大加权系数定义为COEmax,0≤COEmax≤1,图像行数定义为256,列数定义为320,连通线的长度分别为{Len(k)|k=1,2,...,A},其中,k为自变量,A≤320为连通线数量,1≤j≤256。
(4)使用划窗模型MWin对边缘图像Fedge进行全图遍历,求取第j行海天线长度和LENsum(j),1≤j≤256;
(5)将海天线长度和LENsum(j)最大的行定义为当前图像海天线位置j,1≤j≤256;
(6)对当前图像海天线进行滤波,输出第p+1帧海天线滤波值Lineseasky(p+1)。
第p+1帧海天线滤波值Lineseasky(p+1)的计算公式如下:
其中,Num表示滤波系数,Num=100,Lineseasky(p)表示第p帧的海天线滤波值,j表示第p+1帧图像海天线位置,1≤j≤256。
实施例2
图3中,图3a表示简单背景下海天线检测结果,图3b表示复杂背景下海天线检测结果。图像显示和数据仿真表明,本发明提出的海天线自适应检测方法的检测精度是正负1个像素,海天线检测精确高。
本发明提供了一种红外图像的海天线自适应检测方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (6)
1.一种红外图像海天线自适应检测方法,其特征在于,包括以下步骤:
(1)获取原始红外图像Forg;
(2)对原始红外图像进行边缘检测,得到边缘图像Fedge;
(3)建立划窗模型MWin和海天线长度和模型LENsum(j);
(4)使用划窗模型MWin对边缘图像Fedge进行全图遍历,求取不同行的海天线长度和LENsum;
(5)将海天线长度和LENsum最大的行数定义为当前帧红外图像的海天线位置j;
(6)对当前帧第p帧红外图像海天线进行滤波,输出红外图像第p+1帧海天线滤波值Lineseasky(p+1)。
2.根据权利要求1所述的一种红外图像的海天线自适应检测方法,其特征在于,步骤(2)中,使用Canny算子计算原图图像的边缘图像Fedge。
3.根据权利要求1所述的一种红外图像的海天线自适应检测方法,其特征在于,步骤(3)中,使用5行划窗模型对红外图像进行遍历;5行划窗模型MWin定义如下:
MWin={{0,2},{-1,1},{0,1},{1,1},{-2,0},{-1,0},{0,0},{1,0},{2,0},{-1,-1},{0,-1},{1,-1},{0,-2}},
则第j行的海天线长度和模型LENsum(j)计算公式如下:
COEmax为最大加权系数,0≤COEmax≤1,N为红外图像列数,{Len(k)|k=1,2,...,A}为连通线的长度,k为第k个连通线,A为连通线数量。
4.根据权利要求3所述的一种红外图像的海天线自适应检测方法,其特征在于,步骤(3)中,连通线长度{Len(k)|k=1,2,...,A}计算公式如下:
Len(k)=Pixelright(j)-Pixelleft(j),
其中,Pixelleft(j)为第j行海天线的划窗MWin从左向右遍历时连通区域最左边的像素点位置,Pixelright(j)为第j行海天线划窗MWin从左向右遍历时连通区域最右边的像素点位置,1≤j≤M,M表示海天线的总行数。
5.根据权利要求4所述的一种红外图像的海天线自适应检测方法,其特征在于,步骤(3)中,使用5行划窗模型对红外图像进行遍历时,判断连通线是结束的准则是:当划窗模型中5行的最右边像素同时为0时,认为当前连通线结束。
6.根据权利要求5所述的一种红外图像的海天线自适应检测方法,其特征在于,步骤(6)中,第p+1帧海天线滤波值Lineseasky(p+1)的计算公式如下:
其中,Num表示滤波系数,Lineseasky(p)表示当前红外图像第p帧的海天线滤波值,j为步骤(5)中计算得到的当前帧红外图像海天线位置,Num=100。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410505799.XA CN104268877B (zh) | 2014-09-26 | 2014-09-26 | 一种红外图像海天线自适应检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410505799.XA CN104268877B (zh) | 2014-09-26 | 2014-09-26 | 一种红外图像海天线自适应检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104268877A true CN104268877A (zh) | 2015-01-07 |
CN104268877B CN104268877B (zh) | 2017-04-05 |
Family
ID=52160396
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410505799.XA Active CN104268877B (zh) | 2014-09-26 | 2014-09-26 | 一种红外图像海天线自适应检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104268877B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427286A (zh) * | 2015-11-06 | 2016-03-23 | 中国电子科技集团公司第二十八研究所 | 一种基于灰度和梯度分割的红外目标检测方法 |
CN105469390A (zh) * | 2015-11-17 | 2016-04-06 | 哈尔滨工程大学 | 一种基于改进Seam Carving的全景海天线提取方法 |
CN105894516A (zh) * | 2016-04-18 | 2016-08-24 | 江苏科技大学 | 基于海天线旋转后准确确定目标区域的方法 |
CN105931228A (zh) * | 2016-04-15 | 2016-09-07 | 中国科学院上海技术物理研究所 | 一种基于直线查找矩阵的红外图像海天线检测方法 |
CN106096617A (zh) * | 2016-06-08 | 2016-11-09 | 哈尔滨工程大学 | 一种基于改进主动轮廓模型的全景海天线提取方法 |
CN106127728A (zh) * | 2016-06-07 | 2016-11-16 | 电子科技大学 | 一种海空背景下的红外图像海天线连通域检测方法 |
CN106780366A (zh) * | 2016-11-25 | 2017-05-31 | 中国科学院上海技术物理研究所 | 一种红外图像海天线信噪比增强方法 |
CN106952253A (zh) * | 2017-03-14 | 2017-07-14 | 哈尔滨工程大学 | 一种基于行‑信息熵的海天线提取方法 |
CN109978869A (zh) * | 2019-03-29 | 2019-07-05 | 清华大学 | 一种基于灰度共生矩阵和Hough变换的海天线检测方法与系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102054275A (zh) * | 2010-12-13 | 2011-05-11 | 中国科学院长春光学精密机械与物理研究所 | 灰度图像海天线实时检测方法 |
CN102279973A (zh) * | 2010-06-11 | 2011-12-14 | 中国兵器工业第二○五研究所 | 基于高梯度关键点的海天线检测方法 |
-
2014
- 2014-09-26 CN CN201410505799.XA patent/CN104268877B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102279973A (zh) * | 2010-06-11 | 2011-12-14 | 中国兵器工业第二○五研究所 | 基于高梯度关键点的海天线检测方法 |
CN102054275A (zh) * | 2010-12-13 | 2011-05-11 | 中国科学院长春光学精密机械与物理研究所 | 灰度图像海天线实时检测方法 |
Non-Patent Citations (3)
Title |
---|
冯德瀛 等: "红外图像中海天线的检测", 《红外》 * |
夏平 等: "基于形态学处理和最小二乘法的海天线检测算法", 《光学与光电技术》 * |
赵峰民 等: "基于最长曲线法的海天线检测算法", 《兵工自动化》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427286B (zh) * | 2015-11-06 | 2017-12-05 | 中国电子科技集团公司第二十八研究所 | 一种基于灰度和梯度分割的红外目标检测方法 |
CN105427286A (zh) * | 2015-11-06 | 2016-03-23 | 中国电子科技集团公司第二十八研究所 | 一种基于灰度和梯度分割的红外目标检测方法 |
CN105469390A (zh) * | 2015-11-17 | 2016-04-06 | 哈尔滨工程大学 | 一种基于改进Seam Carving的全景海天线提取方法 |
CN105469390B (zh) * | 2015-11-17 | 2018-04-27 | 哈尔滨工程大学 | 一种基于改进Seam Carving的全景海天线提取方法 |
CN105931228A (zh) * | 2016-04-15 | 2016-09-07 | 中国科学院上海技术物理研究所 | 一种基于直线查找矩阵的红外图像海天线检测方法 |
CN105894516A (zh) * | 2016-04-18 | 2016-08-24 | 江苏科技大学 | 基于海天线旋转后准确确定目标区域的方法 |
CN105894516B (zh) * | 2016-04-18 | 2018-08-31 | 江苏科技大学 | 基于海天线旋转后准确确定目标区域的方法 |
CN106127728A (zh) * | 2016-06-07 | 2016-11-16 | 电子科技大学 | 一种海空背景下的红外图像海天线连通域检测方法 |
CN106127728B (zh) * | 2016-06-07 | 2018-08-17 | 电子科技大学 | 一种海空背景下的红外图像海天线连通域检测方法 |
CN106096617A (zh) * | 2016-06-08 | 2016-11-09 | 哈尔滨工程大学 | 一种基于改进主动轮廓模型的全景海天线提取方法 |
CN106096617B (zh) * | 2016-06-08 | 2019-10-08 | 哈尔滨工程大学 | 一种基于改进主动轮廓模型的全景海天线提取方法 |
CN106780366A (zh) * | 2016-11-25 | 2017-05-31 | 中国科学院上海技术物理研究所 | 一种红外图像海天线信噪比增强方法 |
CN106952253A (zh) * | 2017-03-14 | 2017-07-14 | 哈尔滨工程大学 | 一种基于行‑信息熵的海天线提取方法 |
CN109978869A (zh) * | 2019-03-29 | 2019-07-05 | 清华大学 | 一种基于灰度共生矩阵和Hough变换的海天线检测方法与系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104268877B (zh) | 2017-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104268877A (zh) | 一种红外图像海天线自适应检测方法 | |
CN104299229B (zh) | 一种基于时空域背景抑制的红外弱小目标检测方法 | |
CN104657945B (zh) | 复杂背景下多尺度时空联合滤波的红外小目标检测方法 | |
CN102855622B (zh) | 一种基于显著性分析的红外遥感图像海面船只检测方法 | |
CN105574855A (zh) | 云背景下基于模板滤波和虚警抑制的红外小目标检测方法 | |
CN107403433A (zh) | 一种复杂云背景下红外小目标检测方法 | |
CN107403134B (zh) | 基于局部梯度三边的图域多尺度红外弱小目标检测方法 | |
CN103745216B (zh) | 一种基于空域特性的雷达图像杂波抑制方法 | |
CN105654091A (zh) | 海面目标检测方法及装置 | |
CN103679694B (zh) | 一种基于全景视觉的舰船小目标检测方法 | |
CN107329137B (zh) | 基于box样条小波框架的探地雷达地下目标定位方法及系统 | |
CN109711256B (zh) | 一种低空复杂背景无人机目标检测方法 | |
CN105138983B (zh) | 基于加权部件模型和选择性搜索分割的行人检测方法 | |
CN106886747B (zh) | 一种基于扩展小波变换的复杂背景下舰船检测方法 | |
CN107392095A (zh) | 一种基于掩码图像的红外弱小目标检测算法 | |
Wang et al. | An infrared maritime target detection algorithm applicable to heavy sea fog | |
CN105741276A (zh) | 一种船舶水线提取方法 | |
CN104834915A (zh) | 一种复杂云天背景下小红外目标检测方法 | |
CN103646389A (zh) | 一种基于几何模型的sar斜距图像同名点自动提取方法 | |
CN104537688A (zh) | 一种基于背景差分和hog特征的运动目标检测方法 | |
CN104199009B (zh) | 一种基于时域特性的雷达图像杂波抑制方法 | |
CN103413138B (zh) | 一种红外图像序列中点目标检测方法 | |
Lijun et al. | Geo‐information mapping improves Canny edge detection method | |
CN102750675A (zh) | 一种斑点噪声污染图像的非局部均值滤波方法 | |
CN106780545A (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 |