SAR成像算法docx

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

2、.2 SAR信号模型：SAR信号可以分为距离向信号和方位向信号。首先考虑SAR距离向信号。SAR距离像脉冲可表示为：2s rect( )cos 2 fKr (1.2.1)Tr其中，为脉冲持续时间，Kr为距离向昧冲的调频率，fo为中心频率， 以脉神中心为参考原点。任一照射时刻的反射能量脉冲波形和照射区域内地面反射系数gr的卷积，如下所示：SrgrS (1.2.2)考察距雷达Ro处的一个目标点，其后向散射系数的幅度为A，则式(122)中的gr A2Ro.；c，其中c为光速，2&；c为该点的信号延时。所以可知，该点目标的接收信号为：SrArect()cos 2 f02R/ cKr2R c 2Tr(1

3、.2.3)其中， 表示地表散射过程可能引起的首达信号相位改变。现在考虑方位向信号。由于大多数SAR天线在方位面内没有加权，其单程方向图可以近似为一个sinc函数:Pa0.886 sine(1.2.4)bw其中 为斜距平面内测得的与视线的夹角,bw方位向波束宽度精彩文档0.886,/ La，La为方位向天线长度。由于雷达能量的双程传播过程，接收信号的强度由式(1.2.4)平方给出，并且可以表示成方位时间的函数:WaPa2(1.2.5)其中方位时间与的关系是sin所以，点目标的接收信号可以写成Sr2R /cArect () wac cos0 2R c Kr2R , c 2(1.2.6)其中，RR0

4、2 V2 2, Ro为最短距离，c为波束中心穿越时刻，上述信号其实是一个二维信号，它包含了距离同时间和方位向时间，其中题 离向时间又成为快时间，而方位向时间成为慢时间。由子接收信号Sr 包含了雷达载cos2 fo ，在采样之前，载频必 须通过正交解调过程去除。解调后的单个点目标的基带信号可以表示成复数形式：2R C/.2So , Arect(t )Wa c exp j4 f Rc0 exp j Kr 2R . c(127)其中，系数Ao为一个复常数Ao Aexp j 。1.3分辨率1.3.1距离分辨率SAR的距离分辨率仅由雷达发射被形的频带宽度决定。距离分辨率有斜距分辨率r (沿星载SAR与目

5、标的连线方向量度)和地距分辨率g(沿 地量度)之分。在评价SAR的距离分辨率时，一般用地距分辨率。地距分 辨率和斜距分辨率有如下关系：grr； Sin (1.3.1)斜距分辨率为：cr(1.3.2)2B其中，B为雷达发射波形的频带宽度。1.3.2方位分辨率SAR处理之前的方位向分辨率为波束宽度在地面的投影，即Pbw0.886R(1.3.3)该式成为真实孔径雷达分辨率。而以距离为量纲的合成孔径雷达分辨率可成:LaVg2a(134)其中，,a为处理中加窗引入的展宽因子般的，星载SAR情况下VS ,a 1,方位向分辨率可以写成Pa号。这意味着方位向分辨率是天线长度的一半，与距离、速度和波长等因素无关

6、。这是合成孔径雷达系统最显著的特点。1.4 SAR的距离徙动根据RR02 V2 2，瞬时斜距R 随方位时间 而改变，为的双曲函数。该等式表明目标轨迹（以距离为量纲）是方位时间函数。距离采 样间隔为C 2Fr，其中Fr为距离采样率。这意味着在信号存储器中，照射 时间内的目标轨迹经过不同的距离单元，因此称为“距离单元徙动”或者RCM。1002Q0300400500600500100015002C0025003000图2目标轨迹在不同距离上的变化趋势等式RrdfR02 2R0给出了距离多普勒域中的斜距c2f2D f ,VrJ 4Vr f o2等式，其中IRd f近似f的双曲函数。如图2所示，在方位时

7、域中距离双曲函数的弯曲程度随着距离变量的增加而减小。这是由于，从等式V 22R . R02 Vr2 2 Ro 可以看出，距离变量Ro在分母上，因此,2Ro双曲线随着Ro的增加而逐渐张开。1.5 SAR模糊问题1.5.1距离模糊距离模糊是指前后发射周期的一些回波信号会伴随着所期望的发射周期的回波信号同时被雷达接收，在距离向上产生模糊噪声。如图3所示回波窗口图3 SAR的距离模糊原理图距离模糊现象在机载SAR系统中并不严重.因为此时斜距比较小，观 测带回波的最大延时差相对于脉冲重复周期而言是很小的即使第一模糊 区也是远离观测带的其回波能量也将远小于观测带的回波能量甚至可 能会超出波束的照射范围.而

8、对于星载SAR系统，由于斜距比较大，距离 模糊问题必须考虑。1.5.2方位模糊方位模糊主要是由于较低的脉冲重复频率(PRF)造成的。因为目标回波谱是以脉冲重复频率(PRF)为周期重复出现的，在主谱之外的回波信号将折叠到主谱区，如图4所示图4 SAR的方位模糊原理图距离向模糊和方位向模糊取决于脉冲重复频率(PRF)的选择和测绘带的位置。较低的PRF会使方位向模糊增加；较高的 PRF会增加距离向模 糊，或者使测给带宽度受限。故距离模糊和方位模糊是一对相互矛盾的量， 而PRF的选择要综合多种因素折中考虑。2 SAR成像算法2.1 RDA ( Range-Doppler Algorithm )R-D算

