基于高光譜成像技術(shù)的進(jìn)口冰鮮牛肉新鮮度指標(biāo)檢測(cè)

孫宗保，梁黎明，閆曉靜，鄒小波，王天真，劉小裕，李君奎

（江蘇大學(xué)食品與生物工程學(xué)院，江蘇 鎮(zhèn)江 212013）

在我國(guó)，牛肉需求遠(yuǎn)超過國(guó)內(nèi)供給，并且供需差距不斷擴(kuò)大[1]，進(jìn)口牛肉成為解決供需矛盾的重要途徑。冰鮮牛肉因其新鮮、柔嫩等特點(diǎn)深受消費(fèi)者喜愛[2-3]，但牛肉在非冷凍條件下腐敗速度較快，保鮮時(shí)間相對(duì)較短，在一定程度上制約了其在市場(chǎng)上的發(fā)展。進(jìn)口牛肉市場(chǎng)上存在不法商販為了謀取利益以次鮮肉冒充新鮮肉售賣，這不僅損害消費(fèi)者利益，還可能導(dǎo)致食品安全問題的發(fā)生。因此有必要對(duì)不同貯藏時(shí)間冰鮮牛肉的品質(zhì)檢測(cè)進(jìn)行研究。

目前對(duì)肉品品質(zhì)的研究仍主要采用感官評(píng)定、理化指標(biāo)檢測(cè)和微生物檢測(cè)等常規(guī)方法[4]。感官評(píng)價(jià)有主觀性強(qiáng)等缺陷，理化檢測(cè)與微生物檢測(cè)多數(shù)需要復(fù)雜的操作步驟，具有樣品處理繁瑣、檢測(cè)周期長(zhǎng)等問題[5]。為了滿足肉類品質(zhì)快速檢測(cè)與準(zhǔn)確分析的要求，無損檢測(cè)技術(shù)逐步代替常規(guī)檢測(cè)方法運(yùn)用到肉品品質(zhì)安全檢測(cè)中。其中高光譜成像技術(shù)將光譜與二維成像技術(shù)結(jié)合在一起，不僅可以獲得樣本在每個(gè)波段下的圖像信息并用于樣品的外部特征研究，還可以提取圖像上每個(gè)像素點(diǎn)的光譜信息并用于樣品的內(nèi)部品質(zhì)研究[6-8]。高光譜成像技術(shù)在肉品品質(zhì)檢測(cè)方面已有較多應(yīng)用[9]。Kamruzzaman等[10]利用高光譜成像技術(shù)識(shí)別碎牛肉中是否摻假雞肉及其摻假濃度；Crichton等[11]利用高光譜成像技術(shù)檢測(cè)不同貯藏條件下的牛肉的pH值，成功分離出變質(zhì)牛肉；Foca等[12]研究證明高光譜成像技術(shù)可以預(yù)測(cè)干切火腿表明乳酸菌的數(shù)量?，F(xiàn)階段對(duì)牛肉單一或多種指標(biāo)檢測(cè)的研究很多，但在從眾多理化指標(biāo)中篩選出不同貯藏時(shí)間的冰鮮牛肉之間的新鮮度特征指標(biāo)的研究還很少。通過對(duì)特征指標(biāo)含量的預(yù)測(cè)可以進(jìn)一步研究不同貯藏時(shí)間冰鮮牛肉的品質(zhì)差異。

本研究以不同貯藏時(shí)間的冰鮮牛肉為研究對(duì)象，重點(diǎn)研究牛肉貯藏過程中揮發(fā)性鹽基氮（total volatile basic nitrogen，TVB-N）[13]、顏色（L*、a*、b*）、脫氧肌紅蛋白（deoxymyolglobin，DMb）、氧合肌紅蛋白（oxymyoglobin，OMb）和高鐵肌紅蛋白（metmyoglobin，MMb）含量的變化，并找出不同貯藏時(shí)間冰鮮牛肉的新鮮度特征指標(biāo)，然后利用高光譜成像技術(shù)預(yù)測(cè)新鮮度特征指標(biāo)含量，比較不同預(yù)測(cè)模型的預(yù)測(cè)效果。最后利用最佳預(yù)測(cè)模型計(jì)算新鮮度特征指標(biāo)含量，結(jié)合圖像處理方法完成特征指標(biāo)含量的分布可視化[14]。

1 材料與方法

1.1 材料與試劑

