FR 3079309 A1 20190927 FR 1852429 A 20180321 DOMAINE TECHNIQUE La présente invention appartient au domaine de la géolocalisation. Notamment, l'invention concerne un procédé et un système de géolocalisation d'un terminal d'un système de communication sans fil. L'invention s'applique particulièrement bien à la géolocalisation d'objets communicants de type « Internet des Objets » (de l'anglais « Internet Of Things » ou loT) qui ont une propension à se déplacer en groupe. ÉTAT DE LA TECHNIQUE Au cours des dernières années, l'utilisation croissante des systèmes de communication sans fil a naturellement conduit à développer des services basés sur la position géographique d'un objet et utilisés par exemple pour l'aide à la navigation, la gestion du trafic, le suivi d'acheminement de marchandises, etc. Les systèmes de positionnement par satellites comme le GPS (« Global Positioning System ») comptent parmi les techniques de géolocalisation les plus connues. Ces systèmes reposent sur l'exploitation par un terminal récepteur de signaux radio émis par des satellites dédiés. La géolocalisation par satellite est particulièrement précise, mais elle souffre de plusieurs inconvénients, notamment le coût et la consommation électrique induits par l'intégration d'un récepteur GPS dans un objet, ainsi que son manque de performance dans des zones fermées. D'autres techniques existent pour géolocaliser un terminal à partir de signaux échangés avec les stations de base d'un réseau d'accès auquel il est connecté. Dans les réseaux cellulaires, comme par exemple le GSM (« Global System for Mobile Communications »), il est connu d'estimer la position d'un terminal comme étant celle de la station de base à laquelle il est actuellement associé. Cette méthode présente cependant une précision de géolocalisation médiocre puisque la zone de couverture d'une station de base peut atteindre plusieurs kilomètres, voire plusieurs dizaines de kilomètres de rayon. D'autres méthodes consistent à estimer les distances qui séparent un terminal de plusieurs stations de base en calculant des temps d'arrivée, des angles de propagation ou des différences de fréquence pour les signaux 2 échangés entre le terminal et les stations de base. Ces différentes méthodes présentent cependant toutes l'inconvénient de nécessiter du matériel et du logiciel spécifiques. Elles requièrent en effet généralement une synchronisation des terminaux et/ou des différentes stations de base jouant le rôle de points 5 d'observation. En outre, elles sont particulièrement sensibles au phénomène dit de multi-trajet (propagation d'un même signal radio par plusieurs chemins à cause des phénomènes de réflexion, réfraction et diffraction sur les obstacles rencontrés). D'autres méthodes de géolocalisation se basent sur le niveau de 10 puissance reçue (RSSI, de l'anglais « Received Signal Strength Indicator ») d'un signal échangé entre un terminal et une station de base. Ces méthodes sont particulièrement bien adaptées aux systèmes de communication sans fil de type réseaux cellulaires (comme par exemple le GSM) pour lesquels l'information de RSSI est directement disponible car elle est exploitée par le 15 système de communication lui-même. Ces méthodes reposent sur le fait qu'un signal radio est atténué dans l'atmosphère et donc que le niveau de RSSI d'un signal reçu par un récepteur varie en fonction de la distance qui sépare le récepteur de l'émetteur du signal. Ainsi, il est possible de déterminer la position géographique d'un terminal par trilatération en estimant la distance séparant le 20 terminal des différentes stations de base qui l'entourent à partir des niveaux de RSSI mesurés par les stations de base. L'inconvénient d'une telle méthode de géolocalisation par trilatération basée sur les niveaux de RSSI est son manque de précision dû au fait que les nombreux paramètres qui influent sur l'atténuation du signal (obstacles, interférences radio, position et mouvement 25 du terminal, etc.) rendent très complexe la fonction qui définit la distance à partir d'un niveau de RSSI. De nouveaux procédés de géolocalisation basés sur les niveaux de RSSI ont alors été développés. Ces nouveaux procédés reposent sur des techniques d'apprentissage automatique. Concrètement, il s'agit de construire, 30 pendant une phase de calibration, une base de données qui associe à des positions géographiques connues une signature radio correspondant à l'ensemble des niveaux de RSSI mesurés pour un terminal à une position connue pour un ensemble de stations de base du système. Ensuite, pendant 3 une phase de recherche, une signature radio observée pour un terminal situé à une position inconnue est comparée à l'ensemble des signatures de la base de données afin d'estimer la position du terminal à partir de la (ou des) position(s) correspondant à la (aux) signature(s) la (les) plus ressemblante(s). 5 Pour réaliser la phase de calibration, il est connu d'embarquer dans une flotte de véhicules qui parcourent la zone géographique à couvrir des appareils adaptés pour fournir avec précision la position géographique et les niveaux de RSSI pour les stations de base du système de communication en différents points (le terme anglais pour qualifier cette phase est « war- 10 driving »). Plus le nombre de points est important, meilleures seront les performances du procédé de géolocalisation en termes de précision, mais plus la phase de calibration sera longue et coûteuse. De telles méthodes de géolocalisation par apprentissage automatique à partir de signatures radio présentent cependant plusieurs inconvénients, 15 notamment dans le cas où la zone géographique à couvrir est vaste, par exemple s'il faut couvrir tout un pays, voire un continent. Un premier inconvénient concerne la précision de géolocalisation. Il arrive en effet que deux signatures radio correspondant à deux positions géographiques voisines soient significativement différentes, ou bien que deux signatures radio 20 particulièrement ressemblantes correspondent à deux positions géographiques éloignées. L'établissement d'une relation entre une signature radio et une position géographique associée est ainsi rendu difficile à cause du bruit généré par de telles situations. Un autre inconvénient concerne la complexité de l'algorithme d'apprentissage automatique utilisé. En effet, pour obtenir une 25 précision suffisante, il est nécessaire d'insérer un grand nombre d'éléments dans la base de données construite pendant la phase de calibration. De plus, si la zone à couvrir est vaste, il y a un grand nombre de stations de base à prendre en compte. Tout ceci concourt à augmenter considérablement les données en entrée de l'algorithme d'apprentissage automatique, et ainsi à 30 créer des contraintes de capacité et de temps de calcul. 4 EXPOSÉ DE L'INVENTION La présente invention a pour objectif de remédier à tout ou partie des inconvénients de l'art antérieur, notamment ceux exposés ci-avant, en proposant un procédé de géolocalisation d'un terminal d'un système de 5 communication sans fil dont les performances en termes de précision et de complexité sont améliorées en se basant sur l'idée que des terminaux du système peuvent évoluer en groupes. Les méthodes conventionnelles de géolocalisation utilisant une signature radio d'un terminal pour estimer sa position géographique se 10 concentrent en effet uniquement sur des informations liées à une relation existant entre le terminal et un réseau d'accès du système de communication, comme par exemple des niveaux de RSSI ou des temps de propagation de signaux échangés entre le terminal et le réseau d'accès. La méthode de géolocalisation proposée par la présente invention se 15 démarque de l'art antérieur dans la mesure où elle utilise en plus des informations provenant des terminaux du système qui permettent de déterminer par exemple si un terminal dont on cherche à estimer la position géographique s'est déplacé en groupe avec d'autres terminaux pendant une certaine période d'observation. Si tel est le cas, il est par exemple possible 20 d'estimer la position géographique dudit terminal en utilisant des informations relatives à d'autres terminaux du groupe. La présente invention trouve une application particulièrement avantageuse, bien que nullement limitative, dans le domaine du transport de marchandises. Il est par exemple possible d'équiper des palettes utilisées pour 25 le transport de marchandises avec des terminaux d'un système de communication sans fil afin de pouvoir les géolocaliser. De telles palettes ont en effet une certaine propension à être déplacées en groupe. Ainsi, selon un premier aspect, il est proposé par la présente invention un procédé d'estimation de la position géographique d'un terminal, dit 30 « terminal d'intérêt », parmi un panel de terminaux d'un système de communication sans fil. Ledit procédé comporte : 5 une détermination, pour chaque terminal dudit panel, d'une signature temporelle comportant des valeurs mesurées pour ledit terminal pendant une durée d'observation prédéterminée, un calcul de valeurs d'un critère de similarité entre les signatures 5 temporelles des terminaux du panel, une valeur dudit critère de similarité calculée pour deux terminaux étant représentative de la probabilité pour lesdits terminaux de s'être déplacés ensemble ou d'avoir été localisés au même endroit pendant ladite période d'observation, un partitionnement des différents terminaux du panel en différents 10 groupes en fonction des valeurs du critère de similarité, une estimation de la position géographique du terminal d'intérêt en fonction de données de géolocalisation disponibles pour au moins un autre terminal du groupe, dit « groupe d'intérêt », auquel le terminal d'intérêt appartient. 15 Le panel de terminaux peut par exemple être déterminé par une sélection des terminaux du système dont on sait a priori qu'ils sont proches du terminal d'intérêt (par exemple s'ils sont sous la couverture d'au moins une station de base du système de communication sans fil qui couvre également le terminal d'intérêt), ou qu'ils ont une certaine propension à appartenir à un 20 groupe de plusieurs terminaux avec lesquels ils présentent des similarités. Une signature temporelle d'un terminal correspond par exemple à un ensemble de valeurs représentatives des phases de mobilité du terminal pendant ladite période d'observation déterminées à l'aide d'un capteur du terminal comme par exemple un capteur de mouvement. 25 Si deux signatures temporelles pour deux terminaux ont des similarités importantes, cela signifie que les deux terminaux en question se sont probablement déplacés ensemble, ou bien éventuellement qu'ils sont restés au même endroit, pendant la période d'observation. Le partitionnement des différents terminaux du panel en différents 30 groupes permet d'améliorer la géolocalisation d'un terminal d'intérêt dont on cherche à estimer la position géographique. En effet, il est par exemple possible d'estimer la position géographique dudit terminal d'intérêt si l'on 6 cannait avec précision la position géographique d'au moins l'un des terminaux du groupe, par exemple si un des terminaux du groupe est équipé d'un récepteur GPS. Selon un autre exemple, il est possible d'estimer la position géographique dudit terminal en corrélant des positions estimées de plusieurs 5 terminaux du groupe. Dans certains cas, il est possible d'estimer la position géographique dudit terminal d'intérêt en corrélant des données de géolocalisation (par exemple des positions GPS, et/ou des positions estimées, etc.) relatives à plusieurs terminaux de plusieurs groupes différents. Il convient de noter que ces données de géolocalisation peuvent être des informations qui 10 ne correspondent pas directement à une position géographique mais qui permettent de déterminer ou d'estimer une position géographique (par exemple des temps de propagation d'un signal ou des niveaux de RSSI utilisés par des méthodes de trilatération, etc.), Dans des modes particuliers de mise en oeuvre, l'invention peut 15 comporter en outre l'une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons techniquement possibles. Dans des modes particuliers de mise en oeuvre, les valeurs mesurées pour déterminer une signature temporelle d'un terminal comportent: des valeurs mesurées par au moins un capteur dudit terminal pendant la 20 période d'observation, et/ou des signatures radio comportant des valeurs représentatives de la qualité de liens radio existant entre le terminal et une ou plusieurs stations de base du système de communication sans fil à différents instants pendant la période d'observation. 25 Par exemple, une signature temporelle peut comporter des valeurs représentatives des phases de mobilité d'un terminal pendant la période d'observation déterminées par exemple à l'aide d'un capteur de mouvement tel qu'un accéléromètre ou un gyroscope. Alternativement ou en complément, une signature temporelle peut comporter des signatures radio d'un terminal 30 déterminées à différents instants pendant la période d'observation. Si les signatures temporelles de deux terminaux se ressemblent, autrement dit si les valeurs ou les signatures radio qui les composent remplissent certaines conditions de similarité, alors les deux terminaux peuvent probablement être 7 considérés comme appartenant à un même groupe de terminaux qui se sont déplacés ensemble au cours de la période d'observation. Dans des modes particuliers de mise en oeuvre, les valeurs mesurées par au moins un capteur dudit terminal comportent : 5 des valeurs mesurées par un capteur de mouvement représentatives de phases de mobilité dudit terminal (70), et/ou des valeurs mesurées par un capteur de température, et/ou des valeurs mesurées par un capteur de pression atmosphérique, et/ou des valeurs mesurées par un capteur de luminosité, et/ou 10 des valeurs mesurées par un capteur de champ magnétique. En effet, selon le type d'application considérée, différentes données environnementales peuvent être utilisées pour déterminer si des terminaux ont évolué ensemble pendant une période d'observation. Par exemple si les conditions environnementales observées pour deux terminaux distincts sont 15 similaires, cela signifie que lesdits terminaux étaient probablement « co- localisés » pendant la période d'observation. Dans des modes particuliers de mise en oeuvre, une signature temporelle d'un terminal comporte à la fois des valeurs représentatives des phases de mobilité et des signatures radio, et le critère de similarité est 20 déterminé en fonction d'un critère de similarité relatif aux phases de mobilité et d'un critère de similarité relatif aux signatures radio. Dans des modes particuliers de mise en oeuvre, un critère de similarité est déterminé en fonction d'un facteur de pondération affecté à une valeur mesurée représentatif du temps écoulé entre l'instant où ladite valeur est 25 mesurée et l'instant d'estimation de la position géographique du terminal d'intérêt. De telles dispositions permettent d'accorder plus d'importance aux données récentes d'une signature temporelle. En effet, si par exemple deux terminaux se sont déplacés ensemble pendant seulement une partie de la 30 période d'observation, il convient d'accorder plus d'importance aux données récentes afin de déterminer si à la fin de la période d'observation, c'est-à-dire au moment où l'estimation de la position géographique du terminal d'intérêt est 8 faite, les terminaux en question sont ensemble ou non, c'est-à-dire si oui ou non il faut considérer qu'ils appartiennent à un même groupe. Dans des modes particuliers de mise en oeuvre, le partitionnement des différents terminaux du panel est fait de sorte que chaque terminal appartient 5 au plus à un groupe. Un tel partitionnement est par exemple mis en oeuvre par des méthodes de regroupement d'éléments connues sous le terme anglais « hard clustering », pour lesquelles chaque élément appartient à un unique groupe, ou bien il n'appartient à aucun groupe (élément isolé). 10 Dans des modes particuliers de mise en oeuvre, le partitionnement des différents terminaux du panel est fait de sorte que chaque terminal appartient à chaque groupe selon une valeur de probabilité d'appartenance, le groupe d'intérêt correspondant au groupe pour lequel la valeur de probabilité d'appartenance du terminal d'intérêt est la plus grande. 15 Un tel partitionnement est par exemple mise en oeuvre par des méthodes de regroupement d'éléments connues sous le terme anglais « soft clustering », pour lesquelles chaque élément appartient à plusieurs groupes avec un certain degré défini par des méthodes probabilistes. Dans des modes particuliers de mise en oeuvre, la position 20 géographique du terminal d'intérêt est estimée en fonction de données de géolocalisation disponibles pour différents terminaux appartenant à différents groupes et en fonction des valeurs de probabilité d'appartenance du terminal d'intérêt aux différents groupes. Dans des modes particuliers de mise en oeuvre, lesdites données de 25 géolocalisation sont fournies par un terminal sous la forme d'une position géographique obtenue par un système de positionnement dudit terminal. Si certains terminaux du système sont équipés d'un récepteur GPS, il peut être possible de déterminer la position géographique du terminal d'intérêt à l'aide d'une position GPS fourni par l'un des terminaux du groupe d'intérêt 30 auquel le terminal d'intérêt appartient, ou bien à l'aide de plusieurs positions GPS fournies par plusieurs terminaux du groupe d'intérêt, voire à l'aide de plusieurs positions GPS fournies par plusieurs terminaux de plusieurs groupes différents. 9 Dans des modes particuliers de mise en oeuvre, les données de géolocalisation sont disponibles pour un terminal sous la forme d'une signature radio comportant des valeurs représentatives de la qualité de liens radio existant entre ledit terminal et une ou plusieurs stations de base du système de 5 communication sans fil. Dans des modes particuliers de mise en oeuvre, une position géographique d'un terminal du panel est estimée en fonction de la signature radio dudit terminal et d'un ensemble de référence comportant des signatures radio associées à des positions géographiques connues. La position 10 géographique du terminal d'intérêt est alors estimée en fonction de la position géographique estimée d'au moins un autre terminal du groupe d'intérêt. Cet ensemble de référence est par exemple une base de données mémorisant des éléments de référence utilisée par un algorithme d'apprentissage automatique. Chaque élément de référence est alors un 15 couple d'information constitué d'une part par une signature radio, et d'autre part par une position géographique connue associée à cette signature radio. L'algorithme d'apprentissage permet d'estimer la position d'un terminal d'intérêt à partir d'une signature radio dudit terminal d'intérêt, et à l'aide de la base de données de référence. Pour réduire la complexité du procédé de 20 géolocalisation, il peut alors suffire d'estimer la position géographique d'un seul terminal d'un groupe donné et considérer que tous les terminaux du groupe ont la même position géographique. Selon un autre exemple, pour améliorer la précision du procédé de géolocalisation, il peut être avantageux de corréler les positions géographiques 25 estimées de plusieurs terminaux du groupe d'intérêt afin d'affiner la position géographique de chaque terminal du groupe, et notamment celle du terminal d'intérêt. Dans des modes particuliers de mise en oeuvre, une signature radio virtuelle est déterminée pour le groupe d'intérêt en fonction des signatures 30 radio des terminaux appartenant audit groupe d'intérêt, et la position du terminal d'intérêt est estimée en fonction de ladite signature radio virtuelle. Il s'agit par exemple dans ce cas de corréler les signatures radio de plusieurs terminaux du groupe d'intérêt pour déterminer une signature radio 10 virtuelle utilisée pour estimer une position géographique qui est attribuée à tous les terminaux du groupe, et notamment au terminal d'intérêt. Dans des modes particuliers de mise en oeuvre, ledit panel de terminaux est déterminé par une étape de sélection parmi les terminaux du 5 système de communication sans fil de sorte qu'un terminal du panel est couvert par au moins une station de base du système de communication sans fil couvrant également le terminal d'intérêt. Selon un deuxième aspect, la présente invention concerne un système de communication sans fil comportant des terminaux, et un réseau d'accès 10 comportant des stations de base et un serveur relié auxdites stations de base. Ledit système met en oeuvre un procédé de géolocalisation d'un terminal d'intérêt, parmi un panel de terminaux, selon l'un quelconque des modes de mises en oeuvre précédents. Ledit serveur est configuré pour : collecter des valeurs mesurées pour un terminal dudit panel pendant 15 une durée d'observation prédéterminée, déterminer une signature temporelle pour ledit terminal à partir des valeurs mesurées, calculer des valeurs d'un critère de similarité entre les signatures temporelles des terminaux du panel, une valeur dudit critère de similarité 20 calculée pour deux terminaux étant représentative de la probabilité pour lesdits terminaux de s'être déplacés ensemble ou d'avoir été localisés au même endroit pendant ladite période d'observation, partitionner les différents terminaux du panel en différents groupes en fonction des valeurs du critère de similarité, 25 estimer la position géographique du terminal d'intérêt en fonction de données de géolocalisation disponibles pour au moins un autre terminal du groupe, dit « groupe d'intérêt », auquel le terminal d'intérêt appartient. 30 PRÉSENTATION DES FIGURES L'invention sera mieux comprise à la lecture de la description suivante, donnée à titre d'exemple nullement limitatif, et faite en se référant aux figures 1 11 à 7 qui représentent : - Figure 1 : une représentation schématique d'un système de communication sans fil, - Figure 2: les étapes principales d'un procédé conventionnel de 5 géolocalisation, - Figure 3: les étapes principales d'un procédé de géolocalisation selon l'invention, - Figure 4: une représentation schématique de phases de mobilité déterminées par un terminal, 10 - Figure 5: une représentation schématique de la pression atmosphérique mesurée par un terminal, - Figure 6: une représentation schématique de la température mesurée par un terminal, - Figure 7: une représentation schématique de signatures radio 15 déterminées pour un terminal, Dans ces figures, des références identiques d'une figure à une autre désignent des éléments identiques ou analogues. Pour des raisons de clarté, les éléments représentés ne sont pas nécessairement à une même échelle, sauf mention contraire. 20 DESCRIPTION DÉTAILLÉE DE MODES DE RÉALISATION La description ci-après décrit à titre d'exemple et de manière nullement limitative plusieurs modes de mise en oeuvre de l'invention. 1. Le contexte 25 Dans la suite de la description, on se place de manière non limitative dans le cas où des terminaux de type objets connectés d'un système de communication sans fil sont fixés à des palettes de marchandises dans le but de pouvoir les géolocaliser. De telles palettes sont généralement déplacées en groupe, selon la destination des marchandises qu'elles comportent. 30 La figure 1 représente schématiquement un système 60 de communication sans fil, comportant des terminaux 70 et un réseau d'accès 80 comportant par exemple une pluralité de stations de base 81 reliées à un serveur 82. 12 Les terminaux 70 et les stations de base 81 échangent des données sous la forme de signaux radioélectriques. Par « signal radioélectrique », on entend une onde électromagnétique se propageant en espace libre, dont les fréquences sont comprises dans le spectre traditionnel des ondes 5 radioélectriques (quelques hertz à plusieurs centaines de gigahertz). Les échanges de données entre un terminal 70 et les stations de base 81 sont par exemple bidirectionnels. En d'autres termes, un terminal 70 est adapté à émettre des données sur un lien montant à destination des stations de base 81, et à recevoir des données sur un lien descendant depuis lesdites 10 stations de base 81 vers ledit terminal 70. Rien n'exclut cependant, suivant d'autres exemples, d'avoir des échanges unidirectionnels, uniquement sur le lien montant. De nombreuses applications du type loT consistent en de la collecte de données émises par des terminaux 70 et s'accommodent parfaitement d'échanges uniquement sur 15 le lien montant entre chaque terminal 70 et les stations de base 81. Suivant un exemple non limitatif, le système 60 de communication sans fil est à bande ultra étroite. Par « bande ultra étroite » (« Ultra Narrow Band » ou UNB dans la littérature anglo-saxonne), on entend que le spectre fréquentiel instantané des signaux émis par les terminaux 70 est de largeur 20 fréquentielle inférieure à deux kilohertz, voire inférieure à un kilohertz. De telles dispositions sont particulièrement avantageuses en ce que l'émission de tels signaux peut être réalisée avec une consommation électrique très réduite, particulièrement adaptée pour des applications du type loT. Chaque station de base 81 est adaptée à recevoir des messages émis 25 par des terminaux 70 qui se trouvent à sa portée. Chaque message ainsi reçu est par exemple transmis au serveur 82 du réseau d'accès 80, éventuellement accompagné d'autres informations comme un identifiant de la station de base 81 qui l'a reçu, une valeur représentative de la qualité du signal radio transportant le message, la fréquence centrale sur laquelle le message a été 30 reçu, une date à laquelle le message a été reçu, etc. Le serveur 82 traite par exemple l'ensemble des messages reçus des différentes stations de base 81. Pour estimer la position géographique d'un terminal 70, dit « terminal 70i d'intérêt », il est connu par exemple d'utiliser un procédé de géolocalisation 13 basé sur des algorithmes d'apprentissage automatique cherchant à établir une relation entre une signature radio d'un terminal et la position géographique dudit terminal. Le serveur 82 peut notamment être utilisé pour la mise en oeuvre d'un tel procédé de géolocalisation. 5 Dans la suite de la description, on se place à titre d'exemple nullement limitatif dans le cas où l'invention améliore un tel procédé de géolocalisation existant basé sur un algorithme d'apprentissage automatique qui vise à établir une relation entre une signature radio d'un terminal et la position géographique dudit terminal. 10 Il convient de noter cependant qu'il existe d'autres méthodes pour estimer la position d'un terminal 70i d'intérêt, et en particulier des méthodes qui ne sont pas nécessairement basées sur des algorithmes d'apprentissage automatique, comme par exemple des méthodes de trilatération qui permettent de déterminer une position géographique à partir de calculs de différence de 15 temps de propagation d'un signal (« Time Difference Of Arrivai» ou TDOA dans la littérature anglo-saxonne). Le procédé selon l'invention est également applicable à de telles méthodes. Par « position géographique », on entend dans la suite de la description un ensemble de deux coordonnées correspondant respectivement 20 à la latitude et la longitude. Il convient de noter que des variantes pourraient être envisagées pour définir une position géographique. Par exemple, une troisième coordonnée correspondant à l'altitude par rapport au niveau moyen de la mer pourrait également être prise en compte. La signature radio d'un terminal 70 est par exemple déterminée lors de 25 l'émission d'un message dudit terminal 70 à destination du réseau d'accès 80. Par exemple, la signature radio comporte un ensemble de valeurs correspondant à des niveaux de puissance moyenne reçue (RSSI) mesurés par les différentes stations de base (par exemple en décibel) à la réception dudit message. Pour les stations de base 81 qui n'ont pas reçu le message, 30 par exemple parce que le terminal 70 est trop éloigné et ne se trouve pas sous leur couverture radio, une valeur par défaut de -160 dB est par exemple utilisée. Il convient de noter que d'autres méthodes peuvent être envisagées 14 pour déterminer une signature radio. Par exemple, d'autres valeurs représentatives de la qualité des liens radio existant entre un terminal 70 et les stations de base 81 lors de l'émission d'un message pourraient être utilisées pour déterminer une signature radio, comme par exemple un rapport signal sur 5 bruit du signal (SNR ou « Signal on Noise Ratio » dans la littérature anglo- saxonne) ou bien un indicateur de qualité du canal de communication (COI ou « Channel Quality Indicator » dans la littérature anglo-saxonne). Selon un autre exemple, une signature radio peut comporter des temps de propagation d'un signal transportant un message échangé entre un terminal 70 et une station de 10 base 81. Le choix d'une méthode particulière pour déterminer une signature radio ne constitue qu'une variante de mise en oeuvre de l'invention. Pour estimer la position géographique d'un terminal 70i d'intérêt, le procédé de géolocalisation utilise un ensemble d'éléments de référence. Chaque élément de référence comporte une signature radio associée à une 15 position géographique. Pour construire l'ensemble d'éléments de référence, il est par exemple possible d'utiliser certains terminaux 70 du système 60 de communication sans fil, dits « terminaux 70c de calibration », qui possèdent un système de positionnement, par exemple un récepteur GPS, qui permet d'obtenir avec précision la position géographique du terminal 70c. 20 De manière conventionnelle, cet ensemble de référence est par exemple une base de données mémorisée au niveau du serveur 82 et utilisée par un algorithme d'apprentissage automatique. Chaque élément de référence stocké dans cette base de données est un couple d'information constitué d'une part par une signature radio, et d'autre part par une position géographique 25 connue associée à cette signature radio. La signature radio et la position géographique connue d'un élément de référence est obtenu à partir d'un terminal 70c de calibration. L'algorithme d'apprentissage permet alors d'estimer la position d'un terminal 70i d'intérêt à partir d'une signature radio dudit terminal 70i d'intérêt, et à l'aide de la base de données. 30 La figure 2 représente les principales étapes d'un tel procédé 10 de géolocalisation. Une première phase 20 de calibration a pour objectif de réaliser une sorte de cartographie radio de la zone géographique à couvrir. Cette phase 15 comporte : - une étape de détermination 22 de signatures radio associées respectivement à des positions géographiques, - une étape de stockage 24 de l'élément de référence obtenu (association 5 d'une signature radio et d'une position géographique) dans un ensemble 26 d'éléments de référence. L'étape de détermination 22 d'une signature radio consiste par exemple en l'émission par un terminal 70c de calibration d'un message contenant la position géographique courante à destination du réseau d'accès 10 80 du système 60 de communication sans fil. Il est à noter que ce message peut être émis de manière indépendante au procédé 10 de géolocalisation. Par exemple, il peut s'agir d'un message de télérelève d'informations conventionnel incluant la position géographique courante et n'ayant pas pour objectif premier de participer à la phase 20 de 15 calibration du procédé 10 de géolocalisation. Les stations de base 81 du réseau d'accès 80 qui ont reçu le signal contenant ledit message font alors une mesure du niveau RSSI auquel le message a été reçu et transmettent le message et des informations supplémentaires (niveau de RSSI, date de réception du message, identifiant de 20 la station de base, etc.) au serveur 82. Le serveur 82 collecte les niveaux de RSSI mesurés par les différentes stations de base 81 et les inclut dans la signature radio ainsi déterminée pour le terminal 70c de calibration. Il convient de noter que plusieurs méthodes sont possibles pour le 25 choix des stations de base à prendre en compte pour la formation d'une signature radio. Selon un premier exemple, toutes les stations de base du système 60 de communication sans fil sont considérées. Selon un autre exemple, il est possible de limiter le nombre de stations de base à prendre en compte, par exemple en fonction d'une zone géographique limitée à laquelle on 30 s'intéresse. Le serveur 82 procède aussi à l'extraction de la position géographique courante contenue dans le message émis par le terminal 70c de calibration et transmis au serveur 82 par les stations de base 81 qui l'ont reçu. 16 Finalement, le couple d'information formé par la position géographique du terminal 70c de calibration et la signature radio associée est ajouté à l'ensemble 26 des éléments de référence à travers l'étape de stockage 24. Le serveur 82 comporte par exemple un ensemble de moyens matériels et 5 logiciels pour mémoriser les éléments de référence dans une base de données. Dans des modes particuliers de mise en oeuvre, les étapes de la phase 20 de calibration décrites ci-avant sont répétées pendant une durée prédéfinie, par exemple pendant plusieurs jours, voire plusieurs semaines ou 10 plusieurs mois, pour des terminaux 70c de calibration du système 60 de communication sans fil, ceci afin d'obtenir un ensemble 26 d'éléments de référence contenant un nombre conséquent d'informations, c'est-à-dire une cartographie précise de la zone à couvrir. Dans une variante, les étapes de la phase 20 de calibration sont répétées jusqu'à obtenir un nombre prédéfini 15 d'éléments de référence dans l'ensemble 26. Dans un mode préféré de mise en oeuvre, l'ensemble 26 des données de référence est continuellement enrichi par de nouveaux éléments de référence provenant des terminaux 70c de calibration. Ensuite, une phase 40 de recherche comporte : 20 - une étape de détermination 42 d'une signature radio pour un terminal 70i d'intérêt situé à une position géographique à estimer, - une étape d'estimation 44 de la position géographique dudit terminal 70i d'intérêt à partir de la signature radio déterminée et de l'ensemble 26 d'éléments de référence acquis pendant la phase 20 de calibration. 25 L'étape de détermination 42 d'une signature radio consiste par exemple en l'émission d'un message par un terminal 70i d'intérêt à destination du réseau d'accès 80 du système 60 de communication sans fil. Ce message peut être un message quelconque émis de manière indépendante au procédé 10 de géolocalisation du terminal. Par exemple, ce message peut être émis à 30 des fins de télérelève d'informations qui n'ont rien à voir avec la géolocalisation du terminal 70. Alternativement, ce message peut être émis intentionnellement dans le but de géolocaliser le terminal 70. Dans tous les cas, le contenu du message n'a pas nécessairement d'importance dans la phase 40 de recherche. 17 Similairement à ce qui est fait pour la phase 20 de calibration, les stations de base 81 du réseau d'accès 80 qui ont reçu ledit message procèdent à une mesure du niveau de RSSI du signal transportant le message. Le serveur 82 collecte alors les niveaux de RSSI mesurés par les 5 différentes stations de base 81 et les inclut dans la signature radio ainsi déterminée pour le terminal 70i d'intérêt. L'estimation 44 de la position géographique du terminal 70i d'intérêt ayant émis le message est alors effectuée, par exemple à l'aide d'un algorithme d'apprentissage automatique de régression, à partir d'une part de 10 l'ensemble 26 d'éléments de référence, et d'autre part de la signature radio déterminée pour le terminal 70i d'intérêt dont la position géographique doit être estimée. Le procédé 10 de géolocalisation décrit ci-avant présente cependant des performances médiocres en termes de précision de géolocalisation et il est 15 particulièrement complexe (c'est-à-dire qu'il nécessite des capacités et des temps de calcul importants), notamment dans le cas où la zone géographique à couvrir est vaste, par exemple s'il faut couvrir tout un pays, voire un continent. 2. Le principe de l'invention 20 Pour améliorer la précision et/ou pour simplifier la complexité du procédé 10 de géolocalisation, la présente invention se base sur l'idée que les terminaux 70 du système considéré peuvent se déplacer en groupes. L'invention propose ainsi un procédé 10 de géolocalisation d'un terminal 70i d'intérêt qui tire profit des informations fournies par d'autres terminaux 70 qui 25 se sont déplacés conjointement avec le terminal 70i d'intérêt pendant une certaine période. La figure 3 représente les principales étapes d'un mode particulier de mise en oeuvre d'un procédé 10 de géolocalisation selon l'invention. Il est à noter que l'ordre des différentes étapes représentées sur la figure 3 est donné 30 à titre illustratif uniquement, à des fins de clarté de la figure, et ne saurait être considéré comme limitatif de l'invention. Pour le procédé 10 de géolocalisation considéré et présentement décrit à titre d'exemple, la phase 20 de calibration est similaire à celle décrite 18 ci-avant en référence à la figure 2. La phase 40 de recherche comporte en revanche des étapes supplémentaires qui visent à former des groupes de terminaux 70, chaque groupe comportant des terminaux 70 pour lesquels on considère qu'ils se sont déplacés ensemble pendant une certaine période 5 d'observation courant jusqu'à l'instant d'estimation de la position du terminal 70i d'intérêt. La phase 40 de recherche comporte une étape de détermination 42 d'une signature radio pour un terminal 70i d'intérêt dont on souhaite estimer la position géographique. Cette étape est similaire à celle décrite ci-avant en 10 référence à la figure 2. La phase 40 de recherche comporte aussi une étape de détermination 52, pour chaque terminal 70 d'un panel de terminaux 70 du système 60 de communication sans fil, d'une signature temporelle. La signature temporelle d'un terminal 70 comporte des valeurs mesurées pour ledit terminal 70 pendant 15 une période d'observation passée dont la durée est prédéterminée en fonction de l'application considérée. Ces valeurs correspondent à des données liées à l'environnement dans lequel a évolué le terminal 70 ou bien liées au comportement qu'a eu le terminal 70 pendant ladite période d'observation. La phase 40 de recherche comporte aussi une étape de calcul 54 de 20 valeurs d'un critère de similarité entre les signatures temporelles des terminaux 70 du panel. Si deux signatures temporelles pour deux terminaux 70 ont des similarités importantes, cela signifie que les deux terminaux 70 en question se sont probablement déplacés ensemble, ou bien éventuellement qu'ils sont restés au même endroit, pendant la période d'observation. 25 La phase 40 de recherche comporte aussi une étape de partitionnement 56, des différents terminaux 70 du panel en différents groupes en fonction du critère de similarité. Enfin, la phase 40 de recherche comporte une étape d'estimation 44 de la position géographique du terminal 70i d'intérêt en fonction de données de 30 géolocalisation disponibles pour au moins un autre terminal 70 d'un groupe, dit « groupe d'intérêt », auquel le terminal 70i d'intérêt appartient. Les différentes étapes de la phase 40 de recherche sont décrites ci-après en détails. 19 3. Détermination d'une signature temporelle Différentes méthodes sont envisageables pour l'étape de détermination 52 d'une signature temporelle d'un terminal 70 du panel considéré. 5 Dans un mode particulier de mise en oeuvre, la signature temporelle d'un terminal 70 comporte des valeurs représentatives des phases de mobilité dudit terminal 70 pendant la période d'observation. Les phases de mobilité sont par exemple déterminées à l'aide de mesures réalisées par un capteur de mouvement du terminal 70 tel qu'un accéléromètre. 10 On entend par « phase de mobilité » un intervalle de temps pendant lequel le terminal 70 est considéré comme étant en train de se déplacer. Deux phases de mobilité successives sont séparées par une phase d'immobilité, qui correspond à un intervalle de temps pendant lequel le terminal est considéré comme étant immobile. 15 La figure 4 représente schématiquement des phases de mobilité d'un terminal 70 du système 60 de communication sans fil pendant une période d'observation Tobs. La période d'observation Tobs précède un instant to correspondant par exemple à l'instant où la position géographique d'un terminal 70i d'intérêt est estimée par le procédé 10 de géolocalisation. L'instant 20 to correspond généralement à un instant de réception d'un message émis par le terminal 70i d'intérêt, message à partir duquel une signature radio du terminal 70i d'intérêt est déterminée à l'étape 42 de la phase 40 de recherche. La courbe 90 de la figure 4 représente les phases de mobilité du terminal 70. Lorsque la courbe 90 prend une valeur constante non nulle (par 25 exemple la valeur 1), cela signifie que le terminal 70 est en mouvement (il est dans une phase de mobilité). Lorsque la courbe prend une valeur nulle (valeur 0), cela signifie que le terminal 70 est immobile (il est dans une phase d'immobilité). La courbe 90 peut être construite, par exemple par le serveur 82, à 30 l'aide d'informations remontées par le terminal 70 dans des messages émis à destination du réseau d'accès 80. A cet effet, le terminal 70 comporte un capteur de mouvement, par exemple un accéléromètre, et un ensemble de moyens conventionnels 20 configurés de façon logicielle (produit programme d'ordinateur) et ou matérielle (un ou des circuits logiques programmable de type FPGA et/ou un ou des circuits intégrés spécialisés de type ASIC, etc.) pour déterminer des instants de début et des instants de fin respectifs de phases de mobilité du terminal 70. 5 Par exemple, une mesure d'accélération est réalisée de manière récurrente par le capteur (par exemple avec une période comprise entre quelques dizaines de millisecondes et quelques centaines de millisecondes) et si la valeur mesurée est supérieure à un seuil prédéterminé de détection de mouvement, un instant de détection courant est mémorisé et comparé à l'instant de détection 10 précédent. Si les deux instants de détection sont séparés d'une durée inférieure à un seuil prédéterminé de détection d'une nouvelle phase de mobilité, alors il est considéré que les deux instants de détection appartiennent à la même phase de mobilité. Si en revanche les deux instants de détection sont séparés d'une durée supérieure au seuil de détection d'une nouvelle 15 phase de mobilité, alors il est considéré que l'instant de détection précédent correspond à l'instant de fin d'une phase de mobilité, et l'instant de détection courant correspond à l'instant de début d'une nouvelle phase de mobilité. Les instants de début et les instants de fin des phases de mobilité du terminal 70 sont par exemple mémorisés puis envoyés au serveur 82 dans des 20 messages émis de manière récurrente par le terminal 70 à destination du réseau d'accès 80. De tels messages peuvent par exemple être émis périodiquement, ou bien à chaque détection d'une nouvelle phase de mobilité. A partir de la courbe 90, le serveur 82 peut déterminer une signature temporelle pour le terminal 70 pour la période d'observation Tobs, par exemple 25 en échantillonnant un ensemble de valeurs (M1, M2, M3, ..., Md. La période d'échantillonnage est choisie de manière appropriée en fonction des durées typiques des phases de mobilité et des phases d'immobilité du terminal 70 pour l'application considérée. Pour une période d'observation fixe, plus le nombre L de valeurs est grand, plus on a d'information sur les mouvements effectués par 30 le terminal 70, et plus la signature temporelle est précise. En utilisant le même échantillonnage sur la même période d'observation pour tous les terminaux 70 du panel considéré, il est alors possible de comparer les signatures temporelles des différents terminaux 70 21 deux à deux. Si deux terminaux 70 présentent des signatures temporelles similaires cela signifie qu'ils ont eu des phases de mobilité similaires pendant la période d'observation, et en conséquence cela signifie qu'ils se sont probablement déplacés ensemble. 5 Dans un mode particulier de mise en oeuvre, la signature temporelle d'un terminal 70 comporte des valeurs représentatives de la pression atmosphérique à laquelle est soumis ledit terminal 70. Ces valeurs sont par exemple mesurées par un capteur de pression atmosphérique du terminal 70. La figure 5 représente une courbe 91 de l'évolution de la pression 10 atmosphérique à laquelle est soumis un terminal 70 pendant une période d'observation Tobs précédant un instant to d'estimation de la position géographique d'un terminal 70i d'intérêt. La courbe 91 est par exemple construite par le serveur 82 à l'aide de mesures de pression atmosphérique réalisées par le terminal 70 et envoyées 15 au serveur 82 dans des messages émis à destination du réseau d'accès 80 de manière récurrente. A partir de la courbe 91, le serveur 82 peut déterminer une signature temporelle pour le terminal 70 pour la période d'observation Tobs, par exemple en échantillonnant un ensemble de valeurs (Pl, P2, P3, ---, PK). Là encore, la 20 période d'échantillonnage (et donc le nombre K de valeurs échantillonnées pendant la période d'observation Tobs) est choisie de manière appropriée en fonction de l'application considérée. Par exemple, si des palettes de marchandises équipées avec des terminaux 70 sont transportées par avion, il est intéressant de distinguer des 25 phases de transport aérien correspondant aux périodes où la pression atmosphérique est basse (en raison de l'altitude élevée de l'avion), et des phases au sol correspondant aux périodes où la pression atmosphérique est plus élevée. La période d'échantillonnage peut alors être choisie en fonction des durées typiques des phases de transport aérien et des phases au sol. 30 En utilisant le même échantillonnage sur la même période d'observation pour tous les terminaux 70 du panel considéré, il est alors possible de comparer les signatures temporelles des différents terminaux 70 deux à deux. Si deux terminaux 70 présentent des signatures temporelles 22 similaires cela signifie qu'ils ont eu des phases de transport aérien similaires pendant la période d'observation, et en conséquence cela signifie qu'ils se sont probablement déplacés ensemble. Il convient de noter que d'autres méthodes sont possibles pour 5 déterminer une signature temporelle à partir de la courbe 91. Par exemple, il est possible de découper la période d'observation Tobs en différents intervalles de temps, et de déterminer une valeur pour chacun des intervalles, comme par exemple la valeur moyenne prise par la courbe 91 sur un intervalle de temps, ou bien la valeur minimale ou la valeur maximale. 10 Dans un mode particulier de mise en oeuvre, la signature temporelle d'un terminal 70 comporte des valeurs représentatives de la température à laquelle est soumise ledit terminal 70. Ces valeurs sont par exemple mesurées par un capteur de température du terminal 70. La figure 6 représente une courbe 92 de l'évolution de la température 15 à laquelle est soumise un terminal 70 pendant une période d'observation Tobs précédant un instant to d'estimation de la position géographique d'un terminal 70i d'intérêt. La courbe 92 est par exemple construite par le serveur 82 à l'aide de mesures de température réalisées par le terminal 70 et envoyées au serveur 20 82 dans des messages émis à destination du réseau d'accès 80. A partir de la courbe 92, le serveur 82 peut déterminer une signature temporelle pour le terminal 70 pour la période d'observation Tobs, par exemple en échantillonnant un ensemble de valeurs (T1, T2, T3, ..., Tm). Là encore, la période d'échantillonnage est choisie de manière appropriée en fonction de 25 l'application considérée. Par exemple, si des palettes de marchandises équipées avec des terminaux 70 sont transportées par camion frigorifique, il est intéressant de distinguer des phases où la température est basse (palettes localisées dans une chambre froide ou dans un camion frigorifique), et des phases où la 30 température est plus élevée (palettes en train d'être chargées ou déchargées d'un camion frigorifique). La période d'échantillonnage est alors par exemple choisie en fonction des durées typiques des phases où la température est élevée. 23 En utilisant le même échantillonnage sur la même période d'observation pour tous les terminaux 70 du panel considéré, il est alors possible de comparer les signatures temporelles des différents terminaux deux à deux. Si deux terminaux présentent des signatures temporelles similaires 5 cela signifie qu'ils ont été soumis à des changements de température similaires au cours de la période d'observation, et en conséquence cela signifie qu'ils se sont probablement déplacés ensemble. Dans un mode particulier de mise en oeuvre, la signature temporelle d'un terminal 70 comporte des signatures radio déterminées pour ledit terminal 10 70 à différents instants pendant la période d'observation. La figure 7 représente schématiquement des signatures radio SA, SB, Sc, ..., SL déterminées pendant une durée d'observation Tobs précédant un instant to d'estimation de la position géographique d'un terminal 70i d'intérêt. Chaque signature radio est par exemple déterminée à partir de mesures RSSI 15 effectuées par des stations de base 81 du système 60 de communication sans fil lors à la réception de messages émis par le terminal 70 à différents instants pendant la période d'observation Tobs. Il est par exemple possible de déterminer une signature temporelle pour le terminal 70 en découpant la période d'observation Tobs en N intervalles 20 de temps de durée identique, et en définissant une signature radio Si, S2, SN pour chacun desdits intervalles de temps en fonction des signatures radio SA, SB, ..., SL appartenant audit intervalle de temps. Chaque composante d'une signature radio Si, S2, ..., SN est par exemple définie comme la valeur moyenne (ou maximale, minimale, etc.) des valeurs RSSI des composantes 25 correspondantes des signatures radio SA, SB, ..., SL appartenant à l'intervalle de temps correspondant. Par exemple, sur la figure 7, la signature radio S3 correspondant au troisième intervalle de temps de la période d'observation Tobs comporte un ensemble de valeurs (R5513,1, R5513,2, R55/3,Q), et chaque composante R5513,1 est égale à la moyenne des composantes 30 correspondantes RSS/D,i, RSS/E,i, et RSS/Fi des signatures radio SD, SE et SF respectivement. La signature temporelle pour le terminal 70 pendant la période d'observation Tobs est alors l'ensemble des valeurs (Si, S2, S3, ..., SN). A chaque signature radio Si, iE {1..N}, est par exemple associé un temps de 24 mesure ti correspondant à l'instant situé au milieu de l'intervalle de temps utilisé pour déterminer Si. En procédant de même pour tous les terminaux 70 du panel considéré, il est alors possible de comparer les signatures temporelles des 5 différents terminaux deux à deux. Si deux terminaux présentent des signatures temporelles similaires cela signifie qu'ils ont été confrontés à des conditions radio similaires au cours de la période d'observation, et en conséquence cela signifie qu'ils se sont probablement déplacés ensemble. Il convient de noter qu'au lieu d'utiliser des intervalles de temps 10 distincts pour définir les signatures radio S1, S2, ..., SN, il est possible d'utiliser un mécanisme de fenêtre temporelle glissante pour lequel une signature radio SA, SB, ..., SL est impliquée pour plusieurs positions de ladite fenêtre glissante. Chaque signature radio S1, S2, ..., SN de la signature temporelle correspond alors à une convolution entre la position de la fenêtre glissante et les 15 signatures radio SA, SB, ..., SL. Il convient également de noter que d'autres mesures environnementales pourraient être utilisées pour définir la signature temporelle d'un terminal 70, comme par exemple un niveau de luminosité, une intensité de champ magnétique, etc. Le choix d'un type particulier de mesures pour définir 20 une signature temporelle ne constitue qu'une variante de mise en oeuvre de l'invention. Dans des modes particuliers de mise en oeuvre, la signature temporelle d'un terminal 70 est une combinaison de plusieurs ensembles de valeurs correspondant à des mesures environnementales différentes. Par 25 exemple, la signature temporelle d'un terminal 70 comporte à la fois les valeurs (M1, M2, M3, «.», ML) décrites en référence à la figure 4 et les valeurs (Si, S2, S3, ..., SN) décrites en référence la figure 7. 4. Détermination d'un critère de similarité Différentes méthodes sont envisageables pour l'étape de calcul 54 de 30 valeurs d'un critère de similarité entre les signatures temporelles des terminaux 70 du panel considéré. Considérons le cas où la signature temporelle d'un terminal 70 comporte les valeurs (M1, M2, M3, « ML) décrites en référence à la figure 4 25 représentatives des phases de mobilité dudit terminal 70 pendant la période d'observation Tobs. Pour rappel, chaque valeur Mi vaut respectivement 1 ou 0 selon que le terminal 70 est ou non dans une phase de mobilité à un instant ti. Si l'on considère deux terminaux 70 différents, notés A et B, et 5 présentant respectivement les signatures temporelles (M1, M2, ..., Md et (Mi, «.», MO, il est par exemple possible de définir un critère de similarité, notée SIMm(A,B), entre les deux signatures temporelles par les formules suivantes : SIMm(A, B) = 1W - D-1(mi, mi') (1) 10 = e-(to-ti) (2) D -1 (Mi, = C - Mi - Mi' (3) 15 Dans la formule (3), C est une constante définie par exemple en fonction du nombre L et d'un ordre de grandeur souhaité pour SIMm(A,B) (qui prend dans ce cas ses valeurs entre 0 et L x C). Dans la formule (2), ti correspond à l'instant associé aux mesures Mi et Mi'. Le terme West ainsi un facteur de pondération qui permet d'accorder plus 20 d'importance aux mesures récentes. Ceci est avantageux par exemple si deux terminaux se sont déplacés ensemble pendant seulement une partie de la période d'observation. Il convient en effet dans ce cas d'accorder plus d'importance aux données récentes afin de déterminer si à la fin de la période d'observation, c'est-à-dire à l'instant to où l'estimation de la position 25 géographique du terminal 70i d'intérêt est faite, les terminaux en question sont ensemble ou non, c'est-à-dire si oui ou non il faut considérer qu'ils appartiennent à un même groupe. Ainsi, une valeur élevée pour SIMm(A,B) signifie que les signatures temporelles des terminaux A et B sur la période d'observation Tobs sont 30 semblables, ce qui signifie qu'il y a une probabilité forte que les terminaux A et B se soient déplacés ensemble pendant la période d'observation, et plus particulièrement vers la fin de la période d'observation. Considérons maintenant le cas, tel que décrit en référence à la figure 6, où la signature temporelle d'un terminal 70 comporte des valeurs de 26 température, il est par exemple possible de définir un critère de similarité, notée S/MT(A,B), entre les deux signatures temporelles de deux terminaux A et B en utilisant les formules suivantes : simT (A, B) = EnWi - D-1 (Ti, Ti') = e-(t°-ti) 1 D -1 (Ti, = max (I Ti-Till,s at) 10 Le terme sat dans la formule (6) est une constante de faible valeur (par exemple égale à 0,5°C pour l'exemple considéré) destinée à éviter une division par 0 dans la cas où la même température est mesurée pour les deux terminaux A et B à l'instant ti. 15 Considérons maintenant le cas, tel que décrit à la figure 7, où la signature temporelle d'un terminal 70 comporte des signatures radio. Considérons par exemple deux terminaux 70 différents, notés A et B, et présentant respectivement les signatures temporelles (S1, S2, SN) et (S1', S2', ..., SN). Chaque signature radio Si déterminée à un instant ti correspond à 20 un ensemble de niveaux RSSI (RSS/i,i, RSS/1,2, RSS/) mesurés par un nombre Q de stations de base 81. Il est par exemple possible de définir un critère de similarité, notée SIMs(A,B), entre les deux signatures temporelles par les formules suivantes : 25 SI Ms(A, B) = EIY_11/Vi - D-1(Si, Si') (7) e-(to-ti) (8) D -1 (Si, = 1 (9) max (Ey_iIRSSIo-RSSIoll,sat) 30 Dans la formule (9), le terme sat est une constante utilisée pour éviter une division par 0 et ainsi borner la valeur de SIMs(A,B). Dans les formules (3), (6) et (9), la fonction U1 représente l'inverse d'une distance séparant deux composantes d'une signature temporelle. Il 5 (4) (5) (6) 27 convient de noter que d'autres fonctions pourraient être choisies pour définir une telle distance. De la même manière, d'autres fonctions pourraient être choisies pour définir un facteur W1 de pondération en fonction du temps écoulé. Un tel choix ne constitue qu'une variante de mise en oeuvre de l'invention. 5 Dans un mode particulier de mise en oeuvre, la signature temporelle d'un terminal 70 comporte à la fois des valeurs représentatives des phases de mobilité du terminal 70 (telles que celles décrites en référence à la figure 4) et des signatures radio du terminal déterminées à différents instants pendant la période d'observation (telles que celles décrites en référence la figure 7). Un 10 critère de similarité SIM(A,B) peut alors par exemple être défini comme une combinaison convexe d'un critère de similarité SIMm(A,B) pour les phases de mobilité et d'un critère de similarité SIMs(A,B) pour les signatures radio : SIM(A,B) = - SIMm(A,B)+ (1- fi') - SIMs(A,B) (10) " " 15 (t0 -t) - m.b Lii=1 e b2 Dans la formule (11), M est le nombre de valeurs Mi et Mi' dans les signatures temporelles représentatives des phases de mobilité des terminaux A 20 et B, et b est une valeur constante positive. Pour rappel, chaque valeur Mi vaut respectivement 1 ou 0 selon que le terminal 70 est ou non dans une phase de mobilité à un instant ti. Ainsi, le facteur le prend une valeur plus forte si les terminaux A et B ont eu des phases de mobilité similaires en des instants proches de to, et il 25 prend une valeur faible dans le cas contraire. le joue ainsi un rôle de facteur de pondération dans le calcul de SIM(A,B) en accordant plus d'importance à SIMm(A,B) par rapport à SIMs(A,B) lorsque les terminaux A et B ont été en mouvement ensemble en des instants proches de to. Au contraire, le facteur le accorde plus d'importance à SIMs(A,B) par rapport à SIMm(A,B) lorsque les 30 terminaux A et B n'ont pas eu de phases de mobilité similaires en des instants proches de to. Il est ainsi possible d'accorder plus ou moins d'importance aux similarités concernant les phases de mobilité par rapport à celles concernant l'environnement radio selon que les terminaux ont eu ou non des phases de 28 mobilité similaires récemment par rapport à l'instant to d'estimation de la position d'un terminal 70i d'intérêt. 5. Partitionnement des terminaux Différentes méthodes sont envisageables pour l'étape 56 de 5 partitionnement 56 des terminaux 70 du panel considéré en différents groupes en fonction du (ou des) critère(s) de similarité précédemment déterminé(s). L'objectif de cet étape est de déterminer des groupes de terminaux qui ont une probabilité élevée de s'être déplacés ensemble au cours de la période d'observation Tobs précédant l'instant to d'estimation de la position 10 géographique d'un terminal 70i d'intérêt. Il s'agit donc de grouper ensemble des terminaux 70 pour lesquels les valeurs de similarité deux à deux sont particulièrement élevées. Il existe plusieurs méthodes connues pour diviser les terminaux 70 du panel considéré en différents groupes homogènes, en ce sens que les 15 terminaux 70 de chaque groupe partagent des similarités (c'est-à-dire que la valeur du critère de similarité entre deux terminaux d'un même groupe est relativement grande). Par exemple, les méthodes dites « centroïdes » telles que les algorithmes des « k-moyennes » ou « k-médoïdes » définissent des points 20 particuliers de l'espace des mesures considéré, dits « points centraux », qui maximisent la similarité entre un point central et les terminaux 70 du groupe. De telles méthodes associent généralement un terminal 70 à un groupe au plus (méthodes dites de « hard clustering » dans la littérature anglo-saxonne). Si un terminal 70 n'a pas de similarité suffisante avec d'autres terminaux 70 du 25 panel, il n'est associé à aucun groupe et reste isolé. Dans une variante, certaines méthodes de partitionnement associent un terminal 70 à plusieurs groupes selon un certain degré correspondant à une probabilité pour ledit terminal 70 d'appartenir à un groupe (méthodes dites de « soft clustering » dans la littérature anglo-saxonne). 30 Dans les modes de mise en oeuvre où la signature temporelle d'un terminal 70 est une combinaison de deux ensembles de valeurs indépendants, comme par exemple la combinaison de valeurs représentatives de phases de mobilité et de signatures radio, plusieurs méthodes sont possibles pour faire le 29 partitionnement. Selon un premier exemple, le partitionnement est fait en fonction d'un seul critère de similarité prenant en compte les deux ensembles, comme par exemple celle définie ci-avant par la formule (10). 5 Selon un deuxième exemple, un partitionnement peut être fait pour chaque ensemble indépendant selon un critère de similarité qui lui est propre, et un consensus peut être opéré entre les différents partition nements obtenus. Un troisième exemple repose sur une « indépendance conditionnelle » entre des partitionnements effectués selon les critères de similarité qui sont 10 propres à chaque ensemble. Par exemple, un partitionnement est fait pour l'un des deux ensembles à partir d'un critère de similarité S/Mi qui lui est propre. On obtient alors par exemple, pour un terminal A, un vecteur PA = (PA1, PA2, PA3, PAK) où PAi représente la probabilité d'appartenance du terminal A au groupe d'indice j. Puis, les valeurs prises par un critère de similarité S/M2 15 propre à l'autre ensemble sont modifiées en fonction du partitionnement obtenu, par exemple comme suit : siPB Pc = SIM2(A, B) (12) SIM2(A, B) sinon 20 où pc est un seuil prédéterminé. Puis un partitionnement final est fait à partir des valeurs de similarité ainsi modifiées pour le critère S/M2. Les méthodes présentées ci-avant pour faire un partitionnement 56 des terminaux 70 du panel ont été données à titre d'exemples nullement limitatifs. D'autres méthodes pourraient être envisagées, et le choix d'une 25 méthode particulière ne constitue qu'une variante de l'invention. 6. Détermination d'un panel de terminaux Le panel à considérer pour les étapes 52, 54 et 56 peut par exemple être déterminé par une sélection des terminaux 70 du système 60 de communication sans fil dont on sait a priori qu'ils sont proches 30 (géographiquement parlant) du terminal 70i d'intérêt, ou qu'ils ont une certaine propension à appartenir à un groupe de plusieurs terminaux 70. Dans des modes particuliers de mise en oeuvre, ne sont sélectionnés dans le panel que les terminaux 70 qui, lors de l'émission de leur dernier 30 message avant l'instant to, étaient sous la couverture d'au moins une station de base 81 du système 60 de communication sans fil qui couvre également le terminal 70i d'intérêt à l'instant to. Dans des modes particuliers de mise en oeuvre, on note BS 5 l'ensemble des Nb s stations de base 81 qui ont reçu le message émis par le terminal 70i d'intérêt à l'instant to, on note Nm,t(A) le nombre de phases de mobilité d'un terminal A pendant la période d'observation Tobs, et on définit par exemple : s\t 1 10 Pradio (A) = Nb 1-r 1bs E BS (EB#A max ( IRSS/Ams-RSS/Bmsl, sat)) (13) 1 (Nmvt(A)-Nmvt(B))2 (14) e- LB#A Pmvt(A) 2.0-2 Ptot(A) = Pradio(A) - P(A) (15) 15 Ne sont alors sélectionnés dans le panel que les terminaux pour lesquels P0(A) est supérieur à un seuil prédéterminé. Dans la formule (13), RSS/x,bs, bs E BS, est le niveau RSSI mesuré par la station de base bs pour le dernier message émis par un terminal X avant to. 20 Praoio(A) est ainsi une valeur représentative d'une probabilité pour un terminal A de trouver, parmi d'autres terminaux 70 du système 60, des terminaux 70 présentant un environnement radio proche de celui du terminal A. Dans la formule (14), 0- correspond par exemple à l'écart-type des valeurs Nmvt(X). Pmvt(A) est ainsi une valeur représentative d'une probabilité 25 pour un terminal A de trouver, parmi d'autres terminaux 70 du système 60, des terminaux 70 présentant un nombre de phases de mobilité pendant la période Tobs proche de celui du terminal A. Le terme P0(A) définit par la formule (15) représente ainsi une propension pour un terminal A de trouver, parmi d'autres terminaux 70 du 30 système 60, des terminaux 70 qui ont été soumis à des conditions similaires en termes d'environnement radio et en termes de déplacements. Autrement dit, plus P0(A) est grand, et plus la probabilité de pouvoir grouper le terminal A 31 avec d'autres terminaux 70 du système 60 est forte. Il convient de noter que les méthodes présentées ci-avant pour déterminer un panel de terminaux 70 ont été donné à titre d'exemples nullement limitatifs. Aussi, le choix d'une méthode particulière ne constitue 5 qu'une variante de la présente invention. 7. Estimation de la position géographique d'un terminal d'intérêt La phase 40 de recherche du procédé 10 de géolocalisation selon l'invention comporte une étape d'estimation 44 de la position géographique du terminal 70i d'intérêt en fonction de données de géolocalisation disponibles 10 pour au moins un autre terminal 70 d'un groupe, dit « groupe d'intérêt », auquel le terminal 70i d'intérêt appartient. Dans la suite de la description, on considère à titre d'exemple non limitatif que la signature temporelle d'un terminal 70 est une combinaison de valeurs représentatives des phases de mobilité et de signatures radio dudit 15 terminal 70 pendant une période d'observation Tobs précédant un instant to d'estimation de la position du terminal 70i d'intérêt. On considère, de manière avantageuse, que l'instant to correspond à un instant de réception par le serveur 82 d'un message émis par ledit terminal 70i d'intérêt pendant une phase d'immobilité. 20 On considère également qu'un panel de terminaux 70 est déterminé selon l'une quelconque des méthodes décrites ci-avant. On considère également qu'un critère de similarité tel que celui décrit ci-avant par la formule (10) est utilisé pour faire un partitionnement des terminaux 70 du panel. 25 On considère également qu'une position géographique d'un terminal 70 du panel peut être estimée en utilisant par exemple un algorithme conventionnel de géolocalisation tel que celui décrit en référence à la figure 2. Comme cela a déjà été mentionné précédemment, une telle estimation de la position géographique peut être complexe et relativement peu précise, 30 notamment quand la zone géographique à couvrir est vaste et/ou que le nombre de terminaux et/ou de stations de base est élevé. La suite de la description présente à titre d'exemples non limitatifs plusieurs méthodes pour estimer avec plus de précision et/ou avec moins de 32 complexité la position géographique du terminal 70i d'intérêt à l'aide du partitionnement obtenu et de données de géolocalisation fournies par un ou plusieurs terminaux 70 du groupe d'intérêt, voire d'autres groupes. Dans des modes de mise en oeuvre particuliers, les données de 5 géolocalisation sont des positions GPS fournies par des terminaux 70c de calibration équipés d'un récepteur GPS. Selon un premier exemple, le terminal 70i d'intérêt est associé à un seul groupe (méthode dite de « hard clustering »), et au moins un terminal 70c de calibration appartenant à ce groupe a fourni sa position GPS via un 10 message envoyé au réseau d'accès 80 à un instant précédant to depuis lequel il n'a pas bougé. La position géographique du terminal 70i d'intérêt est alors associée à ladite position GPS (ou bien à une moyenne de plusieurs positions GPS disponibles pour le groupe). Selon un deuxième exemple, le terminal 70i d'intérêt est associé à 15 plusieurs groupes avec une certaine probabilité d'appartenance à chaque groupe (méthode dite de « soft clustering »), et une position géographique est attribuée à chaque groupe à partir d'au moins une position GPS fournie par un terminal 70c de calibration dudit groupe (ou éventuellement en moyennant plusieurs positions GPS disponibles pour le groupe). Soit Xk la position 20 géographique attribuée à un groupe k, xj1. Notamment, plusieurs méthodes existent pour déterminer une signature radio pour un terminal 70 : à partir de niveaux RSSI d'un signal échangé entre le terminal 70 et une station de base 81, à partir d'un temps de propagation d'un signal échangé entre le terminal 70 et une station de base 81, 5 etc. Le choix d'une méthode particulière n'est qu'une variante de l'invention. Aussi, plusieurs méthodes existent pour déterminer une signature temporelle pour un terminal 70. Le choix d'un type particulier de mesure (détection de phases de mobilité, valeurs de température, de pression, etc.) n'est qu'une variante de l'invention. 10 De la même manière, plusieurs méthodes existent pour déterminer un critère de similarité entre deux signatures temporelles. Le choix d'une méthode particulière n'est qu'une variante de l'invention. Plusieurs méthodes sont également envisageables pour partitionner les terminaux 70 en plusieurs groupes en fonction d'un critère de similarité, ou 15 bien pour déterminer un panel de terminaux 70 à partitionner. Là encore, le choix d'une méthode particulière n'est qu'une variante de l'invention. Aussi, différents types de données de géolocalisation fournies par certains terminaux 70 du système sont envisageables. Une position GPS d'un terminal 70 a été donnée comme exemple de donnée de géolocalisation, mais 20 d'autres données de géolocalisation sont envisageables. Par exemple, l'adresse MAC (acronyme anglais de « Medium Access Control ») d'un point d'accès Wi-Fi auquel un terminal 70 est connecté peut jouer le rôle de donnée de géolocalisation si la position géographique dudit point d'accès est connue. Le concept de l'invention peut également être appliqué au cas de la 25 géolocalisation d'un seul terminal 70i d'intérêt au cours du temps, indépendamment d'autres terminaux 70 du système 60. Dans le cas par exemple où aucune phase de mobilité du terminal 70i d'intérêt n'est détectée entre deux instants ti et t2 auxquels correspondent respectivement des signatures radio Si et S2, si les deux signatures radio Si et S2 sont similaires, 30 alors elles peuvent être toutes les deux utilisées pour estimer une seule et même position du terminal 70i d'intérêt. Par extension, un groupe de N 2 signatures radio Si similaires et déterminées à des instants ti (l'indice i variant entre 1 et N) entre lesquels le terminal 70i d'intérêt n'a pas bougé peut être 36 utilisé. Une telle approche peut permettre d'améliorer la précision de l'estimation de la position géographique du terminal 70i d'intérêt ainsi que corriger des estimations effectuées dans le passé. Si en revanche des signatures radio ne sont pas similaires, on peut supposer qu'une phase de 5 mobilité a été manquée et il convient d'estimer une position géographique différente pour chaque signature radio. L'invention a été décrite à titre d'exemple en considérant un système 60 de communication sans fil de type loT pour une application dans le domaine de la logistique (transport de palettes équipées avec des terminaux du système 10 60 de communication loT). Rien n'exclut cependant, suivant d'autres exemples, de considérer d'autres système de communication sans fil, comme par exemple des réseaux de téléphonie mobiles, et d'autres applications éventuellement dans d'autres REVENDICATIONS1. Procédé (10) d'estimation de la position géographique d'un terminal, dit « terminal (70i) d'intérêt », parmi un panel de terminaux (70) d'un système (60) de communication sans fil, caractérisé en ce que ledit procédé comporte : une détermination (52), pour chaque terminal (70) dudit panel, d'une signature temporelle comportant des valeurs mesurées pour ledit terminal (70) pendant une durée d'observation prédéterminée, un calcul (54) de valeurs d'un critère de similarité entre les signatures temporelles des terminaux (70) du panel, une valeur dudit critère de similarité calculée pour deux terminaux (70) étant représentative de la probabilité pour lesdits terminaux (70) de s'être déplacés ensemble ou d'avoir été localisés au même endroit pendant ladite période d'observation, un partitionnement (56) des différents terminaux (70) du panel en différents groupes en fonction des valeurs du critère de similarité, une estimation (44) de la position géographique du terminal (70i) d'intérêt en fonction de données de géolocalisation disponibles pour au moins un autre terminal (70) du groupe, dit « groupe d'intérêt », auquel le terminal (70i) d'intérêt appartient. 2. Procédé (10) selon la revendication 1 dans lequel les valeurs mesurées pour déterminer une signature temporelle d'un terminal (70) comportent : des valeurs mesurées par au moins un capteur dudit terminal (70) pendant la période d'observation, et/ou des signatures radio comportant des valeurs représentatives de la qualité de liens radio existant entre le terminal (70) et une ou plusieurs stations de base (81) du système (60) de communication sans fil à différents instants pendant la période d'observation. 3. Procédé (10) selon la revendication 2 dans lequel des valeurs mesurées par au moins un capteur dudit terminal (70) comportent : 38 des valeurs mesurées par un capteur de mouvement représentatives de phases de mobilité dudit terminal (70), et/ou des valeurs mesurées par un capteur de température, et/ou des valeurs mesurées par un capteur de pression atmosphérique, et/ou 5 des valeurs mesurées par un capteur de luminosité, et/ou des valeurs mesurées par un capteur de champ magnétique. 4. Procédé (10) selon la revendication 3 dans lequel une signature temporelle d'un terminal (70) comporte à la fois des valeurs représentatives des phases de mobilité et des signatures radio, et dans lequel le critère de 10 similarité est déterminé en fonction d'un critère de similarité relatif aux phases de mobilité et d'un critère de similarité relatif aux signatures radio. 5. Procédé (10) selon l'une des revendications 1 à 4 dans lequel un critère de similarité est déterminé en fonction d'un facteur de pondération affecté à une valeur mesurée représentatif du temps écoulé entre l'instant où ladite 15 valeur est mesurée et l'instant d'estimation (44) de la position géographique du terminal (70i) d'intérêt. 6. Procédé (10) selon l'une des revendications 1 à 5 dans lequel le partitionnement (56) des différents terminaux (70) du panel est fait de sorte que chaque terminal (70) appartient au plus à un groupe. 20 7. Procédé (10) selon l'une des revendications 1 à 5 dans lequel le partitionnement (56) des différents terminaux (70) du panel est fait de sorte que chaque terminal (70) appartient à chaque groupe selon une valeur de probabilité d'appartenance, le groupe d'intérêt correspondant au groupe pour lequel la valeur de probabilité d'appartenance du terminal (70i) 25 d'intérêt est la plus grande. 8. Procédé (10) selon la revendication 7 dans lequel la position géographique du terminal (70i) d'intérêt est estimée en fonction de données de géolocalisation disponibles pour différents terminaux (70) appartenant à différents groupes et en fonction des valeurs de probabilité d'appartenance 30 du terminal (70i) d'intérêt aux différents groupes. 9. Procédé (10) selon l'une des revendications 1 à 8 dans lequel lesdites données de géolocalisation sont fournies par un terminal (70) sous la forme 39 d'une position géographique obtenue par un système de positionnement dudit terminal. 10. Procédé (10) selon l'une des revendications 1 à 8 dans lequel des données de géolocalisation sont disponibles pour un terminal (70) sous la forme 5 d'une signature radio comportant des valeurs représentatives de la qualité de liens radio existant entre ledit terminal (70) et une ou plusieurs stations de base (81) du système (60) de communication sans fil. 11. Procédé (10) selon la revendication 10 dans lequel une position géographique d'un terminal (70) est estimée en fonction de la signature 10 radio dudit terminal (70) et d'un ensemble (26) de référence comportant des signatures radio associées à des positions géographiques connues, et dans lequel la position géographique du terminal (70i) d'intérêt est estimée en fonction de la position géographique estimée d'au moins un autre terminal (70) du groupe d'intérêt. 15 12. Procédé (10) selon la revendication 10 dans lequel une signature radio virtuelle est déterminée pour le groupe d'intérêt en fonction des signatures radio des terminaux (70) appartenant audit groupe d'intérêt, et la position du terminal (70i) d'intérêt est estimée en fonction de ladite signature radio virtuelle. 20 13. Procédé (10) selon l'une des revendications 1 à 12 dans lequel ledit panel de terminaux (70) est déterminé par une étape de sélection parmi les terminaux (70) du système (60) de communication sans fil de sorte qu'un terminal (70) du panel soit couvert par au moins une station de base (81) du système (60) de communication sans fil couvrant également le terminal 25 (70i) d'intérêt. 14.Système (60) de communication sans fil comportant des terminaux (70) et un réseau d'accès (80) comportant des stations de base (81) et un serveur (82) relié auxdites stations de base (81), ledit système 60 étant caractérisé en ce qu'il met en oeuvre un procédé de géolocalisation d'un terminal (70i) 30 d'intérêt, parmi un panel de terminaux (70), selon l'une des revendications 1 à 13, ledit serveur (82) étant configuré pour : - collecter des valeurs mesurées pour un terminal (70) dudit panel pendant une durée d'observation prédéterminée, 40 déterminer une signature temporelle pour ledit terminal (70) à partir des valeurs mesurées, calculer des valeurs d'un critère de similarité entre les signatures temporelles des terminaux (70) du panel, une valeur dudit critère de 5 similarité calculée pour deux terminaux (70) étant représentative de la probabilité pour lesdits terminaux (70) de s'être déplacés ensemble ou d'avoir été localisés au même endroit pendant ladite période d'observation, partitionner les différents terminaux (70) du panel en différents groupes 10 en fonction des valeurs du critère de similarité, estimer la position géographique du terminal (70i) d'intérêt en fonction de données de géolocalisation disponibles pour au moins un autre terminal (70) du groupe, dit « groupe d'intérêt », auquel le terminal (70i) d'intérêt appartient. 15