[发明专利]一种基于时序的InSAR技术监测地表形变的方法

权利要求书

1.一种基于时序的InSAR技术监测地表形变的方法，其特征在于，包括：

步骤1，SAR数据选取：通过SAR传感器获取SAR数据，形成立体像对；

步骤2，立体像对配准：将立体像对进行精确配准，根据配准多项式将立体像对的配准误差控制在1/8个像元以下；

步骤3，噪声滤波处理：对立体像对进行带通滤波处理；

步骤4，干涉纹图生成：对立体像对进行共轭相乘操作，生成干涉纹图；

步骤5，进行去平地效应操作，得到去平干涉相位图；

步骤6，进行相位解缠，计算影像的真实相位值，获得地表高程形变结果；

步骤7，进行高程及高程坐标系转换，得到国际标准地理参考系中的地表高程形变结果。

2.根据权利要求1所述的方法，其中，所述步骤1，SAR数据选取：通过SAR传感器获取SAR数据，形成立体像对，具体包括：

按照时间顺序获取N景SAR影像，其中，第1景SAR影像为主影像，其余N-1景SAR影像为副影像，形成立体像对。

3.根据权利要求1所述的方法，其中，所述步骤6，进行相位解缠，计算影像的真实相位值，获得地表高程形变结果，具体包括：

步骤6-1，求取相邻像元之间的整周数差；

步骤6-2，计算整周数；

步骤6-3，获得解缠相位结果；

步骤6-4，获得地表高程形变结果。

4.根据权利要求3所述的方法，其中，所述步骤6-1，求取相邻像元之间的整周数差，具体包括：

步骤6-1-1，第1景SAR主影像和其余N-1景SAR副影像顺序组成N-1个SAR干涉像对；

步骤6-1-2，由于第1景和第2景SAR影像拍摄时间相距很近，期间地表形变很小予以忽略，从第3景开始的影像至第1景SAR影像期间地表发生的形变不能忽略，针对每一干涉像对，根据时序建立所述N-1景副影像至第1景主影像的地表高程形变与第s个像元对应的绝对相位以及整周数的第一关系模型：

其中，表示第i个干涉像对中第s个像元对应的绝对相位；k_i(s)(i＝1,2,…,N-1)表示第i个干涉像对中第s个像元对应的整周数；B_i(i＝1,2,…,N-1)表示第i个干涉像对的垂直基线长；Δh^m(s)(m＝2,3,…,N-1)表示第m+1景SAR影像拍摄时相对于第1景SAR影像拍摄时发生的地表高程形变；λ表示波长；

步骤6-1-3，针对每一干涉像对，根据时序建立所述N-1景副影像至第1景主影像的地表高程形变与第s-1个像元对应的绝对相位以及整周数的第二关系模型，所述第s-1个像元表示第s个像元的相邻像元：

步骤6-1-4，由于第s-1个像元与第s个像元是相邻像元，令两者的高程差Δh(s)-Δh(s-1)＝0，结合所述第一关系模型和所述第二关系模型，针对每一干涉像对，根据时序建立第s个像元与第s-1个像元的绝对相位差与整周数差的第三关系模型：

其中，表示第i个干涉像对中第s个像元与第s-1个像元之间的绝对相位差，Δk_i(s,s-1)＝k_i(s)-k_i(s-1)(i＝1,2,…,N-1)，表示第i个干涉像对中第s个像元与第s-1个像元之间的整周数差；

步骤6-1-5，为求取Δk_i(s,s-1)，将第三关系模型转换为最优解模型：

其中，i,j＝1,2,…,N-1且互不相同；

步骤6-1-6，求取最优解模型的最优解，得到相邻像元之间的整周数差Δk_i(s,s-1)。