一种基于压缩感知的ISAR成像脉冲估计算法制造技术

本发明专利技术属于雷达成像技术领域，具体提供了一种基于压缩感知的ISAR成像脉冲估计算法，其实现步骤为：（1）获得部分ISAR原始回波数据；（2）获得各脉冲一维距离像；（3）对各距离像进行PCA融合；（4）对融合结果进行Renyi熵阈值分割；（5）通过CS精确重构定理对分割结果进行脉冲估。本发明专利技术通过图像融合的方法对ISAR距离像进行信号融合及阈值分割，在简单的距离像预处理中，估计出压缩感知超分辨成像为达到最合适的成像效果所需的ISAR原始回波的最小数据量。

【技术实现步骤摘要】
一种基于压缩感知的ISAR成像脉冲估计算法
本专利技术属于雷达成像
，更进一步涉及ISAR图像信号通过压缩感知(CS)技术进行成像所需数据量估计方法，具体为一种基于压缩感知的ISAR成像脉冲估计算法。本专利技术可用于估计ISAR成像所需最小脉冲数目,而达到通过使用最少的数据量对ISAR信号进行压缩感知超分辨成像的目的。
技术介绍
通过对线性系统分析，可以描述雷达回波的性质：输入为雷达的发射脉冲，在过对目标的作用后，得到的输出即为雷达回波。雷达电波对目标(例如飞机等)作用，其一些部件会产生后向散射，一些平面部分也会后向镜面散射，这些构成了雷达回波的主要部分。所以，我们可以用散射点来表示目标，则可以用散射点模型来近似描述目标的分布函数，则其中具有相对较强后向散射的位置即为这些散射点。所以，可以看出，若目标在运动过程中，相对于雷达发生了转动或偏离，那写镜面反射的回波将不能被雷达接受到。由上所述，目标的散射点模型不仅与雷达长有关，同时随视角的变化而变化，经过分析和实验证明，若满足视角变化约为10°，则目标上散射点的位置及强度基本不会发生变化，而对于ISAR目标来说，目标一般需转动3°左右，所以能够考虑该模型。现代军用雷达一般作用的目标均为非机动的，还需对其进行识别跟踪。在采用了宽频带雷达信号后，使距离分辨率有了大幅度的提升，并且宽频带为雷达目标识别提供了很好的基础，较常规窄频带雷达在波形上包含了更加丰富的目标信息。若提高雷达频宽到几百兆赫，则回波能够达到高距离分辨率(HRR)为亚米级信号，且一般来说，HRR回波信号均表现为一维距离像。在考虑ISAR信号稀疏性后，可以对图像进行低数据量超分辨重建，目前已有许多学者验证了基于CS方法对ISAR信号进行成像的可能性以及相较于传统成像的优越性，但同时新的问题有出现了，到底该用多少数据即可重构出ISAR图像，同时保证较好的成像质量？ISAR原始数据需要进行距离向与方位向的二维处理，才能进行成像，针对实际效率要求，分析ISAR数据特性，希望在通过简单的距离向分析后即可给出估计最优成像数据量的方法，以在保证成像质量的同时最大化的减轻雷达系统压力。
技术实现思路
本专利技术的目的在于在保证成像效果及效率的情况下，估计出ISAR稀疏信号的压缩感知超分辨成像所需的最小脉冲数，减少雷达系统及硬件系统的压力，提出了一种基于压缩感知的ISAR成像脉冲估计算法，来估计压缩感知ISAR成像所需最小脉冲数的方法。通过图像融合的方法对ISAR距离像进行信号融合及阈值分割，在简单的距离像预处理中，估计出压缩感知超分辨成像为达到最合适的成像效果所需的ISAR原始回波的最小数据量。为实现上述目的，本专利技术技术方案是：提供一种基于压缩感知的ISAR成像脉冲估计算法，包括对ISAR信号回收系统获得的少量脉冲回波进行信号融合、对融合结果进行信号杂波分割及最终脉冲数估计三个过程，具体实现步骤如下：第一个过程，ISAR距离像PCA图像融合的具体步骤如下：(1)得到ISAR部分原始回波数据的一维距离像；(2)计算得到各距离像协方差；(3)根据所得协方差求得各距离像特征向量，并确定第一主成分所对应的特征向量；(4)据特征向量确定各距离像所应该分配的权重；(5)对各距离像进行加权融合得到融合图像；第二个过程，利用所获得的融合后的距离像进行信号及杂波分割，步骤为：(6)Renyi熵阈值分割；第三个过程，利用所获得的信号及杂波分割结果来估计压缩感知ISAR成像所需最小脉冲数的具体步骤如下：(7)CS精确重构理论：对于一个P稀疏的向量S＝{S1,S2,…,SQ}，Q表示向量长度，若要精确重构S，则观测值M需满足：M≥P*logQ据信号杂波区间计算信号区间内信号点数，通过CS精确重构理论计算最小观测值M，即为成像所需最小脉冲数；上述第一过程ISAR距离像PCA图像融合中步骤(1)中所述的得到ISAR部分原始回波数据的一维距离像，即为通过传统距离-多普勒成像方法，及脉冲压缩方法对部分原始数据进行距离像处理，获得ISAR数据的一维距离像图像大小为N，J为距离像个数；步骤(2)、(3)中所述计算得到各距离像协方差，而后根据所得协方差求得各距离像特征向量，并确定第一主成分所对应的特征向量，其具体方法如下：将步骤(1)中所得的距离像X＝{x1,x2,…,xN}T，求出其协方差矩阵Σ=Σj=1J(Xj-μ)(Xj-μ)T]]>其中，μ为均值向量；求出协方差矩阵Σ的全部N个特征值λ1,λ2,…λN，以及对应的特征向量u1,u2,…uN；其中λ1≥λ2≥…≥λN，第一主成分对应的特征向量即为u1；对于步骤(4)中的融合权重定义如下：假设有2幅图像A、B进行融合，其第一主成分对应的投影系数为xA,xB，则图A的权重分配为xA/(xA+xB)，图B的权重分配为xB/(xA+xB)；上述第二过程中对所获得的融合后的距离像进行信号及杂波分割中的步骤(6)中的Renyi熵阈值分割，按如下步骤进行：(6a)所获得的融合后的距离像Y＝{y1,y2,…,yN}，其中yi为图像中第i点的幅度值，而yi∈G＝{0,1,…,K-1}，K表示幅度最大值，记y(g)为图像中幅度为g的像素点个数，故幅度值g出现的概率h(g)为：h(g)=y(g)N,g=0,1,...,K-1]]>(6b)假设通过阈值k将图像划分成两类——信号和杂波背景，分别为C0和C1，则它们的先验概率如下：P0(k)=Σkh(g)]]>P1(k)=ΣKh(g)]]>则给出图像Y的δ(δ＞0)阶Renyi熵定义如下：Hδ0(k)=11-δlnΣg=0k[h(g)P0(k)]δ]]>Hδ1(k)=11-δlnΣg=k+1K[h(g)P1(k)]δ]]>(6c)则一维Renyi熵阈值方法得到的最优阈值k*为：k*=argmax0kK-1[Hδ0(k)+Hδ1(k)]]]>则所区分的信号与杂波区间二值图像定义为：其中0表示信号区间，K-1表示杂波区间。本专利技术的优点是：本专利技术就问题到底该用多少数据即可重构出ISAR图像，同时保证较好的成像质量，分析ISAR数据特性，给出估计最优成像数据量的方法，通过对原始雷达数据少量脉冲的距离像进行图像融合，并通过合适的信号和背景的阈值分割，压缩感知精确重构理论，估算出最少成像脉冲数，高效率的在保证成像质量的同时最大化的减轻雷达系统压力。以下将结合附图对本专利技术做进一步详细说明。附图说明图1是本专利技术的流程图；图2是本专利技术部分ISAR一维距离像；图3是本专利技术信号与噪声分割示例图；图4是本专利技术前8次脉冲的PCA融合结果图；图5是本专利技术方法ISAR数据距离像及估计所得脉冲数的CS成像结果。具体实施方式本专利技术为一种基于压缩感知的ISAR成像脉冲估计算法，包括对ISAR信号回收系统获得的少量脉冲回波进行信号融合、对融合结果进行信号杂波分割及最终脉冲数估计三个过程，具体实现步骤如下：第一个过程，ISAR距离像PCA图像融合的具体步骤如下：(1)通过本文档来自技高网...

