SAR合成孔径雷达图像点目标仿真报告.docx

SAR合成孔径雷达图像点目标仿真报告 〔 附matlab 代码 〕SAR图像点目标仿真报告徐一凡1 SAR 原理简介合成孔径雷达 〔Synthetic Aperture Radar ，简称 SAR〕是一种高辨论率成像雷达技术；它利用脉冲压缩技术获得高的距离向辨论率，利用合成孔径原理获得高的方位向辨论 率，从而获得大面积高辨论率雷达图像；SAR 回波信号经距离向脉冲压缩后，雷达的距离辨论率雷达发射信号带宽准备：rC ，式中 r 表示雷达的距离辨论率， Br 表示雷达发射信号带宽， C 表示光速；同 2Br 样， SAR回波信号经方位向合成孔径后，雷达的方位辨论率雷达方位向的多谱勒带宽决定： ava，式中 a 表示雷达的方位辨论率， Ba 表示雷达方位向多谱勒带宽， va 表示 BaD，其中 D 为 2 方位向 SAR平台速度；在小斜视角的情形下，方位辨论率近似表示为 a 方位向合成孔径的长度；2 SAR 的几何关系雷达位置和波束在地面掩盖区域的简洁几何模型如图 1所示；此次仿真考虑的是正侧视的条带式仿真，也就是说倾斜角为零， SAR波束中心和 SAR平台运动方向垂直的情形；图 1 雷达数据猎取的几何关系建立坐标系 XYZ如图 2 所示，其中 XOY平面为地平面；SAR平台距地平面高 H，以速度 V 沿 X 轴正向匀速飞行； P 点为 SAR平台的位置矢量，设其坐标为 〔x,y,z〕 ； T 点为目标的位置矢量， 设其坐标为 〔xT,yT,zT〕 ；几何关系， 目标与 SAR平台的斜距为：RPT〔xxT〕2〔yyT〕2〔zzT〕2〔1〕 图可知： y0,zH,zT0 ； 令 x其中 v 为平台速度， s 为慢时间变量， vs；假设 xTvs ，其中 s 表示 SAR平台的 x 坐标为 xT 的时刻；再令 r 标与 SAR的垂直斜距，重写 〔1〕 式为：H2yT2， r 表示目PTR〔s;r〕r2v2〔ss0〕2〔2〕R〔s;r〕 就表示任意时刻 s 时，目标与雷达的斜距；一般情形下， vss0r ，于是通过傅里叶技术开放，可将 〔2〕 式可近似写为：v2R〔s;r〕rv〔ss0〕r〔ss0〕2〔3〕2r222 可见，斜距是 s 和 r 的函数，不同的目标， r 也不一样，但当目标距 SAR较远时，在观测带内，可近似认为r 不变，即 rR0 ；图 2：空间几何关系 〔a〕 正视图 〔b〕 侧视图图 2〔a〕 中， Lsar 表示合成孔径长度，它和合成孔径时间 Tsar 的关系是 LsarvTsar ；〔b〕 中，为雷达天线半功率点波束角，为波束轴线与 Z 轴的夹角，即波束视角， Rmin 为近距点距离， Rmax为远距点距离， W为测绘带宽度，它们的关 系为：RminHtg〔2〕RmaxHtg〔2〕〔4〕 WRmaxRmin3SAR的回波信号模型SAR在运动中以确定的周期 〔1/PRF〕 发射和接收信号， 具体过程如图 3 所示；发射机以 l 的时间发射啁啾脉冲，然后切换天线开关接收回波信号；图 3 雷达发射脉冲串的时序当雷达不处于发射状态时，它接收 3 反射回波；发射和接收回波的时间序列如图 4 所示；在机载情形下，每个回波可以在脉冲发射间隔内直接接收到；但是在星载情形下，于距离过大，某个脉冲的回波要经过 6~10 个脉冲间隔才能接收到；这里仿真为了便利，默认为机载情形；图 4 脉冲雷达的发射与接收周期假设 Tr 为 chirp 信号连续时间，下标 r 表示距离向；PRF为重复频率， PRT为重复周期；等于 1/PRF；接收序列中， n2*R〔s;r〕 表示发射第 i 个脉冲时，目标回波相对于发射 C 序列的延时；雷达的发射序列 数学表达式为式 〔5〕:s〔t〕 np〔tn*PRT〕2tp〔t〕rect ejKrtej2fctTr〕 表示矩形信号， Kr 为距离向的 chirp 信号调频率， fc为载频； 式中， rect〔 雷达回波信号发射信号波形，天线 方向图，斜距，目标 RCS，环境等因素共同准备；如不考虑环境因素，就单点目标雷达回波信号可写成式〔6〕 所示：sr〔t〕nwp〔tnPRTn〕 〔6〕其中，表示点目标的雷达散射截面， w 表示点目标天线方向图双向幅度加权， n 表示载机发射第 n 个脉冲时，电磁波再次回到载机时的延 时 n2*R〔s;r〕 ， 带 入 式 〔6〕 中 得 ：Csr〔t〕nwrect〔4tnPRT2R〔s;r〕/C〕Tr 〔7〕 exp[jKr〔tnPRT2R〔s;r〕/C〕2] exp[-jR〔s;r〕]exp[j2fc〔tnPRTn〕] 式〔7〕 就是单点目标回波信号模型，其中， exp[jKr〔tnPRT2R〔s;r〕/C〕2] 是chirp 重量，它准备距离向辨论率； exp[-j4 它准备方位向辨论率； R〔s;r〕] 为多普勒重量； 对于任意一个脉冲，回波信号可表示为式 〔8〕 所示：sr〔t,s〕A0wr〔2R〔s;r〕/C〕wa〔ssc〕exp{j4f0R〔s;r〕/C}exp{jKr〔2R〔s;r〕/C〕}2我们知道， 于 R〔s;r〕 随慢时间 s 的变化而变化， 所以运算机记录到的回波数据储备形式如图 5 所示：图 5 目标照射时间内， 单个点目标回波能量在信号处理器的二维储备器中的轨迹4 距离徙动及校正依据图 2 可知，在倾斜角为零或很小的时候，目标与雷达的瞬时距离为 R〔s;r〕 ，依据几何关系可知， R〔s;r〕r2v2〔ss0〕2 ，依据泰勒级数开放可得：222v2R〔s;r〕rv〔ss0〕r〔ss0〕2 〔9〕 2r式〔9〕 可知，不同慢时间对应着不同的 R〔s;r〕 ，并且是一个双曲线形式或者近似为一个二次形式；如图 5 所示，同一目标的回波储备在运算 机里不在同始终线上， 存在距离徙动； 从而定义距离徙动量：v2R〔s,r〕〔ss0〕22r 为了进行方位向的压缩， 方位向的回波数据必需在同 一条直线上，也就是说必需校正距离徙动 R〔s,r〕 ；式可知，不同的最近距离 r 对应着不同的 R〔s,r〕 ，因此在时域处理距离徙动会特殊麻烦；因此，对方位向进行傅里叶变换，对距离向不进行变换，得到新的域；于方位向的频率即为多普勒频率，所以这个新的域也称为距离多普勒域；将斜距 R 写成多普勒 fa 的函数，即 R〔fa,r〕 ；众所周知，对最近距离为 r 的点目标 P，回波多普勒 fa 是倾斜角的函数，即 fa2Vsin ，斜距 R〔fa,r〕r/cos ，于是R〔fa,r〕r/cosr/1sin2 r/1〔fa2〕2V 〔11〕1 r 2rfa28V 所以距离多普勒域中的我距离徙动为R〔fa,r〕= 2rfa2 ，可发觉它不随慢时间变换，同一最短距离r 对应着相同大小的距离徙动；因此在距离多普勒域对一个距离徙动校正就是对一组具有相同最短距离的点目标的距离徙动校正，这样可以节省运算量； 为了对距离徙动进行校正，需要得到距离徙动单元，即距离徙动表达在储备单元中的移动数值，距离徙动单元可以表示为 R〔fa,r〕/r ，这个值通常为一个分数，于储备单元都是离散的，所以不同通过在储备单元简洁的移动得到精确的值；为了得到精确的徙动校正值，通常需要进行插值运算；本仿真接受了两种插值方法最近邻点插值和 sinc 插值，下面分别进行介绍；最近邻点插值法的优点是简洁而快速，缺点是不够精确；其中 N 为整数部分， R〔fa,r〕/r=Nn ， n 为小数部分，整数部分徙动可以直接通过平移排除，对于小数部分就通过四舍五入的方法变为 0 或者 1，这样就可以得到较为精确的插值；Sinc 插值原理如下：在基带信号下，卷积核是 sinc 函数18Vh〔x〕sinc〔x〕 插值信号为sin〔x〕 x g〔x〕gd〔i〕sinc〔xi〕 i即为全部输入样本的加权平均；可通过频域来懂得，如图 6 所示，采样信号 gd〔i〕 的频谱 Gd〔f〕 等于以采样率重复的信号频谱；为了重建信号 g〔x〕 ，只需要一个周期频谱，因此需要理想矩形低通滤波器在频域中提取基带频谱所示；已 知该理想滤波器在时域中是 sinc 函数；于频域相乘相当于时域卷积，故插值可以通过与 sinc 核的卷积来实现；图 6 理想低通滤波器怎样对采样信号进行插值5 点目标成像 matlab 仿真距离多普勒算法 距离多普勒算法是在 1976 年至 1978年为民用星载 SAR提出的，它兼顾了成熟、简洁、高效和精确等因素，至今仍是使用最广泛的成像算法；它通过距离和方位上的频域操作，到达了高效的模块化处理要求，同时又具有了一维操作的简便性； 图 7 示意了 RDA的处理流程；这里主要争辩小倾斜角及短孔径下的基本 RDA处理框图； 1.当数据处在方位时域时，可通过快速卷积进行距离压缩；也就是说，距离 FFT 后立刻进行距离向匹配滤波，再利用距离IFFT 完成距离压缩；回波信号为：s0〔t,s〕A0wr[t2R〔s〕/c]wa〔ssc〕 exp{-j4f0R〔s〕/c}exp{jKr〔t-2R〔s〕/c〕} 距离向压缩后的信号为：2src〔t,s〕IFFTt{S0〔ft,s〕H〔ft〕}A0r[t2R〔s〕/c]wa〔ssc〕exp{j4f0R〔s〕/c} ff2H〔ft〕rect{}exp{j}exp{j2ft0}|K|TK 2. 通过方位 FFT 将数据变换至距离多普勒域，多普勒中心频率估量以及大部分后续操作都在该域进行；方位向傅里叶变换后信号为：S1〔t,fs〕FFTs{src〔t,s〕}2Rrd〔fs〕]Wa〔fsfsc〕 c4f0R0fs2 exp{-j}exp{j}cKa A0pr[t 3. 在距离多普勒域进行随距离时间及方位频率变化的 RCM，C该域中同一距离上的一组目标轨迹相互重合； RCMC将距离徙动曲线拉直到与方位 频率轴平行的方向； 这里可以接受最近邻点插值法或者 sinc 插值法，具体插值方法见前面；假设 RCMC插值是精确的，信号变为：2R0〕Wa〔fsfsc〕c4f0R0fs2 exp{-j}exp{j}cKaS2〔t,fs〕A0pr〔t4. 通过每一距离门上的频域匹配滤波实现方位压缩；为进行方位压缩，将 RCMC后的S2〔t,fs〕 乘以频域匹配滤波器 Haz〔fs〕 ；fs2Haz〔fs〕exp{j} KaS3〔t,fs〕S2〔t,fs〕Haz〔fs〕 A0pr〔t2R0/c〕Wa〔fsfsc〕exp{j4f0R0}c 5. 最终通过方位 IFFT 将数据变换回时域，得到压缩后的复图像；复原后的图像为：sac〔t,s〕IFFTs{S3〔t,fs〕} A0pr〔t-2R0/c〕pa〔s〕 exp{-j4f0R0}exp{j2fscs}c图 8 距离多普勒算法流程图Chirp Scaling 算法 距离多普勒算法具有诸多优点，但是距离多普勒算法有两点不足：第一，当用较长的核函数 提高距离徙动校正精度时，运算量较大；其次，二次距离压 缩对方位频率的依靠性问题较难解决，从而限制了其对某些大斜视角和长孔径 SAR的处理精度； Chirp Scaling 算法防止了 RCMC中的插值操作，通过对 Chirp 信号进行频率调制，实现了对该信号的尺度变换或平移； 图 8 显示了 ChirpScaling 算法处理。