FR 2850183 A1 20040723 FR 0300571 A 20030120 La présente invention concerne un procédé et un dispositif de pilotage d'un ensemble de restitution d'un champ acoustique comportant une pluralité d'éléments de restitution, à partir d'une pluralité de signaux acoustiques ou audiophoniques associés chacun à une direction générale de restitution prédéter5 minée, définie par rapport à un point de l'espace donné. Un tel ensemble de signaux est couramment désigné par l'expression " signal multicanal " et correspond à une pluralité de signaux, appelés canaux, transmis en parallèle ou multiplexés les uns avec les autres, chacun destiné à un élément ou un groupe d'éléments de restitution, disposé dans une direction géné10 raie prédéfinie par rapport à un point donné. Par exemple, un système multicanal classique et connu sous le nom de " 5. 1 ITU-R BF 775-1 " et comporte cinq canaux destinés à des éléments de restitutions placés dans cinq directions générales prédéterminées par rapport à un centre d'écoute, définies par les angles 0 , + 30 , -30 , + 1100 et -1100. Une telle disposition correspond donc à la disposition d'un haut-parleur ou un groupe de haut-parleurs devant au centre, un de chaque côté devant à droite et à gauche et un de chaque côté derrière à droite et à gauche. Les signaux de pilotage étant associés chacun à une direction déterminée, l'application de ces signaux à un ensemble de restitution dont les élé20 ments ne répondent pas à la configuration spatiale prédéterminée, entraîne des déformations importantes du champ acoustique restitué. Il existe des systèmes qui intègrent des moyens de retard sur les canaux, afin de compenser au moins partiellement, l'éloignement des éléments de restitution par rapport au centre d'écoute. Ces systèmes ne permettent cepen25 dant pas de prendre en compte la disposition dans l'espace de l'ensemble de restitution. Il apparaît donc qu'aucun procédé ou système existant, ne permet une restitution de bonne qualité à partir d'un signal de type multicanal avec un ensemble de restitution de configuration spatiale quelconque. La présente invention a pour objet de remédier à ce problème, en définissant un procédé et un système de pilotage de l'ensemble de restitution dont la configuration spatiale est quelconque. L'invention a pour objet un procédé de pilotage d'un ensemble de restitution d'un champ acoustique comportant une pluralité d'éléments de restitution associés chacun à une direction générale de restitution prédéterminée définie par rapport à un point donné, pour obtenir un champ acoustique restitué de caractéristiques spécifiques sensiblement indépendantes des caractéristiques intrinsèques de restitution dudit ensemble, caractérisé en ce qu'il comporte: - une étape de détermination de caractéristiques au moins spatiales dudit ensemble de restitution, permettant la détermination de paramètres représentatifs pour au moins un élément dudit ensemble de restitution de sa position dans les trois dimensions de l'espace par rapport audit point donné - une étape de détermination de filtres d'adaptation à partir desdites 10 caractéristiques au moins spatiales dudit ensemble de restitution et desdites directions générales de restitution prédéterminée associées à ladite pluralité de signaux d'entrée d'informations acoustiques; - une étape de détermination d'au moins un signal de pilotage desdits éléments dudit ensemble de restitution par l'application desdits filtres 15 d'adaptation à ladite pluralité de signaux d'entrée d'informations acoustiques; et - une étape de délivrance dudit au moins un signal de pilotage en vue d'une application auxdits éléments de restitution. Suivant d'autres caractéristiques: - ladite étape de détermination de caractéristiques au moins spatiales 20 dudit ensemble de restitution, comporte une sous-étape de saisie permettant de déterminer tout ou partie des caractéristiques dudit ensemble de restitution; - ladite étape de détermination de caractéristiques au moins spatiales dudit ensemble de restitution, comporte une étape de calibrage permettant de délivrer tout ou partie des caractéristiques dudit ensemble de restitution - ladite sous-étape de calibrage comporte pour au moins un des éléments de restitution: - une sous-étape d'émission d'un signal spécifique vers ledit au moins un élément dudit ensemble de restitution; - une sous-étape d'acquisition de l'onde sonore émise en ré30 ponse par ledit au moins un élément; - une sous-étape de transformation desdits signaux acquis en un nombre fini desdits coefficients représentatifs de l'onde sonore émise; et - une sous-étape de détermination de paramètres spatiaux et/ou acoustiques dudit élément à partir desdits coefficients représentatifs de l'onde sonore émise; - ladite sous-étape de calibrage comporte en outre une sous-étape de 5 détermination de la position dans au moins l'une des trois dimensions de l'espace dudit au moins un élément dudit ensemble de restitution; - ladite étape de calibrage comporte une sous-étape de détermination de la réponse en fréquence dudit au moins un élément dudit ensemble de restitution; - ladite étape de détermination de filtres d'adaptation comprend: - une sous-étape de détermination d'une matrice de décodage représentative de filtres permettant la compensation des altérations de restitution dues aux caractéristiques spatiales dudit ensemble de restitution; - une sous-étape de détermination d'une matrice de rayonnement 15 multicanal idéale représentative des directions générales prédéterminées associes à chaque signal d'informations de la pluralité des signaux d'entrée; et - une sous-étape de détermination d'une matrice représentative desdits filtres d'adaptation à partir de ladite matrice de décodage et de ladite matrice de rayonnement multicanal; ladite étape de détermination de filtres d'adaptation comporte une pluralité de sous-étapes de calcul permettant de délivrer un ordre limite de précision spatiale des filtres d'adaptation, une matrice correspondant à une fenêtre spatiale représentative de la répartition dans l'espace de la précision souhaitée lors de la reconstruction du champ acoustique et une matrice représentative du 25 rayonnement de l'ensemble de restitution, ladite sous-étape de calcul de la matrice de décodage étant réalisée à partir des résultats de ces sous-étapes de calcul; - les matrices de décodage, de rayonnement multicanal idéal et d'adaptation sont indépendantes de la fréquence, l'étape de détermination d'au 30 moins un signal de pilotage desdits éléments dudit ensemble de restitution par l'application desdits filtres d'adaptation correspondant à de simples combinaisons linéaires suivies de retard. - ladite étape de détermination de caractéristiques dudit ensemble de restitution permet la détermination de caractéristiques acoustiques dudit ensem- ble de restitution et ledit procédé comporte une étape de détermination de filtres de compensation de ces caractéristiques acoustiques, ladite étape de détermination d'au moins un signal de pilotage comprenant alors une sous-étape d'application desdits filtres de compensation acoustique; - ladite étape de détermination de caractéristiques acoustiques est adaptée pour délivrer des paramètres représentatifs pour au moins un élément de sa réponse en fréquence; ladite étape de détermination d'au moins un signal de pilotage comporte une sous-étape d'ajustement de gain et d'application de retards afin 10 d'aligner temporellement le front d'onde des éléments de restitution en fonction de leur'distance par rapport audit point donné. L'invention a également pour objet un programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé lorsque ledit programme est exécuté par un ordinateur. L'invention a encore également pour objet un support amovible du type comprenant au moins un processeur de traitement et un élément de mémoire non volatile, caractérisé en ce que ladite mémoire comprend un programme comprenant des instructions de code pour l'exécution des étapes du procédé, lorsque ledit processeur exécute ledit programme. L'invention a encore pour objet un dispositif de pilotage d'un ensemble de restitution d'un champ acoustique comportant une pluralité d'éléments de restitution, comportant des moyens d'entrée d'une pluralité de signaux d'entrée d'informations acoustiques associés chacun à une direction générale de restitution prédéterminée définie par rapport à un point donné, caractérisé en ce qu'il 25 comporte en outre: - des moyens de détermination de caractéristiques au moins spatiales dudit ensemble de restitution, permettant la détermination de paramètres représentatifs pour au moins un élément dudit ensemble de restitution de sa position dans les trois dimensions de l'espace par rapport audit point donné; - des moyens de détermination de filtres d'adaptation à partir desdites caractéristiques au moins spatiales dudit ensemble de restitution et des directions générales de restitution prédéterminée associées à ladite pluralité des signaux d'entrée d'informations acoustiques; et - des moyens de détermination d'au moins un signal de pilotage desdits éléments dudit ensemble de restitution par l'application desdits filtres d'adaptation à ladite pluralité de signaux d'entrée d'informations acoustiques. Suivant d'autres caractéristiques de ce dispositif: - lesdits moyens de détermination des caractéristiques au moins spatiales dudit ensemble de restitution comportent des moyens de saisie directe desdites caractéristiques; - il est adapté pour être associé à des moyens de calibrage permettant la détermination des caractéristiques au moins spatiales dudit ensemble de 10 restitution; - lesdits moyens de calibrage comprennent des moyens d'acquisition d'une onde sonore comportant quatre capteurs de pression disposés selon une forme générale de tétraèdre; lesdits moyens de détermination de caractéristiques sont adaptés 15 pour la détermination de caractéristiques acoustiques d'au moins un desdits éléments de restitution dudit ensemble de restitution, ledit dispositif comportant des moyens de détermination de filtres de compensation acoustique à partir desdites caractéristiques acoustiques et lesdits moyens de détermination d'au moins un signal de pilotage étant adaptés pour l'application desdits filtres de compensa20 tion acoustique; lesdits moyens de détermination des caractéristiques acoustiques sont adaptés pour la détermination de la réponse en fréquence desdits éléments de l'ensemble de restitution. L'invention a également pour objet un appareil de traitement de don25 nées audio et vidéo comportant des moyens de détermination d'une pluralité de signaux d'entrée d'informations acoustiques associés chacun à une direction générale de restitution prédéterminée définie par un point donné, caractérisé en ce qu'il comporte en outre un dispositif de pilotage d'un ensemble de restitution - lesdits moyens de détermination d'une pluralité de signaux d'entrée 30 sont formés d'une unité de lecture et de décodage des disques audio et/ou vidéo numériques. L'invention sera mieux comprise à la lecture de la description qui va 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 schéma d'un système de restitution selon l'invention; - la Fig.3 est un organigramme du procédé de l'invention; - la Fig.4 est un schéma de moyens de calibrage mis en oeuvre dans le procédé de l'invention; - la Fig.5 est un organigramme détaillé de l'étape de calibrage; - la Fig.6 est une représentation simplifiée d'un capteur utilisé pour la mise en oeuvre de l'étape de calibrage; - la Fig.7 est un organigramme détaillé de l'étape de détermination de 10 filtres d'adaptation; et - les Figs.8 et 9 sont des schémas de moyens de détermination de signaux de pilotage; et - la Fig.10 est un schéma d'un mode de réalisation d'un dispositif mettant en oeuvre le 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 O et comportant 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,0), o r désigne la distance par rapport à l'origine O, O l'orientation dans le plan vertical et 0 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éep(r,O,0,t), dont la trans25 formée de Fourier temporelle est notée P(r, O, 0,J) ofdésigne la fréquence. L'invention est fondée sur l'utilisation d'une famille de fonctions spatiotemporelles permettant de décrire les caractéristiques de tout champ acoustique. Dans le mode de réalisation décrit, ces fonctions sont des fonctions dites de Fourier-Bessel sphériques de première espèce appelées par la suite fonc30 tions de Fourier-Bessel. Dans une zone vide de sources sonores et vide d'obstacles, les fonctions de Fourier-Bessel sont solutions de l'équation des ondes et constituent une base qui engendre tous les champs acoustiques produits par des sources sonores situées à l'extérieur de cette zone. Tout champ acoustique tridimensionnel s'exprime donc par une combinaison linéaire des fonctions de Fourier-Bessel, selon l'expression de la transformée de Fourier-Bessel inverse qui s'exprime: P(r,0,0,bf)=4 ? ' P, rm(f)jlj](kr)yp( 0,0) 1=0 m=-I Dans cette équation, les termes Pl,m(J) sont, par définition, les coefficients de Fourier-Bessel du champp(r,,0,, t), k= 2, c est la célérité du son dans l'air (340 ms-1), j/(kr)est la fonction de Bessel sphérique de première espèce d'ordre l définie par j, (x)= x J,+112(x) o J,(x) est la fonction de Bessel de première espèce d'ordre v, et ynm(0,ó) est l'harmonique sphérique réelle d'ordre I et 10 de terme m, avec m allant de -I à 1, définie par: 1 PI(CosO) cos(me) pour m>0 Ynlm(O,)= 1 P10(cosO) pour m=0 -1 Plll(cos) sin(me5) pour m En variante, le procédé de l'invention utilise des bases de fonctions s'exprimant comme des combinaisons linéaires, éventuellement infinies, de fonctions de Fourier-Bessel. Sur la figure 2, on a représenté schématiquement un système de restitution dans lequel on met en oeuvre le procédé de l'invention. Ce système comporte un décodeur ou adaptateur 1 pilotant un ensemble de restitution 2 qui comporte une pluralité d'éléments 31 à 3N, tels que des haut-parleurs, des enceintes acoustiques ou toute autre source sonore ou groupe de sources sonores, agencés de manière quelconque dans un lieu d'écoute 4. On place arbitrairement, dans le lieu d'écoute 4, l'origine O du repère que l'on appelle centre 5 de l'ensemble de restitution. L'ensemble des caractéristiques spatiales, acoustiques et électrodynamiques est considéré comme les caractéristiques intrinsèques de l'ensemble de restitution 2. L'adaptateur 1 reçoit en entrée un signal Si de type multicanal comportant des informations acoustiques à restituer et un signal de définition SL compor10 tant des informations représentatives de caractéristiques au moins spatiales de l'ensemble de restitution 2 et notamment permettant la détermination de paramètres représentatifs pour au moins un élément 3, de l'ensemble de restitution 2 de sa position dans les trois dimensions de l'espace par rapport au point donné 5. A l'issue du traitement correspondant au procédé de l'invention, 15 l'adaptateur 1 émet à l'attention de chacun des éléments ou groupes d'éléments 31 à 3N de l'ensemble de restitution 2, un signal sc1 à sCNde pilotage spécifique. Sur la figure 3, on a représenté schématiquement les étapes principales du procédé selon l'invention mis en òuvre avec un système de restitution tel que celui décrit en référence à la figure 2. Ce procédé comporte une étape 10 de détermination de paramètres de fonctionnement, adaptée pour permettre au moins la détermination des caractéristiques spatiales de l'ensemble de restitution 2. L'étape 10 comporte une étape 20 de saisie des paramètres et/ou une étape 30 de calibrage permettant de déterminer et/ou mesurer des caractéristi25 ques de l'ensemble de restitution 2. Dans le mode de réalisation décrit, l'étape 10 comporte également une étape 40 de détermination de paramètres de description des directions générales prédéterminées associées aux différents canaux du signal d'entrée multicanal SI. A l'issue de l'étape 10, des informations relatives au moins aux diffé30 rentes directions générales prédéterminées associées à chacun des canaux d'entrée ainsi qu'à la position dans les trois dimensions de l'espace de chacun des éléments ou groupes d'éléments 3n de l'ensemble de restitution 2, sont déterminées. Ces informations sont utilisées lors d'une étape 50 de détermination des filtres d'adaptation permettant de prendre en compte les caractéristiques spatiales de l'ensemble de restitution 2 afin de définir des filtres d'adaptation du signal d'entrée multicanal à la configuration spatiale spécifique de l'ensemble de restitution 2. Avantageusement, l'étape 10 permet également de déterminer des caractéristiques acoustiques pour tout ou partie des éléments 31 à 3N de l'ensemble de restitution 2. Dans ce cas, le procédé comporte une étape 60 de détermination de 10 filtres de compensation acoustique permettant de compenser l'influence des caractéristiques acoustiques spécifiques des éléments 31 à 3N. Les filtres définis lors des étapes 50 et avantageusement 60, peuvent ainsi être mémorisés, de sorte que les étapes 10, 50 et 60 ne doivent être répétées qu'en cas de modification de la configuration spatiale de l'ensemble de resti15 tution 2 et/ou de la nature du signal d'entrée multicanal. Le procédé comporte ensuite une étape 70 de détermination des signaux de pilotage sc1 à scN destinés aux éléments de l'ensemble de restitution 2, comportant une sous-étape 80 d'application des filtres d'adaptation déterminés lors de l'étape 50 aux différents canaux ci(t) à cQ(t) formant le signal multicanal 20 d'entrée Si et avantageusement, une sous-étape 90 d'application des filtres de compensation acoustique déterminés lors de l'étape 60. Les signaux sc1 à SCN ainsi délivrés, sont appliqués aux éléments 31 à 3N de l'ensemble de restitution 2, afin de restituer le champ acoustique représenté par le signal multicanal d'entrée Si avec une adaptation optimum aux caracté25 ristiques spatiales et avantageusement acoustiques, de l'ensemble de restitution 2. Il apparaît donc que grâce à la mise en oeuvre du procédé de l'invention, les caractéristiques du champ acoustique restitué sont sensiblement indépendantes des caractéristiques intrinsèques de restitution de l'ensemble de 30 restitution 2 et notamment de sa configuration spatiale. On va maintenant décrire plus en détail les étapes principales du procédé de l'invention. Lors de l'étape 20 de saisie de paramètres un opérateur ou un système de mémoire adapté peut spécifier tout ou partie des paramètres de calcul et notamment: - des paramètres x,, exprimés dans le repère sphérique au moyen des 5 coordonnées ri, On, et 0, et représentatifs de la position des éléments 3n par rapport au centre d'écoute 5; et/ou - des paramètres Hf/), représentatifs de la réponse en fréquence des éléments 3n Cette étape 20 est mise en oeuvre au moyen d'une interface de type 10 classique telle qu'un micro-ordinateur ou tout autre moyen approprié. On va maintenant décrire plus en détail l'étape 30 de calibrage ainsi que des moyens de mise en oeuvre de cette étape. Sur la figure 4 on a représenté le détail de moyens de calibrage. Ils comportent un module 91 de décomposition, un module 92 de détermination de 15 réponse impuisionnelle et un module 93 de détermination de paramètres de calibrage. Les moyens de calibrage sont adaptés pour être connectés à un dispositif d'acquisition sonore 100 tel qu'un microphone ou tout autre dispositif adapté, et pour être connectés tour à tour à chaque élément 3, de l'ensemble de resti20 tution 2 afin de prélever des informations sur cet élément. Sur la figure 5, on a représenté le détail d'un mode de réalisation de l'étape 30 de calibrage mise en oeuvre par les moyens de calibrage décrits précédemment et permettant de mesurer des caractéristiques de l'ensemble de restitution 2. Lors d'une sous-étape 32, les moyens de calibrage émettent un signal spécifique un(t) tel qu'une séquence pseudo-aléatoire MLS (Maximum Length Sequence) à l'attention d'un élément 3n. Le dispositif d'acquisition 100 reçoit, lors d'une sous-étape 34, l'onde sonore émise par l'élément 3n en réponse à la réception du signal un(t) et transmet I signaux cpi(t) à cpjXt) représentatifs de l'onde re30 çue au module de décomposition 91. Lors d'une sous-étape 36, le module de décomposition 91 décompose les signaux captés par le dispositif d'acquisition 100 en un nombre fini de coefficients de Fourier-Bessel qlm(t) Par exemple, le dispositif d'acquisition 100 est constitué de 4 capteurs de pression localisés aux 4 sommets d'un tétraèdre de rayon R ainsi que cela est représenté en référence à la figure 6. Les signaux des 4 capteurs de pression sont alors notés cpl(t) à cp4(t). Les coefficients qo0,o(t) à q1,(t) représentatifs du 5 champ acoustique capté se déduisent des signaux cpl(t) à cp4(t) selon les relations suivantes: Qoo(f) = 1 CP(f) +CF2f(f)+CP3(f) + CP4(f) 3 c Qi=-i(f) =- 27rjRf (CPI(f) -CP2(f) + CP3(f) -CP4(f)) o(8) / 2=jRf (cP1(f) + CP2(f) - CP3(f) - CP4(f)) Qo(f = V, 2rjRf 3 c Qii(f) = 8Vr 2 jRf (CPl (f) - CP2(f) - CP3 (f) + CP4(f) Dans ces relations CP(/i à CP4(/) sont les transformées de Fourier de 10 cpl(t) à cp4(t) et Q0,00o( à Ql,i(0) sont les transformées de Fourier de q0,0(t) à ql,l(t). Lorsque ces coefficients sont définis par le module 91, ils sont adressés au module 92 de détermination de réponse. Lors d'une sous-étape 38, le module 92 de détermination de réponse détermine les réponses impulsionnelles hpl,m(t) qui relient les coefficients de Fou15 rier-Bessel ql,,m(t) et le signal émis u,(t). La méthode de détermination dépend du signal spécifique émis. Le mode de réalisation décrit utilise une méthode adaptée aux signaux de type MLS, comme par exemple la méthode de corrélation. La réponse impulsionnelle délivrée par le module 92 de détermination de réponse est adressée au module 93 de détermination de paramètres. Lors d'une sous-étape 39, le module 93 déduit des informations sur des éléments de l'ensemble de restitution. Dans le mode de réalisation décrit, le module 93 de détermination de paramètres détermine la distance r, entre l'élément 3n et le centre 5 à partir de sa réponse hpo,o(t) et de la mesure du temps mis par le son pour se propager de 25 l'élément 3n au dispositif d'acquisition 100, grâce à des méthodes d'estimation de retard sur la réponse hpo,o(t). La direction (0,ón,) de l'élément 3n est déduite par calcul du maximum de la transformée de Fourier sphérique inverse appliquée aux réponses hpo, o(t) à hpl,l(t) prises à l'instant t o hpo,o(t) présente un maximum. Avantageusement, les 30 coordonnées On et ó,f sont estimées sur plusieurs instants, choisis de préférence autour de l'instant o hpo,o(t) présente un maximum. La détermination finale des coordonnées O,, et 0, est obtenue au moyen de techniques de moyennage entre les différentes estimations. Ainsi, dans le mode de réalisation décrit, le dispositif d'acquisition 100 5 est en mesure d'encoder sans ambiguté l'orientation d'une source dans l'espace. En variante, les coordonnées A, et, sont estimées à partir d'autres réponses parmi les hplm(t) disponibles ou sont estimées dans le domaine fréquentiel à partir des réponses HPIm(I), correspondant aux transformées de Fourier des réponses hplm(t). Ainsi l'étape 30 permet de déterminer les paramètres r,, O, et p Dans le mode de réalisation décrit, le module 93 délivre également la fonction de transfert Hf(/) de chaque élément 3n, à partir des réponses hplm(t) issues du module 92 de détermination de réponse. Une première solution consiste à construire la réponse hp'oo(t) corres15 pondant à la sélection de la partie de la réponse hpo,o(t) qui comporte un signal non nul et dénué des réflexions introduites par le lieu d'écoute 4. La réponse en fréquence HW(/) est déduite par transformée de Fourier de la réponse hp'o,o(t) préalablement fenêtrée. La fenêtre peut être choisie parmi les fenêtres classiques de lissage, comme par exemple rectangulaire, Hamming, Hanning, et Blackman. Une seconde solution plus complexe consiste à appliquer un lissage sur le module et avantageusement sur la phase de la réponse en fréquence HPoo(f) obtenue par transformée de Fourier de la réponse hpoo(t). Pour chaque fréquences, le lissage est obtenu par convolution de la réponse HPoo(J) par une fenêtre centrée surf Cette convolution correspond à un moyennage de la ré25 ponse HPoo(/) autour de la fréquences La fenêtre peut être choisie parmi les fenêtres classiques, comme par exemple rectangulaires, triangles et Hamming. Avantageusement, la largeur de la fenêtre varie avec la fréquence. Par exemple, la largeur de la fenêtre peut être proportionnelle à la fréquence f à laquelle est appliquée le lissage. Par rapport à une fenêtre fixe, une fenêtre variable avec la 30 fréquence permet d'éliminer au moins partiellement l'effet de salle dans les hautes fréquences tout en évitant un effet de troncature de la réponse HPoo(j) dans les basses fréquences. Les sous-étapes 32 à 39 sont répétées pour tous les éléments 31 à 3N de l'ensemble de restitution 2. En variante, les moyens de calibrage comportent d'autres moyens d'acquisition d'informations relatives aux éléments 31 à 3N, tels que des moyens de mesure de position au laser, des moyens de traitement du signal mettant en oeuvre des techniques de formation de voies ou tout autre moyen approprié. Les moyens mettant en oeuvre l'étape 30 de calibrage sont constitués par exemple d'une carte électronique ou d'un programme d'ordinateur ou de tout autre moyen approprié. L'étape 40 permet ainsi que cela a été dit précédemment, de déterminer des paramètres décrivant le format du signal multicanal d'entrée et notam10 ment les directions prédéterminées générales associées à chaque canal. Cette étape 40 peut correspondre à une sélection par un opérateur d'un format parmi une liste de formats associés chacun à des paramètres mémorisés, et peut également correspondre à une détection automatique de format effectuée sur le signal multicanal d'entrée. Alternativement, le procédé est adapté 15 pour un unique format de signal multicanal donné. Dans encore un autre mode de réalisation, l'étape 40 permet à unutilisateur de spécifier son propre format en saisissant manuellement les paramètres décrivant les directions associées à chaque canal. Il apparaît que les étapes 20, 30 et 40 formant l'étape 10 de détermi20 nation de paramètres, permettent au moins la détermination de paramètres de positionnement dans l'espace des éléments 3n de l'ensemble de restitution 2 et du format du signal multicanal Si. Sur la figure 7, on a représenté un organigramme détaillé de l'étape 50 de détermination des filtres d'adaptation. Cette étape comporte une pluralité de sous-étapes de calcul et de détermination de matrices représentatives des paramètres déterminés préalablement. Ainsi, lors d'une sous-étape 51, un paramètre L, appelé ordre limite représentatif de la précision spatiale souhaitée lors de l'étape 50 de détermina30 tion des filtres d'adaptation, est déterminé par exemple de la manière suivante: - le plus petit angle amin formé par une paire d'éléments de l'ensemble de restitution 2 est calculé automatiquement au moyen d'une relation trigonométrique, telle que par exemple: an *,n2' = aCos(sinO,,, sinOn2cos(n 1 -0n2) + CoS n1 CoSO,,2) amin = min(anl,,l2) parmi l'ensemble des couples (ni, n2) tels que nl n2; et - ensuite, l'ordre maximum L est déterminé automatiquement comme étant le plus grand entier respectant la relation suivant: L L'étape 50 de détermination de filtres d'adaptation comprend ensuite une sous-étape 52 de détermination d'une matrice W de pondération du champ acoustique. Cette matrice W correspond à une fenêtre spatiale W(rf) représentative de la répartition dans l'espace de la précision souhaitée lors de la reconstruc10 tion du champ. Une telle fenêtre permet de spécifier la taille et la forme de la zone o le champ doit être correctement reconstruit. Par exemple, il peut s'agir d'une boule centrée sur le centre 5 de l'ensemble de restitution. Dans le mode de réalisation décrit, la fenêtre spatiale et la matrice W sont indépendantes de la fréquence. West une matrice diagonale de taille (L+1)2 contenant des coefficients de pondération W1 et dans laquelle chaque coefficient W1 se trouve 21+1 fois à la suite sur la diagonale. La matrice W a donc la forme suivante: 0. o w W=. DW OWL Dans le mode de réalisation décrit, les valeurs prises par les coeffi20 cients W1 sont les valeurs d'une fonction telle qu'une fenêtre de Hamming de taille de 2L+1 évaluée en 1, de sorte que le paramètre W1 est déterminé pour I allant de 0 àL. L'étape 50 comporte ensuite une sous-étape 53 de détermination d'une matrice M représentative du rayonnement de l'ensemble de restitution no25 tamment à partir des paramètres de position x,. La matrice de rayonnement M permet de déduire des coefficients de Fourier-Bessel représentant le champ acoustique qu'émet chaque élément 3, de l'ensemble de restitution en fonction du signal qu'il reçoit. M est une matrice de taille (L+1)2 sur N, constituée d'éléments MImn,, les indices l,m désignant la ligne 12+l+m et n désignant la colonne n. La matrice M a donc la forme suivante: M 0,0,1 M 0,0,2... M 0,0,N M1,-1,1 M1, -1,2...M1,-1,N M],o,1 MI,0,2...... M,O,N M l,l,l M [,1,2...... M l,l,N ML L,-L,lM L M L,-L,N M L,o,1 ML,o,2.... ML,0o,N M L,L,1 M L,L,2......M L,L, N Dans le mode de réalisation décrit, les éléments MImn,,,, s'obtiennent à partir d'un modèle de rayonnement en ondes planes, de sorte que: La matrice M ainsi définie est représentative du rayonnement de 10 I'ensemble de restitution. En particulier, M est représentative de la configuration spatiale de l'ensemble de restitution. Les sous-étapes 51 à 53 peuvent être exécutées séquentiellement ou simultanément. L'étape 50 de détermination de filtres d'adaptation comprend ensuite 15 une sous-étape 54 de prise en compte de l'ensemble des paramètres du système de restitution 2 déterminés précédemment, afin de délivrer une matrice D de décodage représentative de filtres dits de reconstruction. En effet, les éléments Dn, l,m(f) de la matrice D correspondent à des filtres de reconstruction qui, appliqués aux coefficients de FourrierBessel Pl,m(/) 20 d'un champ acoustique connu, permettent de déterminer les signaux de pilotage d'un ensemble de restitution pour reproduire ce champ acoustique. La matrice de décodage D est donc l'inverse de la matrice de rayonnement M. La matrice D est obtenue à partir de la matrice M au moyen de métho25 des d'inversion sous contraintes faisant intervenir des paramètres d'optimisation supplémentaires. Dans le mode de réalisation décrit, l'étape 50 est adaptée pour réaliser une optimisation grâce à la matrice de pondération du champ acoustique W qui permet notamment de réduire la distorsion spatiale dans le champ acoustique reproduit. Cette matrice D est délivrée notamment à partir de la matrice M, selon l'expression suivante: D=(MTWM)-I MTW dans laquelle MTest la matrice transposée conjuguée de M. Dans le mode de réalisation décrit, les matrices M et W sont indépendantes de la fréquence, de sorte que la matrice D est également indépendante de la fréquence. Elle est constituée d'éléments notés D,,,,, organisés de la ma10 nière suivante: D1,o,o Dl,1, 1 D,, D1l,,o D.,., D1,L,-L. D1,L,O'' D,L,L D2,0,o D2,-I D2 l,,o D2,1, ... D2,,L.. "D2.,L,O - D2,L,L DN,O,ODN, 1,-IDN,I,ODN,1,r DN,L,-L *DN,L,O "DN, L,L L'étape 54 permet ainsi de délivrer la matrice D représentative de filtres dits de reconstruction et permettant la reconstruction d'un champ acoustique à partir d'une configuration quelconque de l'ensemble de restitution. Grâce à 15 cette matrice, le procédé de l'invention permet de prendre en compte la configuration de l'ensemble de restitution 2 et notamment de compenser les altérations du champ acoustique dues à sa configuration spatiale spécifique. En variante, les paramètres relatifs à l'ensemble de restitution 2 peuvent être variables en fonction de la fréquence. Par exemple, dans un tel mode de réalisation, chaque élément D,,,,,(,) de la matrice D peut être déterminé en associant à chacun des N signaux de pilotage une fonction de directivité D,(O,óJ) spécifiant à chaque fréquence f l'amplitude, et avantageusement la phase souhaitée sur le signal de pilotage sc. dans le cas d'une onde plane dans la direction (0,ó). On entend par fonction de directivité D,(0,ó> une fonction 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. Dans le mode de réalisation décrit, les fonctions de directivités sont indépendantes de la fréquence et notés D,(O,ó). Ces fonctions de directivité D,(0,5) peuvent être déterminées en spécifiant que certaines grandeurs physiques entre un champ idéal et le même champ reproduit par l'ensemble de restitution respectent des lois prédéterminées. Par exemple, ces grandeurs peuvent être la pression au centre et l'orientation du vecteur vitesse. Dans certains cas, on souhaite que seulement 3 signaux de pilotage soient actifs pour reproduire une onde plane. Les signaux de pilotage actifs, no5 tés scni à SCn3, sont ceux qui alimentent les éléments de restitution dont les directions sont les plus proches de la direction (O,ó) de l'onde plane. Les éléments de restitution actifs, notés 3n1 à 3n3, forment un triangle contenant la direction (O,') de l'onde plane. Dans ce cas, les valeurs des directivités D,1(0,ó) à D,3(,q') associés aux 3 éléments actifs 3n1 à 3n3 sont donnés par: r-z 1] Tf-lr avec sin,,lcosq5r,, sinBn, coso,,, sinO, ,,cosq5.,,0 sin 0cos1 0 X sin 0,l sin, sin 2, sin.,,, sin|, sin 0 sió= 1 COS t4. COS n2 COSP, SnfOf3 COtJn CO = \ 9ccos 0 / Dans cette relation, a correspond au vecteur contenant [D1( 0,ó')...D.3(0,ó)] et les directions (0,1,óql), (0n2, ,2) et (0n3,ó,3) correspondent 15 respectivement aux directions des éléments 3n1, 3n2 et 3n3. On considère que les valeurs des directivités D,(0,ó) correspondant aux éléments de restitution non actifs sont nulles. La relation précédente est répétée pour K directions (Ok,ók) d'ondes planes différentes. Ainsi, chacune des fonctions de directivité D,(0,ó) est fournie 20 sous la forme d'une liste de K échantillons. Chaque échantillon est fournit sous la forme d'un couple {((1Ok,k),D,(Ok,qk))} o (Ok,k) est la direction de l'échantillon k et o Dn(Ok,k) est la valeur de la fonction de directivité associée au signal de pilotage scn pour la direction (Ok, k). Pour chaque fréquence f, les coefficients Dn,,,,,m(f) de chaque fonction 25 de directivité sont déduits des échantillons {((Ok,ók), D(0k,,k))}. Ces coefficients sont obtenus par inversion du processus d'échantillonnage angulaire qui permet de déduire les échantillons de la liste {((O,ók), Dn(Ok,k))} à 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 30 échantillons. Dans d'autres modes de réalisation, les fonctions de directivité sont directement fournies sous la forme de coefficients Dn,im() de type Fourrier- Bessel. Les coefficients Dn,i,m(I) ainsi déterminés sont utilisés pour former la matrice D. L'étape 50 comporte ensuite une étape 55 de détermination d'une matrice de rayonnement multicanale idéale S représentative des directions généra5 les prédéterminées associées à chaque canal du signal multicanal d'entrée SI. La matrice S est représentative du rayonnement d'un ensemble de restitution idéal, c'est à dire respectant parfaitement les directions générales prédéterminées du format multicanal. Chaque élément Sl,m,q(f) de la matrice S permet de déduire les coefficients de Fourier-Bessel Pl, m(/) du champ acoustique idéale10 ment restitué par chaque canal cq(t). La matrice S est déterminée en associant à chaque canal d'entrée cq(t) et avantageusement pour chaque fréquencef, une figure de directivité représentative d'une distribution de sources supposées émettre le signal du canal cq(t). La distribution de sources est donnée sous la forme de coefficients 15 d'harmoniques sphériques Smq(f) Les coefficients Sl,m,q(f) sont rangés dans la matrice S de taille (L+1)2 sur Q, o Q est le nombre de canaux. Dans le mode de réalisation décrit, l'étape de mise en forme associe à chaque canal Cq(t) une source d'onde plane orientée dans la direction (q, bq) correspondant à la direction (Oc,") associée au canal cq(t) dans le format multi20 canal d'entrée. Les coefficients Sl,m,q(f) sont alors indépendants de la fréquence. Ils sont notés Sl,m,q et s'obtiennent par la relation: Sl,m,q=y, (OqIq) Dans d'autres modes de réalisation, la matrice de rayonnement idéale S associe une distribution discrète de sources d'ondes planes à certains canaux 25 pour simuler l'effet d'une ceinture de haut-parleurs. Dans ce cas, les coefficients Sli,m,q sont obtenus par sommation des contributions de chacune des sources élémentaires. Dans encore d'autres modes de réalisation, la matrice de rayonnement idéale S associe certains canaux cq(t) à une distribution continue de sources 30 d'ondes planes décrite par une fonction de directivité Sq(0,ó). Dans ce cas, les coefficients Sl,m,q de la matrice S sont obtenus directement par transformée de Fourier Sphérique de la fonction de directivité Sq(O,ó). Dans ces modes de réalisation, la matrice S est indépendante de la fréquence. Dans d'autres modes de réalisation plus complexes la matrice S associe à certains canaux, une distribution de sources produisant un champ diffus. Dans ce cas, la matrice S varie avec la fréquence. Ces modes de réalisation sont adaptés aux formats multicanaux qui considèrent différemment les canaux avants 5 et arrières. Par exemple, dans les applications destinées à la restitution dans des salles cinéma, les canaux arrières sont souvent destinés à recréer une ambiance diffuse. Dans d'autres modes de réalisation, la matrice S associe à certains canaux des sources sonores dont la réponse n'est pas plate. Par exemple, dans 10 le cas o le format multicanal associe au canal cq(t) une source d'onde plane de réponse en fréquence H(q)(W, les SI,m,q(f) varient avec la fréquence et s'obtiennent par la relation: Si,m,q()=yrl( Oq,q)H(J)W Si le format multicanal associe à certains canaux une superposition 15 des types de distributions de sources pré-citées, les coefficients Simq(P de la matrice de rayonnement s'obtiennent par sommation des coefficients associés à chaque type de distribution de source. Enfin, l'étape 50 comporte une sous-étape 56 de détermination d'une matrice d'adaptation spatiale A correspondant aux filtres d'adaptation à appliquer 20 au signal d'entrée multicanal pour obtenir une restitution optimum prenant en compte la configuration spatiale de l'ensemble de restitution 2. La matrice d'adaptation spatiale A est obtenue à partir des matrices de mise en forme S et de décodage D au moyen de la relation: A =DS La matrice d'adaptation A permet de générer des signaux sai(t) à saN(t) adaptés à la configuration spatiale de l'ensemble de restitution à partir des canaux ci(t) à CQ(t). Chaque élément Anqj) est un filtre spécifiant la contribution du canal cq(t) au signal adapté sa,(t). Grâce à la matrice d'adaptation A, le procédé 30 de l'invention permet la restitution optimum du champ acoustique décrit par le signal multicanal'par un ensemble de restitution de configuration spatiale quelconque. Dans le mode de réalisation décrit, les matrices D et S sont indépendantes de la fréquence et la matrice A également. Dans ce cas, les éléments de 35 la matrice A sont des constantes notées An,q et chacun des signaux adaptés saz(t) à saN(t) s'obtient par de simples combinaisons linéaires des canaux d'entrée cj(t) à CQ(t), le cas échéant suivies de retard ainsi qu'il sera décrit ci- après. Les filtres représentés par la matrice A peuvent être mis en oeuvre sous différentes formes de filtres et lou méthodes de filtrage. Dans le cas o les 5 filtres utilisés sont paramétrés directement avec des réponses en fréquence, les coefficients Anq(j sont directement délivrés par l'étape 50. Avantageusement, l'étape 50 de détermination de filtres d'adaptation comprend une sous-étape 57 de conversion afin de déterminer les paramètres des filtres pour d'autres méthodes de filtrage. Par exemple, les combinaisons de filtrage Anq(f) sont converties en: - des réponses impulsionnelles finies an,q(t) calculées par transformée de Fourier temporelle inverse de An,q(j), chaque réponse impulsionnelle an,q(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 Anvq() avec des méthodes d'adaptation. A l'issue de l'étape 50 les paramètres des filtres d'adaptation Anq L'étape 60 permet ainsi que cela a été dit précédemment, de détermi20 ner les filtres de compensation des caractéristiques acoustiques des éléments de l'ensemble de restitution 2 dans le cas o des paramètres relatifs à ces caractéristiques acoustiques tels que les réponses en fréquence Hn(f), sont déterminés lors de l'étape 10 de détermination des paramètres. La détermination de tels filtres, notés H,'>(J), à partir des réponses en 25 fréquence HW(), peut être réalisée de manière classique en appliquant des méthodes d'inversion de filtres, comme par exemple l'inversion directe, les méthodes de déconvolution, les méthodes Wiener ou d'autres. En fonction des modes de réalisation, la compensation porte uniquement sur l'amplitude de la réponse ou encore sur l'amplitude et la phase. Cette étape 60 permet de déterminer un filtre de compensation pour chaque élément 3n de l'ensemble de restitution 2 en fonction de ses caractéristiques acoustiques spécifiques. De même que précédemment, ces filtres peuvent être mis en oeuvre sous différentes formes de filtres et/ou méthodes de filtrage. Dans le cas o les filtres utilisés sont paramétrés directement avec des réponses en fréquence, les réponses H2(l)(f) sont directement appliquées. Avantageusement, l'étape 60 de détermination de filtres de compensation comprend une sous-étape de conversion afin de déterminer les paramètres des filtres pour d'autres méthodes de fil5 trage. Par exemple, les combinaisons de filtrage Hg(,)(j) sont converties en des réponses impulsionnelles finies h(P)(t) calculées par transformée de Fourier temporelle inverse de Hk(')(f), chaque réponse impulsionnelle he( (t) est échantillonnée puis tronquée à une longueur propre 10 à chaque réponse; ou - des coefficients de filtres récursifs à réponses impulsionnelles infinies calculées à partir des Hnl)(j) avec des méthodes d'adaptation. A l'issue de l'étape 60 les paramètres des filtres de compensation Hn(j)W sont fournis. On va maintenant décrire plus en détail l'étape 70 de détermination de signaux de pilotage. Cette étape 70 comporte une sous-étape 80 d'application des filtres d'adaptation représentés par la matrice A au signal multicanal d'entrée SI correspondant au champ acoustique à restituer. Ainsi que cela a été dit précédemment, 20 les filtres d'adaptation Anq(j) intègrent les paramètres caractéristiques de l'ensemble de restitution 2. Lors de la sous-étape 80, des signaux adaptés sai(t) à saN(t) sont obtenus par l'application des filtres d'adaptation An,q(f) aux canaux ci(t) à CQ(t) du signal Si. Dans le mode de réalisation décrit, la matrice d'adaptation A est indépendante de la fréquence et les coefficients d'adaptation Anq sont appliqués de la façon suivante Q Vn(t) =2Ec4(t) Anq q=I L'adaptation se poursuit par un ajustement des gains et l'application 30 de retards afin d'aligner temporellement les fronts d'onde des éléments 31 à 3N de l'ensemble de restitution 2 par rapport à l'élément le plus éloigné. Les signaux adaptés sai(t) à saN(t) se déduisent des signaux vi(t) à VN(t) selon l'expression: san(t) = rn vn(t - max( rn) n Dans d'autres modes de réalisation, la matrice d'adaptation A varie avec la fréquence et les filtres d'adaptation An,q(f) sont appliqués de la façon suivante: Q Vn (O) =Cq(I) Anq() q=1 avec Cq-) la transformée de Fourier temporelle du canal cq(t) et Vn(I) défini par: Vn(/) = SAn(f) e-27rjrn f /c ln o SAn() est la transformée de Fourier temporelle de san(t). Selon la forme des paramètres des filtres d'adaptation Anq(), chaque filtrage des canaux cq(t) par les filtres d'adaptation An,q(j) peut être réalisé selon des méthodes classiques de filtrage, telles que par exemple: - les paramètres sont directement les réponses en fréquence Anq(j), et le filtrage est effectué dans le domaine fréquentiel, par exemple, à l'aide des 15 techniques usuelles de convolution par blocs; - les paramètres sont directement les réponses impulsionnelles finies an, q(t), et le filtrage est effectué dans le domaine temporel par convolution; ou - les paramètres sont les coefficients de filtres récursifs à réponses impulsionnelles infinies, et le filtrage est effectué dans le domaine temporel au 20 moyen des relations de récurrence. La sous étape 80 se termine par un ajustement des gains et l'application de retards afin d'aligner temporellement les fronts d'onde des éléments 31 à 3N de l'ensemble de restitution 2 par rapport à l'élément le plus éloigné. Les signaux adaptés sai(t) à say(t) se déduisent des signaux v(t) ( à vN(t) se25 Ion l'expression: scn(t) = rn vnt -max( rn)-rn Sur la figure 8, on a représenté la structure de filtrage correspondant à la sous étape 80 d'application des filtres d'adaptation spatiale telle que décrite précédemment. Avantageusement, l'étape 70 comporte une sous-étape 90 de compensation des caractéristiques acoustiques de l'ensemble de restitution. Chaque filtre de compensation H,()(f) est appliqué au signal adapté sa,(t) correspondant afin d'obtenir le signal de pilotage sc,(t) de l'élément 3,, selon la relation SCn(f) = SAn(/) H,(t)() o SC,(f) est la transformée de Fourier temporelle de sc(t) et o SA,(I) est la transformée de Fourier temporelle de sa,(t). L'application des filtres de compensation HYkJ) des caractéristiques acoustiques est décrite en référence à la figure 9. Selon la forme des paramètres de ces filtres, chaque filtrage des signaux sa,(t) peut être réalisé selon des méthodes classiques de filtrage, telles 10 que par exemple: - dans le cas o les paramètres de filtrage sont des réponses en fréquence H,()(f), 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 le cas o les paramètres de filtrage sont des réponses impulsionnelles h(')(t), le filtrage peut être effectué dans le domaine temporel par convolution temporelle; - dans le cas o les paramètres de filtrage sont des coefficients de relations de récurrence, le filtrage peut être réalisé dans le domaine temporel au 20 moyen de filtres récursifs à réponse impulsionnelle infinie. Dans certains modes de réalisation simplifiés, le procédé de l'invention ne compense pas les caractéristiques acoustiques spécifiques des éléments de l'ensemble de restitution. Dans ce cas, l'étape 60 ainsi que la sousétape 90 ne sont pas réalisées et les signaux adaptés sai(t) à sapit) correspondent directe25 ment aux signaux de pilotage sc1 à scN. Par l'application du procédé de l'invention, chaque élément 31 à 3N reçoit donc un signal de pilotage spécifique sc1 à SCN et émet un champ acoustique qui contribue à la reconstruction optimale du champ acoustique à restituer. En effet, le pilotage simultané de l'ensemble des éléments 31 à 3N permet une re30 construction optimale du champ acoustique correspondant au signal multicanal d'entrée par l'ensemble de restitution 2 dont la configuration spatiale est quelconque, soit encore ne correspond pas à une configuration fixe. Par ailleurs, d'autres modes de réalisation du procédé de l'invention peuvent être envisagés et notamment des modes de réalisation inspirés de techniques décrites dans la demande de brevet en France déposée le 28 février 2002, sous le n0 02 02 585. Notamment, l'étape 50 de détermination des filtres d'adaptation spatiale peut prendre en compte de nombreux paramètres d'optimisation tel que: G,(p, représentatif du gabarit de l'élément 3n de l'ensemble de restitution spécifiant la bande de fréquence de fonctionnement de cet élément; - Ni.mn(f, représentatif de la réponse spatio-temporelle de l'élément 10 3n correspondant au champ acoustique produit dans le lieu d'écoute 4 par l'élément 3n, lorsque celui-ci reçoit en entrée un signal impulsionnel; - W(rj), décrivant pour chaque fréquence f considérée une fenêtre spatiale représentative de la répartition dans l'espace de contraintes de reconstruction du champ acoustique, ces contraintes permettant de spécifier la réparti15 tion dans l'espace de l'effort de reconstruction du champ acoustique; - W1çp, décrivant directement sous forme de pondération des coefficients de Fourier-Bessel et pour chaque fréquencef considérée, une fenêtre spatiale représentative de la répartition dans l'espace de contraintes de reconstruction du champ acoustique; - R(f), représentatif, pour chaque fréquences considérée, du rayon de la fenêtre spatiale lorsque celle-ci est une boule; - u(), représentatif, pour chaque fréquences considérée, de la capacité d'adaptation locale souhaitée à l'irrégularité spatiale de la configuration de l'ensemble de restitution; - {(lk, m*)}(f), constituant pour chaque fréquences considérée, une liste de fonctions spatio-temporelles dont la reconstruction est imposée; - L(f), imposant, pour chaque fréquence f considérée, l'ordre limite de détermination de filtres; - RM(t), définissant, pour chaque fréquences considérée, le modèle de 30 rayonnement des éléments 31 à 3N de l'ensemble de restitution 2. Tout ou partie de ces paramètres d'optimisation peuvent intervenir lors de la sous-étape 54 de détermination de la matrice de décodage D. Ainsi, comme cela est décrit dans la demande de brevet en France déposée sous le numéro 02 02 585, les paramètres Nl,,,ff) et RM(W) interviennent dans la sous-étape 53 de détermination de la matrice de rayonnement M, les paramètres W(rj), W1(J), R(j) interviennent dans la sous-étape 52 de détermination de la matrice W, les paramètres {(lk, mk)}(J) interviennent dans une sous-étape supplémentaire dans la détermination d'une matrice F. La matrice de décodage D est alors déterminée 5 lors de la sous-étape 54, pour chaque fréquence f en fonction des matrices M, W et F et des paramètres G,(j) et u(f. Toujours, selon la demande de brevet 02 02 585, le calcul de la matrice D peut être effectué fréquence par fréquence en considérant uniquement les élements actifs pour chaque fréquence considérée. Cette méthode de détermina10 tion de la matrice D fait intervenir le paramètre GW(j) et permet d'exploiter au mieux un ensemble de restitution dont les éléments ont des bandes de fréquences de fonctionnement différentes. Il apparaît que la mise en oeuvre du procédé de l'invention décrit ici est plus efficace et donc plus rapide que les procédés existants et notamment 15 que le procédé décrit dans la demande de brevet français déposée sous le numéro 02 02 585. En effet, pour adapter un signal multicanal comportant Q canaux à un ensemble de restitution comportant N éléments avec une précision spatiale d'ordre L, il apparaît que le procédé de l'invention nécessite QxN filtres 20 d'adaptation au lieu des Q(L+1)2 + (L+ 1)2N filtres nécessaires à la mise en oeuvre du procédé décrit dans la demande de brevet en France déposée sous le n0 02 02 585. Par exemple, l'adaptation d'un signal " 5.1 ITU-R BF 775-1 " à un ensemble de restitution à 5 haut parleurs avec une précision d'ordre 5 nécessite 25 25 filtres au lieu de 360 filtres. Sur la figure 10, on a représenté un schéma d'un mode de réalisation d'un appareil mettant en oeuvre le procédé tel que décrit précédemment. Cet appareil comprend l'adaptateur 1 qui est formé d'une unité 110 délivrant un signal multicanal telle qu'une unité de lecture de disques audio-vidéo 30 dite lecteur DVD 112. Le signal multicanal délivré par l'unité 110 est destiné aux éléments de l'ensemble de restitution 2. Le format de ce signal Si est reconnu automatiquement par l'adaptateur 1 qui est adapté pour lui faire correspondre des paramètres décrivant la direction générale prédéterminée associée à chaque canal du signal SI. Selon l'invention, cet adaptateur 1 intègre également une unité de calcul supplémentaire 114 ainsi que des moyens de saisie d'informations 116. Par exemple, les moyens de saisie 116 sont formés d'une interface infrarouge avec une télécommande ou encore avec un ordinateur et permettent à 5 un utilisateur de déterminer les paramètres définissant les positions dans l'espace des éléments de restitution 31 à 3N. Ces différents paramètres sont utilisés par le calculateur 114 pour déterminer la matrice A définissant les filtres d'adaptation. Ultérieurement, le calculateur 114 applique ces filtres d'adaptation au 10 signal multicanal SI afin de délivrer les signaux de pilotage sc1 à SCN à destination de l'ensemble de restitution 2. Bien entendu, le dispositif mettant en oeuvre l'invention peut prendre d'autres formes, telles qu'un logiciel mis en oeuvre sur un ordinateur ouencore un dispositif complet intégrant des moyens de calibrage ainsi que des moyens de 15 saisie et de détermination des caractéristiques de l'ensemble de restitution plus complet. Ainsi, le procédé peut également être mis en oeuvre sous la forme d'un dispositif dédié a l'optimisation de systèmes de restitution multicanal, extérieur à un décodeur audio-vidéo et associé à celui-ci. Dans ce cas, le dispositif 20 est adapté pour recevoir en entrée un signal multicanal et délivrer en sortie des signaux de pilotage d'éléments d'un ensemble de restitution. Avantageusement, le dispositif est adapté pour être connecté au dispositif 100 d'acquisition nécessaire à l'étape de calibrage et/ou est muni d'une interface permettant de saisir des paramètres, notamment, la position des élé25 ments de l'ensemble de restitution et éventuellement le format multicanal d'entrée. Un tel dispositif d'acquisition 100 peut être relié de manière filaire ou non filaire (radio, infra-rouge) et peut être intégré à un accessoire, comme une télécommande, ou être indépendant. Le procédé peut être mis en oeuvre par un dispositif intégré dans un élément d'une chaîne audio-vidéo chargé du traitement de signaux multicanal, comme par exemple un processeur ou décodeur dit " surround ", un amplificateur audio-vidéo intégrant des fonctions de décodage multicanal ou encore une chaîne audio-vidéo complètement intégrée. Le procédé de l'invention peut également être mis en oeuvre dans une carte électronique ou dans une puce dédiée. Avantageusement, il peut être intégré sous la forme d'un programme dans un processeur de traitement de signal (DSP). Le procédé peut prendre la forme d'un programme informatique destiné à être exécuté par un ordinateur. Le programme reçoit en entrée un signal multicanal et délivre les signaux de pilotage d'un ensemble de restitution éventuellement intégré à cet ordinateur. Par ailleurs, les moyens de calibrage peuvent être réalisés en mettant 10 en oeuvre un procédé différent de celui décrit précédemment, tel que par exemple, un procédé inspiré de techniques décrites dans la demande de brevet en France déposée le 7 mai 2002 sous le numéro 02 05 741. REVENDICATIONS 1. Procédé de pilotage d'un ensemble de restitution (2) d'un champ acoustique comportant une pluralité d'éléments de restitution (3.) à partir d'une pluralité de signaux d'entrée d'informations acoustiques (SI) associés chacun à 5 une direction générale de restitution prédéterminée définie par rapport à un point donné (5) de l'espace, pour obtenir un champ acoustique restitué de caractéristiques spécifiques sensiblement indépendantes des caractéristiques intrinsèques de restitution dudit ensemble (2), caractérisé en ce qu'il comporte: - une étape (10) de détermination de caractéristiques au moins spatia10 les dudit ensemble de restitution (2), permettant la détermination de paramètres représentatifs pour au moins un élément (3n) dudit ensemble de restitution (2) de sa position dans les trois dimensions de l'espace par rapport audit point donné (5); - une étape (50) de détermination de filtres d'adaptation (A) à partir 15 desdites caractéristiques au moins spatiales dudit ensemble de restitution (2) et desdites directions générales de restitution prédéterminée associées à ladite pluralité de signaux d'entrée d'informations acoustiques (SI); - une étape (70) de détermination d'au moins un signal de pilotage desdits éléments dudit ensemble de restitution par l'application desdits filtres 20 d'adaptation à ladite pluralité de signaux d'entrée d'informations acoustiques (SI); et - une étape de délivrance dudit au moins un signal de pilotage en vue d'une application auxdits éléments de restitution (3v). 2. Procédé selon la revendication 1, caractérisé en ce que ladite étape 25 (10) de détermination de caractéristiques au moins spatiales dudit ensemble de restitution (2), comporte une sous-étape (20) de saisie permettant de déterminer tout ou partie des caractéristiques dudit ensemble de restitution (2). 3. Procédé selon l'une quelconque des revendications 1 et 2, caractérisé en ce que ladite étape (10) de détermination de caractéristiques au moins 30 spatiales dudit ensemble de restitution (2), comporte une étape (30) de calibrage permettant de délivrer tout ou partie des caractéristiques dudit ensemble de restitution (2). 4. Procédé selon la revendication 3, caractérisé en ce que ladite sousétape de calibrage (30) comporte pour au moins un des éléments de restitution (3n,) - une sous-étape d'émission (32) d'un signal spécifique (un(t)) vers le5 dit au moins un élément (3h) dudit ensemble de restitution (2); - une sous-étape d'acquisition (34) de l'onde sonore émise en réponse par ledit au moins un élément (3n); - une sous-étape de transformation (36) desdits signaux acquis en un nombre fini desdits coefficients représentatifs de l'onde sonore émise; et - une sous-étape (39) de détermination de paramètres spatiaux et/ou acoustiques dudit élément (3n) à partir desdits coefficients représentatifs de l'onde sonore émise. 5. Procédé selon l'une quelconque des revendications 3 et 4, caractérisé en ce que ladite sous-étape de calibrage (30) comporte en outre une sous15 étape de détermination de la position dans au moins l'une des trois dimensions de l'espace dudit au moins un élément (3n) dudit ensemble de restitution (2). 6. Procédé selon l'une quelconque des revendications 3 à 5, caractérisé en ce que ladite étape de calibrage (30) comporte une sous-étape de détermination de la réponse en fréquence (H,(f)) dudit au moins un élément (3n) dudit 20 ensemble de restitution (2). 7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que ladite étape (50) de détermination de filtres d'adaptation comprend: - une sous-étape (54) de détermination d'une matrice de décodage (D) représentative de filtres permettant la compensation des altérations de restitution 25 dues aux caractéristiques spatiales dudit ensemble de restitution (2); - une sous-étape (55) de détermination d'une matrice de rayonnement multicanal idéale (S) représentative des directions générales prédéterminées associes à chaque signal d'informations de la pluralité des signaux d'entrée (SI) ; et - une sous-étape (56) de détermination d'une matrice (A) représenta30 tive desdits filtres d'adaptation à partir de ladite matrice de décodage (D) et de ladite matrice de rayonnement multicanal (S). 8. Procédé selon la revendication 7, caractérisé en ce que ladite étape (50) de détermination de filtres d'adaptation comporte une pluralité de sousétapes de calcul (51, 52, 53) permettant de délivrer un ordre (L) limite de la préci- sion spatiale des filtres d'adaptation, une matrice (W) correspondant à une fenêtre spatiale représentative de la répartition dans l'espace de la précision souhaitée lors de la reconstruction du champ acoustique et une matrice (M) représentative du rayonnement de l'ensemble de restitution (2) , ladite sous-étape (54) de 5 calcul de la matrice (2) de décodage (D) étant réalisée à partir des résultats de ces sous-étapes de calcul. 9. Procédé selon l'une des revendications 7 ou 8, caractérisé en ce que les matrices de décodage (D), de rayonnement multicanal idéal S et d'adaptation A sont indépendantes de la fréquence, l'étape (70) de détermination d'au 10 moins un signal de pilotage desdits éléments dudit ensemble de restitution par l'application desdits filtres d'adaptation correspondant à de simples combinaisons linéaires suivies de retard. 10. Procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce que ladite étape (10) de détermination de caractéristiques dudit en15 semble de restitution (2) permet la détermination de caractéristiques acoustiques dudit ensemble de restitution (2) et en ce que ledit procédé comporte une étape (60) de détermination de filtres de compensation de ces caractéristiques acoustiques, ladite étape (70) de détermination d'au moins un signal de pilotage comprenant alors une sous- étape (90) d'application desdits filtres de compensation 20 acoustique. 11. Procédé selon la revendication 10, caractérisé en ce que ladite étape (10) de détermination de caractéristiques acoustiques est adaptée pour délivrer des paramètres représentatifs pour au moins un élément (3,) de sa réponse en fréquence (H,(f)). 12. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce que ladite étape (70) de détermination d'au moins un signal de pilotage comporte une sous-étape d'ajustement de gain et d'application de retards afin d'aligner temporellement le front d'onde des éléments de restitution (3,) en fonction de leur distance par rapport audit point donné (5). 13. Programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 12 lorsque ledit programme est exécuté par un ordinateur. 14. Support amovible du type comprenant au moins un processeur de traitement et un élément de mémoire non volatile, caractérisé en ce que ladite mémoire comprend un programme comprenant des instructions de code pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 12, lorsque ledit processeur exécute ledit programme. 15. Dispositif de pilotage d'un ensemble de restitution (2) d'un champ 5 acoustique comportant une pluralité d'éléments de restitution (3n), comportant des moyens (112) d'entrée d'une pluralité de signaux d'entrée d'informations acoustiques (SI) associés chacun à une direction générale de restitution prédéterminée définie par rapport à un point donné (5), caractérisé en ce qu'il comporte en outre: - des moyens (116) de détermination de caractéristiques au moins spatiales dudit ensemble de restitution (2), permettant la détermination de paramètres représentatifs pour au moins un élément (3h) dudit ensemble de restitution (2) de sa position dans les trois dimensions de l'espace par rapport audit point donné (5); - des moyens (114) de détermination de filtres d'adaptation (A) à partir desdites caractéristiques au moins spatiales dudit ensemble de restitution (2) et des directions générales de restitution prédéterminée associées à ladite pluralité des signaux d'entrée d'informations acoustiques (SI); et - des moyens (114) de détermination d'au moins un signal de pilotage 20 (scn) desdits éléments (3n) dudit ensemble de restitution (2) par l'application desdits filtres d'adaptation (A) à ladite pluralité de signaux d'entrée d'informations acoustiques (SI). 16. Dispositif selon la revendication 15, caractérisé en ce que lesdits moyens de détermination des caractéristiques au moins spatiales dudit ensemble 25 de restitution (2) comportent des moyens (116) de saisie directe desdites caractéristiques. 17. Dispositif selon l'une quelconque des revendications 15 et 16, caractérisé en ce qu'il est adapté pour être associé à des moyens de calibrage (91, 92, 93, 100) permettant la détermination des caractéristiques au moins spatiales 30 dudit ensemble de restitution (2). 18. Dispositif selon la revendication 17, caractérisé en ce que lesdits moyens de calibrage comprennent des moyens d'acquisition d'une onde sonore (100) comportant quatre capteurs de pression disposés selon une forme générale de tétraèdre 19. Dispositif selon l'une quelconque des revendications 15 à 18, caractérisé en ce que lesdits moyens de détermination de caractéristiques sont adaptés pour la détermination de caractéristiques acoustiques d'au moins un desdits éléments (3n) dudit ensemble de restitution (2), ledit dispositif comportant 5 des moyens de détermination de filtres de compensation acoustique à partir desdites caractéristiques acoustiques et lesdits moyens de détermination d'au moins un signal de pilotage étant adaptés pour l'application desdits filtres de compensation acoustique. 20. Dispositif selon la revendication 19, caractérisé en ce que lesdits 10 moyens de détermination des caractéristiques acoustiques sont adaptés pour la détermination de la réponse en fréquence (Hn(f") desdits éléments (3.) de l'ensemble de restitution (2). 21. Appareil de traitement de données audio et vidéo comportant des moyens (112) de détermination d'une pluralité de signaux d'entrée d'informations 15 acoustiques (SI) associés chacun à une direction générale de restitution prédéterminée définie par un point donné (5), caractérisé en ce qu'il comporte en outre un dispositif de pilotage d'un ensemble de restitution (2) selon l'une quelconque des revendications 1 à 19. 22. Appareil selon la revendication 21, caractérisé en ce que lesdits 20 moyens de détermination d'une pluralité de signaux d'entrée sont formés d'une unité (112) de lecture et de décodage des disques audio et/ou vidéo numériques.