FR 3011377 A1 20150403 FR 1359515 A 20131001 PROCEDE DE LOCALISATION D'UNE SOURCE SONORE ET ROBOT HUMANOIDE UTILISANT UN TEL PROCEDE L'invention porte sur un procédé pour la localisation spatiale d'une source sonore, ainsi que sur un robot humanoïde pouvant mettre en oeuvre et utiliser un tel procédé. La localisation spatiale d'une source sonore est nécessaire dans des nombreuses applications notamment, mais pas exclusivement, en robotique humanoïde. Un robot peut être qualifié d'humanoïde dès lors qu'il possède certains attributs de l'apparence et des fonctionnalités de l'homme : une tête, un tronc, deux bras, éventuellement deux jambes, etc. Généralement, on demande à un robot humanoïde de pouvoir interagir avec des êtres humains de façon aussi « naturelle » que possible, en s'apercevant de la présence d'une personne, en comprenant son langage, en engageant une conversation avec elle, etc. Une capacité de localisation de sources sonores est très utile, voire nécessaire, pour atteindre ce but. En effet, une telle capacité peut permettre à un robot humanoïde de déterminer la direction dont provient un son et tourner sa tête vers cette direction ; si le son a été produit par une personne, le robot peut alors activer un logiciel de reconnaissance des visages, configurer de manière optimale un système de reconnaissance vocale, suivre « du regard » les déplacements de cette personne, etc. Plusieurs procédés et systèmes de localisation spatiale d'une ou plusieurs sources sonores sont connus de l'art antérieur. Ces procédés et systèmes se basent généralement sur une pluralité de microphones peu ou pas directifs et sur un traitement numérique des signaux captés par lesdits microphones. L'article de J. DiBiase et al. « Robust localization in reverberant rooms » dans l'ouvrage « Microphone Arrays: Signal Processing Techniques and Applications » publié sous la direction de M. S. Brandstein et D. B. Ward 3 0 par Springer-Verlag, 2001, Berlin, Allemagne, décrit trois approches principales à la localisation d'une source sonore. - Une première approche utilise des techniques d'estimation spectrale, basées sur la matrice de corrélation des signaux captés par les 3 0 1 1 3 7 7 2 microphones. Les méthodes basées sur cette approche tendent à être sensibles aux erreurs de modèles et très demandeuses de puissance de calcul. Elles conviennent principalement pour des signaux à bande étroite. - Une deuxième approche se base sur l'estimation des 5 décalages temporels entre les signaux sonores reçus par des paires de microphones (TDOA, pour « Time Difference Of Arrivai », c'est-à-dire « différence de temps d'arrivée »). Ces estimations sont utilisées, avec la connaissance des positions des microphones, pour calculer des courbes hyperboliques, dont l'intersection donne la position de la source. Les décalages 10 temporels peuvent notamment être estimés par la méthode dite PHAT-GCC (pour « PHAse Transform - Generalized Cross-Correlation », ou « transformation de phase - intercorrélation généralisée ») qui exploite le calcul d'une intercorrélation - ou corrélation croisée - entre signaux préalablement « blanchi » par filtrage. La méthode PHAT-GCC est décrite plus en détail dans 15 l'article de Ch. H. Knapp et G. C. Carter « The Generalized Correlation Method for Estimation of Time Delay », IEEE Transaction on Acoustics, Speech and Signal Processing, Vol. ASSP-24, No. 4, août 1976 pp. 320 - 327. Ces méthodes sont légère computationnellement, mais elles ne sont pas robustes au bruit corrélé provenant de sources multiples et sons sujettes aux « faux 20 positifs ». En outre, elles sont peu robustes à la réverbération, à l'exception de la méthode PHAT-GCC. - Une troisième approché consiste à synthétiser un faisceau acoustique orientable en additionnant les signaux captés par les différents microphones auxquels a été appliqué un décalage temporel variable, et à 25 identifier l'orientation du faisceau qui maximise la puissance du signal composite ainsi reçu. Les méthodes basées sur cette approche tendent à être peu robustes à la réverbération et au bruit, sauf certaines variantes qui sont cependant très demandeuses de puissance de calcul. L'article décrit plus particulièrement un procédé combinant la synthèse d'un faisceau acoustique orientable et l'intercorrélation généralisée avec transformation de phase. Cette méthode est dénommée SRP-PHAT (pour « Steered Response Power - PHAse Transform » ou « puissance de réponse dirigée - transformation de phase »). Par rapport à la méthode PHAT-GCC, elle est plus robuste au bruit mais plus sensible à la réverbération. L'invention vise à procurer un procédé de localisation d'une source sonore présentant des meilleures propriétés d'immunité à la fois au bruit et à la réverbération par rapport aux procédés connus de l'art antérieur, tout en étant suffisamment léger du point de vue computationnel pour être mise en oeuvre dans un système embarqué tel qu'un robot humanoïde, en temps réel et en parallèle d'autres tâches. Conformément à l'invention ce but est atteint par un procédé 10 dérivé de la méthode SRP-PHAT précitée et par un robot humanoïde comportant des moyens de mise en oeuvre de ce procédé. Un objet de l'invention est donc un procédé de localisation d'une source sonore comportant les étapes suivantes : a) capter des signaux sonores issus d'une source sonore à 15 localiser au moyen d'un ensemble d'au moins trois microphones ; b) sélectionner au moins trois paires de microphones dudit ensemble et, pour chaque dite paire, calculer une intercorrélation généralisée des signaux sonores captés, ledit calcul étant effectué pour une pluralité de valeurs d'un retard - dit différence de temps interauriculaire - entre lesdits 20 signaux sonores ; c) à partir desdites intercorrélations généralisées, calculer une puissance de réponse dirigée exprimée en fonction d'un vecteur des différences de temps interauriculaires pour chaque dite paire de microphones ; d) déterminer le vecteur des différences de temps 25 interauriculaires qui maximise ladite puissance de réponse dirigée ; et e) estimer une direction de localisation de ladite source sonore en fonction du vecteur des différences de temps interauriculaires déterminé lors de ladite étape d), caractérisé en ce que : 3 0 - lesdites étapes c) et d) sont mises en oeuvre en considérant une pluralité de vecteurs des différences de temps interauriculaires formant un ensemble qui comprend : un premier sous-ensemble de vecteurs compatibles avec des signaux sonores provenant d'une source sonore unique à distance infinie desdits microphones ; et un deuxième sous-ensemble de vecteurs qui ne sont pas compatibles avec des signaux sonores provenant d'une source unique à distance infinie desdits microphones ; en ce que : - chaque vecteur dudit premier sous-ensemble est associé à une direction de localisation de la source sonore unique correspondante, et chaque vecteur dudit deuxième sous-ensemble est associé à la direction de localisation associée à un vecteur dudit premier sous-ensemble qui lui est plus proche selon une métrique prédéfinie ; et en ce que : - la direction estimée lors de ladite étape e) est celle associée au vecteur des différences de temps interauriculaires déterminé lors de ladite étape d). Selon des caractéristiques avantageuses d'un tel procédé : - Les intercorrélations généralisées calculées lors de ladite étape b) peuvent être du type GCC-PHAT, correspondant à l'intercorrélation de signaux préalablement blanchis. - Lors de ladite étape b), on peut sélectionner toutes les paires de microphones dudit ensemble. - Ladite puissance de réponse dirigée peut être donnée par, ou proportionnelle à, la somme des intercorrélations généralisées pour lesdites paires de microphones. - Ledit premier sous-ensemble peut être constitué par les vecteurs dudit ensemble dont les composantes satisfont la condition dite de Chasles : tij-Etjk=tik où les indices i#j#k représentent différent microphones dudit ensemble et 'Ci], tik et tik représentent les différences de temps interauriculaires pour les paires (i,j), (j,k) et (i,k) respectivement. - Ledit deuxième sous-ensemble peut être constitué par les vecteurs dudit ensemble n'appartenant pas audit premier sous-ensemble et dont la distance à un vecteur dudit premier sous-ensemble, déterminée selon ladite métrique prédéfinie et après normalisation, est inférieure à un seuil prédéfini. - Ladite métrique prédéfinie peut être une métrique euclidienne. - Ladite étape e) peut être mise en oeuvre au moyen d'une table de correspondance associant une direction de localisation à chaque valeur dudit ensemble. Un autre objet de l'invention est un robot humanoïde comprenant : - un ensemble d'au moins trois microphones; et - un processeur connecté auxdits microphones pour recevoir en entrée des signaux sonores captés par ces derniers et programmé ou configuré pour mettre en oeuvre un procédé selon l'une des revendications précédentes. Avantageusement, lesdits microphones peuvent être agencés sur une surface supérieure d'une tête dudit robot. 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, respectivement : - la figure 1, un diagramme fonctionnel illustrant la mise en oeuvre d'un procédé selon un mode de réalisation de l'invention ; - la figure 2, la disposition des microphones sur la tête d'un robot humanoïde selon un mode de réalisation de l'invention ; et - Les figures 3A, 4A, 5A, 6A et 7A, des graphiques illustrant les performances d'un procédé selon un mode de réalisation de l'invention ; et - Les figures 3B, 4B, 5B, 6B et 7B, données à titre de comparaison, des graphiques illustrant les performances d'un procédé selon l'art antérieur. La figure 1 illustre très schématiquement la mise en oeuvre d'un procédé selon un mode de réalisation de l'invention au moyen d'un système comprenant un réseau de microphones M1 - M4 et un processeur PR. 3 0 Dans le cas de la figure, le réseau comprend quatre microphones, mais leur nombre peut être inférieur (trois au minimum) ou supérieur ; en outre, bien que les microphones soient représentés alignés, leur agencement spatial peut être quelconque. De préférence, on choisira un agencement minimisant les cas de masquage, c'est-à-dire les cas dans lesquels un obstacle s'interpose entre les microphones et les régions d'espace où peut se situer une source sonore à localiser. Un système réel devrait comprendre également des circuits de conditionnement des signaux issus par les microphones, assurant leur amplification, filtrage, échantillonnage et conversion analogique/numérique. Ces circuits sont omis dans un souci de simplicité et lisibilité de la figure. Dans la suite on considèrera donc que les microphones M1 - M4 fournissent des signaux au format numérique. Le processeur PR est programmé ou configuré pour réaliser les opérations suivantes : 1. Calcul d'une intercorrélation généralisée des signaux issus des différents microphones, et notamment d'une intercorrélation généralisée de type GCC-PHAT. Soient si(t) et si(t) les signaux issus des microphones Mi et Mi (i,j=1 - 4 ; i#j). L'intercorrélation généralisée Rsi,si de ces signaux est définie comme l'intercorrélation entre ces signaux après un filtrage préalable ; il s'agit d'une fonction d'une variable T, homogène à un temps et représentative d'un décalage temporel entre les signaux issus des deux microphones (différence de temps interauriculaire). Avantageusement, l'intercorrélation généralisée est calculée dans le domaine de la fréquence : Rsi,,;(r) = f+(Pi (w) Si(w)crii(w)Si(w)ei" et o (1) où Si,i(co) est la transformée de Fourier du signal si,i(t),« * » est l'opérateur de conjugaison complexe et c(w) est la fonction de transfert d'un filtre. Dans le cas de l'intercorrélation généralisée GCC-PHAT on choisit un filtre « blanchissant » les signaux, c'est-à-dire égalisant l'intensité de toutes leurs composantes spectrales pour ne conserver que l'information de phase : i,PHAT = iSi(60)1 1 (2) Bien entendu, dans une réalisation pratique du procédé les 3 0 transformées de Fourier seront des transformées de Fourier discrètes, calculées notamment par l'algorithme de transformée de Fourier rapide (FFT) et les intégrales seront remplacés par des sommes d'un nombre fini de termes. Ainsi, comme illustré sur la figure 1, les signaux issus des microphones M1 - M4 sont convertis dans le domaine fréquentiel (blocs FT réalisant une transformée de Fourier discrète, notamment en utilisant l'algorithme de transformée de Fourier rapide ou « FFT »), subissent un seuillage pour supprimer les composantes spectrales dont l'intensité est inférieure à une valeur prédéfinie afin d'éviter que des composantes spectrales contenant principalement ou exclusivement du bruit ne soient amplifiées par le filtre blanchissant (bloc de seuillage spectral SCS) puis sont filtrées par ledit filtre blanchissant PHAT - mais d'autres types de filtrage peuvent être utilisés, voir l'article précité de Ch. H. Knapp et G. C. Carter. Ensuite, les signaux filtrés sont multipliés deux à deux, et leurs produits sont reconvertis dans le domaine temporel (blocs FT-1 réalisant une transformée de Fourier discrète inverse, notamment en utilisant l'algorithme de transformée de Fourier rapide ou « FFT »). On obtient ainsi les intercorrélations généralisées des signaux sonores captés par les différents microphones, considérés deux à deux. Il a été choisi de prendre toutes les paires de microphones - au nombre de N(N-1)/2 si N est le nombre de microphones - mais on aurait pu se limiter à considérer trois paires choisies parmi elles, ou toute valeur comprise entre 3 et N(N-1)/2. Le passage par le domaine fréquentiel n'est pas indispensable 2 0 pour calculer les intercorrélations, mais est néanmoins très avantageux. Les calculs étant effectués de façon discrète, chaque intercorrélation généralisée n'est calculée que pour un ensemble discret de valeurs de la différence de temps interauriculaire. 2. Calcul d'une puissance de réponse dirigée SRP (en 25 anglais « Steered Response Power »). Ce calcul est effectué simplement en additionnant les différentes intercorrélations généralisées : SRP (T) = (3) où la somme est effectuée sur M paires de signaux si,si associés à des paires de microphones respectives (31\ffl(N-1)/2). La variable vectorielle i a M 30 composantes, correspondant chacune à la différence de temps interauriculaire pour une dite paire. Plus généralement, on pourrait définir la SRP comme une combinaison linéaire des intercorrélations généralisées. Il convient de noter que toutes les valeurs de T - et donc toutes les combinaisons de valeurs ne sont pas « physiquement possibles ». En effet, si on considère le cas d'une seule source sonore suffisamment éloignée du réseau de microphones pour que les ondes acoustiques parvenant à ces derniers puissent être considérées planes, la direction de localisation de ladite source est complètement identifiée par deux différences de temps interauriculaires. En d'autres termes, deux composantes du vecteur T déterminent de manière univoque (au moins en théorie) les valeurs des autres composantes. 3. Maximisation de la puissance de réponse dirigée L'identification de la direction de localisation de la source sonore se fait en maximisant la fonction SRP(r) par rapport à la variable vectorielle T. Plusieurs méthodes numériques peuvent être utilisées pour résoudre ce problème d'optimisation ; on peut citer à titre d'exemples non limitatifs les algorithmes de gradient et les méthodes de « force brute ». Cette maximisation pourrait se faire sous contrainte, en prenant en compte uniquement les valeurs de T qui sont « physiquement possibles » au sens explicité plus haut, c'est-à-dire compatibles avec des signaux sonores provenant d'une source unique à distance « infinie ». Cette approche est connue, en particulier, de l'article précité de J. H. DiBiase et al. Cependant, en pratique, elle s'avère peu robuste. En effet, à cause du bruit, des effets de masquage (obstacles interposés entre la source et un ou plusieurs microphones) et, surtout, des effets de réverbération, il est fréquent que le maximum absolu de SRP(r) corresponde à une valeur « impossible » de T. Conformément à l'invention, afin d'améliorer la robustesse de l'algorithme de localisation spatiale de la source, la recherche du maximum de SRP(r) ne se fait pas uniquement sur les valeurs « possibles » de T, mais également sur des valeurs théoriquement non admissibles, c'est-à-dire incompatibles avec des signaux sonores provenant d'une source unique à distance « infinie ». Plus précisément, comme indiqué sur la figure 1, la maximisation (bloc fonctionnel « MAX ») de la fonction SRP(r) est effectuée en considérant des valeurs du vecteur i qui forment un ensemble E constitué de deux sous-ensembles : - Un premier sous-ensemble El de vecteurs i compatibles avec des signaux sonores provenant d'une source sonore unique à distance infinie du réseau de microphones. Du point de vue mathématique, ces vecteurs ont des composantes qui satisfont aux conditions dites « de Chasles » tij-Etik=tik (4) où les indices i#j#k représentent différent microphones dudit ensemble et tij, tik et tik représentent les différences de temps interauriculaires 10 pour les paires (i,j), (j,k) et (i,k) respectivement. Ces conditions sont nécessaires, mais pas suffisantes, pour qu'un vecteur i soit compatible avec des signaux sonores provenant d'une source sonore unique à distance infinie. - Un deuxième sous-ensemble E2 de vecteurs i qui ne sont pas compatibles avec des signaux sonores provenant d'une source sonore 15 unique à distance infinie du réseau de microphones. En particulier, ce deuxième sous-ensemble est constitué par les vecteurs dudit ensemble n'appartenant pas audit premier sous-ensemble et dont la distance à un vecteur dudit premier sous-ensemble, déterminée selon une métrique prédéfinie (notamment euclidienne) et après normalisation (chaque vecteur est divisé par 20 sa norme), est inférieure à un seuil prédéfini. La détermination de ce seuil constitue un réglage du procédé, qui peut être effectué expérimentalement. Pour sélectionner les vecteurs i appartenant à l'ensemble E est les attribuer E aux sous-ensembles El et E2 on procède de la façon suivante. Premièrement, il convient d'écrire, sous forme matricielle, le 25 système d'équation qui lie le vecteur unitaire pointant dans la direction de la source, noté x, au vecteur des différencse de temps interauriculaires i : T=AX (5) où A est une matrice Mx3 (on rappelle que M est le nombre de composants du vecteur t). La matrice A n'étant pas carrée, elle ne peut pas 30 être inversée directement. Si les microphones ne sont pas coplanaires, la matrice A est de rang 3. On peut alors choisir trois lignes linéairement indépendants pour construire une matrice carrée inversible À ; on note --c" le vecteur de dimension 3 obtenu en prenant les trois composants de t correspondant auxdits composants indépendants de la matrice A. On a alors : ''c" = Âx (6) et donc x = À-1'i (7) Ensuite, on considère un ensemble de départ E' de vecteurs T, qui correspond à une sphère ou un cube dans II: (car on impose une valeur maximale admissible au module de t ou à chacune de ses composantes) discrétisée (car on ne considère que des valeurs discrètes de ces 10 composantes). Pour chaque vecteur t de E', on calcule sa projection -i dans le sous-espace des vecteurs qui vérifient les relations de Chasles de dimensions 3. Si la différence di=1T--il dépasse un seuil prédéfini c, le vecteur n'est pas retenu. Sinon, on extrait le vecteur --c" de -i (c'est-à-dire le vecteur de ses trois 15 composants linéairement indépendants) et on calcule la distance d2 donnée par : d2 = IIITII-di 11x11 1 (8) L'équation (8) indique qu'en raccourcissant le vecteur --c" de la valeur d2 sans modifier sa direction, ce dernier devient de norme unitaire. 20 Si c/ = -'iclî + cq est inférieur au seuil c, le vecteur t fait donc partie de l'ensemble E (du sous-ensemble El si d=0, du sous-ensemble E2 autrement ; on remarquera toutefois que cette distinction n'est pas opérationnelle, tous les vecteurs de l'ensemble E sont traités de la même façon). 25 L'équation 7 permet d'associer à chaque dit vecteur une direction de localisation de la source. Ces calculs étant très lourds, il est avantageux de les effectuer une seule fois, et de stocker les résultats dans un tableau de correspondance (voir ci-après). 30 4. Identification de la direction de localisation de la source sonore. Après avoir identifié le vecteur te E qui maximise SRP(r) indiqué par tmax, il faut l'associer à une direction de localisation de la source sonore. Si tmaxE E1, cela ne pose pas de problème, car le vecteur des différences de temps interauriculaires est alors compatible avec une et une seule direction de localisation. Si tmaxe E2, on choisit la direction de localisation compatible avec le vecteur appartenant à El qui est le plus proche de tmax selon ladite métrique prédéfinie. Avantageusement, chaque vecteur de l'ensemble E est associé une fois pour toutes à une direction de localisation. Chaque vecteur de l'ensemble et la direction de localisation qui lui est associée sont stockées dans un fichier chargé dans une mémoire du processeur pour former une table de correspondance (référence LUT, de l'anglais « Look-Up Table »), construite de la façon décrite plus haut. Ainsi, la détermination de la direction de localisation de la source à partir de la valeur de tmax se fait par simple lecture de la table. Avantageusement, le fichier contenant la table de correspondance est organisé de la façon suivante (en considérant des vecteurs à six composantes) : L'indice n'est pas explicitement écrit dans le fichier, il découle de l'ordre. Le champ « voisins » contient les indices des vecteurs qui ont des écarts inférieurs ou égaux à 1 (ou à une autre valeur prédéfinie) sur chacune des composantes. Ce champ sert à faciliter la mise en oeuvre d'un algorithme de gradient pour la maximisation de la SRP. Lorsque le processeur démarre, il charge en mémoire le contenu du fichier contenant la table de correspondance. À chaque fois qu'un son est détecté, il calcule les valeurs de SRP(T) pour tous les vecteurs stockés dans la table de correspondance, puis cherche le maximum de ces valeurs, identifie le vecteur i correspondant et lit dans la table de correspondance les valeurs d'azimut et d'élévation qui lui sont associées. En {v Indice 2 {0,3} ( 1 ,(),6,(),(-i,(1) {0,2} azimut, élévation Vecteur (1,7.6,6.6.0) 301 1 3 7 7 12 variante, le processeur pourrait chercher les N>1 plus grands maxima locaux de SRP(T) pour effectuer une multi-localisation (localisation simultanée de plusieurs sources) ; mais cette approche s'avère peu robuste. Le processeur PR, dont le fonctionnement a été illustré à l'aide 5 de la figure 1, peut comprendre un microprocesseur programmé de manière opportune couplé à une ou plusieurs mémoires pour stocker l'ensemble E et la table de correspondance LUT. Ce microprocesseur peut être dédié à la localisation des sources sonores ou accomplir également d'autres tâches ; à la limite, il peut s'agir du processeur unique d'un ordinateur portable ou de bureau 10 ou d'un robot très simple. De même, l'ensemble E et la table de correspondance LUT peuvent être stockés dans une unité de mémoire dédiée ou dans une mémoire centrale. Le processeur peut également comprendre des circuits logiques dédiés, programmables ou pas. La figure 2 représente la tête TRH d'un robot humanoïde 15 adapté pour la mise en oeuvre de l'invention, comprenant un réseau de quatre microphones M1 - M4 agencés sur une surface supérieure de ladite tête. Cet agencement permet d'éviter tout effet de masquage lorsque la source sonore est située à une hauteur supérieure à celle de la tête du robot. Cela convient particulièrement bien à la détection et localisation d'un interlocuteur humain de 20 la part d'un robot humanoïde présentant une taille inférieure à celle de la plupart des êtres humains, par exemple entre 50 et 150 cm. La distance entre deux microphones adjacents peut être comprise, à titre indicatif, entre 3 et 30 CM. Le procédé de l'invention a été testé en équipant de 25 microphones, conformément à la figure 2, un robot « Nao » produit par la société demanderesse. Le robot et une source sonore ont été placés dans une salle ordinaire, présentant une réverbération « normale », avec la source agencée devant le robot, à sa droite, à sa gauche, à l'avant-droite (à un angle de -45° par rapport à un axe arrière-avant) et à l'avant-gauche (à un angle de +45° par rapport à un axe arrière-avant). La source sonore était un équipement de reproduction sonore émettant des signaux vocaux. La localisation a été effectuée en considérant des fenêtres de calcul de 1024 échantillons chacune. Plusieurs essais ont été répétés, et les directions de localisation déterminées par le robot (identifiées par un azimut et un angle d'élévation) ont été regroupées pour former des histogrammes. Les figures 3A/3B correspondent à une source localisée devant le robot (azimut théorique : °). Les figures 4A/4B correspondent à une source localisée à la gauche du robot (azimut nominal : 90°). Les figures 5A/5B correspondent à une source localisée à l'avant-gauche du robot (azimut nominal : 45 °). Les figures 6A/6B correspondent à une source localisée à la droite du robot (azimut nominal : -90°). Les figures 7A/7B correspondent à une source localisée à l'avant- droite du robot (azimut nominal : -45 °). L'élévation nominale n'a pas été mesurée. En effet, dans les applications robotiques visées à titre principal, l'élévation est moins importante que l'azimut. Les figures « A » se rapportent au procédé de l'invention, tel que décrit plus haut. Les figures « B » sont données à titre de comparaison et ont été obtenues au moyen d'une méthode TDOA classique. Dans cette méthode de référence, la localisation est considérée en échec lorsque la valeur maximale de SRP, normalisée, est inférieure à un seuil prédéfini. Plus généralement, même dans une méthode selon l'invention, cette valeur maximale normalisée peut être considérée un indicateur de confiance de la localisation. La normalisation se fait par rapport à : À = >,\1147(si)W(si) sisi; où W(si) et W(si) sont les énergies des signaux issus des microphones i et j 25 blanchis. On peut remarquer que, dans le cas de l'invention (figures « A »), la grande majorité des essais conduit à une estimation satisfaisante, à quelques degrés près, de l'azimut de la source, tandis que la méthode de référence (figures « B ») présente un taux d'échec assez élevé. Il convient de 3 0 souligner que, dans les applications robotiques visées à titre principal, on ne demande pas une précision élevée (une erreur de localisation de quelques degrés est sans incidence), mais plutôt une grande robustesse couplée à une relative simplicité computationnelle.5 REVENDICATIONS1. Procédé de localisation d'une source sonore comportant les étapes suivantes : a) capter des signaux sonores issus d'une source sonore à localiser au moyen d'un ensemble d'au moins trois microphones (M1, M2, M3, M4) ; b) sélectionner au moins trois paires de microphones dudit ensemble et, pour chaque dite paire, calculer une intercorrélation généralisée des signaux sonores captés, ledit calcul étant effectué pour une pluralité de valeurs d'un retard - dit différence de temps interauriculaire - entre lesdits signaux sonores ; c) à partir desdites intercorrélations généralisées, calculer une puissance de réponse dirigée exprimée en fonction d'un vecteur des différences de temps interauriculaires pour chaque dite paire de microphones ; d) déterminer le vecteur des différences de temps interauriculaires qui maximise ladite puissance de réponse dirigée ; et e) estimer une direction de localisation de ladite source sonore en fonction du vecteur des différences de temps interauriculaires déterminé lors de ladite étape d), caractérisé en ce que : - lesdites étapes c) et d) sont mises en oeuvre en considérant une pluralité de vecteurs des différences de temps interauriculaires formant un ensemble (E) qui comprend : un premier sous-ensemble (El ) de vecteurs compatibles avec des signaux sonores provenant d'une source sonore unique à distance infinie desdits microphones ; et un deuxième sous-ensemble (E2) de vecteurs qui ne sont pas compatibles avec des signaux sonores provenant d'une source unique à distance infinie desdits microphones ; en ce que : - chaque vecteur dudit premier sous-ensemble est associé à une direction de localisation de la source sonore unique correspondante, et chaque vecteur dudit deuxième sous-ensemble est associé à la direction de localisation associée à un vecteur dudit premier sous-ensemble qui lui est plus proche selon une métrique prédéfinie ; et en ce que :- la direction estimée lors de ladite étape e) est celle associée au vecteur des différences de temps interauriculaires déterminé lors de ladite étape d). 2. Procédé selon la revendication 1 dans lequel les intercorrélations généralisées calculées lors de ladite étape b) sont du type GCC-PHAT, correspondant à l'intercorrélation de signaux préalablement blanchis. 3. Procédé selon l'une des revendications précédentes dans lequel, lors de ladite étape b), on sélectionne toutes les paires de microphones dudit ensemble. 4. Procédé selon l'une des revendications précédentes dans lequel ladite puissance de réponse dirigée est donnée par, ou proportionnelle à, la somme des intercorrélations généralisées pour lesdites paires de microphones. 5. Procédé selon l'une des revendications précédentes dans lequel ledit premier sous-ensemble est constitué par les vecteurs dudit ensemble dont les composantes satisfont la condition dite de Chasles : tij-Etjk=tik où les indices i#j#k représentent différent microphones dudit ensemble et 'Ci], tik et tik représentent les différences de temps interauriculaires pour les paires (i,j), (j,k) et (i,k) respectivement. 6. Procédé selon l'une des revendications précédentes dans lequel ledit deuxième sous-ensemble est constitué par les vecteurs dudit ensemble n'appartenant pas audit premier sous-ensemble et dont la distance à 3 0 un vecteur dudit premier sous-ensemble, déterminée selon ladite métrique prédéfinie et après normalisation, est inférieure à un seuil prédéfini. 7. Procédé selon l'une des revendications précédentes dans lequel ladite métrique prédéfinie est une métrique euclidienne. 8. Procédé selon l'une des revendications précédentes dans 5 lequel ladite étape e) est mise en oeuvre au moyen d'une table de correspondance (LUT) associant une direction de localisation à chaque valeur dudit ensemble. 9. Robot humanoïde comprenant : 10 - un ensemble d'au moins trois microphones (M1, M2, M3, M4) ; et - un processeur (PR) connecté auxdits microphones pour recevoir en entrée des signaux sonores captés par ces derniers et programmé ou configuré pour mettre en oeuvre un procédé selon l'une des revendications 15 précédentes. 10. Robot humanoïde selon la revendication 9 dans lequel lesdits microphones sont agencés sur une surface supérieure d'une tête (TRH) dudit robot. 20