[发明专利]基于InSAR反演与多影响因子的大范围滑坡形变预测方法

权利要求书

1.一种基于InSAR反演与多影响因子的大范围滑坡形变预测方法，其特征在于：包括如下步骤：

S1.采集目标区域的SAR影像数据，并对SAR影像数据进行InSAR反演处理，得到时序形变数据；

S2.对时序形变数据进行聚类处理，得到若干类别的时序形变数据；

S3.将各类别的时序形变数据分解为周期项形变序列与趋势项形变序列；

S4.提取若干影响因子，计算各影响因子分别与各类别的周期项形变序列之间的关联度，选取出与周期项形变显著相关的影响因子，并将显著相关的影响因子作为目标影响因子；

S5.对各类别的时序形变数据分别建立LSTM模型，具体包括：

同一类别的趋势项形变，将趋势项形变作为变量输入，建立趋势项单变量LSTM模型；

同一类别的周期项形变，将周期项形变与目标影响因子作为变量输入，建立周期项多变量LSTM模型；

S6.对步骤S5中的单变量LSTM模型以及多变量LSTM模型进行网络模型训练，得到训练好的LSTM模型，并使用训练好的LSTM模型对各类别的形变进行预测，得到各类别的周期项形变预测值与趋势项形变预测值；

S7.将各类别的周期项形变预测值与趋势项形变预测值相加，得到各类别的形变量预测结果，并将各类别的形变量预测结果进行合并，得到大范围的滑坡形变预测结果。

2.根据权利要求1所述的基于InSAR反演与多影响因子的大范围滑坡形变预测方法，其特征在于：对SAR影像数据进行InSAR反演处理，具体包括：

S11.对SAR影像数据进行预处理；所述预处理包括数据导入以及数据裁剪；

S12.从SAR影像中选取超级主影像，并将其他所有影像配准至超级主影像；

S13.通过基线估算、连接图生成、差分干涉、干涉图滤波以及相位解缠处理，得到SBAS反演估算用的基础数据；

S14.从基础数据中去除残余地形误差与残余大气相位，对形变信息进行单位转换，得到形变结果；所述形变结果包括形变速率与形变时间序列；

S15.将形变结果转到地理坐标系下，得到目标区域的若干个高相干点，每个高相干点有若干时间段的时序形变数据。

3.根据权利要求1所述的基于InSAR反演与多影响因子的大范围滑坡形变预测方法，其特征在于：对时序形变数据进行聚类处理，具体包括：

S21.在所有时序形变数据中随机确定N个数据分别作为聚类中心；

S22.将每一条时序形变数据分别与N个聚类中心计算余弦相似度值CS，并将该条数据分配到CS值最高的聚类中；

S23.各个聚类中的每条数据分别与本聚类中的其他数据计算CS值，将与所述其他数据的平均CS值最高的数据重新定义为聚类中心；

S24.重复步骤S22-S23进行迭代，直到各聚类中的数据不再变化或迭代到K次为止；

S25.计算各聚类中心数据的CS值，若CS值大于A，则将类别数N加1，并重新聚类；若CS值小于B，则将类别数N减1，并重新聚类；聚类完成后分别输出各类别数据。

4.根据权利要求3所述的基于InSAR反演与多影响因子的大范围滑坡形变预测方法，其特征在于：根据如下公式计算余弦相似度值CS：

其中，与分别代表两条不同的时序数据，x₀,…,x_N-1以及y₀,…,y_N-1分别代表N个时序的形变。

5.根据权利要求1所述的基于InSAR反演与多影响因子的大范围滑坡形变预测方法，其特征在于：步骤S3中，所述分解采用小波分解，所述小波的分解函数为Daubechies，所述小波的分解层数为四层。

6.根据权利要求1所述的基于InSAR反演与多影响因子的大范围滑坡形变预测方法，其特征在于：所述若干影响因子包括降雨量、库水位高程以及库水位变化量。

7.根据权利要求1所述的基于InSAR反演与多影响因子的大范围滑坡形变预测方法，其特征在于：根据如下步骤计算各影响因子分别与各类别的周期项形变序列之间的关联度：

S41.将周期项形变序列取作为将n个影响因子序列取作并对周期项形变序列与影响因子序列进行归一化处理；

S42.计算周期项形变序列与影响因子序列对应元素的绝对差值|x₀(t)-x_i(t)|；其中，x₀(t)为周期项形变序列在t时刻的值，x_i(t)为影响因子序列在t时刻的值；

S43.计算关联系数ξ_0i(t)：

其中，ξ_0i(t)为周期项形变序列与影响因子序列在t时刻的关联系数，表示两序列的最小绝对差值，表示两序列的最大绝对差值，ρ为分辨系数；

S44.计算关联度r_0i：

其中，n为时序序列长度。