CN109669183A

CN109669183A - 一种基于Keystone和时频变换的地球同步轨道SAR运动目标成像处理装置

Info

Publication number: CN109669183A
Application number: CN201811615054.3A
Authority: CN
Inventors: 于泽; 于金栋; 李春升
Original assignee: Beihang University
Current assignee: Beihang University
Priority date: 2017-12-27
Filing date: 2018-12-27
Publication date: 2019-04-23
Anticipated expiration: 2038-12-27
Also published as: CN109669183B

Abstract

本发明公开了一种基于时频变换和多普勒中心校正的地球同步轨道SAR运动目标成像处理装置，包括有，雷达参数输出模块，此模块输出雷达的工作参数及成像处理所需的成像参数；徙动校正及距离压缩模块，该模块根据模块一输出的雷达参数，对回波进行未知目标运动参数下的距离徙动校正，同时进行距离压缩操作；运动目标检测及参数估计模块，该模块根据运动目标目标和静止目标回波时频特性的差异检测运动目标，并根据频谱能量对多普勒中心进行校正，估计运动参数；运动补偿及成像处理模块，根据模块三得到的运动参数对运动信号进行补偿，方位压缩完成成像。

Description

一种基于Keystone和时频变换的地球同步轨道SAR运动目标成像处理装置

技术领域

本发明涉及一种有效的GEO SAR运动目标成像处理装置，具体来说是一种基于Keystone和时频变换的GEO SAR运动目标成像处理装置。

背景技术

合成孔径雷达(Synthetic Aperture Radar，SAR)是一种全天时、全天候的主动式信息获取系统，合成孔径雷达的全链路系统可以分为三个部分，即雷达载体平台上的回波数据获取部分、空间段数据传输部分以及地面段数据处理部分。合成孔径雷达系统通过发射和接收脉冲信号，获取合成孔径雷达原始回波数据。

军事热点地区观察、重要区域监视、地震多发区的地震预警等需要长时间对一固定区域观察的情况下，在轨的星载SAR由于其回归周期长，并不能胜任。近几年，地球同步轨道合成孔径雷达(GEO SAR，Geosynchronous orbit Synthetic Aperture Radar)成了热点，其轨道回归周期小，可以对某地区进行长时间的观测，且其具有很宽的测绘带宽度，能进行大范围观测，波束覆盖地面的时间也长，基于上述特点，GEO SAR在地面物体检测方面有很大前景。

用传统的如RDA算法(Range-Doppler Algorithm，距离多普勒算法)、CSA算法(Chirp Scaling Algorithm，线性调频信号变标算法)、WKA算法(omega-k算法，波束域算法)等算法对动目标进行成像时，其会偏离实际的位置而造成误差，导致对目标位置的错误估计，大部分SAR成像系统针对的是地面静目标，根据计算，动目标具有径向分量速度后，其在成像中的方位向位置会迁移，而具有方位向分速度时，会造成图像的散焦。GEO SAR下的方位向调频率比较低，其成像结果受该值估计精度的影响，且信杂比较低，需要相干积累才能检测出来。

