-
滑坡作为最为常见的地质灾害之一,对经济、社会、人员安全等方面带来巨大损失[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,地形落差大,地形复杂,在河流附近极易孕育出滑坡群。
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)
格网的大气数值产品,采用迭代对流层分解(iterative tropospheric decomposition,ITD)模型从对流层延迟中分离出分层和湍流分量,分别进行空间内插,得到高空间分辨率(默认为90 m)的天顶对流层延迟(zenith troposphere delay,ZTD)改正影像[23]。 InSAR Stacking的基本假设是研究区域内的地表形变是线性的,不存在非线性的变化,其优点是运算高效,可有效削弱随机误差的影响。引入GACOS的ZTD产品削弱了大气误差,残余的大气误差可视为随机误差,通过InSAR Stacking技术可提高信噪比[24]。本文利用最小二乘对解缠相位求平均,干涉相位按时间间隔加权,得到研究区域LOS向年形变速率为[25]:
为平均相位速率; 为解缠后的相位; 为对于干涉图的时间基线;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向的年形变速率。
-
坡向是指坡面法线在水平面上投影的方向,坡度代表地表单元陡缓的程度,通常把坡面垂直高度与水平方向距离的比叫做坡度。坡度和坡向大致能反映一个地形的特征要素,往往滑坡的发生除了必要的外部环境条件,其位移特征与坡度和坡向有着很大的关系。
本文基于滑坡的运动规律,以滑坡的滑动面作为观测其运动的参考面,引入坡度和坡向两个地形因子作为滑动面坐标系的要素。在现实生活中滑动面是不断变化的,本文假设其为一个均匀的滑动面,设滑坡滑动面的坡向轴为
,沿坡向 向下为正,垂直于滑动面的法线方向 向上为正,垂直于面AOC的 轴向外为正构建以滑动面为参考面的独立坐标系。因为滑动面坐标系与站心坐标系(east/north/up,E/N/U)存在几何关系,见图 2。因此可将站心坐标系作为过渡,由LOS的年形变速率得到 、 、 方向的年形变速率。在自然环境下,除了借助外力,一般沿滑动面的法线方向 形变远远小于其他方向,所以本文将 方向的结果视为零[28]。 图 2 滑动面坐标系与站心坐标系的关系
Figure 2. Relationship Between the Sliding-Surface-Defined Coordinate and Topocentric Coordinate Systems
由此得到三轨雷达影像的LOS向年形变速率,沿滑坡向(
)、垂直滑坡向( )的年形变速率计算如下: 为SAR卫星LOS向的年形变速率; 为雷达卫星的入射角; 为雷达卫星的方位角; 、 、 分别为在站心坐标系下N、E、U方向的年形变速率; 为坡度; 为坡向; 、 、 分别为滑坡向、垂直滑坡向、滑坡面法向(即垂直方向)的年形变速率。 -
本文将GACOS辅助下的InSAR Stacking和LiCSBAS技术分别应用于西藏自治区贡觉县雄松乡至沙东乡的金沙江流域进行滑坡群探测与监测,利用2017—2019年Sentinel-1和ALOS-2 3个轨道约2 a的影像对该滑坡群进行探测得到LOS向年形变速率,结果如图 3所示,负值表示远离雷达视线方向,正值表示靠近雷达视线方向。图 4为沿坡向和垂直滑坡向的年形变速率。
由图3、4可以看出,该滑坡群有多处均发生明显形变,两种方法所得结果的量级和滑坡范围基本一致。将两种方法所得结果进行相关性分析,结果如图 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向年形变速率的
在0.92以上,沿坡向和沿垂直滑坡向年形变速率的 在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个剖面线a、b、c,并根据滑坡探测所得结果,将该滑坡群分为
、 、 、 、 、 、 7个区域,见图 6。图 7为两种方法沿剖面线a、b、c获得的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改正后的残余也随之变大。通过以上对比分析更加进一步验证该区域滑坡识别的可靠性。 图 7 沿剖面线a、b、c的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、B、C、D为较活跃滑坡,具有整体失稳的可能。区域A滑坡覆盖范围约为8.5
,通过光学影像可以看出该区域为积层滑坡,在滑坡的前后缘均存在明显的形变,尤其是在滑坡前缘区域沿坡向最大年形变速率为 142 mm/a,可能原因是常年受河水的侵蚀,导致局部发生垮塌,可判定为老滑坡;区域B滑坡覆盖范围约为3 ,沿坡向最大年形变速率为 82 mm/a,通过光学遥感影像可以判断该滑坡为岩质滑坡,该滑坡后缘裂缝明显,前缘上部多级滑动交错,中部滑动明显;区域C和区域D滑坡覆盖范围分别约为12 、2.4 ,同属于积层滑坡,尤其是区域C前后缘均存在明显的形变,该区域属于迎风坡,常年受强烈的太阳光照,坡脚被江水掏蚀严重,沿坡向最大年形变速率为 163 mm/a。区域F属于岩质滑坡,前后缘均形变缓慢,区域E、G为积层滑坡,目前区域E、F、G滑坡仍处于缓慢变形阶段。区域A、B、C、D形变明显,在强降雨江水的侵蚀以及地震等诱发下,可能会整体失稳造成堵江后果。后续应继续对该区域进行InSAR时间序列分析,在条件允许情况下进行现场设备监测,避免不必要的损失。 -
本文对西藏自治区贡觉县雄松乡至沙东乡的金沙江流域范围进行滑坡探测与监测,获取该区域的LOS向、坡向和垂直滑坡向年形变速率,沿坡向最大年形变速率为-163 mm/a,探测结果较可靠。并通过谷歌地球光学遥感影像对滑坡群进行判别,发现了一个包含7处滑坡隐患点(A~G)的滑坡群。该滑坡群形成的原因包括:(1)研究区域海拔为4 000 m左右,距离冰川近,近年来气温上升,冰川融水增加,导致以冰川为驱动力的滑坡数量增加;(2)在2018-10-11和2018-11-03,研究区域上游白格村附近连续两次发生大规模高位滑坡形成堰塞湖,在人工泄洪时,对该滑坡群前缘带来较大影响。探测到7个隐患点具有较明显的形变,区域A、B、C、D为较活跃滑坡,可能会整体失稳造成堵江后果,目前区域E、F、G滑坡仍处于缓慢变形阶段。GACOS辅助下InSAR Stacking技术可有效去除长波段及与地形相关大气延迟误差,具有简单、有效、快速、易于推广应用等优点。
-
摘要: 西藏自治区贡觉县雄松乡至沙东乡金沙江流域作为川藏铁路的必经河流,地形崎岖、地质灾害隐患点多,亟需对该地区隐患点进行全方位的识别。首先,选取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个区域进行实时监测。Abstract:
Objectives The Sichuan-Tibet Railway crosses the Jinsha River specifically the section from Xiongsong Town to Shadong Town, Gongjue County, Tibet Autonomous Region, in which exist strong topography variations and many potential geohazards, and hence it urgent to detect the potential geohazards in this region to ensure the railway safety. Methods In this paper, 61 Sentinel-1 ascending, 53 Sentinel-1 descending and 7 advanced land observing satellite 2 (ALOS-2) ascending images are used to derive the annual mean deformation rates in the satellite radar line of sight (LOS) with two advanced InSAR(interferometric synthetic aperture radar) approaches, namely generic atmospheric correction online service for InSAR(GACOS)assisted InSAR Stacking and LiCSBAS. The three LOS annual mean deformation rate maps are then used to determine 2D surface movements (one along the slope and the other perpendicular to the slope) in the study region. Results The comparison between GACOS assisted InSAR Stacking and LiCSBAS results show that they agree with each other with correlation coefficients over 0.92 for the LOS deformation rates and over 0.85 for the 2D deformation rates, suggesting the reliability of GACOS assisted InSAR Stacking. The maximum annual deformation rate of 163 mm/a can be observed in the slope direction and seven landslides (A, B, C, D, E, F, G) can be clearly identified, which in turn lays the foundation for future real-time monitoring. Based on the detailed analysis of the seven regions by using InSAR and optical interpretation results, it is found that the seven landslides had relatively obvious deformation. Landslides A, B, C and D are active landslides, which might cause the result of river blocking due to overall instability. At present, landslides E, F and G are still in the stage of slow deformation. Conclusions This study also find that GACOS assisted InSAR Stacking can effectively remove long-band and topographic related atmospheric delay error with the help of GACOS, which has the advantages of being simple, effective, fast, and easy to popularize and apply. GACOS assisted InSAR Stacking technology can be used to quickly identify potential landslide hazards. -
Key words:
- the Jinsha River region /
- landslide detection /
- GACOS /
- InSAR stacking /
- LiCSBAS /
- Sentinel-1 /
- ALOS-2
-
表 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 -
[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 -