胡 程 鄧云開(kāi) 田衛(wèi)明 曾 濤

①(北京理工大學(xué)信息與電子學(xué)院雷達(dá)技術(shù)研究所 北京 100081)

②(北京理工大學(xué)衛(wèi)星導(dǎo)航電子信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室 北京 100081)

1 引言

作為一種高精度的形變測(cè)量?jī)x器，地基干涉合成孔徑雷達(dá)(Ground-Based Interferometric Synthetic Aperture Radar,GB-InSAR)已經(jīng)在形變監(jiān)測(cè)領(lǐng)域得到了廣泛的應(yīng)用[1]。GB-InSAR通常是基于差分干涉測(cè)量技術(shù)，通過(guò)對(duì)同一位置、不同時(shí)刻獲取的兩幅SAR圖像進(jìn)行差分干涉處理，基于相位信息來(lái)實(shí)現(xiàn)形變測(cè)量，其一般工作在X或者Ku波段，形變測(cè)量精度可以達(dá)到毫米或者亞毫米量級(jí)。GB-InSAR測(cè)量誤差的主要來(lái)源是大氣相位，由于不同時(shí)刻氣象條件的不同，電磁波在大氣中傳播的速度會(huì)發(fā)生改變，從而導(dǎo)致大氣延時(shí)誤差[2]。

為實(shí)現(xiàn)大氣相位的補(bǔ)償，國(guó)內(nèi)外學(xué)者主要提出了3種解決方法。第1種的補(bǔ)償方法是在觀測(cè)場(chǎng)景內(nèi)建立氣象站，基于大氣折射率模型，利用氣象數(shù)據(jù)(溫度、濕度、大氣壓)來(lái)對(duì)大氣相位進(jìn)行定量估計(jì)[3]；第2種是在場(chǎng)景中人共布設(shè)或者選擇出一些高度穩(wěn)定的特征點(diǎn)，采用空間插值的方法實(shí)現(xiàn)對(duì)整幅圖像的大氣相位的補(bǔ)償[4,5]；第3種則是基于永久散射體(Permanent Scatterer,PS)技術(shù)，根據(jù)大氣相位的空間分布特征，建立描述大氣相位的方程，估計(jì)大氣相位參數(shù)，實(shí)現(xiàn)大氣相位的補(bǔ)償[6]。

在基于PS技術(shù)進(jìn)行大氣相位補(bǔ)償時(shí)，首先需要選擇出未處在形變區(qū)域的PS點(diǎn)，然后對(duì)大氣相位進(jìn)行合理的建模，并建立線性方程組，實(shí)現(xiàn)對(duì)大氣參數(shù)的粗估計(jì)，此后剔除與模型偏差較大的PS點(diǎn)，逐步迭代實(shí)現(xiàn)大氣參數(shù)的準(zhǔn)確估計(jì)。該方法不需要?dú)庀髤?shù)及布設(shè)特征點(diǎn)，可以基于大量的PS點(diǎn)進(jìn)行大氣參數(shù)的估計(jì)，估計(jì)精度較高。一般情況下，大氣在空間上均勻變化，可以將大氣相位建模為隨斜距線性變化的分量[7]。在地形陡峭的山區(qū)，大氣在空間上非均勻變化，可以采用兼顧斜距和高程的多參數(shù)模型[8]。

基于PS技術(shù)的大氣相位補(bǔ)償方法，已經(jīng)在GB-InSAR領(lǐng)域得到了廣泛的應(yīng)用，但還存在一些典型的問(wèn)題。首先，該方法要求采用場(chǎng)景中的非形變PS點(diǎn)建立觀測(cè)方程組；其次，氣象條件一直在隨時(shí)間改變，大氣相位的時(shí)變性很強(qiáng)，在較差天氣條件(降雨、降雪、強(qiáng)風(fēng)等)下，大氣在空間上非均勻變化，導(dǎo)致大氣相位可能表現(xiàn)出復(fù)雜的空變性，無(wú)法建立合理的多參數(shù)模型來(lái)模擬大氣相位。因此，常規(guī)的基于PS技術(shù)的補(bǔ)償方法，在應(yīng)用于時(shí)序GB-InSAR圖像處理時(shí)，還存在較大的改進(jìn)空間。

