Brevet FR2865310A1 — publié le 20050722

PROCEDE DE RESTAURATION DE PARTIELS D'UN SIGNAL SONORE La présente invention se rapporte au domaine des télécommunications et en particulier au domaine du traitement numérique d'un signal sonore et à la représentation harmonique d'un tel signal. En modélisation harmonique des signaux audionumériques, le signal sonore est représenté par un ensemble d'oscillateurs dont les paramètres (fréquence, amplitude, phase) varient lentement au cours du temps. L'analyse harmonique comprend une analyse temps/fréquence à court terme qui permet de déterminer les valeurs de ces paramètres, suivie d'une extraction de pics puis d'un suivi de partiels. Le signal à modéliser est découpé en trames de.Q échantillons (typiquement 2 =1024). Un premier module d'analyse temps/fréquence à court terme (qui effectue typiquement une transformée de Fourier) permet de calculer le spectre du signal à court terme pour chaque trame. Un deuxième module d'extraction de pics permet de ne retenir que les pics les plus pertinents à priori, un critère étant par exemple de ne garder que les pics les plus énergétiques. Un troisième et dernier module cherche à lier les pics entre eux au cours du temps, c'est-à-dire d'une trame à l'autre pour former les partiels. Chaque partiel correspond pendant sa durée de vie à un oscillateur. Ce type d'analyse et de représentation peut être utilisé en particulier lors d'un codage à réduction de débit, lors d'un codage paramétrique (c'est-à-dire un codage qui traite le signal suivant trois aspects: transitoires, sinusoïdes, bruit), lors de la séparation et l'indexation de sources sonores et lors de la restauration de fichiers sonores. Il est couramment admis que la synthèse des partiels est de meilleure qualité en utilisant des techniques d'interpolation des phases proposées par Robert J.McAulay et Thomas F.Quatieri dans l'article "Speech Analysis/Synthesis Based on a Sinusoidal Representation", 1EEE Transaction on Acoustics, Speech and Signal Processing, PP 744-754, 1986 ou encore proposées par Laurent Girin, Sylvain Marchand, Joseph di Martino, axel Rôbel et Geoffroy Peeters dans l'article "Comparing the order of a Polynomial Phase Model for the Synthesis of Quasi-Harmonic Audio Signals", WASPAA, New Paltz, NY, USA, October 2003. Ces techniques permettent la synthèse d'un partiel d'un pic (A, f,., ço.) à un pic (A;+,, f + rp;+,) en calculant toutes les phases intermédiaires à l'aide de polynômes d'ordre 3 ou 5, les fréquences se déduisant par dérivation. Une interpolation d'ordre 3 est utilisée quand seules sont connues les fréquences et les phases de départ et d'arrivée. Une interpolation d'ordre 5 est utilisée quand sont, en outre, connues les variations d'ordre 2 de la phase (équivalentes aux variations à l'ordre 1 de la fréquence puisque par définition la fréquence est la dérivée de la phase). La synthèse d'un partiel entre les pics Pi (Ai, f,rpi) et P+l(Ai+,,f+I, çoi+,) consiste à calculer les valeurs p (n) du partiel entre les trames i et i+l: pi (n) = p(li + n) = Ai (n) cos (Mn , n = 0, e 1 (1) A cette fin, il est connu de calculer toutes les phases intermédiaires par l'une des deux méthodes d'interpolation suivantes. Pour l'interpolation d'ordre 3 selon Mac Aulay et al, la phase est calculée au moyen de l'expression suivante: roi(n)=roi+2DcfnTe+a(nTe)2+ fi(nTe)3 (2) Où Te est la période d'échantillonnage Les deux inconnues a et fi se calculent par résolution d'un système d'équations mettant en jeu (f, roi, f +, , roi+,) . Les fréquences se déduisent par dérivation 27c fi (n) = 21cf +2anTe+3f (nTe)2 (3) Pour l'interpolation d'ordre 5 selon Girin et al, les variations 8f et 8 f+, à l'ordre 1 de la fréquence aux pics P et sont supposées connues. La phase est alors calculée au moyen de l'expression suivante: roi (n)=çoi+27rfnTe+Bf (nTe)2+fi(nTe)3+ y(nTe)4+8(nTe)5 (4) Les trois inconnues)3, 8, y se calculent par résolution d'un système d'équations mettant enjeu (f, f +i, roi, roi+, , 8 f, 8 fi+,) . Les fréquences se déduisent par dérivation: 27r f (n) = 2c f + 8 fi nTe + 3 f (nTe)2 + 4y (nTe)3 + 58 (nTe)4 (5) Pour des raisons diverses, il peut arriver que certains partiels existants dans le signal soient absents, corrompus ou discontinus en sortie d'analyse et/ou en entrée de synthèse. Par exemple, ils peuvent être absents en entrée du décodeur dans une application de diffusion de programmes sonores sur Internet en cas de perte de paquets, ils peuvent être corrompus en cas de perturbations du signal à analyser par un signal parasite (bruit, clic, autre signal, etc), ou ils peuvent être discontinus dans le cas où ils ont une énergie trop faible pour être détectés correctement de façon continue. Il s'avère alors nécessaire de mettre en uvre des techniques de restauration des pics manquants pour pouvoir recréer un signal synthétisé le plus proche possible du signal original. Ceci nécessite de recréer des pics caractérisés chacun par une amplitude, une fréquence et une phase. Les techniques précédentes d'interpolation, connues de l'art antérieur, sont utilisées pour synthétiser les parties correspondant aux pics manquants et restaurer les partiels. Cependant ces techniques d'interpolation connues sont adaptées au court terme, c'est-à-dire pour une période inférieure à 10ms. Pour des durées plus longues, le signal re-synthétisé est souvent éloigné de l'original et des artefacts désagréables peuvent apparaître. En effet, ces techniques assurent une continuité de phase entre les pics existants et les pics restaurés, mais en contre partie elles ne permettent pas de contrôler les fréquences induites données par les équations (3) et (5). Cet effet est d'autant plus marqué que la distance d'interpolation est grande. Un but de l'invention est de proposer une solution alternative au problème de la restauration de la partie manquante et identifiée comme telle d'un partiel, notamment lorsque la partie manquante correspond à des temps longs (supérieurs à 10 ms) pour lesquels les techniques connues sont peu efficaces. Aussi, le problème technique à résoudre par l'objet de la présente invention est de proposer un procédé de restauration de parties manquantes des partiels d'un signal sonore, lors d'une analyse harmonique suivant laquelle le signal sonore est découpé en trames temporelles sur lesquelles est appliquée une analyse temps/fréquence qui fournit des spectres successifs à court terme représentés par des trames fréquentielles d'échantillons, l'analyse consistant en outre à extraire des pics spectraux dans les trames fréquentielles et à les lier entre eux au cours du temps pour former des partiels, ce procédé étant une alternative aux solutions connues. Une solution au problème technique posé consiste, selon la présente invention, en ce que ledit procédé de restauration d'un partiel entre un pic P et un pic P,.+N dont les fréquences co et phases cl) sont connues est tel qu'il comprend les étapes qui consistent: - à estimer la fréquence ci de chacun des pics manquants l';+1 à P+N_1 de ce partiel, à calculer la phase déroulée cp de pic en pic, depuis la phase du pic P jusqu'à celle du pic P+N et pour toutes les fréquences di précédemment estimées, - à calculer l'erreur de phase. errcp entre la phase déroulée cp et la phase connue au même pic P+N, - à corriger chaque phase déroulée cp d'une valeur fonction de l'erreur de phase errcp. Un procédé selon l'invention se différencie des méthodes connues en ce qu'il effectue un contrôle plus fin de la fréquence des pics manquants et un calcul après coup des phases correspondantes pour assurer la continuité avec les phases des pics existants. Ainsi, un procédé selon l'invention permet une re-synthèse sans artefacts des signaux correspondants aux morceaux de partiels manquants contrairement aux méthodes connues préalablement décrites. En outre, de manière avantageuse, un procédé selon l'invention permet une reconstruction du signal plus proche au sens de l'erreur de reconstruction du signal original que celle obtenue par les méthodes connues. Enfin, un procédé selon l'invention présente avantageusement un algorithme à faible complexité. L'invention a en outre pour objet un dispositif de synthèse d'un signal sonore pour la mise en oeuvre d'un procédé de restauration d'un partiel entre un pic P et un pic P+N. Ce dispositif est par exemple un décodeur audio ou un codeur paramétrique adapté pour la mise en oeuvre d'un procédé selon l'invention. L'invention a en outre pour objet un produit programme d'ordinateur chargeable directement dans la mémoire interne d'un dispositif précédent, ou d'un groupe de dispositifs. Ce produit programme d'ordinateur comprend des portions de code logiciel pour l'exécution des étapes d'un procédé selon l'invention lorsque le programme est exécuté sur le dispositif ou le groupe de dispositifs. L'invention a en outre pour objet un support utilisable dans un dispositif précédent ou un groupement de dispositifs et sur lequel est enregistré un produit pro gamme d'ordinateur chargeable direectemeent dans la mémoire interne du dispositif âG1.111111V ou du groupe de dispositifs, comprenant des portions de code logiciel pour l'exécution des étapes d'un procédé selon l'invention, lorsque le programme est exécuté sur le dispositif ou le groupe de dispositifs. D'autres caractéristiques et avantages de l'invention apparaîtront lors de la description qui suit faite en regard de figures annexées données à titre d'exemples non limitatifs. La figure 1 est un organigramme d'un exemple de déroulement d'un procédé selon l'invention. La figure 2 est un schéma d'un exemple d'utilisation d'un procédé selon l'invention. Un procédé selon l'invention se déroule de la façon suivante décrite en regard de l'organigramme de la figure 1. Le procédé 1 consiste à restaurer un partiel entre un pic P et un pic P+N dont les fréquences ça et phases (p sont connues. Dans une première étape 2, le procédé estime la fréquence cô et l'amplitude A de chacun des pics manquants P+1 à P,+N_1. Cette estimation est effectuée par exemple par interpolation ou prédiction linéaire selon des méthodes connues. Soit un partiel constitué d'une succession de pics P(A,,c);,cp,) liés entre eux, connus à des instants iT et caractérisés par: A, l'amplitude du pic au temps iT ca, , la fréquence du pic au temps iT (p. la phase du pic au temps iT, donnée modulo 2n. L'estimation de la fréquence des pics manquants entre les pics P et P+N est effectuée par exemple au moyen d'une interpolation linéaire entre co, et w,+N, ou au moyen d'une prédiction linéaire sur le passé ou sur le futur, décrite par exemple dans l'article "Enhanced Partial Tracking using linear Prediction", de Mathieu Lagrange, Sylvain Marchand, martin Raspaud et Jean-Bernard Rault, Proceedings of the Digital Audio Effects (DAFx) Conference, pp141-146, Queen Mary, University of London, UK, September 2003, ou encore au moyen d'une combinaison pondérée sur le passé ou sur le futur. L'estimation de l'amplitude A des pics manquants est effectuée par exemple au moyen d'une interpolation linéaire entre A, et A,+N, ou au moyen d'une prédiction linéaire sur le passé ou sur le futur ou encore au moyen d'une combinaison pondérée sur le passé ou sur le futur. Dans une deuxième étape 3, le procédé calcule la phase déroulée cp de pic en pic, depuis la phase du pic P jusqu'à celle du pic P+N. Ce calcul est effectué pour chacune des fréquences co précédemment estimées. Soient cf), et a i la phase et la fréquence de départ et {eÎi+1,...,Cili+ N-1} une estimation des fréquences dans l'intervalle à reconstruire. Pour prolonger le partiel entre le pic P et le pic P+N, le procédé déroule la phase selon l'expression suivante: n rPi+n=mod(cpi+E i+j+CAi+j-1 T, 27r n=1,...,N j=1 2 (6) Pour ne pas engendrer de discontinuités nuisibles à la qualité de la re-synthèse, il faut obtenir à l'instant i+N une phase reconstruite cpi+N égale à (pi+N. Les données intervenant dans l'expression (6) précédente étant soient approximées, soient prédites, il est statistiquement impossible d'obtenir cette égalité. Par conséquent, le procédé répartit l'erreur de phase erre) calculée à l'instant i+N entre tous les pics manquants et préalablement reconstruits P+1 à P+N_1 au moyen des étapes suivantes. Dans une troisième étape 4, le procédé calcule l'erreur de phase erre) entre la phase déroulée qPi+N et la phase connue cpi+N au même pic P+N. Ce calcul peut être effectué selon le système d'équations suivant: erre) =(pi+N +27t si Irpi+N (Pi+N + 27t I Dans une quatrième étape 5, le procédé corrige chaque phase déroulée. q i+ n d'une valeur fonction de l'erreur de phase erre) . Typiquement l'erreur de phase calculée à l'instant i+N est répartie de manière uniforme sur chacune des phases déroulées selon l'expression suivante: cp corrigéei+n = mod l cpi+n + errcp N, 27ti] n =1, ..., N -1 (10) La répartition peut ne pas être uniforme et suivre une loi non linéaire par exemple. L'exemple d'utilisation illustré par la figure 2, consiste à restaurer les partiels au moyen d'un procédé 1 selon l'invention lors d'une analyse harmonique d'un signal sonore par exemple lors d'un codage paramétrique. Le signal sonore s(n) est représenté par un ensemble d'oscillateurs dont les paramètres (fréquence, amplitude) varient lentement au cours du temps. De manière classique, l'analyse harmonique comprend une analyse 6 temps/fréquence à court terme qui permet de déterminer les valeurs de ces paramètres, suivie d'une extraction 7 de pics puis d'un suivi 8 de partiels. Une détection 9 de trous dans les partiels précède la mise en oeuvre d'un procédé 1 de restauration des partiels selon l'invention. Les pics reconstruits P,.+n (A;+n i+n ç +n) lors de la mise en oeuvre du procédé 1, sont ensuite traités comme des pics issus de l'analyse harmonique et la synthèse 10 additive du signal correspondant au partiel restauré à partir de ces pics reconstruits peut se faire, par exemple, par une des méthodes connues d'interpolation des phases (ordre 3 ou ordre 5). REVENDICATIONS, 1. Procédé (1) de restauration de partiels d'un signal sonore, lors d'une analyse harmonique suivant laquelle le signal sonore est découpé en trames temporelles sur lesquelles est appliquée une analyse temps/fréquence qui fournit des spectres successifs à court terme représentés par des trames fréquentielles d'échantillons, l'analyse consistant en outre à extraire des pics spectraux dans les trames fréquentielles et à les lier entre eux au cours du temps pour former des partiels, le procédé de restauration d'un partiel entre un pic P et un pic P+N dont les fréquences et phases sont connues est caractérisé en ce qu'il comprend les étapes qui consistent: à estimer (2) la fréquence éà de chacun des pics manquants P+, à P+N_, de ce partiel, 2. Procédé (1) de restauration de partiels d'un signal sonore selon la revendication 1, dans lequel la phase déroulée cp est calculée par la formule: +a)" ; ,+n = mod ç1 + ,+; I+ 1 T, 2rc, n =1,..., N j_l 2 avec cp; et w, = wi la phase et la fréquence du pic P, çp;+N et 6i +N = w;+N la phase 25 et la fréquence du pic P+N 3. Procédé (1) de restauration de partiels d'un signal sonore selon l'une des revendications 1 et 2, dans lequel l'estimation de la fréquence Co des pics manquants P+, à est effectuée au moyen d'une interpolation linéaire entre les fréquences des pics connus P,. et P+N. - à calculer (3) la phase déroulée P de pic en pic, depuis la phase du pic P jusqu'à celle du pic P,.+N et ce, pour toutes les fréquences Co précédemment estimées, - à calculer (4) l'erreur de phase errço entre la phase déroulée et la phase connue au même pic P+N, - à corriger (5) chaque phase déroulée d'une valeur fonction de l'erreur de phase errço. 4. Procédé (1) de restauration de partiels d'un signal sonore selon l'une des revendications 1 et 2, dans lequel l'estimation de la fréquence w des pics manquants P+, à P+,_, est effectuée au moyen d'une prédiction linéaire sur le passé. 5. Procédé (1) de restauration de partiels d'un signal sonore selon l'une des revendications 1 et 2, dans lequel l'estimation de la fréquence tv des pics manquants P+, à P+N_, est effectuée au moyen d'une prédiction linéaire sur le futur. 6. Procédé (1) de restauration de partiels d'un signal sonore selon l'une des revendications 1 et 2, dans lequel l'estimation de la fréquence éo' des pics manquants P+, à P+,_, est effectuée au moyen d'une combinaison pondérée d'une prédiction linéaire sur le passé et d'une prédiction linéaire sur le futur. 7. Procédé (1) de restauration de partiels d'un signal sonore selon l'une des revendications précédentes, qui comprend en outre l'étape qui consiste: - à estimer l'amplitude de chacun des pics manquants P+, à de ce partiel au moyen d'une interpolation linéaire entre les amplitudes A des pics connus Pi et P+N 8. Procédé (1) de restauration de partiels d'un signal sonore selon l'une des revendications 1 à 6, qui comprend en outre l'étape qui consiste: - à estimer l'amplitude de chacun des pics manquants P+, à P+N-I de ce partiel 25 au moyen d'une prédiction linéaire sur le passé. 9. Procédé (1) de restauration de partiels d'un signal sonore selon l'une des revendications 1 à 6, qui comprend en outre l'étape qui consiste: à estimer l'amplitude de chacun des pics manquants P+, à P+N_, de ce partiel 30 au moyen d'une prédiction linéaire sur le futur. 10. Procédé (1) de restauration de partiels d'un signal sonore selon l'une des revendications 1 à 6, qui comprend en outre l'étape qui consiste: - à estimer l'amplitude de chacun des pics manquants à P+,_, de ce partiel au moyen d'une prédiction linéaire sur le passé et d'une prédiction linéaire sur le futur. 11. Procédé (1) de restauration de partiels d'un signal sonore selon l'une des revendications précédentes, dans lequel la correction de phase consiste à répartir uniformément l'erreur de phase errço calculée à l'instant i+N entre tous les pics manquants P+, à P+N_, du partiel. 12. Procédé (1) de restauration de partiels d'un signal sonore selon la revendication précédente, dans lequel la phase corrigée est déterminée par l'équation: rpcorrigéei+n = mod tpi+n + erre N,27r n =1,..., N -1 13. Procédé (1) de restauration de partiels d'un signal sonore selon la revendication précédente, dans lequel l'erreur de phase err0 est déterminée par le système d'équations: errço = çoi+N Oi+N + 27r si I çoi+ N Oi+N + 27rI 14. Dispositif de synthèse d'un signal sonore pour la mise en oeuvre d'un procédé selon l'une des revendications précédentes, caractérisé en ce qu'il comprend: - des moyens d'estimation de la fréquence w de chacun des pics manquants P+, à P+N_, de ce partiel, des moyens de calcul de la phase déroulée 0 de pic en pic, depuis la phase du pic P,. jusqu'à celle du pic P+N et ce, pour toutes les fréquences w précédemment estimées, - des moyens de calcul de l'erreur de phase err0 entre la phase déroulée 0 et la phase connue au même pic P+N, - des moyens de correction de chaque phase déroulée b d'une valeur fonction de l'erreur de phase errrp. 15. Un produit programme d'ordinateur chargeable directement dans la mémoire interne d'un dispositif ou d'un groupe de dispositifs selon la revendication précédente, comprenant des portions de code logiciel pour l'exécution des étapes d'un procédé (1) selon l'une des revendications 1 à 13, lorsque le programme est exécuté sur le dispositif ou le groupe de dispositifs. 16. Support utilisable dans un dispositif ou un groupement de dispositifs selon la revendication 14 et sur lequel est enregistré un produit programme d'ordinateur chargeable directement dans la mémoire interne du dispositif ou du groupe de dispositifs, comprenant des portions de code logiciel pour l'exécution des étapes d'un procédé (1) selon l'une des revendications 1 à 13, lorsque le programme est exécuté sur le dispositif ou le groupe de dispositifs.