譚 潔，陳 嚴(yán)，周衛(wèi)軍，崔浩杰，劉 沛

基于實(shí)驗(yàn)室高光譜數(shù)據(jù)的大圍山森林土壤氧化鐵全量反演①

譚 潔1，陳 嚴(yán)1，周衛(wèi)軍2*，崔浩杰2，劉 沛1

(1湖南農(nóng)業(yè)大學(xué)風(fēng)景園林與藝術(shù)設(shè)計(jì)學(xué)院，長(zhǎng)沙 410128；2湖南農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院，長(zhǎng)沙 410128)

氧化鐵是土壤中含鐵礦物的主體，是土壤發(fā)育和土壤分類最明顯和最有用的指標(biāo)之一。本文以湖南省大圍山森林土壤為研究對(duì)象，通過(guò)實(shí)驗(yàn)室化學(xué)成分測(cè)定和光譜采集，在光譜預(yù)處理及組合變換基礎(chǔ)上，采用相關(guān)性分析篩選土壤氧化鐵全量的敏感波段，并分別建立多元逐步回歸和偏最小二乘回歸反演模型。結(jié)果表明：不同土壤光譜曲線趨勢(shì)基本一致，均形似陡坎，且在420 ～ 580 nm波段，土壤氧化鐵全量與光譜反射率呈負(fù)相關(guān)關(guān)系；不同的光譜數(shù)據(jù)變換方式可以提高光譜與氧化鐵全量的相關(guān)性，Savitzky-Golay(S-G)平滑和去包絡(luò)線相結(jié)合優(yōu)于其他預(yù)處理方法；土壤氧化鐵全量的特征波段主要為392、427、529、523、549、559、565、570、994和1 040 nm，偏最小二乘回歸模型比多元逐步回歸模型具有更好的穩(wěn)定性，適合于快速估算紅黃壤區(qū)森林土壤氧化鐵全量。

土壤光譜；氧化鐵；多元逐步回歸；偏最小二乘回歸

鐵是巖石圈的主要元素之一，是土壤中重要的微量元素，不僅對(duì)植物的生長(zhǎng)發(fā)育和農(nóng)產(chǎn)品品質(zhì)具有極其重要的作用，還對(duì)人體健康有著重要的意義。土壤中的氧化鐵是在地表特定水熱條件下，由硅酸鹽類礦物經(jīng)徹底風(fēng)化作用而形成的重要土壤物質(zhì)。由于土壤氧化鐵的聚積和遷移活動(dòng)受到化學(xué)風(fēng)化、生物循環(huán)等因素影響，反映了土壤的淋溶過(guò)程、風(fēng)化發(fā)育程度以及土壤的地帶性分布特征，常被作為描述土壤發(fā)育和土壤分類的最明顯和最有用的指標(biāo)之一[1-3]。傳統(tǒng)土壤化學(xué)分析方法存在著數(shù)據(jù)獲取難、分析成本高、時(shí)效性差等問題，隨著土壤紅外光譜[4]、漫反射光譜[5]、高光譜[6]等光譜技術(shù)的快速發(fā)展，為快速、準(zhǔn)確、定量獲取土壤氧化鐵全量信息提供了新的科學(xué)途徑。

