FR2882458A1 - Procede de mesure de la gene due au bruit dans un signal audio - Google Patents
Procede de mesure de la gene due au bruit dans un signal audio Download PDFInfo
- Publication number
- FR2882458A1 FR2882458A1 FR0501747A FR0501747A FR2882458A1 FR 2882458 A1 FR2882458 A1 FR 2882458A1 FR 0501747 A FR0501747 A FR 0501747A FR 0501747 A FR0501747 A FR 0501747A FR 2882458 A1 FR2882458 A1 FR 2882458A1
- Authority
- FR
- France
- Prior art keywords
- noise
- signal
- frame
- nob
- coefficients
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 76
- 230000005236 sound signal Effects 0.000 title claims abstract description 48
- 108090000623 proteins and genes Proteins 0.000 title 1
- 238000012360 testing method Methods 0.000 claims abstract description 67
- 230000009467 reduction Effects 0.000 claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims description 37
- 230000003595 spectral effect Effects 0.000 claims description 23
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 230000007480 spreading Effects 0.000 claims description 4
- 230000006870 function Effects 0.000 description 27
- 230000000875 corresponding effect Effects 0.000 description 8
- 238000011156 evaluation Methods 0.000 description 8
- 230000008901 benefit Effects 0.000 description 6
- 238000001514 detection method Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 description 2
- 230000035807 sensation Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/48—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
- G10L25/69—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for evaluating synthetic or decoded voice signals
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Quality & Reliability (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Noise Elimination (AREA)
- Signal Processing Not Specific To The Method Of Recording And Reproducing (AREA)
Abstract
L'invention concerne un procédé de calcul d'une note objective (NOB) de la gêne due au bruit dans un signal audio traité par une fonction de réduction de bruit, ledit procédé comprenant une étape préalable d'obtention d'un signal audio prédéfini de test (x[m]) contenant un signal utile dépourvu de bruit, d'un signal bruité (xb[m]), obtenu en additionnant un signal de bruit prédéfini audit signal de test (x[m]), et d'un signal traité (y[m]), obtenu par application de la fonction de réduction de bruit audit signal bruité (xb[m]), ledit procédé étant caractérisé en ce qu'il inclut une étape (a3, a4) de mesures de sonie de trames dudit signal bruité (xb[m]) et dudit signal traité (y[m]), et de mesures de coefficients de tonalité de trames dudit signal traité (y[m]).
Description
Procédé de mesure de la gêne due au bruit dans un signal audio
La présente invention se situe de manière générale dans les domaines du traitement du signal de parole et de la psychoacoustique. Plus précisément l'invention concerne un procédé et un dispositif d'évaluation objective de la gêne due au bruit dans des signaux audio.
L'invention permet notamment de noter objectivement la gêne due au bruit dans un signal audio traité par une fonction de réduction de bruit.
Dans le domaine de la transmission de signaux audio, une fonction de réduction de bruit, aussi appelée fonction de suppression de bruit ou de débruitage, a pour objectif de réduire le niveau de bruit de fond dans une communication vocale, ou ayant au moins une composante vocale. Elle présente un intérêt spécifique lorsque l'un des interlocuteurs de cette communication est immergé dans un milieu bruité qui nuit fortement à l'intelligibilité de sa voix. Les algorithmes de réduction de bruit sont basés sur une estimation en continu du niveau du bruit de fond à partir du signal incident et d'une détection d'activité vocale permettant de distinguer les périodes de bruit seul de celles avec du signal de parole utile. Un filtrage du signal de parole incident, correspondant au signal de parole bruité, est ensuite effectué de façon à réduire la contribution du bruit déterminée à partir de l'estimée du bruit.
La gêne due à la présence de bruit dans un signal audio traité par une telle fonction de réduction de bruit est évaluée aujourd'hui de manière subjective seulement en se basant sur l'exploitation de résultats de tests mis en oeuvre selon le document "Recommandation UIT-T P.835 (11/2003)". Cette évaluation est faite sur une échelle de type MOS, d'après l'anglais Mean Opinion Score, qui donne une note de un à cinq de la gêne due au bruit, appelée "background noise" dans ce même document.
L'inconvénient majeur de cette technique d'évaluation est la nécessité de mettre en oeuvre des tests subjectifs, cette mise en oeuvre étant très lourde et très coûteuse. En effet chaque contexte particulier, c'est-àdire un type de signal incident associé à un type de bruit et une fonction de réduction de bruit, nécessite de mettre un panel de personnes en situation d'écoute réelle d'échantillons de parole afin de leur demander de noter la gêne due au bruit selon une échelle de type MOS.
C'est pourquoi le développement de méthodes objectives alternatives pouvant compléter ou suppléer les méthodes subjectives est un sujet de grand intérêt. L'illustration la plus frappante de ce phénomène est le modèle de qualité d'écoute, en constante évolution, contenu dans le document "Recommandation UIT-T P.862 (02/2001)". Néanmoins ce modèle ne s'applique pas à l'évaluation de la gêne due au bruit. L'invention concerne en effet des signaux de parole dans lesquels la gêne due au bruit peut être importante, ceci avant ou après traitement de ces signaux par une éventuelle fonction de réduction de bruit.
II est de plus à noter que bien qu'en général l'invention sera utilisée pour évaluer la gêne due au bruit en sortie d'un équipement de communication implémentant une fonction de réduction de bruit, l'invention s'applique aussi aux signaux bruités non traités par une telle fonction. Le cas d'utilisation de l'invention sur un signal audio bruité quelconque est donc un cas particulier du cas plus général d'utilisation de l'invention sur un signal audio traité par une fonction de réduction de bruit.
La présente invention a pour but de résoudre les inconvénients de la technique antérieure en fournissant un procédé et un dispositif de calcul objectif d'une note équivalente à la note subjective telle qu'indiquée dans le document "Recommandation UIT-T P.835", caractérisant la gêne due à la présence de bruit dans un signal audio. Le procédé selon l'invention varie suivant que l'invention est utilisé sur un signal audio bruité quelconque ou sur un signal audio traité par une fonction de réduction de bruit, notamment dans les paramètres de calcul de la note objective selon l'invention. Afin de bien décrire ces deux cas d'utilisation, deux modes de réalisation pouvant aussi être considérés comme deux procédés distincts sont présentés. Cependant le second mode de réalisation, s'appliquant à un signal auclio bruité quelconque, et plus général que le premier mode de réalisation, se déduit aisément de celui-ci.
A cette fin, l'invention propose un procédé de calcul d'une note objective de la gêne due au bruit dans un signal audio traité par une fonction de réduction de bruit, ledit procédé comprenant une étape préalable d'obtention d'un signal audio prédéfini de test contenant un signal utile dépourvu de bruit, d'un signal bruité, obtenu en additionnant un signal de bruit prédéfini audit signal de test, et d'un signal traité, obtenu par application de la fonction de réduction de bruit audit signal bruité, ledit procédé étant caractérisé en ce qu'il inclut une étape de mesures de sonie de trames dudit signal bruité et dudit signal traité, et de mesures de coefficients de tonalité de trames dudit signal traité.
Ce procédé a l'avantage d'une mise en oeuvre simple, immédiate et rapide contrairement aux tests subjectifs. On rappellera ici que l'expression "sonie psychoacoustique" peut être définie comme le caractère de la sensation auditive lié à la pression acoustique et à la structure du son. En d'autres termes, il s'agit de la force sonore d'un son ou d'un bruit en tant que sensation auditive (cf. Office de la langue française, 1988). La sonie est représentée par une échelle de sonie psychoacoustique (en sones). D'autre part, la densité de sonie, encore désignée par "intensité subjective", est une mesure particulière de la sonie.
Selon une caractéristique préférée, ce procédé selon l'invention comprend les étapes de: - Calcul de densités de sonie moyenne -s-y(m)de trames du signal traité, de densités de sonie moyenne respectives Sxe(m parole) et Sr(m _ parole) de trames de signal utile "m_parole" respectivement du signal bruité et du signal traité, de densités de sonie moyenne SY(m_bruit)de trames de bruit "m_bruit" du signal traité, et de coefficients de tonalité aY(m bruit)de trames de bruit "m_bruit" du signal traité, - Calcul d'une note objective de la gêne due au bruit dans le signal traité, à partir desdites densités de sonie moyenne et desdits coefficients de tonalité calculés, et de coefficients de pondération prédéfinis. Selon une caractéristique préférée, l'étape de calcul de densités de sonie moyenne et de coefficients de tonalité est suivie d'une étape de calcul des moyennes SY, SXb parole, SY _ parole, SY _bruit et aY _ brui.' desdites densités de sonie moyenne et desdits coefficients de tonalité sur l'ensemble des trames concernées des signaux correspondants, et la note objective de la gêne due au bruit est calculée selon l'équation suivante: s NOB = cofacteur(i) + (06 r= où facteur(1)= SYbruit
SY
facteur(2)= SY_bruit SY _ parole facteur(3)= Ecart_type (Sxb(m _ parole) SY(m parole)) , l'opérateur "Ecart_type (v(m))" désignant l'écart- type de la variable v sur l'ensemble des trames d'indice m, facteur(4)= aY bruit facteur(5)= Ecart_type (aY(m _ bruit)), et les coefficients col à 06 sont déterminés de manière à obtenir une corrélation maximale entre les données subjectives issues d'une base de données de tests subjectifs et les notes objectives calculées par ledit procédé pour les signaux de tests, bruités et traités correspondants utilisés lors desdits tests subjectifs.
Les coefficients de cette combinaison linéaire ont l'avantage de pouvoir être recalculés si de nouvelles données de tests subjectifs modifient de manière sensible la corrélation précédemment éi:ablie. Ceci permet d'améliorer un modèle objectif alimenté par le procédé selon l'invention, de calcul de la gêne due au bruit dans un signal audio traité par une fonction de réduction de bruit, par une simple reconfiguration des paramètres du procédé.
L'invention concerne aussi un procédé de calcul d'une note objective de la gêne due au bruit dans un signal audio, ledit procédé comprenant une étape préalable d'obtention d'un signal audio prédéfini de test contenant un signal utile dépourvu de bruit, et d'un signal bruité, obtenu en additionnant un signal de bruit prédéfini audit signal de test, ledit procédé étant caractérisé en ce qu'il inclut une étape de mesures de sonie et de mesures de coefficients de tonalité de trames dudit signal bruité.
Ce procédé a les mêmes avantages que le procédé précédent, mais s'applique à un signal audio bruité quelconque.
Selon une caractéristique préférée, ce procédé selon l'invention comporte les étapes de: - Calcul de densités de sonie moyenne Svv(m) de trames du signal bruité, de densités de sonie moyenne Sxb(m _ parole) de trames de signal utile "m_parole" du signal bruité, de densités de sonie moyenne Sxb(m _ bruit) de trames de bruit "m_bruit" du signal bruité, et de coefficients de tonalité axb(m _ bruit) de trames de bruit "m_bruit" du signal bruité, Calcul d'une note objective de la gêne due au bruit dans le signal bruité, à partir desdites densités de sonie rnoyennes et desdits coefficients de tonalité calculés, et de coefficients de pondération prédéfinis.
Selon une caractéristique préférée, l'étape de calcul de densités de sonie moyenne et de coefficients de tonalité est suivie d'une étape de calcul des moyennes Sxb,-S- parole, Sxb _ bruit et axb _ bruit desdites densités de sonie moyenne et desdits coefficients de tonalité sur l'ensemble des trames concernées des signaux correspondants, et en ce que ladite note objective de la gêne due au bruit est calculée selon l'équation suivante: NOB = E wfacteur(i) + cos où i=1 facteur(1)= Sxb= brui Sxb facteur(2)= Sxb_brut Sxb _ parole facteur(3)= axb bruit, facteur(4)= Ecart_type(axb(m_bruit)), l'opérateur "Ecart_type (v(m))" désignant l'écart-type de la variable v sur l'ensemble des trames d'indice m, et les coefficients coi à w5 sont déterminés de manière à obtenir une corrélation maximale entre les données subjectives issues d'une base de données de tests subjectifs et les notes objectives calculées par ledit procédé pour les signaux de tests et les signaux bruités correspondants utilisés lors desdits tests subjectifs.
Comme pour le procédé précédent, les coefficients de cette combinaison linéaire ont l'avantage de pouvoir être recalculés si de nouvelles données de tests subjectifs modifient de manière sensible la corrélation précédemment établie. Ceci permet d'améliorer un modèle objectif alimenté par le procédé selon l'invention, de calcul de la gêne due au bruit dans un signal audio, par une simple reconfiguration des paramètres du procédé.
Selon une caractéristique préférée de ces deux procédés selon l'invention, l'étape de calcul de densités de sonie et de coefficients de tonalité est précédée d'une étape de détection d'activité vocale sur le signal de test, de manière à déterminer si une trame courante du signal bruité, et du signal traité dans le cas du premier procédé, est une trame "m_bruit" contenant seulement du bruit, ou une trame "m_parole" contenant de la parole, dite trame de signal utile.
Cette étape de détection d'activité vocale permet de séparer très simplement les différents types de trames du signal bruité, et du signal traité dans le cas du premier procédé, par l'utilisation du signal de test.
Selon une caractéristique préférée de ces deux procédés selon l'invention, l'étape de calcul de la note objective est suivie d'une étape de calcul d'une note objective sur l'échelle MOS de la gêne due au bruit, calculée selon l'équation suivante: NOB _ MOS = 2; (NOB)-' , dans laquelle les coefficients XI à 24 sont déterminés de manière à ce que ladite nouvelle note objective obtenue caractérise la gêne due au bruit sur l'échelle MOS.
Le fait d'utiliser une fonction polynomiale d'ordre 3 permet d'obtenir une note objective sur l'échelle MOS très proche de la note subjective MOS que donnerait un groupe d'auditeurs dans le cadre d'un test subjectif conforme à la "Recommandation UIT-T P.835".
Selon une caractéristique préférée de ces deux procédés selon l'invention, l'étape de calcul de densités de sonie et de coefficients de tonalité, le calcul de la densité de sonie moyenne Su (m) d'une trame d'indice m quelconque d'un signal audio donné u, comprend les étapes suivantes: fenêtrage, par exemple de type Hanning, de la trame d'indice m et obtention d'une trame fenêtrée u_w[m], - application d'une transformée de Fourier rapide à la trame fenêtrée u_w[m] et obtention d'une trame correspondante U(m,f) dans le domaine fréquentiel, calcul de la densité spectrale de puissance yu (m, f) de la trame U(m,f), -application à la densité spectrale de puissance yu(m,f) d'une conversion de l'axe des fréquences à l'échelle des Barks et obtention d'une densité spectrale de puissance Bu(m,b) sur l'échelle des Barks, - convolution de la densité spectrale de puissance sur l'échelle des Barks, Bu(m,b), avec la fonction d'étalement couramment utilisée en psychoacoustique et obtention d'une densité spectrale étalée sur l'échelle des Barks, Eu(m,b), calibration de la densité spectrale étalée sur l'échelle des Barks, Eu(m, b), par les facteurs respectifs d'échelonnement en puissance et d'échelonnement en sonie couramment utilisés en psychoacoustique, conversion de la grandeur ainsi obtenue sur l'échelle des phones puis conversion sur l'échelle des sones de la grandeur précédemment convertie en phones, et obtention en conséquence d'un nombre B de valeurs de densité de sonie, Su(m,b), de la trarne d'indice m pour la bande critique b, B étant le nombre de bandes critiques considérées dans l'échelle des Barks et l'indice b variant de 1 à B, calcul de la densité de sonie moyenne Su (m) de la trame d'indice m à partir desdites B valeurs de densités de sonie Su (m,b), selon l'équation suivante:
B
Su(m)= E S u (m, b) B b= Selon une caractéristique préférée de ces deux procédés selon l'invention, dans l'étape de calcul de densités de sonie et de coefficients de tonalité, le calcul du coefficient de tonalité a(m) d'une trame d'indice m quelconque d'un signal audio donné u, comprend les étapes suivantes: - fenêtrage, par exemple de type Hanning, de la trame d'indice m et obtention d'une trame fenêtrée u_w[m], -application d'une transformée de Fourier rapide à la trame fenêtrée u_w[m] et obtention d'une trame correspondante U(m,f) dans le domaine fréquentiel, - calcul de la densité spectrale de puissance yu (m, f) de la trame U(m,f), - calcul du coefficient de tonalité a(m) selon l'équation suivante: / ]ive (N-1 f ru(m,f) \I=o N-1 I Yu (m, f) N f=o a(m) _ 60 où * symbolise l'opérateur de multiplication dans l'espace des nombres réels, f représente l'indice fréquentiel de la densité spectrale de puissance, et N désigne la taille de la transformée de Fourier rapide.
L'invention concerne également un équipement de test destiné à évaluer une note objective de la gêne due au bruit dans un signal audio, caractérisé en ce qu'il comporte des moyens adaptés à mettre en oeuvre l'un ou l'autre des procédés selon l'invention.
Selon une caractéristique préférée, l'équipement de test inclut des moyens informatiques et un programme d'ordinateur, ledit programme comportant des instructions adaptées à mettre en oeuvre l'un ou l'autre desdits procédés, lorsqu'il est exécuté par lesdits moyens informatiques.
L'invention concerne encore un programme d'ordinateur sur un support d'informations, comportant des instructions adaptées à a mise en oeuvre de l'un ou l'autre des procédés selon l'invention, lorsque le programme est chargé et exécuté dans un système informatique.
Les avantages de cet équipement de test ou de ce programme d'ordinateur sont identiques à ceux mentionnés plus haut en relation avec les procédés de l'invention.
10*1og10 D'autres caractéristiques et avantages apparaîtront à la lecture de modes de réalisation préférés décrits en référence aux figures dans lesquelles: - la figure 1 représente un environnement de test destiné à calculer une note objective de la gêne due au bruit dans un signal audio traité par une fonction de réduction de bruit, selon un premier mode de réalisation de l'invention, - la figure 2 est un organigramme illustrant un procédé de calcul d'une note objective de la gêne due au bruit dans un signal audio traité par une fonction de réduction de bruit selon un premier mode de réalisation du procédé selon l'invention, - la figure 3 est un organigramme illustrant un procédé de calcul d'une note objective de la gêne due au bruit dans un signal audio selon un second mode de réalisation du procédé selon l'invention, - la figure 4 est un organigramme illustrant le mode de calcul de la densité de sonie moyenne et du coefficient de tonalité d'une trame de signal audio selon l'invention.
Deux modes de réalisation du procédé selon l'invention sont décrits dans la suite, le premier étant appliqué à un signal audio traité par une fonction de réduction de bruit, et le second étant appliqué à un signal audio bruité quelconque. Le principe du procédé selon l'invention est le même dans ces deux modes de réalisation, en particulier le procédé de calcul est exactement le même, mais dans le second mode de réalisation le signal audio traité par une fonction de réduction de bruit est pris égal au signal bruité. Le second mode de réalisation peut en effet être considéré comme un cas particulier du premier mode de réalisation, avec une fonction de réduction de bruit inhibée.
Selon le premier mode de réalisation du procédé l'invention, la gêne due à la présence de bruit dans un signal audio traité par une fonction de réduction de bruit est évaluée de manière objective dans un environnement de test représenté à la figure 1. Un tel environnement de test comprend une source de signaux audio SSA délivrant un signal audio de test x(n) ne contenant que du signal utile, c'est-à-dire dépourvu de bruit, par exemple un signal de parole, et une source de bruit SB délivrant un signal de bruit prédéfini.
Aux fins de test, ce signal de bruit prédéfini est ajouté au signal de test x(n) choisi, comme représenté par l'opérateur d'addition AD. Le signal audio résultant de cette addition de bruit au signal de test x(n) est noté xb(n) et est désigné par l'expression "signal bruité".
Le signal bruité xb(n) constitue alors le signal d'entrée d'un module MRB de réduction de bruit mettant en oeuvre une fonction de réduction de bruit délivrant en sortie un signal audio y(n) désigné par l'expression "signal traité". Le signal traité y(n) est donc un signal audio contenant du signal utile et un bruit résiduel.
Le signal traité y(n) est ensuite délivré à un équipement de test EQT mettant en oeuvre un procédé d'évaluation objective de la gêne due au bruit dans le signal traité, selon l'invention. Typiquement le procédé selon l'invention est implémenté dans l'équipement de test EQT sous la forme d'un programme d'ordinateur. En plus ou en remplacement de moyens logiciels, l'équipement de test EQT comporte éventuellement des moyens matériels électroniques pour implémenter le procédé selon l'invention. Outre le signal y(n), l'équipement de test EQT reçoit en entrée le signal de test x(n) et le signal bruité xb(n).
L'équipement de test EQT délivre en sortie un résultat d'évaluation RES, qui est une note objective NOB_MOS de la gêne due à la présence de bruit dans le signal traité y(n). Le mode de calcul de cette note objective NOB MOS sera décrit plus bas.
Les signaux audio précités x(n), xb(n) et y(n) sont des signaux échantillonnés dans un format numérique, n désignant un échantillon quelconque. Ces signaux sont par exemple supposés échantillonnés à la fréquence d'échantillonnage de 8 kHz (kilo Hertz).
Dans le mode de réalisation décrit et représenté ici, le signal de test x(n) est un signal de parole dépourvu de bruit. Le signal bruité xb(n) représente alors le signal vocal initial x(n) dégradé par un environnement bruité (bruit de fond ou bruit ambiant), et le signal y(n) représente le signal xb(n) après réduction de bruit.
Selon un exemple de mise en oeuvre de l'invention, le signal x(n) est généré dans une chambre anéchoïque. Cependant, le signal x(n) peut être aussi généré dans une pièce "calme" ayant un temps de réverbération "moyen", inférieur à 0,5 seconde.
Le signal bruité xb(n) est obtenu en ajoutant une contribution prédéterminée de bruit au signal x(n). Le signal y(n) est obtenu soit en sortie d'un algorithme de réduction de bruit implanté sur un ordinateur personnel, soit à la sortie d'un équipement réseau réducteur de bruit et dans ce dernier cas, le signal y(n) est prélevé au niveau d'un codeur MIC (modulation par impulsion et codage).
En référence à la figure 2, le procédé de calcul de la note objective NOB MOS de la gêne due au bruit dans le signal traité y(n) selon l'invention est représenté sous la forme d'un algorithme comportant des étapes al à a7.
Dans une première étape a1, les signaux x(n), xb(n) et y(n) sont respectivement découpés en fenêtres temporelles successives appelées trames. Chaque trame de signal, notée m, contient un nombre prédéterminé d'échantillons du signal, l'étape al consiste donc en un changement de cadence de chacun de ces signaux. Les signaux x(n), xb(n) et y(n) passés en cadence trames produisent respectivement les signaux x[m], xb[m], et y[m].
Dans une seconde étape a2, une détection d'activité vocale (DAV) est effectuée sur le signal x[m] de manière à déterminer si chaque trame respective courante d'indice m des signaux xb[m] et y[m], est une trame contenant seulement du bruit, notée "m bruit", ou une trame contenant de la parole, c'est-à-dire du signal utile, et notée "m_parole". Cette détermination se fait par comparaison des signaux xb[m] et y[m] avec le signal de test x[m] dénué de bruit. Chaque trame de silence de x[m] correspond en effet à une trame de bruit pour les signaux xb[m] et y[m], tandis que chaque trame de parole de x[m] correspond à une trame de parole pour les signaux xb[m] et y[m] Comme représenté sur la figure 2, en sortie de l'étape a2, trois types de trames sont sélectionnés à partir des signaux x[m], xb[m] et y[m] : - les trames de parole du signal bruité xb[m], notées xb[m_parole], - les trames de parole du signal traité y[m], notées y[m_parole], - les trames de bruit du signal traité y[m], notées y[m_bruit].
Dans une troisième étape a3, des mesures de sonie sont effectuées sur au moins des ensembles de trames y[m_bruit], y[m_parole], xb[m_parole] issues de l'étape précédente a2, et au moins un ensemble de trames du signal y[m] en sortie de l'étape al. Par exemple si on utilise 8 secondes de signal de test échantillonné à 8kHz, on pourra travailler sur 250 trames y[m] de 256 échantillons de signal y(n). De plus les coefficients de tonalité d'au moins un ensemble de trames y[m_bruit] sont mesurées.
Plus précisément, à cette étape, on calcule les densités de sonie moyennes Sxb(m _ parole), Sy(m _ parole), Sy(m), et Sy(m _ bruit) de respectivement chacune des trames xb[m_parole], y[m_parole], y[m] et y[m_bruit] des ensembles de trames considérés. De même les coefficients de tonalité ay(m_bruit) de chacune des trames y[m_bruit] de l'ensemble considéré de trames y[m_bruit] sont calculés.
Le calcul d'une densité de sonie moyenne Su (m) et d'un coefficient de tonalité a(m) d'une trame d'indice m quelconque d'un signal audio donné u, sera détaillé plus loin en liaison avec la figure 4.
Dans une quatrième étape a4, on calcule les moyennes respectives Sxb_parole, SY_parole, Sy, et SY_bruit des densités de sonie moyenne Sxb(m _ parole) , S- y(m _ parole) , S- y(m) , et Sy(m _bruit) précédemment calculées sur les ensembles respectifs considérés des trames xb[m_parole], y[m_parole], y[m] et y[m_bruit]. La moyenne ay_bru!tdes coefficients de tonalité ay(m _bruit) précédemment calculés sur l'ensemble considéré de trames y[m_bruit] est également calculée.
Dans une cinquième étape a5, on calcule cinq facteurs facteur(i), i étant un entier variant de un à cinq, caractéristiques de la gêne due au bruit dans le signal y(n), selon les formules suivantes: facteur(1)= 'S'Y brui'
SY
facteur(2)= SY _-bru S'Y _ parole facteur(3)= Ecart_type (Sxb(m _ parole) SY(m _ parole) , l'opérateur "Ecart_type (v(m))" désignant l'écart-type de la variable v sur l'ensemble des trames m, facteur(4)= aY bruit, facteu r(5)= Ecart_type (aY(m _ bruit)) . Dans une sixième étape a6, le calcul d'une note objective intermédiaire NOB est obtenue par combinaison linéaire des cinq facteurs calculés à l'étape a5, suivant l'équation suivante: NOB = E cofacteur(i) + C06, i=l où les coefficients col à 0)6 sont des coefficients de pondération prédéfinis. Ces coefficients ont été déterminés de manière à obtenir une corrélation maximale entre les données subjectives issues d'une base de données de tests subjectifs, et les notes objectives NOB calculées par cette combinaison linéaire en utilisant les signaux de tests, bruités et traités x[m], xb[m] et y[m] utilisés lors de ces mêmes tests subjectifs. La base de données de tests subjectifs est par exemple une base de données de notes obtenues avec des groupes d'auditeurs conformément à la "Recommandation UIT-T P.835", dans laquelle ces notes sont appelées notes "background noise".
II est à noter que l'obtention des coefficients de pondération par l'utilisation d'une base de données de tests subjectifs n'est pas indispensable à chaque étape de calcul d'une note objective NOB. En effet, ces coefficients doivent être obtenus préalablement à la première utilisation du procédé, et peuvent être les mêmes pour toutes les utilisations du procédé. Ces coefficients sont néanmoins amenés à évoluer lorsque de nouvelles données subjectives viendront alimenter la base de données de tests subjectifs utilisée.
Enfin dans une dernière étape a7, une note objective NOB_MOS de la gêne due au bruit dans le signal traité y(n) sur l'échelle MOS est calculée en utilisant par exemple une fonction polynomiale d'ordre 3, suivant l'équation suivante: NOB _ MOS = a,, (NOB)-' , où les coefficients X l à 24 sont déterminés de manière à ce que la note objective obtenue NOB_MOS caractérise la gêne due au bruit sur l'échelle MOS, c'est-à-dire sur une échelle de 1 à 5.
Selon un second mode de réalisation du procédé l'invention, la gêne due à la présence de bruit dans un signal audio bruité quelconque est évaluée de manière objective. On utilise le même environnement: de test qu'à la figure 1, mais en ôtant le module MRB de réduction de bruit. La source de signaux audio SSA délivre un signal audio de test x(n) ne contenant que du signal utile, auquel est ajouté un signal de bruit prédéfini généré par la source de bruit SB, pour obtenir en sortie de l'opérateur d'addition AD un signal bruité xb(n).
Le signal de test x(n) et le signal bruité xb(n) sont alors directement envoyés à l'entrée de l'équipement de test EQT mettant en oeuvre un procédé d'évaluation objective de la gêne due au bruit dans le signal bruité xb(n) selon l'invention. Comme dans le premier mode de réalisation, les signaux x(n) et xb(n) sont supposés échantillonnés à la fréquence d'échantillonnage 8 kHz.
L'équipement de test EQT délivre en sortie un résultat d'évaluation RES,qui est une note objective NOB_MOS de la gêne due à la présence de bruit dans le signal bruité xb(n).
En référence à la figure 3, le procédé de calcul de la note objective NOB_MOS de la gêne due au bruit dans le signal bruité xb(n) selon l'invention est représenté sous la forme d'un algorithme comportant des étapes b1 à b7. Ces étapes sont similaires aux étapes al à a7 précédemment décrites dans le premier mode de réalisation, et seront donc un peu moins détaillées. II est en effet à noter que si l'on applique les étapes de calcul a3 à a7 avec le signal y(n) égal au signal xb(n) dans le cas du premier mode de réalisation, on aboutit au deuxième mode de réalisation.
Dans une première étape b1, les signaux x(n) et xb(n) sont découpés en trames x[m] et xb[m] d'indice temporel m.
Dans une seconde étape b2, une détection d'activité vocale est effectuée sur le signal x[m] de manière à déterminer si chaque trame courante d'indice m du signal bruité xb[m] est une trame contenant seulement du bruit, notée "m bruit", ou une trame contenant aussi de la parole, notée "m_parole". Deux types de trames sont donc sélectionnés à partir des signaux x[m] et xb[m] en sortie de l'étape b2: - les trames de parole du signal bruité xb[m], notées xb[m_parole], - et les trames de bruit du signal bruité xb[m], notées xb[m_bruit].
Dans une troisième étape b3, des mesures de sonie sont effectuées sur au moins des ensembles de trames xb[m_bruit] et xb[m_parole] issues de l'étape précédente b2, et au moins un ensemble de trames du signal xb[m] en sortie de l'étape b1. De plus les coefficients de tonalité d'au moins un ensemble de trames xb[m_bruit] sont mesurées.
Plus précisément, à cette étape, on calcule les densités de sonie moyennes Sxb(m), Sxb(m_ parole) et Sbb(m_bruit) de respectivement chacune des trames xb[m], xb[m_parole] et xb[m_bruit] des ensembles de trames considérés. De même les coefficients de tonalité axb(m_bruit) de chacune des trames xb[m_bruit] de l'ensemble considéré de trames xb[m_bruit] sont calculés.
Dans une quatrième étape b4, on calcule les moyennes respectives SXb SXb _ parole, et SXb bruit des densités de sonie moyenne SXb(m), SXb(m _ parole) et Sxb(m _bruit) précédemment calculées sur les ensembles respectifs considérés des trames xb[m], xb[m_parole] et xb[m_bruit]. La moyenne axb_bruit des coefficients de tonalité axb(m__bruit) précédemment calculés sur l'ensemble considéré de trames xb[m_bruit] est également calculée.
Dans une cinquième étape b5, on calcule quatre facteurs facteur(i), i étant un entier variant de un à quatre, caractéristiques de la gêne due au bruit dans le signal bruité xb(n), selon les formules suivantes: facteur(1)= Sxb-bruit SXb facteur(2)= s'ai-bruit SXb_ parole facteur(3)= axb bruit facteur(4)= Ecart_type(avb(m_bruit)), l'opérateur "Ecart_type (v(m))" désignant l'écart-type de la variable v sur l'ensemble des trames m.
Dans une sixième étape b6, le calcul,d'une note objective intermédiaire NOB est obtenue par combinaison linéaire des quatre facteurs calculés à l'étape b5, suivant l'équation suivante: NOB =1 cofacteur (i) +W5, i=l où les coefficients col à (05 sont des coefficients de pondération prédéfinis. Ces coefficients ont été déterminés de manière à obtenir une corrélation maximale entre les données subjectives issues d'une base de données de tests subjectifs, et les notes objectives NOB calculées par cette combinaison linéaire en utilisant les signaux de tests et les signaux bruités x[m] et xb[m] utilisés lors de ces mêmes tests subjectifs. Tout comme pour l'étape a6, l'obtention des coefficients de pondération par l'utilisation d'une base de données de tests subjectifs n'est pas indispensable à chaque étape de calcul d'une note objective NOB.
Enfin dans une dernière étape b7, une note objective NOB_MOS de la gêne due au bruit dans le signal bruité xb(n) sur l'échelle MOS est calculée en utilisant par exemple une fonction polynomiale d'ordre 3, suivant l'équation suivante: NOB MOS = 12i (NOB)'-' , i=1 où les coefficients XI à X4 sont déterminés de manière à ce que la note objective obtenue NOB_MOS caractérise la gêne due au bruit sur l'échelle MOS, c'est-à-dire sur une échelle de 1 à 5.
Le calcul de densité de sonie moyenne et du coefficient de tonalité d'une trame d'un signal audio, utilisé dans les étapes a3 et b3, est maintenant décrit en relation avec la figure 4, selon un mode de réalisation préféré de l'invention.
Le calcul selon l'invention de la densité de sonie moyenne Su (m) d'une trame d'indice m quelconque d'un signal audio donné u[m], comprend les étapes cl à c7 représentées à la figure 4 et exposées ci-après. Le calcul selon l'invention du coefficient de tonalité a(m) d'une trame d'indice m quelconque d'un signal audio donné u[m], comprend les étapes cl, c2, c3 et c8 représentées à la figure 4 et exposées ci-après.
Dans ce qui suit, on considère une trame d'indice m quelconque d'un signal u[m], sachant que tout ou partie des trames du signal considéré subissent le même traitement. Le signal u[m] représente n'importe lequel des signaux x[m], xb[m], ou y[m] définis plus haut.
A la première étape cl, on applique à la trame d'indice m du signal u[m] un fenêtrage, par exemple un fenêtrage de type Hanning, Hamming ou équivalent. On obtient alors une trame fenêtrée u_w[m].
A l'étape suivante c2, on applique à la trame fenêtrée u_w[m], une transformée de Fourier rapide (FFT) et on obtient en conséquence une trame correspondante U(m,f) dans le domaine fréquentiel.
A l'étape suivante c3, on calcule la densité spectrale de puissance yu(m, f) de la trame U(m,f). Un tel calcul est connu de l'homme du métier et ne sera pas, par conséquent, détaillé ici.
A l'issue de l'étape c3, pour le signal y[m_bruit] de l'étape a3 ou le signal xb[m_bruit] de l'étape b3, on passe par exemple à l'étape c8 pour le calcul du coefficient de tonalité, puis à l'étape c4 pour le calcul de la densité de sonie moyenne Su (m) , puisque pour ces deux signaux les deux calculs sont nécessaires. Pour les autres signaux des étapes a3 et b3 on passe à l'étape c4 pour le calcul de la densité de sonie moyenne Su (m) . Il est à noter que le calcul du coefficient de tonalité est indépendant du calcul de la densité de sonie moyenne Su (m) , les deux calculs peuvent donc s'effectuer en parallèle ou l'un après l'autre.
A l'étape c4, on applique à la densité spectrale de puissance yu(m, f) obtenue à l'étape précédente, une conversion de l'axe des fréquences à l'échelle des Barks, et on obtient en conséquence une densité spectrale de puissance, Bu(m,b), sur l'échelle des Barks, appelée aussi spectre de Bark.
Pour une fréquence d'échantillonnage de 8kHz, 18 bandes critiques doivent être considérées. Ce type de conversion est connu de l'homme du métier, le principe de cette conversion Hertz/Bark consiste à additionner toutes les contributions fréquentielles présentes dans la bande critique considérée de l'échelle des Barks.
Ensuite, à l'étape c5, on applique à la densité spectrale de puissance sur l'échelle des Barks, Bu(m,b), une convolution avec la fonction d'étalement couramment utilisée en psychoacoustique, et on obtient en conséquence une densité spectrale étalée sur l'échelle des Barks, notée Eu(m,b). Cette fonction d'étalement a été formulée mathématiquement et une expression possible est: 101oglO(E(b)) =15.81+7. 5*(b+0.474) 17.5*. \/(1+(b+0.474)2) , où E(b) est la fonction d'étalement appliquée à la bande critique b considérée dans l'échelle des Barks et * symbolise l'opérateur de multiplication dans l'espace des nombres réels. Cette étape permet de prendre en compte l'interaction des bandes critiques adjacentes.
A l'étape suivante c6, on convertit la densité spectrale étalée E u (m, b) obtenue précédemment en densités de sonie exprimées en sones.
Pour cela, on opère une calibration de la densité spectrale étalée sur l'échelle des Barks, Eu(m,b), par les facteurs respectifs d'échelonnement en puissance et d'échelonnement en sonie couramment utilisés en psychoacoustique. Le document "Recommandation UIT-T P.862", sections 10.2. 1.3 et 10.2.1.4, donne un exemple d'une telle calibration par les facteurs précités. On convertit ensuite sur l'échelle des phones la grandeur obtenue. La conversion sur l'échelle des phones est effectuée en s'appuyant sur les courbes d'isosonie (courbes de Fletcher) conformément à la norme NF ISO 226 "Lignes isosoniques normales". On effectue alors une conversion sur l'échelle des sones de la grandeur précédemment convertie en phones. La conversion en sones est effectuée conformément à la loi de Zwicker selon laquelle: (N(phone) 40" N(sone) = 2 10 Pour obtenir plus d'information sur la conversion phone/sone, on pourra se reporter au document "PSYCHOACOUSTIQUE, L'oreille récepteur d'information", de E. Zwicker et R. Feldtkeller, édition Masson, 1981.
A l'issue de l'étape c6, on dispose d'un nombre B de valeurs de densité de sonie, Su(m,b), de la trame d'indice m pour la bande critique b, B étant le nombre de bandes critiques considérées dans l'échelle des Barks et l'indice b variant de 1 à B. Enfin, à l'étape c7, on calcule la densité de sonie rnoyenne Su (m) de la trame d'indice m à partir desdites B valeurs de densité de sonie, selon l'équation suivante:
B
Su(m)= ISu(m,b) B b=1 Autrement dit, la densité de sonie moyenne Su(m) selon l'invention d'une trame d'indice m, est donc la moyenne des B valeurs de densité de sonie Su(m,b), de la trame d'indice m pour une bande critique b considérée.
Ces deux dernières étapes c6 et c7 correspondent à une conversion du domaine des Barks vers le domaine des Sones, permettant de calculer une intensité subjective moyenne, c'est-à-dire telle que perçue par l'oreille humaine.
En outre à l'étape c8, le coefficient de tonalité a(rr1) de la trame d'indice m est calculé selon l'équation suivante: /( 1/N N-1 lYu (m, f) f=o N-1 EYu(m,f) N f.o - 60 où * symbolise l'opérateur de multiplication dans l'espace des nombres réels, f représente l'indice fréquentiel de la densité spectrale de puissance, et N désigne la taille de la transformée de Fourier rapide. Ce calcul est effectué selon le principe défini par J. D. Johnston dans son article "Transform coding of audio signais using perceptual noise criteria" du journal "IEEE Journal on selected areas in communications, vol.6, n 2, February 1988".
Le coefficient de tonalité a d'un signal de base est une mesure permettant de montrer si certaines fréquences pures ressortent de ce signal. Il est équivalent à une densité tonale. En effet, plus le coefficient de tonalité a est proche de 0, plus le signal est assimilé à du bruit. A l'inverse, plus le coefficient de tonalité a est proche de 1, plus le signal est à composante * 1og10 a(m) = tonale majoritaire. Un coefficient de tonalité a proche de 1 atteste donc de la présence de signal utile, ou signal de parole.
Claims (13)
1. Procédé de calcul d'une note objective (NOB) de la gêne due au bruit dans un signal audio traité par une fonction de réduction de bruit, ledit procédé comprenant une étape préalable d'obtention d'un signal audio prédéfini de test (x[m]) contenant un signal utile dépourvu de bruit, d'un signal bruité (xb[m]), obtenu en additionnant un signal de bruit prédéfini audit signal de test (x[m]), et d'un signal traité (y[m]), obtenu par application de la fonction de réduction de bruit audit signal bruité (xb[m]), ledit procédé étant caractérisé en ce qu'il inclut une étape (a3, a4) de mesures de sonie de trames dudit signal bruité (xb[m]) et dudit signal traité (y[m]), et de mesures de coefficients de tonalité de trames dudit signal traité (y[m]).
2. Procédé selon la revendication 1, caractérisé en ce qu'il comprend les étapes de: - Calcul (a3) de densités de sonie moyenne -s-y(m) de trames du signal traité (y[m]), de densités de sonie moyenne respectives Sxb(m _ parole) et Sr(m _ parole) de trames de signal utile "m_parole" respectivement du signal bruité (xb[m]) et du signal traité (y[m]), de densités de sonie moyenne Sr(m bruit) de trames de bruit "m_bruit" du signal traité (y[m]), et de coefficients de tonalité ar(m _bruit) de trames de bruit "m_bruit" du signal traité (y[m]), - Calcul (a5, a6) d'une note objective (NOB) de la gêne due au bruit dans le signal traité (y[m]), à partir desdites densités de sonie moyenne et desdits coefficients de tonalité calculés, et de coefficients de pondération prédéfinis.
3. Procédé selon la revendication 2, caractérisé en ce que l'étape de calcul (a3) de densités de sonie moyenne et de coefficients de tonalité est suivie d'une étape de calcul (a4) des moyennes Sy, Sxb parole, SY_ parole,SY_bruit et aY bruit desdites densités de sonie moyenne et desdits coefficients de tonalité sur l'ensemble des trames concernées des signaux correspondants, et en ce que la note objective (NOB) de la gêne due au bruit est calculée selon l'équation suivante: NOB = cofacteur(i) + C06 i=1 où facteur(1)= SY 2rutt Sr facteur(2)= SY -bruit SY _ parole facteur(3)= Ecart_type (Sxb(m _ parole) Sr (m _ parole)) , l'opérateur "Ecart_type (v(m))" désignant l'écart-type de la variable v sur l'ensemble des trames d'indice m, facteur(4)= aY bruit, facteur(5)= Ecart_type (al, (m _ bruit)) , et les coefficients w, à w6 sont déterminés de manière à obtenir une corrélation maximale entre les données subjectives issues d'une base de données de tests subjectifs et les notes objectives (NOB) calculées par ledit procédé pour les signaux de tests, bruités et traités (x[m], xb[m], y[m]) correspondants utilisés lors desdits tests subjectifs.
4. Procédé de calcul d'une note objective (NOB) de la gêne due au bruit dans un signal audio, ledit procédé comprenant une étape préalable d'obtention d'un signal audio prédéfini de test (x[m]) contenant un signal utile dépourvu de bruit, et d'un signal bruité (xb[m]), obtenu en additionnant un signal de bruit prédéfini audit signal de test (x[m]), ledit procédé étant caractérisé en ce qu'il inclut une étape (b3, b4) de mesures de sonie et de mesures de coefficients de tonalité de trames dudit signal bruité (xb[m]).
5. Procédé selon la revendication 4, caractérisé en ce qu'il comprend les étapes de: - Calcul (b3) de densités de sonie moyenne -Sv(m) de trames du signal bruité (xb[m]), de densités de sonie moyenne Sxb(m _ parole) de trames de signal utile "m_parole" du signal bruité (xb[m]), de densités de sonie moyenne Sx- b(m_bruit)de trames de bruit "m_bruit" du signal bruité (xb[m]), et de coefficients de tonalité axb(m_bruit) de trames de bruit "m_bruit" du signal bruité (xb[m]), - Calcul (b5, b6) d'une note objective (NOB) de la gêne due au bruit dans le signal bruité (xb[m]), à partir desdites densités de sonie moyennes et desdits coefficients de tonalité calculés, et de coefficients de pondération prédéfinis.
6. Procédé selon la revendication 5, caractérisé en ce que l'étape de calcul (b3) de densités de sonie moyenne et de coefficients de tonalité est suivie d'une étape de calcul (b4) des moyennes S x b, Sxb_ parole, Sxb _bruit et axb bruit desdites densités de sonie moyenne et desdits coefficients de tonalité sur l'ensemble des trames concernées des signaux correspondants, et en ce que ladite note objective (NOB) de la gêne due au bruit est calculée selon l'équation suivante: NOB = cofacteur(i) + cos où facteur(1)= Sxb-bruit Sxb facteur(2)= Sxb_brui' Sxb _ parole facteur(3) = a bruit, facteur(4)= Ecart type (axb(m _ bruit)), l'opérateur "Ecart type (v(m))" désignant l'écart-type de la variable v sur l'ensemble des trames d'indice m, et les coefficients û1 à (05 sont déterminés de manière à obtenir une corrélation maximale entre les données subjectives issues d'une base de données de tests subjectifs et les notes objectives (NOB) calculées par ledit procédé pour les signaux de tests et les signaux bruités (x[m], xb[m]) correspondants utilisés lors desdits tests subjectifs.
7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que ladite étape de calcul (a3, b3, a4, b4) de densités de sonie et de coefficients de tonalité est précédée d'une étape (a2, b2) de détection d'activité vocale sur le signal de test, de manière à déterminer si une trame courante d'indice m du signal bruité (xb[m]), et du signal traité (y[m]) dans le cas des revendications 1 à 3, est une trame "m bruit" contenant seulement du bruit, ou une trame "m_parole" contenant de la parole, dite trame de signal utile.
8. Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce que l'étape de calcul (a6, b6) de la note objective (NOB) est suivie d'une étape de calcul (a7, b7) d'une note objective sur l'échelle MOS (NOB_MOS) de la gêne due au bruit, calculée selon l'équation suivante: NOB _ MOS = (NOB)-' , i=l dans laquelle les coefficients 2,1 à X4 sont déterminés de manière à ce que ladite nouvelle note objective obtenue (NOB_MOS) caractérise la gêne due au bruit sur l'échelle MOS.
9. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que, dans l'étape de calcul (a3, b3, a4, b4) de densités de sonie et de coefficients de tonalité, le calcul de la densité de sonie moyenne Su(m) d'une trame d'indice m quelconque d'un signal audio donné u, comprend les étapes suivantes: - fenêtrage (cl), par exemple de type Hanning, de la trame d'indice m et obtention d'une trame fenêtrée u_w[m], - application (c2) d'une transformée de Fourier rapide à la trame fenêtrée u_w[m] et obtention d'une trame correspondante U(m,f) dans le domaine fréquentiel, calcul (c3) de la densité spectrale de puissance yu(m,f) de la trame U(m,f), - application (c4) à la densité spectrale de puissance yu(m,f) d'une conversion de l'axe des fréquences à l'échelle des Barks et obtention d'une densité spectrale de puissance Bu(m,b) sur l'échelle des Barks, - convolution (c5) de la densité spectrale de puissance sur l'échelle des Barks, Bu(m,b), avec la fonction d'étalement couramment utilisée en psychoacoustique et obtention d'une densité spectrale étalée sur l'échelle des Barks, Eu(m,b), - calibration (c6) de la densité spectrale étalée sur l'échelle des Barks, Eu(m,b), par les facteurs respectifs d'échelonnement en puissance et d'échelonnement en sonie couramment utilisés en psychoacoustique, conversion de la grandeur ainsi obtenue sur l'échelle des phones puis conversion sur l'échelle des sones de la grandeur précédemment convertie en phones, et obtention en conséquence d'un nombre B de valeurs de densité de sonie, Su(m,b), de la trame d'indice m pour la bande critique b, B étant le nombre de bandes critiques considérées dans l'échelle des Barks et l'indice b variant de 1 à B, - calcul (c7) de la densité de sonie moyenne Su(m) de la trame d'indice m à partir desdites B valeurs de densités de sonie Su(m,b), selon l'équation suivante:
B
Su(m)= 1 ESu(m,b) B b=1
10. Procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce que, dans l'étape de calcul (a3, b3, a4, b4) de densités de sonie et de coefficients de tonalité, le calcul du coefficient de tonalité a(m) d'une trame d'indice m quelconque d'un signal audio donné u, comprend les étapes suivantes: - fenêtrage (cl) , par exemple de type Hanning, de la trame d'indice m et obtention d'une trame fenêtrée u_w[m], - application (c2) d'une transformée de Fourier rapide à la trame fenêtrée u_w[m] et obtention d'une trame correspondante U(m,f) dans le 15 domaine fréquentiel, - calcul (c3) de la densité spectrale de puissance yu (m, f) de la trame U(m,f), - calcul (c8) du coefficient de tonalité a(m) selon l'équation suivante: N-1 fI Yu (m, f) \I=0 i N-I -EYu(m,f) N f=0 \1IN \ 10*1og10 a(m) = 60 où * symbolise l'opérateur de multiplication dans l'espace des nombres réels, f représente l'indice fréquentiel de la densité spectrale de puissance, et N désigne la taille de la transformée de Fourier rapide.
11. Equipement de test destiné à évaluer une note objective de la gêne 25 due au bruit dans un signal audio, caractérisé en ce qu'il comporte des moyens adaptés à mettre en oeuvre un procédé selon l'une quelconque des revendications 1 à 10.
12. Equipement de test selon la revendication 11, caractérisé en ce qu'il inclut des moyens informatiques et un programme d'ordinateur, ledit programme comportant des instructions adaptées à mettre en oeuvre ledit procédé, lorsqu'il est exécuté par lesdits moyens informatiques.
13. Programme d'ordinateur sur un support d'informations, caractérisé en ce qu'il comporte des instructions adaptées à la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à 10, lorsque le programme est chargé et exécuté dans un système informatique.
Priority Applications (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0501747A FR2882458A1 (fr) | 2005-02-18 | 2005-02-18 | Procede de mesure de la gene due au bruit dans un signal audio |
| ES06709505T ES2329932T3 (es) | 2005-02-18 | 2006-02-13 | Procedimiento de medida de la molestia debida al ruido de una señal de audio. |
| PCT/FR2006/050126 WO2006087490A1 (fr) | 2005-02-18 | 2006-02-13 | Procede de mesure de la gene due au bruit dans un signal audio |
| AT06709505T ATE438173T1 (de) | 2005-02-18 | 2006-02-13 | Verfahren zur messung von durch geräusche in einem audiosignal verursachten beeinträchtigungen |
| US11/884,573 US20080267425A1 (en) | 2005-02-18 | 2006-02-13 | Method of Measuring Annoyance Caused by Noise in an Audio Signal |
| EP06709505A EP1849157B1 (fr) | 2005-02-18 | 2006-02-13 | Procede de mesure de la gene due au bruit dans un signal audio |
| DE602006008111T DE602006008111D1 (de) | 2005-02-18 | 2006-02-13 | Verfahren zur messung von durch geräusche in einem audiosignal verursachten beeinträchtigungen |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0501747A FR2882458A1 (fr) | 2005-02-18 | 2005-02-18 | Procede de mesure de la gene due au bruit dans un signal audio |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| FR2882458A1 true FR2882458A1 (fr) | 2006-08-25 |
Family
ID=34981381
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0501747A Pending FR2882458A1 (fr) | 2005-02-18 | 2005-02-18 | Procede de mesure de la gene due au bruit dans un signal audio |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US20080267425A1 (fr) |
| EP (1) | EP1849157B1 (fr) |
| AT (1) | ATE438173T1 (fr) |
| DE (1) | DE602006008111D1 (fr) |
| ES (1) | ES2329932T3 (fr) |
| FR (1) | FR2882458A1 (fr) |
| WO (1) | WO2006087490A1 (fr) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113473314A (zh) * | 2020-03-31 | 2021-10-01 | 华为技术有限公司 | 音频信号处理方法以及相关设备 |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102005052630A1 (de) * | 2005-08-25 | 2007-03-01 | Robert Bosch Gmbh | Verfahren und Vorrichtung zur Bewertung der Lästigkeit von Quietschgeräuschen |
| KR100810077B1 (ko) * | 2006-05-26 | 2008-03-05 | 권대훈 | 표준 청력특성을 이용한 이퀄라이제이션 방법 |
| JP5270006B2 (ja) | 2008-12-24 | 2013-08-21 | ドルビー ラボラトリーズ ライセンシング コーポレイション | 周波数領域におけるオーディオ信号ラウドネス決定と修正 |
| EP2685448B1 (fr) * | 2012-07-12 | 2018-09-05 | Harman Becker Automotive Systems GmbH | Synthèse de son moteur |
| KR101757341B1 (ko) | 2013-01-29 | 2017-07-14 | 프라운호퍼 게젤샤프트 쭈르 푀르데룽 데어 안겐반텐 포르슝 에.베. | 저-복잡도 음조-적응 오디오 신호 양자화 |
| EP3471605B1 (fr) | 2016-06-17 | 2023-08-09 | Predictive Safety Srp, Inc. | Système et procédé de détection d'invalidité |
| CN110688712A (zh) * | 2019-10-11 | 2020-01-14 | 湖南文理学院 | 汽车风振噪声声品质客观烦恼度评价指标及其计算方法 |
| CN116429245B (zh) * | 2023-06-13 | 2023-09-01 | 江铃汽车股份有限公司 | 一种雨刮电机噪声测试方法及系统 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6651041B1 (en) * | 1998-06-26 | 2003-11-18 | Ascom Ag | Method for executing automatic evaluation of transmission quality of audio signals using source/received-signal spectral covariance |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5574824A (en) * | 1994-04-11 | 1996-11-12 | The United States Of America As Represented By The Secretary Of The Air Force | Analysis/synthesis-based microphone array speech enhancer with variable signal distortion |
| FI100840B (fi) * | 1995-12-12 | 1998-02-27 | Nokia Mobile Phones Ltd | Kohinanvaimennin ja menetelmä taustakohinan vaimentamiseksi kohinaises ta puheesta sekä matkaviestin |
| US6446038B1 (en) * | 1996-04-01 | 2002-09-03 | Qwest Communications International, Inc. | Method and system for objectively evaluating speech |
| FI114833B (fi) * | 1999-01-08 | 2004-12-31 | Nokia Corp | Menetelmä, puhekooderi ja matkaviestin puheenkoodauskehysten muodostamiseksi |
| US6490552B1 (en) * | 1999-10-06 | 2002-12-03 | National Semiconductor Corporation | Methods and apparatus for silence quality measurement |
| FI116643B (fi) * | 1999-11-15 | 2006-01-13 | Nokia Corp | Kohinan vaimennus |
| DE60104091T2 (de) * | 2001-04-27 | 2005-08-25 | CSEM Centre Suisse d`Electronique et de Microtechnique S.A. - Recherche et Développement | Verfahren und Vorrichtung zur Sprachverbesserung in verrauschte Umgebung |
| EP1760696B1 (fr) * | 2005-09-03 | 2016-02-03 | GN ReSound A/S | Méthode et dispositif pour l'estimation améliorée du bruit non-stationnaire pour l'amélioration de la parole |
-
2005
- 2005-02-18 FR FR0501747A patent/FR2882458A1/fr active Pending
-
2006
- 2006-02-13 ES ES06709505T patent/ES2329932T3/es active Active
- 2006-02-13 AT AT06709505T patent/ATE438173T1/de not_active IP Right Cessation
- 2006-02-13 WO PCT/FR2006/050126 patent/WO2006087490A1/fr active Application Filing
- 2006-02-13 US US11/884,573 patent/US20080267425A1/en not_active Abandoned
- 2006-02-13 EP EP06709505A patent/EP1849157B1/fr not_active Not-in-force
- 2006-02-13 DE DE602006008111T patent/DE602006008111D1/de not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6651041B1 (en) * | 1998-06-26 | 2003-11-18 | Ascom Ag | Method for executing automatic evaluation of transmission quality of audio signals using source/received-signal spectral covariance |
Non-Patent Citations (3)
| Title |
|---|
| BEERENDS J G ET AL: "Perceptual Evaluation of Speech Quality (PESQ), the new ITU standard for end-to-end speech quality assessment. Part II - Psychoacoustic model", AES: JOURNAL OF THE AUDIO ENGINEERING SOCIETY, June 2001 (2001-06-01), XP002206026 * |
| JOHNSTON J D: "TRANSFORM CODING OF AUDIO SIGNALS USING PERCEPTUAL NOISE CRITERIA", IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS, IEEE INC. NEW YORK, US, vol. 6, no. 2, 1 February 1988 (1988-02-01), pages 314 - 323, XP002003779, ISSN: 0733-8716 * |
| RIX A W ET AL: "The perceptual analysis measurement system for robust end-to-end speech quality assessment", ACOUSTICS, SPEECH, AND SIGNAL PROCESSING, 2000. ICASSP '00. PROCEEDINGS. 2000 IEEE INTERNATIONAL CONFERENCE ON 5-9 JUNE 2000, PISCATAWAY, NJ, USA,IEEE, vol. 3, 5 June 2000 (2000-06-05), pages 1515 - 1518, XP010507639, ISBN: 0-7803-6293-4 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113473314A (zh) * | 2020-03-31 | 2021-10-01 | 华为技术有限公司 | 音频信号处理方法以及相关设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| ES2329932T3 (es) | 2009-12-02 |
| WO2006087490A1 (fr) | 2006-08-24 |
| ATE438173T1 (de) | 2009-08-15 |
| EP1849157A1 (fr) | 2007-10-31 |
| US20080267425A1 (en) | 2008-10-30 |
| DE602006008111D1 (de) | 2009-09-10 |
| EP1849157B1 (fr) | 2009-07-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1849157B1 (fr) | Procede de mesure de la gene due au bruit dans un signal audio | |
| EP2419900B1 (fr) | Procede et dispositif d'evaluation objective de la qualite vocale d'un signal de parole prenant en compte la classification du bruit de fond contenu dans le signal | |
| Li et al. | On the importance of power compression and phase estimation in monaural speech dereverberation | |
| EP2415047B1 (fr) | Classification du bruit de fond contenu dans un signal sonore | |
| EP0002998B1 (fr) | Procédé de compression de données relatives au signal vocal et dispositif mettant en oeuvre ledit procédé | |
| EP1593116B1 (fr) | Procédé pour le traitement numérique différencié de la voix et de la musique, le filtrage de bruit, la création d'effets spéciaux et dispositif pour la mise en oeuvre dudit procédé | |
| JP5542206B2 (ja) | オーディオ・システムの知覚品質を判定する方法およびシステム | |
| CN104919525B (zh) | 用于评估退化语音信号的可理解性的方法和装置 | |
| Santos et al. | Speech dereverberation with context-aware recurrent neural networks | |
| Kumar | Comparative performance evaluation of MMSE-based speech enhancement techniques through simulation and real-time implementation | |
| Kumar | Real-time performance evaluation of modified cascaded median-based noise estimation for speech enhancement system | |
| RU2312405C2 (ru) | Способ осуществления машинной оценки качества звуковых сигналов | |
| FR2894707A1 (fr) | Procede de mesure de la qualite percue d'un signal audio degrade par la presence de bruit | |
| Westhausen et al. | Reduction of subjective listening effort for TV broadcast signals with recurrent neural networks | |
| US9659565B2 (en) | Method of and apparatus for evaluating intelligibility of a degraded speech signal, through providing a difference function representing a difference between signal frames and an output signal indicative of a derived quality parameter | |
| WO2020049263A1 (fr) | Dispositif de rehaussement de la parole par implementation d'un reseau de neurones dans le domaine temporel | |
| EP1229517B1 (fr) | Méthode de reconnaissance de la parole avec normalisation de la variance dépendante du bruit | |
| WO2006032751A1 (fr) | Procede et dispositif d'evaluation de l'efficacite d'une fonction de reduction de bruit destinee a etre appliquee a des signaux audio | |
| Pourmand et al. | Computational auditory models in predicting noise reduction performance for wideband telephony applications | |
| Ghimire | Speech intelligibility measurement on the basis of ITU-T Recommendation P. 863 | |
| JP6232710B2 (ja) | 録音音声の明瞭化装置 | |
| Rahali et al. | Enhancement of noise-suppressed speech by spectral processing implemented in a digital signal processor | |
| Dubey | Non-intrusive objective speech quality assessment using features at single and multiple time-scales | |
| FR3051958A1 (fr) | Procede et dispositif pour estimer un signal dereverbere | |
| CN115101082A (zh) | 语音增强方法、装置、设备、存储介质及程序产品 |