各种SAR成像算法总结(共21页).doc

精选优质文档-倾情为你奉上各种SAR成像算法总结1 SAR成像原理SAR成像处理的目的是要得到目标区域散射系数的二维分布，它是一个二维相关处理过程，通常可以分成距离向处理和方位向处理两个部分。在处理过程中，各算法的区别在于如何定义雷达与目标的距离模型以及如何解决距离方位耦合问题，这些问题直接导致了各种算法在成像质量和运算量方面的差异。一般来说，忽略多普勒频移所引起的距离向相位变化，距离向处理变为一维的移不变过程且相关核已知，即退化为一般的脉冲压缩处理；同时将雷达与目标的距离按2阶Taylor展开并忽略高次项，则方位向处理也是一个一维的移不变过程，并退化为一般的脉冲压缩处理，这就是经典的距离多普勒（Range-Doppler RD）算法的实质。若考虑多普勒频移对距离向相位的影响，同时精确的建立雷达与目标的距离模型，则不论距离向处理还是方位向处理都变为二维的移变相关过程。线性调频尺度变换（Chirp-Scaling CS）算法即在此基础之上将二维数据变换到频域，利用Chirp Scaling原理及频域的相位校正方法，对二维数据进行距离徙动校正处理、距离向及方位向的聚焦处理，最终完成二维成像处理。当方位向数据积累延迟小于全孔径时间（即方位向为子孔径数据）的情况下，方位向处理必须使用去斜（dechirp）处理及频谱分析的方法。在RD和CS算法的基础之上，采用dechirp处理及频谱分析的方法完成方位向处理的算法分别称为频谱分析（SPECAN）算法和扩展CS（Extended Chirp Scaling ECS）算法。1.1 SAR成像原理本节以基本的正侧视条带工作模式为例，对SAR的成像原理进行分析和讨论。正侧视条带SAR的空间几何关系如下图所示。图中，o平面为地平面，o垂直于o平面。SAR运动平台位于S点，其在地面的投影为G点。SAR运动平台的运动方向Sx平行于o，速度大小为。SAR天线波束中心与地面的交点为C，CG与运动方向Sx垂直；S与C的距离为，称为天线波束的方位向宽度，大小为。P为测绘带内的某一点，一般情况下取斜距平面CSP进行分析，称SAR运动的方向Sx为方位向（或方位维），称天线波束指向方向SC为距离向（或距离维）。正侧视条带SAR几何关系示意图假定P的方位向坐标为；在时刻，SAR运动平台S与P的距离为。若当时刻，SAR运动平台位于方位向0点，则当时刻，的表达式为：(1.1)将式(1.1)在附近进行2阶Taylor展开，有： (1.2)假设雷达发射连续的正弦波，即发射信号为：(1.3)其中，为发射正弦波的幅度，为发射信号的载频。发射信号经点目标P散射后，雷达接收机收到的信号为：(1.4)其中：为光速，为复常数，为回波信号相对于发射信号的时间延迟：(1.5)为考虑雷达水平方向增益变化而引入的加权函数。若不考虑雷达天线的加权作用，即令，则式(1.4)变为：(1.6)根据式(1.6)，雷达运动平台相对于点目标的运动将造成回波信号的相位随时间不断变化，从而引起回波瞬时频率的变化，产生多普勒频移。多普勒频移量为：(1.7)将式(1.2)内的代入可得：(1.8)其中：为雷达工作波长，且，为雷达波束中心通过P点的时间。回波信号的瞬时频率为：(1.9)由式(1.9)可知，多普勒频移的存在将使回波信号的瞬时频率在载波频率附近作线性变化。也就是说，由于雷达运动平台匀速直线前进，回波信号在方位向将为线性调频（chirp）信号：(1.10)其中为固定相位项，略去后，式(1.10)可简化为：(1.11)通常为便于对回波信号进行处理，需要将回波信号经频率变换调至较低频率，回波多普勒频率将以为中心变化。中心频率称为偏置频率。因此有：(1.12)式中表示回波信号经变频处理将载频降至偏置频率后的瞬时频率变化。通常称它为点目标回波信号的多普勒频率历史，简称多普勒历史。由式(1.12)可见，多普勒历史是一按负斜率变化的chirp信号，其调频斜率为：(1.13)即点目标回波信号的调频斜率与成正比、与成反比。点目标横过波束的最大距离称为合成孔径长度，其大小与以及方位向波束宽度有关；点目标横过波束的时间称为合成孔径时间。有：(1.14)(1.15)在合成孔径时间里，多普勒频率的变化范围称为多普勒带宽，用表示。由式(1.14)、(1.15)得到的表达式为：(1.16)考虑到对于方位向天线直径为的天线，近似有：(1.17)因此，SAR的方位向理论分辨率为：(1.18)从上述分析可以看出，由于雷达运动平台作等高匀速直线运动，使得目标的回波信号在方位向上具有线性调频特性，对回波信号进行脉冲压缩处理，可以获得方位向的高分辨率。在理想情况下，SAR方位向分辨率与雷达平台的速度、飞行高度、作用距离、雷达工作波长等参数无关，只与天线尺寸有关，为天线方位向口径尺寸的一半，这是SAR的一大特点和优势。1.2 SAR回波信号模型1.1节分析了SAR成像的基本原理，本节推导SAR回波信号的数学模型，给出SAR信号处理的理论基础。chirp信号是SAR系统中最常用的发射信号形式。假设雷达发射的chirp脉冲串为：(1.19)其中，为发射信号的包络，为chirp信号的调频斜率，为发射信号脉宽，为发射信号的中心频率，为脉冲重复周期。则雷达于时刻，接收到斜距为处目标反射的回波信号为： (1.20)其中，为目标的后向散射特性，为方位向的天线方向性函数，为光速。经正交解调后的复信号可以表示为： (1.21)其中，为雷达工作波长。式(1.21)中的两个指数项分别代表方位向的相位调制和距离向发射的相位调制。考虑到相对于雷达发射脉冲而言，和是时间的慢变化函数，可以作如下近似：(1.22)(1.23)同时，将时间分解为快时间分量和慢时间分量之和，即：，(1.24)通过变量置换，可以将转换成二维形式： (1.25)其中，表示对的卷积，表示对的卷积，表示二维卷积。因此，雷达系统接收回波信号的过程，可以看作是地面目标的后向散射特性通过一个线性系统的过程。式(1.25)可简化表示为：(1.26)其中，为线性系统的冲激响应函数： (1.27)式(1.27)可以进一步表示为：(1.28)其中，(1.29)(1.30)则式(1.26)可进一步表示为：(1.31)式(1.29)中，的指数项代表了由于雷达运动平台与目标间相对运动所带来的方位向相位调制。如果对采用式(1.2)所示的2阶Taylor展开方式，则回波的方位向相位为慢时间的2次函数，即一个chirp信号；的冲击函数表达式代表了由于相对运动，回波包络的中心在距离向上的位置发生变化，即距离徙动现象。2 RD算法原理RD算法流程如下图所示，包括距离压缩处理、方位压缩处理两个主要处理步骤，以及作为辅助处理步骤的距离徙动校正处理。由于具有概念简单、易于实现、处理效率高等优点，RD算法成为最经典、最成熟的SAR成像处理算法。RD算法流程RD算法的本质是对采用式(2.2)所示的2阶Taylor展开方式，将距离向处理和方位向处理解耦，分解为两个一维处理分别完成。其中距离向处理利用脉冲压缩技术实现距离向高分辨，方位处理则利用回波中的多普勒信息完成方位高分辨。2.1 RD算法的距离向处理SAR回波信号的表达式为：(2.1)其中，(2.2)(2.3)由于为chirp 信号，距离向处理就是针对完成匹配滤波处理。选取距离向处理参考函数：(2.4)则距离向处理后的信号近似为： (2.5)其中，为距离向处理结果的包络，当发射信号的包络为门函数时：(2.6)为sinc函数：(2.7)其中为发射信号的带宽，为发射脉冲宽度。一般情况下为了获得距离向的高分辨，发射脉冲的带宽很大，此时近似为函数。2.1.1 距离徙动校正处理将距离向处理结束后的信号重写如下： (2.8)由于在不同的慢时间，雷达和目标的距离不同，因此式(2.47)中距离向处理结果包络的最大值随慢时间的变化出现在不同的距离向位置上，这种现象称为距离徙动现象。距离徙动现象的本质是回波信号的方位向和距离向发生耦合，如果要进行精确成像，方位向就需要进行二维相关处理。为了使信号的方位向与距离向解耦，从而简化方位处理，使之变为一维相关处理，就需要在方位向处理之前进行距离徙动校正，使式(2.8)变为如下形式：(2.9)其中，为不随慢时间变化的参考距离。对(2.8)中的斜距按二阶Taylor展开，有：(2.10)式(2.8)可以改写为如下形式： (2.11)其中，为回波方位向多普勒调频斜率。处于不同方位向位置的点目标，其距离徙动变化曲线各不相同。在实际处理过程中，必须针对不同方位向位置逐一进行距离徙动校正处理。为了简化距离徙动校正处理，减小处理量，可以利用方位向回波chirp信号的时频关系： (2.12)使得随方位向频率的表达式与目标所处的方位向位置无关： (2.13)对式(2.12)进行方位向Fourier变换，得到方位向频域信号： (2.14) 两个点目标A，B的距离徙动曲线时域及频域示意图可见将数据变换到方位向频域以后，不同方位向位置的点目标的距离徙动曲线将重合起来。上述过程如上图所示。距离徙动校正处理的实际工作过程一般是针对方位向频域信号，根据式(2.13)由方位向频率计算出的大小，然后对进行相应的距离向移位操作。2.2 RD算法的方位向处理经过距离徙动校正处理的信号可以表示为： (2.15)其中，为距离徙动校正后的参考距离，一般情况下为；为方位向多普勒调频斜率：(2.16)因此是一个在距离向处出现，方位向中心位于，调频斜率为的chirp信号。构造方位向参考函数：(2.17)对进行方位向脉冲压缩处理，处理后的信号为： (2.18)其中，为方位向处理结果的包络，通常情况下也是一个sinc函数。3 SPECAN算法原理SPECAN算法是在RD算法的基础之上发展出来的一种时域和频域混合的SAR成像算法，其距离向处理方法与RD算法相同，而方位向处理方法则与RD算法不同。RD算法的方位向处理采用的是基于相关处理的脉冲压缩算法；而SPECAN则利用方位向信号的线性调频特性，利用去斜（dechirp）处理和频谱分析方法实现方位向的聚焦。3.1 dechirp处理和频谱分析方法对于chirp信号(3.1)构造参考函数(3.2)其中，即的支撑域包含了的支撑域。则与相乘的结果为： (3.3)对进行Fourier变换有：(3.4)式(3.3)过程即称为dechirp处理，式(2.63)的过程称为频谱分析处理。通过观察式(3.2)可见，对于中心时刻位于的chirp信号，经过dechirp及频谱分析处理之后变成了一个在频域中心位于、宽度为的sinc信号。3.2 SPECAN算法的方位向处理SPECAN算法的流程图SPECAN算法正是利用了dechirp及频谱分析的方法来进行数据的方位向处理。其流程如上图所示。SPECAN算法的原理示意如下图所示。运动平台S从S1飞行到S2，飞行距离为一个合成孔径长度，但对应于覆盖地面方位向则为两个合成孔径长度。同一距离上的一族目标的多普勒历程示于图 (b)中，它们是一族chirp信号，每个时刻点接收的信号来自多个目标回波的叠加，而每个频率是由无数个目标的回波信息组成。因为目标的频率是斜线，因此无论在时域还是频域，都不可能把目标分离开。显然，如果选择一个与回波信号频率相反的参考函数（图(c)）对回波信号进行差频处理，就可将回波信号的多普勒历程变为图 (d)所示，即每个目标的频率平行于时间轴，这时就可在频域将信号分开。SPECAN 算法原理示意图下面仍以点目标的回波信号模型来推导SPECAN算法的原理。根据3节的分析，经过与RD算法相同的距离向处理以后的信号为：(3.5)构造相应的参考函数为：(3.6)对回波信号进行dechirp处理，得到： (3.7)对慢时间作方位向Fourier变换，得到最终的成像结果：(3.8)其中，为慢时间频率，方位向处理结果的包络为：(3.9)为点目标回波的持续时间。可见，经过SPECAN算法的距离向和方位向处理后，点目标的处理结果为信号平面上处的一个方位向受天线方向图调制的二维sinc函数，其峰值大小与点目标的后向散射系数有关。4 CS算法原理CS算法利用Chirp Scaling原理，在信号变换到二维频域之前，先初步校正所有距离单元的距离徙动曲线，使之与参考距离处的距离徙动曲线相同。这样的曲线函数仅与方位向有关，并不随距离的变化而变化，因此可以在二维频域通过简单的相位相乘完成距离徙动校正，从而避免了复杂的插值运算，这也正是CS算法与RD算法相比最大的优势所在。CS算法的流程示意图如下图所示。由图中可见，CS算法是以方位向FFT而不是距离向处理开始，并且以方位向IFFT结束，距离向处理则隐含在中间。这种处理流程使得CS算法与RD算法相比，多需要两次数据矩阵转角处理。三次转角处理也是CS算法的一大特点。另外可以看到，在整个处理过程中，CS算法只用到了两种操作：FFT/IFFT和复乘。Chirp Scaling算法流程示意图如图中所示，CS算法共需要进行三次相位因子相乘：第一次相位因子相乘在距离多普勒域进行，目的是进行Chirp Scaling处理，使所有距离单元的距离徙动曲线形状一致，与参考距离处的距离徙动曲线相同；第二次相位因子相乘在二维频域进行，目的是同时完成距离向处理和距离徙动校正，其中距离向处理包括距离压缩和二次距离压缩；第三次相位因子相乘在距离多普勒域进行，目的是补偿Chirp Scaling处理时引入的相位误差，同时完成方位压缩。下面逐一介绍CS算法各个步骤的理论公式及相应的物理意义。4.1 方位向FFT基带回波信号可表示为： (4.1)其中：(4.2)为天线方向图加权，为发射信号包络，为方位向时间，为距离向时间，为发射信号的调频斜率，为目标与雷达的最短斜距。式(2.69)中第一个指数项表示距离向的相位调制，第二个指数项表示方位向的相位调制。根据驻定相位原理，经方位向FFT后，在距离多普勒域的表达式为： (4.3)其中，为复常数，为距离徙动在距离-多普勒域的表示：(4.4)令：(4.5)则有：(4.6)为实际的距离向调频斜率：(4.7)令：(4.8)则有：(4.9)在上面的推导中，式(4.5)定义的称为弯曲因子，由式(4.6)可知，由于不同距离处对弯曲因子的加权不同，因此不同距离处目标的距离徙动曲线也就不同。弯曲因子是用来进行Chirp Scaling处理的关键。式(4.8)定义的称为为距离失真因子，它的存在使得不同目标回波的距离向调频斜率不一致，如果不补偿将导致距离向散焦。对距离失真的补偿就是二次距离压缩处理（SRC）。由此可见，式(4.3)的第一个指数项仍表示距离向相位调制，第二个指数项仍表示方位向相位调制。同时可以看出，目标的距离徙动是随方位向多普勒频率及变化的函数，目标回波的距离向调频斜率也是随以及变化的函数。4.2 Chirp Scaling处理简单说来，Chirp Scaling处理的基本原理是：对目标回波的相位进行微调，使得距离压缩结果在位置上发生偏移，不同距离单元内的目标，其偏移量也不同。借此来调整各距离单元目标的距离徙动曲线，使之与参考距离的距离徙动曲线一致，从而可以对所有目标进行统一的距离徙动校正。构造Chirp Scaling因子： (4.10)其中：(4.12)则对距离处的Chirp Scaling处理可表示为： (4.13)其中，为Chirp Scaling处理引入的相位残差(4.14)为回波信号相位中心的变化轨迹，可表示为：(4.15)也就是说：目标在距离多普勒域的距离徙动曲线变为：(4.16)其变化量为：(4.17)上式说明与距离无关，仅为的函数，其物理意义是：各距离单元内目标的距离徙动曲线形状相同，更加有利于在频域进行距离徙动校正。Chirp Scaling处理中引入的相位残差在方位压缩时较容易补偿，不会损失精度。4.3 距离向处理、距离徙动校正Chirp Scaling处理后，再经距离向FFT，信号在二维频域上的表达式为： (4.18)式(4.18)中，第一个指数项表示方位向相位调制及残差相位；第二个指数项表示距离向chirp调制；第三个指数项则包含了每个点目标的实际距离徙动量。相应的，CS算法的第二个相位因子为： (4.19)式(4.19)中的第一项用来完成距离向处理及二次距离压缩，第二项用来完成距离徙动校正。4.4 相位残差补偿、方位压缩经距离向IFFT后的信号为： (4.20)其中，第一个指数项表示方位向chirp调制，第二个指数项表示进行Chirp Scaling处理时与相乘后引入的相位残差。由此可得，CS算法的第三个相位因子为： (4.21)的第一项用来实现方位向匹配滤波，第二项则用来校正相位残差。4.5 方位向IFFT处理经方位向IFFT处理完成方位压缩后的最终成像结果为。(4.22)其中，及分别为方位向处理及距离向处理后的包络。专心-专注-专业