土壤是由不同理化性質(zhì)的物質(zhì)混合組成的系統(tǒng)，其高光譜特性受到有機(jī)質(zhì)、水分、黏土礦物、氧化鐵、母質(zhì)及土壤質(zhì)地等因素共同作用。氧化鐵在土壤中常處于還原淋溶和氧化淀積的交替過(guò)程，氧化鐵的形態(tài)、顆粒粗細(xì)以及是否以膜狀態(tài)出現(xiàn)，對(duì)土壤顏色產(chǎn)生了很大的影響，從而也顯著影響了土壤光譜特征[7]。Stone等人[8]的研究認(rèn)為氧化鐵的吸收帶主要集中在反射光譜曲線的700 nm和 900 nm處；Hunt和Ashley[9]研究發(fā)現(xiàn)Fe2+在1 000 ～ 1 100 nm 附近可產(chǎn)生一個(gè)最常見的強(qiáng)而寬的譜帶，F(xiàn)e3+形成的吸收帶主要位于 900 nm 附近；Galv?o等人[10]研究巴西熱帶土壤時(shí)發(fā)現(xiàn)土壤中氧化鐵含量與 1 710 nm 處的反射率的相關(guān)性顯著。何挺等[11]的研究發(fā)現(xiàn)土壤氧化鐵含量與1 550 ～ 2 320 nm 波段反射率的相關(guān)性很強(qiáng)，且含量越高，在400 ～ 1 200 nm 波段的反射率越低，利用紅波段和藍(lán)波段平均反射率的對(duì)數(shù)構(gòu)成的土壤氧化鐵指數(shù)，在高光譜圖像解譯中具有應(yīng)用潛力；彭杰等[12]認(rèn)為隨游離氧化鐵含量的增加，350 ～ 570 nm 波段土壤反射率降低，570 ～ 2 500 nm 波段的反射率增加?；趯?duì)土壤氧化鐵的高光譜特征分析，學(xué)者們還采用多元逐步回歸模型[13]、神經(jīng)網(wǎng)絡(luò)模型(BPNN)[14]、偏最小二乘回歸模型(PLSR)[15]、多層感知器神經(jīng)網(wǎng)絡(luò)模型(MLP)[16]、隨機(jī)森林回歸算法(RFR)[17]、組合模型[18]和模糊識(shí)別[19]等方法，分別估算了土壤氧化鐵全量、游離鐵、針鐵礦和赤鐵礦含量。由此可見，土壤氧化鐵全量對(duì)土壤光譜存在明顯影響，可借助高光譜技術(shù)對(duì)土壤氧化鐵全量進(jìn)行估測(cè)。然而，目前針對(duì)森林土壤氧化鐵全量的光譜反演研究相對(duì)較少，且高光譜測(cè)量易受到各種環(huán)境因素的影響，在一定程度上增加了微量元素反演的難度。因此，實(shí)驗(yàn)室的可見光到近紅外區(qū)域高光譜反射率用于預(yù)測(cè)土壤氧化鐵全量的潛力還需要有效深入地挖掘。

本研究以湖南省大圍山為研究區(qū)，利用室內(nèi)高光譜數(shù)據(jù)，通過(guò)運(yùn)用多種平滑和數(shù)學(xué)變換形式對(duì)原始光譜進(jìn)行預(yù)處理，分析土壤氧化鐵全量的高光譜特征和特征波段，厘清多元逐步回歸模型和偏最小二乘模型估算土壤氧化鐵全量的精度差異，明確模型適用性以及最佳預(yù)處理方法，為快速定量獲取紅黃壤區(qū)森林土壤氧化鐵全量提供高光譜模型依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

大圍山(114°2′ ～ 114°12′ E，28°21′ ～ 28°26′ N)位于湖南省瀏陽(yáng)市東北部大圍山鎮(zhèn)與張坊鎮(zhèn)交界處，與江西省銅鼓縣接壤，屬中亞熱帶季風(fēng)濕潤(rùn)氣候，年平均氣溫11 ℃ ～ 17 ℃，年平均降雨量1 200 ～ 2 000 mm，年相對(duì)濕度高于83%。研究區(qū)屬中山地貌，海拔范圍為230 ～ 1 607.9 m，相對(duì)高差1 377.9 m，地質(zhì)單元主要為燕山期花崗巖(γ)和中元古代冷家溪群板頁(yè)巖(Ptln)[20]，山地土壤和植被具有明顯垂直地帶性。

1.2 樣品采集與測(cè)試

本研究采集的樣品為大圍山花崗巖發(fā)育的紅壤、黃紅壤、黃壤、黃棕壤。在海拔177、474、545、720、860、1 078、1 210、1 382 m選擇了8個(gè)具有代表性的土壤剖面設(shè)置樣地，按土壤發(fā)生學(xué)劃分原采集各土壤剖面樣品共計(jì) 26 個(gè)，樣品經(jīng)自然風(fēng)干，檢出肉眼可見的沙礫及植物殘?bào)w，研磨過(guò)100目篩。采用四分法將每個(gè)樣品分為2份，1份采用三酸消解–原子吸收分光光度法測(cè)定土壤氧化鐵全量[21]，另一份用于土壤光譜測(cè)量。

