毛曉康，汪長(zhǎng)城,周玉杉，王會(huì)強(qiáng)

(中南大學(xué) 地球科學(xué)與信息物理學(xué)院，湖南 長(zhǎng)沙 410083)

?

基于最小二乘平差的全極化SAR配準(zhǔn)偏移量估計(jì)方法

毛曉康，汪長(zhǎng)城,周玉杉，王會(huì)強(qiáng)

(中南大學(xué) 地球科學(xué)與信息物理學(xué)院，湖南 長(zhǎng)沙 410083)

全極化合成孔徑雷達(dá)影像(PolSAR)相對(duì)單極化SAR影像有更加豐富的地表信息。為了提高SAR影像偏移量跟蹤技術(shù)估算偏移量的精度，提出一種基于最小二乘平差法的全極化SAR數(shù)據(jù)偏移量估計(jì)方法。首先利用全極化SAR不同極化通道數(shù)據(jù)估算偏移量得到多個(gè)觀測(cè)值，然后通過最小二乘平差法對(duì)多余觀測(cè)值循環(huán)剔除粗差得到最優(yōu)偏移量。文中選取美國(guó)科羅拉多州湖城(Lake city)區(qū)域的兩景JPL無人機(jī)UAVSAR全極化SAR影像進(jìn)行實(shí)驗(yàn)，結(jié)果表明，文中新方法具有良好的粗差探測(cè)和去除功能，相對(duì)于已有研究結(jié)果，在方位向和距離向的偏移量估算精度都有明顯提高，達(dá)到15%～25%。新方法提高了偏移量跟蹤估算偏移量的精度，這對(duì)利用偏移量跟蹤技術(shù)監(jiān)測(cè)滑坡、地震及冰川等有重要的意義。

全極化SAR；最小二乘平差法；偏移量跟蹤；配準(zhǔn)；粗差

合成孔徑雷達(dá)(Synthetic Aperture Radar，SAR)圖像之間利用偏移量跟蹤技術(shù)來研究地物運(yùn)動(dòng)特性一直備受國(guó)內(nèi)外學(xué)者青睞[1]。尤其在研究冰川、地震及滑坡運(yùn)動(dòng)特性，取得了較好的成果[2-3]。目前大多數(shù)研究都是利用單極化SAR影像估算形變偏移量，其中一些偏移量很難區(qū)分是粗差還是形變信息[4-5]?？紤]到全極化SAR影像相對(duì)單極化SAR影像包含更豐富的地物散射信息[6]，在實(shí)際應(yīng)用表現(xiàn)出很多優(yōu)勢(shì)。例如極化干涉測(cè)量(PolInSAR)就能被用來反演森林高度及林下地形[7]。因而在提高偏移量跟蹤技術(shù)估算形變偏移量的精度有很大潛力需要去挖掘。綜合上述兩點(diǎn)，本文目的：①對(duì)比不同極化方式SAR偏移量結(jié)果；②利用全極化SAR影像來提高傳統(tǒng)單極化SAR影像偏移量跟蹤技術(shù)的精度，增強(qiáng)區(qū)分是否為形變信息的可靠性。

傳統(tǒng)的SAR影像偏移量形變監(jiān)測(cè)精度主要受影像空間分辨率與配準(zhǔn)誤差影響[8]。SAR影像空間分辨率是系統(tǒng)固有參數(shù)，難以通過數(shù)據(jù)處理手段提高。但配準(zhǔn)誤差可以在全極化SAR影像匹配中得到改善。不同極化電磁波與地物相互作用存在差異[9]，因此造成不同極化SAR影像配準(zhǔn)誤差不相同。本文首先利用全極化SAR數(shù)據(jù)分析不同極化方式SAR影像配準(zhǔn)偏移量估算的相關(guān)性。然后利用不同極化SAR組合匹配估算偏移量產(chǎn)生多余觀測(cè)值，提出基于最小二乘平差法的全極化SAR偏移量估算方法。最后選取覆蓋美國(guó)科羅拉多州湖城的JPL無人機(jī)UAVSAR全極化SAR影像進(jìn)行實(shí)驗(yàn)驗(yàn)證。

1 不同極化通道SAR偏移量相關(guān)性分析

PolSAR影像在地物的分類分割等相關(guān)領(lǐng)域已得到廣泛應(yīng)用[10]，主要是依據(jù)地表物體對(duì)不同極化方式電磁波的散射特性存在差異[6]。同一地物對(duì)不同極化電磁波反射能力有差異，在不同的地物類型、地表粗糙度及濕度等條件下也會(huì)有所不同[9]。為了分析不同極化SAR精配準(zhǔn)后估算偏移量的相關(guān)性，本文設(shè)計(jì)了下面的實(shí)驗(yàn)。

