啪啪国产视频_一级高清毛片_日韩精品中文字幕在线观看_免费视频一区二区三区四区

用戶名: 密碼: 驗證碼:         注冊
返回首頁

基于InSAR技術的變形監測方法與研究進展

時間:2019-07-25 11:39來源:未知 作者:監測人 點擊:
變形監測是星載InSAR技術應用最為成熟的領域之一。本文首先介紹了InSAR變形監測的基本原理和衛星數據來源;然后對InSAR變形監測方法進行了系統性的分類,分析了D-InSAR、PS-InSAR、SBAS-InSAR、DS-InSAR和MAI等方法的技術特點和適用范圍;進而從應用的角度分析了InSAR技
  

InSAR變形監測方法與研究進展

朱建軍, 李志偉, 胡俊     

中南大學地球科學與信息物理學院, 湖南 長沙 410083

基金項目:國家重點基礎研究發展計劃(973計劃)(2013CB733303);國家自然科學基金(41371335;41404011);高分辨率對地觀測系統重大專項(民用部分)(03-Y20A11-9001-15/16)

第一作者簡介:朱建軍(1962-), 男, 教授, 博士生導師, 研究方向為測量平差與InSAR數據處理。E-mail:zjj@mail.csu.edu.cn

 

星載合成孔徑雷達干涉測量(interferometric synthetic aperture radar, InSAR)因其高精度、高分辨率、全天候等優點已迅速成為常用的大地測量技術之一,旨在通過計算兩次過境時SAR影像的相位差來獲取數字高程模型[1-2]。隨之而來的差分雷達干涉技術(differential interferometric synthetic aperture radar, D-InSAR)則是通過引入外部DEM或三軌/四軌差分實現了地表變形監測。除了地形因素外,時空失相干、大氣、軌道等也是影響DInSAR技術的形變監測精度的主要因素[3-4]。與此同時,衛星硬件的不斷發展,重返周期越來越短會導致微小形變與噪聲之間的混疊,空間分辨率越來越高也會導致噪聲更加復雜,這些都對常規D-InSAR的形變監測提出更高的挑戰。此外,D-InSAR技術獲取的形變是地表真實三維形變在雷達視線方向(line-of-sight, LOS)上的投影,目前多通過除以入射角的余弦值得到垂直向的形變。但是,這需要假設水平方向上無形變[5],很多情況下并不滿足,從而導致形變信號的錯誤解譯。

為了突破DInSAR技術的這些限制,學者們提出了多時相InSAR技術(multi-temporal InSAR, MT-InSAR)來進行高精度的形變監測,如永久散射體(persistent scatterer InSAR, PS-InSAR)、小基線集(small baseline subsets InSAR, SBAS-InSAR)和分布式散射體(distributed scatterer InSAR, DS-InSAR)等。此外,為了彌補D-InSAR或MT-InSAR技術只能獲取LOS形變的缺陷,學者又提出多孔徑InSAR(multi-aperture InSAR, MAI)技術來獲取方位向(即衛星飛行方向)上的形變信息。以上這些InSAR技術已廣泛應用于各類形變監測中,如城市地面沉降、礦區沉降、地震及板塊運動、火山噴發、基礎設施變形、冰川漂移、凍土形變、滑坡等。

本文旨在通過分析以上這些InSAR技術及其在各應用領域的研究現狀,使讀者能夠對InSAR變形監測有更為系統性的理解。本文章節安排如下:第1節介紹InSAR技術與現狀,包括雷達干涉測量原理和SAR衛星現狀;第2節介紹目前主流的幾種InSAR變形監測方法及進展;第3節介紹InSAR變形監測在重點應用領域的研究現狀;第4節則對InSAR變形監測的挑戰進行分析。

1 InSAR技術與現狀1.1 雷達干涉測量原理

