
Citation: | Yu Yang, Pi Yi-ming. Terahertz High-resolution Radar Clutter Measurement and Analysis[J]. Journal of Radars, 2015, 4(2): 217-221. doi: 10.12000/JR14123 |
双基星载合成孔径雷达(Synthetic Aperture Radar, SAR)利用信号收发平台的分置,能够同时获取不同视角的观测数据,在测绘、干涉测量、地面目标识别、自然灾害监测等领域[1,2]具有重要的应用价值。以Tandem-L为代表的新一代双基星载SAR系统应用多通道、数字波束形成(Digital Beam Forming, DBF)等技术,实现方位向高分辨率和距离向宽测绘带SAR(High Resolution Wide Swath SAR, HRWS-SAR)成像,系统的成像能力得到显著改善[3,4]。
HRWS-SAR系统在方位向采用多通道采样降低系统的脉冲重复频率(Pulse Recurrence Frequency, PRF),在不降低方位分辨率的前提下实现系统的宽测绘带成像。HRWS-SAR系统进行成像处理时,若方位向信号为均匀采样可直接采用传统SAR成像处理方法,而受载星平台轨迹约束及应用场景的限制,方位向的非均匀采样更为普遍[5],因而对方位向信号进行重构以获得其均匀采样信号或多普勒谱是HRWS-SAR系统成像处理的一项关键技术。单基星载HRWS-SAR系统的方位向信号重构得到了广泛而深入的研究[6–17]。重排算法[6]和插值算法[7]是两种典型的时域重构算法。重排算法依据接收信号的方位向位置将各通道信号重新排列,只能在特定的PRF得到方位向信号的均匀采样。插值算法则依据文献[7]所提出的周期性非均匀采样信号重构公式通过时域插值得到方位向信号的均匀采样,该算法运算复杂度高且精度依赖于插值核的长度。Krieger等依据广义采样定理提出矩阵求逆算法[8,9],该算法通过线性方程求解从混叠的多通道信号中重构出无模糊的多普勒谱。该算法不需要协方差矩阵等先验信息且易于实现,但在重叠采样时该算法无法进行信号重构,接近重叠采样时重构性能也较差。文献[10–13]则采用不同的方法对矩阵求逆算法进行改进,实现重叠采样附近的高性能信号重构。文献[14,15]基于统一的信号模型,选取不同代价函数进行优化,提出了多种自适应波束形成(DBF)算法:正交投影算法、信号最大化算法、最大化信号模糊噪声比算法、最小均方误差算法等。与矩阵求逆算法相比,DBF类算法运算复杂度较高,并假定各通道的噪声为高斯白噪声且相互独立。文献[16]则对方位向信号重构算法的性能进行了仿真对比分析。文献[17]采用NUFFT直接重构方位向非均匀采样信号的多普勒频谱。
上述方位向信号重构方法理论上均能推广到双基星载HRWS-SAR系统。考虑到矩阵求逆算法易于实现且研究较为广泛,本文研究该算法在双基HRWS-SAR系统中的实现。本文首先将方位照射时间内时变的发射接收距离比近似为一个常数,推导了双基星载HRWS-SAR系统与某个单基系统在方位向多信道间传递函数上的等效性,从而构建了双基星载HRWS-SAR系统的方位向信号模型。而后,提出了适用于一般双基构型星载HRWS-SAR系统方位向信号重构的矩阵求逆算法,并给出了信噪比缩放因子及方位模糊比这两个重构性能指标的计算公式。最后,通过对几种典型双基构型的星载HRWS-SAR系统进行方位向信号重构仿真,验证了矩阵求逆算法在一般双基构型星载HRWS-SAR系统中的适用性。
为保持双基星载HRWS-SAR系统对地面的持续观测,收发平台应置于同一轨道(顺飞模式)或高度相同的平行轨道(平飞模式)。图1给出了一般构型双基星载HRWS-SAR系统的信号收发几何。发射天线和接收天线的最短距离分别为rT0, rR0,接收天线共有M个通道,RXequ为方位向信号重构后的等效接收通道,通道i到等效接收通道RXequ的方位向距离为
不失一般性,假定点目标P在t=0时刻位于接收天线的零多普勒面,容易得到t时刻点P在接收通道i的收发距离和
R(t)=√r2T0+(vt−vtfd)2+√r2R0+(vt−Δxi)2 | (1) |
忽略天线方向图的影响,接收通道i的冲激响应函数可表示为:
hi(t)=exp(−j2πλR(t))=exp(−j2πλ(√r2T0+(vt−vtfd)2+√r2R0+(vt−Δxi)2)) | (2) |
单基星载HRWS-SAR系统可视为rT0=rR0=r0, tfd=0,发射天线为等效接收通道的特殊情况。由式(2)可得此时接收通道i的冲激响应为:
˜hi(t)=exp(−j2πλ(√r20+(vt)2+√r20+(vt−Δxi)2)) | (3) |
令此时单基单通道SAR的冲激响应函数为
˜hmono(t)=exp(−j4πλ√r20+(vt)2) | (4) |
文献[9]对
{˜hmono(t)≈exp(−j4πλr0)⋅exp(−j2πv2t2λr0)˜hi(t)≈exp(−j4πλr0)⋅exp(−jπΔx2i2λr0) ⋅exp(−j2πv2(t−Δxi2v)2λr0) | (5) |
因而
˜hi(t)=exp(−jπΔx2i2λr0)⋅δ(t−Δxi2v)∗˜hmono(t) | (6) |
将式(6)进行Fourier变换,可得单基单通道SAR与接收通道i的系统函数关系:
˜Hi(f)=exp(−jπΔx2i2λr0)⋅exp(−j2πΔxi2vf)⋅˜Hmono(f) | (7) |
单基单通道SAR到接收通道i的传递函数
˜Gi(f)=exp(−jπΔx2i2λr0)⋅exp(−j2πΔxi2vf) | (8) |
基于上述分析,单基星载HRWS-SAR系统的方位向信号生成模型可由图2进行描述。方位向各通道信号在时域为单基单通道SAR信号经对应的相位偏移和时间延迟的结果,而在多普勒域为单基单通道SAR信号与对应的常数相位因子和线性相位因子相乘后的结果。PRF小于多普勒带宽时,方位向各通道信号为时域上的欠采样,将导致多普勒域的频谱混叠。
一般双基构型的星载HRWS-SAR系统中rT0≠rR0, tfd≠0,不能通过对hmono及hi的泰勒级数展开进行2阶相位近似处理得到与式(6)、式(7)相似的表达式。本节将对式(2)进一步处理,推导一般双基构型星载HRWS-SAR系统中
引入点目标P的收发距离比函数
hi(t)=exp(−j2πλ(C(t)√r2R0+(vt)2+√r2R0+(vt−Δxi)2)) | (9) |
等效接收通道RXequ的冲激响应函数为:
hmono(t)=exp(−j2πλ((C(t)+1)√r2R0+(vt)2)) | (10) |
对式(9)、式(10)的泰勒级数展开进行2阶相位近似处理,可得:
{hi(t)≈exp(−j2πλ(C(t)+1)rR0) ⋅exp(−jπλrR0⋅C(t)Δx2iC(t)+1) ⋅exp(−jπλrR0⋅(C(t)+1) ⋅v2(t−Δxi(C(t)+1)v)2)hmono(t)≈exp(−j2πλ(C(t)+1)rR0) ⋅exp(−jπλrR0⋅(C(t)+1)⋅v2t2) | (11) |
在目标照射时间Ta内,
hi(t)=exp(−jπλrR0⋅C0Δx2iC0+1) ⋅δ(t−Δxi(C0+1)v)∗hmono(t) | (12) |
式(12)进行Fourier变换,可得等效接收通道RXequ与接收通道i的系统函数关系:
Hi(f)=exp(−jπλrR0⋅C0Δx2iC0+1)⋅exp(−j2πΔxi(C0+1)vf)⋅Hmono(f) | (13) |
等效接收通道RXequ到接收通道i的传递函数为:
Gi(f)=exp(−jπλrR0⋅C0Δx2iC0+1)⋅exp(−j2πΔxi(C0+1)vf) | (14) |
考察式(8)、式(14),传递函数
Gi(f,Δxi)=˜Gi(f,2(C0+1)Δxi) | (15) |
式(15)表明双基星载HRWS-SAR系统的方位向信号的生成模型可等效为一单基系统模型,只需将图2中的
式(12)、式(13)、式(14)是在
文献[8]给出了单基HRWS-SAR系统方位向信号重构的矩阵求逆算法,由于双基HRWS-SAR系统在方位向信号模型上与单基系统的等效性,矩阵求逆算法在双基系统同样适用。
矩阵求逆算法的理论依据是广义采样定理。在双基HRWS-SAR系统中,通过发射到接收通道i的传递函数,能够得到单基单通道SAR的信号
矩阵求逆算法主要包括3个步骤:
首先,构造系统的传递函数矩阵
G(f)=(G1(f)G1(f+PRF)···G1(f+(M−1)⋅PRF)G2(f)G2(f+PRF)···G2(f+(M−1)⋅PRF)⋮⋮⋱⋮GM(f)GM(f+PRF)···GM(f+(M−1)⋅PRF)) | (16) |
˜S(f)=(S1(f),S2(f),⋯,SM(f))T | (17) |
然后,对
P(f)=(P1(f)P2(f)···PM(f)P1(f+PRF)P2(f+PRF)···PM(f+PRF)⋮⋮⋱⋮P1(f+(M−1)⋅PRF)P2(f+(M−1)⋅PRF)···PM(f+(M−1)⋅PRF)) | (18) |
最后,利用公式
星载HRWS-SAR系统的单通道信号带宽为PRF,采用矩阵求逆算法重构后的方位向信号带宽为
3.2.1 信噪比缩放因子 文献[8]定义信噪比缩放因子表征方位向重构对信噪比的影响,其定义式及计算公式为:
Φbf=SNRin/SNRout(SNRin/SNRout)|PRFuni=M∑j=1E(|Pj(f)|2)N |
(19)
其中,SNRin, SNRout为方位向信号重构前后的信噪比,PRFuni表示方位向均匀采样时的PRF,
Φbf=M∑j=1λj(f)N | (20) |
其中,
3.2.2 方位模糊比 文献[8]详细推导了单基HRWS-SAR系统的方位模糊比(AASR)计算公式,与之类似可推导出双基HRWS-SAR系统的AASR计算公式,可概括为如下4个计算式:
ek(f)=Ak(f)∑mM∑j=1Gj(f+k⋅PRF) ⋅Pj(f+m⋅PRF) | (21) |
eΣ(f)=∞∑k=−∞,k≠0ek(f)=∞∑k=−∞,k≠0(Ak(f)∑mM∑j=1Gj(f+k⋅PRF) ⋅Pj(f+m⋅PRF)) | (22) |
ps=E(|A(f)⋅rect(f/I)|2) | (23) |
AASR=E(|eΣ(f)|2)ps | (24) |
本节对几种典型的双基构型星载HRWS-SAR系统进行方位向信号重构进行仿真,验证矩阵求逆算法的正确性,并分析算法的重构性能。表1列出了系统的方位向系统参数,对于单基系统容易计算出多普勒带宽Bfd=5.61 kHz,照射时间Ta=1.05 s。
参数 | 数值 |
发射天线方位尺寸(m) | 2.4 |
接收通道方位尺寸(m) | 2.4 |
接收通道数目 | 5 |
轨道高度(km) | 600 |
接收天线最短距离(km) | 700 |
方位向速度(m/s) | 7600 |
载波波长(cm) | 3.1 |
表2列出了7种典型的双基构型。表2中采用发射天线和接收天线的零多普勒时间差tfd与轨道距离L对双基构型进行表征:tfd=0时系统工作在平飞模式,L=0时系统工作在顺飞模式,L<0时发射天线位于测绘带近端,反之L>0时接收天线位于测绘带近端,构型Ⅰ中tfd=0且L=0退化为单基系统。
构型编号 | tfd (s) | L (km) |
Ⅰ | 0 | 0 |
Ⅱ | 1 | 0 |
Ⅲ | 10 | 0 |
Ⅳ | 0 | 10 |
Ⅴ | 0 | 100 |
Ⅵ | 0 | –10 |
Ⅶ | 0 | –100 |
将
各接收通道的信噪比设置为20 dB,取PRF=2.0 kHz,对图1所示的点目标P方位向信号重构过程进行仿真,其结果如图4和图5所示。
图4(a)、图4(b)中,由于PRF小于多普勒带宽Bfd,单个接收通道的多普勒频谱有严重的混叠现象,图4(c)则表明采用矩阵求逆算法进行方位向信号重构消除了多普勒频谱中的混叠现象。图4(c)及图5给出了各双基构型下的方位向信号重构后的成像结果,验证了矩阵求逆算法对一般双基构型的适用性,同时也表明重构性能受双基构型的影响,尤其是构型Ⅴ、Ⅶ与其他构型性能差异较为明显。
将PRF设置在区间1.4 kHz≤PRF≤2.8 kHz,采用表1的系统参数和表2的双基构型,利用式(20)、式(24)可以得到方位向信号的重构性能曲线,如图6所示。为了便于图6的分析,首先计算各双基构型下的C0及典型PRF值,计算结果为表3,其中PRFuni为均匀采样时的PRF值,PRFrep1及PRFrep2为PRF范围内出现重叠采样的两个PRF值。构型Ⅰ的PRFuni, PRFrep1和PRFrep2在图6中予以标注。
构型编号 | C0 | PRFuni (kHz) | PRFrep1 (kHz) | PRFrep2 (kHz) |
Ⅰ | 1.0000 | 2.533 | 1.583 | 2.111 |
Ⅱ | 1.0001 | 2.533 | 1.583 | 2.111 |
Ⅲ | 1.0059 | 2.540 | 1.588 | 2.117 |
Ⅳ | 0.9927 | 2.524 | 1.577 | 2.103 |
Ⅴ | 0.9345 | 2.450 | 1.531 | 2.041 |
Ⅵ | 1.0074 | 2.542 | 1.589 | 2.118 |
Ⅶ | 1.0805 | 2.635 | 1.647 | 2.196 |
(1) 均匀采样时PRF=PRFuni,方位向信号的重构性能达到局部最优;重叠采样时PRF=PRFrep,信噪比缩放因子趋向于无穷大,矩阵求逆算法将不能实现信号重构;PRF位于PRFrep附近时,矩阵求逆算法的性能急剧下降。从矩阵理论的观点分析,PRF=PRFrep时,传递函数
(2) C0值的差异直接导致同一PRF时方位向信号重构性能差异。与其他构型相比构型Ⅴ、Ⅶ的C0值差异较大因而构性能差异也大,与图5的结论相一致。根本原因在于与双基HRWS-SAR系统等效的单基系统,通道间的方位向距离为
本文通过分析单基与双基星载HRWS-SAR系统的方位向信号模型,给出了适用于一般双基构型星载HRWS-SAR系统方位向信号重构的矩阵求逆算法,并使用信噪比缩放因子及方位模糊比两个指标分析了该算法的重构性能。本文的分析方法对其他方位向信号重构算法推广到双基星载HRWS-SAR系统具有借鉴意义。双基星载HRWS-SAR系统方位向信号重构的工程应用中,需改进矩阵求逆算法以改善重叠采样附近的重构性能。
[1] |
Oliver C J. Optimum texture estimators for SAR clutter[J]. Journal of Physics D, Applied Physics, 1993, 26(11): 1824-1835.
|
[2] |
Li H C, Hong W, Wu Y R, et al.. On the empiricalstatistical modeling of SAR images with generalized gamma distribution[J]. IEEE Journal of Selected Topics in Signal Processing, 2011, 5(3): 386-397.
|
[3] |
Frery A C, Muller H J, Yanasse C C F, et al.. A model for extremely heterogeneous clutter[J]. IEEE Transactions on Geoscience and Remote Sensing, 1997, 35(3): 648-659.
|
[4] |
Cui S and Datcu M. Coarse to fine patches-based multitemporal analysis of very high resolution satellite images[C]. Proceeding of 6th International Workshop on the Analysis of Multi-Temporal Remote Sensing Images, Trento, Italy, 2011: 85-88.
|
[5] |
Oliver C J and Quegan S. Understanding Synthetic Aperture Radar Images[M]. Norwood, MA, USA: Artech House, 1998: 128-130.
|
[6] |
Nicolas J M. Introduction to second kind statistics: application of logmoments and log-cumulants to second kind statistics: application of logmoments and log-cumulants to analysis of radar images[J]. Traitement du Signal, 2002, 19(3): 139-167.
|
[7] |
Krylov V and Zerubia J. Generalized gamma mixtures for supervised SAR image classification[C]. Proceedings of GraphiCon 2010, Saint-Pertersburg, Russia, 2010: 107-110.
|
[8] |
Tadjudin S and Landgrebe D A. Robust parameter estimation for mixture model[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(1): 439-445.
|
[9] |
Moser G, Zerubia J, and Serpico S B. SAR amplitude probability density function estimation based on a generalized Gaussian model[J]. IEEE Transactions on Image Processing, 2006, 15(6): 1429-1442.
|
[10] |
Szajnowski W. Estimator of log-normal distribution parameters[J]. IEEE Transactions on Aerospace and Electronic Systems, 1977, AES-13(5): 533-536.
|
[11] |
Stacy E W. A generalized of the Gamma distribution[J]. The Annals of Mathematical Statistics, 1962, 33(3): 1187-1192.
|
1. | 叶恺,禹卫东,王伟. 基于矩阵束方法的星载MEB SAR俯仰向DBF处理方法. 电子与信息学报. 2018(11): 2659-2666 . ![]() |
参数 | 数值 |
发射天线方位尺寸(m) | 2.4 |
接收通道方位尺寸(m) | 2.4 |
接收通道数目 | 5 |
轨道高度(km) | 600 |
接收天线最短距离(km) | 700 |
方位向速度(m/s) | 7600 |
载波波长(cm) | 3.1 |
构型编号 | tfd (s) | L (km) |
Ⅰ | 0 | 0 |
Ⅱ | 1 | 0 |
Ⅲ | 10 | 0 |
Ⅳ | 0 | 10 |
Ⅴ | 0 | 100 |
Ⅵ | 0 | –10 |
Ⅶ | 0 | –100 |
构型编号 | C0 | PRFuni (kHz) | PRFrep1 (kHz) | PRFrep2 (kHz) |
Ⅰ | 1.0000 | 2.533 | 1.583 | 2.111 |
Ⅱ | 1.0001 | 2.533 | 1.583 | 2.111 |
Ⅲ | 1.0059 | 2.540 | 1.588 | 2.117 |
Ⅳ | 0.9927 | 2.524 | 1.577 | 2.103 |
Ⅴ | 0.9345 | 2.450 | 1.531 | 2.041 |
Ⅵ | 1.0074 | 2.542 | 1.589 | 2.118 |
Ⅶ | 1.0805 | 2.635 | 1.647 | 2.196 |
参数 | 数值 |
发射天线方位尺寸(m) | 2.4 |
接收通道方位尺寸(m) | 2.4 |
接收通道数目 | 5 |
轨道高度(km) | 600 |
接收天线最短距离(km) | 700 |
方位向速度(m/s) | 7600 |
载波波长(cm) | 3.1 |
构型编号 | tfd (s) | L (km) |
Ⅰ | 0 | 0 |
Ⅱ | 1 | 0 |
Ⅲ | 10 | 0 |
Ⅳ | 0 | 10 |
Ⅴ | 0 | 100 |
Ⅵ | 0 | –10 |
Ⅶ | 0 | –100 |
构型编号 | C0 | PRFuni (kHz) | PRFrep1 (kHz) | PRFrep2 (kHz) |
Ⅰ | 1.0000 | 2.533 | 1.583 | 2.111 |
Ⅱ | 1.0001 | 2.533 | 1.583 | 2.111 |
Ⅲ | 1.0059 | 2.540 | 1.588 | 2.117 |
Ⅳ | 0.9927 | 2.524 | 1.577 | 2.103 |
Ⅴ | 0.9345 | 2.450 | 1.531 | 2.041 |
Ⅵ | 1.0074 | 2.542 | 1.589 | 2.118 |
Ⅶ | 1.0805 | 2.635 | 1.647 | 2.196 |