劉燕德，饒 宇，孫旭東，姜小剛，徐 海，李 雄，王觀田，徐 佳

華東交通大學(xué)機(jī)電與車輛工程學(xué)院，江西 南昌 330013

引 言

糖度(SSC)是影響水果口感的一個(gè)重要因素，SSC差異在2°Brix以上，水果口感差異明顯[1]。近年來(lái)，人們對(duì)口感品質(zhì)十分重視，傳統(tǒng)外部品質(zhì)分選難以滿足口感品質(zhì)分選需求，可見/近紅外光譜技術(shù)作為一種快速無(wú)損檢測(cè)手段，目前已經(jīng)廣泛應(yīng)用于蘋果[2]、梨、桃[3]等水果的內(nèi)部品質(zhì)檢測(cè)。

崔豐娟[4]等研究了蘋果動(dòng)態(tài)在線檢測(cè)的傳輸速度對(duì)模型的精度的影響。曹松濤在他的學(xué)位論文中報(bào)道了對(duì)梨樣品相關(guān)因素對(duì)內(nèi)部品質(zhì)普適模型的穩(wěn)健性影響，樣品的尺寸差異、產(chǎn)地差異、品種差異都是主要影響，通過(guò)對(duì)光譜進(jìn)行一階、二階微分預(yù)處理，基本能夠消除差異影響，提高了模型的預(yù)測(cè)精度。Saranwong[5]等對(duì)比分析了兩種不同近紅外漫反射儀器在芒果品質(zhì)檢測(cè)上的應(yīng)用。Greensill[6]等使用自行研發(fā)的光學(xué)系統(tǒng)，研究了硬皮甜瓜果實(shí)不同果肉位置的SSC取樣準(zhǔn)確度，結(jié)果表明，果實(shí)赤道部位SSC取樣值與光譜建模模型準(zhǔn)確度最高，能代表水果的真實(shí)糖度值。Mizuki[7]等結(jié)合偏最小二乘判別分析，應(yīng)用蘋果分選機(jī)頂部和底部的光譜儀，建立了蘋果內(nèi)部褐變發(fā)生的判別模型。

本工作主要研究了臍橙的尺寸差異對(duì)可見近紅外漫透射及兩種不同漫反射檢測(cè)方式模型的影響，對(duì)比分析了大果、小果及混合果模型的預(yù)測(cè)精度，考察了近紅外光譜模型及多元散射校正預(yù)處理方法對(duì)消除尺寸差異的效果。

1 實(shí)驗(yàn)部分

1.1 材料

實(shí)驗(yàn)所用的臍橙樣品購(gòu)于南昌市某水果市場(chǎng)，將臍橙表面擦凈，在25 ℃的室溫下靜置24個(gè)小時(shí)，依次采集臍橙樣品的橫縱徑、重量、糖度的理化指標(biāo)如表1所示。實(shí)驗(yàn)共采用480個(gè)樣品，按照臍橙的橫縱徑指標(biāo)劃分大小果，其中橫徑范圍在75～90 mm，縱徑在80～95 mm為大果共312個(gè)，橫徑范圍在60～75 mm，縱徑在65～85 mm為小果共168個(gè)，大果建模集為234個(gè)，預(yù)測(cè)集為78個(gè)，小果建模集為126個(gè)，預(yù)測(cè)集42個(gè)，混合果建模集360個(gè)，預(yù)測(cè)集120個(gè)。

表1 臍橙的理化指標(biāo)范圍統(tǒng)計(jì)結(jié)果Table 1 Statistical results of physical andchemical indicators of navel orange

1.2 裝置和光譜采集

