地磅 地基雷達(dá)測量誤差源及提高精度的措施
為了提高地基雷達(dá)系統(tǒng)的監(jiān)測精度,提出一種新的高精度的地基合成孔徑雷達(dá)干涉 GBSAR( Ground Based SAR) 監(jiān)測技術(shù)。融合了 GBSAR 關(guān)鍵技術(shù)以及數(shù)據(jù)處理理論,對影響測量精度的誤差項進(jìn)行了分析研究并從 3 個方面( 雷達(dá)系統(tǒng)、數(shù)據(jù)獲取、數(shù)據(jù)處理) 具體給出了相應(yīng)的精度提高措施。基于變形監(jiān)測系統(tǒng)( IBIS-S) 對 GBSAR 理論研究進(jìn)行實驗驗證。實驗表明 GBSAR 技術(shù)對于目標(biāo)體的監(jiān)測精度較高。最終得到 GBSAR 技術(shù)產(chǎn)生的誤差的途徑主要源自數(shù)據(jù)采集的過程。
干涉合成孔徑雷達(dá) ( Interferometric Synthetic( InSAR) 技術(shù)是一種基于獲取到多幅 SAR 復(fù)圖像的相位特性信息,全天候、高精度的對某一特定地區(qū)進(jìn)行較大范圍的地理地貌以及地表運(yùn)動變化發(fā)生形變的規(guī)模進(jìn)行實時監(jiān)測,因此,InSAR 成為了地形測繪以及地表形變監(jiān)控測量的一項新技術(shù) 。InSAR技術(shù)無論實在對地形形變數(shù)據(jù)采樣頻率以及數(shù)據(jù)采樣的密度方面的精確程度,都遠(yuǎn)遠(yuǎn)高于 GPS 系統(tǒng),因此,InSAR 技術(shù)適合對地形起伏波動較大的山區(qū)
進(jìn)行大面積多點長期的監(jiān)控測量 |
外學(xué)者基于 InSAR 技術(shù)理論,進(jìn)一步衍生出了一套雷達(dá)干涉技術(shù)- 地基合成孔徑 ( GBSAR) 。GBSAR
技術(shù)可以對所監(jiān)控的區(qū)域通過主動探測微波成像技術(shù),得到所監(jiān)控的區(qū)域的二維圖像,并且基于合成孔以及頻率步進(jìn)的理論來對所得圖像的方位以及距離進(jìn)行高空間分辨率的處理。目前,在各類( 大壩、建
筑物、滑坡、冰川以及橋梁) 變形以及位移監(jiān)測的研究中廣泛采用 GBSAR 技術(shù)。該技術(shù)的廣泛應(yīng)用有: 滑坡、冰川、大壩、建筑物和橋梁等變形和位移的監(jiān)測等 。并且與傳統(tǒng)的 GPS 測量方法所得的監(jiān)測結(jié)果進(jìn)行誤差對比分析,得到 GBSAR 技術(shù)在各類變性以及位移監(jiān)測中的監(jiān)控精度較高,滿足監(jiān)測要求 。為了充分發(fā)揮該技術(shù)的優(yōu)勢 拓寬其應(yīng)用范圍,需對各項誤差來源進(jìn)行分析并提出相應(yīng)的改正措施。
與此同時,GBSAR 技術(shù)在地形起伏波動較大的山區(qū),此類地形具有疊掩、陰影面積比例較大,采集到的干涉相位會出現(xiàn)不連續(xù)甚至全是噪的數(shù)字高程模型 DEM ( Digital Elevation Model ) 。同 時,
GBSAR 技術(shù)亦容易受大氣效應(yīng)、時間去相干等因素的影響,解決這一問題的一種思路是融合同一區(qū)域的多角度觀測數(shù)據(jù),從而在某一角度下的幾何畸變區(qū)域可以利用其他角度的數(shù)據(jù)補(bǔ)償 。
GBSAR 理論是由步進(jìn)頻率連續(xù)波 ( SFCW) 技術(shù)、合成孔徑雷達(dá)成像( SAR) 技術(shù)和差分干涉技術(shù) ( DIT)等 3 個關(guān)鍵技術(shù)構(gòu)成: 步進(jìn)頻率連續(xù)波 ( SFCW) 技術(shù)保證了 GBSAR 形變遠(yuǎn)距離測量、合成孔徑雷達(dá)成像( SAR) 技術(shù)保證了 GBSAR 形變大范圍測量、而差分干涉技術(shù)( DIT) 可以實現(xiàn)對形變的高精度測量。
1.雷達(dá)系統(tǒng)的誤差與改正
1.1 相位不穩(wěn)定性誤差與改正
相位穩(wěn)定性是雷達(dá)系統(tǒng)監(jiān)測精度的一個重要指標(biāo),并且雷達(dá)系統(tǒng)監(jiān)測精度受其參照物自身本振特性、信號發(fā)送以及數(shù)據(jù)接受設(shè)備與雷達(dá)系統(tǒng)間的傳輸效率以及傳輸質(zhì)量的影響。
相位的累計干涉表達(dá)式為:
N - 1 | |||||
φ = ∑ ∠{ e | j{ φ | i+1 | -φ(·) ] | } | ( 1) |
i |
i = 1
設(shè) R0 = R2 -R1 為無頻率偏移時的形變值,R'為中心頻率偏移 fD 時的形變值,則由雷達(dá)系統(tǒng)相位偏移造成的距離向偏移為:
R = | R'- | R0 | = - | λ | ( | fD | + | fD | ) | fD | = αf R2 | = | fD | R2 | ( 2) | |
2 | fR | f | ||||||||||||||
2 | R |
式中,λ 為雷達(dá)波長; αf = fD 為頻率偏移比。fR
為了減小相位偏移誤差對雷達(dá)系統(tǒng)監(jiān)測精度的影響,應(yīng)在監(jiān)測過程中,采用同一種設(shè)備,包括雷達(dá)系統(tǒng)的數(shù)據(jù)接收以及信號采集設(shè)備以及數(shù)據(jù)圖像經(jīng)過處理以后的高像素點數(shù)據(jù)圖像。與此同時,在長期、反復(fù)監(jiān)測過程中,發(fā)展具有多波段、多極化以及多個工作模式的地貌成像系統(tǒng),比如可以選擇穩(wěn)定性能優(yōu)異的頻率合成器,基于多梯次的校正方法對其采集到的數(shù)據(jù)圖像進(jìn)行精度校正分析,以此實現(xiàn)可提供質(zhì)量更高、數(shù)據(jù)更精準(zhǔn)的地基 SAR 差分干涉評價數(shù)據(jù)。
1.2 熱噪聲誤差與改正
雷達(dá)系統(tǒng)在發(fā)射信號、接收返回數(shù)據(jù)的電磁波以及數(shù)據(jù)存儲、反饋過程中由雷達(dá)系統(tǒng)特性而自發(fā)產(chǎn)生的信號熱噪聲。
雷達(dá)系統(tǒng)熱噪聲的高斯分布服從均值為 0,標(biāo)準(zhǔn)差為 σn ,并且與系統(tǒng)回波數(shù)據(jù)信號分別獨立統(tǒng)計處理。評估系統(tǒng)熱噪聲的指標(biāo)為統(tǒng)計數(shù)據(jù)的平均值來反映系統(tǒng)特性,因此,I/Q 兩路統(tǒng)一采用 n 表示噪聲。其概率密度函數(shù)為
M N M N
s = ∑∑ Aij exp{ j[θij +φij ]} = ∑∑ Aij exp{ jφ}
i = 1 j = 1 i = 1 j = 1
( 3)
式中,SNRi 為相位干涉通道的信噪比。
系統(tǒng)所產(chǎn)生的熱噪聲在統(tǒng)計學(xué)方面具有彼此獨立、圓高斯的特性,其導(dǎo)致的去相關(guān)效應(yīng)表達(dá)式為:
1
γn = ( 4) 1+SNR1-1 1+SNR2-1
信噪比( SNR) 是指雷達(dá)系統(tǒng)在所涉及的頻帶內(nèi)的輸入端所產(chǎn)生的信號與噪聲功率之間的比值,表達(dá)式為:
P | signal | V | , | ||||||||
SNR = 10lg( | ) dB = 20lg( | signal rms | ) dB | ( 5) | |||||||
Pnoise | Vnoise,rms | ||||||||||
首先回波的實、虛兩部的歸一化因子為: | |||||||||||
1 | MN | ( 6) | |||||||||
n = | ∑∑A2ij +σ2n | ||||||||||
2 | |||||||||||
i = 1 j = 1 |
而 SAR 疊加熱噪聲回波歸一化表達(dá)式為:
1 | M | N | n} | |||
NSr | = | [ ] + | + | |||
n | { ∑∑Aij cos φij | |||||
i = 1 | j = 1 |
N
j | [ ] + | n} | ( | 7 | ) |
n | { ∑∑Aij sin φij | ||||
i = 1 j = 1 |
數(shù)據(jù)獲取中的誤差與改正
地基 GBSAR 視向形變測量誤差主要是由干涉相位誤差所引起的,因此對干涉相位誤差要進(jìn)行長期以及重復(fù)觀測,以確保所監(jiān)測的數(shù)據(jù)結(jié)果能夠?qū)崟r且準(zhǔn)確的反映出監(jiān)測區(qū)域的干涉相位,其表達(dá)式為:
φ = φdisp +φgeom +φatmo +φnoise +δφ+εφ | ( 8) | |||
式中 | , | , | ||
φ 是主從影像計算所得干涉相位 φgeom 為系 | ||||
, | ||||
統(tǒng)設(shè)備安裝時發(fā)生的相位影響 φdisp 為主從影像體 | ||||
, | , | |||
現(xiàn)的變形相位 φatmo 是大氣擾動對相位的影響 φnoise | ||||
, , | 均為相位纏繞對監(jiān)測 | |||
為噪聲對相位的影響 δφ εφ | ||||
系統(tǒng)造成的影響。 | ||||
數(shù)據(jù)獲取中的誤差主要包括平臺偏移誤差和大 | ||||
氣擾動誤差。 | ||||
φ = φ-φdisp = φatmo +φnoise +δφ+εφ | ( 9) | |||
2.1 | 平臺偏移誤差與改正 |
地基 GBSAR 系統(tǒng)監(jiān)測平臺通常設(shè)置在地面、建筑物頂部或者以陸地各種交通工具上,此三類觀
測平臺容易發(fā)生一系列的微小且隨機(jī)性極強(qiáng)的微小的相位偏移,導(dǎo)致雷達(dá)系統(tǒng)的自調(diào)節(jié)判斷其改變監(jiān)測角度以及運(yùn)行軌道來與平臺偏移進(jìn)行匹配,這將嚴(yán)重影響相位數(shù)據(jù)圖片的相干性,使監(jiān)測精度大大的降低,而這種影響,將會給平臺的長期監(jiān)測帶來較為嚴(yán)重的問題。
假設(shè)平臺監(jiān)測到的兩幅復(fù)圖像分別為 I1 、I2 :
I1 = | I1 |·ejφ1 | |||||||||
I | = I ·ejφ2 | ( 10) | |||||||
2 | 2 | ||||||||
在 N×N 大小的匹配窗口區(qū)域內(nèi)形成干涉相位: | |||||||||
I = I ·I | = | I | | I | ej( φ1-φ2) | ( 11) | |||||
1 | 2 | 1 | 2 | ||||||
此種方法稱之為最大信噪比法,即最大分量與 | |||||||||
其他各個分量之和間的比值。 | |||||||||
SNR = | fmax | ( 12) | |||||||
f | , | - f | max | ||||||
∑ | m n |
2.2 大氣擾動誤差與改正
綜上,GBSAR 技術(shù)的監(jiān)測結(jié)果隨著監(jiān)測平臺以及監(jiān)測對象所處的環(huán)境實時動態(tài)的變化而發(fā)生較大范圍的變化。特別是大氣擾動因素對于 GBSAR 技術(shù)監(jiān)測結(jié)果的影響,在對 GBSAR 技術(shù)監(jiān)測雷達(dá)系統(tǒng)所采集到的實時數(shù)據(jù)處理過程中,大氣擾動是較為復(fù)雜的,也是亟待解決的一個關(guān)鍵性問題。本文選取某一個目標(biāo),假設(shè)其是穩(wěn)定,大氣擾動中較為重要的因素為大氣折射指數(shù),而大氣折射指數(shù)具有實時變化的差異性,設(shè)定 φ 為在 t1 和 t2 兩個不同時刻所存在的相位差:
φ( r,t) = φdisp( r,t) +φgeom( r,t) | ( 13) | ||||
式中, | ( r,t) 為監(jiān)測地形的形變相位的實際觀測 | ||||
φdisp | |||||
值; r = | r | | |||||
4πrn | |||||
φ = | [n( t2 ) -n( t1 ) ] | ( 14) | |||
λ | |||||
具有隨機(jī)性與多樣性的大氣擾動無處不在,即使在小尺度規(guī)模的空間上亦極大的影響著檢測精度,目前,對于雷達(dá)監(jiān)測系統(tǒng)實時采集數(shù)據(jù)的過程中大氣擾動級別為厘米級。目前,基于實測的大氣實時變化的氣象數(shù)據(jù)( 溫度、濕度、氣壓) 來補(bǔ)償校正GBSAR 大氣擾動誤差的理論建立大氣擾動模型,即可計算出較為精確的大氣折射率的實時變化情況,使之在補(bǔ)償值的校正下,實現(xiàn)對大氣擾動所產(chǎn)生的
。當(dāng)波長為 λ | 時,距離雷達(dá) rn | ||||
誤差的即使校正 | |||||
處的目標(biāo)點的回波相位表達(dá)式如下 | |||||
[珋,,( ) ] | |||||
K·h( t) = | φ0 r0 t h t | ( 15) | |||
r0 | |||||
式中,可實測溫度 T、相對濕度 H 和總氣壓 P,干氣
壓為 Pd 。
大氣擾動影響的目標(biāo)點相位一階差分校正模型為:
φcorr | ( r,t) = | φ | ( r,t) -[ | ( r | ,t) ]·r( 16) |
φ0 | 0 |
式中,φcorr( r,t) 表示校正后的差分相位 r0 表示視線方向距離。r 表示視線向距離。
大氣影響的目標(biāo)點的二階干涉相位校正模型可得的。因此,為實現(xiàn)高精度位移監(jiān)測,仍需優(yōu)化大氣擾動誤差的改正模型。
3.數(shù)據(jù)處理的誤差
3.1 殘余誤差與改正
雷達(dá)系統(tǒng)監(jiān)測平臺陣列的相位誤差估計均方差 ( ARMSE) 以及幅度誤差估計均方差 ( ARMSE) 分別為:
W | [( ) i | ( | ) i]2 | ||||||||||||
N | ∑ | - | φn | ||||||||||||
1 | φn | ( 17) | |||||||||||||
σp = | ∑ | i = 1 | |||||||||||||
N - 1 n - 2 | W | ||||||||||||||
W | |||||||||||||||
1 | N | ∑[( ρ珋n) i | - ( ρn ) i]2 | ||||||||||||
σa = | ∑ | i = 1 | ( 18) | ||||||||||||
N - 1 n - 2 | W |
式中,W 監(jiān)測采集數(shù)據(jù)的次數(shù),i 為第 i 次監(jiān)測的估計值。在每次實驗中,計算方程均會重新生成幅相誤差,并其在該次監(jiān)測過程中保持不變。
3.2 形變量解算誤差與改正
GBSAR 所采用的 2 維分辨成像技術(shù)是基于控制
信號發(fā)射設(shè)備的發(fā)射軌跡為直線,如圖 1 | 。 | |
所示 |
GBSAR y 方向為軌道方向,長為 L,與觀測目標(biāo)水平高度差為 H,發(fā)射信號的天線的照射俯角度為 θ,觀測范圍為 M。GBSAR 技術(shù)的實際監(jiān)測場景的大小( M) 遠(yuǎn)遠(yuǎn)大于方位向的軌道長度 L,如圖 1 所示,因此在對采集的圖像進(jìn)行處理時,在其方位向要對其進(jìn)行補(bǔ)
零。以保證補(bǔ)零后的方位向信號長度與實際觀測區(qū)域大小實現(xiàn)最優(yōu)匹配,避免了方位向采集的圖像變得過于模糊,從而影響了 GBSAR 圖像質(zhì)量。
在 GBSAR 技術(shù)監(jiān)測過程中的視線向形變真
值為: | ||||||||||||||
( , ) = | ( | , | ) × | 1 | - | 4πf | ( | ) ( | ) | |||||
xi | exp | j | 19 | |||||||||||
S f R | ∑ σ | yi | R( i) | [ | R i | ] | ||||||||
i | c |
式中,R( i) = ( xi ) 2 +( y-yi ) 2 +H2 為觀測目標(biāo)到雷達(dá)系統(tǒng)接受信號的天線的距離,σ( xi ,yi ) 為監(jiān)測區(qū)域目標(biāo)復(fù)散射系數(shù),f 為雷達(dá)監(jiān)測系統(tǒng)設(shè)備發(fā)射信號的頻率,c 為電磁波傳播速度。
4.GBSAR 誤差分析實驗
5.結(jié)論
GBSAR 技術(shù)由于可以獲得很高的監(jiān)測精度,是一種創(chuàng)新的并且得到廣泛應(yīng)用的形變監(jiān)測方法。其
監(jiān)測頻率在使用過程中可根據(jù)實際情況自由設(shè)定并且可以達(dá)到實時監(jiān)測,并且 GBSAR 技術(shù)完善了傳統(tǒng)技術(shù)的缺陷( 星載或機(jī)載 SAR 的失相干嚴(yán)重、時
空分辨率低) 。定性與定量地對 GBSAR 測量精度的誤差影響進(jìn)行分析。實驗表明了 IBIS-S 系統(tǒng)具有較為優(yōu)異的穩(wěn)定性,監(jiān)測系統(tǒng)誤差量級僅僅為亞毫米級; 并且得到大氣擾動因素為對相位形變誤差影響最為關(guān)鍵的因素,GBSAR 誤差隨著視線距離的增大而增加,精度隨著距離的增大而降低。在變形監(jiān)測中 GBSAR 技術(shù)具有實實在在的可實用性,該方法要求數(shù)據(jù)量不大,且適用的區(qū)域類型廣。亦可彌補(bǔ)傳統(tǒng)監(jiān)測技術(shù)的缺點。