時(shí)序合成孔徑雷達(dá)干涉測(cè)量的多級(jí)化解算

胡鳳鳴, 伍吉倉(cāng)

(同濟(jì)大學(xué) 測(cè)繪與地理信息學(xué)院,上海 200092)

上海磁懸浮全長(zhǎng)30 km,是世界第一條商用的高架磁懸浮運(yùn)營(yíng)線.磁懸浮工程在建設(shè)和運(yùn)營(yíng)期間對(duì)沿線的地面沉降控制要求嚴(yán)格,每一跨徑的沉降差異須小于3 mm[1]. 上海浦東地區(qū)第四紀(jì)松散土層厚300余米,極易受地下水開(kāi)采以及工程活動(dòng)的影響誘發(fā)地面沉降[2].為保證磁懸浮的安全運(yùn)營(yíng),需要進(jìn)行高頻高密度高精度的沉降監(jiān)測(cè).傳統(tǒng)的精密水準(zhǔn)測(cè)量由于空間密度較低,測(cè)量時(shí)間間隔大,難以滿足磁懸浮沿線沉降監(jiān)測(cè)需求.

合成孔徑雷達(dá)干涉測(cè)量(interferometric synthetic aperture radar,InSAR)是一種新興的空間對(duì)地觀測(cè)技術(shù),能夠全天候、大范圍、高精度獲取地球表面地形和形變信息.針對(duì)傳統(tǒng)InSAR技術(shù)時(shí)空間失相干以及大氣延遲影響[3]提出的多時(shí)序InSAR(MT-InSAR)可以有效提高測(cè)量精度.研究表明,通過(guò)分析處理高質(zhì)量的相干點(diǎn)可以獲得亞毫米級(jí)的形變監(jiān)測(cè)精度[4].該技術(shù)已經(jīng)被成功應(yīng)用于地殼形變[5]、城市地表[6-7]以及大型結(jié)構(gòu)體的監(jiān)測(cè)中[8-9].X波段的TerraSAR-X/TanDEM-X以及COSMO-Skyed系統(tǒng)獲取的高分辨率SAR數(shù)據(jù)使得該技術(shù)可以獲取更加精細(xì)的地表信息[10-11],這對(duì)于城市區(qū)域內(nèi)大型結(jié)構(gòu)體健康監(jiān)測(cè)具有重要意義.

MT-InSAR技術(shù)僅使用在時(shí)間序列上始終保持高相干性的散射體,依據(jù)干涉相位的時(shí)空特性分離不同的分量,從而獲取高精度的形變速率.因此選取點(diǎn)的質(zhì)量直接影響最終的解算結(jié)果,選取高相干點(diǎn)的方法可以分為基于點(diǎn)和弧段兩類,前者包括歸一化振幅離差[12-13],記為NAD,以DA表示,信噪比[14],相干性[15],相位分析[16],極大似然法[17],相干矩陣特征值分析[18]等,這些方法通過(guò)依據(jù)不同原則分析散射體在時(shí)間序列上相位的穩(wěn)定性來(lái)選取高相干點(diǎn).后者首先進(jìn)行點(diǎn)分析得到較多的備選點(diǎn)集合[19-20],之后分析弧段質(zhì)量以及點(diǎn)之間的連通性得到高相干點(diǎn).選取高質(zhì)量的相干點(diǎn)并保證空間點(diǎn)密度之間保持平衡是MT-InSAR技術(shù)的一大難點(diǎn).寬松的閾值可以得到高密度的觀測(cè)結(jié)果,但是測(cè)量精度隨之降低,反之可以保證精度,但是觀測(cè)結(jié)果卻過(guò)于稀疏.因此需要研究針對(duì)不同質(zhì)量的相干點(diǎn)需采用不同的解算方式.目前的解決方案分為兩步,利用較小閾值選取較高質(zhì)量相干點(diǎn),解算結(jié)果作為一級(jí)參考網(wǎng),保證解算精度;之后擴(kuò)大閾值,構(gòu)成二級(jí)網(wǎng)求解一般質(zhì)量相干點(diǎn),以保證空間點(diǎn)密度[21-22].但是現(xiàn)有的空間點(diǎn)分級(jí)方法沒(méi)有考慮空間點(diǎn)密度的不均勻性,不同級(jí)相干點(diǎn)連接方式以及精度評(píng)價(jià)方法也有待完善,并且計(jì)算效率不高,使用原始分辨率數(shù)據(jù)的大范圍處理通常難以實(shí)施.

