王國梁， 余克強， 成 鍇， 劉 鑫， 王文俊， 李 洪， 郭二虎， 李志偉*

1. 山西農(nóng)業(yè)大學農(nóng)業(yè)工程學院， 山西 太谷 030801 2. 山西農(nóng)業(yè)大學谷子研究所， 山西 長治 046000 3. 西北農(nóng)林科技大學機械與電子工程學院， 陜西 楊凌 712100

引 言

隨著現(xiàn)代人對生活品質的提升以及健康產(chǎn)業(yè)興起， 小米產(chǎn)業(yè)作為山西省“有機旱作”農(nóng)業(yè)發(fā)展的龍頭產(chǎn)業(yè)， 其深加工產(chǎn)品市場價值逐漸被大眾所普遍肯定， 而小米米粉又以其獨特的食味品質優(yōu)勢， 現(xiàn)已成為嬰幼兒輔食佳品。 小米米粉中主要成分是淀粉， 淀粉在一定條件下具有晶體性不可逆的特點， 在加入足量水情況下， 將其混合液加熱、 升溫， 淀粉吸水膨脹、 溶解呈現(xiàn)出這種特性， 而這一現(xiàn)象稱為糊化現(xiàn)象。 其中堿消值是能夠直接反應糊化特性的主要特征指標， 堿消值越低， 糊化溫度越高， 直鏈淀粉含量越高， 而小米米粉口感粘糯性變差[1-2]， 最終影響小米產(chǎn)品的市場價格。 通常利用人工評級、 冷堿糊化法和粘度速測儀(rapid visco analyser, RVA)測定淀粉粘度特征值來評價小米米粉的糊化特性及品質[3]； 而人工評級人為因素大， 結果不可靠； 冷堿糊化法和RVA粘度測定耗時， 破壞樣品理化特性， 并且制備樣品操作程序繁瑣， 成本太高， 因此， 樣品批量處理和實現(xiàn)快速檢測存在一定困難。

高光譜技術是一種通過精確采集被測物每個像素點的連續(xù)光譜信息， 并對被測樣品的感興趣區(qū)域(region of interest， ROI)數(shù)據(jù)提取和處理， 從而反映被測物組成成分和分子結構的無損檢測技術， 具有檢測效率高， 不破壞被測樣品的理化特性， 檢測時間短并具有線上實時分析等特點。 目前高光譜技術已廣泛被應用在土壤和農(nóng)產(chǎn)品的分類、 成分反演等農(nóng)業(yè)檢測相關工作中[3-14]。 Shao等利用高光譜成像結合化學計量學方法對不同種類的米粉進行分類、 判別， 運用競爭性自適應重加權采樣法(competitive adaptive reweighted sampling， CARS)建立最小二乘支持向量機回歸模型， 對品種分類的結果最好[4]； Chen等采用連續(xù)投影算法結合多元線性回歸對谷子蛋白質、 糖類和粗脂肪含量進行近紅外光譜特性分析， 結果表明利用谷子近紅外光譜數(shù)據(jù)能夠反演谷子中蛋白質和糖類含量[5]； Chen等運用隨機蛙跳(random frog， RF)算法提取特征波段結合高斯過程回歸預測木材含水模型效果最好[6]。 當前， 小米米粉的高光譜研究主要集中在品種分類和成分反演等方面[4-5]， 而針對小米米粉糊化特性的研究相關報道仍然較少。

以晉東南地區(qū)“羊肥小米”主產(chǎn)地采集到的小米米粉作為實驗樣本， 運用高光譜儀器采集小米米粉光譜數(shù)據(jù)， RVA儀器測定小米米粉堿消值指標， 結合化學計量方法， 建立快速檢測小米米粉堿消值預測模型， 旨在尋求一種快速、 無損、 低成本預測小米米粉堿消值的方法。 通過探索小米米粉糊化特性與高光譜之間的聯(lián)系， 間接反映小米米粉中直鏈淀粉占比， 為小米米粉品質評級及加工提供參考依據(jù)， 進而為小米米粉堿消值傳感器的開發(fā)提供理論支撐。

