基于RD算法的时频域混合校正距离徙动的SAR成像方法与流程

本发明涉及时频域混合校正距离徙动的SAR成像方法。

背景技术：

在机载SAR对场景成像过程中，由于合成孔径雷达为宽频带发射信号，近场工作的大阵列，点目标回波在距离向与方位向必然会存在耦合现象，具体表现为系统响应在快时间域与慢时间域(坐标)二维平面呈现为曲线，即距离徙动。现考察最为一般的机载SAR斜视模型，该模型下波束射线指向的斜视角为θ₀，雷达合成孔径中心记为A点，它在X轴的位置为X₀(＝R_Btanθ₀)，RB为雷达航线与目标场景中心线的垂直距离，距点目标P的距离为R₀＝(R_Bsecθ₀)，有效合成孔径长度为L。为了得到斜视模型下雷达到目标距离R的近似式，以天线波束中心线指向目标时雷达的位置(X＝X₀)为基准，根据斜视SAR模型的几何关系，可以得到载机沿航线飞行时位于某一坐标X(记为A’点)的瞬时斜距：

由于L＜＜R₀，对式(2-1)在X＝X₀附近作泰勒级数展开，省略(x-x₀)的三次项以上的高次项，可以得到式(2-1)的近似式：

即有

式(2-3)表明，在合成孔径某处(X)到目标的距离与阵列中心A(X＝X₀)到该目标的距离差有(X＝X₀)的一次项和二次项所组成，通常把一次项(即线性项)称为距离走动，将二次项称为距离弯曲。

传统的RD算法对于距离徙动有时域校正和频域校正两种方式。其中，时域校正算法利用系统匹配函数分别对成像场景的距离向和方位向进行匹配滤波，距离徙动的校正是在距离向脉冲压缩时乘以线性相位因子来完成的，这种方法可以很好地校正雷达成像中的线性距离走动，但是不适用于雷达斜视角较大的情况，斜视角较大时距离走动校正不明显，成像目标难以识别。频域校正算法将场景回波信号进行二维傅里叶变换，在频域通过乘以频率函数进行二维去耦和脉冲压缩，再变换到距离-多普勒域对方位向匹配滤波，最后进行方位逆傅里叶变换完成方位压缩，这样就完成了整个频域校正距离徙动的成像算法。但是这种算法仅适用于雷达斜视角较小的情况，并且对距离徙动，尤其是线性距离走动项的校正效果欠佳，导致方位向聚焦效果较差，甚至散焦。

技术实现要素：

本发明的目的是为了解决现有时域校正算法不适用于雷达斜视角较大的情况，斜视角较大时距离走动校正不明显，成像目标难以识别以及频域校正算法仅适用于雷达斜视角较小的情况，并且对距离徙动，尤其是线性距离走动项的校正效果欠佳，导致方位向聚焦效果较差，甚至散焦的问题，而提出基于RD算法的时频域混合校正距离徙动的SAR成像方法。

基于RD算法的时频域混合校正距离徙动的SAR成像方法具体过程为：

步骤一、对雷达数据进行采样，得到原始数据；

步骤二、对原始数据进行距离FFT，将原始数据由时域变换到频域f_r-t_m，得到S₁(f_r,t_m,R₀)；

所述FFT为快速傅立叶变换；为原始数据距离向时域时间；f_r为原始数据距离向频域频率；t_m为原始数据方位向时域时间，R₀为雷达到目标的距离；S₁(f_r,t_m,R₀)为原始数据距离向傅里叶变换后的数据；

步骤三、将S₁(f_r,t_m,R₀)乘以H₁以完成f_r-t_m各点在距离向上的距离徙动补偿一致，得到S₂(f_r,t_m,R₀)；

H₁为线性相位函数；S₂(f_r,t_m,R₀)为S₁(f_r,t_m,R₀)补偿一致后的数据；

步骤四、对步骤三完成的补偿一致后的数据S₂(f_r,t_m,R₀)进行方位FFT，将补偿一致后的数据由时域f_r-t_m变换到二维频域f_r-f_a，得到S₃(f_r,f_a,R₀)；

f_a为原始数据方位向频域频率，S₃(f_r,f_a,R₀)为S₂(f_r,t_m,R₀)方位向傅里叶变换后的数据；

步骤五、将步骤四变换到二维频域f_r-f_a后的数据S₃(f_r,f_a,R₀)乘以H₂在频域上完成距离徙动校正以及距离压缩并进行距离向逆傅里叶变换，得到

H₂为距离徙动校正以及脉冲压缩函数，为S₃(f_r,f_a,R₀)距离徙动校正以及脉冲压缩并进行距离向逆傅里叶变换后的数据；

步骤六、将乘以方位脉压匹配函数H₃并进行方位IFFT，得到SAR图像；

IFFT为逆傅立叶变换，SAR为合成孔径雷达。

本发明的有益效果为：

在相同的雷达参数模型下，如本发明给出的参数模型，传统时域校正方法在斜视角为10°以内校正效果尚佳，斜视角大于15°后距离徙动校正不明显，导致方向压缩效果差，目标混叠难以识别；传统频域校正方法在斜视角8°以内可以校正距离徙动，斜视角大于10°时徙动校正难以实现，导致方位压缩出现散焦，目标难以识别，而时频域混合校正方法可在雷达斜视角60°以内对距离徙动进行良好地校正，保证方位向的聚焦效果，对场景中的目标清晰成像。综上，机载SAR时频域混合校正距离徙动算法对距离徙动的校正、方位聚焦效果以及场景成像效果都要比传统的两种RD算法好。

附图说明

图1为本发明时频域混合校正距离徙动流程图；

图2为传统时域校正方法距离压缩图；

图3为传统时域校正方法方位压缩图；

图4为传统时域校正方法三维成像图；

图5为传统频域校正方法距离压缩图；

图6为传统频域校正方法方位压缩图；

图7为传统频域校正方法三维成像图；

图8为本发明时频域混合校正方法距离压缩图；

图9为本发明时频域混合校正方法方位压缩图；

图10为本发明时频域混合校正方法三维成像图。

具体实施方式

具体实施方式一：结合图1说明本实施方式，本实施方式的基于RD算法的时频域混合校正距离徙动的SAR成像方法具体过程为：

步骤一、对雷达数据进行采样，得到原始数据；采用MATLAB对雷达数据进行采样，得到原始数据；该模型为机载SAR，雷达架设在飞机上；

步骤二、对原始数据进行距离FFT，将原始数据由时域变换到频域f_r-t_m，得到S₁(f_r,t_m,R₀)；

所述FFT为快速傅立叶变换；为原始数据距离向时域时间；f_r为原始数据距离向频域频率；t_m为原始数据方位向时域时间，R0为雷达到目标的距离；S₁(f_r,t_m,R₀)为原始数据距离向傅里叶变换后的数据；

步骤三、将S₁(f_r,t_m,R₀)乘以H₁以完成f_r-t_m各点在距离向上的距离徙动补偿一致，得到S₂(f_r,t_m,R₀)；

H₁为线性相位函数；S₂(f_r,t_m,R₀)为S₁(f_r,t_m,R₀)补偿一致后的数据；

步骤四、对步骤三完成的补偿一致后的数据S₂(f_r,t_m,R₀)进行方位FFT，将补偿一致后的数据由时域f_r-t_m变换到二维频域f_r-f_a，得到S₃(f_r,f_a,R₀)；

f_a为原始数据方位向频域频率，S₃(f_r,f_a,R₀)为S₂(f_r,t_m,R₀)方位向傅里叶变换后的数据；

