多基線極化干涉SAR植被高度反演方法

張 兵,朱建軍,付海強,何永紅,2,賀文杰,孫明星

(1.中南大學 地球科學與信息物理學院，湖南 長沙 410083；2.湖南科技學院 土木與環境工程學院，湖南 永州 425199)

多基線極化干涉SAR植被高度反演方法

張 兵1,朱建軍1,付海強1,何永紅1,2,賀文杰1,孫明星1

(1.中南大學 地球科學與信息物理學院，湖南 長沙 410083；2.湖南科技學院 土木與環境工程學院，湖南 永州 425199)

針對傳統單基線PolInSAR三階段植被高度算法，反演易受到噪聲因素影響導致反演結果精度較低等問題。文中采用幾何結構最優法和高程精度最優法兩種多基線PolInSAR植被高度反演方法，結合瑞典北部E-SAR P波段數據，驗證得出兩種多基線算法植被高度反演的精度以及其在不同植被高度覆蓋區的適用條件。實驗結果表明兩種多基線方法植被高度反演精度相對傳統單基線算法分別提高44.05%和22.11%，且當兩種多基線方法同時反演植被高度時，植被高度大于10 m時采用幾何結構最優法較好；植被高度小于10 m時采用高程精度最優法。

極化干涉SAR;三階段算法;多基線;植被高反演

植被高度是林業資源調查、生態環境管理以及碳匯量評估重要的基礎數據。傳統植被高監測手段因其費時、費力，難以實現大范圍快速監測，無法滿足當前對植被高數據的迫切需求。極化合成孔徑雷達干涉測量(Polarimetric SAR Interferometry, PolInSAR)既有InSAR大范圍、高空間分辨率、全天時、全天候測高能力，又有PolSAR對植被散射空間分布、形狀、方向敏感的特性，為大范圍、快速獲取植被高度提供新的契機。鑒于此，積極開展PolInSAR植被高度反演具有重要的現實意義。

目前PolInSAR植被高度提取應用最廣泛的模型為隨機地體二層散射模型(Random Volume over Ground, RVoG)，該模型建立了InSAR相干性與植被生物物理參數間的關聯[1-3]。2001年Papathanassiou首次建立PolInSAR觀測數據與RVoG模型之間的關聯，將單基線PolInSAR植被高反演問題概括為六維非線性求解問題[3]。但是該方法初值難以確定，容易陷入局部最優問題。針對這一問題，2003年Cloude提出單基線PolInSAR植被高度反演的三階段算法[4]，該方法建立了RVoG假設條件下復相干系數的幾何表達，極大程度簡化了植被高度反演的復雜程度，是目前最為通用的植被高度提取算法。Cloude等[4]認為在復平面單位圓內極化復相干系數呈直線分布，直線擬合的精度直接影響著最終植被高度的估計精度，然而實際情況下存在很多噪聲因素影響復相干系數，其位置也會偏離理想直線呈現區域分布，且此區域隨著誤差增大而增大。當噪聲過大時，會造成擬合直線精度較低，導致地表相位估計錯誤，植被高度反演失敗；此外，垂直向有效波數kz影響著體去相干系數的確定。kz的大小影響著植被高度反演的敏感程度以及所能反演的植被高度范圍，過大或過小都會造成植被高度反演精度的降低。

針對傳統單基線三階段算法直線擬合精度較低及體去相干系數估計精度較差等問題，Lee提出兩種多基線PolInSAR植被參數反演方法[7]，幾何結構最優法及高程精度最優法，在具有多基線數據的情況下實現最優基線的自主選擇，在一定程度上提高了植被高度反演的精度；然而目前已有文獻中并沒有考慮兩種多基線方法在不同植被高度覆蓋區的適用條件，本文以此為基礎，首先通過兩種多基線方法植被高度反演結果與LiDAR測量結果對比，之后根據樣地區域植被高度反演結果對比，分析兩種多基線PolInSAR植被高度反演方法的適用條件，以提高兩種多基線方法的適用性，為后續研究提供基礎。

1 PolInSAR植被高度反演原理

經典三階段算法基于隨機地體二層散射模型(RVoG)，RVoG模型將PolInSAR觀測量與植被參數間的關系概括為[3]

(1)

