留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

星载分布式InSAR基线定标新方法

钱方明 姜挺 楼良盛 陈刚 刘薇

钱方明, 姜挺, 楼良盛, 陈刚, 刘薇. 星载分布式InSAR基线定标新方法[J]. 武汉大学学报 ( 信息科学版), 2020, 45(1): 126-133. doi: 10.13203/j.whugis20180199
引用本文: 钱方明, 姜挺, 楼良盛, 陈刚, 刘薇. 星载分布式InSAR基线定标新方法[J]. 武汉大学学报 ( 信息科学版), 2020, 45(1): 126-133. doi: 10.13203/j.whugis20180199
QIAN Fangming, JIANG Ting, LOU Liangsheng, CHEN Gang, LIU Wei. A New Method of Space-Borne Distributed InSAR Baseline Calibration[J]. Geomatics and Information Science of Wuhan University, 2020, 45(1): 126-133. doi: 10.13203/j.whugis20180199
Citation: QIAN Fangming, JIANG Ting, LOU Liangsheng, CHEN Gang, LIU Wei. A New Method of Space-Borne Distributed InSAR Baseline Calibration[J]. Geomatics and Information Science of Wuhan University, 2020, 45(1): 126-133. doi: 10.13203/j.whugis20180199

星载分布式InSAR基线定标新方法

doi: 10.13203/j.whugis20180199
基金项目: 

国家自然科学基金 41371439

详细信息
    作者简介:

    钱方明, 硕士, 助理研究员, 主要从事InSAR数据处理与几何定标方法研究。122748721@qq.com

  • 中图分类号: P237

A New Method of Space-Borne Distributed InSAR Baseline Calibration

Funds: 

The National Natural Science Foundation of China 41371439

