王 平，周 忠 發(fā)*，廖 娟

(貴州師范大學(xué)喀斯特研究院，貴州 貴陽(yáng) 550001；貴州省遙感中心，貴州 貴陽(yáng) 550001)

基于Freeman分解的喀斯特高原山區(qū)煙田土壤水分反演研究

王 平，周 忠 發(fā)*，廖 娟

(貴州師范大學(xué)喀斯特研究院，貴州 貴陽(yáng) 550001；貴州省遙感中心，貴州 貴陽(yáng) 550001)

為探究SAR技術(shù)在喀斯特高原山區(qū)煙田旱情監(jiān)測(cè)，選取貴州省清鎮(zhèn)市流長(zhǎng)鄉(xiāng)現(xiàn)代烤煙農(nóng)業(yè)基地單元為研究區(qū)，采用SAR目標(biāo)極化分解技術(shù)反演團(tuán)棵期煙田土壤水分。選取Radarsat-2全極化數(shù)據(jù)，通過(guò)Freeman-Durden極化目標(biāo)分解技術(shù)，分解出3種散射機(jī)制，并從目視解譯假彩色合成影像與散射功率兩個(gè)角度綜合分析主要散射機(jī)制，根據(jù)SAR影像散射量中的介電常數(shù)與土壤水分密切相關(guān)的原理，嘗試不同的回歸模型，用主要散射機(jī)制占比反演煙田土壤水分，并用線(xiàn)性回歸檢驗(yàn)法檢驗(yàn)其反演精度。結(jié)果表明：此方法能夠較好地反演喀斯特高原山區(qū)煙田團(tuán)棵期土壤水分，體現(xiàn)出SAR影像信息中單次散射模擬微粗糙表面散射及反演土壤水分的優(yōu)勢(shì)，為喀斯特高原山區(qū)煙田旱情監(jiān)測(cè)提供技術(shù)參考。

喀斯特高原山區(qū)；土壤水分反演；Freeman-Durden分解；Radarsat-2；煙田

0 引言

地表土壤水分的時(shí)空分布信息在水文學(xué)、氣候?qū)W、植物生態(tài)學(xué)等學(xué)科中有著重要地位，同時(shí)也在預(yù)報(bào)旱情、估算農(nóng)作物產(chǎn)量等方面起到重要作用[1]?？緹熓俏覈?guó)重要的經(jīng)濟(jì)作物之一，提高烤煙種植的精準(zhǔn)度，已成為現(xiàn)代烤煙農(nóng)業(yè)發(fā)展的需要[2]。土壤水分是烤煙生長(zhǎng)發(fā)育過(guò)程中生理需水和生態(tài)需水的主要來(lái)源[3]。由于傳統(tǒng)基于點(diǎn)的水分測(cè)量方法只能觀測(cè)稀疏空間點(diǎn)的土壤水分信息，不能全面反映不同尺度上的水量分布和變化規(guī)律[4]，無(wú)法滿(mǎn)足現(xiàn)代山地高效農(nóng)業(yè)大范圍、實(shí)時(shí)、精準(zhǔn)監(jiān)測(cè)的需求。目前，已有眾多學(xué)者通過(guò)各類(lèi)衛(wèi)星遙感反演地表參數(shù)，包括光學(xué)遙感、微波遙感及多傳感器聯(lián)合反演算法，并建立相應(yīng)的經(jīng)驗(yàn)或半經(jīng)驗(yàn)?zāi)Ｐ蚚5-10]。諸多研究表明，C波段(約5 GHz)雷達(dá)系統(tǒng)估計(jì)0～5 cm厚的裸露地表土壤水分可達(dá)到較好的精度[11]，微波遙感對(duì)土壤水分的高敏感性，已在土壤水分研究中發(fā)揮了重要作用。目前在軌運(yùn)行的Radarsat-2提供的多極化、多角度的C波段SAR數(shù)據(jù)對(duì)于地表土壤水分監(jiān)測(cè)具有很大潛力。