其中,ω為各極化狀態，表達目標的散射特性；γ(ω)為不同極化狀態ω對應的極化復相干系數，已知復數；i為復數標識；φ0為地表相位，未知實數；μ(ω)為各極化狀態ω對應的地體散射幅度比，未知實數；γν為完全由植被冠層貢獻的“純”體去相干系數：

(2)

如圖1所示，從幾何角度看，式(1)表示復平面內通過點P和Q的直線，RVoG模型具有線性特征，其相干跡為一條直線，即為相干直線。這條直線與單位圓的交點為地表相位點，可視線段中最遠離地表相位的復相干系數點為體散射復相干系數點。傳統單基線三階段算法分為三個步驟：整體最小二乘直線擬合、地表相位和體去相干系數確定、建立二維查找表反演植被參數[4，9-11]。直線擬合的精度直接影響著地表相位估計的精度。計算中，一般假設體散射通道的地體幅度比為0進行參數反演，當體散射去相干存在誤差時，也會影響著三階段算法的反演精度。

圖1 RVoG模型相干跡

2 多基線PolInSAR植被高度反演

針對傳統單基線三階段算法的不足，本文對比兩種多基線算法原理，可以在具有多基線數據的情況下自主選擇最優基線進行參數解算，得到的參數解算結果最為可靠。

2.1 幾何結構最優法

對任意給定的極化干涉矩陣，在復平面單位圓內均可以形成子區域，可以覆蓋所有復相干系數在單位圓內的位置，這個子區域即為極化干涉矩陣的相干區域。不失一般性，本文定義一個橢圓相干區域[1，5-8]。

SAR傳感器在重軌飛行模式時，獲得基線兩端的散射矢量k1和k2，得到極化干涉相干矩陣T6(6×6)：

(3)

其中，*T表示共軛轉置，T11,T22,Ω12均為(3×3)的矩陣。

重新定義復數矩陣Π：

(4)

對矩陣Π進行舒爾分解，可得

Π=[U3]*T[R][U3]=

(5)

(6)

其中，a代表橢圓長軸長；b代表橢圓短軸長。

如圖2所示，相干邊界的幾何形狀可以作為評判各基線估計精度的標準：橢圓離心率越大，橢圓線性程度越高，參數估計精度越高，反演結果越準確；理想狀態下，橢圓離心率為1，此時相干區域退化為一條直線，即符合RVoG模型假設，顯然這種情況下參數反演結果最為準確。因此，當PolInSAR觀測值存在多條基線時，根據式(6)以及植被高度三階段反演算法，最優植被高度反演結果可得

(7)

其中，N代表觀測值基線條數。

圖2 復相干系數在復平面單位圓內分布

2.2 高程精度最優法

傳統InSAR高度反演結果由干涉相位φInt標準差與垂直向有效波數kz共同確定，因此可以應用高程精度標準差最小的方法確定最優基線：

(8)

(9)

式中：γ為復相干系數，LInt為進行相干性估計時的所用窗口數，δInt為干涉相位標準差。

當雷達回波信號的強度受到非體去相干因素(如時間去相干等)影響而減弱時，干涉相位標準差δInt因此增大，高程反演精度降低。也就是說，高程精度標準差Haccuracy越小，植被高度反演精度越高。

當PolInSAR觀測值存在多條基線時，根據式(8)，式(9)以及植被高度三階段反演算法，最優植被高度反演結果可得如下約束：

(10)

其中，N代表觀測值基線條數。

3 實驗與分析

3.1 實驗數據

為驗證本文采用的兩種多基線PolInSAR植被高度反演方法的可靠性，采用BioSAR 2008項目獲取的瑞典北部全極化數據對本文算法進行實驗，選取的實驗區域內植被以針葉林為主，采用E-SAR P波段數據進行實驗，P波段雷達波長較長，對于植被覆蓋層有較好的穿透能力，更好地獲取植被層垂直結構信息。圖3為實驗區域的影像圖以及Pauli基彩色合成圖。表1、表2分別為基線參數以及SAR影像參數信息表。本文中分別應用兩種最優基線選擇方法選擇基線，同時采用經典三階段算法反演植被高度。采用PD(相位最大分離算法：PDHigh、PDLow)和MCD(相干最大分離算法：opt1、opt2、opt3)等5種極化方式[1，9，16]所對應的復相干系數進行參數解算。

圖3 實驗區域影像圖以及Pauli彩色合成圖

表1 基線參數

表2 SAR影像參數

