CN104730496B - 一种斜视sar地面动目标提取方法 - Google Patents

Abstract

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

Description

一种斜视SAR地面动目标提取方法

技术领域

本发明属于雷达技术领域，它特别涉及合成孔径雷达(SAR)成像技术中的斜视SAR地面动目标提取方法。

背景技术

斜视合成孔径雷达(SAR)可以根据实际应用需求改变波束指向，具有很高的灵活性和机动性，可应用于指定地区的持续观测与侦查，并实现热点地区的重复观测。此外，工作在斜视模式下的雷达平台与目标场景距离更远，增加了平台的隐蔽性和安全性。使用斜视SAR进行地面动目标提取，可以提前发现动目标并得到其位置信息以及运动参数信息。因此，研究斜视SAR地面动目标提取具有重要的实际意义。

然而，斜视SAR地面动目标提取中存在两个难题。首先，动目标位置和运动参数是未知的，这为运动目标回波的距离单元徙动矫正(RCMC)和方位向压缩带来了困难。使用静止场景的RCMC和方位向压缩参数对动目标进行处理会导致动目标在距离向和方位向散焦。另外，由于雷达工作在斜视模式，导致动目标回波的距离向和方位向耦合更加严重，并带来回波的距离单元徙动和多普勒参数的二维空变性。传统的动目标提取方法需要对目标的运动参数进行估计，以进行距离单元徙动矫正。然而，这种方法需要进行二维参数搜索，运算量大。另外，当动目标回波的信噪比较低或回波中存在强的静止目标散射点时，估计结果的误差很大，造成检测结果不准确。在文献“A generalization of DPCA processing formultichannel SAR/GMTI radars”(《IEEE Transactions on Geoscience and RemoteSensing》vol.51,no.1,pp.560-572,2013)中，提出通过多通道DPCA的方法检测动目标。但是，这种方法需要多个接收通道，增加了系统成本和信号处理复杂度。在文献“Groundmoving targets indication in SAR images by symmetric defocusing”(《IEEEGeoscience and Remote Sensing Letters》vol.10,no.2,pp.241–245,2013.)中提出了一种动目标检测方法。该方法首先使用两个不同的滤波器使图像散焦，通过比较图像的散焦程度来检测动目标。但是，通常情况下，无法确定对原始图像散焦处理到何种程度；而且，弱目标会被强散射目标的回波所覆盖。

发明内容

本发明的目的在于克服单通道斜视模式下地面动目标提取中距离方位耦合严重和动目标运动参数未知的问题，提供了一种不需要运动参数估计，效率更高，可广泛应用于杂波抑制、动目标提取等领域的斜视SAR地面动目标提取方法。

本发明的目的是通过以下技术方案来实现的：一种斜视SAR地面动目标提取方法，包括以下步骤：

S1、计算成像区域的静止目标和动目标的距离历史，产生静止目标和动目标回波信号；

S2、将步骤S1得到的回波信号进行距离向傅里叶变换，将回波信号变换到距离频域、方位时域；

S3、在距离频域、方位时域进行一致距离走动校正，去除静止目标空变的一次距离徙动和动目标的一阶距离徙动；

S4、进行距离向残余相位补偿，得到距离徙动校正和距离向压缩后的回波矩阵；

S5、变参数扩展方位向非线性CS，将静止目标的多普勒质心和调频斜率分别进行均衡；

S6、变参数去斜处理：将经过扩展方位向非线性CS后的回波信号分别乘以对应的方位向去因子，并变换到方位向频域；

S7、进行静止场景杂波抑制：将两个回波的幅度相减，去除静止目标回波，实现静止目标杂波抑制；

S8、幅度检测和动目标提取：采用恒虚警处理，进行幅度检测，实现动目标提取。

进一步地，所述的步骤S1具体实现方法为：设M为一地面动目标，其距离向速度为V_r，沿航迹方向上速度为V_c；P为斜视SAR平台，V为平台飞行速度，θ为平台斜视角，波束中心穿越时刻的瞬时斜距记为R₀，方位向时间为t时刻时的瞬时斜距记为R(t)，距离向时间变量记为τ，方位向时间变量记为t；

