引用本文
邵园园, 王永贤, 玄冠涛, 高宗梅, 刘艺, 韩翔, 胡志超. 高光谱成像的猕猴桃货架期快速预测[J]. 光谱学与光谱分析, 2020,40(6): 1940-1946.
SHAO Yuan-yuan, WANG Yong-xian, XUAN Guan-tao, GAO Zong-mei, LIU Yi, HAN Xiang, HU Zhi-chao. Hyperspectral Imaging Technique for Estimating the Shelf-Life of Kiwifruits[J]. Spectroscopy and Spectral Analysis, 2020,40(6): 1940-1946.
Doi:10.3964/j.issn.1000-0593(2020)06-1940-07  
Permissions
Copyright©2020, 《光谱学与光谱分析》期刊社
《光谱学与光谱分析》期刊社 所有
高光谱成像的猕猴桃货架期快速预测
邵园园1,2, 王永贤1, 玄冠涛1,3,*, 高宗梅4, 刘艺1, 韩翔1, 胡志超2,*
1. 山东农业大学机械与电子工程学院, 山东 泰安 271018
2. 农业农村部南京农业机械化研究所, 江苏 南京 210014
3. College of Agriculture, Food and Natural Resources, University of Missouri, Columbia 65211, USA
4. Biological Systems Engineering, Washington State University, Washington 99350, USA
*通讯联系人 e-mail: xuangt@sina.com; zchu369@163.com

作者简介: 邵园园, 1980年生, 山东农业大学机械与电子工程学院副教授 e-mail: syy007@sdau.edu.cn

收稿日期: 2019-05-07
基金: 国家自然科学基金项目(31671632, 31701325)资助
摘要

货架期是影响果蔬品质和供应安全的重要因素, 快速准确预测果蔬货架期已成为消费者、 生产者和管理者共同关注的问题。 猕猴桃含有多种有机物和氨基酸, 具有丰富的营养价值, 深受广大消费者的喜爱。 但由于猕猴桃表面颜色变化不明显, 人们仅凭感官难以准确判断猕猴桃的货架期和质量等级。 采用高光谱成像结合化学计量学方法对不同储存条件下的保鲜猕猴桃进行了货架期预测。 首先采集了4 ℃和(18±2) ℃下保鲜时间为0, 2, 4天各120个猕猴桃样本在400~1 000 nm的高光谱数据, 测定其硬度值和可溶性固形物含量(SSC), 获取猕猴桃切片高光谱图像。 对猕猴桃平均光谱提取并进行Savitzky-Golay卷积平滑预处理后, 通过光谱数据主成分分析(PCA), 发现不同货架期和储存温度的猕猴桃样本在前2个主成分空间形成一定的聚类, 4 ℃下猕猴桃样本出现少量重叠。 为了减少波长变量, 提高运算速度, 使用载荷系数法(XL)与连续投影算法(SPA)选择特征波长。 其中, 4 ℃猕猴桃样本的XL和SPA特征波长分别7个(481, 501, 547, 665, 723, 839, 912 nm)和10个(406, 428, 520, 617, 665, 682, 723, 818, 878和983 nm); (18±2) ℃猕猴桃样本XL特征波长为508, 545, 665, 672, 720, 839和909 nm, SPA特征波长为575, 622, 731, 756, 779, 800, 828, 865, 920和983 nm。 基于3:1的光谱数据集划分三个货架期虚拟等级值1, 2和3, 以全光谱数据、 特征波长为输入, 建立非线性最小二乘支持向量机(LS-SVM)预测模型。 结果表明, 对于4 ℃下3种货架期猕猴桃样本, 在全光谱和XL, SPA特征波长上的预测集准确度分别达到92.2%, 92.2%, 91.1%; (18±2) ℃时预测准确度均为100%。 猕猴桃切片图像PCA分析显示, 除PC5中有部分噪声影响外, 其他主成分图像均能完整反映猕猴桃切片信息, PC2图像可以明显呈现出猕猴桃切片在不同货架期的变化程度。 进一步分析猕猴桃硬度和可溶性固形物含量发现, 随着货架期延长, 猕猴桃可溶性固形物含量逐渐增加, (18±2) ℃时二者存在正相关性, 相关系数为0.557 6。 硬度则随货架期延长逐渐减小, 4和(18±2) ℃下硬度值和货架期之间存在负相关性, 相关系数分别为-0.335 6和-0.562 0。 结合猕猴桃光谱信息, 可以发现猕猴桃光谱反射率与其单个理化指标不成线性关系, 而是多个指标的综合反映。 因此, 采用高光谱成像技术可以全面、 准确、 快速的预测猕猴桃货架期, 为猕猴桃的生产、 销售提供技术指导。

