一种多站雷达ISAR图像配准方法技术

本发明专利技术涉及雷达技术领域，尤其涉及一种多站雷达的ISAR图像配准方法。对多站雷达ISAR像进行配准时，选取其中回波信噪比最大的雷达作为雷达A，并将其经过平动分量补偿后得到的一维像序列作为联合平动分量补偿中的参考模板，其他雷达分别作为雷达B进行配准操作，即可完成多幅ISAR像之间的配准。本发明专利技术中的配准方法，是基于单个回波的，即对雷达B的每个回波分别进行配准。这种配准方式使得不同回波之间的配准过程不会相互干扰，保证了算法良好的抗噪性能。本发明专利技术中的配准方法是在信号域进行，因而ISAR成像之前就已经完成了配准过程。采用上述的配准操作，能够稳健高效地消除或者抑制多种失配因素，包括平动分量差异、雷达系统之间的时频不同步等。

【技术实现步骤摘要】
一种多站雷达ISAR图像配准方法
本专利技术涉及雷达
，尤其涉及一种多站雷达的逆合成孔径(ISAR,Inversesyntheticapertureradar)图像配准方法，该方法能够应用于干涉ISAR三维成像中。
技术介绍
ISAR成像技术是雷达对非合作目标跟踪、观测的重要手段。目前，ISAR成像技术已广泛应用于宽带雷达系统，能够获取空间卫星、空中飞机、海面舰船等目标的电磁散射系数二维分布图像。在二维ISAR成像的基础上，学者们提出了基于多站雷达的干涉ISAR三维成像方法。该方法需要位置不同的两部或者多部雷达同时接收目标反射电磁波并进行ISAR成像，利用同一散射中心在不同雷达ISAR像中的相位之差，即干涉相位，解算得到散射中心的空间坐标。结合ISAR像的二维分辨能力和干涉相位中包含的散射点坐标信息，理论上干涉ISAR成像能够获得目标散射中心的三维分布。干涉ISAR成像在目标侦察、三维结构重构、空间姿态估计等方面有着很大潜力，是宽带成像雷达发展的重要方向之一。ISAR图像配准是干涉ISAR三维成像中的关键步骤，其目标是确保同一散射点在不同雷达ISAR图像中所处的位置精确相同。配准结果影响散射中心提取、关联、相位干涉，从而最终影响散射中心的坐标解算结果，因此该过程需要尽可能精确地完成。实现ISAR图像配准主要有两个思路：基于图像域的配准方法和基于信号域的配准方法。参考文献一(BiaoTian,etal.“ANovelimageregistrationmethodforInISARimagingsystem”.SPIESecurity&Defense2014,Amsterdam,Netherlands.)中给出了基于图像域的ISAR图像配准方法：文中利用两幅ISAR图像之间的相关系数评价配准程度，根据最大相关系数搜索得到最佳的失配量并对其补偿，进而实现配准。该方法计算量较大，且抗噪性较差。参考文献二(ZhangQ,YeoTS,DuG,etal:“EstimationofThree-DimensionalMotionParametersinInterferometricISARImaging”,IEEETransactionsonGeoscienceandRemoteSensing,2004,42,(2),pp.292–300.)中研究了一种基于信号域的配准方法，该方法利用额外的天线对，估计目标相对不同雷达天线的平动分量差异，然后补偿这一差异，进而实现配准。该方法既需要额外的天线对配置，又要求极高的系统同步精度，工程上实现难度大。参考文献三(WenzhenWu,PengjiangHu,ShiyouXuetal.“ImageRegistrationforInISARBasedonJointTranslationalMotionCompensation”,IETRadar,Sonar&Navigation,2017,11(10),pp.1597-1603.)中公开了一种基于联合平动分量补偿的ISAR图像配准方法，是另一种信号域配准方法。该方法对来自两部雷达的信号进行相干处理，在平动分量补偿的过程中实现配准，具有良好的配准精度和抗噪性能。然而其配准过程包含多次迭代，计算效率不高，不便于处理实际雷达系统中的大批量数据。为适应干涉ISAR成像技术的发展，亟需一种精度高、抗噪性能好、计算量小的ISAR图像配准方法。
技术实现思路
本专利技术要解决的技术问题在于提供一种多站雷达ISAR图像配准方法，能够准确、高效、稳定地实现不同雷达ISAR图像之间的配准。该方法在平动分量补偿的过程中，采用联合处理的方法实现不同雷达ISAR图像之间的配准，不需要雷达之间高精度的时钟同步，更加适用于实际工程实践。相比于现有的ISAR图像配准方法，该方法不包含迭代过程，能够快速高效地完成ISAR图像距离向、方位向的配准，计算量小，并且能够保证较高的配准精度。本专利技术采用的技术方案为，一种多站雷达ISAR图像配准方法，包括如下步骤：S1其中一部雷达发射宽带信号照射目标，位置不同的多部雷达以相同的工作参数同时接收目标回波，并进行脉冲压缩和径向速度补偿，得到一维距离像序列；S2计算并对比各雷达一维距离像的信噪比，其中信噪比最大的记为雷达A，其他分别记为雷达B、雷达C、雷达D、…；S3完成雷达A一维像序列的平动分量补偿，其中包络对齐采用基于最大相关准则的对齐方法，初相校正采用基于最小熵准则的相位自聚焦方法；S4雷达B一维距离像序列的联合包络对齐：采用最大相关准则，将雷达B的每个一维距离像，参考同时刻下包络对齐后的雷达A一维距离像实现包络对齐；S5雷达B一维距离像序列的快速联合相位补偿：计算同时刻下，雷达B一维距离像相对雷达A一维距离像的平均相位差，并补偿该相位差，至此完成雷达B一维像序列的平动补偿；S6雷达A、雷达B一维距离像序列分别进行方位压缩，得到配准后的两幅ISAR像；S7将雷达C、雷达D、…等其他雷达的一维像序列分别等效为雷达B一维像序列，进行S4、S5、S6的操作，即可完成多站雷达ISAR图像的配准。上述的配准过程中，雷达A采用现有的距离对齐和相位补偿方法实现平动分量补偿。由于平动补偿方法的数据驱动性，补偿结果中不可避免地存在剩余误差。同理，雷达B、雷达C、雷达D、…等的一维像序列的平动补偿结果中也存在剩余误差。雷达A、雷达B、雷达C、雷达D、…等的一维像序列中不同的剩余误差，造成多站雷达ISAR像之间的失配。本专利技术中的联合包络对齐和快速联合相位补偿方法以多站雷达中回波信噪比最大的那部雷达(雷达A)的一维距离像为参考模板，完成雷达B、雷达C、雷达D、…等多站雷达一维距离像的平动分量补偿，使得两者的剩余平动分量误差相同。本专利技术中的配准方法将雷达A、雷达B、雷达C、雷达D、…等多站雷达的回波进行一一配准，是一种回波到回波的信号域配准方法，最终能够实现良好稳定的ISAR像配准效果。雷达B一维距离像序列的平动分量补偿包含两步——联合包络对齐和快速联合相位补偿，分别对应于上述配准过程的S4、S5。以雷达A、雷达B在同时刻接收到的一对回波为例，说明联合包络对齐和快速联合相位补偿的具体步骤。S4中，联合包络对齐包含以下步骤：S4.1记经过S3进行平动分量补偿后的雷达A一维距离像为：其中σi、RAi分别表示第i个散射点的后向散射系数、第i个散射点相对雷达A的距离，p(·)为归一化的回波包络，fc为载频，c为光速；在雷达A收到sA_comp(t)信号的同一时刻，雷达B收到目标回波，记为sB(t)，经过脉冲压缩，sB(t)可表达为：其中RBi表示第i个散射点到雷达B的距离；S4.2以sA_comp(t)为参考模板，根据最大相关准则，求取sB(t)相对sA_comp(t)的最佳偏移量；sA_comp(t)和sB(t)之间的相关函数可记为：sCorr(Δt)＝|sA_comp(t)·sB(t+Δt)|其中Δt表示时间偏移量；由于sA_comp(t)、sB(t)是同一目标在同一时刻的反射回波，它们之间的差异项为其中RAi-RBi为回波路程差，τm为接收采样窗口时间差；相关函数sCorr(Δt)存在唯一的峰值，记峰值位置为Δtmax；S4.3通过本文档来自技高网...