本文首先闡述了基于弧段選取的多級(jí)化MT-InSAR技術(shù)處理方法,包括基本觀測(cè)量,參考網(wǎng)和附加網(wǎng)的求解和精度評(píng)價(jià)方法,并對(duì)每一級(jí)網(wǎng)絡(luò)中相干點(diǎn)的選取以及求解策略作了詳細(xì)討論.利用C波段的ENVISAT、Sentinel-1A以及X波段的TerraSAR-X數(shù)據(jù)對(duì)上海磁懸浮及其沿線區(qū)域開(kāi)展沉降監(jiān)測(cè),分析不同數(shù)據(jù)的適用性,并與精密水準(zhǔn)數(shù)據(jù)驗(yàn)證,為城市大型結(jié)構(gòu)體健康監(jiān)測(cè)提供參考.

1 多級(jí)化MT-InSAR技術(shù)

1.1 基本觀測(cè)量

在MT-InSAR技術(shù)中,基本觀測(cè)量為雙差相位觀測(cè)值.假設(shè)在時(shí)間序列上有m幅已配準(zhǔn)的合成孔徑雷達(dá)(synthetic aperture radar,SAR)數(shù)據(jù),選取其中一幅作為主影像后可以得到m-1幅干涉相位圖.干涉圖中任意兩個(gè)像元i與j構(gòu)成弧段,去除平地相位后,第k幅干涉圖中弧段的相位差可以表示為[23]

(1)

式中:Δhi,j和Δvi,j分別為點(diǎn)i與j之間的高程改正之差和形變速率差;ni,j,k是相位模糊度;Bcik和tk為垂直基線和時(shí)間基線;Ri為像元點(diǎn)i與衛(wèi)星間的斜距;λ為雷達(dá)波長(zhǎng);θi為入射角;Ci,j為相位常數(shù),與主影像的大氣延遲有關(guān);ei,j,k為相位噪聲.對(duì)式(1)相位解纏后使用最小二乘估計(jì)即可得到參數(shù)解和對(duì)應(yīng)的參數(shù)精度[23].

在基于弧段選取的MT-InSAR技術(shù)中,初始相干點(diǎn)的選取包含大量的噪聲點(diǎn),最終可靠相干點(diǎn)的篩選依賴于選取高質(zhì)量弧段解.弧段質(zhì)量一般采用時(shí)序相干性來(lái)衡量[9].由于模型中采用了線性形變速率的假設(shè),因此加入線性速率的假設(shè)檢驗(yàn)進(jìn)一步篩選弧段.此二者構(gòu)成文獻(xiàn)中的弧段可接受條件[20].

基于本節(jié)得到的弧段解進(jìn)行多級(jí)化相干點(diǎn)構(gòu)網(wǎng)求解.由于城市區(qū)域的相干點(diǎn)之間相關(guān)性較低,因此采用歸一化振幅離差方法進(jìn)行相干點(diǎn)選取,該方法在選取相干點(diǎn)時(shí)僅考慮相干點(diǎn)的時(shí)間序列噪聲,并通過(guò)選取不同的閾值來(lái)控制選取相干點(diǎn)的質(zhì)量,選點(diǎn)的閾值越低,獲取點(diǎn)的精度越高,但是點(diǎn)的數(shù)目越少.數(shù)據(jù)處理的主要流程見(jiàn)圖1.

圖1 多級(jí)化相干點(diǎn)處理流程

Fig.1 Flowchart of processing coherent points with multi-level strategy

1.2 參考網(wǎng)

