Procédé d’estimation de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes (2) mis en œuvre par ordinateur, le procédé comprenant des étapes de : - estimation (202) d’une matrice de covariance de signaux acquis par le réseau d’antennes (2), les signaux étant issus de l’onde électromagnétique, - calcul (204) d’un vecteur propre normalisé ( ) de la matrice de covariance, le vecteur propre étant associé à une valeur propre de la matrice de covariance, - corrélation (205) du vecteur propre normalisé ( ) avec une première table de référence ( ) de sorte à produire un premier spectre de corrélation ( ) comprenant des indices de corrélation respectivement associés à différentes directions d’arrivée, la première table de référence ( ) comprenant des réponses normalisées du réseau d’antennes (2) à une première onde électromagnétique de référence polarisée dans une première direction respectivement associées à différentes directions d’arrivée de la première onde électromagnétique de référence, - corrélation (206) du vecteur propre normalisé ( ) avec une deuxième table de référence ( ) de sorte à produire un deuxième spectre de corrélation ( ) comprenant des indices de corrélation respectivement associés à différentes directions d’arrivée, la deuxième table de référence ( ) comprenant des réponses normalisées du réseau d’antennes (2) à une deuxième onde électromagnétique de référence polarisée dans une deuxième direction respectivement associées à différentes directions d’arrivée de la deuxième onde électromagnétique de référence, la deuxième direction étant différente de la première direction, - construction (312) d’une troisième table de référence ( ) par combinaison linéaire de la première table de référence ( ) et de la deuxième table de référence ( ) avec respectivement une composante de polarisation (P1) de l’onde électromagnétique dans la première direction et une composante de polarisation (P2) de l’onde électromagnétique dans la deuxième direction, - corrélation (314) du vecteur propre normalisé ( ) avec la troisième table de référence ( ) de sorte à produire un troisième spectre de corrélation ( ) comprenant des indices de corrélation associés à différentes directions d’arrivée, - identification (316) d’une direction d’arrivée associée à un indice de corrélation maximal du troisième spectre de corrélation ( ). Procédé selon la revendication 1, comprenant en outre des étapes de : - sélection (210) d’une direction d’arrivée associée à un indice de corrélation maximal du premier spectre de corrélation ( ) ou du deuxième spectre de corrélation ( ), - construction (304) d’une table de recherche en polarisation à partir d'une première réponse normalisée du réseau d’antennes (2) à la première onde électromagnétique de référence, la première réponse normalisée étant associée dans la première table de référence ( ) à la direction d’arrivée sélectionnée, et à partir d’une deuxième réponse normalisée du réseau d’antennes (2) à la deuxième onde électromagnétique de référence, la deuxième réponse normalisée étant associée dans la deuxième table de référence ( ) à la direction d’arrivée sélectionnée, - corrélation (306) du vecteur propre normalisé ( ) avec la table de recherche en polarisation, de sorte à obtenir un quatrième spectre de corrélation ( ) comprenant des indices de corrélation associés à différentes valeurs d’une variable complexe, - ajustement (310) de la composante de polarisation de l’onde électromagnétique dans la deuxième direction (P2) à une valeur complexe (p2) associée à un indice de corrélation maximal du quatrième spectre de corrélation ( ), puis normalisation d’un vecteur de polarisation formé par les composantes de polarisation (P1) de l’onde électromagnétique dans la première direction et la composante de polarisation (P2) de l’onde électromagnétique dans la deuxième direction, la construction (312) de la troisième table de référence ( ) étant mise en œuvre après la normalisation du vecteur de polarisation. Procédé selon la revendication 2, dans lequel le premier spectre de corrélation ( ) comprend un premier indice de corrélation maximal et deuxième spectre de corrélation ( ) comprend un deuxième indice de corrélation maximal, et dans lequel la direction d’arrivée sélectionnée à l’étape de sélection (210) est : - une direction d’arrivée associée au premier indice de corrélation maximal dans le premier spectre de corrélation ( ), lorsque le premier indice de corrélation maximal est supérieur au deuxième indice de corrélation maximal, - une direction d’arrivée associée au deuxième indice de corrélation maximal dans le deuxième spectre de corrélation ( ), lorsque le premier indice de corrélation maximal n’est pas supérieur au deuxième indice de corrélation maximal. Procédé selon l’une des revendications 2 et 3, comprenant la création (302) d’un espace de recherche comprenant différentes valeurs pour la variable complexe, et dans lequel la construction (304) de la table de recherche en polarisation comprend une combinaison linéaire de la première réponse et de la deuxième réponse avec respectivement la composante de polarisation de l’onde électromagnétique dans la première direction et un poids dépendant du module de la deuxième composante de polarisation de l’onde électromagnétique et d’une des différentes valeurs pour la variable complexe, la combinaison linéaire étant répétée pour chacune des différentes valeurs. Procédé selon la revendication 4, comprenant les étapes de : - mise à jour de la direction d’arrivée sélectionnée à la valeur de la direction d’arrivée identifiée, - mise à jour de l’espace de recherche en un nouvel espace de recherche incluant la composante de polarisation (P2) de l’onde électromagnétique dans la deuxième direction du vecteur de polarisation obtenue à l’issue de l’étape de normalisation, puis - répétition des étapes de sélection (303), construction (304) d’une table de recherche en polarisation, corrélation (306) du vecteur propre normalisé ( ) la table de recherche en polarisation, ajustement (310), construction (312) d’une troisième table de référence , corrélation (314) du vecteur propre normalisé ( ) avec la troisième table de référence ( ), et identification (316) d’une direction d’arrivée. Procédé selon la revendication 5, dans lequel l’espace de recherche forme une spirale dans le plan complexe, et la mise à jour de l’espace de recherche comprend une dilatation de la spirale, de sorte que le nouvel espace de recherche forme une nouvelle spirale dans le plan complexe qui est dilatée par rapport à la spirale. Procédé selon l’une des revendications 5 et 6, comprenant en outre une comparaison entre l’indice de corrélation maximal du troisième spectre de corrélation et un seuil prédéfini, ladite répétition n’étant mise en œuvre que si l’indice de corrélation maximal du troisième spectre de corrélation est inférieur au seuil prédéfini. Procédé selon l’une des revendications précédentes, dans lequel la première direction et la deuxième direction sont orthogonales. Mémoire (8) lisible par ordinateur stockant des instructions exécutables par l’ordinateur pour l’exécution des étapes du procédé selon l’une des revendications précédentes. Dispositif d’estimation (4) de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes (2), le dispositif d’estimation (4) comprenant au moins un processeur (6) configuré pour mettre en œuvre le procédé selon l’une des revendications 1 à 8. Système (1) comprenant : - un réseau d’antennes (2) configuré pour acquérir des signaux issus d’une onde électromagnétique, - un dispositif d’estimation (4) de la direction d’arrivée de l’onde électromagnétique au réseau d’antennes (2) sur la base des signaux, le dispositif d’estimation (4) étant conforme à la revendication 10. Porteur mobile, tel qu’un aéronef, comprenant un système (1) selon la revendication 11.