澳大利亞進(jìn)口冰鮮牛腱肉購自上海新臻元食品有限公司。牛肉樣品在冰箱冷鮮層（（4±0.5）℃）進(jìn)行不同時(shí)間（1、15、30、45、60 d和90 d）的冷藏，記為F-1、F-15、F-30、F-45、F-60、F-90。

氧化鎂、硼酸、甲基紅指示劑、溴甲酚綠指示劑、磷酸氫二鈉、磷酸二氫鈉、氯化鉀、無水乙醇（均為分析純） 國(guó)藥集團(tuán)化學(xué)試劑有限公司；鹽酸標(biāo)準(zhǔn)溶液上海安譜實(shí)驗(yàn)科技股份有限公司。

1.2 儀器與設(shè)備

高光譜成像系統(tǒng)主要組成硬件包括：ImSpector V10E CCD攝像頭 芬蘭Specim公司；SC30021A精密電控平移臺(tái)、SC300-1A電子控制箱 北京卓立漢光儀器有限公司；Fiber-Lite DC-950 Illuminator 150 W光纖鹵素?zé)裘绹?guó)Dolan-Jenner公司。

Color Quest XE分光測(cè)色儀 美國(guó)Hunter Lab公司；A25-Digital分散機(jī) 上海歐河機(jī)械設(shè)備有限公司；TGL-16MC離心機(jī) 長(zhǎng)沙湘銳離心機(jī)有限公司；V-5100B分光光度計(jì) 上海元析儀器有限公司。

1.3 方法

1.3.1 樣本的預(yù)處理

實(shí)驗(yàn)前在無菌冷藏室去掉樣本中可見的脂肪、結(jié)締組織等，切成3 cm×3 cm×2 cm表面平整的塊狀樣本，編號(hào)后進(jìn)行指標(biāo)測(cè)定。每個(gè)貯藏時(shí)間下含有30 塊樣本，測(cè)定顏色（L*、a*、b*）并采集高光譜圖像后，再將其利用絞肉機(jī)充分?jǐn)囁楹筮M(jìn)行TVB-N和肌紅蛋白（DMb、OMb和MMb）的測(cè)定。

1.3.2 理化指標(biāo)測(cè)定

1.3.2.1 TVB-N含量的測(cè)定

參考GB 5009.228—2016《食品中揮發(fā)性鹽基氮的測(cè)定》中的半微量定氮法進(jìn)行測(cè)定。

1.3.2.2 顏色的測(cè)定

L*、a*、b*代表物體顏色的色度值，其中L*代表明暗度（黑白），a*代表紅綠色，b*代表黃藍(lán)色。利用分光測(cè)色儀測(cè)定牛肉樣本的L*、a*、b*值。待測(cè)色儀開機(jī)、校正后，對(duì)每個(gè)肉樣在不同肌肉部分測(cè)定3 次取其平均值。

1.3.2.3 肌紅蛋白含量的測(cè)定

參照Krzywicke[15]的方法，首先將肉樣充分?jǐn)囁椋? g碎肉加入20 mL、0.04 mol/L的磷酸鈉緩沖液（pH 6.8）中，然后用分散機(jī)在室溫下10 000 r/min均質(zhì)15 s。將均質(zhì)液在4 ℃冰箱中冷藏1 h后使用離心機(jī)在5 000 r/min、4 ℃條件下離心20 min。最后將上清液過濾，使用分光光度計(jì)分別在525、545、565、572 nm波長(zhǎng)處測(cè)量濾液的吸光度。肌紅蛋白含量按照以下公式計(jì)算：

式中：P1、P2、P3分別為DMb、OMb、MMb的質(zhì)量分?jǐn)?shù)/%；R1、R2、R3分別為吸光度比值A(chǔ)572/A525、A565/A525、A545/A525。

1.3.3 高光譜成像采集

高光譜采集參數(shù)設(shè)置為：曝光時(shí)間45 ms，圖像分辨率618×1 628，平移臺(tái)速率為90 mm/s，光譜分辨率2.8 nm，光譜波長(zhǎng)范圍為431～962 nm，采樣間隔0.858 nm。高光譜數(shù)據(jù)采集后將得到一個(gè)三維圖像數(shù)據(jù)塊。樣本采集過程中的噪聲、暗電流會(huì)對(duì)圖像造成影響，因此對(duì)原始圖像進(jìn)行校正后再提取數(shù)據(jù)。利用ENVI軟件打開牛肉樣本標(biāo)定后的高光譜圖像，采用矩形工具選取每個(gè)樣本中心附近的200×200感興趣區(qū)域（region of interest，ROI），計(jì)算ROI內(nèi)的所有像素點(diǎn)的光譜反射率作為此樣本的光譜數(shù)據(jù)。

