基于深度稀疏學(xué)習(xí)的土壤近紅外光譜分析預(yù)測(cè)模型

王儒敬, 陳天嬌, 汪玉冰, 汪六三, 謝成軍, 張 潔, 李 瑞, 陳紅波

(中國(guó)科學(xué)院 合肥智能機(jī)械研究所, 安徽 合肥 230031)

?

基于深度稀疏學(xué)習(xí)的土壤近紅外光譜分析預(yù)測(cè)模型

王儒敬, 陳天嬌*, 汪玉冰, 汪六三, 謝成軍, 張 潔, 李 瑞, 陳紅波

(中國(guó)科學(xué)院 合肥智能機(jī)械研究所, 安徽 合肥 230031)

提出一種基于深度稀疏學(xué)習(xí)的土壤近紅外光譜分析預(yù)測(cè)模型。首先，使用稀疏特征學(xué)習(xí)方法對(duì)土壤近紅外光譜數(shù)據(jù)進(jìn)行約簡(jiǎn)，實(shí)現(xiàn)土壤近紅外光譜內(nèi)容的稀疏表示；然后采用徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)以稀疏表示特征系數(shù)為輸入，以所測(cè)土壤成分為輸出，分別建立土壤有機(jī)質(zhì)、速效磷、速效鉀的非線性預(yù)測(cè)模型。結(jié)果表明用該模型預(yù)測(cè)土壤有機(jī)質(zhì)的含量是可行的，但對(duì)土壤速效磷和速效鉀含量的預(yù)測(cè)還需對(duì)模型做進(jìn)一步的優(yōu)化。

土壤近紅外光譜； 深度稀疏學(xué)習(xí); 神經(jīng)網(wǎng)絡(luò)模型

1 引 言

近年來(lái)，國(guó)家為了能準(zhǔn)確掌握土壤養(yǎng)分分布信息,每年都要投入幾十甚至上百億元的專項(xiàng)資金 在全國(guó)范圍內(nèi)實(shí)施推廣測(cè)土配方施肥工作，采用網(wǎng)格式定點(diǎn)取樣和實(shí)驗(yàn)室分析檢測(cè)相結(jié)合的方法。然而，傳統(tǒng)的土壤檢測(cè)方法存在著一系列的弊端，不僅需要投入大量的人力物力，檢測(cè)時(shí)間長(zhǎng)，成本高，而且所使用的化學(xué)藥品對(duì)環(huán)境也具有一定的污染性[1]。近紅外光譜分析方法具有操作簡(jiǎn)便、成本低廉、耗時(shí)短、綠色無(wú)污染等優(yōu)點(diǎn)[2]，近年來(lái)，可見(jiàn)/近紅外光譜已被國(guó)內(nèi)外學(xué)者廣泛應(yīng)用于土壤多種養(yǎng)分含量的檢測(cè)中。近紅外光譜分析方法是一種間接校正方法，分析的準(zhǔn)確性和魯棒性取決于所建預(yù)測(cè)模型的質(zhì)量[3]。目前，人們主要使用線性和非線性的多元校正方法進(jìn)行可見(jiàn)/近紅外光譜模型的搭建。在土壤近紅外光譜分析中，建立定量分析模型的線性方法主要有多元線性回歸(MLR)、主成分回歸(PCR)、偏最小二乘(PLS)以及局部權(quán)重的偏最小二乘回歸(LWPLS)等，其中 PLS 是目前主流的線性定量分析建模方法。何勇等[4-5]采用NIRS結(jié)合偏最小二乘回歸方法對(duì)土壤中有機(jī)質(zhì)和三大營(yíng)養(yǎng)元素的含量進(jìn)行了回歸建模分析。朱登勝等[6]用偏最小二乘回歸方法對(duì)經(jīng)簡(jiǎn)單處理土樣的有機(jī)質(zhì)含量和pH 值進(jìn)行了研究。李偉等[7]采用偏最小二乘法和人工神經(jīng)網(wǎng)絡(luò)方法分別建立了0.9 mm 篩的風(fēng)干土壤堿解氮、速效磷和速效鉀含量預(yù)測(cè)的近紅外光譜分析模型, 結(jié)果表明，人工神經(jīng)網(wǎng)絡(luò)方法所建立的校正模型優(yōu)于偏最小二乘法所建立的模型。Mouazen等[8]結(jié)合偏最小二乘回歸和BP神經(jīng)網(wǎng)絡(luò)建立了預(yù)測(cè)模型，其效果優(yōu)于偏最小二乘和主成分回歸模型。這些線性方法在解決非線性的土壤各項(xiàng)指標(biāo)預(yù)測(cè)問(wèn)題上均具有一定的局限性。近年來(lái)，學(xué)者們又開(kāi)始研究將一些非線性回歸或者分類預(yù)測(cè)方法與土壤近紅外光譜數(shù)據(jù)結(jié)合起來(lái)進(jìn)行建模，以充分發(fā)掘待測(cè)樣品組分與采集的光譜數(shù)據(jù)之間的非線性關(guān)系，如神經(jīng)網(wǎng)絡(luò)(BP算法)、支持向量回歸等，這些方法在解決非線性問(wèn)題與模型抗干擾方面比PLS等優(yōu)越[9-10]。神經(jīng)網(wǎng)絡(luò)常用的方法有BP[11-15]、自組織映射[16-19]和徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)[20-22]。

