CN103901422A - 一种水下目标回波几何亮点结构特征提取方法 - Google Patents

一种水下目标回波几何亮点结构特征提取方法 Download PDF

Info

Publication number
CN103901422A
CN103901422A CN201410108098.2A CN201410108098A CN103901422A CN 103901422 A CN103901422 A CN 103901422A CN 201410108098 A CN201410108098 A CN 201410108098A CN 103901422 A CN103901422 A CN 103901422A
Authority
CN
China
Prior art keywords
time
highlight
geometrical
target
target echo
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.)
Pending
Application number
CN201410108098.2A
Other languages
English (en)
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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201410108098.2A priority Critical patent/CN103901422A/zh
Publication of CN103901422A publication Critical patent/CN103901422A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/534Details of non-pulse systems
    • G01S7/536Extracting wanted echo signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/04Systems determining presence of a target
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/539Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明属于水声技术应用领域,特别涉及用于水下沉底或掩埋目标的探测与识别的一种水下目标回波几何亮点结构特征提取方法。获得目标回波的时频分布图像;构造与目标几何亮点时频分布相匹配的结构元;对目标回波信号的时频分布图像进行形态滤波开操作。本方法通过对目标回波信号的时频分布进行图像处理,在时频平面上抑制混响并提取目标的几何亮点结构。本发明的成果可以应用于水下目标探测识别系统中,提取稳定的目标回波信号特征,提高探测准确率。

Description

