[發(fā)明專利]基于D-InSAR技術(shù)與圖像加權(quán)疊加的沉降監(jiān)測方法

在審

• 申請?zhí)枺?/td>: 202110086848.0
• 申請日: 2021-01-22
• 公開（公告）號: CN112857310A
• 公開（公告）日: 2021-05-28
• 發(fā)明（設計）人: 楊學志;李克沖;董張玉
• 申請（專利權(quán)）人: 合肥工業(yè)大學
• 主分類號: G01C5/00
• 分類號: G01C5/00;G01S13/90
• 代理公司: 安徽省合肥新安專利代理有限責任公司 34101
• 代理人: 陸麗莉;何梅生

1.一種基于D-InSAR技術(shù)與圖像加權(quán)疊加的沉降監(jiān)測方法，其特征包括如下步驟：

步驟1.通過公共主影像選擇方法從SAR影像數(shù)據(jù)庫中選取多幅公共主影像；

步驟1.1.對SAR影像數(shù)據(jù)庫中的每幅SAR影像進行SAR數(shù)據(jù)預處理，得到校正后的SAR影像；

步驟1.2.利用式(1)計算第l幅影像的距離正則化的平均值d_l：

式(1)中，l表示SAR影像數(shù)據(jù)庫中的SAR影像的序號；B_i,l表示在影像l的正則化距離計算過程中第i幅校正后的SAR影像與第l幅校正后的SAR影像之間的垂直空間基線，B_c表示垂直空間基線的臨界值；T_i,l表示在第l幅校正后的SAR影像的正則化距離計算過程中第i幅校正后的SAR影像與第l幅校正后的SAR影像之間的時間基線，T_c表示時間基線的臨界值；f_i,l表示在第l幅校正后的SAR影像的正則化距離計算過程中第i幅校正后的SAR影像與第l幅校正后的SAR影像之間的多普勒質(zhì)心頻率差，f_c表示多普勒質(zhì)心頻率差的臨界值；N表示第l幅影像在干涉圖生成過程中參與生成的干涉像對個數(shù)；

步驟1.3.從SAR影像數(shù)據(jù)庫中選取前m幅最小的平均值所對應的SAR影像作為公共主影像；

步驟2.將m幅公共主影像分別與其余M-m幅校正后的SAR影像進行基線估算并生成連接圖；其中，M為SAR影像數(shù)據(jù)庫中SAR影像數(shù)量；

步驟3.根據(jù)數(shù)字高程模型DEM、SAR影像對應的軌道數(shù)據(jù)以及連接圖，對校正后的SAR影像進行處理并生成若干幅干涉圖；

步驟4.對所有干涉圖進行濾波處理并計算其相干性；

步驟5.對每幅濾波處理后的干涉圖進行相位解纏，得到解纏后的干涉圖，再對解纏后的干涉圖剔除相干性差的干涉圖，從而得到剔除后的干涉圖；對剔除后的干涉圖進行軌道精煉和重去平處理，得到重去平后的干涉圖；計算重去平后的干涉圖中相干點的數(shù)量，并選取相干點數(shù)量最高的干涉圖作為參考干涉圖；

步驟6.干涉圖的加權(quán)疊加；

步驟6.1.當重去平后的干涉圖中的相干點的數(shù)量小于參考干涉圖中相干點數(shù)量的一半時，將重去平后的干涉圖的權(quán)重設置為0，否則，設置為1；

步驟6.2.將所有重去平后的干涉圖根據(jù)其對應的權(quán)重進行疊加，從而得到疊加后的相位圖；

步驟7.將所述相位圖進行相位轉(zhuǎn)形變的處理并經(jīng)過地理編碼后生成形變圖，用于監(jiān)測沉降。