发明名称 一站固定式调频连续波双基地SAR成像方法
摘要 本发明公开一种一站固定式调频连续波双基地SAR成像算法,主要解决现有方法不能同时校正方位向距离徙动空变性和方位调频率空变性的问题。其实现步骤是:(1)获取目标的两维时域回波信号,并计算其距离频域方位时间域回波信号;(2)根据距离频域方位时间域回波信号校正方位向距离徙动空变性和方位调频率空变性;(3)对校正后的信号补偿剩余相位,并校正距离向距离徙动的空变性;(4)对距离向距离徙动空变性校正后的信号进行方位向脉冲压缩,并做方位向逆傅里叶变换得到目标时域图像。本发明具有能够同时校正方位向距离徙动空变性和方位调频率空变性以及实现效率高的优点,可用于一站固定式调频连续波双基地SAR的成像。
申请公布号 CN103454632A 申请公布日期 2013.12.18
申请号 CN201310242314.8 申请日期 2013.08.24
申请人 西安电子科技大学 发明人 廖桂生;夏艳杰;徐青;李东;王威;杨科
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一站固定式调频连续波双基地SAR成像方法,包括如下步骤:1)通过静态发射机向目标发出调频连续波的雷达信号,由动态的接收机接收雷达回波信号,并对该回波信号依次进行去调频和去剩余视频相位误差处理,得到目标的两维时域回波信号s(t<sub>r</sub>,t<sub>a</sub>),其中,t<sub>r</sub>是距离时间,t<sub>a</sub>是方位时间;2)将经过预处理的雷达回波信号s(t<sub>r</sub>,t<sub>a</sub>)用f<sub>r</sub>等效代换t<sub>r</sub>,再根据泰勒级数展开法展开化简,得到距离频域方位时间域的回波信号s<sub>1</sub>(f<sub>r</sub>,t<sub>a</sub>),其中,f<sub>r</sub>是距离频率;3)根据目标到雷达接收机和发射机的距离,构造空变性校正函数<img file="FDA00003366264000011.GIF" wi="460" he="84" />以对雷达回波信号中方位向距离徙动的空变性和方位调频率的空变性进行校正,得到校正后的新的距离频域方位时间域的回波信号:s<sub>2</sub>(f<sub>r</sub>,t<sub>a</sub>)=s<sub>1</sub>(f<sub>r</sub>,t<sub>a</sub>)·H<sub>1</sub>(t<sub>a</sub>)其中,<img file="FDA00003366264000012.GIF" wi="335" he="151" />是求解得到的空变性校正函数的系数,v是接收机的运动速度,<img file="FDA00003366264000013.GIF" wi="292" he="135" />是等效波长,f<sub>c</sub>是中心频率,<img file="FDA00003366264000014.GIF" wi="449" he="87" />是经过固定的发射机调制后的新的距离,t<sub>n</sub>是目标的方位时间,d<sub>0R</sub>和d<sub>0T</sub>分别是目标到接收机和发射机的最短距离;4)对空变性校正后的雷达回波信号进行方位向的傅里叶变换,得到两维频域回波信号:G(f<sub>r</sub>,f<sub>a</sub>)=FFT<sub>a</sub>{s<sub>2</sub>(f<sub>r</sub>,t<sub>a</sub>)},其中,f<sub>a</sub>是方位频率;5)根据步骤4)得到的两维频域回波信号G(f<sub>r</sub>,f<sub>a</sub>),构造剩余相位补偿函数H<sub>2</sub>(f<sub>r</sub>,f<sub>a</sub>):<maths num="0001"><![CDATA[<math><mrow><msub><mi>H</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>r</mi></msub><mo>,</mo><msub><mi>f</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><msub><mi>&pi;cd</mi><mrow><mn>0</mn><mi>Rref</mi></mrow></msub><msubsup><mi>f</mi><mi>a</mi><mn>2</mn></msubsup></mrow><mrow><msup><mi>v</mi><mn>2</mn></msup><msubsup><mi>f</mi><mi>c</mi><mn>2</mn></msubsup></mrow></mfrac><msub><mi>f</mi><mi>r</mi></msub><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mfrac><mrow><msub><mi>&pi;cd</mi><mrow><mn>0</mn><mi>Rref</mi></mrow></msub><msubsup><mi>f</mi><mi>a</mi><mn>2</mn></msubsup></mrow><mrow><msup><mi>v</mi><mn>2</mn></msup><msubsup><mi>f</mi><mi>c</mi><mn>3</mn></msubsup></mrow></mfrac><msubsup><mi>f</mi><mi>r</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mi>exp</mi><mo>[</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>a</mi></msub><mrow><mo>(</mo><mfrac><msub><mi>f</mi><mi>r</mi></msub><msub><mi>K</mi><mi>r</mi></msub></mfrac><mo>+</mo><mfrac><msub><mi>d</mi><mn>0</mn></msub><mi>c</mi></mfrac><mo>)</mo></mrow><mo>]</mo><mo>,</mo></mrow></math>]]></maths>其中,c是光速,f<sub>c</sub>是中心频率,K<sub>r</sub>是发射信号的调频率,d<sub>0</sub>=d<sub>0Rref</sub>+d<sub>0Tref</sub>,是去调频时使用的参考距离,d<sub>0Rref</sub>和d<sub>0Tref</sub>分别是场景中心到接收机和发射机的最短距离;6)利用剩余相位补偿函数H<sub>2</sub>(f<sub>r</sub>,f<sub>a</sub>)补偿雷达回波信号中的距离方位耦合相位误差和由雷达接收机的运动引起的多普勒频移误差,得到新的两维频域回波信号G<sub>1</sub>(f<sub>r</sub>,f<sub>a</sub>)为:G<sub>1</sub>(f<sub>r</sub>,f<sub>a</sub>)=G(f<sub>r</sub>,f<sub>a</sub>)·H<sub>2</sub>(f<sub>r</sub>,f<sub>a</sub>);7)对经过剩余相位误差补偿后的雷达回波信号G<sub>1</sub>(f<sub>r</sub>,f<sub>a</sub>)做逆Chirp-z变换,得到去除距离向距离徙动空变性的回波信号s<sub>3</sub>(t<sub>r</sub>,f<sub>a</sub>);8)根据去除距离向距离徙动空变性的回波信号s<sub>3</sub>(t<sub>r</sub>,f<sub>a</sub>)构造方位压缩函数<img file="FDA00003366264000021.GIF" wi="636" he="165" />以对雷达回波信号的方位向脉冲进行压缩,得到脉冲压缩以后的雷达回波信号:s<sub>4</sub>(t<sub>r</sub>,f<sub>a</sub>)=s<sub>3</sub>(t<sub>r</sub>,f<sub>a</sub>)·H<sub>3</sub>(f<sub>a</sub>);9)对脉冲压缩以后的雷达回波信号s<sub>4</sub>(t<sub>r</sub>,f<sub>a</sub>)做方位向逆傅里叶变换,得到时域成像结果,并对该时域成像结果进行插值校正,获得目标图像。
地址 710071 陕西省西安市太白南路2号