用于高分辨率宽测绘带SAR系统的SAR/GMTI处理方法研究

邢孟道 孙光才 李学仕

(西安电子科技大学雷达信号处理国家重点实验室 西安 710071)

用于高分辨率宽测绘带SAR系统的SAR/GMTI处理方法研究

邢孟道 孙光才 李学仕*

(西安电子科技大学雷达信号处理国家重点实验室 西安 710071)

随着科技的不断发展，高分辨宽测绘(High-Resolution and Wide-Swath, HRWS)SAR成像已经越来越受到人们的关注，具有地面动目标指示(Ground Moving Target Indication, GMTI)功能的SAR成像系统因其具有对静止场景进行高分辨成像和动目标检测能力在很多军用和民用领域受到广泛运用。具有HRWS成像能力的沿方位多通道SAR系统，其可以有效地解决高分辨和低脉冲重复频率(Pulse Repetition Frequency, PRF)的矛盾，该矛盾在HRWS成像处理过程中经常遇到。由于方位空域自由度可以用来进行杂波抑制，因此多通道构型具有提供GMTI潜能。该文提出一种新的杂波抑制和动目标的成像方法，使得在低PRF的HRWS系统进行SAR成像的同时可以完成动目标的检测与成像处理，而不需要单独的高PRF系统操作模式。

合成孔径雷达；高分辨宽测绘；动目标检测；动目标成像；同时

1 引言

传统的杂波抑制算法主要包括DPCA(Displaced Phase Center Antenna, DPCA)技术[1,2]，ATI(Along Track Interferometry, ATI)技术[3–5]以及空时自适应处理 (Space-Time Adaptive Processing, STAP)技术[6]。对于多通道SAR系统，传统算法中的绝大部分要求每个通道接收到杂波回波不存在多普勒谱模糊。然而，对于具有GMTI能力的HRWS MC-SAR系统，其单个通道接收到杂波回波是存在多普勒模糊。虽然文献[7]提出的方法可以对HRWS MC-SAR系统进行杂波抑制，但其要求通道数必须大于多普勒模糊数的两倍且DPCA条件必须严格满足，这将明显地增加系统复杂度。本文针对具有GMTI能力的HRWS MC-SAR系统，提出了一种新颖的杂波抑制算法。在每个通道回波存在多普勒谱模糊时，利用方位的CFT (Chirp-Fourier Transform, CFT)[8]操作来对每个通道回波的方位信号进行压缩。因此，可以在距离脉压方位CFT域(RC-CFT domain)中获得每个通道回波的粗聚焦SAR图像。通常，在杂波抑制之前，动目标参数是未知的。为了在杂波抑制过程中对所有的动目标尽可能保持最大能量，假设动目标处在动目标区域的中间位置(动目标区域定义为相邻两个杂波模糊分量之间的区域)。接着，针对所有动目标区域，构造出一系列空域滤波器，同时在杂波方向形成零点。利用构造的空域滤波器，对杂波进行抑制和提取动目标回波信号。

经过杂波抑制之后，动目标成像处理是另外一个重要的问题，其包括动目标参数估计，特别是斜距速度估计。在本文中，提出了一种新的动目标成像算法。在距离频率域对每个通道回波沿方位向进行CFT之后，QRCM (Quadratic Range Cell Migration, QRCM)得到有效地校正并只存在动目标的线性RCM (Range Cell Migration)。在本文中，首先将动目标回波变换到距离频率方位时域。接着，沿距离频率方向构造两组信号和利用文献[9]所提出算法进行基带多普勒中心估计。利用多视互相关的多普勒中心模糊数求解方式[10]来获得准确的动目标斜距速度和对动目标线性RCM进行有效地校正。最后，经过方位逆傅里叶变换之后，可以获得聚焦好的动目标图像。

本文内容安排如下：第2节给出了HRWS MCSAR系统的几何关系；第3节推导了动目标和杂波的粗聚焦图像及复响应；第4节介绍了杂波抑制和动目标检测方法；第5节给出了动目标成像和斜距速度估计方法；第6节对实测的动目标HRWS MCSAR数据进行了处理；第7节对全文进行了总结。

2 HRWS MC-SAR系统的几何关系

