FR 2948484 A1 20110128 FR 0955133 A 20090723 L'invention concerne le traitement de la parole en milieu bruité. Elle concerne notamment, mais de façon non limitative, le traitement des signaux de parole captés par des dispositifs de téléphonie pour véhicules automobiles. Ces appareils comportent un microphone ("micro") sensible captant non seulement la voix de l'utilisateur, mais également le bruit environnant, bruit qui constitue un élément perturbateur pouvant aller, dans certains cas, jusqu'à rendre incompréhensibles les paroles du locuteur. Il en est de même si l'on veut mettre en oeuvre des techniques de reconnaissance vo- cale, car il est très difficile d'opérer une reconnaissance de forme sur des mots noyés dans un niveau de bruit élevé. Cette difficulté liée aux bruits environnants est particulièrement contraignante dans le cas des dispositifs "mains-libres". En particulier, la dis-tance importante entre le micro et le locuteur entraîne un niveau relatif de bruit élevé qui rend difficile l'extraction du signal utile noyé dans le bruit. De plus, le milieu très bruité typique de l'environnement automobile pré-sente des caractéristiques spectrales non stationnaires, c'est-à-dire qui évoluent de manière imprévisible en fonction des conditions de conduite : passage sur des chaussées déformées ou pavées, autoradio en fonction- nement, etc. Certains de ces dispositifs prévoient l'utilisation de plusieurs micros, généralement deux micros, et utilisent la moyenne des signaux captés, ou d'autres opérations plus complexes, pour obtenir un signal avec un niveau de perturbations moindre. En particulier, une technique dite beamforming permet de créer par des moyens logiciels une directivité qui améliore le rapport signal/bruit, mais les performances de cette technique sont très limitées lorsque seulement deux microphones sont utilisés. Par ailleurs, les techniques classiques sont surtout adaptées au filtrage des bruits diffus, stationnaires, provenant des alentours du dispositif et se retrouvant à des niveaux comparables dans les signaux captés par les deux micros. En revanche, un bruit non stationnaire, c'est-à-dire évoluant de manière imprévisible en fonction du temps, ne sera pas discriminé de la parole et ne sera donc pas atténué. Or, dans un environnement automobile ces bruits non stationnaires et directifs sont très fréquents : coup de klaxon, passage d'un scooter, dépassement par une voiture, etc. L'une des difficultés du filtrage de ces bruits non stationnaires tient au fait que leurs caractéristiques temporelles et spatiales sont très proches de celles de la parole, d'où la difficulté d'une part, d'estimer la présence d'une parole (car le locuteur ne parle pas tout le temps) et d'autre part d'extraire le signal utile de parole dans un environnement très bruité tel qu'un habitacle de véhicule automobile. L'un des buts de l'invention est de mettre à profit la structure multimicrophone du dispositif pour opérer une détection spatiale de ces bruits non stationnaires, puis de discriminer, parmi toutes les composantes non stationnaires (ci-après "transients") celles qui sont des composantes de bruit non stationnaires d'avec celles qui sont des composantes de parole, et enfin de traiter le signal capté pour le débruiter de manière efficace tout en minimisant les distorsions introduites par ce traitement. Dans la suite, on appellera "bruit latéral" un bruit non stationnaire directif dont la direction d'arrivée est éloignée de celle du signal utile, et on appellera "cône privilégié" la direction ou secteur angulaire de l'espace où se trouve la source de signal utile (la parole du locuteur) par rapport au ré-seau de micros. Lorsqu'une source sonore se manifestera en dehors du cône privilégié, il s'agira donc d'un bruit latéral, que l'on cherchera à atténuer. Le point de départ de l'invention consiste à associer les propriétés de non- stationnarité temporelle et fréquentielle, d'une part, et de directivité spatiale, d'autre part, pour détecter un type de bruit qu'il est d'ordinaire difficile de discriminer de la parole, puis pour en déduire une probabilité de présence de la parole qui servira à atténuer ce bruit. Plus précisément, l'invention a pour objet un procédé de débruitage d'un signal acoustique bruité capté par une pluralité de microphones d'un dis-positif audio multi-microphone opérant dans un milieu bruité. Le signal acoustique bruité comprend une composante utile de parole issue d'une source de parole directive et une composante parasite de bruit, cette composante de bruit incluant elle-même une composante de bruit latéral non stationnaire directif. Essentiellement, et de façon caractéristique de l'invention, le procédé comporte les étapes suivantes de traitement, exécutées dans le domaine fréquentiel: a) combinaison de la pluralité de signaux captés par la pluralité corres-5 pondante de microphones en un signal combiné bruité ; b) à partir du signal combiné bruité, estimation d'une composante de bruit pseudo-stationnaire contenue dans ce signal combiné bruité ; c) à partir de la composante de bruit pseudo-stationnaire estimée à l'étape b) et du signal combiné bruité, calcul d'une probabilité de pré- 10 sence de transients dans le signal combiné bruité ; d) à partir de la pluralité de signaux captés par la pluralité correspondante de microphones et de la probabilité de présence de transients calculée à l'étape c), estimation d'une direction principale d'arrivée des transients ; 15 e) à partir de la direction principale d'arrivée des transients estimée à l'étape d), calcul d'une probabilité de présence de parole sur un critère spatial, propre à discriminer entre parole utile et bruit latéral parmi les transients ; f) à partir de la probabilité de présence de parole calculée à l'étape e) 20 et du signal combiné bruité, réduction sélective du bruit par application d'un gain variable propre à chaque bande de fréquences et à chaque trame temporelle. Selon diverses formes de mise en oeuvre subsidiaires avantageuses : le traitement de l'étape a) est un traitement de préfiltrage de type fixed 25 beamforming ; le traitement de l'étape d) comprend les sous-étapes successives sui-vantes : dl) partition de l'espace en une pluralité de secteurs angulaires ; d2) pour chaque secteur, évaluation d'un estimateur de direction d'arrivée à partir de la pluralité de signaux captés par la pluralité cor- 30 respondante de microphones ; d3) pondération de chaque estimateur par la probabilité de présence de transients calculée à l'étape c) ; d4) à partir des valeurs d'estimateurs pondérées calculées à l'étape d3), estimation d'une direction principale d'arrivée des transients ; et d5) validation ou invalidation de l'estimation de la direction principale d'arrivée 35 des transients opérée à l'étape d4). - à l'étape d5) l'estimation n'est validée que si la valeur de l'estimateur pondéré correspondant à la direction estimée est supérieure à un seuil prédéterminé, et/ou en l'absence de maximum local de l'estimateur pondéré dans le secteur angulaire d'origine du signal de parole utile, et/ou que si la valeur de l'estimateur est croissante de façon monotone sur une pluralité de trames temporelles successives ; - le procédé comprend en outre une étape de maintien de l'estimation de la direction principale d'arrivée pendant un laps de temps minimal pré-déterminé ; - la probabilité de présence de parole calculée à l'étape e) est soit une probabilité binaire, prenant une valeur 1 ou 0 selon que la direction principale d'arrivée des transients estimée à l'étape d) est située ou non dans le secteur angulaire d'origine du signal de parole utile, soit une probabilité à valeurs multiples, fonction de l'écart angulaire entre la di- rection principale d'arrivée des transients estimée à l'étape d) et la di- rection d'origine du signal de parole utile ; le traitement de l'étape f) est un traitement de réduction sélective du bruit par application d'un gain à amplitude log-spectrale modifié optimisé OM-LSA. 0 On va maintenant décrire un exemple de mise en oeuvre du procédé de l'invention en référence à la figure annexée. La Figure 1 est un schéma par blocs montrant les différents modules et 25 fonctions mis en oeuvre par le procédé de l'invention ainsi que leurs inte- ractions. 0 Le procédé de l'invention est mis en oeuvre par des moyens logiciels, qu'il 30 est possible de décomposer et schématiser par un certain nombre de mo- dules 10 à 24 illustrés Figure 1. Ces traitements sont mis en oeuvre sous forme d'algorithmes appropriés exécutés par un microcontrôleur ou un processeur numérique de signal. Bien que, pour la clarté de l'exposé, ces divers traitements soient présen- 35 tés sous forme de modules distincts, ils mettent en oeuvre des éléments communs et correspondent en pratique à une pluralité de fonctions globalement exécutées par un même logiciel. Le signal que l'on souhaite débruiter est issu d'une pluralité de signaux captés par un réseau de micros (qui, dans la configuration minimale, peut être simplement un réseau de deux micros) disposés selon une configuration prédéterminée. Le réseau de micros capte le signal émis par la source de signal utile (signal de parole), et la différence de position entre les micros induit un en-semble de déphasages et variations d'amplitude dans l'enregistrement des signaux émis par la source de signal utile. Plus précisément, le micro d'indice n délivre un signal : xä(t) = x s(tùzä)+vä(t) où aä est l'atténuation d'amplitude due à la perte d'énergie entre la position de la source sonore s et le micro, 'r,-, est le déphasage entre le signal émis et reçu par le micro et vä représente la valeur du champ de bruit diffus à la position du micro. Dans la mesure où la source est éloignée d'au moins quelques centimè- tres des micros, on pourra faire l'approximation que la source sonore émet une onde plane. Les retards rn pourront alors être calculés à partir de l'anglets, défini comme l'angle entre les médiatrices des couples de micros (n, m) et la direction de référence correspondant à la source s de signal utile. Lorsque le système considéré comporte deux micros dont la médiatrice coupe la source, l'angle Os est nul. Transformée de Fourier des signaux captés par les micros (blocs 10) Le signal dans le domaine temporel x,,(t) issu de chacun des N micros est numérisé, découpé en trames de T points temporels, fenêtré temporelle-ment par une fenêtre de type Nanning, puis la transformée de Fourier rapide FFT (transformée à court terme) X,,(k,l) est calculée pour chacun de ces signaux : X,, (k, l) = a,,.d,, (k) x S(k, l) + Vä (k, l ) avec : 1 étant l'indice de la trame temporelle, k étant l'indice de la bande de fréquences, et fk étant la fréquence centrale de la bande de fréquence indicée par k. Constitution d'un signal combiné partiellement débruité (bloc 12) Les signaux X,, (k,l) peuvent être combinés entre eux par une technique simple de préfiltrage par beamforming du type Delay and Sum qui est appliquée pour obtenir un signal combiné X(k,l) partiellement débruité : X(k,l) = N d,(k)Xn(k,l) nE[1,N] Il est à noter que, concrètement, le nombre de micros étant limité, ce traitement ne procure qu'une faible amélioration du rapport signal/bruit, de l'ordre de 1 dB seulement. Lorsque le système considéré comporte deux micros dont la médiatrice 20 coupe la source, l'angle 0s est nul et il s'agit d'une simple moyenne qui est faite sur les deux microphones. Estimation du bruit pseudo-stationnaire (bloc 14) 25 Cette étape a pour objet de calculer une estimation de la composante de bruit pseudo-stationnaire V(k,l) présente sur le signal X(k,l). Il existe de très nombreuses publications sur ce sujet, l'estimation et la réduction du bruit pseudo-stationnaire étant en effet un problème classique assez bien résolu. Différentes méthodes sont efficaces et utilisables pour 30 obtenir V(k,l), notamment un algorithme d'estimation de l'énergie de la composante de bruit pseudo-stationnaire à moyennage récursif par contrôle des minima (MCRA) comme celui décrit par I. Cohen et B. Berdugo, Noise Estimation by Minima Controlled Recursive Averaging for 6 5 Robust Speech Enhancement, IEEE Signal Processing Letters, Vol. 9, No 1, pp. 12-15, Jan. 2002. Calcul de la probabilité de présence des transients (bloc 16) Les "transients" désignent tous les signaux non-stationnaires, incluant aussi bien la parole utile que les bruits non-stationnaires sporadiques, qui peuvent avoir une énergie équivalente ou parfois supérieure à la parole utile (passage d'un véhicule, sirène, klaxon, parole d'autres personnes 10 etc. ). II est possible de détecter ces transients à l'aide de l'estimation précédemment établie de la composante de bruit pseudo-stationnaire en retranchant cette dernière du signal global X(k,l). On verra plus loin (description détaillée des blocs 18 et 20) la manière 15 dont il est possible de discriminer parmi ces transients entre ceux qui correspondent à la parole utile et ceux qui correspondent à des bruits non-stationnaires et qui ont des caractéristiques similaires à la parole utile. Le traitement opéré par le bloc 16 consiste seulement à calculer une probabilité priansient (k, l) de présence de signaux transients, sans distinction entre 20 parole utile et bruits parasites non-stationnaires. L'algorithme est le suivant : Pour chaque trame 1 et pour chaque bande de fréquence k, Calculer le "Transient to Stationary Ratio" : TSR(k,l)= X(k,1)-V(k,l) V(k,l) (ii) Si TSR(k, l) TSRmax PTransient (k, l) = 1 30 (iv) Si TSR,nin 10 stationnaire du bruit qui est diffuse, les transients sont souvent directifs, c'est-à-dire issus d'une source sonore ponctuelle (comme la bouche du locuteur pour la parole utile, ou le moteur d'une motocyclette pour un bruit latéral). Il est donc judicieux de calculer la direction d'arrivée de ces signaux, qui sera en général bien définie, et de comparer cette direction 15 d'arrivée à l'angle Os correspondant à la direction d'origine parole utile), de manière à déterminer si le signal non-stationnaire considéré est utile ou parasite, et d'effectuer ainsi la discrimination entre parole utile et bruit non-stationnaire. La première étape consiste à estimer la direction d'arrivée du transient. 20 La méthode utilisée ici est basée sur l'utilisation de la probabilité de présence des transients p,.,.�äsieäi(k,l) déterminée par le bloc 18 de la manière exposée plus haut. Plus précisément, on opère une partition de l'espace en secteurs angulaires, chacun correspondant à une direction définie par un angle 8,i e [1,M] 25 (par exemple M=19, avec la collection d'angles {-9o°,-8o°...,o°,...+8o°,+90°}). On notera qu'il n'y a aucun lien entre le nombre N de micros et le nombre Md'angles testés. Par exemple, il est tout à fait possible de tester une dizaine d'angles (M=10) avec un seul couple de micros (N=2). Chaque angle Bi est testé de façon à déterminer celui qui est le plus pro- 30 che de la direction d'arrivée du signal non-stationnaire étudié. Pour ce faire, on considère chaque couple de micros (n, m) et on calcule un estimateur de direction d'arrivée P,,m(Bi,k,l) correspondant, dont le module sera maximal lorsque l'angle Bi testé sera le plus proche de la direction d'arrivée du transient. Cet estimateur peut par exemple s'appuyer sur un calcul d'intercorrélation et prendre la forme : Pä ä, (O , k, 1) = E(X,n (k, l)X (k, l ).e ;2 n l , avec ri = l"'' sin 9 l c 1 étant la distance entre les micros d'indices n et m , et c étant la célérité du son. Une première méthode, classique, consiste à prendre pour estimation de la direction d'arrivée l'angle qui maximise le module de cet estimateur, 10 soit : esta (k, l) = arg max O,,ie[1,M](ei,k,l) Une autre méthode, utilisée ici de façon préférentielle, consiste à pondérer l'estimateur P, ,n, (ei, k, l) par la probabilité de présence de transients 15 prransien(k,l), et définir une nouvelle stratégie de décision. L'estimateur de direction d'arrivée correspondant sera : PNev,,,m (o , k, l) Pn (8;,k,l)x PTYa"sient(k,l 20 L'estimateur peut être moyenné sur les couples de micros (n,m) 1 PNew(ei,k,1) = N(N-1) n PNew,,.,,,(8e,k,l) L'intégration de la probabilité de présence de transients dans l'estimateur 25 de direction d'arrivée présente trois avantages importants : - l'estimation de direction est ciblée sur les parties non-stationnaires du signal (où la probabilité pT,.aäsient(k,l)est proche de 1), dont la direction d'arrivée est bien définie, ce qui rend l'estimation consistante ; - l'estimation de direction est robuste au bruit diffus (où la probabilité 30 pr,.ansient(k,l)est proche de zéro), qui d'ordinaire perturbe les estima- tions de direction d'arrivée ; - la fiabilité de l'estimateur PNew,,.,,,(et, k, l) permet de distinguer plu-sieurs signaux non-stationnaires correspondant à différentes directions et simultanément présents (on verra plus bas que cette distinction peut se faire par bande de fréquences ou par analyse des maxi- ma angulaires locaux sur une même bande de fréquences). Ainsi, si l'on a en même temps un signal de parole utile et un bruit latéral puissant, les deux types de signaux seront détectés, évitant que le signal de parole utile concomitant soit éliminé par erreur dans la suite du processus, même si son énergie est faible. On va maintenant expliciter les règles de décision permettant à partir de PNew - soit de délivrer une estimation 8(k,l) de la direction d'arrivée du transient, - soit d'indiquer qu'aucune estimation de direction d'arrivée ne peut être fournie, si ces règles ne sont pas satisfaites. 1 °) Significativité de PÀ,P1t,(8n,ax,k,l) (Onax étant l'angle qui maximise la valeur PNew (9,,k,l) ) Règle 1 : Une estimation de direction ne peut être fournie que si PNew (emax,k,1) dépasse un seuil donné PMJN , Cette première règle permet de s'assurer que sur la partie (k,l) du signal considéré, la probabilité de présence d'un transient et le niveau d'inter-corrélation sont assez élevés pour que l'estimation soit consistante. 2°) Monotonie de PNew sur l'intervalle 8s - enax; emax (pour alléger les notations, dans la suite les barres de module de PNe14, seront enlevées) Règle 2 : Si 8,nax est en dehors du cône privilégié, une estimation d'angle ne sera validée que siP1eV augmente de façon monotone sur l'intervalle r Les ù emax ; ema,j. 35 Cette deuxième règle analyse le contenu du "cône privilégié", correspondant au secteur angulaire sur lequel est centré la source s et qui présente une étendue angulaire de 80. Ce cône privilégié est défini par les angles 8 tels que 8ûOs û Le "bruit latéral" correspondra à un signal dont la direction d'arrivée est extérieure au cône privilégié, et l'on considérera donc qu'un bruit latéral est présent si 8,nax û 8S dépasse le seuil 80 . Pour valider cette détection d'un bruit latéral, il faut vérifier qu'un signal de parole utile ne se trouve pas simultanément à l'entrée du système. Pour cela, PNew(8lnax,k,l) est confronté aux valeurs de PNe,t,(8;,k,l) obtenues pour d'autres angles, notamment ceux qui appartiennent au cône privilégié. La règle permet ainsi de s'assurer qu'il n'y a pas de maximum local dans le cône privilégié. 3°) Fiabilisation de la détection d'un bruit latéral Règle 3 : Si B,nax est en dehors du cône privilégié pour la première fois sur la trame Z considérée, une estimation d'angle ne sera validée que si PNew(8max,k,l)a1 XPNew(Bmax,k,lû1) , et si PNew(8max,k,1)~a2 X 1 PNew(Bmax,k,i). 5 ie[lù5;1ù1] Si un bruit latéral est détecté, cette troisième règle tient compte des trames précédentes pour éviter les faux déclenchements. Elle ne s'applique qu'à la première trame d'un bruit latéral présumé, et vérifie que PNew(8,nax,k,l) augmente de façon significative par rapport aux données correspondantes obtenues sur les cinq trames précédentes. Les paramètres a1 et a2 sont choisis de manière à correspondre à des situations typiques, proches de la réalité. Si les trois règles 1 à 3 ci-dessus sont vérifiées, l'estimation 8(k,l) de la direction d'arrivée sera donnée par : B(k,l) = 8,ax 4°) Stabilisation de la détection d'un bruit latéral : Les deux dernières règles sont destinées à empêcher les coupures dans la détection d'un bruit latéral. Après une période de détection, elles conti- nuent à maintenir cet état pendant un laps de temps dit de hangover, quand bien même les règles de décision précédentes ne seraient plus vérifiées. Cela permet de détecter les éventuelles périodes à basse énergie d'un bruit non-stationnaire. Règle 4 : Si 9(k, l -1) est en dehors du cône privilégié (trame précédente), Si cpt, S HangoverTime,, (i.e. la période de Hangover n'est pas terminée), et si PNew (9(k, l -1), k, l) est supérieur à un seuil donné P alors 15 l'estimation d'angle est maintenue et cpt, est incrémenté. Règle 5 : Si 0(k, l -1) est en dehors du cône privilégié (trame précédente), Si cpt2 HangoverTime2 et si 1 PNew (8(k, l -1), k, i) est supérieur 20 à un seuil donné P2 alors l'és 6if)à/)on d'angle est maintenue et cpt2 est incrémenté. Si l'une de ces deux dernières règles (Règle n°4 ou n°5) est vérifiée, elle est prioritaire, et il en résulte : 8(k,l) = â(k,l -1) , donc avec correction é- 25 ventuelle de la valeur de 8(k,l), qui ne sera pas égale à 8max mais qui sera maintenue à sa valeur précédente. En résumé, le calcul de 8(k,l) suit trois cas possibles : (i) si la règle n°4 ou n°5 est vérifiée, alors 8(k,l) = 8(k,l -1) ; 30 (ii) dans le cas contraire (ni la règle n°4, ni la règle n°5 n'est vérifiée), si les règles n°1, n°2 et n°3 sont vérifiées, alors 8(k,l) = emax ; (iii) sinon (ni la règle n°4, ni la règle n°5 n'est vérifiée, et l'une au moins des règles n°1, n°2 et n°3 n'est pas vérifiée), alors 8^(k,l) n'est pas défini. 35 Dans une variante, l'estimateurPNe1,, est moyenné sur des paquets de bandes de fréquences K1,K2...,Kp : C. désignant le cardinal de K.. Dans ce cas, l'estimation d'angle 19. n'est pas faite sur chaque bande de fréquences, mais sur chaque paquet Ki de bandes de fréquences. On notera aussi qu'une approche "pleine bande" est possible (p =1, un seul angle étant estimé par trame). On notera enfin que la méthode proposée est compatible avec l'utilisation de micros unidirectionnels. Dans ce cas il sera courant d'utiliser un réseau linéaire (micros alignés et dont les directions privilégiées sont identiques) et orienté vers le locuteur. Dans ce cas la valeur de Os est donc naturellement connue et égale à zéro. Calcul d'une probabilité de présence de parole sur critère spatial (bloc 20) L'étape suivante, caractéristique du procédé de l'invention, consiste à cal-culer une probabilité de présence de parole basée sur l'estimation de di- rection d'arrivée 8(k,l) obtenue de la manière indiquée ci-dessus. Il s'agit d'une probabilité, notée pspa(k,l) , qui a donc pour originalité d'être calculée sur un critère spatial (à partir de 0(k,l) ), et qui permettra de distinguer parmi les signaux non-stationnaires la parole utile des bruits para-sites. Cette probabilité sera ensuite utilisée dans une structure classique de débruitage (bloc 22, décrit ci-après). La probabilité pspa(k,l) peut être calculée de différentes manières, donnant une valeur binaire ou bien des valeurs multiples. On donnera ci-dessous deux exemples de calcul pspa(k,l) , sachant que d'autres lois peuvent être utilisées pour exprimer pspa(k,l) à partir de 8(k,l) . PVew(0,,K1,1) = N(Nû1) Cl nu 1 °) Calcul d'une probabilité pspa (k, i) binaire La probabilité de présence de parole prendra les valeurs '0' ou '1' : - elle sera mise à '0' lorsqu'un bruit latéral, c'est-à-dire un transient pro- venant d'une direction extérieure au cône privilégié, est détecté ; - elle sera mise à '1' lorsque la direction d'arrivée du transient est à l'intérieur du cône privilégié, ou lorsqu'aucune estimation fiable n'a pu être faite sur cette direction. L'algorithme correspondant est le suivant : Si 8(k, i) est à l'intérieur du cône privilégié (9(k, 80 ), alors pspa (k, i) =1 l) - 8s 80 ), alors pspa (k, ) = 0 Si 8(k,l) n'est pas défini, alors pspa (k, l) = 1 2°) Calcul II est possible à valeurs continues dans [0;1] d'une probabilité pspa(k,l) un calcul progressif, par exemple d'utiliser pour pspa(k,l) selon l'algorithme suivant : Si 8(k,l) est à l'intérieur du cône privilégié ( 9(k,l) ù 8s 5 80 ), alors pspa (k, ) =1 8(k, l) ù 8s Si 8(k,l) n'est pas défini, alors pspa (k, l) =1 alors pspa (k, ) =1- Si 8(k,l) est à l'extérieur du cône privilégié ( 8(k,l)ùB0 Réduction de bruit latéral (bloc 22) La probabilité pspa(k,l) de présence de parole calculée au bloc 20, dé-pendant elle-même de la probabilité pT,.a,s,eät (k, l) de présence de tran- sients calculée au bloc 16, va être utilisée comme paramètre d'entrée dans une technique classique de débruitage. On sait que la probabilité de présence de parole est un estimateur crucial pour le bon fonctionnement d'un algorithme de débruitage, car elle sous-tend la bonne estimation du bruit et le calcul d'un gain optimal efficace. On peut avantageusement utiliser une méthode de débruitage de type OM-LSA (Optimally Modified - Log Spectral Amplitude) telle que celle décrite par :1. Cohen, Optimal Speech Enhancement Under Signal Presence Uncertainty Using Log-Spectral Amplitude Estimator, IEEE Signal Processing Letters, Vol. 9, No 4, April 2002. Essentiellement, l'application d'un gain nommé "gain LSA" (Log-Spectral Amplitude) permet de minimiser la distance quadratique moyenne entre le logarithme de l'amplitude du signal estimé et le logarithme de l'amplitude du signal de parole originel. Ce second critère se montre supérieur au premier car la distance choisie est en meilleure adéquation avec le corn- portement de l'oreille humaine et donne donc qualitativement de meilleurs résultats. Dans tous les cas, l'idée essentielle est de diminuer l'énergie des composantes fréquentielles très parasitées en leur appliquant un gain faible, tout en laissant intactes (par l'application d'un gain égal à 1) celles qui le sont peu ou pas du tout. L'algorithme "OM-LSA" (Optimally-Modified Log-Spectral Amplitude) améliore le calcul du gain LSA à appliquer en le pondérant par la probabilité conditionnelle de présence de parole. Dans cette méthode, la probabilité de présence de parole intervient à deux moments importants, pour l'estimation de l'énergie du bruit et pour le calcul du gain final, et la probabilité pspa(k,l) sera utilisée à ces deux ni-veaux. Si l'on note 2.Br,,;t(k,l)l'estimation de la densité spectrale de puissance du bruit, cette estimation est donnée par : ÂB,tuit(k,l)=aB,-u (k,l)•ÂBY�u(k,lû1)+{lûaBY�,t(k,l)1 X(k,l 2 avec : aBY,,it (k, l) = aB + (1û aB) •pspa (k, l ) On peut noter ici que la probabilité pspa(k,l) module le facteur d'oubli dans l'estimation du bruit, qui est mise à jour plus rapidement sur le signal brui-té X(k,l) lorsque la probabilité de parole est faible, ce mécanisme conditionnant entièrement la qualité de 2.8,.,ä t (k, l) . Le gain de débruitage Gom-csA(k,l) est donné par : Pspn (k,l ~ 1ùpspn (k>l) GOM LSA(k,l) = ~GHI(k,l)~ Gmin Gin (k, l) étant un gain de débruitage (dont le calcul dépend de l'estimation du bruit ÂB,.uit) décrit dans l'article précité de Cohen, et G,,,;n étant une constante correspondant au débruitage appliqué lorsque la 15 parole est considérée comme absente. On note ici que la probabilité pspa(k,l) joue un grand rôle dans la détermination du gain GoM-LsA(k,l) . Notamment, lorsque cette probabilité est nulle le gain est égal à et une réduction de bruit maximale est appliquée : si par exemple une valeur de 20 dB est choisie pourG,,,;n , les bruits 20 non-stationnaires précédemment détectés sont atténués de 20 dB. Le signal débruité S(k,l) en sortie du bloc 22 est donné par : S(k,l) = GoM_LsA(k,l).X(k,l) 25 On notera que d'ordinaire une telle structure de débruitage produit un résultat peu naturel et agressif sur les bruits non-stationnaires, qui sont confondus avec la parole utile. L'un des intérêts majeurs de la présente invention est d'éliminer efficacement ces bruits non-stationnaires. Par ailleurs, il est possible d'utiliser dans les expressions ci-dessus une 30 probabilité de présence de parole hybride phybÏld(k,l) , c'est-à-dire calculée à l'aide de pspa(k,l) combinée à une autre probabilité de présence de parole p(k,l), par exemple calculée selon la méthode décrite dans le WO 2007/099222 Al (Parrot SA). II vient : 1610 phybrid (k, l) = min(p(k, l), pspa (k, l)) Cette probabilité hybride permet de bénéficier du repérage des bruits non-stationnaires associé aux petites valeurs de pspa (k, l) , et de compléter l'estimation de la probabilité p,,,,brid(k,l) sur les parties (k,l) où l'estimation de direction d'arrivée 9(k,l) n'a pas été définie (produisant une probabilité pspa(k,l) forcée à la valeur 1 par sécurité). La probabilité hybride phvbrid(k,l) intègre ainsi à la fois les bruits non-stationnaires détectés par pspa(k,l) et les autres bruits (par exemple pseudo-stationnaires) détectés par p(k,l). Reconstitution temporelle du signal (bloc 24) La dernière étape consiste à appliquer au signal S(k,l) une transformée de Fourier rapide inverse iFFT pour obtenir dans le domaine temporel le signal de parole débruité s(t) . REVENDICATIONS1. Un procédé de débruitage d'un signal acoustique bruité capté par une pluralité de microphones d'un dispositif audio multi-microphone opérant dans un milieu bruité, notamment un dispositif téléphonique "mains libres" pour véhicule automobile, le signal acoustique bruité comprenant une composante utile de parole issue d'une source de parole directive et une composante parasite de bruit, cette composante de bruit incluant elle-même une composante de bruit latéral non stationnaire directif, procédé caractérisé en ce qu'il comporte, dans le domaine fréquentiel pour une pluralité de bandes de fréquences définies pour des trames temporelles successives de signal, les étapes de traitement du signal sui-vantes : a) combinaison (12) de la pluralité de signaux captés par la pluralité cor- respondante de microphones en un signal combiné bruité (X(k, l)) ; 15 b) à partir du signal combiné bruité, estimation (14) d'une composante de bruit pseudo-stationnaire (V(k,l)) contenue dans ce signal combi- né bruité ; c) à partir de la composante de bruit pseudo-stationnaire estimée à l'étape b) et du signal combiné bruité, calcul (16) d'une probabilité de 20 présence de transients (pT,.aäS;e,,,(k,l)) dans le signal combiné bruité ; d) à partir de la pluralité de signaux captés par la pluralité correspon- dante de microphones et de la probabilité de présence de transients calculée à l'étape c), estimation (18) d'une direction principale d'arri- vée des transients (9(k, l)) ; 25 e) à partir de la direction principale d'arrivée des transients estimée à l'étape d), calcul (20) d'une probabilité de présence de parole sur un critère spatial (pspa(k,l) ), propre à discriminer entre parole utile et bruit latéral parmi les transients ; f) à partir de la probabilité de présence de parole calculée à l'étape e) 30 et du signal combiné bruité, réduction sélective du bruit (22) par application d'un gain variable propre à chaque bande de fréquences et à chaque trame temporelle. 2. Le procédé de la revendication 1, dans lequel le traitement de l'étape a) est un traitement de préfiltrage de type fixed beamforming. 3. Le procédé de la revendication 1, dans lequel le traitement de l'étape d) comprend les sous-étapes successives suivantes : dl) partition de l'espace en une pluralité de secteurs angulaires ; d2) pour chaque secteur, évaluation d'un estimateur de direction d'arri- vée à partir de la pluralité de signaux captés par la pluralité corres- pondante de microphones ; d3) pondération de chaque estimateur par la probabilité de présence de transients calculée à l'étape c) ; d4) à partir des valeurs d'estimateurs pondérées calculées à l'étape d3), estimation d'une direction principale d'arrivée des transients ; d5) validation ou invalidation de l'estimation de la direction principale d'arrivée des transients opérée à l'étape d4). 4. Le procédé de la revendication 3, dans lequel à l'étape d5) l'estimation n'est validée que si la valeur de l'estimateur pondéré correspondant à la direction estimée est supérieure à un seuil prédéterminé. 5. Le procédé de la revendication 3, dans lequel à l'étape d5) l'estimation n'est validée qu'en l'absence de maximum local de l'estimateur pondéré dans le secteur angulaire d'origine du signal de parole utile. 25 6. Le procédé de la revendication 3, dans lequel à l'étape d5) l'estimation n'est validée que si la valeur de l'estimateur est croissante de façon monotone sur une pluralité de trames temporelles successives. 7. Le procédé de la revendication 3, comprenant en outre une étape de 30 maintien de l'estimation de la direction principale d'arrivée pendant un laps de temps minimal prédéterminé. 8. Le procédé de la revendication 1, dans lequel la probabilité de présence de parole calculée à l'étape e) est une probabilité binaire, prenant 35 une valeur 1 ou 0 selon que la direction principale d'arrivée des transients20 estimée à l'étape d) est située ou non dans le secteur angulaire d'origine du signal de parole utile. 9. Le procédé de la revendication 1, dans lequel la probabilité de pré- sence de parole calculée à l'étape e) est une probabilité à valeurs multiples, fonction de l'écart angulaire entre la direction principale d'arrivée des transients estimée à l'étape d) et la direction d'origine du signal de parole utile. 10. Le procédé de la revendication 1, dans lequel le traitement de l'étape f) est un traitement de réduction sélective du bruit par application d'un gain à amplitude log-spectrale modifié optimisé OM-LSA.