SAR成像算法docx

1、SAR成像1 合成孔径雷达（SAR）1.1 SAR简介合成孔径雷达(SAR)是一种可以全天候、全天时工作的高分辨率成像雷达。它利用天线和目标之间相对运动而形成等效合成孔径，解决了雷达设计中高分辨率与大尺寸天线和短工作波长之间的矛盾，在遥感和国防中潜在着极大的应用价值。星载SAR一般工作在正侧视状态,但在特殊应用中，也会工作在斜视状态。图1给出了星载SAR正侧视模式的空间几何关系。飞行路径在地面上的投影(地面轨迹)方向称为方位方向，而与其垂直的方向称为距离方向。距离向使用脉冲压缩技术实现高分辨率；方位向利用多普勒效应,经过相干处理得到高分辨率。图1 SAR的几何关系1.2 SAR信号模型：SAR

2、信号可以分为距离向信号和方位向信号。首先考虑SAR距离向信号。SAR距离像脉冲可表示为: (1.2.1)其中，为脉冲持续时间，为距离向昧冲的调频率，为中心频率， 以脉神中心为参考原点。任一照射时刻的反射能量脉冲波形和照射区域内地面反射系数的卷积，如下所示: (1.2.2) 考察距雷达处的一个目标点，其后向散射系数的幅度为，则式(1.2.2)中的，其中为光速，为该点的信号延时。所以可知，该点目标的接收信号为:(1.2.3)其中，表示地表散射过程可能引起的首达信号相位改变。现在考虑方位向信号。由于大多数SAR天线在方位面内没有加权，其单程方向图可以近似为一个函数: (1.2.4)其中为斜距平面内测

3、得的与视线的夹角，方位向波束宽度，为方位向天线长度。由于雷达能量的双程传播过程，接收信号的强度由式(1.2.4)平方给出，并且可以表示成方位时间的函数: (1.2.5)其中方位时间与的关系是。所以，点目标的接收信号可以写成: (1.2.6)其中，, 为最短距离，为波束中心穿越时刻，上述信号其实是一个二维信号，它包含了距离同时间和方位向时间，其中题离向时间又成为快时间，而方位向时间成为慢时间。由子接收信号包含了雷达载，在采样之前，载频必须通过正交解调过程去除。解调后的单个点目标的基带信号可以表示成复数形式: (1.2.7)其中，系数为一个复常数。1.3 分辨率1.3.1距离分辨率SAR的距离分辨

4、率仅由雷达发射被形的频带宽度决定。距离分辨率有斜距分辨率（沿星载SAR与目标的连线方向量度）和地距分辨率（沿地量度）之分。在评价SAR的距离分辨率时，一般用地距分辨率。地距分辨率和斜距分辨率有如下关系： (1.3.1)斜距分辨率为： (1.3.2)其中，为雷达发射波形的频带宽度。1.3.2方位分辨率SAR处理之前的方位向分辨率为波束宽度在地面的投影，即 (1.3.3)该式成为真实孔径雷达分辨率。而以距离为量纲的合成孔径雷达分辨率可成: (1.3.4)其中，为处理中加窗引入的展宽因子。一般的，星载SAR情况下，方位向分辨率可以写成。这意味着方位向分辨率是天线长度的一半，与距离、速度和波长等因素无

5、关。这是合成孔径雷达系统最显著的特点。1.4 SAR的距离徙动根据，瞬时斜距随方位时间而改变，为的双曲函数。该等式表明目标轨迹（以距离为量纲）是方位时间函数。距离采样间隔为，其中为距离采样率。这意味着在信号存储器中，照射时间内的目标轨迹经过不同的距离单元，因此称为“距离单元徙动”或者RCM。图2 目标轨迹在不同距离上的变化趋势等式给出了距离多普勒域中的斜距等式，其中近似的双曲函数。如图2所示，在方位时域中距离双曲函数的弯曲程度随着距离变量的增加而减小。这是由于，从等式可以看出，距离变量在分母上，因此，双曲线随着的增加而逐渐张开。1.5 SAR模糊问题1.5.1距离模糊距离模糊是指前后发射周期的

