一种基于子孔径划分的ISAR成像包络对齐方法技术

本发明专利技术一种基于子孔径划分的ISAR成像包络对齐方法，属于ISAR雷达成像技术方法。本发明专利技术的方法利用了运动的平稳连续性，将全孔径分成若干个子孔径，保证了子孔径内部散射系数的近似均匀性，抑制了包络误差的积累，同时每个子孔径单独进行包络误差估计，再利用多项式拟合出最优化的全局包络误差函数，克服了由于大角度高分辨成像时包络漂移误差累积的问题；本发明专利技术的方法区别于最大相关处理方法和整体最优准则的对齐方法，本发明专利技术的方法利用了观测目标连续平稳的运动特性，将全孔径划分为若干个子孔径，对子孔径和全孔径包络误差分别进行线性和高阶多项式建模，有效抑制由于积累转角引起的包络漂移，从而有效提升包络对齐精度。

A method of ISAR imaging envelope alignment based on subaperture partition

The invention relates to a ISAR imaging envelope alignment method based on sub aperture division, which belongs to the ISAR radar imaging technique. The method of the invention utilizes the smooth and steady movement, the full aperture is divided into several sub apertures, ensure the uniformity of the internal approximate scattering coefficient of sub aperture, suppress the envelope error accumulation, and each sub aperture separate envelope error estimation, then the polynomial fitting global envelope error function optimization, overcomes the cumulative envelope drift error angle high resolution imaging problems; alignment method difference between the method of the invention to the maximum processing methods and global optimization criterion, the method of the invention utilizes the motion characteristics of the target continuously and smoothly, the full aperture is divided into several sub aperture, sub aperture and full aperture envelope were linear error and high order polynomial modeling, effectively restrain the accumulation caused by envelope drift angle, so as to effectively enhance the envelope alignment Accuracy.

