CN112986174A - 一种基于近红外光谱的果蔬优选分拣方法、系统及可读存储介质 - Google Patents
一种基于近红外光谱的果蔬优选分拣方法、系统及可读存储介质 Download PDFInfo
- Publication number
- CN112986174A CN112986174A CN202110151497.7A CN202110151497A CN112986174A CN 112986174 A CN112986174 A CN 112986174A CN 202110151497 A CN202110151497 A CN 202110151497A CN 112986174 A CN112986174 A CN 112986174A
- Authority
- CN
- China
- Prior art keywords
- fruit
- vegetable
- vegetables
- fruits
- near infrared
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 235000012055 fruits and vegetables Nutrition 0.000 title claims abstract description 187
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000002329 infrared spectrum Methods 0.000 title claims abstract description 48
- 230000007547 defect Effects 0.000 claims abstract description 75
- 230000003595 spectral effect Effects 0.000 claims abstract description 66
- 238000001228 spectrum Methods 0.000 claims description 44
- 235000013399 edible fruits Nutrition 0.000 claims description 28
- 235000013311 vegetables Nutrition 0.000 claims description 25
- 241000607479 Yersinia pestis Species 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 18
- 241000238631 Hexapoda Species 0.000 claims description 14
- 238000009499 grossing Methods 0.000 claims description 12
- 238000007781 pre-processing Methods 0.000 claims description 11
- 230000001419 dependent effect Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 8
- 238000013144 data compression Methods 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 7
- 201000010099 disease Diseases 0.000 claims description 6
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 5
- 230000003760 hair shine Effects 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 5
- 238000001035 drying Methods 0.000 claims description 3
- 238000010238 partial least squares regression Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 description 16
- 235000007688 Lycopersicon esculentum Nutrition 0.000 description 14
- 240000003768 Solanum lycopersicum Species 0.000 description 14
- 230000009466 transformation Effects 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 239000013598 vector Substances 0.000 description 8
- 238000012795 verification Methods 0.000 description 7
- 241000196324 Embryophyta Species 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 235000021022 fresh fruits Nutrition 0.000 description 4
- 238000010606 normalization Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000000513 principal component analysis Methods 0.000 description 4
- 238000003672 processing method Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010835 comparative analysis Methods 0.000 description 2
- 238000002790 cross-validation Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000556 factor analysis Methods 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000011241 protective layer Substances 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 208000035143 Bacterial infection Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 208000022362 bacterial infectious disease Diseases 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3563—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B07—SEPARATING SOLIDS FROM SOLIDS; SORTING
- B07C—POSTAL SORTING; SORTING INDIVIDUAL ARTICLES, OR BULK MATERIAL FIT TO BE SORTED PIECE-MEAL, e.g. BY PICKING
- B07C5/00—Sorting according to a characteristic or feature of the articles or material being sorted, e.g. by control effected by devices which detect or measure such characteristic or feature; Sorting by manually actuated devices, e.g. switches
- B07C5/34—Sorting according to other particular properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/2431—Multiple classes
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B07—SEPARATING SOLIDS FROM SOLIDS; SORTING
- B07C—POSTAL SORTING; SORTING INDIVIDUAL ARTICLES, OR BULK MATERIAL FIT TO BE SORTED PIECE-MEAL, e.g. BY PICKING
- B07C2501/00—Sorting according to a characteristic or feature of the articles or material to be sorted
- B07C2501/009—Sorting of fruit
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N2021/3595—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/68—Food, e.g. fruit or vegetables
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Immunology (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开的一种基于近红外光谱的果蔬优选分拣方法、系统及可读存储介质，通过光源照射待分拣果蔬，采集待分拣果蔬的光谱信息，并将获取的光谱信息进行预处理，引入果蔬辨别模型；根据果蔬辨别模型得到预测缺陷类别值，将所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率；判断所述偏差率是否大于预设偏差率阈值，根据判断结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，以实现优选分拣。
Description
技术领域
本发明涉及一种果蔬优选分拣方法，更具体的，涉及一种基于近红外光谱的果蔬优选分拣方法、及可读存储介质。
背景技术
如今随着人们生活水平的提高，对高品质消费愈发推崇。价格不再是影响人们购买的主要因素，在购买商品时更多人开始看重商品品质。有数据显示，在我国水果领域约有八九成的水果被用于鲜食，这就要求对这些新鲜水果进行品质区分从而满足人们对高品质新鲜水果的需求，传统的分拣工作由果农依靠人眼识别和手工操作进行，往往会出现分选出的新鲜水果大小不一，美观度不够甚至是有内部缺陷的问题，同时通过果农手工人眼识别和手工操作进行分选判别极易因疲劳造成失误，导致分选准确率极大地降低，且手工操作会给果蔬带来二次损伤和细菌感染；所以消费者在选购这类新鲜果蔬使只能凭借经验、靠运气，常常因为经验不足选到品质不好的果蔬。
为了能够满足消费者能够吃到高品质的果蔬产品，需要开发一款可以对果蔬进行优选分拣的系统，该系统通过获取待分拣果蔬的光谱信息，将所述的光谱信息进行预处理，引入果蔬辨别模型，根据果蔬辨别模型得到的预测缺陷类别值，判别待分拣果蔬中的有物理损伤和有虫害的果蔬，在进行果蔬分拣的过程中，如何实现对果蔬光谱信息进行处理，如何通过建立预测能力好、误差值小的果蔬辨别模型判断果蔬的品质好坏与否都是亟不可待要解决的问题。
发明内容
为了解决上述至少一个技术问题，本发明提出了一种基于近红外光谱的果蔬优选分拣方法、系统和可读存储介质。
本发明第一方面提供了一种基于近红外光谱的果蔬优选分拣方法，包括：
获取待分拣果蔬的光谱信息；
将所述光谱信息进行预处理，引入果蔬辨别模型，得到结果信息；
根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，以实现优选分拣。
本方案中，所述的获取待分拣果蔬的光谱信息，具体为：光源照射待分拣果蔬，果蔬底座中设有导入检测器，导入检测器连接光谱仪采集光谱信息，其中光谱类型为短波近红外光谱，果蔬底座中设有遮光圈，可以有效的减少杂散光干扰。
本方案中，所述的将光谱信息进行预处理，具体为在MATLAB软件中，进行原始光谱数据的平滑、插值、滤波、拟合及提高分辨率运算以及进行平滑去燥、数据压缩及信息提取，导出能量光谱。
本方案中，在MATLAB软件中建立果蔬辨别模型，使用Plsregress函数实现偏最小二乘回归，在偏最小二乘算法模型中，通过迭代算法实现自变量和因变量的分解并回归建立二者关系，其中建立果蔬辨别模型具体为：将数据进行标准化处理，求解出符合要求的主成分，根据求解出的主成分建立光谱数据与缺陷类别值之间的回归，通过迭代一步继续求解主成分直到达到目标值，推导出光谱数据与缺陷类别值的回归表达式，即果蔬辨别模型建立完成。
本方案中，所述的引入果蔬辨别模型，得到结果信息，具体为：果蔬辨别模型通过计算得到预测缺陷类别值，根据预测缺陷类别值判别待分拣果蔬中有物理损伤和有虫害果蔬，其中预测缺陷类别值的计算公式为：
其中，F为果蔬辨别模型预测缺陷类别值，a为模型截距，k为预测集中第k个预测自变量，n为自变量数量，E表示能量光谱谱线强度，η为模型拟合参数。
本方案中，所述的根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，具体为：
根据果蔬辨别模型得到预测缺陷类别值；
将所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率；
判断所述偏差率是否大于预设偏差率阈值；
若大于，则判断为有物理损伤或有病虫害的果蔬。
本方案中，所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率，其中偏差率的计算具体为：
其中Qj和Pj分别表示在果蔬辨别模型判别过程中第j个缺陷类别值的实际值和预测值，z表示样品的个数。
本发明第二方面还提供了一种基于近红外光谱的果蔬优选分拣系统，该系统包括：存储器、处理器，所述存储器中包括基于近红外光谱的果蔬优选分拣方法程序，所述基于近红外光谱的果蔬优选分拣方法程序被所述处理器执行时实现如下步骤：
获取待分拣果蔬的光谱信息；
将所述光谱信息进行预处理，引入果蔬辨别模型，得到结果信息；
根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，以实现优选分拣。
本方案中，所述的获取待分拣果蔬的光谱信息，具体为：光源照射待分拣果蔬，果蔬底座中设有导入检测器，导入检测器连接光谱仪采集光谱信息，其中光谱类型为短波近红外光谱，果蔬底座中设有遮光圈，可以有效的减少杂散光干扰。
本方案中，所述的将光谱信息进行预处理，具体为在MATLAB软件中，进行原始光谱数据的平滑、插值、滤波、拟合及提高分辨率运算以及进行平滑去燥、数据压缩及信息提取，导出能量光谱。
本方案中，在MATLAB软件中建立果蔬辨别模型，使用Plsregress函数实现偏最小二乘回归，在偏最小二乘算法模型中，通过迭代算法实现自变量和因变量的分解并回归建立二者关系，其中建立果蔬辨别模型具体为：将数据进行标准化处理，求解出符合要求的主成分，根据求解出的主成分建立光谱数据与缺陷类别值之间的回归，通过迭代一步继续求解主成分直到达到目标值，推导出光谱数据与缺陷类别值的回归表达式，即果蔬辨别模型建立完成。
本方案中，所述的引入果蔬辨别模型，得到结果信息，具体为：果蔬辨别模型通过计算得到预测缺陷类别值，根据预测缺陷类别值判别待分拣果蔬中有物理损伤和有虫害果蔬，其中预测缺陷类别值的计算公式为：
其中，F为果蔬辨别模型预测缺陷类别值，a为模型截距，k为预测集中第k个预测自变量，n为自变量数量，E表示能量光谱谱线强度，η为模型拟合参数。
本方案中，所述的根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，具体为：
根据果蔬辨别模型得到预测缺陷类别值；
将所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率；
判断所述偏差率是否大于预设偏差率阈值；
若大于，则判断为有物理损伤或有病虫害的果蔬。
本方案中，所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率，其中偏差率的计算具体为：
其中Qj和Pj分别表示在果蔬辨别模型判别过程中第j个缺陷类别值的实际值和预测值，z表示样品的个数。
本发明第三方面还提供一种计算机可读存储介质，所述计算机可读存储介质中包括基于近红外光谱的果蔬优选分拣方法程序，所述基于近红外光谱的果蔬优选分拣方法程序被处理器执行时，实现如上述任一项所述的基于近红外光谱的果蔬优选分拣方法的步骤。
本发明公开的一种基于近红外光谱的果蔬优选分拣方法、系统和可读存储介质，通过获取待分拣果蔬的光谱信息，将所述的光谱信息进行预处理，引入果蔬辨别模型，根据果蔬辨别模型得到的预测缺陷类别值，通过预测缺陷类别值与预设缺陷类别值的比较，判别待分拣果蔬中的有物理损伤和有虫害的果蔬，实现果蔬的优选分拣。
附图说明
图1示出了本发明一种基于近红外光谱的果蔬优选分拣方法的流程图；
图2示出了本发明根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬的方法流程图。
图3示出了本发明一种基于近红外光谱的果蔬优选分拣系统的框图。
具体实施方法
为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
图1示出了本发明一种基于近红外光谱的果蔬优选分拣方法的流程图。
如图1所示，本发明第一方面提供了一种基于近红外光谱的果蔬优选分拣方法，包括：
S102,获取待分拣果蔬的光谱信息；
S104,将所述光谱信息进行预处理，引入果蔬辨别模型，得到结果信息；
S106,根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，以实现优选分拣。
需要说明的是，所述的获取待分拣果蔬的光谱信息，具体为：光源照射待分拣果蔬，果蔬底座中设有导入检测器，导入检测器连接光谱仪采集光谱信息，其中光谱类型为短波近红外光谱，果蔬底座中设有遮光圈，可以有效的减少杂散光干扰。
需要说明的是，所述的将光谱信息进行预处理，具体为在MATLAB软件中，进行原始光谱数据的平滑、插值、滤波、拟合及提高分辨率运算以及进行平滑去燥、数据压缩及信息提取，导出能量光谱，预处理可以减弱甚至消除各种非目标因素对光谱的影响，尽可能的去处去除无关信息变量，从而提高模型的预测能力和精度；光谱处理方法有数据增强算法、平滑、导数、标准正态变量变换、多元散射校正、傅立叶变换、直线差值、矢量归一化、最大最小归一化、小波变换、正交信号校正和净分析信号等方法及这些方法的组合，其中导数方法可以有效消除基线和其他背景的干扰，分辨重叠峰，提高分辨率和灵敏度，但是它同时会引入噪声，降低信噪比，常用的有一阶导数和二阶导数；近红外光谱标准正态变量变换处理办法主要是消除固体颗粒大小，表面散射以及光程变化对漫反射光谱的影响；去趋势算法通常是用于标准正态算法处理后的光谱，用来消除漫反射光谱的基线漂移；小波变换主要作用为光谱去噪平滑、光谱数据压缩和化学信息提取等；正交信号校正可以滤除与浓度无关的光谱信号，以简化模型和提高模型预测能力。
需要说明的是，偏最小二乘法是一种类似于主成分分析的光谱定量分析方法，同时执行光谱数据和浓度数据的分解，并不像主成分分析算法中有单独回归步骤，计算模型中的每一个新荷载向量时，当某个因素所起的作用从原始数据前得分就被替换，缩减了的数据矩阵用来计算下一个载荷向量，过程重复直到所有期望的载荷向量计算完，偏最小二乘法是比较完善的基于因子分析原理的校正方法，可降低噪声对校正模型的影响，比较适用于处理变量多而样本数少的问题。
需要说明的是，在MATLAB软件中建立果蔬辨别模型，使用Plsregress函数实现偏最小二乘回归，在偏最小二乘算法模型中，通过迭代算法实现自变量和因变量的分解并回归建立二者关系，其中建立果蔬辨别模型具体为：将数据进行标准化处理，求解出符合要求的主成分，根据求解出的主成分建立光谱数据与缺陷类别值之间的回归，通过迭代一步继续求解主成分直到达到目标值，推导出光谱数据与缺陷类别值的回归表达式，即果蔬辨别模型建立完成。
需要说明的是，所述的引入果蔬辨别模型，得到结果信息，具体为：果蔬辨别模型通过计算得到预测缺陷类别值，根据预测缺陷类别值判别待分拣果蔬中有物理损伤和有虫害果蔬，其中预测缺陷类别值的计算公式为：
其中，F为果蔬辨别模型预测缺陷类别值，a为模型截距，k为预测集中第k个预测自变量，n为自变量数量，E表示能量光谱谱线强度，η为模型拟合参数。
图2示出了本发明根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬的方法流程图。
根据本发明实施例，所述的根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，还包括：
S202，根据果蔬辨别模型得到预测缺陷类别值；
S204，将所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率；
S206，判断所述偏差率是否大于预设偏差率阈值；
S208，若大于，则判断为有物理损伤或有病虫害的果蔬。
本方案中，所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率，其中偏差率的计算具体为：
其中Qj和Pj分别表示在果蔬辨别模型判别过程中第j个缺陷类别值的实际值和预测值，z表示样品的个数。
需要说明的是，在建立模型之前需要对样品集进行筛选，选取的样品指标范围应恰当合适，同时应该考虑多个品种的果蔬在实验集和验证集的个数和性质是否平均，从而确保模型的准确性以及通用性，有物理损伤或者有虫害果蔬的光谱信息与正常果蔬的光谱信息有不同，但光谱的总体变化趋势基本一致，例如在研究正常番茄与有物理损伤番茄红外漫投射光谱特征对比分析时，由于有物理损伤番茄外部保护层受到破坏，使得透光能力变强，有物理损伤番茄的光谱能量比正常番茄的光谱能量高，由光谱分析可知，两条对比光谱在波长680nm-730nm范围内同时有波谷、波峰，在波长630nm处有物理损伤的番茄比正常番茄多有两个波峰，综上分析，可得出光谱的有效信息集中在689nm-836nm之间，即选用689nm-730nm波长范围进行建模。
需要说明的是，采用近红外光谱建立模型时，样本集按一定比例划分为校正集和验证集，用校正集建立校正模型并用此模型对验证集进行预测，检验模型的性能；将所得果蔬辨别模型根据比较分析各个模型的交互验证均方根、均方根误差以及相关系数等多个参数来进行模型的评估，若模型的交互验证均方根和均方根误差越小，相关系数越趋近于1，则证明模型的预测能力越好。当主乘数增加时，预测集的均方根误差逐渐降低，当主成分增加到一定数值时，均方根误差最小，模型预测能力达到最优效果。光谱变量在果蔬辨别模型中的贡献率越大，对应的回归系数越大，即当正回归系数对象对应的光谱变量越大，果蔬有物理损伤和有虫害的几率越大。
本发明第二方面还提供了一种基于近红外光谱的果蔬优选分拣系统3，该系统包括：存储器31、处理器32，所述存储器中包括基于近红外光谱的果蔬优选分拣方法程序，所述基于近红外光谱的果蔬优选分拣方法程序被所述处理器执行时实现如下步骤：
获取待分拣果蔬的光谱信息；
将所述光谱信息进行预处理，引入果蔬辨别模型，得到结果信息；
根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，以实现优选分拣。
需要说明的是，所述的获取待分拣果蔬的光谱信息，具体为：光源照射待分拣果蔬，果蔬底座中设有导入检测器，导入检测器连接光谱仪采集光谱信息，其中光谱类型为短波近红外光谱，果蔬底座中设有遮光圈，可以有效的减少杂散光干扰。
需要说明的是，所述的将光谱信息进行预处理，具体为在MATLAB软件中，进行原始光谱数据的平滑、插值、滤波、拟合及提高分辨率运算以及进行平滑去燥、数据压缩及信息提取，导出能量光谱，预处理可以减弱甚至消除各种非目标因素对光谱的影响，尽可能的去处去除无关信息变量，从而提高模型的预测能力和精度；光谱处理方法有数据增强算法、平滑、导数、标准正态变量变换、多元散射校正、傅立叶变换、直线差值、矢量归一化、最大最小归一化、小波变换、正交信号校正和净分析信号等方法及这些方法的组合，其中导数方法可以有效消除基线和其他背景的干扰，分辨重叠峰，提高分辨率和灵敏度，但是它同时会引入噪声，降低信噪比，常用的有一阶导数和二阶导数；近红外光谱标准正态变量变换处理办法主要是消除固体颗粒大小，表面散射以及光程变化对漫反射光谱的影响；去趋势算法通常是用于标准正态算法处理后的光谱，用来消除漫反射光谱的基线漂移；小波变换主要作用为光谱去噪平滑、光谱数据压缩和化学信息提取等；正交信号校正可以滤除与浓度无关的光谱信号，以简化模型和提高模型预测能力。
需要说明的是，偏最小二乘法是一种类似于主成分分析的光谱定量分析方法，同时执行光谱数据和浓度数据的分解，并不像主成分分析算法中有单独回归步骤，计算模型中的每一个新荷载向量时，当某个因素所起的作用从原始数据前得分就被替换，缩减了的数据矩阵用来计算下一个载荷向量，过程重复直到所有期望的载荷向量计算完，偏最小二乘法是比较完善的基于因子分析原理的校正方法，可降低噪声对校正模型的影响，比较适用于处理变量多而样本数少的问题。
需要说明的是，在MATLAB软件中建立果蔬辨别模型，使用Plsregress函数实现偏最小二乘回归，在偏最小二乘算法模型中，通过迭代算法实现自变量和因变量的分解并回归建立二者关系，其中建立果蔬辨别模型具体为：将数据进行标准化处理，求解出符合要求的主成分，根据求解出的主成分建立光谱数据与缺陷类别值之间的回归，通过迭代一步继续求解主成分直到达到目标值，推导出光谱数据与缺陷类别值的回归表达式，即果蔬辨别模型建立完成。
需要说明的是，所述的引入果蔬辨别模型，得到结果信息，具体为：果蔬辨别模型通过计算得到预测缺陷类别值，根据预测缺陷类别值判别待分拣果蔬中有物理损伤和有虫害果蔬，其中预测缺陷类别值的计算公式为：
其中，F为果蔬辨别模型预测缺陷类别值，a为模型截距，k为预测集中第k个预测自变量，n为自变量数量，E表示能量光谱谱线强度，η为模型拟合参数。
根据本发明实施例，所述的根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，还包括：
根据果蔬辨别模型得到预测缺陷类别值；
将所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率；
判断所述偏差率是否大于预设偏差率阈值；
若大于，则判断为有物理损伤或有病虫害的果蔬。
本方案中，所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率，其中偏差率的计算具体为：
其中Qj和Pj分别表示在果蔬辨别模型判别过程中第j个缺陷类别值的实际值和预测值，z表示样品的个数。
需要说明的是，在建立模型之前需要对样品集进行筛选，选取的样品指标范围应恰当合适，同时应该考虑多个品种的果蔬在实验集已经验证集的个数和性质是否平均，从而确保模型的准确性以及通用性，有物理损伤或者有虫害果蔬的光谱信息与正常果蔬的光谱信息有不同，但光谱的总体变化趋势基本一致，例如在研究正常番茄与有物理损伤番茄红外漫投射光谱特征对比分析时，由于有物理损伤番茄外部保护层受到破坏，使得透光能力变强，有物理损伤番茄的光谱能量比正常番茄的光谱能量高，由光谱分析可知，两条对比光谱在波长680nm-730nm范围内同时有波谷、波峰，在波长630nm处有物理损伤的番茄比正常番茄多有两个波峰，综上分析，可得出光谱的有效信息集中在689nm-836nm之间，即选用689nm-730nm波长范围进行建模。
需要说明的是，采用近红外光谱建立模型时，样本集按一定比例划分为校正集和验证集，用校正集建立校正模型并用此模型对验证集进行预测，检验模型的性能；将所得果蔬辨别模型根据比较分析各个模型的交互验证均方根、均方根误差以及相关系数等多个参数来进行模型的评估，若模型的交互验证均方根和均方根误差越小，相关系数越趋近于1，则证明模型的预测能力越好；当主乘数增加时，预测集的均方根误差逐渐降低，当主成分增加到一定数值时，均方根误差最小，模型预测能力达到最优效果。光谱变量在果蔬辨别模型中的贡献率越大，对应的回归系数越大，即当正回归系数对象对应的光谱变量越大，果蔬有物理损伤和有虫害的几率越大。
图3示出了本发明一种基于近红外光谱的果蔬优选分拣系统的框图。
本发明第三方面还提供一种计算机可读存储介质，所述计算机可读存储介质中包括基于近红外光谱的果蔬优选分拣方法程序，所述基于近红外光谱的果蔬优选分拣方法程序被处理器执行时，实现如上述任一项所述的基于近红外光谱的果蔬优选分拣方法的步骤。
本发明公开的一种基于近红外光谱的果蔬优选分拣方法、系统和可读存储介质，通过获取待分拣果蔬的光谱信息，将所述的光谱信息进行预处理，引入果蔬辨别模型，根据果蔬辨别模型得到的预测缺陷类别值，通过预测缺陷类别值与预设缺陷类别值的比较，判别待分拣果蔬中的有物理损伤和有虫害的果蔬，实现果蔬的优选分拣。
在本申请所提供的几个实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，如：多个单元或组件可以结合，或可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口，设备或单元的间接耦合或通信连接，可以是电性的、机械的或其它形式的。
上述作为分离部件说明的单元可以是、或也可以不是物理上分开的，作为单元显示的部件可以是、或也可以不是物理单元；既可以位于一个地方，也可以分布到多个网络单元上；可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。
另外，在本发明各实施例中的各功能单元可以全部集成在一个处理单元中，也可以是各单元分别单独作为一个单元，也可以两个或两个以上单元集成在一个单元中；上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：移动存储设备、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
或者，本发明上述集成的单元如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括：移动存储设备、ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。
Claims (10)
1.一种基于近红外光谱的果蔬优选分拣方法，其特征在于，包括：
获取待分拣果蔬的光谱信息；
将所述光谱信息进行预处理，引入果蔬辨别模型，得到结果信息；
根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，以实现优选分拣。
2.根据权利要求1所述的一种基于近红外光谱的果蔬优选分拣方法，其特征在于，所述的获取待分拣果蔬的光谱信息，具体为：光源照射待分拣果蔬，果蔬底座中设有导入检测器，导入检测器连接光谱仪采集光谱信息，其中光谱类型为短波近红外光谱，果蔬底座中设有遮光圈，可以有效的减少杂散光干扰。
3.根据权利要求1所述的一种基于近红外光谱的果蔬优选分拣方法，其特征在于，所述的将光谱信息进行预处理，具体为在MATLAB软件中，进行原始光谱数据的平滑、插值、滤波、拟合及提高分辨率运算以及进行平滑去燥、数据压缩及信息提取，导出能量光谱。
4.根据权利要求1所述的一种基于近红外光谱的果蔬优选分拣方法，其特征在于，还包括，在MATLAB软件中建立果蔬辨别模型，使用Plsregress函数实现偏最小二乘回归，在偏最小二乘算法模型中，通过迭代算法实现自变量和因变量的分解并回归建立二者关系，其中建立果蔬辨别模型具体为：将数据进行标准化处理，求解出符合要求的主成分，根据求解出的主成分建立光谱数据与缺陷类别值之间的回归，通过迭代一步继续求解主成分直到达到目标值，推导出光谱数据与缺陷类别值的回归表达式，即果蔬辨别模型建立完成。
5.根据权利要求1所述的一种基于近红外光谱的果蔬优选分拣方法，其特征在于，所述的引入果蔬辨别模型，得到结果信息，具体为：果蔬辨别模型通过计算得到预测缺陷类别值，根据预测缺陷类别值判别待分拣果蔬中有物理损伤和有虫害果蔬，其中预测缺陷类别值的计算公式为：
其中，F为果蔬辨别模型预测缺陷类别值，a为模型截距，k为预测集中第k个预测自变量，n为自变量数量，E表示能量光谱谱线强度，η为模型拟合参数。
6.根据权利要求1所述的一种基于近红外光谱的果蔬优选分拣方法，其特征在于，所述的根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，具体为：
根据果蔬辨别模型得到预测缺陷类别值；
将所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率；
判断所述偏差率是否大于预设偏差率阈值；
若大于，则判断为有物理损伤或有病虫害的果蔬。
7.一种基于近红外光谱的果蔬优选分拣系统，其特征在于，该系统包括：存储器、处理器，所述存储器中包括基近红外光谱的果蔬优选分拣方法程序，所述基于近红外光谱的果蔬优选分拣方法程序被所述处理器执行时实现如下步骤：
获取待分拣果蔬的光谱信息；
将所述光谱信息进行预处理，引入果蔬辨别模型，得到结果信息；
根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，以实现优选分拣。
8.根据权利要求7所述的一种基于近红外光谱的果蔬优选分拣系统，其特征在于，所述的引入果蔬辨别模型，得到结果信息，具体为：果蔬辨别模型通过计算得到预测缺陷类别值，根据预测缺陷类别值判别待分拣果蔬中有物理损伤和有虫害果蔬，其中预测缺陷类别值的计算公式为：
其中，F为果蔬辨别模型预测缺陷类别值，a为模型截距，k为预测集中第k个预测自变量，n为自变量数量，E表示能量光谱谱线强度，η为模型拟合参数。
9.根据权利要求7所述的一种基于近红外光谱的果蔬优选分拣系统，其特征在于，所述的根据结果信息判别待分拣果蔬中有物理损伤和有虫害的果蔬，具体为：
根据果蔬辨别模型得到预测缺陷类别值；
将所述的预测缺陷类别值与预设缺陷类别值进行比较，得到偏差率；
判断所述偏差率是否大于预设偏差率阈值；
若大于，则判断为有物理损伤或有病虫害的果蔬。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中包括基于近红外光谱的果蔬优选分拣方法程序，所述基于近红外光谱的果蔬优选分拣方法程序被处理器执行时，实现如权利要求1至6中任一项所述的基于近红外光谱的果蔬优选分拣方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110151497.7A CN112986174A (zh) | 2021-02-03 | 2021-02-03 | 一种基于近红外光谱的果蔬优选分拣方法、系统及可读存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110151497.7A CN112986174A (zh) | 2021-02-03 | 2021-02-03 | 一种基于近红外光谱的果蔬优选分拣方法、系统及可读存储介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112986174A true CN112986174A (zh) | 2021-06-18 |
Family
ID=76346559
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110151497.7A Pending CN112986174A (zh) | 2021-02-03 | 2021-02-03 | 一种基于近红外光谱的果蔬优选分拣方法、系统及可读存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112986174A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114871120A (zh) * | 2022-05-26 | 2022-08-09 | 江苏省徐州医药高等职业学校 | 基于图像数据处理的药品确定分拣方法及装置 |
| CN117074353A (zh) * | 2023-08-18 | 2023-11-17 | 广东省农业科学院设施农业研究所 | 一种荔枝果实蒂蛀虫无损检测方法和系统 |
| CN117368146A (zh) * | 2023-12-08 | 2024-01-09 | 苏陀科技(北京)有限公司 | 一种菌丝体蛋白含量的快速检测方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5808305A (en) * | 1996-10-23 | 1998-09-15 | Src Vision, Inc. | Method and apparatus for sorting fruit in the production of prunes |
| US5884775A (en) * | 1996-06-14 | 1999-03-23 | Src Vision, Inc. | System and method of inspecting peel-bearing potato pieces for defects |
| CN1430723A (zh) * | 2000-03-13 | 2003-07-16 | 奥特莱有限公司 | 用可见光/近红外光谱测量和相关水果特性的方法和装置 |
| CN106323909A (zh) * | 2016-09-14 | 2017-01-11 | 江苏大学 | 果蔬品质手持式近红外光谱检测系统及检测方法 |
| US20190340749A1 (en) * | 2017-01-26 | 2019-11-07 | ClariFruit | System and method for evaluating fruits and vegetables |
-
2021
- 2021-02-03 CN CN202110151497.7A patent/CN112986174A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5884775A (en) * | 1996-06-14 | 1999-03-23 | Src Vision, Inc. | System and method of inspecting peel-bearing potato pieces for defects |
| US5808305A (en) * | 1996-10-23 | 1998-09-15 | Src Vision, Inc. | Method and apparatus for sorting fruit in the production of prunes |
| CN1430723A (zh) * | 2000-03-13 | 2003-07-16 | 奥特莱有限公司 | 用可见光/近红外光谱测量和相关水果特性的方法和装置 |
| CN106323909A (zh) * | 2016-09-14 | 2017-01-11 | 江苏大学 | 果蔬品质手持式近红外光谱检测系统及检测方法 |
| US20190340749A1 (en) * | 2017-01-26 | 2019-11-07 | ClariFruit | System and method for evaluating fruits and vegetables |
Non-Patent Citations (7)
| Title |
|---|
| 刘德华等: "基于高光谱成像技术的山楂损伤和虫害缺陷识别研究", 《光谱学与光谱分析》 * |
| 刘燕德等: "鸭梨黑心病可见/近红外漫透射光谱在线检测", 《光谱学与光谱分析》 * |
| 刘燕德等: "黄桃表面缺陷和可溶性固形物光谱同时在线检测", 《农业工程学报》 * |
| 司守奎 等: "《数学建模算法与应用》", 31 May 2015, 国防工业出版社 * |
| 向延菊等: "近红外光谱在果蔬品质定性分析中的应用研究进展", 《食品工业》 * |
| 孙海燕 等: "《数理统计》", 31 October 2016, 北京航空航天大学出版社 * |
| 李龙等: "苹果内外品质在线无损检测分级系统设计与试验", 《农业工程学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114871120A (zh) * | 2022-05-26 | 2022-08-09 | 江苏省徐州医药高等职业学校 | 基于图像数据处理的药品确定分拣方法及装置 |
| CN114871120B (zh) * | 2022-05-26 | 2023-11-07 | 江苏省徐州医药高等职业学校 | 基于图像数据处理的药品确定分拣方法及装置 |
| CN117074353A (zh) * | 2023-08-18 | 2023-11-17 | 广东省农业科学院设施农业研究所 | 一种荔枝果实蒂蛀虫无损检测方法和系统 |
| CN117074353B (zh) * | 2023-08-18 | 2024-05-03 | 广东省农业科学院设施农业研究所 | 一种荔枝果实蒂蛀虫无损检测方法和系统 |
| CN117368146A (zh) * | 2023-12-08 | 2024-01-09 | 苏陀科技(北京)有限公司 | 一种菌丝体蛋白含量的快速检测方法 |
| CN117368146B (zh) * | 2023-12-08 | 2024-03-12 | 苏陀科技(北京)有限公司 | 一种菌丝体蛋白含量的快速检测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112986174A (zh) | 一种基于近红外光谱的果蔬优选分拣方法、系统及可读存储介质 | |
| Gomes et al. | Comparison of different approaches for the prediction of sugar content in new vintages of whole Port wine grape berries using hyperspectral imaging | |
| Ataş et al. | A new approach to aflatoxin detection in chili pepper by machine vision | |
| Suárez et al. | Artificial intelligence in agriculture | |
| Lorente et al. | Visible–NIR reflectance spectroscopy and manifold learning methods applied to the detection of fungal infections on citrus fruit | |
| Wang et al. | Technology using near infrared spectroscopic and multivariate analysis to determine the soluble solids content of citrus fruit | |
| US11682203B2 (en) | Feature extraction method, model training method, detection method of fruit spectrum | |
| Shafiee et al. | Combined data mining/NIR spectroscopy for purity assessment of lime juice | |
| Li et al. | Supervised classification of slightly bruised peaches with respect to the time after bruising by using hyperspectral imaging technology | |
| Zhang et al. | A simple identification model for subtle bruises on the fresh jujube based on NIR spectroscopy | |
| CN106442561A (zh) | 用于检测柑橘表皮缺陷的在线图像采集系统及方法 | |
| Wang et al. | Identification of tomatoes with early decay using visible and near infrared hyperspectral imaging and image‐spectrum merging technique | |
| Zhang et al. | Identification of common skin defects and classification of early decayed citrus using hyperspectral imaging technique | |
| Xu et al. | Research on nondestructive identification of grape varieties based on EEMD‐DWT and hyperspectral image | |
| Nogales-Bueno et al. | Comparative study on the use of anthocyanin profile, color image analysis and near-infrared hyperspectral imaging as tools to discriminate between four autochthonous red grape cultivars from La Rioja (Spain) | |
| Chen et al. | Discrimination of Fresh Tobacco Leaves with Different Maturity Levels by Near‐Infrared (NIR) Spectroscopy and Deep Learning | |
| Font et al. | An image processing method for in-line nectarine variety verification based on the comparison of skin feature histogram vectors | |
| Ropelewska | The application of image processing for cultivar discrimination of apples based on texture features of the skin, longitudinal section and cross-section | |
| CN110208211B (zh) | 一种用于农药残留检测的近红外光谱降噪方法 | |
| CN113751345B (zh) | 一种菠萝无损检测分级装置和方法 | |
| CN110308111A (zh) | 一种应用近红外光谱技术快速预测远安黄茶闷黄时间的方法 | |
| CN110533102A (zh) | 基于模糊推理的单类分类方法以及分类器 | |
| CN106841070A (zh) | 一种白酒真伪鉴定方法及装置 | |
| CN112782148B (zh) | 阿拉比卡和罗伯斯塔两种咖啡豆的快速鉴别方法 | |
| CN114965486A (zh) | 基于高光谱成像对柿子的轻微瘀伤检测方法与装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210618 |
|
| RJ01 | Rejection of invention patent application after publication |