合成孔徑雷達干涉測量技術(InSAR)集合成孔徑雷達技術與干涉測量技術于一體。合成孔徑雷達(SAR)是一種主動式微波遙感,用來記錄地物的散射強度信息及相位信息,前者反映了地表屬性(含水量、粗糙度、地物類型等),后者則蘊含了傳感器與目標物之間的距離信息。干涉的基本原理是同一區域兩次或多次過境的SAR影像的復共軛相乘,來提取地物目標的地形或者形變信息。雷達干涉的模式分為:沿軌道干涉法、交叉軌道干涉法、重復軌道干涉法。其中利用重軌干涉最為常用的方式,以此為例,所得到的干涉相位表達式為[4]  (1)

式中,φflat為平地相位;φtopo為地形相位,該相位可以用來恢復地形信息;φdef為地表形變引起的相位;φatmo為大氣延遲相位;φnoise為觀測噪聲引起的相位。

 將平地相位、地形相位、噪聲相位、大氣相位去除,即可得到地表形變相位,這樣的雷達干涉技術被稱為雷達差分干涉技術(D-InSAR),文獻[2]首次論證該項技術可用于探測厘米級的地表形變。其幾何關系如圖 1所示。

 

 

圖 1 差分干涉測量幾何關系Fig. 1 Geometric relation of differential interferometry

從圖3中可以導出以下關系

  (2)

式中,Δr為所要求解的地面目標點的LOS向形變信息;λ為雷達波長。

 1.2SAR衛星的發展現狀

 1978年,美國發射了L波段SEASAT衛星,在1981、1984年又相繼發射了SIR-A和SIR-B衛星,為最初的研究提供了少量的星載數據。隨后其他國家也相繼研究并研發了不同波段的中低分辨率的SAR衛星。1991年7月、1995年4月歐洲空間局(ESA)相繼發射了C波段ERS-1和ERS-2衛星,采用相同的軌道參數,組成串行飛行模式,為對地觀測提供了大量的數據,極大地推動了InSAR技術的發展。1992年日本空間局發射了L波段JERS-1衛星,長波段的穿透力強,使其在InSAR監測中具有很強的優勢,但JERS-1衛星并不是專門針對干涉測量而設計的,軌道誤差較為嚴重,影響了它的適用性。1995年,加拿大發射了C波段RADARSAR-1衛星,獲取了覆蓋全球的大量數據。21世紀以后,SAR衛星系統的硬件得到更大的改進,專門針對干涉測量的衛星陸續發射成功。2000年美國實施的SRTM計劃,利用航天飛機獲取了全球80%的30m和90m分辨率的DEM數據,對差分干涉測量提供了可靠的外部DEM數據。歐空局2002年3月發射的ERS系列后繼衛星ENVISAT不僅提供了SAR數據,還提供了MERIS水汽數據,對InSAR數據中大氣延遲誤差的改正提供了極大的幫助。2006年1月日本空間局發射了L波段的ALOS-1衛星,因其長波長穿透性強、抗失相干的特性在地震、礦山等方面的監測得到了廣泛的應用。ERS-1/2、ENVISAT、ALOS-1、RADARSAT-1獲取豐富的中低分辨率SAR數據成為21世紀早期InSAR技術發展的主要數據源,大部分的研究成果也是基于這些數據得到。

 

