CN110673145A - 一种基于间断相干的InSAR地表形变监测方法及系统 - Google Patents

Abstract

本发明公开了一种基于间断相干的InSAR地表形变监测方法及系统。该方法包括：获取被监测区域的地表SAR影像；对SAR影像进行配准；采用小基线集法提取SAR影像中的干涉对；对干涉对进行差分干涉处理；提取干涉对中的相干点，相干点包括在时间序列上连续保持高相干性的连续相干点和在时间序列上间断保持高相干性的间断相干点；对相干点进行相位解缠；根据相干点相位解缠后的相位确定相干点的线性形变；从相干点相位解缠后的相位中减去相干点的线性形变相位，得到残余相位；对残余相位进行滤波处理，得到相干点的非线性形变；综合相干点的线性形变和非线性形变，确定被监测区域地表的形变或形变速率。本发明能够对时序上间断相干点位的区域形变进行有效监测。

Description

一种基于间断相干的InSAR地表形变监测方法及系统

技术领域

本发明涉及微波遥感技术领域，特别是涉及一种基于间断相干的InSAR地表形变监测方法及系统。

背景技术

InSAR是一种主动式微波雷达探测技术。它具有大规模，低成本和高精度等技术优势。近年来，它在许多领域显示出巨大的变形监测潜力，如滑坡，冰川运动，地面沉降等。常规的差分干涉测量(Differential SAR Interferomtry,DInSAR)技术通过对两幅影像作差分来计算形变，但是由于DInSAR技术受失相干现象的影响，目前国际上对InSAR技术的研究主要集中在长时间序列DInSAR技术的研究。长时间序列的DInSAR与传统的DInSAR有显著的不同，主要是对长时间范围内相位和幅度变化稳定的点进行分析，利用这些稳定点的相位特征，消除大气延迟成分，获得高精度的地表形变信息，从而实现长时间尺度的形变分析。目前关于长时间序列的DInSAR的研究方法有最小二乘方法、永久散射体方法、小基线集方法(Small Baseline Differential SAR Interferogram，SBAS)和相干目标方法。

Berardino等人(2002)提出的小基线集技术是一种经典的长时间序列的DInSAR方法。通过选择时空基线较小的干涉像对，限制了长基线导致的几何失相关问题，使更多的干涉相对参与计算，增加了时间上的采样，采用最小二乘或者奇异值分解获取地表形变信息，是一种有效的大尺度的形变监测技术。但是这种方法只能监测到时间序列中保持高相干性的区域，如城市区域、裸露岩石和大型人工设施。对于低相干值区域(植被)是无效的，这样不能完整反映形变区域，给后续的形变解译带来困难。由于受到人类活动、地质活动和地质环境的影响，如季节性植被覆盖、土壤湿度和气候变化等，都会突然改变地表的物理特征，降低干涉图的相干性，造成这些像元在时间序列中间断相干，虽然这样的点在时间序列中包含有形变信息，但是在传统的SBAS技术中是不考虑的，以至于在低相干区域点位分布较少，使得这些区域的形变信息无法得到正确的反应。

发明内容

本发明的目的是提供一种基于间断相干的InSAR地表形变监测方法及系统，能够对时序上间断相干点位的区域形变进行有效监测。

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

一种基于间断相干的InSAR地表形变监测方法，包括：

获取被监测区域的地表SAR影像，所述SAR影像包括时间序列上的多个卫星雷达影像；

从所述SAR影像中选取主影像，并根据所述主影像对所述SAR影像进行配准；

采用小基线集法提取所述SAR影像中的干涉对；

对所述干涉对进行差分干涉处理；

提取所述干涉对中的相干点，所述相干点包括连续相干点和间断相干点，所述连续相干点为在时间序列上相干值的均值大于第一预设值的相干点；所述间断相干点为在时间序列上具有多个相干值大于第二预设值且大于所述第二预设值的相干值的数量大于设定值的相干点；

对所述相干点进行相位解缠；

根据所述相干点相位解缠后的相位确定所述相干点的线性形变；

从所述相干点相位解缠后的相位中减去所述相干点的线性形变相位，得到残余相位；

对所述残余相位进行滤波处理，得到所述相干点的非线性形变；

