基于高光譜成像技術(shù)的大曲還原糖含量預(yù)測及其可視化

劉亮,黃丹平*,田建平,黃丹,羅惠波,田穎,徐佳樂,葉建秋

1(四川輕化工大學(xué) 機械工程學(xué)院,四川 宜賓,644000)2(四川輕化工大學(xué) 生物工程學(xué)院,四川 宜賓,644000)

大曲主要以生料小麥為原料,通過自然網(wǎng)羅制曲環(huán)境中的微生物接種發(fā)酵,微生物在曲坯中此消彼長,自然積溫轉(zhuǎn)化并風干而成的一種多酶多菌的微生態(tài)制品[1-3]。大曲是白酒生產(chǎn)過程中的關(guān)鍵材料,大曲的好壞直接影響白酒品質(zhì)[4]。在發(fā)酵過程中,還原糖的含量是大曲質(zhì)量評價的重要指標。因此,實時檢測還原糖含量對生產(chǎn)優(yōu)質(zhì)大曲至關(guān)重要。大曲還原糖含量的常規(guī)檢測方法為化學(xué)方法,需要用化學(xué)試劑對樣品進行復(fù)雜而緩慢的實驗,且會對大曲樣品造成不可逆轉(zhuǎn)的破壞。因此,實現(xiàn)大曲還原糖含量的快速無損檢測具有重要意義。

高光譜成像技術(shù)作為一種新興的快速無損檢測技術(shù)[5]。它將機器視覺與光譜技術(shù)相結(jié)合,采集目標對象的二維圖像和一維光譜信息,從而獲取高分辨率的圖像和光譜數(shù)據(jù)[6]。區(qū)別于一般測量技術(shù)中對單點物質(zhì)含量進行測量的方法,利用高光譜成像技術(shù)獲得的光譜與圖像信息能夠?qū)崿F(xiàn)物質(zhì)含量變化的空間分析[7]。近年來,高光譜成像技術(shù)以其無損、綠色、高效等優(yōu)點逐步引起食品領(lǐng)域研究人員的關(guān)注[8-10],高光譜成像技術(shù)與化學(xué)計量學(xué)方法結(jié)合,被廣泛用于各種物質(zhì)含量檢測及可視化,如冬小麥氮素營養(yǎng)指標[11]、羊肉中飽和脂肪酸含量[12]、水稻籽粒直鏈淀粉含量[13]和馬鈴薯葉片葉綠素含量[14]。以上對物質(zhì)含量檢測方面的研究均取得了較好的結(jié)果,但在大曲物質(zhì)檢測方面,還沒有將高光譜無損技術(shù)應(yīng)用于大曲還原糖含量檢測的相關(guān)研究報道。

因此,本研究以濃香型大曲為研究對象,針對大曲還原糖含量這一重要指標,提出一種基于高光譜成像技術(shù)、圖像處理技術(shù)結(jié)合統(tǒng)計學(xué)方法對不同發(fā)酵時間的大曲還原糖含量進行檢測的方法。研究內(nèi)容如下：首先,采集大曲樣本對應(yīng)的高光譜圖像；使用標準正態(tài)變量校正(standard normal variables,SNV)、卷積平滑(smoothing convolution,SG)、多元散射校正(multiplicative scatter correction,MSC)對大曲樣本的高光譜數(shù)據(jù)進行預(yù)處理；采用主成分分析(principal component analysis,PCA)中的載荷系數(shù)法、連續(xù)投影法(successive projections algorithm,SPA)以及PCA-SPA算法分別提取了樣本的特征波長；建立了基于全波長和特征波長的偏最小二乘回歸(partial least squares regression,PLSR)和最小支持向量機(least squares support vector machine,LS-SVM)模型,預(yù)測大曲的還原糖含量；采用最優(yōu)預(yù)測模型對不同發(fā)酵時間的大曲樣本進行了還原糖含量可視化。

1 材料和方法

1.1 樣品制備

以四川宜賓某酒廠生產(chǎn)的大曲為樣本,該大曲發(fā)酵周期為22 d。為獲取不同發(fā)酵時期的大曲樣品,在發(fā)酵過程中每天從曲房預(yù)設(shè)位置分別取6個樣品,共獲得132個大曲樣品。采集的大曲樣品用于后續(xù)高光譜圖像采集與還原糖含量測定。

1.2 儀器與設(shè)備

對大曲近紅外高光譜數(shù)據(jù)進行采集的系統(tǒng)由芬蘭FX17系列高光譜相機、2個150 W鹵素光源、電控移動平臺、計算機組成。其中FX17系列相機光譜范圍900～1 700 nm,采樣間隔3.58 nm,共采集224個波段,采集方式為推掃式,分辨率640×181像素,采樣曝光時間4.02 ms,掃描速度16.57 mm/s。

1.3 高光譜圖像獲取

