基于成像坐标系优化的中轨星载SAR成像方法李 航 刘文康* 孙光才* 邢孟道 李光伟 费晓燕①(西安电子科技大学雷达信号处理国家重点实验室 西安 710071)②(西安电子科技大学信息感知技术协同创新中心 西安 710071)摘 要：在中轨合成孔径雷达(MEO SAR)成像中，大弯曲轨道以及长合成孔径时间会导致信号产生严重的两维空变。常规方法分别在距离和方位两个方向处理空变，计算复杂度通常比较高。该文研究了大场景中的多普勒调频率的空间分布，并提出将数据变换到一种非正交非线性成像坐标系中进行成像，使中轨SAR信号在该坐标系中满足方位平移不变性，由于不需要对方位空变做额外处理，该成像方法的运算量显著降低。最后通过多普勒线性化处理可以进一步补偿高阶多普勒参数的影响，以实现场景边缘点更精确的聚焦，并校正由非线性坐标系变换引入的方位聚焦位置偏移。最后，在条带模式下仿真2 m分辨率的数据，可以验证所提出算法的有效性。关键词：中轨SAR；两维空变；成像坐标系优化；非正交非线性成像坐标系中图分类号：TN957.52 文献标识码：A 文章编号：2095-283X(2020)05-0856-09DOI: 10.12000/JR20098引用格式：李航, 刘文康, 孙光才, 等. 基于成像坐标系优化的中轨星载SAR成像方法[J]. 雷达学报, 2020, 9(5):856–864. doi: 10.12000/JR20098.Reference format: LI Hang, LIU Wenkang, SUN Guangcai, et al. MEO SAR imaging based on imagingcoordinate system optimization[J]. Journal of Radars, 2020, 9(5): 856–864. doi: 10.12000/JR20098.MEO SAR Imaging Based on Imaging Coordinate System OptimizationLI Hang LIU Wenkang* SUN Guangcai* XING MengdaoLI Guangwei FEI Xiaoyan①(National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China)②(Collaborative Innovation Center of Information Sensing and Understanding,Xidian University, Xi’an 710071, China)Abstract: In the Medium-Earth-Orbit Synthetic Aperture Radar (MEO SAR), the curved trajectory and longsynthetic aperture time lead to a two-dimensional spatial variation in the signals. Traditional methods usuallyprocess the range and azimuth variations separately, and the computational complexities are high. Herein, westudy the Doppler rate distribution across a large scene and propose a non-orthogonal and nonlinear coordinatesystem wherein the MEO SAR signals satisfy the azimuth-shift–invariant property. Thus, the efficiency of theimage formation processor can be significantly improved. The higher-order Doppler parameters are addressedby the Doppler linearization. Then, more precise focusing can be achieved, and the azimuth time-shift causedby the changes in signal distribution is addressed. Finally, the processing results of simulated stripmap-modedata with a 2-m resolution are presented to validate the effectiveness of the proposed algorithm.Key words: Medium-Earth-Orbit Synthetic Aperture Radar (MEO SAR); 2-D space variance; Imagingcoordinate system optimize; Non-orthogonal and nonlinear coordinate system收稿日期：2020-07-08；改回日期：2020-09-15；网络出版：2020-10-09*通信作者： 刘文康 wkliu@stu.xidian.edu.cn； 孙光才 rsandsgc@126.com*Corresponding Author: LIU Wenkang, wkliu@stu.xidian.edu.cn; SUN Guangcai, rsandsgc@126.com基金项目：国家自然科学基金重点项目(61931025)，高等学校学科创新引智计划资助(B18039)Foundation Items: The State Key Program of National Natural Science China (61931025), The 111 Project (B18039)责任主编：李宁 Corresponding Editor: LI Ning第9 卷第5 期雷 达 学 报Vol. 9No. 52 0 2 0年1 0月Journal of Radars Oct. 20201 引言与目前的低地球轨道星载合成孔径雷达(Low-Earth-Orbit Synthetic Aperture Radar, LEOSAR)相比，中地球轨道星载合成孔径雷达(Medium-Earth-Orbit Synthetic Aperture Radar, MEOSAR)具有更大的幅宽；它的广域覆盖范围长达数百公里，重访时间短，只有几分钟，在许多领域中显示出巨大应用潜力。MEO SAR长时间观测目标区域，提供连续的SAR图像，有助于实现动态监视。利用多级带宽设计和波束扫描，MEO SAR可以工作在不同的模式，并实现灵活的分辨率和场景覆盖。由于MEO SAR所在的中轨卫星与地球之间的相对运动较慢，在满足分辨率要求的情况下，合成孔径时间变得更长。而随着合成孔径时间的增加会出现一些问题，例如大气相位误差、斑点噪声等[1–3]。中地球轨道的轨迹曲率也可能影响波束驻留时间和等效卫星速度，从而导致分辨率和信号特性随轨道时间改变，因此波束指向控制对于中高轨曲线轨迹SAR具有重要的意义——能够调整和保持系统标称分辨率，波束扫描速率需要根据轨道位置、视角等进行设计。另外，对于大场景成像，MEO SAR信号方位平移不变性失效，此时回波信号是两维空变的[4,5]。对于MEO SAR信号的两维空变问题，常规的成像算法将其分两维处理。一些成熟的距离空变校正方法，例如chirp scaling[6–8]和omega-K[9,10]算法，经过改进后，同样能用来有效处理弯曲轨迹的SAR数据。经过方位非线性变标[11–13]，重采样[14,15]等一些方位向预处理方法处理之后，信号能够满足方位平移不变性，但是这些方位空变校正处理会引入多次傅里叶变换或者插值操作，明显增加了成像算法的运算量。一些时域算法也适用于弯曲轨迹的星载SAR成像，例如BPAs[16,17]和FFBPAs[18–20]，这些基于时域插值的算法能够直接处理复杂的两维空变，是未来高分辨率、弯曲轨道星载SAR处理的一个潜在选择。通常，尽管频域算法需要一些额外的操作来获得方位平移不变性，但它们的计算复杂度通常小于时域算法。本文为条带模式下的MEO SAR数据的两维空变处理提出了一种有效的处理方法。基于MEOSAR大场景中的多普勒调频率分布，通过非正交非线性坐标系变换将两维空变集中在一维。在推导出的非正交非线性坐标系中，新的方位向对应于多普勒调频率等高线，因此，在该方向上的多普勒调频率空变可以最小化，在获得方位平移不变性的同时避免复杂的方位空变校正处理。然后，直接应用经典星载SAR成像算法来完成信号聚焦。为了实现更精确的方位压缩，采用多普勒重采样以抵消高阶多普勒参数的空变影响，并且消除由非线性坐标系变换引起的方位聚焦偏移。2 信号建模ST 

 a如图1描述了MEO SAR成像几何。表示卫星位置， 表示照射区域的目标， 表示仰角， 表示斜视角， 表示入射角，在数据录取期间，轨道弯曲明显，MEO卫星与场景中目标的相对斜距历程很复杂。在本文中，相对斜距历程使用4阶泰勒展开模型表示，将斜距在多普勒中心时刻扩展为R(ta) = ∥ST∥  R0 +Σ4i=1qi(ta � tc)i (1)ST S T R0tatc qi; i = 1; 2; 3; 4ifdc = �2q1/ fdr = �4q2/其中， 是从卫星到目标的向量， 表示波束中心时刻雷达到目标的中心斜距， 表示方位慢时间， 表示目标波束中心时刻， 表示第项的多普勒系数。多普勒中心与线性分量直接相关，多普勒中心表示为： , 是电磁波波长。多普勒调频率可以用二次项系数表示为，为了实现方位向聚焦，需要准确计算二次项。用将距离模型的2阶系数展开为方位时间与距离的函数q2 = q20+q21;rr+q22;rrr2+q21;tc tc+q22;tctc t2c +q22;rtcrtc(2)类似地，3阶系数可以展开为q3 = q30 + q31;rr + q31;tc tc (3)r = R0 � RrefRref qij;xy qi x yj其中， 是与场景中心相关的距离偏移， 是场景中心斜距， 是关于和的第阶泰勒展开式系数，高阶量一般忽略不计。SZVYXOTqbqa图 1 MEO SAR成像几何示意图Fig. 1 MEO SAR imaging geometry第5 期李 航等：基于成像坐标系优化的中轨星载SAR成像方法8573 最优成像坐标系聚焦算法3.1 最优成像坐标系图2给出了场景中不同目标多普勒调频率在不同成像坐标系中的分布情况，黑色曲线组表示等多普勒调频率线。图2(a)表示原始回波录取坐标系，横轴表示方位时间，纵轴表示距离，在该坐标系中由于卫星速度随时间变化以及轨道的弯曲，多普勒调频率存在沿方位向的非线性变化[21]。本文通过调整场景中的多普勒调频率分布来解决方位空变的问题。图2(b)为非正交坐标系中将数据进行斜置后的多普勒调频率分布，对不同方位时间的脉冲信号进行了线性距离平移操作，此时多普勒调频率方位的方位线性变化量得到了消除，只剩余了二次变化。在图2(c)所示的非线性非正交坐标系，经过对不同时刻的信号进行了非线性的平移，可以使等多普勒调频率线沿水平向分布，信号在方位向上满足平移不变性。从图2(a)到图2(b)和图2(c)的成像坐标系变化能够利用类似距离徙动校正的距离偏移函数来实现，而不同方位时刻的距离徙动量根据多普勒调频率的分布得到。为了最小化多普勒调频率方位向上的线性与非线性变化，距离移动量包含线性项与二次项ΔRNONL (ta) = a1ta + a2t2a (4)则最优成像坐标系变换函数可以表示为HNONL = exp(�j4π (fc + fr)cΔRNONL (ta))(5)fc c fra1 a2其中，定义为载频， 为电磁波的传播速度， 为距离向频率。参数与需要根据多普勒调频率的空变特性求解。在最优成像坐标系变换之后，斜距历程为R′(ta)  R0 +Σ4i=1qi(ta � tc)i + ΔRNONL (ta) (6)在最优成像坐标系变换之后，多普勒如fdc (ta) = � 2(q1 + a1 + 2a2tc + 2 (q2 + a2) (ta � tc))+ O (ta) (7)O (ta) ta fdc (t′c) = fdc (0)ta其中为的高阶项。令，得到新的方位聚焦时间， 的高阶项影响较小可以忽略，那么有t′c =q2tcq2 + a2 q20tcq20 + a2=