综合所述相干点的线性形变和非线性形变，确定被监测区域地表的形变或形变速率。

可选的，所述对所述干涉对进行差分干涉处理，具体包括：

采用轨道数据和地形数据对所述干涉对进行二轨差分干涉处理。

可选的，所述根据所述相干点相位解缠后的相位确定所述相干点的线性形变，具体包括：

对所述相干点构建三角网；

以三角网中位于同一边两端的相干点的形变速率差分为待求解量，以所述相干点相位解缠后的相位为已知量，构建线性形变求解方程，并约束相邻时期的形变速率差，最后求解得到所述相干点的线性形变。

可选的，所述约束相邻时期的形变速率，具体包括：

将相邻时期形变速率的约束确定为-0.5Δv_k-1+Δv_k-0.5Δv_k+1＝0，其中，Δv_k为k时期的形变速率差。

可选的，所述对所述相干点进行相位解缠，具体包括：

采用最小费用流方法对所述相干点进行相位解缠。

本发明还提供了一种基于间断相干的InSAR地表形变监测系统，包括：

地表SAR影像获取模块，用于获取被监测区域的地表SAR影像，所述SAR影像包括时间序列上的多个卫星雷达影像；

影像配准模块，用于从所述SAR影像中选取主影像，并根据所述主影像对所述SAR影像进行配准；

干涉对提取模块，用于采用小基线集法提取所述SAR影像中的干涉对；

差分干涉处理模块，用于对所述干涉对进行差分干涉处理；

相干点提取模块，用于提取所述干涉对中的相干点，所述相干点包括连续相干点和间断相干点，所述连续相干点为在时间序列上相干值的均值大于第一预设值的相干点；所述间断相干点为在时间序列上具有多个相干值大于第二预设值且大于所述第二预设值的相干值的数量大于设定值的相干点；

相位解缠模块，用于对所述相干点进行相位解缠；

线性形变确定模块，用于根据所述相干点相位解缠后的相位确定所述相干点的线性形变；

残余相位确定模块，用于从所述相干点相位解缠后的相位中减去所述相干点的线性形变相位，得到残余相位；

非线性形变确定模块，用于对所述残余相位进行滤波处理，得到所述相干点的非线性形变；

形变确定模块，用于综合所述相干点的线性形变和非线性形变，确定被监测区域地表的形变或形变速率。

可选的，所述差分干涉处理模块，具体包括：

差分干涉处理单元，用于采用轨道数据和地形数据对所述干涉对进行二轨差分干涉处理。

可选的，所述线性形变确定模块，具体包括：

三角网构建单元，用于对所述相干点构建三角网；

线性形变求解单元，用于以三角网中位于同一边两端的相干点的形变速率差分为待求解量，以所述相干点相位解缠后的相位为已知量，构建线性形变求解方程，并约束相邻时期的形变速率差，最终求解得到所述相干点的线性形变。

可选的，相邻时期形变速率的约束为-0.5Δv_k-1+Δv_k-0.5Δv_k+1＝0，其中，Δv_k为k时期的形变速率差。

可选的，所述相位解缠模块，具体包括：

相位解缠单元，用于采用最小费用流方法对所述相干点进行相位解缠。

根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明提供的基于间断相干的InSAR地表形变监测方法及系统在选取被监测区域SAR影像上的相干点时，既选取了在时间序列上连续保持高相干性的连续相干点，还选取了在时间序列上间断保持高相干性的间断相干点；根据这些相干点的相位变化，确定被监测区域地表的形变或形变速率。由于本发明在提取相干点时，考虑了时序间断相干点，使得本发明能够对时序上间断相干点位的区域形变进行有效监测。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例中基于间断相干的InSAR地表形变监测方法的流程示意图；

图2为本发明实施例中间断相干点示意图；

图3为本发明实施例中的干涉序列示意图；

图4为本发明实施例中示例的处理流程图；

图5为本发明实施例中的干涉基线组合图；

图6为本发明实施例中的相干点分布图；

图7为本发明实施例采用本发明提出的方法得到的形变速率图；

图8为传统小基线集方法得到的形变速率图；

图9为本发明实施例中基于间断相干的InSAR地表形变监测系统的结构示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明的目的是提供一种基于间断相干的InSAR地表形变监测方法及系统，能够对时序上间断相干点位的区域形变进行有效监测。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

