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

说明书

本发明提供了一种地表形变监测与预测方法，包括：获取研究区域的InSAR时间序列地表形变数据集，并对其进行处理，获取研究区域的环境因素数据集，并对其进行处理，将处理后的环境因素数据集及处理后的InSAR时间序列地表形变数据集进行合并，得到模型数据集，搭建Autoformer模型，通过模型数据集对Autoformer模型进行训练，调整模型超参数，获取最佳模型，选取研究区域内地表形变变化最大的研究点集，输入最佳模型进行预测预警。本发明提供的地表形变监测与预测方法，无需地面及人工配合，能够充分发挥卫星监测的优点，可对研究区域的全范围进行地表形变监测，且能够实现预测，能够降低事故发生率及财产损失率。

技术领域

本发明涉及地表形变监测技术领域，特别是涉及一种地表形变监测与预测方法。

背景技术

在地表形变的监测手段中，传统的大地测量技术，例如GPS及人工水准仪等容易受到低空间分辨率以及高成本的限制，另外，在对地物的巡检中，传统技术往往需要大量的专业人员及重型器械，监测时间长且监测成本高。

现有技术中还使用多时相干涉合成孔径雷达(MT-InSAR)技术实现大规模、高精度、低成本的监测，其中，典型的MT-InSAR技术之一是永久散射雷达干扰合成孔径雷达(PS-InSAR)，在该方法中，选择了建筑、裸露岩石等时间稳定点，这些点在多幅SAR图像中具有高度稳定性和反射性，但是现有技术中对PS-InSAR进行使用时往往未将其与长期时间序列预测处理分析框架相配合使用，不能实现密集地物地表沉降的预测，而且并未考虑影响地表形变的环境因素。因此，设计一种基于Autoformer及PS-Insar的地表形变监测与预测方法是十分有必要的。

发明内容

本发明的目的是提供一种地表形变监测与预测方法，无需地面及人工配合，能够充分发挥卫星监测的优点，可对研究区域的全范围进行地表形变监测，且能够实现预测，能够降低事故发生率及财产损失率。

为实现上述目的，本发明提供了如下方案：

一种地表形变监测与预测方法，包括如下步骤：

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

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

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

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

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

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

可选的，步骤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点，建立形变反演模型，为：