基于Landsat-8數(shù)據(jù)的土壤顏色預(yù)測方法研究①

曹 振，王昌昆，馬海藝，劉 杰，徐愛愛，張芳芳，楊 穎,3，潘賢章*

基于Landsat-8數(shù)據(jù)的土壤顏色預(yù)測方法研究①

曹 振1,2，王昌昆1,2，馬海藝1,2，劉 杰1，徐愛愛1,2，張芳芳1,2，楊 穎1,2,3，潘賢章1,2*

(1 土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點(diǎn)實(shí)驗(yàn)室(中國科學(xué)院南京土壤研究所)，南京 210008；2 中國科學(xué)院大學(xué)，北京 100049；3 南京信息工程大學(xué)地理科學(xué)學(xué)院，南京 210044)

選擇黑龍江省292個自然風(fēng)干土壤樣品，室內(nèi)測定土壤高光譜反射率，然后依據(jù)三刺激值法計算土壤CIE XYZ色彩空間各顏色分量，用于土壤顏色預(yù)測和驗(yàn)證。同時，提取各樣品采集點(diǎn)位的Landsat-8 OLI原始反射率數(shù)據(jù)，計算歸一化差值植被指數(shù)、歸一化差值水體指數(shù)、歸一化差值濕度指數(shù)、歸一化差值不透水面指數(shù)，并據(jù)此提出建模光譜篩選閾值。進(jìn)一步采用偏最小二乘回歸模型，結(jié)合提取出的遙感光譜進(jìn)行各土壤顏色分量的預(yù)測。結(jié)果顯示：CIE XYZ顏色各分量的驗(yàn)證2分別為0.76、0.76和0.69，RPD分別為1.74、1.76和1.68，表明利用偏最小二乘法建立的模型可以對土壤顏色進(jìn)行近似預(yù)測。不同土地利用類型預(yù)測結(jié)果擬合顯示，耕地土壤顏色各分量預(yù)測效果均優(yōu)于林地和草地。在不同有機(jī)碳含量下分別進(jìn)行土壤顏色建模預(yù)測，當(dāng)有機(jī)碳含量較低時，土壤顏色預(yù)測結(jié)果較好。

土壤顏色；遙感數(shù)據(jù)；高光譜；偏最小二乘回歸；預(yù)測

土壤顏色是土壤的重要屬性，與土壤組成有關(guān)。Stoner和Baumgardner[1]對全美485個土壤研究發(fā)現(xiàn)，土壤顏色主要受土壤有機(jī)質(zhì)和氧化鐵不同程度的影響。章明奎等[2]指出，碳酸鹽巖發(fā)育的土壤其顏色差異性主要決定于土壤氧化鐵礦物類型及有機(jī)質(zhì)含量。張蕾等[3]發(fā)現(xiàn)氧化鐵隨著黏粒在土壤孔隙中的遷移，淀積在土壤結(jié)構(gòu)體表面上，導(dǎo)致某些黏化層結(jié)構(gòu)體表面較其內(nèi)部土壤基質(zhì)的顏色更紅。當(dāng)赤鐵礦含量較高時土壤呈現(xiàn)紅色，而針鐵礦含量較高時呈現(xiàn)黃色[4]。同時，土壤顏色還與淋溶作用有關(guān)，當(dāng)淋溶作用較強(qiáng)時，土壤顏色通常較亮[5-6]。此外，土壤顏色可以指示土壤肥力狀況。

土壤顏色的獲取方式主要有Munsell色卡目視比色法和CIE三刺激值法[7-9]。目視比色法通過人眼匹配土壤樣品和標(biāo)準(zhǔn)色卡來確定土壤顏色，簡單便捷成本低，但是存在主觀性比較強(qiáng)、系統(tǒng)性偏差比較大的問題，難以滿足土壤顏色定量化表達(dá)的需求[10-11]。三刺激值法利用土壤可見光光譜，結(jié)合標(biāo)準(zhǔn)觀察者配色函數(shù)和光源相對光譜功率分布計算土壤顏色(CIE XYZ)[12-14]，具有精度高、抗人為因素干擾強(qiáng)的特點(diǎn)，方便進(jìn)行定量化表達(dá)。雖然目前土壤顏色表達(dá)主要采用Munsell顏色空間，但是由于該系統(tǒng)將顏色劃分為一系列不連續(xù)的色塊，不適合進(jìn)行數(shù)值統(tǒng)計分析[15]。因此，Melville和Atkinson[16]對比土壤顏色表達(dá)的不同方法后，建議使用CIELAB顏色系統(tǒng)。Rossel等[17]建議將CIE XYZ顏色空間作為不同色彩之間轉(zhuǎn)換的媒介。Fernandez和Schulze[7]使用Cary 17D光譜分光儀獲取土壤反射光譜，并采用分光積分方式計算土壤CIE XYZ顏色值，證明相比于目視匹配顏色的方式，這種積分計算土壤顏色的方式更能量化土壤顏色的細(xì)微差異。