则动目标的距离历史R_M(t)表示为：

将上述距离历史对方位向时间t在t＝0处展开到二阶得到：

在上述动目标距离历史R_M(t)中令V_r＝0、V_c＝0，得到静止目标点的距离历史R_S(t)表达式为：

设发射信号为线性调频信号，则经过解调后的动目标回波信号S_M(t,τ)表示为_：

经过解调后的静止目标回波信号S_S(t,τ)表示为：

其中，A为回波复幅度，ω_r[·]和ω_a[·]分别代表距离时间窗和方位时间窗，K_r是发射信号的时间调频斜率，T_a为合成孔径时间，c为光速，f_c为载波频率。

进一步地，所述的步骤S3的具体实现方法为：将方位向时间变量t进行变标，变标后的方位向时间变量记为t_m，其变标算法为：

其中，f_τ为距离向频率变量，t_m为变换后的方位向时间变量；

采用上述变标算法去除静止目标空变的一次距离徙动和动目标的一阶距离徙动，得到变换后的回波矩阵S(t_m,f_τ)。

进一步地，所述的步骤S4的具体实现方法为：将回波矩阵S(t_m,f_τ)乘以距离向残余相位补偿因子ψ_Rcom(t_m,f_τ)，得到距离徙动校正和距离向压缩后的回波矩阵：

进一步地，所述的步骤S5包括以下子步骤：

S51、将同一距离单元的静止场景目标点的多普勒质心和多普勒调频斜率随方位向的空变分别建模为一次函数和二次函数：

f_Dc＝f_Dcref+at_m

其中，f_Dcref为静止场景参考点的多普勒质心，f_Drref为静止场景参考点的多普勒调频斜率，a、b和d分别为拟合系数；

S52、根据上述拟合系数，构造四次滤波和方位向非线性CS函数：

四次滤波函数为：

方位向非线性CS函数为：

其中，G₃，G₄，e₂，e₃，e₄分别为四次滤波和非线性CS系数：

其中

P＝-[d(e₂+f_Drref)²-b²(e₂+f_Drref)]-3f_D3b(a+e₂)

+3G₃e₂bf_Drref(3f_Drrefe₂-2ae₂+f_Drrefa)-3e₃b(e₂-2f_Drref+3^a)

其中，f_D3为多普勒参数三阶项，α为可变参数；

S53、在式(1)中，令参数α分别等于α₁和α₂，得到两组系数，并分别构造两组四次滤波和非线性CS函数，第一组系数表示为：

其中为第一组四次滤波和非线性CS系数，将α＝α₁带入(1)即可得出；

第二组系数表示为：

其中为第二组四次滤波和非线性CS系数，将α＝α₂带入(1)即可得出；

S54、将回波信号S₂(t_m,f_τ)通过距离向IFFT和方位向FFT变换到距离时域方位频域，得到S₂(f_a,τ)；

S55、将回波信号S₂(f_a,τ)分别与和相乘，然后经过方位向IFFT变换到二维时域，得到和

S56、将和分别与和相乘，得到变参数扩展方位向非线性CS后的二维时域回波信号和

进一步地，所述的步骤S6具体实现方法为：将回波信号和分别乘以对应的方位向去斜因子，并变换到方位向频域，得到方位向聚焦的频域回波信号和

其中，FFT_azi{·}为方位向傅里叶变换，和分别为去斜因子：

进一步地，所述的步骤S7具体实现方法为：

进一步地，所述的步骤S8具体实现方法为：采用恒虚警处理，取一个自适应门限值δ，保持输出的虚警率不超过恒定值，进行幅度检测可以实现动目标提取，对于距离门τ₀，若满足：

|S₅(f_a,τ₀)|≥δ

则该距离门存在动目标，则对该动目标进行幅度检测，得到动目标提取结果。