針對(duì)上述問(wèn)題，本文借鑒在場(chǎng)景中布設(shè)若干個(gè)穩(wěn)定的地面控制點(diǎn)，通過(guò)空間維插值估計(jì)大氣相位的思想。首先采用常規(guī)補(bǔ)償方案對(duì)所有的干涉相位圖進(jìn)行大氣相位補(bǔ)償，并分析PS點(diǎn)的相位序列的標(biāo)準(zhǔn)差，設(shè)定門限選擇出穩(wěn)定PS點(diǎn)；然后基于K均值聚類算法 (K-means clustering algorithm,K-means)將穩(wěn)定PS點(diǎn)劃分出一定數(shù)量的子區(qū)域，將每一個(gè)子區(qū)域的中心點(diǎn)設(shè)定為控制點(diǎn)，采用反距離加權(quán)插值算法估計(jì)所有PS點(diǎn)的大氣相位，從而實(shí)現(xiàn)非線性大氣相位的補(bǔ)償。

本文詳細(xì)介紹了補(bǔ)償方法的實(shí)現(xiàn)流程，主要分為穩(wěn)定PS點(diǎn)選擇和空間維插值補(bǔ)償兩部分，然后分別采用改進(jìn)方法和常規(guī)方法對(duì)460幅地基多輸入多輸出 (Multiple-Input Multiple-Output,MIMO)雷達(dá)圖像進(jìn)行了分析，對(duì)比驗(yàn)證了本文方法的有效性。

2 補(bǔ)償方法

2.1 穩(wěn)定PS點(diǎn)選擇

在利用像素點(diǎn)的相位信息進(jìn)行形變測(cè)量時(shí)，差分干涉相位的質(zhì)量直接影響到形變測(cè)量的精度。因此，GB-InSAR差分干涉處理時(shí)，通常需要選擇出一些高質(zhì)量的像素點(diǎn)，即PS點(diǎn)，來(lái)進(jìn)行形變分析。在GB-InSAR領(lǐng)域，廣泛采用幅度離差法來(lái)進(jìn)行PS點(diǎn)的選擇，一個(gè)像素點(diǎn)的時(shí)序幅度序列的標(biāo)準(zhǔn)差與均值之比，即為幅度離差指數(shù)。通過(guò)設(shè)置合理的幅度離差門限，即可以實(shí)現(xiàn)PS點(diǎn)的選擇[9]。

一個(gè)PS點(diǎn)的差分干涉相位 Δφ可以建模為

其中，φdefo為形變相位分量；φatm為兩幅圖像獲取期間，由大氣條件改變所導(dǎo)致的大氣相位分量；φnoi為噪聲相位分量；2kπ為相位模糊度，k為整數(shù)。由于相位周期性的影響，Δφ處在-π～π的范圍內(nèi)[10]。在進(jìn)行大氣相位補(bǔ)償前，要對(duì)干涉相位圖進(jìn)行相位解纏，可以采用非均勻網(wǎng)格下的最小費(fèi)用流算法，下文中采用 Δφ指代解纏相位。

大氣相位φatm可以建模為

其中，λ為信號(hào)的波長(zhǎng)，ΔN表示折射率的變化，其隨著空間r和時(shí)間t發(fā)生變化，L表示信號(hào)的傳輸路徑。一般情況下，大氣的空間同質(zhì)性很好，可以假設(shè)ΔN在空間r上不發(fā)生變化，在時(shí)間t上隨機(jī)變化，因此可以將φatm建模為隨斜距線性變化的分量

其中，β0為常數(shù)分量，β1表示與斜距相關(guān)的線性系數(shù)，R表示雷達(dá)與目標(biāo)點(diǎn)之間的距離?；赑S技術(shù)進(jìn)行大氣相位補(bǔ)償時(shí)，首先建立線性方程組

