周子立+馮雷

摘要：為了滿足消費(fèi)市場(chǎng)對(duì)大米品質(zhì)快速實(shí)時(shí)檢測(cè)的需要，對(duì)未經(jīng)過(guò)粉碎、過(guò)篩等處理的大米，采集了4000cm-1—12500cm-1范圍的近紅外光譜。通過(guò)計(jì)算機(jī)數(shù)據(jù)處理，研究了大米的光譜數(shù)據(jù)，建立了一階微分光譜數(shù)據(jù)與淀粉含量和蛋白質(zhì)含量之間的偏最小二乘回歸法定量分析模型。試驗(yàn)分析表明：淀粉的預(yù)測(cè)相關(guān)系數(shù)為0.912，預(yù)測(cè)標(biāo)準(zhǔn)偏差SEP為0.084，預(yù)測(cè)均方根誤差為RMSEP為0.058；蛋白質(zhì)的預(yù)測(cè)相關(guān)系數(shù)為0.898，SEP為0.095，RMSEP為0.121。表明采用一階微分光譜數(shù)據(jù)分析法可以較好地預(yù)測(cè)經(jīng)簡(jiǎn)單處理后的大米中淀粉含量和蛋白質(zhì)含量，該結(jié)論為日后大米品質(zhì)的快速特性光譜測(cè)量奠定了基礎(chǔ)。

關(guān)鍵詞：近紅外光譜；數(shù)據(jù)處理；偏最小二乘；淀粉；蛋白質(zhì)

中圖分類號(hào)：S123,S15 文獻(xiàn)標(biāo)識(shí)碼：A

Determination of Starch Contents and Protein Contents in Rice Using Spectroscopy

Combined with Data Processing Technique

ZHOU Zili1,FENG Lei2

(1.Zhejiang Institute of Mechanical & Electrical Engineering,Hangzhou 310053,China;

2.College of Biosystems Engineering and Food Science,Zhejiang University,Hangzhou 310029,China)

Abstract:The aim of this study is to investigate the potential of the near infrared spectroscopy technique for non-destructive measurement of rice starch contents and protein contents to meet the needs of consumer markets.Spectral properties of 100 simply treated rice samples were analyzed.The quantitative analysis models were founded to determine the starch contents and protein contents.The correlation coefficients are 0.912 and 0.898,standard errors of prediction are 0.084 and 0.095, and root mean standard errors of prediction are 0.058 and 0.121, respectively.The results show that it can predict starch contents and protein contents commendably with the first derivate spectra.This research lays the basis for field experiment by using infrared spectroscopy technique to fast measurement of rice starch contents and protein contents.

Keywords:near infrared spectroscopy;DP;partial least square;starch; protein

1 引 言(Introduction)

隨著消費(fèi)市場(chǎng)對(duì)食品品質(zhì)要求的不斷提高，對(duì)大米品質(zhì)的快速檢測(cè)已成為國(guó)內(nèi)外學(xué)者研究的一個(gè)重要方向[1,2]，快速獲取大米信息符合消費(fèi)者對(duì)大米品質(zhì)了解的實(shí)際需求。以往的大米成分檢測(cè)分析即耗時(shí)、費(fèi)工、有害又對(duì)周邊環(huán)境產(chǎn)生污染，很難在田間分析測(cè)定[3]。采用光譜技術(shù)結(jié)合計(jì)算機(jī)數(shù)據(jù)處理法因其簡(jiǎn)潔、高效、低成本等優(yōu)點(diǎn)已經(jīng)受到學(xué)者們的重視，此方法還可以對(duì)多組待測(cè)試大米進(jìn)行同時(shí)測(cè)定。隨著近紅外光譜技術(shù)的不斷推廣和應(yīng)用，采用近紅外光譜技術(shù)獲取大米成分信息，再通過(guò)計(jì)算機(jī)數(shù)據(jù)處理技術(shù)對(duì)原始的光譜信息做進(jìn)一步的處理，可以快速得到大米品質(zhì)高低的鑒定結(jié)果，對(duì)保護(hù)消費(fèi)者權(quán)益，防止以次充好，規(guī)范消費(fèi)市場(chǎng)都有積極的作用。在現(xiàn)有大米成分品質(zhì)研究時(shí)大多需要一些烘干、篩選等預(yù)處理，所有研究也都是在室內(nèi)進(jìn)行，同時(shí)對(duì)大米測(cè)試樣品的要求高(粒徑小而均勻)。而對(duì)沒(méi)有進(jìn)行預(yù)處理的大米進(jìn)行測(cè)試時(shí)，其預(yù)測(cè)效果就很差。這說(shuō)明傳統(tǒng)的測(cè)試方法只能對(duì)滿足條件的室內(nèi)測(cè)試有效，不具備實(shí)用性和推廣價(jià)值，如何找到一種實(shí)用有效的針對(duì)田間實(shí)際環(huán)境中檢測(cè)大米的成分含量的方法很有必要。針對(duì)上述的問(wèn)題，本文采用計(jì)算機(jī)數(shù)據(jù)處理技術(shù)中偏最小二乘回歸分析方法定量分析大米淀粉和蛋白含量，并以此建立了一階微分光譜數(shù)據(jù)與淀粉含量和蛋白質(zhì)含量之間的定量分析模型。既滿足了室外檢測(cè)需求，又保證了研究的通用性，值得推廣。

