CN107860740A - 一种基于近红外光谱技术的红茶发酵品质的评价方法 - Google Patents
一种基于近红外光谱技术的红茶发酵品质的评价方法 Download PDFInfo
- Publication number
- CN107860740A CN107860740A CN201711292162.7A CN201711292162A CN107860740A CN 107860740 A CN107860740 A CN 107860740A CN 201711292162 A CN201711292162 A CN 201711292162A CN 107860740 A CN107860740 A CN 107860740A
- Authority
- CN
- China
- Prior art keywords
- fermentation
- quality
- sample
- black tea
- model
- 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
- 238000000855 fermentation Methods 0.000 title claims abstract description 79
- 230000004151 fermentation Effects 0.000 title claims abstract description 79
- 244000269722 Thea sinensis Species 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 42
- 235000006468 Thea sinensis Nutrition 0.000 title claims abstract description 35
- 235000020279 black tea Nutrition 0.000 title claims abstract description 35
- 238000011156 evaluation Methods 0.000 title claims abstract description 34
- 238000002329 infrared spectrum Methods 0.000 title claims abstract description 20
- 230000003595 spectral effect Effects 0.000 claims abstract description 24
- 235000013616 tea Nutrition 0.000 claims abstract description 22
- 238000005516 engineering process Methods 0.000 claims abstract description 10
- 238000002360 preparation method Methods 0.000 claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 14
- 239000000126 substance Substances 0.000 claims description 14
- ADRVNXBAWSRFAJ-UHFFFAOYSA-N catechin Natural products OC1Cc2cc(O)cc(O)c2OC1c3ccc(O)c(O)c3 ADRVNXBAWSRFAJ-UHFFFAOYSA-N 0.000 claims description 13
- 235000005487 catechin Nutrition 0.000 claims description 13
- 238000012216 screening Methods 0.000 claims description 13
- PFTAWBLQPZVEMU-DZGCQCFKSA-N (+)-catechin Chemical compound C1([C@H]2OC3=CC(O)=CC(O)=C3C[C@@H]2O)=CC=C(O)C(O)=C1 PFTAWBLQPZVEMU-DZGCQCFKSA-N 0.000 claims description 12
- 229950001002 cianidanol Drugs 0.000 claims description 12
- IPMYMEWFZKHGAX-UHFFFAOYSA-N Isotheaflavin Natural products OC1CC2=C(O)C=C(O)C=C2OC1C(C1=C2)=CC(O)=C(O)C1=C(O)C(=O)C=C2C1C(O)CC2=C(O)C=C(O)C=C2O1 IPMYMEWFZKHGAX-UHFFFAOYSA-N 0.000 claims description 11
- UXRMWRBWCAGDQB-UHFFFAOYSA-N Theaflavin Natural products C1=CC(C2C(CC3=C(O)C=C(O)C=C3O2)O)=C(O)C(=O)C2=C1C(C1OC3=CC(O)=CC(O)=C3CC1O)=CC(O)=C2O UXRMWRBWCAGDQB-UHFFFAOYSA-N 0.000 claims description 11
- 230000001953 sensory effect Effects 0.000 claims description 11
- IPMYMEWFZKHGAX-ZKSIBHASSA-N theaflavin Chemical compound C1=C2C([C@H]3OC4=CC(O)=CC(O)=C4C[C@H]3O)=CC(O)=C(O)C2=C(O)C(=O)C=C1[C@@H]1[C@H](O)CC2=C(O)C=C(O)C=C2O1 IPMYMEWFZKHGAX-ZKSIBHASSA-N 0.000 claims description 11
- 235000014620 theaflavin Nutrition 0.000 claims description 11
- 229940026509 theaflavin Drugs 0.000 claims description 11
- 235000008118 thearubigins Nutrition 0.000 claims description 9
- 239000000203 mixture Substances 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 7
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 claims description 6
- XABJJJZIQNZSIM-UHFFFAOYSA-N azane;phenol Chemical compound [NH4+].[O-]C1=CC=CC=C1 XABJJJZIQNZSIM-UHFFFAOYSA-N 0.000 claims description 6
- 238000001035 drying Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 238000004445 quantitative analysis Methods 0.000 claims description 4
- ISWSIDIOOBJBQZ-UHFFFAOYSA-N Phenol Chemical compound OC1=CC=CC=C1 ISWSIDIOOBJBQZ-UHFFFAOYSA-N 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 229910021529 ammonia Inorganic materials 0.000 claims description 3
- 229920002678 cellulose Polymers 0.000 claims description 3
- 239000001913 cellulose Substances 0.000 claims description 3
- 230000002860 competitive effect Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000010187 selection method Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 claims description 3
- 238000004108 freeze drying Methods 0.000 claims description 2
- 239000001054 red pigment Substances 0.000 claims description 2
- 238000013515 script Methods 0.000 claims description 2
- 238000004611 spectroscopical analysis Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000013441 quality evaluation Methods 0.000 abstract description 3
- 238000009614 chemical analysis method Methods 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 44
- 230000008859 change Effects 0.000 description 7
- DATAGRPVKZEWHA-YFKPBYRVSA-N N(5)-ethyl-L-glutamine Chemical compound CCNC(=O)CC[C@H]([NH3+])C([O-])=O DATAGRPVKZEWHA-YFKPBYRVSA-N 0.000 description 6
- 230000002596 correlated effect Effects 0.000 description 5
- 230000000875 corresponding effect Effects 0.000 description 5
- 235000013824 polyphenols Nutrition 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 150000008442 polyphenolic compounds Chemical class 0.000 description 4
- 238000012545 processing Methods 0.000 description 3
- 229940026510 theanine Drugs 0.000 description 3
- WMBWREPUVVBILR-UHFFFAOYSA-N GCG Natural products C=1C(O)=C(O)C(O)=CC=1C1OC2=CC(O)=CC(O)=C2CC1OC(=O)C1=CC(O)=C(O)C(O)=C1 WMBWREPUVVBILR-UHFFFAOYSA-N 0.000 description 2
- 244000062793 Sorghum vulgare Species 0.000 description 2
- 238000002835 absorbance Methods 0.000 description 2
- 238000007705 chemical test Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 235000019713 millet Nutrition 0.000 description 2
- 238000007254 oxidation reaction Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 238000003892 spreading Methods 0.000 description 2
- WMBWREPUVVBILR-WIYYLYMNSA-N (-)-Epigallocatechin-3-o-gallate Chemical compound O([C@@H]1CC2=C(O)C=C(C=C2O[C@@H]1C=1C=C(O)C(O)=C(O)C=1)O)C(=O)C1=CC(O)=C(O)C(O)=C1 WMBWREPUVVBILR-WIYYLYMNSA-N 0.000 description 1
- 235000009024 Ceanothus sanguineus Nutrition 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 240000003553 Leptospermum scoparium Species 0.000 description 1
- 235000015459 Lycium barbarum Nutrition 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 150000001413 amino acids Chemical class 0.000 description 1
- 239000003963 antioxidant agent Substances 0.000 description 1
- 230000003078 antioxidant effect Effects 0.000 description 1
- 235000006708 antioxidants Nutrition 0.000 description 1
- 150000001765 catechin Chemical class 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 235000009508 confectionery Nutrition 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002255 enzymatic effect Effects 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 239000000796 flavoring agent Substances 0.000 description 1
- 235000019634 flavors Nutrition 0.000 description 1
- 238000004898 kneading Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004811 liquid chromatography Methods 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 239000002932 luster Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 150000002894 organic compounds Chemical class 0.000 description 1
- 235000020567 organic tea Nutrition 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 239000012071 phase Substances 0.000 description 1
- 239000000049 pigment Substances 0.000 description 1
- -1 polyphenol compound Chemical class 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
-
- 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
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
本发明涉及检测技术领域,尤其涉及一种基于近红外光谱技术的红茶发酵品质的评价方法,包括以下步骤:红茶样本的制备;光谱信息的采集;发酵品质指标的测定,样品集的划分,光谱信息的预处理,特征波长变量筛选,发酵品质预测模型建立和各发酵品质指标预测。本发明样品预处理简单,能准确、快速地鉴定出发酵的品质优劣,该方法可操作性强,避免了人工感官审评中主观意识的影响以及传统化学分析法的繁杂操作,结果客观、准确,提高红茶发酵品质评价的客观性和精准性,为进一步实现自动化、智能化和精准化发酵控制技术提供理论基础,在茶叶品质指标检测中表现出很好的现实意义和应用前景。
Description
技术领域
本发明属于检测技术领域,尤其是涉及一种基于近红外光谱技术的红茶发酵品质的评价方法。
背景技术
红茶是世界上消费量最大的茶类,具有清除体内自由基和抗氧化等保健作用。发酵是儿茶素为主体的多酚类化合物发生酶促氧化反应,生成水溶性的有色氧化产物茶黄素(TFs)、茶红素(TRs)、茶褐素(TBs)和芬香物质,最终形成了红茶特有的色、香、味等品质的化学过程。
大量研究表明,决定红茶发酵品质的主要生化指标有:儿茶素组分(Catechinscontent,C、EC、GC、EGC、CG、ECG、GCG、EGCG等物质的总量)、酚氨比(TP)/游离氨基酸总量(FAA)、TFs/TRs。茶多酚呈“苦”味,而氛基酸呈“甜”味,合适的酚氨比是形成茶汤的“鲜”味特征。儿茶素是TFs和TRs形成的基础,其含量与红茶感官品质呈正相关。茶黄素(TFs)和茶红素(TRs)是决定茶汤亮度和色泽的化学物质,并影响味觉强度、浓度和鲜爽度。此外,发酵中叶面色泽发生显著变化(由绿色到红色,再到褐色),本质是TFs和TRs的绝对含量及组成比例变化的表象,当TFs/TRs为1:9~1:10时为最佳发酵适度点,此时在制品呈现红铜色。国内外红茶生产中,均通过观察这种颜色表象变化,凭借经验判断发酵品质适度,制茶学称之为“看茶做茶”,其具有主观性和不稳定性。
若要精准获取发酵品质信息,必须借助化学仪器持续地检测品质成分,如:分光光度计法是茶色素成分的通用分析方法,样品前处理及检测过程需要耗费大量的时间,且需要具备专业技能的实验人员操作,显然传统的化学分析手段无法在线、实时和快速的获取品质成分信息,此已成为自动化和智能化红茶加工技术实现的主要制约因素。
茶多酚(化学式为C22H18O11)、儿茶素(C15H14O6)、茶氨酸(C7H14N2O3)、茶黄素(C29H24O12)和茶红素(C33H47NO13)以有机化合物形式存在,分子结构复杂,其含量随发酵进程不断变化,具有复杂的光谱吸收特征,这些光谱信息可用于分析和表征内质组分的变化。
因此,探求一种加工过程中发酵品质指标的快速检测方法,对实现加工品质的实时监控技术至关重要。
发明内容
本发明的目的在于,旨在克服人工感官评审和理化检测的不足,提高红茶发酵品质评价的客观性和精准性,为进一步实现自动化、智能化和精准化发酵控制技术提供理论基础,而提供一种基于近红外光谱技术的红茶发酵品质的评价方法。
为此,本发明的上述目的通过以下技术方案来实现:
一种基于近红外光谱技术的红茶发酵品质的评价方法,其特征在于,所述基于近红外光谱技术的红茶发酵品质的评价方法包括以下操作:
(1)红茶样本的制备:样本为发酵过程中不同时序下的条形工夫红茶发酵叶,具体为:发酵总时长为360min,每隔5min取样1次,每次取在制品样本5个;
每次取样的5个在制品样本中:3个用于光谱信息的采集;1个样品经冷冻干燥至足干,用于检测理化成分;1个样品经电磁干燥至足干,用于检测感官评分;
(2)光谱信息的采集:采用光谱分析仪对发酵样本进行漫反射光谱扫描获取的近红外光谱数据;
(3)发酵品质的理化指标测定:测定的每个冷冻干燥样品的酚氨比、儿茶素、茶黄素、茶红素和茶褐素含量;发酵品质的感官评分:测定每个电磁干燥样品的感官评分;
(4)样品集的划分:采用K-S法对样品集进行校正集和预测集的划分,校正集和测试集比例约为3:1;
(5)光谱信息的预处理:对获取的近红外光谱数据在Matlab软件中调用Zscore脚本进行预处理,并通过主成分分析法对预处理后的近红外光谱数据进行降维,获取主成分数据;
(6)特征波长变量筛选:采用联合区间算法、竞争性自适应权重取样算法、随机蛙跳算法和连续投影算法对光谱特征波长变量进行筛选,并基于偏最小二乘法进行建模;然后分别筛选与感官评分、茶黄素、茶红素、茶褐素、儿茶素和酚氨比密切相关的特征波长变量;每样物质在不同波长下吸光度会不同(造成响应光谱的形态是波浪曲线),随着发酵的变化,成分含量高低和分子组成随之变化,会引起吸光度的变化,变量筛选就是要找到与成分相关的有响应关系的特征波长。
对每种方法的筛选结果建立线性模型进行比较,选择预测集中RPD值最高的方法,作为该品质指标的最优变量筛选方法;
(7)发酵品质预测模型建立和各发酵品质指标预测:基于筛选的特征波长变量,采用支持向量机回归法,分别建立各发酵品质指标的非线性定量预测模型,发酵品质预测模型为包含各发酵品质指标的非线性定量预测模型。
在采用上述技术方案的同时,本发明还可以采用或者组合采用以下进一步的技术方案:
优选地,(2)中:光谱分析仪采用SupNIR-1520光谱分析仪,波长范围为1000~1800nm;光谱精度为±1nm。
优选地,(2)中:进行光谱扫描时,每个样本取约50g于样品池中,每个样品进行3次光谱采集,单次扫描32次,取3次采集光谱的均值作为该样品的光谱。
优选地,(7)中:在设计基于支持向量机回归的非线性定量模型时,使用MATLAB软件下的Li-SVRLAB工具箱创建GA-SVR,选用RBF核函数,训练建模函数为trainlsvr,测试函数为simlssvr,以10组主成分数据作为基于支持向量机回归的非线性定量模型的输入,以其对应的各品质成分含量作为网络的输出,采用tunelsssvm函数来优选参数。
优选地,(7)中:发酵品质预测模型评价指标包括:校正集的相关系数Rc、RMSEC、预测集的相关系数Rp和预测均方根误差RMSEP;通常Rc、Rp值越大,RMSEC、RMSEP值越小,则所建发酵品质预测模型的预测性能、泛化性能越好,并将相对标准偏差RPD作为最终评价指标,当RPD>2时,表明发酵品质预测模型具有极好的预测能力;当RPD在1.8-2.0时,表明发酵品质预测模型效果较好,可用于样品的定量分析;当RPD在1.4-1.8时,表明发酵品质预测模型可对样品作粗略预测和关联评估;而RPD<1.4时,表明发酵品质预测模型效果很差,无法应用。
本发明提供一种基于近红外光谱技术的红茶发酵品质的评价方法,具有如下优点:
(1)本发明的红茶发酵品质评价方法基于近红外光谱技术结合变量筛选方法,建立非线性模型可实现对发酵品质指标的量化评价,从而克服人工感官评审和理化检测的不足,提高红茶发酵品质评价的客观性和精准性;
(2)本发明红茶发酵品质评价方法,在模型建立后,可实现红茶加工过程中对发酵品质指标含量的快速在线检测,突破传统取样后进行大量预处理(干燥、固样、理化分析)后,经专门仪器分析或感官评审判断发酵品质程度的缺陷,提高了品质评价的时效性、无损伤性、在线和快速性。
(3)本发明的红茶发酵品质评价方法为进一步实现自动化、智能化和精准化发酵控制技术提供理论基础,在茶叶检测中有较好的应用前景。
附图说明
图1为光谱信息的预处理中通过主成分分析法得到的主成分的分散点图;
图2为基于不同筛选算法建模得到的感官品质的特征波长变量优选图;
图3a为基于连续投影算法建模对感官品质的特征波长变量个数的选择过程图;
图3b为基于连续投影算法建模筛选得到的感官品质的特征波长变量与光谱位置的对应图;
图4为各理化指标的最优特征波长变量结果图。
具体实施方式
参照附图和具体实施例对本发明作进一步详细地描述。
一种基于近红外光谱技术的红茶发酵品质的评价方法,包括以下步骤:
(1)红茶样本的制备
试验于2015年7月在中国最大的有机茶叶加工基地(浙江更香茶厂)实施,茶鲜叶原料采用福鼎大白品种(中国制作红茶的主要茶树品种),嫩度为一芽一叶,萎凋、揉捻、发酵和干燥等工艺条件一致。试验工艺:
1)萎凋:采用自然摊放萎凋方式,匀叶厚度为1-2cm,温度控制在30℃左右,湿度控制在60%左右,萎凋至水分含量为60%左右。
2)揉捻:采用40型揉捻机,按照空揉、轻压、中压和重压各交替方式进行揉捻。
3)发酵:设定发酵温度为27℃、发酵湿度大于90%和翻拌间隔为25min。
4)干燥:干燥分为2次,首先120℃毛火干燥,烘至水分为20%左右;风冷后,进行90℃足火干燥,烘至水分为5%左右。
红茶样本的制备,样本为发酵过程中不同时序的条形工夫红茶发酵叶;具体为发酵过程为360min,每隔5min取样1次,每次取在制品样本5个;
每次取样的5个在制品样本中,3个用于光谱信息的采集;1个样品经冷冻干燥至足干,用于检测理化成分;1个样品经电磁干燥至足干,用于检测感官评分。
(2)光谱信息的采集
采用SupNIR-1520光谱分析仪对发酵样本进行漫反射光谱扫描。波长范围:1000~1800nm;光谱精度:±1nm。每个样本取约50g于样品池中,每个样品进行3次光谱采集(单次扫描32次),取3次采集光谱的均值作为该样品的光谱。
(3)发酵品质指标的测定
样品的理化成分:按照国标法《GB/T 30483-2013茶叶中茶黄素的测定-高效液相色谱法》测定茶色素含量;按照GB/T 31740.2-2015茶制品第2部分:茶多酚测定茶多酚和儿茶素含量;按照GB/T 23193-2008茶叶中茶氨酸的测定高效液相色谱法测定茶氨酸含量;对样品进行冻干和磨粉预处理,利用高效液相色谱仪测定茶黄素、茶红素和茶褐素等理化品质成分含量;
感官品质评分:依据GB/T 23776-2009《茶叶感官审评方法》及密码审评形式,评定各发酵时序下的成品红茶的感官品质,并计算审评总分。
(4)样品集的划分
采用K-S法对样品集进行校正集和预测集的划分,从收集的219个样本集中挑选出146个作为校正集,余下的55个作为预测集,校正集和测试集比例约为3:1,表1为样本集中各发酵品质指标值的分布情况。
表1
aS.N.,样本数.bS.D.,标准偏差.
(5)光谱信息的预处理
对获取的近红外光谱数据在Matlab软件中调用Zscore脚本进行预处理;
进一步,通过主成分分析法(PCA)对预处理后的光谱数据进行降维,获取主成分数据。第一主成分(PC1)得分为56.44%,是发酵样品的NIR光谱的最重要信息,第二主成分(PC2)得分为37.35%,前两个主成分累积贡献率为93.79%,仅6.21%的信息丢失,表明2个主成分能代表原始光谱的绝大部分信息。图1为前二个主成分的得分散点图,有图1可知,同一发酵程度样品的得分分布能够成簇聚集,不同发酵阶段样品可以聚集在不同的区域,反映出不同发酵适度的样品具有差异性,NIR光谱结合主成分分析法能较有效区分不同发酵品质的红茶样品。
(6)特征波长变量筛选
采用联合区间算法(Si)、竞争性自适应权重取样算法(CARS)、随机蛙跳算法(SFLA)和连续投影算法(SPA)对光谱特征波长变量进行筛选,并基于偏最小二乘法(PLS)进行建模,分别得到Si-PLS模型、CARS-PLS模型、SFLA-PLS模型和SPA-PLS模型;然后分别筛选与感官评分、茶黄素、茶红素、茶褐素、儿茶素和酚氨比密切相关的特征波长变量;
对每种方法的筛选结果建立线性模型进行比较,并绘制图2,图2为基于不同筛选算法建模得到的感官品质的特征波长变量优选图。表2为不同变量优选方法所建感官品质评价模型的结果,如表2所示,选择预测集中RPD值最高的方法,作为该品质指标的最优变量筛选方法。
图3a为基于连续投影算法建模对感官品质的特征波长变量个数的选择过程图。感官品质的最优波长筛选模型为SPA-PLS模型,按照被选择的概率降序排名选出9个特征波长变量,分别为:1561、1076、1142、1567、1429、1532、1792、1646和1617nm;并绘制得到图3b。图3b为基于连续投影算法建模筛选得到的感官品质的特征波长变量与光谱位置的对应图;
各理化品质特征波长的优选:
茶黄素相关特征波长变量的筛选模型为SPA-PLS模型,其以投影向量最大的波长为待选波长,并按向量值的大小降序建立由待选波长组成的集合,将集合拆分为多个子集,子集数量分别为1,2…,波长数依次增加,构建各子集的校正模型,以模型的RMSEC最小值对应子集为最优组合,组合内波长为特征波长,选出16个敏感特征波长;
茶红素相关特征波长变量的筛选模型为SPA-PLS模型,按照被选择的概率降序排名筛选出10个敏感特征波长;茶褐素相关特征波长变量的筛选模型为CARS-PLS模型,其随机抽取一定数量的波长构建校正模型,以回归系数绝对值最大的点作为待选波长,并按照系数绝对值大小降序,经过50-100次的循环提出后,建立由待选波长组成的集合,将集合拆分为多个子集,子集数量分别为1,2…,波长数依次增加,构建各子集的校正模型,以模型的RMSEC最小值对应子集为最优组合,组合内波长为特征波长。按照被选择的概率降序排名筛选出33个敏感特征波长;儿茶素总量相关特征波长变量的筛选模型为SPA-PLS模型,按照被选择的概率降序排名筛选出9个极其敏感的特征波长;酚氨比相关特征波长变量的筛选模型为SFLA-PLS模型,按照被选择的概率降序排名筛选出12个极其敏感的特征波长,并绘制得到图4,图4为各理化指标的最优特征波长变量结果图。
表2
(7)模型建立和各发酵品质指标预测
在设计SVR(支持向量机回归法)非线性定量模型时,使用MATLAB软件下的Li-SVRLAB工具箱创建GA-SVR,选用RBF核函数,训练建模函数为trainlsvr,测试函数为simlssvr,以10组主成分数据作为SVR模型的输入,以其对应的各品质成分含量作为网络的输出,采用tunelsssvm函数来优选参数,所建GA-SVR模型结果见表3。表3为各发酵品质指标的GA-SVR模型输出结果。
表3
从表3可知,感官品质、茶褐素和儿茶素的最优变量模型分RPD值均大于3,分别为3.923、3.234和5.462,表明近红外定量分析结果准确、可靠,模型预测性能极高,可用于质量控制。酚氨比的RPD值为2.815(大于2),说明模型具有很好的预测效果,可用于定量分析。茶红素的RPD值为2.223(大于2),表明该模型具有很好的预测性能。茶黄素的RPD值为1.770(介于1.4-1.8之间),表明模型预测能力较差,可对样品作粗略预测和关联评估。
上述具体实施方式用来解释说明本发明,仅为本发明的优选实施例,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改、等同替换、改进等,都落入本发明的保护范围。
Claims (5)
1.一种基于近红外光谱技术的红茶发酵品质的评价方法,其特征在于,所述基于近红外光谱技术的红茶发酵品质的评价方法包括以下操作:
(1)红茶样本的制备:样本为发酵过程中不同时序下的条形工夫红茶发酵叶,具体为:发酵总时长为360 min,每隔5 min取样1次,每次取在制品样本5个;
每次取样的5个在制品样本中:3个用于光谱信息的采集;1个样品经冷冻干燥至足干,用于检测理化成分;1个样品经电磁干燥至足干,用于检测感官评分;
(2)光谱信息的采集:采用光谱分析仪对发酵样本进行漫反射光谱扫描获取的近红外光谱数据;
(3)发酵品质的理化指标测定:测定的每个冷冻干燥样品的酚氨比、儿茶素、茶黄素、茶红素和茶褐素含量;发酵品质的感官评分:测定每个电磁干燥样品的感官评分;
(4)样品集的划分:采用K-S法对样品集进行校正集和预测集的划分,校正集和测试集比例约为3:1;
(5)光谱信息的预处理:对获取的近红外光谱数据在Matlab软件中调用Zscore脚本进行预处理,并通过主成分分析法对预处理后的近红外光谱数据进行降维,获取主成分数据;
(6)特征波长变量筛选:采用联合区间算法、竞争性自适应权重取样算法、随机蛙跳算法和连续投影算法对光谱特征波长变量进行筛选,并基于偏最小二乘法进行建模;然后分别筛选能够预测感官评分、茶黄素、茶红素、茶褐素、儿茶素和酚氨比含量的特征波长变量;
对每种方法的筛选结果建立线性模型进行比较,选择预测集中RPD值最高的方法,作为该品质指标的最优变量筛选方法;
(7)发酵品质预测模型建立和各发酵品质指标预测:基于筛选的特征波长变量,采用支持向量机回归法,分别建立各发酵品质指标的非线性定量预测模型,发酵品质预测模型为包含各发酵品质指标的非线性定量预测模型。
2. 根据权利要求1所述的基于近红外光谱技术的红茶发酵品质的评价方法,其特征在于,(2)中:光谱分析仪采用SupNIR-1520光谱分析仪,波长范围为1000~1800 nm;光谱精度为±1 nm。
3. 根据权利要求1或2所述的基于近红外光谱技术的红茶发酵品质的评价方法,其特征在于,(2)中:进行光谱扫描时,每个样本取约50 g于样品池中,每个样品进行3次光谱采集,单次扫描32次,取3次采集光谱的均值作为该样品的光谱。
4.根据权利要求1所述的基于近红外光谱技术的红茶发酵品质的评价方法,其特征在于,(7)中:在设计基于支持向量机回归的非线性定量模型时,使用MATLAB软件下的Li-SVRLAB工具箱创建GA-SVR,选用RBF核函数,训练建模函数为trainlsvr,测试函数为simlssvr,以10组主成分数据作为基于支持向量机回归的非线性定量模型的输入,以其对应的各品质成分含量作为网络的输出,采用tunelsssvm函数来优选参数。
5. 根据权利要求1所述的基于近红外光谱技术的红茶发酵品质的评价方法,其特征在于,(7)中:发酵品质预测模型评价指标包括:校正集的相关系数Rc、RMSEC、预测集的相关系数Rp和预测均方根误差RMSEP;通常Rc、Rp值越大,RMSEC、RMSEP值越小,则所建发酵品质预测模型的预测性能、泛化性能越好,并将相对标准偏差RPD作为最终评价指标,当RPD > 2时,表明发酵品质预测模型具有极好的预测能力;当RPD在1.8-2.0时,表明发酵品质预测模型效果较好,可用于样品的定量分析;当RPD在1.4-1.8时,表明发酵品质预测模型可对样品作粗略预测和关联评估;而RPD < 1.4时,表明发酵品质预测模型效果很差,无法应用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711292162.7A CN107860740A (zh) | 2017-12-08 | 2017-12-08 | 一种基于近红外光谱技术的红茶发酵品质的评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711292162.7A CN107860740A (zh) | 2017-12-08 | 2017-12-08 | 一种基于近红外光谱技术的红茶发酵品质的评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107860740A true CN107860740A (zh) | 2018-03-30 |
Family
ID=61705208
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711292162.7A Pending CN107860740A (zh) | 2017-12-08 | 2017-12-08 | 一种基于近红外光谱技术的红茶发酵品质的评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107860740A (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109030408A (zh) * | 2018-06-29 | 2018-12-18 | 无锡济民可信山禾药业股份有限公司 | 一种快速的儿茶药材检测方法 |
CN109329501A (zh) * | 2018-10-29 | 2019-02-15 | 江苏大学 | 一种速溶发酵茶粉的智能化制备方法 |
CN109497194A (zh) * | 2018-10-29 | 2019-03-22 | 江苏大学 | 一种发酵茶智能化固态发酵方法 |
CN109668858A (zh) * | 2019-02-14 | 2019-04-23 | 大连理工大学 | 基于近红外光谱检测发酵过程生物量和组分浓度的方法 |
CN110308111A (zh) * | 2019-06-14 | 2019-10-08 | 湖北省农业科学院果树茶叶研究所 | 一种应用近红外光谱技术快速预测远安黄茶闷黄时间的方法 |
CN110308110A (zh) * | 2019-06-14 | 2019-10-08 | 湖北省农业科学院果树茶叶研究所 | 基于最小二乘支持向量机的远安黄茶闷黄时间的无损预测方法 |
CN110530914A (zh) * | 2019-09-26 | 2019-12-03 | 北京农业质量标准与检测技术研究中心 | 土壤重金属检测系统及检测方法 |
CN111044483A (zh) * | 2019-12-27 | 2020-04-21 | 武汉工程大学 | 一种基于近红外光谱的奶油中色素测定方法、系统和介质 |
CN111458308A (zh) * | 2020-06-04 | 2020-07-28 | 杭州岚达科技有限公司 | 一种近红外光谱龙胆鉴别方法和系统 |
CN111480712A (zh) * | 2020-05-26 | 2020-08-04 | 无锡迅杰光远科技有限公司 | 一种茶叶发酵过程全自动品质控制设备 |
CN112085342A (zh) * | 2020-08-14 | 2020-12-15 | 北京农业质量标准与检测技术研究中心 | 农业投入品复合品质评价方法 |
CN112435721A (zh) * | 2020-11-03 | 2021-03-02 | 江南大学 | 一种基于偏最小二乘构建龙井绿茶品质判别模型的方法 |
CN113092406A (zh) * | 2021-04-08 | 2021-07-09 | 晨光生物科技集团股份有限公司 | 一种快速预判辣椒红外观品质变化的方法 |
CN113360836A (zh) * | 2021-06-21 | 2021-09-07 | 湖北民族大学 | 基于回归方程基础上的藤茶品质评审方法 |
CN113504198A (zh) * | 2021-06-15 | 2021-10-15 | 武汉工程大学 | 一种苋菜红色素浓度检测方法、装置及存储介质 |
CN113567390A (zh) * | 2021-06-24 | 2021-10-29 | 贵州省茶叶研究所 | 一种基于近红外光谱技术审评红茶的方法及系统 |
CN113945539A (zh) * | 2021-10-18 | 2022-01-18 | 江西农业大学 | 基于gwas筛选后近红外光谱特征波段预测品质的方法及系统 |
CN117541763A (zh) * | 2024-01-05 | 2024-02-09 | 汉中市黑金茶科技有限公司 | 一种基于图像识别的茶叶发酵度测定方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106172960A (zh) * | 2016-08-24 | 2016-12-07 | 贵州贵茶有限公司 | 一种红茶的发酵工艺 |
CN106212753A (zh) * | 2016-08-22 | 2016-12-14 | 操长安 | 一种红茶的加工工艺 |
CN106212768A (zh) * | 2016-07-08 | 2016-12-14 | 湖南利源隆茶业有限责任公司 | 一种山苍子茶及其制备方法 |
-
2017
- 2017-12-08 CN CN201711292162.7A patent/CN107860740A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106212768A (zh) * | 2016-07-08 | 2016-12-14 | 湖南利源隆茶业有限责任公司 | 一种山苍子茶及其制备方法 |
CN106212753A (zh) * | 2016-08-22 | 2016-12-14 | 操长安 | 一种红茶的加工工艺 |
CN106172960A (zh) * | 2016-08-24 | 2016-12-07 | 贵州贵茶有限公司 | 一种红茶的发酵工艺 |
Non-Patent Citations (9)
Title |
---|
QUANSHENG CHEN等: "Recent developments of green analytical techniques in analysis of tea’s quality and nutrition", 《TRENDS IN FOOD SCIENCE & TECHNOLOGY》 * |
侯英等: "随机蛙跳算法应用于烟草中多酚含量近红外校正模型的建立", 《中国第五届近红外光谱学术会议》 * |
刘洪林: "基于近红外光谱技术(NIRS)对工夫红茶茶色素含量测定", 《分析检测》 * |
吕海鹏: "普洱茶品质的仪器鉴定研究", 《中国博士学位论文全文数据库 农业科技辑》 * |
宁井铭等: "祁门红茶加工中氨基酸和儿茶素快速检测模型建立", 《光学与光谱分析》 * |
张裕莹: "基于光谱方法检测绿茶中儿茶素单体含量的研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑 2017年期》 * |
彭清维等: "基于可见/近红外光谱技术的湄潭翠芽等级判别", 《茶叶科学》 * |
杨杰: "葡萄内部品质的高光谱成像检测研究", 《中国优秀硕士学位论文全文数据库 农业科技辑》 * |
韦玲冬: "基于近红外光谱技术的红茶品质与等级评判方法的研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑 2017年期》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109030408A (zh) * | 2018-06-29 | 2018-12-18 | 无锡济民可信山禾药业股份有限公司 | 一种快速的儿茶药材检测方法 |
CN109329501A (zh) * | 2018-10-29 | 2019-02-15 | 江苏大学 | 一种速溶发酵茶粉的智能化制备方法 |
CN109497194A (zh) * | 2018-10-29 | 2019-03-22 | 江苏大学 | 一种发酵茶智能化固态发酵方法 |
CN109668858A (zh) * | 2019-02-14 | 2019-04-23 | 大连理工大学 | 基于近红外光谱检测发酵过程生物量和组分浓度的方法 |
CN110308111A (zh) * | 2019-06-14 | 2019-10-08 | 湖北省农业科学院果树茶叶研究所 | 一种应用近红外光谱技术快速预测远安黄茶闷黄时间的方法 |
CN110308110A (zh) * | 2019-06-14 | 2019-10-08 | 湖北省农业科学院果树茶叶研究所 | 基于最小二乘支持向量机的远安黄茶闷黄时间的无损预测方法 |
CN110530914A (zh) * | 2019-09-26 | 2019-12-03 | 北京农业质量标准与检测技术研究中心 | 土壤重金属检测系统及检测方法 |
CN111044483A (zh) * | 2019-12-27 | 2020-04-21 | 武汉工程大学 | 一种基于近红外光谱的奶油中色素测定方法、系统和介质 |
CN111480712A (zh) * | 2020-05-26 | 2020-08-04 | 无锡迅杰光远科技有限公司 | 一种茶叶发酵过程全自动品质控制设备 |
CN111458308A (zh) * | 2020-06-04 | 2020-07-28 | 杭州岚达科技有限公司 | 一种近红外光谱龙胆鉴别方法和系统 |
CN112085342A (zh) * | 2020-08-14 | 2020-12-15 | 北京农业质量标准与检测技术研究中心 | 农业投入品复合品质评价方法 |
CN112085342B (zh) * | 2020-08-14 | 2023-09-01 | 北京市农林科学院 | 农业投入品复合品质评价方法 |
CN112435721A (zh) * | 2020-11-03 | 2021-03-02 | 江南大学 | 一种基于偏最小二乘构建龙井绿茶品质判别模型的方法 |
CN113092406A (zh) * | 2021-04-08 | 2021-07-09 | 晨光生物科技集团股份有限公司 | 一种快速预判辣椒红外观品质变化的方法 |
CN113092406B (zh) * | 2021-04-08 | 2022-08-02 | 晨光生物科技集团股份有限公司 | 一种快速预判辣椒红外观品质变化的方法 |
CN113504198A (zh) * | 2021-06-15 | 2021-10-15 | 武汉工程大学 | 一种苋菜红色素浓度检测方法、装置及存储介质 |
CN113360836A (zh) * | 2021-06-21 | 2021-09-07 | 湖北民族大学 | 基于回归方程基础上的藤茶品质评审方法 |
CN113567390A (zh) * | 2021-06-24 | 2021-10-29 | 贵州省茶叶研究所 | 一种基于近红外光谱技术审评红茶的方法及系统 |
CN113945539A (zh) * | 2021-10-18 | 2022-01-18 | 江西农业大学 | 基于gwas筛选后近红外光谱特征波段预测品质的方法及系统 |
CN113945539B (zh) * | 2021-10-18 | 2023-07-04 | 江西农业大学 | 基于gwas筛选后近红外光谱特征波段预测品质的方法及系统 |
CN117541763A (zh) * | 2024-01-05 | 2024-02-09 | 汉中市黑金茶科技有限公司 | 一种基于图像识别的茶叶发酵度测定方法及系统 |
CN117541763B (zh) * | 2024-01-05 | 2024-03-26 | 汉中市黑金茶科技有限公司 | 一种基于图像识别的茶叶发酵度测定方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107860740A (zh) | 一种基于近红外光谱技术的红茶发酵品质的评价方法 | |
Li et al. | Simultaneous determination of six main types of lipid-soluble pigments in green tea by visible and near-infrared spectroscopy | |
Xiaowei et al. | Measurement of total anthocyanins content in flowering tea using near infrared spectroscopy combined with ant colony optimization models | |
Wang et al. | Evaluation of matcha tea quality index using portable NIR spectroscopy coupled with chemometric algorithms | |
Jin et al. | Intelligent evaluation of black tea fermentation degree by FT-NIR and computer vision based on data fusion strategy | |
Yang et al. | Quantitative prediction and visualization of key physical and chemical components in black tea fermentation using hyperspectral imaging | |
Wang et al. | Enhanced quality monitoring during black tea processing by the fusion of NIRS and computer vision | |
Teye et al. | Estimating cocoa bean parameters by FT-NIRS and chemometrics analysis | |
Teye et al. | Cocoa bean and cocoa bean products quality evaluation by NIR spectroscopy and chemometrics: A review | |
Zareef et al. | Prediction of amino acids, caffeine, theaflavins and water extract in black tea using FT-NIR spectroscopy coupled chemometrics algorithms | |
Dong et al. | Prediction of black tea fermentation quality indices using NIRS and nonlinear tools | |
Jiang et al. | Determination of tea polyphenols in green tea by homemade color sensitive sensor combined with multivariate analysis | |
Ouyang et al. | Intelligent evaluation of color sensory quality of black tea by visible-near infrared spectroscopy technology: A comparison of spectra and color data information | |
Cai et al. | Using FTIR spectra and pattern recognition for discrimination of tea varieties | |
US20100297291A1 (en) | Visible/near-infrared spectrum analyzing method and grape fermenting method | |
Li et al. | Monitoring black tea fermentation using a colorimetric sensor array‐based artificial olfaction system | |
CN103278473B (zh) | 白胡椒中胡椒碱及水分含量的测定和品质评价方法 | |
Huang et al. | Qualitative discrimination of Chinese dianhong black tea grades based on a handheld spectroscopy system coupled with chemometrics | |
CN112595692A (zh) | 一种水果总糖含量预测模型的建立方法及水果总糖含量的预测方法 | |
CN106950192A (zh) | 一种基于近红外光谱分析技术的植物蛋白饮料中主要成分含量快速检测的方法 | |
Wang et al. | Development of near‐infrared online grading device for long jujube | |
Hu et al. | Determination of Tibetan tea quality by hyperspectral imaging technology and multivariate analysis | |
CN105973840B (zh) | 一种快速检测酱香型酿酒酒醅理化指标的方法 | |
Wu et al. | Determination of apple varieties by near infrared reflectance spectroscopy coupled with improved possibilistic Gath–Geva clustering algorithm | |
Sheng et al. | Data fusion strategy for rapid prediction of moisture content during drying of black tea based on micro-NIR spectroscopy and machine vision |
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: 20180330 |
|
RJ01 | Rejection of invention patent application after publication |