关键词: 猕猴桃; 货架期; 近红外高光谱技术; 化学计量学
中图分类号:TP391.4 文献标志码:A
Hyperspectral Imaging Technique for Estimating the Shelf-Life of Kiwifruits
SHAO Yuan-yuan1,2, WANG Yong-xian1, XUAN Guan-tao1,3,*, GAO Zong-mei4, LIU Yi1, HAN Xiang1, HU Zhi-chao2,*
1. College of Mechanical and Electrical Engineering, Shandong Agricultural University, Tai’an 271018, China
2. Nanjing Research Institute for Agricultural Mechanization, Ministry of Agriculture and Rural Affairs, Nanjing 210014, China
3. College of Agriculture, Food and Natural Resources, University of Missouri, Columbia 65211, USA
4. Biological Systems Engineering, Washington State University, Washington 99350, USA
*Corresponding authors
Abstract

The shelf-life of fruits and vegetables is an important factor that affects the quality, which is concerned by the consumers, farmers and producers. Kiwifruit contains a variety of organic substances and amino acids, which has rich nutritional value and is deeply loved by consumers. However, due to its own characteristics such as the color characteristics of kiwifruit, it is difficult for consumers to make an accurate judgment on the edible degree of kiwifruit in the shelf-life by sensory evaluation. Therefore, non-destructive testing of the shelf life of fruits and vegetables is vital for agricultural products. In this research, hyperspectral imaging technology with chemometric methods was employed to estimate the shelf-life of kiwifruits which were stored in 4 ℃ and (18±2) ℃ among 3 periods (0, 2, 4 d). The spectral data covering the range of 400~1 000 nm were collected from 720 kiwifruit samples of 3 periods at 4 ℃ and (18±2) ℃. Meanwhile, the firmness and solid soluble content (SSC) of kiwifruits were measured, and the spectral data of kiwifruit slices were collected. The mean spectra (90 kiwifruits in the training set and 30 kiwifruits in prediction set) were extracted from each kiwifruit. Then, principal component analysis (PCA) was implemented for samples stored at different temperatures. Cluster analysis was performed based on PC1, while some overlap phenomenon showed in kiwi samples at 4 ℃. X-loadings of principal component analysis (PCA) and successive projection algorithm (SPA) method were applied to select the effective wavelengths, which are helpful for enhancing computer velocity. Based on X-loadings, 7 wavelengths (481, 501, 547, 665, 723, 839, 912 nm) were selected for samples stored at 4 ℃ and 7 wavelengths (508, 545, 665, 672, 720, 839, 909 nm) were selected for samples stored at (18±2) ℃, respectively. Similarly, for the SPA method, 10 wavelengths (406, 428, 520, 617, 665, 682, 723, 818, 878, 983nm) were selected for samples stored at 4 ℃ and 10 wavelengths (575, 622, 731, 756, 779, 800, 828, 865, 920, 983 nm) were selected for samples stored at (18±2) ℃, respectively. Thereafter, virtual levels (1, 2, 3) were assigned to the samples of 3 periods at 4 ℃ and (18±2) ℃, respectively. Least square-support vector machine (LS-SVM) was used to build classification models on full spectral data, effective wavelengths selected based on PCA and SPA, respectively. The results showed that the accuracy of the predictions reached to 92.2%, 92.2% and 91.1% among 3 periods at 4 ℃ and the accuracy of the predictions reached to 100% among 3 periods at (18±2) ℃, respectively. Also, the firmness and SSC of kiwifruits were measured and analyzed by one-way analysis of variance (ANOVA), the results showed that there was a negative correlation between firmness and shelf-life and the correlation coefficient was -0.335 6 and -0.562 0 at 4 ℃ and (18±2) ℃, respectively. There was a positive correlation between SSC and shelf-life and the correlation coefficient was 0.557 6 at (18±2) ℃. The shelf-life of kiwifruits can be estimated by the firmness index of samples stored atboth 4 ℃ and (18±2) ℃. While the SSC of samples stored at (18±2) ℃ was a significant estimation index. Further, the images of PC1—PC7 can preserve the integrity of the kiwifruit slice surface information, PC2 image can clearly show the degree of kiwifruit slices with different shelf-life. The results of this study indicate that it is feasible to use the hyperspectral imaging technique combined with the chemometric methods to classify the shelf-life of kiwifruits. Meanwhile, this research realized the rapid prediction of the shelf-life of kiwifruits and provided theoretical support for the quality and classification of fruit and vegetable shelf-life. Further, this study help forproviding technical supports for the developed instruments used for real time estimating the shelf-life of fruits and vegetables in further study.

Keyword: Kiwifruit; Shelf-life; Near-Infrared hyperspectral imaging technique; Chemometrics methods
引言

猕猴桃酸甜可口、 质地柔软, 富含人体所需的维生素、 氨基酸、 抗氧化剂和矿物质等多种营养成分, 深受消费者喜爱, 被誉为“水果之王”[1]。 猕猴桃属于呼吸跃变型果实, 采后成熟、 衰老迅速, 极易软化腐烂变质, 货架寿命非常有限[2]。 为此, 超市等零售场所通常使用保鲜膜或保鲜盒储存猕猴桃, 以降低其代谢水平, 延长货架期[3, 4]。 猕猴桃的果品品质与供应安全与其货架期密切相关, 但是, 由于猕猴桃表面颜色变化不明显, 人们仅凭感官难以准确判断猕猴桃的货架期和质量等级。 此外, 不同阶段货架期的猕猴桃放置在一起, 果品品质参差不齐, 严重影响其销售周期和利润。 因此, 快速准确预测保鲜猕猴桃所处货架期成为消费者、 生产者和管理者共同关注的一个技术问题。

传统的货架期人工检测方法操作繁琐、 实验条件要求高[5, 6]。 高光谱成像检测技术融合了样本的空间和光谱信息, 具有简便、 快速、 无损、 准确等优点, 已经成功应用于果蔬、 肉类、 水产等农副产品的货架期预测[7, 8]。 Chaudhry等[9]利用高光谱成像与多变量加速货架期实验法预测芝麻菜叶货架期。 Siripatrawa等[10]采用高光谱成像、 偏最小二乘回归分析方法(partial least squares regression, PLSR)及人工神经网络同时对袋装小香肠的品质和货架期进行了预测和评价; 周莉萍等[11]利用高光谱技术对覆盖保鲜膜的新鲜菠菜叶片货架期进行预测, 建立了偏最小二乘法判别分析(Partial least squares discrimination analysis, PLS-DA)模型。 Taghizadeh等[12]利用高光谱成像对不同包装膜的双孢蘑菇货架期保鲜进行评价, 建立了蘑菇品质指标和平均光谱的线性模型, Wang等[13]基于“大麦克”香蕉不同褐变程度进行货架期预测, 融合样本图像特征和平均光谱信息, 建立了香蕉货架期分类预测模型。 现有猕猴桃高光谱检测主要集中在品质参数、 损伤识别等方面[14, 15, 16, 17], 利用高光谱成像技术预测猕猴桃货架期鲜见报道。

本工作利用近红外高光谱成像技术获取冷藏、 室温条件下不同保鲜时间的猕猴桃高光谱信息, 结合化学计量学方法, 对猕猴桃货架期进行快速预测和判别。

1 实验部分
1.1 猕猴桃样本

2018年9月20日在山东省泰安市水果批发市场购买之后立即运回实验室。 选取形状大小均匀、 无腐烂、 无疤痕、 颜色正常的猕猴桃作为实验样本, 6个猕猴桃一组储存在超市普遍使用的聚对苯二甲酸类塑料(PET)保鲜盒中(图1), 共120盒(720个)猕猴桃样本。 保鲜盒尺寸为21.5 cm×16.5 cm×6 cm, 表层有10个直径为0.5 cm的圆形小孔, 以保证猕猴桃在储存过程中正常的呼吸作用和新陈代谢。 为分析比较冷藏和室温两种条件下猕猴桃货架期, 各选取60盒猕猴桃在4 ℃恒温储存室及(18±2) ℃室温下保鲜。 当猕猴桃储存0 d(购买当天), 储存2 d, 储存4 d时, 分别取出20盒样本共120个猕猴桃进行高光谱数据采集。 其中, 建模集90个, 预测集30个。

图1 猕猴桃样本Fig.1 Kiwifruit samples

1.2 高光谱图像采集

1.2.1 高光谱成像系统

采用GaiaField便携式高光谱系统(双利合谱, 四川, 中国)采集猕猴桃高光谱信息, 系统组成主要包括高光谱成像仪(GaiaField-V10E)、 成像镜头(HSIA-OL23)、 专用光源(HSIA-LS-T-200W)、 标准白板(HSIA-CT-150×150)、 三脚架(HSIA-TP-S)及装有高光谱数据采集软件(SpecView)的专用计算机等。 光谱范围为400~1 000 nm, 光谱分辨率2.8 nm, 入射狭缝宽30 μm, 视场角22° , CCD像素1 394×1 040, 光源对称分布, 入射角度45° 。

1.2.2 图像采集和校准

从保鲜盒取出猕猴桃, 待其恢复至常温状态且表面无水分残留后, 采集光谱数据。 为了得到清晰不失真图像, 系统参数设置为: 相机曝光时间13 ms, 镜头与样本间距离210 mm。 为了消除相机暗电流、 光照不均等对图像的影响, 需要对高光谱成像系统进行黑白校正。 通过遮盖镜头、 扫描标准白板分别获得全黑标定图像Idark和全白图像Iwhite, 根据式(1)获得校正图像[18]。

I0=Iraw-IdarkIwhite-Idark(1)

式中, I0为校准后图像, Iraw为原始高光谱图像。

利用软件ENVI4.6(Environment for Visualizing Images software, Research Systems Inc., Boulder, CO, USA)提取校正图像感兴趣区域(region of interst, ROI)的高光谱数据, 并计算ROI内光谱反射率的平均值。

1.3 理化指标测定

猕猴桃硬度测定使用指针式水果硬度计(GY-1型, 浙江托普云农科技股份有限公司), 测头直径3.5 mm, 下压深度10 mm, 测试速度0.5 mm·s-1, 选取5个位置进行测量, 取平均值作为猕猴桃硬度值。 硬度测定完成后, 使用ST-208B型手动切片机[19]沿猕猴桃中心线进行切片, 获取切片高光谱图像, 随后使用数显折射仪(PAL-1, Atago Co, Tokyo, Japan)进行可溶性固形物(soluble solids content, SSC)测定, 三次测量平均值为猕猴桃可溶性固形物含量测定值。

2 结果与讨论
2.1 猕猴桃光谱特征

为了进一步消除暗电流、 噪声等对高光谱信息的影响, 对猕猴桃平均光谱曲线进行Savitzky-Golay卷积平滑(法)预处理如图2。 图2(a)为不同货架期猕猴桃平均光谱曲线, 可以看出, 不同货架期猕猴桃平均反射率光谱曲线显著不同, 储存2 d的平均光谱反射率最高, 0 d的最低。 进一步分析, 光谱在650 nm附近出现小吸收峰, 是由叶绿素对光的吸收引起。 在780和970 nm附近出现的小吸收峰, 分别为水的O—H三级和二级吸收倍频引造成[20]。 图2(b)为不同温度下猕猴桃平均光谱曲线, 可以看出, (18±2) ℃室温下猕猴桃的光谱反射率要高于4 ℃恒温冷藏时的光谱反射率。

图2 猕猴桃平均光谱曲线
(a): 不同货架期; (b): 不同温度Fig.2 Mean spectral curvesof kiwifruits
(a): Different shelf-life; (b): Different temperature

2.2 PCA分析

主成分分析(principal component analysis, PCA)通过降维和特征提取实现数据集的简化, 每个主成分(principal component, PC)是所有原始变量的线性组合, 几个主成分(PCs)标准正交并在PCA变换中产生[21]。 图3(a)为不同货架期猕猴桃的前2个主成分得分图。 从图中可以看出, PC-1和PC-2解释了98%的光谱数据信息, 3组样本存在明显的聚类现象, 几乎没有重合。 图3(b)为样本PC-1得分图, 反映不同货架期猕猴桃光谱特征随储存温度的变化而变化。 随着货架期的延长, (18±2) ℃下猕猴桃的PC-1得分逐渐高于4 ℃下, 相应的, 室温下不同货架期猕猴桃具有更好的区分度。 为了实现不同温度下猕猴桃货架期的准确定性分析, 还需要建立猕猴桃货架期的预测模型。

图3 主成分得分散点图
(a): 不同货架期; (b): 不同温度Fig.3 Scores scatter plot of the first two PCs
(a): Different shelf-life; (b): Different temperature

2.3 特征波长XL提取

载荷系数法(XL)通过计算各波长所对应的载荷系数, 从而筛选影响样本聚类的特征波长, 前两个主成分PC-1和PC-2的载荷系数曲线如图4所示。

图4 主成分PC-1和PC-2载荷图
(a): 4 ℃; (b): (18±2) ℃Fig.4 Loadings plots of the PC-1 and PC-2 of PCA
(a): 4 ℃; (b): (18±2) ℃

图4(a)为储存温度为4 ℃时猕猴桃的PC-1和PC-2载荷系数曲线, 分析载荷系数局部最值, 可以发现481, 501, 547, 665, 723, 839和912 nm这些波长对不同货架期猕猴桃样聚类起重要作用。 图4(b)为储存温度(18±2) ℃下猕猴桃的PC-1和PC-2载荷系数曲线, 可以看出7个特征波长分别在508, 545, 665, 672, 720, 839和909 nm处。

2.4 特征波长SPA提取

SPA算法是一种特征波长前向选择算法, 通过比较波长投影向量大小, 将最大投影量波长列为有效波长, 并根据校正模型确定最佳特征波长。 以建模集光谱数据为输入, 指定波长范围N最大为30, 选取变量数为10时, 均方根误差(RMSE)最小。 采用SPA方法选择的特征波长如图5所示, 4 ℃下特征波长为406, 428, 520, 617, 665, 692, 723, 818, 878和983 nm, (18±2) ℃下特征波长为575, 622, 731, 756, 779, 800, 828, 865, 920和983 nm。

图5 SPA优选波长分布Fig.5 Distribution of wavelengths selected by SPA

XL, SPA两种方法选取的特征波长如表1所示。 其中, XL方法选取的特征波长数为7个, 两种温度下特征波长位置基本相同, 主要分布在500~900 nm范围内, SPA方法选取的特征波长数为10个, 分布范围在400~1 000 nm, 4 ℃下特征波长在400~600 nm范围的多一些, 与叶绿素有关, 而(18±2) ℃下特征波长更多是在800~100 nm, 由水或氧的窄吸收带引起。

表1 X-loadings和SPA算法选择的特征波长 Table 1 Characteristic wavelength selected by X-loadingsand SPA
2.5 LS-SVM预测模型

最小二乘支持向量机(least square-support vector machine, LS-SVM)作为一种常用的建模方法, 在数据降维及信息提取方面性能良好。 其中, 核函数选择径向基函数(RBF), 模型参数通过网格搜索和十折交叉验证进行优化。 为了建立猕猴桃货架期LS-SVM预测模型, 首先对货架期0, 2和4 d的猕猴桃样本分别赋予虚拟等级值1, 2, 3。 然后, 以全光谱数据和表1特征波长为输入, 建立LS-SVM预测模型, 猕猴桃货架期识别结果如表2所示。

表2 货架期识别结果 Table 2 Identification results for shelf life

由表2可以看出, 4 ℃下时基于全光谱和XL、 SPA特征波长的三个LS-SVM模型总体识别准确率分别为92.2%, 92.2%, 91.1%, 取得了较为满意的预测结果, 但第2天和第4天的猕猴桃样本存在小部分误判, 这也从侧面说明猕猴桃在适宜的低温下储存可以延长猕猴桃货架期。 (18±2) ℃下猕猴桃样本总体识别准确率均为100%, 无误判现象。 进一步分析, 与10个SPA特征波长相比, XL特征波长仅为7个, 可以使模型进一步简化, 提高运算效率, 便于开发便携式货架期检测仪。

2.6 猕猴桃切片图像分析

对猕猴桃切片图像进行PCA分析, 图6为4 ℃和(18±2) ℃的猕猴桃切片PC1—PC7图像。 由图中可以看出, 除PC5中有部分噪声影响外, 其他主成分图像均能完整反映猕猴桃切片信息。 其中, PC2图像反映猕猴桃切片信息最明显, 4 ℃和(18±2) ℃的猕猴桃切片PC2图像在各货架期均呈现出不同程度的内部变化。 从所得切片信息也进一步验证了高光谱成像技术是实现猕猴桃货架期预测的可靠工具。

图6 猕猴桃切片图像主成分分析Fig.6 Principal component images for kiwifruit slice

2.7 猕猴桃理化指标分析

表3和表4为不同货架期猕猴桃硬度、 可溶性固形物含量统计。 由表中可以看出, 随着货架期的延长, 猕猴桃硬度会下降, 这是由于淀粉等多糖物质逐渐降解导致。 相反, 随着储藏时间延长, 猕猴桃逐渐成熟, 可溶性固形物含量不断升高。 与(18±2) ℃室温相比, 4 ℃恒温冷藏条件下, 猕猴桃硬度、 可溶性固形物含量变化缓慢。 进一步采用单因素方差分析(analysis of variance, ANOVA)对猕猴桃硬度和可溶性固形物含量进行统计分。 表3中, 4 ℃的猕猴桃样本, F> F crit, p< 0.05, 说明不同货架期猕猴桃的硬度值存在显著性差异。 (18±2) ℃的猕猴桃样本, F> F crit, p< 0.01, 样本之间的硬度值存在高度显著性差异。 表4中, 4 ℃的猕猴桃样本, F< F crit, p> 0.05, 说明不同货架期猕猴的可溶性固形物含量值无显著性差异。 (18±2) ℃的猕猴桃样本, F> F crit, p< 0.01, 说明不同货架期样本的可溶性固形物含量值存在高度显著性差异。 相关性数据的邓肯式新复极差测验表明, 3个不同货架期下猕猴桃硬度值在0.05水平上存在显著性差异, 4 ℃下样本之间相关系数为-0.336 5, (18±2) ℃下样本之间相关系数为-0.526 0, 呈负相关。 3个不同货架期下猕猴桃可溶性固形物含量值在0.05水平上存在显著性差异, (18±2) ℃下样本之间相关系数0.557 6, 呈正相关。 由此可见, 硬度值均可作为区分4 ℃和(18±2) ℃不同货架期猕猴桃的特征指标, 而可溶性固形物含量仅能作为区分(18±2) ℃不同货架期猕猴桃的理化参数。 结合猕猴桃光谱信息, 可以看出猕猴桃光谱反射率与其单个指标硬度、 可溶性固形物含量不成线性关系, 而是多个指标的综合反映。

表3 不同货架期猕猴桃硬度 Table 3 Firmness of kiwifruits during shelf life
表4 不同货架期猕猴桃可溶性固形物含量 Table 4 SSC of kiwifruits during shelf life
3 结论

(1)采用高光谱成像技术结合化学计量方法, 对不同储存温度的猕猴桃进行货架期快速预测, 4 ℃时基于全光谱和XL、 SPA特征波长的三个LS-SVM模型各获得92.2%, 92.2%, 91.1%的预测准确率, (18±2) ℃时预测准确率均达到100%。 不同货架期猕猴桃切片主成分图像信息也显示出明显变化, 表明高光谱成像技术能够实现猕猴桃货架期快速预测。

(2)猕猴桃理化指标显示, 随着货架期延长, 猕猴桃可溶性固形物含量逐渐增加, (18±2) ℃时二者存在正相关性, 相关系数为0.557 6。 硬度则随货架期延长逐渐减小, 4 ℃和(18±2) ℃下硬度值和货架期之间存在负相关性, 相关系数分别为-0.335 6和-0.562 0。 与理化检测方法相比, 高光谱成像技术可以全面、 准确、 快速的预测猕猴桃货架期, 为猕猴桃的生产、 销售提供技术指导。

参考文献
[1] YU Jing-fen, SHANG Hai-tao, LING Jian-gang, et al(俞静芬, 尚海涛, 凌建刚, 等). Farm Products Processing(农产品加工), 2018, (2): 4. [本文引用:1]
[2] YAN Gen-zhu, WANG Chun-sheng, WANG Hua-rui, et al(阎根柱, 王春生, 王华瑞, 等). Chinese Agricultural Science Bulletin(中国农学通报), 2018, 34(22): 52. [本文引用:1]
[3] Tylewicz U, Fito P J, Castro-Giráldez M, et al. Journal of Food Engineering, 2011, 104(3): 438. [本文引用:1]
[4] Lara M A, Lleó L, Diezma-Iglesias B, et al. Journal of Food Engineering, 2013, 119(2): 353. [本文引用:1]
[5] LIU Ming, PAN Lei-qing, TU Kang, et al(刘明, 潘磊庆, 屠康, 等). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2010, 26(4): 317. [本文引用:1]
[6] WANG Rui-qing, ZHANG Ji-shu, MA Shu-shang, et al(王瑞庆, 张继澍, 马书尚, 等). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2009, 25(4): 243. [本文引用:1]
[7] CHEN Xiao-yu, ZHU Zhi-qiang, ZHANG Xiao-shuan, et al(陈晓宇, 朱志强, 张小栓, 等). Transactions of the Chinese Society of Agricultural Machinery(农业机械学报), 2015, 46(8): 192. [本文引用:1]
[8] Yang Q, Sun D W, Cheng W. Journal of Food Engineering, 2017, 192: 53. [本文引用:1]
[9] Chaudhry M M, Amodio M L, Babellahi F, et al. Journal of Food Engineering, 2018, 238: 122. [本文引用:1]
[10] Siripatrawan U, Makino Y. Meat Science, 2018, 146: 26. [本文引用:1]
[11] ZHOU Li-ping, ZHAO Yan-ru, YU Ke-qiang, et al(周莉萍, 赵艳茹, 余克强, 等). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2017, 37(2): 423. [本文引用:1]
[12] Taghizadeh M, Gowen A, Ward P, et al. Innovative Food Science and Emerging Technologies, 2010, 11(3): 423. [本文引用:1]
[13] Wang N N, Yang Y C, Sun D W, et al. Food Analytical Methods, 2015, 8(5): 1173. [本文引用:1]
[14] Guo W C, Zhao F, Dong J L. Food Analytical Methods, 2016, 9(1): 38. [本文引用:1]
[15] Liu D Y, Guo W C. Food Analytical Methods, 2015, 8(1): 164. [本文引用:1]
[16] GAO Zhen-hong, ZHANG Zhi-wei, HUO Ying-qiu, et al(高振红, 张志伟, 霍迎秋, 等). Transactions of the Chinese Society for Agricultural Machinery(农业机械学报), 2016, 47(6): 222. [本文引用:1]
[17] CHI Qian, WANG Zhuan-wei, YANG Ting-ting(迟茜, 王转卫, 杨婷婷, 等). Transactions of the Chinese Society for Agricultural Machinery(农业机械学报), 2015, 46(3): 235. [本文引用:1]
[18] Li X L, Wei Y Z, Xu J, et al. Postharvest Biology and Technology, 2018, 143: 112. [本文引用:1]
[19] XU Ying-chao, WANG Xiang-you, YIN Xiang, et al(许英超, 王相友, 印祥, 等). Transactions of the Chinese Society for Agricultural Machinery(农业机械学报), 2018, 49(2): 339. [本文引用:1]
[20] CHU Xiao-li(褚小立). Near-infrared Spectroscopy Analytical Technology Practical Hand book(近红外光谱分析技术实用手册). Beijing: China Machine Press(北京: 机械工业出版社), 2016. [本文引用:1]
[21] Zhang C, Guo C, Liu F, et al. Journal of Food Engineering, 2016, 179: 11. [本文引用:1]