FR 3033648 A1 20160916 FR 1500507 A 20150313 1 La présente invention concerne le domaine de la surveillance de lieux comme par exemple les terminaux de transport (aéroports, gares, ports), les sites militaires, les sites industriels, les lieux publics... De tels lieux sont le plus souvent équipés d'un réseau de caméras reliées à un poste de surveillance dans lequel se trouvent un serveur d'enregistrement des images capturées par les caméras et des écrans de visualisation des images capturées par les caméras. Ceci permet à des opérateurs surveillant les écrans de pouvoir détecter des comportements suspects et d'utiliser les images enregistrées pour visualiser une scène qui s'est déroulée dans le passé. L'un des intérêts de disposer d'un réseau de ca- méras réparties dans un lieu est qu'il permet de suivre le cheminement des individus dans le lieu. Ceci suppose toutefois qu'un nombre suffisant d'opérateurs soient présents dans le poste de surveillance, et ce d'autant plus que la fréquentation du lieu est importante. Un but de l'invention est de fournir un moyen permettant de faciliter le suivi d'individus dans un lieu équipé par exemple de caméras. A cet effet, on prévoit, selon l'invention, un procédé de détermination de trajets d'individus dans un lieu ayant des zones équipées de moyens de détection d'individus et dont certaines sont inaccessibles directement depuis d'autres zones, les moyens de détection étant reliés à une unité informatique de traitement incorporant une topographie du lieu et un modèle temporel mettant en relation un temps de parcours et au moins un trajet pos- sible entre les zones, le procédé comprenant les étapes de : - mémoriser dans l'unité informatique des évène- ments initiaux associant chacun une détection des carac- téristiques représentatives d'un individu dans une pre- 3033648 2 mière zone et une information temporelle de détection, - répertorier dans l'unité informatique des évènements postérieurs associant chacun une détection de caractéristiques représentatives d'individu dans une autre 5 des zones du lieu et une information temporelle de détec- tion, - reconstituer dans l'unité informatique un trajet de chaque individu dans le lieu depuis la première zone en agglomérant des évènements postérieurs à chaque 10 évènement initial par comparaison des caractéristiques représentatives de chaque évènement postérieur avec des caractéristiques représentatives de chaque évènement initial, et vérification d'une cohérence des trajets des individus d'une zone à l'autre à partir de la topographie 15 du lieu, du modèle temporel et des informations tempo- relles de détection. La première zone est ici la première zone dans laquelle les caractéristiques représentatives dudit au moins un individu ont été détectées. Il en résulte que la 20 première zone n'est pas forcément la même pour tous les individus. L'information temporelle de détection permet de déterminer à quel moment (heure, minute, seconde) a été capturée l'image dans laquelle les caractéristiques représentatives ont été détectées. Lorsqu'un individu est 25 détecté dans une deuxième zone, on cherche à le recon- naître, c'est-à-dire à déterminer s'il fait partie des individus détectés au préalable dans la première zone. Cette reconnaissance est effectuée par comparaison des caractéristiques représentatives détectées dans la deu- 30 xième zone aux caractéristiques mémorisées. Lors de la reconnaissance, il est tenu compte de la topographie du lieu, c'est-à-dire des trajets possibles entre les zones, et du temps écoulé entre la détection des caractéristiques représentatives dans la première zone et la détec- 35 tion des caractéristiques représentatives dans la deu- 3033648 3 xième zone (si la deuxième zone n'est pas matériellement atteignable depuis la première zone en moins d'un temps prédéterminé, un individu détecté dans la première zone ne peut se trouver dans la deuxième zone qu'une fois ce 5 temps écoulé). La prise en compte de la topographie du lieu et du modèle temporel permet de renforcer la fiabilité de la reconnaissance. De préférence, lorsqu'au moins un évènement pos- térieur, dit résiduel, n'a pu être aggloméré pour former 10 un trajet reconstitué l'intégrant, le procédé comprend les étapes ultérieures de déterminer une probabilité que l'évènement postérieur résiduel suive le dernier évènement postérieur de chaque trajet reconstitué et d'agglomérer l'évènement postérieur résiduel à un des 15 trajets reconstitués en fonction de cette probabilité. Ce cas se présente par exemple lorsque la zone de l'évènement postérieur résiduel ne peut en théorie pas être atteinte sans que l'individu ait été détecté dans une zone intermédiaire située entre la zone de 20 l'évènement postérieur résiduel et la zone du dernier évènement postérieur de chaque trajet reconstitué. Il se peut en effet que la détection n'ait pas fonctionné lors du passage de l'individu dans la zone intermédiaire par exemple parce que le capteur était défaillant à ce mo- 25 ment-là, que l'individu a été masqué par un autre ou que l'individu ne soit pas rentré dans le champ d'une caméra mobile. Le procédé de l'invention permet de tenir compte de ces situations. Avantageusement, lors des étapes ultérieures, à 30 partir de la topographie sont recherchés des chemins pou- vant amener un individu à la zone de l'évènement postérieur résiduel en partant de la zone du dernier évènement postérieur de chaque trajet reconstitué et une densité de probabilités est calculée pour l'ensemble de ces chemins. 35 L'invention a également pour objet un dispositif 3033648 4 pour la mise en oeuvre de ce procédé. D'autres caractéristiques et avantages de l'invention ressortiront à la lecture de la description qui suit de modes de réalisation particuliers non limita- 5 tifs de l'invention. Il sera fait référence aux dessins annexés, parmi lesquels : - la figure 1 est une vue schématique d'un lieu équipé d'un dispositif de suivi pour la mise en 10 oeuvre de l'invention ; - la figure 2 est un vue d'une topographie représentative du lieu et exploitable pour la mise en oeuvre de l'invention. En référence aux figures, l'invention est ici dé- 15 crite en application à la surveillance d'un lieu L, ici sous la forme d'un hangar, ayant des zones ZO, Zl, Z2, Z3, Z4, formant des halls dans le lieu L. La zone ZO est un hall d'entrée du lieu L et la zone Z4 est le hall de sortie du lieu L. La zone ZO communique avec les zones 20 Zl, Z2, Z3 d'une part et avec l'extérieur du lieu L d'autre part. La zone Zl communique avec les zones ZO et Z4. La zone Z2 communique avec les zones ZO et Z4. La zone Z3 communique avec les zones ZO et Z4. La zone Z4 communique avec la zone Z1, la zone Z2 et la zone Z3. On 25 note que : la zone Z4 est inaccessible directement depuis la zone ZO et inversement ; la zone Z3 est inaccessible directement depuis la zone Zl et inversement. Chacune des zones ZO, Zl, Z2, Z3, Z4 est équipée d'au moins une caméra CO, Cl, C2, C3, C4 agencée pour 30 capturer des images des individus évoluant dans les zones ZO, Zl, Z2, Z3, Z4, images qui soient suffisamment définies pour qu'on puisse détecter des caractéristiques représentatives des individus se trouvant sur les images. Ces caractéristiques représentatives comprennent par 35 exemple les vêtements, la coiffure, toutes caractéris- 3033648 5 tiques biométriques dont notamment les lignes du visage. Dans la zone Zl, la caméra Cl est de préférence disposée au voisinage de l'entrée, d'un comptoir d'accueil ou de contrôle d'accès où tout individu pénétrant dans le lieu 5 L est forcé de passer et éventuellement de présenter au moins un document attestant son identité ou une autorisation d'accès : on est donc assuré d'avoir une image de tout individu étant régulièrement entré dans le lieu L. De la même manière, la caméra 07 est de préférence dispo- 10 sée près de la sortie pour capturer des images de tout individu quittant régulièrement le lieu L. Le procédé de l'invention est mis en oeuvre au moyen d'un dispositif de suivi, généralement désigné en 1, comprenant une unité informatique de traitement 2 qui 15 est reliée aux caméras CO, Cl, 02, 03, 04 et qui est agencée pour traiter des données transmises par les caméras CO, Cl, 02, 03, 04. L'unité de traitement 2 a une mémoire 3 contenant une topographie du lieu L (représentée schématiquement sur la figure 2) et un modèle temporel 20 mettant en relation des temps de parcours moyens et des trajets entre les différentes zones ZO, Zl, Z2, Z3, Z4. Le modèle temporel est basé par exemple sur une courbe gaussienne centrée sur un temps de parcours moyen sur chaque trajet. 25 L'unité de traitement exécute un programme infor- matique de suivi d'individus. Le programme analyse les images capturées par les caméras CO, Cl, 02, 03, 04, images qui lui sont transmises au fur et à mesure de leur capture. 30 Pour chaque zone, le programme est agencé pour détecter dans les images qui lui sont transmises des caractéristiques représentatives de chaque individu qui s'y trouve et pour mémoriser ces caractéristiques en association avec une information temporelle de détection, un 35 identifiant d'individu et un identifiant de zone. 3033648 6 L'information temporelle de détection permet de déterminer à quel moment (heure, minute, seconde) a été capturée l'image dans laquelle les caractéristiques représentatives ont été détectées. Avant la mémorisation des carac- 5 téristiques représentatives, il est procédé à une étape de reconnaissance consistant à comparer ces caractéristiques représentatives avec des caractéristiques représentatives préalablement mémorisées pour déterminer si l'individu détecté dans ladite zone a été détecté dans 10 d'autres zones. La comparaison est par exemple réalisée par la mise en oeuvre des techniques d'identification biométrique également appelées techniques de « matching ». Une étape de confirmation de la validation est également effectuée à partir d'une topographie du lieu par vérifi- 15 cation d'une cohérence d'un trajet de l'individu d'une zone à l'autre et du modèle temporel par comparaison des informations temporelles de détection des caractéris- tiques représentatives dans les zones. Dans l'affirmative, les caractéristiques représentatives nou- 20 vellement détectées sont enregistrées en association avec l'identifiant préexistant ; dans la négative, les caractéristiques représentatives nouvellement détectées sont enregistrées en association avec un nouvel identifiant choisi arbitrairement. 25 En théorie, les identifiants sont créés pour les individus nouvellement détectés dans les images capturées dans la zone ZO (il n'y a pas d'entrée dans les autres zones) et les identifiants et les données associées sont supprimées lorsqu'on détecte la sortie des individus par 30 la zone Z4. Il est donc possible de suivre automatiquement un individu dans ses déplacements dans le lieu L. Un exemple particulier va maintenant être décrit pour illustrer le principe du procédé (les individus men- 35 tionnés dans l'exemple qui suit ne sont pas représentés 3033648 7 sur la figure). A un instant donné, la mémoire contient les données suivantes : - identifiant II, avec ses caractéristiques re- 5 présentatives et les informations temporelles de détection en ZO et Z1 ; - identifiant 12, avec ses caractéristiques représentatives et les informations temporelles de détection en ZO, Z2 ; 10 - identifiant 13, avec ses caractéristiques re- présentatives et les informations temporelles de détection en ZO, Z3 ; - identifiant 14, avec ses caractéristiques représentatives et l'information temporelle de 15 détection en ZO. Un individu pénètre dans la zone ZO et apparaît sur les images capturées par la caméra CO. Ses caractéristiques représentatives sont détectées et comparées à celles mémorisées. Comme aucune ne correspond, les carac- 20 téristiques représentatives et l'information temporelle de détection sont mémorisées en association avec un nouvel identifiant 15. On notera que, pour aller plus vite lors de la comparaison, on aurait pu choisir d'utiliser pour la comparaison les caractéristiques des individus 25 dernièrement détectés dans la zone ZO ou une zone commu- niquant avec la zone ZO (dans le présent exemple simplifié tous les individus sont concernés), les autres individus ne pouvant atteindre la zone ZO directement depuis la zone dans laquelle ils ont été dernièrement détectés. 30 Imaginons maintenant que les caractéristiques re- présentatives de l'individu il sont proches de celles de l'individu il et de celles de l'individu 13. Sur des images de la zone Z4 sont détectées des caractéristiques représentatives qui sont comparées à celles mémorisées. 35 La comparaison fait ressortir une liste de trois résul- 3033648 8 tats possibles, à savoir Il, 13, 14, avec un taux de correspondance pour chacun. Il est procédé alors à une confirmation à partir de la topographie du lieu L par vérification d'une cohérence d'un trajet de l'individu d'une 5 zone à l'autre et du modèle temporel par comparaison des informations temporelles de détection des caractéristiques représentatives dans les zones. En l'occurrence : - Il se trouvait dernièrement dans la zone Z1, - 13 se trouvait dernièrement dans la zone Z3 10 d'où il peut accéder directement à la zone Z4, - 14 se trouvait dernièrement dans la zone ZO d'où il ne peut accéder directement à la zone Z4. L'étude de la topographie du lieu L permet de di- 15 minuer la probabilité que les caractéristiques représen- tatives détectées dans la zone Z4 appartiennent à 14, voir d'écarter carrément cette possibilité. Les caractéristiques représentatives ont été dé- tectées dans la zone Z4 à l'instant t. Le temps moyen 20 pour accéder à la zone Z4 est de 60 secondes depuis la zone Z1 et de 30 secondes depuis la zone Z3 d'après le modèle temporel. Or, les caractéristiques représentatives de I1 ont été détectées dernièrement dans la zone Z1 à t-30s et les caractéristiques représentatives de 13 ont été 25 détectées dernièrement dans la zone Z3 à t-45s. Ainsi, la probabilité que il se trouve à l'instant t en zone Z4 est plus faible que la probabilité que 13 se trouve à l'instant t en zone Z4. On notera que le procédé de l'invention peut être 30 mis en oeuvre en temps réel, au fur et à mesure de la cap- ture des images et de leur réception par l'unité de traitement, ou en différé sur la base d'images mémorisées précédemment. Dans le détail, le procédé de l'invention est mis 35 en oeuvre pour la détermination de trajets d'individus 3033648 9 dans le lieu L et comprend les étapes de : - mémoriser dans l'unité informatique des évènements initiaux associant chacun une détection des caractéristiques représentatives d'un individu dans une pre- 5 mière zone et une information temporelle de détection, - répertorier dans l'unité informatique des évènements postérieurs associant chacun une détection de caractéristiques représentatives d'individu dans une autre des zones du lieu et une information temporelle de détec- 10 tion, - reconstituer dans l'unité informatique un trajet de chaque individu dans le lieu depuis la première zone en agglomérant des évènements postérieurs à chaque évènement initial par comparaison des caractéristiques 15 représentatives de chaque évènement postérieur avec des caractéristiques représentatives de chaque évènement initial, et vérification d'une cohérence des trajets des individus d'une zone à l'autre à partir de la topographie du lieu, du modèle temporel et des informations tempo- 20 relies de détection. Ainsi, l'unité informatique de traitement mémo- rise : - pour chaque évènement initial, les caractéristiques biométriques d'un individu détectées dans la zone 25 ZO, un identifiant affecté à cet évènement (ou à cet in- dividu) et une information temporelle de détection (ici la date et l'heure de la détection) ; - pour chaque évènement postérieur, les caractéristiques biométriques d'un individu détectées dans une 30 des zones Zl à Z4 du lieu L, un identifiant affecté à cet évènement et une information temporelle de détection (ici la date et l'heure de la détection). La reconstitution des trajets des individus dans le lieu L depuis la zone ZO est effectuée par un algo- 35 rithme agglomérant des évènements postérieurs à chaque 3033648 10 évènement initial par : - comparaison des caractéristiques représentatives de chaque évènement postérieur avec des caractéristiques représentatives de chaque évènement initial, et 5 - vérification d'une cohérence des trajets des individus d'une zone à l'autre à partir de la topographie du lieu, du modèle temporel et des informations temporelles de détection. Pour effectuer la reconstitution, il est calculé 10 trois probabilités que deux évènements El, E2 soient liés à une même personne : - une probabilité biométrique pbio(Ei,E2) en fonction de la similarité des caractéristiques biométriques, - une probabilité topographique ntopo(ElP E2) déteru 15 minée sur la base de la topographie, - une probabilité temporelle ptemps('1,E2) représentative d'une distribution de probabilité de temps de passage entre les deux zones. Ces probabilités sont utilisées ensuite pour cal- 20 culer une probabilité globale p qui vaut mi \Y P = (Pbio(ElP E2))a (Ptemps(ED E2))e Vtopo(Ei,E2)) dans laquelle les exposants a, p, y permettent de pondérer les probabilités les unes par rapport aux autres. Le procédé est mis en oeuvre en deux temps. 25 Dans une première phase, l'algorithme détermine une matrice permettant d'associer les différents évènements entre eux pour créer des trajets et identifier les trajets les plus probables. Dans cette matrice, chaque ligne représente pour 30 un évènement toutes les associations possibles prenant cet évènement comme évènement antérieur et chacun des évènements comme évènement postérieur. Dans chaque cellule de la ligne, on indiquera 0 lorsque les deux évènements associés sont les mêmes ou la 35 probabilité globale p que l'évènement puisse être associé 3033648 11 à un autre évènement. L'algorithme choisit dans chaque ligne l'association présentant la plus forte probabilité, cette probabilité devant néanmoins être supérieure à un seuil 5 prédéterminé pour être considérée comme vraie par l'algorithme. La ligne et la colonne à l'intersection desquelles se trouve la probabilité sélectionnée sont alors remplies de 0 pour signifier qu'aucune autre association ayant le même évènement de départ et le même évè- 10 nement d'arrivée ne peut être retenue. Ceci est répétée jusqu'à ce qu'il ne soit plus possible de sélectionner des cellules. En l'absence de problème de détection, il est possible d'agglomérer les évènements associés pour former 15 des trajets reconstitués depuis un évènement initial pour tous les individus dont les caractéristiques représentatives ont été détectées dans la zone ZO. En revanche, en cas de problème de détection, il est possible que des évènements postérieurs n'aient pas 20 pu être agglomérés à un évènement initial. On a alors : des trajets reconstitués complets entre la zone ZO et la zone Z4 ; des trajets reconstitués incomplets depuis la zone ZO mais avec des évènements manquants ; et des évènements postérieurs ou associations d'évènements posté- 25 rieurs non rattaché à un évènement initial. On parle alors d'évènement postérieur résiduel pour désigner un évènement seul qui n'a pu être aggloméré dans un trajet reconstitué depuis un évènement initial et de trajet résiduel pour désigner une association d'évènements posté- 30 rieurs qui n'a pu être aggloméré dans un trajet reconsti- tué depuis un évènement initial. Dans une deuxième phase, l'algorithme exécute alors les étapes ultérieures de déterminer la probabilité p que l'évènement postérieur résiduel, noté Er, suive le 35 dernier évènement postérieur de chaque trajet reconsti- 3033648 12 tué, noté Et, et d'agglomérer l'évènement postérieur résiduel Er à un des trajets reconstitués en fonction de cette probabilité. On s'intéresse ici aux trajets reconstitués qui sont incomplets en considérant qu'aucun évène- 5 ment ne manque dans les trajets reconstitués complets. La probabilité que l'évènement postérieur rési- duel suive le dernier évènement postérieur Et d'un trajet reconstitué identifié lors de la première phase pour former un nouveau trajet reconstitué c est le produit de : 10 - la probabilité phi° que les caractéristiques re- présentatives de l'évènement postérieur Et appartiennent au même individu que les caractéristiques représentatives de l'évènement postérieur résiduel Er, - une probabilité ptOPOqu'un individu passe de la 15 zone du dernier évènement postérieur Et à la zone de l'évènement postérieur résiduel Er sans avoir été détecté, - une probabilité e qu'un individu passe de la zone du dernier évènement postérieur du trajet reconsti- 20 tué à la zone de l'évènement postérieur résiduel dans le temps séparant l'information temporelle du dernier évènement postérieur Et et l'information temporelle de l'évènement postérieur résiduel. La probabilité ee° est le produit de : 25 - une probabilité d'erreur de détection pcFTA u représente la probabilité qu'un individu suive le trajet c sans être vu sur une zone de ce trajet connaissant le taux d'erreur de détection du système (on désigne usuellement ce taux par FTA, il est fixé par les concepteurs 30 du système comme un compromis entre le risque de ne pas tenir compte d'une détection et le risque de tenir compte d'une fausse détection), - une probabilité pFlux c qu'un individu passe de la zone du dernier évènement postérieur Et à la zone de 35 l'évènement postérieur résiduel E, compte tenu d'un his- 3033648 13 torique des trajets reconstitués. La probabilité pcF-TA est calculée comme suit pcFTA = (PFTA)N avec Ali le nombre de zones sans évènement sur le chemin c. La probabilité euxse calcule en fonction 5 des probabilités de passage d'un lieu de capture à un autre P1-c2 (par exemple pour la probabilité de passage devant le capteur Cl puis devant le capteur 02) qui s'obtient par observation du flux de personnes évoluant dans le lieu par le passé. On a pcFlux = (Ci dé_ u 10 signant le capteur d'indice i). Finalement, on obtient pctopo = (pcFTA yeux)8 pour pondérer l'importance relative des probabilités en4 Flux topo et pc dans le calcul de la probabilité p°"°. Pour un chemin c comprenant des détections suc- 15 cessives par les capteurs CO, Ci et CN (i variant de 1 à N-1), il est possible de calculer la probabilité qu'un individu ait parcouru celui-ci en un temps t grâce au modèle temporel. Celui-ci comprend des lois de temps entre chaque couple de zones qui communiquent (et donc chaque 20 couple de capteurs susceptibles de détecter successive- ment un même individu sans passer par un autre capteur). Nous pouvons les combiner pour créer une loi par chemin en effet : n rcT=pii_>6.2*...*piN_i_,cN avec * symbolisant un produit de convolution. 25 L'algorithme va alors calculer la densité de pro- babilité pg(Et,Er) que les évènements Et, Er soient liés compte tenu de la topographie et du modèle temporel. La densité de probabilité s'écrit : Pg (Et, Er) = EcEch(Et,Er) K. (Et, Er) lep() avec Ch(Et,Er) l' ensemble des 30 chemins c = Ct -> C2 ... -> Cr entre Et et Er . Il sera ici cal- culé en réalité : pe (Et, Er) = pg (Et, Er) cECh(Et,Er) OÙ pf (Et, Er) = (Et, Er)Y3 (pFTA)y (pFlux)8 Comme il peut y avoir beaucoup de chemins entre Et avec y et ô sont des exposants choisis 3033648 14 et Er (voir une infinité dans le cas d'une topographie avec boucle), l'algorithme réalise une approximation pour accélérer les calculs. L'approximation ici retenue consiste à négliger les trajectoires de probabilités trop 5 faibles. Un seuil e a été sélectionné tel qu'est négligé r Fluxv5 tout chemin pour lequel q)TA\y t. j q)Ft j ---->CP+1 ).1 Co-->---->Cn+i n - Si (317A )y (pFlUX p )S> Co--*--- Cn+1 P C0- n+1 > e alors poser 25 n: =n+1 et Cn := C÷i puis répéter l' étape 3. - Sinon passer 2. Si Cli:+1 = Cr - Calculer 30 Pg ccnp+i (Et, Er) = (p7c: (Et, Er))fl (7,7A)y (ppux)8 Si P:0p (Et, Er) > E , alors pg (Et, Er): = pg (Et, Er) + It+...,(Et, Er) cnp+i - Si (p7c. (Et, Er)) f3 (pcFna)y(ppux)8 > e alors poser n: =n+1 et Cn := CnP+1 et Répéter l'étape 3 . a l'n+1 3033648 15 p+1 - Sinon passer à Cr1+1 4. Retournerpg(Et,Er). Fin de la routine lorsqu'ont été parcourus tous les chemins permettant d'atteindre Er depuis Et. 5 L'algorithme peut dès lors calculer la probabili- té 'g qu'une trajectoire identifiée soit liée à un résidu à partir, d'une part, de la densité de probabilités pe(Et,Er) et, d'autre part, la probabilité Nio. Ainsi, à chaque couple Et, Er, l'algorithme associe une densité de 10 probabilités K. = (Pbio)crPg (Et, Er) - L'algorithme exploite alors les probabilités p.f. de la manière suivante : 1. Créer un tableau dont chaque ligne correspond à un trajet reconstitué et chaque colonne à un résidu et 15 dans lequel chaque cellule contient la probabilité définie précédemment. 2. Sélectionner la cellule de probabilité maximum supérieure à un seuil positif prédéterminé et associer le trajet reconstitué et le résidu de la cellule sé- 20 lectionnée pour former un nouveau trajet reconsti- tué. 3. Si aucune cellule n'est sélectionnable le processus est interrompu. Sinon, l'algorithme revient à l'étape 1 ci-dessus. 25 A titre d'exemple numérique, on applique la deu- xième phase du procédé de l'invention à la topographie représentée en figure 2. On choisit : OE = P = y = 5 1, 30 E = 10-14, PFTA - 5%. L'individu a mis 16 secondes pour se rendre de ZO à Z4. Par hypothèse, 3033648 16 - le modèle temporel prévoit que : Pio ci C4(Et) Er) = Pio C2 C4(Et, Er) = Pi0 C3 C4(Ete Er) = 1,5 .10-6. - les lois de flux prévoient que : 5 topo topo topo = 1/3, PCO Cl - - PCO C2 - PCO C3 = 1. topo topo- PC1 C4 - - PC2 C4 PC3 C4 Au début de l'algorithme, la densité de probabilités pg(Et,Er) = O. On commence le calcul de la densité de probabilités en ZO (capteur CO) et p9(CO) = 1 soit supérieure à E. 10 Ce chemin étant d'intérêt, la routine d'exploration est poursuivie vers la zone Z1 pour calculer la probabilité du chemin ZO vers Zl. Cette probabilité pour le chemin ZO vers Zl est égale à (pcFTAyyr FluxJ8 1 J =-3>e- 15 Comme Z1 n'est pas Z4 et que le chemin est d'intérêt, la routine d'exploration est poursuivie vers la zone Z4 pour calculer la probabilité du chemin ZO vers Zl vers Z4 en suivant la topographie. On arrive finalement à 20 pg (17t,4) = 1/3x1x1,15.10-6X0,05 = 1,9.10-8. CO->C1->C4 Cette valeur étant supérieure à E, on calcule alors : pg (Et, Er) = pg (Et, Er) + Pg (Et, Er) = 0 + 1,9.10_8, CO->C1-$C4 soit 1,9.10-8. Comme il n'est pas possible de reboucler sur Z4, 25 l'algorithme remonte en Z1 ou le seul lien est vers Z4. L'algorithme remonte alors en ZO et recherche un nouveau chemin, qui passe par Z2 puis Z4. La probabilité pour le chemin ZO vers Z2 est égale à (KTA)y(pcF/uxr = _31 > e On cal- cule alors : 30 Pg (Et, Er) = 1/3x1x1,15.10-6x0, 05 = 1,9.10-8 CO->C2-)C4 qui est supérieure au seuil et on calcule pg(Et,Er)= pg (Et,Er) + Pg (Et,Er)+ Pg (Et,Er)=3, 8.10-8. CO-4C1->C4 CO->C2->C4 Comme précédemment, il n'est pas possible de reboucler sur Z4, l'algorithme remonte en Z2 ou le seul 3033648 17 lien est vers Z4. L'algorithme remonte alors en ZO et re- cherche un nouveau chemin, qui passe par Z3 puis Z4. La probabilité pour le chemin ZO vers Z3 est égale à (pcFrit)Y(pcnux)5 = > e 3 5 On calcule alors : /39 (Et,Er) = 1/3x1x1, 15 . 10-EX0, 05 = 1,9.108 CO->C3->C4 qui est supérieure au seuil et on calcule Pg (Et, Er) = Pg (Et,Er) + Pg (Et, Er) + Pg (Et, Er) CO->C1->C4 CO->C2->C4 --FPg (Et,Er)= 5,7.10-8. CO->C3->C4 Comme tous les chemins ont été explorés, la rou- 10 tine est arrêtée et la densité de probabilités pe(Et,Er) vaut 5,7.10-8. C'est ce qui sera reporté dans le tableau. Bien entendu, l'invention n'est pas limitée aux modes de mise en oeuvre décrits mais englobe toute variante entrant dans le champ de l'invention telle que dé- 15 finie par les revendications. En particulier, les moyens de détection peuvent comprendre d'autres types de détecteurs comme un capteur d'empreintes digitales. L'identifiant attribué ici de manière arbitraire 20 peut, selon une variante, correspondre à un identifiant figurant sur un document présenté par l'individu à l'entrée du lieu (document d'identité, carte d'embarquement ou autre autorisation d'accès...). Le modèle temporel peut être basé sur un histo- 25 gramme élaboré à partir des données acquises précédemment ou toute autre loi tirée des données historiques antérieures. Bien qu'ici les zones soient séparées physique- ment les unes des autres, les zones peuvent ne pas être 30 cloisonnées. Une zone correspond alors par exemple à la zone de couverture d'une caméra. Les moyens de détection peuvent comprendre plu- sieurs capteurs biométriques, les capteurs biométriques étant de types différents pour capturer des caractéris- 3033648 18 tiques biométriques de types différents avec par exemple : un détecteur d'empreintes digitales et un détecteur de visage dans une zone d'entrée, un détecteur d'empreintes digitales dans une des zones suivantes et un 5 détecteur de visage dans une autre des zones suivantes. REVENDICATIONS1. Procédé de détermination de trajets d'individus dans un lieu ayant des zones équipées de moyens de détection d'individus et dont certaines sont inaccessibles directement depuis d'autres zones, les moyens de détection étant reliés à une unité informatique de traitement incorporant une topographie du lieu et un modèle temporel mettant en relation un temps de parcours et au moins un trajet entre les zones, le procédé compre- nant les étapes de : - mémoriser dans l'unité informatique des évènements initiaux associant chacun une détection des caractéristiques représentatives d'un individu dans une pre- mière zone et une information temporelle de détection, - répertorier dans l'unité informatique des évènements postérieurs associant chacun une détection de caractéristiques représentatives d'individu dans une autre des zones du lieu et une information temporelle de détec- tion, - reconstituer dans l'unité informatique un trajet de chaque individu dans le lieu depuis la première zone en agglomérant des évènements postérieurs à chaque évènement initial par comparaison des caractéristiques représentatives de chaque évènement postérieur avec des caractéristiques représentatives de chaque évènement initial, et vérification d'une cohérence des trajets des individus d'une zone à l'autre à partir de la topographie du lieu, du modèle temporel et des informations tempo- relies de détection. 2. Procédé selon la revendication 1, dans lequel, lorsqu'au moins un évènement postérieur, dit résiduel, n'a pu être aggloméré dans un trajet, le procédé comprend les étapes ultérieures de déterminer une probabilité que l'évènement postérieur résiduel suive le dernier évène- 3033648 20 ment postérieur d'un trajet reconstitué considéré comme incomplet et d'agglomérer l'évènement postérieur résiduel à un des trajets reconstitués en fonction de cette probabilité. 5 3. Procédé selon la revendication 2, dans lequel la probabilité que l'évènement postérieur résiduel suive le dernier évènement postérieur d'un trajet reconstitué considéré comme incomplet pour former un nouveau trajet reconstitué est le produit de : 10 - une probabilité que les caractéristiques repré- sentative de l'évènement postérieur du trajet reconstitué appartiennent au même individu que les caractéristiques représentatives de l'évènement postérieur résiduel à partir d'un modèle des caractéristiques, 15 - une probabilité qu'un individu passe de la zone du dernier évènement postérieur du trajet reconstitué à la zone de l'évènement postérieur résiduel sans avoir été détecté dans une zone intermédiaire. 4. Procédé selon la revendication 3, dans lequel 20 le produit permettant de calculer la probabilité que l'évènement postérieur résiduel suive le dernier évènement postérieur d'un trajet reconstitué comprend également la probabilité qu'un individu passe de la zone du dernier évènement postérieur du trajet reconstitué à la 25 zone de l'évènement postérieur résiduel dans le temps sé- parant l'information temporelle du dernier évènement postérieur du trajet reconstitué et l'information temporelle de l'évènement postérieur résiduel. 5. Procédé selon la revendication 3, dans lequel 30 la probabilité qu'un individu passe de la zone du dernier évènement postérieur du trajet reconstitué à la zone de l'évènement postérieur résiduel sans avoir été détecté est le produit de : - une probabilité d'erreur de détection, 35 - une probabilité qu'un individu passe de la zone 3033648 21 du dernier évènement postérieur du trajet reconstitué à la zone de l'évènement postérieur résiduel compte tenu d'un historique des trajets reconstitués. 6. Procédé selon la revendication 2, dans lequel, 5 lors des étapes ultérieures, à partir de la topographie sont recherchés des chemins pouvant amener un individu à l'évènement postérieur résiduel en partant du dernier évènement postérieur de chaque trajet reconstitué et une densité de probabilités est calculée pour l'ensemble de 10 ces chemins. 7. Procédé selon la revendication 6, dans lequel un chemin est négligé si la probabilité qu'un individu emprunte ce chemin pour passer de la zone du dernier évènement postérieur du trajet reconstitué à la zone de 15 l'évènement postérieur résiduel compte tenu d'un histo- rique des trajets reconstitués est inférieure à un seuil prédéterminé. 8. procédé selon la revendication 1, dans lequel le modèle temporel est basé sur un temps de parcours 20 moyen sur le trajet. 9. Procédé selon la revendication 1, dans lequel les moyens de détection comprennent au moins un capteur biométrique. 10. Procédé selon la revendication 9, dans lequel 25 les moyens de détection comprennent plusieurs capteurs biométriques, les capteurs biométriques étant de types différents pour capturer des caractéristiques biométriques de types différents. 11. Procédé selon la revendication 13, dans le- 30 quel les moyens de détection de la première zone sont disposés au voisinage d'un poste de contrôle d'accès dans lequel chaque individu doit présenter un document d'identification. 12. Procédé selon la revendication 1, dans lequel 35 les caractéristiques représentatives détectées dans 3033648 22 chaque zone sont mémorisées. 13. Dispositif de suivi d'individus dans un lieu ayant des zones équipées de moyens de détection d'individus et dont certaines sont inaccessibles directe- 5 ment depuis d'autres zones, le dispositif de suivi com- prenant une unité informatique de traitement qui est reliée aux moyens de détection et qui est agencée pour traiter des données transmises par les moyens de détection, l'unité de traitement ayant une mémoire contenant 10 une topographie du lieu et un modèle temporel mettant en relation un temps de parcours et au moins un trajet entre les zones et exécutant un programme de mise en oeuvre du procédé selon l'une quelconque des relations précédentes.