【技术保护点】
1.一种多站雷达ISAR图像配准方法，其特征在于，该方法包括以下步骤：S1其中一部雷达发射宽带信号照射目标，位置不同的多部雷达以相同的工作参数同时接收目标回波，并进行脉冲压缩和径向速度补偿，得到一维距离像序列；S2计算并对比各雷达一维距离像的信噪比，其中信噪比最大的记为雷达A，其他分别记为雷达B、雷达C、雷达D、…；S3完成雷达A一维像序列的平动分量补偿，其中包络对齐采用基于最大相关准则的对齐方法，初相校正采用基于最小熵准则的相位自聚焦方法；S4雷达B一维距离像序列的联合包络对齐：采用最大相关准则，将雷达B的每个一维距离像，参考同时刻下包络对齐后的雷达A一维距离像实现包络对齐；S5雷达B一维距离像序列的快速联合相位补偿：计算同时刻下，雷达B一维距离像相对雷达A一维距离像的平均相位差，并补偿该相位差，至此完成雷达B一维像序列的平动补偿；S6雷达A、雷达B一维距离像序列分别进行方位压缩，得到配准后的两幅ISAR像；S7将雷达C、雷达D、…等其他雷达的一维像序列分别等效为雷达B一维像序列，进行S4、S5、S6的操作，即可完成多站雷达ISAR图像的配准。

