FR 3074305 A1 20190531 FR 1761431 A 20171130 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). É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, la télérelève d'informations en général, 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 échangés entre le terminal et les stations de base. Ces différentes méthodes 2 présentent cependant toutes l'inconvénient de nécessiter du matériel et du logiciel spécifiques. De plus elles requièrent souvent une synchronisation coûteuse des différentes stations de base jouant le rôle de points d'observation. Enfin, elles sont particulièrement sensibles au phénomène dit de 5 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 puissance reçue (RSSI, de l'anglais « Received Signal Strength Indicator ») 10 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 système de communication lui-même. Ces méthodes reposent sur le fait qu'un 15 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 terminal des différentes stations de base qui l'entourent à partir des niveaux de 20 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, mouvement du terminal, etc.) rendent très complexe la fonction qui définit la distance à partir d'un 25 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, pendant une phase de calibration, une base de données qui associe à des 30 positions géographiques connues une signature radio correspondant à l'ensemble des niveaux de RSSI mesurés pour un terminal à la position considérée pour un ensemble de stations de base du système. Ensuite, pendant une phase de recherche, une signature radio observée pour un 3 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 De tels procédés présentent cependant plusieurs inconvénients, notamment dans le cas où la zone géographique à couvrir est très 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 10 soient significativement différentes, ou bien que deux signatures radio 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 15 l'algorithme d'apprentissage automatique utilisé. En effet, pour obtenir une 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 20 données en entrée de l'algorithme d'apprentissage automatique, et ainsi à créer des contraintes de capacité et de temps de calcul. Pour réaliser la phase de calibration, il est connu d'embarquer dans une flotte de véhicules qui parcourent la zone à 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 25 système de communication en différents points (le terme anglais pour qualifier cette phase est « war-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. 30 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. 4 A cet effet, et selon un premier aspect, il est proposé par la présente invention un procédé de géolocalisation d'un terminal, dit « terminal d'intérêt », d'un système de communication sans fil. Le procédé comporte : une détermination, à partir de messages émis par d'autres terminaux, dits 5 « terminaux de calibration », et reçus par des stations de base dudit système de communication sans fil, d'un ensemble d'éléments de référence, chaque élément de référence comportant une signature radio associée à une position géographique d'un terminal de calibration lors de l'émission d'un message, ladite signature radio correspondant à un 10 ensemble de valeurs représentatives de la qualité de liens radio existant entre ledit terminal de calibration et des stations de base, une détermination, à partir d'un message émis par ledit terminal d'intérêt situé à une position géographique à estimer, d'une signature radio pour ledit terminal d'intérêt, 15 une identification de positions géographiques particulières, dites « points d'intérêt », à partir des éléments de référence de l'ensemble d'éléments de référence, - une modification de l'ensemble d'éléments de référence à partir des points d'intérêt identifiés, 20 - une estimation de la position géographique dudit terminal d'intérêt à partir de la signature radio dudit terminal d'intérêt et de l'ensemble d'éléments de référence modifié. L'estimation de la position géographique d'un terminal d'intérêt est par exemple réalisée par un algorithme d'apprentissage automatique qui a pour but 25 d'établir une relation entre une signature radio et une position géographique. Un point d'intérêt correspond à un point prioritaire d'émission de messages où, pour les terminaux du système de communication sans fil considéré, il y aura en moyenne plus de messages émis depuis ou à proximité d'un point d'intérêt que depuis une autre position géographique. 30 Par exemple, un point d'intérêt peut correspondre à une position géographique particulière à laquelle la probabilité de présence d'un terminal du système considéré est plus grande qu'ailleurs. Notamment, dans le cas du transport de palettes de marchandises équipées avec des terminaux et se 5 déplaçant d'entrepôt en entrepôt avec un temps de voyage relativement cours par rapport au temps passé dans un entrepôt, lesdits entrepôts peuvent correspondre à des points d'intérêt. Selon un autre exemple, un point d'intérêt peut correspondre à un 5 point de convergence de terminaux du système. Notamment, des aéroports peuvent correspondre à des points d'intérêt pour des bagages équipés avec des terminaux du système. Les points d'intérêt peuvent être identifiés dynamiquement à partir des éléments de référence fournis par les terminaux de calibration. Dans des 10 modes particuliers de mise en oeuvre, au moins une partie des points d'intérêt peuvent également être connus a priori. La modification de l'ensemble d'éléments de référence correspond par exemple à une phase d'affinement au cours de laquelle des éléments de référence disparates indésirables sont supprimés (par exemple ceux qui sont 15 jugés éloignés des points d'intérêt) et/ou des informations supplémentaires sont ajoutées à des éléments de référence pertinents (par exemple ceux qui sont jugés proches des points d'intérêt), afin de nettoyer et/ou d'enrichir l'ensemble d'éléments de référence. De telles dispositions permettent de gagner en précision de géolocalisation car cette phase d'affinement renforce la 20 relation existant entre une signature radio et une position géographique en tirant profit de l'existence de points d'intérêt auxquels la probabilité de présence d'un terminal est plus grande qu'ailleurs. En outre, la suppression d'éléments disparates contribue à diminuer le nombre d'éléments de référence utilisés par l'algorithme d'apprentissage, et 25 ainsi à en réduire la complexité, ce qui se traduit par un gain en termes de capacité et de temps de calcul. Enfin, ces dispositions peuvent permettre de réduire la taille de la signature radio, par exemple en ne considérant que les stations de base qui sont localisées à proximité des points d'intérêt, ce qui contribue encore 30 davantage à réduire la complexité de l'algorithme d'apprentissage. Dans des modes particuliers de mise en oeuvre, l'invention peut comporter en outre l'une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons techniquement possibles. 6 Dans des modes particuliers de mise en oeuvre, la position géographique du terminal d'intérêt est estimée sous la forme d'une valeur représentative d'une probabilité pour ledit terminal d'intérêt d'être situé en un point d'intérêt identifié. 5 De telles dispositions permettent par exemple de conclure que le terminal d'intérêt est situé au point d'intérêt pour lequel la valeur de probabilité est la plus forte si cette valeur est supérieure à un certain seuil, ou bien de conclure, si cette valeur est inférieure audit seuil, que le terminal d'intérêt n'est situé à proximité d'aucun point d'intérêt (il est par exemple en transit entre deux 10 points d'intérêt). Dans des modes particuliers de mise en oeuvre, l'identification des points d'intérêt comporte : - un découpage d'une zone géographique à couvrir en plusieurs cellules de tailles prédéfinies, 15 - un calcul pour chaque cellule d'une valeur représentative d'un nombre d'éléments de référence situés dans ladite cellule, si la valeur calculée satisfait un critère prédéfini, une définition d'un point d'intérêt en fonction de ladite cellule et/ou des positions géographiques des éléments de référence situés dans ladite cellule. 20 Dans des modes particuliers de mise en oeuvre, l'identification des points d'intérêt comporte une étape de calcul, pour chaque élément de référence, de distances séparant ledit élément de référence et chaque autre élément de référence, ainsi que les étapes suivantes effectuées itérativement sur un ensemble d'éléments de référence temporaire correspondant 25 initialement à l'ensemble d'éléments de référence, jusqu'à ce qu'une condition d'arrêt soit atteinte : un calcul pour chaque élément de référence de l'ensemble d'éléments de référence temporaire, d'une valeur représentative d'un nombre d'éléments de référence situés à une distance inférieure à un seuil prédéfini, dits 30 « éléments voisins », si la plus grande valeur calculée, pour un élément dit « élément d'intérêt », est supérieure à un seuil prédéfini : 7 o une définition d'un point d'intérêt en fonction des positions géographiques de l'élément d'intérêt et/ou de ses éléments voisins, o une suppression de l'élément d'intérêt et de ses éléments voisins de l'ensemble d'éléments de référence temporaire. 5 Dans des modes particuliers de mise en oeuvre, le procédé de géolocalisation comporte en outre une étape de filtrage des points d'intérêt identifiés en fonction de leurs positions géographiques et/ou d'une valeur représentative d'un nombre d'éléments de référence de l'ensemble d'éléments de référence situés à une distance inférieure à un seuil prédéfini. 10 Dans des modes particuliers de mise en oeuvre, la valeur représentative d'un nombre d'éléments de référence est calculée en fonction d'un facteur de pondération affecté à chaque élément de référence et représentatif d'une durée écoulée entre le moment d'émission du message par le terminal de calibration et le moment où les points d'intérêt sont identifiés. 15 De telles dispositions permettent notamment de favoriser des points d'intérêt récents en accordant plus d'importance (c'est-à-dire en affectant une pondération plus forte) aux éléments de référence obtenus à un instant proche de l'instant d'estimation de la positon du terminal d'intérêt. Dans des modes particuliers de mise en oeuvre, la modification de 20 l'ensemble d'éléments de référence comporte les étapes suivantes pour chaque élément de référence de l'ensemble d'éléments de référence : un calcul de la plus petite distance séparant ledit élément de référence et un des points d'intérêt identifiés, si la plus petite distance calculée est inférieure à un seuil prédéfini : 25 o une association audit élément de référence du point d'intérêt correspondant, - sinon : o une suppression dudit élément de référence de l'ensemble d'éléments de référence. 30 Dans des modes particuliers de mise en oeuvre, la modification de l'ensemble d'éléments de référence comporte les étapes suivantes pour chaque élément de référence de l'ensemble d'éléments de référence : 8 un calcul, pour chaque point d'intérêt identifié, d'une valeur représentative de la distance qui sépare ledit élément de référence dudit point d'intérêt, une association audit élément de référence des valeurs ainsi calculées. Dans des modes particuliers de mise en oeuvre, un lien radio est un 5 lien montant vers une station de base du système de communication sans fil. Dans des modes particuliers de mise en oeuvre, une détermination de la signature radio d'un terminal comporte : une émission par le terminal d'un message à destination des stations de base, 10 une mesure, pour chaque station de base, d'une valeur représentative de la qualité du lien radio existant entre ledit terminal et la station de base considérée à partir du message reçu dudit terminal, une formation, par un serveur relié aux stations de base, de la signature radio du terminal à partir des valeurs mesurées. 15 Dans des modes particuliers de mise en oeuvre, un terminal de calibration est équipé d'un système de positionnement et la détermination de la signature radio dudit terminal de calibration comporte en outre : une inclusion, dans le message émis par ledit terminal de calibration, de la position géographique courante dudit terminal de calibration mesurée 20 par le système de positionnement, une extraction, par le serveur, de la position géographique contenue dans ledit message. De telles dispositions permettent notamment de pouvoir facilement enrichir l'ensemble d'éléments de référence avec de nouveaux éléments de 25 référence provenant des terminaux de calibration, ce qui permet de découvrir de nouveaux points d'intérêt pouvant apparaître au fil du temps et/ou d'oublier des points d'intérêt devenus désuets. La phase de calibration consistant à définir ou enrichir l'ensemble d'éléments de référence peut être effectuée une seule fois au début, ou bien 30 elle peut être répétée de manière récurrente (par exemple périodiquement ou bien de manière contrôlée par envoi de messages de contrôle aux terminaux de calibration), ou bien elle peut être effectuée de façon continue (par exemple si les terminaux de calibration émettent régulièrement des messages incluant 9 leur position géographique courante). Grâce à l'existence de certains terminaux spécifiques équipés d'un système de positionnement, il est possible d'automatiser entièrement la phase de calibration qui consiste à définir l'ensemble d'éléments de référence. En 5 outre, la mise en oeuvre du procédé de géolocalisation au sein du système de communication sans fil est particulièrement aisée et peu onéreuse car elle ne nécessite pas de matériel spécifique additionnel au niveau des stations de base ou au niveau de la majorité des terminaux (seuls les terminaux jouant le rôle de terminaux de calibration peuvent nécessiter des modifications 10 matérielles et/ou logicielles). Dans des modes particuliers de mise en oeuvre, la valeur représentative de la qualité du lien radio existant entre un terminal et une station de base est un niveau de puissance reçue d'un signal radio échangé entre ladite station de base et ledit terminal. 15 Dans des modes particuliers de mise en oeuvre, le lien radio est un canal de communication à bande ultra étroite. Un système de communication sans fil à bande ultra étroite ("Ultra Narrow Band" ou UNB dans la littérature anglo-saxonne) est particulièrement adapté pour des applications du type loT. Il permet en effet des 20 communications bas-débit avec une longue portée et des consommations énergétiques basses (système de type LPWAN, pour "Low Power Wide Area Network" en anglais). Dans des modes particuliers de mise en oeuvre, l'estimation de la position géographique d'un terminal d'intérêt à partir de la signature radio dudit 25 terminal d'intérêt et de l'ensemble d'éléments de référence modifié est réalisée à l'aide d'un algorithme d'apprentissage automatique de classification basé sur une technique d'arbres décisionnels. 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 30 comportant des stations de base et un serveur relié auxdites stations de base. Ledit système comporte une base de données mémorisant un ensemble d'éléments de référence formés à partir de messages émis par des terminaux, dits « terminaux de calibration », et reçus par au moins une station de base. 10 Chaque élément de référence comporte une signature radio associée à une position géographique d'un terminal de calibration lors de l'émission d'un message. Ladite signature radio comprend des valeurs représentatives de la qualité de liens radio existant entre ledit terminal de calibration et des stations 5 de base. En outre, le réseau d'accès est configuré pour : - mesurer des valeurs représentatives de la qualité de liens radio existant entre des stations de base et un terminal dont la position géographique doit être estimée, dit « terminal d'intérêt », - former une signature radio correspondant à l'ensemble des valeurs 10 mesurées pour ledit terminal d'intérêt, - identifier des positions géographiques particulières, dites « points d'intérêt », à partir de l'ensemble d'éléments de référence, modifier l'ensemble d'éléments de référence à partir des points d'intérêt identifiés, 15 estimer la position géographique du terminal d'intérêt à partir de la signature radio du terminal d'intérêt et de l'ensemble d'éléments de référence modifié. PRÉSENTATION DES FIGURES 20 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 à 14 qui représentent : - Figure 1 : une représentation schématique d'un système de communication sans fil comportant un dispositif de calibration, 25 - Figure 2 : un diagramme illustrant les étapes principales d'un procédé conventionnel de géolocalisation d'un terminal d'un système de communication sans fil, - Figure 3 : un diagramme illustrant les étapes principales d'un procédé de géolocalisation selon l'invention, 30 - Figure 4 : une représentation schématique d'un mode préféré de réalisation où le dispositif de calibration est un terminal équipé d'un système de positionnement, - Figure 5 : un diagramme illustrant les étapes principales d'un mode préféré 11 de mise en oeuvre de la phase de calibration du procédé de géolocalisation, - Figure 6 : un diagramme illustrant les étapes principales d'un mode préféré de mise en oeuvre de la phase de recherche du procédé de 5 géolocalisation, - Figure 7 : un diagramme illustrant les étapes principales d'un premier mode particulier de mise en oeuvre d'une étape d'identification de points d'intérêt, - Figure 8 : une illustration du premier mode particulier de mise en oeuvre de l'étape d'identification de points d'intérêt, 10 - Figure 9 : un diagramme illustrant les étapes principales d'un deuxième mode particulier de mise en oeuvre de l'étape d'identification de points d'intérêt, - Figure 10 : une illustration du deuxième mode particulier de mise en oeuvre de l'étape d'identification de points d'intérêt, 15 - Figure 11 : un diagramme illustrant les étapes principales d'un premier mode particulier de mise en oeuvre d'une étape de modification d'un ensemble de référence, - Figure 12 : un diagramme illustrant les étapes principales d'un deuxième mode particulier de mise en oeuvre de l'étape de modification de 20 l'ensemble de référence, - Figure 13 : des courbes de répartition d'une distance d'erreur pour un procédé conventionnel de géolocalisation, - Figure 14: des courbes de répartition d'une distance d'erreur pour le procédé de géolocalisation selon l'invention. 25 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 à l'échelle, sauf mention contraire. 30 DESCRIPTION DÉTAILLÉE DE MODES DE RÉALISATION Tel qu'indiqué précédemment, la présente invention vise à estimer la position géographique d'un terminal, dit « terminal d'intérêt », d'un système de communication sans fil à l'aide d'une méthode d'apprentissage basée sur des 12 informations représentatives de la qualité de liens radio existant entre ledit terminal et des stations de base dudit système de communication sans fil. La figure 1 représente schématiquement un système 60 de communication sans fil, comportant plusieurs terminaux 70 et un réseau 5 d'accès 80 comportant plusieurs stations de base 81 reliées et un serveur 82 reliées auxdites stations de base 81. Dans un tel système 60, les communications peuvent généralement être bidirectionnelles, c'est-à-dire que des données peuvent être transmises du réseau d'accès 80 vers un terminal 70 sur un lien radio descendant, ou bien d'un terminal 70 vers le réseau 10 d'accès 80 sur un lien radio montant. Aussi, une mesure de la qualité du lien radio existant entre un terminal 70 et une station de base 81 du réseau d'accès 80 peut par exemple être effectuée par l'une ou l'autre de ces entités. Aussi, il est possible d'introduire dans le système 60 de communication sans fil un ou plusieurs dispositifs 71 de calibration adaptés 15 pour déterminer avec précision la position géographique courante, par exemple avec un système de positionnement tel qu'un récepteur GPS. Ce dispositif 71 de calibration possède en outre des moyens configurés pour permettre la mesure d'une valeur représentative de la qualité du lien radio existant entre lui et les stations de base 81 du réseau d'accès 80. Comme indiqué 20 précédemment, cette mesure peut par exemple être faite par le dispositif 71 de calibration lui-même, sur le lien descendant, ou bien par une station de base 81 du réseau d'accès 81, sur le lien montant. Si la mesure est effectuée par le dispositif 71 de calibration sur le lien descendant, cette mesure peut par exemple être faite sur des signaux radio destinés à d'autres terminaux 70 du 25 système que le dispositif 71 de calibration. Les valeurs représentatives de la qualité de lien radio existant entre un terminal 70 ou un dispositif 71 de calibration et une station de base 81 peuvent être transmises au serveur 82 afin qu'il mette en oeuvre certaines étapes d'un procédé de géolocalisation. 30 La figure 2 représente les deux principales phases d'un procédé 10 de géolocalisation d'un terminal 70 d'intérêt d'un système 60 de communication sans fil basé sur une méthode d'apprentissage. Une première phase 20 de calibration comporte : 13 - une étape de détermination 22 de signatures radio associées respectivement à des positions géographiques connues, chaque signature radio correspondant à un ensemble de valeurs représentatives de la qualité de liens radio existant entre un 5 dispositif 71 de calibration dont la position géographique est connue et une pluralité de stations de base 81 du réseau d'accès 80 du système 60 de communication sans fil, et - une étape de stockage 24 dans un ensemble 26 d'éléments de référence de chaque signature radio déterminée et de la position 10 géographique associée du dispositif 71 de calibration. Par « position géographique », on entend dans la suite de la description un ensemble de deux coordonnées correspondant respectivement à 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 15 troisième coordonnée correspondant à l'altitude par rapport au niveau moyen de la mer pourrait également être prise en compte. En résumé, il s'agit pendant cette phase 20 de calibration de faire une sorte de cartographie radio de la zone géographique considérée. Ensuite, une phase 40 de recherche comporte : 20 - une étape de détermination 42 d'une signature radio pour ledit terminal 70 d'intérêt situé à une position géographique à estimer, - une étape d'estimation 44 de la position géographique dudit terminal 70 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 25 de calibration. Le procédé 10 de géolocalisation selon l'invention introduit une phase 30 supplémentaire, dite « phase 30 d'affinement », qui a pour but de nettoyer l'ensemble 26 d'éléments de référence (c'est-à-dire supprimer des éléments de référence indésirables) et/ou de l'enrichir (c'est-à-dire d'ajouter des 30 informations aux éléments de référence conservés) afin d'améliorer l'estimation 44 de la position géographique d'un terminal 70 de référence. Cette phase 30 d'affinement se base sur l'idée que les terminaux 70 peuvent avoir tendance à passer par des positions géographiques particulières 14 de convergence des terminaux 70 du système 60 (on entend par là des positions géographiques auxquelles un grand nombre de terminaux 70 ont tendance à passer). A de telles positions géographiques particulières, le nombre de messages émis par les terminaux 70 est alors en moyenne plus 5 élevé qu'ailleurs. Dans certains cas, les terminaux 70 peuvent également occuper préférentiellement un ensemble fini de positions géographiques particulières (autrement dit, à un instant donné, un terminal 70 peut statistiquement avoir une probabilité plus forte d'être localisé à l'une de ces positions géographiques particulières plutôt qu'ailleurs). Dans la description, on 10 appelle « point d'intérêt » une telle position géographique particulière. Une telle situation est observée par exemple dans le cas du transport de marchandises se déplaçant d'entrepôt en entrepôt avec un temps de voyage relativement cours par rapport au temps passé dans un entrepôt. Il est par exemple possible d'équiper des palettes utilisés pour le transport de 15 marchandises avec des terminaux 70 du système 60 de communication sans fil afin de pouvoir les géolocaliser. Les différents entrepôts dans lesquels sont successivement stockées les palettes correspondent alors à des points d'intérêt pour le procédé 10 de géolocalisation. Un camion transportant lesdites palettes peut par exemple embarquer un dispositif 71 de calibration permettant 20 de construire l'ensemble 26 d'éléments de référence. Pour estimer la position géographique d'un terminal 70 d'intérêt, plutôt que d'estimer la position du terminal 70 sur l'ensemble de la zone géographique à couvrir, le procédé 10 de géolocalisation selon l'invention estime la position du terminal 70 d'intérêt en fonction des points d'intérêt 25 identifiés. Autrement dit, plutôt que d'utiliser une méthode de régression conventionnelle qui nécessiterait un nombre d'éléments de référence très élevé (et donc des capacités et des temps de calcul très importants) pour obtenir une précision raisonnable, le procédé 10 de géolocalisation selon l'invention utilise un algorithme d'apprentissage automatique de classification qui consiste par 30 exemple à attribuer au terminal 70 d'intérêt une probabilité d'être situé à chacun des points d'intérêt identifiés. La position géographique estimée du terminal 70 d'intérêt peut alors par exemple être définie comme la position géographique du point d'intérêt auquel il a la plus grande probabilité d'être 15 situé. En fonction de la valeur de la probabilité remontée (par exemple si elle est inférieure à un seuil prédéfini), il est également possible de déterminer si le terminal 70 n'est pas situé à un point d'intérêt (par exemple s'il s'agit d'un terminal 70 embarqué sur une palette qui est en transit entre deux entrepôts). 5 La figure 3 représente schématiquement les principales phases d'un procédé 10 de géolocalisation selon l'invention. Outre la phase 20 de calibration et la phase 40 de recherche déjà décrites en référence à la figure 2, le procédé 10 de géolocalisation comporte une phase 30 d'affinement. Cette phase 30 d'affinement comporte une première étape 10 d'identification 32 de points d'intérêt qui permet de définir, à partir de l'ensemble 26 d'éléments de référence acquis pendant la phase 20 de calibration, un ensemble 35 de points d'intérêt. Cette étape sera détaillée ultérieurement en référence aux figures 7 à 10. Ensuite, à l'aide de l'ensemble 35 de points d'intérêt identifiés, une 15 étape de modification 34 de l'ensemble 26 d'éléments de référence fournit un ensemble 36 d'éléments de référence modifié qui est alors utilisé pendant la phase 40 de recherche pour déterminer la position géographique d'un terminal 70 d'intérêt. Cette étape de modification 34 sera détaillée ultérieurement en référence aux figures 10 et 11. 20 Il est en effet avantageux d'identifier des points d'intérêt et d'accorder plus d'importance aux éléments de référence dont la position géographique est proche de celle d'un point d'intérêt. Il est par exemple possible de supprimer les éléments de référence dont la position géographique est éloignée d'un point d'intérêt d'une distance supérieure à un seuil prédéfini. En outre, il est 25 avantageux d'associer à chaque élément de référence le point d'intérêt le plus proche, ou bien une valeur représentative de la distance qui sépare la position géographique dudit élément de référence de chaque point d'intérêt. De telles dispositions permettent en effet de réduire un « bruit » généré par les éléments qui ne sont pas regroupés à proximité des points d'intérêt dans la relation que 30 l'algorithme d'apprentissage automatique cherche à établir entre une signature radio et un point d'intérêt. On entend ici par « distance » une distance calculée de manière conventionnelle entre deux positions géographiques exprimées à l'aide de la 16 longitude et de la latitude. Aussi, on entend par « distance entre deux éléments de référence », la distance entre les positions géographiques desdits deux éléments de référence. Il en sera de même tout au long de la description, sauf mention contraire. 5 Dans la suite de la description, on se place à titre d'exemple et de manière non limitative dans le cas d'un système 60 de communication sans fil à 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 radioélectriques émis par les terminaux est de largeur 10 fréquentielle inférieure à deux kilohertz, voire inférieure à un kilohertz. Par « signal radioélectrique », on entend une onde électromagnétique se propageant via des moyens non filaires, dont les fréquences sont comprises dans le spectre traditionnel des ondes radioélectriques (quelques hertz à plusieurs centaines de gigahertz). De tels systèmes de communication sans fil 15 UNB sont particulièrement adaptés pour des applications du type loT ou du type M2M (acronyme anglo-saxon pour « Machine to Machine »). Dans un tel système 60 de communication sans fil, les échanges de données sont essentiellement monodirectionnels, en l'occurrence sur un lien montant des terminaux 70 vers le réseau d'accès 80 dudit système 60 de 20 communication sans fil. Afin de minimiser les risques de perdre un message émis par un terminal 70, la planification du réseau d'accès 80 est souvent réalisée de telle sorte qu'une zone géographique donnée est couverte simultanément par plusieurs stations de base 81, de telle manière qu'un message émis par un terminal 70 peut être reçu par plusieurs stations de base 25 81. Chaque station de base 81 est adaptée à recevoir des messages 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 30 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é 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. 17 Le serveur 82 peut notamment être utilisé pour la mise en oeuvre du procédé 10 de géolocalisation d'un terminal 70 du système 60. La figure 4 représente schématiquement un mode préféré de réalisation d'un tel système 60 dans lequel certains terminaux 72 du système 5 60 de communication sans fil, dit « terminaux 72 de calibration », 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 72 (par exemple quand celui-ci est situé à un endroit où il peut capter les signaux des satellites dudit système de positionnement). Ces terminaux 72 de calibration jouent le 10 rôle des dispositifs 71 de calibration mentionnés précédemment en plus de leurs fonctions courantes dans le système 60 de communication sans fil. Il convient de noter que d'autres méthodes peuvent être envisagées pour obtenir la position géographique d'un terminal 72 de calibration. Par exemple, la position géographique d'un terminal 72 de calibration peut être 15 déterminée à partir de l'adresse MAC (acronyme anglais de « Medium Access Control ») d'un point d'accès Wi-Fi auquel ledit terminal 72 de calibration est connecté si la position géographique dudit point d'accès est connue. Le choix d'une méthode particulière pour déterminer la position géographique d'un terminal 72 de calibration ne constitue qu'une variante de mise en oeuvre de 20 l'invention. La figure 5 représente les principales étapes d'un mode préféré de mise en oeuvre de la phase 20 de calibration du procédé 10 de géolocalisation. Pour cette phase 20 de calibration, plutôt que de parcourir intentionnellement la zone géographique à couvrir avec un appareil spécifique 25 adapté pour déterminer une position géographique et mesurer des valeurs représentatives de qualité de lien radio à cette position, il est avantageux d'utiliser la capacité éventuelle de certains terminaux du système de remplir ce rôle, c'est le cas des terminaux 72 de calibration. Ainsi, une première étape de la phase 20 de calibration consiste en 30 l'émission 220 par un terminal 72 de calibration d'un message contenant la position géographique courante à destination du réseau d'accès 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 18 s'agir d'un message de télérelève d'information conventionnel incluant la position géographique courante et n'ayant pas pour objectif premier de participer à la phase 20 de 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 5 contenant ledit message font alors une mesure 221 de la qualité du lien radio sur lequel le message a été transmis. Dans un mode préféré de mise en oeuvre, et pour la suite de la description à titre d'exemple non limitatif, la valeur représentative de la qualité du lien radio utilisée est le niveau moyen de puissance reçue (RSSI), par exemple exprimé en décibel, par la station de 10 base 81 pour le signal transportant ledit message. Il est à noter que d'autres valeurs représentatives de la qualité du lien radio pourraient être utilisées, comme par exemple l'atténuation du signal, un rapport signal sur bruit du signal (SNR ou « Signal on Noise Ratio » dans la littérature anglo-saxonne) ou bien un indicateur de qualité du canal de 15 communication (COI ou « Channel Quality Indicator » en anglais). Dans une étape 222 de formation de la signature radio, le serveur 82 collecte alors 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 72 de calibration. Pour les stations de base 81 qui n'ont pas reçu le message, par 20 exemple parce que le terminal 72 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 plusieurs variantes sont possibles pour le choix des stations de base à prendre en compte pour la formation d'une signature radio. Selon un premier exemple, de manière conventionnelle, toutes 25 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 s'intéresse. Le serveur 82 procède aussi à l'extraction 224 de l'information sur la 30 position géographique contenue dans le message émis par le terminal 72 de calibration (pour rappel, ce message a été préalablement transmis au serveur 82 par les stations de base 81 qui l'ont reçu). Finalement, le couple d'information formé par la position géographique 19 du terminal 72 de calibration et la signature radio associée est ajouté à l'ensemble 26 des éléments de référence acquis pendant la phase 20 de calibration. Par exemple, cette étape consiste au stockage 24 du couple d'information formé par la position géographique et la signature radio associée 5 dans une base de données mémorisée au sein du serveur 82. 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 plusieurs mois, pour des terminaux 72 de calibration du système 60 de 10 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 d'éléments de référence dans l'ensemble 26. 15 Dans un mode préféré de mise en oeuvre, l'ensemble 26 d'éléments de référence est continuellement enrichi par de nouveaux éléments de référence provenant des terminaux 72 de calibration. De telles dispositions permettent notamment de découvrir de nouveaux points d'intérêt pouvant apparaître au fil du temps et/ou d'oublier des points d'intérêt devenus désuets. 20 Il convient de noter que, tout comme les terminaux 70 d'intérêt, les terminaux 72 de calibration sont mobiles et ils peuvent par conséquent fournir des données de référence correspondant à différentes positions géographiques. Dans des modes particuliers de mise en oeuvre, l'ensemble 26 25 d'éléments de référence est également continuellement enrichi par des éléments correspondant chacun à une position géographique estimée d'un terminal 70 d'intérêt associée à la signature radio dudit terminal 70 d'intérêt. Il convient également de noter que selon d'autres modes de mise en oeuvre, l'information sur la position géographique courante du terminal 72 de 30 calibration pourrait être émise par un autre système de communication que celui pour lequel sont faites les mesures de RSSI permettant de définir la signature radio dudit terminal 72 de calibration. Par exemple, la position géographique courante du terminal 72 de calibration pourrait être émise au 20 serveur 82 via un réseau de téléphonie mobile du type GSM, UMTS ou LTE, alors que le message pour lequel sont effectuées les mesures de RSSI est émis vers le réseau d'accès 80 du système 60 de communication sans fil UNB. Afin de faciliter, au niveau du serveur 82, l'association de la signature radio 5 déterminée pour le terminal 72 de calibration avec sa position géographique courante, il est envisageable, par exemple, d'horodater le message transportant l'information sur la position géographique courante et le message émis vers le système de communication sans fil UNB. La figure 6 représente les principales étapes d'un mode préféré de 10 mise en oeuvre de la phase 40 de recherche pour le procédé 10 de géolocalisation. Cette phase 40 de recherche est initiée par l'émission 420 d'un message par un terminal 70 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 15 quelconque émis de manière indépendante au procédé 10 de géolocalisation du terminal. Par exemple, ce message peut être émis à 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 20 nécessairement d'importance dans la phase 40 de recherche. 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 421 du niveau de RSSI du signal transportant le message. Dans une étape 422 de formation de la signature radio, le serveur 82 25 collecte alors 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 70 d'intérêt. L'estimation 44 de la position géographique du terminal 70 d'intérêt ayant émis le message est alors effectuée à partir d'une part de l'ensemble 36 30 d'éléments de référence modifié par la phase 30 d'affinage, et d'autre part de la signature radio déterminée pour le terminal 70 d'intérêt dont la position géographique doit être estimée. La figure 7 représente schématiquement les principales étapes d'un 21 premier mode particulier de mise en oeuvre de l'étape d'identification 32 des points d'intérêt. Cette étape d'identification 32 des points d'intérêt est mise en oeuvre par le serveur 82 pendant la phase 30 d'affinement à partir de l'ensemble 26 des éléments de référence acquis pendant la phase 20 de 5 calibration. Une première étape correspond à un découpage 320 de la zone géographique à couvrir en plusieurs cellules dont les dimensions sont prédéfinies. Ce découpage peut correspondre par exemple à un quadrillage régulier de la zone géographique à couvrir, les cellules correspondant alors à 10 des carrés identiques. Selon un autre exemple, les cellules peuvent avoir une forme circulaire avec un diamètre fixe prédéfini et un recouvrement partiel de chaque cellule avec au moins une autre cellule afin de couvrir la totalité de la zone géographique. Selon encore un autre exemple, les cellules peuvent avoir des formes et/ou des dimensions différentes définies par exemple en fonction 15 de certaines informations connues a priori sur la répartition des terminaux 70 sur la zone géographique à couvrir. Le serveur 82 procède ensuite à une étape de calcul 322, pour chaque cellule, d'un nombre Klc d'éléments de référence dont la position géographique est située à l'intérieur de ladite cellule. 20 Si, pour une cellule donnée, le nombre Nlc ainsi calculé est supérieur ou égal à un seuil prédéfini N1, alors un point d'intérêt est identifié pour ladite cellule. Le point d'intérêt peut par exemple être défini comme étant le centre de la cellule. Dans une variante, le point d'intérêt peut être défini comme une moyenne des positions géographiques des éléments de référence situés dans 25 la cellule. Lorsque toutes les cellules ont été prises en compte, l'étape d'identification 32 des points d'intérêt prend fin. Un ensemble 35 des points d'intérêt identifiés est alors disponible et mémorisé au sein du serveur 82. La figure 8 est une illustration du premier mode particulier de mise en 30 oeuvre de l'étape d'identification 32 de points d'intérêt décrit en référence à la figure 7. Une zone 50 géographique est quadrillée par 16 cellules identiques en forme de carré. Chaque cellule est identifiée par un couple « L-C », L étant une lettre de A à D, et C étant un chiffre de 1 à 4. Des éléments 51 de référence 22 acquis pendant la phase 20 de calibration sont représentés sous forme de croix. Chaque croix représente la position géographique associée à un élément 51 de référence acquis dans la zone 50 géographique considérée. Dans l'exemple considéré : 5 - Nc = 0 pour les cellules A-1, D-1, C-2 et B-3, on a Nc = 1 pour les cellules C-1, D-2, A-3, C-4 et D-4, - Nc = 2 pour les cellules A-2, B-2, D-3 et B-4, on a Nc= 5 pour la cellule A-4, - Nc = 6 pour la cellule C-3, et - Nc = 7 pour la cellule B-1. 10 Dans l'exemple considéré, le seuil N1 est fixé à 5, et un point d'intérêt 52 est défini comme étant le centre d'une cellule pour laquelle Nc N1. Chaque point d'intérêt 52 ainsi défini est représenté sur la figure 7 par un point noir. Dans l'exemple considéré, trois points d'intérêt 52 ont été identifié. Ils correspondent respectivement au centre des cellules B-1, C-3 et A-4. 15 Ce premier mode particulier de mise en oeuvre de l'étape d'identification 32 des points d'intérêt est plus ou moins précis selon les dimensions des cellules. Il apparaît par exemple à la figure 8 que le centre de la cellule C-3 est relativement éloigné d'une position plus probable d'un point d'intérêt qui serait situé en bas à droite de la cellule, là où se trouvent la 20 majorité des éléments 51 de référence situés dans ladite cellule. La figure 9 représente schématiquement les principales étapes d'un deuxième mode particulier de mise en oeuvre de l'étape d'identification 32 des points d'intérêt. Pour ce deuxième mode particulier de mise en oeuvre de l'étape 25 d'identification 32 des points d'intérêt, des opérations sont effectuées sur un ensemble 27 d'éléments de référence temporaire qui correspond initialement à l'ensemble 26 d'éléments de référence acquis pendant la phase 20 de calibration. Ainsi, l'ensemble 26 d'éléments de référence acquis pendant la phase 20 de calibration n'est pas modifié par ce deuxième mode de mise en 30 oeuvre de l'étape d'identification 32 des points d'intérêt. Dans une première étape, le serveur 82 procède à un calcul 321 des distances séparant deux éléments de référence de l'ensemble 27. 23 Dans une deuxième étape, pour chaque élément de référence de l'ensemble 27 d'éléments de référence temporaire, le serveur 82 procède à un calcul 323 du nombre Nv d'éléments de référence, dit « éléments voisins », dont la distance à l'élément de référence considéré est inférieure ou égale à un 5 seuil prédéfini D1. Vient ensuite une étape de détermination 325 du nombre Nvmax correspondant à la plus grande valeur parmi les nombres Nv précédemment calculés. Si Nvmax est supérieur ou égal à une valeur seuil prédéfinie N2, le 10 serveur 82 procède alors à une étape de définition 327 d'un point d'intérêt à partir de l'élément de référence pour lequel le nombre de voisins est égal à Nvmax, dit « élément d'intérêt », puis à une étape de suppression 329 dudit élément d'intérêt et de ses éléments voisins dans l'ensemble 27 temporaire. Le point d'intérêt est par exemple défini comme étant la moyenne des positions 15 géographiques de l'élément d'intérêt et de ses éléments voisins. Selon un autre exemple, le point d'intérêt peut être défini simplement comme étant la position géographique de l'élément d'intérêt. Une condition d'arrêt est ensuite évaluée. Tant que ladite condition d'arrêt n'est pas atteinte, on réitère depuis l'étape de calcul 323 du nombre de 20 voisins pour les éléments restant dans l'ensemble 27 temporaire. Si le nombre Nvmax calculé à l'étape 325 est inférieur à N2, alors la condition d'arrêt est directement évaluée. La condition d'arrêt est par exemple atteinte si un nombre prédéfini de points d'intérêt ont été identifiés, ou bien si le nombre Nvmax est inférieur à un 25 seuil prédéfini, ou bien s'il n'y a plus d'éléments dans l'ensemble 27 temporaire. Lorsque la condition d'arrêt est atteinte, l'étape d'identification 32 des points d'intérêt prend fin. Un ensemble 35 des points d'intérêt identifiés est alors disponible et mémorisé au sein du serveur 82. La figure 10 est une illustration du deuxième mode particulier de mise 30 en oeuvre de l'étape d'identification 32 de points d'intérêt décrit en référence à la figure 9. La figure 10 représente une zone 50 géographique et des éléments 51 de référence identiques à ceux représentés dans la figure 8. Des cercles 530, 24 531, 532 en pointillé regroupent respectivement un élément d'intérêt 510, 511, 512 et ses éléments voisins. Dans l'exemple considéré, un premier élément 512 d'intérêt a été déterminé avec 5 éléments voisins. Un point d'intérêt 522 a alors été défini comme étant la moyenne des positions géographiques dudit 5 élément d'intérêt 512 et de ses éléments voisins (situés dans le cercle 532). Puis, dans une deuxième itération, un deuxième élément 511 d'intérêt a été déterminé avec 4 éléments voisins. Un point d'intérêt 521 a alors été défini. Enfin, dans une troisième et dernière itération, un troisième point d'intérêt 520 a été défini à partir d'un élément 510 d'intérêt possédant également 4 éléments 10 voisins. Il convient de noter que dans l'exemple considéré, si plusieurs éléments ont la même valeur Nvmax de voisins, le choix de l'un ou de l'autre comme élément d'intérêt est arbitraire. Dans des modes particuliers de mise en oeuvre, l'étape d'identification 32 des points d'intérêt décrite en référence aux figures 7 et 9 peut en outre 15 comporter une étape d'agrégation et/ou de filtrage des points d'intérêt identifiés. Par exemple, si deux points d'intérêt sont situés à une distance inférieure à un seuil prédéfini, alors ils sont agrégés en un unique point d'intérêt. Selon un autre exemple, si pour un point d'intérêt donné, le nombre d'éléments de l'ensemble 26 d'éléments de référence situés à une distance 20 dudit point d'intérêt inférieure à une distance prédéfinie est inférieur à un seuil prédéfini, alors ledit point d'intérêt est supprimé. Dans des modes particuliers de mise en oeuvre, au lieu de calculer un nombre entier d'éléments de référence (par exemple aux étapes 322 et 323), il est possible d'utiliser une valeur représentative d'un nombre d'éléments de 25 référence calculée en fonction d'un facteur de pondération affecté à chaque élément de référence et représentatif d'une durée écoulée entre le moment d'utilisation des liens radio (dont les niveaux de RSSI mesurés servent à définir une signature radio) et le moment où les points d'intérêt sont identifiés. De telles dispositions permettent notamment de favoriser des points d'intérêt 30 récents en accordant plus d'importance (c'est-à-dire en affectant une pondération plus forte) aux éléments de référence récents, c'est-à-dire aux éléments de référence correspondant à des messages émis récemment par des terminaux 72 de calibration. Par exemple, un facteur de pondération peut 25 être défini pour un élément de référence en fonction du temps par : t w(t) = A - exp( ) 2 0-2 fonction dans laquelle : - A et o- sont des paramètres constants, 5 - t est le temps écoulé depuis le moment où le message émis par un terminal 72 de calibration et utilisé pour définir l'élément de référence a été reçu, - expO est la fonction mathématique exponentielle. Avec une telle définition du facteur de pondération, plus le temps 10 écoulé depuis la réception du message qui a servi à définir l'élément de référence est grand, et plus la valeur du facteur de pondération associé sera petite et donc plus l'importance accordée à l'élément de référence sera faible comparativement à d'autres éléments de référence pour lesquels le message associé a été reçu plus récemment. Il est ainsi possible d'accélérer la 15 découverte de nouveaux points d'intérêt et d'ignorer des points d'intérêt désuets. La figure 11 représente schématiquement les principales étapes d'un premier mode particulier de mise en oeuvre de l'étape de modification 34 de l'ensemble 26 d'éléments de référence à partir de l'ensemble 35 des points 20 d'intérêt identifiés. Dans une première étape, le serveur 82 procède, pour chaque élément de l'ensemble 26 d'éléments de référence, à un calcul 340 de la plus petite distance Dr,* entre ledit élément de référence et un point d'intérêt de l'ensemble 35 des points d'intérêt identifiés. 25 Si cette distance Dmin est inférieure ou égale à un seuil D2 prédéfini, alors le serveur 82 procède à une association 342 entre le point d'intérêt correspondant à la distance Dr,* et l'élément de référence. Autrement dit, l'élément de référence est enrichi avec un identifiant ou label unique du point d'intérêt. L'élément de référence ainsi enrichi comporte alors une signature 30 radio, une position géographique, ainsi qu'un label correspondant au point d'intérêt auquel il est associé. Si par contre Dr,* > D2, alors le serveur 82 procède à la suppression 26 344 de l'élément de référence. Lorsque tous les éléments de référence ont été pris en compte, l'étape de modification 34 s'arrête et un ensemble 36 d'éléments de référence modifié est alors disponible. Cet ensemble 36 comporte uniquement les éléments de 5 référence qui ont pu être associé à un point d'intérêt, c'est-à-dire les éléments de référence suffisamment proches d'un point d'intérêt. L'ensemble 36 contient pour chaque point d'intérêt un certain nombre de signatures radio associées. La figure 12 représente schématiquement les principales étapes d'un deuxième mode particulier de mise en oeuvre de l'étape de modification 34 de 10 l'ensemble 26 d'éléments de référence à partir de l'ensemble 35 des points d'intérêt identifiés. Dans une première étape, le serveur 82 procède, pour chaque élément de référence de l'ensemble 26 d'éléments de référence, à un calcul 341, pour chaque point d'intérêt identifié, d'une valeur représentative de la 15 distance qui sépare la position géographique de l'élément de référence de celle du point d'intérêt. Pour un élément de référence donné, il est possible de définir, par exemple selon une loi de distribution de Laplace, une fonction de « probabilité » associant une valeur comprise entre 0 et 1 à chaque point d'intérêt en fonction de la distance séparant le point d'intérêt et l'élément de 20 référence, de telle sorte que plus l'élément de référence est proche d'un point d'intérêt et plus la valeur de probabilité est grande, et de telle sorte que la somme des valeurs de probabilité soit égale à 1. Le serveur 82 procède alors à une association 343 des valeurs de probabilité ainsi calculées à l'élément de référence. Autrement dit, l'élément de 25 référence est enrichi avec un ensemble de valeurs, chaque valeur représentant une probabilité d'être associé à un point d'intérêt. L'élément de référence ainsi enrichi comporte alors une signature radio, une position géographique, ainsi qu'une valeur de probabilité pour chaque point d'intérêt identifié. Lorsque tous les éléments de référence ont été pris en compte, l'étape 30 de modification 34 s'arrête et un ensemble 36 d'éléments de référence modifié est alors disponible. L'ensemble 36 comporte ainsi pour chaque signature radio d'un élément de référence une probabilité d'être associée à un point d'intérêt identifié. 27 Pendant la phase 40 de recherche, l'estimation 44 de la position géographique d'un terminal 70 d'intérêt est par exemple réalisée à l'aide d'un algorithme d'apprentissage automatique de classification conventionnel de type « arbre décisionnel ». Le fait d'avoir enrichi les éléments de référence présents 5 dans l'ensemble 36 permet en effet de transformer un problème généralement abordé sous la forme d'une régression en un problème de classification. En effet, pendant la phase 30 d'affinement, à chaque point d'intérêt est attribué un label, et la signature radio associée à chaque élément de référence est « classée » selon sa probabilité d'appartenir à tel ou tel label. Pendant la phase 10 40 de recherche, à partir de la signature radio d'un terminal 70 d'intérêt, et à l'aide des éléments de référence présents dans l'ensemble 36 modifié, l'algorithme de classification est alors capable de fournir une probabilité pour ladite signature radio d'appartenir à tel ou tel label. La position géographique du terminal 70 d'intérêt est alors généralement estimée comme étant la 15 position géographique du point d'intérêt correspondant au label ayant la plus grande valeur de probabilité. Si la plus grande valeur de probabilité est inférieure à un seuil prédéfini, il est aussi possible de conclure que le terminal 70 d'intérêt n'est pas localisé à un des points d'intérêt identifiés (auquel cas le terminal 70 d'intérêt est probablement en transit entre deux points d'intérêt). 20 Plusieurs simulations montrent que le fait d'avoir enrichi et/ou nettoyé l'ensemble 26 des éléments de référence dans la phase 30 d'affinement apporte un gain de précision significatif à l'algorithme de classification. Les résultats de ces simulations sont illustrés dans les figures 13 et 14. Pour réaliser ces simulations, un ensemble 26 d'éléments de 25 référence est acquis pendant une phase 20 de calibration, et il est utilisé pour chacune des simulations. D'autre part, des éléments, dits « éléments de test », différents des éléments de référence, sont disponibles. Chaque élément de test comporte une signature radio et une position géographique connue. Pour chaque élément de test, une position géographique est estimée et une distance 30 d'erreur correspondant à la distance entre la position géographique estimée et la position géographique connue est calculée. Une courbe 90 à 96 représentant une répartition de la distance d'erreur parmi tous les éléments de test, ou bien parmi un échantillon seulement des éléments de test, est alors 28 tracée. Dans ces simulations, 200.000 éléments de référence et 200.000 éléments de tests sont utilisés. La taille de la zone géographique considérée correspond environ à celle d'un pays grand comme la France. 5 Pour les simulations dont les résultats sont illustrés à la figure 13, un procédé de géolocalisation conventionnel basé sur un algorithme d'apprentissage automatique de régression est utilisé. Une position géographique est estimée par l'algorithme pour chaque élément de test à partir de la signature radio dudit élément de test et de l'ensemble 26 d'éléments de 10 référence acquis pendant la phase 20 de calibration. La courbe 90 représente la répartition de la distance d'erreur pour tous les éléments de test disponibles. La courbe 91 représente quant à elle la répartition de la distance d'erreur pour un échantillon correspondant à des éléments de test 15 suffisamment proches de points d'intérêt identifiés. Dans l'exemple considéré, des points d'intérêt sont identifiés à partir de l'ensemble 26 d'éléments de référence à l'aide du deuxième mode particulier de mise en oeuvre de l'étape d'identification 32 des points d'intérêt, tel que décrit en référence à la figure 9, et un élément de test est considéré suffisamment proche d'un point d'intérêt si 20 la plus petite distance (Dmin) entre la position géographique connue dudit élément de test et un point d'intérêt identifié est inférieure ou égale à un seuil (D2). Pour ces simulations, la valeur de D2 est fixée à 250 mètres. Il apparaît ainsi sur la figure 13 que l'algorithme de régression est plus précis lorsqu'il estime la position géographique d'un terminal 70 d'intérêt situé à 25 proximité d'un des points d'intérêt identifiés. Cela s'explique intuitivement par le fait que l'algorithme d'apprentissage a une meilleure connaissance de la relation existant entre une signature radio et une position géographique pour des éléments proches d'un point d'intérêt car dans l'ensemble 26 la densité géographique d'éléments de référence est plus forte à proximité des points 30 d'intérêt. La précision de l'algorithme de régression est cependant décevante, par exemple la distance d'erreur est inférieure à 0,5 km seulement dans environ 50% des cas. La figure 14 illustre les résultats obtenus avec un procédé 10 de 29 géolocalisation selon l'invention. Un algorithme d'apprentissage automatique de classification est utilisé pour la phase 40 de recherche, et il utilise un ensemble 36 d'éléments de référence modifié par une phase 30 d'affinement. Dans l'exemple considéré, la phase 30 d'affinement comporte une étape 5 d'identification 32 des points d'intérêt selon le mode particulier de mise en oeuvre décrit en référence à la figure 9 et une étape 34 de modification de l'ensemble 26 d'éléments de référence selon le mode particulier de mise en oeuvre décrit en référence à la figure 11. L'algorithme de classification fournit alors, pour un élément de test considéré, une valeur de probabilité pour ledit 10 élément de test d'être associé à un point d'intérêt identifié. La position géographique estimée de l'élément de test est alors définie comme étant la position géographique du point d'intérêt pour lequel ladite valeur de probabilité est la plus grande. Pour ces simulations, la valeur de D1 est également fixée à 250 mètres, et la valeur du seuil N2 est fixée à 300. 15 La courbe 92 représente la répartition de la distance d'erreur pour tous les éléments de test disponibles. Dans ce cas, la distance d'erreur est inférieure à 0,5 km dans environ 70% des cas. La courbe 93 représente la répartition de la distance d'erreur pour l'échantillon d'éléments de test utilisé pour la courbe 91 correspondant à des 20 éléments de test considérés comme suffisamment proches de points d'intérêt identifiés (distance à un point d'intérêt inférieure ou égale au seuil D2). Pour cet échantillon, la distance d'erreur est inférieure à 0,5 km dans environ 90% des cas. La courbe 94 représente la répartition de la distance d'erreur pour un 25 échantillon d'éléments de test pour lesquels la plus grande valeur de probabilité associée à un point d'intérêt est supérieure ou égale à 0,9. Pour cet échantillon, la distance d'erreur est inférieure à 0,5 km dans environ 87% des cas. Cet échantillon regroupe 56% des éléments de test. Dans ce cas, 32% des éléments de test considérés comme proche d'un point d'intérêt (distance 30 inférieure ou égale au seuil D2) ont été « faussement » écartés de l'échantillon. La courbe 95 représente la répartition de la distance d'erreur pour un échantillon d'éléments de test pour lesquels la plus grande valeur de probabilité associée à un point d'intérêt est supérieure ou égale à 0,7. Pour cet 30 échantillon, la distance d'erreur est inférieure à 0,5 km dans environ 82% des cas. Cet échantillon regroupe 76% des éléments de test. Dans ce cas, 15% des éléments de test considérés comme proche d'un point d'intérêt (distance à un point d'intérêt inférieure ou égale au seuil D2) ont été « faussement » 5 écartés de l'échantillon. La courbe 96 représente la répartition de la distance d'erreur pour un échantillon d'éléments de test pour lesquels la plus grande valeur de probabilité associée à un point d'intérêt est supérieure ou égale à 0,4. Pour cet échantillon, la distance d'erreur est inférieure à 0,5 km dans environ 73% des 10 cas. Cet échantillon regroupe 96% des éléments de test. Dans ce cas, 1,5% des éléments de test considérés comme proche d'un point d'intérêt (distance à un point d'intérêt inférieure ou égale au seuil D2) ont été « faussement » écartés de l'échantillon. Dans les différents cas représentés par les courbes 92 à 96, la 15 distance d'erreur est inférieure à 0,2 km dans environ 50% des cas. Il apparaît ainsi que le procédé 10 de géolocalisation selon l'invention dont les résultats sont illustrés par les courbes 92 à 96 de la figure 14 est nettement plus précis qu'un procédé 10 de géolocalisation conventionnel dont les résultats sont illustrés par les courbes 90 et 91 de la figure 13. 20 Les résultats illustrés à la figure 14 montrent qu'avec le procédé 10 de géolocalisation selon l'invention, il est possible d'estimer avec une bonne précision la position géographique d'un terminal 70 d'intérêt situé à proximité d'un point d'intérêt. Il est également possible de déterminer si le terminal 70 d'intérêt a une forte probabilité d'être en transit entre deux points d'intérêt (par 25 exemple si la plus grande valeur de probabilité retournée est inférieure à un seuil prédéfini). La valeur seuil de probabilité permet de déterminer si la position d'un terminal 70 d'intérêt peut être associée à la position d'un point d'intérêt, ou bien si au contraire il faut considérer que ledit terminal 70 d'intérêt est en transit 30 entre des points d'intérêt. Cette valeur seuil de probabilité peut par exemple être déterminée en fonction d'une courbe reliant la qualité de la prédiction et la probabilité de confiance de cette prédiction (par exemple à partir du point d'inflexion ou indice de Youden d'une courbe ROC (acronyme anglais de 31 « Receiver Operating Characteristic »y Selon un autre exemple, elle peut être choisie de manière arbitraire en fonction de contraintes liées au domaine d'application. Par exemple s'il convient d'avoir une très forte précision de géolocalisation d'un terminal 70 d'intérêt, au risque de ne pas pouvoir estimer 5 la position d'un nombre important de terminaux 70 d'intérêt, il faut prendre une valeur seuil de probabilité grande, par exemple la valeur 0,9 comme pour l'exemple illustré avec la courbe 94. Si au contraire il s'agit d'être capable de géolocaliser une grande proportion de terminaux 70 d'intérêt, au risque de réduire la précision de géolocalisation, alors il faut prendre une valeur seuil de 10 probabilité plus faible, par exemple la valeur 0,4 comme pour l'exemple illustré avec la courbe 95. Pour l'exemple considéré, à savoir un procédé de géolocalisation de palettes de marchandises sur une surface géographique importante correspondant par exemple à un pays ou un continent, une valeur seuil de probabilité de l'ordre de 0,7 semble être un bon compromis. 15 La description ci-avant illustre clairement que, par ses différentes caractéristiques et leurs avantages, la présente invention atteint les objectifs fixés et apportent en outre des avantages supplémentaires. En particulier, le procédé 10 de géolocalisation selon l'invention améliore considérablement les performances en termes de précision par 20 rapport à des procédés conventionnels de géolocalisation basés sur des méthodes d'apprentissage. Ce gain de précision est notamment dû à la phase 30 d'affinement qui permet de renforcer une relation entre une signature radio et une position géographique en tirant profit de l'existence de points d'intérêt auxquels la probabilité de présence d'un terminal 70, 72 est plus grande 25 qu'ailleurs. La phase 30 d'affinement permet également de supprimer des éléments disparates qui ne correspondent pas à un point d'intérêt et qui contribuent à générer un « bruit » qui limite les performances de l'algorithme d'apprentissage. En outre, la suppression d'éléments disparates contribue à diminuer le 30 nombre d'éléments de référence utilisés par l'algorithme d'apprentissage, et ainsi à en réduire la complexité, ce qui se traduit en un gain en termes de capacité et de temps de calcul. De plus, la suppression d'éléments disparates peut contribuer à 32 réduire la taille de la signature radio (c'est-à-dire le nombre de mesures RSSI d'une signature radio) en ne considérant que les stations de base qui sont localisées à proximité des points d'intérêt, ce qui contribue encore davantage à réduire la complexité de l'algorithme d'apprentissage. 5 Le fait de connaître la liste des points d'intérêt pour les terminaux 70, 72 du système 60 peut également être intéressant en soi, par exemple pour connaître la liste des entrepôts utilisés dans le transport de marchandises. Un autre avantage repose sur la possibilité d'attribuer un facteur de pondération à chaque élément de référence en fonction du temps écoulé depuis le moment 10 où un message a été émis par un terminal 72 de calibration pour enrichir l'ensemble 26 d'éléments de référence. Cela permet en effet d'accélérer la découverte de nouveaux points d'intérêt et de favoriser l'oubli de points d'intérêt désuets. Aussi, la mise en place d'un tel procédé 10 de géolocalisation au sein 15 d'un système 60 de communication sans fil comme celui qui a été décrit à titre d'exemple est particulièrement peu onéreuse. En effet, elle ne nécessite pas de matériel spécifique additionnel au niveau de la majorité des terminaux 70 ou des stations de base 81, et grâce à l'existence de certains terminaux spéciaux 72 équipés d'un système de positionnement dans le système 60 elle permet 20 d'automatiser entièrement la phase 20 de calibration. De manière générale, il est à noter que les modes de mise en oeuvre et de réalisation considérés ci-dessus ont été décrits à titre d'exemples non limitatifs, et que d'autres variantes sont par conséquent envisageables. Notamment, les modes de mise en oeuvre de l'étape d'identification 32 25 des points d'intérêt et de modification 34 de l'ensemble 26 d'éléments de référence ont été donnés à titre d'exemple, et d'autres méthodes pourraient être envisagées comme des variantes de l'invention. L'invention a été décrite pour un système 60 de communication radio sans fil à bande ultra étroite adapté à des applications du type loT, mais rien 30 n'empêche de mettre en oeuvre le procédé 10 de géolocalisation selon l'invention pour d'autres types de système de communication comme par exemple des réseaux cellulaires GSM, UMTS, LTE, Wi-Fi, etc. Les mesures représentatives du niveau de qualité de lien radio 33 utilisées pour former une signature radio peuvent être effectuées par des stations de base, sur un lien radio montant, mais aussi par des terminaux, sur un lien radio descendant. Les mesures sont par exemple transmises à un serveur qui réalise alors les autres étapes du procédé de géolocalisation. Selon 5 d'autres exemples, les mesures représentatives du niveau de qualité de lien radio pourraient être effectuées par le serveur lui-même, à partir d'informations reçues des terminaux et/ou des stations de base. Comme indiqué précédemment, la valeur représentative de la qualité de lien radio peut être différente du niveau de RSSI. Cela peut être par 10 exemple, un niveau d'atténuation du signal, un rapport signal sur bruit, ou bien un autre indicateur de qualité de canal radio. Dans l'exemple considéré, un élément de référence correspond à la position géographique obtenue par le récepteur GPS d'un terminal 72 de calibration associée à la signature radio dudit terminal 72 de calibration à ladite 15 position géographique. Dans des modes particuliers de réalisation, l'ensemble 26 d'éléments de référence peut en outre être continuellement enrichi avec des éléments de référence correspondant chacun à une position géographique d'un terminal 70 d'intérêt estimée pendant la phase de recherche et associée à la signature radio dudit terminal 70 d'intérêt situé à cette position. 20 D'autres méthodes peuvent néanmoins être envisagées pour définir l'ensemble 26 d'éléments de référence. Par exemple, la position géographique d'un terminal (72) de calibration peut être estimée pendant la phase de calibration par une méthode conventionnelle de géolocalisation basée sur un algorithme d'apprentissage automatique par régression. Selon encore un autre 25 exemple, la position géographique d'un terminal (72) de calibration peut être estimée pendant la phase de calibration à l'aide d'une autre méthode de géolocalisation, par exemple une méthode de type TDOA (acronyme anglais de « Time Difference Of Arrivai ») basée sur des mesures de temps de propagation d'un signal échangé entre un terminal et différentes stations de 30 base. Pour la phase de recherche, différents algorithmes d'apprentissage automatique de classification peuvent être utilisés pour estimer la position géographique d'un terminal 70 d'intérêt à partir de sa signature radio et de 34 l'ensemble 36 d'éléments de référence modifié. Le fait d'utiliser l'un ou l'autre de ces algorithmes ne constitue qu'une variante de mise en oeuvre de l'invention. L'invention a été décrite en considérant le domaine du transport de 5 marchandises. Rien n'exclut cependant de considérer d'autres domaines d'application pour le procédé 10 de géolocalisation selon l'invention, en particulier des domaines pour lesquels les terminaux à géolocaliser sont susceptibles d'occuper préférentiellement des positions géographiques particulières ou de passer par des positions géographiques particulières 10 auxquelles un grand nombre de terminaux 70 du système 60 ont tendance à passer (géolocalisation de bagages dans des aéroports, géolocalisation de colis, etc.). Dans l'exemple décrit, les points d'intérêt sont identifiés dynamiquement à partir des éléments de référence fournis par des terminaux 15 72 de calibration. Rien n'empêche cependant que dans des modes particuliers de mise en oeuvre, tout ou partie des points d'intérêts soient connus a priori et que des signatures radio associées à ces points d'intérêt soient déterminées à l'aide de dispositifs 71 de calibration pouvant être déplacés intentionnellement sur lesdits points d'intérêt REVENDICATIONS1. Procédé (10) de géolocalisation d'un terminal, dit « terminal (70) d'intérêt », d'un système (60) de communication sans fil, comportant : une détermination, à partir de messages émis par d'autres terminaux, dits « terminaux (72) de calibration », et reçus par des stations de base (81) dudit système (60) de communication sans fil, d'un ensemble (26) d'éléments de référence, chaque élément de référence comportant une signature radio associée à une position géographique d'un terminal de calibration lors de l'émission d'un message, ladite signature radio correspondant à un ensemble de valeurs représentatives de la qualité de liens radio existant entre ledit terminal (72) de calibration et des stations de base (81), une détermination (42), à partir d'un message émis par ledit terminal (70) d'intérêt situé à une position géographique à estimer, d'une signature radio pour ledit terminal (70) d'intérêt, ledit procédé (10) de géolocalisation étant caractérisé en ce qu'il comporte en outre : une identification (32) de positions géographiques particulières, dites « points d'intérêt », à partir des éléments de référence de l'ensemble (26) d'éléments de référence, une modification (34) de l'ensemble (26) d'éléments de référence à partir des points d'intérêt identifiés, une estimation (44) de la position géographique dudit terminal d'intérêt à partir de la signature radio dudit terminal (70) d'intérêt et de l'ensemble (36) d'éléments de référence modifié. 2. Procédé (10) de géolocalisation selon la revendication 1 dans lequel la position géographique du terminal (70) d'intérêt est estimée sous la forme d'une valeur représentative d'une probabilité pour ledit terminal (70) d'intérêt d'être situé en un point d'intérêt identifié. 3. Procédé (10) de géolocalisation selon l'une des revendications 1 à 2 dans lequel l'identification (32) des points d'intérêt comporte : 36 un découpage (320) d'une zone (50) géographique à couvrir en plusieurs cellules de tailles prédéfinies, un calcul (322) pour chaque cellule d'une valeur représentative d'un nombre d'éléments de référence situés dans ladite cellule, 5 si la valeur calculée satisfait un critère prédéfini, une définition (324) d'un point d'intérêt en fonction de ladite cellule et/ou des positions géographiques des éléments de référence situés dans ladite cellule. 4. Procédé (10) de géolocalisation selon l'une des revendications 1 à 2 dans lequel l'identification (32) des points d'intérêt comporte une étape de calcul 10 (321), pour chaque élément de référence, de distances séparant ledit élément de référence et chaque autre élément de référence, ainsi que les étapes suivantes effectuées itérativement sur un ensemble (27) d'éléments de référence temporaire correspondant initialement à l'ensemble (26) d'éléments de référence, jusqu'à ce qu'une condition d'arrêt soit atteinte : 15 un calcul (323) pour chaque élément de référence de l'ensemble (27) d'éléments de référence temporaire, d'une valeur représentative d'un nombre d'éléments de référence situés à une distance inférieure à un seuil prédéfini, dits « éléments voisins », si la plus grande valeur calculée, pour un élément dit « élément 20 d'intérêt », est supérieure à un seuil prédéfini : o une définition (327) d'un point d'intérêt en fonction des positions géographiques de l'élément d'intérêt et/ou de ses éléments voisins, o une suppression (329) de l'élément d'intérêt et de ses éléments voisins de l'ensemble (27) d'éléments de référence temporaire. 25 5. Procédé (10) de géolocalisation selon l'une des revendications 1 à 4 comportant en outre une étape de filtrage des points d'intérêt identifiés en fonction de leurs positions géographiques et/ou d'une valeur représentative d'un nombre d'éléments de référence de l'ensemble (26) d'éléments de référence situés à une distance inférieure à un seuil prédéfini. 30 6. Procédé (10) de géolocalisation selon l'une des revendications 3 à 5 dans lequel la valeur représentative d'un nombre d'éléments de référence est calculée en fonction d'un facteur de pondération affecté à chaque élément de référence et représentatif d'une durée écoulée entre le moment 37 d'émission du message par le terminal (72) de calibration et le moment où les points d'intérêt sont identifiés. 7. Procédé (10) de géolocalisation selon l'une des revendications 1 à 6 dans lequel la modification (34) de l'ensemble (26) d'éléments de référence 5 comporte les étapes suivantes pour chaque élément de référence de l'ensemble (26) d'éléments de référence : un calcul (340) de la plus petite distance séparant ledit élément de référence et un des points d'intérêt identifiés, si la plus petite distance calculée est inférieure à un seuil prédéfini : 10 o une association (342) audit élément de référence du point d'intérêt correspondant, - sinon : o une suppression (344) dudit élément de référence de l'ensemble (26) d'éléments de référence. 15 8. Procédé (10) de géolocalisation selon l'une des revendications 1 à 6 dans lequel la modification (34) de l'ensemble (26) d'éléments de référence comporte les étapes suivantes pour chaque élément de référence de l'ensemble (26) d'éléments de référence : un calcul (341), pour chaque point d'intérêt identifié, d'une valeur 20 représentative de la distance qui sépare ledit élément de référence dudit point d'intérêt, une association (343) audit élément de référence des valeurs ainsi calculées. 9. Procédé (10) de géolocalisation selon l'une des revendications 1 à 8 dans 25 lequel un lien radio est un lien montant vers une station de base (81) du système (60) de communication sans fil. 10. Procédé (10) de géolocalisation selon la revendication 9 dans lequel une détermination (22, 42) de la signature radio d'un terminal (70, 72) comporte 30 une émission (220, 420) par le terminal (70, 72) d'un message à destination des stations de base (81), une mesure (221, 421), pour chaque station de base (81), d'une valeur représentative de la qualité du lien radio existant entre ledit terminal (70, 38 72) et la station de base (81) considérée à partir du message reçu dudit terminal (70, 72), une formation (222, 422), par un serveur (82) relié aux stations de base (81) de la signature radio du terminal (70, 72) à partir des valeurs 5 mesurées. 11. Procédé (10) de géolocalisation selon la revendication 10 dans lequel un terminal (72) de calibration est équipé d'un système de positionnement et la détermination (22) de la signature radio dudit terminal (72) de calibration comporte en outre : 10 une inclusion, dans le message émis par ledit terminal (72) de calibration, de la position géographique courante dudit terminal (72) de calibration mesurée par le système de positionnement, une extraction (223), par le serveur (82), de la position géographique contenue dans ledit message. 15 12. Procédé (10) de géolocalisation selon l'une des revendications 1 à 11 dans lequel la valeur représentative de la qualité du lien radio existant entre un terminal (70, 72) et une station de base (81) est un niveau de puissance reçue d'un signal radio échangé entre ladite station de base (81) et ledit terminal (70, 72). 20 13. Procédé (10) de géolocalisation selon l'une des revendications 1 à 12 dans lequel le lien radio est un canal de communication à bande ultra étroite. 14. Procédé (10) de géolocalisation selon l'une des revendications 1 à 13 dans lequel l'estimation (44) de la position géographique d'un terminal (70) d'intérêt à partir de la signature radio dudit terminal (70) d'intérêt et de 25 l'ensemble (36) d'éléments de référence modifié est réalisée à l'aide d'un algorithme d'apprentissage automatique de classification basé sur une technique d'arbres décisionnels. 15.Système (60) de communication sans fil comportant des terminaux (70, 72) et un réseau d'accès (80) comportant des stations de base (81) et un 30 serveur (82) relié auxdites stations de base (81), ledit système (60) étant caractérisé en ce qu'il comporte une base de données mémorisant un ensemble (26) d'éléments de référence formés à partir de messages émis par des terminaux, dits « terminaux (72) de calibration », et reçus par des 39 stations de base (81), chaque élément de référence comportant une signature radio associée à une position géographique d'un terminal de calibration lors de l'émission d'un message, ladite signature radio comprenant des valeurs représentatives de la qualité de liens radio existant 5 entre ledit terminal (72) de calibration et des stations de base (81), et en ce que le réseau d'accès (80) est configuré pour : mesurer des valeurs représentatives de la qualité de liens radio existant entre des stations de base (81) et un terminal (70) dont la position géographique doit être estimée, dit « terminal (70) d'intérêt », 10 former une signature radio correspondant à l'ensemble des valeurs mesurées pour ledit terminal (70) d'intérêt, identifier des positions géographiques particulières, dites « points d'intérêt », à partir de l'ensemble (26) d'éléments de référence, modifier l'ensemble (26) d'éléments de référence à partir des points 15 d'intérêt identifiés, estimer la position géographique du terminal (70) d'intérêt à partir de la signature radio du terminal (70) d'intérêt et de l'ensemble (36) d'éléments de référence modifié. 20