稻谷脂肪近紅外光譜特征篩選及檢測模型構(gòu)建

李路1，黃漢英1*，李毅1，趙思明2，楊素仙1

1(華中農(nóng)業(yè)大學(xué) 工學(xué)院，湖北 武漢，430070) 2(華中農(nóng)業(yè)大學(xué) 食品科技學(xué)院，湖北 武漢，430070)

摘 要 應(yīng)用近紅外光譜技術(shù)對稻谷脂肪含量進(jìn)行檢測。采集了90個稻谷樣本的漫反射近紅外光譜，運(yùn)用Kennard-Stone法選取校正集及預(yù)測集樣本。對比研究了歸一化、一階導(dǎo)、二階導(dǎo)、一階導(dǎo)+歸一化等4種預(yù)處理方法對模型性能的影響，確定一階導(dǎo)為優(yōu)良預(yù)處理方法。運(yùn)用競爭性自適應(yīng)重加權(quán)采樣技術(shù)篩選出與稻谷脂肪含量檢測相關(guān)的特征波長，再用多元線性回歸對特征波長進(jìn)行優(yōu)選，最終得到30個特征波長。其中最典型的特征波長為1 343、1 489和1 583 nm，反映了稻谷脂肪中大量存在的—CH和—OH基團(tuán)。所建立的基于近紅外光譜分析技術(shù)的稻谷脂肪含量檢測模型具的決定系數(shù)為0.958 9，定標(biāo)標(biāo)準(zhǔn)差RMSEC為0.223 6，相對偏差為5.53%。

關(guān)鍵詞 近紅外光譜；稻谷；脂肪；競爭性自適應(yīng)重加權(quán)采樣；多元線性回歸

脂肪含量是評價稻谷品質(zhì)的重要指標(biāo)[1]。而其傳統(tǒng)測定方法存在過程繁瑣、條件不易控制、藥品消耗大、費(fèi)時費(fèi)力等不足。近紅外光譜分析技術(shù)以其檢測過程簡單迅速、耗材少、無損、分析重現(xiàn)性好、成本低等優(yōu)點(diǎn)，在糧食、蔬果、肉制品等方面的檢測中得到了廣泛的應(yīng)用[2-4]。

目前國內(nèi)外學(xué)者利用近紅外分析技術(shù)做了一些食品檢測方面的研究，BAGCHI等[5]、SUN等[6]和XIE等[7]建立了稻米中蛋白質(zhì)及淀粉含量的檢測模型。SIRISOMBOON等[8]研究了稻谷中黃曲霉素的近紅外檢測方法。張中衛(wèi)等[9]、郭中華等[10]和趙明富等[11]建立了牛奶中脂肪和蛋白質(zhì)含量的數(shù)學(xué)模型。林家永等[12]對不同品種和儲存期限的大米水分、脂肪酸值進(jìn)行檢測。上述研究所建立的數(shù)學(xué)模型對食品主要營養(yǎng)成分的檢測都能取得較好的精度，檢測結(jié)果的決定系數(shù)R2大多在0.8以上，但是模型的建立思路主要集中在2個方面：(1)根據(jù)食品營養(yǎng)成分中特殊化學(xué)基團(tuán)所對應(yīng)的近紅外光譜特征波長建立模型，具有模型簡單，物理意義明確的優(yōu)點(diǎn)，但檢測結(jié)果不穩(wěn)定，精度不高；(2)利用近紅外光譜全波長建立模型，具有適應(yīng)性廣，精度理想的優(yōu)點(diǎn)，但模型復(fù)雜，計算量較大。

本文針對模型性能與復(fù)雜度的矛盾關(guān)系，采用競爭性自適應(yīng)重加權(quán)采樣(competitive adaptive reweighted sampling，CARS)[13]方法，在稻谷近紅外光譜的全波段，篩選出適量的特征波長，然后采用多元線性回歸法建立稻谷脂肪含量的檢測模型，在保證較好模型性能的前提下盡量減小計算量，為快速檢測稻谷主要營養(yǎng)成分提供新的方法。

1 材料與方法

1.1 試驗材料

試驗材料包含A4A/R326、巨風(fēng)A/R257、廣占S/R166、中9A/R591、岡紅1A/R15等46個品種的稻谷，由黃岡農(nóng)科院提供，產(chǎn)自海南省。袋裝封存于陰涼處，收獲至試驗在2個月內(nèi)完成。將部分品種按質(zhì)量1∶1的比例兩兩混合，得到44個混合樣本，共計90個稻谷樣本。