深度學(xué)習(xí)本質(zhì)是復(fù)雜的非線性模型的學(xué)習(xí)。從機(jī)器學(xué)習(xí)的發(fā)展史來(lái)看，深度學(xué)習(xí)的興起代表著機(jī)器學(xué)習(xí)技術(shù)的自然演進(jìn)，自2006年被Hinton等[23]提出以來(lái)，已被成功地應(yīng)用到圖像處理、語(yǔ)音處理、自然語(yǔ)言處理等多個(gè)領(lǐng)域，取得了巨大成功，成為當(dāng)今具有代表性的IT先進(jìn)技術(shù)[24-25]。光譜數(shù)據(jù)量大且冗余，如何有效地表示光譜信號(hào)是近紅外土壤預(yù)測(cè)處理的基礎(chǔ)，而高效的表示又要求用較少的信息(如數(shù)據(jù)壓縮、特征抽取等)來(lái)描述光譜的重要特征。由于土壤近紅外光譜信息獲取的時(shí)間連續(xù)性與空間關(guān)聯(lián)性，土壤樣本近紅外光譜數(shù)據(jù)包含有大量冗余，有必要去除土壤近紅外光譜冗余信息，實(shí)現(xiàn)土壤近紅外光譜內(nèi)容的非線性表示，提高光譜表示能力。作為深度學(xué)習(xí)模型之一的稀疏表示學(xué)習(xí)是一種無(wú)監(jiān)督學(xué)習(xí)方法，它用來(lái)尋找一組“超完備”基向量來(lái)更高效地表示樣本數(shù)據(jù)。數(shù)據(jù)的稀疏表示可以降低數(shù)據(jù)處理的成本，提高壓縮效率[26-27]。對(duì)光譜信號(hào)的稀疏表示涉及冗余字典的選取，MOD、MAD、FOCUSS、K-SVD[28]等字典訓(xùn)練算法已具備利用近紅外光譜完成冗余字典的訓(xùn)練能力。

本文對(duì)室內(nèi)光譜儀測(cè)量的土壤近紅外反射率光譜進(jìn)行研究,嘗試使用土壤近紅外光譜在訓(xùn)練字典深度稀疏特征學(xué)習(xí)下結(jié)合徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)的回歸校正方法來(lái)預(yù)測(cè)土壤的有機(jī)質(zhì)、速效磷和速效鉀含量，研究了土壤近紅外光譜在冗余字典的稀疏分解下結(jié)合徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)的回歸預(yù)測(cè)方法能否顯著提高土壤重要成分含量的預(yù)測(cè)性能。

2 實(shí) 驗(yàn)

2.1土壤樣本的采集

供試土樣來(lái)自安徽省亳州市蒙城縣某農(nóng)場(chǎng)的田塊。采樣時(shí)分別選擇某一小區(qū)，在小區(qū)4個(gè)角落和其對(duì)角線的交點(diǎn)處分別采樣，采樣深度為30 cm，然后將這些樣本混合作為一個(gè)土樣。總共采集了420個(gè)土樣，隨機(jī)選取400個(gè)樣本作為校正集，其余20個(gè)樣本作為驗(yàn)證集，將樣本分為 2 份以備后續(xù)試驗(yàn)使用，其中一份用于化學(xué)分析，另一份用于近紅外分析。

2.2土壤養(yǎng)分的測(cè)量

分別采用常規(guī)化學(xué)分析方法獲得校正集土壤樣本中有機(jī)質(zhì)( OM) 、速效磷( P) 、速效鉀( K) 的測(cè)量結(jié)果，統(tǒng)計(jì)參數(shù)見(jiàn)表1。