图1为本发明实施例中的基于间断相干的InSAR地表形变监测方法的流程示意图，如图1所示，该方法包括以下步骤：

步骤101：获取被监测区域的地表SAR影像，SAR影像包括时间序列上的多个卫星雷达影像；

步骤102：从SAR影像中选取主影像，并根据主影像对SAR影像进行配准；

步骤103：采用小基线集法提取SAR影像中的干涉对；

步骤104：对干涉对进行差分干涉处理；

步骤105：提取干涉对中的相干点，所述相干点包括连续相干点和间断相干点；所述连续相干点为在时间序列上连续保持高相干性的相干点，所述间断相干点为在时间序列上具有多个相干值大于第二预设值且大于第二预设值的相干值的数量大于设定值的相干点；

步骤106：对相干点进行相位解缠；

步骤107：根据相干点相位解缠后的相位确定相干点的线性形变；

步骤108：从相干点相位解缠后的相位中减去相干点的线性形变相位，得到残余相位；

步骤109：对残余相位进行滤波处理，得到相干点的非线性形变；

步骤110：综合相干点的线性形变和非线性形变，确定被监测区域地表的形变或形变速率。

在一实施例中，步骤102中的影像配准方法为：首先以最大化所有干涉像对间的相干性的原则选择参考影像，干涉像对间的相干值依赖像对间的时间基线(T)、空间垂直基线(B_⊥)、多普勒中心基线(F_DC)和热噪声四部分，是四部分的乘积形式，因此从多时相SAR影像中选择主影像必须综合考虑这四部分的分布情况，从中选择使得整体最优的一个影像作为主影像。然后通过卫星轨道参数法对影像进行粗配准，再通过相干系数法对影像进行精配准，在精配准的过程中，需要对影像进行重采样，使每幅图的像元精确配准。

在一实施例中，步骤103中：为了保证干涉相干性，选择时间和空间基线较小的干涉对，但是，不同于传统的SBAS，本发明关注的是间断相干点，因此可以在保证差分干涉相位空间域正确解缠的前提下，放大时空基线，以ERS数据为例，空间基线可以扩大到300m。

在一实施例中，步骤104中的差分干涉方法为：在配准后，依据轨道数据和地形数据对干涉对进行二轨差分干涉，差分干涉过程中，为了减少相位解缠的难度，可采用多视处理。

在一实施例中，步骤105中连续相干点的选取得到可以是：选取相干值的均值大于第一预设值的相干点作为连续相干点，或者可以是选取时序上相干值均大于第一预设值的相干点作为连续相干点，其中第一预设值是人为设定的，可以为0.7或0.8。

在一实施例中，步骤105中所述间断相干点的具体选择方法为：每幅图上选取相干值大于一定阈值(第二预设值)的点，并且这样的点的个数满足一定要求(个数大于设定值)。比如：一个点在每个干涉图中相干性分别为：0.9、0.2、0.3、0.8、0.7，如果选择0.6作为第二预设值，那么有3个图是满足的。如果个数阈值(即设定值)设置为3，那么这个点就满足间断相干点的选取条件，这样，该点就可以选为间断相干点。

在一实施例中，步骤105中识别间断相干点的方法为：假设有N景SAR影像，按照小基线集的干涉组合规则，组成M个干涉对，差分干涉后获得M个相干值图像。传统SBAS方法选取具有高时序相干均值且在时序中保持稳定高相干值的点作为PS点，如图2中的A点。但本发明除了选择A点外，还顾忌了时间序列中间断但保持高相干值的点，如图2中B点所示，这样的点可以通过选定相干值大于这个某个阈值的点得到，在实践应用中，依据SAR数据量，还需合理设置时序中高相干值的个数阈值m(m＜M)。

在一实施例中，步骤106中相位解缠的方法为：采用最小费用流方法进行相位解缠，该算法将相位解缠问题转移到网络规划领域中的最小费用流问题上，用网络规划的方法实现相位解缠。

在一实施例中，步骤107具体包括：

对相干点构建三角网；

以三角网中位于同一边两端的相干点的形变速率差分为待求解量，以相干点相位解缠后的相位为已知量，构建线性形变求解方程，并约束相邻时期的形变速率差，最终求解得到相干点的线性形变。下面对具体求解过程和原理进行解释说明：