本发明的有益效果是：去除了静止目标空变的一阶距离徙动和未知运动参数的动目标一阶距离徙动；并使用变参数扩展方位向非线性CS和相应的去斜处理进行方位向聚焦，并根据两次聚焦后回波的动目标差异进行幅度检测，从而解决了单通道斜视模式下地面动目标提取中距离方位耦合严重和动目标运动参数未知的问题，与现有SAR动目标方法相比，不需要运动参数估计，效率高、成本低，可广泛应用于杂波抑制、动目标提取等领域。

附图说明

图1为本发明的目标提取方法的流程图；

图2为大斜视SAR的几何结构图；

图3为本发明具体实施例中得到的静止目标和动目标方位向频域成像结果，其中，(a)为参数为α₁的方位向扩展非线性CS的成像结果，(b)是参数为α₂的方位向扩展非线性CS的成像结果；

图4为本发明具体实施例中得的去除静止杂波后的动目标提取结果；

图5为图4中三个动目标所在的距离单元的方位向切片。

具体实施方式

本发明主要采用仿真实验的方法进行验证，所有步骤、结论都在Matlab2012上验证正确性。下面就具体实施例和附图对本发明的技术方案作进一步的详细描述。

如图1所示，一种斜视SAR地面动目标提取方法，首先去除动目标和静止场景的一阶RCM。然后在距离频域方位时域进行二阶距离徙动校正并完成距离压缩。接着，采用变参数扩展方位向非线性CS算法将静止目标点的多普勒质心和调频斜率均衡，再使用对应的去斜因子进行方位向聚焦；最后在方位频域距离时域对两幅聚焦后的图像进行幅度相减，去除静止场景回波，通过幅度检测可得到动目标提取结果。具体包括以下步骤：

S1、计算成像区域的静止目标和动目标的距离历史，产生静止目标和动目标回波信号，其具体实现方法为：设M为一地面动目标，其距离向速度为V_r，沿航迹方向上速度为V_c；P为斜视SAR平台，V为平台飞行速度，θ为平台斜视角，波束中心穿越时刻的瞬时斜距记为R₀，方位向时间为t时刻时的瞬时斜距记为R(t)，距离向时间变量记为τ，方位向时间变量记为t，如图2所示；

则动目标的距离历史R_M(t)表示为：

将上述距离历史对方位向时间t在t＝0处展开到二阶得到：

在上述动目标距离历史R_M(t)中令V_r＝0、V_c＝0，得到静止目标点的距离历史R_S(t)表达式为：

设发射信号为线性调频信号，则经过解调后的动目标回波信号S_M(t,τ)表示为：

经过解调后的静止目标回波信号S_S(t,τ)表示为：

其中，A为回波复幅度，ω_r[·]和ω_a[·]分别代表距离时间窗和方位时间窗，K_r是发射信号的时间调频斜率，T_a为合成孔径时间，c为光速，f_c为载波频率。

本实施例的仿真所需的参数如表一所示。

表一

S2、将步骤S1得到的回波信号进行距离向傅里叶变换，将回波信号变换到距离频域、方位时域，得到矩阵记为S(t,f_τ)：S(t,f_τ)＝FFT_range{S(t,τ)}

S3、在距离频域、方位时域进行一致距离走动校正，去除静止目标空变的一次距离徙动和动目标的一阶距离徙动，其具体实现方法为：将方位向时间变量t进行变标，变标后的方位向时间变量记为t_m，其变标算法为：

其中，f_τ为距离向频率变量，t_m为变换后的方位向时间变量；

采用上述变标算法去除静止目标空变的一次距离徙动和动目标的一阶距离徙动，得到变换后的回波矩阵S(t_m,f_τ)。

S4、进行距离向残余相位补偿，得到距离徙动校正和距离向压缩后的回波矩阵，其具体实现方法为：将回波矩阵S(t_m,f_τ)乘以距离向残余相位补偿因子ψ_Rcom(t_m,f_τ)，得到距离徙动校正和距离向压缩后的回波矩阵：

距离向残余相位补偿后，回波信号中静止场景的残余距离徙动已经去除，且静止目标和动目标距离向得到很好的聚焦。

