发明名称 一种基于高精度插值的双基地合成孔径雷达成像方法
摘要 本发明公开了一种基于高精度插值的双基地合成孔径雷达成像方法,它是针对当前的双基地合成孔径雷达后向投影成像方法插值精度低的缺点,通过将高精度的非均匀快速傅里叶变换方法应用到时域成像插值步骤中,实现了高精度双基地合成孔径雷达成像。与现有的方法相比,本发明方法既保持了时域成像方法适用性广、内存开销小、易实施运动补偿等优点,而且能够获得更高质量的成像结果。
申请公布号 CN103529448A 申请公布日期 2014.01.22
申请号 CN201310488448.8 申请日期 2013.10.18
申请人 电子科技大学 发明人 刘喆;张晓玲;于淼
分类号 G01S13/90(2006.01)I;G01S7/41(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 曾磊
主权项 1.一种基于高精度插值的双基地合成孔径雷达成像方法,其特征是它包括如下步骤:步骤1:初始化参数初始化成像系统参数包括:雷达系统工作的信号波长,记做λ,雷达系统发射信号带宽,记做B,雷达系统发射脉冲时宽,记做T<sub>r</sub>,雷达系统采样频率,记做F<sub>s</sub>,雷达系统脉冲重复频率,记做PRF,雷达发射平台速度矢量,记做<img file="FDA0000397752700000011.GIF" wi="82" he="84" />雷达接收平台速度矢量,记做<img file="FDA0000397752700000012.GIF" wi="96" he="84" />雷达发射平台初始位置矢量,记做<img file="FDA0000397752700000013.GIF" wi="99" he="84" />雷达接收平台初始位置矢量,记做<img file="FDA0000397752700000014.GIF" wi="105" he="84" />场景参考点位置矢量,记做<img file="FDA0000397752700000015.GIF" wi="86" he="84" />雷达系统距离向采样点数,记做N<sub>r</sub>,雷达系统方位向采样点数,记做N<sub>a</sub>,光的传播速度,记做C;初始化场景参数:图像距离向像素点间隔,记做d<sub>r</sub>,图像方位向像素点间隔,记做d<sub>a</sub>,图像距离向像素点数,记做s<sub>r</sub>,图像方位向像素点数,记做s<sub>a</sub>;双基地合成孔径雷达原始回波数据,记做<img file="FDA0000397752700000016.GIF" wi="68" he="73" />是一个二维的数据矩阵,时刻雷达原始回波数据<img file="FDA0000397752700000017.GIF" wi="40" he="72" />数据矩阵每列数据是方位向回波信号的采样,每行的数据是逐个单脉冲距离向回波信号的采样;步骤1中所述参数均为已知;步骤2:距离压缩对步骤1中的双基地合成孔径雷达原始回波数据<img file="FDA0000397752700000018.GIF" wi="44" he="72" />采用传统的合成孔径雷达标准距离压缩方法进行压缩,得到距离压缩后的数据<img file="FDA0000397752700000019.GIF" wi="411" he="82" />是一个二维矩阵,τ代表矩阵第一维索引,t代表矩阵第二维索引;步骤3:计算天线相位中心位置利用公式<maths num="0001"><![CDATA[<math><mrow><mover><msub><mi>P</mi><mi>t</mi></msub><mo>-</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mover><msub><mi>P</mi><mrow><mi>t</mi><mn>0</mn></mrow></msub><mo>-</mo></mover><mo>+</mo><mover><msub><mi>V</mi><mi>t</mi></msub><mo>-</mo></mover><mo>&CenterDot;</mo><mi>n</mi><mo>/</mo><mi>PRF</mi></mrow></math>]]></maths>和公式<maths num="0002"><![CDATA[<math><mrow><mover><msub><mi>P</mi><mi>r</mi></msub><mo>-</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mover><msub><mi>P</mi><mrow><mi>r</mi><mn>0</mn></mrow></msub><mo>-</mo></mover><mo>+</mo><mover><msub><mi>V</mi><mi>r</mi></msub><mo>-</mo></mover><mo>&CenterDot;</mo><mi>n</mi><mo>/</mo><mi>PRF</mi></mrow></math>]]></maths>分别计算得到雷达发射、接收平台在第n个PRF时刻的天线相位中心矢量<img file="FDA00003977527000000112.GIF" wi="118" he="84" />以及<img file="FDA00003977527000000113.GIF" wi="147" he="84" />其中,<img file="FDA00003977527000000114.GIF" wi="202" he="84" />分别是雷达发射、接收平台初始位置矢量,<img file="FDA00003977527000000115.GIF" wi="162" he="98" />分别是雷达发射、接收平台速度矢量,PRF是雷达系统脉冲重复频率,n=1,...,N<sub>a</sub>,n表示第n个PRF时刻,N<sub>a</sub>是雷达系统方位向采样点数;步骤4:计算目标位置采用公式<img file="FDA00003977527000000116.GIF" wi="795" he="79" />计算得到图像中像素点P(a,r)的位置矢量<img file="FDA0000397752700000021.GIF" wi="171" he="77" />其中,<img file="FDA0000397752700000022.GIF" wi="55" he="79" />是场景参考点位置矢量,<img file="FDA0000397752700000023.GIF" wi="53" he="65" />和<img file="FDA0000397752700000024.GIF" wi="50" he="65" />是合成孔径雷达成像空间的坐标基,r表示像素点位于图像距离向的第r个位置,r=1,...,s<sub>r</sub>,s<sub>r</sub>为子图像距离向的像素点数,a表示像素点位于图像方位向的第a个位置,a=1,...,s<sub>a</sub>,s<sub>a</sub>为图像方位向的像素点数,d<sub>r</sub>是图像距离向像素点间隔,d<sub>a</sub>是图像方位向像素点间隔;步骤5:计算距离门利用步骤3中得到的雷达发射、接收平台在第n个PRF时刻的天线相位中心矢量<img file="FDA0000397752700000025.GIF" wi="112" he="84" />以及<img file="FDA0000397752700000026.GIF" wi="121" he="83" />以及步骤4中得到的图像中像素点P(a,r)的位置矢量<img file="FDA0000397752700000027.GIF" wi="180" he="77" />采用公式ID(n,a,r)=R(n,a,r)·F<sub>s</sub>/C得到双基地雷达距离史R(n,a,r)对应的距离门ID(n,a,r);其中,F<sub>s</sub>是雷达系统采样频率,C是光的传播速度,R(n,a,r)为双基地合成孔径雷达距离史,由公式<maths num="0003"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>a</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><msub><mrow><mo>|</mo><mo>|</mo><mover><msub><mi>P</mi><mi>t</mi></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>P</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>a</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub><mo>+</mo><msub><mrow><mo>|</mo><mo>|</mo><mover><msub><mi>P</mi><mi>r</mi></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>P</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>a</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub></mrow></math>]]></maths>得到,||·||<sub>2</sub>为L<sub>2</sub>范数,a表示像素点位于图像方位向的第a个位置,a=1,...,s<sub>a</sub>,s<sub>a</sub>为图像方位向的像素点数,r表示像素点位于图像距离向的第r个位置,r=1,...,s<sub>r</sub>,s<sub>r</sub>为图像距离向的像素点数;步骤6:高精度插值依据步骤5中得到的距离门ID(n,a,r),采用一维非均匀快速傅里叶变换方法对由步骤2得到的距离压缩数据矩阵<img file="FDA0000397752700000029.GIF" wi="186" he="82" />进行高精度插值,得到插值后的数据<img file="FDA00003977527000000210.GIF" wi="400" he="82" />其中ID(n,a,r)为距离门,n表示第n个PRF时刻,n=1,...,N<sub>a</sub>,a和r表示像素点位于成像空间中方位向的第a个位置、距离向第r个位置;步骤7:相位补偿将步骤6中得到的插值后的数据<img file="FDA00003977527000000211.GIF" wi="368" he="82" />与补偿相位因子K(n,a,r)相乘,得到相位补偿后的数据A(n,a,r),其中,补偿相位因子K(n,a,r)为<img file="FDA00003977527000000212.GIF" wi="695" he="145" />j为虚数单位(即-1开根),R(n,a,r)为步骤4中得到的距离史,λ为雷达系统工作的信号波长,n表示第n个PRF时刻;步骤8:相干积累利用公式<img file="FDA00003977527000000213.GIF" wi="384" he="124" />将步骤7得到的相位补偿后的数据A(n,a,r)累加起来得到最后的成像结果<img file="FDA0000397752700000031.GIF" wi="141" he="72" />其中,n表示第n个PRF时刻,a和r表示像素点位于成像空间中方位向的第a个位置、距离向第r个位置。
地址 611731 四川省成都市高新区(西区)西源大道2006号