高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

基于高光谱成像的枇杷果实品质检测

,  ,  ,  ,  , 

上一篇

下一篇

吴习宇, 曾凯芳, 郭启高, 等. 基于高光谱成像的枇杷果实品质检测[J]. 西南大学学报（自然科学版）, 2022, 44(9): 52-60. doi: 10.13718/j.cnki.xdzk.2022.09.006
引用本文: 吴习宇, 曾凯芳, 郭启高, 等. 基于高光谱成像的枇杷果实品质检测[J]. 西南大学学报（自然科学版）, 2022, 44(9): 52-60. doi: 10.13718/j.cnki.xdzk.2022.09.006
shu
WU Xiyu, ZENG Kaifang, GUO Qigao, et al. Detection of Loquat Fruit Quality Based on Hyperspectral Imaging[J]. Journal of Southwest University Natural Science Edition, 2022, 44(9): 52-60. doi: 10.13718/j.cnki.xdzk.2022.09.006
Citation: WU Xiyu, ZENG Kaifang, GUO Qigao, et al. Detection of Loquat Fruit Quality Based on Hyperspectral Imaging[J]. Journal of Southwest University Natural Science Edition, 2022, 44(9): 52-60. doi: 10.13718/j.cnki.xdzk.2022.09.006
shu

基于高光谱成像的枇杷果实品质检测

  • 1.

    西南大学 食品科学学院, 重庆 400716

  • 2.

    西南大学 食品贮藏与物流研究中心, 重庆 400716

  • 3.

    西南大学 园艺园林学院, 重庆 400716

  • 基金项目: “十三五”国家重点研发项目子课题(2019YFD1000905); 西南大学博士基金项目(swu118083)
详细信息
    作者简介:

    吴习宇，博士，讲师，主要从事农产品无损检测、农产品贮藏保鲜的研究 .

    通讯作者: 徐丹，博士，教授

  • 中图分类号: Q946;TS255

