[发明专利]一种基于高分辨率的时序InSAR形变监测方法及系统

权利要求书

1.一种基于高分辨率的时序InSAR形变监测方法，其特征在于，包括：

获取被监测区域时间序列上的多个地表SAR图像；

对所述地表SAR图像进行配准，并对各地表SAR图像进行差分干涉处理，得到多个差分干涉图；

采用外部DEM数据对所述差分干涉图进行去高程处理，得到多个第一图像；

对所述第一图像中的PS点进行识别；

对所述第一图像中点密度超过设定阈值的PS点进行抽稀处理；

对抽稀后的PS点构建三角网；

对所述三角网中相邻的PS点的相位进行差分，通过解空间搜索的方法得到相邻的PS点之间的形变速率差和高程改正量差，并根据所述形变速率差和高程改正量差计算得到各PS点的残差相位；

从所述残差相位中分离出大气延迟相位，并对所述大气延迟相位进行插值处理；

从所述第一图像中除去所述大气延迟相位，得到第二图像；以所述第二图像替换所述第一图像，跳转至“对抽稀后的PS点构建三角网”步骤，直至满足设定条件后，停止跳转，得到多次去除大气延迟相位后的图像，记为第三图像；

对各差分干涉图对应的第三图像进行相位解缠；

根据相位解缠后的各第三图像中所有PS点的形变，确定地表的形变速率。

2.根据权利要求1所述的基于高分辨率的时序InSAR形变监测方法，其特征在于，所述对所述第一图像中的PS点进行识别，具体包括：

采用振幅离差法或相干系数法对所述第一图像中的PS点进行识别。

3.根据权利要求2所述的基于高分辨率的时序InSAR形变监测方法，其特征在于，在采用所述振幅离差法对所述第一图像中的PS点进行识别时，选取振幅离差小于0.6的点作为PS点。

4.根据权利要求1所述的基于高分辨率的时序InSAR形变监测方法，其特征在于，所述对所述第一图像中点密度超过设定阈值的PS点进行抽稀处理，具体包括：

确定各所述PS点的点密度；

抽取设定区域内点密度超过设定阈值的PS点中时间序列相干均值相对高的PS点。

5.根据权利要求1或4所述的基于高分辨率的时序InSAR形变监测方法，其特征在于，所述设定阈值为50。

6.根据权利要求1所述的基于高分辨率的时序InSAR形变监测方法，其特征在于，所述从所述残差相位中分离出大气延迟相位，具体包括：

对所述残差相位进行多视处理；

采用空间滤波法对多视处理后的残差相位进行滤波；

对滤波后的残差相位进行插值处理。

7.一种基于高分辨率的时序InSAR形变监测系统，其特征在于，包括：

图像获取模块，用于获取被监测区域时间序列上的多个地表SAR图像；

图像预处理模块，用于对所述地表SAR图像进行配准，并对各地表SAR图像进行差分干涉处理，得到多个差分干涉图；

去高程处理模块，用于采用外部DEM数据对所述差分干涉图进行去高程处理，得到多个第一图像；

PS点识别模块，用于对所述第一图像中的PS点进行识别；

PS点抽稀模块，用于对所述第一图像中点密度超过设定阈值的PS点进行抽稀处理；

三角网构建模块，用于对抽稀后的PS点构建三角网；

相邻点相位差分模块，用于对所述三角网中相邻的PS点的相位进行差分，通过解空间搜索的方法得到相邻的PS点之间的形变速率差和高程改正量差，并根据所述形变速率差和高程改正量差计算得到各PS点的残差相位；

大气延迟相位分离模块，用于从所述残差相位中分离出大气延迟相位，并对所述大气延迟相位进行插值处理；

大气延迟相位去除模块，用于从所述第一图像中除去所述大气延迟相位，得到第二图像；以所述第二图像替换所述第一图像，跳转至所述三角网构建模块，直至满足设定条件后，停止跳转，得到多次去除大气延迟相位后的图像，记为第三图像；

相位解缠模块，用于对各差分干涉图对应的第三图像进行相位解缠；

形变速率确定模块，用于根据相位解缠后的各第三图像中所有PS点的形变，确定地表的形变速率。