參考網(wǎng)的求解精度直接決定了最終結(jié)果的精度,構(gòu)成該網(wǎng)絡(luò)的相干點(diǎn)應(yīng)當(dāng)保證有足夠高的精度.因此采用歸一化振幅離差法選取相干點(diǎn)時(shí)使用盡可能低的閾值.考慮到大氣效應(yīng)的影響,求解弧段不宜過(guò)長(zhǎng)[13],參考網(wǎng)中相干點(diǎn)在空間內(nèi)應(yīng)盡可能均勻分布.依據(jù)這兩個(gè)準(zhǔn)則,參考網(wǎng)相干點(diǎn)選取的閾值使用0.25,是經(jīng)典永久散射體干涉測(cè)量方法中高相干點(diǎn)的選取閾值[12].這一步求解的目標(biāo)是將盡可能多的點(diǎn)連接起來(lái),采用二次構(gòu)網(wǎng)策略.首先使用Delaunay三角網(wǎng)連接所有點(diǎn),篩選弧段后,對(duì)孤立的點(diǎn)集采用漸進(jìn)式星形網(wǎng)構(gòu)成新的弧段.詳細(xì)處理步驟可參見(jiàn)文獻(xiàn)[20].

(2)

式中:B1為構(gòu)成參考網(wǎng)弧段相關(guān)的設(shè)計(jì)矩陣;Δv1為求解的所有弧段形變速率之差;(·)+是求解廣義逆運(yùn)算符;P1為弧段解的權(quán)陣,可采用如下定義:

(3)

參考網(wǎng)的形變速率估計(jì)精度計(jì)算方式為

(4)

式中:σ1為參考網(wǎng)的單位權(quán)中誤差.

當(dāng)測(cè)區(qū)范圍較大或是處理高分辨率數(shù)據(jù)時(shí),法方程系數(shù)矩陣階數(shù)較大,廣義逆求解緩慢或是由于矩陣過(guò)大而難以求解,因此需要減少構(gòu)成參考網(wǎng)的相干點(diǎn)數(shù)目.若使用降低閾值減少相干點(diǎn)的方式會(huì)導(dǎo)致相干點(diǎn)在空間分布不均勻,造成最終結(jié)果的精度在空間不均勻分布.此外,附加網(wǎng)在缺少參考網(wǎng)的相干點(diǎn)的區(qū)域無(wú)法生成弧段,從而造成相干點(diǎn)的遺漏.因此,本文采用基于點(diǎn)密度的分層抽樣方法來(lái)對(duì)空間點(diǎn)進(jìn)行稀疏,在保持相干點(diǎn)均勻分布的同時(shí),減少空間點(diǎn)密度較大區(qū)域的相干點(diǎn)數(shù)目.抽樣后得到的相干點(diǎn)組成參考網(wǎng),余下的點(diǎn)使用附加網(wǎng)進(jìn)行求解.

1.3 附加網(wǎng)

采用寬松閾值下的歸一化振幅離差方法選取相干點(diǎn),使用星形網(wǎng)將這些點(diǎn)與參考網(wǎng)中的點(diǎn)連接構(gòu)成新的弧段.構(gòu)網(wǎng)時(shí)仍然采用漸進(jìn)式策略,每次將待求點(diǎn)與若干個(gè)參考網(wǎng)中的點(diǎn)連接,循環(huán)處理直到該點(diǎn)與參考網(wǎng)之間有弧段連接,否則剔除該點(diǎn).聯(lián)合參考網(wǎng)中相干點(diǎn)的參數(shù)對(duì)新求解的弧段進(jìn)行積分可以得到附加網(wǎng)中相干點(diǎn)的形變速率:

(5)

式中:B2為構(gòu)成附加網(wǎng)弧段相關(guān)的設(shè)計(jì)矩陣;Δv2為求解的所有弧段相對(duì)形變速率之差;P2為權(quán)陣,與P1為相同定義方式;k為拉格朗日函數(shù)的系數(shù)矩陣;v2是附加網(wǎng)點(diǎn)待求解的形變速率;G為限制矩陣,定義如下:

Gv2=v1

(6)

附加網(wǎng)的形變速率的精度計(jì)算方式為

(7)

式中:σ2為附加網(wǎng)的單位權(quán)中誤差,Ql定義如下:

(8)

其中,B11和B22是B2的子矩陣,定義如下:

B2=[B11|B22]

(9)

求解參考網(wǎng)后,附加網(wǎng)可以使用任意多次以獲取更多相干點(diǎn)的信息.求解附加網(wǎng)時(shí),不同相干點(diǎn)的處理是相互獨(dú)立的.當(dāng)區(qū)域較大時(shí),可以將整個(gè)區(qū)域分為獨(dú)立的若干塊,對(duì)每塊區(qū)域內(nèi)的相干點(diǎn)依次處理,從而減少運(yùn)算的內(nèi)存開(kāi)銷,適用于大范圍的數(shù)據(jù)處理.

