肖虹雁，田傳召，岳彩榮

(西南林業(yè)大學(xué)林學(xué)院，云南 昆明 650224)

基于全極化SAR數(shù)據(jù)的云南松后向散射特征分析

肖虹雁，田傳召，岳彩榮

(西南林業(yè)大學(xué)林學(xué)院，云南 昆明 650224)

以昆明市宜良縣花園林場(chǎng)為研究區(qū)，選擇2009年4月23日ALOS PALSAR 1.1等級(jí)的數(shù)據(jù)產(chǎn)品L波段全極化單視復(fù)型數(shù)據(jù)，結(jié)合2007年完成的宜良縣森林資源二類調(diào)查數(shù)據(jù)，分析HH、HV、VH和VV 4種極化方式下云南松后向散射特征，采用SPSS進(jìn)行單因素方差分析法統(tǒng)計(jì)試驗(yàn)結(jié)果。結(jié)果表明，在云南松幼齡林、中齡林和近熟林中，HH極化方式下反射強(qiáng)度值都是最高的，其次是VV極化方式，而HV和VH極化方式下反射值都比較低；云南松幼齡林、中齡林和近熟林與后向散射系數(shù)之間具有極顯著影響，根據(jù)此特性，區(qū)分不同齡組的云南松。

全極化SAR；云南松；后向散射系數(shù)；單因素方差分析

云南松(Pinusyunnanensis)多分布于海拔1000～3000 m的廣大地區(qū)，常形成大面積的純林，也是云南省最重要和面積最大的森林植被類型[1]。大面積分布于滇中高原、金沙江流域及南盤江流域。云南松樹(shù)干通直，木質(zhì)輕軟細(xì)密，是良好的建筑用材；而云南松樹(shù)脂、松花粉、松尖、球果等體現(xiàn)了云南松的非木材資源價(jià)值。

合成孔徑雷達(dá)(SAR)以其自身優(yōu)勢(shì)以及極化組合多樣的特點(diǎn)，成為森林研究的一種有效手段。全極化SAR的極化方式是指發(fā)射和接收電磁波所采用的極化方式，可以有4種組合HH、HV、VH和VV，對(duì)于互易介質(zhì)，HV=VH，通過(guò)它們可以保持圖像空間分辨率同時(shí)可以減小信號(hào)干擾[2]。極化不同，目標(biāo)對(duì)電磁波的響應(yīng)就不同。全極化SAR具有對(duì)地表植被散射體的空間分布和高度很敏感的特性，又具有對(duì)植被散射體的形狀和方向很敏感的特性[3]，在對(duì)植被探測(cè)上具有獨(dú)特的優(yōu)勢(shì)。因此，本文利用ALOS PALSAR全極化SAR數(shù)據(jù)，結(jié)合2007年完成的宜良縣森林資源二類調(diào)查數(shù)據(jù)，分析HH、HV、VH和VV等4種極化方式下云南松后向散射特征。

1 研究區(qū)與數(shù)據(jù)

1.1 研究區(qū)概況

研究區(qū)位于云南省昆明市宜良縣花園林場(chǎng)，地勢(shì)北高南低。海拔1300～2800 m，地貌特點(diǎn)是山地、谷地與壩子相間分布。該地氣候溫和，年平均氣溫16.3 ℃，屬于北亞熱帶季風(fēng)氣候。森林起源以天然林為主，主要優(yōu)勢(shì)樹(shù)種為云南松。調(diào)查林分為云南松純林，群落結(jié)構(gòu)簡(jiǎn)單，各林班平均樹(shù)高為1.2～16 m，平均胸徑為2～31.2 cm，郁閉度為0.2～0.85，其中云南松幼齡林占22%，云南松中齡林占30.3%，云南松近熟林占47.7%。

1.2 數(shù)據(jù)的獲取

1.2.1 SAR數(shù)據(jù) 本試驗(yàn)采用2009年4月23日獲取的一景覆蓋試驗(yàn)區(qū)的ALOS PALSAR 1.1等級(jí)的數(shù)據(jù)產(chǎn)品L波段全極化(HH、HV、VH和VV)單視復(fù)型(SLC)數(shù)據(jù)。圖像中心點(diǎn)坐標(biāo)：北緯24°53′31″、東經(jīng)103°12′43″，其中方位向分辨率為9.37 m，距離向分辨率為3.56 m，像幅為18432×1248。

1.2.2 2007年宜良縣森林資源二類調(diào)查數(shù)據(jù) 2007年宜良縣森林資源二類調(diào)查數(shù)據(jù)[4]中小班因子調(diào)查數(shù)據(jù)，內(nèi)容包括編號(hào)、圖色名稱、鄉(xiāng)鎮(zhèn)、林班號(hào)、小班號(hào)、小班面積、地類、立地類型、優(yōu)勢(shì)樹(shù)種、起源、樹(shù)種組成、齡組、林木蓄積、經(jīng)營(yíng)措施、土壤類型、土壤質(zhì)地、土壤厚度、坡度、坡向、權(quán)屬、海拔等屬性數(shù)據(jù)項(xiàng)。根據(jù)2007年宜良縣森林資源二類調(diào)查數(shù)據(jù)中的齡組可知云南松生長(zhǎng)時(shí)期，而立地類型可以看出云南松生長(zhǎng)環(huán)境。

