1.一种基于Omega-K算法的正交频分复用信号的合成孔径成像方法,其特征在于,包括以下步骤:步骤1、对基于OFDM信号的SAR成像的回波数据s r (τ,η)去掉中心载频后,进行距离向傅里叶变换;对距离向参考信号s 0 (τ)进行傅里叶变换;两个傅里叶变换后信号进行共轭相乘即可得到距离压缩频域信号S ic (f τ ,η);步骤2、对信号S ic (f τ ,η)进行方位向傅里叶变换,得到二维频域信号S ic (f τ ,f η );步骤3、二维频域信号S ic (f τ ,f η )与二维频域滤波器H(f τ ,f η )相乘,得到未校正的距离向方位向高阶耦合相位的二维频域信号S b (f τ ,f η );步骤4、利用Stolt插值校正距离向与方位向之间的高阶耦合相位,得到校正后距离向方位向高阶耦合相位的二维频域信号S l (f τ ′,f η );步骤5、对信号S l (f τ ′,f η )进行二维傅里叶逆变换,得到时域数据s out (τ′,η)作为成像结果。
2.根据权利要求1所述的基于Omega-K算法的正交频分复用信号的合成孔径成像方法,其特征在于,步骤1的具体过程为:距离向雷达参考信号s 0 (τ)表示为:其中,τ为快时间变量,w r 为距离窗函数,R 0 为雷达距离目标的最短距离,c为光速;N为OFDM信号子载波数目;0≤k≤N-1,S k 为第k个子载波所传输的信息; T为雷达脉冲宽度,j为虚部符号, 基于OFDM信号的SAR成像的回波数据s r (τ,η)表示为:其中, v为雷达平台的运动速度,η为慢时间变量;w a 为方位包络窗函数,η c 为波束中心偏离时间;g m 为目标反射系数,f c 为载波频率, 去掉中心载频后,s r (τ,η)表示为:分别对s r (τ,η)和s 0 (τ)进行距离向傅里叶变换并共轭相乘,得到距离压缩后的频域信号S ic (f τ ,η):其中,f τ 为快时间频域变量,sinc为辛格函数,0≤n≤N-1,S n 为第n个子载波所传输的信息,W r 为距离窗函数的频域表达形式。
3.根据权利要求2所述的基于Omega-K算法的正交频分复用信号的合成孔径成像方法,其特征在于,步骤2中,通过驻定相位原理得到二维频域信号S ic (f τ ,f η ),具体为:上式中,S ic (f τ ,f η )的相位为 根据驻定相位原理 则S ic (f τ ,f η )表示为:其中,f η 为慢时间频域变量,Wa 为 方位包络窗函数的频域表达形式; 为多普勒频率中心,t 0 为零多普勒时刻。
4.根据权利要求3所述的基于Omega-K算法的正交频分复用信号的合成孔径成像方法,其特征在于,步骤3中,二维频域滤波器H(f τ ,f η )为:其中,R ref 为参考距离;方位向脉冲压缩后的二维频域信号S b (f τ ,f η )为:
5.根据权利要求4所述的基于Omega-K算法的正交频分复用信号的合成孔径成像方法,其特征在于,步骤4的具体过程为:通过Stolt插值映射方式实现如下校正: 对S b (f τ ,f η )进行上述操作得到S l (f τ ′,f η )的表达式为:其中,f τ ′为Stolt插值后的快时间频域变量。
6.根据权利要求5所述的基于Omega-K算法的正交频分复用信号的合成孔径成像方法,其特征在于,步骤5的具体过程为:在S l (f τ ′,f η )中将 记为常数,S k 为1与-1编码,则时域数据s out (τ′,η)的表达式为:其中,τ′为Stolt插值后的快时间变量。