CN102393520A - 基于目标回波多普勒特性的声纳运动目标成像方法 - Google Patents

基于目标回波多普勒特性的声纳运动目标成像方法 Download PDF

Info

Publication number
CN102393520A
CN102393520A CN2011102877777A CN201110287777A CN102393520A CN 102393520 A CN102393520 A CN 102393520A CN 2011102877777 A CN2011102877777 A CN 2011102877777A CN 201110287777 A CN201110287777 A CN 201110287777A CN 102393520 A CN102393520 A CN 102393520A
Authority
CN
China
Prior art keywords
frequency
signal
data
meant
fft
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
CN2011102877777A
Other languages
English (en)
Other versions
CN102393520B (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.)
Jinzhou Northern Navigation Instrument Co ltd
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 CN 201110287777 priority Critical patent/CN102393520B/zh
Publication of CN102393520A publication Critical patent/CN102393520A/zh
Application granted granted Critical
Publication of CN102393520B publication Critical patent/CN102393520B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

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

Abstract

本发明提供的是一种基于目标回波多普勒特性的声纳运动目标成像方法。对接收基阵中各个基元的回波数据进行带通滤波;对各个频点信号进行带通滤波,获得到各频点的带通数据;对各频点带通滤波后的数据进行正交解调,获得实信号的解析信号形式;对各频点正交解调后的回波数据进行切比雪夫加权和聚焦操作;利用FFT波束形成获得各频点极坐标系下的图像数据;对各个频点FFT波束形成后的图像数据进行插值合并,并利用扇形变换,将极坐标系下的图像转成符合人们视觉习惯的直角坐标系下的图像。本发明采用FFT波束形成,能够利用数字信号处理中的FFT算法实现,满足系统实时成像的要求。

Description

基于目标回波多普勒特性的声纳运动目标成像方法
技术领域
本发明涉及的是一种声纳成像方法,具体地说是基于回波多普勒特性的声纳成像方法。
背景技术
传统的二维声纳成像的方法是时延波束形成,这种方法不能把波束定义在任意方向上,如果想要把波束定义在一定的方向上,要求的ADC(即模拟信号到数字信号的转换器)的采样频率较高,需要的存储量会很大,在硬件实现上会很困难。针对时延波束形成以上的缺点,可以通过内插波束形成来避免高的采样率,如果想得到理想的内插效果,对内插滤波器的要求较高,计算量会很大。
发明内容
本发明的目的在于提供一种能够满足系统实时成像的要求的基于目标回波多普勒特性的声纳运动目标成像方法。
本发明的目的是这样实现的:
(1)对接收基阵中各个基元的回波数据进行带通滤波;
(2)对各个频点信号进行带通滤波,获得到各频点的带通数据;
(3)对各频点带通滤波后的数据进行正交解调,获得实信号的解析信号形式;
(4)对各频点正交解调后的回波数据进行切比雪夫加权和聚焦操作;
(5)利用FFT(即快速傅立叶变换)波束形成获得各频点极坐标系下的图像数据;
(6)对各个频点FFT波束形成后的图像数据进行插值合并,并利用扇形变换,将极坐标系下的图像转成符合人们视觉习惯的直角坐标系下的图像。
本发明还可以包括:
1、所述的对接收基阵中各个基元的回波数据进行带通滤波是指对接收基阵的各个基元的回波数据进行FIR(即有限冲击响应)或者是IIR(即无限冲击响应)带通滤波处理。带通滤波器的下限频率和上限频率要根据发射信号的中心频率、目标运动所引起的回波信号的多普勒频移及接收机的带宽来确定。
2、所述的对各个频点信号进行带通滤波是指高速运动目标的回波信号已经不是单频信号而是一个宽带信号,因此窄带FFT波束形成已经不适用,需要将信号的频率范围分成M段,分别用带通滤波器滤出每段的信号,之后对每段信号分别进行窄带的FFT波束形成。
3、所述的正交解调因为实信号的傅里叶变换具有共轭对称性,因此成像时在对称的方位上会出现目标的虚影,而正交解调实质上是信号的希尔伯特变换,获得信号的解析表达式,解析信号的频谱是单边的,并能够完全描述原始信号,不会丢失任何信息。
4、所述的对正交解调后的回波数据进行切比雪夫加权是指对各个阵元接收信号的幅度进行加权,目的在于改善阵的方向性,降低旁瓣对主瓣的影响。聚焦是指正前方多点聚焦。这种方法适合声纳成像范围不是很大的情况下,设定基阵的尺度为L,目标距离为r,信号波长为λ,当r≤2L2/λ时,目标位于基阵近场,此时回波信号是按照球面波传播的,若直接对回波信号进行FFT波束形成,将会造成成像效果的恶化;此时需要对不同基元接收信号的相位进行补偿,以符合远场平面波的相位要求才能进行FFT波束形成。
5、所述的FFT波束形成的特点是能够用数字信号处理中的FFT算法实现。极坐标系下的目标图像是指图像的纵轴是距离信息R,横轴是目标的方位角度信息θ。
6、所述的插值合并是指对于512点FFT波束形成,扇面的开角为-52.8°-53.1°,将各个频点FFT波束形成后的角度索引经过线性插值到-45°-45°,0.5°间隔的范围上,之后将插值后的各个频点的成像数据在该角度范围上进行对应相加,得到总的输出的图像数据;所述的扇形变换是指极坐标与直角坐标之间的坐标转换,为了和现实的世界建立一种映射,方便人们观察,需要把图像从极坐标系下转换到直角坐标系下。
为实现本发明的目的,需要对基元数据进行带通滤波(步骤1),带通滤波器的下限频率和上限频率要根据发射信号的中心频率、目标运动所引起的回波信号的多普勒频移及接收机的带宽来确定。例如当发射信号频率为f0,发射信号的过零点带宽为Δf,目标运动所造成的回波信号的多普勒频移为ξ(目标背向接收基阵运动时,ξ<0;目标相向接收基阵运动时,ξ>0),则目标回波信号的频率变为f0+ξ,这里设定目标是背离接收基阵运动的,则ξ<0,将各个基于的回波数据经过带通采样及量化处理后信号频率变为f0′,目标回波信号的频率变为f0′+|ξ|,|ξ|要满足|ξ|>Δf,设定接收机的带宽为BW,则带通滤波器的下限频率为f0′+Δf/2,上限频率为f0′+BW/2。
为实现本发明的目的,对各个频点信号进行带通滤波(步骤2),需要将信号的频率范围分成M段,分别用带通滤波器滤出每段的信号,带通滤波器的带宽可以取定为1kHz或2kHz,用Δfbp表示,则第i段带通滤波器的下限频率为f0′+Δf/2+(2*i-1)/2*Δfbp,第i段带通滤波器的上限频率为f0′+Δf/2+(2*i+1)/2*Δfbp
为实现本发明的目的,需要对基元数据进行正交解调(步骤3),所述的正交解调是指将带通滤波后的基元数据与两个正交本振序列cosω0n、sinω0n相乘,再通过数字低通输出信号。
为实现本发明的目的,需要对回波数据进行切比雪夫加权(步骤4)。目的在于改善阵的方向性,降低旁瓣对主瓣的影响。还需要对回波数据进行正前方多点聚焦。以基阵的正前方为参考点,对于各个阵元的回波信号补偿一个相位差,这个相位差就是各个阵元对来自焦点位置的回波信号与参考点竖直位置的声程差d,焦点位置是基阵的正前方的声场,一般是在成像距离等间隔地取,补偿的方式就是把回波信号乘以exp(2πj*d/λ)。
为实现本发明的目的,需要对回波数据进行FFT波束形成(步骤5)。利用数字信号处理中的FFT算法获得距离R和目标方位角度θ间关系的极坐标系下的目标图像。
为实现本发明的目的,需要对各个频点FFT波束形成后的图像数据进行插值合并,对极坐标系下的图像数据进行扇形变换(步骤6)。目的在于和现实的世界建立一种映射,方便人们观察。
本发明的技术方案的核心为:对水下运动目标的回波数据进行成像处理并抑制多普勒的影响;对接收基阵中各个基元的回波数据进行带通滤波,滤除频率为发射信号中心频率的回波信号;再对带通滤波后的数据进行针对各个频点的带通滤波,得到各频点的带通数据,对各频点的带通数据分别进行常规成像处理,即进行正交解调,获得实信号的解析信号形式,对信号进行加权,抑制旁瓣的影响,对近场的回波数据进行聚焦,改善FFT波束形成在近场成像效果恶化的现象,采用能够满足实时成像的要求的FFT快速波束形成方法;对各个频点波束形成后的数据进行插值合并;最后进行扇形变换,将FFT波束形成后的极坐标系下的图像转成符合人们视觉习惯的直角坐标系下的图像。
本发明的特点是当水中体积混响与界面混响较强时,目标的反射被强目标所掩盖,因此常规的成像处理方法很难从强混响和目标散射背景中分辨出运动目标,由于目标是运动的,因此目标体散射的回波信号与混响信号以及静止物体回波信号之间存在着一等的频差,本发明正是基于该物理基础上对混响背景进行抑制,进而突出目标。
附图说明
图1是基于回波多普勒特性的声纳成像方法总流程图;
图2是常规成像处理方法流程图;
图3是正前方单波束聚焦波束形成几何关系图;
图4a是极坐标变换示意图,图4b是直角坐标变换示意图;
图5是扇形变换中插值算法示意图。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
图1给出了基于回波多普勒特性的声纳成像方法流程图,基于回波多普勒特性的声纳成像方法包括以下步骤:(1)带通滤波;(2)各频点带通滤波;(3)各频点常规成像处理;(4)插值合并;(5)扇形变换。由图2可知常规成像处理包括以下步骤:(1)正交变换;(2)加权操作;(3)聚焦操作;(4)基于FFT的波束形成操作。经过带通采样和幅度量化后的各个基元的回波信号首先要经过带通滤波处理,带通滤波器的下限频率和上限频率要根据发射信号的中心频率、目标运动所引起的回波信号的多普勒频移及接收机的带宽来确定。例如当发射信号频率为f0,发射信号的过零点带宽为Δf,目标运动所造成的回波信号的多普勒频移为ξ(目标背向接收基阵运动时,ξ<0;目标相向接收基阵运动时,ξ>0),则目标回波信号的频率变为f0+ξ,这里设定目标是背离接收基阵运动的,则ξ<0,将各个基于的回波数据经过带通采样及量化处理后信号频率变为f0′,目标回波信号的频率变为f0′+|ξ|,|ξ|要满足|ξ|>Δf,设定接收机的带宽为BW,则带通滤波器的下限频率为f0′+Δf/2,上限频率为f0′+BW/2。
对各个频点信号进行带通滤波是指高速运动目标的回波信号已经不是单频信号而是一个宽带信号,可知信号的频率范围是f0′+Δf/2-f0′+BW/2,因此窄带FFT波束形成已经不适用,需要将该频率范围分成M段,分别用带通滤波器滤出每段的信号,带通滤波器的带宽可以取定为1kHz或2kHz,用Δfbp表示,则第i段带通滤波器的下限频率为f0′+Δf/2+(2*i-1)/2*Δfbp,第i段带通滤波器的上限频率为f0′+Δf/2+(2*i+1)/2*Δfbp
实信号x(t)的傅里叶变换具有共轭对称性,即X(f)=X*(-f),其正、负频率的幅度分量是对称的,相位分量是相反的,只需用其正频分量或负频分量就可完全描述x(t),不会丢失任何信息,实信号的解析表达式正好符合此特点,设实信号x(t)的解析表达式为z(t),则:
z(t)=x(t)+j*H[x(t)]                        (1)
其中,H[x(t)]表示x(t)的Hilbert变换。由Hilbert变换的定义可知x(t)的解析表达式为
z(t)=x(t)*cos(ω0t)+j*x(t)*sin(ω0t)          (2)
所以正交解调是只将带通滤波后的基元数据与两个正交本振序列cosω0n、sinω0n相乘,再通过数字低通输出信号。
由于接收基阵是等间隔的直线阵,最常用的幅度加权法是道夫-契比雪夫(Dolph-Chebyshev)加权。它可在要求的主旁瓣比下,获得等旁瓣级;在任意给定的旁瓣级下使主瓣宽度最窄;在给定主瓣宽下使旁瓣级最低。但是幅度加权得到的旁瓣级的降低是以主瓣宽度的增加为代价的,所以不能一味要求主旁瓣比,要在主旁瓣比和主瓣宽度之间做一个适当的折中。
声波在近场不是按照平面规律进行传播,根据经验,远场平面波近似在r≥L2/λ处开始有效,L为基阵尺寸,λ为信号波长。在近场使用FFT波束形成会带来严重的波束性能恶化。在近场,可认为声波是按照球面波规律传播,对于等间隔线阵,可以得到阵元m相对于阵中心由球面波而引起的相位差φm,在波束形成之前,对这个相位差φm进行补偿,就可以得到较好的波束形成效果,φm就是聚焦因子。本发明采用的是如图2所示的正前方单个波束聚集方法,s是参考点,沿声轴方向距离阵中心r0处,近场波阵面为球面,可以利用几何关系计算得到基阵任意阵元与参考点的声程差φ1,φ2,φ3,·,φm,按照波束形成的基本原理:对各阵元接收的信号进行时延或相移补偿,使各阵元信号在正前方方向上形成同相相加,在波束形成时先把这些声程差补偿回来,φ1,φ2,φ3,·,φm被称作聚焦校正因子。
一个N元等间隔直线阵阵元间隔为d,当接收信号为单频或窄带信号时,基阵第i号阵元的输出信号可用复数表示为:
Figure BDA0000094328410000051
Figure BDA0000094328410000052
(i=0,1,·,N-1)          (3)
其中
Figure BDA0000094328410000053
为相邻阵元接收信号间的相位差,θ为信号入射方向与基阵法线之间的夹角,wi为基阵各阵元的幅度加权系数,ω为信号角频率,λ为信号波长。忽略时间因子,记阵元输出为:
Figure BDA0000094328410000054
如果在相邻阵元间插入相移βr,则基阵第i号阵元的输出信号变为:
Figure BDA0000094328410000055
阵输出为:
忽略时间因子,阵输出为:
Br = Σ i = 0 N - 1 x i e - ji β r - - - ( 7 )
如果取
Figure BDA0000094328410000062
(7)式就可以化为:
Br = Σ i = 0 N - 1 x i e - j 2 π N · i · r - - - ( 8 )
式(8)实际上是离散傅立叶变换的形式,因此计算一个等间隔直线阵各波束输出值就等价于计算各阵元的输出信号xi的离散傅立叶变换。
对各个频点波束FFT波束形成后的数据进行插值合并,对于512点FFT波束形成,扇面的开角为-52.8°-53.1°,将各个频点FFT波束形成后的角度索引经过线性插值到-45°-45°,0.5°间隔的范围上,之后将插值后的各个频点的成像数据在该角度范围上进行对应相加,得到总的输出的图像数据;FFT波束形成之后得到的是(θ,r)极坐标系下的图像,为了方便观察需要将极坐标系下的图像转换成直角坐标系下的图像,如图3所示,
r = x 2 + y 2 - - - ( 9 )
θ=arctg(y/x)                        (10)
由公式(9)、(10)可以找到直角坐标系下的像素点n(就是图2中的网格交叉点)在极坐标系下对应像素点m,如果m正好落在极坐标系下的栅格中,采用插值算法求出m点的值,此值就是像素点n的值;如果m落在极坐标系的栅格交叉点上,则极坐标下栅格交叉点的值就是像素点n的值。对于插值算法如图4所示,本发明采用的是最近邻域插值,它输出的像素值等于距离它映射道德位置最近的输入像素值。对于二维图像,该法是取代采样点周围4个相邻像素点中距离最近的1个邻点的灰度值作为改点的灰度值。插值公式:
f(i+u,j+v)=f(i,j)                (11)
其中i,j为非负整数,u,v为[-0.5,0.5)的浮点数,f(i,j)表示源图像的像素值,f(i+u,j+v)表示所求的像素点的像素值。

Claims (9)

1.一种基于目标回波多普勒特性的声纳运动目标成像方法,其特征是包括如下步骤:
(1)对接收基阵中各个基元的回波数据进行带通滤波;
(2)对各个频点信号进行带通滤波,获得到各频点的带通数据;
(3)对各频点带通滤波后的数据进行正交解调,获得实信号的解析信号形式;
(4)对各频点正交解调后的回波数据进行切比雪夫加权和聚焦操作;
(5)利用FFT波束形成获得各频点极坐标系下的图像数据;
(6)对各个频点FFT波束形成后的图像数据进行插值合并,并利用扇形变换,将极坐
标系下的图像转成符合人们视觉习惯的直角坐标系下的图像。
2.根据权利要求1所述的基于目标回波多普勒特性的声纳运动目标成像方法,其特征是:所述的对接收基阵中各个基元的回波数据进行带通滤波是指对接收基阵的各个基元的回波数据进行FIR或者是IIR带通滤波处理。
3.根据权利要求1或2所述的基于目标回波多普勒特性的声纳运动目标成像方法,其特征是:所述的对各个频点信号进行带通滤波是将信号的频率范围分成M段,分别用带通滤波器滤出每段的信号,之后对每段信号分别进行窄带的FFT波束形成。
4.根据权利要求1或2所述的基于目标回波多普勒特性的声纳运动目标成像方法,其特征是:所述的对正交解调后的回波数据进行切比雪夫加权是指对各个阵元接收信号的幅度进行加权,所述聚焦是指正前方多点聚焦。
5.根据权利要求3所述的基于目标回波多普勒特性的声纳运动目标成像方法,其特征是:所述的对正交解调后的回波数据进行切比雪夫加权是指对各个阵元接收信号的幅度进行加权,所述聚焦是指正前方多点聚焦。
6.根据权利要求1或2所述的基于目标回波多普勒特性的声纳运动目标成像方法,其特征是:所述的插值合并是指对于512点FFT波束形成,扇面的开角为-52.8°-53.1°,将各个频点FFT波束形成后的角度索引经过线性插值到-45°-45°,0.5°间隔的范围上,之后将插值后的各个频点的成像数据在该角度范围上进行对应相加,得到总的输出的图像数据;所述的扇形变换是指极坐标与直角坐标之间的坐标转换。
7.根据权利要求3所述的基于目标回波多普勒特性的声纳运动目标成像方法,其特征是:所述的插值合并是指对于512点FFT波束形成,扇面的开角为-52.8°-53.1°,将各个频点FFT波束形成后的角度索引经过线性插值到-45°-45°,0.5°间隔的范围上,之后将插值后的各个频点的成像数据在该角度范围上进行对应相加,得到总的输出的图像数据;所述的扇形变换是指极坐标与直角坐标之间的坐标转换。
8.根据权利要求4所述的基于目标回波多普勒特性的声纳运动目标成像方法,其特征是:所述的插值合并是指对于512点FFT波束形成,扇面的开角为-52.8°-53.1°,将各个频点FFT波束形成后的角度索引经过线性插值到-45°-45°,0.5°间隔的范围上,之后将插值后的各个频点的成像数据在该角度范围上进行对应相加,得到总的输出的图像数据;所述的扇形变换是指极坐标与直角坐标之间的坐标转换。
9.根据权利要求5所述的基于目标回波多普勒特性的声纳运动目标成像方法,其特征是:所述的插值合并是指对于512点FFT波束形成,扇面的开角为-52.8°-53.1°,将各个频点FFT波束形成后的角度索引经过线性插值到-45°-45°,0.5°间隔的范围上,之后将插值后的各个频点的成像数据在该角度范围上进行对应相加,得到总的输出的图像数据;所述的扇形变换是指极坐标与直角坐标之间的坐标转换。
CN 201110287777 2011-09-26 2011-09-26 基于目标回波多普勒特性的声纳运动目标成像方法 Active CN102393520B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110287777 CN102393520B (zh) 2011-09-26 2011-09-26 基于目标回波多普勒特性的声纳运动目标成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110287777 CN102393520B (zh) 2011-09-26 2011-09-26 基于目标回波多普勒特性的声纳运动目标成像方法