在上述的中低分辨率的SAR衛星陸續停止服務后,新一代高分辨率衛星相繼發射,而且這些SAR衛星系統也從單一極化、模式、波段、固定入射角向多極化、多模式、多波段、可變入射角發展。由此,SAR衛星系統迎來了一個全新的時代。2007年德國航空航天局成功發射了X波段TerraSAR-X衛星,隨后在2010年發射了其姊妹星TanDEM-X,構成分布式協同工作模式,可以提供全球高精度的數據及數字高程模型,分辨率最高可達0.25m。意大利空間局設計的X波段COSMO-SkyMed星座包括4顆星,分別于2007年6月、2007年12月、2008年10月、2010年11月成功發射,分辨率高達1m,重訪周期短,在風險預警、災害管理等方面發揮著重要的作用。2007年12月加拿大發射的C波段RADARSAT-2衛星繼承了RADARSAT-1衛星的工作模式,還增加了多極化成像等,增強了其在地形監測領域的觀測能力。2014年4月歐空局成功發射了首顆環境監測衛星C波段Sentinel-1A衛星,并于2016年4月發射了Sentinel-1B衛星,對全球大范圍的地質、環境災害的監測提供了豐富的數據。2014年5月日本成功發射的L波段ALOS-2衛星,空間分辨率最高達到1m,其在大型地質災害(如地震)監測領域具有獨一無二的優勢。與此同時,我國也在積極的研發SAR衛星,并于2016年8月10日在太原成功發射了高分三號(GF-3)衛星,實現了我國衛星SAR影像干涉測量零的突破。

 

現有的SAR衛星系統具體參數如表1所示。在未來,SAR衛星系統的設計將會針對不同的研究領域而進行規劃,例如,歐空局規劃在2020年發射P波段BIOMASS衛星是為森林樹高測繪和生物量反演而設計;德國航天局預計在2022年發射的L波段TanDEM-L衛星,主要針對全球陸表動態變化監測而設計。

 

表1現有SAR衛星系統具體參數Tab.1SpecificparametersofcurrentSARsatellites




 

 

2InSAR變形監測方法與進展

2.1D-InSAR方法

 

D-InSAR是在傳統的InSAR技術上發展起來的方法,主要通過引入外部DEM去除InSAR獲取的干涉圖中的地形相位,進而得到差分干涉圖[2]。而這差分干涉相位中除了地形相位外,還存在平地效應、形變信號、大氣以及噪聲成分。由于該差分干涉相位仍然是纏繞的,因此需對剩余的相位成分進行整周模糊度求解。具體地,通過精確的軌道信息來去除平地效應,利用一定的濾波來降噪,通過外部數據或模型來去除大氣相位[6],通過相應的相位解纏方法來獲取整周模糊度[7]。自Gabriel在1989年首次利用D-InSAR成功獲取了農田區的形變信息后,學者們相繼對DInSAR技術做了一系列的改進,如濾波算法的提出和改進[8-9]、軌道相位的去除[10]、相干性的穩健性估計、電離層的去除[12]等。同時,D-InSAR也從最初的地震、火山等大范圍地表形變監測慢慢發展到局部形變監測,如結合角反射器的礦區、滑坡、油氣田以及城市地表形變。雖然D-InSAR已取得了一定的成果,但因監測精度仍受時空失相干、大氣等因素的影響,因此基于D-InSAR發展起來的MT-InSAR是現今InSAR技術監測變形的研究熱點。

 

2.2PS-InSAR方法

 

文獻[13—14]首次提出了PS-InSAR方法,其基本思想是:第一,利用覆蓋同一研究區的多景單視SAR影像,選取其中一景SAR影像作為主影像,其余SAR影像分別與主影像配準,依據時間序列上的幅度和(或)相位信息的穩定性選取永久散射體(PersistentScatterer,PS)目標;第二,經過干涉和去地形處理,得到基于永久散射體目標的差分干涉相位,并對相鄰的永久散射體目標的差分干涉相位進行再次差分;第三,根據兩次差分后的干涉相位中各個相位成分的不同特性,采用構建形變相位模型和時空濾波或方式估計形變和地形殘余信息。

 

PS-InSAR技術不是針對SAR影像中的所有像元進行數據處理,而是選取在時間上散射特性相對穩定、回波信號較強的PS點作為觀測對象。這些PS點通常包括人工建筑物、燈塔、裸露的巖石以及人工布設的角反射器等[14]。PS點的準確選取可以確保即便在干涉對的時間或空間基線很長的條件下(甚至達到臨界基線),PS點依然呈現出較好的相干性和穩定性。常用PS點選取方法包括振幅離差閾值法[14]、相干系數法[16]、相位分析法[17],以及這些方法的組合等。而形變和地形殘差解算通常采用解空間搜索法[14]、LAMBDA方法[18]和StamPS中的三維解纏法[19]等。

 