可見近紅外漫反射光譜采集平臺(tái)如圖1所示，實(shí)驗(yàn)平臺(tái)搭建在重量分選裝置上的出果口前，考慮實(shí)際應(yīng)用場(chǎng)景，光源發(fā)射與接收探頭安裝在距離果杯15 cm處。采用Ocean Optics公司的QE PRO光譜儀作為實(shí)驗(yàn)的光譜采集裝置，其波長(zhǎng)范圍為350～1 150 nm，光源為1個(gè)12 V、100 W的鹵鎢燈，采用穩(wěn)壓電源供電。光纖探頭分別為環(huán)形發(fā)射與接收和多點(diǎn)發(fā)射與接收漫反射式光纖探頭，均設(shè)計(jì)為將發(fā)射光源與接收光源裝置集成在同一探頭上，通過(guò)鹵鎢燈提供直射光源，從光源接收端口進(jìn)入光纖，由探頭發(fā)射端發(fā)射出光照，光信號(hào)經(jīng)過(guò)樣品內(nèi)部后，經(jīng)過(guò)吸收、反射等，由探頭接收端口接收到漫反射光信號(hào)，經(jīng)過(guò)光纖傳輸并觸發(fā)光譜儀進(jìn)行數(shù)據(jù)保存，每觸發(fā)一次，光譜采集軟件便儲(chǔ)存一條光譜。

圖1 近紅外漫反射在線檢測(cè)裝置1： 光源； 2： 光源接收端口； 3： 光源發(fā)射與接收端口；4： QE PRO光譜儀； 5： 電腦； 6： 臍橙樣品； 7： 果杯Fig.1 Near-infrared diffuse reflection on-line detecting device

開機(jī)預(yù)熱0.5 h，聚四氟乙烯球作為光譜參比，采集臍橙樣品光譜，實(shí)驗(yàn)參數(shù)為： 積分時(shí)間100 ms，速度5個(gè)·s-1，光強(qiáng)100 W。光譜采集軟件采集光譜，分別標(biāo)記面光譜，取平均光譜作為該樣品可見近紅外漫反射光譜。

漫透射光譜采集裝置如圖2所示，采用海洋光學(xué)公司的QE65000光譜儀，波長(zhǎng)范圍為350～1 150 nm，光源為10個(gè)12 V、100 W的鹵鎢燈，均勻分布在傳送帶兩側(cè)，通過(guò)穩(wěn)壓電源供電。光纖探頭位于果杯正下方，光信號(hào)從樣品表面進(jìn)入，經(jīng)過(guò)透射到達(dá)樣品下方的光纖探頭，探頭每接收到光信號(hào)一次，光譜儀觸發(fā)一次，則光譜采集軟件保存一條光譜。果杯設(shè)計(jì)有彈性避光圈與透光孔，避免雜散光與外界干擾光被光纖接收到。

開機(jī)預(yù)熱0.5 h，聚四氟乙烯球作為光譜參比，采集臍橙樣品光譜，實(shí)驗(yàn)參數(shù)為： 積分時(shí)間100 ms，速度5個(gè)·s-1，光強(qiáng)1 000 W。采用光譜采集軟件采集光譜，分別標(biāo)記面光譜，取平均光譜作為該樣品可見近紅外漫反射光譜。

圖2 近紅外漫透射在線檢測(cè)裝置1： 光源； 2： 樣品； 3： 遮光圈；4： 光路； 5： 光纖探頭； 6： 彈跳爪Fig.2 Near-infrared diffuse transmissionon-line detecting device1： Light source； 2： Sample； 3： Shading ring；4： Optical path； 5： Fiber optic probe； 6： Bouncing claw

1.3 糖度含量判定

臍橙糖度含量由折射式數(shù)字糖度計(jì)(ATAGO PAL-α，日本ATAGO公司)測(cè)定，測(cè)量前，用純凈水洗凈并擦干糖度計(jì)，并進(jìn)行標(biāo)定，糖度值為0°Birx，測(cè)量時(shí)，使用一次性紗布將標(biāo)記部分約5 mm深果肉濾出少量果汁，滴于糖度計(jì)的測(cè)試窗口，重復(fù)3次，取3次的平均值作為糖度真實(shí)值，所得值以質(zhì)量分?jǐn)?shù)表示。

2.4 數(shù)據(jù)處理

