一种适用于存在遮挡效应的旋转目标稀疏多普勒成像方法
【技术领域】
[0001] 本发明属于雷达成像技术领域，涉及一种适用于存在遮挡效应的旋转目标稀疏多 普勒成像方法，可用于对存在遮挡效应的旋转目标进行成像。
【背景技术】
[0002] 多普勒成像方法能够提供目标在距离维和横向距离维的精确信息，是实现对3维 旋转目标进行2维成像的一种常用方法。针对不同的旋转目标设计有效的多普勒成像方法 在实际应用中具有重要意义。例如，对太空碎片进行多普勒成像是空间碎片监测系统中的 关键技术之一。目前，随着人类空间活动的日益频繁，空间碎片的数量也在与日俱增，已成 为了外太空中极具危害的一种污染源。它们不仅拥有庞大的数量，还具有很高的运动速度， 通常在近地轨道上能够达到l〇km/ S，从而极易与在轨飞行器发生碰撞，导致大灾难的发生。 有效的太空碎片多普勒成像在防撞指示、碰撞程度预测以及碎片的回收清理等工作中均具 有重要的指导意义。
[0003] 在多普勒成像过程中，目标的雷达回波在整个观测时间段内都被看作位于单个距 离单元内，并且被一系列的雷达发射脉冲序列采样。根据散射中心模型理论，多普勒成像 的关键在于估计目标上主要散射中心的位置和后向散射强度，并且估计精度需要尽可能地 高。仅用多普勒信息的传统单快拍成像方法，其假设在整个成像过程中观测目标上所有散 射中心始终对于雷达是可见的。但在实际中，当目标的自旋轴与雷达视线不一致时，会导 致目标上的部分散射中心在部分观测间隔内对于雷达是不可见的，即遮挡效应，此时的遮 挡效应相当于对目标在方位向进行了非均匀采样，这个因素在已有成像算法中没有被涉及 到。
[0004] 一般而言，雷达成像的过程可看作是由测量场景的电磁散射信号来反推场景信息 的过程。根据雷达成像的数学模型，可以将雷达成像的过程看成是一个逆问题。由于遮挡 效应的存在，测量数据的信息量通常是有限的，无法满足人们对场景信息详细程度的要求。 所以，上述逆问题是病态的，无法通过经典的最小二乘估计方法得以很好地解决。传统的多 普勒成像方法通过对最小二乘估计中不可逆或不稳定的部分作近似处理来实现成像，但代 价是具有一定的主瓣宽度和旁瓣效应，即造成场景图像中的细节信息有所损失。
【发明内容】
[0005] 本发明的目的在于克服上述已有技术的不足，提出一种适用于存在遮挡效应的旋 转目标稀疏多普勒成像方法，可实现对存在遮挡效应的旋转目标进行成像，并且能够自适 应地估计遮挡区域。
[0006] 本发明的实现思路是：将对存在遮挡效应的旋转目标进行多普勒成像的过程看作 是求解一个逆问题的过程，利用目标上的散射中心是稀疏分布的（观测场景具有一定的稀 疏性）以及目标上各个散射中心的后向散射现象随时间缓慢变化（慢变）这两个假设，在 多普勒成像过程中增加相应的稀疏约束，将病态的逆问题转化成良态的2维Fused Lasso 模型来求解，从而实现成像。
[0007] 为达到上述技术目的，本发明采用以下技术方案予以实现。
[0008] -种适用于存在遮挡效应的旋转目标稀疏多普勒成像方法，其特征在于，包括以 下步骤：
[0009] 步骤1，建立存在遮挡效应的旋转目标的时变雷达回波模型；
[0010] 步骤2,对存在遮挡效应的旋转目标的时变雷达回波模型进行稀疏信号表示；
[0011] 步骤3,根据时变雷达回波模型的稀疏信号表示形式，将对存在遮挡效应的旋转目 标的成像过程转化为2维Fused Lasso模型的稀疏优化过程；
[0012] 步骤4,利用凸优化工具包CVX来求解2维Fused Lasso模型的稀疏优化过程，得 到稀疏后向散射强度矩阵B的估计矩阵B ;
[0013] 步骤5,取估计矩阵?中各行的最大值构成等效后向散射强度向量，即实现成像。
[0014] 本发明的有益效果是：将对存在遮挡效应的旋转目标进行多普勒成像的过程看作 是求解逆问题的过程，利用目标上的散射中心是稀疏分布的以及目标上各个散射中心的后 向散射现象具有慢变特性这两个假设，在多普勒成像过程中引入稀疏约束，将病态的逆问 题转化成良态的2维Fused Lasso模型来求解，从而实现对存在遮挡效应的旋转目标的多 普勒成像，同时自适应地估计目标的遮挡区域。
【附图说明】
[0015] 下面结合【附图说明】和【具体实施方式】对本发明作进一步详细说明。
[0016] 图1是本发明的流程图。
[0017] 图2是用本发明对不同遮挡比例下的旋转目标的成像结果图，其中：
[0018] 图2a是本发明对遮挡比例为30%的旋转目标的成像结果图；
[0019] 图2b是本发明对遮挡比例为60%的旋转目标的成像结果图；
[0020] 图2c是本发明对遮挡比例为90%的旋转目标的成像结果图；
[0021] 图2d是本发明对遮挡比例为95%的旋转目标的成像结果图。
[0022] 图3是本发明对不同信噪比下遮挡比例为50%的旋转目标的成像结果图，其中：
[0023] 图3a是本发明对20dB信噪比下遮挡比例为50%的旋转目标的成像结果图；
[0024] 图3b是本发明对15dB信噪比下遮挡比例为50%的旋转目标的成像结果图；
[0025] 图3c是本发明对IOdB信噪比下遮挡比例为50%的旋转目标的成像结果图。
[0026] 图4是电磁仿真的目标模型示意图。
[0027] 图5是电磁仿真目标的一维高分辨距离像随时间的变化图。
[0028] 图6是本发明对电磁仿真数据的成像结果图。
【具体实施方式】
[0029] 参照图1，本发明的适用于存在遮挡效应的旋转目标稀疏多普勒成像方法，包括以 下步骤：
[0030] 步骤1，建立存在遮挡效应的旋转目标的时变雷达回波模型。
[0031] 假设旋转目标由K个散射中心构成，当宽带雷达带宽的距离分辨率大于目标尺寸 时，旋转目标上各个散射中心的雷达子回波都被看作是来自于同一个距离单元。未考虑遮 挡效应时，旋转目