发明名称 多通道多子带滑动聚束模式SAR成像方法
摘要 本发明公开了一种多通道多子带滑动聚束模式SAR成像方法,主要解决步进频率多通道合成孔径雷达在滑动聚束模式下的方位重建和频带合成难度大的问题。其实现过程是:1.对多通道天线进行子孔径划分,并对每个子孔径接收的回波依次进行等效相位中心的误差补偿和距离向脉冲压缩;2.对距离向脉冲压缩后的信号进行多普勒频移操作,将多普勒谱移到多普勒基带;3.计算出空域滤波的权重系数;4.利用空域滤波权重系数和基带信号进行多普勒解模糊;5.对多普勒解模糊后信号进行频带合成,完成二维频谱的重建并成像。本发明避免了高分辨率宽测绘数据的巨大计算量,提高了宽测绘带上的目标高分辨率,可用于机载或星载平台SAR多通道多子带成像。
申请公布号 CN104865571A 申请公布日期 2015.08.26
申请号 CN201510155399.5 申请日期 2015.04.02
申请人 西安电子科技大学 发明人 孙光才;景国彬;张佳佳;邢孟道;保铮
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种多通道多子带滑动聚束模式SAR成像,包括如下步骤:(1)对步进频率多通道合成孔径雷达SF‑MIMO‑SAR的天线进行子孔径划分,并对子孔径接收的回波依次进行多通道天线等效相位中心的误差补偿和距离向脉冲压缩;(2)对距离向脉冲压缩后的信号进行多普勒频移操作,将多普勒谱搬移到多普勒基带;(3)对以上得到的多普勒基带进行方位快速傅里叶FFT处理,即将多普勒基带的第i子孔径信号转化到多普勒频域,并为了减少空域滤波器中的权重系数的计算量,继续对以上方位FFT处理后的结果按如下补偿函数H<sub>3,n</sub>(f<sub>a</sub>,m)进行补偿:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>H</mi><mrow><mn>3</mn><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>a</mi></msub><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mfrac><msubsup><mi>X</mi><mi>nm</mi><mo>&prime;</mo></msubsup><mi>&upsi;</mi></mfrac><mrow><mo>(</mo><msub><mi>f</mi><mi>a</mi></msub><mo>+</mo><msub><mi>f</mi><mi>dc</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>c</mi></msub><mo>,</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000692812280000011.GIF" wi="973" he="162" /></maths>式中,X′<sub>nm</sub>是第n和m两个通道的等效位置,υ是雷达运动速度,f<sub>a</sub>是方位多普勒频率,f<sub>c</sub>是距离频率,t<sub>i</sub>为第i子孔径的方位中心时间;(4)对补偿后的基带信号数据进行多普勒解模糊操作,计算空域滤波的权重系数矩阵<img file="FDA0000692812280000012.GIF" wi="788" he="103" />其中:w<sub>n</sub>(m)是第n个发射通道,第m个接收通道的回波的权重系数;(5)利用权重系数<img file="FDA0000692812280000013.GIF" wi="91" he="86" />与补偿后的基带信号相乘,即进行空域滤波处理,得到多普勒解模糊后的基带信号多普勒谱S<sub>2,n</sub>(t,f<sub>d</sub>):<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><msub><mi>S</mi><mrow><mn>2</mn><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>w</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mn>2</mn><mi>c</mi></mfrac><msub><mi>R</mi><mi>ref</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><msub><mi>W</mi><mi>azi</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>d</mi></msub><mo>+</mo><msub><mi>f</mi><mi>dc</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>c</mi></msub><mo>,</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mfrac><mi>X</mi><mi>&upsi;</mi></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&upsi;</mi></mfrac><msub><mi>R</mi><mi>b</mi></msub><msqrt><msup><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&upsi;</mi></mrow><mi>c</mi></mfrac><msub><mi>f</mi><mi>d</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>f</mi><mi>d</mi></msub><mo>+</mo><msub><mi>f</mi><mi>dc</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>c</mi></msub><mo>,</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000692812280000014.GIF" wi="1387" he="357" /></maths>式中:R<sub>ref</sub>是雷达作用的参考距离,w<sub>r</sub>(·)为距离向的窗函数,w<sub>azi</sub>(·)为方位向的窗函数,f<sub>dc</sub>是多普勒中心频率,f<sub>c</sub>是距离频率,t<sub>i</sub>为第i子孔径的方位中心时间,f<sub>d</sub>为方位多普勒频率;(6)对多普勒解模糊后相邻基带信号进行频带合成,依次在距离频谱上进行小数部分的频移,并将小数部分的频移后相邻基带信号转换到距离频域,进行相干叠加,完成二维频谱重建;(7)对重建后的二维频谱进行成像。
地址 710071 陕西省西安市太白南路2号