表1 土壤化學(xué)成分統(tǒng)計(jì)參數(shù)

Tab.1 Statistical data of soil properties measured by chemical methods

OM/(g·kg-1)速效P/(mg·kg-1)速效K/(mg·kg-1)最小值10.46.3692最大值29.548.24256平均值18.728.35163標(biāo)準(zhǔn)偏差值1.172.786.8

3 土壤近紅外光譜的稀疏表示

土壤近紅外信息的時(shí)間和地理空間位置強(qiáng)相關(guān)，因而光譜結(jié)構(gòu)復(fù)雜，數(shù)據(jù)量大，且包含有大量冗余。稀疏表示可以降低數(shù)據(jù)處理的成本，提高壓縮效率，實(shí)現(xiàn)光譜降噪、壓縮和有效的特征表示。另外，由于基于神經(jīng)網(wǎng)絡(luò)的非線性預(yù)測(cè)模型計(jì)算量很大，所以很有必要對(duì)光譜進(jìn)行一定程度的稀疏降維。

在土壤近紅外光譜信號(hào)的有限維空間Rn中, 近紅外光譜y∈Rn可由一組稀疏基D=[d1d1...dk],dk∈Rn線性表示，y=Dx，如圖1所示。矩陣D中列向量代表稀疏基, 若k>n,則認(rèn)為D是冗余的(過(guò)完備、超完備),稱為冗余字典。當(dāng)訓(xùn)練樣本數(shù)量較大時(shí)，此時(shí)的稀疏編碼非常耗時(shí)，因此學(xué)習(xí)一個(gè)更加緊湊的并且具有稀疏表示能力的字典十分必要。

圖1 稀疏表示

從線性組合角度看，可以將上述問(wèn)題轉(zhuǎn)換為式(1)所示的優(yōu)化問(wèn)題：

(1)

其中T0為稀疏表示系數(shù)中非零分量的數(shù)目的上限。對(duì)已知的光譜信號(hào)做預(yù)處理后(預(yù)處理方法見(jiàn)5.1)，應(yīng)用非監(jiān)督學(xué)習(xí)算法奇異值分解算法(K-SVD)創(chuàng)建過(guò)完備字典，求解如下優(yōu)化問(wèn)題：

(2)

Y=[y1,y2...yN]為已知樣本的光譜信號(hào)集合，xi是某一個(gè)樣本的稀疏表示向量，X=[x1,x2...xN]，因此通過(guò)解式(2)最優(yōu)化問(wèn)題，可以得出土壤近紅外光譜的過(guò)完備訓(xùn)練字典D。

4 深度稀疏學(xué)習(xí)下的土壤近紅外光譜預(yù)測(cè)模型

由于近紅外光譜儀的狀態(tài)、測(cè)量環(huán)境等因素對(duì)光譜的影響大多屬于非線性的，還有一些質(zhì)量參數(shù)和光譜的關(guān)系也是非線性的，所以最直接的解決方法是非線性校正方法。神經(jīng)網(wǎng)絡(luò)在處理高復(fù)雜度的非線性問(wèn)題時(shí)有顯著的優(yōu)越性，常用的方法有BP神經(jīng)網(wǎng)絡(luò)和徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)等。在逼近能力、分類能力和學(xué)習(xí)速率等方面，徑向基函數(shù) (Radial basis function，RBF)網(wǎng)絡(luò)表現(xiàn)出比 BP網(wǎng)絡(luò)更優(yōu)的性能，其結(jié)構(gòu)如圖2所示。

圖2 RBF神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

隱含層是基于高斯徑向基函數(shù)，其形式為:

(3)

輸入層實(shí)現(xiàn)從光譜特征向量x→Ri(x)的非線性映射，輸出層實(shí)現(xiàn)從Ri(x)→y的線性映射，即：

(4)

本文預(yù)測(cè)模型如圖3所示，主要步驟如下：

(1)根據(jù)已知的光譜數(shù)據(jù)集應(yīng)用非監(jiān)督學(xué)習(xí)算法K-SVD創(chuàng)建過(guò)完備字典；

(2)求解校正集在字典上的稀疏系數(shù)；

(3)將校正集的稀疏系數(shù)作為輸入層參數(shù)，與該光譜對(duì)應(yīng)的樣本的成分含量參數(shù)作為輸出層參數(shù)，通過(guò)徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)訓(xùn)練RBF神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型；