2 數(shù)據(jù)處理結(jié)果

本次實(shí)驗(yàn)采用2007~2018年間的31景C波段ENVISAT、53景X波段TerraSAR-X和51景C波段Sentinel-1A影像,詳細(xì)數(shù)據(jù)參數(shù)見(jiàn)表1.不同衛(wèi)星平臺(tái)數(shù)據(jù)均完全覆蓋磁懸浮沿線區(qū)域.由于上海地勢(shì)較為平坦,且樓房眾多,處理中不再使用地面高程數(shù)據(jù)進(jìn)行初始高程相位改正,僅改正平地相位.地面點(diǎn)高程引起的干涉相位在求解中通過(guò)時(shí)間序列方法進(jìn)行分離.

表1 SAR數(shù)據(jù)主要參數(shù)

以高分辨率TerraSAR-X數(shù)據(jù)為例,驗(yàn)證多級(jí)化數(shù)據(jù)方法的有效性.構(gòu)成參考網(wǎng)的相干點(diǎn)采用經(jīng)典的永久散射體[12]判斷方法進(jìn)行選取,DA閾值取為0.25,該閾值可以保證選出的相干點(diǎn)在整個(gè)空間均有分布.選取初始相干點(diǎn)數(shù)目為472 929個(gè),點(diǎn)分布參見(jiàn)圖2a,從圖中可以發(fā)現(xiàn)相干點(diǎn)在磁懸浮沿線呈現(xiàn)明顯的線性分布規(guī)律,磁懸浮周邊房屋區(qū)域也存在著高密度的相干點(diǎn).初始相干點(diǎn)數(shù)目多,難以直接求解,需要降低參考網(wǎng)相干點(diǎn)數(shù)目.

根據(jù)1.2節(jié)中的處理策略,首先使用降低閾值法,將初始選點(diǎn)的閾值更改為0.1,可以得到17 680個(gè)相干點(diǎn),點(diǎn)分布可參見(jiàn)圖2b.從該圖可以發(fā)現(xiàn),相干點(diǎn)在空間的分布很不均勻,磁懸浮沿線幾乎不存在相干點(diǎn).在多級(jí)化數(shù)據(jù)處理中,估算精度隨著分級(jí)的進(jìn)行而逐步遞減,若磁懸浮沿線的相干點(diǎn)只存在于附加網(wǎng)中,那么這些點(diǎn)的精度相對(duì)于整個(gè)測(cè)區(qū)是偏低的.因此降低閾值選取參考網(wǎng)的相干點(diǎn)是不合理的.

使用基于空間點(diǎn)密度的分層抽樣方法來(lái)合理選擇參考網(wǎng)的規(guī)模.基于圖2a的初始相干點(diǎn)集,先計(jì)算相干點(diǎn)的空間密度.取100 m為搜索半徑,對(duì)每個(gè)點(diǎn)統(tǒng)計(jì)搜索窗口內(nèi)的相干點(diǎn)數(shù)目,從而得到各點(diǎn)的空間相干點(diǎn)密度.為簡(jiǎn)化運(yùn)算,將空間點(diǎn)密度量化為等間距的5級(jí)密度,點(diǎn)密度分布見(jiàn)圖3.由圖3可以發(fā)現(xiàn),磁懸浮沿線相干點(diǎn)密度遠(yuǎn)低于周邊房屋區(qū)域.獲取空間點(diǎn)密度后進(jìn)行分層抽樣,依據(jù)空間點(diǎn)密度將空間點(diǎn)分為5層,每層抽樣時(shí)采用等間距隨機(jī)選點(diǎn),空間密度大的區(qū)域采用較大的采樣間隔,反之則采用較小的采樣間隔.在本次試驗(yàn)中,5層抽樣的采樣間隔比為1∶2∶3∶4∶5.設(shè)置目標(biāo)抽樣點(diǎn)數(shù)目為15 000,可以得到5層的采樣間隔依次是14,28,42,56,70.進(jìn)而從472 929個(gè)初始相干點(diǎn)中抽樣得到17 187個(gè)相干點(diǎn),這與降低閾值法得到的相干點(diǎn)數(shù)目相當(dāng),分層抽樣得到的相干點(diǎn)分布見(jiàn)圖2c.對(duì)比圖2b可以發(fā)現(xiàn),相干點(diǎn)在空間的分布更加均勻,并且較好地保留了磁懸浮沿線的相干點(diǎn).