由于土壤光譜及土壤遙感數(shù)據(jù)直接包含紅、綠、藍(lán)等不同波段的地面光譜反射信息，因此，常用于土壤顏色預(yù)測。李怡春等[18]基于土壤高光譜反射率和偏最小二乘回歸(partial least squares regression，PLSR)方法進(jìn)行土壤顏色預(yù)測，并與色系轉(zhuǎn)換法進(jìn)行對比，證明了利用PLSR方法預(yù)測土壤Munsell顏色較色系轉(zhuǎn)換法更優(yōu)。Mattikalli[19]利用Landsat MSS傳感器獲取土壤反射率數(shù)據(jù)并與實(shí)測土壤Munsell顏色進(jìn)行建模擬合，結(jié)果具有很好的相關(guān)性。Escadafal[20]討論了遙感數(shù)據(jù)與土壤顏色的相關(guān)性和相關(guān)研究現(xiàn)狀，認(rèn)為利用遙感數(shù)據(jù)進(jìn)行土壤顏色的建模預(yù)測是可行的。

雖然近年來使用可見–近紅外光譜直接測量土壤顏色的方法研究越來越多，但是利用遙感數(shù)據(jù)進(jìn)行區(qū)域土壤顏色制圖還需要大量研究。如何挑選最能代表裸露地表遙感信息的反射光譜，并基于該反射光譜，進(jìn)行CIE XYZ色系的土壤顏色反演是需要解決的關(guān)鍵問題。因此，本研究基于30 m空間分辨率Landsat-8 OLI原始反射率數(shù)據(jù)，提出遙感指數(shù)的篩選條件，并進(jìn)一步采用PLSR方法進(jìn)行土壤CIE XYZ顏色各分量的預(yù)測研究，以為遙感獲取大范圍土壤顏色提供新途徑。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于黑龍江省，該地區(qū)地勢大致是西北、北部和東南部高，東北、西南部低，西部屬松嫩平原，東北部為三江平原，是我國黑土分布的核心地帶。黑龍江省屬于寒溫帶與溫帶大陸性季風(fēng)氣候，年平均降水量為526.64 mm，年平均氣溫為3℃。該區(qū)分布著多種地貌類型，包括平原、臺地、丘陵和起伏山地；土地利用類型以耕地、林地、草地為主，且耕地以旱地為主；成土母質(zhì)主要為沖積–洪積物、中性沖積物、淺色結(jié)晶巖風(fēng)化物、砂頁巖風(fēng)化物、石灰?guī)r風(fēng)化物以及黏質(zhì)黃土；土壤類型包括草甸土、黑土、暗棕壤、水稻土、黑鈣土、白漿土、沼澤土、風(fēng)沙土和棕色針葉林土。

1.2 土壤樣品采集與制備

根據(jù)地形、母質(zhì)、土地利用類型、土壤類型、交通條件等在研究區(qū)布置采樣點(diǎn)292個(圖1)，于2016年采用五點(diǎn)法采集土壤樣品，采樣時在8 m × 8 m的方形范圍內(nèi)4個頂點(diǎn)和中心位置各采1個土樣，混合后裝袋，取樣深度為0 ～ 15 cm，每個樣品1 kg。土壤樣品自然風(fēng)干后去除植物殘體，將大土塊敲碎，并取出適量樣品研磨過100目(0.149 mm)尼龍篩后備用。

1.3 土壤光譜與有機(jī)碳測定

