FR 2844894 A1 20040326 FR 0211739 A 20020923 I La présente invention concerne un procédé et un système de traitement d'une représentation d'un champ acoustique afin de délivrer une représentation de ce champ acoustique modifiée par un traitement au moins spatial. Actuellement, il existe des procédés et des systèmes qui, à partir 5 d'une représentation d'un champ acoustique sous la forme de coefficients en fonction du temps et de l'espace, permettent de délivrer une représentation modifiée de ces coefficients par des combinaisons linéaires correspondant à un traitement au moins spatial du champ acoustique. Cependant, les traitements possibles sont limités et dans tous les cas 10 prédéterminés. Par exemple, les systèmes regroupés sous le terme " ambisonic " utilisent une représentation d'un environnement sonore sous la forme d'une distribution angulaire de sources virtuelles autour d'un point assimilé à la position d'écoute. Ces systèmes ne permettent d'effectuer que certaines modifications prédéterminées telles qu'une distorsion angulaire sur l'azimut, une focalisation dans une direction privilégiée, ou une rotation par rapport à un axe donné. De plus, ces systèmes sont limités au niveau de leur précision. Ainsi, dans les systèmes existants, des traitements tels que la focalisation dans une 20 direction quelconque ne peuvent être appliqués qu'à des représentations limitées à un niveau de précision correspondant au premier ordre mathématique, c'est à dire de faible qualité. Il apparaît donc qu'aucun des systèmes de l'art antérieur ne permet d'effectuer une modification quelconque d'une représentation d'un champ acous25 tique afin d'obtenir une représentation modifiée, notamment dans le cas d'une représentation initiale de qualité, d'ordre élevé par exemple. Le but de l'invention est de remédier à ce problème en fournissant un procédé et un système de traitement de la représentation d'un champ acoustique permettant d'effectuer des traitements quelconques, présentant une très large 30 diversité L'invention a pour objet un procédé de traitement d'une représentation d'un champ acoustique, afin de délivrer une représentation modifiée par application de ce traitement correspondant à un traitement au moins spatial du champ acoustique, comportant: l'établissement d'une représentation initiale du champ acoustique sous la forme d'un ensemble de coefficients représentatifs de ce champ dans le temps et dans les trois dimensions de l'espace; et - l'établissement d'une représentation modifiée par des combinaisons 5 de filtrage correspondant audit traitement, appliquées auxdits coefficients de ladite représentation initiale, caractérisé en ce qu'il comporte en outre une étape de définition dudit traitement, laquelle comprend: - une étape d'établissement d'un ensemble de fonctions de directivité 10 représentatif dudit traitement, sous la forme d'un ensemble de coefficients correspondant à la décomposition de ces fonctions sur une base d'harmoniques sphériques; et - une étape de détermination desdites combinaisons de filtrage à partir de cette décomposition pour la mise en oeuvre d'une étape d'applications des 15 combinaisons de filtrage ainsi définies auxdits coefficients de la représentation initiale, afin d'obtenir ladite représentation modifiée. Suivant d'autres caractéristiques de l'invention - ladite étape de détermination d'un ensemble de fonctions de directivité comporte une sousétape de détermination directe de tout ou partie dudit 20 ensemble de fonctions de directivité dans l'un ou l'autre d'une pluralité de formats prédéterminés; - ladite étape de détermination d'un ensemble de fonctions de directivité, comporte une sous-étape d'interprétation automatique d'une consigne de traitement permettant de délivrer automatiquement tout ou partie dudit ensemble de 25 fonctions de directivité dans l'un ou l'autre d'une pluralité de formats prédéterminés, à partir de ladite consigne de traitement; - ledit ensemble de fonctions de directivité représentatif du traitement est exprimé selon l'un ou l'autre d'une pluralité de formats prédéterminés et ladite étape d'établissement d'un ensemble de fonctions de directivité comporte une 30 sous-étape de conversion, à partir dudit format prédéterminé, en un ensemble de coefficients correspondant à la décomposition desdites fonctions de directivité sur une base d'harmoniques sphériques; - ladite étape de définition du traitement comprend en outre, une étape d'établissement d'un ensemble de paramètres dudit traitement; ledit ensemble de paramètres comprend tout ou partie des paramètres suivants: - un ensemble de constantes formant des coefficients de pondération et défini en fonction de la nature dudit traitement; et - un paramètre spécifiant l'ordre limite de fonctionnement du procédé et correspondant à la limite mathématique de précision souhaitée; - ladite étape d'établissement d'un ensemble de paramètres comporte une sous-étape de détermination directe de tout ou partie dudit ensemble de paramètres; ladite étape d'établissement d'un ensemble de paramètres comporte une sous-étape d'interprétation automatique d'une consigne de traitement permettant de délivrer automatiquement tout ou partie dudit ensemble de paramètres à partir de ladite consigne de traitement; - ladite étape d'établissement d'un ensemble de paramètres comprend 15 une sous-étape de calcul de tout ou partie dudit ensemble de constantes en fonction de la nature dudit traitement; - ladite étape de détermination desdites combinaisons de filtrage comprend au moins une sous-étape de calcul de combinaisons de filtrage par des combinaisons linéaires desdits coefficients correspondant à la décomposition sur 20 une base d'harmoniques sphériques dudit ensemble de fonctions de directivité représentatif dudit traitement; - ladite sous-étape de calcul par des combinaisons linéaires fait intervenir ledit ensemble de constantes en tant que coefficients de pondération; - il est répété une pluralité de fois et les combinaisons de filtrage résul25 tant des différentes sous-étapes de calcul sont combinées entre elles pour former les combinaisons de filtrage destinées à être appliquées auxdits coefficients de la représentation initiale; - le résultat de ladite sous-étape de calcul forme directement lesdites combinaisons de filtrage destinées à être appliquées auxdits coefficients de la 30 représentation initiale; ledit traitement est un traitement de convolution; - ledit traitement est un traitement de fenêtrage; - ledit traitement comprend au moins une distorsion à appliquer à la représentation initiale du champ acoustique et tout ou partie des fonctions de di- rectivité déterminées lors de l'étape d'établissement d'un ensemble de fonctions de directivité forment un ensemble de couples de distorsion représentatif de ladite distorsion sous la forme d'un ensemble de coefficients correspondant à la décomposition desdits couples de distorsion sur une base d'harmoniques sphéri5 ques; - ledit traitement comprend un premier traitement dit nominal et la ou lesdites distorsions, ladite sous-étape de calcul délivrant des combinaisons de filtrage correspondant audit traitement nominal et ladite étape de détermination de combinaisons de filtrage comprenant en outre au moins une sous-étape de 10 distorsion par modification des combinaisons de filtrage délivrées par ladite sousétape de calcul, à partir desdits couples de distorsion représentatifs de la ou desdites distorsions, afin de délivrer lesdites combinaisons de filtrage destinées à être appliquées auxdits coefficients de la représentation initiale; - ledit ensemble de paramètres comporte un paramètre représentatif de l'effort de 15 distorsion, soit encore spécifiant la ressemblance entre ledit traitement à appliquer et ledit traitement nominal; - ledit traitement comprend uniquement la ou lesdites distorsions, ladite étape de détermination des combinaisons de filtrage comprenant au moins une sous-étape de distorsion délivrant des combinaisons de filtrage uniquement 20 à partir desdits couples de distorsion représentatifs de la ou desdites distorsions, afin de délivrer lesdites combinaisons de filtrage destinées à être appliquées auxdits coefficients de la représentation initiale; - il comporte une pluralité de sous-étapes de distorsion, récursives les unes par rapport aux autres, de sorte que chaque sous-étape de distorsion reçoit 25 en entrée les combinaisons de filtrage délivrées par la sous-étape précédente; - ladite au moins une sous-étape de distorsion comprend une sousétape d'optimisation permettant de déterminer les combinaisons de filtrage correspondant à des modifications à effectuer sur des parties de la représentation initiale qui ne sont pas directement visées par lesdits couples de distorsion; 30 - ledit traitement se décompose en une pluralité de traitements élémentaires et au moins ladite étape de détermination de combinaisons de filtrage est répétée une pluralité de fois, les combinaisons de filtrage résultant des différentes sous-étapes de calcul et/ou de distorsion et correspondant auxdits traitements élémentaires étant arrangées entre elles pour former les combinaisons de filtrage destinées à être appliquées auxdits coefficients de la représentation initiale; - ladite étape de détermination des combinaisons de filtrage comprend une sous-étape de conversion desdites combinaisons de filtrage représentées 5 dans un format prédéterminé pour délivrer des paramètres de filtres adaptés pour le paramétrage de filtres donnés en vue de la réalisation de ladite étape d'application des combinaisons de filtrage; ledit traitement est un traitement spatial et ladite étape d'applications des combinaisons de filtrage comprend des modifications de gain appliquées 10 auxdits coefficients de la représentation initiale; - ledit traitement est un traitement spatial et fréquentiel et ladite étape d'application des combinaisons de filtrage comprend l'application auxdits coefficients de la représentation initiale, de filtres variant en fonction de la fréquence; lesdits ensembles de coefficients représentatifs de la représentation 15 initiale du champ acoustique et de la représentation modifiée sont des ensembles de coefficients dits de Fourier-Bessel. L'invention a également pour objet un dispositif de traitement d'une représentation d'un champ acoustique afin de délivrer une représentation de ce champ acoustique par application de ce traitement correspondant à un traitement 20 au moins spatial du champ acoustique, ledit dispositif recevant en entrée une représentation initiale du champ acoustique sous la forme d'un ensemble de coefficients représentatif de ce champ dans le temps et dans les trois dimensions de l'espace et délivrant en sortie une représentation modifiée par les combinaisons de filtrage correspondant audit traitement appliqué auxdits coefficients de la 25 représentation initiale, caractérisé en ce qu'il comporte des moyens de définition dudit traitement, comprenant: - des moyens d'établissement d'un ensemble de fonctions de directivité représentatif dudit traitement, sous la forme d'un ensemble de coefficients correspondant à la décomposition de ces fonctions de directivité sur une base 30 d'harmoniques sphériques; et - des moyens de détermination desdites combinaisons de filtrage à partir de cette décomposition pour leur mise en oeuvre par des moyens (80) de filtrage recevant en entrée lesdits coefficients de la représentation initiale et délivrant en sortie ladite représentation modifiée. Selon d'autres caractéristique de ce dispositif: - lesdits moyens de définition dudit traitement comportent un module de saisie directe d'informations de fonctionnement du dispositif dans l'un ou l'autre d'une pluralité de formats; - lesdits moyens de définition dudit traitement comportent un module d'interprétation automatique d'une consigne de traitement permettant de délivrer automatiquement des informations de fonctionnement dans l'un ou l'autre d'une pluralité de formats prédéterminés; - lesdites informations de fonctionnement comprennent ledit ensemble 10 de fonctions de directivité représentatif dudit traitement: lesdits moyens de définition dudit traitement comportent un module de conversion adapté pour convertir ledit ensemble de fonctions de directivité de l'un ou l'autre d'une pluralité de formats prédéterminés en un ensemble de coefficients correspondant à la décomposition desdites fonctions de directivité sur une 15 base d'harmoniques sphériques; lesdites informations de fonctionnement comprennent un ensemble de paramètres dudit traitement; - lesdits moyens de définition dudit traitement comportent un module de calcul de tout ou partie dudit ensemble de paramètres en fonction de la nature 20 dudit traitement: - lesdits moyens de détermination desdites combinaisons de filtrage comprennent un module de calcul de combinaisons de filtrage par combinaisons linéaires desdits coefficients correspondant à la décomposition sur une base d'harmoniques sphériques dudit ensemble de directivité représentatif dudit trai25 tement; - ledit ensemble des paramètres dudit traitement comprend un ensemble de constantes et ledit module de calcul par combinaisons linéaires fait intervenir ledit ensemble de constantes en tant que coefficients de pondération; - ledit traitement comprend au moins une distorsion, tout ou partie du30 dit ensemble de fonctions de directivité formant un ensemble de couples de distorsion correspondant à la ou aux distorsions à appliquer à la représentation initiale et lesdits moyens de détermination desdites combinaisons de filtrage comprennent en outre un module de distorsion par modification des combinaisons de filtrage délivrées par ledit module de calcul, à partir desdits couples de distorsion afin de délivrer lesdites combinaisons de filtrage destinées à être appliquées auxdits coefficients de la représentation initiale; ledit traitement comprend uniquement une ou plusieurs distorsions, tout ou partie dudit ensemble de fonctions de directivité formant un ensemble de 5 couples de distorsion correspondant à la ou aux distorsions à appliquer à la représentation initiale et lesdits moyens de détermination desdites combinaisons de filtrage comprennent uniquement un module de distorsion pour délivrer lesdites combinaisons de filtrage uniquement à partir desdits couples de distorsion. L'invention sera mieux comprise à la lecture de la description qui va 10 suivre donnée uniquement à titre d'exemple, et faite en se référant aux dessins annexés sur lesquels: - la Fig.1 est une représentation d'un repère sphérique; - la Fig.2 est un organigramme général du procédé de l'invention; - la Fig.3 est un organigramme détaillé d'un premier mode de réalisa15 tion de l'invention; - la Fig.4 est un schéma d'un mode de réalisation de l'étape d'application du traitement; - la Fig.5 est un organigramme détaillé d'un second mode de réalisation de l'invention; et la Fig.6 est un schéma synoptique d'un dispositif adapté pour la mise en oeuvre du procédé de l'invention. Sur la figure 1, on a représenté un repère sphérique classique, de manière à préciser le système de coordonnées auquel il est fait référence dans le texte. Ce repère est un repère orthonormal, d'origine 0, et comporte trois axes (OX), (OY) et (OZ). Dans ce repère, une position notée x est décrite au moyen de ses coordonnées sphériques (r,0,ç), o r désigne la distance par rapport à l'origine 0, O l'orientation dans le plan vertical et s l'orientation dans le plan horizontal. Dans un tel repère, un champ acoustique est connu si l'on définit en tout point à chaque instant t la pression acoustique notée p(r, O,,t), dont la transformée de Fourier est notée P(r, O, 0,J) o f désigne la fréquence. Le procédé de l'invention se base sur l'utilisation de fonctions spatiotemporelles permettant de décrire un champ acoustique quelconque dans le temps et dans les trois dimensions de l'espace. Dans les modes de réalisation décrits, ces fonctions sont des fonctions 5 dites de Fourier-Bessel sphériques de première espèce, appelées par la suite fonctions de Fourier-Bessel. Dans une zone vide de sources et vide d'obstacles, les fonctions de Fourier-Bessel correspondent aux solutions de l'équation des ondes et constituent une base qui engendre tous les champs acoustiques produits par des sour10 ces situées à l'extérieur de cette zone. Tout champ acoustique tridimensionnel peut donc s'exprimer par une combinaison linéaire des fonctions de Fourier-Bessel, selon l'expression suivante: P(r,9,QbJ)=4îrcE R,m(f)j'ji(kr)yp(0,ç$) P(r,0,q), f) =4 zc],. fij(lr)y[" O, l=0 M=-i Dans cette équation, les termes Pl,m(/) sont définis comme les coefficients de Fourier-Bessel du champ p(r,O,0,t), k= 24, c est la célérité du son C dans l'air (340 ms'1), jKkr) est la fonction de Bessel sphérique de première espèce d'ordre I définie par j, (x) = - J+,112(x) o Jv(x) est la fonction de Bessel de première espèce d'ordre v, et yp(O,S) est l'harmonique sphérique réelle d'ordre I 20 et de terme m, avec m allant de -! à 1, définie par: y"(O,O) = plIml(cos) trg m(O) avec: cos(mrd) pour m>O t pour m>0 trgm(Ob)= 1pour m=0 1 sin( m) pour m temporelle inverse des coefficients de Fourier-Bessel Pl,,m(/). Dans d'autres modes de réalisation, le champ acoustique est décomposé sur une base de fonctions, o chacune des fonctions s'exprime par une combinaison linéaire, éventuellement infinies de fonctions de FourierBessel. Sur la figure 2, on a représenté un organigramme général du procédé de l'invention. De manière générale, le procédé débute par une étape 2 de définition d'un traitement à appliquer à une représentation d'un champ acoustique et correspondant à un traitement au moins spatial de ce champ acoustique. Cette étape 2 de définition du traitement comporte une étape 4 d'établissement d'un ensemble de fonctions de directivité représentatif du traite15 ment. On entend par fonction de directivité une fonction, notée G(Of,ç4, qui associe une valeur réelle ou complexe, éventuellement fonction de la fréquence, ou d'une plage de fréquences, à chaque direction de l'espace. Toute fonction de directivité G(O,05f) se décompose en coefficients 20 d'harmoniques sphériques, notés Gl,m(P, par transformée de Fourier sphérique. Ainsi, la fonction de directivité G(09, V) s'obtient à partir des coefficients de Fourier sphérique Glm(f) par transformée de Fourier sphérique inverse, selon l'expression G(0,0,,f)=yEF Gj.(f) y, (0fJ/0l) 1=0 m=-l Dans le reste du document, les fonctions de directivité peuvent être désignées par leur expression G(9,ç dans l'espace direct ou sous forme de coefficients de Fourier sphérique Gl,m,(f). Cet ensemble de fonctions de directivité est délivré à l'issue de l'étape 4 sous la forme d'un ensemble de coefficients correspondant à la décomposition 30 des fonctions de directivité sur une base d'harmoniques sphériques. Par la suite, une étape 6 de détermination de combinaisons de filtrage est réalisée à partir de la décomposition des fonctions de directivité réalisée à l'étape 4. La notion de filtrage recouvre le filtrage fréquentiel, d'amplitude, de phase ou leur combinaison. Ainsi, l'étape 6 de détermination de combinaison de filtrage délivre des combinaisons de filtrage correspondant au traitement à appliquer sous la forme 5 de paramètres de filtres permettant ultérieurement de mettre en oeuvre ces combinaisons de filtrage à l'aide par exemple d'un ensemble de filtres numériques. En vue de l'application du traitement, le procédé comporte une étape 8 d'établissement d'une représentation initiale du champ acoustique sous la forme d'un ensemble de coefficients ou signaux représentatif de ce champ dans le 10 temps et dans les trois dimensions de l'espace. Notamment, cet ensemble de coefficients est un ensemble de coefficients dit de " Fourier-Bessel " tels que définis en référence à la figure 1. Le procédé comporte ensuite une étape 10 d'application des combinaisons de filtrage, définies lors de l'étape 2 de définition, aux coefficients de la 15 représentation initiale tels que définis lors de l'étape 8. Cette étape 10 permet d'obtenir une représentation modifiée correspondant à une représentation dans le temps et dans les trois dimensions de l'espace du champ acoustique initial modifié par le traitement défini, sous la forme d'un ensemble de coefficients dits de Fourier-Bessel. Sur la figure 3, on a représenté en détail l'organigramme d'un premier mode de réalisation de l'invention. Dans ce mode de réalisation, les traitements appliqués sont des traitements tels que le fenêtrage spatial consistant à modifier le gain de tout ou partie des sources virtuelles dont on considère qu'elles génèrent le champ acousti25 que initial, et la convolution spatiale ou filtrage spatial, consistant à modifier l'étalement angulaire de tout ou partie de ces sources virtuelles, ou tout autre traitement analogue. Ces traitements de fenêtrage et convolution peuvent varier avec la fréquence, ce qui permet d'effectuer des modifications d'éloignement et de spectre 30 de tout ou partie de ces sources virtuelles. Ainsi que cela a été décrit à la figure 2, le procédé débute par une étape 4 de détermination d'un ensemble de fonctions de directivité représentatif du traitement à appliquer. il Dans le cas o le traitement ne comporte qu'un seul traitement tel qu'un traitement de fenêtrage ou de convolution, une seule fonction de directivité, éventuellement variant avec la fréquence, suffit à le représenter. Avantageusement, l'étape 4 comporte une sous-étape 12 de 5 détermination directe de tout ou partie de cet ensemble de fonctions de directivité dans l'un ou l'autre d'une pluralité de formats prédéterminés. En outre ou alternativement, l'étape 4 comporte une sous-étape 14 d'entrée d'une consigne de traitement, laquelle est interprétée automatiquement lors d'une sous-étape 16 d'interprétation permettant de délivrer automatiquement 10 tout ou partie dudit ensemble de fonctions de directivité dans l'un ou l'autre d'une pluralité de formats prédéterminés. Par exemple, le traitement est mis en oeuvre par des moyens logiciels sur un ordinateur et cette consigne de traitement a la forme d'une icône dont la validation entraîne l'exécution par les moyens logiciels d'un traitement donné. Les sous-étapes 12 et 14 correspondent par exemple à une saisie manuelle ou encore à la réception de données de directivité provenant d'une mémoire interne ou d'un dispositif externe. En fonction du format dans lequel est déterminé ledit ensemble de fonctions de directivité, l'étape 4 peut requérir une sous-étape 18 de conversion 20 de l'ensemble de fonctions de directivité à partir d'un format prédéterminé en un ensemble de coefficients correspondant à la décomposition des fonctions de directivité sur une base d'harmoniques sphériques. Pour chaque fonction de directivité Glm,() de l'ensemble de fonctions de directivité, cette sous-étape 18 de conversion peut être réalisée de plusieurs 25 manières. Par exemple, la fonction de directivité Glm(f) est définie à partir de paramètres (rq,, Oq, bq) spécifiant la position de Q sources virtuelles et de paramètres Hq(/) spécifiant la réponse en fréquence de chacune de ces sources. A chaque source virtuelle est associé un modèle de rayonnement de sorte que l'ensemble 30 de ces sources est supposé produire un champ acoustique dont la représentation en coefficients de Fourier-Bessel permet d'associer à ce champ, une fonction de directivité. Par exemple, le modèle de rayonnement associé à chaque source est un modèle de rayonnement en ondes sphériques et la fonction de directivité Gl,m(f) est obtenue lors de la sous-étape 18 grâce aux équations suivantes: Q j2 2rq f Gi,.W H(f) e c M(rq,J) yjm(Oq, %) q=1 rq avec ( (/+k)! j2:órq f. i(rq = 2k! (1-k)! c En variante, ce calcul est réalisé à partir de paramètres (0q, @q) spécifiant les directions des Q sources. Dans ce cas, un modèle de rayonnement en ondes planes est associé à chaque source et la fonction de directivité Gl,m(f) est 10 obtenue grâce à l'équation suivante: Q Gi,.(f)= Hq (f) yP (Oq, q) q=l Dans d'autres modes de réalisation, GI,m(f) est déterminé à partir d'une combinaison de modèles de rayonnement en ondes planes et en ondes sphériques Dans un autre cas, la fonction de directivité est définie à partir de paramètres représentatifs d'une description paramétrique d'une fonction de directivité. Cette description s'effectue dans l'espace direct et associe à chaque direction une valeur pouvant être complexe. Avantageusement, les paramètres peuvent varier avec la fréquence. Dans ce cas, les coefficients de la fonction de directivité Gl,m(f) sont déterminés par la mise en oeuvre d'un modèle paramétrique. Les paramètres du modèle sont placés dans un vecteur;) et le modèle paramétrique de la fonction de directivité s'écrit G(O,,,/f)). Les coefficients de la fonction de directivité Gl,m(f) s'obtiennent par transformée de Fourier sphérique de G(O,q,(ji)) selon 25 l'équation: Igr 217r Gm(f) = jJ y G(O,(, X(f)) y, (0,O) sin dO do Par exemple, le modèle paramétrique peut correspondre à une fonction de directivité présentant un unique lobe. Ce modèle est contrôlé par les pa- ramètres (Oi(f),0f(f) spécifiant la direction du lobe et le paramètre 0(f) spécifiant l'angle d'ouverture du lobe. Ainsi, le modèle paramétrique peut être une fonction de directivité valant 0 pour les directions quelconques et valant 1 pour chaque direction (0,5) véri5 fiant la relation suivante: (010,*01 * (0,O) > CosOS(t" o (i(f),0(f)). (0,ó) désigne le produit scalaire entre les directions (0(0,),010) et (O,"5). Dans ce cas les coefficients GI,,m(f) de la fonction de directivité sont 10 donnés par: G/,m() = 2' y.]l (f1nq1,l (f)) P-,(cos0 (f) - P+I(cosoo(t) 21+1 Dans cette expression P.I(x) = 1. En variante, les paramètres de la description paramétrique de la fonction de directivité peuvent varier avec le temps. Dans ce cas, les paramètres du 15 modèle sont placés dans un vecteur (t) et le modèle paramétrique de la fonction de directivité s'écrit g(0,0, f(t)). Les coefficients de la fonction de directivité gl,m(t) exprimés dans le domaine temporel s'obtiennent alors par transformée de Fourier sphérique de g(0,0, f(t)) selon l'expression r7r 2ir gl,m(t) = j g(O,,'(t)) ym(0,) sin O dO d 20 Les coefficients Gl, m(f) de la fonction de directivité exprimés dans le domaine fréquentiel s'obtiennent alors par transformée de Fourier temporelle des coefficients gl,m(t). Dans encore un autre cas, la fonction de directivité est définie sous la 25 forme d'une liste de Q échantillons. La liste est formée de Q couples {((Oq,q), vq)} o (O,óq) est la direction de l'échantillon q et o Vq est la valeur de la fonction de directivité dans la direction (0q,óq). Avantageusement, la liste d'échantillons peut être saisie pour une pluralité de fréquencesf et est alors notée {((Oq,óq) , Vq)}00. Pour chaque fréquence f de fonctionnement, les coefficients Gl,m() de 30 la fonction de directivité sont déduits de la liste d'échantillons {((Oq, óq), vq)}0). Ces coefficients sont obtenus par inversion du processus d'échantillonnage angulaire qui permet de déduire les échantillons de la liste {((Oq,bq), vq)}(j) à partir d'une fonction de directivité fournie sous forme de coefficients d'harmoniques sphériques. Cette inversion peut prendre différentes formes afin de contrôler l'interpolation entre les échantillons. Dans le mode de réalisation décrit, les échantillons sont placés dans un vecteur Vde la façon suivante. [v, v2... Vq] Une matrice Y est constituée à partir des directions d'échantillonnage de la façon suivante: Yl,o,o YI,1,- 1 Yl,l,O YI,1,1... Y1,L,-L." Y1,L,O Y1,L,L 10 Y210,0 Y2,1,-l Y211,0 Y2zl,l.. Y2,L,-L '.Y2,,o Y2,L,L L YQ,o,oYQ,1,_1YQ,I,OYQ,1,1 YQ,L,-L YQ,L,o '.YQ,L,L o chaque élément Yq,i,m est égal à: Yq, Im=Yn( Oq,q) On rappelle que yr(Oq,oq) désigne l'harmonique sphérique d'ordre I et 15 de terme m pour la direction d'échantillonnage 0q,Oq considérée. On détermine un vecteur G contenant les coefficients Gl,m(j) pour la fréquencefs Le vecteur G peut être déterminé au moyen des relations suivantes: - si le nombre Q d'échantillons à la fréquence f est supérieur au nombre de coefficients Gl,m(f): G= ( yT y)-l yT V; et - si le nombre Q d'échantillons à la fréquence f est inférieur au nombre de coefficients Gl,m(): G= Y ( YY f)- V. Dans le vecteur G, les coefficients Gl,m(f) sont organisés de la manière 25 suivante: [G0,o() G,,_(f) G1,o(o) Gl,,(f)... GL,L(f)... GL,o(f) GLL()]1 En variante, la liste d'échantillons peut être saisie pour une pluralité d'instants t et est alors notée {((Oq, q), vq)}(t). Pour chaque instant t saisi, les coefficients gl,m(t) de la fonction de di30 rectivité exprimés dans le domaine temporel sont déduits de la liste d'échantillons {((fOq,,q), vq)}(t). Ces coefficients sont obtenus par inversion du processus d'échantillonnage angulaire qui permet de déduire les échantillons de la liste {((fOqoq), vq)}(t) à partir d'une fonction de directivité fournie sous forme de coefficients d'harmoniques sphériques. Les coefficients Glm(j) de la fonction de directivité exprimés dans le domaine fréquentiel s'obtiennent alors par transformée de Fourier temporelle des coefficients gl,m(t). Bien entendu, dans le cas o les étapes 12 ou 16 délivrent directement la fonction de directivité sous la forme de coefficients correspondant 10 à la décomposition de cette fonction sur une base d'harmoniques sphériques, la sous-étape 18 de conversion ne la modifie pas. Ainsi, à l'issue de l'étape 4, un ensemble de fonctions de directivité représentatif du traitement à appliquer est délivré sous la forme d'un ensemble de coefficients correspondant à la décomposition de ces fonctions sur une base 15 d'harmoniques sphériques. Par ailleurs, l'étape 2 de définition du traitement comporte également une étape 20 d'établissement d'un ensemble de paramètres du traitement. De même que dans le cas de l'étape 4, cet ensemble de paramètres peut être déterminé directement lors d'une sous-étape 22 de détermination di20 recte etiou à la suite d'une sous-étape 24 d'entrée d'une consigne de traitement et d'une sous-étape 26 d'interprétation de cette consigne de traitement permettant de délivrer automatiquement tout ou partie de l'ensemble de paramètres. Les sous-étapes 22 et 24 correspondent par exemple à une saisie manuelle de certains paramètres ou encore à la réception de données de para25 mètres à partir d'une mémoire interne ou d'un dispositif externe. Notamment, lors de l'étape 20 d'établissement d'un ensemble de paramètres, tout ou partie des paramètres suivants sont définis: - un paramètre L spécifiant l'ordre limite de fonctionnement du procédé et correspondant à la limite mathématique de précision souhaitée; et - un ensemble de constantes notées C[,,m formant des coefficients de pondération définis en fonction de la nature du traitement. En particulier, ces coefficients sont avantageusement appliqués aux fonctions de directivité pour obtenir les combinaisons de filtrage précitées. Certains paramètres peuvent être invariant et être directement intégrés dans le procédé. Le paramètre L permet au procédé de s'adapter au niveau de précision de la représentation initiale et permet notamment au procédé, objet de 5 l'invention, d'appliquer un traitement à une représentation initiale présentant un haut niveau de précision. Ainsi que cela a été décrit précédemment, les constantes C'rm" peuvent être déterminées directement lors de l'étape 20 par un opérateur, à partir d'une mémoire ou d'un dispositif interne. Par exemple, elles sont toutes calculées 10 préalablement, stockées dans une mémoire et sélectionnées en fonction de la nature du traitement. Alternativement, une consigne définissant la nature du traitement est saisie lors de la sous-étape 24 et la sous-étape 26 correspond au calcul de tout ou partie de cet ensemble de constantes en fonction de la nature du traitement et 15 à partir de relations décrites ultérieurement dans le texte. Bien entendu, les constantes cfm,'... peuvent être déterminées à l'aide d'une combinaison des solutions décrites. Les paramètres établis lors de l'étape 20, ainsi que les coefficients correspondants aux fonctions de directivité établis lors de l'étape 4, sont utilisés 20 lors de l'étape 6 de détermination des combinaisons de filtrage. Cette étape 6 de détermination comporte une sous-étape 30 de calcul de combinaisons de filtrage par des combinaisons linéaires des coefficients correspondants à la décomposition sur une base d'harmoniques sphériques de l'ensemble de fonctions de directivité représentatif du traitement. Lors de cette première sous-étape 30 de calcul, des combinaisons de filtrage, notées CI,:'(, sont déterminées à partir des fonctions de directivité Gi,m terme des coefficients de la représentation modifiée et 1", m" désignent l'ordre et le terme des coefficients de la fonction de directivité. Dans le cas o le traitement est un traitement de type fenêtrage, les constantes C[ç C,m,., sont déterminés par les relations: 1l m' 72(l', mIlml + CI,:-m,it,mi,) l"',m(1) ( cImI( il - l,"1) pour m > O pour m = O pour m l',m' (1) Cllm,l,,mI - 1/ ',m' (2) ',lm' (2) I 2 ( Clml,,,m", + c'l,m,l",-m"! l ',m' (2) CIm,l",, O | ( l'm'(2) - cl ',m' (2) & _ \Iml",-m" -CI,m,l",m") pour m" > O pour m": O pour m" O pour m' = O pour m' l',m' (3) Cl m,l,m - Ce traitement correspond à l'application d'un fenêtrage spatial avec filtrage temporel. En effet, les coefficients d'harmoniques sphériques Gi, m() corres20 pondent à la fonction de directivité G(0,ç). De même, les coefficients de Fourier- Bessel P()i,m(f) du champ acoustique initial correspondent à des coefficients d'harmoniques sphériques d'une fonction de directivité P()(0, ó,). Ainsi, le fenêtrage spatial avec filtrage temporel se traduit, pour chaque fréquence f par un fenêtrage de la fonction de directivité P()(1, çJ) par la fonction de directivité 5 G(O,j), ce qui correspond à une multiplication pour chaque direction (0,) de la valeur de ces deux fonctions de directivité, les valeurs de ces deux fonctions de directivité pouvant être complexes. Ainsi que cela a été décrit précédemment, ce traitement appliqué au champ acoustique a pour effet de modifier l'amplitude, l'éloignement et le spectre 10 de tout ou partie de l'ensemble des sources virtuelles dont on considère qu'elles génèrent le champ acoustique. La fonction de directivité G(O,ó) peut être indépendante de la fréquence et est alors noté G(O,). Le traitement est alors uniquement spatial et correspond, pour chaque fréquences, à la multiplication de la fonction de directi15 vité P( (O,,/j) par la fonction de directivité G(0,). Dans le cas o le traitement est un traitement de type convolution, les constantes C i',mV, sont déterminées par la relation ( 1 f ',m (1) _ ', n'(1) r ' (Clm,'(,m" + Cl' m(n>mî) pour m > O C -,,m C01'r'(m) pour m = 0 t ', (C>l-m,)lm"r -C m,-lfm) pour m O l',m' (1) l'_ m' (2) p Cl mr 1"J rma cl,m,lI,0 pour m = O In (Cl m1'$ m" -C m lrn(2)) ppour m O C lml(2)",= Cl,' , (3) pour m' = O t () l', m' (3) -C,m' (3) * l cM'itfM - {rn;,$3),) pourm' CP+ O _11P+k ± (!P (".CP+M;++ -l+m Il-Mil'" X I /C'p+k+mn En effet, les coefficients d'harmoniques sphériques Gl,m() correspondent à des coefficients de Fourier-Bessel d'un champ acoustique produit par un système acoustique linéaire et invariant angulairement, en réponse à une ex10 citation par un champ acoustique de référence tel qu'une onde plane impulsionnelle provenant de la direction (0,0). Une telle réponse est désignée par l'expression " réponse spatio-temporelle ". Un système acoustique linéaire est invariant angulairement lorsque qu'une rotation du champ qu'il reçoit induit la même rotation sur le champ qu'il génère. Ainsi, les coefficients Gum(j) corres15 pondent à la réponse spatiotemporelle d'un tel système acoustique. Les coefficients de Fourier-Bessel P(0I,m(/) du champ acoustique initial sont également considérés comme des coefficients d'harmoniques sphériques d'une fonction de directivité P( (0,bj). Ainsi, la convolution spatiale avec filtrage temporel se traduit, pour chaque fréquence f par une convolution spatiale de la 20 fonction de directivité P(")(O,óO par la fonction de directivité G(O,,), qui est considérée à chaque fréquencefcomme un noyau de convolution directif. Dans le mode de réalisation décrit, la convolution est appliquée en considérant une direction (00,ó), le noyau de convolution G(0,ód) est alors tourné en amenant son pôle nord (direction verticale avec 0 = O et 0 = 0) en (00,0"). Pour 25 obtenir la valeur du champ acoustique modifié P(l(O,ó,) dans la direction (O','), la valeur dans la direction (O',0') du noyau tourné est multipliée avec la valeur de P('i(O,çO dans la direction (0o,ó0) et ces produits sont sommés pour toutes les directions (60,ó0). Ce traitement appliqué au champ acoustique permet, par exemple, de 30 modifier l'étalement angulaire, soit la précision de localisation spatiale, c'est-à- dire la ponctualité des sources virtuelles, dont on considère qu'elles génèrent le champ acoustique initial, pour chaque fréquence f. Ce traitement est l'analogue spatial de la convolution temporelle par une fenêtre, qui a pour effet de modifier l'étalement temporel d'un événement sonore, comme par exemple un effet de moyennage ou de retard. Dans le cas o la fonction de directivité Glm(t) est à symétrie cylindrique, c'est-àdire si Gl,m(I)=0 pour mÉO, les équations ci-dessus se simplifient, et les coefficients c;,m. sont alors déterminés par la relation:,',m' _ 1 l m O /4 cl=, mi' =m - 61 m,, V21 + 1 10 La fonction de directivité G(O,A) peut être indépendante de la fréquence et est alors noté G(09,'). Le traitement est alors uniquement spatial et correspond, pour chaque fréquencef, à la convolution de la fonction de directivité P('(O,OJ) associée au champ initial par la fonction de directivité G(O,ç). Par ailleurs, le procédé de l'invention permet également la réalisation 15 de traitements de rotation. Dans ce cas cependant, il n'est pas fait usage de fonctions de directivité et le procédé requiert des paramètres (0,0,q,), représentatifs d'une rotation du champ acoustique autour de l'origine O. Ainsi, si le traitement est une rotation, les combinaisons de filtrage sont déterminées à partir des paramètres (O0,V,/) représentatifs de la rotation du 20 champ acoustique selon la relation: C,', ' I' R C/,m _ 6 D-m,m avec 8R(D I, + Dl_',m) pour m' > O et m > O V2R(DMo) pour m' > O et m O (D/,m + Dlm,m) pour m'> O et m O DM,m= DO,0 pour m' = Oetm= O (D'om) (Dpour m'= 0O et m O -x/--q9(DOm,o) pour m' Dans cette relation la convention utilisée pour les angles 0, 0 et sf', 10 consiste à effectuer une rotation d'un angle V' autour de l'axe (OZ) puis d'un angle O autour de l'axe (OY) puis d'un angle ó autour de l'axe (OZ). Enfin, le traitement effectué peut également correspondre à un traitement neutre. Dans ce cas, les combinaisons de filtrage sont toutes nulles, sauf pour les valeurs o l'= 1 et m'= m o elles valent 1. Bien entendu, d'autres traitements ainsi que des combinaisons de traitements décrits peuvent être réalisés grâce à la définition de l'ensemble de constantes cf,'.' de pondération. Par exemple, un traitement aléatoire peut être réalisé par la définition aléatoire de l'ensemble de constantes. Avantageusement, la sous-étape 30 de calcul est conduite uniquement pour les constantes C,m. non nulles. Les combinaisons de filtrage CGr;n'(/) ainsi obtenues sont placées dans une matrice C de taille (L'+1)2 sur (L+1)2. Dans le cas o le traitement à appliquer consiste uniquement en une 25 rotation, une convolution, ou un filtrage, la matrice C correspond directement aux combinaisons de filtrage à appliquer et est notée de manière générale F de sorte que les résultats de la sous-étape 30 de calcul forment les combinaisons de filtrage destinées à être appliquées aux coefficients de la représentation initiale. Dans le cas o le traitement correspond à une combinaison de traite30 ments tels que des rotations, filtrages, convolutions etiou traitements quel- conques, la sous-étape 30 est répétée une pluralité de fois et chaque matrice C est successivement mémorisée afin d'être combinée avec les matrices suivantes pour finalement délivrer une matrice F comprenant les coefficients correspondant aux combinaisons de filtrage destinées à être appliquées aux coefficients de la 5 représentation initiale. Ainsi, la matrice F est obtenue par multiplication matricielle des matrices C successives. Chaque élément F,/m(I) de la matrice F représente la réponse en fréquence d'un filtre de traitement spatio-temporel. Différentes formes de filtres et de méthodes de filtrage peuvent être 10 mises en oeuvre lors de l'étape 10 d'application des combinaisons de filtrage à la représentation initiale du champ acoustique. Dans le cas o les filtres sont paramétrés directement avec des réponses en fréquence, les coefficients Flm''() sont directement appliqués. Avantageusement, l'étape 6 de détermination des combinaisons de filtrage comprend 15 une sous-étape 32 de conversion afin de déterminer les paramètres des filtres pour d'autres méthodes de filtrage. Par exemple, les combinaisons de filtrage F,,rm'(j) sont converties en: des réponses impulsionnelles finies f1tmm'(t) calculées par transformée de Fourier temporelle inverse de Fzm'(î), chaque réponse impulsionnelle 20 flIm'(t) est échantillonnée puis tronquée à une longueur propre à chaque réponse; ou - des coefficients de filtres récursifs à réponses impulsionnelles infinies calculées à partir des F,m'(j) avec des méthodes d'adaptation. Les combinaisons de filtrage ainsi définies sont appliquées lors de 25 l'étape 10 d'application des combinaisons de filtrage laquelle transforme la représentation initiale P(im( en une représentation modifiée notée P(7Im représentatifs du traitement au moins spatial défini lors de l'étape 2. Lors de cette étape 10, des coefficients p(7),m(t) correspondant aux 30 coefficients de la représentation modifiée P('), sont déduits des coefficients p(l)',m(t) correspondant aux coefficients de la représentation initiale P(), par l'application des filtres de traitement de réponse en fréquence F/,,m'' de la façon suivante: L 1 1=0 m=-I o P('I,M(f) est la transformée de Fourier de p( l),m(t) et o P(>lm0 est la transfor5 mée de Fourier de p(nl,m(t), coefficients relatifs à la représentation initiale et respectivement à la représentation modifiée. Ainsi, l'étape 10 met en oeuvre les méthodes classiques de filtrage temporel pour l'application des combinaisons de filtrages délivrées par l'étape 2 dans une forme adaptée. Par exemple: - Le filtrage peut être effectué au moyen de méthodes de filtrage dans le domaine fréquentiel, comme par exemple des techniques de convolution par bloc. Dans ce cas les combinaisons de filtrage sont fournies sous la forme des réponses en fréquences F,,m'(/). - Le filtrage peut être effectué dans le domaine temporel par convolution temporelle. Dans ce cas, les combinaisons de filtrage sont fournies sous la forme des réponses impulsionnelles f,,m'(t). - Le filtrage peut être réalisé dans le domaine temporel au moyen de filtres récursifs à réponse impulsionnelle infinie. Dans ce cas, les combinaisons 20 de filtrage sont fournies sous la forme des coefficients des relations de récurrence. Sur la figure 4, on a décrit le cas d'un filtrage par réponse impulsionnelle finie. Ce filtrage requiert initialement la détermination d'une variable Tr,' correspondant au nombre d'échantillons propres à chaque réponse f*,7m(t), ce 25 qui conduit à l'expression de convolution temporelle suivante: L I Tm I P:[t] =E j I'm[ r] plE[t-r] 1=0 m=-I r=O Ces coefficients p('î,m(t) sont représentatifs d'un champ acoustique correspondant au champ acoustique initial après l'application du traitement au moins spatial qui associe à un champ acoustique produit par au moins une 30 source, un autre champ acoustique produit par au moins une source dont la direction, l'éloignement et/ou le gain sont ainsi modifiés. Dans un autre mode de réalisation, le procédé de l'invention permet en outre l'application de traitements de type distorsion correspondant à une distorsion de parties du champ acoustique par des modifications des caractéristiques d'étalement, de position, d'amplitude de tout ou partie des sources virtuelles gé5 nératrices du champ acoustique initial. Ce traitement de distorsion peut varier avec la fréquence, ce qui permet d'effectuer des modifications d'éloignement et de spectre de tout ou partie de ces sources virtuelles. Avantageusement, ce traitement peut être combiné à un des traite10 ments tels que décrits en référence à la figure 3. Sur la figure 5, on a décrit l'organigramme du procédé de l'invention dans un tel mode de réalisation combinant un premier traitement appelé traitement nominal, correspondant à un traitement tel que décrit précédemment, avec un traitement de distorsion spatiale qui se traduit par des distorsions de caracté15 ristiques de fréquence, d'amplitude et/ou de phase de composants de la représentation du champ acoustique initial. Ainsi, dans ce mode de réalisation, l'étape 2 de définition du traitement comprend, de même que précédemment, l'étape 4 d'établissement d'un ensemble de fonctions de directivité comportant la sous-étape 12 de détermination di20 recte et/ou les sous-étapes 14 de saisie et 16 d'interprétation ainsi que la sousétape 18 de conversion permettant de délivrer un ensemble de fonctions de directivité sous la forme d'un ensemble de coefficients correspondants à la décomposition de ces fonctions sur une base d'harmoniques sphériques. Dans ce mode de réalisation, une partie de l'ensemble de fonctions de 25 directivité déterminé lors de l'étape 4 correspond à N couples de distorsion représentatifs d'une distorsion à appliquer. De manière similaire au mode de réalisation décrit précédemment, ces couples de distorsion sont décrits sous la forme d'un ensemble de coefficients, notés {(Himn décomposition sur une base d'harmoniques sphériques. Chaque couple de distorsion est constitué de deux fonctions de directivité qui correspondent aux coefficients de Fourier-Bessel d'un champ acoustique source et d'un champ acoustique cible, l'opération de distorsion appliqué au champ acoustique source tendant à le faire ressembler au champ acoustique cible. Ainsi, l'étape 4 délivre un ensemble de fonctions de directivité correspondant au traitement nominal ainsi qu'à des couples de distorsion représentatifs de la distorsion à appliquer sous la forme d'un ensemble de coefficients correspondant à la décomposition des fonctions de directivité sur une base d'harmoniques sphériques. De même que précédemment, l'étape 2 comporte également l'étape de détermination de l'ensemble de paramètres du traitement par le biais de la sous-étape 22 de détermination directe et/ou des sous-étapes 24 de saisie et 26 d'interprétation. Toutefois, dans ce mode de réalisation, l'étape 20 de détermination de paramètres permet en outre la détermination d'un paramètre 2, pouvant être ajusté entre O et 1, spécifiant la ressemblance entre le traitement final à appliquer et le traitement nominal déterminé par l'étape 4, soit encore l'importance relative de la distorsion à effectuer. Ainsi, si 2 est égal à 1, seule la distorsion est prise en compte et si est égal à 0, seul le traitement déterminé lors de l'étape 4 est pris en compte. De même que dans le mode de réalisation décrit en référence à la figure 3, l'étape 6 de détermination des combinaisons de filtrage comprend la sousétape 30 de calcul réalisée ainsi que décrit précédemment et permettant de 20 délivrer une matrice C comprenant un ensemble de combinaisons de filtrage correspondant à l'application du traitement nominal tel que défini par l'ensemble de fonctions de directivité déterminé à l'étape 4. L'ensemble de couples de distorsion est utilisé dans l'étape 6 de détermination des combinaisons de filtrage lors d'une sous-étape 50 de distorsion 25 par modification des combinaisons de filtrage déterminées à l'issue de la sousétape 30 de calcul. Lors de cette sous-étape 50, les combinaisons de filtrage Ch,In(j), sont modifiées à partir des couples de distorsion notés (Hnt)H,,)}et du paramètre afin de fournir en sortie les combinaisons de filtrage Fr,m(f) représentati30 ves du traitement incluant le traitement nominal et la distorsion et destinées à être appliquées aux coefficients de la représentation initiale. Lors de la sous-étape 50, les coefficients d'harmoniques sphériques Hi,m, n(/ représentatifs des N fonctions de directivité sources sont placés dans une matrice H de taille (L 1)2 sur N de la façon suivante: Ho,o,,(f) Ho,o,2(f)...... Hoo,N(f) H1,-l,1() Hl,-1,2()...... Hl,-i,N(J) Hloa(f) Hl2o,2(f)...... Hl,oN(f) Hl,.l H1,1,2(f)..... H,"N(f) HL,-Ll()HL,-L,2()... .HL,-L,N(f) HL,o, (f) HL,o,2()...... HL,o,N() HL,L,(,) HL,L,2(f)...... HL,L,N(f) De même, les coefficients d'harmoniques sphériques H'lmn() représentatifs des Nfonctions de directivité cibles sont placés dans une matrice I/ de taille (L'+1)2 sur N: H'o,,o(f) H'0,0,2(f)...... H'ooN(f) H'L-1 j(f) H'2,-2,2f)...... H',,,, -N (f) H',,o,1(J) H'1,o,2(f)...... H'l, ,N (f) m'1,,2 (f) ',2()......H'L,,L,N (f) m 'L -.t,, (f')Hm L',-.L',,2 (|) *).....',-.z,,N (f) H'L,Po,() H'L',O,2(f)...... H'L',o,N(f) H'L.,L.,I(f) H'L,,L',2(f)...... H'L',L,N(f) Les coefficients CIm'(/) des combinaisons de filtrage du traitement nominal sont également placés dans une matrice C de taille (L'+1)2 sur (L+1)2: 0C0 0 l' 1 l o C '0 Co,i 0 Cl,-ll_ Cll(o o ci -") r10{ 1,0(j {n 1,0 c lo, c0{ff Co, 1 -I C 'o 1i9 :0Of 0 Cl'll(lf c '. C.0.3..DTD: 220 10 C 2,1C I L' nCL',-L' Co:,O(f Cl',O CL'OU) CL',O_) 0,0 l'- 1 lC'o() l,'! C L '-L' C L' -L' L'-L' L-L' 0,0 ','- C,: 09 l"' '' C-0) C... o... o o CL'_ L.o.. LO0 C L:L03 CL'MV c CL,0O C, C L -L l' * CL,20... LL C ci,_ou)... cioffJ L' L C;o(f 'O C zL L 11 2 1 cL',-L'ff c'LC,-. 03'-'CL,L. 0 C L'" f.. L,Ofl.CL, L,-L% À.. '.. CL' z' Cz, 03'"Lfl. C L, '(f)l L, L L 0 L L' Les combinaisons de filtrage Frm'() représentatives du traitement dans son ensemble sont ensuite déterminées au cours de la sous-étape 50 de 15 distorsion correspondant à une optimisation qui détermine un ensemble de combinaisons de filtrage vérifiant au mieux les couples de distorsion et la ressemblance entre F,,m'() et Ctm'(p. lm V2lm V Cette optimisation peut être réalisée au moyen de méthodes telles que des méthodes de descente de gradient ou d'inversion matricielle et permet également de déterminer, sur les sources virtuelles dont on considère qu'elles génèrent le champ acoustique initial, des modifications qui ne sont pas directement spécifiées par les couples de distorsion. Par exemple, les modifications définies par les couples de distorsion sont interpolées sur les parties voisines des parties directement visées de manière à obtenir les combinaisons de filtrage, ce qui se traduit par des modifications progressives de la représentation initiale du champ acoustique autour des 10 parties directement visées par les couples de distorsion. Le mode de réalisation décrit utilise une méthode d'inversion matricielle qui conduit à la détermination d'une matrice F lors de la sous-étape 50 selon la relation: F=((1-2) C+ fH'HT)((l-) I+ AHHT)15 o Idésigne la matrice identité de taille (L+1)2. * Dans le cas o 2 vaut 1 et o N des modifications entre les parties directement visées par les couples de distorsion et les parties voisines. Par la suite, les combinaisons de filtrage F/,'(1) sont éventuellement converties lors de la sous-étape 32 et sont appliquées lors de l'étape 10 ainsi que cela a été décrit en référence aux figures 3 et 4. Avantageusement, plusieurs distorsions sont appliquées successive5 ment, le procédé comportant alors une pluralité de sous-étapes 50 de distorsion, récursives les unes par rapport aux autres, de sorte que chaque sousétape 50 reçoit en entrée, comme traitement nominal, les combinaisons de filtrage délivrées par la sous-étape précédente. Dans le cas o le paramètre A vaut 0, la distorsion n'est pas appliquée 10 et les combinaisons des filtrages déterminées lors de la sous-étape de calcul 30 correspondent directement aux combinaisons de filtrage représentatives du traitement. Ainsi paramétré, ce mode de réalisation correspond au premier mode de réalisation décrit en référence aux figures 3 et 4. Dans un autre mode de réalisation, seule la distorsion est appliquée et 15 les combinaisons de filtrages sont déterminées uniquement au moyen des couples de distorsion selon la relation: -SiN F=I + (IF -H) (HT H)' HT -SiNÄ(L+1)2: F=H HT (HHT) Dans ces équations, les matrices H et If sont déterminées ainsi que cela a été décrit précédemment. Dans ce cas, la matrice C n'est pas utilisée et la sous-étape 30 n'est pas réalisée. Ainsi, l'étape 2 comporte uniquement la sousétape 50 de distorsion. Ce mode de réalisation est une simplification du mode de réalisation décrit lorsque A vaut 1 et que le traitement nominal est le traitement neutre. Il apparaît donc que le procédé de l'invention permet d'appliquer des traitements quelconques, présentant une très large diversité à la représentation d'un champ acoustique présentant un haut niveau de précision grâce à 30 l'utilisation d'un ensemble de fonctions de directivité représentatif du traitement et représenté sous la forme d'un ensemble de coefficients correspondant à la décomposition de ces fonctions sur une base d'harmoniques sphériques. Dans le cas o le traitement à appliquer peut se décomposer en une pluralité de traitements élémentaires tels que décrits précédemment, les combi- naisons de filtrage du traitement dans son ensemble peuvent être obtenues par multiplication matricielle des matrices correspondant à chacun de la pluralité de traitements, soit encore par arrangement des combinaisons de filtrage délivrées par la sous-étape 30 de calcul et/ou la sous-étape 50 de distorsion. Ce traitement est ainsi répété une pluralité de fois, les résultats de chaque occurrence étant combinés avec ceux de l'occurrence précédente. Sur la figure 6, on a représenté un schéma synoptique d'un dispositif adapté pour la mise en couvre du procédé de l'invention. Ce dispositif 60 permet d'effectuer un traitement au moins spatial 10 d'une représentation d'un champ acoustique afin de délivrer une représentation modifiée par le traitement au moins spatial. Le dispositif 60 comporte des moyens 61 de définition du traitement pour la mise en oeuvre de l'étape 2 du procédé de l'invention. Ces moyens 61 comportent des moyens d'établissement d'un ensemble de fonctions de directivi15 té représentatif du traitement, sous la forme d'un ensemble de coefficients correspondant à la décomposition de ces fonctions de directivité sur une base d'harmoniques sphériques. Avantageusement, ils comportent un module 62 de saisie directe, par un opérateur, à partir d'une mémoire ou d'un dispositif extérieur, d'informations 20 de fonctionnement du dispositif 60 permettant notamment la détermination directe de l'ensemble de fonctions de directivité. Par exemple, ce module 62 de saisie est constitué d'un clavier permettant la saisie directe d'informations ou encore d'une interface avec un dispositif extérieur fournissant directement les informations nécessaires à la mise en ceu25 vre du procédé de l'invention. Le module 62 peut également être mis en couvre pour l'entrée d'une consigne de traitement, les moyens 61 comportent alors également un module 64 d'interprétation automatique de cette consigne, afin de délivrer des informations de fonctionnement du dispositif 60 et notamment l'ensemble de fonctions de di30 rectivité. Les informations entrées grâce au module 62 de saisie et/ou 64 d'interprétation, peuvent en outre être de différents formats prédéterminés et le dispositif 60 comporte un module 66 de conversion des informations entre des formats prédéterminés et des formats sélectionnés pour la mise en oeuvre du procédé de l'invention. Notamment, le module 66 permet de convertir ledit ensemble de fonctions de directivité de l'un ou l'autre d'une pluralité de formats en un ensemble de 5 coefficients correspondant à leur décomposition sur une base d'harmoniques sphériques. Ainsi, les modules 62, 64 et 66 permettent la mise en oeuvre des étapes 4 et 20 du procédé de l'invention, de manière à délivrer un ensemble de fonctions de directivité représentatif du traitement sous la forme des coefficients 10 correspondant à la décomposition de ces fonctions sur une base d'harmoniques sphériques. Dans le cas o le dispositif 60 est utilisé pour la réalisation d'un traitement de distorsion, les modules 62, 64 et 66 délivrent des fonctions de directivité formant des couples de distorsion. Les modules 62, 64 et 66 permettent également la détermination de tout ou partie d'un ensemble de paramètres nécessaire à la réalisation du procédé de l'invention. Notamment, ces modules 62, 64 et 66 permettent la définition de tout ou partie des paramètre suivants: - le paramètre L spécifiant l'ordre limite de fonctionnement du procédé et correspondant à la limite mathématique de précision souhaitée; - un paramètre S définissant la nature du traitement à effectuer; - tout ou partie de l'ensemble de constantes notées Citm formant des coefficients de pondération et définies en fonction de la nature du traitement; 25 et - le paramètre A, spécifiant l'importance d'une distorsion à effectuer. Avantageusement, le dispositif 60 comporte une mémoire permettant de stocker tout ou partie de ces paramètres et notamment l'ensemble de constantes. Par exemple, une telle mémoire comporte une pluralité d'ensembles de 30 constantes, chacun correspondant à un traitement, l'ensemble de constantes mis en oeuvre étant sélectionné dans cette pluralité en fonction du paramètre S représentatif du traitement. Alternativement, les moyens 61 de définition du traitement comprennent un module de calcul adapté pour déterminer tout ou partie de l'ensemble de constantes en fonction du paramètre S représentatif de la nature du traitement. L'ensemble de ces informations de fonctionnement est introduit dans 5 des moyens 70 de détermination de combinaisons de filtrage, compris dans les moyens 61 de définition du traitement, et adaptés pour la mise en oeuvre de l'étape 6 du procédé. Ces moyens 70 comportent un module 72 de calcul par application de combinaisons linéaires permettant la mise en oeuvre de la sousétape 30 de calcul du procédé. Ce module 72 reçoit en entrée l'ensemble de 10 fonctions de directivité afin de délivrer des combinaisons de filtrage par combinaisons linéaires des coefficients correspondants à la décomposition sur une base d'harmoniques sphériques des ensembles de fonctions de directivité représentatives du traitement. Avantageusement, les constantes délivrées par les modules 62, 64 et 15 66, forment des coefficients de pondération et sont utilisées comme coefficients des combinaisons linéaires dans le module 72. Dans le cas o le dispositif 60 met également en oeuvre un traitement de distorsion, il comporte également un module 74 de distorsion par modification des combinaisons de filtrage délivrées par le module 72, mettant en oeuvre la 20 sous-étape 50 du procédé. Ce module reçoit en entrée les couples de distorsion représentatifs de la distorsion à appliquer, ainsi que le paramètre A représentatif de l'effort de distorsion mis en oeuvre ainsi que cela a été décrit précédemment. Dans le cas o seule une distorsion est appliquée, le module 72 n'intervient pas dans la mise en oeuvre du dispositif et uniquement le module 74 25 de distorsion est nécessaire. Ainsi, les moyens 70 délivrent un ensemble de coefficients correspondant à des combinaisons de filtrage représentatives du traitement à appliquer. Eventuellement, ces combinaisons de filtrage sont introduites dans un module 76 de conversion adapté pour délivrer des paramètres de filtres à partir 30 des combinaisons de filtrage en mettant en oeuvre la sous-étape 32 du procédé ainsi que cela est décrit précédemment. Ces paramètres de filtres ou directement les combinaisons de filtrage sont appliqués à des moyens de filtrage 80 qui mettent en oeuvre l'étape 10 du procédé et reçoivent en entrée la représentation initiale du champ acoustique sous forme d'un ensemble de coefficients représentatif de ce champ dans le temps et dans les trois dimensions de l'espace pour appliquer les filtres déterminés afin de délivrer une représentation modifiée correspondant au champ acoustique initial modifié par le traitement défini précédemment. Par exemple, les moyens de filtrage 80 sont des moyens de filtrage numériques appliqués par un processeur de traitement du signal à une représentation numérique du champ acoustique initial. Dans d'autres modes de réalisation, ces moyens de filtrage sont des moyens de filtrage analogiques à gain commandé, appliqués à une représenta1 0 tion analogique du champ acoustique sous la forme d'un ensemble de signaux correspondant chacun à un coefficient de la représentation initiale du champ acoustique. Par exemple, le dispositif 60 peut être réalisé sous la forme d'un ordinateur équipé de moyens logiciels de calcul et de moyens matériels d'acquisition 1 5 et de restitution, de manière à permettre la mise en oeuvre des moyens 62 d'entrée d'informations ainsi que l'acquisition de la représentation initiale du champ acoustique et la sortie de la représentation modifiée du champ acoustique, les différents éléments du dispositif 60 étant réalisés par des combinaisons de moyens logiciels et de moyens matériels. Dans d'autres modes de réalisation, le dispositif 60 correspond à une carte d'acquisition spécifique dédiée uniquement à la réalisation de traitements acoustiques ou encore sous la forme d'un composant programmé pour la mise en oeuvre de ces fonctions uniquement. Avantageusement, le dispositif 60 est intégré dans un système acous25 tique comprenant des moyens d'acquisition d'un champ acoustique et des moyens de représentation de ce champ acoustique sous la forme d'un ensemble de coefficients dans le temps et dans les trois dimensions de l'espace afin de délivrer la représentation initiale du champ acoustique. Un tel système comporte également un ensemble de restitution commandé par des moyens de pilotage 30 recevant en entrée, la représentation modifiée du champ acoustique et pilotant un ensemble de restitution afin d'obtenir la restitution de ce champ acoustique modifié. REVENDICATIONS 1. Procédé de traitement d'une représentation d'un champ acoustique, afin de délivrer une représentation modifiée par application de ce traitement correspondant à un traitement au moins spatial du champ acoustique, compor5 tant: - l'établissement (8) d'une représentation initiale du champ acoustique sous la forme d'un ensemble de coefficients représentatifs de ce champ dans le temps et dans les trois dimensions de l'espace; et l'établissement (10) d'une représentation modifiée par des combinai10 sons de filtrage correspondant audit traitement, appliquées auxdits coefficients de ladite représentation initiale, caractérisé en ce qu'il comporte en outre une étape (2) de définition dudit traitement, laquelle comprend: - une étape (4) d'établissement d'un ensemble de fonctions de directi15 vité représentatif dudit traitement, sous la forme d'un ensemble de coefficients correspondant à la décomposition de ces fonctions sur une base d'harmoniques sphériques; et - une étape (6) de détermination desdites combinaisons de filtrage à partir de cette décomposition pour la mise en oeuvre d'une étape (10) 20 d'application des combinaisons de filtrage ainsi définies auxdits coefficients de la représentation initiale, afin d'obtenir ladite représentation modifiée. 2. Procédé selon la revendication 1, caractérisé en ce que ladite étape (4) de détermination d'un ensemble de fonctions de directivité comporte une sous-étape (12) de détermination directe de tout ou partie dudit ensemble de 25 fonctions de directivité dans l'un ou l'autre d'une pluralité de formats prédéterminés. 3. Procédé selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que ladite étape (4) de détermination d'un ensemble de fonctions de directivité, comporte une sous-étape (16) d'interprétation automatique d'une 30 consigne de traitement permettant de délivrer automatiquement tout ou partie dudit ensemble de fonctions de directivité dans l'un ou l'autre d'une pluralité de formats prédéterminés, à partir de ladite consigne de traitement. 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que ledit ensemble de fonctions de directivité représentatif du traitement est exprimé selon l'un ou l'autre d'une pluralité de formats prédéterminés et en ce que ladite étape (4) d'établissement d'un ensemble de fonctions de directivité comporte une sous-étape (18) de conversion, à partir dudit format prédéterminé, en un ensemble de coefficients correspondant à la décomposition desdites fonc5 fions de directivité sur une base d'harmoniques sphériques. 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que ladite étape de définition (2) du traitement comprend en outre, une étape (20) d'établissement d'un ensemble de paramètres dudit traitement. 6. Procédé selon la revendication 5, caractérisé en ce que ledit en10 semble de paramètres comprend tout ou partie des paramètres suivants: - un ensemble de constantes (c - un paramètre (L) spécifiant l'ordre limite de fonctionnement du procédé et correspondant à la limite mathématique de précision souhaitée. 7. Procédé selon l'une quelconque des revendications 5 ou 6, caractérisé en ce que ladite étape (20) d'établissement d'un ensemble de paramètres comporte une sous-étape (22) de détermination directe de tout ou partie dudit ensemble de paramètres. 8. Procédé selon l'une quelconque des revendications 5 à 7, caractéri20 sé en ce que ladite étape (20) d'établissement d'un ensemble de paramètres comporte une sous-étape (26) d'interprétation automatique d'une consigne de traitement permettant de délivrer automatiquement tout ou partie dudit ensemble de paramètres à partir de ladite consigne de traitement. 9. Procédé selon l'une quelconque des revendications 6 à 8, caractéri25 sé en ce que ladite étape (20) d'établissement d'un ensemble de paramètres comprend une sous-étape de calcul de tout ou partie dudit ensemble de constantes en fonction de la nature dudit traitement. 10. Procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce que ladite étape (6) de détermination desdites combinaisons de filtrage 30 comprend au moins une sous-étape (30) de calcul de combinaisons de filtrage par des combinaisons linéaires desdits coefficients correspondant à la décomposition sur une base d'harmoniques sphériques dudit ensemble de fonctions de directivité représentatif dudit traitement. 11. Procédé selon les revendications 6 et 10 prises ensemble, caractérisé en ce que ladite sous-étape (30) de calcul par des combinaisons linéaires fait intervenir ledit ensemble de constantes en tant que coefficients de pondération. 12. Procédé selon l'une quelconque des revendications 10 ou 11, ca5 ractérisé en ce que le résultat de ladite sous-étape (30) de calcul forme directement lesdites combinaisons de filtrage destinées à être appliquées auxdits coefficients de la représentation initiale. 13. Procédé selon la revendication 12, caractérisé en ce que ledit traitement est un traitement de convolution. 14. Procédé selon la revendication 12, caractérisé en ce que ledit traitement est un traitement de fenêtrage. 15. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce que ledit traitement comprend au moins une distorsion à appliquer à la représentation initiale du champ acoustique et en ce que tout ou partie des 15 fonctions de directivité déterminées lors de l'étape (4) d'établissement d'un ensemble de fonctions de directivité forment un ensemble de couples de distorsion représentatif de ladite distorsion sous la forme d'un ensemble de coefficients correspondant à la décomposition desdits couples de distorsion sur une base d'harmoniques sphériques. 16. Procédé selon les revendications 10 et 15 prises ensemble, caractérisé en ce que ledit traitement comprend un premier traitement dit nominal et la ou lesdites distorsions, ladite sous-étape (30) de calcul délivrant des combinaisons de filtrage correspondant audit traitement nominal et ladite étape (6) de détermination de combinaisons de filtrage comprenant en outre au moins une sous25 étape (50) de distorsion par modification des combinaisons de filtrage délivrées par ladite sous-étape (30) de calcul, à partir desdits couples de distorsion représentatifs de la ou desdites distorsions, afin de délivrer lesdites combinaisons de filtrage destinées à être appliquées auxdits coefficients de la représentation initiale. 17. Procédé selon les revendications 5 et 16 prises ensemble, caractérisé en ce que ledit ensemble de paramètres comporte un paramètre (X.) représentatif de l'effort de distorsion, soit encore spécifiant la ressemblance entre ledit traitement à appliquer et ledit traitement nominal. 18. Procédé selon la revendication 15, caractérisé en ce que ledit traitement comprend uniquement la ou lesdites distorsions, ladite étape (6) de détermination des combinaisons de filtrage comprenant au moins une sous-étape (50) de distorsion délivrant des combinaisons de filtrage uniquement à partir des5 dits couples de distorsion représentatifs de la ou desdites distorsions, afin de délivrer lesdites combinaisons de filtrage destinées à être appliquées auxdits coefficients de la représentation initiale. 19. Procédé selon l'une quelconque des revendications 16 à 18, caractérisé en ce qu'il comporte une pluralité de sous-étapes (50) de distorsion, 10 récursives les unes par rapport aux autres, de sorte que chaque sous-étape (50) de distorsion reçoit en entrée les combinaisons de filtrage délivrées par la sousétape précédente. 20. Procédé selon l'une quelconque des revendications 16 à 19, caractérisé en ce que ladite au moins une sous-étape (50) de distorsion, permet de 15 déterminer les combinaisons de filtrage correspondant à des modifications à effectuer sur des parties de la représentation initiale qui ne sont pas directement visées par lesdits couples de distorsion. 21. Procédé selon l'une quelconque des revendications 10, 16 ou 18, caractérisé en ce que ledit traitement se décompose en une pluralité de traite20 ments élémentaires et en ce qu'au moins ladite étape (6) de détermination de combinaisons de filtrage est répétée une pluralité de fois, les combinaisons de filtrage résultant des différentes sous-étapes (30) de calcul et/ou (50) de distorsion et correspondant auxdits traitements élémentaires étant arrangées entre elles pour former les combinaisons de filtrage destinées à être appliquées auxdits 25 coefficients de la représentation initiale. 22. Procédé selon l'une quelconque des revendications 1 à 21, caractérisé en ce que ladite étape (6) de détermination des combinaisons de filtrage comprend une sous-étape (32) de conversion desdites combinaisons de filtrage représentées dans un format prédéterminé pour délivrer des paramètres de filtres 30 adaptés pour le paramétrage de filtres donnés en vue de la réalisation de ladite étape (10) d'application des combinaisons de filtrage. 23. Procédé selon l'une quelconque des revendications 1 à 22, caractérisé en ce que ledit traitement est un traitement spatial et en ce que ladite étape (10) d'applications des combinaisons de filtrage comprend des modifications de gain appliquées auxdits coefficients de la représentation initiale. 24. Procédé selon l'une quelconque des revendications 1 à 23, caractérisé en ce que ledit traitement est un traitement spatial et fréquentiel et en ce 5 que ladite étape (10) d'application des combinaisons de filtrage comprend l'application auxdits coefficients de la représentation initiale, de filtres variant en fonction de la fréquence. 25. Procédé selon l'une quelconque des revendications 1 à 24, caractérisé en ce que lesdits ensembles de coefficients représentatifs de la représen10 tation initiale du champ acoustique et de la représentation modifiée sont des ensembles de coefficients dits de Fourier-Bessel. 26. Dispositif de traitement d'une représentation d'un champ acoustique afin de délivrer une représentation de ce champ acoustique par application de ce traitement correspondant à un traitement au moins spatial du champ 15 acoustique, ledit dispositif (60) recevant en entrée une représentation initiale du champ acoustique sous la forme d'un ensemble de coefficients représentatif de ce champ dans le temps et dans les trois dimensions de l'espace et délivrant en sortie une représentation modifiée par les combinaisons de filtrage correspondant audit traitement appliqué auxdits coefficients de la représentation initiale, caracté20 risé en ce qu'il comporte des moyens (61) de définition dudit traitement, comprenant: - des moyens (62, 64, 66) d'établissement d'un ensemble de fonctions de directivité représentatif dudit traitement, sous la forme d'un ensemble de coefficients correspondant à la décomposition de ces fonctions de directivité sur 25 une base d'harmoniques sphériques; et - des moyens (70) de détermination desdites combinaisons de filtrage à partir de cette décomposition pour leur mise en couvre par des moyens (80) de filtrage recevant en entrée lesdits coefficients de la représentation initiale et délivrant en sortie ladite représentation modifiée. 27. Dispositif selon la revendication 26, caractérisé en ce que lesdits moyens (61) de définition dudit traitement comportent un module (62) de saisie directe d'informations de fonctionnement du dispositif dans l'un ou l'autre d'une pluralité de formats. 28. Dispositif selon l'une quelconque des revendications 26 ou 27, caractérisé en ce que lesdits moyens (61) de définition dudit traitement comportent un module (64) d'interprétation automatique d'une consigne de traitement permettant de délivrer automatiquement des informations de fonctionnement dans l'un ou l'autre d'une pluralité de formats prédéterminés. 29. Dispositif selon l'une quelconque des revendications 27 ou 28, caractérisé en ce que lesdites informations de fonctionnement comprennent ledit ensemble de fonctions de directivité représentatif dudit traitement. 30. Dispositif selon la revendication 29, caractérisé en ce que lesdits 10 moyens (61) de définition dudit traitement comportent un module (66) de conversion adapté pour convertir ledit ensemble de fonctions de directivité de l'un ou l'autre d'une pluralité de formats prédéterminés en un ensemble de coefficients correspondant à la décomposition desdites fonctions de directivité sur une base d'harmoniques sphériques. 31. Dispositif selon l'une quelconque des revendications 27 ou 28, caractérisé en ce que lesdites informations de fonctionnement comprennent un ensemble de paramètres dudit traitement. 32. Dispositif selon la revendication 31, caractérisé en ce que lesdits moyens (61) de définition dudit traitement comportent un module de calcul de 20 tout ou partie dudit ensemble de paramètres en fonction de la nature dudit traitement. 33. Dispositif selon l'une quelconque des revendications 26 à 32, caractérisé en ce que lesdits moyens (70) de détermination desdites combinaisons de filtrage comprennent un module (72) de calcul de combinaisons de filtrage par 25 combinaisons linéaires desdits coefficients correspondant à la décomposition sur une base d'harmoniques sphériques dudit ensemble de directivité représentatif dudit traitement. 34. Dispositif selon les revendications 31 et 33 prises ensemble, caractérisé en ce que ledit ensemble des paramètres dudit traitement comprend un 30 ensemble de constantes et en ce que ledit module (72) de calcul par combinaisons linéaires fait intervenir ledit ensemble de constantes en tant que coefficients de pondération. 35. Dispositif selon l'une quelconque des revendications 26 à 33, caractérisé en ce que ledit traitement comprend au moins une distorsion, tout ou partie dudit ensemble de fonctions de directivité formant un ensemble de couples de distorsion correspondant à la ou aux distorsions à appliquer à la représentation initiale et en ce que lesdits moyens (70) de détermination desdites combinaisons de filtrage comprennent en outre un module (74) de distorsion par modifica5 tion des combinaisons de filtrage délivrées par ledit module (72) de calcul, à partir desdits couples de distorsion afin de délivrer lesdites combinaisons de filtrage destinées à être appliquées auxdits coefficients de la représentation initiale. 36. Dispositif selon l'une quelconque des revendications 26 à 32, caractérisé en ce que ledit traitement comprend uniquement une ou plusieurs dis10 torsions, tout ou partie dudit ensemble de fonctions de directivité formant un ensemble de couples de distorsion correspondant à la ou aux distorsions à appliquer à la représentation initiale et en ce que lesdits moyens (70) de détermination desdites combinaisons de filtrage comprennent uniquement un module (74) de distorsion pour délivrer lesdites combinaisons de filtrage uniquement à partir 15 desdits couples de distorsion.