合成孔径雷达成像

《合成孔径雷达成像——算法与实现》读书笔记(1—4章)徐一凡第1章概论1.1SAR在遥感邻域运用越来越多的原因雷达自带照射源，在黑夜中同样能出色工作一般雷达所使用电磁波几乎可以无失真地穿透水汽云层物质的光学散射能量与雷达电磁波散射能量不同，二者可以进行互补1.2合成孔径雷达的分类合成孔径雷达分为：条带式合成孔径雷达、扫描式合成孔径雷达、聚束合成孔径雷达、逆合成孔径雷达、双站合成孔径雷达和干涉合成孔径雷达1.3距离徙动合成孔径雷达对大量的回波脉冲进行处理，由于合成孔径内传感器的移动，雷达与目标的距离时间变化，这个变化是引起回波数据多普勒频移的原因，然而这种距离变化同时也导致了存储数据的距离徙动现象，如图1所示：图1距离徙动第2章信号处理基础2.1线性卷积2.1.1连续时间卷积一个信号s(t)通过滤波器h(t)的过程可以理解为卷积运算，输出y(t)可表示为:J(t)=s(t)⑤h(t)=js(u)h(t—u)du=jJ(t)=s(t)⑤h(t)=其理解如下：一个冲击信号导致一个冲击响应，无数个冲击信号导致无数个冲击响应，由于符合线性运算，无数个冲击信号叠加为s(t)，其冲击响应为无数个冲击响应的叠加，即J(t)。连续时间卷积几何计算过程：反褶、移位、重合部分乘积积分。连续时间卷积满足线性时不变特性、交换特性。相关定乂如下:军①，(t)=js(u)h*(u-1)du-s相关不反褶，h(t)取复共轭，相关不可交换。气(t)=七(-t)二维卷积y(t,t)=s(t,t)⑤h(t,t)TOC\o"1-5"\h\z121212=ffs(u,u)h(t-u,t-u)dudu12112212-s-s=ffs(t-u,t-u)h(u,u)dudu11221212-s-s由于SAR信号是二维信号，故滤波形式为二维，但通常可将二维滤波器解耦为两个一维滤波器。h(t,t)=h(t)⑤h(t)1211222.1.2离散时间卷积对于离散信号s(n)和滤波器h(n)，对它们进行离散时间卷积，使用的MATLAB函数为conv(s,h)，其运算过程如图2所示：信号MXXXXXXXX滤波器岚o0O>OO0输出冲)♦。。0。◊00。。对应的时间门0123456789图2离散时间卷积过程从图2中可以看出，对于信号长度为K，滤波器长度为M的离散卷积，其输出y(n)的长度为K+M-1，但是在n=0和1及n=8和9上，信号仅与滤波器系数的子集相乘，输出结果为部分卷积的。二维卷积二维卷积滤波器h(n,n)可表示为h(n,n)=h(n)⑤h(n)，例如12121122「1221一h(n1,n2)=36631221

该式可分解为4（〃i）二和耿2）』221。在MATLAB该式可分解为4（〃i）二上述分解可以显著减少算术运算，同时还节省了存储空间。2.2傅里叶变换离散傅里叶变换中，时域第一个点g（0）对应着零时刻，信号在时域上以1/f等间隔采s样。类似地，频域第一个点G（0）对应的频率为零频，频域采样间隔为N/f，频率样本G（k）s对应的频率为kf/N。s频谱泄露来自非周期信号的周期假设。傅里叶变换的性质：复共轭、线性、尺度变化特性、位移/调制、均值、对称性、Parseval定理、卷积/乘法。补零：某一域（频域或时域）中的序列补零相当于对另一域进行升采样，这使得另一域中的数据量增大，但不会改变序列的信息内容（例如带宽）。二维扭曲和旋转：在某一域中沿一个数据轴进行的扭曲，将导致另一域中沿另一轴的数据扭曲。—M务也）+~~♦。（/1,为+口/1）如一。和）*—某一域中的角度旋转，将导致另一域中同样的角度旋转。时间带宽积（TBP）：信号3dB时宽与3dB带宽的乘积。矩形函数和sinc函数的TBP近似为1.峰值旁瓣比：最大旁瓣功率与峰值功率的比值8528（*-nTJ4一►££心-"）

n=—c»rt=—ao2.3卷积的DFT计算由于DFT和IDFT的循环性质，卷积实际上是周期卷积或循环卷积。通过将序列补零至合适的长度，就可以利用DFT进行线性卷积。注意弃置点、有效点还有部分卷积点。2.4信号采样当对信号进行N点DFT时，输出序列第一个样本的频率为零，最后一个样本的频率为(N-1)f/Nnf。DFT只能“观察”到频谱的基本部分(即基带频谱)。通过观察DFT输出不可能判断出初始信号所处的“实际频段”。2.5平滑窗在脉冲压缩中，使用窗来控制旁瓣，同时尽可能保持高的分辨率。由于旁瓣的降低会导致分辨率的展宽，故需折中考虑。通常在处理中使用Kaiser窗。2.6插值非整数点X上的信号g3)可以通过DFT的平移/调制性质或者频谱补零得到。但是这些方法不够灵活，所以使用sinc插值的方法。g⑴=Zgd(i)h(x-i)ih(x)=sinc(x)注意插值核的归一化。非基带信号的插值方法：1、将信号转到基带；2、将基带滤波器移至信号中心频率处。2.7点目标分析IRW：冲击响应宽度指冲击响应的3dB主瓣宽度，在SAR处理中又称为图像分辨率。PSLR：峰值旁瓣比指最大旁瓣与主瓣的高度比，sinc函数为-13dB，比较大。一维ISLR：ISLR=10lOg{P°ta「Lin}10Pmain峰值位置：指二维样本空间中的冲击响应峰值位置。信号峰值幅度：指峰值点处目标响应幅度。相位：某一方向上测量出的点目标相位特性可能与另一方向不同。第三章线性调频信号的脉冲压缩脉冲压缩是一种频谱扩展方法，用于最小化峰值功率、最大化信噪比以及获得高分辨率目标。3.2线性调频信号3.2.1时域表达s（t）=rect（彳）exp{jKt2}相位是时间的二次函数。"土警=K、砒*|T，时间带宽积TBP*I口3.2.2线性调频脉冲的频谱f、r.f2、G（f）=rect（吉）exp{-"J}包络与矩形包络近似一致，相位为二次相位。KTK图3,2线性调频脉冲的复频谱3.2.3调频信号采样f过采样因子a=，过采样因子应选在1.1-1.4之间。OSIKIT3.2.4频率和时间不连续性补零区的选择会影响频谱相位。3.3脉冲压缩3.3.1脉冲压缩原理如果发射脉冲的持续时间为T，则每一目标的回波数据中占据相同的时间间隔T,故压缩前的可分辨能力为：P'=T1、最短可实现信号的TBP近似为1.2、一定带宽下的最短脉冲可近似为sinc函数，其时域能量分布非常集中。3、sinc函数可以通过对矩形函数进行傅里叶逆变换得到，矩形函数的相位必须与正弦波相位对应（即必须是线性的）。线性调频信号的频谱中，相位包含二次分量，为了得到要求的线性相位，可以与含有二

次共轭相位的类似频谱信号相乘，相乘后的信号相位即是线性的。s(t)=s(t)⑤h(t)3.3.2线性调频信号的时域压缩压缩比为初始信号长度除以压缩脉冲的3dB宽度，接近未压缩脉冲的TBPot—ts(t)=rect(t0)exp{j兀K(t—t)2}h(th(t)=exp{-j兀Kt2}s(t)=s(t)⑤h(t)wTsinc(KT(t-1))outr0脉冲分辨率°-0.886〜1「—IKIT~IKIT压缩比等于未压缩脉冲的TBPO非基带信号s(t)=Texp{-j2兀Kt(t-1)}sinc(KT(t-1))outc003.3.4窗效应Sinc函数的PSLR为-13dB,这个PSLR过高，可能会淹没附近的弱目标，降低PSLR的一种方法是对频域匹配滤波器引入平滑窗，以减少主瓣到旁瓣的能量泄露。窗能平滑频谱，降低PSLR，但是会使主瓣展宽，损失了分辨率°Kaiser窗使PSLR降低至-21dB，分辨率则扩展了1.18倍。3.4匹配滤波器的实现第四章合成孔径的概念4.1概述信号处理器是对雷达照射时间内来自某一特定目标的一组信号进行的处理，这相当于产生了一个具有很长孔径的等效天线。4.2SAR几何关系4.2.1属于定义目标、波束覆盖区域、星下点、雷达轨迹：星下点在地球表面上的移动轨迹。平台速度：指平台沿飞行路径的速度，用匕表示；波束速度：指零多普勒线扫过地面的速度，用^巴表示。方位向、零多普勒面：这是一个垂直于平台速度矢量的包含传感器的平面；这个平面与地面的交线为零多普勒线。最短距离、最近位置、零多普勒时刻：传感器与目标最接近的时刻。波束宽度：在方位向，天线的激励是均匀的，波束宽度近似等于波长除以方位向的天线长度。在俯仰向，波束宽度决定了可成像的“距离带”宽度，俯仰向的辐射方向图通常由不均匀孔径形成。目标轨迹、波束中心穿越时刻、信号空间和图像空间、斜距平面、地距、斜视角、距离横向4.2.2卫星地距几何凡+h凡+h2凡g+h)AR和AG的大小可以认为是一个距离分辨率单元。4.3距离等式R%)=用+吁/⑴地球弯曲几何国直线几何4.4SAR距离向信号发射脉冲：S(T)=W(T)COS{2兀fT+KKT2}以长度量纲表示距离分辨率为：c0.886YC1P=W^铝r2IKIT21KIT某目标点的接收信号：S「(T)=只0SptJ(T—2=A)心(T一2死/切XCOfl^2kfo(T—2Ra/c)+JtKr(T—2&/c)2+W}4.5SAR方位向信号如果天线与散射体不断接近，接收信号的频率则会增加。反之，如果两者不断远离，接收信号的频率则会减小。相干脉冲：具有均匀间隔，每个脉冲用Spui(t)=w,(t)cos{2兀fT+兀Kt2}表示。在机载情况下，每个回波可以在脉冲发射间隔内直接接收到。

在目标被波束中心（图中B点）照射之前，接收信号强度一直不断增加。当波束中心穿过目标之后，在目标被波束方向图的第一个零点（图中的C点）照射到之前，信号强度又逐渐减弱。PaW其中9为斜距平面内测得的与视线的夹角，P旅为方向位波束宽度°.886人/、为方向位天线长度。雷达能量的双程传播过程，接收信号的强度由P（9）的平方给出。a吸何）sjncsjnc2（0・886㈤二。螳）ap：{arctans"m）=Ad"（t-27?（杉/c）wa（T）-Tfc）x86｛2时o（*-迎（切/e）+JtKr（T-2J?（ij）/c）2十破｝多普勒中心频率:多普勒带宽:目标照射时间:方位向调频率:=2衅"@r£_2V；2CO33eTiCKa多普勒中心频率:多普勒带宽:目标照射时间:方位向调频率:=2衅"@r£_2V；2CO33eTiCKa=:d羿）4.6二维信号距离时间也称“快时间”，而方位时间称为“慢时间”，因为距离间隔由距离时间和光速得到，而方位间隔由方位时间和比光速慢得多的波束照射区推移速度得到。解调后的基带信号：8。（"）A剧妇&叫（T-2R（m/e）W7?-％）xexp{“4瓦/（）7?3）/c}exp{九玲（―杰⑴）//}SAR冲击响应:无isp（EE）=吓（"27?（叮）上）吨（叮一%）*弥p〈-j4兀/b&S）/c}exp｛扑K『（t-2*弥p〈-j4兀/b&S）/c}exp4.7SAR分辨率与合成孔径SAR处理前的方位向分辨率为波束宽度在地面上的投影，即：/=R（%）8%=。•塑&方）人该式称为真实孔径雷达分辨率。机载情况下为数百米，星载情况下为数千米。由于平台的运动，方位向信号也受到频率调制，因此，也可以通过匹配滤波器得到高的分辨率。方位向分辨率可以写成0.886乘以带宽的倒数：通常，星载情况下，治-1，方位向分辨率可直接写成^=LJ2，这意味着方位向分辨率近似为天线长度的一半。一般雷达或处理前的SAR方位向分辨率由方位向波束宽度确定，波束宽度则由雷达波长人和天线长度L决定。天线的真实孔径、仅在1—15m，但合成孔径在机载情况下可达数百米，星载数千米。合成孔径在下图中用L表示，它是目标在雷达波束照射期间传感器S所经过的路径长度。该长度决定了某一特定目标所能获得的数据量。合成孔径七表示为：L_性乌_0.886RqA*_0,886A(%)AV；；_cos_Lcos农Vq—瓦％_坚86人合成后的半功率波束宽度为钏=FT，式中的系数2来自于雷达的双程传播。假设一个天线的波束宽度为气，则方位向分辨率为R(nc)^?。方位向压缩比：J1总的压缩比：J=仝屈.您瓦上4.8小结距离向的调制源自发射脉冲