1.2 試驗設(shè)備與樣本光譜采集

近紅外光譜儀：漫反射式Supnir-2720近紅外光譜儀，杭州聚光科技股份有限公司。其測定范圍為1 000～1 799 nm，光譜采樣間隔為1 nm，光譜分辨率為10 nm。為減小誤差，每個樣本掃描3次。稻谷籽粒置于樣本盤內(nèi)，且裝滿、壓實，逐一掃描。

1.3 稻谷脂肪含量測定

按照GB/T 5511—2008《糧油檢驗 糧食中粗脂肪含量測定》中的索氏抽提法[14]測定稻谷的脂肪含量，取3次測定的平均值。

1.4 數(shù)據(jù)處理方法

1.4.1 樣本集劃分

運(yùn)用Kennard-Stone方法[15]對90個稻谷樣本進(jìn)行分組，80%為校正集，20%為驗證集，分別用于模型的建立與驗證。Kennard-Stone方法基于變量之間的歐氏距離，在樣本光譜的特征空間中均勻選取樣本，依次選取歐氏距離最遠(yuǎn)的點(diǎn)，進(jìn)入校正集，留下馬氏距離居中的點(diǎn)在驗證集之中[16]。這樣可使光譜差異較大的樣本全部進(jìn)入校正集，從而在一定程度上避免了校正集樣本分布的不均勻。

1.4.2 近紅外光譜預(yù)處理

預(yù)處理能將光譜的有效信息放大，且將光譜壓縮在相同的范圍進(jìn)行比較，消除了基線漂移、樣品不均勻、光散射、光程變化等對光譜的影響。選擇合適的預(yù)處理方法，對提高模型的檢測能力和精度十分重要。本文運(yùn)用歸一化、一階導(dǎo)、二階導(dǎo)等預(yù)處理方法中的一種或多種的組合對稻谷近紅外光譜做預(yù)處理。然后根據(jù)不同預(yù)處理方法所建立模型的性能指標(biāo)來確定最終的光譜預(yù)處理方案。

1.4.2.1 歸一化

在建立近紅外定量模型前，為了將光譜的有效信息放大，常采用歸一化來進(jìn)行預(yù)處理，使光譜數(shù)據(jù)的方差為1，均值為0，將光譜數(shù)據(jù)放在相同的范圍進(jìn)行比較[17]。本文運(yùn)用Z-score歸一化進(jìn)行預(yù)處理試驗。

1.4.2.2 一階導(dǎo)

導(dǎo)數(shù)法可消除基線偏移、背景干擾出現(xiàn)的譜線重疊，呈現(xiàn)變化明顯的波峰波谷，提供比原始光譜更清晰的光譜輪廓變化[18]。

1.4.2.3 二階導(dǎo)

原始光譜經(jīng)過二階導(dǎo)數(shù)處理后，會變得很尖銳，有利于更好的確定波峰和波谷的位置，可使各吸收峰變得更容易區(qū)分。二階導(dǎo)數(shù)可以消除基線的線性傾斜[19]。

1.4.3 特征波長篩選

稻谷脂肪的近紅外光譜波長點(diǎn)數(shù)為800個，而樣本數(shù)是90個，利用此數(shù)據(jù)進(jìn)行回歸分析，共線性非常嚴(yán)重。利用CARS方法[20]對特征波長進(jìn)行篩選可簡化模型，并提高模型的檢測能力。

1.4.4 建模

為了降低模型的運(yùn)算量，便于后續(xù)研究中在線快速檢測系統(tǒng)的開發(fā)，使用多元線性回歸方法建立基于近紅外光譜分析技術(shù)的稻谷脂肪含量檢測模型，以模型的決定系數(shù)(R2)、定標(biāo)標(biāo)準(zhǔn)差(RMSEC)、相對偏差來評價模型的穩(wěn)定性、檢測能力和優(yōu)劣。當(dāng)R2越趨近于1，RMSEC和相對偏差的值越趨近于0，則模型的穩(wěn)定性和檢測能力越好，在實際中的檢測越準(zhǔn)確。同時，在建模的過程中通過顯著性指標(biāo)進(jìn)一步對特征波長進(jìn)行優(yōu)選，使模型進(jìn)一步簡化。

