基于三维重构的空间目标姿态外推方法技术

本发明专利技术公开了一种基于三维重构的空间目标姿态外推方法，主要解决现有技术在多任务雷达样本个数较少时目标识别率低的问题，其实现过程是：1)对数据进行均匀分帧，并进行ISAR成像；2)对每一幅ISAR图像进行方位定标；3)提取每一幅ISAR图像的散射中心参数集合；4)对ISAR图像集合的散射中心进行关联，得到轨迹矩阵；5)对轨迹矩阵进行奇异值分解，得到投影矩阵和位置矩阵；6)对投影矩阵和位置矩阵进行三维重构，得到散射中心矩阵；7)根据散射中心矩阵，对每一帧内的数据进行外推，得到每一帧外推后的时域实雷达距离像样本，用于后续目标识别。本发明专利技术采用数据外推的方法能显著提高目标识别率，可用于对雷达距离像的识别。

Attitude Extrapolation Method of Space Target Based on Three-dimensional Reconstruction

【技术实现步骤摘要】
基于三维重构的空间目标姿态外推方法
本专利技术属于数据处理
，特别涉及一种小角域内样本的数据外推方法，可用于目标识别。
技术介绍
当目标沿雷达视线方向RLOS的尺寸远大于雷达的距离分辨率时，就不能把目标当作一个“点”来处理。而由于目标占据多个距离分辨单元，相应的雷达回波中呈现出的起伏和尖峰，就反映着目标散射点沿距离方向上的分布，这些回波就被称为高分辨距离像HRRP。近年来，基于HRRP的雷达目标识别技术由于其特有的优势越来越引起各国的关注。然而，在实际情况中，一方面由于雷达的多任务工作模式，导致对目标的回波的非连续性采集；另一方面，由于空间目标的非合作性，空间目标的变轨会导致空间目标的姿态变化，进而影响统计模型的准确性，难以得到理想的识别结果。因此考虑如何在已知角域内少量样本的情况下得到较高的识别结果是亟待解决的问题。西安电子科技大学在其申请的专利技术专利“基于稀疏散射中心提取的距离像数据外推方法”(公开号：105068062B，申请号：201510512213.7)中公开了一种基于稀疏散射中心提取的距离像数据外推方法。该方法通过稀疏求解，得到单幅ISAR图像的散射中心的参数集合；根据目标散射中心参数组合，得到每一帧外推后的复HRRP样本；对每一帧数据进行逆快速傅里叶变换IFFT，得到每一帧外推后的时域实HRRP样本，用于后续的目标识别。该方法存在的不足之处是，只适合扁平目标或者在同一平面旋转的目标，如飞机等,对于具有三维旋转的空间目标，将不再适用。
技术实现思路
本专利技术的目的在于针对上述现有技术的不足，提出一种基于三维重构的空间目标距离像数据外推方法，以在三维旋转的目标的不同姿态下，增多样本个数，提高目标识别率。本专利技术的技术思路是在高分辨率雷达中，将目标看作是由一系列散射中心组成；通过对序列ISAR图像的散射中心提取、关联、重构，估计得到目标的三维散射中心模型，通过数据外推得到更大角域内的复HRRP样本，进而增多样本个数，提高模型参数估计的准确性，最终达到提高目标识别率的目的。其实现步骤包括如下：(1)对获取的复高分辨距离像HRRP进行均匀分段，每一段定义为一帧，并对每一帧进行ISAR成像，得到ISAR像集合；(2)对ISAR像集合里的每一帧ISAR像进行方位定标，得到图像的方位分辨率；(3)设定散射中心的幅度阈值，利用峰值搜索方法获取每一帧ISAR像的散射中心参数集合：其中，表示中参数组合的个数，表示第i个参数组合对应的复强度，表示第i个参数组合对应的位置信息；(4)根据目标散射中心参数集合对散射中心进行关联，得到散射中心轨迹矩阵W；(5)对轨迹矩阵W进行奇异值分解W＝UΣVT，得到投影矩阵位置矩阵其中表示左酉矩阵U的前3列，表示Σ的前3列和前3行，表示右酉矩阵V的前3列；(6)根据投影矩阵和位置矩阵进行散射中心的三维重构，得到散射中心矩阵P0：(7)根据目标散射中心矩阵P0，通过下式外推得到未观测的复HRRP样本其中，f表示雷达回波的频率，和分别表示雷达射线的方位角和俯仰角,exp(·)表示自然指数函数，(P0(1,i),P0(2,i),P0(3,i))表示矩阵P0的第i列，c为光速；(8)对每一段内的数据进行步骤(1)到步骤(7)操作，得到各段外推后的复HRRP样本，对各段复HRRP样本作逆快速傅里叶变换IFFT，并取模值，得到外推后的时域实HRRP样本，用于后续的目标识别。本专利技术与现有技术相比具有以下优点：1.解决了多任务雷达不能长时间获取HRRP样本的问题。现有技术中，多任务雷达间歇性采集HRRP样本导致可用复HRRP样本数较少的问题始终不能较好的解决，本专利技术提出通过目标的二维ISAR图像进行目标三维散射中心的重构，可以较准确的得到目标三维散射中心的参数组合，保证了外推角域后复HRRP样本的准确性，较好的解决了可用HRRP数据较少的问题。2.提高了样本不充足情况下的雷达目标识别率。现有的基于统计建模的HRRP识别方法均是以大量训练数据为前提的，只有在训练样本足够多的情况下，估计得到的统计模型才有可信度，而实际情况中，由于大多数空中目标的非合作性，得到的回波数据往往是有限的，本专利技术提出的基于少量复HRRP样本的三维散射中心估计，根据估计模型进行数据外推的方法可以有效提高识别率。附图说明图1是本专利技术的实现总流程图；图2是本专利技术中进行三维重构的子流程图；图3是用本专利技术外推后的数据和现有数据进行目标识别的识别率比较图。具体实施方式以下结合附图对本专利技术的实施力和效果、作进一步详细描述。参照图1，本专利技术的具体实现步骤如下：步骤1，对数据进行均匀分帧并进行ISAR成像。对雷达获取的复高分辨距离像HRRP进行均匀分帧，将长度为α的连续复HRRP样本定义为一帧数据，其中，B表示雷达带宽，L表示目标横向尺寸，c表示光速，并对每一帧进行ISAR成像，得到ISAR像集合。步骤2，对每一幅ISAR图像进行方位定标。根据雷达观测的跟踪数据得到每一帧数据雷达射线相对于目标的转角β，根据下式求得方位分辨率：式中，dA表示方位分辨率，λ表示雷达发射波长，β表示每一帧数据雷达射线相对于目标的转角。步骤3，提取每一幅ISAR图像的散射中心参数集合。现有提取ISAR图像散射中心参数集合的方法有峰值搜索法和稀疏求解法，本步骤使用但不限于峰值搜索法，其实现步骤如下：(3a)依据经验设定散射中心的幅度阈值TA；(3b)提取每一幅ISAR图像中幅度值高于幅度阈值TA的像素位置(Xp,Yp)，以及该像素的幅度值其中，(Xp,Yp)表示第p个像素的坐标，p＝1…P，P表示像素的总个数；(3c)根据像素位置(Xp,Yp)，按照下式求解散射中心位置(3d)根据散射中心位置和幅度得到散射中心参数集合其中，表示中参数组合的个数。步骤4，对ISAR图像集合的散射中心进行关联，得到轨迹矩阵W。(4a)对散射中心参数集合通过下式进行关联：式中，表示第n帧图像的第i个参数组合对应的位置信息，表示第n帧图像的第i个参数组合对应的复强度，argmin表示求最小值操作，||·||2表示2范数操作；(4b)根据关联之后的散射中心位置获得的轨迹矩阵W如下：式中，表示第j帧图像的第i个参数组合对应的位置信息，j＝1,2,…,n，表示中参数组合的个数，n表示每一组数据中图像帧数。步骤5，对轨迹矩阵W进行奇异值分解。(5a)将轨迹矩阵W分解为3项，即W＝UΣVT，得到左酉矩阵U和右酉矩阵V，以及对角矩阵Σ；(5b)对这两个酉矩阵U，V和对角矩阵Σ按照下式求解，得到投影矩阵位置矩阵其中表示左酉矩阵U的前3列，表示对角矩阵Σ的前3列和前3行，表示右酉矩阵V的前3列。步骤6，根据投影矩阵和位置矩阵进行散射中心的三维重构，得到散射中心矩阵P0：现有三维重构的方法有矩阵分解法和速度场重构法，本步骤使用但不限于矩阵分解法，参照图2，本步骤的具体实现如下：(6a)按照下式求解内积向量G：其中ii(i＝1,2,…,n)表示的第i行，ji(i＝1,2,…,n)表示的第n+i行，gT是一种定义的计算法则，定义如下：gT(ii,ji)＝[ii(1)·ji(1)ii(1)·ji(2)+ii(2)·ji(1)ii(1)·ji(3)+ii(3)·ji(1)ii(2)·ji(2)ii(2)·j本文档来自技高网...

【技术保护点】
1.一种基于三维重构的空间目标姿态外推方法，包括如下：(1)对获取的复高分辨距离像HRRP进行均匀分段，每一段定义为一帧，并对每一帧进行ISAR成像，得到ISAR像集合；(2)对ISAR像集合里的每一帧ISAR像进行方位定标，得到图像的方位分辨率；(3)设定散射中心的幅度阈值，利用峰值搜索方法获取每一帧ISAR像的散射中心参数集合：

【技术特征摘要】
1.一种基于三维重构的空间目标姿态外推方法，包括如下：(1)对获取的复高分辨距离像HRRP进行均匀分段，每一段定义为一帧，并对每一帧进行ISAR成像，得到ISAR像集合；(2)对ISAR像集合里的每一帧ISAR像进行方位定标，得到图像的方位分辨率；(3)设定散射中心的幅度阈值，利用峰值搜索方法获取每一帧ISAR像的散射中心参数集合：其中，表示中参数组合的个数，表示第i个参数组合对应的复强度，表示第i个参数组合对应的位置信息；(4)根据目标散射中心参数集合对散射中心进行关联，得到散射中心轨迹矩阵W；(5)对轨迹矩阵W进行奇异值分解W＝UΣVT，得到投影矩阵位置矩阵其中表示左酉矩阵U的前3列，表示对角矩阵Σ的前3列和前3行，表示右酉矩阵V的前3列；(6)根据投影矩阵和位置矩阵进行散射中心的三维重构，得到散射中心矩阵P0：(7)根据目标散射中心矩阵P0，通过下式外推得到未观测的复HRRP样本其中，f表示雷达回波的频率，和分别表示雷达射线的方位角和俯仰角,exp(·)表示自然指数函数，(P0(1,i),P0(2,i),P0(3,i))表示矩阵P0的第i列，c为光速；(8)对每一段内的数据进行步骤(1)到步骤(7)操作，得到各段外推后的复HRRP样本，对各段复HRRP样本作逆快速傅里叶变换IFFT，并取模值，得到外推后的时域实HRRP样本，用于后续的目标识别。2.根据权利要求1所述的方法，其特征在于，(1)中的均匀分段，每一段长度为其中，B表示雷达带宽，L表示目标横向尺寸，c表示光速。3.根据权利要求1所述的方法，其特征在于，(2)中对ISAR像集合里的每一帧ISAR像进行方位定标，...

【专利技术属性】
技术研发人员：王家东，杜兰，杨栋文，张晓乐，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61