王 冬,韓 平*,吳靜珠,趙麗麗,徐 恒
1. 北京農(nóng)業(yè)質(zhì)量標(biāo)準(zhǔn)與檢測技術(shù)研究中心,北京 100097 2. 全國生態(tài)環(huán)保優(yōu)質(zhì)農(nóng)業(yè)投入品評價(jià)技術(shù)機(jī)構(gòu)(CAQS-TRP-004),北京 100097 3. 北京工商大學(xué)食品安全大數(shù)據(jù)技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室,北京 100048 4. 北京博普特科技有限公司,北京 100193
玉米種子不完善粒(Defective Kernals)指有缺陷或受到損傷但尚有使用價(jià)值的玉米顆粒。 中華人民共和國國家標(biāo)準(zhǔn)GB 1353—2018 《玉米》[1]規(guī)定, 1, 2, 3, 4和5等級玉米中, 不完善粒含量分別為不超過4.0%, 6.0%, 8.0%, 10.0%和15.0%; 中華人民共和國農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)NY/T 519—2002 《食用玉米》[2]規(guī)定, 1, 2和3等級玉米的不完善粒占比不超過5.0%。 熱損傷粒是指發(fā)熱或干燥受熱后籽粒顯著變色或受到損傷的顆粒, 是不完善粒的一種。 玉米干燥過程對溫度有嚴(yán)格的規(guī)定。 中華人民共和國國家標(biāo)準(zhǔn)GB/T 21017—2007《玉米干燥技術(shù)規(guī)范》規(guī)定[3], 食用玉米、 淀粉發(fā)酵工業(yè)用玉米、 飼料用玉米的允許受熱溫度分別為50, 55和60 ℃。 超過規(guī)定的干燥溫度很容易導(dǎo)致玉米種子熱損傷。 熱損傷玉米種子不僅在種子活力方面會受到影響, 而且在營養(yǎng)價(jià)值、 口感等方面亦會降低其價(jià)值。 中華人民共和國國家標(biāo)準(zhǔn)GB/T 5494—2019 《糧油檢驗(yàn) 糧食、 油料的雜質(zhì)、 不完善粒檢驗(yàn)》[4]規(guī)定了不完善粒的檢測方法, 包含電動篩選器法和手篩法。 然而, 目前對熱損傷粒的檢測仍需人工觀察玉米種子表皮、 胚、 胚乳是否有顯著變色, 不同檢驗(yàn)人員對玉米種子不完善粒定義理解存在偏差, 檢驗(yàn)時(shí)把握尺度也有差異, 玉米種子不完善粒檢驗(yàn)過程中很可能出現(xiàn)差異較大的情況[5]; 相對于蟲蝕粒、 病斑粒、 破損粒、 生芽粒、 生霉粒等其他種類不完善粒, 熱損傷粒外觀的改變相對而言不明顯, 故較其他不完善粒更難以準(zhǔn)確分辨。 因此, 針對玉米種子熱損傷粒, 建立一種客觀、 高效、 無損、 快速的檢測方法成為當(dāng)前亟待解決的問題。
隨著化學(xué)計(jì)量學(xué)的發(fā)展和計(jì)算機(jī)技術(shù)的進(jìn)步, 近紅外光譜技術(shù)近幾年飛速發(fā)展, 并在種子品質(zhì)無損快速檢測等領(lǐng)域得以廣泛應(yīng)用[6-7]。 與此同時(shí), 多光譜以及高光譜技術(shù)在提供光譜信息的同時(shí), 還可以提供樣品的空間分布信息[8-9], 從而為種子胚部、 胚乳部的局部光譜數(shù)據(jù)的采集提供了解決方案。 近年來, 一些學(xué)者就多光譜、 高光譜技術(shù)檢測種子品質(zhì)開展了有關(guān)研究, 并取得了一定的進(jìn)展。 Choi等[10]采用高光譜成像研究了羅勒種子溯源, 發(fā)現(xiàn)1 449~1 457, 1 242~1 254, 1 380和1 696 nm對產(chǎn)自新加坡、 巴基斯坦、 越南、 印度的羅勒種子產(chǎn)地溯源有重要作用。 張婷婷等[11]采用可見-近紅外高光譜成像結(jié)合連續(xù)投影算法和遺傳算法測定甜玉米種子電導(dǎo)率, 進(jìn)而對甜玉米種子活力進(jìn)行預(yù)測; 其中, 全譜數(shù)據(jù)經(jīng)多元散射校正預(yù)處理, 采用偏最小二乘回歸所建模型, 校正、 預(yù)測決定系數(shù)分別為0.983和0.974; 采用遺傳算法從全譜853個(gè)變量篩選出25個(gè)關(guān)鍵變量; 采用多元散射校正預(yù)處理, 基于遺傳算法篩選關(guān)鍵變量所建的偏最小二乘回歸模型的校正、 預(yù)測決定系數(shù)分別為0.976和0.973。 Nie等[12]采用近紅外高光譜成像結(jié)合深度卷積神經(jīng)網(wǎng)絡(luò)算法分別研究了雜交秋葵種子和雜交絲瓜種子的分類; 當(dāng)類別數(shù)量從2個(gè)增至6個(gè)時(shí), 雜交秋葵種子深度卷積神經(jīng)網(wǎng)絡(luò)模型的校正、 預(yù)測準(zhǔn)確率不低于98%, 雜交絲瓜種子深度卷積神經(jīng)網(wǎng)絡(luò)模型的校正、 預(yù)測準(zhǔn)確率不低于95%。 Jayme G. A. Barbedo等[13]采用近紅外高光譜成像研究了3個(gè)品種小麥籽粒發(fā)芽的檢測, 結(jié)果表明, 918和1 411 nm和小麥籽粒發(fā)芽損傷具有高度關(guān)聯(lián)性, 將上述兩個(gè)波長數(shù)據(jù)結(jié)合為一個(gè)指數(shù)即可準(zhǔn)確判斷小麥籽粒的發(fā)芽情況。 Collins Wakholi等[14]采用近紅外高光譜成像研究了玉米種子活力的快速檢測, 分別采用線性判別分析、 偏最小二乘判別分析和支持向量機(jī)建立玉米種子活力判別模型, 結(jié)果表明, 支持向量機(jī)模型準(zhǔn)確率最高, 白色、 紫色和黃色玉米種子活力判別準(zhǔn)確率分別為100%, 100%和98%。
糯玉米(waxy corn)又稱黏玉米或蠟質(zhì)型玉米, 其栽培技術(shù)簡單、 周期短, 可用作食材、 牲畜飼料以及工業(yè)原料等, 具有較高的經(jīng)濟(jì)價(jià)值、 營養(yǎng)價(jià)值和加工價(jià)值。 然而, 針對糯玉米種子熱損傷粒的無損快速鑒別目前尚未見報(bào)道。 中華人民共和國國家標(biāo)準(zhǔn)GB/T 22326—2008 《糯玉米》規(guī)定[15], 1, 2和3等級糯玉米中的不完善粒占比不超過6.0%; 中華人民共和國農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)NY/T 524—2002《糯玉米》規(guī)定[16], 1, 2和3等級糯玉米中的不完善粒占比不超過5.0%。 為糯玉米種子熱損傷粒的無損快速鑒別、 種子不完善粒無損快速檢測、 種子質(zhì)量快速篩查等提供具有參考價(jià)值的信息。以糯玉米為例, 采用多光譜成像技術(shù)研究糯玉米熱損傷粒的無損快速檢測, 分別以胚面向上、 胚面向下方式采集對照組和熱損傷組糯玉米種子多光譜成像數(shù)據(jù), 提取胚面向上胚部、 胚乳部單點(diǎn)多光譜數(shù)據(jù)以及胚面向下胚乳部單點(diǎn)多光譜數(shù)據(jù), 并融合胚面向上胚部和胚乳部多光譜數(shù)據(jù); 對各數(shù)據(jù)做基線校正預(yù)處理, 計(jì)算各光譜數(shù)據(jù)樣本標(biāo)準(zhǔn)差, 分析熱損傷對糯玉米種子不同部位多光譜數(shù)據(jù)影響的差異; 建立熱損傷糯玉米種子偏最小二乘判別分析(partial least square-discriminant analysis,PLS-DA)模型, 統(tǒng)計(jì)模型校正、 交互驗(yàn)證數(shù)據(jù)正確率。 作為對比, 采用近紅外光譜儀配合顆粒樣品采樣附件采集糯玉米種子胚面向上、 胚面向下近紅外光譜數(shù)據(jù), 并將糯玉米種子胚面向上及胚面向下近紅外光譜數(shù)據(jù)融合; 做基線校正預(yù)處理后, 計(jì)算各數(shù)據(jù)的各光譜數(shù)據(jù)樣本標(biāo)準(zhǔn)差, 分析熱損傷對糯玉米種子不同部位近紅外光譜數(shù)據(jù)影響的差異; 建立熱損傷糯玉米種子PLS-DA模型, 統(tǒng)計(jì)模型校正、 交互驗(yàn)證數(shù)據(jù)正確率。
京科糯2000玉米種子, 采購于京研益農(nóng)(北京)種業(yè)科技有限公司。 將糯玉米種子分為對照組和熱損傷組(實(shí)驗(yàn)組)各50粒, 對照組置于陰涼干燥處, 熱損傷組樣品放入電熱恒溫鼓風(fēng)干燥箱, 做高溫干燥處理, 模擬糯玉米種子熱損傷過程。
(1)多光譜成像儀: 丹麥Videometer公司, 型號: VideometerLab4。 采集多光譜圖像面積為90 mm×90 mm, 圖像有效像素2 192×2 192, 圖像空間分辨率41 μm·pixel-1; 采用19個(gè)高功率發(fā)光二極管作為光源, 按波長順序依次照射樣品, 光譜波長范圍365~970 nm, 光譜變量數(shù)19個(gè), 以儀器內(nèi)部自帶白色陶瓷片作為光譜參比, 積分球收集樣品反射光信息, 電荷耦合元件(charge-coupled device,CCD)檢測器對樣品反射光信號進(jìn)行檢測。
(2)近紅外光譜儀: 美國VIAVI公司, 型號: NIR1700。 采用鹵鎢燈作為光源, 采用線性漸變分光(linear variable filter,LVF)原理, 波長范圍908~1 676 nm, 光譜中心分辨率10 nm, 光譜變量數(shù)125個(gè), 以聚四氟乙烯白塊作為光譜參比, 光譜累加平均50次, 配合顆粒樣品采樣附件, 采用線陣列銦鎵砷檢測器(InGaAs detector)測量樣品漫反射近紅外光譜數(shù)據(jù)。
(3)電熱恒溫鼓風(fēng)干燥箱: 上?;芫G科學(xué)儀器有限公司, 型號: DGG-9246A, 溫度分辨率0.1 ℃, 恒溫波動度±1.0 ℃。
采用電熱恒溫鼓風(fēng)干燥箱對糯玉米種子熱損傷組樣品進(jìn)行高溫干燥處理, 溫度設(shè)定80.0 ℃, 干燥8 h。 將糯玉米種子胚面向上放置記為GU, 糯玉米種子胚面向下放置記為GD, 分別采用多光譜成像儀、 近紅外光譜儀采集糯玉米種子胚面和胚乳面多光譜成像數(shù)據(jù)和近紅外光譜數(shù)據(jù)。
(1)多光譜成像儀數(shù)據(jù)處理: 采用Videometer多光譜成像儀采集對照組和熱損傷組糯玉米種子多光譜成像數(shù)據(jù)。 其中, 糯玉米種子胚面向上放置(GU)時(shí)分別采集胚部和胚乳部單點(diǎn)多光譜數(shù)據(jù), 胚面向上胚部光譜數(shù)據(jù)記為GUEm, 胚面向上胚乳部光譜數(shù)據(jù)記為GUEn; 糯玉米種子胚面向下放置(GD)時(shí)采集胚乳部單點(diǎn)多光譜數(shù)據(jù), 記為GDEn。 將GUEm和GUEn數(shù)據(jù)初級融合為GUEm-GUEn數(shù)據(jù)。 對各組數(shù)據(jù)做基線校正預(yù)處理, 在此基礎(chǔ)上計(jì)算各光譜數(shù)據(jù)樣本標(biāo)準(zhǔn)差; 分別采用GUEm, GUEn, GDEn和GUEm-GUEn數(shù)據(jù)建立偏最小二乘-判別分析(partial least square-discriminant analysis,PLS-DA)模型, 并討論前3主成分得分分布情況以及各模型判別正確率。 糯玉米種子多光譜采集區(qū)域示意圖如圖1所示。
圖1 糯玉米種子多光譜采集區(qū)域示意圖(a): 胚面向上胚部GUEm(+)、 胚面向上胚乳部GUEn(×);(b): 胚面向下胚乳部GDEn(○)Fig.1 The diagram of the areas for collectingthe spectra of waxy corn seed
(2)近紅外光譜數(shù)據(jù)處理: 采用VIAVI近紅外光譜儀采集對照組和熱損傷組單粒糯玉米種子多光譜成像數(shù)據(jù)。 其中, 糯玉米種子胚面向上放置所采集光譜數(shù)據(jù)記為GUS, 胚面向下放置所采集光譜數(shù)據(jù)記為GDS。 對各組數(shù)據(jù)做基線校正預(yù)處理, 在此基礎(chǔ)上計(jì)算各光譜數(shù)據(jù)的樣本標(biāo)準(zhǔn)差; 分別采用GUS, GDS以及GUS-GDS初級融合數(shù)據(jù)建立PLS-DA模型, 并討論前3主成分得分分布情況以及各模型判別正確率。
圖2是京科糯2000玉米種子對照組(a)和熱損傷組(b)部分樣品。 采用Videometer多光譜成像儀采集樣品的多光譜成像后, 使用儀器自帶圖像處理軟件扣除背景得到圖2。 從圖2可見, 僅用肉眼觀察可見光照片難以將對照樣品和熱損傷樣品加以區(qū)分。
圖2 京科糯2000玉米種子對照組(a)和熱損傷組(b)Fig.2 The control group (a) and the heat-damagedgroup (b) of Jingkenuo 2000 corn seeds
京科糯2000玉米種子對照組和熱損傷組樣品的多光譜成像數(shù)據(jù)如圖3所示, 其中, (a1)為對照組胚面向上胚部光譜(GUEmDZ)、 (a2)為熱損傷組胚面向上胚部光譜(GUEmRSS)、 (b1)為對照 組胚面向上胚乳部光譜(GUEnDZ)、 (b2)為熱損傷組胚面向上胚乳部光譜(GUEnRSS)、 (c1)為對照組胚面向下胚乳部光譜(GDEnDZ)、 (c2)為熱損傷組胚面向下胚乳部光譜(GDEnRSS)。 從圖3可見, 各光譜從365~970 nm的吸光度基本呈上升趨勢, 并普遍存在970 nm處吸光度較高而365 nm處吸光度較低的趨勢。 另一方面, 熱損傷樣品和對照樣品的多光譜數(shù)據(jù)在365~970 nm范圍內(nèi)的光譜輪廓未見明顯區(qū)別。
圖3 糯玉米種子多光譜(a1): GUEmDZ; (a2): GUEmRSS; (b1): GUEnDZ; (b2): GUEnRSS; (c1): GDEnDZ; (c2): GDEnRSSFig.3 The multi-spectra of waxy corn seeds(a1): GUEmDZ; (a2): GUEmRSS; (b1): GUEnDZ; (b2): GUEnRSS; (c1): GDEnDZ; (c2): GDEnRSS
為研究多光譜數(shù)據(jù)中各變量所對應(yīng)吸光度值的變化規(guī)律, 對GUEmDZ, GUEmRSS, GUEnDZ, GUEnRSS, GDEnDZ和GDEnRSS數(shù)據(jù)分別做基線校正預(yù)處理以消除光譜平移對數(shù)據(jù)的影響, 再統(tǒng)計(jì)各光譜數(shù)據(jù)樣本標(biāo)準(zhǔn)差, 如圖4所示。 由于對數(shù)據(jù)做基線校正預(yù)處理, 波長365 nm所對應(yīng)變量數(shù)據(jù)的樣本標(biāo)準(zhǔn)差值為0。 在405~780 nm范圍內(nèi), GUEmDZ樣本標(biāo)準(zhǔn)差較GUEmRSS樣本標(biāo)準(zhǔn)差大, 在850~970 nm范圍內(nèi), GUEmDZ和GUEmRSS樣本標(biāo)準(zhǔn)差規(guī)律不明顯; 在405~970 nm范圍內(nèi), GUEnDZ樣本標(biāo)準(zhǔn)差普遍較GUEnRSS樣本標(biāo)準(zhǔn)差?。?在405~970 nm范圍內(nèi), 除波長690 nm所對應(yīng)的標(biāo)準(zhǔn)差值較為接近外, 其他波長對應(yīng)的GDEnDZ樣本標(biāo)準(zhǔn)差較GDEnRSS的樣本標(biāo)準(zhǔn)差小。 以上數(shù)據(jù)結(jié)果說明, 高溫削弱了糯玉米種子胚的差異, 在光譜上表現(xiàn)為熱損傷組糯玉米種子胚的多光譜數(shù)據(jù)變異相對較小而對照組糯玉米種子胚的多光譜數(shù)據(jù)變異相對較?。?另一方面, 高溫增強(qiáng)了糯玉米種子胚乳的差異, 在光譜上表現(xiàn)為熱損傷組糯玉米種子胚乳的多光譜數(shù)據(jù)變異相對較大而對照組糯玉米種子胚乳的多光譜數(shù)據(jù)變異相對較小。
圖4 糯玉米種子多光譜各變量數(shù)據(jù)樣本標(biāo)準(zhǔn)差○: GUEmDZ; *: GUEmRSS; △: GUEnDZ;▽: GUEnRSS; +: GDEnDZ; ×: GDEnRSSFig.4 The standard deviation of the multi-spectraof waxy corn seeds○: GUEmDZ; *: GUEmRSS; △: GUEnDZ;▽: GUEnRSS; +: GDEnDZ; ×: GDEnRSS
對造成上述結(jié)果的原因分析如下: 種子胚的主要成分是蛋白質(zhì), 蛋白質(zhì)的種類很多, 多以蛋白質(zhì)高級結(jié)構(gòu)(折疊、 螺旋、 疊加及組合等)的形式存在, 這是蛋白質(zhì)具有生物學(xué)功能的重要原因之一。 高溫可導(dǎo)致蛋白質(zhì)發(fā)生變性, 失去了原有的一些高級結(jié)構(gòu), 變性后的蛋白質(zhì)較變性前的蛋白質(zhì)在化學(xué)鍵組成、 分子結(jié)構(gòu)及分子構(gòu)象等方面的差異相對較小, 因此胚對照組光譜變量數(shù)據(jù)樣本標(biāo)準(zhǔn)差普遍高于熱損傷組; 而胚乳的主要成分是以淀粉為主的多糖類物質(zhì), 本研究的高溫環(huán)境導(dǎo)致熱損傷組糯玉米種子胚乳中部分多糖類物質(zhì)化學(xué)鍵斷裂, 形成由多種糖類構(gòu)成的復(fù)雜混合物, 其化學(xué)鍵的組成、 分子結(jié)構(gòu)及分子構(gòu)象較對照組樣品增大, 因此對于胚乳部多光譜數(shù)據(jù), 對照組光譜變量數(shù)據(jù)樣本標(biāo)準(zhǔn)差普遍低于熱損傷組。
對胚面向上胚部光譜(GUEm)PLS-DA模型、 胚面向上胚乳部光譜(GUEn)PLS-DA模型、 胚面向下胚乳部光譜(GDEn)PLS-DA模型、 胚面向上胚部光譜和胚乳部光譜融合數(shù)據(jù)(GUEm-GUEn)PLS-DA模型計(jì)算主成分累積貢獻(xiàn)率, 上述4個(gè)模型的前3主成分累積貢獻(xiàn)率分別為97%, 98%, 95%和93%, 亦即各模型潛3主成分累積貢獻(xiàn)率皆超過90%。 將GUEm, GUEn, GDEn, GUEm-GUEn四模型的前三主成分得分繪制3D散點(diǎn)圖, 分別如圖5(a—d)所示, 其中, “○”表示對照組樣品, “△”表示熱損傷組樣品。 從圖5可見, 四種模型的前三主成分得分的3D散點(diǎn)圖中, 對照組和熱損傷組具有一定的分離趨勢。
圖5 糯玉米種子多光譜PLS-DA判別模型前3主成分得分散點(diǎn)圖(a): GUEm; (b): GUEN; (c): GDEn; (d): GUEm-GUEn; ○: 對照組; △: 熱損傷組Fig.5 The scatter plots of the first 3 principal components of PLS-DA models of the multi-spectra of waxy corn seeds(a): GUEm; (b): GUEN; (c): GDEn; (d): GUEm-GUEn; ○: control group; △: heat-damaged group
分別采用胚面向上胚部光譜數(shù)據(jù)(GUEm)、 胚面向上胚乳部光譜數(shù)據(jù)(GUEn)、 胚面向下胚乳光譜數(shù)據(jù)(GDEn)建立京科糯2000玉米種子熱損傷粒的PLS-DA判別模型, 并建立GUEm和GUEn融合數(shù)據(jù)的熱損傷粒PLS-DA判別模型, 結(jié)果如表1所示。 從表1可見, GUEm, GUEn, GDEn和GUEm-GUEn四模型的校正數(shù)據(jù)正確判別率皆在96%~100%之間, 而交互驗(yàn)證數(shù)據(jù)正確判別率皆在92%~100%之間; 四模型中, GUEm-GUEn模型正確率最高。
表1 京科糯2000玉米種子熱損傷粒PLS-DA正確率Table 1 The accuracy of PLS-DA for heat-damagedkernels of Jingkonuo 2000 corn seeds
采用VIAVI近紅外光譜儀采集京科糯2000玉米種子對照組和熱損傷組樣品的近紅外光譜如圖6所示, 其中, a為對照組胚面向上光譜(GUSDZ)、 b為熱損傷組胚面向上光譜(GUSRSS)、 c為對照組胚面向下光譜(GDSDZ)、 d為熱損傷組胚面向下光譜(GDSRSS)。 從圖6可見, 各光譜呈相似的趨勢, 熱損傷樣品和對照樣品在908~1 676 nm范圍內(nèi)未見明顯區(qū)別。
圖6 糯玉米種子近紅外光譜(a): GUSDZ; (b): GUSRSS; (c): GDSDZ; (d): GDSRSSFig.6 The NIR spectra of waxy corn seeds(a): GUSDZ; (b): GUSRSS; (c): GDSDZ; (d): GDSRSS
為研究近紅外光譜數(shù)據(jù)中各變量所對應(yīng)吸光度值的變化規(guī)律, 對GUSDZ, GUSRSS, GDSDZ和GDSRSS數(shù)據(jù)分別做基線校正預(yù)處理以消除光譜平移對數(shù)據(jù)的影響, 再統(tǒng)計(jì)各光譜數(shù)據(jù)樣本標(biāo)準(zhǔn)差, 如圖7所示。 由于對數(shù)據(jù)做基線校正預(yù)處理, 波長964 nm所對應(yīng)樣本標(biāo)準(zhǔn)差值為0。 從圖7可見: 對于胚面向上數(shù)據(jù), GUSDZ樣本標(biāo)準(zhǔn)差普遍高于GUSRSS; 而對于胚面向下數(shù)據(jù), GDSDZ樣本標(biāo)準(zhǔn)差普遍低于GDSRSS; 該規(guī)律在1 400~1 676 nm范圍內(nèi)更為明顯; 該結(jié)果和多光譜樣本標(biāo)準(zhǔn)差的分析結(jié)果一致。
圖7 糯玉米種子近紅外光譜數(shù)據(jù)的樣本標(biāo)準(zhǔn)差○: GUSDZ; *: GUSRSS; △: GDSDZ; ▽: GDSRSSFig.7 The standard deviation of the near-infraredspectral data of waxy corn seeds○: GUSDZ; *: GUSRSS; △: GDSDZ; ▽: GDSRSS
分別對胚面向上(GUS)、 胚面向下(GDS)以及GUS-GDS融合數(shù)據(jù)建立PLS-DA模型, 并計(jì)算模型主成分累積貢獻(xiàn)率。 結(jié)果表明, GUS, GDS和GUS-GDS模型的前3主成分累積貢獻(xiàn)率分別為99%, 99%和96%, 亦即各模型前3主成分累積貢獻(xiàn)率皆超過95%。 將GUS, GDS和GUS-GDS模型的前三主成分得分繪制3D散點(diǎn)圖, 分別如圖8(a—c)所示, 其中, “○”表示對照組樣品, “△”表示熱損傷組樣品。 從圖8可見, 三個(gè)模型的前三主成分得分的3D散點(diǎn)圖中, 對照組和熱損傷組具有一定的分離趨勢, 其中GUS和GDS模型對照組和熱損傷組的分離趨勢更明顯。
圖8 糯玉米種子近紅外光譜PLS-DA判別模型前3主成分得分散點(diǎn)圖(a): GUS; (b): GDS; (c): GUS-GDS; ○: 對照組; △: 熱損傷組Fig.8 The scatter plots of the first 3 principal components of PLS-DA models of the NIR spectra of waxy corn seeds(a): GUS; (b): GDS; (c): GUS-GDS; ○: control group; △: heat-damaged group
分別采用胚面向上光譜數(shù)據(jù)(GUS)、 胚面向下光譜數(shù)據(jù)(GDS)以及GUS-GDS融合數(shù)據(jù)建立京科糯2000玉米種子熱損傷粒的PLS-DA判別模型, 結(jié)果如表2所示。 從表2可見, 各模型正確率皆為100%。
表2 京科糯2000玉米種子熱損傷粒PLS-DA正確率Table 2 The accuracy of PLS-DA for heat-damagedkernels of Jingkonuo 2000 corn seeds
分別采用多光譜成像儀、 近紅外光譜儀針對京科糯2000玉米種子熱損傷組、 對照組分別采集了多光譜成像數(shù)據(jù)和近紅外光譜數(shù)據(jù)。 針對多光譜成像數(shù)據(jù), 提取了胚面向上胚部、 胚面向上胚乳部單點(diǎn)多光譜數(shù)據(jù)以及胚面向下胚乳部單點(diǎn)多光譜數(shù)據(jù), 并將胚面向上胚部和胚乳部數(shù)據(jù)融合; 針對近紅外光譜數(shù)據(jù), 采集了單粒種子胚面向上、 胚面向下近紅外光譜數(shù)據(jù)并將單粒種子數(shù)據(jù)融合。 對各數(shù)據(jù)做基線校正預(yù)處理后, 計(jì)算樣本標(biāo)準(zhǔn)差, 并分析熱損傷和對糯玉米種子各部位的影響, 結(jié)果表明, 高溫過程對糯玉米種子胚、 胚乳有不同程度的影響, 表現(xiàn)在光譜數(shù)據(jù)樣本標(biāo)準(zhǔn)差上, 熱損傷組胚部多光譜數(shù)據(jù)樣本標(biāo)準(zhǔn)差較對照組普遍較小, 熱損傷組胚向上近紅外光譜變量數(shù)據(jù)樣本標(biāo)準(zhǔn)差較對照組普遍較小, 熱損傷組胚乳部多光譜數(shù)據(jù)樣本標(biāo)準(zhǔn)差較對照組普遍較大, 熱損傷組胚乳部近紅外光譜數(shù)據(jù)樣本標(biāo)準(zhǔn)差較對照組普遍較大, 亦即本研究中采用的近紫外-可見-短波近紅外多光譜數(shù)據(jù)和近紅外光譜數(shù)據(jù)具有相似的規(guī)律。 分別基于多光譜數(shù)據(jù)、 近紅外光譜數(shù)據(jù)建立糯玉米熱損傷粒PLS-DA模型, 結(jié)果表明, 多光譜數(shù)據(jù)PLS-DA的校正數(shù)據(jù)正確率在96%~100%之間, 交互驗(yàn)證數(shù)據(jù)正確率在92%~100%之間, 其中, 基于糯玉米種子胚向上胚部和胚乳部融合的多光譜數(shù)據(jù)所建PLS-DA模型準(zhǔn)確率最高, 校正數(shù)據(jù)正確率為100%, 交互驗(yàn)證數(shù)據(jù)正確率在98%~100%之間; 近紅外光譜數(shù)據(jù)PLS-DA模型的校正數(shù)據(jù)、 交互驗(yàn)證數(shù)據(jù)正確率皆為100%。
本研究表明, 基于近紫外-可見-短波近紅外多光譜數(shù)據(jù)可以對糯玉米種子熱損傷粒進(jìn)行無損快速鑒別, 而基于近紅外光譜數(shù)據(jù)所建糯玉米種子熱損傷粒判別模型具有更高的準(zhǔn)確度。 然而, 跨年度樣品、 不同高溫環(huán)境導(dǎo)致的熱損傷樣品以及其他品種玉米種子的熱損傷粒無損快速鑒別仍需進(jìn)一步深入研究。