[发明专利]一种基于缠绕相位的GBInSAR大气校正方法及系统

权利要求书

1.一种基于缠绕相位的GBInSAR大气校正方法，其特征在于，包括：

获取不同时刻对同一地区进行探测得到的两幅单视复数影像以及由两幅所述单视复数影像差分干涉形成的干涉图；

根据所述干涉图中各个像素点的相干性，选取所述干涉图中的高相干点，具体包括：

采用相干性计算公式：

其中，表示相干性，m和n分别表示滑动窗口的两条边的大小，M和S分别表示两幅单视复数影像，i和j分别表示干涉图块中像素点的两个坐标，M(i,j)表示单视复数影像M中坐标为(i,j)的像素点，S(i,j)表示单视复数影像S中坐标为(i,j)的像素点，S^*为S的矩阵转置，S^*(i,j)表示S^*中坐标为(i,j)的像素点；

选取相干性超过预设阈值的像素点作为高相干点；

获取处于同一个方位向的相邻的所述高相干点之间的实测的相位差值和距离差值，得到相邻相位差和相邻距离差；

根据所述相邻相位差和相邻距离差利用最小二乘回归法计算大气延迟系数；

根据所述大气延迟系数计算大气延迟相位；

根据所述大气延迟相位与所述相邻相位差的差值大小剔除超出预设阈值的所述差值大小所对应的相邻相位差和相邻距离差，得到去粗差后的相邻相位差和相邻距离差；

根据所述去粗差后的相邻相位差和相邻距离差利用最小二乘回归法更新大气延迟系数，得到去粗差后的大气延迟系数；

利用所述去粗差后的大气延迟系数计算粗差改正后的大气延迟相位。

2.根据权利要求1所述的基于缠绕相位的GBInSAR大气校正方法，其特征在于，在所述获取不同时刻对同一地区进行探测得到的两幅单视复数影像以及由两幅所述单视复数影像差分干涉形成的干涉图之后，在所述根据所述干涉图中各个像素点的相干性，选取所述干涉图中的高相干点之前，还包括：

对所述干涉图进行Goldstein滤波处理。

3.根据权利要求1所述的基于缠绕相位的GBInSAR大气校正方法，其特征在于，所述根据所述干涉图中各个像素点的相干性，选取所述干涉图中的高相干点，具体包括：

计算所述干涉图中每个像素的相干性，所述相干性表示通过差分干涉构成所述干涉图中某个像素点的两个像素点之间回波相位的相似程度；

选取相干性超过预设阈值的像素点作为高相干点。

4.根据权利要求2所述的基于缠绕相位的GBInSAR大气校正方法，其特征在于，所述对所述干涉图进行Goldstein滤波处理，具体包括：

利用滑动窗口将所述干涉图划分成多个干涉图块，且任意相邻的两个所述干涉图块的重叠率大于或等于预设比率；

对每个所述干涉图块进行傅里叶变换，得到每个所述干涉图块的频谱图；

对每个所述频谱图进行平滑处理，得到平滑处理的幅值；

利用所述平滑处理的幅值对各个所述频谱图的频谱进行处理，得到处理后的频谱图；

对所述处理后的频谱图进行傅里叶逆变换得到滤波后的干涉图。

5.根据权利要求1所述的基于缠绕相位的GBInSAR大气校正方法，其特征在于，在所述获取处于同一个方位向的相邻的所述高相干点之间的实测的相位差值和距离差值，得到相邻相位差和相邻距离差之后，在所述根据所述相邻相位差和相邻距离差利用最小二乘回归法计算大气延迟系数之前，还包括：

将高于预设距离差阈值的相邻距离差以及对应的相邻相位差数据剔除。