Detection of Loquat Fruit Quality Based on Hyperspectral Imaging

  • 1.

    College of Food Science, Southwest University, Chongqing 400716, China

  • 2.

    Food Storage and Logistics Research Center, Southwest University, Chongqing 400716, China

  • 3.

    College of Horticulture and Landscape, Southwest University, Chongqing 400716, China

  • 摘要: 采用高光谱成像技术(HSI)在可见/近红外(363~1 026 nm)区域检测枇杷果实的可溶性固形物(SSC)和硬度，并判断其成熟度，以实现枇杷果实品质的无损检测和分级分选. 利用蒙特卡洛法(MC)剔除异常样本，基于联合X-Y距离(SPXY)进行建模集和预测集样本的划分，再采用竞争性自适应权重采样算法(CARS)和连续投影算法(SPA)选取特征波长，与全波段光谱(FS)比较，分别建立偏最小二乘回归(PLSR)模型. 结果显示，CARS-PLSR模型更优，CARS提取的SSC特征波长和硬度特征波长分别占总波长的8.52%和5.36%，枇杷果实中SSC和硬度的建模集相关系数Rc分别为0.981 7，0.970 7，预测集相关系数Rp分别为0.918 5，0.742 3，说明CARS能有效地对光谱进行降维，简化了数据处理过程. 枇杷果实SSC和硬度的变化与果实成熟度显著相关，建立判别偏最小二乘法(DPLS)成熟度预测模型，预测集总识别准确率为89.29%. 由此说明，高光谱成像技术可对枇杷品质进行有效检测，为枇杷果实的无损检测和分级分选提供了理论依据.
    Abstract: Hyperspectral imaging technology was used to detect the soluble solid content (SSC), firmness and maturity of loquat fruits in visible and near infrared (363~1 026 nm) region, in order to achieve nondestructive testing and sorting for loquat fruit. Abnormal samples were eliminated by using Monte Carlo (MC) method and the remaining samples were divided into calibration set and prediction set based on joint X-Y distances (SPXY). Then the competitive adaptive reweighted sampling (CARS) and successive projections algorithm (SPA) were used to select characteristic wavelengths, and compared with full spectrum (FS), the partial least squares regression (PLSR) models were established, respectively. The results showed that the CARS-PLSR models were better. The characteristic wavelengths of SSC and firmness extracted by CARS accounted for 8.52% and 5.36% of the total wavelength, respectively. The correlation coefficients Rc of calibration set for SSC and firmness of loquat fruit were 0.981 7 and 0.970 7, and the correlation coefficients Rp of prediction set were 0.918 5 and 0.742 3, respectively, which indicated that CARS effectively reduced the dimension of spectrum and simplified the data processing. There was a significant correlation between the changes of SSC and firmness of loquat fruit with fruit maturity. The discriminant partial least squares (DPLS) model of maturity was established and the total recognition accuracy of prediction set was 89.29%. Therefore, hyperspectral imaging technology can effectively detect the quality and maturity of loquat fruit, also can provide a theoretical basis for nondestructive rapid detection, sorting and grading of loquat fruit.
  • 加载中

  • 图 1  高光谱成像装置及示意图

    下载: 全尺寸图片 幻灯片

    图 2  枇杷果实光谱数据

    下载: 全尺寸图片 幻灯片

    图 3  蒙特卡洛偏最小二乘法检测结果

    下载: 全尺寸图片 幻灯片

    图 4  SSC和硬度直方图

    下载: 全尺寸图片 幻灯片

    图 5  CARS挑选SSC和硬度特征波长过程

    下载: 全尺寸图片 幻灯片

    图 6  SPA挑选SSC和硬度值特征波长过程

    下载: 全尺寸图片 幻灯片

    表 1  枇杷品质指标理化测定结果

    品质指标 建模集(样本数85) 预测集(样本数28)
    最小值 最大值 平均值 标准偏差 最小值 最大值 平均值 标准偏差
    SSC/°Brix 5.90 13.57 8.85 1.53 6.77 11.87 8.40 1.07
    硬度/(kg·cm-2) 4.07 7.17 5.90 0.47 5.53 6.63 5.99 0.26
    下载: 导出CSV

    表 2  不同特征波长下SSC的PLSR模型

    品质指标 特征波长提取方法 波长数 建模集 预测集
    Rc RMSEC Rp RMSEP
    SSC/°Brix FS(全光谱) 1 232 0.905 0 0.647 5 0.780 9 0.617 2
    CARS 105 0.981 7 0.294 2 0.918 5 0.373 8
    SPA 12 0.689 5 1.102 6 0.504 0 0.903 4
    下载: 导出CSV

    表 3  不同特征波长下硬度值的PLSR模型

    品质指标 特征波长提取方法 波长数 建模集 预测集
    Rc RMSEC Rp RMSEP
    硬度/(kg·cm-2) FS(全光谱) 1 232 0.805 1 0.278 6 0.728 8 0.202 8
    CARS 66 0.970 7 0.113 5 0.742 3 0.165 2
    SPA 10 0.690 6 0.409 3 0.630 2 0.298 9
    下载: 导出CSV

    表 4  成熟度预测的样本划分及预测结果

    成熟度 校正集样本数/个 预测集样本数/个 预测集误判数/个 预测集判别正确率/%
    八成熟及以上 33 6 2 66.67
    七成熟 52 22 1 95.45
    合计 85 28 3 89.29
    下载: 导出CSV
  • [1] SULTAN M Z. Loquat (Eriobotrya Japonica Lindl. ), Fruit and Vegetable Phytochemicals: Chemistry and Human Health[M]. 2nd Edition. West Sussex: John Wiley & Sons Ltd, 2018.
    [2] ALOS E, MARTINEZ-FUENTES A, REIG C, et al. Involvement of Ethylene in Color Changes and Carotenoid Biosynthesis in Loquat Fruit (Eriobotrya Japonica Lindl. Cv. Algerie)[J]. Postharvest Biology and Technology, 2019, 149: 129-138. doi: 10.1016/j.postharvbio.2018.11.022
    [3] 邓朝军, 吴琼, 许奇志, 等. 不同成熟度'贵妃'枇杷果实色泽与糖酸含量关系[J]. 热带作物学报, 2016, 37(9): 1747-1751. doi: 10.3969/j.issn.1000-2561.2016.09.015
    [4] 许奇志, 陈秀萍, 邓朝军, 等. 不同成熟度"新白8号"枇杷果实品质特性分析[J]. 中国南方果树, 2019, 48(2): 71-74, 81. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-FRUI201902017.htm
    [5] CAÑETE M L, HUESO J J, PINILLOS V, et al. Ripening Degree at Harvest Affects Bruising Susceptibility and Fruit Sensorial Traits of Loquat (Eriobotrya Japonica Lindl. )[J]. Scientia Horticulturae, 2015, 187: 102-107. doi: 10.1016/j.scienta.2015.03.008
    [6] CUBERO S, ALEIXOS N, MOLTÓ E, et al. Advances in Machine Vision Applications for Automatic Inspection and Quality Evaluation of Fruits and Vegetables[J]. Food and Bioprocess Technology, 2011, 4(4): 487-504. doi: 10.1007/s11947-010-0411-8
    [7] 高迎旺, 耿金凤, 饶秀勤. 果蔬采后内部损伤无损检测研究进展[J]. 食品科学, 2017, 38(15): 277-287. doi: 10.7506/spkx1002-6630-201715044
    [8] 张义志, 王瑞, 张伟峰, 等. 高光谱技术检测农产品成熟度研究进展[J]. 湖北农业科学, 2020, 59(12): 5-8, 12. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-HBNY202012002.htm
    [9] GOWEN A A, O'DONNELL C P, CULLEN P J, et al. Hyperspectral Imaging-an Emerging Process Analytical Tool for Food Quality and Safety Control[J]. Trends in Food Science & Technology, 2007, 18(12): 590-598.
    [10] RAJKUMAR P, WANG N, EIMASRY G, et al. Studies on Banana Fruit Quality and Maturity Stages Using Hyperspectral Imaging[J]. Journal of Food Engineering, 2012, 108(1): 194-200. doi: 10.1016/j.jfoodeng.2011.05.002
    [11] 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(15): 132-140.
    [12] MA T, XIA Y, INAGAKI T, et al. Rapid and Nondestructive Evaluation of Soluble Solids Content (SSC) and Firmness in Apple Using Vis-NIR Spatially Resolved Spectroscopy[J]. Postharvest Biology and Technology, 2021, 173: 111417. doi: 10.1016/j.postharvbio.2020.111417
    [13] 邵园园, 王永贤, 玄冠涛, 等. 基于高光谱成像的肥城桃品质可视化分析与成熟度检测[J]. 农业机械学报, 2020, 51(8): 344-350. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-NYJX202008038.htm
    [14] LI J B, RAO X Q, YING Y B. Detection of Common Defects on Oranges Using Hyperspectral Reflectance Imaging[J]. Computers and Electronics in Agriculture, 2011, 78(1): 38-48. doi: 10.1016/j.compag.2011.05.010
    [15] LIU F, JIN Z L, NAEEM M S, et al. Applying Near-Infrared Spectroscopy and Chemometrics to Determine Total Amino Acids in Herbicide-Stressed Oilseed Rape Leaves[J]. Food and Bioprocess Technology, 2011, 4(7): 1314-1321. doi: 10.1007/s11947-010-0445-y
    [16] 何凤芸, 杨晓东, 周胜灵, 等. 基于FT-MIR的三七粉末掺伪成分的定量检测模型研究[J]. 西南大学学报(自然科学版), 2022, 44(4): 54-61. doi: http://xbgjxt.swu.edu.cn/article/doi/10.13718/j.cnki.xdzk.2022.04.007
    [17] WOLD H. Estimation of Principal Components and Related Models by Iterative Least Squares[J]. Journal of Multivariate Analysis, 1966, 1: 391-420.
    [18] 吴习宇, 祝诗平, 王谦, 等. DPLS和SVM的掺假花椒粉近红外光谱定性鉴别[J]. 光谱学与光谱分析, 2018, 38(8): 2369-2373. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201808012.htm
    [19] 杨晓玉, 刘贵珊, 丁佳兴, 等. 灵武长枣VC含量的高光谱快速检测研究[J]. 光谱学与光谱分析, 2019, 39(1): 230-234. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201901041.htm
  • 加载中
