星载高分辨频率步进SAR成像技术

龙腾 丁泽刚 肖枫 王岩 李喆

引用本文:
Citation:

星载高分辨频率步进SAR成像技术

    作者简介:
    龙 腾(1968–),男,湖北黄陂人。北京理工大学教授、博士生导师。现任北京理工大学副校长、雷达技术研究所所长。研究方向为新体制雷达与实时信息处理。E-mail: longteng@bit.edu.cn;
    丁泽刚(1980–),男,河南信阳人,博士。现担任北京理工大学研究员、博士生导师,主要研究方向为新体制雷达成像机理、成像处理和图像信息提取。E-mail: z.ding@bit.edu.cn;
    肖 枫(1993–),男,内蒙古包头人。北京理工大学雷达技术研究所在读博士研究生,主要研究领域为多模式星载SAR成像技术、SAR运动补偿技术。E-mail: 1120646612@qq.com;
    王 岩(1989–),男,河北沧州人。博士后,现担任北京理工大学副研究员、硕士生导师,主要研究方向为新体制雷达系统、成像、干涉和极化应用。E-mail: yan_wang@bit.edu.cn;
    李 喆(1996–),男,河北任丘人。北京理工大学信息与电子学院硕士研究生,主要研究方向星载合成孔径雷达系统设计与成像。E-mail: 3120180770@bit.edu.cn.
    通讯作者: 丁泽刚 z.ding@bit.edu.cn; 龙腾 longteng@bit.edu.cn
  • 基金项目:

    国家杰出青年科学基金(61625103),国家自然科学基金重点项目(11833001)

  • 中图分类号: TN 958

