[发明专利]融合北斗和InSAR数据的地表沉降监测方法及装置

权利要求书

1.一种融合北斗和InSAR数据的地表沉降监测方法，其特征在于，包括：

获取监测区域内监测点的北斗观测数据和监测区域内永久散射体PS点的合成孔径雷达SAR影像数据；

分别对所述北斗观测数据和所述SAR影像数据进行沉降时序解算，得到监测点的北斗沉降时序和PS点的SAR沉降时序；

将所述PS点中与所述监测点重合的公共点的北斗沉降时序和公共点的SAR沉降时序进行融合，得到公共点的融合沉降时序；

基于所述公共点的融合沉降时序，构建监测区域的沉降场曲面方程；

基于所述公共点的融合沉降时序和所述沉降场曲面方程，对所述PS点中未与所述监测点重合的其他PS点进行时空插值计算，得到其他PS点的融合沉降时序，以基于所述公共点的融合沉降时序和所述其他PS点的融合沉降时序，实现对监测区域的监测。

2.根据权利要求1所述的方法，其特征在于，分别对所述北斗观测数据和所述SAR影像数据进行沉降时序解算，得到监测点的北斗沉降时序和PS点的SAR沉降时序，包括：

基于所述北斗观测数据，利用目标方法求取监测点的原始沉降时序，其中，所述目标方法包括：实时载波相位差分方法和精密单点定位方法；

对所述监测点的原始沉降时序进行第一预处理操作，得到所述监测点的北斗沉降时序；其中，所述第一预处理操作包括：粗差探测处理操作和剔除处理操作。

3.根据权利要求1所述的方法，其特征在于，分别对所述北斗观测数据和所述SAR影像数据进行沉降时序解算，得到监测点的北斗沉降时序和PS点的SAR沉降时序，包括：

在对所述SAR影像数据进行配准的过程中，按照预设稳定性条件选取PS点；

基于所述PS点对所述SAR影像数据进行第二预处理操作，得到PS点的形变相位信息；所述第二预处理操作包括：滤波操作，干涉操作和去地形处理操作；

基于所述PS点的形变相位信息，确定所述PS点的SAR沉降时序。

4.根据权利要求1所述的方法，其特征在于，将所述PS点中与所述监测点重合的公共点的北斗沉降时序和公共点的SAR沉降时序进行融合，得到公共点的融合沉降时序，包括：

确定所述PS点中与所述监测点重合的公共点；

利用卡尔曼滤波算法对公共点的北斗沉降时序和公共点的SAR沉降时序进行融合，得到所述公共点的融合沉降时序。

5.根据权利要求1所述的方法，其特征在于，基于所述公共点的融合沉降时序，构建监测区域的沉降场曲面方程，包括：

基于所述公共点的融合沉降时序，提取所述公共点的属性值；其中，所述属性值包括平面坐标和沉降值；

基于所述公共点的平面坐标和沉降值，进行最小二乘法曲面拟合，得到监测区域的沉降场曲面方程。

6.根据权利要求1所述的方法，其特征在于，基于所述公共点的融合沉降时序和所述沉降场曲面方程，对所述PS点中未与所述监测点重合的其他PS点进行时空插值计算，得到其他PS点的融合沉降时序，包括：

从所述PS点中筛除与所述监测点重合的公共点，得到其他PS点；

从所述PS点的SAR沉降时序中获取所述其他PS点的SAR沉降时序；

基于所述其他PS点的SAR沉降时序和所述公共点的融合沉降时序，预测其他PS点在InSAR监测时刻的融合沉降时序；

基于所述沉降场曲面方程，预测所述其他PS点在非InSAR监测时刻的融合沉降时序；

将所述其他PS点在InSAR监测时刻的融合沉降时序和/或所述其他PS点在非InSAR监测时刻的融合沉降时序确定为其他PS点的融合沉降时序。