1.4 數(shù)據(jù)處理與分析

1.4.1 統(tǒng)計(jì)分析方法

采用偏最小二乘法進(jìn)行定量模型的構(gòu)建。使用SPSS 25.0軟件對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行單因素方差分析，通過Duncan新復(fù)極差法進(jìn)行顯著性分析，P＜0.05，差異顯著；并使用SPSS進(jìn)行Pearson相關(guān)性分析。

1.4.2 光譜預(yù)處理方法

本實(shí)驗(yàn)選用的預(yù)處理方法包括一階導(dǎo)數(shù)（first derivative，1st Der）、二階導(dǎo)數(shù)（second derivative，2nd Der）、均值中心化（mean centering，MC）、多元散射校正（multiplicative scatter correction，MSC）、卷積平滑（savitzky-golay，SG）、標(biāo)準(zhǔn)正態(tài)變量變換（standard normal variate transformation，SNVT）和歸一化（normalization，N）。將未經(jīng)預(yù)處理組分設(shè)置為對(duì)照組。

1.4.3 特征變量提取方法

區(qū)間偏最小二乘（interval partial least squares，iPLS）方法：將整個(gè)光譜區(qū)域分成若干個(gè)等寬子區(qū)間，對(duì)每個(gè)子區(qū)間構(gòu)建偏最小二乘回歸模型，采用留一法交互驗(yàn)證計(jì)算交叉驗(yàn)證均方根誤差（root mean square error of cross validation，RMSECV）值，根據(jù)RMSECV最低的原則確定最佳建模區(qū)間[16]。

競(jìng)爭(zhēng)性自適應(yīng)重加權(quán)算法（competitive adaptive reweighted sampling，CARS）是一種蒙特卡羅（Monte Carlo）結(jié)合偏最小二乘（partial least squares，PLS）回歸系數(shù)的特征波長(zhǎng)選擇方法。它可以有效去除無關(guān)變量，減少共線性變量的影響[17]。

1.4.4 樣本集的劃分

采用Kennard-Stone劃分樣本集，以2∶1的比例劃分訓(xùn)練集和預(yù)測(cè)集，得到120 個(gè)樣本作為訓(xùn)練集，60 個(gè)樣本作為測(cè)試集。

2 結(jié)果與分析

2.1 理化指標(biāo)測(cè)定結(jié)果分析

為獲得較好的新鮮度預(yù)測(cè)性能，利用顯著性和相關(guān)性分析從諸多新鮮度指標(biāo)中篩選出與貯藏時(shí)間相關(guān)性最優(yōu)的特征指標(biāo)，然后利用預(yù)測(cè)模型計(jì)算新鮮度特征指標(biāo)含量，最后結(jié)合圖像處理方法完成特征指標(biāo)含量的分布可視化，以此獲得表征牛肉新鮮度的最優(yōu)可視化圖像。

圖1 不同貯藏時(shí)間冰鮮牛肉的TVB-N含量Fig. 1 TVB-N content of chilled beef during different storage times

從圖1可以看出，冰鮮牛肉的初始TVB-N含量為8.35 mg/100 g。隨著貯藏時(shí)間的延長(zhǎng)，冰鮮牛肉中的TVB-N含量顯著性增加（P＜0.05），與肉品腐敗程度呈正相關(guān)[18-20]。當(dāng)貯藏時(shí)間為45 d時(shí)，TVB-N含量接近15 mg/100 g，貯存至60 d時(shí)，TVB-N含量已經(jīng)超過國(guó)家規(guī)定。

圖2 不同貯藏時(shí)間冰鮮牛肉的顏色參數(shù)值Fig. 2 Color parameters of chilled beef during different storage times

由圖2可知，隨著貯藏時(shí)間的延長(zhǎng)，冰鮮牛肉的L*值變化波動(dòng)較大，貯藏初期（1、15 d）與后期（60、90 d）的L*值差異具有顯著性（P＜0.05）。b*值在貯藏15 d后達(dá)到最大，隨后逐漸降低，這與肉中微生物的繁殖以及水分損失有關(guān)。a*值整體上一直減小，但是在貯藏中期（15、30、45 d）變化不顯著（P＞0.05）。其原因是貯藏時(shí)間的延長(zhǎng)對(duì)氧合肌紅蛋白的形成不利，并且貯藏后期細(xì)菌的大量繁殖，促進(jìn)了高鐵肌紅蛋白的形成[21]，使牛肉a*值下降。

