CN110412586B - 一种基于相控发射和聚束成像的水下目标探查方法 - Google Patents

一种基于相控发射和聚束成像的水下目标探查方法 Download PDF

Info

Publication number
CN110412586B
CN110412586B CN201910611111.9A CN201910611111A CN110412586B CN 110412586 B CN110412586 B CN 110412586B CN 201910611111 A CN201910611111 A CN 201910611111A CN 110412586 B CN110412586 B CN 110412586B
Authority
CN
China
Prior art keywords
sub
azimuth
aperture
distance
mapping
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.)
Active
Application number
CN201910611111.9A
Other languages
English (en)
Other versions
CN110412586A (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.)
Ningbo Bohai Shenheng Technology Co.,Ltd.
Original Assignee
Ningbo Bohai Shenheng Technology Co ltd
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 Ningbo Bohai Shenheng Technology Co ltd filed Critical Ningbo Bohai Shenheng Technology Co ltd
Priority to CN201910611111.9A priority Critical patent/CN110412586B/zh
Publication of CN110412586A publication Critical patent/CN110412586A/zh
Application granted granted Critical
Publication of CN110412586B publication Critical patent/CN110412586B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8997Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using synthetic aperture techniques

Abstract

本发明公开了一种基于相控发射和聚束成像的水下目标探查方法,该方法包括如下步骤:利用相控发射技术对待探查测绘带进行探测,利用俯仰向多发射阵列获取与多个子测绘带一一对应的回波信号;将回波信号从空间域分开并进行解距离模糊,得到不模糊的子测绘带;对不模糊的子测绘带的回波信号进行成像处理,分析存在可疑目标的图像得到待精细查证的子测绘带;对待精细查证的子测绘带进行滑动聚束合成孔径成像,以实现对水下目标的精细查证。

Description

一种基于相控发射和聚束成像的水下目标探查方法
技术领域
本发明属于水下探测领域,具体涉及一种基于相控发射和聚束成像的水下目标探查方法。
背景技术
侧扫声纳亦称“旁侧声纳”或“海底地貌仪”。利用回声测深原理探测海底地貌和水下物体的设备。其换能器阵装在船壳内或拖曳体中,走航时向侧下方发射扇形波束的声脉冲。波束平面垂直于航行方向,沿航线方向束宽很窄,以保持较高分辨率;垂直于航线方向束宽较宽,以保证一定的扫描宽度。工作时发射出的声波投射在海底的区域呈长条形,换能器阵接收来的自照射区各点的反向散射信号,经放大、处理和记录,在记录纸上显示出海底的图像。波束侧扫声纳是传统的单波束侧扫声纳的升级改进,它在沿方位向形成多个平行的波束,最大航速可提高到单波束侧扫的M倍,M为形成波束的个数。但它的方位分辨率会随着距离的增加变差。
合成孔径声呐(Synthetic Aperture Sonar,简称SAS)是一种水下高分辨率的成像设备,其基本原理是利用小孔径基阵的移动来合成虚拟大孔径;之后通过将孔径收集到回波数据以一定的方式高效地相干叠加即可得到方位向的高分辨率。根据小孔径基阵的移动方式的不同,它一般可分为条带式、聚束式、滑动聚束式等成像方式。其中,条带式合成孔径成像的探测效率相比聚束式合成孔径成像要高;而聚束式合成孔径成像通过不断改变波束指向,使波束持续指向同一成像区域,因此可以获得比条带模式更长的合成孔径长度,其方位向分辨率仅取决于波束方位向转角积累和声波波长,且不受方位向天线尺寸的限制,相比于波束指向固定的SAS系统,聚束模式能够获取更高的方位向分辨率;另外聚束式成像可以在一次航行过程中对目标进行多视角成像,有利于目标的检测和识别,然而聚束式成像其计算较复杂,且计算费时,并不适用于大范围的成像识别。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种基于相控发射和聚束成像的水下目标探查方法,其利用俯仰向多发射阵列获取与多个子测绘带一一对应的回波信号,分析回波信号得到不模糊的子测绘带并进行成像处理,得到待精细查证的子测绘带,对待精细查证的子测绘带进行滑动聚束合成孔径成像,以实现对水下目标的精细查证,并提高水下目标的探查效率。
为实现上述目的,按照本发明的一个方面,提供了一种基于相控发射和聚束成像的水下目标探查方法,该方法包括如下步骤:
利用相控发射技术对待探查测绘带进行探测,利用俯仰向多发射阵列获取与多个子测绘带一一对应的回波信号;将回波信号从空间域分开并进行解距离模糊,得到不模糊的子测绘带;
对不模糊的子测绘带的回波信号进行成像处理,分析存在可疑目标的图像得到待精细查证的子测绘带;
对待精细查证的子测绘带进行滑动聚束合成孔径成像,即子孔径回波信号经过斜视ωk算法处理变换到子孔径回波信号的二维频域,进行方位频域Scaling和方位逆傅里叶变换得到方位数据频域,根据子孔径划分时的重叠率进行子孔径拼接得到方位全孔径数据;对方位全孔径数据去旋转和脉压得到方位全孔径脉压数据,对方位全孔径脉压数据进行方位向逆傅里叶变换得到SAS图像,以实现对水下目标的精细查证。
作为本发明的进一步改进,俯仰向多发射阵列为俯仰方向设置的N个天线阵元,第n个天线阵元坐标为(X0,Yn,Z0),
Figure GDA0003024468760000021
其中,X0、Y0和Z0分别为俯仰阵列波束中心的航迹向、距离向和深度向的坐标分量,Yn为第n个天线阵元的距离向的坐标分量,d为天线阵元的阵元间距。
作为本发明的进一步改进,第n个天线阵元接收的回波信号为:
Figure GDA0003024468760000022
其中,A(θl)为第n个天线阵元的方向性函数,θl为与第l个子测绘带对应的俯仰向波束中心的下视角,tk为快时间,tm=kTr,k为整数,Tr为脉冲重复周期,L为子脉冲的总个数,ΔTl为第l个子脉冲的发射时延,Rnl(tm)为第n个天线阵元等效相位中心与第l个子测绘带散射点的瞬时斜距,Tpl为子脉冲的时间宽度,fc为发射载频,γ为调频率。
作为本发明的进一步改进,将瞬时斜距Rnl(tm)近似为:
Rnl(tm)≈Rl(tm)+dsinθl′[n-(N+1)/2]/λ
θl′为与第l个子测绘带散射点对应的下视角,Rl(tm)为俯仰阵列波束中心与第l个子测绘带的瞬时斜距,λ=c/fc为发射波长,
对基带回波信号进行距离向匹配滤波后得到:
Figure GDA0003024468760000031
其中,A′(θl)=A(θl)exp(-j2πfcΔTl);
定义矩阵向量Z(l)为:
Z(l)={exp[-j2πdsinθl′(1-N)/2λ],...,exp[-j2πdsinθl′(N-1)/2λ]}T
取空域滤波的权矢量Wl=[wl1,...,wln,...wlN]T
取子测绘带对应的下视角输出为1,其他模糊下视角输出为0,即
ZT(l)Wl=Hl,l=1…L
其中,Hl=[h1,...,hL]T,hl=1,hi≠l=0;
在tm=MTr时刻,满足:
Figure GDA0003024468760000032
求解权矢量Wl解模糊得到不模糊的子测绘带回波信号为:
Figure GDA0003024468760000033
作为本发明的进一步改进,对不模糊的子测绘带的回波信号进行成像处理的算法包括:CS、RMA和FFBP。
作为本发明的进一步改进,第i个子孔径回波信号表示为:
Figure GDA0003024468760000034
式中:i取值范围为0~N-1,R(ta;Rp)为声纳瞬时位置到任意点P的距离,Rp为任意点P的径向距离,ta和τ分别为方位向时间和距离向时间,Kr为发射信号调频斜率,σn为反射系数,Tsub代表方位向子孔径的时间长度,fc为发射载频。
作为本发明的进一步改进,利用斜视ωk算法得到第i个子孔径回波信号的二维频域表达式为:
Figure GDA0003024468760000041
其中
Figure GDA0003024468760000042
Figure GDA0003024468760000043
分别为距离向和方位向频域,Xp为点P方位向坐标,
Figure GDA0003024468760000044
为参考距离,
Figure GDA0003024468760000045
为点P方位的位置信息,exp(-j2πfaiTsub)为第i个子孔径方位的位置信息,fτ为距离向时间变量τ为对应的频域变量,fa为方位向时间变量ta对应的频域变量,fdc,i为第i个子孔径多普勒中心频率,Va为方位向的运动速度。
作为本发明的进一步改进,用于方位频域Scaling的方位频域Scaling因子为:
Figure GDA0003024468760000046
式中:r为距离轴,
Figure GDA0003024468760000047
为用于消除方位时间位移的参量,
Figure GDA0003024468760000048
为随距离单元变化的参量。
作为本发明的进一步改进,用于子孔径拼接的方位去旋转因子Hscaling
Hde-rotation=exp(jπKrotta_all 2+j2πfdc_allta_all)
式中:ta_all和fdc_all分别为全孔径方位时间和方位频率,j2πfdc_allta_all为消除多普勒频率偏移的参量。
作为本发明的进一步改进,有效的方位多普勒斜率为:
Keff(r)=Krot-Kscl(r)
Figure GDA0003024468760000049
Krot为旋转中心产生的多普勒斜率,Rrot是孔径到旋转中心最短距离,θs为旋转中心的斜视角,
Figure GDA00030244687600000410
Kscl(r)为方位频域Scaling产生的多普勒斜率;
用于脉压的脉压因子为:
Figure GDA00030244687600000411
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
本发明的一种基于相控发射和聚束成像的水下目标探查方法,其利用俯仰向多发射阵列获取与多个子测绘带一一对应的回波信号,分析回波信号得到不模糊的子测绘带并进行成像处理,得到待精细查证的子测绘带,对待精细查证的子测绘带进行滑动聚束合成孔径成像,以实现对水下目标的精细查证,并提高水下目标的探查效率。
本发明的一种基于相控发射和聚束成像的水下目标探查方法,其未锁定可疑目标时,采用远距离和大测绘带的探测模式,对俯仰向多发射阵的回波信号进行距离向解模糊,之后采用条带式合成孔径成像;探测到可疑目标后,将波束指向锁定到目标处,距离向不再进行解模糊处理,直接进行滑动聚束合成孔径成像,实现对水下目标的精细查证。
本发明的一种基于相控发射和聚束成像的水下目标探查方法,通过控制发射阵的方位波束照射区域在水底的移动速度,来控制方位向相干累积的时间和方位分辨率,其成像的面积比聚束SAS大,分辨率高于相同阵列尺寸条带SAS的分辨率,可在分辨率和成像面积间进行权衡。
附图说明
图1为本发明实施例的俯仰向发射阵列的示意图;
图2为本发明实施例的滑动聚束合成孔径成像的示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。下面结合具体实施方式对本发明进一步详细说明。
一种基于相控发射和聚束成像的水下目标探查方法,具体步骤包括:
利用相控发射技术对待探查测绘带进行探测,利用俯仰向多发射阵列获取与多个子测绘带一一对应的回波信号;将回波信号从空间域分开并进行解距离模糊,得到不模糊的子测绘带;
具体地,控制俯仰阵列波束的方向,使其在不同的子脉冲时间照射不同的子测绘带,通过调整参数使得不同的子测绘带的回波时间处于预设时间范围内;可利用航行器或水下拖体艏部的俯仰向多发射阵和舷侧的线阵进行大测绘带的多波束侧扫成像或条带式合成孔径成像;
图1为本发明实施例的俯仰向发射阵列的示意图。如图1所示,俯仰向多发射阵列为俯仰方向设置的N个天线阵元,第n个天线阵元坐标为(X0,Yn,Z0),
Figure GDA0003024468760000061
其中,X0、Y0和Z0分别为俯仰阵列波束中心的航迹向、距离向和深度向的坐标分量,Yn为第n个天线阵元的距离向的坐标分量,d为天线阵元的阵元间距;
脉冲时宽Tp分为L个子脉冲,每个子脉冲时间宽度为Tpl,不同子脉冲时间内利不同方向的俯仰阵列波束照射不同子测绘带,在第l个子脉冲时宽内,俯仰天线阵元利用波束形成发射信号可表示为
Figure GDA0003024468760000062
其中θl为与第l个子测绘带对应的俯仰向波束中心的下视角,λ=c/fc为发射波长,c为光速,fc为发射载频,ΔTl为第l个子脉冲的发射时延,ΔTl=Tp1+Tp2+…+TpL-1,t为全时间,tk为快时间,t=tk+mTr,Tr为脉冲重复周期,m为整数,γ为调频率;
第n个天线阵元接收的基带回波信号可表示为:
Figure GDA0003024468760000063
其中,
Figure GDA0003024468760000064
其为第n个天线阵元的方向性函数,tm=kTr,k为整数。Rnl(tm)为第n个天线阵元等效相位中心与第l个子测绘带散射点的瞬时斜距,
(xl,yl,zl)为第l个子测绘带的散射点坐标,则瞬时斜距Rnl(tm)为
Figure GDA0003024468760000065
其可以近似为:
Rnl(tm)≈Rl(tm)+dsinθl′[n-(N+1)/2]/λ
令H=Z0-zl
Figure GDA0003024468760000066
θl′为与第l个子测绘带散射点对应的下视角,Rl(tm)为俯仰阵列波束中心与第l个子测绘带的瞬时斜距;
对基带回波信号进行距离向匹配滤波后得到:
Figure GDA0003024468760000071
其中,A′(θl)=A(θl)exp(-j2πfcΔTl);
由于不同子测绘带上不同斜距的散射点会同时达到接收阵元,从而导致距离模糊,尽管不同斜距的散射点回波信号在快时间域上是模糊,但是可以利用各回波信号对应的下视角不同将其从空间域分开,利用俯仰数字波束形成解距离模糊,从而得到不模糊的子测绘带。
定义矩阵向量Z(l)为:
Z(l)={exp[-j2πdsinθl′(1-N)/2λ],...,exp[-j2πdsinθl′(N-1)/2λ]}T
取空域滤波的权矢量Wl=[wl1,...,wln,...wlN]T
取子测绘带对应的下视角输出为1,其他模糊下视角输出为0,即
ZT(l)Wl=Hl,l=1…L
其中Hl=[h1,...,hL]T,hl=1,hi≠l=0
在tm=MTr时刻,满足:
Figure GDA0003024468760000072
从而,求解权矢量Wl解模糊得到不模糊的子测绘带回波信号为:
Figure GDA0003024468760000073
对不模糊的子测绘带的回波信号进行成像处理,分析存在可疑目标的图像得到待精细查证的子测绘带;对待精细查证的子测绘带进行滑动聚束合成孔径成像,实现对水下目标进行精细查证。
利用CS、RMA、FFBP等聚焦成像算法对不模糊的子测绘带回波信号进行后续的成像处理;分析存在可疑目标的图像得到待精细查证的子测绘带;
对待精细查证的子测绘带进行滑动聚束合成孔径成像,即子孔径回波信号经过斜视ωk算法处理变换到子孔径回波信号的二维频域,进行方位频域Scaling和方位逆傅里叶变换得到方位数据频域,根据子孔径划分时的重叠率进行子孔径拼接得到方位全孔径数据;对方位全孔径数据去旋转得到方位全孔径脉压数据,对方位全孔径脉压数据进行方位向逆傅里叶变换得到SAS图像,以实现对水下目标的精细查证。
其中,第i个子孔径回波信号可以表示为:
Figure GDA0003024468760000081
式中:i取值范围为0~N-1,R(ta;Rp)为声纳瞬时位置到任意点P的距离,Rp为任意点P的径向距离,ta和τ分别为方位向时间和距离向时间,Kr为发射信号调频斜率,σn为反射系数,Tsub代表方位向子孔径的时间长度。
利用斜视ωk算法得到第i个子孔径回波信号的二维频域表达式为:
Figure GDA0003024468760000082
其中
Figure GDA0003024468760000083
Figure GDA0003024468760000084
分别为距离向和方位向频域,Xp为点P方位向坐标,
Figure GDA0003024468760000085
为参考距离。式中第一项
Figure GDA0003024468760000086
为点P方位的位置信息,第二项exp(-j2πfaiTsub)为第i个子孔径方位的位置信息。正是由于子孔径方位信息的不同,才能实现子孔径方位数据的拼接。SSi(fτ,fa;Rp)是经过ωk算法处理后得到的二维频域表达式,其中,fτ为距离向时间变量τ为对应的频域变量,fa为方位向时间变量ta对应的频域变量,fdc,i为第i个子孔径多普勒中心频率,Va为方位向的运动速度。
图2为本发明实施例的滑动聚束合成孔径成像的示意图。如图2所示,O和O′分别为孔径中心和场景中心点,Orot为旋转中心点,θs为斜视角,R0和Rrot分别为孔径到场景最短距离和孔径到旋转中心最短距离。载机速度为Va,在孔径路径起始点A处,波束中心指向场景边缘点A′,此时波束中心斜视角为θa,水下载体沿着孔径AB航行,声纳波束中心在水底以匀速Vg扫描,到孔径路径终止点B处,波束中心指向场景另一边缘点B′,此时波束中心斜视角为θb。滑动聚束SAS方位带宽随斜视角增大而变宽,但是系统的PRF通常并不能满足方位奈奎斯特采样,因此方位频谱会出现混叠现象。可以将全孔径划分为子孔径处理,消除方位频谱混叠,然后再将子孔径拼接实现全孔径高分辨率成像。
滑动聚束中的子孔径拼接处理是将斜视ωk算法聚焦后,探测待精细查证的子测绘带得到的子孔径方位向信号线性调频后得到方位数据,然后根据子孔径划分时的重叠率进行子孔径拼接。子孔径信号经过斜视ωk算法处理变换到子孔径距离多普勒域,乘以方位频域Scaling因子,恢复方位线性调频,方位频域Scaling因子为:
Figure GDA0003024468760000091
式中:r为距离轴,
Figure GDA0003024468760000092
为用于消除方位时间位移的参量,fdc,i为第i个子孔径多普勒中心频率,
Figure GDA0003024468760000093
为随距离单元变化的参量。
经过方位频域Scaling后,进行方位逆傅里叶变换得到方位数据域,然后根据子孔径划分时的重叠率进行子孔径拼接得到方位全孔径数据。利用方位去旋转因子进行全孔径方位去旋转,消除方位波束围绕旋转中心产生的调频,方位去旋转因子Hscaling
Hde-rotation=exp(jπKrotta_all 2+j2πfdc_allta_all)
式中:ta_all和fdc_all分别为全孔径方位时间和方位频率,j2πfdc_allta_all为消除多普勒频率偏移的参量。
此时有效的方位多普勒斜率为Keff(r)=Krot-Kscl(r),脉压因子为
Figure GDA0003024468760000094
Figure GDA0003024468760000095
Krot为旋转中心产生的多普勒斜率,其中Rrot是孔径到旋转中心最短距离,θs为旋转中心的斜视角,
Figure GDA0003024468760000096
Kscl(r)为方位频域Scaling产生的多普勒斜率;Keff(r)为有效的方位多普勒斜率。
本发明提供了一种利用俯仰向多发射技术和滑动聚束合成孔径成像结合的水下目标探查方法,其中,俯仰向多发射技术可以提高合成孔径声纳的距离向成像的测绘带宽,它将发射脉冲分为多个子脉冲,在不同子脉冲时间内,通过俯仰阵列天线波束形成控制波束照射不同场景,从而得到大测绘带场景;在没有锁定可疑目标时,采用远距离和大测绘带的探测模式,此时对俯仰向多发射阵的回波信号进行距离向解模糊,之后采用条带式合成孔径成像。探测到可疑目标后,将波束指向锁定到目标处,距离向不再进行解模糊处理,直接进行滑动聚束合成孔径成像,实现对水下目标的精细查证。
本发明的滑动聚束合成孔径成像是介于条带模式与聚束模式之间的一种工作模式,它通过控制发射阵的方位波束照射区域在水底的移动速度,来控制方位向相干累积的时间和方位分辨率,其成像的面积比聚束SAS大,分辨率高于相同阵列尺寸条带SAS的分辨率,可在分辨率和成像面积间进行权衡。
假设声纳波束足迹的移动速度为Vf,声纳载体的移动速度为Vp,则滑动聚束模式的方位分辨率比条带模式的方位分辨率高
Figure GDA0003024468760000101
倍(其中通过载体速度和声纳波束控制,使Vp>Vf);当Vf=0时,等效于聚束模式成像,对滑动聚束模式Vf>0,可见滑动聚束的成像区域将大于聚束模式。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种基于相控发射和聚束成像的水下目标探查方法,其特征在于,该方法包括如下步骤:
利用相控发射技术对待探查测绘带进行探测,利用俯仰向多发射阵列获取与多个子测绘带一一对应的回波信号;将回波信号从空间域分开并进行解距离模糊,得到不模糊的子测绘带;
对不模糊的子测绘带的回波信号进行成像处理,分析存在可疑目标的图像得到待精细查证的子测绘带;
对待精细查证的子测绘带进行滑动聚束合成孔径成像,即子孔径回波信号经过斜视ωk算法处理变换到子孔径回波信号的二维频域,进行方位频域Scaling和方位逆傅里叶变换得到方位数据频域,根据子孔径划分时的重叠率进行子孔径拼接得到方位全孔径数据;对方位全孔径数据去旋转和脉压得到方位全孔径脉压数据,对方位全孔径脉压数据进行方位向逆傅里叶变换得到SAS图像,以实现对水下目标的精细查证;
所述俯仰向多发射阵列为俯仰方向设置的N个天线阵元,第n个天线阵元坐标为(X0,Yn,Z0),
Figure FDA0003153667400000011
其中,X0、Y0和Z0分别为俯仰阵列波束中心的航迹向、距离向和深度向的坐标分量,Yn为第n个天线阵元的距离向的坐标分量,d为天线阵元的阵元间距;
第n个天线阵元接收的回波信号为:
Figure FDA0003153667400000012
其中,A(θl)为第n个天线阵元的方向性函数,θl为与第l个子测绘带对应的俯仰向波束中心的下视角,tk为快时间,tm=kTr,k为整数,Tr为脉冲重复周期,L为子脉冲的总个数,ΔTl为第l个子脉冲的发射时延,Rnl(tm)为第n个天线阵元等效相位中心与第l个子测绘带散射点的瞬时斜距,Tpl为子脉冲的时间宽度,fc为发射载频,γ为调频率,c为光速,t为全时间。
2.根据权利要求1所述的一种基于相控发射和聚束成像的水下目标探查方法,其特征在于,将瞬时斜距Rnl(tm)近似为:
Rnl(tm)≈Rl(tm)+dsinθl′[n-(N+1)/2]/λ
θl′为与第l个子测绘带散射点对应的下视角,Rl(tm)为俯仰阵列波束中心与第l个子测绘带的瞬时斜距,λ=c/fc为发射波长,
对基带回波信号进行距离向匹配滤波后得到:
Figure FDA0003153667400000021
其中,A′(θl)=A(θl)exp(-j2πfcΔTl);
定义矩阵向量Z(l)为:
Z(l)={exp[-j2πdsinθl′(1-N)/2λ],...,exp[-j2πdsinθl′(N-1)/2λ]}T
取空域滤波的权矢量Wl=[wl1,...,wln,...wlN]T
取子测绘带对应的下视角输出为1,其他模糊下视角输出为0,即
ZT(l)Wl=Hl,l=1…L
其中,Hl=[h1,...,hL]T,hl=1,hi≠l=0;
在tm=MTr时刻,满足:
Figure FDA0003153667400000022
求解权矢量Wl解模糊得到不模糊的子测绘带回波信号为:
Figure FDA0003153667400000023
i为子孔径回波信号,M为形成波束的个数。
3.根据权利要求1或2所述的一种基于相控发射和聚束成像的水下目标探查方法,其特征在于,所述对不模糊的子测绘带的回波信号进行成像处理的算法包括:
CS、RMA和FFBP。
4.根据权利要求1或2所述的一种基于相控发射和聚束成像的水下目标探查方法,其特征在于,第i个子孔径回波信号表示为:
Figure FDA0003153667400000031
式中:i取值范围为0~N-1,R(ta;Rp)为声纳瞬时位置到任意点P的距离,Rp为任意点P的径向距离,ta和τ分别为方位向时间变量和距离向时间变量,Kr为发射信号调频斜率,σn为反射系数,Tsub代表方位向子孔径的时间长度,fc为发射载频。
5.根据权利要求4所述的一种基于相控发射和聚束成像的水下目标探查方法,其特征在于,利用斜视ωk算法得到第i个子孔径回波信号的二维频域表达式为:
Figure FDA0003153667400000032
其中
Figure FDA0003153667400000033
Figure FDA0003153667400000034
分别为距离向和方位向频域,Xp为点P方位向坐标,
Figure FDA0003153667400000035
为参考距离,
Figure FDA0003153667400000036
为点P方位的位置信息,exp(-j2πfaiTsub)为第i个子孔径方位的位置信息,fτ为距离向时间变量τ为对应的频域变量,fa为方位向时间变量ta对应的频域变量,fdc,i为第i个子孔径多普勒中心频率,Va为方位向的运动速度,θs为斜视角,R0为孔径到场景最短距离。
6.根据权利要求5所述的一种基于相控发射和聚束成像的水下目标探查方法,其特征在于,用于方位频域Scaling的方位频域Scaling因子为:
Figure FDA0003153667400000037
式中:r为距离轴,
Figure FDA0003153667400000038
为用于消除方位时间位移的参量,
Figure FDA0003153667400000039
Figure FDA00031536674000000310
为随距离单元变化的参量,λ为发射波长。
7.根据权利要求6所述的一种基于相控发射和聚束成像的水下目标探查方法,其特征在于,用于子孔径拼接的方位去旋转因子Hscaling
Hde-rotation=exp(jπKrotta_all 2+j2πfdc_allta_all)
式中:ta_all和fdc_all分别为全孔径方位时间和方位频率,j2πfdc_allta_all为消除多普勒频率偏移的参量,Krot为旋转中心产生的多普勒斜率。
8.根据权利要求7所述的一种基于相控发射和聚束成像的水下目标探查方法,其特征在于,有效的方位多普勒斜率为:
Keff(r)=Krot-Kscl(r)
Figure FDA0003153667400000041
Krot为旋转中心产生的多普勒斜率,Rrot是孔径到旋转中心最短距离,θs为斜视角,
Figure FDA0003153667400000042
Kscl(r)为方位频域Scaling产生的多普勒斜率;
用于脉压的脉压因子为:
Figure FDA0003153667400000043
CN201910611111.9A 2019-07-08 2019-07-08 一种基于相控发射和聚束成像的水下目标探查方法 Active CN110412586B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910611111.9A CN110412586B (zh) 2019-07-08 2019-07-08 一种基于相控发射和聚束成像的水下目标探查方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910611111.9A CN110412586B (zh) 2019-07-08 2019-07-08 一种基于相控发射和聚束成像的水下目标探查方法