選取間隔7 d獲取的兩景美國(guó)JPL高分辨率無人機(jī)UAVSAR全極化SAR影像。影像覆蓋美國(guó)科羅拉多州湖城(Lake city)區(qū)域。截取影像中覆蓋Lake city的一塊區(qū)域進(jìn)行實(shí)驗(yàn)，其SAR強(qiáng)度影像如圖1所示。影像具體參數(shù)見表1。湖城位于美國(guó)克羅多河西南部圣胡安山脈下，北緯38°1′53″,西經(jīng)107°18′53″，四周高山環(huán)繞。湖城建筑相對(duì)密集，對(duì)雷達(dá)的信號(hào)反射復(fù)雜，同時(shí)產(chǎn)生陰影也比較多。由于路面對(duì)雷達(dá)信號(hào)發(fā)生鏡面反射，因此在SAR強(qiáng)度影像表現(xiàn)出較弱的信號(hào)。此外，SAR影像的噪聲也比較嚴(yán)重，主要是來自SAR系統(tǒng)成像固有的斑點(diǎn)噪聲[11]。SAR影像中的陰影和噪聲等嚴(yán)重影響配準(zhǔn)偏移量估算的精度[12-13]，甚至?xí)鸫植睢?/p>

表1 JPL無人機(jī)UAVSAR全極化SAR影像具體參數(shù)

圖1 2015-05-03獲取的HH極化方式SAR強(qiáng)度圖像及其地理位置圖

利用傳統(tǒng)的歸一化交叉相關(guān)(Normalized Cross Correlation, NCC)匹配算法，分別對(duì)HH、HV和VV 3種極化SAR影像做精配準(zhǔn)后估算偏移量，結(jié)果見圖2。其中(a)、(b)和(c)分別代表SHH、SHV和SVV在方位向的偏移量，(d)、(e)和(f)分別表示SHH、SHV和SVV距離向的偏移量。由于影像覆蓋區(qū)域?yàn)槌菂^(qū)，且獲取時(shí)間間隔短，因此在不考慮軌道誤差的影響情況下可假設(shè)偏移量為零。本文選取偏移量絕對(duì)值大于3倍中誤差視為粗差，可知圖中不同極化得到的粗差區(qū)域有明顯的差異。用皮爾遜相關(guān)系數(shù)分析[14]統(tǒng)計(jì)3種極化估算的偏移量之間的相關(guān)性。SHH與SHV在方位向偏移量相關(guān)值為0.24，距離向?yàn)?.19；SHH與SVV在方位向?yàn)?.42，距離向?yàn)?.46；SHV與SVV在方位向?yàn)?.25，距離向?yàn)?.28。它們之間的相關(guān)都未超過0.5，因此不同極化SAR影像估算的偏移量之間的相關(guān)性較小。

圖2 不同極化方式SAR影像方位向和距離向偏移量估算結(jié)果比較

2 全極化SAR數(shù)據(jù)偏移量估算方法

2.1 單極化SAR數(shù)據(jù)偏移量估算算法

偏移量跟蹤技術(shù)首先利用兩景相關(guān)的影像中的像元進(jìn)行精配準(zhǔn)，得到兩景影像配準(zhǔn)同名像元的坐標(biāo)偏移量。坐標(biāo)偏移量中包含：①地形起伏，電離層等因素引起的偏移量；②獲取影像的軌道差異引起的偏移量；③地表形變引起的偏移量。由地形起伏引起的偏移量在基線較短和地形相對(duì)平坦的時(shí)候很小[15]，一般忽略不計(jì)?？紤]機(jī)載無人機(jī)平臺(tái)低空飛行，因而不考慮電離層、大氣等因素引起的偏移量。因此坐標(biāo)偏移量Roffset就等于軌道偏移量Rorbit和形變偏移量Rchange之和。

Roffset=Rorbit+Rchange.

(1)

由于影像中絕大部分地表處于穩(wěn)定狀態(tài)，其理論偏移量為零，但是在實(shí)際應(yīng)用中受到軌道偏移量的影響而與理論值產(chǎn)生偏離。因此可以利用穩(wěn)定區(qū)域的偏移量來近似估計(jì)軌道引起的偏移量，即假設(shè)Roffset≈Rorbit。一般情況下，采用二次方程來建立軌道偏移量模型，并利用最小二乘方法擬合模型參數(shù)。

Rorbit=a0+a1·l+a2·p+a3·l·p.

(2)

其中:a0,a1,a2,a3為待定系數(shù)；l,p分別代表像元在SAR影像中的行列數(shù),行代表方位向，列代表距離向。在實(shí)驗(yàn)過程中，利用最小二乘法對(duì)模型參數(shù)進(jìn)行擬合(采樣點(diǎn)為2 000)，模型的內(nèi)符合精度一般作為評(píng)價(jià)偏移量精度的一個(gè)指標(biāo)。