ΔΦ為n個(gè)PS點(diǎn)的解纏相位Δφ1,Δφ2,···,Δφn構(gòu)成的n×1維向量，X為常數(shù)1與n個(gè)PS點(diǎn)的斜距R1,R2,···,Rn構(gòu)成的n×2維矩陣，β為待估計(jì)的2×1維向量，ε為n×1維的隨機(jī)誤差向量，ε1,ε2,···,εn為各PS點(diǎn)的模型誤差相位分量。采用最小二乘法對(duì)β進(jìn)行估計(jì)，可以得到

其中，T表示矩陣的轉(zhuǎn)置。大氣相位的估計(jì)分量為

ΔΦ和ΔΦA(chǔ)PS之間的差值即為補(bǔ)償后相位。首先基于所有的PS點(diǎn)對(duì)β進(jìn)行估計(jì)，然后為提高估計(jì)的精度，剔除一些不可靠PS點(diǎn)。不滿足式(7)準(zhǔn)則的PS點(diǎn)即為不可靠PS點(diǎn)，ΔTatm的取值在0.1～0.2 rad范圍內(nèi)。剔除不可靠PS點(diǎn)后，基于剩余的PS點(diǎn)對(duì)大氣相位參數(shù)進(jìn)行二次估計(jì)[11]

除了最基本的線性斜距模型外，常用的參數(shù)模型還包括高階斜距模型、斜距-方位角模型、斜距-高程模型等，如式(8)-式(10)。其中β0,β1和β2是各模型中待估計(jì)的未知參數(shù)，R,θ和h分別代表目標(biāo)點(diǎn)的斜距、方位角和高程。對(duì)線性方程組式(4)修正后，即可以進(jìn)行大氣相位補(bǔ)償

基于PS技術(shù)的補(bǔ)償方法，其有效補(bǔ)償?shù)臈l件，首先是要求所采用的參數(shù)模型可以準(zhǔn)確模擬大氣相位，其次是迭代估計(jì)參數(shù)時(shí)，可以有效剔除不可靠PS點(diǎn)。采用上述方法對(duì)一幅干涉相位圖進(jìn)行大氣相位補(bǔ)償后，一個(gè)PS點(diǎn)的補(bǔ)償后相位ΔφAPC可以表示為

2.2 空間維插值補(bǔ)償

2.2.1 子區(qū)域劃分

在選擇出穩(wěn)定PS點(diǎn)后，對(duì)于每一幅干涉相位圖，分別進(jìn)行大氣相位補(bǔ)償。由于穩(wěn)定PS點(diǎn)的相位分量中包含噪聲相位和大氣相位，考慮到噪聲相位在干涉相位圖上隨機(jī)變化，不具備空間相關(guān)性，而大氣相位雖然會(huì)在整幅圖像范圍內(nèi)發(fā)生變化，但在較小的距離范圍內(nèi)可以視為一個(gè)常數(shù)。因此如果對(duì)空間上較小距離范圍內(nèi)的所有PS點(diǎn)進(jìn)行相位平均，則噪聲相位可以得到很好的濾除。由于穩(wěn)定PS點(diǎn)是非均勻的分布在整幅圖像范圍內(nèi)，可以基于K-means算法對(duì)穩(wěn)定PS點(diǎn)進(jìn)行簇劃分，將每一個(gè)簇定義為一個(gè)子區(qū)域。每一個(gè)簇的中心點(diǎn)作為控制點(diǎn)，并對(duì)該簇內(nèi)的所有PS點(diǎn)進(jìn)行相位平均，作為當(dāng)前控制點(diǎn)的相位。

K-means算法的實(shí)現(xiàn)原理是對(duì)于給定的樣本集x，按照樣本之間的距離大小，將樣本集劃分為K個(gè)簇，讓各個(gè)簇內(nèi)的點(diǎn)的距離盡可能的小，而讓各個(gè)簇之間的點(diǎn)的距離盡可能的大[12]。假設(shè)將樣本集x劃分為K個(gè)簇(C1,C2,···,CK)，各簇間的平方誤差和E可以表示為

