CN116520335B - 一种多接收阵元合成孔径声呐波数域成像方法 - Google Patents

一种多接收阵元合成孔径声呐波数域成像方法 Download PDF

Info

Publication number
CN116520335B
CN116520335B CN202310783214.XA CN202310783214A CN116520335B CN 116520335 B CN116520335 B CN 116520335B CN 202310783214 A CN202310783214 A CN 202310783214A CN 116520335 B CN116520335 B CN 116520335B
Authority
CN
China
Prior art keywords
distance
array element
representing
frequency domain
receiving array
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
CN202310783214.XA
Other languages
English (en)
Other versions
CN116520335A (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.)
Haiying Deep Sea Technology Co ltd
Original Assignee
Haiying Deep Sea 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 Haiying Deep Sea Technology Co ltd filed Critical Haiying Deep Sea Technology Co ltd
Priority to CN202310783214.XA priority Critical patent/CN116520335B/zh
Publication of CN116520335A publication Critical patent/CN116520335A/zh
Application granted granted Critical
Publication of CN116520335B publication Critical patent/CN116520335B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (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

一种多接收阵元合成孔径声呐波数域成像方法
技术领域
本发明涉及图像处理技术领域,具体地说,涉及一种多接收阵元合成孔径声呐波数域成像方法。
背景技术
目前,合成孔径声呐技术是将常规声呐的方位分辨率提高一个数量级的有效技术。其原理是利用小孔径基阵的匀速直线运动在方位向上虚拟合成一个大孔径基阵,由于这个虚拟的孔径正相关于测绘距离,所以方位分辨率具有与作用距离和工作频率无关的特性。合成孔径成像是合成孔径声呐信号处理的核心内容,是获取高分辨图像和后期干涉信号处理的关键因素。根据线性系统的观点,可以将合成孔径声呐等效为一个线性系统,系统的输入是散射场的二维声散射函数,系统的输出就是合成孔径声呐系统所接收到的目标回波信号,将合成孔径声呐系统中这个接收信号的过程看成一个正向过程,那么对场景目标声散射特性的求取就是一个逆向重构过程,也就是说合成孔径成像的任务就是在合成孔径声呐系统接收到的回波和系统转移函数基础上重构成像场景的二维声散射场。
多接收阵元合成孔径声呐是当前能够同时提高方位分辨率与距离测绘带的有效系统,其在方位向上包含一个发射阵元和多个接收阵元,由于其精确双程斜距历程中包含两个含有根号形式的单程斜距历程,所以利用相位驻留原理不能得到精确且解析的二维频域系统函数,因此必须首先对精确二维频域系统函数采用合理的近似,然后再基于近似的二维频域系统函数设计多接收阵元合成孔径声呐快速成像算法。如Callow H J,Hayes MP,Gough P T.等人在《Electronics Letters》发表的《Wavenumber domainreconstruction of SAR/SAS imagery using single transmitter and multiple-receiver geometry》(译为《利用单发射机和多接收机几何结构对SAR/SAS图像的波数域重建》)中提出一种图像的波数域成像方法。再如黄攀,唐劲松,钟何平等人在《武汉大学学报(信息科学版)》发表的《基于有理函数曲面拟合的InSAS复图像配准新方法》和《干涉合成孔径声呐复图像配准分段曲面拟合法》中,均提出了一种基于曲面拟合的InSAS复图像配准方法,该方法具有拟合精度高,计算量小的优点。
然而往往通过某些近似得到的二维频域系统函数非常复杂,这进一步给后期成像算法的设计增加了难度,为此为简化成像算法设计,同时能够提高成像算法的成像性能,拟对二维频域系统函数进行线性化处理,我们提出了一种多接收阵元合成孔径声呐波数域成像方法。
发明内容
本发明的目的在于提供一种多接收阵元合成孔径声呐波数域成像方法,以解决上述背景技术中提出的问题。
为实现上述技术问题的解决,本发明的目的之一在于,提供了一种多接收阵元合成孔径声呐波数域成像方法,包括如下步骤:
S1、首先计算包含距离空变和距离空不变的二维频域系统函数;
S2、针对每个接收阵元的数据在二维频域进行一致相位补偿;
S3、针对每个接收阵元的数据在二维频域进行距离向坐标变换操作;
S4、在二维频域对所有接收阵元数据融合后进行距离向傅里叶逆变换;
S5、针对融合后的数据,在距离-多普勒域补偿方位走动误差,方位傅里叶变换后得到最终的高分辨结果。
作为本技术方案的进一步改进,所述S1中,计算包含距离空变和距离空不变的二维频域系统函数,其计算公式如下:
其中,下标表示第个接收阵元和发射阵元所组成的第个子系统;表示斜 距;
表示距离空不变相位;表示距离空变 相位;表示第个接收阵元和发射阵元所组成的第个子系统中收发阵元之间的距 离;表示参考距离;表示空变距离;表示线性调频信号中心频率;分别表示对 应于快、慢时间的瞬时频率;均表示中间 参数;分别表示水声声速与声呐平台拖曳速度。
作为本技术方案的进一步改进,所述S2中,针对每个接收阵元的数据在二维频域进行一致相位补偿,主要补偿与距离无关的相位误差,其公式如下:
其中,表示相位补偿函数;表示线性相位项,物理意义在于将 图像平移至参考距离处;表示发射信号的频谱,表示取共轭操作,表示虚数单 位。
作为本技术方案的进一步改进,所述S3中,针对每个接收阵元的数据在二维频域通过插值进行距离向坐标变换操作,即在二维频域采用Stolt插值的方式进行距离向坐标变换操作,将二维波数域内的弧形坐标变换成直线坐标,其变换后的距离向波数域坐标如下:
其中,表示在波数域进行距离向坐标变换后的坐标,表示 收发分置所造成的频率尺度因子。
作为本技术方案的进一步改进,所述S5中,针对融合后的数据在距离-多普勒域补偿方位走动误差,其补偿函数为:
其中,表示方位走动误差补偿函数。
本发明的目的之二在于,提供了一种方法运行计算平台装置,包括处理器、存储器以及存储在存储器中并在处理器上运行的计算机程序,处理器用于执行计算机程序时实现上述的新的多接收阵元合成孔径声呐波数域成像方法的部分计算步骤。
本发明的目的之三在于,提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现上述的新的多接收阵元合成孔径声呐波数域成像方法的部分计算步骤。
与现有技术相比,本发明的有益效果:
该新的多接收阵元合成孔径声呐波数域成像方法中,通过一致相位补偿操作对距离空不变相位进行补偿,降低了后续距离空变误差之间的耦合性,从而能提供更为精确的成像结果。
附图说明
图1为本发明中的流程图;
图2为实施例中的成像处理流程图;
图3为本发明中示例性的多接收阵元合成孔径声呐成像结果示意图;
图4为本发明中示例性的传统方法与本方案方法成像结果的方位剖面图;
图5为本发明中示例性的电子计算机平台装置结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
如图1-图4所示,本实施例提供了一种多接收阵元合成孔径声呐波数域成像方法,其实施流程图如图1所示,主要由以下步骤组成:首先计算一种较为精确的包含距离空变和距离空不变的二维频域系统函数,并针对每个接收阵元的数据在二维频域进行一致相位补偿,然后针对每个接收阵元的数据在二维频域进行距离向坐标变换操作,在二维频域对所有接收阵元数据融合后进行距离向傅里叶逆变换,在距离-多普勒域补偿方位走动误差,方位傅里叶变换后便得到最终的高分辨结果。
本实施例中,上述新的多接收阵元合成孔径声呐波数域成像方法具体包括如下步骤。
S1、首先计算一种较为精确的包含距离空变和距离空不变的二维频域系统函数;
本步骤中,原始的解析二维频域系统函数中含有距离变量,为使得二维频域 系统函数仅距离变量成正比,在计算包含距离空变和距离空不变的二维频域系统函数,针 对原始的二维频域系统函数,采用恒等式与近似式将原始解析 二维频域系统函数中的距离空变项、距离空不变项分离开来,其 计算公式如下:
其中,下标表示第个接收阵元和发射阵元所组成的第个子系统;表示斜 距;表示距离 空不变相位;表示距离空变相位;表 示第个接收阵元和发射阵元所组成的第个子系统中收发阵元之间的距离;表示参 考距离;表示空变距离;表示线性调频信号中心频率;分别表示对应于快、慢时 间的瞬时频率;均表示中间参数; 分别表示水声声速与声呐平台拖曳速度。
S2、针对每个接收阵元的数据在二维频域进行一致相位补偿;
本步骤中,针对每个接收阵元的数据在二维频域进行一致相位补偿,主要补偿与距离无关的相位误差,其公式如下:
其中,表示相位补偿函数;表示线性相位项,物理意义在于将 图像平移至参考距离处;表示发射信号的频谱,表示取共轭操作,表示虚数单 位。如图2所示,M表示接收阵元的总数,1≤m≤M。
S3、针对每个接收阵元的数据在二维频域进行距离向坐标变换操作;
本步骤即在二维频域采用Stolt插值的方式进行距离向坐标变换操作,将二维波数域内的弧形坐标变换成直线坐标,其变换后的距离向波数域坐标如下:
其中,表示在波数域进行距离向坐标变换后的坐标,表示 收发分置所造成的频率尺度因子。
S4、在二维频域对所有接收阵元数据融合后进行距离向傅里叶逆变换。
S5、针对融合后的数据,在距离-多普勒域补偿方位走动误差,方位傅里叶变换后得到最终的高分辨结果;
本步骤中,针对融合后的数据在距离-多普勒域补偿方位走动误差,其补偿函数为:
其中,表示方位走动误差补偿函数。
按照上述步骤对距离坐标为100m、方位坐标为7m理想点目标所产生的多接收阵元回波数据进行成像处理,成像结果如图3所示。
此外,本实施例还通过仿真实验,对本方案的方法、Callow在其发表论文“Wavenumber domain reconstruction of SAR/SAS imagery using single transmitterand multiple-receiver geometry”中提出的方法(以下简称为Callow的方法)以及黄攀在其发表论文“Synthetic aperture imagery for high-resolution imaging sonar”中提出的方法(以下简称为黄攀方法)对仿真数据进行成像后的方位剖面进行重叠对比,如图4所示。从图4的结果中不难发现,本方案的方法相对Callow的方法有较大的提升;而黄攀方法在5.7m处的旁瓣比本方案的方法要高约5.6dB,进一步体现了本方案的方法的优越性。
如图5所示,本实施例还提供了一种方法运行计算平台装置,该装置包括处理器、存储器以及存储在存储器中并在处理器上运行的计算机程序。
处理器包括一个或一个以上处理核心,处理器通过总线与存储器相连,存储器用于存储程序指令,处理器执行存储器中的程序指令时实现上述的新的多接收阵元合成孔径声呐波数域成像方法的部分计算步骤。
可选的,存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现,如静态随时存取存储器(SRAM),电可擦除可编程只读存储器(EEPROM),可擦除可编程只读存储器(EPROM),可编程只读存储器(PROM),只读存储器(ROM),磁存储器,快闪存储器,磁盘或光盘。
此外,本发明还提供一种计算机可读存储介质,计算机可读存储介质存储有计算机程序,计算机程序被处理器执行时实现上述的新的多接收阵元合成孔径声呐波数域成像方法的部分计算步骤。
可选的,本发明还提供了一种包含指令的计算机程序产品,当其在计算机上运行时,使得计算机执行上述各方面新的多接收阵元合成孔径声呐波数域成像方法的部分计算步骤。
本领域普通技术人员可以理解,实现上述实施例的全部或部分步骤的过程可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,程序可以存储于计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的仅为本发明的优选例,并不用来限制本发明,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (3)

1.一种多接收阵元合成孔径声呐波数域成像方法,其特征在于,包括如下步骤:
S1、首先计算包含距离空变和距离空不变的二维频域系统函数;
S2、针对每个接收阵元的数据在二维频域进行一致相位补偿;
S3、针对每个接收阵元的数据在二维频域进行距离向坐标变换操作;
S4、在二维频域对所有接收阵元数据融合后进行距离向傅里叶逆变换;
S5、针对融合后的数据,在距离-多普勒域补偿方位走动误差,方位傅里叶变换后得到最终的高分辨结果;
所述S1中,计算包含距离空变和距离空不变的二维频域系统函数,其计算公式如下:
其中,下标m表示第m个接收阵元和发射阵元所组成的第m个子系统;r表示斜距; 表示距离空不变相位;/>表示距离空变相位;dm表示第m个接收阵元和发射阵元所组成的第m个子系统中收发阵元之间的距离;rc表示参考距离;r'表示空变距离;fc表示线性调频信号中心频率;fτ和ft分别表示对应于快、慢时间的瞬时频率;/> 均表示中间参数;c和v分别表示水声声速与声呐平台拖曳速度;
所述S3中,针对每个接收阵元的数据在二维频域通过插值进行距离向坐标变换操作,即其变换后的坐标如下:
其中,kr表示在波数域进行距离向坐标变换后的坐标,表示收发分置所造成的频率尺度因子。
2.根据权利要求1所述的多接收阵元合成孔径声呐波数域成像方法,其特征在于,所述S2中,针对每个接收阵元的数据在二维频域进行一致相位补偿,其公式如下:
其中,Hm表示相位补偿函数;表示线性相位项,物理意义在于将图像平移至参考距离rc处;P(fτ)表示发射信号的频谱,*表示取共轭操作,j表示虚数单位。
3.根据权利要求1所述的多接收阵元合成孔径声呐波数域成像方法,其特征在于,所述S5中,针对融合后的数据在距离-多普勒域补偿方位走动误差,其补偿函数为:
其中,Hr(ft;r)表示方位走动误差补偿函数。
CN202310783214.XA 2023-06-29 2023-06-29 一种多接收阵元合成孔径声呐波数域成像方法 Active CN116520335B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310783214.XA CN116520335B (zh) 2023-06-29 2023-06-29 一种多接收阵元合成孔径声呐波数域成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310783214.XA CN116520335B (zh) 2023-06-29 2023-06-29 一种多接收阵元合成孔径声呐波数域成像方法

Publications (2)

Publication Number Publication Date
CN116520335A CN116520335A (zh) 2023-08-01
CN116520335B true CN116520335B (zh) 2023-09-05

Family

ID=87390570

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310783214.XA Active CN116520335B (zh) 2023-06-29 2023-06-29 一种多接收阵元合成孔径声呐波数域成像方法

Country Status (1)

Country Link
CN (1) CN116520335B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116953677A (zh) * 2023-09-18 2023-10-27 海底鹰深海科技股份有限公司 基于深度学习的声呐目标识别算法
CN116990791B (zh) * 2023-09-25 2023-11-28 海底鹰深海科技股份有限公司 一种多接收阵合成孔径声纳回波仿真方法
CN117538881B (zh) * 2024-01-10 2024-05-07 海底鹰深海科技股份有限公司 一种声呐水体成像波束形成方法、系统、设备及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772326A (zh) * 2016-12-05 2017-05-31 中国人民解放军91388部队 一种多子阵合成孔径声纳相位误差分析方法
CN106842210A (zh) * 2016-12-05 2017-06-13 中国人民解放军91388部队 一种新的多子阵合成孔径声纳快速成像算法
CN108037497A (zh) * 2018-01-04 2018-05-15 中国人民解放军91388部队 多子阵合成孔径声纳数据的收发合置转换方法
CN116106914A (zh) * 2023-01-09 2023-05-12 中国人民解放军海军工程大学 多接收阵sas六自由度运动误差补偿及成像方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10624612B2 (en) * 2014-06-05 2020-04-21 Chikayoshi Sumi Beamforming method, measurement and imaging instruments, and communication instruments
US11125866B2 (en) * 2015-06-04 2021-09-21 Chikayoshi Sumi Measurement and imaging instruments and beamforming method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772326A (zh) * 2016-12-05 2017-05-31 中国人民解放军91388部队 一种多子阵合成孔径声纳相位误差分析方法
CN106842210A (zh) * 2016-12-05 2017-06-13 中国人民解放军91388部队 一种新的多子阵合成孔径声纳快速成像算法
CN108037497A (zh) * 2018-01-04 2018-05-15 中国人民解放军91388部队 多子阵合成孔径声纳数据的收发合置转换方法
CN116106914A (zh) * 2023-01-09 2023-05-12 中国人民解放军海军工程大学 多接收阵sas六自由度运动误差补偿及成像方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
适用于宽带多子阵SAS成像的非线性CS算法;吕金华;《中国声学学会水声学分会2015年学术会议论文集》;全文 *

Also Published As

Publication number Publication date
CN116520335A (zh) 2023-08-01

Similar Documents

Publication Publication Date Title
CN116520335B (zh) 一种多接收阵元合成孔径声呐波数域成像方法
Chen et al. Motion compensation/autofocus in airborne synthetic aperture radar: A review
CN108051809B (zh) 基于Radon变换的运动目标成像方法、装置及电子设备
CN105759263B (zh) 一种高分辨率大场景下的星载斜视sar雷达成像方法
CN104898119B (zh) 一种基于相关函数的动目标参数估计方法
CN104749570B (zh) 一种移不变机载双基合成孔径雷达目标定位方法
CN105182340B (zh) 一种双基地前视sar运动补偿方法
CN106842210A (zh) 一种新的多子阵合成孔径声纳快速成像算法
Zhang et al. Focusing multireceiver SAS data based on the fourth-order legendre expansion
Zhou et al. A new fast factorized back projection algorithm for bistatic forward-looking SAR imaging based on orthogonal elliptical polar coordinate
CN110954899B (zh) 高海况下海面舰船目标成像方法及装置
WO2021250943A1 (en) Graph-based array signal denoising for perturbed synthetic aperture radar
CN112485783A (zh) 目标探测方法、装置、计算机设备和存储介质
CN113534077B (zh) 一种雷达辐射源威力反演方法、装置及电子设备
CN104931965B (zh) 一种固定发射站双基正侧视sar成像方法
CN116500626B (zh) 一种多接收阵元数据收发合置转换方法
CN114706076A (zh) 一种基于改进距离徙动算法的毫米波近场sar速成像方法
CN110596706A (zh) 一种基于三维图像域投射变换的雷达散射截面积外推方法
Zhang et al. High‐resolution imaging for the multireceiver SAS
CN105549010A (zh) 频域合成孔径雷达成像方法
CN117129994A (zh) 基于相位补偿核gnss-sar改进后向投影成像方法
CN116774178A (zh) 单频mimo弧形阵列雷达的近距频域三维成像方法
CN108802706A (zh) 基于位置标定的调频步进雷达信号目标抽取方法
Gong et al. High resolution 3d InISAR imaging of space targets based on PFA algorithm with single baseline
CN113640793A (zh) 基于mrf的实孔径扫描雷达超分辨成像方法

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
GR01 Patent grant
GR01 Patent grant