1.3 數(shù)據(jù)處理

對(duì)ALOS PALSAR數(shù)據(jù)預(yù)處理過(guò)程包括輻射定標(biāo)、多視化處理和地理編碼。在GAMMA軟件中，首先讀取Level 1.1數(shù)據(jù)，生成的單視復(fù)數(shù)圖像(SLC)再進(jìn)行輻射定標(biāo)，JAXA提供的定標(biāo)常量為-115 dB[5]，將強(qiáng)度圖數(shù)字量化輸出的DN值轉(zhuǎn)換為后向散射系數(shù)圖，便于對(duì)圖像進(jìn)行點(diǎn)目標(biāo)定量度量和分析。其次對(duì)數(shù)據(jù)進(jìn)行多視處理，距離向和方位向視數(shù)分別為1和7，生成強(qiáng)度圖。最后采用數(shù)字高程模型(DEM)數(shù)據(jù)進(jìn)行地理編碼，是將SAR圖像的斜距距離/方位坐標(biāo)變換成地圖幾何坐標(biāo)，處理后得到PALSAR影像。

在ArcGIS中利用2007年宜良縣森林資源二類調(diào)查數(shù)據(jù)，選出宜良縣花園林場(chǎng)中不同齡組的云南松，本文主要提取云南松幼齡林、中齡林以及近熟林的小班數(shù)據(jù)，其中相同齡組的云南松郁閉度和坡度基本一致。

2 結(jié)果與分析

2.1 云南松全極化后向散射特征分析

根據(jù)2007年宜良縣森林資源二類調(diào)查數(shù)據(jù)采集研究區(qū)云南松幼齡林平均樹(shù)高為4 m，平均胸徑為8 cm的255個(gè)樣本數(shù)據(jù)；云南松中齡林平均樹(shù)高為6 m，平均胸徑為12 cm的255個(gè)樣本數(shù)據(jù)以及云南松近熟林平均樹(shù)高為11 m，平均胸徑為18 cm的255個(gè)樣本數(shù)據(jù)。然后分別提取影像中云南松不同齡組的云南松后向散射系數(shù)，并統(tǒng)計(jì)云南松不同極化方式下不同齡組的云南松平均后向散射系數(shù)，見(jiàn)表1。

從表1可以看出，無(wú)論云南松屬于哪種齡組，在HH極化方式下反射強(qiáng)度值最高，其次是VV極化方式，而HV和VH極化方式下反射值都比較低。云南松在幼齡林到中齡林時(shí)期，處于生長(zhǎng)期，冠層增大，對(duì)各極化波的作用表現(xiàn)明顯。因此，在這期間，各齡組云南松平均后向散射系數(shù)都是增強(qiáng)趨勢(shì)。HH極化方式下反射強(qiáng)度值最強(qiáng)，因?yàn)镠H極化方式反應(yīng)的大多是植被冠層葉片的信息，而此時(shí)期云南松冠層增大，含水量較高；同時(shí)云南松形態(tài)變大又削弱了對(duì)VV極化方式下的穿透深度。而從中齡林到近熟林時(shí)期，云南松冠層針葉掉落，形態(tài)變小，因而云南松平均后向散射系數(shù)減小。

表1 云南松各齡組平均后向散射系數(shù)

2.2 云南松全極化后向散射特征方差分析

利用SPSS軟件實(shí)現(xiàn)對(duì)不同極化方式下各齡組云南松255個(gè)樣本進(jìn)行單因素方差分析(One way ANOVA)，并采用Least-signifcant difference(LSD)方法對(duì)各組均值間的變量進(jìn)行多重比較，確定云南松各齡組與后向散射系數(shù)是否顯著。由表2可知，顯著水平小于0.0003，表明云南松各齡組與后向散射系數(shù)之間具有極顯著影響，意味著根據(jù)后向散射系數(shù)可以區(qū)分不同齡組的云南松。通過(guò)LSD方法對(duì)云南松后向散射系數(shù)各組均值間的變量進(jìn)行多重比較，其中因變量后向散射系數(shù)會(huì)隨云南松齡組不同而變化，因此進(jìn)一步確定經(jīng)過(guò)單因素方差分析可以區(qū)分不同齡組的云南松，見(jiàn)表3。

表2 云南松各齡組后向散射系數(shù)的單因素方差分析結(jié)果

表3 云南松各齡組后向散射系數(shù)LSD多重比較

*：(I)齡組中1為云南松幼齡林；2為云南松中齡林；3為云南松近熟林；*為均值差的顯著性水平為0.05。