貴州喀斯特地區(qū)經(jīng)常面臨旱災(zāi)，復(fù)雜的地表與氣候環(huán)境對(duì)該區(qū)域的土壤水分監(jiān)測(cè)帶來(lái)很大的困難，多云雨的氣候使得運(yùn)用傳統(tǒng)的可見(jiàn)光成像難以完整地獲取地物信息。星載SAR數(shù)據(jù)具有全天時(shí)、全天候、分辨率高、信息豐富、穿透能力強(qiáng)等成像優(yōu)勢(shì)，使監(jiān)測(cè)更加及時(shí)[12,13]。SAR回波信號(hào)直接受到地物復(fù)介電常數(shù)、地表粗糙度及植被覆蓋影響，而復(fù)介電常數(shù)與含水量直接相關(guān)，因此可根據(jù)SAR回波信號(hào)分析土壤含水量。Freeman-Durden模型匹配方法是在植被覆蓋情況下，基于雷達(dá)散射回波的物理模型，而不是單純的數(shù)學(xué)推導(dǎo)，避免了由于各種物理?xiàng)l件不同造成的反演模型不合適問(wèn)題[14,15]。因此，通過(guò)Freeman-Durden極化分解技術(shù)，分析出煙田團(tuán)棵期的主導(dǎo)分解分量，并采用其分解分量占比反演土壤水分，從而探索Radarsat-2數(shù)據(jù)在喀斯特高原山區(qū)對(duì)煙田團(tuán)棵期土壤水分的反演能力。

1 研究區(qū)概況及數(shù)據(jù)處理

1.1 研究區(qū)概況

研究區(qū)地處喀斯特高原山區(qū)貴州清鎮(zhèn)流長(zhǎng)現(xiàn)代烤煙農(nóng)業(yè)基地單元(位于106°7′6″～106°29′37″E、26°24′5″～26°45′45″N)，轄流長(zhǎng)、犁倭、紅楓湖3個(gè)鄉(xiāng)(鎮(zhèn))，總面積為489 km2，所選研究面積為20.8 km2，包括4個(gè)自然村。該區(qū)域?yàn)榈湫偷目λ固氐孛玻苑鍏餐莸?、峰叢谷地為主，海拔?00～1 400 m之間，相對(duì)高差較大，地下暗河和落水洞廣布，地表崎嶇破碎并伴隨有石漠化問(wèn)題，導(dǎo)致煙田分布較分散且地塊破碎[16]；該區(qū)域?qū)賮啛釒Ц咴撅L(fēng)濕潤(rùn)氣候區(qū)，年平均氣溫14℃，年均降水量1 150.4 mm，無(wú)霜期290 d左右，全年陰雨天多，日照少，導(dǎo)致光學(xué)遙感數(shù)據(jù)難以獲取。土壤類(lèi)型以黃沙壤、黃壤為主，pH值為5.5～6.5，有機(jī)質(zhì)含量豐富，含鹽和氯較低，土壤質(zhì)地較為黏重但排水良好，較適宜烤煙生長(zhǎng)，因此區(qū)域內(nèi)常年種植烤煙，主要品種為云煙85、云煙87、K326及南江三號(hào)等。近幾年隨著研究區(qū)烤煙種植規(guī)模的逐年遞增，其種植面積占全國(guó)烤煙種植面積的15%左右，居全國(guó)第二[1，17]，種植效益明顯提高。

1.2 數(shù)據(jù)選取

研究采用的數(shù)據(jù)包括遙感數(shù)據(jù)和實(shí)測(cè)數(shù)據(jù)，整個(gè)監(jiān)測(cè)時(shí)期與研究區(qū)內(nèi)烤煙的生長(zhǎng)期基本同步。微波遙感數(shù)據(jù)選擇2014年5月29日Radarsat-2全極化SLC數(shù)據(jù)，標(biāo)稱(chēng)空間分辨率為8 m(圖1，見(jiàn)封3)，圖1中的區(qū)域是根據(jù)研究區(qū)域邊界對(duì)SAR影像裁剪的結(jié)果?？緹焾F(tuán)棵期土壤水分直接影響著烤煙移栽成活率，對(duì)指導(dǎo)煙民田間管理具有重要意義，因此獲取研究區(qū)煙田團(tuán)棵期Radarsat-2數(shù)據(jù),并在同一時(shí)期獲取研究區(qū)0.5 m分辨率航拍圖，用于輔助解譯地面信息。