2 結(jié)果與分析

2.1 樣本近紅外光譜

圖1為90個稻谷樣本的近紅外光譜圖?？梢?，不同樣本的近紅外光譜在總體趨勢上是一致的。但由于物質(zhì)成分含量的微小差異，不同樣本的吸光度略有不同。

圖1 90個稻谷樣本的原始近紅外光譜圖
Fig.1 Original Near-infrared spectra of ninety rice samples

2.2 樣本集劃分結(jié)果

采用Kennard-Stone法選取校正集與驗證集，結(jié)果見表1。驗證集的范圍包含在校正集內(nèi)，分組合理。

表1 Kennard-Stone 分組結(jié)果
Table1 Results of Kennard-Stone

2.3 樣本預(yù)處理

4種預(yù)處理方法后的結(jié)果如圖2和表2所示。由表2可知，運(yùn)用一階導(dǎo)處理后所建模型的決定系數(shù)最大、定標(biāo)標(biāo)準(zhǔn)差最小，因此選用一階導(dǎo)為優(yōu)良的預(yù)處理方法。

圖2 經(jīng)過4種預(yù)處理后的光譜圖
Fig.2 Spectra after 4 preprocessing methods

表2 不同預(yù)處理方法的比較
Table2 Comparison of different pretreatments

2.4 特征波長篩選

對光譜進(jìn)行預(yù)處理后，采用CARS方法對稻谷脂肪的特征波長進(jìn)行多次篩選，其篩選過程如圖3所示。圖3-a為篩選過程中選出變量的變化趨勢，隨著運(yùn)行次數(shù)的增加，保留的變量數(shù)越來越少，且由快到慢呈指數(shù)函數(shù)遞減。圖3-b為波長變量篩選過程中交叉驗證均方差RMSECV的變化趨勢，在1～53次采樣過程中，RMSECV呈現(xiàn)遞減趨勢，表明篩選過程中剔除了與樣本性質(zhì)無關(guān)的變量，此時對應(yīng)的變量數(shù)為34，53次后開始遞增，說明篩選過程中開始剔除了與稻谷脂肪含量相關(guān)的變量，從而導(dǎo)致RMSECV值增大，可見在第53次時，已將無關(guān)變量全部剔除，最后保留波長數(shù)34個。圖3-c中“*”所對應(yīng)的點(diǎn)即為RMSECV最低點(diǎn)，圖3-c中各線表示隨著運(yùn)行次數(shù)增加各波長變量回歸系數(shù)的變化趨勢。特征波長篩選時蒙特卡羅采樣次數(shù)為100，挑選出的34個特征波長見圖4。

圖3 稻谷脂肪特征波長篩選圖
Fig.3 Key wavelengths selection of fat in rice

2.5 稻谷脂肪數(shù)學(xué)模型的建立

運(yùn)用多元線性回歸方法，以34個特征波長建立稻谷脂肪含量的檢測模型，將顯著性不高的特征波長剔除，優(yōu)選特征波長，最終得到30個特征波長，回歸系數(shù)及顯著性見表3。檢測模型的具體形式為：

圖4 稻谷脂肪特征波長的分布情況(o表示特征波長)
Fig.4 Distribution of key wavelengths (shown by “o”) for fat in rice

y=b+∑aixi

(1)

式中：y，稻谷脂肪含量；b，回歸常數(shù)項；ai，各特征波長的回歸系數(shù)；xi，特征波長的吸光值經(jīng)過一階導(dǎo)數(shù)處理后的數(shù)值；i=1 000, 1 001, … , 1 799，波長所對應(yīng)的納米數(shù)。

表3為稻谷脂肪檢測模型的參數(shù)及其顯著性?？梢?，在x1 343、x1 489和x1 583處，回歸系數(shù)的絕對值最大，說明這些特征波長對模型的影響最大，t值相對較大，說明對模型的影響較顯著。其中，1 343 nm反映了C—H基團(tuán)的第二組合頻和O—H基團(tuán)的伸縮振動的二級倍頻吸收帶，1 489 nm反映了O—H基團(tuán)的伸縮振動的二級倍頻，1 583 nm反映了O—H基團(tuán)的伸縮振動的一級倍頻吸收帶。檢測模型的決定系數(shù)R2為0.958 9，定標(biāo)標(biāo)準(zhǔn)差RMSEC為0.223 6，相對偏差為5.53%，說明模型具有較好的穩(wěn)定性和準(zhǔn)確性。

