[发明专利]一种基于InSAR的填海地区机场形变监测方法

权利要求书

1.一种基于InSAR的填海地区机场形变监测方法，其特征在于，包括如下步骤：

S1、根据预设的振幅离差阈值提取填海地区机场所在区域中离散分布的候选点目标；

S2、根据所述离散分布的候选点目标计算每个候选点目标的相位标准差，得到稳定点目标，所述步骤S2包括如下步骤：

S201、利用离散分布的候选点目标构建空间滤波格网，去除空间相关的相位贡献；

S202、利用最邻近差值计算候选点空间不相关的地形残余相位和误差相位，得到地形残余相位和时间相干系数；

S203、根据所述时间相干系数，利用统计法得到候选点目标的时间相干系数阈值，并选出相干点目标；

S204、判断所述时间相干系数是否小于所述时间相干系数阈值，若是，则进入步骤S206，反之，则进入步骤S205；

S205、根据所述候选点空间不相关的误差相位，计算相干点目标的信噪比，并通过信噪比更新空间滤波格网中的网格相位权值，返回步骤S201；

S206、根据判断结果计算每个候选点的相位标准差；

S207、判断所述候选点相位标准差是否小于预设的最小相位标准差阈值，若小于，则进入步骤S3，反之，则去除错误的候选点目标，从而完成对稳定点目标的提取；

S3、将所述稳定点目标的相位进行三维时空解缠，得到解缠相位；

S4、根据所述解缠相位利用最小二乘法得到包含有空间干扰分量的形变相位；

S5、根据所述包含有空间干扰分量的形变相位分离出大气相位和轨道误差相位，得到形变相位；

S6、根据所述形变相位利用地表形变与雷达成像的关系，计算形变相位的权值以及稳定点目标的时间序列形变，从而完成使用时间序列InSAR技术对填海地区机场形变的监测。

2.根据权利要求1所述的基于InSAR的填海地区机场形变监测方法，其特征在于，所述步骤S205中空间滤波格网中的网格相位权值的表达式如下：

其中，φ表示空间滤波格网中的网格相位权值，ADI_i表示振幅离差，n表示候选点的总数，表示初始的权值，φ_i表示第i个候选点的相位值。

3.根据权利要求1所述的基于InSAR的填海地区机场形变监测方法，其特征在于，所述步骤S4中包含有空间干扰分量的形变相位的计算公式如下：

其中，表示包含有空间干扰分量的形变相位，表示干涉相位中去除地形相位后的剩余相位，表示大气相位与轨道误差相位之和，表示噪声相位。

4.根据权利要求1所述的基于InSAR的填海地区机场形变监测方法，其特征在于，所述步骤S5具体为：

根据包含有空间干扰分量的形变相位，采用时间域上的高通滤波和空间域上的低通滤波分离出大气相位和轨道误差相位，得到形变相位。

5.根据权利要求1所述的基于InSAR的填海地区机场形变监测方法，其特征在于，所述步骤S6中的形变相位的权值计算公式如下：

P＝diag(σ¹，σ²，…σ^N)

其中，P表示形变相位的权值,diag(·)表示对角矩阵函数，σ^k表示第k干涉影像对的干涉相位标准差,N表示干涉影像对的数量，表示第x个稳定点目标在第k个干涉影像对上的噪声相位，M表示稳定点目标的数量。

6.根据权利要求1所述的基于InSAR的填海地区机场形变监测方法，其特征在于，所述步骤S6中的候选点目标的时间序列形变v的表达式如下：

v＝(T'PT)^-1T'Pd

其中，T'表示候选点目标的时间基线矩阵的转置矩阵，P表示形变相位的权矩阵，T表示候选点目标的时间基线矩阵，d表示候选点目标的形变量。