采用matlab2012a將光譜數(shù)據(jù)轉(zhuǎn)換為數(shù)據(jù)矩陣，再通過(guò)Unscrambler進(jìn)行PLS建模。建模前先進(jìn)行光譜數(shù)據(jù)預(yù)處理，再使用偏最小二乘法(PLS)分別建立漫透射及兩種不同漫反射式的糖度模型。偏最小二乘法中的相關(guān)系數(shù)r和預(yù)測(cè)集均方根誤差RMSEP、建模集均方根誤差RMSEC計(jì)算公式如式(1)—式(3)所示

(1)

(2)

(3)

2 結(jié)果與討論

2.1 不同檢測(cè)方式的尺寸差異影響分析

近紅外光譜分析是一種應(yīng)用于水果內(nèi)部品質(zhì)無(wú)損檢測(cè)的重要手段，但由于樣品的相關(guān)因素如樣品溫度、大小、產(chǎn)地、品種等的影響，模型的穩(wěn)健性受到較大影響，其中，樣品大小是一個(gè)重要的影響因素。近紅外漫透射和漫反射檢測(cè)方式不同大小樣品光路圖如圖3所示，其中圖3(a)為漫透射檢測(cè)大、小尺寸的臍橙光路圖，在相同光照路徑下，光在大果中的光程要比小果中的光程更長(zhǎng)，光的能量衰減嚴(yán)重，光纖探頭接收到的能量強(qiáng)度更低，而小果由于它的小尺寸，光能快速透過(guò)，光纖探頭接收到的能量衰減少，強(qiáng)度強(qiáng)，攜帶的有效信息多。圖3(b)為漫反射檢測(cè)大、小尺寸的臍橙光路圖，在相同光照路徑下，光在大果中的光程要比小果中的光程要短，可能穿過(guò)的果肉組織不多，攜帶的有效信息少，而在小果中，光程較長(zhǎng)，穿過(guò)的果肉組織較多，攜帶的有效信息也多。

圖3 漫透射(a)和漫反射(b)上大、小尺寸的臍橙光路示意圖

2.2 不同檢測(cè)方式的光譜對(duì)比分析

對(duì)樣品的尺寸進(jìn)行測(cè)量，每個(gè)樣品的最大橫徑測(cè)量5次，取平均值作為樣品的橫徑指標(biāo)，考察尺寸差異對(duì)臍橙糖度在線檢測(cè)精度的影響。由圖4可知，無(wú)論是漫透射檢測(cè)還是漫反射檢測(cè)，小果的光譜能量強(qiáng)度都要比大果強(qiáng)。對(duì)同一樣品的三種不同檢測(cè)方式的光譜進(jìn)行對(duì)比分析。兩種漫反射光譜的能量要比漫透射光譜強(qiáng)，而環(huán)形發(fā)射與接收漫反射光譜又要比多點(diǎn)發(fā)射與接收漫反射光譜的能量強(qiáng)。漫反射光譜在690 nm左右存在吸收峰，漫透射光譜的峰在700和825 nm左右，從兩種檢測(cè)方式的光譜波形上來(lái)看，漫透射光譜要比漫反射光譜多一個(gè)吸收峰，且漫透射光譜集中在600～850 nm范圍內(nèi)，漫反射光譜范圍則在500～900 nm內(nèi)，比漫透射光譜范圍大一點(diǎn)。對(duì)比分析三種不同檢測(cè)方式，漫透射檢測(cè)方式中，光源位于樣品兩側(cè)，接收裝置位于樣品下方，果杯上的橡膠遮光圈能有效消除樣品表面反射光及其他干擾光進(jìn)入接收器； 在漫反射檢測(cè)方式中，光源與接收裝置集成在同一光纖上，環(huán)形發(fā)射與接收的漫反射檢測(cè)方式中，由于它的環(huán)形光源的結(jié)構(gòu)特點(diǎn)，光在樣品表面形成一個(gè)圓形光斑，可以有效的將外界干擾光及樣品表面反射光隔離在光斑外面，達(dá)到克服雜散光的目的； 多點(diǎn)發(fā)射與接收漫反射檢測(cè)方式中，由于它是點(diǎn)光源，點(diǎn)接收方式，所以它克服樣品表面反射光和其他干擾光的能力要比前兩種檢測(cè)方式稍差。