在人工暗室中利用SOC-710(surface optics cor-poration，USA)成像高光譜儀對(duì)土壤樣品進(jìn)行了高光譜圖像采集。儀器光譜范圍為400 ～ 1 050 nm，光譜分辨率為1.3 nm，128個(gè)波段。將土壤樣品平鋪至黑色塑料盒(直徑10 cm，深1 cm)，水平放置在成像高光譜儀的視野中，灰板也放置于同視場(chǎng)角內(nèi)，寬光譜石英鹵素?zé)舸怪闭障蛲寥罉悠?，與土樣表面距離約100 cm。高光譜儀鏡頭光圈設(shè)置為5.6 mm，調(diào)整焦距直至顯示出清晰土壤樣品圖像，曝光時(shí)間為自動(dòng)，采集高光譜數(shù)據(jù)。通過(guò)SRAnal710e軟件將土壤光譜數(shù)據(jù)進(jìn)行輻射定標(biāo)并將能量值轉(zhuǎn)化成真實(shí)反射率，借助ENVI軟件在不同的鋁盒中分別選取靠近鋁盒中心的50個(gè)像素點(diǎn)作為ROI區(qū)域建立波譜庫(kù)。

1.3 高光譜數(shù)據(jù)處理與特征波段提取

將 26個(gè)土壤樣品分別依照土壤氧化鐵全量進(jìn)行聚類分析，并計(jì)算出每一類土壤氧化鐵全量的光譜均值，獲得土壤氧化鐵全量的光譜響應(yīng)范圍。將土壤氧化鐵全量分別與原始光譜反射率、數(shù)學(xué)變換光譜反射率和去包絡(luò)線光譜反射率進(jìn)行逐波段的相關(guān)分析，根據(jù)相關(guān)系數(shù)最值對(duì)應(yīng)波長(zhǎng)選擇特征波段。相關(guān)系數(shù)計(jì)算公式如下：

按照相關(guān)系數(shù)的顯著程度可以分成以下3類：|y|<0.4為低度線性相關(guān)，0.4<|y|<0.7為顯著性相關(guān)，0.7<|y|<1為高度線性相關(guān)[23]。

1.4 光譜預(yù)測(cè)模型與驗(yàn)證

為保證預(yù)測(cè)樣本和驗(yàn)證樣本的不同土壤類型的土壤氧化鐵全量變化范圍一致，將不同土壤類型的土壤，根據(jù)氧化鐵全量的不同，劃分為若干個(gè)小范圍，統(tǒng)計(jì)每個(gè)范圍內(nèi)的樣本數(shù)后，隨機(jī)選取20個(gè)建模樣本，6 個(gè)驗(yàn)證樣本。基于去包絡(luò)線全波段(380 ～ 1 050 nm)高光譜數(shù)據(jù)建立多元線性回歸估測(cè)模型(MLR)；按照建模組與驗(yàn)證組比例分別為0.5:0.5和0.75:0.25，建立了兩種土壤偏最小二乘回歸估測(cè)模型(PLSR)。模型精度評(píng)價(jià)為決定系數(shù)2、校正均方根誤差RMSEC、預(yù)測(cè)均方根RMSEP，2越大，RMSEC越小，RMSEP與RMSEC越接近，表明模型的校正及預(yù)測(cè)精度越高。根據(jù)Rossel和Webster[24]對(duì)殘余預(yù)測(cè)偏差的分析可知，RPD<1.4說(shuō)明模型的預(yù)測(cè)能力很差；1.4≤RPD<2.0，表明模型初步具備定量預(yù)測(cè)可能；RPD≥2.0表明模型具備較好的定量預(yù)測(cè)能力。本研究PLSR 建模和驗(yàn)證利用Unscrambler?X10.4完成，多元線性回歸建模、數(shù)據(jù)分析在 SPSS 22.0 完成，繪圖通過(guò)GraphPad Prism 8完成。

2 結(jié)果與分析

2.1 土壤氧化鐵全量變化特征