【技术实现步骤摘要】
一种基于子孔径划分的ISAR成像包络对齐方法
本专利技术一种基于子孔径划分的ISAR成像包络对齐方法，属于ISAR雷达成像技术方法。
技术介绍
逆合成孔径雷达(ISAR)的关键技术包括平动补偿和转台成像等，其中平动补偿又分为包络误差补偿和相位误差补偿2个步骤。包络误差补偿一般利用包络对齐技术实现。精确的包络对齐是相位误差补偿的基础，又是高分辨成像的前提条件。1980年代初期由Chen等提出包络对齐的互相关法以来，国内外学者又相继提出了多种ISAR包络对齐算法，大致可分为2类：第1类是最大相关处理方法，此类方法运用相邻距离像的相关性估计目标回波的包络徙动量，然后进行包络对齐处理，但其具有较强的噪声敏感性；第2类是整体最优准则处理方法，通常以平均距离像的熵最小或对比度最大为准则，通过优化算法求解各次回波间的包络偏移来进行包络对齐处理，此类方法利用了整体回波距离像的相关性信息，对噪声和散射起伏都具有较好的抑制作用，但一般需对全孔径包络误差进行搜索，相对于第1类方法运算量增大。实质上，两类包络对齐算法都是基于回波数据的相关性进行处理。但由于ISAR目标成像需要利用上百次回波(总转角一般要求大于等于5°)，随着视线角的变化，目标散射点相对雷达几何构型上会发生细微变化，从而引起包络漂移，产生较大的包络误差积累，这种情况随着分辨率的提升影响会更加剧烈，影响传统算法的包络对齐精度，从而影响最终的ISAR成像。由于ISAR目标成像需要利用上百次回波(总转角一般要求大于5°)，随着视线角的变化，目标散射点相对雷达几何构型上会发生细微变化，从而引起包络漂移，产生较大的包络误差积累，这种情况随着分辨率的提升影响会更加剧烈，而传统方法的包络误差估计是全孔径相关方法，无法修正由于大积累转角引起的散射点散射系数起伏带来的全孔径回波包络相关性变差的问题，致使包络对齐误差累积增大，严重影响ISAR的方位向成像，造成方位成像的模糊。
技术实现思路
本专利技术的技术解决问题是：克服现有技术的不足，提出一种基于子孔径划分的ISAR成像包络对齐方法。本专利技术的技术解决方案是：一种基于子孔径划分的ISAR成像包络对齐方法，该方法的步骤包括：(1)将ISAR成像周期所需的全孔径时间分为M个子孔径，每个子孔径长度相同；然后采用最小熵准则获取每个子孔径一阶包络误差导数向量ΔR′m，m＝1,2...M；其中，ΔR′m为第m个子孔径的一阶包络误差导数向量；则全孔径时间内的一阶包络误差导数向量为ξ：ξ＝[ΔR′1,ΔR′2,...ΔR′m,...ΔR′M-1,ΔR′M]T(1)(2)用每个子孔径的中心时刻tm构造时间的序列α＝[t1,t2,...tm,...tM]T，然后采用高阶多项式近似拟合包络误差，构造每个子孔径的包络误差函数为ΔR(tm)：其中，t1为第一个子孔径的中心时刻，t2为第二个子孔径的中心时刻，tm为第m个子孔径的中心时刻，tM为第M个子孔径的中心时刻，aq为多项式系数，Q为包络误差拟合的阶数，q为子孔径的中心时刻tm的阶数，q＝1,2,...Q；对式(2)求导得到每个子孔径的一阶包络误差导数向量为ΔR′(tm)：(3)将步骤(2)得到的公式(3)进行展开，得到根据步骤(1)得到的公式(1)可知，令令多项式系数得到(4)对于ξ＝φβ形式的超定矩阵，利用最小二乘法对包络误差函数中多项式系数β进行估计，得到：(5)将步骤(4)得到的多项式系数β的估计结果带入到步骤(2)得到的公式(2)中，得到每个子孔径的包络误差函数的估计结果再利用公式(2)得到M个子孔径对应的包络误差值ΔR(tm)；(6)根据每个子孔径的包络误差函数值ΔR(tm)，将每个子孔径内的每一个回波信号的所有距离单元向距离零点方向平移ΔR(tm)个距离单元，完成包络对齐。ISAR成像积累转角大于等于5°。所述的ISAR成像一个周期所需的时间为全孔径时间。子孔径数M的确定方法为：首先使得每个子孔径的包络误差小于半个距离单元，然后设定M的初值M0，根据公式(2)获得每个子孔径的包络误差函数若对于每个子孔径都满足：式(8)中，c为光速，fs为采样率；则求得满足公式(8)的最小M0值即为最小的子孔径个数。包络误差拟合的阶数Q满足Q＝M-1。本专利技术与现有技术相比的优点在于：(1)本专利技术的方法利用了运动的平稳连续性，将全孔径分成若干个子孔径，保证了子孔径内部散射系数的近似均匀性，抑制了包络误差的积累，同时每个子孔径单独进行包络误差估计，再利用多项式拟合出最优化的全局包络误差函数，克服了由于大角度高分辨成像时包络漂移误差累积的问题；(2)本专利技术的方法区别于最大相关处理方法和整体最优准则的对齐方法，本专利技术的方法利用了观测目标连续平稳的运动特性，将全孔径划分为若干个子孔径，对子孔径和全孔径包络误差分别进行线性和高阶多项式建模，有效抑制由于积累转角引起的包络漂移，从而有效提升包络对齐精度；(3)由于传统包络对齐方法只着眼于回波数据的相关性，而没有利用观测目标的运动特性，而在实际ISAR成像中，一定的观测时间内，目标特别是空间卫星目标的平动一般较为平稳，回波信号的包络误差调制也是平稳连续变化的，而且包络对齐的精度仅要求在距离分辨单元数量级(一般在米级或者分米级)，所以在包络对齐处理中，根据上述平稳运动特性，将包络误差建模为慢时间的高阶多项式形式，就能够有效克服低信噪比情况下包络误差估计的突跳问题；(4)本专利技术的方法，将包络对齐的全孔径分成若干个子孔径，每个子孔径单独进行包络误差估计，再利用多项式拟合出最优化的全局包络误差函数，克服了由于大角度高分辨成像时包络漂移误差累积的问题；(5)本专利技术的方法将包络对齐的全孔径分成若干个子孔径，每个子孔径单独进行包络误差估计，再利用多项式拟合出最优化的全局包络误差函数，从而实现包络的精确对齐。附图说明图1为相关最大法对齐方法得到的ISAR图像；图2为实施例中采用本申请的方法得到的ISAR图像。具体实施方式一种基于子孔径划分的ISAR成像包络对齐方法，该方法针对成像积累转角大于等于5°时的包络对齐问题，该方法的步骤包括：(1)定义一个ISAR成像周期所需的时间为全孔径时间，将全孔径时间分为M个子孔径，每个子孔径长度相同；用每个子孔径的中心时刻tm构造时间的序列α＝[t1,t2,...tm,...tM]T，t1为第一个子孔径的中心时刻，t2为第二个子孔径的中心时刻，tm为第m个子孔径的中心时刻，tM为第M个子孔径的中心时刻，m＝1,2...M；当M足够大时，即每个子孔径的方位向采样时间足够小，近似认为动目标在子孔径时间内的速度不变，所以可将包络对齐误差可看做为线性，简化了子孔径的包络误差模型；采用最小熵准则获取每个子孔径一阶包络误差导数向量ΔR′m，m＝1,2...M，ΔR′1为第一个子孔径的一阶包络误差导数向量，ΔR′2为第二个子孔径的一阶包络误差导数向量，ΔR′m为第m个子孔径的一阶包络误差导数向量；ΔR′M为第M个子孔径的一阶包络误差导数向量；则全孔径时间内的一阶包络误差导数向量为ξ：ξ＝[ΔR′1,ΔR′2,...ΔR′m,...ΔR′M-1,ΔR′M]T(1)其中，ΔR′m为第m个子孔径的一阶包络误差导数向量；(2)采用高阶多项式近似拟合包络误差，构造每个子孔径的本文档来自技高网...

