一种机载扫描雷达实时超分辨成像方法与流程

1.本发明属于雷达信号处理技术领域，特别涉及一种扫描雷达技术中的实时快速成像技术。


背景技术：

2.机载雷达斜前视高分辨率成像具备全天时、全天候工作的特点，作为一种主动式的遥感手段在地质测绘、灾害监测、军事侦察等领域有着广泛的应用。机载雷达往往采用多普勒波束锐化技术(dbs)对斜前视区成像，其主要利用雷达与目标之间的相对运动所引起的多普勒信息区分目标。但是，在宽幅快扫的成像模式中，基于传统快速傅里叶变换的dbs技术往往因其观测孔径不足而导致方位分辨率受限，其点扩散函数往往呈现宽主瓣、高旁瓣的特点。
3.为进一步提升方位分辨率，文献“super
‑
resolution doppler beam sharpening imaging based on an iterative adaptive approach，remote sensing letters,vol.7,pp.259
–
268,01 2016.”提出了一种基于迭代自适应方法(iaa)谱估计方法用于解决dbs目标散射系数估计问题。但是，由于每次迭代的矩阵求逆运算，该方法的计算复杂度非常高，不利于实时处理。文献“super
‑
resolution doppler beam sharpening method using fast iterative adaptive approach based spectral estimation，journal of applied remote sensing,vol.12,no.1,pp.1,2018.”提出了一种基于加权的快速迭代自适应超分辨率频谱估计方法提出了最小二乘理论以增强dbs的分辨率。然而，该方法确实有效地降低了迭代自适应方法的计算复杂度，但其对于回波的处理方式仍采用批处理形式。这意味着该方法无法在扫描的同时进行超分辨处理，从而导致实时性能的损失。


技术实现要素：

