发明名称 弹载SAR子孔径前斜视高阶非线性调频变标成像方法
摘要 本发明属于雷达信号处理领域,公开了一种弹载SAR子孔径前斜视高阶非线性调频变标成像方法,其利用子孔径数据实现弹载SAR前斜视高阶非线性调频变标成像,可用于机载、弹载平台SAR成像。其主要步骤为:(1)对回波信号进行距离脉压和时域走动校正;(2)将信号变到二维频域进行频域徙动校正和二次压缩并补偿方位向的高次相位;(3)在方位频域引入高阶非线性调频变标扰动因子,校正多普勒调频率和方位向高次项的空变;(4)通过谱分析处理将图像聚焦在方位频域上。本发明解决了距离方位的解耦合以及时域校正走动而带来的方位聚焦深度问题,其方法能够适应不同的场景及高分辨率要求,也可用于地图测绘等领域。
申请公布号 CN103901428A 申请公布日期 2014.07.02
申请号 CN201410122328.0 申请日期 2014.03.28
申请人 西安电子科技大学 发明人 梁毅;李震宇;马倩;邢孟道;王虹现
分类号 G01S13/90(2006.01)I;G06F19/00(2011.01)I 主分类号 G01S13/90(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 1.一种弹载SAR子孔径前斜视高阶非线性调频变标成像方法,其弹载SAR工作在条带模式,弹载平台高度为H,以速度v沿X轴匀速直线飞行,θ<sub>0</sub>为波束射线指向的斜视角,R<sub>0</sub>为波束中心线扫过目标时的斜距,t<sub>m</sub>为方位慢时间,以弹载平台位于O点的时刻为方位慢时间的起点,R<sub>B</sub>为场景中某一点目标与弹载SAR的最近距离;点目标到弹载SAR的瞬时斜距为<maths num="0001"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>R</mi><mn>0</mn></msub><mi>cos</mi><msub><mi>&theta;</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>R</mi><mn>0</mn></msub><mi>sin</mi><msub><mi>&theta;</mi><mn>0</mn></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>vt</mi><mi>m</mi></msub><mo>-</mo><msub><mi>X</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>,</mo></mrow></math>]]></maths>其中X<sub>n</sub>=vt<sub>n</sub>,t<sub>n</sub>为天线波束中心穿越目标点的时刻;其特征在于,包括以下步骤:步骤1,弹载SAR发射的是线性调频信号,则其基频回波信号为S<sub>0</sub>(t<sub>r</sub>,t<sub>m</sub>),其中t<sub>r</sub>为距离快时间,t<sub>m</sub>为方位慢时间;步骤2,对基频回波信号S<sub>0</sub>(t<sub>r</sub>,t<sub>m</sub>)进行距离向傅里叶变换处理,得到距离频域方位时域信号S<sub>1</sub>(f<sub>r</sub>,t<sub>m</sub>),其中f<sub>r</sub>为距离频率;步骤3,对距离频域方位时域信号S<sub>1</sub>(f<sub>r</sub>,t<sub>m</sub>)进行距离向脉冲压缩处理,得到距离脉压信号S<sub>2</sub>(f<sub>r</sub>,t<sub>m</sub>);步骤4,对距离脉压信号S<sub>2</sub>(f<sub>r</sub>,t<sub>m</sub>)进行线性距离走动校正,得到线性距离走动校正后的信号S<sub>3</sub>(f<sub>r</sub>,t<sub>m</sub>);步骤5,采用驻定相位原理,对线性距离走动校正后的信号S<sub>3</sub>(f<sub>r</sub>,t<sub>m</sub>)作方位向傅里叶变换,得到二维频域信号S<sub>4</sub>(f<sub>r</sub>,f<sub>a</sub>),其中f<sub>a</sub>为方位频率;步骤6,对二维频域信号S<sub>4</sub>(f<sub>r</sub>,f<sub>a</sub>)进行距离弯曲校正和二次距离压缩,得到解耦后的二维频域信号S<sub>5</sub>(f<sub>r</sub>,f<sub>a</sub>),然后对解耦后的二维频域信号S<sub>5</sub>(f<sub>r</sub>,f<sub>a</sub>)进行距离向逆傅里叶变换,得到距离时域方位频域信号S<sub>6</sub>(t<sub>r</sub>,f<sub>a</sub>);步骤7,对距离时域方位频域信号S<sub>6</sub>(t<sub>r</sub>,f<sub>a</sub>)在f<sub>a</sub>=0处进行方位向五阶泰勒级数展开,并对方位向四阶、五阶相位进行补偿,得到补偿完方位向四阶、五阶相位的距离时域方位频域信号S<sub>9</sub>(t<sub>r</sub>,f<sub>a</sub>);步骤8,对补偿完方位向四阶、五阶相位的距离时域方位频域信号S<sub>9</sub>(t<sub>r</sub>,f<sub>a</sub>)作方位向逆傅里叶变换,得到补偿完方位向四阶、五阶相位的时域信号S<sub>10</sub>(t<sub>r</sub>,t<sub>m</sub>),<img file="FDA0000483718310000021.GIF" wi="1481" he="202" />其中函数<img file="FDA0000483718310000022.GIF" wi="422" he="140" />B<sub>r</sub>为发射信号的频带,w<sub>a</sub>(·)为方位窗函数,t<sub>r</sub>为距离快时间,t<sub>m</sub>为方位慢时间,c为光速,exp(·)为指数函数,t<sub>n</sub>为天线波束中心穿越目标点的时刻,R<sub>0</sub>为波束中心线扫过目标时的斜距,<img file="FDA0000483718310000025.GIF" wi="57" he="48" />为常数项,f<sub>DR</sub>为多普勒调频率,f<sub>DT</sub>为方位向三次项系数;对f<sub>DR</sub>和f<sub>DT</sub>分别作泰勒级数展开得<maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>f</mi><mi>DR</mi></msub><mo>&ap;</mo><msub><mi>f</mi><mi>DRref</mi></msub><mo>+</mo><msub><mi>k</mi><mi>SAC</mi></msub><msub><mi>t</mi><mi>n</mi></msub><mo>+</mo><msub><mi>k</mi><mi>TAC</mi></msub><msubsup><mi>t</mi><mi>n</mi><mn>2</mn></msubsup></mtd></mtr><mtr><mtd><msub><mi>f</mi><mi>DT</mi></msub><mo>&ap;</mo><msub><mi>f</mi><mi>DTref</mi></msub><mo>+</mo><msub><mi>k</mi><mi>SAT</mi></msub><msub><mi>t</mi><mi>n</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>其中f<sub>DRref</sub>为多普勒调频率的参考值,k<sub>SAC</sub>t<sub>n</sub>为多普勒调频率的一阶空变项,<img file="FDA0000483718310000026.GIF" wi="119" he="66" />为多普勒调频率的二阶空变项,f<sub>DTref</sub>为方位三次相位的参考值,k<sub>SAT</sub>t<sub>n</sub>为方位三次相位的一阶空变项;提取为多普勒调频率的一阶空变项k<sub>SAC</sub>t<sub>n</sub>、多普勒调频率的二阶空变项<img file="FDA0000483718310000027.GIF" wi="144" he="67" />为方位三次相位的一阶空变项k<sub>SAT</sub>t<sub>n</sub>;步骤9,对补偿完方位向四阶、五阶相位的距离时域方位频域信号S<sub>9</sub>(t<sub>r</sub>,f<sub>a</sub>)引入用于消除多普勒调频率的一阶空变项k<sub>SAC</sub>t<sub>n</sub>、多普勒调频率的二阶空变项<img file="FDA0000483718310000028.GIF" wi="119" he="67" />和方位三次相位的一阶空变项k<sub>SAT</sub>t<sub>n</sub>的方位高阶非线性变标扰动因子H<sub>ncs</sub>以消除调频率和相位的误差,再对引入方位高阶非线性变标扰动因子H<sub>ncs</sub>后的回波信号进行方位向逆傅里叶变换,得到消除方位空变后的时域信号S<sub>11</sub>(t<sub>r</sub>,t<sub>m</sub>);步骤10,消除方位空变后的时域信号S<sub>11</sub>(t<sub>r</sub>,t<sub>m</sub>)中含有因引入方位高阶非线性调频变标扰动因子H<sub>ncs</sub>产生的高次相位,对消除方位空变后的时域信号S<sub>11</sub>(t<sub>r</sub>,t<sub>m</sub>)补偿剩余高次相位,其方位剩余高次相位补偿因子H<sub>amz</sub>为<maths num="0003"><![CDATA[<math><mrow><msub><mi>H</mi><mi>amz</mi></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mi>j&pi;</mi><mfrac><msub><mi>K</mi><mi>SAC</mi></msub><mn>3</mn></mfrac><msubsup><mi>t</mi><mi>m</mi><mn>3</mn></msubsup><mo>-</mo><mi>j&pi;</mi><mfrac><msub><mi>K</mi><mi>TAC</mi></msub><mn>6</mn></mfrac><msubsup><mi>t</mi><mi>m</mi><mn>4</mn></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>将消除方位空变后的时域信号S<sub>11</sub>(t<sub>r</sub>,t<sub>m</sub>)与方位剩余高次相位补偿因子H<sub>amz</sub>相乘,得到补偿完方位剩余高次相位的时域信号S<sub>12</sub>(t<sub>r</sub>,t<sub>m</sub>);步骤11,对补偿完方位剩余高次相位的时域信号S<sub>12</sub>(t<sub>r</sub>,t<sub>m</sub>)进行方位向的去斜处理,然后作方位向傅里叶变换,得到方位频域聚焦信号S<sub>13</sub>(t<sub>r</sub>,f<sub>a</sub>),从而将图像聚焦在方位频域上。
地址 710071 陕西省西安市太白南路2号