【技术保护点】
一种基于子孔径划分的ISAR成像包络对齐方法，其特征在于该方法的步骤包括：(1)将ISAR成像周期所需的全孔径时间分为M个子孔径，每个子孔径长度相同；然后采用最小熵准则获取每个子孔径一阶包络误差导数向量ΔR′m，m＝1,2...M；其中，ΔR′m为第m个子孔径的一阶包络误差导数向量；则全孔径时间内的一阶包络误差导数向量为ξ：ξ＝[ΔR'1,ΔR'2,...ΔR'm,...ΔR'M‑1,ΔR'M]

【技术特征摘要】
1.一种基于子孔径划分的ISAR成像包络对齐方法，其特征在于该方法的步骤包括：(1)将ISAR成像周期所需的全孔径时间分为M个子孔径，每个子孔径长度相同；然后采用最小熵准则获取每个子孔径一阶包络误差导数向量ΔR′m，m＝1,2...M；其中，ΔR′m为第m个子孔径的一阶包络误差导数向量；则全孔径时间内的一阶包络误差导数向量为ξ：ξ＝[ΔR'1,ΔR'2,...ΔR'm,...ΔR'M-1,ΔR'M]T(1)(2)用每个子孔径的中心时刻tm构造时间的序列α＝[t1,t2,...tm,...tM]T，然后采用高阶多项式近似拟合包络误差，构造每个子孔径的包络误差函数为ΔR(tm)：其中，t1为第一个子孔径的中心时刻，t2为第二个子孔径的中心时刻，tm为第m个子孔径的中心时刻，tM为第M个子孔径的中心时刻，aq为多项式系数，Q为包络误差拟合的阶数，q为子孔径的中心时刻tm的阶数，q＝1,2,...Q；对式(2)求导得到每个子孔径的一阶包络误差导数向量为ΔR'(tm)：(3)将步骤(2)得到的公式(3)进行展开，得到根据步骤(1)得到的公式(1)可知，令令多项式系数得到

【专利技术属性】
技术研发人员：宋大伟，尚社，孙文锋，罗熹，李栋，温媛媛，王建晓，范晓彦，
申请(专利权)人：西安空间无线电技术研究所，
类型：发明
国别省市：陕西,61