[發明專利]一種基于差頻InSAR的緩變高程反演方法

在審

• 申請號: 202110995142.6
• 申請日: 2021-08-27
• 公開（公告）號: CN113835089A
• 公開（公告）日: 2021-12-24
• 發明（設計）人: 劉波;孔維亞;眭曉虹;袁仕耿;李洋;王沫
• 申請（專利權）人: 中國空間技術研究院;航天東方紅衛星有限公司
• 主分類號: G01S13/90
• 分類號: G01S13/90
• 代理公司: 北京善任知識產權代理有限公司 11650
• 代理人: 張振偉

1.一種基于差頻InSAR的緩變高程反演方法，該方法包括：

S1：對成像后的第一天線的單視復圖像和第二天線的單視復圖像進行配準；

S2：對配準后的第一天線的單視復圖像和第二天線的單視復圖像共軛相乘，取復數據的相位部分即得到干涉相位，所述干涉相位中包含平地相位和高程相位，暫時將平地相位從所述干涉相位中去除，其中計算平地相位所采用的載頻與步驟S3中同頻等效后的頻率一致；

S3：去除平地相位后，對剩余相位補償差頻相位項，將所述剩余相位等效轉化為差頻InSAR中的第一天線的載頻對應的同頻高程相位，第二天線的載頻與步驟S2中計算平地相位所采用的載頻一致；

S4：將所述同頻高程相位中的同頻平地相位加回至高程相位之中，基于同頻InSAR高程反演數學式計算獲得高程。

2.根據權利要求1所述的基于差頻InSAR的緩變高程反演方法，其中，所述步驟S3還包括，當所述同頻高程相位出現相位纏繞時，對所述同頻高程相位進行相位解纏。

3.根據權利要求1所述的基于差頻InSAR的緩變高程反演方法，其中，所述步驟S1還包括，

基于差頻InSAR系統參數計算兩幅單視復圖像在方位向和距離向上的像素級差異并完成粗配準，然后采用相干系數法對兩幅單視復圖像進行亞像素級配準。

4.根據權利要求3所述的基于差頻InSAR的緩變高程反演方法，其中，所述差頻InSAR系統參數為基線長度、平臺高度、平臺速度中的一種或多種。

5.根據權利要求1所述的基于差頻InSAR的緩變高程反演方法，其中，所述步驟S2還包括：根據公式計算平地相位，從而將平地相位從干涉相位中暫時去除，其中R₁和R₂分別為第一天線和第二天線到參考平面某點的距離，λ₁和λ₂分別為第一天線和第二天線的發射波長。

6.根據權利要求1所述的基于差頻InSAR的緩變高程反演方法，其中，在步驟S3中，對剩余相位補償的差頻相位項為

7.根據權利要求1所述的基于差頻InSAR的緩變高程反演方法，其中，步驟S4包括：

基于公式計算獲得高程，

其中，H為平臺飛行高度；r₁為目標至第一天線的距離；B為干涉SAR基線長度；α為基線傾角；λ₂為第二天線的發射波長；△φ′為補償差頻相位項后的高程相位。

8.根據權利要求1所述的基于差頻InSAR的緩變高程反演方法，其中，所述步驟S1還包括：在對成像后的第一天線的單視復圖像和第二天線的單視復圖像進行配準之前，第一天線發射中心頻率為f₀的信號，第二天線發射中心頻率為f₀+Δf的信號，第一天線接收回波并進行SAR成像得到所述第一天線的單視復圖像，第二天線接收回波并進行SAR成像得到所述第二天線的單視復圖像。