汪 玲 孙玲玲 宫 蕊 朱岱寅
(南京航空航天大学电子信息工程学院微波光子与雷达成像教育部重点实验室 南京 211000)
大气折射率的实部和虚部分别表示大气减缓电磁波速率和吸收电磁能量的能力。雷达信号处理时通常认为发射波和散射波都是在真空中传播的,默认大气折射率为1,也就是认为电磁波在大气中的传 播速率等于自由空间光速且忽略大气吸收引起的电磁波衰减。然而,真实的传播过程中,电磁波速率会因为湍流、风场等不可抗因素而产生波动,大气吸收会导致电磁波能量衰减,引起折射率波动,从而导致合成孔径雷达(Synthetic Aperture Radar,SAR)成像误差。
Porcello教授[1](1970年)在对大气传播产生的SAR成像影响的研究中指出低层大气中湍流导致的折射率波动会限制大型天线孔径和阵列的性能。Danklmayer等人[2](2009年)将SAR图像与几乎同时刻获得的天气雷达数据进行比较,得出的结论是50 dBZ的气象雷达图像的反射率可能会导致SAR图像中的可见伪影,此外还讨论了对流层对TerraSAR-X外部定标的影响。对流层和电离层对星载SAR系统的影响的研究结果表明,未校正的电离层延迟差异将会导致图像散焦[3]。大气状况的变化对SAR成像的影响分为两类,部分研究学者将图像伪影归因于大气不均匀性[4–7],部分则认为是由于入射到场景上的功率发生变化而导致了伪影的产生[8]。在大气折射率的时间变化对圆迹SAR聚焦性能的影响[9]中,通过建立相位误差模型推导出了折射率时间变化对圆迹SAR方位向聚焦性能的影响。
上述现有的关于大气传播对雷达成像的影响的研究的局限性在于只是提出了影响因素和影响在图像上的表现形式,没有给出大气折射率波动和图像误差之间的定量解析关系。针对此局限性,本文将Yazici等人[10](2019年)研究工作的基础上对SAR成像误差展开研究,定量分析大气折射率误差(包括电磁波速率误差和大气吸收)对SAR成像的影响,并通过仿真实验来验证理论分析。
本文的结构安排如下:第2节在考虑电磁波速率误差和电磁波衰减的前提下建立回波模型,然后在单位折射率条件下使用滤波反投影(Filtered Back Projection, FBP)算法进行图像重建,得到复折射率影响下的成像点扩散函数误差。第3节定量推导出电磁波速率误差引起的成像位置误差的解析式,分析电磁吸收导致的幅度误差。第4节通过仿真实验来验证误差分析的正确性。第5节总结全文并得出结论。
X为待成像场景中的一散射点,SAR回波的时延和多普勒频移分别为
图1 地形平坦时单基SAR中向量和角度示意图
由式(20)可以得出以下结论:
(1)电磁波在大气中的真实速率越低,即电磁波速率误差∆ c越大,径向和横向位置误差越大。
(2)天线与成像目标的距离R 越大,径向和横向位置误差越大。
(3)擦地角η 和 η⊥越大,径向和横向位置误差越大。
(4)雷达观测角φ 仅影响横向位置误差。观测角越小或斜视角越大,横向位置误差越大。当φ ≈90◦时,即雷达的速度方向与视线方向垂直时,电磁波速率误差引起的横向位置误差可以忽略不计。
(5) α , α⊥的增加,意味着x 处的地形变得陡峭,径向和横向位置误差将会变小,因此,电磁波速率误差导致的平原地区(如沙漠)SAR图像的位置误差大于山区SAR图像的位置误差。
(6)位置误差是依赖于 s 和x 的,这说明位置误差具有时空变化特征。
(7)小场景和合成孔径相对较短的情况下,可以认为式(20)中的角度和距离是不随着慢时间 s和散射点位置 x变化的,可以用场景中心来代替整个场景,合成孔径中心来代表整个合成孔径。如果成像过程中的电磁波速率误差是恒定的,即∆ c为一常数,那么∆ zr和 ∆ zt在整个孔径和场景中也保持恒定。
为了能够准确分析大气吸收的影响,本节中不考虑电磁波速率误差,即假设电磁波在大气中的真实速率等于真空中光速[20,21],点目标T(x)=δ(x−x0)的重建反射率可表示为
由式(30)看出,双基距离越大, a 越小,滤波器的截止频率就越小,频谱衰减越快。
此等效滤波器的频谱特性如图2(a)所示,滤波器使用到的参数分别为中心频率f0=9.6 GHz,双程斜距 R=100 km,复折射率虚部ni=10−4,蓝色竖线为截止频率,棕色曲线的走势表明吸收引起的滤波会引起幅度呈现指数衰减,且此滤波器幅度数量级为1 0−135,也就是说幅度衰减将会非常严重。
对线性调频信号(Linear Frequency Modulation,LFM)使用该滤波器来说明折射率虚部的影响,如图2(b)所示,红色曲线为滤波前幅度为1的LFM信号的时域压缩波形,蓝色曲线为滤波后波形。可以看出,滤波器增加了主瓣宽度和副瓣电平,因此图像分辨率和对比度降低。
由于此滤波器还受双程斜距的影响,因此距离不同时引起重建图像中散射体强度的变化也是不一样的。对LFM信号使用两个距离参数不同的滤波器来滤波,如图2(c)所示,红色和蓝色波形分别为使用 R=99.9 km和R =100 km的滤波器进行滤波后的时域信号波形,可以看到两个信号的幅度是不一样的,双基距离越大,信号幅度越小,因此可以得出大气吸收影响下,目标离雷达越远,其重建反射 率就越弱的结论。
图2 吸收引起的滤波效应
本节中将进行仿真实验以验证前文分析的正确性,考虑两种平坦成像场景,第1个为单点目标,此点目标位于场景中心。第2个场景是合成场景,包含单词“NUAA”,为多点目标,所有散射点反射率均设置为1。使用LFM信号作为发射波形,带宽为62.5 MHz,载波频率为f0=700 MHz。
使用单基线性轨迹SAR分别对单点目标和多点目标进行误差仿真,天线运动轨迹为γ (s)=vs+H,其中 速 度 v =[69.4,0,0] m/s,初始 位置H =[0,0,7300] m,场景中心位于[7100 Rc] m处, Rc的大小反映的是天线和目标之间的距离。
本节中假设电磁波速率等于真空中光速且不存在吸收,即折射率为1,成像结果如图3所示,可以看到点目标和目标“NUAA”在正确位置(图中用红色方框标出)重建,且聚焦良好。
本节中假设电磁波真实传播速率不等于真空中光速,而成像过程中忽略了速率误差,即认为n=1。为了验证电磁波速率误差大小以及天线到目标的距离大小对重建位置误差的影响,此节将在3种配置对单点目标进行误差仿真,配置的参数表1所示,成像结果如图4所示,图4中红色圆圈表示点目标的正确位置。
从这3幅图中可以发现,电磁波速率误差将会导致目标重建位置相对正确位置产生一个恒定偏移,而偏移主要集中在径向视线方向,横向上基本为零。对比图4(a)和图4(b),可以发现图4(b)中的位置误差较大,是由于配置II中的电磁波速率误差比配置I中的电磁波速率误差大。对比图4(a)和图4(c),发现图4(c)中的位置误差较小,这是因为配置III中天线到目标的距离比配置I中的距离小。由以上3幅图的对比可以得出结论:天线与目标的距离不变时,电磁波速率误差越大,重建位置误差越大;电磁波速率误差恒定时,天线与目标的距离越远,图像位置误差越大。此结论与前文折射率误差导致的成像位置误差的理论分析结果一致。
图3 无折射率误差时的单点和多点目标重建图
表1 电磁波速率误差仿真参数设置
图5为各个配置下点目标成像位置误差的理论值和仿真实验测量值曲线,其中图5(a)和图5(b)分别为配置I条件下的小孔径和大孔径成像位置误差曲线,图5(c)和图5(d)为配置II和配置III条件下的大孔径成像位置误差曲线。图中蓝色细实线和蓝色虚点线分别表示径向位置误差的理论值和测量值,红色虚划线和点划线分别表示横向位置误差的理论值和测量值。图5中各配置下的成像位置误差的理论值曲线和实测值曲线基本重合,这说明了电磁波速率误差对SAR成像的影响的理论分析是正确的。图5(a)中的径向位置误差在整个成像过程中保持恒定,这是因为合成孔径时间较短导致的,当合成孔径时间较长时,误差曲线便会随着慢时间缓慢变化,如图5(b)、图5(c)、图5(d)所示。由图5(b)、图5(c)、图5(d)中的径向位置误差的理论值曲线可知,3种配置下径向位置误差的绝对值约为30 m,49 m和16 m,对比配置I和配置II,可得电磁波速率误差越大,径向位置误差越大,对比配置I和配置III,可得天线离目标越远,径向位置误差越大,这与从电磁波速率误差影响下的重建图像(图4)中得到的结论是一致的。
图4 电磁波速率误差影响下的点目标重建图
图5 波速误差导致的SAR图像位置误差的理论值和测量值
本节中假设回波是经过电磁吸收后的信号,即大气折射率为复数,而成像是认为不存在吸收的前提下进行的,以此来验证吸收带来的影响。
图6 电磁波速率误差影响下的多目标场景“NUAA”重建图
首先对点目标进行仿真,假设大气复折射率为n=1+6×10−5i ,将天线到目标的距离Rc的值分别设置为 30 和6 0 km,来验证距离大小对幅度误差的影响,SAR成像结果如图7所示。与图3(a)中的正确重建相比,可以发现目标重建位置是正确的,但是目标能量不再聚焦,而是往径向扩散,目标副瓣振荡明显,这就是大气吸收对图像“滤波”作用的表现。对比图7(a)和图7(b),当天线与成像场景的距离变远时,大气吸收带来的影响将会变大,图像散焦现象更加严重,这与2.2节中得到的大气吸收对距离较远处的散射体影响更大的结论是一致的。
吸收作用下的场景“NUAA”的重建图像如图8所示,可以发现目标在正确的位置重建,但“NUAA”能量沿着径向扩散,图像严重散焦。
本文研究了大气折射率误差对SAR成像的影响,分析得出,电磁波速率误差引起的折射率实部误差会导致散射点在重建图像中的位置误差。推导得出两者之间的定量关系,得知位置误差取决于电磁波速率误差大小以及天线到目标的距离。吸收效应引起的电磁衰减会使点扩散函数产生额外幅度项,对PSF产生“滤波”作用,滤波后的波形呈指数衰减,衰减程度取决于双基距离和吸收系数,将导致散射点重建反射率的变化,表现为重建图像的幅度误差。仿真实验验证了误差理论分析的正确性,电磁波速率误差引起的成像位置误差在仿真中的实测值与理论分析值一致,忽略吸收时的重建图像沿径向散焦,场景反射率发生变化。
本文中所完成的SAR成像的误差定量分析工作有助于SAR图像正确解译,也是后续SAR自聚焦技术研究的支撑,所使用的基本分析方法也可以用于其他因素导致的SAR成像误差分析,除此之外,还可延伸到其他领域,如声呐成像、医学图像伪影分析等。
图7 大气吸收作用下的点目标重建图
图8 大气吸收作用下的多目标场景“NUAA”重建图