其中，∥·∥表示2階范數(shù)，即向量的模。μi是簇Ci的均值向量，可以表示為

其中，|·|表示1階范數(shù)，即簇中點(diǎn)的數(shù)量。

經(jīng)過(guò)K-means劃分后可以得到K個(gè)子區(qū)域，將每一個(gè)子區(qū)域中心點(diǎn)定義為控制點(diǎn)CPKM，其相位為該子區(qū)域內(nèi)所有PS點(diǎn)的相位均值。

2.2.2 大氣相位估計(jì)

基于這K個(gè)控制點(diǎn)CPKM的平均相位，進(jìn)行空間維插值來(lái)估計(jì)所有PS點(diǎn)的大氣相位，可以采用反距離加權(quán)插值算法。反距離加權(quán)插值算法是利用已知點(diǎn)與待插點(diǎn)之間的距離來(lái)定義加權(quán)因子，然后加權(quán)計(jì)算待插點(diǎn)的相位，距離越近加權(quán)比重越大。其計(jì)算公式為

其中，Z(x,y)為插值結(jié)果，(x,y)為待插點(diǎn)的空間坐標(biāo)，Zi為第i(i=1,2,···,n)個(gè)參考點(diǎn)的數(shù)值，|di|表示待插點(diǎn)與第i個(gè)參考點(diǎn)之間的空間距離，μ表示加權(quán)因子的冪指數(shù)，一般取為2[13]。

將這些控制點(diǎn)視為一組離散點(diǎn)，可以基于Delaunay三角剖分法則來(lái)構(gòu)建一個(gè)三角形網(wǎng)絡(luò)。對(duì)于該網(wǎng)絡(luò)內(nèi)的任意一個(gè)三角形，其由3個(gè)控制點(diǎn)作為頂點(diǎn)，且外接圓中不包含其他控制點(diǎn)。在基于反距離加權(quán)插值算法來(lái)估計(jì)所有PS點(diǎn)的大氣相位時(shí)，如果一個(gè)PS點(diǎn)處在某一個(gè)三角形內(nèi)部，則以該三角形的3個(gè)控制點(diǎn)作為參考點(diǎn)，基于式(15)估計(jì)當(dāng)前PS點(diǎn)的大氣相位。如果一個(gè)PS點(diǎn)處在所有三角形的外部，則選擇最近的3個(gè)控制點(diǎn)，同樣基于式(15)來(lái)估計(jì)大氣相位。由于每一個(gè)控制點(diǎn)的相位，均是通過(guò)對(duì)一個(gè)子區(qū)域內(nèi)的所有PS點(diǎn)進(jìn)行相位平均獲得，則上述方法可以視為采用了大量的PS點(diǎn)作為參考點(diǎn)，空間維插值精度可以得到保證。

3 實(shí)驗(yàn)信息

本實(shí)驗(yàn)選定的區(qū)域?yàn)橐宦短扉_(kāi)采礦坑(E118°36′23′′,N40°06′44′′)，其位于河北省遷安市馬蘭莊鎮(zhèn)。圖1(a)所示為場(chǎng)景照片，其中紅色橢圓所示為場(chǎng)景的形變區(qū)域。圖1(b)所示為場(chǎng)景衛(wèi)星圖，礦坑正上方為橢圓形，其長(zhǎng)軸約1050 m，短軸約680 m，其中黃色矩形代表雷達(dá)的布放位置，雷達(dá)的觀測(cè)視角范圍為 60°。礦坑邊坡為典型巖質(zhì)邊坡，基本無(wú)植被覆蓋，最大開(kāi)采深度約400 m，邊坡傾角為38°～47°[14,15]。