2 材料與方法(Materials and methods)

2.1 大米樣本的獲取

從市場(chǎng)上采購(gòu)東北大米作為大米品質(zhì)測(cè)量的樣品。為保證測(cè)量的準(zhǔn)確性和代表性特選取了同一個(gè)廠家在不同日期加工生產(chǎn)的大米，共得到100個(gè)大米樣品。將大米樣品分成三份，先后用于室內(nèi)光譜分析、室外光譜分析和室內(nèi)化學(xué)分析。對(duì)于光譜分析的大米樣品，不需要進(jìn)行任何預(yù)處理，而化學(xué)分析的大米樣品則需要進(jìn)行顆粒直徑大小的統(tǒng)一篩選，費(fèi)時(shí)費(fèi)力。通過(guò)化學(xué)測(cè)量法得到大米樣本的淀粉含量分布范圍為65%-85%；蛋白質(zhì)含量的分布范圍為3.29%-3.96%。

2.2 測(cè)量?jī)x器及方法

實(shí)驗(yàn)使用美國(guó)尼高力(Nicolet)儀器公司的傅立葉變換光譜儀，波數(shù)范圍為：4000cm-1—12500cm-1，分辨率為：0.1—16cm-1。具體測(cè)試方法是：根據(jù)實(shí)驗(yàn)儀器對(duì)環(huán)境溫度、濕度的要求，實(shí)驗(yàn)室溫度應(yīng)保持在22℃左右，濕度保持在45%左右，以保證實(shí)驗(yàn)數(shù)據(jù)的可靠性。先取出光譜儀中的光纖，將光纖固定在一個(gè)支架上，再把大米樣品倒入直徑為15cm的圓形玻璃器皿中后放置在同一支架上。測(cè)試時(shí)，通過(guò)旋轉(zhuǎn)圓形玻璃器皿從0°、90°、180°、270°四個(gè)不同角度對(duì)大米樣品進(jìn)行光譜數(shù)據(jù)的獲取，再通過(guò)運(yùn)行光譜儀自帶的OMNIC6.1軟件進(jìn)行光譜數(shù)據(jù)處理并將處理后的數(shù)據(jù)保存起來(lái)，以便后期對(duì)比驗(yàn)證。實(shí)際測(cè)試時(shí)每個(gè)角度掃描三次，取其平均值作為該大米樣品的光譜數(shù)據(jù)值，這樣做的好處是可以減少圓形玻璃器皿對(duì)光譜測(cè)量的干擾。

2.3 定量分析方法

將100個(gè)大米樣品分成兩部分，70個(gè)樣品作為校正樣品用于測(cè)試數(shù)據(jù)并以此建立分析模型。30個(gè)樣品作為預(yù)測(cè)樣品用于檢驗(yàn)?zāi)Ｐ偷恼_性和測(cè)試效率。大米樣品的選取沒(méi)有任何限制，只要保證其溫度濕度滿足測(cè)試要求即可。

采用光譜技術(shù)結(jié)合數(shù)據(jù)分析技術(shù)測(cè)定大米成分的關(guān)鍵所在是建立定量分析模型。早前的大米品質(zhì)分析中運(yùn)用的多元線性回歸法[4](MLR)、主成分回歸法[5](PCR)等方法。雖然方法簡(jiǎn)單，物理意義明確，但這些方法多受使用變量數(shù)的限制，在建立校正模型時(shí)分析過(guò)程中會(huì)丟失許多光譜數(shù)據(jù)信息?；蚴窃诮⒛Ｐ瓦^(guò)程中并不能保證大米樣品的主成分一定與被測(cè)組分或性質(zhì)相關(guān)。人工神經(jīng)網(wǎng)絡(luò)(ANN)功能強(qiáng)大，但由于計(jì)算量大，而且輸入神經(jīng)元數(shù)不能過(guò)多等因素，使得ANN的使用通常需要結(jié)合主成分分析(PCA)等方法進(jìn)行，這樣做的結(jié)果就是計(jì)算量更大。本文采用偏最小二乘(PLS)方法，經(jīng)過(guò)對(duì)光譜測(cè)試數(shù)據(jù)進(jìn)行處理就可以對(duì)簡(jiǎn)單處理過(guò)的大米樣品做大米淀粉含量和蛋白質(zhì)含量的分析并建立與其光譜吸光度之間的定量分析模型。確定最佳主成分?jǐn)?shù)可以通過(guò)計(jì)算交互驗(yàn)證法(Cross-Validation)結(jié)合預(yù)測(cè)殘差平方和(PRESS)的方法來(lái)實(shí)現(xiàn)。評(píng)定分析模型效果好壞的標(biāo)準(zhǔn)為：好的模型需要有較低的建模和預(yù)測(cè)均方根誤差(RMSEC和RMSEP)，以及較高的相關(guān)系數(shù)(r)。而為了避免建立分析模型時(shí)過(guò)于擬合，則要求預(yù)測(cè)標(biāo)準(zhǔn)偏差(SEP)和建模標(biāo)準(zhǔn)偏差(SEC)盡量接近。

