发明名称 一种合成孔径声纳快速实时成像方法
摘要 本发明提供一种合成孔径声纳快速实时成像方法,该方法根据运动误差查表获取声波从像素点到空间采样点的精确传播延时,并利用阵元单个空间采样点回波数据的反向传播结果与过程图像相干叠加,进行孔径合成,具有成像精确、计算量小、回波数据零库存和便于大规模流水线处理的优点,能够实现对海底的远距离、大场景的高分辨率快速实时成像。
申请公布号 CN103645479A 申请公布日期 2014.03.19
申请号 CN201310712769.1 申请日期 2013.12.20
申请人 中国船舶重工集团公司第七一五研究所 发明人 汲夏;朱必波;丛卫华
分类号 G01S15/89(2006.01)I 主分类号 G01S15/89(2006.01)I
代理机构 杭州九洲专利事务所有限公司 33101 代理人 陈健
主权项 1.一种合成孔径声纳快速实时成像方法,其特征在于:该方法利用阵元单个空间采样点所采集回波数据的反向传播结果,与缓存的过程图像相干叠加,实现快速实时合成孔径处理,其具体步骤为:步骤一:假设声纳平台运动轨迹为直线,预存声信号从阵元波束范围内各像素点到阵元的传播延时表<![CDATA[<math><mrow><mi>&tau;</mi><mo>=</mo><mo>[</mo><msub><mi>&tau;</mi><mrow><mi>k</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>]</mo><mo>&Element;</mo><msup><mi>Q</mi><mrow><msub><mi>L</mi><mi>SA</mi></msub><mo>&times;</mo><mi>N</mi></mrow></msup><mo>,</mo></mrow></math>]]></maths>其中<![CDATA[<math><mrow><msub><mi>&tau;</mi><mrow><mi>k</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>2</mn><msqrt><msubsup><mi>x</mi><mi>k</mi><mn>2</mn></msubsup><mo>+</mo><msup><msub><mi>r</mi><mi>n</mi></msub><mn>2</mn></msup></msqrt><mo>/</mo><mi>c</mi><mo>,</mo></mrow></math>]]></maths>τ<sub>k,n</sub>是像素点到阵元空间采样点的传播延时,x<sub>k</sub>是阵元空间采样点(x<sub>i</sub>,0)到像素点(x<sub>m</sub>,r<sub>n</sub>)的相对方位,即阵元在虚拟孔径中的方位位置,x<sub>k</sub>=x<sub>i</sub>-x<sub>m</sub>,k为相对方位对应的像素点个数,k=i-m,n为像素点与阵元相对距离对应的像素点个数,L<sub>SA</sub>为一个虚拟孔径长度内的方位向像素点个数,N为距离向像素点个数,c为水声声速,Q表示有理数集合;步骤二:参考阵元波束覆盖范围,开辟一个相同大小的过程图像缓存区B,<![CDATA[<math><mrow><mi>B</mi><mo>=</mo><mo>[</mo><msub><mi>b</mi><mrow><mi>k</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>]</mo><mo>&Element;</mo><msup><mi>C</mi><mrow><msub><mi>L</mi><mi>SA</mi></msub><mo>&times;</mo><mi>N</mi></mrow></msup><mo>,</mo></mrow></math>]]></maths>其中C表示复数集合;步骤三:回波数据接收并匹配滤波后,首先将光纤罗经和冗余相位中心算法相结合进行多源运动估计,并通过坐标系转换,计算出载体运动误差在空间采样点i处引入的阵元位置误差(Δx<sub>i</sub>,Δr<sub>i</sub>);步骤四:根据阵元位置误差(Δx<sub>i</sub>,Δr<sub>i</sub>),查表取传声播延时<img file="FDA0000442719890000014.GIF" wi="229" he="65" />作为像素点(x<sub>m</sub>,r<sub>n</sub>)到阵元空间采样点(x<sub>i</sub>,0)的传播延时,其中Δk<sub>i</sub>=round(Δx<sub>i</sub>/δ<sub>x</sub>)、Δn<sub>i</sub>=round(Δr<sub>i</sub>/δ<sub>r</sub>),δ<sub>x</sub>、δ<sub>r</sub>分别为方位向和距离向像素点间隔,round表示取整数;步骤五:根据传播延时,将阵元在空间采样点位置(x<sub>i</sub>,0)接收的回波数据e<sub>i</sub>(t)反向传播到各像素点,并与缓存的过程图像相干叠加,进行孔径合成处理,即<![CDATA[<math><mrow><mi>B</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>B</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>e</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>&tau;</mi><mrow><mi>k</mi><mo>+</mo><msub><mi>&Delta;k</mi><mi>i</mi></msub><mo>,</mo><mi>n</mi><mo>+</mo><msub><mi>&Delta;n</mi><mi>i</mi></msub></mrow></msub><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>步骤六:沿声纳航行反方向,将孔径合成处理结果搬移一个空间采样点间隔,即B(k+q,n)=B(k,n),此时缓存区会自动溢出像素点B(L<sub>SA</sub>+1,n)□B(L<sub>SA</sub>+q,n),自动溢出像素点均已完成虚拟孔径的合成处理,可以输出送显示,其中q为一个空间采样点间隔对应的方位向像素点个数;步骤七:保存缓存区的过程图像,接收下一空间采样点接收的回波数据e<sub>i+1</sub>(t),回到步骤三。
地址 310012 浙江省杭州市西湖区华星路96号