1、第45卷第2 0 期2023年10 月舰船科学技术SHIP SCIENCEAND TECHNOLOGYVol.45,No.20Oct.,2023多子阵合成孔径声呐后向投影自聚焦成像方法姚永红，张旭，宋伟1(1.无锡职业技术学院控制技术学院，江苏无锡2 1412 1；2.中国科学院声学研究所东海研究站，上海2 0 18 15)摘要：针对多子阵合成孔径声呐的非理想运动导致后向投影（BackProjection,BP）成像结果散焦的问题，提出一种结合自聚焦的后向投影成像方法。该方法首先对各脉冲回波数据进行粗栅格BP成像，然后根据图像对比度最优化准则估计出后向投影子图的相位误差，同时计算得到运动误差。
2、最后，对相位补偿及包络搬移后的原始回波数据进行细栅格后向投影自聚焦成像，得到精聚焦图像。仿真结果表明，与传统的BP算法相比，该方法能有效改善运动误差下的多子阵SAS成像效果中显著提高点目标图像的分辨率、峰值旁瓣比及积分旁瓣比指标。关键词：合成孔径声呐；后向投影算法；自聚焦；多子阵；成像方法中图分类号：TB566文章编号：16 7 2-7 6 49(2 0 2 3)2 0-0 17 6-0 6A back projection autofocus imaging method for multi-receiver synthetic aperture sonar(1.Department of
3、Control Technology,Wuxi Institute of Technology,Wuxi 214121,China;2.Shanghai Acoustics Laboratory,Chinese Academy of Sciences,Shanghai 201815,China)Abstract:Aimed at the problem of multi-receiver synthetic aperture sonar(SAS)back projection(BP)image defocus-ing caused by the imperfect motion of the
4、platform,an improved BP imaging method combined with autofocus algorithm isproposed in this paper.Firstly,image with coarser grid for raw data of each receiver is reconstructed by BP algorithm.Secondly,the phase error of each tomographic image is estimated according to the contrast optimization crit
5、erion,and thecorresponding motion error is calculated.Finally,the range-compressed echo data can be corrected by phase compensationand range profile shifts,subsequent another BP combined with autofocus is applied to get a well-focused fine grid image,which is devoid of the previously observable degr
6、adation.The simulation results show that,compared with the conventionalBP algorithm,the proposed method can correct the image degradation,and significantly improve the resolution,peak side-lobe ratio and integral sidelobe ratio of the point target image.Key words:synthetic aperture sonar;back projec
7、tion;autofocus;multi-receiver;imaging method0引言合成孔径声呐（SyntheticAperture Sonar，SA S）是一种水下成像设备，其利用小孔径基阵的运动来合成等效的大孔径，实现方位向的高分辨率成像。与传统侧扫声呐相比，SAS在水下目标探测、识别上具有较大的优势2 。不同于合成孔径雷达，SAS通过多子阵方法来解决探测距离和方位向高分辨率之间的矛盾，进而提高测绘速率。另外，由于水声信号的传播较电磁信号需要更长的时间，使得“停-走-停”的近似在声呐成像中不再有效。多子阵结构及非“停-走-收稿日期：2 0 2 2-0 9-2 6作者简介：姚永红（
8、198 4），男，硕士，高级工程师，研究方向为声呐信号处理。文献标识码：AYAO Yong-hong,ZHANG Xu,SONG Weildoi:10.3404/j.issn.1672-7649.2023.20.033停”传播几何增加了成像算法的复杂性，需对传统的单基算法加以适应性改进，才能将其应用于SAS成像3-4。SAS后向投影算法（BackProjection,BP）是最基本的成像算法，其基本原理是计算各方位时刻平台位置与目标点的双程延时，再据此选取该目标点在各方位时刻对应的回波数据相干累加，进而得出目标函数。BP算法的图像重建过程不存在理论近似，且多子阵结构、非“停-走-停”传播几何及
9、空间采样不均匀等因素对其影响小，是一种精确的时域高分辨成像算法。计算量大的缺点限制了其在实时成像中的应用，但随着电子技术的发展和信号处理能力的不断提升，BP算第45卷法的研究又重新引起了人们的重视5-7 。平台沿理想直线运动是SAS获得高分辨图像的关键。但在实际工作中，海浪、暗涌和平台自身的不规则运动会影响航行轨迹，导致其产生运动误差。高分辨率SAS一般要求运动误差小于波长的1/8，否则会影响图像分辨率及聚焦效果。因此，需要对运动误差进行补偿才能获得聚焦良好的SAS图像8-9。由于运动传感器精度有限，不能满足高分辨率SAS成像要求，因此经过传感器测量运动信息补偿后的回波数据仍会有相当部分的残余
10、运动误差存在，引起接收信号相位误差，导致图像聚焦质量下降；且当残余运动误差引起的距离徙动大于一个距离分辨单元时，便会导致SAS图像二维散焦。解决上述问题的一种思路是用图像自聚焦方法估计出回波中的相位误差和残余距离徙动，并在回波数据中予以补偿后再成像，使成像质量得到一定程度的改善10-11。由于BP算法属于时域成像算法，因此传统适合频域成像算法的自聚焦方法难以与其结合。综合目前的研究成果，可知针对SAS时域成像的自聚焦问题还缺乏相关研究。本文通过分析非“停-走-停”条件下多子阵SAS后向投影算法成像运动误差模型，结合BP算法层析成像的特点及图像对比度最优化准则，提出一种多子阵SAS后向投影自聚焦
11、成像方法，有效解决了残余运动误差下BP成像散焦的问题，增大了BP算法的适用范围。1多子阵SAS后向投影成像模型1.1回波信号模型多子阵SAS数据采集几何模型如图1所示，场景中心点O定义为坐标原点，沿航迹方向定义为X方向。不失一般性，设声呐工作在斜视模式，斜视角为s，正侧视模式时，s=90。假设载体沿着理想航迹以速度匀速直线航行，u表示某方位时刻发射基元的方位向坐标，di为发射基元与第i个接收基元的间距。发射基元相位中心的瞬时坐标为(xa,Ja,za)，H 为航迹与场景的高度差，Rref为作用距离，即孔径中心与场景中心的距离。在成像场景所在坐标系中，发射基元相位中心的瞬时坐标在平台理想运动的情况
12、下可表示为：(xa,ya,za)=u,VRre in(0.)-H2,H对于场景中位置坐标为(x,y,0)的目标，在理想航迹情况下，声信号在发射接收期间的传播路程为：姚永红，等：多子阵合成孔径声呐后向投影自聚焦成像方法Va2+0a-y2+(xa-02+(xa+Ad-0)2-12实际情况下，受海浪、暗涌等因素的影响，声呐载体一般为非理想直线航行。假设其在、z向上分别存在Ay、z 大小的运动误差，令y=ya+Ay，Za=Za+Ay，那么发射基元的瞬时坐标为(xa,ya,za)、接收基元的瞬时坐标为(xa+A d;+v,y a,z a)。从而可得到运动误差情况下的声信号双程传播距离表达式为：R(xa,
13、Ja,x,y)=V(xa-x)2+(Va-y)?+(za)2+V(xa+Ad;+v*-x)2+(a-y)?+(a)?=Ro(xa,ya,x,y)+R(x a+d;+v t*)。(6)(1)式中：R为理想航迹与实际航迹下场景中某点目标P的声信号双程传播距离差。由于对于任意目标而言，R(xa,ya,x,y)与Ro(xa,ya,x,y)都是随接收基元位置:177合成孔径理想航迹vtRRiRT图1多子阵SAS成像几何模型Fig.1 Imaging geometry of the multiple receiver SASRo(Xa,ya,X,y)=RT+RRio式中：Rr为发射换能器到目标的传播距离；
14、RRi为非“停-走-停”情况下第i个接收基元到目标的距离，用坐标表示为：RT=/xa-x)2+(Va-y)?+za?,(RRi=V(xa+Ad;+v-x)?+(Va-y)?+za?。式中，*为声信号从发射到被第i个接收基元接收所经历的时间。声信号发射接收期间的双程传播路程等于目标与发射阵、接收阵的实际距离之和，即ct*=Ro(xa,ya,X,y),式中，C为声速。解式（4）得到：2/2HRrefP(x,y,0)4zX(2)(3)(4)(5)178Xa+d;+v t*变化，所以R也是xa+Ad;+vt*的函数。声呐发射线性调频信号为：s(t)=rect(/T)exp(jnKt2)exp(j2元f
15、et)。其中：t为距离向时间；K为调频斜率；T为信号的脉宽；rect(*)为信号包络。rect对回波信号接收解调及距离脉压后可得s(Xa,T)=0(Xa,x,y)x(x.y)ORo(xa,Ja,x,y)+ARsincTKt-CRo(Xa,Ja,X,y)+AR)exp-j2元.feC式中：(xa,x,y)为目标对于发射基元在xa的声反射系数；Q为成像场景范围。1.2BP成像及运动误差的影响对式（9）所示声呐回波数据进行后向投影处理来重Ro(xa,ya,x,y)建场景图像，重建过程中需乘以expj2元feC用于后续层析图像的相干叠加。各方位脉冲得到的层析图像为：Sm(x,y)=s(xa,T).ex
16、p/j2元fe0(xa,x,y)x(x.y)eRo(xa,a,x,y)+ARsinKRexp-j2元/式中：m为方位向回波脉冲索引，xa与ya都随m而变化；Sm(x,y)为第m个回波脉冲后向投影所得复图像，为了分析方便将R写成R(m)的形式。将合成孔径时间内各方位回波后向投影得到的层析图像相干叠加，得到BP成像结果：Na-1Sm(x,y)=m=0m=0(x.y)eQRo(xa,ya,x,y)+R(m)sincTKexp12元舰船科学技术可知，运动误差R(m)使目标回波在距离向存在偏移，同时引入相位误差。在运动传感器精度不够的情况下，运动误差是未知的，所以一方面相位误差会导(7)致层析图像相干叠
17、加失败，另一方面超出一个距离分辨单元的包络偏移会使重建后的目标位置在距离向出现偏移，两者导致图像的二维散焦。若能通过图像自聚焦技术估计出相位误差，并利用其计算得到相应的运动误差R(m)，则能够改善图像聚焦质量。(8)lo,others。dxdyRo(xa,ya,x,y)CCdxdy。Na-1o(xa,x,y)xCAR(m)dxdy。C第45卷2结合自聚焦的BP成像根据式(11)，当R(m)小于一个距离门宽度时，只需对各脉冲后向投影子图进行相位误差补偿，然后通过相干叠加便能得到聚焦良好的SAS图像：Na-1I(x,y)=2sm(a.y).exp(-jic(m).(9)m=0式中：(m)为第m个脉
18、冲对应的相位误差补偿量，当R(m)0(m)=-2元fc2时，相位误差被完整补偿，图像可C完全聚焦。SAS图像的对比度为图像强度的标准差与均值的比值：1(x,)P-E(x,)P)E(I(x,)式中：E()为求均值；I(x,y)P为图像中坐标为(x,J)点像素的强度。对于给定的SAS图像，信号能量E(I(x,y)是一个常数，因此式(13)可进一步等价为：十(14)X=0y=0(10)式中：N和N,分别为图像在方位向与距离向的像素单元个数。基于图像对比度的自聚焦算法是以图像对比度最优化为标准，对比度越高则聚焦效果越好；而当SAS图像中相位误差为0 时，其对比度则最大。图像自聚焦过程也是相位误差补偿的
19、过程，相位误差补偿精度越高，图像的聚焦效果越好。对比度是(m)的函数，通过求解对比度的最大值问题，即式(14)对6(m)求导，得到对比度最大时的极值点即为相位误差函数的X估计值为：(11)(12)(13)Nx-1Ny-1(15)第45卷式中：A(m)=x=0y=0分析多子阵SAS信号双程传播距离表达式发现，即便声呐载体在y向、z向上的运动误差均小于一个距离单元，R(m)在某些方位位置会大于一个距离单元，若不将该运动误差在BP成像之前的距离脉压信号上予以补偿，则会导致自聚焦失败。针对该问题，先通过粗栅格BP成像以降低距离分辨率，使得R(m)小于一个距离分辨单元，从而可利用对比度最优准则估计出相位
20、误差，同时根据相位误差估计结果计算出相应运动误差R(m)，然后对脉压后的原始回波数据进行相位及包络补偿。最后，对补偿后的回波数据进行细栅格BP成像得到精聚焦的SAS图像。具体实现步骤如下：步骤1设置BP成像栅格距离向宽度为设计值的一倍，将各方位脉压后的回波数据分别进行后向投影得到粗栅格复图像Sm(x,y)；同时初始化迭代次数索引=0，初始化各脉冲相位误差为0，初始化对比度门限为To。步骤2 通过各脉冲后向投影复图像Sm(x,y)及相位误差(m)，求出Ii(x,y),并根据式(14)计算对比度1Cr；按式(15)估计出相位误差函数exp(-j(m)，将各方位后向投影子图Sm(x,y)乘以exp(
21、-j(m)叠加进行更新，并叠加得到新的Ii+1(x,y)。步骤3按式（14）计算图像对比度C+1，同时计算对比度增加的比例T=(Ci+1-C)/CI，将T,与初始化设定的迭代终止门限T。比较，当TiTo时则返回步骤2继续进行下一次迭代，否则送代终止。步骤4对每次迭代所估计出的相位误差求和得到总的相位误差(m)，将N。个方位脉冲对应的相位误差解缠绕处理，并根据Q(m)=-2元fcR(m)/c 计算出各方位脉冲包络偏移量的估计值R(m)。步骤5通过插值对脉压后的各方位原始回波数据进行大小为R(m)包络搬移，同时乘以相位误差补偿量exp(-jd(m)。然后设置BP成像栅格距离向宽度为系统设计值，进行
22、细栅格BP成像，得到各方位脉冲后向投影得到高分辨率复图像Sm(x,y)。步骤6 重复步骤2 与步骤3，当T小于终止门限T。时，步骤2 中的I(x,y)即为最终的聚焦好的BP成像结果。该步骤是对细栅格（距离像素为系统设计值）姚永红，等：多子阵合成孔径声呐后向投影自聚焦成像方法Nx-1Ny-1Sm(x,y)-1(x,y)P1*(x,y)。1+1179.下的BP图像进行自聚焦处理，由于经过前面步骤1步骤5的处理，引起图像距离向散焦的残余距离运动(16)误差得以补偿，故只需对细栅格BP图像进行方位向的相位补偿处理，即步骤2 和步骤3的处理。上述步骤中包含2 次BP成像过程，其中第1次是粗距离分辨率成像
23、，该次成像是为了避免跨距离单元的运动误差对相位估计的影响；第2 次是对回波数据进行运动误差、相位误差补偿后的高分辨率成像。若运动误差R(m)在整个合成孔径时间内都小于一个距离单元，则可以直接进行高分辨率的BP自聚焦成像，以节省成像时间。本文提出的多子阵合成孔径声呐后向投影自聚焦成像方法信号处理流程如图2 所示，其中Pyo为系统设计距离向像素大小。初始化BP成像栅格距离向宽度Py-2Pyo初始化选代次数1-0 初始化对比度门限T初始化相位误差。(0),(1).0(Na-1)为0输入各接收回波距离脉压数据Py-Pyo?Y包络搬移、相位误差补偿BP得到图像域各脉冲对应子图Sm(x,J)各子图乘以估计
24、相位误差函数并叠加：计算对比度:C-23/0,0NY计算对比度增加的比例：TF(CM-C)/C,T,T.?YPy-Pyo?Y1,y)为BP自聚焦最终成像结果图2 多子阵SAS后向投影自聚焦成像信号处理流程图Fig.2 Signal processing flowchart of BP autofocus formultiplereceiverSAS3仿真实验及分析3.1仿真参数仿真参数设置见表1，点目标回波时延的由式（5）给出。仿真场景为理想平地，区域内共设置9个点目标，1个位于场景的坐标原点，另外8 个均匀分布在更新Py-Pyo1-1+14(m)exp(j0(m)=A(m)N各脉冲相位误差函
25、数的估计：Nexp(j0(m)-hexp(-j0(m)并计算得到包络偏移量AR(m):180:以原点为中心、边长为10 m的正方形4条边上。为验证新方法对声呐残余运动误差的补偿效果，在一个合成孔径时间内加入1个周期的侧摆扰动和2 个周期的升沉扰动，扰动形式都为正弦震荡，幅度分别为0.0 2 m、0.01m，运动误差位移曲线如图3所示。表1系统仿真参数Tab.1 Parameters for system simulation参数中心频率/kHz信号带宽/kHz脉冲重复周期/s发射阵宽度/m接收阵宽度/m接收阵个数斜视角/)作用距离/m场景与平台的相对高/m平台速度/ms-0.020.010-0
26、.01-0.02012345678方位时间/s图3运动误差曲线Fig.3Motion errorcurves3.2结果及分析对点目标仿真数据进行成像处理，传统的后向散射算法成像结果如图4（a）所示。可以看出，9个点目标都呈“长条状”，存在严重的散焦。本文所提方法的成像结果如图4（b）所示，其中各点目标聚焦良好。为了衡量本文所述方法的成像效果，取右下角位置点目标周围的6 46 4像素进行放大显示，各步骤对应的结果如图5所示。图5（a）是没有自聚焦的传统BP成像结果，点目标响应在方位向分散于31个像素单元里，在距离向分散于3个像素单元里，总体呈现两维散焦，且方位散焦尤为严重。图5（b）是利用本文方
27、法进行处理，但未对包络偏移进行补偿的结果，点目标响应在方位向分散于4个像素里、距离向上分散于3个像素里，可见经过相位误差补偿后方位向的舰船科学技术聚焦情况改善明显，但距离向散焦情况没有改善。图5（c）是经过本文方法全部步骤处理得到的成像结果，点目标响应基本聚焦在1个像素单元里，可见本文方法能够很好地处理BP成像中残余运动误差的影响，明显改善成像质量。数值110200.20.060.033270100202.4一侧摆运动误差-升沉运动误差第45卷(a)传统的BP成像结果(b)结合自聚焦的BP成像结果图4运动误差下的成像结果Fig.4 Imaging results of simulation d
28、ata with motion error(a)传统BP(c)包络、相位补偿图5场景中某边缘点目标的成像结果对比Fig.5 Comparison of imaging results of an edge point target withdifferent method表2 给出了图3中3幅图像的熵和对比度。由表中数据可直观的看出，相比传统的BP算法，结合自聚焦的BP算法所得图像的熵降低了50.2%、对比度提高了36 9%，进一步证明了本文方法能够显著提高SAS图像的聚焦效果。表2 图像熵及对比度的比较Tab.2 Image entropy and contrast of different
29、 methods本文方法成像方法传统BP仅相位补偿包络、相位补偿图像5.167对比度6.631选取某边缘点目标的两维剖面进行量化分析，为了显示效果，对数据作16 倍内插处理，点目标的方位、距离剖面如图6 所示。其中，虚线为未作自聚焦的传统BP成像结果，实线为本文方法成像结果。由图6 可知，残余运动误差导致传统BP成像结果中点目标方位向主瓣严重展宽、旁瓣抬高，且距离向旁瓣升高；而通过本文方法处理后，两维聚焦效果改善明显。为了能更清楚地比较2 种方法的成像效果，对(b)相位补偿3.5062.57317.10631.133本文方法指标增幅/%-50.2360.9第45卷-10-20-304000-1
30、020-30-400图6 点目标的两维剖面Fig.6Range and azimuth slices of point target图6 中的目标分辨率（IRW）、峰值旁瓣比（PSLR）、积分旁瓣比（ISLR）进行测量，结果如表3所示。可知，经本文方法处理的点目标积分旁瓣比、峰值旁瓣比明显降低，方位向分辨率提高了4.5倍，进一步证明了本文方法能够很好地消除残余运动误差给BP成像带来的影响，提高点目标成像质量。表3两种方法的点目标质量指标对比Tab.3 Point target quality parameters of the two methods成像方法本文方法ISLR/dB-0.20距离
31、向PSLR/dBIRW/cmISLR/dB方位向PSLR/dBIRW/cm4结语本文提出一种用于高分辨率多子阵SAS的后向投影自聚焦成像方法，该方法在BP算法的基础上，以图姚永红，等：多子阵合成孔径声呐后向投影自聚焦成像方法0.传统BP-BP自聚焦，100200方位向（0.0 0 2 34m/采样间隔）(a)点目标方位向剖面.传统BP-BP自聚焦100200距离向（0.0 0 2 34m/采样间隔）(b）点目标距离向剖面传统BP-1.293.754.99-0.8416.40181像对比度大小为准则、通过多次迭代估计相位误差。同时针对超出一个距离分辨单元的距离动误差导致自聚焦效果不理想的问题，提
32、出首先通过在粗分辨率图像下利用自聚焦估计出相位误差和残余运动误差，然后利用估计结果对原始回波数据补偿，最后再进行300400500300400细栅格下的后向投影自聚焦成像以得到聚焦良好的SAS图像。仿真试验表明，新方法能够很好地在BP算法中对残余运动误差及相位误差进行补偿，显著改善目标响应特性，提升成像质量，证明了该方法的有效性和可行性。500参考文献：1】刘纪元.合成孔径声呐技术研究进展1.中国科学院院刊,2019,34(3):283-288.【2 刘昊，刘维，刘纪元.合成孔径成像中掩埋目标的深度误差分析.声学技术,2 0 15,34(3):2 14-2 18.【3徐江，唐劲松，张春华，等.
33、多子阵合成孔径声纳逐点成像算法1.信号处理,2 0 0 3,19(2):157-16 0.【4】吕金华，唐扶光，赵煦，等.一种小斜视多接收阵合成孔径声呐距离多普勒成像算法J.舰船科学技术，2 0 19,41(7):124-129.5 ZHANG X B,DAI X T,YANG B.Fast imaging algorithm forthe multiple receiver synthetic aperture sonarsJ.IET Radar,Sonar&Navigation,2018,12(11):1276-1284.6 ZHANG X B,YANG B,et al.BP algori
34、thm for the multi-re-ceiver SASJ.IET Radar,Sonar&Navigation,2019,13(5):-10.89830-838.-15.23【7】许炎义，钟何平，唐劲松.多接收阵合成孔径声呐后向投影成像快速算法.武汉大学学报，2 0 15,40(10)：140 9-1413.3.75【8 张羽，王朋，刘纪元，等.合成孔径声呐重叠相位中心与惯性9.73导航系统联合估计运动误差算法.兵工学报，2 0 2 1,42(3):-10.73588-597.3.629】马梦博，唐劲松，田振，等.基于惯导系统的多子阵SAS运动补偿.华中科技大学学报，2 0 2 0,48(8)：7 3-7 8.10】陈家瑞，张劲东，班阳阳，等.基于对比度最优准则的反投影自聚焦方法.雷达科学与技术,2 0 15,13(1):2 1-2 6.11】段江涛，黄勇，刘纪元.PGA算法在条带式SAS场景目标成像中的应用.应用声学，2 0 17,36(2)：17 5-18 1.
