Comparative Experiments on Separation Performance of Overlapping Scatterers with Several Tomography Imaging Methods
-
摘要: 层析技术因具有解译城区SAR影像上复杂叠掩场景的能力而备受关注。层析成像包含两个部分:估计散射体在高程向的分布和确定散射体在混叠像元内的真实数目。该文以中科院空天院峨眉数据的机载阵列干涉系统参数为基础,选取了若干代表性的方法,包括OMP, SLIM和MUSIC等层析谱估计方法以及BIC和GLRT等模型定阶方法,进行了模拟叠掩目标的层析反演实验,使用了克拉默-拉奥界和重建成功率来评估实验结果。实验表明:在机载阵列数很有限的条件下,(1)使用2阶统计量反演的高程估计量的方差比单个观测矢量反演结果的方差更小;(2)叠掩散射体间的幅度比、相位差和散射间距会影响层析算法解叠掩的成功率;(3)叠掩散射体间的相位差会使层析算法的高程估计发生偏差。Abstract: The tomographic technique has attracted much attention because of its ability to separate overlapping scatterers in urban Synthetic Aperture Radar (SAR) images. The general method of SAR Tomography (TomSAR) imaging combines the following two aspects: estimating the distribution of the scatterers in the elevation direction and determining the number of strong scatterers in an overlapped pixel. This study applied several sophisticated spectrum estimations (e.g., Orthogonal Matching Pursuit, Sparse Learning via Iterative Minimization and Multiple Signal Classification) and model order selection approaches (e.g., Bayesian information criterion and generalized likelihood ratio test) with highly technical potential to recover the simulated overlapping scatterers. This simulation experiment is based on the parameters of the AIRCAS X-band TomoSAR data from Emei, Sichuan, China. The Cramér-Rao Lower Bound (CRLB) and recovery probability are used to evaluate the performances of different methods for the separation of overlapped scatterers. The experimental results revealed the following: (1) the standard deviation of estimation using second-order statistics is smaller than that of a single observation vector, especially when the number of acquisitions is very limited; (2) the amplitude ratio, phase difference, and elevation spacing between overlapping scatterers will have a significant impact on the different kinds of algorithms; and (3) the phase difference between overlapping scatterers will make the phase center estimation of greedy algorithm or spectrum estimation algorithm biased.
-
1. 引言
合成孔径雷达(Synthetic Aperture Radar, SAR)可以实现全天时、全天候的高分辨率成像[1],星载SAR系统为全球地理遥感、地形测绘、形变监测等应用提供了丰富的高分辨数据。分布式星载SAR系统是一种特殊模式的合成孔径雷达系统,雷达发射机和接收机位于不同平台,这种收发分置的特点带来了很多优势,比如灵活的基线配置,多样的接收模式,可以提供更精细的地物分类信息等,具有单基情况无可比拟的优势。但是收发分置的特点也带来了很多新的挑战,尤其是时间同步、相位同步和波束同步等三大同步问题[2,3]。由于收发雷达不共频率源,导致了时间同步和相位同步问题[4],时间同步问题是指,发射雷达和接收雷达使用不同的晶振产生脉冲信号,导致了触发脉冲不一致,进而导致信号发射时刻与接收窗起始时刻不一致,使距离向的录取回波减少,时间同步误差会在回波中引入相位误差,并且晶振偏差会随时间累积,对回波数据产生较大影响,严重降低雷达图像质量[5,6]。为达到时间同步的要求,Tandem-X SAR系统采用频率一致性指标为10–8量级的晶振,并采用跳跃PRI的方法,确保600 s内,时间同步误差小于5 μs[7]。相位同步问题是指,由于晶振频率的差异,导致了发射信号调制频率和接收信号解调频率不一致,造成解调回波中存在沿方位时间变化的相位误差,该误差在成像后传递到图像域,导致方位向位置偏移、主瓣展宽、积分旁瓣比恶化等问题,并且相位误差最终会传递到干涉相位,造成干涉测高精度损失。TanDEM-X SAR系统采用脉冲交替对传的同步方案解决相位同步的问题。波束同步问题是指,由于收发天线的控制精度不同,使得收发天线无法在同一时刻完全覆盖地面同一区域,导致双基合成天线方向图的增益降低,同时接收天线无法获得全部的目标散射信号,造成有效观测幅宽损失。为保证系统系能,TanDEM-X SAR系统的波束指向精度控制在0.01°,最差情况下导致的天线增益损失小于0.1 dB。
为了消除相位同步带来的影响,可以建立专门的同步链路,文献[8]提出了卫星间连续波对传的同步方法,文献[9]和文献[10]进一步提出了脉冲交替对传的方法,文献[11]提出一种基于信号交换的调频连续波合成孔径雷达相位同步方法。无论采用哪种同步链路方式,目的都是获得相位误差信息,最后在进行数据处理时提取并补偿相位误差,这也是解决分布式卫星SAR相位同步的优选方法[12]。文献[13]建立了一种时间同步误差模型,分析了时间同步误差对成像和干涉测高性能的影响,但是没有分析时间同步和相位同步的联系。时间同步和相位同步来自同一误差源,都是由于收发平台采用了独立的晶振系统来校准系统时钟[14],所以两者之间存在一定的关系。波束同步的研究主要集中在波束同步控制方法,以及波束同步误差对雷达成像性能的影响,文献[15]分析了几种可行的飞行控制策略,同时分析了波束同步对收发平台和姿态控制精度的要求,文献[16]分析了系统波束控制精度需要满足的要求。但是,在不同的卫星姿态控制精度和波束指向精度情况下,波束同步误差对分布式星载SAR系统的综合影响还需要进一步分析。
本文结构安排如下。第2节介绍了时间同步误差模型和波束同步误差模型。第3节基于误差模型分析了时间同步和相位同步导致的相位误差之间的关系,并分析晶振稳定度需要满足的条件。第4节分析了波束同步误差对合成天线方向图、双基图像信噪比、多普勒去相干以及波束覆盖偏差的具体影响。第5节,进行仿真验证,为分布式SAR系统设计提供参考。
2. 误差模型
2.1 时间同步误差模型
收发雷达之间的时间同步误差模型可以表示为常数项、线性项和随机项3部分[13],随机项由频率源的短期稳定性决定,由于随机的特点,该项影响结果不具备积累性,一般只有在频率源稳定度较低或者所需相干积累时间较长时有较大影响[17–20],因此该模型可以表达为:
ε(t)=Δt+βt (1) 其中,
Δt 为常数项,表示初始状态主辅星间雷达累积的定时信号之差;βt 为线性项,由主辅星晶振频率偏差导致,β 可以表示为:β=ΔPRI1−ΔPRI2PRI=Δf1f1−Δf2f2≈Δf1−Δf2f0 (2) 其中,f0表示标称频率,
Δf1 和Δf2 分别表示发射雷达和接收雷达的频率偏移。当存在同步误差时,在以发射信号为参考的时间坐标系中,不同脉冲的回波在接收雷达采样窗中的位置会发生偏移[21]。2.2 波束同步误差模型
卫星姿态误差会影响卫星和地面目标的波束传播距离,同时影响回波多普勒频率,如图1所示,X轴从星下点指向卫星方向,Y轴指向方位向方向,Z轴指向距离向方向,卫星姿态误差可以分解为偏航角、俯仰角和横滚角,其中Rs表示斜距,R0表示最近斜距,P0表示姿态导引后期望的波束中心地面落点,P表示波束中心实际落点,
θy 表示偏航角,θp 表示俯仰角,θr 表示横滚角。可见,偏航角和俯仰角分别导致波束中心线绕X轴和Z轴旋转,进而偏离零多普勒面,均会导致斜视角。横滚角导致波束中心线绕Y轴旋转,所以横滚角引起波束在距离向移动,在距离向产生波束指向误差。根据几何关系可以得到,偏航角和俯仰角引起的斜视角分别为:
θsy=arctan(tanθy⋅sinθl) (3) θsp=arctan(tanθp⋅cosθl) (4) 其中,
θl 为下视角。横滚角引起的距离向波束指向误差为:θsr=θr (5) 3. 时间同步误差分析
主星发射的线性调频脉冲信号可以表示为:
s(τ,t)=rect(τTp)exp(jπkτ2)exp(jπf0τ) (6) 其中,
rect[⋅] 为矩形窗函数,τ 为距离向快时间,t为方位向慢时间,Tp为脉冲宽度,f0为发射信号载频,k为调频率。将时间同步误差模型代入发射信号,经辅星雷达接收解调后含有时间同步误差的回波信号为:sD(τ,t)=rect(τ−Rt/c−ε(t)Tp)⋅rect(t−tc)⋅exp(jπk(τ−Rt/c−ε(t))2)⋅exp(−j2πRt/λ)exp(−j2πf0ε(t)) (7) 其中,c为光速,
λ 为波长,tc为波束中心偏离时间,Rt为雷达和卫星之间的波束传播距离,R(t)=√R12+V21(t−t01)2+√R22+V22(t−t02)2 (8) 其中,R1和R2分别表示收发最近斜距,t01和t02分别表示波束中心穿越时刻,V1和V2分别表示收发卫星的速度。为简化分析,下文中,假设两星同时穿越波束中心时刻且以该时刻为方位时间零点,即t01=t02=0。
以主星发射信号,辅星接收信号为例,采用距离多普勒聚焦算法,推导时间同步误差对SAR成像的影响,假设主星雷达晶振为标准频率源,辅星雷达晶振含有误差项,即时间同步误差全部表现在辅星获得的图像上。
3.1 固定时间误差分析
辅星雷达只存在固定时间误差的情况下,即
ε(t)=Δt (9) 成像处理之后的辅星图像将存在时间同步误差项,具体可以表示为:
σ(τ,t)=sinc(B(τ−(R1+R2)/c−Δt))⋅sinc(Bd⋅t)⋅exp(−j2πR1+R2λ)⋅exp(j2πfact)⋅exp(−j2πf0Δt) (10) 其中,Bd是多普勒带宽,fac是多普勒中心频率。由式(10)可以得到,固定时间同步误差在辅星图像中引入距离向主瓣偏移和固定相位误差,固定相位误差在经过干涉处理之后会产生常量干涉相位误差。主瓣偏移量和干涉相位误差分别为:
Δr=Δt⋅c2 (11) φe=−2πf0Δt (12) 常量相位误差会影响绝对测高精度。
3.2 时间同步与相位同步关系
辅星雷达只存在线性时间误差的情况下,即
ε(t)=βt (13) 对式(7)中带有线性时间误差的相位项进一步推导可以得到
φlinear=−2πβf0t=−2πΔf1−Δf2f0f0t=−2π(Δf1−Δf2)t (14) 由于晶振频率存在偏移,t时刻主星的晶振频率为
f1=f0+Δf1 ,辅星的晶振频率为f2=f0+Δf2 ,根据文献[9]中对脉冲对传相位同步方案的分析,辅星需要补偿的相位偏差为:φc(t)=2π(Δf2−Δf1)t+π(Δf2−Δf1)(τ12+τsys)−πfdτsys+12(nφ2(t+τsys)+nφ2(t+τ)−nφ1(t)−nφ1(t+τ+τsys))+12(φSNR1(t)−φSNR2(t))+Δφsys(t)+Δφant(t) (15) 其中,fd是多普勒频率,
τsys 是同步脉冲的发射延迟,nφ1 和nφ2 分别代表晶振的相位噪声,φSNR1 和φSNR2 分别代表主星和辅星的接收机噪声,Δφsys(t) 代表通道链路的相位噪声,Δφant(t) 代表天线方向图引入的相位噪声。式(15)的第1项
2π(Δf2−Δf1) ,等于式(7)中的线性偏差φlinear ,即随方位时间线性变化的相位误差部分可以通过相位同步方案补偿,补偿后的回波信号经过距离徙动校正后在方位频域的表达形式为:SR−com(τ,fa)=sinc(B(τ−R1+R2c+βfaKa))⋅Wa(fa−fac)⋅exp(−j2π(R1+R2)/λ)⋅exp(jπf2aKa) (16) 其中,Ka为方位向调频率。式(16)中
β(fa/faKaKa) 项会造成距离徙动校正误差,使方位信号弥散在若干距离单元上,导致主瓣展宽。如果要求距离徙动校正误差造成的主瓣展宽小于2%,那么要求未被校正的徙动量小于半个分辨率单元[22],即max(βKa⋅fa)<0.5⋅0.886B (17) 根据时频关系可以得到
βKa⋅Ka⋅Tsar<0.5⋅0.886B⇒β<0.443B⋅Tsar (18) 其中Tsar为合成孔径时间。
满足式(18)的情况下,经方位压缩后的SAR信号为:
σ(τ,t)=sinc(B(τ−R1+R2c))⋅sinc(Bd⋅t)⋅exp(−j2πR1+R2λ) ⋅exp(j2πfact) (19) 此时主要导致主瓣展宽,严重情况可能导致散焦。
4. 波束同步误差分析
理想情况下,SAR系统处于匀速运动状态,但是在实际工作情况下,卫星的运动受到控制系统精度的影响,会产生姿态误差,所以必须保证分布式系统的姿态控制精度[23]。波束同步误差可以分为两个方面,一方面是卫星姿态引起的误差,另一方面是波束指向误差。波束同步误差会导致天线照射范围不一致,导致信噪比损失,多普勒带宽偏移,方位向天线增益降低等问题[24]。同时,由于照射区域的偏差,还会导致距离向重叠幅宽减小。
4.1 对方位向影响
本节首先引入文献[7]中介绍的基于椭圆轨道的双星一发双收工作模式下的姿态导引方式,在此基础上,分析波束同步误差对多普勒去相干、天线增益以及双基图像信噪比损失的影响。
双星有两种姿态导引方式[7],一种是只针对主星进行零多普勒姿态导引,此时,场景内的回波多普勒中心频率不为0。另一种是主辅双星均进行零多普勒姿态导引,此时辅星与主星照射的共同区域等效为正侧视,即场景内多普勒中心频率为0。图2和图3给出了在主辅星顺轨基线从0到3 km情况下,不同姿态导引方式下方位向合成天线方向图的情况,其它仿真参数如表1所示。
图2中,由于仅主星进行零多普勒姿态导引,导致合成后的方位天线方向图指向等效于斜视,即场景回波的多普勒中心频率不为0,产生方位多普勒去相干,导致相干系数发生变化,如图4所示。
考虑到顺轨基线和波束同步误差的综合影响,相干系数的计算方法可以表示为:
γdc=1−|ΔfDop1+ΔfDop2Bd| (20) ΔfDop1 表示由顺轨基线导致的多普勒频率偏移,计算公式为[7]:ΔfDop1≈v⋅dalongλ⋅r (21) 其中,dalong为基线长度,r为平均斜距。
ΔfDop2 表示由波束指向误差和卫星姿态误差导致的多普勒频率偏移,ΔfDop2=2v⋅sin(θsq+Δθa)λ−2v⋅sin(θsq)λ (22) Δθa=θerr+θsy+θsp (23) 其中,
θsq 表示斜视角,Δθa 表示方位向波束同步误差,包括方位向波束指向误差以及偏航角和俯仰角对斜视角的影响部分。图3中,由于对主辅双星均进行零多普勒姿态导引,合成后的方位天线方向图等效于正侧视,此时对多普勒去相干影响很小。考虑到主辅双星之间顺轨基线的影响,合成后的天线方向图幅度会有一定的下降,从仿真中可以看出,在顺轨基线小于3 km的情况下,导致的幅度衰减小于0.05 dB,可以忽略。在不同顺轨基线情况下,主辅星同时进行零多普勒导引的合成天线方向图对频谱加权的计算公式为[7]:
Wbistat(fDop;ΔfDop)=sinc(l2v(fDop−ΔfDop))⋅sinc(l2v(fDop+ΔfDop)) (24) 其中,l为天线长度,v为卫星速度,
ΔfDop 是由顺轨基线导致的多普勒频率偏移。在此基础上,加入波束指向误差和卫星姿态误差的影响,顺轨基线和波束同步误差的综合加权公式可以表示为:
Wbistat(fDop;ΔfDop)=sinc(l2v(fDop−ΔfDop1−ΔfDop2))⋅sinc(l2v(fDop+ΔfDop1+ΔfDop2)) (25) 由于波束同步误差,导致成像后辅星图像信噪比降低,信噪比损失量为[25]:
SNRloss=∫θbw1/2−θbw1/2(sinc2(l⋅sinφ/λ)⋅sinc2(l⋅sin(φ−Δθa)/λ))dφ∫θbw1/2−θbw1/2(sinc2(l⋅sinφ/λ)⋅sinc2(l⋅sinφ/λ))dφ (26) 其中,
φ 为方位向离轴角,θbw1 为方位向波束宽度。4.2 对距离向影响
针对距离向,波束控制误差主要影响主、辅星的重叠幅宽,当主、辅星距离向同时存在相反方向的指向偏差时,距离向的重叠区域(图5和图6中黄色区域)会大幅减小,影响后续的干涉处理。根据几何模型,可以推导出覆盖区域偏差的计算公式:
ΔW=(H+Re)cos(θl+θbw2/2−Δθr)−√(H+Re)2cos2(θl+θbw2/2−Δθr)−2⋅H⋅Re−(H+Re)cos(θl+θbw2/2+Δθr)+√(H+Re)2cos2(θl+θbw2/2+Δθr)−2⋅H⋅Re (27) 其中,H为卫星高度,Re为地球半径,
θl 为下视角,θbw2 为俯仰向波束宽度,Δθr 为距离向波束同步误差。5. 仿真结果
仿真参数如表1所示。
表 1 仿真参数表Table 1. Simulation parameters参数 数值 载频(GHz) 1.30 卫星高度(km) 700 方位向天线长度(m) 6 脉冲宽度(μs) 80 距离向带宽(MHz) 100 PRF(Hz) 1600 卫星有效速度(m/s) 7500 距离向波束指向精度(°) 0.06 方位向波束指向精度(°) 0.05 姿态控制精度(°) 0.02 顺轨基线(km) 3 5.1 时间同步误差
仿真中加入参数为
Δt=100 ns 的固定时间误差,点目标分析结果如图7所示,通过计算可以得到理论上的距离向位置偏移为Δr=Δt⋅c/2=15 m ,仿真得到的对应偏差为15 m,如图7(b)所示。时间同步误差会引起主瓣展宽,如式(18)所示,仿真参数下,当
β<1.2×10−9 时,引起的展宽可以忽略。当β=5×10−8 时,距离向和方位向剖面图如图8所示,由仿真结果可以看到明显的主瓣展宽。5.2 波束同步误差
本节的仿真给出了当顺轨基线为3 km时,波束指向误差和卫星姿态误差对系统性能的影响。图9和图10分别比较了方位向波束指向精度在0.01°到0.05°情况下,对合成天线方向图及相干系数的影响,仿真结果表明,与无误差情况对比,0.05°波束指向误差造成方位天线增益下降小于0.05 dB,造成相干系数减小0.02。
图11和图12分别给出了卫星姿态控制精度在0.005°到0.020°的情况下,对方位向合成天线方向图及相干系数的影响,仿真结果表明,与无误差情况对比,0.020°的姿态控制误差造成方位天线增益下降小于0.02 dB,造成相干系数减小0.01。
不同波束指向和卫星姿态误差情况下,对合成天线方向图及相干系数的综合影响结果如图13和图14所示。波束指向精度为0.05°,卫星姿态控制精度0.02°的情况下,和无误差情况对比,造成方位天线增益下降小于0.1 dB,造成相干系数减小0.03。
不同波束指向和卫星姿态误差情况下,辅星图像信噪比损失情况如表2所示,0.05°的波束指向误差和0.02°的姿态误差导致的辅星图像信噪比损失小于0.1 dB。
表 2 波束同步误差造成的辅星图像信噪比损失(dB)Table 2. Influence of beam synchronization error on SNR (dB)波束指向精度(°) 姿态控制精度(°) 0.01 0.02 0.01 0.021 0.025 0.02 0.028 0.032 0.03 0.036 0.039 0.04 0.046 0.048 0.05 0.055 0.057 0°~0.06°的波束指向误差对距离向重叠幅宽的影响如图15所示。横滚角会在距离向产生波束控制误差,对距离向覆盖偏差产生影响,0.02°的姿态控制误差对重叠幅宽的影响如图16所示。横滚角从0°步进至0.02°,距离向波束指向精度从0°步进至0.06°情况下,两者对距离向重叠幅宽的综合影响如图17所示。上述结果表明,0.02°的姿态控制误差对距离向重叠幅宽的影响为368 m, 0.06°的波束指向误差对距离向重叠幅宽的影响为1100 m,由于波束指向和姿态控制导致的误差均为随机误差,所以综合影响采用几何平均数的形式,两者综合之后对距离向重叠幅宽的影响为1164 m。考虑到轨道控制精度的影响,在重轨工作模式下,两次测量中卫星的位置会在一定范围内偏移,在卫星姿态误差和波束指向误差基础上,系统设计中应综合考虑卫星位置偏移的影响。
由上述分析可知,在仿真精度情况下,波束同步精度对多普勒去相干、天线增益以及双基图像信噪比损失的影响较小,但对距离向重叠幅宽影响较大,因此需要在波位设计中保留一定的幅宽余量,或者根据幅宽需求,对波束同步精度提出要求。不同卫星姿态控制精度和波束指向精度对重叠幅宽的综合影响总结在表3可供系统设计参考。针对不同的卫星系统,可以根据系统参数结合式(27)求得对波束同步误差的限制。
表 3 波束同步误差对重叠幅宽的影响(m)Table 3. Influence of beam synchronization error to overlapping area (m)波束指向精度(°) 姿态控制精度(°) 0.01 0.02 0.01 260 411 0.02 411 520 0.03 582 663 0.04 759 823 0.05 938 991 0.06 1119 1164 6. 结论
本文分析了分布式SAR系统中时间同步误差和波束同步误差的影响,根据固定项和线性项的时间同步误差模型,推导时间同步和相位同步的关系,并分析线性项需要满足的条件。波束同步部分,本文分析了波束指向误差和卫星姿态误差对方位向合成天线方向图、多普勒去相干以及双基图像信噪比损失情况的影响,并根据几何关系推导距离向重叠幅宽的损失量。根据仿真参数,给出实验结果,以便在系统设计中充分考虑同步问题的影响,提出指标要求,为系统设计提供指导。
-
表 1 峨眉数据的参数
Table 1. Parameters of Emei data
参数 数值 波长 λ 0.031 m 斜距 r 1629 m 入射角 θ 47∘ 基线 b 0.2×11 m 瑞利分辨率 ρs 23.78 m -
[1] KNAELL K K and CARDILLO G P. Radar tomography for the generation of three-dimensional images[J]. IEE Proceedings-Radar, Sonar and Navigation, 1995, 142(2): 54–60. doi: 10.1049/ip-rsn:19951791 [2] 廖明生, 魏恋欢, 汪紫芸, 等. 压缩感知在城区高分辨率SAR层析成像中的应用[J]. 雷达学报, 2015, 4(2): 123–129. doi: 10.12000/JR15031LIAO Mingsheng, WEI Lianhuan, WANG Ziyun, et al. Compressive sensing in high-resolution 3D SAR tomography of urban scenarios[J]. Journal of Radars, 2015, 4(2): 123–129. doi: 10.12000/JR15031 [3] REIGBER A and MOREIRA A. First demonstration of airborne SAR tomography using multibaseline L-band data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(5): 2142–2152. doi: 10.1109/36.868873 [4] FORNARO G, SERAFINO F, and SOLDOVIERI F. Three-dimensional focusing with multipass SAR data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(3): 507–517. doi: 10.1109/TGRS.2003.809934 [5] BUDILLON A, EVANGELISTA A, and SCHIRINZI G. SAR tomography from sparse samples[C]. 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 2009: 865–868. doi: 10.1109/IGARSS.2009.5417514. [6] ZHU Xiaoxiang and BAMLER R. Super-resolution power and robustness of compressive sensing for spectral estimation with application to spaceborne tomographic SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(1): 247–258. doi: 10.1109/TGRS.2011.2160183 [7] VRIEZE S I. Model selection and psychological theory: A discussion of the differences between the Akaike information criterion (AIC) and the Bayesian information criterion (BIC)[J]. Psychological Methods, 2012, 17(2): 228–243. doi: 10.1037/a0027127 [8] BUDILLON A and SCHIRINZI G. GLRT based on support estimation for multiple scatterers detection in SAR tomography[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(3): 1086–1094. doi: 10.1109/JSTARS.2015.2494376 [9] ZHU Xiaoxiang, GE Nan, and SHAHZAD M. Joint sparsity in SAR tomography for urban mapping[J]. IEEE Journal of Selected Topics in Signal Processing, 2015, 9(8): 1498–1509. doi: 10.1109/JSTSP.2015.2469646 [10] FERRO-FAMIL L, HUANG Yue, and TEBALDINI S. Polarimetric characterization of 3-D scenes using high-resolution and Full-Rank Polarimetric tomographic SAR focusing[C]. 2016 IEEE International Geoscience and Remote Sensing Symposium, Beijing, China, 2016: 5694–5697. doi: 10.1109/IGARSS.2016.7730487. [11] TEBALDINI S, ROCCA F, D’ALESSANDRO M M, et al. Phase calibration of airborne tomographic SAR data via phase center double localization[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(3): 1775–1792. doi: 10.1109/TGRS.2015.2488358 [12] NANNINI M, SCHEIBER R, HORN R, et al. First 3-D reconstructions of targets hidden beneath foliage by means of polarimetric SAR tomography[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(1): 60–64. doi: 10.1109/LGRS.2011.2160329 [13] 柳祥乐. 多基线层析成像合成孔径雷达研究[D]. [博士论文], 中国科学院研究生院(电子学研究所), 2007.LIU Xiangle. Study on multibaseline tomography synthetic aperture radar[D]. [Ph. D. dissertation], Chinese Acedamy of Sciences, 2007. [14] 王金峰. SAR层析三维成像技术研究[D]. [博士论文], 电子科技大学, 2010.WANG Jinfeng. Study of three-dimentional synthetic aperture radar tomography imaging technology[D]. [Ph. D. dissertation], University of Electronic Science and Technology of China, 2010. [15] WANG Xiao, XU Feng, and JIN Yaqiu. Numerical simulation of tomographic-SAR imaging and object reconstruction using compressive sensing with L1/2-norm regularization[J]. Chinese Science Bulletin, 2014, 59(33): 4600–4607. doi: 10.1007/s11434-014-0554-5 [16] 刘康. 基于高分辨率SAR影像提取建筑物高度的研究[D]. [博士论文], 武汉大学, 2012.LIU Kang. Advanced methods for building height extraction from high resolution SAR images[D]. [Ph. D. dissertation], Wuhan University, 2012. [17] SUN Xilong, YU Anxi, DONG Zhen, et al. Three-dimensional SAR focusing via compressive sensing: The case study of angel stadium[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(4): 759–763. doi: 10.1109/LGRS.2011.2181321 [18] 李震, 张平, 乔海伟, 等. 层析SAR地表参数信息提取研究进展[J]. 雷达学报, 2021, 10(1): 116–130. doi: 10.12000/JR20095LI Zhen, ZHANG Ping, QIAO Haiwei, et al. Advances in information extraction of surface parameters using tomographic SAR[J]. Journal of Radars, 2021, 10(1): 116–130. doi: 10.12000/JR20095 [19] 仇晓兰, 焦泽坤, 彭凌霄, 等. SARMV3D-1.0: SAR微波视觉三维成像数据集[J]. 雷达学报, 2021, 10(4): 485–498. doi: 10.12000/JR21112QIU Xiaolan, JIAO Zekun, PENG Lingxiao, et al. SARMV3D-1.0: Synthetic aperture radar microwave vision 3D imaging dataset[J]. Journal of Radars, 2021, 10(4): 485–498. doi: 10.12000/JR21112 [20] JIAO Zekun, DING Chibiao, QIU Xiaolan, et al. Urban 3D imaging using airborne TomoSAR: Contextual information-based approach in the statistical way[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2020, 170: 127–141. doi: 10.1016/j.isprsjprs.2020.10.013 [21] ZHANG Fubo, LIANG Xingdong, CHENG Ruichang, et al. Building corner reflection in MIMO SAR tomography and compressive sensing-based corner reflection suppression[J]. IEEE Geoscience and Remote Sensing Letters, 2020, 17(3): 446–450. doi: 10.1109/LGRS.2019.2926301 [22] 李杭, 梁兴东, 张福博, 等. 基于高斯混合聚类的阵列干涉SAR三维成像[J]. 雷达学报, 2017, 6(6): 630–639. doi: 10.12000/JR17020LI Hang, LIANG Xingdong, ZHANG Fubo, et al. 3D imaging for array InSAR based on Gaussian mixture model clustering[J]. Journal of Radars, 2017, 6(6): 630–639. doi: 10.12000/JR17020 [23] KAY S M, 罗鹏飞, 张文明, 刘忠, 等译. 统计信号处理基础: 估计与检测理论[M]. 北京: 电子工业出版社, 2011: 451–453.KAY S M, LUO Pengfei, ZHANG Wenming, LIU Zhong, et al. translation. Fundamentals of Statistical Signal Processing Volume I: Estimation Theory[M]. Beijing: Publishing House of Electronics Industry, 2011: 451–453. [24] BAMLER R, EINEDER M, ADAM N, et al. Interferometric potential of high resolution spaceborne SAR[J]. Photogrammetrie - Fernerkundung - Geoinformation, 2009, 5: 407–419. doi: 10.1127/1432-8364/2009/0029 [25] CHEN S, BILLINGS S A, and LUO W. Orthogonal least squares methods and their application to non-linear system identification[J]. International Journal of Control, 1989, 50(5): 1873–1896. doi: 10.1080/00207178908953472 [26] DONOHO D L. Compressed sensing[J]. IEEE Transactions on Information Theory, 2006, 52(4): 1289–1306. doi: 10.1109/TIT.2006.871582 [27] MA Peifeng, LIN Hui, LAN Hengxing, et al. On the performance of reweighted L1 minimization for tomographic SAR imaging[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(4): 895–899. doi: 10.1109/LGRS.2014.2365613 [28] TAN Xing, ROBERTS W, LI Jian, et al. Sparse learning via iterative minimization with application to MIMO radar imaging[J]. IEEE Transactions on Signal Processing, 2011, 59(3): 1088–1101. doi: 10.1109/TSP.2010.2096218 [29] SCHMIDT R. Multiple emitter location and signal parameter estimation[J]. IEEE Transactions on Antennas and Propagation, 1986, 34(3): 276–280. doi: 10.1109/TAP.1986.1143830 [30] MOON T K. The expectation-maximization algorithm[J]. IEEE Signal Processing Magazine, 1996, 13(6): 47–60. doi: 10.1109/79.543975 [31] KRIEGER G, ROMMEL T, and MOREIRA A. MIMO-SAR tomography[C]. The 11th European Conference on Synthetic Aperture Radar, Berlin, Germany, 2016. -