由于海拔高低相差大，氣候條件差異較大，大圍山土壤呈明顯垂直地帶分布特征，低于海拔500 m多為山地紅壤，海拔500 ～ 800 m為山地黃紅壤，海拔800 ～ 1 200 m為山地黃壤，山地黃棕壤則主要分布在海拔1 200 ～ 1 600 m[25]。由圖1和圖2可知，8個(gè)土壤剖面的氧化鐵全量均較高，均值為46.33 ～ 79.61 g/kg，表現(xiàn)為紅壤>黃紅壤>黃壤>黃棕壤，海拔與氧化鐵全量之間呈極顯著負(fù)相關(guān)關(guān)系(= –0.861 5，<0.001，= 26)。這與隨著海拔升高，溫度降低，相對(duì)濕度升高，土壤脫硅富鐵鋁化程度降低，氧化鐵全量隨之不斷減少有關(guān)。剖面土壤氧化鐵全量以剖面紅壤1最高，為73.01 ～ 88.56 g/kg，剖面黃棕壤2最低，僅為46.48 ～ 50.41 g/kg。所有剖面土壤氧化鐵全量均表現(xiàn)出AB層、B層和BC層較A層高的特征。這反映出研究區(qū)雨水豐富，淋溶作用強(qiáng)烈，氧化鐵具有從表土層向下聚集遷移的趨勢(shì)。隨著土壤剖面深度的增加，外界環(huán)境和淋溶作用減弱，土壤的氧化鐵全量逐步下降并趨于穩(wěn)定狀態(tài)，這與前人對(duì)太白山、廬山土壤氧化鐵全量的研究結(jié)論類似[26-27]。

圖1 各土壤剖面氧化鐵全量變化特征

圖2 土壤氧化鐵全量與海拔的相關(guān)性

2.2 土壤光譜曲線特征

根據(jù)樣本的土壤氧化鐵全量，將樣本分為 >70、60 ～ 70、50 ～ 60和 <50 g/kg 4個(gè)區(qū)間，分別計(jì)算樣本在 4 個(gè)區(qū)間內(nèi)的反射率平均值得到平均光譜曲線(圖3)。雖然樣品的氧化鐵全量不同，但高光譜平均反射率曲線均呈上凸的拋物線型。在430 ～ 560 nm，隨著氧化鐵全量的增加，反射率呈現(xiàn)下降的趨勢(shì)，氧化鐵和反射率為負(fù)相關(guān)關(guān)系；在500 ～ 620 nm波段，反射率曲線斜率增加最快，反射率以遞增的速率增加，形似陡坎；在620 ～ 780 nm 波段反射率曲線斜率略微變小，反射率以遞減的速率增加；在波長(zhǎng)780 nm處，反射率出現(xiàn)最大值，隨后波段反射率開始下降；在波長(zhǎng)900 nm附近，高光譜曲線存在一個(gè)微弱的吸收峰，被認(rèn)為是Fe2+、Fe3+的微弱吸收峰[28]。這與土壤氧化鐵全量增加時(shí)，可見光與近紅外部分吸收增強(qiáng)，使土壤光譜反射率朝著藍(lán)波段方向下降，導(dǎo)致土壤光譜反射率降低有關(guān)[7]。此結(jié)果與何挺等[11]、彭杰等[12]、郭穎等[29]認(rèn)為氧化鐵全量與反射光譜率為負(fù)相關(guān)一致。

圖3 不同氧化鐵全量的土壤反射率均值光譜

2.3 土壤氧化鐵全量的光譜特征波段

2.4 土壤氧化鐵全量多元逐步回歸模型

基于去包絡(luò)線全波段高光譜數(shù)據(jù)，建立多元線性回歸估測(cè)模型及誤差分析見表1。391、565、994和1040為對(duì)應(yīng)391、565、994及1 040 nm波段反射率經(jīng)去包絡(luò)線后的數(shù)值。通過(guò)比較2及 RPD值可以看出，模型 D 的決定系數(shù)(2)及相對(duì)分析誤差(RPD)達(dá)到最大，所以模型D建立的多元逐步回歸模型精度相對(duì)較高。

2.5 土壤氧化鐵全量的偏最小二乘回歸模型

