会员体验
专利管家(专利管理)
工作空间(专利管理)
风险监控(情报监控)
数据分析(专利分析)
侵权分析(诉讼无效)
联系我们
交流群
官方交流:
QQ群: 891211   
微信请扫码    >>>
现在联系顾问~
首页 / 专利库 / 反刍动物 / 反刍动物监测系统

反刍动物监测系统

申请号 CN202080029974.6 申请日 2020-05-04 公开(公告)号 CN113727646A 公开(公告)日 2021-11-30
申请人 莱利专利股份有限公司; 发明人 宋翔宇;
摘要 提供了一种用于自动监测反刍动物的系统和方法。所述系统包括3D相机系统,所述3D相机系统从感兴趣区域、特别是腰旁窝获得图像。图像处理器将所述感兴趣区域中的表面曲率确定为时间的函数。基于所述函数达到局部极大值的频率,生成所述动物的健康指示。
权利要求

1.一种用于自动监测反刍动物、特别是比如奶牛等乳畜的系统,所述系统包括:‑3D相机系统,用于在至少预定时间段期间的连续时间点处获得所述动物的至少一个感兴趣区域的多个3D图像,‑连接到所述3D相机系统的控制设备,所述控制设备设置有‑用于处理所获得的多个3D图像的图像处理设备,以及‑输出设备,

所述控制设备被布置为基于经处理的图像确定健康指示,并将所述健康指示输出到所述输出设备,

其中,所述图像处理设备被布置为‑确定所述多个3D图像中的感兴趣区域,以及以下之一:

‑针对每个3D图像计算所述感兴趣区域的曲率值,并将计算出的曲率确定为时间的函数,或者

‑针对每个3D图像测量所述感兴趣区域中的至少一个预定点相对于所述动物的相对位置,并将所述相对位置确定为时间的函数,其中,所述控制设备被布置为确定:‑所述时间函数的局部极值出现的时间点,以及‑所述健康指示,这通过相对于预定标准分析所述时间点来进行。

2.根据权利要求1所述的系统,其中,所述3D相机系统包括飞行时间相机或结构光相机。

3.根据任一项前述权利要求所述的系统,其中,所述感兴趣区域是左侧腰旁窝、特别是所述左侧腰旁窝内的部分区域。

4.根据任一项前述权利要求所述的系统,其中,所述曲率值包括或者是所述感兴趣区域的曲率值的平均值、特别是在跨所述感兴趣区域的多个点上求平均得到的平均值。

5.根据任一项前述权利要求所述的系统,其中,所述控制设备被布置为确定所述局部极值在时间上的频率,并且通过将所确定的频率与预定频率标准进行比较来分析所述时间点。

6.根据权利要求5所述的系统,其中,所述控制设备被布置为通过滤除所述计算出的曲率的时间函数中的具有预定频率范围之外的频率的时间变化来过滤所述函数,特别是所述控制设备被布置为:

‑对所述计算出的曲率函数执行傅立叶变换以构建经变换的频率函数,‑去除所述经变换的频率函数的在所述预定频率范围之外的所有部分,以获得干净的频率函数,以及

‑通过分析所述干净的频率函数、特别是通过确定具有最强信号值的频率值或所述干净的频率函数内的频率信号的平均值来确定收缩频率。

7.根据权利要求6所述的系统,其中,所述预定频率范围是或包括出自每分钟0.5次至4次的频率,包含两端值。

8.根据权利要求5至7中任一项所述的系统,其中,所述频率标准包括当所确定的频率低于预定的频率阈值时生成健康警告。

9.根据权利要求8所述的系统,其中,所述频率阈值是所述动物的历史值、特别是所述动物的历史平均值,

或者其中,所述频率阈值是所述动物的文献值、特别是分别取决于所述动物的品种、年龄、泌乳天数、或饲料类型或饲养方案中的一个或多个的文献值。

10.根据权利要求8或9所述的系统,其中,所述频率阈值取决于所述动物进行的活动类型,例如饲喂、挤奶或休息。

11.根据任一项前述权利要求所述的系统,包括至少一个站,每个站选自由以下各项组成的组:挤奶站、饲养站、饮水站、治疗站、隔离站和选择站,其中,所述3D相机系统包括设置在所述站中的至少一个站中的3D相机。

12.根据权利要求1所述的系统,其中,所述控制设备被布置为基于所确定的时间点来确定所述局部极值之间的时间间隔,并且通过将所确定的时间间隔与预定时间间隔标准进行比较来分析所述时间点。

