唐 彥

(東北林業(yè)大學(xué),黑龍江 哈爾濱 150040)

土壤含鹽量反演的研究

唐 彥

(東北林業(yè)大學(xué),黑龍江 哈爾濱 150040)

運(yùn)用 Hyperion數(shù)據(jù),以黑龍江省大慶市某一實(shí)驗(yàn)區(qū)為例,開展對土壤含鹽量定量提取的研究,通過對圖像預(yù)處理、特征提取、建立BP神經(jīng)網(wǎng)絡(luò)模型(Back Propagation Network)等研究工作,探討反演土壤含鹽量的方法。研究結(jié)果表明:神經(jīng)網(wǎng)絡(luò)模型具有極強(qiáng)的線性和非線性擬合能力,模擬遙感影像特征與土壤鹽分之間比較復(fù)雜的關(guān)系上有很大優(yōu)勢。研究結(jié)果不但為利用 Hyperion數(shù)據(jù)反演土壤含鹽量提供理論依據(jù),而且還為其它地表參數(shù)的反演提供參考。

土壤含鹽量;Hyperion數(shù)據(jù);反演;BP神經(jīng)網(wǎng)絡(luò)模型

土地鹽堿化是土壤退化的一種重要表現(xiàn),大范圍的實(shí)地調(diào)查研究需要投入大量的人力、物力和時間[1],采用遙感技術(shù)對土壤鹽漬化程度進(jìn)行分析,對鹽堿地的綜合治理尤其重要,國內(nèi)外對土壤的鹽漬化問題做了大量的研究[2]。高光譜數(shù)據(jù)提供了連續(xù)窄帶光譜信息,為土壤評價與監(jiān)測提供了強(qiáng)有力的工具,Dehaan采用高光譜數(shù)據(jù)對土壤鹽漬化評價進(jìn)行了研究[3]。Hyperion數(shù)據(jù)是目前少有的星載高光譜數(shù)據(jù)[4],當(dāng)前國內(nèi)外已經(jīng)積極利用該數(shù)據(jù)進(jìn)行地表參數(shù)的定量反演,并且在諸多領(lǐng)域的研究中方法較為成熟,如對植被指數(shù)、植被葉面積指數(shù)、干物質(zhì)積累量、水體污染等方面的研究[5-6],但對土壤鹽分含量的研究還不多見。土壤鹽分含量反演比較復(fù)雜,不是一個簡單的線性問題,BP算法是目前應(yīng)用最為廣泛的神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)算法,可以實(shí)現(xiàn)輸入和輸出間的任意非線性映射和泛化功能[7]。

1 材料與方法

1.1 技術(shù)路線

本研究采用 Hyperion數(shù)據(jù),在對影像進(jìn)行預(yù)處理、特征選擇的基礎(chǔ)上,運(yùn)用BP神經(jīng)網(wǎng)絡(luò)模型反演土壤的含鹽量,探討建立高光譜數(shù)據(jù) Hyperion土壤含鹽量反演模型的過程和方法。技術(shù)路線如圖1所示。

1.2 樣品采集與處理

松嫩平原屬于干旱半干旱氣候區(qū),區(qū)內(nèi)現(xiàn)有鹽堿化土地373萬hm2,是世界上三大蘇打鹽堿地集中分布區(qū)域之一[8]。試驗(yàn)區(qū)設(shè)在肇州與肇源市的西南處,考慮到鹽堿土空間分布的不均勻性,在野外采樣前先對該地區(qū)的影像進(jìn)行了解譯,結(jié)合鹽堿土統(tǒng)計資料確定野外采樣路線。為盡量減小植被等因素對研究結(jié)果的影響,樣品采集時間為4月份,取得厚度為0～5 cm表層土樣200個,并測得采樣點(diǎn)的平面坐標(biāo)與高程數(shù)據(jù)。樣地點(diǎn)的土壤含鹽總量測定采用電導(dǎo)法,測定結(jié)果統(tǒng)計見表1。

圖1 技術(shù)路線

表1 土壤含鹽量樣地統(tǒng)計

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

對 Hyperion數(shù)據(jù)采用目視法,從242個波段中進(jìn)行波段選擇,保留176個波段。從輻射校正、幾何精校正、去噪聲、去條帶等方面對數(shù)據(jù)進(jìn)行預(yù)處理。假定地表為近似朗伯反射體,并且地物反射率與DN值為線性關(guān)系,采用回歸分析法,先對 Hyperion進(jìn)行大氣校正,消除大氣影響,然后使用實(shí)測光譜值進(jìn)行校正;由于研究區(qū)域地勢平坦,采用二次多項(xiàng)式函數(shù)進(jìn)行幾何精校正,校正后誤差控制在0.5個像元以內(nèi)。

1.4 影像與土壤含鹽量的相關(guān)性分析

采用相關(guān)系數(shù)進(jìn)行單波段光譜反射率與土壤鹽分分析,公式為

