CN110441249B - 基于高光谱图像技术的窖泥总酸预测模型建立的方法 - Google Patents
基于高光谱图像技术的窖泥总酸预测模型建立的方法 Download PDFInfo
- Publication number
- CN110441249B CN110441249B CN201910854782.8A CN201910854782A CN110441249B CN 110441249 B CN110441249 B CN 110441249B CN 201910854782 A CN201910854782 A CN 201910854782A CN 110441249 B CN110441249 B CN 110441249B
- Authority
- CN
- China
- Prior art keywords
- total acid
- pit mud
- visible light
- spectral data
- spectral
- 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.)
- Active
Links
- 239000002253 acid Substances 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005516 engineering process Methods 0.000 title claims abstract description 11
- 238000001228 spectrum Methods 0.000 claims abstract description 48
- 238000012549 training Methods 0.000 claims abstract description 29
- 238000012360 testing method Methods 0.000 claims abstract description 28
- 238000002310 reflectometry Methods 0.000 claims abstract description 12
- 238000012216 screening Methods 0.000 claims abstract description 10
- 238000012937 correction Methods 0.000 claims abstract description 6
- 238000005286 illumination Methods 0.000 claims abstract description 3
- 230000003595 spectral effect Effects 0.000 claims description 38
- 230000000694 effects Effects 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000012512 characterization method Methods 0.000 claims description 3
- 239000010453 quartz Substances 0.000 claims description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 9
- 238000001514 detection method Methods 0.000 description 13
- 238000010521 absorption reaction Methods 0.000 description 6
- 238000007781 pre-processing Methods 0.000 description 6
- 238000000855 fermentation Methods 0.000 description 4
- 230000004151 fermentation Effects 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 125000003118 aryl group Chemical group 0.000 description 3
- 238000002329 infrared spectrum Methods 0.000 description 3
- 238000010563 solid-state fermentation Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- SHZIWNPUGXLXDT-UHFFFAOYSA-N ethyl hexanoate Chemical compound CCCCCC(=O)OCC SHZIWNPUGXLXDT-UHFFFAOYSA-N 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 125000000524 functional group Chemical group 0.000 description 2
- FUZZWVXGSFPDMH-UHFFFAOYSA-N hexanoic acid Chemical compound CCCCCC(O)=O FUZZWVXGSFPDMH-UHFFFAOYSA-N 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 230000002378 acidificating effect Effects 0.000 description 1
- 238000013124 brewing process Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 125000003178 carboxy group Chemical group [H]OC(*)=O 0.000 description 1
- 238000000701 chemical imaging Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 150000002148 esters Chemical class 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000003205 fragrance Substances 0.000 description 1
- 229910052736 halogen Inorganic materials 0.000 description 1
- 150000002367 halogens Chemical class 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- -1 hydrogen ions Chemical class 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000004768 lowest unoccupied molecular orbital Methods 0.000 description 1
- 230000000813 microbial effect Effects 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000003918 potentiometric titration Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 235000021419 vinegar Nutrition 0.000 description 1
- 239000000052 vinegar Substances 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
-
- 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/3577—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种基于高光谱图像技术的窖泥总酸预测模型建立的方法，首先采集近红外以及可见光波段的高光谱图像信息，通过ENVI软件提取感兴趣区域(ROI)的平均光谱信息，经过黑白校正以消除光照不均引起的误差，估算出图像中每个像素点的平均光谱反射率，经过不做处理以及SNV、不同建模方法，分别建立基于全波段和特征光谱的定量预测模型，通过比较分析各拟合模型训练集和测试集决定系数和均方根误差的大小，从而筛选出窖泥总酸定量预测的最佳模型，本发明为白酒酿造产业化的转型升级以及数字化、智能化在线实时监测提供一种强有力的技术支持。
Description
技术领域
本发明涉及一种基于高光谱图像技术的窖泥总酸预测模型建立的方法，属于固态发酵指标检测技术领域。
背景技术
固态发酵技术在我国有着悠久的历史，也是我国特有的酿造工艺。在白酒酿造行业中，窖泥是发酵的基础，其中以己酸菌为主的有益微生物菌群，利用糟醅中的营养物质进行产酸产酯，己酸乙酯是浓香型白酒的主体香味成分，因此窖泥在生产过程中对浓香型白酒质量起着决定性作用，评价窖泥好坏就显得比较重要。没有好的窖泥，就不能生产上乘的优质酒。
目前很多酒厂都面临入窖酸度过高的情况，这会严重影响酒糟的正常发酵，且导致出酒率过低。因为参与糖化和发酵的各种酶系，需要在适宜的酸度条件下才能发挥显著作用。因此入窖前，实时在线监测窖泥酸度，以便及时做出工艺调整，对防止窖泥酸化、老化是很有必要的。目前常用的总酸检测方法是电位滴定法，此方法耗时耗力，影响因素众多，导致误差偏大，并且检测结果往往滞后于实际生产。针对此问题，有学者研究用其他方法进行检测，如专利：高光谱图像技术定量检测固态发酵指标分布差异的方法(CN201510641747)，采集镇江香醋醋醅的全波段高光谱图像，筛选最优特征波数，建立含量预测模型；计算每个像素点对应的pH值；绘制含量分布图；计算目标区域所有像素点pH均值、标准差、NU值，以此定量分析指标含量分布图像的分布差异。
该专利的不足：该专利检测的是pH值，即游离氢离子浓度，因此pH与酸性物质的种类及含量息息相关，发酵基质对光谱的特征吸收性能也通过各种酸类物质所携带的发色基团表现，因为一种物质的光学特性一般是通过官能团表现出来的，不同官能团对不同波段光的吸收强弱存在差异，与分子结构有关。经前期预实验以及资料查询发现，总酸主要是通过羧基表现出其光学特性，并且在近红外波段有较强的吸收，而该专利仅采集了可见光区域的图像和光谱信息；未做全光谱与特征光谱的比较，直接筛选特征波长对应的光谱信息作为模型的输入变量，容易漏掉重要信息；未进行预处理去噪，易将非样本信号纳入计算；建模方式也较单一。
发明内容
本发明的目的是提供一种基于高光谱图像技术的窖泥总酸预测模型建立的方法，通过搭建高光谱数据采集系统，采集不同窖龄不同层次窖泥的光谱数据，对其进行校正处理，提取具有代表性的感兴趣区域每个像素点的平均光谱反射率，利用近红外以及可见光、不做处理以及SNV、全光谱以及特征光谱，结合不同的建模方式，建立窖泥总酸的预测模型，根据模型在训练集和测试集数据上的表现性能，筛选最优模型。本发明为白酒酿造产业化的转型升级以及数字化、智能化在线实时监测提供一种强有力的技术支持。
本发明采取的技术方案是：
基于高光谱图像技术的窖泥总酸预测模型建立的方法，包括以下步骤：
步骤一、光谱采样标定：为了消除光照不均匀等因素引起的误差，需要在采样之前，通过标定板对样品光谱数据进行黑白校正；
标定板颜色要与样品接近，否则容易干扰对样品光谱数据的采集，且反射性能要略优于样品的标定板才能实现对样品原始光谱的校正，此外，样品在可见光和近红外区域下的光谱反射性能不同，因此需要选择不同颜色的标定板，才能获取接下来可分析的数据信息。
针对近红外(900～1700nm)和可见光(400～1000nm)波段的差异，在可见光下选择灰色标定板，近红外光下选择米黄色标定板，同时对两个波段的标定分别进行参数设置，所述参数包括打光峰值、曝光频率、时间、平台移动速度；参数设置详见表格。
步骤二、高光谱数据采集：在设置好的参数下，分别采集n份窖泥样本在近红外和可见光波段的全光谱数据，通过ENVI软件进行黑白校正并将原始光谱信息转化成光谱反射率，提取感兴趣区域(ROI)，保证能代表样本绝大部分信息；
全光谱数据指的是近红外波段下的224个波长对应的光谱数据和可见光波段下的448个波长对应的光谱数据，而原始光谱信息指的是未经过任何处理的光谱数据。
步骤三、表征指标含量的测定：本发明参考GBT 12456-2008中pH电位法，对样品总酸进行检测，为提高模型预测效果，总酸的测定区域要与手动选择的感兴趣区域对应；
步骤四、预处理：考虑到数据采集过程中的噪音、基线漂移等因素可能影响模型精度，将步骤二提取出的光谱信息经过SNV预处理后，得到用于后续分析的数据，同时不做处理对照；
步骤五、数据分集：将近红外以及可见光范围内的窖泥样本分别通过K-S算法分为训练集和测试集；
训练集数据主要用于模型的建立，测试集数据主要用于检验模型的预测效果。
步骤六、筛选特征光谱：运用SPA算法分别从近红外以及可见光的全光谱数据中筛选出与总酸表征指标紧密相关的特征光谱，其与全光谱相对应；因可见光和近红外，不做处理的和SNV预处理得到的用于后续分析的光谱数据不同，经过SPA筛选分别得到4组特征波长。
考虑到总酸对光谱的吸收具有选择性，窖泥内部存在多种化合物对光能的吸收，可能会造成相关波长的波谷或波峰的叠加，为提高模型的稳健性和预测效果，需寻找与总酸紧密相关的特征波长，剔除不相关波长。
步骤七、建模：选取近红外区域的前180个波长、可见光区域全部波长的光谱信息，分别基于全光谱和特征光谱对训练集样本建立不同的定量预测模型(PLSR、SVM、BP)，共计24种模型，通过训练集和测试集的决定系数和均方根误差来评价并筛选出最优预测模型，从而建立一种快速预测窖泥水总酸量的方法。
进一步的，训练集和测试集按照一定的比例进行划分，最终近红外训练集样本72个，测试集样本33个；可见光训练集样本72个，测试集样本36个。
本发明的有益效果：
首次将高光谱成像技术与窖泥总酸实测值相结合，从而建立一种新型、快速、高效的窖泥总酸检测方法，为白酒产业数字化在线监控提供强有力的技术支持。
经过预实验以及前期资料查询发现，总酸所含的分子结构在可见光和近红外波段均有吸收。为筛选与总酸密切相关的特征波段，本发明在可见光(400～1000nm)和近红外(900～1700nm)波段下分别采集窖泥总酸光谱信息，用于后续分析，力争建立窖泥总酸预测的最佳模型。
在前期进行采样标定时，本发明针对特定样品，选择与之相近的，并在光谱反射性能上略优于样品的标定板。此外，考虑近红外和可见光波段的差异，本发明在可见光下选择灰色标定板，近红外下选择米黄色标定板。
考虑到玻璃器皿较强的反光性能，本发明优先选用石英材质的器具用于样本光谱数据的采集，提高模型的精度。
为获得窖泥总酸预测的最优模型，后期做了大量分析工作，包括不做处理以及SNV预处理，全光谱以及特征光谱、不同建模方法等。
附图说明
图1最优模型建立流程图；
图2不同窖泥近红外光谱反射率曲线图；
图3不同窖泥可见光光谱反射率曲线图；
图4近红外光谱训练集模型总酸指标检测效果图；
图5近红外光谱测试集模型总酸指标检测效果图；
图6可见光光谱训练集模型总酸指标检测效果图；
图7可见光光谱测试集模型总酸指标检测效果图。
具体实施方式
窖泥总酸定量预测模型建立的工艺流程如图1所示。
(1)窖泥样本采集：样本采自四川省宜宾市某知名酒厂，采样时工人刨开遗漏的酒糟，以黄水为分界线，用铁铲收集同一水平面四周的窖泥，将其用取样袋封好即成一个样本，同一窖池的采样部位包括窖帽、黄水中和窖底泥三层，共3个样。以此方式获得不同窖龄的窖泥总计108份。
(2)高光谱采集系统：采用芬兰SPECIM FX17系列高光谱相机，2组160W的Y型光纤卤素灯作为光源，利用LUMO Scanner软件控制相机和精密电控载物台进行光谱数据采集。
(3)参数设置及光谱数据采集：通过调节电动载物台运行速率，使其与相机图像扫描频率相匹配，设置高光谱近红外相机下的标定板为米黄色，曝光时间为4.02ms，采集频率为50Hz，载物台的移动速度为16.42mm/sec；高光谱可见光相机下的标定板为灰色，曝光时间为8.00ms，采集频率为50Hz，载物台的移动速度为10.79mm/sec。将装好窖泥并铺平的石英器皿放置于载物台中心位置，在400～1000和900～1700波段下分别采集108份窖泥样本的高光谱原始图像信息，通过ENVI软件手动提取像素为100x100的矩形感兴趣区域(ROI)，并进行校正处理转化为每个像素点的平均光谱反射率，从而分别得到近红外波段105条包含224个波长下的光谱数据以及可见光范围108条包含448个波长的光谱数据。为方便观察，本发明以两个窖池的窖泥为例，阐述不同窖龄窖泥在近红外和可见光波段下光谱反射率曲线的差异，选取了6条分别包含近红外全波段(224个波长)以及可见光全波段(448个波长)的光谱反射率曲线，如图2、3所示，图中一个是6年的新窖，一个是30年的老窖，30年的窖泥光谱反射率普遍高于6年的光谱反射率。不同层级的窖泥光谱反射率曲线也不同，说明高光谱可以对不同年份不同层级的窖泥进行识别。
(4)窖泥总酸检测：为提高模型的精度，取样部位与感兴趣区域对应，总酸的测定方法参考GBT 12456-2008中pH电位法。
(5)光谱数据预处理：光谱数据包含大量的噪声信号，需要对其去噪，经前期大量分析工作发现，标准变量变换(SNV)能有效去除高频噪音、防止基线漂移，优化光谱信号，因此本发明采用的光谱预处理方法是SNV。
(6)样本分集：在数据建模前，通常会将数据集划分为训练集和测试集，训练集数据主要用于模型的建立，测试集数据主要用于检验模型的预测效果，根据预测效果筛选最优模型。通过K-S算法将样本数据划分为训练集和测试集。训练集和测试集按照一定的比例进行划分，最终近红外训练集样本72个，测试集样本33个；可见光训练集样本72个，测试集样本36个。
(7)特征光谱的筛选：总酸对光谱的吸收具有选择性，窖泥内部存在多种化合物对光能的吸收，会造成相关波长的波谷或波峰的叠加，为提高模型的稳健性和预测效果，需寻找与总酸紧密相关的特征波长，剔除不相关波长。本发明采用SPA算法筛选特征波长。因不做处理和SNV预处理得到的用于后续分析的光谱数据不一致，经过SPA筛选分别得到4组特征波长，如表1所示。
表1特征波长统计表
(8)预测模型的建立：为比较获得最佳预测模型，将近红外、可见光、不做处理及预处理、全光谱及特征光谱得到的光谱数据分别结合窖泥总酸实测值建立PLSR、LS-SVM、BP三种模型，通过对最终拟合模型的训练集和测试集均方根误差的比较获得最优模型。鉴于经过不同处理的拟合模型较多，本发明以SNV-SPA-PLSR为例，分别阐述基于近红外和可见光的偏最小二乘模型，近红外的训练集和测试集模型的总酸检测效果如图4、5所示，训练集的决定系数R2＝0.9982，均方根误差RMSE＝0.1058；测试集的决定系数R2＝0.9953，均方根误差RMSE＝0.3023。可见光的训练集和测试集模型的总酸检测效果如图6、7所示，训练集的决定系数R2＝0.9997，均方根误差RMSE＝0.0668；测试集的决定系数R2＝0.9946，均方根误差RMSE＝0.2858。说明在所选取的两种建模方式中，基于可见光的SNV-SPA-PLSR模型在训练集上真实值与预测值的拟合效果更好，但模型在测试集上的性能相对较差，基于近红外的SNV-SPA-PLSR模型的稳健性和泛化性能更好，抗干扰能力更强，更适合应用于窖泥总酸的快速检测。经过不同处理的拟合模型在训练集和测试集上的表现性能如表2所示，从表2可以看出，与BP相比，PLSR与SVR两种建模效果更好。近红外波段下，经过SNV预处理后的特征光谱SVR与可见光下的不做处理特征光谱SVR两种模型的泛化性能较好，精度较高，抗干扰能力较一致。在可见光波段下，经过SNV预处理，分别基于全光谱和特征光谱建立的PLSR两种模型在训练集上各真实值与预测值的拟合效果最好，但对模型的预测效果稍差，模型的稳健性欠缺。因此，综合考虑，通过对比分析各拟合模型的检测效果，最终选定基于可见光波段下的不做处理特征光谱的SVR模型做为窖泥总酸定量预测的最佳模型。
表2各拟合模型在训练集和测试集上的表现性能统计表
Claims (2)
1.基于高光谱图像技术的窖泥总酸预测模型的建立方法，其特征在于，包括以下步骤：
步骤一、光谱采样标定：为了消除光照不均匀引起的误差，在采样之前，通过标定板对光谱数据进行标定；标定板颜色与样品接近，且反射性能略优于样品；针对窖泥样品在可见光波段和近红外波段下的光谱反射性能差异，在可见光波段下选择灰色标定板，在近红外波段下选择米黄色标定板；具体地，设置高光谱近红外相机下的标定板为米黄色，设置高光谱可见光相机下的标定板为灰色；同时对近红外波段和可见光波段的标定分别进行参数设置，设置的参数包括打光峰值、曝光频率、曝光时间、电动载物台移动速度；
步骤二、高光谱数据采集：将装好窖泥并铺平的石英器皿放置于电动载物台中心位置，在设置好的参数下，分别采集n份窖泥样本在近红外和可见光波段的高光谱原始图像信息；通过ENVI软件进行黑白校正，并通过ENVI软件手动选择感兴趣区域，计算感兴趣区域每个像素点的平均光谱反射率，分别得到近红外波段和可见光波段的光谱反射率曲线；
步骤三、总酸含量的测定：使用GB/T 12456-2008中的pH电位法，对样品总酸进行检测；为提高模型预测效果，总酸的测定区域与手动选择的感兴趣区域对应；
步骤四、预处理：将步骤二提取出的光谱反射率曲线进行SNV预处理，同时保留未经SNV预处理的光谱反射率曲线，得到用于后续分析的光谱数据；
步骤五、数据分集：将经过SNV预处理的光谱数据以及未经SNV预处理的光谱数据分别通过K-S算法分为训练集和测试集；
步骤六、筛选特征光谱：运用SPA算法分别在经过SNV预处理的光谱数据以及未经SNV预处理的光谱数据中筛选出与总酸表征指标紧密相关的特征光谱；同时保留未筛选特征光谱的全光谱数据；
步骤七、建模：为获得最佳预测模型，分别利用近红外及可见光波段、未经SNV预处理及经SNV预处理、全光谱及特征光谱的光谱数据结合窖泥总酸实测值建立PLSR、LS-SVM、BP三种模型，共计获得24个模型，通过比较各模型训练集和测试集的决定系数和均方根误差，筛选出最优的窖泥总酸预测模型。
2.如权利要求1所述的基于高光谱图像技术的窖泥总酸预测模型的建立方法，其特征在于，对于近红外波段的光谱数据，训练集的样本个数为72个，测试集的样本个数为33个；对于可见光波段的光谱数据，训练集的样本个数为72个，测试集的样本个数为36个。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910854782.8A CN110441249B (zh) | 2019-09-10 | 2019-09-10 | 基于高光谱图像技术的窖泥总酸预测模型建立的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910854782.8A CN110441249B (zh) | 2019-09-10 | 2019-09-10 | 基于高光谱图像技术的窖泥总酸预测模型建立的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110441249A CN110441249A (zh) | 2019-11-12 |
| CN110441249B true CN110441249B (zh) | 2021-12-28 |
Family
ID=68439959
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910854782.8A Active CN110441249B (zh) | 2019-09-10 | 2019-09-10 | 基于高光谱图像技术的窖泥总酸预测模型建立的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110441249B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110823829A (zh) * | 2019-11-21 | 2020-02-21 | 四川长虹电器股份有限公司 | 一种基于sg校准片的光谱数据补偿方法 |
| CN113920113B (zh) * | 2021-12-06 | 2022-03-18 | 深圳市海谱纳米光学科技有限公司 | 一种对高光谱图像进行点对点白参比校正的方法与系统 |
| CN114441468B (zh) * | 2022-02-10 | 2023-06-27 | 合肥工业大学 | 一种基于近红外-中红外融合算法的窖泥品质检测方法 |
| CN117252875B (zh) * | 2023-11-17 | 2024-02-09 | 山东大学 | 基于高光谱图像的医疗图像处理方法、系统、介质及设备 |
| CN117288692B (zh) * | 2023-11-23 | 2024-04-02 | 四川轻化工大学 | 一种酿酒粮食中单宁含量的检测方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103616335A (zh) * | 2013-12-02 | 2014-03-05 | 江苏大学 | 发酵食品产酸菌产酸能力的快速鉴别方法 |
-
2019
- 2019-09-10 CN CN201910854782.8A patent/CN110441249B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103616335A (zh) * | 2013-12-02 | 2014-03-05 | 江苏大学 | 发酵食品产酸菌产酸能力的快速鉴别方法 |
Non-Patent Citations (2)
| Title |
|---|
| 采用近红外光谱定量分析葡萄酒发酵液中总酸含量和pH值;贾柳君 等;《食品与发酵工业》;20170225;第43卷(第2期);191-195 * |
| 高光谱图像技术检测镇江香醋醋醅理化指标分布研究;申婷婷;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20161115(第11期);1-51 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110441249A (zh) | 2019-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110441249B (zh) | 基于高光谱图像技术的窖泥总酸预测模型建立的方法 | |
| CN110441248B (zh) | 基于高光谱图像技术的窖泥水分预测模型建立的方法 | |
| Cozzolino et al. | Analysis of grapes and wine by near infrared spectroscopy | |
| Wu et al. | Monitoring of fermentation process parameters of Chinese rice wine using attenuated total reflectance mid-infrared spectroscopy | |
| Huang et al. | Rapid and nondestructive prediction of amylose and amylopectin contents in sorghum based on hyperspectral imaging | |
| CN101915738A (zh) | 基于高光谱成像技术的茶树营养信息快速探测方法及装置 | |
| CN109766909B (zh) | 基于谱图融合的海岸环境微塑料老化行为解析方法 | |
| CN113554232A (zh) | 一种农作物产量预测方法及系统 | |
| CN112414950B (zh) | 一种协同高光谱与红外技术的绝缘子等值盐密检测方法 | |
| Fernández-Novales et al. | A feasibility study on the use of a miniature fiber optic NIR spectrometer for the prediction of volumic mass and reducing sugars in white wine fermentations | |
| CN102876816A (zh) | 基于多传感器信息融合的发酵过程状态监测与控制方法 | |
| CN111537469A (zh) | 一种基于近红外技术的苹果品质快速无损检测方法 | |
| CN108801975A (zh) | 一种微型化近红外光谱仪检测酒糟成分的光谱预处理方法 | |
| CN113008815A (zh) | 一种基于高光谱图像信息无损检测酸枣仁中总黄酮的方法 | |
| CN109187578A (zh) | 基于高光谱成像技术的猕猴桃表面缺陷快速无损识别方法 | |
| CN103234923A (zh) | 一种黄酒发酵过程中总糖含量的在线监控方法 | |
| CN113514410A (zh) | 一种基于冠层高光谱技术的夏玉米全生育期氮素利用率垂向分布的实时定量监测方法 | |
| CN113686811A (zh) | 一种基于双传感器的光谱数据处理方法 | |
| Wen et al. | Visualization of chlorophyll content distribution in apple leaves based on hyperspectral imaging technology. | |
| CN113267458A (zh) | 甘薯可溶性蛋白质含量定量预测模型的建立方法 | |
| CN109164069B (zh) | 一种果树叶面病害级别的鉴定方法 | |
| CN113049526A (zh) | 一种基于太赫兹衰减全反射的玉米种子水分含量测定方法 | |
| CN112834454A (zh) | 基于近红外高光谱技术的冬枣检测方法 | |
| CN105241824A (zh) | 高光谱图像技术定量检测固态发酵指标分布差异的方法 | |
| CN109827908B (zh) | 一种利用光谱数据判别红富士苹果腐烂程度的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |