[发明专利]融合精密水准与InSAR可靠沉降值的地面沉降监测方法

权利要求书

1.融合精密水准与InSAR可靠沉降值的地面沉降监测方法，其特征在于，包括：

获取两期InSAR影像；

对两期InSAR影像进行差分干涉测量沉降计算，以得到干涉图；

在干涉图的基础上叠加城市不透水面层，得到可靠沉降观测值；

利用精密水准沉降测量结果对所述可靠沉降观测值进行改正与融合，以得到融合沉降观测值；

对融合沉降观测值进行融合点以及面沉降观测值的联合插值计算，以得到计算结果；

根据所述计算结果计算城市地面沉降等值线图。

2.根据权利要求1所述的融合精密水准与InSAR可靠沉降值的地面沉降监测方法，其特征在于，所述对两期InSAR影像进行差分干涉测量沉降计算，以得到干涉图，包括：

匹配覆盖同一地区的两期InSAR影像的对应像元，以得到配准结果；

将配准结果内的像元的复数进行共轭相称，以得到干涉图。

3.根据权利要求1所述的融合精密水准与InSAR可靠沉降值的地面沉降监测方法，其特征在于，所述在干涉图的基础上叠加城市不透水面层，得到可靠沉降观测值，包括：

根据小基线对组合原则对所述干涉图进行干涉滤波，以得到干涉图组合；

对所述干涉图组合去除用外部DEM模拟的地形起伏的相位，以得到差分干涉图；

基于最小费用流法对所述差分干涉图进行缠绕相位的解缠处理，以得到绝对干涉相位；

利用SVD反演算法获取地表形变参数；

利用InSAR主影像卫星轨道数据和外部DEM数据计算R-D坐标系和地理坐标系的映射关系，将两期InSAR影像数据转换到地理坐标系，以得到可靠沉降观测值。

4.根据权利要求3所述的融合精密水准与InSAR可靠沉降值的地面沉降监测方法，其特征在于，所述根据小基线对组合原则对所述干涉图进行干涉滤波，以得到干涉图组合，包括：

根据所述干涉图计算局部坡度以及干涉相位的平均值；

计算所述干涉图的功率谱，依据功率谱设计滤波器，根据小基线对组合原则对所述干涉图进行滤波，形成干涉图组合。

5.根据权利要求1所述的融合精密水准与InSAR可靠沉降值的地面沉降监测方法，其特征在于，所述利用精密水准沉降测量结果对所述可靠沉降观测值进行改正与融合，以得到融合沉降观测值，包括：

制定衡量指标；

基于已知精密观测沉降站点的精密水准沉降测量结果，采用引入地理空间权重矩阵的GM(1，3)预测模型计算未知沉降站点沉降量；

获取目标区域高程变化的干涉影像，将水准测量中选取的水准点通过克里金差值方法生成高程变化的影像，并将目标区域高程变化的干涉影像与生成的高程变化的影像进行叠加和相减，得到处理结果；

采用指数模型对所述处理结果进行回归分析和拟合，以得到改正模型；

利用改正模型对所述可靠沉降观测值进行改正，以得到融合沉降观测值。

6.根据权利要求1所述的融合精密水准与InSAR可靠沉降值的地面沉降监测方法，其特征在于，所述对融合沉降观测值进行融合点以及面沉降观测值的联合插值计算，以得到计算结果，包括：

利用反距离加权平均法、克里金和样条函数中至少一种方法对所述融合沉降观测值进行空间插值计算，以得到计算结果。

7.根据权利要求1所述的融合精密水准与InSAR可靠沉降值的地面沉降监测方法，其特征在于，所述根据所述计算结果计算城市地面沉降等值线图，包括：

对所述计算结果进行采样，以得到采样结果；

对所述采样结果进行提取高程信息，且采用克里金插值对提取高程信息的特征点进行插值处理，以得到栅格数据；

根据所述栅格数据生成等值线，以形成城市地面沉降等值线图。