基于高光谱成像的贵长猕猴桃硬度快速预测

尚静; 冯树南; 谭涛; 吴美芝; 陈海江; 孟庆龙

doi:10.13386/j.issn1002-0306.2022060168

基于高光谱成像的贵长猕猴桃硬度快速预测

尚静^{1, 2,},
冯树南¹,
谭涛¹,
吴美芝¹,
陈海江¹,
孟庆龙^{1, 2, ,}

1.
贵阳学院食品与制药工程学院，贵州贵阳 550005
2.
贵州省农产品无损检测工程研究中心，贵州贵阳 550005

基金项目: 贵州省科技计划项目（黔科合基础[2020]1Y270）；贵阳学院专项资金（GYU-KY-[2022]）。

详细信息

作者简介:
尚静（1988−），女，硕士，讲师，研究方向：农产品无损检测，E-mail：shji0124@163.com

通讯作者:
孟庆龙（1989−），男，博士，副教授，研究方向：农产品无损检测，E-mail：scumql@163.com

中图分类号: O657.3；S663.4
计量
- 文章访问数: 191
- HTML全文浏览量: 60
- PDF下载量: 23
出版历程
- 收稿日期: 2022-06-19
- 网络出版日期: 2023-01-11
- 刊出日期: 2023-03-14

Rapid Prediction for the Firmness of Guichang Kiwifruit by Hyperspectral Imaging

1.
Food and Pharmaceutical Engineering Institute, Guiyang University, Guiyang 550005, China
2.
Research Center of Nondestructive Testing for Agricultural Products of Guizhou Province, Guiyang 550005, China

摘要

摘要: 为了提高猕猴桃硬度预测的效率。应用可见/近红外（390~1030 nm）高光谱成像系统获取贵长猕猴桃的高光谱图像信息，并提取每个样品感兴趣区域的反射光谱，采用标准正态变换消除原始反射光谱中的噪声干扰；基于竞争性自适应重加权（competitive adaptive reweighted sampling，CARS）与连续投影算法筛选特征变量，建立误差反向传播神经网络和多元线性回归（multi linear regression，MLR）模型检测贵长猕猴桃硬度。结果表明：应用CARS从256个变量中筛选了35个特征变量，运算效率提升了近11倍（即运算时间从5.84 s降到了0.54 s）；构建的CARS-MLR检测模型具有较大的r_c=0.95和r_p=0.92，较小的RMSEC=1.65 kg/cm²和RMSEP=1.99 kg/cm²，且RPD值（2.47）大于2，表明CARS-MLR模型具有非常好的检测性能，利用高光谱成像技术以及化学计量学可以实现贵长猕猴桃硬度的快速无损检测。
- 贵长猕猴桃 /
- 硬度 /
- 高光谱成像 /
- 化学计量学 /
- 无损检测
Abstract: The work aimed to improve the prediction efficiency in rapid determination of the firmness of kiwifruit. A Vis/NIR (390~1030 nm) hyperspectral imaging system was applied to obtain the images of Guichang kiwifruit, and the reflective spectrum in the regions of interest on each sample was acquired. Noise from original reflective spectrum was reduced by the standard normal variation method. The competitive adaptive reweighted sampling (CARS) and the successive projection algorithm were applied to select feature variables. Finally an error back propagation neural network and a multi linear regression (MLR) model were constructed to predict the firmness of kiwifruit. A total of 35 feature variables were selected by CARS from 256 variables. The working efficiency of the final prediction model was improved by 11-fold, with the runtime dropped from 5.84 s to 0.54 s. Overall, the CARS-MLR model showed a relatively good detection capability (r_c=0.95, r_p=0.92, RMSEC=1.65 kg/cm², RMSEP=1.99 kg/cm², RPD above 2). This study demonstrated the application potential of the nondestructive hyperspectral imaging technology for fast determination of the firmness of kiwifruit.
- Guichang kiwifruit /
- firmness /
- hyperspectral imaging /
- chemometrics /
- nondestructive detection

HTML全文

贵长猕猴桃产于贵州省修文县，其果肉细嫩而且多浆、果汁酸甜可口，受到广大消费者的喜爱^[1-3]。猕猴桃属于后熟水果，采收后需要贮藏在室温条件下继续熟化，待完全熟透后才可以食用。但是假如食用过早，猕猴桃果会因生硬而影响口感，如果放置时间太久，猕猴桃果就会太软而易腐烂。猕猴桃的硬度是判别其成熟度和品质的重要指标，我国水果硬度的传统测定方法（利用果实硬度计的探头压入果肉中测定）虽然准确性相对较高，却拥有破坏检测对象、效率低、费时等缺陷^[4]；而依靠经验主观判断水果成熟度和品质就无统一标准，导致结果存在随意性，对品质的分选就不准确。因此，在水果品质的分选方面，无论是传统的检测方法还是个人经验主观判断都造成了水果品质参差不齐，难以体现我国水果的市场竞争优势。由此可见，建立一套快速无损的猕猴桃果实的硬度检测方法对于监控其成熟度和品质、指导市场价格是具有实际意义的。

