一种基于多源对地观测的活动性地质灾害探测方法与系统技术方案

本发明专利技术公开了一种基于多源对地观测的活动性地质灾害探测方法与系统，涉及地质灾害探测技术领域，包括如下步骤：获取卫星遥感影像；联合

【技术实现步骤摘要】
一种基于多源对地观测的活动性地质灾害探测方法与系统


[0001]本专利技术涉及地质灾害探测
，特别涉及一种基于多源对地观测的活动性地质灾害探测方法与系统
。

技术介绍

[0002]地质灾害是在地球表面普遍发生的地质现象
,
地质灾害的发生经常造成人员伤亡和经济损失
。
为减少地质灾害风险，查明活动性地质灾害的清单至关重要
。
[0003]针对几十万平方公里的区域，只有基于星载遥感影像的
InSAR
技术和
POT
技术可以以低成本
、
高效率的方式获取广域地表变形，特别是在野外调查具有挑战性或不可行的危险区域，已成为探测地质灾害不可或缺的工具
。
[0004]以往研究大多数采用单一技术手段获取大范围地表形变，但事实上不同阶段的地质灾害具有不同梯度的形变信号，移动速度每年达到米级形变的地质灾害常被忽视
。
尤其是在几十万平方公里的大区域情况下，目前获取形变信息后活动性地质灾害的位置和范围确定依赖于人工解译，耗时耗力
。

技术实现思路

