張靜

摘要 [目的]對土壤有機碳含量進行預測研究。[方法]利用高光譜儀對表層土壤進行光譜測定并且進行光譜數(shù)據(jù)的預處理，通過多元線性逐步回歸（SMLR）和偏最小二乘回歸（PLSR）方法對土壤有機碳含量進行預測，并對2種模型的精度進行比較。[結(jié)果]LSR模型的精度高于SMLR模型。[結(jié)論]偏最小二乘回歸法優(yōu)于多元逐步回歸法，對有機碳的預測具有更好的效果。

關鍵詞 土壤；有機碳；高光譜；多元線性逐步回歸；偏最小二乘回歸

中圖分類號 S126 文獻標識碼 A 文章編號 0517-6611（2018）02-0001-03

Abstract [Objective]To predict soil organic carbon content.[Method]Surface soil was detected by high spectrometer spectrometric and spectral data was treated， through stepwise multiple linear regression（SMLR） and partial leastsquares regression（PLSR） method，soil organic carbon content was predicted ，and the accuracy of the two models was compared.[Result]The accuracy of PLSR model was higher than SMLR model.[Conclusion]PLSR method is better than SMLR method in forecasting organic carbon.

Key words Soil；Organic carbon；Hyperspectral；SMLR；PLSR

土壤是人類生活和農(nóng)業(yè)生產(chǎn)的支撐和基礎，對土壤類型的區(qū)分可以幫助人類因地制宜地利用土壤，同時也給土地利用規(guī)劃、土地評價、農(nóng)業(yè)技術(shù)推廣提供重要的依據(jù)。

伴隨著高光譜遙感技術(shù)的不斷發(fā)展，可利用高光譜技術(shù)提取具有診斷意義的地物光譜特征。土壤光譜反射率是反映土壤理化性質(zhì)的綜合指標，故土壤的光譜特征可用于土壤各種物質(zhì)含量的反演、分類等。Bowers等[1]發(fā)現(xiàn)影響土壤反射特性的主要原因是土壤有機質(zhì)，并且還氧化了土壤有機質(zhì)，發(fā)現(xiàn)能提高土壤的反射系數(shù)。徐彬彬等[2]發(fā)現(xiàn)在可見—近紅外波段，有機質(zhì)對土壤光譜的影響最大。Van Waes等[3]基于一階微分變換結(jié)合偏最小二乘回歸預測模型，成功預測了自然草地和農(nóng)業(yè)草地的土壤有機質(zhì)含量。盧艷麗等[4]利用逐步回歸分析和主成分回歸分析等統(tǒng)計方法對土壤進行了處理，最終建立了東北黑土有機質(zhì)含量預測模型。

目前國內(nèi)外對高光譜土壤有機碳含量研究進行了大量的工作，但是由于不同土壤類型、不同區(qū)域土壤都有其各自的特性，研究特定區(qū)域的土壤就勢在必行。筆者主要是對豐樂河流域的土壤進行研究，基于土壤高光譜反射特性對土壤進行反射率、一階微分和連續(xù)統(tǒng)去除法的數(shù)據(jù)變換，并且通過多元線性逐步回歸模型和偏最小二乘回歸模型對土壤有機碳含量進行預測。

1 材料與方法

1.1 土壤樣品的采集

研究區(qū)為豐樂河流域，為巢湖的重要支流之一。地理坐標為116°45′～117°21′E、31°18′～31°33′N。豐樂河流域西起橫塘崗鄉(xiāng)，東至柿樹崗鄉(xiāng)，南起東河口鎮(zhèn)，北至椿樹鎮(zhèn)，面積為1 500 km2。流域地形以丘陵為主，土地利用類型主要為林地和農(nóng)業(yè)耕作用地，林地主要為松林和果園，農(nóng)業(yè)耕作用地一般為水田和旱地，其土壤類型也較多，有黃褐土、黃棕壤、紅壤、潴育水稻土等。土壤樣品采集自表層0～20 cm，共162個土樣，每個土樣分成2份：1份過20目篩，用于光譜測量；1份過100目篩，用于化學方法測定土樣有機碳含量。光譜測量通過室內(nèi)光譜儀進行，化學分析利用重鉻酸鉀-硫酸亞鐵滴定法進行。

1.2 光譜的測定

光譜測試采用荷蘭Avantes公司生產(chǎn)的型號分別為Avaspec-2048×14、Avaspec-NIR256-2.5的2種地物光譜儀，其中Avaspec-2048×14光譜儀光譜范圍為188～1 170 nm，光譜分辨率為2.4 nm，光譜采樣間隔0.6 nm。型號Avaspec-NIR256-2.5光譜范圍是928～2 528 nm，光譜分辨率為15 nm，光譜采樣間隔6 nm。2種光譜儀所輸出的光譜曲線均設定由10條原始掃描光譜自動平均所獲。對于每一個土壤樣品進行測試前使用標準參考白板進行定標，同時每份土壤從4個不同方向?qū)λ型翗舆M行平均值運算，最終得到土壤樣品的光譜反射率數(shù)據(jù)。