實(shí)驗(yàn)中采用一部MIMO雷達(dá)對(duì)該礦坑進(jìn)行了監(jiān)測(cè)，連續(xù)獲取了460幅雷達(dá)圖像，時(shí)間從2018年10月26日17:30～2018年10月27日11:30。該MIMO雷達(dá)采用16個(gè)發(fā)射天線構(gòu)成兩個(gè)密集子陣列，16個(gè)接收天線構(gòu)成一個(gè)稀疏子陣列，如圖2所示。其工作在Ku波段，波長(zhǎng)為λ=1.86 cm，等效合成孔徑長(zhǎng)度為1.138 m，角分辨率為0.466°[16]。

圖3(a)所示為該露天礦坑在極坐標(biāo)系下的成像結(jié)果，邊坡區(qū)域內(nèi)像素點(diǎn)的幅值主要分布在-30～0 dB的范圍內(nèi)。基于幅度離差法進(jìn)行PS點(diǎn)的選擇，設(shè)置幅度離差門限0.15，幅度門限-25 dB，可篩選出61764個(gè)PS點(diǎn)，如圖3(b)所示。可以看出，坡體上大部分像素點(diǎn)的幅度穩(wěn)定性很高，僅中間的道路上由于雷達(dá)觀測(cè)視角的原因，未能有效選擇出PS點(diǎn)。

圖1 場(chǎng)景信息Fig.1 Scene information

圖2 MIMO雷達(dá)系統(tǒng)照片F(xiàn)ig.2 Photo of the MIMO radar system

以第1幅圖像為主圖像，最后一幅圖像為輔圖像，獲取差分干涉相位圖，如圖4(a)所示，將該相位圖反投到3維地形圖上，如圖4(b)所示。借助于3維地形，可以更加準(zhǔn)確地確定形變區(qū)域的發(fā)生位置。圖4(c)所示為PS點(diǎn)的干涉相位圖，顯然大部分PS點(diǎn)的相位在0 rad左右，而紅色橢圓標(biāo)識(shí)出的一塊區(qū)域，其干涉相位與其他區(qū)域的明顯不同，該區(qū)域?yàn)榘l(fā)生了較大形變的形變區(qū)域。圖4(d)所示為PS點(diǎn)的干涉相位隨其斜距變化的分布圖，如果將大氣相位建模為隨斜距線性變化的相位分量，其中紅色實(shí)線所示為線性大氣相位估計(jì)結(jié)果，顯然存在很大的誤差。因此，對(duì)較長(zhǎng)時(shí)間基線的干涉相位圖進(jìn)行大氣相位補(bǔ)償時(shí)，如果不能合理的剔除形變PS點(diǎn)，將會(huì)嚴(yán)重影響到對(duì)大氣參數(shù)的估計(jì)。在處理時(shí)序GB-InSAR圖像時(shí)，為提高大氣相位補(bǔ)償?shù)臏?zhǔn)確度，進(jìn)行差分干涉的兩幅雷達(dá)圖像之間的時(shí)間基線不宜過(guò)大。

圖3 MIMO雷達(dá)圖像與PS圖Fig.3 MIMO radar image and PS map

圖4 長(zhǎng)時(shí)間基線干涉相位圖Fig.4 Interferometric phase map with long temporal baseline

4 實(shí)驗(yàn)結(jié)果

分析這460幅時(shí)序MIMO雷達(dá)圖像時(shí)，采用相鄰的兩幅圖像構(gòu)成一個(gè)干涉圖像對(duì)，則可以獲取459幅干涉相位圖。每一幅干涉相位圖的時(shí)間基線僅為2～3 min，通常在這么短的時(shí)間內(nèi)，大氣條件的變化很小，相應(yīng)的大氣相位誤差分量也很小。圖5(a)和圖5(b)所示為干涉相位圖A和其相位散點(diǎn)圖，所有PS點(diǎn)的干涉相位均在0 rad左右。在采用線性斜距模型補(bǔ)償大氣相位時(shí)，最大的補(bǔ)償分量不超過(guò)0.05 rad。但在部分時(shí)間段，大氣在空間上不再是均勻變化，導(dǎo)致大氣相位呈現(xiàn)出非線性變化。圖5(c)和圖5(d)所示為干涉相位圖B和其相位散點(diǎn)圖，顯然PS點(diǎn)的干涉相位變化比較復(fù)雜，采用線性斜距模型補(bǔ)償大氣相位時(shí)，最大會(huì)帶來(lái)約0.5 rad的補(bǔ)償誤差，相應(yīng)的形變測(cè)量誤差約為0.74 mm，會(huì)嚴(yán)重影響到MIMO雷達(dá)的形變測(cè)量精度。即使是采用其他模型對(duì)相位圖B進(jìn)行分析，依然無(wú)法有效的補(bǔ)償大氣相位。