endprint

3 試驗(yàn)結(jié)果與分析

依據(jù)上述測(cè)試方式得到的大米樣品光譜圖如圖1所示。圖1中顯示，隨著光譜波數(shù)的不斷增加，大米樣品吸光度數(shù)值逐漸增大，而且在5000cm-1和7000cm-1左右處出現(xiàn)明顯的由于大米成分不同引起的光波吸收峰現(xiàn)象。

圖1 大米光譜圖

Fig.1 The spectrum of the rice

3.1 大米淀粉含量的預(yù)測(cè)

大米淀粉含量是大米的重要組成部分，是衡量大米品質(zhì)質(zhì)量的一個(gè)重要指標(biāo)[6]。由于人為因素造成的誤差，使得在化學(xué)法測(cè)量淀粉含量時(shí)出現(xiàn)三個(gè)大米樣品無(wú)效，在淀粉含量的預(yù)測(cè)過(guò)程中，大米樣品的數(shù)量仍為30個(gè)，但參與校正的樣品總數(shù)則只有67個(gè)。先確定在4050cm-1—10000cm-1對(duì)測(cè)量到的原始光譜進(jìn)行一階微分處理和基線校正，然后建立大米淀粉含量和光譜吸光度之間的PLS模型。當(dāng)主成分?jǐn)?shù)為3時(shí)，PRESS值達(dá)到最小(0.642)。建模效果為：R為0.899；RMSEC為0.086，SEC為0.053。再用此模型對(duì)其余30個(gè)大米樣品進(jìn)行預(yù)測(cè)。預(yù)測(cè)結(jié)果為：R為0.912；RMSEC為0.084，SEC為0.058。可見(jiàn)，大米樣品的淀粉含量定量分析模型有較高的預(yù)測(cè)和校正系數(shù)，比較低的預(yù)測(cè)和校正偏差，而校正和預(yù)測(cè)之間的差距也小，實(shí)踐證明，通過(guò)此模型預(yù)測(cè)大米樣品的淀粉含量是準(zhǔn)確的。由于本文研究的大米樣品未經(jīng)過(guò)任何篩選等預(yù)處理，所建立的大米成分分析模型能夠在實(shí)際應(yīng)用中使用，本研究也更接近實(shí)際。圖2是大米淀粉含量測(cè)試校正樣品與預(yù)測(cè)樣品比較關(guān)系圖。

(a)校正樣本 (b)預(yù)測(cè)樣本

圖 2淀粉含量建模與預(yù)測(cè)情況

Fig.2 The model and prediction of starch content

2.2 大米蛋白質(zhì)含量的預(yù)測(cè)

大米蛋白作為大米主要成分之一，在近紅外光譜照射下可以直接對(duì)應(yīng)分子的振動(dòng)吸收光譜的一定頻帶，所以近紅外光譜能夠定量分析大米蛋白的含量[7]，并對(duì)不同大米樣品的蛋白含量進(jìn)行準(zhǔn)確預(yù)測(cè)。在運(yùn)用PLS建立分析模型前，必須對(duì)異常的大米樣品進(jìn)行剔除，以保證分析模型的可靠性。當(dāng)主成分?jǐn)?shù)為3時(shí)，PRESS值達(dá)到最低，此時(shí)建立的大米蛋白含量與光譜吸光度之間的PLS模型及大米蛋白的定量校正模型最合適。圖3是對(duì)大米樣品蛋白質(zhì)含量分析時(shí)，校正樣品的預(yù)測(cè)殘差圖。圖中顯示大米樣品的預(yù)測(cè)殘差基本上都是以零水平線為中心軸，且分布均勻，另有3個(gè)大米樣品的預(yù)測(cè)殘差值明顯偏離中心軸線，所以在定量分析時(shí)，將這3個(gè)大米樣品按異常樣品剔除。此后，對(duì)正常樣品進(jìn)行光譜基線校正、Savitzky-Golay卷積平移、一階導(dǎo)數(shù)處理等處理，在光譜波數(shù)為4050cm-1—10000cm-1，蛋白質(zhì)測(cè)量值與預(yù)測(cè)值之間的R達(dá)到0.899，SEC為0.098，RMSEC為0.099。并對(duì)校正樣品中大米的蛋白含量與大米吸光度之間的關(guān)系進(jìn)行了分析。同時(shí)，采用PLS模型對(duì)預(yù)測(cè)樣品進(jìn)行預(yù)測(cè)，預(yù)測(cè)結(jié)果的R達(dá)到0.898，SEP為0.095，RMSEP為0.121。二者的SEC和SEP值都比較低而且相關(guān)系數(shù)都比較高，因此可以認(rèn)為利用PLS模型對(duì)大米蛋白含量的預(yù)測(cè)也是比較成功的。從圖4可見(jiàn)，本研究結(jié)果數(shù)值相對(duì)較低，其主要原因大米中含水量較少，且蛋白質(zhì)含量范圍的方差較小，這些都對(duì)蛋白含量預(yù)測(cè)產(chǎn)生影響。

