[發(fā)明專利]閉合魯棒雙基線InSAR相位解纏方法、系統(tǒng)及可讀存儲介質

有效

• 申請?zhí)枺?/td>: 201911147113.3
• 申請日: 2019-11-21
• 公開（公告）號: CN110988876B
• 公開（公告）日: 2023-01-03
• 發(fā)明（設計）人: 袁志輝;尹英豪;邢學敏;陳立福;谷翼策
• 申請（專利權）人: 長沙理工大學
• 主分類號: G01S13/90
• 分類號: G01S13/90;G06V10/762
• 代理公司: 長沙市融智專利事務所(普通合伙) 43114
• 代理人: 龔燕妮

1.一種閉合魯棒雙基線InSAR相位解纏方法，其特征在于，包括：

根據雙基線InSAR系統(tǒng)參數求出各干涉圖對應的高度模糊數，并將它們分解為公因子和整數相乘的形式，所述整數被稱為各干涉圖對應的分解整數；

根據各干涉圖對應的分解整數及由雙基線InSAR系統(tǒng)測得的帶噪聲纏繞干涉相位矢量求出各像素對應的像素截距；

根據各像素對應的像素截距對所有像素進行聚類，從而獲得像素的總類數、各像素所屬的類和每個類對應的類截距；

利用閉合求解公式計算所述每個類對應的類模糊矢量，進而得到各像素對應的像素模糊矢量；

采用類相位濾波法，將各像素對應的帶噪聲纏繞干涉相位矢量點沿投影直線投影到該像素對應的類線段上進行濾波，獲得濾波后的纏繞干涉相位矢量點；所述投影直線是過帶噪聲纏繞干涉相位矢量點的一條直線，其方向向量由各干涉圖對應的相干系數決定；所述類線段是指同類像素在以兩纏繞干涉相位為坐標軸所構成的平面中對應的一條線段；根據各像素對應的像素模糊矢量和濾波后的纏繞干涉相位得到解纏后的絕對干涉相位，并依據解纏后的絕對干涉相位和雙基線InSAR系統(tǒng)的其它參數重建地形高程。

2.根據權利要求1所述的方法，其特征在于，當各干涉圖對應的高度模糊數均為整數時，公因子為所有高度模糊數的最大公約數；

當各干涉圖對應的高度模糊數不全為整數時，首先將所有高度模糊數保留相同的n位小數點位數，接著將所有高度模糊數與10ⁿ相乘，然后求取最大公約數，最后將所得的最大公約數與10ⁿ相除，從而得到公因子；各干涉圖對應的分解整數為各干涉圖對應的高度模糊數除以公因子后所得的整數。

3.根據權利要求1所述的方法，其特征在于，所述類相位濾波法，將各像素對應的帶噪聲纏繞干涉相位矢量點沿投影直線投影到該像素對應的類線段上進行濾波，獲得濾波后的纏繞干涉相位矢量點，具體過程如下：

把平面上通過像素s對應的帶噪聲纏繞干涉相位矢量點且方向向量為[-1/|γ₁|，1/|γ₂|]的投影直線與該像素對應的類線段的交點作為濾波后的纏繞干涉相位矢量點類線段所在的直線經過一個與類截距intercept_c(l)有關的特定點[0，-2π·intercept_c(l)]，且其方向向量為[1/Γ₁，1/Γ₂]；γ₁和γ₂分別為兩干涉圖對應的相干系數；

從而獲得以下兩條點向式直線方程：

經過變換后得：

其中，為自變量，為因變量，分別對應平面上的橫坐標軸和縱坐標軸，為像素s對應的帶噪聲纏繞干涉相位矢量點，和分別為像素s在兩干涉圖上的帶噪聲纏繞干涉相位；intercept_c(l)為像素s對應的類截距，像素s假設屬于第l類，其中l(wèi)表示各像素所屬 的類；

求解上述方程組便獲得濾波后的纏繞干涉相位矢量點其表達式如下：

其中，Γ₁和Γ₂分別表示兩干涉圖對應的分解整數。

4.根據權利要求3所述的方法，其特征在于，當|γ₁|/|γ₂|＝Γ₂/Γ₁時，采用垂直投影；當|γ₁|/|γ₂|＝0時，采用水平投影；當|γ₁|/|γ₂|＝∞時，采用豎直投影。