土壤樣品高光譜采集采用美國ASD FieldSpec4 Hi-Res地物光譜儀，其波段范圍為350 ～ 2 500 nm。土壤樣品進(jìn)行室內(nèi)高光譜測定時，盡可能均勻地將土樣填充于涂黑的樣品盤中。位于樣品兩側(cè)等距位置45°角的斜上方安置兩個50 W的鹵素?zé)?，保證最大限度消除土壤顆粒陰影帶來的測量誤差。ASD探頭視場角為25°，探頭放置于樣品正上方固定高度處，保證入射視場合理地覆蓋于土壤樣品盤中心。進(jìn)行光譜測量前，提前30 min預(yù)熱光譜儀并進(jìn)行暗電流校正和白板校正，光譜測量時每隔10 min進(jìn)行一次白板校正，白板校正使用的是25 cm×25 cm的漫反射標(biāo)準(zhǔn)參考板。對于每個土壤樣品，為保證測量結(jié)果不受土壤表面狀態(tài)的干擾，在依次旋轉(zhuǎn)90°的4個方向上各測量5條光譜曲線，共20條光譜曲線，以這20條光譜曲線的平均值作為樣品的測定值[1]。通過對光譜曲線進(jìn)行拼接校正，獲得土壤樣品在不同波段的光譜反射率值，本研究選取的波段范圍位于可見光部分380 ～ 780 nm。

土壤樣品有機(jī)碳含量基本信息見表1，測定方法參照《土壤農(nóng)業(yè)化學(xué)分析方法》[21]。

(該圖基于國家測繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下載的審圖號為GS(2019)3333的標(biāo)準(zhǔn)地圖制作，底圖無修改)

表1 土壤有機(jī)碳含量統(tǒng)計信息

1.4 土壤顏色獲取

土壤顏色利用CIE XYZ顏色空間模型來表示，CIE XYZ顏色空間是由國際照明委員會在1931年提出的，其中Y代表亮度，X和Z是主光譜的虛擬組成部分。三刺激值XYZ計算公式為[22]：

1.5 遙感數(shù)據(jù)的獲取與處理

Landsat-8遙感影像數(shù)據(jù)集來自于Google Earth Engine(GEE)，采用經(jīng)過大氣校正和正射校正的Landsat-8表面反射率數(shù)據(jù)產(chǎn)品(Surface Reflectance Tier 1)[25]。本研究采用Landsat-8 OLI遙感影像數(shù)據(jù)集中的可見–近紅外波段以及短波紅外波段(共7個波段)進(jìn)行土壤顏色預(yù)測?？紤]到遙感影像由于云、雪等對下墊面遮擋造成的遙感數(shù)據(jù)缺失以及土壤顏色在短期內(nèi)變化幅度很小兩個因素，對使用遙感光譜數(shù)據(jù)的時間進(jìn)行適當(dāng)放寬，在土壤樣品采集年份(2016年)前后各延伸1年，最終選取了黑龍江省2015—2017年遙感影像，并通過質(zhì)量評估波段去除包含云、云陰影、雪等干擾的影像，獲得了所有土壤樣點(diǎn)在1 248景遙感影像上的共9 885條光譜數(shù)據(jù)。

本研究首先提出篩選建模光譜的遙感指數(shù)條件，然后通過遙感指數(shù)構(gòu)建建模光譜數(shù)據(jù)的篩選條件。選用的遙感指數(shù)分別為歸一化差值植被指數(shù)(normalized difference vegetation index, NDVI)、歸一化差值水體指數(shù)(normalized difference water index, NDWI)、歸一化差值濕度指數(shù)(normalized difference moisture index, NDMI)、歸一化差值不透水面指數(shù)(normalized difference impervious surface index, NDISI)[26-27]。各指數(shù)具體計算公式為：

式中：3、4、5、6、10分別與Landsat-8 OLI遙感影像數(shù)據(jù)集波段相對應(yīng)，代表綠色、紅色、近紅外、短波紅外、熱紅外波段。通過對遙感影像光譜反射率數(shù)據(jù)和土壤樣品實(shí)測反射率數(shù)據(jù)的對比分析，確定了篩選光譜約束條件為：

對于每個土壤樣點(diǎn)，滿足所有篩選條件的遙感影像光譜通常不只一條，還需進(jìn)一步選取。先計算各波段光譜平均值，然后計算篩選出來的各條光譜與該平均光譜的歐氏距離，并選擇歐氏距離最小的一條光譜曲線作為對應(yīng)點(diǎn)用來進(jìn)行顏色預(yù)測的光譜數(shù)據(jù)。歐氏距離計算公式為：

1.6 建模與評價方法