采集數(shù)據(jù)時將曲塊碾成曲粉,用勺子取曲粉樣品于石英器皿中,填充至與邊緣平齊的位置,再將石英器皿放置在電控移動平臺的中心位置進行數(shù)據(jù)采集,得到132組大曲樣本的原始高光譜數(shù)據(jù)。

提取的高光譜數(shù)據(jù)經(jīng)過黑白校正,才能轉(zhuǎn)化為表征每個樣品的反射率[15]。此外,黑白校正可以消除樣本形狀的差異以及高光譜相機中存在暗電流引起的噪聲[16]。校正公式如公式(1)所示：

(1)

式中：I,校正后反射光譜率；I0,原始高光譜數(shù)據(jù)；B,暗電流數(shù)據(jù)；W,標準白板數(shù)據(jù)。

校正后選取大曲樣本的感興趣區(qū)域(region of interest,ROI),將ROI的光譜數(shù)據(jù)進行平均處理后作為大曲樣本的原始高光譜數(shù)據(jù),原理如圖1所示。使用圖像處理方法檢測培養(yǎng)皿的邊緣,由邊緣得到培養(yǎng)皿中心坐標；并以中心坐標為圓心,繪制直徑為161像素的圓形掩膜；提取掩膜區(qū)域的光譜數(shù)據(jù),進行平均處理,得到該ROI的光譜反射率。每個樣本提取一個ROI,求得每個ROI內(nèi)所有像素點的平均光譜,得到132組光譜數(shù)據(jù)。

a-高光譜數(shù)據(jù)；b-提取的ROI；c-ROI內(nèi)每個像素點的光譜曲線；d-ROI平均光譜曲線

1.4 大曲還原糖測定

采用直接滴定法測定,此方法常用于食品中還原糖含量的測定。除去樣品中的蛋白質(zhì)后,以亞甲藍為指示劑,在加熱條件下滴定校準的堿性酒石酸銅溶液(用標準還原糖溶液校準),根據(jù)消耗量計算還原糖含量。省略了相關(guān)的溶液制備方法和具體的滴定過程,樣品的還原糖含量按公式(2)計算：

(2)

式中：X,樣品的還原糖含量,g/100 g；m1,還原糖的質(zhì)量,mg；m,樣品的質(zhì)量,g；F,系數(shù),1；10,試樣溶液的體積,mL；250,定容體積,mL；1 000為換算系數(shù)。

2 數(shù)據(jù)分析方法

2.1 光譜數(shù)據(jù)預(yù)處理

由于暗電流、環(huán)境光等因素的影響,采集的原始光譜中摻雜了與樣本無關(guān)的信息。采用合適的光譜預(yù)處理算法可以有效減弱各種因素對還原糖原始光譜數(shù)據(jù)的影響,提高預(yù)測模型的穩(wěn)定性和預(yù)測能力[17-18]。本文采用MSC、SG、SNV共3種預(yù)處理方法。

2.2 特征波長篩選

獲取的大曲高光譜樣本的數(shù)據(jù)量大,波段數(shù)量多,且相鄰波長圖像高度相關(guān),含有很多冗余信息和干擾信息,會影響建模的效率[19]。因此需要采用合適的方法篩選與表征指標高度相關(guān)的特征波長,以增加模型的魯棒性和泛化性[20]。本文采取PCA算法與SPA進行特征波長提取。

2.2.1 PCA

PCA是一種常見的數(shù)據(jù)分析方式,常用于高維數(shù)據(jù)的降維,可用于提取數(shù)據(jù)的主要特征分量。PCA得到的每個主成分都是原始數(shù)據(jù)中所有變量的線性組合。使用PCA載荷系數(shù)提取光譜特征波長步驟如下：首先使用PCA算法計算光譜數(shù)據(jù)的主成分；然后,選擇前幾個貢獻率大的主成分對應(yīng)的載荷因子；最后,根據(jù)閾值選擇載荷因曲線對應(yīng)的波長作為特征波長。具有較大載荷因子的波長對應(yīng)于光譜曲線中的重要峰谷[21]。本文使用PCA提取特征波長的目的是去除數(shù)據(jù)中不相關(guān)或信息含量低的波長,以降低數(shù)據(jù)維度。

2.2.2 SPA算法

SPA算法是一種使矢量空間共線性最小化的前向變量選擇算法,可以將有效的信息從大量的光譜數(shù)據(jù)中篩選出來,找到光譜變量之間共線性最小的特征波長[22]。本文使用SPA提取特征波長的目的是找到光譜變量之間共線性最小的特征波長,優(yōu)化建模條件。

2.2.3 PCA-SPA算法