實(shí)測(cè)數(shù)據(jù)通過(guò)建立野外訓(xùn)練樣方與采集煙田土樣實(shí)現(xiàn)。由于雷達(dá)影像的空間分辨率為8 m，因此在研究區(qū)內(nèi)建立了20個(gè)面積為16 m*16 m的試驗(yàn)樣方(圖1，表1)，圖1中綠色點(diǎn)位置表示樣方所在地。將所有樣方平均分為實(shí)驗(yàn)組、檢驗(yàn)組兩組，樣方標(biāo)準(zhǔn)一致且均勻分布于研究區(qū)域內(nèi)。由于研究區(qū)地形復(fù)雜、地塊破碎，因此在每個(gè)煙田樣方的4角點(diǎn)處用GPS定位，定點(diǎn)順序由樣方4個(gè)角點(diǎn)中任意一個(gè)為起點(diǎn)，順時(shí)針采集GPS定位數(shù)據(jù)。依據(jù)樣方的經(jīng)緯度信息，在Google Earth軟件中解譯出樣方所在地，并展示在SAR影像中(圖1)。為保證實(shí)測(cè)數(shù)據(jù)與衛(wèi)星遙感數(shù)據(jù)的獲取實(shí)現(xiàn)同步或準(zhǔn)同步，在獲取SAR數(shù)據(jù)的當(dāng)天，進(jìn)行煙田土樣采集。選用規(guī)格為5 cm的鋁盒為樣品存儲(chǔ)容器，在樣方內(nèi)沿“S”形路線(xiàn)采集土樣。最后，需要給煙田樣地拍照，詳細(xì)記錄煙田狀況，從而為反演結(jié)果的分析提供依據(jù)。

表1 2014年實(shí)驗(yàn)樣方基本情況Table 1 The basic information of test sample plots in 2014

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

1.3.1 SAR影像Freeman-Durden三分量分解 數(shù)據(jù)處理包括SAR影像的預(yù)處理及實(shí)測(cè)數(shù)據(jù)的處理。全極化雷達(dá)通過(guò)同時(shí)發(fā)射和接收水平與垂直極化電磁波來(lái)獲取地物信息，目標(biāo)極化分解理論是分析與提取極化信息的一種方法[18],可將復(fù)雜地物回波信號(hào)分解為多個(gè)簡(jiǎn)單的標(biāo)準(zhǔn)目標(biāo)回波信號(hào)之和，通過(guò)分析不同散射機(jī)制下的各分量獲得原目標(biāo)信息[19]。由于極化分解技術(shù)與獲取目標(biāo)信號(hào)中的相位有關(guān)，F(xiàn)reeman-Durden目標(biāo)分解可在一定程度上抑制SAR系統(tǒng)產(chǎn)生的相干噪聲，因此為了保持原始數(shù)據(jù)信息的完整性，未對(duì)影像進(jìn)行濾波處理。

極化目標(biāo)分解技術(shù)包括相干目標(biāo)分解和非相干目標(biāo)分解[20]，目前針對(duì)植被覆蓋環(huán)境，F(xiàn)reeman-Durden分解作為非相干分解，將復(fù)雜的地物極化信息簡(jiǎn)化表達(dá)為3種分量，分別通過(guò)對(duì)3種散射機(jī)制進(jìn)行建模，形成由隨機(jī)取向偶極子組成的云狀管層體散射，由一對(duì)不同介電常數(shù)的正交平面構(gòu)成的二次散射，以及適度粗糙表面的布拉格(Bragg)散射，即來(lái)自地表的單次散射、地表與植被層共同構(gòu)成的二次散射以及來(lái)自植被冠層構(gòu)成的體散射，且可用來(lái)初步確定在極化SAR數(shù)據(jù)的后向散射中占主導(dǎo)地位的散射機(jī)制。利用Freeman-Durden目標(biāo)分解方法處理雷達(dá)數(shù)據(jù)，將回波信號(hào)簡(jiǎn)化為3種分量，在研究中認(rèn)為同時(shí)期煙田植被覆蓋與地表粗糙度條件相同，模型中的二次散射能夠表達(dá)出植被與土壤表面結(jié)構(gòu)的共同作用。假設(shè)在3種散射機(jī)制成分互不相干的情況下，總的后向散射模型為三者之和(式(1))[15，20]。

C=Cs+Cd+Cv

(1)

(2)

(3)

(4)