表3 稻谷脂肪檢測模型的參數(shù)及其顯著性
Table 3 Parameters and significances of fat prediction model

利用上述回歸方程預(yù)測18個驗證集的脂肪含量，偏差的絕對值最大為0.716 1，最小為0.003 2，模型驗證的決定系數(shù)R2為0.825 0，校驗標(biāo)準(zhǔn)差RMSEP為0.339 8，相對偏差為8.55%，說明所建模型具有較好的檢測能力。

3 結(jié)論

本文首先采集了90個稻谷樣本的漫反射近紅外光譜，使用Kennard-Stone法選取了校正集與驗證集。然后對比了歸一化、一階導(dǎo)、二階導(dǎo)、一階導(dǎo)+歸一化等4種預(yù)處理方法對模型性能的影響，使用CARS法確定了與稻谷脂肪含量檢測相關(guān)的近紅外特征波長。最后利用多元線性回歸理論對特征波長進(jìn)行優(yōu)選，并建立了基于近紅外光譜分析技術(shù)的稻谷脂肪含量的檢測模型。結(jié)果表明，稻谷紅外光譜的優(yōu)良預(yù)處理方法為一階導(dǎo)，與稻谷脂肪含量檢測相關(guān)的近紅外特征波長為30個，其中最典型的特征波長為：1 343、1 489和1 583 nm。檢測模型的決定系數(shù)R2為0.958 9，定標(biāo)標(biāo)準(zhǔn)差RMSEC為0.223 6，相對偏差為5.53%。將檢測模型用于驗證集，偏差的絕對值最大為0.716 1，最小為0.003 2，模型驗證的決定系數(shù)R2為0.825 0，校驗標(biāo)準(zhǔn)差RMSEP為0.339 8，相對偏差為8.55%，說明所建模型具有較好的穩(wěn)定性和檢測能力。

參考文獻(xiàn)

[1] 黃麗, 柏蕓, 韓文芳, 等. 稻米質(zhì)量對食品安全的影響[J]. 中國糧油學(xué)報, 2013, 28(4): 113-117.

[2] MILICA M P, JASNA S M. Near infrared spectroscopy-advanced analytical tool in wheat breeding, trade, and processing [J]. Food and Bioprocess Technology, 2013, 6: 330-352.

[3] 陳辰, 魯曉翔, 張鵬, 等. 玫瑰香葡萄貯藏期間糖酸品質(zhì)的近紅外檢測[J]. 食品與發(fā)酵工業(yè), 2015, 41(6): 175-180.

[4] 涂斌, 陳志, 彭博, 等. 基于多源光譜特征融合技術(shù)的花生油摻偽檢測[J]. 食品與發(fā)酵工業(yè), 2014, 42(4): 169-173.

[5] BAGCHI T B, SHARMA S G, CHATTOPADHYAY K. Development of NIRS models to predict protein and amylose content of brown rice and proximate compositions of rice bran [J]. Food Chemistry, 2016, 191(1): 20-27.

[6] SUN Cheng-xiao, YU Yong-hong, DUAN Bing-wu, et al. Rapid prediction of rice quality characteristics by near-infrared reflectance spectroscopy for breeding programs [J]. Cereal Chemistry, 2014, 91(3): 270-275.

[7] XIE L H, TANG S Q, CHEN N, et al. Optimisation of near-infrared reflectance model in measuring protein and amylose content of rice flour [J]. Food Chemistry, 2014, 142: 92-100.

[8] SIRISOMBOON C D, PUTTHANG R, SIRISOMBOON P. Application of near infrared spectroscopy to detect aflatoxigenic fungal contamination in rice [J]. Food Control, 2013, 33(1): 207-214.

[9] 張中衛(wèi), 溫志渝, 曾甜玲, 等. 微型近紅外光纖光譜儀用于奶粉中蛋白質(zhì)脂肪的定量檢測研究[J]. 譜學(xué)與光譜分析, 2013, 33(7): 1 796-1 800.

[10] 郭中華, 王磊, 金靈, 等. 基于近紅外透射光譜的乳制品蛋白質(zhì)、脂肪含量檢測[J]. 光電子·激光, 2013(6): 1 163-1 168.