3 種肌紅蛋白的顏色分別為：DMb暗紅色（紫紅色）、OMb鮮紅色、MMb灰棕色（褐色）。根據(jù)謝小雷等[22]研究，OMb經(jīng)氧化后生成MMb，而MMb經(jīng)過高鐵肌紅蛋白還原酶系統(tǒng)的作用可以被還原為OMb，這兩個(gè)反應(yīng)同時(shí)進(jìn)行。由圖3可知，不同貯藏時(shí)間的OMb含量差異最顯著，經(jīng)短暫上升后逐漸下降。MMb含量變化與OMb正好相反，先下降后逐漸上升，導(dǎo)致肉色變暗、a*值下降。這與González-Calvo等[23]研究結(jié)果相符。

圖3 不同貯藏時(shí)間冰鮮牛肉的DMb、OMb和MMb含量Fig. 3 DMb, OMb and MMb contents of chilled beef samples during different storage times

2.2 新鮮度特征指標(biāo)篩選

對(duì)比圖1～3的顯著性分析結(jié)果，發(fā)現(xiàn)不同貯藏時(shí)間的TVB-N含量差異最顯著。如表1所示，TVB-N、DMb和MMb與貯藏時(shí)間呈正相關(guān)，L*、a*、b*和OMb與貯藏時(shí)間呈負(fù)相關(guān)。另外，通過對(duì)比相關(guān)系數(shù)發(fā)現(xiàn)，TVB-N、b*、a*的相關(guān)系數(shù)均大于0.8。因此本研究從顯著差異性、正負(fù)相關(guān)性、相關(guān)系數(shù)3 個(gè)方面綜合考慮，選取TVB-N、a*和b*作為不同貯藏時(shí)間冰鮮牛肉的新鮮度特征指標(biāo)。

表1 貯藏時(shí)間與理化指標(biāo)間的相關(guān)系數(shù)Table 1 Correlation coefficients between storage time and physicochemical indexes

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

圖4 冰鮮牛肉樣本的原始光譜圖（a）與平均光譜圖（b）Fig. 4 Raw Spectra (a) and average spectra (b) of chilled beef samples

從圖4可以看出，不同貯藏時(shí)間牛肉樣本之間的光譜曲線具有相似的趨勢(shì)，但是反射率存在差異。在560 nm波長(zhǎng)附近主要是脫氧肌紅蛋白的吸收峰，700 nm波長(zhǎng)附近的光譜吸收與NH3基團(tuán)的三級(jí)倍頻吸收有關(guān)，800 nm波長(zhǎng)附近的光譜吸收主要與水分子中O—H鍵的二級(jí)倍頻吸收有關(guān)[24-26]。光譜信息中存在著部分與樣本自身性質(zhì)不相關(guān)的冗雜信息，這會(huì)影響模型的穩(wěn)定性與預(yù)測(cè)的準(zhǔn)確度[27-28]。因此采用7 種不同的預(yù)處理方法（1st Der、2nd Der、MC、MSC、SG、SNVT和N）對(duì)原始光譜數(shù)據(jù)進(jìn)行預(yù)處理。

2.4 新鮮度特征指標(biāo)含量預(yù)測(cè)

根據(jù)2.2節(jié)的結(jié)論，將TVB-N、a*和b*作為不同貯藏時(shí)間的冰鮮牛肉的新鮮度特征指標(biāo)。利用PLS、iPLS、CARS-PLS模型對(duì)上述指標(biāo)含量進(jìn)行預(yù)測(cè)。在建立模型之前，采用Kennard-Stone方法進(jìn)行樣本集劃分，且訓(xùn)練集中樣本的每個(gè)特征指標(biāo)含量范圍包含測(cè)試集中樣本對(duì)應(yīng)的指標(biāo)含量，有助于預(yù)測(cè)模型精度的提升。

2.4.1 TVB-N含量的預(yù)測(cè)

2.4.1.1 PLS模型預(yù)測(cè)結(jié)果

如表2所示，當(dāng)原始光譜數(shù)據(jù)經(jīng)過SNVT方法預(yù)處理后，PLS模型的預(yù)測(cè)效果較好，訓(xùn)練集和測(cè)試集的相關(guān)系數(shù)rc和rp分別為0.938 2、0.925 4，RMSECV值和RMSEP值分別為1.46、1.49 mg/100 g。后續(xù)實(shí)驗(yàn)采用SNVT方法進(jìn)行預(yù)處理。

