CN110058247B - 一种合成孔径声呐实时成像的方法 - Google Patents

一种合成孔径声呐实时成像的方法 Download PDF

Info

Publication number
CN110058247B
CN110058247B CN201910248973.XA CN201910248973A CN110058247B CN 110058247 B CN110058247 B CN 110058247B CN 201910248973 A CN201910248973 A CN 201910248973A CN 110058247 B CN110058247 B CN 110058247B
Authority
CN
China
Prior art keywords
sampling
synthetic aperture
echo data
imaging
receiving transducer
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
CN201910248973.XA
Other languages
English (en)
Other versions
CN110058247A (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.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi 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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN201910248973.XA priority Critical patent/CN110058247B/zh
Publication of CN110058247A publication Critical patent/CN110058247A/zh
Application granted granted Critical
Publication of CN110058247B publication Critical patent/CN110058247B/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
    • G01S15/8902Side-looking sonar
    • G01S15/8904Side-looking sonar using synthetic aperture techniques

Landscapes

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

Abstract

本发明公开了一种合成孔径声呐实时成像的方法,本发明不需要预先构建延时表,通过简化时延求解结构,采用分段寻址并参与虚拟孔径合成的处理方式,降低复杂度,减少缓存消耗,加快了运算速度,能够实时完成成像。本发明采用按批次进行孔径处理,数据采集和算法处理同步进行,成像过程不需要再存储全部回波数据,减少硬件开销,提高了算法运算效率。本发明所述技术具有可扩展性,即可根据实际需要,将延时索引计算模块化并行处理,进一步加快运算速度。

Description

一种合成孔径声呐实时成像的方法
技术领域
本发明涉及水声通信中合成孔径声呐技术领域,具体涉及一种合成孔径声呐实时成像的方法。
背景技术
合成孔径声呐(SAS)是一种高分辨率水下成像声呐,其原理是利用多个小尺寸基阵在方位向匀速运动虚拟成大尺寸孔径,从而获得与距离无关的高分辨率图像。SAS被广泛应用在水下环境勘察,沉船搜寻与海底管线测量等领域。随着SAS成像技术的发展和实际应用需求的提升,对成像系统实时性处理的要求日益增加。
常用的SAS波束形成的核心算法为延时求和成像算法,该算法根据球面波聚焦原理,计算合成孔径长度内各接收位置的时间延迟,通过延时叠加来对成像区域中的各像素点进行聚焦成像,从而得到整个目标场景的图像。算法的关键之处在于首先需要求解目标点与接收基阵之间的双程时延,然后根据基阵在运动的过程中,发射阵的水平开角不变,导致成像区域的大小与阵元相对位置关系不变的特点,将每个像素点的延时索引值预先计算好存放在延时表中,供孔径合成处理时直接对成像数据寻址叠加,从而完成延时求和。该算法需要较高的采样率,采样点数多,计算量很大。为提高处理速度,需由多个延时表同时完成寻址叠加,这需要有极大的存储空间。因此当前常用的SAS成像处理方法为:声呐平台将采集到的回波数据预先存储,再统一传送给算法处理机,然后由算法处理机中的后处理软件根据所设计的运行参数产生发射信号的参考信号,对回波数据进行距离向脉冲压缩,并采用基于延时表的成像算法完成孔径合成处理,最后由显控上位机接收成像数据并显示,无法满足实时性成像的要求。
发明内容
本发明的目的在于克服上述技术难以进行实时性处理的不足,提供一种低复杂度、运行效率高,时延短的SAS实时成像方法,在保证高分辨率的前提下能够有效缩短成像时间,可实现合成孔径声呐实时成像。
假设声纳平台的运动轨迹为直线,且空间采样均匀。设阵元尺寸为D,一个合成孔径长度内方位向像素点个数为LA,其距离向像素点个数为Nr,回波数据的总批次数为LINE,LINE的数值应大于LA。
具体实现步骤如下:
步骤(1):对接收基阵采集到的原始回波数据进行采样,采样频率为fs,采样后得到的信号为xk(m);其中,k=0,1,2,...,LINE-1表示采样后回波数据的批次数,m=0,1,2,...,Nr-1表示距离向采样点序号;
步骤(2):对当前第k批次的回波数据进行距离向脉冲压缩;根据SAS系统的运行参数产生发射信号的参考信号h(m),对xk(m)按下式处理得到yk(m);
yk(m)=F-1{F[xk(m)]×F[h*(m)]}
其中,F和F-1分别表示正逆傅里叶变换,h*(m)表示参考信号的复共轭;
步骤(3):计算成像区域内各像素点到接收换能器的相对距离rm
Figure BDA0002011845280000021
其中,Ts表示从测绘带边缘像素点到接收换能器的到达时间,c表示声波在水中传播的速度;
步骤(4):计算中间变量
Figure BDA0002011845280000022
步骤(5):将一个合成孔径长度内方位向像素点均匀分为Q段,按下式求解每段内的方位向像素点到接收换能器的初始相对距离的平方值
Figure BDA0002011845280000023
Figure BDA0002011845280000024
其中,q=0,1,2,..,Q-1表示段号;
步骤(6):由下式求解成像区域内各段像素点到接收换能器的距离索引
Figure BDA0002011845280000025
Figure BDA0002011845280000026
其中,每段内的方位向像素点到接收换能器的距离索引初始值
Figure BDA0002011845280000027
步骤(7):计算对应的延时索引
Figure BDA0002011845280000028
Figure BDA0002011845280000031
其中,int[]表示取整操作;
步骤(8):根据延时索引对yk(m)进行叠加,输出第q段的结果fk q(m);
Figure BDA0002011845280000032
同时判断当前q值,如果q<Q-1,返回步骤(5);否则,进入步骤(9);
步骤(9):对得到的Q次结果求和,输出最终值fk(m)
Figure BDA0002011845280000033
本发明的技术优点主要有:
1、本发明不需要预先构建延时表,通过简化时延求解结构,采用分段寻址并参与虚拟孔径合成的处理方式,降低复杂度,减少缓存消耗,加快了运算速度,能够实时完成成像。
2、本发明采用按批次进行孔径处理,数据采集和算法处理同步进行,成像过程不需要再存储全部回波数据,减少硬件开销,提高了算法运算效率。
3、本发明所述技术具有可扩展性,即可根据实际需要,将延时索引计算模块化并行处理,进一步加快运算速度。
具体实施方式
以下结合实施例对本发明作进一步说明。
假设声纳平台的运动轨迹为直线,且空间采样均匀。设阵元尺寸为D,一个合成孔径长度内方位向像素点个数为LA,其距离向像素点个数为Nr,回波数据的总批次数为LINE,LINE的数值应大于LA。
具体实现步骤如下:
步骤(1):对接收基阵采集到的原始回波数据进行采样,采样频率为fs,采样后得到的信号为xk(m);其中,k=0,1,2,...,LINE-1表示采样后回波数据的批次数,m=0,1,2,...,Nr-1表示距离向采样点序号;
步骤(2):对当前第k批次的回波数据进行距离向脉冲压缩;根据SAS系统的运行参数产生发射信号的参考信号h(m),对xk(m)按下式处理得到yk(m);
yk(m)=F-1{F[xk(m)]×F[h*(m)]}
其中,F和F-1分别表示正逆傅里叶变换,h*(m)表示参考信号的复共轭;
步骤(3):计算成像区域内各像素点到接收换能器的相对距离rm
Figure BDA0002011845280000041
其中,Ts表示从测绘带边缘像素点到接收换能器的到达时间,c表示声波在水中传播的速度;
步骤(4):计算中间变量
Figure BDA0002011845280000042
步骤(5):将一个合成孔径长度内方位向像素点均匀分为Q段,按下式求解每段内的方位向像素点到接收换能器的初始相对距离的平方值
Figure BDA0002011845280000043
Figure BDA0002011845280000044
其中,q=0,1,2,..,Q-1表示段号;
步骤(6):由下式求解成像区域内各段像素点到接收换能器的距离索引
Figure BDA0002011845280000045
Figure BDA0002011845280000046
其中,每段内的方位向像素点到接收换能器的距离索引初始值
Figure BDA0002011845280000047
步骤(7):计算对应的延时索引
Figure BDA0002011845280000048
Figure BDA0002011845280000049
其中,int[]表示取整操作;
步骤(8):根据延时索引对yk(m)进行叠加,输出第q段的结果
Figure BDA00020118452800000410
Figure BDA0002011845280000051
同时判断当前q值,如果q<Q-1,返回步骤(5);否则,进入步骤(9);
步骤(9):对得到的Q次结果求和,输出最终值fk(m)
Figure BDA0002011845280000052

Claims (1)

1.一种合成孔径声呐实时成像的方法,其特征在于,
假设声纳平台的运动轨迹为直线,且空间采样均匀;设阵元尺寸为D,一个合成孔径长度内方位向像素点个数为LA,其距离向像素点个数为Nr,回波数据的总批次数为LINE,LINE的数值应大于LA;
该方法包括以下步骤:
步骤(1):对接收基阵采集到的原始回波数据进行采样,采样频率为fs,采样后得到的信号为xk(m);其中,k=0,1,2,...,LINE-1表示采样后回波数据的批次数,m=0,1,2,...,Nr-1表示距离向采样点序号;
步骤(2):对当前第k批次的回波数据进行距离向脉冲压缩;根据SAS系统的运行参数产生发射信号的参考信号h(m),对xk(m)按下式处理得到yk(m);
yk(m)=F-1{F[xk(m)]×F[h*(m)]}
其中,F和F-1分别表示正逆傅里叶变换,h*(m)表示参考信号的复共轭;
步骤(3):计算成像区域内各像素点到接收换能器的相对距离rm
Figure FDA0002011845270000011
其中,Ts表示从测绘带边缘像素点到接收换能器的到达时间,c表示声波在水中传播的速度;
步骤(4):计算中间变量
Figure FDA0002011845270000012
步骤(5):将一个合成孔径长度内方位向像素点均匀分为Q段,按下式求解每段内的方位向像素点到接收换能器的初始相对距离的平方值
Figure FDA0002011845270000013
Figure FDA0002011845270000014
其中,q=0,1,2,..,Q-1表示段号;
步骤(6):由下式求解成像区域内各段像素点到接收换能器的距离索引
Figure FDA0002011845270000015
Figure FDA0002011845270000016
其中,每段内的方位向像素点到接收换能器的距离索引初始值
Figure FDA0002011845270000021
步骤(7):计算对应的延时索引
Figure FDA0002011845270000022
Figure FDA0002011845270000023
其中,int[]表示取整操作;
步骤(8):根据延时索引对yk(m)进行叠加,输出第q段的结果
Figure FDA0002011845270000026
Figure FDA0002011845270000024
同时判断当前q值,如果q<Q-1,返回步骤(5);否则,进入步骤(9);
步骤(9):对得到的Q次结果求和,输出最终值fk(m)
Figure FDA0002011845270000025
CN201910248973.XA 2019-03-29 2019-03-29 一种合成孔径声呐实时成像的方法 Active CN110058247B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910248973.XA CN110058247B (zh) 2019-03-29 2019-03-29 一种合成孔径声呐实时成像的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910248973.XA CN110058247B (zh) 2019-03-29 2019-03-29 一种合成孔径声呐实时成像的方法

Publications (2)

Publication Number Publication Date
CN110058247A CN110058247A (zh) 2019-07-26
CN110058247B true CN110058247B (zh) 2020-11-17

Family

ID=67318024

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910248973.XA Active CN110058247B (zh) 2019-03-29 2019-03-29 一种合成孔径声呐实时成像的方法

Country Status (1)

Country Link
CN (1) CN110058247B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109239771B (zh) * 2018-08-10 2020-01-31 杭州电子科技大学 一种基于非均匀背景介质的弹性波成像方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103425990A (zh) * 2012-05-15 2013-12-04 中国科学院声学研究所 一种基于纹理特征的合成孔径声纳目标检测方法
CN107367731A (zh) * 2016-05-11 2017-11-21 中国科学院声学研究所 适合非匀速直线航迹的sas成像和运动补偿方法

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005011047A2 (en) * 2003-04-09 2005-02-03 Mano Dorsey Judd Virtual antenna technology (vat) and applications
CN100529795C (zh) * 2006-01-13 2009-08-19 中国科学院声学研究所 用于合成孔径声纳信号实时处理的稳健运动速度跟踪方法
CN101551452B (zh) * 2008-04-01 2012-05-16 中国科学院声学研究所 一种合成孔径声纳联合运动补偿方法及系统
CN101644763B (zh) * 2008-12-12 2012-12-12 中国科学院声学研究所 基于集群处理机的合成孔径声纳实时信号处理方法及系统
CN103593877A (zh) * 2013-11-07 2014-02-19 清华大学 合成孔径声纳图像的仿真方法及系统
CN103645479B (zh) * 2013-12-20 2016-01-27 中国船舶重工集团公司第七一五研究所 一种合成孔径声纳快速实时成像方法
CN104240212B (zh) * 2014-09-03 2017-03-29 西安电子科技大学 基于目标特征的isar图像融合方法
CN104330795B (zh) * 2014-10-20 2017-01-11 北京理工大学 一种基于Keystone变换的地基合成孔径雷达快速成像方法
CN105429922B (zh) * 2015-11-10 2018-10-12 杭州电子科技大学 用于ddo-ofdm系统的基于梳状导频的信道估计法
US10107900B1 (en) * 2016-11-30 2018-10-23 United States Of America As Represented By The Secretary Of The Navy System for calculating displacement of a sonar array
CN107085208B (zh) * 2017-04-20 2022-12-13 中国人民解放军海军工程大学 基于分段曲面拟合的干涉合成孔径声纳复图像配准方法
CN106980117B (zh) * 2017-05-31 2019-08-06 西安电子科技大学 基于重叠孔径的毫米波雷达视频成像方法
CN107632304A (zh) * 2017-08-18 2018-01-26 杭州电子科技大学 一种大斜视条件下的合成孔径声呐成像处理方法
CN108490443B (zh) * 2018-02-24 2020-12-04 中国人民解放军海军工程大学 基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法
CN109239717A (zh) * 2018-08-15 2019-01-18 杭州电子科技大学 一种基于水下无人航行器的合成孔径声呐实时成像方法
CN109444901B (zh) * 2018-11-14 2021-02-26 杭州电子科技大学 一种异构环境下多子阵sas子孔径成像方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103425990A (zh) * 2012-05-15 2013-12-04 中国科学院声学研究所 一种基于纹理特征的合成孔径声纳目标检测方法
CN107367731A (zh) * 2016-05-11 2017-11-21 中国科学院声学研究所 适合非匀速直线航迹的sas成像和运动补偿方法

Also Published As

Publication number Publication date
CN110058247A (zh) 2019-07-26

Similar Documents

Publication Publication Date Title
CN109444901B (zh) 一种异构环境下多子阵sas子孔径成像方法
CN112505710B (zh) 一种多波束合成孔径声呐三维成像算法
JP2001511250A (ja) 多次元ビーム形成装置
CN109581388B (zh) 一种实时三维成像声纳的近场宽视角波束形成方法
CN107966677B (zh) 一种基于空间稀疏约束的圆阵模态域方位估计方法
CN104656073B (zh) 三维成像声纳波束形成方法及在多核处理器上的实现方法
CN110907938B (zh) 一种近场的快速下视合成孔径三维成像方法
CN109884627B (zh) 任意线阵构型的近程毫米波快速三维成像方法
CN109489796A (zh) 一种基于单元辐射法的水下复杂结构辐射噪声源定位识别与声辐射预报方法
CN108519602A (zh) 海底冷泉快速成像系统及其阵列信号处理方法
CN105259557A (zh) 一种多频率发射波束形成方法及应用
CN110058247B (zh) 一种合成孔径声呐实时成像的方法
CN109375227A (zh) 一种解卷积波束形成三维声成像方法
Gough et al. An efficient image reconstruction algorithm for a multiple hydrophone array synthetic aperture sonar
CN114185007B (zh) 基于多核dsp的空时联合降维自适应杂波抑制方法
CN108508446A (zh) 基于海底冷泉成像数据的扇形变换方法
CN111366904A (zh) 一种基于二维频谱相干合并的cfar检测方法
CN110412585B (zh) 一种基于mvdr的下视合成孔径三维成像方法及系统
CN108614234A (zh) 基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法
CN112558070A (zh) 圆周扫描地基sar的频域成像方法及装置
CN110554393B (zh) 基于深度学习的高对比度最小方差成像方法
CN108761433B (zh) 一种使用mimo声纳差合阵处理的高分辨成像方法
CN115390052A (zh) 一种成像声呐实物模拟装置及方法
CN111551943B (zh) 一种基于damas2的稀疏布阵高分辨三维声成像方法及系统
CN110632606B (zh) 一种基于空间重采样的二维快速傅里叶变换三维成像方法

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
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20190726

Assignee: ZHEJIANG REALLIN ELECTRON Co.,Ltd.

Assignor: HANGZHOU DIANZI University

Contract record no.: X2021330000745

Denomination of invention: A real-time imaging method for synthetic aperture sonar

Granted publication date: 20201117

License type: Common License

Record date: 20211111