RD成像算法分析

1、RD成像算法分析摘 要本文是雷达成像原理的第一次课程作业。在此次作业中我完成了对已提供RD算法的Matlab程序的分析，并且根据自己对该算法的理解重新编写了Matlab程序，用它对RADARSAT-1的真实回波数据进行了成像实验。本文分三部分：第一部分对经典RD算法做了一个简单的介绍；第二部分对已提供RD算法的Matlab程序进行了分析；第三部分用自己编写的RD算法的Matlab程序对RADARSAT-1的真实回波数据进行了成像实验，并给出了实验结果。一、 经典RD算法简介RD（Range Doppler）算法是一种非常经典的SAR成像算法。虽然它有着比较长的历史，但它鲜明的特点使得它能“长盛

2、不衰”，在一些SAR系统中，它是一种实际采用的成像算法。经典的RD算法只考虑了相位展开的一次项，它分为三个步骤：距离压缩（Range Compression，RC）、距离迁移校正（Range Cell Migration Correction，RCMC）和方位压缩（Azimuth Compression，AC）。其中RCMC既可以在Range-Doppler域完成，也可以在Frequency-Azimuth Range域完成。它的流程如Error! Reference source not found.所示。回波数据距离向FFT距离匹配函数距离向IFFT方位向FFT方位匹配函数方位向IFFTS

3、AR图像FFTFFT回波数据距离向FFT距离向IFFT方位向FFT方位向IFFTSAR图像RCMC距离匹配函数方位匹配函数FFTFFT线性相位函数（a）Range-Doppler域RCMC（b）Frequency-Azimuth Range域RCMC图 1 经典RD算法流程图二、 RD算法的Matlab程序分析已提供的RD算法Matlab程序是严格基于Error! Reference source not found.（a）中的流程进行编写的。由于该程序并没有特别新颖的地方，所以本文只利用该程序对RD算法作进一步的理解，而将工作重点主要放在对RADARSAT-1真实回波数据的成像上。该程序可以

4、分为六个小部分：设定参数、模拟回波、RC、RCMC、AC以及结果显示。1、 设定参数该程序是在空域中实现SAR成像的，所以距离是基本的度量单位。程序采用的坐标系如Error! Reference source not found.所示。y/mx/mAzimuthRange图 2 程序采用的坐标系其中一些重要变量和表达式的意义见Error! Reference source not found.（若程序中有明确注释，在此不予列出）。表 1 变量和表达式意义变量/表达式意义Xc测绘带中心距离2*X0测绘带宽2*Y0方位向上接收回波的距离2、 模拟回波程序假定了5个点目标，并且它们对电磁波的反射能力

5、是相同的，它们各自回波的线性迭加就构成了仿真的SAR回波数据，其幅值如Error! Reference source not found.所示。图 3 回波幅值图3、 RC进行RC时先产生距离向上所发射的LFM信号，然后对其作FFT得到发射信号的频谱，根据匹配滤波器的原理，它的共轭就是距离向上匹配滤波器的频谱。将匹配滤波器的频谱与回波没距离向的FFT相乘进行匹配滤波，再将结果作IFFT就完成了RC。其流程见Error! Reference source not found.（a）中相应部分。RC后的数据的幅值图如Error! Reference source not found.所示。图 4

6、RC后的幅值图从Error! Reference source not found.可以明显的看出，RC后的图像在距离向上有明显的迁移，所以要进行RCMC。4、 RCMC由于算法只考虑了相位展开中的一线性项，故RCM只是距离的线性函数。程序是通过Sinc插值来补偿RCM的。插值的运算量是比较大的。RCMC后的幅值图如Error! Reference source not found.所示。图 5 RCMC后的幅值图5、 ACAC的步骤和RC是一样的。也是先产生方位向上的LFM信号，然后进行匹配滤波，其流程见Error! Reference source not found.（a）中相应部分，在

7、这里就不多加表述了。完成这一步后就得到了最终的SAR图像，如Error! Reference source not found.所示。图 6 SAR图像三、 RADARSAT-1真实回波数据成像实验如果要用已提供的程序对RADARSAT-1的真实回波数据进行成像，在参数方面要做较大的修改，为了更好的理解RD算法，本文针对RADARSAT-1回波数据的特点重新编写了RD算法的Matlab程序，并在此基础上进行了成像实验。RADARSAT-1是将回波数据进行自动增益控制（AGC）后，以CEOS的格式进行文件存储的。所以先要从CEOS格式文件中读取数据，并进行AGC还原后才能得到真正的回波数据。在不

8、引起混淆的情况下，下文中的回波数据是指RADARSAT-1的经过以上两个步骤后获得的真正回波数据。得到回波数据后，本文是按Error! Reference source not found.（b）所示的流程来实现RD算法的。算法所需的RADARSAT-1的参数是从Ian Cumming等所著Digital Processing of Synthetic Aperture Radar Data一书及其相关Matlab程序中得到的。本文在Matlab R2006b（Matlab 7.3.0）编写和调试的程序，程序文件名为RADARSAT_RDA_Imaging.m。由于程序中注释比较详尽，在此不对

9、程序本身作过多的说明，只对RD算法本身的关键步骤加以分析。Error! Reference source not found.是回波数据经过RC后的幅值图，从图中可以看到一些明亮的条纹，它们表明了在回波数据中有强散射目标。但是这些条纹并不与方位向平行，而是有一定的倾斜，这种倾斜就是由于距离迁移造成的，所以需要进行RCMC。Error! Reference source not found.是RCMC的幅值图，从图中可以看到Error! Reference source not found.中的斜条纹现在已基本与方位向平行了，说明算法中的RCMC是十分有效的。图 7 RC后的幅值图图 8 RCM

10、C后的幅值图在进行AC时，理论上方位向上的匹配滤波器的调频斜率是与距离的关系的，但是对于星载SAR来说这二者的关系不是很明显。所以在一些应用场合下可以用某一固定距离下的匹配滤波器来对方位向进行脉冲压缩。Error! Reference source not found.是用最近距离下的匹配滤波器来对整个实验数据进行AC后的幅值图。为了作一对比，Error! Reference source not found.是用与距离相对应的匹配滤波器进行AC后的幅值图，可以发现Error! Reference source not found.只是在某些细节上稍稍逊色于Error! Reference source not found.。Error! Reference source not found.是将不同数据对应的图像进行拼接，得到的整个测绘区域的SAR图像，从成像质量上看经典RD算法对星载SAR的成像是比较有效的。