1 實驗部分

1.1 樣本制備

樣本材料“羊肥小米”收獲于山西省長治市武鄉(xiāng)縣。 在大田收獲季采集實驗樣本材料過程中， 采用棋盤式取樣法取樣， 取樣點數(shù)1 800個， 單點采樣面積2 m2， 將相鄰5個采樣點采集到的材料混勻， 得到實驗樣本358份經(jīng)清選、 自然晾曬、 去殼、 研磨、 過80目篩， 取得待測小米米粉樣品， 并采用四分法取樣[7]， 被測樣品先用于小米米粉高光譜數(shù)據(jù)采集， 然后進行實驗室堿消值測定。

1.2 儀器與數(shù)據(jù)采集

采用美國Headwall Photonics公司的高光譜成像儀采集光譜數(shù)據(jù)。 該設備工作組件主要包括： 高光譜成像儀(光圈1.4， 焦距25 mm)、 鎢絲燈、 可移動式掃描平臺、 電腦等。 數(shù)據(jù)采集設備參數(shù)為： 光譜通道步長4.715 nm， 采集范圍900～1 700 nm， 共包含170個波段。 為能采集到清晰圖像， 調(diào)整平臺移動速度為16 mm·s-1， 曝光時間0.9 ms， 被測樣品上表面距鏡頭300 mm。

實驗前， 采集暗背景(反射率0.1%)圖像B， 掃描白背景(反射率99.9%)圖像W， 采集樣品的絕對圖像I， 根據(jù)式(1)校正為相對圖像X[8]。 每采集3幅高光譜圖像， 重復一次光譜校正。

(1)

進行高光譜數(shù)據(jù)采集時， 將被測樣品裝入直徑5 cm， 高度1 cm的圓形試驗器皿內(nèi)， 采樣過程中要保證樣品表面處理平整、 緊實， 將被測樣品置于推掃平臺上， 每份樣品推掃三次， 采集光譜數(shù)據(jù)， 同一樣品3次推掃， 取均值、 編號， 按實驗時序保存。

1.3 小米米粉堿消值測定

采用澳大利亞Newport公司生產(chǎn)粘度分析儀， 根據(jù)GB/T 24852—2010方法測得小米米粉堿消值[9]。

1.4 小米米粉高光譜數(shù)據(jù)處理方法

光譜數(shù)據(jù)采集時， 在被測樣品感興趣區(qū)域ROI， 按像素點逐一選擇， 提取光譜數(shù)據(jù)矩陣， 并進行均值運算， 得到每個樣品在各個光譜波段的平均光譜值。 利用Kennard-Stone算法(K-S)將光譜數(shù)據(jù)按照2∶1分為訓練集和預測集[10-11]。 對光譜數(shù)據(jù)全波段、 CARS[4，12]和RF[6，13]算法選擇特征波段后建立偏最小二乘回歸(partial least square regression， PLSR)預測模型， 模型評估利用相關系數(shù)(correlation coefficient，R)值。 為消除噪聲對光譜數(shù)據(jù)的影響， 提高模型預測精度， 進一步對原始光譜數(shù)據(jù)分別運用S-G， 多元散射校正(multiplicative scatter correction， MSC)和S-G+MSC進行預處理[14]， 再利用CARS和RF法選擇特征波段， 建立PLSR模型。 光譜數(shù)據(jù)處理軟件主要有The Unscrambler X 10.1(CAMO， Norway)和MATLAB 2018a(The MathWorks， USA)。

2 結果與討論

2.1 小米米粉堿消值統(tǒng)計分析

稱重(3±0.01) g小米米粉樣品(按12%濕基校準)與(25±0.1) mL水(按12%濕基校準)在鋁盒中混勻， 置于RVA分析儀中進行堿消值檢測， 3次重復測量堿消值， 取每份樣品3次測量均值為該實驗最終值。 小米米粉堿消值統(tǒng)計結果見表1。

表1 小米米粉堿消值統(tǒng)計結果

2.2 小米米粉光譜特性

圖1為小米米粉平均光譜曲線， 由于靠近光譜量程首尾兩端(≤900 nm和≥1 700 nm)， 光譜反射率噪聲較大， 故截取光譜可用范圍為950～1 650 nm， 光譜數(shù)降為148條， 在每幅圖像中的ROI提取光譜數(shù)據(jù)， 經(jīng)均值運算后組成358×148待數(shù)據(jù)處理的光譜數(shù)據(jù)矩陣。 從圖1可看出， 曲線吸收峰值主要集中在980， 1 210和1 450 nm附近， 其中980和1 210 nm處主要受被測樣品淀粉影響， 1 450 nm處受被測樣品中水分子吸收影響[12]。 利用小米米粉全波段繪制成的光譜曲線變化平滑， 但會存在數(shù)據(jù)重疊、 平移和微小噪聲的問題。 為提高精度， 對特征波段提取后的數(shù)據(jù)矩陣進行預處理。

圖1 小米米粉平均光譜曲線

2.3 特征波段的選擇

2.3.1 CARS特征波段選擇

CARS是借鑒達爾文進化論“優(yōu)勝劣汰”的理論思想， 數(shù)據(jù)處理采用該算法進行光譜特征波段快速選擇[10]。 圖2為CARS選擇特征波段的過程和主要參數(shù)的變化曲線， 其中圖2(a)表示被測樣品數(shù)隨采樣次數(shù)增加的變化趨勢； 圖2(b)中可以看出通過交叉檢驗(cross validation， CV)得到隨采樣次數(shù)增加均方根誤差(root mean square error cross validation， RMSECV)的變化曲線， 該曲線在1—26平滑下降， 27—50逐步階梯狀回升， 最小值出現(xiàn)在26， 此點處RMSECV達到最小值； 從圖2(c)中回歸系數(shù)(regression coefficients， RC)， 變量最優(yōu)選擇出現(xiàn)在圖2(c)中藍色星號豎線標注的位置， 在這點處盡可能多的保留了變量信息并選擇出特征波段。 經(jīng)過CARS特征波段選擇， 共有16條特征波段被選出， 分別為964， 1 006， 1 053， 1 166， 1 171， 1 195， 1 213， 1 218， 1 223， 1 246， 1 298， 1 308， 1 350， 1 378， 1 444和1 529 nm。 其中964， 1 006和1 053 nm在O—H三級倍頻附近(980 nm)， 1 166， 1 171， 1 195， 1 213， 1 218， 1 223， 1 246， 1 298， 1 308， 1 350和1 378 nm在C—H二級倍頻及組合頻附近(1 220， 1 360 nm)， 1 444 nm此處O—H鍵與樣品中水分子的吸收峰重合[12](1 450 nm)與淀粉的近紅外特性相應。

圖2 利用CARS算法的特征變量選擇

2.3.2 RF特征波段選擇

RF是以統(tǒng)計學思想為理論基礎， 通過重復計算各變量被選概率來評價變量， 具有不易過擬合優(yōu)點[11]。 RF算法運行前， 先設置相關參數(shù)(迭代參數(shù)、 運行次數(shù)、 初始跳轉變量、 候選變量集被選概率等)， 依據(jù)相關文獻[13]選擇0.4作為特征變量閾值， 運行次數(shù)50次， 將運算后的均值作為變量評價的標準， 每個波長被選擇的可能性(selection probability， SP)如圖3所示。 在SP曲線中， 可以發(fā)現(xiàn)有些特征值具有極高的SP值， 這些峰值與小米米粉堿消值有極高的相關性。 最終選出10條特征波段， 分別為964， 1 053， 1 171， 1 195， 1 218， 1 223， 1 242， 1 270， 1 275和1 279 nm。 采用RF與CARS選出特征波段大部分重疊， 這些特征波段與C—H和O—H鍵的倍頻吸收有關系， 并與淀粉的近紅外特性相應。 執(zhí)行CARS和RF運行結果雖均為隨機， 但為保證波段選擇的有效性， 兩種算法通過多次運算(均設定50次)， 從而達到減少隨機因素的影響[10]。