式中：Cs、Cd與Cv分別為單次散射協(xié)方差矩陣、二次散射協(xié)方差矩陣與體散射平均協(xié)方差矩陣；fs、fd與fv分別為單次、二次和體散射分量；α為產(chǎn)生二次散射的兩平面復(fù)介電常數(shù)之比相關(guān)參數(shù)，β為產(chǎn)生單次散射平面的地表相關(guān)參數(shù)。

1.3.2 實(shí)測(cè)數(shù)據(jù)處理 主要是對(duì)野外采集的土樣處理。先將每個(gè)鋁盒(5cm規(guī)格)貼上標(biāo)簽,將裝有土樣的鋁盒稱(chēng)重后放進(jìn)烘土箱，設(shè)置100℃;待土壤干透后開(kāi)始稱(chēng)重記錄鋁盒與干土總重，并稱(chēng)取烘干后鋁盒重量。土壤水分的計(jì)算采用土壤重量含水量計(jì)算方法(式(5))，通過(guò)計(jì)算每個(gè)煙田樣方中土樣重量含水量的平均值得出各樣方的土壤重量含水量。通過(guò)研究區(qū)氣象站實(shí)時(shí)獲取研究區(qū)2014年5月29日的土壤墑情數(shù)據(jù)(表2)，以便驗(yàn)證反演模型的可靠性。

(5)

式中:θg為實(shí)測(cè)計(jì)算出的土壤重量含水量，MW為土壤總含水重量，MS為土壤干重，M1為烘干前鋁盒及土壤質(zhì)量，M2為烘干后鋁盒及土樣質(zhì)量，MC為烘干后鋁盒重量。

表2 研究區(qū)土壤墑情匯總Table 2 Summary of soil moisture in the study area

2 結(jié)果與分析

2.1 Freeman-Durden分解結(jié)果

應(yīng)用極化分解技術(shù)，對(duì)研究區(qū)頻率為5 GHz的C波段雷達(dá)數(shù)據(jù)應(yīng)用式(2)進(jìn)行極化分解。圖2中，(a)(b)(c)三幅圖分別是應(yīng)用目標(biāo)信息的Freeman-Durden三分量分解技術(shù)得到的研究區(qū)域內(nèi)單次散射、二次散射與體散射的后向散射系數(shù)圖。

圖2 Freeman Durden目標(biāo)分解結(jié)果Fig.2 Freeman Durden target decomposition results

將3種散射機(jī)制結(jié)果用偽彩色合成顯示(圖3，見(jiàn)封3)，其中紅色對(duì)應(yīng)二次散射，綠色對(duì)應(yīng)體散射，藍(lán)色對(duì)應(yīng)單次散射，從合成圖(圖3)中目視解譯出在煙田集中的地區(qū)，藍(lán)色、紅色及紫色(藍(lán)色與紅色疊加效果)明顯分布較多。根據(jù)樣方分布位置，裁剪出煙田樣方，并計(jì)算得到每個(gè)樣方中3種散射機(jī)制的散射功率占比(圖4)。通過(guò)定量分析可以得出，在團(tuán)棵期煙田雷達(dá)影像中，單次散射回波最強(qiáng)，其次為二次散射，體散射最弱。

圖4 三種散射機(jī)制散射功率占比Fig.4 The proportion of scattered power on the three kinds of scattering mechanisms

基于以上散射機(jī)制分布結(jié)果，進(jìn)一步分析其形成原因。由于所獲得的雷達(dá)數(shù)據(jù)成像于煙田團(tuán)棵期，根據(jù)實(shí)地考察及照片的記錄可知，在此時(shí)期煙田內(nèi)主要以煙苗及裸露地表為主；由于SAR數(shù)據(jù)中的C波段波長(zhǎng)較長(zhǎng)，具有較好的穿透性能，部分雷達(dá)信號(hào)可穿透植被層，產(chǎn)生地表的后向回波，因此初步分析出在雷達(dá)后向回波中，地表的后向回波及植被與地表共同作用產(chǎn)生的回波信號(hào)較強(qiáng)，即構(gòu)成單次散射項(xiàng)的直接地表散射和地表與烤煙層共同構(gòu)成的二次散射較強(qiáng)。又由于此時(shí)期烤煙葉片很小，葉長(zhǎng)在13～35 cm之間，葉寬在6～18 cm之間，葉面積指數(shù)(LAI)在0.16～0.22之間，地表散射信息占主導(dǎo)，導(dǎo)致分解的雷達(dá)后向散射總量中單次散射較二次散射大。由于體散射一般主要由森林、部分建筑物及其他結(jié)構(gòu)復(fù)雜的目標(biāo)形成，因此體散射在煙田內(nèi)的散射機(jī)制中所占比例較小。

