Method and apparatus for compressing and decompressing a Higher Order Ambisonics representation

Abstract

Higher Order Ambisonics represents three-dimensional sound independent of a specific loudspeaker set-up. However, transmission of an HOA representation results in a very high bit rate. Therefore compression with a fixed number of channels is used, in which directional and ambient signal components are processed differently. The ambient HOA component is represented by a minimum number of HOA coefficient sequences. The remaining channels contain either directional signals or additional coefficient sequences of the ambient HOA component, depending on what will result in optimum perceptual quality. This processing can change on a frame-by-frame basis.

Claims

1. Method for compressing using a fixed number of perceptual encodings a Higher Order Ambisonics representation of a sound field, denoted HOA, with input time frames of HOA coefficient sequences, said method comprising the following which is carried out on a frame-by-frame basis: for a current frame, estimating a set of dominant directions and a corresponding data set of indices of detected directional signals; separating from the HOA coefficient sequences of said current frame a non-fixed number of directional signals with respective directions contained in said set of dominant direction estimates and with a respective delayed data set of indices of said directional signals, wherein said non-fixed number is smaller than said fixed number, and an ambient HOA component that is represented by a reduced number of HOA coefficient sequences and a corresponding data set of indices of said reduced number of ambient HOA coefficient sequences, which reduced number corresponds to the difference between said fixed number and said non-fixed number; assigning said directional signals and the HOA coefficient sequences of said ambient HOA component to channels the number of which corresponds to said fixed number, wherein for said assigning said delayed data set of indices of said directional signals and said data set of indices of said reduced number of ambient HOA coefficient sequences are used; perceptually encoding said channels of the related frame so as to provide an encoded compressed frame.

2. Method according to claim 1, wherein said non-fixed number of directional signals is determined according to a perceptually related criterion such that: a correspondingly decompressed HOA representation provides a lowest perceptible error which can be achieved with the fixed given number of channels for the compression, wherein said criterion considers the following errors: the modelling errors arising from using different numbers of said directional signals and different numbers of HOA coefficient sequences for the ambient HOA component; the quantization noise introduced by the perceptual coding of said directional signals; the quantization noise introduced by coding the individual HOA coefficient sequences of said ambient HOA component; the total error, resulting from the above three errors, is considered for a number of test directions and a number of critical bands with respect to its perceptibility; said non-fixed number of directional signals is chosen so as to minimize the average perceptible error or the maximum perceptible error so as to achieve said lowest perceptible error.

3. Method according to claim 1, wherein the choice of the reduced number of HOA coefficient sequences to represent the ambient HOA component is carried out according to a criterion that differentiates between the following three cases: in case the number of HOA coefficient sequences for said current frame is the same as for the previous frame, the same HOA coefficient sequences are chosen as in said previous frame; in case the number of HOA coefficient sequences for said current frame is smaller than that for said previous frame, those HOA coefficient sequences from said previous frame are de-activated which were in said previous frame assigned to a channel that is in said current frame occupied by a directional signal; in case the number of HOA coefficient sequences for said current frame is greater than for said previous frame, those HOA coefficient sequences which were selected in said previous frame are also selected in said current frame, and these additional HOA coefficient sequences can be selected according to their perceptual significance or according the highest average power.

4. Method according to claim 1, wherein said assigning is carried out as follows: active directional signals are assigned to the given channels such that they keep their channel indices, in order to obtain continuous signals for said perceptual coding; the HOA coefficient sequences of said ambient HOA component are assigned such that a minimum number of such coefficient sequences is always contained in a corresponding number of last channels; for assigning additional HOA coefficient sequences of said ambient HOA component it is determined whether they were also selected in said previous frame: if true, the assignment of these HOA coefficient sequences to the channels to be perceptually encoded is the same as for said previous frame; if not true and if HOA coefficient sequences are newly selected, the HOA coefficient sequences are first arranged with respect to their indices in an ascending order and are in this order assigned to channels to be perceptually encoded which are not yet occupied by directional signals.

5. Method according to claim 1, wherein O.sub.RED is the number of HOA coefficient sequences representing said ambient HOA component, and wherein parameters describing said assignment are arranged in a bit array that has a length corresponding to an additional number of HOA coefficient sequences used in addition to the number O.sub.RED of HOA coefficient sequences for representing said ambient HOA component, and wherein each o-th bit in said bit array indicates whether the (O.sub.RED+o)-th additional HOA coefficient sequence is used for representing said ambient HOA component.

6. Method according to claim 1, wherein parameters describing said assignment are arranged in an assignment vector having a length corresponding to the number of inactive directional signals, the elements of which vector are indicating which of the additional HOA coefficient sequences of the ambient HOA component are assigned to the channels with inactive directional signals.

7. Method according to claim 1, wherein said separating of the HOA coefficient sequences of said current frame in addition provides parameters which can be used at decompression side for predicting portions of the original HOA representation from said directional signals.

8. Method according to claim 4, wherein said assigning provides an assignment vector, the elements of which vector are representing information about which of the additional HOA coefficient sequences for said ambient HOA component are assigned into the channels with inactive directional signals.

9. Apparatus for compressing using a fixed number of perceptual encodings a Higher Order Ambisonics representation of a sound field, denoted HOA, with input time frames of HOA coefficient sequences, said apparatus carrying out a frame-by-frame based processing and comprising: an estimator which estimates for a current frame a set of dominant directions and a corresponding data set of indices of detected directional signals; a separator which separates from the HOA coefficient sequences of said current frame a non-fixed number of directional signals with respective directions contained in said set of dominant direction estimates and with a respective delayed data set of indices of said directional signals, wherein said non-fixed number is smaller than said fixed number, and an ambient HOA component that is represented by a reduced number of HOA coefficient sequences and a corresponding data set of indices of said reduced number of ambient HOA coefficient sequences, which reduced number corresponds to the difference between said fixed number and said non-fixed number; an assignor which assigns said directional signals and the HOA coefficient sequences of said ambient HOA component to channels the number of which corresponds to said fixed number, thereby obtaining parameters of indices of the chosen ambient HOA coefficient sequences describing said assignment, which can be used for a corresponding re-distribution at a decompression side, wherein for said assigning said delayed data set of indices of said directional signals and said data set of indices of said reduced number of ambient HOA coefficient sequences are used; an encoder which perceptually encodes said channels of the related frame so as to provide an encoded compressed frame.

10. Apparatus according to claim 9, wherein said non-fixed number of directional signals is determined according to a perceptually related criterion such that: a correspondingly decompressed HOA representation provides a lowest perceptible error which can be achieved with the fixed given number of channels for the compression, wherein said criterion considers the following errors: the modelling errors arising from using different numbers of said directional signals and different numbers of HOA coefficient sequences for the ambient HOA component; the quantization noise introduced by the perceptual coding of said directional signals; the quantization noise introduced by coding the individual HOA coefficient sequences of said ambient HOA component; the total error, resulting from the above three errors, is considered for a number of test directions and a number of critical bands with respect to its perceptibility; said non-fixed number of directional signals is chosen so as to minimize the average perceptible error or the maximum perceptible error so as to achieve said lowest perceptible error.

