基于辅助数字高程模型的方位多通道SAR相位失配校正方法

李强 范怀涛

李强, 范怀涛. 基于辅助数字高程模型的方位多通道SAR相位失配校正方法[J]. 雷达学报, 2019, 8(5): 616–623. doi: 10.12000/JR19009
引用本文: 李强, 范怀涛. 基于辅助数字高程模型的方位多通道SAR相位失配校正方法[J]. 雷达学报, 2019, 8(5): 616–623. doi: 10.12000/JR19009
LI Qiang and FAN Huaitao. Channel phase mismatch calibration for multichannel in azimuth SAR imaging aided by digital elevation model[J]. Journal of Radars, 2019, 8(5): 616–623. doi: 10.12000/JR19009
Citation: LI Qiang and FAN Huaitao. Channel phase mismatch calibration for multichannel in azimuth SAR imaging aided by digital elevation model[J]. Journal of Radars, 2019, 8(5): 616–623. doi: 10.12000/JR19009

基于辅助数字高程模型的方位多通道SAR相位失配校正方法

DOI: 10.12000/JR19009
基金项目: 国家重点研发计划项目(2017YFB0502700)
详细信息
    作者简介:

    李 强(1979–),男,山西榆次人,博士,助理研究员,2007年在原总装备部装备学院(现航天工程大学)获得博士学位,现为北京跟踪与通信技术研究所助理研究员。主要研究方向为航天工程总体。E-mail: lenovo1979lq@163.com

    范怀涛(1990–),男,山东潍坊人,博士,助理研究员,2017年在中国科学院大学获得博士学位,现为中国科学院电子学研究所助理研究员。主要研究方向为方位多通道合成孔径雷达信号处理。E-mail: huaitaofan@163.com

    通讯作者:

    范怀涛  huaitaofan@163.com

  • 责任主编:李宁 Corresponding Editor: LI Ning
  • 中图分类号: TN957

Channel Phase Mismatch Calibration for Multichannel in Azimuth SAR Imaging Aided by Digital Elevation Model