[11] 趙明富, 施玉佳, 羅彬彬, 等. 基于近紅外透射光譜的牛奶中脂肪和蛋白質(zhì)含量檢測[J]. 激光雜志, 2014(1): 44-45.

[12] 林家永, 范維燕, 薛雅琳, 等. 稻米儲藏品質(zhì)近紅外光譜快速判定技術(shù)及儀器研發(fā)[J]. 中國糧油學(xué)報, 2011, 26(7): 113-118.

[13] LI Hong-dong, LIANG Yi-zeng, XU Qing-song, et al. Key wavelengths screening using competitive adaptive reweighted sampling method for multivariate calibration[J]. Analytica Chimica Acta, 2009, 648(1): 77-84.

[14] GB 5009.5—2010. 食品安全國家標(biāo)準(zhǔn): 食品中蛋白質(zhì)的測定[S].

[15] 陳鑫. 基于智能算法的近紅外光譜分析預(yù)測模型建立方法研究[D]. 無錫：江南大學(xué), 2013.

[16] 張華秀, 李曉寧, 范偉, 等. 近紅外光譜結(jié)合CARS變量篩選方法用于液態(tài)奶中蛋白質(zhì)與脂肪含量的測定[J]. 分析測試學(xué)報, 2010, 29(5): 430-434.

[17] 王艷陽, 魏永霞. 基于Z比分?jǐn)?shù)的參考作物騰發(fā)量計算方法優(yōu)選[J]. 節(jié)水灌溉, 2015, 29(3): 52-54.

[18] 宋麗華. 花生籽仁蛋白質(zhì)含量近紅外光譜模型的建立及育種應(yīng)用[D]. 保定：河北農(nóng)業(yè)大學(xué), 2011.

[19] 白雁, 余振喜, 孫素琴, 等. 近紅外漫反射光譜技術(shù)測定牛膝中蛻皮甾酮[J]. 中草藥, 2005, 36(9): 115-118.

[20] 劉星, 單楊, 李高陽. CARS結(jié)合PLS-LDA法識別奶牛飼料中土霉素的可行性研究[J]. 包裝與食品機(jī)械, 2012, 30(4): 1-4.

Establishment of a selection and detection model of fat inrice by nearinfrared spectrum characteristics

LI Lu1, HUANG Han-ying 1*, LI Yi1, ZHAO Si-ming2, YANG Su-xian1

1 (College of Engineering, Huazhong Agricultural University, Wuhan 430070, China) 2 (College of Food Science & Technology, Huazhong Agricultural University, Wuhan 430070, China)

ABSTRACT Near Infrared (NIR) spectrum was used to detect fat content in rice. NIR spectra of 90 rice samples were measured. Kennard-Stone method was used to select the calibration set and prediction set samples. The effects of different pretreatment (normalize, first derivative and second derivative methods) have been compared for the accuracy of the models. The best pretreatment method is the first derivative. The competitive self-adaptive weighted sampling technology is used to screen the key wavelengths associated with sample properties. Finally, thirty key wavelengths are selected by Multiple Linear Regression further. The most typical key wavelengths are 1 343 nm, 1 489 nm and 1 583 nm which related to the groups of —CH and —OH in rice fat. The detection model of fat content of rice based on near infrared spectroscopy has higher precision with the coefficient of determination, root mean square error of calibration and relative deviation are 0.958 9, 0.223 6 and 5.53%, respectively.

Key words near infrared spectrum; rice; fat; competitive adaptive reweighted sampling; multiple linear regression

DOI:10.13995/j.cnki.11-1802/ts.014950

引用格式：李路，黃漢英，李毅,等.稻谷脂肪近紅外光譜特征篩選及檢測模型構(gòu)建[J].食品與發(fā)酵工業(yè),2018,44(2):87-91.

LI Lu, HUANG Han-ying, LI Yi,et al.Establishment of a selection and detection model of fat in rice by near infrared spectrum characteristics[J].Food and Fermentation Industries,2018,44(2):87-91.

第一作者：博士，講師(黃漢英碩士為通訊作者，E-mail：hhywmx@mail.hzau.edu.cn)。

基金項目：中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(2662015PY078)；湖北省重大科技創(chuàng)新計劃(2014ABC009)

收稿日期：2017-06-13，改回日期：2017-06-21