FR 3148094 A1 20241025 FR 2304055 A 20230421 Méthode, mise en œuvre par ordinateur, d’apprentissage d’un modèle d’intelligence artificielle destiné à être utilisé pour déterminer une information de localisation d’un récepteur de radionavigation par satellite, la méthode comprenant les étapes de : Recevoir (402) un ensemble de données d’apprentissage comprenant plusieurs ensembles de mesures de radionavigation par satellite associés chacun à une information de positionnement de référence, Déterminer, pour chaque ensemble de mesures, un ensemble de métriques comprenant au moins un ensemble de résidus calculés pour plusieurs sous-ensembles de mesures excluant chacun au moins une mesure de l’ensemble, Pour chaque ensemble de mesures, Déterminer (404) un ensemble de coefficients de pondération de référence, Entrainer (403,405) le modèle d’intelligence artificielle à produire un ensemble de coefficients de pondération à partir des ensembles de métriques des données d’apprentissage, de manière à minimiser une distance entre lesdits coefficients de pondération et les coefficients de pondération de référence, Figure 4 Méthode de localisation GNSS assistée par un modèle d’intelligence artificielle L’invention concerne le domaine des systèmes de positionnement et navigation par satellite connus sous l’acronyme systèmes GNSS (Global Navigation Satellite System). L’invention concerne plus précisément une méthode de localisation d’un système GNSS, assistée par un modèle d’intelligence artificielle. La représente un schéma d’une architecture d’un récepteur GNSS selon l’art antérieur. Un tel récepteur REC a pour fonction d’acquérir et traiter des signaux de radionavigation par satellite ou signaux GNSS transmis par des satellites d’une ou plusieurs constellations parmi les constellations existantes GPS, Glonass, Galileo ou encore Beidou. Les signaux sont transmis dans une ou plusieurs bandes de fréquences, par exemple les bandes L1, L2 ou L5. Le récepteur REC comprend une antenne ANT, un étage d’entrée radiofréquence RF qui assure l’amplification, le filtrage, la transposition en fréquence et la numérisation des signaux analogiques reçus par l’antenne ANT. Il comprend également un étage BB de traitement en bande de base permettant l’acquisition des signaux, leur démodulation, l’extraction des données contenues dans les signaux pour générer des messages GNSS MG et des mesures brutes relatives à ces signaux MB. Le récepteur REC comprend enfin un processeur de navigation NAV configuré pour réaliser un calcul d’une solution de positionnement PVT qui comprend une information de positionnement du récepteur et/ou une information de vitesse de déplacement du récepteur et/ou une information de temps, c’est-à-dire l’horloge du récepteur par rapport à celle de la constellation des satellites. Plus précisément, l’étage de traitement en bande de base BB fournit au processeur de navigation NAV au moins les éléments suivants : Les mesures brutes MB contenant au moins l’une des mesures suivantes: les mesures de pseudo-distances, les mesures Doppler ou taux de pseudo-distance, les mesures de phase des ondes porteuses du signal GNSS et les instants de réception des signaux GNSS. Il peut également fournir des indicateurs de qualité de lien tels que le rapport signal-à-bruit (C/N0), Les messages GNSS MG contenant au moins : les éphémérides ou paramètres permettant d’estimer les positions et vitesses des satellites ; les paramètres de corrections d’erreurs d’horloge des satellites ; des éléments caractérisant l’état des satellites et l’état de leurs mesures, les paramètres d’un modèle d’estimation des délais ionosphériques. L’invention concerne la mise en œuvre du processeur de navigation NAV et plus précisément l’algorithme permettant d’estimer les informations de positionnement PVT. Le processeur de navigation NAV reçoit en entrée les mesures brutes MB et les messages GNSS MG et produit, typiquement, pour chacun des satellites pour lequel un signal est reçu des informations parmi : une mesure de pseudo-distance, une mesure de phase de porteuse, une mesure de décalage Doppler. Les mesures de pseudo-distance ou de phase peuvent être utilisées pour déterminer la position ainsi que le décalage d’horloge du récepteur, tandis que la mesure de décalage Doppler est utilisée pour la détermination de la vitesse et de la dérive d’horloge du récepteur. Par exemple, si l’on considère une mesure de pseudo-distance entre un récepteur de coordonnées et un satellite SVi de coordonnées , cette mesure respecte la relation suivante : 2.1 Où est la distance géométrique entre le récepteur et le satellite SVi, 2 . 2 est un décalage de distance résultant d’un décalage d’horloge entre le récepteur et la constellation (tous les satellites appartenant à une même constellation GNSS étant parfaitement synchronisés, le décalage est le même pour tous les satellites de la constellation) converti en mètres en le multipliant avec la vitesse de la lumière .c 2 . 3 est une erreur de mesure causée par de multiples phénomènes tels que : le bruit de mesure du récepteur, de l’ordre de quelques cm à quelques mètres, des retards de propagation dans l’ionosphère et la troposphère, de l’ordre de quelques cm à quelques mètres, des erreurs de position des satellites, des effets liés à la propagation du signal. Le fonctionnement du système GNSS est tel que les observations en provenance de différents satellites sont collectées simultanément au cours d’une période de temps appellee « époque » ( epoch = transmissions simultanées des signaux satellitaires) à l’issue de laquelle un ensemble de mesures est collectée. Le calcul de la solution de position ( ) dépend du mode de fonctionnement du récepteur. On distingue en effet un mode de poursuite et un mode dit “single-epoch”. Le mode de poursuite est un mode de fonctionnement du récepteur pour lequel une solution a déjà été calculée lors d’une époque précédente k- 1. Dans ce cas, la solution à l’instant k est obtenue par la mise à jour de la solution k -1 en utilisant les mesures collectées lors de l’époque k. Ce calcul est généralement effectué par un filtre de Kalman ou l’un de ses nombreux dérivés (EKF, UKF, etc.). Le mode dit « single-epoch » est un mode de fonctionnement du récepteur pour lequel aucune solution précédente n’est disponible. Ce mode, a priori moins performant que le mode de poursuite, trouve néanmoins un intérêt pour certaines applications telles que : l’initialisation de l’algorithme de poursuite, l’utilisation dans un contexte de fusion avec d’autres modalités de mesures telles que les centrales inertielles, où les solutions calculées sont filtrées à nouveau par un autre algorithme, dans certaines applications (ex : IoT pour Internet of Things) où le calcul est effectué à distance (dans un serveur distant) et pour lesquelles différentes contraintes (consommation, débit) ne permettent qu’un transfert ponctuel des mesures depuis les nœuds/tags capteurs à localiser. L’invention concerne plus particulièrement, mais pas uniquement, le mode dit « single-epoch » en ce qu’il vise le calcul d’une solution de positionnement pour une époque donnée à l’instant k sans bénéfice d’une solution antérieure. En mode dit single epoch , le calcul d’une solution de position ( ) comprenant quatre inconnues peut donc être déterminé à partir de l’ensemble des mesures collectées à l’époque k en utilisant l’équation 2.1. Puisque le nombre de mesures N (typiquement, de 4 à 50) est généralement supérieur au nombre d’inconnues qui est au moins égal à 4, la solution doit être donc calculée au sens d’un critère de minimisation des erreurs affectant les mesures – encore appelés résidus - tel que les moindres carrés. Plus précisément, le nombre d’inconnus est égal à 4 lorsque tous les signaux reçus sont émis par des satellites de la même constellation (par exemple la constellation GPS). Si les signaux sont émis par des satellites d’au moins deux constellations (par exemple GPS et GALILEO) alors le nombre d’inconnues est égal à 5. De façon générale, le nombre d’inconnues est égal à 3+m, où m est le nombre de constellations de satellites envisagé. La relation 2.4 donne l’équation des moindres carrés permettant de déterminer la solution de position X en recherchant la solution qui minimise une somme pondérée des résidus élevés au carré. 2. 4 Les résidus correspondent à la difference entre la prediction de la mesure selon le modèle 2.1 pour la solution X d’une part, et la mesure réalisée d’autre part, où x,y,z et δ sont les solutions de l’équation 2.4. 2 . 5 Et définit le coefficient de pondération associé à la mesure # i , c’est-à-dire « l’importance » que celle-ci aura dans le calcul de la solution. Le terme définit le coefficient de pondération ou poids associé à chaque mesure . Plus celui-ci est élevé, plus la mesure associée aura une contribution élevée sur la solution calculée. Un poids nul signifie qu’une mesure n’a aucune influence sur la solution calculée et est donc équivalent à l’absence de cette mesure. Dans le cas où le poids est non nul, son importance est évaluée relativement aux poids des autres mesures, sa valeur seule ne permet donc pas de définir son importance. Le même type d’équation que l’équation 2.5 peut être établie en utilisant les mesures de phase plutôt que celles de pseudo-distance, les deux grandeurs étant identiques dans leur modèle, à un coefficient près. Le même type d’équation peut aussi être établi pour calculer une solution de vitesse ( ), incluant les composantes de vitesse ainsi que la dérive d’horloge du récepteur, à partir des mesures de doppler. Dans la suite du document, l’invention est décrite pour la recherche d’une solution de position à partir de mesures de pseudo-distances mais elle peut ainsi être étendue directement à la recherche d’une solution de position à partir de mesures de phase ou à la recherche d’une solution de vitesse à partir de mesures de Doppler en appliquant le même algorithme. Il existe plusieurs stratégies de pondérations pour la résolution des moindres carrés donnée par l’équation 2.4. Une solution possible consiste à calculer le poids en fonction de l’erreur probable associée à chaque mesure en utilisant le modèle de l’équation 2.6, c’est-à-dire l’inverse de la variance de l’erreur probable sur la mesure d’indice i. 2.6 Cette stratégie de pondération trouve en effet des fondements théoriques (estimation Bayésienne), qui permettent de la justifier sous certaines hypothèses. Dans ce cas, l’erreur probable est l’écart type associé à la distribution de l’erreur de la mesure , c’est-à-dire la distribution du terme aléatoire de l’équation 2.1. Outre sa justification théorique, cette approche permet de relier la valeur exprimée en mètres, à des phénomènes physiques mesurables. Ainsi, la distribution des erreurs de mesure peut être établie empiriquement par des campagnes de test, ou à partir d’autres éléments connus (précision de positionnement des satellites, propagation dans l’ionosphère, bruit du récepteur etc.). Typiquement, les valeurs sont comprises entre 1m et 20m, mais peuvent varier en fonction de l’intensité du signal reçu ou de l’élévation du satellite par exemple. Les techniques reliant différentes métriques annexes (rapport signal à bruit C/N0, élévation du satellite etc.) à la précision supposée de la mesure sont appelées pondération stochastique. Par exemple, la variance de l’erreur peut être obtenue au moyen du modèle empirique donné dans le document [1] et repris à l’équation 2.7. 2.7 Dans laquelle et représentent respectivement l’angle d’élévation ainsi que l’intensité du signal reçu du satellite SVi, et et sont deux termes empiriques. D’autres modèles empiriques peuvent être utilisés en remplacement de celui donné ci-dessus. Si ces formules empiriques sont largement utilisées, elles n’en demeurent pas moins difficiles à régler, et une modélisation plus fine de la distribution de chacune des mesures reçues permettrait une amélioration des performances. Cependant, cette modélisation plus fine peut s’avérer extrêmement complexe en pratique, tant le nombre de phénomènes en jeu sont nombreux (état des couches atmosphériques, type de récepteur, précision de positionnement du satellite, présence d’obstacles dans l’environnement du récepteur etc.) et difficiles à prédire. Par ailleurs, afin de garantir l’optimalité, et donc une bonne précision, des solutions calculées, il est important que l’ensemble des mesures prises en compte pour établir cette solution soient conformes au modèle utilisé a priori, c’est-à-dire dans une plage d’erreur attendue. Or, plusieurs phénomènes peuvent significativement dégrader la qualité d’une mesure tels que : un dysfonctionnement du satellite, une attaque – par exemple de type brouillage ou leurrage- à l’encontre des signaux GNSS, la présence d’obstacles (par exemple des bâtiments en environnement urbain) bloquant le trajet direct de l’onde, conduisant le récepteur à ne mesurer qu’un trajet réfléchi qui aura parcouru une distance potentiellement très supérieure à celle attendue. Ces mesures biaisées doivent donc être détectées comme telles, et exclues du calcul de la solution ou bien se voir assigner un poids nul, ce qui est équivalent. On parle alors de détection de fautes. Il existe donc un besoin pour une méthode de localisation qui permette de mieux prendre en compte les incertitudes liées aux modèles théoriques de calculs des coefficients de pondération lors de la résolution de la solution par les moindres carrés. Il existe aussi un besoin pour une méthode de localisation qui permette d’éliminer de façon fiable les mesures biaisées du calcul de la solution de positionnement. Les schémas de surveillance d'intégrité existants, tels que la surveillance d'intégrité autonome du récepteur connue sous l’acronyme RAIM pour « Receiver Autonomous Integrity Monitoring », déterminent s'il y a une défaillance dans une mesure de satellite en examinant la cohérence d'un ensemble de mesures redondantes. Une façon de faire est d'utiliser la méthode de séparation de solution. La méthode de séparation de solution pour les récepteurs RAIM est basée sur le calcul de la différence entre une solution de navigation complète qui est rendue en utilisant tous les N satellites visibles et un ensemble de sous-solutions de navigation qui sont rendues en utilisant N-i satellites visibles. En calculant l'ensemble de sous-solutions de navigation, l’algorithme RAIM suppose qu'une seule défaillance de satellite se produit à la fois. Cependant, avec l'introduction de nouvelles constellations (par exemple Galileo, BeiDou) et l'utilisation continue des constellations existantes (par exemple les systèmes de positionnement mondiaux (GPS et GLONASS), il est plus probable qu'il puisse y avoir plusieurs défaillances simultanées de satellite à un moment donné. De plus, les défaillances de constellations entières devront également être prises en compte par les futurs schémas de surveillance d'intégrité. En réponse à la possibilité de plusieurs défaillances se produisant en même temps, la surveillance d'intégrité autonome avancée du récepteur ou ARAIM « Advanced Receiver Autonomous Integrity Monitoring » a été développée. La technologie ARAIM est basée sur la méthode de séparation de solution, mais elle a été modifiée pour prendre en compte l’apparition de plusieurs défaillances simultanées et les défaillances de constellation. Pour chaque défaillance à surveiller, une sous-solution de navigation qui n'inclut pas les mesures associées à la défaillance est créée. Par exemple, si deux défaillances simultanées simples doivent être surveillées, alors un ensemble de sous-solutions basées sur l'élimination de toutes les combinaisons possibles de deux satellites doit être créé. Cependant, l'augmentation du nombre de satellites visibles et des probabilités plus élevées de défaillances simultanées (comme prévu avec les nouvelles constellations) peut considérablement augmenter le nombre de sous-solutions à créer. Cela a un impact important sur les exigences de calcul de l'algorithme, ce qui se traduit par des puces plus coûteuses pour l’implémentation des récepteurs. De nombreuses techniques de détection et d’exclusion de fautes ( FDE, Fault detection and exclusion en anglais) ont été développées, celles-ci peuvent être classées en deux catégories principales, selon que la faute est analysée dans l’espace des mesures ( Residual-Based RAIM, notées ici RB) comme décrit dans [3] ou dans l’espace des solutions ( Solution- Separation RAIM, notées ici SS ) comme décrit dans [2]. Dans le cas des méthodes RB , une première solution est calculée en utilisant toutes les mesures, puis un test statistique est effectué sur l’ensemble des résidus afin de déterminer si une faute est présente parmi les mesures. En effet, si toutes les mesures respectent le modèle attendu (par exemple modèle Gaussien, d’un écart type donné) alors la somme des résidus au carré est donnée par la relation suivante 3.1 3.1 La somme des résidus suit une loi du chi-2 et cette hypothèse peut être testée en comparant la valeur obtenue à un seuil T. 3.2 Dans le cas des méthodes SS, une première solution est calculée en utilisant toutes les mesures, ainsi que N solutions basées chacune sur N -1 mesures en excluant le satellite SV i . Ensuite, chacune des solutions est comparée (en norme) par rapport à la solution , 3.3 Si aucune d’elles ne s’écarte plus d’un certain seuil T prédéterminé, 3.4 alors cela signifie qu’aucune mesure ne comporte de faute. Dans les deux cas, si le test est non concluant, cela signifie qu’il existe au moins une faute parmi les mesures. Afin d’exclure la faute, on peut construire N sous ensembles de N -1 mesures et recommencer l’opération. Si l’un des tests est satisfait, alors la totalité des mesures associées au sous-ensemble sont considérées comme valides et la solution calculée peut être retenue. Si en revanche aucun des tests n’est satisfait, l’opération est répétée en construisant des sous-ensembles de N -2 mesures qui seront ensuite testés individuellement, et ainsi de suite. Il est à noter qu’à chaque itération correspondant à un nombre h de fautes testées, le nombre de sous-ensembles à construire pour chaque hypothèse de nombre de fautes h, est donné par la relation : 3.5 Ce nombre peut rapidement devenir considérable lorsque le nombre de fautes testées devient important. Par exemple, un récepteur standard multi-constellations peut typiquement obtenir 40 mesures, et, en environnement urbain, le nombre de fautes peut fréquemment se situer entre 3 et 10, conduisant à des nombres de sous-ensembles à tester compris entre 9880 et 847660528 ce qui peut rapidement devenir prohibitif en cout de calcul. Il est à noter que toutes ces approches ont en commun de reposer sur une décision (via un test statistique), celle d’inclure ou non certaines mesures, mais ne modifient pas les poids associés à chacune d’elles, sinon par un forçage à 0, ce qui est équivalent à une exclusion. D’autres approches connues sont basées sur une modification des pondérations telles que la méthode proposée dans la référence [5]. Cette approche se distingue par l’assignation de coefficients de pondération non nuls à des mesures pourtant détectées comme étant affectées d’une faute, ici principalement causée par du multi-trajet. L’idée est de préserver un équilibre entre d’une part, la réduction du poids des mesures fautives dans le calcul de la solution et, d’autre part, la dégradation de la dilution géométrique de la précision causée par l’exclusion des mesures fautives. La dilution géométrique de la précision concerne un problème d’imprécision des mesures due à des facteurs géométriques. Autrement dit, les positions relatives du récepteur et des satellites agissent de façon significative sur la précision de la position. Par exemple, pour une précision de mesure identique, des satellites « plutôt alignés » dans le ciel conduiront à une moins bonne précision de positionnement que des satellites répartis de manière homogène dans le ciel. Ce phénomène, appelé dilution géométrique de la précision ou simplement dilution de la précision peut être quantifié précisément à partir des positions même approximatives (>1km) du récepteur et des satellites. Il se représente souvent sous la forme de termes multiplicatifs dans les trois directions de l’espace, qui permettent de traduire l’erreur de mesure en erreur de position dans ces mêmes directions. Par exemple, une dilution horizontale (HDOP) de la précision de 2, signifie qu’une erreur de mesure de 1m (supposée identique pour l’ensemble des satellites) se traduit par une erreur de positionnement horizontale de 2m. La méthode propose dans [5] consiste à : Détecter toutes les mesures fautives, Réduire itérativement le poids des mesures fautives jusqu’à ce que un nombre maximum d’itération est atteinte, ou la dilution géométrique de la précision, calculée en prenant en compte la nouvelle pondération, dépasse un certain seuil, Les mesures non fautives utilisent les formules classiques de l’état de l’art telles que celle de l’équation 2.7 par exemple. Une autre approche proposée dans [8] s’intéresse au réglage des formules de pondération des mesures dans le cas de multi-trajets avec une notion d’apprentissage (sous la forme d’une régression linéaire des valeurs), cependant elle ne considère pas comme entrées pour le calcul des pondérations des informations conjointes (c.-a-d. reposant sur des calculs combinant plusieurs signaux, comme les résidus par exemple). La référence [9] décrit l’utilisation d’un réseau de neurones convolutionnel (CNN) pour classifier les signaux GNSS selon trois catégories : vue directe ( line-of- sight , LOS), non vue directe ( non-line of sight , NLOS) et multi-trajets ( multipath , MP). Cependant cet algorithme utilise comme entrées directement les sorties des corrélateurs du ranging processor , qui sont des données internes et généralement non accessibles à l’extérieur du récepteur, et n’utilise pas d’entrées conjointes, c’est-à-dire mettant en jeu plusieurs signaux issus de différents satellites simultanément. De plus, les sorties associées à chaque signal/satellite, sont de nature discrète. En synthèse, il existe donc un besoin pour améliorer les méthodes de l’art antérieur en particulier dans la prise en compte plus précise des erreurs et incertitudes d’origines diverses qui impactent les mesures de radionavigation par satellite. L’invention vise à proposer une nouvelle méthode de localisation qui a pour objectif d’améliorer la précision de positionnement ou de vitesse grâce à une meilleure pondération des mesures. L’invention exploite un modèle d’intelligence artificielle qui est entrainé pour apprendre à déterminer, à partir de mesures conjointes réalisées sur l’ensemble des signaux reçus au sein d’une époque, un ensemble de coefficients de pondération permettant un calcul plus précis d’une information de localisation via une résolution au sens des moindres carrés. L’invention a pour objet une méthode, mise en œuvre par ordinateur, d’apprentissage d’un modèle d’intelligence artificielle destiné à être utilisé pour déterminer une information de localisation d’un récepteur de radionavigation par satellite, la méthode comprenant les étapes de : Recevoir un ensemble de données d’apprentissage comprenant plusieurs ensembles de mesures de radionavigation par satellite associés chacun à une information de positionnement de référence, Déterminer, pour chaque ensemble de mesures, un ensemble de métriques comprenant au moins un ensemble de résidus calculés pour plusieurs sous-ensembles de mesures excluant chacun au moins une mesure de l’ensemble, Pour chaque ensemble de mesures, Déterminer un ensemble de coefficients de pondération de référence, Entrainer le modèle d’intelligence artificielle à produire un ensemble de coefficients de pondération à partir des ensembles de métriques des données d’apprentissage, de manière à minimiser une distance entre lesdits coefficients de pondération et les coefficients de pondération de référence, les coefficients de pondération étant destinés à pondérer un ensemble de résidus, égaux à une différence entre une mesure et une information de positionnement prédite, lors d’un calcul de prédiction d’une information de positionnement par minimisation de la somme des résidus élevés au carré et pondérés par les coefficients de pondération. Selon un aspect particulier de l’invention, chaque coefficient de pondération de référence est associé à un satellite et est fonction d’un résidu calculé comme la différence entre une pseudo-distance associée à une mesure de radionavigation émise par ce satellite et une pseudo-distance de référence calculée à partir de l’information de positionnement de référence. Selon un aspect particulier de l’invention, l’information de positionnement de référence est fournie par un moyen de positionnement pris parmi : un système inertiel ou un système GNSS à haute précision ou une combinaison des deux systèmes. Selon un aspect particulier de l’invention, le modèle d’intelligence artificielle est un réseau de neurones artificiels, par exemple un réseau de neurones récurrents. Selon un aspect particulier de l’invention, ensemble de mesures des données d’apprentissage comprend en outre un ensemble d’indicateurs de qualité du signal reçu pour chaque mesure. L’invention a aussi pour objet une méthode de localisation d’un récepteur de radionavigation par satellite comprenant les étapes de : Recevoir un ensemble de N mesures de radionavigation par satellite, N étant un entier au moins égal à 4, Déterminer, à partir des mesures reçues, un ensemble de métriques comprenant au moins un ensemble de résidus calculés pour plusieurs sous-ensembles de mesures comprenant chacun au plus N-1 mesures, Exécuter une phase d’inférence du modèle d’intelligence artificielle entrainé au moyen de la méthode d’apprentissage selon l’invention, à partir dudit ensemble de métriques pour déterminer un ensemble de coefficients de pondérations, Déterminer une information de localisation à partir desdites mesures de radionavigation par satellite et des coefficients de pondération en recherchant la valeur de l’information de localisation qui minimise la somme des résidus élevés au carré, associés aux N mesures, pondérés par les coefficients de pondération. Selon un aspect particulier de l’invention, les mesures de radionavigation par satellite sont des mesures de pseudo-distance et l’information de localisation est une position d’un récepteur. Selon un aspect particulier de l’invention, un résidu est défini par la différence entre une mesure de pseudo-distance et une prédiction de cette mesure. Selon un aspect particulier de l’invention, les mesures de radionavigation par satellite sont des mesures de phases ou de Doppler et l’information de localisation est une position ou une vitesse d’un récepteur. Selon un aspect particulier de l’invention, l’ensemble de métriques comprend en outre un ensemble d’indicateurs de qualité du signal reçu pour chaque mesure. L’invention a aussi pour objet un récepteur de signaux de radionavigation par satellite comprenant un étage de réception de signaux GNSS, un étage de traitement en bande de base et un processeur de navigation configuré pour exécuter les étapes de la méthode de localisation selon l’invention. L’invention a aussi pour objet un programme d'ordinateur comportant des instructions pour l'exécution de la méthode selon l’invention, lorsque le programme est exécuté par un processeur. L’invention a aussi pour objet un support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution de la méthode selon l’invention, lorsque le programme est exécuté par un processeur. D’autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés suivants. représente un schéma d’un récepteur de radionavigation par satellite selon l’art antérieur, représente un schéma d’un récepteur de radionavigation par satellite selon l’invention, représente un organigramme détaillant les étapes principales d’une méthode de détermination d’une information de localisation d’un récepteur GNSS selon un mode de réalisation de l’invention, représente un organigramme détaillant les étapes d’une méthode d’apprentissage d’un modèle d’intelligence artificielle destiné à être utilisé pour déterminer une information de localisation, selon un mode de réalisation de l’invention, représente un organigramme détaillant les étapes d’une méthode de prédiction d’une information de localisation à partir du modèle d’intelligence artificielle entrainé par la méthode de la , selon un mode de réalisation de l’invention, représente une matrice de résidus destinée à alimenter un moteur d’intelligence artificielle selon un premier exemple de réalisation, représente une matrice de résidus destinée à alimenter un moteur d’intelligence artificielle selon un second exemple de réalisation. La description de l’invention est, par la suite, faite pour le cas d’une détermination de position à partir de mesures de pseudo-distances, mais l’invention s’applique de façon similaire pour déterminer une information de vitesse ou de temps à partir de mesures de pseudo-distances, de phases ou de Doppler. La schématise l’architecture fonctionnelle d’un récepteur de signaux de radio-navigation par satellite selon un mode de réalisation de l’invention. Le récepteur comporte un premier étage de réception REC qui comporte une antenne, un étage radiofréquence RF et un étage de traitement en bande de base BB similaires à ceux déjà décrits à la . L’invention peut être implémentée dans un processeur numérique de navigation NAV, comme représenté à la . Alternativement, les mesures fournies par l’étage de réception REC peuvent être transmises à un serveur distant dans lequel est implémentée l’invention. L’étage de réception REC est apte à recevoir différents signaux de radio-navigation émis par différents satellites émettant chacun sur une ou plusieurs fréquences et appartenant à une même constellation ou à des constellations différentes. Les signaux sont émis durant une époque. L’étage de réception REC fournit en sortie des mesures de pseudo-distances (ou de phase ou de Doppler) ainsi que différents indicateurs de qualité de chaque signal reçu, par exemple un indicateur de rapport signal à bruit C/N0 ou un indicateur d’instant d’accrochage ou « lock time » en anglais. Un premier module EXT d’extraction de métriques reçoit les mesures et les indicateurs et fournit en sortie un ensemble de métriques caractérisant les signaux. Cet ensemble de métriques comprend au moins un ensemble de métriques conjointes, c’est-à-dire des métriques calculées à partir de plusieurs signaux. Les métriques conjointes sont par exemple un ensemble de résidus organisés dans une matrice MR comme représentée à la , les résidus étant calculés à partir de chaque sous-ensemble de N-1 satellites parmi N satellites en visibilité du récepteur et pour lesquels des signaux ont été reçus. Optionnellement, l’ensemble de métriques comprend également des métriques par lien qui correspondent ou sont calculées à partir des indicateurs de lien fournis par l’étage de réception REC. Un deuxième module IA d’intelligence artificielle est configuré pour fournir en sortie un ensemble de coefficients de pondération pour chacune des mesures. Par exemple, le module IA est implémenté au moyen d’un réseau de neurones artificiel entrainé au moyen de techniques d’apprentissage machine. Un troisième module RMC de résolution des moindres carrés détermine une solution de positionnement à partir des mesures de pseudo-distances fournies par l’étage de réception REC et des coefficients de pondération fournis par le module IA en appliquant l’équation 2.4. La schématise, sur un organigramme, les étapes principales de l’invention. A l’étape 301, le modèle d’intelligence artificielle utilisé est initialisé, en particulier ses paramètres sont initialisés en vue d’apprendre une fonction f qui permet de calculer les coefficients de pondération associés à des mesures . A l’étape 302, le modèle d’intelligence artificielle est entrainé afin d’ajuster ses paramètres pour apprendre à réaliser la fonction f. L’entrainement est réalisé sur un nombre K de jeux de mesures par exemple associés à k positions connues. L’entrainement est réalisé de manière que les coefficients de pondération obtenues soient le plus proche possible de coefficients de référence . Le nombre K de jeux de mesures utilisé est tel que K.N soit supérieur au nombre de paramètres P ( K.N > P ) et est par exemple compris entre 10000 et 10000000. N est le nombre de mesures réalisées à partir de signaux émis par des satellites différents. Avantageusement, les coefficients de pondération de référence sont calculés en utilisant l’équation 2.6 à partir d’une solution de position de référence connue qui est par exemple fournie par un récepteur à haute précision qui peut être muni d’une centrale inertielle ou tout autre moyen indépendant des signaux GNSS. Les coefficients de pondération de référence correspondent à des résidus définis par l’équation 4.1 4.1 Les poids de référence sont calculés en utilisant un modèle empirique, par exemple celui donné à l’équation 2.7. Les paramètres sont ajustés par exemple de manière à minimiser une fonction de cout C qui est fonction d’une distance entre les coefficients de pondération obtenus et les coefficients de pondération de référence. La fonction de cout C est par exemple donnée par la relation 4.2 : 4.2 A l’étape 303, le modèle d’intelligence artificielle entrainé à l’étape 302 est ensuite utilisé pour prédire, à partir de nouvelles mesures, un ensemble de coefficients de pondération qui sont utilisés pour calculer une solution de positionnement à partir de l’équation des moindres carrés 2.4. On décrit à présent plus en détail chaque étape de l’invention. Le modèle d’intelligence artificielle utilisé peut être un réseau de neurones artificiel, par exemple un réseau de neurones récurrents ou tout autre modèle qui peut être entrainé pour apprendre à exécuter une fonction particulière. Selon un mode de réalisation particulier de l’invention, le modèle d’intelligence artificielle est un réseau de neurones récurrents ou LSTM-NN pour « Long Short Term Memory Neural Network ». Le réseau comporte, par exemple deux couches cachées de respectivement 893 et 517 neurones. La couche d’entrée comporte autant de neurones que de métriques fournies par le module d’extraction de métriques EXT et la couche de sortie comporte un nombre N de neurones égal au nombre de coefficients de pondération à générer qui est lui-même égal au nombre maximum de mesures associées à des satellites différents qu’un récepteur peut recevoir. Ce paramètre peut être défini a priori. Les fonctions d’activation des couches du réseau de neurones sont par exemple, une fonction tangente hyperbolique tanh ou une fonction d’activation ReLU. Les paramètres du réseau de neurones, en particulier les coefficients synaptiques sont initialisés à une valeur prédéterminée, par exemple 0 ou 1 ou une valeur aléatoire entre 0 et 1. Le nombre de couches (1 à 100), le nombre de neurones (5 à 1000 par couche) de même que le nombre de paramètres P = 1000 à 100 000 000 peut varier. La « fonction d’activation » de chaque couche peut être choisie parmi des exemples tels que ReLU, tanh, ou sigmoid par exemple. Le nombre de couches, de même que le nombre de neurones et le choix de la fonction d’activation peuvent eux même faire partie des paramètres pouvant être optimisés lors de l’étape d’apprentissage 302. D’autres types de réseaux peuvent également être utilisés tels que des réseaux de type totalement connectés « fully connected » ou convolutionnels « CNN » (convolutionnal-NN) ou encore des réseaux GRU (« Gate Recurrent Unit ») par exemple. Les données d’entrée du réseau de neurones sont normalisées entre 0 et 1 en divisant chacune des valeurs par une valeur prédéterminée (exemple 10, 100 ou 1000) ou par la valeur maximale des valeurs calculées sur toutes les données d’entrée. La représente, sur un organigramme, les étapes de mise en œuvre d’une méthode d’apprentissage d’un modèle d’intelligence artificielle selon un mode de réalisation de l’invention. La méthode reçoit en entrée les paramètres déterminés lors de la phase d’initialisation représentée par l’étape 401 sur la . L’étape 402 consiste à générer une base de données d’apprentissage pour réaliser l’apprentissage du modèle. Les données d’apprentissage sont constituées de plusieurs ensembles de test, chaque ensemble étant constitué de plusieurs mesures de radionavigation par satellite et d’une solution de positionnement de référence associée . La base de données peut être construite en collectant des mesures réelles au moyen d’un récepteur et en mesurant la position de référence à partir d’un système très performant ( ), par exemple un système inertiel ou un système GNSS à haute précision ou une combinaison des deux systèmes. Dans ce cas, le décalage d’horloge du récepteur (non fourni par le système de référence, car celui-ci est propre au récepteur utilisé pour la collecte) peut être calculé en utilisant la relation 5.1 (on omet l’indice k pour une meilleure lisibilité). 5.1 Selon une variante de réalisation, les données d’apprentissage comprennent en plus des mesures de pseudo-distances, d’autres métriques associées à chacun des signaux reçus tels que, par exemple, le rapport signal à bruit C/N0 ou la durée de verrouillage du signal. Selon une autre variante de réalisation, la base de données d’apprentissage est générée par simulation. Les données d’apprentissage sont organisées sous forme de métriques identiques à celle produites par le module d’extraction de métriques, c’est-à-dire une matrice de résidus complétée éventuellement par des métriques spécifiques à chaque lien. A l’étape 403, pour chaque ensemble de mesures de la base de données d’apprentissage, une phase d’inférence du réseau de neurones est exécutée pour calculer les coefficients de pondération à partir des paramètres du réseau. A l’étape 404, un ensemble de coefficients de pondération de référence sont déterminés. Avantageusement, ils sont calculés à partir de la position de référence et des mesures de références en utilisant, par exemple l’équation 2.6 ou une équation similaire. Par exemple, les coefficients de pondération de référence peuvent être calculés au moyen de l'une des formules suivantes: , A l’étape 405, le réseau de neurones est entrainé de sorte à ajuster ses paramètres en fonction : des coefficients de pondération calculés : , des coefficients de pondération de référence : . Plus précisément, l’ajustement des paramètres est par exemple réalisé au moyen de : une phase de propagation des données d’apprentissage dans le réseau de neurones depuis les entrées du réseau vers ses sorties pour calculer des coefficients de pondération , un calcul d’erreur entre la sortie calculée et le résultat à obtenir via une fonction de cout, une phase de rétropropagation de l’erreur qui vise à ajuster les paramètres ou coefficients synaptiques pour chaque connexion entre deux neurones à l’aide d’un algorithme de rétropropagation qui vise à minimiser la fonction de cout. L’algorithme de rétropropagation est par exemple basé sur une descente de gradient ou tout autre algorithme équivalent. La fonction de cout est une fonction de distance entre les coefficients de pondération calculés et les coefficients de pondération de référence de sorte à faire converger les coefficients de pondération calculés vers les coefficients de référence. Par exemple, la fonction de cout suivante peut être utilisée en remplacement de celle donnée à l’équation 4.2 : 5.2 Toute autre fonction de cout traduisant un écart entre les coefficients calculés par l’algorithme et les coefficients de référence peut être utilisée. Une fois le modèle entrainé, il peut être utilisé pour prédire une information de positionnement. La schématise les étapes de mise en œuvre d’une méthode de localisation selon un mode de réalisation de l’invention. A l’étape 501, on reçoit un nouvel ensemble de mesures acquises par un récepteur de radio-navigation, par exemple un ensemble de mesures de pseudo-distances correspondant à une même époque GNSS. Le nombre N de mesures réalisées peut varier typiquement de l’ordre de 4 à 50 mesures. Le récepteur peut recevoir des signaux satellites de différentes constellations (GPS, Galiléo Glonass, Beidou, Qzs etc.). Le récepteur peut recevoir plusieurs signaux à des fréquences différentes (bandes L1, L2, L5, E1, E2 etc.) pour un même satellite. Ainsi, un seul satellite peut conduire à la réception de plusieurs signaux par exemple, si toutefois le récepteur dispose de cette capacité (récepteur multi-bandes). Le récepteur peut produire, pour chaque signal reçu, plusieurs types de mesures parmi : Des mesures de pseudo-distance, de phase de la porteuse, des mesures Doppler. Des mesures de rapport signal à bruit C/N0, des mesures de temps d’accrochage du récepteur. Les mesures de porteuses ou de Doppler peuvent être utilisées à la place des mesures de pseudo-distances. Dans ce cas la matrice de résidus est calculée sur les mesures de phase ou sur les mesures de Doppler. Les mesures de porteuses ou de Doppler peuvent être utilisées en plus des mesures de pseudo-distances. Dans ce cas, plusieurs matrices de résidus peuvent être calculées et présentées en entrée de l’algorithme. Dans une variante de réalisation, une sélection est opérée au sein des mesures. Certaines mesures peuvent être exclues d’emblée des traitements car supposées non-fiables. Un test préliminaire de validité peut s’appuyer, par exemple, sur une élévation minimale du satellite i : , avec choisi entre 1 et 15 deg par exemple (i.e. la mesure est conservée si l’élévation du satellite est supérieure à la valeur minimale ) ; ou sur des critères de niveau de signal : . A l’étape 502, on extrait des mesures un ensemble de métriques. Dans un mode de réalisation, l’étape 502 consiste à construire une matrice de résidus M telle qu’illustrée à la . Une telle matrice est de dimensions N par N, où N est le nombre de mesures associées à des satellites différents. Chaque ligne de la matrice correspond à un sous ensemble de mesures dans lequel on a exclu un satellite. On a donc N sous ensembles de N-1 mesures. Pour chaque sous ensemble on calcule N-1 résidus en utilisant les équations 2.4 et 2.5 et en utilisant un modèle empirique pour la détermination des coefficients de pondération. Sur chaque ligne de la matrice M on a donc N-1 valeurs de résidus avec j l’indice du satellite et Xi correspond à l’ensemble de mesures d’indice i. Les résidus sont calculés en se limitant aux N-1 mesures de l’ensemble Xi. Le terme diagonal de la matrice peut être remplacé par une valeur arbitraire prise dans la plage [1 1000] par exemple ou peut correspondre au calcul d’un résidu pour la mesure exclue. Afin que la matrice conserve une taille constante de dimensions NxN (ce qui facilite l’implémentation) si le nombre de mesures reçues est inférieur à N , alors tous les coefficients manquants sont affectés d’une valeur prédéfinie . D’autres indicateurs peuvent être ajoutés à la matrice M tels que les indicateurs de qualité de chaque lien ( , LT i ), ou l’élévation du satellite . On peut également ajouter des indicateurs statistiques dérivant des indicateurs de qualité du lien. Par exemple, on peut adjoindre les statistiques empiriques telles que la valeur moyenne et l’écart type calculées sur une fenêtre glissante de M échantillons précédent l’époque courante, ainsi que le nombre d’échantillons utilisés pour calculer ces statistiques empiriques ( WS i ) : , Dans le cas ou d’autres informations relatives à chaque lien sont utilisées en plus des résidus, celles-ci peuvent être avantageusement ajoutées à la matrice de résidus comme illustré sur la formant une matrice M de caractéristiques composée d’une matrice de résidus MR ainsi qu’une matrice de métrique par lien ML. Dans tous les cas, les métriques déterminées à l’étape 502 doivent être de même nature que celles utilisées pour entrainer le réseau de neurones lors de la phase d’apprentissage. L’utilisation de métriques conjointes telles que des résidus calculés pour différents sous-ensembles de mesures permet de fournir au modèle d’intelligence artificielle des informations plus complètes que les seules métriques de qualité de lien individuelles pour chaque mesure. A partir de l’ensemble de ces métriques, le modèle d’intelligence artificielle peut apprendre à identifier quelle mesure est susceptible d’être entachée d’une erreur et de l’affecter alors d’un coefficient de pondération de valeur plus faible que les autres mesures. Si le nombre de métriques obtenues à l’étape 502, qui dépend du nombre de mesures réalisées à l’étape 501, est inférieur au nombre de métriques utilisé lors de la phase d’apprentissage (correspondant au nombre d’entrées du réseau de neurones), alors les entrées supplémentaires du réseau de neurones sont complétées par des valeurs nulles ou choisies arbitrairement. A l’étape 503, le moteur d’intelligence artificielle entrainé lors de la phase d’apprentissage est exécuté pour déterminer un ensemble de coefficients de pondération à partir des métriques calculées à l’étape 502 qui lui sont fournies en entrée. Les coefficients de pondération obtenus permettent de pondérer avec des valeurs plus élevées les mesures identifiées comme fiables et avec des valeurs plus faibles les mesures identifiées comme erronées. La caractérisation des mesures comme fiable ou erronée est possible car le réseau de neurones a été entrainé pour générer les meilleurs coefficients possibles dans un nombre important de situations différentes couvertes par les données d’apprentissage. En particulier les données d’apprentissage couvrent, avantageusement, de nombreux cas d’application dont des scénarii pour lesquels certaines mesures sont entachées d’erreurs. Dans une variante de réalisation, le moteur d’intelligence artificielle peut être entrainé pour produire une sortie qui n’est pas directement un coefficient de pondération tel que défini dans l’équation 2.4, mais relié à celui-ci via une fonction g comme par exemple : 5.4 Dans ce cas, lors de la phase d’entrainement, les coefficients de référence sont transformés de la même façon en utilisant la même fonction g. A l’étape 504, une solution de positionnement est calculée à partir des coefficients de pondération déterminés à l’étape 503 par le réseau de neurones. La solution est calculée à partir de l’équation 2.4 appliquée aux mesures et aux coefficients calculés à l’étape 503. Différents algorithmes peuvent être utilisés pour résoudre l’équation 2.4, tels que les moindres carrés récursifs, l’algorithme de Levenberg-Marquardt ou la descente de gradient. Lorsque les mesures sont issues de satellites appartenant à plusieurs constellations différentes, il faut estimer un décalage d’horloge par constellation. La définition de la solution devient donc et le calcul du résidu de l’équation 2.5 devient l’équation 5.6: 5.6 L’équation 5.6 utilise le décalage approprié , où C désigne la constellation de la mesure considérée. Le calcul de la solution peut s’appuyer sur une solution précédente (époque k -1) qui peut être utilisée pour initialiser la résolution de l’équation 2.4 et ainsi accélérer le traitement. Le calcul de la solution peut également être effectué par des algorithmes de filtrage (type filtre de Kalman), qui effectuent une mise à jour de la solution obtenue à l’époque k -1, (i.e ) pour calculer la solution à l’époque courante (i.e ). Dans le cas où l’étape 503 ne produit pas des coefficients de pondération en sortie tels que définis à l’équation 2.4 (cf variante de l’étape 503 et équation 5.4), ceux-ci seront tout d’abord transformés par une fonction g( ) avant d’être utilisés. La méthode de localisation décrite à la permet de fournir une information de positionnement comprenant au moins l’une des indications parmi une position, une vitesse ou une référence temporelle d’horloge. La méthode de localisation peut être exécutée dans un processeur de navigation NAV d’un récepteur GNSS du type décrit à la . Alternativement, la méthode de localisation peut être exécutée dans un serveur distant à partir de mesures réalisées par un récepteur GNSS puis transmises vers ce serveur. De façon générale, les méthodes d’apprentissage et de localisation selon l’invention peuvent être mises en œuvre à l’aide de composants matériels et/ou logiciels. Les éléments logiciels peuvent être disponibles en tant que produit programme d’ordinateur sur un support lisible par ordinateur, support qui peut être électronique, magnétique, optique ou électromagnétique. Les éléments matériels peuvent être disponibles tous ou en partie, notamment en tant que circuits intégrés dédiés (ASIC) et/ou circuits intégrés configurables (FPGA) et/ou en tant que circuits neuronaux ou en tant que processeur de signal numérique DSP et/ou en tant que processeur graphique GPU, et/ou en tant que microcontrôleur et/ou en tant que processeur général par exemple. Méthode, mise en œuvre par ordinateur, d’apprentissage d’un modèle d’intelligence artificielle destiné à être utilisé pour déterminer une information de localisation d’un récepteur de radionavigation par satellite, la méthode comprenant les étapes de : Recevoir (402) un ensemble de données d’apprentissage comprenant plusieurs ensembles de mesures de radionavigation par satellite associés chacun à une information de positionnement de référence, Déterminer, pour chaque ensemble de mesures, un ensemble de métriques comprenant au moins un ensemble de résidus calculés pour plusieurs sous-ensembles de mesures excluant chacun au moins une mesure de l’ensemble, Pour chaque ensemble de mesures, Déterminer (404) un ensemble de coefficients de pondération de référence, Entrainer (403,405) le modèle d’intelligence artificielle à produire un ensemble de coefficients de pondération à partir des ensembles de métriques des données d’apprentissage, de manière à minimiser une distance entre lesdits coefficients de pondération et les coefficients de pondération de référence, les coefficients de pondération étant destinés à pondérer un ensemble de résidus, égaux à une différence entre une mesure et une information de positionnement prédite, lors d’un calcul de prédiction d’une information de positionnement par minimisation de la somme des résidus élevés au carré et pondérés par les coefficients de pondération. Méthode d’apprentissage d’un modèle d’intelligence artificielle selon la revendication 1 dans laquelle chaque coefficient de pondération de référence est associé à un satellite et est fonction d’un résidu calculé comme la différence entre une pseudo-distance associée à une mesure de radionavigation émise par ce satellite et une pseudo-distance de référence calculée à partir de l’information de positionnement de référence. Méthode d’apprentissage d’un modèle d’intelligence artificielle selon l’une quelconque des revendications précédentes dans laquelle l’information de positionnement de référence est fournie par un moyen de positionnement pris parmi : un système inertiel ou un système GNSS à haute précision ou une combinaison des deux systèmes. Méthode d’apprentissage d’un modèle d’intelligence artificielle selon l’une quelconque des revendications précédentes dans laquelle le modèle d’intelligence artificielle est un réseau de neurones artificiels, par exemple un réseau de neurones récurrents. Méthode d’apprentissage d’un modèle d’intelligence artificielle selon l’une quelconque des revendications précédentes dans laquelle chaque ensemble de mesures des données d’apprentissage comprend en outre un ensemble d’indicateurs de qualité du signal reçu pour chaque mesure. Méthode de localisation d’un récepteur de radionavigation par satellite comprenant les étapes de : Recevoir (501) un ensemble de N mesures de radionavigation par satellite, N étant un entier au moins égal à 4, Déterminer (502), à partir des mesures reçues, un ensemble de métriques comprenant au moins un ensemble de résidus calculés pour plusieurs sous-ensembles de mesures comprenant chacun au plus N-1 mesures, Exécuter (503) une phase d’inférence du modèle d’intelligence artificielle entrainé au moyen de la méthode d’apprentissage selon l’une quelconque des revendications précédentes, à partir dudit ensemble de métriques pour déterminer un ensemble de coefficients de pondérations, Déterminer (504) une information de localisation à partir desdites mesures de radionavigation par satellite et des coefficients de pondération en recherchant la valeur de l’information de localisation qui minimise la somme des résidus élevés au carré, associés aux N mesures, pondérés par les coefficients de pondération. Méthode selon la revendication 6 dans laquelle les mesures de radionavigation par satellite sont des mesures de pseudo-distance et l’information de localisation est une position d’un récepteur. Méthode selon la revendication 7 dans laquelle un résidu est défini par la différence entre une mesure de pseudo-distance et une prédiction de cette mesure. Méthode selon la revendication 6 dans laquelle les mesures de radionavigation par satellite sont des mesures de phases ou de Doppler et l’information de localisation est une position ou une vitesse d’un récepteur. Méthode selon l’une quelconque des revendications 6 à 9 dans laquelle l’ensemble de métriques comprend en outre un ensemble d’indicateurs de qualité du signal reçu pour chaque mesure. Récepteur de signaux de radionavigation par satellite comprenant un étage de réception de signaux GNSS (RF), un étage de traitement en bande de base (BB) et un processeur de navigation (NAV) configuré pour exécuter les étapes de la méthode de localisation selon l’une quelconque des revendications 6 à 10. Programme d'ordinateur comportant des instructions pour l'exécution de la méthode selon l’une quelconque des revendications 1 à 10, lorsque le programme est exécuté par un processeur. Support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution de la méthode selon l’une quelconque des revendications 1 à 10, lorsque le programme est exécuté par un processeur.