11. Apparatus according to claim 9, wherein the choice of the reduced number of HOA coefficient sequences to represent the ambient HOA component is carried out according to a criterion that differentiates between the following three cases: in case the number of HOA coefficient sequences for said current frame is the same as for the previous frame, the same HOA coefficient sequences are chosen as in said previous frame; in case the number of HOA coefficient sequences for said current frame is smaller than that for said previous frame, those HOA coefficient sequences from said previous frame are de-activated which were in said previous frame assigned to a channel that is in said current frame occupied by a directional signal; in case the number of HOA coefficient sequences for said current frame is greater than for said previous frame, those HOA coefficient sequences which were selected in said previous frame are also selected in said current frame, and these additional HOA coefficient sequences can be selected according to their perceptual significance or according the highest average power.

12. Apparatus according to claim 9, wherein said assigning is carried out as follows: active directional signals are assigned to the given channels such that they keep their channel indices, in order to obtain continuous signals for said perceptual coding; the HOA coefficient sequences of said ambient HOA component are assigned such that a minimum number of such coefficient sequences is always contained in a corresponding number of last channels; for assigning additional HOA coefficient sequences of said ambient HOA component it is determined whether they were also selected in said previous frame: if true, the assignment of these HOA coefficient sequences to the channels to be perceptually encoded is the same as for said previous frame; if not true and if HOA coefficient sequences are newly selected, the HOA coefficient sequences are first arranged with respect to their indices in an ascending order and are in this order assigned to channels to be perceptually encoded which are not yet occupied by directional signals.

13. Apparatus according to claim 9, wherein O.sub.RED is the number of HOA coefficient sequences representing said ambient HOA component, and wherein parameters describing said assignment are arranged in a bit array that has a length corresponding to an additional number of HOA coefficient sequences used in addition to the number O.sub.RED of HOA coefficient sequences for representing said ambient HOA component, and wherein each o-th bit in said bit array indicates whether the (O.sub.RED+o)-th additional HOA coefficient sequence is used for representing said ambient HOA component.

14. Apparatus according to claim 9, wherein parameters describing said assignment are arranged in an assignment vector having a length corresponding to the number of inactive directional signals, the elements of which vector are indicating which of the additional HOA coefficient sequences of the ambient HOA component are assigned to the channels with inactive directional signals.

15. Apparatus according to claim 9, wherein said separating of the HOA coefficient sequences of said current frame in addition provides parameters which can be used at decompression side for predicting portions of the original HOA representation from said directional signals.

16. Apparatus according to claim 12, wherein said assigning provides an assignment vector, the elements of which vector are representing information about which of the additional HOA coefficient sequences for said ambient HOA component are assigned into the channels with inactive directional signals.

17. Digital audio signal that is compressed according to the method of claim 1.

18. Digital audio signal according to claim 17, which includes an assignment parameters bit array as defined in claim 13.

19. Digital audio signal according to claim 17, which includes an assignment vector.

20. Method for decompressing a Higher Order Ambisonics representation compressed according to the method of claim 1, said decompressing comprising: perceptually decoding a current encoded compressed frame so as to provide a perceptually decoded frame of channels; re-distributing said perceptually decoded frame of channels, using said data set of indices of directional signals and said data set of indices of the chosen ambient HOA coefficient sequences, so as to recreate the corresponding frame of directional signals and the corresponding frame of the ambient HOA component; re-composing a current decompressed frame of the HOA representation from said frame of directional signals and from said frame of the ambient HOA component, using said data set of indices of detected directional signals and said set of dominant direction estimates, wherein directional signals with respect to uniformly distributed directions are predicted from said directional signals, and thereafter said current decompressed frame is re-composed from said frame of directional signals, said predicted signals and said ambient HOA component.

21. Method according to claim 20, wherein said prediction of directional signals with respect to uniformly distributed directions is performed from said directional signals using said received parameters for said predicting.

22. Method according to claim 20, wherein in said re-distribution, instead of the data set of indices of detected directional signals and the data set of indices of the chosen ambient HOA coefficient sequences, a received assignment vector is used, the elements of which vector are representing information about which of the additional HOA coefficient sequences for said ambient HOA component are assigned into the channels with inactive directional signals.

23. Apparatus for decompressing a Higher Order Ambisonics representation compressed according to the method of claim 1, said apparatus comprising: a decoder which perceptually decodes a current encoded compressed frame so as to provide a perceptually decoded frame of channels; a re-distributor which re-distributes said perceptually decoded frame of channels, using said data set of indices of detected directional signals and said data set of indices of the chosen ambient HOA coefficient sequences, so as to recreate the corresponding frame of directional signals and the corresponding frame of the ambient HOA component; a re-composer which re-composes a current decompressed frame of the HOA representation from said frame of directional signals and from said frame of the ambient HOA component, using said data set of indices of detected directional signals and said set of dominant direction estimates, wherein directional signals with respect to uniformly distributed directions are predicted from said directional signals, and thereafter said current decompressed frame is re-composed from said frame of directional signals, said predicted signals and said ambient HOA component.

24. Apparatus according to claim 23, wherein said prediction of directional signals with respect to uniformly distributed directions is performed from said directional signals using said received parameters for said predicting.

25. Apparatus according to claim 23, wherein in said re-distribution, instead of the data set of indices of detected directional signals and the data set of indices of the chosen ambient HOA coefficient sequences, a received assignment vector is used, the elements of which vector are representing information about which of the additional HOA coefficient sequences for said ambient HOA component are assigned into the channels with inactive directional signals.

Description

BRIEF DESCRIPTION OF DRAWINGS

(1) Exemplary embodiments of the invention are described with reference to the accompanying drawings, which show in:

(2) FIG. 1 block diagram for the HOA compression;

(3) FIG. 2 estimation of dominant sound source directions;

(4) FIG. 3 block diagram for the HOA decompression;

(5) FIG. 4 spherical coordinate system;

(6) FIG. 5 normalised dispersion function υ.sub.N(Θ) for different Ambisonics orders N and for angles θ∈[0,π].

DESCRIPTION OF EMBODIMENTS

(7) A. Improved HOA Compression

(8) The compression processing according to the invention, which is based on EP 12306569.0, is illustrated in FIG. 1 where the signal processing blocks that have been modified or newly introduced compared to EP 12306569.0 are presented with a bold box, and where ‘’ (direction estimates as such) and ‘C’ in this application correspond to ‘A’ (matrix of direction estimates) and ‘D’ in EP 12306569.0, respectively. For the HOA compression a frame-wise processing with non-overlapping input frames C(k) of HOA coefficient sequences of length L is used, where k denotes the frame index. The frames are defined with respect to the HOA coefficient sequences specified in equation (45) as
C(k):=[c((kL+1)T.sub.s)c((kL+2)T.sub.s)c((k+1)LT.sub.s)],  (1)

(9) where T.sub.s indicates the sampling period.

(10) The first step or stage 11/12 in FIG. 1 is optional and consists of concatenating the non-overlapping k-th and the (k−1)-th frames of HOA coefficient sequences into a long frame {tilde over (C)}(k) as
{tilde over (C)}(k):=[C(k−1)C(k)],  (2)

(11) which long frame is 50% overlapped with an adjacent long frame and which long frame is successively used for the estimation of dominant sound source directions. Similar to the notation for {tilde over (C)}(k), the tilde symbol is used in the following description for indicating that the respective quantity refers to long overlapping frames. If step/stage 11/12 is not present, the tilde symbol has no specific meaning.

(12) In principle, the estimation step or stage 13 of dominant sound sources is carried out as proposed in EP 13305156.5, but with an important modification. The modification is related to the determination of the amount of directions to be detected, i.e. how many directional signals are supposed to be extracted from the HOA representation. This is accomplished with the motivation to extract directional signals only if it is perceptually more relevant than using instead additional HOA coefficient sequences for better approximation of the ambient HOA component. A detailed description of this technique is given in section A.2.

