亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        多視線向D-InSAR三維地表形變解算中GPS約束定權(quán)

        2018-08-31 05:46:02王志偉劉國林陶秋香于勝文周傳義
        測繪通報 2018年8期
        關(guān)鍵詞:視線觀測方向

        王志偉,劉國林,陶秋香,于勝文,王 珂,周傳義

        (山東科技大學測繪科學與工程學院,山東 青島 266590)

        合成孔徑雷達干涉測量(interferometric synthetic aperture radar,InSAR)技術(shù)以其高精度、高空間分辨率的地表形變監(jiān)測能力,目前已經(jīng)被廣泛應用于監(jiān)測各種因素導致的地表形變[1-7]。隨著越來越多的衛(wèi)星發(fā)射升空,同一地區(qū)可以獲取多種不同傳感器SAR影像,進而可以獲取不同雷達視線向地表形變監(jiān)測結(jié)果(多視線向D-InSAR形變量)。將多視線向D-InSAR地表形變量作為觀測量,以衛(wèi)星軌道方位角和入射角等參數(shù)構(gòu)建系數(shù)矩陣,建立三維地表形變解算模型,可以獲取真實三維地表形變[12]。然而,由于受到不同誤差因素的影響,多視線向D-InSAR地表形變監(jiān)測結(jié)果的精度不盡相同。因此,基于多視線向D-InSAR形變量聯(lián)合平差,以獲取高精度三維形變解算結(jié)果尤為重要。其研究關(guān)鍵點是如何合理確定各觀測量的權(quán)重。

        本文提出一種基于GPS約束的多視線向D-InSAR三維地表形變解算中定權(quán)方法。該方法以高精度GPS形變監(jiān)測結(jié)果作為約束,通過對多視線向D-InSAR觀測值進行評價確定其權(quán)重。本文以2009年意大利L’Aquila地震三維地表形變解算為例,采用GPS約束的定權(quán)方法實現(xiàn)三維形變解算。試驗結(jié)果表明,本文提出的GPS約束定權(quán)的多視線向D-InSAR三維形變解算方法可獲取更可靠的解算結(jié)果,可以解決不同精度觀測量聯(lián)合平差權(quán)重確定存在的定權(quán)難、定權(quán)適用性差、精度低等突出問題。

        1 多視線向D-InSAR三維地表形變GPS約束定權(quán)方法

        1.1 多視線向D-InSAR三維地表形變解算

        SAR衛(wèi)星的成像幾何關(guān)系如圖1所示。AZ向表示衛(wèi)星飛行方向在地面上的投影,GR向表示衛(wèi)星視線向在地面的投影,α表示衛(wèi)星方位角,θ表示SAR影像入射角,U、E和N分別代表大地測量坐標系中的垂直、東西和南北方向,dU,dE和dN分別代表對應方向上的形變信息。

        圖1 SAR衛(wèi)星三維成像示意圖

        R=Bd

        (1)

        式中,B為三維形變模型系數(shù)矩陣,具體可以表示為

        (2)

        根據(jù)最小二乘原理,可以解算出三維形變信息

        d=(BTPB)-1BTPR

        (3)

        1.2 基于GPS約束的定權(quán)方法

        針對多視線向D-InSAR三維地表形變解算中存在的定權(quán)問題,本文提出以GPS精密形變觀測量為基準數(shù)據(jù),對多視線向D-InSAR觀測量進行精度評價,綜合估算出軌道和地形誤差及噪聲等因素引起的誤差。由于GPS測站數(shù)量有限,無法逐像元進行精度評價。顧及軌道和地形及噪聲等因素引起的誤差對形變監(jiān)測結(jié)果影響的特點,建立聯(lián)合擬合模型并計算擬合參數(shù)。根據(jù)確定出的擬合模型估算多視線向D-InSAR每個像元觀測量的方差,以計算出的方差的倒數(shù)確定權(quán)矩陣,建立加權(quán)最小二乘三維地表形變解算模型。具體步驟如下:

        (1) 將GPS三維形變觀測量換算到D-InSAR視線方向,如下

        (4)

        (2) 根據(jù)雙觀測值之差標準差計算公式,計算出GPS測站對應D-InSAR觀測量方差

        (5)

        (3) 根據(jù)步驟(2)計算出的雙觀測值之差標準差及對應像元相干值,按照下式計算方差

        (6)

        式中,γi為第i個像元相干系數(shù);E(γGPS)為GPS觀測站對應的SAR影像像元相干系數(shù)均值。

        (4) 以步驟(3)估算出的D-InSAR像元的方差為依據(jù),充分考慮軌道和地形誤差對形變量影響的特點,建立方差與像元點的位置和高程聯(lián)合擬合模型

        (7)

        式中,等號右側(cè)前6項為軌道誤差項;x表示第i個像元方位向坐標;y表示第i個像元距離向坐標;后兩項表示地形誤差項,z表示坐標為(x,y)像元的高程值,ai(i=1,2,…,5)和bj(j=1,2)分別表示待估參數(shù)。

        估算出上式中系數(shù)參數(shù)后便可以確定出擬合模型,按照上述模型估算每個像元D-InSAR形變監(jiān)測結(jié)果方差。

        (5) 根據(jù)D-InSAR形變監(jiān)測結(jié)果方差,按照下式定權(quán)

        (8)

        按照式(8)計算出權(quán)陣P之后,按照式(3)進行三維地表形變解算。

        2 L’Aquila地震多視線向D-InSAR地表形變監(jiān)測

        2.1 試驗區(qū)與數(shù)據(jù)源

        2009年4月6日意大利中部城市L’Aquila發(fā)生了Mw6.3級地震,整個意大利中部地區(qū)均有震感。本文以該地震震中區(qū)域為研究區(qū)驗證定權(quán)方法的有效性。考慮干涉對的選取原則(足夠短的垂直基線和時間基線),選取了4個干涉像對SAR數(shù)據(jù),分別用P1、P2、P3和P4表示,影像覆蓋范圍如圖2所示。同期的GPS數(shù)據(jù)來源于文獻[13]。研究區(qū)范圍內(nèi)共有48個點,如圖2所示,其中黑點表示用來估算聯(lián)合擬合模型參數(shù)的GPS數(shù)據(jù),白點用來驗證解算精度。

        表1 選取的SAR數(shù)據(jù)基本信息

        圖2 研究區(qū)概況

        2.2 數(shù)據(jù)處理

        以震前、后獲取的SAR影像分別作為主、輔影像,采用雙軌D-InSAR技術(shù)對數(shù)據(jù)進行處理,對主、輔影像進行配準、重采樣和共軛相乘處理,獲取干涉相位圖。采用30 m分辨率的SRTM DEM數(shù)據(jù)去除地形相位的影響。采用Goldstein濾波方法對干涉圖進行噪聲濾除,獲取最終的干涉圖。采用最小費用流的方法進行相位解纏。解纏之后的干涉圖只包含形變相位、大氣相位和殘余噪聲相位。其中,大氣相位對干涉圖的影響主要體現(xiàn)在對流層延遲方面,往往跟地形高程有一定的相關(guān)性,通過GAMMA軟件的相關(guān)算法可以一定程度上去除大氣誤差相位。

        圖3為地理編碼后的形變圖。從圖中可以看出,多視線向D-InSAR形變監(jiān)測結(jié)果在形變趨勢和量級上基本一致,震中和斷裂帶周邊形變?yōu)樨撝当硎驹搮^(qū)域地面下沉(呈現(xiàn)深色)。多視線向D-InSAR形變監(jiān)測結(jié)果在雷達視線向上下沉的一致性說明拉奎拉地震是以地表下沉為主。根據(jù)地質(zhì)調(diào)查的結(jié)果顯示,拉奎拉地區(qū)所在的阿爾卑斯地區(qū)斷層走向為西北-東南方向,大小在40~50之間的西南傾向,正斷層。

        圖3 不同軌道SAR影像獲取的形變圖

        3 L’Aquila地震三維形變解算

        3.1 三維地表形變解算

        獲取了多視線向D-InSAR地表形變信息,按照第一部分介紹的數(shù)據(jù)處理方法進行三維地表形變解算。值得注意的是兩種數(shù)據(jù)形變監(jiān)測結(jié)果具有不同的空間分辨率,本文將ASAR數(shù)據(jù)監(jiān)測結(jié)果重采樣至ALOS數(shù)據(jù)分辨率。入射角和方位角信息見表1。用來計算擬合參數(shù)的GPS數(shù)據(jù)具體分布如圖2中白色點所示。解算出的三維地表形變結(jié)果如圖4所示。

        圖4 三維地表形變解算結(jié)果

        從解算出的三維地表形變結(jié)果中可以看出,在震中主形變區(qū)內(nèi),地表發(fā)生了東南方向的水平錯動。在斷層周邊的三維形變場中表現(xiàn)出了比較明顯的形變特征。在東西方向形變場的下半部分,地表向西發(fā)生了移動,上側(cè)發(fā)生了向東的移動,表現(xiàn)出比較明顯的形變錯動。在南北方向下半部分,地表向北發(fā)生了移動,上側(cè)發(fā)生了向南的移動,同樣表現(xiàn)出比較明顯的形變錯動。在垂直方向的右上方,地表表現(xiàn)為上升的運動,而在下方表現(xiàn)為沉降,垂直方向上也發(fā)生了明顯了錯動。

        3.2 精度驗證

        為了說明本文提出定權(quán)方法的有效性,通過解算出的三維結(jié)果與真實GPS數(shù)據(jù)進行對比。GPS點的位置如圖2白點位所示。將解算結(jié)果與GPS數(shù)據(jù)的作差值比較。圖5給出了對比結(jié)果,曲線越接近于0說明解算結(jié)果越準確。

        圖5 結(jié)果對比

        從圖5可以看出,相對于南北方向來說,東西方向和垂直方向2個方向的縱坐標取值范圍小,這說明東西和垂直向的解算精度比南北方向高(縱軸表示與GPS監(jiān)測出的三維形變的差值)。導致這一結(jié)果的原因有兩點:①衛(wèi)星近南北方向飛行,導致雷達視線向?qū)δ媳狈较蛐巫儾幻舾校诮馑氵^程中觀測值中的噪聲被放大,從而導致南北方向解算精度差[17];②GPS監(jiān)測結(jié)果出垂直向精度比水平方向監(jiān)測精度要差[18],這會對GPS約束改進定權(quán)方法解算結(jié)果產(chǎn)生影響。

        為了進一步說明本文提出定權(quán)方法的有效性,表2計算解算出的三維地表形變與收集的GPS三維地表形變對比的均方根誤差(RMSE)及統(tǒng)計了對比的最大差值。

        4 結(jié) 論

        本文針對多視向D-InSAR三維地表形變解算定權(quán)問題,提出了一種基于GPS約束的多視線向D-InSAR三維地表形變監(jiān)測定權(quán)方法。以2009年拉奎拉地震實測數(shù)據(jù)為例進行了可行性驗證,可以得出以下結(jié)論:

        (1) 兩種數(shù)據(jù)監(jiān)測出的研究區(qū)形變結(jié)果基本一致,但兩種數(shù)據(jù)獲取的時間間隔、數(shù)據(jù)軌道和波段會導致監(jiān)測結(jié)果在形變位置、范圍、形變量等方面略有不同。

        表2 對比結(jié)果定量統(tǒng)計

        (2) 從解算出的三維地表形變分布圖來看,垂直方向和東西方向的解算結(jié)果相對較高,南北方向解算結(jié)果受噪聲影響較大。

        (3) 相對于東西方向來說,垂直方向的均方誤差較大,解算精度差一些,這表明GPS低精度的垂直方向形變監(jiān)測結(jié)果會對三維形變解算精度產(chǎn)生影響。

        (4) 通過收集到的GPS監(jiān)測數(shù)據(jù)對解算結(jié)果進行精度分析與驗證,對比發(fā)現(xiàn)解算出的東西方向均方根誤差為0.036 4 m,南北方向均方根誤差為0.513 5 m,垂直方向均方根誤差為0.439 m。

        綜上所述,利用本文所提出的GPS約束定權(quán)方法能正確確定多視線向D-InSAR觀測量的權(quán),且計算簡單。采用GPS約束改進定權(quán)方法能提高三維解算精度。

        猜你喜歡
        視線觀測方向
        觀測到恒星死亡瞬間
        軍事文摘(2023年18期)2023-11-03 09:45:42
        2022年組稿方向
        計算機應用(2022年2期)2022-03-01 12:33:42
        要去就去視線盡頭的山
        2021年組稿方向
        計算機應用(2021年4期)2021-04-20 14:06:36
        2021年組稿方向
        計算機應用(2021年1期)2021-01-21 03:22:38
        你吸引了我的視線
        天測與測地VLBI 測地站周圍地形觀測遮掩的討論
        可觀測宇宙
        太空探索(2016年7期)2016-07-10 12:10:15
        高分辨率對地觀測系統(tǒng)
        太空探索(2015年8期)2015-07-18 11:04:44
        當代視線
        當代貴州(2015年19期)2015-06-13 09:42:32
        狠干狠爱无码区| 久久人人玩人妻潮喷内射人人 | 亚洲av美国av产亚洲av图片| 欧美成人午夜精品久久久| 日本午夜免费福利视频| 欧美中文在线观看| 尤物AV无码色AV无码麻豆| 日本成人在线不卡一区二区三区| 麻豆视频av在线观看| 蜜桃视频免费进入观看| 精品国产青草久久久久福利| 亚洲综合无码无在线观看| 日韩毛片在线| 抖射在线免费观看视频网站| 日韩精品国产精品亚洲毛片| 艳妇臀荡乳欲伦交换h在线观看| 国产精品无码专区av在线播放| 久久福利青草精品免费| A阿V天堂免费无码专区| 午夜精品人妻中字字幕| 偷拍色图一区二区三区| 乱中年女人伦| 久久久天堂国产精品女人| 国产高清无码在线| 国产精品亚洲综合色区丝瓜| 成人自拍三级在线观看| 国产内射爽爽大片| 色拍自拍亚洲综合图区| 久久亚洲sm情趣捆绑调教| 亚洲中文久久久久无码| 日本老熟女一区二区三区| 亚洲综合精品亚洲国产成人| 性猛交╳xxx乱大交| 欧美丰满大乳高跟鞋| 亚洲在战AV极品无码| 中文文精品字幕一区二区| 九九在线中文字幕无码| 男同gay毛片免费可播放| 中文字幕国产精品中文字幕| 国产免费一区二区三区在线视频| 黄片视频免费观看蜜桃|