【技术保护点】
一种基于压缩感知的ISAR成像脉冲估计算法，其特征在于：包括对ISAR信号回收系统获得的少量脉冲回波进行信号融合、对融合结果进行信号杂波分割及最终脉冲数估计三个过程，具体实现步骤如下： 第一个过程，ISAR距离像PCA图像融合的具体步骤如下： （1）得到ISAR部分原始回波数据的一维距离像； （2）计算得到各距离像协方差； （3）根据所得协方差求得各距离像特征向量，并确定第一主成分所对应的特征向量； （4）据特征向量确定各距离像所应该分配的权重； （5）对各距离像进行加权融合的到融合图像； 第二个过程，利用所获得的融合后的距离像进行信号及杂波分割，步骤为： （6）Renyi熵阈值分割； 第三个过程，利用所获得的信号及杂波分割结果来估计压缩感知ISAR成像所需最小脉冲数的具体步骤如下： （7）CS精确重构理论：对于一个P稀疏的向量S＝{S1,,S1,,…,SQ}，Q表示向量长度，若要精确重构S，则观测值M需满足： M≥P*logQ 据信号杂波区间计算信号区间内信号点数，通过CS精确重构理论计算最小观测值M，即为成像所需最小脉冲数。

【技术特征摘要】
1.一种基于压缩感知的ISAR成像脉冲估计算法，其特征在于：包括对ISAR信号回收系统获得的少量脉冲回波进行信号融合、对融合结果进行信号杂波分割及最终脉冲数估计三个过程，具体实现步骤如下：第一个过程，ISAR距离像PCA图像融合的具体步骤如下：(1)得到ISAR部分原始回波数据的一维距离像；(2)计算得到各距离像协方差；(3)根据所得协方差求得各距离像特征向量，并确定第一主成分所对应的特征向量；(4)据特征向量确定各距离像所应该分配的权重；(5)对各距离像进行加权融合得到融合图像；第二个过程，利用所获得的融合后的距离像进行信号及杂波分割，步骤为：(6)Renyi熵阈值分割；第三个过程，利用所获得的信号及杂波分割结果来估计压缩感知ISAR成像所需最小脉冲数的具体步骤如下：(7)CS精确重构理论：对于一个P稀疏的向量S＝{S1，S2，…，SQ}，Q表示向量长度，若要精确重构S，则观测值M需满足：M≥P*logQ据信号杂波区间计算信号区间内信号点数，通过CS精确重构理论计算最小观测值M，即为成像所需最小脉冲数。2.如权利要求1所述的一种基于压缩感知的ISAR成像脉冲估计算法，其特征在于：所述第一过程ISAR距离像PCA图像融合中步骤(1)中所述的得到ISAR部分原始回波数据的一维距离像，即为通过传统距离-多普勒成像方法，及脉冲压缩方法对部分原始数据进行距离像处理，获得ISAR数据的一维距离像图像大小为N，J为距离像个数；步骤(2)、...

【专利技术属性】
技术研发人员：张小华，焦李成，白婷，马晶晶，马文萍，王爽，侯彪，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61