對(duì)于坐標(biāo)偏移量的計(jì)算，核心是尋找主從影像同名像點(diǎn)。尋找同名像點(diǎn)的算法有很多，一般利用基于強(qiáng)度信息的相關(guān)系數(shù)法。相關(guān)系數(shù)計(jì)算公式如下[16]：

(3)

當(dāng)相關(guān)系數(shù)取得最大值時(shí)的位置為最佳匹配點(diǎn)。主從影像的同名像點(diǎn)的位置用slave(l′,p′),master(l,p)表示，按式(4)計(jì)算偏移量。

Rorbit=slave(l′,p′)-master(l,p).

(4)

2.2 全極化SAR數(shù)據(jù)偏移量估算方法

(5)

(BTPB)X=BTPL.

(6)

解得

X=(BTPB)-1BTPL.

(7)

X的中誤差為

(8)

3 實(shí)驗(yàn)驗(yàn)證

為了驗(yàn)證本文算法，利用圖1對(duì)應(yīng)的兩景JPL無人機(jī)UAVSAR全極化SAR影像進(jìn)行實(shí)驗(yàn)。方位向和距離向的偏移量估算結(jié)果如圖3、圖4所示。圖3(a)～(f)分別為SHH，SHV，SVV，SHH+SVV，SHH+SHV+SVH+SVV及平差方法得到的方位向偏移量估算結(jié)果。圖4(a)～(f)分別為上述方法在距離向偏移量估算結(jié)果。從圖3和圖4可以看出，本文提出的方法的結(jié)果(圖3(f)和圖4(f))中偏移量粗差明顯少于單極化SAR影像估算結(jié)果，尤其是在距離向比較明顯(見圖4中的黑色圓框所示區(qū)域)。

圖3 不同組合方位向形變偏移量

圖4 不同組合的距離向偏移量

另外，為了分析不同的匹配窗口大小對(duì)結(jié)果的影響，分別利用4種不同大小的匹配窗口64×32(方位向×距離向)、96×48、128×64及160×80進(jìn)行實(shí)驗(yàn)驗(yàn)證。統(tǒng)計(jì)不同方法在不同大小匹配窗口條件下的偏移量估算結(jié)果的中誤差，結(jié)果如表2所示。從表2可知，隨著窗口增大，不同方法在方位向及距離向偏移量估算精度均隨之提高，然而結(jié)果的細(xì)節(jié)信息會(huì)隨著匹配窗口的增大而損失；本文方法相對(duì)單極化SAR結(jié)果在不同的匹配窗口條件下精度均提高15%～25%，且在匹配窗口較小時(shí)精度改善更為明顯。

表2 不同極化組合偏移量中誤差

本文方法利用小的匹配窗口估算偏移量精度與單極化SAR利用大的匹配窗口估算偏移量的精度相當(dāng)。例如本文方法選擇96×48的匹配窗口偏移量精度：方位向5.654 0 cm；距離向1.990 4 cm。傳統(tǒng)單極化SAR只能選擇128×64的匹配窗口才能達(dá)到與本文方法相近的精度，方位向5.615 4 cm；距離向1.959 5 cm。同理其他窗口也有類似規(guī)律。詳見表2中虛線框和實(shí)線框分別對(duì)應(yīng)的結(jié)果。

4 結(jié) 論

1)利用偏移量跟蹤算法處理不同極化方式SAR影像，估算的偏移量存在差異，尤其是粗差區(qū)域更明顯。用皮爾遜相關(guān)系數(shù)分析得它們之間偏移量的相關(guān)性都在0.5以下。

2)利用本文提出的基于最小二乘平差的全極化SAR配準(zhǔn)偏移量估計(jì)方法，處理高分辨率全極化機(jī)載SAR數(shù)據(jù)。結(jié)果表明，本文方法能夠較好地剔除單極化SAR估算偏移量中的粗差。對(duì)比傳統(tǒng)單極化SAR處理的結(jié)果，本文方法偏移量估計(jì)精度提高了15%～25%。

3)比較不同大小的匹配窗口估算偏移量，本文方法對(duì)比傳統(tǒng)方法，精度均得到提高，且在匹配窗口較小時(shí)精度改善更為明顯。傳統(tǒng)方法一般通過增大匹配窗口大小，提高偏移量估算精度，而本文方法利用小的匹配窗口(96×48)估算的精度與傳統(tǒng)方法選擇大的匹配窗口(128×64)的精度相當(dāng)。因此，本文方法在保證精度不變的情況下，能降低匹配窗口的大小，不僅減少大量的計(jì)算，而且還能保留更多細(xì)節(jié)信息。

