[发明专利]一种斜视SAR地面动目标提取方法

摘要

本发明公开了一种斜视SAR地面动目标提取方法，包括以下步骤S1、计算静止目标和动目标回波信号；S2、进行距离向傅里叶变换；S3、在距离频域、方位时域进行一致距离走动校正；S4、进行距离向残余相位补偿；S5、变参数扩展方位向非线性CS；S6、将经过扩展方位向非线性CS后的回波信号分别乘以对应的方位向去因子，并变换到方位向频域；S7、将两个回波的幅度相减，去除静止目标回波，实现静止目标杂波抑制；S8、幅度检测和动目标提取。本发明解决了单通道斜视模式下地面动目标提取中距离方位耦合严重和动目标运动参数未知的问题，不需要运动参数估计，效率更高，可广泛应用于杂波抑制、动目标提取等领域。

主权项

一种斜视SAR地面动目标提取方法，其特征在于，包括以下步骤：S1、计算成像区域的静止目标和动目标的距离历史，产生静止目标和动目标回波信号；具体实现方法为：设M为一地面动目标，其距离向速度为Vr，沿航迹方向上速度为Vc；P为斜视SAR平台，V为平台飞行速度，θ为平台斜视角，波束中心穿越时刻的瞬时斜距记为R0，方位向时间为t时刻时的瞬时斜距记为R(t)，距离向时间变量记为τ，方位向时间变量记为t；则动目标的距离历史RM(t)表示为：RM(t)=(R0cosθ-Vrt)2+[R0sinθ-(Vt-Vct)]2]]>将上述距离历史对方位向时间t在t＝0处展开到二阶得到：RM(t)≈R0-(Vrcosθ+(V-Vc)sinθ)t+[Vrsinθ+(V-Vc)cosθ]22R0t2]]>在上述动目标距离历史RM(t)中令Vr＝0、Vc＝0，得到静止目标点的距离历史RS(t)表达式为：RS(t)=(R0cosθ)2+(R0sinθ-Vt)2]]>设发射信号为线性调频信号，则经过解调后的动目标回波信号SM(t,τ)表示为：SM(t,τ)=Aωr(τ-2RM(t)c)ωa(tTa)×exp{jπKr(τ-2RM(t)c)2}exp{-j4πfcRM(t)c}]]>经过解调后的静止目标回波信号SS(t,τ)表示为：SS(t,τ)=Aωr(τ-2RS(t)c)ωa(tTa)×exp{jπKr(τ-2RS(t)c)2}exp{-j4πfcRS(t)c}]]>其中，A为回波复幅度，ωr[·]和ωa[·]分别代表距离时间窗和方位时间窗，Kr是发射信号的时间调频斜率，Ta为合成孔径时间，c为光速，fc为载波频率；S2、将步骤S1得到的回波信号进行距离向傅里叶变换，将回波信号变换到距离频域、方位时域；S3、在距离频域、方位时域进行一致距离走动校正，去除静止目标空变的一次距离徙动和动目标的一阶距离徙动；具体实现方法为：将方位向时间变量t进行变标，变标后的方位向时间变量记为tm，其变标算法为：t=fcfτ+fctm]]>其中，fτ为距离向频率变量，tm为变换后的方位向时间变量；采用上述变标算法去除静止目标空变的一次距离徙动和动目标的一阶距离徙动，得到变换后的回波矩阵S(tm,fτ)；S4、进行距离向残余相位补偿，得到距离徙动校正和距离向压缩后的回波矩阵；具体实现方法为：将回波矩阵S(tm，fτ)乘以距离向残余相位补偿因子ψRcom(tm,fτ)，得到距离徙动校正和距离向压缩后的回波矩阵：S2(tm,fτ)=S(tm,fτ)·ψRcom(tm,fτ)=S(tm,fτ)·exp{jπfτ2Kr-j2π(V2sin2θR0c)tm2fτ};]]>S5、变参数扩展方位向非线性CS，将静止目标的多普勒质心和调频斜率分别进行均衡；包括以下子步骤：S51、将同一距离单元的静止场景目标点的多普勒质心和多普勒调频斜率随方位向的空变分别建模为一次函数和二次函数：fDc＝fDcref+atmfDr=fDrref+btm+dtm2]]>其中，fDcref为静止场景参考点的多普勒质心，fDrref为静止场景参考点的多普勒调频斜率，a、b和d分别为拟合系数；S52、根据上述拟合系数，构造四次滤波和方位向非线性CS函数：四次滤波函数为：H1F(fa)=exp{jπ(G3fa3+G4fa4)}]]>方位向非线性CS函数为：HANCS(tm)=exp{jπ(e2tm2+e3tm3+e4tm4)}]]>其中，G3，G4，e2，e3，e4分别为四次滤波和非线性CS系数：e2=-2aα+(2α-1)fDrrefe3=2b(e2+a)(e2+fDrref)-fD3(a+e2)2-e2[b(2e2+a+fDrref)-fD3(a+e2)]3(fDrref-a)2e4=Q/4-(fDrref-a)fDrref3e2G4a-fDrrefG3=b(2e2+a+fDrref)-fD3(a+e2)3(fDrref-a)2e2fDrrefG4=P/6-Q(a-fDrref)/4(fDrref-a)2e2fDrref2(e2+fDrref)---(1)]]>其中P=-[d(e2+fDrref)2-b2(e2+fDrref)]-3fD3b(a+e2)+3G3e2bfDrref(3fDrrefe2-2ae2+fDrrefa)-3e3b(e2-2fDrref+3a)]]> Q=-3fD3b+3G3e2bfDrref2-3e3b]]>其中，fD3为多普勒参数三阶项，α为可变参数；S53、在式(1)中，令参数α分别等于α1和α2，得到两组系数，并分别构造两组四次滤波和非线性CS函数，第一组系数表示为：H1F,α1(fa)=exp{jπ(G3,α1fa3+G4,α1fa4)}]]>HANCS,α1(tm)=exp{jπ(e2,α1tm3+e4,α1tm4+e4,α1tm4)}]]>其中为第一组四次滤波和非线性CS系数，将α＝α1带入(1)即可得出；第二组系数表示为：H1F,α2(fa)=exp{jπ(G3,α2fa3+G4,α2fa4)}]]>HANCS,α2(tm)=exp{jπ(e2,α2tm2+e3,α2tm3+e4,α2tm4)}]]>其中为第二组四次滤波和非线性CS系数，将α＝α2带入(1)即可得出；S54、将回波信号S2(tm,fτ)通过距离向IFFT和方位向FFT变换到距离时域方位频域，得到S2(fa,τ)；S55、将回波信号S2(fa,τ)分别与和相乘，然后经过方位向IFFT变换到二维时域，得到和S56、将和分别与和相乘，得到变参数扩展方位向非线性CS后的二维时域回波信号和S6、变参数去斜处理：将经过扩展方位向非线性CS后的回波信号分别乘以对应的方位向去斜因子，并变换到方位向频域；S7、进行静止场景杂波抑制：将两个回波的幅度相减，去除静止目标回波，实现静止目标杂波抑制；S8、幅度检测和动目标提取：采用恒虚警处理，进行幅度检测，实现动目标提取。