More Information
    Author Bio:

    QIAN Fangming, master, assistant researcher, specializes in InSAR data processing and geometry calibration method. E-mail: 122748721@qq.com

  • 摘要: 为了解决星载分布式干涉合成孔径雷达三维基线的高精度定标问题,提出了一种基线定标新方法。首先,该方法将辅雷达多普勒方程与辅雷达距离改化方程相结合,形成基线定标模型,再基于地面控制点,采用最小二乘迭代解算基线三个轴向误差。与基于辅雷达距离改化方程的基线定标模型相比,该基线定标模型优化了法方程系数矩阵条件数,减小了矩阵的病态程度,提高了定标精度。然后,根据基线定标实际上是利用多个地面控制点交会出基线矢量,以及增大交会角可提高交会精度的原理,给出了控制点布设策略,即在近、远距端两侧布设控制点。与均匀布设控制点方式相比,该方式能显著提高定标精度和稳健性。最后,利用星载分布式干涉合成孔径雷达仿真数据进行基线定标实验,先后开展基于不同控制点数量、不同控制点精度、不同控制点布设策略的基线定标实验。实验结果表明,该方法能实现三维基线高精度定标,对跨轨、顺轨、径向基线的定标精度能达到毫米级,可为星载分布式干涉合成孔径雷达基线定标提供一种新途径。
  • 图  1  主雷达APC坐标系中表示的基线矢量示意图

    Figure  1.  Diagram of Baseline Vector in Main Radar APC Coordinate System

    图  2  基线定标原理示意图

    Figure  2.  Diagram of Baseline Calibration Principle

    图  3  控制点子条带位置示意图

    Figure  3.  Diagrams of Control Points Sub-band Location

    图  4  控制点子条带不同位置时定标误差统计

    Figure  4.  Calibration Error Statistics for Sub-bands of Control Points at Different Positions

    图  5  控制点子条带不同位置时标出值与设定值互差

    Figure  5.  Difference Between the Calibrated Value and the Set Value for Sub-bands of Control Points at Different Positions

    表  1  控制点均匀分布的定标实验结果

    Table  1.   Results of Calibration Experiments with Evenly Distributed Control Points

    控制点数量/个 基线分量 系统误差/cm μ/cm σ/cm |δ|/cm εp/cm εh/cm
    20 BX -5.00 -4.48 7.95 0.52 355.31 75.23
    BY -5.00 -5.29 5.60 0.29
    BZ 5.00 5.46 6.99 0.46
    60 BX -5.00 -5.29 4.06 0.29 354.89 75.40
    BY -5.00 -5.10 3.34 0.10
    BZ 5.00 4.75 3.57 0.25
    100 BX -5.00 -4.77 2.81 0.23 355.19 75.23
    BY -5.00 -5.05 2.43 0.05
    BZ 5.00 5.20 2.47 0.20
    140 BX -5.00 -4.85 2.65 0.15 355.33 75.21
    BY -5.00 -4.88 2.13 0.12
    BZ 5.00 5.13 2.33 0.13
    180 BX -5.00 -5.05 2.25 0.05 355.44 75.26
    BY -5.00 -5.13 2.02 0.13
    BZ 5.00 4.95 1.98 0.05
    下载: 导出CSV

    表  2  不同控制点精度的定标实验结果

    Table  2.   Results of Calibration Experiments with Control Points of Different Accuracy

    控制点精度/m 基线分量 系统误差/cm μ/cm σ/cm |δ|/cm εp/cm εh/cm
    2.0 BX -5.00 -5.45 4.59 0.45 357.35 76.22
    BY -5.00 -5.65 22.38 0.65
    BZ 5.00 4.61 4.03 0.39
    1.0 BX -5.00 -5.36 4.14 0.36 354.88 75.42
    BY -5.00 -5.33 11.18 0.33
    BZ 5.00 4.69 3.64 0.31
    0.5 BX -5.00 -5.31 4.06 0.31 354.88 75.40
    BY -5.00 -5.16 5.58 0.16
    BZ 5.00 4.73 3.57 0.27
    0.1 BX -5.00 -5.27 4.08 0.27 354.89 75.39
    BY -5.00 -5.03 1.10 0.03
    BZ 5.00 4.76 3.59 0.24
    下载: 导出CSV

    表  3  控制点子条带不同位置的定标结果

    Table  3.   Calibration Results for Sub-bands of Control Points at Different Positions/cm

    子条带位置 基线分量 系统误差 μ σ |δ| εp εh
    距离向中部 BX -5.00 -3.79 22.26 1.21 356.19 75.37
    BY -5.00 -4.74 3.57 0.26
    BZ 5.00 6.06 19.58 1.06
    距离向宽度1/3、2/3处 BX -5.00 -5.35 6.35 0.35 355.04 75.38
    BY -5.00 -4.87 3.44 0.13
    BZ 5.00 4.69 5.58 0.31
    近、远距端 BX -5.00 -5.17 2.28 0.17 355.12 75.30
    BY -5.00 -4.87 3.44 0.13
    BZ 5.00 4.85 2.00 0.15
    下载: 导出CSV
  • [1] 张永俊.星载分布式InSAR系统的误差分析与DEM精度提高方法研究[D].长沙: 国防科学技术大学, 2011

    Zhang Yongjun. Study on Error Analysis and DEM Precision Improvement Methods of Spaceborne Distributed InSAR[D]. Changsha: National University of Defense Technology, 2011
    [2] 楼良盛, 汤晓涛, 牛瑞.基于卫星编队InSAR空间同步对系统性能影响的分析[J].武汉大学学报·信息科学版, 2007, 32(10):892-894 http://ch.whu.edu.cn/CN/abstract/abstract1999.shtml

    Lou Liangsheng, Tang Xiaotao, Niu Rui. Analysis of Influence on System Performance from InSAR Space Synchronization Based on Formation Satellites Based InSAR[J]. Geomatics and Information Science of Wuhan University, 2007, 32(10):892-894 http://ch.whu.edu.cn/CN/abstract/abstract1999.shtml
    [3] Tu J, Gu D F, Wu Y, et al. Error Modeling and Analysis for InSAR Spatial Baseline Determination of Satellite Formation Flying[J]. Mathematical Problems in Engineering, 2012(4):346-360 http://d.old.wanfangdata.com.cn/OAPaper/oai_doaj-articles_205ad9312f31c6cdc4e5b80e8ef2c6d1
    [4] Small D. Baseline Modelling for ERS-1 SAR Interferometry[C]. International Geoscience and Remote Sensing Symposium, Piscataway, New Jersey, USA, 1993 https://www.researchgate.net/publication/224665930_Baseline_modelling_for_ERS-1_SAR_interferometry
    [5] Knedlik S, Loffeld O. Baseline Estimation and Prediction Referring to the SRTM[C]. International Geoscience and Remote Sensing Symposium, Piscataway, New Jersey, USA, 2002 Baseline Estimation and Prediction Referring to the SRTM
    [6] Zink M, Geudtner D. Calibration of the Interferometric X-SAR System on SRTM[C]. International Geoscience and Remote Sensing Symposium, Piscataway, New Jersey, USA, 1999 https://www.researchgate.net/publication/3803258_Calibration_of_the_interferometric_X-SAR_system_on_SRTM
    [7] Costantini M, Minati F, Quagliarini A, et al. SAR Interferometric Baseline Calibration Without Need of Phase Unwrapping[C]. International Geoscience and Remote Sensing Symposium, Piscataway, New Jersey, USA, 2004 https://www.researchgate.net/publication/4110474_SAR_interferometric_baseline_calibration_without_need_of_phase_unwrapping
    [8] 刘国祥.利用雷达干涉测量技术监测区域地表形变[M].北京:测绘出版社, 2006

    Liu Guoxiang. Monitoring Surface Deformation Area by Radar Interferometry[M]. Beijing: Surveying and Mapping Press, 2006
    [9] 何秀凤, 何敏. InSAR对地观测数据处理方法与综合测量[M].北京:科学出版社, 2012:63-65

    He Xiufeng, He Min. InSAR Earth Observation Data Processing Method and Comprehensive Measurement[M] Beijing: Science Press, 2012:63-65
    [10] Gonzalez J H, Bachmann M, Krieger G, et al. Development of the TanDEM-X Calibration Concept: Analysis of Systematic Errors[J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(2): 716–726 doi:  10.1109/TGRS.2009.2034980
    [11] Bachmann M, Gonzalez J H, Krieger G, et al. Calibration of the Bistatic TanDEM-X Interferometer[C]. European Conference on Synthetic Aperture Radar, Berlin, German, 2012
    [12] Antony J W, Gonzalez J H, Schwerdt M, et al. Results of the TanDEM-X Baseline Calibration[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2013, 6(3): 1 495–1 501 doi:  10.1109/JSTARS.2013.2258328
    [13] 陈刚, 汤晓涛, 钱方明.星载InSAR立体基线定标方法[J].武汉大学学报·信息科学版, 2014, 39(1):37-41 http://ch.whu.edu.cn/CN/abstract/abstract2860.shtml

    Chen Gang, Tang Xiaotao, Qian Fangming. A Three-Dimensional Baseline Calibration Method of Spaceborne InSAR[J]. Geomatics and Information Science of Wuhan University, 2014, 39(1): 37-41 http://ch.whu.edu.cn/CN/abstract/abstract2860.shtml
    [14] Zhao S Y, Gu D F, Yi B, et al. Error Analysis for the Baseline Estimation and Calibration of Distributed InSAR Satellites[C]. International Geoscience and Remote Sensing Symposium, Piscataway, New Jersey, USA, 2016
    [15] 楼良盛, 汤晓涛, 黄启来.基于卫星编队InSAR方位向基线影响分析[J].武汉大学学报·信息科学版, 2007, 32(1):59-61 http://ch.whu.edu.cn/CN/abstract/abstract1799.shtml

    Lou Liangsheng, Tang Xiaotao, Huang Qilai. Influence of Azimuth Base-Line on Formation-Flying Satellites Based InSAR[J]. Geomatics and Information Science of Wuhan University, 2007, 32(1): 59-61 http://ch.whu.edu.cn/CN/abstract/abstract1799.shtml
  • [1] 黄攀, 唐劲松, 钟何平.  干涉合成孔径声呐复图像配准分段曲面拟合法 . 武汉大学学报 ( 信息科学版), 2021, 46(8): 1259-1264. doi: 10.13203/j.whugis20190306
    [2] 杨梦诗, 廖明生, 史绪国, 张路.  联合多平台InSAR数据集精确估计地表沉降速率场 . 武汉大学学报 ( 信息科学版), 2017, 42(6): 797-802. doi: 10.13203/j.whugis20140924
    [3] 李宁, 汪长城, 彭星, 付海强, 李鹏飞, 熊文秀.  一种双基线极化干涉SAR的建筑物高度估计方法 . 武汉大学学报 ( 信息科学版), 2016, 41(10): 1333-1338. doi: 10.13203/j.whugis20140675
    [4] 王利朋, 刘成龙, 刘胜, 何林烜.  坐标参数化的平面基准转换方法 . 武汉大学学报 ( 信息科学版), 2016, 41(10): 1409-1413. doi: 10.13203/j.whugis20140416
    [5] 孙文川, 暴景阳, 金绍华, 肖付民, 张志伟.  一种多波束换能器横摇角度偏差二次校准方法 . 武汉大学学报 ( 信息科学版), 2016, 41(11): 1440-1444. doi: 10.13203/j.whugis20140481
    [6] 李晶, 李冬海, 赵拥军.  利用角度和时差的单站外辐射源定位方法 . 武汉大学学报 ( 信息科学版), 2015, 40(2): 227-232.
    [7] 徐亚明, 周建国, 张鹏.  一种WiFi无线接入点位置估计新算法 . 武汉大学学报 ( 信息科学版), 2014, 39(3): 279-282. doi: 10.13203/j.whugis20120051
    [8] 师芸.  加乘性混合误差模型参数估计方法及其应用 . 武汉大学学报 ( 信息科学版), 2014, 39(9): 1033-1037. doi: 10.13203/j.whugis20130355
    [9] 陈刚, 汤晓涛, 钱方明.  星载InSAR立体基线定标方法 . 武汉大学学报 ( 信息科学版), 2014, 39(1): 37-41.
    [10] 李永生, 张景发, 李振洪, 罗毅.  利用短基线集干涉测量时序分析方法监测北京市地面沉降 . 武汉大学学报 ( 信息科学版), 2013, 38(11): 1374-1377.
    [11] 闫利, 胡晓斌.  利用GA求解卫星影像的空间后方交会 . 武汉大学学报 ( 信息科学版), 2013, 38(11): 1286-1289.
    [12] 柯灏, 赵建虎, 张红梅.  短潮位序列系统误差的探测及修复方法研究 . 武汉大学学报 ( 信息科学版), 2012, 37(7): 843-846.
    [13] 周志伟, 鄢子平, 刘苏, 李振洪.  永久散射体与短基线雷达干涉测量在城市地表形变中的应用 . 武汉大学学报 ( 信息科学版), 2011, 36(8): 928-931.
    [14] 黄登山, 杨敏华, 胥海威, 姚学恒.  利用最优估计理论进行多光谱与全色影像融合 . 武汉大学学报 ( 信息科学版), 2011, 36(9): 1039-1042.
    [15] 肖金群, 汪长城, 李志伟, 杨亚夫.  基于配准偏移量的InSAR基线估计 . 武汉大学学报 ( 信息科学版), 2010, 35(10): 1236-1239.
    [16] 李国平, 陈娇娜, 黄丁发, 郭洁.  地基GPS水汽实时监测系统及其气象业务应用 . 武汉大学学报 ( 信息科学版), 2009, 34(11): 1328-1331.
    [17] 童小华, 杜道生.  地图数字化中基于点位坐标的统一平差模型 . 武汉大学学报 ( 信息科学版), 2006, 31(3): 194-198.
    [18] 罗红霞, 龚健雅.  UFCLS线性光谱混合分析法在遥感图像分类识别中的效果分析 . 武汉大学学报 ( 信息科学版), 2004, 29(7): 615-618.
    [19] 邓德祥, 吴章华, 胡志雄.  一种用于影像镶嵌的快速最小二乘序贯算法 . 武汉大学学报 ( 信息科学版), 1997, 22(1): 29-31.
    [20] 陶本藻.  论最小二乘拟合 . 武汉大学学报 ( 信息科学版), 1980, 5(1): 27-34.
  • 加载中
图(5) / 表(3)
计量
  • 文章访问数:  1174
  • HTML全文浏览量:  224
  • PDF下载量:  160
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-12-20
  • 刊出日期:  2020-01-05

星载分布式InSAR基线定标新方法

doi: 10.13203/j.whugis20180199
    基金项目:

    国家自然科学基金 41371439

    作者简介:

    钱方明, 硕士, 助理研究员, 主要从事InSAR数据处理与几何定标方法研究。122748721@qq.com

  • 中图分类号: P237

摘要: 为了解决星载分布式干涉合成孔径雷达三维基线的高精度定标问题,提出了一种基线定标新方法。首先,该方法将辅雷达多普勒方程与辅雷达距离改化方程相结合,形成基线定标模型,再基于地面控制点,采用最小二乘迭代解算基线三个轴向误差。与基于辅雷达距离改化方程的基线定标模型相比,该基线定标模型优化了法方程系数矩阵条件数,减小了矩阵的病态程度,提高了定标精度。然后,根据基线定标实际上是利用多个地面控制点交会出基线矢量,以及增大交会角可提高交会精度的原理,给出了控制点布设策略,即在近、远距端两侧布设控制点。与均匀布设控制点方式相比,该方式能显著提高定标精度和稳健性。最后,利用星载分布式干涉合成孔径雷达仿真数据进行基线定标实验,先后开展基于不同控制点数量、不同控制点精度、不同控制点布设策略的基线定标实验。实验结果表明,该方法能实现三维基线高精度定标,对跨轨、顺轨、径向基线的定标精度能达到毫米级,可为星载分布式干涉合成孔径雷达基线定标提供一种新途径。

English Abstract

钱方明, 姜挺, 楼良盛, 陈刚, 刘薇. 星载分布式InSAR基线定标新方法[J]. 武汉大学学报 ( 信息科学版), 2020, 45(1): 126-133. doi: 10.13203/j.whugis20180199
引用本文: 钱方明, 姜挺, 楼良盛, 陈刚, 刘薇. 星载分布式InSAR基线定标新方法[J]. 武汉大学学报 ( 信息科学版), 2020, 45(1): 126-133. doi: 10.13203/j.whugis20180199
QIAN Fangming, JIANG Ting, LOU Liangsheng, CHEN Gang, LIU Wei. A New Method of Space-Borne Distributed InSAR Baseline Calibration[J]. Geomatics and Information Science of Wuhan University, 2020, 45(1): 126-133. doi: 10.13203/j.whugis20180199
Citation: QIAN Fangming, JIANG Ting, LOU Liangsheng, CHEN Gang, LIU Wei. A New Method of Space-Borne Distributed InSAR Baseline Calibration[J]. Geomatics and Information Science of Wuhan University, 2020, 45(1): 126-133. doi: 10.13203/j.whugis20180199
  • 星载分布式干涉合成孔径雷达(interferometric synthetic aperture radar,InSAR)系统是指利用多颗满足一定编队构形的合成孔径雷达(synthetic aperture radar,SAR)卫星对地面同一地区进行一发多收、自发自收、多发多收或交替发射各自接收式干涉成像,获取SAR复影像对的InSAR测量系统[1-2]。星载分布式InSAR已成为星载InSAR的发展趋势,与重复轨道InSAR相比,其回波之间具有较好的相干性,数据可用率高、产品精度高;与单平台双天线InSAR相比,其工程实现难度较低,通过编队构形切换可支撑多任务多模式。

    InSAR测量中的干涉基线是指两颗雷达对同一地面目标观测时,主辅雷达天线相位中心(antenna phase centre,APC)的连线[3],简称基线。基线是干涉测量的基础,InSAR测量系统的模糊高度和相位的高程敏感度都与基线长度直接相关。此外,基线矢量误差是影响高程精度的主要因素之一,若星载分布式InSAR系统的模糊高度为45 m,雷达波长为3.2 cm(X波段),在主雷达视线方向上1 cm的基线误差会带来14 m左右的高程误差。基线定标的任务是利用地面控制信息(控制点、数字高程模型数据等)确定基线的系统误差,从而确定真实基线矢量。分布式InSAR系统的设计指标中一般要求基线定标的精度在毫米级。基线定标受到了国内外许多学者的关注,文献[4-9]通过建立切向、交向、径向坐标系,对重复轨道、单平台双天线或分布式InSAR基线进行了定标,此类方法中大多数没有直接给出定标方法的精度,而是利用标定后基线计算地面点位置,通过地面点的高程和平面定位精度来评估基线定标算法,且只实现了二维基线(跨轨、径向基线或基线长度、基线方位角)定标,其中何秀凤等[9]提出的基于卡尔曼滤波和配准参数的基线估计方法,以精轨参数计算的基线为真值,该方法得到的二维基线估计精度为分米级。文献[10-12]在TanDEM系统的定标中将基线误差分解为平行基线误差和垂直基线误差,采用不同波位获取的相邻两景数据联合解算二维基线误差,定标后的基线误差优于2 mm,该方法对数据获取要求较高,且需要较大面积的定标场支持,其定标场选在了澳大利亚北部,范围400 km×800 km。通常情况下星载InSAR的基线是三维的,分布式InSAR采用绕飞构型,两星沿航向存在较长的距离(德国TanDEM-X系统沿航向距离最大可达到500 m[10]),因此有必要对顺轨基线进行定标。陈刚等[13]开展了星载InSAR基线定标方法研究,提出了一种星载InSAR立体基线定标方法,并利用仿真数据进行了实验,顺轨基线的精度在分米级,其他两个方向基线的精度在厘米级,该方法定标精度较低,无法满足星载分布式InSAR系统毫米级定标精度要求。

    针对以上方法的不足,本文提出了一种新的星载分布式InSAR三维基线定标方法。首先从改善法方程系数矩阵状态的角度出发,将辅雷达距离改化方程与辅雷达多普勒方程组合,形成新的基线定标模型,给出了基线误差解算方法;然后通过对基线定标原理进行分析,提出了一种在近、远距端两侧布设控制点的策略;最后通过仿真数据实验验证了新基线定标模型和控制点布设策略的正确性和合理性。

    • 基线矢量可以表示为辅雷达APC相对于主雷达APC的偏移,因此以主雷达APC为原点建立空间直角坐标系,就可以用3个轴向的坐标表示基线矢量。参考卫星轨道数据处理中的卫星中心坐标系[2]定义建立主雷达APC坐标系,如图 1所示。A1为主雷达APC,与坐标原点O重合;A2为辅雷达APC,$P\left({{X_d}, {Y_d}, {Z_d}} \right)$为地面点;R1R2分别为主、辅雷达斜距;B为基线矢量;S为主雷达APC在地固系的位置矢量;V为主雷达速度矢量。主雷达速度矢量方向定义为Y轴,称为顺轨方向;Y轴与S矢量确定平面的法向量方向定义为X轴,称为跨轨方向;X轴和Y轴的正交方向为Z轴,称为径向。这样对于每一方位时刻,基线矢量在主雷达APC坐标系中可分解为顺轨基线BY、跨轨基线BX和径向基线BZ

      图  1  主雷达APC坐标系中表示的基线矢量示意图

      Figure 1.  Diagram of Baseline Vector in Main Radar APC Coordinate System

    • 分布式InSAR对地定位模型可以用主、辅雷达距离和多普勒方程表示,在这4个方程中辅雷达距离方程和辅雷达多普勒方程包含基线参数,根据三角形余弦定理可对辅雷达距离方程进行改化[13],将辅雷达距离改化方程和多普勒方程组成基线定标模型,在主雷达APC坐标系下,该模型可表示为:

      $$\left\{ {\begin{array}{*{20}{l}} {{F_R} = R_1^2 + B_X^2 + B_Y^2 + B_Z^2 - 2{B_X}{X_d} - }\\ {\begin{array}{*{20}{c}} {}&{2{B_Y}{Y_d} - 2{B_Z}{Z_d} - {{\left( {{R_1} - \frac{{\lambda {\phi _i}}}{{2\rho {\rm{ \mathit{ π} }}}}} \right)}^2} = 0} \end{array}}\\ {{F_D} = {V_{{X_2}}}\left( {{B_X} - {X_d}} \right) + {V_{{Y_2}}}\left( {{B_Y} - {Y_d}} \right) + }\\ {\begin{array}{*{20}{c}} {}&{{V_{{Z_2}}}\left( {{B_Z} - {Z_d}} \right) + \frac{{\lambda {R_2}{f_{{d_2}}}}}{2} = 0} \end{array}} \end{array}} \right.$$ (1)

      式中,$\left({{B_X}, {B_Y}, {B_Z}} \right)$为基线矢量;λ为雷达波长;ϕi为地面点的绝对干涉相位;ρ表示收发模式,ρ=1表示一发双收模式,ρ=2表示自发自收模式;$\left({{V_{{X_2}}}, {V_{{Y_2}}}, {V_{{Z_2}}}} \right)$为辅雷达速度矢量;${f_{{d_2}}}$为辅雷达多普勒中心频率。式(1)为非线性方程,若想实现各参数的平差,需对各参数求导,将其线性化。线性化后,式(1)可写成:

      $$\left\{ {\begin{array}{*{20}{l}} {{F_R} = {F_{{R_0}}} + {a_0}{\rm{\Delta }}{B_X} + {a_1}{\rm{\Delta }}{B_Y} + {a_2}{\rm{\Delta }}{B_Z} = 0}\\ {{F_D} = {F_{{D_0}}} + a_0^{'}{\rm{\Delta }}{B_X} + a_1^{'}{\rm{\Delta }}{B_Y} + a_2^{'}{\rm{\Delta }}{B_Z} = 0} \end{array}{\rm{}}} \right.$$ (2)

      式中,未知数系数分别为:

      $$\left\{ {\begin{array}{*{20}{l}} {{a_0} = \frac{{\partial {F_R}}}{{\partial {B_X}}} = 2\left( {{B_X} - {X_d}} \right)}\\ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{l}} {{a_1} = \frac{{\partial {F_R}}}{{\partial {B_Y}}} = 2\left( {{B_Y} - {Y_d}} \right)}\\ {{a_2} = \frac{{\partial {F_R}}}{{\partial {B_Z}}} = 2\left( {{B_Z} - {Z_d}} \right)} \end{array}}&{} \end{array}}\\ {a_0^{'} = \frac{{\partial {F_D}}}{{\partial {B_X}}} = {V_{{X_2}}}}\\ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{l}} {a_1^{'} = \frac{{\partial {F_D}}}{{\partial {B_Y}}} = {V_{{Y_2}}}}\\ {a_2^{'} = \frac{{\partial {F_D}}}{{\partial {B_Z}}} = {V_{{Z_2}}}} \end{array}}&{}&{} \end{array}} \end{array}} \right.$$ (3)

      常数项为:

      $${F_0} = \left[ {\begin{array}{*{20}{l}} {{F_{{R_0}}}}\\ {{F_{{D_0}}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {R_1^2 + B_{{x_0}}^2 + B_{{Y_0}}^2 + B_{{z_0}}^2 - 2{B_{{x_0}}}{X_d} - }\\ {2{B_{{Y_0}}}{Y_d} - 2{B_{{z_0}}}{Z_d} - {{\left( {{R_1} - \frac{{\lambda {\phi _i}}}{{2\rho {\rm{ \mathit{ π} }}}}} \right)}^2}}\\ {{V_{{X_2}}}\left( {{B_{{X_0}}} - {X_d}} \right) + {V_{{Y_2}}}\left( {{B_{{Y_0}}} - {Y_d}} \right) + }\\ {{V_{{Z_2}}}\left( {{B_{{Z_0}}} - {Z_d}} \right) + \frac{{\lambda {R_2}{f_{{d_2}}}}}{2}} \end{array}} \right]$$ (4)

      式中,$\left({{B_{{X_0}}}, {B_{{Y_0}}}, {B_{{Z_0}}}} \right)$为基线矢量初值。

      当有nn≥2)个地面控制点,则可以采用最小二乘法求解,按照式(2)列出误差方程组。实际计算时,在式(2)第一个方程中,Yd的绝对值相对于XdZd较小(若主雷达在零多普勒面内成像的条件下,Yd值为0),系数a1的绝对值相对于a0a2较小;在式(2)第二个方程中,系数a'1为速度在Y轴方向也就是卫星飞行方向的分量,一般为7.6 km/s左右,它相对其他两个轴方向的分量要大得多,这与在第一个方程中的情况正好相反。因此,由两个方程组成的定标模型与只有一个方程的定标模型[13]相比,法方程系数矩阵条件数更优,矩阵的病态程度更小,最小二乘解算精度更高。

    • 布设控制点的常规方法是在场景内均匀布设,为了提高定标模型平差解算的精度和稳定性,一般应尽量多布设控制点,但在实际情况下,若场景内地形复杂,在某些区域布放控制点难度很大,另外由于成本和时间等因素,可能无法进行大量控制点的布设,因此如何能以较少的控制点实现较高的定标精度是值得研究的问题。除了改进模型外,优化控制点的布设策略也是一种解决问题的途径。图 2为基线定标原理示意图,通过干涉方式可获取高精度斜距R2R3,在已知地面点P1P2点坐标后,以斜距R2R3为半径做圆弧,两个圆弧交于A2,从而求解出基线矢量。在其他误差条件不变情况下,要想提高交会精度,就必须拉大P1P2点之间距离,即增加$\angle {P_1}{A_2}{P_2}$的角度,因此,当采用多个控制点进行基线定标时,控制点应尽量布设在近距端和远距端附近。

      图  2  基线定标原理示意图

      Figure 2.  Diagram of Baseline Calibration Principle

    • 为了验证算法的正确性,本文利用InSAR仿真数据进行定标实验,仿真数据的优势是可以设定与实际情况相似的误差值,通过比较定标出的误差值与设定值之间的偏差,来确定定标算法的精度。仿真所使用的主要参数如下:轨道高度为538.22 km;主星多普勒中心频率为-7.12 Hz;辅星多普勒中心频率为-75.31 Hz;脉冲重复频率为3 490 Hz;主星速度为7 656.55 m/s;波长为0.03 m;分辨率为3.20 m×1.61 m(方位向×距离向);地面场景大小为30 km×30 km,场景大地高范围4.22~397.78 m。首先对场景进行电磁建模,得到雷达回波数据;然后进行成像处理、干涉处理,得到基线参数和干涉相位数据。基线定标中主要误差源为控制点地面坐标误差、干涉相位误差、斜距误差[14],原始仿真数据无误差,为了与实际情况接近,分别加入误差。实际中控制点地面坐标的误差由图像上刺点误差和数字表面模型误差决定,一般在亚米级,因此给控制点地面坐标(3个坐标分量)分别加入随机误差${{\rm{e}}_1}\sim N\left({0, \sigma _1^2} \right)$,σ1=0.3 m。InSAR高程测量相位误差最大值在30°左右[15],因此干涉相位加入随机误差${e_2}\sim N\left({0, \sigma _2^2} \right)$,σ2=30°。斜距误差由设备内延迟误差和大气延迟误差组成,最大误差在3 m左右,因此斜距加入随机误差${e_3}\sim N\left({0, \sigma _3^2} \right)$,σ3=3 m。参考国外TanDEM系统定标前的基线误差为厘米级,同时存在少量的随机误差影响,因此给跨轨基线BX、顺轨基线BY、径向基线BZ分别增加了-5.00 cm、-5.00 cm、5.00 cm的常数误差(系统误差),并且给每个分量加入随机误差${e_4}\sim N\left({0, \sigma _4^2} \right)$,σ4=0.1 cm。

      由于主要误差源属于随机误差,因此定标结果会在真值附近随机波动,为了抑制随机误差的影响,需要进行多次定标实验。从精度评估角度来说,应该采用准确度作为衡量基线定标方法精度的指标,准确度指数学期望与真值的接近程度,常用偏差表示,这里准确度是指多次定标标出误差的数学期望与加入系统误差值的偏差。同时将标准差作为另外一个判据,标准差是指随机变量与其数学期望的接近程度,它表征算法受偶然误差影响的程度。标准差越小,说明单次定标结果与数学期望值越接近,算法的稳健性也越好,少量的定标实验即可满足精度要求;标准差越大,说明算法的稳健性较差,需要大量的定标实验方能满足定标要求。InSAR基线定标的主要目的是提高数字高程模型(digital elevation model,DEM)产品精度,因此可将基线定标后生成DEM精度与基线定标前生成DEM精度进行比较,来验证各种条件下定标方法的正确性。

    • 本实验设计了5种控制点布设方案,方位向×距离向分别为20(5×4)、60(10×6)、100(10×10)、140(14×10)、180(15×12)个控制点在场景中均匀分布。通过像坐标量测将控制点转刺到图像上,利用以上5种控制点布设方案进行基线定标实验。

      实验中,每种控制点布设方案都进行了200次定标,迭代终止条件设置为前后两次计算得到的每个轴向基线改正数之差的绝对值都小于0.1 mm,每次定标的迭代次数为2次,5种控制点布设方案的法方程系数矩阵条件数为3.75×104左右,统计均值μ和标准差σ记录在表 1中,同时计算出准确度的绝对值|δ|。在场景中均匀选取100(10×10)个地面点(与控制点不重合)作为检查点,验证定标前后DEM的精度,定标前由于各种误差源的存在,生成的DEM平面位置精度为55.22 m,高程精度为48.94 m,定标后DEM平面位置精度εp和高程精度εh记录在表 1中。

      表 1  控制点均匀分布的定标实验结果

      Table 1.  Results of Calibration Experiments with Evenly Distributed Control Points

      控制点数量/个 基线分量 系统误差/cm μ/cm σ/cm |δ|/cm εp/cm εh/cm
      20 BX -5.00 -4.48 7.95 0.52 355.31 75.23
      BY -5.00 -5.29 5.60 0.29
      BZ 5.00 5.46 6.99 0.46
      60 BX -5.00 -5.29 4.06 0.29 354.89 75.40
      BY -5.00 -5.10 3.34 0.10
      BZ 5.00 4.75 3.57 0.25
      100 BX -5.00 -4.77 2.81 0.23 355.19 75.23
      BY -5.00 -5.05 2.43 0.05
      BZ 5.00 5.20 2.47 0.20
      140 BX -5.00 -4.85 2.65 0.15 355.33 75.21
      BY -5.00 -4.88 2.13 0.12
      BZ 5.00 5.13 2.33 0.13
      180 BX -5.00 -5.05 2.25 0.05 355.44 75.26
      BY -5.00 -5.13 2.02 0.13
      BZ 5.00 4.95 1.98 0.05

      表 1中可以看出:(1)顺轨基线准确度绝对值随控制点数量的增加先减小后增加,其他两个方向基线标准差与准确度绝对值随控制点数量的增加而减小,标准差减小的幅度更大。(2)5种控制点布设方案基线三轴定标结果的准确度绝对值大多在毫米级,有的甚至达到亚毫米级,也就是多次定标标出误差的数学期望与误差设定值(真值)之间的偏差很小,说明了基于辅雷达距离改化方程和多普勒方程的基线定标模型是正确的,且具有较高的定标精度。(3)控制点数量对于定标精度具有重要影响。随着控制点数量增多,跨轨基线、径向基线定标准确度绝对值逐渐减小,顺轨基线定标准确度绝对值保持在1 mm左右;3个轴向基线的标准差逐渐减小,算法的稳健性逐渐提升。因此在实际定标中,应尽可能多布设控制点来提高定标精度和稳健性。(4)通过基线定标可显著提高DEM精度,DEM平面精度可提高一个数量级,DEM高程精度可提高两个数量级。由于存在其他误差源,基线定标后,生成的DEM平面精度为3.55 m左右,高程精度为0.75 m左右,不同控制点数量情况下,精度差别在毫米级。

    • 由于地面控制点获取途径不同,精度会存在较大的差异,因此给控制点地面坐标分别加入随机误差2.0 m、1.0 m、0.5 m、0.1 m后进行定标实验,分析控制点精度对定标结果影响,在场景中均匀选取60(10×6)个控制点,其他误差源加入的误差大小与§2.1中的相同,将结果记录在表 2中。

      表 2  不同控制点精度的定标实验结果

      Table 2.  Results of Calibration Experiments with Control Points of Different Accuracy

      控制点精度/m 基线分量 系统误差/cm μ/cm σ/cm |δ|/cm εp/cm εh/cm
      2.0 BX -5.00 -5.45 4.59 0.45 357.35 76.22
      BY -5.00 -5.65 22.38 0.65
      BZ 5.00 4.61 4.03 0.39
      1.0 BX -5.00 -5.36 4.14 0.36 354.88 75.42
      BY -5.00 -5.33 11.18 0.33
      BZ 5.00 4.69 3.64 0.31
      0.5 BX -5.00 -5.31 4.06 0.31 354.88 75.40
      BY -5.00 -5.16 5.58 0.16
      BZ 5.00 4.73 3.57 0.27
      0.1 BX -5.00 -5.27 4.08 0.27 354.89 75.39
      BY -5.00 -5.03 1.10 0.03
      BZ 5.00 4.76 3.59 0.24

      表 2可以看出,随着地面控制点精度提高,3个轴向基线的准确度绝对值和标准差逐渐减小,定标精度和算法的稳健性逐渐提升。利用定标后基线生成DEM,其平面精度和高程精度基本稳定,不同控制点精度情况下,精度差别在毫米级。因此在实际定标中,应尽可能采用高精度控制点。

    • 在场景中选取两个子条带,子条带的距离向宽为3 km,方位向长度与场景方位向长度一致,每个子条带布设30个控制点,控制点在子条带内均匀分布,共布设60个控制点。如图 3所示,三角形为控制点,矩形框为场景的边界,水平方向为方位向,竖直方向为距离向,虚线为子条带距离向边界。图 3(a)表示两个子条带位于距离向中部(不重叠);图 3(b)表示两个子条带分别位于场景距离向宽度1/3、2/3处;图 3(c)表示两个子条带分别位于近、远距端。利用以上3种控制点布设方式进行基线定标实验,各种误差源加入的误差大小与§2.1中的相同。

      图  3  控制点子条带位置示意图

      Figure 3.  Diagrams of Control Points Sub-band Location

      图 4为3种控制点布设方式基线定标结果统计示意图,每个正方形表示一次定标标出的基线分量误差值。图 5为基线三轴标出值与设定值互差示意图。从图 4图 5中可以看出,随着控制点子条带在距离向间距变大,标出三轴误差的均值与系统误差设定值之间的偏差变小,单次标出的跨轨基线、径向基线误差值与误差均值的偏差也变小,单次标出的顺轨基线误差值与误差均值的偏差变化量很小。

      图  4  控制点子条带不同位置时定标误差统计

      Figure 4.  Calibration Error Statistics for Sub-bands of Control Points at Different Positions

      图  5  控制点子条带不同位置时标出值与设定值互差

      Figure 5.  Difference Between the Calibrated Value and the Set Value for Sub-bands of Control Points at Different Positions

      将定标实验的均值μ和标准差σ记录在表 3中,同时计算出准确度的绝对值|δ|,可见标准差与准确度绝对值随着控制点子条带在距离向间距变大而减小,跨轨基线、径向基线的标准差与准确度绝对值减小幅度较顺轨基线的大。利用定标后基线生成DEM,其平面精度和高程精度基本稳定,不同控制点子条带位置情况下,精度差别在毫米级。

      表 3  控制点子条带不同位置的定标结果

      Table 3.  Calibration Results for Sub-bands of Control Points at Different Positions/cm

      子条带位置 基线分量 系统误差 μ σ |δ| εp εh
      距离向中部 BX -5.00 -3.79 22.26 1.21 356.19 75.37
      BY -5.00 -4.74 3.57 0.26
      BZ 5.00 6.06 19.58 1.06
      距离向宽度1/3、2/3处 BX -5.00 -5.35 6.35 0.35 355.04 75.38
      BY -5.00 -4.87 3.44 0.13
      BZ 5.00 4.69 5.58 0.31
      近、远距端 BX -5.00 -5.17 2.28 0.17 355.12 75.30
      BY -5.00 -4.87 3.44 0.13
      BZ 5.00 4.85 2.00 0.15

      比较表 1表 3可知,60个控制点两侧布设于近、远距端方式的定标结果优于均匀布设方式,接近140个控制点均匀布设方式的定标结果。由此可以推出,若要达到相同的定标精度,前者布设方式所需的控制点数量明显少于后者。控制点布设策略对于定标精度同样具有重要影响,控制点在近、远距端两侧布设的方式能显著提高定标精度和稳健性,这与理论分析一致,是一种较好的布设策略。

    • 本文提出了一种星载分布式InSAR基线定标新方法,该方法基于辅雷达距离改化方程和多普勒方程的定标模型,可采用均匀或近、远距端两侧布设控制点的策略,在条件允许的情况下优先选择后一种策略,同时应尽可能采用高精度地面控制点,通过最小二乘平差迭代计算确定基线误差,基线定标精度能达到毫米级。本文提出的方法与国外TanDEM系统中的定标方法相比,精度相当,增加了对顺轨基线的定标,从而可实现三维基线定标,所需场地面积较小,一景数据范围(约30 km×30 km)即可满足要求,控制点的布设方式相对灵活,兼顾了定标精度与地面控制条件,因此具有一定优势。本文提出的方法可为星载分布式InSAR基线定标提供一种新途径。

参考文献 (15)

目录

    /

    返回文章
    返回