步骤五、将步骤四变换到二维频域f_r-f_a后的数据S₃(f_r,f_a,R₀)乘以H₂在频域上完成距离徙动校正以及距离压缩并进行距离向逆傅里叶变换，得到

H₂为距离徙动校正以及脉冲压缩函数，为S₃(f_r,f_a,R₀)距离徙动校正以及脉冲压缩并进行距离向逆傅里叶变换后的数据；

步骤六、将乘以方位脉压匹配函数H₃并进行方位IFFT，得到SAR图像；

IFFT为逆傅立叶变换，SAR为合成孔径雷达。

具体实施方式二：本实施方式与具体实施方式一不同的是：所述步骤三中将S₁(f_r,t_m,R₀)乘以H₁以完成f_r-t_m各点在距离向上的距离徙动补偿一致，得到S₂(f_r,t_m,R₀)；

H₁为线性相位函数；S₂(f_r,t_m,R₀)为S₁(f_r,t_m,R₀)补偿一致后的数据；

具体过程为：

将原始数据从域变换到f_r-t_m域得S₁(f_r,t_m,R₀)，然后将S₁(f_r,t_m,R₀)乘以线性相位函数H₁，得到S₂(f_r,t_m,R₀)

式中，△R(t_m)＝-(Vsinθ₀)t_m，

j为虚数，j²＝-1；△R(t_m)为雷达到目标的距离差，△R(t_m)为随t_m变化的变量；C为光速3×10⁸m/s；f_c为雷达信号的频率；V为雷达速度；θ₀为雷达斜视角，0°≤θ₀≤90°。

其它步骤及参数与具体实施方式一相同。

具体实施方式三：本实施方式与具体实施方式一或二不同的是：所述步骤五中将步骤四变换到二维频域f_r-f_a后的数据S₃(f_r,f_a,R₀)乘以H₂在频域上完成距离徙动校正以及距离压缩并进行逆傅里叶变换，得到

H₂为距离徙动校正以及脉冲压缩函数，为S₃(f_r,f_a,R₀)距离徙动校正以及脉冲压缩并进行逆傅里叶变换后的数据；

具体过程为：

将S₃(f_r,f_a,R₀)乘以H₂在频域上完成距离徙动校正以及距离压缩并进行距离向逆傅里叶变换，得到

γ是线性调频率，为雷达信号频带宽度与信号脉冲宽度的比值；λ是雷达信号波长。

其它步骤及参数与具体实施方式一或二相同。

具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：所述步骤六中将乘以方位脉压匹配函数H₃并进行方位IFFT，得到SAR图像；

IFFT为逆傅立叶变换，SAR为合成孔径雷达；

具体过程为：

将乘以方位脉压匹配函数H₃并进行方位IFFT，得到SAR图像，方位脉压匹配函数H₃为

f_aM为方位向多普勒频率的最大值；

其它步骤及参数与具体实施方式一至三之一相同。

采用以下实施例验证本发明的有益效果：

实施例一：

本实施例基于RD算法的时频域混合校正距离徙动的SAR成像方法具体是按照以下步骤制备的：

现提供一种参数模型：

表7-1时频域混合校正距离徙动算法的基本参数值

点目标的仿真结果如图2、3、4、5、6、7、8、9、10；图2为传统时域校正方法距离压缩图；图3为传统时域校正方法方位压缩图；图4为传统时域校正方法三维成像图；图5为传统频域校正方法距离压缩图；图6为传统频域校正方法方位压缩图；图7为传统频域校正方法三维成像图；图8为本发明时频域混合校正方法距离压缩图；图9为本发明时频域混合校正方法方位压缩图；图10为本发明时频域混合校正方法三维成像图。从图8和图2，图5的对比结果可明显得出本发明方法距离徙动校正效果优于传统方法。从图9和图3，图6的对比结果可明显得出本发明方法方位聚焦效果优于传统方法。从图10和图4，图7的对比结果可明显得出本发明方法场景成像效果要比传统方法好。

本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。