离包络和方位相位2维空变，严重影响成像质量[3,4]。因此，对于俯冲段大斜视SAR子孔径成像两维空变校正方法的研究具有重要的意义。针对俯冲SAR大斜视成像2维空变校正，文献[5]提出一种曲线轨迹高分辨俯冲大斜视波数域成像算法，但是成像中需要复杂的插值操作，时间复杂度高，不适于实时处理。文献[6]提出通过方位分块处理来实现沿方位空变的RCM和多普勒空变参数校正，得到曲线轨迹俯冲模式下的2维图像。但是在非线性轨迹下，方位分块可能会导致数据的不连续，在图像拼接中出现误差；同时分块处理也会增加运算量，不利于实时实现。文献[7]提出基于方位非线性变标(Nonlinear Chirp Scaling, NCS)来校正多普勒调频率的空变，但未考虑到平台存在3维加速度的情况。文献[8]也采用NCS来补偿空变的多普勒调频率，虽然考虑了加速度，但未考虑3次项系数的空变性，在方位幅宽较大的情况下，残留的3次空变相位限制了方位聚焦深度。对于高速机动平台SAR成像，采用子孔径处理更具优势。针对俯冲SAR大斜视子孔径成像问题，本文首先建立俯冲段瞬时斜距模型，随后分析斜视角沿距离向的空变和3维加速度对2维频谱的影响；然后通过预处理来降低距离方位的耦合及补偿加速度对2维频谱的影响，在此基础上提出基于2级频域滤波(Two-stage Frequency Filtering, Ts FF)的2维空变校正新方法以改善方位聚焦深度，适用于大斜视宽幅成像；最后通过点目标仿真验证了本文方法的有效性。2    俯冲大斜视回波信号分析2.1  俯冲段瞬时斜距模型OXnYeZsg= (vn; ve; vs)g= (an; ae; as)vs°= (vx; 0; vz)= (ax; ay; az)tm= 0对SAR运动平台，通过惯导系统可得到其在地理坐标系(即图1中表示北天东方向的 坐标系 ) 中 的 3 维 速 度 和 3 维 加 速 度。在传统的SAR成像中，平台沿北向或偏离北向一定角度的方向水平飞行，天向速度 为零，此时数据录取平面为速度方向和波束射线在零多普勒面内的投影组成的平面。但在俯冲成像中，由于3维速度的存在，平台沿偏离正北方向的方向做俯冲运动；由于3维加速度的存在，平台偏离理想运动轨迹作曲线运动，如图1所示。此时，数据录取平面发生旋转。图1中，假设在子孔径时间内平台从M点以初速度 和加速度 运动到N点，Q点为方位慢时间 时刻平台所在位置，原点O为Q点在地面的投影点。以XOY平面为成像平面，平台位于Q点P1D1xnX 0Y 0Z时，波束中心射线(Line Of Sight, LOS)与地面交点为P, D为场景中与P有相同X轴坐标E的任意一点， 和 为点P和D在零多普勒面内的投影点，T为场景中与D点沿X轴相距 的任意一点。俯冲成像时，数据录取平面为速度方向和波束中心射线组成的平面，此时的成像坐标系可等效为图1中的坐标系，斜视角可等效为速度方向和波束中心射线方向夹角的余角。这种成像坐标系的旋转使传统的成像方法不再实用，因此需要研究基于该模型下的成像算法。µ®RsRdµs µd H0R定义斜距 为平台运动至MN之间任意位置时到T点的斜距， 为对应的斜视角(LOS与YOZ平面的夹角)，H为雷达平台高度， 为方位侧视角(LOS在XOY平面的投影与X轴正向的夹角)。假定平台位于Q点时到P点和D点的斜距分别为 和 ，对应的斜视角分别为 和 ，平台高度为 。以任意一点C为例，平台在该点到T点的斜距 可用泰勒级数展开到4阶项，即R (tm; Rd; xn) = ¯¯¯¡¡!¯¯¯=4Xj=0kj(Rd; xn) tjm(1)kj(Rd; xn)=(1=j!)¡djR (tm; Rd; xn) =dtjm¢¯¯tm=0tm= 0 k1k2k3k4其中，表示斜距在 处的第j 项展开式系数。其中为线性距离徙动项(Linear RCM, LRCM),  为QRCM项， 和 为高阶距离徙动项(High-orderRCM, HRCM)。假设发射信号为线性调频信号(Linear Fre-quency Modulation, LFM)信号，将T点的回波信号解调后变换到距离频域为Ss (fr; tm; Rd; xn)= Wr(fr) wa(tm) exp µ¡jfr2Kr¶¢exp µ¡j4c(fc+ fr) R (tm; Rd; xn)¶ (2)c Krfcfr其中， 为光速， 为调频率， 和 分别表示雷  图 1 俯冲SAR大斜视成像几何模型2622 电  子  与  信  息  学  报 第 40 卷Wr(¢)wa(¢)达中心载频和距离频率， 为距离窗的频域形式， 方位窗函数。式(2)中第1个指数项表示距离调制项，第2个指数项表示距离方位耦合项。2.2  距离向空变的斜视角对距离包络的影响分析在图1中，定义擦地角¯为LOS与XOY平面的夹角，从而平台位于C点时到成像平面中任意一点的斜视角可表示为sin µ = cos ¯ cos ® =pR2¡H 2=R¢cos ® (3)在传统的大斜视成像中，通常先以场景中心为参考进行走动校正来降低2维耦合[9]。但由于斜视角沿距离向的空变，导致走动校正后存在线性残留量，这个残留量对距离向包络会带来很大的影响，而现有文献常忽略了这部分的影响[10]。下面对残余线性项对距离包络空变性影响进行分析，仿真参数如表1所示，线性残余量对距离包络的影响如图2所示。从图2中可看出，在距离向单边幅宽大于210 m的情况下，线性残余量会超过1/2距离分辨率，根据单边包络误差不可忽略的准则，此时的距离向幅宽不能大于420 m。此外，2阶项的残余量远小于分辨率，在后续处理中可将其忽略。因此，本文将在考虑斜视角线性空变的基础上研究俯冲SAR大斜视2维成像算法。2.3  3维加速度对2维频谱影响分析kj(Rd; xn) =kj( ; Rd; xn)¡kj( ; Rd)下面分析加速度的存在对2维频谱的影响。可将式(1)中各阶项系数分离成两部分：，其中第1部分可看作平台匀速运动时轨迹斜距历程，第2部分反映了加速度对斜距的影响。加速度不会影响线性项，但会直接影响高次项系数，尤其是方位带宽，具体地Ba= 4Ta¢ jk2( ; Rd; xn)¡k2( ; Rd)j=Kr (4)Tak2( ;Rd) k2( ; Rd; xn)其中， 为子孔径时间。从式(4)可以看出，若和 取值相反，会导致方位带宽超过采样率，造成方位频谱混叠；若二者取值接近，会导致方位频谱位于零频处[11]。因此，为了保证后续处理的有效性，需要考虑到加速度对频域的影响。3    俯冲SAR子孔径两级频域滤波成像算法在分析了斜视角沿距离向的空变对距离包络的影响和3维加速度对2维频谱的影响的基础上，本节研究了俯冲SAR子孔径成像算法。首先通过预处理进行走动校正及补偿加速度对2维频谱的影响；然后采用第1级频域滤波(First-stage FrequencyFiltering, Fs FF)校正距离包络沿方位向的空变来实现RCM统一校正；通过Fs FF可消除大部分RCM的空变，最后再采用第2级频域滤波(Second-stageFrequency Filtering, Ss FF)校正方位相位的空变来实现统一聚焦。3.1  预处理在俯冲大斜视SAR成像中，较大的线性走动造成了距离向和方位向严重的耦合，同时加速度的存在会影响2维频谱。首先以场景中心点P点为参考点构造走动校正函数和加速度补偿函数，可表示为Hp(fr; tm) = exp µj4c(fc+ fr)hk1(Rs; 0) tm+4Xi=2ki( ; Rs) timi!(5)k1(Rs其中，; 0)表示参考点P点的线性项系数。Ss (fr; tm; Rd; xn将式(5)和)相乘，可得预处理后的信号相位表达式为Á (fr; tm; Rd; xn)=¡j4c(fc+ fr)"A0(Rd; xn)+A1(Rd; xn) tm+4Xj=2Aj(Rd; xn) tjm#(6)A0(Rd; xn) = k0(Rd; xn) A1(Rd; xn) = k1(Rd; xn)¡k1(Rs; )Aj(Rd; xn)=kj(Rd; xn; )¡kj(Rs; ) ; j =2;其中， 为常数项系数，为残余线性项系数，表 1  仿真参数载频 15 GHz 脉冲重复频率 2.5 k Hz距离带宽 200 MHz 高度 5 km脉冲宽度 10 μs 场景中心斜距 10 km方位侧视角 28° 平台速度 (145, –38, –35) m/s加速度 (1.2, 0.5, –0.8) m/s2合成孔径时间 0.6 s  图 2 距离幅宽与RCM残余量关系曲线图第 1 1 期 党彦锋等：俯冲段大斜视SAR子孔径成像二维空变校正方法 2623xn= 0为弯曲项和高阶项系数。为方便后续处理，将式(7)中各阶项系数在 处泰勒级数展开，可得Ai(Rd; xn) = Ai(Rd; 0) +rAi(Rd; xn) xn= Ai0+ Ai1xn; i = 1; 2; 3 (7)由于4阶项系数的空变性很小，可将其忽略，A4(Rd; xn) = A40即 。下面采用Fs FF抑制距离包络的空变来实现RCM统一校正。3.2  基于Fs FF的距离包络空变校正采用级数反演法(Method of Series Reverse,MSR)将信号变换到2维频域[12]，其相位表达式为(fr; fa; Rd; xn) =¡4 A0(Rd; xn)c(fc+ fr) +c4A2(Rd; xn) (fc+ fr)·fa+2A1(Rd; xn)c(fc+ fr)¸2+c2A3(Rd; xn)16A23(Rd; xn) (fc+ fr)2·fa+2A1(Rd; xn)c(fc+ fr)¸3+c3£9A32(Rd; xn)¡4A2(Rd; xn) A4(Rd; xn)¤256A25(Rd; xn) (fc+ fr)3·fa+2A1(Rd; xn)c(fc+ fr)¸4(8)为实现距离包络空变校正，同时防止引入多余的相位使后续的方位向处理复杂化，本文引入Fs FF因子，可表示为HF(fr; fa) = exp µj µp1c22 (fc+ fr)3frf3a+ p2c34 (fc+ fr)4frf4a¶¶ (9)p1p2其中， 和 为待定量。将式(8)与式(9)中的相位部分相加，并变换到方位时域，将方位时域相位信号整理后可表示为Á (fr; ta; Rd; xn)=¡4 fcc4Xi=0Ai(Rd; xn) tim¡4 frc4Xi=0Bi(Rd; xn; p1; p2) tim(10)Bi(Rd; xn; p1; p2)xntmxnt2mx2nt2mC1(p1; p2; Rd) = 0 C2(p1; p2; Rd) = 0其 中 ， 第 1 项 为 空 变 的 方 位 相 位 ， 第 2 项 中为各阶RCM项系数。为消除包络空变，将第2项中各阶项系数按照式(7)展开，可得一系列关于 和 的各阶耦合量。一般地，和 的存在对包络空变性影响最大。令这两项的系数满足 和 ，联立方程式，可得p1=¡¡A112A202+ 4A10A11A20A21+ A102A212¢12A10A207A113³3A10A11A21+ 2A112A20´p2=¡¡A10A212+ 2A11A20A21¢24A10A207A113³3A10A11A21+ 2A112A20´9>>>>>>=>>>>>>;(11)fr= 0将式(11)代入式(9)，然后将式(9)中的相位部分与式(8)相加，就可在2维频域消除包络的空变性。然后将频域信号在 处泰勒级数展开，可得(fr; fa) = j [Áaz(fa; Rd; xn) + 'rcc (f a; Rd) f r+'src (f a; Rd) f2r¤(12)其中，Áaz (f a; Rd; xn)为空变的方位相位，' rcc (f a; Rd)和'src (f a; Rd)分别为距离弯曲项和2次脉冲压缩项相位。由于其不存在方位空变，可以统一校正，以场景中心点P点为参考RCC和SRC因子，可表示为Hc(fr; fa) = exp [¡j ('rcc(fa; Rs) fr+'src (f a; Rs) f2r¢¤(13)将式(13)中相位部分与式(12)相加就可实现距离包络空变校正。经Fs FF处理已经极大地弱化了多普勒参数随方位位置的空变，但是由于相位对空变的敏感性高于包络，尤其是在高频段。为保证聚焦质量，本文提出Ss FF来校正方位相位的空变以实现方位统一聚焦。3.3  基于Ss FF的方位相位空变校正Rd0 = Rd+ xnsin µ在时域走动校正后，距离单元相同方位位置不同的点目标将出现在不同的距离单元内，且其位置满足 ，这导致多普勒调频率和方位向高次项的空变，方位向平移不变性不再满足[2]。根据之前的分析，距离向处理后的方位相位可表示为Á (tm; Rd0; xn) =¡4¸4Xi=0Ai(Rd0; xn) tim(14)A1(Rd0; xn)=4可以看出式(14)中各阶方位相位系数均具有方位空变性。特别地， 表示残留多普勒中心项，该项对于小斜视窄波束成像影响较小，但在大斜视宽幅成像中(残留多普勒中心项远大于 )，它的存在将会严重影响方位聚焦深度，因此在后续处理中需要考虑。具体地，各阶空变的方位相位误差为¢ m =[ m(Rd0 ; xn)¡ m(Rd0 ; 0)] [(¢fa+ Bsub=2)m¡(¢fa¡Bsub=2)m] ; m = 2; 3; 4 (15)m(Rd0; xn) Bsub其中， 为各阶相位频域表达式， 为2624 电  子  与  信  息  学  报 第 40 卷俯冲段大斜视SAR子孔径成像二维空变校正方法党彦锋①②      梁   毅*①②      别博文①②      丁金闪①②      张玉洪③①(西安电子科技大学雷达信号处理国家重点实验室   西安   710071)②(西安电子科技大学信息感知技术协同创新中心   西安   710071)③(西安电子科技大学电子工程学院   西安   710071)摘   要：俯冲合成孔径雷达(SAR)常采用大斜视子孔径成像来满足平台机动性和处理实时性；而大斜视模式下距离方位的严重耦合、斜视角沿距离向的空变及3维速度和加速度的存在会导致距离包络和方位相位的2维空变，严重影响成像质量。针对这些问题，该文提出一种新的两级频域滤波算法，在预处理后，首先引入第1级频域滤波因子以校正距离包络沿方位位置的空变以实现距离徙动统一校正，然后采用第2级频域滤波校正方位多普勒参数的空变实现方位相位统一聚焦。仿真实验验证了该算法的有效性。关键词：合成孔径雷达；子孔径；俯冲大斜视；两级频域滤波算法；2维空变校正中图分类号：TN957.52 文献标识码：A 文章编号：1009-5896(2018)11-2621-09DOI: 10.11999/JEIT180021Two-dimension Space-variance Correction Approach for Diving HighlySquinted SAR Imaging with Sub-apertureDANG Yanfeng①②      LIANG Yi①②      BIE Bowen①②     DING Jinshan①②      ZHANG Yuhong③①(National Key Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China)②(Collaborative Innovation of Information Sensing and Understanding, Xidian University, Xi’an 710071, China)③(School of Electronic Engineering, Xidian University, Xi’an 710071, China)Abstract: The diving SAR usually adopts the highly squinted mode and sub-aperture to satisfy themaneuvering and real-time processing. However, the existence of severe range-azimuth coupling, range-dependent squint angle and three-dimension velocity and acceleration leads to the space variance of rangeenvelope and azimuth phase, which makes imagery unfocused seriously. To solve these problems, a Two-stageFrequency Filtering Algorithm (Ts FFA) is proposed. After preprocessing, the First-stage Frequency Filtering(Fs FF) factor is first introduced to correct azimuth-dependent Range Cell Migration (RCM) and realize theunified RCM correction. Furthermore, the Second-stage Frequency Filtering (Ss FF) factor is adopted toequalize azimuth-dependent Doppler parameters and realize unified azimuth phase focused. Simulation resultsare presented to validate the effectiveness of the proposed approach.Key words: SAR; Sub-aperture; Diving highly squinted; Two-stage Frequency Filtering Algorithm (Ts FFA);Two-Dimension (2-D) space-variance correction1    引言将合成孔径雷达(SAR)与高速平台结合并运用于制导技术已成为当前研究热点[1]。为保证平台机动性，需以大斜视模式成像；为保证观测实时性，需采用子孔径进行相干处理，以简化处理流程，减小存储量和计算量，实现快视成像[2]。理想情况下的SAR平台是沿航迹方向匀速直线飞行，垂直方向无速度分量。但高速无人机SAR等高速机动平台具有灵活多样的飞行轨迹，如地形回避、战术规避及向下俯冲等，此时回波特性要比一般的机载SAR复杂。本文就俯冲状态下的SAR信号进行研究，在俯冲大斜视模式下，垂向速度和3维加速度的存在导致距离方位严重耦合，传统的方位平移不变性假设不再成立；同时存在斜视角沿距离向的空变，这  收稿日期：2018-01-05；改回日期：2018-06-13；网络出版：2018-07-20*通信作者： 梁毅　yliang@xidian.edu.cn基金项目：海洋公益性行业科研专项经费(201405001)Foundation Item: The Public Science and Technology ResearchFunds Projects of Ocean (201405001)¢fa子孔径成像时的方位向多普勒带宽， 为任意方位位置处点目标的多普勒频率与多普勒中心频率之差。以第2节中参数为例，可得方位相位误差如图3所示。可以看出，2阶相位和3阶相位误差随方位位=4=4置的增大而逐渐大于 ，而4阶相位误差小于，因此进行方位相位空变校正时需考虑残留多普勒中心项、2阶和3阶相位空变，忽略4次相位的空变。jk (xn)j<jk (x0)j<jk (x¡n)jTsub同样地，在具体公式推导前，首先用图4所示的时频分布线对方位相位补偿过程进行说明。在采用子孔径成像时，场景中各点数据支撑区在时域是重叠的，在频域是分开的；如图4(a)所示，假设场景中心和两侧等距处各有一个点目标，由于各点方位位置不同，对应的相位曲率也不同，且满足，图中 为子孔径录取时间，亦等于时频分布线在时间轴的投影。通过分析发现，首先进行高次非空变相位补偿可缩短数据支撑区的宽度，进而减少补零点数(后续会进行分析)，提高实时性；因此，构造如图4(b)所示的高阶项补偿函数，在方位时域对高次相位进行补偿。对补偿后的经FFT后，相频分布线具有不同的T0subjk (x0)jTsub ext曲率(为时域曲率的倒数)，且其在频率轴上的投影是错开的，如图4(c)所示，图中 为补偿高阶非空变相位后的数据支撑区(小于原数据支撑区)。为了实现方位相位统一聚焦，在方位多普勒域引入Ss FF因子(如图4(d)所示)，该因子可将相频分布线校正至相同的曲率。如图4(e)所示，将校正后的信号经IFFT后，不同点的时频分布线具有相同的斜率 ，但此时时频分布线在时间轴的投影超出原数据支撑区，这种位置的偏移会在聚焦是产生虚假点，因此需要将支撑区扩展来消除这些虚假点，图中 为扩展后的数据支撑区；最后构造De-ramp因子实现方位统一聚焦。下面进行具体的公式推导。  图 3 方位空变相位误差随方位位置变化曲线  图 4 方位相位空变统一校正示意图第 1 1 期 党彦锋等：俯冲段大斜视SAR子孔径成像二维空变校正方法 2625首先进行高阶非空变相位项补偿，构造高阶项补偿函数为Hh(tm; Rd0 )= exp·j4¸ ¡A30(Rd0 ) t3m+A40(Rd0 ) t4m¢¸(16)将补偿完高阶非空变相位的信号经FFT到方位频域，根据之前分析，要消除调频率和3阶项系数的空变，需要构建Ss FF因子，表示为HS(fa) = exp£j¡q1f3a+ q2f4a¢¤(17)q1q2其中 和 为未知量。将式(17)与方位频域信号相乘，然后经IFFT到方位时域，得ss (tr; tm; Rd0; xn) = sinc ·Brµtr¡2Rd0c¶¸ wa(tm)¢exp (¡jÁ (tm; Rd0; xn)) (18)trBrÁ(tm; Rd0 ;xn)xntmxnt2mx2nt2mD1(q1; q2; Rd0 ) = 0 D2(q1;q2; Rd0 ) = 0其中， 为距离快时间， 为距离向带宽，为方位时域相位。将方位相位按照式(7)展开，同样可得一系列关于 和 的各阶耦合量，由于和 的存在对方位相位的空变性影响最大[11]。令这两项的系数满足 和，联立方程组，可得q1=¡2C213C202C11; q2=¡4C2123C112C203(19)Ai(Rd0 ; 0) +rAi(Rd0 ; xn) xn= Ci0+ Ci1xn(i = 1; 2)其中，，将式(19)代入式(17)，然后与信号方位频域表达式相乘就可消除方位相位的空变性。最后构造方位统一聚焦因子实现方位相位空变校正，具体地可表示为Hd(tm; Rd0 ) = exp£j A¡t2m; t3m; t4m; q1; q2¢¤(20)A¡t2m; t3m; t4m; q1; q2¢其中， 表示与方位位置无关，可统一补偿的方位调制项。从而聚焦后的图像为s S (tr; fa; Rd0) = sinc ·Brµtr¡2Rd0c¶¸¢sinc ·1Bsub¡fa¡b1xn¡b2x2n¢¸(21)b1b2其中， 反映点目标的方位位置， 会造成目标方位位置的偏差，需要进行几何校正来得到正确的点目标位置。整个处理流程如图5所示。Nr£Na5NrNalog2£Nr=a¤6NrNa下面分析本文方法的浮点运算量。假设距离与方位采样点数为 。对于每次的FFT/IFFT运算来说，其浮点运算量为 ；复乘运算量为 。从流程图5可知，所提方法共有2次距离向的FFT/IFFT操作，5次方位向的FFT/IFFT操作，5次复乘。因此，总的运算量为C = NrNa£[10£log2Nr+ 25£log2Na+ 30] (22)4    仿真数据结果及分析为验证本文所提2维空变校正方法的有效性，这里通过仿真数据来说明，仿真参数如表1所示。在地面场景中沿雷达波束视线方向和垂直于波束视线方向放置一大小为1 km×1 km的3×3矩形点阵，如图6所示，图中点1和点5为方位空变最剧烈的点，点2和点4位于同一距离单元，以场景中心点3为参考点。首先验证Fs FF校正包络空变的有效性，选取点2、点3和点4的校正结果进行分析。图7(a1)—7(a3)为 采 用 参 考 算 法 [1 1]进 行 包 络 空 变 校 正 结 果 ，图7(b1)—7(b3)为本文算法校正结果。可以看出，采用参考算法校正后的边界点存在包络线越距离单元的现象，会直接影响最终的方位聚焦效果；而采用本文方法进行统一RCM校正后，包络线不存在越距离单元的现象，验证了本文方法的有效性。然后在Fs FF的基础上验证Ss FF校正相位空变  图 5 俯冲大斜视子孔径成像算法流程图  图 6 点目标位置示意图2626 电  子  与  信  息  学  报 第 40 卷、点3和点5的校正结果进行分析。为了更清晰地反映边缘点成像效果，图8(a1)—8(a3)给出了采用参考算法[8]方位聚焦结果剖面图，可以看出，对于边缘点1和5，由于忽略了残留多普勒中心和3次相位系数随方位位置的空变性，造成第1旁瓣的抬高，进而使峰值旁瓣比损失，最终影响方位分辨率。图8(b1)—8(b3)给出了采用本文方法方位聚焦结果剖面图，可以看出，在考虑了残留多普勒中心项和3次相位系数随方位位置的空变的基础上，可得到聚焦效果良好的成像结果。进一步地，计算点1、点3和点5的峰值旁瓣比、积分旁瓣比和方位分辨率，结果如表2所示。可以看出参考算法所得指标相比理论值(峰值旁瓣比–13.26 d B，积分旁瓣比–9.8 d B，方位分辨率1 m)偏差较大，而本文方法所得指标与理论指标基本吻合，从而进一步证明本文方法有效性。5    结束语本文针对俯冲段大斜视SAR子孔径成像中存在的2维空变问题，提出一种两级频域滤波校正空变的方法。根据平台的实际运动状态，建立了俯冲大斜视SAR成像几何模型，并分析了斜视角沿距离向的空变对距离包络的影响及3维加速度对2维频谱的影响。在预处理实现线性走动校正和加速度补偿后，针对距离包络空变，在2维频域引入第1级滤波因子来校正其空变性，从而实现距离徙动统一校正；针对方位相位空变，考虑了残留多普勒中心项，并详细分析了多普勒调频率和3次项系数的空变性，然后根据子孔径数据特征在方位频域引入第2级滤波因子来校正其空变性，最后基于谱分析实现方位统一聚焦。整个算法只包含FFT和复乘运算，不涉及插值处理，易于工程实时实现。仿真数据处理验证了本文方法的有效性。表 2  性能指标测量结果(未加窗)性能指标参考方法 本文方法点1 点3 点5 点1 点3 点5峰值旁瓣比(d B) –3.46 –13.23 –5.63 –13.23 –13.25 –13.21积分旁瓣比(d B) –4.24 –9.77 –5.92 –9.77 –9.83 –9.76方位分辨率(m) 1.58 1.05 1.46 1.03 1.01 1.04  图 7 包络空变校正结果图第 1 1 期 党彦锋等：俯冲段大斜视SAR子孔径成像二维空变校正方法 2627参 考 文 献保铮, 邢孟道, 王彤. 雷达成像技术[M]. 北京: 电子工业出版社, 2005: 123–182.[1]李震宇, 梁毅, 邢孟道, 等. 弹载合成孔径雷达大斜视子孔径频域相位滤波成像算法[J].  电子与信息学报,  2015,  37(4):953–960. doi: 10.11999/JEIT140618.LI Zhenyu, LIANG Yi, XING Mengdao, et al. A frequencyphase  filtering  imaging  algorithm  for  highly  squintedmissile-borne Synthetic Aperture Radar with subaperture[J].Journal of Electronics & Information Technology, 2015,37(4): 953–960. doi: 10.11999/JEIT140618.[2]LI  Dong,  LIN  Huan,  LIU  Hongqing,  et al.  Focusimprovement  for  high-resolution  highly  squinted  SARimaging based on 2-D spatial-variant linear and quadraticRCMs  correction  and  azimuth-dependent  Dopplerequalization[J]. IEEE Journal of Selected Topics in AppliedEarth Observations and Remote Sensing,  2017,  10(1):168–183. doi: 10.1109/JSTARS.2016.2569561.[3]HUAI Yuanyuan, LIANG Yi, DING Jinshan, et al. Aninverse extended Omega-K algorithm for SAR raw datasimulation with trajectory deviations[J]. IEEE Geoscienceand Remote Sensing Letters, 2016, 13(6): 826–830. doi:10.1109/LGRS.2016.2548240.[4]LI Zhenyu, XING Mengdao, XIONG Tao, et al. A modifiedequivalent range model and wavenumber-domain imaging[5]approach for high-resolution high-squint SAR with curvedtrajectory[J]. IEEE Transactions on Geoscience and RemoteSensing,  2017,  55(7):  3271–3734.  doi:  10.1109/TGRS.2017.2678763.ZENG Tao, LI Yinghe, DING Zegang, et al. Subapertureapproach based on azimuth-dependent range cell migrationcorrection and azimuth focusing parameter equalization formaneuvering high-squint-mode SAR[J]. IEEE Transactionson Geoscience and Remote Sensing,  2015,  52(12):6718–6734. doi: 10.1109/TGRS.2015.2447393.[6]肖忠源,  徐华平,  李春生.  基于俯冲模型的频域距离走动校NLCS-SAR成像算法[J].  电子与信息学报,  2013,  35(5):1090–1096. doi: 10.3724/SP.J.1146.2012.01207.XIAO Zhongyuan, XU Huaping, and LI Chunsheng. NLCSSAR  imaging  algorithm  with  range-walk  correction  infrequency  domain  based  on  dive  model[J].  Journal ofElectronics &  Information Technology,  2013,  35(5):1090–1096. doi: 10.3724/SP.J.1146.2012.01207.[7]周松, 包敏, 周鹏, 等. 基于方位非线性变标的弹载SAR下降段成像算法[J]. 电子与信息学报, 2011, 33(6): 1420–1426. doi:10.3724/SP.J.1146.2010.01124.ZHOU Song, BAO Min, ZHOU Peng, et al. An imagingalgorithm for missile-borne SAR with downward movementbased on azimuth nonlinear chirp scaling[J]. Journal ofElectronics &  Information Technology,  2011,  33(6):[8]  图 8 点目标方位聚焦结果剖面图2628 电  子  与  信  息  学  报 第 40 卷