6、一些回波信号会伴随着所期望的发射周期的回波信号同时被雷达接收，在距离向上产生模糊噪声。如图3所示。图3 SAR的距离模糊原理图距离模糊现象在机载SAR系统中并不严重因为此时斜距比较小，观测带回波的最大延时差相对于脉冲重复周期而言是很小的即使第一模糊区也是远离观测带的其回波能量也将远小于观测带的回波能量甚至可能会超出波束的照射范围而对于星载SAR系统，由于斜距比较大，距离模糊问题必须考虑。1.5.2方位模糊方位模糊主要是由于较低的脉冲重复频率（PRF）造成的。因为目标回波谱是以脉冲重复频率（PRF）为周期重复出现的，在主谱之外的回波信号将折叠到主谱区，如图4所示。图4 SAR的方位模糊原理图距离

7、向模糊和方位向模糊取决于脉冲重复频率(PRF)的选择和测绘带的位置。较低的PRF会使方位向模糊增加；较高的PRF会增加距离向模糊，或者使测给带宽度受限。故距离模糊和方位模糊是一对相互矛盾的量，而PRF的选择要综合多种因素折中考虑。2 SAR成像算法2.1 RDA（Range-Doppler Algorithm）R-D算法基于匹配滤波的原理，将SAR成像中的二维联合处理简化为两个一维的级联R-D算法的参考函数选择为接收信号频谱的复共轭，时域上是接收信号的逆时复共轭。R-D算法的实现步骤为，先对每个回波信号进行距离向压缩，然后在R-D域中对距离向压缩后的数据进行距离徙动校正，在大斜视角情况下再进行

8、二次距离压缩，最后进行方位向压缩.2.1.1距离向压缩根据SAR的成像原理，得到SAR回波信号 (2.1.1)式中，是距离向的快时间变量，是方位向的慢时间变量，是矩形窗信号，是卫星与地面目标的距离，是线性调频脉冲的调频斜率，是点目标的后向散射系数，是雷达天线增益由(2.1.1)式可以得到回波信号的距离向频谱 (2.1.2)式中，是距离向频率，是发射脉冲的宽度，是距离向带宽，距离向压缩就是对(2.1.2)式进行匹配滤波，滤波函数为 (2.1.3)其时域形式为 (2.1.4)滤波后的频谱为 (2.1.5)滤波后的时域信号为 (2.1.6)可见，距离压缩后的信号仍然是距离向和方位向的二维信号，距离向

9、和方位向的耦合仍然没有解除。2.1.2距离徙动从(2.1.6)式可以看出，经距离压缩后不同的点目标响应出现在不同的距离向上，这是由距离徙动造成的。根据SAR的多普勒历程，有 (2.1.7)式中，表示位于波束照射中心的目标与雷达之间的距离，为照射到目标的时刻，为多普勒中心频率，为多普勒调频斜率，是距离徙动。由于时间一带宽积较大，依据驻定相位定理，驻定相位点与多普勒频率之间的锁定关系为 (2.1.8)于是时域距离徙动曲线在频域内表示成 (2.1.9)式中，是方位向频率。为了消除距离徙动引起的距离向和方位向耦合，必须做距离徙动校正。同一合成孔径内，不同方位但同一距离的目标点的距离弯曲具有相同的频域形

10、式，因此弯曲校正一般在频域内进行。对于距离走动校正，理想情况在时域校正地球自转引起的距离走动，在频域校正斜视产生的距离走动。当距离走动不大时，时域校正简单有效但是很小的斜视角会产生大的距离走动，超越聚焦深度，导致时域校正失效。因此，为了提高运算效率，扩展算法的适用范围，我们在频域内进行距离徙动校正。2.1.3距离徙动校正设雷达波束中心照射到目标的时刻为时间零点，结合(2.1.9)式和(2.1.5)式，得到二维频谱 (2.1.10)式中，。因为，雷达载波频率通常比带宽大得多，于是可以认为。对中的在处做Taylor级数展开并做二阶近似得 (2.1.11)经过这两个近似处理后的距离压缩信号二维频谱为