图( 6) 表( 4)
计量
  • 文章访问数:  2457
  • HTML全文浏览数:  2457
  • PDF下载数:  167
  • 施引文献:  0
出版历程
  • 收稿日期:  2021-12-10
  • 刊出日期:  2022-09-20

基于高光谱成像的枇杷果实品质检测

    通讯作者: 徐丹，博士，教授
    作者简介: 吴习宇，博士，讲师，主要从事农产品无损检测、农产品贮藏保鲜的研究
  • 1. 西南大学 食品科学学院, 重庆 400716
  • 2. 西南大学 食品贮藏与物流研究中心, 重庆 400716
  • 3. 西南大学 园艺园林学院, 重庆 400716
  • 收稿日期:  2021-12-10
基金项目:  “十三五”国家重点研发项目子课题(2019YFD1000905); 西南大学博士基金项目(swu118083)

摘要: 采用高光谱成像技术(HSI)在可见/近红外(363~1 026 nm)区域检测枇杷果实的可溶性固形物(SSC)和硬度，并判断其成熟度，以实现枇杷果实品质的无损检测和分级分选. 利用蒙特卡洛法(MC)剔除异常样本，基于联合X-Y距离(SPXY)进行建模集和预测集样本的划分，再采用竞争性自适应权重采样算法(CARS)和连续投影算法(SPA)选取特征波长，与全波段光谱(FS)比较，分别建立偏最小二乘回归(PLSR)模型. 结果显示，CARS-PLSR模型更优，CARS提取的SSC特征波长和硬度特征波长分别占总波长的8.52%和5.36%，枇杷果实中SSC和硬度的建模集相关系数Rc分别为0.981 7，0.970 7，预测集相关系数Rp分别为0.918 5，0.742 3，说明CARS能有效地对光谱进行降维，简化了数据处理过程. 枇杷果实SSC和硬度的变化与果实成熟度显著相关，建立判别偏最小二乘法(DPLS)成熟度预测模型，预测集总识别准确率为89.29%. 由此说明，高光谱成像技术可对枇杷品质进行有效检测，为枇杷果实的无损检测和分级分选提供了理论依据.