(4)求解驗(yàn)證集在字典上的稀疏系數(shù)，輸入訓(xùn)練好的RBF網(wǎng)絡(luò),可得對(duì)驗(yàn)證集樣本某一成分含量的預(yù)測(cè)值。

5 實(shí)驗(yàn)結(jié)果及其分析

5.1光譜數(shù)據(jù)的采集

將采集的土壤樣本在室內(nèi)攤開(kāi)，風(fēng)干、碾碎后，將土壤樣本過(guò)1 mm的篩，然后用于光譜實(shí)驗(yàn)。

圖3 光譜冗余字典稀疏表示下的RBF神經(jīng)網(wǎng)絡(luò)回歸預(yù)測(cè)

實(shí)驗(yàn)中，采用鹵素?zé)粽丈渫寥辣砻妫寥赖穆瓷涔庾V通過(guò)光纖被可見(jiàn)光/近紅外光譜儀(美國(guó)Ocean Optics公司的USB4000光譜儀和日本Hamamstu公司的TG-cooled NIR-2 光譜儀)接收，USB4000光譜儀的波長(zhǎng)范圍為350～1 000 nm，分辨率為1.5～2.3 nm；TG-cooled NIR-2光譜儀的波長(zhǎng)范圍為1 100～2 200 nm，最大光學(xué)分辨率為8 nm。

測(cè)量前根據(jù)儀器測(cè)量要求，調(diào)整好測(cè)量高度和角度，將被測(cè)土樣放在培養(yǎng)皿，每一位置的測(cè)量次數(shù)為5次，每個(gè)樣本的掃描次數(shù)為15次。平均處理后，將其按照l(shuí)g(1/R)轉(zhuǎn)換為吸光度，然后將數(shù)據(jù)導(dǎo)出進(jìn)行譜圖預(yù)處理。使用目前應(yīng)用較廣泛的Savitzky-Golay卷積平滑法和多元散射校正進(jìn)一步對(duì)光譜數(shù)據(jù)進(jìn)行預(yù)處理。一個(gè)土壤樣本的光譜數(shù)據(jù)共有381個(gè)點(diǎn)，圖4為校正集部分土壤樣本的近紅外光譜圖。圖中橫坐標(biāo)是波長(zhǎng)(nm), 縱坐標(biāo)為吸光度lg(1/R)。訓(xùn)練時(shí)每個(gè)樣本去除原始光譜邊緣不平滑區(qū)域，只取波長(zhǎng) 600～2 000 nm的光譜信息使用。

圖4 土壤的紅外光譜

5.2光譜稀疏度對(duì)模型預(yù)測(cè)的影響

由于光譜數(shù)據(jù)中噪聲的不確定性, 實(shí)驗(yàn)設(shè)置T0參數(shù)在6～15范圍內(nèi)分別對(duì)有機(jī)質(zhì)、速效K和速效P模型預(yù)測(cè)結(jié)果進(jìn)行觀察。對(duì)校正集采取交叉驗(yàn)證, 可利用交叉驗(yàn)證標(biāo)準(zhǔn)差(SECV)來(lái)確定模型最佳預(yù)測(cè)結(jié)果下的稀疏度T0，圖5表示了土壤成分有機(jī)質(zhì)、速效P和速效K在不同稀疏度下的模型預(yù)測(cè)結(jié)果。當(dāng)稀疏度為10時(shí), 有機(jī)質(zhì)、速效P和速效K模型的SECV基本均具有局部極小點(diǎn)。

5.3不同字典表示下的土壤成分預(yù)測(cè)

將校正集樣本由小波字典、傅立葉字典以及訓(xùn)練字典下經(jīng)過(guò)稀疏表示提取的稀疏系數(shù)作為RBF神經(jīng)網(wǎng)絡(luò)的輸入變量分別建立預(yù)測(cè)模型, 再通過(guò)調(diào)整隱含層的節(jié)點(diǎn)數(shù)來(lái)優(yōu)化網(wǎng)絡(luò)結(jié)構(gòu)。針對(duì)不同的土壤成分，經(jīng)多次實(shí)驗(yàn)后, 得到土壤不同成分網(wǎng)絡(luò)模型的最佳隱含層節(jié)點(diǎn)數(shù),設(shè)定目標(biāo)誤差，其他參數(shù)設(shè)置為默認(rèn)值。訓(xùn)練字典結(jié)合RBF網(wǎng)絡(luò)模型的參數(shù)值如表2所示。用3種不同的字典結(jié)合RBF神經(jīng)網(wǎng)絡(luò)建模方法建立土壤有機(jī)質(zhì)、速效磷和速效鉀的含量預(yù)測(cè)模型，對(duì)驗(yàn)證集的20個(gè)樣本進(jìn)行預(yù)測(cè)，實(shí)驗(yàn)結(jié)果如表3～5所示。