將26個(gè)土壤樣本分成建模組與驗(yàn)證組的比例分別為0.5:0.5(樣本形式A)和0.75:0.25(樣本形式B)，基于去包絡(luò)線高光譜數(shù)據(jù)建立了2種偏最小二乘反演模型。由表2可知，2種樣本形式的 RPD 值均大于1.4，表明樣本形式A、B模型都能較好地估測(cè)土壤氧化鐵全量。雖然樣本形式A的建模組的決定系數(shù)R和校正均方根誤差 RMSEC與樣本形式B相接近，分別為0.930 6和3.643 4，但樣本形式A的驗(yàn)證組決定系數(shù)2僅為0.630 0，驗(yàn)證組預(yù)測(cè)均方根誤差RMSEP為7.740 0，表明該模型雖然建模效果好，但預(yù)測(cè)能力較差。樣本形式 B 的建模組的決定系數(shù)2達(dá)到0.935 8，校正均方根誤差RMSEC為3.419 7，驗(yàn)證組的2達(dá)到0.981 8，預(yù)測(cè)均方根誤差RMSEP為3.155 0，表明該模型不僅建模效果好，而且預(yù)測(cè)能力也很強(qiáng)?？梢?，樣本形式 B 的建模效果優(yōu)于A，且形式 B 的預(yù)測(cè)效果也更佳，能夠更好地估測(cè)土壤氧化鐵全量。

2.6 模型的驗(yàn)證與優(yōu)選

從多元逐步回歸模型和偏最小二乘回歸模型建模結(jié)果對(duì)比可見(表3)，偏最小二乘回歸模型樣本形式B精度更高，其驗(yàn)證組的決定系數(shù)(2)為0.981 8，預(yù)測(cè)均方根誤差RMSEP為3.155 0，預(yù)測(cè)相對(duì)分析誤差RPD為4.235 2。由兩種模型預(yù)測(cè)值和實(shí)測(cè)值擬合圖(圖5)也可看出，偏最小二乘回歸模型B的數(shù)據(jù)組相比多元逐步回歸模型，建模組和驗(yàn)證組的預(yù)測(cè)效果更好，該結(jié)論和表2的分析結(jié)果一致。

圖4 高光譜變換的相關(guān)性曲線

表 1 基于去包絡(luò)線的多元逐步回歸模型

表 2 基于去包絡(luò)線的偏最小二乘模型精度

表 3 建模結(jié)果對(duì)比

(A.多元逐步回歸模型預(yù)測(cè)值和實(shí)測(cè)值擬合；B.偏最小二乘回歸模型預(yù)測(cè)值和實(shí)測(cè)值擬合)

3 結(jié)論

1)不同森林土壤氧化鐵全量和不同土壤類型光譜曲線在可見光–近紅外波波段(400 ～ 1 000 nm)內(nèi)變化趨勢(shì)基本一致，均形似陡坎。在420 ～ 580 nm波段，氧化鐵全量與光譜反射率呈負(fù)相關(guān)關(guān)系。

2)土壤原始反射率經(jīng)過(guò)S-G平滑，并去包絡(luò)線和微分處理后，光譜曲線特征差異明顯，其中去包絡(luò)線的處理反演效果最好。土壤氧化鐵高光譜的特征波段主要為392、427、529、523、549、559、565、570、994和1 040 nm。

3)以去包絡(luò)線的土壤光譜反射率為自變量建立的土壤氧化鐵全量多元逐步回歸模型=–115.015×565–275.853×994+231.346×1040–76.7665×391+227.603精度較好，建模組與驗(yàn)證組的比例為0.75:0.25優(yōu)于比例為0.5:0.5的偏最小二乘回歸模型。其中，建模組與驗(yàn)證組的比例為0.75:0.25的偏最小二乘反演模型的2最高，為0.935 8；RMSEC最小，為3.419 7；RPD最高，為4.235 2，是較為適用于紅黃壤區(qū)森林土壤氧化鐵全量估測(cè)的高光譜模型。

[1] Schwertmann U, Cornell R M.Iron oxides in the laboratory: Preparation and characterization[M].Wiley‐VCH Verlag GmbH: 2000.

