发明名称 一种基于Keystone变换的地基合成孔径雷达快速成像方法
摘要 本发明公开了一种基于Keystone变换的地基合成孔径雷达快速成像方法。使用本发明能够适应GB SAR形变监测系统目标区域距离跨度大、多普勒中心空变性严重、方位分辨率随距离扩大等成像处理问题,并能满足系统实时性需求。发明是一种基于Keystone变换的方位Dechirp成像处理算法,相比于现有的经典SAR成像算法,本发明充分考虑了GB SAR信号模型特点,更适用于GB SAR实时系统成像处理任务。
申请公布号 CN104330795B 申请公布日期 2017.01.11
申请号 CN201410559229.9 申请日期 2014.10.20
申请人 北京理工大学 发明人 胡程;曾涛;毛聪;龙腾;田卫明;毛二可
分类号 G01S13/90(2006.01)I;G06F19/00(2011.01)I 主分类号 G01S13/90(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 付雷杰;杨志兵
主权项 一种基于楔石形Keystone变换的地基合成孔径雷达快速成像方法,其特征在于,包括如下步骤:步骤1,对雷达回波信号进行距离向处理:将回波信号进行距离压缩,在距离频域对距离压缩后的信号进行楔石形Keystone变换,完成距离徙动校正;步骤2,方位向处理:步骤2.1,将完成距离徙动校正后的信号变换到距离多普勒域,将多普勒域等效为方位角正弦域t‑sinθ,然后对信号进行方位分块:将信号沿整个sinθ轴划分为N<sub>sub</sub>个宽度为Δ<sub>sinθ</sub>的子块,各子块的中心位于(sinθ)<sub>i</sub>,i=1,2,…,N<sub>sub</sub>,其中,N<sub>sub</sub>为正整数,Δ<sub>sinθ</sub>为满足下式的最大值:sinθ=Δ<sub>sinθ</sub>〃N<sub>sub</sub><maths num="0001"><math><![CDATA[<mrow><msub><mi>&Delta;</mi><mrow><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi></mrow></msub><mo>&le;</mo><msqrt><mrow><mfrac><mrow><mn>2</mn><msub><mi>&lambda;</mi><mi>c</mi></msub><msub><mi>&rho;</mi><mi>min</mi></msub></mrow><msup><mi>L</mi><mn>2</mn></msup></mfrac><mo>+</mo><mn>4</mn><mo>|</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi><msubsup><mo>|</mo><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mn>2</mn></msubsup></mrow></msqrt><mo>-</mo><mn>2</mn><mo>|</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi><msub><mo>|</mo><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mo>|</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi><msubsup><mo>|</mo><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mn>2</mn></msubsup><mo>)</mo><mi>L</mi></mrow><msub><mi>&rho;</mi><mi>min</mi></msub></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0001073213140000011.GIF" wi="1214" he="166" /></maths>其中,λ<sub>c</sub>为发射信号载波对应的波长,ρ<sub>min</sub>为最小成像距离,L是雷达合成孔径长度,|sinθ|<sub>max</sub>表示场景内全部目标的方位角正弦值的绝对值|sinθ|中的最大值;步骤2.2,将各子块的信号变换到方位时域,各子块在二维时域分别与相应的参考信号进行去倾Dechirp处理,其中,参考信号d<sub>a</sub>为<maths num="0002"><math><![CDATA[<mrow><msub><mi>d</mi><mi>a</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><msub><mi>x</mi><mi>n</mi></msub><mo>;</mo><msub><mrow><mo>(</mo><mrow><mi>sin</mi><mi>&theta;</mi></mrow><mo>)</mo></mrow><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>r</mi><mi>e</mi><mi>c</mi><mi>t</mi><mrow><mo>(</mo><mfrac><msub><mi>x</mi><mi>n</mi></msub><mi>L</mi></mfrac><mo>)</mo></mrow><mi>exp</mi><mo>&lsqb;</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mi>&lambda;</mi><mi>c</mi></msub></mfrac><mfrac><msubsup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mo>(</mo><mi>sin</mi><mi>&theta;</mi><mo>)</mo><mo>)</mo></mrow><mi>i</mi><mn>2</mn></msubsup><mrow><mi>t</mi><mo>&CenterDot;</mo><mi>c</mi><mo>/</mo><mn>2</mn></mrow></mfrac><msubsup><mi>x</mi><mi>n</mi><mn>2</mn></msubsup><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001073213140000012.GIF" wi="1125" he="174" /></maths>其中,t为距离快时间,x<sub>n</sub>是方位采样位置,c是光速;步骤2.3,将各子块去倾Dechirp处理结果在二维时域相加,然后经方位快速傅里叶变换获得伪极坐标系下的SAR图像。
地址 100081 北京市海淀区中关村南大街5号