一种合成孔径声纳成像的自聚焦处理方法技术

本发明专利技术涉及一种合成孔径声纳成像的自聚焦处理方法，主要包括以下步骤：（１）每次回波信号经过复基带解调以及匹配滤波处理得到宽带解析信号ｒ↓［ｐ］（ｔ）；（２）相邻回波解析信号之间相关处理；（３）通过求复数相角，得到相邻回波解析信号之间的相位误差估计；（４）通过空间位置累积处理获得运动误差引起的累积相位误差；（５）利用获得的相位误差序列｛＊↓［ｐ］｝进行补偿来实现回波数据校正；（６）补偿后数据应用相应的图像重构算法，则可获得聚焦的合成孔径声纳图像。本发明专利技术的有益的效果是：它不需要迭代运算，并具有较好的实时运算性能和较好的图像宏观校正能力。并且进一步通过α加权削弱强目标散射点和噪声对运动误差估计的影响，从而拓宽了应用条件。

【技术实现步骤摘要】

本专利技术涉及合成孔径声纳成像的信号处理，主要是，用于获得更加精细的海底地貌测绘和沉物探测的声成像。
技术介绍
合成孔径声纳(Synthetic Aperture Sonar，简称SAS)是一种高分辨力声学成像设备。它是利用小孔径基阵匀速直线运动，在运动过程中相继发射和接收信号，并将一段运动历程中接收信号进行相干迭加，从而得到等效于实际物理声阵几倍到几十倍的合成虚拟孔径阵，获得很高的分辨力和空间增益。然而在实际航行过程中这个条件并不满足，使得合成孔径图像产生强度失真和几何失真。实现运动误差补偿的自聚焦处理方法是合成孔径声纳的一项关键技术。合成孔径声纳的估计和校正方法通常分为两类一类为采用传感器的硬件补偿，另一类为自聚焦算法补偿。由于传感器精度难于满足合成孔径声纳的要求，合成孔径声纳一般都需要采用自聚焦算法来实现运动误差补偿。自聚焦算法通常又可分为数据域自聚焦算法(如冗余相位中心自聚焦算法)和图像域自聚焦算法(如相位梯度自聚焦)算法两种。国外合成孔径声纳所采用的自聚焦算法主要是冗余相位中心(Redundance Phase Center，简称RPC)自聚焦算法，它是利用接收回波信号的冗余性来处理的，需要采用冗余子阵设计或牺牲空间采样来实现。经典的相位梯度自聚焦(Phase Gradient Autofocus，简称PGA)算法，需要进行迭代处理才具有良好的估计性能，同时经典自聚焦处理方法需要选择强目标散射点。
技术实现思路
本专利技术的目的是为了克服上述方法的不足，而提供。本专利技术解决其技术问题所采用的技术方案。这种合成孔径声纳成像的自聚焦处理方法，主要包括以下步骤 (1)、每次回波信号经过复基带解调以及匹配滤波处理得到宽带解析信号rp(t)；(2)、相邻回波解析信号之间相关处理；相邻两个发射脉冲回波的宽带解析信号rp(t)、rp-1(t)的互相关运算为rrp=∫0Trp*(t)rp-1(t)dt----(5)]]>式中*表示复共轭；(3)、通过求复数相角，得到相邻回波解析信号之间的相位误差估计Δφ^p=Arg{rrp}----(6)]]>(4)、通过空间位置累积处理获得运动误差引起的累积相位误差ϵ^p=ϵ^p-1+Δφ^p;ϵ^0≡0----(7)]]>(5)、利用获得的相位误差序列 进行补偿来实现回波数据校正RpIdeal(f)=Rp(f)exp(jϵ^p)----(8)]]>式中，Rp(f)为rp(t)的傅立叶变换；(6)、补偿后数据应用相应的图像重构算法，则可获得聚焦的合成孔径声纳图像。本专利技术的有益的效果是它是一种利用合成孔径声纳混响中冗余信息的数据域相关处理的运动误差估计方法。采用以均自聚焦算法为核心的、通过测量宽带解析信号复相关的相位角进行合成孔径声纳运动误差估计和补偿，从而获得合成孔径声纳成像的自聚焦处理。它不需要迭代运算，并具有较好的实时运算性能和较好的图像宏观校正能力。并且进一步通过α加权削弱强目标散射点和噪声对运动误差估计的影响，使得切变平均自聚焦算法既适应于空间缓变场景的宽波束SAS处理；也适用于存在强目标散射点的场景，从而拓宽了应用条件。附图说明图1本专利技术实施例2的α加权切变平均自聚焦处理方法的流程图；具体实施方式下面结合附图和实施例对本专利技术作进一步介绍实施例1一、基于混响估计的运动误差模型对于一个混响限制(混响占主导地位，而噪声相对较弱)慢变化场景海底模型，假定海底的各个散射体的尺寸远小于波长，可以把它们看作“点源”目标，在这种近似下可以不考虑频率和波束方向对散射的影响，这样混响模型可看为很多点源目标回波的线性之和。如果在测绘区域内共有Q个散射体，则第p个发射脉冲的回波可描述为ep(t)=Σq=1Qsp(t)⊗h(t,τpq)+np(t)----(1)]]>式中表示卷积操作，sp(t)是宽带发射信号，h(t，τpq)为点扩展函数(包括介质中传播函数，散射体散射响应以及基阵响应函数)，np(t)为噪声项。在没有运动误差时，回波信号经过复基带解调以及匹配滤波处理得到解析信号rp(t)，表示为rpIdeal(t)=|rpIdeal(t)|exp(jψpIdeal(t))----(2)]]>由于混响场中散射体服从均匀分布，接收回波为各个散射体散射信号之和，故而ψpIdeal(t)为内均匀分布的随机变量。假设在声学换能器基阵发生摆动误差xp，而散射体的主体近似不变，则对于场景中每一个散射体的回波延时误差都等于εp≈2xp/c，第p个回波可表示为rp(t)=rpIdeal(t-ϵp)exp(-j2πf0ϵp)]]>=|rpIdeal(t-ϵp)|exp(jψpIdeal(t-ϵp))exp(-j2πf0ϵp)----(3)]]>故而获得测量相位为ψp(t)=Arg{rp(t)}=mod2π----(4)]]>式中Arg{·}为取复数相位操作。如果我们能够对每一个ψp(t)进行去卷绕处理，则获得一个常量相移2πf0εp和一个随机变量ψpIdeal(t-εp)，进一步通过在时间t上做平均处理就可获得常量相移，从而求出误差时延估计 二、合成孔径声纳成像的切变平均自聚焦处理方法(1)每次回波信号经过复基带解调以及匹配滤波处理得到宽带解析信号rp(t)。(2)相邻回波解析信号之间相关处理。利用(4)式中单个脉冲回波在时间t上做平均处理获得的常量误差相位2πf0εp一般大于2π从而会导致相位模糊问题。一个可行的解决这种相位卷绕的方法就是估计相邻两个发射回波信号之间的相位差，其核心是一种相关处理。相邻两个发射脉冲回波的宽带解析信号rp(t)、rp-1(t)的互相关运算为rrp=∫0Trp*(t)rp-1(t)dt----(5)]]>式中*表示复共轭。(3)通过求复数相角，得到相邻回波解析信号之间的相位误差估计Δφ^p=Arg{rrp}----(6)]]>(4)通过空间位置累积处理获得运动误差引起的累积相位误差ϵ^p=ϵ^p-1+Δφ^p;ϵ^0≡0----(7)]]>(5)利用获得的相位误差序列 进行补偿就可实现回波数据校正RpIdeal(f)=Rp(f)exp(jϵ^p)----(8)]]>式中，Rp(f)为rp(t)的傅立叶变换。(6)补偿后数据应用相应的图像重构算法，则可获得聚焦的合成孔径声纳图像。实施例2α加权切变平均自聚焦处理方法为了提高误差相位变化量估计的精度，K.本文档来自技高网...