圖3 利用RF算法的光譜波段被選概率

2.4 小米米粉堿消值的PLSR回歸模型建立

表2列出采用全波段(RAW)和特征波段(CARS、 RF提取)建立的小米米粉堿消值偏最小二乘回歸模型。 從結果看出用RAW建立的PLSR，Rp值最高為0.77； 采用CARS、 RF提取特征波段建立的PLSR，Rp值分別為0.72和0.7，Rp值降低； 究其原因在于這兩種方法只提取到了部分與O—H鍵和C—H鍵光譜特性有關的光譜數(shù)據(jù)[12-13]， 而導致其他相關數(shù)據(jù)信息的丟棄。 通過采用CARS、 RF提取特征波段建立的PLSR，Rp值與運用RAW建立的PLSR接近， 這反映了采用CARS和RF建模具有可行性， 可以對小米米粉高光譜反射率進行堿消值的預測。

表2 不同特征波段提取方法PLSR模型結果

為了提高模型預測精度， 基于表2運算結果， 對光譜數(shù)據(jù)預處理， 從表3中可看出， MSC處理后光譜數(shù)據(jù)[14]Rp值達到0.83， 精度最高。 因此， 對原始光譜數(shù)據(jù)采用MSC預處理所得到的回歸模型， 數(shù)據(jù)集的擬合程度較高， 訓練集與校正集的擬合情況如圖4所示。

表3 不同預處理方法對原始數(shù)據(jù)PLSR模型分析結果

圖4 采用MSC預處理下的訓練集和預測集擬合情況

對MSC預處理的光譜數(shù)據(jù)分別進行二次CARS和RF特征波段提取， 從表4中可看出， 兩種算法各選擇10條特征波段， 前者選出964， 978， 992， 1 053， 1 166， 1 218， 1 246， 1 270， 1 275和1 553 nm， 后者選出978， 992， 1 053， 1 171， 1 195， 1 218， 1 242， 1 246， 1 270和1 275 nm， 從以上數(shù)據(jù)結果可以看出經(jīng)MSC預處理后提取的特征波段數(shù)量減少， 但與未進行預處理的回歸模型相比，Rp值變化不大， 這也說明CARS和RF具有一定的穩(wěn)定性， 可以作為小米米粉高光譜反射率預測堿消值的參考方法。

表4 采用MSC預處理下的不同特征波長提取方法PLSR預測模型結果

從以上數(shù)據(jù)處理結果看， 在不破壞小米米粉內(nèi)部理化特性的情況下， 采用MSC預處理光譜數(shù)據(jù)， 建立快速檢測小米米粉堿消值預測模型， 能夠反演小米米粉堿消值。 從而反映小米米粉中直鏈淀粉的占比， 進而為小米米粉品質評級、 加工及小米米粉堿消值傳感器的開發(fā)提供參考依據(jù)和數(shù)據(jù)支撐。

3 結 論

采用高光譜技術可以實現(xiàn)小米米粉堿消值快速檢測， 進一步為小米米粉品質評級及加工提供參考依據(jù)， 主要結論為：

(1)采用全光譜波段、 CARS和RF提取特征波段建立的小米米粉堿消值偏最小二乘回歸模型， 結果表明三者Rp值較為相近， 全波段最高Rp為0.77， 同時采用CARS和RF建模具有可行性， 能夠利用小米米粉高光譜反射率反演其堿消值。

(2)采用MSC預處理原始數(shù)據(jù)并作出PLSR模型分析，Rp值達到0.83。 對經(jīng)MSC預處理的數(shù)據(jù)集， 再次作CARS和RF特征波段提取， 建立回歸模型。 實驗結果能夠為小米米粉品質評級、 加工及堿消值傳感器的開發(fā)提供參考依據(jù)和數(shù)據(jù)支撐。