11、 (2.1.12)(2.1.12)式中，的二次项就是距离徙动引起的距离向上的线性调频为了消除线性调频，需要进行二次距离压缩(SRC)理想的二次距离压缩采用频域参考函数进行滤波，去掉(2.1.12)式中的最后指数项，但实现起来比较困难为简化处理，一种方法是取多普勒中心频率作为方位向的参考频率，取时域参考函数的频域函数进行二次距离压缩，消除距离向的扩展 (2.1.13)式中，为驻定相位点，就是距离徙动部分也可采用其他不同的简化方法，确定不同的参考函数对经过SRC后的二维频谱做距离向的傅立叶反变换，得 (2.1.14)接着使用插值运算和数据阵移动，使得 (2.1.15)进行频域距离徙动校正 (2.1

12、.16)式中，2.1.4 方位向压缩经过距离徙动矫正后的是一个方位向的频谱函数，用方位向时域参考函数对应个频域参考函数进行匹配滤波，完成方位压缩。 (2.1.17)对上式进行傅里叶反变换，得到时域目标点的散射系数 (2.1.18)即在方位向上的时刻和距离向延迟后时刻所在对应的目标点的像素值这样就实现了对目标点的成像。RD 算法流程图如图5所示。图5 RD算法流程图2.2 CSA（Chirp -Scaling Algorithm）在SAR成像中，距离-多普勒(R-D)算法是为民用SAR开发的第一个成像处理算法。由于它等兼顾成熟、简单、高效和精确等因素，至今仍是使用最广的一种算法。但在一定条件下，

13、该算法存在两点不足：首先，当用较长的核函数提高距离徙动矫正(RCMC)精度时，运算量较大。其次，二次距离压缩(SRC)对方位频率的依赖性问题较难解决，从而限制了其对某些大斜视角和长孔径SAR的处理精度。其基本思路是,在信号变换到二维波数域之前校正所用距离单元的距离徙动曲线,使之与参考距离相同,这样就可以在二维波数域通过简单的相位相乘完成距离徙动校正,从而避免了插值运算。在CSA中,Chirp Scaling原理是指线性调频信号与一个具有适当相关调频率的调频信号(称为Chirp Scaling因子)相乘,结果仍然是一个调频信号,只是相位中心和调频率发生微小的变化。因此,采用新的调频率进行距离压缩

14、后,信号的位置产生位移,修正了不同距离上目标距离徙动曲线的微小差别,使得所有距离徙动曲线具有相同的形状,这样就可以进行统一的距离徙动校正和,实现精确成像。算法流程图如下所示。图6 CS算法流程图2.2.1 Chirp Scaling原理 Chirp Scaling充分利用了发射信号为线性调频信号的性质。对于一个调频斜率为、中心为的线性调频信号，若把它乘上另一调频斜率为、中心为的线性调频信号，则有以下结果， (2.2.1)我们得到一个新的中心，和一个残余相位。通过选取适当的参数，实现对徙动的校正。2.2.2 方位向FFT首先将回波信号在方位向进行FFT，使信号数据转换到距离一多普勒域，在距离压缩

15、之前，CS方法利用了距离向线性调频信号的特性。根据驻留相位原理，对回波信号进行方位向FFT，信号在R-D域表示为 (2.2.2)式中，、为包含SRC参数在内的综合调频斜率，C为弯曲因子，表示信号轨迹与多普勒频率之间的关系。为距离徙动在距离一多普勒域的表示。 (2.2.3) (2.2.4)2.2.3 Chirp Scaling相位相乘我们需要校正的徙动为 (2.2.5)它显然与距离有关。我们对做如下的Chirp Scaling处理 (2.2.6) (2.2.7)式中为一固定的参考距离，一般可取在成像测绘带的中心;为所对应的回波中心。 (2.2.8) (2.2.9)如果对做距离向压缩，则其峰值所对

16、应的徙动为 (2.2.10)该徙动已与距离无关，而与参考距离上的一致。2.2.4 距离向FFT将Chirp Scaling相位修正后的信号做距离向FFT，回波信号转换为二维频域信号： (2.2.11) 图7 Chirp Scaling示意 (2.2.12)式中，为距离向频谱偏移。考虑式中的相位项，第一项对应横向聚焦;第二项包含，对应距离压缩和SRC;第三项是的线性项，对应目标的正确位置和徙动。2.2.5 距离徙动校正、距离压缩及二次距离压缩通过乘以一个可以完成距离徙动校正及距离聚焦处理(包括多普勒域内的二次距离压缩及补偿线性调频率在Chirp Scaling中的变化)。(2.2.13) (2.2.14) 式(2.2.13)中的第一项完成了二次距离压缩及距离向聚焦处理，第二项完成了距离徙动的校正。2.2.6 距离向IFFT所有距离向调制的相位补偿后，在距离向做IFFT，得到。2.2.7 方位滤波及残余相位消除进行横向聚焦处理，在上乘上下面的横向聚焦参考函数 (2.2.15)得到 (2.2.16)2.2.8 方位向IFFTCS算法的最后一步是对信号做方位向IFFT (2.2.17)完成了成像处理。式中为一复常数，和分别为距离和横向上的目标冲击响应，前者与发射