tc (8)其中，用
 = q20/ (q20 + a2)表示方位时域尺度变化因子。方位时移会引起目标聚焦位置的变化，可近似表示为ΔRshift = �(1 �
)[(q1 + a1) tc � a2t2c](9)那么，新的斜距历程变为R′(ta) = R′0 +Σ4i=1ki(ta � t′c)i (10)R′0 其中， 为新的聚焦距离，具体表达式为R′0 R0 +(a1 � a2q1q20)t′c +(a2 � a22q20)(t′c)2 (11)ki; i = 1; 2; ···; 4 t′是在新的方位聚焦时间 c下新的距离系数。具体表达式为k1=q1 + a1k2 = q20 + a2 + +q21;rrnew + q22;rrr2new+(1aq21;tc� q21;r(a1 � q1a2q20)� 3a2q30q20+(q22;tcr � 2(a1 � q1a2q20)q22;rr � 3q31;rq20)rnew)t′c+(1
2 q22;tctc� q21;ra2 � 3a2q31;tcq20+3q31;ra2q20(a1 � a2q1q20)� 1
q22;tc r(a1 � q1a2q20))(t′c)2k3 = q30+q31;rrnew +(1
q31;tc� q31;r(a1 � a2q1q20))t′ck4 = q49>>>>>>>>>>>>>>=>>>>>>>>>>>>>>;(12)新的方位空变的距离系数也在式(13)中给出k21;tc =1
q21;tc� q21;r(a1 � a2q1q20)� 3a2q30q20k22;tctc = q22;tctc +(2q22;tctck20� q21;r � 3q31;tcq20+q1q22;tc rq20)a2 +(3q31;r � q22;tc r) a1a2q20� q22;tc ra19>>=>>;(13)k2 k21;tc k22;tctca1 a2为了使距离单元偏移(Range Cell Migration, RCM)的方位空变最小，在参考距离处进行最优成像坐标系变换后，假设在方位向2阶系数为常数，令与为零，可以得到最优坐标系转换函数中的参数与858 雷 达 学 报第9 卷a2 =q20 (q22;tctcq21;r � q22;tc rq21;tc )[(2q22;tctc� q20q21;r � 3q31;tc + q1q22;tc r) q21;r + (3q31;r � q22;tc r) q21;tc ]a1 =q21;tcq21;r+q21;tc+q21;rq1 � 3q30q21;rq20a29>>=>>;(14)k22;rtc k31;tc 另外，由和造成的空变相位误差相对较小，可以利用3.2节提出的多普勒重采样将其处理。3.2 聚焦算法在上述最优成像坐标系转换的基础上，本文提出了一种新的MEO SAR聚焦算法，如图3所示。该算法首先利用最优成像坐标系转换函数处理多普勒调频率空变的问题；接着，利用距离非线性频调变标算法(Nonlinear Chirp Scaling Algorithm,N C S A ) [ 2 2 ]或距离徙动算法( R a n g e M i g r a t i o nAlgorithm, RMA)[23]补偿距离徙动；然后，用多普勒重采样补偿多普勒参数且校正最优成像坐标系转换造成的聚焦时移；最后利用几何形变校正来处理算法带来的聚焦位置偏移，从而获得最终的聚焦图像。在最优成像坐标系变换步骤中，没有考虑高阶距离系数的变化。三次距离系数的方位空变可能会导致非中心目标的旁瓣不对称，二次距离系数中的交叉项可能会导致场景角落的目标聚焦深度下降。为了实现更精确的方位聚焦，在最优成像坐标系变换之后，数据变换到方位时域之前，需要校正高阶多普勒参数和距离-方位耦合的空变分量。文献[21]提出使用与距离相关的多普勒重采样实现精确的方位聚焦。本文提出使用基于文献[21]的改进多普勒重采样方法实现多普勒线性化与方位时移校正。此改进的多普勒重采样方法仅多普勒重采样函数与文献[21]不同，其余多普勒重采样步骤请参考文献[21]。在多普勒域，方位匹配滤波器可以被表达为Hac (r; fa)=exp (�j2πφ) (15)其中，φ在式(16)中给出φ =(8k2;tc=0(fa +2a1)2+2k3032k32;tc=0(fa +2a1)3+3(9k230� 4k20k40)512k520)(fa +2a1)4k2;tc=0 = k20 + k21;rr + k22;rrr29>>>>>>>>>>>>=>>>>>>>>>>>>;(16)根据式(14)中的剩余空变和式(8)中目标的方位偏移因子，构造一个改进的多普勒重采样函数fan=
(fa + p21 (rnew) f2a + p31f3a)(17)p21 (rnew) p31 其中， 与是多普勒参数关于时间的导数，其表达式在式(18)中给出。p21 (rnew) =k22;rt′crnew8k22 (rnew; 0); p31 = �2k31;t′c32k320(18)根据式(18)在多普勒域中对数据进行重采样，能够消除残余的多普勒调制项，并校正方位聚焦偏移。成像坐标系变换引入的几何畸变，可以在方位压缩后通过在距离频域乘以线性校正函数式(19)来校正(a) 原始数据坐标系(a) Original coordinate system(b) 非正交坐标系(b) Non-orthogonal coordinate system(c) 非线性非正交坐标系(c) Non-orthogonal nonlinear coordinate systemRtcR′tcR″tc图 2 多普勒调频率平面坐标系Fig. 2 Doppler rate plane回波信号距离向FFT距离向NCS多普勒重采样方位IFFT几何校正聚焦成像最优成像坐标系变换方位压缩图 3 聚焦算法流程图Fig. 3 Flowchart of the proposed imaging algorithm第5 期李 航等：基于成像坐标系优化的中轨星载SAR成像方法859HGDC = exp(j4πfr (R0;new � R0)c)= exp(j4πfrc((
 �
k1k20)tc+(
 �
2k20)t2c))(19)完成几何校正后，即可获得良好的聚焦图像。4 仿真分析在本节中，将使用表1中的MEO SAR参数进行仿真实验，以验证所提出的方法。轨道高度为13000 km，在100 km的场景中布置了11 × 11个点目标(如图4所示)，能够实现2 m方位分辨率。为了评估本文所提方法对场景中不同点目标的聚焦效果，挑选3个典型的点目标A, B, C进行了分析，并且将本文方法与文献[21]中联合时间多普勒重采样算法(Joint Time-Doppler Resampling, JTDR)、常规的文献[22]中非线性频调变标算法(NCSA)进行对比分析。图5给出了在坐标系变换前和变换后，目标斜距与多普勒调频率的关系。通过对坐标系变换前后的斜距系数进行拟合，可以对比多普勒调频率的变化情况。在坐标系变换之前，相同斜距但具有不同波束中心时间的目标的多普勒调频率分布在较宽的范围内。在坐标系变换之后，所有目标在多普勒调频率-距离平面内严格分布在一条线上，说明此时方位向上多普勒调频率是恒定的。图6给出了文献[22]中NCSA算法得到的仿真结果，只有中心目标聚焦良好。图7给出了文献[21]中JTDR算法得到的仿真结果，选定的3个目标都聚焦良好，方位峰值旁瓣比(Peak SideLobe Ratio, PSLR)低于–13.26 dB，方位积分旁瓣比(Integrated Sidelobe Ratio, ISLR)低于–10.14 dB(具体数值见表2)，因此，文献[21]中的算法对于MEO SAR聚焦是有效的。图8给出了本文所提算法得到的仿真结果，选表 1 仿真参数Tab. 1 Simulation parameters类型名称值轨道参数轨道高度(km) 13000偏心率0倾角(°) 90近地点幅角(°) 0雷达参数载频(GHz) 5.2带宽(MHz) 105PRF(Hz) 830斜视角(°) 0入射角(°) 40合成孔径时间(s) 40.1地面距离/方位分辨率(m) 2/2场景参数方位场景幅宽(km) 100距离场景幅宽(km) 100-50 0 50A BC方位 (km)-50050 12345678910111 2 3 4 5 6 7 8 9 10 11距离 (km)图 4 仿真场景目标分布Fig. 4 Arrangement of simulated targets(a) 坐标系变换之前(a) Before coordinate system transformation-40 -20 0 20 40Range (km)-12.45-12.40-12.35-12.30-12.25-12.20-12.15Doppler rate (Hz/s)(b) 坐标系变换之后(b) After coordinate system transformationDoppler rate (Hz/s)-50 0 50Range (km)-10.80-10.75-10.70-10.65-10.60-10.55-10.50图 5 目标斜距与多普勒调频率的关系Fig. 5 Relationships between Doppler rates and ranges of the simulated targets860 雷 达 学 报第9 卷(a) 目标A仿真结果(a) Simulation result of target A(b) 目标B仿真结果(b) Simulation result of target B(c) 目标C仿真结果(c) Simulation result of target C(d) 目标A方位向剖面图(d) Azimuth profile of target A(e) 目标B方位向剖面图(e) Azimuth profile of target B(f) 目标C方位向剖面图(f) Azimuth profile of target CAmplitude (dB) Range (samples)Amplitude (dB)Amplitude (dB)5004003002001000Range (samples)500400300200100Range (samples)500400300200100-5-10-15-20-25-30-35-400-5-10-15-20-25-30-35-400-5-10-15-20-25-30-35-400 100 200 300 400 500 01000 2000 3000 40001000 2000 3000 4000 0 1000 2000 3000 40001000 2000 3000 4000Azimuth (samples) Azimuth (samples) Azimuth (samples)Azimuth (samples) Azimuth (samples) Azimuth (samples)100 200 300 400 500图 6 文献[22]中NCS算法的点目标仿真结果Fig. 6 Simulation results using the NCS algorithm in Ref. [22](a) 目标A仿真结果(a) Simulation result of target A(b) 目标B仿真结果(b) Simulation result of target B(c) 目标C仿真结果(c) Simulation result of target C(d) 目标A方位向剖面图(d) Azimuth profile of target A(e) 目标B方位向剖面图(e) Azimuth profile of target B(f) 目标C方位向剖面图(f) Azimuth profile of target CAmplitude (dB) Range (samples)Amplitude (dB)Amplitude (dB)5004003002001000Range (samples)500400300200100Range (samples)500400300200100-5-10-15-20-25-30-35-400-5-10-15-20-25-30-35-400-5-10-15-20-25-30-35-400 100 200 300 400 500Azimuth (samples) Azimuth (samples) Azimuth (samples)Azimuth (samples) Azimuth (samples) Azimuth (samples)100 200 300 400 5000 100 200 300 400 500100 200 300 400 5000 100 200 300 400 500100 200 300 400 500图 7 文献[21]中的JTDR算法的点目标仿真结果Fig. 7 Simulation results using the JTDR algorithm in Ref. [21]第5 期李 航等：基于成像坐标系优化的中轨星载SAR成像方法861定的3个目标都聚焦良好，方位PSLR低于–13.25 dB，方位ISLR低于–10.06 dB(具体数值见表2)。因此，本文所提算法对于MEO SAR聚焦是有效的，并且聚焦效果与文献[21]中的JTDR算法基本一致。M NK下面针对本文算法的计算复杂度进行仿真分析。假设信号距离向有个采样点，方位向有个采样点。用距离向的非线性频调变标(NCSA)实现距离单元徙动校正。所提方法包括4个距离向的FFT/IFFT，两个方位向的FFT/IFFT, 4个相位乘法，1个方位向插值。sinc插值核的的长度是。总的计算复杂度是4NM2log2M + 2MN2log2N + 4MN + KMN= (2log2M + log2N + K + 4)MN (20)图9给出了3种算法运算量随场景大小的变化关系，可以看出本文算法的运算量介于改进的NCS算法和JTDR算法之间。这主要是多普勒重采样操作引入了较大的运算量，事实上如果适当放宽对场景角落点的聚焦深度要求，不进行重采样操作，那么本文算法的运算量将接近改进的NCS算法。可见，本文所提算法实现良好聚焦效果的同时，运算量比参考算法低。5 结论本文提出一种基于最优成像坐标系变换的MEOSAR数据处理算法。相对于分别处理两维空变的常规处理算法，本文提出一种非正交非线性坐标系转换，在该坐标系中信号满足方位平移不变性，因此可以采用类似L E O S A R的处理算法进行MEO表 2 文献[21]中JTDR算法与本文所提算法仿真PSLR及ISLR数值Tab. 2 Compare of simulated values of PSLR and ISLR using the JTDR algorithm in Ref. [21] and the proposed method目标坐标文献[21]中JTDR算法本文所提算法PSLR ISLR PSLR ISLR方位距离方位距离方位距离方位距离A(6,6) –13.26 –13.28 –10.17 –10.06 –13.25 –13.28 –10.08 –10.02B(6,11) –13.26 –13.29 –10.14 –10.07 –13.26 –13.28 –10.06 –10.03C(1,11) –13.32 –13.28 –10.21 –10.06 –13.27 –13.28 –10.09 –10.02(a) 目标A仿真结果(a) Simulation result of target A(b) 目标B仿真结果(b) Simulation result of target B(c) 目标C仿真结果(c) Simulation result of target C(d) 目标A方位向剖面图(d) Azimuth profile of target A(e) 目标B方位向剖面图(e) Azimuth profile of target B(f) 目标C方位向剖面图(f) Azimuth profile of target CAmplitude (dB) Range (samples)Amplitude (dB)Amplitude (dB)5004003002001000Range (samples)500400300200100Range (samples)500400300200100-5-10-15-20-25-30-35-400-5-10-15-20-25-30-35-400-5-10-15-20-25-30-35-400 100 200 300 400 500Azimuth (samples) Azimuth (samples) Azimuth (samples)Azimuth (samples) Azimuth (samples) Azimuth (samples)100 200 300 400 5000 100 200 300 400 500100 200 300 400 5000 100 200 300 400 500100 200 300 400 500图 8 本文所提算法的点目标仿真结果Fig. 8 Simulation results using the proposed method862 雷 达 学 报第9 卷SAR的数据处理。与常规方法相比，由于简化方位空变的校正处理，因此能够显著降低算法的运算量。另外，针对波束扫描模式(如滑动聚束/TOPS模式等)，需要采用子孔径处理、或者两步法方位解模糊处理，本文提出的新的空变处理方法，经过改进后结合预处理，也能够应用于其他模式的数据处理中。参 考 文 献李春升, 于泽, 陈杰. 高分辨率星载SAR成像与图像质量提升方法综述[J]. 雷达学报, 2019, 8(6): 717–731. doi: 10.12000/JR19085.LI Chunsheng, YU Ze, and CHEN Jie. Overview oftechniques for improving high-resolution spaceborne SARimaging and image quality[J]. Journal of Radars, 2019, 8(6):717–731. doi: 10.12000/JR19085.[1]SUN Guangcai, XING Mengdao, WANG Yong, et al. A 2-Dspace-variant chirp scaling algorithm based on the RCMe q u a l i z a t i o n a n d s u b b a n d s y n t h e s i s t o p r o c e s sgeosynchronous SAR data[J]. IEEE Transactions onGeoscience and Remote Sensing, 2014, 52(8): 4868–4880.doi: 10.1109/TGRS.2013.2285721.[2]SUN Guangcai, JIANG Xiuwei, XING Mengdao, et al.Focus improvement of highly squinted data based onazimuth nonlinear scaling[J]. IEEE Transactions onGeoscience and Remote Sensing, 2011, 49(6): 2308–2322.doi: 10.1109/TGRS.2010.2102040.[3]HU Bin, JIANG Yicheng, ZHANG Shunsheng, et al.Generalized omega-k algorithm for geosynchronous SARimage formation[J]. IEEE Geoscience and Remote SensingLetters, 2015, 12(11): 2286–2290. doi: 10.1109/LGRS.2015.2470516.[4]李春升, 杨威, 王鹏波. 星载SAR成像处理算法综述[J]. 雷达学报, 2013, 2(1): 111–122. doi: 10.3724/SP.J.1300.2013.20071.LI Chunsheng, YANG Wei, and WANG Pengbo. A reviewof spaceborne SAR algorithm for image formation[J].Journal of Radars, 2013, 2(1): 111–122. doi: 10.3724/SP.J.1300.2013.20071.[5]WONG F W and YEO T S. New applications of nonlinearchirp scaling in SAR data processing[J]. IEEE Transactionson Geoscience and Remote Sensing, 2001, 39(5): 946–953.doi: 10.1109/36.921412.[6]RANEY R K, RUNGE H, BAMLER R, et al. PrecisionSAR processing using chirp scaling[J]. IEEE Transactionson Geoscience and Remote Sensing, 1994, 32(4): 786–799.doi: 10.1109/36.298008.[7]SUN Guangcai, XING Mengdao, LIU Yan, et al. ExtendedNCS based on method of series reversion for imaging ofhighly squinted SAR[J]. IEEE Geoscience and RemoteSensing Letters, 2011, 8(3): 446–450. doi: 10.1109/LGRS.2010.2084562.[8]SHIN H S and LIM J T. Omega-K algorithm for spacebornespotlight SAR imaging[J]. IEEE Geoscience and RemoteSensing Letters, 2012, 9(3): 343–347. doi: 10.1109/LGRS.2011.2168380.[9]SHIN H S and LIM L T. Omega-k algorithm for airbornespatial invariant bistatic spotlight SAR imaging[J]. IEEETransactions on Geoscience and Remote Sensing, 2009,47(1): 238–250. doi: 10.1109/TGRS.2008.2002954.[10]ZHANG Tianyi, DING Zegang, TIAN Weiming, et al. A 2-D nonlinear chirp scaling algorithm for high squint GEOSAR imaging based on optimal azimuth polynomialcompensation[J]. IEEE Journal of Selected Topics inApplied Earth Observations and Remote Sensing, 2017,10(12): 5724–5735. doi: 10.1109/JSTARS.2017.2765353.[11]HUANG Lijia, QIU Xiaolan, HU Donghui, et al. Mediumearth-orbit SAR focusing using range doppler algorithmwith integrated two-step azimuth perturbation[J]. IEEEGeoscience and Remote Sensing Letters, 2015, 12(3):626–630. doi: 10.1109/LGRS.2014.2353674.[12]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.[13]LI Dexin, WU Manqing, SUN Zaoyu, et al. Modeling andp r o c e s s i n g o f t w o - d i m e n s i o n a l s p a t i a l - v a r i a n tgeosynchronous SAR data[J]. IEEE Journal of SelectedTopics in Applied Earth Observations and Remote Sensing,2015, 8(8): 3999–4009. doi: 10.1109/JSTARS.2015.2418814.[14]汪丙南, 向茂生. 地球同步轨道圆迹SAR三维分辨特性分析[J].雷达学报, 2012, 1(3): 314–322. doi: 10.3724/SP.J.1300.2012.20044.WANG Bingnan and XIANG Maosheng. Three-dimensionalresolution analysis for geosynchronous circular SAR[J].Journal of Radars, 2012, 1(3): 314–322. doi: 10.3724/SP.J.1300.2012.20044.[15]20 40 60 80 100场景大小 (km)00.51.01.52.02.53.0运算量 (×1011)本文算法NCSAJTDR图 9 不同算法运算量对比Fig. 9 Computation comparison using different methods第5 期李 航等：基于成像坐标系优化的中轨星载SAR成像方法863CHEN Jianlai, SUN Guangcai, WANG Yong, et al. ATSVD-NCS algorithm in Range-Doppler domain forgeosynchronous synthetic aperture Radar[J]. IEEEGeoscience and Remote Sensing Letters, 2016, 13(11):1631–1635.[16]TANG Shiyang, LIN Chunhui, ZHOU Yu, et al. Processingof long integration time spaceborne SAR data with curvedorbit[J]. IEEE Transactions on Geoscience and RemoteSensing, 2018, 56(2): 888–904. doi: 10.1109/TGRS.2017.2756109.[17]LI Zhuo, LI Chunsheng, YU Ze, et al. Back projectionalgorithm for high resolution GEO-SAR image formation[C].2011 IEEE International Geoscience and Remote SensingSymposium, Vancouver, Canada, 2011: 336–339.[18]RAN Lei, LIU Zheng, LI Tao, et al. An adaptive fastfactorized back-projection algorithm with integrated targetdetection technique for high-resolution and high-squintspotlight SAR imagery[J]. IEEE Journal of Selected Topicsin Applied Earth Observations and Remote Sensing, 2018,11(1): 171–183. doi: 10.1109/JSTARS.2017.2771503.[19]ZHANG Lei, LI Haolin, QIAO Zhijun, et al. A fast BPalgorithm with wavenumber spectrum fusion for highresolutionspotlight SAR imaging[J]. IEEE Geoscience andRemote Sensing Letters, 2014, 11(9): 1460–1464.[20]LIU Wenkang, SUN Guangcai, XIA Xianggen, et al. Amodified CSA based on joint Time-Doppler resampling forMEO SAR stripmap mode[J]. IEEE Transactions onGeoscience and Remote Sensing, 2018, 56(6): 3573–3586.doi: 10.1109/TGRS.2018.2802545.[21]HUANG Lijia, QIU Xiaolan, HU Donghui, et al. Focusing ofmedium-earth-orbit SAR with advanced nonlinear chirpscaling algorithm[J]. IEEE Transactions on Geoscience andRemote Sensing, 2011, 49(1): 500–508. doi: 10.1109/TGRS.2010.2053211.[22]LIU Wenkang, SUN Guangcai, XIA Xianggen, et al. Highlysquinted MEO SAR focusing based on extended Omega-Ka l g o r i t h m a n d m o d i f i e d j o i n t t i m e a n d D o p p l e rresampling[J]. IEEE Transactions on Geoscience andRemote Sensing, 2019, 57(11): 9188–9200. doi: 10.1109/TGRS.2019.2925385.[23]作 者 简 介李　航(1996–)，女，山西吕梁人，西安电子科技大学电子工程学院博士生。主要研究方向为星载合成孔径雷达成像、海洋微波遥感观测等。E-mail: hli_xidian@163.com刘文康(1994–)，男，河南周口人，西安电子科技大学电子工程学院博士生。主要研究方向为中轨星载合成孔径雷达成像、高轨星载合成孔径雷达成像、动目标成像处理等。E-mail: wkliu@stu.xidian.edu.cn孙光才(1984–)，男，湖北孝感人，博士，华山特聘教授。2012年在西安电子科技大学电子工程学院获得博士学位，现担任西安电子科技大学电子工程学院副教授。主要研究方向为新体制雷达成像、运动目标检测成像。E-mail: rsandsgc@126.com邢孟道(1975–)，男，浙江嵊州人，博士，教授，2002年在西安电子科技大学电子工程学院获得博士学位，现担任西安电子科技大学电子工程学院教授。研究方向为雷达探测、雷达成像、运动目标检测成像。E-mail: xmd@xidian.edu.cn864 雷 达 学 报第9 卷