3.2 實驗結果

如圖4所示，為傳統單基線三階段算法、高程精度最優法以及幾何結構最優法植被高度反演結果，整體看，傳統單基線三階段算法相較其他兩種方法反演結果呈現低估。為了對植被高度反演結果進行定量分析，以實驗區域LiDAR測量結果作為真值進行精度評定。在實驗區域內選取22塊樣地對反演的植被高度值進行數據分析，分別計算三種方法植被高度反演方法的均方根誤差，傳統單基線三階段算法為5.29 m，高程精度最優法為4.12 m，幾何結構最優法為2.96 m。高程精度最優法相對傳統單基線三階段算法精度提高了22.11%，幾何結構最優法提高了44.05%。分析原因：①幾何結構最優法將相干區域的線性程度作為判斷標準，自主選擇最優基線進行植被高度反演。相干邊界橢圓離心率越大，相干區域越趨于直線，說明復相干系數噪聲水平越低，噪聲對植被高反演貢獻越小。②高程精度最優法，基于相干性的先驗統計特性，考慮了噪聲對干涉SAR測高精度的影響，可以在多基線數據中選擇測高精度最高的基線進行植被高度反演。

圖4 植被高度反演結果對比

如圖5所示，植被高度小于10 m，高程精度最優法相較幾何結構最優法植被高度反演結果更為準確；植被高度大于10 m時，幾何結構最優法相較高程精度最優法植被高度反演結果更為準確。分析原因在于：在低植被區，由于P波段雷達穿透能力較強，各極化通道相位中心相對集中且貼近于地面，地表相位估計精度差別不大；高程精度最優法通過衡量干涉測量精度，選擇體散射通道高程精度最優的基線進行植被參數反演，在一定程度上提高了植被高度反演的精度；高植被區不同極化通道相位中心差異較大，復相干系數在復平面單位圓內分布較為離散，體散射通道高程精度不再對植被高度反演起決定性因素，地表相位估計精度對植被高度反演精度影響較大。幾何結構最優法可以選取直線擬合最為準確的基線，其地表相位估計最為準確，地表相位偏差對于植被高度反演的貢獻最小，因而反演的植被高度最為準確。

圖5 樣地區域植被高度反演結果對比

4 結 論

1)傳統單基線三階段植被高度反演方法在植被高度解算過程中可能因為去相干因素以及非去相干因素影響而導致反演結果失敗。在具有多基線數據的情況下，本文對比了兩種多基線PolInSAR植被高度反演方法：①幾何結構最優法；②高程精度最優法，并應用BioSAR2008項目獲得的E-SAR P波段數據對兩種多基線方法進行實驗驗證，理論和實驗結果證明多基線PolInSAR植被高度反演方法行之有效且結果較傳統單基線三階段算法分別提高了44.05%和22.11%。

2)針對兩種多基線算法反演結果，得出如下結論：兩種多基線方法聯合反演植被高度時，植被高度小于10 m時采用高程精度最優法；植被高度大于10 m采用幾何結構最優法。

[1] TREUHAFT R N, MADSEN S N, MOGHADDAM M, et al. Vegetation characteristics and underlying topography from interferometric radar[J]. Radio Science, 1996, 31(31):1449-1486.

[2] CLOUDE S R, PAPATHANASSIOU K P. Polarimetric SAR Interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 1998,36(5):1551-1565.

[3] PAPATHANASSIOU K P, CLOUDE S R. Single-baseline Polarimetric SAR Interferometry[J]. Radio Science, 2000, 35(1):141-177.

[4] CLOUDE S R, PAPATHANASSIOU K P. Three-stage Inversion Process for Polarimetric SAR Interferometry[J]. Journal of Remote Sensing, 2002,6(6):424-429.

[5] SEYMOUR S, CUMMING I G. Maximum likelihood estimation for SAR interferometry. Proceeding of IEEE Geoscience and Remote Sensing Symposium, IGARSS’94, Pasadena, USA, 1994.

[6] CLOUDE S R. Polarisation: Applications in Remote Sensing[M].London: Oxford University Press,2009.

[7] Lee S K, Kugler F, Papathanassiou K, et al. MULTIBASELINE POLARIMETRIC SAR INTERFEROMETRY FOREST HEIGHT INVERSION APPROACHES[C]// ESA POLinSAR Workshop. 2011.