1.3 光譜數(shù)據(jù)預處理

1.3.1 異常樣本剔除。在獲取樣本土樣和測量樣本時，由于各因素的影響，會出現(xiàn)部分誤差甚至錯誤，因此在進行數(shù)據(jù)處理前，需對樣本進行檢驗，判定和剔除異常樣本。該研究采用箱形圖對異常樣本進行剔除，剔除后樣本共156個。

1.3.2 光譜曲線平滑。

光譜曲線平滑是為了剔除獲取的原始光譜曲線的高頻噪聲，以及曲線上的細小“毛刺”現(xiàn)象，得到一條較為平滑的曲線。該研究使用移動平均法，對獲取的原始光譜反射率數(shù)據(jù)進行平滑去噪處理。移動平均法是對包括待平滑光譜數(shù)據(jù)點及前后各若干個數(shù)據(jù)點進行平均，用平均值取代該點原來的光譜值[5]。圖1為平滑后的曲線圖。

1.3.3 微分變換。

一階微分變換下的光譜曲線值反映了原始光譜反射率曲線的斜率值及其變化情況，獲得相應的曲線拐點。光譜一階微分的最大值所對應的波段是光譜反射率增加速度最大的波段位置，而最小值則是反映了光譜反射率減少速度最大的波段位置。對于一個離散光譜，其波長k點處的一階微分[6]為：

對光譜曲線進行分變換，主要是為了擴大重疊混合部分光譜的差異，以便更好地辨別樣品之間的光譜特征[8]。另外光譜微分變換可以降低因土壤組成成分含量不同而引起的對土壤光譜特征的影響[7]，并且還可以減弱由大氣散射和吸收作用帶來的對地物光譜特征的影響。如圖2所示，在1 000～2 300 nm波段之間，原始光譜中的細節(jié)吸收谷和吸收峰也得到了凸顯。

1.3.4 連續(xù)統(tǒng)去除。連續(xù)統(tǒng)去除法即去包絡線法。光譜曲線由一條接連不斷的折線段包裹（這個折線段即稱為包絡線），將包絡線看作背景，去掉包絡線后可以壓制光譜背景，獲得光譜特征吸收帶，提取光譜特征信息。光譜特征吸收帶的計算公式是：

式中，Rc、Rcr、R分別是光譜特征吸收帶、包絡線和原始光譜反射率；λ為波長。

連續(xù)統(tǒng)去除法可顯著地突出光譜曲線的反射和吸收特征，并將其歸一到同一個光譜背景上，在與其他光譜曲線的特征數(shù)值的比較上具有很大的優(yōu)勢[9]。如圖3所示，原始數(shù)據(jù)經(jīng)過去包絡線處理后，3個水汽吸收帶在1 400、1 900、2 200 nm處的吸收谷更加明顯，還可以看出在可見光 500 nm附近有2個又窄又深的吸收谷，800～1 100 nm 波段處還有一個又寬又淺的吸收谷。

1.4 建模方法與精度評價

該研究采用多元線性逐步回歸模型（SMLR）和偏最小二乘回歸模型（PLSR）分別建立土壤有機碳含量預測模型。SMLR是按規(guī)則對輸入的自變量進行篩選，找出最優(yōu)組合建立回歸方程，它要求自變量與因變量之間為線性關系，并且自變量之間不存在嚴重的共線性，因變量符合正態(tài)分布；選擇特征波段數(shù)據(jù)作為預測模型的因變量。PLSR 在進行預算前先對數(shù)據(jù)進行降維，在數(shù)據(jù)壓縮時，兼顧了目標變量矩陣的作用[10]。另外PLSR對樣本數(shù)量的要求較低，樣本個數(shù)小于自變量個數(shù)時，也可建立回歸模型。PLSR突出的特點是逐步提取光譜數(shù)據(jù)中的成分，并且還考慮了對目標變量的最大解釋。

模型精度評價主要按照以下幾個參數(shù)進行：均方根誤差RMSE，決定系數(shù)R2，以及統(tǒng)計量F所對應的顯著性水平Sig.。RMSE越小越好，R2越大越好，Sig.<0.05即代表顯著性水平高，模型有效。

2 結(jié)果與分析

2.1 有機碳的多元線性逐步回歸模型

該研究將156個土壤樣本分為建模集97個，驗證集59個。對土樣進行相關分析選取特征波段，分別將反射率、一階微分以及連續(xù)統(tǒng)去除這3個光譜指標的特征波段作為多元逐步回歸模型的自變量，土壤有機碳含量作為因變量，建立有機碳含量反演模型。表1為多元逐步回歸模型的檢驗結(jié)果。觀察發(fā)現(xiàn)3種模型精度都不高，可能的原因是土壤中有機碳成分復雜，其光譜特性很難用幾個波段解釋，只關注相關性高的波段，而忽略了整個波段都有一定的相關性。相比較之下，一階微分變換下R2最大，RMSE最小，模型效果最好，可能是因為其擴大了重疊部分的光譜差異。