本研究使用了一種新興的特征波長提取算法。該算法由PCA和SPA組成,稱為PCA-SPA算法[23]。使用該算法時先用PCA通過主成分的載荷系數(shù)提取光譜特征波長,然后在此基礎(chǔ)上使用SPA算法進行特征波長的二次提取。使用PCA-SPA算法的目的是確定還原糖的特征波長,構(gòu)建更穩(wěn)定、更簡單的檢測模型。

2.3 數(shù)學(xué)模型的建立與評價

2.4 還原糖含量可視化

為直觀測定指標的含量分布情況,為大曲樣品的ROI中還原糖的含量進行可視化處理。提取大曲樣本高光譜圖像每個像元的光譜反射率,依據(jù)檢測模型計算每個像元的還原糖含量。由于大曲還原糖含量在0.8～10 g/100 g,為了將不同發(fā)酵期間的大曲樣本感興趣區(qū)域內(nèi)各像素點的還原糖含量用偽彩色圖顯示,需要將各像素點的大曲還原糖含量數(shù)值映射到0～255灰度圖像的范圍期間。然后對灰度圖像進行偽彩色處理形成還原糖含量的可視化分布圖。還原糖含量的可視化彩色分布圖直觀的顯示出還原糖含量的二維分布,可以更好的觀察不同發(fā)酵時期的大曲還原糖含量及其分布情況。

3 結(jié)果與分析

3.1 大曲還原糖含量的變化

本研究采用Kennard-Stone(KS)算法將132個大曲光譜數(shù)據(jù)分為訓(xùn)練集(100)和預(yù)測集(32),如表1所示。

表1 大曲樣品還原糖含量統(tǒng)計表

3.2 大曲樣品的光譜特征

由于光譜曲線的首尾波長噪聲較大,為保證數(shù)據(jù)的準確性,去掉首尾10個波長的數(shù)據(jù)。圖2顯示了大曲在不同發(fā)酵時期的204個波段的平均光譜反射率曲線。

圖2 原始光譜圖

在波長范圍內(nèi),光譜曲線表現(xiàn)出相同的吸收帶和反射峰,但不同帶的反射率不同,這主要是由于大曲在不同發(fā)酵時期的物質(zhì)組成存在差異。近紅外波長范圍內(nèi)的吸收峰都是由O—H、C—H、N—H等共價鍵的吸收引起的。O—H作用于980 nm左右的波長,這是因為水的作用；C—H作用于1 210 nm附近的波長,這歸功于脂肪；N—H作用于1 442 nm附近的波長,這是由于蛋白質(zhì)[12]。

3.3 基于全波長光譜預(yù)測大曲還原糖的含量

基于原始光譜(無處理)和3種預(yù)處理的204個波長,分別建立PLSR和LS-SVM預(yù)測模型,建模效果如表2所示。

表2 PLSR和LS-SVM模型全光譜建模效果

圖3 SNV預(yù)處理結(jié)果

3.4 模型優(yōu)化

3.4.1 基于SPA提取特征波長的大曲還原糖預(yù)測模型

對SNV算法預(yù)處理后的204個波段的樣本光譜進行SPA特征波長提取,將提取后的特征波長進行PLSR與LS-SVM建模預(yù)測的結(jié)果如表3所示。

表3 基于SPA提取特征波長的光譜建模效果

3.4.2 基于PCA載荷因子提取特征波長的大曲還原糖預(yù)測模型

對SNV算法預(yù)處理后的204個波段的樣本光譜進行PCA運算,前3個主成分的貢獻率如圖4所示,其中還原糖的前3個主成分貢獻率之和達到99.96%,后續(xù)主成分的累計貢獻率提高幅度逐步變小,因此選擇前3個主成分的載荷因子提取特征波長。

圖4 PCA中主成分的貢獻率

利用PCA光譜數(shù)據(jù)主成分,根據(jù)前3個主成分對應(yīng)的載荷因子提取還原糖的特征波長,根據(jù)設(shè)置的載荷因子的閾值,可得到不同數(shù)目的特征波長,并進行PLSR和LS-SVM建模,建模效果如表4所示。隨著載荷閾值的增加,提取的特征波長數(shù)量會快速的減少,而后減少得較慢,最后不發(fā)生變化,這表明PCA載荷因子對變量的選擇是一個從粗略到精細的過程。當載荷因子的閾值達到0.060時,選擇的特征波長數(shù)目變化趨于平緩,且此時建模效果較差,預(yù)測精度皆在0.86及以下,說明閾值設(shè)置過高,在提取特征波長的過程中刪去了一部分與還原糖相關(guān)的波段。說明當載荷因子的閾值小于0.060時,算法濾除與大曲還原糖無關(guān)的光譜信息。

表4 基于PCA提取特征波長的光譜PLSR建模效果

3.4.3 基于PCA-SPA提取特征波長的大曲還原糖預(yù)測模型

