FR 2763774 A1 19981127 FR 9706191 A 19970521 PROCEDE ET DISPOSITIF DE REDUCTION D'ECHO ACOUSTIQUE MULTIVOIES ET DE SPATIALISATION SONORE La présente invention concerne un procédé et un dispositif de réduction d'écho acoustique multivoies, et de restitution du relief sonore ou "spatialisation" sonore. Elle trouve une application particulièrement importante dans les systèmes de transmission numérique de signaux sonores sur plusieurs canaux de transmission entre un site local et un ou plusieurs sites distants, chaque site comprenant plu- sieurs microphones et plusieurs haut-parleurs. De tels systèmes sont dits "multivoies" à la prise de son (le son est capté par plusieurs microphones), à la transmission (sur plusieurs canaux) et à la restitution du son (au moyen de plusieurs haut-parleurs). A titre d'exemple non limitatif, la présente invention sera ici décrite dans son application aux salles de visioconférence, aussi bien pour la transmission entre deux salles, dans une configuration dite "point à point" (o l'image de la salle distante est affichée sur un écran, ou bien plusieurs points de vue de la salle distante sont affichés sur plusieurs écrans) que pour la transmission dans un réseau comportant plus de deux salles, dans une configu- ration dite "multipoint" (o les images des salles distantes sont affichées sur plusieurs écrans). Dans une salle comportant N microphones et M haut- parleurs, il existe NxM chemins d'écho acoustique. D'une part, les systèmes classiques d'annulation d'écho acoustique ne permettent pas de réaliser un contrôle de l'écho à un coût raisonnable dans un tel contexte. D'autre part, associer pour la transmission chaque microphone d'une salle locale à un haut-parleur d'une salle distante, afin d'obte- nir la meilleure spatialisation possible du son, multiplie- rait le nombre de canaux de transmission, ce qui entraîne- rait un coût de transmission prohibitif pour un développe- ment commercial. On connait, outre les systèmes classiques d'annula- tion d'écho acoustique, des systèmes de réduction d'écho acoustique fondés sur des variations de niveau des signaux sonores. Un inconvénient majeur généralement rencontré dans de tels systèmes est la limitation de l'interactivité, c'est-à-dire une baisse de qualité de réception, notamment une variation sensible du niveau sonore en réception, dans des situations dites de "double parole", c'est-à-dire en cas de présence simultanée d'un signal sonore local utile et d'un signal sonore utile provenant de la salle distante, dans une salle locale donnée. La présente invention vise à fournir un procédé et un dispositif de réduction d'écho acoustique multivoies et de spatialisation, qui permettent de réduire l'écho tout en préservant l'interactivité et qui, tout en utilisant un nombre relativement faible de canaux de transmission, assurent la restitution du relief sonore. Dans son applica- tion à la visioconférence, la présente invention permet de rapprocher la réunion à distance d'une situation de communi- cation naturelle. Dans ce but, l'invention propose notamment un procédé de réduction d'écho acoustique et de spatialisation, dans un système de transmission numérique de signaux sonores sur P canaux de transmission, entre un site local et au moins un site distant, comportant chacun N microphones et M haut-parleurs, N, M, P étant des entiers, les valeurs de N et M pouvant différer selon les sites, suivant lequel: (a) on calcule un cumulant pour chaque signal microphonique xi(n) du site local, i étant un entier compris entre 1 et N et n désignant le rang temporel des échantil- ions, et on calcule un cumulant pour chaque signal de haut- parleur zj(n) du site local, j étant un entier compris entre 1 et M (la notion de cumulant est définie plus loin); puis, pour tout i, 1 5 i N, (b) on calcule une première atténuation G.i:(i,n) pour le signal microphonique xi(n) du site local à partir d'un rapport entre des cumulants obtenus précédemment; (c) on modifie la première atténuation Gmic(i,n) de façon à obtenir une deuxième atténuation G',mi(i,n) définie comme suit: G'lij(i,n) = S1(G.jó(i,n)) o Sl(Gjmi(i,n)) = s si Gj,(i, n) S s, S1(G.ic(i,n)) = G.mic(i,n) si s qui dans ledit premier cas, est égale à la deuxième atténua- tion G'.ij(i,n) et dans ledit second cas, est égale à la deuxième atténuation G'. c(i,n) dans le calcul de laquelle on augmente le seuil minimal s d'une valeur prédéterminée; (f) on calcule une quatrième atténuation F(i,n) à partir d'un rapport entre des cumulants de signaux micropho- niques; (g) on modifie la quatrième atténuation r(i,n) de façon à obtenir une cinquième atténuation r'(i,n) définie comme suit: r'(i, n) = S2(r(i,n)) o S2(r(i,n)) = s' si r(i,n) s s', S2(r(i,n)) = r(i,n) si s' o ak i(n) désigne des coefficients de codage réels prédé- terminés et gi(n) désigne les coefficients de pondération obtenus précédemment; puis, pour tout entier J, 1 S j S M: (k) on calcule, pour le signal de haut-parleur zj(n) du site distant, une sixième atténuation GHp(j,n) à partir de cumulants calculés pour chaque signal Yk(n) transmis en provenance du site local; (Q) on modifie la sixième atténuation GHP(j,n) de façon à obtenir un coefficient de pondération kX(n) défini comme suit: X%(n) = S3(GHp(j,n)) O S3(GHP(j, n)) = s* si GHP(j,n) s s*, S3(G.p(j,n)) = GHP(j,n) si s* 1, s* étant un seuil minimal prédéterminé strictement inférieur à 1; (m) on détermine le signal de haut-parleur zj(n) du site distant à partir d'une combinaison linéaire des signaux transmis yk(n) pondérés, définie comme suit: P ze(n) = Xk(n).E yjk(n).yk(n) k=l o yJk(n) désigne des coefficients de décodage réels prédé- terminés, et o XA(n) désigne les coefficients de pondéra- tion obtenus précédemment; et (n) on émet sur le jéme haut-parleur du site distant le signal de haut-parleur z3(n) ainsi obtenu. Les opérations ci-dessus sont effectuées de façon similaire dans toutes les salles du réseau considéré, à la fois en émission et en réception. Dans toute la suite, pour simplifier, on considère une salle donnée dite "locale" et une ou plusieurs salles "distantes", et on ne décrit que les opérations effectuées dans la salle locale à l'émission et les opérations effectuées dans la ou les salle(s) distan- te(s) à la réception, bien que des opérations similaires soient également effectuées dans la ou les salle(s) distan- te(s) à l'émission et dans la salle locale à la réception, respectivement. Dans le même but que décrit plus haut, la présente invention propose également un dispositif de réduction d'écho acoustique et de spatialisation, dans un système de transmission numérique de signaux sonores sur P canaux de transmission, entre un site local et au moins un site distant, comportant chacun N microphones et M haut- parleurs, N, M, P étant des entiers, les valeurs de N et M pouvant différer selon les sites, caractérisé en ce qu'il comprend: un module de codage des signaux sonores, recevant en entrée N signaux numériques xi(n) provenant respectivement des N microphones du site local, o i est un entier compris entre 1 et N, et n est un entier qui désigne le rang temporel des échantillons, ce module de codage fournissant en sortie P signaux numériques yk(n), o k est un entier compris entre 1 et P. définis par la formule suivante: N yk(n) = aki(n).li(n).xi(n) i=l o aki(n) désigne des coefficients de codage réels prédéter- minés, et o pi(n) désigne des coefficients de pondération qui dépendent de cumulants des signaux reçus par les microphones du site local et de cumulants des signaux émis par les haut- parleurs du site local, les P signaux yk(n) étant transmis respectivement sur les P canaux de transmission entre le site local et le site distant; et un module de décodage des signaux sonores, recevant en entrée les P signaux yk(n) et fournissant en sortie M signaux numériques ze(n) à émettre respectivement par les M haut- parleurs de chaque site, o j est un entier compris entre 1 et M, définis par la formule suivante: ze(n) = k3(n).E yjk(n).yk(n) k=l o yjk(n) désigne des coefficients de décodage réels prédé- terminés, et o Aj(n) désigne des coefficients de pondération qui dépen- dent des signaux yk(n). D'autres caractéristiques et avantages de la présente invention apparaîtront à la lecture de la descrip- tion détaillée qui suit de modes particuliers de réalisa- tion, donnés à titre d'exemples non limitatifs. La descrip- tion se réfère aux dessins annexés, dans lesquels: - la figure 1 représente de façon schématique et simplifiée une salle de visioconférence classique du type à plusieurs microphones et plusieurs haut-parleurs, dans laquelle la disposition des principaux éléments est donnée à titre d'exemple non limitatif; - la figure 2 est un organigramme général du procédé de la présente invention, dans un mode particulier de réalisation; - la figure 3 représente de façon schématique le dispositif de la présente invention, dans un mode particu- lier de réalisation; - la figure 4 représente de façon schématique un module de codage compris dans le dispositif de la présente invention, dans un mode particulier de réalisation; et - la figure 5 représente de façon schématique un module de décodage compris dans le dispositif de la présente invention, dans un mode particulier de réalisation. Dans toute la suite, on décrit la présente invention dans son application à un réseau d'au moins deux salles de visioconférence. Cette application fournit un exemple non limitatif de système de transmission numérique de signaux sonores entre un site local et un ou plusieurs sites distants comportant chacun plusieurs microphones et plu- sieurs haut-parleurs. Les signaux sonores considérés ici sont généralement des signaux de parole. On suppose dans la présente description que dans une salle donnée, qu'il s'agisse d'une salle locale ou distante, les microphones et les haut-parleurs sont choisis de telle façon que le couplage entre haut-parleurs et microphones est négatif. Cela signifie que la valeur maximale, en dB, de la fonction de transfert entre l'entrée des amplificateurs des haut-parleurs et la sortie des préamplificateurs auxquels sont connectés les microphones est négative. L'énergie d'un signal issu du préamplificateur d'un microphone et provenant d'un haut-parleur de la même salle est donc inférieure à l'énergie du signal en entrée de l'amplificateur de ce haut- parleur. Le couplage négatif assure la stabilité de la boucle acoustique; en particulier, il n'y a pas d'effet Larsen. D'autre part, cela permet de discriminer entre l'écho et la parole distante utile, puisque le niveau moyen de l'écho est toujours inférieur au niveau moyen de la parole distante utile. L'hypothèse du couplage négatif est vérifiée dans la plupart des salles de visioconférence classiques. On appelle "signal microphonique à un instant donné" l'ensemble des signaux sonores captés par un microphone à cet instant. On appelle "signal de haut-parleur à un instant donné" le signal sonore émis par un haut-parleur à cet instant. On rappelle, en relation avec la figure 1, qu'une salle de visioconférence 10 classique du type à plusieurs microphones et plusieurs haut-parleurs comprend plusieurs microphones 12, plusieurs haut-parleurs 14, et au moins un écran 16, en plus du mobilier courant (généralement table(s) et siège(s)). Pour plus de clarté, on a représenté seulement un microphone 12, deux haut-parleurs 14, un écran 16 et une table 18. On considère dans toute la suite une salle locale ayant N microphones et M haut-parleurs, les entiers N et M n'étant pas nécessairement égaux, et une ou plusieurs salles distantes, selon le mode de réalisation envisagé, ayant également N microphones et M haut-parleurs. Les entiers N et M peuvent être différents d'une salle du réseau à une autre. Les salles de visioconférence considérées ne comportent généralement pas plus de 12 microphones et 16 haut-parleurs. On décrit d'abord le procédé de la présente inven- tion, dans un mode particulier de réalisation, à l'aide de la figure 2. Au cours d'une première étape, on calcule des grandeurs statistiques appelées cumulants, sur les signaux de chaque microphone et haut-parleur de la salle locale. Ces grandeurs décrivent les distributions des signaux micropho- niques et de haut-parleur. A titre d'exemple non limitatif, on peut les estimer soit en large bande, soit sur certaines bandes de fréquence. On rappelle que les cumulants sont les dérivées partielles successives de la seconde fonction caractéristi- que d'un vecteur aléatoire, prises à l'origine. En particu- lier, on montre que pour deux signaux a et b, si a = b = x alors le cumulant d'ordre 2 des signaux a et b, appelé dans ce cas auto-cumulant d'ordre 2 du signal x, représente l'énergie du signal x, et si a e b alors le cumulant d'ordre 2 des signaux a et b, appelé dans ce cas inter-cumulant d'ordre 2 des signaux a et b, représente la fonction d'intercorrélation des signaux a et b. Dans un mode particulier de réalisation du procédé de la présente invention, on choisit comme cumulant pour les signaux microphoniques et les signaux de haut-parleur l'auto-cumulant d'ordre 2, c'est-à- dire l'énergie de ces signaux. On note xi(n) le signal capté par le ièm' microphone d'une salle donnée, i étant un entier compris entre 1 et N, et n désignant le rang temporel des échantillons. On note zj(n) le signal émis par le jèm'e haut-parleur d'une salle donnée, j étant un entier compris entre 1 et M. On note Cum2(xi(n)) et Cum2(zj(n)) les énergies respectives des signaux microphoniques et de haut-parleur. Cette première étape est indiquée par le chiffre de référence 1000 sur la figure 2. Dans une variante non détaillée ici, on peut choisir comme cumulants des inter-cumulants d'ordre 2, c'est-à-dire des fonctions d'intercorrélation des signaux microphoniques et de haut-parleur. Puis, comme indiqué en 1001 sur la figure 2, à partir de ces grandeurs statistiques, on calcule pour chaque signal microphonique de la salle locale une première attenuation Gmc(i,n) permettant de réduire l'écho en fonction des signaux de haut-parleurs de la salle locale. Dans le mode particulier de réalisation mentionné précédemment, o les cumulants choisis sont les énergies, on calcule la première atténuation G. ó(i,n) à partir de la formule suivante: Cum-(xi(n)) Geic(i,n) = (1) M C. E Cum2(zj(n)) j=l o C est une constante réelle strictement supérieure à 1. En présence de signaux de haut-parleurs, cette constante permet d'atténuer plus rapidement les signaux microphoniques, ce qui permet de réduire l'écho plus rapidement. Le choix de C est lié en particulier à des caractéristiques acoustiques de la salle et des transducteurs utilisés. A titre d'exemple non limitatif et dans un contexte acoustique particulier, on peut choisir C = 8. On peut calculer les cumulants ainsi que la première atténuation Gió(i, n) soit à chaque échantillon, soit à une cadence inférieure à la fréquence d'échantillonnage, soit par trame de plusieurs échantillons, à titre d'exemples non limitatifs. Dans ce dernier cas, on peut par exemple calculer une valeur moyenne des cumulants et de la première atténuation sur une trame donnée et attribuer cette valeur moyenne à chaque échantillon de cette trame. Cependant, on n'applique pas la première atténuation G.,c(i,n) en l'état au signal microphonique correspondant x,(n). Pour éviter une saturation des signaux microphoni- ques, on modifie la première atténuation G, ,(i,n) de façon qu'elle ait une valeur maximale égale à 1 et qu'elle ait comme valeur minimale l'atténuation souhaitée, qui dépend du niveau de l'écho acoustique. Comme indiqué en 1002 sur la figure 2, on obtient ainsi pour chaque signal microphonique de la salle locale une atténuation modifiée appelée deuxième atténuation, notée G'm.ic(i,n) et définie par: G'.iC(i,n) = Si(Gc(i,n)) (2) o Sl(Gmic(i,n)) = s si Gic(i,n) S s, Sl(G.ji(i,n)) = Gmic(i,n) si s S1(Gió(i,n)) = 1 si Gmió(i,n) 2 1. Le seuil minimal s est l'atténuation souhaitée précitée. Dans l'exemple précité o C = 8, pour une fonction de transfert entre haut- parleurs et microphones ayant une valeur moyenne de -20 dB sur l'ensemble du domaine des fréquences mises en oeuvre, un seuil minimal avantageux est s - 0,1, ce qui correspond, en dB, à 20 loge0s - -20 dB, o log10 désigne le logarithme décimal. Cependant, on n'applique pas non plus la deuxième atténuation G'mie(i,n) en l'état au signal microphonique correspondant xi(n). Sinon, dans le cas de la double parole, c'est-à-dire en présence de parole locale et distante simultanées, le signal de parole local subirait la même atténuation que l'écho, ce qui limiterait l'interactivité. On effectue pour cela un ajustement supplémentaire des atténuations à appliquer aux signaux microphoniques. A partir du rapport de cumulants (1) calculé précédemment, on détermine, microphone par microphone, pour tout entier i compris entre 1 et N, si le signal microphoni- que xi(n) est un signal d'écho seul ou un signal de parole locale seule, suivant un premier cas, ou si xi(n) est un signal de double parole, suivant un second cas. Dans un mode particulier de réalisation, si la première atténuation G.ic(i,n) a une valeur très inférieure à 1, cela signifie qu'on est en présence d'un signal d'écho seul, étant donné qu'on a supposé que le couplage entre haut-parleurs et microphones est négatif. Si G..,(i,n) est très supérieur à 1, on est en présence d'un signal de parole locale seule. Si Gmic(i,n) a une valeur intermédiaire, on est en présence de double parole. Cette détermination du type de chaque signal microphonique peut être effectuée à l'aide de toute autre méthode jugée appropriée. En fonction de la détermination du type de chaque signal microphonique xi(n), on calcule une troisième atténuation G"mió(i,n), comme indiqué en 1003 sur la figure 2, à partir de la deuxième atténuation G'.jó(i,n), comme suit. Dans le premier cas précité, o le signal microphoni- que xi(n) est un signal d'écho seul ou un signal de parole locale seule, la troisième atténuation est G",,,(i,n) = G'., (i,n). Dans le second cas précité, o le signal micro- phonique xi(n) est un signal de double parole, on "relache" l'atténuation, c'est-à-dire que dans le calcul de la deuxième atténuation G'mic(i,n), on remplace la valeur initiale du seuil minimal s par une valeur plus élevée, dans une proportion telle que l'écho soit perceptible, mais non gênant. Dans l'exemple précité o C = 8 et s = 0,1, on peut par exemple choisir d'augmenter le seuil s lorsque le rapport de cumulants (1) est supérieur à -12 dB; on peut choisir de le faire passer à s = 0,7, soit environ -3 dB. Ce relâchement de l'atténuation, effectué uniquement pour les microphones qui captent un signal de parole locale, contri- bue notamment à préserver l'interactivité. Les étapes décrites précédemment concernent des opérations effectuées séparément pour chaque microphone de la salle locale. Afin de prendre en compte dans chacune de ces opérations la présence des autres microphones de la salle locale, on calcule pour chaque signal microphonique xi(n) une atténuation additionnelle, appelée quatrième atténuation et notée r(i,n), comme indiqué en 1004 sur la figure 2. La quatrième atténuation F(i,n) tient compte des signaux microphoniques xq(n), o q est un entier différent de i et compris entre 1 et N. Elle est calculée à partir d'un rapport entre des cumulants de signaux microphoniques. Dans un mode particulier de réalisation, on choisit comme cumulant pour les signaux microphoniques xi(n) l'auto- cumulant d'ordre 2, c'est-à-dire l'énergie de ces signaux, et on calcule la quatrième atténuation r(i,n) pour chaque signal microphonique xi(n) à partir de la formule suivante: (Cum2(xi(n)))N- r(i,n) =(3) N n Cum2(xq(n)) q=l qwi Comme on l'a indiqué plus haut, on peut estimer les auto-cumulants des signaux microphoniques, d'une part, soit en large bande, soit sur certaines bandes de fréquence, et d'autre part, soit à chaque échantillon, soit à une cadence inférieure à la fréquence d'échantillonnage, soit par trame de plusieurs échantillons. Dans le rapport de cumulants (3), il est particuliè- rement avantageux d'utiliser au dénominateur le produit des cumulants des signaux microphoniques autres que le iêm' signal microphonique. Néanmoins, en variante, on peut utiliser la somme des mêmes cumulants. Cependant, on n'applique pas la quatrième atténua- tion r(i,n) en l'état au signal microphonique correspondant xi(n). Pour éviter une saturation des signaux microphoni- ques, on modifie la quatrième atténuation r(i,n) suivant une étape analogue à l'étape de modification de la première atténuation G.i=(i,n), décrite plus haut. Comme indiqué en 1005 sur la figure 2, on obtient ainsi pour chaque signal microphonique de la salle locale une cinquième atténuation r'(i, n), définie par: r'(i,n) = S2(r(i,n)) (4) o S2(r(i,n)) = s' si r(i,n) 5 s', S2(r(i,n)) = r(i,n) si s' S2(r(i,n)) = 1 si r(i,n) > 1. Le seuil minimal s' est l'atténuation souhaitée, qui dépend principalement du niveau de diaphonie perceptible, et dans une moindre mesure, du niveau d'écho. La diaphonie est ici caractérisée par des interférences perceptibles dans le signal microphonique xi(n), entre la parole qui provient du locuteur s'exprimant devant le i*' microphone de la salle locale et la parole qui provient de locuteurs s'exprimant devant les autres microphones de la salle locale et qui est également captée par le i*è' microphone. Dans l'exemple précité o C = 8 et s = 0,1 (ou s = 0,7 en cas de double parole), un seuil s' avantageux est s' = 0,125, soit environ -18 dB. On calcule ensuite le produit des troisième et cinquième atténuations G",ió(i,n) et r'(i,n) obtenues précédemment, de façon à obtenir pour chaque signal micro- phonique xi(n) une atténuation globale G*.ic(i,n) définie par: G*mic(i,n) = G"mió(i,n).r'(i,n). On modifie cette atténuation globale de façon à obtenir pour chaque signal microphonique xi(n), comme indiqué en 1006 sur la figure 2, un coefficient de pondéra- tion gi(n) défini comme suit: gi(n) = S4(G*mic(i,n)) (5) o S4(G*mij(i,n)) = s" si G*,jc(i,n) s s" et S4(G* i(in)) = G*mic(i,n) si s" en cas de double parole) et s' = 0,125, un seuil s" avanta- geux est s" = 0,1. On applique à chaque signal microphonique xi(n) le coefficient de pondération gi(n) correspondant. En variante aux étapes 1005 et 1006, on peut tout d'abord calculer l'atténuation rF"(i,n) puis l'appliquer au signal microphonique correspondant xi(n), de façon à obtenir un signal microphonique atténué; on calcule ensuite l'atténuation G"oic(i,n) et on l'applique au signal micropho- nique atténué. On suppose que P canaux sont prévus pour la trans- mission numérique entre la salle locale et la ou les salles distantes. Dans le mode particulier de réalisation décrit, l'étape suivante du procédé de la présente invention consiste à transmettre sur chaque canal de transmission un signal yk(n), k étant un entier compris entre 1 et P. comme indiqué en 1007 sur la figure 2, sous forme d'une combinai- son linéaire des signaux microphoniques xi(n) pondérés, définie comme suit: N yk(n) = aki(n).pi(n).xi(n) (6) i=l o aki(n) désigne des coefficients de codage réels prédé- terminés et pi(n) désigne les coefficients de pondération obtenus précédemment. Les coefficients gi(n) varient au cours du temps et contribuent notamment à réduire l'écho tout en préservant l'interactivité. Ils permettent en outre de privilégier le ou les microphone(s) actif(s) et de favoriser la spatialisa- tion du son. Les coefficients ak (n), 1 s k s P. 1 s i s N. forment une matrice de codage à P lignes et N colonnes. Le coefficient situé à l'intersection de la k0ol ligne et de la io0èm colonne, noté akOiO(n), o k0 est un entier compris entre 1 et P et i0 est un entier compris entre 1 et N, représente la contribution du io0m' microphone au k0o' canal de transmis- sion. Ainsi, chaque ligne de la matrice de codage fournit les contributions des N microphones à un canal de transmis- sion donné, et chaque colonne fournit les contributions d'un microphone donné aux différents canaux de transmission. La somme des coefficients de chaque colonne vaut 1. Ces contributions, c'est-à-dire les coefficients aSi(n), sont inférieures ou égales à 1 et sont choisies de façon appropriée selon le contexte acoustique et vidéo. Les coefficients aki(n) dépendent notamment de la configuration du réseau de salles de visioconférence. Par exemple, ils diffèrent selon que la configuration est du type point à point ou multipoint (ces notions ont été définies en introduction). Les coefficients a1,(n) dépendent également de la prise de vue effectuée dans les salles de visioconfé- rence du réseau, de façon à assurer la concordance entre les signaux sonores et les images fournies par les écrans des différentes salles. Chaque microphone est généralement associé à une région de la salle locale. Chaque canal transmet générale- ment les contributions de microphones voisins d'une même région. A titred'exemple non limitatif, dans le cas o N = 6 (six microphones) et P = 2 (deux canaux de transmission), la matrice de codage peut être choisie comme suit: 1 1 1/2 1/2 0 0 0 0 1/2 1/2 1 1 On a choisi dans cet exemple d'associer les premier et deuxième microphones à une région "gauche", les deuxième et troisième microphones à une région "centre" et les cinquième et sixième microphones à une région "droite", en fonction de leur emplacement dans la salle. Les signaux des premier et deuxième microphones sont transmis sur le premier canal (a,, = a,,2 = 1). Les signaux des cinquième et sixième microphones sont transmis sur le deuxième canal (a2 5 = a2,6 = 1). Les signaux des troisième et quatrième microphones sont transmis à 50 % sur le premier canal et à 50 % sur le deuxième canal (al 3 = al,4 = a2,3 = a24 = 1/2). Les contributions des cinquième et sixième micropho- nes au premier canal sont considérées comme négligeables, étant donné l'emplacement sensiblement éloigné de ces microphones par rapport à la région "gauche" (al.5 = a,6 = O). De même, les contributions des premier et deuxième microphones au deuxième canal sont considérées comme négligeables, étant donné l'emplacement sensiblement éloigné de ces microphones par rapport à la région "gauche" (a2.1 = a2 2 =). On considère dans ce qui suit la réception dans une salle distante donnée. L'étape suivante du procédé consiste à calculer, pour chaque signal de haut-parleur ze(n) de cette salle distante, une sixième atténuation GHp(j,n), comme indiqué en 1008 sur la figure 2, à partir de cumulants calculés pour chaque signal yk(n) transmis en provenance de la salle locale considérée précédemment. Dans un mode particulier de réalisation, la sixième atténuation GHp(j,n) est issue d'un rapport entre des auto- cumulants d'ordre 2 (énergies) et des inter-cumulants d'ordre 2 (fonctions d'intercorrélation) des signaux transmis yk(n). De même que pour les cumulants mentionnés plus haut, le calcul peut être effectué par bande de fréquence ou en large bande. De même que pour les atténuations calculées à l'émission, on applique à la sixième atténuation GHp(J,n) un seuil minimal, noté s* et strictement inférieur à 1, et un seuil maximal de 1, de façon que cette atténuation ait pour valeur maximale 1 et pour valeur minimale l'atténuation souhaitée en fonction de la diaphonie et de l'écho. Comme indiqué en 1009 sur la figure 2, on obtient ainsi un coefficient de pondération kX(n) défini comme suit: kX(n) = S3(Gp(j,n)) (7) O S3(G.p(j,n)) = s* si GHp(j,n) s s*, S3(GHp(j,n)) = GHp(j,n) si s* S3(GHp(j,n)) = 1 si GHp(j,n) 2 1. A titre d'exemple non limitatif, dans le cas o P = 2 (deux canaux de transmission) et M = 3 (trois haut- parleurs), les coefficients de pondération peuvent être obtenus comme suit: X2(n) = S3(GHP(2,n)), o 2.Cum2(y1(n)).Cum2(y2(n)) Gp(2,n) = Cum2(yl(n))2 + Cum2(y2(n))2 o Cum2(yl(n)) est un auto-cumulant d'ordre 2 qui représente l'énergie du signal transmis yl(n), et o Cum2(y2(n)) est un auto-cumulant d'ordre 2 qui représente l'énergie du signal transmis y2(n). Dans le même exemple, on peut choisir kl(n) = X3(n) = I + s* - k2(n). Ces expressions se généralisent aisément à un nombre plus élevé de canaux de transmission et de haut-parleurs. Dans l'exemple précité o C = 8, s = 0,1 (ou s = 0,7 en cas de double parole) et s' = 0,125, un seuil s* avanta- geux est s* = s' = 0,125, soit environ -18 dB. Comme indiqué en 1010 sur la figure 2, on détermine ensuite chaque signal de haut-parleur z](n) de la salle distante, à partir d'une combinaison linéaire des signaux transmis yk(n) pondérés, définie comme suit: P z](n) = Xj(n).E y k(n).yk(n) (8) k=l1 o yjk(n) désigne des coefficients de décodage réels prédéterminés et X3(n) désigne les coefficients de pondéra- tion obtenus précédemment. Les coefficients kX(n) varient au cours du temps et contribuent notamment à améliorer la restitution du relief sonore et à réduire la diaphonie éventuelle, en particulier en cas de transmission sur un nombre P de canaux inférieur au nombre M de haut-parleurs. Ils contribuent également à réduire l'écho acoustique en limitant le nombre de haut- parleurs actifs à tout moment. Les coefficients yjk(n), 1 s j s M, 1 s k s P, forment une matrice de décodage à P lignes et M colonnes. Ils dépendent notamment de la disposition des haut-parleurs et du mode d'affichage des images dans les salles de visio- conférence du réseau. A la différence des coefficients de codage aki(n), les coefficients yjk(n) sont indépendants de la configuration point à point ou multipoint du réseau. L'étape suivante du procédé consiste à émettre, au moyen du jè" hautparleur de la salle distante, pour tout entier j compris entre 1 et M, le signal de haut-parleur z](n) obtenu à partir de la formule (8). Dans un mode particulier de réalisation, on transmet en outre sur les P canaux de transmission les produits yj.k(n).kJ(n), 1 s j S M, 1 s k ce signal doit être affecté avant émission par le haut- parleur correspondant de la salle distante. Cela revient à effectuer une détection d'activité vocale et à transmettre l'information correspondante. Dans ce mode de réalisation, on adapte la matrice de codage en fonction des locuteurs actifs et on calcule la matrice de décodage à l'émission, à partir de la matrice de codage et des caractéristiques de restitution sonore. Lorsque le nombre de microphones actifs, tel que déduit de la valeur de la cinquième atténuation r'(i,n) ou de la troisième atténuation G"., c(i,n), est inférieur ou égal au nombre de canaux de transmission, chaque signal issu d'un microphone actif est aiguillé vers un des canaux. La matrice de codage est alors constituée uniquement de 1 et de 0 (les 1 correspondant aux microphones actifs). * Dans ce mode de réalisation, la matrice de décodage est la transposée de la matrice de codage. On peut transmettre les M x P coefficients yjk(n) au moyen d'un faible nombre de bits. A titre d'exemple non limitatif, trois bits par coefficient peuvent suffire. On peut par exemple effectuer un codage de type logarithmique o chaque coefficient, exprimé en dB, est quantifié de façon uniforme sur un nombre de niveaux donné. A titre d'exemple non limitatif, on peut prévoir huit niveaux. Dans un mode particulier de réalisation, on diffuse en outre, dans la salle locale, les signaux à transmettre Yk(n), au moyen des haut- parleurs locaux, en plus des signaux reçus en provenance de la ou des salle(s) distan- te(s) et à un niveau plus faible que les signaux reçus. Cette diffusion d'un retour local amène les locuteurs à parler moins fort ou plus loin du microphone, au cas o ils parlaient initialement à un niveau très élevé ou trop près du microphone sans s'en rendre compte. Cette diffusion contribue également à créer un environnement acoustique commun dans une salle donnée, ce qui est particulièrement avantageux pour des locuteurs situés à des emplacements éloignés d'une même salle. On décrit dans ce qui suit le dispositif de la présente invention, dans un mode particulier de réalisation, à l'aide de la figure 3. Ce dispositif s'applique à un système de transmis- sion numérique de signaux sonores du même type que celui considéré en relation avec le procédé, à savoir, un système mettant en oeuvre P canaux de transmission entre une salle locale et au moins une salle distante, qui comportent chacune N microphones et M haut-parleurs, les valeurs de N et M pouvant différer selon les salles. Comme le montre la figure 3, le dispositif comprend un module 20 de codage des signaux sonores. Le module 20 reçoit en entrée N signaux numériques xi(n) qui proviennent respectivement des N microphones de la salle locale. Il fournit en sortie P signaux numériques yk(n) définis par la formule (6) donnée plus haut. Les P signaux yk(n) sont transmis respectivement sur les P canaux de transmission dans le réseau 22. Le dispositif comprend également un module 24 de décodage des signaux sonores. Le module 24 reçoit en entrée les P signaux yk(n) et fournit en sortie M signaux numéri- ques ze(n), définis par la formule (8) donnée plus haut, à émettre respectivement par les M haut-parleurs de chaque salle. Dans le mode particulier de réalisation représenté sur la figure 4, le module 20 de codage des signaux sonores comprend un module 26 pour calculer un cumulant pour chaque signal microphonique xi(n) de la salle locale et un cumulant pour chaque signal de haut-parleur zj(n) de la salle locale. Le module 20 comprend également un module 28, connecté en sortie du module 26, pour calculer la première atténuation G,ic(i,n), décrite plus haut, pour chaque signal microphonique xi(n) de la salle locale, à partir d'un rapport entre des cumulants fournis par le module 26 de calcul de cumulants. Dans un mode particulier de réalisa- tion, le module 28 fournit en sortie des valeurs de G.îó(i,n) telles que définies par la formule (1) donnée plus haut. Un module 30 pour modifier chaque première atténua- tion G. i(i,n) est connecté en sortie du module 28. Il fournit en sortie la deuxième atténuation G'oIc(i,n) définie par la formule (2) donnée plus haut. Le module 20 comprend de plus un module 32, connecté en sortie du module 26 de calcul de cumulants, pour détermi- ner, pour chaque signal microphonique xi(n), à partir des cumulants de signaux microphoniques et de haut-parleurs fournis par le module 26, si ce signal microphonique est un signal d'écho seul ou un signal provenant uniquement de la salle locale, suivant un premier cas, ou si ce signal microphonique présente des composantes provenant de la salle locale et d'autres composantes provenant de la salle distante, suivant un second cas, comme on l'a décrit plus haut. Un module 34 pour calculer la troisième atténuation G".i(i,n) définie plus haut est connecté en sortie des modules 30 et 32. Le module 20 comprend en outre un module 36, connecté en sortie du module 26 de calcul de cumulants, pour calculer la quatrième atténuation r(i,n) à partir d'un rapport entre des cumulants de signaux microphoniques. Dans un mode particulier de réalisation, le module 36 fournit en sortie des valeurs de r(i,n) telles que définies par la formule (3) donnée plus haut. Un module 38 pour modifier chaque quatrième atténua- tion r(i,n) est connecté en sortie du module 36. Il fournit en sortie la cinquième atténuation r' (i,n) définie par la formule (4) donnée plus haut. Un module 40 pour calculer et modifier le produit des troisième et cinquième atténuations G"mic(i,n) et r'(i,n) est connecté en sortie des modules 34 et 38. Le module 40 fournit en sortie les coefficients de pondération i(n) définis par la formule (5) donnée plus haut. Dans le mode particulier de réalisation représenté sur la figure 5, le module 24 de décodage des signaux sonores comprend un module 42 pour calculer, pour chaque signal de haut-parleur ze(n) de la salle distante considé- rée, la sixième atténuation GHP(j,n) définie plus haut, à partir de cumulants calculés pour chaque signal transmis yk(n). Ces cumulants peuvent être fournis par le module 26 de calcul de cumulants représenté sur la figure 4. Un module 44 pour modifier chaque sixième atténua- tion GHP(j,n) est connecté en sortie du module 42. Il fournit en sortie les coefficients de pondération kX(n) définis par la formule (7) donnée plus haut. REVENDICATIONS 1. Procédé de réduction d'écho acoustique et de spatialisation sonore, dans un système de transmission numérique de signaux sonores sur P canaux de transmission, entre un site local et au moins un site distant, comportant chacun N microphones et M haut-parleurs, N, M, P étant des entiers, les valeurs de N et M pouvant différer selon les sites, suivant lequel: (a) on calcule un cumulant pour chaque signal microphonique xi(n) du site local, i étant un entier compris entre 1 et N et n désignant le rang temporel des échantil- lons, et on calcule un cumulant pour chaque signal de haut- parleur zj(n) du site local, j étant un entier compris entre 1 et M; puis, pour tout i, 1 5 i S N, (b) on calcule une première atténuation Gjic(i,n) pour le signal microphonique xi(n) du site local à partir d'un rapport entre des cumulants obtenus précédemment; (c) on modifie la première atténuation Gmic(i,n) de façon à obtenir une deuxième atténuation G'.1ó(i,n) définie comme suit: G'.ic(i,n) = Sl(G. ic(i,n)) o Sl(G.mi(i,n)) = s si G.mic(i,n) s, Sl(G.ó(i,n)) = G.. i(i,n) si s 1, s étant un seuil minimal prédéterminé strictement inférieur à 1; (d) à partir des cumulants de signaux microphoniques et de haut-parleurs calculés précédemment, on détermine si le signal microphonique xi(n) est un signal d'écho seul ou un signal provenant uniquement du site local, suivant un premier cas, ou si le signal microphonique xi(n) présente des composantes provenant du site local et d'autres composantes provenant du site distant, suivant un second cas; (e) on calcule une troisième atténuation G"mic(i,n), qui dans ledit premier cas, est égale à la deuxième atténua- tion G'mic(i,n) et dans ledit second cas, est égale à la deuxième atténuation G',ic(i,n) dans le calcul de laquelle on augmente le seuil minimal s d'une valeur prédéterminée; (f) on calcule une quatrième atténuation r(i,n) à partir d'un rapport entre des cumulants de signaux micropho- niques; (g) on modifie la quatrième atténuation r(i,n) de façon à obtenir une cinquième atténuation r' (i,n) définie comme suit: F'(i,n) = S2(r(i,n)) o S2(r(i,n)) = s' si r(i,n) s 1, s' étant un seuil minimal prédéterminé strictement inférieur à 1; (h) on calcule le produit des troisième et cin- quième atténuations G"mic(i,n) et r' (i,n) obtenues précédem- ment, de façon à obtenir une atténuation globale G*.ió(i,n) définie par: G*mi.c(i,n) = G".ió(i,n).r'(i,n) (i) on modifie l'atténuation globale G*mic(i,n) de façon à obtenir un coefficient de pondération Pi(n) défini comme suit: PI(n) = S4(G*.ó(i,n)) O S4(G*.mic(i,n)) = s" si G*.ji(i, n) à 1; (j) on transmet sur chaque canal de transmission un signal yk(n), k étant un entier compris entre 1 et P. sous forme d'une combinaison linéaire des signaux microphoniques xi(n) pondérés, définie comme suit: N yk(n) = ú akl(n).Pi(n).xi(n) i=l o ak,.(n) désigne des coefficients de codage réels prédé- terminés et pi(n) désigne les coefficients de pondération obtenus précédemment; puis, pour tout entier j, 1 _ j _ M: (k) on calcule, pour le signal de haut-parleur zj(n) du site distant, une sixième atténuation GHp(j,n) à partir de cumulants calculés pour chaque signal Yk(n) transmis en provenance du site local; (Q) on modifie la sixième atténuation GHp(j,n) de facon à obtenir un coefficient de pondération Aj(n) défini comme suit: XA(n) = S3(GHp(j,n)) O S3(GHp(j,n)) = s* si GHp(j,n) 5 s*, S3(GHp(j,n)) = GHp(j,n) si s* 1, s* étant un seuil minimal prédéterminé strictement inférieur à 1; (m) on détermine le signal de haut-parleur zt(n) du site distant à partir d'une combinaison linéaire des signaux transmis yk(n) pondérés, définie comme suit: p zj(n) = kj(n).E yj.k(n).yk(n) k=l o yjk(n) désigne des coefficients de décodage réels prédé- terminés, et o kA(n) désigne les coefficients de pondéra- tion obtenus précédemment; et (n) on émet sur le jr" haut- parleur du site distant le signal de haut-parleur zj(n) ainsi obtenu. 2. Procédé selon la revendication 1, suivant lequel, à l'étape (a), on choisit comme cumulant pour les signaux microphoniques xi(n) et les signaux de haut-parleur z](n) l'énergie de ces signaux, et à l'étape (b), on calcule la première atténuation Gj ó(i,n) pour chaque signal microphoni- que xi(n) à partir de la formule suivante: Cum2(xi(n)) Gmic(i,n) = M C. E Cum2(z3(n)) j=1 o Cum2(.) désigne l'énergie et o C est une constante réelle prédéterminée strictement supérieure à 1. 3. Procédé selon la revendication 1, suivant lequel, à l'étape (a), on choisit comme cumulant pour les signaux microphoniques xi(n) l'énergie de ces signaux, et à l'étape (f), on calcule la quatrième atténuation F(i,n) pour chaque signal microphonique xi(n) à partir de la formule suivante: (Cum2(xi(n)))- r(i,n) = N F Cum2(xq(n)) q=l qui o Cum2(.) désigne l'énergie et o q est un entier. 4. Procédé selon la revendication 1, 2 ou 3, suivant lequel on transmet en outre sur les P canaux les produits Yj.k(n).kj(n), 1 S j s M, 1 s k s P. à la même cadence que les signaux sonores. 5. Procédé selon l'une quelconque des revendications 1 à 4, suivant lequel on diffuse en outre sur le site local les signaux à transmettre Yk(n), au moyen des haut-parleurs locaux, en plus des signaux reçus en provenance du site distant et à un niveau plus faible que ces derniers. 6. Dispositif de réduction d'écho acoustique et de spatialisation sonore, dans un système de transmission numérique de signaux sonores sur P canaux de transmission, entre un site local et au moins un site distant, comportant chacun N microphones et M haut-parleurs, N, M, P étant des entiers, les valeurs de N et M pouvant différer selon les sites, caractérisé en ce qu'il comprend: des moyens (20) de codage des signaux sonores, recevant en entrée N signaux numériques xi(n) provenant respectivement des N microphones du site local, o i est un entier compris entre 1 et N, et n est un entier qui désigne le rang temporel des échantillons, lesdits moyens (20) de codage fournissant en sortie P signaux numériques yk(n), o k est un entier compris entre 1 et P. définis par la formule suivante: N yk(n) = aki(n).pi(n).xi(n) i=l o akl(n) désigne des coefficients de codage réels prédéter- minés, et o P3(n) désigne des coefficients de pondération qui dépendent de cumulants des signaux reçus par les microphones du site local et de cumulants des signaux émis par les haut- parleurs du site local, les P signaux Yk(n) étant transmis respectivement sur les P canaux de transmission entre le site local et le site distant; et des moyens (24) de décodage des signaux sonores, recevant en entrée les P signaux yk(n) et fournissant en sortie M signaux numériques zl(n) à émettre respectivement par les M haut- parleurs de chaque site, o j est un entier compris entre 1 et M, définis par la formule suivante: P z3(n) = kj(n).E Yj.k(n).yk(n) k=l o yj.k(n) désigne des coefficients de décodage réels prédéterminés, et o Xj(n) désigne des coefficients de pondération qui dépen- dent des signaux yk(n). 7. Dispositif selon la revendication 6, caractérisé en ce que lesdits moyens (20) de codage des signaux sonores comprennent: des moyens (26) pour calculer un cumulant pour chaque signal microphonique xi(n) du site local et un cumulant pour chaque signal de haut-parleur z](n) du site local; des moyens (28) pour calculer une première atténua- tion Gmi(i,n) pour chaque signal microphonique xi(n) du site local à partir d'un rapport entre des cumulants fournis par lesdits moyens (26) de calcul de cumulants; des moyens (30) pour modifier chaque première atténuation Gjic(i,n) de façon à obtenir, pour tout i, 1 s i N, une deuxième atténuation G'mic(i,n) définie comme suit: G',ic(i,n) = Sl(G. ,(i,n)) o Sl(Gmic(i,n)) = s si Gmic(i,n) s s, S1(G.ió(i,n)) = Gmic(i,n) si s et dans ledit second cas, est égale à la deuxième atténua- tion G'mic(i,n) dans le calcul de laquelle on augmente le seuil minimal s d'une valeur prédéterminée; des moyens (36) pour calculer, pour tout i, 1 5 i N, une quatrième atténuation r(i,n) à partir d'un rapport entre des cumulants de signaux microphoniques; des moyens (38) pour modifier chaque quatrième atténuation r(i,n) de façon à obtenir une cinquième atténua- tion r'(i,n) définie comme suit: r'(i,n) = S2(r(i,n)) o S2(r(i,n)) = s' si r(i,n) 1, s' étant un seuil minimal prédéterminé strictement inférieur à 1; et des moyens (40) pour calculer et modifier, pour tout i, 1 s i tions G".i(i,n) et r'(i,n) obtenues précédemment, de façon à obtenir lesdits coefficients de pondération gi(n), définis comme suit: Pl(n) = S4(G*miC(i,n)) o G*.li(i,n) = G"lmic(i,n).r'(i,n), o S4(G*,ij(i,n)) = s" si G*jic(i,n) 1, s* étant un seuil minimal prédéterminé strictement inférieur à 1. 9. Dispositif selon la revendication 7, caractérisé en ce que lesdits moyens (28) pour calculer chaque première atténuation G.mic(i,n) fournissent en sortie: Cum2(xi (n)) Gmic(i,fn) = M C. E Cum2(zj(n)) j=1 o Cum2(.) désigne le cumulant qui représente l'énergie et o C est une constante réelle prédéterminée strictement supérieure à 1. 10. Dispositif selon la revendication 7, caractérisé en ce que lesdits moyens (36) pour calculer chaque quatrième atténuation r(i,n) fournissent en sortie: (Cum2(xi(n)))"-: r(i,n) = N n Cum2(xq(n)) q=l qxi o Cum2(.) désigne le cumulant qui représente l'énergie et o q est un entier.