[2] 章明奎.土壤氧化鐵的形態(tài)在浙江紅壤和黃壤分類中的意義[J].浙江農(nóng)業(yè)大學(xué)學(xué)報(bào), 1990, 16(1): 42–50.

[3] Terra F S, Demattê J A M, Viscarra Rossel R A.Proximal spectral sensing in pedological assessments: Vis-NIR spectra for soil classification based on weathering and pedogenesis[J].Geoderma, 2018, 318: 123–136.

[4] Todorova M, Atanassova S, Sitaula B, et al.Application of NIRS as a rapid and alternative method for prediction of heavy metals content in soil [J].Agricultural Science and Technology, 2012, 4(4): 440-444.

[5] Gupta A, Vasava H B, Das B S, et al.Local modeling approaches for estimating soil properties in selected Indian soils using diffuse reflectance data over visible to near-infrared region[J].Geoderma, 2018, 325: 59–71.

[6] 史舟, 王乾龍, 彭杰, 等.中國(guó)主要土壤高光譜反射特性分類與有機(jī)質(zhì)光譜預(yù)測(cè)模型[J].中國(guó)科學(xué): 地球科學(xué), 2014, 44(5): 978–988.

[7] Fontes M P F, Carvalho I A Jr.Color attributes and mineralogical characteristics, evaluated by radiometry, of highly weathered tropical soils[J].Soil Science Society of America Journal, 2005, 69(4): 1162–1172.

[8] 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.

[9] Hunt G R, Ashley R P.Spectra of altered rocks in the visible and near infrared [J].Economic Geology, 1979, 74(7): 1613–1629.

[10] Galv?o L S, Pizarro M A, Epiphanio J C N.Variations in reflectance of tropical soils: Spectral-chemical composition relationships from AVIRIS data[J].Remote Sensing of Environment, 2001, 75(2): 245–255.

[11] 何挺, 王靜, 程燁, 等.土壤氧化鐵光譜特征研究[J].地理與地理信息科學(xué), 2006, 22(2): 30–34.

[12] 彭杰, 向紅英, 周清, 等.土壤氧化鐵的高光譜響應(yīng)研究[J].光譜學(xué)與光譜分析, 2013, 33(2): 502–506.

[13] de Aldana B R V, Criado B G, Ciudad A G, et al.Estimation of mineral content in natural grasslands by near infrared reflectance spectroscopy[J].Communications in Soil Science and Plant Analysis, 1995, 26(9/10): 1383–1396.

[14] 魏昌龍, 趙玉國(guó), 鄔登巍, 等.基于光譜分析的土壤游離鐵預(yù)測(cè)研究[J].土壤, 2014, 46(4): 678–683.

[15] Liu Y L, Chen Y Y.Estimation of total iron content in floodplain soils using VNIR spectroscopy - a case study in the Le'an River floodplain, China[J].International Journal of Remote Sensing, 2012, 33(18): 5954–5972.

[16] Magendran T, Sanjeevi S.Statistical and neural network analysis of hyperspectral radiometric data to characterise hematite of Singbhum iron ore belt, India[C]//SPIE Asia-Pacific Remote Sensing.Proc SPIE 9263, Multispectral, Hyperspectral, and Ultraspectral Remote Sensing Technology, Techniques and Applications V, Beijing, China.2014, 9263: 92630M.

[17] 秦倩如, 齊雁冰, 吳娟, 等.基于高光譜的土壤游離鐵隨機(jī)森林模型估算研究[J].土壤通報(bào), 2018, 49(6): 1286–1293.

[18] 謝文, 趙小敏, 郭熙, 等.基于組合模型的廬山森林土壤有效鐵光譜間接反演研究[J].土壤學(xué)報(bào), 2017, 54(3): 601–612.

[19] 高開強(qiáng), 徐泮林, 趙曉旭.基于模糊模式識(shí)別的土壤含鐵量高光譜估測(cè)[J].測(cè)繪與空間地理信息, 2020, 43(1): 127–129,132.

[20] 趙振華, 劉偉, 趙亞輝, 等.大圍山第四紀(jì)冰川地質(zhì)遺跡特征及成因探討[J].國(guó)土資源導(dǎo)刊, 2008, 5(6): 42–45.

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