LSD法多重比較法中，(I)齡組為比較基準(zhǔn)齡組，(J)齡組為比較齡組。通過(guò)LSD法多重比較表明：齡組1與齡組2之間存在顯著性差異；齡組2與齡組1和齡組3之間存在顯著性差異；齡組3與齡組2之間存在顯著性差異；齡組1即是云南松幼齡林，齡組2即是云南松中齡林，齡組3即是云南松近熟林。因此，根據(jù)云南松的后向散射系數(shù)可以把不同齡組的云南松區(qū)分開(kāi)來(lái)。

3 結(jié)論與討論

本文通過(guò)研究宜良縣花園林場(chǎng)地區(qū)天然林——云南松不同齡組在不同極化方式下L波段后向散射系數(shù)相關(guān)分析，同時(shí)進(jìn)行單因素方差分析。結(jié)果表明，在云南松幼齡林、中齡林和近熟林中，HH極化方式下反射強(qiáng)度值都是最高的，其次是VV極化方式，而HV和VH極化方式下反射值都比較低；通過(guò)對(duì)云南松全極化后向散射特征進(jìn)行方差分析，可知云南松各齡組與后向散射系數(shù)之間具有極顯著影響，根據(jù)其顯著性，可以把云南松幼齡林、中齡林和近熟林區(qū)分開(kāi)來(lái)。

由于電磁波在經(jīng)過(guò)植被層后，其回波的能量會(huì)有所改變，而且同一地物對(duì)于不同的極化方式結(jié)果也不同，反映在雷達(dá)圖像上就是后向散射系數(shù)的差異。一般認(rèn)為HV極化方式下為植被冠層信息，HH極化方式下為近地表信息[6]。粗壯樹(shù)木的樹(shù)干與地面都可能發(fā)生偶次散射，特別是對(duì)HH、VV分量回波很強(qiáng)，所以HH、VV的后向散射系數(shù)相對(duì)較大。而HV和VH相比HH和VV的滲透能力弱，回波信號(hào)強(qiáng)度低，所以HV、VH的后向散射系數(shù)相對(duì)較小。但HV圖像能給出地物目標(biāo)較細(xì)微的結(jié)構(gòu)，地物界限更分明。對(duì)樹(shù)木的種類、HV圖像有較好的識(shí)別[7]。

[1]吳興華.云南松次生林生長(zhǎng)研究[J].林業(yè)建設(shè)，2010(3)：47-50.

[2]王曉蘭，戴博偉，楊汝良.全極化合成孔徑雷達(dá)多視圖像的極化特征分析[J].電子與信息學(xué)報(bào)，2002(4)：507-513.

[3]Sagues L，Lopez-Sanchez J M，et al.Polarimetric radar interferometry for improved mine detection and surface clutter rejection[J].IEEE Transactions on Geosciences and Remote Sensing，2001，39(6)：1271-1278.

[4]云南省林業(yè)調(diào)查規(guī)劃院.2007年宜良縣森林資源二類調(diào)查報(bào)告[Z].昆明，2007.

[5]SHIMADA M，ISOGUCHI O，TADONO T.PLASAR radiometric calibration and geometric calibration[J].IEEE Transactions on Geoscience and Remote Sensing，2009，47(12)：3915-3932.

[6]張慶云，劉國(guó)林，陶秋香，等.基于極化干涉SAR的植被高度反演算法對(duì)比分析[J].山東科技大學(xué)學(xué)報(bào)，2014(6)：77-83.

[7]郭華東.機(jī)載雷達(dá)遙感應(yīng)用實(shí)驗(yàn)研究[M].北京：科學(xué)出版社，1992：81-85.

Analysis on Backscattering Characteristics of thePinusyunnanensisBased on Full-polarization SAR Data

XIAO Hong-yan，TIAN Chuan-zhao，YUE Cai-rong

(ForestryCollege，SouthwestForestryUniversity，Kunming650224，Yunnan，China)

Using the garden forest farm in the Yiliang County of Kunming as study area，the full-polarization single complex type of L-band data is ALOS PALSAR Level 1.1 on April 23，2009.The backscattering characteristics of thePinusyunnanensis analyzed. Using SPSS software conducts experimental results by one-way ANOVA statistical analysis.The results indicates that the backscattering characteristics of thePinusyunnanensisand different ages(young，middle-aged and over mature forests) have a significant effect.According to these characteristics，thePinusyunnanensiscan be distinguished from different ages ground.

full-polarization SAR；Pinusyunnanensis；backscattering coefficient；one-way analysis of variance

2014-09-24；

2014-10-23

國(guó)家自然基金(31260156)

肖虹雁(1989—)，女，福建三明人，西南林業(yè)大學(xué)林學(xué)院碩士研究生，從事資源與環(huán)境遙感研究。E-mail：742929194@qq.com。

岳彩榮(1964—)，男，云南建水人，西南林業(yè)大學(xué)林學(xué)院教授、博士、博士生導(dǎo)師，從事遙感與GIS的教學(xué)與研究。E-mail：cryue@163.com。

10.13428/j.cnki.fjlk.2015.03.018

S791.257

A

1002-7351(2015)03-0085-03