2.2 反演模型建立

根據(jù)Freeman-Durden三分量分解結(jié)果可知，單次散射在研究區(qū)內(nèi)的散射機(jī)制中占主導(dǎo)。由于植被覆蓋地區(qū)的單次散射直接反映的是地表土壤的散射量，地表散射量中包含的土壤介電常數(shù)信息又與土壤水分密切相關(guān)，因此，嘗試多種回歸方程的擬合，建立實(shí)測(cè)樣方土壤水分與分解的單次散射占比之間的關(guān)系(表3)，從表中可以看出，三次函數(shù)模型相關(guān)性系數(shù)R2最高，達(dá)到0.7893。從三次函數(shù)模型分析結(jié)果(圖5)中可以看出兩者相關(guān)性很好，充分證明了利用土壤介電常數(shù)信息作為媒介，將SAR影像中地表散射回波信息與地表土壤水分建立關(guān)系，能夠很好地解釋土壤水分。

表3 土壤重量含水量與單次散射占比回歸分析Table 3 Regression analysis of soil gravimetric water content and single scattering proportion

圖5 三次函數(shù)反演模型Fig.5 Three function regression model

2.3 反演模型精度分析

運(yùn)用線(xiàn)性回歸檢驗(yàn)法，對(duì)以上模型進(jìn)行精度驗(yàn)證，選取研究區(qū)內(nèi)10個(gè)檢驗(yàn)組土壤水分作為模型精度的驗(yàn)證數(shù)據(jù)，根據(jù)Freeman-Durden三分量分解技術(shù)將SAR影像中檢驗(yàn)組煙田的單次散射占比帶入三次函數(shù)模型計(jì)算出土壤水分反演值，圖6中點(diǎn)表示每個(gè)土壤水分反演值對(duì)應(yīng)的土壤水分實(shí)測(cè)值，虛線(xiàn)是模型的回歸結(jié)果，實(shí)線(xiàn)表示“1∶1”線(xiàn)。

圖6 土壤水分反演值與實(shí)測(cè)值關(guān)系Fig.6 Relationship between soil moisture inversion and actual measured value

土壤水分實(shí)測(cè)值與反演值間的擬合度R2=0.7881，模擬直線(xiàn)在y軸的截距為0.0006，由圖5可知，10個(gè)隨機(jī)檢驗(yàn)樣本中的個(gè)別檢驗(yàn)樣本數(shù)據(jù)與土樣的整體趨勢(shì)有偏差，煙田樣方個(gè)別反演值與實(shí)測(cè)值偏差較大。根據(jù)檢驗(yàn)組樣方的記錄，主要是在樣方選取過(guò)程中雖已考慮煙田標(biāo)準(zhǔn)一致問(wèn)題，但在喀斯特地區(qū)，地塊破碎，烤煙種植的壟距與株距存在一定差距，極化分解出的三種散射機(jī)制間存在細(xì)微且不同程度的相干影響，導(dǎo)致極化分解的單次散射強(qiáng)度有差異，因此擬合線(xiàn)與“1∶1”線(xiàn)之間存在偏移，但整體基本滿(mǎn)足精度要求，可以認(rèn)為利用Freeman-Durden三分量分解獲得的單次散射機(jī)制占比反演土壤水分的結(jié)果與實(shí)測(cè)結(jié)果擬合較好，反演精度較高，能較好地反映煙田土壤水分狀況。

3 結(jié)論與討論

