EP1563485A1 - Procede de traitement de donnees sonores et dispositif d'acquisition sonore mettant en oeuvre ce procede - Google Patents

Procede de traitement de donnees sonores et dispositif d'acquisition sonore mettant en oeuvre ce procede

Info

Publication number
EP1563485A1
EP1563485A1 EP03782553A EP03782553A EP1563485A1 EP 1563485 A1 EP1563485 A1 EP 1563485A1 EP 03782553 A EP03782553 A EP 03782553A EP 03782553 A EP03782553 A EP 03782553A EP 1563485 A1 EP1563485 A1 EP 1563485A1
Authority
EP
European Patent Office
Prior art keywords
distance
sound
components
point
signals
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
EP03782553A
Other languages
German (de)
English (en)
Other versions
EP1563485B1 (fr
Inventor
Jérôme DANIEL
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Orange SA
Original Assignee
France Telecom SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=32187712&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=EP1563485(A1) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by France Telecom SA filed Critical France Telecom SA
Publication of EP1563485A1 publication Critical patent/EP1563485A1/fr
Application granted granted Critical
Publication of EP1563485B1 publication Critical patent/EP1563485B1/fr
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10H—ELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H1/00—Details of electrophonic musical instruments
    • G10H1/0091—Means for obtaining special acoustic effects
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10H—ELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
    • G10H1/00—Details of electrophonic musical instruments
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S2400/00—Details of stereophonic systems covered by H04S but not provided for in its groups
    • H04S2400/15—Aspects of sound capture and related signal processing for recording or reproduction
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S2420/00—Techniques used stereophonic systems covered by H04S but not provided for in its groups
    • H04S2420/11—Application of ambisonics in stereophonic audio systems

Definitions

  • the present invention relates to the processing of sound data.
  • Techniques relating to the propagation of a sound wave in three-dimensional space involving in particular a simulation and / or a specialized sound reproduction, implement audio signal processing methods applied to the simulation of acoustic and psycho-acoustic phenomena .
  • Such processing methods provide for spatial encoding of the acoustic field, its transmission and its spatial reproduction on a set of loudspeakers or on headphones of a stereophonic headset.
  • a first category of treatments relates to processes for synthesizing room effects, or more generally environmental effects. From a description of one or more sound sources (signal emitted, position, orientation, directivity, or other) and based on a room effect model (involving room geometry, or even an acoustic perception desired), we calculate and describe a set of elementary acoustic phenomena (direct waves, reflected or diffracted), or even a macroscopic acoustic phenomenon (reverberant and diffuse field), allowing to translate the spatial effect at the level of a listener located at a chosen point of auditory perception, in three-dimensional space. We then calculate a set of signals typically associated with reflections
  • Secondary sources active by re-emission of a main wave received, having a spatial position attribute
  • a late reverberation decorrelated signals for a diffuse field
  • a second category of process concerns the positional or directional rendering of sound sources. These methods are applied to signals determined by a method of the first category described above.
  • Methods according to this second category make it possible to obtain signals to be broadcast over loudspeakers or headphones, in order to finally give a listener the auditory impression of sound sources placed at respective predetermined positions, around the auditor.
  • the methods according to this second category are qualified as "creators of three-dimensional sound images", because of the distribution in three-dimensional space of the feeling of the position of the sources by a listener.
  • Methods according to the second category generally include a first step of spatial encoding of elementary acoustic events which produces a representation of the sound field in three-dimensional space. In a second step, this representation is transmitted or stored for deferred use. In a third decoding step, the decoded signals are delivered to speakers or headphones of a playback device.
  • the present invention falls rather into the second category mentioned above. It concerns in particular the spatial encoding of sound sources and a specification of the three-dimensional sound representation of these sources. It applies as well to an encoding of "virtual" sound sources (applications where sound sources are simulated such as games, a spatial conference, or other), as an "acoustic" encoding of a natural sound field, when taking sound through one or more three-dimensional networks of microphones.
  • Ambisonic encoding which will be described in detail later, consists in representing signals relating to one or more sound waves in a base of spherical harmonics (in spherical coordinates notably implying an elevation angle and an azimuthal angle, characterizing a direction of the sound (s)).
  • the components representing these signals and expressed in this base of spherical harmonics are also a function, for waves emitted in the near field, of a distance between the sound source emitting this field and a point corresponding to the origin of the harmonic base. spherical. More specifically, this distance dependence is expressed as a function of the sound frequency, as will be seen below.
  • This ambisonic approach offers a large number of possible functionalities, in particular in terms of simulation of virtual sources, and, in general, has the following advantages:
  • the representation of acoustic phenomena is scalable: it offers a spatial resolution which can be adapted to different situations. Indeed, this representation can be transmitted and used as a function of rate constraints during the transmission of the encoded signals and / or of limitations of the restitution device;
  • the ambisonic representation is flexible and it is possible to simulate a rotation of the sound field, or even, at restitution, to adapt the decoding of ambisonic signals to any restitution device, of various geometries.
  • the encoding of virtual sources is essentially directional.
  • the encoding functions amount to calculating gains which depend on the incidence of the sound wave expressed by the spherical harmonic functions which depend on the elevation angle and the azimuth angle in spherical coordinates.
  • the loudspeakers, at restitution are far away. he this results in a distortion (or curvature) of the shape of the reconstructed wave fronts.
  • the components of the sound signal in the base of spherical harmonics, for a near field in fact also depend on the distance from the source and on the sound frequency.
  • these components can be expressed mathematically in the form of a polynomial, the variable of which is inversely proportional to the aforementioned distance and to the sound frequency.
  • the ambisonic components in the sense of their theoretical expression, are divergent in the low frequencies and, in particular, tend towards infinity when the sound frequency decreases towards zero, when they represent a sound in near field emitted by a source located at a finite distance.
  • This mathematical phenomenon is known, in the field of ambisonic representation, already for order 1, by the term "bass boost", in particular by: - .A.GERZON, "General Metatheory of Audi tory
  • This phenomenon becomes particularly critical for high spherical harmonic orders involving high power polynomials.
  • spatial aliasing a processing artefact that appears because of the distance between the microphones, unless you choose a virtual microphone mesh tight in space, which weighs down the processing;
  • this technique is difficult to transpose to a real case of sensors to be placed on a network, in the presence of a real source, at acquisition;
  • the three-dimensional sound representation is implicitly subject to a fixed radius of the rendering device because the ambisonic decoding must be done, here, on a network of loudspeakers of the same dimensions as the initial microphone network, this document does not offering no way to adapt encoding or decoding to other sizes of rendering devices.
  • this document presents a horizontal network of sensors, which supposes that the acoustic phenomena which are taken into account here, propagate only in horizontal directions, which excludes any other direction of propagation and which, therefore, does not represent the physical reality of an ordinary sound field.
  • An object of the present invention is to provide a method for processing, by encoding, transmission and reproduction, any type of sound field, in particular the effect of a sound source in the near field.
  • Another object of the present invention is to provide a method allowing the encoding of virtual sources, not only in direction, but also in distance, and to define a decoding adaptable to any rendering device.
  • Another object of the present invention is to provide a robust processing method for sounds of all sound frequencies (including low frequencies), especially for taking sound from natural acoustic fields using three-dimensional networks of microphones.
  • the present invention provides a process for processing sound data, in which: a) signals representative of at least one sound propagating in three-dimensional space and coming from a source located at a first distance are coded from a reference point, to obtain a representation of the sound by components expressed in a base of spherical harmonics, of origin corresponding to said reference point, and b) a compensation for a near field effect is applied to said components by filtering which is a function of a second distance defining substantially, for a sound reproduction by a reproduction device, a distance between a reproduction point and a hearing perception point.
  • said source being far from the reference point
  • a filter is applied, the coefficients of which, each applied to a component of order m, are expressed analytically in the form of the inverse of a polynomial of power m, the variable of which is inversely proportional to the sound frequency and to said second distance, to compensate for a near field effect at the level of the restitution device.
  • said source being a virtual source provided at said first distance
  • the numerator is a polynomial of power m, the variable of which is inversely proportional to the sound frequency and to said first distance, to simulate a field effect close to the virtual source, and
  • the denominator is a polynomial of power m, the variable of which is inversely proportional to the sound frequency and to said second distance, to compensate for the effect of the near field of the virtual source in the low sound frequencies.
  • the coded and filtered data in steps a) and b) are transmitted to the restitution device with a parameter representative of said second distance.
  • the restitution device comprising means for reading a memory medium, the coded and filtered data are stored on a memory medium intended to be read by the restitution device. in steps a) and b) with a parameter representative of said second distance.
  • an adaptation filter is applied to the coded and filtered data, the coefficients of which are a function of said seconds and third distances.
  • the coefficients of this adaptation filter, each applied to a component of order m are expressed analytically in the form of a fraction, of which:
  • the numerator is a polynomial of power m, the variable of which is inversely proportional to the sound frequency and to said second distance,
  • the denominator is a polynomial of power m, the variable of which is inversely proportional to the sound frequency and to said third distance.
  • step b for the implementation of step b), provision is made: - for components of even order m, digital audio filters in the form of a cascade of cells of order two; and
  • digital audio filters in the form of a cascade of cells of order two and an additional cell of order one.
  • the coefficients of a digital audio filter, for a component of order m are defined from the numerical values of the roots of said polynomials of power m.
  • the above-mentioned polynomials are Bessel polynomials.
  • a microphone comprising a network of acoustic transducers arranged substantially on the surface of a sphere whose center corresponds substantially to said reference point, in order to obtain said signals representative of at least one sound. propagating in three-dimensional space.
  • a global filter is applied in step b) on the one hand, to compensate for a near field effect as a function of said second distance and, on the other hand, to equalize the signals coming from the transducers to compensate for a directivity weighting of said transducers.
  • a number of transducers is provided as a function of a chosen total number of components to represent the sound in said base of spherical harmonics.
  • step a) a total number of components is chosen from the base of spherical harmonics to obtain, at the restitution, a region of the space around the point of perception in which the restitution of the sound is faithful and whose dimensions are increasing with the total number of components.
  • a reproduction device comprising a number of loudspeakers at least equal to said total number of components.
  • a reproduction device comprising at least a first and a second loudspeaker arranged at a chosen distance from a listener, for this listener, information is obtained on the feeling expected from the position in the space of sound sources located at a predetermined reference distance from the listener for the application of a technique called "binaural synthesis” or "transaural", and
  • step b) The compensation of step b) is applied with said reference distance substantially as a second distance.
  • a reproduction device comprising at least a first and a second loudspeaker arranged at a chosen distance from a listener, - information is obtained for this listener about feeling the position in the space of sound sources located at a predetermined reference distance from the auditor, and
  • an adaptation filter is applied to the coded and filtered data in steps a) and b), the coefficients of which are a function of the second distance and substantially of the reference distance.
  • the restitution device comprises a helmet with two earphones for the respective ears of the listener
  • steps a) and b) are applied for respective signals intended to supply each earpiece, with, as first distance, respectively a distance separating each ear from a position of a source to be restored in the reproduction space.
  • a matrix system is formed in steps a) and b) comprising at least:
  • the reproduction device upon return the reproduction device comprises a plurality of loudspeakers arranged substantially at the same distance from the point of auditory perception, and
  • a matrix system comprising said matrix result of compensated components and a predetermined decoding matrix, specific to the restitution device, and
  • the present invention also relates to a sound acquisition device, comprising a microphone provided with a network of acoustic transducers arranged substantially on the surface of a sphere.
  • the device further comprises a processing unit arranged for:
  • a filtering which is a function, on the one hand, of a distance corresponding to the radius of the sphere and, on the other hand, of a reference distance.
  • the filtering carried out by the processing unit consists, on the one hand, in equalizing, as a function of the radius of the sphere, the signals coming from the transducers to compensate for a weighting of directivity of said transducers and, on the other hand, in compensate for a near field effect as a function of said reference distance.
  • FIG. 1 schematically illustrates a system of acquisition and creation, by simulation of virtual sources, of sound signals, with encoding, transmission, decoding and restitution by a spatialized restitution device
  • FIG. 4 illustrates a representation by a three-dimensional metric in a coordinate system of spherical coordinates, of spherical harmonics Y ⁇ n of different orders;
  • - Figure 5 is a diagram of the variations of the module of radial functions j a (kr), which are spherical Bessel functions, for successive values of order m, these radial functions intervening in the ambisonic representation of a field of sound pressure;
  • - Figure 6 shows the amplification due to the near field effect for different successive orders m, in particular in the low frequencies;
  • FIG. 7 schematically shows a playback device comprising a plurality of speakers
  • FIG. 8 shows schematically the parameters involved in the ambisonic encoding, with a directional encoding, as well as a distance encoding according to one invention
  • FIG. 11A shows a reconstruction of the near field with compensation, within the meaning of the present invention, for a spherical wave in the horizontal plane;
  • FIG. 11B represents the initial wavefront, coming from a source S;
  • FIG. 12 schematically represents a filtering module for adapting the received and pre-compensated ambisonic components to the encoding for a distance of reference R as a second distance, to a reproduction device comprising a plurality of loudspeakers arranged at a third distance R 2 from a point of auditory perception;
  • - Figure 13A schematically shows the arrangement of a sound source M, during playback, for a listener using a playback device applying binaural synthesis, with a source emitting in the near field;
  • - Figure 13B schematically shows the steps of encoding and decoding with near field effect in the context of the binaural synthesis of Figure 13A which is combined encoding / decoding ambisonic;
  • - Figure 14 schematically shows the processing of signals from a microphone comprising a plurality of pressure sensors arranged on a sphere, for illustrative purposes, by ambisonic encoding, equalization and near field compensation within the meaning of the invention.
  • FIG. 1 shows by way of illustration a global sound spatialization system.
  • a module simulating a virtual scene defines a sound object as a virtual source of a signal, for example monophonic, of position chosen in three-dimensional space and which defines a direction of sound. There may also be specifications for the geometry of a virtual room, to simulate reverberation of sound.
  • a processing module 11 applies management of one or more of these sources with respect to a listener (definition of a virtual position of the sources with respect to this listener). It implements a room effect processor to simulate reverberations or the like by applying common delays and / or filtering.
  • the signals thus constructed are transmitted to a module 2a for spatial encoding of the elementary contributions of the sources.
  • a natural sound recording can be carried out as part of a sound recording by one or more.
  • microphones arranged in a chosen manner in relation to the real sources (module 1b).
  • the signals picked up by the microphones are encoded by a module 2b.
  • the acquired and encoded signals can be transformed according to an intermediate representation format (module 3b), before being mixed by the module 3 with the signals generated by the module la and encoded by the module 2a (from virtual sources).
  • the mixed signals are then transmitted, or even memorized on a medium, for later playback (arrow TR). They are then applied to a decoding module 5, for playback on a playback device 6 comprising speakers.
  • the decoding step 5 can be preceded by a step for manipulating the sound field, for example by rotation, by means of a processing module 4 provided upstream of the decoding module 5.
  • the rendering device can be in the form of a multiplicity of loudspeakers, arranged for example on the surface of a sphere in a three-dimensional (peripheric) configuration to ensure, during the rendering, in particular a feeling of a direction sound in three-dimensional space.
  • an auditor generally places at the center of the sphere formed by the network of loudspeakers, this center corresponding to the point of auditory perception mentioned above.
  • the loudspeakers of the reproduction device can be arranged in a plane (two-dimensional panoramic configuration), the loudspeakers being arranged in particular on a circle and the listener usually placing themselves in the center of this circle.
  • the restitution device can be in the form of a "surround" type device (5.1).
  • the reproduction device can be in the form of a headset with two headphones for a binaural synthesis of the restored sound, which allows the listener to feel a direction of the sources in three-dimensional space, as will be seen in detail below.
  • a reproduction device with two loudspeakers, for feeling in three-dimensional space can also be in the form of a transaural reproduction device, with two loudspeakers arranged at a chosen distance from a listener.
  • FIG. 2 We now refer to FIG. 2 to describe a spatial encoding and a decoding for a three-dimensional sound reproduction, of elementary sound sources.
  • the signal from a source 1 to N is transmitted to a spatial encoding module 2, as well as its position (real or virtual). Its position can be defined as well in terms of incidence (direction of the source seen by the listener) as in terms of distance between this source and a listener.
  • the plurality of signals thus encoded allows to obtain a multi-channel representation of a global sound field.
  • the encoded signals are transmitted (arrow TR) to a sound reproduction device 6, for sound reproduction in three-dimensional space, as indicated above with reference to FIG. 1.
  • the set of weighting factors B ⁇ n which are implicitly a function of the frequency, thus describe the pressure field in the zone considered. For this reason, these factors are called "spherical harmonic components" and represent a frequency expression of the sound (or of the pressure field) in the base of the spherical harmonics Y n .
  • spherical harmonics The angular functions are called “spherical harmonics" and are defined by:
  • the spherical harmonics form an orthonormee base where the scalar products between harmonic components and, in general between two functions F and G, are respectively defined by:
  • Spherical harmonics are real bounded functions, as shown in Figure 4, as a function of the order m and the indices n and ⁇ .
  • the dark and light parts correspond respectively to the positive and negative values of the spherical harmonic functions.
  • the radial functions j m (kr) are spherical Bessel functions, the modulus of which is illustrated for some values of the order m in Figure 5.
  • the ambisonic representation of the sound is however less satisfactory as one moves away from the origin O. This effect becomes critical in particular for high sound frequencies (of short wavelength) . It is therefore advantageous to obtain a number of ambisonic components which is as large as possible, which makes it possible to create a region of space around the point of perception, in which the reproduction of the sound is faithful and whose dimensions are increasing with the total number of components.
  • an ambisonic system takes into account a subset of spherical harmonic components, as described above.
  • a system of order M when it takes into account ambisonic components of index m ⁇ M.
  • the reproduction device comprises loudspeakers arranged on the surface of a sphere ("periphery"), it is in principle possible to use as many harmonics as there are loudspeakers.
  • the pressure signal S designates the pressure signal carried by a plane wave and picked up at point O corresponding to the center of the sphere in FIG. 3 (origin of the base in spherical coordinates).
  • the incidence of the wave is described by the azimuth ⁇ and the elevation ⁇ .
  • the expression of the components of the field associated with this plane wave is given by the relation:
  • the encoding produces signals which differ from the original signal only by a real, finite gain, which corresponds to a purely directional encoding (relation [A3]);
  • the additional filter Fjjjf '( ⁇ ) encodes distance information by introducing, in the expression of the ambisonic components, complex amplitude relationships which depend on the frequency, as expressed in the relation [A5].
  • this additional filter is of the "integrator" type, with an increasing and diverging amplification effect (unbounded) as the sound frequencies decrease towards zero.
  • a reproduction device comprises a plurality of loudspeakers HPi, arranged at the same distance R, in the example described, from an auditory perception point P.
  • R distance
  • the point M corresponds to the position of a source (real or virtual) located at the first distance p, stated above, from the reference point O.
  • a near field pre-compensation is introduced at the very stage of encoding, this compensation involving filters of the form
  • coefficients of this compensation filter are increasing with the frequency of the sound and, in particular, tending towards zero, for the low frequencies.
  • this pre-compensation carried out from the encoding, ensures that the transmitted data are not divergent for the low frequencies.
  • the distance p then represents a distance between a near virtual source M and the point O representing the origin of the spherical base of FIG. 3.
  • a first near field simulation filter is thus applied to simulate the presence of a virtual source at the distance p described above.
  • the terms of the coefficient of this filter diverge in the low frequencies ( Figure 6) and, on the other hand, the distance p above will not necessarily represent the distance between the high- speakers of a playback device and a point P of perception (figure 7).
  • a pre-compensation is applied to the encoding, putting in
  • the near-field pre-compensation of the loudspeakers (placed at the distance R), at the encoding stage, can be combined with a simulated near-field effect of a virtual source placed at a distance p.
  • a total filter is finally brought into play resulting, on the one hand, from the near field simulation, and, on the other hand, from the near field compensation, the coefficients of this filter being able to be expressed analytically by the relation:
  • the total filter given by the relation [Ail] is stable and constitutes the "distance encoding" part in the spatial ambisonic encoding according to the invention, as represented in FIG. 8.
  • the coefficients of these filters correspond to functions of monotonic transfer of the frequency, which tend towards the value 1 in • high frequencies and towards the value (R / p) m in low frequencies.
  • R between an auditory perception point and the HPi loudspeakers is effectively of the order of one or a few meters.
  • total filters are also provided (near field compensation and, if necessary, simulation of a near field) H m ⁇ pc ' c ' (co) which are applied to the ambisonic components, according to their order m, to carry out the encoding of the distance, as represented in FIG. 8.
  • H m ⁇ pc ' c ' (co) which are applied to the ambisonic components, according to their order m, to carry out the encoding of the distance, as represented in FIG. 8.
  • FIG. 11B the propagation of the initial sound wave has been represented from a near-field source situated at a distance p from a point in the acquisition space which corresponds, in the reproduction space , at point P of Figure 7 of auditory perception.
  • the listeners symbolized by diagrammed heads
  • an advantageous method for defining a digital filter from the analytical expression of this filter in the analog domain in continuous time consists of a
  • ⁇ p / c (c being the acoustic speed in the environment, typically 340 m / s in the air).
  • the bilinear transform consists in presenting, for a sampling frequency f s , the relation [Ail] in the form:
  • X m7q are the q successive roots of the Bessel polynomial
  • q l and are expressed in table 1 below, for different orders m, in the respective forms of their real part, their modulus (separated by a comma) and their value (real) when m is odd.
  • Table 1 values P e X m , q ⁇ ⁇ X m, q (and R e ⁇ . X m, m ⁇ when m is odd) of a Bessel polynomial calculated using the MATLAB ⁇ calculation software.
  • Digital filters are thus produced in the form of an infinite impulse response, easily configurable as shown above. It should be noted that an implementation in the form of a finite impulse response can be envisaged and consists in calculating the complex spectrum of the transfer function from the analytical formula, then to deduce a finite impulse response by inverse Fourier transform. We then apply a convolution operation for filtering.
  • FIG. 8 a modified ambisonic representation is defined (FIG. 8), by adopting as transmissible representation of the signals expressed in the frequency domain, in the form:
  • R is a reference distance with which a compensated near field effect is associated and c is the speed of sound (typically 340 m / s in air).
  • c is the speed of sound (typically 340 m / s in air).
  • modified ambisonics has the same scalability properties (schematically represented by transmitted data "surrounded” near the arrow TR in figure 1) and obeys the same rotational transformations of the field (module 4 in figure 1) as the usual ambisonic representation .
  • the decoding operation is adaptable to any restitution device, of radius R 2 , different from the reference distance R above.
  • type filters are applied H ( ⁇ ) ⁇ as described above, but with distance parameters R and R 2 , instead of p and R.
  • R / c is to be memorized (and / or transmit) between encoding and decoding.
  • the filter module shown there is provided for example in a processing unit of a rendering device.
  • the ambisonic components received were pre-compensated for encoding for a reference distance Ri as the second distance.
  • the reproduction device comprises a plurality of loudspeakers arranged at a third distance R 2 from an auditory perception point P, this third distance R 2 being different from the aforementioned second distance R_.
  • the invention also makes it possible to mix several ambisonic representations of sound fields (real and / or virtual sources), the reference distances R of which are different (if necessary with infinite reference distances and corresponding to distant sources).
  • the reference distances R are different (if necessary with infinite reference distances and corresponding to distant sources).
  • we will filter a pre-compensation of all these sources at the smallest reference distance, before mixing the signals ambisonics, which allows the reproduction to obtain a correct definition of the sound relief.
  • the distance encoding with near field pre-compensation is advantageously applied in combination with the focusing processing.
  • the wave emitted by each loudspeaker is defined by a prior processing of "re-encoding" of the ambisonic field at the center of the restitution device, as follows.
  • decoding consists in comparing the original ambisonic signals received by the restitution device, in the form:
  • the number of loudspeakers is greater than or equal to the number of ambisonic components to be decoded and the decoding matrix D is expressed, as a function of the re-encoding matrix C, in the form:
  • the matrixing operation is preceded by a filtering operation which compensates for the near field on each component B mn ⁇ and which can be implemented in digital form, as described above, with reference to the relation [A14] .
  • the "re-encoding" matrix C is specific to the restitution device. Its coefficients can be determined initially by configuration and sound characterization of the reproduction device reacting to a predetermined excitation.
  • the decoding matrix D is also specific to the restitution device. Its coefficients can be determined by the relation [B8]. Using the previous notation where B is the matrix of precompensated ambisonic components, these can be transmitted to the restitution device in matrix form B with:
  • the restitution device then decodes the data received in matrix form B (column vector of the transmitted components) by applying the decoding matrix D to the pre-compensated ambisonic components, to form the signals Si intended to supply the loudspeakers HPi, with:
  • FIG. 13A a listener having a headset with two headphones of a binaural synthesis device is represented.
  • the two ears of the listener are arranged at respective points 0 L (left ear) and 0 R (right ear) in space.
  • the center of the listener's head is arranged at point O and the radius of the listener's head is a value.
  • a sound source must be heard hearing at a point M in space, located at a distance r from the center of the listener's head (and respectively at distances r R from the right ear and r L from the ear left).
  • the direction of the source placed at point M is defined by the vectors F, FR and f j _.
  • binaural synthesis is defined as follows.
  • Each listener has his own ear shape.
  • the perception of a sound in space by this listener is done by learning, from birth, according to the shape of the ears (in particular the shape of the flags and the dimensions of the head) specific to this listener.
  • the perception of a sound in space is manifested inter alia by the fact that the sound reaches one ear, before the other ear, which results in a delay ⁇ between the signals to be emitted by each earpiece of the listening device.
  • restitution applying binaural synthesis.
  • the restitution device is initially configured, for the same listener, by scanning a sound source around its head, at the same distance R from the center of its head. It will thus be understood that this distance R can be considered as a distance between a "restitution point" as stated above and a hearing perception point (here the center O of the listener's head).
  • the index L is associated with the signal to be restored by the earpiece attached to the left ear and the index R is associated with the signal to be restored by the earpiece attached to the right ear.
  • a delay is applied to the initial signal S for each channel intended to produce a signal for a separate listener.
  • These delays ⁇ L and ⁇ R are a function of a maximum delay X MAX which here corresponds to the ratio a / c where a, as indicated previously, corresponds to the radius of the listener's head and c to the speed of sound.
  • these delays are defined as a function of the difference in distance from point O (center of the head) to point M (position of the source whose sound is to be reproduced, in FIG. 13A) and of each ear at this point
  • respective gains S and -R / are also applied to each channel, which are a function of a ratio of the distances from point O to point M and of each ear to point M.
  • Respective modules applied to each channel 2 L and 2 R encode the signals of each channel, in an ambisonic representation, with NFC near field pre-compensation (for "jVear Field Compensation") within the meaning of present invention.
  • the signals originating from the source M are transmitted to the restitution device comprising ambisonic decoding modules, for each channel, 5 L and 5 R .
  • the restitution device comprising ambisonic decoding modules, for each channel, 5 L and 5 R .
  • an ambisonic encoding / decoding is applied, with near field compensation, for each channel (left earpiece, right earpiece) in the reproduction with binaural synthesis (here of "B-FORMAT" type), in split form.
  • the near field compensation is carried out, for each channel, with a distance r L and r R as the first distance between each ear and the position M of the sound source to be restored.
  • a microphone 141 comprises a plurality of transducer capsules, capable of picking up acoustic pressures and reproducing electrical signals Si, ..., S N.
  • CAPi capsules are arranged on a sphere of predetermined radius r (here, a rigid sphere, such as a ping-pong ball for example). The capsules are spaced evenly on the sphere. In practice, we choose the number N of capsules as a function of the order M desired for the ambisonic representation.
  • the near-field pre-compensation can be applied not only for the simulation of virtual source, as indicated above, but also to the acquisition and, more generally, by combining the pre-compensation of field close to all types of treatment involving ambisonic representation.
  • EQ ra is an equalizing filter which compensates for a weighting m which is linked to the directivity of the capsules and which also includes diffraction by the rigid sphere.
  • this equalization filter is not stable and an infinite gain is obtained at very low frequencies.
  • the spherical harmonic components themselves, are not of finite amplitude when the sound field is not limited to a propagation of plane waves, i.e. from from distant sources, as we saw earlier.
  • the signals Si to S N are recovered from the microphone 141. If necessary, a pre-equalization of these signals is applied by a processing module 142.
  • the module 143 makes it possible to express these signals in the ambisonic context, in the form matrix.
  • Module 144 applies the filter of relation [C7] to the components ambisonics expressed as a function of the radius r of the sphere of the microphone 141.
  • the near field compensation is performed for a reference distance R as the second distance.
  • the signals encoded and thus filtered by the module 144 can be transmitted, if necessary, with the parameter representative of the reference distance R / c.
  • the near field compensation within the meaning of the present invention can be applied to all types of processing involving an ambisonic representation.
  • This near field compensation makes it possible to apply the ambisonic representation to a multiplicity of sound contexts where the direction of a source and advantageously its distance must be taken into account.
  • the possibility of representing sound phenomena of all types (near or far fields) in the ambisonic context is ensured by this pre-compensation, due to the limitation to finite real values of the ambisonic components.
  • the present invention is not limited to the embodiment described above by way of example; it extends to other variants.
  • the near field pre-compensation can be integrated, in the encoding, as much for a near source as for a far source.
  • the distance p expressed above will be considered to be infinite, without substantially modifying the expression of the filters H m given above.
  • processing using room effect processors which generally provide decorrelated signals usable for modeling the late diffuse field (late reverberation) can be combined with near-field pre-compensation.
  • the encoding principle within the meaning of the present invention can be generalized to radiation models other than monopolar sources (real or virtual) and / or loudspeakers.
  • any form of radiation in particular a source spread out in space
  • any form of radiation can be expressed by integration of a continuous distribution of point elementary sources.
  • the present invention applies to all types of sound spatialization systems, in particular for "virtual reality” type applications (navigation in virtual scenes in three-dimensional space, games with three-dimensional sound spatialization, "chat” type conversations with sound on the Internet), interface sonifications, audio editing software for recording, mixing and playing music, but also for acquisition, from the use of three-dimensional microphones, for taking musical or cinematographic sound, or for transmitting sound environment on the Internet, for example for sound-activated "WebCam”.
  • “virtual reality” type applications novigation in virtual scenes in three-dimensional space, games with three-dimensional sound spatialization, "chat” type conversations with sound on the Internet
  • interface sonifications audio editing software for recording, mixing and playing music, but also for acquisition, from the use of three-dimensional microphones, for taking musical or cinematographic sound, or for transmitting sound environment on the Internet, for example for sound-activated "WebCam”.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Stereophonic System (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

La présente invention concerne le traitement de données sonores. Selon l'invention : a) on code des signaux représentatifs d'un son se propageant dans l'espace tridimensionnel et issu d'une source située à une première distance (rho) d'un point de référence, pour obtenir une représentation du son par des composantes exprimées dans une base d'harmoniques sphériques, d'origine correspondant audit point de référence, b) et on applique à ces composantes une compensation d'un effet de champ proche par un filtrage qui est fonction d'une seconde distance (R) définissant, pour une restitution du son, une distance entre un point de restitution (HPi) et un point (P) de perception auditive où se place généralement un auditeur.

Description

Procédé de traitement de données sonores et dispositif d'acquisition sonore mettant en œuyre ce procédé
La présente invention concerne le traitement de données sonores .
Des techniques relatives à la propagation d'une onde sonore dans l'espace tridimensionnel, impliquant notamment une simulation et/ou une restitution sonores spécialisées, mettent en œuvre des procédés de traitement du signal audio appliqués à la simulation de phénomènes acoustiques et psycho-acoustiques. De tels procédés de traitement prévoient un encodage spatial du champ acoustique, sa transmission et sa reproduction spatialisée sur un ensemble de hauts-parleurs ou sur des écouteurs d'un casque stéréophonique.
Parmi les techniques de son spatialisé, on distingue deux catégories de traitements complémentaires l'une de l'autre mais qui sont généralement mise en œuvre, l'une et l'autre, au sein d'un même. système.
D'une part, une première catégorie de traitements concerne les procédés de synthèse d'effet de salle, ou plus généralement d'effets environnementaux. A partir d'une description d'une ou plusieurs sources sonores (signal émis, position, orientation, directivité, ou autre) et en se basant sur un modèle d'effet de salle (impliquant une géométrie de salle, ou encore une perception acoustique souhaitée), on calcule et l'on décrit un ensemble de phénomènes acoustiques élémentaires (ondes directes, réfléchies ou diffractées) , ou encore un phénomène acoustique macroscopique (champ réverbéré et diffus) , permettant de traduire l'effet spatial au niveau d'un auditeur situé à un point choisi de perception auditive, dans l'espace tridimensionnel. On calcule alors un ensemble de signaux associés typiquement aux réflexions
(sources " secondaires" , actives par ré-émission d'une onde principale reçue, ayant un attribut de position spatiale) et/ou associés à une réverbération tardive (signaux décorrélés pour un champ diffus) .
D'autre part, une seconde catégorie de procédés concerne le rendu positionnel ou directionnel de sources sonores. Ces procédés sont appliqués à des signaux déterminés par un procédé de la première catégorie décrite ci-avant
(impliquant des sources primaires et secondaires) en fonction de la description spatiale (position de la source) qui leur est associée. En particulier, de tels procédés selon cette seconde catégorie permettent d'obtenir des signaux à diffuser sur des hauts-parleurs ou écouteurs, pour finalement donner à un auditeur l'impression auditive de sources sonores placées à des positions respectives prédéterminées, autour de l'auditeur. Les procédés selon cette seconde catégorie sont qualifiés de " créateurs d ' images sonores tridimensionnelles ", du fait de la répartition dans l'espace tridimensionnel du ressenti de la position des sources par un auditeur. Des procédés selon la seconde catégorie comportent généralement une première étape d'encodage spatial des événements acoustiques élémentaires qui produit une représentation du champ sonore dans l'espace tridimensionnel. Dans une seconde étape, cette représentation est transmise ou stockée pour un usage différé. Dans une troisième étape, de décodage, les signaux décodés sont délivrés sur des hauts-parleurs ou des écouteurs d'un dispositif de restitution.
La présente invention s'inscrit plutôt dans la seconde catégorie précitée. Elle concerne en particulier l'encodage spatial de sources sonores et une spécification de la représentation sonore tridimensionnelle de ces sources. Elle s'applique aussi bien à un encodage de sources sonores " virtuelles" (applications où des sources sonores sont simulées telles que des jeux, une conférence spatialisée, ou autres) , qu'un encodage " acoustique" d'un champ sonore naturel, lors d'une prise de son par un ou plusieurs réseaux tridimensionnels de microphones.
Parmi les techniques envisageables de spatialisation du son, l'approche "amb isonique " est préférée. L'encodage ambisonique, qui sera décrit en détail plus loin, consiste à représenter des signaux relatifs à une ou plusieurs ondes sonores dans une base d'harmoniques sphériques (en coordonnées sphériques impliquant notamment un angle d'élévation et un angle azimutal, caractérisant une direction du ou des sons) . Les composantes représentant ces signaux et exprimées dans cette base d'harmoniques sphériques sont aussi fonction, pour les ondes émises en champ proche, d'une distance entre la source sonore émettant ce champ et un point correspondant à 1 ' origine de la base des harmonique sphériques. Plus particulièrement, cette dépendance de la distance s'exprime en fonction de la fréquence sonore, comme on le verra plus loin.
Cette approche ambisonique offre un grand nombre de fonctionnalités possibles, notamment en terme de simulation de sources virtuelles, et, de manière générale, présente les avantages suivants :
- elle traduit, de façon rationnelle, la réalité des phénomènes acoustiques et apporte un rendu auditif spatial réaliste, convaincant et immersif ;
- la représentation des phénomènes acoustiques est scalable : elle offre une résolution spatiale qui peut être adaptée à différentes situations. En effet, cette représentation peut être transmise et exploitée en fonction de contraintes de débit lors de la transmission des signaux encodés et/ou de limitations du dispositif de restitution ;
- la représentation ambisonique est flexible et il est possible simuler une rotation du champ sonore, ou encore, à la restitution, d'adapter le décodage des signaux ambisoniques à tout dispositif de restitution, de géométries diverses .
Dans l'approche ambisonique connue, l'encodage des sources virtuelles est essentiellement directionnel. Les fonctions d'encodage reviennent à calculer des gains qui dépendent de l'incidence de l'onde sonore exprimée par les fonctions harmoniques sphériques qui dépendent de 1 ' angle d'élévation et de l'angle azimutal en coordonnées sphériques. En particulier, au décodage, on suppose que les hauts-parleurs, à la restitution, sont lointains. Il en résulte une distorsion (ou une incurvation) de la forme des fronts d'onde reconstruits. En effet, comme indiqué ci-avant, les composantes du signal sonore dans la base des harmoniques sphériques, pour un champ proche, dépendent en fait aussi de la distance de la source et de la fréquence sonore. Plus précisément, ces composantes peuvent s'exprimer mathématiquement sous la forme d'un polynôme dont la variable est inversement proportionnelle à la distance précitée et à la fréquence sonore. Ainsi, les composantes ambisoniques, au sens de leur expression théorique, sont divergentes dans les basses fréquences et, en particulier, tendent vers l'infini quand la fréquence sonore décroît vers zéro, lorsqu'elles représentent un son en champ proche émis par une source située à une distance finie. Ce phénomène mathématique est connu, dans le domaine de la représentation ambisonique, déjà pour l'ordre 1, par le terme de "bass boost " , notamment par : - .A.GERZON, " General Metatheory of Audi tory
Localisation" , preprint 3306 of the 92nd AES Convention, 1992, page 52.
Ce phénomène devient particulièrement critique pour des ordres d'harmoniques sphériques élevées impliquant des polynômes de puissance élevée.
On connaît par :
SONTACCHI et HÔLDRICH, "Furtizer Investigations on 3D Sound Fields using Distance Coding" (Proceedings of the COST G- 6 Conférence on Digital Audio Effects (DAFX-01) , Limerick, Irlande, 6-8 Décembre 2001) , une technique pour prendre en compte une incurvation des fronts d'ondes au sein d'une représentation proche d'une représentation ambisonique, dont le principe consiste à :
- appliquer un encodage ambisonique (d'ordre élevé) aux signaux issus d'une prise de son virtuelle (simulée) , de type WFS (pour "Wave Field Synthesis ") ;
- et reconstruire le champ acoustique sur une zone d'après ses valeurs sur une frontière de zone, se fondant ainsi sur le principe de HUYGENS-FRESNEL.
Cependant, la technique présentée dans ce document, bien que prometteuse du fait qu'elle utilise une représentation ambisonique à un ordre élevé, pose un certain nombre de problèmes : - les ressources informatiques nécessaires pour le calcul de toutes les surfaces permettant d'appliquer le principe de HUYGENS-FRESNEL, ainsi que les temps de calcul nécessaires, sont excessifs ;
- des artefacts de traitement dits d' "aliasing spatial " apparaissent à cause de la distance entre les microphones, à moins de choisir un maillage de microphone virtuels serré dans l'espace, ce qui alourdit les traitements ;
- cette technique est difficilement transposable à un cas réel de capteurs à disposer en réseau, en présence d'une source réelle, à l'acquisition ;
- à la restitution, la représentation sonore tridimensionnelle est implicitement assujettie à un rayon figé du dispositif de restitution car le décodage ambisonique doit se faire, ici, sur un réseau de hauts- parleurs de mêmes dimensions que le réseau de microphones initial, ce document ne proposant aucun moyen d'adapter l'encodage ou le décodage à d'autres tailles de dispositifs de restitution.
Surtout, ce document présente un réseau horizontal de capteurs, ce qui suppose que les phénomènes acoustiques dont on tient compte, ici, ne se propagent que dans des directions horizontales, ce qui exclut toute autre direction de propagation et qui, donc, ne représente pas la réalité physique d'un champ acoustique ordinaire.
De façon plus générale, les techniques actuelles ne permettent pas de traiter de façon satisfaisante tout type de sources sonores, notamment en champ proche, mais plutôt des sources sonores lointaines (ondes planes) , ce qui correspond à une situation restrictive et artificielle dans nombreuses applications.
Un objet de la présente invention est de fournir un procédé pour traiter, par encodage, transmission et restitution, un type quelconque de champ sonore, en particulier l'effet d'une source sonore en champ proche.
Un autre objet de la présente invention est de fournir un procédé permettant l'encodage de sources virtuelles, non seulement en direction, mais aussi en distance, et de définir un décodage adaptable à un dispositif de restitution quelconque.
Un autre objet de la présente invention est de fournir un procédé de traitement robuste pour des sons de toutes fréquences sonores (y compris les basses fréquences) , notamment pour la prise de son de champs acoustiques naturels à l'aide de réseaux tridimensionnels de microphones .
A cet effet, la présente invention propose un procédé de traitement de données sonores, dans lequel : a) on code des signaux représentatifs d'au moins un son se propageant dans l'espace tridimensionnel et issu d'une source située à une première distance d'un point de référence, pour obtenir une représentation du son par des composantes exprimées dans une base d'harmoniques sphériques, d'origine correspondant audit point de référence, et b) on applique auxdites composantes une compensation d'un effet de champ proche par un filtrage qui est fonction d'une seconde distance définissant sensiblement, pour une restitution du son par un dispositif de restitution, une distance entre un point de restitution et un point de perception auditive.
Dans un premier mode de réalisation, ladite source étant lointaine du point de référence,
- on obtient des composantes d'ordres successifs m pour la représentation du son dans ladite base d'harmoniques sphériques, et
- on applique un filtre dont les coefficients, appliqués chacun à une composante d'ordre m, s'expriment analytiquement sous la forme de 1 ' inverse d ' un polynôme de puissance m, dont la variable est inversement proportionnelle à la fréquence sonore et à ladite seconde distance, pour compenser un effet de champ proche au niveau du dispositif de restitution.
Dans un second mode de réalisation, ladite source étant une source virtuelle prévue à ladite première distance,
- on obtient des composantes d'ordres successifs m pour la représentation du son dans ladite base d'harmoniques sphériques, et
- on applique un filtre global dont les coefficients, appliqués chacun à une composante d'ordre m, s'expriment analytiquement sous la forme d'une fraction, dont :
- le numérateur est un polynôme de puissance m, dont la variable est inversement proportionnelle à la fréquence sonore et à ladite première distance, pour simuler un effet de champ proche de la source virtuelle, et
- le dénominateur est un polynôme de puissance m, dont la variable est inversement proportionnelle à la fréquence sonore et à ladite seconde distance, pour compenser l'effet du champ proche de la source virtuelle dans les basses fréquences sonores.
Préférentiellement, on transmet au dispositif de restitution les données codées et filtrées aux étapes a) et b) avec un paramètre représentatif de ladite seconde distance.
En complément ou en variante, le dispositif de restitution comportant des moyens de lecture d'un support mémoire, on mémorise sur un support mémoire destiné à être lu par le dispositif de restitution les données codées et filtrées aux étapes a) et b) avec un paramètre représentatif de ladite seconde distance.
Avantageusement, préalablement à une restitution sonore par un dispositif de restitution comportant une pluralité de hauts-parleurs disposés à une troisième distance dudit point de perception auditive, on applique aux données codées et filtrées un filtre d'adaptation dont les coefficients sont fonction desdites seconde et troisième distances.
Dans une réalisation particulière, les coefficients de ce filtre d'adaptation, appliqués chacun à une composante d'ordre m, s'expriment analytiquement sous la forme d'une fraction, dont :
- le numérateur est un polynôme de puissance m, dont la variable est inversement proportionnelle à la fréquence sonore et à ladite seconde distance,
- et le dénominateur est un polynôme de puissance m, dont la variable est inversement proportionnelle à la fréquence sonore et à ladite troisième distance.
Avantageusement, pour la mise en œuvre de l'étape b) , on prévoit : - pour des composantes d'ordre m pair, des filtres audionumériques sous la forme d'une cascade de cellules d'ordre deux ; et
- pour des composantes d'ordre m impair, des filtres audionumériques sous la forme d'une cascade de cellules d'ordre deux et une cellule supplémentaire d'ordre un. Dans cette réalisation, les coefficients d'un filtre audionumérique, pour une composante d'ordre m, sont définis à partir des valeurs numériques des racines desdits polynômes de puissance m.
Dans une réalisation particulière, les polynômes précités sont des polynômes de Bessel.
A 1 'acquisition des signaux sonores, on prévoit avantageusement un microphone comportant un réseau de transducteurs acoustiques agencés sensiblement sur la surface d'une sphère dont le centre correspond sensiblement audit point de référence, pour obtenir lesdits signaux représentatifs d'au moins un son se propageant dans l'espace tridimensionnel.
Dans cette réalisation, on applique à l'étape b) un filtre global pour, d'une part, compenser un effet de champ proche en fonction de ladite seconde distance et, d'autre part, égaliser les signaux issus des transducteurs pour compenser une pondération de directivité desdits transducteurs .
Préférentiellement, on prévoit un nombre de transducteurs fonction d'un nombre total choisi de composantes pour représenter le son dans ladite base d'harmoniques sphériques .
Selon une caractéristique avantageuse, on choisit à l'étape a) un nombre total de composantes dans la base des harmoniques sphériques pour obtenir, à la restitution, une région de l'espace autour du point de perception dans laquelle la restitution du son est fidèle et dont les dimensions sont croissantes avec le nombre total de composantes .
Préférentiellement, on prévoit en outre un dispositif de restitution comportant un nombre de haut-parleurs au moins égal audit nombre total de composantes .
En variante, dans le cadre d'une restitution avec synthèse binaurale ou transaurale :
- on prévoit un dispositif de restitution comportant au moins un premier et un second haut-parleur disposés à une distance choisie d'un auditeur, on obtient, pour cet auditeur, une information de ressenti attendu de la position dans l'espace de sources sonores situées à une distance de référence prédéterminée de l'auditeur pour l'application d'une technique dite de "synthèse binaurale " ou "transaurale " , et
- on applique la compensation de l'étape b) avec ladite distance de référence sensiblement en tant que seconde distance.
Dans une variante où l'on introduit une adaptation au dispositif de restitution à deux écouteurs :
- on prévoit un dispositif de restitution comportant au moins un premier et un second haut-parleur disposés à une distance choisie d'un auditeur, - on obtient, pour cet auditeur, une information de ressenti de la position dans l'espace de sources sonores situées à une distance de référence prédéterminée de l'auditeur, et
- préalablement à une restitution sonore par le dispositif de restitution, on applique aux données codées et filtrées aux étapes a) et b) un filtre d'adaptation dont les coefficients sont fonction de la seconde distance et sensiblement de la distance de référence.
En particulier, dans le cadre d'une restitution avec synthèse binaurale :
- le dispositif de restitution comporte un casque à deux écouteurs pour les oreilles respectives de l'auditeur,
- et préférentiellement, séparément pour chaque écouteur, on applique le codage et le filtrage des étapes a) et b) pour des signaux respectifs destinés à alimenter chaque écouteur, avec, en tant que première distance, respectivement une distance séparant chaque oreille d'une position d'une source à restituer dans l'espace de restitution.
Préférentiellement, on met en forme, aux étapes a) et b) , un système matriciel comportant au moins :
- une matrice comportant lesdites composantes dans la base des harmoniques sphériques, et - une matrice diagonale dont les coefficients correspondent à des coefficients de filtrage de l'étape b), et on multiplie lesdites matrices pour obtenir une matrice résultat de composantes compensées .
De préférence, à la restitution - le dispositif de restitution comporte une pluralité de haut-parleurs disposés sensiblement à une même distance du point de perception auditive, et
- pour décoder lesdites données codées et filtrées aux étapes a) et b) et former des signaux adaptés pour alimenter lesdits haut-parleurs :
* on forme un système matriciel comportant ladite matrice résultat de composantes compensées et une matrice de décodage prédéterminée, propre au dispositif de restitution, et
* on obtient une matrice comportant des coefficients représentatifs des signaux d'alimentation des hauts- parleurs par multiplication de la matrice résultat par ladite matrice de décodage.
La présente invention vise aussi un dispositif d'acquisition sonore, comportant un microphone muni d'un réseau de transducteurs acoustiques disposés sensiblement sur la surface d'une sphère. Selon l'invention le dispositif comporte en outre une unité de traitement agencée pour :
- recevoir des signaux émanant chacun d'un transducteur,
- appliquer auxdits signaux un codage pour obtenir une représentation du son par des composantes exprimées dans une base d'harmoniques sphériques, d'origine correspondant au centre de ladite sphère,
- et appliquer auxdites composantes un filtrage qui est fonction, d'une part, d'une distance correspondant au rayon de la sphère et, d'autre part, d'une distance de référence . Préférentiellement, le filtrage effectué par l'unité de traitement consiste, d'une part, à égaliser, en fonction du rayon de la sphère, les signaux issus des transducteurs pour compenser une pondération de directivité desdits transducteurs et, d'autre part, à compenser un effet de champ proche en fonction de ladite distance de référence.
D'autres avantages et caractéristiques de l'invention apparaîtront à la lecture de la description détaillée ci- après et à l'examen des figures qui l'accompagnent, sur lesquelles :
- la figure 1 illustre schématiquement un système d'acquisition et création, par simulation de sources virtuelles, de signaux sonores, avec encodage, transmission, décodage et restitution par un dispositif de restitution spatialisé,
- la figure 2 représente plus précisément un encodage de signaux définis à la fois en intensité et par rapport à la position d'une source dont ils sont issus, - la figure 3 illustre les paramètres en jeu dans la représentation ambisonique, en coordonnées sphériques ;
- la figure 4 illustre une représentation par une métrique tridimensionnelle dans un repère de coordonnées sphériques, d'harmoniques sphériques Y^n de différents ordres ;
- la figure 5 est un diagramme des variations du module de fonctions radiales ja (kr) , qui sont des fonctions de Bessel sphériques, pour des valeurs d'ordre m successives, ces fonctions radiales intervenant dans la représentation ambisonique d'un champ de pression acoustique ; - la figure 6 représente l'amplification due à l'effet de champ proche pour différents ordres successifs m, en particulier dans les basses fréquences ;
- la figure 7 représente schématiquement un dispositif de restitution comportant une pluralité de hauts-parleurs
HPi , avec le point (référencé P) de perception auditive précité, la première distance précitée (référencée p) et la seconde distance précitée (référencée R) ;
- la figure 8 représente schématiquement les paramètres mis en jeu dans l'encodage ambisonique, avec un encodage directionnel, ainsi qu'un encodage de distance selon 1 ' invention ;
- la figure 9 représente des spectres d'énergie des filtres de compensation et de champ proche simulés pour une première distance d'une source virtuelle p = 1 m et une pré-compensation de hauts-parleurs situés à une seconde distance R = 1,5 m ;
- la figure 10 représente des spectres d'énergie des filtres de compensation et de champ proche simulés pour une première distance de la source virtuelle p = 3 m et une pré-compensation de hauts-parleurs situés à une distance R = 1,5 m ;
- la figure 11A représente une reconstruction du champ proche avec compensation, au sens de la présente invention, pour une onde sphérique dans le plan horizontal ;
- la figure 11B, à comparer avec la figure 11A, représente le front d'onde initial, issu d'une source S ;
- la figure 12 représente schématiquement un module de filtrage pour adapter les composantes ambisoniques reçues et pré-compensées à l'encodage pour une distance de référence R en tant que seconde distance, à un dispositif de restitution comportant une pluralité de hauts-parleurs disposés à une troisième distance R2 d'un point de perception auditive ; - la figure 13A représente schématiquement la disposition d'une source sonore M, à la restitution, pour un auditeur utilisant un dispositif de restitution appliquant une synthèse binaurale, avec une source émettant en champ proche ; - la figure 13B représente schématiquement les étapes d'encodage et de décodage avec effet de champ proche dans le cadre de la synthèse binaurale de la figure 13A à laquelle est combiné un encodage/décodage ambisonique ; - la figure 14 représente schématiquement le traitement des signaux issus d'un microphone comportant une pluralité de capteurs de pression agencés sur une sphère, à titre illustrâtif, par encodage ambisonique, égalisation et compensation de champ proche au sens de l'invention.'
On se réfère tout d'abord à la figure 1 qui représente à titre illustratif un système global de spatialisation sonore. Un module la de simulation d'une scène virtuelle définit un objet sonore comme une source virtuelle d'un signal, par exemple monophonique, de positiqn choisie dans l'espace tridimensionnel et qui définit une direction du son. Il peut être prévu en outre des spécifications de la géométrie d'une salle virtuelle, pour simuler une réverbération du son. Un module de traitement 11 applique une gestion d'une ou plusieurs de ces sources par rapport à un auditeur (définition d'une position virtuelle des sources par rapport à cet auditeur) . Il met en œuvre un processeur d'effet de salle pour simuler des réverbérations ou autres en appliquant des retards et/ou des filtrages usuels. Les signaux ainsi construits sont transmis à un module 2a d'encodage spatial des contributions élémentaires des sources .
Parallèlement, une prise de son naturelle peut être effectuée dans le cadre d'un enregistrement sonore par un ou plusieurs . microphones disposés de façon choisie par rapport aux sources réelles (module 1b) . Les signaux captés par les microphones sont encodés par un module 2b. Les signaux acquis et encodés peuvent être transformés selon un format de représentation intermédiaire (module 3b), avant d'être mixés par le module 3 aux signaux générés par le module la et encodés par le module 2a (issu des sources virtuelles) . Les signaux mixés sont ensuite transmis, ou encore mémorisés sur un support, en vue d'une restitution ultérieure (flèche TR) . Ils sont ensuite appliqués à un module de décodage 5, en vue de la restitution sur un dispositif de restitution 6 comportant des hauts-parleurs . Le cas échéant, l'étape de décodage 5 peut être précédée d'une étape de manipulation du champ sonore, par exemple par rotation, grâce à un module de traitement 4 prévu en amont du module de décodage 5.
Le dispositif de restitution peut se présenter sous la forme d'une multiplicité de hauts-parleurs, agencés par exemple à la surface d'une sphère dans une configuration tridimensionnelle (périphonique) pour assurer, à la restitution, notamment un ressenti d'une direction du son dans l'espace tridimensionnel. A cet effet, un auditeur se place généralement au centre de la sphère formée par le réseau de haut-parleurs, ce centre correspondant au point de perception auditive cité ci-avant. En variante, les hauts-parleurs du dispositif de restitution peuvent être agencés dans un plan (configuration panoramique bidimensionnelle) , les hauts-parleurs étant disposés en particulier sur un cercle et l'auditeur se plaçant habituellement au centre de ce cercle. Dans une autre variante, le dispositif de restitution peut se présenter sous la forme d'un dispositif de type "surround " (5.1). Enfin, dans une variante avantageuse, le dispositif de restitution peut se présenter sous la forme d'un casque à deux écouteurs pour une synthèse binaurale du son restitué, qui permet à l'auditeur de ressentir une direction des sources dans l'espace tridimensionnel, comme on le verra plus loin de façon détaillée. Un tel dispositif de restitution à deux hauts-parleurs, pour un ressenti dans l'espace tridimensionnel, peut se présenter aussi sous la forme d'un dispositif de restitution transaurale, à deux hauts-parleurs disposés à une distance choisie d'un auditeur.
On se réfère maintenant à la figure 2 pour décrire un encodage spatial et un décodage pour une restitution sonore tridimensionnelle, de sources sonores élémentaires. On transmet à un module d'encodage spatial 2 le signal issu d'une source 1 à N, ainsi que sa position (réelle ou virtuelle) . Sa position peut être aussi bien définie en terme d'incidence (direction de la source vue de l'auditeur) qu'en terme de distance entre cette source et un auditeur. La pluralité des signaux ainsi encodés permet d'obtenir une représentation multi-canale d'un champ sonore global. Les signaux encodés sont transmis (flèche TR) à un dispositif de restitution sonore 6, pour une restitution sonore dans l'espace tridimensionnel, comme indiqué ci-avant en référence à la figure 1.
On se réfère maintenant à la figure 3 pour décrire ci- après la représentation ambisonique par des harmoniques sphériques dans l'espace tridimensionnel, d'un champ acoustique. On considère une zone autour d'une origine 0 (sphère de rayon R) exempte de source acoustique. On adopte un système de coordonnées sphériques dans lequel chaque vecteur F dès 1 ' origine 0 à un point de la sphère est décrit par un azimut θr, une élévation δr et un rayon r (correspondant à la distance à l'origine 0).
Le champ de pression p(r) à l'intérieur de cette sphère (r < R où R est le rayon de la sphère) peut s'écrire dans le domaine fréquentiel comme une série dont les termes sont les produits pondérés de fonctions angulaires y^n (θ,δ) et de fonction radiale jm(kr) qui dépendent ainsi d'un terme de propagation où k=2πf/c, où f est la fréquence sonore et c est la vitesse du son dans le milieu de propagation.
Le champ de pression s'exprime alors par
L ' ensemble des facteurs de pondération B ιn , qui sont implicitement fonction de la fréquence, décrivent ainsi le champ de pression dans la zone considérée. Pour cette raison, ces facteurs sont appelés "composantes harmoniques sphériques " et représentent une expression fréquentielle du son (ou du champ de pression) dans la base des harmoniques sphériques Y n .
Les fonctions angulaires sont appelées "harmoniques sphériques " et sont définies par :
[A2]
ou
Pmn isinδ) sont des fonctions de Legendre de degré m et d'ordre n ; δp,q est le symbole de Krônecker (égal à 1 si p=q et 0, sinon)
Les harmoniques sphériques forment une base orthonormee où les produits scalaires entre composantes harmoniques et, de façon générale entre deux fonctions F et G, sont respectivement définies par :
γυ 1mn Ym'n' } = δmm δnnyδσσκ . [ A ' 2 ]
4π = -~- F(θ,δ)G(θ,δ),dΩ(θ,δ)
Les harmoniques sphériques sont des fonctions réelles bornées, comme représenté sur la figure 4, en fonction de l'ordre m et des indices n et σ. Les parties sombres et claires correspondent respectivement aux valeurs positives et négatives des fonctions harmoniques sphériques. Plus l'ordre m est élevé et plus la fréquence angulaire (et donc la discrimination entre fonctions) est élevée. Les fonctions radiales jm (kr) sont des fonctions de Bessel sphériques, dont le module est illustré pour quelques valeurs de l'ordre m dans la figure 5.
On peut donner une interprétation de la représentation ambisonique par une base d'harmoniques sphériques comme suit. Les composantes ambisoniques de même ordre m expriment finalement des "dérivées " ou des "moments " d'ordre m du champ de pression au voisinage de l'origine O (centre de la sphère représentée sur la figure 3) .
En particulier, BQQ = W décrit la grandeur scalaire de la
pression, tandis que = Y, B Q = Z sont liés aux gradients de pression (ou encore à la vélocité particulaire), à l'origine O. Ces quatre premières composantes , X, Y et Z sont obtenues lors d'une prise de son naturelle à l'aide de microphones omnidirectifs (pour la composante W d'ordre 0) et bidirectifs (pour les trois autres composantes suivantes) . En utilisant un plus grand nombre de transducteurs acoustiques, un traitement approprié, notamment par égalisation, permet d'obtenir d'avantage de composantes ambisoniques (ordres m plus élevés supérieurs à 1) .
En prenant en compte des composantes supplémentaires d'ordre plus élevé (supérieur à 1), donc en augmentant la résolution angulaire de la description ambisonique, on accède à une approximation du champ de pression sur un voisinage plus large au regard de la longueur d'onde de l'onde sonore, autour de l'origine O. On comprendra ainsi qu'il existe une relation étroite entre la résolution angulaire (ordre des harmoniques sphériques) et la portée radiale (rayon r) qui peut être représentée. En bref, lorsque l'on s'écarte spatialement du point d'origine O de la figure 3, plus le nombre de composantes ambisoniques est élevé (ordre M élevé) et meilleure est la représentation du son par 1 ' ensemble de ces composantes ambisoniques . On comprendra aussi que la représentation ambisonique du son est toutefois moins satisfaisante au fur et à mesure que l'on s'éloigne de l'origine O. Cet effet devient critique en particulier pour des fréquences sonores élevées (de longueur d'onde courte) . On a donc intérêt à obtenir un nombre de composantes ambisoniques qui soit le plus grand possible, ce qui permet de créer une région de l'espace autour du point de perception, dans laquelle la restitution du son est fidèle et dont les dimensions sont croissantes avec le nombre total de composantes .
On décrit ci-après une application à système d'encodage/transmission/restitution d'un son spatialisé. En pratique, un système ambisonique prend en compte un sous-ensemble de composantes harmoniques sphériques, comme décrit ci-avant. On parle d'un système d'ordre M lorsque celui-ci prend en compte des composantes ambisoniques d'indice m < M. Lorsqu'il s'agit d'une restitution par un dispositif de restitution à hauts-parleurs, on comprendra que si ces hauts-parleurs sont disposés dans un plan horizontal, seules les harmoniques d'indice m=n sont exploitées. En revanche, lorsque le dispositif de restitution comporte des hauts-parleurs disposés sur la surface d'une sphère ( "périphonie" ) , on peut en principe exploiter autant d'harmoniques qu'il existe de haut- parleurs .
On désigne par la référence S le signal de pression porté par une onde plane et capté au point O correspondant au centre de la sphère de la figure 3 (origine de la base en coordonnées sphériques). L'incidence de l'onde est décrite par l'azimut θ et l'élévation δ. L'expression des composantes du champ associé à cette onde plane est donnée par la relation :
Bmn - S.Ymn(θ,δ) [A3]
Pour encoder (simuler) une source en champ proche à une
distance p de l'origine O, on applique un filtre Ft {pic) m pour "incurver" la forme des fronts d'onde, en considérant qu'un champ proche émet, en première approximation, une onde sphérique. Les composantes encodées du champ deviennent :
et l'expression du filtre précité F}f ' est donnée par la relation :
jωPic)-" wι où ω = 2πf est la pulsation de l'onde, f étant la fréquence du son.
Ces deux dernières relations [A4] et [A5] montrent finalement que, aussi bien pour une source virtuelle (simulée) que pour une source réelle en champ proche, les composantes du son dans la représentation ambisonique s'expriment mathématiquement (en particulier analytiquement) sous la forme d'un polynôme, ici de Bessel, de puissance m et dont la variable (c/2jωp) est inversement proportionnelle à la fréquence sonore.
Ainsi, on comprendra que :
- dans le cas d'une onde plane, l'encodage produit des signaux qui ne diffèrent du signal d'origine que d'un gain réel, fini, ce qui correspond à un encodage purement directionnel (relation [A3] ) ;
- dans le cas d'une onde sphérique (source en champ proche) , le filtre supplémentaire Fjjjf '(ω) encode l'information de distance en introduisant, dans l'expression des composantes ambisoniques, des rapports d'amplitudes complexes qui dépendent de la fréquence, comme exprimé dans la relation [A5] .
Il est à noter que ce filtre supplémentaire est de type "intégrateur" , avec un effet d'amplification croissant et divergent (non-borné) au fur et à mesure que les fréquences sonores décroissent vers zéro. La figure 6 montre, pour chaque ordre m, une augmentation du gain en basses fréquences (ici la première distance p = lm) . Il s'agit donc de filtres instables et divergents lorsque l'on cherche à les appliquer à des signaux audio quelconques. Cette divergence est d'autant plus critique pour les ordres m de valeur élevée.
On comprendra en particulier, à partir des relations [A3] , [A4] , et [A5] , que la modélisation d'une source virtuelle en champ proche présente des composantes ambisoniques divergentes en basses fréquences, de façon particulièrement critique pour des ordres m élevés, comme représenté sur la figure 6. Cette divergence, dans les basses fréquences, correspond au phénomène de "bass boost" énoncé ci-avant. Il se manifeste aussi en acquisition sonore, pour des sources réelles.
Pour cette raison notamment, l'approche ambisonique, en particulier pour des ordres m élevés, n'a pas connu, dans l'état de la technique, une application concrète (autre que théorique) dans le traitement du son. On comprend en particulier qu'une compensation du champ proche est nécessaire pour respecter, à la restitution, la forme des fronts d'ondes encodés dans- la représentation ambisonique. En se référant à la figure 7, un dispositif de restitution comporte une pluralité de hauts-parleurs HPi, disposés à une même distance R, dans l'exemple décrit, d'un point de perception auditive P. Sur cette figure 7 :
- chaque point où se situe un haut-parleur HPi correspond à un point de restitution énoncé ci-avant,
- le point P est le point de perception auditive énoncé ci-avant,
- ces points sont séparés de la seconde distance R énoncée ci-avant, tandis que sur la figure 3 décrite ci-avant :
- le point 0 correspond au point de référence, énoncé ci- avant, qui forme l'origine de la base des harmoniques sphériques ,
- le point M correspond à la position d'une source (réelle ou virtuelle) située à la première distance p, énoncée ci- avant, du point de référence O.
Selon l'invention, on introduit une pré-compensâtion du champ proche au stade même de l'encodage, cette compensation mettant en jeu des filtres de la forme
analytique et qui s'appliquent aux composantes ambisoniques B^n précitées. Selon l'un des avantages que procure l'invention, l'amplification Fm ' ( )) dont l'effet apparaît sur la figure 6 est compensée par l'atténuation du filtre
appliqué dès l'encodage . En particulier, les
coefficients de ce filtre de compensation sont croissants avec la fréquence du son et, en particulier, tendent vers zéro, pour les basses fréquences.
Avantageusement, cette pré-compensation, effectuée dès l'encodage, assure que les données transmises ne sont pas divergentes pour les basses fréquences.
Pour indiquer la signification physique de la distance R qui intervient dans le filtre de compensation, on considère, à titre illustratif, une onde plane réelle, initiale, à l'acquisition des signaux sonores. Pour simuler un effet de champ proche de cette source lointaine, on applique le premier filtre de la relation
[A5] , comme indiqué dans la relation [A4] . La distance p représente alors une distance entre une source virtuelle proche M et le point O représentant 1 ' origine de la base sphérique de la figure 3. On applique ainsi un premier filtre de simulation de champ proche pour simuler la présence d'une source virtuelle à la distance p décrite ci-avant. Néanmoins, d'une part, comme indiqué ci-avant, les termes du coefficient de ce filtre divergent dans les basses fréquences (figure 6) et, d'autre part, la distance p précitée ne représentera pas forcément la distance entre les hauts-parleurs d'un dispositif de restitution et un point P de perception (figure 7). Selon l'invention, on applique une pré-compensation, à l'encodage, mettant en
jeu un filtre de type comme indiqué ci-avant, ce qui permet, d'une part, de transmettre des signaux bornés, et, d'autre part, de choisir la distance R, dès l'encodage, pour la restitution du son à partir des hauts- parleurs HPi, tel que représenté sur la figure 7. En particulier, on comprendra que si l'on a simulé, à l'acquisition, une source virtuelle placée à la distance p de l'origine O, à la restitution (figure 7), un auditeur placé au point P de perception auditive (à une distance R des hauts-parleurs HPi) ressentira, à l'audition, la présence d'une source sonore S, placée à la distance p du point de perception P et qui correspond à la source virtuelle simulée lors de l'acquisition.
Ainsi, la pré-compensation du champ proche des hauts- parleurs (placés à la distance R) , au stade de l'encodage, peut être combinée à un effet de champ proche simulé d'une source virtuelle placée à une distance p. A l'encodage, on met finalement en jeu un filtre total résultant, d'une part, de la simulation du champ proche, et, d'autre part, de la compensation du champ proche, les coefficients de ce filtre pouvant s'exprimer analytiquement par la relation :
Le filtre total donné par la relation [Ail] est stable et constitue la partie " encodage de distance " dans l'encodage ambisonique spatial selon 1 ' invention, tel que représenté sur la figure 8. Les coefficients de ces filtres correspondent à des fonctions de transfert monotones de la fréquence, qui tendent vers la valeur 1 en • hautes fréquences et vers la valeur (R/p)m en basses fréquences. En se référant à la figure 9, les spectres d'énergie des filtres H^FC^P' 'c'R l >ό)(ω) traduisent l'amplification des composantes encodées, dues à l'effet de champ de la source virtuelle (placée ici à une distance p = 1 m) , avec une pré-compensation du champ des hauts-parleurs (placés à une distance R = 1,5 m). L'amplification en décibels est donc positive lorsque p < R (cas de la figure 9) et négative quand p > R (cas de la figure 10 où p = 3 m et R = 1, 5 m) .
Dans un dispositif de restitution spatialisée, la distance
R entre un point de perception auditive et les haut- parleurs HPi est effectivement de l'ordre de un ou quelques mètres .
En se référant à nouveau à la figure 8, on comprendra que, outre les paramètres de direction θ et δ habituels, on transmettra une information sur les distances qui interviennent à l'encodage. Ainsi, les fonctions angulaires correspondant aux harmoniques sphériques Y^n sont conservées pour 1 ' encodage directionnel .
Toutefois, au sens de la présente invention, on prévoit en outre des filtres totaux (compensation de champ proche et, le cas échéant, simulation d'un champ proche) Hm ^p c' c'(co) qui sont appliqués aux composantes ambisoniques, en fonction de leur ordre m, pour réaliser l'encodage de la distance, comme représenté sur la figure 8. Un mode de réalisation de ces filtres dans le domaine audionumérique sera décrit en détail plus loin.
On remarquera, en particulier que ces filtres peuvent être appliqués dès même l'encodage de distance (r) et avant même l'encodage de direction (θ,δ) . On comprendra ainsi que les étapes a) et b) ci-avant peuvent être rassemblées en une même étape globale, ou même être interverties (avec un encodage de distance et filtrage de compensation, suivis d'un encodage de direction) . Le procédé selon 1 ' invention ne se limite donc pas à une mise en œuvre successive dans le temps des étapes a) et b) .
La figure 11A représente une visualisation (vue de dessus) d'une reconstruction d'un champ proche avec compensation, d'une onde sphérique, dans le plan horizontal (avec les mêmes paramètres de distance que ceux de la figure 9) , pour un système d'ordre total M = 15 et une restitution sur 32 hauts-parleurs . Sur la figure 11B, on a représenté la propagation de l'onde sonore initiale à partir d'une source en champ proche située à une distance p d'un point de l'espace d'acquisition qui correspond, dans l'espace de restitution, au point P de la figure 7 de perception auditive. On remarque sur la figure 11A que les auditeurs (symbolisés par des têtes schématisées) peuvent localiser la source virtuelle en un même lieu géographique situé à la distance p du point de perception P sur la figure 11B. On vérifie bien ainsi que la forme du front d'onde encodé est respectée après décodage et restitution. Toutefois, on constate sensiblement des interférences à droite du point P tel que représenté sur la figure 11A qui sont dues au fait que le nombre de hauts-parleurs (donc de composantes ambisoniques prises en compte) n'est pas suffisant pour restituer parfaitement le front d'ondes en jeu sur toute la surface délimitée par les haut-parleurs .
Dans ce qui suit, on décrit, à titre d'exemple, l'obtention d'un filtre audionumérique pour la mise en œuvre du procédé au sens de l'invention.
Comme indiqué ci-avant, si l'on cherche à simuler un effet de champ proche, compensé dès l'encodage, on applique aux composantes ambisoniques du son un filtre de la forme :
HNFC{Plc,R ,c) ω) [AII]
De l'expression de la simulation d'un champ proche donné par la relation [A5] , il apparaît que pour des sources lointaines (p = ∞) , la relation [Ail] devient simplement :
Il apparaît donc de cette' dernière relation [A12] que le cas où la source à simuler émet en champ lointain (source lointaine) n'est qu'un cas particulier de l'expression générale du filtre formulée dans la relation [Ail] .
Dans le domaine des traitements audionumériques, un procédé avantageux pour définir un filtre numérique à partir de l'expression analytique de ce filtre dans le domaine analogique à temps continu consiste en une
"transformée bilinéaire " .
On exprime d'abord la relation [A5] sous la forme d'une transformée de Laplace, ce qui correspond à :
où τ = p/c (c étant la vitesse acoustique dans le milieu, typiquement 340 m/s dans l'air).
La transformée bilinéaire consiste à présenter, pour une fréquence d'échantillonnage fs, la relation [Ail] sous la forme :
f-2 _•« H--,"*-1 +έfz"2 bXl) l1 +-,<*"+1>2z-1
«•« - π q= α °^ + α '^,z-, + α yz 2 - ÛQ D '/2 , + α^ o '/2 z-ι [AI4]
si m est impair et
si m est pair,
où z est défini par p = 2fs par rapport à la relation
[A13] précédente,
et avec :
et
Re(Xm,q X m,q x2 = 1 +2 + a a'
où α = 4fs R/c pour x=a et α = 4fs p/c pour x=b
Xm7q sont les q racines successives du polynôme de Bessel
q=l et sont exprimées dans le tableau 1 ci-après, pour différents ordres m, sous les formes respectives de leur partie réelle, leur module (séparés par une virgule) et leur valeur (réelle) lorsque m est impair.
Tableau 1 : valeurs Pe X m,q\ < Xm,q (et Re\.Xm,m\ lorsque m est impair) d'un polynôme de Bessel calculées à l'aide du logiciel de calcul MATLAB©.
On implémente ainsi les filtres numériques, à partir des valeurs du tableau 1, en prévoyant des cascades de cellules d'ordre 2 (pour m pair), et une cellule supplémentaire (pour m impair) , à partir des relations [A14] données ci-avant.
On réalise ainsi des filtres numériques sous une forme de réponse impulsionnelle infinie, aisément paramétrable comme montré ci-avant. Il est à noter qu'une implémentation sous une forme de réponse impulsionnelle finie peut être envisagée et consiste à calculer le spectre complexe de la fonction de transfert à partir de la formule analytique, puis à en déduire une réponse impulsionnelle finie par transformée de Fourier inverse. On applique ensuite une opération de convolution pour le filtrage.
Ainsi, en introduisant cette pré-compensation du champ proche à l'encodage, on définit une représentation ambisonique modifiée (figure 8) , en adoptant comme représentation transmissible des signaux exprimés dans le domaine fréquentiel, sous la forme :
βσ(R / c) = 1 Bσ **A151
F mR l c(ω)
Comme indiqué ci-avant, R est une distance de référence à laquelle est associé un effet de champ proche compensé et c est la vitesse du son (typiquement 340 m/s dans l'air) . Cette représentation . ambisonique modifiée possède les mêmes propriétés de scalabilité (schématiquement représentée par des données transmises "entourées" près de la flèche TR de la figure 1) et obéit aux mêmes transformations de rotation du champ (module 4 de la figure 1) que la représentation ambisonique habituelle.
On indique ci-après les opérations à mettre en œuvre pour le décodage des signaux ambisoniques reçus.
On indique tout d'abord que l'opération de décodage est adaptable à un dispositif de restitution quelconque, de rayon R2, différent de la distance de référence R ci- avant. A cet effet, on applique des filtres de type H (ω) ι tels que décrits plus haut, mais avec des paramètres de distance R et R2, au lieu de p et R. En particulier, il est à noter que seul le paramètre R/c est à mémoriser (et/ou transmettre) entre l'encodage et le décodage .
En se référant à la figure 12, le module de filtrage qui y est représenté est prévu par exemple dans une unité de traitement d'un dispositif de restitution. Les composantes ambisoniques reçues ont été pré-compensées à l'encodage pour une distance de référence Ri en tant que seconde distance. Toutefois, le dispositif de restitution comporte une pluralité de hauts-parleurs disposés à une troisième distance R2 d'un point de perception auditive P, cette troisième distance R2 étant différente de la seconde distance précitée R_ . Le module de filtrage de la figure
12, sous la forme H^FC^ ' C'R2 ' c) (ω) , adapte alors, à la réception des données, la pré-compensation à la distance Rx pour une restitution à la distance R2. Bien entendu, comme indiqué ci-avant, le dispositif de restitution reçoit aussi le paramètre Ri/c .
Il est à noter que 1 ' invention permet en outre de mixer plusieurs représentations ambisoniques de champs sonores (sources réelles et/ou virtuelles) , dont les distances de référence R sont différentes (le cas échéant avec des distances de référence infinies et correspondant à des sources lointaines) . Préférentiellement, on filtrera une pré-compensation de toutes ces sources à une distance de référence la plus petite, avant de mélanger les signaux ambisoniques, ce qui permet à la restitution d'obtenir une définition correcte du relief sonore.
Dans le cadre d'un traitement dit de "focalisation sonore " avec, à la restitution, un effet d'enrichissement sonore pour une direction choisie de l'espace (à la manière d'un projecteur lumineux éclairant dans une direction choisie en optique) , impliquant un traitement matriciel de focalisation sonore (avec pondération des composantes ambisoniques), on applique avantageusement l'encodage de distance avec pré-compensation de champ proche de façon combinée au traitement de focalisation.
Dans ce qui suit, on décrit un procédé de décodage ambisonique, avec compensation du champ proche des hauts- parleurs, à la restitution.
Pour reconstruire un champ acoustique encodé suivant le formalisme ambisonique, à partir des composantes B^n et en utilisant des hauts-parleurs d'un dispositif de restitution qui prévoit un emplacement "idéal " d'un auditeur qui correspond au point de restitution P de la figure 7, l'onde émise par chaque haut-parleur est définie par un traitement préalable de "ré-encodage " du champ ambisonique au centre du dispositif de restitution, comme suit.
Dans ce contexte de " ré- encodage " , on considère dans un premier temps et pour simplification que les sources émettent en champ lointain. En se référant à nouveau à la figure 7, l'onde émise par un haut-parleur d'indice i et d'incidence (θi et δi) est alimenté par un signal Si. Ce haut-parleur participe à la
reconstruction de la composante Bmn , par sa contribution
Le vecteur Ci des coefficients d'encodage associés aux hauts-parleurs d'indice i s'exprime par la relation :
*00 (#ι A) ,- = A) [Bl]
Le vecteur S des signaux émanant de l'ensemble des N hauts-parleurs est donné par l'expression :
La matrice d'encodage de ces N hauts-parleurs (qui correspond finalement à une matrice de "ré -encodage ") , s ' exprime par la relation :
C= [c_ C2 ...CN] [B3] où chaque terme ci représente un vecteur selon la relation [Bl] ci-avant.
Ainsi, la reconstruction du champ ambisonique B' est définie par la relation :
La relation [B4] définit ainsi une opération de ré- encodage, préalable à la restitution. Finalement, le décodage, en tant que tel, consiste à comparer les signaux ambisoniques originaux et reçus par le dispositif de restitution, sous la forme :
^00 R+1
B = B 11 [B5]
B mn
au signaux ré-encodés B , pour définir la relation générale :
B' = B [B6] Il s'agit, en particulier, de déterminer les coefficients d'une matrice de décodage D, qui vérifie la relation :
S ≈ D.B [B7]
De préférence, le nombre de hauts-parleurs est supérieur ou égal au nombre de composantes ambisoniques à décoder et la matrice de décodage D s'exprime, en fonction de la matrice de ré-encodage C, sous la forme :
D = CT . (C.CT)~ [B8]
où la notation Cτ correspond à la transposée de la matrice C.
Il est à noter que la définition d'un décodage vérifiant des critères différents par bandes de fréquences est possible, ce qui permet d'offrir une restitution optimisée en fonction des conditions d'écoute, notamment pour ce qui concerne la contrainte de positionnement au centre O de la sphère de la figure 3, lors de la restitution. A cet effet, on prévoit avantageusement un filtrage simple, en égalisation fréquentielle par paliers, à chaque composante ambisonique .
Toutefois, pour obtenir une reconstruction d'une onde originellement encodée, il faut corriger l'hypothèse de champ lointain pour les hauts-parleurs, c'est-à-dire exprimer 1 ' effet de leur champ proche dans la matrice de ré-encodage C ci-avant et inverser ce nouveau système pour définir le décodeur. A cet effet, en supposant une concentricité des hauts-parleurs (disposés à une même distance R du point P de la figure 7) , tous les hauts-
parleurs ont un même effet de champ proche F^ ' (ω) , sur
chaque composante ambisonique du type R'^„ . En introduisant les termes de champ proche sous la forme d'une matrice diagonale, la relation [B4] ci-avant devient :
B'= Diag[ FR lc(ω) FR l c(ω) F mR l c(ω) F mR l c(ω) ]) . C.S [B9]
La relation [B7] ci-avant devient
Ainsi, l'opération de matriçage est précédée par une opération de filtrage qui compense le champ proche sur chaque composante B mn ι et qui peut être mise en œuvre sous forme numérique, comme décrit ci-avant, en référence à la relation [A14] .
On retiendra qu'en pratique, la matrice C de "ré-encodage" est propre au dispositif de restitution. Ses coefficients peuvent être déterminés initialement par paramétrage et caracterisation sonore du dispositif de restitution réagissant à un excitation prédéterminée. La matrice de décodage D est, elle aussi, propre au dispositif de restitution. Ses coefficients peuvent être déterminés par la relation [B8] . En reprenant la notation précédente où B est la matrice des composantes ambisoniques précompensées, ces dernières peuvent être transmises au dispositif de restitution sous forme matricielle B avec :
Le dispositif de restitution décode ensuite les données reçues sous forme matricielle B (vecteur colonne des composantes transmises) en appliquant la matrice de décodage D aux composantes ambisoniques pré-compensées, pour former les signaux Si destinés à alimenter les haut- parleurs HPi, avec :
En se référant à nouveau à la figure 12, si une opération de décodage doit être adaptée à un dispositif de restitution de rayon R2 différent de la distance de référence R*-., un module d'adaptation préalable au décodage proprement dit et décrit ci-avant permet de filtrer chaque composante ambisonique Bmn , pour l'adapter à un dispositif de restitution de rayon R2. L'opération de décodage proprement dite est effectuée ensuite, comme décrit ci-avant, en référence à la relation [Bll] . On décrit ci-après une application de l'invention à la synthèse binaurale.
On se réfère à la figure 13A sur laquelle un auditeur disposant d'un casque à deux écouteurs d'un dispositif de synthèse binaurale est représenté. Les deux oreilles de 1 ' auditeur sont disposées à des points respectifs 0L (oreille gauche) et 0R (oreille droite) de l'espace. Le centre de la tête de l'auditeur est disposé au point O et le rayon de la tête de l'auditeur est de valeur a. Une source sonore doit être perçue auditivement à un point M de l'espace, situé à une distance r du centre de la tête de l'auditeur (et respectivement à des distances rR de l'oreille droite et rL de l'oreille gauche). Par ailleurs, la direction de la source placée au point M est définie par les vecteurs F , FR et fj_ .
De façon générale, la synthèse binaurale se définit comme suit.
Chaque auditeur a une forme d'oreille qui lui est propre. La perception d'un son dans l'espace par cet auditeur se fait par apprentissage, depuis la naissance, en fonction de la forme des oreilles (notamment la forme des pavillons et les dimensions de la tête) propre à cet auditeur. La perception d'un son dans l'espace se manifeste entre autres par le fait que le son parvient à une oreille, avant l'autre oreille, ce qui se traduit par un retard τ entre les signaux à émettre par chaque écouteur du dispositif de restitution appliquant la synthèse binaurale . Le dispositif de restitution est paramétré initialement, pour un même auditeur, en balayant une source sonore autour de sa tête, à une même distance R du centre de sa tête. On comprendra ainsi que cette distance R peut être considérée comme une distance entre un "point de restitution" comme énoncé ci-avant et un point de perception auditive (ici le centre O de la tête de 1 'auditeur) .
Dans ce qui suit, l'indice L est associé au signal à restituer par l'écouteur accolé à l'oreille gauche et 1 ' indice R est associé au signal à restituer par l'écouteur accolé à l'oreille droite. En se référant à la figure 13B, on applique au signal initial S un retard pour chaque voie destinée à produire un signal pour un écouteur distinct. Ces retards τL et τR sont fonction d-'un retard maximum XMAX qui correspond ici au rapport a/c où a, comme indiqué précédemment, correspond au rayon de la tête de l'auditeur et c à la vitesse du son. En particulier, ces retards sont définis en fonction de la différence de distance du point O (centre de la tête) au point M (position de la source dont le son est à restituer, sur la figure 13A) et de chaque oreille à ce point M. Avantageusement, on applique en outre des gains respectifs S et -R/ à chaque voie, qui sont fonction d'un rapport des distances du point O au point M et de chaque oreille au point M. Des modules respectifs appliqués à chaque voie 2L et 2R encodent les signaux de chaque voie, dans une représentation ambisonique, avec pré-compensation de champ proche NFC (pour "jVear Field Compensation" ) au sens de la présente invention. On comprendra ainsi que, par la mise en œuvre du procédé au sens de la présente invention, on peut définir les signaux issus de la source M, non seulement par leur direction (angles azimutaux Θ et ΘR et angles d'élévation δL et δR) , mais aussi en fonction de la distance séparant chaque oreille rL et rR de la source M. Les signaux ainsi encodés sont transmis au dispositif de restitution comportant des modules de décodage ambisonique, pour chaque voie, 5L et 5R. Ainsi, on applique un encodage/décodage ambisonique, avec compensation de champ proche, pour chaque voie (écouteur gauche, écouteur droit) dans la restitution avec synthèse binaurale (ici de type "B-FORMAT") , sous forme dédoublée. La compensation de champ proche s'effectue, pour chaque voie, avec comme première distance p une distance rL et rR entre chaque oreille et la position M de la source sonore à restituer.
On décrit ci-après une application de la compensation au sens de l'invention, au contexte de l'acquisition sonore en représentation ambisonique.
On se réfère à la figure 14 sur laquelle un microphone 141 comporte une pluralité de capsules transductrices, capables de capter des pressions acoustiques et restituer des signaux électriques Si, ...,SN. Les capsules CAPi sont agencées sur une sphère de rayon r prédéterminé (ici, une sphère rigide, telle qu'une balle de ping-pong par exemple) . Les capsules sont espacées d'un pas régulier sur la sphère. En pratique, on choisit le nombre N de capsules en fonction de l'ordre M désiré pour la .représentation ambisonique.
On indique ci-après, dans le contexte d'un microphone comportant des capsules agencées sur une sphère rigide, comment compenser l'effet de champ proche, dès l'encodage dans le contexte ambisonique. On montrera ainsi que la pré-compensation du champ proche peut s'appliquer non seulement pour la simulation de source virtuelle, comme indiqué ci-avant, mais aussi à l'acquisition et, de façon plus générale, en combinant la pré-compensation de champ proche à tous types de traitements impliquant une représentation ambisonique .
En présence d'une sphère rigide (susceptible d'introduire une diffraction des ondes sonores reçues) , la relation [Al] donnée, ci-avant devient :
pr& = ("t) [ci]
Les dérivées des fonctions de Hankel sphériques h" obéissent à la loi de récurrence :
(2m+ ï)h~ (x) = m h~_ι (x) - (m +1) A~ +1 (x) [C2]
On déduit les composantes ambisoniques R^κ du champ initial à partir du champ de pression à la surface de la sphère, en mettant en œuvre des opérations de projection et d'égalisation données par la relation :
Dans cette expression, EQra est un filtre égaliseur qui compense une pondération m qui est liée à la directivité des capsules et qui inclut en outre la diffraction par la sphère rigide.
L'expression de ce filtre EQm est donnée par la relation suivante :
EQm = (kr)2h- (kr)rm+1 [C4] w m,
Les coefficients de ce filtre d'égalisation ne sont pas stables et on obtient un gain infini en très basses fréquences. D'ailleurs, il convient de noter que les composantes harmoniques sphériques, elles-mêmes, ne sont pas d'amplitude finie lorsque le champ sonore n'est pas limité à une propagation d'ondes planes, c'est-à-dire issues de sources lointaines, comme on l'a vu précédemment .
Par ailleurs, si, plutôt que de prévoir des capsules encastrées dans une sphère solide, on prévoit des capsules de type cardioïdes, avec une directivité en champ lointain donnée par l'expression : G(θ) = + (l -a)cosθ [C5]
En considérant ces capsules montées sur un support "transparent acoustiquement " , le terme de pondération à compenser devient :
Wm = jm {c jm(kr) - j(l -α)jm kr)) [C6]
Il apparaît encore que les coefficients d'un filtre d'égalisation correspondant à l'inverse analytique de cette pondération donnée par la relation [C6] sont divergents pour les très basses fréquences.
De façon générale, on indique que pour tout type de directivité de capteurs, le gain du filtre EQm pour compenser la pondération m liée à la directivité des capteurs est infini pour les basses fréquences sonores. En se référant à la figure 14, on applique avantageusement une pré-compensation de champ proche dans l'expression même du filtre d'égalisation EQm, donnée par la relation :
Ainsi, les signaux Si à SN sont récupérés du microphone 141. Le cas échéant, on applique une pré-égalisation de ces signaux par un module de traitement 142. Le module 143 permet d'exprimer ces signaux dans le contexte ambisonique, sous forme matricielle. Le module 144 applique le filtre de la relation [C7] aux composantes ambisoniques exprimées en fonction du rayon r de la sphère du microphone 141. La compensation de champ proche s'effectue pour une distance de référence R en tant que seconde distance. Les signaux encodés et ainsi filtrés par le module 144 peuvent être transmis, le cas échéant, avec le paramètre représentatif de la distance de référence R/c.
Ainsi, il apparaît dans les différents modes de réalisation liés respectivement à la création d'une source virtuelle en champ proche, à l'acquisition de signaux sonores issues de sources réelles, ou même à la restitution (pour compenser un effet de champ proche des hauts-parleurs) , que la compensation de champ proche au sens de la présente invention peut s'appliquer à tous types de traitements faisant intervenir une représentation ambisonique. Cette compensation de champ proche permet d'appliquer la représentation ambisonique à une multiplicité de contextes sonores où la direction d'une source et avantageusement sa distance doivent être prises en compte. De plus, la possibilité de la représentation de phénomènes sonores de tous types (champs proches ou lointains) dans le contexte ambisonique est assurée par cette pré-compensation, du fait de la limitation à des valeurs réelles finies des composantes ambisoniques.
Bien entendu, la présente invention ne se limite pas à la forme de réalisation décrite ci-avant à titre d'exemple ; elle s'étend à d'autres variantes. Ainsi, on comprendra que la pré-compensation de champ proche peut être intégrée, à l'encodage, autant pour une source proche que pour une source lointaine. Dans ce dernier cas (source lointaine et réception d'ondes planes) , la distance p exprimée ci-avant sera considérée comme infinie, sans modifier de façon substantielle l'expression des filtres Hm donnée ci-avant. Ainsi, le traitement utilisant des processeurs d'effet de salle qui fournissent en général des signaux décorrélés utilisables pour modeliser le champ diffus tardif (réverbération tardive) peut être combiné à une pré-compensâtion de champ proche. On peut considérer que ces signaux sont de même énergie et correspondent à une part de champ diffus correspondant à la composante omnidirective W = BQQ (figure 4) . On peut alors construire les diverses composantes harmoniques sphériques (avec un ordre M choisi) en appliquant une correction de gain pour chaque composante ambisonique et on applique une compensation de champ proche des hauts-parleurs (avec une distance de référence R séparant les haut-parleurs du point de perception auditive comme représenté sur la figure 7) .
Bien entendu, le principe d'encodage au sens de la présente invention est generalisable à des modèles de rayonnement autres que des sources monopolaires (réelles ou virtuelles) et/ou des hauts-parleurs . En effet, toute forme de rayonnement (notamment une source étalée dans l'espace) peut être exprimée par intégration d'une distribution continue de sources élémentaires ponctuelles. En outre, dans le contexte de la restitution, il est possible d'adapter la compensation de champ proche à tout contexte de restitution. A cet effet, il peut être prévu de calculer des fonctions de transfert (ré-encodage des composantes harmoniques sphériques de champ proche pour chaque haut-parleur, compte tenu d'une propagation* réelle dans la salle où le son est restitué) , ainsi qu'une inversion de ce ré-encodage pour redéfinir le décodage.
On a décrit ci-avant un procédé de décodage dans lequel on appliquait un système matriciel faisant intervenir les composantes ambisoniques. Dans une variante, il peut être prévu un traitement généralisé par transformées de Fourier rapides (circulaire ou sphérique) pour limiter les temps de calcul et les ressources informatiques (en terme de mémoire) nécessaires au traitement de décodage.
Comme indiqué ci-avant en référence aux figures 9 et 10, on constate que le choix d'une distance de référence R par rapport à la distance p de la source en champ proche introduit une différence de gain pour différentes valeurs de la fréquence sonore . On indique que le procédé d'encodage avec pré-compensation peut être couplé à une compression audionumérique permettant de quantifier et d'ajuster le gain pour chaque sous-bande fréquentielle.
Avantageusement, la présente invention s'applique à tous types de systèmes de spatialisation sonore, notamment pour des applications de type "réalité virtuelle " (navigation dans des scènes virtuelles dans l'espace tridimensionnel, jeux avec spatialisation sonore tridimensionnelle, conversations de type "chat " sonorisées sur le réseau Internet), à des sonifications d'interfaces, à des logiciels d'édition audio pour enregistrer, mixer et restituer de la musique, mais aussi à l'acquisition, à partir d'usage de microphones tridimensionnels, pour la prise de son musicale ou cinématographique, ou encore pour la transmission d'ambiance sonore sur Internet, par exemple pour des "WebCam" sonorisées.

Claims

REVENDICATIONS
1. Procédé de traitement de données sonores, dans lequel : a) on code des signaux représentatifs d'au moins un son se propageant dans l'espace tridimensionnel et issu d'une source située à une première distance (p) d'un point de référence (O) , pour obtenir une représentation du son par des composantes (Bran σ) exprimées dans une base d'harmoniques sphériques, d'origine correspondant audit point de référence (O) , b) et on applique auxdites composantes (B^) une compensation d'un effet de champ proche par un filtrage qui est fonction d'une seconde distance (R) définissant sensiblement, pour une restitution du son par un dispositif de restitution, une distance entre un point de restitution (HPi) et un point (P) de perception auditive.
2. Procédé selon la revendication 1, dans lequel, ladite source étant lointaine du point de référence (O) , - on obtient des composantes d'ordres successifs m pour la représentation du son dans ladite base d'harmoniques sphériques , et
- on applique un filtre (l/Fm) dont les coefficients, appliqués chacun à une composante d'ordre m, s'expriment analytiquement sous la forme de l'inverse d'un polynôme de puissance m, dont la variable est inversement proportionnelle à la fréquence sonore et à ladite seconde distance (R) , pour compenser un effet de champ proche au niveau du dispositif de restitution.
3. Procédé selon la revendication 1, dans lequel, ladite source étant une source virtuelle prévue à ladite première distance (p) ,
- on obtient des composantes d'ordres successifs m pour la représentation du son dans ladite base d'harmoniques sphériques , et
- on applique un filtre global (Hm) dont les coefficients, appliqués chacun à une composante d'ordre m, s'expriment analytiquement sous la forme d'une fraction, dont : - le numérateur est un polynôme de puissance m, dont la variable est inversement proportionnelle à la fréquence sonore et à ladite première distance (p) , pour simuler un effet de champ proche de la source virtuelle, et - le dénominateur est un polynôme de puissance m, dont la variable est inversement proportionnelle à la fréquence sonore et à ladite seconde distance (R) , pour compenser 1 ' effet du champ proche de la source virtuelle dans les basses fréquences sonores.
4. Procédé selon l'une des revendications précédentes, dans lequel on transmet au dispositif de restitution les données codées et filtrées aux étapes a) et b) avec un paramètre représentatif de ladite seconde distance (R/c) .
5. Procédé selon l'une des revendications 1 à 3, dans lequel, le dispositif de restitution comportant des moyens de lecture d'un support mémoire, on mémorise sur un support mémoire destiné à être lu par le dispositif de restitution les données codées et filtrées aux étapes a) et b) avec un paramètre représentatif de ladite seconde distance (R/c) .
6. Procédé selon l'une des revendications 4 et 5, dans lequel, préalablement à une restitution sonore par un dispositif de restitution comportant une pluralité de hauts-parleurs disposés à une troisième distance (R2) dudit point de perception auditive (P) , on applique aux données codées et filtrées un filtre d'adaptation (Hm (Ri/c,R2/c) -, ς ont ιes coefficients sont fonction desdites seconde (Ri) et troisième distances (R2) .
7. Procédé selon la revendication 6, dans lequel les coefficients dudit filtre d'adaptation (Hm <R1/c'R2/c) ) , appliqués chacun à une composante d'ordre m, s'expriment analytiquement sous la forme d'une fraction, dont :
- le numérateur est un polynôme de puissance m, dont la variable est inversement proportionnelle à la fréquence sonore et à ladite seconde distance (R) , - et le dénominateur est un polynôme de puissance m, dont la variable est inversement proportionnelle à la fréquence sonore et à ladite troisième distance (R2) .
8. Procédé selon l'une des revendications 2, 3 et 7, dans lequel, pour la mise en œuvre de l'étape b) , on prévoit :
- pour des composantes d'ordre m pair, des filtres audionumériques sous la forme d'une cascade de cellules d'ordre deux ; et
- pour des composantes d'ordre m impair, des filtres audionumériques sous la forme d'une cascade de cellules d'ordre deux et une cellule supplémentaire d'ordre un.
9. Procédé selon la revendication 8 , dans lequel les coefficients d'un filtre audionumérique, pour une composante d'ordre m, sont définis à partir des valeurs numériques des racines desdits polynômes de puissance m.
10. Procédé selon l'une des revendications 2, 3, 7, 8 et 9, dans lequel lesdits polynômes sont des polynômes de Bessel.
11. Procédé selon l'une des revendications 1, 2 et 4 à 10, dans lequel on prévoit un microphone comportant un réseau de transducteurs acoustiques agencés sensiblement sur la surface d'une sphère dont le centre correspond sensiblement audit point de référence (O) , pour obtenir lesdits signaux représentatifs d'au moins un son se propageant dans l'espace tridimensionnel.
12. Procédé selon la revendication 11, dans lequel on applique à l'étape b) un filtre global pour, d'une part, compenser un effet de champ proche en fonction de ladite seconde distance (R) et, d'autre part, égaliser les signaux issus des transducteurs pour compenser une pondération de directivité desdits transducteurs.
13. Procédé selon l'une des revendications 11 et 12, dans lequel on prévoit un nombre de transducteurs fonction d'un nombre total choisi de composantes pour représenter le son dans ladite base d'harmoniques sphériques.
14. Procédé selon l'une des revendications précédentes, dans lequel on choisit à l'étape a) un nombre total de composantes dans la base des harmoniques sphériques pour obtenir, à la restitution, une région de l'espace autour du point de perception (P) dans laquelle la restitution du son est fidèle et dont les dimensions sont croissantes avec le nombre total de composantes .
15. Procédé selon la revendication 14 , dans lequel on prévoit un dispositif de restitution comportant un nombre de haut-parleurs au moins égal audit nombre total de composantes .
16. Procédé selon l'une des revendications 1 à 5 et 8 à 13, dans lequel :
- on prévoit un dispositif de restitution comportant au moins un premier et un second haut-parleur disposés à une distance choisie d'un auditeur, on obtient, pour cet auditeur, une information de ressenti de la position dans l'espace de sources sonores situées à une distance de référence prédéterminée (R) de l'auditeur, et
- on applique la compensation de l'étape b) avec ladite distance de référence sensiblement en tant que seconde distance.
17. Procédé selon l'une des revendications 1 à 3 et 8 à 13, prises en combinaison avec l'une des revendications 4 et 5, dans lequel : - on prévoit un dispositif de restitution comportant au moins un premier et un second haut-parleur disposés à une distance choisie d'un auditeur, on obtient, pour cet auditeur, une information de ressenti de la position dans l'espace de sources sonores situées à une distance de référence prédéterminée (R2) de l'auditeur, et
- préalablement à une restitution sonore par le dispositif de restitution, on applique aux données codées et filtrées aux étapes a) et b) un filtre d'adaptation (Hm (R/c'R2 c) ) dont les coefficients sont fonction de la seconde distance (R) et sensiblement de la distance de référence (R2) .
18. Procédé selon l'une des revendications 16 et 17, dans lequel :
- le dispositif de restitution comporte un casque à deux écouteurs pour les oreilles respectives de l'auditeur, et
- séparément pour chaque écouteur, on applique le codage et le filtrage des étapes a) et b) pour des signaux respectifs destinés à alimenter chaque écouteur, avec, en tant que première distance (p) , respectivement une distance (rR,r-**.) séparant chaque oreille d'une position (M) d'une source à restituer.
19. Procédé selon l'une des revendications précédentes, dans lequel on met en forme, aux étapes a) et b) , un système matriciel comportant au moins :
- une matrice (B) comportant lesdites composantes dans la base des harmoniques sphériques, et - une matrice diagonale (Diag(l/Fm) ) dont les coefficients correspondent à des coefficients de filtrage de l'étape b), et on multiplie lesdites matrices pour obtenir une matrice résultat de composantes compensées ( B ) .
20. Procédé selon la revendication 19, dans lequel :
- le dispositif de restitution comporte une pluralité de haut-parleurs disposés sensiblement à une même distance (R) du point de perception auditive (P) , et
- pour décoder lesdites données codées et filtrées aux étapes a) et b) et former des signaux adaptés pour alimenter lesdits haut-parleurs :
* on forme un système matriciel comportant ladite matrice résultat ( B ) et une matrice de décodage (D) prédéterminée, propre au dispositif de restitution, et
* on obtient une matrice (S) comportant des coefficients représentatifs des signaux d'alimentation des hauts-parleurs par multiplication de la matrice des composantes compensées ( B ) par ladite matrice de décodage (D) .
21. Dispositif d'acquisition sonore, comportant un microphone muni d'un réseau de transducteurs acoustiques disposés sensiblement sur la surface d'une sphère, caractérisé en ce qu'il comporte en outre une unité de traitement agencée pour :
- recevoir des signaux émanant chacun d'un transducteur, - appliquer auxdits signaux un codage pour obtenir une représentation du son par des composantes {B^) exprimées dans une base d'harmoniques sphériques, d'origine correspondant au centre de ladite sphère (O) , - et appliquer auxdites composantes (Bπuι σ) un filtrage qui est fonction, d'une part, d'une distance correspondant au rayon de la sphère (r) et, d'autre part, d'une distance de référence (R) .
22. Dispositif selon la revendication 21, caractérisé en ce que ledit filtrage consiste, d'une part, à égaliser, en fonction du rayon de la sphère, les signaux issus des transducteurs pour compenser une pondération de directivité desdits transducteurs et, d'autre part, à compenser un effet de champ proche en fonction d'une distance de référence choisie (R) , définissant sensiblement, pour une restitution du son, une distance entre un point de restitution (HPi) et un point (P) de perception auditive .
EP03782553A 2002-11-19 2003-11-13 Procede de traitement de donnees sonores et dispositif d'acquisition sonore mettant en oeuvre ce procede Expired - Lifetime EP1563485B1 (fr)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR0214444A FR2847376B1 (fr) 2002-11-19 2002-11-19 Procede de traitement de donnees sonores et dispositif d'acquisition sonore mettant en oeuvre ce procede
FR0214444 2002-11-19
PCT/FR2003/003367 WO2004049299A1 (fr) 2002-11-19 2003-11-13 Procede de traitement de donnees sonores et dispositif d'acquisition sonore mettant en oeuvre ce procede

