一种能够去除突跳和漂移误差的改进ISAR距离对准方法技术

本发明专利技术公开了一种能够去除突跳和漂移误差的改进ISAR距离对准方法，主要用于解决全局最小熵距离对准算法的突跳误差和成块距离漂移问题，包括顺序执行的以下步骤：步骤一、采用全局最小熵算法对ISAR脉压数据进行距离对准；步骤二、计算距离对准后每条回波与相邻一维距离像的归一化相关系数；步骤三、通过分析归一化相关系数判断回波数据是否为突跳脉冲；步骤四、通过突跳脉冲与相邻脉冲互相关将成块距离漂移的脉冲对齐；步骤五、对步骤四中获得的对准数据进行图像处理，得到距离对准图像的边缘轨迹；步骤六、利用Radon变换计算步骤五两条边缘轨迹的倾斜角；步骤七、根据步骤六得到的倾斜角对步骤四对准数据进行倾斜校正，获得ISAR高精度距离对准结果。

【技术实现步骤摘要】

本专利技术属于雷达成像技术，特别是一种ISAR成像技术。
技术介绍
逆合成孔径雷达(Inverse-Synthetic-Aperture-Radar，简称ISAR)距离对准是消除目标相对于雷达平动造成的相邻回波在距离向上的错位，对准后，目标上同一散射点的信号在不同的回波脉冲中都位于同一距离单元中。典型的算法有：相邻包络互相关法、时域和频域积累互相关法、包络最大修正峰度法、全局对准算法和全局最小熵算法等。全局最小熵算法是以平均距离像(Averagerangeprofile，简称ARP)熵最小为原则作为对准依据，通过迭代的方法使其对准效果更优。这类算法对脉冲样本之间的相关性要求不高，同时可以较好地抑制突跳和漂移误差，尤其当遇到回波有间断或回波突跳的情况时该算法具有较好的鲁棒性。由于孔径时间内目标相对雷达的转动，回波相关性有较为复杂的变化，并且相关性总体呈下降趋势，该变化对全局最小熵算法对准准确性有一定影响，会造成脉冲突跳以及成块脉冲距离漂移的问题，对准的漂移误差会造成ISAR经过相位梯度自聚焦算法(PhaseGradientAutofocus，简称PGA)处理后成像出现散焦和重影现象。
技术实现思路
专利技术目的：为了克服现有技术的不足，本专利技术提供一种能够去除突跳和漂移误差的改进ISAR距离对准方法，用于解决现有距离对准算法在目标运动复杂情况下的突跳和成块漂移误差造成的成像散焦和重影问题。技术方案：一种能够去除突跳和漂移误差的改进ISAR距离对准方法，包括顺序执行的以下步骤：步骤一：采用全局最小熵算法对ISAR脉压数据进行距离对准获得最小熵对准数据；步骤二：计算距离对准后每条回波数据与相邻一维距离像的归一化相关系数；步骤三：通过分析归一化相关系数判断回波数据是否为突跳脉冲；步骤四：通过突跳脉冲与相邻脉冲互相关将成块距离漂移的脉冲对齐；步骤五：对步骤四中获得的对准数据进行图像处理，得到距离对准图像的两条边缘轨迹；步骤六：利用Radon变换计算步骤五中两条边缘轨迹的倾斜角；步骤七：根据步骤六得到的倾斜角对步骤四中获得的对准数据进行倾斜校正，获得改进ISAR高精度距离对准结果。进一步，本专利技术中解决距离对准突跳和成块漂移误差的过程：所述步骤三中：先计算归一化相关系数的平均值，将每条回波数据的归一化相关系数与上述平均值比较，若归一化相关系数小于一定范围，判定为突跳脉冲；所述步骤四中：通常在最小熵距离对准算法中为实现亚单元对准精度，采用频域补零实现时域插值。为不破坏对准精度，首先根据对准数据频谱特性补零，补零倍数一般为8倍，然后返回到距离时域计算突跳脉冲与相邻脉冲互相关函数最大时对应的偏移量，最后将偏移量除以相应补零倍数构造频移因子，实现该突跳脉冲与下一个突跳脉冲间所有漂移脉冲的对准。所述步骤五中：利用canny算子提取去除成块漂移对准图像边缘，再利用形态学图像闭运算处理填充边缘孔洞，再提取边缘得到两条倾斜的边缘轨迹；所述步骤六中：利用Radon变换计算两条边缘轨迹的倾斜角，取两倾斜角平均值，即为去除成块距离漂移的距离对准图像倾斜角度；所述步骤七中：利用步骤六所得平均倾斜角度计算倾斜斜率，对步骤四所得去除距离漂移的对准数据进行倾斜校正。有益效果：本专利技术将归一化相关系数作为脉冲相关性能评价标准检测突跳脉冲，对于降低成像质量的突跳脉冲数据，不是采用剔除的方式，而是利用与其前一个相邻脉冲互相关计算偏移量将成块漂移脉冲对准，解决了传统ISAR距离对准突跳和漂移误差；考虑到成块漂移脉冲对准后引入的距离对准图像倾斜问题，首次提出利用图像处理方法对距离对准结果做倾斜校正；经过成块漂移脉冲对准和倾斜校正后的距离对准数据利用PGA相位补偿后的成像结果聚焦效果有了很大改善，消除了重影。附图说明图1为本专利技术流程图；图2为Radon变换原理说明图；图3为传统距离对准算法处理结果；图4为归一化相关系数值；图5为成块距离漂移对准后结果；图6是图像处理结果；图6(a)为canny算子提取边缘图像，图6(b)是闭运算后二值图，图6(c)是提取的两条边缘轨迹；图7是Radon变换检测结果图；图7(a)为上边缘轨迹Radon变换结果图，图7(b)为下边缘轨迹Radon变换结果图；图8是倾斜校正后对准结果；图9是经过PGA相位补偿后成像结果图；图9(a)为传统全局最小熵对准后成像结果图，图9(b)为成块漂移对准后成像结果图，图9(c)为倾斜校正后成像结果图。具体实施方式下面结合附图对本专利技术做更进一步的解释。如图1所示为本专利技术的基于归一化相关系数和图像处理的ISAR距离对准改进方法的流程图。包含的处理步骤可以归为两个阶段。第一价段为预处理，将ISAR脉压数据经过传统全局最小熵算法进行距离对准。第二阶段可以分为两个部分，第一部分为去除突跳和成块距离漂移，第二部分为基于图像处理的倾斜校正。首先对去除突跳和成块距离漂移问题进行说明。在最小熵距离对准后，利用(1)式计算每条回波与其相邻的脉冲回波即相邻一维距离像的归一化相关系数：cor(n)=snTsn-1Σm=1Msn2(m)Σm=1Msn-12(m),1≤n≤N---(1)]]>其中sn和sn-1为两个相邻脉冲回波，N为脉冲数，M为距离单元数，m为距离单元序号，n为回波序号，T代表转置。在已有归一化相关系数的基础上判断是否为突跳脉冲，先利用(2)式计算归一化相关系数的平均值：cor_ave=1NΣn=1Ncor(n)---(2)]]>若cor(n)<cor_ave-Δ，则sn为突跳脉冲，Δ视具体数据取值，数据质量好即脉冲之间相关性很强(相关系数的平均值很接近1)的情况下δ值小一点，可取0.02-0.04，否则Δ应大点，取0.1左右。将最小熵对准数据傅里叶变换至距离向频域，根据频谱对称性补零后逆傅里叶变换至距离时域。将内插处理后突跳脉冲与前一个相邻脉冲互相关，求互相关函数取最大值对应的距离偏移量，互相关函数为下列(3)式：R(s)=∫Si*(r)Si-1(r+s)dr[∫|Si(r)|2dr·∫|Si-1(r)|2dr·]1/2---(3)]]>其中Si(r)和Si-1(r)为突跳脉冲对应的一维距离像及其相邻一维距离像，r为距离偏移量，i表示脉冲序号，表示Si(r)取共轭。用每个突跳脉冲计算得到的偏移量除以相应补零倍数构造频移因子，利用离散傅里叶变换性质频移特性通过频域乘以线性相位因子实现亚距离单元的时域移位，将该突跳脉冲与下一个突跳脉冲之间所有成块距离漂移脉冲对准。接下来对基于图像处理的倾斜校正进行说明。在去除突跳和成块距离漂移之后，对准结果因为偏移量积累会出现距离倾斜，如果不进行倾斜校正将会降低成像质量。倾斜校正需要得到距离对准边缘轨迹，再利用Radon变换检测倾斜角。首先通过canny算子检测目标距离对准边缘。canny算子是一种性能优良的边缘检测算子，由下列基本步骤组成：用高斯滤波器平滑距离对准图像；用一阶偏导的有限差分来计算梯度值和方向角；对梯度幅值进行非极大值抑制；用双阈值处理和连接分析来检测并连接边缘。然后利用形态学闭运算填充边本文档来自技高网...

