专利摘要
专利摘要
本发明基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法，获取N个子阵的回波，对单个子阵的回波数据进行二维快速傅立叶变换，将回波数据变换到二维波数域，并求解单个子阵回波的二维波数谱解析解；对各子阵的方位向波数谱进行延展；矫正各个子阵的方位向时间偏移；消除数据基线长度对二维波数谱的影响，同时完成参考函数相乘和距离向脉冲压缩；基于非均匀傅里叶变换进行Stolt映射及距离向IFFT；将所得结果数据进行相干叠加；对叠加后的距离多普勒域中的数据进行方位向的IFFT实现成像。本发明构建了多接收子阵SAS在非停走停假设下几何模型，适用于多接收子阵SAS成像的精确的二维波数谱解析解，通过仿真实验和湖试数据成像结果验证了其正确性和有效性。
权利要求
1.基于解析解及NUFFT的多子阵合成孔径声纳ωk成像方法，其特征在于：
步骤1，获取N个子阵的回波；
步骤2，对单个子阵的回波数据进行二维快速傅立叶变换，将回波数据变换到二维波数域，并求解单个子阵回波的二维波数谱解析解；
步骤3，对各子阵的方位向波数谱进行延展；
步骤4，矫正各个子阵的方位向时间偏移；
步骤5，消除数据基线长度对二维波数谱的影响，同时完成参考函数相乘和距离向脉冲压缩；
步骤6，基于非均匀傅里叶变换进行Stolt映射；
步骤7，将步骤6所得结果数据进行相干叠加；
步骤8，对叠加后的距离多普勒域中的数据进行方位向的IFFT实现成像；
所述步骤1获取的第i个子阵的回波为
s
其中，i＝1,2,3……N，A
2.根据权利要求1所述的基于解析解及NUFFT的多子阵合成孔径声纳ωk成像方法，其特征在于，所述步骤2中第i个子阵的回波的二维波数谱解析解为：
其中，
其中h
3.根据权利要求2所述的基于解析解及NUFFT的多子阵合成孔径声纳ωk成像方法，其特征在于，所述步骤3，对方位向波数谱进行复制延展得到：
其中，K
4.根据权利要求3所述的基于解析解及NUFFT的多子阵合成孔径声纳ωk成像方法，其特征在于，所述步骤4，矫正各个子阵的方位向时间轴的方法包括：
将步骤4所得结果乘以
5.根据权利要求4所述的基于解析解及NUFFT的多子阵合成孔径声纳ωk成像方法，其特征在于，所述步骤5，消除数据基线长度对二维波数谱的影响，同时完成参考函数相乘和距离向脉冲压缩的方法包括：
将步骤4所得结果乘以
6.根据权利要求5所述的基于解析解及NUFFT的多子阵合成孔径声纳ωk成像方法，其特征在于，所述步骤6对步骤5所得结果进行基于NUFFT的Stolt映射，得到经过该映射后回波变化到距离多普勒域：
7.根据权利要求6所述的基于解析解及NUFFT的多子阵合成孔径声纳ωk成像方法，其特征在于，所述步骤7对步骤6所得的各个子阵的结果进行相干叠加。
8.根据权利要求7所述的基于解析解及NUFFT的多子阵合成孔径声纳ωk成像方法，其特征在于，对所述步骤8对步骤7所得的结果进行方位向逆傅立叶变换，变换后的结果即为成像结果。
说明书
技术领域
本发明属于成像算法技术领域，具体涉及一种基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法。
背景技术
一般来说，合成孔径声纳(SAS，Synthetic Aperture Sonar)采用多个接收子阵以提高测绘速率，通过增加发射阵元的数量来提高探测距离。多接收子阵SAS在提高测绘速率时，由于平台运动速度增大，而水声的传播速度较低，使得采用偏置相位中心(DPC，Displaced Phase Center)近似将收发分置模式下的回波数据转化为收发合置模式下的回波进行处理时，需要在非停走停假设下予以重新考虑。针对这一问题，对非停走停假设下的DPC近似乘以一个相位项来补偿DPC近似引入的误差，但在方位高分辨率的要求下该方法仍不能满足需要。可以在非停走停模式下通过等效的基线长度对DPC模型进行修正，但该方法对近距离的目标成像时有一定的散焦问题。近年来双基地合成孔径雷达(Bistatic SAR)回波的二维频谱计算取得了较大的进展，采用瞬时多普勒波数(IDW，Instantaneous Dopplerwavenumber)和半双基地角(Half Quasi-Bistatic Angle，HQBA)得到了回波的二维波数谱，在顺轨Bistatic SAR几何模型下得到了二维波数谱中HQBA的精确解析解，并分别采用ωk算法(ωkA，ωk Algorithm)和距离多普勒算法(RDA，Range Doppler Algorithm)对仿真数据进行了成像。经比较，该方法得到的解析解在一定条件下优于LBF(Loffeld’sBistatic Formula)方法、MSR(Method of Series Reversion)方法以及DMO(Dip MoveOut)方法得到的解。以上二维波数谱解析解都是对Bistatic SAR而言的，由于较低的声速、较高的相对带宽等问题使得SAS的处理过程实际上有别于SAR，目前已有借鉴双基地雷达的二维频谱的解析解并利用MSR方法，采用方位谱扩展的RDA算法进行成像的方法，但是该算法仅适用于窄带、窄测绘带等条件下。
发明内容：
为了克服上述背景技术的缺陷，本发明提供一种基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法。
为了解决上述技术问题本发明的所采用的技术方案为：
基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法：
步骤1，获取N个子阵的回波；
步骤2，对单个子阵的回波数据进行二维快速傅立叶变换，将回波数据变换到二维波数域，并求解单个子阵回波的二维波数谱解析解；
步骤3，对各子阵的方位向波数谱进行延展；
步骤4，矫正各个子阵的方位向时间偏移；
步骤5，消除数据基线长度对二维波数谱的影响，同时完成参考函数相乘和距离向脉冲压缩；
步骤6，基于非均匀傅里叶变换进行Stolt映射及距离向IFFT；
步骤7，将步骤6所得结果数据进行相干叠加；
步骤8，对叠加后的距离多普勒域中的数据进行方位向的IFFT实现成像。
较佳地，步骤1中子阵i的回波可表示为
si(τ，η)＝A0wr(τ-Ri/c)wa(η)exp{-j2πf0Ri/c}exp{jπμ(τ-Ri/c)
其中A0为常数，wr(·)为距离向包络，wa(·)为方位向包络，τ为快时间，Ri为发射阵发射脉冲到达接收子阵i时历经的距离，c为声速，η为慢时间，f0为发射的线性调频信号的载波频率，μ为线性调频信号的调频率。在非停走停假设下， 其中xT为发射阵发射时刻所在的方位位置，xp为目标方位位置，r为目标零多普勒距离；Δdi为发射阵相位中心与第i个接收子阵相位中心的距离，V为合成孔径声纳运动速度，τi
较佳地，步骤2中第i个子阵的回波的二维波数谱解析解为：
其中， fτ为发射信号基带频率，f0为发射信号载频，c为波速，fa为方位向多普勒频率，μ为发射信号调频率，Wr(kr)表示发射信号谱包络形状，wa(kx)表示方位谱包络形状，A为常量，RB为目标的零多普勒距离，xn为目标方位向位置，β0是收发分置半角β在发射波束中心扫过条带区域中心点目标时的接收子阵的数据基线长度hic和条带区域中心距离RB＝RBc时的值，收发分置半角β的解析解为
其中hi为子阵i的数据基线长度。
较佳地，步骤3，对方位向波数谱进行复制延展得到：
其中，Kx_M为方位频谱延展后的方位向波数，方位频谱延展的倍数为不小于由方位向多普勒带宽与脉冲重复频率的比值的整数。
较佳地，步骤4，矫正各个子阵的方位向时间轴的方法包括：
将步骤4所得结果乘以 将各个子阵的方位向时间校正到设定的方位向0时刻，校正后回波的二维波数谱解析解为
较佳地，步骤5，消除数据基线长度对二维波数谱的影响，同时完成参考函数相乘和距离向脉冲压缩的方法包括：
将步骤4所得结果乘以
得到回波的二维波数谱解析解为
较佳地，步骤6对步骤5所得结果进行基于NUFFT的Stolt映射，得到经过该映射后回波变化到距离多普勒域： 其中A0为常量，Br为发射信号带宽。
较佳地，步骤7对步骤6所得的各个子阵的结果进行相干叠加。
较佳地，对步骤8对步骤7所得的结果进行方位向逆傅立叶变换，变换后的结果即为成像结果。
本发明的有益效果在于：本发明构建了多接收子阵SAS在非停走停假设下的几何模型，推导了适用于多接收子阵SAS成像的精确的二维波数谱解析解，提出了一种与与该二维波数谱相结合的、基于NUFFT的多接收子阵SAS的ωkA成像算法，最后通过仿真实验和湖试数据成像结果验证了其正确性和有效性。
附图说明
图1为本发明实施例的流程图；
图2为本发明实施例多子阵SAS成像几何模型；
图3为本发明实施例斜距平面上点目标分布示意图；
图4为本发明实施例RBc＝85m时成像结果；
图5为本发明实施例RBc＝45m时成像结果；
图6为本发明实施例所述方法进行湖试实测数据成像结果；
图7为采用RDA成像算法进行湖试实测数据成像结果。
具体实施方式
下面结合附图和实施例对本发明做进一步的说明，如图1所示，包括如下步骤：
步骤1：获取N个子阵的回波，回波的几何模型如图2所示。
图2为多子阵SAS几何模型。为斜距平面，x轴为方位轴，r轴为距离轴。声纳沿x轴正方向以速度V运动，发射阵在前接收子阵在后，发射阵及各接收子阵孔径宽度均为La且紧密相邻。θT为发射阵到目标的视角，RT为发射阵到目标的距离，θR为接收阵到目标的视角，RR为接收子阵到目标的距离，条带成像区域中心距离为RB，β为收发分置半角， 设定发射阵波束中心扫过条带成像区域中心点目标P(Xn，Yn)的时刻为方位向0时刻。
在非停走停假设下发射阵发射的信号经过目标反射回接收阵期间，接收阵将发生移动，定义发射阵发射时刻所在的位置与第i个接收子阵接收时刻所在的位置之间的间距hi为“数据基线”的长度。设数据基线中点为Mi，Δdi是发射阵相位中心与第i个接收子阵相位中心之间的距离，若发射阵发射的信号被第i个接收子阵接收到的时间为τi
其中xp是目标方位位置，r为目标零多普勒距离，c为声速。
步骤2对单个子阵的回波数据进行二维快速傅立叶变换，将回波数据变换到二维波数域，并求解单个子阵回波的二维波数谱解析解。
由于在非停走停假设下hi是变化的，因此对τi
在非停走停假设下，单个子阵回波的二维波数谱的解析解为
其中 fr为发射信号基带频率，f0为发射信号载频，c为波速，fa为方位向多普勒频率，μ为发射信号调频率。Wr(kr)表示发射信号谱包络形状，Wa(kx)表示方位谱包络形状，A为常量。
β的解析解为：
选择发射波束中心扫过条带区域中心点目标时的接收子阵的数据基线长度hic和条带区域中心距离RBc来对β近似，设此时β值为β0。对(2)式在RB＝RBc处做1阶泰勒级数展开，即可得第i个子阵的回波的二维波数谱解析解为
步骤3对方位向波数谱进行延展。
对多子阵SAS来说，单个子阵的方位谱存在严重的多普勒频谱卷绕，如果直接处理，则卷绕的频谱将严重影响成像质量。因此先需要对方位谱进行延展以解决多普勒频谱卷绕问题，方位频谱延展的倍数为大于方位向多普勒带宽与脉冲重复频率比值的最小整数。延展后回波的二维波数谱解析解为
其中Kx_M为方位频谱延展后的方位向波数。
步骤4各个子阵的方位向时间轴的矫正。
将(10)乘以 将各个子阵的方位向时间校正到设定的方位向0时刻。校正后回波的二维波数谱解析解为
步骤5参考函数相乘。
将(11)式乘以 H2中的第一项实现距离向脉压，第二项消除数据基线对二维波数谱的影响，并且以条带区域中心距离RBc为参考完成一致相位去除。参考函数相乘后回波的二维波数谱解析解为
步骤6：基于非均匀傅里叶变换进行Stolt映射。
则经过Stolt映射后的二维波数谱为
G(Kx，KR)＝AWr(Kr)Wa(Kx)exp{-j(RB-RBc)Ky-jKx_MXn)} (12)
即完成 得映射，相应的频域的映射关系为 其中fτ′为映射后的新的频率轴上的频率，fτ′各频点之间的间隔是非均匀的，此时使用NUFFT直接完成Stolt插值以及距离向傅立叶逆变换(IFFT)的计算。计算完成后回波数据变换到距离多普勒域，
步骤7：将步骤6所得结果数据的相干叠加，相干叠加等价于提高了方位向的采样频率，用于消除多普勒频谱混叠的影响。根据傅立叶变换是线性变换的特点，该叠加也可以在方位向IFFT之后进行。
步骤8：方位向逆傅立叶变换实现成像。对叠加后的距离多普勒域中的数据进行方位向的IFFT从而实现成像。
仿真数据成像结果如图4、图5所示，实测数据成像结果如图6所示。设发射信号载频为150kHz，带宽40kHz，脉冲宽度20ms，发射脉冲重复间隔为0.2s，平台的运动速度为2.5m/s，发射和接收阵元孔径均为0.04m，接收子阵个数为25个，条带区域距离向宽度为30m，条带中心距离为RBc。成像区域中心点有1个目标T3，其它目标相对于中心目标在距离上有±3m的偏移，在方位上有±1m的偏移，如图3所示。
图4给出了采用仿真的多子阵数据进行成像后的结果。其中图4、图5是本算法结果，计算过程中NUFFT的计算误差设定为10
图6和图7给出了ChiSAS-150多子阵SAS系统在千岛湖湖底采集的数据进行成像的结果，通过与参考算法的成像结果对比可知，在细节上图6好于图7，说明本文提出的算法在在实际应用中有较好的适用性。
本实施例以基于HQBA的二维波数谱位基础，首先进行方位谱延展，然后完成方位时间校正、消除数据基线长度对二维波数谱的影响、参考函数相乘以及距离向脉冲压缩，随后完成基于NUFFT(NUFFT，Nonuniform FFT)的Stolt插值，最后对插值后多子阵数据的相干叠加并进行方位向逆傅立叶变换，实现最终成像。
应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法专利购买费用说明
Q：办理专利转让的流程及所需资料
A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。
1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。
2：按规定缴纳著录项目变更手续费。
3：同时提交相关证明文件原件。
4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。
Q:专利著录项目变更费用如何缴交
A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式
Q：专利转让变更，多久能出结果
A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。
动态评分
0.0