S5、经过上述步骤，静止场景和动目标点目标在距离向得到了很好的聚焦，然而，静止场景的多普勒质心和调频斜率沿方位向是空变的，再变参数扩展方位向非线性CS，将静止目标的多普勒质心和调频斜率分别进行均衡，包括以下子步骤：

S51、将同一距离单元的静止场景目标点的多普勒质心和多普勒调频斜率随方位向的空变分别建模为一次函数和二次函数：

f_Dc＝f_Dcref+at_m

其中，f_Dcref为静止场景参考点的多普勒质心，f_Drref为静止场景参考点的多普勒调频斜率，a、b和d分别为拟合系数；

S52、根据上述拟合系数，构造四次滤波和方位向非线性CS函数：

四次滤波函数为：

方位向非线性CS函数为：

其中，G₃，G₄，e₂，e₃，e₄分别为四次滤波和非线性CS系数：

其中

P＝-[d(e₂+f_Drref)²-b²(e₂+f_Drref)]-3f_D3b(a+e₂)

+3G₃e₂bf_Drref(3f_Drrefe₂-2ae₂+f_Drrefa)-3e₃b(e₂-2f_Drref+3a)

其中，f_D3为多普勒参数三阶项，α为可变参数，通过选择不同的α值，可以使动目标与静止目标区分开；

S53、在式(1)中，令参数α分别等于α₁和α₂，得到两组系数，并分别构造两组四次滤波和非线性CS函数，第一组系数表示为：

其中为第一组四次滤波和非线性CS系数，将α＝α₁带入(1)即可得出；

第二组系数表示为：

其中为第二组四次滤波和非线性CS系数，将α＝α₂带入(1)即可得出；

S54、将回波信号S₂(t_m,f_τ)通过距离向IFFT和方位向FFT变换到距离时域方位频域，得到S₂(f_a,τ)；

S55、将回波信号S₂(f_a,τ)分别与和相乘，然后经过方位向IFFT变换到二维时域，得到和

S56、将和分别与和相乘，得到变参数扩展方位向非线性CS后的二维时域回波信号和

S6、变参数去斜处理：将经过扩展方位向非线性CS后的回波信号分别乘以对应的方位向去因子，并变换到方位向频域，具体实现方法为：将回波信号和分别乘以对应的方位向去斜因子，并变换到方位向频域，得到方位向聚焦的频域回波信号和

其中，FFT_azi{·}为方位向傅里叶变换，和分别为去斜因子：

图3为本实施例中得到的静止目标和动目标方位向频域聚焦成像结果，其中，(a)为参数为α₁的方位向扩展非线性CS的成像结果，(b)是参数为α₂的方位向扩展非线性CS的成像结果。图中圈出的目标点M1，M2和M3是三个动目标，其距离向速度分别为10m/s和3m/s，-7m/s，切距离向速度为7m/s和3.4m/s,-5m/s。

S7、进行静止场景杂波抑制：经过方位向去斜处理后，在方位向频域，回波和中的静止场景目标点得到了聚焦，且聚焦位置相同；而动目标是散焦的，而且两个回波中的动目标的位置和散焦程度不同。因此，将两个回波的幅度相减，去除静止目标回波，实现静止目标杂波抑制：

S8、幅度检测和动目标提取：采用恒虚警处理，进行幅度检测，实现动目标提取：取一个自适应门限值δ，保持输出的虚警率不超过恒定值，进行幅度检测可以实现动目标提取，对于距离门τ₀，若满足：

|S₅(f_a,τ₀)|≥δ

则该距离门存在动目标，则对该动目标进行幅度检测，得到动目标提取结果，如图4、图5所示。从图中可以看出，本发明提供的方法可以很好的实现斜视SAR动目标提取。

本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

Claims (4)

1.一种斜视SAR地面动目标提取方法，其特征在于，包括以下步骤：