13.一种用于确定反刍动物、特别是奶牛的健康指示的方法,所述方法使用根据前述权利要求中任一项所述的系统并且包括以下步骤:‑在至少预定时间段期间的连续时间点处获得所述动物的至少一个感兴趣区域的多个

3D图像,

‑处理所获得的多个3D图像,‑确定所述多个3D图像中的感兴趣区域,以下之一:

‑针对每个3D图像计算所述感兴趣区域的曲率值,并将计算出的曲率确定为时间的函数,或者

‑针对每个3D图像测量所述感兴趣区域中的至少预定点相对于所述动物的相对位置,并将所述相对位置确定为时间的函数,以及确定:

‑所述时间函数的局部极值出现的时间点,以及‑所述健康指示,这通过相对于预定标准分析所述时间点来进行。

14.一种用于监测反刍动物的健康的方法,所述方法包括:‑通过根据权利要求1至12中任一项所述的系统来自动确定瘤网胃收缩频率和瘤胃填充,以及

基于所确定的瘤网胃收缩频率和所确定的瘤胃填充分别与预定的频率值、预定的瘤胃填充值的比较,生成健康警告。

说明书全文

反刍动物监测系统

[0001] 本发明涉及一种用于自动监测反刍动物、特别是比如奶牛等乳畜的系统,所述系统包括:3D相机系统,用于在至少预定时间段期间的连续时间点处获得所述动物的至少一
个感兴趣区域的多个3D图像;连接到所述3D相机系统的控制设备,所述控制设备设置有用
于处理所获得的多个3D图像的图像处理设备;以及输出设备。
[0002] 这样的动物监测系统是已知的。例如,文件WO 12/138290 A1披露了一种用于确定乳畜瘤胃的肠道填充水平的系统。3D相机系统确定瘤胃三角形的深度或体积。还已知基于
3D相机的身体状况评分系统。
[0003] 然而,肠道填充水平或身体状况评分对于反刍动物的健康状态的重要性有限。
[0004] 因此,需要实用、简单且可靠的系统和方法来基于动物的3D图像确定更多和/或更好的关于反刍动物的健康指示信息。
[0005] 本发明至少部分地实现了该目标,并且在一个方面提供了一种用于自动监测反刍动物、特别是比如奶牛等乳畜的系统,所述系统包括:3D相机系统,用于在至少预定时间段
期间的连续时间点处获得所述动物的至少一个感兴趣区域的多个3D图像;连接到所述3D相
机系统的控制设备,所述控制设备设置有用于处理所获得的多个3D图像的图像处理设备;
以及输出设备,所述控制设备被布置为基于经处理的图像确定健康指示,并将所述健康指
示输出到所述输出设备,其中,所述图像处理设备被布置为:确定所述多个3D图像中的感兴
趣区域,以及以下之一:针对每个3D图像计算所述感兴趣区域的曲率值,并将计算出的曲率
确定为时间的函数,或者针对每个3D图像测量所述感兴趣区域中的至少预定点相对于所述
动物的相对位置,并将所述相对位置确定为时间的函数,其中,所述控制设备被布置为确定
所述时间函数的局部极值出现的时间点以及所述健康指示(这通过相对于预定标准分析所
述时间点来进行)。
[0006] 本文使用以下见解:已知瘤网胃运动性可以提供有价值的健康信息。瘤网胃通过周期性收缩和放松来执行消化功能。在本文中,最初的收缩从网胃开始并通过瘤胃。这些收
缩混合和循环食糜。这些收缩可以通过3D相机观察到。这里要注意,仅确定曲率的单个值
(该值比如可以从根据单个图像确定瘤胃区域的深度或体积或身体状况评分的现有技术系
统导出)无法产生有关收缩的任何信息。
[0007] 注意,瘤胃运动性往往由兽医通过听诊(听身体声音)和颤动(通过感觉身体来诊察)来确定。这需要兽医在场,对于连续监测而言这是不切实际的。因此,在实践中,仅在怀
疑存在某种健康问题时才使用瘤胃运动性。因此,由此排除了在很早的阶段确定健康问题
或其前兆的可能性。另外,实际测试表明,人类观察填充满的瘤胃的收缩比使用根据本发明
的系统所能实现的更困难,这因此提供了可靠性方面的优势。由于是人工评估,因此兽医的
运动性确定也可能产生错误和主观性。而且这也会给动物或畜群带来压力或其他扰动。这
本身可能是导致真实瘤网胃运动性偏差的原因,使得兽医可能由于其自身的评估引起的压
力等而产生错误评估。由于这些原因,期望具有如本发明所提供的用于确定健康指示的非
接触式系统。
[0008] 在该方法中,可以将函数的局部极值确定为特别是计算出的曲率值或相对位置在时间上的局部极大值。这里的位置极大值被认为是最接近3D相机的位置。“局部”在数学意
义上被认为是“相对于变量”,不一定与空间位置有关。
[0009] 图像是在连续时间点处获得的。这些时间点不需要以规律的时间间隔分布,但它们必需是已知的以便允许对图像进行时间分析。此外,还应该在至少预定时间段期间获得
图像,同样以便允许进行合理的时间分析。该时间范围取决于预期的重现间隔,该重现间隔
可能取决于动物的类型等。在实践中,如果图像是在至少2分钟期间、并且更优选地在至少5
分钟期间获得的,则是优选的。至少对于奶牛,这允许观察足够数量的收缩,以可靠地确定
健康指示值。
[0010] 3D相机可以是固定的,也可以是可移动的,只要该相机被布置为观察感兴趣区域即可。此外,考虑到动物可能会在图像获取期间移动,该3D相机可能设置有具有这样的视野
的镜头,使得对于动物相对于3D相机的所有可能位置,该感兴趣区域都将包括在图像中。特
别地,3D相机可以被布置为可通过连接到控制设备的所提供的移动器控制装置而指向感兴
趣区域。然后,控制设备可能能够例如基于图像处理主动地将3D相机指向感兴趣区域。然后
可以使用更小的视野,这进而允许更高的准确度。
[0011] 实际上,确定多个3D图像中的感兴趣区域可以有利地受益于某种形式的跟踪,尤其是当动物移动时。例如,图像处理设备被布置为识别图像中的特定点,并使用该点来跟踪
和定位连续图像以进行比较。基于这样的跟踪,相机的移动器控制装置可以被布置为将3D
相机定位成使得所识别的点总是在图像中的相同位置。
[0012] 其他特殊的和本发明的实施例在从属权利要求中以及在说明书的当前以下部分中进行了描述。
[0013] 在实施例中,3D相机系统包括飞行时间相机或结构光相机。这种类型的3D相机提供高达30Hz或更高的相对较高的帧率。这确保不会错过图像之间感兴趣区域的收缩或其他
运动,同时在所获得的3D图像中可以很容易地考虑动物的运动。然而,不一定需要具有如此
高的帧率,这样的高帧率会导致计算数量相对较多。比如立体相机或激光扫描仪等其他3D
相机类型也是可能的。优选的3D相机的示例是“SwissRanger”3D相机或等效物、以及
“Kinect”型相机设备。在所有这些中,相机的提供3D信息的特性在本发明中被用于提供动
物的感兴趣区域的曲率数据。
[0014] 在实施例中,感兴趣区域是左侧腰旁窝,也称为瘤胃三角形。这是反刍动物左侧的或多或少的三角形区域,其提供了相对清晰的瘤网胃视野。它可以很容易被自动识别,例如
通过使用三角形并匹配相对较高局部曲率的区域来进行模板匹配。以下网址中描述了类似
的方法:https://nl.mathworks.com/help/images/ref/normxcorr2.html,并且参考了此
来源以获得数学细节。特别地,感兴趣区域是所述左侧腰旁窝内的部分区域,即,不是整个
瘤胃三角形而仅是其较小或较大的部分。已经发现,也可以通过查看腰旁窝的一部分来收
集有用的信息,这限制了确定和计算的数量。
[0015] 在实施例中,所述曲率值包括或者是所述感兴趣区域的曲率值的平均值、特别是在跨所述感兴趣区域的多个点上求平均得到的平均值。该曲率值可以根据任何期望的方法
来确定,只要它是表面弯曲程度的指示即可。根据该实施例,可以通过分析完整的或部分感
兴趣区域来确定该值。也可以分析所述完整的或部分区域中的多个点,同样只要获得弯曲
程度的良好指示即可。另外需要注意的是,曲率值的符号不能忽略,因为凸感兴趣区域应该
给予正值,而凹感兴趣区域应该给予负值,反之亦然,只要这两者是不同的即可。在感兴趣
区域的外观在凸与凹之间变化的情况下,这可以防止无意义的双峰值。
[0016] 计算表面(感兴趣区域)曲率的有用示例是,通过对与表面上点的邻居的协方差矩阵进行主成分分析。这给出了矩阵的三个特征值λ0、λ1、λ2,点的曲率σ可以通过σ=λ0/(λ0+
λ1+λ2)来估计,其中,λ0是最小的特征值。然后,可以通过对感兴趣区域中的多个点(比如点
矩阵)的表面曲率求平均(或相加等)来确定表面曲率。结果证明,在感兴趣区域中的某个固
定且易于识别的点周围的表面点网格或栅格就足够了,如上文已经说明的。该固定点可以
是例如用于跟踪的具有最高互相关系数的表面点。然而,这并不排除其他估计表面曲率的
方法,比如通过计算感兴趣区域或其部分/窗口中的所有表面法向量的变化。
[0017] 在确定表面曲率(可以针对每个图像完成此步骤)之后,结果是作为时间的函数的曲率。在分析该函数时,确定局部极值、特别是局部极大值。在本文中,如果在某个值的时间
之前和之后的预定时间窗口(比如几秒,例如3s)中该值是最高的,则将该值视为局部极值
是有利的。这可以防止噪声尖峰等过多地干扰。
[0018] 在实施例中,所述控制设备被布置为确定所述局部极值、特别是局部极大值在时间上的频率,并且通过将所确定的频率与预定频率标准进行比较来分析所述时间点。极值、
特别是极大值将以或多或少的规律的频率发生。通过确定该频率或频带,可以获得关于动
物健康的有价值的信息,并给出健康指示。
[0019] 在实施例中,所述控制设备被布置为通过滤除所述计算出的曲率的时间函数中的具有预定频率范围之外的频率的时间变化来过滤所述函数。通过这种措施,可以有效地去
除不相关的变化,比如由动物运动/位移或噪声引起的变化。本文使用以下见解:相关收缩
发生在动物特定的频带内。特别地但非限制性实施例,所述控制设备被布置为对所述计算
出的曲率函数执行傅立叶变换以构建经变换的频率函数,然后去除所述经变换的频率函数
的在所述预定频率范围之外的所有部分,以获得干净的频率函数,以及通过分析所述干净
的频率函数来确定收缩频率。傅里叶变换是众所周知的从时间信号中去除频率分量的有效
方法。然而,可以应用实现类似效果的其他滤波方式,比如带通滤波器。对干净的频率函数
的分析可以通过确定具有最强信号值的频率值或干净的频率函数中剩余信号内的频率信
号的平均值来完成,或者通过给出指示收缩频率值的任何其他分析来完成。另外但仅可选
地,控制设备可以被布置为对干净的频率函数执行逆傅立叶变换,以获得干净的曲率函数。
可以显示该干净的曲率函数以供视觉检查。
[0020] 特别地,所述预定频率范围是或包括出自每分钟0.5次至4次的频率,包含两端值。对于奶牛,其被证实是有用的频带,在该频带之外没有表示有意义的收缩的信号。对于其他
物种,不同的频带可能是适用的。
[0021] 在实施例中,所述频率标准包括当所确定的频率低于预定的频率阈值时生成健康警告。已经发现,对于正常动物,瘤网胃收缩发生的频率具有至少一定的值,并且还发现,对
于患有临床甚至亚临床疾病的动物,该频率会降低。这允许确定健康警告,比如注意列表上
的条目、给农场主或兽医的消息、可听或可视警报等。所有这些都用于引起对特定动物的注
意以进行进一步的健康检查,因为发现所述动物的某个值不同寻常。
[0022] 在实施例中,所述频率阈值是所述动物的历史值、特别是所述动物的历史平均值,或者所述频率阈值是所述动物的文献值、特别是分别取决于所述动物的品种、年龄、泌乳天
数、或饲料类型或饲养方案中的一个或多个的文献值。这允许根据动物的具体情况微调频
率阈值,以便防止不必要的健康警告。例如,当动物主要由粗饲料饲养时,收缩频率通常高
于当由富含精饲料的膳食饲养时的收缩频率,于是可以相应地调整阈值。这里要注意的是,
文献值可以是起始点,而可选地通过深度学习等方式对历史值进行处理可以产生具有更高
可靠性的阈值。
[0023] 特别地,反刍动物是反刍动物群中的一个,并且所述频率阈值取决于所述反刍动物群的平均频率,更特别地,其中所述平均值是瞬时平均值或移动平均值。这确保可以考虑
适用于整个畜群的变化,而不会引起警报(或健康警告)。如果在比如冬季时在室内停留和
饲养数月后,允许畜群在特定的一天外出吃草,则可能会出现全畜群的瘤网胃频率变化,这
是正常的。因此,为所有动物发出警报是无意义的。但是,如果对于某头动物而言这种变化
明显不同,则仍然可能会发出警报。
[0024] 更一般地,认为实际上在所有情况下,都应该为相对于畜群平均值的变化发出警报,至少要调查这种偏差的可能原因。在许多情况下,与畜群平均频率的偏差是对于疾病、
比如特别是代谢性疾病的良好指示值。
[0025] 例如,频率阈值的值可以低于或低于畜群平均频率某一绝对值或相对值,比如高于或低于畜群平均频率值0.5Hz或15%。但是取决于动物的类型、畜群的大小、期望的准确
度等,当然可以选择其他数量。
[0026] 在实施例中,所述频率阈值取决于所述动物进行的活动类型,例如饲喂、挤奶或休息。本文使用以下见解:收缩频率取决于动物短时间内参与的活动类型。例如,如果动物正
在进食,则该频率通常会较高,即,在预期频带内的高侧。相反,在压力情况下或休息时,该
频率会较低。
[0027] 在实施例中,系统被布置为基于所获得的多个图像或者至少其中的绝对或相对变化来确定瘤胃填充水平。基本上,可以应用与现有技术文件WO 12/138290A1中相同的方法,
该文件中的描述这种方法的部分被视为并入。简而言之,该文件中披露的方法和系统确定
瘤胃三角形的平均凹度或凸度,并由此确定瘤胃填充值。注意,在本发明的当前实施例中,
可以并且有利的是使用多于一个图像,以便确定曲率的平均值并因此确定瘤胃填充的平均
值。这确保了正确地确定真正的瘤胃填充值,而不受或几乎不受突然运动、呼吸等的影响。
[0028] 通过这样确定瘤胃填充值,系统能够甚至更准确地确定何时可以给出健康警告或警报。特别地,控制设备被布置为在发生以下情况之一时发出健康警告:
[0029] ‑所确定的频率在时间上降低,所确定的瘤胃填充值在时间上增加;
[0030] ‑在制度改变后,所确定的频率与畜群平均频率值相比变化更大;
[0031] ‑在制度改变后,动物的所确定的频率没有恢复到在饲养制度改变之前的该动物的频率;或者
[0032] ‑在制度改变后,动物的所确定的频率确实恢复到其原始值,但其恢复速度慢于畜群平均值至少预定量。
[0033] 这四种可能性中的第一种是瘤胃运动性受抑制的迹象,这是代谢性疾病、特别是瘤胃酸中毒的迹象,应予治疗。
[0034] 第二到第四种可能性与对代谢可塑性的量化和对制度变化的适应有关。这种制度变化特别是饲养制度的变化,但也可能是(畜棚)气候、畜群密度或任何农场管理实践的变
化。在每种情况下,都建议更密切地监测,因为所指示的变化是代谢问题的迹象。在本文中,
应当注意,迄今为止,已经通过血液值和/或牛奶的实验室测试确定或量化了这种代谢可塑
性。本发明人已经发现,通过本系统和方法可获得类似有效的健康警告系统,该系统是非侵
入性的,并且因此对于动物而言不适感较小。
[0035] 也可以使用所确定的频率以及可选的瘤胃填充值作为对其他指示值的检查。例如,经常使用产奶量作为健康警告标准。然而,产奶量下降本身并不指示生病,并且健康警
告仅基于产奶量可能会过于经常导致误报的健康警告或出现过晚的警告,并导致健康或产
量永久性变差。通过本发明,所确定的瘤网胃收缩频率过高或过低(无论是在绝对意义上还
是相对于群体平均值)、或者所确定的频率在瘤胃填充值上升时下降、根据前文提及的其他
组合的等等可能不仅本身会引起健康警告,而且实际上也可能反过来,即,如果产奶量下降
并且正常地引起健康警告,而所确定的频率保持稳定,则控制单元会否决该健康警告。于
是,控制单元被布置为相应地给出或不给出健康警告。类似地,如果瘤胃填充和/或乳脂浓
度降低,但动物的所确定的瘤网胃收缩频率保持稳定,则无需健康警告,并且控制单元也应
相应地布置。然而,可能会给出饲养警告,因为后一种情况与饲料(摄入量)不足有关,而不
是与动物的(急性)疾病有关。在所有情况下,系统应当包括本领域已知的用于确定比如产
奶量或乳脂浓度等附加量的任何装置。
[0036] 综上所述,本发明提供了一种用于监测反刍动物的健康的有利方法,所述方法包括:通过根据本发明的系统自动确定瘤网胃收缩频率和瘤胃填充,以及基于所确定的瘤网
胃收缩频率和瘤胃填充分别与预定的频率值、预定的瘤胃填充值的比较,生成健康警告。在
本文中,相应的预定值可以是时间平均值或畜群平均值、文献值等。这种方法的一大优点
是,使健康警告基于这些参数可以加快对严重问题的检测,因为当存在与代谢相关的问题
时,瘤胃收缩和瘤胃填充首先发生变化,然后才是身体状况评分、产奶量和奶品质的变化。
现在将继续描述进一步的系统实施例。
[0037] 在实施例中,该系统包括至少一个站,每个站选自由以下各项组成的组:挤奶站、饲养站、饮水站、治疗站、隔离站和选择站,其中,3D相机系统包括设置在所述站中的至少一
个站中的3D相机。这些站非常适合放置3D相机系统,因为动物将在一个位置停留相对较长
时间。特别是挤奶站非常适合,因为甚至可以基于预期的产奶量和其他历史数据预测动物
将在此停留多长时间。平均而言,动物将在挤奶站停留至少5分钟,这为如上文已经陈述地
那样确定频率提供了非常合适的时间窗口。此外,在其他站也可以使用3D相机系统进行此
类测量。即使不是每个实例都产生可靠的测量值(例如因为动物停留在合适位置的时间过
短),本发明的系统也允许在许多实例中监测动物,使得可以获得对该动物的健康的良好见
解,并且可以在出现问题时尽早地发出警告。这里要注意的是,在每个站,除了测量瘤网胃
收缩之外,还可以执行一些其他动作,比如挤奶、饲养、治疗或隔离动物。有利地,隔离站联
接到挤奶站或饲养站,使得在针对站中的特定动物给出健康警告的情况下,控制设备可以
将动物送到隔离站以等待由农场主或兽医进行健康检查。
[0038] 在实施例中,控制设备被布置为基于所确定的时间点来确定局部极值之间的时间间隔,并且通过将所确定的时间间隔与预定时间间隔标准进行比较来分析时间点。在上文
中,系统使用基于频率的分析来确定健康指示或警告。然而,也可以基于局部极值之间的时
间间隔进行分析。实际上,在经历从时域到频域的变换之后,这种分析将归结为与基于频率
的分析相同。然而,在某些情况下,仅分析各种时间间隔可能更简单,比如通过确定局部极
大值之间的平均时间间隔。于是,这样的平均时间间隔应该在一定时间间隔范围内,该时间
间隔范围可能取决于动物、活动等。与频率相关分析系统中的特殊实施例有关的所有特征
也适用于基于时间间隔的分析系统。
[0039] 在第二方面,本发明涉及一种用于确定反刍动物、特别是奶牛的健康指示的方法,所述方法使用根据本发明第一方面的系统并且包括以下步骤:在至少预定时间段期间的连
续时间点处获得所述动物的至少一个感兴趣区域的多个3D图像,处理所获得的多个3D图
像,确定所述多个3D图像中的感兴趣区域,以下之一:针对每个3D图像计算所述感兴趣区域
的曲率值,并将计算出的曲率确定为时间的函数,或者针对每个3D图像测量所述感兴趣区
域中的至少预定点相对于所述动物的相对位置,并将所述相对位置确定为时间的函数,以
及确定时间函数的局部极值出现的时间点和健康指示(这通过相对于预定标准分析所述时
间点来进行),和/或确定所述函数的幅度和健康指示(这通过相对于预定幅度标准分析所
确定的幅度来进行)。由于这是与本发明的系统方面对应的方法,因此这里只须说明针对与
本发明的系统方面的特殊实施例有关的措施提及的所有特殊特征和优点也适用于方法方
面就足够了。
[0040] 这里要强调的是,本方法不是诊断方法,尽管其可以用于诊断方法中。毕竟,该方法仅有助于确定动物是否有问题。其无法做出关于给出健康警告的动物存在何种疾病(如
果有的话)的诊断。
[0041] 现在将参考一个或多个示例性且非限制性实施例以及附图来阐明本发明,在附图中:
[0042] ‑图1非常示意性地示出了根据本发明的系统1,以及
[0043] ‑图2示出了原始曲率与时间的示例性曲线图、以及经平滑的曲率对时间的函数。
[0044] 图1非常示意性地示出了根据本发明的系统1,该系统用于确定奶牛2的健康指示,并且包括3D相机3、具有图像处理设备5的控制设备4、具有乳杯7的挤奶机器人6和具有传感
器9的饲料槽8。呈左侧腰旁窝形式的感兴趣区域用附图标记10指示,并且子区域或窗口用
附图标记11指示。
[0045] 在所示的实施例中,非常示意性地示出了通过挤奶机器人6利用乳杯7(这里示出了其中一个乳杯)给奶牛2挤奶的挤奶站。一旦挤奶过程开始,即在识别奶牛并决定她将被
挤奶之后,控制设备4将能够粗略估计该奶牛将在挤奶站花费的时间。为此,其可以使用标
准最少时间、该奶牛的历史挤奶时间、或甚至是基于生产和挤奶间隔的估计挤奶时间,如本
领域本身已知的。对于几乎每次挤奶,此时间至少为5分钟,通常长达8或9分钟。在奶牛不被
挤奶的情况下,她会被驱赶到外面,并且没有足够的时间进行有意义的测量。
[0046] 可替代地,该站是饲养站,通过饲料槽8示意性地指示,该饲料槽具有传感器9,该传感器通过奶牛鼻子的按压而指示奶牛开始进食。大多数挤奶站也会有饲料槽,但饲料槽
也可能是独立的系统。另外的其他替代性方案可以是饮水站(或饮水槽)、治疗站等。而且,
大多数站都具有动物标识(这里未示出),通过该标识,可以识别动物并且可以单独调整设
置(挤奶、饲养、治疗)。此外,还可以通过这样的动物识别来触发3D相机系统,例如,在系统
应仅监测特定动物的情况下。
[0047] 当3D相机被打开时,其开始对动物的感兴趣区域(下文中:ROI)进行成像,这里感兴趣区域指示为左侧腰旁窝(由虚线指示)。为了能够获得足够可靠的数据,帧速率至少为
每秒一帧,但优选地至少为其十倍。为了确保在成像过程中ROI在相机的视野内,可以在相
机上提供广角镜头,使得当动物在(挤奶、饲养、...)站处或中时,ROI将在视野内。还可以提
供用于基于图像处理设备5对图像中的ROI的识别来移动相机的马达。以这种方式跟踪ROI
本身是已知技术。其优点是ROI可以形成图像的较大部分,并且因此可以以更高的分辨率成
像。
[0048] 所获得的3D图像由图像处理设备5处理,这将在下面进一步阐明。该处理的结果是作为时间的函数的ROI的曲率值。由控制设备5分析该函数,并且应用一个或多个标准来确
定奶牛2的健康指示。在健康指示引起警报等的情况下,控制设备4可以将奶牛2加入注意列
表、发出可听或可视警报、在站2、6之后隔离奶牛2等。然后,奶牛2将由农场主或兽医进一步
诊察。
[0049] 在使用该系统时和在该方法中,所获得的3D图像形成ROI的3D表示。为了限制计算的数量,将在下面描述,可以将ROI限制到左侧腰旁窝的仅一部分,比如限制到本示例中的
子区域或窗口11,尽管这并非是必要的。
[0050] 代替确定该ROI的深度或体积(如在现有技术系统中为确定瘤胃填充所做的那样),本发明确定ROI的变化曲率值。这基于以下见解:自然过程影响瘤胃形状,使得瞬时瘤
胃填充值的预测值似乎是受限的,但相反,对瘤胃或左侧腰旁窝区域的曲率的时间变化分
析的预测值似乎是有意义的。
[0051] 针对每个所获得的3D图像,通过图像处理(比如通过识别左侧腰旁窝的左上角)并重新定位图像/调整图像大小来跟踪感兴趣(子)区域。然后,计算曲率值。这可以通过多种
方式完成,只要其以系统的方式表达ROI的弯曲程度即可。现在将简要阐述一个示例。
[0052] 曲率的确定将限于子区域11。针对该子区域11,曲率值如下确定。
[0053] 首先,从点的最近邻居计算协方差矩阵。
[0054]
[0055] 其中,k是邻近点的数量,pi是第i个邻近点的位置向量,并且 是邻近点的质心的位置向量。所得到的协方差矩阵C将是具有3个特征值的3×3矩阵。表面曲率σ可以通过以下
等式估计
[0056]
[0057] 其中,λ0、λ1、λ2是协方差矩阵C的特征值,其中λ0是最小的特征值。
[0058] 然后,可以将所得到的曲率值σ绘制为时间的函数。这是在示例性的图2中完成的,作为变化有点大的曲线。注意,x轴表示图像编号,在这种情况下,帧速率为30Hz,或1800张
图像/分钟。显然,虽然在曲线中可辨别出粗略的“节拍”,但很难从中提取出有意义的信息。
然而,已经意识到可以有效地消除各种噪声源。一方面,可以去除动物运动以及变化的光照
条件,这些光照条件原则上是一次性变化而不是规律的变化。另外,可以消除比预期收缩频
率快得多的规律的运动,比如呼吸。后者通常在每分钟约25到50次呼吸之间,这比瘤网胃收
缩高一个数量级。在这种情况下,上述“噪声信号”借助通过信号频率分解信号的方式来去
除,其中,“通过”频率范围在每分钟0.5到4.0次收缩之间,而丢弃其余频率。为此,构建信号
的傅立叶变换,应用该通过频率范围以及低于和高于该范围的截止频率范围,并构建逆傅
立叶变换以重新生成曲率‑时间函数。该函数也绘制在图2中,作为平滑曲线。在平滑曲线
中,或多或少均匀间隔的峰值清晰可见,而且这些峰值的频率也清晰地在预期范围内。注
意,经平滑的曲率函数仅用于视觉检查,而经傅立叶变换的函数则用作计算和监测的依据。
[0059] 在经平滑的曲线中,峰值或极大值已用点指示。因此,对于瘤网胃收缩,平均频率约为每分钟9/(5800/1800)=2.8次收缩。这是被测奶牛的正常频率,在这种情况下是初产
的健康荷斯坦和瑞士棕奶牛,因此对于该特定奶牛,可以给出“健康”指示,并且不需要给出
健康警报。
[0060] 然而,可能的是,基于历史测量结果,某头特定奶牛通常具有更高或更低的值。在这种情况下,可以给出“仍然健康,但请检查”的健康指示,即,可以进行一些诊察,但不是很
紧急,或者可以更密切地监测该奶牛。还可能的是,计算出的收缩频率实际上低于预定阈
值,比如当在饲养期间所确定的频率为每分钟大约1到2之间时。这样的值在休息期间并不
少见,但在饲养期间通常应该更高。因此,在这种情况下,可以由控制设备发出“急需检查”
的健康指示,或警报发声等。
[0061] 在所示的示例中,在挤奶站中提供了3D相机系统。在该站内,每天会给动物挤奶几次,比如2‑4次/天。如果更频繁地监测动物是有利的,因为例如测量可能会由于时间太短而
无法得出有意义的信息、或者由于运动太剧烈等而失败。因此,如果3D相机系统在比如饲养
站、供水站、甚至隔间或“休息站”等其他位置包括一个或多个附加相机,则是有利的。注意,
可以将所确定的频率与对应地调整的阈值频率(比如在隔间中时的较低阈值频率)进行比
较。
[0062] 注意,3D相机生成3D图像,这些图像表示每个像素都结合了关于到相机的距离的信息的动物的2D图像。在上述实施例中,在感兴趣区域中计算曲率,并且基于对曲率的时间
分析得出结论。还可能的是,具有图像处理设备的控制设备计算感兴趣区域中的特定点与
相机之间的距离。该点也会以与曲率变化相同的频率朝向相机移动并再次移回。换句话说,
该点相对于相机的相对位置以相同的相关频率变化。因此,本发明——系统和方法两者也
在它们(被布置为)测量感兴趣区域(左侧腰旁窝或其子区域)中的固定点的相对位置时起
作用,并且分析极值、特别是固定点距离相机最近的时间点。
[0063] 有利的是,以足够的精度和准确度确定固定点,以便防止人为因素或简单的错误测量。此外有利的是,固定点在图像中可容易地识别。在本文中,确定腰旁窝的边界是有帮
助的,该腰旁窝是动物左侧可相对容易识别的三角形。然后,可以由图像处理软件确定固定
点位于如此确定的感兴趣区域的边界的相对位置处,比如在几何中心或任何其他位置。然
后,通过使用区域的边界的相对位置来确定动物的基本位置是相对简单的,即使在事后而
不是实时进行分析时也是如此(比如使用肋骨和脊椎,这另外地表明确定呼吸率是有利的,
以便减去对应的信号)。这些边界整体的运动算作动物整体的位移,并且对于收缩没有意
义。在减去这些、或以其他方式产生影响的信号之后,可以确定固定点的真实相对运动,并
且可以针对该相对运动或多或少地复制剩余的分析。
[0064] 上述实施例仅用于帮助说明本发明,而不以任何方式限制本发明。本发明的范围而是由所附权利要求确定。