CN112262408A - 用于从一距离处表征活体标本的方法和装置 - Google Patents

用于从一距离处表征活体标本的方法和装置 Download PDF

Info

Publication number
CN112262408A
CN112262408A CN201980028070.9A CN201980028070A CN112262408A CN 112262408 A CN112262408 A CN 112262408A CN 201980028070 A CN201980028070 A CN 201980028070A CN 112262408 A CN112262408 A CN 112262408A
Authority
CN
China
Prior art keywords
image
living specimen
living
specimen
map
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.)
Pending
Application number
CN201980028070.9A
Other languages
English (en)
Inventor
I·阿马特·罗尔丹
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.)
Non Contact Animal Measurement Co ltd
Original Assignee
Non Contact Animal Measurement Co ltd
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 Non Contact Animal Measurement Co ltd filed Critical Non Contact Animal Measurement Co ltd
Publication of CN112262408A publication Critical patent/CN112262408A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • G06T7/0014Biomedical image inspection using an image reference approach
    • AHUMAN NECESSITIES
    • A01AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01KANIMAL HUSBANDRY; CARE OF BIRDS, FISHES, INSECTS; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
    • A01K29/00Other apparatus for animal husbandry
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • G06T7/73Determining position or orientation of objects or cameras using feature-based methods
    • G06T7/74Determining position or orientation of objects or cameras using feature-based methods involving reference images or patches
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/44Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/103Static body considered as a whole, e.g. static pedestrian or occupant recognition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10024Color image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10048Infrared image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20036Morphological image processing
    • G06T2207/20041Distance transform
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20036Morphological image processing
    • G06T2207/20044Skeletonization; Medial axis transform
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20061Hough transform
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30168Image quality inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation
    • G06T2207/30188Vegetation; Agriculture

Abstract

本发明涉及一种用于从一距离处表征活体标本的方法,所述方法包括采集活体标本的图像和分割所述图像,从而提供分割图像;测量至所述图像的若干部分的距离,从而提供若干距离测量值,和选择所述活体标本的所述分割图像中所包含的那些距离测量值的子组;处理所述分割图像和所述若干距离测量值,所述若干距离测量值通过以下方式指代所述分割图像内所包含的不同位置:表征所述活体标本的形状,表征所述活体标本在所述分割图像所包含的所述距离测量值内的深度;和比较所述形状分析图和所述深度剖面分析图。如果所述比较的结果包括于给定范围之内,则进一步确定所述活体标本的参数,包括姿势参数,和/或所述解剖基准点的定位或校正,和/或身体尺寸参数,和/或表示所述活体标本的身体图。

Description