解缠后的差分干涉相位采用下式表示：

式中，k为干涉图序列号，

为残余高程引起的相位值，是地物在雷达视线方向移动引起的相位值，

为基线残余误差相位，为大气延迟引起的相位，

为失相关噪声相位。(x,y)为像元的行列号。λ为波长，t_j为第j时刻的时间，v_j为第j时刻的速率，β_k为系数矩阵，V为速率矩阵，

为k幅图的形变量。

为卫星到地面点的距离，θ为入射角，B_⊥为垂直基线，Δh为高程。

对选择的间断相干点构建delaunry三角网，将三角网中的每个arc(边)的两端点(x,y)和(x',y')进行差分，差分结果如下式所示：

对应的矩阵形式为：

ΔΦ为M个已知的解缠相位，

为求解参数。假设存在8景SAR数据，组成14个干涉对，干涉对组合信息如下表所示。

表一8景SAR数据的干涉对组合信息

主影像	副影像	主影像	副影像	主影像	副影像	主影像	副影像
								1	2	3	4	4	7	6	8
1	3	3	5	5	6	7	8
								2	3	3	7	5	7
2	5	4	5	6	7

连续相干点对应的观测方程如式4所示

如果假定时间序列中形变为线性形变，即Δv₁＝Δv₂＝…Δv_N-1，此时只有两个待求参数，采用最小二乘即可求解。但是实践中这样的假设很难满足，不同时期的形变类型不一样，即使为线性形变，不同时期的形变速率也不一样。为此对相邻时期的形变速率差添加约束条件，如式5所示：

-0.5Δv_k-1+Δv_k-0.5Δv_k+1＝0 (5)

由于系数矩阵奇异，采用最小二乘可以获得无数解，因此不能采用最小二乘方法求解。需要对矩阵A进行SVD分解求取参数。

下面对间断相干点的分析处理进行解释说明，以图3为例，不是所有的ΔΦ都参与计算，而是待求点和参考点(参考点和待求点为三角网同一边两端的点，选取其中一点作为参考点，另一点为待求点)在时间序列中同时保持高相干性的ΔΦ。最终参与计算的观测量如图3中的Int1、Int2、Int4、Int5、Int7、Int8、Int9、Int10、Int11、Int13标签所示，对应的观测方程如式6：

相干点虽然在时间序列中保持高的相干性，但是不同时期的相干值不一样，对应的解缠的干涉相位也存在差异，而且可能存在解缠错误的情况。目前大都采用L₂范获得参数的最小二乘解，但是最小二乘方法在具有相同分布且独立的情况是最优的，对于其他情况，得到的参数值严重偏移模型值，同时，最小二乘方法容易受到数据中的粗差值的影响。基于这样的情况，本发明采用L₁范求解，如下：

假定A满秩，采用下式进行迭代求解：

其中P为权矩阵，采用下式确定

在求得每个弧段ΔV后，选定一个已知的稳定区域的点，使其作为一个参考点，可求得每个点的线性形变速率V。

在一实施例中，步骤108、步骤109对非线性形变的求解方法为：将解缠后的相位

减去线性形变相位

即可得到残余相位

在残余相位中，利用大气和非线性形变在时间域和空间域的特性，可以通过滤波的方法把他们进行分离。大气相位在时间上是一个白噪声，而非线性形变在空间上相关长度较小，在时间域具有低频特征。因此，把各差分干涉图的残余相位的均值

作为主影像的大气相位，

在时间域上的低通部分作为非线性形变的估计值，

在时间域上的高通和空间域上的平滑可以作为从影像上的大气相位，把从影像的大气相位加到主影像的大气相位就可以获得对应每个干涉图的大气相位。

下面以示例的方式对本发明进行说明，如图4所示，具体实施步骤如下：

步骤一：数据采用2012-2013年期间郑州市的17景TerraSAR-X数据，采用2013年01月11日的数据作为主影像，进行影像配准。

步骤二：选择时间基线小于77天和垂直基线小于120米的干涉组合，共组成26个干涉对，见图5。

步骤三：DInSAR差分干涉数据处理采用GAMMA DIFF软件包，采用SRTM DEM去除高程相位，多视视数采用21×21，多视后的像元空间分辨率40×40m。