English Abstract

    全文HTML

  • 开放科学(资源服务)标志码(OSID):

  • 枇杷(Eriobotrya japonica Lindl.)属于蔷薇科(Rosaceae)枇杷属常绿多年生植物，原产于中国南方及与云南省交界的老挝、越南等东南亚国家[1]. 枇杷果实含有丰富的矿物质、糖类和有机酸等营养成分，根据果肉颜色可分为红肉枇杷和白肉枇杷两种类型[2]. 枇杷属于呼吸非跃变型果实，成熟度低的枇杷糖分低、含酸量高，鲜食口感偏差; 八成熟以上果实含糖量持续增加，酸度降低[3-4]. 当前，枇杷采收后的品质分级主要通过肉眼可观测和可触摸的外观指标(颜色、大小、硬度、质量等)来进行[5]，但人的主观判断容易受很多因素的影响，效率和准确性都较差. 采用理化实验的方法可以对枇杷果实的糖度、酸度、硬度等进行较高精度的检测，但这些检测方法均为破坏性实验，需要耗费大量的化学试剂且费时费力. 而采用光学无损检测，建立数学模型后，仅通过机器扫描就能定量分析果实中的营养指标含量，这为我国未来水果成熟度的准确预测以及果品质量分级制度的建立带来了可能性.

    高光谱成像技术(hyperspectral imaging，HSI)利用很多窄的电磁波波段的电磁光谱以成像的形式获取物体特性的有关数据，把传统的二维成像技术和光谱技术有机地结合在一起，能同时分析样品的光谱信息和相应的空间信息[6-7]，具有“图谱合一”的特点[8]. 高光谱成像获取的原始图像是三维的，是一系列光波波长处的光学图像，图像像素的横坐标轴和纵坐标轴分别用x和y表示，光谱的波长信息以λ(z轴)表示[9]. 当前，一些学者研究了高光谱技术对香蕉[10-11]、苹果[12]、桃子[13]等水果品质及成熟度的无损检测，取得了较好的效果，但关于枇杷果实品质与成熟度的高光谱成像检测研究鲜有报道.

    本文以枇杷为研究对象，在可见-近红外(363~1 026 nm)波长区域内获得高光谱信息，建立果实光谱信息与可溶性固形物(soluble solid content，SSC)、硬度、成熟度回归模型，旨在对HSI用于枇杷果实品质无损检测和成熟度的预测潜力进行评估.

1.   材料与方法

    1.1.   样品的制备

  • 以西南大学果树学重点实验室枇杷种质资源圃的红肉枇杷“金华1号”品种为实验材料，采摘大小均匀、无疤痕、无病虫害、成熟度七成(经验判断)以上的果实，采后2 h内运至实验室，放置于7℃的冷藏柜中保存，于3 d内完成高光谱图像获取和理化实验分析，试验样本总计115个.

    • 1.2.   仪器

  • CT-3型质构仪，美国Brookfield公司. PAL-1型数字折射仪，日本ATAGO CO.，LTD.. 可见/近红外高光谱成像系统(363~1 026 nm): 包括1个线性扫描成像仪(Imspector V10E，Spectral Imaging Ltd.)，1个高性能的EMCCD相机(Raptor EM285CL)，1个由步进电机驱动的移动平台(Isuzu Optics Corp.)，2个光纤卤素灯(IT 3900，150W)，1台处理数据的计算机和配套图像获取软件(Isuzu Optics Corp.). 系统放置在封闭的箱体内，如图 1.

    • 1.3.   图像采集及校正

  • 将待测的枇杷样本从冷藏柜中取出，放置于20℃室温环境下1 h后进行高光谱图像采集. 参数设置: 光谱分辨率为0.5 nm，样本与镜头之间距离为20 cm，曝光时间为500 ms，载物台移动速度为0.65 mm/s，图像的分辨率为1 632×1 232像素.

    为消除光照对图像的影响，使用配套图像处理软件对采集的图像进行黑白校正，校正图像(R)的计算公式：

    式中，R0为原始图像，W为白标定图像，B为黑标定图像.

    黑白校正后的枇杷高光谱图像需要去除背景仅保留样本区域，或选取图像中的感兴趣区域(ROI)，因此，在ENVI 5.1软件中进行ROI提取，并将ROI内所有光谱信息的平均值作为对应反射光谱值，计算公式：

    式中，I为ROI区域的平均光谱，m为ROI区域中像素点的个数，Ii为第i个像素的光谱，共1 232个光谱点的光谱数据.

    • 1.4.   理化检测

    1.4.1.   SSC

  • 去掉样品的果皮，将果肉挤压出汁液于PAL-1型数字折射仪上读数，重复取3次枇杷果肉的汁液，3次读数的平均值作为该样品的SSC值.

    • 1.4.2.   硬度

  • 将果实赤道面果皮小心剥除后，放置于质构仪的测试平板上，探头直径5 mm，测试速度1.0 mm/s，两次压缩停顿时间为5 s，深度5.0 mm，触发点负载0.2 N，重复3次，由质地特征曲线得到果肉硬度.

    • 1.5.   数据处理

  • 光谱数据处理使用Matlab 2018a软件.

    • 1.5.1.   异常值剔除

  • 由于环境、仪器、人为操作等因素造成的实验误差，应对光谱异常和成分含量值异常的样本进行剔除. 蒙特卡洛方法(Monte Carlo，MC)具有同时检测光谱异常值和理化性质异常值的优点，将所有样本作为建模集建立模型，根据交叉验证均方根误差(RMSECV)最小原则确定最佳主成分数，通过计算各样本预测误差均值和标准差，结合3σ判据方法，判定具有较高预测误差均值和标准差的样本为异常样本[14].

    • 1.5.2.   样本集划分

  • 联合X-Y距离(samples set partitioning based on joint X-Y distances，SPXY)样本划分方法能够在对样本间距离计算的同时考虑光谱向量和浓度向量[15]，因此本文采用SPXY对枇杷样本的SSC，硬度，成熟度按3∶1进行建模集和预测集划分，选择样本时计算其样本间距离dxy(i，j)：

    式中，dx (i，j)和dy(i，j)分别为仅以光谱和仅以浓度为特征参数统计的各样本间距离，并通过除以各自最大值，从而使两个特征参量在样本选择时权重相同.

    • 1.5.3.   特征波长选取

  • 获得的高光谱图像中共有1 232个波长变量，冗余较多，可从所有波段中选择可分性好的波段子集，降低数据维度. 因此，对样本进行偏最小二乘回归分析，利用竞争性自适应权重采样算法(competitive adaptive reweighted sampling，CARS)、连续投影算法(successive projections algorithm，SPA)进行特征波长的选取.

    CARS算法是一种结合蒙特卡洛采样与偏最小二乘(Partial least squares regression，PLSR)模型回归系数的特征变量选择方法，模仿达尔文理论中的“适者生存”的原则. CARS算法中，每次通过自适应加权采样保留PLSR模型中回归系数绝对值大的波长点，去掉权值较小的波长点，再利用交叉验证优选出RMSECV最小的波长子集[13]. SPA算法是前向特征变量选择方法，利用向量的投影分析，通过将波长投影到其他波长上，比较投影向量大小，以投影向量最大的波长为待选波长，然后基于矫正模型选择最终的特征波长，SPA选择的是含有最少冗余信息及最小共线性的变量组合[16].

    • 1.5.4.   模型的建立与评价

  • 在光谱的定量分析模型中，PLSR是一种处理两数据块之间关系的数学方法和常用的化学计量学算法，能够同时做到回归建模和数据降维，并结合了典型相关分析和多元线性回归分析、主成分分析[17]. PLSR将n个样品m个组分的浓度矩阵Y=(yi，j)n×m和仪器测定的n个样品p个波长点处吸光度矩阵X=(xi，j) n×p分解为特征向量形式Y = UQ + F和X = TP + E，其中U，T分别为浓度特征因子矩阵和吸光度特征因子矩阵，Q，P分别为浓度载荷阵和吸光度载荷阵，F，E分别为浓度残差阵和吸光度残差阵. 而后建立PLS回归模型:

    式中，Ed为随机误差矩阵，B为d维对角回归系数矩阵. 对未知待测样品，如果吸光度为x，则其浓度(y)可以求解为：

    判别偏最小二乘法(Discriminant partial least squares，DPLS)是基于判别分析基础上的偏最小二乘法，它用类别信息矩阵代替了偏最小二乘法回归模型中的浓度矩阵[18]. 本实验采用PLSR和DPLS进行建模.

    SSC和硬度回归模型建立后，将建模集留一法交叉验证相关系数Rc和均方根误差RMSEC作为回归模型性能的主要评价指标，以独立预测集相关系数Rp，预测集均方根误差(RMSEP)作为模型预测性能主要评价指标[19]. 成熟度分类模型以预测集的判别准确率作为评价指标. 建模集的内部留一法交叉验证，使用式(6)和式(7)计算相关系数和均方根误差，$\hat{y}_{c_i} $为用标准方法测定的结果，yci为通过光谱数据测量及数学模型交互验证预测的结果，nc为建模集样品数量; 外部预测集验证使用式(8)和式(9)计算相关系数和均方根误差，$\hat{y}_{p_i} $为用标准方法测定的结果，ypi为通过光谱数据测量及数学模型预测的结果，np为预测集样品数量.

