Brevet FR2927442A1 — publié le 20090814

- 1- " Procédé pour déterminer un taux de transformation local d'un objet d'intérêt." La présente invention se rapporte à un procédé pour déterminer un taux de transformation local d'un objet d'intérêt présent dans une zone d'observation, dans lequel procédé on acquiert une pluralité d'images de la zone d'observation au moyen d'une caméra. Elle trouve une application particulièrement intéressante mais non limitative dans le domaine de io l'analyse comportementale par vision artificielle. Plus précisément, la présente invention peut s'appliquer à des systèmes dans lesquels on dispose d'au moins une caméra disposée en hauteur dans un endroit stratégique d'un magasin où l'on est capable de visualiser des déplacements de personnes. Le flux vidéo de la caméra est traité de façon notamment à 15 réaliser des analyses comportementales. Un but de la présente invention est d'améliorer la disposition d'objets dans les rayons d'un magasin. La présente invention a aussi pour but un nouvel outil d'analyse comportementale. Un autre but de l'invention est un procédé pour estimer un taux de transformation d'un produit de manière 20 précise et locale. On atteint au moins l'un des buts précités avec un procédé pour déterminer un taux de transformation local d'un objet d'intérêt présent dans une zone d'observation. Cet objet d'intérêt peut être physiquement présent ou 25 simplement représenté, par exemple sous forme d'une image, d'un panneau publicitaire,... Dans ce procédé, on acquiert une pluralité d'images de la zone d'observation au moyen d'une caméra. Selon l'invention, pour chaque image acquise on réalise les étapes suivantes : - détection de zones de pixels se détachant du fond fixe de l'image, 30 - détection d'individus en déplacement dans la zone d'observation, - pour chaque individu détecté, détermination de la position occupée par cet individu, et - pour l'ensemble des positions déterminées, on minimise la distance globale entre les positions ainsi déterminées dans l'image en cours et 35 des positions déterminées lors d'une image précédente. - 2- Avantageusement, pour un intervalle de temps donné, on détermine la trajectoire de chaque individu ayant traversé la zone d'observation, et on détermine un taux de transformation comme étant le rapport entre le nombre de trajectoires ainsi déterminées et une valeur d'intérêt qui est fonction dudit objet d'intérêt. L'homme du métier comprendra aisément que ce rapport peut être défini dans un sens (nombre de trajectoires/valeur d'intérêt) ou dans le sens inverse (valeur d'intérêt/nombre de trajectoires). Avec le procédé selon l'invention, on peut déterminer de façon très précise un taux de transformation local : la zone d'observation est une zone io qui est judicieusement définie autour notamment de l'objet d'intérêt. Cette zone d'observation peut avantageusement correspondre au champ de vision de la caméra. En d'autres termes, pour réaliser le taux de transformation, on réalise un processus de suivi de déplacements, puis on en déduit des trajectoires. 15 Ensuite on compare le nombre de ces trajectoires avec une valeur d'intérêt. Selon l'invention, cette valeur d'intérêt peut être le nombre d'individus dont la trajectoire présente un temps d'arrêt de durée prédéterminée devant l'objet d'intérêt. Cette valeur d'intérêt représente en fait l'intérêt que les passants ont pu avoir sur l'objet d'intérêt. Le fait de marquer un temps d'arrêt 20 signifie que l'individu a porté un certain intérêt à l'objet d'intérêt. A titre d'exemple non limitatif, l'objet d'intérêt peut être : - un outil de communication tel qu'un écran plasma, un élément d'affichage quelconque...placé dans un lieu fermé comme un magasin ou non (à l'extérieur) ; 25 - une vitrine de présentation de produit, en interne ou en externe ; - l'entrée d'un espace tel qu'un magasin : la valeur d'intérêt pouvant alors être le nombre de personnes ayant franchi une porte d'entrée. Selon l'invention, lorsque plusieurs objets d'intérêt sont disponibles dans la zone d'observation ou ailleurs, cette valeur d'intérêt peut également 30 être le nombre d'objets d'intérêt ayant été détectés par un appareil de détection de sortie d'objets d'intérêt. Avantageusement cet appareil de détection peut être un lecteur infrarouge présent dans les caissiers de magasins et relié à un serveur de gestion de stock par exemple. En fait, dans ce cas, la valeur d'intérêt représente le nombre d'objets d'intérêt vendus par 35 le magasin dans l'intervalle de temps donné. On peut aussi imaginer la détermination du rapport entre le nombre d'individus dont la trajectoire présente un temps d'arrêt de durée prédéterminée devant l'objet d'intérêt et le nombre d'objets d'intérêt ayant été détectés par l'appareil de détection de sortie d'objets d'intérêt. La présente invention est notamment basée sur un processus de suivi de déplacements d'individus dans la zone d'observation. Le processus de suivi fait intervenir une détection d'individus image par image. Plus deux détections sont proches, plus elles sont susceptibles de correspondre à un même io individu. Ainsi la trajectoire d'un individu est formée par un ensemble de positions déterminées à travers des images successives et appariées ensemble selon le critère de la minimisation globale. Selon une caractéristique avantageuse de l'invention, la minimisation de 15 la distance globale consiste à réaliser les étapes suivantes : pour chaque position d'individu déterminée, on évalue la distance par rapport aux autres positions d'individus, définissant ainsi une matrice de distances que l'on convertit en matrice de similitudes, puis on normalise des vecteurs propres de cette matrice de similitudes. En d'autres termes, on réalise en fait une 20 méthode à base de décomposition en valeurs singulières. La normalisation des vecteurs propres de cette matrice de similitudes permet d'identifier le meilleur appariement global de l'ensemble des positions. Selon l'invention, la distance entre deux positions est définie en prenant 25 en compte : - une distance géométrique 3D : chaque position est identifiée dans un repère 3D fixe, il est donc possible d'évaluer la distance en mètres séparant deux positions ; - une distance géométrique 2D avec prise en compte de recouvrement 30 chaque position est identifiée dans le plan de l'image, il est donc possible d'évaluer la distance en pixels séparant deux positions dans cette représentation, on peut notamment utiliser la distance dite de Hausdorf, et/ou utiliser des critères de recouvrement de positions ; - une distance colorimétrique : pour chaque position, il est possible 35 d'établir une signature à partir de ses histogrammes de couleur ; la - 4- mesure de corrélation entre les signatures de deux positions permet d'en déduire une distance d'apparence des deux positions ; - une cohérence temporelle : la date de détection de deux positions peut être intégrée dans la mesure de distance pour apparier de préférence deux positons proches dans le temps. L'invention est notamment remarquable par le fait qu'on effectue un appariement global sur l'ensemble des individus détectées, en utilisant des notions riches de distance entre individus pouvant être basées sur des notions de position géométrique, dimensions, signature d'apparence, cohérence io temporelle, . Selon une caractéristique avantageuse de l'invention, l'étape b) de détection des individus comprend les étapes suivantes : - détection de têtes des individus dans chaque zone de pixels, et - dénombrement des individus dans chaque zone de pixels. 15 En particulier, l'étape de détection de têtes comprend les étapes suivantes - application d'un filtre de Canny sur l'image par rapport aux zones de pixels de façon à générer une image de contours des zones de pixels, 20 - application d'une transformée en distance de façon à élaborer une carte de distances, -réalisation d'une convolution entre la carte de distance et un modèle de tête (� template �) de façon à obtenir une carte de convolution, et 25 -application d'un algorithme de ligne de partage des eaux (� watershed �) de façon à détecter des têtes d'individus. L'utilisation d'algorithme de ligne de partage des eaux pour la détection de têtes d'individus permet d'obtenir les positions les plus probables de têtes dans l'image sans avoir recours à des valeurs arbitraires 30 (notamment de valeurs de seuils) ce qui limite les risques de non détection et garantit l'unicité d'une solution locale. Avantageusement, le dénombrement consiste à identifier simultanément le nombre d'individus présents et leur position respective. En d'autres termes, on réalise une segmentation de groupes. En effet, dans la 35 plupart des circonstances, les images proviennent d'une caméra ayant parfois un angle de visée oblique par rapport au sol, provoquant ainsi de nombreuses situations d'occlusion entre individus, il est alors difficile d'utiliser une recherche individu par individu. Pour cela, la présente invention propose donc une solution cherchant à détecter et dénombrer un groupe (identification simultanée du nombre de personnes présentes et de leur position respective), et non à détecter chaque individu séparément. Selon une caractéristique de l'invention, l'identification des individus est effectuée à partir d'un modèle d'apparence à trois ellipses imitant un individu. io Ce modèle d'apparence est optimisé selon les critères suivants : - minimisation du nombre d'individus, - maximisation du nombre de pixels différents du fond et appartenant aux silhouettes, - minimisation du nombre de pixels différents du fond et 15 n'appartenant pas aux ellipses, - minimisation du nombre de pixels appartenant au fond et aux ellipses, et - minimisation de la distance des contours d'un modèle 2D/3D. Ensuite, l'optimisation est réalisée au moyen d'un algorithme itératif 20 de descente de gradient. Selon une caractéristique avantageuse de l'invention, on répartit la convergence de l'algorithme itératif sur des images successives. La convergence de l'algorithme itératif vers un dénombrement optimal des individus peut prendre un temps considérable. Avec la présente invention, 25 on se contente d'une solution partielle approchée sur une image au moment de passer à l'image suivante; cette solution approchée servant de point de départ pour l'image suivante. Selon une variante avantageuse de l'invention, l'étape b) de détection des individus peut comprendre les étapes suivantes : 30 à partir d'un modèle 3D de la zone d'observation, d'un modèle 3D d'individu, et des caractéristiques liées à l'appareil d'acquisition de la pluralité d'images d'observation, on réalise une projection perspective dans l'image en cours de chaque surface élémentaire de la zone d'observation de façon à obtenir une projection 2D des modèles d'individu ; on calcule la correspondance des modèles 2D sur chaque zone de pixels de façon à obtenir une carte de densité de probabilité de présence ; et on réalise une recherche de maxima par l'algorithme de ligne de partage des eaux de façon à déterminer les individus détectés. De préférence, les caractéristiques de l'appareil d'acquisition comprennent le champ de vision, le positionnement et la distorsion. D'une façon générale, selon l'invention, préalablement à l'étape a), on réalise une étape d'extraction d'un fond de l'image de façon à déterminer un io modèle de fond fixe. L'extraction du fond fixe de l'image est obtenue en moyennant les valeurs des pixels de l'image au cours du temps. Selon l'invention, pour tenir compte des faibles variations des valeurs des pixels, on peut modéliser la valeur de chaque pixel de fond par une distribution de probabilités. Cette distribution peut être représentée par une 15 gaussienne, ou plus sobrement par une valeur moyenne et deux valeurs min et max. Par ailleurs, afin de minimiser autant que possible l'influence des objets n'appartenant pas au fond, avant de moyenner les valeurs des pixels de l'image, il est possible d'effectuer une détection instantanée de 20 mouvement par soustraction d'images successives et d'appliquer un seuil afin d'obtenir un masque correspondant aux pixels qui ne seront pas moyennés...DTD: D'autres avantages et caractéristiques de l'invention apparaîtront à 25 l'examen de la description détaillée d'un mode de mise en oeuvre nullement limitatif, et des dessins annexés, sur lesquels : La figure 1 est une vue schématique générale d'un exemple de mise en oeuvre du procédé selon l'invention, La figure 2 est une vue schématique illustrant le positionnement de 3o plusieurs individus à deux instants différents, La figure 3 est une vue schématique illustrant l'appariement des positions de la figure 2, La figure 4 est une vue générale d'un organigramme illustrant différentes étapes d'un procédé selon l'invention, - 7- Les figures 5 et 6 sont des images provenant d'une caméra de surveillance sans et avec des blobs apparents, La figure 7 est une vue d'un organigramme détaillé d'un algorithme de détection de têtes d'individus, La figure 8 illustre un modèle d'apparence à trois ellipses selon l'invention, La figure 9 est un schéma simplifié illustrant différentes étapes d'une variante d'algorithme de détection d'individus, et La figure 10 est une vue de l'algorithme décrit sur la figure 9. i0 Sur la figure 1 on voit un rayon 1 d'un magasin équipé d'une caméra de surveillance 2 disposée en hauteur. Cette caméra de surveillance est fixée au plafond et est disposée en oblique de façon à filmer un objet d'intérêt 3 dans un endroit stratégique tel qu'un carrefour de plusieurs 15 allées dans le magasin. Cette caméra peut tout aussi bien être disposée ailleurs qu'au plafond. Par exemple, elle peut être disposée en hauteur mais pas forcément au plafond. Sa position n'est pas nécessairement oblique, elle peut être verticale ou autre. On définit une zone d'observation 4 qui est une surface autour de l'objet d'intérêt 3 et comprise dans le champ de vision de 20 la caméra 2. La présente invention a pour but de déterminer un taux de transformation entre le nombre de personnes ayant traversé la zone d'observation et le nombre de personnes ayant acheté ensuite l'objet d'intérêt. Cet objet peut être un produit quelconque dans un magasin. Plusieurs objets d'intérêt 3a peuvent être disposés dans un rayon de 25 magasin à proximité dudit objet d'intérêt 3 placé en exposition. Sur la figure 1, on voit que la caméra de surveillance 2 est reliée à un serveur distant 5 destiné à mettre en oeuvre le procédé selon l'invention. Ce serveur distant 5 peut également être relié à un terminal 13 d'une caisse. D'une façon générale, lorsque le magasin comporte plusieurs caisses, le 30 serveur distant 5 est relié à tous les terminaux de caisse. Le terminal 13 est également relié à un lecteur infrarouge 14 destiné à lire le code de chaque produit acheté. Par ce moyen, le serveur distant 5 peut enregistrer et comptabiliser le nombre produits, notamment d'objets d'intérêt 3a achetés en une journée, une semaine, un mois ou un an par exemple. Selon un 35 autre mode de réalisation, on peut prévoir un serveur 5 relié à une ou plusieurs caméras 2, ce serveur 5 ayant un rôle d'unité de traitement t pour les images provenant des caméras 2. On peut alors prévoir plusieurs serveurs 5, chacun relié à une ou plusieurs caméras d'une part, et également relié à un serveur externe de gestion 5bis. Par ailleurs, plusieurs terminaux de caisses 13 peuvent être reliés à un autre serveur externe 13bis qui peut être spécifiquement dédié au magasin. Le serveur externe 5bis peut se connecter au serveur externe 13bis de façon à récupérer des données afin de déterminer des taux de transformation. Lorsqu'on souhaite réaliser une analyse comportementale sur un mois io par exemple. On place l'objet d'intérêt 3 dans la zone d'observation. On détecte les trajectoires de toute personne traversant la zone d'observation durant ce mois d'analyse. On peut éventuellement isoler de ces trajectoires, celles que l'on considère pertinentes : c'est-à-dire les trajectoires pour lesquels les personnes correspondantes se sont arrêtées par exemple une 15 dizaine de secondes devant l'objet d'intérêt 3. Par exemple, en un mois d'observation, 2500 personnes se sont arrêtées au moins dix secondes devant l'objet d'intérêt 3. Durant la même période, 125 produits 3a se sont vendus dans le magasin. On peut donc en déduire un taux de transformation égal à 2500/125, soit 20. On peut aussi définir plutôt le taux 20 de transformation comme étant le rapport 125/2500, soit 0.05. Dans ce cas le taux de transformation est de 5 pourcent. On va maintenant décrire le processus de suivi des individus. Le suivi des individus par la caméra permet d'établir des flux dans la zone d'observation. 25 La figure 2 illustre la problématique de l'appariement. On distingue en b21-b25 des positions d'individus pour une image à l'instant t+1, et en a21-a24 des positions d'individus pour une image à l'instant t. On a donc déterminé avec précision la position de chaque individu dans la zone d'observation pour chaque instant. Par exemple pour la position a22, il 30 s'agit de savoir si sa position suivante est la position b21 ou la position b23. Avec le procédé selon l'invention, on réalise une minimisation de la distance globale. En fait on normalise la matrice de correspondance. Sur la figure 3, on voit que l'appariement est optimal par minimisation des distances de façon globale : a22 est apparié à b21, et a24 est apparié à b23. 35 On va maintenant décrire un peu plus en détail le processus de détection des individus. Sur la figure 4, on voit un organigramme d'un procédé selon l'invention mis en oeuvre dans le serveur distant 5. Cet organigramme a pour entrées une image acquise par la caméra de surveillance 2 et la zone d'observation 4. A l'étape ai, on extrait le fond de l'image acquise de façon à obtenir un modèle de fond fixe. L'algorithme d'extraction du fond peut être de type itératif sur notamment plusieurs images ou un flux d'image vidéo. io L'approche la plus simple pour extraire le fond consiste à moyenner les valeurs des pixels au cours du temps. La contribution à la moyenne d'un objet en mouvement sera d'autant plus faible que son mouvement est rapide par rapport à la cadence des images. Plus précisément, l'extraction de fond selon l'invention conduit à une 15 modélisation de niveaux de gris, gradients et orientations de gradients moyens et de leurs écarts types. A titre d'exemple, sur la figure 5 on voit une image acquise par la caméra de surveillance 2. On distingue des individus à l'intérieur et hors de la zone d'observation 4 en pointillée. Sur la figure 6, on a réalisé l'étape ai de sorte que les éléments en mouvements 20 hors fond fixe sont symbolisés par des blobs 12. A l'étape a2, on réalise une détection de zone en comparant l'image et le fond, et en extrayant des zones de pixels, appelées � blobs �, n'appartenant pas au fond. A l'étape bi, on réalise une détection de têtes consistant à rechercher 25 des positions possibles de têtes d'individus par analyse de contours. Pour ce faire on reçoit en entrée les blobs précédemment détectés et ladite image acquise. Le résultat de l'étape bi est un ensemble d'hypothèses de position de têtes. Puis, à l'étape b2, on réalise un dénombrement des individus. Plus 30 précisément, on réalise une identification du nombre d'individus dans chaque blob et une estimation de leur position. Sur la figure 7, on voit un peu plus en détail un algorithme de détection de têtes de l'étape bi. En entrée, on distingue l'image acquise, les 35 blobs se détachant du fond, et un modèle de tête appelé � template �. La -10 - première étape consiste à appliquer un filtre de Canny pour extraire les contours des têtes d'individus sur l'image acquise de façon à obtenir une image de contours dans la zone des blobs. Ensuite, on réalise une transformée en distance. Il s'agit d'un calcul d'une carte de distance entre contours. C'est un calcul rapide et stable de la corrélation du � template � de tête avec les contours. La sortie de cette étape est une carte de distances. Ensuite, on réalise une convolution entre la carte de distances et le � template � de façon à obtenir une carte de convolution. La convolution io comprend un calcul de la corrélation du � template � de tête avec les contours. Enfin, on applique un algorithme de ligne de partage des eaux (� watershed � en langue anglaise) pour localiser et quantifier des maxima de corrélation, et déterminer des maxima de probabilité de présence de 15 têtes dans l'image acquise. On obtient en sortie une hypothèse sur la position des têtes. Sur la figure 8, on voit un exemple de modèle d'apparence d'un individu. Ce modèle est utilisé à l'étape b2. En particulier, la détection d'hypothèses de têtes de l'étape bi génère un grand nombre de fausses 20 détections. Pour aboutir à l'interprétation complète de la scène, on cherche à identifier les individus à partir d'un modèle d'apparence à trois ellipses de la figure 13 (silhouettes pleines), en optimisant les critères suivants : - minimiser le nombre de personnes, - maximiser le nombre de pixels différents du fond et 25 appartenant aux silhouettes, - minimiser le nombre de pixels différents du fond et n'appartenant pas aux silhouettes, - minimiser le nombre de pixels appartenant au fond et aux silhouettes, et/ou 30 - minimiser la distance des contours d'un modèle 2D/3D. La convergence vers la solution optimale s'effectue par des méthodes itératives de type descente de gradient. L'évaluation de cette solution à partir de la superposition des modèles d'apparence permet de gérer naturellement les cas d'occlusion. 2927442 - 11- Sur la figure 9, on voit un autre exemple de procédé de détection d'individus. Dans cette variante, on fait l'hypothèse d'une taille standard de personne. La zone d'observation est considérée dans ses coordonnées réelles 3D, une "fenêtre de recherche" 6, en fait un parallélépipède 8 de la 5 taille d'une personne, est déplacé dans cette zone suivant un maillage régulier. Pour chaque position de cette fenêtre de recherche, la projection 7 du parallélépipède dans l'image acquise est calculée, et la vraisemblance de présence d'une personne à cette position est évaluée en fonction du remplissage de ce parallélépipède par des pixels des blobs. Le remplissage io est représenté sur le graphe 9. Cette variante permet donc de gérer la perspective de la scène, ainsi que les défauts de la caméra tels que la distorsion, ce paramètre étant pris en compte dans le calcul de la projection des parallélépipèdes. Cette variante permet également de prendre en compte le masquage 15 par des objets : si l'on sait qu'une zone correspond à un objet derrière lequel des personnes seraient susceptibles d'être partiellement masquées, l'estimation de la probabilité de présence est corrigée en fonction du nombre de pixels de l'objet masquant situés dans la fenêtre de recherche. Sur l'image 11 de la figure 9, après lissage, des maxima locaux de 20 vraisemblance sont calculés par une méthode de ligne de partage des eaux sur la vraisemblance inversée; chaque bassin est considéré comme une détection possible. Afin d'éviter des détections multiples d'un même individu (l'une des détections correspondant à son tronc, et est interprétée comme une 25 personne plus éloignée), les détections obtenues à l'étape précédente sont classées et considérées par ordre croissant de distance les séparant de la caméra de surveillance. On valide les détections de la plus proche à la plus lointaine. Pour qu'une détection soit validée, on prévoit à ce que le 3o parallélépipède soit rempli par des pixels d'un blob au-delà d'un taux prédéterminé. Chaque détection validée efface du blob les pixels contenus dans son parallélépipède, afin d'éviter les détections multiples. Sur la figure 10, on voit un exemple d'organigramme décrivant cette variante. En entrée, on prend en compte la zone d'observation 3D, les 35 caractéristiques de la caméra de surveillance (champ, position, 2927442 -12 - distorsion,...) et un modèle d'individu, de façon à réaliser la projection perspective 7 de chaque surface élémentaire de la zone d'observation. On obtient un modèle 2D des modèles 3D d'individu. Le calcul 9 de correspondance des modèles (le remplissage) entre une image des blobs et 5 le modèle 2D des modèles permet d'établir une carte de densité de probabilité de présence de chaque individu. En 11, on recherche donc les maxima par ligne de partage des eaux (� watershed �) de façon à obtenir les individus détectés. Bien sûr, l'invention n'est pas limitée aux exemples qui viennent io d'être décrits et de nombreux aménagements peuvent être apportés à ces exemples sans sortir du cadre de l'invention. REVENDICATIONS 1. Procédé pour déterminer un taux de transformation local d'un objet d'intérêt présent dans une zone d'observation, dans lequel procédé on acquiert une pluralité d'images de la zone d'observation au moyen d'une caméra, caractérisé en ce que pour chaque image acquise on réalise les étapes suivantes : a) détection de zones de pixels se détachant du fond fixe de l'image, b) détection d'individus en déplacement dans la zone d'observation, io c) pour chaque individu détecté, détermination de la position occupée par cet individu, et d) pour l'ensemble des positions déterminées, on minimise la distance globale entre les positions ainsi déterminées dans l'image en cours et des positions déterminées lors d'une image précédente, 15 et en ce que pour un intervalle de temps donné, on détermine la trajectoire de chaque individu ayant traversé la zone d'observation, et on détermine un taux de transformation comme étant le rapport entre le nombre de trajectoires ainsi déterminées et une valeur d'intérêt qui est fonction dudit objet d'intérêt. 20 2. Procédé selon la revendication 1, caractérisé en ce que ladite valeur d'intérêt est le nombre d'individus dont la trajectoire présente un temps d'arrêt de durée prédéterminée devant ledit objet d'intérêt. 3. Procédé selon la revendication 1, dans lequel plusieurs objets d'intérêt sont 25 disponibles, ce procédé est caractérisé en ce que ladite valeur d'intérêt est le nombre d'objets d'intérêt ayant été détectés par un appareil de détection de sortie d'objets d'intérêt. 4. Procédé selon la revendication 2 et 3, caractérisé en ce qu'on détermine en 30 outre un taux de transformation comme étant le rapport entre le nombre d'individus dont la trajectoire présente un temps d'arrêt de durée prédéterminée devant l'objet d'intérêt et le nombre d'objets d'intérêt ayant été détectés par l'appareil de détection de sortie d'objets d'intérêt. 15- 14 - 5. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la minimisation de la distance globale consiste à réaliser les étapes suivantes : pour chaque position d'individu déterminée, on évalue la distance par rapport aux autres positions d'individus, définissant ainsi une matrice de distances que l'on convertit en matrice de similitudes, puis on normalise des vecteurs propres de cette matrice de similitudes. 6. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la distance entre deux positions est définie en prenant en compte 10 une distance géométrique 3D. 7. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la distance entre deux positions est définie en prenant en compte une distance géométrique 2D. 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la distance entre deux positions est définie en prenant en compte une distance colorimétrique. 20 9. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la distance entre deux positions est définie en prenant en compte une cohérence temporelle. 10. Procédé selon l'une quelconque des revendications précédentes, 25 caractérisé en ce que l'étape b) de détection des individus comprend les étapes suivantes : - détection de têtes des individus dans chaque zone de pixels, et - dénombrement des individus dans chaque zone de pixels. 30 11. Procédé selon la revendication 10, caractérisé en ce que l'étape de détection de têtes comprend les étapes suivantes : - application d'un filtre de Canny sur l'image par rapport aux zones de pixels de façon à générer une image de contours des zones de pixels, - application d'une transformée en distance de façon à élaborer une carte 35 de distances, 5 2927442 - 15 - - réalisation d'une convolution entre la carte de distance et un modèle de tête (template) de façon à obtenir une carte de convolution, et - application d'un algorithme de ligne de partage des eaux (watershed) de façon à détecter des têtes d'individus. 12. Procédé selon la revendication 10 ou 11, caractérisé en ce que le dénombrement consiste à identifier simultanément le nombre d'individus présents et leur position respective. io 13. Procédé selon la revendication 12, caractérisé en ce que l'identification des individus est effectuée à partir d'un modèle d'apparence à trois ellipses imitant un individu. 14. Procédé selon la revendication 13, caractérisé en ce que le modèle 15 d'apparence est optimisé selon les critères suivants : - minimisation du nombre d'individus, - maximisation du nombre de pixels différents du fond et appartenant aux silhouettes, - minimisation du nombre de pixels différents du fond et n'appartenant 20 pas aux ellipses, - minimisation du nombre de pixels appartenant au fond et aux ellipses, et - minimisation de la distance des contours d'un modèle 2D/3D. 25 15. Procédé selon la revendication 14, caractérisé en ce que l'optimisation est réalisée au moyen d'un algorithme itératif de descente de gradient. 16. Procédé selon la revendication 15, caractérisé en ce qu'on répartit la convergence de l'algorithme itératif sur plusieurs images successives. 30 17. Procédé selon la revendication 1, caractérisé en ce que l'étape b) de détection des individus comprend les étapes suivantes : à partir d'un modèle 3D de la zone d'observation, d'un modèle 3D d'individu, et des caractéristiques liées à l'appareil d'acquisition de ladite pluralité 35 d'images d'observation, on réalise une projection perspective dans l'image en 20- 16- cours de chaque surface élémentaire de la zone d'observation de façon à obtenir une projection 2D des modèles d'individu ; on calcule la correspondance des modèles 2D sur chaque zone de pixels de façon à obtenir une carte de densité de probabilité de présence ; et on réalise une recherche de maxima par l'algorithme de ligne de partage des eaux de façon à déterminer les individus détectés. 18. Procédé selon la revendication 17, caractérisé en ce que les caractéristiques de l'appareil d'acquisition comprennent le champ de vision, le io positionnement et la distorsion. 19. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que préalablement à l'étape a), on réalise une étape d'extraction d'un fond de l'image de façon à déterminer un modèle de fond 15 fixe. 20. Procédé selon la revendication 19, caractérisé en ce que l'extraction du fond fixe de l'image est obtenue en moyennant les valeurs des pixels de l'image au cours du temps. 21. Procédé selon la revendication 20, caractérisé en ce qu'on modélise la valeur de chaque pixel de fond par une distribution de probabilités. 22. Procédé selon la revendication 20 ou 21, caractérisé en ce qu'avant de 25 moyenner les valeurs des pixels de l'image, on réalise une détection instantanée de mouvement par soustraction d'images successives et on applique un seuil afin d'obtenir un masque correspondant aux pixels qui ne seront pas moyennés. 30