留言板

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

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

利用GACOS辅助下InSAR Stacking对金沙江流域进行滑坡监测

张成龙 李振洪 余琛 宋闯 肖儒雅 彭建兵

张成龙, 李振洪, 余琛, 宋闯, 肖儒雅, 彭建兵. 利用GACOS辅助下InSAR Stacking对金沙江流域进行滑坡监测[J]. 武汉大学学报 ● 信息科学版, 2021, 46(11): 1649-1657. doi: 10.13203/j.whugis20200675
引用本文: 张成龙, 李振洪, 余琛, 宋闯, 肖儒雅, 彭建兵. 利用GACOS辅助下InSAR Stacking对金沙江流域进行滑坡监测[J]. 武汉大学学报 ● 信息科学版, 2021, 46(11): 1649-1657. doi: 10.13203/j.whugis20200675
ZHANG Chenglong, LI Zhenhong, YU Chen, SONG Chuang, XIAO Ruya, PENG Jianbing. Landslide Detection of the Jinsha River Region Using GACOS Assisted InSAR Stacking[J]. Geomatics and Information Science of Wuhan University, 2021, 46(11): 1649-1657. doi: 10.13203/j.whugis20200675
Citation: ZHANG Chenglong, LI Zhenhong, YU Chen, SONG Chuang, XIAO Ruya, PENG Jianbing. Landslide Detection of the Jinsha River Region Using GACOS Assisted InSAR Stacking[J]. Geomatics and Information Science of Wuhan University, 2021, 46(11): 1649-1657. doi: 10.13203/j.whugis20200675

利用GACOS辅助下InSAR Stacking对金沙江流域进行滑坡监测

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

国家自然科学基金 41941019

国家自然科学基金 41804005

中央高校基本科研业务费专项资金 300102260301/087

中央高校基本科研业务费专项资金 300102260404/087

详细信息

Landslide Detection of the Jinsha River Region Using GACOS Assisted InSAR Stacking

Funds: 

The National Natural Science Foundation of China 41941019

The National Natural Science Foundation of China 41804005

the Fundamental Research Funds for the Central Universities, CHD 300102260301/087

the Fundamental Research Funds for the Central Universities, CHD 300102260404/087

