EP3043710A1 - Procédé de quantification du comportement d'un animal - Google Patents
Procédé de quantification du comportement d'un animalInfo
- Publication number
- EP3043710A1 EP3043710A1 EP14761892.0A EP14761892A EP3043710A1 EP 3043710 A1 EP3043710 A1 EP 3043710A1 EP 14761892 A EP14761892 A EP 14761892A EP 3043710 A1 EP3043710 A1 EP 3043710A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- mesh
- animal
- labeled
- animals
- batch
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 241001465754 Metazoa Species 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000036544 posture Effects 0.000 claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 230000000007 visual effect Effects 0.000 claims description 29
- 238000005457 optimization Methods 0.000 claims description 9
- 241000894007 species Species 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 abstract description 4
- 241000699666 Mus <mouse, genus> Species 0.000 description 30
- 230000006399 behavior Effects 0.000 description 23
- 238000004422 calculation algorithm Methods 0.000 description 17
- 238000004458 analytical method Methods 0.000 description 11
- 241000699670 Mus sp. Species 0.000 description 8
- 238000010171 animal model Methods 0.000 description 7
- 210000000988 bone and bone Anatomy 0.000 description 6
- 230000033001 locomotion Effects 0.000 description 5
- 230000007935 neutral effect Effects 0.000 description 5
- 241000700159 Rattus Species 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 241000283984 Rodentia Species 0.000 description 2
- 210000000038 chest Anatomy 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 238000012896 Statistical algorithm Methods 0.000 description 1
- 241000251539 Vertebrata <Metazoa> Species 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000004197 pelvis Anatomy 0.000 description 1
- 230000001144 postural effect Effects 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
- A61B5/1104—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb induced by stimuli or drugs
- A61B5/1105—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb induced by stimuli or drugs of laboratory animals, e.g. activity
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
- A01K1/00—Housing animals; Equipment therefor
- A01K1/02—Pigsties; Dog-kennels; Rabbit-hutches or the like
- A01K1/03—Housing for domestic or laboratory animals
- A01K1/031—Cages for laboratory animals; Cages for measuring metabolism of animals
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
- A01K29/00—Other apparatus for animal husbandry
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
- A61B5/0015—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by features of the telemetry system
- A61B5/002—Monitoring the patient using a local or closed circuit, e.g. in a room or building
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/40—Animals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/42—Evaluating a particular growth phase or type of persons or animals for laboratory research
Definitions
- the subject of the invention is that of methods for quantifying the behavior of animals, and in particular laboratory animals.
- Biological research uses laboratory animals, such as rodents and, in particular, rats and mice. These animals are used, for example, for testing active molecules, or for identifying the physiological role of genes or target proteins.
- rodents such as rodents and, in particular, rats and mice.
- mice These animals are used, for example, for testing active molecules, or for identifying the physiological role of genes or target proteins.
- proteomics studies we measure the phenotype, that is, the behavioral and functional consequences of genetic mutations (phenotyping).
- This procedure takes place over a short period of time that does not exceed about ten minutes per day. Thus, often the behaviors that would be relevant are not observed because occurring outside the observation period.
- Such systems use a two-dimensional image analysis algorithm (2D), able to extract the contour of the animal at each sampling time.
- 2D two-dimensional image analysis algorithm
- 3D posture three-dimensional in what follows, and the evolution of this 3D posture during time.
- the present invention aims to meet the aforementioned need for quantifying the behavior of a laboratory animal.
- the subject of the invention is a method for quantifying the behavior of an animal of a given species, characterized in that it comprises the steps:
- acquisition method comprising acquiring a measurement of an external surface of the body of the observed animal from a plurality of frames obtained simultaneously by means of a plurality of imaging devices;
- the labeled envelope template for optimizing the parameters of a labeled envelope template, the labeled envelope template associating a labeled mesh and a control volume including labeled mesh deformation cages, the labeled mesh including control points associated with tags features of a skeletal model of animals of said species, the optimization step consisting of an adjustment of the labeled mesh to the measurement, so as to generate an optimal labeled mesh;
- the method comprises one or more of the following characteristics, taken separately or in any technically possible combination:
- the steps are iterated for each animal of a batch of animals to obtain an overall quantity representative of the behavior of the animals of said batch of animals;
- the method comprises a step of constructing a labeled envelope model from a visual mesh, by selecting certain vertices of said visual grid as control points, each control point being associated with a characteristic point of the model skeleton and subdividing the labeled mesh thus obtained into portions, each portion being associated with a deformation cage;
- the invention also relates to a system comprising means for storing a labeled envelope model and a skeleton model of the animals of the species, and a control unit programmed to execute a method according to the preceding method. .
- the subject of the invention is also a batch of animals, characterized in that each animal in the batch has been selected to be part of the batch on the basis of the value of an individual quantity representative of its behavior obtained by the implementation of a process according to the preceding method.
- the subject of the invention is also a kit, characterized in that it associates a batch of animals and an information recording medium comprising at least one overall quantity representative of the behavior of the animals of said batch of animals, obtained by implementation of a method according to the preceding method.
- the subject of the invention is also a method of using a batch of animals conforming to the preceding batch, in the screening of the activity of a pharmaceutical molecule, characterized in that it comprises a contacting step. of each animal in the batch with the molecule, a step of implementing a method according to the preceding method, to obtain an overall quantity representative of the behavior of the animals of said batch, a step of comparing the overall quantity obtained and a overall reference quantity, and a step of identifying the activity of the molecule on the basis of the result of the comparison step.
- the reference quantity used in the comparison step is an overall quantity previously associated with said batch of animals, or a global quantity measured on a batch of control animals.
- the invention also relates to an information recording medium, characterized in that it comprises instructions for executing a method according to the preceding quantization method, when the instructions are executed by an electronic calculator. .
- FIG. 1 is a block representation of the system for quantifying the behavior of a laboratory animal
- FIG. 2 is a schematic representation in block form of the method for quantifying the behavior of a laboratory animal implemented by the system of FIG. 1;
- Figure 3 is a label envelope model, consisting of a labeled mesh and a set of cages.
- a 3D posture corresponds to a representation of the three-dimensional morphology of the skeleton of an animal.
- the 3D posture designates a location, in position and rotation in the 3D space, anatomical groups relative to each other. Indeed, it is better not to go down to a too fine level of resolution in the structure of the skeleton, but to consider only groups of bones that are substantially fixed relative to each other whatever the posture adopted by the animal: skull, thorax, abdomen, tail, paws, vertebrae, ribs, femur, etc.
- a posture corresponds to the three-dimensional morphology of the skeleton of a mouse and a pose corresponds to the three-dimensional morphology of the outer surface (skin and / or fur) of the body of a mouse.
- an envelope of reference Env associates a Mesh mesh, presenting a high resolution, and a volume of control VC, constituting a mesh of lower resolution containing said mesh Mesh.
- Mesh mesh is the outer surface of the body of a mouse.
- the mesh is a high resolution surface, consisting of triangular facets adjacent to each other.
- a mesh has 10,000 vertices and 10,000 facets.
- the mesh can be defined by the positions of each of its vertices.
- the control volume VC results from the juxtaposition of a set of contiguous hexahedral Cm cages, controlling the deformation of the mesh Mesh.
- Each cage is defined by two faces Qn opposite characteristics.
- Each face is entirely defined by six degrees of freedom, which correspond to three position coordinates and three rotation coordinates.
- the Cm cages subdivide the Mesh mesh into several Meshm portions. Each Meshm portion is associated with a cage Cm of the VC control volume.
- the Cm cages must subdivide the Mesh mesh in a morphologically adapted manner. For example, a cage must be associated with the portion corresponding to the head of the animal, another one, the thorax, another one, the pelvis, etc.
- the mesh Mesh of the external surface of the body of a mouse is associated with nineteen cages Cm, defined by nineteen faces Qn.
- certain points of the mesh Mesh are characteristic points PCq whose utility is presented below.
- An elementary transformation on a cage Cm makes it possible to simply deform the associated Meshm mesh portion.
- the propagation of a deformation of a cage on the associated portion of mesh is performed for example by an interpolation algorithm.
- Each cage Cm being defined by a reduced number of degrees of freedom (six degrees of freedom), the laying of an animal is described by a hundred or so parameters, compared to the positions of the 10,000 vertices of the mesh Mesh underlying.
- the system 10 for quantifying the behavior of a laboratory animal comprises an off-line part 12 and an on-line part 14, as well as a reference database 16 and a database of data. analysis 18.
- the reference database 16 includes a skeleton template MSq, a envelope template labeled MEnvE and a pose space E, constructed from a set of EnvRef reference envelopes.
- the analysis database 18 comprises a table Til storing the instantaneous 3D postures, Pil (t), for each mouse i of a mouse lot I, over a determined period of time T.
- the off-line portion 12 of the system 10 includes an information module 22 of the reference database 16; and a module 24 for analyzing the data contained in the analysis database 18.
- the line portion 14 of the system 10 includes an acquisition module 32 generating a MeshV visual mesh of the outer surface of the body of an observed animal; an envelope model adjustment module 34 labeled MEnvE on the MeshV visual mesh in the pose space E, to obtain an optimized MeshEO labeled mesh; and a 3D posture reconstruction module 36, P (t), from the MeshEO optimized tag mesh and the MSq backbone model.
- the acquisition module 32 comprises a plurality of cameras 42 fitted to a transparent box 44, in which the mouse i of the batch I to be observed (Si,) is placed.
- eight optical video cameras having a resolution of 640/480 pixels, are used and arranged so that at any time t, the mouse can be seen by at least two cameras 42.
- the cameras 42 are synchronized with each other and calibrated (in particular in order to determine an enlargement parameter of each camera) so as to deliver, at each sampling instant t, a set of Ik frames.
- the module 32 comprises a means 46 for determining a so-called MeshV visual mesh corresponding to the external surface of the body of the mouse observed, from the set of Ik frames, as described for example in the 1994 article.
- Other means for determining MeshV visual mesh can be envisaged, such as depth camera devices (for example the "Microsoft Kinect").
- the means 46 can also implement an algorithm according to that described in paragraph 4.2 of reference 1 above.
- the MeshV visual mesh is a surface composed of triangular facets. For example, it has 10,000 vertices and 10,000 triangular facets.
- Intelligence module
- the information module 22 of the reference database 16 includes a first user interface 52, a means 54 for importing a skeleton template MSq, a means 56 for producing a envelope template labeled MEnvE, means for defining reference envelopes EnvRef and means 60 for determining the posespace E.
- the first user interface 52 allows a user to access the means 54, 56 or 58 to inform the database 16.
- the means 54 for importing a skeleton model MSq makes it possible to store a model of the skeleton MSq of the studied mice.
- each bone is a solid characterized by six degrees of freedom and the articulation between two bones is modeled by an angular spring.
- the model of the skeleton MSq corresponds to a neutral 3D posture of the animal, which corresponds to a neutral pose of the external surface of the animal. Skeletal features are associated with control points on the outer surface of the body of the animal.
- the means 56 for producing a labeled envelope template makes it possible to construct an envelope model labeled MEnvE from a MeshVO visual mesh corresponding to a neutral pose of a mouse.
- This MeshVO visual mesh is obtained at the output of the acquisition module 32.
- the means 56 makes it possible to label a group of vertices of the visual mesh
- the means 56 makes it possible to subdivide the mesh labeled MeshE into different MeshEm portions and to associate each portion with a cage Cm of deformation of a control volume.
- An envelope model labeled MEnvE is thus determined by association of a cage Cm and a portion MeshEm.
- the means 58 for defining reference envelopes EnvRef makes it possible to adjust the MeshE mesh of the envelope model labeled MEnvE to different MeshV visual meshes obtained at the output of the acquisition module 32. These visual meshes correspond to poses which, for the user, are characteristic of the behavior of a mouse. As many EnvRef reference envelopes are obtained.
- the means 60 for determining the posespace E also called the latent space, makes it possible to select a set of reference envelopes, EnvRef, and to combine them to define the posespace E. All the possible transformations on the cages of the control volume of an envelope, that is to say all the deformations of the underlying mesh, do not correspond to poses that can take the observed animal.
- the space E thus corresponds to the only transformations of the cages which correspond to a possible pose of the animal.
- the space E is a subspace of all the deformations of the cages.
- the space E has a reduced number of parameters relative to the number of degrees of freedom of all the cages of the control volume of the envelope. This reduction in the number of parameters reflects the constraints linking the degrees of freedom of the cages to describe a possible pose of the animal.
- the algorithm implemented by the means 60 is for example described in reference 2 above.
- the adjustment module 34 comprises an optimization algorithm capable of adjusting, in real time, in the pose space E, the envelope model labeled MEnvE, with a MeshV visual mesh obtained at the output of the acquisition module 32 .
- the adjustment means 34 makes it possible to obtain an optimized labeled mesh
- the reconstruction module 36 is capable of implementing a reconstruction algorithm making it possible, from the optimized MeshEO labeled mesh and the MSq skeleton model, to determine, in real time, the 3D posture, P N (t), of the mouse. S N observed.
- an inverse kinematic algorithm is implemented by which the MSq skeleton model is deformed taking into account the position of the PCq control points of the MeshEO optimized tag mesh.
- the analysis module 24 includes a second user interface 62 enabling a user to perform different queries on the different tables Til stored in the analysis database 18 and to analyze the extracted data.
- the different cameras 42 are calibrated, among other things, to take account of their respective magnification parameters.
- the first interface 52 is used to allow a user to enter the reference database 16.
- a step 1 10 the skeleton model MSq is stored in the database 16 using the means 54.
- a visual mesh is obtained by placing a reference mouse, characteristic of the studied species, in the box 44 and using the acquisition module 32.
- the visual mesh at the output of the acquisition module 32 corresponds to a neutral pose, the user selects this MeshVO visual mesh.
- control points PCq selects certain vertices of this MeshV visual mesh as control points PCq.
- Each control point is associated with a characteristic point of the MSq skeleton model.
- a mesh labeled MeshE is obtained.
- the operator sub-divides the mesh labeled MeshE into MeshEm portions, and associates with each portion a cage Cm of a control volume.
- this subdivision takes into account the skeleton model MSq, in particular groups of bones which are substantially fixed relative to one another, from one pose to another of the animal.
- An envelope model labeled MEnvE is obtained by associating with each cage Cm a MeshEm portion. Such a model is shown in Figure 3 for a mouse.
- a step 130 the database 16 is filled with EnvRef reference envelopes.
- the means 58 of the acquisition module 32 is again used with a reference mouse, so as to obtain different mesh meshV visual corresponding to poses which, according to the user, are representative of the poses that can take a mouse.
- the mesh of the envelope model labeled MEnvE is then deformed to be adjusted to each MeshV visual mesh obtained.
- An algorithm similar to that used by the adjustment module 34 is used by the means 58. During this adjustment, the cages of the envelope model labeled MEnvE are deformed. Each deformed envelope thus produced is saved in the database 16 as an EnvRef reference envelope.
- the method according to the invention provides a step 140 of generating a posespace E corresponding to the subset of the deformations of the cages having a meaning in terms of laying a mouse.
- the user selects a set of reference envelopes EnvRef in the reference database 16. For example 12 EnvRef reference envelopes are selected.
- This set of reference envelopes constitutes a kind of basis for the E-poses space.
- the details of the procedure for generating such a space are presented in paragraph 4 of reference 2.
- an intelligence phase of the database 18 can be started.
- a mouse S N to be observed is placed in the box 44.
- the algorithm of the determination means 46 of the acquisition module 32 is executed in order to obtain, in real time, a MeshVil (t) visual mesh, from at least two of the Ik frames acquired. at the current time t by the different cameras 42.
- the optimization algorithm of the adjustment module 34 is executed to adjust the labeled mesh of the envelope model labeled MEnvE on the MeshVil (t) visual mesh, while remaining in the space E (process called "morphing").
- a distance is calculated between the meshVil (t) visual mesh and the deformed labeled mesh.
- the distance is a function of the individual distances separating a vertex from the MeshVil (t) visual mesh, from the corresponding vertex of the deformed labeled mesh.
- a minimization procedure is used to find the values of the parameters in the space E that minimize this distance.
- an optimized mesh MeshEOil (t) is obtained for the mouse S N at time t.
- a modeling of the displacement of a control point with respect to the corresponding point of the skeleton is used in order to take account of a displacement and / or elasticity of the skin of the mouse in the pose studied with respect to the neutral pose.
- this modeling takes the form of a term proportional to the distance between the control point on the mesh and the corresponding point on the skeleton.
- the coefficient of proportionality is negative so as to penalize in the minimization procedure of the reconstruction algorithm, too large deviations between the control point of the external surface of the animal and the corresponding points of the skeleton of the animal.
- a step 180 the posture Pil (t) is stored in the analysis database 18, in a table Til labeled with the identifier i of the observed mouse and the identifier I of the batch to which it belongs.
- the various steps 150 to 180 of the method are repeated for the study of another mouse i + 1 of the same batch I or a mouse j of another batch 1 + 1, so as to lead to the creation of a table for monitoring the instant postures of these other mice.
- Step by step, the analysis database 18 is filled in for each mouse of a batch and for several batches.
- This morphological information can advantageously be obtained in parallel by equipping several boxes 44 with acquisition, optimization and reconstruction modules, 32, 34, 36.
- a user through the interface 62 of the module 24, makes adapted queries on the various tables Til of the database 18 in order to perform a quantitative study of the postures of the studied mice. For example statistical algorithms are executed to determine at least one individual size G M characteristic of the behavior of a mouse S ih or at least a global magnitude Gl characteristic of the overall behavior of the different mice of a batch I.
- the comparison of the individual magnitude Gi or global G1 obtained over a characteristic time window can be compared with a reference variable GRef.
- This reference quantity GRef can correspond to an individual or global quantity which is obtained on the batch of the studied mice, but during a previous period (so as to study the evolution in the time of this size), or on a batch of control mice belonging to the same species as the mice studied, but not exposed to the active molecule for example.
- the animals studied belong to the same strain of animals, a strain being defined, within a species, by a set of reproducible characters with little variability.
- These hereditary genetic characteristics are expressed (phenotype) or not (silent character), and common to all animals of said strain.
- the common character may be a motor expressed character (motor phenotype) characterized by a particular value of at least an individual magnitude Gi characteristic of the behavior of each animal of said strain.
- the adjustment procedure of a reference envelope is performed on a measurement of the outer surface of the mouse, which is a 3D mesh
- d other measurements of the outer surface of the mouse could be used.
- the acquisition module extracts a plurality of 2D contours from the outer surface of the observed animal.
- the adjustment module is then adapted to align, by deformation, the 3D mesh of the labeled envelope model on these 2D contours.
- the measurement of the external surface may alternatively come from observation means such as stereoscopic vision, depth camera, etc.
- the present method has been successfully implemented, since it has shown the possibility of performing 3D postural encoding, via a 3D registration of the position of the bone elements of an articulated skeleton from an observation of the animal.
- this method makes it possible to automate and reduce the observation time of the animals by making it automatic and more reliable. This makes it possible to meet the requirements of respect for better animal welfare.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Environmental Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Husbandry (AREA)
- Biodiversity & Conservation Biology (AREA)
- Physiology (AREA)
- Zoology (AREA)
- Clinical Laboratory Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medicinal Chemistry (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Computer Networks & Wireless Communication (AREA)
- Processing Or Creating Images (AREA)
Abstract
Ce procédé comporte les étapes : - d'acquisition d'une mesure (MeshV) d'une surface externe de l'animal observé, à partir de trames (Ik) obtenues au moyen d'un imageur (42); - d'optimisation des paramètres d'un modèle d'enveloppe étiquetée (MEnvE), associant un maillage étiqueté (MeshE) et un volume de contrôle (VC) comportant des cages (Cm) de déformation du maillage étiqueté, le maillage étiqueté comportant des points de contrôle (PCq) associés à des points caractéristiques d'un modèle du squelette (MSq) des animaux, l'étape d'optimisation consistant en un ajustement du maillage étiqueté (MeshE) à la mesure, (MeshV) pour générer un maillage étiqueté optimal (MeshEO); - de reconstruction consistant à calculer une posture (Pil(t)) à partir des points de contrôle (PCq) du maillage étiqueté optimal (MeshEO) et du modèle du squelette (MSq); et, - d'association d'une pluralité de postures pour obtenir une grandeur individuelle (Gi) représentative du comportement de l'animal.
Description
Procédé de quantification du comportement d'un animal
L'invention a pour domaine celui des procédés de quantification du comportement d'animaux, et, en particulier, d'animaux de laboratoire.
La recherche biologique a recours à des animaux de laboratoire, tels que des rongeurs et, en particulier, des rats et des souris. Ces animaux sont utilisés par exemple pour le test de molécules actives, ou pour identifier le rôle physiologique de gènes ou de protéines cibles. Dans le cadre d'études portant sur la génétique ou la protéomique, on mesure le phénotype, c'est-à-dire les conséquences comportementales et fonctionnelles de mutations génétiques (phénotypage).
Au cours de ce type de test, est suivie l'évolution du comportement d'un animal exposé à la molécule active par rapport au comportement d'un animal n'ayant pas été exposé.
Actuellement, la caractérisation du comportement d'un animal s'effectue d'abord par la simple observation du comportement de l'animal par un expérimentateur.
Une telle manière de procéder est longue et fastidieuse, en particulier lorsqu'il s'agit de suivre un lot de plusieurs animaux afin de réaliser une étude statistique. C'est essentiellement cette contrainte pratique qui limite le nombre des animaux dans les lots de test actuellement utilisés.
Cette manière de procéder, essentiellement qualitative, n'est pas objective. De nombreux biais sont introduits. Par exemple, l'expérimentateur peut introduire un biais lorsqu'il cherche à décrire le comportement d'un animal en mouvement, le mouvement d'un rongeur étant particulièrement rapide et erratique.
Cette manière de procéder s'effectue sur un laps de temps court qui n'excède pas quelque dizaine de minutes par jour. Ainsi, bien souvent, les comportements qui seraient pertinents ne sont pas observés, car survenant hors de la période d'observation.
Il a récemment été proposé plusieurs systèmes d'imagerie permettant de suivre de manière automatique et en permanence le déplacement d'un animal dans sa cage.
De tels systèmes utilisent un algorithme d'analyse d'image à deux dimensions (2D), propre à extraire le contour de l'animal à chaque instant d'échantillonnage.
Une information qui présenterait un grand intérêt dans le cadre de l'étude du comportement d'un animal de laboratoire réside dans la morphologie spatiale de son squelette, dénommée posture 3D tridimensionnelle dans ce qui suit, et l'évolution de cette posture 3D au cours du temps.
L'article de Lionel REVERET et al, "Kameleon : masse de données automatiques internes et externes pour l'étude des structures squelettiques des petits vertébrés »,
volume 28/3-4-201 1 page 309-328 du journal « traitement du signal », [Référence 1 dans le présent document] présente un algorithme de reconstruction permettant de déterminer la posture 3D instantanée d'un rat à partir de l'observation, au moyen de caméras standards et cinéradiographiques, d'un animal portant des marqueurs réels opaques aux rayons X. Chaque marqueur réel est collé sur la peau ou la fourrure de l'animal en un point de contrôle correspondant à un point caractéristique, d'extrémité ou saillant, d'un os du squelette d'un rat.
L'article de Estelle DUVEAU et al, « Cage-based Motion Recovery using Manifold Learning», second international conférence on 3D imaging, modelling, processing, visualization and transmission, Zurich Octobre 2012 [Référence 2 dans le présent document] présente un algorithme d'ajustement automatique d'un maillage 3D de référence de la surface externe du corps d'un individu humain, sur un maillage 3D bruité résultant d'une mesure directe de la surface externe du corps d'un individu humain. Lors de la procédure d'ajustement, les déformations du maillage 3D de référence sont calculées à partir de la position 3D d'un maillage de contrôle, de faible résolution, formant un volume de contrôle associé au maillage 3D de référence. L'évolution temporelle des paramètres d'ajustement du volume de contrôle permet un suivi du mouvement de l'individu.
La présente invention a pour but de répondre au besoin précité de quantification du comportement d'un animal de laboratoire.
Pour cela, l'invention a pour objet un procédé de quantification du comportement d'un animal d'une espèce donnée, caractérisé en ce qu'il comporte les étapes :
d'acquisition, consistant à acquérir une mesure d'une surface externe du corps de l'animal observé, à partir d'une pluralité de trames obtenues simultanément au moyen d'une pluralité de dispositifs imageurs ;
d'optimisation des paramètres d'un modèle d'enveloppe étiquetée, le modèle d'enveloppe étiquetée associant un maillage étiqueté et un volume de contrôle comportant des cages de déformation du maillage étiqueté, le maillage étiqueté comportant des points de contrôle associés à des points caractéristiques d'un modèle du squelette des animaux de ladite espèce, l'étape d'optimisation consistant en un ajustement du maillage étiqueté à la mesure, de manière à générer un maillage étiqueté optimal ;
de reconstruction, consistant à calculer une posture tridimensionnelle instantanée à partir des points de contrôle du maillage étiqueté optimal et du modèle du squelette ; et, - d'association d'une pluralité de postures tridimensionnelles instantanées de l'animal, déterminées en itérant les étapes d'acquisition, d'optimisation et de
reconstruction à des instants différents, pour obtenir une grandeur individuelle (Gi) représentative du comportement de l'animal.
Suivant des modes particuliers de réalisation, le procédé comporte une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles :
- les étapes sont itérées pour chaque animal d'un lot d'animaux pour obtenir une grandeur globale représentative du comportement des animaux dudit lot d'animaux ;
- il comporte une étape de comparaison de la grandeur avec une grandeur de référence ;
- le procédé comporte une étape de construction d'un modèle d'enveloppe étiqueté à partir d'un maillage visuel, en sélectionnant certains sommets dudit maillage visuel en tant que points de contrôle, chaque point de contrôle étant associé à un point caractéristique du modèle du squelette et en subdivisant le maillage étiqueté ainsi obtenu en portions, chaque portion étant associée à une cage de déformation ;
- la mesure d'une surface externe du corps de l'animal est un maillage visuel. L'invention a également pour objet un système comportant des moyens de stockage d'un modèle d'enveloppe étiquetée et d'un modèle du squelette des animaux de l'espèce, et une unité de commande programmée pour exécuter un procédé conforme au procédé précédent.
L'invention a également pour objet un lot d'animaux, caractérisé en ce que chaque animal du lot a été sélectionné pour faire partie du lot sur la base de la valeur d'une grandeur individuelle représentative de son comportement obtenue par la mise en œuvre d'un procédé conforme au procédé précédent.
L'invention a également pour objet un kit, caractérisé en ce qu'il associe un lot d'animaux et un support d'enregistrement d'informations comportant au moins une grandeur globale représentative du comportement des animaux dudit lot d'animaux, obtenue par mise en œuvre d'un procédé conforme au procédé précédent.
L'invention a également pour objet un procédé d'utilisation d'un lot d'animaux conforme au lot précédent, dans le criblage de l'activité d'une molécule pharmaceutique, caractérisé en ce qu'il comporte une étape de mise en contact de chaque animal du lot avec la molécule, une étape de mise en œuvre d'un procédé conforme au procédé précédent, pour obtenir une grandeur globale représentative du comportement des animaux dudit lot, une étape de comparaison de la grandeur globale obtenue et d'une
grandeur globale de référence, et une étape d'identification de l'activité de la molécule sur la base du résultat de l'étape de comparaison.
De préférence, la grandeur de référence utilisée à l'étape de comparaison est une grandeur globale associée précédemment audit lot d'animaux, ou une grandeur globale mesurée sur un lot d'animaux témoins.
L'invention a également pour objet, un support d'enregistrement d'informations, caractérisé en ce qu'il comporte des instructions pour l'exécution d'un procédé conforme au procédé de quantification précédent, lorsque les instructions sont exécutées par un calculateur électronique.
L'invention et ses avantages seront mieux compris à la lecture de la description qui va suivre, donnée uniquement à titre illustratif et non limitatif, d'un exemple de réalisation, faite en se référant aux dessins annexés sur lesquels :
- la figure 1 est une représentation sous forme de blocs du système de quantification du comportement d'un animal de laboratoire ;
- la figure 2 est une représentation schématique sous forme de blocs du procédé de quantification du comportement d'un animal de laboratoire mis en œuvre par le système de la figure 1 ; et,
- la figure 3 est un modèle d'enveloppe étiquetée, constituée par un maillage étiqueté et un ensemble de cages.
Dans ce qui suit, une posture 3D correspond à une représentation de la morphologie tridimensionnelle du squelette d'un animal.
Plus particulièrement, la posture 3D désigne un repérage, en position et rotation dans l'espace 3D, de groupes anatomiques les uns par rapport aux autres. En effet, il est préférable de ne pas descendre à un niveau de résolution trop fin dans la structure du squelette, mais de ne considérer que des groupes d'os qui sont sensiblement fixes les uns par rapport aux autres quelle que soit la posture adoptée par l'animal : crâne, thorax, abdomen, queue, pattes, vertèbres, côtes, fémur, etc.
Dans ce qui suit, une posture correspond à la morphologie tridimensionnelle du squelette d'une souris et une pose correspond à la morphologie tridimensionnelle de la surface externe (peau et/ou fourrure) du corps d'une souris.
De manière générale, une enveloppe de référence Env associe un maillage Mesh, présentant une résolution élevée, et un volume de contrôle VC, constituant un maillage de plus faible résolution contenant ledit maillage Mesh.
Le maillage Mesh correspond à la surface externe du corps d'une souris. Le maillage est une surface de haute résolution, constituée de facettes triangulaires adjacentes les unes aux autres. Par exemple un maillage comporte 10 000 sommets et 10 000 facettes. Le maillage peut être défini par les positions de chacun de ses sommets.
Le volume de contrôle VC résulte de la juxtaposition d'un ensemble de cages Cm hexaédriques contigues, contrôlant la déformation du maillage Mesh. Chaque cage est définie par deux faces Qn caractéristiques opposées. Chaque face est entièrement définie par six degrés de liberté, qui correspondent à trois coordonnées de position et à trois coordonnées de rotation.
Les cages Cm subdivisent le maillage Mesh en plusieurs portions Meshm. Chaque portion Meshm est associée à une cage Cm du volume de contrôle VC. Les cages Cm doivent subdiviser le maillage Mesh d'une manière adaptée sur le plan morphologique. Par exemple une cage doit être associée à la portion correspondant à la tête de l'animal, une autre, au thorax, une autre encore, au bassin, etc.
Comme représenté sur la figure 3, le maillage Mesh de la surface externe du corps d'une souris est associé à dix-neuf cages Cm, définies par dix-neuf faces Qn. Comme cela est également représenté sur la figure 3, certains points du maillage Mesh sont des points caractéristiques PCq dont l'utilité est présentée ci-dessous.
Une transformation élémentaire sur une cage Cm (déplacement et/ou rotation d'une ou de ses deux faces caractéristiques) permet de déformer simplement la portion de maillage Meshm associée. La propagation d'une déformation d'une cage sur la portion de maillage associée s'effectue par exemple par un algorithme d'interpolation.
Chaque cage Cm étant définie par un nombre réduit de degré de liberté (six degrés de liberté), la pose d'un animal est décrite par une centaine de paramètres environ, par comparaison aux positions des 10 000 sommets du maillage Mesh sous- jacent.
La manipulation de ce nombre réduit de paramètres permet un ajustement rapide du maillage Mesh de l'enveloppe, au maillage visuel dérivé d'un animal observé. SYSTEME
En se référant à la figure 1 , le système 10 de quantification du comportement d'un animal de laboratoire comporte une partie hors ligne 12 et une partie en ligne 14, ainsi qu'une base de donnée de référence 16 et une base de données d'analyse 18.
La base de données de référence 16 comporte un modèle de squelette MSq, un modèle d'enveloppe étiquetée MEnvE et un espace des poses E, élaboré à partir d'un ensemble d'enveloppes de référence EnvRef.
La base de données d'analyse 18 comporte une table Til stockant les postures 3D instantanée, Pil(t), pour chaque souris i d'un lot I de souris, sur une période de temps déterminée T.
La partie hors ligne 12 du système 10 comporte un module 22 de renseignement de la base de données de référence 16 ; et un module 24 d'analyse des données contenues dans la base de données d'analyse 18.
La partie en ligne 14 du système 10 comporte un module 32 d'acquisition générant un maillage visuel MeshV de la surface externe du corps d'un animal observé ; un module 34 d'ajustement du modèle d'enveloppe étiquetée MEnvE sur le maillage visuel MeshV dans l'espace des poses E, pour obtenir un maillage étiqueté optimisé MeshEO ; et un module 36 de reconstruction d'une posture 3D, P(t), à partir du maillage étiqueté optimisé MeshEO et du modèle de squelette MSq.
Module d'acquisition
Le module 32 d'acquisition comporte une pluralité de caméras 42 équipant une boîte 44 transparente, dans laquelle est placée la souris i du lot I à observer (Si,).
Par exemple, huit caméras vidéo optiques, ayant une résolution de 640/480 pixels, sont utilisées et disposées de manière à ce qu'à chaque instant t, la souris puisse être vue par au moins deux caméras 42.
Les caméras 42 sont synchronisées entre elles et calibrées (en particulier afin de déterminer un paramètre d'agrandissement de chaque caméra) de manière à délivrer, à chaque instant t d'échantillonnage, un ensemble de trames Ik.
Le module 32 comporte un moyen 46 de détermination d'un maillage dit visuel MeshV correspondant à la surface externe du corps de la souris observée, à partir de l'ensemble de trames Ik, comme décrit par exemple dans l'article de 1994. de A. LAURENTINI : The visual hull concept for silhouette-based image understanding (IEEE Transactions on Pattern Analysis and Machine Intelligence 16, 2, 150-162). D'autres moyens de détermination du maillage visuel MeshV peuvent être envisagés, tels que des dispositifs de caméra de profondeur (par exemple la « Microsoft Kinect »). Le moyen 46 peut aussi mettre en œuvre un algorithme conforme à celui décrit au paragraphe 4.2 de la référence 1 ci-dessus.
Le maillage visuel MeshV est une surface constituée de facettes triangulaires. Il comporte par exemple 10 000 sommets et 10 000 facettes triangulaires.
Module de renseignement
Le module 22 de renseignement de la base de données de référence 16 comporte une première interface utilisateur 52, un moyen 54 d'importation d'un modèle de squelette MSq, un moyen 56 d'élaboration d'un modèle d'enveloppe étiquetée MEnvE, un moyen 58 de définition d'enveloppes de référence EnvRef et un moyen 60 de déterminer l'espace des poses E.
La première interface utilisateur 52 permet à un utilisateur d'accéder aux moyens 54, 56 ou 58 pour renseigner la base de données 16.
Le moyen 54 d'importation d'un modèle de squelette MSq permet de stocker un modèle du squelette MSq des souris étudiées. Dans le modèle du squelette MSq utilisé, chaque os est un solide caractérisé par six degrés de liberté et l'articulation entre deux os est modélisée par un ressort angulaire. Le modèle du squelette MSq correspond à une posture 3D neutre de l'animal, celle-ci correspondant à une pose neutre de la surface externe de l'animal. Des points caractéristiques du squelette sont associés à des points de contrôle sur la surface externe du corps de l'animal.
Le moyen 56 d'élaboration d'un modèle d'enveloppe étiquetée permet de construire un modèle d'enveloppe étiquetée MEnvE à partir d'un maillage visuel MeshVO correspondant à une pose neutre d'une souris.
Ce maillage visuel MeshVO est obtenu en sortie du module d'acquisition 32.
Le moyen 56 permet d'étiqueter un groupe de sommets du maillage visuel
MeshVO en tant que points de contrôle PCq associés au point caractéristique du squelette MSq. Un maillage étiqueté MeshE est ainsi obtenu.
Le moyen 56 permet de subdiviser le maillage étiqueté MeshE en différentes portions MeshEm et d'associer chaque portion à une cage Cm de déformation d'un volume de contrôle.
Un modèle d'enveloppe étiquetée MEnvE est ainsi déterminé par association d'une cage Cm et d'une portion MeshEm.
Le moyen 58 de définition d'enveloppes de référence EnvRef permet d'ajuster le maillage MeshE du modèle d'enveloppe étiquetée MEnvE à différents maillages visuels MeshV obtenus en sortie du module d'acquisition 32. Ces maillages visuels correspondent à des poses qui, pour l'utilisateur, sont caractéristiques du comportement d'une souris. Autant d'enveloppes de référence EnvRef sont obtenues.
Le moyen 60 de détermination de l'espace des poses E, aussi dénommé espace latent, permet de sélectionner un ensemble d'enveloppes de référence, EnvRef, et à les combiner pour définir l'espace des poses E.
Toutes les transformations possibles sur les cages du volume de contrôle d'une enveloppe, c'est-à-dire toutes les déformations du maillage sous-jacent, ne correspondent pas à des poses que peut prendre l'animal observé.
L'espace E correspond ainsi aux seules transformations des cages qui correspondent à une pose possible de l'animal. L'espace E est un sous-espace de l'ensemble des déformations des cages.
L'espace E présente un nombre réduit de paramètres par rapport au nombre de degrés de liberté de l'ensemble des cages du volume de contrôle de l'enveloppe. Cette réduction du nombre des paramètres traduit les contraintes liant les degrés de liberté des cages pour qu'elles décrivent une pose possible de l'animal.
L'algorithme mis en œuvre par le moyen 60 est par exemple décrit dans la référence 2 ci-dessus.
Module d'ajustement
La module d'ajustement 34 comporte un algorithme d'optimisation propre à ajuster, en temps réel, dans l'espace des poses E, le modèle d'enveloppe étiquetée MEnvE, à un maillage visuel MeshV obtenu en sortie du module d'acquisition 32.
Les détails de l'algorithme d'optimisation mis en œuvre sont présentés au point 5 de la référence 2 mentionnée ci-dessus.
Le moyen d'ajustement 34 permet d'obtenir un maillage étiqueté optimisé
MeshEO.
Module de reconstruction
Le module de reconstruction 36 est propre à mettre en œuvre un algorithme de reconstruction permettant à partir du maillage étiqueté optimisé MeshEO et du modèle de squelette MSq, de déterminer, en temps réel, la posture 3D, PN(t), de la souris SN observée.
En particulier, est mis en œuvre un algorithme dit de cinématique inverse par lequel le modèle de squelette MSq est déformé en tenant compte de la position des points de contrôle PCq du maillage étiqueté optimisé MeshEO.
Les détails d'un exemple d'un tel algorithme de reconstruction sont donnés au point 5.2 de la référence 1 ci-dessus, en substituant aux marqueurs externes utilisés dans cette référence, les points de contrôle PCq du maillage étiqueté optimisé MeshEO.
Module d'analyse
Le module d'analyse 24 comporte une seconde interface utilisateur 62 permettant à un utilisateur d'effectuer différentes requêtes sur les différentes tables Til stockées dans la base de données d'analyse 18 et d'analyser les données extraites.
PROCEDE
Dans une étape 100, les différentes caméras 42 sont calibrées, entre autre pour tenir compte de leurs paramètres d'agrandissement respectifs.
La première interface 52 est utilisée pour permettre à un utilisateur de renseigner la base de données de référence 16.
Dans une étape 1 10, le modèle de squelette MSq est stocké dans la base de données 16 en utilisant le moyen 54.
Dans une étape 120 de construction d'un modèle d'enveloppe étiqueté MEnvE, un maillage visuel est obtenu en plaçant une souris de référence, caractéristique de l'espèce étudiée, dans la boîte 44 et en utilisant le module d'acquisition 32. Lorsque le maillage visuel en sortie du module d'acquisition 32 correspond à une pose neutre, l'utilisateur sélectionne ce maillage visuel MeshVO.
En utilisant le moyen 56, l'utilisateur sélectionne ensuite certains sommets de ce maillage visuel MeshV en tant que points de contrôle PCq. Chaque point de contrôle est associé à un point caractéristique du modèle de squelette MSq. Un maillage étiqueté MeshE est obtenu.
L'opérateur subdivise le maillage étiqueté MeshE en portions MeshEm, et associe à chaque portion une cage Cm d'un volume de contrôle. Avantageusement cette subdivision s'effectue en tenant compte du modèle de squelette MSq, en particulier, des groupes d'os qui sont sensiblement fixes les uns par rapport aux autres, d'une pose à l'autre de l'animal.
Un modèle d'enveloppe étiquetée MEnvE est obtenu en associant à chaque cage Cm une portion MeshEm. Un tel modèle est représenté sur la figure 3 pour une souris.
Dans une étape 130, la base de données 16 est renseignée avec des enveloppes de référence EnvRef.
Pour ce faire, le moyen 58 du module d'acquisition 32 est à nouveau utilisé avec une souris de référence, de manière à obtenir différents maillages visuels MeshV correspondant à des poses qui, selon l'utilisateur, sont représentatifs des poses que peut prendre une souris.
Le maillage du modèle d'enveloppe étiquetée MEnvE est alors déformé pour être ajusté à chaque maillage visuel MeshV obtenu. Un algorithme similaire à celui utilisé par
le module d'ajustement 34 est utilisé par le moyen 58. Lors de cet ajustement, les cages du modèle d'enveloppe étiquetée MEnvE sont déformées. Chaque enveloppe déformée ainsi réalisée est sauvegardée dans la base de données 16 en tant qu'enveloppe de référence EnvRef.
Le procédé selon l'invention prévoit une étape 140 de génération d'un espace des poses E correspondant au sous-ensemble des déformations des cages ayant un sens en termes de pose d'une souris.
Pour créer l'espace E, l'utilisateur sélectionne un ensemble d'enveloppes de référence EnvRef dans la base de données de référence 16. Par exemple 12 enveloppes de référence EnvRef sont sélectionnées.
Cet ensemble d'enveloppes de référence constitue une sorte de base de l'espace des poses E. Le détail de la procédure de génération d'un tel espace est présenté dans le paragraphe 4 de la référence 2.
Une fois réalisées ces différentes étapes de configuration du système 10, une phase de renseignement de la base de données 18 peut être lancée.
Une souris SN à observer est placée dans la boîte 44.
Dans une étape 150 d'acquisition, l'algorithme du moyen de détermination 46 du module d'acquisition 32 est exécuté pour obtenir, en temps réel, un maillage visuel MeshVil(t), à partir d'au moins deux des trames Ik acquises à l'instant courant t par les différentes caméras 42.
Dans une étape 160 d'optimisation, l'algorithme d'optimisation du module d'ajustement 34 est exécuté pour ajuster le maillage étiqueté du modèle d'enveloppe étiqueté MEnvE sur le maillage visuel MeshVil(t), tout en restant dans l'espace E (processus dit de « morphing »).
Plus précisément, au cours de cet ajustement, les valeurs des paramètres dans l'espace E du modèle d'enveloppe étiqueté MEnvE sont modifiées. Cela entraîne une déformation du maillage étiqueté sous-jacent.
Une distance est calculée entre le maillage visuel MeshVil(t) et le maillage étiqueté déformé. La distance est une fonction des distances individuelles séparant un sommet du maillage visuel MeshVil(t), du sommet correspondant du maillage étiqueté déformé.
Une procédure de minimisation est utilisée pour rechercher les valeurs des paramètres dans l'espace E qui minimisent cette distance.
A la fin de l'exécution de l'algorithme du module 34, un maillage étiqueté optimisé MeshEOil(t) est obtenu pour la souris SN à l'instant t.
Dans une étape 170 de reconstruction, l'algorithme de reconstruction du module
36 est exécuté pour obtenir, en temps réel, une posture 3D, Pil(t), en déformant le modèle
de squelette MSq, contenu dans la base de données 16, de manière à l'adapter au maillage étiqueté optimisé MeshEOil(t). Cette adaptation s'effectue en tenant compte des points de contrôle PCq que comporte le maillage MeshEOil(t).
En particulier, est utilisée une modélisation du déplacement d'un point de contrôle par rapport au point correspondant du squelette afin de tenir compte d'un déplacement et/ou d'une élasticité de la peau de la souris dans la pose étudiée par rapport à la pose neutre. Dans le mode de réalisation actuellement envisagé, cette modélisation prend la forme d'un terme proportionnel à la distance entre le point de contrôle sur le maillage et le point correspondant sur le squelette. Le coefficient de proportionnalité est négatif de manière à pénaliser dans la procédure de minimisation de l'algorithme de reconstruction, les écarts trop importants entre le point de contrôle de la surface externe de l'animal et les points correspondant du squelette de l'animal.
Ainsi, une posture 3D à l'instant courant t est obtenue pour la souris SN. Cette posture est référencée par Pil(t).
Dans une étape 180, la posture Pil(t) est stockée dans la base de données d'analyse 18, dans une table Til étiquetée avec l'identifiant i de la souris observée et l'identifiant I du lot auquel elle appartient.
Les différentes étapes 150 à 180 ayant conduit à la détermination d'une posture 3D sont répétées à un instant ultérieur t+At de manière à obtenir, de la même souris Sy, la posture 3D à cet autre instant Pil(t+ At). La table Til de la base de données d'analyse est ainsi renseignée périodiquement avec la posture instantanée de la souris.
Les différentes étapes 150 à 180 du procédé sont répétées pour l'étude d'une autre souris i+1 du même lot I ou d'une souris j d'un autre lot 1+1 , de manière à conduire à la création d'une table de suivi des postures instantanées de ces autres souris.
De proche en proche, la base de données d'analyse 18 est renseignée pour chaque souris d'un lot et pour plusieurs lots.
Ces informations morphologiques peuvent avantageusement être obtenues parallèlement en équipant plusieurs boîtes 44 de modules d'acquisition, d'optimisation et de reconstruction, 32, 34, 36.
Dans une étape d'analyse 190, un utilisateur, au travers de l'interface 62 du module 24, effectue des requêtes adaptées sur les différentes tables Til de la base de données 18 afin de réaliser une étude quantitative des postures des souris étudiées. Par exemple des algorithmes statistiques sont exécutés pour déterminer au moins une grandeur individuelle GM caractéristique du comportement d'une souris Sih ou au moins une grandeur globale Gl caractéristique du comportement global des différentes souris d'un lot I.
En particulier, la comparaison de la grandeur individuelle Gi ou globale Gl obtenue sur une fenêtre de temps caractéristique peut être comparée à une grandeur de référence GRef. Cette grandeur de référence GRef peut correspondre à une grandeur individuelle ou globale qui est obtenue sur le lot des souris étudiés, mais au cours d'une période antérieure (de manière à étudier l'évolution dans le temps de cette grandeur), ou sur un lot de souris témoins appartenant à la même espèce que les souris étudiées, mais n'ayant pas été exposées à la molécule active par exemple.
Avantageusement, les animaux étudiés appartiennent à une même souche d'animaux, une souche étant définie, au sein d'une espèce, par un ensemble de caractères reproductibles avec peu de variabilité. Ces caractère génétiques héréditaires sont exprimés (phénotype) ou non (caractère silencieux), et communs à tous les animaux de ladite souche. En particulier, le caractère commun peut être un caractère exprimé moteur (phénotype moteur) se caractérisant par une valeur particulière d'au moins une grandeur individuelle Gi caractéristique du comportement de chaque animal de ladite souche.
En variante, alors que, dans le mode de réalisation préféré décrit ci-dessus en détail, la procédure d'ajustement d'une enveloppe de référence s'effectue sur une mesure de la surface externe de la souris qui est un maillage 3D, d'autres mesures de la surface externe de la souris pourraient être utilisées. Par exemple, le module d'acquisition extrait une pluralité de contours 2D de la surface externe de l'animal observé. Le module d'ajustement est alors adapté pour aligner, par déformation, le maillage 3D du modèle d'enveloppe étiquetée sur ces contours 2D.
La mesure de la surface externe peut provenir en variante de moyens d'observation du type vision stéréoscopique, caméra de profondeur, etc.
Le présent procédé a été mise en œuvre avec succès, puisqu'il a montré la possibilité de réaliser l'encodage postural 3D, via un repérage 3D de la position des éléments osseux d'un squelette articulé à partir d'une observation de l'animal.
De manière générale, ce procédé permet d'automatiser et de réduire le temps d'observation des animaux en le rendant automatique et plus fiable. Ceci permet de répondre aux exigences du respect d'un meilleur bien-être animal.
Claims
1 . - Procédé de quantification du comportement d'un animal d'une espèce donnée, caractérisé en ce qu'il comporte les étapes :
- d'acquisition (150), consistant à acquérir une mesure (MeshV) d'une surface externe du corps de l'animal observé, à partir d'une pluralité de trames (Ik) obtenues simultanément au moyen d'une pluralité de dispositifs imageurs (42) ;
d'optimisation (160) des paramètres d'un modèle d'enveloppe étiquetée (MEnvE), ledit modèle d'enveloppe étiquetée associant un maillage étiqueté (MeshE) et un volume de contrôle (VC) comportant des cages (Cm) de déformation dudit maillage étiqueté, ledit maillage étiqueté comportant des points de contrôle (PCq) associés à des points caractéristiques d'un modèle du squelette (MSq) des animaux de ladite espèce, ladite étape d'optimisation consistant en un ajustement dudit maillage étiqueté (MeshE) à ladite mesure,(MeshV) de manière à générer un maillage étiqueté optimal (MeshEO) ;
- de reconstruction (170), consistant à calculer une posture tridimensionnelle instantanée (Pil(t)) à partir des points de contrôle (PCq) du maillage étiqueté optimal (MeshEO) et du modèle du squelette (MSq) ; et,
d'association d'une pluralité de postures tridimensionnelles instantanées de l'animal, déterminées en itérant les étapes d'acquisition, d'optimisation et de reconstruction à des instants différents, pour obtenir une grandeur individuelle (Gi) représentative du comportement de l'animal.
2. - Procédé selon la revendication 1 , dans lequel les étapes sont itérées pour chaque animal (i) d'un lot (I) d'animaux pour obtenir une grandeur globale (GI) représentative du comportement des animaux dudit lot d'animaux.
3. - Procédé selon la revendication 2, comportant une étape de comparaison de la grandeur (Gi, GI) avec une grandeur de référence (GRef).
4.- Procédé selon l'une quelconque des revendications précédentes, comportant une étape (120) de construction d'un modèle d'enveloppe étiqueté (MEnvE) à partir d'un maillage visuel, en sélectionnant certains sommets dudit maillage visuel en tant que points de contrôle, chaque point de contrôle étant associé à un point caractéristique du modèle du squelette (MSq) et en subdivisant le maillage étiqueté (MeshE) ainsi obtenu en portions, chaque portion étant associée à une cage (Cm) de déformation.
5. - Procédé selon l'une quelconque des revendications 1 à 4, dans lequel la mesure d'une surface externe du corps de l'animal est un maillage visuel (MeshV).
6. - Système, caractérisé en ce qu'il comporte des moyens de stockage d'un modèle d'enveloppe étiquetée (MEnvE) et d'un modèle du squelette (MSq) des animaux de ladite espèce, et une unité de commande programmée pour exécuter un procédé conforme à l'une quelconque des revendications 1 à 5.
7. - Lot d'animaux, caractérisé en ce qu'un comportement de chaque animal dudit lot est caractérisé par au moins une grandeur individuelle (Gi) représentative de comportement dudit animal, qui est obtenue par la mise en œuvre du procédé selon l'une quelconque des revendications 1 à 5 et qui correspond à une grandeur de référence (GRef) associée audit lot.
8.- Kit, caractérisé en ce qu'il associe un lot d'animaux conforme à la revendication 7 et un support d'enregistrement d'informations, ledit support d'enregistrement d'informations comportant ladite grandeur de référence (GRef) associée audit lot.
9.- Procédé d'utilisation d'un kit conforme à la revendication 8 dans le criblage de l'activité d'une molécule pharmaceutique, caractérisé en ce qu'il comporte une étape de mise en contact de chaque animal du lot dudit kit avec la molécule, une étape de mise en œuvre du procédé selon la revendication 2 pour obtenir une grandeur globale (GI) représentative du comportement des animaux dudit lot après mise en contact avec la molécule, une étape de comparaison de la grandeur globale obtenue avec une grandeur globale de référence (GRef), et une étape d'identification de l'activité de la molécule sur la base du résultat de l'étape de comparaison.
10. - Procédé selon la revendication 9, dans lequel la grandeur globale de référence (GRef) utilisée à l'étape de comparaison est la grandeur de référence stockée sur le support d'enregistrement d'informations dudit kit ou une grandeur globale mesurée sur un lot d'animaux témoins.
1 1 . - Support d'enregistrement d'informations, caractérisé en ce qu'il comporte des instructions pour l'exécution d'un procédé conforme à l'une quelconque des
revendications 1 à 5, lorsque les instructions sont exécutées par un calculateur électronique.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1358835A FR3010623B1 (fr) | 2013-09-13 | 2013-09-13 | Procede de quantification du comportement d'un animal |
| PCT/EP2014/069415 WO2015036495A1 (fr) | 2013-09-13 | 2014-09-11 | Procédé de quantification du comportement d'un animal |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| EP3043710A1 true EP3043710A1 (fr) | 2016-07-20 |
Family
ID=50231251
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| EP14761892.0A Withdrawn EP3043710A1 (fr) | 2013-09-13 | 2014-09-11 | Procédé de quantification du comportement d'un animal |
Country Status (3)
| Country | Link |
|---|---|
| EP (1) | EP3043710A1 (fr) |
| FR (1) | FR3010623B1 (fr) |
| WO (1) | WO2015036495A1 (fr) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022142458A1 (fr) | 2020-12-31 | 2022-07-07 | 武汉联影生命科学仪器有限公司 | Procédé, appareil et système de balayage d'animal |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2446853A1 (fr) * | 2001-05-15 | 2002-11-21 | Psychogenics Inc. | Systemes et procedes de controle informatique du comportement |
-
2013
- 2013-09-13 FR FR1358835A patent/FR3010623B1/fr not_active Expired - Fee Related
-
2014
- 2014-09-11 EP EP14761892.0A patent/EP3043710A1/fr not_active Withdrawn
- 2014-09-11 WO PCT/EP2014/069415 patent/WO2015036495A1/fr not_active Ceased
Non-Patent Citations (1)
| Title |
|---|
| See references of WO2015036495A1 * |
Also Published As
| Publication number | Publication date |
|---|---|
| FR3010623A1 (fr) | 2015-03-20 |
| WO2015036495A1 (fr) | 2015-03-19 |
| FR3010623B1 (fr) | 2015-10-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110662484B (zh) | 用于全身测量结果提取的系统和方法 | |
| Panagiotou et al. | Generating elevation surface from a single RGB remotely sensed image using deep learning | |
| CN111462120A (zh) | 一种基于语义分割模型缺陷检测方法、装置、介质及设备 | |
| JP7518490B2 (ja) | 圃場の育種データ収集装置、育種における特徴解析装置、圃場の育種データ収集方法、プログラム、および記録媒体 | |
| CN116519106B (zh) | 一种用于测定生猪体重的方法、装置、存储介质和设备 | |
| Jiang et al. | Cotton3DGaussians: Multiview 3D Gaussian Splatting for boll mapping and plant architecture analysis | |
| CN117275613A (zh) | 一种面膜配方的调配方法 | |
| CN116229001A (zh) | 一种基于空间熵的城市三维数字地图生成方法及系统 | |
| Bernabé et al. | Portability and acceleration of deep learning inferences to detect rapid earthquake damage from vhr remote sensing images using intel openvino toolkit | |
| CN117912056A (zh) | 基于计算机视觉的养殖鱼类全周期识别方法 | |
| US20230169662A1 (en) | System and method for generating a morphological atlas of an embryo | |
| Li et al. | Dessie: Disentanglement for articulated 3d horse shape and pose estimation from images | |
| Yin et al. | Automated follicle counting system (AFCS) using YOLO-based object detection algorithm and its application in the POI model | |
| Ramachandran et al. | Tiny Criss-cross network for segmenting paddy panicles using aerial images | |
| Wang et al. | Using UAVRS and deep learning to conduct resource surveys of threatened Tibetan medicinal plants in the Qinghai-Tibet Plateau | |
| Rodrıguez-Malpica et al. | Tree crown prediction of spruce tree using a machine-learning-based hybrid image processing method | |
| EP3043710A1 (fr) | Procédé de quantification du comportement d'un animal | |
| CN114264355A (zh) | 重量检测方法、装置、电子设备及存储介质 | |
| Zhao et al. | Automatic waterfowl and habitat detection using drone imagery and deep learning | |
| CN116012528B (zh) | 一种大白菜卷曲性状量化的三维重建方法 | |
| Mirka et al. | A photogrammetric approach to the estimation of distance to animals in camera trap images | |
| Zhou et al. | High-throughput wheat seedling phenotyping via UAV-based semantic segmentation and ground sample distance driven pixel-to-area mapping | |
| Li et al. | Object-Centric 3D Gaussian Splatting for Strawberry Plant Reconstruction and Phenotyping | |
| Au et al. | ARISGAN: Extreme super-resolution of arctic surface imagery using generative adversarial networks | |
| Siatou et al. | Adaption of Imaging Techniques for Monitoring Cultural Heritage Objects |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
| 17P | Request for examination filed |
Effective date: 20160316 |
|
| AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
| AX | Request for extension of the european patent |
Extension state: BA ME |
|
| DAX | Request for extension of the european patent (deleted) | ||
| STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN |
|
| 18D | Application deemed to be withdrawn |
Effective date: 20200603 |