【技术保护点】
一种能够去除突跳和漂移误差的改进ISAR距离对准方法，其特征在于：包括顺序执行的以下步骤：步骤一、采用全局最小熵算法对ISAR脉压数据进行距离对准获得最小熵对准数据；步骤二、计算距离对准后每条回波数据与相邻一维距离像的归一化相关系数；步骤三、通过分析归一化相关系数判断回波数据是否为突跳脉冲；步骤四、通过突跳脉冲与相邻脉冲互相关将成块距离漂移的脉冲对齐；步骤五、对步骤四中获得的对准数据进行图像处理，得到距离对准图像的两条边缘轨迹；步骤六、利用Radon变换计算步骤五中两条边缘轨迹的倾斜角；步骤七、根据步骤六得到的倾斜角对步骤四中获得的对准数据进行倾斜校正，获得改进ISAR距离对准结果。

【技术特征摘要】
1.一种能够去除突跳和漂移误差的改进ISAR距离对准方法，其特征在于：包括顺序执行的以下步骤：步骤一、采用全局最小熵算法对ISAR脉压数据进行距离对准获得最小熵对准数据；步骤二、计算距离对准后每条回波数据与相邻一维距离像的归一化相关系数；步骤三、通过分析归一化相关系数判断回波数据是否为突跳脉冲；步骤四、通过突跳脉冲与相邻脉冲互相关将成块距离漂移的脉冲对齐；步骤五、对步骤四中获得的对准数据进行图像处理，得到距离对准图像的两条边缘轨迹；步骤六、利用Radon变换计算步骤五中两条边缘轨迹的倾斜角；步骤七、根据步骤六得到的倾斜角对步骤四中获得的对准数据进行倾斜校正，获得改进ISAR距离对准结果。2.根据权利要求1所述的一种能够去除突跳和漂移误差的改进ISAR距离对准方法，其特征在于：所述步骤三中，先计算归一化相关系数的平均值，将每条回波数据的归一化相关系数与上述平均值比较，若归一化相关系数与归一化相...

【专利技术属性】
技术研发人员：朱岱寅，左潇丽，聊蕾，许至威，杨鸣东，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏;32