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

摘要

本发明公开了一种结合路径跟踪的多基线InSAR高程重建方法，1）获取干涉相位质量图；2）选择相位质量值最大的非边界点，用最大似然算法对非边界点进行最大似然估计得估计高程值；3）将起始像素点4邻域中未解缠的点，按相位质量从大到小排列作为邻接列；4）取邻接列中质量最好的点，以该点为中心3*3的窗口内计算该窗口内已解缠像素点间的被估计点的高程值，若最大高程差小于阈值，则该点保持连续；否则该点为突变点，并将该点4邻域中未解缠的点加入邻接列；5）判断邻接列是否为空，若非空，则继续执行步骤4）；若为空，则算法结束，得到重建的高程图。该方法有效消除了ML算法中的毛刺现象，提高了多基线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)；若为空，则算法结束，得到重建的高程图。