圖3 蛋白質(zhì)預(yù)測(cè)殘差圖

Fig.3 Protein prediction residual plot

(a)校正樣本(b)預(yù)測(cè)樣本

圖4 蛋白質(zhì)含量建模與預(yù)測(cè)情況

Fig.4The model and prediction of protein content

4 結(jié)論(Conclusion)

(1)在4000cm-1—12500cm-1的光譜波數(shù)范圍內(nèi)對(duì)未經(jīng)任何篩選處理的大米進(jìn)行了測(cè)試，通過(guò)對(duì)測(cè)試得到的原始光譜數(shù)據(jù)的處理，得到了大米樣品的光譜特性。結(jié)果表明，隨著光譜波數(shù)的增加，吸光度數(shù)值也逐漸增大，并且在5000cm-1和7000cm-1附近產(chǎn)生明顯的吸收峰，這主要由大米水分引起的。

(2)在4050cm-1—10000cm-1的波數(shù)范圍內(nèi)，對(duì)測(cè)試的大米樣品光譜數(shù)據(jù)進(jìn)行一階微分處理和基線校正處理后，應(yīng)用PLS建立了大米淀粉含量和吸光度之間的定量校正模型。校正模型的R、RMSEC和SEC分別為0.899、0.086和0.053。采用此模型對(duì)其余30個(gè)樣本進(jìn)行預(yù)測(cè)，其預(yù)測(cè)樣本的R、SEP和RMSEP分別為0.912、0.084和0.058。這表明應(yīng)用光譜技術(shù)結(jié)合數(shù)據(jù)處理技術(shù)預(yù)測(cè)大米樣品中的淀粉含量是可行的。由于本文測(cè)試的大米樣品未經(jīng)過(guò)任何預(yù)先處理，所以結(jié)果更接近實(shí)際情況。

(3)在4050cm-1—10000cm-1的波數(shù)范圍內(nèi)，對(duì)測(cè)試的大米樣品光譜數(shù)據(jù)經(jīng)基線校正、一階導(dǎo)數(shù)處理和Savitzky-Golay卷積平移后，應(yīng)用PLS建立了大米樣品的蛋白含量與吸光度之間的定量校正模型。校正模型的R、SEC和RMSEC分別為0.899、0.098和0.099。采用此模型對(duì)預(yù)測(cè)大米樣品的蛋白含量進(jìn)行預(yù)測(cè)，預(yù)測(cè)的R、SEP和RMSEP為0.898、0.095和0.121。這表明應(yīng)用光譜技術(shù)結(jié)合數(shù)據(jù)處理技術(shù)可準(zhǔn)確地預(yù)測(cè)未經(jīng)過(guò)預(yù)先篩選處理的大米蛋白含量。

(4)本文介紹的光譜技術(shù)結(jié)合數(shù)據(jù)處理技術(shù)分析研究大米淀粉和蛋白含量，通過(guò)分析模型的建立可以準(zhǔn)確的預(yù)測(cè)出未經(jīng)任何篩選處理的大米樣品淀粉和蛋白含量，這在農(nóng)副產(chǎn)品的成分品質(zhì)鑒別中有其十分重要的實(shí)際作用和推廣價(jià)值。

參考文獻(xiàn) (References)

[1] 胡濤.模式識(shí)別在大米品質(zhì)檢測(cè)中的應(yīng)用[J].農(nóng)業(yè)網(wǎng)絡(luò)信息,2010,10:44-46.

[2] 陳建偉.基于機(jī)器視覺(jué)技術(shù)大米品質(zhì)檢測(cè)[J].糧食與油脂,2008,5:39-41.

[3] 何劍飛,陳召桂.不同品種的大米品質(zhì)評(píng)價(jià)分析[J].福建稻麥科技,2009,27(4):31-33.

[4] 張嫻,等.正交信號(hào)校正應(yīng)用于多元線性回歸建模的研究[J].光譜學(xué)與光譜分析,2011,31(12):3228-3231.

[5] 劉旭華,等.有監(jiān)督主成分回歸法在近紅外光譜定量分析中的應(yīng)用研究[J].光譜學(xué)與光析,2009,11:2959-2961.

[6] 寧希燁,劉亞偉,劉曉峰.大米淀粉特性及修飾作用對(duì)其食用功能的影響[J].糧食與飼料工業(yè),2011,3:35-38.

[7] 王遠(yuǎn)宏,常若葵,張偉玉.基于近紅外光譜的大米蛋白質(zhì)含量的研究[J].農(nóng)產(chǎn)品加工學(xué)刊,2009,7:30-32.

作者簡(jiǎn)介:

周子立(1963-),男,本科,高級(jí)工程師,副教授.研究領(lǐng)域:計(jì)算機(jī)信息處理及計(jì)算機(jī)應(yīng)用技術(shù).

