光纤光谱技术对猕猴桃品质及成熟度的无损检测

光纤光谱技术对猕猴桃品质及成熟度的无损检测 Nondestructive detection for kiwifruit quality and maturity by optical fiber spectroscopy technology 1 2 first-author 尚　静（1988-），女，河北张家口人，硕士，讲师，2017年于四川大学获得硕士学位，现为贵阳学院农产品无损检测工程研究中心主要成员，主要从事基于光谱技术的农产品无损检测方面的研究。E-mail： shji0124@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848661&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848663&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848662&type= 22.01333618 32.00401306 尚　静 （1988-），女，河北张家口人，硕士，讲师，2017年于四川大学获得硕士学位，现为贵阳学院农产品无损检测工程研究中心主要成员，主要从事基于光谱技术的农产品无损检测方面的研究。E-mail： shji0124@163.com shji0124@163.com 1 2 1 2 2 corresp E-mail： eileen_zyool@sohu.com E-mail： eileen_zyool@sohu.com 张　艳（1977-），女，贵州遵义人，博士后，教授，分别于2004年、2007年在四川大学获得硕士、博士学位，现为贵阳学院农产品无损检测工程研究中心负责人，主要从事农产品无损检测方面的研究。E-mail： eileen_zy001@sohu.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848664&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848667&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848665&type= 22.01333618 32.00402832 张　艳 （1977-），女，贵州遵义人，博士后，教授，分别于2004年、2007年在四川大学获得硕士、博士学位，现为贵阳学院农产品无损检测工程研究中心负责人，主要从事农产品无损检测方面的研究。E-mail： eileen_zy001@sohu.com eileen_zy001@sohu.com 贵阳学院 食品与制药工程学院，贵州 贵阳 550005 贵阳学院 食品与制药工程学院，贵州 贵阳 550005 Food and Pharmaceutical Engineering Institute， Guiyang University， Guiyang 550005， China Food and Pharmaceutical Engineering Institute， Guiyang University， Guiyang 550005， China 贵阳学院 农产品无损检测工程研究中心，贵州 贵阳 550005 贵阳学院 农产品无损检测工程研究中心，贵州 贵阳 550005 Research Center of Nondestructive Testing for Agricultural Products， Guiyang University， Guiyang 550005， China Research Center of Nondestructive Testing for Agricultural Products， Guiyang University， Guiyang 550005， China 猕猴桃可溶性固形物含量（SSC）和硬度是评价其品质的关键参数，同时也是判别其成熟度的重要指标。为探究基于光纤光谱技术预测猕猴桃SSC、硬度和成熟度的可行性并寻求最佳预测模型。首先，采用光纤光谱（200~1 000 nm）采集系统获取不同成熟期“贵长”猕猴桃的反射光谱，并测定SSC和硬度的参考值。接着，基于全光谱和参考值构建偏最小二乘回归（PLSR）和主成分回归（PCR）预测模型。然后，应用连续投影算法（SPA）和竞争性自适应重加权算法（CARS）选取特征波长，构建简化的多元线性回归（MLR）和误差反向传播（BP）网络预测模型。最后，通过偏最小二乘判别分析（PLS-DA）和简化的K近邻（SKNN）算法，构建预测猕猴桃成熟度检测模型。结果表明：CARS-BP模型对SSC的预测性能最优，其预测集决定系数 <math id="M1"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848674&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848673&type= 3.72533321 3.38666677 =0.90，预测集均方根误差（RMSEP）和剩余预测偏差（RPD）分别为0.64和3.22；CARS-MLR对硬度的预测性能相对最优，其 <math id="M2"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848674&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848673&type= 3.72533321 3.38666677 =0.83，RMSEP和RPD分别为1.67和2.47；PLS-DA模型对猕猴桃成熟度的检测性能最优，其正确识别率高达100%。该研究为水果品质和成熟度的无损检测提供重要指导。 The soluble solids content （SSC） and firmness of kiwifruit are two important indices for evaluating its quality and distinguishing its maturity. In this paper， we explore the feasibility of predicting the SSC， firmness， and maturity of kiwifruit using optical fiber spectroscopy technology and of finding the best prediction model. First， an optical fiber spectroscopy （200~1 000 nm） acquisition system was used to collect the reflectance spectra of the different maturity stages of ‘Guichang’ kiwifruit. Simultaneously， the reference values of the SSC and firmness were measured. Two methods， namely， partial least-squares regression （PLSR） and principal components regression （PCR）， were employed to establish the models on the basis of the full spectra and the reference values. Then， multiple linear regression （MLR） and an error back-propagation （BP） network were applied to build simplified models on the basis of the selected characteristic variables from the full wavelengths using the methods of successive projection algorithm （SPA） and competitive adaptive reweighted sampling （CARS）. Finally， partial least-squares discrimination analysis （PLS-DA） and a simplified K nearest neighbor （SKNN） algorithm were applied to build models for predicting the maturity of kiwifruit. The results showed that， for the SSC， the CARS-BP model had the best prediction ability （ <math id="M3"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848678&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848677&type= 3.72533321 3.38666677 =0.90， RMSEP=0.64， RPD=3.22）， and for the firmness， the CARS-MLR model had the best prediction ability （ <math id="M4"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848678&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848677&type= 3.72533321 3.38666677 =0.83， RMSEP=1.67， RPD=2.47）. The PLS-DA model had the best detection ability， and the maturity discrimination accuracy was up to 100%. These results can provide important guidance for the nondestructive prediction of the quality and maturity of fruits. 光纤光谱 猕猴桃 可溶性固形物含量 硬度 成熟度 无损检测 optical fiber spectroscopy kiwifruit soluble solids content firmness maturity nondestructive detection 1 引　言 1 “贵长”猕猴桃拥有果肉细嫩，肉质多浆，果汁丰富，清甜爽口，酸甜适中的独特品质，作为后熟型水果之一，实际采收时为了延长其贮藏期，经常采摘还未成熟的果实，但是如果过早采摘，果实就会特别生硬，这将会使口感受到影响，果实也易受冷害；若过晚采摘，果实十分柔软就会因为易腐烂而难以贮藏 ［ 1 1 ］ 。猕猴桃的可溶性固形物含量（Soluble Solids Content，SSC）和硬度是衡量其品质和成熟度的关键指标，因此，开展猕猴桃SSC、硬度和成熟度的无损检测对于指导其采收时间、采后储藏和加工具有重要的意义。 水果SSC和硬度的传统检测方法是采用生化检测手段，然而这种方法最大的缺点是破损样本 ［ 2 2 ］ 。基于光谱技术的无损检测方法具有诸多优势，比如分析速度快、无污染、无损伤等，所以，近年来，这种方法已被广泛地用于水果品质的快速检测，引起国内外研究学者的颇多关注 ［ 3 3 - 12 12 3-12 ］ 。在预测SSC方面，董金磊等人 ［ 13 13 ］ 基于特征光谱构建的连续投影算法-误差反向传播（Successive Projection Algorithm-error Back Propagation，SPA-BP）网络模型对猕猴桃SSC的预测性能相对最优，其预测集相关系数（ R P ）和预测集均方根误差（Root Mean Squares Errors for Prediction，RMSEP）分别为0.92和0.77。詹白勺等人 ［ 14 14 ］ 基于竞争性自适应重加权算法（Competitive Adaptive Reweighted Sampling，CARS）提取的特征光谱构建的非线性最小二乘支持向量机（LS-SVM）模型对库尔勒香梨SSC的预测效果相对最优，其预测集决定系数（ <math id="M5"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 ），RMSEP和剩余预测偏差（Residual Predictive Deviation，RPD）分别为0.85，0.29，2.59。在预测硬度方面，郭文川等 ［ 15 15 ］ 通过高光谱成像技术结合人工神经网络构建的全光谱BP网络模型拥有相对最优的预测效果，其 R P 和RMSEP分别为0.86和0.93。冯迪等 ［ 16 16 ］ 利用高光谱成像技术提取了可以同时预测苹果SSC和硬度的最佳波长。综上所述，利用光纤光谱技术检测猕猴桃SSC、硬度和成熟度是可行的。 本文采用光纤光谱采集系统获取不同成熟期“贵长”猕猴桃的反射光谱，为确定预测猕猴桃SSC、硬度和成熟度最佳模型，先基于全光谱和参考值构建偏最小二乘回归和主成分回归预测模型；然后分别应用竞争性自适应重加权算法和连续投影算法对全波段光谱数据进行降维，并基于提取的特征光谱构建简化的多元线性回归和误差反向传播网络预测模型；最后通过偏最小二乘判别分析和简化的 K 近邻算法，构建猕猴桃成熟度预测模型，选取最优模型，以期为水果品质和成熟度无损检测装备的开发提供理论基础。 2 试验材料、仪器与方法 1 2.1　材料 2 不同成熟期（未熟期、半熟期、成熟期、过熟期）“贵长”猕猴桃于2019年9月17日~2019年10月19日分4批次采自贵州省修文县龙关口猕猴桃果园，样品从不同的果树上随机采摘，每批次采摘50个无病虫害且无机械损伤的样品，共计200个。样品采摘后立即运到实验室，在实验室温度为（22±2）℃下进行实验，用柔软的纸品轻轻地擦掉猕猴桃表面的尘土等杂物，依次对其编号后采集光谱并测定SSC和硬度。 2.2　仪器 2 实验中使用的仪器主要有光纤光谱采集系统（蔚海光学仪器（上海）有限公司）、ATAGO PAL- α 型数字手持袖珍折射仪（日本ATAGO公司）、GY-4型数显果实硬度计（杭州绿博仪器有限公司）、TD4Z-WS型台式低速离心机（湖南湘立科学仪器有限公司）以及JYZ-V911榨汁机（九阳股份有限公司）。 如 图1 图1 所示，光纤光谱采集系统主要由：光谱仪（QEPro，波长范围为198.2~1 006.4 nm，分辨率为2.84~3 nm）、光纤（R600-7-VIS-125F，直径为600 μm）、卤钨灯光源（HL-2000，波长范围为360~2 400 nm）、反射探头支架（RPH-1）、适配器（RPH-ADP）、标准反射白板（WS-1）和高性能计算机（Lenovo）等构成。 2.3　光谱数据采集 2 实验根据测试标准反射白板的光谱强度，调节积分时间和平均次数等参数，并保证测量猕猴桃样品时的参数一致。其中，系统的积分时间为110 ms，扫描平均次数为8，滑动平均宽度为1。将待测猕猴桃紧贴在RPH-1表面上（反射探头通过RPH-ADP固定在RPH-1上，RPH-1表面距离猕猴桃约1 cm），光谱采集部位为猕猴桃赤道位置，5次采集的平均值作为该样品的光谱数据。 2.4　SSC和硬度的测定 2 实验使用数显果实硬度计测定猕猴桃样品的硬度，硬度计的探针直径长7.9 mm，探头压入果子的深度约10 mm，计量单位：kg/cm 2 。测量时先将猕猴桃测量部位去皮（测量部位为猕猴桃赤道位置），然后将探头正对果肉处，缓慢匀速压入，至刻线处止，每个猕猴桃分别测定3个位置，3个位置的平均值为该样品的参考值。 实验使用数字手持袖珍折射仪测定猕猴桃的SSC，折射仪的测量范围：0.0~85 °Brix，测量精度：±0.2 °Brix。先将猕猴桃榨汁后再离心（离心机的转速为3 000 r/min，离心时间为5 min），再将猕猴桃汁涂抹在折光棱镜的镜面上，连续按测量键按钮多次，当最后液晶显示屏3次显示值一致时记录该值，作为该样品的参考值。 2.5　建模及模型评价方法 2 实验基于全波段光谱数据和参考值分别构建偏最小二乘回归（Partial Least Square Regression，PLSR）和主成分回归（Principal Components Regression，PCR）预测模型，并采用连续投影算法（Successive Projection Algorithm，SPA）和竞争性自适应重加权算法（Competitive Adaptive Reweighted Sampling，CARS）选取特征波长，进而基于提取的特征光谱和参考值分别构建简化的多元线性回归（Multiple Linear Regression，MLR）和误差反向传播（Error Back Propagation，BP）网络预测模型，并通过偏最小二乘判别分析（Partial Least Square Discrimination Analysis，PLS-DA）和简化的K近邻（Simplified K Nearest Neighbor，SKNN）算法，构建猕猴桃成熟度预测模型，最终确定最优预测模型。 将校正集的决定系数（ <math id="M6"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 ）及其均方根误差（Root Mean Squares Errors for Calibration，RMSEC），预测集的决定系数（ <math id="M7"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 ）及其均方根误差（Root Mean Squares Errors for Prediction，RMSEP），剩余预测偏差（Residual Predictive Deviation，RPD）作为模型性能的评价指标。通常一个好的预测模型拥有较大的 <math id="M8"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 ， <math id="M9"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 和RPD值，以及较小的RMSEC和RMSEP值。RPD值位于1.5~2之间意味着模型的预测性能较差，RPD值位于2~2.5之间表明模型可以粗略地完成定量预测，当RPD值位于2.5~3之间甚至大于3时说明模型具有非常好的预测性能 ［ 17 17 ］ 。 <math display="block" id="M10"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>C</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>=</mo><mn mathvariant="normal">1</mn><mo>-</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>C</mi></mrow></msub></mrow></munderover></mstyle><mrow><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">l</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mrow><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>C</mi></mrow></msub></mrow></munderover></mstyle><mrow><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848696&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848695&type= 50.54600143 19.64266586 ， <math display="block" id="M11"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>=</mo><mn mathvariant="normal">1</mn><mo>-</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>P</mi></mrow></msub></mrow></munderover></mstyle><mrow><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">p</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mrow><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>P</mi></mrow></msub></mrow></munderover></mstyle><mrow><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848698&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848697&type= 50.46133423 19.64266586 ， <math display="block" id="M12"><mi mathvariant="normal">R</mi><mi mathvariant="normal">M</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">E</mi><mi mathvariant="normal">C</mi><mo>=</mo><mroot><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>C</mi></mrow></msub></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>C</mi></mrow></msub></mrow></munderover></mstyle><mrow><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">l</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848701&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848699&type= 58.16600037 10.58333302 ， <math display="block" id="M13"><mi mathvariant="normal">R</mi><mi mathvariant="normal">M</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">E</mi><mi mathvariant="normal">P</mi><mo>=</mo><mroot><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>P</mi></mrow></msub></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mi>P</mi></mrow></msub></mrow></munderover></mstyle><mrow><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">p</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi></mrow></msub><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848704&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848703&type= 58.08132935 10.58333302 ， <math display="block" id="M14"><mi mathvariant="normal">R</mi><mi mathvariant="normal">P</mi><mi mathvariant="normal">D</mi><mo>=</mo><mfrac><mrow><mi mathvariant="normal">S</mi><mi mathvariant="normal">D</mi></mrow><mrow><mi mathvariant="normal">R</mi><mi mathvariant="normal">M</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">E</mi><mi mathvariant="normal">P</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848706&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848705&type= 26.67000008 7.61999989 ， 其中： n C 和 n P 为校正集和预测集中的猕猴桃数， y act 和 y mean 为样品参考值及其平均值， y cal 和 y pre 为校正集和预测集中SSC和硬度预测值，SD表示预测集中参考值的标准偏差。 2.6　数据分析 2 实验采用OceanView（Ocean Optics，Copyright 2013）软件采集光谱，利用MATLAB R2016b软件对光谱处理、划分样本集、建模。 3 结果与讨论 1 3.1　SSC和硬度测定结果 2 表1 表1 给出了不同成熟期猕猴桃SSC和硬度的统计数据， 表1 表1 可表明，猕猴桃SSC值随着其逐渐成熟而增大，硬度值随着其逐渐成熟而减小。在构建回归预测模型之前，需要基于采集的光谱和测定的参考值将所有猕猴桃样本划分为校正样本集和预测样本集，校正集中参考值的范围要比预测集的范围宽是划分样本集的评判标准。实验采用样本集划分方法（Sample set Partitioning based on joint x - y distances，SPXY） ［ 18 18 ］ 将所有猕猴桃按照3：1的比例划分为校正样本集和预测样本集。校正样本集和预测样本集中猕猴桃SSC和硬度统计数据见 表2 表2 ，从 表2 表2 可见，校正集中SSC和硬度的范围要比预测集中的范围宽，这样划分的样本集有助于后期构建较优的预测模型。 参数 成熟度 范围 平均值 标准差 4 SSC/（°Brix） 未熟期 4.80~5.57 5.23 0.16 半熟期 6.00~7.20 6.57 0.28 成熟期 7.10~9.00 8.11 0.50 过熟期 9.50~12.20 10.81 0.63 4 硬度/ （kg·cm -2 ） 未熟期 20.39~26.19 22.51 1.55 半熟期 15.91~18.15 17.10 0.51 成熟期 13.80~16.88 15.39 0.62 过熟期 9.10~15.68 12.36 1.86 参数 样本集 范围 平均值 标准差 2 SSC/（°Brix） 校正集（150） 4.80~12.20 8.28 2.03 预测集（50） 4.93~11.40 8.29 2.07 2 硬度/ （kg·cm -2 ） 校正集（150） 9.10~26.19 17.44 4.26 预测集（50） 10.13~25.36 18.15 4.13 3.2　反射光谱及预处理 2 由于原始光谱的首末两头含有较多的噪声，因此剔除前10个和后10个波段，选择波段206.33~999.06 nm为有效光谱区域，该区域共有1 024个波段。为了提升预测模型的精确度和稳定性，进一步采用标准正态变换（Standard Normal Variation，SNV）对原始光谱数据进行预处理 ［ 19 19 ］ 。 图2 图2 分别给出了猕猴桃样本的原始反射光谱（a）以及经过SNV（b）预处理后的相对反射光谱。 1,2; 3.3　基于全光谱建立的预测模型 2 实验首先基于全光谱分别构建了PLSR和PCR预测模型，建模结果见 表3 表3 。从 表3 表3 可见，构建的PCR预测模型对猕猴桃SSC和硬度检测性能均优于PLSR预测模型，其中PCR预测模型对SSC的检测性能相对较优，其 <math id="M15"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 和RMSEP分别为0.88和0.70，RPD高达2.94，表明该预测模型具有相对较好的预测性能。 2 指标 2 模型 2 潜变量数 2 校正集 3 预测集 决定系数 均方根误差 决定系数 均方根误差 剩余偏差 2 SSC PLSR 8 0.88 0.70 0.88 0.71 2.92 PCR 16 0.83 0.82 0.88 0.70 2.94 2 硬度 PLSR 10 0.91 1.26 0.78 1.89 2.18 PCR 9 0.77 2.05 0.83 1.66 2.49 3.4　特征光谱的提取 2 由于采集的反射光谱包含较多的波段（1 024个），且波段之间存在较大的相关性及冗余信息，为提高预测模型的检测效率和稳定性，对全波段光谱进行特征波长的提取。 3 . 4 . 1　SPA提取特征波长 3 利用SPA提取特征波长时，通常依据均方根误差（RMSE）的最小值来选定最优的特征波长数。 在预测SSC时，RMSE随SPA中有效波长数的变化情况如 图3 图3 （a）所示，表明RMSE随有效波长数的增加而减小，当有效波长数大于5时，RMSE减小的趋势不明显，因此将这5个波长作为预测SSC的特征波长， 图3 图3 （b）给出了提取的特征波长分布图，选取的特征波长按重要性排列分别为312.69，230.83，397.63，255.24和886.38 nm。 在预测硬度时， 图4 图4 （a）给出了RMSE随SPA中有效波长数的变化情况，当有效波长数大于6时，RMSE减小的趋势不明显，进而将这6个波长作为预测硬度的特征波长，特征波长分布如 图4 图4 （b）所示，选取的特征波长按重要性排列分别为337.63，398.42，233.27，251.99，235.72和885.63 nm。 3 . 4 . 2　CARS提取特征波长 3 应用CARS提取特征波长时，其蒙特卡洛采样次数设定为50次，利用5折交叉验证方法计算所构建的偏最小二乘（PLS）模型的交叉验证均方根误差（RMSECV）。 图5 图5 给出了应用CARS提取预测SSC特征波长结果， 图5 图5 （a）表明第29次采样获得的RMSECV值最小，该最优波长子集包含29个特征波长，提取的特征波长分布情况见 图5 图5 （b）。 图6 图6 给出了应用CARS提取预测硬度特征波长的结果，表明第27次采样获得的RMSECV值最小，该最优波长子集包含37个特征波长， 图6 图6 （b）给出了提取的特征波长分布情况。 3.5　基于特征光谱建立的预测模型 2 分别将SPA和CARS提取的特征光谱作为MLR和BP网络模型的自变量，校正集和预测集中猕猴桃SSC和硬度参考值作为因变量，构建预测不同成熟期猕猴桃SSC和硬度检测模型。MLR和BP模型对猕猴桃SSC和硬度的预测结果见 表4 表4 。 2 指标 2 模型 2 特征波长提取方法 2 特征 波长数 2 校正集 3 预测集 决定 系数 均方根 误差 决定 系数 均方根 误差 剩余 偏差 4 SSC 2 MLR SPA 5 0.72 1.07 0.89 0.68 3.03 CARS 29 0.91 0.61 0.88 0.70 2.94 2 BP SPA 5 0.75 1.01 0.84 0.81 2.54 CARS 29 0.90 0.66 0.90 0.64 3.22 4 硬度 2 MLR SPA 6 0.62 2.62 0.84 1.64 2.52 CARS 37 0.91 1.27 0.83 1.67 2.47 2 BP SPA 6 0.66 2.49 0.84 1.62 2.55 CARS 37 0.90 1.33 0.82 1.75 2.36 从 表4 表4 可以明显看出，对于预测SSC而言，CARS-BP预测模型具有相对较好的校正性能和预测性能，其 <math id="M16"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 和 <math id="M17"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 均为0.90，RMSEC和RMSEP分别为0.66和0.64，其中RPD=3.22大于3进一步说明该预测模型对不同成熟度猕猴桃SSC的预测效果非常好。 对于预测硬度而言，虽然CARS-MLR预测模型的预测性能稍微劣于SPA-MLR和SPA-BP预测模型，但其校正性能明显好于这两种模型，因此，CARS-MLR预测模型的性能相对较好，其 <math id="M18"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 和 <math id="M19"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 分别为0.91和0.83，RMSEC和RMSEP分别为1.27和1.67，另外剩余预测偏差RPD=2.47接近于2.5说明CARS-MLR预测模型可以用于预测猕猴桃硬度。 对比 表3 表3 和 表4 表4 可以看出，CARS-BP预测模型对猕猴桃SSC的预测性能相对最优，而CARS-MLR预测模型对猕猴桃硬度的预测性能相对最优。另外，采用CARS从全光谱1 024个波段中分别提取了29个和37个特征波长作为预测猕猴桃SSC和硬度的特征光谱，大大提升了预测模型的检测效率。其中， 图7 图7 给出了CARS-BP预测模型对SSC的预测结果，CARS-MLR预测模型对硬度的预测结果如 图8 图8 所示。 3.6　成熟度预测 2 将不同成熟度猕猴桃样本的光谱数据作为自变量，样本的类别信息（未熟期、半熟期、成熟期、过熟期）作为因变量，构建预测猕猴桃成熟度检测模型。 表5 表5 给出了SKNN和PLS-DA模型对猕猴桃成熟度的预测结果。从表中可见PLS-DA模型的性能要优于SKNN模型，其正确识别率达到了100%。 2 模型 2 潜变 量数 2 校正集（150） 2 预测集（50） 误判数 识别率/% 误判数 识别率/% SKNN 8 1 99.3 0 100 PLS-DA 6 0 100 0 100 4 结　论 1 本研究以不同成熟期“贵长”猕猴桃为研究对象，探究基于光纤光谱技术预测猕猴桃SSC、硬度和成熟度的可行性并寻求最佳预测模型。主要结论如下： （1）基于采集的全光谱数据，以及测定的SSC和硬度参考值，分别构建基于全波段光谱的PLSR和PCR预测模型，得出PCR模型对SSC和硬度检测性能均优于PLSR模型，其中PCR模型对SSC的预测性能相对较好， <math id="M20"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 =0.88，RMSEP为0.70，RPD为2.94。 （2）为实现光谱数据的降维，应用SPA和CARS选取特征波长，得出相比于SPA特征波长提取方法，CARS更适合光谱数据的降维，采用CARS从1 024个全波段光谱中分别提取了29和37个特征波长用于预测SSC和硬度，明显地提升了预测模型的检测效率；而且构建的CARS-BP模型对SSC的预测性能相对最优，其 <math id="M21"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 =0.90，RMSEP为0.64，RPD高达3.22，而CARS-MLR对硬度的预测性能相对最优，其 <math id="M22"><msubsup><mrow><mi>R</mi></mrow><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848785&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848782&type= 4.14866638 3.97933316 =0.83，RMSEP和RPD分别为1.67和2.47。 （3）通过PLS-DA和SKNN判别分析方法，构建猕猴桃成熟度预测模型，得出PLS-DA模型对猕猴桃成熟度的检测性能相对最优，其正确识别率高达100%。 综上所述，采用光纤光谱技术对猕猴桃SSC、硬度和成熟度的无损检测是可行的，为水果品质和成熟度无损检测装备的开发提供理论基础。 参考文献 张承 ， 李 明 ， 龙友华 ， 等 . 采前喷施壳聚糖复合膜对猕猴桃软腐病的防控及其保鲜作用 ［J］. 食品科学 ， 2016 ， 37 （ 22 ）： 274 - 281 . 采前喷施壳聚糖复合膜对猕猴桃软腐病的防控及其保鲜作用 ZHANG CH ， LI M ， LONG Y H ， et al . Control of soft rot of kiwifruit by pre-harvest application of chitosan composite film and tis effect on quality-improving and fresh-keeping of fruit ［J］. Food Science ， 2016 ， 37 （ 22 ）： 274 - 281 . （in Chinese） Control of soft rot of kiwifruit by pre-harvest application of chitosan composite film and tis effect on quality-improving and fresh-keeping of fruit 程丽娟 ， 刘贵珊 ， 万国玲 ， 等 . 可见/近红外高光谱成像技术对长枣中葡萄糖含量的无损检测 ［J］. 发光学报 ， 2019 ， 40 （ 8 ）： 1055 - 1063 . 可见/近红外高光谱成像技术对长枣中葡萄糖含量的无损检测 CHENG L J ， LIU G S ， WAN G L ， et al . Non-destructive detective of glucose contect in lingwu jujube by Vis /NIR hyperspectral imaging technology ［J］. Chinese Journal of Luminescence ， 2019 ， 40 （ 8 ）： 1055 - 1063 . （in Chinese） Non-destructive detective of glucose contect in lingwu jujube by Vis /NIR hyperspectral imaging technology 郭文川 ， 王铭海 ， 谷静思 ， 等 . 近红外光谱结合极限学习机识别贮藏期的损伤猕猴桃 ［J］. 光学 精密工程 ， 2013 ， 21 （ 10 ）： 2720 - 2727 . 近红外光谱结合极限学习机识别贮藏期的损伤猕猴桃 GUO W CH ， WANG M H ， GU J S ， et al . Nondestructive detection on firmness of peaches based on hyperspectral imaging and artificial neural networks ［J］. Opt. Precision Eng. ， 2013 ， 21 （ 10 ）： 2720 - 2727 . （in Chinese） Nondestructive detection on firmness of peaches based on hyperspectral imaging and artificial neural networks 邵园园 ， 王永贤 ， 玄冠涛 ， 等 . 高光谱成像的猕猴桃货架期快速预测 ［J］. 光谱学与光谱分析 ， 2020 ， 40 （ 6 ）： 1940 - 1946 . 高光谱成像的猕猴桃货架期快速预测 SHAO Y Y ， WANG Y X ， XUAN G T ， et al . Hyperspectral imaging technique for estimating the shelf-life of kiwifruits ［J］. Spectroscopy and Spectral Analysis ， 2020 ， 40 （ 6 ）： 1940 - 1946 . （in Chinese） Hyperspectral imaging technique for estimating the shelf-life of kiwifruits 刘燕德 ， 马奎荣 ， 孙旭东 ， 等 . 梨和苹果糖度在线检测通用数学模型研究 ［J］. 光谱学与光谱分析 ， 2017 ， 37 （ 7 ）： 2177 - 2183 . 梨和苹果糖度在线检测通用数学模型研究 LIU Y D ， MA K R ， SUN X D ， et al . The fruits soluble solids content detection online using universal mathematical model ［J］. Spectroscopy and Spectral Analysis ， 2017 ， 37 （ 7 ）： 2177 - 2183 . （in Chinese） The fruits soluble solids content detection online using universal mathematical model HU W ， SUN D ， 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 . Rapid monitoring 1-MCP-induced modulation of sugars accumulation in ripening ‘hayward’ kiwifruit by Vis/NIR hyperspectral imaging MA T ， LI X ， INAGAKI T ， et al . Noncontact evaluation of soluble solids content in apples by near-infrared hyperspectral imaging ［J］. Journal of Food Engineering ， 2018 ， 224 ： 53 - 61 . Noncontact evaluation of soluble solids content in apples by near-infrared hyperspectral imaging ZHANG D ， XU L ， WANG Q ， et al . The optimal local model selection for robust and fast evaluation of soluble solid content in melon with thick peel and large size by VIS-NIR spectroscopy ［J］. Food Analytical Methods ， 2018 ， 12 ： 136 - 147 . The optimal local model selection for robust and fast evaluation of soluble solid content in melon with thick peel and large size by VIS-NIR spectroscopy XIE C ， CHU B ， HE Y . Prediction of banana color and firmness using a novel wavelengths selection method of hyperspectral imaging ［J］. Food Chemistry ， 2018 ， 245 ： 132 - 140 . Prediction of banana color and firmness using a novel wavelengths selection method of hyperspectral imaging YU X ， LU H ， WU D . Development of deep learning method for predicting firmness and soluble solid content of postharvest korla fragrant pear using VIS/NIR hyperspectral reflectance imaging ［J］. Postharvest Biology and Technology ， 2018 ， 141 ： 39 - 49 . Development of deep learning method for predicting firmness and soluble solid content of postharvest korla fragrant pear using VIS/NIR hyperspectral reflectance imaging ZHANG D ， ZHAN B ， 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 ： 1 - 9 . Determination of soluble solids content in oranges using visible and near infrared full transmittance hyperspectral imaging with comparative analysis of models ZHAO Y ， ZHANG C ， ZHU S ， et al . Shape induced reflectance correction for non-destructive determination and visualization of soluble solids content in winter jujubes using hyperspectral imaging in two different spectral ranges ［J］. Postharvest Biology and Technology ， 2020 ， 161 ： 1 - 11 . Shape induced reflectance correction for non-destructive determination and visualization of soluble solids content in winter jujubes using hyperspectral imaging in two different spectral ranges 董金磊 ， 郭文川 . 采后猕猴桃可溶性固形物含量的高光谱无损检测 ［J］. 食品科学 ， 2015 ， 36 （ 16 ）： 101 - 106 . 采后猕猴桃可溶性固形物含量的高光谱无损检测 DONG J L ， GUO W CH . Nondestructive detection of soluble solid content of postharvest kiwifruits based on hyperspectral imaging technology ［J］. Food Science ， 2015 ， 36 （ 16 ）： 101 - 106 . （in Chinese） Nondestructive detection of soluble solid content of postharvest kiwifruits based on hyperspectral imaging technology 詹白勺 ， 倪君辉 ， 李 军 . 高光谱技术结合CARS算法的库尔勒香梨可溶性固形物定量测定 ［J］. 光谱学与光谱分析 ， 2014 ， 34 （ 10 ）： 2752 - 2757 . 高光谱技术结合CARS算法的库尔勒香梨可溶性固形物定量测定 ZHAN B SH ， NI J H ， LI J . Hyperspectral technology combined with CARS algorithm to quantitatively determine the SSC in korla fragrant pear ［J］. Spectroscopy and Spectral Analysis ， 2014 ， 34 （ 10 ）： 2752 - 2757 . （in Chinese） Hyperspectral technology combined with CARS algorithm to quantitatively determine the SSC in korla fragrant pear 郭文川 ， 董金磊 . 高光谱成像结合人工神经网络无损检测桃的硬度 ［J］. 光学 精密工程 ， 2015 ， 23 （ 6 ）： 1530 - 1537 . 高光谱成像结合人工神经网络无损检测桃的硬度 GUO W CH ， DONG J L . Nondestructive detection on firmness of peaches based on hyperspectral imaging and artificial neural networks ［J］. Opt. Precision Eng. ， 2015 ， 23 （ 6 ）： 1530 - 1537 . （in Chinese） . Nondestructive detection on firmness of peaches based on hyperspectral imaging and artificial neural networks 冯 迪 ， 纪建伟 ， 张 莉 ， 等 . 基于高光谱成像提取苹果糖度与硬度最佳波长 ［J］. 发光学报 ， 2017 ， 38 （ 6 ）： 799 - 806 . 基于高光谱成像提取苹果糖度与硬度最佳波长 FENG D ， JI J W ， ZHANG L ， et al . Optimal wavelengths extraction of apple brix and firmness based on hyperspectral imaging ［J］. Chinese Journal of Luminescence ， 2017 ， 38 （ 6 ）： 799 - 806 . （in Chinese） Optimal wavelengths extraction of apple brix and firmness based on hyperspectral imaging 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 ， 46 （ 2 ）： 99 - 118 . Nondestructive measurement of fruit and vegetable quality by means of NIR spectroscopy： a review 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 . A method for calibration and validation subset partitioning 倪力军 ， 张立国 . 基础化学计量学及其应用 ［M］. 上海 ： 华东理工大学出版社 ， 2011 . 基础化学计量学及其应用 NI L J ， ZHANG L G . Basic Chemometrics and Its Application ［M］. Shanghai ： East China University of Science and Technology Press ， 2011 . （in Chinese） Basic Chemometrics and Its Application 10.37188/OPE.20212905.1190 O433.4 TP391.4 A 1004-924X（2021）05-1190-09 尚静,孟庆龙,黄人帅等.光纤光谱技术对猕猴桃品质及成熟度的无损检测[J].光学精密工程,2021,29(05):1190-1198. 尚静,孟庆龙,黄人帅等.光纤光谱技术对猕猴桃品质及成熟度的无损检测[J].光学精密工程,2021,29(05):1190-1198. SHANG Jing,MENG Qing-long,HUANG Ren-shuai,et al.Nondestructive detection for kiwifruit quality and maturity by optical fiber spectroscopy technology[J].Optics and Precision Engineering,2021,29(05):1190-1198. SHANG Jing,MENG Qing-long,HUANG Ren-shuai,et al.Nondestructive detection for kiwifruit quality and maturity by optical fiber spectroscopy technology[J].Optics and Precision Engineering,2021,29(05):1190-1198. 国家自然科学基金资助项目（No.61505036）；贵州省基础研究计划资助项目（No.黔科合基础［2020］1Y270）；贵阳市财政支持贵阳学院学科建设与研究生教育项目资助项目（No.SY-2020）；大学生创新创业项目资助（No.PX-62058） 2020-08-13 2020-10-08 2021-05-15 2021-06-24 光学精密工程 5 29