CN113552089A - 一种红茶萎凋程度数字化的无损检测方法 - Google Patents
一种红茶萎凋程度数字化的无损检测方法 Download PDFInfo
- Publication number
- CN113552089A CN113552089A CN202110765625.7A CN202110765625A CN113552089A CN 113552089 A CN113552089 A CN 113552089A CN 202110765625 A CN202110765625 A CN 202110765625A CN 113552089 A CN113552089 A CN 113552089A
- Authority
- CN
- China
- Prior art keywords
- withering
- spectrum
- withered
- sample
- degree
- 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
- 244000269722 Thea sinensis Species 0.000 title claims abstract description 53
- 235000006468 Thea sinensis Nutrition 0.000 title claims abstract description 41
- 235000020279 black tea Nutrition 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000009659 non-destructive testing Methods 0.000 title claims abstract description 11
- 238000001228 spectrum Methods 0.000 claims abstract description 46
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 39
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 8
- 238000013528 artificial neural network Methods 0.000 claims description 31
- 230000003595 spectral effect Effects 0.000 claims description 23
- 238000012549 training Methods 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 10
- 238000000513 principal component analysis Methods 0.000 claims description 10
- 238000009499 grossing Methods 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 5
- 238000012795 verification Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000004891 communication Methods 0.000 claims description 4
- 238000010438 heat treatment Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000003860 storage Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 16
- 238000004458 analytical method Methods 0.000 abstract description 6
- 235000013616 tea Nutrition 0.000 description 12
- 238000003062 neural network model Methods 0.000 description 10
- 238000010238 partial least squares regression Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 238000012545 processing Methods 0.000 description 7
- 238000010200 validation analysis Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 210000004027 cell Anatomy 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 210000002569 neuron Anatomy 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 3
- 235000013305 food Nutrition 0.000 description 3
- 238000012886 linear function Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 description 2
- 238000004497 NIR spectroscopy Methods 0.000 description 2
- 230000004913 activation Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000000855 fermentation Methods 0.000 description 2
- 230000004151 fermentation Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 230000001953 sensory effect Effects 0.000 description 2
- 238000010561 standard procedure Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 235000009024 Ceanothus sanguineus Nutrition 0.000 description 1
- 241000037488 Coccoloba pubescens Species 0.000 description 1
- 240000003553 Leptospermum scoparium Species 0.000 description 1
- 235000015459 Lycium barbarum Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 235000013361 beverage Nutrition 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001599 direct drying Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000796 flavoring agent Substances 0.000 description 1
- 235000019634 flavors Nutrition 0.000 description 1
- 230000007407 health benefit Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 235000003715 nutritional status Nutrition 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010998 test method 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/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
- 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/01—Arrangements or apparatus for facilitating the optical investigation
-
- 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/01—Arrangements or apparatus for facilitating the optical investigation
- G01N2021/0106—General arrangement of respective parts
- G01N2021/0112—Apparatus in one mechanical, optical or electronic block
Landscapes
- Physics & Mathematics (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)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种红茶萎凋程度数字化的无损检测方法，包括以下步骤：S1、萎凋叶样品光谱的采集，将处于不同萎凋时间节点的萎凋叶样品，分别通过连接智能手机的微型近红外光谱仪采集，进而得到不同萎凋时间节点下的样品光谱；所述萎凋时间节点包括过轻萎凋、轻萎凋、适度萎凋Ⅰ、适度萎凋Ⅱ、重萎凋、过重萎凋；S2、样品萎凋程度数字化参数的检测，随着萎凋的进行，检测不同萎凋时间节点萎凋叶的萎凋程度数字化参数。本发明将近红外光谱应用于萎凋叶含水率的快速无损检测，提出的基于扩展输入空间的Elman模型，具有更好的定量预测分析能力，相对于传统的线性PLS模型和非线性BP模型具有更强的预测能力。
Description
技术领域
本发明属于食品加工技术领域，尤其涉及一种红茶萎凋程度数字化的无损检测方法。
背景技术
红茶是全球最受欢迎的饮品之一，具有独特的风味和多种健康益处。红茶的初制涉及多个过程，包括萎凋、揉捻、发酵和干燥等四道工序。萎凋是红茶初加工中的关键工序，叶内水分和品质生化成分发生剧烈变化，为红茶品质的形成奠定重要基础，对红茶加工和生产具有重要意义。
在萎凋的过程中，鲜叶适当散失水分，降低鲜叶细胞表面张力，使芽叶由脆变软，增强可塑性、柔软性等物理特性，为揉捻工序获得一定细胞破碎率和初步做形作业创造必要条件，所以生产中将萎凋叶的含水率是判别萎凋适度的关键指标，现有研究认为当鲜叶的含水率达到60％-64％为适度萎凋。目前，对于萎凋程度的传统判断方法仍依靠人的感官经验判断，容易受到人为因素的影响，而传统的含水率检测方法操作复杂、检测周期长，难以实现萎凋叶含水率的快速、准确、便捷、无损地检测。近红外光谱法(NIRS)是一种先进的无损检测技术，已广泛应用于茶叶的定性和定量的分析，如茶的等级分类、发酵度预测、感官质量确定、主要成分等。但是，先前的研究中使用的大部分都是基于实验室条件下的台式近红外光谱仪，这种光谱仪价格高、体积大且不可携带，在实验室以外的生产中难以使用。此外，这种台式光谱仪结构复杂、需要连接电脑且操作技能要求高，难以成为茶厂和农民的常规生产工具。
随着计算机、智能手机和光学仪器的快速普及发展，已经开发了通过无线蓝牙技术连接到智能手机的微型近红外光谱仪，可利用智能手机中的应用程序操控光谱仪完成光谱数据的采集、存储和传输，再通过微型计算机对数据完成分析、处理和显示。由于微型计算机和微型近红外光谱仪价格低、便携性且对使用环境
要求低，基于智能手机的微型近红外光谱系统得到了研究人员的青睐。微型近红外光谱仪已成功用于茶树的营养状况诊断、茶叶的质量评估、红茶的品味评估。然而，通过使用微型近红外光谱仪结合微型计算机进行红茶萎调中含水率在线定量分析的研究暂无。
发明内容
为了解决上述现有技术方案中的不足，本发明提供了一种红茶萎凋程度数字化的无损检测方法，该方法实现了红茶萎凋叶含水率快速无损的在线检测。
具体技术方案如下：
一种红茶萎凋程度数字化的无损检测方法，包括以下步骤：
S1、萎凋叶样品光谱的采集，将处于不同萎凋时间节点的萎凋叶样品，分别通过连接智能手机的微型近红外光谱仪采集，进而得到不同萎凋时间节点下的样品光谱；所述萎凋时间节点包括过轻萎凋、轻萎凋、适度萎凋Ⅰ、适度萎凋Ⅱ、重萎凋、过重萎凋；
S2、样品萎凋程度数字化参数的检测，随着萎凋的进行，检测不同萎凋时间节点萎凋叶的萎凋程度数字化参数；
S3、对采集的萎凋叶样品光谱进行预处理和主成分分析，降低采集光谱时萎凋叶样品表面漫反射对光谱信息造成的干扰，提取光谱数据主要信息；
S4、Elman神经网络的建立，将S3进行预处理后的光谱数据融合，通过扩展输入空间的方法，按照5：1随机划分为校正集和预测集，校正集与S2步骤的不同萎凋时间节点下萎凋叶相对应的萎凋程度数字化参数用于训练Elman神经网络的红茶萎凋水分定量预测模型，预测集作为外部验证样本，进行萎凋程度数字化参数的预测；根据预测集样本所预测的萎凋程度数字化参数与实际值的相关系数分析红茶萎凋水分定量预测模型的性能；
S5、不同萎凋程度数字化参数的快速检测，通过无线蓝牙技术连接到智能手机的微型近红外光谱仪，利用智能手机中的应用程序操控微型近红外光谱仪完成萎凋叶样品光谱数据的采集和存储，通过USB将存储的光谱数据传输到PC端，通过MATLAB软件建立水分预测模型；
S6、将PC端建立好的水分预测模型部署到微型计算机Raspberry Pi上独立运行，通过微型近红外光谱仪、OLED显示屏与微型计算机Raspberry Pi通信后，光谱数据即可实时传输至已部署好水分预测算法模型的微型计算机Raspberry Pi上，通过模型运行计算后，萎凋叶水分预测结果即可通过OLED显示屏实时显示。
优选地，所述S1步骤中，从红茶萎凋开始至萎凋过度后1小时止，共18h；期间每隔1h测取该萎凋程度下萎凋叶的光谱数据。
优选地，所述微型近红外光谱仪的近红外波段范围为900-1700nm，分辨率为10-12nm；光谱仪的尺寸为75mm×58mm×26.5mm，重约77g，智能手机通过无线蓝牙连接到微型近红外光谱仪，此种光谱仪具有低成本且极易便携的特点；微型计算机Raspberry Pi为4B主板，8GB内存，其尺寸为85.60mm×56mm×21mm，重约45g，数据处理性能好且易便携；OLED显示屏为0.96寸，分辨率128*64，控制芯片为SSD1306，通信方式为IIC。
优选地，所述S1步骤中光谱采集的条件为将每个萎凋程度下萎凋叶样品轻轻按压均匀的表面并扫描3次，然后将平均光谱用作样品的代表性光谱数据。
优选地，所述S2步骤中描述红茶萎凋程度数字化的指标是萎凋叶的含水率，所述含水率采用水分测定仪的加热计算差值测定，所述含水率采用GB 5009.3—2016进行测定。
优选地，所述S3步骤中，对近红外光谱采用平滑预处理，对预处理后的光谱进行PCA主成分分析。
优选地，所述PCA主成分分析的主成分数为10。
优选地，所述Elman神经网络结构包括2个输入层、9个隐藏层、1个输出层，它比普通的神经网络在结构上多一个承接层，承接层用于承接整个网络上一次的状态，可看作神经网络的记忆延时单元；承接层的增益因子取[0,1]，用来确保Elman神经网络的收敛性。
本发明的有益效果是：
1.创新性地将近红外光谱应用于萎凋叶含水率的快速无损检测。在萎凋的过程中，萎凋叶内部发生了一系列的物理变化和化学变化，在此过程中，鲜叶会不断失去水分而变得萎蔫柔软，水分的变化幅度较大，而近红外光谱可以描述有机分子中-OH基团的响应，能表征样品的内部信息。所以通过近红外光谱结合萎凋叶含水率，基于Elman神经网络建立红茶萎凋程度的判别模型，为红茶萎凋程度提供一种科学准确、适合在线、快速无损检测的判别方法；
2.创新性地将光谱数据进行平滑预处理和主成分分析多重分析，提出了数据扩展输入空间的思想，并应用扩展输入空间建立了红茶萎凋程度的判别的Elman神经网络模型。本发明中便携式的红茶萎凋程度数字化的快速无损检测方法所建立的Elman神经网络萎凋叶含水率预测模型，在获取萎凋叶样品的近红外光谱后，将原始光谱经过光滑预处理和主成分分析后的数据共同用作模型输入，最大程度地去除了噪声干扰，通过扩展输入空间的方法增加了模型的泛化性，将扩展输入空间后的数据与对应萎凋叶的含水率参与模型训练中，因此该发明提出的基于扩展输入空间的Elman模型具有更好的定量预测分析能力，相对于传统的线性PLS模型和非线性BP模型具有更强的预测能力；
3.发明了一套基于智能手机通过无线蓝牙技术结合微型近红外光谱仪、微型计算机Raspberry Pi(树莓派)进行萎凋叶叶含水率的快速显示系统。本发明中所使用的通过无线蓝牙技术连接到智能手机的微型近红外光谱仪与微型计算机Raspberry Pi(树莓派)通信后，光谱数据即可实时传输至已部署好水分预测模型的微型计算机Raspberry Pi(树莓派)上，通过已部署好的模型运行计算后，即可将萎凋叶水分预测结果通过OLED显示屏实时显示，从而实现红茶萎凋程度数字化的定量快速无损检测。此套萎凋叶含水率快速无损检测设备具有价格低、可便携、易操作，可在线快速检测萎凋叶含水率的优点，可成为茶厂和茶农加工茶叶的常规生产工具。
附图说明
图1是本发明使用的萎凋叶数字化水分预测系统；
图2是本发明基于Elman神经网络的红茶萎凋含水率定量预测模型的算法流程图；
图3是本发明Elman神经网络的结构图；
图4是本发明实施例1中经过微型近红外光谱仪采集的原始光谱图像(a)以及经过平滑预处理之后的图像(b)；
图5是本发明实施例1中基于Elman的神经网络模型外部验证集实测值与预测值散点图(a)，外部验证集预测误差图(b)；
图6是本发明对照实施例1中基于PLS线性模型的主成分贡献率(a)，外部验证集实测值与预测值散点图(b)；
图7是本发明对照实施例1中基于BP神经网络非线性模型外部验证集实测值与预测值散点图(a)，外部验证集预测误差图(b)；
图中，1-Raspberry Pi；2-智能手机；3-样品池；4-微型近红外光谱仪；5-OLED显示屏。
具体实施方式
下面将各个附图和具体的实施例相结合进行详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
实施例1
一种红茶萎凋程度数字化的无损检测方法，包括如下步骤：
S1.萎凋叶样品光谱的采集，将处于不同萎凋程度节点(过轻萎凋、轻萎凋、适度萎凋Ⅰ、适度萎凋Ⅱ、重萎凋、过重萎凋等)的萎凋叶样品，分别通过连接智能手机2的微型近红外光谱仪4(NIR-S-G1；InnoSpectra公司，中国台湾)采集样品池3，波长范围:900-1700nm，光谱采样间隔:1.5nm，光谱分辨率:10nm，波长精度:±1nm，镜头紧贴样品池3中样品表面，利用智能手机2中的应用程序操控微型近红外光谱仪4完成萎凋光谱数据的采集和存储，通过串行总线(USB)将微型近红外光谱仪4中存储的光谱传输到PC端，进而得到不同萎凋程度节点下的光谱数据；
S2.样品萎凋程度数字化参数的检测，随着萎凋的进行，检测不同萎凋时间节点萎凋叶的萎凋程度数字化参数；
S3.对采集的样品光谱依次进行SNV平滑预处理和PCA主成分分析，降低采集光谱时样品表面漫反射对光谱信息造成的干扰和提取光谱数据主要信息；
S4.Elman神经网络的建立，将S3进行平滑预处理与主成分分析后的光谱数据融合，通过扩展输入空间的方法后，按照5：1随机划分为校正集和预测集，校正集与S2步骤的不同萎凋时间下萎凋叶相对应的含水率指标用于训练Elman神经网络的红茶萎凋含水率定量预测模型，验证集作为外部验证样本，进行萎凋程度指标的预测，根据预测集样本所预测的含水率与实际值的相关系数分析预测模型的性能；
S5.不同萎凋程度样品含水率的快速检测，通过无线蓝牙技术连接到智能手机2的微型近红外光谱仪4，通过MATLAB软件写入预测模型，近而实现红茶萎凋程度数字化的定量快速无损检测。
S6.将电脑PC端的含水率预测模型部署到微型计算机Raspberry Pi 1(树莓派)上独立运行，与OLED显示屏5通信后将萎凋叶含水率预测结果实时显示，从而实现红茶萎凋程度数字化的定量快速无损检测。
S5.不同萎凋程度含水率的快速检测。将预测集作为外部验证样本，进行萎凋程度指标的预测，根据预测集样本相关系数分析预测模型的性能。
S6.不同萎凋程度的数字化实时显示。将电脑PC端建立好的含水率预测模型部署到微型计算机Raspberry Pi 1(树莓派)上独立运行，通过微型近红外光谱仪4与微型计算机Raspberry Pi 1(树莓派)通信后，光谱即可实时传输至已部署好含水率预测模型的微型计算机Raspberry Pi 1(树莓派)上，通过模型运行计算后，萎凋叶含水率预测结果即可通过OLED显示屏5实时显示，从而实现红茶萎凋程度数字化的定量快速无损检测。
本发明使用的Elman神经网络模型的建立过程如下：
Elman模型的输入为u(t)＝[v(t),p(v(t))]，其中原始光谱经过SNV平滑预处理后为v(t)，v(t)经过主成分分析后为p(v(t))；
故存在一对一映射，W[v(t)]＝NN[v(t),h(v(t))]+ε。其中，NN(·)表示多层前馈Elman神经网络。ε为误差近似值，存在任意εN＞0,|ε|≤εN。W(v(t))是使用水分测定仪测定每一时间光谱节点下萎凋叶的含水率值，由此构建基于Elman神经网络的红茶萎凋水分定量预测模型框图如图1所示，u(t)＝[v(t),P(v(t))]和
分别是Elman神经网络模型的输入和输出。为简易起见，用y(t)和
分别表示W(v(t))和
Elman神经网络的结构图如图3所示。
然后，所提出的数学模型表示为：
x(k)＝f(ω2xc(k)+ω1u(k-1)) 1-(1)
xc(k)＝αxc(k-1)+x(k-1) 1-(2)
其中，α是上下文节点的增益因子。0≤α＜1用来确保Elman神经网络的收敛性，k＝2,3,…,n.n是样本数。f(·)和g(·)分别为S函数和线性函数。Elman神经网络的输出表示为：
目标函数定义为：
通过以下公式更新神经网络的权重：
其中
是将输入层中的节点j连接到隐藏层中的节点q的权重；
是将上下文层中的节点j连接到隐藏层中的节点l的权重；
是将隐藏层中的节点j连接到输出层中的节点i的权重。η是正数，代表权重的学习因子。
E相对于
的导数可以得出：
此处，
i＝1,2,…,m；j＝1,2,…,n.
同样，E关于
的导数可以计算为：
此处
可以通过以下公式计算：
当不考虑xc(k)对连接权重
的依赖性时，可以将1-(9)改写为：
E相对于
的导数可以得出：
此处，
q＝1,2,…,r.
本实施例采用的Elman神经网络各层神经元的数量是通过试凑法给定的，其输入节点、隐藏节点和输出节点数分别为2、9、1。隐藏层和输出层的激活函数分别为单极性S函数和线性函数。Levenberg-Marquardt算法用于训练Elman神经网络。
本实施例采用如图1所示的光谱采集获取系统，包括微型计算机Raspberry Pi 1(树莓派)、智能手机2、样品池3、微型近红外光谱仪4和OLED显示屏5，智能手机2、微型近红外光谱仪4采集系统通过蓝牙连接进行信号传输，微型近红外光谱仪4作为红茶萎凋实验光谱取样系统，该光谱仪型号是NIR-S-G1(InnoSpectra公司，中国台湾)，波长范围:900-1700nm，光谱采样间隔:1.5nm，光谱分辨率:10nm，波长精度:±1nm。光谱数据采集完后以xsl格式保存至智能手机2，微型计算机Raspberry Pi 1(树莓派)用来部署Elman神经网络预测模型，型号为：4B主板，8GB内存，其尺寸为85.60mm×56mm×21mm，重约45g。OLED显示屏5为0.96寸，分辨率128*64，控制芯片为SSD1306，通信方式为IIC。
利用微型近红外光谱仪4采集萎凋样品的光谱数据，具体的光谱采集获取过程如下：将样品放置于样品池3中，通过蓝牙连接的智能手机2客户端发出采集信号控制微型近红外光谱仪4，微型近红外光谱仪4镜头紧贴样品池3中样品表面并扫描3次，然后将平均光谱用作样品的代表性光谱数据，进而得到不同萎凋程度节点下的光谱数据。样本光谱采集完之后，光谱数据会自动存储在智能手机2端，将以.xsl格式保存。光谱采集完毕后按照图2所示算法流程进行萎凋叶含水率预测模型的建立。
在实验数据处理的硬件方面，本实施例选用的是华为-Mate40 Pro智能手机2，光谱采集完后，将智能手机2内的数据通过USB数据线传输至电脑PC端，本次样本数据集基于MATLAB软件运行，然后将建立好的Elman神经网络预测模型部署到微型计算机RaspberryPi 1(树莓派)中，完成微型近红外光谱仪4、OLED显示屏5与Raspberry Pi 1的通信，从而实现对不同程度的萎凋叶实现含水率数字化的快速在线检测。
本实施例以萎凋过程中茶鲜叶的变化为研究对象，品种为云南大叶种，嫩度为1芽2叶。萎凋装置为萎凋槽，槽内温度控制在30℃，相对湿度控制在50％。从萎凋0h开始取样至萎凋过度，共18h，萎凋过程中每间隔1h取一次样品，总共18个萎凋时间节点，每个时间节点取20组样品进行光谱采集，并使用水分测定仪(MA35M-000230V1,Sartorious)测定每一时间节点下萎凋叶的含水率，该过程每个时间节点做3次重复，取每个时间节点萎凋叶的三次含水率的平均值作为该时间节点下样品对应的含水率。萎凋18h后，萎凋叶含水率由80.01％降至57.68％，共采集18个时间节点下的360组萎凋样品光谱。
将360个样品光谱进行平滑预处理、主成分分析、扩展输入空间，以起到减少数据的噪声，提取光谱主要信息，提高所建立模型的泛化能力及鲁棒性，接着按照5:1的比例随机分为校正集和预测集，即校正集305个样品，预测集55个样品。
萎凋叶含水率的Elman模型建立：
根据萎凋叶光谱及含水率数据集建立Elman模型，以下是模型的具体参数设定及说明：
训练次数及各层神经元的数量：理想情况下神经网络的训练终止条件是达到误差足够小，所需要对模型进行多少次的训练。训练次数过少或者过多均会对模型造成干扰，达不到预设的误差标准。在本实施例中，通过多次试验方法来确定，通过观察发现当训练次数为6000时，模型基本可以达到误差稳定。各层神经元的数量是通过试凑法给定的，其输入节点、隐藏节点和输出节点数分别为2、9、1。隐藏层和输出层的激活函数分别为单极性S函数和线性函数。Levenberg-Marquardt算法用于训练Elman神经网络，同时将学习率设置为0.0001。
为了能对所建模型做出客观评价，本实施例选用预测集中55个未参与Elman模型训练的萎凋叶光谱数据作为模型的验证集，以其相关系数(R)、均方根误差(RMSEP)、RPD(定义为标准偏差(SD)与RMSEP之比)作为模型预测性能的评价指标。
在萎凋叶的光谱数据集中，校正集中165个萎凋叶光谱数据经过6000次训练后，对所得模型进行外部预测集测试。提出模型的外部预测集相关系数Rp为0.9615，外部预测集均方根误差RMSE为0.3958，RPD为16.1132，表明本发明提出的模型对萎凋叶的含水率具有准确的预测能力。如图5(a)所示，我们展示了55个外部验证集预测值与实测值的散点图，从图中可以看出，模型在55个从未在训练集出现过的样本时，所得出预测值与实际值基本重合，且外部验证集相对误差均小于0.014，如图5(b)所示，说明模型具有较强的鲁棒性和泛化能力。
对照实施例1
将本发明提出的基于扩展输入空间的Elman神经网络模型和传统的线性模型偏最小二乘回归(PLS)、非线性模型BP神经网络模型的预测效果进行比较。关于萎凋叶含水率的偏最小二乘回归(PLS)与BP神经网络模型分析如下：
1.PLS线性模型
PLS作为典型的线性回归模型，可有效解决多重共线性问题。如图6(a)所示，表示
模型误差值与选取不同主成分数之间的关系，当主成分数NPC为10时，模型的RMSEV值低，所得模型的R p值为0.8755，RMSEP值为1.7751，RPD为2.5491，相对误差值均小于0.0840。表明建立的PLS模型对萎凋叶的含水率具有较好的预测功能。预测值与实测值散点分布如图6-(b)所示。
2.BP神经网络模型
BP神经网络模型作为经典的非线性模型，通过对样本模式的学习模拟信息之间的内在机制，实现输入与输出之间的高度非线性映射。当训练次数为1000时，模型基本可以达到误差稳定。各层神经元的数量是通过试凑法给定的，其输入节点、隐藏节点和输出节点数分别为1、4、1。同时将学习率设置为0.0001。模型的RMSEV值低，所得模型的Rp值为0.94951，RMSEP值为1.320，RPD为4.4732，相对对误差值均小于0.0541。表明建立的BP模型对萎凋叶的含水率具有较好的预测功能。预测值与实测值散点分布如图7(a)。
3.模型对比与分析
将线性模型PLS、非线性模型BP神经网络模型与本发明提出的Elman模型做对比，其结果如表1所示：
表1不同预测模型的结果对比
从预测集结果可以看出，三种模型的RPD值均大于2，表明对于萎凋叶含水率的预测均具有良好的预测能力，且基于扩展输入空间的Elman模型具有最大Rp和RPD值，是三种模型中的最优模型。此外由图5(a)、图6(b)、图7(a)模型外部预测集实际值与预测值可看出基于扩展输入空间的Elman模型相比其余两种模型有较强的鲁棒性和泛化能力，因此本发明可以实现便携式的红茶萎凋叶含水率的检测。
对照实施例2
以国标法检测结果为对照，对快速水分测定仪、台式近红外光谱仪、便携式近红外光谱仪等三种方式不同萎凋程度含水率测定效果进行比较，结果见表2。从表2可以看出：基于微型计算机Raspberry Pi 1(树莓派)和智能手机2操控的便携式近红外光谱仪检测含水率的时间最短，仅33s，同时含水率检测的精准度较其他2种显著高，达96.78％，并且检测不受条件限制，可应用于不同类型的作业场地，即证明了本发明获得的红茶萎凋程度数字化的快速无损检测方法可行、科学。
表2不同含水率检测方式结果对比
注：国标法采用GB 5009.3-2016《食品安全国家标准食品中的水分的测定》中的直接干燥法对不同萎凋程度的含水率进行测定；快速水分测定仪设定130℃、称取3.0g样品对不同萎凋程度的含水率进行测定。含水率检测精准度为不同方式检测出的不同萎凋程度叶含水率与真实含水率比值的均值。
以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项发明技术思想的范围内，进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。
Claims (8)
1.一种红茶萎凋程度数字化的无损检测方法，其特征在于，包括以下步骤：
S1、萎凋叶样品光谱的采集，将处于不同萎凋时间节点的萎凋叶样品，分别通过连接智能手机的微型近红外光谱仪采集，进而得到不同萎凋时间节点下的样品光谱；所述萎凋时间节点包括过轻萎凋、轻萎凋、适度萎凋Ⅰ、适度萎凋Ⅱ、重萎凋、过重萎凋；
S2、样品萎凋程度数字化参数的检测，随着萎凋的进行，检测不同萎凋时间节点萎凋叶的萎凋程度数字化参数；
S3、对采集的萎凋叶样品光谱进行预处理和主成分分析，降低采集光谱时萎凋叶样品表面漫反射对光谱信息造成的干扰，提取光谱数据主要信息；
S4、Elman神经网络的建立，将S3进行预处理后的光谱数据融合，通过扩展输入空间的方法，按照5：1随机划分为校正集和预测集，校正集与S2步骤的不同萎凋时间节点下萎凋叶相对应的萎凋程度数字化参数用于训练Elman神经网络的红茶萎凋水分定量预测模型，预测集作为外部验证样本，进行萎凋程度数字化参数的预测；根据预测集样本所预测的萎凋程度数字化参数与实际值的相关系数分析红茶萎凋水分定量预测模型的性能；
S5、不同萎凋程度数字化参数的快速检测，通过无线蓝牙技术连接到智能手机的微型近红外光谱仪，利用智能手机中的应用程序操控微型近红外光谱仪完成萎凋叶样品光谱数据的采集和存储，通过USB将存储的光谱数据传输到PC端，通过MATLAB软件建立水分预测模型；
S6、将PC端建立好的水分预测模型部署到微型计算机Raspberry Pi上独立运行，通过微型近红外光谱仪、OLED显示屏与微型计算机Raspberry Pi通信后，光谱数据即可实时传输至已部署好水分预测算法模型的微型计算机Raspberry Pi上，通过模型运行计算后，萎凋叶水分预测结果即可通过OLED显示屏实时显示。
2.根据权利要求1所述的一种红茶萎凋程度数字化的无损检测方法，其特征在于，所述S1步骤中，从红茶萎凋开始至萎凋过度后1小时止，共18h；期间每隔1h测取该萎凋程度下萎凋叶的光谱数据。
3.根据权利要求1所述的一种红茶萎凋程度数字化的无损检测方法，其特征在于，所述微型近红外光谱仪的近红外波段范围为900-1700nm，分辨率为10-12nm；微型计算机Raspberry Pi为4B主板，OLED显示屏控制芯片为SSD1306，通信方式为IIC。
4.根据权利要求1所述的一种红茶萎凋程度数字化的无损检测方法，其特征在于，所述S1步骤中光谱采集的条件为将每个萎凋程度下萎凋叶样品轻轻按压均匀的表面并扫描3次，然后将平均光谱用作样品的代表性光谱数据。
5.根据权利要求1所述的一种红茶萎凋程度数字化的无损检测方法，其特征在于，所述S2步骤中描述红茶萎凋程度数字化的指标是萎凋叶的含水率，所述含水率采用水分测定仪的加热计算差值测定。
6.根据权利要求1所述的一种红茶萎凋程度数字化的无损检测方法，其特征在于，所述S3步骤中，对近红外光谱采用平滑预处理，对预处理后的光谱进行PCA主成分分析。
7.根据权利要求6所述的一种红茶萎凋程度数字化的无损检测方法，其特征在于，所述PCA主成分分析的主成分数为10。
8.根据权利要求1所述的一种红茶萎凋程度数字化的无损检测方法，其特征在于，所述Elman神经网络结构包括2个输入层、9个隐藏层、1个输出层，它比普通的神经网络在结构上多一个承接层，承接层用于承接整个网络上一次的状态，可看作神经网络的记忆延时单元；承接层的增益因子取[0,1]。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110765625.7A CN113552089A (zh) | 2021-07-06 | 2021-07-06 | 一种红茶萎凋程度数字化的无损检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110765625.7A CN113552089A (zh) | 2021-07-06 | 2021-07-06 | 一种红茶萎凋程度数字化的无损检测方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113552089A true CN113552089A (zh) | 2021-10-26 |
Family
ID=78102792
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110765625.7A Pending CN113552089A (zh) | 2021-07-06 | 2021-07-06 | 一种红茶萎凋程度数字化的无损检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113552089A (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009014700A (ja) * | 2007-01-31 | 2009-01-22 | Osaka Univ | 緑茶の品質予測方法 |
| CN104931453A (zh) * | 2015-06-12 | 2015-09-23 | 湖北省农业科学院果树茶叶研究所 | 一种基于近红外光谱技术预测绿茶摊青叶含水量的方法 |
| CN106290230A (zh) * | 2016-07-26 | 2017-01-04 | 安徽农业大学 | 一种近红外光谱结合化学成分的红茶萎凋程度的判别方法 |
| CN106568738A (zh) * | 2016-10-31 | 2017-04-19 | 湖北省农业科学院果树茶叶研究所 | 一种近红外光谱快速判定不同质量等级茶鲜叶的方法 |
| CN106568739A (zh) * | 2016-10-31 | 2017-04-19 | 湖北省农业科学院果树茶叶研究所 | 一种近红外光谱快速判定不同海拔高度茶鲜叶的方法 |
| US20190162658A1 (en) * | 2017-11-24 | 2019-05-30 | Oil Crops Research Institute, Chinese Acadamy Of Agricultural Sciences | Method for detecting multivariate adulteration of edible oil based on near-infrared spectroscopy |
| CN111524092A (zh) * | 2020-03-16 | 2020-08-11 | 中国农业科学院茶叶研究所 | 一种红茶萎凋程度指标的无损检测方法 |
-
2021
- 2021-07-06 CN CN202110765625.7A patent/CN113552089A/zh active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009014700A (ja) * | 2007-01-31 | 2009-01-22 | Osaka Univ | 緑茶の品質予測方法 |
| CN104931453A (zh) * | 2015-06-12 | 2015-09-23 | 湖北省农业科学院果树茶叶研究所 | 一种基于近红外光谱技术预测绿茶摊青叶含水量的方法 |
| CN106290230A (zh) * | 2016-07-26 | 2017-01-04 | 安徽农业大学 | 一种近红外光谱结合化学成分的红茶萎凋程度的判别方法 |
| CN106568738A (zh) * | 2016-10-31 | 2017-04-19 | 湖北省农业科学院果树茶叶研究所 | 一种近红外光谱快速判定不同质量等级茶鲜叶的方法 |
| CN106568739A (zh) * | 2016-10-31 | 2017-04-19 | 湖北省农业科学院果树茶叶研究所 | 一种近红外光谱快速判定不同海拔高度茶鲜叶的方法 |
| US20190162658A1 (en) * | 2017-11-24 | 2019-05-30 | Oil Crops Research Institute, Chinese Acadamy Of Agricultural Sciences | Method for detecting multivariate adulteration of edible oil based on near-infrared spectroscopy |
| CN111524092A (zh) * | 2020-03-16 | 2020-08-11 | 中国农业科学院茶叶研究所 | 一种红茶萎凋程度指标的无损检测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cozzolino | Use of infrared spectroscopy for in-field measurement and phenotyping of plant properties: instrumentation, data analysis, and examples | |
| CN102879353B (zh) | 近红外检测花生中蛋白质组分含量的方法 | |
| CN108152235B (zh) | 一种联合土壤室内外光谱的重金属含量反演方法 | |
| CN109253985B (zh) | 基于神经网络的近红外光谱识别古筝面板用木材等级的方法 | |
| CN112098358B (zh) | 基于四元数卷积神经网络的近红外光谱并行融合定量检测方法 | |
| CN110118742A (zh) | 基于光谱分析的陆表植被冠层生态水含量遥感反演方法 | |
| CN106841083A (zh) | 基于近红外光谱技术的芝麻油品质检测方法 | |
| CN110672546A (zh) | 一种基于便携式近红外光谱仪的酒糟模型建模方法 | |
| CN108169165B (zh) | 基于太赫兹光谱和图像信息融合的麦芽糖混合物定量分析方法 | |
| CN102876816A (zh) | 基于多传感器信息融合的发酵过程状态监测与控制方法 | |
| Power et al. | From the laboratory to the vineyard—evolution of the measurement of grape composition using NIR spectroscopy towards high-throughput analysis | |
| Ning et al. | Rapid evaluation of soil fertility in tea plantation based on near-infrared spectroscopy | |
| Shen et al. | Rapid and real-time detection of moisture in black tea during withering using micro-near-infrared spectroscopy | |
| CN114818985B (zh) | 基于中心锚点三元组优化伪孪生网络的茶叶品质评价方法 | |
| CN114676636A (zh) | 一种综合植被和生境特征的草原区土壤水分快速反演方法 | |
| Puttipipatkajorn et al. | Development of calibration models for rapid determination of moisture content in rubber sheets using portable near-infrared spectrometers | |
| CN108693139A (zh) | 电子烟油理化指标的近红外光谱预测模型建立方法及应用 | |
| CN110567894A (zh) | 一种含水溶液pH值快速检测方法和检测装置 | |
| CN103234923A (zh) | 一种黄酒发酵过程中总糖含量的在线监控方法 | |
| Gong et al. | Spectral reflectance reconstruction from red-green-blue (RGB) images for chlorophyll content detection | |
| CN108007917B (zh) | 希尔伯特法建立稻株中氮素含量拉曼光谱测量模型方法 | |
| CN113552089A (zh) | 一种红茶萎凋程度数字化的无损检测方法 | |
| CN108007916B (zh) | 希尔伯特黄法建立稻株氮含量的共聚焦显微拉曼测量模型 | |
| Hu et al. | Nondestructive discrimination of internal defects in jujube (Huizao) of Xinjiang based on visible and near-infrared spectroscopy | |
| CN111220565B (zh) | 一种基于cpls的红外光谱测量仪器标定迁移方法 |
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: 20211026 |
|
| RJ01 | Rejection of invention patent application after publication |