【技术特征摘要】
1.一种多站雷达ISAR图像配准方法，其特征在于，该方法包括以下步骤：S1其中一部雷达发射宽带信号照射目标，位置不同的多部雷达以相同的工作参数同时接收目标回波，并进行脉冲压缩和径向速度补偿，得到一维距离像序列；S2计算并对比各雷达一维距离像的信噪比，其中信噪比最大的记为雷达A，其他分别记为雷达B、雷达C、雷达D、…；S3完成雷达A一维像序列的平动分量补偿，其中包络对齐采用基于最大相关准则的对齐方法，初相校正采用基于最小熵准则的相位自聚焦方法；S4雷达B一维距离像序列的联合包络对齐：采用最大相关准则，将雷达B的每个一维距离像，参考同时刻下包络对齐后的雷达A一维距离像实现包络对齐；S5雷达B一维距离像序列的快速联合相位补偿：计算同时刻下，雷达B一维距离像相对雷达A一维距离像的平均相位差，并补偿该相位差，至此完成雷达B一维像序列的平动补偿；S6雷达A、雷达B一维距离像序列分别进行方位压缩，得到配准后的两幅ISAR像；S7将雷达C、雷达D、…等其他雷达的一维像序列分别等效为雷达B一维像序列，进行S4、S5、S6的操作，即可完成多站雷达ISAR图像的配准。2.一种根据权利要求1所述多站雷达ISAR图像配准方法，其特征在于：S4中，联合包络对齐包含以下步骤：S4.1记经过S3进行平动分量补偿后的雷达A一维距离像为：其中σi、RAi分别表示第i个散射点的后向散射系数、第i个散射点相对雷达A的距离，p(·)为归一化的回波包络，fc为载频，c为光速；在雷达A收到sA_comp(t)信号的同一时刻，雷达B收到目标回波，记为sB(t)，经过脉冲压缩，sB(t)可表达为：其中RBi表示第i个散射点...

【专利技术属性】
技术研发人员：吴文振，张月，徐世友，陈曾平，
申请(专利权)人：中国人民解放军国防科技大学，
类型：发明
国别省市：湖南,43