圖5 不同稀疏度下的模型預(yù)測(cè)結(jié)果。(a)有機(jī)質(zhì)；(b)速效磷；(c)速效鉀。

Fig.5 Prediction results under different sparsity. (a) OM. (b) Soil available P. (c) Soil available K.

表2 RBF網(wǎng)絡(luò)模型參數(shù)值

表3 土壤有機(jī)質(zhì)在不同字典表示下的建模方法的結(jié)果

Tab.3 Results of different modeling methods for soil available OM

字典R2RMSEP小波字典0.76940.3582傅里葉字典0.79210.2783訓(xùn)練字典0.87370.1420

表4 土壤速效磷在不同字典表示下的建模方法的結(jié)果

Tab.4 Results of different modeling methods for soil available P

字典R2RMSEP小波字典0.521810.8283傅里葉字典0.56829.8429訓(xùn)練字典0.68395.7320

表5 土壤速效鉀在不同字典表示下的建模方法的結(jié)果

Tab.5 Results of different modeling methods for soil available K

字典R2RMSEP小波字典0.581014.6385傅立葉字典0.538913.9453訓(xùn)練字典0.66729.4376

如表3～5所示, 有機(jī)質(zhì)、速效磷和速效鉀的基于訓(xùn)練字典稀疏表示下結(jié)合RBF神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)模型的相關(guān)系數(shù)(R2：0.873 7,0.683 9，0.667 2)最大，預(yù)測(cè)集標(biāo)準(zhǔn)差(RMSEP：0.142 0, 5.732 0，9.437 6)最小,表明訓(xùn)練字典相比于現(xiàn)有的小波字典和傅立葉字典對(duì)土壤近紅外光譜具有更好的表示能力。

5.4不同網(wǎng)絡(luò)模型下的土壤成分預(yù)測(cè)

將校正集樣本由訓(xùn)練字典下經(jīng)過(guò)稀疏表示提取的稀疏系數(shù)作為BP神經(jīng)網(wǎng)絡(luò)和RBF神經(jīng)網(wǎng)絡(luò)的輸入變量分別建立預(yù)測(cè)模型，再通過(guò)調(diào)整隱含層的節(jié)點(diǎn)數(shù)來(lái)優(yōu)化網(wǎng)絡(luò)結(jié)構(gòu)。用訓(xùn)練字典結(jié)合RBF神經(jīng)網(wǎng)絡(luò)建模方法建立土壤有機(jī)質(zhì)、速效磷和速效鉀的含量預(yù)測(cè)模型，對(duì)驗(yàn)證集樣本進(jìn)行預(yù)測(cè)，實(shí)驗(yàn)結(jié)果如表6～8所示。

表6 土壤有機(jī)質(zhì)在不同網(wǎng)絡(luò)模型下的結(jié)果對(duì)比

Tab.6 Results of different modeling methods for soil available OM

網(wǎng)絡(luò)模型R2RMSEPBP神經(jīng)網(wǎng)絡(luò)0.84290.2728RBF神經(jīng)網(wǎng)絡(luò)0.87370.1420

表7 土壤速效磷在不同網(wǎng)絡(luò)模型下的結(jié)果對(duì)比

Tab.7 Results of different modeling methods for soil available P

網(wǎng)絡(luò)模型R2RMSEPBP神經(jīng)網(wǎng)絡(luò)0.65729.1029RBF神經(jīng)網(wǎng)絡(luò)0.68395.7320

如表6～8所示, 有機(jī)質(zhì)、速效磷和速效鉀的基于訓(xùn)練字典稀疏表示下結(jié)合RBF神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)模型的相關(guān)系數(shù)(R2：0.873 7, 0.683 9，0.6672)最大(圖6)，預(yù)測(cè)集標(biāo)準(zhǔn)差(RMSEP：0.1420, 5.7320，9.4376)最小,表明RBF網(wǎng)絡(luò)對(duì)土壤近紅外光譜具有更好的預(yù)測(cè)能力。同時(shí)，有機(jī)質(zhì)的預(yù)測(cè)模型相關(guān)度為0.873 7，預(yù)測(cè)集標(biāo)準(zhǔn)差為0.142 0, 表明該方法預(yù)測(cè)土壤有機(jī)質(zhì)含量是可行的。但對(duì)土壤速效磷和速效鉀含量的預(yù)測(cè)并不是特別理想, 還有待進(jìn)一步研究。