圖4 三種檢測(cè)方式的大小果光譜對(duì)比Fig.4 Spectral comparison of fruit with differentsizes using three detection methods

2.3 不同檢測(cè)方式下尺寸預(yù)測(cè)模型建立

不同樣品尺寸在某一波長(zhǎng)范圍光下具有一定的衰減規(guī)律，實(shí)驗(yàn)選用波長(zhǎng)范圍為900～936 nm建立近紅外光譜尺寸預(yù)測(cè)模型?？梢?近紅外光譜波長(zhǎng)范圍320～500 nm是可見光波段，受顏色變化影響較大； 500～900 nm是光譜的有效信息集中波段，用來(lái)建立糖度預(yù)測(cè)模型； 936～1 150 nm波段噪聲和無(wú)效信息較多，故選用900～936 nm建立尺寸預(yù)測(cè)模型。不同檢測(cè)方式下，尺寸模型預(yù)測(cè)散點(diǎn)圖如圖5所示。由圖5可知，漫透射檢測(cè)方式下尺寸模型的預(yù)測(cè)集相關(guān)系數(shù)為0.60，預(yù)測(cè)集均方根誤差3.95°Brix，截距B=81.27； 多點(diǎn)發(fā)射與接收漫反射檢測(cè)方式下尺寸模型的預(yù)測(cè)集相關(guān)系數(shù)為0.97，預(yù)測(cè)集均方根誤差1.46°Brix，截距B=94.3； 環(huán)形發(fā)射與接收漫反射檢測(cè)方式下尺寸模型的預(yù)測(cè)集相關(guān)系數(shù)為0.96，預(yù)測(cè)集均方根誤差1.73°Brix，截距B=70.73。由不同的檢測(cè)方式下尺寸預(yù)測(cè)散點(diǎn)圖可知，漫透射檢測(cè)方式的尺寸模型預(yù)測(cè)精度較低，兩種漫反射檢測(cè)方式的尺寸模型預(yù)測(cè)精度較高，可以有效的預(yù)測(cè)樣品的尺寸大小。

圖5 不同檢測(cè)方式下尺寸模型預(yù)測(cè)結(jié)果散點(diǎn)圖(a)： 漫透射檢測(cè)方式尺寸預(yù)測(cè)散點(diǎn)圖；(b)： 多點(diǎn)發(fā)射與接收漫反射檢測(cè)方式尺寸預(yù)測(cè)散點(diǎn)圖；(c)： 環(huán)形發(fā)射與接收漫反射檢測(cè)方式尺寸預(yù)測(cè)散點(diǎn)圖Fig.5 Scatter plot of dimensional model predictionresults under different detection modes

2.4 糖度預(yù)測(cè)模型建立

王動(dòng)民(本刊34卷9期)利用MSC方法對(duì)光譜進(jìn)行預(yù)處理，基于不同波段的光譜數(shù)據(jù)會(huì)得到不同的校正結(jié)果，進(jìn)而影響所得定標(biāo)模型的可靠性。分別建立三種不同檢測(cè)方式的不同果型的PLS模型，模型預(yù)測(cè)結(jié)果如表2所示，通過(guò)預(yù)測(cè)集相關(guān)系數(shù)Rp來(lái)評(píng)價(jià)模型的預(yù)測(cè)能力。由表2可知，三種不同的近紅外光譜檢測(cè)方式下，小果的PLS模型預(yù)測(cè)能力要強(qiáng)于大果和混合果，其中漫透射方式的Rp為0.76，多點(diǎn)發(fā)射與接收方式的Rp為0.72，環(huán)形發(fā)射與接收檢測(cè)方式的Rp為0.72。在相同光照路徑下，光在大果中的光程要比小果的長(zhǎng)，衰減嚴(yán)重，接收到的光譜能量強(qiáng)度低，因此攜帶的有效信息比小果光譜的要少很多，所以小果的模型預(yù)測(cè)精度相對(duì)較高。