GEO SAR对运动目标的检测和成像处理中，由于轨道高度位于地球同步轨道，覆盖范围大的同时引入了大量的杂波，因此信号难以检测。对于运动目标，在传统模式下，方位向信号难以相干积累导致了散焦的现象，因此针对检测方面而言有非相干方法，如J.Carretero＝Moya的《Application of the radon transform to detect smalltargets in sea clutter》中的Radon变换、B.D.Carlson《Search radar detection andtrack with the Hough transform》的Hough变换等基于幅度进行积累的检测方法，另一类为如WVD算法(Wigner-Ville Distribution，Wigner-Ville分布)、FrFT算法(FractionalFourier Transform,分数阶傅里叶变换)等时频处理方法，充分利用了信号的相位信息，进行相干积累来达到高的增益，方便低信杂比下信号的检测。对于该类信号的检测处理往往和参数估计密不可分，检测利用了信号的特性，而特性又反映了目标的参数。2002年HONG-BO SUN发表了《Application of the Fractional Fourier Transform to Moving TargetDetection in Airborne SAR》，将分数阶傅里叶变换用到了雷达运动目标的检测中，与常用的时频变换WVD进行了对比，并提出多目标CLEAN方法有效防止弱小目标的漏检。而在GEOSAR的应用场景下，距离徙动会使信号在低信杂比的环境中更难以检测，必须对距离徙动进行校正。Daiyin Zhu于2007年提出的《A Keystone Transform Without Interpolationfor SAR Ground Moving-Target Imaging》中采用频域处理来代替时域的插值有效降低了Keystone变换的运算量，尤其适用于大数据量的处理工作。Hang Yi在2009年论文《Imagingand Locating Multiple Ground Moving Targets Based on Keystone Transform andFrFT for Single Channel SAR System》也采用了Keystone和FrFT的方法，但文中未分析同一距离门多个目标的情况，仅靠FrFT不能完成这种情况的定位处理。

针对现有检测和参数估计方法的不足，本发明提出一种利用时频变换和频谱能量相结合的方法来估计完整的多普勒参数，同时为了在低信杂比下检测目标，利用了非参数化距离徙动校正方法，获取了良好的成像结果。

发明内容

本发明提出了一种基于Keystone和时频变换的GEO SAR运动目标成像处理装置，该装置针对GEO SAR运动目标非合作、信杂比低的特点，利用不依赖于参数估计的快速Keystone变换方法校正距离徙动，方便在后续的方位向脉冲压缩中获取高的能量积累；接着利用FrFT没有交叉项干扰的特点，估计运动目标信号的多普勒参数，使用频谱能量方法辅助确定同一距离门下多个目标的多普勒中心频率，方便后续的定位处理。最后根据估计出的多普勒参数进行运动的补偿，最终形成免受运动影响的图像。

本发明采用GEO SAR下的运动目标模型，长合成孔径时间使得目标积累效果更好，有利于信杂比的提升；Keystone变换在未知运动参数下可对整个场景进行距离徙动校正，尤其适合GEO SAR大场景的情况，采用基于频域变标思想的快速Keystone变换，降低运算量，提高运算效率；利用分数阶傅里叶变换这种时频分析方法检测运动目标，由于FrFT本身对多目标情况下的中心频率估计效果不佳，本发明采用基于频谱能量估计的多普勒中心频率测量方法，有效弥补了FrFT的缺点。

本发明的一种基于Keystone和时频变换的GEO SAR运动目标成像处理装置，该装置利用雷达参数输出模块来生成所需雷达参数，其特征在于还包括有：徙动校正及距离压缩模块、运动目标检测及参数估计模块和运动补偿及成像处理模块；

雷达参数输出模块，用于生成雷达工作参数f₁＝{Na,Nr,PRF,F_s,B_r,T_r,K_r,λ,V,t_d,t_center,η,τ,f_η,f_τ}；

徙动校正及距离压缩模块，用于对原始回波S_origin进行距离向傅里叶变换并进行距离向匹配滤波完成脉冲压缩，接着在距离频率和方位时间域进行变标处理，去掉方位时间和距离频率的一阶耦合，也就完成了非参数化的距离徙动工作；

运动目标检测及参数估计模块，用于接受经距离压缩和徙动校正的数据S_rc，对距离门信号进行分数阶傅里叶变换分析，进行信号的检测，用频谱能量发辅助估计多普勒参数，将各距离门中信号的多普勒参数作为输出结果；

运动补偿及成像处理模块，用于补偿运动，方位压缩形成最终图像。该模块接收运动目标检测及参数估计模块的参数，针对存在运动目标的距离门中的单个运动目标进行补偿成像，最后将所有目标在图像域叠加，形成完整的图像。

