一种基于波程差补偿的InISAR图像配准方法

周子铂 蒋李兵 王壮

引用本文:
Citation:

一种基于波程差补偿的InISAR图像配准方法

    作者简介: 周子铂(1994–),男,江西九江人。2016年在西安电子科技大学获得学士学位,现在国防科技大学自动目标识别重点实验室攻读硕士研究生,主要研究方向为空间信息获取与处理技术。E-mail: zibo_travel@163.com;蒋李兵(1982–),男,江苏启东人。2014年在国防科技大学获得博士学位,现任国防科技大学电子科学学院自动目标识别重点实验室讲师,主要研究方向为电磁散射建模、微波成像和雷达图像解译。E-mail: jianglibing@nudt.edu.cn;王 壮(1973–),男,陕西西安人。2001年在国防科技大学获得博士学位,现任国防科技大学自动目标识别重点实验室教授,博士生导师,主要研究方向为雷达信息处理、空间目标监视、目标识别。E-mail: wangzhuang@nudt.edu.cn.
    通讯作者: 王壮, wangzhuang@nudt.edu.cn
  • 基金项目:

    国家部委基金

  • 中图分类号: TN95

Image Registration Based on Wave Path Difference Compensation for InISAR

    Corresponding author: Wang Zhuang, wangzhuang@nudt.edu.cn
  • Fund Project: The National Ministries Foundation

    CLC number: TN95

  • 摘要: 逆合成孔径雷达(ISAR)图像配准是干涉逆合成孔径雷达(InISAR)成像领域一个关键的课题,可以实现同一散射点在不同ISAR图像中的对齐,以便于后续的ISAR图像干涉处理。该文分析了ISAR图像失配准的原因,即散射点到不同天线之间的波程差,并据此提出一种基于波程差补偿的方法来实现不同天线ISAR图像之间的精确配准。首先通过调频傅里叶变换估计目标相对于雷达的转速;进一步根据波程差与目标转动角速度的关系构建补偿相位消除散射点到不同天线间的波程差,并通过2维傅里叶变换获得配准之后的ISAR图像。最后利用干涉处理获得目标真实的3维结构。该文方法可以在回波域通过波程差补偿实现ISAR图像配准,配准之后的各散射点在图像中的位置相同;而经过相关法配准之后的ISAR图像中的各散射点之间有一个像素单元的错位,即该方法的配准效果更精确。此外,基于相关法的图像配准方法耗时达到万秒级,而基于该方法的ISAR图像配准时间仅为秒级,即该方法计算效率更高。最终的InISAR 3维成像结果中,该方法的散射点坐标重构误差为0.3034,而基于相关法的成像结果的误差(45.8529)远大于此。因此,基于所提出方法的InISAR 3维成像结果精度更高。
  • 图 1  L型三天线构型

    Figure 1.  Configuration of L-shape three-antenna InISAR imaging system

    图 2  目标平面运动模型

    Figure 2.  Projection of motion model in XOY plane

    图 3  ISAR成像平面转台模型

    Figure 3.  Turntable model of ISAR imaging

    图 4  强散射距离单元块的选择方法

    Figure 4.  The selection of strong scattering range cell block

    图 5  InISAR 3维成像的算法流程

    Figure 5.  Flowchart of the rotational angular velocity estimation algorithm

    图 6  旋转角速度估计结果

    Figure 6.  Estimation results of the Doppler chirp rate

    图 7  角速度w估计的鲁棒性

    Figure 7.  Estimation of w for different SNR

    图 8  经过配准之后的1维距离像

    Figure 8.  Range profiles without mismatch

    图 9  第50次脉冲的1维距离像比较

    Figure 9.  The comparison of range profiles of the 50th pulse

    图 10  未经配准的目标ISAR图像

    Figure 10.  ISAR images without image registration

    图 11  经过配准之后的目标ISAR图像

    Figure 11.  ISAR images after image registration based on the proposed method

    图 12  基于波程差补偿配准方法得到的干涉3维成像结果

    Figure 12.  InISAR imaging result through the proposed method

    图 13  图像配准过程中相关系数分布

    Figure 13.  Distribution of correlation coefficients during registration

    图 14  相关法配准之后的ISAR图像

    Figure 14.  ISAR image after registration based on correlation coefficient

    图 15  基于相关法图像配准最终3维成像结果

    Figure 15.  InISAR imaging result through the method based on correlation efficient

    图 16  天宫一号光学图像

    Figure 16.  Picture of TianGong

    图 17  天宫一号仿真模型

    Figure 17.  Scatter point model of TianGong

    图 18  天宫一号模型InISAR成像结果

    Figure 18.  The imaging result of TianGong model

    表 1  雷达和目标的参数

    Table 1.  Parameters of experiment

    雷达参数 数值 目标参数 数值
    信号载频fc 10 GHz 散射点数 8
    信号带宽BW 1 GHz 目标高度 100 km
    距离分辨率 0.15 m 转动角速度 0.0112 rad/s
    脉冲时宽Tp 0.1 ms 最大尺寸 5.7 m
    脉冲重复频率PRF 100 Hz
    脉冲数量 500
    下载: 导出CSV

    表 2  基于波程差补偿方法配准的ISAR图像中各散射点位置

    Table 2.  Locations of scatterers after image registration based on the proposed method

    散射点标号 天线A 天线B 天线C
    1 (247,218) (247,218) (247,218)
    2 (279,229) (279,229) (279,229)
    3 (229,243) (229,243) (229,243)
    4 (260,248) (260,248) (260,248)
    5 (253,267) (253,267) (253,267)
    6 (242,269) (242,269) (242,269)
    7 (256,286) (256,286) (256,286)
    8 (244,295) (244,295) (244,295)
    下载: 导出CSV

    表 3  相关法配准之后ISAR图像中各散射点位置

    Table 3.  Locations of scatterers after image registration based on correlation method

    散射点标号 A图像 B图像 C图像
    1 (247,218) (246,218) (247,218)
    2 (279,229) (278,229) (278,230)
    3 (229,243) (229,243) (229,243)
    4 (260,248) (260,247) (260,248)
    5 (253,267) (252,266) (253,267)
    6 (242,269) (242,269) (241,270)
    7 (256,286) (256,286) (256,286)
    8 (244,295) (243,295) (244,295)
    下载: 导出CSV
  • [1] Xing M, Wu R, and Bao Z. High resolution ISAR imaging of high speed moving targets[J]. IEE Proceedings-Radar,Sonar and Navigation, 2005, 152(2): 58–67. DOI: 10.1049/ip-rsn:20045084
    [2] Chen C C and Andrews H C. Target-motion-induced radar imaging[J]. IEEE Transactions on Aerospace and Electronic Systems, 1980, AES-16(1): 2–14. DOI: 10.1109/TAES.1980.308873
    [3] Martorella M, Stagliano D, Salvetti F, et al. 3D interferometric ISAR imaging of noncooperative targets[J]. IEEE Transactions on Aerospace and Electronic Systems, 2014, 50(4): 3102–3114. DOI: 10.1109/TAES.2014.130210
    [4] Tian B, Li N, Liu Y, et al.. A novel image registration method for InISAR imaging system[C]. Proceedings of SPIE 9252, Millimetre Wave and Terahertz Sensors and Technology VII, Amsterdam, Netherlands, 2014: 92520U.
    [5] Tian B, Zou J W, Xu S Y, et al. Squint model interferometric ISAR imaging based on respective reference range selection and squint iteration improvement[J]. IET Radar,Sonar&Navigation, 2015, 9(9): 1366–1375. DOI: 10.1049/iet-rsn.2015.0195
    [6] Nasirian M and Bastani M H. A novel model for three-dimensional imaging using interferometric ISAR in any curved target flight path[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(6): 3236–3245. DOI: 10.1109/TGRS.2013.2271875
    [7] 王俊, 文亚亚, 禹娟, 等. 一种基于运动目标ISAR像序列的三维重构方法[J]. 系统仿真学报, 2013, 25(4): 809–816. DOI: 10.16182/j.cnki.joss.2013.04.040Wang Jun, Wen Ya-ya, Yu Juan, et al. Three dimensional reconstruction method for moving target from ISAR sequences[J]. Journal of System Simulation, 2013, 25(4): 809–816. DOI: 10.16182/j.cnki.joss.2013.04.040
    [8] Staglianò D, Lischi S, Massini R, et al.. Soft 3D-ISAR image reconstruction using a dual interferometric radar[C]. 2015 IEEE Radar Conference (RadarCon), Arlington, VA, USA, 2015: 572–576. DOI: 10.1109/RADAR.2015.7131063.
    [9] Wu W Z, Hu P J, Xu S Y, et al. Image registration for InISAR based on joint translational motion compensation[J]. IET Radar,Sonar&Navigation, 2017, 11(10): 1597–1603. DOI: 10.1049/iet-rsn.2017.0140
    [10] 王青松, 瞿继双, 黄海风, 等. 联合实、复相关函数的干涉SAR图像配准方法[J]. 测绘学报, 2012, 41(4): 563–569Wang Qingsong, Qu Jishuang, Huang Haifeng, et al. A method based on integrating real and complex correlation function for InSAR image coregistration[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(4): 563–569
    [11] 曾琪明, 解学通. 基于谱运算的复相关函数法在干涉复图像配准中的应用[J]. 测绘学报, 2004, 33(2): 127–131. DOI: 10.3321/j.issn:1001-1595.2004.02.007Zeng Qi-ming and Xie Xue-tong. A FFT-based complex correlation function method applied to interferometric complex image corregistration[J]. Acta Geodaetica et Cartographica Sinica, 2004, 33(2): 127–131. DOI: 10.3321/j.issn:1001-1595.2004.02.007
    [12] Stone H S, Orchard M T, Chang E C, et al. A fast direct Fourier-based algorithm for subpixel registration of images[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(10): 2235–2243. DOI: 10.1109/36.957286
    [13] 刘承兰, 高勋章, 贺峰, 等. 一种基于相位校正的InISAR图像配准新方法[J]. 国防科技大学学报, 2011, 35(5): 116–122. DOI: 10.3969/j.issn.1001-2486.2011.05.022Liu Cheng-lan, Gao Xun-zhang, He Feng, et al. A novel method for image registration in InISAR imaging based on phase correction[J]. Journal of National University of Defense Technology, 2011, 35(5): 116–122. DOI: 10.3969/j.issn.1001-2486.2011.05.022
    [14] Zhang Q, Yeo T S, Du G, et al. Estimation of three-dimensional motion parameters in interferometric ISAR imaging[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(2): 292–300. DOI: 10.1109/TGRS.2003.815669
  • [1] 路满宋红军罗运华 . 基于调频连续波信号的圆弧式合成孔径雷达成像方法. 雷达学报, 2016, 5(4): 425-433. doi: 10.12000/JR16007
    [2] 种劲松周晓中 . 合成孔径雷达图像海洋内波探测研究综述. 雷达学报, 2013, 2(4): 406-421. doi: 10.3724/SP.J.1300.2013.13012
    [3] 蒋彦雯邓彬王宏强秦玉亮庄钊文 . 基于双频联合处理的太赫兹InISAR成像方法. 雷达学报, 2018, 7(1): 139-146. doi: 10.12000/JR17109
    [4] 周叶剑张磊王虹现邢孟道 . 空间轨道目标的逆合成孔径雷达成像质量分析. 雷达学报, 2017, 6(1): 17-24. doi: 10.12000/JR16136
    [5] 钟雪莲向茂生郭华东陈仁元 . 机载重轨干涉合成孔径雷达的发展. 雷达学报, 2013, 2(3): 367-381. doi: 10.3724/SP.J.1300.2013.13005
    [6] 吴一戎 . 多维度合成孔径雷达成像概念. 雷达学报, 2013, 2(2): 135-142. doi: 10.3724/SP.J.1300.2013.13047
    [7] 张月婷丁赤飚王宏琦胡东辉 . 分析圆柱油罐目标合成孔径雷达图像的新方法. 雷达学报, 2012, 1(2): 190-195. doi: 10.3724/SP.J.1300.2012.20020
    [8] 赵雨露张群英李超纪奕才方广有 . 视频合成孔径雷达振动误差分析及补偿方案研究. 雷达学报, 2015, 4(2): 230-239. doi: 10.12000/JR14153
    [9] 詹学丽王岩飞王超李和平 . 一种用于合成孔径雷达的数字去斜方法. 雷达学报, 2015, 4(4): 474-480. doi: 10.12000/JR14117
    [10] 林世斌李悦丽严少石周智敏 . 平地假设对合成孔径雷达时域算法成像质量的影响研究. 雷达学报, 2012, 1(3): 309-313. doi: 10.3724/SP.J.1300.2012.20035
    [11] 王岩飞刘畅詹学丽韩松 . 无人机载合成孔径雷达系统技术与应用. 雷达学报, 2016, 5(4): 333-349. doi: 10.12000/JR16089
    [12] 李海英张珊珊李世强张华春 . 环境一号C 卫星合成孔径雷达相干性分析. 雷达学报, 2014, 3(3): 320-325. doi: 10.3724/SP.J.1300.2014.13060
    [13] 曾操梁思嘉王威徐青 . 基于频率步进信号的旋转式合成孔径雷达成像方法. 雷达学报, 2014, 3(4): 401-408. doi: 10.3724/SP.J.1300.2014.14043
    [14] 曾涛 . 双基地合成孔径雷达发展现状与趋势分析. 雷达学报, 2012, 1(4): 329-341. doi: 10.3724/SP.J.1300.2012.20093
    [15] 杨汝良戴博伟李海英 . 极化合成孔径雷达极化层次和系统工作方式. 雷达学报, 2016, 5(2): 132-142. doi: 10.12000/JR16013
    [16] 金添 . 叶簇穿透合成孔径雷达增强成像方法. 雷达学报, 2015, 4(5): 503-508. doi: 10.12000/JR15114
    [17] 张珂殊潘洁王然李光祚王宁吴一戎 . 大幅宽激光合成孔径雷达成像技术研究. 雷达学报, 2017, 6(1): 1-10. doi: 10.12000/JR16152
    [18] 禹卫东杨汝良邓云凯赵凤军雷宏 . HJ-1-C 卫星合成孔径雷达载荷的设计与实现. 雷达学报, 2014, 3(3): 256-265. doi: 10.3724/SP.J.1300.2014.14020
    [19] 张文彬邓云凯王宇 . 星地双基合成孔径雷达聚束模式快速BP算法. 雷达学报, 2013, 2(3): 357-366. doi: 10.3724/SP.J.1300.2013.13031
    [20] 何峰杨阳董臻梁甸农 . 曲线合成孔径雷达三维成像研究进展与展望. 雷达学报, 2015, 4(2): 130-135. doi: 10.12000/JR14119
  • 加载中
图(18)表(3)
计量
  • 文章访问数:  213
  • HTML浏览量:  97
  • PDF下载量:  151
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-03
  • 录用日期:  2018-11-06
  • 网络出版日期:  2018-11-21
  • 刊出日期:  2018-12-28

一种基于波程差补偿的InISAR图像配准方法

    通讯作者: 王壮, wangzhuang@nudt.edu.cn
    作者简介: 周子铂(1994–),男,江西九江人。2016年在西安电子科技大学获得学士学位,现在国防科技大学自动目标识别重点实验室攻读硕士研究生,主要研究方向为空间信息获取与处理技术。E-mail: zibo_travel@163.com;蒋李兵(1982–),男,江苏启东人。2014年在国防科技大学获得博士学位,现任国防科技大学电子科学学院自动目标识别重点实验室讲师,主要研究方向为电磁散射建模、微波成像和雷达图像解译。E-mail: jianglibing@nudt.edu.cn;王 壮(1973–),男,陕西西安人。2001年在国防科技大学获得博士学位,现任国防科技大学自动目标识别重点实验室教授,博士生导师,主要研究方向为雷达信息处理、空间目标监视、目标识别。E-mail: wangzhuang@nudt.edu.cn
  • 国防科技大学ATR重点实验室   长沙   410073
基金项目:  国家部委基金

摘要: 逆合成孔径雷达(ISAR)图像配准是干涉逆合成孔径雷达(InISAR)成像领域一个关键的课题,可以实现同一散射点在不同ISAR图像中的对齐,以便于后续的ISAR图像干涉处理。该文分析了ISAR图像失配准的原因,即散射点到不同天线之间的波程差,并据此提出一种基于波程差补偿的方法来实现不同天线ISAR图像之间的精确配准。首先通过调频傅里叶变换估计目标相对于雷达的转速;进一步根据波程差与目标转动角速度的关系构建补偿相位消除散射点到不同天线间的波程差,并通过2维傅里叶变换获得配准之后的ISAR图像。最后利用干涉处理获得目标真实的3维结构。该文方法可以在回波域通过波程差补偿实现ISAR图像配准,配准之后的各散射点在图像中的位置相同;而经过相关法配准之后的ISAR图像中的各散射点之间有一个像素单元的错位,即该方法的配准效果更精确。此外,基于相关法的图像配准方法耗时达到万秒级,而基于该方法的ISAR图像配准时间仅为秒级,即该方法计算效率更高。最终的InISAR 3维成像结果中,该方法的散射点坐标重构误差为0.3034,而基于相关法的成像结果的误差(45.8529)远大于此。因此,基于所提出方法的InISAR 3维成像结果精度更高。

English Abstract

    • 雷达成像是雷达信号处理领域一个非常重要的课题,逆合成孔径雷达(Inverse Synthetic Aperture Radar, ISAR)成像可以有效提升雷达的空间检测和目标识别能力[13]。ISAR图像是目标沿雷达视线方向和目标运动轨迹方向所成平面的投影,不同观测角度下得到的ISAR图像相差比较大,而且ISAR图像不能完全反映目标真实结构,无法提供目标俯仰向的信息。 随着空间监测对于高分辨图像的需求,3维成像技术日趋成为雷达成像领域热点课题[4,5]

      干涉技术首先是在合成孔径雷达(Synthetic Aperture Radar, SAR)成像中被应用的,考虑到ISAR成像为SAR成像的逆过程,自然也就可以将干涉技术引用到ISAR成像领域,发展成为干涉逆合成孔径雷达(Interferometric Inverse Synthetic Aperture Radar, InISAR)3维成像技术[6,7]。InISAR 3维成像系统包含有单天线结构、双天线结构、三天线结构和阵列多天线结构。考虑到单天线和双天线系统结构对于目标运动轨迹的特殊要求以及天线阵列复杂的硬件结构和高昂成本,本文主要研究经典的L型三天线结构。通过使用3部天线,此种3维干涉成像技术可以通过CLEAN算法提取出ISAR图像中的强散射点,并通过同一基线上的两幅ISAR图像的干涉处理得到该散射点的干涉相位[8]。由于干涉相位里面包含有散射点到各雷达天线之间的波程差,经过计算可以得到散射点的3维分布[9]。图像配准是InISAR 3维成像技术里面相当重要的一环,目的是将3幅ISAR图像中的每个散射点对应到相同的像素单元,以确保后续的干涉处理是在相同的散射点之间进行的。

      近些年,很多国内外学者发表了ISAR图像配准和InISAR 3维成像相关的算法和应用的论文。InSAR领域基于相关系数的图像配准方法首先被引用到ISAR图像的配准上[1012],但是SAR和ISAR观测的目标场景有很大的区别,ISAR图像中目标所占的像素单元只是很小的一部分,所以基于相关系数的配准方法并不适用于ISAR图像的配准。文献[13]提出一种基于相位校正的图像配准方法,该方法通过在回波域进行处理,将3部天线的ISAR图像利用各自的参考距离分别完成聚焦,这样在回波就克服了图像失配准的本质原因。该方法对于雷达的测距和测角精度要求比较高,在低信噪比情况下雷达测距和测角精度比较差时,图像配准和3维成像精度较差。在文献[14]中,张群等人提出一种多天线结构3维成像技术,在互相垂直的方向上布置多部天线,利用较长基线对的天线对目标测角,较短的基线对的天线对目标成像,通过长基线测得的目标的较精确的转角对短基线的回波进行补偿,实现ISAR图像的配准,这种成像系统的天线结构相当复杂,给天线布阵带来极大的挑战,并且硬件成本太高。

      本文提出一种基于回波域波程差补偿的图像配准方法,使用调频傅里叶变换估计得到目标的旋转角速度,并根据目标旋转角速度构造波程差补偿相位对雷达回波进行相位补偿,最终经过距离维和方位维的距离压缩获得已经配准的ISAR图像。

      本文其余部分的结构如下:第2节分析InISAR成像系统的信号模型。本文提出的基于波程差补偿的图像配准方法在第3节进行讨论,在第4节给出实验验证并且与基于相关系数的图像配准方法的性能进行比较。第5节为本文的结论。

    • L型三天线构型的干涉3维成像的结构如图1所示,图中包含两个坐标系,即雷达坐标系(A-XYZ)和目标坐标系(O-xyz)。3部天线在雷达坐标系下的位置分别为A(0, 0, 0), B(L, 0, 0), C(0, 0, L),其中L表示成像系统的基线长度,而且基线AB和AC互相垂直。且仅仅A天线为发射天线,A, BC 3部天线都接收目标回波。

      图  1  L型三天线构型

      Figure 1.  Configuration of L-shape three-antenna InISAR imaging system

      目标坐标系原点为目标的几何中心点O,在雷达坐标系下的坐标为 $O\left( {{X_0},{Y_0},{Z_0}} \right)$ ,参考中心点(即目标几何中心点)到3部天线的距离分别为 ${R_{OA}}$ , ${R_{OB}}$ , ${R_{OC}}$ 。点P为目标上的第k个散射点,在目标坐标系下的坐标为 $\left( {{x_k},{y_k},z{}_k} \right)$ , P点到3部天线的距离分别为 ${R_{Ak}}$ , ${R_{Bk}}$ , ${R_{Ck}}$ 。此外,目标沿AO, BOCO方向飞行的速度分别为 ${V_{Ak}}$ , ${V_{Bk}}$ ${V_{Ck}}$ $\hat t$ 为目标成像快时间, ${t_m}$ 为成像慢时间。

      此外,设雷达发射线性调频信号,信号载频为 ${f\!_{\rm{c}}}$ ,带宽为BW,脉冲持续时间为 ${T_{\rm{p}}}$ ,脉冲重复间隔为 ${T_{\rm{r}}}$ ,成像总时间为 ${T_1}$ ,第k个散射点的散射强度为 ${\sigma _k}$ 。发射信号的调频斜率为 $\gamma $ ,雷达波长为 $\lambda $ ,光速为 $c$ 。所有3部雷达均使用A天线测得的参考距离 ${R_{{\rm{ref}}A}}{\rm{ = }}{R_{AO}}$ 进行解线频调处理。设 $ {R_{\Delta Ak}} =$ $ {R_{Ak}} - {R_{{\rm{ref}}A}}$ , ${R_{\Delta Bk}} = {R_{Bk}} - {R_{{\rm{ref}}A}}$ , ${R_{\Delta Ck}} = {R_{Ck}}$ $ - {R_{{\rm{ref}}A}}$ ,然后雷达A, BC经过解线频调处理之后的回波为:

      $ \begin{align} {S_{{\rm{if}}A}}\left( {{\overset{\frown} t} ,{t_m}} \right) =& \sum {\sigma \!_k}{\rm{rect}}\left( {\frac{{\hat t - 2{R_{Ak}}/c}}{{{T_{\rm{p}}}}}} \right)\\ & \cdot \exp \left\{ {{\rm j}\frac{{4{{π}} \gamma }}{c}\left( {\hat t - \frac{{2{R_{{\rm{ref}}A}}}}{c}} \right){R_{\Delta Ak}}} \right\}\\ & \cdot \exp \left\{ {{\rm j}4{{π}} \left( {\frac{{{R_{\Delta Ak}}}}{\lambda } - \frac{{\gamma R_{\Delta Ak}^2}}{{{c\,^2}}}} \right)} \right\} \end{align} $

      $ \begin{align} {S_{{\rm{if}}B}} & \left( {\mathord{\overset{\frown} t} ,{t_m}} \right) \\ \quad &= \sum {\sigma \!_k}{\rm{rect}}\left( {\frac{{\hat t - \left( {{R_{Ak}} + {R_{Bk}}} \right)/c}}{{{T_{\rm{p}}}}}} \right)\\ \quad & \quad\cdot \exp \left\{ {{\rm j}\frac{{2{{π}} \gamma }}{c}\left( {\hat t - \frac{{2{R_{{\rm{ref}}A}}}}{c}} \right)\left( {{R_{\Delta Ak}} + {R_{\Delta Bk}}} \right)} \right\} \\ \quad &\quad\cdot \exp \left\{ {\rm j}2{{π}} \frac{{\left( {{R_{\Delta Ak}} + {R_{\Delta Bk}}} \right)}}{\lambda }\right. \\ \quad &\quad \left.- {\rm j}{{π}} \frac{{\gamma {{\left( {{R_{\Delta Ak}} + {R_{\Delta Bk}}} \right)}^2}}}{{{c^2}}} \right\} \end{align} $

      $ \begin{align} {S_{{\rm{if}}C}} & \left({\mathord{\overset{\frown} t} ,{t_m}} \right) \\ \quad & = \sum {\sigma \!_k}{\rm{rect}}\left( {\frac{{\hat t - \left( {{R_{Ak}} + {R_{Ck}}} \right)/c}}{{{T_{\rm{p}}}}}} \right)\\ \quad & \quad\cdot \exp \left\{ {{\rm j}\frac{{2{{π}} \gamma }}{c}\left( {\hat t - \frac{{2{R_{{\rm{ref}}A}}}}{c}} \right)\left( {{R_{\Delta Ak}} + {R_{\Delta Ck}}} \right)} \right\} \\ \quad & \quad \cdot \exp \left\{ {\rm j}2{{π}} \frac{{\left( {{R_{\Delta Ak}} + {R_{\Delta Ck}}} \right)}}{\lambda } \right.\\ \quad &\quad\left. - {\rm j}{{π}} \frac{{\gamma {{\left( {{R_{\Delta Ak}} + {R_{\Delta Ck}}} \right)}^2}}}{{{c^2}}} \right\} \end{align} $

      干涉逆合成孔径雷达成像技术是高分辨2维ISAR成像和干涉技术的结合,因此,高分辨的ISAR成像是InISAR 3维成像的前提条件。经过文献研究可知,ISAR 2维成像技术已经比较成熟,根据传统的平动补偿方法和脉冲压缩方法可以得到目标比较清晰的2维ISAR图像。在3部天线都完成2维ISAR成像之后,需要通过图像配准方法将同一散射点在3幅不同ISAR图像当中对齐到同一像素单元,以确保后续的干涉处理过程是对相同散射点进行干涉处理,以解算出各散射点的空间位置。根据ISAR成像理论,A, B, C 3部天线所成的ISAR复图像分别为:

      $ \begin{align} {S\!_A}\left( {f,{f\!_m}} \right) =& \sum\limits_k {{{\sigma \!_k}}{T_{\rm{p}}}{T_1}} \sin\!{\rm c}\left[ {{T_{\rm{p}}}\left( {{f\!_i} - \frac{{2\gamma }}{c}{R_{\Delta Ak}}} \right)} \right]\\ & \cdot \sin\!{\rm c}\left[ {{T_1}\left( {{f\!_m} - \frac{{2{V_{Ak}}}}{\lambda }} \right)} \right]\\ & \cdot \exp \left( { - {\rm j}\frac{{4{{π}} {R_{\Delta Ak0}}}}{\lambda }} \right) \end{align} $

      $ \begin{align} {S\!_B}\left( {f,{f\!_m}} \right) =& \sum\limits_k {{{\sigma \!_k}}{T_{\rm{p}}}{T_1}} \\ & \cdot\sin\!{\rm c}\left[ {{T_{\rm{p}}}\left( {{f\!_i} - \frac{\gamma }{c}\left( {{R_{\Delta Ak}} + {R_{\Delta Bk}}} \right)} \right)} \right]\\ & \cdot \sin\!{\rm c}\left[ {{T_1}\left( {{f\!_m} - \frac{{\left( {{V_{Ak}} + {V_{Bk}}} \right)}}{\lambda }} \right)} \right]\\ & \cdot \exp \left( { - {\rm j}\frac{{2{{π}} \left( {{R_{\Delta Ak0}} + {R_{\Delta Bk0}}} \right)}}{\lambda }} \right) \end{align} $

      $ \begin{align} {S\!_C}\left( {f,{f\!_m}} \right) =& \sum\limits_k {{{\sigma \!_k}}{T_{\rm{p}}}{T_1}}\\ & \cdot \sin\!{\rm c}\left[ {{T_{\rm{p}}}\left( {{f\!_i} - \frac{\gamma }{{\rm c}}\left( {{R_{\Delta Ak}} + {R_{\Delta Ck}}} \right)} \right)} \right]\\ & \cdot \sin\!{\rm c}\left[ {{T_1}\left( {{f\!_m} - \frac{{\left( {{V_{Ak}} + {V_{Ck}}} \right)}}{\lambda }} \right)} \right]\\ & \cdot \exp \left( { - {\rm j}\frac{{2{{π}} \left( {{R_{\Delta Ak0}} + {R_{\Delta Ck0}}} \right)}}{\lambda }} \right) \end{align} $

      根据式(4)和式(5)可知,同一散射点在A, B两幅ISAR图像当中的位置有一定的偏移:

      $ \begin{align} \Delta f &= \frac{\gamma }{c}\left( {{R_{\Delta Ak0}} - {R_{\Delta Bk0}}} \right) = \frac{\gamma }{c}\left( {{R_{\Delta Ak}} - {R_{\Delta Bk}}} \right)\Big| {_{{t_m} = 0}} \\ &= \frac{\gamma }{c}\left( {{R_{Ak}} - {R_{Bk}}} \right)\Big| {_{{t_m} = 0}} \end{align} $

      $ \begin{align} \Delta {f\!_m} & = \frac{{{V_{Ak}} - {V_{Bk}}}}{\lambda } \\ & = \frac{{\left( {{R_{Ak}} - {R_{Bk}}} \right)\left| {_{{t_m} = {T_1}} - \left( {{R_{Ak}} - {R_{Bk}}} \right)\Big| {_{{t_m} = 0}}} \right.}}{{\lambda {T_1}}} \end{align} $

      根据式(7)和式(8),代入X波段宽带成像雷达典型系统参数 $\rm{BW} = 1\;GHz$ , $c = 3 \times {10^8}$ , $\lambda = 0.03\;{\rm m}$ 以及散射点不同时间到雷达距离,可得到此时A, B天线ISAR图像在距离向和方位向偏移的距离单元数:

      $ \begin{align} \Delta {{{N}}_f} = & \frac{{\rm{BW}}}{c}\left( {{R_{Ak}} - {R_{Bk}}} \right)\Big| {_{{t_m} = 0}} \\ & < 6.67 \times {10^{ - 5}} < < 1 \end{align} $

      $ \begin{align} \Delta {N_{{f_m}}} & = \frac{{\left( {{R_{Ak}} - {R_{Bk}}} \right)\Big| {_{{t_m} = {T_1}} - \left( {{R_{Ak}} - {R_{Bk}}} \right)\Big| {_{{t_m} = 0}} } }}{\lambda } \\ & \approx 16.67 > 1 \end{align} $

      由式(7)和式(8)可知,ISAR图像距离向的失配主要由成像起始时刻目标( ${t_m} = 0$ )到不同雷达的波程差 $\left(\left( {{R_{\Delta Ak0}} - {R_{\Delta Bk0}}} \right)\left| {_{{t_m} =0}}\right.\right) $ 引起,而方位向的失配则是由成像时间内积累的波程差 $\left(\left( {{R_{Ak}} - {R_{Bk}}} \right)\left| {_{{t_m} = {T_1}}}\right.\right.$ $ \left. {- \left( {{R_{Ak}} - {R_{Bk}}} \right)\left| {_{{t_m} = 0}}\right. } \right)$ 引起。根据式(9)和式(10)的计算结果可知,ISAR图像在距离向的失配量比较小,但是方位向的失配量已经达到十几个分辨单元。综上所述,当3部天线均采用A天线测得的到目标的距离进行解线频调处理时,ISAR图像之间的失配准现象比较严重。

      因此,在回波域将由天线位置不同而引起的波程差补偿掉,就可以实现各天线ISAR图像的配准。在完成图像配准之后,对两幅ISAR图像进行干涉处理,进而解算出目标中各散射点的3维位置关系。

      $ \begin{align} \Delta {\varphi _{{{AB}}}} & = - {\varphi _A} + {\varphi _B} = \frac{{2{{π}} }}{\lambda }\left( {{R_{\Delta Ak0}} - {R_{\Delta Bk0}}} \right) \\ & = \frac{{2{{π}} }}{\lambda }\frac{{2L\left( {{X_0} + {x_k}} \right) - {L^2}}}{{{R_{Ak}} + {R_{Bk}}}} \end{align} $

      $ \begin{align} \Delta {\varphi _{{{AC}}}} & = - {\varphi _A} + {\varphi _C} = \frac{{2{{π}} }}{\lambda }\left( {{R_{\Delta Ak0}} - {R_{\Delta Ck0}}} \right) \\ & = \frac{{2{{π}} }}{\lambda }\frac{{2L\left( {{Y_0} + {y_k}} \right) - {L^2}}}{{{R_{Ak}} + {R_{Bk}}}} \end{align} $

    • 从式(9)和式(10)可见,当雷达B和雷达C的回波都用雷达A测得的参考距离 ${R_{{\rm{ref}}A}} = {R_{AO}}$ 进行解线频调处理并成像时,B图像和A图像以及C图像和A图像之间的失配准现象比较严重。但是,用雷达B和雷达C到目标的距离分别对各自回波进行解线频调处理时,即 ${R_{{\rm{ref}}B}} = {R_{BO}}$ , ${R_{{\rm{ref}}C}} = {R_{CO}}$ ,图像之间的失配准现象会得到较好的改善。这样可得 $R'\!\!_{\Delta Bk} = {R_{Bk}} - {R_{BO}}$ , $R'\!\!_{\Delta Ck} = {R_{Ck}} - {R_{CO}}$ 。最终,距离向和方位向的偏移量为:

      $ \begin{align} \Delta {f'}\!\!_m &= \frac{1}{\lambda }\Big[ {\left( {{R_{\Delta Ak}} - {R'}\!\!_{\Delta Bk}} \right)\Big| {_{{t_m} = {\rm{ }}{T_1}}} - \left( {{R_{\Delta Ak}} - {R'}\!\!_{\Delta Bk}} \right)\Big| {_{{t_m} = {\rm{ }}0}} } \Big] \\ & = \frac{1}{\lambda }\left\{ {\bigr[ {\left( {{R_{Ak}} - {R_{Bk}}} \right) - \left( {{R_{AO}} - {R_{BO}}} \right)} \bigr]\left| {_{{t_m} = {\rm{ }}{T_1}}} \right. - \left[ {\left( {{R_{Ak}} - {R_{Bk}}} \right) - \left( {{R_{AO}} - {R_{BO}}} \right)} \right]\left| {_{{t_m} = {\rm{ }}0}} \right.} \right\} \\ &< \frac{1}{\lambda }\bigr\{ {\left[ {\left( {{R_{Ak}} - {R_{Bk}}} \right) - \left( {{R_{AO}} - {R_{AO}}} \right)} \right]\left| {_{{t_m} = {\rm{ }}{T_1}}} \right. - \left[ {\left( {{R_{Ak}} - {R_{Bk}}} \right) - \left( {{R_{AO}} - {R_{AO}}} \right)} \right]\left| {_{{t_m} = {\rm{ }}0}} \right.} \bigr\} \\ & = \frac{1}{\lambda }\bigr[ {\left( {{R_{\Delta Ak}} - {R_{\Delta Bk}}} \right)\left| {_{{t_m} = \; {\rm{ }}{T_1}}} \right. - \left( {{R_{\Delta Ak}} - {R_{\Delta Bk}}} \right)\left| {_{{t_m} = {\rm{ }}0}} \right.} \bigr] \end{align} $

      $ \begin{align} \Delta {f'} =& \frac{\gamma }{c}\left( {{R_{\Delta Ak0}} - {R'}\!\!_{\Delta Bk0}} \right) = \frac{\gamma }{c}\bigr[ {\left( {{R_{Ak}} - {R_{Bk}}} \right){\rm{ - }}\left( {{R_{AO}} - {R_{BO}}} \right)} \bigr]\left| {_{t = {\rm{ }}0}} \right. \\ <& \frac{\gamma }{c}\bigr[ {\left( {{R_{Ak}} - {R_{Bk}}} \right){\rm{ - }}\left( {{R_{AO}} - {R_{AO}}} \right)} \bigr]\left| {_{{t_m} = {\rm{ }}0}} \right. = \frac{\gamma }{c}\left[ {{R_{Ak}} - {R_{Bk}}} \right]\left| {_{{t_m} = {\rm{ }}0}} \right. \end{align} $

      由式(13)和式(14)可知,采用3部雷达各自到目标的距离对各自雷达回波进行解线频调处理并成像,便可以有效减小距离向和方位向的偏移。因此,雷达B和雷达C需要补偿的波程差分别为:

      $\Delta R_m^B = {R_{AO}} - {R_{BO}}$

      $\Delta R_m^C = {R_{AO}} - {R_{CO}}$

      图2描述了远场目标在XOY平面的运动模型,目标在XOY平面以 ${w^{AB}}$ 做匀速转动。图中 ${O_1}$ ${O_2}$ 分别为不同时刻目标的中心位置。 ${O'}$ 为基线AB的中点,并且直线 ${O_1}{O'}$ 垂直于X轴。P为目标上任意一散射点, ${O'}$ 到散射点P的连接线与直线 ${O_1}{O'}$ 之间的夹角为 $\theta \left( {{t_m}} \right)$ ,此处 ${t_m}$ 为成像慢时间, $0 \le {t_m} \le {T_1}$ 。散射点P到雷达A和雷达B的距离分别为 ${R_1}\left( {{t_m}} \right)$ ${R_2}\left( {{t_m}} \right)$ ,均为慢时间 ${t_m}$ 的函数。平面YOZ内的运动情况和XOY平面内的运动情况类似,且其转动角速度为 ${w^{AC}}$

      图  2  目标平面运动模型

      Figure 2.  Projection of motion model in XOY plane

      图2中散射点与雷达位置之间的几何关系,可计算得到散射点P到雷达A和雷达B的距离差为:

      $ \begin{align} {R_1}\left( {{t_m}} \right) - {R_2}\left( {{t_m}} \right) &= L\sin \bigr( {\theta \left( {{t_m}} \right)}\bigr)\\ &=L \bullet \sin \left( {\theta \left( {{t_0}} \right) + {w^{AB}}{t_m}} \right) \end{align} $

      同理,参考距离差:

      $ \begin{align} \Delta R_m^B &= {R_{A{O_1}}}\left( {{t_m}} \right) - {R_{B{O_1}}}\left( {{t_m}} \right) \\ &= L\sin \bigr( {\theta \left( {{t_m}} \right){\rm{ - }}\theta \left( {{t_0}} \right)} \bigr) \approx L \bullet {w^{AB}}{t_m} \end{align} $

      其中, ${R_{A{O_1}}}\left( {{t_m}} \right)$ ${R_{B{O_1}}}\left( {{t_m}} \right)$ 分别为雷达A和雷达B的参考距离。同样,雷达A和雷达C之间的参考距离差为 $\Delta R_m^C \approx L \bullet {w^{AC}}{t_m}$

      经过上述分析可知,ISAR图像配准问题可由波程差的估计问题转化为目标角度转动参数的估计问题,可以通过自发自收的雷达A对目标的转动角速度进行估计。

      经过平动补偿之后,目标可被视为转台目标,具体转台模型如图3所示。对于相对雷达做匀速转动的目标,目标上散射点到雷达距离为:

      图  3  ISAR成像平面转台模型

      Figure 3.  Turntable model of ISAR imaging

      $ \begin{align} {R_{Ak}} &=\! {R_{AO}} + {R_{\Delta Ak}} = {R_{AO}} + {R_{Ok}}\cos \left( {{\theta _0} + w{t_m}} \right) \\ & = \!{R_{AO}} \!+\! \sqrt {\!x_k^2 \!+\! y_k^2} \left( {\cos {\theta _0}\cos w{t_m} \!-\! \sin {\theta _0}\sin w{t_m}} \right) \\ & =\! {R_{AO}} + {y_k}\cos w{t_m} - {x_k}\sin w{t_m} \end{align} $

      雷达A的经过剩余视频相位和包络斜置项消除之后的1维距离像为:

      $ \begin{align} {S_{{A}}}\left( {f,{t_m}} \right) = & \sum {} {A_{\rm p}}{T_{\rm p}}\sin\!{\rm c}\left( {{T_{\rm p}}\left( {f{\rm{ + }}\frac{{2\gamma }}{c}{R_{\Delta Ak}}} \right)} \right)\\ & \cdot \exp \left( { - {\rm j}\frac{{4{{π}} }}{\lambda }{R_{\Delta Ak}}} \right){\rm rect}\left( {\frac{{{t_m}}}{{{T_1}}}} \right) \end{align} $

      其中, ${R_{\Delta Ak}} = {y_k}\cos w{t_m} - {x_k}\sin w{t_m}$ ,因此1维距离像的相位项为:

      $ {\rm Phase} = \exp \left\{ { -{\rm j}\frac{{4{{π}} \left( {{y_k}\cos w{t_m} - {x_k}\sin w{t_m}} \right)}}{\lambda }} \right\} $

      (21)

      因此,当成像时间内目标转角比较小的情况下,回波多普勒频率为:

      $ \begin{align} {f_{{\rm{dci}}}} & = \frac{1}{{2{{π}} }}\frac{{d{\rm Phase}}}{{d{t_m}}} = - \frac{{2w}}{\lambda }\left( {{y_k}\sin w{t_m} + {x_k}} \right) \\ & \approx - \frac{{2w}}{\lambda }{x_k} - \frac{{2{y_k}{w^2}}}{\lambda }{t_m} \end{align} $

      由式(22)可见,回波在慢时间域为线性调频信号,且调频率为:

      $k = - \frac{{2{y_0}{w^2}}}{\lambda }$

      式(23)表明,在给定雷达波长 $\lambda $ 和散射点的纵坐标 ${y_k}$ 之后,目标回波在多普勒域的调频斜率与目标的转动角速度成一一对应关系。换句话说目标旋转角速度的估计问题就转化为回波多普勒域的调频率估计问题。因此,调频傅里叶变换可以作为目标回波慢时间域调频率k估计的有效方法。回波强散射点的1维方位像的调频傅里叶变换的定义式为:

      $ F\left( {f,\gamma } \right) \!=\! \int {{S_{AA}}\left( {f,{t_m}} \right)\exp \left( {{\rm j}2{{π}} \bullet\frac{1}{2}{\gamma _m}{t_m}\!\!\!^2} \right){\rm d}{t_m}} $

      其中, ${S_{AA}}\left( {f,{t_m}} \right)$ 表示雷达A回波强散射点所在的单元的1维方位像,且调频傅里叶变换过程中需要补偿的调频斜率的范围为:

      ${\gamma _m} = \frac{{2{y_k}{w^2}}}{\lambda } \in \left\{ {\frac{{2{y_k}w_{\min }^2}}{\lambda },\frac{{2{y_k}w_{\max }^2}}{\lambda }} \right\}$

      $\left. \begin{array}{l} H = - \sum\limits_{m = 1}^M {\sum\limits_{n = 1}^N {P\left( {m,n} \right)\ln P} } \left( {m,n} \right)\\ P\left( {m,n} \right) = \frac{{\left| {F\left( {m,n} \right)} \right|}}{{\displaystyle\sum\limits_{m = 1}^M {\displaystyle\sum\limits_{n = 1}^N {\left| {F\left( {m,n} \right)} \right|} } }} \end{array} \right}$

      在每一个补偿的调频斜率 ${\gamma _m}$ 下,计算调频傅里叶变换结果 $F\left( {f,\gamma } \right)$ 的熵,以熵最小为准则选择最佳补偿的调频斜率为回波多普勒域调频斜率的估计值。设调频傅里叶变换的结果为 $F\left( {m,n} \right)$ ,则熵的计算式如式(26)所示。进一步为了估计结果的稳定性,此处选择强散射点所在距离单元的前后多个距离单元(调频斜率变化不大)的强散射单元块的1维方位像作为整体对目标回波多普勒频率的调频斜率进行估计。具体的强散射单元块的选择方法如图4所示。

      图  4  强散射距离单元块的选择方法

      Figure 4.  The selection of strong scattering range cell block

      其中,强散射单元块由2n+1个1维方位像组成,分别为强散射点单元所在的1维方位像,以及其前后各n个距离单元的1维方位像。具体的强散射单元块的组成如图4中阴影部分所示。此处选择强散射单元块遵循两个原则:(1)充分利用强散射单元的峰值的旁瓣信息;(2)同一个强散射单元块不同时包含两个强散射单元。

      在各个强散射单元块估计得到调频斜率之后,再根据调频斜率与散射点距离 ${y_k}$ 成线性关系的特点,对多个强散射单元块估计得到的调频斜率进行最小二乘拟合,得到的直线的斜率K与转动角速度成唯一对应关系:

      $w = \sqrt {\left| K \right|\lambda /2} $

    • 进一步根据上述推导得到的需要补偿的波程差,可以构造雷达B和雷达C回波分别需要补偿的相位为:

      $ {\rm Pha}_m^B =\exp \left\{ { -{\rm j}2{{π}} \gamma \hat t\Delta R_m^B/c} \right\}\exp \left\{ { - {\rm j}2{{π}} \Delta R_m^B/\lambda } \right\} $

      (28)

      $ {\rm Pha}_m^C= \exp \left\{ { - {\rm j}2{{π}} \gamma \hat t\Delta R_m^C/c} \right\} \exp \left\{ { - {\rm j}2{{π}} \Delta R_m^C/\lambda } \right\} $

      (29)

      然后,对雷达B和雷达C经过解线频调处理的接收回波进行相位补偿:

      $S_{{{\rm{if}}B}}\!\!\!\!\!'\ \left( {\mathord{\overset{\frown} t} ,{t_m}} \right) = {S_{{\rm{if}}B}}\left( {\hat t,{t_m}} \right) \bullet {\rm{Pha}}_m^B$

      $S_{{{\rm{if}}C}}\!\!\!\!'\ \left( {\mathord{\overset{\frown} t} ,{t_m}} \right) = {S_{{\rm{if}}C}}\left( {\hat t,{t_m}} \right) \bullet {\rm{Pha}}_m^C$

      最后,对雷达B和雷达C的经过相位补偿之后回波进行成像处理,最终雷达B和雷达C所成的ISAR图像与雷达A的ISAR图像是完全配准的。因为ISAR图像配准仅仅是InISAR 3维成像的中间步骤,所以在经过上述方法获得已经配准的ISAR图像之后,利用干涉技术来提取目标的干涉相位以及各散射点的坐标信息:

      $ \begin{align} {\Delta '}{\varphi _{AB}} &= - {\varphi _A} + {\varphi _B} = \frac{{2{{π}} }}{\lambda }\left( {{R_{\Delta Ak0}} - {R'}\!\!_{\Delta Bk0}} \right) \\ &= \frac{{2{{π}} }}{\lambda }\frac{{2L{x_k}}}{{{R_{Ak}} + {R_{Bk}}}} \approx \frac{{2{{π}} }}{\lambda }\frac{{L{x_k}}}{{{R_{AO}}}} \end{align} $

      $ \begin{align} {\Delta '}{\varphi _{AC}} &= - {\varphi _A} + {\varphi _C} = \frac{{2{{π}} }}{\lambda }\left( {{R_{\Delta Ak0}} - {R'}\!\!_{\Delta Ck0}} \right) \\ &= \frac{{2{{π}} }}{\lambda }\frac{{2L{z_k}}}{{{R_{Ak}} + {R_{Bk}}}} \approx \frac{{2{{π}} }}{\lambda }\frac{{L{z_k}}}{{{R_{AO}}}} \end{align} $

      综上所述,整个InISAR 3维成像算法的流程如图5所示。

      图  5  InISAR 3维成像的算法流程

      Figure 5.  Flowchart of the rotational angular velocity estimation algorithm

    • 为了验证所提出方法的有效性,本部分就所提出的方法进行仿真实验,并与基于相关系数的ISAR图像配准算法进行比较。InISAR成像系统结构如图1所示,基线长度L=10 m。雷达和目标参数如表1所示。

      雷达参数 数值 目标参数 数值
      信号载频fc 10 GHz 散射点数 8
      信号带宽BW 1 GHz 目标高度 100 km
      距离分辨率 0.15 m 转动角速度 0.0112 rad/s
      脉冲时宽Tp 0.1 ms 最大尺寸 5.7 m
      脉冲重复频率PRF 100 Hz
      脉冲数量 500

      表 1  雷达和目标的参数

      Table 1.  Parameters of experiment

    • 目标回波多普勒域调频斜率 ${\gamma _{\rm{comp}}}$ 估计结果如图6所示,蓝色线为调频傅里叶变换之后用最小熵法搜索得到的各距离单元块的调频斜率估计结果,红色线为最小二乘拟合的结果。计算红线的斜率即可计算得到目标的旋转角速度: $w = \sqrt { - \lambda \left| K \right|/2} = 0.{\rm{0108}}$ 图7所示为不同信噪比条件下目标旋转角速度的估计效果,由图可见,在信噪比较高时,估计效果比较稳定,受噪声影响较小。

      图  6  旋转角速度估计结果

      Figure 6.  Estimation results of the Doppler chirp rate

      图  7  角速度w估计的鲁棒性

      Figure 7.  Estimation of w for different SNR

    • 图8所示为各天线经过波程差补偿的回波进行距离压缩得到的目标的1维距离像。显然,3部天线得到的3幅1维距离像极其相似,强散射点基本位于相同的距离单元,忽略噪声的影响,可以认为距离向上基本实现配准。另外,图9给出了雷达A, BC回波的第50次脉冲对齐之后1维距离像的比较结果,峰值点表示散射点位置,并且峰值对齐效果较好,可见散射点在距离向基本实现精确对齐。

      图  8  经过配准之后的1维距离像

      Figure 8.  Range profiles without mismatch

      图  9  第50次脉冲的1维距离像比较

      Figure 9.  The comparison of range profiles of the 50th pulse

      对目标回波的1维距离像进行方位向压缩可以得到目标的2维ISAR图像,3部天线分别各自对目标所成的未进行配准的ISAR图像如图10所示。由于3部天线观测目标的视角不同,且仅仅雷达A发射电磁波,从而雷达B, C接收到的回波的波程与雷达A接收到的回波的波程不同,从而3幅图像会有一定的差别。由图10可见,3幅ISAR图像在距离向没有明显的平移,但是在方位向有明显的错位,即3幅ISAR图像各散射点之间未对齐。考虑到ISAR图像间的错配主要是因为目标上各散射点到不同雷达具有一定的波程差,因此提出一种基于波程差补偿的ISAR图像配准方法。经过所提出的回波域的波程差补偿方法的处理之后,3部天线所成的ISAR图像如图11所示。观察所成图像可以发现无论是在距离向还是方位向,3幅ISAR图像之间均未发生明显的错位。通过峰值检测进一步提取配准之后的ISAR图像当中的强散射点的位置,如表2所示。括号中前一个数为散射点所在的方位向位置,后一个数为散射点在距离向的位置。显然,各强散射点在3幅ISAR图像当中的位置都是一样的,图像间的失配准被消除,即A, B, C 3幅图像完成配准。

      图  10  未经配准的目标ISAR图像

      Figure 10.  ISAR images without image registration

      图  11  经过配准之后的目标ISAR图像

      Figure 11.  ISAR images after image registration based on the proposed method

      散射点标号 天线A 天线B 天线C
      1 (247,218) (247,218) (247,218)
      2 (279,229) (279,229) (279,229)
      3 (229,243) (229,243) (229,243)
      4 (260,248) (260,248) (260,248)
      5 (253,267) (253,267) (253,267)
      6 (242,269) (242,269) (242,269)
      7 (256,286) (256,286) (256,286)
      8 (244,295) (244,295) (244,295)

      表 2  基于波程差补偿方法配准的ISAR图像中各散射点位置

      Table 2.  Locations of scatterers after image registration based on the proposed method

    • 因为ISAR图像配准仅仅是干涉ISAR 3维成像的中间步骤,在完成ISAR图像配准之后还需要进行干涉处理,并根据式(32)和式(33)求解得到目标上各散射点的x坐标和z坐标,y坐标根据各散射点径向距离即可得到。最终,目标的干涉ISAR 3维成像结果如图12所示,其中蓝色“*”点表示目标的真实结构,红色圆圈“○”表示使用波程差补偿方法得到的3维成像结果。其中,红色圆圈与对应的蓝色“*”点的 $\left( {x,y,z} \right)$ 3个坐标的差的平均值为0.3034。可见,提出的方法可以很好地实现L型天线结构中3部天线所成ISAR图像的配准,并得到比较好的3维成像结果。

      图  12  基于波程差补偿配准方法得到的干涉3维成像结果

      Figure 12.  InISAR imaging result through the proposed method

    • 按传统RD成像方法对回波依次进行距离向脉冲压缩,包络对齐,初相校正,方位向脉冲压缩得到目标的ISAR图像。然后分别对BC的图像逐步进行距离向和方位向的平移,并计算相关系数。该仿真实验在设备参数为Intel(R) Core(TM) i7 7700 CPU @ 3.6 GHz and 16 GB RAM计算机的MATLAB环境中进行,相关系数计算结果如图13所示。其中A, B相关系数的计算用时为15916 s, A, C相关系数的计算用时为15740 s;当B图像距离向平移量为0.8,方位向平移量为18.3时,AB图像的相关系数达到最大;当C图像距离向平移量为0.7,方位向平移量为10.8时,AC图像的相关系数达到最大。在实验过程中,为了达到十分之一个像素单元的配准精度,需要对原始图像进行10倍插值处理,另外,还需要搜索距离向和方位向的最佳平移量,这将极大地增加相关法的计算量,因此基于相关系数的图像配准方法的计算时间较长。根据BC图像的最佳平移量对BC图像进行平移,平移过后的ISAR图像如图14所示,从左至右依次为A, B, C图像。同样,可以使用峰值检测的方法提取出3幅图像中的强散射点,如表3所示。根据表3可知,有的散射点在3幅图像中的位置并不是同一个像素单元,8个强散射点并没有完全对齐,可见相关法的配准精度相对较差。因为ISAR图像中目标所占的区域相对较小,大部分像素单元都是被噪声占据,目标散射信息对相关系数的贡献值相对较小,所以两幅图像的相关系数更有可能被噪声所影响,尤其是在信噪比较低的观测条件下。最终导致基于相关系数的ISAR图像配准方法的配准精度较差,容易形成错配现象。

      图  13  图像配准过程中相关系数分布

      Figure 13.  Distribution of correlation coefficients during registration

      图  14  相关法配准之后的ISAR图像

      Figure 14.  ISAR image after registration based on correlation coefficient

      散射点标号 A图像 B图像 C图像
      1 (247,218) (246,218) (247,218)
      2 (279,229) (278,229) (278,230)
      3 (229,243) (229,243) (229,243)
      4 (260,248) (260,247) (260,248)
      5 (253,267) (252,266) (253,267)
      6 (242,269) (242,269) (241,270)
      7 (256,286) (256,286) (256,286)
      8 (244,295) (243,295) (244,295)

      表 3  相关法配准之后ISAR图像中各散射点位置

      Table 3.  Locations of scatterers after image registration based on correlation method

      目标的干涉ISAR 3维图像结果如图15所示,其中,红色圆圈与对应的蓝色“*”点的3个坐标的差的平均值为45.8529。可见基于相关法的图像配准精度比较差,导致最终的InISAR 3维图像与目标的真实结构相差甚远。综上所述,基于波程差补偿的方法通过在回波域进行波程差补偿消除散射点到不同天线之间的波程差进而实现ISAR图像配准,再进行干涉处理获得目标的3维结构,最终的成像结果表明本文方法的有效性。

      图  15  基于相关法图像配准最终3维成像结果

      Figure 15.  InISAR imaging result through the method based on correlation efficient

    • 为进一步验证所提出算法的稳定性,此部分用“天宫一号”卫星模型按照上述算法流程进行InISAR成像仿真实验。具体目标结构以及仿真模型如图16图17所示。

      图  16  天宫一号光学图像

      Figure 16.  Picture of TianGong

      图  17  天宫一号仿真模型

      Figure 17.  Scatter point model of TianGong

      根据上述算法流程,对“天宫一号”散射点模型进行InISAR成像,成像结果如图18所示,其中蓝色“*”点表示目标的真实结构,红色圆圈“○”表示使用波程差补偿方法得到的3维成像结果。因为ISAR图像配准只是InISAR成像的中间步骤,因此这里直接采用最终的InISAR成像结果来衡量算法性能。由成像结果在XOZ平面的投影可见,目标真实结构的投影与成像结果的投影比较一致,成像结果在XOZ平面与真实结构吻合得比较好;再根据YOZ平面的结构可知,干涉结果在y坐标方向上也基本吻合,误差较小,所以,总体而言最终的3维成像质量比较好。可见,所提出的算法对于较复杂的散射点模型也能得出较好的成像效果,由此可见所提出的基于波程差补偿的ISAR图像配准方法可以很好地实现ISAR图像配准,并获得质量较好的InISAR 3维成像结果。

      图  18  天宫一号模型InISAR成像结果

      Figure 18.  The imaging result of TianGong model

    • 本文提出一种基于波程差补偿的ISAR图像配准方法,通过在回波域的相位补偿消除各散射点到不同天线之间的波程差,最终完成ISAR图像的配准。实验结果证明,相比于基于相关系数的ISAR图像配准方法而言,本文方法可以在回波域直接完成ISAR图像的配准,并且配准精度较高,计算负担较轻。最终获得的目标真实的3维成像结果也证实该ISAR图像配准方法更适合被应用于InISAR 3维成像,并进一步通过复杂目标的InISAR 3维成像结构验证了所提出方法的稳定性。

参考文献 (14)

目录

    /

    返回文章
    返回