目前,PS-InSAR技術已在多個城市的高分辨率地面沉降監測中得到廣泛應用,特別是城市重點基礎設施的高分辨率形變監測。通過對比同期的水準和GPS測量數據,證實了PS-InSAR技術具有較高的可靠性(mm級的精度)[20],因此對掌握城市地面沉降動態規律以及分析地面沉降的成因等具有重要意義。

 

然而PS-InSAR技術也存在自身的缺陷。首先,其通常要求覆蓋同一區域的SAR影像數目較多(通常要求大于25景),便于保證模型解算的可靠性。其次,PS-InSAR技術由于是基于大量PS點的迭代回歸或網平差計算,運算效率不高,因此不適合大范圍地區高分辨率的地面沉降監測。

 

2.3SBAS-InSAR方法

 

SBAS-InSAR是一種基于多主影像的InSAR時間序列方法,只利用時空基線較短的干涉對來提取地表形變信息,文獻[21]提出并利用該技術進行意大利南部CampiFlegrei火山口和Naples市的時間序列形變估計,證明了該方法具有能夠準確估計形變的能力。SBAS-InSAR技術克服了PS-InSAR因選取一幅影像作為公共主影像而引起的部分干涉圖相干性較差的不足,同時降低了對SAR數據的需求量,運算效率較高。

 

SBAS-InSAR技術是以多主影像的干涉對為基礎,基于高相干點恢復研究區域的時間序列形變信息,其原理如下:首先對覆蓋某個地區的不同時間段的多景SAR影像計算時間空間基線,選擇恰當的時空基線閾值選取干涉對;然后對選取的干涉對進行差分干涉處理并進行相位解纏;最后根據自由組合的干涉圖形成子集的情況,對所有干涉圖組成相位方程采用最小二乘法或者SVD方法進行形變參數的估計。在實際處理中會采用時空濾波的方法去除大氣延遲的影像分離出非線性形變,估算的低頻形變與此非線性形變的總和即為整個研究區的形變信息。

 

針對SBAS-InSAR技術,各國學者展開了許多研究,在高相干點選取方面主要有相位穩定性選點法[22]、振幅離差指數選點法[13]、空間相干性[23]等方法;在形變觀測的數學模型選取方面,主要有線性模型[23]、二階多項式形變模型[21]、周期性形變項[24]等;在參數估算方面主要有最小二乘法[23]、SVD法[21]、基于L1范數的解算方法[26]等。

 

為了提高高質量相干點的密度、減弱誤差的影像、改善形變參數估計結果的精度,國內外學者提出了一些融合PS-InSAR和SBAS-InSAR優點的算法,如文獻[23]提出了基于最小二乘網絡的平差法,利用少量的SAR數據獲取干涉圖中高相干點的地表形變信息。文獻[27]應用PS-InSAR和SBAS-InSAR融合技術于Eyjafjallajokull火山的形變監測中,獲取了該區域的時間序列形變結果。文獻[28]提出了時域相干點法(TCP-InSAR),對形變參數、軌道誤差同時進行估計,提高了監測精度。

2.4DS-InSAR方法

 

自2011年文獻[29]發布第二代永久散射體技術SqueeSAR以來,時序InSAR領域的研究熱點逐漸轉向對分布式目標(distributedscatterer,DS)的探索。與PS目標的物理屬性不同,DS是指在雷達分辨率單元內沒有任何散射體的后向散射占據統治地位的點目標[30]。事實上,利用DS-InSAR進行形變監測的概念早在SBAS技術和StamPS技術中已經建立[21,27],當時學界更多的是以相干、非相干目標加以區分,而弱化了PS與DS之間的物理界限。因此,從數據處理的角度來說,以SBAS為代表的小基線技術與SqueeSAR技術及其變種有本質不同,但它們又同時強調DS目標的信噪比,因而均屬于相位濾波。

 