本发明采用GEO SAR下的运动目标模型，长合成孔径时间使得目标积累效果更好，有利于信杂比的提升；Keystone变换在未知运动参数下可对整个场景进行距离徙动校正，尤其适合GEOSAR大场景的情况，本文采用基于频域变标思想的快速Keystone变换，降低运算量，提高运算效率；利用分数阶傅里叶变换这种时频分析方法检测运动目标，由于FrFT本身对多目标情况下的中心频率估计效果不佳，本发明采用基于频谱能量估计的多普勒中心频率测量方法，有效弥补了FrFT的缺点。

本发明的优点在于：

(1)针对地球同步轨道SAR的低信杂比特点，利用分数阶傅里叶域的搜索功能实现信号的检测；

(2)为适应地球同步轨道SAR的数据量大的特点，采用快速Keystone变换非参数化校正距离徙动，计量量小，计算效率高

(3)弥补了传统分数阶傅里叶变换中对多目标多普勒中心频率估计的误差，采用频谱能量法辅助估计参数，估计更加准确，具有好的鲁棒性。

附图说明

图1是本发明的基于Keystone和时频变换的GEO SAR运动目标成像处理装置结构框图。

图2是分数阶傅里叶变换域检测结果图。

图3是本发明的基于Keystone和时频变换的GEO SAR运动目标成像处理装置的流程图。

图4是运动目标示意图。

图5A是距离压缩后的二维时域结果。

图5B是Keystone变换后的距离徙动结果。

图5C是经运动补偿后的成像结果。

图5D是同一距离门多目标成像结果图。

图6A是成像结果方位向切片图。

图6B是成像结果方位向切片插值图。

图7A是成像结果距离向切片图。

图7B是成像结果距离向切片插值图。

具体实施方式

下面将结合附图和实施例对本发明做进一步的详细说明。

参见图1、图3所示，本发明的一种基于Keystone和时频变换的GEO SAR运动目标成像处理装置，该装置由雷达参数输出模块、徙动校正及距离压缩模块、运动目标检测及参数估计模块和运动补偿及成像处理模块构成。本发明装置是对二维原始回波信息S_origin生成的运动目标成像进行处理。

(一)雷达参数输出模块

在本发明中，雷达参数输出模块用于产生雷达的工作参数，记为f₁＝{Na,Nr,PRF,F_s,B_r,T_r,K_r,λ,V,t_d,t_center,η,τ,f_η,f_τ}。该雷达参数输出模块的输入参数为二维原始回波信息S_origin。

N_a是发射脉冲总数；

N_r是距离采样点总数；

PRF是脉冲重复频率；

F_s是距离向采样率；

B_r是雷达工作带宽；

T_r是发射脉冲持续时间；

K_r是距离向调频率；

λ是工作波长；

V是雷达平台的飞行速度；

t_d为采样延时；

t_center为仿真中心时刻；

η为方位向时刻，且η＝[uη₁,uη₂,…uη_m,…,uη_Na]^T，m表示雷达发射的脉冲个数的标识号，uη₁表示第一个方位向时刻，uη₂表示第二个方位向时刻，uη_m表示第m个方位向时刻，且uη_Na表示最后一个方位向时刻，上角标T表示转置符号；

f_η为方位向频率，且m′表示雷达发射的频点个数的标识号，表示第一个方位频域频点，表示第二个方位频域频点，表示第m′个方位频域频点，且表示最后一个方位频域频点，在本发明中，m′与m选取的个数相同；

τ为距离向时刻，且τ＝[uτ₁,uτ₂,…uτ_k,…,uτ_Nr]，k表示距离时间采样数的标识号，uτ₁表示第一个距离向时刻，uτ₂表示第二个距离向时刻，uτ_k表示第k个距离向时刻，且(即从接收窗开启时间起的第k个距离向采样时刻)，uτ_Nr表示最后一个距离向时刻；

f_τ为距离向频率，且f_τ＝[uf_τ,1,uf_τ,2,…uf_τ,k′,…,uf_τ,Nr]，k′表示距离频率采样数的标识号，uf_τ,1表示第一个距离向频率，uf_τ,2表示第二个距离向频率，uf_τ,k′表示第k′个距离向频率，且uf_τ,Nr表示最后一个距离向频率；k′与k选取的个数相同。

