[发明专利]一种地表形变监测与预测方法

权利要求书

1.一种地表形变监测与预测方法，其特征在于，包括如下步骤：

步骤1：获取研究区域的InSAR时间序列地表形变数据集；

步骤2：对获取的InSAR时间序列地表形变数据集进行处理；

步骤3：获取研究区域的环境因素数据集；

步骤4：对环境因素数据集进行处理，并将处理后的环境因素数据集及处理后的InSAR时间序列地表形变数据集进行合并，得到模型数据集；

步骤5：搭建Autoformer模型，通过模型数据集对Autoformer模型进行训练，调整模型超参数，获取最佳模型；

步骤6：选取研究区域内地表形变变化最大的研究点集，输入最佳模型进行预测预警。

2.根据权利要求1所述的地表形变监测与预测方法，其特征在于，步骤1中，获取研究区域的InSAR时间序列地表形变数据集，具体为：

获取覆盖研究区域的n幅Sentinel-1A的SAR影像，根据时空基线参数及多普勒中心频率，选取其中一幅SAR影像作为主影像，其余作为辅影像，对主影像及辅影像进行精确配准及干涉处理，生成N-1幅差分干涉图，获取外部DEM数据及卫星轨道数据，使用外部DEM数据与影像覆盖范围进行配准，通过卫星轨道数据进行轨道偏移矢量矫正去除平面效应，并将DEM数据采样到SAR影像的坐标系中，通过振幅离差法识别PS点，设定合理的振幅离差阈值，提取研究区域PS候选点，构建像素点网络，其中，振幅离差法通过评估影像像元的信噪比来确认PS点，相位的信噪比可用振幅离差指数来衡量，即：

式中，σ_A及μ_A分别为SAR影像的振幅标准差和振幅平均值，设定D_A的阈值为x，当该位置振幅离差指数小于x时，则判断为PS点，反之，则判断不是PS点，建立形变反演模型，为：

ψ_x,i＝W{φ_D,x,i+φ_A,x,i+Δφ_S,x,i+Δφ_θ,x,i+Δφ_N,x,i}

式中，W为相位缠绕算子，φ_D,x,i为由地面相位引起的相位，φ_A,x,i为大气延迟误差相位，Δφ_S,x,i为卫星轨道误差相位，Δφ_θ,x,i为残余地形相位，Δφ_N,x,i为配准误差，即热噪声失相关引起的噪声相位，建立形变反演模型，进行相邻PS点相位参数估计，解算残余高程和线性形变及大气相位参数，并测试模型解算的各参数与像素点干涉相位的拟合情况，证明模型解算值的准确性，进行大气相位移除后，重新选取PS候选点，再次进行形变分析，最后对残差进行分解，得到非线性形变信息，获取时间序列地表形变数据集，经过地理编码后得到InSAR时间序列地表形变数据集。

3.根据权利要求2所述的地表形变监测与预测方法，其特征在于，步骤2中，对获取的InSAR时间序列地表形变数据集进行处理，具体为：

将获取的InSAR时间序列地表形变数据集按照研究区域进行裁剪和平滑化处理，获取时间序列数据集长度，并进行文件格式转换，确认研究区域的关键点，即形变程度最大的点集合，得到处理后的InSAR时间序列地表形变数据集。

4.根据权利要求3所述的地表形变监测与预测方法，其特征在于，步骤3中，获取研究区域的环境因素数据集，具体为：

获取NetCDF格式数据的大气驱动数据集、地表温度分析数据集、土壤湿度数据集、土壤温度分析数据集及土壤相对湿度数据集，组成环境因素数据集。

5.根据权利要求4所述的地表形变监测与预测方法，其特征在于，步骤4中，对环境因素数据集进行处理，并将处理后的环境因素数据集及处理后的InSAR时间序列地表形变数据集进行合并，得到模型数据集，具体为：

通过基于Python语言的Pandas库及Numpy库将NetCDF格式数据的环境因素数据集转化为csv格式的环境因素数据集，并抽取环境因素数据集中的环境因素参数，将环境因素数据集按照研究区域经纬度网格划分，并对抽取的环境因素参数进行数据格式标准化，完成标准化后，对环境因素数据集进行数据格式审查，去除不符合常规格式的数据噪声，去除完毕后，对环境因素数据集进行异常值补充，使用各环境因素参数的平均值对环境因素数据集的空值进行补充，得到处理过的环境数据因素数据集，计算研究区域内的各研究点与采集环境因素参数的环境因素气象站的距离，每个研究点选取距其最近的环境因素气象站采集的环境因素参数作为参考，对各研究点进行聚类分簇处理，合并环境因素数据集及InSAR地表形变时间序列数据集，得到模型数据集。