Publications (2)

Publication Number Publication Date
CN102393520A true CN102393520A (zh) 2012-03-28
CN102393520B CN102393520B (zh) 2013-07-31

Family

ID=45860878

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110287777 Active CN102393520B (zh) 2011-09-26 2011-09-26 基于目标回波多普勒特性的声纳运动目标成像方法

Country Status (1)

Country Link
CN (1) CN102393520B (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103954967A (zh) * 2014-05-19 2014-07-30 么彬 一种图像声纳用近场声纳图像快速成像方法
CN105520752A (zh) * 2014-10-21 2016-04-27 爱飞纽医疗机械贸易有限公司 波束成形设备、超声成像设备以及波束成形方法
CN105572656A (zh) * 2015-12-25 2016-05-11 中国船舶重工集团公司七五〇试验场 动态水声目标回波模拟的多普勒频移模拟器及模拟方法
CN108508446A (zh) * 2018-03-28 2018-09-07 青岛海洋地质研究所 基于海底冷泉成像数据的扇形变换方法
CN108519602A (zh) * 2018-03-28 2018-09-11 青岛海洋地质研究所 海底冷泉快速成像系统及其阵列信号处理方法
CN108683972A (zh) * 2014-10-20 2018-10-19 索尼公司 声音处理系统
CN108700651A (zh) * 2015-12-01 2018-10-23 超声成像 成像方法、实施所述方法的装置、计算机程序以及计算机可读存储介质
CN109218238A (zh) * 2018-10-18 2019-01-15 电子科技大学 一种基于希尔伯特变换的实信号多普勒频移方法
CN109283536A (zh) * 2018-09-01 2019-01-29 哈尔滨工程大学 一种多波束测深声呐水体成像波束形成算法
CN109459753A (zh) * 2017-10-16 2019-03-12 内蒙古工业大学 天气雷达数据坐标转换快速插值方法
CN111751890A (zh) * 2020-06-11 2020-10-09 上海微波技术研究所(中国电子科技集团公司第五十研究所) 一种能够抑制背景的地雷探测系统及方法
CN114089342A (zh) * 2022-01-19 2022-02-25 宁波博海深衡科技有限公司武汉分公司 单波束声呐水下测绘方法、系统、存储介质及终端
CN115902853A (zh) * 2023-03-09 2023-04-04 中国科学院深海科学与工程研究所 一种适用高速海底测绘的合成接收孔径聚焦波束形成方法
CN109709562B (zh) * 2018-11-08 2023-08-29 中国船舶重工集团公司七五0试验场 一种探测声纳的目标分辨数据处理方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101581785A (zh) * 2008-05-15 2009-11-18 中国科学院声学研究所 用于小型水下载体的三维前视声像声纳系统和方法
CN101825709A (zh) * 2009-12-08 2010-09-08 中国科学院声学研究所 一种水下高分辨率侧视声学成像系统与方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101581785A (zh) * 2008-05-15 2009-11-18 中国科学院声学研究所 用于小型水下载体的三维前视声像声纳系统和方法
CN101825709A (zh) * 2009-12-08 2010-09-08 中国科学院声学研究所 一种水下高分辨率侧视声学成像系统与方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
孙大军 等: "跳频脉冲信号在水下沉底小目标探测中的应用", 《哈尔滨工程大学学报》 *
郭元曦 等: "并行处理技术在高分辨率成像声纳中的应用", 《数据采集与处理》 *

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103954967B (zh) * 2014-05-19 2018-02-02 北京海卓同创科技有限公司 一种图像声纳用近场声纳图像快速成像方法
CN103954967A (zh) * 2014-05-19 2014-07-30 么彬 一种图像声纳用近场声纳图像快速成像方法
US11172292B2 (en) 2014-10-20 2021-11-09 Sony Corporation Voice processing system
CN108683972A (zh) * 2014-10-20 2018-10-19 索尼公司 声音处理系统
CN108683972B (zh) * 2014-10-20 2022-08-09 索尼公司 声音处理系统
CN105520752A (zh) * 2014-10-21 2016-04-27 爱飞纽医疗机械贸易有限公司 波束成形设备、超声成像设备以及波束成形方法
CN105520752B (zh) * 2014-10-21 2019-11-05 爱飞纽医疗机械贸易有限公司 波束成形设备、超声成像设备以及波束成形方法
CN108700651A (zh) * 2015-12-01 2018-10-23 超声成像 成像方法、实施所述方法的装置、计算机程序以及计算机可读存储介质
CN108700651B (zh) * 2015-12-01 2022-03-29 声科影像有限公司 成像方法、实施所述方法的装置以及计算机可读存储介质
CN105572656A (zh) * 2015-12-25 2016-05-11 中国船舶重工集团公司七五〇试验场 动态水声目标回波模拟的多普勒频移模拟器及模拟方法
CN109459753B (zh) * 2017-10-16 2022-10-11 内蒙古工业大学 天气雷达数据坐标转换快速插值方法
CN109459753A (zh) * 2017-10-16 2019-03-12 内蒙古工业大学 天气雷达数据坐标转换快速插值方法
CN108519602A (zh) * 2018-03-28 2018-09-11 青岛海洋地质研究所 海底冷泉快速成像系统及其阵列信号处理方法
CN108508446A (zh) * 2018-03-28 2018-09-07 青岛海洋地质研究所 基于海底冷泉成像数据的扇形变换方法
CN109283536A (zh) * 2018-09-01 2019-01-29 哈尔滨工程大学 一种多波束测深声呐水体成像波束形成算法
CN109283536B (zh) * 2018-09-01 2023-02-14 哈尔滨工程大学 一种多波束测深声呐水体成像波束形成方法
CN109218238B (zh) * 2018-10-18 2020-12-01 电子科技大学 一种基于希尔伯特变换的实信号多普勒频移方法
CN109218238A (zh) * 2018-10-18 2019-01-15 电子科技大学 一种基于希尔伯特变换的实信号多普勒频移方法
CN109709562B (zh) * 2018-11-08 2023-08-29 中国船舶重工集团公司七五0试验场 一种探测声纳的目标分辨数据处理方法
CN111751890A (zh) * 2020-06-11 2020-10-09 上海微波技术研究所(中国电子科技集团公司第五十研究所) 一种能够抑制背景的地雷探测系统及方法
CN114089342A (zh) * 2022-01-19 2022-02-25 宁波博海深衡科技有限公司武汉分公司 单波束声呐水下测绘方法、系统、存储介质及终端
CN115902853A (zh) * 2023-03-09 2023-04-04 中国科学院深海科学与工程研究所 一种适用高速海底测绘的合成接收孔径聚焦波束形成方法
CN115902853B (zh) * 2023-03-09 2023-05-12 中国科学院深海科学与工程研究所 一种适用高速海底测绘的合成接收孔径聚焦波束形成方法

Also Published As

Publication number Publication date
CN102393520B (zh) 2013-07-31

Similar Documents

Publication Publication Date Title
CN102393520B (zh) 基于目标回波多普勒特性的声纳运动目标成像方法
CN112505710B (zh) 一种多波束合成孔径声呐三维成像算法
US7914455B2 (en) Multi stage beamforming
CN110426707B (zh) 基于轨道角动量的涡旋sar成像方法及成像系统
CN101995574B (zh) 一种近场聚焦波束形成定位法
US7450470B2 (en) High resolution images from reflected wave energy
JPH10506800A (ja) 超音波画像における調節可能な周波数走査方法および装置
CN106680825B (zh) 一种声学阵列成像系统与方法
CN103018739B (zh) 一种校正多通道幅相误差的微波三维成像方法
CN110412587B (zh) 一种基于解卷积的下视合成孔径三维成像方法及系统
CN104597447B (zh) 一种子孔径SAR大斜视改进Omega‑K成像方法
CN102323581B (zh) 斜视聚束合成孔径雷达的成像方法
CN109298422B (zh) 一种用于穿墙雷达的合成孔径成像优化处理方法、装置
CN108008387B (zh) 一种机载阵列天线下视三维成像方法
CN110244300B (zh) 基于球体模型和fenlcs算法的弹载sar平飞段高分辨率成像方法
CN108132466A (zh) 一种机载阵列天线下视三维成像方法和系统
CN104765033B (zh) 利用互相关函数抑制步进频率成像中距离旁瓣的方法
CN105741236A (zh) 超声系统图像广义消旁瓣方法
EP2317335B1 (en) Improved beamforming method for analysing signals received by a transducer arrray, and relative detection system
Ponnle et al. Suppression of grating lobe artifacts in ultrasound images formed from diverging transmitting beams by modulation of receiving beams
CN108761433B (zh) 一种使用mimo声纳差合阵处理的高分辨成像方法
CN111239747A (zh) 一种基于解卷积的声纳高分辨低旁瓣二维成像方法
JP4114838B2 (ja) 超音波映像システム
CN105467387B (zh) 一种用于大孔径mimo阵列的低运算量混和波束形成方法
Su et al. A novel ionospheric clutter mitigation method through time-frequency image processing based on ridgelet analysis

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
ASS Succession or assignment of patent right

Owner name: JINZHOU NORTH NAUTICAL INSTRUMENT CO., LTD.

Free format text: FORMER OWNER: HARBIN ENGINEERING UNIV.

Effective date: 20150716

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20150716

Address after: 121219, No. 1, the Yellow Sea Avenue, Songshan New District, Liaoning, Jinzhou

Patentee after: JINZHOU NORTHERN NAVIGATION INSTRUMENT CO.,LTD.

Address before: 150001 Heilongjiang, Nangang District, Nantong street,, Harbin Engineering University, Department of Intellectual Property Office

Patentee before: HARBIN ENGINEERING University

PP01 Preservation of patent right

Effective date of registration: 20190925

Granted publication date: 20130731

PP01 Preservation of patent right
PD01 Discharge of preservation of patent

Date of cancellation: 20220925

Granted publication date: 20130731

PD01 Discharge of preservation of patent