用于从一距离处表征活体标本的方法和装置
技术领域
本发明整体涉及自动测量方法和系统的领域。特别地,本发明涉及一种用于从一距离处(即,远程地或以非接触方式)表征活体标本(诸如牲畜动物)的方法和装置。表征包括活体标本的尺寸参数(包括取向、尺寸和姿势等)的计算,和/或活体标本的3D表示。
在本文档中,“图”应被理解为多种空间关系或一系列的特征,或其中不同信息相关联的(一维、二维或多维)图表。因此,图可为一系列的身体尺寸和取向,或身体温度在不同位置的关系。这特别地适用于形状分析图、深度剖面分析图和身体图。
背景技术
用于活体标本的远程表征的方法和/或装置在本领域是已知的。
例如,本发明的相同申请人的EP3158289涉及一种用于目标(诸如猪或其它牲畜动物)的自动参数计算的方法和装置。该方法包括:通过二维相机在场景中采集至少一个目标的二维图像;识别该采集二维图像内的目标;考虑到目标和二维相机之间的距离,通过第一装置计算目标在该采集和分割二维图像中的像素尺寸;以及利用目标的所述计算像素尺寸和先验模型,通过第二装置计算若干参数,该若干参数至少包括目标的尺寸、维度、身体部分维度、身体特征、重量和/或体积,其中所述先验模型包括将代表若干目标(200)的不同部分、轮廓或形状(先前预先使用二维相机采集)与所述若干目标的若干参数相联系的信息。
US-5474085提供了一种用于牲畜的远程感测的方法和设备(利用热成像图像感测系统),以确定监视区域中动物的数量、重量、位置、温度、畜体pH等。将包括区域像素的热成像图像发送至微型计算机中的数字化板,图像在其中将被转换成数字阵列。然后,数字通过软件进行解释以可判读形式提供期望的信息。
US-5412420公开了一种系统,该系统测量动物(诸如奶牛)的三维表型特性。该系统利用激光雷达相机的大量调制激光光束对动物的每平方英寸的大约100个点进行测量。每个激光束测量强度、水平、竖直和深度维度,并且通过组合这些测量值,该系统构成了关于该动物的极其准确的三维图像。该系统通过组合动物上所选点的测量值,计算关于该动物的形态的期望表型测量值。然后,该系统将每个动物的测量值存储于计算机数据库中以备使用。该系统还存储由动物标记组成的光强度图像,将该光强度图像与其它存储图像相比较。
US-A1-20150302241公开了用于改善对象在工业环境中的健康和福利的系统和方法。系统可包括相机和联接至相机的处理器,该相机布置成观察对象的一个或多个特征,该处理器配置成分析从相机获得的一个或多个图像,以从对象的一个或多个图像提取一个或多个特征,并对一个或多个特征进行分析,或替换嵌套其中的特征以预测对象状态的结果结果。特别地,该系统可配置成基于分析而生成诊断信号(例如,结果结果、发烧、乳腺炎、病毒、细菌感染、发情期等)。
文献“利用热成像和KINECT传感器的黑牛身体形状和温度测量(Black cattlebody shape and temperature measurement using thermography and KINECT sensor)”介绍了一种黑牛身体形状和温度测量系统。由于本文献的作者指出,重要的是在生长过程期间周期性地评估日本黑牛的质量,所以不仅是黑牛的重量和尺寸,并且姿势、形状和温度也需作为主要评估标准进行追踪。在该研究中,KINECT传感器和热相机获得了身体形状和其温度。整个系统进行校准以在共同坐标系中操作。点云数据从不同角度获得并且在计算机中进行重建。还捕获了热数据。通过考虑到肉牛的取向,将点云数据和热信息进行组合。所收集信息用于评估和估计牛的状况。
然而,这些现有技术文献均不允许执行快速(秒级以下方案)和自动测量来可靠、可重现和准确地估计活体标本的3D取向和/或姿势参数,和/或均不允许当活体标本在农场或在其自然环境中自由地移动时计算其身体图。
发明内容
根据第一方面,本发明提出了一种用于从一距离处表征活体标本(优选地,牲畜动物,诸如猪、公牛、母牛、羊、肉鸡、鸭或雏鸡等)的方法,同时该动物在农场或在其自然环境中自由地移动。应当注意,该方法适用于具有复杂形状的任何目标的表征。
该方法包括a)经由图像采集单元(诸如相机)采集活体标本的一个图像和通过处理单元进一步分割该采集图像,从而提供分割图像;b)通过遥测单元(在图像采集单元的给定距离处)测量至采集图像的若干部分的距离,从而提供若干距离测量值,和选择活体标本的分割图像中所包含的那些距离测量值的子组;和c)通过处理单元(等同于或不同于其它处理单元)处理分割图像和所述若干距离测量值,这些若干距离测量值指代分割图像内所包含的不同位置。
根据所提出方法,所述步骤c)包括表征活体标本的形状,估计活体标本的深度,和比较所述预先表征的结果,以获得质量参数/估计指示(活体标本的身体部分或解剖基准得以实际地测量并适当地定位,或需要被找出的更佳估计)。
即,如果比较的结果包括在给定范围之内,那么意味着所执行的测量是正确的,该方法还可确定活体标本的一些参数(例如,姿势参数(诸如活体标本的身体的深度和/或弯曲的取向)、解剖基准点的定位或校正、身体尺寸参数等等),和/或还可表示活体标本的身体图(优选为3D图)。相反地,如果比较的结果包括在所述给定范围之外,那么意味着所执行的测量是不正确的,例如因为活体标本在采集图像时移动。该方法还可包括使现有步骤a)至步骤c)重复,从而获得新深度剖面分析图和新形状分析图。可选地,如果结果包括于该范围之外,那么可选择什么都不做并且表示活体标本的身体图,该身体图将具有相关联错误。
优选地,形状的表征通过实施一种算法而执行,该算法在分割图像内至少计算以下项的一者或多者:活体标本的质心、活体标本在分割图像内相对于基准点的取向、和/或活体标本的特定身体部分(通过将活体标本的解剖基准点定位于分割图像内)。形状表征的结果提供了形状分析图。
深度的表征还优选地通过实施一种算法而执行,该算法至少在分割图像所包含的距离测量值内,通过将活体标本的解剖基准点定位于该距离测量值内来计算活体标本的特定身体部分。深度表征的结果提供了一个深度剖面分析图(其可提供一个以上)。
应当注意,执行用于表征的步骤的次序为不相关的。此外,两种表征可同时进行。
此外,根据所提出方法,图像采集单元(例如,RGB相机或热相机或该两种相机)和遥测单元(例如,激光雷达或飞行时间(TOF)系统)被校准。优选地,该两种单元布置于它们之间的给定距离处,并且特别地附接至共同支撑件。
在一个实施例中,该方法还估计图像采集单元和活体标本的相对位置的三维信息的一部分以获得一些额外参数,诸如:图像采集单元和活体标本之间的至少一个角度的平均值、活体标本的形状的弯曲程度或平面度、图像采集单元相对于地面的高度或图像采集单元相对于活体标本高度的高度、和/或图像采集单元的光学轴线相对于地面的角度。
在一个实施例中,活体标本的取向通过经由二乘拟合函数、高斯模型、主分量分析(PCA)、最小面积矩形、霍夫变换或与二维傅里叶变换的主轴线的相对性等将分割图像拟合成椭圆进行计算。
当身体部分在形状分析图中进行计算的情况下,该身体部分可通过圆形霍夫变换进行计算,该圆形霍夫变换计算包含分割图像内的大腿或股部的一部分的半径。另选地,身体部分可通过二阶多项式函数进行计算,该二阶多项式函数通过拟合质心和取向轴线周围的抛物线而检测分割图像内的活体标本的尾部。
另外,为改善形状分析图,可执行其它计算。例如,在一个实施例中,计算分割图像内的活体标本的轮廓,在极坐标上编码该计算轮廓并且还将傅里叶变换函数应用至所述极坐标,提供若干傅里叶系数,这些傅里叶系数的模量为旋转不变量并且其论证包含旋转信息。
在另一个实施例中,分割图像可编码为提供若干系数的图像矩,例如:统计矩、中心矩或Hu矩,这些系数为形似傅里叶变换的表示。然而,这种操作可适用于分割区域、轮廓或轮廓子组。
在另一个实施例中,计算活体标本的轮廓,并且距离度量还包括基于欧几里得、测地或城市街区等的距离度量而在该计算轮廓内进行计算。
在另一个实施例中,分割图像的活体标本的轮廓通过骨架化函数进行计算,从而提供活体标本的骨架的图像。任选地,还可计算所述骨架内的分支点和端点以估计不同身体部分的解剖位置。
在又一个实施例中,计算分割图像的距离变换。
步骤a)可包括活体标本的若干图像在不同时间段的采集,使得可捕获活体标本的不同姿势。在这种情况下,对于每个采集图像,获得了一系列距离测量值。
在后者中,关于每次采集所获得的信息可进行整合/组合,使得一系列成对深度剖面分析图和形状分析图得以获得。然后,该方法还可包括对每对图指定分数和选择具有最高分数的一对图。另选地,该方法还可匹配所有采集内的解剖基准点并且累积深度剖面分析图和解剖基准点的不同片段以计算活体标本的三维重建,或甚至该方法可计算每次采集的身体图并且累积每个身体图的所有信息,计分至扩展(或改进)身体图中。
在活体标本的身体图得以表示的情况下,该身体图可用于基于图像采集单元的特征而计算身体图的特性(例如,色彩、温度等),或计算通过预先校准的额外装置(例如,高分辨率热相机、光谱性质)所测量的额外特性。
可发生的是,采集图像包括一个以上的活体标本。在这种情况下,在一个实施例中,所提出方法可计算并比较图像中所包括的每个活体标本的形状分析图和深度剖面分析图,使得在单次采集中可表征一个图像所包括的所有标本。
根据另一个方面,本发明还提出一种用于从一距离处表征活体标本的装置。该装置包括图像采集单元、第一处理单元、遥测单元和第二处理单元;该图像采集单元用以采集一个或多个活体标本的一个或多个图像;该第一处理单元用以分割该采集图像,从而提供分割图像;该遥测单元用以测量至采集图像的若干部分的距离,从而提供若干距离测量值,并且用以测量活体标本的分割图像所包含的那些距离测量值的子组;该第二处理单元配置成处理分割图像和所述若干距离测量值,所述若干距离测量值指代分割图像内所包含的不同位置。
优选地,图像采集单元和遥测单元布置于相同支撑件内的给定距离处。
第一处理单元和第二处理单元可为独立单元或相同单元。
根据所提出装置,第二处理单元适于并配置成实施根据本发明的第一方面的方法。另外,图像采集单元和遥测单元被校准。
图像采集系统可为RGB相机和/或热相机,该RGB相机具有红色通道中的扩展NIR。遥测单元可为旋转激光雷达、扫描激光雷达、多个激光雷达、飞行时间(TOF)传感器、TOF相机,或是具有或不具有基于单点或多点检测的移动部件的任何其它遥测装置。
附图说明
从以下参考附图对实施例的详细描述中,将更充分地理解前述和其他优点和特征,这些附图必须以说明性和非限制性的方式考虑,其中:
图1和图2为两个流程图,示出了用于从一距离处表征活体标本的方法的两个实施例;
图3示出了三个不同选项,这些选项可单独地或组合地用于获得形状分析图;
图3A示出了三次采集,其中尾部检测为深度剖面的最小距离,并且这与形状分析一致,即在中心列中具有与深度剖面最接近的质心;
图3B示出了相同三次采集,其中尾部检测为深度剖面的最小距离,并且这与形状分析一致,即在中心列中具有与深度剖面最平行的质心;
图3C示出了相同三次采集,其中尾部检测为深度剖面的最小距离,并且这与形状分析一致,即将尾部定位在中心列的右侧;
图4A示出了根据一个实施例来执行质心和估计取向的计算;
图4B示出了两个抛物线在校正取向之后如何拟合至轮廓的左极端(虚线)和右极端(实线);
图5示出了基于欧几里得度量作为轮廓线的二值图像的距离变换的表示轮廓;刻度虚线示出了分割图像的边界;对于较远离任何边界的那些点,图像得分更高;
图6A示出了图像骨架化;
图6B三角形示出端点并且圆圈示出分支点;
图6C直线标记了作为解剖基准点的第一估计的前脚、横跨部、臀部和后脚之间的联系;
图6D额外线标记了与其它基准点的联系,诸如头部和尾部(白色圆圈)以及中间身体宽度;
图7示出了形状分析图的匹配和激光雷达扫描的定位是如何允许了解多个图像点的特定深度;
图8顶行示出了根据一个实施例通过旋转激光雷达采集的分割图像和遥测测量值的采集序列;中间行示出了图像中的距离测量值和水平像素位置;底行示出了距离和像素位置至真实空间的变换;
图9示出了所测量深度剖面分析图的角度校正,以测量特定身体特性并估计测量的质量;
图10A示出了校正后的深度剖面的内插和由形状和深度剖面分析所计算的身体部分的注释,其中可以比较两种分析的解剖位置的计算;以及
图10B示出了具有分割图像的二值轮廓的原始图像,内插有深度剖面分析图和如通过图10A的剖面分析所计算的身体部分;
图11为猪的身体图的实例,其中解剖基准点位于三维轴线内;质心为(0,0,0)并且身体的不同部分(诸如尾部或头部)映射于真实空间中,从而示出例如以厘米为单位的坐标;
图12A和图12B为以两个不同矩所采集的相同猪的两个图像,在(A)的右侧可看到栅栏,而在(B),墙壁是唯一背景;
图12C和图12D为分割图像的骨架化,身体部分通过形状分析和图6的方法来获得;
图13A示出了图12所示的两次采集的重叠二值轮廓和基准点以及身体部分;图13B示出了通过平移和旋转的规范化空间;并且图13C示出了基于基准点的空间变换。
图14A和图14B分别示出了从图8左列和中间列的剖面分析的轮廓、深度剖面、基准点;
图14C示出了在图像(B)上具有校正坐标的基准点和深度剖面的重叠;以及
图14D示出了两个深度剖面分析图的基准点、轮廓和累积的重叠;
图15示出了身体部分或基准点的解剖关系;
图16示出了角度φ和θ的明确关系的投影顶视图和正视图;在该图中,还可观察到,图像能够从任何角度和距离进行采集;
图17A为公牛的图像分割,其中犄角的尖端、犄角和嘴部之间的中心点得以检测,并且基准线与该对点关联地描绘以构建形状图,在其中圆点示出了(图17D)的深度剖面得以测量的图像位置;
图17B示出了基于分支点(圆圈)和端点(三角形)的骨架化和检测的形状分析;
图17C为图17B的放大图以示出分支点和端点的特定位置;以及
图17D为犄角的基部和头部的顶部处的深度剖面;
图18示出了头部的形状分析图以定位犄角并计算总长度;
图19A示出了原始图像;
图19B示出了从精细边缘的所勾画的内轮廓;以及
图19C示出了作为形状分析图的基础的骨架化、分支点和端点;
图20A示出了作为轮廓的分割图像、作为白色圆圈的质心、和作为垂直于如虚线所示取向的轴线的深度剖面分析图的图像位置;以及
图20B为真实空间中的所获得深度剖面分析图;
图21A示出了作为轮廓的分割图像、作为白色圆圈的质心、和作为穿过头部和尾部的迹线(如形状分析图的虚线或特定点所示)的深度剖面的图像位置;以及
图21B为真实空间中的所获得深度剖面分析图;
图22为视频图像,其中热图像和小圆点的重叠分割区域示出了激光雷达测量值;
图23示出了TOF图像的距离(左)和霍夫变换的计算(右)以用于尾部和肩部检测,通过扫描身体和方法2中额外解剖点的计算脊柱描记作为中间点;
图24示出了身体图的处理的一个实例,该处理允许从其它系统(即,热相机)提取额外特征。
具体实施方式
本发明提供了一种用于执行活体标本的自动测量以表征这些活体标本的方法和装置。
图1根据一个实施例以图形示出了所提出方法的流程图。根据该实施例,该方法经由图像采集单元(诸如热相机或RGB相机)采集活体标本(在该特定情况下为猪,非限制性地,可表征任何活体标本)的一个图像,并且还分割该采集图像,从而提供分割图像(步骤a)。同时或之后,该方法经由遥测单元测量至所述采集图像的若干部分的距离,从而提供若干距离测量值,并且选择猪的分割图像中所包含的那些距离测量值的子组(步骤b)。然后,分割图像和距离测量值进行处理(步骤c)。在该特定实施例中,处理步骤包括经由计算形状分析图的算法表征猪的形状(步骤c1);经由计算深度剖面分析图的算法表征猪的深度(步骤c2)。最后,该方法涉及形状分析图和深度剖面分析图的比较(步骤c3)。比较的结果/分数可用于决定猪的参数是否可以足够质量进行计算,和/或身体图(参见图2)(优选为3D表示)是否可以足够质量进行计算,或该方法是否必须进行重复、校正或停止。
采集图像优选地为任何类型的二维图像(例如,灰度、色彩、热量,或色彩和热量两者)。可使用将采集图像转换成分割图像的任何分割方法。分割图像由处理一个图像(例如,灰度、色彩、热量,或其组合)并且将该图像的像素分成两类得到:(1)包含于猪中的像素和(2)未包含于猪中的像素。
分割图像可以不同方式进行编码:(1)二值图,其中猪内所包含的像素设定成最大值并且猪内未包含的像素设定成最小值;(2)二值轮廓,其中猪的边缘内所包含的像素设定成最大值并且猪内未包含的像素设定成最小值;(3)矢量,其中边界的位置以矢量来设定。
遥测单元配置成测量至少两个点的距离,这两个点包含于分割图像内。距离测量值可通过不同方法来获得。例如,遥测单元可通过旋转激光雷达来实施,该旋转激光雷达具有10Hz的自转速度(100ms完成对角度和距离的完全读取)和小于一度的分辨率。图像采集单元激光雷达的预先校准或热相机至可视或近红外相机(其然后校准至激光雷达)的校准允许构建表格,该表格用于将激光雷达坐标(即,角度和所测量距离)变换成图像坐标(即,二维图像的行和列)。另选地,具有仅检测激光雷达波长的特定光学滤波器的专用相机可用于图像坐标和激光雷达信息的精确定位。另选地,遥测单元可通过具有TOF技术的新型相机来实施,这些新型相机向二维图像提供距离。速度超过10帧每秒,并且在一些情况下其可达到1000帧每秒。图像采集单元和TOF传感器或相机的预先校准允许找出图像采集单元的像素和TOF传感器或相机的像素之间的关系。
遥测单元和图像采集单元的校准可通过一对加热电阻器来执行,该对加热电阻器以至平面的两个任意深度位于该平面上。在这种情况下,采集单元为热相机,该热相机以该采集单元平行于平面的方式进行定位,并且热电阻器位于热图像的竖直轴线的中心周围。相对于激光雷达坐标(对于旋转激光雷达,其为角度和距离),旋转激光雷达以每个热电阻器的距离d1和d2相对于初始值dL(例如,2米)进行测量的方式进行调整。随着采集图像中的位置改变了距离dL,该操作对于不同距离dL进行重复。该过程允许构建点表格,该点表格与像素位置和测量距离相关。然后,回归模型构建成使任何激光雷达坐标(角度和距离)相关于采集图像和分割图像中的特定(x,y)位置。
在另一个实例中,对于图像采集单元为热相机并且遥测单元为TOF相机或传感器的特定情况,校准如之前进行,但考虑更多个点并且不仅依赖于旋转激光雷达的扫描线。
其它校准方法也为可能的。例如,图像采集单元(由一个RGB相机和一个热相机组成,该RGB相机具有红色通道中的NIR扩展)和基于旋转激光雷达的遥测单元可一起被校准。
用以计算形状分析图的形状表征包括对以下项的计算:猪的质心、猪在分割图像内相对于基准点的取向,和/或猪的特定身体部分(通过将猪的解剖基准点定位于分割图像内)。应当注意,仅需要上文所指示方法的一者来计算形状分析图。然而,其组合也是可能的。图3A至图3C示出了三种计算的实施例。
为表征猪的形状,猪由分割图像来限定。猪的形状为分割图像的形状。采集图像和分割图像可表达为一系列位置以构建二值图、二值轮廓或轮廓的多点近似值。因此,任一种格式的分割图像s(x,y)可表达如下:
Figure BDA0002740273020000081
其中x,y分别为数字图像的列和行。
在一个实施例中,为计算质心,猪的形状通过图像矩来表征:沿循这种格式,然后可以根据标准公式计算任何图像矩Mnk
Figure BDA0002740273020000082
像素的数量等于矩M00,质心等于(M10/M00,M01/M00)。
这些矩可直接地扩展至中心矩(其为平移不变量)。然后,平移不变矩还可扩展至尺度不变量,并且此类尺度不变量还可通过公知的现有技术计算扩展至旋转不变量(Hu矩不变量)。该组的矩允许计算可与特定形状相关联的特性特征,诸如从特定视角(或取向角度)所见的猪形状。
这些矩还可琐细地扩展至多个维度(例如,3D)以表征3D形状:
Figure BDA0002740273020000091
Figure BDA0002740273020000092
其中x,y,z分别为数字体积的列、行和深度。
为计算取向,分割图像可通过最小二乘拟合、高斯模型、主分量分析、霍夫变换等拟合成椭圆。拟合椭圆的取向、高斯分布的取向、第一主分量的角度或霍夫线的平均取向是用以快速且可靠地估计目标取向的方法。
根据一个实施例(参见图3C),为计算特定身体部分,可使用霍夫变换。霍夫变换可以多种形式来实施。特别地,圆形霍夫变换允许识别圆形区域的半径范围。这可用于区分猪的头部和尾部。由于尾部为圆的,其可拟合至较大圆圈。例如,通过将分割图像考虑为如图3C所示的二值轮廓,圆形霍夫变换可设定成检测具有高敏感性的圆圈并且可仅要求多个点来拟合圆圈。圆圈的范围可通过下述估计:1)大腿半径RH应为分割图像的竖直尺寸的约1/3;2)半径的搜索范围然后设定为RH+/-50%的RH。然后,选择5个圆圈中具有最大评分的较大圆圈,得到居中于尾部部分的圆圈。
霍夫分析可通过半圆形霍夫变换进行扩展,并且可获得半圆圈的拟合,这对于获得尾部和头部差异将是更可靠的。其还扩展至椭圆形状以拟合猪的身体或头部的中心部分或其它部分。广义霍夫变换是以不同尺度和角度拟合多个特定形状来匹配形状的另一种方法。可以以等同方式来使用类似的形状匹配方法。
以类似方式,猪的尾部可通过将二阶多项式拟合至由质心和取向角度所限定的轴线进行检测。图4A示出了质心和取向轴线,图4B校正了取向并且将抛物线(二阶多项式)拟合于两个界限周围。底部通过具有最接近于质心的顶点的抛物线来检测。
形状分析图可通过其它若干计算策略来完善。例如,利用傅里叶分析;在这种情况下,猪的轮廓可在极坐标中进行编码然后进行傅里叶变换。这提供了若干傅里叶系数,这些傅里叶系数的模量为旋转不变量并且其论证包含旋转信息。
图5示出了可使用的另一策略。在这种情况下,分割图像根据分割图像内的任何点至最近边界点的距离进行评分。距离度量可以是多种多样的:欧几里得度量、城市街区度量或任何其它距离度量。
另一策略是通过计算分割图像的骨架化函数而计算猪的轮廓。图像骨架为与其边界等距的该形状的薄型版本。骨架通常强调形状的几何性质和拓扑性质,诸如其连接性、拓扑结构、长度、方向和宽度。连同其点至形状边界的距离一起,骨架还可用作形状的表示(它们包含用以重建该形状所需的所有信息)。然后,分支点和端点可用于估计不同身体部分的解剖位置。
应当注意,用以计算形状分析图的这些互补策略可在它们之间组合地使用。
返回参考图3和图4,这些图示出了形状分析如何允许通过适当“尾部检测”来识别头部至尾部的取向(其为用以使图像点关联于身体部分的基本步骤)的实例。形状的额外表征允许使其它图像点关联于其它身体部分。例如,具有高边界分数(根据距离变换所计算)或邻近质心轴线的分支点可关联于如图5所示的横跨部(肩部)和臀部。尾部检测还在横跨部和臀部之间进行区分。足部通过位于分割图像的底部部分处的端点来确定,并且几乎垂直于由横跨部和臀部(或质心轴线)以及如图6C所示的底部部分处的端点所绘制的线。
关于用以计算深度剖面分析图的猪的深度表征,该过程包括在分割图像中所包含的距离测量值内,通过将活体标本的解剖基准点定位于距离测量值内来计算活体标本的特定身体部分。
该步骤可分成两个主要部分:(1)至图像点的遥测距离的定位通过如上文所描述的预先校准来实现并且允许计算具有深度信息的图像坐标;和(2)图像点和至身体部分的距离的关系。
具有分割图像内所包含的深度信息的图像坐标提供了如图7所示的深度剖面;较大圆点为落入分割图像之内的圆点,而较小圆点处于分割图像之外。形状分析图的结果允许使图像点(和深度信息)关联于身体部分或基准点。例如,从特定基准点(诸如质心、横跨部、臀部、尾部),或从解剖基准所获得的其它基准点)获得的深度信息和其它形状分析(诸如距离变换或骨架化)。
图像采集单元和遥测单元的特定对准允许获得相关3D信息,该3D信息与身体尺寸和所测量身体部分或基准点的较高可靠性有关。例如,旋转激光雷达与质心轴线或沿着图像中的横跨部和臀部所限定的身体线的对准允许扫描重要身体特征以获得特定深度信息。图7示出,当对准时,激光雷达扫描质心轴线,深度剖面得到多峰值曲线。
图8示出了由不同激光雷达对准由横跨部点和臀部点所限定的轴线所采集的图像序列的实际实例。在顶行处,可看出虚线形式的三个二值轮廓(或廓影);轮廓之内的白色十字示出了通过旋转激光雷达所采集的距离测量值的图像位置;圆圈示出最接近于上文所解释的骨架、分支点和额外计算所限定的横跨部位置和臀部位置的激光雷达测量值。中间行以水平像素示出了所测量距离和图像位置。底行示出了从图像位置至真实空间的转换测量值。
通过分析真实空间中剖面曲线,可以确认是否已实际测量并适当地定位基准点(诸如身体部分或解剖基准),或是否可找出更佳估计。
图8的左列示出了其中旋转激光雷达靠近由横跨部和臀部所限定的轴线对猪进行扫描的情形,因为其实际上非常接近地经过了图6标记为臀部的点。然而,旋转激光雷达的角度偏离由横跨部和臀部所限定的轴线。在该上下文中,深度剖面分析图似乎缺少重要身体特征。图8的右列示出了其中旋转激光雷达远离由横跨部和臀部所限定的轴线对猪进行扫描的情形,其中不可能获得猪的形状的任何深度信息。图8的中间列示出了其中旋转激光雷达沿由横跨部和臀部所限定的轴线对猪进行扫描的情形。如在真实空间的剖面深度可看出,测量值包含关于后腿的信息。猪的后腿为猪的重要部分,并且猪的后腿的厚度在许多情况中为相关特征。例如,伊比利猪的后腿的尺寸对估计整猪的市场价值是重要的。图9示出了靠近由真实空间的横跨部和臀部所限定的轴线的测量剖面,该轴线通过估计剖面的线性分量进行校正。这允许通过估计正确图像中从横跨部和臀部之间的最远位置至臀部和尾部之间的最近位置的深度差值来测量大腿的厚度。在臀部和尾部之间的最近位置发现的该负峰值可认为是新基准点,称为“最大臀部”,该基准点是后腿的厚度最大的位置。该“最大臀部”的位置可设定为在测量中必须满足约束条件,以验证整个取向估计、执行任何身体部分测量或存储深度信息。
当扫描通过关键点(例如,穿过分割图像上通过形状分析所估计的横跨部点和臀部点)时,深度剖面分析图还可用于估计基准点或身体部分的精确位置。期望深度剖面或样条内插的“先验”模型可用于该目的。图10示出了图9的深度剖面测量值的样条内插,倾斜方形示出了通过形状分析所计算的横跨部位置和臀部位置(如图8所示),新身体部分根据内插深度剖面进行重新计算,其中:横跨部位置为靠近头部侧的最小距离,最大臀部为靠近尾部侧的最小距离,并且臀部为最大臀部和朝向尾部侧靠近质心的最大点之间的转折点。
以更普遍方式,可以通过将深度测量值参考到基准点或身体部分来完成图像点和至身体部分的距离的关系,从而允许组合一系列测量值的信息以具有如图11所示的动物或复杂目标的更完整3D图片,即,身体图。这可以通过以下方式实现:(1)从如图6D所示的分割图像的形状分析获得基准点或身体部分;(2)计算规范化基准空间;(3)基于此类基准点计算至规范化基准空间的空间变换;(3)将此类空间变换应用于所有采集;(4)将所有深度剖面累积至规范化基准空间。
图12示出了用以从分割图像的形状分析获得基准点或身体部分的第一步骤。图13A示出了二值轮廓和基准点在分割图像的坐标中的重叠。第一规范化空间可通过二值轮廓至质心的直接平移进行计算。然后通过校正如上文所计算的取向而进行点的旋转。结果呈现于图13B中。关于在类似距离所获得的两次采集,标度无需进行校正,但在一些其它采集中,标度可需要进行校正。可以通过以不同手段(多项式函数、局部加权映射或分段),基于本领域已知的基准点对构建空间变换,来实现坐标的第二细化。三阶多项式映射呈现于图13C中。
取决于两次采集之间的基准点的一致性,重叠的采集可能被拒绝并且可请求另一次采集。在参照图8所示的形状分析得出的特定基准点时,可将可接受的采集限制为读取满足某些预期规则的深度信息剖面,其中在使旋转激光雷达扫描通过所估计横跨部和臀部位置方面,中间列显示了可接受的采集并且深度剖面包含与大腿和肩部的厚度有关的期望峰值。在该线上,可要求旋转激光雷达进行多次扫描以捕获足够的深度和形状信息,以根据关于动物、人或复杂目标的已知信息计算动物或复杂目标的取向和尺寸。
当使用3D相机和分割图像时,上述解释也适用。例如,热相机和TOF相机可用于估计动物的身体取向和尺寸。热相机可用于生成分割图像,该分割图像根据上文所描述的过程进行处理。TOF相机将提供多个点的深度信息,但分割图像的形状分析将提供必需身体情况以接受采集。以类似方式,TOF相机允许执行单次采集的多线扫描,并且这可显著地加快整体采集时间。
图15示出了形状分析图和/或图像坐标与深度信息的关系,以获得身体部分或基准点,例如,如图所示的横跨部、大腿的起始点和尾部。
在一个实施例中,所提出方法还包括估计图像采集单元和活体标本的相对位置的三维信息的一部分以获得图像采集单元和猪之间的至少一个角度(θ或φ)的平均值(参见图16)(例如,通过计算深度剖面分析图的线性近似值的斜率的反正切值)。另外,还可获得猪的形状的弯曲度或平面度。平面度或弯曲度可通过扩展图9所示的方法进行估计,在其中可以拟合线性函数以估计φ角度。并且,该原理可例如通过二阶或三阶的多项式扩展至任何形状。调整后的R平方系数可用于评估二次函数相比于线性模型是否更好地拟合。当二次函数更可能拟合时,意味着动物为弯曲的并且测量需进行重复或适当地校正。伊比利猪相比于其它类型的猪为高肌肉动物,它通常会弯曲其身体并且采用保护性的形状。很多情况下必须考虑到这一点,以在表征中克服这种错误源。
还可获得图像采集单元相对于地面的高度或图像采集单元相对于猪的高度的高度。在第一种情况下,额外遥测单元还可提供额外距离测量手段以估计相对高度,图像采集单元和遥测单元在该相对高度处操作。在第二种情况下,由于至少一个距离通过遥测单元进行测量并且分割图像直接地关联于距离测量值,所以可以估计动物高度。总动物高度可计算如下:(1)如上文所描述进行取向校准之后的分割轮廓的垂直范围;(2)通过距离和竖直视场或校准的关系来转换像素数量的计算。如果旋转激光雷达配置成竖直地扫描或遥测单元提供距离的2D图像,那么利用基准点或身体部分,将可以提取3D坐标并且将高度计算为坐标点之间的距离。以类似方式,可以从基准点或身体部分估计高度,例如上文所描述的“最大臀部”至也如上文所描述的后脚。然后,另一种校准方法,可根据视场和距离的关系、将像素的数量从竖直旋转激光雷达或距离的2D图像转换为坐标距离。
甚至,可获得图像采集单元的光学轴线相对于地面的角度。
在一个实施例中,所提出方法还允许计算不同身体部分或基准点之间的关系以获得猪的身体分析。所有基准点或身体部分可用于构建如图6所示的猪的简化图。校准允许在真实空间中计算3D坐标的任何图像点,这允许直接估计复杂目标、动物和人的取向和尺寸。头部至尾部作为总动物长度,横跨部至尾部作为身体长度,横跨部至臀部作为短身体长度,臀部至后脚作为大腿长度,横跨部至前脚作为前腿长度,中间身体顶部和底部作为动物宽度。图9还示出了如何从深度剖面分析图估计大腿厚度,该大腿厚度为伊比利猪的重要特征。另外,对于横跨部、臀部或最大臀部的定位所做出的校正可提供更可靠或更感兴趣的尺寸和取向测量。另外,还可进行区域测量和体积测量。大腿的区域测量可通过仅将分割图像的区域保持在臀部点之外来实现。通过添加深度剖面分析图信息,还可产生大腿的体积估计。类似地,可通过将横跨部之外的区域作为基准点来获得身体区域和体积。
所有上下文数据(诸如采集角度Φ、最小距离、最大距离和平均距离、不同尺寸,和不同尺寸(诸如长度、区域或体积)之间的关系)可用于生成猪的一系列特征。
在本文档中,猪的形状充当复杂形状的一般实例。因为其它动物(诸如牛、雏鸡、肉鸡、公牛、母牛、羊)作为牲畜动物,所以将特别地适合该方法。也可在这些参考条件下对人类进行建模,并且复杂物体可能需要特定的适应,因为复杂物体是一个广义术语。然而,可根据这种方法直接地适应遵循具有明显基准点的图案的对象,这些基准点不只是简单的方形、三角形或圆形,该方法结合了分割图像的形状分析和深度信息,该深度信息通过已适当校准的遥测装置(激光雷达、旋转激光雷达、扫描激光雷达、TOF相机,或提供1D或2D距离序列的任何装置)获得。
图17至图20示出了其它牲畜(诸如斗牛或肉鸡雏鸡)的不同实例。在这些图中,示出了用以构建身体图或身体图的一部分的类似方法。例如,图17A示出了解剖基准点如何用于构建斗牛的头部的基准图以测量犄角的距离,该犄角为用于确立斗牛的价值的重要特征。如图17B所示,深度剖面分析图还可用于计算头部取向并且准确地计算此类距离,或甚至计算犄角的长度。相同骨架化过程的应用在图19中示出用于肉鸡雏鸡。在一些情况下,可能重要的是准确地测量肉鸡的宽度,如图20所示。其它努力可要求在图21所示的其它轴线上进行准确测量,或甚至组合两个轴线或2D表面的信息。
关于图17A至图17D,公牛的犄角可检测为头部内或嘴部之上的左顶位置和右顶位置。因为嘴部具有明确不同的色彩外观或温度,所以嘴部质心可通过对所采集图像进行色彩分析或热分析来检测。根据一个实施例,通过参考形状分析图中的分支点和端点,腿部也可通过类似于猪的形状分析进行测量(参见图17B)和检测。通过拟合圆形霍夫变换或二次函数,尾部也可以如上文所描述的与猪相同的方式进行检测。
头部取向可以如上文所描述的类似于猪的方式通过深度剖面分析进行估计。
通过考虑到头部取向并校正由(x,y)距离和深度引起的图像失真,,可成功地计算犄角的尖端之间的距离测量值。如图18所示,可计算出犄角的总长度的额外信息。头部的对称性质可用于定位特定基准点或轴线。
肉鸡或雏鸡也可适应于所提出方法。例如,图19示出了肉鸡(19A),其分割图像和轮廓(19B)和骨架化(19C)(其中分支点和端点的检测用以构建形状分析图)。
质心和垂直于取向的轴线可用作基准以在短轴线上获得深度剖面。类似地,从分支点所获得的头部和尾部可用于识别长轴线并在其它方向上获得深度剖面信息。顶部的图像或TOF相机的使用允许两个深度剖面从相同采集的计算。
此外,分割图像之外的点可用于计算肉鸡的高度。
深度剖面分析图的实例呈现为线性测量值,但捕获2D图像的TOF相机或旋转激光雷达的若干扫描线(这些扫描线在解剖学上映射至身体图)的累积允许执行其它计算,诸如拟合椭圆。当将这种方法相比于以卷尺和量具实际取得的动物高度时,线性测量为更类似的测量。然而,这种方法不限于线性测量,并且TOF信息可用于拟合表面。
另外,解剖点还可通过图像信息来协助。例如,肉鸡的头部相比于身体为较温暖的,所以该特征可用于直接地定位头部。以类似方式,头部相比于其它身体部分通常为较高的,这可以通过遥测或图像位置加以利用。
接下来,详细介绍所提出方法的不同实例:
—实例1:侧视图或倾斜视图的作为形状分析图的质心
热相机、视频相机和激光雷达已校准。因此,该方法包括步骤a):以热相机采集图像并且通过温度阈值分割一头猪。然后,步骤b),该方法包括以旋转激光雷达测量至极坐标(ρ,Φ)中的若干点的距离和使极坐标中的激光雷达测量值相关于图像内的特定像素位置(x,y)。在步骤c1),猪的质心计算为作为中心圆点的分割图像(x0,y0)的质量中心,如图22A所示。在步骤c2),该方法找出如图7、图8、图9和图10A所示的局部最小距离、局部最大距离和转折点,并且计算尾部、最大臀部(局部最小距离)、大腿端部(转折点)和横跨部(局部最小距离)的深度分析图。最后,在步骤c3),该方法检查任何深度点是否在质心附近行进或分割图像内所包含的深度点是否处于距离y0,例如,|y-y0|<30。如果这是事实,那么通过深度分析图所检测的解剖点可接受为正确的。
—实例2:侧视图或倾斜视图的由形状分析图所检测的多个解剖点
热相机、视频相机和激光雷达已校准。该方法包括步骤a):以热相机采集图像并且通过温度阈值分割一头猪。然后,在步骤b),该方法包括以旋转激光雷达测量至极坐标(ρ,Φ)中的若干点的距离和使极坐标中的激光雷达测量值相关于图像内的特定像素位置(x,y)。在步骤c1),该方法通过计算霍夫变换而计算分割图像的质心以定位尾部。如果所检测圆圈的中心处于具有质心的距离范围内,那么该方法还计算如图6所示的分割图像的骨架,从而检测分支点和端点。尾部附近的分支点为最大臀部解剖点的近似值。质心或最大臀部的类似高度处的其它侧的分支点为肩部点,质心之下和尾部附近的最低端点为后腿,质心之下和与尾部相对的足部为前腿足部(还最靠近于肩部位置以与头部进行区分,此时猪正嗅探地面)。这得到如图6C所示的简图。更复杂图可如图11、图12、图13和图14所示进行计算。在步骤c2),该方法找出如图7、图8、图9和图10A所示的局部最小距离、局部最大距离和转折点,并且计算尾部、最大臀部(局部最小距离)、大腿端部(转折点)和横跨部(局部最小距离)的深度分析图。最后,在步骤c3),该方法检查分割图像内的激光雷达测量值根据形状分析是否在肩部和最大臀部附近交叉,并且还检查两个图中常见的所有解剖点是否处于期望位置附近或偏移。如果这些条件是事实,那么通过深度分析图所检测的解剖点可接受为正确的。
—实例3:倾斜视图的热相机和TOF相机
热相机和TOF相机已校准。该方法包括步骤a):以热相机采集图像并且通过温度阈值分割一头猪。然后,在步骤b),该方法包括以TOF相机测量至若干点的距离,计算真实空间中的(rx,ry,rz)位置,和使TOF测量值相关于图像内的特定像素位置(x,y)。然后,在步骤c1),该方法经由霍夫变换而计算分割图像的质心和取向以定位尾部。如果所检测圆圈的中心处于具有质心的距离范围内,那么该方法还计算如图6所示的分割图像的骨架,从而检测分支点和端点。尾部附近的分支点为最大臀部解剖点的近似值。质心或最大臀部的类似高度处其它侧的分支点为肩部点,质心之下和尾部附近的最低端点为后腿,质心之下和与尾部相对的足部为前腿足部(还最靠近于肩部位置以与头部进行区分,此时猪正嗅探地面)。这得到如图6C所示的简图。更复杂图可如图11、图12、图13和图14所示进行计算。然后,在步骤c2,该方法通过在质心附近和以类似于形状分析的取向挑选TOF图像线而提取深度剖面,或在分割图像的表面内执行深度剖面的2D分析,换句话讲,分析分割图像内所包含的(rx,ry,rz)点。接下来,该方法找出如图7、图8、图9和图10A所示的局部最小距离、局部最大距离和转折点,并且计算尾部、最大臀部(局部最小距离)、大腿端部(转折点)和横跨部(局部最小距离)的深度分析图。还可以将动物拟合至表面模板并且定位不同解剖点。最后,在步骤c3),该方法检查两个图中常见的所有解剖点是否处于期望位置附近或偏移。如果这些条件是事实,那么通过深度分析图所检测的解剖点可接受为正确的。
—实例4:鸟瞰图和倾斜视图的热相机和TOF相机
校准热相机和TOF相机。该方法包括步骤a):以热相机采集图像并且通过温度阈值分割一头猪。然后,在步骤b),该方法包括以TOF相机测量至若干点的距离,计算真实空间中的(rx,ry,rz)位置,和使TOF测量值相关于图像内的特定像素位置(x,y)。在步骤c1),该方法利用霍夫变换而执行形状分析以检测肩部和尾部。尾部以多种形式与肩部进行区分,例如,超出肩部的区域(头部)远远大于尾部(仅尾部)。另选地,如图23中方法2所示,因为尾部的梯度远远高于头部,所以轮廓分析允许直接检测尾部,。从肩部中心至尾部中心对图像进行扫描允许确定脊柱。形状分析图包括例如肩部位置、头部和脊柱点,如图23中方法1所示;并且还可扩展至方法2(如果要求额外解剖点的话)。在步骤c2),该方法计算解剖点的高度。最后,在步骤c3),该方法检查形状分析图的所有解剖点是否处于正确高度或高于阈值。
—实例5:鸟瞰图和倾斜视图的TOF相机
校准TOF相机。该方法包括步骤a):以TOF相机采集图像并且以相比于背景的距离阈值分割。然后,在步骤b),该方法包括以TOF相机测量至若干点的距离,计算真实空间中的(rx,ry,rz)位置,和使TOF测量值相关于图像内的特定像素位置(x,y)。在这种情况下,步骤c1等同于实例4。在步骤c2),深度分析图确保所有分割区域都位于地面的给定高度之上。最后,在步骤c3),如果找出形状分析的所有点,那么这意味着它们处于正确距离,因为其为分割(步骤a)的先决条件。此外,还可以包括其它计算(诸如计算脊柱的rx点、ry点、rz点的曲率),和向此类曲率给出特定公差。
—实例6:用于额外特征的TOF图像和热图像的组合(用于身体图)。如果TOF相机和热相机校准,那么不同身体部分处热特征的额外计算可如图24所示进行计算。
还提供了一种用于远程表征活体标本的装置。该装置主要包括所述图像采集单元,分割装置、所引用遥测单元和处理装置,该处理装置用以处理所描述不同信息/数据以允许对活体标本或复杂目标进行表征。该装置还可包括存储器以存储所处理的不同测量值或信息。
本发明可以硬件、软件、固件或其任何组合来实施。如果以软件来实施,那么功能可存储于计算机可读介质中或作为一个或多个指令或代码在计算机可读介质中进行编码。
计算机可读介质包括计算机存储介质。存储介质可为可由计算机访问的任何可用介质。通过示例和非限制的方式,此类计算机可读介质可包括RAM、ROM、EEPROM、CD-ROM或其它光盘存储装置、磁盘存储装置或其它磁性存储装置,或可用于承载或存储指令或数据结构形式的期望程序代码并且可由计算机访问的任何其它介质。如本文所用,磁盘和光盘包括压缩光盘(CD)、激光盘、光学光盘、数字通用光盘(DVD)、软盘和蓝光光盘,其中磁盘通常以磁性方式复制数据,而光盘则通过光学方式以激光复制数据。上述的组合还应包括于计算机可读介质的范围内。任何处理器和存储介质可驻留于ASIC中。ASIC可驻留于用户终端中。另选地,处理器和存储介质可作为分立部件驻留于用户终端中。
如本文所用,计算机程序产品包括计算机可读介质,该计算机可读介质包括所有形式的计算机可读介质,但此类介质被视为非法定的瞬时传播信号的除外。
本发明的范围限定于下述权利要求书中。