[8] FLYNN T, TABB M, CARANDE R. Coherence region shape extraction for vegetation parameter estimation in polarimetric SAR interferometry[C]// Geoscience and Remote Sensing Symposium, 2002. IGARSS '02. 2002 IEEE International. 2002:2596-2598 vol.5.

[9] 伍雅晴,朱建軍,付海強,等. 引入PD極化相干最優的三階段植被高度反演算法[J]. 測繪通報,2016(5):32-35，76.

[10] 許麗穎,李世強,鄧云凱,等. 基于極化干涉SAR反演植被高度的改進三階段算法[J]. 雷達學報,2014(1):28-34.

[11] 李廷偉,梁甸農,朱炬波. 極化干涉SAR森林高度反演綜述[J]. 遙感信息,2009(3):85-91.

[12] 李哲,陳爾學,王建. 幾種極化干涉SAR森林平均高反演算法的比較評價[J]. 遙感技術與應用,2009(5):611-616.

[13] 楊震,楊汝良. 極化合成孔徑雷達干涉技術[J]. 遙感技術與應用,2001(3):139-143.

[14] 陳爾學,李增元,龐勇,等. 基于極化合成孔徑雷達干涉測量的平均樹高提取技術[J]. 林業科學,2007(4):66-70，145.

[15] CLOUDE S R, CHEN H,GOODENOUGH D G. Forest height estimation and validation using Tandem-X polinsar. 2013 IEEE International Geoscience and Remote Sensing Symposium-IGARSS, 2013,07:1889-1892.

[16] 付海強. 基于測量平差理論的PolInSAR植被垂直結構提取模型與方法[D].長沙：中南大學,2014.

[17] 史磊. 多基線POLINSAR植被高度與植被下地形估計[D].武漢：武漢大學,2013.

[18] 宋桂萍. 極化干涉SAR植被高度反演算法研究[D].長沙：中南大學,2013.

[19] 羅環敏. 基于極化干涉SAR的森林結構信息提取模型與方法[D].成都：電子科技大學,2011.

[20] TREUHAFT R N, SIQUEIRA P R. Vertical Structure of Vegetated Land Surfaces from Interferometric and Polarmetric Radar[J]. Radio Science, 2000, 35(1):141-177.

[21] 郭華東. 雷達對地觀測理論與應用[M]. 北京:科學出版社,2000.

[責任編輯：李銘娜]

Multi-baseline PolInSAR vegetation height inversion method

ZHANG Bing1,ZHU Jianjun1,FU Haiqiang1,HE Yonghong1,2,HE Wenjie1,SUN Mingxing1

(1. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China; 2. School of Civil and Environmental Engineering, Hunan University of Science and Engineering, Yongzhou 425199, China)

This paper aims at the problems that the traditional single-baseline three stage inversion process is easily affected by noise factors, which will lead to a less precise inversion result. Two kinds of multi-baseline PolInSAR vegetation height inversion algorithms are adopted, which are the optimal method of geometric structure and the optimal method of vertical accuracy. Combined with the northern Sweden E-SAR P band data, the precision of the two algorithms for multi-baseline vegetation height inversion and their applied conditions in different vegetation areas is verified and obtained. The experimental result shows that the vegetation height inversion precision of the two methods has each increased by 44.05% and 22.11%, comparing with the traditional single-baseline algorithm. Besides, when the two methods works at the same time, it is better to adopt the optimal method of geometric structure when the vegetation height is more than 10m, while the optimal method of vertical accuracy is better when it is less than 10m.

polarimetric SAR interferometry；three-stage algorithm；multi-baseline；vegetation height inversion

2016-08-08

國家自然科學基金資助項目(41531068；41274010；41371335)；國家973計劃項目(2013CB733303)；國家863計劃項目(2012AA121301)；湖南省研究生科研創新項目(150140004)；湖南省自然科學基金資助項目(14JJ2131)；中南大學中央高校基本科研業務費專項資金資助項目(2016zzts083)

張 兵(1992-)，男，碩士研究生.

朱建軍(1962-)，男，教授，博士生導師.

著錄：張兵,朱建軍,付海強,等.多基線極化干涉SAR植被高度反演方法[J].測繪工程，2017，26(9)：23-27，31.

10.19349/j.cnki.issn1006-7949.2017.09.005

P228

A

1006-7949(2017)09-0023-05