(13) The estimation provides a data set .sub.DIR,ACT(k) ⊂ {1, . . . , D} of indices of directional signals that have been detected as well as the set .sub.Ω,ACT(k) of corresponding direction estimates. D denotes the maximum number of directional signals that has to be set before starting the HOA compression.

(14) In step or stage 14, the current (long) frame {tilde over (C)}(k) of HOA coefficient sequences is decomposed (as proposed in EP 13305156.5) into a number of directional signals X.sub.DIR(k−2) belonging to the directions contained in the set .sub.Ω,ACT(k) and a residual ambient HOA component C.sub.AMB(k−2). The delay of two frames is introduced as a result of overlap-add processing in order to obtain smooth signals. It is assumed that X.sub.DIR(k−2) is containing a total of D channels, of which however only those corresponding to the active directional signals are non-zero. The indices specifying these channels are assumed to be output in the data set .sub.DIR,ACT(k−2). Additionally, the decomposition in step/stage 14 provides some parameters ζ(k−2) which are used at decompression side for predicting portions of the original HOA representation from the directional signals (see EP 13305156.5 for more details). In step or stage 15, the number of coefficients of the ambient HOA component C.sub.AMB(k−2) is intelligently reduced to contain only O.sub.RED+D−N.sub.DIR,ACT(k−2) non-zero HOA coefficient sequences, where N.sub.DIR,ACT(k−2)=|.sub.DIR,AcT(k−2)| indicates the cardinality of the data set .sub.DIR,ACT(k−2), i.e. the number of active directional signals in frame k−2. Since the ambient HOA component is assumed to be always represented by a minimum number O.sub.RED of HOA coefficient sequences, this problem can be actually reduced to the selection of the remaining D−N.sub.DIR,ACT(k−2) HOA coefficient sequences out of the possible O−O.sub.RED ones. In order to obtain a smooth reduced ambient HOA representation, this choice is accomplished such that, compared to the choice taken at the previous frame k−3, as few changes as possible will occur.

(15) In particular, the three following cases are to be differentiated: a) N.sub.DIR,ACT(k−2)=N.sub.DIR,ACT(k−3): In this case the same HOA coefficient sequences are assumed to be selected as in frame k−3. b) N.sub.DIR,ACT(k−2)<N.sub.DIR,ACT(k−3): In this case, more HOA coefficient sequences than in the last frame k−3 can be used for representing the ambient HOA component in the current frame. Those HOA coefficient sequences that were selected in k−3 are assumed to be also selected in the current frame. The additional HOA coefficient sequences can be selected according to different criteria. For instance, selecting those HOA coefficient sequences in C.sub.AMB(k−2) with the highest average power, or selecting the HOA coefficients sequences with respect to their perceptual significance. c) N.sub.DIR,ACT(k−2)>N.sub.DIR,ACT(k−3): In this case, less HOA coefficient sequences than in the last frame k−3 can be used for representing the ambient HOA component in the current frame. The question to be answered here is which of the previously selected HOA coefficient sequences have to be deactivated. A reasonable solution is to deactivate those sequences which were assigned to the channels i∈.sub.DIR,ACT(k−2) at the signal assigning step or stage 16 at frame k−3.

(16) For avoiding discontinuities at frame borders when additional HOA coefficient sequences are activated or deactivated, it is advantageous to smoothly fade in or out the respective signals.

(17) The final ambient HOA representation with the reduced number of O.sub.RED+N.sub.DIR,ACT(k−2) non-zero coefficient sequences is denoted by C.sub.AMB,RED(k−2). The indices of the chosen ambient HOA coefficient sequences are output in the data set .sub.AMB,ACT(k−2).

(18) In step/stage 16, the active directional signals contained in X.sub.DIR(k−2) and the HOA coefficient sequences contained in C.sub.AMB,RED(k−2) are assigned to the frame Y(k−2) of I channels for individual perceptual encoding. To describe the signal assignment in more detail, the frames X.sub.DIR(k−2), Y(k−2) and C.sub.AMB,RED(k−2) are assumed to consist of the individual signals X.sub.DIR,d(k−2), d∈ {1, . . . , D}, y.sub.i(k−2), i∈ {1, . . . , I} and C.sub.AMB,RED,o(k−2), o∈ {1, . . . , O} as follows:

(19) $\begin{array}{cc}{X}_{\mathrm{DIR}}\left(k-2\right)=\left[\begin{array}{c}{x}_{\mathrm{DIR},1}\left(k-2\right)\\ x_{\mathrm{DIR},2}\left(k-2\right)\\ \mathrm{.Math.}\\ x_{\mathrm{DIR},D}\left(k-2\right)\end{array}\right],& \left(3\right)\\ C_{\mathrm{AMB},\mathrm{RED}}\left(k-2\right)=\left[\begin{array}{c}{c}_{\mathrm{AMB},\mathrm{RED}1}\left(k-2\right)\\ c_{\mathrm{AMB},\mathrm{RED},2}\left(k-2\right)\\ \mathrm{.Math.}\\ c_{\mathrm{AMB},\mathrm{RED},O}\left(k-2\right)\end{array}\right],& \\ Y\left(k-2\right)=\left[\begin{array}{c}{y}_1\left(k-2\right)\\ y_2\left(k-2\right)\\ \mathrm{.Math.}\\ y_I\left(k-2\right)\end{array}\right].& \end{array}$

(20) The active directional signals are assigned such that they keep their channel indices in order to obtain continuous signals for the successive perceptual coding. This can be expressed by
y.sub.d(k−2)=x.sub.DIR,d(k−2) for all d∈.sub.DIR,ACT(k−2).  (4)

(21) The HOA coefficient sequences of the ambient component are assigned such the minimum number of O.sub.RED coefficient sequences is always contained in the last O.sub.RED signals of Y(k−2), i.e.
y.sub.D+o(k−2)=c.sub.AMB,RED,o(k−2) for 1≦o≦O.sub.RED.  (5)

(22) For the additional D−N.sub.DIR,ACT(k−2) HOA coefficient sequences of the ambient component it is to be differentiated whether or not they were also selected in the previous frame: a) If they were also selected to be transmitted in the previous frame, i.e. if the respective indices are also contained in data set .sub.AMB,ACT(k−3), the assignment of these coefficient sequences to the signals in Y(k−2) is the same as for the previous frame. This operation assures smooth signals y.sub.i(k−2), which is favourable for the successive perceptual coding in step or stage 17. b) Otherwise, if some coefficient sequences are newly selected, i.e. if their indices are contained in data set .sub.AMB,ACT(k−2) but not in data set .sub.AMB,ACT(k−3), they are first arranged with respect to their indices in an ascending order and are in this order assigned to channels i .Math. .sub.DIR,ACT(k−2) of Y(k−2) which are not yet occupied by directional signals. This specific assignment offers the advantage that, during a HOA decompression process, the signal re-distribution and composition can be performed without the knowledge about which ambient HOA coefficient sequence is contained in which channel of Y(k−2). Instead, the assignment can be reconstructed during HOA decompression with the mere knowledge of the data sets .sub.AMB,ACT(k−2) and .sub.DIR,ACT(k).

