Brevet FR2947931A1 — publié le 20110114

LOCALISATION DE SOURCES La présente invention se rapporte de manière générale au comptage automatique du nombre de sources présentes dans un mélange, et à la détermination de la direction d'arrivée de signaux émis par une pluralité de sources, et plus particulièrement au cas où le nombre de sources est inconnu a priori. L'invention trouve des applications, en particulier, dans la réalisation de codeurs/décodeurs audio multipistes, notamment ceux de type � MPEG Surround � pour créer, à partir d'une piste audio stéréo, une piste audio comportant plus de deux canaux, et/ou pour modifier les caractéristiques spatiales en simulant un déplacement des différentes sources ayant émis les signaux enregistrés. Elle s'applique aussi à la séparation de sources automatique pour réduire les perturbations de type bruit, écho, interférences liées aux prises de son éloignées des sources, notamment dans le cas d'audioconférences de groupe. Une autre application est la localisation dans une piste stéréo de la direction d'arrivée de certaines sources spécifiques afin de les supprimer, par exemple les sources correspondant aux voix afin de produire une piste musicale pour Karaoké à partir de la piste originale. Dans un espace E dans lequel un nombre N de sources émettent chacune un signal s;, la séparation aveugle desdites sources consiste à estimer la direction d'arrivée des signaux s; et le nombre N de sources à partir d'un jeu comprenant un nombre M d'observations. Chaque observation est obtenue à l'aide d'un capteur qui enregistre le signal parvenu jusqu'en un point de l'espace où se situe ledit capteur. Le signal ainsi enregistré résulte alors du mélange et de la propagation dans l'espace E des signaux s;, et se trouve donc affecté de différentes perturbations propres au milieu traversé comme par exemple le bruit, la réverbération, les interférences, etc. Les procédés connus de séparation de sources peuvent être classés en deux catégories principales. La première regroupe l'ensemble des méthodes pour lesquelles le nombre de sources N doit être connu a priori. Or, cette information n'est pas toujours disponible. Quand bien même cette dernière est accessible, à chaque première utilisation et ensuite à chaque changement du nombre de sources au cours du temps, il est nécessaire de paramétrer à nouveau lesdites méthodes. Cela rend donc particulièrement complexe l'utilisation de ces dernières dans le cadre de traitements automatisés. De plus, la connaissance du nombre de sources conditionne la précision du résultat obtenu. La rapidité à laquelle évolue le nombre de sources dans l'environnement, ainsi que le délai de modification dudit paramètre dans cette première catégorie de méthodes sont donc susceptibles de dégrader leur fiabilité. Une deuxième catégorie de méthodes, parmi lesquelles on peut penser à la méthode décrite dans le document � A Robust Method to Count and Locate Audio Sources in a Stereophonic Linear Anechoic Mixture � - Arberet, S.; Gribonval, R.; Bimbot, F.- Acoustics, Speech and Signal Processing, 2007. ICASSP 2007. IEEE International Conference on , vol.3, no., pp.11l-745-III-748, 15-20 April 2007�, propose de détecter, dans les observations M, des zones temps-fréquence pour lesquelles une seule source est dominante, c'est-à-dire une source dont l'énergie du signal s; correspondant est supérieure à la somme des énergies des autres signaux s;. On détermine alors les délais et les différences d'amplitude observés entre les capteurs deux à deux pour calculer les indices dits binauraux de chaque source. Les directions d'arrivée des signaux s; sont alors calculées, puis le nombre de sources est alors déduit en comptant le nombre de directions d'arrivée correspondant à une source active. Cependant, les méthodes de cette catégorie ne peuvent être utilisées que pour des environnements dits anéchoïques dans lesquelles les observations ne sont affectées d'aucune réverbération. En effet, ces méthodes reposent sur l'hypothèse qu'un indice binaural est invariant quelque soit la fréquence des signaux s;. Cette hypothèse s'avère valide dans le cas d'environnement anéchoïques mais n'est plus vérifiée dans un environnement affecté de réverbération. En outre, la réverbération rend délicate la détection des zones temps-fréquence où une seule source est dominante active. Il en découle que cette limitation exclut d'utiliser ces méthodes avec des observations réalisées dans des environnements échoïques, parmi lesquels on peut citer, pour les signaux acoustiques, les salles de réunion, les restaurants et plus généralement l'ensemble des lieux courants dans lesquels les prises de son sont réalisées. La présente invention vise à répondre à un problème additionnel. Les méthodes connues de détermination d'une direction d'arrivée d'un signal provenant d'une source unique dans l'espace se basent sur une différence de phase çp(t,w) du signal au niveau de plusieurs capteurs répartis dans l'espace. Pour cela, elles proposent de calculer une probabilité p de présence de la source en fonction du temps et de la direction d'arrivée, reposant par exemple sur la fonction de coût suivante : arg[R(t, w)] ù oYr P(r,t) = f p\\ G(w) dw où la fonction P est une loi de probabilité, R le rapport entre les signaux reçus par les capteurs et G(o)) est une fonction de pondération permettant d'accorder plus de poids à certaines fréquence en fonction de la configuration de l'espace. Or la différence de phase arg[R(t, w)] est une fonctionnelle discontinue car cette dernière est définie modulo 27, et présente généralement des sauts de 2rr en fonction de la fréquence, tandis que la fonction GOT est linéaire et continue. Il existe cependant des techniques dites de � phase unwrapping � en anglais pour rendre cette différence de phase linéaire, mais ces méthodes de correction de la phase ne sont pas performantes, notamment en présence de bruit ou de multiples sources. L'invention propose des moyens, décrits ci-après et notamment illustrés à la figure 3, pour résoudre ce problème additionnel. La présente invention vise à améliorer la situation. Un premier aspect de la présente invention propose un procédé de détermination d'un premier nombre de signaux émis par des sources et de leurs directions d'arrivée, au cours d'une plage temporelle subdivisée en trames, dans un espace, à partir de la connaissance d'au moins deux observations. Le procédé comporte les étapes suivantes : /a/ pour chaque trame, on calcule, à partir des observations, pour chaque direction d'un ensemble de directions de l'espace, une première probabilité de présence d'une des sources ; /b/ pour un deuxième nombre de trames d'une fenêtre temporelle, on calcule une deuxième probabilité de présence d'une des sources dans chaque direction d'un ensemble de directions de l'espace en fonction de la première probabilité; /c/ on recherche les directions pour lesquelles il existe un maximum local de la deuxième probabilité, lesdites directions correspondant chacune à la direction d'arrivée d'un des signaux. L'invention propose ainsi une méthode apte à traiter une pluralité de sources. Le procédé est fonctionnel, y compris lorsque le nombre de sources est supérieur au nombre d'observations. Ce procédé peut être automatisé. En particulier, il n'est pas nécessaire a priori d'indiquer au procédé le nombre de sources. On peut prévoir qu'au cours de l'étape /c/, seules les directions correspondant à un maximum local supérieur à un seuil sont considérées comme des directions d'arrivée d'un des signaux. Ce filtrage par seuil permet en outre d'améliorer la robustesse du procédé face aux fausses détections de direction d'arrivée, alors qu'aucune source ne correspond. Dans le cas par exemple où ce paramètre est inconnu, le procédé peut comporter une étape /e/ où l'on détermine le premier nombre de sources en dénombrant le nombre de maxima locaux de la deuxième probabilité. Ainsi, il est possible de fournir cette information en sortie du procédé par le calcul, sans connaissance a priori de ce nombre. Dans un mode de réalisation de la présente invention, la première probabilité est calculée en réalisant les étapes suivantes : /d/ on détermine les différences de phase entre les observations ; /e/ on obtient la première probabilité (p), en calculant une fonction de coût (p) basée sur le niveau de corrélation dans le domaine vectoriel complexe du rapport (R) entre les observations (x,, x2) et un modèle théorique de ce rapport . En particulier, la première probabilité peut être calculée, en définissant une fonction de pondération de plages de fréquences (G(w)) et à l'aide de l'expression mathématique suivante : p(r;t) =1P 'arg(R(t w)e 1°`)" G(w) dw La première probabilité ainsi obtenue ne subit ainsi aucun saut de 10 fréquence, puisque le niveau de corrélation est calculé dans le domaine complexe. Dans un mode de réalisation de la présente invention, on calcule pour chaque direction de l'espace la deuxième probabilité en déterminant le maximum de la première probabilité sur la fenêtre temporelle. Alternativement, 15 on peut calculer la deuxième probabilité en : • déterminant, pour chaque trame temporelle, la direction d'arrivée dominante z(t)=argmaxp(r,t) ; • calculant une probabilité de dominance m(t)=max p(r,t) • appliquant, pour obtenir la deuxième probabilité, et pour une 20 valeur (0) définissant un degré de lissage, l'expression mathématique suivante : p(T)=>m(t)I Z-e 25 sont susceptibles de se déplacer dans l'espace. Cela permet notamment d'adapter le procédé selon l'invention aux caractéristiques des sources. Un deuxième aspect de la présente invention propose un dispositif comprenant des moyens adaptés pour mettre en oeuvre un procédé de détermination de la direction d'arrivée selon le premier aspect de la présente invention. Un troisième aspect de la présente invention propose un dispositif de décodage audio comprenant un dispositif d'estimation de la direction d'arrivée selon le deuxième aspect de la présente invention. A titre d'exemple, le dispositif de décodage peut décoder des flux stéréo sans information auxiliaire en flux à la norme � MPEG surround � mixés en sur 5 + 1 canaux. Un quatrième aspect de la présente invention propose un programme d'ordinateur comportant des instructions pour la mise en oeuvre du procédé selon le premier aspect de la présente invention lorsque ce programme est exécuté par un processeur. D'autres aspects, buts et avantages de l'invention apparaîtront à la lecture de la description d'un de ses modes de réalisation. L'invention sera également mieux comprise à l'aide des dessins, sur lesquels - la figure 1 illustre un espace comportant une pluralité de sources ; - la figure 2 illustre les étapes principales d'un procédé de détermination de la direction d'arrivée des signaux dans l'espace à partir des observations, mis en oeuvre selon un mode de réalisation de la présente invention ; - la figure 3 illustre les étapes principales d'un mode de réalisation du calcul de la première probabilité p ; - la figure 4 montre, par un schéma de principe, un dispositif d'estimation de la direction d'arrivée de signaux, selon un mode de réalisation de la présente invention ; - la figure 5 illustre un dispositif de décodage, selon un mode de réalisation de la présente invention.. Dans la présente description, et comme illustré sur la figure 1, on considère trois sources 10a, 10b, 10c, émettant respectivement un signal s,, S2, s3, dans l'espace E. Les sources peuvent se déplacer au cours du temps. Par rapport à un point A de l'espace E, et à un instant donné, le signal si a une direction d'arrivée e1, le signal s2 une direction d'arrivée e2, le signal s3 une direction d'arrivée e3 . Sur une plage temporelle P, subdivisée en trames T, des observations xi(t) et x2(t) ont été réalisées respectivement en un point 01 et un point 02 de l'espace E. La présente description fait état, à titre illustratif uniquement, d'un nombre N de sources égal à 3 et d'un nombre M d'observation égal à 2. Il est aisé de l'appliquer à toute autre combinaison de nombre N de sources et de nombre M d'observations, M étant supérieur ou égal à 2. Dans un espace E échoïque, les observations xi(t) et x2(t) peuvent être modélisées par un mélange convolutif bruité composé des signaux si, s2, s3. En posant {afi(k)} les coefficients de la réponse impulsionelle du filtre séparant la ?me source du je' capteur, bi(t) le bruit additif diffus indépendant des sources, * étant le symbole de la convolution, les observations xi(t) et x2(t) peuvent être modélisées ainsi : N N (t)=~~a~i(k)s.(tùk)+nj(t)=(ai,.*s,.)(t)+bi(t) j=1,...,M. Dans le domaine fréquentiel, par transformation temps-fréquence, en posant u(k) une fenêtre apodisante de type fenêtre de Nanning par exemple, on obtient alors: L N Xj(t,û)=u(k)x.,(t+k)e °=A,(û)si(t,w)+BB(û) k=0 où Afi(w) est la transformée temps-fréquence de afi(k), Bj(w) la transformée 20 temps-fréquence de b;(t), et Sj(t,w) la transformée temps-fréquence à court-terme de s;(t). Les signaux si, s2, s3 ont pour caractéristique qu'il existe pour chacun d'eux au moins une trame tp de la plage temporelle P au cours de laquelle l'énergie dudit signal est supérieur à la somme de l'énergie des autres 25 signaux. La source à l'origine dudit signal est alors dite dominante au cours de cette trame tp. L'expression mathématique suivante traduit cette caractéristique: t,, ={tipP i s,(t,w) o s,,(t,w)} avec card(tp)>-1. Il résulte de cette caractéristique, que pour l'ensemble des trames tp, les observations X1(tp, w) et X2(tp, w) dans le domaine temps-fréquence peuvent 30 être approchées par l'expression mathématique suivante : XX(tn,o)~A~n(co)sn(tl,,w)+BJ(co) où B'i est la somme de Bj et des résidus des autres sources non dominantes. La figure 2 illustre les étapes principales d'un procédé de détermination de la direction d'arrivée des signaux dans l'espace E à partir des observations, mis en oeuvre selon un mode de réalisation de la présente invention. Dans une première étape 50, à partir des M observations, on calcule pour toutes les trames T, une première probabilité p de présence d'une des sources pour chaque direction d'arrivée. On peut, de manière équivalente, calculer ladite probabilité p pour une direction exprimée à l'aide d'un vecteur ou pour une direction exprimée à l'aide d'un angle. De même, il est possible de limiter, en fonction de l'espace E, l'intervalle pour lequel on va chercher à obtenir la première probabilité p. Par exemple, il peut être souhaitable de limiter les calculs aux directions comprises dans un cône donné de l'espace. En outre, il est possible de définir un pas de résolution spatiale avec lequel on va balayer l'espace E, pour limiter le nombre de première probabilité p effectivement calculées. Dans une deuxième étape 60, à partir des premières probabilités p, on calcule une deuxième probabilité p de présence d'une des sources en fonction de la direction d'arrivée, pour un sous-ensemble composé d'un nombre g de trames T formant une fenêtre temporelle F. Le nombre g peut être compris entre 1 et le nombre de trames T compris dans la plage temporelle P. Dans le cas où les sources sont susceptibles de se déplacer dans l'espace E, le choix du nombre g de trames T de la fenêtre F est fonction de la vitesse de déplacement des sources. De manière générale, plus la vitesse des sources est susceptible d'être importante, plus le nombre g sera petit. En conséquence, le deuxième nombre g de trames T de la fenêtre temporelle F peut être choisi de sorte à être inversement proportionnel à la vitesse à laquelle les sources sont susceptibles de se déplacer dans l'espace E. Si les sources se déplacent à des vitesses différentes, on pourra par exemple considérer soit la vitesse minimale, soit la vitesse maximale ou encore la vitesse moyenne. On recherche, dans une troisième étape 70, les directions pour lesquelles il existe un maximum local, qualifié dans ce contexte de pic, de la deuxième probabilité �, lesdites directions correspondant chacune à la direction d'arrivée d'un des signaux. A titre d'exemple non limitatif, le maximum local peut être obtenu en lissant, dans un premier temps, la deuxième probabilité � à l'aide par exemple d'un filtre passe-bas de premier ordre, puis en recherchant les directions pour lesquelles la dérivée première de la deuxième probabilité � donne une valeur nulle et pour lesquelles la dérivée seconde de la deuxième probabilité � donne une valeur négative. On peut alors encore éventuellement calculer l'énergie de la deuxième probabilité � pour une plage de valeur autour desdits maxima locaux et éliminer les maxima locaux pour lesquels cette énergie est inférieure à un seuil prédéterminé. Dans un mode de réalisation de la présente invention, le procédé comporte une étape optionnelle 80 au cours de laquelle on détermine le nombre N de sources en dénombrant le nombre de maxima locaux de la deuxième probabilité p. La figure 3 illustre les étapes principales d'un mode de réalisation du calcul de la première probabilité p, mis en oeuvre selon un mode de réalisation de la présente invention. Pour calculer la première probabilité p à l'étape 50, pour une trame T donnée, on détermine, dans une étape 52, la différence de phase ço(t,w) entre les observations X1 et X2 correspondantes. La différence de phase çp(t,w) est supposée linéaire, ce qui est vérifiée y compris dans un espace E échoïque. En effet, la différence de phase théorique q(tp,w) entre les observations X1 et X2 correspondantes, pour les trames tp où une des sources est dominante. A ces instants, la différence de phase q(tp,w) entre les points 01 ou 02 peut être modélisée par l'expression mathématique suivante : X(t,w) çp(tp,w)=argR(tp,w)=arg P =r w+h(tp,w) XZ (tp w) où R est le rapport entre les observations X1 et X2, et h est la somme totale du bruit causé par le bruit additif, le résidus des sources non dominantes et la réverbération. La différence de phase q(tp,w) est donc linéaire en fonction de la fréquence w, tant que h(tp,w) ne dégrade pas la linéarité de çp(tp,w). Pour une source donnée émettant un signal Si dont la direction d'arrivée est notée 8;, ledit signal Si arrive à l'un des deux points 01 ou 02 avec un décalage temporel z par rapport à l'autre point 01 ou 02. Le décalage temporelle z peut être estimé à dcos(9) avec d la distance entre les points 01 c ou 02 et c la célérité du signal Si dans l'espace E. Il en découle que la détermination de la différence de phase ço(t,w) permet l'obtention de la direction d'arrivée 8; pour chaque source. Dans une étape 54, on obtient la première probabilité p en calculant le niveau de corrélation dans le domaine vectoriel complexe d'une fonction de coût portant sur le rapport R entre les observations X1 et X2. La probabilité p peut être obtenue en appliquant la formule mathématique suivante : arg(R(t,w)é-'�' p(r;t) = f p G(w) dw Dans un premier mode de réalisation, la deuxième probabilité �, pour une direction donnée, calculée au cours de la deuxième étape 60, peut quant elle être obtenue en identifiant le maximum de la première probabilité p sur la 20 fenêtre temporelle F, ce qui peut se traduire par l'expression mathématique suivante : ,u(r) = max p(r,t) . tEF Dans un deuxième mode de réalisation, la deuxième probabilité �, pour une direction donnée, calculée au cours de la deuxième étape 60, peut quant elle être obtenue en déterminant pour chaque trame temporelle la direction 25 d'arrivée dominante égale à la position du maximum sur les directions étudiées î(t) = argmax p(r,t) et en calculant la probabilité de cette dominance égale à la r valeur de ce maximum m(t)=maxp(r,t). Ainsi, la deuxième probabilité est donnée par l'histogramme pondéré calculé à partir de l'ensemble des directions d'arrivées dominantes et leur où la fonction p est une loi de probabilité, par exemple une loi de Poisson ou de Gauss, et G(o) est une fonction de pondération permettant d'accorder plus 15 de poids à certaines fréquence en fonction de la configuration de l'espace. probabilité de dominance. u(r)=m(t) i r-e La figure 4 montre, sur un schéma de principe, un dispositif 100 d'estimation de la direction d'arrivée de signaux, selon un mode de réalisation de la présente invention. Le dispositif 100 est notamment adapté à la mise en oeuvre du procédé selon l'invention. Il comporte une entrée 100 lui permettant de recevoir les observations x1, x2 des signaux S1, S2, S3, sur la plage temporelle P. Une unité de transformation Temps-Fréquence 106, par exemple une unité adaptée à la mise en oeuvre d'une transformée de Fourrier rapide communément appelée � FFT �, permet par la suite de travailler sur les observations dans le domaine fréquentiel, les observations notées x1, x2 dans le domaine temporel étant notées classiquement X1, X2 dans le domaine fréquentiel. II comporte une unité de calcul de la direction d'arrivée 110. Cette dernière est reliée à l'unité de transformation Temps-Fréquence 106. Elle est adaptée au calcul de la première probabilité p de présence d'une des sources S pour chaque direction de l'espace E à partir des observations X1, X2. Le dispositif 100 comprend une unité de regroupement temporel 125, coopérant avec l'unité de calcul de la direction d'arrivée 110. Cette unité de regroupement 125 est adaptée au calcul de la deuxième probabilité p de présence d'une des sources pour chaque direction de l'espace, en fonction de la première probabilité p, et sur la fenêtre temporelle F. Le dispositif 100 comporte une unité d'identification 130, coopérant avec l'unité de regroupement temporel 125, adaptée à identifier des directions pour lesquelles il existe un maximum local de la deuxième probabilité �. L'unité d'identification 130 est reliée à la sortie 140 du dispositif 100 de sorte à pouvoir délivrer les directions identifiées correspondant aux directions d'arrivée 9Z des signaux Si. Le dispositif peut également comporter des moyens de comptage 135 pour délivrer sur la sortie 140 le premier nombre N de sources en dénombrant le nombre de maxima locaux de la deuxième probabilité �. Le dispositif peut également comporter des moyens de paramétrage 120 adaptés à modifier, au niveau de l'unité de regroupement temporel 125, le deuxième nombre g de trames T de manière inversement proportionnelle à la vitesse à laquelle les sources sont susceptibles de se déplacer dans l'espace E. La figure 5, illustre par un synoptique, un dispositif de décodage audio, selon un mode de réalisation de la présente invention. Un tel dispositif est par exemple conçu pour notamment créer des flux à la norme � MPEG surround � mixés en 5+1 canaux à partir d'un flux stéréo sans information auxiliaire. Le dispositif de décodage 210 reçoit en entrée des observations x1,..., xN, typiquement un signal stéréo issu du codeur de type AAC par exemple et contenant des signaux Si émis par une pluralité de sources. Le dispositif de décodage comporte un dispositif 100 d'estimation de la direction d'arrivée de signaux selon l'invention, recevant lui aussi les observations x1,..., xN. A partir des informations fournies par le dispositif 100, le dispositif de décodage audio comporte les moyens de traitement 215 nécessaires pour générer des flux multiples spatialisés sur une sortie 220 à partir des directions d'arrivée des signaux et éventuellement du nombre de sources. REVENDICATIONS1. Procédé de détermination d'un premier nombre (N) de signaux (Si) émis par des sources (10) et de leurs directions d'arrivée (9Z ), au cours d'une plage temporelle (P) subdivisée en trames (T), dans un espace (E), à partir de la connaissance d'au moins deux observations (xi, x2), caractérisé en ce qu'il comporte les étapes suivantes : /a/ pour chaque trame (T), on calcule (50), à partir des observations (xi, x2), pour chaque direction d'un ensemble de directions de l'espace (E), une première probabilité (p) de présence d'une des sources ; /b/ pour un deuxième nombre (g) de trames (T) d'une fenêtre temporelle (F), on calcule (60) une deuxième probabilité (�) de présence d'une des sources (10) dans chaque direction d'un ensemble de directions de l'espace (E) en fonction de la première probabilité (p) ; /c/ on recherche (70) les directions pour lesquelles il existe un maximum local de la deuxième probabilité (�), lesdites directions correspondant chacune à la direction d'arrivée (9z) d'un des signaux (Si) ; 2. Procédé selon la revendication 1, dans lequel au cours de l'étape /c/, seules les directions correspondant à un maximum local supérieur à un seuil sont considérées comme des directions d'arrivée (9z) d'un des signaux (Si). 3. Procédé selon l'une quelconque des revendications précédentes comportant en outre une étape /e/ où l'on détermine (80) le premier nombre (N) de sources (10) en dénombrant le nombre de maxima locaux de la deuxième probabilité (�). 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel la première probabilité (p) est calculée (50) en réalisant les étapes suivantes : /d/ on détermine (52) les différences de phase (ço(t,w)) entre les observations (xi, x2) ; /e/ on obtient (54) la première probabilité (p), en calculant une fonction de coût (p) basée sur le niveau de corrélation dans le domaine vectoriel complexe du rapport (R) entre les observations (xi, x2) et un modèle théorique de ce rapport. 5. Procédé selon la revendication 4, dans lequel la première probabilité (p) est calculée, en définissant une fonction de pondération de plages de fréquences (G(w)) et à l'aide de l'expression mathématique suivante : arg(R(t,w)é-''r p(r;t)=$p G(w) dw 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel on calcule (60) pour chaque direction de l'espace (E) la deuxième probabilité (�) en déterminant le maximum de la première probabilité (p) sur la 20 fenêtre temporelle (F). 7. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel on calcule (60) la deuxième probabilité (Il) en : • déterminant, pour chaque trame temporelle, la direction d'arrivée 25 dominante 2(t)= argmax p(r,t) • calculant une probabilité de dominance m(t)=maxp(r,t) ; • appliquant, pour obtenir la deuxième probabilité, et pour une valeur (c) définissant un degré de lissage, l'expression mathématique suivante : 30 u(r)=Dm(t) z-e 8. Dispositif (100) d'estimation de la direction d'arrivée (9Z ), au cours d'une plage temporelle (P) subdivisée en trames (T), dans un espace (E), d'un premier nombre (N) de signaux (Si) émis respectivement par un même premier nombre (N) de sources (10), ledit dispositif étant adapté à la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce qu'il comporte: • une entrée (105) pour recevoir au moins deux observations (xi, x2) desdits signaux (Si) sur la plage temporelle (P) ; • une unité de calcul de la direction d'arrivée (110), apte à recevoir les observations (xi, x2), ladite unité (110) étant adaptée au calcul d'une première probabilité (p) de présence d'une des sources (S) pour chaque direction d'un ensemble de directions de l'espace (E); • une unité de regroupement temporel (125), coopérant avec l'unité de calcul de la direction d'arrivée (110), adaptée au calcul d'une deuxième probabilité (�) de présence d'une des sources (10) pour chaque direction d'un ensemble de directions de l'espace (E), en fonction de la première probabilité (p), et sur une fenêtre temporelle (F) composée d'un deuxième nombre (g) de trames (T); • une unité d'identification (130), coopérant avec l'unité de regroupement temporel (125), adaptée à identifier des directions pour lesquelles il existe un maximum local de la deuxième probabilité (�), lesdites directions correspondant chacune à la direction d'arrivée (9z) d'un des signaux (Si) d'un des signaux (Si) et étant délivrées sur une sortie (140). 9. Dispositif selon la revendication 8, comportant en outre des moyens de comptage (135) pour délivrer sur la sortie (140) le premier nombre (N) de sources (10) en dénombrant le nombre de maxima locaux de la deuxième probabilité (�). 10. Dispositif selon l'une quelconque des revendications 8 à 9, comportant des moyens de paramétrage (120) adaptés à modifier, au niveau de l'unité de30regroupement temporel (125), le deuxième nombre (g) de trames (T) de manière inversement proportionnelle à la vitesse à laquelle les sources (10) sont susceptibles de se déplacer dans l'espace (E). 11. Dispositif de décodage audio comportant une entrée pour recevoir au moins deux observations (xi, x2) de signaux (Si) émis par des sources (10), ledit dispositif comportant un dispositif d'estimation de la direction d'arrivée (9z) selon l'une quelconque des revendications 8 à 10. 12. Programme d'ordinateur comportant des instructions pour la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 7 lorsque ce programme est exécuté par un processeur.