CN103714542A - 低分辨率高频声纳图像中目标亮点提取方法 - Google Patents

低分辨率高频声纳图像中目标亮点提取方法 Download PDF

Info

Publication number
CN103714542A
CN103714542A CN201310726052.2A CN201310726052A CN103714542A CN 103714542 A CN103714542 A CN 103714542A CN 201310726052 A CN201310726052 A CN 201310726052A CN 103714542 A CN103714542 A CN 103714542A
Authority
CN
China
Prior art keywords
sonar image
pixel
bright spot
target
image
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
Application number
CN201310726052.2A
Other languages
English (en)
Other versions
CN103714542B (zh
Inventor
姜可宇
蔡志明
王超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Naval University of Engineering PLA
Original Assignee
Naval University of Engineering PLA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Naval University of Engineering PLA filed Critical Naval University of Engineering PLA
Priority to CN201310726052.2A priority Critical patent/CN103714542B/zh
Publication of CN103714542A publication Critical patent/CN103714542A/zh
Application granted granted Critical
Publication of CN103714542B publication Critical patent/CN103714542B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明属于声纳图像处理技术领域,具体涉及一种低分辨率高频声纳图像中目标亮点提取方法。它包括以下步骤:确定目标亮点在声纳图像中距离维和角度维上的最大尺寸;计算声纳图像中各像素点的信号强度;计算声纳图像中各像素点的背景强度;逐点计算声纳图像中每个像素点的信混比,并据此确定各像素点的声纳图像背景;将原始声纳图像与声纳图像背景相减后,对图像进行逐点扫描,得到目标亮点图像。该提取方法简单、有效,提取结果准确、适用性高,解决了低分辨率高频声纳图像中小目标检测、辨识难的问题。

Description

低分辨率高频声纳图像中目标亮点提取方法
技术领域
本发明属于声纳图像处理技术领域,具体涉及一种低分辨率高频声纳图像中目标亮点提取方法。
背景技术
高频成像声纳对水下物体进行回波成像,在声纳工作频率不是足够高的情况下,声纳图像的分辨率较低,对于尺度在2m以下的海底物体来说,不能形成明显的声影区,给操作员发现类似小目标造成了很大困难。对于高分辨率的声纳图像,有很多图像分割的方法来提取目标的形状和大小,进而识别水下目标,但对于低分辨率的声纳图像,目前尚没有见到目标亮点提取方法的公开报道。
混响是主动声纳的固有干扰,对于高频成像声纳来说,海底混响构成了声纳图像的混响背景,一般起伏不大,当目标信混比足够强时,目标亮点会从声纳图像背景中显现出来,人眼很容易辨别。但在目标信混比不是足够强的情况,受混响背景的干扰,人眼很难辨别其中的小目标亮点,也无法继续进行目标的自动检测和跟踪等处理。
声纳图像预处理包括各种去噪和图像增强等技术,其根本目的是最大限度地消除各类噪声和畸变干扰对图像中目标区域的影响,突出图像中符合预定特征的目标区域。低分辨率的声纳图像一般没有斑点噪声,因此多数图像滤波方法都不适用于此类图像的处理。
发明内容
本发明的目的就是为了解决上述背景技术存在的不足,提供一种能够将低分辨率声纳图像中的小目标清晰提取出来的低分辨率高频声纳图像中目标亮点提取方法。
本发明采用的技术方案是:一种低分辨率高频声纳图像中目标亮点提取方法,包括以下步骤:
(1)、确定声纳图像中目标亮点在距离维和角度维上的最大尺寸;
(2)、计算声纳图像中每个像素点oi,j的信号强度
Figure BDA0000446267150000023
:在像素点oi,j左右宽度为角度维上最大尺寸的范围内,计算n组连续m个像素点的图像幅度值的平均值,取n个平均值中的最大值作为像素点oi,j的信号强度;
(3)、计算声纳图像中每个像素点oi,j的背景强度
Figure BDA0000446267150000024
:以声纳图像中目标亮点在距离维和角度维上的最大尺寸为间隔,取像素点oi,j周围8个角度像素点对应图像幅度值的均值作为像素点oi,j的背景强度;
(4)、逐点计算声纳图像中的每个像素点的信混比;
(5)、根据信混比大小确定每个像素点对应的声纳图像背景;
(6)、将各像素点的原始声纳图像与声纳图像背景相减,得到各像素点的目标亮点图像,对声纳图像进行逐点扫描,得到整体目标亮点图像。
进一步地,所述步骤(1)中目标亮点在距离维上的最大尺寸L的算法为:
L = [ A Δd ]
式中:L为目标亮点在距离维上的最大尺寸,A为目标的最大尺寸,△d为声纳图像的距离分辨长度,[·]表示向下取整。
进一步地,所述步骤(1)中目标亮点在角度维上的最大尺寸W的算法为:
W = K { [ 180 A βπD ] + } ,
式中:W为目标亮点在角度维上的最大尺寸,A为目标的最大尺寸,β为相邻波束间隔角度,D为目标与声纳基阵之间的距离,[·]表示向下取整。
进一步地,所述步骤(2)中n的取值为n=2W-m+1,m的取值为m=[W/2],其中[·]表示向下取整。
进一步地,所述步骤(4)中信混比的算法为
Figure BDA0000446267150000025
更进一步地,所述步骤(5)中每个像素点对应的声纳图像背景IR(i,j)的算法为:
I R ( i , j ) = R O i , j S O i , j / R O i , j &GreaterEqual; Q I ( i , j ) S O i , j / R O i , j < Q ,
式中:IR(i,j)为每个像素点对应的声纳图像背景,I(i,j)为每个像素点对应的原始声纳图像,Q为信混比阈值。
本发明依据特定目标亮点的尺寸,结合信混比大小,有效将目标亮点从声纳图像背景中提取出来,为后续目标检测、跟踪、辨识等提供条件。该提取方法简单、有效,提取结果准确、适用性高,解决了低分辨率高频声纳图像中小目标检测、辨识难的问题。
附图说明
图1为本发明声纳图像目标亮点提取流程图。
图2为本发明声纳图像中像素点信号强度计算模板。
图3为本发明声纳图像中像素点背景强度计算模板。
图4为本发明实施例中的原始声纳图像。
图5为本发明实施例中提取的目标亮点图像。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细说明,便于清楚地了解本发明,但它们不对本发明构成限定。
如图1所示,本发明低分辨率高频声纳图像中目标亮点提取方法包括以下步骤:
S101:确定目标亮点在声纳图像中距离维和角度维上的最大尺寸。
根据声纳图像的距离分辨长度△d、相邻波束间隔角度β以及目标最大尺寸A和距离D,可以确定目标在声纳图像中距离维和角度维上的最大尺寸。
目标亮点在距离维上的最大尺寸L通过最大像素点数来计算,即
L = [ A &Delta;d ] - - - ( 1 )
目标亮点在角度维上的最大尺寸W通过最大像素点数来计算,即
W = K { [ 180 A &beta;&pi;D ] + } , - - - ( 2 )
式中,[·]表示向下取整,K表示声纳由波束数据转换为图像时的角度维插值比例,计算出来的L和W最终以像素点的个数来衡量。
当声纳图像的距离分辨长度△d不是足够小,且角度维的内插比例K比较大时,声纳图像中的目标亮点形状往往表现为横条状。
S102:根据目标亮点在角度维上最大尺寸,计算声纳图像中各像素点的信号强度。
在角度维上,目标亮点中心像素对应幅度值最大,往两边幅度值逐渐减小至背景像素幅度平均值水平,因此目标亮点中各像素点的信号强度可以表示为以像素点为中心,在像素点左右宽度为W的范围内,计算多组连续若干个像素点的图像幅度值的平均值,取多个平均值中的最大值作为像素点的信号强度。由于目标亮点在角度维的尺寸一般较大,为了把目标亮点完整地呈现出来,将一定范围内连续多个像素点幅度的最大平均值作为信号强度的估计值是比较合理的,像素点信号强度的具体计算方法如下。
如图2所示,以像素点oi,j为待考察的中心像素,在中心像素角度维左右宽度为W的范围内,计算n组连续m个像素点的图像幅度值的平均值,取n个平均值中的最大值作为像素点的信号强度,左右宽度W的范围均不包含本身像素点oi,j,则该范围内的像素点的个数为2W+1,故n的取值为n=2W-m+1,即像素点oi,j的信号强度定义如下:
S O i , j = max 0 &le; p &le; 2 W - m + 1 1 m &Sigma; k = j - W + p j - W + p + m - 1 I ( i , k ) - - - ( 3 )
式中,m=[W/2],[·]表示向下取整,I(i,k)表示像素点oi,k的图像幅度值,以m个像素点的宽度范围为一个窗口,窗口每左移或右移一个像素点距离,就能计算出一个平均值,故p为连续m个像素点的宽度移动范围,
S103:根据目标亮点在距离维和角度维上的最大尺寸,计算声纳图像中各像素点的背景强度。
根据步骤S101计算得到目标亮点的距离维和角度维最大像素点数分别为L和W,目标背景强度的计算方法如下:
如图3所示,像素点oi,j为待考察的中心像素,将中心像素周围8个方向上的B1、B2、…、B8八个点作为目标背景强度的取样点,取样点像素与中心像素在角度维和距离维上的间隔分别为W和L,该范围同样不包含本身像素点,则待考察像素点oi,j的背景强度
Figure BDA0000446267150000053
定义为:像素点oi,j周围8个角度像素点对应图像幅度值的均值,以公式表示即为
R O i , j = 1 8 &Sigma; k = 1 8 I B k - - - ( 4 )
S104:逐点计算声纳图像中每个像素点的信混比,并根据信混比大小确定像素点的声纳图像背景。
以声纳图像中的每个像素点oi,j为待考察的中心像素,根据S102步的信号强度计算方法计算当前像素点的信号强度
Figure BDA0000446267150000054
,根据S103步的背景强度计算方法计算当前像素点的背景强度
Figure BDA0000446267150000055
,则每个像素点的信混比定义为
Figure BDA0000446267150000056
声纳图像中像素点对应的声纳图像背景IR(i,j)定义如下:
I R ( i , j ) = R O i , j S O i , j / R O i , j &GreaterEqual; Q I ( i , j ) S O i , j / R O i , j < Q - - - ( 5 )
式中,IR(i,j)为声纳图像中各像素点对应的声纳图像背景,I(i,j)为声纳图像中各像素点的原始声纳图像,Q为信混比阈值,取值一般在3~4之间。信混比阈值越大,提取的目标亮点越少,反之提取目标亮点越多。
S105:将声纳图像中各像素点的原始声纳图像I(i,j)与声纳图像背景IR(i,j)相减,得到各像素点的目标亮点图像IT(i,j),各像素点的目标亮点图像IT(i,j)定义如下:IT(i,j)=I(i,j)-IR(i,j)   (6)
最后对声纳图像进行逐点扫描,得到整体目标亮点图像IT
实施例1:据声纳基阵约300米处有一圆柱体小目标,最大尺度为2m。根据式(1)和(2),目标亮点在距离维和角度维上的最大尺度分别为4和11,信混比阈值取为3.2。按照本发明所述方法从如图4所示的原声纳图像中提取的目标亮点图像如图5所示,图中共提取出3个疑似目标亮点,真实目标亮点根据连续多幅图像来进一步辨识。
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (6)

1.一种低分辨率高频声纳图像中目标亮点提取方法,其特征在于,包括以下步骤:
(1)、确定声纳图像中目标亮点在距离维和角度维上的最大尺寸;
(2)、计算声纳图像中每个像素点oi,j的信号强度
Figure FDA0000446267140000013
:在像素点oi,j左右宽度为角度维上最大尺寸的范围内,计算n组连续m个像素点的图像幅度值的平均值,取n个平均值中的最大值作为像素点oi,j的信号强度;
(3)、计算声纳图像中每个像素点oi,j的背景强度
Figure FDA0000446267140000014
:以声纳图像中目标亮点在距离维和角度维上的最大尺寸为间隔,取像素点oi,j周围8个角度像素点对应图像幅度值的均值作为像素点oi,j的背景强度;
(4)、逐点计算声纳图像中的每个像素点的信混比;
(5)、根据信混比大小确定每个像素点对应的声纳图像背景;
(6)、将各像素点的原始声纳图像与声纳图像背景相减,得到各像素点的目标亮点图像,对声纳图像进行逐点扫描,得到整体目标亮点图像。
2.根据权利要求1所述低分辨率高频声纳图像中目标亮点提取方法,其特征在于:所述步骤(1)中目标亮点在距离维上的最大尺寸L的算法为:
L = [ A &Delta;d ]
式中:L为目标亮点在距离维上的最大尺寸,A为目标的最大尺寸,△d为声纳图像的距离分辨长度,[·]表示向下取整。
3.根据权利要求1所述低分辨率高频声纳图像中目标亮点提取方法,其特征在于:所述步骤(1)中目标亮点在角度维上的最大尺寸W的算法为:
W = K { [ 180 A &beta;&pi;D ] + } ,
式中:W为目标亮点在角度维上的最大尺寸,A为目标的最大尺寸,β为相邻波束间隔角度,D为目标与声纳基阵之间的距离,[·]表示向下取整。
4.根据权利要求1所述低分辨率高频声纳图像中目标亮点提取方法,其特征在于:所述步骤(2)中n的取值为n=2W-m+1,m的取值为m=[W/2],其中[·]表示向下取整。
5.根据权利要求1所述低分辨率高频声纳图像中目标亮点提取方法,其特征在于:所述步骤(4)中信混比的算法为
Figure FDA0000446267140000022
6.根据权利要求1所述低分辨率高频声纳图像中目标亮点提取方法,其特征在于:所述步骤(5)中每个像素点对应的声纳图像背景IR(i,j)的算法为:
I R ( i , j ) = R O i , j S O i , j / R O i , j &GreaterEqual; Q I ( i , j ) S O i , j / R O i , j < Q ,
式中:IR(i,j)为每个像素点对应的声纳图像背景,I(i,j)为每个像素点对应的原始声纳图像,Q为信混比阈值。
CN201310726052.2A 2013-12-25 2013-12-25 低分辨率高频声纳图像中目标亮点提取方法 Expired - Fee Related CN103714542B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310726052.2A CN103714542B (zh) 2013-12-25 2013-12-25 低分辨率高频声纳图像中目标亮点提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310726052.2A CN103714542B (zh) 2013-12-25 2013-12-25 低分辨率高频声纳图像中目标亮点提取方法

Publications (2)

Publication Number Publication Date
CN103714542A true CN103714542A (zh) 2014-04-09
CN103714542B CN103714542B (zh) 2016-09-14

Family

ID=50407485

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310726052.2A Expired - Fee Related CN103714542B (zh) 2013-12-25 2013-12-25 低分辨率高频声纳图像中目标亮点提取方法

Country Status (1)

Country Link
CN (1) CN103714542B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111123274A (zh) * 2019-12-27 2020-05-08 苏州联视泰电子信息技术有限公司 一种水下声呐成像系统目标检测方法
CN111325696A (zh) * 2020-03-03 2020-06-23 杭州瑞利海洋装备有限公司 一种基于正态分布区间估计的水声图像混响抑制方法
CN113538504A (zh) * 2020-04-16 2021-10-22 杭州海康威视数字技术股份有限公司 目标匀速移动场景中的跟踪目标显示方法和电子设备
CN115150390A (zh) * 2022-06-27 2022-10-04 山东信通电子股份有限公司 一种图像显示方法、装置、设备及介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010112426A1 (en) * 2009-04-03 2010-10-07 Tele-Rilevamento Europa - T.R.E. S.R.L. Process for identifying statistically homogeneous pixels in sar images acquired on the same area
CN102034109A (zh) * 2009-12-08 2011-04-27 中国科学院声学研究所 一种基于统计特性的合成孔径声纳目标检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010112426A1 (en) * 2009-04-03 2010-10-07 Tele-Rilevamento Europa - T.R.E. S.R.L. Process for identifying statistically homogeneous pixels in sar images acquired on the same area
CN102034109A (zh) * 2009-12-08 2011-04-27 中国科学院声学研究所 一种基于统计特性的合成孔径声纳目标检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王兴梅: "水下声呐图像目标分割方法的研究及应用", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
王超 等: "一种探雷声纳图像分割方法及应用", 《指挥控制与仿真》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111123274A (zh) * 2019-12-27 2020-05-08 苏州联视泰电子信息技术有限公司 一种水下声呐成像系统目标检测方法
CN111123274B (zh) * 2019-12-27 2021-12-28 苏州联视泰电子信息技术有限公司 一种水下声呐成像系统目标检测方法
CN111325696A (zh) * 2020-03-03 2020-06-23 杭州瑞利海洋装备有限公司 一种基于正态分布区间估计的水声图像混响抑制方法
CN111325696B (zh) * 2020-03-03 2023-07-07 杭州瑞利海洋装备有限公司 一种基于正态分布区间估计的水声图像混响抑制方法
CN113538504A (zh) * 2020-04-16 2021-10-22 杭州海康威视数字技术股份有限公司 目标匀速移动场景中的跟踪目标显示方法和电子设备
CN113538504B (zh) * 2020-04-16 2023-07-14 杭州海康威视数字技术股份有限公司 目标匀速移动场景中的跟踪目标显示方法和电子设备
CN115150390A (zh) * 2022-06-27 2022-10-04 山东信通电子股份有限公司 一种图像显示方法、装置、设备及介质
CN115150390B (zh) * 2022-06-27 2024-04-09 山东信通电子股份有限公司 一种图像显示方法、装置、设备及介质

Also Published As

Publication number Publication date
CN103714542B (zh) 2016-09-14

Similar Documents

Publication Publication Date Title
CN109633633B (zh) 一种基于分段分类增强处理的生命信号增强方法
CN104851097B (zh) 基于目标形状与阴影辅助的多通道sar‑gmti方法
CN102596050B (zh) 超声波成像装置及超声波成像方法
CN109143195B (zh) 一种基于全kl散度的雷达目标检测方法
CN104834915B (zh) 一种复杂云天背景下小红外目标检测方法
CN105319537A (zh) 基于空间相关性的航海雷达同频干扰抑制方法
CN104268505A (zh) 基于机器视觉的布匹疵点自动检测识别装置及方法
CN103197302B (zh) 一种适用于穿墙雷达成像的目标位置提取方法
CN103714542B (zh) 低分辨率高频声纳图像中目标亮点提取方法
CN103839257B (zh) 一种广义高斯k&i的sar图像变化检测方法
CN112541481A (zh) 一种基于深度学习的对海探测雷达目标检测方法
CN106707258A (zh) 一种非高斯背景下微动目标多参数估计方法
CN109188430B (zh) 一种基于地面监视雷达系统的目标提取方法
CN106156758B (zh) 一种sar海岸图像中海岸线提取方法
CN107942329A (zh) 机动平台单通道sar对海面舰船目标检测方法
CN104268877A (zh) 一种红外图像海天线自适应检测方法
CN104569923B (zh) 基于速度约束的Hough变换快速航迹起始方法
CN109100697B (zh) 一种基于地面监视雷达系统的目标凝聚方法
CN103500453A (zh) 基于伽玛分布和邻域信息的sar图像显著性区域检测方法
CN111323757B (zh) 一种航海雷达目标检测方法及装置
CN103473779B (zh) 图像中条纹干扰的检测方法和装置
CN110837079A (zh) 一种基于雷达的目标检测方法及装置
CN107369163A (zh) 一种基于最佳熵双阈值分割的快速sar图像目标检测方法
CN106780545A (zh) 一种结合时空域的弱小目标能量累积增强方法
CN109828275A (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160914

Termination date: 20181225