步骤四：选择时序相干均值大于0.25，且相干值高于0.25的个数大于18的点参与计算。共选择704684个点，占总像元数的71.53％，点位分布如图6所示。

步骤五：对选择的相干点进行相位解缠。

步骤六：对解缠后的相位运用本发明提出的方法解算线性形变。

步骤七：对残余相位滤波分离大气和非线性形变。

步骤八：得到最终的形变速率，见图7。

由图5可见干涉对网络中没有孤立的干涉对，保证了结果的连通性。

由图6可见本发明提出的方法识别的点，不仅在市区而且在郊区的农田、森林和黄河沿岸有茂密的分布。而图8中传统SBAS方法识别的点主要分布在郑州市市区，在其他区域有稀疏的点位分布。

由图7可见本发明提出的方法的处理结果较之图8中的传统小基线方法更能反应地表的真实形变，尤其是在农田，森林等低相干区域，如图下方红圈处。

本发明通过选择在长时间序列中的间断高相干点，运用小基线集InSAR技术，能有效获取低相干区域，如森林，农田等植被覆盖区域的地表形变，提高了地表形变监测的精度。

本发明还提供了一种基于间断相干的InSAR地表形变监测系统，如图9所示，该系统包括：

地表SAR影像获取模块901，用于获取被监测区域的地表SAR影像，SAR影像包括时间序列上的多个卫星雷达影像；

影像配准模块902，用于从SAR影像中选取主影像，并根据主影像对SAR影像进行配准；

干涉对提取模块903，用于采用小基线集法提取SAR影像中的干涉对；

差分干涉处理模块904，用于对干涉对进行差分干涉处理；

相干点提取模块905，用于提取所述干涉对中的相干点，所述相干点包括连续相干点和间断相干点，所述连续相干点为在时间序列上相干值的均值大于第一预设值的相干点；所述间断相干点为在时间序列上具有多个相干值大于第二预设值且大于所述第二预设值的相干值的数量大于设定值的相干点；

线性形变确定模块907，用于根据相干点相位解缠后的相位确定相干点的线性形变；

残余相位确定模块908，用于从相干点相位解缠后的相位中减去相干点的线性形变相位，得到残余相位；

非线性形变确定模块909，用于对残余相位进行滤波处理，得到相干点的非线性形变；

形变确定模块910，用于综合相干点的线性形变和非线性形变，确定被监测区域地表的形变或形变速率。

在一实施例中，差分干涉处理模块904，具体包括：

差分干涉处理单元，用于采用轨道数据和地形数据对干涉对进行二轨差分干涉处理。

在一实施例中，线性形变确定模块907，具体包括：

三角网构建单元，用于对相干点构建三角网；

线性形变求解单元，用于以三角网中位于同一边两端的相干点的形变速率差为待求解量，以相干点相位解缠后的相位为已知量，构建线性形变求解方程，并约束相邻时期的形变速率差，最终求解得到相干点的线性形变。

在一实施例中，相位解缠模块，具体包括：

相位解缠单元906，用于采用最小费用流方法对相干点进行相位解缠。

本发明提供的基于间断相干的InSAR地表形变监测方法及系统在选取被监测区域SAR影像上的相干点时，既选取了在时间序列上连续保持高相干性的连续相干点，还选取了在时间序列上间断保持高相干性的间断相干点；根据这些相干点的相位变化，确定被监测区域地表的形变或形变速率。由于本发明在提取相干点时，考虑了时序间断相干点，使得本发明能够对时序上点位间断相干的区域的形变进行有效监测。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种基于间断相干的InSAR地表形变监测方法，其特征在于，包括：

获取被监测区域的地表SAR影像，所述SAR影像包括时间序列上的多个卫星雷达影像；

从所述SAR影像中选取主影像，并根据所述主影像对所述SAR影像进行配准；

采用小基线集法提取所述SAR影像中的干涉对；

对所述干涉对进行差分干涉处理；

提取所述干涉对中的相干点，所述相干点包括连续相干点和间断相干点，所述连续相干点为在时间序列上相干值的均值大于第一预设值的相干点；所述间断相干点为在时间序列上具有多个相干值大于第二预设值且大于所述第二预设值的相干值的数量大于设定值的相干点；

