[发明专利]一种结合路径跟踪的多基线InSAR高程重建方法

权利要求书

1.一种结合路径跟踪的多基线InSAR高程重建方法，其特征在于，包括如下步骤；

1)获取干涉相位质量图

获取干涉合成孔径雷达SAR图像中最长基线对应的二维干涉图，即缠绕相位图，对二维干涉图进行预处理，获得干涉图像中反映干涉像元质量的像元的相位质量值，即干涉相位质量图，指导相位展开路径，干涉相位质量图由干涉图直接获得，相位质量值ρ_m,n的计算公式为：

公式(1)中，是以像素点(m,n)为中心的局部窗口内缠绕相位值，k×k为局部窗口的尺寸；

2)从步骤1)得到的干涉相位质量图中，选择相位质量值最大的非边界点，使用最大似然算法对该非边界点进行最大似然估计，得到估计高程值，作为路径跟踪策略的起始像素点，最大似然率和联合最大似然率F_ML(φ|h)的计算公式如下：

公式(2)和公式(3)中，κ为高度灵敏度，h为对应的地表高程值，代表第i幅干涉图的缠绕相位γ_i为第i幅干涉图的相干系数；

3)在步骤2)中获得已得到估计高程值的起始像素点，将起始像素点4邻域中未解缠的点，按相位质量从大到小排列作为邻接列；

4)取出邻接列中质量最好的点，在以该点为中心3*3的窗口内，计算该窗口内已解缠像素点间的被估计点的高程值若最大高程差h_diff小于阈值h_th，则认为该点保持连续；否则认为该点为突变点，并将该点4邻域中未解缠的点加入邻接列；

5)判断邻接列是否为空，若非空，则继续执行步骤4)；若为空，则算法结束，得到重建的高程图。

2.根据权利要求1所述的一种结合路径跟踪的多基线InSAR高程重建方法，其特征在于，步骤1)中，利用生成相位质量图，将单基线相位解缠算法中的路径跟踪策略引入到多基线相位解缠算法中。

3.根据权利要求1所述的一种结合路径跟踪的多基线InSAR高程重建方法，其特征在于，步骤4)中，所述的被估计点的高程值计算公式为：

公式(4)和公式(5)中，β是归一化参数；y∈N_x表示y为像素点x邻域的像素点集合，一般取8连通域；表示以x点为中心的3*3矩阵中已解缠像素点的平均高程值，是一个标量；表示待解缠像素点的可能高程值，是一个向量；阈值h_th设置为最长基线对应的高度模糊数，公式(4)为迭代递推的估计运算。