圖6 RBF模型預(yù)測(cè)值和測(cè)定值之間的關(guān)系。(a)有機(jī)質(zhì)；(b)速效磷；(c)速效鉀。

Fig.6 Correlativity between the predicted data of the RBF model and measured data. (a) OM. (b) Soil available P. (c) Soil available K.

表8 土壤速效鉀在不同網(wǎng)絡(luò)模型下的結(jié)果對(duì)比

Tab.8 Results of different modeling methods for soil available K

網(wǎng)絡(luò)模型R2RMSEPBP神經(jīng)網(wǎng)絡(luò)0.531012.3241RBF神經(jīng)網(wǎng)絡(luò)0.66729.4376

6 結(jié) 論

本文基于深度稀疏學(xué)習(xí)特征和神經(jīng)網(wǎng)絡(luò)方法，利用近紅外光譜數(shù)據(jù)分別建立了土壤有機(jī)質(zhì)、速效磷、速效鉀預(yù)測(cè)模型,并對(duì)其預(yù)測(cè)結(jié)果進(jìn)行驗(yàn)證。結(jié)果表明：對(duì)土壤有機(jī)質(zhì)、速效磷和速效鉀驗(yàn)證集樣本的預(yù)測(cè)相關(guān)度分別為0.873 7、0.683 9和0.6672，均方根誤差分別為0.142 0、5.732 0和9.437 6，相對(duì)誤差絕對(duì)平均值分別為2.18%、4.09 %和5.38%。這說(shuō)明采用本文方法對(duì)土壤有機(jī)質(zhì)進(jìn)行預(yù)測(cè)是可行的，針對(duì)土壤速效鉀和速效磷含量的預(yù)測(cè)需對(duì)該模型進(jìn)一步優(yōu)化和研究。隨著土壤近紅外光譜檢測(cè)研究的深入，光譜數(shù)據(jù)勢(shì)必呈幾何級(jí)數(shù)增長(zhǎng)。在分析理解這些海量數(shù)據(jù)時(shí)，該模型的健壯性受限，需進(jìn)一步研究可以更為深層次刻畫(huà)近紅外吸收光譜參數(shù)與樣品含量化學(xué)測(cè)定值之間的非線性關(guān)系預(yù)測(cè)模型。

[1] 隋方功, 李俊良. 土壤農(nóng)化分析實(shí)驗(yàn) [M]. 北京: 中國(guó)農(nóng)業(yè)大學(xué)出版社, 2004. SUI F G, LI J L.AgrochemicalSoilAnalysisExperiment[M]. Beijing: China Agriculture Press, 2004. (in Chinese)

[2] 陸婉珍. 現(xiàn)代近紅外光譜分析技術(shù) [M]. 2版. 北京: 中國(guó)石化出版社, 2007. LU W Z.ModernNearInfraredSpectroscopyAnalyticatTechnology[M]. 2nd ed. Beijing: China Petrochemical Press, 2007. (in Chinese)

[3] 褚小立, 陸婉珍. 近五年我國(guó)近紅外光譜分析技術(shù)研究與應(yīng)用進(jìn)展 [J]. 光譜學(xué)與光譜分析, 2014, 34(10): 2595-2605. CHU X L, LU W Z. Research and application progress of near infrared spectroscopy analytical technology in China in the past five years [J].Spectrosc.Spect.Anal., 2014, 34(10):2595-2605. (in Chinese)

[4] HE Y, HUANG M, GARCA A,etal.. Prediction of soil macronutrients content using near-infrared spectroscopy[J].Comput.Electron.Agric., 2007, 58(2):144-153.

[5] 鮑一丹, 何勇, 方慧, 等. 土壤的光譜特征及氮含量的預(yù)測(cè)研究 [J]. 光譜學(xué)與光譜分析, 2001, 27(1):62-65. BAO Y D, HE Y, FANG H,etal.. Spectral characterization and N content prediction of soil with different particle size and moisture content [J].Spectrosc.Spect.Anal., 2001, 27(1):62-65. (in Chinese)

