CN112485304A - 一种基于电导法的牛奶检测方法及*** - Google Patents
一种基于电导法的牛奶检测方法及*** Download PDFInfo
- Publication number
- CN112485304A CN112485304A CN202011466288.3A CN202011466288A CN112485304A CN 112485304 A CN112485304 A CN 112485304A CN 202011466288 A CN202011466288 A CN 202011466288A CN 112485304 A CN112485304 A CN 112485304A
- Authority
- CN
- China
- Prior art keywords
- data
- current
- sample
- milk
- fitting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 66
- 239000008267 milk Substances 0.000 title claims abstract description 61
- 235000013336 milk Nutrition 0.000 title claims abstract description 61
- 210000004080 milk Anatomy 0.000 title claims abstract description 61
- 238000001514 detection method Methods 0.000 title claims abstract description 36
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 24
- 238000000605 extraction Methods 0.000 claims abstract description 13
- 230000004927 fusion Effects 0.000 claims abstract description 9
- 238000007637 random forest analysis Methods 0.000 claims description 28
- 239000013598 vector Substances 0.000 claims description 23
- 208000004396 mastitis Diseases 0.000 claims description 21
- 238000012549 training Methods 0.000 claims description 12
- 238000011156 evaluation Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 6
- 210000001082 somatic cell Anatomy 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 5
- 238000012216 screening Methods 0.000 claims description 5
- 238000002790 cross-validation Methods 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 239000011521 glass Substances 0.000 claims description 2
- 230000000295 complement effect Effects 0.000 claims 1
- 238000010276 construction Methods 0.000 abstract description 2
- 238000011550 data transformation method Methods 0.000 abstract description 2
- 238000001228 spectrum Methods 0.000 abstract description 2
- 239000000523 sample Substances 0.000 description 36
- 241000283690 Bos taurus Species 0.000 description 15
- 235000013365 dairy product Nutrition 0.000 description 10
- 238000003066 decision tree Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 239000007788 liquid Substances 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000013361 beverage Nutrition 0.000 description 1
- 238000009614 chemical analysis method Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 244000144980 herd Species 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006651 lactation Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/06—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a liquid
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Evolutionary Computation (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明为一种基于电导法的牛奶检测方法及***,牛奶样本特征的提取过程,只需要得到单一数据量:电流值,又根据得到的单一数据量的各种特性及其分类目标,提出了电流的数据变换方法,对原始的电流数据进行特征构造,即降低了原始电流数据的冗余性,充分提取和丰富了采集的数据特征。相比于其他高精度方法如获取光谱、超声波值的特征获取方式更为简单,易得,是针对样品牛奶的直接信息获取,所获得的融合特征配合智能算法,能提高分类精度或数据检测结果精度。
Description
技术领域
本发明涉及液体成分检测技术领域,尤其涉及一种基于电导法的牛奶检测方法及***。
背景技术
牛奶作为人类的第二大饮品,乳汁的质量问题一直是人们关注的焦点。奶牛乳腺炎一直是影响牛奶质量与产量的主要因素之一,患有奶牛乳腺炎的奶牛泌乳量会减少,乳汁稀薄,乳汁中出现凝块,乳汁质量下降,甚至于丧失泌乳能力和死亡。
传统的牛奶检测方法例如化学分析法检测步骤繁琐,试剂又破坏了牛奶的成分结构,红外线检测和超声波检测的成本昂贵、对环境要求高。
目前而言,广泛用于鉴定牛奶检测以及牛奶的健康管理是DHI(Dairy HerdImprovement)中心。然而,DHI检测设备少而贵。一份牛奶从牛场送出到检测报告出具,往往耗时一个月之久。
现有电导法确定奶牛乳腺炎方法是选取每次挤奶过程中出现的10个最高电导率之平均值作为电导率的稳定值,当将诊断临界条件为8000-9000μmho/cm之间时,亚临床乳腺炎测诊率可达到89%。实验条件较为苛刻,不易操作,且测诊率不高。
因此,设计一种快速、准确、低成本的基于电导法的牛奶检测方法及***去解决上述技术问题是必要的。
发明内容
本发明的目的是,提供一种基于电导法的牛奶检测方法及***,该方法可以快速、低成本、高精度的进行牛奶检测,可用以分类奶牛乳腺炎。该***简单,便捷。
为解决上述技术问题,本发明采用的技术方案是:
一种基于电导法的牛奶检测方法及***,该方法包括以下内容:
采集不同牛奶样品中一段时间内产生的电导信息作为样本数据,即电流信息。
利用样本数据进行数学特征提取,分别取每个样本数据的最大值、最小值、方差值、平均值和响应曲线的面积值,作为样本数据的A类数据特征;
对样本数据用直线拟合包裹法进行特征提取,即以电压作为纵坐标、电流作为横坐标按照单个个体分别拟合样本数据的数据点,获得拟合之后的斜率和与横坐标电流的截距,再同时通过包裹法获得两条包裹平行线与横坐标电流的截距;
利用交流电特性对样本数据进行拟合处理,获得每个样本的拟合电流曲线,满足三角函数公式
X1sin(X2t+X3)+X4;
式中,X1为幅度、X2为初相、X3为角频率、X4为偏距、t为时间,保留电流的幅度、初相、角频率、偏距这四个参数信息;
将上述通过直线拟合包裹法获得的参数信息及交流电特性获得的参数信息作为样本数据的B类数据特征;
利用SVM-RFE(支持向量机结合后向搜索策略算法)对样本数据进行特征筛选,获得C类数据特征;
电流信息作为单一变量,采集间隔0.02s,连续采集1万个数据,数据量比较大。本发明降低了智能算法模型的电流数据的输入数据量,提升了特征数量,上述三类数据特征作为特征矩阵输入到智能算法模型中,三类特征相互补充,去除了一定的冗余,得到当前牛奶中体细胞数。可以用于后续奶牛乳腺炎分类。
所述智能算法模型为改进深度随机森林算法,该算法的具体过程是:
为A类、B类和C类数据特征赋予权值,构成的融合特征矩阵,输入到改进深度随机森林算法中,进入多粒度扫描阶段;使用一个长度为L的采样窗口对融合特征矩阵进行滑动采样,得到子样本向量;
上述每个子样本向量通过级联模块中的完全随机森林和随机森林分别生成一个表征向量H(x),再分别给两个输出向量进行加权处理,然后叠加生成一个特征向量Zi,也作为后续深度森林的输入;
赋予权值的处理方式为改进的critic权重法
森林输出所述加权处理方式为AHP和熵值法组合权重法;
在模型训练中,对每层深度随机森林进行交叉验证,以防止数据过拟合,停止规则是正确率对比上一层不再提升,获得最终的预测模型。
所述赋予权值的过程,采取了改进的critic权重法,更能体现对比强度和指标冲突性,提高融合效果。
首先进行归一化,
其第f个指标的客观权重Wf为:
所述的随机森林和完全随机森林的输出均为:
式中,H(x)为最终结果,ht为第t个决策树预测结果,Y为真实值,I为度量函数,M为决策树的数量。
对深度随机森林进行训练时,采用以下指标进行综合评判,误差越小,系数越高,表明深度随机森林的性能越好,当各性能指标均满足设定条件时才能停止训练,各评判标准的计算方法如下:
(1)平均绝对误差MAE计算公式:
(2)均方根误差RMSE计算公式:
(3)决定系数R2计算公式:
(4)校正后的决定系数Ra 2计算公式:
一个基于电导法的牛奶检测方法及***,包括:
密封的玻璃容器,用以盛放待测牛奶样品;
两个电极,放置在待测牛奶样品中;
电导分析仪,电导分析仪连接两个电极,用于检测待测牛奶样品中的电流;
数据特征提取模块,用于电导分析仪输入的数据进行数据处理和特征提取;
乳腺炎预测模块,用于对提取的数据特征进行分析预测乳腺炎分类;
数据特征提取模块和乳腺炎预测模块加载在终端中,终端与电导分析仪连接;,
终端用以分析实时采集的牛奶电导数据。
所述终端可以为单片机也可以为电脑端、智能手机端等,电导分析仪由接口总线连接至电脑端。
与现有技术相比,本发明的有益效果是:
本发明牛奶样本特征的提取过程,只需要得到单一数据量:电流值,又根据得到的单一数据量的各种特性及其分类目标,提出了电流的数据变换方法,对原始的电流数据进行特征构造,即降低了原始电流数据的冗余性,充分提取和丰富了采集的数据特征。相比于其他高精度方法如获取光谱、超声波值的特征获取方式更为简单,易得,是针对样品牛奶的直接信息获取,所获得的融合特征配合智能算法,能提高分类精度或数据检测结果精度。
在多粒度扫描阶段,可以有效的处理有关联的数据,将A、B、C三类数据相互融合,得到新的特征向量,明显提高后期模型效果。又针对形成的新的特征向量,采用改进critic权重法给已提取特征增加权重,使特征向量在输入森林前更具有客观性,在一定程度上对特征的重要性进行了梳理,让有效特征的作用效果更强。而通过级联模块中的级联森林产生的输出通过AHP和熵值法组合权重法加权后,不仅利用了先验知识,主观随机性控制在一定范围内,实现主客观赋权中的中正、客观方面,得到的结果更加精确。综上,针对牛奶电流样本对深度随机森林改进了两点:一、在多粒度对深度随机森林提出了新的权重配比模式,在样本特征选择方面,没有改变Bagging的随机又放回抽样的特性,但是改变了权重相等的模式,而是采取了Boosting的优点,针对每一个样本特征向量去配比权重,有利于提高了分类的准确率。二、针对非最后一层森林的输出采取AHP和熵值法组合权重法先对其进行加权处理后,再作为下一层森林的输入,即结合专家经验和客观因素对森林的输出进行权重配比,进一步提高智能算法模型的准确度。为工作人员提供了实时的乳腺炎评估建议,可以及时发现奶牛问题,及时更改饲料和药物配比,有效的减少了损失。
本发明牛奶检测牛奶检测方法及***首先保留了待测牛奶的完好性,不会破坏牛奶的原有结构,不造成资源浪费;其次与现用的DHI检测方法相比检测速度快,操作简单,价格更便宜,可以在任何一个牛场中进行,数据易得,也可以在DHI中心对样本进行初步检测。
牛奶检测方法及***通过电极采集了牛奶的电流信息,通过对电流特征的数据梳理和特征选择,然后输入已经训练好的智能算法中进行分析,通过跟设定好的分类标准进行比较,若超过阈值则发送信息至牛场相关工作人员,进行有效预警,实现了牛场实时便捷的智能监测,监测所需变量少,速度高,成本低,准确率高。乳腺炎可根据DHI中心的评估报告分为四类:A类属于非临床型乳腺炎(体细胞数未超过50万/ml),B类属于轻度乳腺炎(体细胞在50-100万/ml),C类属于重度临床型乳腺炎(体细胞在100-500万/ml),需要进行单独饲喂和挤奶,D类属于机型全身性乳腺炎(体细胞超过500万/ml),需要兽医检查是否有并发症,并采取有效措施。
附图说明书
图1是本发明基于电导法的牛奶数据检测装置的结构示意图;
图2是本发明基于电导法的牛奶数据检测方法的流程示意图;
图3(a)是本发明基于电导法的牛奶数据检测方法中电直线拟合包裹结果图;
图3(b)是本发明基于电导法的牛奶数据检测方法中交流电特性对样本电流数据进行拟合处理结果图;
图4是本发明改进的多粒度扫描结构示意图;
图5是改进的级联森林结构示意图;
图6是不同算法对比的预测值与实际值的残差图。
具体实施方式
下面结合附图对本发明进行进一步详细说明:
牛奶液体本身具有一定的阻抗,患有奶牛乳腺炎的奶牛泌乳量会减少,乳汁稀薄,乳汁中出现凝块,其阻抗数据会发生一定变化,本发明先通过探针加载牛奶样品上的电压信号并在信号线上捕获电流信号来分析。电导分析仪可以通过接口总线连接PC端快速记录采样电流和电压信号。
牛奶的阻抗值会受温度影响。所以应先牛奶样品储存在38℃的恒温水箱中静置5-8分钟。再取样并快速测量电信号,并记录温度。这在一定程度上消除了由温度引起的实验误差。
牛奶的电导率与探针之间的液体电导率有关,即***液体表面的探针的高度和角度。在实验中,在探头上标记,使探头与液体表面之间的距离高度和角度相同,这确保了探头和液体之间的接触面积和位置在每次测量期间不变。
通过上述方法采集到实际电导数据作为样本,进行数据预处理。所述的数据预处理是对数据进行清洗。舍弃其中的不合理数据。结果数据的90%作为训练集,10%作为测试集。最后开始进行特征提取,分为三个部分:
第一部分通过数理知识,提取电流响应数据的最大值、最小值、方差值、平均值和响应曲线的面积值,作为电流数据的A类数据特征;
第二部分提取数据的数据点显示在以为电压纵坐标,电流为横坐标的坐标系中。利用直线拟合出所有数据点的中心线并记录其斜率和与横坐标电流的截距作为数据特征。设置中心线产生两条平行线去包裹95%以上的数据点,同时也剔除了某些异常点,在一定程度上消除了因个人操作产生噪声点。再保留两条平行线与横坐标电流的截距,将其也作为数据特征,最终得到三条关于数据点几何特性的直线,记录三条直线与横坐标电流的截距作为数据特征。
再利用电学特性,由于电压和电流信号都是交流信号,阻抗的值和幅度与相位有关,因此通过三角函数拟合原始信息,得到电流与时间曲线特征。
两者综合得到B类数据特征,即包括三条直线与横坐标电流的截距、拟合直线的斜率、电流幅度、初相、角频率、偏距。
第三部分利用多分类SVM-RFE算法(支持向量机结合后向搜索策略)对样本数据进行筛选得到重要特征,将样本中的原始电流电压数据点作为多分类SVM-RFE算法的输入,多分类结果作为输出,获得C类数据特征。
直线拟合方法为:先利用直线拟合的方法在横坐标为电流,纵坐标为电压的笛卡尔坐标系中观察和拟合所有数据点,记录拟合直线的斜率和截距作为数据特征。
进一步的,拟合线不能完全代表电流数据点的横向特征,因此使用一组平行线来包围95%以上的数据点,同时也剔除了某些异常点,在一定程度上消除了因个人操作产生噪声点。再保留两条线的截距也作为数据特征。
进一步的,电压和电流信号都是交流信号,阻抗的值和幅度与相位有关,因此通过三角函数公式
X1sin(X2t+X3)+X4
并保留幅度、初相、角频率、偏距4个参数信息作为数据特征。
利用SVM-RFE进行特征选择:利用数据挖掘对所有数据点进行清洗和初筛,形成预选特征集合;然后采用SVM-RFE对预选特征集合进一步筛选。
下面进行训练深度森林模型,本发明把原有的多粒度特征扫描进一步优化为critic加权的多粒度扫描,其中每一个扫描结构都由一个随机森林和一个完全随机森林组成。
进一步的,通过采样窗口(定长为K)对所有的特征进行步长为L的滑动采样,采样后会得到S=(P-K)/L+1个K维特征子样本,再由相应的完全随机森林和随机森林去训练特征子样本生成对应C维的表征向量H(x),这样会得到2*S*C维的特征向量Zi,让它作为训练级联随机森林的输入。其中长度K、步长L、所选取的特征向量个数P、森林输出H(x)的维数C、决策树的数量M均统一使用默认值。
扫描后的特征向量其量纲、数量级都不同,且特别是每个类的特征之间数据相关性较强,针对这个问题,本发明采取了改进的critic权重法,基本思想是基于评价指标间的辨别力和冲突性来确立指标的客观权数。在原有critic基础上,首先对特征向量进行归一化,消除特征向量数据之间的量级问题。
其第j个指标的客观权重Wj为
本发明中使用CART的决策树构成森林,并选定基尼指数作为***属性的评价指标。基尼指数系数越小,其子样本的“差异性”越小,相应的该属性值作为***属性时的效果就越理想,所以在选取属性值时,会选取具有最小基尼指数属性值,作为最优的***属性和最有***属性值。针对某样本S,其总方差σ计算方式如下:
其中,μ表示预测结果的均值,yd表示第d个样本的预测结果。
明确了要***的属性后(例如属性为V,某个属性值为u),可以将原有的数据集S根据属性值u划分为两个部分S1和S2,然后计算基尼指数Gain_σ,如下:
Gain-σV,u(S)=σ(S1)+σ(S2)
相应的,针对属性V,可以根据上述属性值划分的方式,计算相应部分的基尼指数,选取最佳的基尼指数,就可以得到属于属性V的最优方案:
对于整个样本集X,可以通过上述方式得到所有属性的最优分类方案,最终得到整个样本集X的最优方案:
在执行完CART决策树后,对森林中的所有树的输出结果进行整理,其中每个森林的输出结果为:
H(x)为一个森林的最终输出,ht为第t个决策树结果,Y为真实结果,I为度量函数,M为决策树的数量。
Ht(xs)=H(x)θs
Ht(xs)表示第t个决策树加综合权重之后的森林输入,θs表示AHP和熵值法组合权重法的综合权重。
设θ1和θ2作为AHP和熵权法得到的两组权重,为了使综合权重θs更好的贴合θ1和θ2,由于预期结果是可知的,属于有监督学习,所以基于最小二乘法原理得到θ1、θ2的值使随机森林的输出值H(x)最接近预期目标。
本发明中,对级联森林进行了优化,将其中的两个完全随机森林和两个随机森林产生的预测值分别与AHP和熵值法组合权值相乘后进行整合,得到2*S*C维的特征向量Z1作为级联森林的输入。
Z1=[H1(x1),H1(x2)…H1(xs),H2(x1),H2(x2)…H2(xs)]T
其中级联森林的每一层包括两个随机森林和两个完全随机森林,每个森林计算方法同上。
生成一个增强型的二维特征向量,这样每层生成一个四维的特征向量,再与之前的原始特征向量融合在一起得到:
用作下一级联森林的输入,以此类推,直至到达训练停止规则。整个多粒度级联森林是自适应调节的,层数不固定,而且在每一级都进行交叉验证去避免过拟合。一级结束后做一次性能测试,然后再继续生成下一级,当扩展一个新的级后,整个级联的性能将在验证集上进行估计,同时调整综合权重AHP部分进行调整,如果没有显着的性能增益,训练过程将终止并选定最优模型和权重参数。深度森林算法默认为3折交叉验证,对最后一层4个级联森林的输出结果H(x)取平均值,最后输出4个C维的数据H(x),然后在这个C维的数据里面取最大的,作为预测值。输出为
采用以下指标进行综合评判预测值的准确性:
(1)平均绝对误差MAE计算公式:
(2)均方根误差RMSE计算公式:
(3)决定系数R2计算公式:
(4)校正后的决定系数Ra 2计算公式:
实例:
选取了5500份牛奶电导数据进行验证。先进行数据筛选(选取完整数据)之后选取得到5000份,样本数据数量为5000,按照前述的特征提取过程进行获取ABC类特征,用于训练智能算法模型,其中90%用以训练集模型使用,剩余10%作为测试集使用。
本实施例中改进的多粒度扫描结构如图4所示,改进的级联森林结构如图5所示,***的超参数设置如表1:
表1
其测试结果如表2,从上述4个指标来对***作评价,得到***的数据梳理和特征提取的有效性,也证明了改进的深度森林算法相比于传统机器学习算法更加有效准确的结论,证明了本发明方法的可靠性和实用性。
表2
首先通过电导数据采集单元,实时快速的提取电流数据,同时通过电导分析仪和数据接口传入PC端,先进行了数据特征提取,再输入改进深度随机森林建立的奶牛乳腺炎预测模型,对乳腺炎进行分类,并对照牛场专家人员的查验,实时有效的得出相关成分报告。实现了奶牛乳腺炎实时在线检测,速度快、准确率高,对牛场的管理和奶牛的健康提供了依据和保障。
本发明未述及之处适用于现有技术。
Claims (8)
1.一种基于电导法的牛奶检测方法,该检测方法包括以下内容:
采集不同牛奶样品中一段时间内产生的电导信息作为样本数据,即电流信息;
利用样本数据进行数学特征提取,分别取每个样本电流数据的最大值、最小值、方差值、平均值和响应曲线的面积值,作为电流数据的A类数据特征;
对样本数据用直线拟合包裹法进行特征提取,即以电压作为纵坐标、电流作为横坐标按照单个个体分别拟合所有样本数据的数据点,获得拟合之后的斜率和与横坐标电流的截距,再同时通过包裹法获得两条包裹平行线与横坐标电流的截距;
利用交流电特性对样本数据进行拟合处理,获得每个样本的拟合电流曲线,满足三角函数公式
X1sin(X2t+X3)+X4;
式中,X1为幅度、X2为初相、X3为角频率、X4为偏距,保留电流的幅度、初相、角频率、偏距这四个参数信息;
将上述通过直线拟合包裹法获得的参数信息及交流电特性获得的参数信息作为电流数据的B类数据特征;
利用SVM-RFE对样本数据进行特征筛选,获得C类数据特征。
2.根据权利要求1所述的检测方法,其特征在于,将上述三类数据特征作为特征矩阵输入到智能算法中,三类特征相互补充,得到当前牛奶中体细胞数;所述智能算法为改进深度随机森林算法,该算法的具体过程是:
将A类、B类和C类数据特征赋予权值,构成的融合特征矩阵,输入到改进深度随机森林算法中,进入多粒度扫描阶段;使用一个长度为L的采样窗口对融合特征矩阵进行滑动采样,得到子样本向量;
上述每个子样本向量通过级联模块中的完全随机森林和随机森林分别生成一个表征向量H(x),再分别给两个向量进行加权处理,然后叠加生成一个特征向量Zi,也作为训练深度森林的输入;
赋予权值处理方式为改进的critic权重法
森林输出所述加权处理方式为AHP和熵值法组合权重法;
在模型训练中,对每层深度随机森林进行交叉验证,以防止数据过拟合,停止规则是正确率对比上一层不再提升,获得最终的预测模型。
4.根据权利要求1所述的检测方法,其特征在于,直线拟合法获得的拟合线不能完全代表电流数据点的横向特征,使用包裹法剔除异常点,包裹法为:一组平行线来包围95%以上的数据点,再保留两条平行线的截距也作为数据特征。
6.一种基于电导法的牛奶数据检测***,其特征在于,该***包括:
密封的玻璃容器,用以盛放待测牛奶样品;
两个电极,放置在待测牛奶样品中;
电导分析仪,电导分析仪连接两个电极,用于检测待测牛奶样品中的电流、电压;
数据特征提取模块,用于电导分析仪输入的数据进行数据处理和特征提取;
乳腺炎预测模块,用于对提取的数据特征进行分析预测乳腺炎分类,
数据特征提取模块和乳腺炎预测模块加载在终端中,终端与电导分析仪连接,
终端用以分析实时采集的牛奶电导数据。
7.根据权利要求6所述的检测***,其特征在于,所述终端为单片机、电脑端、智能手机端,电导分析仪由接口总线连接至电脑端。
8.根据权利要求6所述的检测***,其特征在于,所述数据特征提取模块包括以下内容:
取每个样本电流数据的最大值、最小值、方差值、平均值和响应曲线的面积值,作为电流数据的A类数据特征;
对样本电流数据用直线拟合包裹法进行特征提取,即以电压作为纵坐标、电流作为横坐标按照单个个体分别拟合所有样本数据的数据点,获得拟合之后的斜率和与横坐标电流的截距,再同时通过包裹法获得两条包裹平行线与横坐标电流的截距;
利用交流电特性对样本电流数据进行拟合处理,获得每个样本的拟合电流曲线,满足三角函数公式
X1sin(X2t+X3)+X4;
式中,X1为幅度、X2为初相、X3为角频率、X4为偏距,保留电流的幅度、初相、角频率、偏距这四个参数信息;
将上述通过直线拟合包裹法获得的参数信息及交流电特性获得的参数信息作为样本数据的B类数据特征;
利用SVM-RFE对样本数据进行特征筛选,获得C类数据特征;
上述A类数据特征、B类数据特征、C类数据特征构成最终提取的数据特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011466288.3A CN112485304B (zh) | 2020-12-14 | 2020-12-14 | 一种基于电导法的牛奶检测方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011466288.3A CN112485304B (zh) | 2020-12-14 | 2020-12-14 | 一种基于电导法的牛奶检测方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112485304A true CN112485304A (zh) | 2021-03-12 |
CN112485304B CN112485304B (zh) | 2022-08-09 |
Family
ID=74916719
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011466288.3A Expired - Fee Related CN112485304B (zh) | 2020-12-14 | 2020-12-14 | 一种基于电导法的牛奶检测方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112485304B (zh) |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5416417A (en) * | 1991-05-08 | 1995-05-16 | S.A.E. Afikim | Method and system for investigating mastitis of a cow by measuring electrical conductivity of milk |
US6307362B1 (en) * | 1998-02-17 | 2001-10-23 | Agricultural Instruments Canada Ltd. | Somatic cell analyser |
CN201828539U (zh) * | 2010-09-30 | 2011-05-11 | 李延辉 | 一种奶牛***炎检测器 |
CN103025890A (zh) * | 2010-04-06 | 2013-04-03 | 卡里斯生命科学卢森堡控股 | 疾病的循环生物标志物 |
CN203203956U (zh) * | 2013-05-06 | 2013-09-18 | 北京倍特双科技发展有限公司 | 奶牛乳腺炎检测装置 |
CN104271759A (zh) * | 2011-09-22 | 2015-01-07 | 林纳格生物科学股份有限公司 | 作为疾病信号的同种型谱的检测 |
CN104330443A (zh) * | 2014-09-30 | 2015-02-04 | 山东国家农产品现代物流工程技术研究中心 | 一种预测冷链肉类微生物数量的方法与*** |
CN205607860U (zh) * | 2016-04-12 | 2016-09-28 | 中国农业科学院农业信息研究所 | 一种牛奶体细胞数测量装置 |
CN108363972A (zh) * | 2018-02-05 | 2018-08-03 | 西北农林科技大学 | 一种基于介电频谱技术的牛乳体细胞数快速检测方法 |
CN109711549A (zh) * | 2018-12-27 | 2019-05-03 | 中国农业大学 | 一种基于遗传算法优化bp神经网络的奶牛***炎检测方法 |
CN111007220A (zh) * | 2019-12-17 | 2020-04-14 | 浙江农林大学 | 一种生猪养殖污水水质敏捷监测的方法 |
CN111177650A (zh) * | 2019-12-18 | 2020-05-19 | 国网浙江省电力有限公司绍兴供电公司 | 一种配电网电能质量监测与综合评估***及方法 |
CN111505058A (zh) * | 2020-03-05 | 2020-08-07 | 艾普康(香港)有限公司 | ***炎检测分析仪 |
CN111638249A (zh) * | 2020-05-31 | 2020-09-08 | 天津大学 | 基于深度学习的含水率测量方法及其在油井开采中的应用 |
CN111812158A (zh) * | 2020-07-10 | 2020-10-23 | 桂林优利特医疗电子有限公司 | 一种尿液电导率、渗透压的变频检测方法 |
CN112014461A (zh) * | 2019-05-28 | 2020-12-01 | 滁州市康普检测设备有限公司 | 一种用于牛奶多指标分析的检测***和方法 |
-
2020
- 2020-12-14 CN CN202011466288.3A patent/CN112485304B/zh not_active Expired - Fee Related
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5416417A (en) * | 1991-05-08 | 1995-05-16 | S.A.E. Afikim | Method and system for investigating mastitis of a cow by measuring electrical conductivity of milk |
US6307362B1 (en) * | 1998-02-17 | 2001-10-23 | Agricultural Instruments Canada Ltd. | Somatic cell analyser |
CN103025890A (zh) * | 2010-04-06 | 2013-04-03 | 卡里斯生命科学卢森堡控股 | 疾病的循环生物标志物 |
CN201828539U (zh) * | 2010-09-30 | 2011-05-11 | 李延辉 | 一种奶牛***炎检测器 |
CN104271759A (zh) * | 2011-09-22 | 2015-01-07 | 林纳格生物科学股份有限公司 | 作为疾病信号的同种型谱的检测 |
CN203203956U (zh) * | 2013-05-06 | 2013-09-18 | 北京倍特双科技发展有限公司 | 奶牛乳腺炎检测装置 |
CN104330443A (zh) * | 2014-09-30 | 2015-02-04 | 山东国家农产品现代物流工程技术研究中心 | 一种预测冷链肉类微生物数量的方法与*** |
CN205607860U (zh) * | 2016-04-12 | 2016-09-28 | 中国农业科学院农业信息研究所 | 一种牛奶体细胞数测量装置 |
CN108363972A (zh) * | 2018-02-05 | 2018-08-03 | 西北农林科技大学 | 一种基于介电频谱技术的牛乳体细胞数快速检测方法 |
CN109711549A (zh) * | 2018-12-27 | 2019-05-03 | 中国农业大学 | 一种基于遗传算法优化bp神经网络的奶牛***炎检测方法 |
CN112014461A (zh) * | 2019-05-28 | 2020-12-01 | 滁州市康普检测设备有限公司 | 一种用于牛奶多指标分析的检测***和方法 |
CN111007220A (zh) * | 2019-12-17 | 2020-04-14 | 浙江农林大学 | 一种生猪养殖污水水质敏捷监测的方法 |
CN111177650A (zh) * | 2019-12-18 | 2020-05-19 | 国网浙江省电力有限公司绍兴供电公司 | 一种配电网电能质量监测与综合评估***及方法 |
CN111505058A (zh) * | 2020-03-05 | 2020-08-07 | 艾普康(香港)有限公司 | ***炎检测分析仪 |
CN111638249A (zh) * | 2020-05-31 | 2020-09-08 | 天津大学 | 基于深度学习的含水率测量方法及其在油井开采中的应用 |
CN111812158A (zh) * | 2020-07-10 | 2020-10-23 | 桂林优利特医疗电子有限公司 | 一种尿液电导率、渗透压的变频检测方法 |
Non-Patent Citations (6)
Title |
---|
C. KAMPHUIS: "Decision-tree induction to detect clinical mastitis with automatic milking", 《COMPUTERS AND ELECTRONICS IN AGRICULTURE》 * |
J.M. LUKAS: "A novel method of analyzing daily milk production and electrical conductivity to predict disease onset", 《JOURNAL OF DAIRY SCIENCE》 * |
崔传金: "基于电参数与神经网络的奶牛乳腺炎检测方法", 《农业机械学报》 * |
曹利军: "高压脉冲细胞破壁条件奶牛乳腺炎的电检测技术研究", 《中国优秀硕士学位论文全文数据库 农业科技辑》 * |
李迪星: "颜色空间融合色差法的牛奶细菌含量检测研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
王杰: "电导法牛奶细菌总数检测***设计", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112485304B (zh) | 2022-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ebrahimie et al. | A large-scale study of indicators of sub-clinical mastitis in dairy cattle by attribute weighting analysis of milk composition features: highlighting the predictive power of lactose and electrical conductivity | |
CN113053535A (zh) | 一种医疗信息预测***及医疗信息预测方法 | |
Jiang et al. | Identification of tomato maturity based on multinomial logistic regression with kernel clustering by integrating color moments and physicochemical indices | |
CN112101459B (zh) | 一种基于近红外光光谱特征的动物骨质识别方法及*** | |
CN113310930A (zh) | 高温灭菌奶、巴氏杀菌奶和掺加高温灭菌奶的巴氏杀菌奶的光谱鉴定方法 | |
CN113310936A (zh) | 四种高温灭菌商品牛奶的快速鉴定方法 | |
Mavi et al. | Mango ripeness classification system using hybrid technique | |
Devi et al. | IoT-deep learning based prediction of amount of pesticides and diseases in fruits | |
Loddo et al. | Automatic monitoring cheese ripeness using computer vision and artificial intelligence | |
CN102749420B (zh) | 一种适宜凝胶型蛋白质加工的花生品质测定及其评价方法 | |
CN113310934A (zh) | 骆驼奶中掺加奶牛奶及其掺加比例的快速鉴定方法 | |
CN112485304B (zh) | 一种基于电导法的牛奶检测方法及*** | |
CN113324943A (zh) | 牦牛奶及其掺加奶牛奶的快速鉴别模型 | |
CN112666114A (zh) | 利用光谱识别水牛奶与马奶的方法 | |
CN110991469B (zh) | 一种水果可溶性固形物在线检测方法及*** | |
Gao et al. | Mass detection of walnut based on X‐ray imaging technology | |
Guo et al. | Characterization of different processes lemon slice using electronic tongue | |
CN108562622B (zh) | 一种基于介电特性技术的生鲜羊乳菌落总数快速检测方法 | |
Sanchez et al. | Determination of Sugar Apple Ripeness via Image Processing Using Convolutional Neural Network | |
CN113324942A (zh) | 原料牛奶、高温灭菌奶和掺加高温灭菌奶的原料牛奶的快速鉴定模型 | |
CN113310933A (zh) | 原料水牛奶保存天数的光谱鉴定方法 | |
CN113310938A (zh) | 一种巴氏杀菌的水牛鲜奶和奶牛鲜奶的快速鉴定方法 | |
CN106501325A (zh) | 一种优化传感数据及模式识别对辐照食品的实时快检方法 | |
CN112801172A (zh) | 一种模糊模式识别的白菜农药残留定性分析方法 | |
Rukmana et al. | Mastitis Detection System in Dairy Cow Milk based on Fuzzy Inference System using Electrical Conductivity and Power of Hydrogen Sensor Value |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220809 |