(23) Advantageously, this assigning operation also provides the assignment vector γ(k)∈.sup.D−N.sup.DIR,ACT(k−2), whose elements γ.sub.o(k), o=1, . . . , D−N.sub.DIR,ACT(k−2), denote the indices of each one of the additional D−N.sub.DIR,ACT(k−2) HOA coefficient sequences of the ambient component. To say it differently, the elements of the assignment vector γ(k) provide information about which of the additional O−O.sub.RED HOA coefficient sequences of the ambient HOA component are assigned into the D−N.sub.DIR,ACT(k−2) channels with inactive directional signals. This vector can be transmitted additionally, but less frequently than by the frame rate, in order to allow for an initialisation of the re-distribution procedure performed for the HOA decompression (see section B). Perceptual coding step/stage 17 encodes the I channels of frame Y(k−2) and outputs an encoded frame {hacek over (Y)}(k−2).

(24) For frames for which vector γ(k) is not transmitted from step/stage 16, at decompression side the data parameter sets .sub.DIR,ACT(k) and .sub.AMB,ACT(k−2) instead of vector γ(k) are used for the performing the re-distribution.

(25) A.1 Estimation of the Dominant Sound Source Directions

(26) The estimation step/stage 13 for dominant sound source directions of FIG. 1 is depicted in FIG. 2 in more detail. It is essentially performed according to that of EP 13305156.5, but with a decisive difference, which is the way of determining the amount of dominant sound sources, corresponding to the number of directional signals to be extracted from the given HOA representation. This number is significant because it is used for controlling whether the given HOA representation is better represented either by using more directional signals or instead by using more HOA coefficient sequences to better model the ambient HOA component.

(27) The dominant sound source directions estimation starts in step or stage 21 with a preliminary search for the dominant sound source directions, using the long frame {tilde over (C)}(k) of input HOA coefficient sequences. Along with the preliminary direction estimates {tilde over (Ω)}.sub.DOM.sup.(d)(k), 1≦d≦D, the corresponding directional signals {tilde over (x)}.sub.DOM.sup.(d)(k) and the HOA sound field components {tilde over (C)}.sub.DOM,CORR.sup.(d)(k), which are supposed to be created by the individual sound sources, are computed as described in EP 13305156.5. In step or stage 22, these quantities are used together with the frame {tilde over (C)}(k) of input HOA coefficient sequences for determining the number {tilde over (D)}(k) of directional signals to be extracted. Consequently, the direction estimates {tilde over (Ω)}.sub.DOM.sup.(d)(k), {tilde over (D)}(k)<d≦D, the corresponding directional signals{tilde over (x)}.sub.DOM.sup.(d)(k), and HOA sound field components {tilde over (C)}.sub.DOM,CORR.sup.(d)(k) are discarded. Instead, only the direction estimates {tilde over (Ω)}.sub.DOM.sup.(d)(k), 1≦d≦{tilde over (D)}(k) are then assigned to previously found sound sources.

(28) In step or stage 23, the resulting direction trajectories are smoothed according to a sound source movement model and it is determined which ones of the sound sources are supposed to be active (see EP 13305156.5). The last operation provides the set .sub.DIR,ACT(k) of indices of active directional sound sources and the set .sub.Ω,ACT(k) of the corresponding direction estimates.

(29) A.2 Determination of Number of Extracted Directional Signals

(30) For determining the number of directional signals in step/stage 22, the situation is assumed that there is a given total amount of I channels which are to be exploited for capturing the perceptually most relevant sound field information. Therefore the number of directional signals to be extracted is determined, motivated by the question whether for the overall HOA compression/decompression quality the current HOA representation is represented better by using either more directional signals, or more HOA coefficient sequences for a better modelling of the ambient HOA component. To derive in step/stage 22 a criterion for the determination of the number of directional sound sources to be extracted, which criterion is related to the human perception, it is taken into consideration that HOA compression is achieved in particular by the following two operations: reduction of HOA coefficient sequences for representing the ambient HOA component (which means reduction of the number of related channels); perceptual encoding of the directional signals and of the HOA coefficient sequences for representing the ambient HOA component.

(31) Depending on the number M, 0≦M≦D, of extracted directional signals, the first operation results in the approximation

(32) $\begin{array}{cc}\tilde{C}\left(k\right)\approx {\tilde{C}}^{\left(M\right)}\left(k\right)& \left(6\right)\\ :={\tilde{C}}_{\mathrm{DIR}}^{\left(M\right)}\left(k\right)+{\tilde{C}}_{\mathrm{AMB},\mathrm{RED}}^{\left(M\right)}\left(k\right),& \left(7\right)\\ \mathrm{where}{\tilde{C}}_{\mathrm{DIR}}^{M)}\left(k\right)\text{:}=\underset{d=1}{\overset{M}{.Math.}}{\tilde{C}}_{\mathrm{DOM},\mathrm{CORR}}^{\left(d\right)}\left(k\right)& \left(8\right)\end{array}$
denotes the HOA representation of the directional component consisting of the HOA sound field components {tilde over (C)}.sub.DOM,CORR.sup.(d)(k), 1≦d≦M, supposed to be created by the M individually considered sound sources, and {tilde over (C)}.sub.AMB,RED.sup.(M)(k) denotes the HOA representation of the ambient component with only I−M non-zero HOA coefficient sequences.

(33) The approximation from the second operation can be expressed by

(34) $\begin{array}{cc}\tilde{C}\left(k\right)\approx {\hat{\tilde{C}}}^{\left(M\right)}\left(k\right)& \left(9\right)\\ :={\hat{\tilde{C}}}_{\mathrm{DIR}}^{\left(M\right)}\left(k\right)+{\hat{\tilde{C}}}_{\mathrm{AMB},\mathrm{RED}}^{\left(M\right)}\left(k\right)& \left(10\right)\end{array}$
where {tilde over (Ĉ)}.sub.DIR.sup.(M)(k) and {tilde over (Ĉ)}.sub.AMB,RED.sup.(M)(k) denote the composed directional and ambient HOA components after perceptual decoding, respectively.

(35) Formulation of Criterion

(36) The number {tilde over (D)}(k) of directional signals to be extracted is chosen such that the total approximation error
{tilde over (Ê)}.sup.(M)(k):={tilde over (C)}(k)−{tilde over (Ĉ)}.sup.(M)(k)  (11)
with M={tilde over (D)}(k) is as less significant as possible with respect to the human perception. To assure this, the directional power distribution of the total error for individual Bark scale critical bands is considered at a predefined number Q of test directions Ω.sub.q, q=1, . . . , Q, which are nearly uniformly distributed on the unit sphere. To be more specific, the directional power distribution for the b-th critical band, b=1, . . . , B, is represented by the vector
.sup.(M)(k,b):=[.sub.1.sup.(M)(k,b).sub.2.sup.(M)(k,b) . . . .sub.Q.sup.(M)(k,b)].sup.T,  (12)

(37) whose components .sub.q.sup.(M)(k,b) denote the power of the total error {tilde over (Ê)}.sup.(M)(k) related to the direction Ω.sub.q, the b-th Bark scale critical band and the k-th frame. The directional power distribution .sup.(M)(k,b) of the total error {tilde over (Ê)}.sup.(M)(k) is compared with the directional perceptual masking power distribution
.sub.MASK(k,b):=[.sub.MASK,1(k,b).sub.MASK,2(k,b) . . . .sub.MASK,Q(k,b)].sup.T  (13)