Claims (15)

1.一种用于从一距离处表征活体标本的方法,所述方法包括:
a)经由图像采集单元采集至少一个活体标本的一个图像和通过处理单元进一步分割所采集的图像,从而提供分割图像;
b)通过遥测单元测量至所述采集图像的若干部分的距离,从而提供若干距离测量值,和选择所述活体标本的分割图像中所包含的那些距离测量值的子组,
其中所述图像采集单元和所述遥测单元进行校准;
c)通过处理单元处理所述分割图像和所述若干距离测量值,所述若干距离测量值通过以下方式指代所述分割图像内所包含的不同位置:
c1)通过实施算法而表征所述活体标本的形状,所述算法在所述分割图像内至少计算以下项的一者或多者:
·所述活体标本的质心,和/或
·所述活体标本在所述分割图像内相对于基准点的取向,和/或
·通过将所述活体标本的解剖基准点定位于所述分割图像内确定的所述活体标本的特定身体部分,
所述形状表征的结果提供了形状分析图;和
c2)通过实施算法而表征所述活体标本的深度,所述算法在所述分割图像中所包含的距离测量值内至少计算:
·通过将所述活体标本的解剖基准点定位于所述距离测量值内确定的所述活体标本的特定身体部分,从而提供至少一个深度剖面分析图;和
c3)比较所述形状分析图和所述深度剖面分析图,其中:
—如果所述比较的结果包括于给定范围之内,那么所述方法还包括确定所述活体标本的参数、包括姿势参数,和/或所述解剖基准点的定位或校正,和/或身体尺寸参数,和/或表示所述活体标本的身体图;
—如果所述比较的所述结果包括于所述给定范围之外,那么所述方法还包括使步骤a)至步骤c)重复,从而获得新深度剖面分析图和新形状分析图;和
其中步骤a)和步骤b)在任何时间执行。
2.根据权利要求1所述的方法,其中步骤a)包括在不同时间段采集所述活体标本的若干图像,捕获所述活体标本的不同姿势,和对于每个采集图像获得一系列距离测量值。
3.根据前述权利要求中任一项所述的方法,还包括估计所述图像采集单元和所述活体标本的相对位置的三维信息的一部分,以获得下述参数的至少两者:所述图像采集单元和所述活体标本之间的至少一个角度的平均值,所述活体标本的所述形状的弯曲度或平面度,所述图像采集单元与地面高度的相对高度或所述图像采集单元与所述活体标本高度的相对高度,和/或所述图像采集单元的光学轴线相对于所述地面的角度。
4.根据权利要求1所述的方法,其中所述活体标本的所述取向经由二乘拟合函数、高斯模型、主分量分析、最小面积矩形或霍夫变换,通过将所述分割图像拟合成椭圆进行计算。
5.根据权利要求1所述的方法,其中所述形状分析图的所述身体部分通过以下项进行计算:
—圆形霍夫变换,所述圆形霍夫变换计算所述分割图像内的大腿或股部的半径;或
—二阶多项式函数,所述二阶多项式函数通过拟合所述质心和取向轴线周围的抛物线而检测所述分割图像内的所述活体标本的尾部。
6.根据权利要求1或5所述的方法,其中所述形状分析图的所述身体部分还通过以下方式进行计算:
—计算所述分割图像内的所述活体标本的轮廓,在极坐标上编码所述计算轮廓并且进一步将傅里叶变换函数应用至所述极坐标,提供若干傅里叶系数,所述傅里叶系数的模量为旋转不变量并且其论证包含旋转信息;和/或
—计算所述活体标本的轮廓并且进一步基于包括欧几里得、测地或城市街区的距离度量来计算所述计算轮廓内的距离度量;和/或
—计算所述分割图像的距离变换;和/或
—通过从所述分割图像计算骨架化函数而计算所述活体标本的轮廓,提供所述活体标本的骨架的图像,和进一步可选地计算所述骨架内的分支点和端点以估计不同身体部分的解剖位置;和/或
—计算所述分割图像的图像矩。
7.根据权利要求2所述的方法,还包括:
—组合每次采集所获得的信息,获得一系列成对深度剖面分析图和形状分析图;
—对每对图指定分数;和
—选择具有最高分数的所述对。
8.根据权利要求2所述的方法,还包括:
—组合每次采集所获得的信息,获得一系列成对深度剖面分析图和形状分析图;
—对每对图指定分数;和
—匹配所有采集内的解剖基准点并且累积所述深度剖面分析图和解剖基准点的不同片段以计算所述活体标本的三维重建,或
—计算每次采集的身体图并且累积每个身体图的所有信息,并且计分至扩展身体图中。
9.根据权利要求1所述的方法,其中所述活体标本的身体图得以表示,所述身体图进一步被用于基于所述图像采集单元的特征,包括色彩和/或温度,而计算所述身体图的特性,和/或计算由预先校准的额外装置所测量的额外特性。
10.根据权利要求1所述的方法,其中步骤a)的所述采集图像包括两个或更多个活体标本,并且所述方法包括为图像中包括的每个活体标本计算并比较所述形状分析图和所述深度剖面分析图,从而在同一采集中评估所述两个或更多个活体标本。
11.根据权利要求1、8或9所述的方法,其中所述身体图为所述活体标本的3D表示。
12.根据前述权利要求中任一项所述的方法,其中所述活体标本为牲畜动物,包括猪、公牛、母牛、羊、肉鸡、鸭或雏鸡。
13.一种用于从一距离处表征活体标本的装置,所述装置包括:
—图像采集单元,所述图像采集单元配置成采集至少一个活体标本的至少一个图像;
—第一处理单元,所述第一处理单元配置成分割所述采集图像,从而提供分割图像;
—遥测单元,所述遥测单元配置成测量至所述采集图像的若干部分的距离,从而提供若干距离测量值,并且配置成测量所述活体标本的所述分割图像所包含的那些距离测量值的子组,
其中所述图像采集单元和所述遥测单元被校准;
—第二处理单元,所述第二处理单元配置成通过执行根据权利要求1的步骤c)来处理分割图像和若干距离测量值,所述若干距离测量值指代所述分割图像内所包含的不同位置。
14.根据权利要求13所述的装置,其中所述图像采集系统包括RGB相机和/或热相机,所述RGB相机具有红色通道中的扩展NIR。
15.根据权利要求13或14所述的装置,其中所述遥测单元包括旋转激光雷达、扫描激光雷达、多个激光雷达、飞行时间TOF传感器或TOF相机。
CN201980028070.9A 2018-02-26 2019-02-14 用于从一距离处表征活体标本的方法和装置 Pending CN112262408A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP18382110.7 2018-02-26
EP18382110.7A EP3531373B1 (en) 2018-02-26 2018-02-26 A method and device for the characterization of living specimens from a distance
PCT/EP2019/053654 WO2019162180A1 (en) 2018-02-26 2019-02-14 A method and device for the characterization of living specimens from a distance

