2018年10月11日上午7:00,位于川藏交界带的金沙江发生山体滑坡,形成堰塞湖,地理位置为东经98°42′17.98″,北纬31°4′56.41″。为分析金沙江区域滑坡前后的地表形变情况,基于Sentinel1数据,以DInSAR方法对此次滑坡事件进行干涉测量。
图1 图中红圈位置为主要滑坡点
此次下载了金沙江沿岸的Sentinel1数据,地理范围为:(96.591962,32.697645,98.431692,33.317714,101.853589,27.681709,103.8221,29.226344,104.944335,28.302296,101.743205,25.431675,98.744446,27.812658)。同时手动下载并拼接对应区域dem数据(因为ENVI自动下载DEM时报错)。数据下载链接:https://scihub.copernicus.eu/s3/#/home下sentinel数据或https://vertex.daac.asf.alaska.edu/,https://dds.cr.usgs.gov/srtm/version2_1/SRTM3/Eurasia/下dem数据。数据下载太慢可尝试:在mundi网络里,下载matelink文件,然后用jdownloader下载,速度就快多了。
按照数据处理流程对金沙江沿岸的Sentinel1数据进行处理。主要包括以下几个步骤:
DEM数据导入:涉及格式转换问题,参考http://blog.sciencenet.cn/blog-791897-1083342.html
基线估算
DInSAR工作流
结果分析: 经处理,得到金沙江沿岸地质形变图,可为灾后防范提供预警。