a 初始選取相干點(diǎn)

b 降低選點(diǎn)閾值法

c 基于點(diǎn)密度的分層抽樣法

圖2 參考網(wǎng)相干點(diǎn)分布

Fig.2 Distribution of coherent points in the reference network

圖3 空間點(diǎn)五級(jí)密度分布

求得參考網(wǎng)中相干點(diǎn)的參數(shù)后,擴(kuò)大DA閾值至0.4,選取中等質(zhì)量相干點(diǎn),并利用附加網(wǎng)求解分層抽樣余下的相干點(diǎn)以及中等質(zhì)量的相干點(diǎn),為保證高質(zhì)量相干點(diǎn)的精度,將附加網(wǎng)分為兩次進(jìn)行,分層抽樣余下的相干點(diǎn)構(gòu)成一級(jí)附加網(wǎng),中等質(zhì)量相干點(diǎn)構(gòu)成二級(jí)附加網(wǎng),各級(jí)網(wǎng)獲取的相干點(diǎn)數(shù)目見(jiàn)表2.

表2 各級(jí)網(wǎng)絡(luò)獲取的相干點(diǎn)數(shù)目

Tab.2 The number of accepted coherent points in each network

參考網(wǎng)一級(jí)附加網(wǎng)二級(jí)附加網(wǎng)初選點(diǎn)數(shù)17 155470 1392 210 867接受點(diǎn)數(shù)15 813435 2341 711 388通過(guò)率/%92.292.677.4

基于弧段解的時(shí)間序列分析方法與選點(diǎn)策略無(wú)關(guān),在數(shù)據(jù)處理的流程中通過(guò)網(wǎng)絡(luò)連通性分析來(lái)去除噪聲點(diǎn),以保證最終結(jié)果的可靠性.初始選點(diǎn)質(zhì)量越高,經(jīng)網(wǎng)絡(luò)連通性分析后可接受的相干點(diǎn)越多.由表2可以發(fā)現(xiàn),參考網(wǎng)與一級(jí)附加網(wǎng)中的相干點(diǎn)是通過(guò)同一DA閾值選取的,因此有著相似的通過(guò)率,二級(jí)附加網(wǎng)的相干點(diǎn)由于選點(diǎn)閾值較為寬松,噪聲點(diǎn)多,因此可接受相干點(diǎn)的通過(guò)率低于前者.

Sentinel-1A與ENVISAT數(shù)據(jù)的處理使用類似的多級(jí)化處理方式.不同數(shù)據(jù)集得到的磁懸浮沿線及其周邊區(qū)域的視線向形變速率如圖4所示.為了盡可能突出沉降結(jié)果的細(xì)節(jié)信息,所有結(jié)果均保持SAR圖像的原始分辨率,未做多視處理.

不同衛(wèi)星數(shù)據(jù)在磁懸浮軌道及其周邊區(qū)域得到的相干點(diǎn)數(shù)目依次為:ENVISAT 129 182個(gè),Sentinel-1A 356 663個(gè),TerraSAR-X 2 165 435個(gè).由圖4b和4c可知,Sentinel-1A與TerraSAR-X得到的結(jié)果中,磁懸浮沿線均存在大量呈線性排列的相干點(diǎn),而在ENVISAT的結(jié)果中(圖4a),磁懸浮沿線相干點(diǎn)極少.Sentinel-1A與ENVISAT均為中等分辨率C波段數(shù)據(jù),且Sentinel-1A數(shù)據(jù)的方位向分辨率低于ENVISAT數(shù)據(jù).由此可以得知磁懸浮沿線上相干點(diǎn)的數(shù)目與距離向分辨率有關(guān),而與方位向分辨率關(guān)系較弱.ENVISAT數(shù)據(jù)由于距離向分辨率較低,在磁懸浮沿線難以形成有效的點(diǎn)目標(biāo),不能獲取磁懸浮軌道上足夠的形變信息.

對(duì)比表1中不同傳感器的測(cè)繪幅寬,盡管TerraSAR-X數(shù)據(jù)分辨率較高,但是單景數(shù)據(jù)覆蓋范圍僅為Sentinel-1A數(shù)據(jù)的1/30.ENVISAT數(shù)據(jù)與Sentinel-1A數(shù)據(jù)有著相似的空間分辨率,但是Sentinel-1A數(shù)據(jù)的覆蓋范圍是ENVISAT數(shù)據(jù)的5倍,且相干點(diǎn)的密度更高.因此Sentinel-1A數(shù)據(jù)的可用性更好.

此外,選取一塊形變顯著的區(qū)域(圖4中白色方框區(qū)域),詳細(xì)研究磁懸浮周邊區(qū)域沉降對(duì)軌道的影響.不同衛(wèi)星數(shù)據(jù)獲取的該區(qū)域沉降速率見(jiàn)圖5.圖5中表示沉降速率的色條與圖4一致,白線表示磁懸浮軌道位置.從圖5中可以發(fā)現(xiàn),在2007~2011年間,磁懸浮周邊區(qū)域沉降較大,最大年沉降速率大于20 mm·年-1(圖5a);2011~2015年間,該區(qū)域趨于穩(wěn)定,但是仍然存在沉降速率大于15 mm·年-1的點(diǎn)(圖5b);2016~2018年間,該區(qū)域大多數(shù)點(diǎn)保持穩(wěn)定,沉降速率小于10 mm·年-1(圖5c).不同數(shù)據(jù)得到的結(jié)果中沿軌道的線相干點(diǎn)密度有很大差異,

a ENVISAT

b TerraSAR-X

c Sentinel-1A

Fig.4 Subsidence velocity along the line of sight within the maglev zone

但是這些點(diǎn)在不同時(shí)段均保持穩(wěn)定,年平均沉降速率小于3 mm·年-1.因此周邊區(qū)域的沉降并未影響磁懸浮軌道的安全.

a ENVISAT

b TerraSAR-X

c Sentinel-1A

圖5 沉降較大區(qū)域放大圖

Fig.5 Magnified view of the selected area with large subsidence

3 結(jié)果驗(yàn)證

利用實(shí)測(cè)水準(zhǔn)數(shù)據(jù)驗(yàn)證本文多級(jí)MT-InSAR結(jié)果的可靠性.精密水準(zhǔn)包含2009~2014年間5期數(shù)據(jù),與TerraSAR-X數(shù)據(jù)獲取的時(shí)間段一致.該區(qū)域共26個(gè)水準(zhǔn)點(diǎn)(BM1~BM26),位于磁懸浮軌道以及周邊區(qū)域,具體位置見(jiàn)圖6.

圖6 磁懸浮區(qū)域水準(zhǔn)點(diǎn)分布

Fig.6 Distribution of levelling points within the maglev zone

水準(zhǔn)測(cè)量基于地面單個(gè)點(diǎn),而SAR圖像中的單個(gè)像元反映了地面某一塊區(qū)域的綜合信息,二者不能直接對(duì)比.為了獲取MT-InSAR結(jié)果中與水準(zhǔn)點(diǎn)對(duì)應(yīng)的沉降信息,以每個(gè)水準(zhǔn)點(diǎn)為中心,選取30 m×30 m的搜索窗口,以搜索窗口內(nèi)所有相干點(diǎn)的加權(quán)平均值作為與水準(zhǔn)點(diǎn)對(duì)應(yīng)的結(jié)果.MT-InSAR結(jié)果表明磁懸浮軌道與周邊區(qū)域呈現(xiàn)不同的沉降信息,在使用搜索窗口時(shí)必須將軌道與周邊區(qū)域相干點(diǎn)分開(kāi)處理.因此水準(zhǔn)數(shù)據(jù)的驗(yàn)證分為兩組,磁懸浮軌道上點(diǎn)以及軌道周邊區(qū)域的點(diǎn).

