[發明專利]一種InSAR干涉相位真值確定及差分干涉測量方法

有效

• 申請號: 201810600536.5
• 申請日: 2018-06-12
• 公開（公告）號: CN108919264B
• 公開（公告）日: 2020-12-01
• 發明（設計）人: 趙超英;康亞;王寶行;劉曉杰
• 申請（專利權）人: 長安大學
• 主分類號: G01S13/90
• 分類號: G01S13/90
• 代理公司: 西安恒泰知識產權代理事務所 61216
• 代理人: 張明

本發明公開了一種InSAR干涉相位真值確定及差分干涉測量方法，用于確定InSAR差分干涉圖中每一個像素點的觀測值的解纏誤差，以獲得InSAR差分干涉圖干涉相位真值，本發明通過擬準檢定法將解纏誤差與其他相對較小誤差分離，并將解纏誤差估計、去除；本發明提供的InSAR干涉相位真值確定及差分干涉測量方法，可以解決由于干涉相位真值估計的不準確的問題，從而提高SBAS?InSAR的監測精度。

技術領域

本發明涉及InSAR數據處理領域，具體涉及一種InSAR干涉相位真值確定及差分干涉測量方法。

背景技術

合成孔徑雷達差分干涉測量技術(InSAR)是一種可以有效獲取地表形變的技術。目前InSAR技術憑借其全天候、無接觸、高精度等特點，已經廣泛地應用于地震、火山、地面沉降、滑坡等多類地物的監測。為使用SBAS-InSAR技術達到較高的監測精度，需要對InSAR觀測值各類誤差進行估計并減弱。其誤差可以分為：大氣延遲，時間去相干，DEM誤差，解纏誤差等。其中除了解纏誤差外，其余幾種誤差可以通過時序分析或者外部數據有效地減弱。

目前針對解纏誤差的識別與消除主要通過構建最小閉合環的方式進行，理論上可以采用相位自動檢測與校正的方法進行解纏誤差的改正，但是在實際情況下，解纏干涉圖中存在很多局部誤差，這些誤差加上噪聲等因素造成的誤差，校正非常復雜也非常耗時，因此目前通用的方法為人機交互手動校正，而隨著SAR數據量的不斷增加，人工識別改正解纏誤差將耗費大量的人力與物力，因此現有技術采用SBAS-InSAR技術時，由于干涉相位真值估計的不準確，導致SBAS-InSAR的監測精度低。

發明內容

本發明的目的在于提供一種修正InSAR解纏誤差方法，用以解決現有技術中由于干涉相位真值估計的不準確，導致SBAS-InSAR的監測精度低的問題。

為了實現上述任務，本發明采用以下技術方案：

一種InSAR干涉相位真值確定方法，用于獲得InSAR差分干涉圖中一個像素點觀測值的真值，所述的方法包括：

步驟1、通過建立觀測方程求取當前像素點的m個觀測值各自的閉合殘差，根據閉合殘差的大小從所有觀測值中選擇n+k個初選觀測值，n為觀測方程中設計矩陣的秩，m1，nm,k＝[0,1,…,m-n]；

步驟2、將m個觀測值按照各自真誤差絕對值的大小從小到大進行排序，獲得觀測值序列；

若|dΔ|_max＝||Δ_n+k|-|Δ_n+k-1||，將步驟1中選擇的n+k個初選觀測值作為復選觀測值，執行步驟3；否則，k＝k+1，從所述觀測值序列中選擇前n+k個觀測值作為初選觀測值，重新執行步驟2；

其中，dΔ為所述觀測值序列中相鄰觀測值的真誤差絕對值之間的差值，Δ_n+k為觀測值序列中第n+k個觀測值的真誤差，Δ_n+k-1為觀測值序列中第n+k-1個觀測值的真誤差；

步驟3、根據每一個觀測值的評價參數絕對值的大小從m個觀測值中選取y個觀測值，獲得一個y×1維的觀測值向量L_y；

其中，第i個觀測值的評價參數i＝(1,2,…,m)，Δ_i為第i個觀測值的真誤差，Δ_z為第z個觀測值的真誤差，z＝(1,2,…,n+k)；

步驟4、利用式VII計算當前像素點觀測值的解纏誤差向量

式VII中，L_r為r×1維的觀測值向量，r＝m-y，R為m×1維的觀測值向量L的平差因子矩陣；

步驟5、利用式X獲得當前像素點觀測值的真值L'：