频域合成孔径雷达成像方法
【技术领域】
[0001] 本发明涉及合成孔径雷达技术领域，尤其一种频域合成孔径雷达成像方法。
【背景技术】
[0002] 合成孔径雷达(SAR)是一种工作于主动探测方式的微波成像遥感设备。SAR装载在 飞机、卫星等运动平台上，随着平台的移动发射探测信号并接收目标的回波信号，通过将一 系列接收到的信号进行综合处理使其等效于一个大孔径天线发射和接收到的信号的方式 来实现高分辨率的目标探测，这一技术称为孔径合成技术。
[0003] 合成孔径雷达具有探测距离远、成像区域宽、分辨率高，可全天候昼夜工作。在资 源普查、测绘、海洋环境调查、灾害监测、军事侦察等领域具有广泛的应用。
[0004] 雷达发射的电磁波为球面波，在低精度成像时，可近似为平面波进行成像处理。而 在高精度成像时，不能忽略电磁波存在的波前弯曲。在波束指向偏离正侧视时也会存在波 前与平台运动方向不平行的问题，上述问题统称为距离迀移。
[0005] 通常，为了保证雷达的高精度成像性能，需要在成像处理中加进距离迀移校正 (RCMC)处理。目前，主要的成像方法包括距离多普勒域成像方法(RD方法）、波数域成像方法 (ω -k方法）、Chirp Scaling方法(CS方法）AD以及ω -k等方法具有处理简单、精度高等特 点，然而由于在处理的过程中需要进行插值，带来了插值误差，为了减少插值带来的误差则 需要采用更复杂的插值算法，又带来了运算量增加的问题。因此，在高精度的SAR成像处理 中，通常采样CS方法。CS方法避免了插值运算，但是带来了处理步骤增多、需要多次转置存 储、处理模型存在近似误差等问题，增加了处理的复杂性和局限性。
[0006] 在不断的实践过程中，申请发现目前已有的SAR成像方法中，存在着或者精度低， 或者运算处理复杂、运算量大等问题。因此，本领域存在对进一步简化处理和提高精度的 SAR成像方法的需要。
【发明内容】
[0007] (一)要解决的技术问题
[0008] 鉴于上述技术问题，本发明提供了一种频域合成孔径雷达成像，以满足对降低计 算复杂性、降低运算量、提高成像精度等要求。
[0009] (二)技术方案
[0010] 本发明频域合成孔径雷达成像包括:步骤A，对合成孔径雷达接收的回波信号依次 进行下变频和解调后，采样得到信号s(ta，t r)，对信号s(ta，tr)进行两维傅里叶变换，得到S (fa，fr)，其中，1为方位向变量，tr为距离向变量，fa为方位向变量对应的频率变量；f r为距 离向变量对应的频率变量;步骤B，将信号S (fa，fr)与雷达系统工作参考信号HRF (fa，fr)相 乘，得到信号步骤C，将信号按照设定的插值算法进行Stolt变换，得到 信号S 2(fa，f/)，其中，f/为新的距离频率变量；以及步骤D:将信号S2(f a，f/)进行二维逆 傅里叶变换后，与设定的插值算法对应的插值误差补偿函数相乘，得到高精度合成孔径雷 达图像SI(ta，tr) 〇
[0011] (三)有益效果
[0012] 从上述技术方案可以看出，本发明频率合成孔径雷达成像具有以下有益效果：
[0013] (1)采用频域精确、无误差的成像算法，计算过程中的插值运算带来的误差采用插 值误差补偿方法加以控制，可避免插值计算误差对SAR成像的影响，实现高精度成像；
[0014] (2)插值误差补偿函数由插值算法确定，可以选用便于实现的插值算法，大大简化 了插值运算量；
[0015] (3)实现SAR成像的处理算法及流程简单，可大量节省转置处理的时间和对转置存 储器的需求，降低了成像算法实现的复杂性。
[0016] 从上述特点中可以看出，本发明给出的SAR高精度成像方法，对处理误差进行了有 针对性的补偿、成像精度高，计算简单、运算量少，简化了处理流程、占用的转置处理时间和 存储空间少，具有适合于计算机通过数值计算软件来实现、以及采用数字信号处理器(DSP) 以及FPGA芯片硬件电路实现等优势，具有较好的应用前景。
【附图说明】
[0017] 图1是根据本发明实施频域合成孔径雷达成像方法的流程图。
【具体实施方式】
[0018] 本发明采用频域精确、无误差的成像算法，计算过程中的插值运算带来的误差采 用插值误差补偿方法加以控制，以降低计算复杂性、降低运算量和成像精度。
[0019] 为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照 附图，对本发明进一步详细说明。
[0020] 在本发明的一个示例性实施例中，提供了一种频域合成孔径雷达成像方法。如图1 所示，本实施例频域合成孔径雷达成像方法包括：
[0021] 步骤Α:对合成孔径雷达接收的回波信号依次进行下变频和解调后，采样得到信号 8(^，^)，在距离向补零后进行两维傅里叶变换，得到3(匕，匕）；
[0022] 假设合成孔径雷达发射的线性调频信号为：
[0024] 其中：tr表示距离时间变量，W( ·)表示发射信号的包络，fQ是发射信号的载频，Kr 是发射信号的线性调频率，其中j为虚数单位
[0025] 雷达接收到的回波信号经过下变频、解调后的信号为S(ta，tr)，其中ta为方位向变 量，tr为距呙向变量。
[0026] 为了避免后续频域插值引起模糊，将信号S(ta，tr)在距离向补零至原来信号长度 的2~4倍以上，之后进行距离向和方位向二维傅里叶变换，得到信号S(f a，fr)，其中，fa为方 位向变量对应的频率变量;fr为距离向变量对应的频率变量。
[0027] 步骤B:将信号S(fa，fr)与雷达系统工作参考信号HRF(f a，fr)相乘，得到信号Si(fa， fr)；
[0028]雷达系统工作参考信号为：
[0030] 上式（2)中，c为光速，v为雷达平台的运动速度，rref为进行成像处理的参考距离， 通常选为雷达成像观测区域距离向中心位置与雷达所在位置的距离。
[0031] 将S(fa，fr)与参考函数相乘，得到如下结果：
[0032] Sl(fa,fr)=S(fa,fr) * HRF(fa,fr) (3)
[0033] 步骤C:将信号SKfufr)按照设定的插值算法进行stolt变换，得到信号S2(fa， f·/ )，其中，f/为新的距离频率变量；
[0034] 进行stolt变换，将原来的距离频率变量fr映射为新的距离频率变量f/，如下所 示：
[0038]在通过数字信号处理方式进行上式的stolt变换时，会遇到计算的信号取值在离 散采样点之间的问题，需要对信号进行插值。
[0039]设定的插值算法是：首先在信号采样点之间插入m-1个零值点（相当于m倍增采 样），之后用插值算法对应的卷积函数与插入零值点的信号进行卷积运算从而完成插值。其 中，最邻域插值算法对应长度为m的方波卷积函数，线性插值算