没有合适的资源?快使用搜索试试~ 我知道了~
利用GACOS辅助下InSAR Stacking对金沙江流域进行滑坡监测.docx
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 1 下载量 40 浏览量
2022-11-29
17:43:14
上传
评论
收藏 973KB DOCX 举报
温馨提示
试读
13页
利用GACOS辅助下InSAR Stacking对金沙江流域进行滑坡监测.docx
资源推荐
资源详情
资源评论
滑坡作为最为常见的地质灾害之一,对经济、社会、人员安全等方面带来巨大损失
[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)年形变速率,进而提取出二维的
形变场,并针对研究区域开展了滑坡探测、监测和对比分析。
1. 研究区概况和研究数据
本文研究区域位于西藏自治区贡觉县雄松乡(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 升
轨影像用于研究。
2. 数据处理
2.1 干涉处理
本文采用 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
剩余12页未读,继续阅读
资源评论
- xingwenglianbo2023-07-14发现一个超赞的资源,赶紧学习起来,大家一起进步,支持!
罗伯特之技术屋
- 粉丝: 3591
- 资源: 1万+
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功