FR 3061560 A1 20180706 FR 1601897 A 20161229 1 Procédé d'ajustement optimal de bornes d'erreurs GIVE ou de calcul optimal des variances de résidus de points IGP d'une grille ionosphérique de correction d'un système SBAS et système SBAS de mise en oeuvre dudit procédé La présente invention concerne un premier procédé d'ajustement optimal de bornes d'erreurs des estimations des retards ionosphériques, appelés GIVE (en anglais « Grid lonospheric Vertical Error »), participant à assurer l'intégrité des services de navigation aéronautique critiques (en anglais « Safety of Life), fournie par un système d'augmentation spatial (en anglais « Satellite Based Augmentation System »). Ces corrections ainsi que leurs bornes d'erreurs associées sont élaborées pour des points, appelés IGP, d'une grille ionosphérique de correction du système SBAS. L'invention concerne également un deuxième procédé de calcul optimal des écarts-type appelés o-mops, permettant de modéliser l'erreur d'interpolation, réalisée entre les IGP et le point de percement du signal reçu par l'utilisateur, et l'erreur de la fonction de rabattement du délai vertical (en anglais « vertical delay ») vers le délai oblique (en anglais « slant delay »), lesdits écarts-type o-MoPs intervenant dans le calcul des GIVE finaux, et étant associées aux corrections d'erreur ionosphérique élaborées pour des points IGP d'une grille ionosphérique de correction d'un système SBAS pour assurer l'intégrité d'un service GNSS. L'invention concerne également un système SBAS d'augmentation à base de satellite(s) configuré pour la mise en oeuvre du premier procédé 25 et/ou du deuxième procédé de l'invention. Un système SBAS d'augmentation à base de satellite(s), par exemple le système EGNOS (en anglais European Geostationary Navigation Overlay Service), est un système destiné à augmenter les performances d'un système global de navigation par satellites GNSS, par exemple le système 30 GPS (en anglais « Global Positioning System »), en fournissant une précision améliorée et en garantissant un niveau d'intégrité rendant l'utilisation du système GPS compatible des opérations aériennes. De manière connue, les corrections et les bornes d'erreurs associées à cette utilisation, sont élaborées par des centres de calcul, faisant partie du 35 système SBAS et utilisant des données récupérées à partir de récepteurs 3061560 2 GNSS ou de stations spécifiques déployés de manière optimisée dans la zone de service à augmenter, ces stations étant appelées des stations RIMS (en anglais « Ranging and Integrity Monitoring Station ») dans le cas d'un système de type EGNOS (en anglais « European Geostationary Navigation 5 Overlay Services »). Le système SBAS élabore des corrections permettant de compenser les principales sources d'erreur du positionnement GPS à savoir : .- la correction de l'orbite et de l'horloge des satellites augmentés ; .- et également la correction de l'erreur ionosphérique liée au retard du 10 signal GPS causé par la traversée dudit signal dans l'ionosphère. L'ensemble des corrections et de leurs bornes d'erreur associées sont envoyées aux utilisateurs, plus précisément les récepteurs GNSS embarqués des terminaux utilisateurs, par l'intermédiaire de satellites géostationnaires émettant des signaux standardisés qui pourront ensuite, 15 être décodés et utilisés par les utilisateurs équipés pour l'augmentation SBAS. Généralement, un utilisateur équipé pour l'augmentation SBAS dispose d'un récepteur GPS spécifique, capable de traquer et décoder les messages des satellites géostationnaires SBAS, conformément aux normes définis par le standard applicable RTCA-D0229, toutes versions comprises 20 jusqu'à la version E actuelle, et le standard ICAO SARPs, toutes versions comprises jusqu'à la 6ième version actuelle. Concernant les erreurs d'horloge et d'orbite, les éphémérides directement reçues au travers des signaux GPS sont généralement de très bonne qualité et ne nécessitent pas forcément de correction. Cependant, c'est surtout l'intégrité garantie sur ces données qui ne suffit pas à l'aviation civile. Le système SBAS doit donc observer, corriger et surtout garantir avec le niveau d'intégrité requis les informations d'orbite et d'horloge qui seront envoyées aux utilisateurs finaux. Le système SBAS a donc en particulier la capacité de détecter les erreurs et de les réparer, ou si cela ne s'avère pas possible d'avertir l'utilisateur dans un temps suffisamment court de l'existence de ces erreurs non réparées, pour garantir sa sécurité. Concernant le retard ionosphérique, le problème est tout autre. En effet, non corrigées, ces erreurs peuvent peser plusieurs dizaines de mètres sur l'erreur de positionnement GPS. Bien que les signaux envoyés par les 35 satellites GPS contiennent un modèle de correction de ces retards en 3061560 3 fonction de l'heure et de la position des points de percement IPP (en anglais lonospheric Pierce Point) des signaux reçus par l'utilisateur, ce modèle s'avère peu fiable, en particulier lorsque l'activité électromagnétique de l'ionosphère est forte et instable, notamment dans les zones équatoriales ou 5 polaires. Ce modèle portant, le nom de son inventeur Klobuchar, ne peut donc servir pour garantir l'intégrité que pour des opérations aériennes ne nécessitant que peu de précision,». Pour les approches précises lesdites procédures étant dénommées classiquement en anglais « Non Precision Approach / En route en vue de l'atterrissage d'un avion, ce modèle est donc 10 insuffisant compte-tenu du niveau de précision requis. Le système SBAS doit donc élaborer une estimation de ces retards ou délais au-dessus de la zone de couverture en définissant un modèle plus précis et actuel du contenu électronique de l'ionosphère en terme de retard de traversée de l'ionosphère des signaux GNSS. Ce modèle, plus précis et 15 représentatif en temps réel des retard ionosphériques, est réalisé au travers d'un maillage ou réseau discret virtuel de points appelés IGP (en anglais « lonospheric Grid Point »), situés à 350 km d'altitude, l'ionosphère étant assimilée à une couche mince concentrant le principal de sa densité électronique à cette altitude. 20 Ces délais de traversée calculés sur les IGPs sont établis pour des signaux verticaux, et sont interpolés de façon normalisée par le terminal utilisateur en fonction de la localisation de ses IPP par rapport aux différents points IGP du maillage. Une fonction de rabattement, standardisée également, est ensuite utilisée en fonction de l'élévation par rapport à la 25 sphère terrestre des différentes lignes de vues des satellites visibles depuis l'antenne GNSS du terminal utilisateur pour assurer l'allongement du délai en fonction de l'obliquité des signaux traqués. Les délais et leur borne d'erreur sur chacun des points IGP du maillage observables par le système SBAS sont calculés à partir des 30 mesures bi-fréquence collectées par les stations de télémétrie et de contrôle d'intégrité réparties sur l'ensemble de la zone de couverture du système SBAS, lesdites stations étant dénommées RIMS dans le cas du système EGNOS et WRS (en anglais « Wide-area Reference Station ») dans le cas du système WAAS (en anglais « Wide area Augmentation System ») par exemple. 3061560 4 Un algorithme dédié réalise l'estimation des délais et fournit une première borne d'erreur simplifiée appelée « Basic GIVE » (en anglais « Basic Grid lonospheric Vertical Error »). Un certain nombre de termes sont ensuite calculés et ajoutés afin de conférer à la borne d'erreur finale assurant 5 l'intégrité, dénommé le GIVE (en anglais « Grid lonospheric Vertical Error »), la confiance requise pour être conforme à la sécurité relative aux opérations aériennes précises. D'après le standard du système GPS, le GIVE est déterminé suivant l'expression : 10 GIVE = 3,29 x 0G1VE dans lequel o-ân,E désigne une variance corrigée de la variance fournie par le filtre estimateur du GIVD, désignée par o-, DASICGIVE Ici et suivant un exemple d'algorithme dédié, la variance corrigée o-lvE est calculée suivant l'expression : 15 ,.2 ___ ,2 ,2 + ,r2 er.2 + ''USERsmooTHiNG + ,T2 uGIVE = uBASICGIvE + ',MOPS , ',GEOMETRYscALE , -DRIFT dans laquelle ,r2 '--BASICGIVE est la variance fournie par le filtre estimateur du GIVD, oHOpS est la variance des résidus obtenus par la différence entre 20 l'UISD estimé à partir des IGP avec la méthode d'interpolation du système MOPS, et le STEC (« Slant Total Electron Content ») obtenu directement du point de mesure de la station considérée, ,2 uGEOMETRYSCALE est une variance représentative de la qualité de l'observation du ciel autour de l'IGP à estimer, 25 ,.2 ''USERSMOOTHING est une variance de protection des utilisateurs contre un biais généré par le filtrage de leurs mesures en présence de gradient sur le délai ionosphérique, ,.2 ''DRIFT est une variance qui permet d'assurer le maintien de la protection offerte par le GIVE jusqu'au prochain rafraîchissement de celui- 30 Ci. Une fois le GIVE final élaboré pour chaque IGP, il est vérifié de nouveau que lorsqu'on applique les corrections ionosphériques aux différents IPP mesurés par les stations RIMS que l'on compare par rapport à la valeur directement mesurée à leur niveau, les différences entre la valeur de délai 3061560 5 ionosphérique mesuré à l'IPP et la valeur reconstruite par interpolation sont correctement couvertes par leur UIVE (User lonospheric Vertical Error) associé, résultant de l'interpolation (selon la méthode d'interpolation imposée par le standard) des GIVE des IGP situés autour des IPP considérés. C'est 5 sur ce point précis que se concentre le premier procédé de l'invention. Avant de décrire le procédé actuel de vérification d'intégrité des informations de correction associées aux IGP, des définitions et des illustrations préliminaires d'un certain nombre de grandeurs utilisés par ledit procédé actuel sont décrites ci-dessous. 10 Suivant la Figure 1 un point de percement IPP, désigné par la référence numérique 2, correspond à l'intersection d'une ligne de vue 4 LoS (en anglais « Line of Sight ») entre un satellite GNSS 8, ici un satellite GPS, et un récepteur utilisateur 10 avec l'ionosphère 12, assimilée à une couche 14 infiniment mince par la norme RTCA-D0229, toutes versions comprises 15 jusqu'à la version E actuelle, et la norme ICAO SARPs, et située à 350 km d'altitude. Le récepteur utilisateur 108, ici embarqué à bord d'un aéronef, et le satellite 8 GPS font partie d'un système 14 GNSS. De la même façon qu'un récepteur utilisateur possède à un instant donné un IPP par ligne de vue du ou de chacun des satellites GNSS en 20 visibilité, suivant la Figure 2 chaque station d'observation 22 RIMS d'un système SBAS 24 possède un nombre donné d'IPP mesurés sur une fenêtre temporelle donnée. Ici, sur la Figure 2, un seul point de percement d'observation IPP 26 à un instant donné est représenté comme l'intersection d'une ligne de vue 28 entre un satellite GPS 30 et un terminal récepteur 32 25 d'une station d'observation 34, faisant partie de l'ensemble des stations d'observations 22 du système SBAS 24, avec l'ionosphère 12, assimilée à la couche 14 infiniment mince. Les mesures fournies par les stations d'observation RIMS 22 étant en bi-fréquence aux fréquences L1 (C/A) et L2 (P(Y)), contrairement aux 30 terminaux utilisateurs qui ne disposent actuellement que de la fréquence L1 (C/A) pour les usages civils, le délai ionosphérique en chacun des IPP d'observation peut être calculé à partir de la différence faites sur les « pseudo-distances », désignées respectivement par C1 et P2, entre L1 et L2, après retranchement par estimation des biais inter-fréquences générés 35 par le terminal récepteur GPS 32 de la RIMS 34 de la Figure 2, et par le 3061560 6 satellite 30 également. Au niveau du satellite 30, le biais correspond à un biais de transmission entre les signaux de fréquences différentes. Au niveau récepteur, il s'agit d'un biais provoqué par la chaîne de traitement du terminal récepteur entre chaque fréquence. 5 La mesure du délai ionosphérique par cette méthode repose sur le fait que la vitesse de propagation d'un signal électromagnétique tel que le signal GPS dans l'ionosphère dépend de sa fréquence, ici notée f. Plus précisément, le délai ionosphérique, ici noté d, est au premier ordre proportionnel au contenu électronique du chemin parcouru STEC (en anglais 10 « Slant Total Electronic Content ») et inversement proportionnel au carré de la fréquence, ce qui peut être formulé par l'expression : d STEC P2 - Cl = (y - 1) x STEC IFB"tellite IFBstation (équation #1) 25 avec: f2 Pour information, la fréquence L2 n'est que partiellement utilisée par le système SBAS qui ne fait que traquer le signal en lui-même, car le contenu du signal porté par cette fréquence L2 est à usage militaire exclusif 15 et donc inaccessible pour des utilisations civiles du fait de puissantes méthodes de cryptage qui y sont introduites. Cependant, pour la mesure du retard ionosphérique, cette utilisation est amplement suffisante. Les avions civils ne sont pas quant à eux équipés de récepteurs capables de traquer L2 ne serait-ce que partiellement du fait de l'absence de 20 standard encadrant son utilisation, et surtout du fait que L2 n'est pas une fréquence civile à la base. Ainsi la différence des pseudo-distances est régie par la relation: fi, 2 ^ - Y- , le rapport des fréquences L1 et L2 permettant de ramener fL2 le délai ionosphérique à la fréquence de référence L1; ^ STEC, le délai ionosphérique de la ligne de vue réelle, donc fonction 30 de l'élévation (en anglais Slant TEC, avec TEC "Total Electron Content") ; wBsatellite du satellite et du récepteur. Ce calcul est répété pour toutes les lignes de vues de chaque RIMS déployée dans toute la zone de service du système SBAS. Ainsi on obtient 35 un nombre important d'IPP de mesures situés aux alentours et dans des et IFBstation sont respectivement les biais inter-fréquence 3061560 7 voisinages de tous les points de maillage IGP observables dans la zone de service sur lesquels le STEC est mesuré. Puis, chaque délai ionosphérique de la ligne de vue réelle STEC est « redressé » à la verticale en devenant un VTEC (V pour vertical) suivant 5 une relation du type : STEC = Fpp. VTEC (équation #2) dans laquelle Fpp désigne une fonction de rabattement du délai ionosphérique en fonction de l'élévation de la ligne de vue. On pourra noter ici que ce facteur vaut 1 quand l'élévation de la ligne 10 de vue est à 90°, donc verticale. A cette étape du calcul, le choix de la loi de rabattement est libre. On n'utilise donc pas nécessairement la loi de rabattement définie dans le standard pour un terminal utilisateur dans le cadre du service de navigation. On peut définir la loi que l'on juge la plus efficace et la plus adaptée pour l'obtention du meilleur GIVD possible. 15 Ensuite, un algorithme est utilisé pour calculer un délai ionosphérique estimé sur chaque IGP, appelé alors délai vertical du réseau maillé ionosphérique GIVD (en anglais « Grid lonosphere Vertical Delay »). Suivant la Figure 3, un exemple de collecte de mesures, effectuées par une ou plusieurs station RIMS, en un ensemble ou nuage de points IPP 20 contenus dans une zone d'influence ou voisinage 42 d'un IGP prédéterminé d'un réseau maillé ionosphérique 44, est illustré. Ici, le maillage illustré est un maillage de points IGP, défini par le standard RTCA-D0229, toutes versions comprises jusqu'à la version E actuelle, et le standard ICAO SARPs, avec des mailles de forme carrée. Dans cet exemple, le point IGP 46 dont la zone 25 d'influence 42 est illustrée est entouré de huit points IGP immédiatement voisins 48, 50, 52, 54, 46, 58, 60, 62 et forme le point commun de quatre mailles adjacentes 64, 66, 66, 68 qui l'entourent et constituent l'essentiel de sa zone d'influence 42. Les points IPP d'observation, représentés par des étoiles et situés dans la zone d'influence 42 définie par un réglage dans 30 l'algorithme, portent et sont associés chacun une mesure de VTEC. Ainsi, chaque IPP situé dans cette zone d'influence contribue à l'élaboration du GIVD de l'IGP considéré avec des poids variables selon la qualité de leur mesure et leur distance par rapport à l'IGP concerné. Cette même opération est répétée sur tous les IGP de la zone de 35 service du SBAS ayant suffisamment d'IPP dans leur zone d'influence pour 3061560 8 estimer correctement une valeur de GIVD et un écart-type appelé Basic GIVE ou GIVE de départ. Une fois les GIVD obtenus pour chacun des IGP de la zone SBAS de service, suivant la Figure 4 on calcule l'UIVD (en anglais « User lonospheric 5 Vertical Delay » sur chaque IPP situé dans les facettes ou cellules ou mailles contenant l'IPP que l'on cherche à évaluer, par exemple l'IPP 76, en interpolant le GIVD des quatre IGP, ici désignés par 78, 80, 82, 84, définissant la facette dans laquelle l'interpolation peut être réalisée. Il est rappelé que le standard du SBAS utilisé impose la loi 10 d'interpolation qui utilise une méthode de transformation bilinéaire, pondérée en fonction de la distance des IGP intervenant par rapport à l'IPP. Ainsi, l'UIVD obtenu de l'IPP considéré permet d'observer l'écart obtenu entre le délai issu directement de la mesure sur l'IPP et l'UIVD issu de l'estimation du GIVD des IGP voisins et de leur interpolation. Une fois 15 l'UIVD obtenu, on le rabat avec la fonction de rabattement définie dans le standard pour obtenir l'UISD : UISD = Fpp standard- UIVD (équation #3) En effet, ici on se place dans les mêmes conditions qu'un utilisateur et on veut déterminer les erreurs qui vont être générées par l'utilisation de la loi 20 d'interpolation et la fonction de rabattement standardisées. Pour les bornes d'erreur ionosphériques GIVE et suivant le Figure 5, exactement la même méthode est appliquée. Une fois interpolé au niveau de l'IPP 76, le GIVE devient UIVE (en anglais User lonospheric Vertical Error). Ensuite, comme pour l'UIVD, l'UIVE est rabattu selon l'élévation de la 25 ligne de vue associée à l'IPP considéré pour obtenir l'UISE (en anglais User lonospheric Slanted Error) : UISE = Fpp standard- UIVE (équation #4) Pour la vérification de l'intégrité garantie au niveau de l'utilisateur, on calcule les délais et les bornes sous leur forme verticale et non rabattue. 30 L'erreur de rabattement sera intégrée à la mesure de l'IPP lors de sa conversion de STEC vers VTEC. Suivant la méthode actuelle de vérification d'intégrité, on vérifie si le GIVE de chaque IGP suffit bien à borner l'erreur mesurée à partir de chaque IPP collecté par une station RIMS et ayant hérité de son couple UIVD et 3061560 9 UIVE, obtenu par interpolation des GIVD et GIVE des IGP constituant la cellule dans laquelle l'interpolation à l'IPP a été réalisée. Pour cela, un IGP courant ayant été choisi parmi l'ensemble des IGP à parcourir, ici l'IGP désigné par la référence numérique 92, et suivant la 5 Figure 6, on considère l'ensemble des IPP rattachés à au moins une cellule contenant l'IGP considéré dans leur interpolation. Sur la Figure 6, l'ensemble 94 des IPP d'observation qui seront utilisés pour la vérification de l'intégrité GIVE assurée par le GIVE de l'IGP testé se trouvent à l'intérieur des quatre cellules 96, 98, 100, 102 représentées. L'un 10 des IPP est d'ailleurs mis en évidence avec son VTEC mesuré. Il servira par la suite d'exemple dans l'application de la méthode de contrôle d'intégrité. Pour chaque IPP, on calcule l'écart entre la valeur de délai ionosphérique vertical interpolée à son niveau, c'est-à-dire l'UIVD et la valeur de délai vertical directement mesuré VTEC sur cet IPP. 15 On rappelle que le VTEC est obtenu à partir du STEC mesuré que l'on multiplie par l'inverse de la fonction de rabattement : STEC; VTEC. = Fpp standard (équation #6) On vérifie alors que l'erreur mesurée est bien couverte par l'UIVE avec le niveau d'intégrité requis, c'est-à-dire par exemple 10-7, ce qui correspond à 20 5,33 a pour une distribution gaussienne, pour chaque IPPi , i désignant un indice de parcours des IPP contenus dans les cellules immédiatement voisines de l'IGP testé et définies au préalable. L'UIVE est interprété dans le standard comme une variance avec un niveau de confiance ramené à 99.9%, c'est-à-dire le carré de 3.29a. Dans 25 notre cas, nous considérons un écart-type, donc 1 a. Nous l'appellerons donc : auwE. Tout d'abord, une innovation de chaque IPPi désignée par stdUIVDerrori , est calculée suivant l'expression : 30 IVTECi - UIVDi I stdUIVDerrori = 6vTEci2 + 6invEi2 (équation #7) avec : VTEC;, le délai ionosphérique vertical mesuré à l'IPP 3061560 10 UIVDi, le délai ionosphérique vertical interpolé à partir des GIVD des IGP entourant l'IPP i ; -- avTEci, l'écart-type du bruit de mesure à l'IPP oeUIVEil'écart-type issu des GIVE des IGP ayant participé à 5 l'interpolation de l'UIVD de l'IPPi concerné. Il est à remarquer que les valeurs Cl et P2 servant à l'élaboration de la mesure de VTECi à l'IPPi sont bruitées et généralement filtrées en amont. la prise en compte du bruit résiduel est donc nécessaire pour limiter les fausses alarmes. 10 Parmi l'ensemble des innovations calculées, on regarde si la valeur maximum des obtenue dépasse la valeur de seuil prédéterminée 5.33 qui correspond au niveau d'intégrité de 10-7 : stdUIVDerrori > 5.33 15 Si l'innovation est supérieure au seuil autorisé, alors un essai de réajustement du GIVE de l'IGP testé est entrepris en faisant gonfler (en anglais « inflate ») le premier rapport stdUIVDerroripar un facteur de 5.33 réajustement, par exemple ici égal à 1,1 lorsqu'une marge de 10% est souhaitée, pour obtenir un deuxième rapport « gonflé », désigné par K et 20 ayant pour expression : stdUIVD error; K = 5.33 x 1.1 La variance du GIVE réajusté, désignée par aGIVE(K) , s'écrit donc : acIvE (K) = K. aGivE 25 Cependant, la principale limitation de cette opération réside dans le fait qu'elle n'assure pas une correction parfaitement ajustée car on ne prend pas ici en compte les autres IGP ayant servi à l'obtention de l'UIVDi de l'IPPi considéré. 30 Afin d'éviter tout risque de non respect du niveau d'intégrité, lorsque la condition stdUIVDerrori > 5.33 est réalisée, condition ici associée au niveau 3061560 11 d'intégrité de 99,99999%, soit un risque de 10-7, alors l'utilisation de l'IGP considéré est interdite sans essayer de corriger son GIVE. Cela a pour conséquence de potentiellement réduire la capacité d'un terminal utilisateur à utiliser l'ensemble des satellites qui lui sont visibles pour 5 réaliser, s'il s'agit d'un avion, son approche vers la piste d'atterrissage et donc d'accroître le risque d'empêcher le terminal utilisateur d'atteindre le niveau de service nécessaire pour effectuer une tâche prédéterminée, une approche précise dans le cas d'un avion. De plus, si cet empêchement survient pendant une approche précise 10 dans le cas d'un avion, ce dernier n'aura de choix alors que d'arrêter son approche. L'avion remet alors les gaz, s'éloigne de la piste et se met en attente au-dessus de l'aéroport en attendant que la tour de contrôle lui propose une autre procédure d'approche, par exemple une approche dite « non précise ». Cela pourrait conduire à une gestion plus complexe du trafic 15 aérien au-dessus de la zone impactée. Une première invention vise à remédier aux inconvénients décrits ci-dessus. Le premier problème technique que résout la première invention est de proposer un procédé d'ajustement optimal des bornes d'erreurs de GIVE 20 de points IGP d'une grille ionosphérique de correction d'une zone de service d'un système SBAS, et un système de mise en oeuvre dudit procédé, qui garantissent à la fois la tenue de l'exigence d'intégrité d'un service SBAS prédéterminé et l'utilisation d'un nombre maximal d'IGP ayant leur GIVE ajustés pour une meilleure disponibilité du service SBAS. 25 A cet effet, la première l'invention a pour premier objet un procédé d'ajustement optimal de bornes d'erreurs de correction GIVE d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système d'augmentation à base de satellites SBAS, la grille de correction ionosphérique étant structurée en un réseau maillé de cellules 30 de forme polygonale prédéterminée, les cellules du réseau maillé correspondant projectivement à la zone de service SBAS et ayant pour sommets des point IGP du premier ensemble. Le procédé d'ajustement optimal des bornes d'erreur comprend les étapes consistant à : .- calculer par au moins un calculateur électronique, à partir d'informations 35 prédéterminées de corrections d'erreurs ionosphériques des points IGP du 3061560 12 premier ensemble, et de mesures de pseudo-distances d'un deuxième ensemble de points de percement IPP d'observation et de contrôle contenus dans les cellules de la grille ionosphérique, pour chaque point de percement d'observation IPPi du deuxième ensemble, une innovation désignée par 5 stdUIVDerrori , suivant l'expression : stdUIVDerrori = IVTEC; - UIVD; I icrvrEci2+ cruivEi2 dans laquelle : .* VTECi désigne le délai ionosphérique vertical mesuré à l'IPPi, . UIVDi désigne le délai ionosphérique vertical interpolé à partir des GIVDj des IGPj de la maille de rang m entourant l'IPPi concerné; 10 -* avTEci désigne l'écart-type du bruit de mesure à l'IPPi, * crtnvEi désigne l'écart-type issu des GIVE des IGPj de la maille de - rang m ayant participé à l'interpolation de l'UIVDi de l'IPPi concerné ; puis .- lorsque l'innovation de l'IPPi parcouru dans la maille actuelle de rang m est strictement supérieur à un seuil théorique MaxThd correspondant à niveau 15 de confiance et d'intégrité prédéterminé requis par le service SBAS, calculer un incrément de variance AuwEi2 qu'il faudrait rajouter au 6uwEi2 sur l'IPPi considéré pour qu'il recouvre le niveau d'intégrité requis suivant l'équation : AuivEi2= (3+ 6uivEi2) X (Kfact2 - 1) dans laquelle : 206uwEi désigne l'écart-type issu des GIVE des IGPj de la maille de rang m ayant participé à l'interpolation de l'UIVDi de l'IPPi concerné ; F3 et Kfact2 désignent respectivement un premier terme et un deuxième terme, Le premier terme F3 étant déterminé selon l'équation: = CYVTECi2 25 où VTECi désigne le délai ionosphérique vertical mesuré à l'IPPi considéré, et le deuxième terme Kfact2 étant déterminé selon l'équation : Kfact2 = (SafMarg x MaxThd)2 .SafMarg désignant la marge d'intégrité en pourcentage du seuil 30 maximum toléré que l'on souhaite générer, supérieur à 0 et strictement stdUIVDerrori 2 3061560 13 inférieur à 1, et configurée en fonction d'une marge de garantie d'intégrité en valeur relative prédéterminée, notée X et exprimée en pourcentage, suivant l'expression : SafMarg = (1 - X). Le procédé d'ajustement optimal étant caractérisé en ce qu'il 5 comprend en outre les étapes consistant à, lorsque l'innovation de l'IPPi parcouru dans la maille actuelle de rang m est strictement supérieur à un seuil théorique MaxThd : .- déterminer des valeurs d'incréments de variance AGivEk2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré par une méthode 10 d'interpolation inverse utilisant une méthode de Moindres Carrés de sorte à répartir l'incrément de variance AuivEi2 de l'IPPi suivant la relation : NIGP àUIVEi2 = wk- àGIVEk2 k=1 dans laquelle NIGP est le nombre d'IGPs de la maille de forme prédéterminée utilisés par le calcul de l'UIVDi et UIVEi de l'IPPi considéré ; 15 k désigne un indice de numérotation dans la maille contenant l'IPPi des IGP de ladite maille, wk désignent les poids respectifs des IGPk, k variant de 1 à NiGp, calculés par application du standard GNSS suivant une fonction du standard GNSS qui dépend de la distance entre l'IPPi et l'IGPk , la somme des poids 20 wk étant égale à un ; puis .- pour chaque IGPk de la maille m à laquelle l'IPPi considéré, mettre à jour la variance aGivEk dudit IGPk en remplaçant la valeur actuelle crGivEkdu GIVEk par une nouvelle valeur égale à _/ GGIVEk2 AGivEk2 - Suivant des modes particuliers de réalisation, le procédé d'ajustement 25 optimal des bornes d'erreurs de correction ionosphérique GIVE comprend l'une ou plusieurs des caractéristiques suivantes : .- l'étape de détermination des valeurs d'incréments de variance AGivEk2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré est mise en oeuvre par une méthode de Moindre Carré Classique dans 30 laquelle un vecteur Yi de répartition du AuivEi2 de l'IPPi considéré sur les IGPk associés est calculé suivant l'équation : Yi = Hit. (Hi. Hit)-1. -uivEi2 dans laquelle : 3061560 14 -* Hi = [wl. w2 --- wNIGP] désigne un vecteur-ligne à NIGP composantes pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk, k variant de 1 à NIGP, correspondant au poids d'un IGPk obtenu avec la méthode de calcul d'interpolation directe définie par le standard pour l'IPPi considéré ; AGivEi2 2 5 yi AGIVE2 désigne un vecteur-colonne à NIGP composantes de la GIVENIGP 2 répartition de l'incrément de variance 41\7E2 à répartir sur les NIGP IGPk lorsque la maille de l'IPPi considéré comporte NIGP IGPs, comme résultat de la résolution de l'équation ; -* AUIVEi2 désigne l'incrément de variance à répartir ; 10 .- l'étape de détermination des valeurs d'incréments de variance AGIVEk 2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré est mise en oeuvre par une méthode de moindre carré ajustée par un niveau de confiance 1-p dans laquelle un vecteur Yi de répartition du AuwEi2 de l'IPPi considéré sur les IGPk associés est calculé suivant l'équation : 15 Yi = TiCALS Q. Hi t. (Hi. Hit)- 1. A -uivEi2 dans laquelle : 3,GIVE12 2 .* Yi = AGIVE2 désigne un vecteur-colonne à NIGP composantes de la A GIVENIGP 2 répartition de l'incrément de variance AuivEi2 à répartir sur les NIGP IGPk lorsque la maille de l'IPPi considéré comporte NIGP IGP, comme résultat de la résolution de l'équation ; 20 -* Hi = [wl. w2 --- wNIGP] désigne un vecteur-ligne à NIGP composantes pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk, k variant de 1 à NIGP, correspondant au poids d'un IGPk obtenu avec la méthode de calcul d'interpolation directe définie par le standard pour l'IPPi considéré ; et .* TCALSQ I 25 tv,,a2 0 0 TicALsQ = 1 0 tv2,0,2 0 G(p)2 0 0 tv4,a2 avec tvk,OE désignant le facteur de Student associé l'IGPk de rang k fonction de vk le degré de liberté et fonction de la probabilité de confiance a , ledit est une matrice carré d'inflation de rang NIGP définie par : 3061560 15 facteur de Student tvk,a permettant de borner une erreur avec un niveau de confiance (1-p) prédéterminé et requis par le service SBAS, la probabilité de confiance a étant lié au niveau de confiance par la relation a = 1- p , et le degré de liberté vk étant égal ici au nombre d'IPP ayant généré un résidu 5 pour l'IGPk concerné moins un ; et G (p) désignant la valeur obtenue lorsque l'hypothèse gaussienne est appliquée pour la même probabilité de confiance (1-p) que celle utilisée pour le calcul de a, c'est-à-dire la limite du facteur de Student tyk,c, lorsque le degré de liberté vk tend vers l'infini ; et -* AUIVE; 2 désigne l'incrément de variance à répartir ; 10 .- l'étape de détermination des valeurs d'incréments de variance A G IVEk 2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré est mise en oeuvre par une méthode de moindre carré pondéré par un niveau de confiance niv_conf dans laquelle un vecteur Yi de répartition du AuivEi2 de l'IPPi considéré sur les IGPk associés est calculé suivant l'équation : 15 Yi = Pi. Hit. (Hi. Pi. Hit)-1. AuivE,2 dans laquelle : AGIVE12 2 Yi = AGIVE2 [ désigne un vecteur-colonne à NiGp composantes de la AGIVENIGP 2 répartition de l'incrément de variance AuwEi2 à répartir sur les NiGp IGPk lorsque la maille de l'IPPi considéré comporte NiGp IGP, comme résultat de la résolution de l'équation ; 20 -* Hi = [w1 w2 --- wNIGPI désigne un vecteur-ligne à NiGp composantes pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk, k variant de 1 à NiGp, correspondant au poids d'un IGPk obtenu avec la méthode de calcul d'interpolation directe définie par le standard pour l'IPPi considéré ; et .* Pi est une matrice carré de pondération diagonale de rang NiGp définie 25 par : Pi = (;11. 0 CINIGP 0 q2 0 0 t 2 avec qk = (e'v associé à l'IGPk, k variant de 1 à NiGp, t vk,a désignant le facteur de Student associé l'IGPk de rang k fonction de vk le degré de liberté 3061560 16 et fonction de la probabilité de confiance a , ledit facteur de Student tyk,« permettant de borner une erreur avec un niveau de confiance (1-p) prédéterminé et requis par le service SBAS, la probabilité de confiance a étant lié au niveau de confiance par la relation a = 1- , et le degré de 5 liberté vk étant égal ici au nombre d'IPP ayant généré un résidu pour l'IGPk concerné moins un ; et G(p) désignant la valeur obtenue lorsque l'hypothèse gaussienne est appliquée pour la même probabilité de confiance (1-p) que celle utilisée pour le calcul de a, c'est à dire la limite du facteur de Student tyk,,,, lorsque le degré de liberté vk tend vers l'infini ; et 10 .* AUIVEi 2 désigne l'incrément de variance à répartir ; .- l'étape de détermination des valeurs d'incréments de variance AGivEk2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré est mise en oeuvre par une méthode de moindre carré pondéré par un niveau de confiance 1-p dans laquelle un vecteur Yi de répartition du AinvEi2 de 15 l'IPPi considéré sur les IGPk associés est calculé suivant l'équation : 2 (tNbiPP MoyEN,a)2 Yi = Hit. (Hi. Hit)-1. A -UIVEi - G(p) dans laquelle : AGIvEi2 2 .* Yi = AcivE2 désigne un vecteur-colonne à NiGp composantes de la AGIVENIOP 2 répartition de l'incrément de variance A lorsque la maille de l'IPPi considéré comporte NiGp IGP, comme résultat de 20 la résolution de l'équation ; .* Hi = [w1 w2 --- WNIGP] désigne un vecteur-ligne à NiGp composantes pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk, k variant de 1 à NIGP, correspondant au poids d'un IGPk obtenu avec la méthode de calcul d'interpolation directe définie par le standard pour l'IPPi considéré ; et 25 -* tNbIPP MoyEN,a désigne un facteur de Student fonction de la probabilité de confiance a , ledit facteur de Student tN,,IpP MOYEN, c, permettant de borner une erreur avec un niveau de confiance (1-p) prédéterminé et requis par le service SBAS, la probabilité de confiance a étant lié au niveau de confiance par la relation a = 1 , et le degré de liberté Nbipp MOYEN étant un 30 nombre calculé à partir de la somme pondérée du nombre d'IPP au à répartir sur les NiGp IGPk 3061560 17 voisinage des IGPk ayant servis à l'obtention de l'UIVDi et UIVEi de l'IPPi testé, suivant la formule : vk=1 NIGP (Nbipp k - 1) Zi NbIPP MOYEN - NIGP Lik=1 W v k avec : 5 .* wk, le poids de l'IGPk obtenu par application de la méthode d'interpolation définie par le standard RTCA DO-229D sur l'IPPi concerné pour l'obtention de son UIVDi et UIVEi ; .* Nbipp k correspond au nombre d'IPP situé autour de l'IGPk dans le voisinage constitué de la réunion de toutes les cellules contenant cet IGPk ; 10 et .* G (p) désignant la valeur obtenue lorsque l'hypothèse gaussienne est appliquée pour la même probabilité de confiance (1-p) que celle utilisée pour a, c'est-à-dire la limite du facteur de Student tvivo, lorsque le degré de liberté vk tend vers l'infini ; et 15 -* AuniE12 désigne l'incrément de variance à répartir ; .- la forme de chaque maille est une forme triangulaire et le nombre NiGp des IGP sommets de chaque maille et égale à 3, ou la forme de chaque maille est la forme d'un quadrilatère, de préférence comprise parmi les formes d'un rectangle, d'un carré et d'un losange, et le nombre NiGp des 20 IGP sommets de chaque maille et égale à 4. La première invention a également pour deuxième objet un système SBAS d'augmentation à base de satellites pour augmenter les performances d'un système global de navigation par satellites GNSS, le système SBAS comprenant des terminaux utilisateurs de service(s) SBAS, un ou plusieurs 25 satellites d'augmentation des satellites du système GNSS et de diffusion de messages d'information aux terminaux utilisateurs, une ou plusieurs stations d'observation RIMS munis de récepteurs GNSS, et un ou plusieurs ou calculateurs. Le système SBAS est configuré pour ajuster de manière optimale des bornes d'erreurs de correction ionosphérique, appelés GIVE 30 finaux, d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service du système SBAS, la grille de correction ionosphérique étant structurée en un réseau maillé de cellules de forme polygonale prédéterminée, les cellules du réseau maillé correspondant 3061560 18 projectivement à la zone de service SBAS et ayant pour sommets des point IGP du premier ensemble. Le système SBAS est caractérisé en ce que le ou les calculateurs électroniques sont configurés pour : à partir d'informations prédéterminées de correction d'erreur 5 ionosphériques des points IGP du premier ensemble, et de mesures de pseudo-distances d'un deuxième ensemble de points de percement IPP d'observation et de contrôle contenus dans les cellules de la grille ionosphérique, pour chaque point de percement d'observation IPPi du deuxième ensemble, une innovation désignée par stdUIVDerrori , suivant 10 l'expression : stdUIVDerrori = 1VTECi - UIVDil VCrVTECi2 CrUIVEi2 dans laquelle : .*VTECi désigne le délai ionosphérique vertical mesuré à l'IPP UIVD; désigne le délai ionosphérique vertical interpolé à partir des GIVDj des IGPj de la maille de rang m entourant l'IPPi concerné; 15 .* GwEci désigne l'écart-type du bruit de mesure à l'IPPi, -* CrUIVEi désigne l'écart-type issu des GIVE des IGPj de la maille de rang m ayant participé à l'interpolation de l'UIVDi de l'IPPi concerné ; puis .- lorsque l'innovation de l'IPPi parcouru dans la maille actuelle de rang m est strictement supérieur à un seuil théorique MaxThd correspondant à un 20 niveau de confiance ou d'intégrité (1-p) prédéterminé exigé par le service SBAS, calculer un incrément de variance AuivEi2 qu'il faudrait rajouter au CrUIVEi 2 sur l'IPPi considéré pour qu'il recouvre le niveau d'intégrité requis suivant l'équation : AuivEi 2 = 6UWEj2) X (Kfact2 - 1) dans laquelle : -* aUIVEi désigne l'écart-type issu des GIVE des IGPj de la maille de rang m 25 ayant participé à l'interpolation de l'UIVDi de l'IPPi concerné ; .* 13 et Kfact2 désignent respectivement un premier terme et un deuxième terme, le premier terme f3 étant déterminé selon l'équation: 13 = avTEci2 où VTECi désigne le délai ionosphérique vertical mesuré à l'IPPi considéré, etle deuxième terme Kfact2 étant déterminé selon l'équation : stdUIVDerrori2 Kfact2= (SafMarg x MaxThd)2 30 .*SafMarg désignant la marge d'intégrité en pourcentage du seuil maximum toléré que l'on souhaite générer, supérieure à 0 et strictement inférieure à 1, et configurée en fonction d'une marge de garantie d'intégrité en valeur 3061560 19 relative prédéterminée, notée X et exprimée en pourcentage, suivant l'expression : SafMarg = (1 - X) ; puis lorsque l'innovation de l'IPPi parcouru dans la maille actuelle de rang m est strictement supérieur à un seuil théorique MaxThd : 5 . - déterminer des valeurs d'incréments de variance AGivEk2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré par une méthode d'interpolation inverse utilisant une méthode de Moindres Carrés de sorte à répartir l'incrément de variance àuivEi2 de l'IPPi suivant la relation : NIGP UIVE 2 = Wk. AGIVEk2 k=1 dans laquelle 10 .* NIGP est le nombre d'IGPs de la maille de forme prédéterminée utilisés par le calcul de l'UIVDi et UIVEi de l'IPPi considéré ; .* k désigne un indice de numérotation dans la maille contenant l'IPPi des IGP de ladite maille, .* wk désignent les poids respectifs des IGPk, k variant de 1 à NIGP, calculés 15 par application du standard GNSS suivant une fonction du standard GNSS qui dépend de la distance entre l'IPPi et l'IGPk , la somme des poids wk étant égale à un ; puis .- pour chaque IGPk de la maille m à laquelle appartient l'IPPi considéré, mettre à jour la variance GIVEk dudit IGPk en remplaçant la valeur actuelle 20 aGivEkdu GIVEk par une nouvelle valeur égale à -la v/- GIVEk2 AGIvEk2 Une deuxième invention concerne un deuxième procédé de calcul optimal des écarts-type appelés o-mops qui permettant de modéliser l'erreur d'interpolation, réalisée entre les IGP et le point de percement du signal reçu par l'utilisateur, ainsi que l'erreur de la fonction de rabattement du délai 25 vertical (en anglais « vertical delay ») vers le délai oblique (en anglais « slant delay »), ladite fonction de rabattement dépendant de l'élévation satellite et de l'antenne du terminal GNSS). Les écarts-type 0MOPS interviennent dans le calcul des GIVE finaux et constituent un des termes ajoutés aux Basic GIVE pour l'obtention du GIVE final. Ainsi l'écart-type o-mops contribue t'il à conférer 30 à la borne d'erreur finale assurant l'intégrité, c'est-à-dire le GIVE final, la confiance requise pour être conforme à la sécurité élevée de services SBAS prédéterminés comme ceux relatifs aux opérations aériennes précises. Ce terme est particulièrement important car il permet de recentrer l'erreur 3061560 20 obtenue du fait des limitations induites par les imprécisions provoquées par les lois de calcul imposées par le standard. Suivant la méthode actuelle de calcul des écarts-type appelés YMOPS, pour chaque point IGP un voisinage V de points de percement entourant 5 l'IGP est défini. Pour chaque point de percement de mesure IPPi contenu dans ce voisinage V, un résidu noté Resi , ramené à la verticale est calculé suivant l'équation : STECi - UISDi Resi = Fpp standard (équation #8) 10 dans lequel : ^ STECi désigne le délai ionosphérique de la ligne de vue réelle, donc fonction de l'élévation ; ^ UISDi désigne le délai ionosphérique vertical interpolé suivant le standard à partir des GIVDk des IGPk de la maille ou cellule entourant 15 l'IPPi ; -- Fpp standard désigne la fonction de rabattement standard du délai ionosphérique en fonction de l'élévation de la ligne de vue. Le résidu Resi est ramené à la verticale en raison du fait que les délais et leur borne d'erreur au niveau des IGPk sont toujours formulés en vertical 20 pour permettre à l'utilisateur de les adapter selon l'élévation de leurs propres lignes de vues. Le résidu Resi du point de percement de mesure IPPi est ensuite pondéré suivant la qualité de sa mesure et la distance de l'IPPi par rapport à l'IGP considéré : 25 pi = f(6vTEc2, Distance) (équation #9) Une fois le résidu Resi obtenu pour chaque IPPi situé dans le voisinage V contenant l'IGP considéré, une erreur moyenne, issue des différents résidus, est évaluée. Cette valeur d'erreur peut être calculée de 30 plusieurs façons différentes, le but étant d'obtenir une valeur la plus représentative possible du biais généré par l'utilisation des méthodes d'interpolation et de rabattement définies par le standard afin de « re- 3061560 21 centrer » la borne d'erreur CIVE et garantir le niveau d'intégrité requis suivant une loi gaussienne centrée. Dans le système EGNOS (en anglais « European Geostationary Navigation Service »), ce terme est un écart-type désigné et appelé aMOPS (l'acronyme MOPS étant le standard définissant le système EGNOS). Cette valeur d'écart-type ou de RMS (en anglais Root Mean Square) est calculée comme étant la racine carrée de la moyenne des carrés des résidus pondérés suivant l'équation : Nipp OEMOPS 7, X - Resi2 P Pi i=1 10 (équation #10) avec : . P égal à la somme des poids normalisés suivant la formule : P = EN. i" !i où Nipp désigne le nombre de point de percement de mesure 1=1 p IPPi contenu dans le voisinage V entourant l'IGP considéré. 15 La principale limitation liée à cette méthode de calcul actuellement utilisée est la mise en place de la pondération, car chaque IPPi contribue différemment en fonction de la distance d'éloignement de l'IGP et de la qualité de la mesure. De plus, étant donné qu'au moins trois IGP, quatre dans le cas d'une maille rectangulaire ou carrée, sont en réalité intervenus 20 dans le calcul de l'UISDi de chaque IPPi, de l'information a été perdue et on ne maîtrise pas parfaitement la quantité d'erreur que l'on injecte dans le résidu une fois pondéré. Ensuite, le aMOPS est dilaté suivant le nombre d'IPPi ayant contribué à son estimation par un facteur de Student permettant de définir l'indice de 25 confiance de la valeur mesurée par rapport à une hypothèse gaussienne : (tvk,05.33 .99999995) OEMOPS X aMOPS avec : 3061560 22 tyk,0.999999955 le facteur de Student en fonction de vk, le degré de liberté (ici c'est le nombre d'IPP) permettant de borner une erreur avec une probabilité de 10-7 requise par le standard aéronautique ; et .- 5.33 qui est la valeur obtenue lorsque l'hypothèse gaussienne est 5 appliquée, pour la même probabilité que celle définie ci-dessus, c'est à dire 0,9999999 de confiance. La deuxième invention vise à remédier aux inconvénients précités, en particulier la principale limitation mentionnée. Le deuxième problème technique que résout la deuxième invention 10 est de proposer un procédé d'ajustement optimal des écarts-type crMOps de points IGP d'une grille ionosphérique de correction d'une zone de service d'un système SBAS, et un système de mise en oeuvre dudit procédé, qui soient indépendants de la configuration ou de la répartition spatiale des points de percement de mesure IPP situés autour des IGP d'interpolation. 15 A cet effet, la deuxième invention a pour troisième objet un procédé de calcul optimal des variances des résidus d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système SBAS, la grille de correction ionosphérique étant structurée en un réseau maillé de cellules de forme polygonale prédéterminée, les cellules du 20 réseau maillé correspondant projectivement à la zone de service SBAS et ayant pour sommets des point IGP du premier ensemble. Le procédé de calcul optimal des variances des résidus des IGP comprend les étapes consistant à : .- calculer par au moins un calculateur électronique, à partir d'informations 25 prédéterminées de correction d'erreur ionosphériques des points IGP du premier ensemble, et de mesures de pseudo-distances d'un deuxième ensemble de points de percement IPP d'observation et de contrôle contenus dans les cellules de la grille ionosphérique, pour chaque point de percement d'observation IPPi du deuxième ensemble, un résidu, désigné par Resi , 30 ramené à la verticale, suivant l'équation : STECi - UISDi Resi = r pp standard dans laquelle : STECi désigne le délai ionosphérique de la ligne de vue réelle mesuré par une station d'observation et fonction de l'élévation de ladite ligne de vue, UISDi désigne le délai ionosphérique vertical interpolé suivant le standard du terminal utilisateur à partir des GIVDk des IGPk de la maille 3061560 23 entourant l'IPPi , et Fpp standard désigne la fonction de rabattement standard du délai ionosphérique en fonction de l'élévation de la ligne de vue, puis le carré du résidu . Le procédé de calcul optimal des variances des résidus des IGP est caractérisé en ce qu'il comprend en outre les étapes consistant à : 5 pour chaque point de percement d'observation IPPi du deuxième ensemble, .- déterminer des valeurs d'incréments de variance A -Resk2 à attribuer respectivement aux IGPk de la maille m de l'IPPi considéré par une méthode d'interpolation inverse utilisant une méthode de Moindres Carrés dans laquelle un vecteur Yi de répartition du carré du résidu Resi2 de l'IPPi 10 considéré est calculé suivant l'équation Yi = Hit. (Hi. Hit)-1. Resi2 dans laquelle : .* H = [w1 w2 wNiGp] désigne un vecteur-ligne à NiGp composantes pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk , k variant de 1 à NIGP , correspondant au poids d'un IGPk obtenu avec la méthode de 15 calcul d'interpolation directe définie par le standard pour l'IPPi considéré ; ARes12 2 Yi = ARes2 désigne un vecteur-colonne à NiGp composantes de la AResNIGP2 répartition du carré du résidu Resi2 à répartir sur les NiGp IGPk lorsque la maille de l'IPPi considéré comporte NiGp IGP, comme résultat de la résolution de l'équation ; 20 .* Res? désigne le carré du résidu à répartir, puis .- déterminer un coefficient de pondération pi comme inverse d'un poids, ce coefficient de pondération étant représentatif de la qualité de la mesure du STECi redressé en VTECi par la division avec la fonction de rabattement de l'IPPi considéré, et exprimé suivant la formule suivante 25 pi = a VTEC12 2 dans laquelle CrVTEC NOMINAL -* aVTECi2 désigne le bruit de mesure de la pseudo-distance mesuré par le terminal de mesure pour l'IPPi considéré, et -* aVTECi NOMINAL2 désigne un bruit de mesure nominal de référence du terminal de mesure ; 30 .- ensuite, après avoir calculé tous les Yi et pi correspondant au premier ensemble des points de percement IPPi observé, pour chaque IGP, calculer une variance de résidu non dilaté suivant l'équation : 3061560 24 aMOPS = Nipp 1 1 2 rs X - X AResipp i p i=1 avec : .* Ps, la somme des poids normalisés 1 /pi de chaque résidu disponibles pour l'IGP considéré suivant : Ps = p p 1pi .* Nipp étant le nombre des IPP ayant généré un résidu pour l'IGP concerné, 5 -* AResipp i 2 le résidu interpolé en inverse à l'IGP considéré pour l'IPPi faisant partie des IPP ayant généré un résidu pour l'IGP concerné. Suivant des modes particuliers de réalisation, le procédé de calcul optimal des variances des résidus d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système 10 SBAS comprend l'une ou plusieurs des caractéristiques suivantes : .- le coefficient de pondération pi est modulé par l'inverse du poids wi de l'IPPi suivant l'IGP tel que défini par le standard MOPS DO-229 et s'écrit 1 ci i OEVTECi 2 VTEC NOMINAL suivant l'expression : pi = w- x 2 , et l'écart-type aMOPS de l'IGP considéré est calculé suivant l'expression 15 MOPS = j(K. P. Kt)-1. K. P. Xt dans laquelle : .* K désigne le vecteur unitaire ligne de dimension Nipp K = [1 1 ... 1], 1 0 - Pi .* P désigne la matrice diagonale des poids P = : : , et 1 PNipp .* X désigne le vecteur colonne des résidus : X = I AResNipP ARes 12 2 20 .- le coefficient de pondération pi est modulé par l'inverse du poids wi center norm de l'IPPi suivant l'IGP l'expression : 1 arrEci2 Pi = dans laquelle le poids wi center norm wi center norm ciVTEC NOMINAL2 est calculé selon l'équation : Si w1 = O, alors l'IPP n'est pas utilisé pour le calcul du o-mops de l'IGP 25 considéré Si O MaxThd Par exemple, MaxThd, est pris égal à 5,33 sous des hypothèses de distribution gaussienne de l'erreur, ce qui correspond à un niveau de -\/aVTECi2 aUIVEi2 3061560 34 confiance exigé par le service SBAS égal à (1-10-7), soit une confiance de 99,99999%. Ce seuil théorique pourra être modifié si on souhaite modifier le niveau de probabilité requis et sera noté par la suite G (p) avec (p) = er f -1(1- x , erf-10 désignant la fonction réciproque de la fonction 5 erreur et (1-p) désignant le niveau de confiance requis. Il est à noter ici que le seuil maximum autorisé peut être calculé avec une marge, évitant ainsi les cas jugés trop limites. Si cette condition n'est pas vérifiée, c'est-à-dire si l'innovation de l'IPPi actuel est inférieure ou égale au seuil théorique prédéterminé MaxThd , 10 cela veut dire que les IGPj de la maille actuelle de rang m ayant servi à l'interpolation de l'UIVDi de l'IPPi considéré ont un GIVEj suffisant pour garantir le niveau d'intégrité requis. Dans ce cas, les GIVEj de la maille actuelle restent inchangés et une quatrième étape 216 de test de la fin de la deuxième boucle locale de parcours des IPPi contenus dans la maille de 15 rang m actuel est exécutée. Si cette condition est vérifiée, c'est-à-dire si l'innovation de l'IPPi actuel dépasse le seuil MaxThd, alors, la quantité unrEi qu'il faudrait rajouter au auwEi sur l'IPPi considéré pour qu'il recouvre le niveau d'intégrité requis est calculée dans une cinquième étape 218 suivant l'équation : 20 AUIVEi 2 = (13 + CTUIVEi x (Kfact2 - 1) (équation #12) dans laquelle : -- auivEi désigne l'écart-type issu des GIVE des IGPj de la maille de rang m ayant participé à l'interpolation de l'UIVDi de l'IPPi concerné ; 25 .- 0 et Kfact2 désignent respectivement un premier terme et un deuxième terme. Le premier terme s'écrit selon l'équation: = aVTECi 2 où VTECi désigne le délai ionosphérique vertical mesuré à l'IPPi considéré. 30 le deuxième terme Kfact2 s'écrit selon l'équation : 1 REVENDICATIONS1. Procédé d'ajustement optimal de bornes d'erreurs de correction 5 ionosphérique GIVE d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système d'augmentation à base de satellites SBAS, la grille de correction ionosphérique étant structurée en un réseau maillé de cellules de forme polygonale prédéterminée, les cellules du réseau maillé correspondant 10 projectivement à la zone de service SBAS et ayant pour sommets des point IGP du premier ensemble, le procédé d'ajustement optimal des bornes d'erreur comprenant les étapes consistant à : .- calculer (212) par au moins un calculateur électronique, à partir 15 d'informations prédéterminées de corrections d'erreurs ionosphériques des points IGP du premier ensemble, et de mesures de pseudo-distances d'un deuxième ensemble de points de percement IPP d'observation et de contrôle contenus dans les cellules de la grille ionosphérique, pour chaque point de percement d'observation IPPi du deuxième ensemble, une innovation 20 désignée par stdUIVDerrori , suivant l'expression : IVTECi - UIVDi I stdUIVDerrori = dans laquelle : .* VTECi désigne le délai ionosphérique vertical mesuré à l'IPPI, . UIVDi désigne le délai ionosphérique vertical interpolé à partir des GIVDj des IGPj de la maille de rang m entourant l'IPPi concerné; 25 .* OvrEci désigne l'écart-type du bruit de mesure à l'IPPi, -* cruivEi désigne l'écart-type issu des GIVE des IGPj de la maille de rang m ayant participé à l'interpolation de l'UIVDi de l'IPPi concerné ; puis .- lorsque l'innovation de l'IPPi parcouru dans la maille actuelle de rang m est strictement supérieur à un seuil théorique MaxThd correspondant 30 à niveau de confiance et d'intégrité prédéterminé requis par le service SBAS, calculer (218) un incrément de variance AuwEi2 qu'il faudrait rajouter au CrUIVEi 2 sur l'IPPi considéré pour qu'il recouvre le niveau d'intégrité requis suivant l'équation : \/CrVTECi2 CrUIVEi2 3061560 60 àUIVEi2 = (13 + 6UIVEj2) X (Kfact2 - 1) dans laquelle : -- auivE; désigne l'écart-type issu des GIVE des IGPj de la maille de rang m ayant participé à l'interpolation de l'UIVDi de l'IPPi concerné ; 13 et Kfact2 désignent respectivement un premier terme et un deuxième terme, Le premier terme p étant déterminé selon l'équation: 2 I3 = CiVTECi où VTECi désigne le délai ionosphérique vertical mesuré à l'IPPi considéré, et le deuxième terme Kfact2 étant déterminé selon l'équation : Kfaa 2 = (SafMarg x MaxThd)2 .SafMarg désignant la marge d'intégrité en pourcentage du seuil maximum toléré que l'on souhaite générer, supérieure à zéro et strictement inférieur à 1, et configurée en fonction d'une marge de garantie d'intégrité en valeur relative prédéterminée, notée X et exprimée en pourcentage, suivant l'expression : SafMarg = (1 - X) ; le procédé d'ajustement optimal étant caractérisé en ce qu'il comprend en outre les étapes (220, 222) consistant à, lorsque l'innovation de l'IPPi parcouru dans la maille actuelle de rang m est strictement supérieur à un 20 seuil théorique MaxThd : . - déterminer (220) des valeurs d'incréments de variance AGIvEk2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré par une méthode d'interpolation inverse utilisant une méthode de Moindres Carrés de sorte à répartir l'incrément de variance AuwEi2 de l'IPPi suivant la 25 relation : NIGP AutvEi2= >Wk. AGIVEk2 k=1 dans laquelle NIGP est le nombre d'IGPs de la maille de forme prédéterminée utilisés par le calcul de l'UIVDi et UIVEi de l'IPPi considéré ; stdUIVDerrori2 3061560 61 ^ k désigne un indice de numérotation dans la maille contenant l'IPPi des IGP de ladite maille, ^ wk désignent les poids respectifs des IGPk, k variant de 1 à NIGP, calculés par application du standard GNSS suivant une fonction du standard 5 GNSS qui dépend de la distance entre l'IPPi et l'IGPk , la somme des poids wk étant égale à un ; puis .- pour chaque IGPk de la maille m à laquelle l'IPPi considéré, mettre à jour (222) la variance o-GivEk dudit IGPk en remplaçant la valeur actuelle aGivEkdu GIVEk par une nouvelle valeur égale à _/ v/acivEk2+ AGIvEk2 10 2. Procédé d'ajustement optimal de bornes d'erreurs de correction ionosphérique GIVE selon la revendication 1, dans lequel l'étape de détermination (220) des valeurs d'incréments de variance AGivEk2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré est mise en oeuvre 15 par une méthode de Moindre Carré Classique dans laquelle un vecteur Yi de répartition du AuwEi2 de l'IPPi considéré sur les IGPk associés est calculé suivant l'équation : Yi Hit. (Hi. Hit)-f. A dans laquelle .- Hi = [w1 w2 WNIGp] désigne un vecteur-ligne à NIGP composantes pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk, k variant de 1 à NIGP, correspondant au poids d'un IGPk obtenu avec la méthode de calcul d'interpolation directe définie par le standard pour l'IPPi considéré ; AGRTE,2 Yi = '4IVE2 2 désigne un vecteur-colonne à NIGP composantes AGIVENIGP 2 de la répartition de l'incrément de variance AuwEi2 à répartir sur les NIGP IGPk lorsque la maille de l'IPPi considéré comporte NIGP IGPs, comme résultat de la résolution de l'équation ; AuivE,2 désigne l'incrément de variance à répartir. 3061560 62 3. Procédé d'ajustement optimal de bornes d'erreurs de correction ionosphérique GIVE selon la revendication 1, dans lequel l'étape de détermination (220) des valeurs d'incréments de variance AGIvEk2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré est mise en oeuvre 5 par une méthode de moindre carré ajustée par un niveau de confiance 1-p dans laquelle un vecteur Yi de répartition du AuivEi2 de l'IPPi considéré sur les IGPk associés est calculé suivant l'équation : Yi = TiCALSQ Hit. (Hi. Hit)-1. A dans laquelle AGIVE1 2 2 10 yi =GivE2 désigne un vecteur-colonne à 1\11Gp composantes AGIVENIGP 2 de la répartition de l'incrément de variance LluivEi2 à répartir sur les NIGP IGPk lorsque la maille de l'IPPi considéré comporte NiGp IGP, comme résultat de la résolution de l'équation ; Hi = w2 --- wNIGP] désigne un vecteur-ligne à NiGp composantes 15 pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk, k variant de 1 à NIGP, correspondant au poids d'un IGPk obtenu avec la méthode de calcul d'interpolation directe définie par le standard pour l'IPPi considéré ; et TCALSQ I 20 t 2 vi,a CALSQ = 1 0 tv2,a2 G(p)2 0 0 2 tv4,a avec tvk,a désignant le facteur de Student associé l'IGPk de rang k fonction de vk le degré de liberté et fonction de la probabilité de confiance a , ledit facteur de Student tvk,a permettant de borner une erreur avec un niveau de 25 confiance (1-p) prédéterminé et requis par le service SBAS, la probabilité de confiance a étant lié au niveau de confiance par la relation a = 1 - p , et le degré de liberté vk étant égal ici au nombre d'IPP ayant généré un résidu pour l'IGPk concerné moins un ; et est une matrice carré d'inflation de rang NiGp définie par : 0 0 o 3061560 63 G (p) désignant la valeur obtenue lorsque l'hypothèse gaussienne est appliquée pour la même probabilité de confiance (1-p) que celle utilisée pour le calcul de a, c'est-à-dire la limite du facteur de Student tyk,a lorsque le degré de liberté vk tend vers l'infini ; et 5 3iinvE12 désigne l'incrément de variance à répartir. 4. Procédé d'ajustement optimal de bornes d'erreurs de correction ionosphérique GIVE selon la revendication 1, dans lequel l'étape de détermination (220) des valeurs d'incréments de variance AGivEk2 à attribuer 10 respectivement aux IGPk de la maille de l'IPPi considéré est mise en oeuvre par une méthode de moindre carré pondéré par un niveau de confiance niv_conf dans laquelle un vecteur Yi de répartition du AuwEi2 de l'IPPi considéré sur les IGPk associés est calculé suivant l'équation : Yi = Pi. Hit. (Hi. Pi. Hit)-1. A -uivEi2 15 dans laquelle AGIvEi2 AGIVE22 désigne un vecteur-colonne à NiGp composantes AGIVENIGP 2 de la répartition de l'incrément de variance AuivE2 à répartir sur les NiGp IGPk lorsque la maille de l'IPPi considéré comporte NiGp IGP, comme 20 résultat de la résolution de l'équation ; Hi = [w1 w2 wNiGPi désigne un vecteur-ligne à NiGp composantes pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk, k variant de 1 à NIGP, correspondant au poids d'un IGPk obtenu avec la méthode de calcul d'interpolation directe définie par le standard pour l'IPPi 25 considéré ; et .- Pi est une matrice carré de pondération diagonale de rang NiGp définie par : q1 0 0 Pi= ° (12 0 0 0 CINIGP 3061560 64 l 2 avec qk = (-tvk associé à l'IGPk, k variant de 1 à NIGP, G(p) tvivec désignant le facteur de Student associé l'IGPk de rang k fonction de vk le degré de liberté et fonction de la probabilité de confiance a , ledit facteur de Student tvke permettant de borner une erreur avec un niveau de 5 confiance (1-p) prédéterminé et requis par le service SBAS, la probabilité de confiance a étant lié au niveau de confiance par la relation a = 1 - p , et le degré de liberté vk étant égal ici au nombre d'IPP ayant généré un résidu pour l'IGPk concerné moins un ; et G(p) désignant la valeur obtenue lorsque l'hypothèse gaussienne est 10 appliquée pour la même probabilité de confiance (1-p) que celle utilisée pour le calcul de a, c'est-à-dire la limite du facteur de Student tvio,' lorsque le degré de liberté vk tend vers l'infini ; et -- AuivEi2 désigne l'incrément de variance à répartir. 15 5. Procédé d'ajustement optimal de bornes d'erreurs de correction ionosphérique GIVE des IGP selon la revendication 1, dans lequel dans lequel l'étape de détermination (220) des valeurs d'incréments de variance AcivEk2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré est mise en oeuvre par une méthode de moindre carré pondéré par un niveau 20 de confiance 1-p dans laquelle un vecteur Yi de répartition du Au1vEi2 de l'IPPi considéré sur les IGPk associés est calculé suivant l'équation : 2 [ Yi = Hit. (Hi.Hityi. A rbiPP MoyEN,a) UIVEi 2 - G(p) dans laquelle AGIVE12 25 Yi = GivE22 désigne un vecteur-colonne à NiGp composantes LIGIvENIGp2 de la répartition de l'incrément de variance AuivEi2 à répartir sur les NiGp IGPk lorsque la maille de l'IPPi considéré comporte NiGp IGP, comme résultat de la résolution de l'équation ; Hi = [w1 w2 --- WNIGP] désigne un vecteur-ligne à NiGp composantes 30 pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk, k variant de 1 à NIGP, correspondant au poids d'un IGPk obtenu avec la 3061560 65 méthode de calcul d'interpolation directe définie par le standard pour l'IPPi considéré ; et tNbIPP MOYEN,C1 désignant un facteur de Student fonction de la probabilité de confiance a , ledit facteur de Student tNbIPP MOYEN, a permettant de borner 5 une erreur avec un niveau de confiance (1-p) prédéterminé et requis par le service SBAS, la probabilité de confiance a étant lié au niveau de confiance par la relation a = 1 p , et le degré de liberté Nbipp MOYEN étant un nombre calculé à partir de la somme pondérée du nombre d'IPP au voisinage des IGPk ayant servis à l'obtention de l'UIVDi et UIVEi de l'IPPi 10 testé, suivant la formule : NbIPP MOYEN = NIGP w, Zik=i K avec : wk, le poids de l'IGPk obtenu par application de la méthode d'interpolation définie par le standard RTCA DO-229D sur l'IPPi concerné 15 pour l'obtention de son UIVDi et UIVEi ; Nbipp k correspond au nombre d'IPP situé autour de l'IGPk dans le voisinage constitué de la réunion de toutes les cellules contenant cet IGPk ; et G (p) désignant la valeur obtenue lorsque l'hypothèse gaussienne est 20 appliquée pour la même probabilité de confiance (1-p) que celle utilisée pour le calcul de a, c'est-à-dire la limite du facteur de Student tvkx, lorsque le degré de liberté vk tend vers l'infini ; et AuivEi2 désigne l'incrément de variance à répartir. 25 6. Procédé d'ajustement optimal de bornes d'erreurs de correction ionosphérique GIVE selon l'une quelconque des revendications 1 à 5, dans lequel La forme de chaque maille est une forme triangulaire et le nombre NIGP des IGP sommets de chaque maille et égale à 3, ou 30 La forme de chaque maille est la forme d'un quadrilatère, de préférence comprise parmi les formes d'un rectangle, d'un carré et d'un losange, et le nombre NIGP des IGP sommets de chaque maille et égale à 4. EkN=IGiP wk. (Nbipp k - 1) 3061560 66 7. Système SBAS d'augmentation à base de satellites pour augmenter les performances d'un système global de navigation par satellites GNSS, le système SBAS comprenant des terminaux utilisateurs de service(s) SBAS, un ou plusieurs satellites d'augmentation des satellites du système 5 GNSS et de diffusion de messages d'information aux terminaux utilisateurs, une ou plusieurs stations d'observation RIMS munis de récepteurs GNSS, et un ou plusieurs ou calculateurs, le système SBAS étant configuré pour ajuster de manière optimale des bornes d'erreurs de correction ionosphérique, appelés GIVE finaux, d'un 10 premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service du système SBAS, la grille de correction ionosphérique étant structurée en un réseau maillé de cellules de forme polygonale prédéterminée, les cellules du réseau maillé correspondant projectivement à la zone de service SBAS et ayant pour sommets des point 15 IGP du premier ensemble, et le système SBAS étant caractérisé en ce que le ou les calculateurs électroniques sont configurés pour : .- à partir d'informations prédéterminées de correction d'erreur ionosphériques des points IGP du premier ensemble, et de mesures de 20 pseudo-distances d'un deuxième ensemble de points de percement IPP d'observation et de contrôle contenus dans les cellules de la grille ionosphérique, pour chaque point de percement d'observation IPPi du deuxième ensemble, une innovation désignée par stdUIVDerrori , suivant l'expression : stdUIVDerrori = IVTEC; - UIVDi I -NICrVTECi2 CrUIVE12 25 dans laquelle : .* VTECi désigne le délai ionosphérique vertical mesuré à l'IPP .* UIVD; désigne le délai ionosphérique vertical interpolé à partir des GIVDj des IGPj de la maille de rang m entourant l'IPPi concerné; o-vrEci désigne l'écart-type du bruit de mesure à l'IPPi, 30 -* CruivE, désigne l'écart-type issu des GIVE des IGPj de la maille de rang m ayant participé à l'interpolation de l'UlVDi de l'IPPi concerné ; puis .- lorsque l'innovation de l'IPPi parcouru dans la maille actuelle de rang m est strictement supérieur à un seuil théorique MaxThd correspondant 3061560 67 à un niveau de confiance ou d'intégrité (1-p) prédéterminé exigé par le service SBAS, calculer un incrément de variance AinvEi2 qu'il faudrait rajouter au auivEi2 sur l'IPPi considéré pour qu'il recouvre le niveau d'intégrité requis suivant l'équation : 5 AUIVEi2 - (P CiUIVEi 2) X (Kfact2 - 1) dans laquelle : aCJIVEi désigne l'écart-type issu des GIVE des IGPj de la maille de rang m ayant participé à l'interpolation de l'UIVDi de l'IPPi concerné ; R et Kfact2 désignent respectivement un premier terme et un 10 deuxième terme, Le premier terme fi étant déterminé selon l'équation: 13 = avTEci2 où VTECi désigne le délai ionosphérique vertical mesuré à l'IPPi considéré, et le deuxième terme Kfact2 étant déterminé selon l'équation : 15 Kfact2= (SafMarg x MaxThd)2 .SafMarg désignant la marge d'intégrité en pourcentage du seuil maximum toléré que l'on souhaite générer, supérieure à zéro et strictement inférieure à 1, et configurée en fonction d'une marge de garantie d'intégrité en valeur relative prédéterminée, notée X et exprimée en pourcentage, 20 suivant l'expression : SafMarg = (1 - X) ; puis lorsque l'innovation de l'IPPi parcouru dans la maille actuelle de rang m est strictement supérieur à un seuil théorique MaxThd : . - déterminer des valeurs d'incréments de variance LXGIVEk2 à attribuer respectivement aux IGPk de la maille de l'IPPi considéré par une méthode 25 d'interpolation inverse utilisant une méthode de Moindres Carrés de sorte à répartir l'incrément de variance AuwEi2 de l'IPPi suivant la relation : NIGP AuivEi2= >Wk- àGIVEk2 k=1 dans laquelle NiGp est le nombre d'IGPs de la maille de forme prédéterminée utilisés par le calcul de l'UIVDi et UIVEi de l'IPPi considéré ; stdUIVDerrori 2 3061560 68 k désigne un indice de numérotation dans la maille contenant l'IPPi des IGP de ladite maille, wk désignent les poids respectifs des IGPk, k variant de 1 à NIGP, calculés par application du standard GNSS suivant une fonction du standard 5 GNSS qui dépend de la distance entre l'IPPi et l'IGPk , la somme des poids wk étant égale à un ; puis .- pour chaque IGPk de la maille m à laquelle appartient l'IPPi considéré, mettre à jour la variance GIVEk dudit IGPk en remplaçant la valeur actuelle aGIVEkdu GIVEk par une nouvelle valeur égale à 10 8. Procédé de calcul optimal des variances des résidus d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système SBAS, la grille de correction ionosphérique étant structurée en un réseau maillé de cellules de forme polygonale 15 prédéterminée, les cellules du réseau maillé correspondant projectivement à la zone de service SBAS et ayant pour sommets des point IGP du premier ensemble, le procédé de calcul optimal des variances des résidus des IGP comprenant les étapes consistant à : 20 .- calculer (312) par au moins un calculateur électronique, à partir d'informations prédéterminées de correction d'erreur ionosphériques des points IGP du premier ensemble, et de mesures de pseudo-distances d'un deuxième ensemble de points de percement IPP d'observation et de contrôle contenus dans les cellules de la grille ionosphérique, pour chaque point de 25 percement d'observation IPPi du deuxième ensemble, un résidu, désigné par Resi , ramené à la verticale, suivant l'équation : Resi = STECi - UISDi Fpp standard dans lequel : STECi désigne le délai ionosphérique de la ligne de vue réelle 30 mesuré par une station d'observation et fonction de l'élévation de ladite ligne de vue; UISDi désigne le délai ionosphérique vertical interpolé suivant le standard du terminal utilisateur à partir des GIVDk des IGPk de la maille entourant l'IPPi ; ,./6GivEk2+ AGivEk2 3061560 69 Fpp standard désigne la fonction de rabattement standard du délai ionosphérique en fonction de l'élévation de la ligne de vue, puis le carré du résidu ; le procédé de calcul optimal des variances des résidus des IGP étant 5 caractérisé en ce qu'il comprend en outre les étapes consistant à, pour chaque point de percement d'observation IPPi du deuxième ensemble, . - déterminer (314) des valeurs d'incréments de variance A -Resk2 à attribuer respectivement aux IGPk de la maille m de l'IPPi considéré par une 10 méthode d'interpolation inverse utilisant une méthode de Moindres Carrés dans laquelle un vecteur Yi de répartition du carré du résidu Resi2 de l'IPPi considéré est calculé suivant l'équation : Yi = Hit. (Hi. Hit)-1. Resi2 15 dans laquelle -* H = [w1 w2 --- wNIGP] désigne un vecteur-ligne à NIGP composantes pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk , k variant de 1 à NIGP , correspondant au poids d'un IGPk obtenu avec la méthode de calcul d'interpolation directe définie par le standard pour l'IPPi 20 considéré ; ARes12 2 * Yi = 0Res2 désigne un vecteur-colonne à NiGp composantes AResNIGP2 de la répartition du carré du résidu Resi2 à répartir sur les %FI IGPk lorsque la maille de l'IPPi considéré comporte NiGp IGP, comme résultat de la résolution de l'équation ; 25 .* Req désigne le carré du résidu à répartir, puis .- déterminer (316) un coefficient de pondération pi comme inverse d'un poids, ce coefficient de pondération étant représentatif de la qualité de la mesure du STECi redressé en VTECi par la division avec la fonction de rabattement de l'IPPi considéré, et exprimé suivant la formule suivante : 30 2 aVTECi aVTEC NOMINAL 2 Pi = 3061560 70 dans laquelle 6vrEci2 désigne le bruit de mesure de la pseudo-distance mesuré par le terminal de mesure pour l'IPPi considéré, et 5 crwrEci NOMINAL2 désigne un bruit de mesure nominal de référence du terminal de mesure ; ensuite, après avoir calculé tous les Yi et pi correspondant au premier ensemble des points de percement IPPi observé, pour chaque IGP, calculer (324) une variance de résidu non dilaté 10 suivant l'équation : Nipp 1 1 2 - X X AResipp Ps aMOPS = p' I avec : Ps, la somme des poids normalisés 1/pi de chaque résidu Nipp 1 disponibles pour l'IGP considéré suivant : Ps = Ei=i - Pi 15 Nipp étant le nombre des IPP ayant généré un résidu pour l'IGP concerné, 3'Resipp 2 le résidu interpolé en inverse à l'IGP considéré pour l'IPPi faisant partie des IPP ayant généré un résidu pour l'IGP concerné. 20 9. Procédé de calcul optimal des variances des résidus d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système SBAS selon la revendication 8, dans lequel .- le coefficient de pondération pi est modulé par l'inverse du poids wi de l'IPPi suivant l'IGP tel que défini par le standard MOPS DO-229 et s'écrit 2 2 25 suivant l'expression : pi = - x avTEci et aMOPS = P. Kt)-1. K. P. Xt dans laquelle : 30 Wi avrEc NOMINAL .- l'écart-type amops de l'IGP considéré est calculé suivant l'expression .* K désigne le vecteur unitaire ligne de dimension Nipp : K = [1 1 ... 1], 3061560 71 i [ - ... 0 pi .* P désigne la matrice diagonale des poids P = : '.. : , et 1 PNipp ARes 12 .* X désigne le vecteur colonne des résidus : X = I 2 AResNipp 10. Procédé de calcul optimal des variances des résidus d'un premier 5 ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système SBAS selon la revendication 8, dans lequel .- le coefficient de pondération pi est modulé par l'inverse du poids 0 - - - - wi center norm ayTEC NOMINAL2 10 est calculé selon l'équation : Wi center norm de l'IPPi suivant l'IGP l'expression : 1 X aVTE Ci 2 Pi = dans laquelle le poids wi center norm Si wi = 0, alors l'IPP n'est pas utilisé pour le calcul du o-mops de l'IGP considéré Si 0 11. Procédé de calcul optimal des variances des résidus d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système SBAS selon l'une quelconque des revendications 8 à 10, comprenant une étape supplémentaire consistant à : 5 .- pour chaque IGPk du deuxième ensemble, calculer (326) une variance de résidu dilaté o-Ifilops à partir de la variance de résidu non dilaté amops suivant l'équation : ft ka\ CifilOPS = C(13) - aMOPS G avec dk = (tv(kp'a)) le coefficient de dilatation associé à l'IGPk, k 10 variant de 1 à NIGP, tvk,a désignant le facteur de Student associé l'IGPk de rang k fonction de vk le degré de liberté et fonction de la probabilité de confiance a , ledit facteur de Student tyk,o, permettant de borner une erreur avec un niveau de confiance (1-p) prédéterminé et requis par le service SBAS, la probabilité de 15 confiance a étant lié au niveau de confiance par la relation a = 1 - p , et le degré de liberté vk étant égal ici au nombre d'IPP ayant généré un résidu pour l'IGPk concerné moins un ; et G (p) désignant la valeur obtenue lorsque l'hypothèse gaussienne est appliquée pour la même probabilité de confiance a, c'est à dire la limite du 20 facteur de Student tyk,a lorsque le degré de liberté vk tend vers l'infini. 12. Procédé de calcul optimal simplifié des variances des résidus d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système SBAS, la grille de correction 25 ionosphérique étant structurée en un réseau maillé de cellules de forme polygonale prédéterminée, les cellules du réseau maillé correspondant projectivement à la zone de service SBAS et ayant pour sommets des point IGP du premier ensemble, le procédé de calcul optimal des variances des résidus des IGP 30 comprenant les étapes consistant à : .- calculer par au moins un calculateur électronique, à partir d'informations prédéterminées de correction d'erreur ionosphériques des points IGP du premier ensemble, et de mesures de pseudo-distances d'un deuxième ensemble de points de percement IPP d'observation et de contrôle 3061560 73 contenus dans les cellules de la grille ionosphérique, pour chaque point de percement d'observation IPPi du deuxième ensemble, un résidu, désigné par Resi , ramené à la verticale, suivant l'équation : STECi - UISDi Resi = Fpp pp standard 5 dans lequel : STECi désigne le délai ionosphérique de la ligne de vue réelle mesuré par une station d'observation et fonction de l'élévation de ladite ligne de vue; UISDi désigne le délai ionosphérique vertical interpolé suivant le 10 standard du terminal utilisateur à partir des GIVDk des IGPk de la maille entourant l'IPPi ; -- Fpp standard désigne la fonction de rabattement standard du délai ionosphérique en fonction de l'élévation de la ligne de vue, puis le carré du résidu ; 15 le procédé de calcul optimal des variances des résidus des IGP étant caractérisé en ce qu'il comprend en outre les étapes consistant à, pour chaque point de percement d'observation IPPi du deuxième ensemble, . - déterminer des valeurs d'incréments de variance A -Resk2 à attribuer 20 respectivement aux IGPk de la maille m de l'IPPi considéré par une méthode d'interpolation inverse utilisant une méthode de Moindres Carrés dans laquelle un vecteur Yi de répartition du carré du résidu Resi2 de l'IPPi considéré est calculé suivant l'équation : Yi = Hit. (Hi. Hit)-1. Resi2 25 dans laquelle -*H = w2 --- wNIGP1 désigne un vecteur-ligne à NiGp composantes pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk , k variant de 1 à NIGP , correspondant au poids d'un IGPk obtenu avec la 30 méthode de calcul d'interpolation directe définie par le standard pour l'IPPi considéré ; 3061560 74 Y ARes12 Yi = ARes22 désigne un vecteur-colonne à NiGp composantes - - - AREesNIGP2 de la répartition du carré du résidu Resi2 à répartir sur les NiGp IGPk lorsque la maille de l'IPPi considéré comporte NiGp IGP, comme résultat de la résolution de l'équation ; 5 .* Res? désigne le carré du résidu à répartir, puis .- déterminer un coefficient de pondération pi comme inverse d'un poids, ce coefficient de pondération étant représentatif de la qualité de la mesure du STECi redressé en VTECi par la division avec la fonction de rabattement de l'IPPi considéré, et exprimé suivant la formule suivante : 10 1 avTEci2 pi = X 2 Wi aVTEC NOMINAL dans laquelle avTEci2 désigne le bruit de mesure de la pseudo-distance mesuré par le terminal de mesure pour l'IPPi considéré, 15 wi désigne le poids de l'IPPi suivant l'IGP tel que défini par le standard DO-229 ; et -- 6VTECi2 désigne un bruit de mesure nominal de référence du terminal de mesure ; ensuite, après avoir calculé tous les Yi et pi correspondant au premier 20 ensemble des points de percement IPPi observé, pour chaque IGP, faire un tri des résidus calculés par la méthode d'interpolation inverse en enlevant les résidus calculés à partir des IPPi dont les qualités de mesures sont insuffisantes par rapport à un seuil de qualité prédéterminé, et 25 garder les Nipp jil résidus calculés, puis calculer une variance de résidu simplifié non dilaté suivant l'équation : CYMOPS_simp = 1 Nipp_fil 2 X i=1 AResipp ipp_fil avec : 3061560 75 ^ Nipp fil le nombre des IPP ayant généré un résidu pour l'IGP concerné et dont la qualité de mesure est suffisante ; 3^ `Resipp i 2 est le résidu interpolé en inverse à l'IGP considéré pour l'IPPi faisant partie des IPP triés ayant généré un résidu pour l'IGP concerné. 5 13. Procédé de calcul optimal simplifié des variances des résidus d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système SBAS selon la revendication 10, dans lequel 10 pour chaque IGPk du deuxième ensemble, calculer une variance de résidu simple dilaté creops_simp à partir de la variance de résidu simple non dilaté amops_simp suivant l'équation : c1 tvkji , "MOPS_simp = (p)) - aMOPS_simp tv p)a 15 avec dk = PL-) associé à l'IGPk, k variant de 1 à NiGp, G( tvivot désignant le facteur de Student associé l'IGPk de rang k fonction de vk le degré de liberté et fonction de la probabilité de confiance a , ledit facteur de Student tvk,,, permettant de borner une erreur avec un niveau de confiance (1-p) prédéterminé et requis par le service SBAS, la probabilité de 20 confiance a étant lié au niveau de confiance par la relation cc = 1- p , et le degré de liberté vk étant égal ici au nombre d'IPP ayant généré un résidu pour l'IGPk concerné moins un ; et G (p) désignant la valeur obtenue lorsque l'hypothèse gaussienne est appliquée pour la même probabilité de confiance a, c'est à dire la limite du 25 facteur de Student tvk,a lorsque le degré de liberté vk tend vers l'infini. 14. Procédé de calcul optimal des variances des résidus d'un premier ensemble de points IGP d'une grille de correction ionosphérique d'une zone de service d'un système SBAS selon l'une quelconque des revendications 8 30 à 11, dans lequel La forme de chaque maille est une forme triangulaire et le nombre NiGp des IGP sommets de chaque maille est égal à 3, ou 3061560 76 La forme de chaque maille est la forme d'un quadrilatère, de préférence comprise parmi les formes d'un rectangle, d'un carré et d'un losange, et le nombre NiGp des IGP sommets de chaque maille est égal à 4. 5 15. Système SBAS d'augmentation à base de satellites pour augmenter les performances d'un système global de navigation par satellites GNSS, le système SBAS comprenant des terminaux utilisateurs de service(s) SBAS, un ou plusieurs satellites d'augmentation des satellites du système 10 GNSS et de diffusion de messages d'information aux terminaux utilisateurs, une ou plusieurs stations d'observation RIMS munis de récepteurs GNSS, et un ou plusieurs ou calculateurs, le système SBAS étant configuré pour calculer de manière optimale des variances des résidus d'un premier ensemble de points IGP d'une grille 15 de correction ionosphérique d'une zone de service d'un système SBAS, la grille de correction ionosphérique étant structurée en un réseau maillé de cellules de forme polygonale prédéterminée, les cellules du réseau maillé correspondant projectivement à la zone de service SBAS et ayant pour sommets des point IGP du premier ensemble, 20 le système SBAS étant caractérisé en ce que le ou les calculateurs électroniques sont configurés pour : .- à partir d'informations prédéterminées de correction d'erreur ionosphériques des points IGP du premier ensemble, et de mesures de pseudo-distances d'un deuxième ensemble de points de percement IPP 25 d'observation et de contrôle contenus dans les cellules de la grille ionosphérique, pour chaque point de percement d'observation IPPi du deuxième ensemble, un résidu, désigné par Resi , ramené à la verticale, suivant l'équation : Resi = STECi - UISDi Fpp standard 30 dans lequel : STECi désigne le délai ionosphérique de la ligne de vue réelle mesuré par une station d'observation et fonction de l'élévation de ladite ligne de vue; 3061560 77 UISDi désigne le délai ionosphérique vertical interpolé suivant le standard du terminal utilisateur à partir des GIVDk des IGPk de la maille entourant l'IPPi ; Fpp standard désigne la fonction de rabattement standard du délai 5 ionosphérique en fonction de l'élévation de la ligne de vue, puis calculer le carré du résidu ; et ensuite pour chaque point de percement d'observation IPPi du deuxième ensemble, . - déterminer des valeurs d'incréments de variance A -Resk2 à attribuer 10 respectivement aux IGPk de la maille m de l'IPPi considéré par une méthode d'interpolation inverse utilisant une méthode de Moindres Carrés dans laquelle un vecteur Yi de répartition du carré du résidu Resi2 de l'IPPi considéré est calculé suivant l'équation : Yi = Hit. (Hi. Hit)-1. Resi2 15 dans laquelle -* H = [w1 w2 --- WNIGP1 désigne un vecteur-ligne à NiGp composantes pour réaliser l'interpolation de l'IPPi considéré, chaque composante wk , k variant de 1 à NIGP , correspondant au poids d'un IGPk obtenu avec la 20 méthode de calcul d'interpolation directe définie par le standard pour l'IPPi considéré ; [3`Res12 Yi = 3`Res22 désigne un vecteur-colonne à NiGp composantes AResNIGP2 de la répartition du carré du résidu Resi2 à répartir sur les NiGp IGPk lorsque la maille de l'IPPi considéré comporte NiGp IGP, comme résultat de la 25 résolution de l'équation ; .* Res? désigne le carré du résidu à répartir, puis .- déterminer un coefficient de pondération pi comme inverse d'un poids, ce coefficient de pondération étant représentatif de la qualité de la mesure du STECi redressé en VTECi par la division avec la fonction de 30 rabattement de l'IPPi considéré, et exprimé suivant la formule suivante : 3061560 78 2 aVTECi pl = 2 "VTEC NOMINAL dans laquelle 6vTEci2 désigne le bruit de mesure de la pseudo-distance mesuré par le terminal de mesure pour l'IPPi considéré, et 5 aVTECi 2 désigne un bruit de mesure nominal de référence du terminal de mesure ; ensuite, après avoir calculé tous les Yi et pi correspondant au premier ensemble des points de percement IPPi observé, pour chaque IGP, calculer une variance de résidu non dilaté suivant 10 l'équation : CYMOPS = - X - X ,AResipp Ps p 1=1 i Nipp 1 1 2 avec - Ps, la somme des poids normalisés 1/pi de chaque résidu disponibles pour l'IGP considéré suivant : Ps = EN»PP 1=1 pi 15 Nipp étant le nombre des IPP ayant généré un résidu pour l'IGP concerné, AResipp i 2 le résidu interpolé en inverse à l'IGP considéré pour l'IPPi faisant partie des IPP ayant généré un résidu pour l'IGP concerné. 20