FR 3050346 A1 20171020 FR 1653447 A 20160419 1 DOMAINE TECHNIQUE La présente invention appartient au domaine de la surveillance d'aéronefs. La présente invention concerne plus particulièrement la détection de messages Mode S en squitter étendu (« Mode S extended squitter » dans la littérature anglo-saxonne), dits « messages 1090 ES », émis par des aéronefs dans le contexte de systèmes de surveillance ADS-B (« Automatic Dependent Surveillance-Broadcast » dans la littérature anglo-saxonne). La présente invention trouve une application particulièrement avantageuse, bien que nullement limitative dans le cas où plusieurs messages Mode S d'aéronefs différents, dont un ou plusieurs messages 1090 ES, sont susceptibles d'être reçus simultanément. ÉTAT DE LA TECHNIQUE Actuellement, les aéronefs émettent des messages Mode S qui sont reçus par des stations réceptrices au sol. Le nombre d'aéronefs susceptibles d'être présents, à un instant donné, dans la zone géographique desservie par une station réceptrice au sol est généralement assez limité. Par conséquent, la probabilité de collision entre messages Mode S émis par des aéronefs différents, c'est-à-dire la probabilité pour qu'au moins deux messages Mode S soient reçus en même temps par une station réceptrice au sol, est assez faible. Cependant, il est envisagé de nos jours d'effectuer la réception de messages Mode S émis par des aéronefs, en particulier des messages 1090 ES, au moyen de satellites en orbite basse (« Low Earth Orbit » ou LEO dans la littérature anglo-saxonne) ou en orbite moyenne (« Medium Earth Orbit » ou MEO dans la littérature anglo-saxonne). Du fait de son positionnement à une altitude élevée par rapport aux altitudes de vol des aéronefs, un tel satellite est adapté à desservir une zone géographique de dimensions beaucoup plus importantes que celles d'une zone géographique desservie par une station réceptrice au sol. La zone géographique desservie par un satellite est donc susceptible d'être survolée, à un instant donné, par un nombre d'aéronefs beaucoup plus important qu'une zone géographique desservie par une station réceptrice au sol. Il en résulte que la probabilité de collision au niveau d'un satellite est supérieure à la probabilité de collision au niveau d'une station réceptrice au sol. 3050346 2 La probabilité de collision au niveau d'un satellite peut en outre être très élevée (de l'ordre de 50% ou davantage), et ce d'autant plus qu'un même aéronef est susceptible d'émettre un grand nombre de messages Mode S sur des intervalles de temps relativement courts. Par exemple, un même aéronef 5 émet des messages 1090 ES de manière récurrente avec une période moyenne de 0.5 secondes (en pratique comprise entre 0.4 et 0.6 secondes), et les messages 1090 ES ne représentent qu'un type particulier des messages Mode S susceptibles d'être émis par des aéronefs. La probabilité de collision étant généralement assez faible dans le cas 10 d'une station réceptrice au sol, le protocole de communication Mode S n'a pas été conçu pour permettre de distinguer facilement des messages Mode S reçus en même temps, mais principalement pour permettre la détection et le décodage des messages Mode S avec des moyens à la fois simples et robustes, exploitant de simples mesures de puissance. 15 Toutefois, de tels moyens de détection et de décodage conventionnels ne permettent pas de détecter des messages Mode S dès lors qu'au moins deux messages Mode S sont reçus simultanément, et ne sont donc pas adaptés à la collecte de messages 1090 ES par un satellite, ou dans le cas d'une station réceptrice au sol desservant une zone géographique dense, 20 comme par exemple un aéroport. EXPOSÉ DE L'INVENTION La présente invention a pour objectif de remédier à tout ou partie des limitations des solutions de l'art antérieur, notamment celles exposées ci-avant, en proposant une solution qui permette d'améliorer la capacité de détection 25 d'un message 1090 ES en présence d'une collision avec un autre message Mode S, de type 1090 ES ou d'un autre type. A cet effet, et selon un premier aspect, l'invention concerne un procédé de détection de messages Mode S en squitter étendu, dit « message 1090 ES », dans un signal global mesuré dans une bande de fréquences 30 d'émission de messages 1090 ES par des aéronefs, comportant l'obtention, à partir du signal global, d'un signal comportant deux voies en quadrature de phase, dit « signal IQ », et la sélection d'une ou de plusieurs fréquences candidates pour la détection d'un message 1090 ES. Ledit procédé comporte 3050346 3 en outre, pour chaque fréquence candidate : - le recalage fréquentiel du signal IQ avec un signal de référence en fonction de la fréquence candidate considérée, ledit signal de référence comportant un préambule de message 1090 ES et un 5 indicateur de format de message 1090 ES, - la mesure de ressemblance entre le signal de référence et le signal 10, de sorte à obtenir plusieurs valeurs d'un signal de ressemblance associées respectivement à différentes positions dudit signal de référence par rapport audit signal IQ pour la 10 fréquence candidate considérée, un message 1090 ES étant détecté lorsqu'une valeur du signal de ressemblance vérifie un critère prédéterminé de détection. Ainsi, le procédé de détection repose notamment sur l'obtention d'un signal IQ à partir du signal global, c'est-à-dire une représentation complexe du 15 signal global en amplitude et en phase. Grâce à un tel signal 10, qui préserve l'information de phase du signal global là où les procédés de détection selon l'art antérieur ne conservaient que l'information d'amplitude/puissance, il est possible d'améliorer la détection d'un message 1090 ES en cas de collision avec un autre message Mode S, de type 1090 ES ou d'un autre type. 20 En particulier, il est alors possible d'exploiter la dimension fréquentielle pour distinguer des messages Mode S reçus en même temps. En effet, bien que les messages Mode S soient théoriquement tous émis sur la fréquence 1090 MHz, il s'avère en pratique que la fréquence réelle sur laquelle est reçu un message Mode S peut être différente de 1090 MHz et varier légèrement 25 d'un message Mode S à un autre. Notamment, l'imprécision tolérée pour la génération de la fréquence 1090 MHz étant de l'ordre de ± 1 MHz, la fréquence réelle sur laquelle est émis un message Mode S est en pratique comprise entre 1089 MHz et 1091 MHz. En outre, la fréquence réelle sur laquelle un message Mode S est reçu peut, par effet Doppler, différer de la fréquence sur laquelle 30 ledit message Mode S a été émis. En particulier, dans le cas où ce message Mode S est reçu par un satellite en orbite LEO, la variation fréquentielle introduite par effet Doppler peut être de l'ordre de ± 25 kHz. Ces perturbations, à savoir l'imprécision sur la génération de la 3050346 4 fréquence 1090 MHz et/ou l'effet Doppler deviennent cependant avantageuses dans un contexte de collisions dans la mesure où elles introduisent une diversité fréquentielle susceptible de permettre de distinguer des messages Mode S reçus en même temps. Il est par conséquent possible de considérer 5 plusieurs fréquences candidates pour la détection d'un message 1090 ES, autour de leur fréquence 1090 MHz théorique d'émission. La détection comporte alors, pour chaque fréquence candidate, le recalage fréquentiel du signal IQ avec un signal de référence, et la mesure de la ressemblance entre le signal de référence et le signal 10, pour plusieurs 10 positions dudit signal de référence par rapport audit signal 10, de sorte à obtenir des valeurs d'un signal de ressemblance. Le signal de référence reprend certains champs des messages Mode S tels que définis par le protocole de communication Mode S. En particulier, selon le protocole de communication Mode S, un message Mode S débute par 15 un préambule et est suivi par un indicateur de format. Dans le protocole de communication Mode S, le préambule est constant et le même pour tous les messages Mode S, y compris pour les messages 1090 ES. Par contre, plusieurs formats sont possibles pour les messages Mode S, et un indicateur de format égal à 17 (DF17) indique que le format utilisé est celui des messages 20 1090 ES. Le signal de référence utilisé comporte le préambule et l'indicateur de format de message 1090 ES (DF17). En présence d'un message 1090 ES dans le signal 10, le signal de ressemblance comporte en principe un optimum local lorsque le signal de référence est synchronisé avec le début dudit message 1090 ES. On comprend 25 donc qu'il est possible de détecter ce message 1090 ES à partir du signal de ressemblance. Dans des modes particuliers de mise en oeuvre, le procédé de détection de messages 1090 ES peut comporter en outre l'une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les 30 combinaisons techniquement possibles. Dans des modes particuliers de mise en oeuvre, la sélection de fréquences candidates comporte l'analyse d'un spectre fréquentiel du signal 10, une fréquence candidate correspondant à une fréquence pour laquelle la 3050346 5 valeur du spectre fréquentiel du signal IQ vérifie un critère prédéterminé de détection. Dans des modes particuliers de mise en oeuvre, le procédé de détection de messages 1090 ES comporte le recalage en phase du signal IQ 5 avec le signal de référence en fonction d'au moins une phase candidate pour la détection d'un message 1090 ES. Dans des modes particuliers de mise en oeuvre, les fréquences candidates sont des fréquences choisies arbitrairement et le signal de détection correspond à une fonction sur un espace tridimensionnel dont les 10 dimensions correspondent respectivement aux fréquences candidates, aux instants candidats et à des phases candidates. Dans des modes particuliers de mise en oeuvre, le procédé de détection de messages 1090 ES comporte l'extraction de bits inclus dans chaque message 1090 ES détecté. 15 Dans des modes particuliers de mise en oeuvre, le procédé de détection de messages 1090 ES comporte, pour chaque message 1090 ES détecté, l'extraction d'un identifiant ICAO de l'aéronef ayant émis ledit message 1090 ES détecté, et les messages 1090 ES suivants émis par ledit aéronef sont détectés au moyen d'un signal de référence augmenté d'un identifiant 20 d'émetteur égal à l'identifiant ICAO dudit aéronef. De telles dispositions permettent d'améliorer la détection de messages 1090 ES émis par un aéronef pour lequel on a préalablement détecté un message 1090 ES. En effet, du fait que le signal de référence pour cet aéronef est plus long que lors de la détection initiale, un gain de traitement plus 25 important sera introduit lors de la mesure de ressemblance. Dans des modes particuliers de mise en oeuvre, chaque message 1090 ES comportant des bits de données et des bits de détection d'erreur, l'extraction de bits d'un message 1090 ES détecté comporte une étape d'extraction initiale de bits et de calcul d'une valeur de confiance pour chaque 30 bit de données extrait, et comporte en outre les étapes suivantes qui sont itérées jusqu'à ce qu'un critère prédéterminé d'arrêt soit vérifié : - calcul, à partir des bits de données extraits, de bits de détection d'erreur théoriques, 3050346 6 - comparaison des bits de détection d'erreur théoriques avec les bits de détection d'erreur extraits du message 1090 ES détecté, - lorsque les bits de détection d'erreur théoriques sont différents des bits de détection d'erreur extraits : modification d'au moins le bit 5 extrait présentant la plus faible valeur de confiance parmi les bits extraits qui n'ont pas encore été modifiés. Dans des modes particuliers de mise en oeuvre, l'extraction des bits inclus dans chaque message 1090 ES détecté est réalisée de manière conjointe pour l'ensemble des messages 1090 ES détectés. 10 Dans des modes particuliers de mise en oeuvre, chaque message 1090 ES comportant des bits de données et des bits de détection d'erreur, l'extraction de bits comporte les étapes suivantes qui sont itérées jusqu'à ce qu'un critère prédéterminé d'arrêt soit vérifié : - extraction de bits pour chaque message 1090 ES détecté présent 15 dans le signal 10, - calcul, pour chaque message 1090 ES, de bits de détection d'erreur théoriques à partir des bits de données extraits, - comparaison, pour chaque message 1090 ES, des bits de détection d'erreur théoriques avec les bits de détection d'erreur 20 extraits, les bits extraits pour un message 1090 ES étant considérés comme corrects lorsque les bits de détection d'erreur théoriques sont égaux aux bits de détection d'erreur extraits et comme erronés lorsque les bits de détection d'erreur théoriques sont différents des bits de détection d'erreur extraits, 25 - lorsque les bits extraits sont corrects pour au moins un message 1090 ES et erronés pour au moins un autre message 1090 ES : o reformation de chaque message 1090 ES pour lequel les bits extraits sont corrects, à partir desdits bits extraits, o suppression de chaque message 1090 ES reformé de sorte 30 à obtenir un signal IQ dépourvu des messages 1090 ES pour lesquels les bits extraits sont corrects. Selon un second aspect, la présente invention concerne un produit programme d'ordinateur comportant un ensemble d'instructions de code de 3050346 7 programme qui, lorsqu'elles sont exécutées par un processeur, configurent ledit processeur pour mettre en oeuvre un procédé de détection de messages 1090 ES selon l'un quelconque des modes de mise en oeuvre de l'invention. Selon un troisième aspect, la présente invention concerne un dispositif 5 de détection de messages 1090 ES dans un signal global mesuré dans une bande de fréquences d'émission de messages 1090 ES par des aéronefs, comportant des moyens configurés pour mettre en oeuvre les étapes d'un procédé de détection de messages 1090 ES selon l'un quelconque des modes de mise en oeuvre de l'invention. 10 Selon un quatrième aspect, la présente invention concerne un satellite comportant un dispositif de mesure d'un signal global mesuré dans une bande de fréquences d'émission de messages 1090 ES par des aéronefs, et un dispositif de détection de messages 1090 ES l'un quelconque des modes de réalisation de l'invention. 15 Selon un cinquième aspect, la présente invention concerne un système de détection de messages 1090 ES, comportant un satellite comportant un dispositif de mesure d'un signal global mesuré dans une bande de fréquences d'émission de messages 1090 ES par des aéronefs, ledit système comportant en outre des moyens, répartis entre ledit satellite et une 20 ou plusieurs stations de traitement au sol, configurés pour mettre en oeuvre les étapes d'un procédé de détection de messages 1090 ES selon l'un quelconque des modes de mise en oeuvre de l'invention. PRÉSENTATION DES FIGURES L'invention sera mieux comprise à la lecture de la description suivante, 25 donnée à titre d'exemple nullement limitatif, et faite en se référant aux figures qui représentent : - Figure 1 : une représentation schématique d'un exemple de réalisation d'un système de détection de messages 1090 ES, - Figure 2: un diagramme illustrant les principales étapes d'un 30 exemple de mise en oeuvre d'un procédé de détection de messages 1090 ES, - Figure 3 : une représentation temporelle d'un exemple de signal de référence pour la détection de messages 1090 ES, 3050346 8 - Figure 4: un diagramme illustrant les principales étapes d'un premier mode particulier de mise en oeuvre d'une étape d'extraction de bits de messages 1090 ES détectés, - Figure 5: un diagramme illustrant les principales étapes d'un 5 second mode particulier de mise en oeuvre d'une étape d'extraction de bits de messages 1090 ES détectés. Dans ces figures, des références identiques d'une figure à une autre désignent des éléments identiques ou analogues. Pour des raisons de clarté, les éléments représentés ne sont pas à l'échelle, sauf mention contraire. 10 DESCRIPTION DÉTAILLÉE DE MODES DE RÉALISATION La figure 1 représente schématiquement un exemple de réalisation d'un système 10 de détection de messages 1090 ES émis par des aéronefs 40. Tel qu'illustré par la figure 1, le système 10 de détection comporte au moins un satellite 20 et au moins une station 30 de traitement au sol. Le 15 satellite 20 est de préférence placé sur une orbite non-géosynchrone, telle qu'une orbite LEO ou une orbite MEO. Le satellite 20 comporte un dispositif de mesure (non représenté sur les figures) adapté à mesurer un signal global dans une bande de fréquences d'émission de messages 1090 ES par des aéronefs 40, c'est-à-dire une bande 20 de fréquences comportant la fréquence 1090 MHz qui correspond à la fréquence théorique d'émission des messages 1090 ES par les aéronefs 40. De manière générale, le système 10 de détection comporte des moyens, répartis entre ledit satellite 20 et une ou plusieurs stations 30 de traitement au sol, configurés pour mettre en oeuvre les étapes d'un procédé 50 25 de détection de messages 1090 ES, qui sera décrit ci-après. En d'autres termes, suivant le mode de réalisation considéré : - toutes les étapes du procédé 50 de détection sont mises en oeuvre par le satellite 20, auquel cas la station 30 de traitement au sol se contente de recevoir dudit satellite 20 le résultat de la détection, 30 - toutes les étapes du procédé 50 de détection sont mises en oeuvre par une ou plusieurs stations 30 de traitement au sol, auquel cas le satellite 20 se contente de transmettre à une station 30 de traitement au sol le signal global mesuré, 3050346 9 - le satellite 20 met en oeuvre une partie des étapes du procédé 50 de détection, et la ou les étapes non mises en oeuvre par le satellite 20 sont mises en oeuvre par une ou plusieurs stations 30 de traitement au sol. 5 Dans la suite de la description, on se place de manière non limitative dans le cas où toutes les étapes du procédé 50 de détection sont mises en oeuvre par un dispositif de détection (non représenté sur les figures) embarqué dans le satellite 20. Le dispositif de détection comporte par exemple un ou plusieurs 10 processeurs et des moyens de mémorisation (disque dur magnétique, mémoire électronique, disque optique, etc.) dans lesquels est mémorisé un produit programme d'ordinateur, sous la forme d'un ensemble d'instructions de code de programme à exécuter pour mettre en oeuvre les différentes étapes du procédé 50 de détection de messages 1090 ES. Alternativement ou en 15 complément, le dispositif de détection peut comporter un ou des circuits logiques programmables (FPGA, PLD, etc.), et/ou un ou des circuits intégrés spécialisés (ASIC), et/ou un ensemble de composants électroniques discrets, etc., adaptés à mettre en oeuvre tout ou partie desdites étapes du procédé de détection de messages 1090 ES. 20 En d'autres termes, le dispositif de détection est configuré de façon logicielle (produit programme d'ordinateur spécifique) et/ou matérielle (FPGA, PLD, ASIC, composants électroniques discrets, etc.) pour mettre en oeuvre les différentes étapes du procédé 50 de détection de messages 1090 ES. La figure 2 représente les principales d'un exemple de mise en oeuvre 25 d'un procédé 50 de détection de messages 1090 ES. Tel qu'illustré par la figure 2, le procédé 50 de détection comporte tout d'abord une étape 51 d'obtention, à partir du signal global, d'un signal comportant deux voies en quadrature de phase, dit « signal IQ », et une étape 52 de sélection d'une ou de plusieurs fréquences candidates pour la détection 30 d'un message 1090 ES. Le procédé 50 de détection comporte ensuite, pour chaque fréquence candidate, des étapes de : - 53 recalage fréquentiel du signal IQ avec un signal de référence en fonction de la fréquence candidate considérée, 3050346 10 - 54 mesure de ressemblance entre le signal de référence et le signal IQ de sorte à obtenir, pour la fréquence candidate considérée, plusieurs valeurs d'un signal de ressemblance associées respectivement à différentes positions dudit signal de 5 référence par rapport audit signal IQ. Le procédé 50 de détection de messages 1090 ES comporte également une étape 55 de détection d'un message 1090 ES lorsqu'une valeur du signal de ressemblance vérifie un critère prédéterminé de détection. Il est à noter que le procédé 50 de détection peut se contenter de 10 détecter les messages 1090 ES, sans extraire les bits inclus dans les messages 1090 ES détectés. Dans des modes préférés de mise en oeuvre, et tel qu'illustré par la figure 2, le procédé 50 de détection peut cependant comporter une étape 56 d'extraction des bits inclus dans les messages 1090 ES détectés. 15 On décrit à présent des exemples non limitatifs de mise en oeuvre des différentes étapes du procédé 50 de détection de messages 1090 ES. A) Etape 51 d'obtention de signal IQ L'étape 51 d'obtention de signal IQ peut mettre en oeuvre toute méthode connue de l'homme de l'art permettant d'obtenir deux voies I et Q en 20 quadrature de phase. Le signal IQ est par exemple un signal dans lequel la fréquence 1090 MHz d'émission des messages 1090 ES a été ramenée sur une fréquence FO qui correspond à la fréquence nulle (bande de base) ou à une fréquence intermédiaire non nulle prédéterminée. En considérant que K messages 1090 ES sont présents dans le signal 25 global, reçus d'aéronefs différents, alors le signal IQ est de préférence un signal numérique qui peut par exemple être exprimé sous la forme suivante : K n j- 2-rrfkn+(j)k (1) r[n] 1Ak - Sk -5k - e Fs + o[n] k=1 Fs expression dans laquelle : - Fs correspond à la fréquence d'échantillonnage du signal 10, 30 - r[n] correspond au signal IQ à l'instant n/Fs, - sk correspond au message 1090 ES reçu de l'aéronef de rang k, 1 3050346 11 L expression dans laquelle L correspond à une valeur seuil positive prédéterminée. Si le critère prédéterminé de détection est vérifié pour la fréquence fek, alors cela signifie que le signal IQ est susceptible de comporter un message 1090 ES sur la fréquence fek qui est sélectionnée comme 10 fréquence candidate. Le cas échéant, la phase de la valeur Dfreq[fek] correspond à une estimation (pek de la phase dudit message 1090 ES détecté, dite « phase candidate » dudit message 1090 ES détecté. En pratique, le spectre fréquentiel Dfreq est de préférence calculé pour un ensemble de fréquences prédéterminées, dites « fréquences d'analyse ». 15 Le cas échéant, une fréquence candidate sélectionnée correspond à l'une des fréquences d'analyse. Il est cependant possible, dans des modes particuliers de mise en oeuvre, d'effectuer une interpolation du spectre fréquentiel 1Dfreq12 au voisinage d'une fréquence d'analyse ayant permis de vérifier le critère prédéterminé de détection. De telles dispositions permettent d'améliorer la 20 précision de la sélection des fréquences candidates. L'interpolation du spectre fréquentiel IDfreq12 repose par exemple sur une fonction polynomiale d'ordre 2. C) Etape 53 de recalage fréquentiel Pour chaque fréquence candidate sélectionnée, le procédé 50 de détection comporte alors une étape 53 au cours de laquelle le signal IQ est 25 recalé en fréquences avec un signal de référence, en fonction de la fréquence candidate considérée. Par « recalage fréquentiel », on entend que le signal IQ et le signal de référence sont ramenés sur la même fréquence. Dans la suite de la description, on se place de manière non limitative dans le cas où le recalage fréquentiel, pour une fréquence candidate fek, 30 correspond à ramener en bande de base un message 1090 ES reçu sur ladite fréquence candidate fek, c'est-à-dire à ramener ledit message 1090 ES sur la fréquence nulle. Pour la fréquence candidate fek, le signal IQ est par exemple recalé en fréquences selon l'expression suivante : 3050346 13 - j2-rrfek-n (3) rek [ri] = r[n] - e Fs De préférence, si l'étape 52 de sélection de fréquences candidates a également fourni une phase candidate (pek du message 1090 ES détecté sur la fréquence candidate fek, le signal IQ est recalé selon l'expression suivante 5 -j. 211-fek n ±(Pek (4) rek [ri] = r[n] - e Fs En d'autres termes, le signal IQ est également recalé en phase avec le signal de référence en fonction de la phase candidate (pek. De manière plus générale, si les phases candidates n'ont pas été obtenues au cours de l'étape 52 de sélection de fréquences candidates, par 10 exemple du fait que les fréquences candidates sont des fréquences choisies arbitrairement autour de la fréquence F0, il est possible de considérer plusieurs phases candidates pour chaque fréquence candidate, et de déterminer plusieurs signaux IQ recalés en phase pour chaque fréquence candidate, en fonction des différentes phases candidates sélectionnées. 15 D) Etape 54 de mesure de ressemblance avec le signal de référence Après recalage fréquentiel en fonction d'une fréquence candidate, on mesure la ressemblance entre un signal de référence et le signal 10, pour plusieurs positions dudit signal de référence par rapport audit signal 10. Lesdites positions correspondent respectivement à plusieurs instants candidats 20 pour la détection d'un message 1090 ES dans le signal 10. Le protocole de communication Mode S définit plusieurs champs à l'intérieur d'un message Mode S. En particulier, selon le protocole de communication Mode S, un message Mode S débute par un préambule et est suivi par un indicateur de format (« Downlink Format » ou DF). 25 Dans le protocole de communication Mode S, le préambule est constant et le même pour tous les messages Mode S, y compris pour les messages 1090 ES. Par contre, plusieurs formats sont possibles pour les messages Mode S, et la durée d'un message Mode S et l'interprétation des champs présents 30 dans un message Mode S peuvent dépendre du format utilisé. L'indicateur de format permet donc de connaître la durée dudit message Mode S qui, suivant 3050346 14 le format utilisé, peut comporter 56 microsecondes (ps) ou 112 ps d'informations après le préambule, ce qui correspond à respectivement 56 bits ou 112 bits. L'indicateur de format permet également de savoir comment interpréter les bits extraits des différents champs du message Mode S. 5 Dans le cas d'un message 1090 ES, l'indicateur de format est égal à 17 (DF17) et ledit message 1090 ES comporte 112 ps d'informations après le préambule, c'est-à-dire 112 bits. Le tableau ci-après donne les différents champs d'un message 1090 ES (format DF17), tels que définis par le protocole de communication 10 Mode S, selon leur ordre d'apparition dans un message 1090 ES. Champ préambule DF CA/A F AA ME PI Durée / 8 ps 5 bits 3 bits 24 bits 56 bits 24 bits Bits Ainsi, un message 1090 ES comporte, après le préambule, un champ DF de 5 bits qui correspond à l'indicateur de format, égal à 17 (DF17) dans le cas d'un message 1090 ES. Un message 1090 ES comporte ensuite : 15 - un champ CA/AF de 3 bits, utilisé pour transmettre des informations de contrôle complémentaires, - un champ AA de 24 bits, qui comporte un identifiant d'émetteur du message 1090 ES, qui correspond à l'identifiant ICAO (« International Civil Aviation Organization ») de l'aéronef 40 qui 20 émet le message 1090 ES, - un champ ME de 56 bits qui comporte des informations relatives notamment à la position ou à l'altitude de l'aéronef 40 qui émet le message 1090 ES, - un champ PI de 24 bits, qui comporte des bits de détection d'erreur 25 permettant de détecter, lors de l'extraction, des erreurs dans les 88 bits de données extraits correspondant aux bits des champs DF, CA/AF, AA et ME. Le signal de référence correspond au début d'un message 1090 ES, et comporte, outre le préambule, un indicateur de format de message 1090 ES, 30 c'est-à-dire un indicateur de format égal à 17 (DF17). 3050346 15 La figure 3 représente schématiquement le signal de référence en bande de base. Tel qu'illustré par la figure 3, le préambule tel que défini par le protocole de communication Mode S comporte 4 impulsions carrées, chacune 5 de durée 0.5 ps, qui débutent respectivement à 0 ps, 1 ps, 3.5 ps et 4.5 ps. Selon le protocole de communication Mode S, chaque bit est encodé par un symbole de durée Ts égale à 1 ps. Si le bit est égal à 0, on émet un symbole 50 qui est nul sur la première moitié de la durée Ts, constant et positif sur la seconde moitié de la durée Ts. Si le bit est égal à 1, on émet un symbole 10 51 qui est positif et constant sur la première moitié de la durée Ts, nul sur la seconde moitié de la durée Ts. Tel qu'illustré sur la figure 3, le signal de référence comporte, après le préambule, un indicateur de format égal à 17, c'est-à-dire correspondant à la séquence de bits « 1 0 0 0 1 » qui est encodée par la séquence de symboles 15 « 51 50 50 50 51 ». Avec un tel signal de référence, on comprend que, si le signal IQ comporte un message 1090 ES qui débute à un instant candidat ôek, le signal de ressemblance doit comporter un optimum local lorsque le signal de référence est positionné audit instant candidat ôek, puisque le signal IQ doit, 20 sur la durée dudit signal de référence, ressembler audit signal de référence. Il est à noter que toute fonction permettant de mesurer la ressemblance entre deux signaux différents, comportant chacun plusieurs échantillons, peut être mise en oeuvre au cours de l'étape 54 de mesure de ressemblance. Dans la suite de description, on se place dans le cas où la 25 mesure de ressemblance comporte une corrélation entre le signal de référence et le signal 10, qui est calculée pour plusieurs positions du signal de référence par rapport au signal 10, c'est-à-dire pour plusieurs instants candidats. On se place tout d'abord dans le cas où le signal IQ rek est recalé à la fois en fréquence et en phase en fonction d'une fréquence candidate et d'une 30 phase candidate, selon l'expression (4) ci-dessus, par exemple suite à la détection de fréquences candidates par analyse spectrale du signal 10. Dans un tel cas, le signal de ressemblance Dtiming, pour une fréquence candidate fek et un instant candidat ôek, est par exemple calculé selon l'expression suivante : 3050346 16 (M-1 Dtiming, fek, (pek (Elek ) = real rek [m + Elek ] - sref [m] 0-1=o expression dans laquelle real(x) correspond à la valeur réelle du nombre complexe x, et M correspond au nombre d'échantillons du signal de référence sref. Si l'on désigne par OSF le facteur de sur-échantillonnage par rapport à la 5 durée Ts (1 ps), c'est-à-dire le nombre d'échantillons que comporte le signal IQ par durée Ts, alors le nombre d'échantillons M du signal de référence sref, de durée égale à (8 + 5) ps, est égal à 13.05F. Dans le cas, par exemple, où les fréquences candidates sont des fréquences choisies arbitrairement sans analyse spectrale du signal 10, alors il 10 est possible de considérer, pour chaque fréquence candidate, plusieurs phases candidates possibles. Dans un tel cas, le signal IQ rek recalé est par exemple donné par l'expression (3) ci-dessus, et le signal de ressemblance Dtiming comporte, pour chaque fréquence candidate fek et chaque instant candidat ôek considéré, plusieurs valeurs associées respectivement à différentes phases 15 candidates (pek, calculées par exemple selon l'expression suivante : (m-i Dti min g (fek, ôek, (Pek ) rr1=0 Dans certains cas, il est possible de considérer un signal de référence plus long que le signal de référence illustré par la figure 3. En particulier, certains champs ne varient pas d'un message 1090 ES à un autre émis par un 20 même aéronef 40. C'est le cas du champ AA qui comporte l'identifiant ICAO de l'aéronef 40 qui a émis le message 1090 ES, l'identifiant ICAO d'un aéronef 40 étant un identifiant unique et permanent dudit aéronef 40. Par conséquent, il est possible, dans des modes préférés de mise en oeuvre, pour un aéronef 40 donné pour lequel on a préalablement détecté un message 1090 ES et extrait 25 l'identifiant ICAO dudit message 1090 ES détecté, de rajouter au signal de référence un champ AA égal à l'identifiant ICAO dudit aéronef 40 pour la détection des messages 1090 ES suivants émis par cet aéronef 40. Le cas échéant, le champ CA/AF, entre le champ DF et le champ AA, est par exemple forcé à zéro. Le champ AA comportant 24 bits, il est possible d'améliorer 30 fortement le gain de traitement introduit par la mesure de ressemblance lors de (5) real rek [m + Elekl- sref [m] - e-i(Pek (6) 3050346 17 la détection des messages 1090 ES suivants émis par cet aéronef 40. E) Etape 55 de détection de message 1090 ES Ainsi, à l'issue de l'étape 54 de mesure de ressemblance, on dispose de valeurs de signal de ressemblance associées aux différentes fréquences 5 candidates, et pour chaque fréquence candidate, aux différents instants candidats et, le cas échéant, aux différentes phases candidates. Au cours de l'étape 55 de détection de message 1090 ES, on évalue les différentes valeurs du signal de ressemblance, et un message 1090 ES est détecté lorsqu'une valeur du signal de ressemblance vérifie un critère 10 prédéterminé de détection. Dans le cas où une analyse spectrale du signal IQ a été effectuée au cours de l'étape 52 de sélection de fréquences candidates, on vérifie, pour chaque fréquence candidate susceptible de comporter un message 1090 ES, si la valeur du signal de ressemblance pour l'un des instants candidats permet de 15 confirmer la détection d'un message 1090 ES pour la fréquence candidate considérée. Par exemple, dans le cas de l'expression (5), le critère de détection est considéré comme vérifié par la valeur du signal de ressemblance à un instant candidat ôek si ladite valeur correspond à un maximum local du signal de ressemblance suivant la dimension des instants candidats, et si en 20 outre ladite valeur est supérieure à une valeur seuil prédéterminée. Lorsque le critère de détection est vérifié pour une valeur du signal de ressemblance, les paramètres de détection du message 1090 ES détecté, c'est-à-dire les estimations de la fréquence, de la phase et de l'instant de réception dudit message 1090 ES, correspondent à : 25 - la fréquence candidate et à la phase candidate, déterminées au cours de l'étape 53 de sélection, associées à la valeur dudit signal de ressemblance pour laquelle le critère de détection a été vérifié, - l'instant candidat associé à la valeur dudit signal de ressemblance pour laquelle le critère de détection a été vérifié. 30 Dans le cas où les fréquences candidates sont choisies arbitrairement, et dans le cas où plusieurs phases candidates sont considérées, le signal de ressemblance, par exemple donné par l'expression (6) ci-dessus, correspond à une fonction sur un espace tridimensionnel dont les dimensions correspondent 3050346 18 respectivement aux fréquences candidates, aux phases candidates et aux instants candidats. Par exemple, dans le cas de l'expression (6), le critère de détection est considéré comme vérifié par une valeur du signal de ressemblance pour un triplet donné (fréquence candidate, phase candidate, 5 instant candidat) si ladite valeur correspond à un maximum local du signal de ressemblance, et si en outre ladite valeur est supérieure à une valeur seuil prédéterminée. La détection de maximum local du signal de détection met par exemple en oeuvre des méthodes conventionnelles de traitement de volume (par exemple filtres d'ordre, analyse différentielle (laplacien, hessien), etc.). 10 Lorsque le critère de détection est vérifié pour une valeur du signal de ressemblance, le triplet associé donne les différents paramètres de détection du message 1090 ES détecté, c'est-à-dire des estimations de la fréquence, de la phase et de l'instant de réception dudit message 1090 ES. Ainsi, à l'issue de l'étape 55 de détection, un ou plusieurs messages 15 1090 ES peuvent avoir été détectés et, pour chaque message 1090 ES détecté, on peut également estimer un triplet de paramètres de détection correspondant à des estimations de la fréquence, de la phase et de l'instant de réception dudit message 1090 ES détecté. F) Etape 56 d'extraction de bits de données 20 L'étape 56 d'extraction vise à extraire les bits des messages 1090 ES détectés, c'est-à-dire à décoder chacun des messages 1090 ES détectés. Plusieurs méthodes d'extraction des bits des messages 1090 ES détectés peuvent être envisagées, et on comprend que le choix d'une méthode d'extraction particulière ne constitue qu'une variante d'implémentation de 25 l'invention. Notamment, il est possible soit d'extraire séparément les bits des différents messages 1090 ES détectés, soit d'extraire conjointement lesdits bits des messages 1090 ES détectés. F.1) Extraction séparée des bits des messages 1090 ES détectés Tel qu'indiqué ci-dessus, à l'issue de l'étape 55 de détection, on 30 dispose, pour chaque message 1090 ES détecté, d'un triplet de paramètres de détection qui correspondent à des estimations de la fréquence, de la phase et de l'instant de réception dudit message 1090 ES détecté. Pour chaque message 1090 ES détecté, il est donc possible de 3050346 19 calculer un signal IQ recalé en fréquence et en phase, à l'intérieur duquel on connaît l'instant de début du préambule du message 1090 ES détecté et les instants de transitions entre symboles. Pour extraire les bits du message 1090 ES détecté de rang k, pour 5 lequel le triplet de paramètres de détection est (fdk, (pdk, ôdk), on calcule par exemple un signal d'extraction selon l'expression suivante : ( r2Trfdk 11 Fs-Fek j rk [n] = real r[n]-e Les premiers échantillons des différents symboles du message 1090 ES détecté correspondent aux échantillons rk[b'dk + m.OSF], m étant un 10 entier positif ou nul, expression dans laquelle ô'dk correspond au début du champ DF du message 1090 ES détecté, égal à (Eiclk + 8.05F). Pour extraire les différents bits, il est par exemple possible de corréler les échantillons du signal d'extraction rk, qui correspondent à un même symbole, avec les symboles SO et 51 décrits en référence à la figure 3. Par exemple, il est 15 possible de calculer, pour chaque symbole du signal d'extraction rk, des valeurs COk et Cl k selon les expressions suivantes : OSF-1 OSF-1 COk (M) = rk [5' dk +m - OSF + il - SO[i] = rk dk m - OSF + il - SO [i] i=o i=OSF/2 OSF-1 OSF/2-1 C1k(m)= rk dk +m - OSF + il - S1[i] = rk dk +m - OSF + il - S1[i] i=o i=o Par conséquent, si COk(m) est supérieur à Cl k (m) alors le symbole de 20 rang m est considéré comme étant un symbole SO, de sorte qu'un bit égal à 0 est extrait. Dans le cas contraire (Ci k (m) > COk(m)), le symbole de rang m est considéré comme étant un symbole 51, de sorte qu'un bit égal à 1 est extrait. Les valeurs COk(m) et Ci k(m) peuvent également être utilisées pour calculer une valeur de confiance Vck(m) pour le bit extrait pour le symbole de 25 rang m. Par exemple, une telle valeur de confiance Vck(m) peut être calculée selon l'expression suivante : Vck(m) = IC0k(m) - C1 k(m)I D'autres méthodes d'extraction séparée de bits peuvent être 3050346 20 cependant mises en oeuvre, et on comprend que le choix d'une méthode particulière ne constitue qu'une variante d'implémentation de l'invention. F.2) Extraction conjointe des bits des messages 1090 ES détectés Dans des modes particuliers de mise en oeuvre, l'extraction des bits 5 inclus dans chaque message 1090 ES détecté est réalisée de manière conjointe pour l'ensemble des messages 1090 ES détectés, par exemple par une estimation au sens du maximum de vraisemblance (« Joint Maximum Likelihood Sequence Estimation » ou JMLSE dans la littérature anglo-saxonne) ou une estimation au sens des moindres carrés, etc. 10 En effet, à l'issue de l'étape 55 de détection, on dispose du nombre de messages 1090 ES présents dans le signal 10, mais également des différents paramètres de détection pour chaque message 1090 ES détecté. En estimant par ailleurs, selon toute méthode connue de l'homme de l'art, l'amplitude de chaque message 1090 ES détecté, il est alors possible de calculer un signal IQ 15 théorique pour chaque combinaison possible de séquences de bits possibles pour les différents messages 1090 ES détectés. La combinaison de séquences de bits pour les messages 1090 ES détectés qui permet de maximiser la vraisemblance, ou de minimiser l'écart quadratique moyen entre le signal IQ théorique et le signal IQ obtenu à partir du signal global, donne alors 20 directement tous les bits de tous les messages 1090 ES détectés. F.3) Amélioration des performances de l'extraction de bits On décrit ci-après, en référence aux figures 4 et 5, deux exemples de méthodes permettant d'améliorer les performances de l'extraction, qui peuvent être considérées seules ou en combinaison. Ces méthodes s'appuient 25 notamment sur le fait que chaque message 1090 ES comporte des bits de données et des bits de détection d'erreur (« Cyclic Redundancy Check » ou CRC dans la littérature anglo-saxonne). La figure 4 représente schématiquement les principales étapes d'un mode préféré de mise en oeuvre de l'étape 56 d'extraction de bits des 30 messages 1090 ES détectés, qui s'appuie notamment sur le calcul de valeurs de confiance associées respectivement aux différents bits extraits. Tel qu'illustré par la figure 4, l'extraction de bits comporte tout d'abord, pour chaque message 1090 ES détecté, une étape 570 d'extraction initiale de 3050346 21 bits, au cours de laquelle on calcule des valeurs de confiance associées respectivement aux différents bits extraits. L'extraction de bits comporte ensuite, pour chaque message 1090 ES détecté, les étapes suivantes qui sont itérées jusqu'à ce qu'un critère prédéterminé d'arrêt soit vérifié : 5 - 571 calcul, à partir des bits de données extraits, de bits de détection d'erreur théoriques, - 572 comparaison des bits de détection d'erreur théoriques avec les bits de détection d'erreur extraits du message 1090 ES détecté, - lorsque les bits de détection d'erreur théoriques sont différents des 10 bits de détection d'erreur extraits : 573 modification d'au moins le bit extrait présentant la plus faible valeur de confiance parmi les bits extraits qui n'ont pas encore été modifiés. L'étape 571 est conventionnelle et vise à calculer des bits de détection d'erreur théoriques à partir des bits de données extraits au niveau du dispositif 15 de détection, selon la même méthode de calcul que celle utilisée à l'émission pour calculer les bits de détection d'erreur inclus dans le message 1090 ES. Les bits de détection d'erreur théoriques sont ensuite comparés avec les bits de détection d'erreur extraits du message 1090 ES détecté, au cours de l'étape 572 de comparaison. 20 Lorsque les bits de détection d'erreur théoriques sont différents des bits de détection d'erreur extraits du message 1090 ES détecté considéré (référence 572a sur la figure 4), cela signifie qu'au moins un bit extrait (bit de données ou bit de détection d'erreurs) est erroné. Le cas échéant, la plus faible valeur de confiance calculée est en principe associée au bit extrait qui a le plus 25 de chances d'être celui qui est erroné. En modifiant la valeur du bit extrait ayant la valeur de confiance la plus faible (étape 573 de modification), il est alors possible de calculer à nouveau de bits de détection d'erreur théoriques (étape 571 de calcul) à partir des bits de données extraits, et de les comparer aux bits de détection d'erreur extraits (étape 572 de comparaison) afin de 30 vérifier si cette modification a permis de supprimer les erreurs. Si tel n'est pas le cas, il est possible de modifier la valeur du bit extrait non encore modifié ayant la valeur de confiance la plus faible, et ainsi de suite jusqu'à ce que ce que le critère d'arrêt soit vérifié. 3050346 22 Lorsque, au cours de l'étape 572 de comparaison, les bits de détection d'erreur théoriques sont égaux aux bits de détection d'erreur extraits du message 1090 ES détecté (référence 572b sur la figure 4), alors les bits extraits sont considérés comme corrects, et le critère d'arrêt est considéré 5 comme vérifié. Dans un tel cas, les bits de données extraits correspondent à ceux, obtenus éventuellement après modification d'un ou de plusieurs desdits bits de données extraits, ayant permis d'obtenir des bits de détection d'erreur théoriques égaux aux bits de détection d'erreur extraits. De préférence, le critère d'arrêt comporte également un nombre 10 maximal de modifications autorisées, ledit nombre maximal de modifications autorisées étant par exemple inférieur à quinze. Ainsi, dans l'exemple non limitatif illustré par la figure 4, l'étape 56 d'extraction de bits comporte une étape 574 de détermination si le nombre de bits modifiés a atteint le nombre maximal de modifications autorisées. Si c'est le cas (référence 574b sur la 15 figure 4), le critère d'arrêt est vérifié et l'extraction des bits prend fin même s'il est considéré que les bits extraits comportent encore des erreurs. Si le nombre de bits modifiés n'a pas encore atteint le nombre maximal de modifications autorisées (référence 574a sur la figure 4), alors la valeur du bit non encore modifié ayant la valeur de confiance la plus faible est modifiée, et ainsi de suite 20 jusqu'à ce que ce que le critère d'arrêt soit vérifié. La figure 5 représente schématiquement les principales étapes d'un autre mode préféré de mise en oeuvre de l'étape 56 d'extraction de bits des messages 1090 ES détectés. Tel qu'illustré par la figure 5, l'extraction de bits comporte les étapes 25 suivantes qui sont itérées jusqu'à ce qu'un critère prédéterminé d'arrêt soit vérifié : - 580 extraction de bits de données et de bits de détection d'erreur pour chaque message 1090 ES détecté présent dans le signal 10, - 581 calcul, pour chaque message 1090 ES, de bits de détection 30 d'erreur théoriques à partir des bits de données extraits, - 582 comparaison, pour chaque message 1090 ES, des bits de détection d'erreur théoriques avec les bits de détection d'erreur extraits, 3050346 23 - lorsque les bits extraits sont corrects pour au moins un message 1090 ES et erronés pour au moins un autre message 1090 ES : o 583 reformation de chaque message 1090 ES pour lequel les bits extraits sont corrects, à partir desdits bits extraits, 5 o 584 suppression de chaque message 1090 ES reformé de sorte à obtenir un signal IQ dépourvu des messages 1090 ES pour lesquels les bits extraits étaient corrects. Au cours de l'étape 580, on extrait les bits de données et les bits de détection d'erreur de chaque message 1090 ES détecté présent dans le signal 10 10. Au cours de la première itération, tous les messages 1090 ES détectés sont encore présents dans le signal 10. Par contre, au cours d'une seconde itération, le signal IQ ne comporte plus tous les messages 1090 ES initialement détectés, puisqu'un ou plusieurs messages 1090 ES ont été supprimés du signal 10. Par conséquent, à partir de la seconde itération et pour toutes les 15 éventuelles itérations suivantes, l'étape 580 ne porte que sur l'extraction des bits inclus dans les messages 1090 ES initialement détectés qui n'ont pas encore été reformés et supprimés dans le signal 10. Pour chaque message 1090 ES dont les bits ont été extraits, des bits de détection d'erreur théoriques sont calculés (étape 581 de calcul), à partir 20 des bits de données extraits au niveau du dispositif de détection, et lesdits bits de détection d'erreur théoriques sont comparés (étape 582 de comparaison) avec les bits de détection d'erreur extraits pour le message 1090 ES considéré. Lorsque les bits extraits sont corrects pour au moins un message 1090 ES et erronés pour au moins un autre message 1090 ES (référence 582a 25 sur la figure 5), alors chaque message 1090 ES pour lequel les bits extraits sont corrects est reformé (étape 583 de reformation), par le dispositif de détection, selon le protocole de communication 1090 ES. La reformation du message 1090 ES tient compte en outre des différents paramètres de détection obtenus pour ce message 1090 ES, afin d'obtenir un message 1090 ES 30 reformé qui correspond au message 1090 ES tel qu'il apparaît théoriquement dans le signal 10, compte tenu des bits extraits et des différents paramètres de détection. Chaque message 1090 ES reformé est ensuite soustrait (étape 584 de suppression) au signal 10. Le signal IQ ainsi obtenu, théoriquement 3050346 24 dépourvu des messages 1090 ES pour lesquels les bits extraits étaient corrects, est ensuite utilisé pour itérer à nouveau les étapes décrites ci-dessus. Lorsque les bits extraits sont corrects pour tous les messages 1090 ES initialement détectés dans le signal IQ (référence 582b sur la figure 5 5), le critère d'arrêt est alors considéré comme vérifié. De préférence, le critère d'arrêt comporte également un nombre maximal d'itérations autorisées, par exemple inférieur à quinze. Ainsi, dans l'exemple non limitatif illustré par la figure 5, l'étape 56 d'extraction de bits comporte une étape 585 de détermination si le nombre d'itérations effectuées, 10 c'est-à-dire le nombre de suppressions successives de messages 1090 ES reformés, a atteint le nombre maximal d'itérations autorisées. Si c'est le cas (référence 585b sur la figure 5), le critère d'arrêt est vérifié et l'extraction des bits prend fin même s'il est considéré que les bits extraits comportent encore des erreurs pour certains messages 1090 ES initialement détectés. Si le 15 nombre d'itérations effectuées n'a pas encore atteint le nombre maximal d'itérations autorisées (référence 585a sur la figure 5), alors il est possible de reformer et de supprimer des messages 1090 ES si les bits extraits sont corrects pour au moins un message 1090 ES détecté et erronés pour au moins un autre message 1090 ES détecté. 20 De manière plus générale, il est à noter que les modes de mise en oeuvre et de réalisation considérés ci-dessus ont été décrits à titre d'exemples non limitatifs, et que d'autres variantes sont par conséquent envisageables. Notamment, l'invention a été décrite en considérant principalement un système de détection de messages 1090 ES comportant au moins un satellite. 25 Rien n'exclut, suivant d'autres exemples, d'avoir un système de détection dépourvu de satellites et d'intégrer, dans une station réceptrice au sol, un dispositif de détection de messages 1090 ES selon l'un quelconque des modes de réalisation de l'invention. En effet, l'invention trouve également une application particulièrement avantageuse dans le cas d'une station réceptrice 30 au sol desservant une zone géographique dense, comme par exemple un aéroport, pour laquelle la probabilité de collision entre messages 1090 ES peut également être importante. Il est également possible, suivant d'autres exemples, d'avoir un 3050346 25 système de détection comportant à la fois un ou plusieurs satellites, pour la réception de messages 1090 ES depuis une orbite terrestre, et une ou plusieurs stations réceptrices au sol, pour la réception de messages 1090 ES au sol sans avoir à passer par l'un des satellites, intégrant un dispositif de 5 détection selon l'un quelconque des modes de réalisation de l'invention. REVENDICATIONS1 - Procédé (50) de détection de messages Mode S en squitter étendu, dit « message 1090 ES », dans un signal global mesuré dans une bande de fréquences d'émission de messages 1090 ES par des aéronefs (40), caractérisé en ce qu'il comporte l'obtention (51), à partir du signal global, d'un signal comportant deux voies en quadrature de phase, dit « signal IQ », et la sélection (52) d'une ou de plusieurs fréquences candidates pour la détection d'un message 1090 ES, ledit procédé comportant en outre, pour chaque fréquence candidate : - le recalage fréquentiel (53) du signal IQ avec un signal de référence en fonction de la fréquence candidate considérée, ledit signal de référence comportant un préambule de message 1090 ES et un indicateur de format de message 1090 ES, - la mesure de ressemblance (54) entre le signal de référence et le signal 10, de sorte à obtenir plusieurs valeurs d'un signal de ressemblance associées respectivement à différentes positions dudit signal de référence par rapport audit signal IQ pour la fréquence candidate considérée, un message 1090 ES étant détecté (55) lorsqu'une valeur du signal de ressemblance vérifie un critère prédéterminé de détection. 2 - Procédé (50) selon la revendication 1, dans lequel la sélection (52) de fréquences candidates comporte l'analyse d'un spectre fréquentiel du signal 10, une fréquence candidate correspondant à une fréquence pour laquelle la valeur du spectre fréquentiel du signal IQ vérifie un critère prédéterminé de détection. 3 - Procédé (50) selon l'une des revendications précédentes, comportant le recalage en phase du signal IQ avec le signal de référence en fonction d'au moins une phase candidate pour la détection d'un message 1090 ES. 4 - Procédé (50) selon la revendication 1, dans lequel les fréquences candidates sont des fréquences choisies arbitrairement et le signal de détection correspond à une fonction sur un espace tridimensionnel dont les dimensions correspondent respectivement aux fréquences candidates, aux instants candidats et à des phases candidates pour la détection d'un 3050346 27 message Mode S. 5 - Procédé (50) selon l'une des revendications précédentes, comportant l'extraction (56) de bits inclus dans chaque message 1090 ES détecté. 6 - Procédé (50) selon la revendication 5, comportant, pour chaque message 5 1090 ES détecté, l'extraction d'un identifiant ICAO de l'aéronef (40) ayant émis ledit message 1090 ES détecté, et dans lequel les messages 1090 ES suivants émis par ledit aéronef sont détectés au moyen d'un signal de référence augmenté d'un identifiant d'émetteur égal à l'identifiant ICAO dudit aéronef. 10 7 - Procédé (50) selon l'une des revendications 5 à 6, dans lequel, chaque message 1090 ES comportant des bits de données et des bits de détection d'erreur, l'extraction de bits d'un message 1090 ES détecté comporte une étape (570) d'extraction initiale de bits et de calcul d'une valeur de confiance pour chaque bit de données extrait, et comporte en 15 outre les étapes suivantes qui sont itérées jusqu'à ce qu'un critère prédéterminé d'arrêt soit vérifié : - (571) calcul, à partir des bits de données extraits, de bits de détection d'erreur théoriques, - (572) comparaison des bits de détection d'erreur théoriques avec 20 les bits de détection d'erreur extraits du message 1090 ES détecté, - lorsque les bits de détection d'erreur théoriques sont différents des bits de détection d'erreur extraits : (573) modification d'au moins le bit extrait présentant la plus faible valeur de confiance parmi les bits extraits qui n'ont pas encore été modifiés. 25 8 - Procédé (50) selon l'une des revendications 5 à 7, dans lequel, chaque message 1090 ES comportant des bits de données et des bits de détection d'erreur, l'extraction de bits comporte les étapes suivantes qui sont itérées jusqu'à ce qu'un critère prédéterminé d'arrêt soit vérifié : - (580) extraction de bits pour chaque message 1090 ES détecté 30 présent dans le signal IQ, - (581) calcul, pour chaque message 1090 ES, de bits de détection d'erreur théoriques à partir des bits de données extraits, - (582) comparaison, pour chaque message 1090 ES, des bits de 3050346 28 détection d'erreur théoriques avec les bits de détection d'erreur extraits, les bits extraits pour un message 1090 ES étant considérés comme corrects lorsque les bits de détection d'erreur théoriques sont égaux aux bits de détection d'erreur extraits et 5 comme erronés lorsque les bits de détection d'erreur théoriques sont différents des bits de détection d'erreur extraits, - lorsque les bits extraits sont corrects pour au moins un message 1090 ES et erronés pour au moins un autre message 1090 ES : o (583) reformation de chaque message 1090 ES pour lequel 10 les bits extraits sont corrects, à partir desdits bits extraits, o (584) suppression de chaque message 1090 ES reformé de sorte à obtenir un signal IQ dépourvu des messages 1090 ES pour lesquels les bits extraits sont corrects. 9 - Produit programme d'ordinateur caractérisé en ce qu'il comporte un 15 ensemble d'instructions de code de programme qui, lorsqu'elles sont exécutées par un processeur, configurent ledit processeur pour mettre en oeuvre un procédé (50) de détection de messages 1090 ES selon l'une des revendications précédentes. 10 - Dispositif de détection de messages 1090 ES dans un signal global 20 mesuré dans une bande de fréquences d'émission de messages 1090 ES par des aéronefs (40), caractérisé en ce qu'il comporte des moyens configurés pour mettre en oeuvre les étapes d'un procédé (50) de détection selon l'une des revendications 1 à 8. 11 - Satellite (20) comportant un dispositif de mesure d'un signal global 25 mesuré dans une bande de fréquences d'émission de messages 1090 ES par des aéronefs (40), caractérisé en ce qu'il comporte un dispositif de détection de messages 1090 ES selon la revendication 10. 12 - Système (10) de détection de messages 1090 ES, comportant un satellite (20) comportant un dispositif de mesure d'un signal global mesuré dans 30 une bande de fréquences d'émission de messages 1090 ES par des aéronefs (40), caractérisé en ce que ledit système comporte des moyens, répartis entre ledit satellite (20) et une ou plusieurs stations (30) de traitement au sol, configurés pour mettre en oeuvre les étapes d'un 3050346 29 procédé (50) de détection selon l'une des revendications 1 à 8.