表2 不同尺寸樣品集的不同檢測(cè)方式PLS建模集預(yù)測(cè)結(jié)果Table 2 PLS modeling prediction results for different size sample sets by using different detection methods

糖度預(yù)測(cè)模型公式如式(4)所示。

(4)

其中，n為參與建模的光譜變量數(shù)，β為糖度相關(guān)能量譜強(qiáng)度，b為模型的截距。

2.5 臍橙糖度最佳模型建立和預(yù)測(cè)

通過(guò)以上多種模型的對(duì)比分析，漫透射檢測(cè)方式模型預(yù)測(cè)效果最好，多元散射校正預(yù)處理后，模型預(yù)測(cè)效果有所優(yōu)化，可以消除尺寸差異帶來(lái)的影響。采用多元散射校正預(yù)處理漫透射光譜，再利用偏最小二乘(PLS)建立漫透射及漫反射的預(yù)測(cè)模型如圖6所示。由圖6可知，經(jīng)過(guò)MSC校正的漫透射光譜預(yù)測(cè)模型散點(diǎn)距離擬合線較為集中，預(yù)測(cè)精度較高，建模相關(guān)系數(shù)為0.88，預(yù)測(cè)相關(guān)系數(shù)為0.84，預(yù)測(cè)集均方根誤差RMSEP為0.64，截距b=14.05?；貧w系數(shù)與近紅外漫透射大小果光譜曲線對(duì)比如圖7所示，光譜變量在PLS定量模型中貢獻(xiàn)率越大，回歸系數(shù)越大。所用的光譜變量與回歸系數(shù)加權(quán)求和再加上截距，及PLS模型預(yù)測(cè)值，實(shí)現(xiàn)了對(duì)樣品的糖度值的預(yù)測(cè)。

在漫反射檢測(cè)方式中，先進(jìn)行尺寸分選，區(qū)分出大小果后，再進(jìn)行大小果的糖度分選方式，也可以消除尺寸差異帶來(lái)的影響。

圖6 漫透射MSC-PLS建模散點(diǎn)圖Fig.6 Diffuse transmission MSC-PLS modeling scatter plot

3 結(jié) 論

研究了樣品尺寸差異對(duì)臍橙漫透射及漫反射在線無(wú)損檢測(cè)精度的影響。分析了臍橙的尺寸差異在漫透射和漫反射裝置上的不同影響以及臍橙三種不同檢測(cè)方式的近紅外光譜特性。分別建立了三種檢測(cè)方式下的大果、小果、混合果、經(jīng)MSC預(yù)處理后的混合果糖度PLS模型和尺寸PLS預(yù)測(cè)模型，其中小果的糖度PLS模型預(yù)測(cè)效果均最佳，在漫透射糖度預(yù)測(cè)模型中，經(jīng)MSC預(yù)處理后，混合果預(yù)測(cè)模型精度有所提高，預(yù)測(cè)效果要優(yōu)于小果，模型的預(yù)測(cè)相關(guān)系數(shù)為0.84，預(yù)測(cè)均方根誤差為0.64°Brix； 在不同檢測(cè)方式下的尺寸預(yù)測(cè)模型中，漫透射檢測(cè)方式要比兩種漫反射檢斥方式的尺寸模型預(yù)測(cè)精度低，在漫反射檢測(cè)方式下，先采用尺寸模型進(jìn)行尺寸進(jìn)行分選，然后再使用糖度模型進(jìn)行糖度分選，最終實(shí)現(xiàn)不同檢測(cè)方式下不同尺寸臍橙的同時(shí)在線分選。論文研究為不同果型的水果糖度同時(shí)在線無(wú)損檢測(cè)和分選提供了參考和理論支持。

圖7 回歸系數(shù)與漫透射大小果光譜對(duì)比Fig.7 Comparison of regression coefficientand diffuse transmission size