Brevet FR3051958A1 — publié le 20171201

Procédé et dispositif pour estimer un signal déréverbéré. DOMAINE DE L'INVENTION La présente invention est relative aux procédés et 5 dispositifs pour estimer la phase instantanée d'un signal acoustique déréverbéré et pour estimer un signal déréverbéré. ARRIERE PLAN DE L'INVENTION Lorsqu'un signal acoustique originel est émis dans 10 un milieu réverbérant puis capté par un microphone, le microphone capte un signal réverbéré dépendant du milieu réverbérant. Dans la suite, on entend par « signal acoustique anéchoïque » le signal acoustique originel non réverbéré 15 par un milieu. Un signal acoustique anéchoïque peut parfois être directement enregistré par un microphone, par exemple lorsque le signal acoustique originel est émis dans une chambre anéchoïque. Toutefois, dans des conditions courantes 20 d'enregistrement un microphone enregistre un signal acoustique réverbéré qui est un signal constitué du signal acoustique originel reçu directement, mais également des réflexions du signal acoustique originel sur les éléments réverbérant du milieu, par exemple les murs d'une salle. 25 Une forte réverbération acoustique du milieu peut être particulièrement gênante puisqu'elle dégrade la qualité du son enregistré et réduit l'intelligibilité de la parole et sa reconnaissance par des machines. Pour résoudre ce problème, on connait des procédés 30 et dispositifs qui permettent de reconstruire l'amplitude d'un signal déréverbéré à partir d'un signal acoustique réverbéré par un milieu. Dans la présente demande, on entend par « signal 3051958 2 déréverbéré » une estimation du signal acoustique originel, ou signal anéchoïque, obtenue par traitement analogique ou digital d'un signal acoustique réverbéré enregistré par un microphone. A titre d'exemple, le document US2016035367 décrit un procédé de déreverbération qui permet de reconstruire un signal déréverbéré à partir d'un signal acoustique réverbéré par un milieu en calculant l'amplitude du signal déréverbéré dans plusieurs bandes de fréquences. Il existe un besoin pour améliorer encore les performances de tels procédés en estimant de manière plus précise les caractéristiques du signal déréverbéré à partir d'un signal acoustique réverbéré enregistré par un microphone. La présente invention vient ainsi améliorer cette situation. OBJETS ET RESUME DE L'INVENTION A cet effet, l'invention a pour premier objet un procédé pour estimer une phase instantanée de signal 20 acoustique déréverbéré, le procédé comprenant les étapes suivantes : (a) on mesure un signal acoustique réverbéré par une propagation dans un milieu, (b) on estime au moins une fréquence instantanée 25 lissée du signal réverbéré ainsi qu'un taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré, (c) on calcule au moins une fréquence instantanée de signal déréverbéré à partir de ladite fréquence 30 instantanée lissée du signal acoustique réverbéré, du taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré et d'un facteur d'influence du milieu, ledit facture d'influence étant 3051958 3 fonction d'un temps de réverbération dudit milieu, (d) on détermine au moins une phase instantanée de signal déréverbéré en intégrant la fréquence instantanée de signal déréverbéré au cours du temps. 5 Grâce à ces dispositions, on peut obtenir une estimation de la phase du signal déréverbéré proche de la phase du signal anéchoïque. Cette phase du signal déréverbéré peut être utilisée dans des procédés aval, en particulier en remplacement de la phase du signal 10 réverbéré, pour fournir des résultats de meilleurs qualité, par exemple pour fournir un signal vocal déréverbéré plus intelligible ou encore pour permettre une localisation d'une source sonore plus précise que les procédé de l'état de l'art. 15 Dans des modes de réalisation préférés de l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes : pour calculer ladite au moins une fréquence instantanée de signal déréverbéré, on détermine un facteur de correction en multipliant le taux de variation au cours du temps de la fréquence instantanée lissée du signal réverbéré par le facteur d'influence du milieu, en particulier on additionne ledit facteur de correction avec ladite fréquence instantanée lissée du 25 signal acoustique réverbéré ; - le facteur d'influence du milieu est donné par : 1 min(t, Th) R(t) = 28 + 1 - e2Smin(t,Th) où ô et Th sont respectivement un facteur d'amortissement et une durée d'une décroissance 30 exponentielle p(t)=e-st1[0,Thi de la réponse impulsionnelle du milieu ; 3051958 4 - on estime le facteur d'amortissement ô et la durée de la réponse impulsionnelle Th à partir d'un temps de réverbération mesuré dans le milieu, notamment un temps de réverbération RT60, par exemple en calculant 5 ô = 3.1og(10)/RT60 et Th = 1.5RT60 ; - pour chaque bande fréquentielle k d'une pluralité de N bandes fréquentielles, on estime une fréquence instantanée lissée du signal réverbéré dans ladite bande fréquentielle k ainsi 10 qu'un taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré, on calcule une fréquence instantanée de signal déréverbéré dans ladite bande fréquentielle k à partir de ladite fréquence instantanée du signal acoustique 15 réverbéré, du taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré et du facteur d'influence du milieu, et on détermine une phase instantanée de signal déréverbéré dans ladite bande fréquentielle k en intégrant 20 la fréquence instantanée de signal déréverbéré dans la bande fréquentielle k au cours du temps ; - pour estimer une fréquence instantanée lissée du signal réverbéré pour chaque bande fréquentielle k de la pluralité de N bandes fréquentielles, on calcule une transformée de Fourier locale discrète du signal acoustique réverbéré en utilisant des fonctions fenêtres w(n) avec n entre 0 et N-1, et pour déterminer la phase instantanée de signal déréverbéré dans chaque bande fréquentielle k de ladite pluralité de N bandes fréquentielles, on intègre la fréquence instantanée de signal déréverbéré au cours du temps en ajoutant un terme de correction fonction des fonctions fenêtres w(n) ; 3051958 5 - pour estimer une fréquence instantanée lissée du signal réverbéré pour chaque bande fréquentielle k de la pluralité de N bandes fréquentielles, on applique un algorithme de vocodeur réassigné utilisant une transformée 5 de Fourier locale discrète du signal acoustique réverbéré puis on applique un filtre de lissage temporelle; - au cours d'une étape de calibration, on mesure un signal acoustique de référence réverbéré par une propagation dans le milieu et on détermine le facteur 10 d'influence du milieu à partir dudit signal acoustique de référence, en particulier on détermine un temps de réverbération dudit milieu. L'invention a également pour objet un procédé pour reconstruire un signal déréverbéré dans lequel, 15 on mesure un signal acoustique réverbéré par une propagation dans un milieu, on détermine un spectre d'amplitude de signal déréverbéré pour une pluralité de N bandes fréquentielles à partir du signal acoustique réverbéré, 20 on estime une phase instantanée de signal déréverbéré pour chaque bande fréquentielle k de ladite pluralité de N bandes fréquentielles au moyen d'un procédé tel que décrit ci-avant, et on reconstruit un signal déréverbéré à partir dudit 25 spectre d'amplitude de signal déréverbéré et desdites phases instantanées de signal déréverbéré. L'invention a enfin pour objet un dispositif pour estimer une phase instantanée de signal acoustique déréverbéré, comprenant : 30 (a)des moyens de mesure pour capter au moins un signal acoustique réverbéré par une propagation dans un milieu, (b)des moyens pour estimer au moins une fréquence 3051958 6 instantanée lissée du signal réverbéré ainsi qu'un taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré, (c)des moyens pour calculer au moins une fréquence 5 instantanée de signal déréverbéré à partir de ladite fréquence instantanée lissée du signal acoustique réverbéré, du taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré et d'un facteur d'influence du milieu, ledit facture d'influence 10 étant fonction d'un temps de réverbération dudit milieu, (d)des moyens pour déterminer au moins une phase instantanée de signal déréverbéré en intégrant la fréquence instantanée lissée de signal déréverbéré au cours du temps. BREVE DESCRIPTION DES DESSINS 15 D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante d'un de ses modes de réalisation, donné à titre d'exemple non limitatif, en regard des dessins joints. Sur les dessins : 20 - la figure 1 est une vue schématique illustrant la réverbération du son dans une pièce lorsqu'un sujet parle de façon que ses paroles soient captées par un dispositif selon une forme de réalisation de l'invention, - la figure 2 est un schéma de principe du 25 dispositif de la figure 1, et - la figure 3 est un organigramme d'un procédé de reconstruction de signal déréverbéré selon un mode de réalisation de l'invention, mettant en particulier en oeuvre un procédé d'estimation d'une phase instantanée de signal 30 déréverbéré selon une forme de réalisation de l'invention. DESCRIPTION PLUS DETAILLEE Sur les différentes figures, les mêmes références désignent des éléments identiques ou similaires. 3051958 7 L'invention a pour but d'estimer une phase instantanée de signal acoustique déréverbéré à partir d'une mesure d'un signal acoustique réverbéré par une propagation dans un milieu 7, par exemple une pièce dans un bâtiment 5 telle que schématisée sur la figure 1. L'invention permet ainsi de traiter les signaux acoustiques captés par un dispositif électronique 1 pourvu d'un microphone 2. Le dispositif électronique 1 peut être par exemple un téléphone dans l'exemple représenté, ou un 10 ordinateur ou autre. Lorsqu'un son est émis dans le milieu 7, par exemple par une personne 3, ce son se propage jusqu'au microphone 2 selon divers trajets 4, soit directs, soit après réflexion sur une ou plusieurs parois 5, 6 du milieu 15 7. Comme représenté sur la figure 2, le dispositif électronique 1 peut comporter par exemple une unité centrale électronique 8 telle qu'un processeur ou autre, reliée au microphone 2 et à divers autre éléments, incluant 20 par exemple un haut-parleur 9, un clavier 10, un écran 11. L'unité centrale électronique 8 peut communiquer avec un réseau externe 12, par exemple un réseau téléphonique. L'invention permet au dispositif électronique 1 d'estimer une phase instantanée de signal acoustique 25 déréverbéré. Dans une première application de grand intérêt, la phase instantanée de signal déréverbéré peut être utilisée pour reconstruire un signal déréverbéré à partir d'un signal acoustique réverbéré. 30 Pour cela, on mesure tout d'abord un signal acoustique réverbéré par une propagation dans le milieu. Puis, on détermine un spectre d'amplitude de signal déréverbéré pour une pluralité de N bandes fréquentielles à 3051958 8 partir du signal acoustique réverbéré. De nombreux procédés permettant de déterminer un spectre d'amplitude de signal déréverbéré à partir d'un signal acoustique réverbéré sont connus de l'art antérieur. 5 Ces procédés consistent par exemple à estimer un spectre de réverbération à partir du signal acoustique réverbéré puis à soustraire ledit spectre de réverbération du signal acoustique réverbéré. On connait ainsi des procédés de détermination d'un 10 spectre d'amplitude de signal déréverbéré utilisant : - la prédiction à long-terme tel que décrit dans l'article "Suppression of late reverberation effect on speech signal using long-term multiple-step linear prediction" de K. Kinoshita, M. Delcroix, T. Nakatani, et 15 M. Miyoshi paru dans IEEE Transactions on Audio, Speech, and Language Processing, vol. 17, no. 4, pp. 534-545, Mai 2009, - une modélisation stochastique de la réponse impulsionelle du milieu tel que décrit dans "A new method 20 based on spectral subtraction for speech dereverberation" de K. Lebart et J.M. Boucher, paru dans ACUSTICA, vol. 87, no. 3, pp. 359-366, 2001, ou encore - des réseaux de neurones profonds tel que décrit dans "Speech dereverberation for enhancement and 25 recognition using dynamic features constrained deep neural networks and feature adaptation" de X. Xiao, S. Zhao, D.H. Ha Nguyen, X. Zhong, D.L. Jones, E.S. Chang, et H. Li, paru dans EURASIP Journal on Advances in Signal Processing, vol. 2016, no. 1, pp. 1-18, 2016. 30 Dans ces procédés de l'art antérieur, on reconstitue ensuite un signal déréverbéré à partir du spectre d'amplitude de signal déréverbéré obtenu et de la phase du signal réverbéré. 3051958 9 Il existe toutefois un besoin pour améliorer encore la qualité et l'intelligibilité du signal déréverbéré obtenu par ce procédé. A cette fin, selon l'invention, on détermine par 5 ailleurs, à partir du signal acoustique réverbéré une phase instantanée de signal déréverbéré pour chaque bande fréquentielle k de la pluralité de N bandes fréquentielles au moyen d'un procédé tel que décrit ci-après. Puis, on reconstruit un signal déréverbéré à partir 10 du spectre d'amplitude de signal déréverbéré et de la phase estimée à l'aide de la méthode selon l'invention. On obtient de cette manière un signal déréverbéré reconstruit de qualité nettement supérieure. La phase instantanée de signal déréverbéré 15 déterminée par le procédé selon l'invention peut également avoir d'autres utilisations que la reconstruction du signal déréverbéré et peut par exemple être utilisée pour améliorer la qualité et la précision d'un algorithme de localisation de source sonore tel que connu par ailleurs 20 dans la littérature. Il est connu que le milieu réverbérant peut être modélisé par un modèle stochastique en définissant une réponse impulsionnelle h(t) de la forme : h(t) = b(t)p(t) (1) 25 dans laquelle b(t)-N(0,0-2) est un bruit blanc avec une distribution gaussienne centrée de variance 0-2 et p(C)= e-el[uli] est une décroissance exponentielle de la réponse impulsionnelle du milieu où ô et Th sont respectivement un facteur d'amortissement et une durée de la réponse 30 impulsionnelle du milieu. Un tel modèle stochastique est par exemple exposé dans la thèse de J.D. Polack, « La transmission de 3051958 10 l'énergie sonore dans les salles », soutenue à l'Université du Maine en 1988. Le facteur d'amortissement ô et la durée de la réponse impulsionnelle Th peuvent être déterminés à partir 5 d'un temps de réverbération mesuré dans le milieu. Un temps de réverbération couramment utilisé est le temps de réverbération à 60dB dénoté RT60. Le temps de réverbération à 60dB correspond au temps requis pour que la courbe de décroissance d'énergie (Energy Decay Curve) EDC 10 diminue de 60 dB. Le temps de réverbération à 60dB peut par exemple être défini par la méthode d'intégration inverse de Manfred R. Schroeder (New Method of Measuring Reverberation Time, The Journal of the Acoustical Society of Ameri.ca, 37(3) :409, 15 1965) par la courbe de décroissance d'énergie (Energy Decay Curve) EDC(n)=ENk_hnh(k)2 où h est la réponse impulsionnelle d'un milieu de longueur Nh et n est un indice temporel, par exemple un nombre d'échantillons obtenus par un échantillonnage de pas temporel constant, n étant compris 20 entre 1 et Nh. RT60 est alors le temps à l'indice temporel n requis pour que EDC(n) diminue de 60 dB. Des valeurs typiques du temps de réverbération RT60 sont par exemple des valeurs entre 0.4 s et 2 s. Bien que le temps de réverbération RT60 soit le plus 25 couramment utilisé, on peut également utiliser un autre temps de réverbération caractéristique du milieu 7. On peut ensuite calculer le facteur d'amortissement du milieu ô à partir du temps de réverbération RT60 par la formule ô = 3.1og(10)/RT60 . 30 La durée de la réponse impulsionnelle Th peut également être définie à partir du temps de temps de réverbération comme Th = 1.5RT60. 3051958 11 Toutefois, le facteur d'amortissement du milieu ô et la durée de la réponse impulsionnelle Th peuvent également être calculés par d'autres méthodes connues de l'art antérieur. 5 A partir du modèle statistique donné par l'équation (1), le signal acoustique réverbéré peut être relié au signal acoustique anéchoïque par l'équation de convolution : y(t) = (h *s)(t) (2) 10 où y(t) est le signal acoustique réverbéré et s(t) est le signal acoustique anéchoïque. La phase instantanée du signal réverbéré peut par ailleurs être exprimée en fonction de la transformée de Hilbert du signal réverbéré comme : 15(Prev = ami-2n (y(t)) (3) y(t) où (pre'(t) est la phase instantanée du signal réverbéré et .9(0 est la transformée de Hilbert du signal réverbéré. Sachant que l'on peut relier la fréquence 20 instantanée du signal réverbéré à la phase instantanée du signal réverbéré par l'expression : frev(0= dePrev(t) 27r dt (t) ( 4) Il est alors possible d'exprimer la fréquence instantanée du signal anéchoïque en fonction de l'espérance 25 mathématique de la fréquence instantanée du signal réverbéré à partir des équations (1) à (4) comme : f + f + 1-e E[frev(t)] ( min(t,Th) 28min(t,T h)) ( 5 où f(t) est la fréquence instantanée du signal anéchoïque estimée à l'instant t, E[fre'(t)] est l'espérance 30 mathématique de la fréquence instantanée du signal 3051958 12 réverbéré à l'instant t et f est le taux de variation au cours du temps de la fréquence instantanée du signal réverbéré. L'espérance mathématique de la fréquence 5 instantanée du signal réverbéré à l'instant t ne peut pas être mesurée mais peut être approximée par un lissage temporelle de la fréquence instantanée du signal réverbéré mesurée. On peut ainsi estimer une fréquence instantanée 10 d'un signal déréverbéré en fonction d'une fréquence instantanée du signal réverbéré à partir des équations (1) à (5) comme : min,Th) i(t) = frev(t) (-218 + 1-e2S(i(t,T h)) (6) où f(t) est la fréquence instantanée du signal 15 déréverbéré estimée à l'instant t, f",'(0 est une fréquence instantanée lissée du signal réverbéré à l'instant t et f est le taux de variation au cours du temps de la fréquence instantanée lissée du signal réverbéré. L'équation (6) permet d'estimer une fréquence instantanée du signal 20 déréverbéré en fonction de la fréquence instantanée lissée du signal réverbéré, du taux de variation au cours du temps de la fréquence instantanée et d'un facteur d'influence du milieu R(t) est donné par I?(t) =:L_f min(t,Th) (7) 28 1-e 25 On peut ainsi réécrire l'équation (6) comme : /(0 = frev(t) f R(t) ( 8) Une phase instantanée du signal déréverbéré 0(0 peut par la suite être déterminée par intégration temporelle comme : 30 0(0 = 2if fot (t)cit (P (0) (9) 2 8min(t,Th) 3051958 13 où 0(0) est une phase originelle du signal déréverbéré. La fréquence et la phase du signal déréverbéré estimées au moyen des équations (6) à (9) sont donc des 5 estimations de la fréquence et de la phase du signal acoustique originel, ou signal anéchoïque. Les tests réalisés par les inventeurs indiquent que ces estimations sont particulièrement bonnes puisqu'elles conduisent à un signal déréverbéré de qualité nettement 10 supérieure à l'état de l'art. Un procédé pour estimer une phase instantanée de signal acoustique déréverbéré selon l'invention comporte ainsi les étapes suivantes : (a) une étape de mesure, au cours de laquelle on 15 mesure le signal acoustique réverbéré par une propagation dans un milieu, (b) une étape d'estimation au cours de laquelle on estime au moins une fréquence instantanée lissée du signal réverbéré ainsi qu'un taux de variation au cours du temps 20 de ladite fréquence instantanée lissée du signal réverbéré, (c) une étape de calcul au cours de laquelle on calcule au moins une fréquence instantanée de signal déréverbéré à partir de ladite fréquence instantanée lissée du signal acoustique réverbéré, du taux de variation au 25 cours du temps de ladite fréquence instantanée lissée du signal réverbéré et d'un facteur d'influence du milieu, ledit facteur d'influence étant fonction d'un temps de réverbération dudit milieu, (d) une étape de détermination au cours de laquelle on 30 détermine au moins une phase instantanée de signal déréverbéré en intégrant la fréquence instantanée de signal déréverbéré au cours du temps. (a) Etape de mesure : 3051958 14 Au cours de cette étape, le microphone 2 capte un signal acoustique réverbéré par une propagation dans le milieu 7, par exemple lorsque le locuteur 3 parle. Ce signal est échantillonné et stocké dans le processeur 8 ou 5 une mémoire annexe (non représentée). Comme indiqué ci-avant, le signal capté y(t) est une convolution du signal anéchoïque s(t) émis (parole) avec la réponse impulsionnelle h(t) du milieu entre le locuteur 3 et le microphone 2, 10 (b) Etape d'estimation: Au cours de cette étape, on estime au moins une fréquence instantanée lissée du signal réverbéré ainsi qu'un taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré à partir 15 du signal acoustique réverbéré mesuré à l'étape (a). Dans un mode de réalisation de l'invention, on peut déterminer la fréquence instantanée lissée du signal réverbéré en mesurant tout d'abord la fréquence instantanée du signal réverbéré puis en lissant ladite fréquence 20 instantanée par exemple par un lissage temporel au moyen d'un filtre de Savitzky-Golay. La fréquence instantanée du signal réverbéré peut être déterminée, de manière générale, par une transformation de Fourier du signal. 25 Dans une variante de réalisation, pour chaque bande fréquentielle k d'une pluralité de N bandes fréquentielles, on peut estimer une fréquence instantanée du signal réverbéré dans ladite bande fréquentielle k ainsi qu'un taux de variation au cours du temps de ladite fréquence 30 instantanée du signal réverbéré. Pour cela, on peut par exemple appliquer un algorithme de vocodeur réassigné utilisant une transformée de Fourier locale discrète du signal acoustique réverbéré. 3051958 15 Un tel algorithme de vocodeur réassigné est par exemple décrit dans le document "Estimation of frequency for AM/FM models using the phase vocoder framework" de M. Betser, P. Collen, G. Richard, et B. David, paru dans IEEE 5 Transactions on Signal Processing, vol. 56, no. 2, pp. 505- 517, Feb 2008. Pour mettre en oeuvre un tel algorithme, on peut en particulier calculer une transformée de Fourier locale discrète du signal acoustique réverbéré en utilisant des 10 fonctions fenêtres w(n) avec n entre 0 et N-1. Une telle transformée de Fourier locale discrète du signal acoustique réverbéré peut être mise en oeuvre avec des fonctions fenêtres w(n) de taille N et des trames temporelles séparées par des sauts de R échantillons de 15 signal. Le signal acoustique réverbéré étant échantillonné avec une fréquence f' par exemple 16 kHz, on obtient ainsi N fréquences discrètes fk =kLs, k E [0, N - 1] et N f trames N temporelles. N vaut par exemple 256, 512 ou 1024. R vaut 20 par exemple la moitié ou le quart de N. L'estimation de la, ou des, fréquence(s) instantanée(s) du signal réverbéré peut typiquement être fait sur un nombre NI, de trames, par exemple une centaine de trames, correspondant à au moins quelques secondes de 25 signal selon les paramètres d'analyse choisis. Les trames peuvent avoir une durée individuelle de 10 à 100 ms, notamment de l'ordre de 32 ms. Les trames peuvent se recouvrir mutuellement, par exemple avec un taux de recouvrement de l'ordre de 50 % entre trames successives. 30 Une fois les fréquences instantanées du signal réverbéré estimées, on peut alors les lisser par un algorithme de lissage temporel comme indiqué ci-dessus pour obtenir les fréquences instantanées lissées du signal 3051958 16 réverbéré. (c) Etape de calcul : Au cours de cette étape, on calcule l'équation (8) ci- dessus f(t) = f ' v (t) + f R (t) pour estimer une fréquence 5 instantanée du signal déréverbéré. Dans la variante de réalisation dans laquelle on a estimé une fréquence instantanée lissée du signal réverbéré pour chaque bande fréquentielle k d'une pluralité de N bandes fréquentielles, on peut alors plus précisément 10 calculer une fréquence instantanée de signal déréverbéré F(774k) dans chaque bande fréquentielle k et pour chaque trame temporelle m. Plus précisément, on calcule alors la fréquence instantanée de signal déréverbéré P(774k) à partir de la 15 fréquence instantanée lissée du signal acoustique réverbéré de ladite bande fréquentielle k, du taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré et du facteur d'influence du milieu R(t). Ce calcul utilise également l'équation (8) qui est 20 appliquée de manière indépendante à chaque bande fréquentielle k, c'est-à-dire en remplaçant f(t) par P(174k). Pour estimer la fréquence instantanée du signal déréverbéré f(t) ou P(m,k), on détermine tout d'abord un facteur de correction iR(t) en multipliant le taux de 25 variation au cours du temps f de la fréquence instantanée lissée du signal réverbéré par le facteur d'influence du milieu R(t) = 1/(28) + min(t, Th)/(1 - exp(28min(t, Th))) . Puis, on additionne le facteur de correction iR(t) avec la fréquence instantanée lissée du signal acoustique 30 réverbéré selon l'équation (8). Le facteur d'influence du milieu R(t) peut être déterminé préalablement dans une étape de calibration 3051958 17 préliminaire. Au cours de cette étape de calibration préliminaire, on mesure un signal acoustique de référence réverbéré par une propagation dans le milieu et on détermine le facteur 5 d'influence du milieu à partir dudit signal acoustique de référence. Pour cela on peut par exemple déterminer un temps de réverbération dudit milieu par des méthodes connues par ailleurs, par exemple le temps de réverbération RT60 tel 10 que décrit ci-avant, et en déduire le facteur d'amortissement ô et la durée de la réponse impulsionnelle Th. Le signal acoustique de référence peut être un signal acoustique réverbéré par le milieu à partir d'un signal 15 originel connu du dispositif. Toutefois, la détermination du facteur d'influence du milieu peut également être réalisée « à l'aveugle », c'est-à-dire à partir d'un signal réverbéré enregistré suite à un signal originel arbitraire. 20 Avantageusement, on peut utiliser une pluralité de signaux acoustiques de référence qui correspondent à une pluralité respective de cas de figures différents (locuteurs différents, positions différentes, milieux 7 différents). Le nombre de signaux acoustiques de référence 25 peut être de plusieurs centaines, voire plusieurs milliers. Dans un mode de réalisation particulier de l'invention, le signal acoustique de référence peut être constitué par le signal acoustique réverbéré utilisé par le procédé selon l'invention, de sorte que la détermination du 30 facteur d'influence du milieu est alors réalisée directement lors de la mise en oeuvre du procédé d'estimation de la phase instantanée et sans nécessiter d'étape de calibration préliminaire. 3051958 18 La détermination du facteur d'influence du milieu peut également être réalisée de façon répétitive, de façon que le dispositif 1 s'adapte par exemple à des changements de locuteurs 3, à des déplacements de locuteurs 5 3, à des déplacements du dispositif 1 ou d'autres objets dans le milieu 7. (d) Etape de détermination : Enfin, au cours de cette étape on détermine la phase instantanée du signal déréverbéré O(t) par une 10 intégration temporelle de la fréquence instantanée déréverbéré comme indiqué dans l'équation (9). Cette intégration temporelle peut être réalisée à partir d'une phase originel du signal déréverbéré 0(0). Dans la plupart des cas, le signal déréverbéré 15 peut-être supposé avoir une phase égale à la phase du signal réverbéré à l'origine, de sorte que l'on peut par exemple prendre 0(0) = (Prev(0). Ceci s'applique en particulier dans le cas dans lequel le signal enregistré est précédé par un silence, de sorte que la réverbération est nulle à 20 l'origine. Ici également, en variante, on peut déterminer une phase instantanée de signal déréverbéré 0(n4k) dans chaque bande fréquentielle k de la pluralité de N bandes fréquentielles et pour chaque trame temporelle m en 25 intégrant la fréquence instantanée de signal déréverbéré de ladite bande fréquentielle k au cours du temps, c'est-à-dire en la sommant sur les trames temporelles m. Lorsque, pour estimer une fréquence instantanée lissée du signal réverbéré pour chaque bande fréquentielle 30 k de la pluralité de N bandes fréquentielles, on a calculé une transformée de Fourier locale discrète du signal acoustique réverbéré en utilisant des fonctions fenêtres w(n) avec n entre 0 et N-1, il est nécessaire de tenir 3051958 19 compte desdites fonctions fenêtres w(n) pour le calcul de la phase instantanée du signal anéchoïque (p(0. On a ainsi: mR mR (13(m, k) = (-) + arg r(k,f )) fs fs où 94mR ---) est la phase de Hilbert telle que définie par 5 l'équation (3) pour la trame temporelle d'indice m, cre(n4k) est la phase du signal anécoïque et F(k,f) est un facteur de correction lié aux fonctions fenêtres w(n) qui peut par exemple s'écrire : N-1 [27r(f - fk)n Tri(fsn)21) r(k,f) w(n) exp n=0 L'intégration temporelle des fréquences instantanées 10 déterminées pour le signal déréverbéré peut alors s'écrire comme une somme sur les trames temporelles cie(m, k) =(m - 1, k) + R 7(k, f. ((m - 1)1'? 2n-r. (m, k) Js Js + arg F (-1F* fs ))/ où k) est la fréquence instantanée de signal déréverbéré pour la bande fréquentielle k et pour la trame temporelle m et F* dénote le complexe conjugué du facteur 15 de correction F lié aux fonctions fenêtres w(n). De manière analogue au cas, détaillé ci-avant, dans lequel une seule fréquence instantanée lissée est déterminée, on peut par exemple initialiser I(),k) pour chaque bande fréquentielle k avec la valeur cre''(0,k), 20 c'est-à-dire considérer la réverbération nulle à l'origine. Les tests effectués montrent que l'utilisation de la phase estimée du signal déréverbéré dans les algorithmes de 3051958 20 reconstruction de signal déréverbéré et de localisation de sources, au lieu de la phase du signal réverbéré qui est traditionnellement utilisée, offre une qualité et une intelligibilité de signal déréverbéré 5 nettement améliorée ainsi qu'une meilleure localisation des sources sonores. Ainsi par exemple des tests ont montré un accroissement de 10 dB du rapport signal sur réverbération ( « signal-to-reverberant ratio » ou SRR) et une diminution 10 de 5 dB de la distance Cepstrale ( « cepstral distance » OU CD) ce qui correspond respectivement à un gain important de déreverbération et une réduction significative de la distorsion . REVENDICATIONS1. Procédé pour estimer une phase instantanée de 5 signal acoustique déréverbéré, le procédé comprenant les étapes suivantes : (a) on mesure un signal acoustique réverbéré par une propagation dans un milieu, (b) on estime au moins une fréquence instantanée 10 lissée du signal réverbéré ainsi qu'un taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré, (c) on calcule au moins une fréquence instantanée de signal déréverbéré à partir de ladite fréquence 15 instantanée lissée du signal acoustique réverbéré, du taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré et d'un facteur d'influence du milieu, ledit facture d'influence étant fonction d'un temps de réverbération dudit milieu, 20 (d) on détermine au moins une phase instantanée de signal déréverbéré en intégrant la fréquence instantanée de signal déréverbéré au cours du temps. 2. Procédé selon la revendication 1, dans lequel, pour calculer ladite au moins une fréquence instantanée de 25 signal déréverbéré, on détermine un facteur de correction en multipliant le taux de variation au cours du temps de la fréquence instantanée lissée du signal réverbéré par le facteur d'influence du milieu, en particulier dans lequel on additionne ledit 30 facteur de correction avec ladite fréquence instantanée lissée du signal acoustique réverbéré. 3. Procédé selon la revendication 1 ou 2, dans lequel, le facteur d'influence du milieu est donné par : 3051958 22 1 min(t, Th) R(t) = 28 + 1 - e2Smin(t,Th) où ô et Th sont respectivement un facteur d'amortissement et une durée d'une décroissance exponentielle p(t)=e-st1[0,Thi de la réponse impulsionnelle du milieu. 4. Procédé selon la revendication 3, dans lequel, on estime le facteur d'amortissement ô et la durée de la réponse impulsionnelle Th à partir d'un temps de réverbération mesuré dans le milieu, notamment un temps de réverbération RT60, par exemple en calculant ô = 3 .log(10) /RT60 et Th = 1.5RT60 . 5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel, pour chaque bande fréquentielle k d'une pluralité de N bandes fréquentielles, on estime une fréquence instantanée lissée du signal réverbéré dans ladite bande fréquentielle k ainsi 15 qu'un taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré, on calcule une fréquence instantanée de signal déréverbéré dans ladite bande fréquentielle k à partir de ladite fréquence instantanée lissée du signal acoustique 20 réverbéré, du taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré et du facteur d'influence du milieu, et on détermine une phase instantanée de signal déréverbéré dans ladite bande fréquentielle k en intégrant 25 la fréquence instantanée de signal déréverbéré dans la bande fréquentielle k au cours du temps. 6. Procédé selon la revendication 5, dans lequel, pour estimer une fréquence instantanée lissée du signal réverbéré pour chaque bande fréquentielle k de la 30 pluralité de N bandes fréquentielles, on calcule une transformée de Fourier locale discrète du signal acoustique 3051958 23 réverbéré en utilisant des fonctions fenêtres w(n) avec n entre 0 et N-1, et pour déterminer la phase instantanée de signal déréverbéré dans chaque bande fréquentielle k de ladite 5 pluralité de N bandes fréquentielles, on intègre la fréquence instantanée de signal déréverbéré au cours du temps en ajoutant un terme de correction fonction des fonctions fenêtres w(n). 7. Procédé selon la revendication 6, dans lequel, 10 pour estimer une fréquence instantanée lissée du signal réverbéré pour chaque bande fréquentielle k de la pluralité de N bandes fréquentielles, on applique un algorithme de vocodeur réassigné utilisant une transformée de Fourier locale discrète du signal acoustique réverbéré puis on applique un filtre de lissage temporelle. 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel au cours d'une étape de calibration, on mesure un signal acoustique de référence réverbéré par une propagation dans le milieu et on détermine le facteur d'influence du milieu à partir dudit signal acoustique de référence, en particulier on détermine un temps de réverbération dudit milieu. 9. Procédé pour reconstruire un signal déréverbéré dans lequel, on mesure un signal acoustique réverbéré par une propagation dans un milieu, on détermine un spectre d'amplitude de signal déréverbéré pour une pluralité de N bandes fréquentielles à partir du signal acoustique réverbéré, on estime une phase instantanée de signal déréverbéré pour chaque bande fréquentielle k de ladite pluralité de N bandes fréquentielles au moyen d'un procédé selon l'une quelconque des revendications 5 à 8, et 3051958 24 on reconstruit un signal déréverbéré à partir dudit spectre d'amplitude de signal déréverbéré et desdites phases instantanées de signal déréverbéré. 10. Dispositif pour estimer une phase instantanée 5 de signal acoustique déréverbéré, comprenant : (a)des moyens de mesure (2) pour capter au moins un signal acoustique réverbéré par une propagation dans un milieu, (b)des moyens (8) pour estimer au moins une 10 fréquence instantanée lissée du signal réverbéré ainsi qu'un taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré, (c)des moyens (8) pour calculer au moins une fréquence instantanée de signal déréverbéré à partir de 15 ladite fréquence instantanée du signal acoustique réverbéré, du taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré et d'un facteur d'influence du milieu, ledit facture d'influence étant fonction d'un temps de réverbération dudit milieu, 20 (d)des moyens (8) pour déterminer au moins une phase instantanée de signal déréverbéré en intégrant la fréquence instantanée de signal déréverbéré au cours du temps.