通過(guò)對(duì)喀斯特高原山區(qū)煙田SAR影像進(jìn)行Freeman-Durden目標(biāo)極化分解，分析出主導(dǎo)散射機(jī)制，并進(jìn)一步用主導(dǎo)散射機(jī)制占比反演煙田土壤水分。得出以下結(jié)論：1)通過(guò)采用極化目標(biāo)分解技術(shù)，根據(jù)散射相關(guān)矩陣特征值，對(duì)SAR影像中的煙田進(jìn)行Freeman-Durden三分量分解，可將混合復(fù)雜的煙田散射回波簡(jiǎn)化。2)三分量分解公式充分表達(dá)了地表、地表與植被、植被3種地物結(jié)構(gòu)，并能夠分析出其主導(dǎo)散射機(jī)制為單次散射，用主導(dǎo)散射機(jī)制占比能較精確地反演煙田土壤水分。與通過(guò)數(shù)學(xué)建模過(guò)程反演地表土壤水分相比，不需要任何地面測(cè)量數(shù)據(jù)，直接從雷達(dá)散射回波的物理基礎(chǔ)分析地表覆蓋狀態(tài)，避免了反演模型的不確定性。3)嘗試用不同回歸模型建立反演關(guān)系，最終選擇擬合度最高的三次函數(shù)模型，并使用線(xiàn)性回歸檢驗(yàn)法驗(yàn)證模型精度，驗(yàn)證模型的擬合度達(dá)到0.7，反演結(jié)果與2014年研究區(qū)氣象站監(jiān)測(cè)結(jié)果相符，證實(shí)該方法的反演精度較好。4)首次在貴州喀斯特高原山區(qū)運(yùn)用SAR目標(biāo)極化分解技術(shù)進(jìn)行煙田土壤水分的反演實(shí)驗(yàn)，充分體現(xiàn)出SAR目標(biāo)極化分解的單次散射模擬喀斯特高原山區(qū)微粗糙表面散射的優(yōu)勢(shì)。

由于喀斯特高原山區(qū)煙田在旺盛期與成熟期地表復(fù)雜，因此本研究未對(duì)其他時(shí)期的煙田進(jìn)行試驗(yàn)，下一步將結(jié)合多時(shí)相SAR影像探究合適的極化目標(biāo)分解方法，從而更好地反映不同時(shí)期的煙田旱情，推進(jìn)SAR目標(biāo)極化分解技術(shù)在喀斯特高原山區(qū)反演煙田土壤水分的研究，從而滿(mǎn)足現(xiàn)代山地高效農(nóng)業(yè)大范圍、實(shí)時(shí)、精準(zhǔn)監(jiān)測(cè)的需求。

[1] 周忠發(fā)，李波，賈龍浩.合成孔徑雷達(dá)技術(shù)在喀斯特山區(qū)煙草種植定量監(jiān)測(cè)應(yīng)用探討[J].測(cè)繪通報(bào)，2012(S1)：246-248.

[2] 中國(guó)農(nóng)業(yè)科學(xué)院烤煙研究所.中國(guó)烤煙栽培學(xué)[M].上海：上?？茖W(xué)技術(shù)出版社，1987.121.

[3] 田國(guó)良.土壤水分的遙感監(jiān)測(cè)方法[J].環(huán)境遙感,1991(2)：89-98.

[4] 隋立春.主動(dòng)式雷達(dá)遙感[M].北京：測(cè)繪出版社，2009.153-154.

[5] 奧利弗.合成孔徑雷達(dá)圖像理解[M].北京：電子工業(yè)出版社，2009.273-293.

[6] 郭華東.中國(guó)雷達(dá)遙感圖像分析[M].北京：科學(xué)出版社，1999.1-35.

[7] 李震，廖靜娟.合成孔徑雷達(dá)地表參數(shù)反演模型與方法[M].北京：科學(xué)出版社，2011.221-255.

[8] 賈明權(quán).水稻微波散射特性研究及參數(shù)反演[D].成都：電子科技大學(xué)，2013.55-59.

[9] 楊虎，施建成.FY-3微波成像儀地表參數(shù)反演研究[J].遙感技術(shù)與應(yīng)用,2005(1)：194-200.

[10] 楊沈斌，李秉柏，申雙和,等.基于ENVISAT ASAR數(shù)據(jù)的水稻遙感監(jiān)測(cè)[J].江蘇農(nóng)業(yè)學(xué)報(bào),2008，24(1)：33-38.

[11] 任鑫.多極化、多角度SAR土壤水分反演算法研究[D].北京：中科院遙感應(yīng)用研究所,2003.16-20.

[12] 張艷寧，李映.SAR圖像處理的關(guān)鍵技術(shù)[M].北京：電子工業(yè)出版社，2014.340-347.