【技术保护点】
一种合成孔径声纳成像的自聚焦处理方法，其特征在于：主要包括以下步骤：（１）、每次回波信号经过复基带解调以及匹配滤波处理得到宽带解析信号ｒ↓［ｐ］（ｔ）；（２）、相邻回波解析信号之间相关处理；相邻两个发射脉冲回波的宽带解析信号 ｒ↓［ｐ］（ｔ）、ｒ↓［ｐ－１］（ｔ）的互相关运算为：ｒｒ↓［ｐ］＝∫↓［０］↑［Ｔ］ｒ↓［ｐ］↑［＊］（ｔ）ｒ↓［ｐ－１］（ｔ）ｄｔ（５）式中＊表示复共轭；（３）、通过求复数相角，得到相邻回波解析信号之间的相位 误差估计：Δ＊↓［ｐ］＝Ａｒｇ｛ｒｒ↓［ｐ］｝（６）（４）、通过空间位置累积处理获得运动误差引起的累积相位误差：＊↓［ｐ］＝＊↓［ｐ－１］＋Δ＊↓［ｐ］；＊↓［０］≡０（７）（５）、利用获得的相位误差 序列｛＊↓［ｐ］｝进行补偿来实现回波数据校正：Ｒ↓［ｐ］↑［ｌｄｅａｌ］（ｆ）＝Ｒ↓［ｐ］（ｆ）ｅｘｐ（ｊ＊↓［ｐ］）（８）式中，Ｒ↓［ｐ］（ｆ）为ｒ↓［ｐ］（ｔ）的傅立叶变换；（６）、补偿后数据应用相应的图像重 构算法，则可获得聚焦的合成孔径声纳图像。...

【技术特征摘要】

【专利技术属性】
技术研发人员：丛卫华，张峰山，谭咏青，何卫彬，董彦涛，朱必波，章灵祥，何永军，曹海林，施国全，付菊英，王欣，杨晓帆，
申请(专利权)人：中国船舶重工集团公司第七一五研究所，
类型：发明
国别省市：86[中国|杭州]