Publications (2)

Publication Number Publication Date
EP1563485A1 true EP1563485A1 (fr) 2005-08-17
EP1563485B1 EP1563485B1 (fr) 2006-03-29

Family

ID=32187712

Family Applications (1)

Application Number Title Priority Date Filing Date
EP03782553A Expired - Lifetime EP1563485B1 (fr) 2002-11-19 2003-11-13 Procede de traitement de donnees sonores et dispositif d'acquisition sonore mettant en oeuvre ce procede

Country Status (13)

Country Link
US (1) US7706543B2 (fr)
EP (1) EP1563485B1 (fr)
JP (1) JP4343845B2 (fr)
KR (1) KR100964353B1 (fr)
CN (1) CN1735922B (fr)
AT (1) ATE322065T1 (fr)
AU (1) AU2003290190A1 (fr)
BR (1) BRPI0316718B1 (fr)
DE (1) DE60304358T2 (fr)
ES (1) ES2261994T3 (fr)
FR (1) FR2847376B1 (fr)
WO (1) WO2004049299A1 (fr)
ZA (1) ZA200503969B (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9460726B2 (en) 2010-03-26 2016-10-04 Dolby Laboratories Licensing Corporation Method and device for decoding an audio soundfield representation for audio playback
US10582329B2 (en) 2016-01-08 2020-03-03 Sony Corporation Audio processing device and method

Families Citing this family (76)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10328335B4 (de) * 2003-06-24 2005-07-21 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Wellenfeldsyntesevorrichtung und Verfahren zum Treiben eines Arrays von Lautsprechern
US20050271216A1 (en) * 2004-06-04 2005-12-08 Khosrow Lashkari Method and apparatus for loudspeaker equalization
KR101304797B1 (ko) * 2005-09-13 2013-09-05 디티에스 엘엘씨 오디오 처리 시스템 및 방법
US8059824B2 (en) * 2006-03-13 2011-11-15 France Telecom Joint sound synthesis and spatialization
FR2899424A1 (fr) * 2006-03-28 2007-10-05 France Telecom Procede de synthese binaurale prenant en compte un effet de salle
US8180067B2 (en) * 2006-04-28 2012-05-15 Harman International Industries, Incorporated System for selectively extracting components of an audio input signal
US7876903B2 (en) * 2006-07-07 2011-01-25 Harris Corporation Method and apparatus for creating a multi-dimensional communication space for use in a binaural audio system
US8036767B2 (en) 2006-09-20 2011-10-11 Harman International Industries, Incorporated System for extracting and changing the reverberant content of an audio input signal
DE602007011955D1 (de) * 2006-09-25 2011-02-24 Dolby Lab Licensing Corp Ür mehrkanal-tonwiedergabesysteme mittels ableitung von signalen mit winkelgrössen hoher ordnung
DE102006053919A1 (de) * 2006-10-11 2008-04-17 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Vorrichtung und Verfahren zum Erzeugen einer Anzahl von Lautsprechersignalen für ein Lautsprecher-Array, das einen Wiedergaberaum definiert
JP2008118559A (ja) * 2006-11-07 2008-05-22 Advanced Telecommunication Research Institute International ３次元音場再生装置
JP4873316B2 (ja) * 2007-03-09 2012-02-08 株式会社国際電気通信基礎技術研究所 音響空間共有装置
EP2094032A1 (fr) * 2008-02-19 2009-08-26 Deutsche Thomson OHG Signal audio, procédé et appareil pour coder ou transmettre celui-ci et procédé et appareil pour le traiter
EP2250821A1 (fr) * 2008-03-03 2010-11-17 Nokia Corporation Appareil de capture et de rendu d'une pluralité de canaux audio
EP2154910A1 (fr) * 2008-08-13 2010-02-17 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Appareil de fusion de flux audio spatiaux
ES2425814T3 (es) 2008-08-13 2013-10-17 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Aparato para determinar una señal de audio espacial convertida
GB0815362D0 (en) 2008-08-22 2008-10-01 Queen Mary & Westfield College Music collection navigation
US8819554B2 (en) * 2008-12-23 2014-08-26 At&T Intellectual Property I, L.P. System and method for playing media
EP2205007B1 (fr) * 2008-12-30 2019-01-09 Dolby International AB Procédé et appareil pour le codage tridimensionnel de champ acoustique et la reconstruction optimale
GB2467534B (en) 2009-02-04 2014-12-24 Richard Furse Sound system
CN102318373B (zh) 2009-03-26 2014-09-10 松下电器产业株式会社 解码装置、编解码装置及解码方法
KR20140010468A (ko) * 2009-10-05 2014-01-24 하만인터내셔날인더스트리스인코포레이티드 오디오 신호의 공간 추출 시스템
JP5672741B2 (ja) * 2010-03-31 2015-02-18 ソニー株式会社 信号処理装置および方法、並びにプログラム
US20110317522A1 (en) * 2010-06-28 2011-12-29 Microsoft Corporation Sound source localization based on reflections and room estimation
US9313599B2 (en) 2010-11-19 2016-04-12 Nokia Technologies Oy Apparatus and method for multi-channel signal playback
US9055371B2 (en) * 2010-11-19 2015-06-09 Nokia Technologies Oy Controllable playback system offering hierarchical playback options
US9456289B2 (en) 2010-11-19 2016-09-27 Nokia Technologies Oy Converting multi-microphone captured signals to shifted signals useful for binaural signal processing and use thereof
EP2541547A1 (fr) 2011-06-30 2013-01-02 Thomson Licensing Procédé et appareil pour modifier les positions relatives d'objets de son contenu dans une représentation ambisonique d'ordre supérieur
US9338572B2 (en) * 2011-11-10 2016-05-10 Etienne Corteel Method for practical implementation of sound field reproduction based on surface integrals in three dimensions
KR101282673B1 (ko) 2011-12-09 2013-07-05 현대자동차주식회사 음원 위치 추정 방법
US8996296B2 (en) * 2011-12-15 2015-03-31 Qualcomm Incorporated Navigational soundscaping
CN104137248B (zh) 2012-02-29 2017-03-22 应用材料公司 配置中的除污及剥除处理腔室
EP2645748A1 (fr) * 2012-03-28 2013-10-02 Thomson Licensing Procédé et appareil de décodage de signaux de haut-parleurs stéréo provenant d'un signal audio ambiophonique d'ordre supérieur
WO2013150341A1 (fr) 2012-04-05 2013-10-10 Nokia Corporation Appareil de capture d'élément audio spatial flexible
US9288603B2 (en) * 2012-07-15 2016-03-15 Qualcomm Incorporated Systems, methods, apparatus, and computer-readable media for backward-compatible audio coding
EP2688066A1 (fr) * 2012-07-16 2014-01-22 Thomson Licensing Procédé et appareil de codage de signaux audio HOA multicanaux pour la réduction du bruit, et procédé et appareil de décodage de signaux audio HOA multicanaux pour la réduction du bruit
US9473870B2 (en) 2012-07-16 2016-10-18 Qualcomm Incorporated Loudspeaker position compensation with 3D-audio hierarchical coding
US9761229B2 (en) * 2012-07-20 2017-09-12 Qualcomm Incorporated Systems, methods, apparatus, and computer-readable media for audio object clustering
US9479886B2 (en) 2012-07-20 2016-10-25 Qualcomm Incorporated Scalable downmix design with feedback for object-based surround codec
KR101676634B1 (ko) * 2012-08-31 2016-11-16 돌비 레버러토리즈 라이쎈싱 코오포레이션 오브젝트―기반 오디오를 위한 반사된 사운드 렌더링
US10203839B2 (en) 2012-12-27 2019-02-12 Avaya Inc. Three-dimensional generalized space
US9838824B2 (en) 2012-12-27 2017-12-05 Avaya Inc. Social media processing with three-dimensional audio
US9892743B2 (en) 2012-12-27 2018-02-13 Avaya Inc. Security surveillance via three-dimensional audio space presentation
US9301069B2 (en) * 2012-12-27 2016-03-29 Avaya Inc. Immersive 3D sound space for searching audio
US9736609B2 (en) 2013-02-07 2017-08-15 Qualcomm Incorporated Determining renderers for spherical harmonic coefficients
US9685163B2 (en) * 2013-03-01 2017-06-20 Qualcomm Incorporated Transforming spherical harmonic coefficients
WO2014162171A1 (fr) 2013-04-04 2014-10-09 Nokia Corporation Appareil de traitement audiovisuel
US9706324B2 (en) 2013-05-17 2017-07-11 Nokia Technologies Oy Spatial object oriented audio apparatus
US9716959B2 (en) 2013-05-29 2017-07-25 Qualcomm Incorporated Compensating for error in decomposed representations of sound fields
US9369818B2 (en) 2013-05-29 2016-06-14 Qualcomm Incorporated Filtering with binaural room impulse responses with content analysis and weighting
EP2824661A1 (fr) * 2013-07-11 2015-01-14 Thomson Licensing Procédé et appareil de génération à partir d'une représentation dans le domaine des coefficients de signaux HOA et représentation dans un domaine mixte spatial/coefficient de ces signaux HOA
DE102013013378A1 (de) * 2013-08-10 2015-02-12 Advanced Acoustic Sf Gmbh Aufteilung virtueller Schallquellen
WO2015054033A2 (fr) 2013-10-07 2015-04-16 Dolby Laboratories Licensing Corporation Système et procédé de traitement audio spatial
EP2866475A1 (fr) 2013-10-23 2015-04-29 Thomson Licensing Procédé et appareil pour décoder une représentation du champ acoustique audio pour lecture audio utilisant des configurations 2D
US9922656B2 (en) 2014-01-30 2018-03-20 Qualcomm Incorporated Transitioning of ambient higher-order ambisonic coefficients
EP2930958A1 (fr) * 2014-04-07 2015-10-14 Harman Becker Automotive Systems GmbH Génération d'un champ d'ondes sonores
US10770087B2 (en) 2014-05-16 2020-09-08 Qualcomm Incorporated Selecting codebooks for coding vectors decomposed from higher-order ambisonic audio signals
JP6388551B2 (ja) * 2015-02-27 2018-09-12 アルパイン株式会社 複数領域音場再現システムおよび方法
DE102015008000A1 (de) * 2015-06-24 2016-12-29 Saalakustik.De Gmbh Verfahren zur Schallwiedergabe in Reflexionsumgebungen, insbesondere in Hörräumen
EP3402223B1 (fr) 2016-01-08 2020-10-07 Sony Corporation Dispositif et procédé de traitement audio, et programme
WO2017119318A1 (fr) * 2016-01-08 2017-07-13 ソニー株式会社 Dispositif et procédé de traitement audio, et programme
AU2017305249B2 (en) 2016-08-01 2021-07-22 Magic Leap, Inc. Mixed reality system with spatialized audio
WO2018064528A1 (fr) * 2016-09-29 2018-04-05 The Trustees Of Princeton University Navigation ambisonique dans des champs sonores à partir d'un réseau de microphones
WO2018081829A1 (fr) * 2016-10-31 2018-05-03 Google Llc Codage audio par projection
FR3060830A1 (fr) * 2016-12-21 2018-06-22 Orange Traitement en sous-bandes d'un contenu ambisonique reel pour un decodage perfectionne
US10182303B1 (en) * 2017-07-12 2019-01-15 Google Llc Ambisonics sound field navigation using directional decomposition and path distance estimation
US10764684B1 (en) 2017-09-29 2020-09-01 Katherine A. Franco Binaural audio using an arbitrarily shaped microphone array
US10721559B2 (en) 2018-02-09 2020-07-21 Dolby Laboratories Licensing Corporation Methods, apparatus and systems for audio sound field capture
EP3759508B1 (fr) * 2018-03-02 2023-09-13 Sonitor Technologies AS Système et procédé d'émetteur et de récepteur de positionnement acoustique
US10771913B2 (en) 2018-05-11 2020-09-08 Dts, Inc. Determining sound locations in multi-channel audio
CN110740416B (zh) * 2019-09-27 2021-04-06 广州励丰文化科技股份有限公司 一种音频信号处理方法及装置
CN110740404B (zh) * 2019-09-27 2020-12-25 广州励丰文化科技股份有限公司 一种音频相关性的处理方法及音频处理装置
WO2021138517A1 (fr) 2019-12-30 2021-07-08 Comhear Inc. Procédé pour fournir un champ sonore spatialisé
CN111537058B (zh) * 2020-04-16 2022-04-29 哈尔滨工程大学 一种基于Helmholtz方程最小二乘法的声场分离方法
US11743670B2 (en) 2020-12-18 2023-08-29 Qualcomm Incorporated Correlation-based rendering with multiple distributed streams accounting for an occlusion for six degree of freedom applications
CN113791385A (zh) * 2021-09-15 2021-12-14 张维翔 一种三维定位方法及系统

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS53114201U (fr) * 1977-02-18 1978-09-11
US4731848A (en) * 1984-10-22 1988-03-15 Northwestern University Spatial reverberator
JP2569872B2 (ja) * 1990-03-02 1997-01-08 ヤマハ株式会社 音場制御装置
JP3578783B2 (ja) * 1993-09-24 2004-10-20 ヤマハ株式会社 電子楽器の音像定位装置
US5745584A (en) * 1993-12-14 1998-04-28 Taylor Group Of Companies, Inc. Sound bubble structures for sound reproducing arrays
GB9726338D0 (en) * 1997-12-13 1998-02-11 Central Research Lab Ltd A method of processing an audio signal
US7231054B1 (en) * 1999-09-24 2007-06-12 Creative Technology Ltd Method and apparatus for three-dimensional audio display
US7340062B2 (en) * 2000-03-14 2008-03-04 Revit Lawrence J Sound reproduction method and apparatus for assessing real-world performance of hearing and hearing aids
EP1275272B1 (fr) * 2000-04-19 2012-11-21 SNK Tech Investment L.L.C. Prise de son ambiant multi-canal et techniques de reproduction qui preservent les harmoniques spatiales en trois dimensions

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO2004049299A1 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9460726B2 (en) 2010-03-26 2016-10-04 Dolby Laboratories Licensing Corporation Method and device for decoding an audio soundfield representation for audio playback
US10582329B2 (en) 2016-01-08 2020-03-03 Sony Corporation Audio processing device and method

Also Published As

Publication number Publication date
BR0316718A (pt) 2005-10-18
CN1735922B (zh) 2010-05-12
FR2847376A1 (fr) 2004-05-21
KR100964353B1 (ko) 2010-06-17
FR2847376B1 (fr) 2005-02-04
JP2006506918A (ja) 2006-02-23
KR20050083928A (ko) 2005-08-26
EP1563485B1 (fr) 2006-03-29
ZA200503969B (en) 2006-09-27
WO2004049299A1 (fr) 2004-06-10
DE60304358D1 (de) 2006-05-18
BRPI0316718B1 (pt) 2021-11-23
CN1735922A (zh) 2006-02-15
DE60304358T2 (de) 2006-12-07
US20060045275A1 (en) 2006-03-02
ATE322065T1 (de) 2006-04-15
AU2003290190A1 (en) 2004-06-18
US7706543B2 (en) 2010-04-27
JP4343845B2 (ja) 2009-10-14
ES2261994T3 (es) 2006-11-16

Similar Documents

Publication Publication Date Title
EP1563485B1 (fr) Procede de traitement de donnees sonores et dispositif d&#39;acquisition sonore mettant en oeuvre ce procede
EP1600042B1 (fr) Procede de traitement de donnees sonores compressees, pour spatialisation
EP1992198B1 (fr) Optimisation d&#39;une spatialisation sonore binaurale a partir d&#39;un encodage multicanal
WO2006075077A2 (fr) Procede et dispositif d’individualisation de hrtfs par modelisation
EP1586220B1 (fr) Procede et dispositif de pilotage d&#39;un ensemble de restitution a partir d&#39;un signal multicanal
EP3475943B1 (fr) Procede de conversion et d&#39;encodage stereophonique d&#39;un signal audio tridimensionnel
JP5611970B2 (ja) オーディオ信号を変換するためのコンバータ及び方法
WO2003073791A2 (fr) Procédé et dispositif de pilotage d&#39;un ensemble de restitution d&#39;un champ acoustique
WO2005069272A1 (fr) Procede de synthese et de spatialisation sonores
EP3400599B1 (fr) Encodeur ambisonique ameliore d&#39;une source sonore a pluralite de reflexions
EP3025514B1 (fr) Spatialisation sonore avec effet de salle
FR3065137A1 (fr) Procede de spatialisation sonore
Ifergan et al. On the selection of the number of beamformers in beamforming-based binaural reproduction
WO2005096268A2 (fr) Procede de traitement de donnees sonores, en particulier en contexte ambiophonique
EP3384688B1 (fr) Décompositions successives de filtres audio
FR3040253B1 (fr) Procede de mesure de filtres phrtf d&#39;un auditeur, cabine pour la mise en oeuvre du procede, et procedes permettant d&#39;aboutir a la restitution d&#39;une bande sonore multicanal personnalisee
FR3069693B1 (fr) Procede et systeme de traitement d&#39;un signal audio incluant un encodage au format ambisonique
EP3484185A1 (fr) Modelisation d&#39;ensemble de fonctions de transferts acoustiques propre a un individu, carte son tridimensionnel et systeme de reproduction sonore tridimensionnelle
EP3449643B1 (fr) Procédé et système de diffusion d&#39;un signal audio à 360°
FR3114209A1 (fr) Systeme de reproduction de sons avec virtualisation du champ reverbere

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20050530

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PT RO SE SI SK TR

AX Request for extension of the european patent

Extension state: AL LT LV MK

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

DAX Request for extension of the european patent (deleted)
AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PT RO SE SI SK TR

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT;WARNING: LAPSES OF ITALIAN PATENTS WITH EFFECTIVE DATE BEFORE 2007 MAY HAVE OCCURRED AT ANY TIME BEFORE 2007. THE CORRECT EFFECTIVE DATE MAY BE DIFFERENT FROM THE ONE RECORDED.

Effective date: 20060329

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

Ref country code: IE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

Ref country code: NL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: FRENCH

REF Corresponds to:

Ref document number: 60304358

Country of ref document: DE

Date of ref document: 20060518

Kind code of ref document: P

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060629

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060629

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060629

GBT Gb: translation of ep patent filed (gb section 77(6)(a)/1977)

Effective date: 20060712

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060829

NLV1 Nl: lapsed or annulled due to failure to fulfill the requirements of art. 29p and 29m of the patents act
REG Reference to a national code

Ref country code: IE

Ref legal event code: FD4D

REG Reference to a national code

Ref country code: ES

Ref legal event code: FG2A

Ref document number: 2261994

Country of ref document: ES

Kind code of ref document: T3

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20061130

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20061130

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20070102

BERE Be: lapsed

Owner name: FRANCE TELECOM

Effective date: 20061130

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060630

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20071130

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20061113

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20071130

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060930

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20060329

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20101217

Year of fee payment: 8

REG Reference to a national code

Ref country code: FR

Ref legal event code: ST

Effective date: 20120731

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: FR

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20111130

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: IT

Payment date: 20221020

Year of fee payment: 20

Ref country code: GB

Payment date: 20221021

Year of fee payment: 20

Ref country code: ES

Payment date: 20221201

Year of fee payment: 20

Ref country code: DE

Payment date: 20221020

Year of fee payment: 20

REG Reference to a national code

Ref country code: DE

Ref legal event code: R071

Ref document number: 60304358

Country of ref document: DE

REG Reference to a national code

Ref country code: ES

Ref legal event code: FD2A

Effective date: 20231124

REG Reference to a national code

Ref country code: GB

Ref legal event code: PE20

Expiry date: 20231112

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GB

Free format text: LAPSE BECAUSE OF EXPIRATION OF PROTECTION

Effective date: 20231112

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: ES

Free format text: LAPSE BECAUSE OF EXPIRATION OF PROTECTION

Effective date: 20231114

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GB

Free format text: LAPSE BECAUSE OF EXPIRATION OF PROTECTION

Effective date: 20231112

Ref country code: ES

Free format text: LAPSE BECAUSE OF EXPIRATION OF PROTECTION

Effective date: 20231114