[6] 朱登勝, 吳迪, 宋海燕, 等. 應(yīng)用近紅外光譜法測(cè)定土壤的有機(jī)質(zhì)和pH值 [J]. 農(nóng)業(yè)工程學(xué)報(bào), 2008, 24(6): 196-199. ZHU D S, WU D, SONG H Y,etal.. Determination of organic matter contents and pH values of soil using near infrared spectroscopy [J].Trans.CSAE, 2008, 24(6):196-199. (in Chinese)

[7] 李偉, 張書(shū)慧, 張倩, 等. 近紅外光譜法快速測(cè)定土壤堿解氮、速效磷和速效鉀含量 [J]. 農(nóng)業(yè)工程學(xué)報(bào), 2007, 23(1):55-59. LI W, ZHANG S H, ZHANG Q,etal.. Rapid prediction of available N, P and K content in soil using near-infrared reflectance spectroscopy[J].Trans.CSAE, 2007, 23(1):55-59. (in Chinese)

[8] MOUAZEN A M, KUANG B, DE BAERDEMAEKER J,etal.. Comparison among principal component, partial least squares and back propagation neural network analyses for accuracy of measurement of selected soil properties with visible and near infrared spectroscopy [J].Geoderma, 2010, 158(1-2):23-31.

[9] 許祿, 邵學(xué)廣. 化學(xué)計(jì)量學(xué)方法 [M]. 2版. 北京: 科學(xué)出版社, 2004. XU L, SHAO X G.ChemometricMethods[M]. 2nd ed. Beijing: Science Press, 2004. (in Chinese)

[10] 褚小立. 化學(xué)計(jì)量學(xué)方法與分子光譜分析技術(shù) [M]. 北京: 化學(xué)工業(yè)出版社, 2011. CHU X L.MolecularSpectroscopyAnalyticalTechnologyCombinedwithChemometricsanditsApplications[M]. Beijing: Chemical Industry Press, 2011. (in Chinese)

[11] 蔣宗禮. 人工神經(jīng)網(wǎng)絡(luò)導(dǎo)論 [M]. 北京: 高等教育出版社, 2001:8. JIANG Z L.IntroductiontoAritificialNeuralNetworks[M]. Beijing: Higher Education Press, 2001:8. (in Chinese)

[12] HAYKIN S. 神經(jīng)網(wǎng)絡(luò)與機(jī)器學(xué)習(xí) [M]. 申富饒, 徐燁, 鄭俊, 等, 譯. 3版. 北京: 機(jī)械工業(yè)出版社, 2011:3. HAYKIN S.NeuralNetworksandLearningMachines[M]. SHEN F R, XU Y, ZHENG J,etal.. Trans. 3rd ed. Beijing: Mechanical Industry Press, 2011:3. (in Chinese)

[13] 牛曉穎, 邵利敏, 趙志磊, 等. 基于BP-ANN的草莓品種近紅外光譜無(wú)損鑒別方法研究 [J]. 光譜學(xué)與光譜分析, 2012, 32(8):2095-2099. NIU X Y, SHAO L M, ZHAO Z L,etal.. Nondestructive discrimination of strawberry varieties by NIR and BP-ANN [J].Spectrosc.Spect.Anal., 2012, 32(8):2095-2099. (in Chinese)

[14] QU Z H, WANG L H. Prediction of lignin content of manchurian walnut by BP neural network and near-infrared spectroscopy [J].Adv.Mater.Res., 2011, 267:991-994.

[15] 鄭立華, 李民贊, 潘孌, 等. 基于近紅外光譜技術(shù)的土壤參數(shù)BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè) [J]. 光譜學(xué)與光譜分析, 2008, 28(5):1160-1164. ZHENG L H, LI M Z, PAN L,etal.. Estimation of soil organic matter and soil total nitrogen based on NIR spectroscopy and BP neural network [J].Spectrosc.Spect.Anal., 2008, 28(5):1160-1164. (in Chinese)

[16] VESANTO J, ALHONIEMI E. Clustering of the self-organizing map [J].IEEETrans.NeuralNetw., 2000, 11(3): 586-600.

[17] KOHONEN T. Essentials of the self-organizing map [J].NeuralNetw., 2013, 37:52-65.

[18] CRéPUT J C, HAJJAM A, KOUKAM A,etal.. Self-organizing maps in population based metaheuristic to the dynamic vehicle routing problem [J].J.Comb.Optim., 2012, 24(4):437-458.

