徐 璐， 陳奕云, 2, 3*， 洪永勝， 魏 鈺， 郭 龍， Marc Linderman

1. 武漢大學(xué)資源與環(huán)境科學(xué)學(xué)院， 湖北 武漢 430079 2. 自然資源部數(shù)字制圖與國(guó)土信息應(yīng)用重點(diǎn)實(shí)驗(yàn)室， 湖北 武漢 430079 3. 土壤與農(nóng)業(yè)可持續(xù)發(fā)展國(guó)家重點(diǎn)實(shí)驗(yàn)室， 江蘇 南京 210008 4. 華中農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院， 湖北 武漢 430070 5. Geographical and Sustainability Sciences, The University of Iowa, Iowa City, IA 52246, USA

引 言

土壤有機(jī)碳(soil organic carbon, SOC)是土壤的重要組成， 在碳循環(huán)與全球氣候變化、 土壤健康與糧食安全等方面起著至關(guān)重要的作用[1-2]。 SOC信息的高效獲取對(duì)于可持續(xù)土壤資源管理與政策制定具有重要意義[3-4]。 因此， 探索快速且精確的SOC含量估測(cè)方法成為了一個(gè)熱點(diǎn)問(wèn)題。 隨著遙感技術(shù)的發(fā)展， 可見近紅外光譜技術(shù)逐漸被用于SOC含量估測(cè)研究[5]。 與傳統(tǒng)的化學(xué)方法相比， 此方法具有快速高效、 低成本、 無(wú)污染的特點(diǎn)。

已有研究表明， 利用非成像的可見近紅外光譜技術(shù)能夠較好地估測(cè)SOC含量[6]。 但是這些研究大多需要對(duì)野外采集到的土壤進(jìn)行風(fēng)干、 研磨、 過(guò)篩等預(yù)處理， 以降低土壤粗糙度對(duì)土壤屬性光譜建模的不利影響。 自然狀況下的土壤顆粒因團(tuán)聚作用形成不同直徑的土壤團(tuán)聚體， 呈現(xiàn)非定向空間分布， 具有一定程度的隨機(jī)粗糙度[7]。 這會(huì)導(dǎo)致入射光在土壤樣本的不同位置發(fā)生不同程度的散射、 吸收和反射； 同時(shí)， 土壤顆粒大小及顆粒間的陰影部分也會(huì)對(duì)光譜反射率造成影響。 顆粒粒徑越大， 顆粒間陰影部分面積越大， 光譜反射率越低。 這些特性導(dǎo)致高粗糙度的土壤樣本表面具有復(fù)雜的光譜散射與反射特征。

非成像光譜技術(shù)僅能夠獲取土壤某個(gè)點(diǎn)處的反射光譜信息， 無(wú)法同時(shí)獲取土壤光譜的空間連續(xù)分布信息， 且容易受到土壤粗糙度對(duì)光譜測(cè)量的影響， 因此難以成功估測(cè)高粗糙度樣本的SOC含量。 土壤樣本經(jīng)過(guò)風(fēng)干、 研磨、 過(guò)篩的預(yù)處理后， 其顆粒大小與形狀相對(duì)均勻， 從而降低了粗糙度， 使得同一個(gè)土壤樣本的表面各點(diǎn)具有近似一致的光譜反射特征。 這為非成像光譜數(shù)據(jù)實(shí)現(xiàn)SOC含量的精確估算提供了保障[8]。 但這一預(yù)處理過(guò)程需要耗費(fèi)大量時(shí)間、 人力、 物力， 具有較高的成本， 時(shí)效性差， 難以實(shí)現(xiàn)大尺度的野外土壤光譜測(cè)量及SOC含量的實(shí)時(shí)監(jiān)測(cè)。 因此， 開展能夠克服土壤粗糙度影響的SOC光譜反演具有重要的科研應(yīng)用價(jià)值。

成像光譜技術(shù)融合了光譜技術(shù)和成像技術(shù)， 能夠同時(shí)獲取土壤樣本中每個(gè)像素點(diǎn)的光譜信息及其空間分布信息。 已有少數(shù)研究將成像光譜技術(shù)成功應(yīng)用于低粗糙度樣本的SOC估算的報(bào)道[9]。 但是成像光譜技術(shù)估算高粗糙度樣本SOC含量的潛力仍有待驗(yàn)證。

因此， 以研磨前后不同粗糙度下的土壤樣本為研究對(duì)象， 提取高光譜影像中的土壤光譜信息， 基于5種光譜預(yù)處理策略， 利用偏最小二乘回歸(PLSR)和支持向量回歸(SVR)算法建立SOC估測(cè)模型， 并與利用非成像光譜建立的SOC估測(cè)模型性能進(jìn)行對(duì)比， 考察成像光譜技術(shù)估測(cè)高粗糙度下土壤樣本SOC含量的能力， 并分析其估測(cè)機(jī)理， 為土壤有機(jī)碳儲(chǔ)量實(shí)時(shí)監(jiān)測(cè)提供新技術(shù)。

1 實(shí)驗(yàn)部分

1.1 研究區(qū)概況及樣本采集

研究區(qū)位于美國(guó)愛(ài)荷華州東南部。 該地區(qū)的海拔位于237～266 m之間， 坡度為0°～9.7°， 年平均氣溫為9.88 ℃， 年平均降雨量為903 mm。 土壤類型主要為粉質(zhì)粘壤土、 粘壤土和粉砂壤土。 土地利用類型以耕地為主。 該研究區(qū)主要種植玉米、 大豆、 黑麥等農(nóng)作物。

圖1 研究區(qū)地理位置及采樣點(diǎn)分布情況

2015年10月， 在研究區(qū)利用網(wǎng)格采樣的策略共采集了60個(gè)土壤樣本(圖1)。 土壤樣本采樣間隔為130 m。 每個(gè)樣本是由五個(gè)子樣本組成的復(fù)合樣本(子樣本采集于1 m×1 m網(wǎng)格的四個(gè)頂點(diǎn)及中心點(diǎn))。 采樣前， 將土壤表面的植被、 雜草和其他雜物去除。 采樣時(shí)， 將土壤樣本按照采樣順序裝入帶有標(biāo)簽的塑料袋中并密封。 采樣后， 將采集到土壤樣本帶回實(shí)驗(yàn)室進(jìn)行光譜測(cè)量及化學(xué)分析。

1.2 光譜測(cè)量及化學(xué)分析

在實(shí)驗(yàn)室中將采集的土壤樣本進(jìn)行自然風(fēng)干。 首先測(cè)量風(fēng)干后未研磨的土壤樣本的成像光譜與非成像光譜。 之后將所有樣本進(jìn)行研磨和過(guò)篩(2 mm)預(yù)處理， 再分別進(jìn)行成像與非成像光譜測(cè)量。 最后， 將樣本繼續(xù)研磨、 過(guò)篩(0.15 mm)， 用于測(cè)量SOC含量。

采用Cubert UHD185型成像光譜儀測(cè)量土壤樣本研磨前后的成像光譜。 該成像光譜儀是由德國(guó)Cubert Gmbh公司研發(fā)， 具有實(shí)時(shí)成像的特點(diǎn)， 光譜范圍為450～998 nm， 光譜分辨率為4 nm[10]。 光譜測(cè)量在暗室中進(jìn)行， 以鹵素?zé)糇鳛槲ㄒ还庠矗?入射光與土壤樣本成45°夾角， 與樣本相距30 cm。 UHD185型成像光譜儀被固定于土壤樣本正上方12 cm處。 測(cè)量土壤樣本光譜前， 需要進(jìn)行白板測(cè)量以校正光譜儀。

采用ASD FieldSpec3便攜式光譜儀測(cè)量土壤樣本研磨前后的非成像光譜。 整個(gè)光譜測(cè)量環(huán)境與成像光譜測(cè)量環(huán)境相同， 在暗室中進(jìn)行。 鹵素?zé)糇鳛槲ㄒ还庠矗?與土壤樣本成45°夾角， 并與樣本相距30 cm。 光纖放置在樣本正上方， 距離樣本12 cm。 在測(cè)量光譜之前利用白板校正光譜儀。 每個(gè)樣本測(cè)量10次光譜， 取平均值作為最終的樣本光譜。 光譜范圍為350～2 500 nm， 光譜分辨率為1 nm。

采用Vario MACRO cube元素分析儀測(cè)量SOC含量[11]。 測(cè)量前， 需要對(duì)儀器進(jìn)行氣密性檢驗(yàn)。 測(cè)量時(shí)， 將用錫箔紙包好的樣本按順序放入元素分析儀中測(cè)量SOC含量。 元素分析儀中氧氣和氦氣的氣瓶減壓閥的輸出壓力分別調(diào)至0.20～0.22和0.20～0.23 MPa， 燃燒管溫度設(shè)置為1 150 ℃， 還原管溫度設(shè)置為850 ℃。

1.3 光譜預(yù)處理

測(cè)量得到的單個(gè)樣本成像光譜影像大小為50×50個(gè)像素點(diǎn)(圖2)。 為了降低周邊測(cè)量環(huán)境的影響(如盛放土壤樣本的培養(yǎng)皿等)， 提取紅色框中14×20共280個(gè)像素點(diǎn)的平均光譜作為土壤樣本的成像光譜數(shù)據(jù)。 由于儀器固有特性以及實(shí)際測(cè)量誤差， 導(dǎo)致450～489 nm范圍內(nèi)信噪比較低， 故選擇490～998 nm范圍內(nèi)的成像光譜數(shù)據(jù)用于建模分析。 成像光譜與非成像光譜數(shù)據(jù)波段范圍與光譜分辨率不一致。 為了有效對(duì)比兩者數(shù)據(jù)估算SOC的潛力， 本文選擇490～998 nm范圍內(nèi)非成像光譜數(shù)據(jù)進(jìn)行建模， 并且將非成像光譜反射率重采樣為4 nm。 此外， 通過(guò)原始成像與非成像光譜進(jìn)行去包絡(luò)線(CR)、 吸光度變換(AB)、 S-G平滑(SG)、 標(biāo)準(zhǔn)正態(tài)變換(SNV)、 多元散射校正(MSC)5種光譜預(yù)處理后， 對(duì)原始光譜及預(yù)處理后光譜開展了SOC估測(cè)建模。

圖2 單個(gè)土壤樣本研磨前及研磨后成像光譜影像

1.4 模型構(gòu)建及驗(yàn)證

本研究中采用偏最小二乘回歸(PLSR)算法和支持向量回歸(SVR)算法建立不同預(yù)處理后成像與非成像光譜的SOC含量估測(cè)模型。 PLSR是一種典型的線性回歸算法， 由Wold提出， 被廣泛應(yīng)用于光譜分析研究。 該方法將原始可見近紅外光譜波段投影到一個(gè)新的空間， 形成新的特征變量， 利用新的特征變量建立回歸模型。 它集合了多元線性回歸、 主成分分析和典型相關(guān)分析優(yōu)點(diǎn)， 能夠有效降低高光譜波段的多重共線性和信息冗余。 SVR是利用核函數(shù)將原始可見近紅外光譜投影至高維線性空間尋找最優(yōu)估算結(jié)果的算法， 對(duì)高維數(shù)據(jù)具有良好的魯棒性。 常見的核函數(shù)包括線性核函數(shù)、 高斯核函數(shù)、 多項(xiàng)式核函數(shù)和徑向基核函數(shù)。 其中基于徑向基核函數(shù)建立的SVR模型為非線性模型， 效果最佳。 因此， 選擇徑向基核函數(shù)的SVR算法建立SOC估算模型。

根據(jù)濃度梯度法， 土壤樣本按2∶1比例劃分為建模集和驗(yàn)證集。 建模集被用來(lái)建立PLSR和SVR估測(cè)模型， 而驗(yàn)證集被用來(lái)評(píng)價(jià)模型性能。 評(píng)價(jià)指標(biāo)包括決定系數(shù)(R2)和均方根誤差(RMSE)。R2越大， RMSE越小， 模型性能越好。

2 結(jié)果與討論

2.1 SOC數(shù)據(jù)分布特征

表1描述了研究區(qū)土壤樣本全集、 建模集和驗(yàn)證集的SOC含量的統(tǒng)計(jì)特征。 樣本全集有60個(gè)樣本， 建模集有40個(gè)樣本， 驗(yàn)證集有20個(gè)樣本。 樣本全集與建模集中SOC含量范圍介于1.21%～2.97%之間， 驗(yàn)證集中SOC含量最大值為2.83%， 最小值為1.24%。 樣本全集、 建模集和驗(yàn)證集的均值皆為2.22%， 中位數(shù)為2.33%。 樣本全集、 建模集、 驗(yàn)證集的標(biāo)準(zhǔn)差分別為0.43， 0.44和0.42； 變異系數(shù)分別為0.19， 0.20和0.19， 屬于中等變異。 三個(gè)樣本集的偏度系數(shù)分別為-0.72， -0.71和-0.82， 峰度系數(shù)分別為-0.28， -0.29和0.04， 為峰值較小的非正態(tài)分布。 根據(jù)SOC的統(tǒng)計(jì)特征可知， 建模集和驗(yàn)證集的SOC 含量特征與樣本全集的SOC含量統(tǒng)計(jì)特征相似。 這表明建模集和驗(yàn)證集可以代表樣本全集， 能夠被用來(lái)建立和驗(yàn)證SOC估測(cè)模型。

表1 SOC含量統(tǒng)計(jì)性特征描述表

2.2 光譜特征

圖3展示了在490～998 nm范圍內(nèi)， 研磨前后土壤成像光譜與非成像光譜平均反射率特征。 其中， 成像光譜反射率是所有樣本中所有像素點(diǎn)的平均光譜反射率。 圖3中d,b曲線分別代表了土壤樣本研磨前后成像光譜的反射率。 土壤樣本研磨前的成像光譜反射率在490～974 nm范圍內(nèi)不斷增加， 在974～998 nm之間開始下降。 土壤樣本研磨后的成像光譜反射率的變化趨勢(shì)與研磨前的光譜反射率變化趨勢(shì)相似。 土壤樣本研磨后的光譜反射率在490～978 nm范圍內(nèi)不斷增加， 在978～998 nm范圍內(nèi)不斷降低。 在490～998 nm范圍內(nèi)， 土壤樣本研磨前的光譜反射率總是低于研磨后的光譜反射率， 且隨著光譜波長(zhǎng)不斷增加， 反射率變化越來(lái)越大。 圖3中c,a曲線分別代表了土壤樣本研磨前后非成像光譜反射率。 從圖中可以看出， 非成像光譜反射率在490～998 nm范圍內(nèi)不斷增加， 無(wú)下降趨勢(shì)。

圖3 土壤樣本研磨前后成像光譜(d, b)與非成像光譜(c, a)平均反射率

土壤樣本研磨前的成像光譜與非成像光譜反射率總是低于研磨后的光譜反射率。 隨著光譜波長(zhǎng)增加， 土壤樣本研磨前后光譜反射率差異逐漸變大， 土壤研磨前后非成像光譜反射率差異小于成像光譜反射率的差異。

2.3 SOC與光譜相關(guān)性分析

圖4展示了樣本研磨前后不同粗糙度條件下土壤成像光譜、 非成像光譜與SOC的皮爾遜相關(guān)系數(shù)曲線。 圖4中d,b曲線分別代表了土壤樣本研磨前、 后成像光譜與SOC的相關(guān)系數(shù)。 研磨前成像光譜反射率與SOC的相關(guān)系數(shù)絕對(duì)值介于0.49～0.70之間， 在606 nm處相關(guān)系數(shù)絕對(duì)值最大。 研磨后成像光譜反射率與SOC的相關(guān)系數(shù)絕對(duì)值范圍為0.75～0.85， 在666 nm處相關(guān)系數(shù)絕對(duì)值最大。

圖4 土壤樣本研磨前后成像光譜(d, b)、 非成像光譜(c, a)與SOC相關(guān)系數(shù)圖

對(duì)于非成像光譜， 研磨前光譜與SOC的相關(guān)系數(shù)最大絕對(duì)值為0.56， 最小絕對(duì)值為0.38。 598 nm處光譜與SOC的相關(guān)性最強(qiáng)。 研磨后光譜與SOC的相關(guān)系數(shù)最大值為0.84， 最小值為0.75。 668 nm處光譜與SOC的相關(guān)性最強(qiáng)。

總體來(lái)說(shuō)， 土壤樣本研磨前后成像光譜與非成像光譜與SOC的相關(guān)性皆為負(fù)相關(guān)。 在490～998 nm范圍內(nèi)， 土壤樣本研磨前， 成像光譜與SOC的相關(guān)性明顯強(qiáng)于非成像光譜與SOC的相關(guān)性； 土壤樣本研磨后， 成像光譜與SOC的相關(guān)性稍微強(qiáng)于非成像光譜與SOC的相關(guān)性。

2.4 模型性能對(duì)比

表2展示了兩種粗糙度下的土壤樣本(研磨前后)， 基于成像光譜與非成像光譜， 利用其原始光譜及5種預(yù)處理光譜， 建立的PLSR和SVR模型性能對(duì)比表。

表2 土壤樣本研磨前后成像光譜與非成像光譜在不同預(yù)處理下PLSR和SVR模型性能表

續(xù)表2

從表中可以看出， 土壤樣本研磨前后， 不同光譜預(yù)處理下， 基于成像光譜建立的PLSR和SVR模型性能總是優(yōu)于基于非成像光譜建立的估測(cè)模型。

總的來(lái)說(shuō)， 對(duì)于成像光譜， 土壤樣本研磨前后建立的估測(cè)模型皆能被用于估測(cè)SOC含量。 土壤樣本研磨前建立 PLSR模型性能總是優(yōu)于樣本研磨之后建立的PLSR模型， 而SVR模型性能正好相反。 對(duì)于非成像光譜， 僅有土壤樣本研磨后建立的PLSR和SVR模型能夠成功估測(cè)SOC含量， 但估測(cè)精度總是低于成像光譜建立的估測(cè)模型。

2.5 粗糙度對(duì)光譜反射率的影響

土壤在不同粗糙度下的光譜反射率特征不同。 整體表現(xiàn)為粗糙度大， 反射率低， 反之， 粗糙度小， 反射率高。 土壤樣本研磨前粗糙度高， 研磨后粗糙度低。 本研究探討了應(yīng)用不同光譜儀測(cè)量時(shí)， 土壤樣本在處于研磨前后的這兩種粗糙度狀態(tài)下的反射率特征。

成像與非成像光譜儀測(cè)量得到的研磨前的土壤光譜平均反射率總是小于研磨后的土壤光譜平均反射率。 這是由于研磨前， 土壤顆粒在團(tuán)聚作用的驅(qū)動(dòng)下形成不同大小的團(tuán)聚體， 具有非定向空間分布， 導(dǎo)致較高的粗糙度。 測(cè)量光譜時(shí)， 光散射現(xiàn)象明顯， 導(dǎo)致較低的反射率。

土壤樣本研磨前后， 成像光譜平均反射率均小于非成像光譜平均反射率。 這是因?yàn)槌上窆庾V可以獲取到土壤樣本陰影部分與非陰影部分像素點(diǎn)的反射率。 陰影部分土壤反射率低， 故而成像光譜平均反射率低于非成像光譜平均反射率。 研磨后土壤樣本陰影部分面積減少， 因此研磨后的成像光譜反射率略微低于非成像光譜反射率。 這說(shuō)明成像光譜數(shù)據(jù)能夠更加真實(shí)地反映土壤粗糙度對(duì)土壤光譜反射率的影響。

2.6 粗糙度對(duì)SOC光譜響應(yīng)波段的影響

不同的土壤粗糙度下， 可見近紅外范圍內(nèi)的光譜波段在估測(cè)SOC含量時(shí)具有不同的重要性。 變量投影重要性(variable importance in projection， VIP)指示每個(gè)光譜波段在估測(cè)SOC含量中的重要性。 圖5展示了樣本研磨前后這兩種粗糙度狀態(tài)下土壤成像光譜與非成像光譜的VIP圖。 其中VIP>1， 光譜波段被視為重要波段， 表明在估測(cè)SOC含量時(shí)起關(guān)鍵作用。 采用成像光譜技術(shù)估測(cè)SOC含量時(shí)， 樣本研磨前重要波段為614～766與902～998 nm。 其中， 766～902 nm范圍內(nèi)的光譜波段在樣本研磨前為非重要波段， 在樣本研磨后為重要波段。 利用非成像光譜建模時(shí)， 767～871 nm范圍內(nèi)的光譜波段在樣本研磨前為非重要波段， 在樣本研磨后為重要波段。 樣本研磨后重要波段為686～998 nm。 樣本研磨前后， 成像光譜與非成像光譜波段重要性曲線相似， 再一次證明了成像光譜技術(shù)被用于估測(cè)SOC含量的可行性。

圖5 土壤樣本研磨前后成像光譜(d, b)與非成像光譜(c, a)VIP圖

2.7 成像光譜技術(shù)的優(yōu)勢(shì)

將土壤樣本研磨前后的成像光譜與非成像光譜數(shù)據(jù)用于SOC含量的估測(cè)， 并對(duì)比了不同數(shù)據(jù)源及光譜預(yù)處理對(duì)PLSR和SVR模型性能的影響。 結(jié)果顯示成像光譜數(shù)據(jù)在土壤樣本研磨前后具有較為優(yōu)秀的表現(xiàn)。 而使用非成像光譜數(shù)據(jù)所構(gòu)建的模型僅能用于研磨后土壤樣本SOC含量估測(cè)， 這是由于研磨后， 降低了粗糙度， 增強(qiáng)了非成像光譜與SOC的相關(guān)性， 從而提高了PLSR模型的估測(cè)能力。 土壤樣本研磨前后， 成像光譜與SOC的相關(guān)性較強(qiáng)， 說(shuō)明成像光譜數(shù)據(jù)能夠克服土壤粗糙度的影響， 可以被用來(lái)建立良好的SOC估測(cè)模型。

此外， 利用成像光譜建立的估測(cè)模型有潛力被推廣應(yīng)用于衛(wèi)星遙感數(shù)據(jù)。 國(guó)內(nèi)現(xiàn)有的高光譜遙感衛(wèi)星有高分五號(hào)、 珠海一號(hào)等。 高分五號(hào)獲取的高光譜衛(wèi)星遙感影像在可見近紅外光譜范圍內(nèi)共有150個(gè)波段， 光譜分辨率為4.28 nm[4]。 珠海一號(hào)獲取的高光譜衛(wèi)星遙感影像在400～1 000 nm范圍內(nèi)共有32個(gè)波段， 光譜分辨率為3～8 nm[12]。 本研究獲取的成像光譜數(shù)據(jù)與高分五號(hào)、 珠海一號(hào)高光譜遙感影像具有相似的波段數(shù)量與光譜分辨率。 可以進(jìn)一步探索成像光譜與衛(wèi)星遙感影像光譜間的傳遞性， 從而將近地成像光譜SOC估測(cè)模型遷移應(yīng)用于衛(wèi)星遙感影像中， 為大尺度的SOC含量估測(cè)提供新型手段。

3 結(jié) 論

以美國(guó)愛(ài)荷華州土壤樣本為對(duì)象， 在實(shí)驗(yàn)室內(nèi)測(cè)量了不同粗糙度下土壤樣本的成像及非成像光譜反射率， 基于5種光譜預(yù)處理算法， 利用PLSR和SVR方法建立了SOC含量估測(cè)模型， 探究了利用成像光譜技術(shù)估測(cè)高粗糙度的土壤樣本的SOC含量的潛力。 結(jié)果表明成像光譜技術(shù)能夠克服土壤粗糙度的影響， 成功預(yù)測(cè)高粗糙度樣本的SOC含量。 成像光譜技術(shù)能夠同時(shí)獲取土壤空間維度與光譜維度信息， 增強(qiáng)了光譜與SOC的相關(guān)性， 因此具有良好的估測(cè)能力。 成像光譜技術(shù)作為一個(gè)新型光譜測(cè)量技術(shù)， 為野外光譜測(cè)量和大尺度SOC含量估測(cè)提供了技術(shù)支持。