GEO SAR处理器接收到的回波信号S_origin是二维数据，采用矩阵形式表示为：

n是一个变量，取值范围1至N_a，表示为方位向第n个点；

m是一个变量，取值范围1至N_r，表示为距离向第m个点；

表示沿方位向上采集到的第1个方位时刻的距离向上的第1个采样点；

表示沿方位向上采集到的第1个方位时刻的距离向上的第2个采样点；

表示沿方位向上采集到的第1个方位时刻的距离向上的第m个采样点；

表示沿方位向上采集到的第1个方位时刻的距离向上的第N_r个采样点；

表示沿方位向上采集到的第2个方位时刻的距离向上的第1个采样点；

表示沿方位向上采集到的第2个方位时刻的距离向上的第2个采样点；

表示沿方位向上采集到的第2个方位时刻的距离向上的第m个采样点；

表示沿方位向上采集到的第2个方位时刻的距离向上的第N_r个采样点；

表示沿方位向上采集到的第n个方位时刻的距离向上的第1个采样点；

表示沿方位向上采集到的第n个方位时刻的距离向上的第2个采样点；

表示沿方位向上采集到的第n个方位时刻的距离向上的第m个采样点；

表示沿方位向上采集到的第n个方位时刻的距离向上的第N_r个采样点；

表示沿方位向上采集到的第N_a个方位时刻的距离向上的第1个采样点；

表示沿方位向上采集到的第N_a个方位时刻的距离向上的第2个采样点；

表示沿方位向上采集到的第N_a个方位时刻的距离向上的第m个采样点；

表示沿方位向上采集到的第N_a个方位时刻的距离向上的第N_r个采样点。

(二)距离徙动校正及距离压缩模块

在本发明中，徙动校正及距离压缩模块一方面用于接收雷达参数输出所述f₁，另一方面接收二维原始回波信息S_origin，然后进行快速Keystone变换完成距离徙动校正。所述徙动校正及距离压缩模块的处理方法如下：

步骤2-1，距离压缩；

是将二维原始回波信息S_origin进行距离向傅里叶变换和距离向匹配滤波，得到距离向匹配滤波后的方位时间—距离频域信，记为号；

A₀表示雷达后向散射系数。

rect表示矩形函数形式。

ω_a表示方位向窗函数形式。

η_c为波束中心穿越时刻。

j表示虚部。

f₀为雷达载波频率。

c表示光速。

R(η)表示斜距历程。

在本发明中，若对公式(1)直接进行方位向傅里叶逆变换，则得到距离压缩后的信号，记为S_rc(τ,η)，即公式(2)。从所述S_rc(τ,η)中可以看出一维脉压后的信号存在着距离徙动，无法将能量集中于一个距离门中。

p_r表示距离脉冲压缩后的幅度形式。

经公式(2)处理后的距离压缩后的二维时域信号，如图5A所示。

假设运动目标只有平动，不存在高阶运动，由基本的GEO SAR与平动目标之间的几何关系，得到距离历程R(η)为：

R₀为零时刻雷达与目标的距离。

v_a为方位向速度。

v_c为径向速度。

a_c为径向加速度。

由距离历程R(η)和方位向时刻η之间的关系简化公式(3)后为：

步骤2-2，Keystone变换；

在方位时域—距离多普勒域进行Keystone变换，Keystone变标为：

具体操作为沿着距离向对每个固定距离向频率的距离门信号进行方位时间的变标处理，变标前方位时间为η，变标后方位时间为t_a，首先选择一距离门，此时距离向频率f_τ是一个固定值，在该距离门内设定一组均匀采样的时间为：

t_a＝[ut_a,1,ut_a,2,…ut_a,k,…,ut_a,Na] (6)

ut_a,1表示第一个采样时刻，ut_a,2表示第二个采样时刻，ut_a,k表示第k个采样时刻，k表示采样时刻标识，ut_a,Na标识最后一个采样时刻。