S1、计算成像区域的静止目标和动目标的距离历史，产生静止目标和动目标回波信号；具体实现方法为：设M为一地面动目标，其距离向速度为V_r，沿航迹方向上速度为V_c；P为斜视SAR平台，V为平台飞行速度，θ为平台斜视角，波束中心穿越时刻的瞬时斜距记为R₀，方位向时间为t时刻时的瞬时斜距记为R(t)，距离向时间变量记为τ，方位向时间变量记为t；

则动目标的距离历史R_M(t)表示为：

$R_M\left(t\right)=\sqrt{{(R_{0}cos\mathrm{\θ}-V_{r}t)}^2+{\[R_0\sin \mathrm{\θ}-(Vt-V_{c}t)\]}^2}$

将上述距离历史对方位向时间t在t＝0处展开到二阶得到：

$R_M\left(t\right)\≈R_0-(V_{r}cos\mathrm{\θ}+(V-V_c\left)sin\mathrm{\θ}\right)t+\frac{{\[V_{r}sin\mathrm{\θ}+(V-V_c)cos\mathrm{\θ}\]}^2}{2R_0}{t}^2$

在上述动目标距离历史R_M(t)中令V_r＝0、V_c＝0，得到静止目标点的距离历史R_S(t)表达式为：

$R_S\left(t\right)=\sqrt{{\left(R_{0}cos\mathrm{\θ}\right)}^2+{(R_{0}sin\mathrm{\θ}-Vt)}^2}$

设发射信号为线性调频信号，则经过解调后的动目标回波信号S_M(t,τ)表示为：

$\begin{array}{c}{S}_M(t,\mathrm{\τ})={\mathrm{A\ω}}_r(\mathrm{\τ}-\frac{2R_M\left(t\right)}{c}){\mathrm{\ω}}_a\left(\frac{t}{T_a}\right)\\ \×\exp \left\{{\mathrm{j\πK}}_r{(\mathrm{\τ}-\frac{2R_M\left(t\right)}{c})}^2\right\}\exp \{-j4{\mathrm{\πf}}_c\frac{R_M\left(t\right)}{c}\}\end{array}$

经过解调后的静止目标回波信号S_S(t,τ)表示为：

$\begin{array}{c}{S}_S(t,\mathrm{\τ})={\mathrm{A\ω}}_r(\mathrm{\τ}-\frac{2R_S\left(t\right)}{c}){\mathrm{\ω}}_a\left(\frac{t}{T_a}\right)\\ \×\exp \left\{{\mathrm{j\πK}}_r{(\mathrm{\τ}-\frac{2R_S\left(t\right)}{c})}^2\right\}\exp \{-j4{\mathrm{\πf}}_c\frac{R_S\left(t\right)}{c}\}\end{array}$

其中，A为回波复幅度，ω_r[·]和ω_a[·]分别代表距离时间窗和方位时间窗，K_r是发射信号的时间调频斜率，T_a为合成孔径时间，c为光速，f_c为载波频率；

S2、将步骤S1得到的回波信号进行距离向傅里叶变换，将回波信号变换到距离频域、方位时域；

S3、在距离频域、方位时域进行一致距离走动校正，去除静止目标空变的一次距离徙动和动目标的一阶距离徙动；具体实现方法为：将方位向时间变量t进行变标，变标后的方位向时间变量记为t_m，其变标算法为：

$t=\frac{f_c}{f_{\mathrm{\τ}}+f_c}{t}_m$

其中，f_τ为距离向频率变量，t_m为变换后的方位向时间变量；

采用上述变标算法去除静止目标空变的一次距离徙动和动目标的一阶距离徙动，得到变换后的回波矩阵S(t_m,f_τ)；

S4、进行距离向残余相位补偿，得到距离徙动校正和距离向压缩后的回波矩阵；具体实现方法为：将回波矩阵S(t_m，f_τ)乘以距离向残余相位补偿因子ψ_Rcom(t_m,f_τ)，得到距离徙动校正和距离向压缩后的回波矩阵：