基于精密光学检测仪器的研发和化学计量学的飞速发展，以及近红外光谱技术与高光谱成像技术具有的可实现无损检测的优势^[5]。国内外科研学者应用近红外光谱技术以及高光谱成像技术开展了关于草莓^[6]、猕猴桃^[7-8]、橙子^[9]、樱桃^[10-11]、葡萄^[12]和蓝莓^[13]等诸多水果内部品质的分析与研究。高升等^[14]采用可见/近红外漫反射光谱技术实现了红提葡萄硬度的无损检测，获得的预测集相关系数：0.93，均方根误差：8.36 N。刘卉等^[15]利用近红外光谱技术结合偏最小二乘法实现了冷藏期间“华优”猕猴桃硬度的无损检测，其校正集和预测集相关系数分别为0.96和0.85。而高光谱成像技术是一种集光谱信息与图像信息于一体的无损检测技术，能全面反映水果的内外部品质信息^[16]。Xie等^[17]使用高光谱成像结合特征变量筛选方法构建了检测香蕉硬度的偏最小二乘模型，获得的预测集的决定系数仅为0.76。Li等^[18]基于高光谱成像构建了评价李果实品质的多种检测模型。邵园园等^[19]采用高光谱成像技术结合人工神经网络实现了不同成熟度肥城桃可溶性固形物含量和硬度的高光谱检测。但是，基于高光谱成像以及化学计量学检测贮藏期内贵长猕猴桃硬度的研究与分析还未见报道，而且高光谱技术应用在水果品质无损检测领域面临的主要问题是数据量冗余。因此，对获取的高光谱数据进行降维处理尤为必要。

本文通过高光谱成像以及化学计量学实现贵长猕猴桃硬度的检测。运用竞争性自适应重加权与连续投影算法筛选特征变量，基于特征变量和测得的硬度参考值构建检测贵长猕猴桃硬度的误差反向传播神经网络和多元线性回归模型，依据预测效果得到最优模型。以期为水果成熟度以及内部品质的快速检测提供一定的理论基础。

1. 材料与方法

1.1 材料与仪器

贵长猕猴桃　来自贵州省修文县龙关口猕猴桃果园。从不同的果树上采摘并挑选没有病虫害且没有损伤的猕猴桃果子，猕猴桃果子摘取后马上运送到实验室，并贮藏在温度为20±2 ℃的实验室环境下，每隔3 d取40个猕猴桃样品进行测量，共进行4次实验，总计160个样品。每次实验前用软纸小心地除掉猕猴桃果子表面的灰尘，逐一编号获取高光谱图像并测定硬度。

TA.XT.Plus质构仪　英国SMS公司；GaiaField-F-V10高光谱成像系统　江苏双利合谱科技有限公司。图1展示了高光谱成像系统结构框图，其中，光谱波长范围：390~1030 nm，CCD相机的曝光时间：12.6 ms，猕猴桃样品距离镜头大约：40 cm，仪器其他参数见参考文献[20]。

图 1 高光谱成像系统框图

Figure 1. Schematic diagram of hyperspectral imaging system

下载: 全尺寸图片幻灯片

1.2 实验方法

1.2.1 高光谱图像获取与校正

将猕猴桃样品按编号顺序依次放在电动平移台上，挨个扫描样品，获取猕猴桃样品的高光谱图像。为消除因光源强度分布不均匀和猕猴桃样品的形状差异等带来的噪声干扰，采集完所有猕猴桃样品的高光谱图像后，在相同的系统参数下获取全白以及全黑的标定图像，在SpecVIEW软件中对所有猕猴桃样品的原始高光谱图像进行校正，公式如下：

${{\rm{R}}_{{\text{ref}}}}{{ = }}\frac{{{{\rm{R}}_{{\text{original}}}} - {{\rm{R}}_{{\text{black}}}}}}{{{{\rm{R}}_{{\text{white}}}} - {{\rm{R}}_{{\text{black}}}}}}$

(1)

式中：R_ref为校正后的图像；R_original为原始的图像；R_white和R_black分别为全白和全黑的标定图像。

1.2.2 猕猴桃硬度测定方法

在完成猕猴桃高光谱图像采集后，采用质构仪测定猕猴桃硬度，质构仪的探头直径：2 mm，探头压入猕猴桃果实的深度：10 mm，仪器测试时的速度：1 mm/s，测前和测后的速度都是2 mm/s。硬度测定前先将猕猴桃样品测定部位去皮，每个样品分别测定三个位置（如图2所示），分别将每一个猕猴桃三次测定值的平均值作为此样品的硬度参考值。

图 2 硬度测量位置

Figure 2. Firmness measured position

下载: 全尺寸图片幻灯片

1.2.3 建模方法

1.2.3.1 反射光谱提取及预处理

采用ENVI 5.4软件从校正后的高光谱图像中选取猕猴桃样品硬度测量的三个位置（如图2所示）作为感兴趣区域（regions of interest，ROIs），并获取反射光谱，将ROIs区域的反射率平均值作为此猕猴桃样品的原始的反射光谱。通过前期探究多种光谱预处理方法（标准正态变换、多元散射校正、二阶导数等）对原始反射光谱的预处理效果，得出标准正态变换是相对较优的预处理方法^[20-21]，因此本文采用标准正态变换（standard normal variation，SNV）进一步消除反射光谱中的噪声信号^[22]，计算公式如下：

$\left\{ \begin{gathered} {{\rm{Z}}_{{\rm{ij}}}} = \frac{{{{\rm{x}}_{{\rm{ij}}}} - \overline {{{\rm{x}}_{\rm{i}}}} }}{{{{\rm{S}}_{\rm{i}}}}} \\ {{\rm{\overline x}} _{\rm{i}}} = \frac{1}{{\rm{p}}}\sum\limits_{{\rm{j}} = 1}^{\rm{p}} {{{\rm{x}}_{{\rm{ij}}}}} \\ {{\rm{S}}_{\rm{i}}} = \sqrt {\frac{1}{{{\rm{p}} - 1}}\sum\limits_{{\rm{j}} = 1}^{\rm{p}} {{{({{\rm{x}}_{{\rm{ij}}}} - {{{\rm{\overline x}} }_{\rm{i}}})}^2}} } \\ \end{gathered} \right.$

(2)

式中：Z_ij代表标准化的光谱；x_ij代表原始的光谱； ${\overline {\rm{x}} _i}$ 是原始光谱平均值；S_i表示标准偏差，i=1，2，···，n；j=1，2，···，p；n为猕猴桃样本个数；p为光谱点个数。

1.2.3.2 特征变量筛选

运用竞争性自适应重加权（competitive adaptive reweighted sampling，CARS）与连续投影算法（successive projection algorithm，SPA）筛选特征变量，两种算法详见参考文献[20]。

1.2.3.3 预测模型及评价方法

实验基于误差反向传播（error back propagation，BP）神经网络^[23]与多元线性回归（multi linear regression，MLR）分别构建了检测贵长猕猴桃硬度的快速无损检测模型。MLR回归方程如下：

$\left. \begin{gathered} {{\rm{Y}}_1} = {\beta _0}+{\beta _1}{{\rm{X}}_{11}}+{\beta _2}{{\rm{X}}_{12}}+\cdots+{\beta _{\rm{m}}}{{\rm{X}}_{1{\rm{m}}}}+{\varepsilon _1} \\ {{\rm{Y}}_2} = {\beta _0}+{\beta _1}{{\rm{X}}_{21}}+{\beta _2}{{\rm{X}}_{22}}+\cdots+{\beta _{\rm{m}}}{{\rm{X}}_{2{\rm{m}}}}+{\varepsilon _2} \\ \cdots \\ {{\rm{Y}}_{\rm{n}}} = {\beta _0}+{\beta _1}{{\rm{X}}_{{\rm{n}}1}}+{\beta _2}{{\rm{X}}_{{\rm{n}}2}}+\cdots+{\beta _{\rm{m}}}{{\rm{X}}_{{\rm{nm}}}}+{\varepsilon _{\rm{n}}} \\ \end{gathered} \right\}$

(3)

式中：Y为硬度值；X_nm表示特征变量；β₀，β₁，β₂，···，β_m为待定系数；ε₁，ε₂，ε₃，···，ε_n为n个独立且服从同一正态分布的随机误差变量；m为特征变量数；n为猕猴桃样本个数。

通过建模集与预测集的相关系数r_c与r_p，均方根误差RMSEC与RMSEP，以及剩余预测偏差RPD衡量模型性能优劣的关键指标。若RPD<1.4，证明构建的检测模型无法完成预测；1.4≤RPD<1.8，则模型仅可粗略地预测；1.8≤RPD<2.0，证明构建的模型可完成相对较好的预测；RPD≥2，证明模型可以完成非常好的预测^[24-25]。

${\text{RPD}} = \frac{{{\text{SD}}}}{{{\text{RMSEP}}}}$

(4)

式中：SD代表预测集中猕猴桃样本硬度参考值的标准偏差。

1.3 数据处理

利用软件Spectral SENS采集所有猕猴桃样品的高光谱图像，采集获得的猕猴桃高光谱图像数据在ENVI 5.4和MATLAB R2016b软件中进行处理与分析。

2. 结果与分析

2.1 硬度测定结果

在建立无损检测模型时，通过获取的反射光谱以及测定的硬度参考值将所有贵长猕猴桃样品分成建模集和预测集。通过光谱-理化值共生距离方法^[26]从160个猕猴桃样品选择120个样品用于建模集，其余的40个用于预测集。表1列出了建模集和预测集中贵长猕猴桃硬度参考值统计结果，从表1明显得出，建模集中硬度的参考值范围大于预测集中硬度的参考值范围，这样分成的建模集和预测集有利于构建准确、稳定的检测模型。

表 1 猕猴桃硬度值

Table 1. Values of kiwifruit firmness

样本集	样本数	硬度（kg/cm²）
样本集	样本数	最小值	最大值	平均值	标准差
建模集	120	1.99	18.80	10.54	5.31
预测集	40	2.02	18.76	13.27	4.92

下载: 导出CSV

| 显示表格

2.2 猕猴桃反射光谱及预处理

由于原始反射光谱中包含一些噪声，为了提升贵长猕猴桃无损检测模型的精准度及稳定性，运用SNV进一步消除原始的反射光谱中的噪声，所有猕猴桃样品的原始的反射光谱以及经过SNV预处理后的相对反射光谱如图3所示。从图3可以看出，波长675 nm附近的吸收峰是由叶绿素的吸收引起的，体现了猕猴桃的颜色信息，波长980 nm附近的吸收峰是由水分的吸收引起的，反映了猕猴桃中的水分含量。对比图3a与图3b可以发现，预处理后的相对反射光谱曲线要比原始的反射光谱曲线相对平滑些，证明SNV预处理去掉了部分噪声以及背景干扰信息。

图 3 猕猴桃的反射光谱

注：a.原始的光谱；b.SNV预处理后的光谱。

Figure 3. Reflectance spectra of kiwifruit

下载: 全尺寸图片幻灯片

2.3 特征变量筛选

2.3.1 采用SPA筛选特征变量

应用SPA筛选特征变量时，基于均方根误差来筛选出相对较优的特征变量。图4a给出了均方根误差随着SPA中变量数的改变规律，从图4a可以明显看出，当模型中包含的变量个数小于22时，均方根误差随着变量数的增多而减小，而当模型中包含的变量个数大于22时，均方根误差趋于稳定，因此将这22个变量（占256个总变量的8.59%）确定为特征变量，这22个变量分别为：427.5、439.2、451.0、500.7、522.2、541.4、604.4、673.2、703.0、773.3、821.6、847.2、873.0、878.2、885.9、888.5、924.9、932.7、993.2、1004.0、1012.0、1020.0，单位均为nm。图4b给出了采用SPA筛选的22个特征变量分布情况。

图 4 SPA筛选特征变量结果

注：a.均方根误差随着变量数的改变规律；b.SPA筛选的特征变量。

Figure 4. Results of feature variables selected by SPA

下载: 全尺寸图片幻灯片

2.3.2 采用CARS筛选特征变量

运用CARS筛选特征变量时，依据交叉验证均方根误差的最小值来筛选出相对较优的特征变量。图5a和图5b分别给出了交叉验证均方根误差和采样变量数随着采样次数的改变规律，从图5a和图5b可以观察到，第21次采样获得的交叉验证均方根误差最小（2.28），该变量集中包括了35个特征变量（占256个总变量的13.67%），这35个变量分别为：427.5、429.8、448.6、455.7、484.1、601.9、616.6、648.5、651.0、670.7、673.2、685.6、695.5、758.2、763.2、773.3、796.1、803.8、806.3、829.3、847.2、862.7、865.3、873.0、896.3、909.3、922.3、932.7、940.6、982.6、987.9、990.5、1004、1012、1017，单位均为nm。图5c给出了应用CARS筛选的35个特征变量分布情况。对比图5c与图4b可以看出，CARS筛选的特征变量相对集中，主要位于波段427.5~484.1、601.9~695.5、758.2~940.6 nm以及982.6~1027 nm之间，而SPA筛选的特征变量要离散一些。

图 5 CARS筛选特征变量的结果

注：a.交叉验证均方根误差随着采样次数的改变规律；b.采样变量数随着采样次数的改变规律；c.CARS筛选的特征变量。

Figure 5. Results of feature variables selected by CARS

下载: 全尺寸图片幻灯片

2.4 基于特征光谱的建模结果

基于CARS和SPA筛选出的特征变量以及建模集和预测集中硬度的参考值，分别构建检测贵长猕猴桃硬度的BP神经网络和MLR快速无损检测模型。表2列出了基于特征变量构建的检测模型对贵长猕猴桃硬度的检测结果，从表2可以得出，基于SPA筛选的特征变量构建的SPA-BP和SPA-MLR两种无损检测模型仅可粗略地预测贵长猕猴桃硬度，其1.4≤RPD<1.8；而基于CARS筛选的特征变量构建的CARS-BP和CARS-MLR两种无损检测模型可以较好地实现贵长猕猴桃硬度的预测，其RPD分别为2.02和2.47，表明CARS特征变量筛选方法比SPA方法相对较好。对比CARS-BP和CARS-MLR检测模型可以看出，构建的CARS-MLR检测模型拥有最大的r_c=0.95和r_p=0.92，以及较小的RMSEC=1.65 kg/cm²和RMSEP=1.99 kg/cm²，且RPD（2.47）大于2，相关文献报道的猕猴桃硬度定量检测模型的r_p为0.85^[15]，小于本研究构建的CARS-MLR检测模型的r_p，进一步证实本研究构建的猕猴桃硬度检测模型具有更好的预测效果。同时，运用CARS从256个变量中筛选出35个特征变量，大大提高了检测模型的运算效率，即对比基于全光谱与特征变量构建的检测模型运算效率，得出基于特征变量构建的模型的运算效率比全光谱模型的运算效率提升了近11倍（运算时间从5.84 s降到了0.54 s）。综上所述，基于特征变量构建的CARS-MLR检测模型对检测贵长猕猴桃硬度拥有非常好的预测性能，且运算效率得到了较明显提升。图6展示了CARS-MLR模型对贵长猕猴桃硬度的检测结果。

表 2 猕猴桃硬度BP和MLR检测模型

Table 2. BP and MLR detection models for the firmness of kiwifruit

模型	特征变量筛选方法	变量数	建模集		预测集
模型	特征变量筛选方法	变量数	r_c	RMSEC （kg/cm²）	r_p	RMSEP （kg/cm²）	RPD
BP	SPA	22	0.80	3.19	0.79	2.97	1.66
BP	CARS	35	0.93	2.11	0.87	2.43	2.02
MLR	SPA	22	0.81	3.08	0.80	2.95	1.67
MLR	CARS	35	0.95	1.65	0.92	1.99	2.47

下载: 导出CSV

| 显示表格

图 6 猕猴桃硬度参考值与预测值

Figure 6. Reference and predicted values of the firmness of kiwifruit

下载: 全尺寸图片幻灯片

3. 结论与讨论

本研究采用高光谱成像技术结合化学计量学，得到了一种能够快速、无损预测贵长猕猴桃硬度的方法。其中，为消除原始反射光谱中的部分噪声和背景干扰信息，运用SNV方法对原始光谱进行预处理，得出预处理后的相对反射光谱曲线要比原始的反射光谱曲线相对平滑些。基于SPA和CARS筛选特征变量，得出构建的SPA-BP和SPA-MLR两种无损检测模型仅可粗略地预测贵长猕猴桃硬度，其1.4≤RPD<1.8，而构建的CARS-BP和CARS-MLR两种无损检测模型可以较好地实现贵长猕猴桃硬度的预测，其RPD分别为2.02和2.47，证明CARS特征变量筛选方法比SPA方法相对较好。基于特征变量构建的CARS-MLR检测模型拥有较大的r_c=0.95与r_p=0.92，以及较小的RMSEC=1.65 kg/cm²和RMSEP=1.99 kg/cm²，且RPD大于2，证明CARS-MLR模型可以非常好的预测贵长猕猴桃硬度。相对于文献报道的基于全变量构建的猕猴桃硬度无损检测模型^[15]，本研究采用CARS从全变量中筛选出35个特征变量，大大提高了模型的运算效率，即CARS-MLR模型的运算效率比全变量模型的运算效率提升了近11倍（运算时间从5.84 s降到了0.54 s）。为水果成熟度以及内部品质的快速检测提供一定的理论基础。此次实验样本仅采摘于一个猕猴桃果园，为了提高模型的适用性和稳定性，后期将进一步扩充和完善实验。

图 1 高光谱成像系统框图

Figure 1. Schematic diagram of hyperspectral imaging system

下载: 全尺寸图片幻灯片

图 2 硬度测量位置

Figure 2. Firmness measured position

下载: 全尺寸图片幻灯片

图 3 猕猴桃的反射光谱

注：a.原始的光谱；b.SNV预处理后的光谱。

Figure 3. Reflectance spectra of kiwifruit

下载: 全尺寸图片幻灯片

图 4 SPA筛选特征变量结果

注：a.均方根误差随着变量数的改变规律；b.SPA筛选的特征变量。

Figure 4. Results of feature variables selected by SPA

下载: 全尺寸图片幻灯片

图 5 CARS筛选特征变量的结果

注：a.交叉验证均方根误差随着采样次数的改变规律；b.采样变量数随着采样次数的改变规律；c.CARS筛选的特征变量。

Figure 5. Results of feature variables selected by CARS

下载: 全尺寸图片幻灯片

图 6 猕猴桃硬度参考值与预测值

Figure 6. Reference and predicted values of the firmness of kiwifruit

下载: 全尺寸图片幻灯片

表 1 猕猴桃硬度值

Table 1 Values of kiwifruit firmness

样本集	样本数	硬度（kg/cm²）
样本集	样本数	最小值	最大值	平均值	标准差
建模集	120	1.99	18.80	10.54	5.31
预测集	40	2.02	18.76	13.27	4.92

下载: 导出CSV

表 2 猕猴桃硬度BP和MLR检测模型

Table 2 BP and MLR detection models for the firmness of kiwifruit

模型	特征变量筛选方法	变量数	建模集		预测集
模型	特征变量筛选方法	变量数	r_c	RMSEC （kg/cm²）	r_p	RMSEP （kg/cm²）	RPD
BP	SPA	22	0.80	3.19	0.79	2.97	1.66
BP	CARS	35	0.93	2.11	0.87	2.43	2.02
MLR	SPA	22	0.81	3.08	0.80	2.95	1.67
MLR	CARS	35	0.95	1.65	0.92	1.99	2.47

下载: 导出CSV

参考文献(26)

[1]	王金华, 杜超, 梁晨, 等. 贵长猕猴桃多糖提取工艺及体外抗氧化功能[J]. 食品科学,2016,37(20):19−23. [WANG J H, DU C, LIANG C, et al. Extraction and antioxidant activity of polysaccharides from Guichang kiwifruit[J]. Food Science,2016,37(20):19−23. doi: 10.7506/spkx1002-6630-201620004
[2]	张承, 李明, 龙友华, 等. 采前喷施壳聚糖复合膜对猕猴桃软腐病的防控及其保鲜作用[J]. 食品科学,2016,37(22):274−281. [ZHANG C, LI M, LONG Y H, et al. Control of soft rot in kiwifruit by pre-harvest application of chitosan composite coating and its effect on preserving and improving kiwifruit quality[J]. Food Science,2016,37(22):274−281. doi: 10.7506/spkx1002-6630-201622042
[3]	靳政时, 牛犇, 刘瑞玲, 等. 干燥方式对猕猴桃果干品质的影响[J]. 食品工业科技,2022,43(24):62−71. [JIN Z S, NIU B, LIU R L, et al. Effects of drying methods on the quality of dried kiwifruit[J]. Science and Technology of Food Industry,2022,43(24):62−71.
[4]	曹念念, 刘强, 彭菁, 等. 基于近红外光谱技术的黄桃脆片可溶性固形物和硬度定量检测方法[J]. 食品机械,2021,37(3):51−57. [CAO N N, LIU Q, PENG J, et al. Study on quantitative detection of soluble solids and firmness of yellow peach chips by near-infrared spectroscopy[J]. Food Machinery,2021,37(3):51−57.
[5]	胡逸磊, 姜洪喆, 周宏平, 等. 水果成熟度近红外光谱及高光谱成像无损检测研究进展[J]. 食品工业科技,2021,42(20):377−383. [HU Y L, JIANG H Z, ZHOU H P, et al. Research progress on nondestructive detection of fruit maturity by near infrared spectroscopy and hyperspectral imaging[J]. Science and Technology of Food Industry,2021,42(20):377−383. doi: 10.13386/j.issn1002-0306.2020070074
[6]	DEVASSY B M, GEORGE S. Estimation of strawberry firmness using hyperspectral imaging: A comparison of regression models[J]. Journal of Spectral Imaging,2021,10(a3):1−14.
[7]	许丽佳, 陈铭, 王玉超, 等. 高光谱成像的猕猴桃糖度无损检测方法[J]. 光谱学与光谱分析,2021,41(7):2188−2195. [XU L J, CHEN M, WANG Y C, et al. Study on non-destructive detection method of kiwifruit sugar content based on hyperspectral imaging technology[J]. Spectroscopy and Spectral Analysis,2021,41(7):2188−2195.
[8]	HU W H, SUN D W, BLASCO J. Rapid monitoring 1-MCP-induced modulation of sugars accumulation in ripening 'Hayward' kiwifruit by Vis/NIR hyperspectral imaging[J]. Postharvest Biology and Technology,2017,125:168−180. doi: 10.1016/j.postharvbio.2016.11.001
[9]	ZHANG H L, ZHAN B S, PAN F, et al. Determination of soluble solids content in oranges using visible and near infrared full transmittance hyperspectral imaging with comparative analysis of models[J]. Postharvest Biology and Technology,2020,163(111148):1−9.
[10]	LI X L, WEI Y Z, XU J, et al. SSC and pH for sweet assessment and maturity classification of harvested cherry fruit based on NIR hyperspectral imaging technology[J]. Postharvest Biology and Technology,2018,143:112−118. doi: 10.1016/j.postharvbio.2018.05.003
[11]	PULLANAGARI R R, LI M. Uncertainty assessment for firmness and total soluble solids of sweet cherries using hyperspectral imaging and multivariate statistics[J]. Journal of Food Engineering,2020,289(110177):1−12.
[12]	蔡正云, 吴龙国, 王菁, 等. 宁夏赤霞珠葡萄水分含量的高光谱无损检测研究[J]. 食品工业科技,2017,38(2):79−88. [CAI Z Y, WU L G, WANG J, et al. Non-destructive determination of moisture composition in Ningxia wine grapes based on visible near-infrared hyperspectral imaging technique[J]. Science and Technology of Food Industry,2017,38(2):79−88. doi: 10.13386/j.issn1002-0306.2017.02.007
[13]	LEIVA-VALENZUELA G A, LU R, AGUILERA J M. Assessment of internal quality of blueberries using hyperspectral transmittance and reflectance images with whole spectra or selected wavelengths[J]. Innovative Food Science & Emerging Technologies,2014,24:2−13.
[14]	高升, 王巧华, 施行, 等. 便携式红提葡萄多品质可见/近红外检测仪设计与试验[J]. 农业机械学报,2021,52(2):308−319. [GAO S, WANG Q H, SHI X, et al. Design and test of portable red globe grape extraction multi-quality visible/near infrared detector[J]. Transactions of the Chinese Society for Agricultural Machinery,2021,52(2):308−319. doi: 10.6041/j.issn.1000-1298.2021.02.030
[15]	刘卉, 郭文川, 岳绒. 猕猴桃硬度近红外漫反射光谱无损检测[J]. 农业机械学报,2011,42(3):145−149. [LIU H, GUO W C, YUE R. Non-destructive detection of kiwifruit firmness based on near-infrared diffused spectroscopy[J]. Transactions of the Chinese Society for Agricultural Machinery,2011,42(3):145−149.
[16]	马佳佳, 王克强. 水果品质光学无损检测技术研究进展[J]. 食品工业科技,2021,42(23):427−437. [MA J J, WANG K Q. Research progress of optical nondestructive testing technology for fruit quality[J]. Science and Technology of Food Industry,2021,42(23):427−437.
[17]	XIE C Q, CHU B Q, HE Y. Prediction of banana color and firmness using a novel wavelengths selection method of hyperspectral imaging[J]. Food Chemistry,2018,245:132−140. doi: 10.1016/j.foodchem.2017.10.079
[18]	LI B, COBO-MEDINA M, LECOURT J, et al. Application of hyperspectral imaging for nondestructive measurement of plum quality attributes[J]. Postharvest Biology and Technology,2018,141:8−15. doi: 10.1016/j.postharvbio.2018.03.008
[19]	邵园园, 王永贤, 玄冠涛, 等. 基于高光谱成像的肥城桃品质可视化分析与成熟度检测[J]. 农业机械学报,2020,51(8):344−350. [SHAO Y Y, WANG Y X, XUAN G T, et al. Visual detection of SSC and firmness and maturity prediction for Feicheng peach by using hyperspectral imaging[J]. Transactions of the Chinese Society for Agricultural Machinery,2020,51(8):344−350. doi: 10.6041/j.issn.1000-1298.2020.08.038
[20]	MENG Q L, SHANG J, HUANG R S, et al. Determination of soluble solids content and firmness in plum using hyperspectral imaging and chemometric algorithms[J]. Journal of Food Process Engineering,2021,44(1):1−9.
[21]	孟庆龙, 尚静, 黄人帅, 等. 基于主成分回归的猕猴桃可溶性固形物无损检测[J]. 包装工程,2021,42(3):19−24. [MENG Q L, SHANG J, HUANG R S, et al. Nondestructive detection for soluble solids content of kiwifruits based on principal component regression[J]. Packaging Engineering,2021,42(3):19−24.
[22]	倪力军, 张立国. 基础化学计量学及其应用[M]. 上海: 华东理工大学出版社, 2011: 67−68. NI L J, ZHANG L G. Basic chemometrics and its application[M]. Shanghai: East China University of Science and Technology Press, 2011: 67−68.
[23]	郭文川, 董金磊. 高光谱成像结合人工神经网络无损检测桃的硬度[J]. 光学精密工程,2015,23(6):1530−1537. [GUO W C, DONG J L. Nondestructive detection on firmness of peaches based on hyperspectral imaging and artificial neural networks[J]. Optics and Precision Engineering,2015,23(6):1530−1537. doi: 10.3788/OPE.20152306.1530
[24]	NICOLAÏ B, BEULLENS K, BOBELYN E, et al. Nondestructive measurement of fruit and vegetable quality by means of NIR spectroscopy: A review[J]. Postharvest Biology and Technology,2007:99−118.
[25]	于雷, 章涛, 朱亚星, 等. 基于IRIV算法优选大豆叶片高光谱特征波长变量估测SPAD值[J]. 农业工程学报,2018,34(16):148−154. [YU L, ZHANG T, ZHU Y X, et al. Determination of soybean leaf SPAD value using characteristic wavelength variables preferably selected by IRIV algorithm[J]. Transactions of the Chinese Society of Agricultural Engineering,2018,34(16):148−154. doi: 10.11975/j.issn.1002-6819.2018.16.019
[26]	GALVÃO R K H, ARAUJO M C U, JOSÉ G E, et al. A method for calibration and validation subset partitioning[J]. Talanta,2005,67(4):736−740. doi: 10.1016/j.talanta.2005.03.025

施引文献(10)

期刊类型引用(5)

1.	刘丽，杨天意，董聪颖，石彩云，司鹏，魏志峰，高登涛. 米良1号猕猴桃贮藏过程中糖酸含量估测模型的构建. 果树学报. 2025(01): 216-226 . 百度学术
2.	王加龙，马坤，高鹏，朱金芳，张平，黄凡. 基于可见-近红外光谱技术研发便携式贝贝南瓜品质无损检测仪. 食品科学. 2025(06): 254-262 . 百度学术
3.	翟依晨，张林，祝佳佳，龙艳，杨鑫睿，赵志博. 猕猴桃软腐病无损高光谱诊断技术研究. 山地农业生物学报. 2025(02): 78-86 . 百度学术
4.	闫鼎，张义志，程森，蔡宪杰，董祥洲，杨悦章，岳耀稳，王大彬，林润英. 烤后烟叶不同部位高光谱特征分析及判别模型构建. 湖北农业科学. 2024(08): 140-146 . 百度学术
5.	金科，郭志强，曾云流，丁港. 结合三维卷积神经网络和Haar小波滤波器的高光谱猕猴桃品质分类. 激光与光电子学进展. 2023(20): 49-58 . 百度学术

其他类型引用(5)

资源附件(0)

图(6) / 表(2)

计量

文章访问数: 191
HTML全文浏览量: 60
PDF下载量: 23
被引次数: 10

1. 材料与方法
1.1 材料与仪器
1.2 实验方法
1.2.1 高光谱图像获取与校正
1.2.2 猕猴桃硬度测定方法
1.2.3 建模方法
1.2.3.1 反射光谱提取及预处理
1.2.3.2 特征变量筛选
1.2.3.3 预测模型及评价方法
1.3 数据处理
2. 结果与分析
2.1 硬度测定结果
2.2 猕猴桃反射光谱及预处理
2.3 特征变量筛选
2.3.1 采用SPA筛选特征变量
2.3.2 采用CARS筛选特征变量
2.4 基于特征光谱的建模结果
3. 结论与讨论

1. 材料与方法
1.1 材料与仪器
1.2 实验方法
1.2.1 高光谱图像获取与校正
1.2.2 猕猴桃硬度测定方法
1.2.3 建模方法
1.2.3.1 反射光谱提取及预处理
1.2.3.2 特征变量筛选
1.2.3.3 预测模型及评价方法
1.3 数据处理
2. 结果与分析
2.1 硬度测定结果
2.2 猕猴桃反射光谱及预处理
2.3 特征变量筛选
2.3.1 采用SPA筛选特征变量
2.3.2 采用CARS筛选特征变量
2.4 基于特征光谱的建模结果
3. 结论与讨论

参考文献(26)

施引文献

资源附件(0)

基于高光谱成像的贵长猕猴桃硬度快速预测

作者简介: 尚静（1988−），女，硕士，讲师，研究方向：农产品无损检测，E-mail：shji0124@163.com

通讯作者: 孟庆龙（1989−），男，博士，副教授，研究方向：农产品无损检测，E-mail：scumql@163.com

计量

出版历程