公式(6)经过公式(5)的变标处理对应的变标前的位置为：

对应着公式(7)中的位置在原信号中进行插值，找到公式(7)中采样时刻对应的信号，组成新的采样时刻即完成插值操作，在每一个距离门内均进行这样的插值操作，每次根据距离门的位置变换f_τ的值，则最终完成对整个场景的Keystone变换。

步骤2-3，快速Keystone变换；

在本发明中，快速Keystone变换是时域插值操作转换为频域的相乘。具体操作如下：

构建滤波器

第一滤波器

第二滤波器

第三滤波器

第四滤波器

所述变标后的结果记为：

FFT表示沿方位向傅里叶变换。

IFFT表示沿方位向逆傅里叶变换。

按照公式(12)逐步对输入的方位时间—距离频率信号进行滤波操作和傅里叶变换、逆傅里叶变换，最后对公式(12)的结果进行方位向傅里叶逆变换即可得到经过距离压缩和徙动校正的二维时域，记为S_rcmc(τ,η)：

在本发明中，利用徙动校正及距离压缩模块对运动目标进行距离压缩和非参数化的距离徙动校正，能够在未知目标运动参数情况下，将能量集中于同一距离门，处理结果如图5B所示。

(三)运动目标检测及参数估计模块

在本发明中，运动目标检测及参数估计模块用于接收经距离压缩和徙动校正的二维时域信号S_rcmc(τ,η)，然后检测每个距离门中的信号并进行参数的估计，具体操作方法如下：

步骤3-1，分数阶傅里叶变换；

针对运动目标和杂波不同特点，利用分数阶傅里叶变换对距离门内的信号进行处理。分数阶傅里叶变换操作的输入数据是该距离门内的时域信号，下面给出分数阶傅里叶变换的操作方式：

F_p(u)为分数阶傅里叶变换的结果；

f(η)为输入的一维时间数据，在本发明中，所述f(η)是从S_rcmc(τ,η)中选取的一距离门的数据；

d为积分形式；dη为对方位向时刻的积分；

K_p(u,η)是分数阶傅里叶变换的核函数；所述K_p(u,η)为：

其中核函数系数α为旋转角度，且有p为分数阶傅里叶变换的阶数；δ为狄拉克函数，u为分数阶傅里叶变换的位置，n为自然数；

联合公式(14)和公式(15)，得到完整核函数F_p(u)为：

结合公式(4)和公式(13)，得出一距离门的信号f(η)为：

使用分数阶傅里叶变换对公式(17)进行处理，得到分数阶傅里叶域信号F_p(u)：

步骤3-2，分数阶傅里叶域参数搜索；

扫描分数阶傅里叶变换域和转角所形成的能量图，如附图2所示的三维图。搜索图中能量最大点(尖峰)，估计尖峰所在的位置，获得二维坐标(α,u)。

当公式(18)的相位中参数满足公式(19)时，F_p(u)具有最大的能量聚集，幅值为最大。则公式(4)中的a,b取值为：

若将公式(17)视为线性调频波的形式，则其对应的多普勒中心频率f_d和多普勒调频率f_r为：

获取分数阶傅里叶域u和转角α以及对应的|F_p(u)|所形成的三维图中最大尖峰值的位置为：

(α₀,u₀)＝argmax|F_p(u)| (21)

步骤3-3，多普勒补偿；

由公式(21)的位置反演信号的多普勒参数及目标的运动参数，即联立公式(20)、公式(21)可得该距离门信号的多普勒参数；

联立公式(3)、公式(4)、公式(19)和公式(21)可得出目标的运动参数。

在本发明中，但在信号检测和成像处理中只需要估计出多普勒中心频率f_d、多普勒调频率f_r或公式(19)的a、b的值，即可完成补偿。

步骤3-4，频谱能量法估计多普勒中心；

采用频谱能量法估计信号多普勒中心参数，该部分输入数据为分数阶傅里叶变换所处理的同一距离门的信号S_rcmc(τ,η)。首先对所述S_rcmc(τ,η)进行方位向傅里叶变换，然后求绝对值的平方。

E(f_η)＝|FFT_η[S_rcmc(τ_k,η)]|² (22)

构建多普勒能量滤波器：

在本发明中，除以外的称为其他。

采用频域相乘的形式来代替时域卷积，有：

F_out(f_η)＝IFFT{FFT[E(f_η)]FFT[F_pb(f_η)]} (24)

取F_out(f_η)中的降交零点即为功率等分频点，即信号的多普勒中心。

步骤3-5，CLEAN算法处理多目标；

当同一距离门存在多个目标的信号时，采取CLEAN算法，具体如下：

采用步骤3-1～3-4估计参数后，在最高尖峰对应的转角下采用窄带滤波的方式在分数阶傅里叶域滤除目标，得到：

p₀为最高尖峰对应的分数阶傅里叶变换的阶数，为对信号进行以p₀为阶数的分数阶傅里叶逆变换，为对信号进行以p₀为阶数的分数阶傅里叶变换。H_p(u)是分数阶域的窄带滤波器，剔除已检测出的运动目标信号，并将所述窄带滤波器H_p(u)滤除的信号保存以便运动补偿及成像处理模块使用，记为f_res(η)为剔除检测出的信号后的剩余信号。将按同样的逆分数阶傅里叶变换操作，得到检测到的信号的时域信号f_detect(η)。

步骤3-6，在同一距离门信号中，不断重复步骤3-5直到检测到所有目标位置；

步骤3-7，沿距离向逐距离门采取步骤3-1～步骤3-6，即可完成所有运动目标的检测和参数估计。

(四)运动补偿及成像处理模块

在本发明中，运动补偿及成像处理模块用于接收运动目标检测及参数估计模块中分数阶傅里叶域对应的已检测出的运动目标信号f_detect(η)，且已知其多普勒中心频率f_d和多普勒调频率f_r。采用匹配滤波方式进行方位向脉冲压缩，而后将所有运动目标的成像结果叠加即得最终图像。

在本发明中，构建校正运动的匹配滤波器：

以单个运动目标成像结果为：

s_single(η)＝IFFT{FFT[f_detect(η)]H_a(f_η)} (27)

在本发明中，将同一距离门中检测并校正的单目标成像结果叠加起来，如图5D所示，同时逐距离门计算，最后将所有单成像结果叠加起来，形成最终图像，如图5C所示。

i为单个运动目标的个数的标识号。

本发明的一种基于Keystone和时频变换的GEO SAR运动目标成像处理装置，对于信号的处理为(参见图3所示)：二维回波信号的距离向傅里叶变换→距离向匹配滤波→距离向逆傅里叶变换完成距离压缩→Keystone变换→选择一距离门信号→进行分数阶傅里叶变换同时采取基于频谱能量估计的方法估计多普勒中心频率→根据变换结果判断有无运动目标→若存在目标，求运动参数进行定位、成像，若无目标则处理下一距离门信号→处理完所有距离门后即得运动目标图像。

实施例1

仿真实施参数以表1所示，本发明方法主要研究运动目标的检测与成像，因此设置场景大小为2km(方位向)×2km(距离向)，为了让雷达运行轨迹尽量符合直线，设置观测点在赤道附近且为了防止高分辨率带来的复杂的空变性而选择了低分辨率来进行实验。

表1雷达参数

参数	取值
		轨道倾角(度)	60
轨道偏心率	0
		轨道半长轴(米)	42164431.345
近心点角距(度)	0
		升交点赤经(度)	91
过近心点时刻(秒)	0
		仿真中心时刻t<sub>center</sub>(s)	46300
合成孔径时间(s)	25
		波长λ(m)	0.24
脉宽T<sub>r</sub>(μs)	50
		信号带宽B<sub>r</sub>(MHz)	8
采样率F<sub>s</sub>(MHz)	16
		重频PRF(Hz)	120