式中:r為光譜反射率與土壤鹽分的單相關(guān)系數(shù),i為波段序號,Rni為第n個土壤樣本第i波段的光譜反射率,Ri為N個土壤樣本在i波段的光譜反射率平均值,yn為實(shí)測的第n個土樣的鹽分,y 為實(shí)測的N個土樣鹽分的平均值,N是土壤樣本的數(shù)目。

根據(jù)土壤含鹽量與單波段相關(guān)性分析(見圖2)可以得出:土壤含鹽量在可見光和近紅外波段相關(guān)系數(shù)比較大,其中31、34波段間相關(guān)系數(shù)最高,達(dá)到0.695;相關(guān)性比較好的波段集中在500～900 nm之間;對該數(shù)據(jù)分別提取了主成分分析(K-L)、獨(dú)立成分分析、最大噪聲主分量變換(MNF)后,將特征圖像與土壤含水量的實(shí)測值進(jìn)行了相關(guān)性分析,結(jié)果表明:經(jīng)過降維映射后,相關(guān)性并沒有明顯的提高。以上試驗(yàn)結(jié)論為選擇最佳變量反演土壤含鹽量提供依據(jù)。

圖2 單波段與土壤含鹽量的相關(guān)系數(shù)散點(diǎn)圖

1.5 BP網(wǎng)絡(luò)模型設(shè)計

本次研究的樣本總數(shù)為200個,為確保精度,構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)模型的層數(shù)為4層,雙隱含層的節(jié)點(diǎn)數(shù)根據(jù)前人的經(jīng)驗(yàn)公式確定[9]式中:n為隱含層節(jié)點(diǎn)數(shù),ni為輸入節(jié)點(diǎn)數(shù),n0為輸出節(jié)點(diǎn)數(shù),a為1～10之間的常數(shù)。

輸入層節(jié)點(diǎn)數(shù)為176,為了克服少數(shù)異常值帶來的干擾有必要在建模和統(tǒng)計分析之間檢驗(yàn),剔除樣本中的異常值。構(gòu)建的網(wǎng)絡(luò)模型如圖3所示。

圖3 4層神經(jīng)網(wǎng)絡(luò)模型

選擇3種傳遞函數(shù):Log-sigmoid、Tan-sigmoid和線性函數(shù),一般而言,輸入層和隱層多采用非線性傳遞函數(shù),輸出層采用線性函數(shù),以保證輸出的范圍。

BP網(wǎng)絡(luò)的訓(xùn)練算法主要有traingd、traingdm、traingdx、trainrp、traincgf、traincgp 、traincgb、trainscg、trainbfg、trainoss、trainlm 等[9]。測試結(jié)果以計算檢驗(yàn)樣本實(shí)測值和模擬值的誤差均方根RRMSE和相關(guān)系數(shù)表示。誤差均方根 RRMSE公式為

其中:Rmeas為檢驗(yàn)樣本的實(shí)測值,Rmod為模型的預(yù)測值。

通過上述設(shè)計,采用MA TLAB6.5建立神經(jīng)網(wǎng)絡(luò)模型,神經(jīng)網(wǎng)絡(luò)模型構(gòu)成后,對BP網(wǎng)絡(luò)訓(xùn)練并對網(wǎng)絡(luò)進(jìn)行仿真,得到土壤含鹽量的反演圖。圖4(a)表示了反演的土壤含鹽量的部分圖斑,圖4(b)表示了以隱含層和輸出層的傳遞函數(shù)為purelin,訓(xùn)練函數(shù)為trainbr時,部分圖斑的反演圖。研究區(qū)域的土壤含鹽量反演前后比較見圖5,其中圖5(a)為研究區(qū)域,圖5(b)為土土壤含鹽量的反演圖。

1.6 BP神經(jīng)網(wǎng)絡(luò)模型性能分析

1.6.1 選擇不同傳遞函數(shù)和訓(xùn)練函數(shù)

取目標(biāo)誤差為0.001。當(dāng)隱含層的傳遞函數(shù)為Log-sigmoid、輸出層為 purelin函數(shù)時,以 trainlm函數(shù)的收斂速度最快,隨著訓(xùn)練次數(shù)的增大,誤差均方根逐漸減小,次數(shù)達(dá)到8時,誤差均方根約為0,該模型的相關(guān)系數(shù)為0.724,數(shù)據(jù)與模型的擬合程度最好。當(dāng)隱含層傳遞函數(shù)為 Tan-sigmoid、輸出層的傳遞函數(shù)為purelin時,訓(xùn)練函數(shù)為trainlm收斂的速度最快,達(dá)到了要求的目標(biāo)誤差,誤差均方根在訓(xùn)練7次后幾乎為0,模型的相關(guān)系數(shù)達(dá)到了0.733,擬合程度最好。當(dāng)隱含層的傳遞函數(shù)為purelin、輸出層的傳遞函數(shù)為purelin時,以trainbr的擬合程度最好。

1.6.2 與其它反演模型的精度分析

以隱含層為log-singmid,輸出層為purelin函數(shù)時,訓(xùn)練函數(shù)為trainlm時,隱含層節(jié)點(diǎn)數(shù)為20時的模型為代表,計算該模型實(shí)際精度,公式為[5]