表2 不同預(yù)處理方法建立TVB-N含量PLS模型預(yù)測(cè)結(jié)果Table 2 Predictive results of PLS model for TVB-N content with different preprocessing methods

2.4.1.2 iPLS模型預(yù)測(cè)結(jié)果

首先將預(yù)處理后的618 個(gè)變量分成10～30 個(gè)區(qū)間，然后對(duì)每一區(qū)間建立PLS模型。當(dāng)i=10時(shí)，主成分?jǐn)?shù)為8，RMSECV值最小為0.851 7 mg/100 g，如圖5a所示。此時(shí)iPLS篩選的最優(yōu)子區(qū)間為第4個(gè)，圖5b顯示了最優(yōu)子區(qū)間的位置，對(duì)應(yīng)波長(zhǎng)范圍587.36～639.17 nm，一共62 個(gè)變量。將這些變量組合成新矩陣并建立PLS預(yù)測(cè)模型，預(yù)測(cè)結(jié)果如圖5c、d所示。與全光譜變量建立的預(yù)測(cè)模型結(jié)果相比較，訓(xùn)練集和測(cè)試集的識(shí)別率均有不同程度的下降。這是因?yàn)椋m然iPLS在提取變量過程中去除了大量的變量，但提取出來的變量是連續(xù)的，相鄰或者相互之間可能會(huì)存在高度相關(guān)性[29]。

圖5 TVB-N含量的iPLS模型的最佳主成分（a）、最優(yōu)子區(qū)間（b）、訓(xùn)練集（c）和測(cè)試集（d）散點(diǎn)圖Fig. 5 iPLS model for TVB-N content: optimum PCs (a), optimum subintervals (b), and scatter plots of training (c) and test set (d)

2.4.1.3 CARS-PLS模型預(yù)測(cè)結(jié)果

圖6 TVB-N含量的CARS-PLS模型的變量選擇Fig. 6 CARS-PLS model results of TVB-N content

利用CARS篩選波長(zhǎng)時(shí)，采樣次數(shù)設(shè)置為100，運(yùn)行結(jié)果如圖6和圖7所示。圖6a為不同采樣次數(shù)下選取的變量數(shù)量，其隨采樣次數(shù)的增加逐漸減??；由圖6b可知，RMSECV值隨著采樣次數(shù)的增加先減小后增大，當(dāng)采樣次數(shù)為35時(shí)，RMSECV值最??；圖6c表示波長(zhǎng)選擇過程中不同波長(zhǎng)變量對(duì)應(yīng)的回歸系數(shù)的變化路徑，*表示RMSECV值最小的位置，此時(shí)共篩選得到65 個(gè)波長(zhǎng)。將這65 個(gè)特征波長(zhǎng)變量組合后建立PLS預(yù)測(cè)模型，結(jié)果如圖7所示，rc和rp分別為0.965 8和0.963 7，RMSECV值和RMSEP值分別為1.06、1.12 mg/100 g。與全光譜建模結(jié)果相比，預(yù)測(cè)精度有所提高。說明通過CARS方法去除無關(guān)變量后，可以簡(jiǎn)化模型，降低模型的計(jì)算量，提高運(yùn)算效率并提高模型預(yù)測(cè)精度[30-31]。

圖7 TVB-N含量的CARS-PLS模型的最佳主成分?jǐn)?shù)（a）、訓(xùn)練集（b）和測(cè)試集（c）散點(diǎn)圖Fig. 7 CARS-PLS model results of TVB-N: optimum PCs (a), and scatter plots of training (b) and test set (c)

2.4.2a*和b*值的預(yù)測(cè)結(jié)果

2.4.2.1 PLS模型預(yù)測(cè)結(jié)果

表3 a*和b*值的PLS模型預(yù)測(cè)結(jié)果Table 3 Predictive results of PLS models for a* and b* values

如表3所示，對(duì)照組對(duì)a*的預(yù)測(cè)結(jié)果最好，說明利用原始光譜數(shù)據(jù)建立a*的PLS預(yù)測(cè)模型效果最佳，rc和rp分別為0.930 1、0.914 8，RMSECV值與RMSEP值分別為0.85、0.90；MSC是預(yù)測(cè)b*值的最佳預(yù)處理方法，rc和rp分別為0.912 8和0.902 0，RMSECV值與RMSEP值分別為1.04和1.24。

2.4.2.2 iPLS、CARS-PLS模型預(yù)測(cè)結(jié)果