Publications (2)

Publication Number Publication Date
CN110412586A CN110412586A (zh) 2019-11-05
CN110412586B true CN110412586B (zh) 2021-08-20

Family

ID=68360626

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910611111.9A Active CN110412586B (zh) 2019-07-08 2019-07-08 一种基于相控发射和聚束成像的水下目标探查方法

Country Status (1)

Country Link
CN (1) CN110412586B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111965640B (zh) * 2020-07-01 2022-11-18 西安电子科技大学 一种大斜视中高轨sar二维波束扫描方法
CN112379375B (zh) * 2020-10-26 2023-05-02 中国科学院空天信息创新研究院 基于双天线sar频分加码分斜视两发两收成像方法及装置
CN114089342A (zh) * 2022-01-19 2022-02-25 宁波博海深衡科技有限公司武汉分公司 单波束声呐水下测绘方法、系统、存储介质及终端

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5626621B2 (ja) * 2010-03-19 2014-11-19 株式会社 日立産業制御ソリューションズ フェイズドアレイ合成開口ソナーシステム
CN103576157A (zh) * 2012-07-18 2014-02-12 中国科学院声学研究所 一种基于多维波形编码的合成孔径声纳成像方法及系统
CN104678393B (zh) * 2015-01-30 2017-01-18 南京航空航天大学 斜视滑动聚束sar的子孔径波数域成像方法
CN109975815B (zh) * 2019-03-22 2021-09-07 宁波博海深衡科技有限公司 一种水下目标多波束声纳探测系统及方法