使用緩沖區(qū)法分離磁懸浮軌道上點(diǎn),選取10 m緩沖區(qū)半徑,在磁懸浮沿線建立緩沖區(qū),位于緩沖區(qū)內(nèi)部的點(diǎn)視為磁懸浮軌道相干點(diǎn).圖7表示磁懸浮軌道上水準(zhǔn)點(diǎn)與MT-InSAR相干點(diǎn)的形變速率對(duì)比結(jié)果,二者結(jié)果差異的平均值是-0.020 4 mm·年-1,標(biāo)準(zhǔn)差是0.628 3 mm·年-1.圖8表示磁懸浮周邊區(qū)域內(nèi)水準(zhǔn)點(diǎn)與MT-InSAR相干點(diǎn)的形變速率對(duì)比結(jié)果,二者差異的均值為-0.414 4 mm·年-1,標(biāo)準(zhǔn)差是1.170 8 mm·年-1.MT-InSAR的結(jié)果中,磁懸浮軌道上點(diǎn)與水準(zhǔn)測(cè)量結(jié)果符合度較高.而磁懸浮周邊區(qū)域中相干點(diǎn)難以與水準(zhǔn)點(diǎn)直接對(duì)應(yīng),因此符合度較低.

圖7 磁懸浮軌道上水準(zhǔn)點(diǎn)與MT-InSAR形變速率對(duì)比

Fig.7 Comparison between deformation velocities derived from levelling and from MT-InSAR along the maglev track

圖8 磁懸浮周邊區(qū)域水準(zhǔn)點(diǎn)與MT-InSAR形變速率對(duì)比

Fig.8 Comparison of deformation velocity between levelling and MT-InSAR nearby the maglev track

選取4個(gè)水準(zhǔn)點(diǎn)BM6、7、11、12,對(duì)比水準(zhǔn)與MT-InSAR獲取的時(shí)間序列累積形變量,如圖9所示.由圖可以發(fā)現(xiàn),MT-InSAR得到的形變時(shí)間序列與水準(zhǔn)結(jié)果有著良好的一致性.并且水準(zhǔn)結(jié)果表明BM11與BM6在2013年10月至2014年期間出現(xiàn)了明顯的抬升現(xiàn)象,這與MT-InSAR獲取的結(jié)果有著相似的形變趨勢(shì).2013年開(kāi)始出現(xiàn)的抬升現(xiàn)象與上海市實(shí)行的限制地下水開(kāi)采以及地下水回灌有關(guān)[24].

a 水準(zhǔn)點(diǎn)11

b 水準(zhǔn)點(diǎn)6

c 水準(zhǔn)點(diǎn)12

d 水準(zhǔn)點(diǎn)7

圖9 水準(zhǔn)點(diǎn)與MT-InSAR時(shí)間序列形變量對(duì)比

Fig.9 Comparison between time series displacement derived from levelling and from MT-InSAR

4 結(jié)語(yǔ)

本文介紹了一種多級(jí)化的MT-InSAR數(shù)據(jù)處理方法,提出了一種基于相干點(diǎn)空間密度的抽樣方法,綜合考慮選取相干點(diǎn)的密度以及空間分布,得到優(yōu)化的分級(jí)MT-InSAR處理方法.并利用該方法獲取了2007~2018年間的上海磁懸浮及其周邊區(qū)域沉降信息,不同數(shù)據(jù)結(jié)果均顯示磁懸浮軌道保持穩(wěn)定,而周邊區(qū)域存在不同程度的沉降,隨著時(shí)間的累積,沉降速率是逐漸變小的.采用精密水準(zhǔn)數(shù)據(jù)對(duì)沉降監(jiān)測(cè)結(jié)果進(jìn)行了精度驗(yàn)證,兩者在時(shí)間以及空間上吻合度較高.MT-InSAR技術(shù)得到的形變速率可以滿足高精度的地面形變監(jiān)測(cè)要求.

此外,本文還分析了不同SAR影像數(shù)據(jù)分辨率對(duì)于磁懸浮沿線相干點(diǎn)獲取的影響,發(fā)現(xiàn)相干點(diǎn)的密度與SAR數(shù)據(jù)的距離向分辨率有著密切關(guān)系,而方位向分辨率影響較小.對(duì)比上一代C波段中等分辨率ENVISAT數(shù)據(jù),開(kāi)源的Sentinel-1A數(shù)據(jù)有著更大的覆蓋范圍,并且憑借較高的距離向分辨率,可以獲取更多的相干點(diǎn),在城市地表以及大型結(jié)構(gòu)體形變檢測(cè)中有著突出的優(yōu)勢(shì).