FR 3005356 A1 20141107 FR 1301029 A 20130503 05 3 56 PROCEDE ET SYSTEME DE DETERMINATION D'UNE VITESSE PAR RAPPORT A L'AIR D'UN AERONEF L'invention porte sur un procédé pour déterminer une vitesse par rapport à l'air (par exemple vitesse « vraie » TAS, de l'anglais « True Air Speed », et/ou vitesse conventionnelle ou corrigée CAS, de l'anglais « Calibrated Air Speed ») d'un aéronef, ainsi que sur un système pour la mise en oeuvre d'un tel procédé. L'invention porte également sur l'utilisation de ce procédé pour détecter une défaillance d'une chaîne anémométrique (basée sur une sonde Pitot ou Pitot-statique ou un LIDAR et de capteurs de pression associés) d'un aéronef, pour pallier une telle défaillance voire pour assurer le pilotage automatique de l'aéronef. L'invention s'applique principalement aux avions, mais peut être utile aussi pour d'autres aéronefs tels que des hélicoptères ou des planeurs. La CAS est un paramètre fondamental de la conduite du vol, car elle représente la capacité d'un avion à se maintenir en vol. La CAS est une information indispensable au fonctionnement des systèmes de pilotage automatique. De manière conventionnelle, la vitesse conventionnelle est mesurée par des chaînes de capteurs et sondes embarqués, de type « Pitotou Pitot-statique » (qui fournissent également des informations d'altitude) ou LIDAR et leurs capteurs de pression associés. En cas de défaillance sur une chaîne de mesure, une information d'altitude peut être obtenue de manière indépendante, par exemple grâce à un système de navigation satellitaire (GNSS). Par contre, il n'existe pas à ce jour de solution satisfaisante pour déterminer la vitesse de l'aéronef par rapport à l'air dans une telle éventualité. Les constructeurs d'avions indiquent les procédures à suivre dans ce cas, pour permettre le pilotage manuel de l'aéronef. Le document EP 2 348 285 divulgue la possibilité de calculer la vitesse par rapport à l'air d'un aéronef à partir de la connaissance de sa vitesse par rapport au sol - obtenue grâce à un système de navigation satellitaire - et de la vitesse du vent. En cas de défaillance des chaînes anémométriques, cette dernière information peut être estimée à partir de données statistiques ou être communiquée par une station au sol. Cette solution est peu satisfaisante car l'estimée de la valeur de vitesse du vent risque d'être entachée d'incertitudes importantes. L'article de A.C. in't Veld et al. « Real-time Wind Profile Estimation using Airborne Sensors », AIAA Guidance Navigation and Control Conference, 8 - 11 août 2011, Portland, Oregon, Etats-Unis (AIAA 2011-6662) divulgue un procédé d'estimation d'un profil de vitesse du vent en fonction de l'altitude dans une région de contrôle terminale (TMA, « Terminal Maneuvering Area ») par filtrage de Kalman de données acquises par des capteurs embarqués sur des aéronefs et transmis selon le protocole ADS-B. Ce procédé ne vise pas à pallier la défaillance ou l'absence des chaînes anémométriques, mais à connaître la vitesse locale du vent à des altitudes différentes de celle à laquelle évolue l'aéronef. L'invention vise à remédier aux inconvénients précités de l'art antérieur et plus particulièrement à procurer un procédé et un système permettant de déterminer de manière fiable la vitesse d'un aéronef par rapport à l'air, même en cas de défaillance des chaînes embarquées. Ce procédé et ce système permettent également de vérifier le bon fonctionnement de ces chaînes, et de détecter leur défaillance éventuelle. Un objet de l'invention permettant d'atteindre ce but est un procédé de détermination d'une vitesse par rapport à l'air d'un aéronef, dit aéronef assisté, comprenant les étapes suivantes : a) déterminer une position dudit aéronef assisté ; b) mesurer une vitesse par rapport au sol dudit aéronef assisté ; c) recevoir une pluralité de messages en provenance d'une pluralité d'autres aéronefs, dits assistants, chaque dit message contenant : une première information, indicative d'une position d'un dit aéronef assistant, et une deuxième information, indicative d'une vitesse du vent à ladite position ; d) estimer une vitesse du vent à la position dudit aéronef assisté par interpolation des valeurs de vitesse du vent aux positions desdits aéronefs assistants obtenues lors de ladite étape c) ; et e) calculer une vitesse vraie TAS dudit aéronef assisté par différence vectorielle entre sa vitesse par rapport au sol, mesurée à ladite étape b), et la vitesse du vent estimée lors de ladite étape d). Le procédé peut comprendre également l'étape suivante : f) calculer, à partir de la vitesse vraie calculée lors de ladite étape e), au moins un paramètre choisi parmi une vitesse conventionnelle ou corrigée CAS et un nombre de Mach. Selon des modes de réalisation particuliers du procédé de l'invention : - Lesdites étapes a) et b) peuvent être mises en oeuvre au moyen d'au moins un dispositif embarqué sur ledit aéronef assisté et choisi parmi : un récepteur GNSS et une centrale inertielle. - Lesdits messages peuvent être diffusés par lesdits aéronefs assistants. En particulier, lesdits messages peuvent être diffusés par lesdits aéronefs assistants et reçus par ledit aéronef assistés selon le protocole ADS- B. - Ladite deuxième information peut comprendre une mesure de vitesse par rapport au sol et une mesure de vitesse par rapport à l'air dudit avion assistant, dont peut être déduite ladite vitesse du vent. - Au moins lesdites étapes d) et e) peuvent être réalisées par un processeur embarqué sur ledit aéronef assisté. Ladite étape d) peut être mise en oeuvre par Kriegeage. Lesdites étapes c) et d) peuvent comprendre : la réception d'un nombre n prédéfini desdits messages, l'extraction des premières et 25 deuxièmes informations contenues dans ces messages et leur stockage dans une file d'attente de longueur n ; le calcul d'une première estimation de ladite vitesse du vent par interpolation à partir des informations stockées dans ladite file d'attente ; et lors de la réception ultérieure d'autres dits messages, le stockage des informations contenues dans ladite file d'attente et la suppression 30 concomitante de ladite file d'attente d'information plus anciennes, selon le principe du premier entré - premier sorti, ou FIFO ; et le calcul d'une nouvelle estimation de ladite vitesse du vent par interpolation à partir des informations stockées dans la file d'attente ainsi modifiée. - Le procédé peut comprendre également la détermination d'un paramètre d'incertitude de ladite vitesse par rapport à l'air dudit aéronef assisté. - Lesdits messages peuvent contenir également une troisième information représentative de la pression au niveau de la mer en correspondance de la position de chaque dit aéronef assistant, auquel cas le procédé peut comporter également une étape d'estimation de la pression au niveau de la mer en correspondance de la position dudit aéronef assisté par interpolation desdites troisièmes informations. Un autre objet de l'invention est un procédé de vérification du fonctionnement d'une chaîne anémométrique embarquée sur un aéronef, dit aéronef assisté, comportant les étapes suivantes : A) mesurer une vitesse par rapport à l'air dudit aéronef assisté au moyen de ladite chaîne anémométrique ; B) déterminer ladite vitesse par rapport à l'air dudit aéronef, ainsi que son paramètre d'incertitude, par un procédé tel que décrit ci-dessus ; et C) considérer que ladite chaîne anémométrique est défaillante lorsque la différence entre ladite vitesse mesurée et ladite vitesse déterminée excède ledit paramètre d'incertitude. Encore un autre objet de l'invention est l'utilisation d'un procédé de détermination d'une vitesse par rapport à l'air d'un aéronef tel que décrit ci-dessus en cas de défaillance d'une chaîne anémométrique embarquée sur ledit aéronef assisté. Encore un autre objet de l'invention est un procédé de pilotage automatique d'un aéronef, dit assisté, comportant les étapes suivantes : i) déterminer une vitesse par rapport à l'air dudit aéronef assisté ainsi qu'un paramètre d'incertitude de ladite vitesse au moyen d'un procédé tel que décrit ci-dessus ; ii) corriger les limites supérieure et inférieure de vitesse par rapport à l'air admises pour le vol dudit aéronef en tenant compte dudit paramètre d'incertitude ; et iii) utiliser un système de pilotage automatique pour maintenir la vitesse par rapport à l'air dudit aéronef assisté à l'intérieur desdites limites supérieure et inférieure corrigées. Encore un autre objet de l'invention est un système de détermination d'une vitesse par rapport à l'air d'un aéronef, dit aéronef assisté par un procédé tel que décrit ci-dessus, comprenant : - un dispositif de détermination de la position et de la vitesse par rapport au sol dudit aéronef ; - un récepteur pour recevoir une pluralité de messages en provenance d'une pluralité d'autres aéronefs, dits assistants, chaque dit message contenant : une première information, indicative d'une position d'un dit aéronef assistant, et une deuxième information, indicative d'une vitesse du vent à ladite position ; et - un processeur configuré ou programmé pour estimer une vitesse du vent à la position dudit aéronef assisté par interpolation des valeurs de vitesse du vent aux positions desdits aéronefs assistants et pour calculer une vitesse par rapport à l'air dudit aéronef assisté par différence vectorielle entre sa vitesse par rapport au sol et la vitesse estimée du vent. Un tel système peut comprendre également une chaîne anémométrique pour déterminer une vitesse du vent ou une vitesse de l'aéronef par rapport à l'air, ainsi qu'un émetteur pour diffuser un message contenant : une première information, indicative d'une position dudit aéronef, et une deuxième information, indicative d'une vitesse du vent à ladite position. Ledit récepteur, et le cas échéant ledit émetteur, peuvent être 25 respectivement un récepteur et un émetteur ADS-B. D'autres caractéristiques, détails et avantages de l'invention ressortiront à la lecture de la description faite en référence aux dessins annexés donnés à titre d'exemple et qui représentent : - La figure 1, une illustration générale d'un procédé selon 30 l'invention ; - La figure 2, l'étape d'interpolation de la vitesse du vent ; et - Les figures 3A et 3B, respectivement, la partie émettrice et la partie réceptrice d'un système selon un mode de réalisation de l'invention. Comme illustré sur la figure 1, un procédé selon l'invention permet à un aéronef AA, qu'on appellera par la suite « aéronef assisté », de calculer une estimation de la vitesse du vent qu'il subit, et d'en déduire sa vitesse par rapport à l'air, à l'aide de messages M1 _ M4 qu'il reçoit d'autres aéronefs (dits « assistants ») AST1 - AST4. Chacun de ces messages, M1, contient essentiellement deux informations : une première information - datée - de position (latitude, longitude et altitude à un instant donné) - de l'aéronef assistant ASTi, Pi (Pià P4 dans l'exemple illustré sur la figure), que ledit aéronef assistant détermine de manière conventionnelle à l'aide d'un système de navigation satellitaire (GNSS, de l'anglais « Global Navigation Satellite System », par exemple GPS) ou d'une centrale inertielle, et une deuxième information relative à la vitesse du vent subi par ledit aéronef assistant. Cette deuxième information peut être directement une mesure de la vitesse du vent (en intensité et direction) Wi. ; c'est le cas pour les messages M1 - M3, émis par 15 les aéronefs assistants ASTI, AST2 et. AST3 En variante, comme dans le cas du message M4 émis par l'aéronef assistant AST4, il peut s'agir d'une mesure conjointe de vitesse par rapport au sol GS4, et de vitesse « vraie » par rapport à l'air, TAS4. Dans ce cas, la vitesse du vent subi par l'aéronef AST4 peut être déterminée par différence entre ces deux vecteurs 171/4 = GS4 - TAS4. La 20 deuxième information est déterminée par chaque aéronef assistant à l'aide de ses chaînes anémométriques et/ou anémobarométriques (sondes Pitot-statiques ou LIDARs et capteurs associés) et d'un système de navigation satellitaire ou inertiel. Les messages Mi peuvent contenir d'autres champs 25 d'information et notamment un identifiant de l'aéronef émetteur et éventuellement un champ de validité des données transmises (ce dernier n'est pas nécessaire si l'on prévoit que la transmission ne s'effectue que lorsque les données sont considérées valides). Les messages peuvent être diffusés directement par les 30 aéronefs, par exemple au moyen du protocole ADS-B, connu de l'art antérieur. D'autres protocoles peuvent également être utilisés, prévoyant par exemple le passage par une station relais au sol et/ou l'établissement d'une liaison entre l'aéronef émetteur (assistant) et l'aéronef récepteur (assisté). Les données transmises par les aéronefs assistants réalisent un échantillonnage du champ vectoriel représentant la vitesse et direction du vent dans une région de l'atmosphère qui entoure l'aéronef assisté AA. Cela est représenté (en deux dimensions : longitude LONG et latitude LAT, toutes deux exprimées en degrés « ° » ; plus souvent, en réalité, l'altitude sera également prise en considération) sur la figure 2, qui montre une situation où neuf aéronefs assistants ASTI - AST9 communiquent leurs vitesses (vectorielles) du vent W1- r179 à un aéronef assisté AA, situé à la position É ; sur la figure 2, chaque aéronef assistant est représenté par un cercle d'autant plus grand que l'information fournie est récent. A bord de l'aéronef assisté AA, un processeur estime la vitesse vectorielle locale du vent, W communiquées W1- W9. De manière préférée, l'interpolation est effectuée en quatre dimensions : trois dimensions spatiales et une dimension temporelle. A titre de comparaison, le vecteur WM représente la vitesse du vent mesurée par une chaîne anémométrique de l'aéronef AA ; on peut constater que l'erreur d'estimation est particulièrement faible. L'interpolation peut être calculée, de manière préférentielle, par l'algorithme de Kriegeage, connu en soi. Cependant, d'autres méthodes peuvent également être utilisées, et notamment une interpolation linéaire, une interpolation par spline cardinale, une approximation par plus proche voisin, etc. Avantageusement, l'algorithme d'interpolation (Kriegeage ou autre) peut être alimenté en données d'entrée par une file d'attente « premier entrant - premier sorti » (FIFO, de l'anglais « First In - First Out ») de longueur constante, gérée de la façon suivante : - lors d'une étape d'initialisation, les données reçues en entrée alimentent la file d'attente jusqu'à ce qu'elle atteigne une longueur définie (N éléments ou données) ; - l'algorithme d'interpolation est appliqué à ces N éléments pour donner une première estimation de la vitesse du vent ; ->EsT par interpolation des valeurs - à chaque fois qu'une nouvelle donnée est acquise, elle est intégrée à la file, tandis que l'élément le plus ancien est supprimé ; l'algorithme d'interpolation est appliqué au nouvel ensemble de N données ainsi obtenu. Le processeur connait également la vitesse de l'aéronef assisté par rapport au sol, GS, qui est déterminée à l'aide d'un récepteur GNSS et/ou d'une centrale inertielle, ou encore qui est mesurée par une station radar au sol et transmise au processeur par une liaison radio. La vitesse vraie par rapport à l'air est donc calculée par différence vectorielle entre la vitesse au sol et la vitesse du vent estimée : TAS Le plus souvent, la connaissance de la TAS n'est pas suffisante, et il est nécessaire de calculer la vitesse conventionnelle ou corrigée CAS, ce qui est connu en soi. Par exemple, on peut procéder de la façon suivante : 1) Déterminer la TAS scalaire : TAS = -NI TAS -TAS où - est l'opérateur produit scalaire. 2) Mesurer la température à l'altitude courante de l'avion ou la calculer à partir de l'atmosphère standard et d'un delta ISA préférentiellement fourni par un organisme météo. Cette température est notée T. rT 2) Calculer le nombre de Mach : M = TAS / yoù y = 1,4 et r est la constante des gaz parfaits (r= 287,053 u.s.i. - unités du système international). 3) calculerqcp = (1+ 0,2M2)3'5 -1. ; cette grandeur représente la différence entre la pression dynamique et la pression statique, rapportée à la pression au niveau de la mer. 4) calculer un facteur de correction qui vaut cf = ( 288,15 288'15 - 0 ,0019 812 Zn ) 5,256114 ) si l'altitude de pression Zp - exprimée en pieds (ft) - est inférieure à l'altitude de la tropopause Tropo_Alt (préférentiellement fournie par un organisme météo ou à défaut entrée dans un Tropo _ Alt-Zp 20804,945 système de gestion de vol, FMS) et cf = 0,223353e sinon GS W EST 5) calculer la CAS (scalaire) en noeuds (kts) en appliquant la formule suivante: 2 CAS = 661.4748\i5(1+ qcp x cf) 7 -1 Dans certains modes de réalisation on pourra s'arrêter à l'étape 2, c'est-à-dire au calcul du nombre de Mach. Avantageusement, la précision de l'estimation est connue et définie par un paramètre d'incertitude 6 tel qu'il existe une probabilité de N% (N étant une valeur prédéfinie inférieure à 100, par exemple 95% ou 99%) que la différence entre la valeur réelle de la CAS et sa valeur estimée se trouve dans une plage [-G, -Fcs]. La valeur du paramètre d'incertitude a peut être fournie par l'algorithme d'interpolation (Kriegeage ou autre), et/ou être déterminé expérimentalement. Le calcul du nombre de Mach et/ou de la CAS nécessite la connaissance de l'altitude de pression Zp de l'aéronef assistée. Ce paramètre peut être connu grâce aux capteurs embarqués, même en présence d'une défaillance sur une voie empêchant la mesure de la vitesse par rapport à l'air. D'autres modes de réalisation permettent de récupérer cette information même si elle n'est pas fournie par les capteurs embarqués. Une première possibilité consiste à calculer Zp à partir de la hauteur par rapport au sol, mesurée par GNSS, de la pression locale au niveau 2 0 de la mer, fournie par exemple par un organisme météorologique via une station au sol, et du modèle d'évolution de la pression avec l'altitude de l'atmosphère standard. Une deuxième possibilité consiste à procéder par interpolation, comme pour la détermination de la vitesse du vent. Dans ce cas, les aéronefs 25 assistants transmettent également une troisième information représentative de leur pression locale au niveau de la mer ; cette information est reçue par l'aéronef assistée et utilisée pour calculer, par interpolation (Kriegeage ou autre), la pression au niveau de la mer à la latitude et longitude dudit aéronef assisté. La troisième information peut être constituée directement par une 30 valeur de pression (« QNH »), ou par des grandeurs permettant de déterminer cette valeur via un calcul (par exemple, une hauteur par rapport au sol et une pression locale perçue par l'aéronef). Il existe plusieurs applications possibles d'un procédé selon l'invention tel que décrit ci-dessus. Premièrement, il peut être utilisé pour pallier la défaillance des chaînes anémométriques embarquées et permettre un pilotage manuel, voire automatique, de l'aéronef dans des conditions de sécurité maximale dans une telle éventualité. En ce qui concerne le pilotage automatique, le paramètre d'incertitude a de l'estimation de la CAS peut être utilisé pour appliquer une limitation supplémentaire au domaine de vol. Soient VIS et VMO, respectivement, la limite basse et haute du domaine de vol, on peut régler le système de pilotage automatique pour qu'il maintienne la vitesse (CAS) de l'aéronef, estimée conformément à l'invention, dans l'intervalle [V''.; V x ] avec V .' = Vis + 6 et V.- = VMO . Cela assure avec une probabilité de N% que l'aéronef reste dans son domaine de vol. Le procédé de l'invention peut également servir à vérifier le bon fonctionnement des chaînes anémométriques. En effet, si la vitesse CASM mesurée par lesdits capteurs diffère de la vitesse estimée conformément à l'invention de plus que a (en valeur absolue), cela signifie qu'il y une probabilité de N% que la valeur mesurée soit erronée. Dans ce cas, avantageusement, seule la vitesse estimée est utilisée pour le pilotage (automatique ou manuel) de l'appareil, et le ou les pilotes sont avertis de ce fait. Les figures 3A et 3B montrent un schéma fonctionnel d'un système de mise en oeuvre du procédé décrit ci-dessus. Dans un souci de clarté, ce système a été décomposé en un sous-système émetteur (fig. 3A), utilisé par un aéronef assistant, et un sous-système récepteur (fig. 3B), utilisé par un aéronef assisté. En général, un aéronef doit pouvoir fonctionner aussi bien en tant qu'assistant que d'assisté, par conséquent il comportera les deux sous-systèmes. Par ailleurs, ces derniers peuvent partager certains éléments (indiqués par les mêmes références sur les deux figures). Le sous-système émetteur de la figure 3A comprend un 30 processeur PROC, un système de gestion de vol FMS, une chaîne anémométrique ou anémobarométrique DAB, un récepteur GNSS ou une station inertielle RSN et un émetteur, de préférence de type ADS-B, EM. Le processeur PROC reçoit de la chaîne anémobarométrique DAB une mesure vectorielle de la vitesse par rapport à l'air CAS ainsi qu'une valeur d'altitude h ; il reçoit du dispositif RSN une information de position P et de vitesse par rapport au sol GS ; en outre, il reçoit du système de gestion de vol FMS une information DeltaISA, d'écart entre la température mesurée à l'extérieur de l'aéronef et la température de l'atmosphère standard à la même altitude. Le processeur PROC utilise ces éléments d'information pour calculer la vitesse du vent de la manière suivante : h et DeltaISA permettent de convertir la vitesse conventionnelle en vitesse vraie TAS, puis la vitesse du vent West calculée par différence vectorielle entre la vitesse au sol et la TAS : W = GS -TAS . La vitesse du vent ainsi calculée et la position obtenue directement par le dispositif RSN sont ensuite communiquées à l'émetteur EM, qui les insère dans un message Mi qu'il transmet en modalité diffusion. La transmission peut être inhibée si une défaillance de la chaîne anémobarométrique est détectée, indiquant que la vitesse du vent déterminée n'est pas correcte ou fiable. En variante, un champ « de validité » du message peut être utilisé pour indiquer si l'information transmise est correcte et fiable ou pas. Le sous-système récepteur de la figure 3B comprend, lui aussi, un processeur PROC (qui peut être le même que celui du sous-système émetteur, ou pas), un système de gestion de vol FMS, une chaîne anémométrique ou anémobarométrique DAB, un récepteur GNSS ou une station inertielle RSN, un récepteur, de préférence de type ADS-B, REC, un système de pilotage automatique P/A et/ou un système d'affichage D/U. Le récepteur REC reçoit une pluralité de messages M1 ...Mn émis par les sous-systèmes émetteurs d'autres aéronefs et en extrait les informations de position P1 ...Pn et de vitesse du vent W>n (cette dernière information pouvant être remplacée, comme expliqué plus haut, par des valeurs de vitesse par rapport au sol et de vitesse vraie par rapport à l'air), qui sont communiquées au processeur. Ce dernier utilise ces informations, avec les valeurs de position et de vitesse par rapport au sol obtenues par le dispositif RSN, pour calculer la vitesse TAS en appliquant le procédé décrit plus haut. Si les messages reçus comprennent un champ de validité, ce dernier est vérifié et seules les informations considérées comme valables sont utilisées. Le processeur reçoit également une valeur DeltaISA du système de gestion de vol TAS et une mesure d'altitude h de la chaîne anémobarométrique DAB, qui permettent de calculer la vitesse corrigée CAS ; il calcule également le paramètre d'incertitude a. En outre, il peut recevoir de la DAB une valeur mesurée de la vitesse corrigée, CASM, et la comparer à celle estimée - et au paramètre d'incertitude 6 - de manière à détecter une éventuelle défaillance de la chaîne anémobarométrique. En fonction du résultat de cette comparaison, le processeur transmet au système de pilotage automatique P/A et/ou à un système d'affichage D/U : - soit la vitesse corrigée mesurée CASM, si la DAB n'a pas été identifié comme défaillante ; - soit, dans le cas contraire, la vitesse corrigée estimée, ainsi que son paramètre d'incertitude 6 et un identifiant MD signifiant que le pilotage doit se faire en « mode dégradé », sur la base d'une vitesse estimée et non mesurée. En variante, le processeur peut communiquer avec une station au sol qui effectue les calculs, et principalement l'interpolation permettant de 2 0 déterminer la vitesse du vent. REVENDICATIONS1. Procédé de détermination d'une vitesse par rapport à l'air d'un aéronef, dit aéronef assisté (AA), comprenant les étapes suivantes : a) déterminer une position (P) dudit aéronef assisté ; b) mesurer une vitesse par rapport au sol (GS) dudit aéronef assisté ; c) recevoir une pluralité de messages (M1 - M4) en provenance d'une pluralité d'autres aéronefs, dits assistants (ASTI - AST4), chaque dit message contenant : une première information - indicative d'une position d'un dit aéronef assistant, et une deuxième information (W1 - 17173; GS4, TAS4 ), indicative d'une vitesse du vent à ladite position ; d) estimer une vitesse du vent (WEST) à la position dudit aéronef assisté par interpolation des valeurs de vitesse du vent aux positions desdits aéronefs assistants obtenues lors de ladite étape c) ; et e) calculer une vitesse vraie TAS dudit aéronef assisté par différence vectorielle entre sa vitesse par rapport au sol, mesurée à ladite étape b), et la vitesse du vent estimée lors de ladite étape d). 2. Procédé selon la revendication 1 comprenant également l'étape suivante : f) calculer, à partir de la vitesse vraie calculée lors de ladite étape e), au moins un paramètre choisi parmi une vitesse conventionnelle ou corrigée CAS et un nombre de Mach. 3. Procédé selon l'une des revendications précédentes dans lequel lesdites étapes a) et b) sont mises en oeuvre au moyen d'au moins un dispositif (RSN) embarqué sur ledit aéronef assisté et choisi parmi : un récepteur GNSS et une centrale inertielle. 4) Procédé selon l'une des revendications précédentes dans lesquels lesdits messages sont diffusés par lesdits aéronefs assistants. 5) Procédé selon la revendication 4 dans lequel lesdits messages sont diffusés par lesdits aéronefs assistants et reçus par ledit aéronef assistés selon le protocole ADS-B. 6) Procédé selon l'une des revendications précédentes dans lequel ladite deuxième information comprend une mesure de vitesse par rapport au sol ( GS4) et une mesure de vitesse par rapport à l'air (TAS4) dudit avion assistant (AST4), dont peut être déduite ladite vitesse du vent. 7) Procédé selon l'une des revendications précédentes dans lequel au moins lesdites étapes d) et e) sont réalisées par un processeur embarqué (PROC) sur ledit aéronef assisté. 8) Procédé selon l'une des revendications précédentes dans lequel ladite étape d) est mise en oeuvre par Kriegeage. 9. Procédé selon l'une des revendications précédentes dans lequel lesdites étapes c) et d) comprennent : - la réception d'un nombre N prédéfini desdits messages, l'extraction des premières et deuxièmes informations contenues dans ces messages et leur stockage dans une file d'attente de longueur N ; - le calcul d'une première estimation de ladite vitesse du vent par interpolation à partir des informations stockées dans ladite file d'attente ; et - lors de la réception ultérieure d'autres dits messages, le stockage des informations contenues dans ladite file d'attente et la suppression concomitante de ladite file d'attente d'information plus anciennes, selon le principe du premier entré - premier sorti, ou FIFO ; et le calcul d'une nouvelle estimation de ladite vitesse du vent par interpolation à partir des informations stockées dans la file d'attente ainsi modifiée. 10) Procédé selon l'une des revendications précédentes comprenant également la détermination d'un paramètre d'incertitude (a) de ladite vitesse par rapport à l'air dudit aéronef assisté. 11) Procédé selon l'une des revendications précédentes dans lequel lesdits messages contiennent également une troisième information représentative de la pression au niveau de la mer en correspondance de la position de chaque dit aéronef assistant, le procédé comportant également une étape d'estimation de la pression au niveau de la mer en correspondance de la position dudit aéronef assisté par interpolation desdites troisièmes informations. 12) Procédé de vérification du fonctionnement d'une chaîne 10 anémométrique (DAB) embarquée sur un aéronef, dit aéronef assisté, comportant les étapes suivantes : A) mesurer une vitesse par rapport à l'air (CASsAB) dudit aéronef assisté au moyen de ladite chaîne anémométrique ; B) déterminer ladite vitesse par rapport à l'air (CAS) dudit 15 aéronef, ainsi que son paramètre d'incertitude (a), par un procédé selon la revendication 9 ; et C) considérer que ladite chaîne anémométrique est défaillante lorsque la différence entre ladite vitesse mesurée et ladite vitesse déterminée excède ledit paramètre d'incertitude. 20 13) Utilisation d'un procédé selon l'une des revendications 1 à 11 en cas de défaillance d'une chaîne anémométrique embarquée sur ledit aéronef assisté. 25 14) Procédé de pilotage automatique d'un aéronef, dit assisté, comportant les étapes suivantes : i) déterminer une vitesse par rapport à l'air dudit aéronef assisté ainsi qu'un paramètre d'incertitude de ladite vitesse au moyen d'un procédé selon la revendication 12 ; 30 ii) corriger les limites supérieure et inférieure de vitesse par rapport à l'air admises pour le vol dudit aéronef en tenant compte dudit paramètre d'incertitude ; etiii) utiliser un système de pilotage automatique (P/A) pour maintenir la vitesse par rapport à l'air dudit aéronef assisté à l'intérieur desdites limites supérieure et inférieure corrigées. 15) Système de détermination d'une vitesse par rapport à l'air d'un aéronef, dit aéronef assisté par un procédé selon l'une des revendications 1 à 11, comprenant : - un dispositif (RSN) de détermination de la position et de la vitesse par rapport au sol dudit aéronef ; - un récepteur (REC) pour recevoir une pluralité de messages (M1...Me) en provenance d'une pluralité d'autres aéronefs, dits assistants, chaque dit message contenant : une première information (Pi ...P'n) indicative d'une position d'un dit aéronef assistant, et une deuxième information (W1 ... Win), indicative d'une vitesse du vent à ladite position ; et - un processeur (PROC) configuré ou programmé pour estimer une vitesse du vent à la position dudit aéronef assisté par interpolation des valeurs de vitesse du vent aux positions desdits aéronefs assistants et pour calculer une vitesse par rapport à l'air dudit aéronef assisté par différence vectorielle entre sa vitesse par rapport au sol et la vitesse estimée du vent. 16) Système selon la revendication 15, comprenant également une chaîne anémométrique (DAB) pour déterminer une vitesse du vent ou une vitesse de l'aéronef par rapport à l'air, ainsi qu'un émetteur (EM) pour diffuser un message (Mi) contenant : une première information, indicative d'une position dudit aéronef, et une deuxième information, indicative d'une vitesse du vent à ladite position. 17) Système selon l'une des revendications 15 et 16 dans lequel ledit récepteur, et le cas échéant ledit émetteur, sont respectivement un 30 récepteur et un émetteur ADS-B.