各种SAR成像算法总结.doc
33页
各种SAR成像算法总结1SAR成像原理SAR成像处理的目的是要得到目标区域散射系数的二维分布，它是一个二维相关处理过程，通常可以分成距离向处理和方位向处理两个部分在处理过程中，各算法的区别在于如何定义雷达与目标的距离模型以及如何解决距离－方位耦合问题，这些问题直接导致了各种算法在成像质量和运算量方面的差异一般来说，忽略多普勒频移所引起的距离向相位变化，距离向处理变为一维的移不变过程且相关核已知，即退化为一般的脉冲压缩处理；同时将雷达与目标的距离按2阶Taylor展开并忽略高次项，则方位向处理也是一个一维的移不变过程，并退化为一般的脉冲压缩处理，这就是经典的距离多普勒(Range-DopplerRD)算法的实质若考虑多普勒频移对距离向相位的影响，同时精确的建立雷达与目标的距离模型，则不论距离向处理还是方位向处理都变为二维的移变相关过程线性调频尺度变换(Chirp-ScalingCS)算法即在此基础之上将二维数据变换到频域，利用ChirpScaling原理及频域的相位校正方法，对二维数据进行距离徙动校正处理、距离向及方位向的聚焦处理，最终完成二维成像处理当方位向数据积累延迟小于全孔径时间(即方位向为子孔径数据)的情况下，方位向处理必须使用去斜(dechirp)处理及频谱分析的方法。
在RD和CS算法的基础之上，采用dechirp处理及频谱分析的方法完成方位向处理的算法分别称为频谱分析(SPECAN)算法和扩展CS(ExtendedChirpScalingECS)算法1.1SAR成像原理本节以基本的正侧视条带工作模式为例，对SAR的成像原理进行分析和讨论正侧视条带SAR的空间几何关系如下图所示图中，ao卩平面为地平面，oy垂直于ao卩平面SAR运动平台位于S点，其在地面的投影为G点SAR运动平台的运动方向Sx平行于卩，速度大小为vOSAR天线波束中心与地面的交点a为C，CG与运动方向Sx垂直；S与C的距离为R，€BSB称为天线波束的方s12位向宽度，大小为PP为测绘带内的某一点，一般情况下取斜距平面CSP进a行分析，称SAR运动的方向Sx为方位向(或方位维)，称天线波束指向方向SC为距离向(或距离维)xvR⑴G测绘带BvCRsa正侧视条带SAR几何关系示意图假定P的方位向坐标为X；在t时刻，SAR运动平台S与P的距离为R(t)若当t„0时亥【」，SAR运动平台位于方位向0点，则当t时亥【」，R(t)的表达式为：R(t)„(R2+(vt-X)2(1.1)Wsa将式(1.1)在t„X/V附近进行2阶Taylor展开，有:a
c发射信号s(t)经点目标P散射后，雷达接收机收到的信号s(t)为：s(t)…Re^KAej€c(t„a)F(x),(1.4)r其中：c为光速，K为复常数，a为回波信号相对于发射信号的时间延迟：a…2R(t)/c(1.5)F(x)为考虑雷达水平方向增益变化而引入的加权函数若不考虑雷达天线的加权作用，即令F(x)…1，则式(1.4)变为:s(t)…Reej€c(t„a)(1.6)根据式(1.6)，雷达运动平台相对于点目标的运动将造成回波信号的相位随时间不断变化，从而引起回波瞬时频率的变化，产生多普勒频移多普勒频移量f(t)为:f(t)…一1d(€a)dt„1dI2€R(t)2兀dtIc(1.7)(t)1d2兀dt1dJ2€R2兀dt 也就是说，由于雷达运动平台匀速直线前进，回波信号s(t)在r方位向将为线性调频(chirp)信号：(1.#)s(t)二Re a由式(1.14)、(1.15)得到B的表达式为:aB二1fIT二丝T二Mla2€d」s,Rs,(1.16)考虑到对于方位向天线直径为D的天线，近似有：a(1.17)因此，SAR的方位向理论分辨率™为:a(1.#)€vvD€———a—aaB2vP/,2aaa从上述分析可以看出，由于雷达运动平台作等高匀速直线运动，使得目标的回波信号在方位向上具有线性调频特性，对回波信号进行脉冲压缩处理，可以获得方位向的高分辨率在理想情况下，SAR方位向分辨率与雷达平台的速度、飞行高度、作用距离、雷达工作波长等参数无关，只与天线尺寸有关，为天线方位向口径尺寸的一半，这是SAR的一大特点和优势1.2SAR回波信号模型1.1节分析了SAR成像的基本原理，本节推导SAR回波信号的数学模型,给出SAR信号处理的理论基础chirp信号是SAR系统中最常用的发射信号形式假设雷达发射的chirp脉冲串s(t)为：t(1.19)s(t)-艺p(t-nPRT)tn—-„其中，p(t)—s(t)cos00s(t)为发射信号的包络，k为chirp信号的调频斜率，t为发射信号脉宽，0rsf为发射信号的中心频率，PRT为脉冲重复周期。 则雷达于时刻t，接收到斜距0为R(t)处目标反射的回波信号s(t)为：tR(t)+2R(t)s(t)—a•Wt-st-ra•c_t•c_(1.20)其中，a为目标的后向散射特性,W()为方位向的天线方向性函数，c为a艺a・WtR(t)t—pt—nPRT—2R(t)an—-„•c_•c_(1.#)光速1.#)(1.#)•s(t-nPRT—2R(t)/c)•0s(t)经正交解调后的复信号s(t)可以表示为:rs(t)„艺,・Wan„—g(1.21)exp“—j+R(t)'exp“jpt—nPRT—其中，尢为雷达工作波长式(1.21)中的两个指数项分别代表方位向的相位调制和距离向发射的相位调制考虑到相对于雷达发射脉冲而言，W(t)和R(t)是时间t的慢变化函数，可以a作如下近似:W(t)„W(nPRT)aa(1.22)R(t)„R(nPRT)(1.23)同时，将时间t分解为快时间分量—和慢时间分量t之和，即:at„—+1，t„k•PRTaa(1.24)通过变量置换，可以将s(t)转换成二维形式:s(—,t)„Q・W(t)exp“—jjR(t)”™aaI九a-t5(t)5—a2R(t)”™{s(—)exp[jp(—)]}—0(1.25)„,•W(t)exp‘—jJR(t)'5——aa•九(5(t)s(—)exp[jp(—)]}a02R(t)a—c其中，™表示对t的卷积，™表示对—的卷积，™表示二维卷积。 ta—a因此，雷达系统接收回波信号的过程，可以看作是地面目标的后向散射特性通过一个线性系统的过程式(1.25)可简化表示为:s(—,t)„，™h(—,t)aa(1.26)其中，h(—,t)为线性系统的冲激响应函数:ah(T,t)…W(t)expaaa4€-j亍R(t),aT2R(t)a—c'sG)exp[j‘(T)]T0(1.27)(1.28)(1.29)(1.30)(1.31)式(1.27)可以进一步表示为：h(T,t)…h(t,t)'h(t)a1aT2其中，h(T,t)…W(t)exp•-jR(t)|sT-2R(ta)1aaa>,ac丿h(t)…sG)exp{j‘(t)}20则式(1.26)可进一步表示为：s(T,t)…◎'h(T,t)'h(T)a1aT2式(1.29)中，h(t,t)的指数项代表了由于雷达运动平台与目标间相对运动所1a带来的方位向相位调制如果对R(t)采用式(1.2)所示的2阶Taylor展开方式，a则回波的方位向相位为慢时间t的2次函数，即一个chirp信号；h(t,t)的冲击a1a函数表达式代表了由于相对运动，回波包络的中心在距离向上的位置发生变化，即距离徙动现象。 2RD算法原理RD算法流程如下图所示，包括距离压缩处理、方位压缩处理两个主要处理步骤，以及作为辅助处理步骤的距离徙动校正处理由于具有概念简单、易于实现、处理效率高等优点，RD算法成为最经典、最成熟的SAR成像处理算法距离徙动校正处理方位向FFT方位向IFFT方位向参考函数成像结果RD算法流程RD算法的本质是对R(t)采用式(2.2)所示的2阶Taylor展开方式，将距离向a处理和方位向处理解耦，分解为两个一维处理分别完成。