[0005]本专利技术的目的在于针对上述现有技术的不足，提供一种基于多源对地观测的活动性地质灾害探测方法，以解决现有技术中移动速度每年达到米级形变的地质灾害常被忽视，且获取形变信息后活动性地质灾害的位置和范围确定依赖于人工解译，耗时耗力的问题
。
[0006]本专利技术具体提供如下技术方案：一种基于多源对地观测的活动性地质灾害探测方法，包括以下步骤：
[0007]获取若干卫星遥感影像；
[0008]对所述卫星遥感影像进行优化处理，获得毫米到米级形变信息；
[0009]联合
Enhanced InSAR Stacking
技术
、SAR POT、OPOT
和
DEM
差分技术对所述卫星遥感影像进行处理，获得毫米到米级形变信息；
[0010]利用
MCD
方法统计识别毫米到米级形变信息中的运动像素；
[0011]利用
DBSCAN
方法对所述运动像素进行处理，圈定活动性地质灾害位置和边界，获得活动性地质灾害范围；
[0012]采集所述活动性地质灾害范围内的平均坡度，利用所述平均坡度将活动性地质灾害划分为活动性滑坡或地面沉降
。
[0013]优选的，所述卫星遥感影像包括：多轨道多波段的合成孔径雷达影像
、
多期高空间分辨率的卫星遥感影像和多期高空间分辨率的
DEM
影像
。
[0014]优选的，所述对所述卫星遥感影像进行优化处理，获得毫米到米级形变信息，包括如下步骤：
[0015]联合
Enhanced InSAR Stacking
技术
、SAR POT、OPOT
和
DEM
差分技术对所述卫星遥
感影像进行处理，获得毫米到米级形变信息
。
[0016]优选的，使用
Enhanced InSAR Stacking
技术对所述卫星遥感影像进行处理，包括如下步骤：
[0017]采取
InSAR Stacking
与大气改正算法的策略减少卫星遥感影像中非平稳信号带来的估计偏差；
[0018]利用减少估计偏差后的卫星遥感影像和辅助数据生成干涉图；
[0019]通过自适应滤波抑制干涉图中的随机噪声；
[0020]对滤波后的干涉图利用最小费用流算法
MCF
进行相位解缠；
[0021]校正相位解缠后干涉图中的大气延迟，并利用最小二乘法估计校正相位解缠后干涉图中每个点的平均速度，获得最终的地表形变信息
。
[0022]优选的，使用
SAR POT、OPOT
技术对所述卫星遥感影像进行处理，包括如下步骤：
[0023]所述
SAR POT
与
OPOT
通过计算搜索窗口中合成孔径雷达影像块之间的相关性，对多时相影像中移动目标进行追踪，捕捉到整个像素与亚像素级的偏移量，获得最终的地表形变信息
。
[0024]优选的，使用
DEM
差分技术对所述卫星遥感影像进行处理，包括如下步骤：
[0025]所述
DEM
差分技术为两期卫星遥感影像的
DEM
影像精配准后作差，获得最终的地表形变信息
。
[0026]优选的，所述毫米到米级形变信息包括：
Sentinel
‑1升轨
LOS
向年形变速率
、Sentinel
‑1降轨
LOS
向年形变速率
、Sentinel
‑1升轨
Range
向形变量
、Sentinel
‑1降轨
Range
向形变量
、Sentinel
‑2获取的东西向形变和
Sentinel
‑2获取的南北向形变
。
[0027]优选的，所述利用
DBSCAN
方法对所述运动像素进行处理，圈定活动性地质灾害位置和边界，包括如下步骤：
[0028]使用
DBSCAN
方法将所述
MCD
确定的运动像素进行聚类；
[0029]将聚类后具有足够密度的区域划分为簇，将所述簇的密度阈值参数
MinPts
设置为3；
[0030]将所述密度阈值参数
MinPts
大于等于3的簇定义为活动性地质灾害；
[0031]将所述簇的参数
Eps
设置为2个像素大小，将距离小于2个像素的簇进行连接；
[0032]利用凸包算法对设置参数后
、
且符合条件的簇进行边界的绘制，获得最终的活动性地质灾害位置和边界
。
[0033]优选的，所述利用所述平均坡度将活动性地质灾害划分为活动性滑坡和地面沉降，包括如下步骤：
[0034]设置所述平均坡度的阈值为5°
，将平均坡度大于阈值的活动性地质灾害划分为活动性滑坡，将平均坡度小于阈值的活动性地质灾害；
[0035]在活动性滑坡速度
V<100mm/yr
时，划定所述活动性滑坡为一级滑坡；
[0036]在活动性滑坡速度
V>100mm/yr
时，划定所述活动性滑坡为二级滑坡；
[0037]其中，所述二级滑坡的速度大于一级滑坡的速度
。
[0038]本专利技术还提供一种基于多源对地观测的活动性地质灾害探测系统，包括：
[0039]图像获取模块，用于获取若干卫星遥感影像；对卫星遥感影像进行优化处本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种基于多源对地观测的活动性地质灾害探测方法，其特征在于，包括以下步骤：获取若干卫星遥感影像；对所述卫星遥感影像进行优化处理，获得毫米到米级形变信息；利用
MCD
方法统计识别毫米到米级形变信息中的运动像素；利用
DBSCAN
方法对所述运动像素进行处理，圈定活动性地质灾害位置和边界，获得活动性地质灾害范围；采集所述活动性地质灾害范围内的平均坡度，利用所述平均坡度将活动性地质灾害划分为活动性滑坡或地面沉降
。2.
如权利要求1所述的一种基于多源对地观测的活动性地质灾害探测方法，其特征在于，所述若干卫星遥感影像包括：多轨道多波段的合成孔径雷达影像
、
多期高空间分辨率的卫星遥感影像和多期高空间分辨率的
DEM
影像
。3.
如权利要求2所述的一种基于多源对地观测的活动性地质灾害探测方法，其特征在于，所述对所述卫星遥感影像进行优化处理，获得毫米到米级形变信息，包括如下步骤：联合
Enhanced InSAR Stacking
技术
、SAR POT、OPOT
和
DEM
差分技术对所述卫星遥感影像进行处理，获得毫米到米级形变信息
。4.
如权利要求3所述的一种基于多源对地观测的活动性地质灾害探测方法，其特征在于，使用
Enhanced InSAR Stacking
技术对所述卫星遥感影像进行处理，包括如下步骤：采取
InSAR Stacking
与大气改正算法的策略减少卫星遥感影像中非平稳信号带来的估计偏差；利用减少估计偏差后的卫星遥感影像和辅助数据生成干涉图；通过自适应滤波抑制干涉图中的随机噪声；对滤波后的干涉图利用最小费用流算法
MCF
进行相位解缠；校正相位解缠后干涉图中的大气延迟，并利用最小二乘法估计校正相位解缠后干涉图中每个点的平均速度，获得最终的地表形变信息
。5.
如权利要求3所述的一种基于多源对地观测的活动性地质灾害探测方法，其特征在于，使用
SAR POT、OPOT
技术对所述卫星遥感影像进行处理，包括如下步骤：所述
SAR POT
与
OPOT
通过计算搜索窗口中合成孔径雷达影像块之间的相关性，对多时相影像中移动目标进行追踪，捕捉到整个像素与亚像素级的偏移量，获得最终的地表形变信息
。6.
如权利要求3所述的一种基于多源对地观测的活动性地质灾害探测方法，其特征在于，使用
DEM
差分技术对所述卫星遥感影像进行处理，包括如下步骤：所述
DEM
差分技术将两期卫星遥感影像的
DEM
影像精配准后作差，获得最终的地表形变信息
...

【专利技术属性】
技术研发人员：李振洪，陈博，宋闯，余琛，张成龙，杜建涛，
申请(专利权)人：长安大学，
类型：发明
国别省市：