[發明專利]一種融合不同合成孔徑雷達干涉測量的地面沉降監測方法

有效

• 申請號: 202110347015.5
• 申請日: 2021-03-31
• 公開（公告）號: CN113008202B
• 公開（公告）日: 2023-02-21
• 發明（設計）人: 鄭曉慧;張鵬;夏錦;邊云云;李程
• 申請（專利權）人: 中鐵上海設計院集團有限公司;南京智慧巖土工程技術研究院有限公司;南京工業大學
• 主分類號: G01C5/00
• 分類號: G01C5/00;G01S13/90
• 代理公司: 南京瑞弘專利商標事務所(普通合伙) 32249
• 代理人: 沈廉

1.一種融合不同合成孔徑雷達干涉測量的地面沉降監測方法，其特征在于該監測方法的步驟如下：

步驟1.獲得N幅合成孔徑雷達SAR圖像數據，構成SAR圖像數據集合，所述N大于等于20；所述圖像數據為SAR對同一個地表觀測場景進行成像并完成了圖像配準的數據；

步驟2.對步驟1中獲得的SAR圖像數據集合進行永久散射體合成孔徑雷達干涉測量PS-InSAR處理，解算出地面沉降監測結果，獲得滿足相干性閾值C為篩分條件的PS點；

步驟3.整理所有PS點的經緯度坐標(x,y)和每個反映PS點計算結果的可靠度或置信度對應相干值c，采用重心法求得PS點空間分布的重心坐標G（a,b）；

步驟4.將重心坐標G（a,b）作為DBSCAN聚類算法的初始聚類中心點，對所有PS點的經緯度坐標(x,y)進行聚類分析，將所有PS點的空間分布分為密集區和稀疏區兩個部分；

步驟5.求得PS點密集區域離散點的最小幾何邊界，建立不規則三角網格并以三角網格表面最大邊長D為約束條件對最小幾何邊界進行修正，修正后的最小幾何邊界即為PS點分布密集區和稀疏區的分界線L；

步驟6.對步驟1中獲得的SAR圖像數據集合進行短基線集合成孔徑雷達干涉測量SBAS-InSAR處理，解算出地面沉降監測結果，將SBAS-InSAR方法獲得的地面沉降監測結果的柵格數據轉換為矢量數據；

步驟7.根據分界線L，將稀疏區內的PS-InSAR監測成果數據更新為對應區域的SBAS-InSAR方法的矢量數據，獲得最終整個研究區域的地面沉降監測結果。

2.根據權利要求1所述的一種融合不同合成孔徑雷達干涉測量的地面沉降監測方法，其特征在于：所述步驟2中相干性閾值C≥0.9。

3.根據權利要求1所述的一種融合不同合成孔徑雷達干涉測量的地面沉降監測方法，其特征在于：步驟3中所述的重心法為：

其中，a,b為重心坐標的橫坐標和縱坐標；x,y為每一個PS點的橫坐標和縱坐標；c為每一個PS點對應的相干值。

4.根據權利要求1所述的一種融合不同合成孔徑雷達干涉測量的地面沉降監測方法，其特征在于：步驟5中所述最大邊長D根據檢測點離散型程度取4000～5000m。