對(duì)459幅干涉相位圖進(jìn)行相位解纏，并采用線性斜距模型進(jìn)行大氣相位補(bǔ)償，然后計(jì)算每一個(gè)PS點(diǎn)的補(bǔ)償后相位序列的標(biāo)準(zhǔn)差。圖6(a)所示為所有PS點(diǎn)的相位標(biāo)準(zhǔn)差圖，顯然形變PS點(diǎn)的標(biāo)準(zhǔn)差很大，設(shè)置0.3 rad的標(biāo)準(zhǔn)差門限，篩選出4341個(gè)PS點(diǎn)。這些PS點(diǎn)的分布情況圖如圖6(b)所示，除紅色橢圓形變區(qū)域內(nèi)的形變PS點(diǎn)外，還有少量的形變區(qū)域外的噪聲PS點(diǎn)也被選擇，如紫色橢圓所標(biāo)識(shí)出的部分點(diǎn)。

經(jīng)過(guò)相位標(biāo)準(zhǔn)差篩選后，得到57423個(gè)穩(wěn)定PS點(diǎn)，設(shè)定每個(gè)子區(qū)域中的PS點(diǎn)平均數(shù)量為200個(gè)，則可以劃分出287個(gè)子區(qū)域。以圖5(c)所示的干涉相位圖B為例來(lái)說(shuō)明本文所提方法，圖7(a)中紅色方形點(diǎn)標(biāo)識(shí)出了每一個(gè)子區(qū)域中心，即287個(gè)控制點(diǎn)，以干涉相位圖B作為背景，圖7(b)為局部放大圖，其中每一個(gè)黑色多邊形，代表各個(gè)子區(qū)域中PS點(diǎn)的最小外接多邊形。之后對(duì)于各個(gè)子區(qū)域中的PS點(diǎn)進(jìn)行相位平均，獲取控制點(diǎn)的平均相位，采用反距離加權(quán)法估計(jì)所有PS點(diǎn)的大氣相位，如圖7(c)所示，可以很直觀地看出，估計(jì)圖中PS點(diǎn)的相位隨其斜距呈現(xiàn)出非線性變化，和圖5(c)高度相似。圖7(d)所示則為補(bǔ)償結(jié)果，所有PS點(diǎn)的補(bǔ)償后相位均在0 rad左右，空變性的大氣相位得到了有效的補(bǔ)償。

圖5 短時(shí)間基線干涉相位圖A和BFig.5 Interferometric phase maps of A and B with short temporal baselines

圖6 形變PS點(diǎn)選擇Fig.6 Selection of deformation PS

圖7 非線性大氣相位補(bǔ)償Fig.7 Non-linear atmospheric phase compensation

