[发明专利]一种InSAR地形复杂区域干涉相位解缠方法

权利要求书

1.一种InSAR地形复杂区域干涉相位解缠方法，其特征在于，包括：

步骤1，获取目标区域的SAR影像数据；

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

步骤3，计算整周数；

步骤4，获得解缠相位结果。

2.根据权利要求1所述的方法，其中，所述步骤1，获取SAR影像数据，具体包括：

按照时间顺序获取N景SAR影像，其中，第1景SAR影像为主影像，其余N-1景SAR影像为副影像，对N-1对主副影像进行配准。

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

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

步骤2-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影像拍摄时发生的地表高程形变；λ表示波长；

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

步骤2-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个像元之间的整周数差；

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

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

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

4.根据权利要求1所述的方法，其中，所述步骤3，计算整周数k，具体包括：

步骤3-1，确定相干系数图：令第i个干涉像对的主副影像，即第1景和第i+1景SAR影像分别为S₁和S₂，建立第i个干涉像对的相干系数图，表示为：其中，E()表示期望值；*表示共轭；表示S2的共轭复数；

步骤3-2，相干系数图分割：在相干系数图中标记残差点并确定所述残差点在所述相干系数图中的分布，确定所述残差点在所述相干系数图上对应的相关系数值，以所有残差点中最大的相关系数值作为分割数据的阈值，将所述相干系数图和对应的缠绕相位图分割成高质量区域和低质量区域；

步骤3-3，高质量区域的相位解缠；

步骤3-4，低质量区域的相位解缠。

5.根据权利要求4所述的方法，其中，所述步骤3-3，高质量区域的相位解缠，具体包括：

步骤3-3-1，获取缠绕相位图的高质量区域中任意一点的GPS信息；

步骤3-3-2，利用所述GPS信息反演得到所述点对应的整周数；

步骤3-3-3，对高质量区域进行整周数积分，得到高质量区域每个像元对应的整周数。