馮 雷(1971-),男,博士,副教授.研究領(lǐng)域:光譜技術(shù),數(shù)據(jù)挖掘,協(xié)同(云)計(jì)算,移動(dòng)互聯(lián).

endprint

3 試驗(yàn)結(jié)果與分析

依據(jù)上述測(cè)試方式得到的大米樣品光譜圖如圖1所示。圖1中顯示，隨著光譜波數(shù)的不斷增加，大米樣品吸光度數(shù)值逐漸增大，而且在5000cm-1和7000cm-1左右處出現(xiàn)明顯的由于大米成分不同引起的光波吸收峰現(xiàn)象。

圖1 大米光譜圖

Fig.1 The spectrum of the rice

3.1 大米淀粉含量的預(yù)測(cè)

大米淀粉含量是大米的重要組成部分，是衡量大米品質(zhì)質(zhì)量的一個(gè)重要指標(biāo)[6]。由于人為因素造成的誤差，使得在化學(xué)法測(cè)量淀粉含量時(shí)出現(xiàn)三個(gè)大米樣品無(wú)效，在淀粉含量的預(yù)測(cè)過(guò)程中，大米樣品的數(shù)量仍為30個(gè)，但參與校正的樣品總數(shù)則只有67個(gè)。先確定在4050cm-1—10000cm-1對(duì)測(cè)量到的原始光譜進(jìn)行一階微分處理和基線校正，然后建立大米淀粉含量和光譜吸光度之間的PLS模型。當(dāng)主成分?jǐn)?shù)為3時(shí)，PRESS值達(dá)到最小(0.642)。建模效果為：R為0.899；RMSEC為0.086，SEC為0.053。再用此模型對(duì)其余30個(gè)大米樣品進(jìn)行預(yù)測(cè)。預(yù)測(cè)結(jié)果為：R為0.912；RMSEC為0.084，SEC為0.058?？梢?jiàn)，大米樣品的淀粉含量定量分析模型有較高的預(yù)測(cè)和校正系數(shù)，比較低的預(yù)測(cè)和校正偏差，而校正和預(yù)測(cè)之間的差距也小，實(shí)踐證明，通過(guò)此模型預(yù)測(cè)大米樣品的淀粉含量是準(zhǔn)確的。由于本文研究的大米樣品未經(jīng)過(guò)任何篩選等預(yù)處理，所建立的大米成分分析模型能夠在實(shí)際應(yīng)用中使用，本研究也更接近實(shí)際。圖2是大米淀粉含量測(cè)試校正樣品與預(yù)測(cè)樣品比較關(guān)系圖。

(a)校正樣本 (b)預(yù)測(cè)樣本

圖 2淀粉含量建模與預(yù)測(cè)情況

Fig.2 The model and prediction of starch content

2.2 大米蛋白質(zhì)含量的預(yù)測(cè)

大米蛋白作為大米主要成分之一，在近紅外光譜照射下可以直接對(duì)應(yīng)分子的振動(dòng)吸收光譜的一定頻帶，所以近紅外光譜能夠定量分析大米蛋白的含量[7]，并對(duì)不同大米樣品的蛋白含量進(jìn)行準(zhǔn)確預(yù)測(cè)。在運(yùn)用PLS建立分析模型前，必須對(duì)異常的大米樣品進(jìn)行剔除，以保證分析模型的可靠性。當(dāng)主成分?jǐn)?shù)為3時(shí)，PRESS值達(dá)到最低，此時(shí)建立的大米蛋白含量與光譜吸光度之間的PLS模型及大米蛋白的定量校正模型最合適。圖3是對(duì)大米樣品蛋白質(zhì)含量分析時(shí)，校正樣品的預(yù)測(cè)殘差圖。圖中顯示大米樣品的預(yù)測(cè)殘差基本上都是以零水平線為中心軸，且分布均勻，另有3個(gè)大米樣品的預(yù)測(cè)殘差值明顯偏離中心軸線，所以在定量分析時(shí)，將這3個(gè)大米樣品按異常樣品剔除。此后，對(duì)正常樣品進(jìn)行光譜基線校正、Savitzky-Golay卷積平移、一階導(dǎo)數(shù)處理等處理，在光譜波數(shù)為4050cm-1—10000cm-1，蛋白質(zhì)測(cè)量值與預(yù)測(cè)值之間的R達(dá)到0.899，SEC為0.098，RMSEC為0.099。并對(duì)校正樣品中大米的蛋白含量與大米吸光度之間的關(guān)系進(jìn)行了分析。同時(shí)，采用PLS模型對(duì)預(yù)測(cè)樣品進(jìn)行預(yù)測(cè)，預(yù)測(cè)結(jié)果的R達(dá)到0.898，SEP為0.095，RMSEP為0.121。二者的SEC和SEP值都比較低而且相關(guān)系數(shù)都比較高，因此可以認(rèn)為利用PLS模型對(duì)大米蛋白含量的預(yù)測(cè)也是比較成功的。從圖4可見(jiàn)，本研究結(jié)果數(shù)值相對(duì)較低，其主要原因大米中含水量較少，且蛋白質(zhì)含量范圍的方差較小，這些都對(duì)蛋白含量預(yù)測(cè)產(chǎn)生影響。

