技术领域
[0001] 本
发明属于雷达技术领域，更进一步涉及雷达成像技术领域中的高
分辨率宽测绘带
合成孔径雷达(Synthetic Aperture Radar，SAR)地面运动目标成像方法。
背景技术
[0002] 合成孔径雷达因具有全天时、全天候、远
距离成像的能
力，在战略防御、地形测绘等军用和民用领域得到了广泛应用。其中，利用SAR获取运动目标成像结果已成为现代雷达的研究热点之一。
[0003] 同时实现高分辨宽测绘带成像是雷达成像技术的追求目标。然而，传统星载SAR系统受最小天线面积限制，无法同时获取高分辨宽测绘带图像。结合数字波束形成技术，方位多通道SAR系统通过发射低脉冲重复
频率(Pulse Repetition Frequency,PRF)
信号可克服此限制，但是低PRF将导致回波信号产生严重的多普勒模糊。针对地面静止目标，现有多种方法可实现多普勒模糊抑制从而实现高分辨宽测绘带SAR成像。但是，对于地面运动目标来说，由于其运动参数未知，无法直接精确构造运动目标的模糊重构函数，从而导致处理性能下降。
[0004] 针对以上运动目标成像存在的问题，各国学者提出了多种方法。Li等人提出通过搜索空时平面内的运动目标区，得到运动目标所有可能的谱分量，但并未给出相应的判断准则。Baumgartner等人根据最大信杂噪比准则搜索运动目标的速度，但由于该方法需要对所有可能的运动目标速度都进行一次全孔径成像操作，因此运算量会大大增加。Yang和Zhang等人假设运动目标谱分量位于杂波谱分量中间，尽量可能地正确约束运动目标谱分量，虽然该方法可以避免繁重的搜索操作，但当运动目标的运动方向偏离该导向方向时，该方法将无法正确对运动目标的谱分量进行约束，从而无法对模糊分量进行抑制并导致运动目标增益的损失。
发明内容
[0005] 本发明针对上述高分辨宽测绘带SAR地面运动目标成像方法面临的问题，提出了一种基于稀疏空间谱估计的高分辨宽测绘带SAR地面运动目标成像方法。与传统方法不同，本发明没有采用搜索或近似假设，并能以较小的运算量实现对地面运动目标的高
质量成像。
[0006] 为实现上述目的，本发明的主要步骤如下：
[0007] (1)对各通道接收的回波信号进行方位向傅里叶变换；
[0008] (2)根据稀疏空间谱估计方法，利用某一多普勒单元回波估计得到运动目标的速度；
[0009] (3)根据估计得到的运动目标速度，构造运动目标的
频谱重构函数，得到运动目标的无模糊全带宽信号；
[0010] (4)根据估计得到的运动目标速度，对无模糊的全带宽回波信号进行传统高分辨运动目标成像，得到地面运动目标高分辨率宽测绘带图像。
[0011] 本发明与现有的技术相比具有以下优点：
[0012] 本发明利用稀疏空间谱估计的方法直接对运动目标的速度进行精确估计，从而避免了传统搜索操作引入的巨大运算量，另外运动目标速度的精确估计在保证运动目标的多普勒频谱精确重构的同时，克服了传统方法带来的运动目标谱分量增益的损失。
附图说明
[0014] 图2为在斜距平面下的SAR系统观测几何示意图；
[0015] 图3为利用回波信号的某一多普勒单元进行空间谱估计的结果；
[0016] 图4为利用传统方法得到的运动目标仿真成像结果；
[0017] 图5为利用本发明方法的仿真成像结果。
具体实施方式
[0018] 参照附图1，本发明的具体实施步骤如下：
[0019] 步骤1，对各通道接收的回波信号进行方位向傅里叶变换。
[0020] 假设SAR系统共有M个沿航向均匀分布的接收通道，如图2所示，x轴为平台运动速度方向，y轴为斜距方向，雷达到地面目标的斜距矢量与零多普勒面的夹
角φ称为锥角。不失一般性，本发明假设SAR系统为正侧视成像。对于方位多通道高分辨宽测绘带SAR系统来说，各通道接收回波均可等效为参考接收通道(如第1个通道)接收回波的时延。为保证宽距离测绘带，SAR系统采用低PRF，这样回波信号将发生多普勒模糊。考虑加性白噪声，第m个通道接收的回波信号进行方位向傅里叶变换后可写为
[0021]
[0022] 其中，m＝1,2,…,M，M为方位接收通道个数，τ为距离时间，fp为系统PRF，fd∈[-fp/2,fp/2]为多普勒频率，vs为平台运动速度，a(τ,fd)为第1个通道的接收回波，Δxm为第m个通道相对第1个通道的距离，L＝(N-1)/2，N是回波信号的多普勒模糊数，nm(τ,fd)为第m个通道的噪声。地面运动目标回波的多普勒频率与锥角φ的关系如下[0023]
[0024] 其中，λ为
载波频率，vt为运动目标的垂直航向速度，
[0025] 忽略τ，式(1)可用矢量形式表示为
[0026]
[0027] 其中，T
[0028] s(τ,fd)＝[s1(τ,fd),s2(τ,fd),…,sM(τ,fd)]， (4)T
[0029] a(τ,fd)＝[a(τ,fd-L·fp),a(τ,fd-(L-1)·fp),…,a(τ,fd+L·fp)]， (5)[0030] P(fd)＝[p-L(fd),p-L+1(fd),…,pL(fd)]， (6)
[0031]
[0032] n(τ,fd)＝[n1(τ,fd),n2(τ,fd),…,nM(τ,fd)]T， (8)[0033] [□]T表示矩阵转置。由此可知，各通道接收的多普勒模糊回波信号可等效看作为来自不同方向的信号
叠加。与地面固定杂波不同，由于运动目标的速度vt未知，所以无法精确获知阵列流形矩阵P(fd)，从而无法直接利用已有方法对回波信号进行频谱重构。
[0034] 步骤2，根据稀疏空间谱估计方法，利用某个多普勒单元回波信号估计得到运动目标的速度。
[0035] 根据式(2)所示的多普勒频率与锥角的关系，本发明将运动目标速度估计问题转换为信号的
波达方向(direction-of-arrival,DOA)估计问题，也即空间谱估计问题。根据SAR系统参数可计算出多普勒模糊数N，也即待估计DOA信号的个数。更进一步，由于运动目标空间谱在空间域上为稀疏分布，因此上述DOA估计问题可转换为以下稀疏问题[0036]
[0037] 其中，
[0038] Y＝[s(τ1,fd),…,s(τK,fd)]， (10)
[0039]
[0040]
[0041]
[0042] N＝[n(τ1,fd),…,n(τK,fd)]， (14)
[0043] Y为各通道接收的快拍信号，K为快拍个数，在实际中，可以利用相邻的K个距离单元回波替代。 为所有可能方向的回波信号构成的阵列流形，φq(q＝1,…,Q)的取值范围为雷达波束的照射范围。 为Q×1维稀疏向量，若第n(n＝1,…,N)个信号的方向角为φq，则 的第q个值为a(τ,fd-ln·fp)，其中ln＝n-L-1，否则其值为零。一般来说，Q□M＞N。由此可将代价函数写为
[0044]
[0045] 其中，||·||f为Frobenius范数， 下标q表示向量的第q个值，λ为稀疏规则化参数。求解上述优化问题即可得到 然后根据波峰
位置得到信号的波达方向，再根据下式计算出运动目标的速度
[0046]
[0047] 其中，fd为所选取的多普勒单元，φq为估计得到的波达方向。在实际操作中，为了减小噪声的影响，可以选择多个多普勒单元估计得到多个vt，然后将其平均求得最终的vt。
[0048] 式(15)所示的优化问题为凸优化问题，现已有多种成熟的求解方法。需要说明的是，对于SAR系统来说，M和N的值都较小，而且在实际操作中，仅需要对式(15)进行若干几次求解，因此相比传统的搜索方法，本发明的运算量将大大降低。
[0049] 步骤3，根据上一步估计得到的运动目标速度，构造运动目标的多普勒频谱重构函数，得到运动目标的无模糊全带宽信号。
[0050] 利用步骤2估计得到的运动目标速度，即可正确构造出运动目标回波信号的阵列流形P(fd)，然后可采用传统的模糊重构方法得到无模糊的运动目标全带宽信号。
[0051] 步骤4，根据估计得到的运动目标速度，对无模糊的全带宽回波信号进行传统高分辨运动目标成像，得到地面运动目标高分辨率宽测绘带图像。
[0052] 利用步骤2估计得到的运动目标速度，对回波信号进行距离压缩和距离徙动校正，然后再进行方位压缩，得到地面运动目标的高分辨宽测绘带SAR图像。
[0053] 下面结合仿真数据实验对本发明的效果做进一步的说明。
[0054] 1、仿真条件：
[0055] 星载方位多通道SAR系统的仿真参数如下表所示。多普勒模糊次数约为3，由于运动目标存在垂直航向速度，其多普勒频谱将偏离杂波谱约682Hz。
[0056]
[0057] 2、仿真数据包实验分析：
[0058] 首先对各通道回波信号进行方位向傅里叶变换，然后选取某一多普勒单元回波进行稀疏空间谱估计。图3给出了利用第1360个多普勒单元回波进行空间谱估计的结果，相邻的10个距离单元作为快拍数。可以看出，利用本发明的方法可以准确估计得到运动目标的空间谱，由此估计得到的目标垂直航向速度为10.53m/s，对应的多普勒偏移约为677Hz，与真实值非常接近。图4给出了假设运动目标谱偏离杂波谱fp/2所构造的频谱重构函数得到的运动目标成像结果的方位剖面图，可以看出成像结果存在较为严重的方位模糊。图5(a)给出了利用本发明所述的方法得到的运动目标成像结果，图5(b)为成像结果的局部放大图，可以看出运动目标的多普勒模糊分量得到了有效的抑制，并实现了很好的聚焦。