(38) due to the original HOA representation Ĉ(k). Next, for each test direction Ω.sub.q and critical band b the level of perception .sub.q.sup.(M)(k,b) of the total error is computed. It is here essentially defined as the ratio of the directional power of the total error {tilde over (Ê)}.sup.(M)(k) and the directional masking power according to

(39) $\begin{array}{cc}{\widetilde{ℒ}}_q^{\left(M\right)}\left(k,b\right)\text{:}=\max \left(0,\frac{{\widehat{\tilde{\mathrm{��}}}}_q^{\left(M\right)}\left(k,b\right)}{{\widetilde{\mathrm{��}}}_{\mathrm{MASK},q}\left(k,b\right)}-1\right).& \left(14\right)\end{array}$

(40) The subtraction of ‘1’ and the successive maximum operation is performed to ensure that the perception level is zero, as long as the error power is below the masking threshold. Finally, the number {circumflex over (D)}(k) of directionals signals to be extracted can be chosen to minimise the average over all test directions of the maximum of the error perception level over all critical bands, i.e.,

(41) $\begin{array}{cc}\tilde{D}\left(k\right)=\underset{M}{\mathrm{argmin}}\frac{1}{Q}\underset{q=1}{\overset{Q}{.Math.}}\underset{b}{\max }{\widetilde{ℒ}}_q^{\left(M\right)}\left(k,b\right).& \left(15\right)\end{array}$

(42) It is noted that, alternatively, it is possible to replace the maximum by an averaging operation in equation (15).

(43) Computation of the Directional Perceptual Masking Power Distribution

(44) For the computation of the directional perceptual masking power distribution {tilde over (∵)}.sub.MASK(k,b) due to the original HOA representation {tilde over (C)}(k), the latter is transformed to the spatial domain in order to be represented by general plane waves {tilde over (υ)}.sub.q(k) impinging from the test directions Ω.sub.q, q=1, . . . , Q. When arranging the general plane wave signals {tilde over (υ)}.sub.q(k) in the matrix {tilde over (υ)}(k) as

(45) $\begin{array}{cc}\tilde{V}\left(k\right)=\left[\begin{array}{c}{\tilde{v}}_1\left(k\right)\\ {\tilde{v}}_2\left(k\right)\\ \mathrm{.Math.}\\ {\tilde{v}}_Q\left(k\right)\end{array}\right],& \left(16\right)\end{array}$

(46) the transformation to the spatial domain is expressed by the operation
{tilde over (V)}(k)=Ξ.sup.T{tilde over (C)}(k),  (17)

(47) where Ξ denotes the mode matrix with respect to the test direction Ω.sub.q, q=1, . . . , Q, defined by
Ξ:=[S.sub.1S.sub.2 . . . S.sub.Q]∈.sup.O×Q  (18)
with S.sub.q:=[S.sub.0.sup.0(Ω.sub.q)S.sub.−1.sup.−1(Ω.sub.q)S.sub.−1.sup.0(Ω.sub.q)S.sub.−1.sup.1(Ω.sub.q)S.sub.−2.sup.−2(Ω.sub.q) . . . S.sub.N.sup.N(Ω.sub.q)].sup.T∈.sup.o.  (19)

(48) The elements .sub.MASK(k,b) of the directional perceptual masking power distribution .sub.MASK(k,b), due to the original HOA representation {tilde over (C)}(k), are corresponding to the masking powers of the general plane wave functions {tilde over (υ)}.sub.q(k) for individual critical bands b.

(49) Computation of Directional Power Distribution

(50) In the following two alternatives for the computation of the directional power distribution .sup.(M)(k,b) are presented: a. One possibility is to actually compute the approximation {tilde over (Ĉ)}.sup.(M)(k) of the desired HOA representation {tilde over (C)}(k) by performing the two operations mentioned at the beginning of section A.2. Then the total approximation error {tilde over (Ê)}.sup.(M)(k) is computed according to equation (11). Next, the total approximation error {tilde over (Ê)}.sup.(M)(k) is transformed to the spatial domain in order to be represented by general plane waves {tilde over (Ŵ)}.sub.q.sup.(M)(k) impinging from the test directions Ω.sub.q, q=1, . . . , Q. Arranging the general plane wave signals in the matrix {tilde over (Ŵ)}.sup.(M)(k) as

(51) $\begin{array}{cc}{\widehat{\stackrel{\sim }{W}}}^{\left(M\right)}\left(k\right)=\left[\begin{array}{c}{\hat{\tilde{w}}}_1^{\left(M\right)}\left(k\right)\\ {\hat{\tilde{w}}}_2^{\left(M\right)}\left(k\right)\\ \mathrm{.Math.}\\ {\hat{\tilde{w}}}_{\mathrm{Q1}}^{\left(M\right)}\left(k\right)\end{array}\right],& \left(20\right)\end{array}$
the transformation to the spatial domain is expressed by the operation
{tilde over (Ŵ)}.sup.(M)(k)=Ξ.sup.T{tilde over (Ê)}.sup.(M)(k).  (21)

(52) The elements .sub.q.sup.(M)(k,b) of the directional power distribution .sup.(M)(k,b) of the total approximation error {tilde over (Ê)}.sup.(M)(k) are obtained by computing the powers of the general plane wave functions {tilde over (Ŵ)}.sub.q.sup.(M)(k), q=1, . . . , Q, within individual critical bands b. b. The alternative solution is to compute only the approximation {tilde over (C)}.sup.(M)(k) instead of {tilde over (Ĉ)}.sup.(M)(k). This method offers the advantage that the complicated perceptual coding of the individual signals needs not be carried out directly. Instead, it is sufficient to know the powers of the perceptual quantisation error within individual Bark scale critical bands. For this purpose, the total approximation error defined in equation (11) can be written as a sum of the three following approximation errors:
{tilde over (E)}.sup.(M)(k):={tilde over (C)}(k)−{tilde over (C)}.sup.(M)(k)  (22)
{tilde over (Ê)}.sub.DIR.sup.(M)(k):={tilde over (C)}.sub.DIR.sup.(M)(k)−{tilde over (Ĉ)}.sub.DIR.sup.(M)(k)  (23)
{tilde over (Ê)}.sub.AMB,RED.sup.(M)(k):={tilde over (C)}.sub.AMB,RED.sup.(M)(k)−{tilde over (Ĉ)}.sub.AMB,RED.sup.(M)(k),  (24) which can be assumed to be independent of each other. Due to this independence, the directional power distribution of the total error {tilde over (Ê)}.sup.(M)(k) can be expressed as the sum of the directional power distributions of the three individual errors {tilde over (E)}.sup.(M)(k), {tilde over (Ê)}.sub.DIR.sup.(M)(k) and {tilde over (Ê)}.sub.AMB,RED.sup.(M)(k).

(53) The following describes how to compute the directional power distributions of the three errors for individual Bark scale critical bands: a. To compute the directional power distribution of the error {tilde over (E)}.sup.(M)(k), it is first transformed to the spatial domain by
{tilde over (W)}.sup.(M)(k)=Ξ.sup.T{tilde over (E)}.sup.(M)(k),  (25) wherein the approximation error {tilde over (E)}.sup.(M)(k) is hence represented by general plane waves {tilde over (w)}.sub.q.sup.(M)(k) impinging from the test directions Ω.sub.q, q=1, . . . , Q, which are arranged in the matrix {tilde over (W)}.sup.(M)(k) according to

