发明名称 一种基于级数反演的机载圆迹环扫SAR成像方法
摘要 本发明涉及一种基于级数反演的机载圆迹环扫SAR成像方法,采用的技术方案是:a.首先对回波信号做二维FFT得到二维频谱;b.得到二维频谱对距离徙动项,距离调制项,残余相位项和距离调制项进行补偿;c.作距离IFFT,将信号变换至距离多普勒域并乘以方位匹配滤波函数完成方位脉冲压缩;d.最后对信号作方位IFFT就可以得到聚焦后的SAR图像。在平台飞行高度较低、俯仰角不大、积累时间较长的情况下能够获得较好的成像效果。
申请公布号 CN102928841A 申请公布日期 2013.02.13
申请号 CN201210423924.3 申请日期 2012.10.30
申请人 西安电子科技大学 发明人 邢孟道;廖轶;张磊
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 代理人
主权项 1.一种基于级数反演的机载圆迹环扫SAR成像方法,雷达平台在距离目标所在平面高度为H<sub>C</sub>的平面上以半径r<sub>a</sub>做圆形轨迹运动,平台运动的角速度为ω,平台旋转运动的方位向慢时间(方位向慢时间表示相邻两波束发射间隔时间)为η,雷达平台运动时,波束的指向始终与飞行速度方向垂直,目标点P离雷达平台的运动圆形轨迹在目标点所在平面的投影的中心点距离为r<sub>p</sub>,雷达发射信号为线性调频LFM信号,脉冲宽度为T<sub>p</sub>,调频率为γ,σ<sub>p</sub>为散射系数,a<sub>r</sub>(t)和a<sub>a</sub>(η)分别为雷达线性调频LFM信号的距离向时域包络函数和方位向时域包络函数,c为光速,λ为雷达线性调频LFM信号的中心频率所对应的波长,f<sub>c</sub>为载频,f<sub>τ</sub>为距离向频率,f<sub>η</sub>为多普勒频率;其特征在于成像方法按以下步骤完成:a).计算<maths num="0001"><![CDATA[<math><mrow><msub><mi>R</mi><mi>cen</mi></msub><mo>=</mo><msqrt><msubsup><mi>H</mi><mi>c</mi><mn>2</mn></msubsup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mi>p</mi></msub><mo>-</mo><msub><mi>r</mi><mi>a</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>;</mo></mrow></math>]]></maths>b).计算<maths num="0002"><![CDATA[<math><mrow><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>r</mi><mi>a</mi></msub><msub><mi>r</mi><mi>p</mi></msub><msup><mi>&omega;</mi><mn>2</mn></msup></mrow><msub><mi>R</mi><mi>cen</mi></msub></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>k</mi><mn>4</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><msub><mi>r</mi><mi>a</mi></msub><msub><mi>r</mi><mi>p</mi></msub><msup><mi>&omega;</mi><mn>4</mn></msup></mrow><mrow><mn>12</mn><msub><mi>R</mi><mi>cen</mi></msub></mrow></mfrac><mo>-</mo><mfrac><mrow><msubsup><mi>r</mi><mi>a</mi><mn>2</mn></msubsup><msubsup><mi>r</mi><mi>p</mi><mn>2</mn></msubsup><msup><mi>&omega;</mi><mn>4</mn></msup></mrow><mrow><mn>4</mn><msubsup><mi>R</mi><mi>cen</mi><mn>3</mn></msubsup></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>c).用FFT求a<sub>r</sub>(t)的时域向频域变换得距离向频率包络A<sub>r</sub>(f<sub>τ</sub>),A<sub>r</sub>(f<sub>τ</sub>)=FFT[a<sub>r</sub>(t)];d).用FFT求a<sub>a</sub>(η)的时域向频域变换得方位向频率包络A<sub>a</sub>(f<sub>η</sub>),A<sub>a</sub>(f<sub>η</sub>)=FFT[a<sub>a</sub>(η)];e).计算<img file="FSA00000797081600014.GIF" wi="772" he="152" /><img file="FSA00000797081600015.GIF" wi="593" he="132" /><img file="FSA00000797081600016.GIF" wi="294" he="113" /><img file="FSA00000797081600017.GIF" wi="705" he="133" /><img file="FSA00000797081600018.GIF" wi="1076" he="135" /><img file="FSA00000797081600019.GIF" wi="619" he="135" />f).计算<img file="FSA000007970816000110.GIF" wi="1604" he="114" />g).计算S<sub>A</sub>(f<sub>τ</sub>,f<sub>η</sub>)=σ<sub>p</sub>A<sub>r</sub>(f<sub>τ</sub>)A<sub>a</sub>(f<sub>η</sub>)·exp(jθ<sub>c</sub>);h).计算相位补偿项<maths num="0004"><![CDATA[<math><mrow><msub><mi>H</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>&tau;</mi></msub><mo>,</mo><msub><mi>f</mi><mi>&eta;</mi></msub><mo>,</mo><msub><mi>R</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mi>j</mi><mfrac><mrow><mi>&pi;</mi><mo>&CenterDot;</mo><mi>c</mi></mrow><mrow><mn>4</mn><msub><mi>k</mi><mn>2</mn></msub></mrow></mfrac><msup><mrow><mo>(</mo><mfrac><msub><mi>f</mi><mi>&eta;</mi></msub><msub><mi>f</mi><mi>c</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mi>j</mi><mfrac><mrow><mn>3</mn><mi>&pi;</mi><msub><mi>k</mi><mn>4</mn></msub><msup><mi>c</mi><mn>3</mn></msup></mrow><mrow><mn>64</mn><msubsup><mi>k</mi><mn>2</mn><mn>4</mn></msubsup></mrow></mfrac><msup><mrow><mo>(</mo><mfrac><msub><mi>f</mi><mi>&eta;</mi></msub><msub><mi>f</mi><mi>c</mi></msub></mfrac><mo>)</mo></mrow><mn>4</mn></msup><mo>-</mo><mi>j</mi><mfrac><mrow><mi>&pi;</mi><mo>&CenterDot;</mo><mi>c</mi><mo>&CenterDot;</mo><msubsup><mi>f</mi><mi>&eta;</mi><mn>2</mn></msubsup></mrow><mrow><mn>4</mn><msub><mi>k</mi><mn>2</mn></msub><msubsup><mi>f</mi><mi>c</mi><mn>3</mn></msubsup></mrow></mfrac></mtd></mtr><mtr><mtd><mo>+</mo><mi>j</mi><mfrac><mrow><mn>3</mn><mi>&pi;</mi><msub><mi>k</mi><mn>4</mn></msub><msup><mi>c</mi><mn>3</mn></msup><msubsup><mi>f</mi><mi>&eta;</mi><mn>4</mn></msubsup></mrow><mrow><mn>32</mn><msubsup><mi>k</mi><mn>2</mn><mn>4</mn></msubsup><msubsup><mi>f</mi><mi>c</mi><mn>5</mn></msubsup></mrow></mfrac><mo>+</mo><mi>j</mi><mfrac><mrow><mi>&pi;</mi><mo>&CenterDot;</mo><mi>c</mi><mo>&CenterDot;</mo><msubsup><mi>f</mi><mi>&eta;</mi><mn>2</mn></msubsup><msubsup><mi>f</mi><mi>&tau;</mi><mn>3</mn></msubsup></mrow><mrow><mn>4</mn><msub><mi>k</mi><mn>2</mn></msub><msubsup><mi>f</mi><mi>c</mi><mn>3</mn></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>c</mi></msub><mo>+</mo><msub><mi>f</mi><mi>&tau;</mi></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mo>+</mo><mi>j</mi><mfrac><mrow><mi>&pi;</mi><msub><mi>k</mi><mn>4</mn></msub><msup><mi>c</mi><mn>3</mn></msup><msubsup><mi>f</mi><mi>&eta;</mi><mn>4</mn></msubsup></mrow><mrow><mn>64</mn><msubsup><mi>k</mi><mn>2</mn><mn>4</mn></msubsup></mrow></mfrac><mo>[</mo><mn>1</mn><mo>/</mo><mrow><mo>(</mo><msub><mi>f</mi><mi>c</mi></msub><mo>+</mo><msub><mi>f</mi><mi>&tau;</mi></msub><mo>)</mo></mrow><mo>-</mo><mn>1</mn><mo>/</mo><msubsup><mi>f</mi><mi>c</mi><mn>3</mn></msubsup><mo>+</mo><mn>3</mn><msub><mi>f</mi><mi>&tau;</mi></msub><mo>/</mo><msubsup><mi>f</mi><mi>c</mi><mn>4</mn></msubsup><mo>-</mo><mn>6</mn><msubsup><mi>f</mi><mi>&tau;</mi><mn>2</mn></msubsup><mo>/</mo><msubsup><mi>f</mi><mi>c</mi><mn>5</mn></msubsup><mo>]</mo></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>i).计算S<sub>A</sub>(f<sub>τ</sub>,f<sub>η</sub>)=S<sub>A</sub>(f<sub>τ</sub>,f<sub>η</sub>)·H<sub>l</sub>(f<sub>τ</sub>,f<sub>η</sub>,R<sub>0</sub>);j).计算距离IFFT,将信号变换至距离多普勒域,G<sub>η</sub>(t,f<sub>η</sub>)=IFFT[S<sub>A</sub>(f<sub>τ</sub>,f<sub>η</sub>)];k).计算方位匹配滤波函数,<maths num="0005"><![CDATA[<math><mrow><msub><mi>H</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>&eta;</mi></msub><mo>,</mo><msub><mi>R</mi><mi>cen</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mi>&pi;</mi><mo>&CenterDot;</mo><mi>c</mi><mo>&CenterDot;</mo><msubsup><mi>f</mi><mi>&eta;</mi><mn>2</mn></msubsup></mrow><mrow><msub><mrow><mn>4</mn><mi>k</mi></mrow><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>R</mi><mi>cen</mi></msub><mo>)</mo></mrow><msub><mi>f</mi><mi>c</mi></msub></mrow></mfrac><mo>+</mo><mi>j</mi><mfrac><mrow><mi>&pi;</mi><mo>&CenterDot;</mo><msub><mi>k</mi><mn>4</mn></msub><mrow><mo>(</mo><msub><mi>R</mi><mi>cen</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>c</mi><mn>3</mn></msup><msubsup><mi>f</mi><mi>&eta;</mi><mn>4</mn></msubsup></mrow><mrow><mn>64</mn><msubsup><mi>k</mi><mn>2</mn><mn>4</mn></msubsup><mrow><mo>(</mo><msub><mi>R</mi><mi>cen</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msubsup><mi>f</mi><mi>c</mi><mn>3</mn></msubsup></mrow></mfrac><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>l).进行方位脉冲压缩,G<sub>η</sub>(t,f<sub>η</sub>)=H2(f<sub>η</sub>,R<sub>cen</sub>)G<sub>η</sub>(t,f<sub>η</sub>);m).计算方位IFFT就可以得到聚焦后的SAR图像,G(t,η)=IFFT[G<sub>η</sub>(t,f<sub>η</sub>)]。
地址 710071 陕西省西安市太白路2号