EP 2 000 002 B1 valeur de retard à chaque sous- bande de rang k, une valeur de retard correspondant au retard de propagation entre l’oreille gauche et l’oreille droite d’un auditeur pour différentes positions. [0086] Ainsi, à partir d’un filtre HRTF, il est possible de calculer de façon automatique les gains et les délais de retard appliqués en sous- bande. A partir de la résolution fréquentielle du banc de filtre HRTF, on associe à chacune des sous- bandes SB$_{k}$ une valeur de retard correspondant au retard de propagation entre l’oreille gauche et l’oreille droite d’un auditeur pour différentes positions. [0087] Ainsi, à partir d’un filtre HRTF, on peut calculer de façon automatique les gains et les délais de retard à appliquer en sous- bande. [0088] A partir de la résolution fréquentielle du banc de filtre, on associe à chacune des bandes une valeur réelle. A titre d’exemple non limitatif, il est possible à partir du module du filtre HRTF, de calculer, pour chaque sous- bande, la moyenne du module du filtre HRTF précité. Une telle opération est similaire à une analyse en bande d’octave ou de Bark des filtres HRTF. De même, on détermine le retard à appliquer pour les canaux indirects, c’est- à- dire les valeurs de retard qui sont applicables plus particulièrement aux canaux dont le retard n’est pas minimum. Il existe de nombreuses méthodes pour déterminer de manière automatique les retards interauraux encore désignés ITD pour « I nteraural Time Difference » et qui correspondent aux retards entre l’oreille gauche et l’oreille droite, pour différentes positions de l’auditeur. On peut utiliser, à titre d’exemple non limitatif, la méthode du seuil décrite par S. Busson dans la thèse de doctorat de l’Université de la Méditerranée Est- Marseille II, 2006, intitulée «Individualisation d’indices acoustiques pour la synthèse binaurale». Le principe des méthodes d’estimation du retard interaural de type seuil est de déterminer le temps d’arrivée, ou encore le retard initial de l’onde sur l’oreille droite Td et sur l’oreille gauche Tg. Le retard interaural est donné par la relation ITD seuil = Td - Tg. [0089] La méthode la plus courante estime le temps d’arrivée comme l’instant où le filtre temporel HRIR dépasse un seuil donné. Par exemple le temps d’arrivée peut correspondre au temps pour lequel la réponse du filtre HRIR atteint 10 % de son maximum. [0090] Un exemple de mise en oeuvre spécifique dans le domaine transformé PQMF sera maintenant donné ci- après. [0091] D’une manière générale, on indique que l’application d’un gain dans le domaine PQMF complexe consiste à multiplier la valeur de chaque échantillon du signal en sous- bande, représenté par une valeur complexe, par la valeur de gain formée par un nombre réel. [0092] En effet, il est bien connu que l’usage d’un domaine transformé PQMF complexe, permet d’appliquer les gains en s’affranchissant des problèmes de repliement de spectre engendrés par le sous- échantillonnage inhérent aux bancs de filtres. Chaque sous- bande SB$_{k}$ de chaque canal se voit ainsi affectée d’un gain déterminé. [0093] En outre, l’application d’un retard dans le domaine transformé PQMF consiste au moins, pour chaque échantillon du signal en sous- bande, représenté par une valeur complexe, à introduire une rotation dans le plan complexe par multiplication de cet échantillon par une valeur exponentielle complexe fonction du rang de la sous- bande considérée, du taux de sous- échantillonnage dans la sous- bande considérée et d’un paramètre de retard lié à la différence de retard interaural d’un auditeur. [0094] La rotation dans le plan complexe est alors suivie d’un retard temporel pur de l’échantillon après rotation. Ce retard temporel pur est une fonction de la différence du retard interaural d’un auditeur et du taux de sous échantillonnage dans la sous- GLYPH<3> bande considérée. [0095] De manière pratique, on indique que les retards précités sont appliqués sur les signaux résultants c’est- à- dire les signaux égalisés et en particulier sur les sous- GLYPH<3> ensembles de ces signaux ou canaux qui ne bénéficient pas d’une trajectoire directe. [0096] En particulier, la rotation est effectuée sous la forme d’une multiplication complexe par une valeur exponentielle de la forme :  et par un retard pur implémenté par une ligne à retard, par exemple réalisant l’opération :  [0097] Dans les relations précédentes : - exp est la fonction exponentielle ; - j est tel que j*j = -1 ; exp (-j*pi*(k+ 0,5)*d/M) y(k,n) = x(k, n-D) 8