[19] BALLABIO D, VASIGHI M. A MATLAB toolbox for self organizing maps and supervised neural network learning strategies [J].Chemom.Intell.Lab.Syst., 2012, 118:24-32.

[20] HUANG G B, SARATCHANDRAN P, SUNDARARAJAN N. A generalized growing and pruning RBF (GGAP-RBF) neural network for function approximation [J].IEEETrans.NeuralNetw., 2005, 16(1):57-67.

[21] MELAGRAKI G, AFANTITIS A, MAKRIDIMA K,etal.. Prediction of toxicity using a novel RBF neural network training methodology[J].J.Mol.Model., 2006, 12(3):297-305.

[22] MELAGRAKI G, AFANTITIS A, MAKRIDIMA K,etal.. Prediction of toxicity using a novel RBF neural network training methodology [J].J.Mol.Model., 2006, 12(3):297-305.

[23] HINTON G E, SALAKHUTDINOV R R. Reducing the dimensionality of data with neural networks[J].Science, 2006, 313(5786):504-507.

[24] 余凱, 賈磊, 陳雨強(qiáng), 等. 深度學(xué)習(xí)的昨天、今天和明天 [J]. 計(jì)算機(jī)研究與發(fā)展, 2013, 50(9):1799-1804. YU K, JIA L, CHEN Y Q,etal.. Deep learning: yesterday, today, and tomorrow [J].J.Comput.Res.Dev., 2013, 50(9):1799-1804. (in Chinese)

[25] 劉建偉, 劉媛, 羅雄麟. 深度學(xué)習(xí)研究進(jìn)展 [J]. 計(jì)算機(jī)應(yīng)用研究, 2014, 31(7):1921-1930. LIU J W, LIU Y, LUO X L. Research and development on deep learning [J].Appl.Comput.Res., 2014, 31(7):1921-1930. (in Chinese)

[26] DONOHO D L. Compressed sensing [J].IEEETrans.Inf.Theory, 2006, 52(4):1289-1306.

[27] ZHANG C M, YIN Z K, CHEN X D,etal.. Signal overcomplete representation and sparse decomposition based on redundant dictionaries [J].Chin.Sci.Bull., 2005, 50(23):2672-2677.

[28] AHARON M, ELAD M, BRUCKSTEIN A. rmK-SVD: an algorithm for designing overcomplete dictionaries for sparse representation [J].IEEETrans.SignalProcess., 2006, 54(11):4311-4322.

王儒敬(1964-)，男，安徽亳州人，博士，研究員，2004年于中國(guó)科技大學(xué)獲博士學(xué)位,主要從事農(nóng)業(yè)智能系統(tǒng)的理論、方法與技術(shù)的研究。

E-mail: rjwang@iim.ac.cn陳天嬌(1989-)，女，安徽宿州人，碩士,2014年于安徽大學(xué)獲得碩士學(xué)位，主要從事模式識(shí)別等方面的研究。

E-mail: 708877137@qq.com

Soil Near-infrared Spectroscopy Prediction Model Based on Deep Sparse Learning

WANG Ru-jing, CHEN Tian-jiao*, WANG Yu-bing, WANG Liu-san, XIE Cheng-jun, ZHANG Jie, LI Rui, CHEN Hong-bo

(InstituteofIntelligentMachines,ChineseAcademyofSciences,Hefei230031,China)
*CorrespondingAuthor,E-mail: 708877137@qq.com

This paper presents a soil near-infrared spectroscopy prediction model based on sparse representation and radial basis function neural. The model first makes the soil near-infrared large spectroscopy data to be sparse, then the model uses radial basis function neural network with sparse representation coefficients as input and the measured soil composition value by chemical methods as output to establish effective nonlinear predictive model of soil organic matter, available phosphorus and potassium respectively. The results show that the model is feasible to predict soil organic matter content, but the model needs to be further optimized on the soil phosphorus or potassium effective prediction.

near-infrared spectroscopy; deep sparse learning; neural network model

2016-06-29；

2016-08-23

中國(guó)科學(xué)院科技服務(wù)網(wǎng)絡(luò)計(jì)劃(KFJ-EW-STS-069)； 國(guó)家自然科學(xué)基金(31671586)資助項(xiàng)目 Supported by Science and Technology Service Network Initiative of Chinese Science Academy(KFJ-EW-STS-069); National Natural Science foundation of China(31671586)

1000-7032(2017)01-0109-08

O235

A

10.3788/fgxb20173801.0109