圖3 蛋白質(zhì)預(yù)測(cè)殘差圖

Fig.3 Protein prediction residual plot

(a)校正樣本(b)預(yù)測(cè)樣本

圖4 蛋白質(zhì)含量建模與預(yù)測(cè)情況

Fig.4The model and prediction of protein content

4 結(jié)論(Conclusion)

(1)在4000cm-1—12500cm-1的光譜波數(shù)范圍內(nèi)對(duì)未經(jīng)任何篩選處理的大米進(jìn)行了測(cè)試，通過(guò)對(duì)測(cè)試得到的原始光譜數(shù)據(jù)的處理，得到了大米樣品的光譜特性。結(jié)果表明，隨著光譜波數(shù)的增加，吸光度數(shù)值也逐漸增大，并且在5000cm-1和7000cm-1附近產(chǎn)生明顯的吸收峰，這主要由大米水分引起的。

(2)在4050cm-1—10000cm-1的波數(shù)范圍內(nèi)，對(duì)測(cè)試的大米樣品光譜數(shù)據(jù)進(jìn)行一階微分處理和基線校正處理后，應(yīng)用PLS建立了大米淀粉含量和吸光度之間的定量校正模型。校正模型的R、RMSEC和SEC分別為0.899、0.086和0.053。采用此模型對(duì)其余30個(gè)樣本進(jìn)行預(yù)測(cè)，其預(yù)測(cè)樣本的R、SEP和RMSEP分別為0.912、0.084和0.058。這表明應(yīng)用光譜技術(shù)結(jié)合數(shù)據(jù)處理技術(shù)預(yù)測(cè)大米樣品中的淀粉含量是可行的。由于本文測(cè)試的大米樣品未經(jīng)過(guò)任何預(yù)先處理，所以結(jié)果更接近實(shí)際情況。

(3)在4050cm-1—10000cm-1的波數(shù)范圍內(nèi)，對(duì)測(cè)試的大米樣品光譜數(shù)據(jù)經(jīng)基線校正、一階導(dǎo)數(shù)處理和Savitzky-Golay卷積平移后，應(yīng)用PLS建立了大米樣品的蛋白含量與吸光度之間的定量校正模型。校正模型的R、SEC和RMSEC分別為0.899、0.098和0.099。采用此模型對(duì)預(yù)測(cè)大米樣品的蛋白含量進(jìn)行預(yù)測(cè)，預(yù)測(cè)的R、SEP和RMSEP為0.898、0.095和0.121。這表明應(yīng)用光譜技術(shù)結(jié)合數(shù)據(jù)處理技術(shù)可準(zhǔn)確地預(yù)測(cè)未經(jīng)過(guò)預(yù)先篩選處理的大米蛋白含量。

(4)本文介紹的光譜技術(shù)結(jié)合數(shù)據(jù)處理技術(shù)分析研究大米淀粉和蛋白含量，通過(guò)分析模型的建立可以準(zhǔn)確的預(yù)測(cè)出未經(jīng)任何篩選處理的大米樣品淀粉和蛋白含量，這在農(nóng)副產(chǎn)品的成分品質(zhì)鑒別中有其十分重要的實(shí)際作用和推廣價(jià)值。

參考文獻(xiàn) (References)

[1] 胡濤.模式識(shí)別在大米品質(zhì)檢測(cè)中的應(yīng)用[J].農(nóng)業(yè)網(wǎng)絡(luò)信息,2010,10:44-46.

[2] 陳建偉.基于機(jī)器視覺(jué)技術(shù)大米品質(zhì)檢測(cè)[J].糧食與油脂,2008,5:39-41.

[3] 何劍飛,陳召桂.不同品種的大米品質(zhì)評(píng)價(jià)分析[J].福建稻麥科技,2009,27(4):31-33.

[4] 張嫻,等.正交信號(hào)校正應(yīng)用于多元線性回歸建模的研究[J].光譜學(xué)與光譜分析,2011,31(12):3228-3231.

[5] 劉旭華,等.有監(jiān)督主成分回歸法在近紅外光譜定量分析中的應(yīng)用研究[J].光譜學(xué)與光析,2009,11:2959-2961.

[6] 寧希燁,劉亞偉,劉曉峰.大米淀粉特性及修飾作用對(duì)其食用功能的影響[J].糧食與飼料工業(yè),2011,3:35-38.

[7] 王遠(yuǎn)宏,常若葵,張偉玉.基于近紅外光譜的大米蛋白質(zhì)含量的研究[J].農(nóng)產(chǎn)品加工學(xué)刊,2009,7:30-32.

作者簡(jiǎn)介:

周子立(1963-),男,本科,高級(jí)工程師,副教授.研究領(lǐng)域:計(jì)算機(jī)信息處理及計(jì)算機(jī)應(yīng)用技術(shù).

馮 雷(1971-),男,博士,副教授.研究領(lǐng)域:光譜技術(shù),數(shù)據(jù)挖掘,協(xié)同(云)計(jì)算,移動(dòng)互聯(lián).