本研究采用PLSR模型對土壤顏色進(jìn)行建模預(yù)測。PLSR算法是一種多元線性回歸建模方法，在建模過程中集成了主成分分析、典型相關(guān)分析和線性回歸分析方法的特點(diǎn)，能夠有效處理變量間存在多重相關(guān)性而樣本量較少的問題[28]。主成分?jǐn)?shù)選取對PLSR建模尤為重要，主成分?jǐn)?shù)過少，沒有足夠信息進(jìn)入模型擬合，會導(dǎo)致擬合不充分，降低模型的預(yù)測精度；主成分?jǐn)?shù)過多，很多無用信息會進(jìn)入模型，容易出現(xiàn)過擬合，同樣會降低模型預(yù)測精度。本文利用赤池信息量準(zhǔn)則值(akaike information criterion, AIC)來選擇模型的最佳主成分?jǐn)?shù)：

AIC=lnRMSE+2(5)

式中：是建模集樣本數(shù)量；RMSE為驗(yàn)證均方根誤差；為對應(yīng)的主成分?jǐn)?shù)，以AIC最小時的主成分?jǐn)?shù)作為PLSR模型的最優(yōu)主成分?jǐn)?shù)[29]。評價建模精度時使用驗(yàn)證決定系數(shù)2、RMSE和RPD(測定值標(biāo)準(zhǔn)偏差與標(biāo)準(zhǔn)預(yù)測誤差的比值)：

式中：SD表示數(shù)據(jù)集實(shí)測值標(biāo)準(zhǔn)差。已有研究結(jié)果表明，當(dāng)RPD≥2.50時，模型可以精確預(yù)測；當(dāng)2.50>RPD≥2.00時，模型可以很好地預(yù)測；當(dāng)2.00>RPD≥1.80時，模型可以較好地預(yù)測；當(dāng)1.80>RPD≥1.40時，模型可以近似地預(yù)測；當(dāng)1.40>RPD≥1.00時，模型可以區(qū)分變量值的大?。划?dāng)RPD<1.00時，模型不適合預(yù)測[29]。

建模過程中數(shù)據(jù)集被隨機(jī)劃分成70% 的建模集(204個土樣)和30% 的驗(yàn)證集(88個土樣)，并選擇Landsat-8遙感產(chǎn)品數(shù)據(jù)集可見–近紅外、短波紅外共7個波段作為建模變量，因變量分別為CIE XYZ三種顏色分量。根據(jù)AIC標(biāo)準(zhǔn)，PLSR方法預(yù)測CIE XYZ顏色各分量的最佳主成分?jǐn)?shù)參數(shù)設(shè)定為5、5、6，其余參數(shù)為默認(rèn)值。

文中涉及的光譜值處理、色系值等計算在Python 3.6中完成，PLSR建模采用基于Python語言編制的scikit-learn機(jī)器學(xué)習(xí)庫[30]完成，圖表由Excel 2013、ArcGIS 10.2制作完成。

2 結(jié)果

2.1 土壤顏色實(shí)測結(jié)果統(tǒng)計特征

實(shí)測土壤樣品CIE XYZ顏色結(jié)果如圖2所示，各顏色分量呈近似正態(tài)分布。其中，X值可以指示人眼對紅色色光的感應(yīng)強(qiáng)弱，其均值為9.11，標(biāo)準(zhǔn)差為2.52，最大值為14.63，最小值為4.57，總體跨度為10.06；Y值可以指示人眼對綠色色光的感應(yīng)強(qiáng)弱，同時反映著土壤顏色亮度大小，其均值為9.01，標(biāo)準(zhǔn)差為2.47，最大值為14.48，最小值為4.59，總體跨度為9.89，表明土壤顏色總體亮度偏低；Z值可以指示人眼對藍(lán)色色光的感應(yīng)強(qiáng)弱，其均值為7.33，標(biāo)準(zhǔn)差為1.73，最大值為11.57，最小值為4.33，總體跨度為7.24?？傮w來看，CIE XYZ顏色各分量跨度不大，均值偏低，這是因?yàn)橥寥罉悠繁旧砉庾V反射率較低，加之樣品中有機(jī)碳含量豐富，導(dǎo)致土壤顏色偏暗。

2.2 土壤顏色預(yù)測結(jié)果

2.2.1 光譜篩選對土壤顏色預(yù)測的影響 光譜篩選前后土壤顏色預(yù)測結(jié)果對比如圖3所示。基于篩選后的建模光譜進(jìn)行土壤顏色的建模預(yù)測，結(jié)果如圖3A所示。在X顏色分量上，PLSR預(yù)測結(jié)果2、RMSE、RPD分別為0.76、0.94和1.74，說明該模型可以近似地預(yù)測X顏色分量；在Y顏色分量上，PLSR預(yù)測結(jié)果2、RMSE、RPD分別為0.76、0.91和1.76，說明該模型可以近似地預(yù)測Y顏色分量；在Z顏色分量上，PLSR預(yù)測結(jié)果2、RMSE、RPD分別為0.69、0.66和1.68，說明該模型也可以近似地預(yù)測Z顏色分量?？傮w來看，基于PLSR方法進(jìn)行建模預(yù)測3種土壤顏色分量是可行的，均可實(shí)現(xiàn)近似預(yù)測。為了比較光譜篩選前后土壤顏色預(yù)測精度，本文對光譜篩選前遙感數(shù)據(jù)進(jìn)行建模預(yù)測。圖3B顯示了進(jìn)行光譜篩選前的土壤顏色的建模預(yù)測結(jié)果，可見土壤顏色CIE XYZ各分量PLSR預(yù)測結(jié)果RPD分別為1.26、1.27和1.26，均低于光譜篩選后的建模預(yù)測結(jié)果，說明進(jìn)行本文光譜篩選能夠有效提升土壤顏色的建模預(yù)測精度。

圖2 土壤樣品實(shí)測CIE XYZ顏色頻率分布圖

圖3 PLSR預(yù)測值與實(shí)測值對比圖

2.2.2 土地利用對土壤顏色預(yù)測的影響 為了探究土地利用類型對土壤顏色預(yù)測精度的影響，按照不同土地利用類型分別對預(yù)測結(jié)果進(jìn)行擬合分析，結(jié)果如圖4所示。在土壤顏色CIE XYZ各顏色分量上，耕地類型的預(yù)測結(jié)果均優(yōu)于草地和林地類型的預(yù)測結(jié)果。其中，耕地類型土壤顏色預(yù)測RPD分別為2.42、2.45和2.21，這表明模型能夠?qū)Ω赝寥李伾M(jìn)行很好的預(yù)測；林地類型土壤顏色預(yù)測RPD分別為1.69、1.69和1.42，這表明模型能夠近似預(yù)測林地土壤顏色；草地類型土壤顏色預(yù)測RPD分別為0.98、1.01和1.28，表明模型對于草地土壤顏色的預(yù)測效果較差?？梢?，在不同土地利用類型下土壤顏色的預(yù)測結(jié)果具有明顯的差異性。這是由于耕地土壤反射光譜較為一致，而林地、草地由于地表覆蓋、腐殖質(zhì)積累等差異較大，導(dǎo)致反射光譜中混雜有較多的噪聲信息。此外，本研究中林地和草地的樣本數(shù)量相對偏少，這也會導(dǎo)致土壤顏色的預(yù)測結(jié)果出現(xiàn)這種差異性。

2.2.3 有機(jī)碳含量對土壤顏色預(yù)測的影響 為了探究不同有機(jī)碳含量對土壤顏色建模預(yù)測結(jié)果的影響，進(jìn)一步按照不同有機(jī)碳含量將驗(yàn)證集土壤樣本劃分成高值區(qū)、中值區(qū)和低值區(qū)3個子集，并分別對3個子集進(jìn)行建模預(yù)測，結(jié)果如圖5所示。其中，高值區(qū)子集有機(jī)碳含量范圍為28.25 ～ 88.92 g/kg，共30個土壤樣品；中值區(qū)子集有機(jī)碳含量范圍為21.00 ～ 28.19 g/kg，共29個土壤樣品；低值區(qū)子集有機(jī)碳含量范圍為7.42 ～ 20.82 g/kg，共29個土壤樣品。分別對有機(jī)碳含量高、中和低值區(qū)3個子集的土壤顏色建模預(yù)測結(jié)果進(jìn)行驗(yàn)證，結(jié)果顯示，在CIE X顏色分量上，預(yù)測土壤顏色驗(yàn)證RPD分別為1.39、1.50和1.77；在CIE Y顏色分量上，預(yù)測土壤顏色驗(yàn)證RPD分別為1.35、1.55和1.82；在CIE Z顏色分量上，預(yù)測土壤顏色驗(yàn)證RPD分別為1.02、1.97和2.07?？傮w來說，當(dāng)土壤有機(jī)碳含量較低時，各顏色分量驗(yàn)證RPD較高，預(yù)測效果較好。

圖4 不同土地利用類型土壤顏色預(yù)測值與實(shí)測值對比圖

圖5 不同有機(jī)碳含量土壤顏色預(yù)測值與實(shí)測值對比圖

3 討論

本文基于土壤反射光譜，依據(jù)三刺激值法計算土壤CIE XYZ色彩空間各顏色分量值，并結(jié)合遙感數(shù)據(jù)進(jìn)行土壤顏色建模預(yù)測。首先定義了光譜數(shù)據(jù)篩選方法對遙感光譜數(shù)據(jù)進(jìn)行篩選，然后采用PLSR方法對土壤顏色各顏色分量進(jìn)行建模預(yù)測，結(jié)果顯示，CIE XYZ顏色各分量的驗(yàn)證RPD分別為1.74、1.76和1.68，表明本方法建立的模型可以對土壤顏色進(jìn)行可近似預(yù)測。土壤顏色的預(yù)測受多種因素影響，本文從光譜篩選、土地利用類型和有機(jī)碳含量3個方面進(jìn)行討論。

從光譜篩選過程來看，本研究通過確定遙感指數(shù)NDVI、NDWI、NDMI、NDISI的閾值，進(jìn)行建模光譜數(shù)據(jù)的篩選。其中，NDVI閾值主要消除植被信息，NDWI、NDMI閾值主要消除水分信息，NDISI閾值主要消除不透水地面信息。對比光譜篩選前后的預(yù)測結(jié)果，可以看見土壤顏色預(yù)測結(jié)果RPD由光譜篩選前的1.26、1.27和1.26分別提升至1.74、1.76和1.68，這說明本文提出的光譜篩選方法能夠有效提升土壤顏色的預(yù)測精度，通過篩選獲得的光譜，可以在最大程度上排除植被等外部參數(shù)的干擾。實(shí)際上，衛(wèi)星過境與地面純裸露時間不一定正好匹配，即便是裸露地表也或多或少存在諸如腐爛秸稈、稀疏植被、動植物殘體等的覆蓋，不可能獲得像在實(shí)驗(yàn)室條件下制作非常干凈的裸土光譜，因此，進(jìn)一步消除植被等外部影響的干擾，提升遙感數(shù)據(jù)預(yù)測土壤顏色的精度是有可能的。

從土地利用類型影響來看，對比不同土地利用類型下土壤顏色的預(yù)測結(jié)果，發(fā)現(xiàn)存在差異，耕地預(yù)測效果要優(yōu)于林地、草地。耕地類型土壤顏色CIE XYZ各分量值預(yù)測RPD分別為2.42、2.45和2.21，這表明模型能夠?qū)Ω赝寥李伾M(jìn)行很好的預(yù)測。這是因?yàn)楦赝寥婪瓷涔庾V較為一致，而林地、草地由于地表覆蓋、腐殖質(zhì)積累等差異較大，反射光譜中混雜有較多的噪聲信息。此外，由于本研究中林地和草地的樣本數(shù)量相對偏少，在土壤顏色建模預(yù)測的過程中偏差相對較大，這也會導(dǎo)致土壤顏色的預(yù)測結(jié)果較差。

從有機(jī)碳含量影響來看，Stoner和Baumgar-dner[1]研究表明，土壤顏色主要受土壤有機(jī)碳和氧化鐵組分影響。本研究采集的土壤樣品分布于黑龍江省全境，土壤類型以草甸土、黑土、暗棕壤和水稻土居多，有機(jī)碳含量豐富。土壤樣品顏色主要受有機(jī)碳含量影響。本研究中，不同有機(jī)碳含量土壤顏色預(yù)測結(jié)果表明，有機(jī)碳含量較低時，土壤顏色預(yù)測結(jié)果較好。出現(xiàn)這種結(jié)果可能與土壤本身反射光譜的能力有關(guān)，當(dāng)有機(jī)碳含量較低時，土壤顏色較亮，反射光譜能力較強(qiáng)，反射光譜中包含更多的有效信息，預(yù)測結(jié)果較好。

本研究中土壤樣品全部來自于黑龍江省，具有一定的區(qū)域性，土樣采集涵蓋當(dāng)?shù)氐湫屯恋乩妙愋秃屯寥李愋?，土壤顏色主要受有機(jī)碳影響。對于其他受不同土壤顏色主導(dǎo)因素影響的土壤如紅壤、紫色土等，是否具有相同結(jié)果還需進(jìn)一步驗(yàn)證。

4 結(jié)論

本研究基于土壤實(shí)測反射光譜和Landsat-8遙感數(shù)據(jù)進(jìn)行土壤顏色建模預(yù)測，并提出了建模光譜數(shù)據(jù)篩選方法。采用PLSR建模預(yù)測CIE XYZ顏色各分量，驗(yàn)證結(jié)果RPD分別為1.74、1.76和1.68，說明按本方法建立的模型可以對X、Y、Z各分量進(jìn)行近似預(yù)測，且耕地土壤顏色預(yù)測效果要優(yōu)于林地和草地。分別對有機(jī)碳含量高、中和低值區(qū)土壤顏色進(jìn)行預(yù)測，各顏色分量的驗(yàn)證RPD隨著有機(jī)碳含量的降低而升高，說明當(dāng)有機(jī)碳含量較低時，土壤顏色預(yù)測結(jié)果較好。

[1] Stoner E R, Baumgardner M F. Characteristic variations in reflectance of surface soils[J]. Soil Science Society of America Journal, 1981, 45(6): 1161–1165.

[2] 章明奎, 姚玉才, 邱志騰, 等. 中國南方碳酸鹽巖發(fā)育土壤的成土特點(diǎn)與系統(tǒng)分類[J]. 浙江大學(xué)學(xué)報(農(nóng)業(yè)與生命科學(xué)版), 2019, 45(1): 54–65.

[3] 張蕾, 張鳳榮, 靳東升, 等. 黃土高原淋溶土黏粒、氧化鐵與顏色的關(guān)系及發(fā)生學(xué)解釋——以山西土系調(diào)查的31個黏化層為例[J]. 土壤學(xué)報, 2021, 58(4): 876–886.

[4] Hu X F, Du Y, Guan C L, et al. Color variations of the quaternary red clay in Southern China and its paleoclimatic implications[J]. Sedimentary Geology, 2014, 303: 15–25.

[5] 邱堋星, 葉飛. 武夷山土壤性狀及其垂直變化規(guī)律[J]. 東北林業(yè)大學(xué)學(xué)報, 2008, 36(1): 55–57.

[6] 熊毅, 李慶逵. 中國土壤[M]. 2版. 北京: 科學(xué)出版社, 1987.

[7] Fernandez R N, Schulze D G. Calculation of soil color from reflectance spectra[J]. Soil Science Society of America Journal, 1987, 51(5): 1277–1282.

[8] Torrent J, Barrón V. Laboratory measurement of soil color: Theory and practice[M]//Soil Color. Madison, WI, USA: Soil Science Society of America, 2015: 21–33.

[9] 陳劍科, 袁大剛, 晏昭敏, 等. 測色儀與中國標(biāo)準(zhǔn)土壤色卡測定土壤顏色比較——以川中丘陵區(qū)為例[J]. 土壤學(xué)報, 2019, 56(1): 78–89.

[10] Moritsuka N, Kawamura K, Tsujimoto Y, et al. Comparison of visual and instrumental measurements of soil color with different low-cost colorimeters[J]. Soil Science and Plant Nutrition, 2019, 65(6): 605–615.

[11] 岳智慧, 黃強(qiáng), 肖理, 等. 土壤顏色由CIE向Munsell系統(tǒng)的定量轉(zhuǎn)換[J]. 光譜學(xué)與光譜分析, 2019, 39(9): 2842–2846.

[12] 荊其誠. 色度學(xué)[M]. 北京: 科學(xué)出版社, 1979.

[13] 吳豪翔. 土壤光譜的顏色計算及其應(yīng)用意義[J]. 科技通報, 1989, 5(2): 32–34.

[14] 徐彬彬. 土壤剖面的反射光譜研究[J]. 土壤, 2000, 32(6): 281–287.

[15] Hartemink A E, Minasny B. Towards digital soil morphometrics[J]. Geoderma, 2014, 230/231: 305–317.

[16] Melville M D, Atkinson G. Soil colour: Its measurement and its designation in models of uniform colour space[J]. Journal of Soil Science, 1985, 36(4): 495–512.

[17] Viscarra Rossel R A, Minasny B, Roudier P, et al. Colour space models for soil science[J]. Geoderma, 2006, 133(3/4): 320–337.

[18] 李怡春, 潘愷, 王昌昆, 等. 基于PLSR的土壤顏色預(yù)測方法及其與色系轉(zhuǎn)換法的對比研究[J]. 土壤學(xué)報, 2018, 55(6): 1411–1421.

[19] Mattikalli N M. Soil color modeling for the visible and near-infrared bands of Landsat sensors using laboratory spectral measurements[J]. Remote Sensing of Environment, 1997, 59(1): 14–28.

[20] Escadafal R. Remote sensing of soil color: Principles and applications[J]. Remote Sensing Reviews, 1993, 7(3/4): 261–279.

[21] 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京: 中國農(nóng)業(yè)科技出版社, 2000.

[22] ASTM Intl. E308-18, Standard practice for computing the colors of objects by using the CIE system[S]. West Conshohocken, Philadclphia: American Society for Testing and Materials, 2018.

[23] 胡威捷, 湯順青, 朱正芳. 現(xiàn)代顏色技術(shù)原理及應(yīng)用: 光學(xué)工程[M]. 北京: 北京理工大學(xué)出版社, 2007.

[24] 徐彬彬, 戴昌達(dá). 土壤顏色的分光光度計算法[J]. 土壤學(xué)報, 1979, 16(4): 397–408.

[25] Gorelick N, Hancher M, Dixon M, et al. Google earth engine: Planetary-scale geospatial analysis for everyone[J]. Remote Sensing of Environment, 2017, 202: 18–27.

[26] 徐涵秋. 福建省長汀縣河田盆地區(qū)近35年來地表裸土變化的遙感時空分析[J]. 生態(tài)學(xué)報, 2013, 33(10): 2946– 2953.

[27] Xu H Q. Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery[J]. International Journal of Remote Sensing, 2006, 27(14): 3025–3033.

[28] 司守奎, 孫璽菁. 數(shù)學(xué)建模算法與應(yīng)用[M]. 北京: 國防工業(yè)出版社, 2011.

[29] Viscarra Rossel R A, McGlynn R N, McBratney A B. Determining the composition of mineral-organic mixes using UV-vis-NIR diffuse reflectance spectroscopy[J]. Geoderma, 2006, 137(1/2): 70–82.

[30] Pedregosa F, Varoquaux G, Gramfort A, et al. Scikit-learn: Machine learning in python[J]. Journal of Machine Learning Research, 2011, 12: 2825–2830.

Soil Color Prediction Method Based on Landsat-8 Data

CAO Zhen1, 2, WANG Changkun1, 2, MA Haiyi1, 2, LIU Jie1, XU Aiai1, 2, ZHANG Fangfang1, 2, YANG Ying1, 2, 3, PAN Xianzhang1, 2*

(1 State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China; 2 University of Chinese Academy of Sciences, Beijing 100049, China; 3 School of Geographical Sciences, Nanjing University of Information Science and Technology, Nanjing 210044, China)

A total of 292 air-dried soil samples were used from Heilongjiang Province, and the hyper-spectral reflectance data were also measured in the laboratory. Meanwhile, color components of all soil samples were calculated in the CIE XYZ color space according to the three-stimulus method to predict and validate soil color. Then, Landsat-8 OLI original reflectance data of each soil sample site were extracted to calculate the normalized difference vegetation index, normalized difference water index, normalized difference moisture index and normalized difference impervious surface index. Based on these remote sensing indexes, the threshold value of modeling spectra screening is proposed. Combined with the extracted remote sensing spectra, the partial least squares regression model was used to predict soil color components. The results showed that for the color components of CIE X, Y and Z, the validation2value were 0.76, 0.76 and 0.69, and the RPD value were 1.74, 1.76 and 1.68, indicating that the model established by partial least squares can be used to predict soil color. The prediction of different land use types showed that cultivated land soil color was better in prediction than those of forest and grassland. Soil color was also predicted under different organic carbon contents, showing the prediction was better in lower organic carbon content.

Soil color; Remote sensing; Hyper-spectral; Partial least squares regression; Prediction

曹振, 王昌昆, 馬海藝, 等. 基于Landsat-8數(shù)據(jù)的土壤顏色預(yù)測方法研究. 土壤, 2022, 54(1): 152–160.

S159.9

A

10.13758/j.cnki.tr.2022.01.020

國家重點(diǎn)研發(fā)計劃項(xiàng)目(2020YFC1807401，2018YFC1800104)和中國科學(xué)院野外站聯(lián)盟項(xiàng)目(KFJ-SW-YW035-3)資助。

(panxz@issas.ac.cn)

曹振(1995—)，男，河南信陽人，碩士研究生，主要從事土壤遙感應(yīng)用研究。E-mail: caozhen@issas.ac.cn