$\begin{array}{c}{S}_2(t_m,f_{\mathrm{\τ}})=S(t_m,f_{\mathrm{\τ}})\·{\mathrm{\ψ}}_{Rcom}(t_m,f_{\mathrm{\τ}})\\ =S(t_m,f_{\mathrm{\τ}})\·\exp \{j\mathrm{\π}\frac{f_{\mathrm{\τ}}^2}{K_r}-j2\mathrm{\π}\left(\frac{V^2{\sin }^2\mathrm{\θ}}{R_{0}c}\right)t_m^2{f}_{\mathrm{\τ}}\}\end{array};$

S5、变参数扩展方位向非线性CS，将静止目标的多普勒质心和调频斜率分别进行均衡；包括以下子步骤：

S51、将同一距离单元的静止场景目标点的多普勒质心和多普勒调频斜率随方位向的空变分别建模为一次函数和二次函数：

f_Dc＝f_Dcref+at_m

$f_{Dr}=f_{Drref}+{\mathrm{bt}}_m+{\mathrm{dt}}_m^2$

其中，f_Dcref为静止场景参考点的多普勒质心，f_Drref为静止场景参考点的多普勒调频斜率，a、b和d分别为拟合系数；

S52、根据上述拟合系数，构造四次滤波和方位向非线性CS函数：

四次滤波函数为：

$H_{1F}\left(f_a\right)=\exp \left\{j\mathrm{\π}(G_3{f}_a^3+G_4{f}_a^4)\right\}$

方位向非线性CS函数为：

$H_{ANCS}\left(t_m\right)=\exp \left\{j\mathrm{\π}(e_2{t}_m^2+e_3{t}_m^3+e_4{t}_m^4)\right\}$

其中，G₃，G₄，e₂，e₃，e₄分别为四次滤波和非线性CS系数：