實(shí)際精度達(dá)到80.69%。為更進(jìn)一步比較,將高光譜數(shù)據(jù)單波段與土壤含鹽量進(jìn)行相關(guān)性分析,得到相關(guān)系數(shù)最大的波段,并建立BP神經(jīng)網(wǎng)絡(luò)模型,并分別建立一次(L IN)、二次(QUA)、三次(CUB)曲線經(jīng)驗(yàn)統(tǒng)計模型,分別計算他們的實(shí)際精度,表2列出了不同的模型反演土壤含鹽量的實(shí)際精度和模型的相關(guān)系數(shù),從表2中可以看出,BP神經(jīng)網(wǎng)絡(luò)模型的實(shí)際精度高于傳統(tǒng)的經(jīng)驗(yàn)統(tǒng)計模型的實(shí)際精度。

表2 不同反演模型的精度分析

2 結(jié) 論

本文運(yùn)用 Hyperion數(shù)據(jù)對土壤含鹽量定量研究,將BP神經(jīng)網(wǎng)絡(luò)模型運(yùn)用到高光譜數(shù)據(jù)對研究地區(qū)土壤含鹽量的反演中。通過研究得出:

1)采用高光譜的 Hyperion數(shù)據(jù),采用BP神經(jīng)網(wǎng)絡(luò)進(jìn)行土壤含鹽量提取是可行的;另外,J.Farifteh[10]等人曾對PLSR和人工神經(jīng)網(wǎng)絡(luò)方法進(jìn)行了對比分析,所以本文建議也可以利用 PLSR(偏最小二乘回歸法)土壤含鹽量預(yù)測模型,進(jìn)行含鹽量的預(yù)測。

2)神經(jīng)網(wǎng)絡(luò)模型具有極強(qiáng)的線性和非線性擬合能力,模擬遙感影像特征與土壤鹽分之間比較復(fù)雜的關(guān)系上有很大優(yōu)勢,但很難全面解釋神經(jīng)網(wǎng)絡(luò)作出決策或者產(chǎn)生輸出的過程,因?yàn)樯窠?jīng)網(wǎng)絡(luò)缺乏洞察數(shù)據(jù)集特性的解釋能力。有關(guān)神經(jīng)網(wǎng)絡(luò)在土壤含鹽量的反演的研究工作還有待于進(jìn)一步深入。

[1]安永清.內(nèi)蒙古河套灌區(qū)土壤鹽堿化遙感監(jiān)測方法研究[J].遙感技術(shù)與應(yīng)用,2008,23(3):316-319.

[2]熊文成,邵蕓.氯化鈉鹽土壤介電虛部特性的初步研究[J].遙感學(xué)報,2006,10(2):280-282.

[3]TAYLOR G,DEHAAN R.Salinity mapping with hyperspectral imagery[EB/OL].http//www.geology.unsw.edu.au/research/Remote Sensing,2000.

[4]王芳.基于 Hyperion高光譜數(shù)據(jù)的城市植被脅迫評價[J].應(yīng)用生態(tài)學(xué)報,2007,18(6):1286-1292.

[5]王秀珍,黃敬峰,李云梅,等.水稻葉面積指數(shù)的高光譜遙感估算模型[J].遙感學(xué)報,2004,8(1):81-88.

[6]HANSEN P M,SCHJOERRING J K.Reflectance measurement of canopy biomass and nitrogen status in w heat crops using no rmalized difference vegetation indices and partial least squares regression〔J〕.Remote Sensing of Environment.2003,86:542-553.

[7]飛思科技產(chǎn)品研發(fā)中心.matlab6.5輔助神經(jīng)網(wǎng)絡(luò)分析與設(shè)計[M].北京:電子工業(yè)出版社,2003.

[8]張巍,馮玉杰.松嫩平原不同鹽漬土條件下藍(lán)藻群落的生態(tài)分布[J].生態(tài)學(xué)雜志,2008,27(5):718-722.

[9]葛哲學(xué),孫志強(qiáng).神經(jīng)網(wǎng)絡(luò)理論與MA TLABR2007實(shí)現(xiàn)[M].北京:電子工業(yè)出版社,2007:54-56.

[10]J.FARIFTEH,F.VAN DER M EER,C.A TZBERGER,et al..Quantitative analysisof salt-affected soil reflectance[J].Remote Sensing of Environment,2007(110):59-78.

Study on inversion of soil salinity

TANG Yan
(Northeast Forestry University,Harbin 150040,China)

One experimental area in Daqing city in Heilongjiang p rovince is taken as an examp le to perfo rm the quantitative inversion of soil salinity using Hyperion data in this paper.The inversion method of soil salinity using Hyperion data is discussed by the image p rep rocessing,the feature extraction and the establishment of BP neural netwo rk model.It gives a lot of help in soil suveying system and p romoting the development in quantitative retrievalof soil salinity.M eanw hile,thismodel p rovides reference for solving other non-linear p roblem s.

soil salinity;Hyperion data;inversion;BP neural network model

TP751

A

1006-7949(2010)06-0065-03

2010-06-28

唐 彥(1957-),女,副教授.

[責(zé)任編輯張德福]