对于正侧视模式，HRWS SAR系统的几何关系如图1(a)所示，其中X轴为方位方向，其表示平台的速度方向，Y轴指向左边，Z轴方向背离地球中心。雷达平台以恒定前向速度v沿方位向运动，θ为俯仰角，φ为锥角，Wg为测绘带宽度，Rb为地距。通常，把整个天线沿方位分割成N个子孔径，其如图1(b)所示且以7个通道为例，即N=7。中间的子孔径天线发射宽波束的线性调频信号，所有的子孔径天线同时接收回波信号。每相邻两个子孔径天线沿方位的距离为2d。HRWS MC-SAR系统中高的方位分辨率由单个子孔径天线尺寸决定的，而利用整个天线的长度来抑制多普勒谱模糊。通过补偿一个常数相位，每一组分开的发射和接收通道可以等效于一个对应于等效相位中心(Effective Phase Center, EPC)的通道自发自收。对于EPC位置，相邻的通道之间的间距为d。令dn=n.d为第n通道和参考通道之间的EPC位置距离(n=1，2，...，N)。在本文中，将第1个通道作为参考通道。

在HRWS MC-SAR 系统数据的录取过程中，假设动目标以恒定速度运动。图1(c)为正侧视模式下平台与动目标在斜距平面上的几何关系，其中vr与va分别表示动目标在径向方向和沿航向方向的速度。t表示方位慢时间，R0为雷达平台与动目标之间(或是成像场景中杂波与雷达平台之间)的最短斜距，RT,n(t)为第n通道与动目标之间的瞬时斜距。在雷达的观测时间Ta内，动目标由a点位置运动到b点位置。本文主要讨论HRWS MC-SAR系统中有关动目标的杂波抑制和成像处理。图2 所示为本文所提出算法的流程图，其主要步骤包括杂波抑制、参数估计和成像处理。此外，Wj为提取动目标回波的权值矢量，ST,j为动目标粗聚焦的图像，为估计多普勒参数构造的两组数据。

3 动目标和杂波的粗聚焦图像及复响应

CFT作为FT的2阶推广，比原始的FT具有更好处理线性调频信号的能力。信号s(t)的CFT可以表示为[8]：

其中φ0为常数相位，f1为CFT频率谱的平移，ka为分析信号采用的线性调频斜率，fa为CFT的频率。

图1 HRWS MC-SAR系统的几何关系Fig. 1 Geometry of an MC-HRWS SAR system

图2 本文所提HRWS MC-SAR中杂波抑制和动目标成像算法流程图Fig. 2 Flowchart of the proposed clutter suppression and moving target imaging algorithm for an MC-HRWS SAR system

在式(1)中，CFT可以认为是信号s(t)经过一个2次多项式相位补偿之后的FT。式(1)对应的CFT逆变换可以表示为：

一个点目标方位信号的CFT频率谱仅仅分布在几个频率单元，这与傅里叶频率谱存在很大区别。此外，CFT频率fa和锥角之间的关系可以表示为图3(a)所示为一个点目标方位信号的FT频率谱和CFT频率谱，其中多普勒带宽大于PRF。从图中可以看到模糊的FT频率谱分布到整个多普勒频率区域而相应的CFT频率谱仅仅分布在一些CFT频率单元。然而，一个点目标CFT频率谱的模糊特性意味着位置模糊，如图3(b)所示。图3(b)为杂波的模糊CFT频率谱和无模糊CFT频率谱之间关系，其中宽的虚线代表杂波无模糊的CFT频率谱，宽实线为杂波模糊的CFT频率谱。图3(b)是以两个点目标(杂波)为例。

图4所示为动目标方位信号的CFT频谱，其中黑的宽虚线表示杂波无模糊的CFT频率谱，蓝色细虚线为动目标无模糊CFT频率谱。绿色宽的虚线和红色细虚线分别表示杂波和动目标模糊的CFT频率谱。此外，动目标区域定义为两个杂波模糊分量之间的区域，图4所示为第1个动目标区域和第l个动

目标区域。动目标CFT频率谱位于动目标区域。每个通道杂波回波的粗聚焦图像位于所有CFT频率单元及其相应的模糊分量。在杂波抑制和动目标检测过程中，要求所有杂波方向得到抑制而动目标区域方向能量尽可能地保持。

图3 杂波回波方位信号的CFTFig. 3 CFT for the azimuth signal of clutter echoes

图4 杂波回波方位信号的CFTFig. 4 CFT for the azimuth signal of clutter echoes

4 杂波抑制和动目标检测

在文献[11]中，利用多个多普勒方向约束方式来抑制多普勒模糊分量和重构方位多普勒无模糊信号，其在多普勒模糊分量方向上形成零点。受多个多普勒方向约束方式的启发，为了在CFT频率域中对杂波进行抑制，零点可以在杂波方向及其相应的模糊方向上形成。此外，动目标潜在区域的能量需要得到保持(在CFT频率域中，由于动目标位于动目标区域如图4所示，动目标潜在方向表示动目标区域方向)。因此，杂波抑制的最优准则可以表示为：

