Brevet FR2851877A1 — publié le 20040903

Procédé de mesure de fonctions de transfert acoustiques associées à la morphologie d'un individu L'invention concerne une mesure de fonctions de transfert 5 acoustiques associées à la morphologie d'un individu, pour une synthèse binaurale ou transaurale de son spatialisé. Un individu, lorsqu'il est plongé, en tant qu'auditeur, dans un champ sonore, perçoit en trois dimensions les 10 différentes sources d'une scène sonore avec leurs informations de position, de niveau, d'étendue et de déplacement. La localisation auditive d'une source sonore repose sur 15 les différences entre les signaux reçus aux entrées des conduits auditifs gauche et droit de l'auditeur. Ces différences temporelles, d'intensité et spectrales sont interprétées par le cerveau de l'auditeur, qui en déduit, en corrélation avec des informations visuelles et par un 20 phénomène d'apprentissage depuis la naissance, la position de la source dans l'espace. Pour simuler la perception spatiale d'une source sonore, il faut reconstruire le signal dit "binaural" qui 25 atteindrait les entrées des conduits auditifs de l'auditeur si cet auditeur était réellement en présence de cette source. Cette construction d'un signal binaural synthétique est désignée par les termes "synthèse binaurale". Toutefois, le trajet d'une onde acoustique issue d'une telle source sonore jusqu'aux conduits auditifs de l'auditeur peut varier suivant la morphologie de l'auditeur, notamment suivant les dimensions de sa tête, de la distance entre ses deux conduits auditifs, ou 5 autres. On comprendra ainsi que la perturbation qu'induit un auditeur sur la propagation d'une onde acoustique est propre à la morphologie de l'auditeur et que cette perturbation permet à l'auditeur de localiser la provenance d'un son dans l'espace. Un individu possède 10 donc, pour chacune de ses oreilles, une fonction de transfert définissant cette perturbation, dite IIHRTFfV (pour "Heaci Related Transfer Function") et qui est propre à cet individu. Les fonctions de transfert HRTFs, pour l'oreille droite et l'oreille gauche, sont habituellement mesurées, sur un individu et pour cet individu, en faisant déplacer angulairement (suivant un angle azimutal et un angle d'élévation) un haut-parleur autour du centre de la tête 20 de l'individu, en disposant un microphone en chaque oreille de l'individu. Les microphones mesurent donc les perturbations des ondes sonores qu'émet le haut-parleur, pour différentes directions de propagation, et, de là, on estime les HRTFs de cet individu. Bien entendu, il peut 25 être prévu, en variante, de placer une source acoustique au voisinage de l'oreille de l'individu, ainsi que des points récepteurs (par exemple des microphones) aux positions auxquelles on souhaite mesurer les fonctions de transfert. Cette variante découle de l'application du 30 principe de réciprocité de la propagation des ondes acoustiques. Ce procédé de mesure des HRTFs, pour un individu et dit à cet effet "d'individualisation des HRTFs", est fastidieux et dure quelques heures. L'individu doit rester immobile 5 pendant cette durée et supporter un casque équipé de transducteurs électroacoustiques pour chaque oreille. On connaît notamment par le document: "Boundary element method calculation of individual head10 related transfer function. I. Rigid model calculation", B.F.G. Katz, J. Acoust. Soc. Am. 110(5), 2440-2448, Nov.2001, un procédé d'obtention des HRTFs d'un individu par une simulation, en modélisant la tête de l'individu et en 15 calculant les perturbations qu'un tel modèle générerait sur une onde acoustique. Toutefois, la modélisation d'une tête humaine, avec ses reliefs, mène à des calculs très longs pour obtenir les 20 HRTFs. Par ailleurs, de tels procédés de modélisation ne sont fiables que pour une bande de fréquences restreinte, typiquement entre 1000 et 5000 Hz environ, car des parties du corps de l'individu, autres que sa tête, sont aussi susceptibles de perturber la propagation d'une onde 25 acoustique et de telles perturbations ne sont pas prises en compte dans ce document. La présente invention vient améliorer la situation. L'un des buts de la présente invention est de fournir un procédé de mesure des HRTFs sur une gamme de fréquences étendue, préférentiellement d'environ 20 Hz à au moins kHz. Un autre but de la présente invention est de fournir un 5 procédé de mesure des HRTFs, de mise en oeuvre simple et rapide. Elle propose à cet effet un procédé de mesure de fonctions de transfert associées à la morphologie d'un individu, 10 pour une spatialisation du son, dans lequel: a) on mesure au moins les dimensions de la tête et du torse d'un individu, b) on définit, en fonction des mesures de l'étape a), un premier ellipsode de dimensions correspondant sensiblement aux dimensions de la tête de l'individu et un second ellipsode de dimensions correspondant sensiblement aux dimensions du torse de l'individu, et c) on estime: - une première perturbation d'une onde acoustique par au 20 moins lesdits premier et second ellipsodes, pour au moins une fréquence de l'onde acoustique comprise dans une gamme de basses fréquences, pour obtenir un premier jeu de fonctions de transfert dans les basses fréquences acoustiques, - et une seconde perturbation d'une onde acoustique par au moins ledit premier ellipsode, pour au moins une fréquence de l'onde acoustique comprise dans une gamme de hautes fréquences, pour obtenir un second jeu de fonctions de transfert dans les hautes fréquences acoustiques. La Demanderesse a en effet constaté que, pour une gamme de basses fréquences typiquement entre environ 20 Hz et environ 2000 Hz, le torse de l'individu joue un rôle considérable dans l'estimation des HRTFs et que, pour une 5 gamme de hautes fréquences typiquement entre environ 2 kHz et au moins 10 kHz, un simple ellipsode peut représenter fidèlement la tête de l'individu, sans nécessairement prendre en compte son relief précis. Ainsi, selon l'un des avantages que procure la présente invention, deux ellipsodes seulement (représentant la tête et le torse) modélisent correctement les parties du corps d'un individu pour obtenir les HRTFs de cet individu, et ce, pour une gamme de fréquences étendue 15 couvrant pratiquement la totalité du spectre des fréquences audibles. On comprendra ainsi que des calculs par simulation des perturbations d'une onde acoustique pour estimer les HRTFs sont beaucoup plus simples et plus rapides que ceux prévus dans les procédés de l'état de la 20 technique. On comprendra aussi que les HRTFs peuvent être obtenus directement par mesures des perturbations acoustiques que génère un mannequin comportant simplement deux ellipsodes représentant la tête et le torse d'un individu, les axes de ces ellipsodes étant réglés aux 25 dimensions de l'individu, ce qui évite à cet individu de subir le procédé de mesure classique "d'individualisation des HRTFs", décrit ci-avant. L'estimation des HRTFs à l'étape c) peut être effectuée ainsi par simulation ou par mesure directe sur un mannequin. Préférentiellement, à l'étape a), on repère au moins une position d'une oreille de l'individu et à l'étape b), on définit cette position sur le premier ellipsode, de sorte qu'à l'étape c), on estime les fonctions de transfert de l'individu pour cette position de l'oreille. Avantageusement, à l'étape a) précitée, on mesure en outre les dimensions du cou de l'individu et, à l'étape b), on définit en outre un cylindre de dimensions correspondant 10 sensiblement aux dimensions du cou de l'individu, de sorte qu'à l'étape c), on estime ladite première perturbation dans la gamme des basses fréquences, à la fois par ce cylindre et par les premier et second ellipsodes précités. Avantageusement, à l'étape a), on mesure en outre les dimensions d'au moins une oreille de l'individu, ainsi que la position de cette oreille sur la tête de l'individu. A l'étape b) : - on définit au moins une portion de cône, de dimensions correspondant sensiblement aux dimensions de l'oreille, et - on reporte cette portion de cône sur le premier ellipsode, de sorte qu'à l'étape c), on estime ladite seconde perturbation dans la gamme des hautes fréquences, à la fois par au moins cette portion de cône et par le premier ellipsode. On comprendra que le terme "portion de cône" désigne de façon générale aussi bien une portion de cylindre, ou autre. On peut obtenir ainsi les HRTFs de l'individu jusqu'à une fréquence d'environ 13 kHz. On indique qu'il est possible, bien entendu, de choisir les bornes des gammes de hautes fréquences et de basses 10 fréquences, selon l'application visée. Ainsi, en variante d'une gamme de basses fréquences entre 20 Hz et 2 kHz, il est possible de viser la gamme de 20 Hz à 8 kHz, en modélisant le torse, le cou et la tête de l'individu comme indiqué ciavant. Pour des hautes fréquences entre 8 kHz 15 et 13 kHz, on modélise dans cette variante la tête et une ou deux oreilles de l'individu. On indique en outre qu'il suffit en principe de modéliser une seule oreille de l'individu, compte tenu de la 20 symétrie par rapport au plan médian entre les deux oreilles. Dans un premier mode de réalisation, on obtient les HRTFs par une modélisation du buste de l'individu, à partir d'au 25 moins deux ellipsodes comme indiqué ci-avant, par simulation des perturbations. Dans ce mode de réalisation, à l'étape c) cl) on affecte une impédance acoustique choisie au moins 30 auxdits premier et second ellipsodes, et c2) on calcule donc par simulation lesdites première et seconde perturbations. L'impédance acoustique choisie à l'étape cl) est, de 5 préférence, infinie. En variante, elle peut être choisie proche de l'impédance acoustique du corps humain, mais les tests effectués avec une impédance infinie fournissent déjà des résultats satisfaisants. Préférentiellement, à l'étape c2) c21) on applique au moins un maillage auxdits premier et second ellipsodes, et c22) on calcule par éléments finis lesdites première et seconde perturbations. Pour une estimation plus fine des fonctions de transfert, ou encore pour réduire le temps de calcul de la simulation par éléments finis dans les basses fréquences, à l'étape c21), on applique avantageusement: - un premier maillage pour le calcul de la première perturbation dans les basses fréquences, et - un second maillage pour le calcul de la seconde perturbation dans les hautes fréquences, les dimensions d'une maille du second maillage étant 25 choisies plus petites que les dimensions d'une maille du premier maillage. Dans ce premier mode de réalisation o l'on estime les fonctions de transfert HRTFs par simulation des 30 perturbations, on place une source acoustique "virtuelle" au voisinage de la position précitée de l'oreille, de préférence à une distance d'environ 2 millimètres de cette position, au-dessus du premier ellipsode. On applique ensuite, pour le calcul des HRTFs, le principe précité de réciprocité, en calculant donc les perturbations, dans 5 l'espace autour du modèle, que subissent les ondes acoustiques issues de la source virtuelle. On indique que l'application du principe de réciprocité pour le calcul par simulation des HRTFs a un grand intérêt puisque l'on évite de réitérer un calcul numérique pour chaque position 10 de récepteur. Dans un second mode de réalisation, on obtient les HRTFs par une mesure des première et seconde perturbations que génère un mannequin représentant le buste d'un individu, 15 comme indiqué ci-avant. Dans ce second mode de réalisation, à l'étape c) cl) on prévoit donc un mannequin comportant un premier ellipsode représentant la tête d'un individu et un second 20 ellipsode, solidaire du premier ellipsode et représentant le torse d'un individu, ainsi que des moyens de réglage des axes desdits ellipsodes, c2) on place au moins un transducteur électroacoustique à une position choisie par rapport au mannequin, et c3) on mesure lesdites première et seconde perturbations d'une onde acoustique que produit ladite source sonore. A ce titre, la présente invention vise aussi un mannequin pour la mise en òuvre du procédé selon ce second mode de 30 réalisation, le mannequin comportant au moins: - un premier ellipsode représentant la tête d'un individu, - un second ellipsode représentant le torse d'un individu, des moyens de solidarisation du premier ellipsode au second ellipsode, des moyens de fixation d'au moins un transducteur électroacoustique en une position choisie sur ledit premier ellipsode, et - des moyens de réglage des axes respectifs desdits premier et second ellipsodes. Préférentiellement, la position choisie correspond sensiblement à la position d'une oreille sur la tête de 15 l'individu et le transducteur électroacoustique peut être aussi bien un microphone qu'un haut-parleur, comme indiqué ci-avant. Toutefois, dans une réalisation préférée, on place un 20 microphone, en tant que transducteur, avantageusement à environ 2 mm de la position de l'oreille, tandis qu'un haut-parleur, en tant que source acoustique réelle, balaye l'espace autour du mannequin. Avantageusement, le mannequin peut comporter d'autres parties de tailles réglables et représentant le corps de l'individu, notamment le cou et au moins une oreille de l'individu. Avantageusement, ces parties du mannequin telles qu'au moins les premier et second ellipsodes représentant la il tête et le torse peuvent être revêtues d'une mousse de polymère telle que du latex pour imiter la texture de la peau humaine. Dans un mode de réalisation préféré de l'étape a) du procédé au sens de la présente invention, al) on obtient une première prise de vue de l'individu, de profil, avec un étalon de mesure, et a2) on obtient une seconde prise de vue de l'individu, de 10 face, avec ledit étalon de mesure. Préférentiellement, pour la vue de profil, all) on repère sur la première prise de vue les positions respectives d'un sourcil et d'une oreille de l'individu, a12) on trace une droite passant par l'oreille et le sourcil, ainsi que par un premier point situé à l'arrière de la tête, a13) on détermine une droite médiane et perpendiculaire du segment défini par les positions du sourcil et dudit 20 premier point, et a14) on repère deux intersections de ladite droite médiane perpendiculaire avec un contour du visage de l'individu sur ladite première prise de vue, lesdites deux intersections correspondant sensiblement au sommet du 25 crâne de l'individu et au menton de l'individu, de sorte qu'à l'étape b), on définit un premier axe du premier ellipsode entre lesdites deux intersections et un second axe du premier ellipsode perpendiculaire et médian dudit premier axe sur ladite première prise de vue. 30 Préférentiellement, pour la vue de face, a21) on repère sur la seconde prise de vue les positions respectives des oreilles de l'individu, a22) on détermine une droite médiane et perpendiculaire du segment défini par lesdites positions des oreilles, et a23) on repère deux intersections de ladite droite médiane avec un contour du visage de l'individu sur ladite seconde prise de vue, lesdites deux intersections correspondant sensiblement au sommet du crâne de l'individu et au menton de l'individu, de sorte qu'à l'étape b), on définit un troisième axe du premier ellipsode entre lesdites deux intersections et un quatrième axe du premier ellipsode perpendiculaire et médian dudit troisième axe sur ladite seconde prise de vue. Comme on utilise le même étalon de mesure entre les vues de face et de profil, on fait avantageusement concider les premier et troisième axes précités et le premier ellipsode est finalement défini par lesdits premier, 20 second et quatrième axes. Pour le torse, préférentiellement a31) on repère sur la seconde prise de vue une position haute du sommet des épaules et une position basse des 25 coudes de l'individu, d'une part, et les limites du bras gauche et du bras droit, d'autre part, a32) on repère sur la première prise de vue les positions respectives de l'arrière du dos et de l'avant de la poitrine de l'individu, a33) on définit un volume délimité par les positions repérées aux étapes a31) et a32), et a34) on ajuste les axes du second ellipsode, d'une part, en orientant dans l'espace ledit ellipsode et, d'autre part, en recherchant un recouvrement volumique maximal avec le volume défini à l'étape a33). Avantageusement, on obtient alors un second ellipsode respectant l'orientation du torse, lorsque l'individu se tient dans une position naturelle. Par exemple, le second ellipsode peut suivre le galbe général du dos et/ou 10 d'autres formes particulières. D'autres caractéristiques et avantages de l'invention apparaîtront à l'examen de la description détaillée ciaprès, et des dessins annexés sur lesquels: - les figures lA et lB sont respectivement des vues de face et de profil de la tête d'un individu; - les figures 2A et 2B sont respectivement des vues de face et de profil du buste d'un individu, en particulier du torse; - les figures 3A et 3B sont respectivement des vues de face et de profil du buste d'un individu, en particulier du cou; - les figures 4A et 4B sont respectivement des vues de face et de profil de la tête d'un individu, en particulier 25 d'une oreille de l'individu; - les figures 5A à 5D représentent une modélisation du buste de l'individu, respectivement vue de dessus, de profil, de face et en perspective; - les figures 6A à 6D représentent une modélisation de la 30 tête et d'une oreille de l'individu, vue respectivement de dessus, de profil, de derrière et en perspective; et - la figure 7 représente schématiquement un mannequin, au sens de la présente invention, modélisant le buste d'un individu. On se réfère tout d'abord aux figures lA à 4B pour décrire, simplement à partir de deux photographies d'un individu (l'une prise de face et l'autre prise de profil), un procédé pour extraire de ces prises de vues des caractéristiques géométriques de "primitives" définissant 10 la morphologie de l'individu. On indique qu'un recouvrement volumique maximal entre ces primitives et les formes reconnues sur les photographies est recherché. On entend par "primitives" les formes de base, telles que des ellipsodes ou des cylindres, qui assurent une 15 modélisation appropriée de la tête, du torse, du cou, ou d'autres parties du corps d'un individu. Dans une première étape, on repère dans des représentations bidimensionnelles (les photographies des figures lA à 4B) les caractéristiques des primitives (dimensions d'axe 20 d'ellipsodes, de cylindre, ou autres) et, dans une seconde étape, on génère, à partir de ces caractéristiques de primitives, des modèles tridimensionnels. On indique que l'individu est coiffé d'un bonnet BON 25 (figure 2A) de manière à ce que ses cheveux ne perturbent pas les mesures anthropométriques effectuées sur la base desdites prises de vue. Pour ce qui concerne l'obtention des caractéristiques 30 géométriques de l'ellipsode représentant la tête de l'individu, en se référant aux figures lA et 1B, cet ellipsode doit être contenu au maximum dans le volume délimité par le segment [AB] (largeur de la tête), le segment [CD] (hauteur de la tête) et le segment [EF] (longueur de la tête). Dans le plan (OXY) de la figure lA, la projection de l'ellipsode, correspondant à une ellipse, est tangente à l'entrée des oreilles A et B, ainsi qu'au haut C et au bas D de la tête. L'axe interaural défini par la droite (CD) 10 partage l'ellipse en deux parties sensiblement égales. Dans le plan (OYZ) de la figure 1B, le point G représente la position du canal auditif et le point H représente la position du haut du sourcil. La droite (GH) passant par 15 ces deux positions coupe le contour arrière de la tête pour former le point I. La médiatrice du segment [IH] permet de placer la droite qui va contenir le grand axe de l'ellipse dans le plan de la figure 1B. Cette droite coupe le contour de la tête pour former le segment [JKI qui 20 correspond alors au grand axe de l'ellipse. La médiatrice du segment [JK] permet de positionner la droite qui contient le petit axe de l'ellipse et qui est parallèle à la droite (IH) . L'intersection de cette droite avec les contours de la tête forme un segment [LM] qui correspond 25 au petit axe de l'ellipse. On notera sur les figures lA et lB la présence d'un étalon de mesure ETA qui permet d'établir une corrélation entre les deux prises de vues, de face et de profil, et pour 30 déterminer ainsi le troisième axe de l'ellipsode, en comparant par exemple la dimension du grand axe de l'ellipse de la figure lB, avec la dimension du grand axe de l'ellipse de la figure lA. En se référant maintenant aux figures 2A et 2B, pour 5 définir les caractéristiques de l'ellipsode du torse, on définit les points suivants: - M qui correspond à la position du haut des épaules; - N qui correspond à la position des coudes; - Q qui correspond à l'avant de la poitrine; 10 P qui correspond à l'arrière du dos; - K qui correspond à la limite du bras droit; et - L qui correspond à la limite du bras gauche. L'ellipsode représentant le torse doit être contenu dans 15 le volume délimité par les segments [MN], [KL] et [PQI. Les points ainsi définis permettent de déterminer les dimensions de l'ellipsode, ainsi que son orientation et son positionnement dans l'espace, comme indiqué ci-avant en référence aux figures lA et 1B pour ce qui concerne 20 l'ellipsode représentant la tête de l'individu. Pour ce qui concerne la primitive représentant le cou d'un individu, on indique qu'une forme cylindrique est préférentiellement choisie. En se référant aux figures 3A 25 et 3B, les contours du cou forment deux droites parallèles (a) et (b) . Les points A, B, C et D de la figure 3B forment respectivement: - le haut de l'arrière du cou, - le haut de l'avant du cou, - le bas de l'arrière du cou, et - le bas de l'avant du cou. Les segments [AC] et [BD] sont' parallèles et préférentiellement de même longueur. Préférentiellement encore, l'ellipse correspondant à la base du cylindre 5 ainsi défini comporte un premier axe [AB] dans la direction (OZ) et un second axe d dans la direction [OX]. Pour ce qui concerne le pavillon auditif (ou oreille), on repère d'abord sa position sur la tête de l'individu. En 10 se référant aux figures 4A et 4D, la primitive représentant le pavillon auditif est définie par exemple par une portion de cône (ou une portion de cylindre, dans un cas particulier simplifié). Les points R, S, T et U représentent respectivement le haut de l'oreille, le bas 15 de l'oreille, l'arrière de l'oreille et l'avant de l'oreille. Ces points permettent de définir les dimensions de la primitive du pavillon auditif. Ainsi, dans un mode de réalisation préférée, on choisit de 20 s'appuyer sur des photographies pour obtenir les données anthropométriques de l'individu. On indique toutefois qu'en variante, les mesures de grandeurs anthropométriques peuvent être effectuées à l'aide d'un capteur de positions monté sur la tête de l'individu, tel que le capteur " POLH MUS" D. Pour obtenir les primitives tridimensionnelles représentant ensemble le buste de l'individu, on utilise préférentiellement un logiciel de modélisation tridimensionnelle, par exemple 3DS Max , qui permet de créer facilement des primitives géométriques de tous types. On détermine alors les caractéristiques géométriques des primitives précitées, tant en dimensions qu'en positions dans l'espace. Le modèle tridimensionnel peut être ainsi obtenu en se basant uniquement sur quelques mesures anthropométriques de l'individu. On entend par les termes "prise de vue", indiqués ciavant, aussi bien une photographie qu'une succession d'images, telle qu'une séquence animée. Ainsi, par une 10 mise en correspondance visuelle, on peut vérifier par rotation du modèle tridimensionnel son adéquation avec la morphologie de l'individu filmé par exemple avec une rotation de la caméra. En variante, la génération du modèle tridimensionnel du buste de l'individu peut être automatisée. On indique en effet que l'étape d'extraction des paramètres anthropométriques et des caractéristiques géométriques des primitives modélisant le buste de l'individu peut se faire 20 en utilisant un traitement d'image, par exemple par détection de contours. On projette alors dans un modéliseur 3D un plan de face et un plan de profil représentant les photographies de l'individu, prises de face et de profil. Pour chaque primitive représentant la tête, le cou, le torse et le pavillon auditif d'un l'individu, on indique une généralisation ci-après de la modélisation tridimensionnelle. Pour la tête: - on centre une primitive ellipsode sur la position relevée sur les photographies ci-avant, du centre de l'ellipsode; - on écrase la primitive sur la vue de profil selon l'axe 5 Y, de façon à atteindre la dimension du petit axe de l'ellipsode en vue de profil; - on écrase la primitive sur la vue de profil selon l'axe Z, de façon à atteindre la dimension du grand axe de l'ellipsode en vue de profil; - on effectue une rotation de l'ellipsode selon l'axe X, cette rotation étant d'un angle qui est calculé à partir des données de position et d'orientation du petit axe et du grand axe de l'ellipsode vu de profil comme indiqué ci-avant; et - on écrase l'ellipsode sur la vue de face selon l'axe X, de façon à atteindre la dimension du segment [AB], comme représenté sur la figure 1A. Pour le cou - on créé une primitive cylindre centrée sur un point choisi de la figure 3A; - on procède à un écrasement de la primitive cylindre selon l'axe X, de façon à atteindre la taille du segment d; - on écrase la primitive cylindre selon l'axe Z, de façon à atteindre la taille du segment [AB] ; et - on effectue une rotation du cylindre selon l'axe X, d'un angle qui est calculé à partir des données de position et d'orientation des points A et B. On indique que la position et l'orientation dans l'espace de la primitive représentant le torse s'appuie sur le même procédé ci-avant de modélisation de la tête par une primitive ellipsode. Par ailleurs, la position et 5 l'orientation dans l'espace de la portion de cône représentant le pavillon auditif s'appuient sensiblement sur le même procédé d'obtention de la primitive cylindre représentant le cou de l'individu. En outre, la position du pavillon auditif de l'individu est repérée sur les photographies de face et de profil, pour être reportée sur l'ellipsode représentant la tête de l'individu. On remarquera sur les prises de vues des figures lA à 4B la présence d'un même étalon de mesures ETA pour corréler les dimensions anthropométriques mesurées sur les prises de vues de face et de profil, et d'une partie du corps de l'individu représentée à l'autre. 20 On indique. que la Demanderesse a conçu un programme informatique pour automatiser l'extraction des paramètres anthropométriques à partir des photographies de face et de profil d'un individu, par reconnaissance de forme, ainsi 25 qu'une automatisation de la création des primitives tridimensionnelles, par positionnement et orientation dans l'espace des caractéristiques anthropométriques relevées. A ce titre, la présente invention vise aussi un tel programme d'ordinateur. On indique que ce programme peut 30 consister, en variante, en une implémentation de routines additionnelles ("plugins") d'un programme déjà existant, comme 3DS Max. En se référant maintenant aux figures 5A à 5D, les 5 primitives simples (ellipsodes, cylindres, portions de cônes) obtenues précédemment sont maillées pour un calcul acoustique par éléments finis de frontière(méthode dite ".BEMR," pour "Boundary Element Method") . Une fois les primitives représentant le modèle géométrique tridimensionnel créées, on génère différents maillages selon la gamme de fréquences sur laquelle on souhaite calculer les réponses acoustiques. On fonction de la fréquence maximale de calcul pour les HRTFs (valeur haute des gammes de fréquences basses et hautes, précitées), on 15 vérifie que la taille des éléments du maillage (éléments dits "mailles") est préférentiellement inférieure à 343/(6xfmax), o la valeur 343 correspond à la vitesse du son exprimée en mètres par seconde et fmax correspond à la fréquence maximale de calcul (2 kHz pour la gamme des 20 basses fréquences précitée et 13 kHz pour la gamme des hautes fréquences précitée). Ce choix correspond à une contrainte du calcul BEM, selon laquelle on utilise préférentiellement au moins six mailles par longueur d'onde. Comme indiqué ci-avant, il est préférable de 25 choisir un maillage de basse résolution (faible nombre de mailles) pour les basses fréquences. En revanche, un maillage de haute résolution est choisi pour des fréquences élevées (préférentiellement supérieures à 2 kHz). Préférentiellement, on utilise le principe de réciprocité, c'est-à-dire que l'on place une source acoustique virtuelle au niveau de l'oreille du modèle. Préférentiellement, cette source est placée à environ 2 mm 5 de la position de l'entrée du conduit auditif. Préférentiellement, un raffinement local du maillage est choisi, sur une zone de 2 à 3 cm autour de la position du conduit auditif. La taille maximale des éléments du maillage, sur cette zone, est ainsi inférieure à la 10 distance de 2 mm indiqué ci-avant. On applique ensuite les maillages grâce à un logiciel de calcul par éléments finis de frontière, tel que le logiciel VnoiseO. On importe ensuite les positions (sous 15 forme de fichier informatique) pour lesquelles on souhaite calculer les HRTFs. Ces positions correspondent à des emplacements de microphones virtuels disposés sur une sphère autour de l'individu modélisé. On place la source sonore virtuelle près du conduit auditif, modélisé et 20 maillé comme indiqué ci-avant. On utilise ensuite un algorithme de calcul tel que BICG-STAB(L) (pour "Bi Conjugate Gradient Stabilized, order L"I), couramment utilisé pour ce type d'application. Préférentiellement, on fixe l'ordre L=32, pour 2000 itérations. 25 Le maillage choisi est de-type triangulaire et les données du maillage (position des noeuds du maillage) sont traités par un programme de calcul basé sur la technique des éléments finis de frontière, avantageusement sur une bande 30 de fréquences logarithmiques (espacement en sixièmes d'octave), pour obtenir une résolution proche de la résolution auditive. Les résultats sont ensuite interpolés en amplitude sur une échelle linéaire, par une méthode dite de "spline". On indique simplement que le calcul des réponses fréquentielles sur une échelle linéaire, à partir 5 de la simulation des perturbations acoustiques sur une échelle logarithmique, est mené par interpolation, pour obtenir finalement les réponse impulsionnelles binaurales propres à l'individu modélisé. Le calcul des fonctions de transfert sur une échelle logarithmique permet de réduire efficacement le temps de calcul des HRTFs par modélisation. En effet, ce choix d'une échelle logarithmique permet de calculer moins de fréquences élevées qu'un calcul avec une échelle linéaire. 15 Or, les fréquences élevées sont ici les plus coteuses en temps de calcul. Comme on choisit ainsi deux maillages différents pour deux gammes de fréquences différentes, il est donc possible de choisir aussi différentes parties du buste de l'individu, modélisées, pour calculer les 20 fonctions de transfert binaurales HRTF. Préférentiellement, on choisit pour les basses fréquences au moins la tête et le torse de l'individu et, pour les hautes fréquences, au moins la tête de l'individu. Dans une réalisation préférée, on considère comme primitive, dans la gamme de fréquences 20 Hz à 2000 Hz, la primitive ellipsodale représentant la tête de l'individu, la primitive ellipsodale représentant son torse, et le cylindre représentant son cou. Dans la gamme de fréquences 30 de 2000 Hz à 13000 Hz, on considère simplement la primitive ellipsodale représentant la tête, ainsi qu'une portion de cône représentant le pavillon auditif de l'individu. Ainsi, en se référant aux figures 5A à 5D, le maillage précité est effectué pour l'ellipsode T représentant le torse, pour l'ellipsode H représentant la 5 tête et pour le cylindre C représentant le cou de l'individu. Simplement, la position du pavillon auditif O est repérée sur l'ellipsode H, de manière à y positionner la source sonore virtuelle, pour la simulation. On obtient ainsi les fonctions de transfert HRTFs pour la gamme de 10 fréquences de 20 Hz à 2000 Hz. En se référant aux figures 6A à 6D, la portion de cône O représentant l'oreille de l'individu est maintenant prise en compte, avec l'ellipsode H représentant sa tête. On 15 obtient alors, en prenant en compte ces dernières formes primitives O et H, les fonctions de transfert HRTFs pour la gamme de fréquences de 2000 Hz à 13000 Hz. On indique que le maillage (non représenté sur les 20 figures) qui est appliqué au modèle des figures 6A à 6D est plus fin que le maillage appliqué au modèle des figures 5A à 5D, comme décrit ci-avant. Une fois que ces fonctions de transfert pour les deux 25 gammes de fréquences basses et hautes ont été calculées, on recompose les fonctions de transfert HRTFs, pour l'oreille droite et pour l'oreille gauche et pour l'ensemble de la gamme de fréquences de 20 Hz à 13000 Hz. Dans une variante du choix des gammes de fréquences, il peut être prévu, pour une gamme de fréquences de 20 Hz à 8000 Hz, de tenir compte dans la modélisation de l'ellipsode H représentant la tête de l'individu, du cylindre C représentant le cou, de l'ellipsode T représentant le torse et de la portion de cône O 5 représentant le pavillon auditif. Dans cette variante, pour la gamme de hautes fréquences de 8000 Hz à 13000 Hz, on considère l'ellipsode H modélisant la tête et la portion de cône O modélisant le pavillon auditif de l'individu. On effectue ensuite une concaténation des fonctions de transfert HRTFs sur l'ensemble de la gamme de fréquences de 20 Hz à 13000 Hz. De préférence, on procède ensuite à une normalisation des HRTFs ainsi calculées et concaténées. Par exemple, les HRTFs sont normalisées par une détermination de la pression acoustique qui est normalement présente, par exemple à la position du centre de la tête de l'individu, 20 mais en l'absence du modèle représentant l'individu. Préférentiellement, cette étape est menée en utilisant un logiciel de calcul tel que MATLABO. Dans ces conditions, on obtient finalement les deux 25 fonctions de transfert HRTFs (oreille droite et oreille gauche) pour un individu. La Demanderesse a constaté systématiquement des résultats satisfaisants en mode statique (synthèse binaurale avec les HRTFs obtenues par simulation sans prise en compte des mouvements de la tête 30 de l'auditeur), et des résultats encore plus satisfaisants en mode dynamique (avec un capteur de positions pour mesurer les mouvements de la tête de l'auditeur, par exemple à l'aide d'un capteur de type "HEAD TRACKER"). Ainsi, dans un premier mode de réalisation, on définit un 5 modèle géométrique de la morphologie d'un auditeur individuel, basé sur des primitives géométriques simples, telles que des sphères, des ellipsodes, des cylindres, des cônes, ou autres. Ces primitives géométriques sont utilisées pour une modélisation acoustique d'une gamme de 10 basses fréquences et d'une gamme de moyennes ou hautes fréquences des fonctions de transfert binaurales HRTFs de cet individu. Dans ce premier mode de réalisation, cette modélisation est menée à l'aide d'un programme de calcul par éléments finis de frontière. Le modèle obtenu permet 15 alors de décrire le comportement acoustique du buste de l'individu auditeur. On indique que la définition des primitives géométriques est effectuée de telle manière que les caractéristiques perceptives des HRTFs soient modélisées de façon efficace dans la gamme de fréquences 20 de 20 Hz, jusqu'à 13 kHz couvrant ainsi presque l'ensemble de la gamme audible (20 Hz à 16 kHz). Cette modélisation, permettant d'obtenir une localisation sonore très efficace, peut être paramétrée de façon très simple néanmoins, puisqu'un jeu restreint de paramètres 25 anthropométriques individualise correctement le modèle de chaque auditeur. Ainsi, à partir de mesures anthropométriques simples, soit à l'aide de photographies, soit à l'aide de mesures directes à partir d'un mètre ou d'un pied à coulisse, par exemple, on évalue les effets 30 acoustiques du buste de l'individu sur deux plages de fréquences, le torse de l'individu étant en particulier modélisé pour la détermination des fonctions de transfert dans les basses fréquences. En effet, la Demanderesse a constaté que le torse de l'individu n'a pas une contribution importante pour les moyennes et hautes fréquences comprises entre 2 kHz et 13 kHz. On se réfère maintenant à la figure 7 pour décrire un mannequin pour la mise en oeuvre du second mode de réalisation précité du procédé selon l'invention. Le 10 mannequin comporte globalement deux ellipsodes H et T, modélisant respectivement la tête et le torse de l'individu. Ces deux ellipsodes H et T sont solidarisés l'un à l'autre par une partie cylindrique C modélisant le cou de l'individu. On indique que chaque partie du 15 mannequin, telle que l'ellipsode H, l'ellipsode T ou le cylindre C, peut se présenter sous la forme d'un assemblage de coques rigides, reliées les unes aux autres par des moyens de réglage des distances entre coques. De tels moyens de réglage peuvent se présenter sous la forme 20 de vis sans fin reliées à un bouton rotatif, en contact par exemple avec une lame ressort. Avantageusement, ces coques rigides sont revêtues d'une mousse de polymère, telle que du latex pour imiter la texture de la peau humaine. Dans l'exemple représenté sur la figure 7, l'ellipsode H représentant la tête de l'individu comporte trois réglages RH1, RH2 et RH3 pour ajuster les trois axes de l'ellipsode en conformité avec les mesures anthropométriques réalisées sur l'individu (par exemple en comparant les vues de face et de profil de l'individu, comme décrit ci-avant en référence aux figures lA et 1B). De même, l'ellipsode T représentant le torse de l'individu comporte trois moyens de réglage RT1, RT2 et RT3 pour ajuster les dimensions des axes de l'ellipsode T 5 en conformité avec les mesures anthropométriques préalablement effectuées (comme décrit ci-avant en référence aux figures 2A et 2B). Dans l'exemple représenté sur la figure 7, le cylindre C représentant le cou de l'individu peut être réglé en hauteur et en largeur, à 10 partir des moyens de réglage RC1 et RC2, respectivement. Bien entendu, un troisième moyen de réglage (non représenté) peut être prévu pour ajuster les axes de l'ellipse correspondant à la base du cylindre. Le pavillon auditif O de l'individu est modélisé par une portion de cône dont la hauteur est réglable par les moyens de réglage RO. On indique qu'il peut être prévu en outre un moyen de réglage des dimensions de la base du cône (non représenté). Par ailleurs, par exemple à l'aide 20 d'une liaison rotule, la portion de cône O peut être disposée dans un cône de débattement ax, en vue de définir précisément la position P de l'entrée du conduit auditif sur le mannequin modélisant l'individu. Cette position du point P peut être déduite des mesures anthropométriques 25 décrites ci-avant en référence aux figures 4A et 4B. C'est sensiblement en cette position P que l'on place ensuite (par exemple à partir d'un casque équipé de transducteurs électroacoustiques) une source sonore ou un microphone pour mesurer les fonctions de transfert HRTFs, sur le 30 mannequin. Ces fonctions de transfert correspondent parfaitement, selon les essais effectués par la Demanderesse, aux fonctions de transfert propres à l'individu modélisé par le mannequin. Bien entendu, la présente invention ne se limite pas à la forme de réalisation décrite ci-avant à titre d'exemple mais elle s'étend à d'autres variantes. Ainsi, on comprendra que les formes géométriques simples des primitives ci-avant sont susceptibles de variantes. 10 Par exemple, dans une réalisation moins sophistiquée, il peut être prévu de modéliser le torse d'un individu, simplement par une demie sphère, ou encore la tête d'un individu par une sphère, ou encore le pavillon auditif par une demie ellipse. En revanche, dans une réalisation plus 15 sophistiquée, il peut être prévu, par exemple pour une tête humaine, de modéliser en outre les pommettes de l'individu par des portions de sphère, ou encore le nez par une portion de cône. Par ailleurs, les gammes de fréquences décrites ci-avant sont susceptibles de variantes, selon la nature ou la complexité des primitives choisies. Les prises de vues ci-avant peuvent être obtenues par 25 photographies de l'individu, ou encore par filmage de l'individu. Dans une autre variante, ces prises de vues peuvent consister encore en un relevé tridimensionnel direct des lignes de niveau du buste d'un individu, notamment à partir d'un balayage laser sur le buste, cette 30 technique de balayage d'un faisceau laser pour obtenir une topographie 3D étant connue en soi. Revendications 1. Procédé de mesure de fonctions de transfert acoustiques associées à la morphologie d'un individu, pour une spatialisation sonore, dans lequel: a) on mesure au moins les dimensions de la tête et du torse d'un individu, b) on définit, en fonction des mesures de l'étape a) , un premier ellipsode de dimensions correspondant 10 sensiblement aux dimensions de la tête de l'individu et un second ellipsode de dimensions correspondant sensiblement aux dimensions du torse de l'individu, et c) on estime: - une première perturbation d'une onde acoustique par au 15 moins lesdits premier et second ellipsodes, pour au moins une fréquence de l'onde acoustique comprise dans une gamme de basses fréquences, pour obtenir un premier jeu de fonctions de transfert dans les basses fréquences acoustiques, - et une seconde perturbation d'une onde acoustique par au moins ledit premier ellipsode, pour au moins une fréquence de l'onde acoustique comprise dans une gamme de hautes fréquences, pour obtenir un second jeu de fonctions de transfert dans les hautes fréquences acoustiques. 25 2. Procédé selon la revendication 1, caractérisé en ce qu'à l'étape a), on repère au moins une position d'une oreille de l'individu, en ce qu'à l'étape b), on définit ladite position sur le 30 premier ellipsode, et en ce qu'à l'étape c), on estime les fonctions de transfert acoustiques pour ladite position de l'oreille. 3. Procédé selon l'une des revendications 1 et 2, 5 caractérisé en ce qu'à l'étape a), on mesure en outre les dimensions du cou de l'individu, en ce qu'à l'étape b), on définit un cylindre de dimensions correspondant sensiblement aux dimensions du cou de l'individu, et en ce qu'à l'étape c) , on estime ladite première perturbation dans la gamme des basses fréquences, à la fois par ledit cylindre et lesdits premier et second ellipsodes. 4. Procédé selon l'une des revendications 1 à 3, caractérisé en ce qu'à l'étape a), on mesure en outre les dimensions d'au moins une oreille de l'individu, ainsi que sa position sur la tête de l'individu, en ce qu'à l'étape b) : - on définit au moins une portion de cône, de dimensions correspondant sensiblement aux dimensions de ladite oreille de l'individu, et - on reporte ladite portion de cône sur le premier ellipsode, et en ce qu'à l'étape c), on estime ladite seconde perturbation dans la gamme des hautes fréquences, à la fois par au moins ladite portion de cône et ledit premier ellipsode. 5. Procédé selon l'une des revendications 1 à 4, caractérisé en ce qu'à l'étape c) : cl) on affecte une impédance acoustique choisie au moins auxdits premier et second ellipsodes, et c2) on calcule, par simulation, lesdites première et seconde perturbations. 6. Procédé selon la revendication 5, caractérisé en ce que ladite impédance acoustique choisie est infinie. 7. Procédé selon l'une des revendications 5 et 6, 10 caractérisé en ce qu'à l'étape c2) : c21) on applique au moins un maillage auxdits premier et second ellipsodes, et c22) on calcule par éléments finis lesdites première et seconde perturbations. 8. Procédé selon la revendication 7, caractérisé en ce qu'à l'étape c21), on applique: - un premier maillage pour le calcul de la première perturbation dans les basses fréquences, et - un second maillage pour le calcul de la seconde perturbation dans les hautes fréquences, et en ce que les dimensions d'une maille du second maillage sont plus petites que les dimensions d'une maille du premier maillage. 9. Procédé selon l'une des revendications 5 à 8, caractérisé en ce que l'étape de calcul c2) est effectuée sur une bande de fréquences logarithmiques, et en ce que l'étape c2) est suivie d'une interpolation sur une échelle 30 linéaire. 10. Procédé selon l'une des revendications 5 à 9, caractérisé en ce qu'à l'étape c2), on simule la présence d'une source acoustique au voisinage de ladite position de l'oreille, lesdites perturbations étant calculées dans un espace entourant les premier et second ellipsodes. 11. Procédé selon la revendication 10, caractérisé en ce que ladite source acoustique est simulée à environ 2 millimètres de ladite position de l'oreille. 12. Procédé selon l'une des revendications 1 à 4, caractérisé en ce qu'à l'étape c) : cl) on prévoit un mannequin comportant un premier ellipsode représentant la tête d'un individu et un second 15 ellipsode, solidaire du premier ellipsode et représentant le torse d'un individu, ainsi que des moyens de réglage des axes desdits ellipsodes, c2) on place au moins une source sonore à une position choisie par rapport au mannequin, et c3) on mesure lesdites première et seconde perturbations d'une onde acoustique que produit ladite source sonore. 13. Procédé selon la revendication 12, caractérisé en ce 25 qu'à l'étape cl), on repère en outre la position de l'oreille de l'individu sur le premier ellipsode du mannequin, et en ce que l'on estime lesdites perturbations à environ 2 millimètres de ladite position de l'oreille, audessus de la surface du premier ellipsode. 14. Procédé selon l'une des revendications précédentes, caractérisé en ce que ladite gamme de basses fréquences s'étend sensiblement entre environ 20 Hz et environ 2000 Hz. 15. Procédé selon l'une des revendications précédentes, caractérisé en ce que ladite gamme de hautes fréquences s'étend sensiblement entre environ 2 kHz et au moins kHz. 16. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'à l'étape a) : ai) on obtient une première prise de vue de l'individu, de profil, avec un étalon de mesure, et a2) on obtient une seconde prise de vue de l'individu, de face, avec ledit étalon de mesure. 17. Procédé selon la revendication 16, caractérisé en ce qu'à l'étape a) : ail) on repère sur la première prise de vue les positions respectives d'un sourcil et d'une oreille de l'individu, a12) on trace une droite passant par l'oreille et le sourcil, ainsi que par un premier point situé à 25 l'arrière de la tête, a13) on détermine une droite médiane et perpendiculaire du segment défini par les positions du sourcil et dudit premier point, et a14) on repère deux intersections de ladite droite médiane 30 avec un contour du visage de l'individu sur ladite première prise de vue, lesdites deux intersections correspondant sensiblement au sommet du crâne de l'individu et au menton de l'individu, et en ce qu'à l'étape b), on définit un premier axe du premier ellipsode entre lesdites deux intersections et un 5 second axe du premier ellipsode perpendiculaire et médian dudit premier axe sur ladite première prise de vue. 18. Procédé selon l'une des revendications 16 et 17, caractérisé en ce qu'à l'étape a) : a21) on repère sur la seconde prise de vue les positions respectives des oreilles de l'individu, a22) on détermine une droite médiane perpendiculaire du segment défini par lesdites positions des oreilles, et a23) on repère deux intersections de ladite droite médiane avec un contour du visage de l'individu sur ladite seconde prise de vue, lesdites deux intersections correspondant sensiblement au sommet du crâne de l'individu et au menton de l'individu, et en ce que, à l'étape b), on définit un troisième axe du premier ellipsode entre lesdites deux intersections et un quatrième axe du premier ellipsode perpendiculaire et médian dudit troisième axe sur ladite seconde prise de vue. 19. Procédé selon la revendication 18, prise en combinaison avec la revendication 17, caractérisé en ce que ledit troisième axe correspond sensiblement au premier axe, et en ce que le premier ellipsode est défini par 30 lesdits premier, second et quatrième axes. 20. Procédé selon l'une des revendications 16 à 19, caractérisé en ce qu'à l'étape a) a3l) on repère sur la seconde prise de vue une position haute du sommet des épaules et une position basse des 5 coudes de l'individu, d'une part, et les limites du bras gauche et du bras droit, d'autre part, a32) on repère sur la première prise de vue les positions respectives de l'arrière du dos et de l'avant de la poitrine de l'individu, a33) on définit un volume délimité par les positions repérées aux étapes a31) et a32), et a34) on ajuste les axes du second ellipsode, d'une part, en orientant dans l'espace ledit ellipsode et, d'autre part, en recherchant un recouvrement 15 volumique maximal avec le volume défini à l'étape a33). 21. Mannequin pour la mise en oeuvre du procédé selon la revendication 12, caractérisé en ce qu'il comporte au 20 moins: - un premier ellipsode représentant la tête d'un individu, - un second ellipsode représentant le torse d'un individu, - des moyens de solidarisation du premier ellipsode au second ellipsode, - des moyens de fixation d'au moins un transducteur électroacoustique en une position choisie sur ledit premier ellipsode, et - des moyens de réglage des axes respectifs desdits premier et second ellipsodes. 22. Mannequin selon la revendication 21, caractérisé en ce que lesdits moyens de solidarisation sont de forme sensiblement cylindrique et représentent sensiblement le 5 cou d'un individu, et en ce qu'ils comportent des moyens de réglage en hauteur et en largeur du cylindre. 23. Mannequin selon l'une des revendications 21 et 22, caractérisé en ce que le premier ellipsode comporte au 10 moins une portion de cône représentant sensiblement une oreille d'un individu, ainsi que des moyens de réglage au moins en hauteur du cône. 24. Mannequin selon l'une des revendications 21 à 23, caractérisé en ce que lesdits premier et second ellipsodes au moins sont revêtus d'une mousse de polymère.