$\left\{\begin{array}{c}{e}_2=-2a\mathrm{\α}+(2\mathrm{\α}-1)f_{Drref}\\ e_3=\frac{2b(e_2+a)(e_2+f_{Drref})-f_{D3}{(a+e_2)}^2-e_2\[b(2e_2+a+f_{Drref})-f_{D3}(a+e_2)\]}{3{(f_{Drref}-a)}^2}\\ e_4=\frac{Q/4-(f_{Drref}-a)f_{Drref}^3{e}_2{G}_4}{a-f_{Drref}}\\ G_3=\frac{b(2e_2+a+f_{Drref})-f_{D3}(a+e_2)}{3{(f_{Drref}-a)}^2{e}_2{f}_{Drref}}\\ G_4=\frac{P/6-Q(a-f_{Drref})/4}{{(f_{Drref}-a)}^2{e}_2{f}_{Drref}^2(e_2+f_{Drref})}\end{array}\right.---\left(1\right)$

其中

$\begin{array}{c}P=-\[d{(e_2+f_{Drref})}^2-b^2(e_2+f_{Drref})\]-3f_{D3}b(a+e_2)\\ +3G_3{e}_2{\mathrm{bf}}_{Drref}(3f_{Drref}{e}_2-2{\mathrm{ae}}_2+f_{Drref}a)-3e_{3}b(e_2-2f_{Drref}+3a)\end{array}$

$Q=-3f_{D3}b+3G_3{e}_2{\mathrm{bf}}_{Drref}^2-3e_{3}b$

其中，f_D3为多普勒参数三阶项，α为可变参数；

S53、在式(1)中，令参数α分别等于α₁和α₂，得到两组系数，并分别构造两组四次滤波和非线性CS函数，第一组系数表示为：

$H_{1F,{\mathrm{\α}}_1}\left(f_a\right)=\exp \left\{j\mathrm{\π}(G_{3,{\mathrm{\α}}_1}{f}_a^3+G_{4,{\mathrm{\α}}_1}{f}_a^4)\right\}$

$H_{ANCS,{\mathrm{\α}}_1}\left(t_m\right)=\exp \left\{j\mathrm{\π}(e_{2,{\mathrm{\α}}_1}{t}_m^3+e_{4,{\mathrm{\α}}_1}{t}_m^4+e_{4,{\mathrm{\α}}_1}{t}_m^4)\right\}$

其中为第一组四次滤波和非线性CS系数，将α＝α₁带入(1)即可得出；

第二组系数表示为：

$H_{1F,{\mathrm{\α}}_2}\left(f_a\right)=\exp \left\{j\mathrm{\π}(G_{3,{\mathrm{\α}}_2}{f}_a^3+G_{4,{\mathrm{\α}}_2}{f}_a^4)\right\}$

$H_{ANCS,{\mathrm{\α}}_2}\left(t_m\right)=\exp \left\{j\mathrm{\π}(e_{2,{\mathrm{\α}}_2}{t}_m^2+e_{3,{\mathrm{\α}}_2}{t}_m^3+e_{4,{\mathrm{\α}}_2}{t}_m^4)\right\}$

其中为第二组四次滤波和非线性CS系数，将α＝α₂带入(1)即可得出；

S54、将回波信号S₂(t_m,f_τ)通过距离向IFFT和方位向FFT变换到距离时域方位频域，得到S₂(f_a,τ)；

S55、将回波信号S₂(f_a,τ)分别与和相乘，然后经过方位向IFFT变换到二维时域，得到和

S56、将和分别与和相乘，得到变参数扩展方位向非线性CS后的二维时域回波信号和

S6、变参数去斜处理：将经过扩展方位向非线性CS后的回波信号分别乘以对应的方位向去斜因子，并变换到方位向频域；

S7、进行静止场景杂波抑制：将两个回波的幅度相减，去除静止目标回波，实现静止目标杂波抑制；

S8、幅度检测和动目标提取：采用恒虚警处理，进行幅度检测，实现动目标提取。

2.根据权利要求1所述的斜视SAR地面动目标提取方法，其特征在于，所述的步骤S6具体实现方法为：将回波信号和分别乘以对应的方位向去斜因子，并变换到方位向频域，得到方位向聚焦的频域回波信号和

$S_{4,{\mathrm{\α}}_1}(f_a,\mathrm{\τ})={\mathrm{FFT}}_{azi}\{S_{3,{\mathrm{\α}}_1}(t_m,\mathrm{\τ})\·{\mathrm{\φ}}_{Der,{\mathrm{\α}}_1}\left(t_m\right)\}$

$S_{4,{\mathrm{\α}}_2}(f_a,\mathrm{\τ})={\mathrm{FFT}}_{azi}\{S_{3,{\mathrm{\α}}_2}(t_m,\mathrm{\τ})\·{\mathrm{\φ}}_{Der,{\mathrm{\α}}_2}\left(t_m\right)\}$

其中，FFT_azi{·}为方位向傅里叶变换，和分别为去斜因子：

${\mathrm{\φ}}_{Der,{\mathrm{\α}}_1}\left(t_m\right)=\exp \{-j\mathrm{\π}(e_{2,{\mathrm{\α}}_1}+f_{Drref})t_m^2\}$

${\mathrm{\φ}}_{Der,{\mathrm{\α}}_2}\left(t_m\right)=\exp \{-j\mathrm{\π}(e_{2,{\mathrm{\α}}_2}+f_{Drref})t_m^2\}.$

3.根据权利要求2所述的斜视SAR地面动目标提取方法，其特征在于，所述的步骤S7具体实现方法为：

$S_5(f_a,{\mathrm{\τ}}_0)=\left|S_{4,{\mathrm{\α}}_1}(f_a,{\mathrm{\τ}}_0)\right|-\left|S_{4,{\mathrm{\α}}_2}(f_a,{\mathrm{\τ}}_0)\right|.$

4.根据权利要求3所述的斜视SAR地面动目标提取方法，其特征在于，所述的步骤S8具体实现方法为：采用恒虚警处理，取一个自适应门限值δ，保持输出的虚警率不超过恒定值，进行幅度检测可以实现动目标提取，对于距离门τ₀，若满足：

|S₅(f_a,τ₀)|≥δ

则该距离门存在动目标，则对该动目标进行幅度检测，得到动目标提取结果。