在式(3)中，当动目标存在时，RX(fa)为在CFT频率单元fa的多通道输出对应的协方差矩阵，即RX(fa)可以表示为：

在杂波的抑制过程中，杂波及其模糊分量方向是已知的且需要形成零点。接着，可以得到导向矢量矩阵：

式(5)中，对于CFT频率单元fa,T，l(fa)=T.c，l(fa)为动目标方向(通常，在动目标参数估计之前，动目标方向是未知的。为了在杂波抑制过程中对所有动目标保持高的信号能量，认为动目标位于动目标区域的中间位置)，ac,l(fa)为杂波方向及其模糊分量方向。通常，动目标方向需要尽可能保持最大能量和杂波方向及其模糊分量方向需要保持尽可能的最小能量，由此可以得到Wl(fa)为在CFT频率谱中提取位于第l个动目标区域动目标的权值矢量。式(3)所示最优化问题的最优解可以表示为：

包含动目标的多通道SAR系统回波信号为Sc(fa)+ST(fa)。假设动目标的CFT频率谱位于动目标的第l个区域位置，如图5所示。杂波抑制过程可以表示为：

图5所示为杂波抑制过程，其中在成像场景中包含一个动目标。从图中看到在杂波方向及其模糊分量方向上形成零点。经过杂波抑制之后，可以提取出动目标回波。

当动目标多普勒中心fdc满足动目标的CFT频率谱和杂波的CFT频率谱处在相同位置区域，其将在杂波抑制过程中同样也得到抑制。由于v其中vr为盲速。注意到当系统的PRF越低，更多的斜距速度将位于盲速位置。

图5 本文所提出的杂波抑制算法处理过程Fig. 5 Clutter suppression processing of the proposed approach

图6 ISTAP方式的杂波抑制处理和动目标信号提取Fig. 6 Clutter suppression processing and moving target signal extraction of the ISTAP approach

图6所示为ISTAP[12–14]方式的杂波抑制处理过程和动目标信号提取。在图6中注意到当动目标的1个信号分量需要从1个多普勒频率单元提取出来时，所有杂波方向和动目标模糊分量方向需要形成零点。经过在每个多普勒频率单元提取动目标信号分量和重排之后，可以得到动目标无模糊的多普勒谱。由于动目标多普勒谱的模糊数和杂波一样，因此所需的通道数至少应为杂波谱模糊数的两倍。然而，对于本文所提出的杂波抑制算法，需要的通道数只需大于杂波的模糊数。

5 动目标成像和斜距速度估计

在这一小节中，将对动目标的成像算法进行讨论，其包括动目标斜距速度估计。经过杂波抑制之后，对提取的动目标回波进行距离FT和方位逆FT。动目标回波可以表示为：

图7 利用提取的动目标回波构造的两组信号Fig. 7 Construction of two groups signals for the extracted moving target echo

其中Mamb为多普勒中心模糊数。基于多视互相关的多普勒中心解模糊算法[10]，可以得到等效的斜距速度

线性RCM校正函数可以构造为：

经过距离逆傅里叶变换和方位傅里叶变换，聚焦的动目标图像可以表示为：

从式(14)中，注意到利用方位全孔径回波信号可以获得方位全孔径分辨率。对于几何分辨率，其方位分辨率为其中Ba为方位多普勒带宽，ka为方位调频率，距离分辨率为

6 实测的动目标HRWS MC-SAR数据

在这一小节中，将讨论在实测的机载HRWS MC-SAR数据中动目标的处理结果。在实验过程中，录取得到实测的机载X波段的HRWS MCSAR数据，其包括动目标回波，主要的系统参数如表1所示。在系统的工作过程中，中间通道发射线性调频信号，所有通道接收场景回波信号。由于多普勒带宽为1000 Hz和多通道系统的PRF为200 Hz，对于杂波回波的方位多普勒信号其多普勒模糊数近似为5。此外，在数据的录取过程中成像场景包括多个动目标。

对于获得的HRWS MC-SAR回波，利用文献[11]提出的HRWS MC-SAR成像算法来重构杂波方位多普勒无模糊信号。进而得到杂波聚焦的成像结果，如图8所示。从图8中，可以看到两个动目标及模糊分量信号，其分别由红色实圆圈和绿色虚圆圈标出。通常，利用传统的HRWS MC-SAR成像算法[11]进行处理，动目标无法取得很好聚焦效果。为了对动目标进行聚焦处理，利用本文所提出的杂波抑制和动目标成像处理算法来提取和重新聚焦动目标回波。

