CN113324941A - 原料牛奶保存时间的快速鉴定方法 - Google Patents
原料牛奶保存时间的快速鉴定方法 Download PDFInfo
- Publication number
- CN113324941A CN113324941A CN202110503813.2A CN202110503813A CN113324941A CN 113324941 A CN113324941 A CN 113324941A CN 202110503813 A CN202110503813 A CN 202110503813A CN 113324941 A CN113324941 A CN 113324941A
- Authority
- CN
- China
- Prior art keywords
- milk
- raw milk
- samples
- model
- mid
- 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.)
- Withdrawn
Links
- 235000020185 raw untreated milk Nutrition 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004321 preservation Methods 0.000 title claims abstract description 17
- 235000013336 milk Nutrition 0.000 claims abstract description 30
- 239000008267 milk Substances 0.000 claims abstract description 30
- 210000004080 milk Anatomy 0.000 claims abstract description 30
- 238000012360 testing method Methods 0.000 claims abstract description 24
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 20
- 238000001514 detection method Methods 0.000 claims abstract description 16
- 238000007781 pre-processing Methods 0.000 claims abstract description 15
- 238000012549 training Methods 0.000 claims abstract description 12
- 230000002159 abnormal effect Effects 0.000 claims abstract description 6
- 238000012795 verification Methods 0.000 claims abstract description 4
- 239000000523 sample Substances 0.000 claims description 25
- 238000001228 spectrum Methods 0.000 claims description 16
- 238000012706 support-vector machine Methods 0.000 claims description 12
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 11
- 238000012216 screening Methods 0.000 claims description 10
- 238000007637 random forest analysis Methods 0.000 claims description 8
- 238000002834 transmittance Methods 0.000 claims description 8
- 238000002835 absorbance Methods 0.000 claims description 7
- 238000010521 absorption reaction Methods 0.000 claims description 7
- 239000007788 liquid Substances 0.000 claims description 7
- 235000021243 milk fat Nutrition 0.000 claims description 7
- 102000014171 Milk Proteins Human genes 0.000 claims description 5
- 108010011756 Milk Proteins Proteins 0.000 claims description 5
- 235000021239 milk protein Nutrition 0.000 claims description 5
- 239000007787 solid Substances 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 4
- 239000013307 optical fiber Substances 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 230000004069 differentiation Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 11
- 238000004458 analytical method Methods 0.000 abstract description 8
- 230000003595 spectral effect Effects 0.000 abstract description 7
- 238000004476 mid-IR spectroscopy Methods 0.000 abstract description 6
- 238000000105 evaporative light scattering detection Methods 0.000 abstract 1
- 235000013365 dairy product Nutrition 0.000 description 10
- GUBGYTABKSRVRQ-QKKXKWKRSA-N Lactose Natural products OC[C@H]1O[C@@H](O[C@H]2[C@H](O)[C@@H](O)C(O)O[C@@H]2CO)[C@H](O)[C@@H](O)[C@H]1O GUBGYTABKSRVRQ-QKKXKWKRSA-N 0.000 description 4
- 239000004615 ingredient Substances 0.000 description 4
- 239000008101 lactose Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000010183 spectrum analysis Methods 0.000 description 4
- 235000013305 food Nutrition 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 239000002994 raw material Substances 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 241000283690 Bos taurus Species 0.000 description 2
- 108010010803 Gelatin Proteins 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000000701 chemical imaging Methods 0.000 description 2
- 238000013145 classification model Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 229920000159 gelatin Polymers 0.000 description 2
- 239000008273 gelatin Substances 0.000 description 2
- 235000019322 gelatine Nutrition 0.000 description 2
- 235000011852 gelatine desserts Nutrition 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 235000015097 nutrients Nutrition 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 235000018102 proteins Nutrition 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 102000004169 proteins and genes Human genes 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 241001440840 Mikania micrantha Species 0.000 description 1
- 230000001580 bacterial effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003115 biocidal effect Effects 0.000 description 1
- 235000013351 cheese Nutrition 0.000 description 1
- 235000020186 condensed milk Nutrition 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 235000020247 cow milk Nutrition 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003640 drug residue Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 235000015243 ice cream Nutrition 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000813 microbial effect Effects 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 230000035764 nutrition Effects 0.000 description 1
- 235000016046 other dairy product Nutrition 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 description 1
- 235000013343 vitamin Nutrition 0.000 description 1
- 229940088594 vitamin Drugs 0.000 description 1
- 229930003231 vitamin Natural products 0.000 description 1
- 239000011782 vitamin Substances 0.000 description 1
- 235000020334 white tea Nutrition 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/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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明属于奶品分析技术领域，具体涉及原料牛奶保存时间的快速鉴定技术。本发明的领域与中红外光谱分析领域有关。本发明的主要步骤为：获得不同保存时间的原料牛奶样品中红外光谱数据，对原始中红外光谱进行预处理，去除异常值，将预处理后的数据集分为训练集和测试集，对训练集进行光谱波段选择，在训练集上构建预测模型，使用测试集和验证集对模型进行评估。所述的检测应用：用模型预测原料牛奶是否保存时间超过1天。本发明具有测试速度快、对样品无损、可同时大批量检测等突出优点。
Description
技术领域
本发明属于奶品分析技术领域，具体涉及一种原料牛奶保存时间的快速鉴定技术。本发明的领域与中红外光谱领域相关。
背景技术
牛奶含有蛋白质、脂肪、乳糖和各种维生素及矿物质等一百多种人体需要的营养元素，是一种营养丰富而全面的理想食品。原料牛奶(原料奶)可以制成种类众多的乳制品，包括液体乳类、乳粉类、炼乳类、乳脂肪类、乳冰淇淋类、干酪类、其它乳制品类，乳制品销售额大，涉及生活方方面面，是一种重要的国民食品来源。原料奶位于奶业产业链的最上游，其质量安全将直接影响到乳品的质量与安全，这就要求进行奶制品生产加工的原料奶必须来自健康奶牛，且感官正常，营养成分均衡，微生物指数低，无抗生素、药物残留，不掺假。原料奶在4℃条件下，随着保存时间延长，其酸度上升，pH值下降，乳糖含量下降，当保存时间超过3天时，酸度大于18°T，不符合国家原料奶标准[1]。原料奶的品质监控一直是企业生产环节中的重要一环，传统方法一般通过感官评价、测定酸度或者菌落总数[2]等进行原料奶品质监控，这些方法费时废力，又或需要专业人士来操作。李亮等人通过将近红外光谱技术结合进原料牛奶的酸度、pH、乳糖含量的检测中，用这些指标来判断原料奶的腐败程度[1]。
中红外光谱(MIR)分析是近年来快速发展起来的一种快速、无损、无公害、可多组分同时分析的现代技术，且广泛应用于乳用动物，特别是奶牛的生产性能测定，全国有30多家DHI中心，均应用中红外光谱技术进行检测。目前中红外技术已经应用在生鲜奶的酸度[3]、pH值[4]等保存天数变化相关指标的检测中，具备了直接用中红外测定奶样保存天数的条件，故而可以利用中红外光谱建立原料牛奶保存时间的快速鉴定技术。
发明内容
本发明的目的在于克服现有技术存在的缺陷，提供一种基于中红外光谱分析的原料牛奶保存时间的快速鉴定方法。为了确定最佳的预处理和建模算法组合，本发明对光谱数据使用了包括不处理在内的5种预处理方法，结合2种建模方法，共建立了10个原料牛奶保存时间鉴定模型。并且通过对光谱数据进行Pearson相关性检验和相关性的显著性分析筛选出建模使用的特征光谱。所建立的模型在测试集上的准确率为1.00，在验证集上的准确率为0.90。
本发明的技术方案如下所述：
一种基于中红外光谱分析的原料牛奶保存天数的快速鉴定方法，所述的方法包括以下步骤：
1)选取奶样：分别采集不同保存天数的原料牛奶作为检测样本；
2)采集中红外光谱:利用乳成分检测仪对步骤1)的检测样本进行扫描，通过与其相连的计算机输出样本对应的透光率，得到样本光谱图；
3)对采集的原始中红外光谱数据进行预处理，将光谱数据由透射率(T)转换成吸光度(A)，去除异常值；
4)数据集的划分:将数据集分为训练集和测试集；
5)建模波段的选择：筛选两种牛奶样本的显著差异波段，并去除水的吸收区域；
6)模型的建立与筛选：以训练集奶样的中红外光谱作为输入值，以原料牛奶不同保存天数的类别作为输出值，使用不同光谱预处理方法和不同建模算法在训练集上构建模型，遵循准确性、特异性、灵敏度和AUC等指标较高的原则对模型进行评估及筛选，选择最优的数据预处理和建模算法组合构建模型；
7)最优模型的验证与应用：另取不同保存天数的原料牛奶样本，使用筛选出的最优模型对样本进行鉴别，评估其应用性能；
其中：
步骤2)中采集中红外光谱时，将不同保存天数的原料牛奶样本分别倒入直径3.5cm，高9cm的圆柱形采样管中，保证液面高度大于6cm，然后将其在42℃水浴锅中水浴15-20min，再将固体光纤探头伸到液体中吸样检测；
步骤3)中根据A＝log10(1/T)将透射率(T)转换为吸光度(A)，使用马氏距离和乳脂乳蛋白的百分含量去除异常值，其中，保留光谱马氏距离≤3、乳脂和乳蛋白百分含量在平均值±3.5个标准差范围内的数据；
步骤5)中使用的筛选差异波段的方法为Pearson相关性检验和相关性的显著性检验，去除的水吸收区域为3587.94-2970.66cm-1和1716.81-1543.2cm-1，选择的建模波段是1769.075-1946.327cm-1和2254.591-2817.172cm-1两个波段；
步骤6)中使用的光谱预处理方法为一阶微分(Diff)、标准正态变量变换(SNV)、多元散射校正(MSC)和Savitzky-Golag卷积平滑(简称SG卷积平滑)，使用的建模算法为随机森林(RF)和支持向量机(SVM)。
与现有技术相比，本发明有益效果如下所述：
本发明的发明点在于：选择建模波段是1769.075-1946.327cm-1和2254.591-2817.172cm-1两个波段，预处理方法为一阶微分，建模方法为支持向量机。
本发明实用性强，比近红外检测技术适用性更广。
本发明操作方便、快捷、省时，可以快速得到检测结果，相较菌落总数计数等方法，本发明可以节约大量时间，符合生鲜奶在生产环节上要求及时和迅速的特点。
附图说明
图1：本发明建模波段的平均光谱。即二类原料牛奶在建模波段的平均吸光度图。附图标记说明：图1中横坐标为光谱波数，纵坐标为吸光度；图标实线0类表示保存时间一天的原料牛奶，虚线1类表示保存时间超过一天牛奶；图1(a)为建模波段的全部平均光谱图，图1(b)、1(c)分别为放大的局部光谱图，其中图1(b)的波段范围是3587.94-2970.66cm-1，图1(c)的波段范围是1716.81-1543.2cm-1。
图2：本发明测试集的混淆矩阵。附图标记说明：ROC曲线可以衡量模型在测试集的性能，横坐标为假阳性率，纵坐标为真阳性率，AUC为ROC曲线下与坐标轴围成的面积，取值范围在0.5和1之间，若AUC越接近1.0，表明本发明的方法的真实性越高。
图3：本发明测试集分类概率。附图标记说明：横坐标为预测概率，纵坐标为预测的类别，例如图中左下方圆点表示此样本被分为0类的概率为0.615，即判定为正确分类；图中的所有样本均被正确分类，且大部分概率为0.840-1.000。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
本发明的实施方式所举的实施例中未注明具体条件的试验方法，按照常规方法和条件进行，或按照制造厂商所建议的条件操作(如产品使用说明书)。
实施例1：模型的建立
仪器与设备：选用Delta公司生产的LactoScopeTMFT-B多功能乳品成分分析仪(按产品使用说明书操作)。
具体步骤如下：
(1)采集奶样和测定中红外光谱
从奶牛场采集原料牛奶样本30个，共获得120个检测样本，所有原料牛奶储存在4℃条件下，连续4天测定对所有原料牛奶的中红外光谱(MIR)。将奶样分别倒入直径3.5cm，高9cm的圆柱形样本管中，保证液面高度大于6cm，然后将其在42℃水浴锅中水浴15-20min，再将固体光纤探头伸到液体中吸样检测，通过其软件得到检测样本的透光率。
(2)数据预处理
对所有原料牛奶样本的中红外光谱计算出马氏距离
保留光谱马氏距离≤3，乳脂和乳蛋白百分含量在3.5个标准差
范围内的数据，表1为此过程的样本量变化统计，无异常样本，得到保存时间1天原料牛奶样本30个和保存时间2-4天原料牛奶样本90个，将其按分层抽样法分为训练集(n＝88：保存时间1天原料牛奶21个、保存时间2-4天原料牛奶67个)和测试集(n＝22：保存时间1天原料牛奶5个、保存时间2-4天原料牛奶17个)。
在建模过程中，0代表保存时间1天原料牛奶类，1代表保存时间2-4天原料牛奶类。表2为两类原料牛奶常规乳成分的描述性统计，由表中可得，两类原料牛奶脂肪、蛋白、乳糖和总固形物等常规乳成分均没有显著差异(P>0.05)。结果见表1和表2。
表1剔除异常值时的样本量变化
表2两类原料牛奶常规乳成分的描述性统计
说明：每个参数数值皆用平均值±平均值的标准差表示。其中a，b：在不同的原料牛奶类型中，相同的乳成分中具有不同上标的均值存在显著差异(P<0.05)。
(3)鉴定模型的建立
将光谱数据由透光率(T)转化为吸光度(A)，并去除水的吸收区域，对光谱数据进行Pearson相关性检验[5]，并对相关性进行显著性分析，最终选择1769.075-1946.327cm-1和2254.591-2817.172cm-1两个波段进行建模。图1为建模波段的光谱。
(4)模型的建立
将数据集分为训练集(n＝88)、测试集(n＝22)和验证集(n＝10)。
分别采用一阶微分(Diff)[6]、标准正态变量变换(SNV)[7]、多元散射校正(MSC)[8]和SG卷积平滑[9]对光谱数据进行预处理，同时也与不使用预处理的数据进行比较。
使用随机森林(RF)[10]和支持向量机(SVM)[11]算法利用训练集数据建立分类模型，并对测试集中的样本进行预测。在不同预处理下，RF和SVM算法的建模结果如下表所示。
表3不同预处理下RF和SVM的建模结果
(5)最优模型的筛选和确定
在此鉴别模型中，准确率为正确判断占所有判断的概率，其值越接近1越好；灵敏度表示二分类中某一类被正确分类的比例，其值越接近1越好；特异性表示二分类中另一类被正确分类的概率，其值越接近1越好；AUC是ROC曲线下的面积，其直观地反映了ROC曲线表达的分类能力，AUC＝1代表此分类器为完美分类器，0.5<AUC<1代表此分类器优于随即分类器，0<AUC<0.5则代表此分类器差于随即分类器。由表3可知，只有一阶微分预处理和用支持向量机建模的方法的各项指标达到最优，因此选择一阶微分与支持向量机的组合建立的模型为最优模型。
利用选择的最优分类模型，预测测试集的22个检测样本。以ROC曲线衡量模型在测试集的性能，以假阳性率为横坐标，真阳性率为纵坐标，得到的ROC曲线如图2所示。AUC为ROC曲线下与坐标轴围成的面积，其取值范围在0.5和1之间，AUC越接近1.0，表明方法的真实性越高。由图2可知，本实施例中测试集的AUC为1，说明模型在测试集上的分类效果很好。
图3为测试集中类别分类的概率，例如图3中左下方的绿点表示此样本被分为0类的概率为0.615，且为正确分类。由图3可知，测试集中的所有样本均被正确分类，且大部分概率为0.840-1.000。
实施例2：本发明模型的应用
采用实施例1的测定光谱、数据预处理等技术，对10个不同保存天数原料牛奶样本进行测定中红外光谱和处理，结果如表4所示。
表4模型应用结果
本发明利用筛选出的最优模型对原料牛奶保存天数进行快速鉴定，结果错判样本为1个，即将1类生奶判定到0类中，准确率达到0.90。
本发明检测时间短(只需两分钟)，所需样本量少(20mL)，并且可以同时大批量检测原料牛奶，判断出原料牛奶是否保存时间超过1天。
主要参考文献
[1]李亮，近红外光谱技术在原料奶掺假快速检测及新鲜度检测中的应用研究[D].西北农林科技大学,2010；
[2]卜海峰.生鲜牛乳微生物的检测与分析[J].农村科技,2011(01):59-60；
[3]Calamari L,Gobbi L,Bani P.Improving the prediction ability of FT-MIR spectroscopy to assess titratable acidity in cow's milk.Food Chem.2016Feb1；192:477-84.doi:10.1016/j.foodchem.2015.06.103.Epub 2015Jun 30.PMID:26304375；
[4]Marchi M D,Fagan C C,O'Donnell C P,et al.Prediction of coagulationproperties,titratable acidity,and pH of bovine milk using mid-infraredspectroscopy.[J].Journal of Dairy Science,2009,92(1):423-432；
[5]Pearson K.Determination of the coefficient of correlation[J].Science,1909,30(757):23-25；
[6]张浩等，基于近红外光谱结合机器学习算法检测食用明胶品种溯源的研究[J/OL].河南农业大学学报:1-10[2021-04-29]；
[7]吕美蓉等，光谱数据预处理对潮间带沉积物氮LSSVM模型的影响研究[J].光谱学与光谱分析,2020,40(08):2409-2414；
[8]张浩等，基于近红外光谱结合机器学习算法检测食用明胶品种溯源的研究[J/OL].河南农业大学学报:1-10[2021-04-29]；
[9]马云强等，基于高光谱成像的不同花期薇甘菊特征分析及监测模型研究[J].云南大学学报(自然科学版),2021,43(02):290-298；
[10]韩兆迎等，基于SVM与RF的苹果树冠LAI高光谱估测[J].光谱学与光谱分析,2016,36(03):800-805；
[11]陈书媛等，基于高光谱成像技术的白茶储藏年份判别[J/OL].食品工业科技:1-13[2021-04-29]。
Claims (1)
1.一种原料牛奶保存时间的快速鉴定方法，其特征包括以下步骤：
1)选取奶样：分别采集不同保存天数的原料牛奶作为检测样本；
2)采集中红外光谱:利用乳成分检测仪对步骤1)的检测样本进行扫描，通过与其相连的计算机输出样本对应的透光率，得到样本光谱图；
3)对采集的原始中红外光谱数据进行预处理，将光谱数据由透射率转换成吸光度，去除异常值；
4)数据集的划分:将数据集分为训练集和测试集；
5)建模波段的选择：筛选两种牛奶样本的显著差异波段，并去除水的吸收区域；
6)模型的建立与筛选：以训练集奶样的中红外光谱作为输入值，以原料牛奶不同保存天数的类别作为输出值，使用不同光谱预处理方法和不同建模算法在训练集上构建模型，以准确性、特异性、灵敏度和AUC指标对模型进行评估和筛选，以最优的数据预处理和建模算法组合构建模型；
7)最优模型的验证与应用：另取不同保存天数的原料牛奶样本，以最优模型对样本进行鉴别，评估其应用性能；
其中：
步骤2)中采集中红外光谱时，将不同保存天数的原料牛奶样本分别倒入直径3.5cm，高9cm的圆柱形采样管中，保证液面高度大于6cm，然后将其在42℃水浴锅中水浴15-20min，再将固体光纤探头伸到液体中吸样检测；
步骤3)中根据A＝log10(1/T)将透射率(T)转换为吸光度(A)，使用马氏距离和乳脂乳蛋白的百分含量去除异常值，其中，保留光谱马氏距离≤3、乳脂和乳蛋白百分含量在平均值±3.5个标准差范围内的数据；
步骤5)中使用的筛选差异波段的方法为Pearson相关性检验和相关性的显著性检验，去除的水吸收区域为3587.94-2970.66cm-1和1716.81-1543.2cm-1，选择的建模波段是1769.075-1946.327cm-1和2254.591-2817.172cm-1两个波段；
步骤6)中使用的光谱预处理方法为用一阶微分、标准正态变量变换、多元散射校正和Savitzky-Golag卷积平滑，使用的建模算法为随机森林和支持向量机，以一阶微分和支持向量机构建预测模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110503813.2A CN113324941A (zh) | 2021-05-10 | 2021-05-10 | 原料牛奶保存时间的快速鉴定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110503813.2A CN113324941A (zh) | 2021-05-10 | 2021-05-10 | 原料牛奶保存时间的快速鉴定方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113324941A true CN113324941A (zh) | 2021-08-31 |
Family
ID=77415105
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110503813.2A Withdrawn CN113324941A (zh) | 2021-05-10 | 2021-05-10 | 原料牛奶保存时间的快速鉴定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113324941A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114184573A (zh) * | 2021-11-01 | 2022-03-15 | 华中农业大学 | 牛奶中κ-酪蛋白的中红外快速批量检测方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109142269A (zh) * | 2018-07-26 | 2019-01-04 | 江苏大学 | 一种冰鲜牛肉不同储存时间的快速鉴别方法 |
| CN110441251A (zh) * | 2019-08-09 | 2019-11-12 | Oppo（重庆）智能科技有限公司 | 奶粉剩余保质期的检测方法、装置及计算机可读存储介质 |
| CN112525850A (zh) * | 2020-10-01 | 2021-03-19 | 华中农业大学 | 奶牛奶、马奶、骆驼奶、山羊奶和水牛奶的光谱指纹识别方法 |
-
2021
- 2021-05-10 CN CN202110503813.2A patent/CN113324941A/zh not_active Withdrawn
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109142269A (zh) * | 2018-07-26 | 2019-01-04 | 江苏大学 | 一种冰鲜牛肉不同储存时间的快速鉴别方法 |
| CN110441251A (zh) * | 2019-08-09 | 2019-11-12 | Oppo（重庆）智能科技有限公司 | 奶粉剩余保质期的检测方法、装置及计算机可读存储介质 |
| CN112525850A (zh) * | 2020-10-01 | 2021-03-19 | 华中农业大学 | 奶牛奶、马奶、骆驼奶、山羊奶和水牛奶的光谱指纹识别方法 |
Non-Patent Citations (2)
| Title |
|---|
| CHAO DU ET AL.: ""Genetic Analysis of Milk Production Traits and Mid-Infrared Spectra in Chinese Holstein Population"", 《ANIMALS》 * |
| 张爱武 等: ""p值统计量建模独立性的高光谱波段选择方法"", 《红外与激光工程》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114184573A (zh) * | 2021-11-01 | 2022-03-15 | 华中农业大学 | 牛奶中κ-酪蛋白的中红外快速批量检测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Woodcock et al. | Better quality food and beverages: the role of near infrared spectroscopy | |
| Pu et al. | Ripeness classification of Bananito Fruit (Musa acuminata, AA): A comparison study of visible spectroscopy and hyperspectral imaging | |
| CN107024450A (zh) | 一种基于近红外光谱技术鉴别不同品牌和段数奶粉的方法 | |
| CN107219184A (zh) | 一种应用于产地溯源的肉类鉴别方法及装置 | |
| Melado-Herreros et al. | Postharvest ripeness assessment of ‘Hass’ avocado based on development of a new ripening index and Vis-NIR spectroscopy | |
| CN105044024A (zh) | 一种基于近红外光谱技术对葡萄果实进行无损检测的方法 | |
| CN110987865A (zh) | 一种基于近红外光谱检测无花果品质的方法 | |
| CN109211829A (zh) | 一种基于SiPLS的近红外光谱法测定大米中水分含量的方法 | |
| CN113310936A (zh) | 四种高温灭菌商品牛奶的快速鉴定方法 | |
| Saad et al. | Quality analysis prediction and discriminating strawberry maturity with a hand-held Vis–NIR spectrometer | |
| CN113310930A (zh) | 高温灭菌奶、巴氏杀菌奶和掺加高温灭菌奶的巴氏杀菌奶的光谱鉴定方法 | |
| Sun et al. | Non-destructive detection of blackheart and soluble solids content of intact pear by online NIR spectroscopy | |
| CN113324941A (zh) | 原料牛奶保存时间的快速鉴定方法 | |
| CN113310937A (zh) | 高温灭菌牛奶、巴氏杀菌奶牛鲜奶及奶粉复原牛奶的快速鉴定方法 | |
| CN113310929A (zh) | 高温灭菌奶中掺加豆粉及其掺加比例的光谱鉴定方法 | |
| CN110231302A (zh) | 一种快速测定奇亚籽粗脂肪含量的方法 | |
| CN113324940A (zh) | 特优优质奶、高蛋白特色奶、高乳脂特色奶和普通奶的光谱分级方法 | |
| CN112213281A (zh) | 一种基于透射近红外光谱快速测定淡水鱼新鲜度综合评价的方法 | |
| CN113324943A (zh) | 牦牛奶及其掺加奶牛奶的快速鉴别模型 | |
| CN105675538A (zh) | 一种胡麻饼养分的检测方法 | |
| CN113310933A (zh) | 原料水牛奶保存天数的光谱鉴定方法 | |
| Wang Xue et al. | Monitoring model for predicting maize grain moisture at the filling stage using NIRS and a small sample size. | |
| CN110231306A (zh) | 一种无损、快速测定奇亚籽蛋白质含量的方法 | |
| CN110609011A (zh) | 单籽粒玉米种子淀粉含量近红外高光谱检测方法及系统 | |
| CN113310928A (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 | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20210831 |
|
| WW01 | Invention patent application withdrawn after publication |