[1] STROZZI T, LUCKMAN A, MURRAY T, et al. Glacier motion estimation using SAR offset-tracking procedures[J]. Geoscience and Remote Sensing, IEEE Transactions on. 2002,40(11): 2384-91.

[2] 李佳，李志偉，汪長(zhǎng)城，等．SAR偏移量跟蹤技術(shù)估計(jì)天山南依內(nèi)里切克冰川運(yùn)動(dòng)[J]．地球物理學(xué)報(bào)，2013，56(4)：1226-1236.

[3] HU J, LI Z W, DING X L, et al. Two-dimensional co-seismic surface displacements field of the Chi-Chi earthquake inferred from SAR image matching[J]. Sensors. 2008, 21;8(10):6484-95.

[5] SINGLETON A, LI Z, HOEY T, et al. Evaluating sub-pixel offset techniques as an alternative to D-InSAR for monitoring episodic landslide movements in vegetated terrain[J]. Remote Sensing of Environment. 2014,5;147:133-44.

[6] ZEBKER H. Polarisation: Applications in Remote Sensing[J]. Physics Today. 2010,1;63(10):53-4.

[7] 付海強(qiáng), 汪長(zhǎng)城, 朱建軍，等. 一種改進(jìn)的PolInSAR PCT方法反演植被垂直結(jié)構(gòu)[J]. 測(cè)繪工程, 2014,23(11):56-61.

[8] 陳強(qiáng),羅容,楊瑩輝,等.利用SAR影像配準(zhǔn)偏移量提取地表形變的方法與誤差分析[J].測(cè)繪學(xué)報(bào), 2015,44(3):301-308.

[9] 王超.全極化合成孔徑雷達(dá)圖像處理[M].北京：科學(xué)出版社, 2008.

[10] 王興玲,胡德勇,唐宏, 等. 基于Bayes決策的機(jī)載全極化SAR圖像滑坡信息提取[J].國(guó)土資源遙感, 2014,26(2):121-127.

[11] LEE J S, POTTIER E. Polarimetric radar imaging: from basics to applications[M]. CRC press, 2009.

[12] 劉國(guó)祥,丁曉利,李志林, 等.星載SAR復(fù)數(shù)圖像的配準(zhǔn)[J].測(cè)繪學(xué)報(bào), 2001,30(1):60-66.

[13] ZITOVA B, FLUSSER J. Image registration methods: a survey [J]. Image and vision computing, 2003,21(11):977-1000.

[14] DUTILLEUL P, STOCKWELL J D, FRIGON D. The Mantel Test versus Pearson's Correlation Analysis: Assessment of the Differences for Biological and Environmental Studies[J]. Journal of Agricultural Biological & Environmental Statistics, 2000,5(2):131-150.

[15] STROZZI T, KOURAEV A, WIESMANN A, et al. Estimation of Arctic glacier motion with satellite L-band SAR data [J]. Remote Sensing of Environment. 2008,112(3):636-45.

[16] YOO J C, HAN T H. Fast Normalized Cross-Correlation[J]. Circuits Systems & Signal Processing, 2009,82(2):144-156.

[責(zé)任編輯：劉文霞]

An offset-tracking method based on least square adjustment of full polarimetric

SARMAO Xiaokang，WANG Changcheng，ZHOU Yushan，WANG Huiqiang

(School of Geosciences and Info-Physics Central South University,Changsha 410083,China)

Comparing with single-channel SAR image, the full polarimetric synthetic aperture radar image (PolSAR) has more abundant terrain information. To improve the offset estimation accuracy of SAR images, the paper proposes an offset-tracking method based on least square adjustment of full polarimetric SAR data. Firstly, the multiple offset observes are obtained of different polarimetric channel SAR images. Then the least square adjustment is used to eliminate the errors and get the optimal results. Finally, two JPL UAVSAR full polarimetric SAR images are selected to cover Lake City, Colorado, USA for validating the proposed method. The experimental results show that, compared with the single-channel SAR image, the offset errors of the proposed method decrease significantly. The precision of offsets in range and azimuth improves by 15%～25%. The proposed method improves the accuracy of offset tracking. It is significant to use the offset-tracking method for monitoring deformation of landslide, earthquake and glacier.

full polarimetric SAR; least squares adjustment; offset tracking; registration；error

10.19349/j.cnki.issn1006-7949.2017.03.005

2016-03-04

國(guó)家自然科學(xué)基金資助項(xiàng)目(41371335)

毛曉康(1990-)，男，碩士.

P237

A

1006-7949(2017)03-0023-05

引用著錄：毛曉康，汪長(zhǎng)城，周玉杉，等.基于最小二乘平差的全極化SAR配準(zhǔn)偏移量估計(jì)方法[J].測(cè)繪工程，2017，26(3)：23-27.