在动目标回波的提取过程中，首先进行杂波抑制，其过程如图9所示。图9(a)所示为录取的HRWS MC-SAR系统单通道回波在RC-CFT域的信号形式，其中两个动目标的回波被标出来，且杂波回波具有高的幅度能量值。接着，得到杂波的抑制结果，如图9(b)所示，杂波得到抑制，动目标回波得到保留。为了进一步说明杂波抑制性能，提取出Line 1和Line 2在图9(a)和图9(b)中的信号进行比

较。图9(c)所示为Line 1所标出信号的比较结果，其中红色虚线为杂波抑制之前的信号，黑色实线为杂波抑制之后的信号。图9(d)所示为Line 2所标出信号的比较结果。从图9(c)和图9(d)中注意到在假设动目标能量不变的情况下杂波可以抑制15 dB，这说明本文所提出的杂波抑制算法是有效的。

表1 实测数据的系统参数Tab. 1 System parameters for the real data

图8 利用参考文献[12–14]提出的方位多普勒信号重构算法进行重构之后场景回波的成像结果Fig. 8 Imaging result for the scene echoes using the azimuth Doppler signal reconstruction approach proposed in Refs. [12–14]

图9 RC-CFT域杂波抑制过程Fig. 9 Clutter suppression process in the RC-CFT domain

经过杂波抑制之后，可以提取出动目标回波。图10(a)所示为提取出来在CR-CFT域中第1个动目标(Target 1)回波信号，图10所示为第1个动目标的成像过程。图10(b)所示为该提取动目标回波在距离脉压和方位慢时间域的信号形式，其存在明显的RCM。经过参数估计和线性RCM校正之后，回波如图10(c)所示，估计得到动目标斜距速度为7.2063 m/s。图10(c)中，动目标的回波校正到一条水平直线上，这说明估计得到的斜距速度是准确的。图10(d)所示为聚焦的动目标图像。类似地，第2个动目标(Target 2)的成像处理过程如图11所示，且聚焦好的动目标图像如图11(d)所示。因此，对于HRWS MC-SAR系统，本文所提出的动目标成像算法是有效的。

图10 第1个动目标的成像过程Fig. 10 Imaging process of the first moving target

7 结束语

利用沿方位的多通道自由度进行多普勒模糊分量抑制，允许多通道雷达系统获取高分辨宽测绘带成像。同样，杂波抑制使得人们可以进行动目标检测和成像处理。在本文中，主要研究了具有GMTI功能的HRWS MC-SAR系统。HRWS MCSAR系统中每个通道杂波回波存在多普勒模糊，传统的杂波抑制算法，如DPCA, EDPCA及基于STAP的杂波抑制算法，不能取得很好的杂波抑制性能。因此，针对HRWS MC-SAR系统，本文提出了一种稳健的杂波抑制算法。经过杂波抑制之后，利用基于CFAR和GLRT检测器可以将动目标回波提取出来，进而完成对动目标的聚焦处理。本文所提出的算法仅要求通道数大于杂波谱的多普勒模糊数。与此同时，所提出算法不仅适用于均匀的DPC采样，同样也适用于非均匀采样。然而，本

文所提方法对于盲速动目标是无能为力的，与常规STAP方法一样，在空时平面中，盲速动目标的空时谱分布与静止背景杂波是重合的，在抑制杂波的同时，动目标能量也会被抑制掉。最后，实测数据处理结果验证了本文所提算法的有效性。

图11 第2个动目标的成像过程Fig. 11 Imaging process of the sceond moving target

[1]Gierull C H and Sikaneta I C. Raw data based two-aperture SAR ground moving target indication[C]. Proceedings 2003 IEEE IGARSS’03, Toulouse, France, 2003, 2: 1032–1034.

[2]Wang H S C. Mainlobe clutter cancellation by DPCA for space-based Radars[C]. Aerospace Applications Conference, Creasted Butte, CO, USA, 1991: 124–128.

[3]Shen Chiu. Application of fractional Fourier transform to moving target indication via along-track interferometry[J].EURASIP Journal on Applied Signal Processing, 2005(20): 3293–3303.

[4]Breit H, Eineder M, Holzner J,et al.. Traffic monitoring using SRTM along-track interferometry[C]. Proceedings 2003 IEEE IGARSS’03, Toulouse, France, 2003, 2: 1187–1189.

[5]Suchandt S, Runge H, Breit H,et al.. Automatic extraction of traffic flows using TerraSAR-X along-track interferometry[J].IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(2): 807–819.

[6]Klemm R. Introduction to space-time adaptive processing[J].Electronics&Communication Engineering Journal, 1999, 11(1): 5–12.

[7]Zhang S, Xing M, Xia X,et al.. A novel moving target imaging algorithm for multi-channel in azimuth HRWS SAR based on local maximum-likelihood minimum entropy in range frequency[J].IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(11): 5333–5347.

[8]Xia X G. Discrete chirp-Fourier transform and its application to chirp rate[J].IEEE Transactions on Signal Processing, 2000, 48(11): 3122–3133.

[9]Madsen S N. Estimating the Doppler centroid of SAR data[J].IEEE Transactions on Aerospace and Electronic Systems, 1989, 25(2): 134–140.

[10]Cumming I G, Kavanagh P F, and Ito M R. Resolving the Doppler ambiguity for spaceborne synthetic aperture radar[C]. Proceedings of IGARSS, Zurich, Switzerland, 1986: 1639–1643.

[11]Zhang S, Xing M, Xia X,et al.. Multi-channel HRWS SAR imaging based on range-variant channel-calibration and multi-Doppler-direction restriction ambiguity suppression[J].IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(9): 4306–4326.

[12]Delphine C and Ishuwa S. “Optimum GMTI processing for space-based SAR/GMTI systems-Theoretical derivation”[C]. 2010 8th European Conference on Synthetic Aperture Radar, Aachen, Germany, 2010: 390–393.

[13]Cerutti-Maori D and Sikaneta I. Optimum GMTI processing for spacebased SAR/GMTI systems-Simulation results[C]. 2010 8th European Conference on Synthetic Aperture Radar, Aachen, Germany, 2010: 653–656.

[14]Cerutti-Maori D, Gierull C H, and Ender J H G. Optimum SAR/GMTI processing and its application to the radar satellite RADARSAT-2 for traffic monitoring[J].IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(10): 3868–3881.

邢孟道(1975–)，男，浙江嵊州人，西安电子科技大学教授，博士生导师，主要研究方向为雷达成像、目标识别和天波超视距雷达信号处理。

E-mail: xmd@xidian.edu.cn

孙光才(1984–)，男，湖北孝感汉川人，西安电子科技大学讲师，主要研究方向为多通道波束指向SAR成像和SAR动目标成像。

E-mail: rsandsgc@126.com

李学仕(1988–)，男，黑龙江哈尔滨人，西安电子科技大学信号与信息处理专业博士研究生，主要研究方向为SAR成像及SAR动目标成像。

E-mail: lixueshi@stu.xidian.edu.cn

Study on SAR/GMTI Processing for High-resolution Wide-swath SAR System

Xing Meng-dao Sun Guang-cai Li Xue-shi
(National Laboratory of Radar Signal Processing,Xidian University,Xi’an710071,China)

In recent years, there have been many studies on Synthetic Aperture Radar (SAR), and High-Resolution and Wide-Swath (HRWS) imaging has attracted increasing attention. High-resolution images of a stationary scene enabling the detection of moving targets can be provided by an SAR system combined with Ground Moving Target Indication (GMTI), which is widely used in many civilian applications. A multichannel in azimuth HRWS SAR system has been presented to effectively deal with conflicts between high resolution and low Pulse Repetition Frequency (PRF), which are encountered during HRWS SAR imaging. Because the spatial degrees of freedom in the azimuth can be employed to suppress a clutter, the multichannel configuration offers significant potential for GMTI capability. In this study, a new clutter suppression and moving target imaging method is proposed, which can be operated in an HRWS SAR system without the need for a separate high PRF operation mode.

Synthetic Aperture Radar (SAR); High-Resolution Wide-Swath (HRWS); Moving targets detection; Moving targets imaging; Simultaneous

TN957

A

2095-283X(2015)-04-0375-11

10.12000/JR15096

邢孟道, 孙光才, 李学仕. 用于高分辨率宽测绘带SAR系统的SAR/GMTI处理方法研究[J]. 雷达学报, 2015, 4(4): 375–385.

10.12000/JR15096.

Reference format:Xing Meng-dao, Sun Guang-cai, and Li Xue-shi. Study on SAR/GMTI processing for highresolution wide-swath SAR system[J].Journal of Radars, 2015, 4(4): 375–385. DOI: 10.12000/JR15096.

2015-08-15；

2015-09-10

*通信作者： 李学仕 lixueshi@stu.xidian.edu.cn

国家自然科学基金(61222108, 61301292, 61101245)和中央高校基本科研业务费专项资金(K5051302058, K5051302046)

Foundation Items: The National Natural Science Foundation of China (61222108, 61301292, 61101245); Fundamental Research Funds for the Central Universities (K5051302058, K5051302046)