Funds: The National Key R&D Program of China (2017YFB0502700)
More Information
  • 摘要: 在方位多通道合成孔径雷达(SAR)系统中,进行非均匀采样重建之前,由于通道特性不一致导致的幅度相位差异必须进行校正,以避免图像中出现“鬼影”虚假目标,影响图像判读。方位多通道SAR工作过程中,平台偏航和俯仰导致的通道相位失配具有方位时变和距离空变特点。目前基于平台姿态信息的通道相位失配校正方法均未考虑地形高程起伏带来的影响。该文提出一种新的方位多通道SAR相位失配校正方法,基于辅助数字高程模型(DEM)信息和平台姿态信息,获得更加精确的场景下视角,在地形起伏较大的场景显著提高了通道间相位失配估计精度。针对提出的算法,开展仿真实验,针对虚假目标抑制效果开展定量评估。同时选取场景高程起伏较大场景开展了机载飞行试验数据处理,并对实验结果进行分析,验证算法的有效性。

     

  • 合成孔径雷达(Synthetic Aperture Radar, SAR)具备全天候、全天时、远距离、高分辨率对地成像能力,已经成为现代微波遥感领域的重要技术手段之一[1]。获取陆海环境更高的分辨能力和更大的测绘范围一直是各种遥感手段不断突破的方向。方位向多通道体制以空间维度采样的增加消除分辨率和幅宽对时间维度采样率的矛盾制约[2],目前已经成为高分辨率宽幅SAR成像的主流技术手段之一,并成功应用于国内外多颗在轨SAR卫星,包括高分三号[3]、TerraSAR-X[4], RADARSAT-2[5]和ALOS-2[6]

    通道失配校正和多通道信号重建是方位多通道SAR信号处理的两个关键技术环节。多通道信号重建旨在对空间上非均匀分布的多个通道采样信号进行联合处理,得到等效于单通道SAR系统的空间维均匀采样信号[7]。在进行非均匀采样数据重建处理之前,多个通道接收回波的幅度、相位、延时等特性需要校正一致,以避免重建图像中出现干扰图像判读的“鬼影”虚假目标[8]。针对方位多通道SAR系统失配校正问题,德国宇航中心、西安电子科技大学、上海交通大学、中国科学院电子学研究所等国内外研究机构相继进行了大量研究,提出了多种处理方案[911]

    通道相位失配校正主要分为基于内定标系统的标定方法和基于回波数据的估计方法[12]。由于实际SAR系统中通道失配来源复杂,完全依靠内定标系统不仅代价高昂而且系统十分庞大[13]。基于回波数据的通道失配估计方法受到国内外学者的广泛研究[14]。通常来讲,通道间特性失配源于系统特性不理想,与成像场景无关,目前已经发展的通道失配校正方法也是基于此假设。但是,虽然失配特性不一致源于SAR系统特性,其造成的影响与场景特性却是紧密相关的,特别是起伏地形高程。目前发展的通道失配校正方法均未考虑地形起伏的影响,基于平地假设的校正方法在地形起伏明显的场景中,失配校正性能将受到很大影响。当前,高精度数字高程模型(Digital Elevation Model, DEM)重建技术经过国内外学者多年研究,已经逐渐成熟并得到广泛应用[15]。本文详细分析了方位向多通道SAR系统存在偏航、俯仰和地形起伏条件下的通道失配特点,给出运用辅助DEM信息进行补偿的解析表达式,针对仿真实验数据和机载4通道SAR实验数据开展处理工作,验证所提方法的有效性。

    方位多通道SAR信号无模糊重建的基础是各天线子孔径沿航迹方向均匀排布。但是,在实际工作过程中,受高空气流和侧风的影响,天线基线常常偏离平台的飞行轨迹,平台姿态往往存在偏航、俯仰和横滚[16],如图1所示。在星载SAR系统中,地球自转等效于系统中存在一个固定偏航角。

    图  1  机载SAR系统飞行姿态示意图
    Figure  1.  The attitude of airborne SAR system

    其中,平台横滚角带来的影响主要是回波的接收增益发生变化,场景目标到多个通道天线相位中心的斜距变化是相同的,因此不会导致通道间的相位误差。天线偏航和俯仰不仅会导致波束在方位向指向发生改变,成像区域改变,而且会导致子孔径天线相位中心位置发生偏移[16]。天线相位中心偏移导致的相位差异在不同通道中是不同的,由此导致通道间的相位误差。

    图2中,Y 轴指向平台前进方向,Z轴垂直地面向上,X轴形成右手直角坐标系,平台飞行高度H, Ln表示天线子孔径n到参考通道的顺轨间隔,场景目标T坐标为(x, y, h)。以方位向两通道SAR系统为例进行分析,方位向前端通道称为主通道,后端通道称为辅通道,偏移前相位中心表示为P,偏移后相位中心表示为¯P。这里以主通道为参考通道,分析平台偏航和俯仰给辅通道带来的影响,结论同样适用于任意通道数目的多通道SAR系统。窄波束条件下SAR系统符合波束中心近似假设,图3给出平台姿态变化在斜距方向导致相位中心变化示意图。

    图  2  天线相位中心偏移示意图
    Figure  2.  Movement of SAR antenna phase center
    图  3  2维横截面图
    Figure  3.  2-D cross-sectional view

    假设θyaw表示偏航角,θpitch表示俯仰角,如图3所示,则偏航和俯仰导致的相位中心在3维空间的变化量可以表示为

    Δx=LnsinθyawΔz=Lnsinθpitch}
    (1)

    目标T和校正前后辅通道的天线相位中心位置

    T:(x,h)P:(0,H)¯P:(Δx,H+Δz)}
    (2)

    其中,xh分别表示目标T的地距向和高度向的位置坐标,H表示平台飞行高度。由成像几何关系,天线相位中心偏移导致的相位失配可以表示为

    φ=2πλ((xΔx)2+(H+Δzh)2x2+(Hh)2)
    (3)

    其中,λ表示载波波长。在星载SAR和机载SAR系统中,条件ΔxxΔzHh均有效成立,可以得到近似不等式如式(4)

    x2+(Hh)22xΔx+2(Hh)ΔzΔx2+Δz2
    (4)

    在上述假设条件下,将式(3)中的双曲形式泰勒级数展开至2次项可以得到

    φ2πλ(xΔx+(Hh)Δzx2+(Hh)2)
    (5)

    基于成像几何关系,可以得到下视角θlook的三角函数关系如式(6)

    cosθlook=(Hh)x2+(Hh)2sinθlook=xx2+(Hh)2}
    (6)

    综上,将式(6)代入式(5)可以得到偏航角和俯仰角导致的通道间相位偏差

    φ=2πLnλ(sinθlooksinθyaw+cosθlooksinθpitch)
    (7)

    在机载SAR系统工作过程中,偏航和俯仰不断发生变化,通道间相位失配具有方位时变特性。同时,从式(7)可以看出,相位失配随场景地物对应的下视角变化,呈现距离向空变特点,因此式(7)可以重写为

    φ(t,h)=2πLnλ(sin(θlook(h))sin(θyaw(t))+cos(θlook(h))sin(θpitch(t)))
    (8)

    在地形起伏明显的区域,场景地物对应的下视角与平地假设下不同。如果采用简单的平地假设,针对存在高程起伏的地物,通道间相位失配补偿将出现错误,导致最终的重建图像中出现虚假目标。如图4所示,假设平台飞行高度为3000 m,场景中高程地物高出参考水平面300 m,建筑物顶端视角θ1为50°。如果不考虑场景地形起伏,建筑物顶端斜距对应的视角θ2约为44.42°。在子孔径间隔0.5 m的C波段机载SAR中,假设机载系统存在5°偏航角和3°俯仰角,采用平地假设进行相位失配校正,造成的偏差约为6.5°,在背景环境较弱的条件下会出现干扰虚假目标。如果在校正相位失配时考虑地形起伏影响,可以准确校正姿态导致的通道间相位失配,有效抑制重建后图像中的虚假目标。

    图  4  高程引起的成像叠掩示意图
    Figure  4.  Imaging overlap caused by topograghic relief

    基于辅助DEM的方位多通道SAR通道失配补偿流程如图5所示,具体处理步骤如下:

    图  5  基于DEM辅助的方位多通道SAR通道失配补偿方法流程图
    Figure  5.  The flowchart of the proposed channel phase mismatch calibration method based on the aided DEM

    (1) 根据惯性测量单元(IMU)记录的3维速度信息,针对多通道回波数据进行1阶运动补偿,补偿之后认为参考通道处于理想直线航迹上;

    (2) 在距离频域和方位时域进行2维自适应校正幅度误差,校正表达式如式(9)所示,其中n表示对通道2进行校正的次数,fr表示距离频率,t表示方位时间;

    Sn+12(fr,t)=Sn2(fr,t)|S1(fr,t)|dt|Sn2(fr,t)|dtSn+22(fr,t)=Sn+12(fr,t)|S1(fr,t)|dfr|Sn+12(fr,t)|dfr}
    (9)

    (3) 进行距离压缩;

    (4) 采用方位互相关方法[11]进行通道间延时误差校正;

    (5) 在机载系统中,飞行姿态具有方位向时变特点,场景高程起伏也具有方位空变特点,采用方位向子孔径处理方式,假设子孔径内方位向高程不变,利用辅助DEM信息计算场景地物对应的下视角,并结合IMU记录的平台姿态数据校正距离空变、方位时变相位误差;

    (6) 采用子空间投影方法[13]校正残余的非空变、非时变相位误差;

    (7) 采用滤波器组方法[7]进行多通道信号重建,得到无模糊的单通道信号;

    (8) 最后采用调频变标(CS)算法进行单通道数据成像处理,得到无模糊SAR图像。

    为了验证本文算法有效性,分别开展仿真数据处理和机载多通道SAR试验数据处理。首先,参照实际飞行试验中采用的方位多通道机载SAR系统设置仿真实验,实验参数由表1给出。

    表  1  仿真机载实验参数
    Table  1.  System parameters of the simulated airborne experiment
    系统参数数值
    通道数目4
    飞行高度(m)3000
    平台速度(m/s)120
    载频(GHz)5.4
    发射信号带宽(MHz)210
    脉冲重复频率(Hz)150
    多普勒带宽(Hz)384
    下载: 导出CSV 
    | 显示表格

    假设平台存在5°偏航角和3°俯仰角。在场景中设置5个点目标,等间距分布,位于中间位置的点目标高程为0,其他4个点目标分别设置不同的高度,如图6所示。

    图  6  场景点目标间距及高程设置
    Figure  6.  Five points set in the simulated scene

    图7展示了仿真实验数据处理结果。图7(a)给出的是未进行通道失配校正的多通道信号重建结果。由于通道失配的影响,场景中5个点目标均存在明显的虚假目标,分别分布在真实目标方位向前后位置。由于虚假目标的位置与目标的最近斜距有关,对比A, B和C点不难发现位于距离向不同位置点目标的虚假目标与真实目标的间距不同。同时,对比B点和E点可以发现,由于高程的影响,位于相同地距位置的点目标在斜距图像中出现在不同距离门中。

    图  7  仿真实验结果
    Figure  7.  Experimental results of the simulation

    图7(b)中给出平地假设条件下通道失配校正结果。可以发现,没有高程的B点虚假目标得到有效抑制。由于未考虑地形起伏影响,场景中A, C, D, E 4个点目标仍然存在可见的虚假目标。同时,高程较高的C点和E点,虚假目标强度明显高于高程较低的A点和D点。图7(c)中给出考虑地形起伏条件下通道失配校正结果。可以发现,场景中A, B, C, D, E 5个点目标对应的虚假目标均得到有效抑制。仿真实验中虚假目标抑制的定量评估结果在表2中给出。

    表  2  仿真实验定量评估结果(dB)
    Table  2.  Quantitative assessment result of the simulation (dB)
    目标未校正未考虑高程起伏考虑高程起伏
    A–5–34–53
    B–7–50–50
    C–5–28–53
    D–6–37–51
    E–5–27–52
    下载: 导出CSV 
    | 显示表格

    中国科学院电子学研究所自主研制C波段方位4通道机载SAR系统,于2017年7月在舟山开展飞行试验,主要系统参数在表3中给出。系统采用1发4收的工作方式,整个天线阵面发射信号,方位向均匀分成4个通道接收回波。图8给出飞行过程中偏航和俯仰角度。

    表  3  C波段方位向4通道机载SAR实验参数
    Table  3.  Experimental parameters of the C-band azimuth four channel airborne SAR system
    系统参数数值
    通道数目4
    平台高度(m)3000
    平台速度(m/s)129.64
    载频(GHz)5.4
    发射信号带宽(MHz)210
    方位向天线尺寸(m)0.624
    脉冲重复频率(Hz)137.19
    下载: 导出CSV 
    | 显示表格
    图  8  试验过程中的平台姿态信息
    Figure  8.  The attitude information recorded by the IMU

    图9展示了方位4通道机载SAR试验数据处理结果。图9(a)给出的是未进行通道失配校正的多通道信号重建结果。由于通道失配的影响,场景中存在明显的虚假目标,干扰图像的正常判读。特别是左岸煤炭码头停靠的干散货船,散射强度较大,对应的虚假目标非常明显。所选区域在方位向临海,虚假目标在平静的海面上尤其明显。图9(b)中给出平地假设条件下通道失配校正结果。可以发现,虚假目标的影响大大减弱。但是由于存在高程的山地影响,图像中强散射的目标仍然存在可见的虚假目标。图9(c)中给出考虑地形起伏条件下通道失配校正结果。可以发现,场景中虚假目标均得到有效抑制。

    图  9  机载SAR数据处理结果
    Figure  9.  Imaging results of the airborne SAR experimental data

    方位向多通道星载SAR是实现高分辨率宽幅成像的重要技术手段之一,机载方位多通道SAR是验证系统和算法有效性的必要手段。基于外部DEM辅助信息,本文提出一种适用于高程起伏条件下,方位多通道条带SAR系统的通道相位失配校正方法。仿真实验表明算法能够精确补偿平台姿态导致的通道相位失配。将本文算法应用于机载多通道SAR实验数据处理,在高程起伏区域获取了良好的通道一致性校正和重建处理结果,验证了本文算法的有效性。

  • 图  1  机载SAR系统飞行姿态示意图

    Figure  1.  The attitude of airborne SAR system

    图  2  天线相位中心偏移示意图

    Figure  2.  Movement of SAR antenna phase center

    图  3  2维横截面图

    Figure  3.  2-D cross-sectional view

    图  4  高程引起的成像叠掩示意图

    Figure  4.  Imaging overlap caused by topograghic relief

    图  5  基于DEM辅助的方位多通道SAR通道失配补偿方法流程图

    Figure  5.  The flowchart of the proposed channel phase mismatch calibration method based on the aided DEM

    图  6  场景点目标间距及高程设置

    Figure  6.  Five points set in the simulated scene

    图  7  仿真实验结果

    Figure  7.  Experimental results of the simulation

    图  8  试验过程中的平台姿态信息

    Figure  8.  The attitude information recorded by the IMU

    图  9  机载SAR数据处理结果

    Figure  9.  Imaging results of the airborne SAR experimental data

    表  1  仿真机载实验参数

    Table  1.   System parameters of the simulated airborne experiment

    系统参数数值
    通道数目4
    飞行高度(m)3000
    平台速度(m/s)120
    载频(GHz)5.4
    发射信号带宽(MHz)210
    脉冲重复频率(Hz)150
    多普勒带宽(Hz)384
    下载: 导出CSV

    表  2  仿真实验定量评估结果(dB)

    Table  2.   Quantitative assessment result of the simulation (dB)

    目标未校正未考虑高程起伏考虑高程起伏
    A–5–34–53
    B–7–50–50
    C–5–28–53
    D–6–37–51
    E–5–27–52
    下载: 导出CSV

    表  3  C波段方位向4通道机载SAR实验参数

    Table  3.   Experimental parameters of the C-band azimuth four channel airborne SAR system

    系统参数数值
    通道数目4
    平台高度(m)3000
    平台速度(m/s)129.64
    载频(GHz)5.4
    发射信号带宽(MHz)210
    方位向天线尺寸(m)0.624
    脉冲重复频率(Hz)137.19
    下载: 导出CSV
  • [1] GEBERT N, KRIEGER G, and MOREIRA A. High resolution wide swath SAR imaging with digital beam-forming-performance analysis, optimization, system design[C]. Proceedings of the 6th European Conference on Synthetic Aperture Radar, Dresden, Germany, 2006: 341–344.
    [2] ZHANG Lei, XING Mengdao, QIU Chengwei, et al. Adaptive two-step calibration for high resolution and wide-swath SAR imaging[J]. IET Radar, Sonar & Navigation, 2010, 4(4): 548–559. doi: 10.1049/iet-rsn.2008.0158
    [3] SUN Jili, YU Weidong, and DENG Yunkai. The SAR payload design and performance for the GF-3 mission[J]. Sensors, 2017, 17(10): E2419. doi: 10.3390/s17102419
    [4] KIM J H, YOUNIS M, PRATS-IRAOLA P, et al. First spaceborne demonstration of digital beamforming for azimuth ambiguity suppression[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(1): 579–590. doi: 10.1109/TGRS.2012.2201947
    [5] WILLIAMS D, LEDANTEC P, CHABOT M, et al. RADARSAT-2 image quality and calibration update[C]. Proceedings of 10th European Conference on Synthetic Aperture Radar, Berlin, Germany, 2014: 1–4.
    [6] SHIMADA M. ALOS-2 Science program[C]. 2013 IEEE International Geoscience and Remote Sensing Symposium, Melbourne, Australia, 2013: 2400–2403. doi: 10.1109/IGARSS.2013.6723303.
    [7] KRIEGER G, GEBERT N, and MOREIRA A. Unambiguous SAR signal reconstruction from nonuniform displaced phase center sampling[J]. IEEE Geoscience and Remote Sensing Letters, 2004, 1(4): 260–264. doi: 10.1109/LGRS.2004.832700
    [8] ZHANG Shuangxi, XING Mengdao, XIA Xianggen, et al. Multichannel HRWS SAR imaging based on range-variant channel calibration and multi-Doppler-direction restriction ambiguity suppression[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(7): 4306–4327. doi: 10.1109/TGRS.2013.2281329
    [9] LI Jianbing, LI Xiaoping, LIN Mingfu, et al. Maximum-likelihood-based Doppler centroid estimation algorithm for MC-HRWS SAR system[J]. Electronics Letters, 2014, 50(22): 1630–1631. doi: 10.1049/el.2014.2723
    [10] GUO Xiaojiang, GAO Yesheng, WANG Kaizhi, et al. Improved channel error calibration algorithm for azimuth multichannel SAR systems[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(7): 1022–1026. doi: 10.1109/LGRS.2016.2561961
    [11] FENG Jin, GAO Canguan, ZHANG Yi, et al. Phase mismatch calibration of the multichannel SAR based on azimuth cross correlation[J]. IEEE Geoscience and Remote Sensing Letters, 2013, 10(4): 903–907. doi: 10.1109/LGRS.2012.2227107
    [12] 刘艳阳, 李真芳, 杨桃丽, 等. 一种单星方位多通道高分辨率宽测绘带SAR系统通道相位偏差时域估计新方法[J]. 电子与信息学报, 2012, 34(12): 2913–2919. doi: 10.3724/SP.J.1146.2012.00562

    LIU Yanyang, LI Zhenfang, YANG Taoli, et al. A novel channel phase bias estimation method for spaceborne along-track multi-channel HRWS SAR in time-domain[J]. Journal of Electronics &Information Technology, 2012, 34(12): 2913–2919. doi: 10.3724/SP.J.1146.2012.00562
    [13] 张磊, 全英汇, 邢孟道, 等. 一种子空间投影的高分辨宽测绘带SAR成像通道均衡方法[J]. 电子与信息学报, 2010, 32(1): 1–6. doi: 10.3724/SP.J.1146.2008.01821

    ZHANG Lei, QUAN Yinghui, XING Mengdao, et al. An SSP based channel calibration for high-resolution and wide-swath SAR imagery[J]. Journal of Electronics &Information Technology, 2010, 32(1): 1–6. doi: 10.3724/SP.J.1146.2008.01821
    [14] LIU Yanyang, LI Zhenfang, WANG Zhibin, et al. On the baseband Doppler centroid estimation for multichannel HRWS SAR imaging[J]. IEEE Geoscience and Remote Sensing Letters, 2014, 11(12): 2050–2054. doi: 10.1109/LGRS.2014.2318511
    [15] ZINK M, BACHMANN M, BRAUTIGAM B, et al. TanDEM-X: The new global DEM takes shape[J]. IEEE Geoscience and Remote Sensing Magazine, 2014, 2(2): 8–23. doi: 10.1109/MGRS.2014.2318895
    [16] 范怀涛, 张志敏, 李宁. 基于特征分解的方位向多通道SAR相位失配校正方法[J]. 雷达学报, 2018, 7(3): 346–354. doi: 10.12000/JR17012

    FAN Huaitao, ZHANG Zhimin, and LI Ning. Channel phase mismatch calibration for multichannel in azimuth SAR imaging based on Eigen-structure method[J]. Journal of Radars, 2018, 7(3): 346–354. doi: 10.12000/JR17012
  • 期刊类型引用(1)

    1. 姜文,梁伟,吴一戎. 机载顺轨干涉SAR时变交轨基线校正算法. 国外电子测量技术. 2020(06): 49-54 . 百度学术

    其他类型引用(1)

  • 加载中
图(9) / 表(3)
计量
  • 文章访问数: 2686
  • HTML全文浏览量: 866
  • PDF下载量: 143
  • 被引次数: 2
出版历程
  • 收稿日期:  2019-01-18
  • 修回日期:  2019-03-18
  • 网络出版日期:  2019-10-01

目录

/

返回文章
返回