SqueeSAR的技術要點是:①通過同質點選取算法增強時序InSAR協方差矩陣的估計精度,并同時輔助PS與DS目標的分離;②通過相位優化算法從協方差矩陣中恢復時序SAR影像的相位。在第1個步驟中,其前提條件是相同SAR影像質地的像素具有相同相位中心,因此在時序統計推斷的框架下,選取具有相同SAR統計分布的像素參與平均不僅可以提升相位信噪比,還能保留圖像的空間分辨率[32]。相比之下,SBAS多采用常規多視處理或空間(自適應)濾波,是一種以犧牲空間分辨率為代價換取相位質量提升的方法[29]。在復雜形變特征環境下,這類方法極易損失細節,并伴隨形變非形變區域的誤讀。在第2個步驟中,SqueeSAR在樣本協方差矩陣服從復Wishart分布的基礎之上,采用極大似然估計方法得到優化后的時序相位。由于似然估計量無解析表達式,需采用非線性優化“擠壓”待估參數,便成為“Squeeze”一詞命名的由來。值得注意的是,SqueeSAR一方面運用了所有干涉對的信息,另一方面卻并未從濾波后的干涉相位中直接提取時序相位,這與大多數時序InSAR技術采用相位三角關系從濾波后的干涉相位中直接解算時序相位有本質區別[29]。在優化DS之后,與PS目標一起融入傳統PS-InSAR數據處理框架就可以獲得精度更高、空間分辨率增強的時序形變產品。

 

SqueeSAR技術也有其缺陷。在同質選點算法中,使用的KS假設檢驗功效低,特別在小樣本條件下易使選取的同質集合中包含許多異質點,而且計算效率低下[32]。因此,近年來國內外學者紛紛提出了計算更為高效或估計精度更高的統計推斷方法,具有代表性的是AD檢驗[33]、似然比檢驗[32]、自適應檢驗和置信區間估計等[34]。在相位優化算法中,計算效率是該方法的主要詬病之處,使得大范圍的精細監測幾乎無法施展,備擇的方案是文獻[35]提出的協方差矩陣奇異值分解方法。在優化可靠性方面,可采用抗差性更強的估計量,如M估計量等[36]。

 

上述改進算法與SqueeSAR技術一起構成了當今DS-InSAR的雛形,旨在提高計算效率、增加觀測密度和改善估計精度。廣義的,除基于PS目標之外的時序InSAR技術都或多或少地應用了DS的概念,而聯合PS-InSAR和DS-InSAR或許是提高最終產品精度的最佳方案[27,29]。

 

頂一下
(0)
0%
踩一下
(0)
0%
------分隔線----------------------------
最新評論 查看所有評論
發表評論 查看所有評論
請自覺遵守互聯網相關的政策法規,嚴禁發布色情、暴力、反動的言論。
評價:
表情:
用戶名: 密碼: 驗證碼:
發布者資料
監測人 查看詳細資料 發表文章 加為好友 用戶等級:注冊會員 注冊時間:1970-01-01 08:01 最后登錄:2023-04-19 18:04
推薦內容
精彩推薦
精彩推薦

關于本站|免責聲明|廣告服務|網站地圖|網站搜索|友情鏈接|聯系我們|RSS訂閱

Copyright © 2008-2023 www.www.mocsil.com All Rights Reserved
監測人 www.www.mocsil.com 安全監測行業最大的門戶網站
稿件聯系E-mail: admin@www.mocsil.com 廣告聯系:QQ:32334427[請詳細說明來意] 電話:13687374912【非誠勿擾】
免責聲明:所有論文及相關資料均源于網上的共享資源及期刊共享,或者來源于會員的共享,監測人網站只是一個資源共享和交流的非盈利平臺。
如有侵犯您論文的版權或其他有損您利益的行為,請聯系批評與指出   湘ICP備14001333號-1

湘公網安備 43011102001622號