一种水下目标回波几何亮点结构特征提取方法
技术领域
本发明属于水声技术应用领域,特别涉及用于水下沉底或掩埋目标的探测与识别的一种水下目标回波几何亮点结构特征提取方法。
背景技术
水下目标探测与识别的关键技术之一是从目标的声纳回波信号中提取出稳定的目标特征,从而实现与海底沉积层、礁石和水下生物等目标的区分。根据对目标声散射的分析,无论声波以何种方式照射目标,目标棱角产生的几何散射总是存在。由几何散射形成的几何亮点结构可以反映目标的几何形状与尺度信息,是目标回波的固有特征。但是当目标处于沉底状态时,强烈的海底混响会淹没目标回波,因此如何有效的抑制混响并提取目标特征一直是制约水下沉底目标探测与识别的难点问题。
目前水下目标回波信号特征提取方法主要根据目标回波与混响不同的时频特性。其中,目标几何亮点的形成服从线性声学规律,当主动声呐发射LFM脉冲时,目标几何亮点具有规则的时频分布形式。而根据混响的点散射模型,形成混响的各散射点的幅度与相位都是随机的,导致混响的时频分布也是随机的。根据这一特点,Wigner-Ville分布、小波包、Hilbert-Huang变换与分数阶傅里叶变换等时频分析方法被用于提取目标回波信号特征。然而这些方法得到的目标回波与混响在时频平面上的混合分布,对混响的抑制能力较弱,特征形式对声波入射角、频率以及水底底质变化敏感,当目标回波被混响淹没时无法识别出目标几何亮点结构。
发明内容
本发明的目的在于提供一种提高探测准确率的水下目标回波几何亮点结构特征提取方法。
本发明的内容是这样实现的:
(1)获得目标回波的时频分布图像:
当主动声纳发射线性调频脉冲时,单个几何亮点在时频平面上的分布投影为有限长线段,令单个几何亮点的信号形式为x(t)=exp[j2π(f0t+mt2/2)],f0为信号起始频率,m为调频斜率,Wigner-Ville分布的时频分布为,
W x ( t , f ) = ∫ - ∞ ∞ e j 2 π ( f 0 + mt ) e - j 2 πfτ dτ = δ [ f - ( f 0 + mt ) ] ;
其中,t与f分别代表各个时频单元对应的时间与频率,两个几何亮点自项之间会存在交叉项,
Wcross(t,f)=2πδ(f-fm)exp(2πfdt);
其中,fm=(f1+f2)/2,fd=f1-f2,f1与f2分别是两个几何亮点在t时刻的频率;
(2)构造与目标几何亮点时频分布相匹配的结构元:
选择具有三角形包络的序列作为结构元,结构元的长度为目标回波几何亮点的时频分布长度的一半,对目标回波的WVD进行Radon变换,搜索目标几何亮点的实际调频斜率,x′与y′表示旋转后的坐标轴,则旋转角度为θ的Radon变换Rθ(x′)为,
R θ ( x ′ ) = ∫ - ∞ ∞ f ( x ′ cos θ - y ′ sin θ , x ′ sin θ + y ′ cos θ ) dy ′
当旋转角度为θmax时Rθ(x′)出现最大值,那么目标几何亮点的调频斜率即为cos(θmax),将结构元调整到此斜率下到与目标几何亮点时频分布的最佳匹配;
(3)对目标回波信号的时频分布图像进行形态滤波开操作:
使用结构元对获得目标回波的WVD二值图像进行形态学开操作。
本发明的有益效果在于:本方法通过对目标回波信号的时频分布进行图像处理,在时频平面上抑制混响并提取目标的几何亮点结构。本发明的成果可以应用于水下目标探测识别系统中,提取稳定的目标回波信号特征,提高探测准确率。
附图说明
图1为本发明的流程图。
图2为双几何亮点的时频分布图像。
图3为经过形态滤波的双几何亮点的时频分布图像。
具体实施方式
下面结合附图对本发明做进一步描述。
为了克服现有水下目标回波信号特征提取方法无法抑制混响干扰的问题,本发明提供了一种水下目标回波几何亮点结构提取方法,该方法通过对目标回波信号的时频分布进行图像处理,在时频平面上抑制混响并提取目标的几何亮点结构。本发明的成果可以应用于水下目标探测识别系统中,提取稳定的目标回波信号特征,提高探测准确率。本发明解决其技术问题所采用的技术方案包含以下流程:
1、获得目标回波信号的时频分布图像;
2、为了匹配目标几何亮点在时频平面上幅度鳍状分布特征,使用具有三角形包络的序列作为结构元;为了匹配目标几何亮点时频分布斜率特征,对目标回波的时频分布进行Radon变换,搜索对应Radon变换最大值旋转角度计算结构元的斜率;
3、对目标回波时频分布图像进行形态滤波开操作,消除混响在时频平面上造成的干扰。
本发明的一维几何回波时域信号在二维时频域中的几何分布特征作为图像特征进行特征提取,通过对目标回波信号的时频分布进行形态滤波,在时频平面上抑制混响并提取目标的几何亮点结构。设计了一种与目标回波时频分布的幅度特征、斜率特征相匹配的结构元。为了匹配目标几何亮点在时频平面上幅度鳍状分布特征,设计了具有三角形包络的序列作为结构元。为了匹配目标几何亮点时频分布斜率特征,对目标回波的时频分布进行Radon变换,搜索对应Radon变换最大值旋转角度计算结构元的斜率。对目标回波信号的时频分布图像进行形态滤波开操作以消除混响在时频平面上造成的干扰。
为了帮助对本发明内容的理解,下面结合附图及实施例对本发明的具体实施方式进行说明。
流程1:获得目标回波的时频分布图像。
根据目标回波的产生机理,当主动声纳发射线性调频脉冲时,单个几何亮点在时频平面上的分布投影为有限长线段,线段的斜率与发射信号调频斜率一致。令单个几何亮点的信号形式为x(t)=exp[j2π(f0t+mt2/2)],f0为信号起始频率,m为调频斜率,以Wigner-Ville分布(WVD)为例其时频分布为,
W x ( t , f ) = ∫ - ∞ ∞ e j 2 π ( f 0 + mt ) e - j 2 πfτ dτ = δ [ f - ( f 0 + mt ) ] - - - ( 1 )
其中,t与f分别代表WVD中各个时频单元对应的时间与频率。式(1)说明单个几何亮点的WVD为冲激谱。此外,由于实际信号长度一定是有限长的,因而其WVD呈鳍状。在WVD中,两个几何亮点自项之间会存在交叉项。
Wcross(t,f)=2πδ(f-fm)exp(2πfdt)  (2)
其中,fm=(f1+f2)/2,fd=f1-f2,f1与f2分别是两个几何亮点在t时刻的频率。式(2)说明交叉项的WVD呈周期起伏状。在本发明中,交叉项与混响及噪声的抑制原理一致,因此本实施例中只以抑制交叉项为例说明本发明的实施方法。双几何亮点的WVD示意图如附图2所示。
流程2:构造与目标几何亮点时频分布相匹配的结构元。
由于目标回波几何亮点的时频分布为鳍状,同时考虑工程实现的难度,选择具有三角形包络的序列作为结构元。序列长度需要考虑目标回波信号的信混比,混响较强时序列的长度也要加长,以达到更好的混响抑制效果。本实施例中结构元的长度为目标回波几何亮点的时频分布长度的一半。
结构元斜率决定了结构元与目标几何亮点时频分布的匹配程度,而几何亮点的调频斜率与发射信号基本一致,但考虑到脉冲展宽现象,实际的几何亮点调频斜率要略小于发射信号的调频斜率。对目标回波的WVD进行Radon变换,搜索目标几何亮点的实际调频斜率。二维图像f(x,y)的Radon变换是将坐标轴旋转一定角度后,沿纵轴方向对图像进行积分。设x′与y′表示旋转后的坐标轴,则旋转角度为θ的Radon变换Rθ(x′)定义为
R θ ( x ′ ) = ∫ - ∞ ∞ f ( x ′ cos θ - y ′ sin θ , x ′ sin θ + y ′ cos θ ) dy ′
如果当旋转角度为θmax时Rθ(x′)出现最大值,那么目标几何亮点的调频斜率即为cos(θmax),将结构元调整到此斜率下可以达到与目标几何亮点时频分布的最佳匹配。
流程3:对目标回波信号的时频分布图像进行形态滤波开操作。
在通过流程1获得目标回波的WVD二值图像的基础上,使用流程2设计的结构元对该图像进行形态学开操作。形态学开操作是由两个基本的形态学操作组成的,先进行的是腐蚀操作,然后进行的是膨胀操作。对目标回波时频分布图像进行腐蚀可以消除图像中与结构元形态特征不符的成分,如交叉项与混响的时频分布,保留几何亮点自项,但长度缩短,而膨胀会恢复信号自项长度。附图3所示是对附图2双几何亮点的时频分布图像进行形态滤波开操作的结果,由于腐蚀已经消除了混响与交叉项,因此经过形态滤波开操作后,时频平面上只存在有目标几何亮点,达到了本发明的目的。

Claims (1)

1.一种水下目标回波几何亮点结构特征提取方法,其特征在于:
(1)获得目标回波的时频分布图像:
当主动声纳发射线性调频脉冲时,单个几何亮点在时频平面上的分布投影为有限长线段,令单个几何亮点的信号形式为x(t)=exp[j2π(f0t+mt2/2)],f0为信号起始频率,m为调频斜率,Wigner-Ville分布的时频分布为,
W x ( t , f ) = ∫ - ∞ ∞ e j 2 π ( f 0 + mt ) e - j 2 πfτ dτ = δ [ f - ( f 0 + mt ) ] ;
其中,t与f分别代表各个时频单元对应的时间与频率,两个几何亮点自项之间会存在交叉项,
Wcross(t,f)=2πδ(f-fm)exp(2πfdt);
其中,fm=(f1+f2)/2,fd=f1-f2,f1与f2分别是两个几何亮点在t时刻的频率;
(2)构造与目标几何亮点时频分布相匹配的结构元:
选择具有三角形包络的序列作为结构元,结构元的长度为目标回波几何亮点的时频分布长度的一半,对目标回波的WVD进行Radon变换,搜索目标几何亮点的实际调频斜率,x′与y′表示旋转后的坐标轴,则旋转角度为θ的Radon变换Rθ(x′)为,
R θ ( x ′ ) = ∫ - ∞ ∞ f ( x ′ cos θ - y ′ sin θ , x ′ sin θ + y ′ cos θ ) dy ′
当旋转角度为θmax时Rθ(x′)出现最大值,那么目标几何亮点的调频斜率即为cos(θmax),将结构元调整到此斜率下到与目标几何亮点时频分布的最佳匹配;
(3)对目标回波信号的时频分布图像进行形态滤波开操作:
使用结构元对获得目标回波的WVD二值图像进行形态学开操作。
CN201410108098.2A 2014-03-21 2014-03-21 一种水下目标回波几何亮点结构特征提取方法 Pending CN103901422A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410108098.2A CN103901422A (zh) 2014-03-21 2014-03-21 一种水下目标回波几何亮点结构特征提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410108098.2A CN103901422A (zh) 2014-03-21 2014-03-21 一种水下目标回波几何亮点结构特征提取方法

Publications (1)

Publication Number Publication Date
CN103901422A true CN103901422A (zh) 2014-07-02

Family

ID=50992872

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410108098.2A Pending CN103901422A (zh) 2014-03-21 2014-03-21 一种水下目标回波几何亮点结构特征提取方法

Country Status (1)

