CN108760655A - A kind of apple sense of taste profile information method for visualizing - Google Patents

A kind of apple sense of taste profile information method for visualizing Download PDF

Info

Publication number
CN108760655A
CN108760655A CN201810435827.3A CN201810435827A CN108760655A CN 108760655 A CN108760655 A CN 108760655A CN 201810435827 A CN201810435827 A CN 201810435827A CN 108760655 A CN108760655 A CN 108760655A
Authority
CN
China
Prior art keywords
taste
apple
spectral
image
sample
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.)
Granted
Application number
CN201810435827.3A
Other languages
Chinese (zh)
Other versions
CN108760655B (en
Inventor
刘晶晶
刘付龙
王晴晴
韩晓菊
门洪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeast Electric Power University
Original Assignee
Northeast Dianli University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeast Dianli University filed Critical Northeast Dianli University
Priority to CN201810435827.3A priority Critical patent/CN108760655B/en
Publication of CN108760655A publication Critical patent/CN108760655A/en
Application granted granted Critical
Publication of CN108760655B publication Critical patent/CN108760655B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00Features of devices classified in G01N21/00
    • G01N2201/12Circuits of general importance; Signal processing
    • G01N2201/129Using chemometrical methods

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种苹果味觉图谱信息可视化方法,包括如下步骤:S1、基于电子舌苹果样本味觉数据的获取;S2、基于高光谱技术苹果样本数据的获取与预处理;S3、特征波段的选取;S4、苹果味觉可视化模型的构建;S5、在GS‑SVM模型建立之后,将10个波段下的预测光谱值按点输入,求出各点光谱值下的味觉预测输出值记为K(i,j),将味觉值转化为对应的实际色彩值,将苹果味觉信息可视化呈现。本发明利用苹果样本在高光谱单点差异上的优点与电子舌味觉整体信息检测技术相结合,从而实现了苹果味觉图谱信息的可视化。

The invention discloses a method for visualizing apple taste map information, comprising the following steps: S1, acquisition of taste data of apple samples based on electronic tongue; S2, acquisition and preprocessing of apple sample data based on hyperspectral technology; S3, selection of characteristic bands ; S4, the construction of the apple taste visualization model; S5, after the GS-SVM model is established, the predicted spectral values under the 10 bands are input by points, and the taste predicted output values under the spectral values of each point are calculated as K(i , j), the taste value is converted into the corresponding actual color value, and the apple taste information is visualized. The invention combines the advantages of apple samples in hyperspectral single-point difference with the electronic tongue taste overall information detection technology, thereby realizing the visualization of apple taste map information.

Description

一种苹果味觉图谱信息可视化方法A method for visualization of apple taste map information

技术领域technical field

本发明涉及苹果味觉分析领域,具体涉及一种苹果味觉图谱信息可视化方法。The invention relates to the field of apple taste analysis, in particular to a method for visualizing apple taste map information.

背景技术Background technique

中国是苹果生产大国,其产量占苹果总产量的65%。苹果果实富含矿物质和维生素,可溶性大,易被人体吸收,故有“活水”之称,是人们经常食用的水果之一。苹果的味觉信息作为反映其品质的重要因素之一,它影响着大多数消费者是否购买时的选择,准确、高效的检测并表征味觉信息对苹果的育种、种植、储存等具有现实指导意义。传统的理化检测方法无法反映苹果的味觉感官信息,而由于鉴评员心理及周围环境等因素影响导致最常用的人工感官鉴评结果不够客观。基于此,SA-402B型电子舌作为单味觉智能仿生检测系统以其客观、精准等优点逐渐取代传统检测方法在味觉信息方面的应用。对苹果样本进行抽样检测、批次处理,实现样本整体味觉信息的检测,并不能反映出各种味觉信息在样本空间上的分布情况。China is a big apple producing country, and its output accounts for 65% of the total apple output. Apple fruit is rich in minerals and vitamins. It is highly soluble and easily absorbed by the human body. Therefore, it is called "living water" and is one of the fruits that people often eat. As one of the important factors reflecting its quality, the taste information of apples affects the choice of most consumers when purchasing. Accurate and efficient detection and characterization of taste information has practical guiding significance for apple breeding, planting and storage. Traditional physical and chemical testing methods cannot reflect the taste and sensory information of apples, and the most commonly used artificial sensory evaluation results are not objective enough due to factors such as the assessor's psychology and the surrounding environment. Based on this, the SA-402B electronic tongue, as a single-taste intelligent bionic detection system, gradually replaces the application of traditional detection methods in taste information with its advantages of objectivity and precision. Sampling and batch processing of apple samples to detect the overall taste information of the sample cannot reflect the distribution of various taste information in the sample space.

发明内容Contents of the invention

为解决上述问题,本发明提供了一种苹果味觉图谱信息可视化方法。In order to solve the above problems, the present invention provides a method for visualizing apple taste map information.

为实现上述目的,本发明采取的技术方案为:In order to achieve the above object, the technical scheme that the present invention takes is:

一种苹果味觉图谱信息可视化方法,包括如下步骤:A method for visualizing apple taste map information, comprising the steps of:

S1、基于电子舌苹果样本味觉数据的获取S1. Acquisition of taste data of apple samples based on electronic tongue

S11、选取同一批次、同一地点采摘的果形正、大小均匀、无缺陷或污染物的苹果,按国标GB/T10651-2008选出90个优等果品;S11. Select apples picked in the same batch and at the same location with positive shape, uniform size, no defects or pollutants, and select 90 high-quality fruits according to the national standard GB/T10651-2008;

S12、将样品编号后置于常温20±2℃,相对湿度55±5%贮藏24小时,保持温湿度不变;S12. Store the sample at room temperature 20±2°C and relative humidity 55±5% for 24 hours after numbering, keeping the temperature and humidity unchanged;

S13、将每个苹果样本清洗、削皮、榨汁、过滤取上层清液40mL分别置于2个纯净测量杯中,按编号顺序依次经电子舌的检测得到30秒处90×3维糖度、酸味、咸味三种味觉数据;测试开始前,传感器先在正负极清洗溶液中清洗90s,结束后在参比溶液中清洗120s,在另一种参比溶液中继续清洗120s,传感器平衡归零30s;达到平衡后,开始进行检测,测试时间为30s,每次测量结束后自动进入清洗步骤;S13. Wash, peel, squeeze, and filter each apple sample to take 40 mL of the supernatant and place them in two pure measuring cups, and obtain 90 × 3-dimensional sugar content in 30 seconds through the detection of the electronic tongue in sequence according to the numbering order. Sour and salty taste data; before the test, the sensor was cleaned in the positive and negative cleaning solutions for 90s, and then in the reference solution for 120s, and then in another reference solution for 120s, and the sensor balance returned to zero. Zero 30s; after the balance is reached, the test starts, the test time is 30s, and the cleaning step is automatically entered after each measurement;

S2、基于高光谱技术苹果样本数据的获取与预处理S2. Acquisition and preprocessing of apple sample data based on hyperspectral technology

S21、通过下式采用黑白标定方法对高光谱图像进行校正,以消除噪声的影响:S21. Correct the hyperspectral image by using the black-and-white calibration method through the following formula to eliminate the influence of noise:

式中,Rd-暗图像,Rw-白板的漫反射图像,Rs-苹果样本原始的漫反射光谱图像,R-校正后的漫反射光谱图。In the formula, R d - dark image, R w - diffuse reflectance image of whiteboard, R s - original diffuse reflectance spectrum image of apple sample, R - corrected diffuse reflectance spectrum image.

S22、对黑白标定后的图像进行矢量化处理,获得图像描述曲线,圈定目标图像后,生成特定感兴趣掩膜图像,控制样本影像处理区域;S22. Carry out vectorization processing on the black-and-white calibrated image, obtain an image description curve, delineate the target image, generate a specific mask image of interest, and control the sample image processing area;

S23、对掩膜处理后的图像进行影像裁剪,掩膜后的光谱区域除感兴趣处外其他光谱数值均为0。S23. Perform image cropping on the masked image, and the spectral values of the masked spectral region are all 0 except for the region of interest.

S24、获取预处理后的90个苹果光谱图像,其中,由于苹果为类球形水果,故在用高光谱分选仪进行照射样品时将一个苹果样本划分为4个面,并在每个面上分别取编号1-5的5个感兴趣区域,每个感兴趣区域大小为300像素点,求取每个样品5个表面整体的光谱均值,最终得到90×256维的苹果数据;S24. Obtain 90 apple spectral images after preprocessing, wherein, since apples are spherical fruits, an apple sample is divided into 4 surfaces when a hyperspectral sorter is used to irradiate the sample, and each surface Take 5 regions of interest numbered 1-5 respectively, each region of interest is 300 pixels in size, calculate the average value of the spectra of the 5 surfaces of each sample, and finally obtain 90×256 dimensional apple data;

S3、特征波段的选取S3. Selection of characteristic bands

分别选取甜、酸、咸3种基础味觉相对应的敏感波段,采用变异系数确定最佳敏感波段数量,然后采用灰色关联度(GRA)的方法确定苹果样品的味觉信息与光谱图像之间的关联程度;The sensitive bands corresponding to the three basic tastes of sweet, sour and salty were respectively selected, the optimal number of sensitive bands was determined by the coefficient of variation, and then the correlation between the taste information of the apple sample and the spectral image was determined by the gray relational degree (GRA) method degree;

S4、苹果味觉可视化模型的构建S4. Construction of Apple Taste Visualization Model

采用SVM预测各个味觉信息离散点的浓度值,其中超平面函数、RBF核函数及回归函数公式(2)、(3)、(4)分别如下:Using SVM to predict the concentration value of each discrete point of taste information, the hyperplane function, RBF kernel function and regression function formulas (2), (3) and (4) are as follows:

K(x,x′)=exp(g Px-x′P2) (2)K(x, x')=exp(g Px-x'P 2 ) (2)

f(x)=wφ(x)+b (3)f(x)=wφ(x)+b (3)

式中,w-超平面的法向向量,φ(x)-非线性映射函数,b-偏置量,g-宽度系数;In the formula, w-the normal vector of the hyperplane, φ(x)-non-linear mapping function, b-bias, g-width coefficient;

S5、味觉可视化呈现S5. Taste visual presentation

在GS-SVM模型建立之后,将10个波段下的预测光谱值按点输入,求出各点光谱值下的味觉预测输出值记为K(i,j),将味觉值转化为对应的实际色彩值,将苹果味觉信息可视化呈现。After the GS-SVM model is established, the predicted spectral values under the 10 bands are input point by point, and the taste prediction output value under the spectral values of each point is calculated as K(i, j), and the taste value is converted into the corresponding actual value. Color value, which visualizes the taste information of apples.

优选地,所述灰色关联度法具体包括如下步骤:首先,将多量纲数据转化成统一的无量纲形式,定义参考数列为个味觉信息数列,比较数列为380-1038nm各波段下的光谱信息数列来进行无量纲处理,以此消除味觉与光谱信息量纲的差异;然后,求取各单一味觉信息数列与全波段下的光谱信息数列的灰色关联系数;最后,求取各波段的灰色关联度。Preferably, the gray relational degree method specifically includes the following steps: first, convert the multidimensional data into a unified dimensionless form, define a reference sequence as a taste information sequence, and compare the sequence as a spectrum information sequence under each band of 380-1038nm to perform dimensionless processing, so as to eliminate the difference in dimension between taste and spectral information; then, calculate the gray correlation coefficient between each single taste information sequence and the spectral information sequence under the full band; finally, calculate the gray correlation degree of each band .

优选地,所述步骤S4基于遗传算法(GA)和网格式搜索(GS)的方法来对参数c和g进行寻优,在遗传算法中,最大的遗传代数为100,初始种群数为20,参数c的搜索范围是0到100,g为0到100;在网格式搜索方法中,以0.5为间隔进行参数寻优,参数c和g的搜索范围是2-10到210Preferably, said step S4 optimizes the parameters c and g based on genetic algorithm (GA) and grid search (GS). In genetic algorithm, the maximum genetic algebra is 100, and the initial population number is 20. The search range of the parameter c is from 0 to 100, and the value of g is from 0 to 100; in the grid search method, the parameters are optimized at an interval of 0.5, and the search range of the parameters c and g is from 2 -10 to 2 10 .

优选地,所述步骤S4中在建立模型的过程中,在进行光谱数据和味觉信息映射时,将每个波段下的光谱值进行均值求取后再进行模型的建立;其中利用Kennard-Stone方法将2/3个光谱数据样本选为训练集,1/3个光谱数据选为预测集,基于GS-SVM和GA-SVM实现光谱-味觉信息的可视化分析。Preferably, in the process of establishing the model in the step S4, when performing spectral data and taste information mapping, the spectral values under each band are averaged before establishing the model; wherein the Kennard-Stone method is used Select 2/3 spectral data samples as the training set and 1/3 spectral data as the prediction set, and realize the visual analysis of spectral-taste information based on GS-SVM and GA-SVM.

本发明利用苹果样本高光谱单点差异上的优点与电子舌味觉整体信息检测技术相结合,从而实现了苹果味觉图谱信息的可视化,从而为苹果味觉的分析提供了一种较为准确的分析方法。The present invention combines the advantages of apple sample hyperspectral single-point difference with the electronic tongue taste overall information detection technology, thereby realizing the visualization of apple taste map information, thereby providing a relatively accurate analysis method for apple taste analysis.

附图说明Description of drawings

图1为本发明实施例中的光谱图像描述曲线。Fig. 1 is a spectrum image description curve in an embodiment of the present invention.

图2为本发明实施例中的矢量掩膜图像。Fig. 2 is a vector mask image in an embodiment of the present invention.

图3为本发明实施例中的感兴趣区域光谱曲线。Fig. 3 is a spectral curve of a region of interest in an embodiment of the present invention.

图4为本发明实施例中的变异系数变化曲线图Fig. 4 is the variation curve diagram of coefficient of variation in the embodiment of the present invention

图5为本发明实施例中的参数寻优过程示意图;FIG. 5 is a schematic diagram of a parameter optimization process in an embodiment of the present invention;

其中,(a)为遗传算法参数寻优过程;(b)为网格式搜索参数寻优过程。Among them, (a) is the genetic algorithm parameter optimization process; (b) is the grid search parameter optimization process.

图6为本发明实施例中的基于GS-SVM的味觉可视化结果图。Fig. 6 is a graph showing the visualization results of taste sensation based on GS-SVM in the embodiment of the present invention.

具体实施方式Detailed ways

为了使本发明的目的及优点更加清楚明白,以下结合实施例对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。In order to make the objects and advantages of the present invention clearer, the present invention will be further described in detail below in conjunction with the examples. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.

实施例Example

本实施例选取同一批次、同一地点采摘的果形正、大小均匀、无缺陷或污染物的阿克苏糖心苹果为研究对象,按国标GB/T10651-2008选出优等果品90个;包括如下步骤:In this example, Aksu candied apples picked in the same batch and at the same location are selected as the research object, with fruit shape, uniform size, and no defects or pollutants, and 90 high-quality fruits are selected according to the national standard GB/T10651-2008; including the following steps :

S1、基于电子舌苹果样本味觉数据的获取S1. Acquisition of taste data of apple samples based on electronic tongue

将样品编号后置于常温20℃(±2℃),相对湿度55%(±5%)贮藏24小时,保持温湿度基本不变。按编号顺序依次测量90个样本,将每个苹果样本清洗、削皮、榨汁、过滤取上层清液40mL分别置于2个纯净测量杯中待测。测试开始前,传感器先在正负极清洗溶液中清洗90s,结束后在参比溶液中清洗120s,在另一种参比溶液中继续清洗120s,传感器平衡归零30s。达到平衡后,开始进行检测,测试时间为30s,每次测量结束后自动进入清洗步骤。经电子舌的检测得到30秒处90×3维三种味觉(糖度、酸味、咸味)数据。After numbering the sample, store it at room temperature at 20°C (±2°C) and relative humidity at 55% (±5%) for 24 hours, keeping the temperature and humidity basically unchanged. Measure 90 samples sequentially according to the numbering order, wash, peel, squeeze, and filter each apple sample to get 40mL of supernatant and place them in 2 pure measuring cups for testing. Before the test starts, the sensor is cleaned in the positive and negative electrode cleaning solution for 90s, and then in the reference solution for 120s, then in another reference solution for 120s, and the sensor balance is reset to zero for 30s. After the balance is reached, the detection is started, the test time is 30s, and the cleaning step is automatically entered after each measurement. The 90×3-dimensional three-dimensional taste (sugar, sour, salty) data at 30 seconds is obtained through the detection of the electronic tongue.

S2、基于高光谱技术苹果样本数据的获取与预处理S2. Acquisition and preprocessing of apple sample data based on hyperspectral technology

高光谱图像采集过程中,因光源在各波段下强度分布的差异性和相机暗电流噪声的影响,会夹杂部分噪声信息。这些噪声信息会影响高光谱图像的质量,进而影响高光谱图像定性或定量分析模型的精度和稳定性。因此采用黑白标定方法对高光谱图像进行校正,以消除噪声的影响,公式(1)所示。During the hyperspectral image acquisition process, due to the difference in the intensity distribution of the light source in each band and the influence of the dark current noise of the camera, some noise information will be included. These noise information will affect the quality of hyperspectral images, and then affect the accuracy and stability of qualitative or quantitative analysis models of hyperspectral images. Therefore, the black and white calibration method is used to correct the hyperspectral image to eliminate the influence of noise, as shown in formula (1).

式中,Rd-暗图像,Rw-白板的漫反射图像,Rs-苹果样本原始的漫反射光谱图像,R-校正后的漫反射光谱图。In the formula, R d - dark image, R w - diffuse reflectance image of whiteboard, R s - original diffuse reflectance spectrum image of apple sample, R - corrected diffuse reflectance spectrum image.

对黑白标定后的图像进行矢量化处理,获得如图1所示图像描述曲线,圈定目标图像后,生成如图2所示特定感兴趣掩膜图像,控制样本影像处理区域。由于获取的图像边缘包含较多的光谱噪声会增加后期数据处理难度,故对掩膜处理后的图像进行影像裁剪,掩膜后的光谱区域除感兴趣处外其他光谱数值均为0。Vectorize the black and white calibrated image to obtain the image description curve as shown in Figure 1. After delineating the target image, generate a specific mask image of interest as shown in Figure 2 to control the sample image processing area. Since the edge of the acquired image contains a lot of spectral noise, which will increase the difficulty of data processing in the later stage, the image after mask processing is cropped, and the spectral values of the masked spectral region are all 0 except for the region of interest.

获取预处理后的90个苹果光谱图像。其中,由于苹果为类球形水果,故在用高光谱分选仪进行照射样品时将一个苹果样本划分为4个面,并在每个面上分别取编号1-5的5个感兴趣区域(每个感兴趣区域大小大约为300像素点)。每个截面上5个感兴趣区域的平均光谱值如图3所示,求取每个样品5个表面整体的光谱均值,最终得到90×256维的苹果数据。S3、特征波段的选取Get 90 preprocessed apple spectral images. Among them, since apples are spherical fruits, an apple sample is divided into 4 surfaces when a hyperspectral sorter is used to irradiate the sample, and 5 regions of interest (numbered 1-5) are taken on each surface ( Each region of interest is approximately 300 pixels in size). The average spectral values of the 5 regions of interest on each section are shown in Figure 3, and the average spectral values of the 5 surfaces of each sample are calculated to obtain the 90×256-dimensional apple data. S3. Selection of characteristic bands

光谱图像反应苹果内部品质的差异,电子舌通过苹果样品中游离出的官能团得到味觉信息。因此采用灰色关联度(GRA)的方法确定两者之间的关联程度,分别选取甜、酸、咸3种基础味觉相对应的敏感波段。在进行灰色关联分析之前,采用变异系数确定最佳敏感波段数量。综合变异系数越大表明变量间的相关性越低,研究以2为间隔的特征波段个数对变异系数的影响,由图4可以看出所示,特征波段个数在10处取得了最高的变异系数。其中灰色关联度法中,首先,将多量纲数据转化成统一的无量纲形式,定义参考数列为个味觉信息数列,比较数列为380-1038nm各波段下的光谱信息数列来进行无量纲处理,以此消除味觉与光谱信息量纲的差异。然后,求取各单一味觉信息数列与全波段下的光谱信息数列的灰色关联系数。最后,求取各波段的灰色关联度。结果如表1所示,各单一味觉信息前10个关联度值较高的特征波段。Spectral images reflect the differences in the internal quality of apples, and the electronic tongue obtains taste information through the free functional groups in apple samples. Therefore, the method of gray relational degree (GRA) is used to determine the degree of correlation between the two, and the sensitive bands corresponding to the three basic tastes of sweet, sour and salty are selected respectively. Before gray correlation analysis, the coefficient of variation was used to determine the optimal number of sensitive bands. The larger the comprehensive coefficient of variation, the lower the correlation between variables. To study the influence of the number of characteristic bands with an interval of 2 on the coefficient of variation, it can be seen from Figure 4 that the number of characteristic bands achieved the highest value at 10. coefficient of variation. Among them, in the gray relational degree method, firstly, the multi-dimensional data is converted into a unified dimensionless form, the reference sequence is defined as a taste information sequence, and the comparison sequence is the spectral information sequence under each band of 380-1038nm for dimensionless processing. This eliminates the difference in the dimensions of gustatory and spectral information. Then, calculate the gray correlation coefficient between each single taste information series and the spectral information series under the full band. Finally, the gray correlation degree of each band is calculated. The results are shown in Table 1, the top 10 characteristic bands with higher correlation value for each single taste information.

表1 各味觉前10特征波段Table 1 The top 10 characteristic bands of each taste

S4、苹果味觉可视化模型的构建S4. Construction of Apple Taste Visualization Model

光谱与味觉信息数据具有离散性,且二者之间呈现局部线性与非线性关系,而SVM主要思想是建立一个回归超平面作为决策面,利用核函数将多维数据映射到高维空间,使其尽可能线性,以解决原始数据的局部非线性,最终使集合中所有数据到超平面的距离最近。基于此本文采用SVM预测各个味觉信息离散点的浓度值,其中超平面函数、RBF核函数及回归函数公式(2)、(3)、(4)分别如下。SVR建模过程结果依赖于参数c和g,正确有效的参数选择对支持向量机具有良好的回归性能。因此,本文基于遗传算法(GA)和网格式搜索(GS)的方法来对参数c和g进行寻优,如图5所示。在遗传算法中,最大的遗传代数为100,初始种群数为20,参数c的搜索范围是0到100,g为0到100。在网格式搜索方法中,以0.5为间隔进行参数寻优,参数c和g的搜索范围是2-10到210。在建立模型的过程中,因实验所检测为整体样品的味觉值,故在进行光谱数据和味觉信息映射时,将每个波段下的光谱值进行均值求取后再进行模型的建立。其中利用Kennard-Stone方法将2/3个光谱数据样本选为训练集,1/3个光谱数据选为预测集,基于GS-SVM和GA-SVM实现光谱-味觉信息的可视化分析。结果如表2显示GA-SVM的参数寻优过程,遗传算法快速的在达到最佳的适应度Mse为0.01667条件下筛选出了建模最优参数c为6.3859,g为83.0159,GS-SVM的参数寻优过程,网格式搜索在达到最佳的适应度Mse为0.0037346条件下筛选出了建模最优参数c为22.6274,g为0.0019531,因此选择GS-SVM作为苹果味觉值求取的模型。Spectral and taste information data are discrete, and there is a local linear and nonlinear relationship between them, and the main idea of SVM is to establish a regression hyperplane as a decision surface, and use the kernel function to map multidimensional data to high-dimensional space, making it Be as linear as possible to solve the local nonlinearity of the original data, and finally make the distance between all the data in the set and the hyperplane the shortest. Based on this, this paper uses SVM to predict the concentration value of each discrete point of taste information, in which the hyperplane function, RBF kernel function and regression function formulas (2), (3), and (4) are as follows. The result of the SVR modeling process depends on the parameters c and g, and the correct and effective parameter selection has good regression performance for the support vector machine. Therefore, this paper optimizes the parameters c and g based on genetic algorithm (GA) and grid search (GS), as shown in Figure 5. In the genetic algorithm, the maximum genetic algebra is 100, the initial population number is 20, the search range of parameter c is 0 to 100, and g is 0 to 100. In the grid search method, parameter optimization is carried out at an interval of 0.5, and the search range of parameters c and g is 2 -10 to 2 10 . In the process of building the model, because the taste value detected by the experiment is the taste value of the whole sample, when mapping the spectral data and taste information, the average value of the spectral value under each band is calculated before the model is established. The Kennard-Stone method is used to select 2/3 of the spectral data samples as the training set, and 1/3 of the spectral data as the prediction set. Based on GS-SVM and GA-SVM, the visual analysis of spectral-taste information is realized. The results are shown in Table 2, which shows the parameter optimization process of GA-SVM. The genetic algorithm quickly screened out the optimal modeling parameters c as 6.3859 and g as 83.0159 under the condition that the optimal fitness Mse was 0.01667. The GS-SVM In the parameter optimization process, the grid search screened out the optimal modeling parameters c as 22.6274 and g as 0.0019531 under the condition that the optimal fitness Mse was 0.0037346, so GS-SVM was selected as the model for calculating the apple taste value.

K(x,x′)=exp(g Px-x′P2) (2)K(x, x')=exp(g Px-x'P 2 ) (2)

f(x)=wφ(x)+b (3)f(x)=wφ(x)+b (3)

式中w-超平面的法向向量,φ(x)-非线性映射函数,b-偏置量,g-宽度系数In the formula, w-the normal vector of the hyperplane, φ(x)-non-linear mapping function, b-bias, g-width coefficient

表2 c与g寻优过程的结果值Table 2 Result values of c and g optimization process

S5、味觉可视化呈现S5. Taste visual presentation

在GS-SVM模型建立之后,将10个波段下的预测光谱值按点输入,求出各点光谱值下的味觉预测输出值记为K(i,j),将味觉值转化为对应的实际色彩值,将苹果味觉信息可视化呈现,图6为甜味、酸味、咸味单点味觉分布图,图例从深蓝色到砖红色味道逐渐变浓。由图可以看出,阿克苏糖心苹果主要以甜味、酸味为主,咸味很淡,且甜味最浓且分布最广。甜味主要分布在赤道线两侧,且中轴线部分较甜;酸味的分布位置与甜味相似;咸味则是在赤道截面较为均匀且分布更为广泛。After the GS-SVM model is established, the predicted spectral values under the 10 bands are input point by point, and the taste prediction output value under the spectral values of each point is calculated as K(i, j), and the taste value is converted into the corresponding actual value. Color value, which visualizes the taste information of apples. Figure 6 is a single-point taste distribution map of sweetness, sourness, and saltiness. The taste of the legend gradually becomes stronger from dark blue to brick red. It can be seen from the picture that Aksu candied apples are mainly sweet and sour, with a very light salty taste, and the sweetest and most widely distributed. The sweet taste is mainly distributed on both sides of the equator line, and the central axis is sweeter; the distribution position of the sour taste is similar to the sweet taste; the salty taste is more uniform and widely distributed in the equator section.

以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。The above is only a preferred embodiment of the present invention, it should be pointed out that for those of ordinary skill in the art, without departing from the principle of the present invention, some improvements and modifications can also be made, and these improvements and modifications should also be It is regarded as the protection scope of the present invention.

Claims (4)

1.一种苹果味觉图谱信息可视化方法,其特征在于,包括如下步骤:1. an apple taste map information visualization method, is characterized in that, comprises the steps: S1、基于电子舌苹果样本味觉数据的获取S1. Acquisition of taste data of apple samples based on electronic tongue S11、选取同一批次、同一地点采摘的果形正、大小均匀、无缺陷或污染物的苹果,按国标GB/T10651-2008选出90个优等果品;S11. Select apples picked in the same batch and at the same location with positive shape, uniform size, no defects or pollutants, and select 90 high-quality fruits according to the national standard GB/T10651-2008; S12、将样品编号后置于常温20±2℃,相对湿度55±5%贮藏24小时,保持温湿度不变;S12. Store the sample at room temperature 20±2°C and relative humidity 55±5% for 24 hours after numbering, keeping the temperature and humidity unchanged; S13、将每个苹果样本清洗、削皮、榨汁、过滤取上层清液40mL分别置于2个纯净测量杯中,按编号顺序依次经电子舌的检测得到30秒处90×3维糖度、酸味、成味三种味觉数据;测试开始前,传感器先在正负极清洗溶液中清洗90s,结束后在参比溶液中清洗120s,在另一种参比溶液中继续清洗120s,传感器平衡归零30s;达到平衡后,开始进行检测,测试时间为30s,每次测量结束后自动进入清洗步骤;S13. Wash, peel, squeeze, and filter each apple sample to take 40 mL of the supernatant and place them in two pure measuring cups, and obtain 90 × 3-dimensional sugar content in 30 seconds through the detection of the electronic tongue in sequence according to the numbering order. Three kinds of taste data: sour taste and mature taste; before the test starts, the sensor is cleaned in the positive and negative electrode cleaning solution for 90s, after the end of the test, it is washed in the reference solution for 120s, and in another reference solution for 120s, and the sensor is balanced. Zero 30s; after the balance is reached, the test starts, the test time is 30s, and the cleaning step is automatically entered after each measurement; S2、基于高光谱技术苹果样本数据的获取与预处理S2. Acquisition and preprocessing of apple sample data based on hyperspectral technology S21、通过下式采用黑白标定方法对高光谱图像进行校正,以消除噪声的影响:S21. Correct the hyperspectral image by using the black-and-white calibration method through the following formula to eliminate the influence of noise: 式中,Rd-暗图像,Rw-白板的漫反射图像,RS-苹果样本原始的漫反射光谱图像,R-校正后的漫反射光谱图。In the formula, R d - dark image, R w - diffuse reflectance image of whiteboard, R S - original diffuse reflectance spectrum image of apple sample, R - corrected diffuse reflectance spectrum image. S22、对黑白标定后的图像进行矢量化处理,获得图像描述曲线,圈定目标图像后,生成特定感兴趣掩膜图像,控制样本影像处理区域;S22. Carry out vectorization processing on the black-and-white calibrated image, obtain an image description curve, delineate the target image, generate a specific mask image of interest, and control the sample image processing area; S23、对掩膜处理后的图像进行影像裁剪,掩膜后的光谱区域除感兴趣处外其他光谱数值均为0。S23. Perform image cropping on the masked image, and the spectral values of the masked spectral region are all 0 except for the region of interest. S24、获取预处理后的90个苹果光谱图像,其中,由于苹果为类球形水果,故在用高光谱分选仪进行照射样品时将一个苹果样本划分为4个面,并在每个面上分别取编号1-5的5个感兴趣区域,每个感兴趣区域大小为300像素点,求取每个样品5个表面整体的光谱均值,最终得到90×256维的苹果数据;S24. Obtain 90 apple spectral images after preprocessing, wherein, since apples are spherical fruits, an apple sample is divided into 4 surfaces when a hyperspectral sorter is used to irradiate the sample, and each surface Take 5 regions of interest numbered 1-5 respectively, each region of interest is 300 pixels in size, calculate the average value of the spectra of the 5 surfaces of each sample, and finally obtain 90×256 dimensional apple data; S3、特征波段的选取S3. Selection of characteristic bands 分别选取甜、酸、咸3种基础味觉相对应的敏感波段,采用变异系数确定最佳敏感波段数量,然后采用灰色关联度(GRA)的方法确定苹果样品的味觉信息与光谱图像之间的关联程度;The sensitive bands corresponding to the three basic tastes of sweet, sour and salty were respectively selected, the optimal number of sensitive bands was determined by the coefficient of variation, and then the correlation between the taste information of the apple sample and the spectral image was determined by the gray relational degree (GRA) method degree; S4、苹果味觉可视化模型的构建S4. Construction of Apple Taste Visualization Model 采用SVM预测各个味觉信息离散点的浓度值,其中超平面函数、RBF核函数及回归函数公式(2)、(3)、(4)分别如下:Using SVM to predict the concentration value of each discrete point of taste information, the hyperplane function, RBF kernel function and regression function formulas (2), (3) and (4) are as follows: K(x,x′)=exp(g Px-x′P2) (2)K(x, x')=exp(g Px-x'P 2 ) (2) f(x)=ωφ(x)+b (3)f(x)=ωφ(x)+b (3) 式中,w-超平面的法向向量,φ(x)-非线性映射函数,b-偏置量,g-宽度系数;In the formula, w-the normal vector of the hyperplane, φ(x)-non-linear mapping function, b-bias, g-width coefficient; S5、味觉可视化呈现S5. Taste visual presentation 在GS-SVM模型建立之后,将10个波段下的预测光谱值按点输入,求出各点光谱值下的味觉预测输出值记为K(i,j),将味觉值转化为对应的实际色彩值,将苹果味觉信息可视化呈现。After the GS-SVM model is established, the predicted spectral values under the 10 bands are input point by point, and the taste prediction output value under the spectral values of each point is calculated as K(i, j), and the taste value is converted into the corresponding actual value. Color value, which visualizes the taste information of apples. 2.如权利要求1所述的一种苹果味觉图谱信息可视化方法,其特征在于,2. a kind of apple taste map information visualization method as claimed in claim 1, is characterized in that, 所述灰色关联度法具体包括如下步骤:首先,将多量纲数据转化成统一的无量纲形式,定义参考数列为个味觉信息数列,比较数列为380-1038nm各波段下的光谱信息数列来进行无量纲处理,以此消除味觉与光谱信息量纲的差异;然后,求取各单一味觉信息数列与全波段下的光谱信息数列的灰色关联系数;最后,求取各波段的灰色关联度。The gray relational degree method specifically includes the following steps: first, convert the multidimensional data into a unified dimensionless form, define a reference sequence as a taste information sequence, and compare the sequence as a spectral information sequence under each band of 380-1038nm to carry out infinite Then, the gray correlation coefficients between each single taste information sequence and the spectral information sequence under the full band are calculated; finally, the gray correlation degree of each band is calculated. 3.如权利要求1所述的一种苹果味觉图谱信息可视化方法,其特征在于,所述步骤S4基于遗传算法(GA)和网格式搜索(GS)的方法来对参数c和g进行寻优,在遗传算法中,最大的遗传代数为100,初始种群数为20,参数c的搜索范围是0到100,g为0到100;在网格式搜索方法中,以0.5为间隔进行参数寻优,参数c和g的搜索范围是2-10到2103. a kind of apple taste map information visualization method as claimed in claim 1, is characterized in that, described step S4 carries out optimization to parameter c and g based on the method of genetic algorithm (GA) and grid search (GS) , in the genetic algorithm, the maximum genetic algebra is 100, the initial population number is 20, the search range of the parameter c is 0 to 100, and the g is 0 to 100; in the grid search method, the parameters are optimized at an interval of 0.5 , the search range of parameters c and g is 2 -10 to 2 10 . 4.如权利要求1所述的一种苹果味觉图谱信息可视化方法,其特征在于,所述步骤S4中在建立模型的过程中,在进行光谱数据和味觉信息映射时,将每个波段下的光谱值进行均值求取后再进行模型的建立;其中利用Kennard-Stone方法将2/3个光谱数据样本选为训练集,1/3个光谱数据选为预测集,基于GS-SVM和GA-SVM实现光谱-味觉信息的可视化分析。4. a kind of apple taste map information visualization method as claimed in claim 1, is characterized in that, in the process of building a model in described step S4, when carrying out spectral data and taste information mapping, the The average value of the spectral value is calculated before the model is established; 2/3 of the spectral data samples are selected as the training set by the Kennard-Stone method, and 1/3 of the spectral data are selected as the prediction set, based on GS-SVM and GA- SVM realizes visual analysis of spectral-taste information.
CN201810435827.3A 2018-04-28 2018-04-28 Apple taste map information visualization method Active CN108760655B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810435827.3A CN108760655B (en) 2018-04-28 2018-04-28 Apple taste map information visualization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810435827.3A CN108760655B (en) 2018-04-28 2018-04-28 Apple taste map information visualization method

Publications (2)

Publication Number Publication Date
CN108760655A true CN108760655A (en) 2018-11-06
CN108760655B CN108760655B (en) 2021-03-30

Family

ID=64009312

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810435827.3A Active CN108760655B (en) 2018-04-28 2018-04-28 Apple taste map information visualization method

Country Status (1)

Country Link
CN (1) CN108760655B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109685134A (en) * 2018-12-21 2019-04-26 龙口盛福达食品有限公司 A kind of baby milk intelligence gustatory analysis method based on unsupervised discretization
CN110763698A (en) * 2019-10-12 2020-02-07 仲恺农业工程学院 A hyperspectral citrus leaf disease identification method based on characteristic wavelengths
CN111289446A (en) * 2020-03-30 2020-06-16 天津工业大学 Method and system for detecting component concentration of complex solution
CN111289516A (en) * 2020-03-26 2020-06-16 中国农业大学 Method and device for detecting amino acid content of plant leaves
CN112285054A (en) * 2020-09-08 2021-01-29 塔里木大学 Apple identity recognition model establishing method and system and identity recognition method and system
EP4070076A4 (en) * 2019-12-04 2024-01-17 Profileprint Pte. Ltd. SYSTEM AND METHOD FOR NON-DESTRUCTIVE RAPID FOOD PROFILING USING ARTIFICIAL INTELLIGENCE

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060203247A1 (en) * 2001-07-09 2006-09-14 Mitsui Mining & Smelting Co., Ltd. Apparatus and method for evaluating the interior quality of fruits and vegetables
CN1995987A (en) * 2007-02-08 2007-07-11 江苏大学 Non-destructive detection method and device for agricultural and animal products based on hyperspectral image technology
CN101021478A (en) * 2007-02-13 2007-08-22 江西农业大学 Laser inducing fluorescent high spectral image detecting method and device for fruit quality
CN201021922Y (en) * 2007-02-13 2008-02-13 江西农业大学 A device for non-destructive testing of fruit quality
CN101832926A (en) * 2010-03-19 2010-09-15 江南大学 Method for performing apple powder materialization non-destructive inspection by using hyper-spectral image technique
CN101853392A (en) * 2010-04-21 2010-10-06 河海大学 Band Selection Method for Remote Sensing Hyperspectral Imagery Based on Conditional Mutual Information
CN105738311A (en) * 2015-12-01 2016-07-06 昆明睿意铂科技股份有限公司 Apple sweetness non-damage quick detection method based on near-infrared spectrum technology
CN105866043A (en) * 2016-06-06 2016-08-17 东北电力大学 Method for detecting apple sour through hyperspectral technology

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060203247A1 (en) * 2001-07-09 2006-09-14 Mitsui Mining & Smelting Co., Ltd. Apparatus and method for evaluating the interior quality of fruits and vegetables
CN1995987A (en) * 2007-02-08 2007-07-11 江苏大学 Non-destructive detection method and device for agricultural and animal products based on hyperspectral image technology
CN101021478A (en) * 2007-02-13 2007-08-22 江西农业大学 Laser inducing fluorescent high spectral image detecting method and device for fruit quality
CN201021922Y (en) * 2007-02-13 2008-02-13 江西农业大学 A device for non-destructive testing of fruit quality
CN101832926A (en) * 2010-03-19 2010-09-15 江南大学 Method for performing apple powder materialization non-destructive inspection by using hyper-spectral image technique
CN101853392A (en) * 2010-04-21 2010-10-06 河海大学 Band Selection Method for Remote Sensing Hyperspectral Imagery Based on Conditional Mutual Information
CN105738311A (en) * 2015-12-01 2016-07-06 昆明睿意铂科技股份有限公司 Apple sweetness non-damage quick detection method based on near-infrared spectrum technology
CN105866043A (en) * 2016-06-06 2016-08-17 东北电力大学 Method for detecting apple sour through hyperspectral technology

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘晶晶 等: "高光谱图像技术检测苹果外部损伤", 《中国食品学报》 *
宫元娟 等: "寒富苹果品质无损检测光谱信息在线分析", 《沈阳农业大学学报》 *
张偲敏 等: "支持向量回归机的参数择优算法", 《中国科技信息》 *
马惠玲 等: "基于高光谱成像的苹果品种快速鉴别", 《农业机械学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109685134A (en) * 2018-12-21 2019-04-26 龙口盛福达食品有限公司 A kind of baby milk intelligence gustatory analysis method based on unsupervised discretization
CN110763698A (en) * 2019-10-12 2020-02-07 仲恺农业工程学院 A hyperspectral citrus leaf disease identification method based on characteristic wavelengths
EP4070076A4 (en) * 2019-12-04 2024-01-17 Profileprint Pte. Ltd. SYSTEM AND METHOD FOR NON-DESTRUCTIVE RAPID FOOD PROFILING USING ARTIFICIAL INTELLIGENCE
CN111289516A (en) * 2020-03-26 2020-06-16 中国农业大学 Method and device for detecting amino acid content of plant leaves
CN111289516B (en) * 2020-03-26 2021-10-08 中国农业大学 Method and device for detecting amino acid content in plant leaves
CN111289446A (en) * 2020-03-30 2020-06-16 天津工业大学 Method and system for detecting component concentration of complex solution
CN112285054A (en) * 2020-09-08 2021-01-29 塔里木大学 Apple identity recognition model establishing method and system and identity recognition method and system

Also Published As

Publication number Publication date
CN108760655B (en) 2021-03-30

Similar Documents

Publication Publication Date Title
CN108760655A (en) A kind of apple sense of taste profile information method for visualizing
CN111968080B (en) A detection method for internal and external quality of Feicheng peach based on hyperspectral and deep learning
US11454594B2 (en) Method for identifying frostbite condition of grain seeds using spectral feature wavebands of seed embryo hyperspectral images
CN108875913B (en) Tricholoma matsutake rapid nondestructive testing system and method based on convolutional neural network
CN110781889B (en) Deep learning-based nondestructive testing method for total sugar content in blueberry fruits
CN110349153A (en) The image method tobacco based on multi-feature fusion detection method of rate containing stalk
WO2020186746A1 (en) Electronic tongue-based soy sauce taste evaluation method and device, medium and calculation apparatus
CN106841053B (en) A kind of rapid identification method of jujube variety
CN109827957A (en) A kind of rice leaf SPAD value estimating and measuring method based on computer vision and system
CN105181606B (en) The method that cane sugar content is distributed in peanut is detected based on high light spectrum image-forming technology
CN114140695A (en) A prediction method and system for nitrogen diagnosis and quality index determination of tea tree based on UAV multispectral remote sensing
CN110516668A (en) A method and device for detecting honey adulteration based on hyperspectral imaging technology
CN116482317B (en) Real-time monitoring methods, systems, equipment and media for lake water nutrient status
Mu et al. Non‐destructive detection of blueberry skin pigments and intrinsic fruit qualities based on deep learning
Patel et al. Rapid assessment of some physical parameters of mangoes using monochrome computer vision
CN106841167A (en) The lossless detection method of garden stuff pesticide residue
Yu et al. Nondestructive testing of external defects on Nanguo pear
CN104568639B (en) A kind of determination method and apparatus of sugar degree
CN114965346A (en) Kiwi fruit quality detection method based on deep learning and hyperspectral imaging technology
CN115545383A (en) Binocular stereoscopic vision detection method for crisp appearance quality of strawberries
CN104198397A (en) Method for detecting nutrient content of greenhouse crops under nitrogen, phosphorus and potassium interaction
CN117288692A (en) Method for detecting tannin content in brewing grains
CN106442397A (en) Rice near infrared spectrum model optimizing method based on spectrum proportional deduction
CN106124447A (en) A kind of based on the method for soluble solid content in near-infrared spectral analysis technology detection Fructus Fragariae Ananssae
CN104574341B (en) A kind of determination method and apparatus of sugar degree

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
GR01 Patent grant
GR01 Patent grant