[发明专利]基于平面近似模型的多基线相位解缠绕方法

权利要求书

1.一种基于平面近似模型的多基线相位解缠绕方法，其特征在于，所述方法包括如下步骤：

步骤1，从多基线InSAR系统中获取两个基线分别对应的干涉相位图，记为第一干涉相位图和第二干涉相位图，两幅干涉相位图的大小相同，且为m×n，其中，m、n分别为大于零的正整数；

步骤2，根据第一干涉相位图中每个像素的缠绕相位，确定第一干涉相位图中水平方向相邻像素间的水平缠绕相位梯度，记为水平缠绕第一相位梯度矩阵，其大小为m×(n-1)，并确定第一干涉相位图中垂直方向相邻像素间的垂直缠绕相位梯度，记为垂直缠绕第一相位梯度矩阵，其大小为(m-1)×n；

根据第二干涉相位图中每个像素的缠绕相位，确定第二干涉相位图中水平方向相邻像素间的水平缠绕相位梯度，记为水平缠绕第二相位梯度矩阵，其大小为m×(n-1)，并确定第二干涉相位图中垂直方向相邻像素间的垂直缠绕相位梯度，记为垂直缠绕第二相位梯度矩阵，其大小为(m-1)×n；

步骤3，设定预设窗口大小，根据预设窗口大小选定给定像素位置对应的估计窗口；

步骤4，给定像素位置(i，j)，1≤i≤m，1≤j≤n-1；

根据给定像素位置(i，j)，以及给定像素位置(i，j)对应的估计窗口，估计得到第一干涉相位图的水平方向模糊数梯度矩阵，其大小为m×(n-1)；以及估计得到第二干涉相位图的水平方向模糊数梯度矩阵，其大小为m×(n-1)；

步骤5，给定像素位置(i，j)，1≤i≤m-1，1≤j≤n；

根据给定像素位置(i，j)，以及给定像素位置(i，j)对应的估计窗口，估计得到第一干涉相位图的垂直方向模糊数梯度矩阵，其大小为(m-1)×n；以及估计得到第二干涉相位图的垂直方向模糊数梯度矩阵，其大小为(m-1)×n；

步骤6，将第一干涉相位图在水平方向的模糊数梯度矩阵中的元素以及干涉相位图1在垂直方向的模糊数梯度矩阵中的元素带入设定的优化模型，得到第一干涉相位图的模糊数矩阵1；

将第二干涉相位图在水平方向的模糊数梯度矩阵的元素以及第二干涉相位图在垂直方向的模糊数梯度矩阵的元素带入所述设定的优化模型，得到第二干涉相位图的模糊数矩阵2；

步骤7，给定像素位置(i，j)，1≤i≤m，1≤j≤n；

根据第一干涉相位图在像素位置(i，j)处对应的缠绕相位和第一干涉相位图在像素位置(i，j)处的模糊数，计算得到第一干涉相位图在像素位置(i，j)处的绝对相位；

根据第二干涉相位图在像素位置(i，j)处对应的缠绕相位和第二干涉相位图在像素位置(i，j)处的模糊数，计算得到第二干涉相位图在像素位置(i，j)处的绝对相位；

步骤8，令像素位置(i，j)遍历集合{(i，j)|1≤i≤m，1≤j≤n}中的所有元素，重复执行步骤7，从而得到第一干涉相位图的绝对相位矩阵，以及第二干涉相位图的绝对相位矩阵。

2.根据权利要求1所述的一种基于平面近似模型的多基线相位解缠绕方法，其特征在于，步骤2具体包括如下子步骤：

(2a)获取第一干涉相位图中像素位置(i，j)处的缠绕相位获取第二干涉相位图中像素位置(i，j)处的缠绕相位像素位置(i，j)处的像素值为像素位置(i，j)处的缠绕相位；

(2b)记第一干涉相位图的水平缠绕第一相位梯度矩阵为中第i行第j列元素

记第二干涉相位图的水平缠绕第二相位梯度矩阵为中第i行第j列元素

(2c)记第一干涉相位图的垂直缠绕第一相位梯度矩阵为中第i行第j列元素

记第二干涉相位图的垂直缠绕第二相位梯度矩阵为中第i行第j列元素