FR 3081662 A1 20191129 FR 1855876 A 20180628 Procédé pour une restitution sonore spatialisée d'un champ sonore audible sélectivement dans une sous-zone d'une zone [001] L'invention se place dans le domaine de l'audio spatialisée et du contrôle de champs sonores. Le procédé a pour but de reproduire différents champs sonores dans une même zone, en utilisant le même système de haut-parleurs. Plus précisément, le procédé cherche à rendre un champ sonore audible dans une première sous-zone et inaudible dans une deuxième sous-zone. Le procédé est notamment basé sur la maximisation du contraste de niveau sonore entre les deux sous-zones de la zone. [002] Actuellement, les méthodes de contrôle de champs sonores permettent de réaliser un contraste important du niveau sonore entre les différentes zones d'écoute. Toutefois, elles ne prennent pas en compte le champ sonore induit dans la sous-zone dans laquelle le champ sonore est à rendre inaudible. En effet, le champ sonore cible à reproduire dans cette sous-zone est caractérisé par une pression acoustique nulle. Ainsi, en ne tenant pas compte du champ sonore induit, le champ sonore réellement reproduit dans cette sous-zone n'est jamais caractérisé par une pression acoustique nulle, bien que le contraste de niveau sonore entre les deux sous-zones soit très important. On connaît notamment le document « Personal Sound Zones : Delivering interface-free audio to multiple listeners », T. Betlehem, W. Zhang, M. A. Poletti, et T. Abhayapala, IEEE, n° 1053-5888/15, p. 81-91, mars 2015. De telles méthodes ne permettent pas d'obtenir un champ sonore homogène sur l'ensemble de la sous-zone dans lequel le champ sonore ne doit pas être reproduit. [3] La présente invention vient améliorer la situation. [4] A cet effet, elle propose un procédé assisté par des moyens informatiques, pour une restitution sonore spatialisée à partir d'un réseau de haut-parleurs, en vue d'une diffusion d'un champ sonore choisi, audible sélectivement dans une première sous-zone prédéterminée d'une zone couverte par les haut-parleurs, ladite zone comportant en outre une deuxième sous-zone dans laquelle on cherche à rendre inaudible ledit champ sonore choisi, dans lequel, les haut-parleurs étant alimentés par des signaux de commande respectifs pour émettre chacun un signal audio en continu, le procédé comprend itérativement et en continu : - estimation d'une pression acoustique dans la deuxième sous-zone, au moins en fonction des signaux de commande respectifs des haut-parleurs, et d'un poids initial respectif des signaux de commande des haut-parleurs, - calcul d'une erreur entre la pression acoustique estimée et une pression acoustique cible, souhaitée dans la deuxième sous-zone, 2 - et, calcul et application de poids respectifs aux signaux de commande des haut-parleurs, en fonction de ladite erreur, le calcul de la pression acoustique dans la deuxième sous-zone étant mis en oeuvre à nouveau en fonction des signaux de commande respectifs, pondérés, des haut-parleurs. 5 [5] Le procédé ne se base pas uniquement sur le contraste sonore souhaité entre le niveau sonore des deux sous-zones. Cela permet, à chaque itération, de prendre en compte le champ sonore effectivement restitué dans la deuxième sous-zone. Ainsi, à chaque nouvelle itération, les signaux de commande des haut-parleurs sont pondérés peu à peu, de manière à ce que le champ sonore dans la deuxième sous-zone, caractérisé par sa 10 pression acoustique, atteigne la valeur du champ sonore cible, caractérisé par sa pression acoustique cible. [6] Selon une réalisation, le procédé comporte préalablement une détermination de distances entre au moins un point de la deuxième sous-zone et des positions respectives de haut-parleurs, pour en déduire des fonctions de transfert acoustiques respectives des haut-parleurs en ledit point, 15 ledit point correspondant ainsi à une position de microphone virtuel, l'estimation de la pression acoustique dans la deuxième sous-zone étant alors effectuée en tenant compte d'un poids initial appliqué à chaque haut-parleur, et de la fonction de transfert propre à chaque haut-parleur en le microphone virtuel dans la deuxième sous-zone, et comportant en outre une détermination de distances entre au moins un point de la première sous-zone et des 20 positions respectives de haut-parleurs, pour en déduire des fonctions de transfert acoustiques respectives des haut-parleurs en ledit point de la première sous-zone, ledit point correspondant ainsi à une position d'un deuxième microphone virtuel. [7] Selon une variante, le point de la deuxième sous-zone peut correspondre à un microphone réel. 25 [8] La création d'un microphone virtuel dans chacune des sous-zones permet d'estimer les fonctions de transferts de chaque couple haut-parleur/microphone virtuel, de sorte que la pression acoustique est estimée de manière précise. 30 [009] Selon une réalisation, on prévoit une pluralité de points au moins dans la deuxième sous-zone formant les positions respectives d'une pluralité de microphones virtuels, pour estimer une pluralité de pressions acoustiques respectives dans la deuxième sous-zone en tenant compte du poids respectif appliqué à chaque haut-parleur et des fonctions de transfert propres à chaque haut-parleur en chaque microphone virtuel. 3 [10] De cette manière, la pression acoustique est estimée en une pluralité de points de la deuxième sous-zone. Cela permet d'appliquer des poids à chaque haut-parleur en tenant compte des écarts de pressions acoustiques pouvant survenir en différents points de la deuxième sous-zone. L'estimation de la pression acoustique est donc réalisée de manière homogène et précise sur l'ensemble de la zone. La précision du 5 procédé en est alors augmentée de manière considérable. [11] Selon une réalisation, la pression cible dans la deuxième sous-zone est évaluée en fonction d'un contraste sonore souhaité entre la première sous-zone et la deuxième sous-zone, et comportant une estimation de : 10 - la pression acoustique dans la première sous-zone en tenant compte du poids respectif appliqué à chaque haut-parleur et de la fonction de transfert propre à chaque haut-parleur en le au moins un deuxième microphone virtuel, - et de là, le contraste sonore en fonction de la pression acoustique estimée dans la première sous-zone et la pression acoustique estimée dans la deuxième sous-zone. 15 [12] Ainsi, la pression acoustique cible dans la deuxième sous-zone est déterminée en tenant compte, à chaque itération, de la pression acoustique dans chacune des sous-zones. Le niveau de contraste désiré est alors contrôlé à chaque itération, ce qui a l'avantage de prendre en compte d'éventuels effets de masquage perceptif fréquentiel du système auditif. Le but n'est plus alors d'uniquement maximiser le contraste sonore 20 entre les deux sous-zones, mais bien de contrôler le contraste sonore souhaité en fonction des pressions acoustiques effectivement estimées dans chaque sous-zone. [13] Selon une réalisation, à chaque itération, un coefficient d'atténuation, défini par la différence entre le contraste sonore souhaité et le contraste sonore estimé entre la première sous-zone et la deuxième sous-zone, 25 est appliqué à la pression acoustique dans la deuxième sous-zone, la pression cible dans la seconde sous-zone prenant alors la valeur de la pression acoustique atténuée de la deuxième sous-zone. [14] Ainsi, le procédé ne cherche pas à produire une pression acoustique nulle dans la deuxième sous-zone, 30 mais au contraire à spécifier le champ sonore cible, caractérisé par sa pression acoustique, dans la deuxième sous-zone. L'efficacité du procédé en est augmentée. De plus, dans le cas d'une pluralité de microphones virtuels, le coefficient d'atténuation est défini en chaque point correspondant aux microphones virtuels. Le procédé permet alors d'obtenir un contraste sonore homogène entre les deux sous-zones, en tous points des deux sous-zones. 4 [015] Selon une réalisation, le procédé comporte préalablement, pour chaque auditeur : - obtenir des données de géolocalisation d'un auditeur dans la zone au moyen d'un capteur de position, - définir la première sous-zone comme correspondant auxdites données de géolocalisation ainsi que la 5 position dudit au moins un deuxième microphone virtuel dans ladite première sous-zone, - définir la deuxième sous-zone comme étant complémentaire de la première sous-zone ainsi que la position dudit au moins un premier microphone virtuel dans ladite deuxième sous-zone. [16] Le procédé permet donc de personnaliser la disposition des deux sous-zones en fonction de la position 10 d'un auditeur. [17] Selon une réalisation, le procédé comporte en outre, pour chaque auditeur, - obtenir répétitivement des données de géolocalisation de l'auditeur dans la zone au moyen d'un capteur de position, et en cas de changement des données de géolocalisation : 15 > définir dynamiquement la première sous-zone comme correspondant auxdites données de géolocalisation ainsi que la position dudit au moins un deuxième microphone virtuel dans ladite première sous-zone, > définir dynamiquement la deuxième sous-zone comme étant complémentaire de la première sous-zone ainsi que la position dudit au moins un premier microphone virtuel dans ladite deuxième sous-zone, 20 distances entre ledit au moins un premier microphone virtuel et les positions respectives des haut-parleurs étant calculées dynamiquement pour en déduire les fonctions de transfert acoustiques respectives des haut-parleurs en ledit au moins un premier microphone virtuel, et l'estimation de la pression acoustique dans la deuxième sous-zone étant alors effectuée en tenant compte du poids respectif appliqué à chaque haut-parleur, 25 - chaque poids respectif appliqué à un haut-parleur comprenant un premier terme correspondant à un poids respectif calculé à une itération précédente et un deuxième terme fonction de l'erreur, un coefficient d'oubli atténuant le terme correspondant au poids calculé à l'itération précédente étant alors appliqué à chaque poids respectif. 30 [018] Ainsi, il est possible de définir dynamiquement les deux sous-zones en fonction de la position d'un auditeur. Plus précisément, la première sous-zone, dans laquelle le champ sonore est audible, suit perpétuellement l'auditeur. L'application d'un poids comprenant un coefficient d'oubli permet de ne pas relancer le procédé à chaque déplacement de l'auditeur, ce qui renforce grandement son efficacité. les 5 [19] Selon une réalisation, le procédé comporte en outre, itérativement et en continu : - estimation de la pression acoustique dans la première sous-zone, en fonction desdits signaux de commande respectifs des haut-parleurs en le au moins un deuxième microphone virtuel et du poids respectif des signaux de commande des haut-parleurs, 5 - calcul d'une dégradation entre la pression acoustique estimée dans la première sous-zone et une pression acoustique cible, souhaitée dans la première sous-zone, ladite pression acoustique cible étant préalablement définie, - et, application de poids respectifs aux signaux de commande des haut-parleurs, en fonction de ladite dégradation, afin que la pression acoustique dans la première sous-zone atteigne la pression cible, 10 - le calcul de la pression acoustique dans la première sous-zone et dans la deuxième sous-zone étant respectivement mis en oeuvre à nouveau en fonction des signaux de commande respectifs, pondérés, des haut-parleurs. [20] Le champ sonore dans la première sous-zone est également contrôlé en continu. De cette manière, le 15 champ sonore cible, caractérisé par la pression acoustique cible dans la sous-zone, peut être continuellement atteint dans la première sous-zone. Les dégradations sonores sont compensées à chaque itération. [21] Selon une réalisation, le nombre d'itérations du procédé, pour que la pression acoustique de la première sous-zone atteigne la pression acoustique cible, souhaitée dans la première sous-zone, est compris entre 10 et 20 1000, de préférence entre 200 et 800, et de préférence entre 400 et 500. [22] Par « atteindre la pression acoustique cible », on entend que l'erreur entre la pression acoustique de la première sous-zone et la pression acoustique cible, souhaitée dans la première sous-zone est négligeable. Notamment, l'erreur peut être définie à partir de données psycho-acoustiques d'utilisateurs. Ces données 25 psycho-acoustiques peuvent être obtenues par des mesures visant à déterminer un seuil de perception pour des utilisateurs. Notamment, le seuil de perception peut évoluer en fonction des signaux sonores constituant la scène sonore ou de la fréquence des signaux sonores. [23] Il apparaît que le procédé présente une très grande efficacité, compte tenu du petit nombre d'itérations 30 avant que le procédé ne converge. [24] Selon une réalisation, le champ sonore choisi est un premier champ sonore, un deuxième champ sonore choisi étant en outre diffusé, audible sélectivement dans la deuxième sous-zone, et dans lequel on cherche à rendre le deuxième champ sonore inaudible dans la première sous-zone, les haut-parleurs étant alimentés par 6 des premiers signaux de commande respectifs pour émettre chacun un signal audio en continu correspondant au premier champ sonore choisi, et par des deuxièmes signaux de commande respectifs pour émettre chacun un signal audio en continu correspondant au deuxième champ sonore choisi, les étapes du procédé selon l'invention étant appliquées à la première sous-zone, de manière à rendre 5 inaudible le deuxième champ sonore dans la première sous-zone. [025] Ainsi, le procédé peut également être mis en oeuvre dans le cas où deux champs sonores sont diffusés dans une même zone. Le procédé peut également être mis en oeuvre pour plus de deux champs sonores. 10 [026] Selon une réalisation, chaque sous-zone comprend au moins un microphone virtuel et deux haut- parleurs, et de préférence chaque sous-zone comprend au moins une dizaine de microphones virtuels et au moins une dizaine de haut-parleurs. [027] L'invention vise également un système de restitution sonore spatialisée à partir d'un réseau de haut- 15 parleurs, en vue d'une diffusion d'un champ sonore choisi, audible sélectivement dans une première sous- zone prédéterminée d'une zone couverte par les haut-parleurs, ladite zone comportant en outre une deuxième sous-zone dans laquelle on cherche à rendre inaudible ledit champ sonore choisi, comprenant une unité de traitement adaptée pour le traitement et la mise en oeuvre du procédé selon l'invention. 20 [028] D'autres avantages et caractéristiques de l'invention apparaîtront à la lecture de la description détaillée ci-après d'exemples de réalisation de l'invention, et à l'examen des dessins annexés sur lesquels : - la figure 1 représente un exemple de système selon un mode de réalisation de l'invention, - les figures 2a et 2b illustrent, sous la forme d'un ordinogramme, les principales étapes d'un mode de réalisation particulier du procédé, 25 - la figure 3 illustre de façon schématique un mode de réalisation dans lequel les deux sous-zones sont définies dynamiquement en fonction des données de géolocalisation d'un auditeur. [029] Les modes de réalisation décrits en référence aux figures peuvent être combinés. 30 [030] La figure 1 illustre schématiquement un système SYST selon un exemple de réalisation. Le système SYST comprend un réseau de haut-parleurs HP comprenant N haut-parleurs (HP1,...,HPN), avec N au moins égal à 2, et de préférence au moins égal à 10. Le réseau de haut-parleurs HP couvre une zone comprenant une première sous-zone SZ1 et une deuxième sous-zone 5Z2. Les haut-parleurs HP sont alimentés par des signaux de commande respectifs pour émettre chacun un signal audio en continu, en vue d'une diffusion sonore 7 spatialisée d'un champ sonore choisi. Le champ sonore choisi est à rendre audible dans une des sous-zones, et à rendre inaudible dans l'autre sous-zone. Par exemple, le champ sonore choisi est audible dans la première sous-zone SZ1. Le champ sonore choisi est à rendre inaudible dans la deuxième sous-zone SZ2. Les haut-parleurs peuvent être définis par leur position dans la zone. 5 [31] Selon un exemple de réalisation, chaque sous-zone SZ peut être définie par les données de géolocalisation d'un auditeur. Notamment, les données de géolocalisation d'un auditeur U sont obtenues au moyen d'un capteur de position POS. Il est alors possible de définir, en fonction des données de géolocalisation de l'auditeur, la première sous-zone SZ1, dans laquelle l'auditeur Ul entend le champ sonore 10 choisi. La sous-zone SZ1 a par exemple des dimensions prédéfinies. Notamment, la première sous-zone peut correspondre à un périmètre de quelques dizaines de centimètres à quelques dizaines de mètres, dont l'auditeur Ul est le centre. La deuxième sous-zone 5Z2, dans laquelle le champ sonore choisi est à rendre inaudible, peut être définie comme la sous-zone complémentaire de la première sous-zone SZ1. L'auditeur U2, présent dans la deuxième sous-zone 5Z2 n'entend donc pas le champ sonore choisi. 15 [32] Selon une variante de réalisation, la définition des sous-zones n'est pas axée sur la position des auditeurs U mais est prédéfinie et constante. Cela peut notamment être le cas dans des musées, où un champ sonore différent est diffusé en association avec chaque tableau. Un auditeur entend alors uniquement le champ sonore associé au tableau devant lequel il se trouve. Le périmètre autour du tableau représente alors la 20 première sous-zone, tandis que le reste du musée représente la deuxième sous-zone. [33] La zone est en outre couverte par des microphones MIC. Plus précisément, chaque sous-zone SZ comprend au moins un microphone MIC. Dans un exemple de réalisation, la zone est couverte par un réseau de M microphones MIC, avec M au moins égal à 2 et de préférence au moins égal à 10. Par exemple, la 25 première sous-zone SZ1 est couverte par les microphones MICi à MICmi et la deuxième sous-zone 5Z2 est couverte par les microphones MICmi+i à MICm. Dans une réalisation particulière, les microphones MIC sont des microphones virtuels. Dans la suite de la description, le terme « microphone MIC » est utilisé, les microphones pouvant être réels ou virtuels. Les microphones MIC sont repérés par leur position dans chacune des sous-zones SZ1, 5Z2. 30 [34] Dans un exemple de réalisation dans lequel le système comprend des microphones virtuels, des points sont définis dans chacune des sous-zones. Les points sont repérés par leurs positions respectives. Un microphone virtuel MIC est associé à chacun des points, de sorte que la position de chaque microphone virtuel est celle du point auquel il est associé. 8 [35] Dans un exemple de réalisation, les points sont définis une seule fois, la position des microphones virtuels MIC étant constante. Selon une autre réalisation, les points peuvent être redéfinis perpétuellement, de sorte que la position des microphones virtuels évolue. Il en est de même dans le cas où des microphones réels 5 sont utilisés. [36] Le système SYST comprend en outre une unité de traitement TRAIT apte à mettre en oeuvre les étapes du procédé. L'unité de traitement TRAIT comprend notamment une mémoire, formant un support de stockage d'un programme d'ordinateur comprenant des portions de code pour la mise en oeuvre du procédé décrit ci- 10 après en référence aux figures 2 et 3. L'unité de traitement TRAIT comprend en outre un processeur PROC apte à exécuter les portions de code du programme d'ordinateur. [37] L'unité de traitement TRAIT reçoit en entrée la position des microphones MIC, les données de géolocalisation des auditeurs U, les positions de chaque haut-parleur HP, les signaux audio à reproduire 15 S(Ui...Us) destinés respectivement aux utilisateurs U1 à Us et le champ sonore cible Pti, Pt2 à atteindre dans chaque sous-zone. A partir de ces données, l'unité de traitement TRAIT calcul le filtre FILT à appliquer aux signaux S(Ui...Us) afin de restituer les champs sonores cibles Pti, Pt2 dans les sous-zones. L'unité de traitement TRAIT délivre en sortie les signaux filtrés S(HPi...HPN) à diffuser respectivement sur les haut-parleurs HP1 à HPN 20 [38] Le système SYST trouve de nombreuses applications. Par exemple, le système SYST peut être destiné à des musées. Par exemple, des haut-parleurs HP sont disposés dans le musée et diffusent des signaux sonores correspondant notamment à des commentaires ou explications des oeuvres d'art. Devant chaque oeuvre d'art peut être définie une première sous-zone, le reste du musée formant la deuxième sous-zone dans laquelle les 25 commentaires correspondant à cette oeuvre d'art sont à rendre inaudible. Dans cette application, l'unité de traitement TRAIT peut être intégrée à un serveur du musée. [39] Une autre application peut être celle du visionnage de contenus multimédia, par exemple pour des particuliers visionnant un film. Ainsi, l'unité de traitement TRAIT peut être installée dans un boîtier de 30 téléviseur. Une première sous-zone dans laquelle les signaux sonores du film sont audibles peut être un périmètre autour du téléviseur, et la deuxième sous-zone dans laquelle les signaux sonores sont à rendre inaudibles peut être le reste de l'habitation. 9 [40] En variante, dans le cas où plusieurs auditeurs sont présents et souhaitent entendre des signaux sonores différents, par exemple dans différentes langues, plusieurs sous-zones peuvent être formées en fonction des positions de chacun des utilisateurs. Ainsi, une première sous-zone est formée autour de chaque utilisateur, le reste de la maison formant la deuxième sous-zone. Chaque utilisateur entend ainsi les signaux sonores dans la 5 langue choisie tandis que ces signaux sonores ne sont pas audibles dans le reste de l'habitation. [41] Les figures 2a et 2b illustrent les principales étapes du procédé selon l'invention. Les étapes du procédé sont mises en oeuvre par l'unité de traitement TRAIT de manière continue et en temps réel. 10 [042] Le procédé tel que décrit en référence aux figures 2a et 2b a pour but de rendre inaudible le champ sonore choisi dans l'une des sous-zones, par exemple dans la deuxième sous-zone SZ2. Le procédé est basé sur une estimation de pressions acoustiques dans chacune des sous-zones, de manière à appliquer un niveau de contraste sonore souhaité entre les deux sous-zones. A chaque itération, les signaux audio S(U1, U2) sont filtrés en fonction des pressions acoustiques estimées et du niveau de contraste sonore pour obtenir les signaux 15 de commande S(HPi...HPN) à diffuser sur les haut-parleurs. [043] Selon l'exemple de réalisation, non limitatif, illustré sur les figures 2a, 2b, les sous-zones SZ1, 5Z2 sont définies par la position d'au moins un auditeur Ul, U2. 20 [044] A l'étape SO, les données de géolocalisation d'au moins un des deux auditeurs Ul, U2 sont déterminées, par exemple au moyen d'un capteur de position POS. A partir de ces données de géolocalisation, les deux sous-zones SZ1, 5Z2 sont définies. Par exemple, la première sous-zone correspond aux données de géolocalisation du premier auditeur Ul. La première sous-zone SZ1 est par exemple définie comme étant une zone de quelques dizaines de centimètres à quelques dizaines de mètres de circonférence, dont le premier 25 auditeur Ul est le centre. La deuxième sous-zone 5Z2 peut être définie comme étant complémentaire de la première sous-zone SZ1. [045] A l'étape Si, le réseau de microphones MIC est défini. Par exemple, dans le cas de microphones virtuels, au moins un point, et de préférence une pluralité de points, chacun associé à une position dans l'une 30 des sous-zones SZ1, 5Z2 est défini. Notamment, chaque sous-zone SZ1, 5Z2 comprend au moins un point. De préférence, chaque sous-zone SZ1, 5Z2 comprend une dizaine de points. Les positions des points dans les sous-zones SZ1, 5Z2 définissent les positions des microphones MIC virtuels dans chaque sous-zone. 10 [46] A l'étape S2, la position de chaque haut-parleur HP est déterminée Notamment, la zone comprend un réseau de haut-parleurs comprenant au moins deux haut-parleurs HP. De préférence, le réseau de haut-parleurs comprend une dizaine de haut-parleur HP. Les haut-parleurs HP peuvent être répartis dans la zone de manière à ce que l'intégralité de la zone soit couverte par les haut-parleurs. 5 [47] A l'étape S3, une distance entre chaque couple de haut-parleur HP et de microphone MIC est calculée. Cela permet de pouvoir calculer chacune des fonctions de transfert Ftransf, pour chaque couple haut-parleur HP/microphone MIC, à l'étape S4. [Pti 10 [048] Autrement dit, le champ sonore cible peut être défini comme un vecteur Pt(to, n) = pour Pt21' l'ensembles des microphones MIC, à chaque instant n pour une pulsation co = 2n-f, f étant la fréquence. Les microphones MICi à MICm sont disposés aux positions xmic = MIC] et capturent un ensemble de pressions acoustiques regroupés dans le vecteur P(co, n). 15 [049] Le champ sonore est reproduit par les haut-parleurs (HP1,...,HPN) fixes et ayant comme position respective xiip = [H131, , HP] . Les haut-parleurs (HP1,...,HPN) sont pilotés par un ensemble de poids regroupés dans le vecteur q(co, n) = [qi (co, n), , qN (co, n)] T. L'exposant T est l'opérateur de transposition. 20 [050] Le trajet de propagation du champ sonore entre chaque couple de haut-parleur HP et microphone MIC peut être défini par un ensemble de fonctions de transferts G(co, n) assemblées dans la matrice [G11(60, n) - G11,(W, n) G(co, n) = I [Gmi (co, n) - - - GAIN(CO,n) [51] Avec les fonctions de transfert définies comme étant : 25 Gnit = jpck e- jkRmi, avec Ri,' la distance entre un couple haut-parleur et microphone, k le nombre 4n-Rmi d'onde, p la masse volumique de l'air et c la célérité du son. [52] A l'étape S5, les pressions acoustiques P1 et P2 sont déterminées respectivement dans la première sous-zone SZ1 et dans la deuxième sous-zone SZ2. 30 11 [53] Selon un exemple de réalisation, la pression acoustique P1 dans la première sous-zone SZ1 peut être la pression acoustique issue des signaux diffusés par les haut-parleurs dans la première sous-zone. La pression acoustique P2 dans la deuxième sous-zone, dans laquelle les signaux sonores sont à rendre inaudibles, peut correspondre à la pression acoustique induite issue des signaux diffusés par les haut-parleurs alimentés par les 5 signaux de commande associés à la pression P1 induite dans la première sous-zone. [54] Les pressions acoustiques P1, P2 sont déterminées à partir des fonctions de transfert Ftransf calculées à l'étape S4, et d'un poids initial appliqué aux signaux de commande de chaque haut-parleur. Le poids initial appliqué aux signaux de commande de chacun des haut-parleurs est égal à zéro. Puis, le poids appliqué aux 10 signaux de commande tend à varier à chaque itération, tel que décrit ci-après. [55] Selon cet exemple de réalisation, les pressions acoustiques P1, P2 comprennent chacun l'ensemble des pression acoustiques déterminées en chacune des positions des microphones virtuels. Ainsi, la pression acoustique estimée dans les sous-zones est plus représentative. Cela permet d'obtenir un résultat homogène en 15 sortie de procédé. [56] En variante, une pression acoustique déterminée en une seule position P1 , P2 est respectivement estimée pour la première sous-zone SZ1 et pour la deuxième sous-zone 5Z2. Cela permet de limiter le nombre de calculs, et donc de diminuer le temps de traitement et par conséquent la réactivité du système. 20 [57] Plus précisément, les pressions acoustiques P1 , P2 dans chacune des sous-zones peuvent être P rassemblées sous la forme d'un vecteur défini comme - p(co, n) = [ Il = G(co, n)q(co, n) P2 [58] A l'étape S6, les niveaux sonores L1 et L2 sont déterminés respectivement dans la première sous-zone 25 SZ1 et dans la deuxième sous-zone 5Z2. Les niveaux sonores L1 et L2 sont déterminés en chaque position des microphones MIC. Cette étape permet de convertir les valeurs des pressions acoustiques estimées P1, P2 en des valeurs mesurables en décibels. De cette manière, le contraste sonore entre la première et la deuxième sous-zone peut être calculé. A l'étape S7, un niveau de contraste sonore souhaité Cc entre la première sous-zone et la deuxième sous-zone est défini. Par exemple, le contraste sonore souhaité Cc entre la première sous- 30 zone SZ1 et la deuxième sous-zone SZ2 est préalablement défini par un concepteur en fonction du champ sonore choisi et/ou de la perception sonore des auditeurs Ul, U2 dans chaque sous-zone. [59] Plus précisément, le niveau sonore L pour un microphone peut être défini par 12 L = 20 10g10 (1.1), avec po la pression acoustique de référence, c'est-à-dire le seuil de perception. Po [60] Ainsi, le niveau sonore moyen dans une sous-zone peut être défini comme - ( pH pi, L = 10 10g10 avec PH la transposée conjuguée du vecteur de pressions acoustiques dans la 7) , , 5 sous-zone et M le nombre de microphones dans cette sous-zone. [61] A partir du niveau sonore L1, L2 dans les deux sous-zones, il est possible de calculer le contraste sonore estimé C entre les deux sous-zones : C = L1 - L2. 10 [062] A l'étape S8, la différence entre le contraste sonore estimé entre les deux sous-zones et le contraste sonore souhaité Cc est calculée. A partir de cette différence, un coefficient d'atténuation peut être calculé. Le coefficient d'atténuation est calculé et appliqué à la pression acoustique estimée P2 dans la deuxième sous-zone à l'étape S9. Plus précisément, un coefficient d'atténuation est calculé et appliqué à chacune des pressions acoustiques estimées P2 en chacune des positions des microphones MIC de la deuxième sous-zone 15 SZ2. La pression acoustique cible Pt2 dans la deuxième sous-zone prend alors la valeur de la pression acoustique P2 atténuée de la deuxième sous-zone. [63] Mathématiquement, la différence C. entre le contraste sonore estimé C et le contraste sonore souhaité Cc peut être calculée comme suit C. = C - Cc = L1 - L2 - Cc On peut alors calculer le coefficient c 20 d'atténuation = 10. [64] Ce coefficient est déterminé par l'amplitude de la pression acoustique à donner à chaque microphone pour que le niveau sonore dans la deuxième sous-zone soit homogène. Lorsque le contraste est équivalent à celui correspondant au contraste sonore souhaité Cc pour un microphone dans la deuxième sous-zone, alors 25 C. 5:---, 0 donc 5:---, 1. Cela signifie que la pression acoustique estimée en ce microphone correspond à la valeur de pression cible dans la deuxième sous-zone. [65] Lorsque la différence entre le contraste sonore estimé C et le contraste sonore souhaité Cc est négative C. 0, la pression acoustique en ce point est trop faible. Elle doit donc être augmentée pour correspondre au contraste sonore souhaité dans la deuxième sous-zone. 5 [067] Le principe est donc d'utiliser le champ de pression présent dans la deuxième sous-zone qui est induit par la pression acoustique dans la première sous-zone, puis d'atténuer ou d'amplifier les valeurs individuelles de pressions acoustiques estimées en chaque microphone, de sorte à ce qu'elles correspondent au champ sonore cible dans la deuxième sous-zone sur l'ensemble des microphones. Pour tous les microphones, on définit le vecteur : e, [el,...,e,,,,1T. 10 [068] Ce coefficient est calculé à chaque itération et peut donc évoluer. Il peut donc être écrit sous la forme e(n). 15 [069] En variante, dans le cas où une unique pression acoustique P2 est estimée pour la deuxième sous-zone SZ2, un seul coefficient d'atténuation est calculé et appliqué à la pression acoustique P2. [70] Les coefficients d'atténuation sont calculés de manière à répondre au critère de contraste défini par le concepteur. Autrement dit, le coefficient d'atténuation est défini de sorte que la différence entre le contraste 20 sonore entre les deux sous-zones SZ2 et le contraste sonore souhaité Cc soit proche de zéro. [71] Les étapes S10 à S12 permettent de définir la valeur des pressions acoustiques cibles Pti, Pt2 dans la première et la deuxième sous-zone SZ1, SZ2. 25 [072] L'étape S10 comprend l'initialisation des pressions acoustiques cibles Pti, Pt2, respectivement dans la première et la deuxième sous-zone SZ1, SZ2. Les pressions acoustiques cibles Pti, Pt2 caractérisent le champ sonore cible à diffuser dans les sous-zones. La pression acoustique cible Pti dans la première sous-zone SZ1 est définie comme étant une pression cible Pti, choisie par le concepteur. Plus précisément, la pression cible Ptidans la première sous-zone SZ1 est supérieure à zéro, de sorte que le champ sonore cible soit audible dans 30 cette première sous-zone. La pression acoustique cible Pt2 dans la seconde sous-zone est initialisée à zéro. Les pressions cibles Pti, Pt2 sont ensuite transmises à l'unité de traitement TRAIT à l'étape S11, sous la forme d'un vecteur Pt. [073] A chaque itération, on attribue de nouvelles valeurs de pressions cibles aux pressions cibles Pti, Pt2 déterminées à l'itération précédente. Cela correspond à l'étape S12. Plus précisément, la valeur de la pression 14 cible Ptidans la première sous-zone est celle définie à l'étape S10 par le concepteur. Le concepteur peut modifier cette valeur à tout moment. La pression acoustique cible Pt2 dans la deuxième sous-zone prend la valeur de la pression acoustique P2 atténuée (étape S9). Cela permet, à chaque itération, de redéfinir le champ sonore cible à restituer dans la deuxième sous-zone, en tenant compte de la perception de l'auditeur et des 5 signaux de commande des haut-parleurs. Ainsi, la pression acoustique cible Pt2 de la deuxième sous-zone n'est égale à zéro que lors de la première itération. En effet, dès lors que les haut-parleurs diffusent un signal, un champ sonore est perçu dans la première sous-zone, mais également dans la deuxième sous-zone [074] Mathématiquement, la pression cible Pt2 dans la deuxième sous-zone est calculée comme suit. 10 [75] A la première itération, Pt2 est égale à zéro : Pt2 (0) = O. [76] A chaque itération, la pression acoustique P2 estimée dans la deuxième sous-zone est calculée. Cette pression acoustique correspond à la pression acoustique induite dans la seconde sous-zone par le rayonnement 15 des haut-parleurs dans la première sous-zone. Ainsi, à chaque itération on a: P2 (co, n) = G2 (co, n)q(co, n), avec G2 (co, n) la matrice de fonctions de transfert dans la deuxième sous-zone à l'itération n. [77] La pression cible Pt2 à l'itération n + 1 peut donc être calculée comme Pt2 (n + 1) = .(n) X P2. 20 [078] A l'étape S13, l'erreur entre la pression cible Pt2 et la pression estimée P2 dans la deuxième sous-zone est calculée. L'erreur est due au fait qu'un pas d'adaptation y est appliqué de manière à ce que la pression cible Pt2 ne soit pas atteinte immédiatement. La pression cible Pt2 est atteinte au bout d'un certain nombre d'itérations du procédé. Cela permet de minimiser les ressources en calcul nécessaires pour atteindre la pression cible Pt2 dans la deuxième sous-zone SZ2. Cela permet en outre d'assurer la stabilité de 25 l'algorithme. De la même manière, le pas d'adaptation y est également choisi de sorte que l'erreur calculée à l'étape S13 ait une petite valeur, afin de stabiliser le filtre. [079] L'erreur calculée à l'étape S13 sert à calculer des poids à appliquer aux signaux de commande des haut-parleurs HP. Plus précisément, à la première itération les poids sont initialisés à zéro (étape S14). Chaque 30 haut-parleur diffuse un signal de commande non pondéré. Puis, à chaque itération, la valeur des poids varie en fonction de l'erreur (étape 515). Les haut-parleurs diffusent alors un signal de commande pondéré, qui peut être différent à chaque nouvelle itération. Cette modification des signaux de commande explique notamment que les pressions acoustiques Pti, Pt2 estimées dans les sous-zones SZ1, 5Z2 puissent être différentes à chaque itération. 15 [80] Mathématiquement, les poids q peuvent être calculés comme suit. q(n + 1) = q(n)(1 - py(n)) + AH(n)(G(n)q(n) - Pt(n)), avec y le pas d'adaptation pouvant varier à chaque itération et y(n) correspondant, dans un exemple de réalisation, à un 5 facteur d'oubli pouvant varier dans le temps et éventuellement être défini arbitrairement. Afin de garantir la stabilité du filtre, il est avantageux d'éviter que le pas d'adaptation y ne soit supérieur à l'inverse de la plus grande valeur propre de GHG. [81] Les filtres FILT à appliquer aux haut-parleurs sont alors déterminés à l'étape S16. Un filtre par haut- 10 parleur HP est par exemple calculé. Il peut donc y avoir autant de filtres que de haut-parleurs. Le type de filtres appliqués à chaque haut-parleur comprend par exemple une transformée de Fourier inverse. [82] Les filtres sont alors appliqués aux signaux audio à reproduire S(Ui, U2) qui ont été obtenus à l'étape S17. L'étape S17 est une étape d'initialisation, mise en oeuvre uniquement à la première itération du procédé. 15 Les signaux audio à reproduire S(Ui, U2) sont destinés respectivement à chacun des auditeurs U1, U2 en fonction de la sous-zone dans laquelle ils se trouvent. A l'étape S18, les filtres FILT sont appliqués aux signaux S(Ui, U2), en vue d'obtenir N signaux de commande S(HP1,...,HPN) filtrés à diffuser respectivement par les haut-parleurs HP1, HPN à l'étape S19. Les signaux de commande S(H131,...,HPN) sont diffusés respectivement par chaque haut-parleur (H131,...,HPN) du réseau de haut-parleurs à l'étape S20. De manière 20 générale, les haut-parleurs HP diffusent les signaux de commande en continu. [83] Puis, à chaque itération, les filtres FILT sont calculés en fonction des signaux S(H131,...,HPN) filtrés à l'itération précédente et diffusés par les haut-parleurs, tels que perçus par le réseau de microphones. Les filtres FILT sont appliqués aux signaux S(Ui, U2) pour obtenir de nouveaux signaux de commande S(H131,...,HPN) à 25 diffuser respectivement sur chaque haut-parleur du réseau de haut-parleurs. [84] Le procédé est alors relancé à partir de l'étape S5 dans laquelle les pressions acoustiques P1, P2 des deux sous-zones SZ1, 5Z2 sont estimées. 30 [85] Du fait des étapes décrites ci-dessus, le nombre d'itérations pour lesquelles le procédé converge, c'est-à-dire pour que la pression acoustique estimée P2 dans la deuxième sous-zone 5Z2 soit sensiblement égale à la pression acoustique cible Pt2 dans la deuxième sous-zone, est compris entre 10 et 1000, de préférence entre 200 et 800 et encore plus de préférence entre 400 et 500. Autrement dit, cela correspond au nombre 16 d'itérations pour que l'erreur soit égale à zéro ou sensiblement égale à zéro. Par « sensiblement égale », on entend qu'une marge d'erreur faible puisse exister. La marge d'erreur permet tout de même de rendre inaudible le champ sonore choisi dans la deuxième sous-zone. Lorsque l'erreur est égale, ou sensiblement égale, à zéro, les poids calculés et appliqués aux signaux de commande des haut-parleurs sont les mêmes que 5 ceux calculés et appliqués à l'itération précédente. [086] Le nombre d'itérations est ainsi relativement faible, ce qui permet un temps de convergence rapide. D'autre part, la mise en oeuvre d'un procédé itératif permet de limiter les ressources en calcul nécessaires. 10 [087] Bien entendu, d'autres exemples de réalisation selon l'invention sont possibles. Notamment, il a été décrit en référence aux figures 2a, 2b que les deux sous-zones SZ1, SZ2 sont définies une seule fois. Plus précisément, les données de géolocalisation d'un auditeur permettent de définir la première sous-zone SZ1, dans laquelle le champ sonore choisi est audible, et par la même de définir le réseau de microphones MIC de la première sous-zone. La deuxième sous-zone SZ2 est alors définie comme étant complémentaire de la 15 première sous-zone SZ1, quelles que soient les données de géolocalisation des autres auditeurs U. [88] Dans une variante de réalisation illustrée sur la figure 3, les données de géolocalisation de l'auditeur Ul situé dans la première sous-zone SZ1 sont obtenues répétitivement, de sorte que la première sous-zone SZ1, ainsi que le réseau de microphones de la première sous-zone SZ1, sont définis dynamiquement à chaque 20 déplacement de l'auditeur de la première sous-zone SZ1. De la même manière, la deuxième sous-zone SZ2 est redéfinie dynamiquement comme étant complémentaire de la première sous-zone SZ1, quelles que soient les données de géolocalisation des autres auditeurs U. [89] Sur la figure 3, la première sous-zone SZ1 est définie par la position al de l'auditeur Ul (représentés en 25 traits plein). Le réseau de microphones MIC est défini de manière à couvrir la première sous-zone SZ1. La deuxième sous-zone 5Z2 est complémentaire de la première sous-zone SZ1. La flèche F illustre un déplacement de l'auditeur Ul vers une position a2. La première sous-zone SZ1 est alors redéfinie autour de l'auditeur Ul (en traits pointillés). Le réseau de microphones MIC est redéfini de manière à couvrir la nouvelle première sous-zone SZ1. Le reste de la zone représente la nouvelle deuxième sous-zone 5Z2. Ainsi, 30 la première sous-zone SZ1 définie initialement par la position al de l'auditeur se trouve dans la deuxième sous-zone 5Z2. [090] Selon cette réalisation, le procédé est relancé, à chaque itération, depuis l'étape Si. Toutefois, la pression cible Pt2 souhaitée dans la deuxième sous-zone est estimée en fonction des poids appliqués aux 17 signaux de commande lors des itérations précédentes. Le procédé comprend donc une mémoire des poids calculés aux itérations précédentes, ce qui peut être problématique lorsque les sous-zones sont redéfinies dynamiquement. Une étape supplémentaire peut alors être prévue. Cette étape comprend la détermination d'un coefficient d'oubli, à appliquer aux poids à l'étape S15. Ce coefficient d'oubli permet d'atténuer le terme 5 correspondant au poids calculé à l'itération précédente. De cette manière, le procédé peut tourner en continu, tout en limitant l'erreur due à la redéfinition dynamique des sous-zones. D'autre part, le facteur d'oubli permet de régulariser le système de sorte à éviter une divergence du système. [91] Comme décrit en référence à l'étape 515 de la figure 2, les poids sont calculés comme suit : 10 q(n + 1) = q(n)(1 - ,uy(n)) + ,uGH(n)(G(n)q(n) - Pt(n)), où y est le facteur d'oubli. [92] Selon une autre réalisation, le procédé comprend une étape de calcul d'une dégradation entre la pression acoustique P1 estimée dans la première sous-zone SZ1 et la pression acoustique cible Pti souhaitée dans la première sous-zone SZ1. Le calcul des poids à l'étape 515 prend en compte le résultat du calcul de la 15 dégradation, de manière à ce que la pression P1 dans la première sous-zone SZ1 reste égale à la pression cible Pti souhaitée. La pression cible Pti peut rester constante durant toute la mise en oeuvre du procédé. En variante, la pression acoustique cible Pti peut être modifiée par le concepteur à tout moment à l'étape S11. Le procédé est alors relancé à partir de l'étape S5 d'estimation des pressions acoustiquesPti, Pt2 dans les sous-zones, en fonction des signaux de commande pondérés. 20 [93] Selon une autre variante de réalisation, le champ sonore choisi est un premier champ sonore, au moins un deuxième champ sonore choisi étant diffusé par le réseau de haut-parleurs HP. Le deuxième champ sonore choisi est audible dans la deuxième sous-zone et est à rendre inaudible dans la première sous-zone. Les haut-parleurs sont alimentés par les premiers signaux de commande pour émettre chacun un signal audio en continu 25 correspondant au premier champ sonore choisi, et sont également alimentés par des deuxièmes signaux de commande pour émettre chacun un signal audio en continu correspondant au deuxième champ sonore choisi. Les étapes du procédé telles que décrites ci-avant peuvent être appliquées à la première sous-zone SZ1, de sorte que le deuxième champ sonore choisi soit rendu inaudible dans la première sous-zone SZ1. 30 [094] Dans cet exemple de réalisation, la définition des première et deuxième sous-zones SZ1, 5Z2 peut être effectuée de manière indépendante pour chaque auditeur U présent dans la zone. Par exemple, une première sous-zone SZ1 est définie en fonction de la position d'un premier auditeur Ul, la deuxième sous-zone 5Z2 étant définie comme étant complémentaire de la première sous-zone SZ1, quelles que soient les données de géolocalisation de chaque autre auditeur U présent dans la zone. La définition des sous-zones est définie de la 18 même manière pour chacun des auditeurs U présents dans la zone. Ainsi, en se plaçant du point de vue d'un auditeur U à la fois, la zone comprend uniquement deux sous-zones SZ1, SZ2. De manière plus globale, la zone peut être décomposée en autant de sous-zones, dans chacune desquelles un champ sonore différent peut être audible, chaque sous-zone étant définie par la position d'un auditeur dans la zone. Le procédé est alors 5 mis en oeuvre en parallèle pour chacune des sous-zones, afin que le champ sonore choisi diffusé dans une sous-zone ne soit pas audible dans le reste de la zone. 1095] Selon un autre exemple de réalisation, les première et deuxième sous-zones ne sont pas complémentaires. Par exemple, dans une zone, REVENDICATIONS1. Procédé assisté par des moyens informatiques, pour une restitution sonore spatialisée à partir d'un réseau de haut-parleurs, en vue d'une diffusion d'un champ sonore choisi, audible sélectivement dans une première sous-zone (SZ1) prédéterminée d'une zone couverte par les haut-parleurs (HP1, HPN), ladite zone comportant en outre une deuxième sous-zone (SZ2) dans laquelle on cherche à rendre inaudible ledit champ sonore choisi, dans lequel, les haut-parleurs (H131,...,HPN) étant alimentés par des signaux de commande respectifs (S(HP1,...,HPN)) pour émettre chacun un signal audio en continu, le procédé comprend itérativement et en continu : - estimation d'une pression acoustique (P2) dans la deuxième sous-zone (SZ2), au moins en fonction des signaux de commande respectifs (S(HP1,...,HPN)) des haut-parleurs (HP1,...,HPN), et d'un poids initial respectif des signaux de commande (S(HP1,...,HPN)) des haut-parleurs (HP1,...,HPN), - calcul d'une erreur entre la pression acoustique estimée (P2) et une pression acoustique cible (Pt2), souhaitée dans la deuxième sous-zone (SZ2), - et, calcul et application de poids respectifs aux signaux de commande (S(HP1,...,HPN) ) des haut-parleurs, en fonction de ladite erreur, le calcul de la pression acoustique (P2) dans la deuxième sous-zone (SZ2) étant mis en oeuvre à nouveau en fonction des signaux de commande respectifs (S(HP1,...,HPN)), pondérés, des haut-parleurs (HP1,...,HPN). 2. Procédé selon la revendication 1, comportant préalablement une détermination de distances entre au moins un point de la deuxième sous-zone (SZ2) et des positions respectives de haut-parleurs (HP1,...,HPN),pour en déduire des fonctions de transfert acoustiques respectives des haut-parleurs (HP1,...,HPN) en ledit point, ledit point correspondant ainsi à une position de microphone virtuel (MIC), l'estimation de la pression acoustique (P2) dans la deuxième sous-zone (SZ2) étant alors effectuée en tenant compte d'un poids initial appliqué à chaque haut-parleur (H131,...,HPN),et de la fonction de transfert propre à chaque haut-parleur (HP1,...,HPN) en le microphone virtuel (MIC) dans la deuxième sous-zone (SZ2), et comportant en outre une détermination de distances entre au moins un point de la première sous-zone (SZ1) et des positions respectives de haut-parleurs (HP1,...,HPN), pour en déduire des fonctions de 20 transfert acoustiques respectives des haut-parleurs (HP1,...,HPN) en ledit point de la première sous-zone (SZ1), ledit point correspondant ainsi à une position d'un deuxième microphone virtuel (MIC). 5 3. Procédé selon la revendication 2, dans lequel on prévoit une pluralité de points au moins dans la deuxième sous-zone (SZ2) formant les positions respectives d'une pluralité de microphones virtuels (MIC), pour estimer une pluralité de pressions acoustiques (P2) respectives dans la deuxième sous-zone (SZ2) en tenant compte du poids respectif appliqué à chaque haut-parleur (HP1,...,HPN) et des fonctions de transfert propres à chaque haut-parleur (HP1,...,HPN) en chaque microphone virtuel (MIC). 10 4. Procédé selon l'une des revendications 2 ou 3, dans lequel la pression cible (Pt2) dans la deuxième sous-zone (SZ2) est évaluée en fonction d'un contraste sonore souhaité (Ce) entre la première sous-zone (SZ1) et la deuxième sous-zone (SZ2), et comportant une estimation de : - la pression acoustique (P1) dans la première sous-zone (SZ1) en tenant compte du poids respectif 15 appliqué à chaque haut-parleur (HP1,...,HPN) et de la fonction de transfert propre à chaque haut- parleur (HP1,...,HPN) en le au moins un deuxième microphone virtuel (MIC), - et de là, le contraste sonore (C) en fonction de la pression acoustique estimée (P1) dans la première sous-zone (SZ1) et la pression acoustique estimée dans la deuxième sous-zone (P2). 20 5. Procédé selon la revendication 4, dans lequel, à chaque itération, un coefficient d'atténuation, défini par la différence entre le contraste sonore souhaité (Ce) et le contraste sonore estimé (C) entre la première sous-zone (SZ1) et la deuxième sous-zone (SZ2), est appliqué à la pression acoustique (P2) dans la deuxième sous-zone (SZ2), la pression cible (Pt2) dans la deuxième sous-zone (SZ2) prenant alors la valeur de la pression acoustique 25 atténuée de la deuxième sous-zone. 6. Procédé selon l'une des revendications 2 à 5, comportant préalablement, pour chaque auditeur : - obtenir des données de géolocalisation d'un auditeur dans la zone au moyen d'un capteur de position, - définir la première sous-zone (SZ1) comme correspondant auxdites données de géolocalisation ainsi 30 que la position dudit au moins un deuxième microphone virtuel (MIC) dans ladite première sous-zone (SZ1), - définir la deuxième sous-zone (SZ2) comme étant complémentaire de la première sous-zone (SZ1) ainsi que la position dudit au moins un premier microphone virtuel (MIC) dans ladite deuxième sous-zone (SZ2). 21 7. Procédé selon l'une des revendications 2 à 6, comportant en outre, pour chaque auditeur : - obtenir répétitivement des données de géolocalisation de l'auditeur dans la zone au moyen d'un capteur de position, et en cas de changement des données de géolocalisation : 5 > définir dynamiquement la première sous-zone (SZ1) comme correspondant auxdites données de géolocalisation ainsi que la position dudit au moins un deuxième microphone virtuel (MIC) dans ladite première sous-zone (SZ1), > définir dynamiquement la deuxième sous-zone (SZ2) comme étant complémentaire de la première sous-zone (SZ1) ainsi que la position dudit au moins un premier microphone virtuel (MIC) dans 10 ladite deuxième sous-zone (SZ2), les distances entre ledit au moins un premier microphone virtuel (MIC) et les positions respectives des haut-parleurs (HP1,...,HPN) étant calculées dynamiquement pour en déduire les fonctions de transfert acoustiques respectives des haut-parleurs (HP1,...,HPN) en ledit au moins un premier microphone virtuel (MIC), et 15 l'estimation de la pression acoustique (P2) dans la deuxième sous-zone (SZ2) étant alors effectuée en tenant compte du poids respectif appliqué à chaque haut-parleur (HP1,...,HPN), chaque poids respectif appliqué à un haut-parleur (HP1,...,HPN) comprenant un premier terme correspondant à un poids respectif calculé à une itération précédente et un deuxième terme fonction de l'erreur, un coefficient d'oubli atténuant le terme correspondant au poids calculé à l'itération précédente 20 étant alors appliqué à chaque poids respectif. 8. Procédé selon l'une des revendications 2 à 7, comportant en outre, itérativement et en continu : - estimation de la pression acoustique (P1) dans la première sous-zone, en fonction desdits signaux de commande respectifs des haut-parleurs (S(HP1,...,HPN)) en le au moins un deuxième 25 microphone virtuel (MIC) et du poids respectif des signaux de commande des haut-parleurs, - calcul d'une dégradation entre la pression acoustique estimée (P1) dans la première sous-zone (SZ1) et une pression acoustique cible (Pti), souhaitée dans la première sous-zone, ladite pression acoustique cible (Pti) étant préalablement définie, - et, application de poids respectifs aux signaux de commande des haut-parleurs (S(HP1,...,HPN)), 30 en fonction de ladite dégradation, afin que la pression acoustique (P1) dans la première sous-zone atteigne la pression cible (Pti), le calcul de la pression acoustique (P1, P2) dans la première sous-zone (SZ1) et dans la deuxième sous-zone (SZ2) étant respectivement mis en oeuvre à nouveau en fonction des signaux de commande respectifs, pondérés, des haut-parleurs. 22 9. Procédé selon l'une des revendications précédentes, dans lequel le nombre d'itérations du procédé, pour que la pression acoustique de la première sous-zone atteigne la pression acoustique cible, souhaitée dans la première sous-zone, est compris entre 10 et 1000, de préférence entre 200 et 800, et de préférence entre 5 400 et 500. 10. Procédé selon l'une des revendications 2 à 9, dans lequel le champ sonore choisi est un premier champ sonore, un deuxième champ sonore choisi étant en outre diffusé, audible sélectivement dans la deuxième sous-zone (SZ2), et dans lequel on cherche à rendre le deuxième champ sonore inaudible dans la première 10 sous-zone (SZ1), les haut-parleurs (HP1,...,HPN) étant alimentés par des premiers signaux de commande respectifs (S(HP1,...,HPN)) pour émettre chacun un signal audio en continu correspondant au premier champ sonore choisi, et par des deuxièmes signaux de commande respectifs (S(HP1,...,HPN)) pour émettre chacun un signal audio en continu correspondant au deuxième champ sonore choisi, les étapes du procédé selon l'une des revendications 1 à 9 étant appliquées à la première sous-zone, de 15 manière à rendre inaudible le deuxième champ sonore dans la première sous-zone. 11. Procédé selon l'une des revendications précédentes, dans lequel chaque sous-zone (SZ1, SZ2) comprend au moins un microphone virtuel (MIC) et deux haut-parleurs (HP1,...,HPN), et de préférence chaque sous-zone (SZ1, SZ2) comprend au moins une dizaine de microphones virtuels (MIC) et au moins une dizaine 20 de haut-parleurs (HP1,...,HPN). 12. Système de restitution sonore spatialisée à partir d'un réseau de haut-parleurs (HP1,...,HPN),en vue d'une diffusion d'un champ sonore choisi, audible sélectivement dans une première sous-zone (SZ1) prédéterminée d'une zone couverte par les haut-parleurs (HP1,...,HPN), ladite zone comportant en outre 25 une deuxième sous-zone (5Z2) dans laquelle on cherche à rendre inaudible ledit champ sonore choisi, caractérisé en ce qu'il comprend une unité de traitement (TRAIT) adaptée pour le traitement et la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 11. 13. Support de stockage d'un programme d'ordinateur, chargeable dans une mémoire associée à un 30 processeur (PROC), et comprenant des portions de code pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à 11 lors de l'exécution dudit programme par le processeur.