Spaceborne High-resolution Stepped-frequency SAR Imaging Technology

    Corresponding author: DING Zegang, z.ding@bit.edu.cn; LONG Teng, longteng@bit.edu.cn
  • Fund Project: The National Science Foundation for Distinguished Young Scholars (61625103), The National Natural Science Foundation of China (11833001)

    CLC number: TN 958

  • 摘要: 星载合成孔径雷达(SAR)是一种2维高分辨率微波成像雷达。它通过发射大带宽信号实现距离向高分辨,通过合成孔径技术实现方位向高分辨。随着人们对分辨率需求的不断提升,星载SAR正朝着分米级分辨率发展。一方面,受限于现有器件水平,可以通过频率步进技术实现大带宽信号发射,需要研究高精度子带拼接技术、子带间幅相误差对成像的影响与补偿技术;另一方面,受限于有限的波束宽度,可以使系统工作在聚束模式或滑聚模式实现长合成孔径,此时需研究轨道弯曲、“Stop-go”假设误差、电离层与对流层传输误差等非理想因素对成像的影响与补偿技术。因此,该文详细介绍了频率步进信号时序设计与子带拼接,研究星载高分辨率频率步进SAR成像算法与非理想因素补偿方法,最后给出成像算法的仿真验证和性能分析。
  • 图 1  子带并发

    Figure 1.  Simultaneously transmitted subpulses

    图 2  脉内子带串发

    Figure 2.  Sequential subpulses transmitted in PRT

    图 3  脉间子带串发

    Figure 3.  Sequential subpulses transmitted interpulse

    图 4  高分辨率星载SAR子带交替串发设计结果

    Figure 4.  Design result of subpulses transmitted alternatively for high resolution spaceborne SAR

    图 5  轨道弯曲引起的等效速度空变

    Figure 5.  Equivalent speed variance caused by curved orbit

    图 6  “Stop-go”误差引起的方位向偏移

    Figure 6.  Azimuth migration caused by stop-go error

    图 7  星载高分辨率SAR 100 s工作期间典型折射率下斜距误差的变化图

    Figure 7.  Range error variance of spaceborne high resolution SAR in typical refractive index during 100 s operation

    图 8  电离层导致的二次相位误差

    Figure 8.  Quadratic phase error caused by ionosphere error

    图 9  高分辨率滑动聚束星载SAR成像算法总流程

    Figure 9.  Flow chart of high resolution sliding spotlight spaceborne SAR imaging algorithm

    图 10  信号子带拼接前后对比

    Figure 10.  Stepped frequency subband signal and simulation result of synthetic bandwidth

    图 11  加入幅相误差的子带拼接结果

    Figure 11.  Compressed synthetic signal with amplitude errors and phase errors

    图 12  加入对流层误差的成像结果

    Figure 12.  Imaging result with troposphere error

    图 13  加入不同电离层误差的成像结果

    Figure 13.  Imaging results with different ionosphere errors

    图 14  加入不同大小DEM误差的成像结果

    Figure 14.  Imaging results with different DEM errors

    图 15  高分辨率滑动聚束星载SAR点阵目标成像结果图

    Figure 15.  Imaging result of dot targets with high resolution sliding spotlight spaceborne SAR

    表 1  高分辨率宽测绘带滑动聚束典型参数

    Table 1.  Typical parameters for sliding spotlight SAR with high resolution and wide swath

    参数数值
    距离向发射带宽2.4 G,单个脉冲400 M
    PRF3000~4000 Hz,分段变重频
    采样延迟分段变采样延迟
    聚束因子1/40
    下载: 导出CSV

    表 2  典型大气参数

    Table 2.  Typical atmosphere parameters

    参数数值
    场景大气压1.014 MPa
    相对湿度20%
    场景气温10°C
    场景海拔100 m
    下载: 导出CSV

    表 3  高分辨滑动聚束成像评估结果

    Table 3.  Evaluation results of high resolution sliding spotlight imaging

    点目标距离向方位向距离向分辨率 (m)方位向分辨率 (m)
    PSLR (dB)ISLR (dB)PSLR (dB)ISLR (dB)
    a–25.72–19.46–23.13–18.220.150.16
    b–25.42–19.46–22.97–19.010.150.16
    c–25.70–19.45–24.76–19.590.150.16
    d–25.66–19.49–21.77–18.690.150.16
    e–25.37–19.23–21.27–18.430.150.16
    下载: 导出CSV
  • [1] CUMMING I G and WONG F H. Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation[M]. Norwood: Artech House, 2005: 1–9.
    [2] 张澄波. 综合孔径雷达[M]. 北京: 科学出版社, 1989: 1–5.ZHANG Chengbo. Synthetic Aperture Radar[M]. Beijing: Science Press, 1989: 1–5.
    [3] 袁孝康. 星载合成孔径雷达导论[M]. 北京: 国防工业出版社, 2003: 1–5.YUAN Xiaokang. Introduce to the Spaceborne Synthetic Aperture Radar[M]. Beijing: National Defend Industry Press, 2003: 1–5.
    [4] WILEY C A. Synthetic aperture radars[J]. IEEE Transactions on Aerospace and Electronic Systems, 1985, AES-21(3): 440–443. doi: 10.1109/TAES.1985.310578
    [5] TSUNODA S I, PACE F, STENCE J, et al. Lynx: A high-resolution synthetic aperture radar[C]. SPIE 3704, Radar Sensor Technology IV, Orlando, USA, 1999: 1–4. doi: 10.1117/12.354602.
    [6] ENDER J H G and BRENNER A R. PAMIR-a wideband phased array SAR/MTI system[J]. IEE Proceedings-Radar, Sonar and Navigation, 2003, 150(3): 165–172. doi: 10.1049/ip-rsn:20030445
    [7] WERNINGHAUS R and BUCKREUSS S. The TerraSAR-X mission and system design[J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(2): 606–614. doi: 10.1109/tgrs.2009.2031062
    [8] 邓云凯, 陈倩, 祁海明, 等. 一种基于频域子带合成的多发多收高分辨率SAR成像算法[J]. 电子与信息学报, 2011, 33(5): 1082–1087. doi: 10.3724/SP.J.1146.2010.01067DENG Yunkai, CHEN Qian, QI Haiming, et al. A high-resolution imaging algorithm for MIMO SAR based on the sub-band synthesis in frequency domain[J]. Journal of Electronics &Information Technology, 2011, 33(5): 1082–1087. doi: 10.3724/SP.J.1146.2010.01067
    [9] LORD R T and INGGS M R. High resolution SAR processing using stepped-frequencies[C]. 1997 IEEE International Geoscience and Remote Sensing Symposium, Remote Sensing - A Scientific Vision for Sustainable Development, Singapore, 1997: 490–492. doi: 10.1109/IGARSS.1997.615924.
    [10] WU Yuan, SUN Guangcai, YANG Chun, et al. Processing of very high resolution spaceborne sliding spotlight SAR data using velocity scaling[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(3): 1505–1518. doi: 10.1109/TGRS.2015.2481923
    [11] PRATS-IRAOLA P, SCHEIBER R, RODRIGUEZ-CASSOLA M, et al. On the processing of very high resolution spaceborne SAR data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(10): 6003–6016. doi: 10.1109/TGRS.2013.2294353
    [12] BELCHER D P. Theoretical limits on SAR imposed by the ionosphere[J]. IET Radar, Sonar & Navigation, 2008, 2(6): 435–448. doi: 10.1049/iet-rsn:20070188
    [13] 王沛. 高分辨率星载合成孔径雷达关键技术研究与验证[D]. [博士论文], 中国科学院大学, 2018: 47–49.WANG Pei. Study on key technology and verification of spaceborne high resolution synthetic aperture radar[D]. [Ph.D. dissertation], University of Chinese Academy of Sciences, 2018: 47–49.
    [14] VILLANO M, KRIEGER G, and MOREIRA A. Staggered-SAR for high-resolution wide-swath imaging[C]. 2012 IET International Conference on Radar Systems, Glasgow, UK, 2012: 1–6. doi: 10.1049/cp.2012.1600.
    [15] 秦显平. 星载GPS低轨卫星定轨理论及方法研究[D]. [博士论文], 解放军信息工程大学, 2009: 1–3.QIN Xianping. Research on precision orbit determination theory and method of low earth orbiter based on GPS technique[D]. [Ph.D. dissertation], The PLA Information Engineering University, 2009: 1–3.
    [16] ULANDER L M H, HELLSTEN H, and STENSTROM G. Synthetic-aperture radar processing using fast factorized back-projection[J]. IEEE Transactions on Aerospace and Electronic Systems, 2003, 39(3): 760–776. doi: 10.1109/TAES.2003.1238734
    [17] DING Zegang, GUO Yansu, GAO Wenbin, et al. A range grating lobes suppression method for stepped-frequency SAR imagery[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(12): 5677–5687. doi: 10.1109/JSTARS.2016.2593711
    [18] DING Zegang, GAO Wenbin, LIU Jingyun, et al. A novel range grating lobe suppression method based on the stepped-frequency SAR image[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(3): 606–610. doi: 10.1109/LGRS.2014.2352676
  • [1] 唐江文邓云凯王宇赵硕李宁 . 高分辨率滑动聚束SAR BP成像及其异构并行实现. 雷达学报, 2017, 6(4): 368-375. doi: 10.12000/JR16053
    [2] 折小强仇晓兰韩冰雷斌 . 一种基于变换域的滑动聚束SAR 调频率估计方法. 雷达学报, 2014, 3(4): 419-427. doi: 10.3724/SP.J.1300.2014.14008
    [3] 马兵强 . 滑动聚束FMCW-SAR的子孔径波数域成像算法. 雷达学报, 2013, 2(3): 319-325. doi: 10.3724/SP.J.1300.2013.13044
    [4] 韦维朱岱寅吴迪 . 基于尺度变换原理的SAR波数域成像算法. 雷达学报, 2020, 9(2): 354-362. doi: 10.12000/JR19112
    [5] 郑建成王党卫马晓岩宣泽平丰孝兵 . 基于自旋运动的高速弹头成像方法. 雷达学报, 2013, 2(3): 300-308. doi: 10.3724/SP.J.1300.2013.13070
    [6] 曾操梁思嘉王威徐青 . 基于频率步进信号的旋转式合成孔径雷达成像方法. 雷达学报, 2014, 3(4): 401-408. doi: 10.3724/SP.J.1300.2014.14043
    [7] 向寅张凯胡程 . 基于NUFFT的调频步进频高分辨成像与目标识别算法. 雷达学报, 2015, 4(6): 639-647. doi: 10.12000/JR15083
    [8] 杨利超高悦欣邢孟道盛佳恋 . 基于广义keystone和频率变标的微波光子ISAR高分辨实时成像算法. 雷达学报, 2019, 8(2): 215-223. doi: 10.12000/JR18120
    [9] 吴玉峰叶少华冯大政 . 基于方位相位编码的脉内聚束SAR成像方法. 雷达学报, 2018, 7(4): 437-445. doi: 10.12000/JR17114
    [10] 陈世阳黄丽佳俞雷 . 基于改进sinc插值的变PRF采样聚束SAR成像. 雷达学报, 2019, 8(4): 527-536. doi: 10.12000/JR18095
    [11] 惠叶白雪茹 . 基于RID序列的微动目标高分辨三维成像方法. 雷达学报, 2018, 7(5): 548-556. doi: 10.12000/JR18056
    [12] 廖明生魏恋欢汪紫芸TimoBalz张路 . 压缩感知在城区高分辨率SAR层析成像中的应用. 雷达学报, 2015, 4(2): 123-129. doi: 10.12000/JR15031
    [13] 李宁汪玲张弓 . 多基ISAR 舰船侧视及俯视高分辨率成像方法. 雷达学报, 2012, 1(2): 163-170. doi: 10.3724/SP.J.1300.2012.20021
    [14] 陈潇翔邢孟道 . 基于空变运动误差分析的微波光子超高分辨SAR成像方法. 雷达学报, 2019, 8(2): 205-214. doi: 10.12000/JR18121
    [15] 李春升于泽陈杰 . 高分辨率星载SAR成像与图像质量提升方法综述. 雷达学报, 2019, 8(6): 717-731. doi: 10.12000/JR19085
    [16] 张文彬邓云凯王宇 . 星地双基合成孔径雷达聚束模式快速BP算法. 雷达学报, 2013, 2(3): 357-366. doi: 10.3724/SP.J.1300.2013.13031
    [17] 罗绣莲徐伟郭磊 . 捷变PRF技术在斜视聚束SAR中的应用. 雷达学报, 2015, 4(1): 70-77. doi: 10.12000/JR14149
    [18] 闫敏韦顺军田博坤张晓玲师君 . 基于稀疏贝叶斯正则化的阵列SAR高分辨三维成像算法. 雷达学报, 2018, 7(6): 705-716. doi: 10.12000/JR18067
    [19] 赵团邓云凯王宇李宁王翔宇 . 基于扇贝效应校正的改进滑动Mosaic全孔径成像算法. 雷达学报, 2016, 5(5): 548-557. doi: 10.12000/JR16014
    [20] 周超伟李真芳王跃锟解金卫 . 联合多方位角调频率估计的星载SAR三维成像方法. 雷达学报, 2018, 7(6): 696-704. doi: 10.12000/JR18094
  • 加载中
图(15)表(3)
计量
  • 文章访问数:  1016
  • HTML浏览量:  586
  • PDF下载量:  241
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-08-30
  • 录用日期:  2019-12-01
  • 网络出版日期:  2019-12-30
  • 刊出日期:  2019-12-28

星载高分辨频率步进SAR成像技术

    通讯作者: 丁泽刚 z.ding@bit.edu.cn; 龙腾 longteng@bit.edu.cn
    作者简介:
    龙 腾(1968–),男,湖北黄陂人。北京理工大学教授、博士生导师。现任北京理工大学副校长、雷达技术研究所所长。研究方向为新体制雷达与实时信息处理。E-mail: longteng@bit.edu.cn;
    丁泽刚(1980–),男,河南信阳人,博士。现担任北京理工大学研究员、博士生导师,主要研究方向为新体制雷达成像机理、成像处理和图像信息提取。E-mail: z.ding@bit.edu.cn;
    肖 枫(1993–),男,内蒙古包头人。北京理工大学雷达技术研究所在读博士研究生,主要研究领域为多模式星载SAR成像技术、SAR运动补偿技术。E-mail: 1120646612@qq.com;
    王 岩(1989–),男,河北沧州人。博士后,现担任北京理工大学副研究员、硕士生导师,主要研究方向为新体制雷达系统、成像、干涉和极化应用。E-mail: yan_wang@bit.edu.cn;
    李 喆(1996–),男,河北任丘人。北京理工大学信息与电子学院硕士研究生,主要研究方向星载合成孔径雷达系统设计与成像。E-mail: 3120180770@bit.edu.cn
  • 北京理工大学信息与电子学院雷达技术研究所 北京 100081
  • 北京理工大学嵌入式实时信息处理技术北京市重点实验室 北京 100081
基金项目:  国家杰出青年科学基金(61625103),国家自然科学基金重点项目(11833001)

摘要: 星载合成孔径雷达(SAR)是一种2维高分辨率微波成像雷达。它通过发射大带宽信号实现距离向高分辨,通过合成孔径技术实现方位向高分辨。随着人们对分辨率需求的不断提升,星载SAR正朝着分米级分辨率发展。一方面,受限于现有器件水平,可以通过频率步进技术实现大带宽信号发射,需要研究高精度子带拼接技术、子带间幅相误差对成像的影响与补偿技术;另一方面,受限于有限的波束宽度,可以使系统工作在聚束模式或滑聚模式实现长合成孔径,此时需研究轨道弯曲、“Stop-go”假设误差、电离层与对流层传输误差等非理想因素对成像的影响与补偿技术。因此,该文详细介绍了频率步进信号时序设计与子带拼接,研究星载高分辨率频率步进SAR成像算法与非理想因素补偿方法,最后给出成像算法的仿真验证和性能分析。

English Abstract

    • 合成孔径雷达[1-3]是一种2维高分辨率微波成像雷达,它利用脉冲压缩技术实现距离向高分辨,利用合成孔径技术实现方位向高分辨。自1951年合成孔径的概念被Wiley[4]提出以来,合成孔径雷达一直在向着高分辨、宽测绘带的方向发展。目前典型的机载SAR系统已实现分米级甚至厘米级分辨率成像,其中美国Sandia实验室和General Atomics公司联合开发的Lynx SAR可实现0.1 m分辨率成像[5],德国Frounhofer高频物理与雷达技术研究所研制的PAMIR可实现0.05 m分辨率成像[6]。星载SAR系统目前也发展到了分米级分辨率成像的水平,其中美国FIA系统实现了优于0.3 m分辨率成像,德国TerraSAR-X实现了方位向0.16 m分辨率成像[7]

      高分辨率能增强SAR图像的地物表征能力,但也会给星载SAR的系统设计和信号处理带来新的问题。一方面,星载SAR的距离向分辨率取决于雷达信号带宽。受限于现有的器件水平,过大的信号带宽会给系统的发射、接收和处理带来极大困难。为降低硬件压力,可以使用频率步进技术,发射多个不同载频的子带信号,再使用子带拼接技术将子带回波合成为宽带回波[8,9]。目前世界上很多高分辨系统都采用了频率步进体制,以PAMIR系统为例,该系统最初能通过5个子带实现1.8 GHz的信号带宽[6],在2011年系统升级后其带宽进一步提升到3.6 GHz。然而在实际工程中,频率步进体制将带来脉冲重复频率增加、子带间能量泄露与子带间幅相误差等问题。这些问题会导致距离向测绘带宽度减小、脉冲压缩质量下降。因此,频率步进信号的时序设计以及精确的子带拼接方法是实现距离高分辨的关键技术。

      另一方面,星载SAR的方位向分辨率取决于合成孔径时间。采用聚束模式或滑聚模式能够在波束宽度受限的情况下,通过波束控制的方式增大合成孔径时间,从而突破条带模式方位分辨率的限制。然而,由于星载SAR轨道弯曲严重、电波传输环境复杂,长合成孔径时间成像会受到诸多非理想因素的影响,例如成像参数空变[10]、“Stop-go”模型误差[11]、对流层传输延迟[11]与电离层导致的色散[12]。这些非理想因素大多存在时空变化的特性,它们会在信号中引入相位误差并在合成孔径时间内积累,最终导致距离徙动校正失败或图像散焦。因此,非理想因素的分析和补偿也是实现星载SAR高分辨成像所必须解决的问题。

      针对上述问题,本文针对星载高分辨频率步进SAR的特点,介绍了频率步进信号的时序设计方法和子带拼接原理,给出了一种高分辨率成像算法与非理想因素补偿方法,最后进行了算法仿真和性能分析验证。本文结构如下:第2节介绍频率步进信号时序设计和子带拼接技术原理;第3节给出了星载高分辨率频率步进SAR成像算法与非理想因素补偿方法;第4节给出了成像仿真和性能分析验证;第5节对全文进行总结。

    • 频率步进技术能够实现大带宽信号发射,但同时也会引入脉冲重复频率增加、子带间能量泄露与子带间幅相误差等问题,而合理的信号时序设计能够在一定程度上减轻上述问题的影响。根据发射和接收信号的子带时序关系,频率步进SAR系统可分为3种基本工作模式:子带并发、脉内子带串发和脉间子带串发[13]

      (1) 子带并发模式

      各子带被调制到不同载频上,在同一时刻通过不同天线孔径发射出去,接收时使用N个独立接收通道接收对应的N个子带并进行解调处理,从而获得各子带的基带信号。子带并发模式的时序关系如图1所示,其中Tp为脉冲宽度,Bi为单子带带宽,B为总带宽。这种模式的优点是能节省发射脉冲的时间,不需要因为子带分割发射而增加系统脉冲重复频率,不影响系统时序设计。缺点是多子带同时发射或接收会导致射频兼容性问题以及能量泄露问题,进而影响子带拼接效果。

      图  1  子带并发

      Figure 1.  Simultaneously transmitted subpulses

      (2) 脉内子带串发模式

      各子带信号以一个大于脉冲宽度的时间间隔在一个PRT内顺序发射,接收时采用N子带接收机同时接收。脉内子带串发模式的工作时序如图2所示,其中ΔT为子带间的发射时间间隔。这种模式由于在发射子带间设置了保护时间,因此可以避免发射子带之间的射频兼容性问题,但由于各子带回波仍然只能通过滤波的方式在频域区分,接收时子带间的频谱泄漏问题仍然存在。

      图  2  脉内子带串发

      Figure 2.  Sequential subpulses transmitted in PRT

      (3) 脉间子带串发模式

      该模式每个PRT仅发射/接收1个脉冲,其工作时序如图3所示。由于发射/接收子带在时间上完全分离,因此可以完全避免子带之间发射和接收的干扰。然而,脉间串发方式会增加系统所需PRF,进而导致距离向测绘带宽度减小。

      图  3  脉间子带串发

      Figure 3.  Sequential subpulses transmitted interpulse

      对于子带数目较多的系统,上述3种方法难以在降低子带间干扰的同时保证测绘带宽度。此时可以将脉内子带串发和脉间子带串发结合,即使用交替串发模式。具体来说,可将所有子带分为两个或多个组合,组合内子带频谱互不相邻,可采用脉内子带串发的方式发射,不同组合则采用脉间串发的方式发射。通过这种方式,交替串发模式能够在避免子带间干扰的同时维持可接受的距离向带宽。

      以2.4 GHz的总带宽为例,若单子带带宽为400 MHz,则一共需要6个子带。若采用脉间串发模式,系统的PRF会增大为原来的6倍,测绘带宽度将降低为原来的1/6。若采用交替串发模式,将各子带按照中心频率的高低编号为1–6,可如图4所示,在奇数PRT中发射子带1, 3和5,在偶数PRT中发射子带2, 4和6。由于每个PRT内仅有3个互不相邻的子带,因此仅需设置3个通道进行回波接收和AD采样,且3个通道间不存在互相干扰的问题。此外,由于PRF仅增加了1倍,距离向测绘带宽仍然可以接受。

      图  4  高分辨率星载SAR子带交替串发设计结果

      Figure 4.  Design result of subpulses transmitted alternatively for high resolution spaceborne SAR

    • 频率步进信号的子带拼接处理可分为频域拼接[8]和时域拼接[9]两种方法。这两种方法都能够实现良好的拼接效果,因此可以根据后续处理算法选择合适的方法。由于频域拼接原理直观,因此论文采用频域拼接方法。

      频域拼接方法的思路是,利用距离向信号脉压后相位为常数的特点,先将各子带信号进行脉压,然后变换到频域进行频谱拼接。在频域合成一个宽带信号脉压频谱后,再对频谱拼接处的相位跳变进行补偿,最终可得到宽带信号的脉压结果。

      设第n个频点解调后的回波信号为

      $ \begin{split} {S_n}\left( t \right) =\, & {\rm{rect}}\left( {\frac{{{t_1} - {t_0} - \dfrac{{{T_{\rm{p}}}}}{2} - {t_{\rm{d}}}\left( n \right)}}{{{T_{\rm{p}}}}}} \right)\\ & \cdot \exp \left( {{\rm{j}}{\rm{\pi }}{K_{\rm{r}}}{{\left( {{t_1} - {t_0} - \frac{{{T_{\rm{p}}}}}{2} - {t_{\rm{d}}}\left( n \right)} \right)}^2}} \right)\\ & \cdot\exp \left( { - {\rm{j}}2{\rm{\pi }}{f_{\rm{c}}}\left( n \right)\frac{{2R}}{\rm c}} \right) \end{split} $

      其中,t1为距离向时间轴,Tp是脉冲宽度,fc(n)是第n个频点的载频,Kr是距离向调频率,R是目标斜距,t0=2(RRmin)/c, Rmin是起始采样距离,c是光速,td(n)是第n个子带与第1个子带间的延迟。

      对回波信号进行傅里叶变换,得到子带信号频谱Sn(f )

      $ \begin{split} {S_n}\left( f \right) =\,& {\rm{rect}}\left( {\frac{f}{B}} \right)\exp \left( { - {\rm{j}}{\rm{\pi }}\frac{{{f^2}}}{{{K_{\rm{r}}}}}} \right)\\ & \cdot\exp \left( { - {\rm{j}}2{\rm{\pi }}f\left( {{t_0}{\rm{ + }}\frac{{{{\rm{T}}_{\rm{p}}}}}{2} + {t_{\rm{d}}}\left( n \right)} \right)} \right)\\ & \cdot \exp \left( { - {\rm{j}}2{\rm{\pi }}{f_{\rm{c}}}\left( n \right)\frac{{2R}}{\rm c}} \right) \end{split} $

      为进行脉冲压缩与时域对齐,将子带回波频谱与式(3)相乘

      $ \begin{split} {H_{{\rm{rc}}}}\left( f \right)=\,&{\rm{ rect}}\left( {\frac{f}{{\Delta f}}} \right)\exp \left( {{\rm{j}}{\rm{\pi }}\frac{{{f^2}}}{{{K_{\rm{r}}}}}} \right)\\ & \cdot\exp \left( {{\rm{j}}2{\rm{\pi }}f\left( {\frac{{{T_{\rm{p}}}}}{2} + {t_{\rm{d}}}\left( n \right)} \right)} \right) \end{split} $

      得到脉压后子带频谱

      $ {S_n}\left( f \right) = {\rm{rect}}\!\left(\! {\frac{f}{B}}\! \right)\exp \left(\! { - {\rm{j}}2{\rm{\pi }}f{t_0}} \right)\exp \!\left( { - {\rm{j}}2{\rm{\pi }}{f_{\rm{c}}}\left( n \right)\frac{{2R}}{\rm c}} \!\right) $

      为了适应合成后的信号带宽,脉压后需对信号进行升采样。升采样倍数与子带数相等,升采样后信号形式不变。

      f0为第1个子带的载频,对于第n个频点,其载频为

      $ {f_{\rm c}}\left( n \right) = {f_0} + \left( {n - 1} \right)\Delta f $

      为了恢复各子带的频率步进关系,在进行频谱叠加之前需进行频移操作,各子带的频率偏移量fshift(n)为

      ${f_{{\rm{shift}}}}\left( n \right) = \left( {n - \frac{{N + 1}}{2}} \right)\Delta f$

      其中,N为子带数量。

      频移后第n个频点的频域表达式为

      $ \begin{split} {S_n}\left( f \right) =\,& {\rm{rect}}\left( {\frac{{f - {f_{{\rm{shift}}}}\left( n \right)}}{B}} \right)\\ & \cdot\exp \Bigr\{ { - {\rm{j}}2{\rm{\pi }}\Bigr[ {f - {f_{{\rm{shift}}}}\left( n \right)} \Bigr]{t_0}} \Bigr\}\\ &\cdot \exp \left\{ { - {\rm{j}}2{\rm{\pi }}{f_{\rm{c}}}\left( n \right)\frac{{2R}}{\rm c}} \right\} \end{split} $

      将式(5)和式(6)代入式(7)得

      $ \begin{split} {S_n}\left( f \right) =\,& {\rm{rect}}\left( {\frac{{f - \left( {n - \dfrac{{N + 1}}{2}} \right) \cdot \Delta f}}{B}} \right)\\ &\cdot \exp \left\{ { - {\rm{j}}2{\rm{\pi }}f\frac{{2\left( {R - {R_{\min }}} \right)}}{\rm c}} \right\}\\ & \cdot \exp \left\{ { - {\rm{j}}2{\rm{\pi }}\left( {{f_0} + \frac{{N - 1}}{2} \cdot \Delta f} \right)\frac{{2R}}{\rm c}} \right\} \\ & \cdot \exp \left\{ { - {\rm{j}}2{\rm{\pi }}\left( {n - \frac{{N + 1}}{2}} \right) \cdot \Delta f \cdot \frac{{2{R_{\min }}}}{\rm c}} \right\} \\ \end{split}$

      式(8)最后1个指数项表明,频移后各子带频谱相位之间仍存在跳变。这会导致拼接信号脉压结果出现栅瓣。因此,为保证子带拼接时相位的连续性,需要对跳变相位进行补偿。补偿项可以表示为

      $ \begin{split} & {\varphi _{{\rm{comp}}}}\left( n \right) = \exp \left\{ {{\rm{j}}2{\rm{\pi }}\left( {n - \frac{{N + 1}}{2}} \right) \cdot \Delta f \cdot \frac{{2{R_{\min }}}}{\rm c}} \right\},\\ & \quad\qquad n = 1,2,··· ,N\\[-10pt] \end{split} $

      补偿后的频域表达式为

      $ \begin{split} {S_n}\left( f \right) =\,& {\rm{rect}}\left( {\frac{{f - \left( {n - \dfrac{{N + 1}}{2}} \right) \cdot \Delta f}}{B}} \right)\\ &\cdot \exp \left\{ { - {\rm{j}}2{\rm{\pi }}f\frac{{2\left( {R - {R_{\min }}} \right)}}{\rm c}} \right\}\\ &\cdot \exp \left\{ { - {\rm{j}}2{\rm{\pi }}\left( {{f_0} + \frac{{N - 1}}{2} \cdot \Delta f} \right)\frac{{2R}}{\rm c}} \right\}\!\! \end{split} $

      将补偿后的各子带信号在频域叠加求和,即可得到带宽为N·B的宽带信号

      $ \begin{split} S\left( f \right) =\,& \sum\limits_{n = 1}^N {{S_n}\left( f \right) = } {\rm{rect}}\left( {\frac{f}{{NB}}} \right) \cdot \exp \left\{ { - {\rm{j}}2{\rm{\pi }}f{t_0}} \right\} \\ &\cdot \exp \left\{ { - {\rm{j}}2{\rm{\pi }}\left( {{f_0} + \frac{{N - 1}}{2} \cdot \Delta f} \right)\frac{{2R}}{\rm c}} \right\}\\[-17pt] \end{split} $

      需要说明的是,上述分析中假设不同子带的斜距是相同的。事实上,脉间串发时卫星的位置变化不可忽略,由此导致的斜距误差需要在子带拼接时进行补偿。另外,子带间幅相误差也会对子带拼接造成影响,相关分析将在第4节给出。

    • 星载SAR方位向分辨率由合成孔径时间决定。在波束宽度受限的情况下,传统条带模式难以满足高分辨成像需求,因此系统需要采用滑动聚束模式。将滑动聚束SAR波足速度与卫星速度之比定义为聚束因子A,那么滑动聚束SAR的分辨率可表示为

      ${\rho _{\rm{a}}} = \frac{{{D_{\rm{a}}}}}{2}A$

      其中,Da为方位向天线尺寸。

      以前面论证的子带交替串发系统为例,其距离向信号带宽为2.4 GHz。设方位向天线尺寸为10 m,那么若要在方位向实现与距离向匹配的高分辨率,聚束因子需设计在1/40左右,此时单点多普勒带宽为1400 Hz。考虑星下点回波时隙保护约束、发射信号时隙保护约束以及交替串发带来的PRF提升,可确定PRF范围为3000~4000 Hz。

      在高分辨滑动聚束模式下,成像幅宽大,斜距历程变化大。如果采用固定PRF和固定采样延迟,会大大增加回波窗口长度,进而导致波位设计失败。这一问题需要通过变重频[14]方法解决,即在整个孔径过程中多次变换PRF和起始采样时刻,以适应回波延迟的变化。为使距离向点数对齐,系统接收回波时长为固定值。为确保回波接收完整,一般取整个合成孔径内的最大值为系统接收回波时长。

      综上,高分辨率滑动聚束典型参数总结如表1所示。

      参数数值
      距离向发射带宽2.4 G,单个脉冲400 M
      PRF3000~4000 Hz,分段变重频
      采样延迟分段变采样延迟
      聚束因子1/40

      表 1  高分辨率宽测绘带滑动聚束典型参数

      Table 1.  Typical parameters for sliding spotlight SAR with high resolution and wide swath

    • 由于星载SAR轨道弯曲严重、电波传输环境复杂,因此长合成孔径时间成像时会受到诸多非理想因素的影响,例如轨道弯曲、“Stop-go”假设误差、对流层误差和电离层误差。

      (1) 轨道弯曲

      由于星载SAR几何关系的复杂性,星载SAR等效速度随方位空变、随高程空变的现象不能忽略,如图5所示。传统基于双曲线斜距模型的成像方法由于没有考虑等效速度的方位空变性,此时不再适用。文献[9]把超高分辨率情况下轨道弯曲导致的误差视为运动误差,并提出一种基于运动补偿的成像方法。另外,考虑到卫星上GPS定轨精度可达到厘米级[15],在构建投影网格时也能将DEM考虑在内,因此采用时域和其快速算法成像也是一种有效的处理方法。

      图  5  轨道弯曲引起的等效速度空变

      Figure 5.  Equivalent speed variance caused by curved orbit

      (2) “Stop-go”假设误差[11]

      在传统SAR成像算法中,通常假设平台在发射脉冲到接收回波期间是静止的,这种假设被称为“Stop-go”假设。当星载SAR分辨率提高到分米级时,“Stop-go”假设引入的误差不可忽略。

      在慢时间方面,平台发射脉冲和接收脉冲的位置不同,因此目标在图像中会产生方位向偏移。由于不同距离门对应的延迟不同,上述方位向偏移会随着目标斜距而变化,如图6所示。“Stop-go”误差可以在距离多普勒域中通过式(13)进行校正

      图  6  “Stop-go”误差引起的方位向偏移

      Figure 6.  Azimuth migration caused by stop-go error

      $ {H_{{\rm{slow}}}}\left( {{f_{\rm{a}}},r} \right) = \exp \left( {{\rm{j}} \cdot 2{\rm{\pi }} \cdot \frac{r}{\rm c} \cdot {f_{\rm{a}}}} \right) $

      其中,r为距离轴,fa为方位向频率轴。

      在快时间方面,由于发射脉冲期间平台的运动,脉冲本身会产生频偏,频偏量等于瞬时多普勒频率。当多普勒带宽较大时,这一偏移会导致距离徙动校正失败,因此需要在2维频域进行校正。

      (3) 对流层误差

      大气层中的对流层是影响电磁波传输的主要因素,它处于大气层的低层,高度约为10~14 km,电磁波在其中传输会产生明显的折射现象。对米级分辨率星载SAR来说,对流层折射现象对SAR成像的影响可以忽略。然而当雷达分辨率达到分米级时,该折射现象将引起合成孔径时间内的斜距误差,进而造成距离徙动校正失败和方位向散焦。在表2所示的典型大气参数下,当雷达工作时间为100 s时,电磁波传输路径延迟变化如图7所示。由图7可见,场景中心点在10 s合成孔径时间内斜距误差0.07 m,场景边缘点斜距误差在合成孔径时间内达到1.6 m,远超过1个距离单元,因此这样的斜距误差必须补偿。若方位向波束宽度内目标的斜距误差差异可以忽略,可以在回波域通过乘以相位exp(j4πΔR(ta)/λ)补偿,其中ΔR(ta)为随时间变化的斜距误差。

      参数数值
      场景大气压1.014 MPa
      相对湿度20%
      场景气温10°C
      场景海拔100 m

      表 2  典型大气参数

      Table 2.  Typical atmosphere parameters

      图  7  星载高分辨率SAR 100 s工作期间典型折射率下斜距误差的变化图

      Figure 7.  Range error variance of spaceborne high resolution SAR in typical refractive index during 100 s operation

      (4) 电离层误差

      地球的电离层高度约为60~2000 km,这也是现有所有星载SAR的运行轨道高度范围,因此星载SAR的发射和接收回波均要经过电离层。对于分辨率较低的星载SAR,电离层延迟引起的色散效应非常小,所以很少考虑其对成像的影响。然而,当分辨率达到分米级,信号带宽达到2.4 GHz时,电离层的散射效应对信号脉压的影响不能再忽略。根据卡普曼(Chapman)模型,电磁波在电离层中传播的延迟Δt可表示为

      $\Delta t = \frac{K}{{{\rm{c}}{f^2}\cos \alpha }} \cdot {\rm{TEC}}$

      其中,K为常数,c为光速,f为信号频率,α为雷达下视角,TEC为电离层的积分电子含量。

      电磁波传输延迟随频率的变化表明,可将电离层理解为一种色散介质,根据信号传播理论并考虑星载SAR信号的双程传播现象,由电离层色散效应导致的距离频域相位误差$ \Delta \varphi $

      $ \Delta \varphi = - 4{\rm{\pi }}\frac{K}{{{\rm{c}}f\cos \alpha }} \cdot {\rm{TEC}} $

      由式(15)可以看到,电离层误差会在距离频域中引入随频率变化的相位误差,将相位误差进行泰勒展开,其二次项及高次项会导致图像散焦。表1所示参数对应的二次误差相位随TEC变化曲线如图8所示。由图8可见在TEC超过28 TECU时电离层带来的相位误差就会超过π/4,从而造成散焦。因此对于2.4 GHz带宽信号而言,当TEC达到28 TECU以上时,需要对电离层误差进行补偿。

      图  8  电离层导致的二次相位误差

      Figure 8.  Quadratic phase error caused by ionosphere error

      对于电离层误差,可以采用中国相关科研机构发布的TEC测量数据来补偿相关附加频域相位,还可以采用自聚焦的手段估计TEC,进而补偿电离层影响。

      考虑到利用频域算法解决等效速度空变问题过于复杂[10],论文拟采用基于BP算法[16]的星载高分辨频率步进SAR成像。算法流程图如图9所示,首先进行对流层误差的补偿;接着在进行电离层误差补偿的同时完成子带拼接;再采用插值方法将方位向分段变重频信号恢复成均匀信号;进而进行距离向脉压;最后根据外部DEM实现空变的地形精确BP成像。

      图  9  高分辨率滑动聚束星载SAR成像算法总流程

      Figure 9.  Flow chart of high resolution sliding spotlight spaceborne SAR imaging algorithm

    • 首先根据表1所示参数进行频率步进信号子带拼接的仿真。仿真采用子带交替串发模式,通过6个子带拼接得到2.4 GHz的总带宽,实现距离向高分辨率成像。子带信号与合成后宽带信号的频谱如图10(a)所示,可见单子带带宽为400 MHz而6子带合成后带宽达到2.4 GHz;子带信号脉压结果和合成后宽带信号脉压结果如图10(b)所示,旁瓣抑制采用–25 dB泰勒窗,可见合成后信号主瓣变窄,分辨率提升。

      图  10  信号子带拼接前后对比

      Figure 10.  Stepped frequency subband signal and simulation result of synthetic bandwidth

      需要说明的是,串发模式中各个子带收发位置并不一样,因此子带信号之间存在一定的斜距误差,且该误差在斜视角最大时(即刚开始工作时)达到最大。在表1所示参数下进行计算,可知此时沿轨道方向各个子带收发位置最大差异为1.9 m,对场景中目标斜距差为0.9 m,此时的斜距误差无法忽略。斜距误差引入的相位误差可以在各子带分别乘以相位exp(j4πΔR(ta)/λ)进行补偿,而包络误差则需要在距离频域乘以相应的时移因子进行校正。斜距误差主要导致的是子带间相位不连续,因此要保证子带拼接效果,斜距补偿精度需在次波长级。

      此外,如第2节所述,实际系统通道间的不一致性会导致子带间存在幅相误差。为便于分析,在6个子带中取出两个子带,单独加入0 dB, 10 dB, 20 dB的幅度误差和0 rad, π/4 rad, π/2 rad的相位误差,经过子带拼接和脉冲压缩后,结果如图11所示。由图11可见,子带间的幅度误差会造成主瓣展宽、旁瓣提升和波形变化,而相位误差会造成旁瓣不对称。因此,在实际多通道系统中需要结合定标或者一些幅相误差估计和栅瓣抑制手段[17,18]来避免此类问题。

      图  11  加入幅相误差的子带拼接结果

      Figure 11.  Compressed synthetic signal with amplitude errors and phase errors

    • 首先使用表1所示参数对典型参数下的对流层误差进行点目标仿真,成像结果如图12所示。由图12可见,对流层误差能够引起明显的方位向散焦。此外,随着场景的增大,方位向边缘目标的对流层误差空变性也会增大,当合成孔径时间内的斜距误差变化超过1个距离单元时,会导致距离徙动校正失败。因此,在高分辨成像过程中,必须采取适当措施对对流层引入的斜距误差进行补偿。

      图  12  加入对流层误差的成像结果

      Figure 12.  Imaging result with troposphere error

      为分析电离层对成像的影响,分别对电子含量为0 TECU, 20 TECU, 40 TECU和60 TECU的情况进行仿真,成像结果如图13所示。由图13可见,在未补偿电离层误差的情况下,信号色散效应明显,图像有严重的主瓣展宽和旁瓣升高。因此,为了实现良好的聚焦效果,需采用自聚焦等手段进行电离层相位补偿,且补偿精度优于10 TECU。一般来说,自聚焦算法精度与信号信噪比、图像内容等因素相关。对于信噪比为0 dB的点目标信号来说,最大对比度算法的补偿精度能达到3 TECU以内,但对于相同信噪比的真实场景来说,算法性能会有所下降。

      图  13  加入不同电离层误差的成像结果

      Figure 13.  Imaging results with different ionosphere errors

      在分米级分辨率下,BP算法的成像网格中需要包含DEM信息,而DEM的精度也会对成像质量产生影响。图14展示了高程误差0 m, 10 m, 20 m, 30 m时的成像结果,可见当存在DEM误差时,图像方位向会有散焦现象。要实现良好的聚焦效果,DEM的精度需达到10 m以上。

      图  14  加入不同大小DEM误差的成像结果

      Figure 14.  Imaging results with different DEM errors

      最后,进行高分辨率滑动聚束SAR的BP成像仿真。由于大场景数据量巨大,这里只仿真了2 km幅宽的场景,场景中设置3×3点阵,点目标两维间距为1 km。对仿真场景进行BP成像时,两维均加–25 dB的泰勒窗,成像结果如图15所示。

      图  15  高分辨率滑动聚束星载SAR点阵目标成像结果图

      Figure 15.  Imaging result of dot targets with high resolution sliding spotlight spaceborne SAR

      为评估该成像处理算法的有效性,取点阵中4个角点以及中心点对成像指标进行评估,结果如表3所示。由评估结果可见该成像算法能够满足高分辨率的成像要求,成像指标良好。

      点目标距离向方位向距离向分辨率 (m)方位向分辨率 (m)
      PSLR (dB)ISLR (dB)PSLR (dB)ISLR (dB)
      a–25.72–19.46–23.13–18.220.150.16
      b–25.42–19.46–22.97–19.010.150.16
      c–25.70–19.45–24.76–19.590.150.16
      d–25.66–19.49–21.77–18.690.150.16
      e–25.37–19.23–21.27–18.430.150.16

      表 3  高分辨滑动聚束成像评估结果

      Table 3.  Evaluation results of high resolution sliding spotlight imaging

    • 本文针对星载高分辨频率步进SAR成像的特点,分析了频率步进信号时序设计方法和子带拼接原理,给出了一种基于BP算法的高分辨星载SAR成像处理流程,分析了轨道弯曲、“Stop-go”假设误差、对流层误差、电离层误差和地形起伏等非理想因素的影响和补偿方法,并通过计算机仿真对所提方法进行了仿真验证。相关分析有助于高分辨率星载SAR成像处理和工程实现。

参考文献 (18)

目录

    /

    返回文章
    返回