FR 3000557 A1 20140704 FR 1262998 A 20121231 Système de cartographies de déplacements 1. Domaine de l'invention L'invention a pour objet l'analyse de déplacements dans des espaces délimités. Plus particulièrement, l'invention porte sur un procédé et un système de cartographie de déplacements dans un espace délimité. 2. Art antérieur On estime à plusieurs centaines de millions le nombre de personnes disposant d'un téléphone intelligent (smartphone) dans le monde et ce chiffre ne cesse de croitre. Ces téléphones intelligents embarquent de nombreuses fonctionnalités et disposent de plusieurs systèmes de réception et de transmission de données sans fils. Par ailleurs, il peut être utile, pour de nombreux acteurs privés ou institutionnels, de connaître les déplacements effectués par des clients ou des visiteurs dans un lieu donné. Ces déplacements, une fois connus, permettent par exemple d'optimiser les surfaces de vente, de prévoir des flux de déplacement de groupes de personnes, etc. Ainsi, des entreprises ont eu l'idée de créer des systèmes de cartographie de déplacement en utilisant les téléphones intelligents en possession des utilisateurs. Des publications antérieures se rapportent à des systèmes de repérage ou de surveillance des mobiles, mais beaucoup exigent l'utilisation de matériel ou de logiciel particulier sur l'appareil mobile. Dans le domaine des achats, le brevet US 5,295,064 décrit une méthode de suivi d'acheteurs en ajoutant du matériel de géolocalisation pour caddies, qui non seulement rend les chariots beaucoup plus cher et n'est pas capable de suivre les consommateurs qui n'utilisent pas de chariot. Dans le domaine plus général de suivi des appareils mobiles, il y a également un certain nombre de publications antérieures, mais encore une fois ils nécessitent tous une certaine interaction avec l'appareil mobile. Ainsi, dans le brevet US 6,716,101, un signal d'appel est envoyé à un appareil mobile et sa réponse est utilisée pour calculer sa position. Cette méthode oblige à la transmission d'un signal complémentaire à destination du terminal, ce qui peut interférer avec d'autres communications. Dans le brevet US 6,246,882, un procédé d'utilisation d'appareils mobiles spécialisés qui transmettent un signal de position est décrit. Cette technique ne fournit une indication utile sur le mouvement d'un groupe de personnes que si un grand nombre de personnes de ce groupe ont des dispositifs spécialisés tels que décrits dans ce brevet. Pour sa part, la demande de brevet EP1779133 décrit un système qui ne nécessite normalement pas l'utilisation d'un matériel ou d'un signal spécifique du côté de l'utilisateur. En revanche, cette demande de brevet nécessite la transmission, par le terminal mobile, d'un identifiant, ce qui pose des problèmes de confidentialité des données. De plus, le réseau utilisé en pratique est le réseau GSM ou UMTS. Par ailleurs, cette demande de brevet ne décrit pas concrètement comment le système est implémenté. Les éléments décrits dans cette demande ne permettent pas en effet de mettre en oeuvre une technique réellement bien définie. On apprend de cette demande que le système comprend de multiples capteurs, qui sont disposés dans les surfaces à cartographiées. Le téléphone transmet son identifiant à ces capteurs, lesquels déterminent la position et la direction du téléphone en analysant les différences de phases des signaux reçus par le téléphone. Ainsi, les systèmes de l'art antérieur ne proposent pas de solution simple à mettre en oeuvre et qui garantisse l'anonymat des utilisateurs et la sécurité des données recueillies. 3. Résumé de l'invention L'invention ne présente pas ces inconvénients de l'art antérieur. Plus précisément, l'invention permet de disposer d'une solution simple et performante d'obtention de données de parcours de personnes dans un espace déterminé. Plus particulièrement, l'invention se rapporte à un procédé détermination de données représentatives d'un déplacement d'un terminal de communication, appelé particule, au sein d'un espace géographique prédéterminé, dit espace à cartographier. Le procédé caractérisé comprend : une phase d'obtention de données représentatives de puissances d'émissions d'au moins un signal de transmission de données au sein dudit espace à cartographier, délivrant un ensemble de données capturées ; une phase de filtrage des données dudit ensemble de données capturées en fonction d'au moins un paramètre de filtrage prédéterminé, délivrant un ensemble de données filtrées ; une phase de calcul, à l'aide dudit ensemble de données filtrées, de données représentatives d'au moins un déplacement effectué par ladite particule au sein dudit espace à cartographier. Ainsi l'invention permet de définir les parcours des terminaux (particules) dans un espace à cartographier de manière simple, sans nécessiter de dispositif particulier. En effet, en se basant sur l'obtention d'une puissance de signal émis, l'invention offre une solution simple au problème susmentionné. Selon un mode de réalisation particulier ladite phase d'obtention de données comprend une étape préalable de calibration d'un réseau de capteurs disposés au sein dudit espace à cartographier, comprenant, pour au moins un capteur une étape de calcul d'une puissance de réception de référence. Ainsi, en calculant une puissance de réception de référence, on édifie un référentiel de captation de données de puissance qui permet de capter n'importe quel signal de façon simple. Selon un mode de réalisation particulier ladite phase de filtrage des données comprend : une étape de tentative d'identification, en fonction dudit au moins un paramètre de filtrage prédéterminé, d'au moins une donnée représentative d'au moins une particule parasite au sein desdites données capturées ; et lorsque ladite étape de tentative d'identification délivre un résultat positif : une étape de suppression desdites données liées audit au moins une particule parasite, délivrant lesdites données filtrées. Ainsi, on élimine facilement des données non nécessaires ou jugées non pertinentes. Selon une caractéristique particulière ledit au moins un paramètre de filtrage prédéterminé appartient au groupe comprenant : une durée d'émission de signal ; un nombre d'occurrence de données identifiant une particule ; une puissance d'émission de signal. Selon un mode de réalisation particulier ladite phase de calcul, à l'aide dudit ensemble de données filtrées, de données représentatives d'au moins un déplacement effectué par ladite particule au sein dudit espace à cartographier comprend : une étape d'obtention de la liste des positions connues de ladite particule, ladite liste comprenant au moins deux positions successives ; pour chaque couple de positions successives de ladite liste, une étape de calcul d'un sous-chemin parcouru entre les deux positions de ce couple de positions successives, délivrant une liste de sous chemins ; une étape de concaténation des sous-chemins de la liste de sous chemins pour former un chemin global. Ainsi, on calcul des sous chemins entre les positions connues ou estimées de la particule. Ces sous chemins permettent par la suite d'obtenir un chemin global. Selon une caractéristique particulière ladite étape d'obtention de la liste des positions connues de ladite particule comprend : une étape de détermination d'un ensemble de zones probables de positionnement (PPZ) de ladite particule en fonction : desdites puissances d'émission dudit ensemble de données filtrées associées aux capteurs d'un réseau de capteurs ; d'une représentation matricielle dudit espace à cartographier ; une étape de détermination d'au moins deux positions probables de ladite particule, à au moins deux moments différents, en fonction d'une loi statistique et/ou de contraintes associées audit espace à cartographier. Dans d'autres modes de réalisation, l'invention se rapporte également aux capteurs et aux dispositifs permettant la mise en oeuvre des procédés préalablement décrits. Dans un autre mode de réalisation, l'invention se rapporte également à un système pour la mise en oeuvre de l'invention telle que décrite dans les procédés précédents. Il est entendu que des procédés indépendants de captures, par les capteurs, des données brutes et de traitement, par un serveur de traitement, des données brutes capturées font intégralement partie de l'invention. Selon un mode de réalisation particulier, l'invention se rapporte à un dispositif de détermination de données représentatives d'un déplacement d'un terminal de communication, appelé particule, au sein d'un espace géographique prédéterminé, dit espace à cartographier. Selon l'invention, il comprend : des moyens d'obtention de données représentatives de puissances d'émissions d'au moins un signal de transmission de données au sein dudit espace à cartographier, délivrant un ensemble de données capturées ; des moyens de filtrage des données dudit ensemble de données capturées en fonction d'au moins un paramètre de filtrage prédéterminé, délivrant un ensemble de données filtrées ; des moyens de calcul, à l'aide dudit ensemble de données filtrées, de données représentatives d'au moins un déplacement effectué par ladite particule au sein dudit espace à cartographier. Dans un autre mode de réalisation, l'invention se rapporte à un système de 15 détermination de données représentatives d'un déplacement d'un terminal de communication, appelé particule, au sein d'un espace géographique prédéterminé, dit espace à cartographier. Ce système comprend : au moins un capteur de puissance d'émissions d'au moins un signal de transmission de données au sein dudit espace à cartographier, délivrant un ensemble de données 20 capturées ; au moins un réseau de communication par lequel les données capturées sont transmises ; au moins un dispositif de détermination, connecté audit réseau de communication et comprenant : 25 des moyens d'obtention desdites données capturées ; des moyens de filtrage des données dudit ensemble de données capturées en fonction d'au moins un paramètre de filtrage prédéterminé, délivrant un ensemble de données filtrées ; des moyens de calcul, à l'aide dudit ensemble de données filtrées, de données représentatives d'au moins un déplacement effectué par ladite particule au sein dudit espace à cartographier. Selon une implémentation préférée, les différentes étapes des procédés selon l'invention sont mises en oeuvre par un ou plusieurs logiciels ou programmes d'ordinateur, comprenant des instructions logicielles destinées à être exécutées par un processeur de données d'un module relais selon l'invention et étant conçu pour commander l'exécution des différentes étapes des procédés. En conséquence, l'invention vise aussi un programme, susceptible d'être exécuté par un ordinateur ou par un processeur de données, ce programme comportant des instructions pour commander l'exécution des étapes d'un procédé tel que mentionné ci-dessus. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée, ou dans n'importe quelle autre forme souhaitable. L'invention vise aussi un support d'informations lisible par un processeur de données, et comportant des instructions d'un programme tel que mentionné ci-dessus. Le support d'informations peut être n'importe quelle entité ou dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen d'enregistrement magnétique, par exemple une disquette (floppy disc) ou un disque dur. D'autre part, le support d'informations peut être un support transmissible tel qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet. Alternativement, le support d'informations peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé en question. Selon un mode de réalisation, l'invention est mise en oeuvre au moyen de composants logiciels et/ou matériels. Dans cette optique, le terme "module" peut correspondre dans ce document aussi bien à un composant logiciel, qu'à un composant matériel ou à un ensemble de composants matériels et logiciels. Un composant logiciel correspond à un ou plusieurs programmes d'ordinateur, un ou plusieurs sous-programmes d'un programme, ou de manière plus générale à tout élément d'un programme ou d'un logiciel apte à mettre en oeuvre une fonction ou un ensemble de fonctions, selon ce qui est décrit ci-dessous pour le module concerné. Un tel composant logiciel est exécuté par un processeur de données d'une entité physique (terminal, serveur, passerelle, set-top-box, routeur, etc) et est susceptible d'accéder aux ressources matérielles de cette entité physique (mémoires, supports d'enregistrement, bus de communication, cartes électroniques d'entrées/sorties, interfaces utilisateur, etc). De la même manière, un composant matériel correspond à tout élément d'un ensemble matériel (ou hardware) apte à mettre en oeuvre une fonction ou un ensemble de fonctions, selon ce qui est décrit ci-dessous pour le module concerné. Il peut s'agir d'un composant matériel programmable ou avec processeur intégré pour l'exécution de logiciel, par exemple un circuit intégré, une carte à puce, une carte à mémoire, une carte électronique pour l'exécution d'un micrologiciel (firmware), etc. Chaque composante du système précédemment décrit met bien entendu en oeuvre ses propres modules logiciels. Les différents modes de réalisation mentionnés ci-dessus sont combinables entre eux pour la mise en oeuvre de l'invention. 4. Figures D'autres caractéristiques et avantages de l'invention apparaîtront plus clairement à la lecture de la description suivante d'un mode de réalisation préférentiel, donné à titre de simple exemple illustratif et non limitatif, et des dessins annexés, parmi lesquels : la figure 1 illustre un système de manière général ; les figure 2a, 2b et 2c sont des illustrations de la trilatération statistique de l'invention ; la figure 3 illustre la structure physique d'un capteur d'un réseau de capteur ; la figure 4 illustre la structure d'un serveur de traitement selon l'invention ; la figure 5 illustre un arbre de représentation des chemins les plus fréquents. 5. Description d'un mode de réalisation Dans ce mode de réalisation, on cherche à simplifier et à automatiser tant la capture des données que leurs exploitations ultérieures. Ceci est possible d'une part en segmentant les processus (obtention de données brutes puis corrélation de ces données brutes) et d'autre part en mettant en oeuvre les optimisations présentées ci-après. 5.1. Rappel du principe de l'invention Comme explicité, le système de l'invention offre de multiples avantages. Pour débuter, ce système ne nécessite pas de réglage particulier. Il est uniquement nécessaire de positionner des capteurs à de multiples endroits de l'espace à cartographier. Le système se calibre ensuite de manière automatique, sans nécessiter d'intervention humaine. Ensuite le système possède la capacité de s'auto-synchroniser, ce qui n'est pas envisagé par les systèmes antérieurs. En effet, l'une des problématiques associée à l'utilisation de multiples capteurs tient à la synchronisation de leurs horloges. Il est bien connu que les composants utilisés pour construire des dispositifs électroniques présentent des différences relativement importantes. Ainsi, même en réglant précisément l'heure de deux capteurs à un temps TO, des différences plus ou moins importantes sont observées entre les heures de ces deux capteurs au bout de quelques temps. Or, la synchronisation de ces capteurs est importante pour pouvoir déterminer la position d'un dispositif de manière assez précise. Les systèmes de positionnement par satellite GPS est soumis à des problématiques similaires. Cependant, à la différence du GPS, et pour des questions de coûts, le système objet de la présente demande de brevet ne peut pas mettre en oeuvre une horloge atomique dans chaque capteur. Par ailleurs, pour assurer la confidentialité des données obtenues, le système ne se base pas sur l'utilisation de données personnelles de l'utilisateur. Au contraire, dès l'origine, les données de l'utilisateur sont remplacées par des données anonymes. Enfin, le système comprend au moins une entité de traitement, laquelle assure l'organisation et la gestion des données. Plus précisément, en lien avec la figure 1, le système objet de la demande de brevet est mis en oeuvre de la façon suivante : Une première phase consiste en la calibration des capteurs du réseau de capteurs une fois qu'ils ont été positionnés. Une seconde phase consiste en la synchronisation initiale de ces capteurs. Une troisième phase consiste en la capture des données en provenance des dispositifs (ordinateurs, terminaux mobiles, tablettes) qui circulent dans l'espace à cartographier. Ces données sont anonymisées dès leurs captures. Il s'en suit une phase de synchronisation en ligne et un stockage des données brutes obtenues. La phase suivante consiste traiter ces données brutes pour permettre leurs interprétations. En d'autres termes, le système comprend d'une part des moyens d'obtention de données individuelles, brutes, qui sont anonymes (il s'agit du réseau de capteurs) et stockées dans un entrepôt de données brutes et des moyens de traitement et de corrélations de ces données individuelles pour former des données consolidées (c'est-à-dire représentatives des déplacements des particules) aptes à être interprétées. Plus particulièrement, en relation avec la figure 1, le système met en oeuvre un procédé détermination de données représentatives d'un déplacement d'un terminal de communication (T1), appelé particule, au sein d'un espace géographique prédéterminé, dit espace à cartographier (Sc). Le procédé de l'invention comprend : une phase d'obtention de données représentatives de puissances d'émissions (P1, P2, P3, etc.) d'au moins un signal de transmission de données au sein dudit espace à cartographier, délivrant un ensemble de données capturées, stockées dans une base de données (RDB), auquel le réseaux de capteurs (C1, C2, etc.) est connecté par l'intermédiaire d'un réseau de communication (NTWK) ; une phase de filtrage des données, par un dispositif de filtrage (PostT), dudit ensemble de données capturées en fonction d'au moins un paramètre de filtrage prédéterminé, délivrant un ensemble de données filtrées (ADB) stocké dans une base de données ; une phase de calcul, par un serveur de traitement (CLC), à l'aide dudit ensemble de données filtrées, de données représentatives d'au moins un déplacement effectué par ladite particule au sein dudit espace à cartographier, elles aussi stockées dans une base de données. Il est clair que d'autres modes de réalisations peuvent être envisagées sans qu'ils soient exclus de la portée de l'invention. Par exemple, la connexion à un réseau de communication ou l'utilisation de multiples bases de données peut être évitée en fonction des modes de réalisation. De même, il n'est pas obligatoire de disposer d'un serveur de traitement CLC et d'un dispositif de filtrage particulier PostT. Dans un mode de réalisation particulier de l'invention, un capteur particulier du réseau de capteur (appelé capteur « maitre ») peut tout à fait remplir ces fonctions. Auquel cas les données sont produites « sur site », c'est-à-dire au sien même de l'espace à cartographier. Le système décrit dans l'invention peut être utilisé dans de nombreuses applications dont certaines sont données ci-dessous à titre illustratif et nullement limitatif : Procédé utilisé pour localiser les détenteurs de téléphones mobiles présents dans un espace de vente, par exemple à l'intérieur d'un magasin, d'un centre commercial, etc. Procédé utilisé pour localiser et analyser les flux des visiteurs des grands centres administratifs, dans les gares, dans les aéroports ; Procédé utilisé pour localiser, analyser les zones d'affluence et les points de blocage dans les centres hospitaliers ; Procédé utilisé pour mesurer les files d'attente et les flux dans les salles d'exposition de musée. 5.2. Description d'un mode de réalisation 5.2.1. Capture des données brutes. La première phase du processus consiste en la capture de données brute. Cette capture permet d'obtenir des données de base qui seront traitées par la suite. Dans ce mode de réalisation de l'invention, les données brutes sont les suivantes : identifiant d'une particule (d'un terminal), heure de capture et puissance du signal capturé. Cependant, pour pouvoir obtenir ces données, plusieurs opération préalables doivent ou peuvent être effectuées en fonction des conditions d'utilisation et de paramètres liés aux capteurs. Les optimisations suivantes peuvent ou non être utilisées, en fonction du mode de réalisation de l'invention. Par exemple, si les horloges des capteurs sont parfaitement synchrones, il n'est pas nécessaire de réaliser une synchronisation. De même, si la puissance de réception à un mètre est connue, il n'est pas nécessaire de réaliser un auto calibrage des capteurs du réseau de capteurs. 5.2.1.1. Calibration du réseau de capteurs. Comme explicité, la première phase consiste à réaliser une auto-calibration des capteurs du réseau de capteurs. Ce calibrage est rendu nécessaire par le fait que la position d'un terminal est déterminée par sa distance avec le capteur qui enregistre sa position. Ainsi, si un terminal T1 est considéré comme étant à une distance dl d'un premier capteur Cl, une distance d2 d'un deuxième capteur C2 et une distance d3 d'un troisième capteur C3, il est possible de déterminer sa position dans l'espace à cartographier (voir par la suite). Cette technique générale est bien connue, mais a été améliorée dans le cadre de l'invention. Dans ce mode de réalisation, la détermination de la distance d entre un capteur et un terminal est réalisée au niveau du capteur lui-même. Or, le capteur ne mesure pas de distance mais capteur une puissance de signal émis par le terminal Tl. Il est donc nécessaire de transformer cette puissance en distance. La puissance de réception au niveau du récepteur, Pr, est définie par la fonction f suivante : Pr= f(Po, a, d) Avec : Po est la puissance de réception à une distance de 1 mètre ; «a » est une constante représentant l'atténuation ; et « d » la distance entre le récepteur et l'émetteur. Plus précisément, la formule est la suivante : Pr= P0-10. a.log(d). Pour déterminer la distance d'un terminal à un récepteur, les inventeurs ont donc eu l'idée d'inverser cette fonction pour déterminer « d », au niveau de chaque capteur. Cependant, pour que cela soit possible, il est nécessaire de connaître la valeur de Po et de a. « a » possède une valeur généralement admise pour les milieux fermés qui est donc utilisée. En revanche, il est nécessaire de déterminer Po. Une première méthode peut consister à mesurer manuellement cette valeur sur le lieu d'installation (dans l'espace à cartographier). Cette technique donne le résultat attendu mais elle est longue et compliquée particulièrement lorsque les capteurs sont installés en hauteur. Par voie de conséquence, cette technique est couteuse et non optimale. Une deuxième technique, objet de la présente demande de brevet, consiste à réaliser une auto calibration (ou auto-calibrage). Selon l'invention, l'auto-calibrage est réalisé grâce à un processus au cours duquel les capteurs écoutent et émettent alternativement sur de courtes périodes de temps. De façon aléatoire, un capteur CA aura écouté un capteur CB et émis vers le capteur CB et réciproquement. Au bout d'un certain temps, le capteur CA aura écouté au moins une fois tous les autres capteurs. Par exemple, dans le cas de six capteurs, le capteur CA aura écouté les capteurs {CB, CC, CD CE, CF} au moins une fois. Au cours de ce processus d'écoute, le capteur CA connaît la puissance des signaux reçus de la part des capteurs {CB, CC, CD CE, CF}. Il dispose donc d'au moins une mesure PAB de la puissance, qui est la puissance mesurée par le capteur CA du signal reçu du capteur CB. Or, dans le cadre de l'installation des capteurs, une donnée a été fournie aux capteurs : il s'agit de la distance qui sépare ce capteur d'au moins un autre capteur présent dans la zone à cartographier. En effet, au cours de l'installation des capteurs dans la zone à cartographier, l'utilisateur sait quelles sont les distances qui séparent les capteurs. Il peut donc fournir cette information aux capteurs. Ainsi, la résolution de la relation entre la puissance et la distance entre le capteur CA et le capteur CB (soit PABePoA, a, DAB)), permet de déterminer POA, la puissance de réception à 1 mètre du capteur CA de façon automatique. Les autres capteurs {CB, CC, CD CE, CF} agissent de même pour obtenir leur propre puissance de réception à 1 mètre. Les données de calibrage de l'ensemble des capteurs peuvent donc être déterminées de façon automatique sans aucune intervention humaine et notamment sans utilisation de matériel de calibrage complémentaire. 5.2.1.2. Synchronisation des Captures de données brutes. Une autre problématique est liée à la synchronisation des données enregistrées au sein des capteurs. Cette synchronisation est nécessaire pour la corrélation ultérieure des données individuelles. Les données individuelles sont enregistrées dans une structure de données de type fichier à plat. Cette structure de données de type « fichier à plat » est appelée fichier de capture. Un fichier de capture recense l'heure et la puissance de réception des émissions d'une particule. Une particule est par exemple un terminal d'un utilisateur se déplaçant dans l'espace à cartographier. Comme cela est explicité par la suite, lors de la phase de traitement des données, les différentes positions d'une particule dans l'espace à cartographier sont déterminées en corrélant les informations des différents capteurs par trilatération. Malheureusement, l'heure système d'un capteur donné n'est généralement pas identique à celle d'un autre capteur. Il en résulte que, pour une même particule, à un instant donné, les heures qui sont enregistrées sur les fichiers de capture de deux capteurs différents peuvent être distinctes de quelques secondes, rendant en pratique ces données inutilisables. Il faut détecter et corriger cet écart en synchronisant les captures. Afin de synchroniser les captures d'un même réseau de capteurs, il convient de désigner un capteur du réseau comme le capteur source dont l'heure est correcte et de décaler les capteurs adjacents pour s'aligner sur l'heure de ce capteur source. Par itération, on synchronise l'ensemble des capteurs. Dans ce mode de réalisation de l'invention, une méthode statistique est employée pour synchroniser les heures de capture. Cette méthode statistique est basée sur le fait qu'un grand nombre de données sont enregistrées. On suppose également, dans ce mode de réalisation, que chaque capteur a une horloge interne qui dévie de l'heure des autres capteurs de façon constante. Ainsi, par exemple, la déviation de l'heure indiquée par l'horloge HA du capteur CA est constante et la déviation de l'heure indiquée par l'horloge HB du capteur CB est constante. Ainsi, l'horloge HA a une déviation constante Dev A et l'horloge HB a une déviation constante Dev B. La relation entre les deux horaires des capteurs CA et CB est donc définie par la fonction : HB=Y + X x HA, avec : X = Dey B ± Dey A. Dans ce mode de réalisation, on suppose donc que la déviation entre deux horloges données est une fonction affine. On suppose en plus que le coefficient directeur de cette fonction est égal à un (X=1) ce qui correspond à une déviation identique de l'horloge HA et de l'horloge HB. Ceci n'est bien sûr qu'une supposition de ce mode de réalisation. D'autres modes de réalisation peuvent supposer une formulation différente. Ces hypothèses posées, la méthode mise en oeuvre consiste à sélectionner des particules de synchronisation (particules à occurrence diverses qui servent d'étalon de mesure) et à calculer des matrices des écarts entre les heures de captures des deux capteurs pour sélectionner dans cette matrice l'écart apparaissant le plus souvent. Les particules de synchronisation sont sélectionnées de la façon suivante : on définit des groupes de particules en fonction du nombre de captures qui sont considérées comme représentatives. On définit les groupes avec peu de captures ( 10 et 100) ; dans ces groupes déterminés, on sélectionne au hasard un nombre n prédéterminé de particules. Pour chaque particule ainsi sélectionnée, on dispose d'une liste de captures à partir d'un capteur CA et une liste de capture à partir d'un capteur CB. Chaque liste (une par capteur) comprend une pluralité d'enregistrements, un enregistrement d'une liste comprenant un triplet composé d'un identifiant de la particule, une heure de capture et une puissance de réception. À partir de deux listes (par exemple les listes issues des capteurs CA et CB), une matrice M est construite, comprenant les coefficients Mij = ai - tir avec : « a » et « b » les heures des capteurs. Plus précisément, « a » correspond à l'heure de capture de la particule par le capteur CA et « b » correspond à l'heure de capture de la particule par le capteur CB. Une fois cette matrice M construite, elle est parcourue et l'écart entre les horloges des capteurs CA et CB, pour cette particule, est déterminé comme étant la valeur Mu qui apparait le plus souvent dans la matrice M. On obtient ainsi, pour chaque particule sélectionnée, un écart calculé. Les inventeurs ont identifié qu'il était important de réaliser une matrice des différences car comme les horloges des capteurs ne sont pas synchronisées et que l'on ne sait pas quand les captures des capteurs sont faites les unes par rapport aux autres (il n'est pas assuré par exemple que la première capture du capteur CA soit celle qui corresponde à la première capture du capteur CB), il faut rechercher un écart à l'aveugle. La meilleure façon de rechercher l'écart à l'aveugle est justement de calculer tous les écarts possibles (dans la matrice) et d'identifier celui qui apparaît le plus souvent. L'écart final entre les deux horloges des capteurs CA et CB est celui apparaissant le plus souvent sur l'ensemble des particules. Ainsi, on obtient l'écart probable (Ep) entre CA et CB. Cet écart probable est par la suite soustrait aux captures du capteur CB (On considère que le capteur CA est le capteur de référence). La même procédure est appliquée aux autres capteurs à partir d'un capteur de référence (dans l'exemple il s'agit de CA). Lorsqu'un capteur de référence et un capteur courant ne partagent pas de données communes (car ils sont physiquement trop éloignés et que les données obtenues ne sont pas pertinentes), un autre capteur de référence est utilisé.30 5.2.1.3. Anonymisation des données. Afin de pouvoir déterminer la position et la parcours d'une particule donnée il est nécessaire d'identifier cette particules. Plusieurs types de données peuvent être utilisés pour pouvoir identifier une particule. En revanche, il est nécessaire d'assurer le caractère anonyme de ces identifiants. Dans ce mode de réalisation de l'invention, on choisit de récupérer les adresses MAC des terminaux (qui sont émises par les terminaux lorsqu'ils cherchent à s'accrocher à un réseau, par exemple un réseau de type wifi). Une particule correspond donc à un terminal et on utilise comme identifiant l'adresse MAC de celui-ci. Pour garantir l'anonymat, on met en oeuvre un « Salted hashing ». Pour assurer le caractère anonyme des données capturées, il s'agit d'ajouter au principe de hachage par SHA-1, un « salt » ou mot secret. En d'autres termes, dans ce mode de réalisation, l'adresse MAC est hashée une première fois, puis, on ajoute à la signature obtenue un mot secret auquel on applique à nouveau un hashage. La donnée d'identification enregistrée est donc une donnée anonyme. Il n'est pas possible d'identifier l'utilisateur auquel cette donnée est associée. La donnée d'identification est rendue anonyme dès l'origine, à savoir avant l'enregistrement du triplet composé de l'identifiant de la particule, de l'heure de capture et de la puissance de réception. 5.2.2. Filtrage des données. La deuxième phase du processus consiste à effectuer un filtrage des données capturées. En effet, parmi l'ensemble des données capturées, nombreuses sont celles qui ne sont pas pertinentes vis-à-vis de l'objet de l'invention. Il s'agit de données parasites émises par des dispositifs parasites. Ainsi, au moment du post-traitement (ou à un autre moment plus pertinent en fonction de mode de réalisation) les données sont filtrées. De façon automatique, sont exclues les données parasites provenant : des émetteurs fixes préalablement déterminés. Ces émetteurs sont repérés par la réception d'un beacon (paquet de données émis par un émetteur fixe). Chaque donnée provenant d'un de ces émetteurs fixes est éliminée de l'analyse. Les émetteurs fixes exclus sont par exemple les routeurs WiFi. Un émetteur fixe émet en permanence mais à une forme de communication différente des émetteurs permanents car c'est un fournisseur d'accès pour les autres communicants WiFi. des émetteurs dits « permanents » repérés en fonction de la durée d'émission des ondes. Si cette durée est supérieure à un laps de temps déterminé, les données seront exclues de l'analyse. Les émetteurs permanents peuvent être : les téléphones des vendeurs, les ordinateurs, tablettes, imprimantes sans fil, télévisions connectées en démonstration dans l'espace étudié, etc. Cette liste n'est pas exhaustive. Pour que ces émetteurs permanents ne soient pas pris en compte dans le calcul, ils sont donc filtrés en fonction du temps de présence. Dans ce mode de réalisation, un émetteur dont la durée de présence est supérieure à deux ou trois heures est considéré comme permanent et est exclu. Cette valeur est bien entendue variable en fonction d'une part de la surface à cartographier et d'autre part en fonction du nombre d'émetteurs permanents. En fonction du mode de réalisation, les émetteurs fixes sont filtrés au moment du « post-traitement » alors que les émetteurs permanents sont filtrés lors du « traitement de données ». Ceci est variable et dépend d'une part de l'autonomie des capteurs du réseau de capteurs (sont-ils capables de prendre des décisions) et d'autre part de la configuration générale du système (a-t-on par exemple un traitement centralisé des données). 5.2.3. Traitement des données brutes. Comme indiqué préalablement, le traitement des données consiste à corréler les données brutes obtenues et à leur donner un sens. En effet, il ressort des techniques précédemment décrites dans les documents de l'art antérieur, par exemple, qu'un des postulats implicites sur lesquels ces techniques antérieures se basent est que la localisation de l'utilisateur peut être connue en permanence. Or, il est clair que cela ne peut pas être le cas, et ce pour de nombreuses raisons (par exemple, de par la présence d'interférences ou dans le cas ou des dispositifs sont ajoutés à des chariots, parce que l'utilisateur peut se déplacer dans la surface de vente sans son chariot). Ainsi, il est nécessaire de prévoir des méthodes qui d'une part permettent d'assurer une localisation partielle des particules dans l'espace à cartographier et d'autre part permettent de définir un réel chemin parcouru par la particule alors même que l'on ne dispose à son sujet que d'informations incomplètes.30 5.2.3.1. Positionner les particules par trilatération La méthode de positionnement par trilatération est une technique connue qui nécessite au moins trois informations de distance pour déterminer une localisation. En revanche, les inventeurs ont mis au point une technique qui ne nécessite qu'une information partielle pour l'obtention de la position probabiliste de la particule. En effet, les données brutes récupérées lors de la capture des données brutes comprennent la puissance du signal qui définit la distance entre le capteur et la particule. Cette méthode est illustrée par les figures 2a, 2b et 2c. Quand, à un instant donné, une particule est vue par au moins trois capteurs, par trilatération, la position de la particule est déterminée. Les inventeurs ont mis en place une méthode permettant de définir la position d'une particule repérée par un, deux, trois ou plus de capteurs. Pour ce faire, la méthode de trilatération est complétée par l'usage de données complémentaires formées d'une part par des contraintes systémiques liées à l'espace à cartographier (telles que la position d'une cloison par exemple, d'une zone infranchissable, etc.) et d'autre part par un pavage de l'espace à cartographier. En d'autres termes, préalablement à l'application de la trilatération, l'espace à cartographier est découpé en cases. Une matrice représente donc l'espace à cartographier. Au sein de cette matrice, un certain nombre de cases sont « indisponibles » : il s'agit par exemple des cases correspondant à des emplacements de murs ou d'étalages (quand il s'agit d'un magasin). Ainsi, en lieu et place des données qui peuvent manquer pour le positionnement d'une particule par trilatération, les inventeurs ont eu l'idée d'utiliser les données de la carte elle-même pour positionner la particule. Plus précisément, la détermination de la position est mise en oeuvre de la façon suivante : On appelle « PPZ » les zones de présence déterminées. À un instant donné, chaque émetteur possède une liste de n-uplets comprenant une puissance associée à chaque particule; : Capteur 1 = { P11; P12; P13 ; ... ; P1n1 } Capteur 2 = { P21 ; P22; P23 ; ... ; P2n2}30 Capteur k = Pk1 ; Pk2 ; Pk3 ; ; Pknk Pour chaque capteur la puissance des émissions de la particule vers ce capteur est fonction de la distance les séparant, D : = Po, - 10 * a * logio (D) + Où P. est l'atténuation à 1 mètre de la puissance de réception, a est le facteur d'atténuation de cette puissance, et X' est une variable aléatoire de loi Normale possédant une moyenne nulle et un écart-type a. Ce dernier terme dans l'équation permet de prendre en compte les variations dans les atténuations dues à l'hétérogénéité du milieu de transmission (présence d'obstacles, de murs, etc.). Pour un capteur l'estimation de la distance entre ce capteur et une particule est effectué en considérant un échantillon de n (n>1) puissances {Pi1 ; Pi2 ; ; Pin} mesurées à un moment donné. Pour cela, il est nécessaire d'estimer à partir de ces puissances une moyenne : 1 Pm = -n* 1Pij j=1 et une variance, qui est la même que celle de la variable aléatoire X, : n 1 i=0 - Po2 S0-2 - n - 1 À partir de cette estimation, il est possible de calculer un intervalle de confiance à 95% de la variable aléatoire X, : [ xcrinf ,xusupi, xcrinf xusup == * (t951 est le quantile d'ordre 95 à n-1 degré de liberté, et est utilisé afin d'estimer l'intervalle de confiance sur des échantillons de petites tailles). Il devient alors possible d'estimer un intervalle pour la distance séparant la particule du capteur i. En effet, en inversant la formule ci-dessus fournissant la puissance à partir de la distance, il est possible de borner cette distance D : P .-Prri+Xianf entre une distance minimale : Dinf = 10( 10*a ) - Po i-Prn-FXS0-UP et une distance maximale :Dsup = 10( io*a La distance minimum et la distance maximum au capteur permet de déterminer un anneau dans lequel est positionnée la particule. Par ailleurs, le plan (l'espace à cartographier) est découpé en «°cases°». Les anneaux définis par les distances min et max définissent alors un ensemble de cases probables pour chacune des captures. Chaque case correspond à une information de présence représentée par un bit : - Case = { 0 } si la particule ne peut être présente dans la case ; - Case = { 1} si la particule peut être dans la case. Avec un capteur, un ensemble de cases représentant une partie de l'anneau permet de positionner la particule de façon probabiliste. Avec 2 ou 3 capteurs, l'intersection des cases probables permet de définir précisément la position de la particule. Plus la taille des cases est petite plus le positionnement sera précis. Il en est de même pour le nombre de capteurs. Les schémas des figures 2a, 2b et 2c représentent trois cas de positionnement de particules avec un, deux et trois capteurs. Dans la figure 2a, un seul capteur est utilisé. Une zone de présence (PPZ1), est constituée d'un arc comprenant un certain nombre de cases probables. Dans la figure 2b, on a deux zones de présences (PPZ1, PPZ2) qui permettent la détermination d'une première intersection (PPZ1+2) correspondant à un certain nombre de cases dans laquelle la particule est susceptible de se trouver. Dans la figure 2c, on a successivement la détermination de trois zones de présence (PPZ1, PPZ2 et PPZ3) qui conduisent à la détermination d'une deuxième intersection (PPZ1+2+3) plus réduite. Une fois ces zones de présence déterminées, les conditions particulières de l'espace à cartographier sont appliquées (comme par exemple les cloisons, les espaces de stockages, les murs, etc.), pour limiter encore les possibilités de présence de la particule. Puis une loi normale est appliquée sur les cases restantes pour déterminer la position probable de la particule (voir section suivante). 5.2.3.2. Déterminer des chemins Une fonctionnalité importante du système consiste à déterminer les chemins qui sont empruntés par les particules et donc à cartographier l'espace. Cette cartographie est effectuée sur la base des positions successives des particules dans l'espace. Plus particulièrement, à l'issue de la phase de trilatération, phase dans laquelle les positions approximatives des particules sont obtenues, il est nécessaire d'une part d'affiner les positions relatives et d'autre part de rechercher les différents chemins parcourus. On note que cette détermination est réalisée dans un milieu dit « discret », c'est-à-dire que ce milieu est découpé en zones (qui sont des cases d'une matrice). Cette détermination n'est pas réalisée dans un milieu « continu », qui nécessite une connaissance en temps réel de la position des particules, ce qui n'est pas le cas dans ce mode de réalisation. 5.2.3.2.1. Positionnement des particules (détermination des positions successives) Chaque particule est reliée à un ou plusieurs zones de présence « PPZ » qui sont référencés par un temps. Du fait de la procédure décrite ci-haut, la zone de présence d'une particule à un moment donné peut être constituée de plusieurs cases. Il est nécessaire pour la reconstitution des chemins de disposer d'un seul point (ou une seule case) à chaque instant. Pour cela de présence est réduite, en sélectionnant parmi les cases potentielles, celle où la particule a la plus grande probabilité de se trouver. Étant donné que la zone de présence est le résultat d'intersection d'anneaux, que dans chacun de ces anneaux la probabilité d'être au centre de l'anneau est beaucoup plus importante que celles d'être aux bords de celui-ci (du fait de l'utilisation d'une loi Normale pour le calcul de ces anneaux), la réduction de la zone de présence est effectuée en choisissant la case la plus proche du « centre » de la zone. Une fois le nombre de points de présence réduit, la méthode d'obtention du chemin consiste à relier l'ensemble de ces points entre eux pour définir le chemin de la particule (effectué dans la phase suivante). À l'issue de cette phase de positionnement des particules, on obtient une liste des positions connues qui comprend des enregistrements comprenant chacun un temps et une position, repérée par une case dans la matrice de l'espace à cartographier. 5.2.3.2.2. Détermination du (sous) chemin emprunté entre deux positions À l'issue de la phase précédente de détermination des positions successives des positions des particules, il est nécessaire d'estimer, pour chaque particule, les sous-chemins empruntés entre les positions connues pour obtenir un chemin « global » de la particule entre son entrée dans l'espace à cartographier et sa sortie de l'espace à cartographier. Cette obtention du chemin global est réalisée par une suite de calculs des sous-chemins empruntés entre les suites de positions connues. On a ainsi une itération pour l'obtention d'un chemin global qui met en oeuvre des calculs de sous-chemins entre deux positions successives connues. La méthode générale est la suivante et comprend : une étape d'obtention de la liste des positions connues de la particule (cette liste des positions connues comprend, on le rappelle, des enregistrements comprenant chacun un temps et une position, repérée par une case dans la matrice de l'espace à cartographier, la liste étant triée par heure croissante) ; pour chaque couple de positions connues d'une particule, une étape de calcul d'un sous-chemin parcouru entre ces deux positions. Par exemple, pour une liste comprenant trois positions successives (pl, p2 et p3), on calcule : le sous-chemin sc/ entre la position pl et la position p2 ; puis le sous-chemin sc2 entre la position p2 et la position p3. une étape de concaténation des sous-chemins pour former un chemin global. Dans cette méthode générale, la détermination des sous-chemins est l'étape la plus consommatrice en termes de temps et de ressources. Pour mettre en oeuvre cette étape de calcul d'un sous-chemin, les inventeurs ont eu l'idée d'utiliser deux méthodes différentes. Ces deux méthodes sont des méthodes elles-mêmes itératives. Entre une position de départ D et une position d'arrivée A (les deux positions étant respectivement identifiées par une case de départ D et une case d'arrivée A), le sous-chemin est déterminé par le calcul de plusieurs positions successives. Le nombre de positions calculées dépend de la longueur du sous-chemin, et donc dépend du nombre de cases de la matrice qui séparent la case de départ D de la case d'arrivée A. Dans un premier mode de réalisation, pour déterminer un sous chemin emprunté, un algorithme de découverte du plus court chemin est implémenté. Il existe plusieurs algorithmes de ce type dans la littérature (par exemple l'algorithme de Dijkstra). Dans ce mode de réalisation de l'invention, l'algorithme utilisé est A* (A Star). Cet algorithme permet de parcourir un plus court chemin entre différents points de présence. Pour effacer les effets de bord (superposition de trajectoires simultanées par exemple) qui ne correspondent que rarement à des comportements humains, A* (Astar) arrondi est utilisé. Lors du calcul des chemins des particules, chaque case est pondérée (avec une valeur donnée) en fonction du nombre de passage à chaque instant dans cette case. Une particule passant par une case à un instant donné, augmentera la pondération de cette case. Lors de la détermination de la prochaine case à utiliser pour leurs chemins, les particules ne peuvent emprunter les cases dont le poids dépasse un certain seuil (Su). De plus, le poids associé aux cases décroit en fonction du temps (la valeur donnée à cette case diminue en fonction du temps). Cette procédure permet de simuler la densité de passage par chaque case et d'obtenir des trajectoires de particules qui s'évitent. Les points de départ et d'arrivée, pour une particule, sont attribués en définissant l'entrée / la sortie de l'espace le plus proche du premier / dernier point de présence connu. Dans un deuxième mode de réalisation, pour déterminer un sous chemin emprunté, une méthode de « Flooding » est utilisée. Le flooding est utilisée dans le monde des jeux vidéo. L'objectif originel de cette méthode est de permettre à plusieurs objets en mouvement, ayant chacun sa position propre, d'atteindre une destination commune (Exemple : plusieurs entités dispersées voulant atteindre un même point où se trouve une autre entité). Le flooding consiste à noter les cases d'un environnement en fonction de la position de cette case par rapport à une certaine destination. Les cases adjacentes à la destination auront une note de 1. Les cases adjacentes aux cases notées 1 auront une note de 2, etc. Le flooding permet de déterminer la prochaine case de déplacement à chaque instant. La case de déplacement suivante est choisie en sélectionnant la case dont la note est la plus petite autour de la case courante, comme décrit dans le tableau ci-dessous: 5 5 5 5 5 6 7 8 9 9 4 4 4 4 5 6 8 D 8 3 3 3 5 6 6 7 8 2 2 2 6 5 5 6 7 8 1 1 1 6 4 5 6 7 8 1 1 3 4 5 6 7 8 1 1 1 2 3 4 5 6 7 8 Avec A la case d'Arrivée D la case de Départ (ou case courante) Et le poids de chaque case Les cases gris clair du tableau représentent des obstacles et les cases dont le poids est en gras définissent le trajet utilisé entre le départ D et l'arrivée A. À la différence de l'algorithme utilisé traditionnellement, le flooding est appliqué par les inventeurs avec n couples (départ, arrivée) simultanément, ce qui permet d'obtenir des trajectoires de l'ensemble des particules simultanément en prenant en compte les interactions entre ces particules. Ainsi, dans ce mode de réalisation, on met en oeuvre un processus itératif complexe, permettant de déterminer pas à pas les différents chemins empruntés par les particules. Dans le monde des jeux vidéo, plusieurs particules ont une seule destination (il y a un seul point destination est recherché par plusieurs entités). L'implémentation du flooding réalisée par les inventeurs associe une destination à une unique particule. De ce fait, il est nécessaire d'associer à chaque couple (case de départ, case d'arrivée) une matrice de distance. De plus, il est nécessaire d'enregistrer dans une matrice supplémentaire, la 'capacité' de chaque case (appelée matrice de Capacité). La capacité d'une case représente le nombre maximal de particules pouvant s'y trouver simultanément. Toutes les cases de cette dernière carte sont initialisées avec une valeur fonction de la taille de la carte représentée par la matrice (la capacité d'une case de 1 mètre * 1 mètre est de 2, par exemple, car il n'est pas possible d'y loger plus de deux piétons simultanément). De plus, chaque particule posée sur sa case de départ décroit la capacité de ladite case de 1. L'ensemble des matrices de distance par couple (départ, arrivée), ainsi que la matrice des capacités sont maintenues et modifiées à chaque itération de la façon suivante : Pour chaque particule non encore traitée à cette itération, choisir parmi les cases adjacentes à la case courante, la case de coût minimum, et dont la capacité n'a pas encore atteint 0. La particule se déplace alors sur cette case, accroit de 1 la capacité de la case qu'elle vient de quitter, et décroit de 1 la capacité de la case sur laquelle elle vient d'atterrir. Si pour une certaine itération, une particule ne peut pas se déplacer (toutes les cases adjacentes sont interdites), alors la particule reste immobile pour cette itération et ne modifie donc pas la carte des capacités. L'objectif étant de prendre en compte pour chaque particule l'environnement, et les possibles interactions avec d'autres particules pour atteindre sa destination. Le processus est ainsi itéré jusqu'à avoir déterminé l'ensemble des chemins de toutes les particules pour une période temporelle donnée. La période temporelle considérée peut être la journée, ou une fraction de journée. Bien entendu, ces modes de réalisations sont combinables et/ou superposables. Il peut ainsi être intéressant de réaliser deux calculs avec ces deux modes de réalisation pour un même sous chemin et de déterminer quel est le sous chemin le plus efficient à l'issue de ces deux calculs. Par ailleurs, en fonction de la topologie du terrain, il peut être intéressant d'appliquer plutôt un calcul pour un sous chemin donné et un autre calcul pour le sous chemin suivant. 5.3. Autres caractéristiques et avantages. 5.3.1. Caractéristiques complémentaires. 5.3.1.1. Multi-étages Pour reconstituer des chemins entre différents étages d'un espace, des « transit » sont déterminés entre les étages, c'est-à-dire des points d'entrées et de sorties reliant des étages d'un magasin. Ces points de transit permettent de dessiner des chemins entre différents étages : Du point de départ D au point d'arrivée A, la particule sera incitée à prendre le transit entre les deux étages et un seul chemin sera déterminé. 5.3.1.2. Multi-visites Un visiteur qui vient plusieurs fois dans la journée n'est pas récurrent pour autant. Il convient donc de déterminer les multi-visites, c'est-à-dire quand une particule vient à plusieurs reprises dans une même journée dans l'espace étudié. Si deux PPZ1 et PPZ2 successifs d'un même visiteur ont des intervalles de temps supérieurs à une durée seuil, alors les PPZ avant PPZ1 et ceux après PPZ2 seront séparés en deux visites. Ainsi, deux chemins seront calculés et correspondront à deux visites différentes du même visiteur. 5.3.1.3. Analyse des chemins les plus utilisés Afin de déterminer quels sont les chemins les plus utilisés, les inventeurs ont eu l'idée de superposer l'ensemble de chemins des particules. Un arbre des différentes possibilités de chemin est établi. La première branche de l'arbre représente les entrées dans l'espace, la dernière les sorties. Chaque branche de l'arbre cumule un poids de 1 à chaque passage d'une particule par cette branche. D'un point de vue algorithmique, un noeud d'entrée et un noeud de sortie sont ajoutés avec des poids de 0. Les poids permettent donc de déterminer le chemin le plus parcouru de l'entrée à la sortie de l'espace en prenant le poids le plus lourd à chaque pas. Le principe est schématisé dans l'arbre présenté en relation avec la figure 5 (Avec E : Entrée ; P : Porte ; S : Sortie). Par itération, pour trouver le 2ème chemin le plus utilisé, on retire le dernier poids des branches du chemin le plus utilisé de toutes les branches de l'arbre correspondant à ce chemin. On peut ainsi déterminer le 2ème, le 3ème, etc. chemins les plus utilisés en utilisant des itérations successives du même algorithme. D'autres méthodes de détermination du chemin le plus utilisé peuvent également être mises en oeuvre. 5.3.2. Mises en oeuvre par des dispositifs. La présentation du système dans lequel l'invention prend place a principalement été décrite en termes de méthode. Cependant, le système précédemment décrit repose sur la mise en oeuvre de dispositifs que l'on peut également décrire. La figure 3 illustre une structure simplifiée d'un capteur Cl d'un réseau de capteurs selon un mode de réalisation de l'invention. Par exemple, le capteur comprend une mémoire 31 constituée d'une mémoire tampon, une unité de traitement 32, équipée par exemple d'un microprocesseur, et pilotée par le programme d'ordinateur 33, mettant en oeuvre un procédé de capture de données. À l'initialisation, les instructions de code du programme d'ordinateur 33 sont par exemple chargées dans une mémoire avant d'être exécutées par le processeur de l'unité de traitement 32. L'unité de traitement 32 reçoit en entrée au moins un signal de transmission de données par exemple émis à partir d'un émetteur sans fil. Le microprocesseur de l'unité de traitement 32 met en oeuvre les étapes du procédé de capture de données, selon les instructions du programme d'ordinateur 33, pour se calibrer, réaliser un filtrage de données, et rendre des données capturées anonyme, lesquelles données sont transmises à un dispositif du réseau auquel le capteur est lui-même connecté. Pour cela, le capteur comprend, outre la mémoire tampon 31, des moyens de calibration, qui sont par exemple un module logiciel de calibration ou un processeur prévu ç cet effet, des moyens de capture de données, qui sont par exemple un émetteur transmetteur sans fil (type carte wifi ou GSM) ou un module logiciel prévu à cet effet. Ces moyens utilisent l'interface d'entrée I. Le capteur Cl comprend aussi : des moyens d'anonymisation des données capturées, qui peuvent être des composants de sécurité ou un module logiciel prévu à cet effet ; des moyens de transmission de ces données par l'interface de sortie T. Ces moyens sont pilotés par le microprocesseur de l'unité de traitement 32. La figure 4 illustre une structure simplifiée d'un serveur de traitement CLC selon un mode de réalisation de l'invention. Par exemple, le serveur de traitement comprend une mémoire 41 constituée d'une mémoire tampon, une unité de traitement 42, équipée par exemple d'un microprocesseur, et pilotée par le programme d'ordinateur 43, mettant en oeuvre un procédé de traitement de données conforme à ce qui a été présenté préalablement. À l'initialisation, les instructions de code du programme d'ordinateur 43 sont par exemple chargées dans une mémoire avant d'être exécutées par le processeur de l'unité de traitement 42. L'unité de traitement 42 reçoit en entrée des données d'indexation d'applications logicielles. Le microprocesseur de l'unité de traitement 42 met en oeuvre les étapes d'un procédé de contrôle d'une application logicielle, selon les instructions du programme d'ordinateur 43, pour réaliser des calculs de chemins empruntés par les différentes particules qui ont traversées l'espace à cartographier. Pour cela, le serveur de traitement CLC comprend, outre la mémoire tampon 41, des moyens de filtrage des données, se présentant sous la forme d'un module logiciel ou matériels de type FCPGA, des moyens de synchronisations temporelles des données en fonction des capteurs qui ont capturées ces données, se présentant sous la forme d'un module logiciel ou matériels de type FCPGA, se présentant sous la forme d'un module logiciel ou matériels de type FCPGA, des moyens de trilatération des données synchronisée pour déterminer des zones de présences probables des particules, se présentant sous la forme d'un module logiciel ou matériels de type FCPGA, des moyens de positionnement des particules, se présentant sous la forme d'un module logiciel ou matériels de type FCPGA, des moyens de calcul de sous chemin entre des successions de deux positions connues des particules afin de déterminer des chemins globaux des particules, se présentant sous la forme d'un module logiciel ou matériels de type FCPGA. Ces moyens sont pilotés par le microprocesseur de l'unité de traitement 42.10 REVENDICATIONS1. Procédé de détermination de données représentatives d'un déplacement d'un terminal de communication, appelé particule, au sein d'un espace géographique prédéterminé, dit espace à cartographier, procédé caractérisé en ce qu'il comprend : une phase d'obtention de données représentatives de puissances d'émissions d'au moins un signal de transmission de données au sein dudit espace à cartographier, délivrant un ensemble de données capturées ; une phase de filtrage des données dudit ensemble de données capturées en fonction d'au moins un paramètre de filtrage prédéterminé, délivrant un ensemble de données filtrées ; une phase de calcul, à l'aide dudit ensemble de données filtrées, de données représentatives d'au moins un déplacement effectué par ladite particule au sein dudit espace à cartographier. 2. Procédé selon la revendication 1, caractérisé en ce que ladite phase d'obtention de données comprend une étape préalable de calibration d'un réseau de capteurs disposés au sein dudit espace à cartographier, comprenant, pour au moins un capteur une étape de calcul d'une puissance de réception de référence. 3. Procédé selon la revendication 1, caractérisé en ce que ladite phase de filtrage des données comprend : une étape de tentative d'identification, en fonction dudit au moins un paramètre de filtrage prédéterminé, d'au moins une donnée représentative d'au moins une particule parasite au sein desdites données capturées ; et lorsque ladite étape de tentative d'identification délivre un résultat positif : une étape de suppression desdites données liées audit au moins une particule parasite, délivrant lesdites données filtrées. 4. Procédé selon la revendication 1, caractérisé en ce que ledit au moins un paramètre de filtrage prédéterminé appartient au groupe comprenant : une durée d'émission de signal ; un nombre d'occurrence de données identifiant une particule ; une puissance d'émission de signal. 5. Procédé selon la revendication 1, caractérisé en ce que ladite phase de calcul, à l'aide dudit ensemble de données filtrées, de données représentatives d'au moins un déplacement effectué par ladite particule au sein dudit espace à cartographier comprend : une étape d'obtention de la liste des positions connues de ladite particule, ladite liste comprenant au moins deux positions successives ; pour chaque couple de positions successives de ladite liste, une étape de calcul d'un sous-chemin parcouru entre les deux positions de ce couple de positions successives, délivrant une liste de sous chemins ; une étape de concaténation des sous-chemins de la liste de sous chemins pour former un chemin global ; 6. Procédé selon la revendication 5, caractérisé en ce que ladite étape d'obtention de la liste des positions connues de ladite particule comprend : une étape de détermination d'un ensemble de zones probables de positionnement (PPZ) de ladite particule en fonction : desdites puissances d'émission dudit ensemble de données filtrées associées aux capteurs d'un réseau de capteurs ; d'une représentation matricielle dudit espace à cartographier ; une étape de détermination d'au moins deux positions probables de ladite particule, à au moins deux moments différents, en fonction d'une loi statistique et/ou de contraintes associées audit espace à cartographier. 7. Dispositif de détermination de données représentatives d'un déplacement d'un terminal de communication, appelé particule, au sein d'un espace géographique prédéterminé, dit espace à cartographier, caractérisé en ce qu'il comprend : des moyens d'obtention de données représentatives de puissances d'émissions d'au moins un signal de transmission de données au sein dudit espace à cartographier, délivrant un ensemble de données capturées ; des moyens de filtrage des données dudit ensemble de données capturées en fonction d'au moins un paramètre de filtrage prédéterminé, délivrant un ensemble de données filtrées ; des moyens de calcul, à l'aide dudit ensemble de données filtrées, de données représentatives d'au moins un déplacement effectué par ladite particule au sein dudit espace à cartographier. 8. Système de détermination de données représentatives d'un déplacement d'un terminal de communication, appelé particule, au sein d'un espace géographique prédéterminé, dit espace à cartographier, caractérisé en ce qu'il comprend : au moins un capteur de puissance d'émissions d'au moins un signal de transmission de données au sein dudit espace à cartographier, délivrant un ensemble de données capturées ; au moins un réseau de communication par lequel les données capturées sont transmises ; au moins un dispositif de détermination, connecté audit réseau de communication et comprenant : des moyens d'obtention desdites données capturées ; des moyens de filtrage des données dudit ensemble de données capturées en fonction d'au moins un paramètre de filtrage prédéterminé, délivrant un ensemble de données filtrées ; des moyens de calcul, à l'aide dudit ensemble de données filtrées, de données représentatives d'au moins un déplacement effectué par ladite particule au sein dudit espace à cartographier. 9. Produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou stocké sur un support lisible par ordinateur et/ou exécutable par un microprocesseur, caractérisé en ce qu'il comprend des instructions de code de programme pour la mise en oeuvre du procédé de détermination selon l'une quelconque des revendications 1 à 6.