endprint

3 試驗(yàn)結(jié)果與分析

依據(jù)上述測(cè)試方式得到的大米樣品光譜圖如圖1所示。圖1中顯示，隨著光譜波數(shù)的不斷增加，大米樣品吸光度數(shù)值逐漸增大，而且在5000cm-1和7000cm-1左右處出現(xiàn)明顯的由于大米成分不同引起的光波吸收峰現(xiàn)象。

圖1 大米光譜圖

Fig.1 The spectrum of the rice

3.1 大米淀粉含量的預(yù)測(cè)

大米淀粉含量是大米的重要組成部分，是衡量大米品質(zhì)質(zhì)量的一個(gè)重要指標(biāo)[6]。由于人為因素造成的誤差，使得在化學(xué)法測(cè)量淀粉含量時(shí)出現(xiàn)三個(gè)大米樣品無(wú)效，在淀粉含量的預(yù)測(cè)過(guò)程中，大米樣品的數(shù)量仍為30個(gè)，但參與校正的樣品總數(shù)則只有67個(gè)。先確定在4050cm-1—10000cm-1對(duì)測(cè)量到的原始光譜進(jìn)行一階微分處理和基線校正，然后建立大米淀粉含量和光譜吸光度之間的PLS模型。當(dāng)主成分?jǐn)?shù)為3時(shí)，PRESS值達(dá)到最小(0.642)。建模效果為：R為0.899；RMSEC為0.086，SEC為0.053。再用此模型對(duì)其余30個(gè)大米樣品進(jìn)行預(yù)測(cè)。預(yù)測(cè)結(jié)果為：R為0.912；RMSEC為0.084，SEC為0.058?？梢?jiàn)，大米樣品的淀粉含量定量分析模型有較高的預(yù)測(cè)和校正系數(shù)，比較低的預(yù)測(cè)和校正偏差，而校正和預(yù)測(cè)之間的差距也小，實(shí)踐證明，通過(guò)此模型預(yù)測(cè)大米樣品的淀粉含量是準(zhǔn)確的。由于本文研究的大米樣品未經(jīng)過(guò)任何篩選等預(yù)處理，所建立的大米成分分析模型能夠在實(shí)際應(yīng)用中使用，本研究也更接近實(shí)際。圖2是大米淀粉含量測(cè)試校正樣品與預(yù)測(cè)樣品比較關(guān)系圖。

(a)校正樣本 (b)預(yù)測(cè)樣本

圖 2淀粉含量建模與預(yù)測(cè)情況

Fig.2 The model and prediction of starch content

2.2 大米蛋白質(zhì)含量的預(yù)測(cè)

大米蛋白作為大米主要成分之一，在近紅外光譜照射下可以直接對(duì)應(yīng)分子的振動(dòng)吸收光譜的一定頻帶，所以近紅外光譜能夠定量分析大米蛋白的含量[7]，并對(duì)不同大米樣品的蛋白含量進(jìn)行準(zhǔn)確預(yù)測(cè)。在運(yùn)用PLS建立分析模型前，必須對(duì)異常的大米樣品進(jìn)行剔除，以保證分析模型的可靠性。當(dāng)主成分?jǐn)?shù)為3時(shí)，PRESS值達(dá)到最低，此時(shí)建立的大米蛋白含量與光譜吸光度之間的PLS模型及大米蛋白的定量校正模型最合適。圖3是對(duì)大米樣品蛋白質(zhì)含量分析時(shí)，校正樣品的預(yù)測(cè)殘差圖。圖中顯示大米樣品的預(yù)測(cè)殘差基本上都是以零水平線為中心軸，且分布均勻，另有3個(gè)大米樣品的預(yù)測(cè)殘差值明顯偏離中心軸線，所以在定量分析時(shí)，將這3個(gè)大米樣品按異常樣品剔除。此后，對(duì)正常樣品進(jìn)行光譜基線校正、Savitzky-Golay卷積平移、一階導(dǎo)數(shù)處理等處理，在光譜波數(shù)為4050cm-1—10000cm-1，蛋白質(zhì)測(cè)量值與預(yù)測(cè)值之間的R達(dá)到0.899，SEC為0.098，RMSEC為0.099。并對(duì)校正樣品中大米的蛋白含量與大米吸光度之間的關(guān)系進(jìn)行了分析。同時(shí)，采用PLS模型對(duì)預(yù)測(cè)樣品進(jìn)行預(yù)測(cè)，預(yù)測(cè)結(jié)果的R達(dá)到0.898，SEP為0.095，RMSEP為0.121。二者的SEC和SEP值都比較低而且相關(guān)系數(shù)都比較高，因此可以認(rèn)為利用PLS模型對(duì)大米蛋白含量的預(yù)測(cè)也是比較成功的。從圖4可見(jiàn)，本研究結(jié)果數(shù)值相對(duì)較低，其主要原因大米中含水量較少，且蛋白質(zhì)含量范圍的方差較小，這些都對(duì)蛋白含量預(yù)測(cè)產(chǎn)生影響。

