FR 3136285 A1 20231208 FR 2205306 A 20220602 Procédé de localisation d'une unité de navigation Ce procédé comporte : - l'établissement (116) de biais de mesure pour un instant t k qui permettent de diminuer l'erreur entachant des mesures brutes de pseudo-distances qui séparent l'unité de navigation d’au moins trois balises, et - l'estimation (122) d'une position externe de l'unité de navigation à l'instant t k à partir : - de positions de balises, - de mesures corrigées des pseudo-distances calculées pour l'instant t k à l’aide des biais de mesure établis, et - de la précédente position externe de l'unité de navigation estimée pour l'instant précédent t k-1 . L’établissement des biais de mesures est réalisé sans prendre en compte la position externe de l'unité de navigation. Fig. 4 Procédé de localisation d'une unité de navigation L’invention concerne un procédé de localisation d'une unité de navigation. L’invention concerne également un support d’enregistrement d’informations et une unité de navigation pour la mise en œuvre de ce procédé. Les procédés de localisation permettent de localiser une unité de navigation dans un environnement comportant des balises qui émettent des messages. Par exemple, ces balises sont des satellites en orbite autour de la terre. Lorsque les balises sont les satellites, l'unité de navigation est connue sous le terme de segment utilisateur du système GNSS (« Global Navigation Satellite System »). Le GNSS américain est connu sous l’acronyme GPS ("Global Positioning System"). Typiquement, les procédés connus de localisation comportent les étapes suivantes : a) la réception des messages émis par au moins trois balises, puis b) l'extraction de mesures brutes de pseudo-distances qui séparent l'unité de navigation de chacune des balises à partir des instants de réception de ces messages mesurés par l'unité de navigation et de la vitesse de propagation de ces messages dans l'espace, puis c) l'estimation d'une position de l'unité de navigation à partir des mesures brutes extraites. Pour améliorer la précision de la position estimée, il est connu de corriger les mesures brutes de pseudo-distance avant qu’elles soient utilisées pour estimer la position de l’unité de navigation. En effet, les signaux émis par les balises peuvent être fortement bruités. Pour corriger les mesures brutes, les documents suivants ont proposés des procédés qui utilisent la précédente position estimée pour l’unité de navigation : - la demande de brevet FR3066027A1, et - A. Giremus, J.-Y. Tourneret et V. Calmettes : “A particle filtering approach for joint detection/ estimation of multipath effects on GNSS measurements”, IEEE Trans. Signal Processing, 2007. Ces procédés sont avantageux car ils permettent d’améliorer la précision de la localisation par rapport aux procédés de localisation qui ne corrige pas les mesures brutes des pseudo-distances en fonction de la précédente position de l’unité de navigation. Toutefois, en pratique, ces procédés sont peu robustes et risquent de diverger au fil du temps. En effet, une erreur dans la correction des mesures brutes des pseudo-distances introduit une erreur dans la position estimée. Or, lors de l’itération suivante du procédé, c’est la position estimée avec une erreur qui est utilisée pour corriger les nouvelles mesures brutes des pseudo-distances. L'erreur commise sur l'estimation de la position peut alors accroître, lors de l'itération suivante, l’erreur commise lors de la correction des mesures brutes des pseudo-distances. Ainsi, il peut apparaître une rétroaction qui amplifie l'erreur au lieu de la compenser. L’invention vise à remédier à cet inconvénient en proposant un procédé de localisation utilisant la position de l’unité de navigation pour corriger les mesures brutes des pseudo-distances tout en étant plus robuste. Elle a donc pour objet un procédé de localisation d'une unité de navigation dans un environnement comportant des balises qui émettent des messages, ce procédé comportant l'exécution des étapes suivantes pour chaque instant t k d'une suite temporelle d'instants {t 0 ; ...; t k-1 ; t k ; ...} : a) la réception des messages émis par au moins trois balises, puis b) l'extraction de mesures brutes de pseudo-distances qui séparent l'unité de navigation de chacune de ces au moins trois balises à partir des instants de réception de ces messages mesurés par l'unité de navigation et de la vitesse de propagation de ces messages dans l'espace, c) l'estimation d'une position interne de l'unité de navigation à l'instant t k à partir des mesures brutes extraites pour cette instant t k et des positions des balises, d) pour plusieurs balises, l'établissement d'un biais de mesure pour l'instant t k qui, lorsqu'il est retranché à la mesure brute extraite à partir de l'instant de réception du message émis par cette balise, permet de diminuer l'erreur entachant cette mesure brute, ce biais étant établi à partir de l'écart entre : - une pseudo-distance calculée qui sépare cette balise de l'unité de navigation, cette pseudo-distance calculée étant calculée à partir d'une position de cette balise et de la position interne de l'unité de navigation estimée pour l'instant t k , et - une pseudo-distance mesurée qui sépare cette balise de l'unité de navigation, cette pseudo-distance mesurée étant calculée à partir de la mesure brute de la pseudo-distance extraite pour cette balise pour l'instant t k et sans prendre en compte la position interne de l'unité de navigation, e) pour chaque balise pour laquelle un biais de mesure a été établi pour l'instant t k , le calcul d'une mesure corrigée de la pseudo-distance à l'instant t k en ajoutant le biais établi pour cette balise à une mesure brute de la pseudo-distance pour cette balise extraite à l'instant t k , dans lequel le procédé comporte également l'étape suivante : f) l'estimation d'une position externe de l'unité de navigation à l'instant t k à partir : - des positions des balises pour lesquelles une mesure corrigée de la pseudo-distance à l'instant t k a été calculée, - des mesures corrigées des pseudo-distances calculées pour l'instant t k , et - de la précédente position externe de l'unité de navigation estimée pour l'instant précédent t k-1 , et dans lequel les étapes b), c) et d) sont réalisées sans prendre en compte la position externe de l'unité de navigation. Les modes de réalisation de ce procédé peuvent comprendre une ou plusieurs des caractéristiques suivantes : 1) - le procédé comporte, après l'étape b) et avant l'étape c), la sélection, en fonction de la position externe de l'unité de navigation à l'instant précédent t k-1 , d'un nombre limité de mesures brutes extraites parmi l'ensemble des mesures brutes extraites lors de l'étape b) pour l'instant t k , puis - lors des étapes c) et d) seuls les mesures brutes ainsi sélectionnées sont utilisées et les mesures brutes non-sélectionnées sont ignorées. 2) chaque biais de mesure pour l'instant t k est également établi en minimisant l'écart entre ce biais de mesure pour l'instant t k et le biais de mesure établi pour l'instant précédent t k-1 . 3) chaque balise est un satellite en orbite autour de la terre. 4) - après l'étape b) et avant l'étape c), le procédé comporte une étape de pré-correction des mesures brutes extraites, cette étape de pré-correction comportant pour chaque balise ayant émis un message à partir duquel est extrait une pseudo-distance : - l'établissement d'une erreur causée par la propagation du message de cette balise à travers la ionosphère ou la troposphère en utilisant pour cela un modèle prédéterminé paramétré par la position interne de l'unité de navigation estimée pour l'instant précédent t k-1 , puis - la soustraction de l'erreur établie à la mesure brute extraite pour obtenir une mesure brute pré-corrigée, puis - lors des étapes c) et d), ce sont les mesures brutes pré-corrigées qui sont utilisées. L’invention a également pour objet un support d’enregistrement d’informations, lisible par un microprocesseur, comportant des instructions exécutables par ce microprocesseur, dans lequel ce support comporte des instructions non transitoires pour l’exécution du procédé de localisation ci-dessus, lorsque ces instructions sont exécutées par le microprocesseur. L’invention a également pour objet une unité de navigation apte à se déplacer dans un environnement comportant des balises qui émettent des messages, cette unité de navigation comportant un microprocesseur configuré pour exécuter les étapes suivantes pour chaque instant t k d'une suite temporelle d'instants {t 0 ; ...; t k-1 ; t k ; ...} : a) la réception des messages émis par au moins trois balises, puis b) l'extraction de mesures brutes de pseudo-distances qui séparent l'unité de navigation de chacune de ces au moins trois balises à partir des instants de réception de ces messages mesurés par l'unité de navigation et de la vitesse de propagation de ces messages dans l'espace, c) l'estimation d'une position interne de l'unité de navigation à l'instant t k à partir des mesures brutes extraites pour cette instant t k et des positions des balises, d) pour plusieurs balises, l'établissement d'un biais de mesure pour l'instant t k qui, lorsqu'il est retranché à la mesure brute extraite à partir de l'instant de réception du message émis par cette balise, permet de diminuer l'erreur entachant cette mesure brute, ce biais étant établi à partir de l'écart entre : - une pseudo-distance calculée qui sépare cette balise de l'unité de navigation, cette pseudo-distance calculée étant calculée à partir d'une position de cette balise et de la position interne de l'unité de navigation estimée pour l'instant t k , et - une pseudo-distance mesurée qui sépare cette balise de l'unité de navigation, cette pseudo-distance mesurée étant calculée à partir de la mesure brute de la pseudo-distance extraite pour cette balise pour l'instant t k et sans prendre en compte la position interne de l'unité de navigation, e) pour chaque balise pour laquelle un biais de mesure a été établi pour l'instant t k , le calcul d'une mesure corrigée de la pseudo-distance à l'instant t k en ajoutant le biais établi pour cette balise à une mesure brute de la pseudo-distance pour cette balise extraite à l'instant t k , dans lequel le microprocesseur est également configuré pour exécuter l'étape suivante : f) l'estimation d'une position externe de l'unité de navigation à l'instant t k à partir : - des positions des balises pour lesquelles une mesure corrigée de la pseudo-distance à l'instant t k a été calculée, - des mesures corrigées des pseudo-distances calculées pour l'instant t k , et - de la précédente position externe de l'unité de navigation estimée pour l'instant précédent t k-1 , et dans lequel les étapes b), c) et d) sont réalisées sans prendre en compte la position externe de l'unité de navigation. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif et faite en se référant aux dessins sur lesquels : - la est une illustration schématique d'un système de localisation d'un véhicule comportant une unité de navigation ; - la est une illustration schématique de l’unité de navigation du système de la , - la est une illustration schématique de différents modules logiciels mis en œuvre dans l’unité de navigation de la ; - la est un organigramme d'un procédé de localisation à l'aide de l’unité de navigation de la ; - la est un organigramme d'un procédé de sélection de mesures mis en œuvre dans le procédé de la . Dans ces figures, les mêmes références sont utilisées pour désigner les mêmes éléments. Dans la suite de cette description, les caractéristiques et fonctions bien connues de l'homme du métier ne sont pas décrites en détail. Dans cette description, des exemples détaillés de modes de réalisation sont d'abord décrits dans le chapitre I en référence aux figures. Ensuite, dans le chapitre II suivant, des variantes de ces modes de réalisation sont introduites. Enfin, les avantages des différents modes de réalisation sont présentés dans un chapitre III. Chapitre I : Exemples de modes de réalisation Par la suite, un mode de réalisation détaillé de l'invention est décrit dans le contexte particulier de la localisation d'un véhicule à la surface de la terre. Dans ce contexte particulier, les balises qui émettent des messages à partir desquels le véhicule peut être localisé sont des satellites en orbite autour de la terre. La représente un véhicule 2 capable de se déplacer à la surface de la terre. A cet effet, par exemple, le véhicule 2 comporte des roues ou des chenilles. Le véhicule 2 est aussi équipé de moyens 4 de propulsion tel qu'un moteur qui entraîne les roues ou les chenilles. Le véhicule 2 est équipé d'un système 6 de géolocalisation de ce véhicule. Ce système 6 est apte à déterminer la position, l'orientation et la vitesse du véhicule 2 dans un repère terrestre R T . Ici, le repère terrestre R T est fixé sans aucun degré de liberté à la terre. Le repère R T comporte trois axes typiquement orthogonaux entre eux. Un repère mobile R b est également fixé sans aucun degré de liberté au véhicule 2. Ce repère R b comporte trois axes orthogonaux entre eux notés respectivement x b , y b et z b . Ici, la position du véhicule 2 dans le repère R T est exprimée par la latitude L, la longitude λ et l'altitude h de l'origine du repère R b . L'orientation du véhicule 2 est exprimée par l'angle ψ de lacet (« yaw angle » en anglais), l'angle θ de tangage (« pitch angle » en anglais) et l'angle φ de roulis (« roll angle » en anglais) du repère R b par rapport au repère R T . La position, l'orientation et la vitesse déterminées par le système 6 sont délivrées sur une sortie 7. Typiquement, le véhicule 2 comporte un poste 8 de pilotage pour guider ou assister au guidage du véhicule 2 vers une destination prédéfinie. Le poste 8 est raccordé à la sortie 7. Le poste 8 peut être un poste de pilotage manuel et/ou automatique. Dans le cas d'un poste de pilotage manuel, la position, l'orientation et la vitesse déterminées sont transmises à une interface homme-machine pour assister un être humain dans le pilotage des moyens 4 de propulsion. Dans le cas d'un poste de pilotage automatique, la position, l'orientation et la vitesse déterminées sont automatiquement converties en commandes de pilotage des moyens 4 de propulsion, puis transmises automatiquement à ces moyens 4 de propulsion. Le système 6 comporte une unité 10 de navigation et une unité 12 de mesure inertielle. L'unité 10 de navigation est capable de déterminer sa position dans le repère R T à partir de messages émis par les satellites. Par exemple, ici, les satellites sont les satellites la constellation GPS ("Global Positioning System"). Dans ce contexte, l'unité 10 est connue sous l'acronyme GNSS (« Global Navigation Satellite System ») et correspond au segment utilisateur de ce système GNSS. A partir des signaux satellitaires qu'elle reçoit, l'unité 10 génère des signaux représentatifs de la position et de la vitesse du véhicule dans le repère R T . L'unité 12 est connue sous l'acronyme IMU (« Inertial Measurement Unit »). L’unité 12 comporte notamment un accéléromètre triaxe 14 et un gyromètre triaxe 16. Grâce à ces capteurs, l'unité 12 est capable de mesurer la variation de l'orientation du véhicule 2. Pour déterminer la position, l'orientation et la vitesse du véhicule 2 à partir des mesures des unités 10 et 12, le système 6 comporte un calculateur électronique programmable 20. Ce calculateur 20 est apte à acquérir les mesures des unités 10 et 12 et, à partir de ces mesures, à déterminer la position, l'orientation et la vitesse du véhicule 2 dans le repère R T . Le calculateur 20 comporte un microprocesseur 22 et une mémoire 24. La mémoire 24 comporte notamment les instructions d'un module logiciel 26 apte à déterminer la position, l'orientation et la vitesse du véhicule 2 à partir des mesures des unités 10 et 12 lorsqu'il est exécuté par le microprocesseur 22. Le module 26 implémente notamment un algorithme de fusion qui établit, à partir d'une précédente estimation de la position, de l'orientation et de la vitesse du véhicule 2 et de nouvelles mesures des unités 10 et 12 acquises depuis cette précédente estimation, une nouvelle estimation de la position, de l'orientation et de la vitesse du véhicule 2. L'algorithme de fusion établit également des marges d'erreur sur chaque nouvelle estimation. De tels algorithmes de fusion sont bien connus de l'homme du métier. Par exemple, une présentation de l'état de l'art en la matière peut être trouvée dans la thèse suivante : S. Godha, “ Performance Evaluation of Low Cost MEMS-Based IMU Integrated With GPS for Land Vehicle Navigation Application ”, PhD report, 2006. Par la suite, cette thèse est désignée par l'expression "Godha2006".. Typiquement, cet algorithme de fusion implémente un ou plusieurs filtres de Kalman. Ici, le module 26 implémente une architecture connue sous le terme de "en boucle fermée" ("closed loop integration scheme" ou "closed loop approach" en anglais). La représente plus en détail les principaux composants de l'unité 10. L'unité 10 comporte successivement : - une antenne 30 pour recevoir les signaux GNSS émis par les satellites, - un étage 32 d’entrée, connu sous le terme anglais de "RF front end", raccordé à l'antenne 30 et qui assure l’amplification, le filtrage, la transposition en fréquence ainsi que la numérisation des signaux GNSS émis, - un étage 34 de traitement en bande de base, connu sous le terme anglais de "Baseband processor" qui reçoit les signaux traités par l'étage 32 d'entrée, et - un processeur électronique 36 de navigation, connu sous le terme anglais de "Navigation processor" qui traite les données transmises par l'étage 34 pour localiser l'unité 10. L'étage 34 acquière les signaux traités par l'étage 32, les démodule et extrait des signaux démodulés les données suivantes : - des mesures brutes Mb relatives à ces signaux, et - des messages émis Mn par les satellites. Les messages émis par les satellites sont connus sous le terme "Messages GNSS". Typiquement, les mesures brutes Mb contiennent pour chaque satellite ayant émis un signal GNSS reçu par l'unité 10, au moins une mesure brute Mb k,i de la pseudo-distance qui sépare l'unité 10 de ce satellite. Par exemple, cette mesure brute Mb k,i de pseudo-distance est calculée ici à l'aide de la relation suivante : Mb k,i = (tr k,i -te k,i ) x c, où: - l'indice i est un identifiant du satellite, - l'indice k est un identifiant d'un instant t k défini plus loin, - Mb k,i est la mesure brute de la pseudo-distance entre l'unité 10 et le satellite identifié par l'indice i, - tr k,i est l'instant de réception du message mesuré par l'étage 34, - te k,i est l'instant d'émission du message par le satellite i, cette instant te k,i étant ici contenu dans le message émis par le satellite, - c est la célérité de la lumière, et - "x" est le symbole qui représente l’opération de multiplication. Éventuellement, les mesures brutes peuvent contenir d'autres mesures pour chaque satellite comme : - une mesure Doppler ou décalages de fréquence, - une mesure de phase de l'onde porteuse du signal GNSS émis, et - l'instant de réception du signal GNSS émis. Ces autres mesures peuvent être utilisées pour estimer, en plus de la position de l'unité 10, d'autres grandeurs physiques comme la vitesse ou l'erreur d'horloge de l'unité 10. Par la suite, lorsque d'autres grandeurs physiques que la position de l'unité 10 sont estimées, il est considéré que ces autres grandeurs physiques sont estimées de façon conventionnelle. Ainsi, par la suite, seul le traitement des mesures brutes de pseudo-distances est décrit en détail. Les messages Mn contiennent pour chaque satellite ayant émis un signal GNSS reçu par l'unité 10 : des éphémérides ou des paramètres permettant d’estimer les positions et vitesses du satellite. Les messages Mn contiennent aussi éventuellement pour chaque satellite d'autres données comme : - des paramètres de corrections d’erreurs d’horloge du satellite, - des éléments caractérisant l’état du satellite et l’état de sa mesure, - des paramètres d’un modèle d’estimation des délais de propagation à travers l'ionosphère et la troposphère. Le processeur 36 comporte typiquement un microprocesseur 38 et un mémoire 40. La mémoire 40 contient les instructions d'un logiciel 42 de localisation qui localise l'unité 10 à partir des signaux émis par les satellites lorsqu'il est exécuté par le microprocesseur 38. Dans ce mode de réalisation, localiser l'unité 10 signifie calculer la position, la vitesse et l'erreur d'horloge de l'unité 10, c'est-à-dire calculer une solution, connue sous l'acronyme PVT (Position, vitesse et temps), pour l'unité 10. Plus précisément, le logiciel 42 met à jour la solution PVT de l'unité 10 à chaque instant t k d'une suite temporelle d'instants {t 0 ; ...; t k-1 ; t k ; ...}. Par exemple, cette suite temporelle est un suite périodique d'instants t k et la période entre deux instants t k successifs est égale à la période d’échantillonnage des signaux GNSS. Par la suite, la solution PVT délivrée par l'unité 10 au calculateur 20 à l'instant t k est notée PVT k . Dans ce texte, l'indice k est également utilisé pour identifier les nouvelles mesures, les nouveaux messages et les nouvelles grandeurs physiques utilisées pour estimer la solution PVT k . Ainsi, les nouvelles mesures brutes Mb et les nouveaux messages Mn obtenus à partir du signal GNSS d'un satellite et utilisés pour estimer la solution PVT k sont notés, respectivement, Mb k,i et Mn k,i , où l'indice i est un identifiant de ce satellite. Pour estimer la solution PVT k , les mesures Mb k,i et les messages Mn k,i d'au moins trois et, dans ce mode de réalisation, d'au moins quatre satellites sont utilisés. En pratique le nombre de satellites dont les messages GNSS sont reçus par l'unité 10 est souvent supérieur à cinq ou huit. Ainsi, l'indice i est ici un nombre entier qui varie de 1 à plus de 4, 5 ou 8. La représente l'architecture du logiciel 42 de localisation. Le logiciel 42 comporte un module 50 d'estimation de la position des satellites à partir des mesures Mb k,i et des données contenues dans les messages Mn k,i . Plus précisément, ici, pour chaque instant t k , le module 50 estime les nouvelles positions, les nouvelles vitesses et les nouvelles erreurs d’horloges des satellites à partir des mesures Mb k,i et des messages Mn k,i reçus depuis l'instant précédent t k-1 . Autrement dit, le module 50 estime la solution PVT de chaque satellite i pour l'instant t k . Cette solution PVT du satellite i pour l'instant t k est notée PVTs k,i par la suite. On note PVTs k , le vecteur dont la i-ième coordonnée est égale à la solution PVTs k,i . De nombreuses méthodes sont connues pour estimer la solution PVTs k,i à partir des mesures Mb k,i et des messages Mn k,i . Par exemple, dans un premier temps, les instants d’émission des signaux GNSS sont estimés à partir : - de la mesure des instants tr k de réception de ces signaux GNSS par l'unité 10, et - de la mesure des pseudo-distances entre les satellites et l’antenne 30. Dans un second temps, les positions, vitesses et erreurs d’horloges des satellites sont estimées en exécutant un algorithme qui utilise des paramètres fournis dans les messages Mn k,i . Par exemple, cet algorithme est l'un de ceux proposés dans les articles suivants : - Global navigation satellite system Glonass ICD Navigational radiosignal in bands L1, L2 , Edition 5.1, Russian Institute of Space Device Engineering, 2008, pages 41-49 - Navstar GPS Space Segment/Navigation User Interfaces , IS-GPS-200, Revision D, ARINC, March 2006, pages 83-98 Ensuite le logiciel 42 comporte : - un processus 52 qui établit des biais c2 k,i de mesure qui sont pratiquement indépendant de la précédente solution PVT k-1 , et - un processus 54 qui estime la solution PVT k en utilisant les biais de mesures c2 k,i pour corriger les mesures brutes Mb k,i et ainsi améliorer la précision de la solution PVT k . Un biais de mesure est une valeur qui, lorsqu'elle est retranchée à une mesure brute permet de diminuer l'erreur entachant cette mesure brute. Le processus 52 reçoit en entrée les mesures Mb k,i , les données contenues dans les messages Mn k,i et les solutions PVTs k,i estimées par le module 50. Ce processus 52 fournit en sortie les biais c2 k,i établis. Ce processus 52 contient trois modules qui sont exécutés les uns après les autres : un module 56 de correction et de sélection des mesures, un module 58 d'estimation d'une position interne P1 k et un module 60 d’estimation des biais c2 k,i . Le processus 54 d’estimation de la solution PVT k a pour données d’entrée : les solutions PVTs k,i , les messages Mn k,i , les mesures brutes Mb k,i et les biais c2 k,i établis par le processus 52. Ce processus 54 se divise en deux modules exécutés l'un après l'autre : un module 62 de correction et de sélection des mesures Mb k,i qui délivre des mesures pré-corrigées Mc2 k,i et un module 64 d'estimation d'une solution externe PVT2 k . La solution externe est délivrée au calculateur 20. Ainsi, la solution PVT2 k et la solution PVTk sont les mêmes. Le processus 54 comporte aussi un soustracteur 66 qui retranche les biais c2 k,i aux mesures pré-corrigées Mc2 k,i délivrées par le module 62. Des exemples de réalisation des différents modules du logiciel 42 et le fonctionnement de l'unité 10 est maintenant décrit en référence au procédé de la . Pour chaque instant t k pour lequel une nouvelle solutions PVT k doit être délivrée au calculateur 20, l'unité 10 exécute les étapes suivantes. Lors d'une étape 100, l'étage 32 reçoit de nouveaux signaux GNSS émis par au moins trois satellites. Ces signaux GNSS sont qualifiés de "nouveau" dans le sens où ils sont reçus après les signaux GNSS utilisés pour estimer la précédente solution PVT k-1 . Typiquement, lors de l'étape 100, l'étage 32 reçoit les signaux GNSS émis par plus de quatre et habituellement plus de huit satellites. Lors d'une étape 102, l'étage 34 extrait, des nouveaux signaux GNSS reçus, les mesures brutes Mb k,i et les messages émis Mn k,i pour chacun des satellites. Ensuite, lors d'une étape 104, le module 50 estime les nouvelles solutions PVTs k,i à partir des mesures Mb k,i et des messages Mn k,i reçus. Lors d'une étape 110, le module 56 pré-corrige les mesures Mb k,i pour obtenir des mesures pré-corrigées Mc1 k,i . Ici, les mesures pré-corrigées Mc1 k,i sont obtenues à partir des mesures Mb k,i , des donnes contenues dans les messages Mn k,i , des solutions PVTs k,i estimée par le module 50 et de la précédente position interne P1 k-1 . Typiquement, les mesures pré-corrigées Mc1 k,i sont obtenues à l'aide de la relation suivante : où c1 k,i est la correction estimée lors de cette étape. La correction c1 k,i est établie en utilisant un modèle prédéterminé m1 des perturbations introduites par le fait que les signaux GNSS traversent la ionosphère et la troposphère. Ici, ce modèle m1 est paramétré par les solutions PVTs k,i et la position P1 k-1 . Par exemple, ici le modèle m1 est le modèle défini dans l'article suivant: Collins, J.P., Assessment and Development of a Tropospheric Delay Model for Aircraft Userw of the Global Positioning System , Technical Report No. 203, University of New Brunswick, 1999, pages 96-97 . Dans ce cas, la correction c1 k,i est calculée à l'aide de la relation (2) suivante : où : - l’angle θ est l’angle d’élévation, c'est-à-dire l'angle entre le plan tangent à la surface de la Terre au point défini par la position P1 k-1 et la droite reliant le satellite à ce point de la Terre, et - A(P1 k-1 ) est une fonction définie par la relation suivante : où : - h est la hauteur orthométrique de l'unité 10 obtenue à partir de la position P1 k-1 , et - pour chaque ligne de cette relation, le terme après le point virgule indique pour quelle plage de hauteurs orthométriques cette ligne doit être utilisée pour calculer A(P1 k-1 ), - m signifie l’unité de longueur mètre. La hauteur orthométrique h est estimée à partir la position P1 k-1 et d’un système géodésique. Lors d'une étape 112, le module 56 sélectionne les mesures pré-corrigées Mc1 k,i qui sont les seules à être utilisées dans la suite du processus 52. Le but de cette sélection est d’éliminer des mesures Mc1 k,i considérées comme étant trop erronées afin de ne pas les prendre en compte lors de l'estimation de la position P1 k . Cela permet d’améliorer la précision. Un exemple détaillé de mode de réalisation de cette étape 112 est décrit plus loin en référence à la . Les étapes suivantes du processus 52 utilisent uniquement les mesures sélectionnées lors de l'étape 112 pour calculer les biais c2 k,i . Les mesures non-sélectionnées ne sont pas utilisées pour calculer les biais c2 k,i de sorte que les biais c2 k,i pour lesquels l'indice i correspond à une mesure non-sélectionnée, ne sont pas mis à jour par le processus 52 pour l'instant t k . Autrement dit, pour les indices i qui correspondent à des mesures non-sélectionnées, le biais c2 k,i est pris égal à c2 k-1,i . De façon similaire, pour les indices i qui correspondent à des mesures non-sélectionnées, la variance du biais c2 k,i n'est pas mise à jour. Elle est donc égale à la variance du biais c2 k-1,i . Pour les étapes suivantes du processus 52, l'indice i correspond uniquement aux indices i des mesures sélectionnées. Lors d'une étape 114, le module 58 estime la position P1 k à partir des mesures pré-corrigées Mc1 k,i sélectionnées lors de l'étape 112, des solutions PVTs k,i estimées lors de l'étape 104 et de la position précédente P1 k-1 . Ici, le module 58 estime la position P1 k par minimisation de la fonction de coût f 1 suivante : où : - x k , y k , z k sont les coordonnées spatiales à déterminer de la position P1 k de l'unité 10, - x t,k est l'erreur d'horloge à déterminer de l'unité 10 exprimée sous la forme d'une distance ajoutée à la pseudo-distance calculée, - X est le vecteur de coordonnées x k , y k , z k et x t,k , - xs k,i , ys k,i et zs k,i sont les coordonnées spatiales du satellite i issues de la solution PVTs k,i estimée par le module 50, - Mc1 k,i est la i-ème mesure pré-corrigée lors de l'étape 110, - w i est un poids affecté à la i-ème mesure Mc1 k,i , et - Ipmin est l'ensemble qui contient les indices i des mesures sélectionnées lors de l'étape 112. La valeur de chaque poids w i est prédéterminée en fonction, par exemple, de l’élévation du satellite i, du rapport signal à bruit mesuré pour le satellite i, de la vitesse à laquelle varie la distance entre l'unité 10 et ce satellite. Dans une variante simplifiée, les poids w i sont constants et, par exemple, tous égaux à un. Les valeurs des coordonnées du vecteur X qui minimisent la fonction de coût f 1 sont estimées en utilisant un algorithme de Gauss-Newton par exemple, ou tout autre algorithme d’optimisation tels que les algorithmes de Levenberg-Marquardt, la descente de gradient ou autre qui sont généralement des algorithmes itératifs. Les valeurs initiales des coordonnées du vecteur X sont choisies, par exemple, égale à zéro. Ici, la recherche des valeurs des coordonnées du vecteur X qui minimisent la fonction f 1 est arrêtée : - après un nombre prédéterminé d'itérations, par exemple compris entre deux et vingt, ou - lorsque la mise à jour du vecteur X devient faible, c'est-à-dire que la norme de la différence entre le vecteur X obtenu à l'issue de l'itération courante et le vecteur X obtenu à l'issue de la précédente itération, est inférieure à un seuil prédéterminé. Ce seuil prédéterminé est, par exemple, compris entre 10 -4 ou 10 -10 . Lors d'une étape 116, le module 60 établit les valeurs des biais c2 k,i à partir des mesures pré-corrigées Mc1 k,i , des solutions PVTs k,i et de la position P1 k . Ici, les biais de mesure c2 k,i sont estimés en minimisant une fonction de coût, notée f 2, dépendant des biais de mesure c2 k-1,i établis pour l’instant précédent t k-1 . Par exemple, la fonction de coût f 2 est définie par la relation suivante : où : - Mc1 k est le vecteur formé par les mesures Mc1 k,i pour i appartenant à l'ensemble Ipmin, - c2 k-1 est le vecteur formé par les biais c2 k-1,i pour i appartenant à l'ensemble Ipmin. La fonction r(Mc1 k , P1 k ) est un vecteur de résidus de mesures entre Mc1 k et les mesures estimées en fonction de P1 k . La i-ème composante r i du vecteur de résidus est définie par la relation suivante : . Dans ce texte, le symbole désigne la norme au carrée du vecteur A issue du produit scalaire définie par la matrice définie positive R -1 . Cette norme est définie par la relation suivante : où, dans ce texte, le symbole " T " désigne l'opération transposée, et le symbole " -1 " désigne l'opération d'inversion d'une matrice. La matrice R c est la matrice de covariance des résidus du vecteur r(Mc1 k , P1 k ). Par exemple, ici, la matrice R c est diagonale et ses éléments sont tous égaux à 25 m². La matrice P c,k est la matrice de covariance des biais c2 k,i . Elle est mise à jour à pour chaque instant t k selon la relation suivante : Une fois les biais c2 k,i établis, l'exécution du processus 52 est terminée et l'exécution du processus 54 débute. Lors d'une étape 118, le module 62 pré-corrige les mesures Mb k,i pour obtenir des secondes mesures pré-corrigées Mc2 k,i . Ici, les mesures pré-corrigées Mc2 k,i sont obtenues à partir des mesures Mb k,i , des données contenues dans les messages Mn k,i , des solutions PVTs k,i estimées par le module 50 et de la précédente solution externe PVT2 k-1 . Ici, l'étape 118 est identique à l'étape 110 sauf que la position P2 k-1 contenue dans la solution PVT2 k-1 est utilisée à la place de la position P1 k-1 . Une fois les étapes 114 et 118 terminées, lors d'une étape 120, des mesures corrigées Mc3 k,i sont calculées à l’aide de la relation suivante : Ensuite lors d’une étape 122, le module 64 estime la solution PVT2 k à partir des solutions des satellites PVTs k,i , des mesures corrigées Mc3 k,i et la précédente solution PVT2 k-1 . Dans cet exemple de réalisation, la solution PVT2 k est estimée par minimisation de la fonction de coût f 3 suivante : où : - X est le variable permettant de minimiser la fonction f 3 , elle contient la solution PVT2 k,i et donc les éléments suivants : X T = (x, y, z, v x , v y , v z , x t , dx t ), avec : - x, y et z sont les coordonnées de la position P2 k de l’unité 10 dans le repère R T , - v x , v y et v z sont les composantes du vecteur vitesse V2 k de l’unité 10 dans le repère R T , - x t est le biais d’horloge de l’unité 10 par rapport au temps GPS, - dx t est la dérive d’horloge de l’unité 10 par rapport au temps GPS. - h(X) est la fonction de mesure dont la i-ème composante est définie par la relation suivante : - Mc3 k est le vecteur dont la i-ième composante est Mc3 k,i , - PVT2 k-1 est le vecteur PVT2 k estimé pour l’instant t k-1 , - R Mc3 est la matrice de covariance des mesures Mc3 k,i , elle est diagonale et chaque élément diagonal est, ici, fixé à 25m², - f est la fonction définie par la relation suivante : avec : - I n x p est la matrice identité de taille n par p, - 0 n x p est la matrice nulle de taille n par p, - dt est la période d’échantillonnage, ici elle est égale à t k -t k-1 , - P PVT2 est une estimation de la matrice de covariance de la solution PVT2 k , elle est obtenue à l’aide de la relation suivante : où : - P PVT2,k-1 est la matrice de covariance de la précédente solution PVT2 k-1 , et - Q est une matrice symétrique et positive dont les valeurs sont données par la relation suivante : Par exemple, la minimisation de la fonction de coût f 3 est effectuée avec la méthode de Gauss-Newton. Lors d’une étape 124, le module 64 met à jour l’estimation P PVT2 de la matrice de covariance de la solution PVT2 k à l’aide de la relation suivante : où : - H est la matrice jacobienne de la fonction h(X) définie ci-dessus lors de l’étape 122, et - P PVT2 est l’estimation de la matrice de covariance calculée lors de l’étape 122. Pour l’instant initial t 0 , la solution PVT2 -1 est définie comme étant le vecteur nul et la matrice de covariance P PVT2,-1 est la matrice diagonale définie par la relation suivante : Un exemple de mode de réalisation de l'étape 112 de sélection des mesures est maintenant décrit en détail en référence au procédé de la . Lors d’une opération 130, le module 56 calcule un score s1 0,i pour chaque satellite à l'aide de la relation suivante : où : - cn0 i (k-j) est la mesure du rapport signal à bruit pour le satellite i à l’instant t k-j , et - N est un nombre entier prédéterminé supérieur à un et, de préférence, supérieur à dix ou vingt. Typiquement, la mesure cn0 i (k-j) est réalisée lors de l'étape 102 par l'étage 34. Ensuite, lors d’une opération 132, le module 56 sélectionne les M meilleures mesures au sens du score s1 0,i . Par exemple, M est égal à huit. Enfin, un processus de détection et d’exclusion de fautes basé sur une approche par test de sous-ensembles est utilisé pour sélectionner le meilleur sous-ensemble de mesures parmi les M meilleures mesures. Pour limiter les calculs, le nombre maximal de fautes recherchées est limité à Mex=2. Plus précisément, la sélection du meilleur sous ensemble est effectuée comme suit. Lors d’une opération 134, le module 56 génère tous les sous-ensembles E p possibles de M -Mex mesures parmi les M meilleures mesures, où l'indice p est un identifiant du sous-ensemble. Ensuite, pour chacun des sous-ensembles E p générés, le module 56 exécute les opérations 136 à 140 suivantes. Lors de l’opération 136, le module 56 calcule de la position P1 k,p issue des mesures du sous-ensemble E p par minimisation de la fonction de coût f 1 suivante : Cette fonction f 1 est identique à la fonction f 1 utilisée lors de l’étape 114 sauf que l’ensemble Ip est l'ensemble qui contient les indices i des mesures contenues dans le sous-ensemble E p . La position qui minimise cette fonction de coût f 1 est trouvée en appliquant, par exemple, le même algorithme que celui mis en œuvre lors de l’étape 114. Lors d’une opération 138, le module 56 calcule les métriques suivantes pour chaque sous-ensemble E p : un résidu r p , une dilution géométrique GDOP p et une sensibilité maximale de l’erreur de position par rapport au résidu r p notée SlopeMax p . Le résidu r p est égal à la norme du vecteur r(Mc1 k,p , P1 k,p ), où : - Mc1 k,p sont les mesures contenues dans le sous-ensemble E p , et - r(…) est la même fonction que celle définie à l’étape 116. La dilution géométrique est plus connue sous le terme anglais de "Geometric dilution of precision". Ici elle est calculée en prenant en compte la position de chacun des satellites dont sont issues les mesures Mc1 k,i du sous-ensemble E p et de la position P1 k,p . Le paramètre SlopeMax p est calculé comme défini dans les pages 105 à 136 du document suivant : Salos Andres, C. D, Integrity monitoring applied to the reception of GNSS signals in urban environments , PhD Thesis, Institut National Polytechnique de Toulouse, 2012. Lors d’une opération 140, le module 56 calcule pour chaque ensemble E p un score s1 1,p à l'aide de la relation suivante : où α 1 , α 2 et α 3 sont des coefficients de pondération prédéterminés. Lors d’une opération 142, le module 56 sélectionne le sous-ensemble E p qui a le score s1 1,p le plus faible. Le sous-ensemble sélectionné est noté E pmin . Lors d’une opération 144, le module 56 valide le meilleur sous-ensemble E pmin , en comparant les métriques r pmin , GDOP pmin et SlopeMax pmin calculées pour le sous-ensemble E pmin , à des seuils prédéterminés, respectivement, Thr0, ThGDOP et ThSlopeMax. Si chacune des métriques r pmin , GDOP pmin et SlopeMax pmin est inférieure à son seuil respectif, alors le sous-ensemble E pmin est validé. Dans le cas contraire, la détection de faute est dite impossible, aucune mesure n'est validée. Lorsque aucune mesure n'est validée, les biais c2 k,i sont considérés comme inchangés et pris égaux aux biais c2 k-1,i . Seulement si le meilleur sous-ensemble E pmin a été validé, lors d’une opération 146, le module 56 calcule la cohérence des mesures restantes par rapport au sous-ensemble E pmin sélectionné. Les mesures restantes sont celles qui n'appartiennent pas au sous-ensemble E pmin . La cohérence est testée en comparant à un seuil prédéterminé Thr1 chaque résidu des mesures restantes. Le résidu de chaque mesure restante est calculé à l'aide de la relation suivante : où : - la fonction h() est la même fonction que celle mise en œuvre lors de l’étape 122, et - les indices i correspondent ici uniquement aux indices i des mesures restantes. Une mesure restante est cohérente avec les mesures du sous-ensemble E pmin si le résidu r i calculé pour cette mesure restante est inférieur à un seuil Thr1 prédéterminé. Ensuite, toutes les mesures restantes cohérentes avec les mesures du sous-ensemble E pmin , sont ajoutées dans le sous-ensemble E pmin et donc utilisées pour la suite du processus 52. A l’inverse, les mesures restantes incohérentes, ne sont pas utilisées pour la suite du processus 52. Chapitre II : Variantes : Variantes de l'unité de navigation : L'unité 10 de navigation a été décrite dans le cas particulier où elle est utilisée au sein d'un système de localisation comportant en plus l'unité 12 de mesure inertielle. Toutefois, l'unité 10 peut être utilisée dans tout système de localisation et, en particulier, des systèmes de localisation ne comportant pas d'unité de mesure inertielle ou comportant en plus ou à la place d'autres capteurs comme un odomètre. Les balises peuvent être des satellites d’une ou plusieurs constellations telles que GPS ("Global Positioning System"), Glonass, Galileo et Beidou. Ce qui a été décrit dans le cas particulier où les balises sont des satellites, peut être appliqué à tout système de localisation et, en particulier, à des systèmes de localisation où les balises ne sont pas des satellites. Par exemple, ce qui a été décrit ici peut aussi être mis en œuvre dans le cas où les balises sont des balises fixes sur terre qui émettent des messages. Par exemple, de telles balises fixes peuvent être installées dans un centre commercial pour localiser des piétons ou être installées à proximité d'un téléphérique pour localiser la position d'une télécabine. Ces balises fixes sont par exemple des balises ultra-large bande ou UWB ("Ultra Wide Band"). Le terme "balise ultra-large bande" ou "balise UWB" désigne ici un émetteur-récepteur qui utilise une large bande de fréquences pour envoyer et recevoir les messages. Une "large" bande de fréquences est une bande de fréquences dont la largeur est supérieure à 0,2f c , où f c est la fréquence médiane de cette bande de fréquences. Typiquement, une large bande de fréquences a une largeur supérieure à 250 MHz ou même supérieure à 400 MHz. A la place d'émettre des ondes électromagnétiques pour transmettre les messages, les balises peuvent émettre des ultrasons. Dans ce cas, la vitesse de propagation des messages dans l'espace est égale à la vitesse du son et non pas à la célérité de la lumière. Variantes du procédé de localisation : Pour mesurer une pseudo-distance, il n'est pas nécessaire que les messages émis par les balises comportent l'instant d'émission de ce message. Par exemple, ce n'est pas nécessaire si les horloges des balises et de l'unité de navigation sont synchronisées temporellement entre elles. Dans ce cas, par exemple, les balises émettent les messages à des instants te k connus de l'unité de navigation de sorte que ces instants te k n'ont pas besoin d’être contenus dans les messages émis. Selon une autre variante, les instants d'émission te k,i sont calculés par l'unité de navigation, par exemple, à partir d'un instant d'émission d'un message par l'unité de navigation et de l'instant de réception de la réponse à ce message émise par la balise. L'étape 110 de pré-correction de la mesure brute peut être réalisée différemment. En particulier, il existe d'autres modèles des perturbations introduites par l'ionosphère ou la troposphère utilisables à la place du modèle STANAG. En particulier, le modèle utilisé peut aussi être paramétré par des données contenues dans les messages Mn k,i . Dans le cas où les balises ne sont pas des satellites, un autre modèle des perturbations introduites par l'environnement extérieure dans lequel se déplace l'unité de navigation peut être utilisé. Dans un mode de réalisation simplifiée, l'étape 110 de pré-correction des mesures brutes est omise. En effet, dans certains contextes, comme le cas où les balises sont fixées à la surface de la terre, une telle pré-correction peut être inutile. L'étape 112 de sélection des mesures corrigées peut être réalisées différemment. Par exemple, l'une des méthodes décrites dans les articles suivants sous le nom anglais de "Fault Detection and Exclusion" peut être utilisée à la place de la méthode décrite dans le chapitre I : - H. Kuusniemi, G. Lachapelle : “GNSS Signal Reliability Testing in Urban and Indoor Environments”, Proceedings of the 2004 National Technical Meeting of The Institute of Navigation, San Diego, CA, pp. 210-224, 2004 - N. Zhu : “GNSS Propagation Channel Modeling in Constrained Environments: Contribution to the Improvements of the Geolocation Service Quality”, Ph.D thesis, University of Lille FR, 2018 Dans une autre variante de l'étape 112, au lieu d'utiliser la solution PVT2 k-1 pour effectuer la sélection, c'est la solution PVT1 k-1 qui est utilisée. D'autres méthodes de sélection du meilleur sous-ensemble de mesures parmi les M meilleures mesures au sens du score s1 0,i peuvent être utilisées. D’autres scores que les scores s1 0,i et s1 1,p peuvent aussi être utilisés pour sélectionner les meilleures mesures. Dans un mode de réalisation simplifié, l'étape 112 de sélection est omise. D’autres modes de réalisation de l’étape 114 d’estimation de la position P1 k,i sont possibles. Par exemple, en variante, c'est une solution interne PVT1 k,i qui est estimée et pas seulement la position interne P1 k,i . La vitesse interne V1 k,i peut être estimée à partir des mesures Doppler extraites par l’étage 34. Dans ce dernier cas, les différentes variantes décrites ci-dessous de l’étape 122 d’estimation de la solution PVT2 k,i sont applicables à l’étape 114. En variante, lors de l'étape 116 les données contenues dans les messages Mn k,i peuvent aussi être prises en compte pour établir les valeurs des biais c2 k,i . D'autres fonctions de coût sont utilisables à la place de la fonction de coût f 2 pour établir les biais c2 k,i . Par exemple, dans un mode de réalisation simplifié, les biais c2 k,i sont établis en minimisant la fonction de coût f 2.1 suivante : Les biais c2 k,i peuvent également être établis en utilisant un filtre de Kalman au lieu de minimiser une fonction de coût. Il existe encore d’autres méthodes pour estimer les biais c2 k,i . Par exemple, le processus 52 est remplacé par N processus 52 (p) qui s’exécutent en parallèle les uns des autres, où l’exposant (p) est un identifiant du processus. L’exposant (p) est un nombre entier compris entre 1 et N. Chaque processus 52 (p) fournit en sortie des biais respectifs notés c2 k,i (p) . Ici, les biais de mesures c2 k,i (p) sont établis par chaque processus 52 (p) en exécutant les étapes 110, 112, 114 et 116 précédemment décrites mais avec des paramètres différents de ceux utilisés par les autres processus 52 (p) . Dans cet exemple, la fonction f 2 est paramétrée par les matrices R c et P c,k-1 et chaque processus 52 (p) correspond à des valeurs spécifiques des coefficients de ces matrices R c et P c,k-1 . Par exemple, le processus 52 (1) utilise des matrices R c (1) et P c,k-1 (1) identiques aux matrices R c et P c,k-1 utilisées dans le cas du procédé de la . Le processus 52 (2) utilise des matrices, respectivement R c (2) et P c,k-1 (2) , dont les coefficients sont par exemple dix fois inférieures aux coefficients correspondants des matrices R c (1) et P c,k-1 (1) . A l’inverse, le processus 52 (3) utilise des matrices, respectivement R c (3) et P c,k-1 (3) , dont les coefficients sont par exemple dix fois supérieures aux coefficients correspondants des matrices R c (1) et P c,k-1 (1) . Ensuite, le meilleur jeu de biais c2 k,i (p) est sélectionné parmi les N jeux de biais c2 k,i (p) établis par les processus 52 (p) . Chaque biais c2 k,i est alors pris égal au biais c2 k,i (p) du meilleur jeu de biais c2 k,i (p) . Par exemple, le critère utilisé pour sélectionner le meilleur jeu de biais c2 k,i (p) est la norme de l’innovation notée nino (p) et chaque biais c2 k,i est égal au biais c2 k,i (p) du jeu de biais c2 k,i (p) qui possède la plus faible valeur nino (p) . Par exemple, la norme nino (p) du jeu de biais c2 k,i (p) est définie par la relation suivante : où : - Mc2 k est le vecteur dont les coefficients sont les mesures corrigées Mc2 k,i , - c2 k (p) est le vecteur dont les coefficient sont les biais c2 k,i (p), - PVT2 k-1 est le vecteur PVT2 k estimé pour l’instant t k-1 , - les fonctions h() et f() sont les mêmes que celles utilisées dans la fonction f 3 précédemment décrite. Après l’étape 116 d’établissement des biais c2 k,i et avant l’étape 120 de correction des mesures Mb k,i avec les biais c2 k,i , le procédé peut comporter d’autres étapes. Par exemple, le procédé comporte une étape de test sur la mise à jour des biais de mesure c2 k,i . Cette étape de test permet de détecter si la mise à jour des biais c2 k,i est conforme à un indicateur de confiance tel que l'écart type attendu. Par exemple, le test effectué est le test défini par la relation suivante pour chacun des indices i de l’ensemble I pmin : où : - r i (Mc1 k,i , P1 k ) est la même fonction que celle précédemment définie lors de l’étape 116, et - la notation (M) ii désigne le i-ème élément diagonal de la matrice M, - R c et P c sont les mêmes matrices que celles définies lors de l’étape 116, - α est un coefficient prédéterminé fixé par l'utilisateur. Si pour chacun des indices i de l’ensemble I pmin , le test ci-dessus est satisfait, alors la mise à jour est acceptée. Dans le cas contraire la mise à jour est supprimée est les valeurs des biais c2 k,i sont inchangées, c’est-à-dire qu’elles sont prises égales aux valeurs des biais c2 k-1,i . D’autres indicateurs de confiance peuvent être établis quant à la précision des biais c2 k,i établis. Par exemple, dans le cas où les biais c2 k,i sont établis par un filtre de Kalman, la covariance associée à chaque biais c2 k,i estimée est un tel indicateur de précision. Après l’étape 116 d’établissement des biais c2 k,i et avant l’étape 120 de correction des mesures Mb k,i avec les biais c2 k,i , le procédé peut aussi comporter une étape de filtrage des biais c2 k,i , par exemple à l’aide d’une moyenne glissante ou d’un filtre de Kalman, afin de les rendre leur estimation plus précise. C’est ensuite les biais filtrés qui sont utilisés pour obtenir les mesures corrigées Mc3 k,i . L'étape 118 peut être exécutée avant l'une des étapes 110, 112 et 114. En variante, lors de l'étape 118, les mesures pré-corrigées Mc2 k,i sont simplement prises égales aux mesures pré-corrigées Mc1 k,i . En variante l'étape 118 est omise. Entre les étapes 118 et 120, en variante, le module 62 sélectionne les mesures pré-corrigées Mc2 k,i qui sont ensuite les seuls à être utilisées dans la suite du processus 54. Par exemple, le module 62 réalise cette étape de sélection de la même manière que l’étape 112 ou ses variantes. De nombreuses variantes de l’étape 122 d’estimation de la solution PVT2 k,i sont possibles. Par exemple, la solution PVT2 k,i peut aussi être estimée en prenant en compte les données contenues dans les messages Mn k,i . Le calcul de la solution PVT2 k décrit utilise, entre autres données d’entrée, une matrice R Mc3 définissant la précision de chacune des mesures. Cette matrice R Mc3 peut ne peut pas être constante. Par exemple, elle est ajustée en fonction de la correction de biais appliquée à l’étape 118 ainsi que de la confiance accordée à cette correction, afin de mieux refléter l’incertitude de chaque mesure à l’issue de cette correction. De nombreux autres algorithmes peuvent être utilisés pour déterminer une position, une solution PT ou une solution PVT à partir des données d’entrée. Par exemple, des filtres de Kalman, des filtres particulaires ou des algorithmes d’intelligence artificielle peuvent aussi être utilisés. De plus, comme illustré dans le mode de réalisation du chapitre I, le calcul d’une position, d’une solution PT ou d’une solution PVT peut également se faire à partir d’un sous-ensemble des signaux d’entrée disponible. Autrement dit, certains signaux peuvent être délibérément exclus du calcul s’il est déterminé qu’ils présentent une mesure fortement erronée. Les algorithmes de sélection des signaux satellitaires sont nombreux et connus de l’homme de l’art. Dans un mode de réalisation simplifié, l'unité 10 délivre seulement sa position. Sa vitesse et son erreur d'horloge ne sont pas calculés ou ne sont pas délivrés au calculateur 20. Lorsque l'erreur d'horloge n'est pas calculée, le nombre minimum de satellites nécessaires pour déterminer la position de l'unité 10 est égal à trois. De plus, dans ce cas, les mesures brutes peuvent contenir seulement les pseudo-distances. Les autres mesures, comme l'effet Doppler, sont omises. Chapitre III : Avantages des modes de réalisation décrits : Puisque les biais sont établis en prenant en compte la position interne de l'unité de navigation, les biais sont plus précis que ceux calculés sans prendre en compte la position de l'unité de navigation. De plus, puisque les biais sont établis sans prendre en compte la position externe de l'unité de navigation, la stabilité et la robustesse du procédé de localisation est améliorée. En effet, une erreur lors de l'établissement des biais peut se propager sur l'estimation de la position externe. Par contre, l'erreur sur l'estimation de la position externe ne peut pas se propager sur l'établissement du biais. Ainsi, la stabilité et la robustesse du procédé de localisation sont accrues sans pour autant détériorer la précision. Le fait de sélectionner, à l'aide de la position externe à l'instant précédent t k-1 , les balises dont les messages sont ensuite traités pour établir les biais de mesure, permet d'améliorer la fiabilité et la précision de cette sélection car la position externe est plus précise que la position interne. Cela permet donc d'établir plus précisément les biais de mesure et donc d'obtenir une position externe plus précise à l'instant t k . De plus, cela ne compromet pas la robustesse du procédé de localisation car la position externe intervient seulement sur la sélection des messages traités. En particulier, la position externe n'intervient pas sur le contenus de ces messages, leurs instants de réception et lors des traitements exécutés pour établir les biais de mesure. Le fait de minimiser l'écart entre les biais de mesure pour les instants t k et t k-1 augmente la précision de la localisation. L'utilisation du modèle paramétré pour estimer l'erreur causée par la propagation des messages des balises à travers la ionosphère ou la troposphère permet d'augmenter la précision de la localisation. Procédé de localisation d'une unité de navigation dans un environnement comportant des balises qui émettent des messages, ce procédé comportant l'exécution des étapes suivantes pour chaque instant t k d'une suite temporelle d'instants {t 0 ; ...; t k-1 ; t k ; ...} : a) la réception (100) des messages émis par au moins trois balises, puis b) l'extraction (102) de mesures brutes de pseudo-distances qui séparent l'unité de navigation de chacune de ces au moins trois balises à partir des instants de réception de ces messages mesurés par l'unité de navigation et de la vitesse de propagation de ces messages dans l'espace, c) l'estimation (114) d'une position interne de l'unité de navigation à l'instant t k à partir des mesures brutes extraites pour cette instant t k et des positions des balises, d) pour plusieurs balises, l'établissement (116) d'un biais de mesure pour l'instant t k qui, lorsqu'il est retranché à la mesure brute extraite à partir de l'instant de réception du message émis par cette balise, permet de diminuer l'erreur entachant cette mesure brute, ce biais étant établi à partir de l'écart entre : - une pseudo-distance calculée qui sépare cette balise de l'unité de navigation, cette pseudo-distance calculée étant calculée à partir d'une position de cette balise et de la position interne de l'unité de navigation estimée pour l'instant t k , et - une pseudo-distance mesurée qui sépare cette balise de l'unité de navigation, cette pseudo-distance mesurée étant calculée à partir de la mesure brute de la pseudo-distance extraite pour cette balise pour l'instant t k et sans prendre en compte la position interne de l'unité de navigation, e) pour chaque balise pour laquelle un biais de mesure a été établi pour l'instant t k , le calcul (120) d'une mesure corrigée de la pseudo-distance à l'instant t k en ajoutant le biais établi pour cette balise à une mesure brute de la pseudo-distance pour cette balise extraite à l'instant t k , caractérisé en ce que le procédé comporte également l'étape suivante : f) l'estimation (122) d'une position externe de l'unité de navigation à l'instant t k à partir : - des positions des balises pour lesquelles une mesure corrigée de la pseudo-distance à l'instant t k a été calculée, - des mesures corrigées des pseudo-distances calculées pour l'instant t k , et - de la précédente position externe de l'unité de navigation estimée pour l'instant précédent t k-1 , et dans lequel les étapes b), c) et d) sont réalisées sans prendre en compte la position externe de l'unité de navigation. Procédé selon la revendication 1, dans lequel : - le procédé comporte, après l'étape b) et avant l'étape c), la sélection (112), en fonction de la position externe de l'unité de navigation à l'instant précédent t k-1 , d'un nombre limité de mesures brutes extraites parmi l'ensemble des mesures brutes extraites lors de l'étape b) pour l'instant t k , puis - lors des étapes c) et d) seuls les mesures brutes ainsi sélectionnées sont utilisées et les mesures brutes non-sélectionnées sont ignorées. Procédé selon l'une quelconque des revendications précédentes, dans lequel chaque biais de mesure pour l'instant t k est également établi (116) en minimisant l'écart entre ce biais de mesure pour l'instant t k et le biais de mesure établi pour l'instant précédent t k-1 . Procédé selon l'une quelconque des revendications précédentes, dans lequel chaque balise est un satellite en orbite autour de la terre. Procédé selon la revendication 4, dans lequel : - après l'étape b) et avant l'étape c), le procédé comporte une étape (110) de pré-correction des mesures brutes extraites, cette étape de pré-correction comportant pour chaque balise ayant émis un message à partir duquel est extrait une pseudo-distance : - l'établissement d'une erreur causée par la propagation du message de cette balise à travers la ionosphère ou la troposphère en utilisant pour cela un modèle prédéterminé paramétré par la position interne de l'unité de navigation estimée pour l'instant précédent t k-1 , puis - la soustraction de l'erreur établie à la mesure brute extraite pour obtenir une mesure brute pré-corrigée, puis - lors des étapes c) et d), ce sont les mesures brutes pré-corrigées qui sont utilisées. Support (40) d’enregistrement d’informations, lisible par un microprocesseur, comportant des instructions exécutables par ce microprocesseur, caractérisé en ce que ce support comporte des instructions non transitoires pour l’exécution d’un procédé de localisation conforme à l’une quelconque des revendications précédente, lorsque ces instructions sont exécutées par le microprocesseur. Unité (10) de navigation apte à se déplacer dans un environnement comportant des balises qui émettent des messages, cette unité de navigation comportant un microprocesseur (38) configuré pour exécuter les étapes suivantes pour chaque instant t k d'une suite temporelle d'instants {t 0 ; ...; t k-1 ; t k ; ...} : a) la réception des messages émis par au moins trois balises, puis b) l'extraction de mesures brutes de pseudo-distances qui séparent l'unité de navigation de chacune de ces au moins trois balises à partir des instants de réception de ces messages mesurés par l'unité de navigation et de la vitesse de propagation de ces messages dans l'espace, c) l'estimation d'une position interne de l'unité de navigation à l'instant t k à partir des mesures brutes extraites pour cette instant t k et des positions des balises, d) pour plusieurs balises, l'établissement d'un biais de mesure pour l'instant t k qui, lorsqu'il est retranché à la mesure brute extraite à partir de l'instant de réception du message émis par cette balise, permet de diminuer l'erreur entachant cette mesure brute, ce biais étant établi à partir de l'écart entre : - une pseudo-distance calculée qui sépare cette balise de l'unité de navigation, cette pseudo-distance calculée étant calculée à partir d'une position de cette balise et de la position interne de l'unité de navigation estimée pour l'instant t k , et - une pseudo-distance mesurée qui sépare cette balise de l'unité de navigation, cette pseudo-distance mesurée étant calculée à partir de la mesure brute de la pseudo-distance extraite pour cette balise pour l'instant t k et sans prendre en compte la position interne de l'unité de navigation, e) pour chaque balise pour laquelle un biais de mesure a été établi pour l'instant t k , le calcul d'une mesure corrigée de la pseudo-distance à l'instant t k en ajoutant le biais établi pour cette balise à une mesure brute de la pseudo-distance pour cette balise extraite à l'instant t k , caractérisé en ce que le microprocesseur (38) est également configuré pour exécuter l'étape suivante : f) l'estimation (122) d'une position externe de l'unité de navigation à l'instant t k à partir : - des positions des balises pour lesquelles une mesure corrigée de la pseudo-distance à l'instant t k a été calculée, - des mesures corrigées des pseudo-distances calculées pour l'instant t k , et - de la précédente position externe de l'unité de navigation estimée pour l'instant précédent t k-1 , et dans lequel les étapes b), c) et d) sont réalisées sans prendre en compte la position externe de l'unité de navigation..