实施例1中运动目标场景示意图如图4所示：图中圆形、正方形、菱形、三角形目标分别为位于不同位置，具有不同速度的目标，五角星形目标为四个目标运动的终点，运动参数如表2所示。

表2运动目标参数

经本发明方法处理所得的相关指标结果如表3所示：

表3运动目标结果

目标名称	定位位置	方位向速度	距离向速度
				目标1	(96.73,20.43)	3.4467m/s	2.8270m/s
目标2	(-300.2,-388.2)	0.5911m/s	2.0149m/s
				目标3	(447,388.2)	5.1273m/s	0.0153m/s
目标4	(403.6,-61.29)	0.8621m/s	-2.6278m/s

利用估计的运动参数进行成像，方位向切片如图6A所示，图6B为图6A的插值结果。距离向切片如图7A所示，图7B为图7A的插值结果。经校正后的运动目标成像结果成像指标均符合要求。图7A中的分辨率为斜距分辨率，经几何投影之后，与设定分辨率相同。

综上结果，本发明在GEO SAR分辨率为50m的情况下，定位误差在100m内，测速误差在0.65m/s内，且能处理位于同一距离门和位于不同距离门的多运动目标，形成集群运动目标的检测与成像，成像结果符合要求。

Claims

1.一种基于Keystone和时频变换的地球同步轨道SAR运动目标成像处理装置，该装置利用雷达参数输出模块来生成所需雷达参数，其特征在于还包括有：徙动校正及距离压缩模块、运动目标检测及参数估计模块和运动补偿及成像处理模块；

2.根据权利要求1所述的一种基于Keystone和时频变换的地球同步轨道SAR运动目标成像处理装置，其特征在于徙动校正及距离压缩模块的处理方法为：

步骤2-1，距离压缩；

步骤2-2，Keystone变换；在方位时域—距离多普勒域进行Keystone变换，Keystone变标为

步骤2-3，快速Keystone变换；所述变标后的结果记为对输入的方位时间—距离频率信号进行滤波操作和傅里叶变换、逆傅里叶变换，然后进行方位向傅里叶逆变换即可得到经过距离压缩和徙动校正的二维时域

是将二维原始回波信息S_origin进行距离向傅里叶变换和距离向匹配滤波，得到距离向匹配滤波后的方位时间—距离频域信，记为假设运动目标只有平动，不存在高阶运动，由基本的GEO SAR与平动目标之间的几何关系，得到距离历程R(η)为由距离历程R(η)和方位向时刻η之间的关系得到R(η)＝R₀+aη+bη²，a＝-v_c，

3.根据权利要求1所述的一种基于Keystone和时频变换的地球同步轨道SAR运动目标成像处理装置，其特征在于：运动目标检测及参数估计模块用于接收经距离压缩和徙动校正的二维时域信号S_rcmc(τ,η)，然后检测每个距离门中的信号并进行参数的估计，具体操作方法如下：

步骤3-1，分数阶傅里叶变换；

步骤3-2，分数阶傅里叶域参数搜索；

步骤3-3，多普勒补偿；

步骤3-4，频谱能量法估计多普勒中心；

步骤3-5，CLEAN算法处理多目标；

4.根据权利要求1所述的一种基于Keystone和时频变换的地球同步轨道SAR运动目标成像处理装置，其特征在于：运动补偿及成像处理模块用于接收运动目标检测及参数估计模块中分数阶傅里叶域对应的已检测出的运动目标信号f_detect(η)，且已知其多普勒中心频率f_d和多普勒调频率f_r；采用匹配滤波方式进行方位向脉冲压缩，而后将所有运动目标的成像结果叠加即得最终图像

5.根据权利要求1所述的一种基于Keystone和时频变换的地球同步轨道SAR运动目标成像处理装置，其特征在于：在GEO SAR分辨率为50m的情况下，定位误差在100m内，测速误差在0.65m/s内，且能处理位于同一距离门和位于不同距离门的多运动目标，形成集群运动目标的检测与成像。