Also Published As

Publication number Publication date
CN110412586A (zh) 2019-11-05

Similar Documents

Publication Publication Date Title
CN108398677B (zh) 三坐标连续波一维相扫无人机低空目标检测系统
CN110412586B (zh) 一种基于相控发射和聚束成像的水下目标探查方法
EP2909652B1 (en) Interferometric inverse synthetic aperture radar and method
CN109477892B (zh) 用于环境检测的车辆雷达
JPH045156B2 (zh)
CN107390181B (zh) 一种基于多波束扫描的雷达高分辨成像方法
CN109975815B (zh) 一种水下目标多波束声纳探测系统及方法
Hansen Synthetic aperture sonar technology review
CN114910915A (zh) 一种侧扫声呐的水下目标多模式成像方法
JP5634238B2 (ja) レーダ装置
US9075142B2 (en) Device and method for selecting signal, and radar apparatus
Baumgartner et al. Multi-channel SAR for ground moving target indication
US20200271782A1 (en) Aft-Looking Sonar
CN110456362B (zh) 一种基于脉冲对发射的目标声成像及测速方法与系统
Saeedi Feasibility study and conceptual design of missile-borne synthetic aperture radar
US6650274B1 (en) Radar imaging system and method
JP5035782B2 (ja) スプリットビーム方式合成開口レーダ
KR20190012069A (ko) 지상배열안테나 시스템 및 지상배열안테나 시스템에서 상공의 탐지 영역의 이미지를 획득하는 방법
Pirkani et al. Doppler Beam Sharpening for MIMO and Real Aperture Radars at mm-wave and Sub-THz Maritime Sensing
CN115015925A (zh) 基于改进匹配追踪的机载阵列雷达超分辨前视成像方法及装置
Fontana et al. 3D ISAR/SAR imaging using multichannel real data
Pan et al. Shallow-water wideband low-frequency synthetic aperture sonar for an autonomous underwater vehicle
Ponce et al. Multiple-input multiple-output circular SAR
Bai et al. DOA Estimation Using Shipborne HFSWR Based on Virtual Synthetic Array
GB2104753A (en) Radars

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
TA01 Transfer of patent application right

Effective date of registration: 20210524

Address after: 12-1-1, No.1, building 1, East District, Ningbo new material innovation center, Ningbo hi tech Zone, 315000, Zhejiang Province

Applicant after: Ningbo Bohai Shenheng Technology Co.,Ltd.

Address before: No.4, fozuling Third Road, Donghu New Technology Development Zone, Wuhan City, Hubei Province 430205

Applicant before: WUHAN YUANHAI BOCHUANG TECHNOLOGY Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant