FR 2939945 A1 20100618 FR 0807060 A 20081216 PIROCEDE DE CORRECTION DE L'ALTITUDE BAROMETRIQUE D'AERONEFS MIS EN OEUVRE DANS UN SYSTEME DE CONTROLE DU TRAFIC AERIEN La présente invention concerne un procédé de correction de l'altitude barométrique d'aéronefs mis en oeuvre dans un système de contrôle du trafic aérien. Elle concerne notamment les systèmes de pistage d'aéronefs sur la base de multiples capteurs. 10 Le contrôle du trafic aérien, ou ATC selon l'acronyme de Air Traffic Contrcl, permet à des contrôleurs aériens d'assurer l'exécution sûre, rapide et efficace des vols d'aéronefs dans l'espace aérien sous surveillance. Son rôle est principalement de prévenir les collisions entre les aéronefs et le sol 15 ou les véhicules, ainsi que les collisions en vol entre aéronefs. II consiste aussi à accélérer et ordonner la circulation aérienne, à fournir aux aéronefs les av s et renseignements utiles à l'exécution sûre et efficace du vol, tels que les informations météorologiques, l'information sur l'état des moyens au sol de navigation, l'information sur le trafic. Il consiste enfin à fournir un 20 service d'alerte pour prévenir les organismes appropriés lorsque les aéronefs ont besoin de l'aide des organismes de secours et de sauvetage, et à prêter à ces organismes le concours nécessaire. Les données nécessaires à la gestion du contrôle du trafic aérien proviennent essentiellement d'une pluralité de capteurs. Parmi ces capteurs, 25 des radars primaires de surveillance ou PSR selon l'acronyme de Primary Surveillance Radars, fournissent des échos de cibles, via des informations d'azimut, de distance et de l'instant de détection. Des données sont aussi fournies par des radars secondaires de surveillance ou SSR selon l'acronyme de Secondary Surveillance Radars ; 30 un SSR envoie des signaux à des cibles. Lorsque ces cibles sont des aéronefs équipés de transpondeurs, ces derniers envoient en retour des informations, qui sont reçues par le SSR. Les informations reçues par le SSR comprennent la distance et l'azimut, ainsi que les informations complémentaires fournies par le transpondeur selon son type. Ainsi un 35 transpondeur compatible avec le mode A/C fournit une donnée d'identification de l'aéronef (mode A), et une donnée d'altitude barométrique5 (mode C). L'ensemble des données fournies par un SSR permet donc une identification des aéronefs évoluant dans son champ de vision, ainsi que la connaissance de leurs positions respectives dans un repère tridimensionnel. II existe encore des transpondeurs équipés d'un mode perfectionné ou Mode S, fournissant sur requête les mêmes données que le mode A/C, en plus d'une identification univoque de l'aéronef codée sur 24 bits, et de moyens de communication bilatérale de données diverses. D'autres données de localisation d'aéronefs et objets au sol sont fournies par des capteurs de multilatération ou MLAT, consistant en une ~o pluralité d'antennes omnidirectionnelles disséminées sur le sol, recevant des signaux envoyés par un aéronef afin de le localiser. Ces signaux peuvent être non sollicités, ou bien envoyés par l'aéronef en réponse de signaux provenant de radars. Des calculs portant sur les différences entre les temps de réception de ces signaux entre au moins deux antennes permettent de 15 déterminer la position de l'aéronef. Ce type de capteur est largement utilisé pour la surveillance de mouvements au sol dans des zones d'aéroport et dans les zones d'approche. Une autre technique procédant d'une manière similaire est utilisée pour la localisation d'aéronefs en dehors de zones d'approche d'aéroports, basée sur les signaux reçus par une pluralité 20 d'ante-mes omnidirectionnelles disposées sur des zones géographiques beaucoup plus larges ; cette technique porte le nom de WAM, selon l'acroryme pour Wide Area Multilateration. II existe aussi des systèmes de surveillance automatique et dépendante des aéronefs, connus sous les acronymes d'ADS pour 25 Automatic Dependent Surveillance. Un premier système d'ADS, dénommé ADS-C selon l'acronyme pour Automatic Dependent Surveillance û Contract, est utilisé essentiellement pour les zones dans lesquelles peu de capteurs d'autres types sont disponibles, par exemple les zones désertiques ou océaniques. Les aéronefs qui en sont équipés radiodiffusent vers un satellite 30 de communication des données relatives à leur position, par exemple déterminée par un calculateur embarqué sur la base de données fournies par un récepteur de géo-localisation par satellite ou GPS selon l'acronyme pour Global Positioning System, et / ou par une centrale inertielle. D'autres données sont diffusées, telles que des données sur la route envisagée par 35 l'aéronef, les vitesses de l'aéronef par rapport à l'air et au sol, des données météorologiques (force et direction du vent, température, etc.). La diffusion de ces données par l'aéronef peut se faire de manière périodique, ou bien en réaction à certains évènements, ou encore en cas de situation d'urgence. Les données sont ensuite transmises par le satellite et sont reçues par une antenne dédiée. Un second système d'ADS est dénommé ADS-B selon l'acronyme pour Automatic Dependent Surveillance û Broadcast. Les aéronefs qui en sont équipés radiodiffusent les données décrites précédemment en référence à l'ADN-C, directement vers une antenne au sol. La transmission se fait de manière périodique, à une fréquence bien plus élevée, d'au moins une émission par seconde. Toutes les données en provenance de tous les capteurs précités, après avoir parfois subi un prétraitement au niveau des capteurs mêmes, sont recueillies par un central de contrôle de trafic aérien, où elles sont traitées, avant d'être diffusées sous forme synthétique aux contrôleurs aériens via des interfaces d'affichage. Il est primordial pour les contrôleurs de pouvoir procéder au pistage des aéronefs présents dans l'espace aérien. Le pistage d'un aéronef consiste à suivre sa piste, c'est-à-dire la représentation de l'évolution de l'aéronef dans l'espace aérien, en projection sur le plan stéréographique situé au niveau de la mer. Une autre donnée d'importance majeure dans la localisation d'un aéronef, donc son pistage, est l'altitude. Les données fournies par un PSR par exemple, qui sont prises en compte dans la détermination de la position d'un aéronef, ne contiennent aucune composante liée à l'altitude : une détection d'un PSR peut correspondre à une multitude d'altitudes différentes ; réciproquement sans la donnée d'altitude, la cible détectée correspond à une multitude de projections possibles sur la surface terrestre. Les données fournies par un SSR contiennent une composante liée à l'altitude : il s'agit de l'altitude barométrique fournie par le mode C du transpondeur. Néanmoins l'altitude barométrique transmise par le mode C et issue d'un alticodeur, est donnée en référence à la pression standard au niveau de la mer dans des conditions normales de pression et de température, à savoir : 1013,25 hectoPascal pour une température de 15° C. La pression atmosphérique pouvant varier significativement d'une zone géographique à une autre zone géographique relativement proche, la donnée de l'altitude barométrique mode C utilisée dans les calculs doit être assortie d'une forte variance. L'incertitude ou la variance quant à l'estimation de l'altitude d'un aéronef se répercute sur les calculs de position, ainsi que sur tous les calculs mis en oeuvre dans un système de contrôle du trafic aérien, comme par exemple les calculs réalisés dans le but de déclencher des alertes en cas de risques de collisions entre aéronefs. Il est donc souhaitable d'appliquer une correction à l'altitude barométrique des radars mode A/C, afin d'améliorer la précision de la position, et de réduire la variance. ~o En outre, il est utile pour les contrôleurs de connaître les altitudes géométriques des aéronefs, ou bien des altitudes barométriques corrigées en fonction des conditions météorologiques dans les zones dans lesquelles les aéronefs se situent. Une altitude barométrique corrigée est fournie en référence au QNH : le QNH prend en considération la valeur de pression 15 actuelle au niveau de la mer. Le central de contrôle du trafic aérien n'effectue une correction des valeurs d'altitude barométrique en fonction du QNH que sur la zone dans laquelle il est situé, c'est-à-dire typiquement sur la zone d'aéroport, en fonctior de données fournies par des services météorologiques. Ces 20 données sont mises à jour avec une fréquence très faible (par exemple : une fois par jour) : non seulement la pression dans la zone d'aéroport peut être amenée à varier de manière significative dans le temps, entre deux mises à jour de la donnée de pression atmosphérique par le service météorologique, mais encore elle peut varier significativement entre la zone d'aéroport et une 25 zone distante de l'aéroport, dans la zone d'intérêt. Pour cette raison, la correction d'altitude barométrique n'est effectivement réalisée que dans les zones d'aéroport, pour des aéronefs évoluant à de faibles altitudes, c'est-à-dire à des altitudes inférieures à l'altitude de transition. Dans une zone d'aéroport, l'altitude barométrique est prise en compte par le filtre en position, 30 et considérée comme une altitude géométrique avec une variance associée augmentée d'une valeur forfaitaire visant à compenser le maximum de variatioi de la correction barométrique. Ainsi, en l'absence de correction de l'altitude barométrique, c'est-à- dire pour les aéronefs situés dans des zones extérieures aux zones 35 d'aéroport, la donnée d'altitude est biaisée, plus ou moins fortement, en fonction des conditions atmosphériques. Par conséquent, l'évaluation de la précision du pistage est rendue aléatoire suivant que les conditions atmosphériques soient favorables ou défavorables, c'est-à-dire plus ou moins proches des conditions atmosphériques standard. En outre, dans des conditions atmosphériques favorables, c'est-à-dire dans des zones où la pressio-i atmosphérique est proche de la pression atmosphérique standard, la composante en altitude est surlissée, ce qui dégrade la précision du filtre en posi.:ion, en particulier lorsque la cible est manoeuvrante. Un but de la présente invention est de pallier les inconvénients précités, en proposant notamment un procédé apte à appliquer aux altitudes barométriques des aéronefs présents dans une zone de surveillance, une correction appropriée en fonction des conditions atmosphériques environnant chacun des aéronefs, afin de réduire la variance sur les valeurs d'altitude et ainsi la précision des calculs, sans que des données spécifiques ne soient fournies par des services météorologiques. A cet effet, l'invention a pour objet un procédé de correction de l'altitude barométrique d'aéronefs, caractérisé en ce qu'il comprend les 20 étapes suivantes : • une étape préalable de découpage du plan stéréographique correspondant à une zone de surveillance en une grille bidimensionnelle comprenant une pluralité de cellules ; • une première étape de démarrage d'une fenêtre temporelle 25 d'une durée déterminée, de détection de rapports de type ADS-B, pour chaque rapport de type ADS-B reçu pendant la fenêtre temporelle de détection : • une étape d'acquisition et de mémorisation des couples de 30 valeurs d'altitude barométrique et altitude géométrique (C ; H) issus des rapports ADS-B associés aux cellules correspondantes ; après la fin de la fenêtre temporelle de détection : • une étape de démarrage d'une nouvelle fenêtre temporelle de 35 la durée déterminée ; • une étape de calcul de valeurs de pression et de température (Po ; To) au niveau de la mer pour toutes les cellules, à partir de tous les couples (C ; H) acquis pour ces cellules au cours d'une période déterminée, et de calcul d'une correction de l'altitude barométrique à partir des valeurs de température et de pression au niveau de la mer calculées, pour toutes les cibles détectées par au moins un capteur, à la verticale des cellules. Dans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus, peut être caractérisé en ce cue le calcul des valeurs (Po ; To) de pression et de température au niveau de la mer réalisé lors de l'étape de calcul au terme de la fenêtre temporelle est opéré sur la base de deux couples de valeurs (C; ; H;) et (Ci ; Hi) acquis au cours de la fenêtre temporelle, selon les formules suivantes : Postandard na(Cj. H1 - C.H. ) 1- TOstandard (Cl ù Ci) et L 0 = Tostandard ù P \ Ci 1ù PO Ostandard ttPo / i Dans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus peut être caractérisé en ce que ledit au moins un capteur est un radar secondaire mode A/C, la correction de l'altitude barométrique réalisée lors de l'étape de calcul pour une détection en provenance d'un radar secondaire de surveillance fournissant une altitude barométrique corrigée selon la formule suivante : T P ' ùT ' 0 T C( Ostandard )n + Ostandard 1 ù ( Ostandard Tostandard PO a To P0 ù i i \\ H._ Tostandard a PO ù Postandard nPa Ccorr Dans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus peut être caractérisé en ce que les couples de valeurs (C ; H) mémorisés lors de l'étape d'acquisition et de mémorisation sont supprimés s'ils correspondent à des détections à des niveaux d'altitude inférieurs à un premier seuil déterminé et / ou supérieurs à un second seuil déterminé. Dans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus peut être caractérisé en ce que les deux couples de valeurs (C; ; H;) et (C; ; Hj) sont respectivement les deux couples correspondant aux deux détections de plus basse altitude géométrique H. Dans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus peut être caractérisé en ce que l'étape de calcul comprend en outre une étape de vérification de la cohérence des valeurs des couples (C ; H) mémorisés pour une cellule donnée, un couple de valeurs (C ; H) étant considéré comme cohérent si la différence entre la valeur corrigée de l'altitude barométrique C et l'altitude géométrique H est inférieure à un seuil déterminé ou incohérent dans le cas contraire ; un ou plusieurs couples de valeurs considérés comme incohérents entraînant le choix de deux couples (C; ; H;) et (Cj ; Hi) différents. Dans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus peut être caractérisé en ce que des valeurs de pression et de température au niveau de la mer (Po ; To) sont calculées lors de l'étape de calcul, pour toutes les paires de couples de valeurs (C ; H) d'une cellule donnée, l'étape de calcul comprenant un calcul supplémentaire fournissant une moyenne des valeurs (Po ; To) calculées et appliquant à la cellule la valeur moyenne du couple (Po ; To). Dans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus peut être caractérisé en ce que l'étape de calcul comprend un calcul supplémentaire des valeurs de pression et de température au niveau de la mer (Po ; To) pour des cellules pour lesquelles aucun couple de valeurs (Po ; To) n'est attribué et juxtaposées à au moins un nombre déterminé de cellules adjacentes pour lesquelles des valeurs (Po ; To) ont été déterminées, le calcul attribuant aux cellules la moyenne barycentrique des valeurs (Po ; To) de toutes les cellules adjacentes. Dans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus peut être caractérisé en ce que le nombre déterminé de cellules adjacentes pour lesquelles des io valeurs (Po ; To) ont été déterminées est égal à 4. Dans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus peut être caractérisé en ce qu'aucune valeur du couple (Po ; To) n'est calculé pour une cellule donnée 15 lors de l'étape de calcul si le nombre de couples de valeurs (C ; H) mémorisés pour cette cellule est inférieur à un nombre prédéterminé. Clans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus peut être caractérisé en 20 ce que le nombre prédéterminé de couples de valeurs (C ; H) mémorisés est égal à 5. Cans un mode de réalisation de l'invention, le procédé de correction de l'altitude barométrique d'aéronefs décrit ci-dessus peut être caractérisé en 25 ce que les valeurs de pression et de température au niveau de la mer (Pon ; To,n) pour une cellule donnée, telles que calculées lors de l'étape de calcul pour urne fenêtre temporelle n, subissent un lissage avec les valeurs (Pon_, ; Toä_,) déterminées pour la cellule pour une fenêtre temporelle précédente n-1, le lissage attribuant au couple de valeurs (Po ; To) la valeur a.(Pon-, ; Ton-,) 30 + (1 - a).(Pon ; Ton), a étant un coefficient compris entre 0 et 1. L invention peut enfin avoir pour objet un système de contrôle du trafic aérien caractérisé en ce qu'il met en oeuvre un procédé de correction de l'altitude d'aéronefs tel que décrit ci-dessus selon divers modes de 35 réalisations de l'invention. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description, donnée à titre d'exemple, faite en regard des dessins annexés qui représentent : la figure 1, un schéma synoptique de la structure générale d'un système de contrôle du trafic aérien et des différents capteurs impliqués ; la figure 2, un exemple de courbe représentant la correction apportée à une altitude barométrique pour estimer une altitude géométrique, et la répercussion de la correction sur la position d'un aéronef donnée par un radar mode A/C ; la figure 3, un exemple du découpage en cellules du plan stéréographique correspondant à la zone de surveillance dans le cadre du contrôle du trafic aérien ; la figure 4, un ordinogramme représentant le procédé 20 d'estimation d'altitude d'aéronefs selon l'invention ; la figure 5, un diagramme représentant un nuage de couples de valeurs d'altitudes barométriques et géométriques reçues au cours d'une fenêtre temporelle. 25 La figure 1 présente un schéma synoptique de la structure générale d'un système de contrôle du trafic aérien. Un aéronef 100 est en vol, il est la cible d'un radar primaire de surveillance PSR 101 monté avec un radar secondaire de surveillance SSR 102 dans cet exemple. Une antenne 30 omnidirectionnelle 103 est dédiée à l'ADS-B. Un centre de contrôle du trafic aérien 110 recueille toutes les données en provenance des différents capteur:;. Le PSR 101 émet une onde vers l'aéronef, et capte un écho de l'aéronef 100 lorsque celui-ci se trouve dans son champ de vision. L'analyse de l'onde réfléchie permet une estimation de la distance entre l'aéronef 100 35 et le PSR 101, ainsi que l'azimut. Le SSR émet une onde, lorsque celle-ci 10 15 est cap':ée par l'aéronef 100, le transpondeur de ce dernier, dans l'exemple de la figure compatible avec le mode NC, émet à son tour un signal comprenant la distance et l'azimut, ainsi que la donnée d'identification de l'aéronef, et la donnée d'altitude barométrique. Dans l'exemple de la figure, l'aéronef 100 est équipé d'un système ADS-B. Il émet donc périodiquement un signal, ou rapport ADS-B, contenant des informations complémentaires fournies par exemple par un récepteur embargaé de géo-positionnement par satellite. Parmi ces données complémentaires, se trouvent les coordonnées géodésiques de l'aéronef 100, selon le système connu sous l'acronyme de WGS 84, ainsi que son altitude barométrique. Ainsi chaque rapport ADS-B comprend, entre autres données, l'altitude géométrique H de l'aéronef 100 et son altitude barométrique C. Le taux d'équipement ADS-B des avions de compagnies aériennes est aujourd'hui tel qu'il permet de couvrir les principales zones d'intérêt du trafic aérien. Un but de la présente invention est de mémoriser des couples de données d'altitudes barométriques et d'altitudes géométriques par zones géographiques et par altitudes géométriques ou niveaux de vol, afin de permettre le calcul, à partir d'au moins deux mesures récentes de tels couples, de la température et de la pression au niveau de la mer. Ce calcul permet de réaliser une cartographie météorologique de la zone d'intérêt. Par exemple, un découpage en cellules du plan stéréographique ramené au niveau moyen de la mer peut être réalisé. Par une pluralité de calculs, opérés sur la base de couples (C ; H) obtenus au cours du temps de divers aéronefs, il est possible d'associer à chacune des cellules des valeurs de pression et de température au niveau du sol. Un couple de pression et de température au niveau du sol, que nous noterons par la sjite (Po ; To), peut être déterminé pour une cellule donnée, dès lors qu'au moins deux couples (C ; H) ont été obtenus d'aéronefs à la verticale de cette cellule. Le calcul des valeurs (Po ; To) peut être réalisé à partir de la formule suivante, selon le modèle d'atmosphère de l'Organisation de l'Aviation Civile Internationale (OACI) combinant la relation d'état des gaz parfaits et l'hypothèse hydrostatique : H = To C( Postandard ) n + TOstandard Tetandard Po a • H est l'altitude géométrique, • C est l'altitude barométrique, • -fo est la température au niveau de la mer, avec Tostandard = 288,15 K, • Po est la pression au niveau de la mer, avec Postandard = 1013,25 hPa, • a est la chute de la température avec l'altitude ; ici a = 0,0065 °C.m-1, • n = 5,2561 Ainsi pour une paire de couples (C ; H) notés (C; ; H;) et (Ci ; Hi) pour des relevés à la verticale d'une même cellule donnée, la pression et la température au niveau de la mer peuvent s'écrire : Postandard Pä = ù ù (TOstandard ) ä To ii (1),où: na(Cj.H. -C.Hj.) (2), 1 TOstandard (C; To = TOstandard ~P -P \\ Ostandard 0 Ostandard a nPo Pù P \ C, 1 ù 0 Ostandard nPo A partir du moment où des valeurs de pression et de température au niveau de la rrer Po et To ont été déterminées pour une cellule donnée, c'est-à-dire initialisoes et validées comme explicité ci-après en référence aux figures 4 et 5, une correction de l'altitude barométrique pourra être effectuée par exemple pour toutes les détections radar mode NC identifiées à la verticale de cette cellule, comme explicité ci-après en référence à la figure 3. II est bien ertendu que d'autres types de détections identifiées à la verticale de cette cellule, par exemple en provenance d'autres capteurs tels que le MLAT ou le VVAM, pourront de la même manière subir une correction de l'altitude barométrique. La figure 2 présente un exemple de courbe représentant la correct on apportée à une altitude barométrique pour estimer une altitude géométrique, et la répercussion de la correction sur la position d'un aéronef 100 non représenté sur cette figure, donnée par un radar mode A/C. Dans cet exemple, une détection réelle d'un radar secondaire de surveillance SSR non représenté sur la figure, mais dont l'antenne est située à l'origine O du repère tridimensionnel orthonormé (O ; x ; y ; z), est représentée par un point 201 ; les coordonnées du point 201 dans le repère orthonormé sont x, y et C. Le point 201 est à une distance p du SSR, et est relevé à un azimut Az, représenté sur la figure en faisant coïncider l'axe des abscisses du repère orthonormé avec l'azimut de référence du SSR. En outre le point 201 est relevé à une altitude barométrique C. Ces données sont fournies dans le rapport mode A/C, ainsi cette altitude C est donnée en référence à la pression standard de 1013,25 hPa. En admettant que pour la cellule (non représentée sur la figure) à la verticale de laquelle se situe le point 201, des valeurs Po et To ont été déterminées, alors l'altitude barométrique C peut subir une correction. L'altitude barométrique après correction est notée Ccorr. est calculé selon la relation (1) décrite plus haut : C = To C.! Postandard )ä +Tostandard 1 _ (Tostandard )ä corr, \\ l Tostandard Po a To II est alors possible de reconstruire une détection SSR virtuelle, représentée par un point virtuel 202, situé de même à la distance p du SSR, à l'azimut Az, et dont l'altitude est l'altitude barométrique corrigée Ccorr. L'abscisse et l'ordonnée du point virtuel 202 sont notées respectivement xcorr et ycorr. Les trois coordonnées du point 202 diffèrent de celles du point 201, et leur précision est améliorée puisque l'altitude barométrique Ccorr correspond au QNH local. Ainsi tous les calculs réalisés sur la base de la positior de l'aéronef 100 sont réalisés à partir de la position corrigée 30 correspondant au point virtuel 202, ce qui évite de propager dans ces calculs l'erreur introduite par l'estimation de l'altitude barométrique en référence à la pressiol de 1013,25 hPa. Ccorr 25 La figure 3 présente un exemple du découpage en cellules du plan stéréographique correspondant à la zone de surveillance dans le cadre du contrôle du trafic aérien selon l'invention. Une grille 301 peut par exemple compr{ ndre des cellules de forme carrée, à l'instar d'une cellule 302 représentée. Les cellules adjacentes à la cellule 302 sont représentées par les références 303 à 310. Dans un exemple de réalisation de l'invention, chaque cellule peut avoir un côté de 10 Nm. Dans un premier exemple, une détection de SSR peut se situer à la verticale de la cellule 302, pour la piste d'un aéronef 100 n'étant pas équipé d'un système ADS-B. L'altitude barométrique donnée par le rapport du mode C peut alors être corrigée en se basant sur les valeurs de pression et de température Po et To déterminées pour cette cellule. Si les valeurs de Po et To n'ont pas pu être déterminées pour cette cellule 302, alors aucune correction n'est appliquée. Avantageusement, si les valeurs de Po et To n'ont pas pi. être déterminées pour la cellule 302, mais que des valeurs de pression et de température au niveau de la mer ont pu être déterminées pour des cellules adjacentes 303 à 309, alors des valeurs Po et To peuvent être attribuées par interpolation à la cellule 302. Par exemple, si des valeurs ont été déterminées pour au moins 4 des cellules 303 à 309, alors des valeurs Po et To peuvent être calculées comme le barycentre de toutes les valeurs attribuées aux cellules adjacentes concernées. La figure 4 présente un ordinogramme représentant le procédé d'estimation d'altitude d'aéronefs selon l'invention. A un instant to, une fenêtre temporelle de durée déterminée démarre, par la réinitialisation d'un compteur de temps. Cette durée peut par exemple durer 20 minutes. Au cours de cette fenêtre temporelle, des rapports ADS-B peuvent être reçus. A la réception d'un rapport ADS-B en provenance d'un aéronef, une première étape 401 a pour but de vérifier si des données d'altitude géométrique et d'altitude barométrique sont comprises dans le rapport. Si tel n'est pas le cas, alors aucune exploitation n'est faite des données reçues. Si de telles données sont présentes dans le rapport, alors une étape 402 suivante a pour but de déterminer la cellule de la grille bidimensionnelle 301 à la verticale de laquelle se trouve l'aéronef dont est issu le rapport ADS-B, au moment de son émission. Lors d'une étape 403 suivante, le couple (C ; H) nouvellement acquis est mémorisé pour la cellule considérée, la mémoire associée à cette cellule est ainsi augmentée d'un couple de valeurs supplémentaires, en plus de couples de valeurs éventuellement acquis précédemment au cours de la fenêtre temporelle courante ou bien au cours de fenêtres temporelles précédentes. Avantageusement, lors d'une étape, non représentée dans la figure, préalable à l'étape de mémorisation 403, il peut être vérifié que le io couple de valeurs nouvellement acquis ne présente pas une aberration, pour cause par exemple d'un dysfonctionnement quelconque. Par exemple, il peut être vérifié que la différence entre les données comprises dans le dernier rapport ADS-B et les données comprises dans un ou des rapports ADS-B précédents pour le même aéronef, n'est pas plus importante qu'un seuil 15 déterminé. Un couple présentant une telle aberration peut alors être définitivement discriminé afin de ne pas générer d'erreurs dans les calculs ultérieurs. Cette discrimination peut tout aussi bien être réalisée lors d'une étape ultérieure. Un certain nombre de couples (C ; H) provenant de différents 20 rapports ADS-B pour un ou une pluralité d'aéronefs sont ainsi mémorisés au cours de la fenêtre temporelle. Lorsque la fenêtre temporelle courante est terminée, une nouvelle fenêtre temporelle est démarrée. En parallèle, les données acquises et mémorisées pendant la période qui vient de se terminer, sont traitées lors d'une étape 404. L'étape 25 404 a pour but de calculer les valeurs Po et To pour toutes les cellules pour lesquelles ce calcul est réalisable, sur la base des relations (2) et (3) explicitées plus haut, et selon une méthode décrite ci-après sur la base d'un exemple, en référence à la figure 5. Avantageusement, les calculs réalisés à cette étape pour une cellule donnée peuvent n'être effectués qu'à la 30 condition qu'au moins un nombre déterminé de couples (C ; H) soit disponible pour cette même cellule. En effet, d'une part il est nécessaire d'avoir au moins deux couples (C ; H) pour pouvoir déterminer les conditions de pression et de température au niveau de la mer ; d'autre part, il est préférable d'avoir un plus grand nombre de couples à disposition, afin 35 d'évite ' de fonder des calculs sur des valeurs isolées, par exemple en provenance d'un unique aéronef, pour se prévaloir plus efficacement contre des aberrations ayant pour cause un éventuel dysfonctionnement, par exemple des satellites de géo-positionnement. Avantageusement, le nombre minimum de couples nécessaires à la réalisation des calculs peut par exemple être fixé à 5. Une fois que les calculs des conditions de pression et de température au niveau de la mer ont été réalisés pour toutes les cellules possibles, les corrections d'altitude barométrique pourront être appliquées pour tous les aéronefs pour lesquels par exemple des détections de radar de io surveillance secondaire SSR 102 auront été reçues à la verticale de ces cellules:. La figure 5 présente un diagramme représentant un nuage de couples de valeurs d'altitudes barométriques et géométriques reçues au 15 cours d'une fenêtre temporelle. Dans un souci de clarté, la figure se limite à une cellule 302 donnée. Bien sûr, la description ci-après s'applique pour toutes les cellules de la grille bidimensionnelle correspondant à l'espace aérien d'intérêt. En outre, la description ci-après est faite en référence à une fenêtre temporelle 20 d'acqui:3ition des rapports ADS-B, démarrant à l'instant t et se terminant à l'instant t + At. Soient (Pon_, ; Ton_,) les valeurs de pression et de température au niveau de la mer associées à la cellule 302 à l'instant t, issues d'un calcul réalisé à la fin de la fenêtre temporelle de détection précédente ; soient (Pon ; Ton) l'es valeurs de pression et de température au niveau de la mer 25 associées à la cellule 302 à l'instant t + At, c'est-à-dire à l'issue de la fenêtre temporelle de détection courante. Au cours de la fenêtre temporelle de détection, 6 couples de valeurs (C ; H) sont acquis dans l'exemple de la figure. Ces valeurs sont notées (C, ; H,) à (C6 ; H6) dans l'ordre chronologique de leur acquisition. 30 Dans un mode de réalisation, seuls les couples (C ; H) acquis au cours de la période temporelle de détection sont utilisés pour le calcul de (Pon ; Ton)• Dans un mode alternatif de réalisation, des valeurs acquises lors d'une ou plusieurs fenêtres temporelles précédentes peuvent être utilisées. Avantageusement, chaque couple de valeurs (C ; H) peut être associé à la 35 date de réception du rapport ADS-B correspondant, et un compteur de temps permet d'éliminer les couples de valeurs ayant atteint un âge déterminé, afin qu'elles ne soient plus utilisées dans les calculs. Un premier couple de valeurs (C ; H) est utilisé comme référence. Par exemple, le couple du niveau d'altitude géométrique le plus bas est utilisé comme référence. Dans l'exemple de la figure, il s'agit du couple de valeurs (C2 ; H2). Avantageusement, les couples de valeurs (C ; H) correspondant à des altitudes géométriques inférieures à un premier seuil d'altitude déterminé, peuvent être éliminés afin de ne pas être utilisés dans les calculs, les altitudes les plus basses étant sujettes à des phénomènes tels que les conditions de l'atmosphère modèle OACI ne peuvent plus y être appliquées en approximation. De la même manière, des couples de valeurs (C ; H) correspondant à des altitudes géométriques supérieures à un second seuil d'altitude déterminé peuvent être éliminés. Un second couple de valeurs (C ; H) est utilisé en combinaison avec le premier couple de valeurs, soit (C2 ; H2) dans l'exemple de la figure, afin de calculer le couple de valeurs (Pon ; Ton). Le second couple de valeurs (C ; H) peut par exemple être choisi comme le second couple de niveau d'altitude géométrique le plus bas, soit le couple de valeurs (C3 ; H3) dans l'exemple de la figure. Avantageusement, il peut être vérifié que les autres couples de valeurs (C ; H) associés à la cellule 302 présentent des valeurs cohérentes. Cette vérification peut par exemple être réalisée en appliquant à chaque valeur d'altitude barométrique C de chaque couple la correction imposée par le couple de valeurs (Pon ; Ton), et en vérifiant que chaque valeur Ccorr résultant de la correction est suffisamment proche de la valeur H correspondante. Si tel n'est pas le cas, il peut alors par exemple être prévu de choisir deux autres couples de valeurs, par exemple (C5 ; H5) et (C, ; H,), pour procéder à un nouveau calcul de (Pon ; Ton) et à une nouvelle vérification avec les autres couples de valeurs (C ; H), et de procéder ainsi jusqu'à ce qu'une telle vérification s'avère satisfaisante. Dans un mode alternatif de réalisation, il peut être envisagé de procéder à différents calculs de (Pon ; Ton) pour toutes les paires de couples possibles par combinaison des couples de valeurs (C ; H) associés à la cellule, puis de procéder à une moyenne des valeurs de pression et de température ainsi déterminées pour obtenir le couple (Pon ; Ton), ou encore il est possible de ne conserver les valeurs (Pon ; Ton) déterminées pour toutes les paires de couples de valeurs (C ; H) possibles qu'à la condition qu'elles soient suffisamment proches des valeurs (Pon-1 ; Ton-,) et / ou qu'elles soient suffisamment proches entre elles. Ces vérifications permettent notamment d'identifier et de s'affranchir d'éventuels effets dépressionnaires pouvant se produire à la verticale de la cellule 302, ou bien d'erreurs plus rares d'origine avionique pouvant conduire à certains rapports ADS-B contenant des valeurs erronées. Dans toutes les situations précitées pouvant conduire à l'impossibilité d'effectuer un calcul du couple de valeurs (Pon ; Ton), il peut être envisagé de n'appliquer aucune correction aux altitudes barométriques des détections faites à la verticale de cette cellule, ou bien d'appliquer les valeurs. (Pon-1 ; Ton-1), assorties d'une augmentation de la variance sur l'amplitude de la correction. Avantageusement les valeurs (Pon ; Ton) peuvent être lissées par 15 rapport aux valeurs (Pon-1 ; Ton-1), par exemple en attribuant au couple de valeurs (Po ; To) un couple de valeurs défini par : a.(Pon-1 ; Ton-1) + (1 - a).(Pon ; Ton) ; a étant compris entre 0 et 1, et typiquement entre 0,3 et 0,4. REVENDICATIONS1. Procédé de correction de l'altitude barométrique d'aéronefs (100), caractérisé en ce qu'il comprend les étapes suivantes : une étape préalable de découpage du plan stéréographique correspondant à une zone de surveillance en une grille bidimensionnelle (301) comprenant une pluralité de cellules (302, 304) ; • une première étape de démarrage d'une fenêtre temporelle d'une io durée déterminée, de détection de rapports de type ADS-B ; pour chaque rapport de type ADS-B reçu pendant la fenêtre temporelle de détection : • une étape (403) d'acquisition et de mémorisation des couples de 15 valeurs d'altitude barométrique et altitude géométrique (C ; H) issus des rapports ADS-B associés aux cellules correspondantes ; après la fin de la fenêtre temporelle de détection : • une étape de démarrage d'une nouvelle fenêtre temporelle de la 20 durée déterminée ; • une étape de calcul (404) de valeurs de pression et de température (Po ; To) au niveau de la mer pour toutes les cellules, à partir de tous les couples (C ; H) acquis pour ces cellules au cours d'une période déterminée, et de calcul d'une correction de l'altitude 25 barométrique à partir des valeurs de température et de pression au niveau de la mer calculées, pour toutes les cibles détectées par au moins un capteur, à la verticale des cellules (302). 2. Procédé de correction de l'altitude barométrique d'aéronefs selon la 30 revendication 1, caractérisé en ce que le calcul des valeurs (Po ; To) de pression et de température au niveau de la mer réalisé lors de 'étape de calcul (404) au terme de la fenêtre temporelle est opéré sur a base de deux couples de valeurs (Ci ; H;) et (Ci ; Hi) acquis au cours de la fenêtre temporelle, selon les formules suivantes : na(C~H ùC Hi) 1ù POstandard , et 5TOstandard (C j ù Ci ) / \ H, ù TOstandard PO ù 1 Ostandard a nPo ( P ùP \ 1 ù 0 Ostandard nP0 3. Procédé de correction de l'altitude barométrique d'aéronefs selon la revendication 1, caractérisé en ce que ledit au moins un capteur est un radar secondaire mode A/C, la correction de l'altitude barométrique réalisée lors de l'étape de calcul (404) pour une détection en provenance d'un radar secondaire de surveillance fournissant une altitude barométrique corrigée selon la formule suivante : 19 T Ci 1 ù (TOstandard )ä \ l TO )/ i TO c( POstandard )ä + TOstandard TOstandard PO a 10 15 20 25 4. Procédé de correction de l'altitude barométrique d'aéronefs selon la revendication 1, caractérisé en ce que les couples de valeurs (C ; H) mémorisés lors de l'étape d'acquisition et de mémorisation (403) sont supprimés s'ils correspondent à des détections à des niveaux d'altitude inférieurs à un premier seuil déterminé et / ou supérieurs à un second seuil déterminé. 5. Procédé de correction de l'altitude barométrique d'aéronefs selon l'une des revendications 2 à 4, caractérisé en ce que les deux couples de valeurs (Ci ; Hi) et (Cj ; Hj) sont respectivement les deux couples correspondant aux deux détections de plus basse altitude géométrique H. 6. Procédé de correction de l'altitude barométrique d'aéronefs selon la revendication 5, caractérisé en ce que l'étape de calcul (404) comprend en outre une étape de vérification de la cohérence desvaleurs des couples (C ; H) mémorisés pour une cellule (302) donnée, un couple de valeurs (C ; H) étant considéré comme cohérent si la différence entre la valeur corrigée de l'altitude barométrique C et l'altitude géométrique H est inférieure à un seuil déterminé ou incohérent dans le cas contraire ; un ou plusieurs couples de valeurs considérés comme incohérents entraînant le choix de deux couples (Ci ; H;) et (Ci ; Hi) différents. 7. Procédé de correction de l'altitude barométrique d'aéronefs selon la revendication 1, caractérisé en ce que des valeurs de pression et de température au niveau de la mer (Po ; To) sont calculées lors de l'étape de calcul (404), pour toutes les paires de couples de valeurs (C ; H) d'une cellule donnée (302), l'étape de calcul (404) comprenant un calcul supplémentaire fournissant une moyenne des valeurs (Po ; To) calculées et appliquant à la cellule (302) la valeur moyenne du couple (Po ; To). 8. Procédé de correction de l'altitude barométrique d'aéronefs selon la revendication 1, caractérisé en ce que l'étape de calcul (404) comprend un calcul supplémentaire des valeurs de pression et de température au niveau de la mer (Po ; To) pour des cellules pour lesquelles aucun couple de valeurs (Po ; To) n'est attribué et juxtaposées à au moins un nombre déterminé de cellules adjacentes pour lesquelles des valeurs (Po ; To) ont été déterminées, le calcul attribuant aux cellules la moyenne barycentrique des valeurs (Po ; To) de toutes les cellules adjacentes. 9. Procédé de correction de l'altitude barométrique d'aéronefs selon la -evendication 8, caractérisé en ce que le nombre déterminé de cellules adjacentes pour lesquelles des valeurs (Po ; To) ont été déterminées est égal à 4. 10. Procédé de correction de l'altitude barométrique d'aéronefs selon l'une quelconque des revendications précédentes, caractérisé en ce qu'aucune valeur du couple (PO ; TO) n'est calculé pour une celluledonnée (302) lors de l'étape de calcul (404) si le nombre de couples de valeurs (C ; H) mémorisés pour cette cellule (302) est inférieur à un nombre prédéterminé. 11. Procédé de correction de l'altitude barométrique d'aéronefs selon la revendication 10, caractérisé en ce que le nombre prédéterminé de couples de valeurs (C ; H) mémorisés est égal à 5. 12. Procédé de correction de l'altitude barométrique d'aéronefs selon l'une quelconque des revendications précédentes, caractérisé en ce que les valeurs de pression et de température au niveau de la mer (Pon ; Ton) pour une cellule (302) donnée, telles que calculées lors de l'étape de calcul (404) pour une fenêtre temporelle n, subissent un lissage avec les valeurs (Pon_1 ; T0n_1) déterminées pour la cellule (302) pour une fenêtre temporelle précédente n-1, le lissage attribuant au couple de valeurs (Po ; To) la valeur a.(Pon-1 ; Ton-1) + (1 - a).(Pon ; Ton), a étant un coefficient compris entre 0 et 1. 13. Système de contrôle du trafic aérien caractérisé en ce qu'il met en oeuvre un procédé de correction de l'altitude d'aéronefs selon l'une quelconque des revendications précédentes.