在研究中,分別將大曲的光譜數(shù)據(jù)進行了PCA和SPA算法的特征波長提取。單獨使用SPA算法時,盡管降低了特征波長點建的共線性,但也會保留下一些無信息、甚至?xí)a(chǎn)生干擾的波長點。同樣,當在單獨使用PCA算法進行特征波長提取時,也會留下一些信息量較少的波長。因此,本實驗中將PCA算法與SPA算法聯(lián)合,實現(xiàn)兩者之間的優(yōu)勢互補,建立基于PCA-SPA特征波段選取的PLSR大曲還原糖含量預(yù)測模型。

使用PCA-SPA算法對光譜數(shù)據(jù)進行特征波長提取。設(shè)定PCA的載荷因子閾值為0.045,經(jīng)過2種算法提取特征,得到26個特征波長點,分布如圖5。圖中藍色光譜曲線為132個大曲樣本的平均光譜反射率曲線,橙色標記點對應(yīng)選擇的波段位置。還原糖的特征波長主要分布在短波近紅外區(qū),且基本分布在光譜曲線的波峰波谷附近。

圖5 PCA-SPA算法提取的特征波長

對基于SNV-PCA-SPA方法提取的特征波長進行PLSR與LS-SVM建模分析。LS-SVM預(yù)測模型的訓(xùn)練集和預(yù)測集的決定系數(shù)分別為0.863 4和0.802 1,訓(xùn)練集和預(yù)測集的均方根誤差分別為0.662 2和1.119 4 g/100 g。PLSR模型預(yù)測大曲還原糖的效果最好,訓(xùn)練集和預(yù)測集的決定系數(shù)分別為0.983 2和0.922 7,訓(xùn)練集和預(yù)測集的均方根誤差分別為0.057 3和0.455 6 g/100 g,光譜波段數(shù)減少了87%。PLSR模型預(yù)測集的擬合效果如圖6所示。

圖6 預(yù)測集的擬合效果圖

3.4.4 最優(yōu)大曲還原糖預(yù)測模型

經(jīng)過上述實驗分析,本研究確定了SNV-PCA-SPA-PLSR為最優(yōu)的大曲還原糖含量的預(yù)測模型。首先對基于原始光譜和3種預(yù)處理的204個波長,建立PLSR和LS-SVM預(yù)測模型,得出經(jīng)SNV預(yù)處理算法的PLSR建模效果最佳；為減少建模的復(fù)雜度以及提升建模精度,分別使用PCA載荷系數(shù)法以及SPA算法對光譜數(shù)據(jù)進行特征提取,經(jīng)建模分析,得出最佳模型為PCA-PLSR以及SPA-PLSR,建模所用光譜數(shù)據(jù)分別減少70%和44%；為進一步提高建模效果,提出PCA-SPA算法進行特征提取,經(jīng)建模驗證,PLSR模型精度及準確度皆有一定的提高,建模所用光譜數(shù)據(jù)減少87%。綜上所述,本研究的最優(yōu)大曲還原糖含量的預(yù)測模型為SNV-PCA-SPA-PLSR模型。

3.5 大曲還原糖含量可視化

在確定了SNV-PCA-SPA-PLSR為最優(yōu)的大曲還原糖含量的預(yù)測模型后,使用最優(yōu)模型對大曲還原糖含量進行可視化分析。不同發(fā)酵時期的大曲樣本的還原糖含量可視化如圖7所示,其中不同的顏色代表不同的還原糖含量值,深紅色區(qū)域表示還原糖含量值較高,深藍色區(qū)域表示還原糖含量值較低。隨著還原糖含量的增加,顏色逐漸由深藍色變?yōu)樯罴t色。ROI中分布不均還原糖的含量可以通過顏色清晰區(qū)分,與單點檢測技術(shù)相比,高光譜成像可以實現(xiàn)大曲任意局部區(qū)域的可視化。

a-第1天；b-第5天；c-第9天；d-第13天；e-第17天；f-第22天

由圖7可知,不同發(fā)酵階段的大曲還原糖含量明顯不同,隨著發(fā)酵時間的進行,還原糖含量不斷降低,顏色逐漸由紅變藍。發(fā)酵前期有較多紅色區(qū)域,因為發(fā)酵剛開始時,微生物大量繁殖,產(chǎn)出大量還原糖；隨著發(fā)酵進行,還原糖含量逐漸下降,這是因為酵母菌等其他菌種的快速增殖消耗了部分還原糖[25]。高光譜成像可以直觀顯示大曲還原糖分布,為判定大曲發(fā)酵狀態(tài),調(diào)節(jié)曲房環(huán)境提供了依據(jù)。

4 結(jié)論

在進一步研究中,應(yīng)擴大樣本獲取范圍,采集不同發(fā)酵批次的大曲樣本,擴大樣本數(shù)量,增大樣本的還原糖范圍,使模型的預(yù)測性更高,應(yīng)用更加廣泛。