CN115600117A - 一种畜牧粪便成分检测与识别方法 - Google Patents

一种畜牧粪便成分检测与识别方法 Download PDF

Info

Publication number
CN115600117A
CN115600117A CN202211616295.6A CN202211616295A CN115600117A CN 115600117 A CN115600117 A CN 115600117A CN 202211616295 A CN202211616295 A CN 202211616295A CN 115600117 A CN115600117 A CN 115600117A
Authority
CN
China
Prior art keywords
value
component detection
sequence
resistor
kth
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
Application number
CN202211616295.6A
Other languages
English (en)
Other versions
CN115600117B (zh
Inventor
何爽
张丽
张田田
赵国强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jining Rencheng District Animal Husbandry And Veterinary Development Center Jining Rencheng District Animal Disease Prevention And Control Center Jining Rencheng District Animal Health And Quarantine Center
Original Assignee
Jining Rencheng District Animal Husbandry And Veterinary Development Center Jining Rencheng District Animal Disease Prevention And Control Center Jining Rencheng District Animal Health And Quarantine Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jining Rencheng District Animal Husbandry And Veterinary Development Center Jining Rencheng District Animal Disease Prevention And Control Center Jining Rencheng District Animal Health And Quarantine Center filed Critical Jining Rencheng District Animal Husbandry And Veterinary Development Center Jining Rencheng District Animal Disease Prevention And Control Center Jining Rencheng District Animal Health And Quarantine Center
Priority to CN202211616295.6A priority Critical patent/CN115600117B/zh
Publication of CN115600117A publication Critical patent/CN115600117A/zh
Application granted granted Critical
Publication of CN115600117B publication Critical patent/CN115600117B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G01N27/06Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a liquid
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/483Physical analysis of biological material
    • G01N33/487Physical analysis of biological material of liquid biological material
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/483Physical analysis of biological material
    • G01N33/487Physical analysis of biological material of liquid biological material
    • G01N33/48707Physical analysis of biological material of liquid biological material by electrical means
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A40/00Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
    • Y02A40/70Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in livestock or poultry

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Chemical & Material Sciences (AREA)
  • Mathematical Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Data Mining & Analysis (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Food Science & Technology (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Computing Systems (AREA)
  • Software Systems (AREA)
  • Mathematical Optimization (AREA)
  • Hematology (AREA)
  • Urology & Nephrology (AREA)
  • Medicinal Chemistry (AREA)
  • Pure & Applied Mathematics (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Computational Linguistics (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Treatment Of Sludge (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

本发明涉及数据处理技术领域,具体涉及一种畜牧粪便成分检测与识别方法,该方法获取畜牧粪便的成分检测过程中畜牧粪便的电阻、测量设备在畜牧粪便中的应力和畜牧粪便的pH值,以得到对应的特征矩阵;根据连续多个成分检测过程的特征矩阵将测量误差相似的两个成分检测过程作为一个匹配对,得到多个匹配对;计算每个匹配对的测量设备误差状态指数,以得到由电阻、应力、pH值和测量设备误差状态指数构成的特征向量;利用特征向量训练LSTM预测网络,以根据预测值对畜牧粪便中检测的实时TK含量进行优化补偿。本发明利用训练好的LSTM预测网络对实时检测的畜牧粪便的TK含量进行优化补偿,使得最终的检测结果更加接近真实值。

Description

一种畜牧粪便成分检测与识别方法
技术领域
本发明涉及数据处理技术领域,具体涉及一种畜牧粪便成分检测与识别方法。
背景技术
各种畜牧在养殖过程中常常会产生大量的粪尿,而这些粪尿中存在大量农作物生长所必需的氮磷钾等营养成分和大量有机质。目前畜牧粪便作为有机肥在农业生产中已得到了广泛应用,但如果不合理使用畜牧粪便或连续过量使用,不仅起不到应有的效果,还会影响农作物生长和人类健康,造成生态环境的二次污染,因此对于畜牧粪便的成分检测非常重要。
目前,对于畜牧粪便的成分检测的常用方法是:通过检测畜牧粪便的电阻率,选择畜牧种类和肥料成分种类,将检测到的电阻率、选择的畜牧种类和肥料成分种类作为输入来选择预测估算模型,利用预测估算模型对输入的数据进行分析,以根据分析结果计算畜牧粪便中的肥料成分含量。但是考虑到检测电阻率的传感器,随着长时间的使用,其稳定性和灵敏度都会下降,或者传感器的探针没有很完美的紧贴被测粪便,进而导致传感器数据存在较大误差,进而使得成分检测结果也存在较大误差。
发明内容
为了解决上述技术问题,本发明的目的在于提供一种畜牧粪便成分检测与识别方法,所采用的技术方案具体如下:
在畜牧粪便的成分检测过程中,根据每个采样时刻的畜牧粪便的电阻、测量设备在畜牧粪便中的应力和畜牧粪便的pH值,分别对应得到电阻序列、应力序列和pH值序列;根据所述电阻序列、应力序列和pH值序列构建所述成分检测过程中的特征矩阵;
获取连续多个所述成分检测过程中的特征矩阵,基于明氏距离计算任意两个所述特征矩阵之间的相似程度,基于所述相似程度利用KM算法将测量误差相似的两个所述成分检测过程作为一个匹配对,得到多个匹配对;根据当前匹配对所对应两个所述成分检测过程的上一个所述成分检测过程的所述特征矩阵,计算测量设备误差状态指数,将当前匹配对所对应两个所述成分检测过程中每个采样时刻的电阻、应力、pH值和所述测量设备误差状态指数构成对应采样时刻的特征向量;
获取每个匹配对的所述测量设备误差状态指数,得到连续多个所述成分检测过程对应每个采样时刻的所述特征向量,利用所述特征向量训练LSTM预测网络,以获取实时采样时刻下实时电阻、实时应力和实时pH值所对应的所述测量设备误差状态指数的预测值,利用预测值对畜牧粪便中检测的实时TK含量进行优化补偿,得到实际TK含量。
进一步的,所述特征矩阵的构建方法,包括:
计算所述电阻序列中第K个电阻与前面第N个电阻之间的差值绝对值,K和N都为正整数,K大于N,以第K个电阻为底数、对应差值绝对值为幂指数得到的结果作为第K个电阻的第一误差因子,组成一维的第一误差因子序列;
计算所述应力序列中第K个应力与前面第N个应力之间的差值绝对值,以第K个应力为底数、对应差值绝对值为幂指数得到的结果作为第K个应力的第二误差因子,组成一维的第二误差因子序列;
计算所述pH值序列中第K个pH值与前面第N个pH值之间的差值绝对值,以第K个pH值为底数、对应差值绝对值为幂指数得到的结果作为第K个pH值的第三误差因子,组成一维的第三误差因子序列;
将所述电阻序列中的第K个电阻所对应的前N个电阻作为一个子集合,获取子集合中的最大电阻和最小电阻,分别计算第K个电阻与最大电阻之间的第一电阻差值绝对值、第K个电阻与最小电阻之间的第二电阻差值绝对值,以第二电阻差值绝对值为分子、第一电阻差值绝对值为分母得到对应的比值,将该比值作为第K个采样时刻的测量质量下降因子,组成一维的测量质量下降因子序列;其中第K个电阻是在第K个采样时刻采集的;
将所述第一误差因子、所述第二误差因子、所述第三误差因子和所述测量质量下降因子序列组成所述成分检测过程中的4维特征矩阵。
进一步的,所述基于明氏距离计算任意两个所述特征矩阵之间的相似程度的方法,包括:
分别计算两个所述特征矩阵中每个维度下的明氏距离,结合所有维度的明氏距离得到所述相似程度。
进一步的,所述测量设备误差状态指数的获取方法,包括:
基于所述成分检测过程的先后顺序,分别获取当前匹配对所对应的两个所述成分检测过程相邻上一个所述成分检测过程的所述特征矩阵,将两个相邻的上一个所述成分检测过程的所述特征矩阵之间的相似程度作为当前匹配对的所述测量设备误差状态指数。
进一步的,所述利用预测值对畜牧粪便中检测的实时TK含量进行优化补偿,得到实际TK含量的方法,包括:
获取实时TK含量与预测值之间的第一乘积;基于历史TK含量获取历史TK含量的中值,获取历史TK含量的中值与第一优化值之间的第二乘积,第一优化值与预测值之和为1,将第一乘积和第二乘积之间的均值作为实时TK含量。
本发明实施例至少具有如下有益效果:(1)获取成分检测过程中畜牧粪便的电阻、测量设备在畜牧粪便中的应力和畜牧粪便的pH值,根据连续测量过中的数据得到每个成分检测过程对应测量设备测量误差的特征矩阵,基于特征矩阵的相似性将测量误差相似的两个成分检测过程作为一个匹配对,根据匹配对中的两个成分检测过程的上一个成分检测过程的特征矩阵之间的差异,分析对应匹配对的测量设备误差状态指数,将测量设备误差状态指数作为对应成分检测过程中每个采样时刻的测量设备的测量误差指标,以保证了LSTM预测网络的训练数据集的有效性和准确性,进而利用时序上每个采样时刻的电阻、应力、pH值和对应的测量设备误差状态指数对LSTM预测网络进行训练,使得LSTM预测网络的预测结果更加准确。
(2)利用训练好的LSTM预测网络获取实时采样时刻下的测量设备误差状态指数,利用测量设备误差状态指数对实时检测的畜牧粪便的TK含量进行优化补偿,使得最后的检测结果更加接近真实值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案和优点,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它附图。
图1为本发明一个实施例提供的一种畜牧粪便成分检测与识别方法的步骤流程图。
具体实施方式
为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例,对依据本发明提出的一种畜牧粪便成分检测与识别方法,其具体实施方式、结构、特征及其功效,详细说明如下。在下述说明中,不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外,一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。
下面结合附图具体的说明本发明所提供的一种畜牧粪便成分检测与识别方法的具体方案。
本方案所针对的具体场景为:畜牧粪便的自动化成分检测流程,此流程包括送料、成分检测、导料和堆垛。
请参阅图1,其示出了本发明一个实施例提供的一种畜牧粪便成分检测与识别方法的步骤流程图,该方法包括以下步骤:
步骤S001,在畜牧粪便的成分检测过程中,根据每个采样时刻的畜牧粪便的电阻、测量设备在畜牧粪便中的应力和畜牧粪便的pH值,分别对应得到电阻序列、应力序列和pH值序列;根据电阻序列、应力序列和pH值序列构建成分检测过程中的特征矩阵。
具体的,对于畜牧粪便成分检测,一般是在批量的回收粪便过程中,对其进行抽样或连续采样检测,例如在传送带安装设备,探针持续采样。
首先,在畜牧粪便的成分检测过程中,获取每个采样时刻的畜牧粪便的电阻,以构成电阻序列。
具体的,电化学法是基于检测物在电化学反应过程中产生的电流、电阻或电位的变化进行检测的方法,在传统实践中,保证准确的电化学法的核心要义是电极材料和探针的选择,电极材料和探针能否紧贴被测粪便或粪便的固液性状决定了测量设备的稳定性和灵敏度等。故利用测量设备采集畜牧粪便的电阻,此处得到的是等效电阻,其中畜牧粪便的电阻的获取方法是公知技术,本方案中不再赘述。
优选的,本发明实施例中使用常规的插入粪便的接触式探头作为测量设备,每测量完成一个得到一个数值,数值单位为kΩ,因此将畜牧粪便的成分检测过程中每个采样时刻的电阻r构成电阻序列R。
然后,在畜牧粪便的成分检测过程中,获取每个采样时刻下测量检测畜牧粪便时,测量设备探入过程中受到的各向应力,构成应力序列。
具体的,因电阻测量时由于粪便过干裂开、撞击到异物、内部溏心等原因,不能保证粪便与测量设备的探头完全接触,进而会出现电阻数据的采集误差。为了提高粪便成分检测的准确性,引入测量设备采集时的各向应力。
为了测量电阻,测量设备的探针性状如同Π形,且是一种直杆状的较细的插入式结构,通过上下导轨移动在完成穿透粪便的过程后具有一定的应力,使得探头能够紧贴粪便,因此通过确定测量设备的应力,以根据应力差异表征测量设备的采集误差。
鉴于粪便质地柔软,基于应力传感器获取所有的轴向应力分量和径向应力分量,将最大的应力分量作为对应采样时刻的测量设备在畜牧粪便中的应力p,进而得到畜牧粪便的成分检测过程中的应力序列P。
最后,在畜牧粪便的成分检测过程中,获取每个采样时刻的畜牧粪便的pH值,以构成pH值序列。
因粪便的来源不同,动物的消化系统特性的差异,导致pH指标不同,虽然不影响其正常施肥使用,但酸碱度能够体现粪便的金属盐和氨氮的理化特性,以便加以区分,故在畜牧粪便的成分检测过程中,利用测量设备采集每个采样时刻对应畜牧粪便的pH值,记为f,构成pH值序列F。
进一步的,为了能够准确评估畜牧粪便的测量设备的误差情况,通过对采集完成电阻序列R、应力序列P和pH值序列F进行计算与分析,能够得到畜牧粪便的测量设备的误差趋势。
计算电阻序列中第K个电阻与前面第N个电阻之间的差值绝对值,K和N都为正整数,K大于N,以第K个电阻为底数、对应差值绝对值为幂指数得到的结果作为第K个电阻的第一误差因子,组成一维的第一误差因子序列。
作为一个示例,令N为10,则第一误差因子的计算公式为:
Figure 353253DEST_PATH_IMAGE001
其中,
Figure 756552DEST_PATH_IMAGE002
为电阻序列中第K个电阻的第一误差因子;
Figure 687599DEST_PATH_IMAGE003
为电阻序列中的第K个电 阻;
Figure 628879DEST_PATH_IMAGE004
为绝对值函数;
Figure 303574DEST_PATH_IMAGE005
为电阻序列中第K个电阻之前的第10个电阻。
需要说明的是,K
Figure 512445DEST_PATH_IMAGE006
10,也即是从电阻序列中的第10个电阻开始计算第一误差因 子,由于前期畜牧粪便不会出现测量设备的测量误差,因此默认第10个电阻之前的电阻没 有误差,进而将计算当前采样时刻的电阻和前面第10次采样时刻的电阻的差值,差值越大, 说明当前采样时刻对应电阻的测量误差越大。
同理,计算应力序列中第K个应力与前面第N个应力之间的差值绝对值,K和N都为正整数,K大于N,以第K个应力为底数、对应差值绝对值为幂指数得到的结果作为第K个应力的第二误差因子,组成一维的第二误差因子序列。
作为一个示例,令N为10,则第二误差因子的计算公式为:
Figure 598082DEST_PATH_IMAGE007
其中,
Figure 511811DEST_PATH_IMAGE008
为应力序列中第K个应力的第二误差因子;
Figure 475350DEST_PATH_IMAGE009
为应力序列中的第K个应 力;
Figure 853242DEST_PATH_IMAGE004
为绝对值函数;
Figure 126092DEST_PATH_IMAGE010
为应力序列中第K个应力之前的第10个应力。
需要说明的是,K
Figure 41964DEST_PATH_IMAGE006
10,也即是从应力序列中的第10个应力开始计算第二误差因 子,由于前期畜牧粪便不会出现测量设备的测量误差,因此默认第10个应力之前的应力没 有误差,进而将计算当前采样时刻的应力和前面第10次采样时刻的应力的差值,差值越大, 说明当前采样时刻对应应力的测量误差越大。
同理,计算pH值序列中第K个pH值与前面第N个pH值之间的差值绝对值,K和N都为正整数,K大于N,以第K个pH值为底数、对应差值绝对值为幂指数得到的结果作为第K个pH值的第三误差因子,组成一维的第三误差因子序列。
作为一个示例,令N为10,则第二误差因子的计算公式为:
Figure 792882DEST_PATH_IMAGE011
其中,
Figure 970486DEST_PATH_IMAGE012
为pH值序列中第K个pH值的第三误差因子;
Figure 679816DEST_PATH_IMAGE013
为pH值序列中的第K个pH 值;
Figure 614143DEST_PATH_IMAGE004
为绝对值函数;
Figure 434331DEST_PATH_IMAGE014
为pH值序列中第K个pH值之前的第10个pH值。
需要说明的是,K
Figure 475231DEST_PATH_IMAGE006
10,也即是从pH值序列中的第10个pH值开始计算第三误差因 子,由于前期畜牧粪便不会出现测量设备的测量误差,因此默认第10个pH值之前的pH值没 有误差,进而将计算当前采样时刻的pH值和前面第10次采样时刻的pH值的差值,差值越大, 说明当前采样时刻对应pH值的测量误差越大。
随着采样次数的增加,测量设备的测量数据的误差程度也会随之增加,相对应检 测的畜牧粪便的检测质量也会逐渐下降,因此基于电阻序列中每个采样时刻下的电阻差异 计算对应采样时刻下的测量质量下降因子,具体方法为:将电阻序列中的第K个电阻所对应 的前N个电阻作为一个子集合,获取子集合中的最大电阻和最小电阻,分别计算子集合中第 K个电阻与最大电阻之间的第一电阻差值绝对值、第K个电阻与最小电阻之间的第二电阻差 值绝对值,以第二电阻差值绝对值为分子、第一电阻差值绝对值为分母得到对应的比值,由 于电阻序列中的第K个元素是在第K个采样时刻下采集的,因此将该比值作为第K个采样时 刻的测量质量下降因子
Figure 89883DEST_PATH_IMAGE015
,同理K
Figure 714768DEST_PATH_IMAGE016
10,进而得到一维的测量质量下降因子序列。
因畜牧粪便的检测质量时的测量设备的误差分析需要从多方面的数据来作支撑, 单一方面的数据不能得到有力支持,因此结合电阻的第一误差因子序列、应力的第二误差 因子序列、pH值的第三误差因子序列和测量质量下降因子序列构成畜牧粪便的成分检测过 程中的特征矩阵
Figure 869806DEST_PATH_IMAGE017
,该矩阵包括4个维度,一个维度表示一行,第一维度为第一误差因子,第 二维度为第二误差因子,第三维度为第三误差因子,第四维度为测量质量下降因子,其中该 特征矩阵更加能够综合体现对应畜牧粪便的成分检测过程中因测量设备误差导致的畜牧 粪便检测结果的变化情况。
步骤S002,获取连续多个成分检测过程中的特征矩阵,基于明氏距离计算任意两个特征矩阵之间的相似程度,基于相似程度利用KM算法将测量误差相似的两个成分检测过程作为一个匹配对,得到多个匹配对;根据当前匹配对所对应两个成分检测过程的上一个成分检测过程的特征矩阵,计算测量设备误差状态指数,将当前匹配对所对应两个成分检测过程中每个采样时刻的电阻、应力、pH值和测量设备误差状态指数构成对应采样时刻的特征向量。
具体的,成分检测系统是面向于批次成分检测中所有的成分检测过程,因此利用步骤S001的方法,获取连续多个成分检测过程中的特征矩阵U。
由于同一批次中不同成分检测过程对畜牧粪便的连续测量都是不同的,因此会出现在相同连续测量段内的特征矩阵是不同,故对任意两个成分检测过程所对应的特征矩阵计算其差异性。
由于特征矩阵是由多个不同方面的因素构成的,存在量纲,故此处采用明可夫斯基距离来表征两个特征矩阵之间的差异性,其中,明可夫斯基距离是欧式距离的推广,是对多个距离度量公式的概括性的表述,则基于明氏距离得到任意两个特征矩阵之间的相似程度的计算公式为:
Figure 280059DEST_PATH_IMAGE018
Figure 78994DEST_PATH_IMAGE019
其中,
Figure 941908DEST_PATH_IMAGE020
为批次内成分检测过程A和成分检测过程B之间对应特征矩阵的 相似程度;n为矩阵的维度;
Figure 87587DEST_PATH_IMAGE021
为成分检测过程A和成分检测过程B对应特征矩阵的第n维度 之间的明氏距离;
Figure 352347DEST_PATH_IMAGE022
为成分检测过程A的特征矩阵中第n维度的第i个数值;
Figure 325113DEST_PATH_IMAGE023
为成分检测 过程B的特征矩阵中第n维度的第i个数值;
Figure 675323DEST_PATH_IMAGE024
为常数。
明氏距离越大,说明这两个成分检测过程中的测量设备误差情况越不相同,对应两个成分检测过程之间对应特征矩阵的相似程度越小,当相似程度趋于0时,可认为这两个成分检测过程中的测量设备误差情况不相同。
由上述计算出的相似程度可以得知同一批次下不同的成分检测过程的畜牧粪便的测量设备误差情况是不同的,且不同成分检测过程之间的差异性都比较大,对于不同测量误差情况的畜牧粪便,可以使用KM算法将测量误差相似的成分检测过程分配为相互对照的匹配对。
将任意两个成分检测过程之间的相似程度D投入到KM算法中,将相似程度最大的两个成分检测过程进行匹配,得到畜牧粪便的测量设备误差最为接近的两个成分检测过程,将其作为一个匹配对,进而将多个成分检测过程进行了配对,得到了多个匹配对。
所配对的成分检测过程的畜牧粪便的测量设备误差是相同,若出现孤立的样本,则将其在后续的数据处理中视为一个匹配对,但匹配对中是两个同样的孤立样本。
基于KM算法,根据每个成分检测过程中特征矩阵对多个成分检测过程进行匹配的目的是找到测量结果相似的数据,从而对其分析测量设备误差状态,即两个匹配对之间若测量记录中误差较大,又因此匹配成一对,则很有可能是探头出现了误差。
根据匹配对中的两个成分检测过程的特征矩阵建立测量设备误差状态指数W,则测量设备误差状态指数W的构建方法为:以一个匹配对为例,基于成分检测过程的先后顺序,分别获取当前匹配对所对应的两个成分检测过程相邻上一个成分检测过程的特征矩阵,进而将两个相邻的上一个成分检测过程的特征矩阵之间的相似程度作为当前匹配对的测量设备误差状态指数。
需要说明的是,获取当前匹配对中两个成分检测过程相邻上一个成分检测过程的特征矩阵,进而计算两者之间的相似程度的目的是为了比较两个相邻上一个成分检测过程的特征矩阵对应测量设备误差是否一致,若不一致,则认为当前匹配对中两个成分检测过程相邻上一个成分检测过程的测量状态虽然不同,但最后导致当前匹配对中两个成分检测过程的测量状态又相同,则认为探头出现了较明显的因插入状态受到影响而导致的误差。
步骤S003,获取每个匹配对的测量设备误差状态指数,得到连续多个成分检测过程对应每个采样时刻的特征向量,利用特征向量训练LSTM预测网络,以获取实时采样时刻下实时电阻、实时应力和实时pH值所对应的测量设备误差状态指数的预测值,利用预测值对畜牧粪便中检测的实时TK含量进行优化补偿,得到实际TK含量。
具体的,利用步骤S002的方法,获取每个匹配对所对应的测量设备误差状态指数W,且每个匹配对所对应的两个成分检测过程共用一个测量设备误差状态指数W,进而能够得到同一批次下连续多个成分检测过程对应每个采样时刻的特征向量。
将时序上的特征向量作为LSTM预测网络的训练集,利用训练集训练LSTM预测网络,且以特征向量中的电阻、应力和pH值作为LSTM预测网络的输入数据、特征向量中的测量设备误差状态指数作为LSTM预测网络的输出数据。由于LSTM预测网络是一种端到端、序列到序列的深度神经网络,基于这样的训练集可以让LSTM预测网络基于电阻、应力和pH值的连续变化特征预测对应的测量设备误差状态指数。
需要说明的是,LSTM预测网络是公知技术,本方案中不再赘述。
利用步骤S001的采样方法,能够获取实时采样时刻下的实时电阻、实时应力和实时pH值,将实时电阻、实时应力和实时pH值输入训练好的LSTM预测网络能够得到实时采样时刻下对应的测量设备误差状态指数的预测值。
基于实时电阻能够得到畜牧粪便中检测的实时TK含量,进而利用预测值对实时TK含量进行优化补偿,以得到实际TK含量。
作为一个示例,对于牛而言,其畜牧粪便中的TK含量的计算公式为TK=0.479/r+ 0.437;首先将实时电阻代入TK含量的计算公式中能够得到实时TK含量,然后基于历史计算 的历史TK含量获取历史TK含量的中值,利用历史TK含量的中值和测量设备误差状态指数的 预测值对实时采样时刻的实时TK含量进行优化补偿,得到实际TK含量,其优化补偿公式为:
Figure 890272DEST_PATH_IMAGE025
,其中,
Figure 9538DEST_PATH_IMAGE026
为实际TK含量,
Figure 624977DEST_PATH_IMAGE027
为实时
Figure 462483DEST_PATH_IMAGE027
含量,
Figure 12282DEST_PATH_IMAGE028
为 预测值,
Figure 454896DEST_PATH_IMAGE029
为历史TK含量的中值。
需要说明的是,预测值越大,说明检测的实时TK含量的误差越大,畜牧粪便的识别 结果越不准确,进而需要对实时TK含量进行更多的补偿,反之,预测值越小,说明检测的实 时TK含量的误差越小,畜牧粪便的识别结果越准确,进而对实时TK含量的补偿越少,因此预 测值与实际TK含量之间呈正相关关系;历史TK含量的中值
Figure 300623DEST_PATH_IMAGE029
越大,预测值对中值
Figure 625425DEST_PATH_IMAGE029
的补偿越少,对应得到的实际TK含量越接近真实值。
综上所述,本发明实施例提供了一种畜牧粪便成分检测与识别方法,该方法获取畜牧粪便的成分检测过程中畜牧粪便的电阻、测量设备在畜牧粪便中的应力和畜牧粪便的pH值,以得到对应的特征矩阵;根据连续多个成分检测过程的特征矩阵将测量误差相似的两个成分检测过程作为一个匹配对,得到多个匹配对;计算每个匹配对的测量设备误差状态指数,以得到由电阻、应力、pH值和测量设备误差状态指数构成的特征向量;利用特征向量训练LSTM预测网络,以根据预测值对畜牧粪便中检测的实时TK含量进行优化补偿。本发明利用训练好的LSTM预测网络对实时检测的畜牧粪便的TK含量进行优化补偿,使得最终的检测结果更加接近真实值。
需要说明的是:上述本发明实施例先后顺序仅仅为了描述,不代表实施例的优劣。且上述对本说明书特定实施例进行了描述。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种畜牧粪便成分检测与识别方法,其特征在于,该方法包括以下步骤:
在畜牧粪便的成分检测过程中,根据每个采样时刻的畜牧粪便的电阻、测量设备在畜牧粪便中的应力和畜牧粪便的pH值,分别对应得到电阻序列、应力序列和pH值序列;根据所述电阻序列、应力序列和pH值序列构建所述成分检测过程中的特征矩阵;
获取连续多个所述成分检测过程中的特征矩阵,基于明氏距离计算任意两个所述特征矩阵之间的相似程度,基于所述相似程度利用KM算法将测量误差相似的两个所述成分检测过程作为一个匹配对,得到多个匹配对;根据当前匹配对所对应两个所述成分检测过程的上一个所述成分检测过程的所述特征矩阵,计算测量设备误差状态指数,将当前匹配对所对应两个所述成分检测过程中每个采样时刻的电阻、应力、pH值和所述测量设备误差状态指数构成对应采样时刻的特征向量;
获取每个匹配对的所述测量设备误差状态指数,得到连续多个所述成分检测过程对应每个采样时刻的所述特征向量,利用所述特征向量训练LSTM预测网络,以获取实时采样时刻下实时电阻、实时应力和实时pH值所对应的所述测量设备误差状态指数的预测值,利用预测值对畜牧粪便中检测的实时TK含量进行优化补偿,得到实际TK含量。
2.如权利要求1所述的一种畜牧粪便成分检测与识别方法,其特征在于,所述特征矩阵的构建方法,包括:
计算所述电阻序列中第K个电阻与前面第N个电阻之间的差值绝对值,K和N都为正整数,K大于N,以第K个电阻为底数、对应差值绝对值为幂指数得到的结果作为第K个电阻的第一误差因子,组成一维的第一误差因子序列;
计算所述应力序列中第K个应力与前面第N个应力之间的差值绝对值,以第K个应力为底数、对应差值绝对值为幂指数得到的结果作为第K个应力的第二误差因子,组成一维的第二误差因子序列;
计算所述pH值序列中第K个pH值与前面第N个pH值之间的差值绝对值,以第K个pH值为底数、对应差值绝对值为幂指数得到的结果作为第K个pH值的第三误差因子,组成一维的第三误差因子序列;
将所述电阻序列中的第K个电阻所对应的前N个电阻作为一个子集合,获取子集合中的最大电阻和最小电阻,分别计算第K个电阻与最大电阻之间的第一电阻差值绝对值、第K个电阻与最小电阻之间的第二电阻差值绝对值,以第二电阻差值绝对值为分子、第一电阻差值绝对值为分母得到对应的比值,将该比值作为第K个采样时刻的测量质量下降因子,组成一维的测量质量下降因子序列;其中第K个电阻是在第K个采样时刻采集的;
将所述第一误差因子、所述第二误差因子、所述第三误差因子和所述测量质量下降因子序列组成所述成分检测过程中的4维特征矩阵。
3.如权利要求2所述的一种畜牧粪便成分检测与识别方法,其特征在于,所述基于明氏距离计算任意两个所述特征矩阵之间的相似程度的方法,包括:
分别计算两个所述特征矩阵中每个维度下的明氏距离,结合所有维度的明氏距离得到所述相似程度。
4.如权利要求1所述的一种畜牧粪便成分检测与识别方法,其特征在于,所述测量设备误差状态指数的获取方法,包括:
基于所述成分检测过程的先后顺序,分别获取当前匹配对所对应的两个所述成分检测过程相邻上一个所述成分检测过程的所述特征矩阵,将两个相邻的上一个所述成分检测过程的所述特征矩阵之间的相似程度作为当前匹配对的所述测量设备误差状态指数。
5.如权利要求1所述的一种畜牧粪便成分检测与识别方法,其特征在于,所述利用预测值对畜牧粪便中检测的实时TK含量进行优化补偿,得到实际TK含量的方法,包括:
获取实时TK含量与预测值之间的第一乘积;基于历史TK含量获取历史TK含量的中值,获取历史TK含量的中值与第一优化值之间的第二乘积,第一优化值与预测值之和为1,将第一乘积和第二乘积之间的均值作为实时TK含量。
CN202211616295.6A 2022-12-16 2022-12-16 一种畜牧粪便成分检测与识别方法 Active CN115600117B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211616295.6A CN115600117B (zh) 2022-12-16 2022-12-16 一种畜牧粪便成分检测与识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211616295.6A CN115600117B (zh) 2022-12-16 2022-12-16 一种畜牧粪便成分检测与识别方法

Publications (2)

Publication Number Publication Date
CN115600117A true CN115600117A (zh) 2023-01-13
CN115600117B CN115600117B (zh) 2023-06-20

Family

ID=84854286

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211616295.6A Active CN115600117B (zh) 2022-12-16 2022-12-16 一种畜牧粪便成分检测与识别方法

Country Status (1)

Country Link
CN (1) CN115600117B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101975800A (zh) * 2010-09-06 2011-02-16 中国农业大学 畜禽粪便肥料成分检测方法及其检测装置
US20150342224A1 (en) * 2008-04-30 2015-12-03 Xyleco, Inc. Processing biomass
CN110570914A (zh) * 2019-08-19 2019-12-13 甘肃方正节能科技服务有限公司 测量沼液营养成分含量的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150342224A1 (en) * 2008-04-30 2015-12-03 Xyleco, Inc. Processing biomass
CN101975800A (zh) * 2010-09-06 2011-02-16 中国农业大学 畜禽粪便肥料成分检测方法及其检测装置
CN110570914A (zh) * 2019-08-19 2019-12-13 甘肃方正节能科技服务有限公司 测量沼液营养成分含量的方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
LEI ZHANG 等: "Metagenomic sequencing analysis of bacterial biodiversity in the treatment ofsurplus sludge and chicken manure mixtures by black soldier fly larvae(Hermetia illucens)", 《RESEARCH SQUARE》, pages 1 - 8 *
NISHA PILLAI 等: "Artificial Intelligence Models for Zoonotic Pathogens: A Survey", 《MICROORGANISMS》, pages 1 - 20 *
何雯菁 等: "畜禽粪便中有害物质分析方法及控制措施研究进展", 《畜牧业环境》, pages 30 - 33 *
姜山 等: "烟台地区主要畜禽粪便养分含量分析及其植物毒性评价", 《黑龙江畜牧兽医》, pages 49 - 53 *
薛超友 等: "沼液中乙酸含量比色测定法研究", 《农业资源与环境学报》, vol. 35, no. 5, pages 476 - 482 *

Also Published As

Publication number Publication date
CN115600117B (zh) 2023-06-20

Similar Documents

Publication Publication Date Title
Mottram et al. A wireless telemetric method of monitoring clinical acidosis in dairy cows
US7676346B2 (en) System and a method for observing and predicting a physiological state of an animal
Bonfatti et al. Genetic parameters of measures and population-wide infrared predictions of 92 traits describing the fine composition and technological properties of milk in Italian Simmental cattle
Brunel et al. Estimating age at maturation and energy-based life-history traits from individual growth trajectories with nonlinear mixed-effects models
CN114707786A (zh) 一种基于共现性网络的河流生态系统健康评价方法
van Breukelen et al. Genetic parameters for repeatedly recorded enteric methane concentrations of dairy cows
CN117828480A (zh) 一种基于深度学习的肉牛发情状态识别与评估方法
CN115600117A (zh) 一种畜牧粪便成分检测与识别方法
EP1540336B1 (en) A system and a method for observing and predicting physiological state of an animal
CN115812626A (zh) 一种种猪性能测定系统
Vaughan et al. Spring nitrogen fertilizer recommendation models for dryland hard red winter wheat
Hammes et al. Using 13C in cattle hair to trace back the maize level in the feeding regime—a field test
CN108225617A (zh) 传感器标定装置及标定方法
CN112560336A (zh) 畜牧物联网采集体温校准方法
CN111223002A (zh) 一种玉米区域干物质产量或青储产量评估方法和系统
Steibel Phenomics in Animal Breeding
Cai et al. Performance indices for tests of soil nutrient status: Extractable phosphorus
CN111238613B (zh) 一种鹅的智能动态无应激称重方法
Dakhlan et al. Use of nonlinear regression in predicting body weight of female Saburai goat
US20240003863A1 (en) Methods for measuring organic carbon content or agricultural soils
Ryan et al. Exploring definitions of daily enteric methane emission phenotypes for genetic evaluations using a population of indoor-fed multi-breed growing cattle with feed intake data
Gholipour et al. Genetic correlation estimates between direct and maternal genetic effects and estimates of the proportion of maternal permanent environmental variance to phenotypic variance for some important traits in Iranian fowls: A meta-analysis
Hong et al. Portable soil ph sensor using isfet electrode
Gholipour et al. Journal of Poultry Sciences and Avian Diseases
Aruna et al. Bayesian approach for evaluation of lactation curve in cross bred cattle based on bimonthly test day milk yield

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