More Information
  • 摘要: 西藏自治区贡觉县雄松乡至沙东乡金沙江流域作为川藏铁路的必经河流,地形崎岖、地质灾害隐患点多,亟需对该地区隐患点进行全方位的识别。首先,选取61景哨兵一号(Sentinel-1)升轨影像、53景Sentinel-1降轨影像和7景陆地观测技术卫星2号(advanced land observing satellite 2,ALOS-2)升轨影像对研究区域进行滑坡探测与监测。然后,利用合成孔径雷达(interferometric synthetic aperture radar,InSAR)通用型大气改正在线服务(generic atmospheric correction online service for InSAR,GACOS)辅助干涉影像堆叠技术(InSAR Stacking)的方法,获取研究区域雷达视线(line of sight, LOS)方向的InSAR年形变平均速率图,并结合3个轨道的结果提取出沿坡向和垂直滑坡向的平均速率图。最后,与LiCSBAS时间序列分析包的结果进行比较,发现两者具有高度一致性,LOS向年形变速率图像的相关系数在0.92以上,沿坡向和沿垂直滑坡向年形变速率的相关系数在0.85以上,证明了GACOS辅助下InSAR Stacking结果的可靠性。此外,还发现研究区域内沿坡向最大年形变速率为-163 mm/a;结合InSAR形变结果与光学遥感影像解译,可将该滑坡群分为A~G 7个区域进行实时监测。
  • 图  1  研究区域与SAR卫星影像覆盖范围

    Figure  1.  Study Region with the Coverages of SAR Imagery

    图  2  滑动面坐标系与站心坐标系的关系

    Figure  2.  Relationship Between the Sliding-Surface-Defined Coordinate and Topocentric Coordinate Systems

    图  3  InSAR LOS向年形变速率图

    Figure  3.  Annual Deformation Rates in the LOS Direction Revealed by InSAR

    图  4  沿坡向OA和垂直滑坡向OB的年形变速率

    Figure  4.  Annual Deformation Rate in theOA and OB Directions

    图  5  GACOS辅助下InSAR Stacking和LiCSBAS技术LOS向年形变速率相关性分析

    Figure  5.  Correlation Analysis Between GACOS Assisted InSAR Stacking and LiCSBAS Derived LOS Annual Deformation Rates

    图  6  研究区域剖面线abc和滑坡区域A~G

    Figure  6.  Profiles a, b, c in the Study Area and Landslide Areas A-G

    图  7  沿剖面线abc的GACOS辅助下InSAR Stacking和LiCSBAS技术获得的LOS向年形变速率对比

    Figure  7.  Comparison Between GACOS Assisted InSAR Stacking and LiCSBAS Derived LOS Annual Deformation Rates Along Profiles a, b, c

    图  8  滑坡区域A~G的光学遥感影像

    Figure  8.  Landslide Areas A-G Revealed by Optical Imagery

    表  1  ALOS-2升轨影像干涉对

    Table  1.   ALOS-2 Ascending Interferograms

    序号 主影像 从影像 垂直基线/m 时间基线/d
    1 20170724 2017112 50 126
    2 20170724 2018052 -226 308
    3 20170724 2018072 -128 364
    4 20170724 2018111 50 476
    5 20170724 2019052 -109 672
    6 20170724 2019072 -102 728
    7 20171127 2018052 -277 182
    8 20171127 2018072 -178 238
    9 20171127 2018111 1 350
    10 20171127 2019052 -159 546
    11 20171127 2019072 -153 602
    12 20180528 2018072 98 56
    13 20180528 2018111 277 168
    14 20180528 2019052 117 364
    15 20180528 2019072 124 420
    16 20180723 2018111 178 112
    17 20180723 2019052 19 308
    18 20180723 2019072 25 364
    19 20181112 2019052 -159 196
    20 20181112 2019072 -153 252
    21 20190527 2019072 6 56
    下载: 导出CSV
  • [1] 许强, 李为乐, 董秀军, 等. 四川茂县叠溪镇新磨村滑坡特征与成因机制初步研究[J]. 岩石力学与工程学报, 2017, 36(11): 2 612-2 628 https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201711002.htm

    Xu Qiang, Li Weile, Dong Xiujun, et al. The Xinmo Village Landslide on June 24, 2017 in Maoxian County, Sichuan Province: Characteristics and Failure Mechanism[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(11): 2 612-2 628 https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201711002.htm
    [2] Ouyang C J, Zhao W, Xu Q, et al. Failure Mechanisms and Characteristics of the 2016 Catastrophic Rockslide at Su Village, Lishui, China[J]. Landslides, 2018, 15(7): 191-100 http://www.researchgate.net/profile/Dalei_Peng/publication/324588780_Failure_mechanisms_and_characteristics_of_the_2016_catastrophic_rockslide_at_Su_village_Lishui_China/links/5adff9afa6fdcc29358fe239/Failure-mechanisms-and-characteristics-of-the-2016-catastrophic-rockslide-at-Su-village-Lishui-China.pdf
    [3] 许强, 郑光, 李为乐, 等. 2018年10月和11月金沙江白格两次滑坡-堰塞堵江事件分析研究[J]. 工程地质学报, 2018, 26(6): 1 534- 1 551 https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201806016.htm

    Xu Qiang, Zheng Guang, Li Weile, et al. Study on Successive Landslide Damming Events of Jinsha River in Baige Village on Octorber 11 and November 3, 2018 [J]. Journal of Engineering Geology, 2018, 26 (6): 1 534-1 551 https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201806016.htm
    [4] 许强, 董秀军, 李为乐. 基于天-空-地一体化的重大地质灾害隐患早期识别与监测预警[J]. 武汉大学学报·信息科学版, 2019, 44(7): 957-966 doi:  10.13203/j.whugis20190088

    Xu Qiang, Dong Xiujun, Li Weile. Early Identification, Monitoring and Early Warning of Major Geolhazards Based on Space-Air-Ground[J]. Geomatics and Information Science of Wuhan University, 2019, 44(7): 957-966 doi:  10.13203/j.whugis20190088
    [5] 许强. 对地质灾害隐患早期识别相关问题的认识与思考[J]. 武汉大学学报·信息科学版, 2020, 45(11): 1 651-1 659 doi:  10.13203/j.whugis20200043

    Xu Qiang. Understanding and Consideration of Related Issues in Early Identification of Potential Geohazards[J]. Geomatics and Information Science of Wuhan University, 2020, 45(11): 1 651-1 659 doi:  10.13203/j.whugis20200043
    [6] Dai Keren, Li Zhenhong, Xu Qiang, et al. Entering the Era of Earth Observation-Based Landslide Warning Systems: A Novel and Exciting Framework[J]. IEEE Geoscience and Remote Sensing Magazine, 2020, 8(1): 136-153 doi:  10.1109/MGRS.2019.2954395
    [7] Ferretti A, Prati C, Rocca F. Permanent Scatterers in SAR Interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(1): 8-20 doi:  10.1109/36.898661
    [8] Canuti P, Casagli N, Ermini L, et al. Landslide Activity as a Geoindicator in Italy: Significance and New Perspectives from Remote Sensing[J]. Environmental Geology, 2004, 45(7): 907-919 doi:  10.1007/s00254-003-0952-5
    [9] Metternicht G, Hurni L, Gogu R. Remote Sensing of Landslides: An Analysis of the Potential Contribution to Geo-Spatial Systems for Hazard Assessment in Mountainous Environments[J]. Remote Sensing of Environment, 2005, 98(2): 284-303 http://www.cabdirect.org/abstracts/20053194375.html
    [10] Van Westen C J, Castellanos E, Kuriakose S L. Spatial Data for Landslide Susceptibility, Hazard, and Vulnerability Assessment: An Overview[J]. Engineering Geology, 2008, 102(3): 112-131 http://www.sciencedirect.com/science/article/pii/s0013795208001786
    [11] Casagli N, Catani F, Del Ventisette C, et al. Monitoring, Prediction, and Early Warning Using Ground-Based Radar Interferometry[J]. Landslides, 2010, 7(3): 291-301 doi:  10.1007/s10346-010-0215-y
    [12] Martha T R, Kerle N, Jetten V, et al. Characterising Spectral, Spatial and Morphometric Properties of Landslides for Semi-Automatic Detection Using Object-Oriented Methods[J]. Geomorphology, 2010, 116(1): 24-36 http://www.researchgate.net/profile/CJ_Westen/publication/222579232_Characterising_spectral_spatial_and_morphometric_properties_of_landslides_for_semi-automatic_detection_using_object-oriented_methods/links/02bfe50d1e7b641a49000000
    [13] Guzzetti F, Mondini A C, Cardinali M, et al. Landslide Inventory Maps: New Tools for an Old Problem[J]. Earth - Science Reviews, 2012, 112(1): 42-66 http://core.ac.uk/download/pdf/82520789.pdf
    [14] Ye X, Kaufmann H, Guo X F. Landslide Monitoring in the Three Gorges Area Using D-InSAR and Corner Reflectors[J]. Photogrammetric Engineering & Remote Sensing, 2004, 70(10): 1 167-1 172 http://essential.metapress.com/content/k3v513t523111536/fulltext.pdf?page=1
    [15] 李振洪, 宋闯, 余琛, 等. 卫星雷达遥感在滑坡灾害探测和监测中的应用: 挑战与对策[J]. 武汉大学学报·信息科学版, 2019, 44(7): 967-979 doi:  10.13203/j.whugis20190098

    Li Zhenhong, Song Chuang, Yu Chen, et al. Application of Satellite Radar Remote Sensing in Landslide Detection and Monitoring: Challenges and Solutions[J]. Geomatics and Information Science of Wuhan University, 2019, 44(7): 967-979 doi:  10.13203/j.whugis20190098
    [16] Tofani V, Segoni S, Agostini A, et al. Technical Note: Use of Remote Sensing for Landslide Studies in Europe[J]. Natural Hazards and Earth System Sciences, 2013, 13(2): 299-309 doi:  10.5194/nhess-13-299-2013
    [17] Riddick S N, Schmidt D A, Deligne N I. An Analysis of Terrain Properties and the Location of Surface Scatterers from Persistent Scatterer Interferometry[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2012, 73: 50-57 doi:  10.1016/j.isprsjprs.2012.05.010
    [18] LiZ H, Fielding E J, Cross P, et al. Interferometric Synthetic Aperture Radar Atmospheric Correction: GPS Topography-Dependent Turbulence Mode l[J]. Journal of Geophysical Research: Solid Earth, 2006, 111(B2): B02404 doi:  10.1029/2005JB003711/pdf
    [19] Gong W, Meyer F J, Liu S Z, et al. Temporal Filtering of InSAR Data Using Statistical Parameters from NWP Models[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(7): 4 033- 4 044 doi:  10.1109/TGRS.2015.2389143
    [20] Torres R, Snoeij P, Geudtner D, et al. GMES Sentinel-1 Mission[J]. Remote Sensing of Environment, 2012, 120: 9-24 doi:  10.1016/j.rse.2011.05.028
    [21] Werner C, Wegmüller U, Strozzi T, et al. Gamma SAR and Interferometric Processing Software[C]. The 2000 ERS - Envisat Symposium, Gothenburg, Sweden, 2000
    [22] Yagüe-Martínez N, Prats-Iraola P, Rodríguez González F, et al. Interferometric Processing of Sentinel-1 TOPS Data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(4): 2 220- 2 234 doi:  10.1109/TGRS.2015.2497902
    [23] Yu C, Penna N T, Li Z H. Generation of RealTime Mode High-Resolution Water Vapor Fields from GPS Observations[J]. Journal of Geophysical Research: Atmospheres, 2017, 122(3): 2 008-2 025 doi:  10.1002/2016JD025753
    [24] Wright T, Parsons B, Fielding E. Measurement of Interseismic Strain Accumulation across the North Anatolian Fault by Satellite Radar Interferometry[J]. Geophysical Research Letters, 2001, 28(10): 2 117-2 120 doi:  10.1029/2000GL012850
    [25] Xiao R Y, Yu C, Li Z H, et al. General Survey of Large-Scale Land Subsidence by GACOS-Corrected InSAR Stacking: Case Study in North China Plain[J]. Proceedings of the International Association of Hydrological Sciences, 2020, 382: 213-218 doi:  10.5194/piahs-382-213-2020
    [26] Morishita Y, Lazecky M, Wright T, et al. LiCSBAS: An Open-Source InSAR Time Series Analysis Package Integrated with the LiCSAR Automated Sentinel-1 InSAR Processor[J]. Remote Sensing, 2020, 12(3): 424 doi:  10.3390/rs12030424
    [27] 李永生, 张景发, 李振洪, 等. 时序InSAR离散相干点相位解缠误差检查与校正方法研究[J]. 武汉大学学报· 信息科学版, 2014, 39(10): 1 199-1 203 http://ch.whu.edu.cn/article/id/3095

    Li Yongsheng, Zhang Jingfa, Li Zhenhong, et al. Detection and Correction of Phase Unwrapping Errors in InSAR Time Series Analysis with Discrete Coherent Points[J]. Geomatics and Information Science of Wuhan University, 2014, 39(10): 1 199-1 203 http://ch.whu.edu.cn/article/id/3095
    [28] Hu X, Lu Z, Pierson T C, et al. Combining InSAR and GPS to Determine Transient Movement and Thickness of a Seasonally Active Low-Gradient Translational Landslide[J]. Geophysical Research Letters, 2018, 45(3): 1 453-1 462 doi:  10.1002/2017GL076623
  • [1] 厉飞, 闫庆武, 邹雅婧, 刘保丽.  利用夜间灯光POI的城市建成区提取精度研究——以珞珈一号01星和NPP/VIIRS夜间灯光影像为例 . 武汉大学学报 ● 信息科学版, 2021, 46(6): 825-835. doi: 10.13203/j.whugis20190266
    [2] 刘权毅, 詹庆明, 李建松, 杨晨, 刘稳.  珞珈一号夜间灯光影像在建设用地提取中的应用:以武汉市为例 . 武汉大学学报 ● 信息科学版, 2021, 46(1): 30-39. doi: 10.13203/j.whugis20190376
    [3] 柳林, 宋豪峰, 杜亚男, 冯光财, 刘清瑶, 孙敏.  联合哨兵2号和Landsat 8估计白格滑坡时序偏移量 . 武汉大学学报 ● 信息科学版, 2021, 46(10): 1461-1470. doi: 10.13203/j.whugis20200596
    [4] 高梦瑶, 许才军, 刘洋.  青藏高原西北缘时序InSAR对流层延迟改正方法评估 . 武汉大学学报 ● 信息科学版, 2021, 46(10): 1548-1559. doi: 10.13203/j.whugis20210254
    [5] 沈维青, 张丰, 刘仁义, 杜震洪.  利用SARA结合高分一号数据反演气溶胶光学厚度 . 武汉大学学报 ● 信息科学版, 2020, 45(6): 914-922. doi: 10.13203/j.whugis20190034
    [6] 陆会燕, 李为乐, 许强, 董秀军, 代聪, 王栋.  光学遥感与InSAR结合的金沙江白格滑坡上下游滑坡隐患早期识别 . 武汉大学学报 ● 信息科学版, 2019, 44(9): 1342-1354. doi: 10.13203/j.whugis20190086
    [7] 唐尧, 王立娟, 马国超, 贾虎军, 李洪梁, 尹恒.  基于“高分+”的金沙江滑坡灾情监测与应用前景分析 . 武汉大学学报 ● 信息科学版, 2019, 44(7): 1082-1092. doi: 10.13203/j.whugis20190048
    [8] 李振洪, 宋闯, 余琛, 肖儒雅, 陈立福, 罗慧, 戴可人, 葛大庆, 丁一, 张宇星, 张勤.  卫星雷达遥感在滑坡灾害探测和监测中的应用:挑战与对策 . 武汉大学学报 ● 信息科学版, 2019, 44(7): 967-979. doi: 10.13203/j.whugis20190098
    [9] 王磊, 陈锐志, 李德仁, 蔚保国, 伍蔡伦.  珞珈一号低轨卫星导航增强系统信号质量评估 . 武汉大学学报 ● 信息科学版, 2018, 43(12): 2191-2196. doi: 10.13203/j.whugis20180413
    [10] 徐文, 邵俊, 喻文勇, 冯鹏.  陆地观测卫星数据中心:大数据挑战及一种解决方案 . 武汉大学学报 ● 信息科学版, 2017, 42(1): 7-13. doi: 10.13203/j.whugis20150172
    [11] 张艳, 王涛, 冯伍法, 王淑香.  “天绘一号”海岛(礁)影像稀少控制下的定位技术研究 . 武汉大学学报 ● 信息科学版, 2016, 41(5): 617-623. doi: 10.13203/j.whugis20140220
    [12] 康一飞, 王树根, 韩飞飞, 孙明伟.  资源一号02C影像条带噪声去除的改进矩匹配方法 . 武汉大学学报 ● 信息科学版, 2015, 40(12): 1582-1587. doi: 10.13203/j.whugis20150717
    [13] 陈王丽, 孙 涛, 秦前清, 石 岩, 张 巍, 林立宇.  利用嫦娥一号ccd相机邻轨下视影像生成月面dem . 武汉大学学报 ● 信息科学版, 2015, 40(9): 1169-1175. doi: 10.13203/j .whu g is20140041
    [14] 陈王丽, 孙涛, 陈喆3马国锐, 秦前清.  利用光流配准进行嫦娥一号CCD多视影像超分辨率重建 . 武汉大学学报 ● 信息科学版, 2014, 39(9): 1103-1108. doi: 10.13203/j.whugis20130289
    [15] 曹建峰, 胡松杰, 黄勇, 刘磊.  嫦娥二号卫星日地拉格朗日L_2点探测轨道定轨分析 . 武汉大学学报 ● 信息科学版, 2013, 38(9): 1029-1033.
    [16] 严韦, 魏二虎, 刘经南.  ΔVLBI用于“嫦娥一号”地月转移轨道段定轨及EOP解算 . 武汉大学学报 ● 信息科学版, 2012, 37(8): 960-962.
    [17] 何雨帆, 王家松, 李远平, 陈建荣.  基于北斗一号的近地卫星天基测控技术及应用 . 武汉大学学报 ● 信息科学版, 2012, 37(4): 441-444.
    [18] 魏二虎, 史青, 严韦, 李雪川.  ERP精度对“嫦娥一号”差分VLBI定位精度的影响 . 武汉大学学报 ● 信息科学版, 2011, 36(11): 1324-1327.
    [19] 陈明, 唐歌实, 曹建峰, 张宇.  嫦娥一号绕月探测卫星精密定轨实现 . 武汉大学学报 ● 信息科学版, 2011, 36(2): 212-217.
    [20] 刘山洪, 鄢建国, 杨轩, 叶茂, 郭茜, 王波, 李斐.  天问一号拓展任务对火星低阶重力场解算的潜在贡献分析 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20210035
  • 加载中
图(8) / 表(1)
计量
  • 文章访问数:  77
  • HTML全文浏览量:  29
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-12-13
  • 刊出日期:  2021-11-05

利用GACOS辅助下InSAR Stacking对金沙江流域进行滑坡监测

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

    国家自然科学基金 41941019

    国家自然科学基金 41804005

    中央高校基本科研业务费专项资金 300102260301/087

    中央高校基本科研业务费专项资金 300102260404/087

    作者简介:

    张成龙,博士生,主要从事InSAR地质灾害早期识别研究。chenglongzhang136@163.com

    通讯作者: 李振洪,博士,教授。zhenhong.li@chd.edu.cn
  • 中图分类号: P237

摘要: 西藏自治区贡觉县雄松乡至沙东乡金沙江流域作为川藏铁路的必经河流,地形崎岖、地质灾害隐患点多,亟需对该地区隐患点进行全方位的识别。首先,选取61景哨兵一号(Sentinel-1)升轨影像、53景Sentinel-1降轨影像和7景陆地观测技术卫星2号(advanced land observing satellite 2,ALOS-2)升轨影像对研究区域进行滑坡探测与监测。然后,利用合成孔径雷达(interferometric synthetic aperture radar,InSAR)通用型大气改正在线服务(generic atmospheric correction online service for InSAR,GACOS)辅助干涉影像堆叠技术(InSAR Stacking)的方法,获取研究区域雷达视线(line of sight, LOS)方向的InSAR年形变平均速率图,并结合3个轨道的结果提取出沿坡向和垂直滑坡向的平均速率图。最后,与LiCSBAS时间序列分析包的结果进行比较,发现两者具有高度一致性,LOS向年形变速率图像的相关系数在0.92以上,沿坡向和沿垂直滑坡向年形变速率的相关系数在0.85以上,证明了GACOS辅助下InSAR Stacking结果的可靠性。此外,还发现研究区域内沿坡向最大年形变速率为-163 mm/a;结合InSAR形变结果与光学遥感影像解译,可将该滑坡群分为A~G 7个区域进行实时监测。

English Abstract

张成龙, 李振洪, 余琛, 宋闯, 肖儒雅, 彭建兵. 利用GACOS辅助下InSAR Stacking对金沙江流域进行滑坡监测[J]. 武汉大学学报 ● 信息科学版, 2021, 46(11): 1649-1657. doi: 10.13203/j.whugis20200675
引用本文: 张成龙, 李振洪, 余琛, 宋闯, 肖儒雅, 彭建兵. 利用GACOS辅助下InSAR Stacking对金沙江流域进行滑坡监测[J]. 武汉大学学报 ● 信息科学版, 2021, 46(11): 1649-1657. doi: 10.13203/j.whugis20200675
ZHANG Chenglong, LI Zhenhong, YU Chen, SONG Chuang, XIAO Ruya, PENG Jianbing. Landslide Detection of the Jinsha River Region Using GACOS Assisted InSAR Stacking[J]. Geomatics and Information Science of Wuhan University, 2021, 46(11): 1649-1657. doi: 10.13203/j.whugis20200675
Citation: ZHANG Chenglong, LI Zhenhong, YU Chen, SONG Chuang, XIAO Ruya, PENG Jianbing. Landslide Detection of the Jinsha River Region Using GACOS Assisted InSAR Stacking[J]. Geomatics and Information Science of Wuhan University, 2021, 46(11): 1649-1657. doi: 10.13203/j.whugis20200675
  • 滑坡作为最为常见的地质灾害之一,对经济、社会、人员安全等方面带来巨大损失[1-2]。川藏铁路沿线及邻区地形地貌、地质构造、气候水文等各方面极为复杂,尤其是在大江大河流域,多为高原深切峡谷,峡谷之间的相对高差大多为2 000~3 000 m,最大可达5 000 m。金沙江流域地质灾害频发,2018-10-11和2018-11-03在西藏江达县波罗乡白格村附近金沙江流域发生两次大规模高位滑坡,对周围居民以及河道造成严重损害[3],因此对金沙江流域进行滑坡隐患定量全方位的探测,具有重要的社会和工程应用价值。

    近年来,越来越多的学者致力于滑坡隐患点的探测和监测研究,传统的野外调查无法探测一些高位和隐蔽性的滑坡,随着高精度对地观测技术的发展,尤其是光学遥感、雷达遥感、机载激光雷达(light detection and ranging,LiDAR)等技术的出现与普及,基于多平台、多技术、多源数据的集成技术在滑坡隐患点的探测和监测方面的研究日益增多。文献[4-5]提出利用天-空-地一体化的三查体系提高中国的地质灾害识别能力。文献[6]提出基于对地观测技术的地质滑坡预警系统,利用星载对地观测技术进行广域的滑坡隐患点探测,然后通过数值模拟和野外调查方法筛选出高风险滑坡体,在重点区域安装北斗/全球导航卫星系统(global navigation satellite system,GNSS)、裂缝计、倾斜计、雨量计等传感器实时监测滑坡体,进而作早期预警,并通过与政府和社区的合作确保预警系统发挥作用。文献[7-14]对滑坡区域勘察、滑坡的探测和监测、滑坡的地质灾害评估以及滑坡的危险性评估进行了研究。随着对地观测技术的迅猛发展,其在滑坡识别、测绘、监测和灾害评估等方面的应用也日益增加。

    合成孔径雷达干涉测量(interferometric synthetic aperture radar,InSAR)作为一项极具潜力的空间对地观测技术,在大范围形变监测方面具有突出优势[15-16]。InSAR在滑坡探测与监测应用中也存在以下问题:(1)相干点密度低,空间分布严重不均匀。高植被、高海拔落差、高冰雪覆盖等因素都会导致InSAR相干点密度低、空间分布不均匀,从而影响到形变量的解算精度[17]。(2)大气延迟误差的干扰。多变的气象和巨大的地形起伏可导致明显的大气延迟误差,直接限制InSAR提取的形变量精度[18-19]。(3)高海拔落差研究区地质灾害具有多样性和复杂性的特点,往往需要从三维方向观测其形变,而传统InSAR获取的形变只是在一维视线向上的投影,难以反映真实的灾害变形情况。

    位于川藏铁路沿线的金沙江流域具有高海拔落差、相干点密度低、空间分布严重不均匀等特点,对雷达遥感影像处理和分析提出了很高的要求。本文以金沙江流域为研究区域,利用合成孔径雷达通用型大气改正在线服务(generic atmospheric correction online service for InSAR,GACOS)辅助下干涉影像堆叠技术(InSAR Stacking)和LiCSBAS时间序列分析包获取雷达视线方向(line of sight,LOS)年形变速率,进而提取出二维的形变场,并针对研究区域开展了滑坡探测、监测和对比分析。

    • 本文研究区域位于西藏自治区贡觉县雄松乡(98.92°E,30.49°N)至沙东乡(98.9°E,30.61°N)的金沙江流域(见图 1),从谷歌影像目视解译发现可能存在滑坡。本文选取南北约15 km、东西约6 km的矩形区域进行滑坡探测实验,黄色矩形为研究区域,蓝色矩形为哨兵一号(Sentinel-1)升轨影像覆盖范围,绿色矩形为Sentinel-1降轨影像覆盖范围,红色矩形为陆地观测技术卫星2号(advanced land observing satellite 2,ALOS-2)升轨影像覆盖范围。研究区域海拔最高为4 608 m,最低为2 594 m,地形落差大,地形复杂,在河流附近极易孕育出滑坡群。

      图  1  研究区域与SAR卫星影像覆盖范围

      Figure 1.  Study Region with the Coverages of SAR Imagery

      Sentine-l系列卫星是欧空局(European Space Agency,ESA)由哥白尼方案资助的对地观测卫星。Sentinel-1A卫星和Sentinel-1B卫星运行在同一轨道平面上,两颗卫星星座昼夜工作,并通过哥白尼访问中心提供免费、开放、连续的产品。该卫星工作在C波段,有条带模式(strip mapmode,SM)、干涉宽幅模式(interferometric wideswath,IW)、超宽幅模式(extra wideswath,EW)、波模式(wave mode,WV)4种工作模式。本文选择了IW模式下垂直发射垂直接收的同向极化方式(vertical vertical,VV)得到的单视复数影像(single look complex,SLC)。无论天气如何,Sentinel-1单颗卫星采用渐进式扫描地形观测模式(terrain observation by progressive scans,TOPS),每12 d就可以获得全球范围内的图像[20]。本文选取2017-10-08—2019-10-10共61景Sentinel-1升轨影像和2017-07-23—2019-07-25共53景Sentinel-1降轨影像用于分析。

      ALOS-2卫星携带L波段合成孔径雷达(synthetic aperture radar,SAR)传感器,穿透力强,具有聚光模式(Spotlight)、SM模式、扫描模式(ScanSAR)3种观测模式。由于观察模式和研究区域的不同,重访时间也从几个月到两周不等,与C波段相比,L波段由于雷达波长较长,能够更好地穿透植被获取地面信息,便于在高海拔、高植被覆盖地区保持较好的相干性,因此,本文选取SM3模式下水平发射水平接收的同向极化方式(horizontal horizontal,HH)得到的SLC,2017-07-24—2019-07-22 7景ALOS-2升轨影像用于研究。

    • 本文采用Gamma处理1级单视复数格式的SAR影像[21],为了抑制斑点噪声,Sentinel-1影像在距离向和方位向采用10×2的多视因子。ALOS-2影像在距离向和方位向采用2×5的多视因子,与其他传统的SAR影像数据处理不同,Sentinel-1在TOPS模式下获取数据时,由于天线会发生旋转,需要更严格的图像配准,配准的方位角精度应优于0.001像素,以避免相邻Burst之间的相位跳跃[22]

      利用30 m分辨率的航天飞机雷达地形测绘使命(shuttle radar topography mission,SRTM)数字高程模型(digital elevation model,DEM)(https://earthexplorer.usgs.gov)削弱地形的影响,Sentinel-1影像设置的时间基线为90 d,空间基线为150 m,由于ALOS-2数据偏少,因此本文对ALOS-2数据进行全干涉,得到ALOS-2数据干涉对组合21对(见表 1),Sentinel-1升轨影像干涉对组合212对,Sentinel-1降轨影像干涉对组合188对。

      表 1  ALOS-2升轨影像干涉对

      Table 1.  ALOS-2 Ascending Interferograms

      序号 主影像 从影像 垂直基线/m 时间基线/d
      1 20170724 2017112 50 126
      2 20170724 2018052 -226 308
      3 20170724 2018072 -128 364
      4 20170724 2018111 50 476
      5 20170724 2019052 -109 672
      6 20170724 2019072 -102 728
      7 20171127 2018052 -277 182
      8 20171127 2018072 -178 238
      9 20171127 2018111 1 350
      10 20171127 2019052 -159 546
      11 20171127 2019072 -153 602
      12 20180528 2018072 98 56
      13 20180528 2018111 277 168
      14 20180528 2019052 117 364
      15 20180528 2019072 124 420
      16 20180723 2018111 178 112
      17 20180723 2019052 19 308
      18 20180723 2019072 25 364
      19 20181112 2019052 -159 196
      20 20181112 2019072 -153 252
      21 20190527 2019072 6 56
    • GACOS系统基于高分辨率欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts,ECMWF)0.1°×0.1°格网的大气数值产品,采用迭代对流层分解(iterative tropospheric decomposition,ITD)模型从对流层延迟中分离出分层和湍流分量,分别进行空间内插,得到高空间分辨率(默认为90 m)的天顶对流层延迟(zenith troposphere delay,ZTD)改正影像[23]

      InSAR Stacking的基本假设是研究区域内的地表形变是线性的,不存在非线性的变化,其优点是运算高效,可有效削弱随机误差的影响。引入GACOS的ZTD产品削弱了大气误差,残余的大气误差可视为随机误差,通过InSAR Stacking技术可提高信噪比[24]。本文利用最小二乘对解缠相位求平均,干涉相位按时间间隔加权,得到研究区域LOS向年形变速率为[25]

      Vmean=i=1NφiΔTi=1NΔT2 ]]>

      ,Vmean为平均相位速率;φi为解缠后的相位;ΔT为对于干涉图的时间基线;N为干涉图的个数。

    • LiCSBAS是一个与LiCSAR系统集成的开放源InSAR时间序列分析软件包[16],本文利用LiCSBAS解算InSAR雷达LOS位移时间序列。尽管Sentinel‐1影像的空间和时间基线较短,但由于研究区域的高植被、高积雪覆盖以及地形起伏出现严重的失相关和解缠误差,在进行时序分析之前需要对解缠后的干涉图进行质量检查,修改或删除[26]。解算步骤如下:首先利用GACOS对流层延迟产品对干涉对进行有效的改正,剔除解缠误差,对整体干涉对进行质检。以SAR影像的有效像元率(每对干涉图上有效像元数/有效像元总数)和平均相干性作为阈值,对不满足条件的干涉对进行剔除,Sentinel-1升轨影像干涉对中22个干涉对被剔除,Sentinel-1降轨影像干涉对中10个干涉对被剔除。然后对干涉图进行相位环闭合的检验[27],逐一计算干涉对环闭合相位均方根(root mean square,RMS),为了防止RMS满足阈值的干涉对中仍然存在解缠误差的像元,该方法统计了每个像元环闭合相位误差大于阈值(1.5 rad)的个数,作为掩模像元的标准。最后在满足条件的情况下,得到研究区域的LOS向的年形变速率。

    • 坡向是指坡面法线在水平面上投影的方向,坡度代表地表单元陡缓的程度,通常把坡面垂直高度与水平方向距离的比叫做坡度。坡度和坡向大致能反映一个地形的特征要素,往往滑坡的发生除了必要的外部环境条件,其位移特征与坡度和坡向有着很大的关系。

      本文基于滑坡的运动规律,以滑坡的滑动面作为观测其运动的参考面,引入坡度和坡向两个地形因子作为滑动面坐标系的要素。在现实生活中滑动面是不断变化的,本文假设其为一个均匀的滑动面,设滑坡滑动面的坡向轴为OA,沿坡向OA向下为正,垂直于滑动面的法线方向OC向上为正,垂直于面AOCOB轴向外为正构建以滑动面为参考面的独立坐标系。因为滑动面坐标系与站心坐标系(east/north/up,E/N/U)存在几何关系,见图 2。因此可将站心坐标系作为过渡,由LOS的年形变速率得到OAOBOC方向的年形变速率。在自然环境下,除了借助外力,一般沿滑动面的法线方向OC形变远远小于其他方向,所以本文将OC方向的结果视为零[28]

      图  2  滑动面坐标系与站心坐标系的关系

      Figure 2.  Relationship Between the Sliding-Surface-Defined Coordinate and Topocentric Coordinate Systems

      由此得到三轨雷达影像的LOS向年形变速率,沿滑坡向(OA)、垂直滑坡向(OB)的年形变速率计算如下:

      VLOS=[sinθcosα-sinθcosαcosθ]VNVEVU ]]>
      VNVEVU=cosβcosγsinγsinβcosγcosβsinγ-cosγsinβsinγ-sinβ0cosβVAVBVC ]]>
      VLOS=abcVAVBVCTa=cosθcosβsinγ+sinθsinαcosβcosγ+sinθcosαsinβb=-cosθcosα+sinθsinαsinγc=cosθsinβsinγ+sinθsinαsinβcosγ-sinθcosαcosβ ]]>

      ,VLOS为SAR卫星LOS向的年形变速率;θ为雷达卫星的入射角;α为雷达卫星的方位角;VNVEVU分别为在站心坐标系下N、E、U方向的年形变速率;β为坡度;γ为坡向;VAVBVC分别为滑坡向、垂直滑坡向、滑坡面法向(即垂直方向)的年形变速率。

    • 本文将GACOS辅助下的InSAR Stacking和LiCSBAS技术分别应用于西藏自治区贡觉县雄松乡至沙东乡的金沙江流域进行滑坡群探测与监测,利用2017—2019年Sentinel-1和ALOS-2 3个轨道约2 a的影像对该滑坡群进行探测得到LOS向年形变速率,结果如图 3所示,负值表示远离雷达视线方向,正值表示靠近雷达视线方向。图 4为沿坡向和垂直滑坡向的年形变速率。

      图  3  InSAR LOS向年形变速率图

      Figure 3.  Annual Deformation Rates in the LOS Direction Revealed by InSAR

      图  4  沿坡向OA和垂直滑坡向OB的年形变速率

      Figure 4.  Annual Deformation Rate in theOA and OB Directions

      由图34可以看出,该滑坡群有多处均发生明显形变,两种方法所得结果的量级和滑坡范围基本一致。将两种方法所得结果进行相关性分析,结果如图 5所示。

      图  5  GACOS辅助下InSAR Stacking和LiCSBAS技术LOS向年形变速率相关性分析

      Figure 5.  Correlation Analysis Between GACOS Assisted InSAR Stacking and LiCSBAS Derived LOS Annual Deformation Rates

      图 5可知,两者结果具有一致性,LOS向年形变速率的R2在0.92以上,沿坡向和沿垂直滑坡向年形变速率的R2在0.85以上。ALOS-2影像所得结果的标准差为1.3 mm/a,Sentinel-1A升轨影像所得结果的标准差为5.4 mm/a,Sentinel-1A降轨影像所得结果的标准差为6.8 mm/a,而由图 1可知,沿坡向年形变速率标准差为7.2 mm/a,沿垂直滑坡向年形变速率标准差为7.4 mm/a,这些均证明了GACOS辅助下InSAR Stacking结果的可靠性。

      为了对结果进行更加详细和直观的分析,在形变量较大的滑坡上随机选取3个剖面线abc,并根据滑坡探测所得结果,将该滑坡群分为ABCDEFG 7个区域,见图 6图 7为两种方法沿剖面线abc获得的LOS向年形变速率对比。由图 7可知,两种方法所得结果相关性都在0.9以上,尤其是Sentinel-1升轨和ALOS-2影像的结果相关性都在0.95以上。对于ALOS-2影像,两种方法计算结果的RMS均在2 mm/a以下,Sentinel-1升轨影像的两种方法解算结果的RMS均在2~4 mm/a之间;相较于其他两轨道的对比结果,Sentinel-1降轨影像的两种方法解算结果的RMS较大,剖线a的RMS为8.5 mm/a,剖线b的RMS为4 mm/a,剖线c的RMS为3.7 mm/a,这可能是由于降轨影像采集的时间在白天,而升轨影像采集的时间在晚上,白天大气水汽时空变化可能比较大,扰动型误差占比较大,而GACOS改正后的残余也随之变大。通过以上对比分析更加进一步验证该区域滑坡识别的可靠性。

      图  6  研究区域剖面线abc和滑坡区域A~G

      Figure 6.  Profiles a, b, c in the Study Area and Landslide Areas A-G

      图  7  沿剖面线abc的GACOS辅助下InSAR Stacking和LiCSBAS技术获得的LOS向年形变速率对比

      Figure 7.  Comparison Between GACOS Assisted InSAR Stacking and LiCSBAS Derived LOS Annual Deformation Rates Along Profiles a, b, c

      图 8为滑坡区域A~G的光学遥感影像,通过谷歌地球将滑坡区域标记,目视解译发现7个区域均有滑坡迹象,该滑坡群在相同的监测时间段内表现出不同的区分性。

      图  8  滑坡区域A~G的光学遥感影像

      Figure 8.  Landslide Areas A-G Revealed by Optical Imagery

      图 8可知,区域ABCD为较活跃滑坡,具有整体失稳的可能。区域A滑坡覆盖范围约为8.5km2,通过光学影像可以看出该区域为积层滑坡,在滑坡的前后缘均存在明显的形变,尤其是在滑坡前缘区域沿坡向最大年形变速率为-142 mm/a,可能原因是常年受河水的侵蚀,导致局部发生垮塌,可判定为老滑坡;区域B滑坡覆盖范围约为3km2,沿坡向最大年形变速率为-82 mm/a,通过光学遥感影像可以判断该滑坡为岩质滑坡,该滑坡后缘裂缝明显,前缘上部多级滑动交错,中部滑动明显;区域C和区域D滑坡覆盖范围分别约为12km2、2.4km2,同属于积层滑坡,尤其是区域C前后缘均存在明显的形变,该区域属于迎风坡,常年受强烈的太阳光照,坡脚被江水掏蚀严重,沿坡向最大年形变速率为-163 mm/a。区域F属于岩质滑坡,前后缘均形变缓慢,区域EG为积层滑坡,目前区域EFG滑坡仍处于缓慢变形阶段。区域ABCD形变明显,在强降雨江水的侵蚀以及地震等诱发下,可能会整体失稳造成堵江后果。后续应继续对该区域进行InSAR时间序列分析,在条件允许情况下进行现场设备监测,避免不必要的损失。

    • 本文对西藏自治区贡觉县雄松乡至沙东乡的金沙江流域范围进行滑坡探测与监测,获取该区域的LOS向、坡向和垂直滑坡向年形变速率,沿坡向最大年形变速率为-163 mm/a,探测结果较可靠。并通过谷歌地球光学遥感影像对滑坡群进行判别,发现了一个包含7处滑坡隐患点(A~G)的滑坡群。该滑坡群形成的原因包括:(1)研究区域海拔为4 000 m左右,距离冰川近,近年来气温上升,冰川融水增加,导致以冰川为驱动力的滑坡数量增加;(2)在2018-10-11和2018-11-03,研究区域上游白格村附近连续两次发生大规模高位滑坡形成堰塞湖,在人工泄洪时,对该滑坡群前缘带来较大影响。探测到7个隐患点具有较明显的形变,区域ABCD为较活跃滑坡,可能会整体失稳造成堵江后果,目前区域EFG滑坡仍处于缓慢变形阶段。GACOS辅助下InSAR Stacking技术可有效去除长波段及与地形相关大气延迟误差,具有简单、有效、快速、易于推广应用等优点。

参考文献 (28)

目录

    /

    返回文章
    返回