Publications (1)

Publication Number Publication Date
CN112262408A true CN112262408A (zh) 2021-01-22

Family

ID=61283149

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201980028070.9A Pending CN112262408A (zh) 2018-02-26 2019-02-14 用于从一距离处表征活体标本的方法和装置

Country Status (16)

Country Link
US (1) US11587350B2 (zh)
EP (1) EP3531373B1 (zh)
JP (1) JP7336691B2 (zh)
KR (1) KR20200127016A (zh)
CN (1) CN112262408A (zh)
AU (1) AU2019223705A1 (zh)
BR (1) BR112020017434A2 (zh)
CA (1) CA3092209A1 (zh)
DK (1) DK3531373T3 (zh)
ES (1) ES2921128T3 (zh)
HU (1) HUE059001T2 (zh)
MX (1) MX2020008861A (zh)
PH (1) PH12020551508A1 (zh)
PL (1) PL3531373T3 (zh)
PT (1) PT3531373T (zh)
WO (1) WO2019162180A1 (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110910393B (zh) * 2018-09-18 2023-03-24 北京市商汤科技开发有限公司 数据处理方法及装置、电子设备及存储介质
WO2020261403A1 (ja) * 2019-06-26 2020-12-30 日本電気株式会社 身長推定装置、身長推定方法及びプログラムが格納された非一時的なコンピュータ可読媒体
EP3756458A1 (en) * 2019-06-26 2020-12-30 Viking Genetics FmbA Weight determination of an animal based on 3d imaging
CN111145240A (zh) * 2019-11-18 2020-05-12 西宁市动物疫病预防控制中心(挂西宁市畜牧兽医站牌子) 一种基于3d相机的活体西门塔尔牛体尺在线测量方法
CN113316803A (zh) * 2019-12-20 2021-08-27 博医来股份公司 使用历史校正的统计分析来校正医学图像的分割
CN113915740B (zh) * 2020-07-08 2023-12-22 海信空调有限公司 一种空调器和控制方法
US11892574B2 (en) * 2020-07-31 2024-02-06 GM Global Technology Operations LLC Dynamic lidar to camera alignment
CN111986168B (zh) * 2020-08-07 2024-03-15 中国农业大学 T形分支的双耳根检测与最佳耳根测温帧检测方法及系统
CN112927282A (zh) * 2021-01-25 2021-06-08 华南农业大学 一种基于机器视觉的畜禽脚参数自动测量方法
CN113632749B (zh) * 2021-09-07 2022-11-22 华中农业大学 一种基于动物行为的猪身体状态监测方法
CN113808156B (zh) * 2021-09-18 2023-04-18 内蒙古大学 一种户外牛体体尺检测方法与装置
CN114243925B (zh) * 2021-12-21 2024-02-09 国网山东省电力公司淄博供电公司 基于智能融合终端的台区配变态势感知方法及系统
CN114838665B (zh) * 2022-07-04 2022-09-02 江西农业大学 一种基于黑皮鸡枞菌的尺寸原位测量方法
CN116363141B (zh) * 2023-06-02 2023-08-18 四川省畜牧科学研究院 一种妊娠母猪智能体型评价装置及系统

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5412420A (en) 1992-10-26 1995-05-02 Pheno Imaging, Inc. Three-dimensional phenotypic measuring system for animals
US5474085A (en) 1994-02-24 1995-12-12 University Of Prince Edward Island Remote thermographic sensing of livestock
IL150915A0 (en) * 2002-07-25 2003-02-12 Vet Tech Ltd Imaging system and method for body condition evaluation
US6974373B2 (en) * 2002-08-02 2005-12-13 Geissler Technologies, Llc Apparatus and methods for the volumetric and dimensional measurement of livestock
WO2008130906A1 (en) * 2007-04-17 2008-10-30 Mikos, Ltd. System and method for using three dimensional infrared imaging to provide psychological profiles of individuals
NL1034292C2 (nl) * 2007-08-27 2009-03-02 Maasland Nv Systeem en werkwijze voor het beheren van een groep dieren.
US9684956B2 (en) * 2008-12-03 2017-06-20 Delaval Holding Ab Arrangement and method for determining a body condition score of an animal
AU2010219406B2 (en) * 2010-05-19 2013-01-24 Plf Agritech Pty Ltd Image analysis for making animal measurements
WO2014083433A2 (en) 2012-12-02 2014-06-05 Agricam Ab Systems and methods for predicting the outcome of a state of a subject
EP2957861A1 (en) 2014-06-17 2015-12-23 Expert Ymaging, SL Device and method for automated parameters calculation of an object
US20170196196A1 (en) * 2016-01-08 2017-07-13 Leo Trottier Animal interaction devices, systems and methods
EP3646287A1 (en) * 2017-06-29 2020-05-06 The GSI Group, Llc Regression-based animal weight estimation

Also Published As

Publication number Publication date
PT3531373T (pt) 2022-08-16
BR112020017434A2 (pt) 2020-12-22
RU2020131518A (ru) 2022-03-29
CA3092209A1 (en) 2019-08-29
PH12020551508A1 (en) 2021-09-06
AU2019223705A1 (en) 2020-10-01
US20210004577A1 (en) 2021-01-07
ES2921128T3 (es) 2022-08-18
JP2021515323A (ja) 2021-06-17
EP3531373A1 (en) 2019-08-28
PL3531373T3 (pl) 2022-08-08
EP3531373B1 (en) 2022-05-04
JP7336691B2 (ja) 2023-09-01
WO2019162180A1 (en) 2019-08-29
MX2020008861A (es) 2020-10-28
KR20200127016A (ko) 2020-11-09
HUE059001T2 (hu) 2022-09-28
US11587350B2 (en) 2023-02-21
DK3531373T3 (da) 2022-07-11

Similar Documents

Publication Publication Date Title
CN112262408A (zh) 用于从一距离处表征活体标本的方法和装置
Wang et al. A portable and automatic Xtion-based measurement system for pig body size
CA2744146C (en) Arrangement and method for determining a body condition score of an animal
Shuai et al. Research on 3D surface reconstruction and body size measurement of pigs based on multi-view RGB-D cameras
US8340402B2 (en) Device and method for detecting a plant
JP6777948B2 (ja) 健康状態推定装置
Wang et al. Automated calculation of heart girth measurement in pigs using body surface point clouds
CN109636779A (zh) 识别家禽体积体尺的方法、装置及存储介质
JP2019211364A (ja) 動物体の体重推定装置及び体重推定方法
Azouz et al. Development of a teat sensing system for robotic milking by combining thermal imaging and stereovision technique
CN112288793A (zh) 畜类个体背膘检测方法、装置、电子设备及存储介质
Hansen et al. Non-intrusive automated measurement of dairy cow body condition using 3D video
US9171228B2 (en) Method and system for estimating a similarity between two binary images
CN113920138A (zh) 一种基于rgb-d相机的奶牛体尺检测装置及其检测方法
CN116071417B (zh) 一种基于Azure Kinect的羊只体尺体重获取系统及方法
CN112528772A (zh) 基于热红外的奶牛跛行运动特征检测方法及装置
RU2785149C2 (ru) Способ и устройство для получения характеристик живых особей на расстоянии
US11341661B2 (en) Method and apparatus for registering live medical image with anatomical model
JP7136344B2 (ja) カメラの較正方法、カメラ及びプログラム
Du et al. Automatic heart girth measurement for cattle based on deep learning
Song et al. Methods of three-dimensional reconstruction and body size measurement of cattle based on multi-view binocular camera
GR1010400B (el) Μεθοδος και συστημα μη επεμβατικης μετρησης μεγεθους ψαριων ιχθυοκαλλιεργειας
Ben Azouz Development of teat sensing system for automated milking

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