[22] 王濤, 喻彩麗, 張楠楠, 等.基于去包絡(luò)線和連續(xù)投影算法的棗園土壤電導(dǎo)率光譜檢測(cè)研究[J].干旱地區(qū)農(nóng)業(yè)研究, 2019, 37(5): 193–199,217.

[23] 魏建榕.經(jīng)濟(jì)物理學(xué)中的金融數(shù)據(jù)分析：統(tǒng)計(jì)與建模[D].上海: 復(fù)旦大學(xué), 2013.

[24] Rossel R A V, Webster R.Predicting soil properties from the Australian soil visible-near infrared spectroscopic database [J].European Journal of Soil Science, 2012, 63(6): 848–860.

[25] 羅蘭花, 王翠紅, 謝紅霞, 等.大圍山花崗巖風(fēng)化物發(fā)育土壤抗蝕性垂直分異[J].水土保持研究, 2018, 25(1): 62–65, 71.

[26] 雷梅, 常慶瑞, 馮立孝, 等.太白山土壤特性及氧化鐵發(fā)生學(xué)特征[J].地理研究, 2001, 20(1): 83–90.

[27] 曹慶.廬山北坡土壤發(fā)生特性與系統(tǒng)分類[D].南京: 南京師范大學(xué), 2012.

[28] 龔紹琦, 王鑫, 沈潤(rùn)平, 等.濱海鹽土重金屬含量高光譜遙感研究[J].遙感技術(shù)與應(yīng)用, 2010, 25(2): 169–177.

[29] 郭穎, 郭治興, 劉佳, 等.亞熱帶典型區(qū)域水稻土氧化鐵高光譜反演——以珠江三角洲為例[J].應(yīng)用生態(tài)學(xué)報(bào), 2017, 28(11): 3675–3683.

[30] Baumgardner M F, Silva L F, Biehl L L, et al.Reflectance properties of soils [J].Advances in Agronomy, 1986, 38: 1–44.

Inversion of Iron Oxide Contents in Forest Soils of Dawei Mountains Using Laboratory Hyperspectral Data

TAN Jie1, CHEN Yan1, ZHOU Weijun2*, CUI Haojie2, LIU Pei1

(1 College of Landscape Architecture and Art Design, Hunan Agricultural University, Changsha 410128, China; 2College of Resources and Environment, Hunan Agricultural University, Changsha 410128, China)

Iron oxide is the main body of iron-bearing minerals in the soil and is one of the most obvious and useful indicators of soil development and soil classification.In this paper, forest soils in Dawei Mountain of Hunan Province were collected, iron oxide contents in soils were determined respectively by conventional chemical method and by hyperspectral inversion with models of multiple stepwise regression and partial least squares regression inversion established by screening sensitive bands after spectral preprocessing and combinatorial transformation.The results showed that soil spectral curves with different iron oxide contents all were in steep-hill shape in the whole band, iron oxide content was negatively correlated with spectral reflectance within 420-580 nm, different spectral transformation could improve the correlation, and the combination of Savitzky-Golay (S-G) smoothing and de-embedding lines was superior to other pretreatment methods in inversion.The characteristic bands of iron oxides were 392, 427, 529, 523, 549, 559, 565, 570, 994 and 1 040 nm.Partial least squares regression model had better stability than multiple stepwise regression model, and is suitable for rapid estimation of iron oxide contents in forest red and yellow soils.

Soil spectra; Iron oxide; Multiple stepwise regression; Partial least squares regression

S155.1

A

10.13758/j.cnki.tr.2021.04.025

譚潔, 陳嚴(yán), 周衛(wèi)軍, 等.基于實(shí)驗(yàn)室高光譜數(shù)據(jù)的大圍山森林土壤氧化鐵全量反演.土壤, 2021, 53(4): 858–864.

國(guó)家自然科學(xué)基金項(xiàng)目(41771272)資助。

(wjzh0108@163.com)

譚潔(1979—)，女，湖南湘潭人，博士，副教授，主要從事土地/土壤環(huán)境過(guò)程及模擬研究。E-mail: 37080689@qq.com