系列文章目录
《SAR学习笔记》
目录
前言
上篇文章主要介绍了SAR成像中所用到的四个核心公式,本文将基于此介绍各种SAR成像算法,首先介绍后向投影算法(BPA)。
一、算法原理
前面介绍了SAR成像的本质是距离(快时间)、方位(慢时间)二维脉冲压缩。距离维(快时间)的脉冲压缩可以由发射的脉冲波形设计对应的匹配滤波器得以实现,对于方位维(慢时间维)脉冲压缩,BPA算法用了比较笨但是非常有效的方法,就是直接沿着距离徙动曲线做相干累加,具体步骤为沿着距离徙动曲线做相位补偿然后相干累加。因此,BPA算法核心就是对探测区域的每一个点,沿着距离徙动曲线,做相干累加。
二、算法步骤
根据BPA算法原理,对于的实现步骤如下图所示:
2.1 回波数据获取
对探测区域以脉冲重复间隔发射连续脉冲信号,并接收探测区域的回波信号。不失一般性,探测区域内任意一个散射目标的回波信号(已经下变频处理)为:
其中
其中为快时间采样间隔,为脉冲重复周期。
2.2 距离脉冲压缩
利用发射的脉冲信号设计距离向匹配滤波器,并用该滤波器对接收的每个脉冲信号匹配滤波,实现距离维的脉冲压缩。为了提高计算效率,通常在频域上处理。得到距离脉压后的信号为:
2.3 方位脉冲压缩
-
场景网格剖分
探测区域离散化处理:对探测的区域进行网格剖分,将探测区域剖分为个网格点
-
徙动曲线确定
对任意网格点,对应回波信号为:
其对应的距离徙动曲线为:
在已知情况下,该距离曲线是确定的,按照该距离徙动曲线在回波数据中找对应的徙动曲线位置,由于是离散采样点,是连续的徙动曲线,这样会导致,在不同慢时间下,找不到对应的离散快时间点。为减少这方面的误差,一般通过距离维的插值处理找到网格点对应距离徙动曲线上的回波信号为:
事实上,沿对应距离徙动曲线上的回波信号除了目标点的信号能量贡献,还包含其他目标回波在该曲线上的能量分布:
但其他目标的多普勒相位随慢时间的变化曲线与的不一致,所以这部分影响会在后续相位补偿相干累加中能够得到抑制。
-
多普勒相位补偿
平台的运动一方面使得距离维发生距离徙动现象,即不同方位(慢时间)同样的网格剖分点的距离不同。另一方方面平台的运动使得不同方位时间上同一个目标的脉冲回波叠加了一个瞬时距离变换引起的多普勒信号,由于每个脉冲时间内时间很短,一般不考虑快时间内的多普勒效应,慢时间的时间间隔较大,多普勒效应不能忽略,此时慢时间(方位)上多普勒效应就在距离曲线上得到体现,为了实现方位维的脉冲压缩,需要沿着距离先动曲线补偿对应的多普勒相位实现方位维的脉冲压缩。补偿的多普勒相位由精确的斜距信息得到:,沿距离徙动曲线补偿相位后的回波信号为:
-
相干累加成像
补偿距离徙动曲线位置处的多普勒相位,此时多普勒曲线上信号相位相干,将其累加即可得到:
2.4 SAR成像
可以看出得到的累加结果包含了的散射面积,通过方位脉冲压缩操作即可得到回波信号中其他散射点的散射面积,形成最终的SAR图像。
三、性能分析
3.1 计算效率
运算效率主要看进行复数乘法的次数。包含运算的步骤有:距离脉冲压缩、距离维插值运算以及距离徙动曲线相位补偿。相干累加只涉及复数加法运算,不做统计。
距离脉冲压缩:共进行了M次脉冲压缩,每次脉冲压缩进行了一次点FFT,一次点相乘,一次点IFFT,运算次数:
插值运算:假设在确定徙动曲线过程的插值方式采用频域补零的插值方式(不同网格点插值方式一致),插值倍数为,则插值运算为次点FFT,次点IFFT,运算次数:
相位补偿:共进行次点方位向相位补偿,运算次数:
3.2 适用场合
由于BPA实现方位向脉冲压缩过程并未用到近似处理(离散导致徙动曲线不在采样点上对于SAR成像算法是共性误差,不做对比考虑)。所以BPA能够用于包括大斜视场景,以及条带、聚束等模式中,唯一的缺点就是计算开销很大。
3.3 仿真结果
距离压缩结果
距离压缩插值前后对比
插值后的徙动曲线更平滑。
插值前后SAR图像:
插值后SAR图像能量旁瓣得到抑制,定位精度更高。
代码见SAR+BPA算法
总结
本节主要介绍后向投影(BPA)的SAR成像算法,在不考虑计算资源的情况下,该算法定位性能最佳,成像质量最好,适用范围最广。