4.为解决上述技术问题，本发明提出一种机载扫描雷达实时超分辨成像方法，不仅有效地改善了传统方法的方位分辨率，而且显著降低了其复杂度和内存占用率，有利于机载雷达的高品质连续实时成像。
5.本发明采用的技术方案为：一种机载扫描雷达实时超分辨成像方法，包括：
6.a1、建立斜前视矢量卷类积精确回波模型；
7.a2、构造描述天线方向图和多普勒导向矩阵与雷达回波信号的可解析最小二乘代价函数，将斜前视成像分辨率提升问题转化为优化估计问题；
8.a3、针对首个窗口数据，采用基于正则化的方法进行目标估计；
9.a4、针对剩下的窗口数据，采用updating和downdating滑窗递归，实现斜前视区域目标散射系数的实时更新。
10.步骤a1所述斜前视矢量卷类积精确回波模型为：
11.y＝aσ e
12.a＝d
⊙
h
13.其中，y表示方位回波矩阵，表示复数集合，a表示导向矩阵，d表示多普勒矩阵，h表示天线方向图矩阵，σ是目标的散射系数，e是加性高斯噪声矩阵，
⊙
表示点积。
14.y通过下式离散采样得到：
[0015][0016]
其中，τ和t分别表示距离和方位时间变量，σ(x,y)和h(
·
)分别表示目标后向散射系数和天线方向图响应函数，sinc(
·
)表示脉冲压缩响应函数，φ为感兴趣的观测区域，b为带宽,λ为波长。
[0017]
步骤a2所述的代价函数表达式为：
[0018][0019]
其中，y
1:l
表示首次滑窗窗口的回波矢量，a
1:l
表示首次滑窗窗口对应的导向矩阵。
[0020]
步骤a3首个窗口数据估计得到的目标散射系数最优解为：
[0021][0022]
步骤a4具体为：
[0023]
a41、引入矩阵
[0024][0025]
其中i表示单位矩阵；
[0026]
a42、downdating的结果可表示为
[0027][0028]
其中
[0029][0030][0031][0032]
a43、当前的目标散射系数由downdating结果通过updating递归计算得出
[0033][0034]
其中
[0035][0036][0037][0038]
a44、对于n>l重复步骤a42
‑
a43，实现当前第n时刻的在线递归更新结果
[0039]
本发明的有益效果：本发明通过建立斜前视矢量卷类积精确回波模型，突破斜前
视角分辨率提升机理性瓶颈；然后构造描述天线方向图和多普勒导向矩阵与雷达回波信号的可解析最小二乘代价函数，将斜前视成像分辨率提升问题转化为优化估计问题；最后，通过updating和downdating滑窗递归，实现斜前视区域目标散射系数的实时更新。相对于传统方法，提升了方位分辨率，同时显著降低了其复杂度和内存占用率。
附图说明
[0040]
图1为本发明的方法流程图；
[0041]
图2为本发明实施例提供的机载雷达斜前视超分辨成像方法的几何构型；
[0042]
图3为本发明方法与传统方法的距离剖面一维对比结果；
[0043]
图4为本发明的成像结果与现有技术的对比；
[0044]
其中，图4(a)为实波束的二维成像结果，图4(b)为dbs方法的二维成像结果，图4(c)为ls方法的二维成像结果，图4(d)为本发明方法的二维成像结果。
具体实施方式
[0045]
为便于本领域技术人员理解本发明的技术内容，下面结合附图对本发明内容进一步阐释。
[0046]
如图1所示，本发明的方法包括以下步骤：
[0047]
步骤一：建立机载雷达斜前视回波精确模型
[0048]
采用如图2所示机载雷达斜前视成像几何模型，选取如表1所示的雷达仿真系统参数。信噪比设置为snr＝20db。设机载雷达平台以恒定速度v＝30m/s沿着x轴方向匀速飞行，雷达天线以角速度ω＝90
°
/s扫描斜前视区域。平台的飞行高度为h＝1000m平台和目标p的初始斜距为r0＝3000m。在t时刻，平台与目标的距离历史r(t)表示为
[0049][0050]
其中，表示目标俯仰角，θ表示目标方位角。
[0051]
表1仿真参数
[0052]
参数数值工作频段f
c
30.75ghz时宽t
p
40mhz带宽b2μs扫描区域
‑
25
°
～
‑5°
扫描速度ω90
°
/s波束宽度δθ4
°
脉冲重复频率prf2khz初始斜距r03000m机载平台飞行速度v30m/s
[0053]
雷达通过发射具有较大时间带宽积的线性调频信号，实现了高距离分辨率。经过距离预处理，即脉冲压缩和徙动校正，回波信号可以写为，
[0054][0055]
其中，τ和t分别表示距离和方位时间变量，σ(x,y)表示处于观测区域φ内坐标(x,y)处的目标后向散射系数，h(
·
)表示天线方向图响应函数，sinc(
·
)表示脉冲压缩响应函数，φ为感兴趣的观测区域，c表示电磁波传播速度，带宽为b＝40mhz,波长为λ＝0.01m。
[0056]
仿真设定回波的方位采样点数为n＝1024，对于某一距离单元，回波的离散矩阵形式可以由上式导出为
[0057]
y＝aσ e
[0058]
a＝d
⊙
h
[0059]
其中，表示方位回波，表示复数集合，表示导向矩阵，由多普勒矩阵和天线方向图矩阵组成。仿真中设定的方位目标采样点数k＝670。是目标的散射系数，是加性高斯噪声矩阵，
⊙
表示点积。
[0060]
本发明通过建立斜前视矢量卷类积精确回波模型，突破斜前视角分辨率提升机理性瓶颈。
[0061]
实波束回波的点目标剖面结果如附图3中实线所示，二维目标成像结果附图4(a)所示。
[0062]
步骤二：针对首个窗口数据处理，采用基于正则化的方法进行目标估计
[0063]
设滑动窗口长度与天线方向图长度相等，即l＝120。定义滑窗窗口内的回波和导向矩阵的矢量形式分别为
[0064][0065][0066]
其中，y
n
‑
l 1:n
表示包含从n
‑
l 1时刻到n时刻的滑动窗口长度为120的离散回波矢量，y(n)表示y
n
‑
l 1:n
中的第n个元素，a
n
‑
l 1:n
表示滑窗窗口内的导向矩阵，a(n)表示矩阵a
n
‑
l 1:n
的第n行矢量，此外，添加正则化项γ抑制噪声。则目标散射估计问题可以被重建为一个正则化约束下的最优化问题
[0067][0068]
||
·
||2表示矩阵的二范数运算；
[0069]
对于首次递归，令n＝l，根据上述优化问题构建代价函数
[0070][0071]
其中，表示对σ的估计值，表示首次滑窗窗口的回波矢量，表示首次滑窗窗口对应的导向矩阵。上式是对于首次目标散射的凸优化问
题，最优解可以表示为
[0072][0073]
其中，γ表示正则化参数，用以权衡噪声和方位分辨率。
[0074]
步骤三：基于downdating和updating的滑动窗口数据在线递归更新
[0075]
在仿真递归更新算法计算时，选取正则化参数为γ＝0.1，对于n>l＝120的方位回波反演可以通过结合downdating和updating递归分别实现和和的递归更新，最终得到从到递归更新算法。
[0076]
首先，引入矩阵
[0077][0078]
其中，表示恒等于，上标h表示共轭转置，i表示单位矩阵。
[0079]
downdating的结果可表示为
[0080][0081]
其中
[0082][0083][0084][0085]
这里的计算采用a
n
来表示a(n)；上标中的l表示downdating递归时的维度；
[0086]
当前的目标散射系数即下一次更新的结果，可以由downdating结果通过updating递归计算得出
[0087][0088]
其中
[0089][0090][0091][0092]
上标中的l
‑
1表示updating递归时的维度；即updating递归与downdating递归的维度不相同。
[0093]
对于n＝l 1,l 2,...，依次遍历上述迭代公式直到波束最终方位，得到第n次在线递归更新时目标散射系数使用本发明提出的滑窗递归实时成像方法对目标进行递归求解，对于一个点目标的距离剖面的超分辨结果如附图3所示(三角型线)，其中虚线表示传统dbs方法的结果，星号线表示传统ls方法的结果。可以看到，本发明方法的分辨率明显优于传统方法。
[0094]
如图4(b)所示为传统dbs方法的二维成像结果，如图4(c)所示为ls方法的二维成像结果。与dbs方法相比，其分辨率有所提高。但是，该方法具有较高的复杂度和较差的实时
性能。图4(d)给出了本发明方法的成像结果，该结果为多个滑动窗口结果的叠加显示。可以看出，本发明方法的结果明显优于传统的dbs方法，并与ls方法的分辨率性能相近。在实际应用中，本发明方法在算法复杂度，实时性能和内存占用方面均优于传统ls方法。
[0095]
综上所述，本发明方法相较于传统成像方法，有效地提升了方位分辨率，同时显著降低了算法复杂度和成像处理器的内存占用率。
[0096]
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。