表4 a*和b*值的iPLS、CARS-PLS and iCARS-PLS模型結(jié)果Table 4 Results of iPLS, CARS-PLS and iCARS-PLS models for a* and b * values

利用iPLS和CARS方法分別對(duì)經(jīng)過上述最佳預(yù)處理方法處理后的數(shù)據(jù)進(jìn)行特征提取，然后將變量組合成新數(shù)據(jù)矩陣并建立a*和b*值PLS預(yù)測(cè)模型，結(jié)果如表4所示。相比于iPLS和全光譜建立的PLS模型結(jié)果，CARS-PLS模型的預(yù)測(cè)精度有非常明顯的提升。此時(shí)對(duì)a*值和b*值進(jìn)行預(yù)測(cè)時(shí)，rc分別為0.949 5、0.964 2，rp分別為0.949 4、0.942 3，RMSECV值分別為0.71、0.76，RMSEP值分別為0.72、0.77。

2.5 新鮮度特征指標(biāo)含量分布可視化

利用CARS-PLS模型實(shí)現(xiàn)TVB-N、a*和b*值的可視化分布，便于更加直觀地分析牛肉的品質(zhì)變化。

如圖8所示，隨著貯藏時(shí)間的延長(zhǎng)，圖像顏色由藍(lán)色逐漸轉(zhuǎn)為黃色。圖8e、f基本為綠色和黃色，與圖8a～d差異明顯。對(duì)應(yīng)的牛肉樣本TVB-N含量均值大于15 mg/100 g，牛肉新鮮度明顯降低。

圖8 不同貯藏時(shí)間牛肉中TVB-N的可視化分布圖Fig. 8 TVB-N distribution maps of beef samples at different storage time pints

在圖9a、b中，黃色像素點(diǎn)和紅色像素點(diǎn)居多，a*值較大；圖9c紅色像素點(diǎn)基本消失，以黃色和淺藍(lán)色像素點(diǎn)為主；圖9d～f中像素點(diǎn)由淺藍(lán)色逐漸向深藍(lán)色轉(zhuǎn)變，表明a*值隨著貯藏時(shí)間的延長(zhǎng)逐漸下降。

圖9 不同貯藏時(shí)間牛肉中a*值的可視化分布圖Fig. 9 a* Distribution maps of beef samples at different storage time points

圖10 反映了b*值變化情況。隨著貯藏時(shí)間的延長(zhǎng)，圖像顏色整體由黃色向藍(lán)色轉(zhuǎn)變。其中，圖10a、b中黃色像素點(diǎn)增多，b*值增大，而圖10b～e中黃色像素點(diǎn)不斷減少，藍(lán)色像素點(diǎn)不斷增多，表明貯藏15 d后b*值開始減小。

圖10 不同貯藏時(shí)間牛肉中b*值的可視化分布圖Fig. 10 b* Distribution maps of beef samples at different storage time points

因此，通過新鮮度特征指標(biāo)的可視化分布圖可以更直觀地看出冰鮮牛肉在貯藏過程中新鮮度的變化趨勢(shì)，幫助消費(fèi)者以及生產(chǎn)者判斷牛肉的新鮮程度。

3 結(jié) 論

本研究以不同貯藏時(shí)間的冰鮮牛肉為研究對(duì)象，采集樣本的高光譜圖像、測(cè)定并分析樣本的TVB-N、顏色參數(shù)（L*、a*、b*）、DMb、OMb和MMb含量變化，結(jié)果表明TVB-N、a*、b*值與貯藏時(shí)間呈現(xiàn)較好的相關(guān)性，選取了TVB-N、a*和b*作為不同貯藏時(shí)間的冰鮮牛肉的新鮮度特征指標(biāo)。同時(shí)利用光譜數(shù)據(jù)對(duì)特征指標(biāo)含量進(jìn)行預(yù)測(cè)，最佳預(yù)測(cè)模型均為CARS-PLS。此時(shí)，模型的rp分別為0.963 7、0.949 4、0.942 3，RMSEP分別為1.12 mg/100 g、0.72、0.77。利用CARS-PLS模型結(jié)合圖像處理繪制了特征指標(biāo)含量可視化分布圖，可以直觀地看出冰鮮牛肉在貯藏過程中新鮮度的變化趨勢(shì)。結(jié)果表明，高光譜成像技術(shù)結(jié)合化學(xué)計(jì)量學(xué)方法可以實(shí)現(xiàn)冰鮮牛肉新鮮度指標(biāo)的快速檢測(cè)及其分布可視化。