对所述相干点进行相位解缠；

根据所述相干点相位解缠后的相位确定所述相干点的线性形变；

从所述相干点相位解缠后的相位中减去所述相干点的线性形变相位，得到残余相位；

对所述残余相位进行滤波处理，得到所述相干点的非线性形变；

综合所述相干点的线性形变和非线性形变，确定被监测区域地表的形变或形变速率。

2.根据权利要求1所述的基于间断相干的InSAR地表形变监测方法，其特征在于，所述对所述干涉对进行差分干涉处理，具体包括：

采用轨道数据和地形数据对所述干涉对进行二轨差分干涉处理。

3.根据权利要求1所述的基于间断相干的InSAR地表形变监测方法，其特征在于，所述根据所述相干点相位解缠后的相位确定所述相干点的线性形变，具体包括：

对所述相干点构建三角网；

以三角网中位于同一边两端的相干点的形变速率差分为待求解量，以所述相干点相位解缠后的相位为已知量，构建线性形变求解方程，并约束相邻时期的形变速率差，最后求解得到所述相干点的线性形变。

4.根据权利要求3所述的基于间断相干的InSAR地表形变监测方法，其特征在于，所述约束相邻时期的形变速率，具体包括：

将相邻时期形变速率的约束确定为-0.5Δv_k-1+Δv_k-0.5Δv_k+1＝0，其中，Δv_k为k时期的形变速率差。

5.根据权利要求1所述的基于间断相干的InSAR地表形变监测方法，其特征在于，所述对所述相干点进行相位解缠，具体包括：

采用最小费用流方法对所述相干点进行相位解缠。

6.一种基于间断相干的InSAR地表形变监测系统，其特征在于，包括：

地表SAR影像获取模块，用于获取被监测区域的地表SAR影像，所述SAR影像包括时间序列上的多个卫星雷达影像；

影像配准模块，用于从所述SAR影像中选取主影像，并根据所述主影像对所述SAR影像进行配准；

干涉对提取模块，用于采用小基线集法提取所述SAR影像中的干涉对；

差分干涉处理模块，用于对所述干涉对进行差分干涉处理；

相干点提取模块，用于提取所述干涉对中的相干点，所述相干点包括连续相干点和间断相干点，所述连续相干点为在时间序列上相干值的均值大于第一预设值的相干点；所述间断相干点为在时间序列上具有多个相干值大于第二预设值且大于所述第二预设值的相干值的数量大于设定值的相干点；

相位解缠模块，用于对所述相干点进行相位解缠；

线性形变确定模块，用于根据所述相干点相位解缠后的相位确定所述相干点的线性形变；

残余相位确定模块，用于从所述相干点相位解缠后的相位中减去所述相干点的线性形变相位，得到残余相位；

非线性形变确定模块，用于对所述残余相位进行滤波处理，得到所述相干点的非线性形变；

形变确定模块，用于综合所述相干点的线性形变和非线性形变，确定被监测区域地表的形变或形变速率。

7.根据权利要求6所述的基于间断相干的InSAR地表形变监测系统，其特征在于，所述差分干涉处理模块，具体包括：

差分干涉处理单元，用于采用轨道数据和地形数据对所述干涉对进行二轨差分干涉处理。

8.根据权利要求6所述的基于间断相干的InSAR地表形变监测系统，其特征在于，所述线性形变确定模块，具体包括：

三角网构建单元，用于对所述相干点构建三角网；

线性形变求解单元，用于以三角网中位于同一边两端的相干点的形变速率差分为待求解量，以所述相干点相位解缠后的相位为已知量，构建线性形变求解方程，并约束相邻时期的形变速率差，最终求解得到所述相干点的线性形变。

9.根据权利要求8所述的基于间断相干的InSAR地表形变监测系统，其特征在于，相邻时期形变速率的约束为-0.5Δv_k-1+Δv_k-0.5Δv_k+1＝0，其中，Δv_k为k时期的形变速率差。

10.根据权利要求6所述的基于间断相干的InSAR地表形变监测系统，其特征在于，所述相位解缠模块，具体包括：

相位解缠单元，用于采用最小费用流方法对所述相干点进行相位解缠。

Images