Geolocation Error Transfer Model and Trajectory Calibration Method for Airborne SAR Considering Motion Compensation Residual Error
-
摘要: 机载合成孔径雷达(SAR)定位误差不仅受载机位置/速度测量误差、系统时间误差等的影响,还与运动补偿残余误差有关。然而现有机载SAR定位模型很少考虑运动补偿误差的影响。该文针对实际中普遍存在的含运动误差和载机航迹测量误差的情况,结合运动补偿和频域成像算法,推导了机载SAR图像定位误差传递模型,阐明了运动补偿残余误差影响下航迹测量误差对定位偏差的影响方式,并基于该模型给出了载机航迹测量误差的标定方法。仿真实验验证了该定位误差传递模型的正确性,相比于不考虑运动补偿残余误差的定位模型,得到了更高精度的航迹测量误差标定结果,证明了该方法的优越性。Abstract: Airborne Synthetic Aperture Radar (SAR) location error is affected by the position/speed measurement error of the aircraft, system time error, etc., and also related to the residual error of motion compensation. However, the existing airborne SAR location model rarely considers the effect of residual motion error. Considering that motion and trajectory measurement errors are common in practice, this paper derives a location error transfer model of an airborne SAR image based on the motion compensation and frequency-domain imaging algorithms. The proposed model clarifies the influence of trajectory measurement error on location deviation when residual motion error exists and provides a method of error calibration measurement. The simulation experiments validate the correctness of the proposed location error transfer model. The present method obtains a more accurate error calibration measurement result than the location error model that does not consider the residual motion error, proving the superiority of the proposed model.
-
1. 引言
由于脉冲雷达发射的大功率信号在城市环境下容易干扰电台等其他通讯设备,可以借助通信信号来探测目标。近年来出现的多载波调制(Multi-Carrier Modulation, MCM)技术[1]引起了人们的关注,由此发展而来的正交频分复用(Orthogonal Frequency Division Multiplex, OFDM)[2]技术采用多路正交子载波进行信号调制,具有频率分集和波形分集的潜力,在通信、雷达一体化的发展背景下有着重要的研究价值。
相位编码OFDM信号具有多普勒高分辨力,用于雷达动目标检测时可以更精确地估计目标速度。不同的相位编码信号具有不同的自相关特性。霍夫曼编码可以降低信号自相关函数的整体旁瓣水平,但其包络峰均比(Peak-to-Mean Envelope Power Ratio, PMEPR)明显增大[3]。研究表明,barker码序列具有较好的综合性能[4],本文将以13位barker码作为相位编码序列。
最近10年,国内外在OFDM信号特性与波形设计上取得了较好的研究成果[5–7],因此OFDM雷达的信号处理问题受到了研究者的关注。张卫等人[8]利用Keystone变换在信号子载波域、快时间域和慢时间域进行联合解耦合处理,解决了目标的跨距离-多普勒单元走动问题,进而可估计出匀速运动下的多目标参数信息,但计算量较大;Lellouch等人[9]利用回波与发射信号的载频相位信息得到了点目标距离及径向速度估计,但要求目标在一个距离门内运动,即不发生越距离单元走动现象。除此以外,还需进一步研究OFDM雷达回波处理中面临的一些特殊问题,如速度补偿和多普勒解模糊。
在信号处理领域,最大似然估计是一种渐进有效估计量,但是对于多测量的非线性模型而言其计算量较大,不利于实际应用。为了提高计算效率,本文借鉴文献[10]中MIMO雷达信号处理的思路,结合OFDM信号多载波正交结构的特点,对信号进行通道分离,形成多通道信号。通过相关处理得到不同子载波上的距离像;利用Keystone变换进行速度补偿并解多普勒模糊,对同一载波的相同距离单元进行脉冲多普勒处理,得到每个子载波对应的多普勒频谱;进一步在子载波域作相参积累,得到距离-多普勒2维谱。通过谱峰搜索和CLEAN技术[11]的运用,从中提取出峰值位置对应的时延和多普勒参数。将其作为初值,结合观测数据的似然函数,利用牛顿迭代法获得更精确的参数估计,本文称之为近似最大似然估计。近似最大似然估计量可构成复合假设检验的重要环节,提升检测器的目标检测性能。论文组织结构如下:第2节给出了相位编码OFDM信号的回波模型;第3节给出了目标距离和速度参数的最大似然估计模型;为了提高目标运动参数估计的运算效率,第4节提出一种基于通道分离的近似最大似然估计算法;第5节利用仿真实验验证了算法的性能;第6节是总结。
2. 相位编码OFDM信号模型
设雷达发射相位编码OFDM信号为:
sT(t)=ej2πfctN−1∑n=0u(t−nTr)rect(t−nTrTr) (1) 其中,
rect(t)={1,0≤t≤10,else 表示窗函数,fc 为雷达信号的工作频率,N是脉冲个数,Tr 是脉冲重复周期,u(t) 为相位编码OFDM信号的复包络:u(t)=K−1∑k=0M−1∑m=0ak,mej2πkΔftrect(t−mtctc) (2) 其中,K是载波个数,第k个子载波上的相位编码序列
ak,m 由M个码元组成,tc 为码元宽度,B是信号总带宽,为了满足子载波间的正交性,子载波间隔Δf=B/K=1/tc 。对目标散射的回波进行下变频处理,可获得相应的OFDM基带信号:
sr(t)=A0e−j2πfctsT(t−τ) (3) 其中,
A0 为目标散射强度,τ=2(R−vt)c= τ0−2vtc 是匀速运动目标对应于t 时刻的时延。τ0=2Rc 表示t=0 时刻的目标距离时延,v 是目标的径向速度。若满足条件2vtc≪tc ,则式(3)可以简化为:sr(t)≈N−1∑n=0K−1∑k=0M−1∑m=0A0ak,mrect(t−nTr−τ0Tr)⋅rect(t−nTr−τ0−mtctc)e−j2π(fc+kΔf)τ0⋅ej2πkΔftej2π(fd+fdk)te−j2πkΔfnTr (4) 其中,子载波间的多普勒频差
fdk=2vkcΔf 。记全时间t=nTr+˜t ,其中快时间˜t∈[0,Tr) ,则式(4)表示为:\begin{align} {s_{\rm r}}(\tilde t,n) =& \sum\limits_{k = 0}^{K - 1} \sum\limits_{m = 0}^{M - 1} {A_0}{a_{k,m}} {\rm{rect}}\left(\frac{{\tilde t - {\tau _0}}}{{{T_{\rm r}}}}\right)\\ & \cdot{\rm{rect}}\left( \frac{{\tilde t - {\tau _0} - m{t_{\rm c}}}}{{{t_{\rm c}}}}\right){{\rm{e}}^{{\rm{ - j}}2{\rm{{{π}} }}({f_{\rm c}} + k\Delta f){\tau _0}}}\\ & \cdot {{\rm{e}}^{{\,\rm{j}}2{\rm{{{π}} }}k\Delta f\tilde t} \;{{\rm{e}}^{{\,\rm{j}}2{\rm{{{π}} }}({f_{\rm d}} + {f_{{\rm d}k}}) (\tilde t + n{T_{\rm r}})}} \end{align} (5) 其中,第2个指数项表示不同子载波的回波信号具有不同的频率偏移;由第3个指数项可以看出,子载波分别与慢时间和快时间相耦合。
在加性高斯白噪声背景下,观测信号可以表示为:
y(˜t,n)=sr(˜t,n)+w(˜t,n) (6) 其中,
E(w(˜t,n))=0 ,Var(w(˜t,n))=σ2 。3. 最大似然估计
针对式(6)在快时间域采样,采样时刻为
˜t=mtc+pTs ,其中,Ts=1KΔf ,p∈[0,K−1] ,m∈[0,M−1] 。由于fdk˜t≪1 ,式(5)进一步化简sr(p,m,n)=K−1∑k=0A0ak,mrect(mtc+pTs−τ0Tr)⋅rect(pTs−τ0tc)e−j2π(fc+kΔf)τ0⋅ej2π(fd+kΔf)(mtc+pTs)ej2π(fd+fdk)nTr (7) 令
A=A0e−j2πfcτ0 ,则采样信号为:Y=AS+W (8) 其中,
S=[s0,···,sN−1]T ,Y=[y0,···,yN−1]T yn=[y(˜t0,n),···,y(˜tMK−1,n)]T sn=K−1∑k=0ak,m ej2π(fd+fdk)nTr ⋅e−j2πkΔfτ0[ej2π(fd+kΔf)˜t0,···, ej2π(fd+kΔf)˜tMK−1]T 记待估参数向量
u=[A,v,R]T ,噪声的协方差矩阵R=σ2I,则观测数据Y对应的似然函数为:p(Y|u)=N−1∏n=01πMK|R|−1e−(yn−Asn)HR−1(yn−Asn) (9) 式(8)是关于A的条件线性模型,其最大似然估计为:
ˆA=N−1∑n=0sHnR−1yn/N−1∑n=0sHnR−1sn (10) 将其代入式(9),则距离和速度参数的最大似然估计为:
[ˆR,ˆv]=argmin[R,v]N−1∑n=0(yn−ˆAsn)HR−1(yn−ˆAsn)=argmax[R,v]|N−1∑n=0sHnR−1yn|2/|N−1∑n=0sHnR−1yn|2N−1∑n=0sHnR−1snN−1∑n=0sHnR−1sn=(σ2KN)−1argmax[R,v]|N−1∑n=0K−1∑p=0M−1∑m=0K−1∑k=0y(p,m,n)⋅a∗k,mej2πkΔfτ0e−j2π(fd+kΔf)(mtc+pTs)⋅e−j2π(fd+fdk)nTr|2 (11) 针对该模型,通过由粗到精的网格搜索可以得到参数的估计值。记速度搜索范围为
[−vmax,vmax] ,距离搜索范围为[0,Rmax] ,则搜索次数对应为U=ceil(2vmax/Δv)+1 和V=ceil(Rmax/ΔR)+1 。其中,ceil(x) 表示不小于x 的最小整数,Δv 和ΔR 表示搜索步长。算法计算复杂度为O(UVKNlog2(KN)) 。令
∑K−1k=0|ak,m|2=1 ,由文献[12]可知,待估参数的Cramer-Rao下限为:CRLB(v)=c2σ2N/c2σ2N(32(πfcTrA0)2(32(πfcTrA0)2⋅[N2(N−1)(2N−1)/6−N2(N−1)2/4]) (12) CRLB(R)=3c2σ2/(8(πΔfA0)2N(K2−1)) (13) 4. 基于通道分离的近似最大似然估计
直接利用式(11)求最大似然估计时计算量较大,在此考虑提高运算效率的参数估计方法。式(11)中相位项
e−j2πfdknTr 和ej2πkΔf(mtc+pTs) 的子载波k分别与慢时间和快时间相耦合。针对不同的子载波,多脉冲联合处理时目标将在距离单元和多普勒单元上走动。鉴于此,文献[8]利用Keystone变换在信号子载波域、快时间域和慢时间域进行联合解耦合处理,解决了目标的跨距离多普勒单元走动问题,计算量约为3K3Nlog2(KN) 。事实上,式(11)中的相位项ej2πkΔfpTs 体现了不同子载波的回波信号具有不同的频率偏移,因此可以将多载波正交结构的OFDM雷达信号进行分离,通过多通道接收的方式增大距离分辨单元,避免目标的跨距离单元走动和3维Keystone变换。当|fd+fdk|<Δf/2 时,对式(5)进行通道分离,即用参考信号ej2πkΔf˜t 与回波混频并经过低通滤波处理,得到各子载波通道上的信号x(˜t,n,k)=M−1∑m=0A0ak,mrect(˜t−τ0Tr)⋅rect(˜t−τ0−mtctc)e−j2π(fc+kΔf)τ0⋅ej2πfd˜tej2π(fd+fdk)nTr (14) 其中,快时间域的采样时刻
˜t=mtc+pTs 。用相位编码信号作为参考对式(14)作相关处理,得到第i个距离单元上的信号
zk(˜t,n,i)=M−1∑m,l=0A0ak,ma∗k,lrect(˜t−τ0Tr−|τ0−itc|)⋅rect(˜t−τ0˜t−|τ0−(i+l−m)tc|)⋅e−j2π(fc+kΔf)τ0ej2πfd˜tej2π(fd+fdk)nTr (15) 此时,信号的距离分辨单元
ΔR=c/(2Δf) 相对较大,目标不易产生跨距离单元走动。当最大多普勒频差小于多普勒分辨单元,即((2v)/c)(K−1)Δf< 1/(NTr) 时,子载波和慢时间之间的耦合可以忽略,本文称之为低速运动,否则,称之为高速运动。此时,可利用Keystone变换消除子载波k与慢时间n之间的耦合,即(fd+fdk)n=fdh ,实现多普勒维对齐。解耦合后的信号gk(˜t,h,i)=zk(˜t,fcfc+kΔfh,i)≈M−1∑m,l=0A0ak,ma∗k,lrect(˜t−τ0Tr−|τ0−itc|)⋅rect(˜t−τ0˜t−|τ0−(i+l−m)tc|)⋅e−j2π(fc+kΔf)τ0ej2πfd˜tej2πfdhTr (16) 经Keystone变换后的信号在慢时间域可能出现多普勒速度模糊现象(
fd=˜fd+r/Tr ,其中|˜fd|<1/(2Tr) ,折叠因子r∈{−8,−7,···,7,8} ),则式(16)改为:gk(˜t,h,i)≈M−1∑m,l=0A0ak,ma∗k,lrect(˜t−τ0Tr−|τ0−itc|)⋅rect(˜t−τ0˜t−|τ0−(i+l−m)tc|)⋅e−j2π(fc+kΔf)τ0ej2πfd˜tej2π˜fdhTr⋅ej2πrfcfc+kΔfh (17) 经过折叠因子补偿后的信号
fk(˜t,h,i)=gk(˜t,h,i)e−j2πrfcfc+kΔfh (18) 对式(18)进行脉冲多普勒处理,即关于慢时间作相参积累,可得第k个子载波上的多普勒频谱
Fk(˜t,u,i)=M−1∑m,l=0A0ak,ma∗k,lrect(˜t−τ0Tr−|τ0−itc|)⋅rect(˜t−τ0˜t−|τ0−(i+l−m)tc|)⋅sin(πs)sin(πs/N)ejπN−1Nse−j2π(fc+kΔf)τ0⋅ej2πfd˜t (19) 其中,
s=u−˜fdNTr ,u∈[0,N−1] 。进一步,关于k作相参积累可得距离-多普勒谱:
F(˜t,u,i)=M−1∑m,l=0A0ak,ma∗k,lrect(˜t−τ0Tr−|τ0−itc|)⋅rect(˜t−τ0˜t−|τ0−(i+l−m)tc|)⋅sin(πs)sin(πs/N)sin(πp)sin(πp/K)⋅ejπN−1NsejπK−1Kpe−j2πfcτ0ej2πfd˜t (20) 其中,
p=q+KΔfτ0 ,q∈[0,K−1] 。上述基于通道分离和Keystone相结合的方法是对式(11)模型的简化处理。为了进一步提高算法的估计精度,将得到的参数值作为初值,根据式(11)进行牛顿迭代,得到近似最大似然估计量。算法计算量约为
KNlog2(KN)+Nlog2N+2Klog2K+K2N 。上述模型均是以单目标为例进行讨论的。在多目标情况,可以利用CFAR检测器在距离-多普勒谱域进行谱峰搜索,根据最大峰值对应的参数估计值重构相应的子信号
b(˜t,n)=K−1∑k=0ej2πkΔfˆτ0e−j2π(ˆfd+kΔf)˜te−j2π(ˆfd+ˆfdk)nTr (21) 利用CLEAN技术从原信号中减去重构的子信号,然后对剩余信号继续上述的操作,直到CFAR检测不出峰值为止。算法处理流程如图1所示。
5. 仿真实验
仿真参数设置如表1所示。根据表中数据并结合前文的分析可知,目标产生速度模糊的阈值是
vuna=15m/s ,出现跨距离单元走动的速度临界值是vt≈4.68m/s ,速度分辨率Δvr=0.06m/s ,且满足(2vTr)/c≪tc 。由于|fd+fdk|≪Δf/2 ,多普勒频偏导致的子载频间串扰可忽略不计。在高斯白噪声背景下,2个目标距离为R1=10.0km ,R2=1.8 km ,R3=17.0km ,径向速度v1=20m/s ,v2=24m/s ,v3= 34 m/s,对应的目标散射强度分别为A01=5 ,A02=2 ,A03=1 。虚警率Pfa=10−3 ,蒙特卡洛仿真100次。表 1 仿真参数Table 1. Simulation parameters参数 数值 工作频率 fc (GHz) 5 带宽B (MHz) 64 脉冲重复周期 Tr (ms) 2 脉冲数N 250 载波数K 64 信号经过通道分离和相关处理后,在子载波-多普勒平面的投影如图2所示。由图2可以看出,目标的多普勒频移随子载波的变化而变化,因此不能直接进行子载波域的相参积累。
经过Keystone变换和CLEAN处理后,相应子载波-多普勒平面的投影如图3所示。可见,多普勒频移与子载波之间的耦合得到校正。
对Keystone变换的数据在子载波域进行相参积累,获得信号的距离-多普勒2维谱如图4所示。其中图4(a)是聚焦于第1个目标的结果,图4(b)是剔除前两个目标后聚焦于第3个目标的结果。结果表明,本文所述的补偿方法能在快时间域、慢时间域以及子载波域将目标的回波能量积累起来,有利于目标检测和后续的参数估计。
下面考察近似最大似然估计量的性能。当
v=20m/s ,SNRi 在0~20 dB之间变化时,分别利用文献[8]所述的联合Keystone变换法、通道分离与Keystone变换组合法以及本文所提的近似最大似然估计法进行速度估计,估计量的均方根误差(RMSE)曲线如图5(a)所示。图5(b)给出了各种算法对应的计算量随脉冲数的变化情况。由图5(a)可知,文献[8](点划线)与基于通道分离和Keystone相结合的方法(实线)的速度估计的RMSE曲线变化趋势类似。在相同的估计精度下,与上述两种方法相比,本文基于通道分离的近似最大似然估计法(点线)的输入信噪比
SNRi 改善约4 dB。同时,本文算法得到的RMSE非常接近Cramer-Rao下限。由图5(b)可知,在相同脉冲数下,本文算法的计算量较文献[8]的算法大幅降低。由于使用了牛顿迭代法,与基于通道分离和Keystone相结合的方法相比,计算量有所增加,但仍然远小于文献[8]中算法的计算量。因此,本文提出的基于通道分离的近似最大似然估计算法在估计精度和计算复杂度上具有综合优势。6. 总结
本文将OFDM通信信号应用在雷达动目标探测中,在通信、雷达一体化的发展背景下有着重要的应用前景。本文参考多载频MIMO雷达通道分离得到目标高分辨距离信息的方法,将OFDM信号的多载波正交结构与脉冲多普勒处理相结合,并借助Keystone变换解决了多普勒偏移问题。为了得到更好的估计精度,利用牛顿迭代法对似然函数进行优化,得到了基于通道分离的近似最大似然估计方法。仿真结果验证了算法的综合性能。今后还可针对机动目标相干化处理以及参数估计问题展开研究,扩展OFDM雷达的应用范围。
-
表 1 定位误差仿真参数
Table 1. Simulation parameters of location error
参数 数值 载波频率(GHz) 30 信号带宽(MHz) 200 脉冲持续时间(μs) 1 脉冲重复频率(Hz) 1000 方位向天线尺寸(m) 2 平台速度(m/s) 70 斜视角(°) 0 目标点1坐标(m) (330,0,0) 目标点2坐标(m) (400,0,0) 目标点3坐标(m) (470,0,0) 表 2
{\boldsymbol{y}} 方向定位误差Table 2. Location error in y direction
目标点 实际定位误差 公式算得误差 两误差的偏差 目标点1 (m) –4.5500 –4.5254 –0.0246 目标点2 (m) –1.6100 –1.6270 0.0170 目标点3 (m) 1.3300 1.3030 0.0270 表 3
{\boldsymbol{x}} 方向定位误差Table 3. Location error in x direction
目标点 实际定位误差 公式算得误差 两误差的偏差 目标点1 (m) 3.4980 3.4192 0.0788 目标点2 (m) 1.2034 1.0431 0.1603 目标点3 (m) –0.4437 –0.6014 0.1577 表 4 测量误差标定结果
Table 4. Calibration results of measurement error
方向 测量航迹 真实航迹 标定出的测量误差 标定航迹 x方向 {\eta ^2} - 3\eta - 1 {\eta ^2} - 3.0161\eta - 0.6755 {\eta ^2} + 0.0161\eta - 0.3245 z方向 {\eta ^2} - 3\eta + 449 {\eta ^2} + 450 - 3.0113\eta - 0.7585 {\eta ^2} + 0.0113\eta + 449.7585 表 5 测量误差标定结果
Table 5. Calibration results of measurement error
方向 测量航迹 真实航迹 标定出的测量误差 标定航迹 x方向 {\eta ^2} - 3\eta - 1 {\eta ^2} - 3.0161\eta - 0.{\rm{2751}} {\eta ^2} + 0.0161\eta - 0.{\rm{7249}} z方向 {\eta ^2} - 3\eta + 449 {\eta ^2} + 450 - 3.0{\rm{222}}\eta - 0.{\rm{3882}} {\eta ^2} + 0.0{\rm{222}}\eta + 449.{\rm{3882}} -
[1] 郑波浪. 机载高分辨率合成孔径雷达运动补偿研究[D]. [硕士论文], 中国科学院电子学研究所, 2006.ZHENG Bolang. A research on motion compensation of high-resolution airborne SAR[D]. [Master dissertation], Institute of Electronics, Chinese Academy of Sciences, 2006. [2] CANTALLOUBE H M J and NAHUM C E. Airborne SAR-efficient signal processing for very high resolution[J]. Proceedings of the IEEE, 2013, 101(3): 784–797. doi: 10.1109/JPROC.2012.2232891 [3] 苗慧. 机载SAR定位精度的研究[D]. [博士论文], 中国科学院电子学研究所, 2007.MIAO Hui. Research on airborne SAR geolocation accuracy[D]. [Ph. D dissertation], Institute of Electronics, Chinese Academy of Sciences, 2007. [4] 高祥武, 黄广民, 杨汝良. 机载SAR目标快速定位方法和定位精度分析[J]. 现代雷达, 2004, 26(9): 4–7. doi: 10.3969/j.issn.1004-7859.2004.09.002GAO Xiangwu, HUANG Guangmin, and YANG Ruliang. Study on a fast target location method for airborne SAR and location precision analysis[J]. Modern Radar, 2004, 26(9): 4–7. doi: 10.3969/j.issn.1004-7859.2004.09.002 [5] 彭代强, 林幼权, 杜鹏飞. 机载SAR图像快速经纬度计算及精度分析[J]. 现代雷达, 2010, 32(3): 48–52, 79. doi: 10.3969/j.issn.1004-7859.2010.03.013PENG Daiqiang, LIN Youquan, and DU Pengfei. A fast algorithm for latitude and longitude calculation of airborne SAR and its location precision analysis[J]. Modern Radar, 2010, 32(3): 48–52, 79. doi: 10.3969/j.issn.1004-7859.2010.03.013 [6] LIU Junbin, QIU Xiaolan, HUANG Lijia, et al. Curved-path SAR geolocation error analysis based on BP algorithm[J]. IEEE Access, 2019, 7: 20337–20345. doi: 10.1109/ACCESS.2019.2897361 [7] XING Mengdao, JIANG Xiuwei, WU Renbiao, et al. Motion compensation for UAV SAR based on raw radar data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(8): 2870–2883. doi: 10.1109/TGRS.2009.2015657 [8] DE MACEDO K A C and SCHEIBER R. Precise topography- and aperture-dependent motion compensation for airborne SAR[J]. IEEE Geoscience and Remote Sensing Letters, 2005, 2(2): 172–176. doi: 10.1109/LGRS.2004.842465 [9] MAO Yongfei, XIANG Maosheng, WEI Lideng, et al. Error analysis of SAR motion compensation[C]. 2012 IEEE International Conference on Imaging Systems and Techniques Proceedings, Manchester, UK, 2012: 377–380. doi: 10.1109/IST.2012.6295562. [10] 李芳芳, 仇晓兰, 孟大地, 等. 机载双天线InSAR运动补偿误差的影响分析[J]. 电子与信息学报, 2013, 35(3): 559–567. doi: 10.3724/SP.J.1146.2012.00850LI Fangfang, QIU Xiaolan, MENG Dadi, et al. Effects of motion compensation errors on performance of airborne dual-antenna InSAR[J]. Journal of Electronics &Information Technology, 2013, 35(3): 559–567. doi: 10.3724/SP.J.1146.2012.00850 [11] PRATS P, REIGBER A, MALLORQUI J J, et al. Efficient detection and correction of residual motion errors in airborne SAR interferometry[C]. IGARSS 2004. 2004 IEEE International Geoscience and Remote Sensing Symposium, Anchorage, USA, 2004: 992–995. doi: 10.1109/IGARSS.2004.1368576. [12] REIGBER A, PRATS P, and MALLORQUI J J. Refined estimation of time-varying baseline errors in airborne SAR interferometry[J]. IEEE Geoscience and Remote Sensing Letters, 2006, 3(1): 145–149. doi: 10.1109/LGRS.2005.860482 [13] 刘云龙, 李焱磊, 周良将, 等. 一种机载SAR快速几何精校正算法[J]. 雷达学报, 2016, 5(4): 419–424. doi: 10.12000/JR16064LIU Yunlong, LI Yanlei, ZHOU Liangjiang, et al. A fast precise geometric calibration method for airborne SAR[J]. Journal of Radars, 2016, 5(4): 419–424. doi: 10.12000/JR16064 [14] REIGBER A, ALIVIZATOS E, POTSIS A, et al. Extended wavenumber-domain synthetic aperture radar focusing with integrated motion compensation[J]. IEE Proceedings-Radar, Sonar and Navigation, 2006, 153(3): 301–310. doi: 10.1049/ip-rsn:20045087 [15] MOREIRA A, MITTERMAYER J, and SCHEIBER R. Extended chirp scaling algorithm for air- and spaceborne SAR data processing in stripmap and ScanSAR imaging modes[J]. IEEE Transactions on Geoscience and Remote Sensing, 1996, 34(5): 1123–1136. doi: 10.1109/36.536528 [16] 孟大地. 机载合成孔径雷达运动补偿算法研究[D]. [博士论文], 中国科学院电子学研究所, 2006.MENG Dadi. Research on motion compensation algorithm for airborne SAR[D]. [Ph. D dissertation], Institute of Electronics, Chinese Academy of Sciences, 2006. [17] YANG Mingdong, ZHU Daiyin, and SONG Wei. Comparison of two-step and one-step motion compensation algorithms for airborne synthetic aperture radar[J]. Electronics Letters, 2015, 51(14): 1108–1110. doi: 10.1049/el.2015.1350 [18] FORNARO G, FRANCESCHETTI G, and PERNA S. On center-beam approximation in SAR motion compensation[J]. IEEE Geoscience and Remote Sensing Letters, 2006, 3(2): 276–280. doi: 10.1109/LGRS.2005.863391 [19] 曾乐天, 邢孟道, 陈士超. 基于窄波束和平地假设的运动补偿方向研究[J]. 电子与信息学报, 2014, 36(10): 2464–2468. doi: 10.3724/SP.J.1146.2013.01671ZENG Letian, XING Mengdao, and CHEN Shichao. The research on the direction of motion compensation according to the narrow beam and flat earth hypothesis[J]. Journal of Electronics &Information Technology, 2014, 36(10): 2464–2468. doi: 10.3724/SP.J.1146.2013.01671 [20] 杨鸣冬, 俞翔, 朱岱寅. 基于距离子带的机载SAR高精度多级空变运动补偿[J]. 航空学报, 2018, 39(2): 321557. doi: 10.7527/S1000-6893.2017.21557YANG Mingdong, YU Xiang, and ZHU Daiyin. High-precision space-variant motion compensation with multi-level processing for airborne SAR based on subswath[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(2): 321557. doi: 10.7527/S1000-6893.2017.21557 [21] 柳俊斌. 弱导航信息下的SAR自主定位与航迹修正技术研究[D]. [硕士论文], 中国科学院空天信息创新研究院, 2020.LIU Junbin. Research on SAR autonomous positioning and track correction technology under weak navigation information[D]. [Master dissertation], Aerospace Information Research Institute, Chinese Academy of Sciences, 2020. [22] 孙立军, 刘鑫, 刘兴春. 关于差分GPS精度的分析[C]. 第十届东北三省测绘学术与信息交流会论文集, 漠河, 2009.SUN Lijun, LIU Xin, and LIU Xingchun. Analysis on the accuracy of differential GPS-DGPS[C]. The 10th Academic and Information Exchange Meeting of Surveying and Mapping in Northeast China, Mohe, China, 2009. [23] BAMLER R and EINEDER M. Accuracy of differential shift estimation by correlation and split-bandwidth interferometry for wideband and delta-k SAR systems[J]. IEEE Geoscience and Remote Sensing Letters, 2005, 2(2): 151–155. doi: 10.1109/LGRS.2004.843203 -