Country Link
CN (1) CN103901422A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106886017A (zh) * 2017-01-11 2017-06-23 浙江大学 基于双频识别声呐的水下目标空间位置计算方法
CN107272005A (zh) * 2017-05-27 2017-10-20 西北工业大学 基于可靠声路径下目标回波到达时延和到达角度的主动定位方法
CN108508483A (zh) * 2018-06-07 2018-09-07 中国人民解放军海军工程大学 一种滩涂近地表掩埋目标的地震波探测方法和系统
CN108593767A (zh) * 2018-01-24 2018-09-28 天津大学 一种浅滩掩埋管道超声探测回波模型的建立方法
CN111323784A (zh) * 2020-03-25 2020-06-23 哈尔滨工程大学 一种近场分布式声源定位方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040042531A1 (en) * 2002-08-28 2004-03-04 Orhan Arikan System and method for CDMA communications
CN1588124A (zh) * 2004-08-20 2005-03-02 清华大学 利用多普勒变换的被动测速测距法及装置
US20090067286A1 (en) * 2007-09-12 2009-03-12 Schlumberger Technology Corporation Dispersion extraction for acoustic data using time frequency analysis
CN101900810A (zh) * 2010-07-15 2010-12-01 哈尔滨工程大学 多潜器载体探测端声纳信息融合方法
CN102305940A (zh) * 2011-05-24 2012-01-04 中国石油集团川庆钻探工程有限公司 流体因子提取方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040042531A1 (en) * 2002-08-28 2004-03-04 Orhan Arikan System and method for CDMA communications
CN1588124A (zh) * 2004-08-20 2005-03-02 清华大学 利用多普勒变换的被动测速测距法及装置
US20090067286A1 (en) * 2007-09-12 2009-03-12 Schlumberger Technology Corporation Dispersion extraction for acoustic data using time frequency analysis
CN101900810A (zh) * 2010-07-15 2010-12-01 哈尔滨工程大学 多潜器载体探测端声纳信息融合方法
CN102305940A (zh) * 2011-05-24 2012-01-04 中国石油集团川庆钻探工程有限公司 流体因子提取方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
何劲等: ""随机线性调频步进雷达波形设计及成像算法研究"", 《电子与信息学报》 *
尚海燕等: ""基于时频形态学滤波的能量积累检测"", 《电子与信息学报》 *
祝海华: ""基于时频分布的数字形态学的信号检测"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
赵学云等: ""基于Radon-WVD变换的编队目标架次识别"", 《电子与信息学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106886017A (zh) * 2017-01-11 2017-06-23 浙江大学 基于双频识别声呐的水下目标空间位置计算方法
CN106886017B (zh) * 2017-01-11 2020-02-07 浙江大学 基于双频识别声呐的水下目标空间位置计算方法
CN107272005A (zh) * 2017-05-27 2017-10-20 西北工业大学 基于可靠声路径下目标回波到达时延和到达角度的主动定位方法
CN107272005B (zh) * 2017-05-27 2020-04-03 西北工业大学 基于可靠声路径下目标回波到达时延和到达角度的主动定位方法
CN108593767A (zh) * 2018-01-24 2018-09-28 天津大学 一种浅滩掩埋管道超声探测回波模型的建立方法
CN108508483A (zh) * 2018-06-07 2018-09-07 中国人民解放军海军工程大学 一种滩涂近地表掩埋目标的地震波探测方法和系统
CN111323784A (zh) * 2020-03-25 2020-06-23 哈尔滨工程大学 一种近场分布式声源定位方法
CN111323784B (zh) * 2020-03-25 2023-01-03 哈尔滨工程大学 一种近场分布式声源定位方法

Similar Documents

Publication Publication Date Title
CN103901422A (zh) 一种水下目标回波几何亮点结构特征提取方法
CN111580048B (zh) 一种利用单矢量水听器的宽带声源深度估计方法
CN102981146A (zh) 一种单矢量水听器被动定位方法
CN107179535A (zh) 一种基于畸变拖曳阵的保真增强波束形成的方法
CN103529441B (zh) 一种被动合成孔径目标信号检测和分辨方法及系统
CN106019263B (zh) 基于多亮点回波模型的目标径向速度测量方法
CN103076604B (zh) 一种基于频散特征的低频水声脉冲信号距离的测量方法
JP2015148452A (ja) 合成開口レーダ装置
CN103076594A (zh) 一种基于互相关的水声脉冲信号双阵元定位的方法
CN102999473A (zh) 一种线性调频信号的检测与参数估计方法
CN103076590A (zh) 一种基于频率预估的水声脉冲信号的定位方法
CN105158734B (zh) 一种基于阵不变量的单矢量水听器被动定位方法
CN104820218A (zh) 一种基于频域自相关的浅海海底单参数反演方法
CN105629220A (zh) 一种基于单水听器的深海水声被动测距方法
CN104076404A (zh) 运用多通道相干抑制地磁背景噪声的磁异常探测方法
CN104678384A (zh) 一种波束域的声压差互相关谱分析水下目标速度估计方法
CN105204022A (zh) 海面风场的反演方法及装置
CN103487796B (zh) 一种利用水声信道统计不变特征实现被动测距的方法
CN105353340B (zh) 一种双层圆柱阵水下被动目标检测方法
CN103809170A (zh) 一种多普勒频偏修正方法及系统
CN108802725A (zh) 一种浅层穿透雷达合成孔径成像方法
CN104062663A (zh) 一种多波束海底浅地层剖面探测设备
CN103105603A (zh) 一种x波段测波雷达海流反演预处理方法
CN110346802A (zh) 基于计算水声信道参数的水下目标探测方法
Turgut et al. Measured depth-dependence of waveguide invariant in shallow water with a summer profile

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20140702

RJ01 Rejection of invention patent application after publication