9、法基于匹配滤波的原理，将SAR成像中的二维联合处理简化为两个一维的级联R-D算法的参考函数选择为接收信号频谱的复共轭，时域上是接收信号的逆时复共轭。R-D算法的实现步骤为，先对每个回波信号进行距离向压缩，然后在R-D域中对距离向压缩后的数据进行距离徙动校正，在大斜视角情况下再进行二次距离压缩，最后进行方位向压缩2.1.1距离向压缩根据SAR的成像原理，得到SAR回波信号2j4 r s / j K 2r s /cg s,Wa s a 2r s /c ee(2.1.1)式中，是距离向的快时间变量，s是方位向的慢时间变量，a t是矩形窗信号，r s是卫星与地面目标的距离，K是线性调频脉冲的调频斜率，

10、 是点目标的后向散射系数，Wa s是雷达天线增益由(2.1.1)式可以得到回波信号的距离向频谱G s, fWa s e j4 r s /1/2 j sgn K /4 j 4 f r s /c j f 2/Kee j e j (2.1.2)式中，f是距离向频率,p是发射脉冲的宽度，Br K p是距离向带宽,1/2 j sgn K /4 j fee/K(2.1.3)Br/2距离向压缩就是对(2.1.2)式进行匹配滤波，滤波函数为其时域形式为(2.1.4)滤波后的频谱为R fsRcD2fRfD4fR s2(2.1.5)滤波后的时域信号为has,BrWasine Br 2r s /c (2.1.6)可

11、见，距离压缩后的信号仍然是距离向和方位向的二维信号，距离向 和方位向的耦合仍然没有解除。2.1.2距离徙动从(2.1.6)式可以看出，经距离压缩后不同的点目标响应出现在不同的距离向上，这是由距离徙动造成的。根据SAR的多普勒历程，有2r sRcfD s Sc /2fR s sc /4 (2.1.7)式中，Rc表示位于波束照射中心的目标与雷达之间的距离，Sc为照射到目标的时刻，fD为多普勒中心频率，fR为多普勒调频斜率，R r s Rc是距离徙动。由于时间一带宽积较大，依据驻定相位定理,驻定相位点？与多普勒频率之间的锁定关系为? s scffD / fR(2.1.8)于是时域距离徙动曲线在频域内

12、表示成R fsD2fRfsfD2fsfD(2.1.9)4fR式中，fs是方位向频率。为了消除距离徙动引起的距离向和方位向耦 合，必须做距离徙动校正。同一合成孔径内，不同方位但同一距离的目标 点的距离弯曲具有相同的频域形式，因此弯曲校正一般在频域内进行。对 于距离走动校正，理想情况在时域校正地球自转引起的距离走动，在频域 校正斜视产生的距离走动。当距离走动不大时，时域校正简单有效但是 很小的斜视角会产生大的距离走动，超越聚焦深度，导致时域校正失效。因此，为了提高运算效率，扩展算法的适用范围，我们在频域内进行距离 徙动校正。2.1.3距离徙动校正设雷达波束中心照射到目标的时刻氏为时间零点，结合(2

13、.1.9)式和(2.1.5)式，得到二维频谱aWae j4 RJ e j4f Rc/cCPl1/2 jesgn p /4e2 2 fcfs/fRfcfej2fDfs/fRe jfDDDf / fcfR(2.1.10)式中，P1 fR 1 f / fc。因为f Br/2，雷达载波频率fc通常比带 宽Br大得多，于是可以认为P1 fR。对ejfcfs/fRfcf中的fc / fc f在f 0处做Taylor级数展开并做二阶近似得2 2fc / fc ffc / fc f / fcf / fc (2.1.11)经过这两个近似处理后的距离压缩信号二维频谱为B fs,fHa fs, f T2 fs, f

14、 r ?Wafse j4 Rc/,1/2 j sgn fR /4 j 2 fR eefD ?2 2 (2.1.12)fD fs/ f Rjfs fD /fRj2 f 2 R3 r ?ee(2.1.12)式中，f的二次项就是距离徙动引起的距离向上的线性调频为了消除线性调频，需要进行二次距离压缩(SRC)理想的二次距离压缩采j f 2 f 2/ f 2 f用频域参考函数进行滤波T2 fs, f e s R，去掉(2.1.12)式中的最后 指数项，但实现起来比较困难为简化处理，一种方法是取多普勒中心频率 fD作为方位向的参考频率，取时域参考函数T2e j fc2 2/fD的频域函数j f2fD/ f

