留言板

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

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

利用高分三号影像监测依嘎冰川表面运动

王群 张蕴灵 范景辉 傅宇浩

王群, 张蕴灵, 范景辉, 傅宇浩. 利用高分三号影像监测依嘎冰川表面运动[J]. 武汉大学学报 ● 信息科学版, 2020, 45(3): 460-466. doi: 10.13203/j.whugis20190036
引用本文: 王群, 张蕴灵, 范景辉, 傅宇浩. 利用高分三号影像监测依嘎冰川表面运动[J]. 武汉大学学报 ● 信息科学版, 2020, 45(3): 460-466. doi: 10.13203/j.whugis20190036
WANG Qun, ZHANG Yunling, FAN Jinghui, FU Yuhao. Monitoring the Motion of the Yiga Glacier Using GF-3 Images[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 460-466. doi: 10.13203/j.whugis20190036
Citation: WANG Qun, ZHANG Yunling, FAN Jinghui, FU Yuhao. Monitoring the Motion of the Yiga Glacier Using GF-3 Images[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 460-466. doi: 10.13203/j.whugis20190036

利用高分三号影像监测依嘎冰川表面运动

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

国家重大科技专项 07-Y30B10-9001-14/16

中欧科技合作“龙计划”四期 32365

中国地质调查局地质调查项目 DD20190515

详细信息
    作者简介:

    王群, 硕士, 助理工程师, 主要从事InSAR技术应用与研究。wangq0723@163.com

    通讯作者: 范景辉, 博士, 教授级高级工程师。jhfan2004@qq.com
  • 中图分类号: P237

Monitoring the Motion of the Yiga Glacier Using GF-3 Images

Funds: 

The National Science and Technology Major Project of China 07-Y30B10-9001-14/16

ESA-MOST China Dragon 4 Cooperation 32365

Geological Survey Program of China Geological Survey DD20190515

More Information
    Author Bio:

    WANG Qun, master, assistant engineer, specializes in interferometric synthetic aperture radar. E-mail: wangq0723@163.com

    Corresponding author: FAN Jinghui, PhD, professor. E-mail: jhfan2004@qq.com
  • 摘要: 冰川运动是冰川及冰川湖危险性评估的重要依据。高分三号作为中国首颗民用星载合成孔径雷达卫星,因其全天时、全天候以及高空间分辨率的特点在冰川运动监测中有着重要优势。利用5景精细条带2成像模式的高分三号影像,应用偏移量跟踪技术监测了念青唐古拉山依嘎冰川在2017-11—2018-03期间4个时段的运动分布情况。通过与高分三号影像时间分布类似的两景Sentinel-1影像的偏移量跟踪结果进行对比,并基于研究区基岩运动速度应为0的假设,计算了高分三号各个时段结果基岩区的速度残差,评估了高分三号影像在冰川表面运动监测中的适用性。高分三号监测结果表明,4个时段依嘎冰川表面具有相似的流速分布,其表面速度最大处均分布在冰川表面高程变化最剧烈的中部,且与基于两景Sentinel-1影像的偏移量跟踪结果在总体强弱态势上吻合;各个时段高分三号监测结果基岩区速度残差的均方根误差分别为1.4 cm/d、2.0 cm/d、1.7 cm/d、2.3 cm/d,验证了利用高分三号影像监测依嘎冰川表面运动的可靠性。以上分析结果表明,高分三号影像可作为冰川表面运动监测的常规数据来源之一。由于高分三号同时具有高空间分辨率和高性价比的特点,因此它能为中国冰川遥感监测发挥独特作用。
  • 图  1  依嘎冰川地理位置

    Figure  1.  Geographic Location of the Yiga Glacier

    图  2  依嘎冰川水平面二维速率图

    Figure  2.  2D Velocity on Horizontal Plane of the Yiga Glacier

    图  3  高分三号影像依嘎冰川监测结果剖面线速率及高程变化图

    Figure  3.  Velocity Profiles of GF-3 Monitoring Results of the Yiga Glacier and Elevation Variation Along the Profile

    图  4  高分三号(2017-12-16—2018-01-14)与Sentinel-1(2017-12-10—2018-01-15)影像监测结果剖面线速率图

    Figure  4.  Velocity Profiles of Monitoring Results of GF-3(2017-12-16—2018-01-14)and Sentinel-1(2017-12-10—2018-01-15) Images

    表  1  本文使用的高分三号数据主要参数

    Table  1.   Main Parameters of GF-3 Images Used in This Paper

    主影像 从影像 垂直基线/m 时间基线/d
    2017-11-17 2017-12-16 -1 766.8 29
    2017-12-16 2018-01-14 2 195.8 29
    2018-01-14 2018-02-12 -1 177.1 29
    2018-02-12 2018-03-13 412.8 29
    下载: 导出CSV

    表  2  本文使用的Sentinel-1数据主要参数

    Table  2.   Main Parameters of Sentinel-1 Images Used in This Paper

    主影像 从影像 垂直基线/m 时间基线/d
    2017-12-10 2018-01-15 37.1 36
    下载: 导出CSV

    表  3  研究区基岩区速率残差的均方根误差/(cm∙d-1)

    Table  3.   The RMSE of Velocity in the Bedrock Region of Offset Tracking Techniques/(cm∙d-1)

    影像对 A B C 均值
    2017-11-17—2017-12-16 1.1 1.7 1.3 1.4
    2017-12-16—2018-01-14 1.8 2.3 1.8 2.0
    2018-01-14—2018-02-12 1.4 1.5 2.2 1.7
    2018-02-12—2018-03-13 1.8 3.1 2.0 2.3
    2017-12-10—2018-01-15 1.6 1.8 1.3 1.6
    下载: 导出CSV
  • [1] Vaughan D G, Comiso J, Allison I, et al. Observations:Cryosphere, Climate Change 2013:The Physical Science Basis[M]. New York:Cambridge University Press, 2013
    [2] 姚晓军, 刘时银, 孙美平, 等. 20世纪以来西藏冰湖溃决灾害事件梳理[J].自然资源学报, 2014, 29(8):1377-1390 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zrzyxb201408010

    Yao Xiaojun, Liu Shiyin, Sun Meiping, et al. Study on the Glacial Lake Outburst Flood Events in Tibet Since the 20th Century[J]. Journal of Natural Resources, 2014, 29(8):1377-1390 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zrzyxb201408010
    [3] 谢自楚, 刘潮海.冰川学导论[M].上海:上海科学普及出版社, 2010

    Xie Zichu, Liu Chaohai. Introduction to Glaciology[M]. Shanghai:Shanghai Popular Science Press, 2010
    [4] Euillades L D, Euillades P A, Riveros N C, et al. Detection of Glaciers Displacement Time-Series Using SAR[J]. Remote Sensing of Environment, 2016, 184:188-198 doi:  10.1016/j.rse.2016.07.003
    [5] Mouginot J, Scheuchl B, Rignot E. Mapping of Ice Motion in Antarctica Using Synthetic-Aperture Radar Data[J].Remote Sensing, 2012, 4(12):2753-2767 http://d.old.wanfangdata.com.cn/OAPaper/oai_doaj-articles_15d8440d5c23ce5ff6c9d8049b9838bf
    [6] Yan S, Run Z, Liu G, et al. Deriving Ice Motion Patterns in Mountainous Regions by Integrating the Intensity-Based Pixel-Tracking and Phase-Based D-InSAR and MAI Approaches:A Case Study of the Chongce Glacier[J]. Remote Sensing, 2016, 8(7):611-624 doi:  10.3390/rs8070611
    [7] Goldstein R M, Engelhardt H, Kamb B, et al.Satellite Radar Interferometry for Monitoring Ice Sheet Motion:Application to an Antarctic Ice Stream[J]. Science, 1993, 262(5139):1525-1530 doi:  10.1126/science.262.5139.1525
    [8] Gray A L, Mattar K E, Vchon P W, et al. InSAR Results from the RADARSAT Antarctic Mapping Mission Data: Estimation of Glacier Motion Using a Simple Registration Procedure[C]. IEEE International Geoscience and Remote Sensing Symposium, Seattle, WA, USA, 1998
    [9] Luckman A, Quincey D, Bevan S. The Potential of Satellite Radar Interferometry and Feature Tracking for Monitoring Flow Rates of Himalayan Glaciers[J].Remote Sensing of Environment, 2007, 111(2/3):172-181 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=7a56657b93478efa36e69c8085b48df7
    [10] 周春霞, 邓方慧, 艾松涛, 等.利用DInSAR的东南极极记录和达尔克冰川冰流速提取与分析[J].武汉大学学报·信息科学版, 2014, 39(8):940-944 http://ch.whu.edu.cn/CN/abstract/abstract3048.shtml

    Zhou Chunxia, Deng Fanhui, Ai Songtao, et al. Determination of Ice-Flow Velocity at the Polar Record Glacier and Dalk Glacier Using DInSAR[J]. Geomatics and Information Science of Wuhan University, 2014, 39(8):940-944 http://ch.whu.edu.cn/CN/abstract/abstract3048.shtml
    [11] Li J, Li Z W, Wu L X, et al. Deriving a Time Series of 3D Glacier Motion to Investigate Interactions of a Large Mountain Glacial System with Its Glacial Lake:Use of Synthetic Aperture Radar Pixel Offset-Small Baseline Subset Technique[J]. Journal of Hydrology, 2018, 559:596-608 doi:  10.1016/j.jhydrol.2018.02.067
    [12] Mouginot J, Rignot E, Scheuchl B, et al. Comprehensive Annual Ice Sheet Velocity Mapping Using Landsat-8, Sentinel-1, and RADARSAT-2 Data[J]. Remote Sensing, 2017, 9(4):364-383 doi:  10.3390/rs9040364
    [13] 王思胜, 江利明, 孙永玲, 等.基于ALOS PALSAR数据的山地冰川流速估算方法比较:以喀喇昆仑地区斯克洋坎力冰川为例[J].国土资源遥感, 2016, 28(2):54-61 http://d.old.wanfangdata.com.cn/Periodical/chxb200904010

    Wang Sisheng, Jiang Liming, Sun Yongling, et al. Evaluation of Methods for Deriving Mountain Glacier Velocities with ALOS PALSAR Images:A Case Study of SKyang Glacier in Central Karakoram[J]. Remote Sensing for Land and Resources, 2016, 28(2):54-61 http://d.old.wanfangdata.com.cn/Periodical/chxb200904010
    [14] 张晓博, 赵学胜, 葛大庆, 等.利用新型C波段雷达卫星研究南伊内里切克冰川运动特征[J].武汉大学学报·信息科学版, 2019, 44(3):429-435 http://ch.whu.edu.cn/CN/abstract/abstract6385.shtml

    Zhang Xiaobo, Zhao Xuesheng, Ge Daqing, et al. Motion Characteristics of the South Inilchek Glacier Derived from New C-Band SAR Satellite[J]. Geomatics and Information Science of Wuhan University, 2019, 44(3):429-435 http://ch.whu.edu.cn/CN/abstract/abstract6385.shtml
    [15] 张庆君.高分三号卫星总体设计与关键技术[J].测绘学报, 2017, 46(3):269-277 http://d.old.wanfangdata.com.cn/Periodical/chxb201703002

    Zhang Qingjun. System Design and Key Technologies of the GF-3 Satellite[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(3):269-277 http://d.old.wanfangdata.com.cn/Periodical/chxb201703002
    [16] 刘时银, 郭万钦, 许君利, 等.中国第二次冰川编目数据集(版本1.0)[R].寒区旱区科学数据中心, 兰州, 中国, 2014

    Liu Shiyin, Guo Wanqing, Xu Junli, et al.The Second Glacier Inventory Dataset of China (Version 1.0)[R]. Cold and Arid Regions Science Data Center, Lanzhou, China, 2014
    [17] Wang H, Yang J, Mouche A, et al. GF-3 SAR Ocean Wind Retrieval:The First View and Preliminary Assessment[J]. Remote Sensing, 2017, 9(7):694-705 doi:  10.3390/rs9070694
    [18] Grohmann C H. Evaluation of TanDEM-X DEMs on Selected Brazilian Sites:Comparison with SRTM, ASTER GDEM and ALOS AW3D30[J]. Remote Sensing of Environment, 2018, 212:121-133 doi:  10.1016/j.rse.2018.04.043
    [19] Strozzi T, Luckman A, Murray T, et al. Glacier Motion Estimation Using SAR Offset-Tracking Procedures[J]. IEEE Transactions on Geoscience and Remote Sensing, 2002, 40(11):2384-2391 doi:  10.1109/TGRS.2002.805079
    [20] Haemmig C, Huss M, Keusen H, et al. Hazard Assessment of Glacial Lake Outburst Floods from Kyagar Glacier, Karakoram Mountains, China[J]. Annals of Glaciology, 2014, 55(66):34-44 doi:  10.3189/2014AoG66A001
    [21] Li J, Li Z, Zhu J, et al. Deriving Surface Motion of Mountain Glaciers in the Tuomuer-Khan Tengri Mountain Ranges from PALSAR Images[J].Global and Planetary Change, 2013, 101:61-71 doi:  10.1016/j.gloplacha.2012.12.004
    [22] 邓方慧, 周春霞, 王泽民, et al.利用偏移量跟踪测定Amery冰架冰流汇合区的冰流速[J].武汉大学学报·信息科学版, 2015, 40(7):901-906 http://ch.whu.edu.cn/CN/abstract/abstract3296.shtml

    Deng Fanghui, Zhou Chunxia, Wang Zemin, et al. Ice-Flow Velocity Derivation of the Confluence Zone of the Amery Ice Shelf Using Offset-Tracking Method[J]. Geomatics and Information Science of Wuhan University, 2015, 40(7):901-906 http://ch.whu.edu.cn/CN/abstract/abstract3296.shtml
    [23] Wegmueller U, Werner C, Strozzi T, et al. Ionospheric Electron Concentration Effects on SAR and InSAR[C]. IEEE International Geoscience and Remote Sensing Symposium, Denver, CO, USA, 2006
    [24] GAMMA Remote Sensing. Differential Interferometry and Geocoding Software: DIFF and GEO(Geocoding and Image Registration)[R]. GAMMA Remote Sensing AG, Bern, Switzerland, 2008
    [25] 王群, 范景辉, 周伟, 等. DEM辅助偏移量跟踪技术的山地冰川运动监测研究[J].国土资源遥感, 2018, 30(3):167-173 http://d.old.wanfangdata.com.cn/Periodical/gtzyyg201803023

    Wang Qun, Fan Jinghui, Zhou Wei, et al. Research on the DEM-Assisted Offset Tracking Technique Applied to Glaciers Movement Monitoring[J]. Remote Sensing for Land and Resources, 2018, 30(3):167-173 http://d.old.wanfangdata.com.cn/Periodical/gtzyyg201803023
    [26] Wang Q, Fan J, Zhou W, et al. 3D Surface Velocity Retrieval of Mountain Glacier Using an Offset Tracking Technique Applied to Ascending and Descending SAR Constellation Data:A Case Study of the Yiga Glacier[J]. International Journal of Digital Earth, 2019, 12(6):614-624 doi:  10.1080/17538947.2018.1470690
    [27] 李佳, 李志伟, 汪长城, 等. SAR偏移量跟踪技术估计天山南依内里切克冰川运动[J].地球物理学报, 2013, 56(4):1226-1236 http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201304017

    Li Jia, Li Zhiwei, Wang Changcheng, et al. Using SAR Offset Tracking Approach to Estimate Surface Motion of the South Inylchek Glacier in Tianshan[J].Chinese Journal of Geophysics, 2013, 56(4):1226-1236 http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201304017
    [28] 闫世勇.山地冰川表面运动雷达遥感监测方法研究[D].北京: 中国科学院大学, 2013 http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y2382133

    Yan Shiyong. Research on Extraction of Alpine Glacier Surface Movement by SAR Remote Sensing[D]. Beijing: University of Chinese Academy of Sciences, 2013 http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y2382133
  • [1] 周中正, 许才军, 刘洋, 熊露雲.  岗纳楼冰川表面流速时空变化特征提取及分析 . 武汉大学学报 ● 信息科学版, 2022, 47(2): 226-233. doi: 10.13203/j.whugis20200133
    [2] 熊俊麟, 范宣梅, 窦向阳, 杨莹辉.  藏东南然乌湖流域雅弄冰川流速季节性变化 . 武汉大学学报 ● 信息科学版, 2021, 46(10): 1579-1588. doi: 10.13203/j.whugis20210104
    [3] 柳林, 宋豪峰, 杜亚男, 冯光财, 刘清瑶, 孙敏.  联合哨兵2号和Landsat 8估计白格滑坡时序偏移量 . 武汉大学学报 ● 信息科学版, 2021, 46(10): 1461-1470. doi: 10.13203/j.whugis20200596
    [4] 王晓文, 刘巧, 张波, 张瑞, 刘国祥.  近40 a昆仑山口37号冰川冰舌滑塌多源遥感监测与分析 . 武汉大学学报 ● 信息科学版, 2020, 45(11): 1687-1696. doi: 10.13203/j.whugis20200214
    [5] 张波, 张瑞, 刘国祥, 刘巧, 蔡嘉伦, 于冰, 符茵, 李志林.  基于SAR影像的贡巴冰川末端冰湖年际变化监测及溃决规律分析 . 武汉大学学报 ● 信息科学版, 2019, 44(7): 1054-1064. doi: 10.13203/j.whugis20190087
    [6] 张晓博, 赵学胜, 葛大庆, 刘斌.  利用新型C波段雷达卫星研究南伊内里切克冰川运动特征 . 武汉大学学报 ● 信息科学版, 2019, 44(3): 429-435. doi: 10.13203/j.whugis20160538
    [7] 刘国祥, 张波, 张瑞, 蔡嘉伦, 符茵, 刘巧, 于冰, 李志林.  联合卫星SAR和地基SAR的海螺沟冰川动态变化及次生滑坡灾害监测 . 武汉大学学报 ● 信息科学版, 2019, 44(7): 980-995. doi: 10.13203/j.whugis20190077
    [8] 王密, 杨博, 李德仁, 龚健雅, 皮英冬.  资源三号全国无控制整体区域网平差关键技术及应用 . 武汉大学学报 ● 信息科学版, 2017, 42(4): 427-433. doi: 10.13203/j.whugis20160534
    [9] 郭丽, 李金岭, 童锋贤, 王广利, 黄飞, 刘庆会, 郑鑫, 郑为民.  同波束VLBI技术对嫦娥三号巡视器的高精度相对定位 . 武汉大学学报 ● 信息科学版, 2016, 41(8): 1125-1130. doi: 10.13203/j.whugis20140439
    [10] 孙维君, 闫明, 艾松涛, 朱国才, 王泽民, 刘雷保, 徐跃通, 任贾文.  北极新奥尔松地区Austre Lovénbreen冰川温度变化特征 . 武汉大学学报 ● 信息科学版, 2016, 41(1): 79-85. doi: 10.13203/j.whugis20150302
    [11] 邓方慧, 周春霞, 王泽民, 鄂栋臣, 张辛.  利用偏移量跟踪测定Amery冰架冰流汇合区的冰流速 . 武汉大学学报 ● 信息科学版, 2015, 40(7): 901-906. doi: 10.13203/j.whugis20130653
    [12] 周春霞, 邓方慧, 艾松涛, 王泽民, 鄂栋臣.  利用DInSAR的东南极极记录和达尔克冰川冰流速提取与分析 . 武汉大学学报 ● 信息科学版, 2014, 39(8): 940-944. doi: 10.13203/j.whugis20130116
    [13] 江敏, 钟敏, 张子占, 冯伟.  冰川均衡调整对全球海水质量变化的影响 . 武汉大学学报 ● 信息科学版, 2013, 38(10): 1188-1191.
    [14] 艾松涛, 王泽民, 鄂栋臣, 闫明.  利用GPS的北极冰川运动监测与分析 . 武汉大学学报 ● 信息科学版, 2012, 37(11): 1337-1340.
    [15] 董晓燕, 丁晓利, 李志伟, 王永哲.  一种新的SAR像素偏移量估计流程及其在同震形变监测中的应用 . 武汉大学学报 ● 信息科学版, 2011, 36(7): 789-792.
    [16] 肖金群, 汪长城, 李志伟, 杨亚夫.  基于配准偏移量的InSAR基线估计 . 武汉大学学报 ● 信息科学版, 2010, 35(10): 1236-1239.
    [17] 张小红, 鄂栋臣.  用PPP技术确定南极Amery冰架的三维运动速度 . 武汉大学学报 ● 信息科学版, 2005, 30(10): 909-912.
    [18] 张世强, 卢健, 刘时银.  利用TM高光谱图像提取青藏高原喀喇昆仑山区现代冰川边界 . 武汉大学学报 ● 信息科学版, 2001, 26(5): 435-440.
    [19] 王清华, 鄂栋臣, 陈春明.  中山站至A冰穹考察及沿线GPS复测结果分析 . 武汉大学学报 ● 信息科学版, 2001, 26(3): 200-204,231.
    [20] 陈宗信, 虞家钧, 汤明宝.  对河谷、冰川、黄土三种地貌表示法的研究 . 武汉大学学报 ● 信息科学版, 1959, 3(3): 58-71.
  • 加载中
图(4) / 表(3)
计量
  • 文章访问数:  1694
  • HTML全文浏览量:  285
  • PDF下载量:  131
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-18
  • 刊出日期:  2020-03-05

利用高分三号影像监测依嘎冰川表面运动

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

    国家重大科技专项 07-Y30B10-9001-14/16

    中欧科技合作“龙计划”四期 32365

    中国地质调查局地质调查项目 DD20190515

    作者简介:

    王群, 硕士, 助理工程师, 主要从事InSAR技术应用与研究。wangq0723@163.com

    通讯作者: 范景辉, 博士, 教授级高级工程师。jhfan2004@qq.com
  • 中图分类号: P237

摘要: 冰川运动是冰川及冰川湖危险性评估的重要依据。高分三号作为中国首颗民用星载合成孔径雷达卫星,因其全天时、全天候以及高空间分辨率的特点在冰川运动监测中有着重要优势。利用5景精细条带2成像模式的高分三号影像,应用偏移量跟踪技术监测了念青唐古拉山依嘎冰川在2017-11—2018-03期间4个时段的运动分布情况。通过与高分三号影像时间分布类似的两景Sentinel-1影像的偏移量跟踪结果进行对比,并基于研究区基岩运动速度应为0的假设,计算了高分三号各个时段结果基岩区的速度残差,评估了高分三号影像在冰川表面运动监测中的适用性。高分三号监测结果表明,4个时段依嘎冰川表面具有相似的流速分布,其表面速度最大处均分布在冰川表面高程变化最剧烈的中部,且与基于两景Sentinel-1影像的偏移量跟踪结果在总体强弱态势上吻合;各个时段高分三号监测结果基岩区速度残差的均方根误差分别为1.4 cm/d、2.0 cm/d、1.7 cm/d、2.3 cm/d,验证了利用高分三号影像监测依嘎冰川表面运动的可靠性。以上分析结果表明,高分三号影像可作为冰川表面运动监测的常规数据来源之一。由于高分三号同时具有高空间分辨率和高性价比的特点,因此它能为中国冰川遥感监测发挥独特作用。

English Abstract

王群, 张蕴灵, 范景辉, 傅宇浩. 利用高分三号影像监测依嘎冰川表面运动[J]. 武汉大学学报 ● 信息科学版, 2020, 45(3): 460-466. doi: 10.13203/j.whugis20190036
引用本文: 王群, 张蕴灵, 范景辉, 傅宇浩. 利用高分三号影像监测依嘎冰川表面运动[J]. 武汉大学学报 ● 信息科学版, 2020, 45(3): 460-466. doi: 10.13203/j.whugis20190036
WANG Qun, ZHANG Yunling, FAN Jinghui, FU Yuhao. Monitoring the Motion of the Yiga Glacier Using GF-3 Images[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 460-466. doi: 10.13203/j.whugis20190036
Citation: WANG Qun, ZHANG Yunling, FAN Jinghui, FU Yuhao. Monitoring the Motion of the Yiga Glacier Using GF-3 Images[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 460-466. doi: 10.13203/j.whugis20190036
  • 受全球气候变暖影响,中国大部分冰川呈现出加速退缩趋势[1],且近年来冰崩、冰川湖溃决洪水等冰川灾害时有发生,规模和频率都有所增加[2]。冰川运动作为冰川的特征之一,是研究冰川物质平衡和全球气候变化的重要参数之一,也是冰川及冰川湖危险性评估的重要依据,监测冰川运动变化对研究气候变化和冰川灾害预警有重要意义[3]

    目前,现场实地监测和遥感监测等多种技术手段都可用于监测冰川运动速度的变化[4-5]。但是冰川往往分布在人迹罕至、环境恶劣的地区,对实地观测人员的生命安全有着极大的威胁,并且实地监测费时费力,监测范围和频次受限[6]。随着遥感技术的发展,特别是星载合成孔径雷达(synthetic aperture radar,SAR)的出现,冰川区的SAR影像可以全天时、全天候获取,有效克服一些冰川区多云雨天气条件的限制,提高影像获取频率;还可以根据SAR影像的强度信息和相位信息,应用偏移量跟踪技术、差分干涉雷达测量技术(differential interferometric synthetic aperture radar,D-InSAR)、多孔径干涉测量技术(multiple aperture interferometry,MAI)等多种技术手段获取大范围、高精度的冰川表面运动情况,成为冰川运动监测的重要手段之一[7-14]

    目前,可用于监测冰川表面运动的主要在轨SAR卫星包括Sentinel-1、TerraSAR-X、COSMO-SkyMed星座、RadarSAT-2等,国内外学者也利用以上SAR卫星数据在冰川表面运动监测研究中取得了较多应用成果。高分三号(GF-3)卫星是首颗中国自主研制的C频段多极化SAR卫星,具有高分辨率、大成像幅宽、多成像模式等特点,与Sentinel-1等免费数据相比具有高分辨率的优势,与TerraSAR-X、COSMO-SkyMed、RadarSAR-2等商业数据相比具有性价比高的优势,其在海洋、气象、水利、减灾等多个领域均有成功应用案例[15]。但利用高分三号影像监测冰川表面运动的研究尚属空白,故本文以发掘高分三号在中国冰川遥感监测中的作用为目标,利用覆盖念青唐古拉山依嘎冰川的5景高分三号影像,采用基于SAR影像强度信息的偏移量跟踪技术,提取了依嘎冰川表面运动分布情况,分析了其运动特征,并与Sentinel-1及前人研究结果对比,评估了高分三号监测冰川表面运动的适用性。

    • 依嘎冰川位于中国西藏自治区那曲地区嘉黎县念青唐古拉山脉西段,地处易贡藏布一级支流尼都藏布上游北侧,具体位置和范围见图 1。根据2014年中国科学院寒区旱区环境与工程研究所发布的《中国第二次冰川编目》,依嘎冰川面积约为46.3 km2,高差接近3 000 m[16]。依嘎冰川末端与尼都藏布相接,其运动极可能堵塞河道,形成堰塞湖,威胁下游众多村落的安全,监测其运动情况具有重要的意义。

      图  1  依嘎冰川地理位置

      Figure 1.  Geographic Location of the Yiga Glacier

    • 高分三号卫星于2016-08-10在中国太原卫星发射基地发射,是国家高分辨率对地观测系统重大专项中的唯一一颗民用微波遥感卫星,同时也是中国自主研制的首颗C波段多极化卫星,重复周期29 d[17]。本文收集了研究区2017-11—2018-03之间的5景升轨高分三号影像,成像模式为精细条带2模式,斜距向分辨率为2.25 m,方位向分辨率为4.78 m,极化方式为HH极化,使用的高分三号数据主要参数见表 1;同时收集了与本文高分三号影像时间分布类似的两景升轨Sentinel-1数据,成像模式为方位波束扫描(terrain observation by progressive scans,TOPS)模式,斜距向分辨率为2.33 m,方位向分辨率为13.96 m,极化方式为VV极化,使用的Sentinel-1数据主要参数见表 2;采用30 m分辨率的AW3D30数字表面模型(digital surface model,DSM)数据辅助SAR影像配准和地理编码,其高程精度优于5 m[18]

      表 1  本文使用的高分三号数据主要参数

      Table 1.  Main Parameters of GF-3 Images Used in This Paper

      主影像 从影像 垂直基线/m 时间基线/d
      2017-11-17 2017-12-16 -1 766.8 29
      2017-12-16 2018-01-14 2 195.8 29
      2018-01-14 2018-02-12 -1 177.1 29
      2018-02-12 2018-03-13 412.8 29

      表 2  本文使用的Sentinel-1数据主要参数

      Table 2.  Main Parameters of Sentinel-1 Images Used in This Paper

      主影像 从影像 垂直基线/m 时间基线/d
      2017-12-10 2018-01-15 37.1 36
    • 本文使用偏移量跟踪技术对覆盖依嘎冰川的高分三号和Sentinel-1影像进行处理。偏移量跟踪技术可以计算两景SAR影像获取的时间段内冰川在距离向和方位向的运动分布情况[19]。一般来说,偏移量跟踪技术的精度可以达到SAR影像像元分辨率的1/10以上[20],以10 m分辨率的高分三号数据为例,其结果精度优于1 m。

      偏移量跟踪技术的关键算法是归一化互相关算法,归一化互相关算法计算出的SAR影像总偏移量一般是地形、电离层、轨道以及冰川运动引起的偏移量的总和[9, 21-22]。在研究区地形起伏不大且SAR影像对空间基线较短时,地形起伏引起的偏移量往往可以忽略[21],而电离层引起的偏移量与研究区的纬度和SAR影像的波长相关[23]。本文研究区处于低纬度区域,SAR数据为C波段,且与电离层变化的空间尺度相比冰川面积很小,电离层影响可以忽略。但研究区地形陡峭,影像对的平均垂直基线在1 000 m以上,地形引起的偏移量不可忽略。故本文使用基于外部数字高程模型(digital elevation model,DEM)辅助的SAR影像配准方法,该方法利用外部DEM和SAR影像的轨道信息建立起主从影像的配准查找表,将研究区地形引起的偏移量计算到查找表内,并且对根据初始查找表重采样获得的主影像与从影像实施互相关配准,据此进一步精化查找表,降低由轨道定位不准确引起的偏移量误差,在空间基线较大时仍能精确地配准两幅SAR影像,提高地形起伏地区偏移量跟踪结果的精度[24-25]

    • 本文使用的5景高分三号影像,除2017-11-17和2018-03-13两个时间段接收的影像具有相似的覆盖范围外,其余3景SAR影像的覆盖范围均有较大差别。故首先依据依嘎冰川在SAR影像中的具体位置,将5景SAR影像适当裁剪为同一覆盖范围,然后利用外部DEM辅助SAR影像配准。配准完成后,设定搜索窗口对SAR影像执行基于影像强度信息的互相关计算,互相关系数最大处即为像元的偏移量。本文将搜索窗口设为150×150像元,可以较好地保留结果的局部细节。对Sentinel-1数据的搜索窗口设为256×64像元。

    • 对高分三号以及Sentinel-1影像应用偏移量跟踪技术获取了依嘎冰川地距向和方位向的位移分布,将地距向和方位向位移合成并转化为冰川日均流速,从而得到了两个数据集的依嘎冰川水平面二维速率场,如图 2所示。图 2中,GF-3表示高分三号影像,S1表示Sentinel-1影像。

      图  2  依嘎冰川水平面二维速率图

      Figure 2.  2D Velocity on Horizontal Plane of the Yiga Glacier

      根据高分三号监测结果图(图 2(a)~2(d))可知:(1)整体来看,在可监测到的冰川运动范围内,冰川末端到冰川上端,冰川表面速度呈现出先减速、后增加到基本稳定,再到显著增加、显著减少至稳定的趋势。其中,冰川最大速度主要分布在冰川中部区域,在4个监测时段内,中部冰川区的运动速率均大于1 m/d。(2)由于缺少明显的地物特征(如裂隙、冰碛物等),依嘎冰川积累区表面较为平滑,偏移量跟踪技术适用性降低,导致大部分积累区无法监测出形变结果。(3)与近似同时期Sentinel-1影像结果相比(图 2(e)),两组冰川流速结果的强弱态势在空间上高度吻合,最大流速均分布在冰川中部,但在细节上略有不同。相比15 m空间分辨率的Sentinel-1流速结果,5 m分辨率的高分三号流速监测结果细节更加丰富。(4)在4个监测时段内,依嘎冰川表面运动分布具有高度一致性,与文献[26]利用COSMO-SkyMed反演出的依嘎冰川运动结果进行对比,在冰川中部坡度陡变处的流速显著增加、显著减小变化情况一致,两种结果显示出的冰川运动分布总体吻合。

      为定量化分析依嘎冰川的速率变化情况,沿冰川主冰流线取一条剖面线(图 2(a)中黑色曲线,起点为冰川末端),统计了沿剖面线高分三号结果在4个时段的冰川速率变化及高程变化(见图 3);将高分三号监测结果欠采样至Sentinel-1相同的分辨率,并对比剖面线上高分三号与Sentinel-1的速率分布(见图 4)。

      图  3  高分三号影像依嘎冰川监测结果剖面线速率及高程变化图

      Figure 3.  Velocity Profiles of GF-3 Monitoring Results of the Yiga Glacier and Elevation Variation Along the Profile

      图  4  高分三号(2017-12-16—2018-01-14)与Sentinel-1(2017-12-10—2018-01-15)影像监测结果剖面线速率图

      Figure 4.  Velocity Profiles of Monitoring Results of GF-3(2017-12-16—2018-01-14)and Sentinel-1(2017-12-10—2018-01-15) Images

      图 3可以看出,剖面线上4个时段依嘎冰川具有相似的速率分布,经统计,各个时段剖面线上冰川的平均运动速率依次为50 cm/d、49 cm/d、48 cm/d和49 cm/d,十分接近,表明2017-11—2018-03依嘎冰川运动速率稳定;与高程联合分析可知,距剖面线起点的4 000~4800 m是冰川表面高程变化最剧烈的区域,此处冰川运动速率急剧增加,到5 500 m处,高程变化稳定,冰川速率也下降并趋于稳定。以上表明,地面高程变化与冰川速率的变化密切相关。在距剖面线起点约600 m处,速率出现陡降,除与局部坡度变缓有关外,也与此处山谷走向变化、冰流受山体阻挡有关。由图 4可知,高分三号和Sentinel-1剖面线上速率整体变化态势较为吻合,在距剖面线起点约4 500 m的坡度陡变处,两种结果的运动速率都显著增加。但Sentinel-1结果速率分布相比高分三号波动较大,且最大速率低于高分三号结果,其原因应为Sentinel-1空间分辨率较低导致的信息采样不足,一方面降低了结果的信噪比,另一方面也难以反映出局部的运动速率高值。

      将本文监测的2018-01-14—2018-02-12期间的剖面线上冰川平均运动速度与文献[26]采用升降轨COSMO-SkyMed影像计算出的2017年1、2月份的三维平均运动速度进行比较[26],两者具有相似的变化趋势:冰舌末端到高程陡变区,两种结果都呈现出速度先增加后减少至稳定的特征;高程变化最剧烈的区域,均是速度变化最大的区域,也均是冰川最大速度的分布位置;高程陡变区至冰川上部,冰川速率均表现为稳定的速度;本文剖面线上速度最大值为2.08 m/d,文献[26]监测的三维平均运动速度最大值为2.4 m/d,两者最大速度值较为接近,但本文剖面线速度总体小于文献[26]估计的速度。一方面是因为本文获取的是冰川表面的二维速率分量,是文献[26]三维速率在水平面的分量,从而出现剖面线上速率较小的现象,另一方面也可能与剖面线的位置以及获取速率的时间没有完全吻合有关系。相比文献[26]仅涉及了约1个月内的两个时段冰川表面运动状况,未能对冰川表面运动时序特征开展分析,本文获得了冰川2017-11—2018-03共4个时段的冰川表面运动时序结果,对依嘎冰川的认识更为深入。

    • 偏移量跟踪技术的主要误差来源是SAR影像的配准、地形等[27]。基于外部DEM辅助的配准方法可有效降低地形引起的偏移量,提高配准精度。本文使用30 m分辨率的数字表面模型(digital surface model, DSM)辅助SAR影像配准,高分三号影像对配准精度最终优于1/15像元,针对本文使用的高分三号影像,每个影像对的配准误差不超过0.4 m,对应的冰川表面运动速率配准误差约为1.3 cm/d;Sentinel-1影像对配准精度优于1/30像元,配准误差不超过0.5 m,对应冰川表面运动速率配准误差不超过1.4 cm/d。在假设研究区基岩的运动速率应为0的基础上,偏移量跟踪技术的结果精度可以通过统计基岩区的均方根误差(root mean square error, RMSE)来评价[28]。通过在研究区基岩区选取3块区域(图 2(a)中红色矩形ABC区域),分别计算了区域内4个时段速率的均方根误差,计算结果如表 3所示。

      表 3  研究区基岩区速率残差的均方根误差/(cm∙d-1)

      Table 3.  The RMSE of Velocity in the Bedrock Region of Offset Tracking Techniques/(cm∙d-1)

      影像对 A B C 均值
      2017-11-17—2017-12-16 1.1 1.7 1.3 1.4
      2017-12-16—2018-01-14 1.8 2.3 1.8 2.0
      2018-01-14—2018-02-12 1.4 1.5 2.2 1.7
      2018-02-12—2018-03-13 1.8 3.1 2.0 2.3
      2017-12-10—2018-01-15 1.6 1.8 1.3 1.6

      高分三号4个时段及Sentinel-1影像对的基岩区速率残差的平均均方根误差依次为1.4 cm/d、2.0 cm/d、1.7 cm/d、2.3 cm/d、1.6 cm/d, 仅占剖面线上冰川平均运动速率的2.8%、4.1%、3.5%、4.7%、3.2%。其中,分辨率较低的Sentinel-1结果的平均均方根误差与高分三号接近,其原因在于Sentinel-1影像对的配准精度更高,时间基线更长。由以上结果可知,两种数据的均方根误差远小于依嘎冰川的运动速率,验证了本文偏移量跟踪技术结果的可靠性。

    • 本文利用5景高分三号影像,应用偏移量跟踪技术对依嘎冰川2017-11—2018-03期间的表面运动速度进行了监测,得出以下结论:

      1)4个时段的依嘎冰川表面速度变化不明显,速度空间分布具有一致性。从冰川末端到冰川上端,冰川表面速度呈现出先减速后增加到基本稳定,再到显著增加、显著减少至稳定的趋势,其最大速率分布在冰川中部,速率高于1 m/d,与高程联合分析发现,此处也是冰川表面高程变化最大的区域。

      2)与以往研究结果及两景Sentinel-1的偏移量跟踪结果进行对比,并计算高分三号4个时段基岩区速率残差的均方根误差,本文结果与以往学者的研究成果及Sentinel-1的结果十分吻合,且基岩区的均方根误差小于3.1 cm/d,远小于依嘎冰川的运动速率,验证了本文结果的可靠性。

      3)高分三号可以全天时、全天候获取影像,具有高分辨率、大幅宽等优势,能有效应用到冰川表面运动监测中,可作为冰川表面运动定期监测的重要数据来源之一。考虑到其高空间分辨率和高性价比的优势,高分三号必将在中国冰川遥感监测中发挥独特作用。

      本文研究表明国产SAR卫星数据可用于定期监测冰川运动变化情况,但本文研究区尚无实测数据验证,下一步工作将结合地基SAR等实地测量仪器、星地联合验证结果的可靠性。高分三号的业务化数据获取能力与现有成熟星载SAR系统相比仍有不足,随着中国更多星载SAR系统的发射运行,国产SAR数据的规模化应用条件将逐步得到满足。

参考文献 (28)

目录

    /

    返回文章
    返回