采用本文所提非線性補(bǔ)償方法和常規(guī)的線性斜距模型補(bǔ)償方法，對(duì)這459幅相位圖分別分析，然后從不同斜距處選出5個(gè)幅度離差指數(shù)最小的PS點(diǎn)作為參考點(diǎn)，來(lái)對(duì)比說(shuō)明本文方法的有效性。圖8(a)所示為采用本文方法獲取的累積相位圖，由459幅補(bǔ)償后的干涉相位圖累加獲取，該幅圖像可以用來(lái)反映整個(gè)監(jiān)測(cè)周期內(nèi)，場(chǎng)景中PS點(diǎn)的相位變化情況。形變區(qū)域PS點(diǎn)的相位變化達(dá)到了約-3 rad，非形變區(qū)域PS點(diǎn)的相位變化在±1 rad范圍內(nèi)，直觀上很難看出非形變區(qū)域PS點(diǎn)的相位是否存在空變性。圖8(a)中紅色方形點(diǎn)標(biāo)識(shí)出的點(diǎn)1～點(diǎn)5為非形變區(qū)域的參考點(diǎn)。圖8(b)所示為這5個(gè)參考點(diǎn)的時(shí)序相位變化曲線，點(diǎn)1～點(diǎn)5的相位變化趨勢(shì)高度相似，這是由于大氣相位具有較高的空間相關(guān)性。圖8(c)和圖8(b)所示分別為采用本文方法和常規(guī)方法對(duì)點(diǎn)1～點(diǎn)5的補(bǔ)償結(jié)果。常規(guī)方法的補(bǔ)償結(jié)果中，點(diǎn)1的相位隨著時(shí)間逐漸增大，這說(shuō)明了大氣相位殘余誤差隨著時(shí)間逐漸積累了起來(lái)。本文方法的補(bǔ)償結(jié)果中，各點(diǎn)的相位序列變化更為隨機(jī)。

圖9(a)和圖9(b)所示為采用常規(guī)方法和本文方法獲取的累積相位圖，為了更直觀地對(duì)比場(chǎng)景中PS點(diǎn)的相位變化情況，兩幅圖像中均忽略了形變區(qū)域?？梢院苤庇^地看出，常規(guī)方法的補(bǔ)償相位圖中，PS點(diǎn)的相位依然呈現(xiàn)出明顯的空變性，尤其是圖像中500～600 m范圍內(nèi)的區(qū)域，該區(qū)域PS點(diǎn)的相位達(dá)到了1 rad，相應(yīng)的形變量為1.48 mm。由于大氣相位補(bǔ)償誤差隨著時(shí)間積累起來(lái)的原因，這部分區(qū)域很容易被誤認(rèn)為在雷達(dá)監(jiān)測(cè)周期內(nèi)，也出現(xiàn)了形變。

圖9 累積相位圖Fig.9 Cumulative phase maps

5 結(jié)束語(yǔ)

本文提出了一種適用于時(shí)序GB-InSAR圖像非線性大氣相位補(bǔ)償?shù)母倪M(jìn)方法，解決了常規(guī)方法中多參數(shù)模型無(wú)法準(zhǔn)確模擬大氣相位，從而無(wú)法對(duì)部分干涉相位圖進(jìn)行有效補(bǔ)償?shù)膯?wèn)題。采用本文方法與常規(guī)方法，對(duì)一露天礦坑的460幅時(shí)序地基MIMO雷達(dá)圖像分別進(jìn)行了處理，對(duì)比了參考點(diǎn)的時(shí)序相位曲線和累積相位圖，驗(yàn)證了本文方法的有效性。

本文方法還存在著一些不足。首先本文方法要求先采用常規(guī)方法對(duì)時(shí)序干涉相位圖進(jìn)行補(bǔ)償，基于PS點(diǎn)的補(bǔ)償后相位序列的標(biāo)準(zhǔn)差大小來(lái)選擇穩(wěn)定PS點(diǎn)，但相位標(biāo)準(zhǔn)差門限是人為設(shè)定的；其次如果部分PS點(diǎn)僅在少量圖像中發(fā)生了形變，這些PS點(diǎn)的相位標(biāo)準(zhǔn)差可能較小，無(wú)法有效地篩選出；最后本文方法目前是應(yīng)用于對(duì)時(shí)序GB-InSAR圖像的事后大氣相位補(bǔ)償，實(shí)際形變監(jiān)測(cè)中更關(guān)注的是實(shí)時(shí)處理，需要進(jìn)一步研究實(shí)時(shí)形變處理中的非線性大氣相位補(bǔ)償方法。