(54) $\begin{array}{cc}{\stackrel{\sim }{W}}^{\left(M\right)}\left(k\right)=\left[\begin{array}{c}{\tilde{w}}_1^{\left(M\right)}\left(k\right)\\ {\tilde{w}}_2^{\left(M\right)}\left(k\right)\\ \mathrm{.Math.}\\ {\tilde{w}}_{\mathrm{Q1}}^{\left(M\right)}\left(k\right)\end{array}\right].& \left(26\right)\end{array}$ Consequently, the elements .sub.q.sup.(M)(k,b) of the directional power distribution .sup.(M)(k,b) of the approximation error {tilde over (E)}.sup.(M)(k) are obtained by computing the powers of the general plane wave functions {tilde over (w)}.sub.q.sup.(M)(k), q=1, . . . , Q, within individual critical bands b. b. For computing the directional power distribution .sub.DIR.sup.(M)(k,b) of the error {tilde over (Ê)}.sub.DIR.sup.(M)(k), it is to be borne in mind that this error is introduced into the directional HOA component {tilde over (C)}.sub.DIR.sup.(M)(k) by perceptually coding the directional signals {tilde over (x)}.sub.DOM.sup.(d)(k), 1≦d≦M. Further, it is to be considered that the directional HOA component is given by equation (8). Then for simplicity it is assumed that the HOA component {tilde over (C)}.sub.DOM,CORR.sup.(d)(k) is equivalently represented in the spatial domain by O general plane wave functions {tilde over (υ)}.sub.GRID,o.sup.(d)(k), which are created from the directional signal {tilde over (x)}.sub.DOM.sup.(d)(k) by a mere scaling, i.e.
{tilde over (υ)}.sub.GRID,o.sup.(d)(k)=α.sub.o.sup.(d)(k){tilde over (x)}.sub.DOM.sup.(d)(k),  (27) where α.sub.o.sup.(d)(k), o=1, . . . , O, denote the scaling parameters. The respective plane wave directions {tilde over (Ω)}.sub.ROT,o.sup.(d)(k), o=1, . . . , O, are assumed to be uniformly distributed on the unit sphere and rotated such that {tilde over (Ω)}.sub.ROT,1.sup.(d)(k) corresponds to the direction estimate {tilde over (Ω)}.sub.DOM.sup.(d)(k). Hence, the scaling parameter α.sub.1.sup.(d)(k) is equal to ‘1’. When defining Ξ.sub.GRID.sup.(d)(k) to be the mode matrix with respect to the rotated directions {tilde over (Ω)}.sub.ROT,o.sup.(d)(k), o=1, . . . , O, and arranging all scaling parameters α.sub.o.sup.(d)(k) in a vector according to
α.sup.(d)(k):=[1α.sub.2.sup.(d)(k)α.sub.3.sup.(d)(k) . . . α.sub.0.sup.(d)(k)].sup.T∈.sup.o,  (28) the HOA component {tilde over (C)}.sub.DOM,CORR.sup.(d)(k) can be written as
{tilde over (C)}.sub.DOM,CORR.sup.(d)(k)=Ξ.sub.GRID.sup.(d)(k)α.sup.(d)(k){tilde over (x)}.sub.DOM.sup.(d)(k)  (29) Consequently, the error {tilde over (Ê)}.sub.DIR.sup.(M)(k) (see equation (23)) between the true directional HOA component
{tilde over (C)}.sub.DIR.sup.(M)(k)=Σ.sub.d=1.sup.M{tilde over (C)}.sub.DOM,CORR.sup.(d)(k)  (30) and that composed from the perceptually decoded directional signals {tilde over ({circumflex over (x)})}.sub.DOM.sup.(d)(k), d=1, . . . M, by

(55) $\begin{array}{cc}\begin{array}{c}{\hat{\tilde{C}}}_{\mathrm{DIR}}^{\left(M\right)}\left(k\right)=\underset{d=1}{\overset{M}{.Math.}}{\hat{\tilde{C}}}_{\mathrm{DOM},\mathrm{CORR}}^{\left(d\right)}\left(k\right)\\ :=\underset{d=1}{\overset{M}{.Math.}}{\Xi }_{\mathrm{GRID}}^{\left(d\right)}\left(k\right){\alpha }^{\left(d\right)}\left(k\right){\hat{\tilde{x}}}_{\mathrm{DOM}}^{\left(d\right)}\left(k\right)\left(32\right)\end{array}& \left(31\right)\end{array}$ can be expressed in terms of the perceptual coding errors
{tilde over ({circumflex over (e)})}.sub.DOM.sup.(d)(k):={tilde over (x)}.sub.DOM.sup.(d)(k)−{tilde over ({circumflex over (x)})}.sub.DOM.sup.(d)(k)  (33) in the individual directional signals by
{tilde over ({circumflex over (E)})}.sub.DIR.sup.(M)(k)=Σ.sub.d=1.sup.MΞ.sub.GRID.sup.(d)(k)α.sup.(d)(k){tilde over ({circumflex over (e)})}.sub.DOM.sup.(d)(k).  (34) The representation of the error {tilde over (Ê)}.sub.DIR.sup.(M)(k) in the spatial domain with respect to the test directions Ω.sub.q, q=1, . . . Q, is given by

(56) 0$\begin{array}{cc}{\widehat{\stackrel{\sim }{W}}}_{\mathrm{DIR},q}^{\left(M\right)}\left(d\right)=\underset{d=1}{\overset{M}{.Math.}}\underset{\underset{=:{\beta }^{\left(d\right)}\left(k\right)}{︸}}{{\Xi }^T{\Xi }_{\mathrm{GRID}}^{\left(d\right)}\left(k\right){\alpha }^{\left(d\right)}\left(k\right)}{\hat{\tilde{e}}}_{\mathrm{DOM}}^{\left(d\right)}\left(k\right).& \left(35\right)\end{array}$ Denoting the elements of the vector β.sup.(d)(k) by β.sub.q.sup.(d)(k), q=1, . . . , Q, and assuming the individual perceptual coding errors {tilde over (ê)}.sub.DOM.sup.(d)(k), d=1, . . . M, to be independent of each other, it follows from equation (35) that the elements .sub.DIR,q.sup.(M)(k,b) of the directional power distribution .sub.DIR.sup.(M)(k,b) of the perceptual coding error {tilde over (Ê)}.sub.DIR.sup.(M)(k) can be computed by
.sub.DIR,q.sup.(M)(k,b)=Σ.sub.d=1.sup.M(β.sub.q.sup.(d)(k)).sup.2{tilde over (σ)}.sub.DIR,d.sup.2(k,b).  (36) {tilde over (σ)}DIR,d.sup.2(k,b) is supposed to represent the power of the perceptual quantisation error within the b-th critical band in the directional signal {tilde over ({circumflex over (x)})}.sub.DOM.sup.(d)(k). This power can be assumed to correspond to the perceptual masking power of the directional signal {tilde over (x)}.sub.DOM.sup.(d)(k). c. For computing the directional power distribution AMB,RED.sup.(M)(k,b) of the error {tilde over (Ê)}.sub.AMB,RED.sup.(M)(k) resulting from the perceptual coding of the HOA coefficient sequences of the ambient HOA component, each HOA coefficient sequence is assumed to be coded independently. Hence, the errors introduced into the individual HOA coefficient sequences within each Bark scale critical band can be assumed to be uncorrelated. This means that the inter-coefficient correlation matrix of the error {tilde over (Ê)}.sub.AMB,RED.sup.(M)(k) with respect to each Bark scale critical band is diagonal, i.e.
{tilde over (Σ)}.sub.AMB,RED.sup.(M)(k,b)=diag({tilde over (σ)}.sub.AMB,RED,1.sup.2(M)(k,b),{tilde over (σ)}.sub.AMB,RED,2.sup.2(M)(k,b), . . . ,{tilde over (σ)}.sub.AMB,RED,O.sup.2(M)(k,b)).  (37) The elements {tilde over (σ)}.sub.AMB,RED,o.sup.2(M)(k,b), o=1, . . . , O, are supposed to represent sent the power of the perceptual quantisation error within the b-th critical band in the o-th coded HOA coefficient sequence in {tilde over (Ĉ)}.sub.AMB,RED.sup.(M)(k). They can be assumed to correspond to the perceptual masking power of the o-th HOA coefficient {tilde over (C)}.sub.AMB,RED.sup.(M)(k). The directional power distribution of the perceptual coding error {tilde over (Ê)}.sub.AMB,RED.sup.(M)(k) is thus computed by
.sub.AMB,RED.sup.(M)(k,b)=diag(Ξ.sup.T{tilde over (Σ)}.sub.AMB,RED.sup.(M)(k,b)Ξ).  (38)

