[發(fā)明專利]一種基于纏繞相位的GBInSAR大氣校正方法及系統(tǒng)

有效

• 申請?zhí)枺?/td>: 201910982029.7
• 申請日: 2019-10-16
• 公開（公告）號: CN110658521B
• 公開（公告）日: 2021-07-27
• 發(fā)明（設計）人: 楊紅磊;劉杰;彭軍還
• 申請（專利權）人: 中國地質大學（北京）
• 主分類號: G01S13/90
• 分類號: G01S13/90;G01S7/40
• 代理公司: 北京高沃律師事務所 11569
• 代理人: 程江濤

1.一種基于纏繞相位的GBInSAR大氣校正方法，其特征在于，包括：

獲取不同時刻對同一地區(qū)進行探測得到的兩幅單視復數(shù)影像以及由兩幅所述單視復數(shù)影像差分干涉形成的干涉圖；

根據(jù)所述干涉圖中各個像素點的相干性，選取所述干涉圖中的高相干點，具體包括：

采用相干性計算公式：

其中，表示相干性，m和n分別表示滑動窗口的兩條邊的大小，M和S分別表示兩幅單視復數(shù)影像，i和j分別表示干涉圖塊中像素點的兩個坐標，M(i,j)表示單視復數(shù)影像M中坐標為(i,j)的像素點，S(i,j)表示單視復數(shù)影像S中坐標為(i,j)的像素點，S^*為S的矩陣轉置，S^*(i,j)表示S^*中坐標為(i,j)的像素點；

選取相干性超過預設閾值的像素點作為高相干點；

獲取處于同一個方位向的相鄰的所述高相干點之間的實測的相位差值和距離差值，得到相鄰相位差和相鄰距離差；

根據(jù)所述相鄰相位差和相鄰距離差利用最小二乘回歸法計算大氣延遲系數(shù)；

根據(jù)所述大氣延遲系數(shù)計算大氣延遲相位；

根據(jù)所述大氣延遲相位與所述相鄰相位差的差值大小剔除超出預設閾值的所述差值大小所對應的相鄰相位差和相鄰距離差，得到去粗差后的相鄰相位差和相鄰距離差；

根據(jù)所述去粗差后的相鄰相位差和相鄰距離差利用最小二乘回歸法更新大氣延遲系數(shù)，得到去粗差后的大氣延遲系數(shù)；

利用所述去粗差后的大氣延遲系數(shù)計算粗差改正后的大氣延遲相位。

2.根據(jù)權利要求1所述的基于纏繞相位的GBInSAR大氣校正方法，其特征在于，在所述獲取不同時刻對同一地區(qū)進行探測得到的兩幅單視復數(shù)影像以及由兩幅所述單視復數(shù)影像差分干涉形成的干涉圖之后，在所述根據(jù)所述干涉圖中各個像素點的相干性，選取所述干涉圖中的高相干點之前，還包括：

對所述干涉圖進行Goldstein濾波處理。

3.根據(jù)權利要求1所述的基于纏繞相位的GBInSAR大氣校正方法，其特征在于，所述根據(jù)所述干涉圖中各個像素點的相干性，選取所述干涉圖中的高相干點，具體包括：

計算所述干涉圖中每個像素的相干性，所述相干性表示通過差分干涉構成所述干涉圖中某個像素點的兩個像素點之間回波相位的相似程度；

選取相干性超過預設閾值的像素點作為高相干點。

4.根據(jù)權利要求2所述的基于纏繞相位的GBInSAR大氣校正方法，其特征在于，所述對所述干涉圖進行Goldstein濾波處理，具體包括：

利用滑動窗口將所述干涉圖劃分成多個干涉圖塊，且任意相鄰的兩個所述干涉圖塊的重疊率大于或等于預設比率；

對每個所述干涉圖塊進行傅里葉變換，得到每個所述干涉圖塊的頻譜圖；

對每個所述頻譜圖進行平滑處理，得到平滑處理的幅值；

利用所述平滑處理的幅值對各個所述頻譜圖的頻譜進行處理，得到處理后的頻譜圖；

對所述處理后的頻譜圖進行傅里葉逆變換得到濾波后的干涉圖。

5.根據(jù)權利要求1所述的基于纏繞相位的GBInSAR大氣校正方法，其特征在于，在所述獲取處于同一個方位向的相鄰的所述高相干點之間的實測的相位差值和距離差值，得到相鄰相位差和相鄰距離差之后，在所述根據(jù)所述相鄰相位差和相鄰距離差利用最小二乘回歸法計算大氣延遲系數(shù)之前，還包括：

將高于預設距離差閾值的相鄰距離差以及對應的相鄰相位差數(shù)據(jù)剔除。