FR 2962834 A1 20120120 FR 1002950 A 20100713 Procédé de détection d'une cible dans des images stéréoscopiques par apprentissage et classification statistique à partir d'une loi de probabilité La présente invention concerne un procédé de détection et de reconnaissance de cibles dans des images stéréoscopiques. Elle peut notamment s'appliquer dans tout domaine utilisant un couple de cameras stéréoscopiques pour détecter des objets. II peut s'agir par exemple d'applications de surveillance, d'assistance ou encore de sécurité embarquées à bord d'un véhicule. La problématique générale est de rechercher dans des images la présence de cibles de types divers, qu'il s'agisse d'objets ou de personnes, les cibles présentant certaines caractéristiques conformes à un modèle. Par exemple, il peut s'agir d'un modèle paramétrique, comme un rapport entre largeur et hauteur qui doit valoir une valeur À donnée ou encore d'un modèle tridimensionnelle de type CAO. Une telle méthode de détection de cibles basée sur un modèle devient difficile à mettre en oeuvre en cas de grande variabilité d'apparence des cibles. Par exemple, l'apparence d'une personne peut varier grandement en fonction de sa posture ou de son habillement. La méthode peut même devenir impossible à mettre en oeuvre. Par exemple, l'exploitant d'un parking aura d'immenses difficultés à détecter des camions dès lors qu'il ne dispose pas des modèles CAO des différents types de camion. Dans ces cas où la modélisation des cibles s'avère difficile voire impossible, une solution connue consiste à mettre en oeuvre une étape d'apprentissage statistique OFF-LINE, c'est-à-dire préalablement à l'exploitation du système de détection, et une étape de classification ONLINE, c'est-à-dire simultanément à l'exploitation du système de détection. En fait, l'étape de classification fait partie intégrante du processus de détection: si l'on considère un système pour détecter des piétons, une détection a lieu dès lors qu'une cible a été classifiée "piéton". L'étape préalable d'apprentissage statistique consiste à apprendre à reconnaître des cibles grâce à un algorithme, qui extrait automatiquement les paramètres les plus pertinents des cibles, afin de les discriminer par rapport aux autres éléments pouvant être présents sur les images. II s'agit en fait de constituer des modèles statistiques de données extraites d'une collection d'images "typiques" de cibles. Ces modèles statistiques sont utilisés plus tard lors de l'étape simultanée de classification. L'étape simultanée de classification est effectuée en temps réel sur les images les plus récemment fournies par les caméras. II s'agit de comparer de nouvelles données extraites des images "réelles" aux modèles statistiques constitués durant l'étape d'apprentissage à partir d'images "typiques". Ainsi, des systèmes permettent d'ores et déjà de détecter et reconnaître des objets ou des personnes, fixes ou mobiles, grâce à des couples d'images fournis par des cameras calibrées formant une tête stéréoscopique, par exemple deux caméras disposées horizontalement. Un tel système calcule tout d'abord, pour chaque couple d'images, une carte de disparité représentant le décalage entre l'image de gauche et l'image de droite. Plus exactement, la disparité est la différence en position de pixel entre deux images pour un même point de la scène observé. Par triangulation, cet écart permet de calculer la coordonnée z des pixels de l'image et donc d'obtenir des informations de profondeur (3D) sur la scène observée. Parfois représentée par des niveaux de gris, une telle carte de disparité est aussi communément appelée image de disparité. Puis, un tel système modélise l'apparence des objets présents sur les images lors d'un processus d'apprentissage statistique. Ce processus se base sur un ensemble de descripteurs calculés dans l'image, comme les niveaux de gris, les données RGB (Red-Green-Blue), les dérivées successives du signal, des convolutions par un ensemble de filtres spécifiques ou encore des histogrammes. Enfin, les flux vidéo fournis par les deux caméras permettent de calculer une carte des positions 2D estimées des pixels pour chacune des deux caméras gauche et droite. Cette information est importante pour distinguer les objets en mouvement. Elle permet de mieux segmenter les objets, notamment lorsque la caméra est fixe ou encore lorsque le mouvement des objets est suffisamment différent de celui de la caméra, comme par exemple un piéton qui traverse la rue devant une voiture en marche portant les caméras. Par exemple, l'article "Improved Multi-Person Tracking with Active Occlusion Handling" (A. Ess, K. Schindler, B. Leibe, L. van Gool, ICRA Workshop on People Detection and Tracking, May 2009) décrit une telle méthode de détection, de reconnaissance et même de suivi d'objets par utilisation d'une tête stéréoscopique. Elle met en oeuvre plusieurs étapes pour intégrer les informations de luminance, de position et de profondeur décrites précédemment. Cette méthode est basée sur la détection préalable de zones d'intérêt et sur la représentation de ces zones par un dictionnaire de motifs élémentaires. L'étape d'apprentissage comporte une étape d'apprentissage d'un "Codebook of Local Appearance", qui est un dictionnaire de motif visuels élémentaires que l'on peut rencontrer sur les objets, et une étape d'apprentissage des "Implicit Shape Models", qui sont les positions relatives de ces motifs élémentaires sur les objets. Lors de l'étape de classification en vue de détecter des objets, un premier détecteur recherche, dans les images et à partir du dictionnaire, des zones d'intérêt susceptibles de contenir des objets, puis un deuxième détecteur recherche des objets dans les zones d'intérêt. Enfin, un mécanisme de vote permet la reconnaissance des objets. Un inconvénient majeur de cette méthode est qu'elle repose sur plusieurs étapes successives de pré-détection, d'abord pré-détection de zones d'intérêt, puis d'objets dans les zones d'intérêt et puis seulement reconnaissance des objets par classification. Il peut en résulter un nombre important de non-détections, car ces étapes successives de pré-détections reviennent d'une certaine manière à un mécanisme de "tout ou rien" : si une étape en amont donne un résultat négatif, les étapes en aval ne sont même pas mises en oeuvre, alors qu'elles auraient pu s'avérer fructueuses en "corrigeant" la non-détection de l'étape amont. Et si, pour tenter d'atténuer ce mécanisme de "tout ou rien", on augmente le nombre de zones d'intérêt, alors on assiste à une véritable explosion des temps de calcul. L'invention a notamment pour but de limiter les non-détections en évitant d'enchaîner de multiples phases de pré-détection durant la phase de classification. Pour cela, l'invention propose notamment une méthode d'apprentissage statistique des objets permettant de fusionner à bas niveau les informations de luminance, de position et de profondeur sur les objets. A cet effet, l'invention a pour objet un procédé de détection d'une cible présente dans au moins deux images d'une même scène capturées simultanément par des caméras distinctes. Le procédé comporte une étape préalable d'apprentissage de cibles dans des conditions de mise au point et, dans des conditions d'exploitation, une étape simultanée de classification d'objets présents sur les images, la cible étant réputée détectée dès lors qu'un objet est classifié comme étant l'une des cibles apprises pendant l'étape d'apprentissage. L'étape d'apprentissage des cibles inclut une étape de modélisation des données X correspondant à une zone d'intérêt dans les images par une loi de distribution P telle que P(X) = P(X2d,X3d,XT) = P(X2d)P(X3d)P(XT) où X2d sont les données de luminance dans la zone d'intérêt, X3d sont les données de profondeur dans la zone d'intérêt et XT sont les données de position dans la zone d'intérêt. Dans un mode de réalisation préférentiel, l'étape de modélisation peut comporter une étape de calcul, en chaque pixel de la zone d'intérêt dans une des deux images, des valeurs dx et dy des dérivées de la luminance de pixel suivant deux directions orthogonales, P(X2d) pouvant être une loi de Von Mises-Fisher sur la sphère unité de R3 d'indice de concentration k2d et de vecteur moyen ,u2d, P(X2d) décrivant la distribution des vecteurs unitaires normaux aux plans d'équation z=dx.x+dy.y correspondant à tous les pixels de la zone d'intérêt. Dans un mode de réalisation préférentiel, l'étape de modélisation peut également comporter une étape de calcul, en chaque pixel de la zone d'intérêt dans une des deux images, des valeurs mx et my des dérivées de la position de pixel suivant deux directions orthogonales, P(XT) pouvant être une loi de Von Mises-Fisher sur la sphère unité de i3 d'indice de concentration kT et de vecteur moyen,uT, P(XT) décrivant la distribution des vecteurs unitaires normaux aux plans d'équation z=mx.x+my.y correspondant à tous les pixels de la zone d'intérêt. Dans un mode de réalisation préférentiel, l'étape de modélisation peut également comporter une étape de calcul, en chaque pixel de la zone d'intérêt dans l'image de disparité correspondant aux deux images, des valeurs px et py des dérivées de la profondeur de pixel suivant deux directions orthogonales, P(X3d) pouvant être une loi de Von Mises-Fisher sur la sphère unité de 3 d'indice de concentration k3d et de vecteur moyen ,u3d, P(X3d) décrivant la distribution des vecteurs unitaires normaux aux plans d'équation z=px.x+py.y correspondant à tous les pixels de la zone d'intérêt. Avantageusement, un vecteur unitaire ayant pour coordonnées (a,/3) dans la sphère unité de J3, la loi de Von Mises-Fischer peut être a- périodique et/ou /3-périodique de période 2 où n entier. n Dans un autre mode de réalisation, P(X3d) peut être une loi normale d'écart-type 03d et de moyenne ,u3d, P(X3d) décrivant la distribution 15 de la profondeur de pixel dans la zone d'intérêt. Par exemple, la zone d'intérêt peut être choisie par utilisation d'un classifieur de type boosting. 20 La présente invention a encore pour principal avantage qu'elle ne fait pas augmenter le taux de fausses alarmes pour un même taux de détection. De plus, présentant une complexité raisonnable, elle peut être mise en oeuvre même sur les matériels les plus légers. 25 D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la description qui suit faite en regard de dessins annexés qui représentent : les figures la et 1 b, par une image 2D et une image 3D 30 respectivement, un exemple d'image 2D en niveaux de gris ainsi que son interprétation, selon l'invention, comme une surface 3D; - la figure 2, par un diagramme de séquence, un exemple de mise en oeuvre de l'invention. La présente invention propose de fusionner les données au plus bas niveau, afin de permettre un apprentissage conjoint des informations 2D de luminance X2d , des informations 3D de profondeur X3d et des informations temporelles de position XT. Pour cela, une image est créée, image dont chaque pixel contient l'ensemble des données X de base calculées à partir des informations provenant des trois sources. La fusion des données est réalisée par utilisation d'une loi P de probabilité jointe entre les données X2d, X3d et XT que l'on suppose indépendantes, c'est-à-dire telles que : P(X) = P(X2d , X3d ,XT ) = P(X2d )P(X3d )P(XT Ainsi, les données X issues d'une zone d'intérêt rectangulaire quelconque 15 de l'image peuvent être modélisées par cette loi P fusionnant les informations de luminance X2d, de profondeur X3d et de position XT. La description qui suit correspond à un exemple de mise en 20 oeuvre de l'invention à partir d'une image gauche et d'une image droite fournies par deux caméras disposées horizontalement. Dans ce mode de réalisation, les données de luminance, de profondeur et de mouvement des pixels d'une zone d'intérêt rectangulaire de l'image peuvent avantageusement être représentées par un descripteur peu couteux à 25 calculer, à savoir en vecteurs de caractéristiques X2dX3dXT prêts à être comparés. En effet, les données du signal de luminance et de mouvement des pixels, qui sont toutes deux des données vectorielles 2D, peuvent être l'une et l'autre modélisées par une loi de probabilité sur données sphériques 30 3D. D'une part, le signal de luminance peut avantageusement être représenté par ses dérivées par rapport à deux directions orthogonales du plan de l'image. A titre d'illustration, l'image 2D en niveaux de gris illustrée par la figure la, chaque niveau de gris traduisant une valeur de luminance 35 d'un pixel, peut être interprétée comme une surface, comme illustrée par la figure lb, dont les hauteurs dans une direction Z sont données par la luminance des pixels. Ainsi, en un point quelconque de la surface peut être construit le plan tangent à la surface, plan d'équation z=a.x+b.y où a=dx est la valeur de la dérivée de la luminance par rapport à une direction X au point considéré et b=dy est la valeur de la dérivée de la luminance par rapport à une direction Y au point considéré. De ce plan tangent peut être facilement déduit le vecteur unitaire normal à la surface au point considéré, qui est orthogonal au plan. Ainsi, l'information de luminance d'un pixel peut être décrite par ce premier vecteur unitaire. Par exemple, le calcul de la dérivée du signal de luminance peut se faire aisément par masque de convolution 2D, dans l'image droite par exemple. Il s'agit là d'un traitement classique en traitement d'images. D'autre part, de manière similaire, le déplacement d'un pixel entre deux images qui se suivent dans le temps peut avantageusment lui aussi être représenté par ses dérivées par rapport à deux directions orthogonales du plan de l'image, ainsi que permettre de construire un plan tangent en un point quelconque d'une pseudo-surface, plan d'équation z=a.x+b.y où a=mx est la valeur de la dérivée de la position par rapport à une direction X au point considéré et b=my est la valeur de la dérivée de la position par rapport à une direction Y au point considéré, ainsi que permettre de construire le vecteur unitaire normal à la pseudo-surface au point considéré. Ainsi, l'information de position d'un pixel peut être décrite par ce deuxième vecteur unitaire. Par exemple, le calcul du déplacement des pixels entre un instant t-1 et un instant peut se faire aisément par corrélation, dans l'image droite par exemple. II s'agit là d'un traitement classique en traitement d'images. Il faut noter que l'invention permet de modéliser, dans une même loi, à la fois la direction des vecteur 2D, mais aussi leurs magnitudes. Enfin, toujours de manière similaire, la profondeur d'un pixel peut avantageusement elle aussi être représentée par ses dérivées par rapport à deux directions orthogonales du plan de l'image et être décrite par un vecteur unitaire tangent à un plan d'équation z=a.x+b.y avec a=px est la valeur de la dérivée de la profondeur par rapport à une direction X au pixel considéré et b=py est la valeur de la dérivée de la profondeur par rapport à une direction Y au pixel considéré. II faut noter que plusieurs opérations préalables sur les images gauche et droite sont effectuées: rectification des images gauche et droite pour le calcul de la carte de disparité, calcul de la disparité entre l'image droite et gauche par corrélation et enfin calcul de la carte de profondeur de l'image droite. II s'agit là de traitements classiques en traitement d'images. Ainsi, les données de luminance, de profondeur et de position peuvent être transformées de manière à être chacune d'entre elles représentée par un point 3D sur la sphère unité, c'est-à-dire par trois couples d'angles (a,/3) pour les données provenant des trois sources. Or, la distribution de ces points sur la sphère unité peut être 10 efficacement modélisée par une loi P de Von Mises-Fisher d'indice de concentration k réel et de vecteur moyen p sur la sphère unité : P(X /k,,u) = k eckuTX> 4,r sinh(k) 15 où sinh désigne la fonction sinus hyperbolique et ,ur désigne le vecteur transposé de p. Si l'on pose X = (a, /3) et ,u = (a,1, /3,u) , la loi P est donc donnée par: P(X / k' 'u) = k 4,z sinh(k) e 20 Au total, la loi de probabilité P représentant les données d'une des sources peut donc être modélisée par trois valeurs scalaires: l'indice k réel et les deux angles ap et . Par ailleurs, P est a-périodique et fi-périodique de période 27r . 25 Cependant, une telle loi P est uni-modale et donc mal adaptée aux ensembles de points présentant plusieurs modes ou points de concentration. Pour surmonter cet inconvénient et plutôt que de mettre classiquement en oeuvre une mixture de lois très coûteuse en temps de calcul sur l'ensemble d'une image, la présente invention propose avantageusement d'utiliser 30 plusieurs lois de périodes différentes. Ceci permet de réduire drastiquement les temps de calcul. En effet, il est possible de définir des lois de Von Mises- k(cos Côp)xcos(/8)+sin(fp )sin(ap)xsin(f)sin(a)+sin(fl )cos(ap)xsin(f)cos(a)) Fisher a-périodiques ou fi-périodiques de période 27r où n entier. Par n exemple, la loi de Von Mises Fisher d'a-période 2 est donnée par: n P(X / k, ,u, n) = cst ek(cos(,Qf,)xcos(fl)+sin(f~, )sin(f)xsin(na,)sin(na)+sin(fl, )sin(8)xcos(na,)cos(na)) où cst désigne une constante réelle. En supposant les lois de périodes différentes indépendantes entre elles, une loi de distribution jointe utilisant N loi de périodes différentes, où N entier, pour décrire les pixels d'une zone d'intérêt rectangulaire de l'image, peut être la suivante: 10 N P(X /k'u) = nP(Xj /k;'uj,ni) où k et ,u désignent des vecteurs rassemblant l'ensemble des paramètres des N lois de distribution. 15 Il faut noter que, si les données X3d sont de mauvaise qualité, P(X3d) peut être une loi normale d'écart-type a3d et de moyenne U3d définie par: P(X3d a3d , P3d) = 27Lla3d 1 (3d-P3d)z 20 La présente invention propose également une méthode pour choisir les zones les plus pertinentes pour la reconnaissance de cibles, c'est-à-dire les zones d'intérêt rectangulaires devant être décrites par une loi jointe comme explicité précédemment. Pour cela, l'invention propose 25 avantageusement un classifieur de type "boosting", ainsi qu'une métrique pertinente entre de tels descripteurs. Les techniques de type "boosting" regroupent des méthodes d'apprentissage qui permettent de construire un classifieur dit fort à partir de plusieurs classifieurs dits faible. C'est un algorithme itératif qui choisit de 30 façon optimale à chaque boucle un nouveau classifieur faible. Le classifieur faible réalise une partition de l'espace d'entrée. La valeur renvoyé par le classifieur est fonction de la zone dans laquelle se trouve le point représenté par les données X2d, X3d et XT. Un score de classification indique quelle classe d'objet est la plus vraisemblable en fonction des données. A l'issue de l'étape d'apprentissage, un ensemble de descripteurs 5 de luminance, de profondeur, et de position est obtenu à partir d'une collection d'exemples d'images représentant des objets O dont les classes sont connues. A chaque exemple est attaché un poids w; qui varie au cours de l'algorithme. Les paramètres des lois de probabilité P(X2d /O, w1) , P(X3d / O, w;) , et P(XT / O, wt) des données sachant les classes sont 10 calculés et pondérées par w;. Les votes du classifieur faible sont calculés en fonction des différentes partitions de l'espace qu'il réalise et en fonction du poids des exemples. Durant l'étape de classification, une fois le vecteur de la zone de l'image calculé, la log-vraisemblance L(Z/Oj) est calculée pour les 15 différents objets O. Puis les probabilités P(O; /Z) sont calculées par la règle de Bayes. Avantageusement, la zone d'intérêt choisie peut être celle associée à la classe j pour laquelle la probabilité P(O; /Z) est la plus grande. Le classifieur renvoie le score calculé a l'apprentissage. Les scores de l'ensemble des classifieurs/descripteurs locaux sont sommés pour produire 20 le score final de la zone d'intérêt. La figure 2 illustre par un diagramme un exemple de séquence d'étapes permettant de mettre en oeuvre l'invention à partir de l'image 25 gauche et de l'image droite fournies par les deux caméras disposées horizontalement. Tout d'abord, les images gauche et droite peuvent être rectifiées lors d'une étape 1. Puis, la disparité entre l'image droite et gauche, la carte de profondeur de l'image droite, la dérivée dû signal de luminance et le déplacement des pixels peuvent être calculés lors d'une étape 2. Puis, lors 30 d'une étape 3, à chaque pixel i de l'image droite peut être associé un vecteur contenant la profondeur du pixel i (z;), la dérivée du signal de luminance au pixel i suivant les lignes et les colonnes (Dx; et Dy;) de l'image et le déplacement du pixel i entre les instants t-1 et t (mx; et my;). Enfin, lors d'une étape 4, des images intégrales peuvent être calculées pour permettre des sommations rapides. Les étapes 1, 2, 3 et 4 sont détaillées dans la suite de la présente demande. Comme explicité précédemment, pour modéliser une zone Z de l'image selon l'invention, l'ensemble des données X2d , X3d et XT de ladite zone Z est modélisé par une loi de probabilité. Puis, lors du processus de détection, les données issues d'une zone de l'image sont comparées aux lois de probabilité calculées pour chaque classe d'objet. La probabilité que l'ensemble des données X. des pixels de la zone Z soit issu de la loi de probabilité associée à un objet O est donnée par: 10 P(Z / O) = [J P(Xi / O) = [J (P(X2dl / O)P(X3dl / O)P(XTi / O)) Le logarithme de cette probabilité est donnée par: 15 L(Z / O) =1n(fJ P(Xi / O)) _ I ln(P(Xi / O) = ln(P(X2di / O)P(X3di / O)P(XTi / 0) L(Z / O) = ln(P(X2di / O)) + E ln(P(X3di / O)) + E 1n(P(XTi / O)) 1 i I 20 Avec les formules suivantes pour chaque loi de probabilité : 1 ('3d -PZ )2 O /, O 1 2 ~3d2 P(X3d l 63d >1"3d ) = e 22r ja k° P X / k° o _ k2d e(e/I? tx2d ) ( 2d 2d ~2d) 47r sinh(k d ) k° O O kT e(kr°dixr) P(XT l kT 'uT) = sinh(k°) 25 k° ln(P(X2di / O)) = 1 ln( k2d = 47r sinh(k2d ) o ) + k0 2dp2dr ox2dl k° = N x ln( k2d 0 1n(P(XTi / O)) = 1 1n( kT o) + E ko~orx,i sinh(kT) 0 =Nxln( kT o )+k°,u°tExTi 42r sinh(kT ) Au final, le calcul de L(Z/O1) nécessite le calcul des sommes des vecteurs X2di et XTi pour les sources 2D et de mouvement, des valeurs scalaires z; et z;2 pour la source 3D sur la zone en question, puis le calcul du produit scalaire entre les vecteurs sommes obtenus et un vecteur de paramètres pour chaque classe d'objet à considérer. Puis, les données d2d , dT et dm issues respectivement du calcul de la dérivée, de la carte de profondeur et de mouvement sont transformées afin d'être comparées facilement aux différents modèles d'objets. Au départ ces données sont sous la forme: d2d= dx Z2d Mx sinh(k2°)) + kzâf~2dr xzat ZT Elles sont transformées en X2d = d2d d2d XT = et x3d = 2d , donc: dT Z3d dT z dx .Nidx2 +dy2 +Z 2 2d dy Vdx2 +dy2 +Z 2 2d Z2d -\idx2 + dy2 + Z 2 2d COS(fl2d ) sin(/32d) sin(a2d ) Sin(,62d) Cos(a2d ) X2d = Mx .\/mx2+my2 + ZT2 My i mx2 + my2 + ZT2 ZT cos(/3T) sin(aT) sin(aT ) sin(aT) cos(aT) X3d = Vmx2 + my2 +ZT2 jdzx2 +dzy2 +Z3d2cos(f33d ) sin(f3d) sin(a3d ) sin(f3d ) cos(a3d ) dzx jdzx2 + dzy2 + Z3d2 My jdzx2 + dzy2 + Z3d2 Z3 d II faut noter que ces transformations sont très rapides a effectuer et qu'elles différent peu des opérations classiques en traitement d'image. En effet, dx +dÿ correspond au calcul de la magnitude du gradient au carré. Lors de la transformation, une constante Z2d est ajoutée à la magnitude au carré avant de calculer sa racine. Cette constante permet de paramétrer plus finement la transformation en fonction des valeurs de magnitude que l'on rencontre typiquement dans les images. Elle permet d'adapter la transformation à des valeurs de magnitudes trop fortes ou trop faibles. Ensuite, les images intégrales de ces données sont calculées pour permettre leurs sommations rapides. Le calcul de la somme d'un vecteur de taille n sur une zone de taille quelconque de l'image se fait ainsi à coût fixe, ne nécessitant que n°x°2 additions, n soustractions et 4°x°n accès aux données des images. Il faut noter que, dans le cas d'un calcul classique, le nombre d'opérations dépend de la surface de la zone et devient prohibitif dès 10°x°10 pixels. Durant la phase de détection, l'image est parcourue par une zone d'intérêt dont les positions et les tailles sont fonction du plan du sol. Une boite d'intérêt 3D est définie dans la scène. Cette boite 3D va être successivement déplacée sur un ensemble de positions candidates à contenir le ou les objets recherchés. Pour chaque position, la projection de la boite 3D est calculée et un vecteur décrivant le contenu de la zone 2D est calculé. Puis un test de reconnaissance est effectué sur ces données. Le calcul de la fonction de vraisemblance entre des descripteurs locaux et les modèles résultant de l'apprentissage s'effectue par l'intermédiaire de la sommation des vecteurs descriptifs sur la zone en question. Ces sommes sont calculées par la technique des images intégrales qui permettent de calculer efficacement des sommes de valeurs sur une aire 2D. Enfin, la zone est ajoutée à la liste des zones positives ou elle est rejetée. L'invention décrite précédemment a encore pour principal avantage de fournir un modèle à la fois plus fidèle et plus simple des données extraites des différentes images, qu'il s'agisse de la luminance, de la profondeur ou de la position des pixels. La loi de probabilité utilisée pour les vecteurs 2D, à savoir le signal image et le mouvement dans l'image, permet la modélisation d'une source par seulement trois paramètres, à savoir deux angles et un paramètre de concentration, tout en étant plus performante qu'une loi jointe classique linéaire ou circulaire. De plus, pour une détection rapide, elle fusionne les informations 3D, 2D et de mouvement au plus bas niveau dans les descripteurs locaux, les descripteurs les plus pertinents étant ensuite choisis par l'algorithme de boosting lors de la phase d'apprentissage. Par ailleurs, les classifieurs y sont organisés en cascade dont les étages sont de plus en plus complexes. Cette technique permet d'accélérer la recherche, car les tests les plus complexes ne sont appliqués qu'aux zones de l'image les plus ambigües. L'invention propose ainsi un processus très simple et très régulier, qui peut être efficacement implémenté tant sur des ordinateurs du marché que sur des architectures dédiées. Enfin, l'invention décrite précédemment est applicable à toutes sortes d'applications de reconnaissance de cibles, dès lors qu'elles sont basées sur l'utilisation de capteurs vidéos stéréoscopiques et l'utilisation d'une méthode de classification statistique comportant une phase préalable d'apprentissage suivie d'une phase simultanée de détection. Il peut s'agir par exemple d'applications de surveillance, comme la surveillance d'une zone à l'avant ou à l'arrière d'un véhicule pour éviter les collisions, la surveillance d'une zone d'un bâtiment, l'aide aux personnes handicapées ou encore la navigation de robots mobiles autonomes.25 REVENDICATIONS1. Procédé de détection d'une cible présente dans au moins deux images d'une même scène capturées simultanément par des caméras distinctes, le procédé comportant: - dans des conditions de mise au point, une étape préalable d'apprentissage de cibles, et; - dans des conditions d'exploitation, une étape simultanée de classification d'objets présents sur les images, la cible étant réputée détectée dès lors qu'un objet est classifié comme étant l'une des cibles apprises pendant l'étape d'apprentissage; le procédé étant caractérisé en ce que l'étape d'apprentissage des cibles inclut une étape de modélisation des données X correspondant à une zone d'intérêt dans les images par une loi de distribution P telle que P(X) = P(X2d,X3d,XT) = P(X2d)P(X3d)P(XT) où X2d sont les données de luminance dans la zone d'intérêt, X3d sont les données de profondeur dans la zone d'intérêt et XT sont les données de position dans la zone d'intérêt. 2. Procédé selon la revendication 1, caractérisé en ce que, l'étape de modélisation comportant une étape de calcul, en chaque pixel de la zone d'intérêt dans une des deux images, des valeurs dx et dy des dérivées de la luminance de pixel suivant deux directions orthogonales, P(X2d) est une loi de Von Mises-Fisher sur la sphère unité de t3 d'indice de concentration k2d et de vecteur moyen P2d, P(X2d) décrivant la distribution des vecteurs unitaires normaux aux plans d'équation z=dx.x+dy.y correspondant à tous les pixels de la zone d'intérêt. 3. Procédé selon la revendication 1, caractérisé en ce que, l'étape de modélisation comportant une étape de calcul, en chaque pixel de la zone d'intérêt dans une des deux images, des valeurs mx et my des dérivées de la position de pixel suivant deux directions orthogonales, P(XT) est une loi de Von Mises-Fisher sur la sphère unité de R3 d'indice de concentration kT et de vecteur moyen PT, P(XT) décrivant la distribution des vecteurs unitaires normaux aux plans d'équation z=mx.x+my.y correspondant à tous les pixels de la zone d'intérêt.17 4. Procédé selon la revendication 1, caractérisé en ce que, l'étape de modélisation comportant une étape de calcul, en chaque pixel de la zone d'intérêt dans l'image de disparité correspondant aux deux images, des valeurs px et py des dérivées de la profondeur de pixel suivant deux directions orthogonales, P(X3d) est une loi de Von Mises-Fisher sur la sphère unité de R3 d'indice de concentration k3d et de vecteur moyen )U3d, P(X3d) décrivant la distribution des vecteurs unitaires normaux aux plans d'équation z=px.x+py.y correspondant à tous les pixels de la zone 1 o d'intérêt. 5. Procédé selon l'une quelconque des revendications 2, 3 ou 4, caractérisé en ce que, un vecteur unitaire ayant pour coordonnées (a,/3) dans la sphère unité de la loi de Von Mises-Fischer est a-périodique et/ou $- 15 périodique de période 2z où n entier. n 6. Procédé selon la revendication 1, caractérisé en ce que P(X3d) est une loi normale d'écart-type 03d et de moyenne /u3d, P(X3d) décrivant la distribution de la profondeur de pixel dans la zone d'intérêt. 7. Procédé selon la revendication 1, caractérisé en ce que la zone d'intérêt est choisie par utilisation d'un classifieur de type boosting.