(57) B. Improved HOA Decompression

(58) The corresponding HOA decompression processing is depicted in FIG. 3 and includes the following steps or stages. In step or stage 31 a perceptual decoding of the I signals contained in {hacek over (Y)}(k−2) is performed in order to obtain the I decoded signals in Ŷ(k−2).

(59) In signal re-distributing step or stage 32, the perceptually decoded signals in Ŷ(k−2) are re-distributed in order to recreate the frame {circumflex over (X)}.sub.DIR(k−2) of directional signals and the frame Ĉ.sub.AMB,RED(k−2) of the ambient HOA component. The information about how to re-distribute the signals is obtained by reproducing the assigning operation performed for the HOA compression, using the index data sets .sub.DIR,ACT(k) and .sub.AMB,ACT(k−2). Since this is a recursive procedure (see section A), the additionally transmitted assignment vector γ(k) can be used in order to allow for an initialisation of the re-distribution procedure, e.g. in case the transmission is breaking down.

(60) In composition step or stage 33, a current frame Ĉ(k−3) of the desired total HOA representation is re-composed (according to the processing described in connection with FIG. 2b and FIG. 4 of EP 12306569.0 using the frame {circumflex over (X)}.sub.DIR(k−2) of the directional signals, the set .sub.DIR,ACT(k) of the active directional signal indices together with the set .sub.Ω,ACT(k) of the corresponding directions, the parameters ζ(k−2) for predicting portions of the HOA representation from the directional signals, and the frame {tilde over (C)}.sub.AMB,RED(k−2) of HOA coefficient sequences of the reduced ambient HOA component. Ĉ.sub.AMB,RED(k−2) corresponds to component {circumflex over (D)}.sub.A(k−2) in EP 12306569.0, and .sub.Ω,ACT(k) and .sub.DIR,ACT(k) correspond to A.sub.{circumflex over (Ω)}(k) in EP 12306569.0, wherein active directional signal indices are marked in the matrix elements of A.sub.{circumflex over (Ω)}(k). I.e., directional signals with respect to uniformly distributed directions are predicted from the directional signals ({circumflex over (X)}.sub.DIR(k−2)) using the received parameters (ζ(k−2)) for such prediction, and thereafter the current decompressed frame (Ĉ(k−3)) is re-composed from the frame of directional signals ({circumflex over (X)}.sub.DIR(k−2)), the predicted portions and the reduced ambient HOA component (Ĉ.sub.AMB,RED(k−2).

(61) C. Basics of Higher Order Ambisonics

(62) Higher Order Ambisonics (HOA) is based on the description of a sound field within a compact area of interest, which is assumed to be free of sound sources. In that case the spatiotemporal behaviour of the sound pressure p(t,x) at time t and position x within the area of interest is physically fully determined by the homogeneous wave equation. In the following a spherical coordinate system as shown in FIG. 4 is assumed. In the used coordinate system the x axis points to the frontal position, the y axis points to the left, and the z axis points to the top. A position in space x=(r,θ,φ).sup.T is represented by a radius r>0 (i.e. the distance to the coordinate origin), an inclination angle θ ∈ [0,π] measured from the polar axis z and an azimuth angle φ ∈ [0,2π] measured counter-clockwise in the x−y plane from the x axis. Further, (.Math.).sup.T denotes the transposition.

(63) It can be shown (see E. G. Williams, “Fourier Acoustics”, volume 93 of Applied Mathematical Sciences, Academic Press, 1999) that the Fourier transform of the sound pressure with respect to time denoted by .sub.t(.Math.) i.e.
P(ω,x)=.sub.t(p(t,x))=∫.sub.−∞.sup.∞p(t,x)e.sup.−iωtdt,  (39)

(64) with ω denoting the angular frequency and i indicating the imaginary unit, can be expanded into a series of Spherical Harmonics according to
P(ω=kc.sub.s,r,θ,φ)=Σ.sub.n=0.sup.NΣ.sub.m=−n.sup.nA.sub.n.sup.m(k)j.sub.n(kr)S.sub.n.sup.m(θ,φ).  (40)

(65) In equation (40), c.sub.s denotes the speed of sound and k denotes the angular wave number, which is related to the angular frequency ω by

(66) $k=\frac{\omega }{c_s}.$
Further, j.sub.n(.Math.) denote the spherical Bessel functions of the first kind and S.sub.n.sup.m(θ,φ) denote the real valued Spherical Harmonics of order n and degree m, which are defined in below section C.1. The expansion coefficients A.sub.n.sup.m(k) are depending only on the angular wave number k. In the foregoing it has been implicitly assumed that sound pressure is spatially band-limited. Thus the series of Spherical Harmonics is truncated with respect to the order index n at an upper limit N, which is called the order of the HOA representation.

(67) If the sound field is represented by a superposition of an infinite number of harmonic plane waves of different angular frequencies ω arriving from all possible directions specified by the angle tuple (θ,φ), it can be shown (see B. Rafaely, “Plane-wave Decomposition of the Sound Field on a Sphere by Spherical Convolution”, Journal of the Acoustical Society of America, vol. 4(116), pages 2149-2157, 2004) that the respective plane wave complex amplitude function C(ω,θ,φ) can be expressed by the following Spherical Harmonics expansion
C(ω=kc.sub.s,θ,φ)=Σ.sub.n=0.sup.NΣ.sub.m=−n.sup.nC.sub.n.sup.m(k)S.sub.n.sup.m(θ,φ),  (41)

(68) where the expansion coefficients C.sub.n.sup.m(k) are related to the expansion coefficients A.sub.n.sup.m(k) by
A.sub.n.sup.m(k)=4πi.sup.nC.sub.n.sup.m(k).  (42)

(69) Assuming the individual coefficients C.sub.n.sup.m(ω=kc.sub.s) to be functions of the angular frequency ω, the application of the inverse Fourier transform (denoted by .sup.−1(.Math.)) provides time domain functions

(70) $\begin{array}{cc}{c}_n^m\left(t\right)={ℱ}_t^{-1}\left(C_n^m\left(\omega /c_s\right)\right)=\frac{1}{2\pi }{\int }_{-\infty }^{\infty }{C}_n^m\left(\frac{\omega }{c_s}\right)e^{i\omega t}d\omega & \left(43\right)\end{array}$

(71) for each order n and degree m, which can be collected in a single vector c(t) by
c(t)=[c.sub.0.sup.0(t)c.sub.1.sup.−1(t)c.sub.1.sup.0(t)c.sub.1.sup.1(t)c.sub.2.sup.−2(t)c.sub.2.sup.−1(t)c.sub.2.sup.0(t)c.sub.2.sup.1(t)c.sub.2.sup.2(t) . . . c.sub.N.sup.N−1(t)c.sub.N.sup.N(t)].sup.T.  (44)

(72) The position index of a time domain function c.sub.n.sup.m(t) within the vector c(t) is given by n(n+1)+1+m. The overall number of elements in vector c(t) is given by O=(N+1).sup.2.

(73) The final Ambisonics format provides the sampled version of c(t) using a sampling frequency f.sub.s as
={c(T.sub.s),c(2T.sub.s),c(3T.sub.s),c(4T.sub.s), . . . }  (45)

(74) where T.sub.s=1/f.sub.s denotes the sampling period. The elements of c(lT.sub.s) are here referred to as Ambisonics coefficients. The time domain signals c.sub.n.sup.m(t) and hence the Ambisonics coefficients are real-valued.

(75) C.1 Definition of Real-Valued Spherical Harmonics

(76) The real-valued spherical harmonics S.sub.n.sup.m(θ,φ) are given by

(77) $\begin{array}{cc}{S}_n^m\left(\theta ,\varphi \right)=\sqrt{\frac{\left(2n+1\right)}{4\pi }\frac{\left(n-.Math.m.Math.\right)!}{\left(n+.Math.m.Math.\right)!}}{P}_{n,.Math.m.Math.}\left(\cos \theta \right){\mathrm{trg}}_m\left(\varphi \right)& \left(46\right)\\ \mathrm{with}{\mathrm{trg}}_m\left(\varphi \right)=\{\begin{array}{cc}\sqrt{2}\cos \left(m\varphi \right)& m>0\\ 1& m=0\\ -\sqrt{2}\sin \left(m\varphi \right)& m<0\end{array}.& \left(47\right)\end{array}$

(78) The associated Legendre functions P.sub.n,m(x) are defined as

(79) $\begin{array}{cc}{P}_{n,m}\left(x\right)={\left(1-x^2\right)}^{\frac{m}{2}}\frac{d^m}{{\mathrm{dx}}^m}{P}_n\left(x\right),m\ge 0& \left(48\right)\end{array}$

(80) with the Legendre polynomial P.sub.n(x) and, unlike in the above-mentioned Williams article, without the Condon-Shortley phase term (−1).sup.m.

(81) C.2 Spatial Resolution of Higher Order Ambisonics

(82) A general plane wave function x(t) arriving from a direction Ω.sub.0=(θ.sub.0,φ.sub.0).sup.T is represented in HOA by
c.sub.n.sup.m(t)=x(t)S.sub.n.sup.m(Ω.sub.0), 0≦n≦N,|m|≦n.  (49)

(83) The corresponding spatial density of plane wave amplitudes c(t, Ω):=.sub.t.sup.−1(C(ω,Ω)) is given by

(84) $\begin{array}{cc}c\left(t,\Omega \right)=\underset{n=0}{\overset{N}{.Math.}}\underset{m=-n}{\overset{n}{.Math.}}{c}_n^m\left(t\right)S_n^m\left(\Omega \right)& \left(50\right)\\ =x\left(t\right)\frac{\left[\underset{n=0}{\overset{N}{.Math.}}\underset{m=-n}{\overset{n}{.Math.}}{S}_n^m\left({\Omega }_0\right)S_n^m\left(\Omega \right)\right]}{v_N\left(\Theta \right)}.& \left(51\right)\end{array}$

(85) It can be seen from equation (51) that it is a product of the general plane wave function x(t) and of a spatial dispersion function υ.sub.N(Θ), which can be shown to only depend on the angle Θ between Ω and Ω.sub.0 having the property
cos Θ=cos θ cos θ.sub.0+cos(φ−φ.sub.0)sin θ sin θ.sub.0.  (52)

(86) As expected, in the limit of an infinite order, i.e., N.fwdarw.∞, the spatial dispersion function turns into a Dirac delta δ(.Math.),

(87) $\begin{array}{cc}i.e.\underset{N.fwdarw.\infty }{\lim }{v}_N\left(\Theta \right)=\frac{\delta \left(\Theta \right)}{2\pi }.& \left(53\right)\end{array}$

(88) However, in the case of a finite order N, the contribution of the general plane wave from direction ω.sub.0 is smeared to neighbouring directions, where the extent of the blurring decreases with an increasing order. A plot of the normalised function υ.sub.N(Θ) for different values of N is shown in FIG. 5.

(89) It should be pointed out that for any direction Ω the time domain behaviour of the spatial density of plane wave amplitudes is a multiple of its behaviour at any other direction. In particular, the functions c(t,Ω.sub.1) and c(t,Ω.sub.2) for some fixed directions Ω.sub.1 and Ω.sub.2 are highly correlated with each other with respect to time t.

(90) C.3 Spherical Harmonic Transform

(91) If the spatial density of plane wave amplitudes is discretised at a number of O spatial directions Ω.sub.o, 1≦o≦O, which are nearly uniformly distributed on the unit sphere, O directional signals c(t,Ω) are obtained. Collecting these signals into a vector as
c.sub.SPAT(t):=[c(t,Ω.sub.1) . . . c(t,Ω.sub.0)].sup.T,  (54)

(92) by using equation (50) it can be verified that this vector can be computed from the continuous Ambisonics representation d(t) defined in equation (44) by a simple matrix multiplication as
c.sub.SPAT(t)=Ψ.sup.Hc(t),  (55)

(93) where (.Math.) indicates the joint transposition and conjugation, and Ψ denotes a mode-matrix defined by
Ψ:=[S.sub.1 . . . S.sub.0]  (56)
with
S.sub.0:=[S.sub.0.sup.0(Ω.sub.0)S.sub.1.sup.−1(Ω.sub.0)S.sub.1.sup.0(Ω.sub.0)S.sub.1.sup.1(Ω.sub.0) . . . S.sub.N.sup.N−1(Ω.sub.0)S.sub.N.sup.N(Ω.sub.0)].  (57)

(94) Because the directions Ω.sub.0 are nearly uniformly distributed on the unit sphere, the mode matrix is invertible in general. Hence, the continuous Ambisonics representation can be computed from the directional signals c(t,Ω.sub.0) by
c(t)=Ψ.sup.−Hc.sub.SPAT(t).  (58)

(95) Both equations constitute a transform and an inverse transform between the Ambisonics representation and the spatial domain. These transforms are here called the Spherical Harmonic Transform and the inverse Spherical Harmonic Transform.

(96) It should be noted that since the directions Ω.sub.0 are nearly uniformly distributed on the unit sphere, the approximation
Ψ.sup.H≈Ψ.sup.−1  (59)

(97) is available, which justifies the use of Ψ.sup.−1 instead of Ψ.sup.H in equation (55).

(98) Advantageously, all the mentioned relations are valid for the discrete-time domain, too.

(99) The inventive processing can be carried out by a single processor or electronic circuit, or by several processors or electronic circuits operating in parallel and/or operating on different parts of the inventive processing.