15、c2fRT2 f e进行二次距离压缩，消除距离向的扩展B fs，f Hawa fs ej4 Rc/1/2 j sgn fR /4 j2 fpfs/fRjfs fD 2/ fR j2 f 2 R r s（2.1.13）e ee式中，? fs fD / fR为驻定相位点，r ?fD? 土 s2就是距2 2离徙动部分也可采用其他不同的简化方法，确定不同的参考函数对经过SRC后的二维频谱B fs, f做距离向的傅立叶反变换，得B fs,BrWa fs ej41/2 jesgn fR/4e jfD /fRsinc Br(2.1.14)接着使用插值运算和数据阵移动，使得si ncBr2 Rc r /csi

16、 nc Br2Rc / c (2.1.15)进行频域距离徙动校正B fs, ,?02 (2.1.16)Dfj4 RJ 1/2 j sgn fR /4 j fs fD /fR .Dco /BrWa fs e c fR eesin cBr2Rc /c式中， s 2r ? / c2.1.4方位向压缩经过距离徙动矫正后的 B fs, ,?是一个方位向的频谱函数，用0j2 fs fDS2,2方位向时域参考函数Ts s e对应个频域参考函数Ts fsej fs f，fD进行匹配滤波，完成方位压缩。B fs，BrWafs e j4 Rc/?Ts0(2.1.17)- 1/2 j sgn fR /4 ._ _.

17、 .ejsincBr2Rc /c对上式进行傅里叶反变换，得到时域目标点的散射系数?Qcsinc Br2Rc/c csinc Bas1/2 (2.1.18)j4 R /1/2；f !A/Q BaBrWa s e j41/2 ej sgn fR /4即在方位向上的s时刻和距离向延迟2Rc.fc后时刻所在对应的目标点的像素值这样就实现了对目标点的成像。RD算法流程图如图5所示图5 RD算法流程图2.2 CSA (Chirp -Sealing Algorithm)在SAR成像中，距离-多普勒(R-D)算法是为民用SAR开发的第个成像处理算法。由于它等兼顾成熟、简单、高效和精确等因素，至今 仍是使用最广

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

19、和调频率发生微小的变化。因此，采用新的调频率进行距离压缩 后,信号的位置产生位移，修正了不同距离上目标距离徙动曲线的微小差 别，使得所有距离徙动曲线具有相同的形状，这样就可以进行统一的距离 徙动校正和，实现精确成像。算法流程图如下所示。距离多普勒域距离多 普勒域二维频域图6 CS算法流程图2.2.1 Chirp Scali ng 原理Chirp Scali ng充分利用了发射信号为线性调频信号的性质。对于一个调频斜率为Km、中心为1的线性调频信号，若把它乘上另一调频斜率为 KmCs、中心为2的线性调频信号，贝U有以下结果,exp j KmexpKmCsexp j Km 1 Cs2 .newJ2

20、2(2.2.1)我们得到一个新的中心new1 Cs 21 Cs和一个残余相位2K Cs 12K m1 Cs。通过选取适当的参数，实现对徙动的校正。222方位向FFT首先将回波信号在方位向进行 FFT,使信号数据转换到距离一多普勒域,在距离压缩之前，CS方法利用了距离向线性调频信号的特性。根据驻留相位 原理，对回波信号进行方位向FFT,信号在R-D域表示为t, ,rCG -? m2v22Rf ft,rcexp j Km ft；r2Rf ft,rcexp jJ(2.2.2)式中,12vRf ft；rr 1 Cs ftr rCs ftKm为包含SRC参数在内的综合调频斜率,C为弯曲因子，表示信号轨迹

21、与多普勒频率之间的关系。Rf ft;r为距离徙动在距离一多普勒域的表示。Km ft；r1 _1 厂1 Ksm_ft/ 1K2(2.2.3)Kft cKsrc(2.2.4)tR ft；r r Cs (2.2.5)它显然与距离r有关。我们对S .做如下的Chirp Scali ng 处理S2 ft, ,rSi ft, ,r Hi ft, ,r（2.2.6）2Hi ft, ,r exp j Km 论牯 Csref（227）2式中ref 2rref 1Cs ft为一固定的参考距离，一般可取在成像测绘C带的中心；ref为rref所对应的回波中心。1(2.2.9)Csft, ; r CG m exoj K

22、m ft；ref 1 Cs- rc2rref Cs.4 r _exp j (2.2.8)exp jftCS2ref如果对S2做距离向压缩，则其峰值所对应的徙动为R ft；rrref Cs（2.2.io）该徙动已与距离无关，而与参考距离上的一致。2.2.4距离向FFT将Chirp Scali ng相位修正后的信号做距离向FFT，回波信号转换为二维频域信号：rft,fr CG2丄m2v2expm 1 CsgeXP j Km 1 Cf f. 4 rexp j-Km 1 Cs j rrref Csfc(2211)图 7 Chirp Scali ng示意ft2 1ref K mc2r rref (2.2.12)c式中，f为距离向频谱偏移。考虑式中的相位项，第一项对应横向聚焦第二项包含f2，对应距离压缩和SRC第三项是f的线性项，对应目标的正确位置和徙动rrefCs。2.2.5距离徙动校正、距离压缩及二次距离压缩S2通过乘以一个H2可以完成距离徙动校正及距离聚焦处理 （包括多普勒域内