圖3 蛋白質(zhì)預(yù)測(cè)殘差圖

Fig.3 Protein prediction residual plot

(a)校正樣本(b)預(yù)測(cè)樣本

圖4 蛋白質(zhì)含量建模與預(yù)測(cè)情況

Fig.4The model and prediction of protein content

4 結(jié)論(Conclusion)

(1)在4000cm-1—12500cm-1的光譜波數(shù)范圍內(nèi)對(duì)未經(jīng)任何篩選處理的大米進(jìn)行了測(cè)試，通過(guò)對(duì)測(cè)試得到的原始光譜數(shù)據(jù)的處理，得到了大米樣品的光譜特性。結(jié)果表明，隨著光譜波數(shù)的增加，吸光度數(shù)值也逐漸增大，并且在5000cm-1和7000cm-1附近產(chǎn)生明顯的吸收峰，這主要由大米水分引起的。

(2)在4050cm-1—10000cm-1的波數(shù)范圍內(nèi)，對(duì)測(cè)試的大米樣品光譜數(shù)據(jù)進(jìn)行一階微分處理和基線校正處理后，應(yīng)用PLS建立了大米淀粉含量和吸光度之間的定量校正模型。校正模型的R、RMSEC和SEC分別為0.899、0.086和0.053。采用此模型對(duì)其余30個(gè)樣本進(jìn)行預(yù)測(cè)，其預(yù)測(cè)樣本的R、SEP和RMSEP分別為0.912、0.084和0.058。這表明應(yīng)用光譜技術(shù)結(jié)合數(shù)據(jù)處理技術(shù)預(yù)測(cè)大米樣品中的淀粉含量是可行的。由于本文測(cè)試的大米樣品未經(jīng)過(guò)任何預(yù)先處理，所以結(jié)果更接近實(shí)際情況。

(3)在4050cm-1—10000cm-1的波數(shù)范圍內(nèi)，對(duì)測(cè)試的大米樣品光譜數(shù)據(jù)經(jīng)基線校正、一階導(dǎo)數(shù)處理和Savitzky-Golay卷積平移后，應(yīng)用PLS建立了大米樣品的蛋白含量與吸光度之間的定量校正模型。校正模型的R、SEC和RMSEC分別為0.899、0.098和0.099。采用此模型對(duì)預(yù)測(cè)大米樣品的蛋白含量進(jìn)行預(yù)測(cè)，預(yù)測(cè)的R、SEP和RMSEP為0.898、0.095和0.121。這表明應(yīng)用光譜技術(shù)結(jié)合數(shù)據(jù)處理技術(shù)可準(zhǔn)確地預(yù)測(cè)未經(jīng)過(guò)預(yù)先篩選處理的大米蛋白含量。

(4)本文介紹的光譜技術(shù)結(jié)合數(shù)據(jù)處理技術(shù)分析研究大米淀粉和蛋白含量，通過(guò)分析模型的建立可以準(zhǔn)確的預(yù)測(cè)出未經(jīng)任何篩選處理的大米樣品淀粉和蛋白含量，這在農(nóng)副產(chǎn)品的成分品質(zhì)鑒別中有其十分重要的實(shí)際作用和推廣價(jià)值。

參考文獻(xiàn) (References)

[1] 胡濤.模式識(shí)別在大米品質(zhì)檢測(cè)中的應(yīng)用[J].農(nóng)業(yè)網(wǎng)絡(luò)信息,2010,10:44-46.

[2] 陳建偉.基于機(jī)器視覺(jué)技術(shù)大米品質(zhì)檢測(cè)[J].糧食與油脂,2008,5:39-41.

[3] 何劍飛,陳召桂.不同品種的大米品質(zhì)評(píng)價(jià)分析[J].福建稻麥科技,2009,27(4):31-33.

[4] 張嫻,等.正交信號(hào)校正應(yīng)用于多元線性回歸建模的研究[J].光譜學(xué)與光譜分析,2011,31(12):3228-3231.

[5] 劉旭華,等.有監(jiān)督主成分回歸法在近紅外光譜定量分析中的應(yīng)用研究[J].光譜學(xué)與光析,2009,11:2959-2961.

[6] 寧希燁,劉亞偉,劉曉峰.大米淀粉特性及修飾作用對(duì)其食用功能的影響[J].糧食與飼料工業(yè),2011,3:35-38.

[7] 王遠(yuǎn)宏,常若葵,張偉玉.基于近紅外光譜的大米蛋白質(zhì)含量的研究[J].農(nóng)產(chǎn)品加工學(xué)刊,2009,7:30-32.

作者簡(jiǎn)介:

周子立(1963-),男,本科,高級(jí)工程師,副教授.研究領(lǐng)域:計(jì)算機(jī)信息處理及計(jì)算機(jī)應(yīng)用技術(shù).

馮 雷(1971-),男,博士,副教授.研究領(lǐng)域:光譜技術(shù),數(shù)據(jù)挖掘,協(xié)同(云)計(jì)算,移動(dòng)互聯(lián).

endprint