[13] 韓桂紅.干旱區(qū)鹽漬地極化雷達(dá)土壤水分反演研究[D].烏魯木齊：新疆大學(xué)，2013.10-13.

[14] 劉偉，施建成，王建明.極化分解技術(shù)在估算植被覆蓋地區(qū)土壤水分變化中的應(yīng)用[J].遙感信息，2005(4)：3-6.

[15] LEE J S,POTTIER E.極化雷達(dá)成像基礎(chǔ)與應(yīng)用[M].北京：電子工業(yè)出版社，2013.133-153.

[16] 高貴龍，鄧自民，熊康寧，等.喀斯特的呼喚與希望[M].貴陽(yáng)：貴州科技出版社，2003.8-20.

[17] 賈龍浩，周忠發(fā)，李波.高分辨率SAR在喀斯特山地?zé)煵萆L(zhǎng)建模中的應(yīng)用探討[J].中國(guó)烤煙科學(xué)，2013，34(5)：104-112.

[18] 徐星歐，舒寧，李磊.Freeman-Durden分解方法提取雷達(dá)影像含水量信息[J].計(jì)算機(jī)工程與應(yīng)用，2010，46(27)：17-19，24.

[19] 王慶，曾琪明，廖靜娟.基于特征向量分解與散射機(jī)制判別指數(shù)的全極化SAR圖像地物提取與分類(lèi)[J].遙感信息,2012(2)：9-14.

[20] FREEMAN A,DURDEN S L.A three-component scattering modle for polarrimetic SAR Data[J].IEEE Transactions on Geoscience and Remote Sensing，1998，36(3)：963-973.

Study on Soil Moisture Retrieval of Tobacco Field in Karst Plateau Mountainous Area Based on Freeman Decomposition

WANG Ping,ZHOU Zhong-fa,LIAO Juan

(InstituteofKarst,GuizhouNormalUniversity，Guiyang550001；NationalRemoteSensingCenterGuizhouBranch，Guiyang550001，China)

Research purposes is to explore drought monitoring with SAR technology in karst plateau mountainous area.The modern tobacco agriculture base in Liuchang towns of Qingzhen City as study area,which is located in Guizhou Province.In this paper,Radarsat-2 full polarization data are selected,the soil moisture content of the tobacco field at rosette stage is retrieved by the SAR polarization decomposition technique.Three kinds of scattering mechanisms are obtained by means of Freeman-Durden polarimetric target decomposition technique.And it synthetically analyzes the main scattering mechanisms from two angles:visual interpretation of false color composite image and the scattering power.According to the principle of the dielectric constant in the scattering parameters of the SAR image and the soil moisture is closely related,different regression models are tried to inversed soil moisture in tobacco field by the main scattering mechanism.The inversion accuracy is tested by linear regression test method.The results shows that this method can be used to inverse the soil moisture content in the tobacco field in the mountain area of karst plateau.The advantages of single scattering in SAR images for simulating the scattering of micro rough surfaces and inversed soil moisture are also reflected.This method provides a technical reference for the tobacco field drought monitoring in karst plateau mountainous area.

karst plateau mountainous area;soil moisture retrieval;Freeman-Durden decomposition;Radarsat-2; tobacco field

2015-08-08；

2015-10-29

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)課題“人為干預(yù)下喀斯特山地石漠化的演變機(jī)制與調(diào)控”(2012CB723202)；貴州省科技計(jì)劃“喀斯特山區(qū)SAR遙感平臺(tái)監(jiān)測(cè)與識(shí)別關(guān)鍵技術(shù)研究與應(yīng)用”(黔科合GY字〔2013〕3062)；貴州省重大應(yīng)用基礎(chǔ)研究項(xiàng)目“喀斯特石漠化生態(tài)修復(fù)及生態(tài)經(jīng)濟(jì)系統(tǒng)優(yōu)化調(diào)控研究-巖土類(lèi)型格局”(黔科合JZ字[2014]200201)

王平(1991-)，女，碩士研究生，研究方向?yàn)榈乩硇畔⑾到y(tǒng)與遙感。*通訊作者E-mail：fa6897@163.com

10.3969/j.issn.1672-0504.2016.02.014

TP79;S152.7

A

1672-0504(2016)02-0072-05