2.   结果与讨论

    2.1.   高光谱图像预处理

  • 采用ENVI 5.1软件对进行黑白校正后的枇杷高光谱图像进行ROI提取及光谱范围选择. 在枇杷高光谱图像中心部位提取像素400×400的正方形区域为光谱感兴趣区域，每个样本的光谱值以ROI中光谱响应平均值来估算. 分析光谱曲线可知，所有枇杷样本的光谱反射率趋势基本保持一致，在420 nm和680 nm附近有较强光谱吸收特征，这主要是由枇杷中类胡萝卜素和叶绿素吸收光谱所致，980 nm附近的吸收与水分有关(图 2).

    • 2.2.   异常样本剔除

  • 采用蒙特卡洛交叉验证法剔除异常值，以全部样本分别建立SSC和硬度的交叉验证模型，计算样本预测误差均值和标准差，散点图见图 3. 结合3σ判据方法对样本中的异常样本进行识别，采用蒙特卡洛交叉验证法建立SSC模型，异常样本有69号，70号，76号，88号，92号，96号; 采用蒙特卡洛交叉验证法建立硬度模型，异常样本有4号，21号，28号，58号，88号，92号. 为尽可能保证样本的多样性，因此仅剔除共有的88号和92号异常样本. 剩余样本进行交叉验证，SSC和硬度交叉检验模型的相关系数Rcv分别从0.892 3，0.851 3提高到0.955 2，0.859 5. 后续用于建模的总样本数为113个.

    • 2.3.   枇杷品质指标理化测定结果及样本划分

  • 经测定，113个“金华1号”枇杷样本的SSC值范围在5.90~13.57 °Brix之间，硬度范围在4.07~7.17 kg/cm2之间，直方图见图 4. SSC值和硬度值分布呈近似正态分布，说明用于实验的样本选择是合适的. 采用SPXY方法将所有样本进行样本集划分，结果见表 1. 可以看出，用于回归模型建立的建模集包含了两个品种枇杷SSC和硬度的最大值与最小值，分布范围较广，样本集划分合理.

    • 2.4.   特征波段选取

  • 将光谱数据进行预处理后，采用CARS算法提取特征波长建立回归模型. 图 5a为CARS算法选择SSC特征波长的过程，运行次数为39时，选取的波长数为105个，占波长总量的8.52%. 同样方法进行硬度特征波长的选择(图 5b)，当运行次数为46时，硬度值选取的波长数为66个，占波长总量的5.36%. 图 6a为SPA算法优选的SSC特征波长分布，选择的波长数范围N=10~84，共选出12个光谱波长，占波长总量的0.97%; 图 6b为SPA算法优选的硬度值特征波长分布，选择的波长数范围N=10~84，共选出10个光谱波长，占波长总量的0.81%.

    • 2.5.   PLSR模型的建立

  • 建立全光谱(FS)和CARS，SPA算法优选特征波长的PLSR模型，对其建模效果进行比较分析，结果见表 2和表 3.

    由表 2可知，SSC的CARS-PLSR模型与FS-PLSR模型和SPA-PLSR模型相比，Rc值和Rp值均为最大，分别为0.981 7和0.918 5; RMSEC，RMSEP均为最小，分别为0.294 2和0.373 8，具有较好的校正性能和预测性能. 由表 3可知，硬度值的CARS-PLSR模型预测效果也是3个模型中最好的，Rc值和Rp值分别为0.970 7和0.742 3; RMSEC和RMSEP分别为0.113 5和0.165 2. SPA-PLSR模型的预测性能最差，这可能是由于全波长有1 232个光谱点，SPA算法以投影向量最大的波长为待选波长，选出的通常为最小共线性的变量组合，选出的特征波长数较少，而枇杷果实的化学成分较为复杂，导致丢失了一些重要的光谱信息.

    • 2.6.   成熟度预测

  • 果实SSC和可溶性糖值变化较大的时间段是在七成熟与八成熟之间. 邓朝军等[3]、许奇志等[4]研究了不同成熟度枇杷的品质特性，发现七成熟至八成熟可溶性糖值迅速上升. 七成熟的枇杷质地较硬，着色淡黄色，SSC和可溶性糖值较低; 八成熟及以上的枇杷则质地较软，着色深黄，香味浓郁，SSC和可溶性糖值显著升高. 采摘时根据样本的外观，对所有样本的成熟度进行标记，并在测定SSC后进行成熟度标记的复核，其中八成熟及以上的枇杷样本有39个，七成熟的枇杷样本有74个. 采用SPXY法按表 4中的数量进行校正集和预测集样本划分; 光谱采用CARS法进行特征波长选择; 类别信息矩阵中，将八成熟及以上的枇杷成熟度属性赋值为1，七成熟的枇杷成熟度属性赋值为2[18].

    以建模集85个样本建立DPLS定性判别模型，独立预测集28个样本进行模型的验证，预测集样本中有6个八成熟及以上的样本，有两个样本的类别预测值为2，即有两个样本被错判为七成熟; 预测集样本中有22个七成熟的样本，有1个样本的类别预测值为1，即有1个七成熟的样本被错判为八成熟及以上. 采用DPLS算法对枇杷果实成熟度进行预测，预测集28个样本中八成熟及以上的枇杷预测准确率为66.67%，七成熟的枇杷预测准确率为95.45%. 所有样本判错数为3，总的预测准确率为89.29%(表 4).

3.   结论

  • 采用高光谱成像技术，开展了枇杷品质检测方法的研究，通过光谱与图像信息相结合的方法实现了枇杷果实SSC和硬度的预测以及成熟度判别，主要结论如下：

    ① 将枇杷果实的光谱信息和理化值经过蒙特卡洛交叉验证法剔除了2个异常值后，所有样本的SSC和硬度交叉检验模型的相关系数Rcv分别从0.892 3，0.851 3提高到0.955 2，0.859 5.

    ② CARS算法进行波段选取可以有效去除全光谱波段中的无关信息和冗余信息，降低模型的复杂度，提高模型的预测能力. 采用SPXY-CARS-PLSR建立回归模型并进行外部预测集数据验证，SSC的建模集相关系数Rc值和预测集相关系数Rp值分别为0.981 7和0.918 5; RMSEC，RMSEP分别为0.294 2和0.373 8; 硬度的建模集相关系数Rc值和预测集相关系数Rp值分别为0.970 7和0.742 3; RMSEC，RMSEP分别为0.113 5和0.165 2，与FS-PLSR模型和SPA-PLSR模型相比，CARS-PLSR模型具有较好的校正性能和预测性能.

    ③ 采用DPLS建立枇杷果实成熟度判别模型并进行外部预测集数据验证，28个预测集样本共有3个判错，总的判别正确率为89.29%. 八成熟及以上的样本数量较少，还需在后续的研究中增加样本.

参考文献 (19)

目录

版权所有：西南大学学报（自然科学版）

地址：重庆市北碚区西南大学 电话：023-68254576, 68252540, 68366165,68366917 电子邮箱：lkxb@swu.edu.cn

本系统由北京仁和汇智信息技术有限公司开发

/

返回文章
返回