多元逐步回歸模型的有機碳含量反演模型表達如下。

圖4是多元線性逐步回歸模型下土壤有機碳含量實測值與預測值的散點圖。樣點越接近1∶1線，說明模型預測的精度越好。從圖4可以看出，3種數(shù)據(jù)下的散點圖緊湊性差不多，一階微分變換下的R2最大。綜合考慮模型的精度穩(wěn)

定性和預測的實際精度，一階微分變換下的多元逐步回歸所得的模型效果最好。

2.2 有機碳的偏最小二乘回歸模型

與構(gòu)建SMLR的建模集和驗證集相同，但PLSR是全波段參與建模，數(shù)據(jù)量龐大，它可對龐大且有多重共線性的土壤高光譜數(shù)據(jù)進行壓縮，并提取出與因變量相關的主成分。

從表2可以看出，3種偏最小二乘回歸模型的建模精度都很高，反射率、一階微分、連續(xù)統(tǒng)去除法的R2相差不多，驗證模型也相差不多。圖6為3種偏最小二乘回歸模型下的預測值與實測值的比較，可以看出，樣點幾乎都分布于1∶1線的周圍，緊湊度各有差異，相對而言，一階微分變換下的散點圖比反射率、連續(xù)統(tǒng)去除法的散點圖稍好，R2值也稍大，故綜合考慮，一階微分變換下的偏最小二乘回歸模型效果最好。

表2與表1相比，明顯表2的建模精度高于表1。圖5與圖4相比，圖5的3幅圖整體精度也高于圖4。由此對比可發(fā)現(xiàn)，PLSR優(yōu)于SMLR，預測值更逼近于實測值。

3 結(jié)語

該研究對土壤樣品進行異常樣本的剔除，避免影響有機碳含量預測精度。運用反射率、一階微分以及連續(xù)統(tǒng)去除3種光譜數(shù)據(jù)進行多元逐步回歸（SMLR）和偏最小二乘回歸分析（PLSR），發(fā)現(xiàn)一階微分變換下的SMLR和PLSR效果最好。當研究的土壤樣本中某一物質(zhì)組分含量較低且無明顯的吸收峰時，進行微分變換后，易找到相關性高的波段，使得該方法在反演土壤成分含量時應用較多；當土壤光譜信息重疊時，微分變換還可以擴大光譜之間的特征差異，故在土壤屬性的反演和類型區(qū)分上較多使用。對于所建模型結(jié)果，PLSR精度總體上比SMLR要高一些，尤其是建模精度。PLSR解決了自變量之間的多重共線性問題，能提供更可靠的建模結(jié)果，對于樣本個數(shù)小于自變量個數(shù)的情況下比較有利。驗證集的預測值和實測值散點圖也表明了整體性PLSR優(yōu)于SMLR，因此PLSR具有快速預測有機碳含量的潛力。

參考文獻

[1]BOWERS S，HANKS R J.Reflection of radiant energy from soils[J].Soil science， 1965， 100（2）：130-138.

[2] 徐彬彬， 季耿善，朱永豪.中國陸地背景和土壤光譜反射特性的地理分區(qū)的初步研究[J].環(huán)境遙感，1991，6（2）：142-151.

[3] VAN WAES C，MEATDAGH I，LOOTENS P，et al.Possibilities of near infrared reflectance spectroscopy for the prediction of organic carbon concentrations in grassland soil[J].Journal of agricultural science，2005，143（6）：487-492.

[4] 盧艷麗，白由路，楊俐蘋，等.基于主成分回歸分析的土壤有機質(zhì)高光譜預測與模型驗證[J].植物營養(yǎng)與肥料學報，2008，14（6）：1076-1082.

[5] 張芳，熊黑鋼，龍?zhí)遥?實測反射率與影像反射率對土壤堿化預測的對比分析[J].光譜學與光譜分析，2011，31（1）：227-232.

[6] 王欣.近紅外分析中光譜預處理方法的研究與應用進展[J].科技資訊，2013（15）：2.

[7] CLOUTIS E A.Hypersoectral geological remote sensing： Evaluation of analytical techniques[J].International journal of remote sensing，1996，17（12）：2215-2242.

[8] 于士凱，姚艷敏，王德營，等.基于高光譜的土壤有機質(zhì)含量反演研究[J].中國農(nóng)學通報，2013，29（23）：146-152.

[9] 謝伯承，薛緒掌，劉偉東，等.基于包絡線法對土壤光譜特征的提取及其分析[J].土壤學報，2005，42（1）：171-175.

[10] 紀文君.基于野外viS-NIR高光譜的土壤屬性預測及田間水分影響去除研究[D].杭州：浙江大學，2014：67-72.