The present invention relates to coding and decoding digital signals in applications that transmit or store multimedia signals such as audio (speech and/or sound) signals or video signals.
To offer mobility and continuity, modern and innovative multimedia communication services must be able to function under a wide variety of conditions. The dynamism of the multimedia communication sector and the heterogeneous nature of networks, access points, and terminals have generated a proliferation of compression formats.
The present invention relates to optimization of the “multiple coding” techniques used when a digital signal or a portion of a digital signal is coded using more than one coding technique. The multiple coding may be simultaneous (effected in a single pass) or non-simultaneous. The processing may be applied to the same signal or to different versions derived from the same signal (for example with different bandwidths). Thus, “multiple coding” is distinguished from “transcoding”, in which each coder compresses a version derived from decoding the signal compressed by the preceding coder.
One example of multiple coding is coding the same content in more than one format and then transmitting it to terminals that do not support the same coding formats. In the case of real-time broadcasting, the processing must be effected simultaneously. In the case of access to a database, the coding could be effected one by one, and “offline”. In these examples, multiple coding is used to code the same signal with different formats using a plurality of coders (or possibly a plurality of bit rates or a plurality of modes of the same coder), each coder operating independently of the others.
Another use of multiple coding is encountered in coding structures in which a plurality of coders compete to code a signal segment, only one of the coders being finally selected to code that segment. That coder may be selected after processing the segment, or even later (delayed decision). This type of structure is referred to below as a “multimode coding” structure (referring to the selection of a coding “mode”). In these multimode coding structures, a plurality of coders sharing a “common past” code the same signal portion. The coding techniques used may be different or derived from a single coding structure. They will not be totally independent, however, except in the case of “memoryless” techniques. In the (routine) situation of coding techniques using recursive processing, the processing of a given signal segment depends on how the signal has been coded in the past. There is therefore some coder interdependency, when a coder has to take account in its memories of the output from another coder.
The concept of “multiple coding” and conditions for using such techniques have been introduced in the various contexts referred to above. The complexity of implementation may prove insurmountable, however.
For example, in the situation of content servers that broadcast the same content with different formats adapted to the access conditions, networks, and terminals of different clients, this operation becomes extremely complex as the number of formats required increases. In the case of real-time broadcasting, as the various formats are coded in parallel, a limitation is rapidly imposed by the resources of the system.
The second use referred to above relates to multimode coding applications that select one coder from a set of coders for each signal portion analyzed. Selection requires the definition of a criterion, the more usual criteria aiming to optimize the bit rate/distortion trade-off. The signal being analyzed over successive time segments, a plurality of codings are evaluated in each segment. The coding with the lowest bit rate for a given quality or the best quality for a given bit rate is then selected. Note that constraints other than those of bit rate and distortion may be used.
In such structures, the coding is generally selected a priori by analyzing the signal over the segment concerned (selection according to the characteristics of the signal). However, the difficulty of producing a robust classification of the signal for the purposes of this selection has led to the proposal for a posteriori selection of the optimum mode after coding all the modes, although this is achieved at the cost of high complexity.
Intermediate methods combining the above two approaches have been proposed with a view to reducing the computation cost. Such strategies are less than the optimum, however, and offer worse performance than exploring all the modes. Exploring all the modes or a major portion of the modes constitutes a multiple coding application that is potentially highly complex and not readily compatible a priori with real-time coding, for example.
At present, most multiple coding and transcoding operations take no account of interaction between formats and between the format and its content. A few multimode coding techniques have been proposed but the decision as to the mode to use is generally effected a priori, either on the signal (by classification, as in the SMV coder (selectable mode vocoder), for example, or as a function of the conditions of the network (as in adaptive multirate (AMR) coders, for example).
Various selection modes are described in the following documents, in particular decision controlled by the source and decision controlled by the network:
“An overview of variable rate speech coding for cellular networks”, Gersho, A.; Paksoy, E.; Wireless Communications, 1992. Conference Proceedings, 1992 IEEE International Conference on Selected Topics, 25-26 Jun. 1992 Page(s): 172-175;
“A variable rate speech coding algorithm for cellular networks”, Paksoy, E.; Gersho, A.; Speech Coding for Telecommunications, 1993. Proceedings, IEEE Workshop 1993, Page(s): 109-110; and
“Variable rate speech coding for multiple access wireless networks”, Paksoy E.; Gersho A.; Proceedings, 7th Mediterranean Electrotechnical Conference, 12-14 Apr. 1994 Page(s): 47-50 vol. 1.
In the case of a decision controlled by the source, the a priori decision is made on the basis of a classification of the input signal. There are many methods of classifying the input signal.
In the case of a decision controlled by the network, it is simpler to provide a multimode coder whose bit rate is selected by an external module rather than by the source. The simplest method is to produce a family of coders each of fixed bit rate but with different coders having different bit rates and to switch between those bit rates to obtain a required current mode.
Work has also been done on combining a plurality of criteria for a priori selection of the mode to be used; see in particular the following documents:
“Variable-rate for the basic speech service in UMTS” Berruto, E.; Sereno, D.; Vehicular Technology Conference, 1993 IEEE 43rd, 18-20 May 1993 Page(s): 520-523; and
“A VR-CELP codec implementation for CDMA mobile communications” Cellario, L.; Sereno, D.; Giani, M.; Blocher, P.; Hellwig, K.; Acoustics, Speech, and Signal Processing, 1994, ICASSP-94, 1994 IEEE International Conference, Volume: 1, 19-22 Apr. 1994 Page(s): I/281-I/284 vol. 1.
All multimode coding algorithms using a priori coding mode selection suffer from the same drawback, related in particular to problems with the robustness of a priori classification.
For this reason techniques have been proposed using an a posteriori decision as to the coding mode. For example, in the following document:
“Finite state CELP for variable rate speech coding” Vaseghi, S. V.; Acoustics, Speech, and Signal Processing, 1990, ICASSP-90, 1990 International Conference, 3-6 Apr. 1990 Page(s): 37-40 vol. 1,
the coder can switch between different modes by optimizing an objective quality measurement with the result that the decision is made a posteriori as a function of the characteristics of the input signal, the target signal-to-quantization noise ratio (SQNR), and the current status of the coder. A coding scheme of this kind improves quality. However, the different codings are carried out in parallel and the resulting complexity of this type of system is therefore prohibitive.
Other techniques have been proposed combining an a priori decision and closed loop improvement. In the document:
“Multimode variable bit rate speech coding: an efficient paradigm for high-quality low-rate representation of speech signal” Das, A.; DeJaco, A.; Manjunath, S.; Ananthapadmanabhan, A.; Huang, J.; Choy, E.; Acoustics, Speech, and Signal Processing, 1999. ICASSP '99 Proceedings, 1999 IEEE International Conference, Volume: 4, 15-19 Mar. 1999 Page(s): 2307-2310 vol. 4,
the proposed system effects a first selection (open loop selection) of the mode as a function of the characteristics of the signal. This decision may be effected by classification. Then, if the performance of the selected mode is not satisfactory, on the basis of an error measurement, a higher bit rate mode is applied and the operation is repeated (closed loop decision).
Similar techniques are described in the following documents:
“Variable rate speech coding for UMTS” Cellario, L.; Sereno, D.; Speech Coding for Telecommunications, 1993. Proceedings, IEEE Workshop, 1993 Page(s): 1-2.
“Phonetically-based vector excitation coding of speech at 3.6 kbps” Wang, S.; Gersho, A.; Acoustics, Speech, and Signal Processing, 1989. ICASSP-89., 1989 International Conference, 23-26 May 1989 Page(s): 49-52 vol. 1.
“A modified CS-ACELP algorithm for variable-rate speech coding robust in noisy environments” Beritelli, F.; IEEE Signal Processing Letters, Volume: 6 Issue: 2, February 1999 Page(s): 31-34.
An open loop first selection is effected after classification of the input signal (phonetic or voiced/non-voiced classification), after which a closed loop decision is made:
either over the complete coder, in which case the whole speech segment is coded again;
or over a portion of the coding, as in the above references preceded by an asterisk (*), in which case the dictionary to be used is selected by a closed loop process.
All of the work referred to above seeks to solve the problem of the complexity of the optimum mode selection by the total or partial use of an a priori selection or preselection that avoids multiple coding or reduces the number of coders to be used in parallel.
However, no prior art technique has ever been proposed that reduces coding complexity.
The present invention seeks to improve on this situation.
To this end it proposes a multiple compression coding method in which an input signal feeds in parallel a plurality of coders each including a succession of functional units with a view to compression coding of said signal by each coder.
The method of the invention includes the following preparatory steps:
a) identifying the functional units forming each coder and one or more functions implemented by each unit;
b) marking functions that are common from one coder to another; and
c) executing said common functions once and for all for at least some of the coders in a common calculation module.
In an advantageous embodiment of the invention, the above steps are executed by a software product including program instructions to this effect. In this regard, the present invention is also directed to a software product of the above kind adapted to be stored in a memory of a processor unit, in particular a computer or a mobile terminal, or in a removable memory medium adapted to cooperate with a reader of the processor unit.
The present invention is also directed to a compression coding aid system for implementing the method of the invention and including a memory adapted to store instructions of a software product of the type cited above.
Other features and advantages of the invention become apparent on reading the following detailed description and examining the appended drawings, in which:
a is a diagram of the application context of the present invention, showing a plurality of coders disposed in parallel;
b is a diagram of an application of the invention with functional units shared between a plurality of coders disposed in parallel;
c is a diagram of an application of the invention with functional units shared in multimode coding;
d is a diagram of an application of the invention to multimode trellis coding;
a is a diagram of the main functional units of a TDAC coder;
b is a diagram of the format of the bit stream coded by the
a is a diagram of the main functional units of an MPEG-1 (layer I and II) coder;
b is a diagram of the format of the bit stream coded by the
Refer first to
For simplicity, all the coders in the
Some functional units BFi are sometimes identical from one mode (or coder) to another; others differ only at the level of the layers that are quantized. Usable relations also exist when using coders from the same coding family employing similar models or calculating parameters linked physically to the signal.
The present invention aims to exploit these relations to reduce the complexity of multiple coding operations.
The invention proposes firstly to identify the functional units constituting each of the coders. The technical similarities between the coders are then exploited by considering functional units whose functions are equivalent or similar. For each of those units, the invention proposes:
to define “common” operations and to effect them once only for all the coders; and
to use calculation methods specific to each coder and in particular using the results of the aforementioned common calculations. These calculation methods produce a result that may be different from that produced by complete coding. The object is then in fact to accelerate the processing by exploiting available information supplied in particular by the common calculations. Methods like these for accelerating the calculations are used in techniques for reducing the complexity of transcoding operations, for example (known as “intelligent transcoding” techniques).
b shows the proposed solution. In the present example, the “common” operations cited above are effected once only for at least some of the coders and preferably for all the coders in an independent module MI that redistributes the results obtained to at least some of the coders or preferably to all the coders. It is therefore a question of sharing the results obtained between at least some of the coders C0 to CN (this is referred to below as “mutualization”). An independent module MI of the above kind may form part of a multiple compression coding aid system as defined above.
In an advantageous variant, rather than using an external calculation module MI, the existing functional unit or units BF1 to BFn of the same coder or a plurality of separate coders are used, the coder or coders being selected in accordance with criteria explained later.
The present invention may employ a plurality of strategies which may naturally differ according to the role of the functional unit concerned.
A first strategy uses the parameters of the coder having the lowest bit rate to focus the parameter search for all the other modes.
A second strategy uses the parameters of the coder having the highest bit rate and then “downgrades” progressively to the coder having the lowest bit rate.
Of course, if preference is to be given to a particular coder, it is possible to code a signal segment using that coder and then to reach coders of higher and lower bit rate by applying the above two strategies.
Of course, criteria other than the bit rate can be used to control the search. For some functional units, for example, preference may be given to the coder whose parameters lend themselves best to efficient extraction (or analysis) and/or coding of similar parameters of the other coders, efficacy being judged according to complexity or quality or a trade-off between the two.
An independent coding module not present in the coders but enabling more efficient coding of the parameters of the functional unit concerned for all the coders may also be created.
The various implementation strategies are particularly beneficial in the case of multimode coding. In this context, shown in
In this particular case of multimode coding, a variant of the present invention represented in
A more sophisticated variant of the multimode structure based on the division into functional units described above is described next with reference to
Thus each coding mode is derived from the combination of operating modes of the functional units: functional unit 1 has N1 operating modes, functional unit 2 has N2, and so on up to unit P. The combination of the NN=N1×N2× . . . ×Np possible combinations is therefore represented by a trellis with NN branches defining, end-to-end, a complete multimode coder with NN modes. Some branches of the trellis may be eliminated a priori to define a tree having a reduced number of branches. A first particular feature of this structure is that, for a given functional unit, it provides a common calculation module for each output of the preceding functional unit. These common calculation modules carry out the same operations, but on different signals, since they come from different previous units. The common calculation modules of the same level are advantageously mutualized: the results from a given module usable by the subsequent modules are supplied to those subsequent modules. Secondly, partial selection following the processing of each functional unit advantageously enables the elimination of branches offering the lowest performance against the selected criterion. Thus the number of branches of the trellis to be evaluated may be reduced.
One advantageous application of this multimode trellis structure is as follows.
If the functional units are liable to operate at respective different bit rates using respective parameters specific to said bit rates, for a given functional unit, the path of the trellis selected is that through the functional unit with the lowest bit rate or that through the functional unit with the highest bit rate, according to the coding context, and the results obtained from the functional unit with the lowest (or highest) bit rate are adapted to the bit rates of at least some of the other functional units through a focused parameter search for at least some of the other functional units, up to the functional unit with the highest (respectively lowest) bit rate.
Alternatively, a functional unit of given bit rate is selected and at least some of the parameters specific to that functional unit are adapted progressively, by focused searching:
up to the functional unit capable of operating at the lowest bit rate; and
up to the functional unit capable of operating at the highest bit rate.
This generally reduces the complexity associated with multiple coding.
The invention applies to any compression scheme using multiple coding of multimedia content. Three embodiments are described below in the field of audio (speech and sound) compression. The first two embodiments relate to the family of transform coders, to which the following reference document relates:
“Perceptual Coding of Digital Audio”, Painter, T.; Spanias, A.; Proceedings of the IEEE, Vol. 88, No 4, April 2000.
The third embodiment relates to CELP coders, to which the following reference document relates:
“Code Excited Linear Prediction (CELP): High quality speech at very low bit rates” Schroeder M. R.; Atal B. S.; Acoustics, Speech, and Signal Processing, 1985. Proceedings. 1985 IEEE International Conference, Page(s): 937-940.
A summary of the main characteristics of these two coding families is given first.
Transform or Sub-Band Coders
These coders are based on psycho-acoustic criteria and transform blocks of the signal in the time domain to obtain a set of coefficients. The transforms are of the time-frequency type, one of the most widely used transforms being the modified discrete cosine transform (MDCT). Before the coefficients are quantized, an algorithm assigns bits so that the quantizing noise is as inaudible as possible. Bit assignment and coefficient quantization use a masking curve obtained from a psycho-acoustic model used to evaluate, for each line of the spectrum considered, a masking threshold representing the amplitude necessary for a sound at that frequency to be audible.
a unit 21 for effecting the time/frequency transform on the input digital audio signal so;
a unit 22 for determining a perceptual model from the transformed signal;
a quantizing and coding unit 23 operating on the conceptual model; and
a unit 24 for formatting the bit stream to obtain a coded audio stream stc.
Analysis by Synthesis Coders (CELP Coding)
In coders of the analysis by synthesis type, the coder uses the synthesis model of the reconstructed signal to extract the parameters modeling the signals to be coded. Those signals may be sampled at a frequency of 8 kilohertz (kHz) (300-3400 hertz (Hz) telephone band) or at higher frequency, for example at 16 kHz for broadened band coding (bandwidth from 50 Hz to 7 kHz). Depending on the application and the required quality, the compression ratio varies from 1 to 16. These coders operate at bit rates from 2 kilobits per second (kbps) to 16 kbps in the telephone band and from 6 kbps to 32 kbps in the broadened band.
In the
Decoding is much less complex than coding. The decoder can obtain the quantizing index of each parameter from the bit stream generated by the coder after demultiplexing. The signal can then be reconstructed by decoding the parameters and applying the synthesis model.
The three embodiments referred to above are described below, beginning with a transform coder of the type shown in
The first embodiment relates to a “TDAC” perceptual frequency domain coder described in particular in the published document US-2001/027393. A TDAC coder is used to code digital audio signals sampled at 16 kHz (broadened band signals).
Dynamic bit assignment (in functional unit 44) is based on a masking curve for each band calculated from the decoded and dequantized version of the spectral envelope (functional unit 42). This makes bit assignment by the coder and the decoder compatible. The normalized MDCT coefficients in each band are then quantized (in functional unit 45) by vector quantizers using size-interleaved dictionaries consisting of a union of type II permutation codes. Finally, referring to
This coder is able to operate at several bit rates and it is therefore proposed to produce a multiple bit rate coder, for example a coder offering bits rates of 16, 24 and 32 kbps. In this coding scheme, the following functional units may be pooled between the various modes:
MDCT (functional unit 41);
voicing detection (functional unit 47,
calculation, quantization and entropic coding of the spectral envelope (functional unit 43); and
calculation of a masking curve coefficient by coefficient and of a masking curve for each band (functional unit 42).
These units account for 61.5% of the complexity of the processing performed by the coding process. Their factorization is therefore of major interest in terms of reducing complexity when generating a plurality of bit streams corresponding to different bit rates.
The results from the above functional units already yield a first portion common to all the output bit streams that contain the bits carrying information on voicing, tonality and the coded spectral envelope.
In a first variant of this embodiment, it is possible to carry out the bit assignment and quantization operations for each of the output bit streams corresponding to each of the bit rates considered. These two operations are carried out in exactly the same way as is usually done in a TDAC coder.
In a second, more advanced variant, shown in
bit assignment (functional unit 44); and
coefficient quantization (functional units 45—i, see below).
In
For the bit assignment and quantization functional units, the strategy employed consists in exploiting the results from the bit assignment and quantization functional units obtained for the bit stream (0), at the lowest bit rate D0, to accelerate the operation of the corresponding two functional units for the K−1 other bit streams (k) (1≦k<K). A multiple bit rate coding scheme that uses a bit assignment functional unit for each bit stream (with no factorization for that unit) but mutualizes some of the subsequent quantization operations may also be considered.
The multiple coding techniques described above are advantageously based on intelligent transcoding to reduce the bit rate of the coded audio stream, generally in a node of the network.
The bit streams k (0≦k<K) are classified in increasing bit rate order (D0<D1< . . . <DK−1) below. Thus bit stream 0 corresponds to the lowest bit rate.
Bit Assignment
Bit assignment in the TDAC coder is effected in two phases. Firstly, the number of bits to assign to each band is calculated, preferably using the following equation:
is a constant,
Each of the values obtained is rounded off to the nearest natural integer. If the total bit rate assigned is not exactly equal to that available, a second phase effects an adjustment, preferably by means of a succession of iterative operations based on a perceptual criterion that adds bits to or removes bits from the bands.
Accordingly, if the total number of bits distributed is less than that available, bits are added to the bands showing the greatest perceptual improvement, as measured by the variation of the noise-to-mask ratio between the initial and final band assignments. The bit rate is increased for the band showing the greatest variation. In the contrary situation where the total number of bits distributed is greater than that available, the extraction of bits from the bands is the dual of the above procedure.
In the multiple bit rate coding scheme corresponding to the TDAC coder, it is possible to factorize certain operations for the assignment of bits. Thus the first phase of determination using the above equation may be effected once only based on the lowest bit rate D0. The phase of adjustment by adding bits may then be effected continuously. Once the total number of bits distributed reaches the number corresponding to a bit rate of a bit stream k (k=1, 2 . . . , K−1), the current distribution is considered to be that used for quantizing normalized coefficient vectors for each band of that bit stream.
Coefficient Quantization
For coefficient quantization, the TDAC coder uses vector quantization employing size-interleaved dictionaries consisting of a union of type II permutation codes. This type of quantization is applied to each of the vectors of the MDCT coefficients over the band. This kind of vector is normalized beforehand using the dequantized value of the spectral envelope over that band. The following notation is used:
C(bi,di) is the dictionary corresponding to the number of bits bi and the dimension di;
N(bi,di) is the number of elements in that dictionary;
CL(bi,di) is the set of its leaders; and
NL(bi,di) is the number of leaders.
The quantization result for each band i of the frame is a code word mi transmitted in the bit stream. It represents the index of the quantized vector in the dictionary calculated from the following information:
the number Li in the set CL(bi,di) of the leaders of the dictionary C(bi,di) of the quantized leader vector {tilde over (Y)}q(i) nearest a current leader {tilde over (Y)}(i);
the rank ri of Yq(i) in the class of the leader {tilde over (Y)}q(i); and
the combination of signs signq(i) to be applied to Yq(i) (or to {tilde over (Y)}q(i)).
The following notation is used:
Y(i) is the vector of the absolute values of the normalized coefficients of the band i;
sign(i) is the vector of the signs of the normalized coefficients of the band i;
{tilde over (Y)}(i) is the leader vector of the vector Y(i) cited above obtained by ordering its components in decreasing order (the corresponding permutation is denoted perm(i)); and
Yq(i) is the quantized vector of Y(i) (or “the nearest neighbor” of Y(i) in the dictionary C(bi,di)).
Below, the notation α(k) with an exponent k indicates the parameter used in the processing effected to obtain the bit stream of the coder k. Parameters without this exponent are calculated once and for all for the bit stream 0. They are independent of the bit rate (or mode) concerned.
The “interleaving” property of the dictionaries referred to above is expressed as follows:
C(bi(0),di)⊂ . . . ⊂C(bi(k−1),di)⊂C(bi(k),di) . . . ⊂C(bi(k−1),di)
also with:
CL(bi(0),di)⊂ . . . ⊂CL(bi(k−1),di)⊂CL(bi(k),di) . . . ⊂CL(bi(k−1),di)
CL(bi(k),di)\CL(bi(k−1),di) is the complement of CL(bi(k−1),di) in CL(bi(k),di). Its cardinal is equal to NL(bi(k),di)−NL(bi(k−1),di).
The code words mi(k) (with 0≦k<K ), which are the results of quantizing the vector of the coefficients of the band i for each of the bit streams k, are obtained as follows.
For the bit stream k=0, the quantizing operation is effected conventionally, as is usual in the TDAC coder. It produces the parameters signq(0)(i), Li(0) and ri(0) used to construct the code word mi(0). The vectors {tilde over (Y)}(i) and sign(i) are also determined in this step. They are stored in memory, together with the corresponding permutation perm(i), to be used, if necessary, in subsequent steps relating to the other bit streams.
For the bit streams 1≦k<K, an incremental approach is adopted, from k=1 to k=K−1, preferably using the following steps:
If (bi(k)=bi(k−1)) then:
If not, i.e. if (bi(k)>bi(k−1)):
The MPEG-1 Layer I&II coder shown in
Starting from this coding scheme, in one application of the invention a multiple bit rate coder may be constructed by pooling the following functional units (see
Bank of analysis filters 61;
Determination of scaling factors 67;
FFT calculation 65; and
Masking threshold determination 64 using a psycho-acoustic model.
The functional units 64 and 65 already supply the signal-to-mask ratios (arrows SMR in
In the embodiment shown in
Bit Assignment
In the MPEG-1 Layer I&II coder, bit assignment is preferably effected by a succession of interactive steps, as follows:
Step 0: Initialize to zero the number of bits bi for each of the sub-bands i (0≦i<M).
Step 1: Update the distortion function NMR(i) (noise-to-mask ratio) over each of the sub-bands NMR(i)=SMR(i)−SNR(bi), where SNR(bi) is the signal-to-noise ratio corresponding to the quantizer having a number of bits bi and SMR(i) is the signal-to-mask ratio supplied by the psycho-acoustic model.
Step 2: Increment the number of bits bi
where ε is a positive integer value depending on the band, generally taken as equal to 1.
Steps 1 and 2 are iterated until the total number of bits available, corresponding to the operational bit rate, has been distributed. The result of this is a bit distribution vector (b0,b1, . . . ,bM−1).
In the multiple bit rate coding scheme, these steps are pooled with a few other modifications, in particular:
the output of the functional unit consisting of K bit distribution vectors (b0(k),b1(k), . . . ,b(M−1)(k) (0≦k≦K−1) , a vector (b0(k),b1(k), . . . ,b(M−1)(k) is obtained when the total number of bits available corresponding to the bit rate Dk of the bit stream k has been distributed, in the iteration of steps 1 and 2; and
the iteration of steps 1 and 2 is stopped when the total number of bits available corresponding to the highest bit rate DK−1 has been totally distributed (the bit streams are in order of increasing bit rate).
Note that the bit distribution vectors are obtained successively from k=0 up to k=K−1. The K outputs of the bit assignment functional unit therefore feed the quantization functional units for each of the bit streams at the given bit rate.
The final embodiment concerns coding multimode speech using the a posteriori decision 3GPP NB-AMR (Narrow-Band Adaptive Multi-Rate) coder, which is a telephone band speech coder conforming to the 3GPP standard. This coder belongs to the well-known family of CELP coders, the theory of which is described briefly above, and has eight modes (or bit rates) from 12.2 kbps to 4.75 kbps, all based on the algebraic code excited linear prediction (ACELP) technique.
In a first variant, only mutualization of identical functional units is exploited (the results of the four codings are then identical to those of the four codings in parallel).
In a second variant, the complexity is reduced further. The calculations of functional units that are not identical for certain modes are accelerated by exploiting those of another mode or of a common processing module (see below). The results with the four codings mutualized in this way are then different from those of the four codings in parallel.
In a further variant, the functional units of these four modes are used for multimode trellis coding, as described above with reference to
The four modes (7.4; 6.7; 5.9; 5.15) of the 3GPP NB-AMR coder are described briefly next.
The 3GPP NB-AMR coder operates on a speech signal band-limited to 3.4 kHz, sampled at 8 kHz and divided into frames of 20 ms (160 samples). Each frame contains four 5 ms subframes (40 samples) grouped two by two into 10 ms “supersubframes” (80 samples). For all the modes, the same types of parameters are extracted from the signal but with variants in terms of the modeling and/or quantization of the parameters. In the NB-AMR coder, five types of parameters are analyzed and coded. The line spectral pair (LSP) parameters are processed once per frame for all modes except the 12.2 mode (and thus once per supersubframe). The other parameters (in particular the LTP delay, adaptive excitation gain, fixed excitation and fixed excitation gain) are processed once per subframe.
The four modes considered here (7.4; 6.7; 5.9; 5.15) differ essentially in terms of the quantization of their parameters. The bit assignment of these four modes is summarized in table 1 below.
These four modes (7.4; 6.7; 5.9; 5.15) of the NB-AMR coder use exactly the same modules, for example preprocessing, linear prediction coefficient analysis and weighted signal calculation modules. The preprocessing of the signal is low-pass filtering with a cut-off frequency of 80 Hz to eliminate DC components combined with division by two of the input signals to prevent overflows. The LPC analysis comprises windowing submodules, autocorrelation calculation submodules, Levinson-Durbin algorithm implementation submodules, A(z)→LSP transform submodules, submodules for calculating LSPi non-quantized parameters for each subframe (i=0, . . . , 3) by interpolation between the LSP of the past frame and those of the current frame, and inverse LSPi→Ai(z) transform submodules.
Calculating the weighted speech signal consists in filtering by the perceptual weighting filter (Wi(z)=Ai(z/γ1)/Ai(z/γ2) where Ai(z) is the non-quantized filter of the subframe of index i, γ1=0.94 and γ2=0.6).
Other functional units are the same for only three of the modes (7.4; 6.7; 5.9). For example, the open loop LTP delay search effected on the weighted signal once per supersubframe for these three modes. For the 5.15 mode, it is effected only once per frame, however.
Similarly, if the four modes used first order predictive weighted vectorial MA (moving average) quantization of with suppressed average and Cartesian product of the LSP parameters in the normalized frequency domain, the LSP parameters of the 5.15 kbps mode are quantized on 23 bits and those of the other three modes on 26 bits. Following transformation into the normalized frequency domain, the “split VQ” vector quantization per Cartesian product of the LSP parameters splits the 10 LSP parameters into three subvectors of size 3, 3 and 4. The first subvector composed of the first three LSP is quantized on 8 bits using the same dictionary for the four modes. The second subvector composed of the next three LSP is quantized for the three high bit rate modes using a dictionary of size 512 (9 bits) and for the 5.15 mode using half of that dictionary (one vector in two). The third and final subvector composed of the last four LSP is quantized for the three high bit rate modes using a dictionary of size 512 (9 bits) and for the lower bit rate mode using a dictionary of size 128 (7 bits). The transformation into the normalized frequency domain, the calculation of the weight of the quadratic error criterion and the moving average (MA) prediction of the LSP residue to be quantized are exactly the same for the four modes. Because the three high bit rate modes use the same dictionaries to quantize the LSP, they can share, in addition to the same vector quantization module, the inverse transform (to revert from the normalized frequency domain to the cosine domain), as well as the calculation of the LSPQi quantized for each subframe (i=0, . . . , 3) by interpolation between the quantized LSP of the past frame and those of the current frame, and finally the inverse transform LSPQi→AQi(z)
Adaptive and fixed excitation closed loop searches are effected sequentially and necessitate calculation beforehand of the impulse response of the weighted synthesis filter and then of target signals. The impulse response (Ai(z/γ1)/[AQi(z)Ai(z/γ2)]) of the weighted synthesis filter is exactly the same for the three high bit rate modes (7.4; 6.7; 5.9). For each subframe, the calculation of the target signal for adaptive excitation depends on the weighted signal (independently of the mode), the quantized filter AQi(z) (which is exactly the same for the three modes) and the past of the subframe (which is different for each subframe other than the first subframe). For each subframe, the target signal for fixed excitation is obtained by subtracting from the preceding target signal the contribution of the filtered adaptive excitation of that subframe (which is different from one mode to the other except for the first subframe of the first three modes).
Three adaptive dictionaries are used. The first dictionary, used for the even subframes (i=0 and 2) of the 7.4; 6.7; 5.9 modes and for the first subframe of the 5.15 mode, includes 256 fractional absolute delays of ⅓ resolution in the range [19+ 1/3.84+⅔] and of entire resolution in the range [85.143]. Searching in this absolute delay dictionary is focused around the delay found in open loop mode (interval of ±5 for the 5.15 mode or ±3 for the other modes). For the first subframe of the 7.4; 6.7; 5.9 modes, the target signal and the open loop delay being identical, the result of the closed loop search is also identical. The other two dictionaries are of differential type and are used to code the difference between the current delay and the entire delay Ti−1 closest to the fractional delay of the preceding subframe. The first differential dictionary on five bits, used for the odd subframes of the 7.4 mode, is of ⅓ resolution about the entire delay Ti−1 in the range [Ti−1−5 +⅔, Ti−1+4 +⅔]. The second differential dictionary on four bits, which is included in the first differential dictionary, is used for the odd subframes of the 6.7 and 5.9 modes and for the last three subframes of the 5.15 mode. This second dictionary is of entire resolution about the entire delay Ti−1 in the range [Ti−1−5, Ti−1+4] plus a resolution of ⅓ in the range [Ti−1−1+⅔, Ti−1+⅔].
The fixed dictionaries belong to the well-known family of ACELP dictionaries. The structure of an ACELP directory is based on the interleaved single-pulse permutation (ISPP) concept, which consists in dividing the set of L positions into K interleaved tracks, the N pulses being located in certain predefined tracks. The 7.4, 6.7, 5.9 and 5.15 modes use the same division of the 40 samples of a subframe into five interlaced tracks of length 8, as shown in Table 2a. Table 2b shows, for the 7.4, 6.7 and 5.9 modes, the bit rate of the dictionary, the number of pulses and their distribution in the tracks. The distributions of the two pulses of the 5.15 mode of the ACELP dictionary with nine bits is even more constrained.
The adaptive and fixed excitation gains are quantized on seven or six bits (with MA prediction applied to the fixed excitation gain) by conjoint vector quantization minimizing the CELP criterion.
Multimode Coding With a Posteriori Decision Exploiting Only Mutualization of Identical Functional Units
An a posteriori decision multimode coder may be based on the above coding scheme, pooling the functional units indicated below.
Referring to
pre-processing (functional unit 81);
analyzing the linear prediction coefficients (windowing and calculating the autocorrelations 82, executing the Levinson-Durbin algorithm 83; A(z)→LSP transform 84, interpolating the LSP and inverse transformation 862);
calculating the weighted input signal 87;
transforming the LSP parameters into the normalized frequency domain, calculating the weight of the quadratic error criterion for vector quantization of the LSP, MA prediction of the LSP residue, vector quantization of the first three LSP (in the functional unit 85).
Thus the cumulative complexity for all these units is divided by four.
For the three highest bit rate modes (7.4, 6.7 and 5.9), there are effected:
vector quantization of the last seven LSP (once per frame) (in functional unit 85 in
open loop LTP delay search (twice per frame) (functional unit 88);
quantized LSP interpolation (861) and inverse transformation to the filters AQi (for each subframe); and
calculation of the impulse response 89 of the weighted synthesis filter (for each subframe).
For these units, the calculations are no longer effected four times but only twice, once for the three highest bit rate modes and once for the low bit rate mode. Their complexity is therefore divided by two.
For the three highest bit rate modes, it is also possible to mutualize for the first subframe the calculation of the target signals for the fixed excitation (functional unit 91 in
Advanced a Posteriori Decision Multimode Coding
Non-identical functional units can be accelerated by exploiting those of another mode or a common processing module. Depending on the constraints of the application (in terms of quality and/or complexity), different variants may be used. A few examples are described below. It is also possible to rely on intelligent transcoding techniques between CELP coders.
Vector Quantization of the Second LSP Subvector
As in the TDAC coder embodiment, interleaving certain dictionaries can accelerate the calculations. Accordingly, as the dictionary of the second LSP subvector of the 5.15 mode is included in that of the other three modes, the quantization of that subvector Y by the four modes can be advantageously combined:
Step 1: Search for nearest neighbor Y1 in the smallest dictionary (corresponding to half the large dictionary)
Y1 quantizes Y for the 5.15 mode
Step 2: Search for the nearest neighbor Yh in the complement in the large dictionary (i.e. in the other half of the dictionary)
Step 3: Test if the nearest neighbor of Y in the 9-bit dictionary is Y1 (“Flag=0”) or Yh (“Flag=1”)
“Flag=0”: Y1 also quantizes Y for the 7.4, 6.7 and 5.9 modes
“Flag=1”: Yh quantizes Y for the 7.4, 6.7 and 5.9 modes
This embodiment gives an identical result to non-optimized multimode coding. If quantization complexity is to be reduced further, we can stop at step 1 and take Y1 as the quantized vector for the high bit rate modes if that vector is deemed sufficiently close to Y. This simplification can therefore yield a result different from an exhaustive search.
Open Loop LTP Search Acceleration
The 5.15 mode open loop LTP delay search can use search results for the other modes. If the two open loop delays found over the two supersubframes are sufficiently close to allow differential coding, the 5.15 mode open loop search is not effected. The results of the higher modes are used instead. If not, the options are:
to effect the standard search; or
to focus the open loop search on the whole of the frame around the two open loop delays found by the higher modes.
Conversely, the 5.15 mode open loop delay search may also be effected first and the two higher mode open loop delay searches focused around the value determined by the 5.15 mode.
In a third and more advanced embodiment shown in
There are therefore P=4 functional units and 2×3×4×2=48 possible combinations. In this particular embodiment the high bit rate of functional unit 2 (LTP bit rate 26 bits/frame) is not considered. Other choices are possible, of course.
The multiple bit rate coder obtained in this way has a high granularity in terms of bit rates with 32 possible modes (see Table 3b). However, the resulting coder cannot interwork with the NB-AMR coder cited above. In Table 3b, the modes corresponding to the 5.15, 5.90 and 6.70 bit rates of the NB-AMR coder are shown in bold, the exclusion of the highest bit rate of the functional unit LTP eliminating the 7.40 bit rate.
This coder having 32 possible bit rates, five bits are necessary for identifying the mode used. As in the previous variant, functional units are mutualized. Different coding strategies are applied to the different functional units.
For example, for functional unit 1 including LSP quantization, preference is given to the low bit rate, as mentioned above, and as follows:
the first subvector made up of the first three LSP is quantized on 8 bits using the same dictionary for the two bit rates associated with this functional unit;
the second subvector made up of the next three LSP is quantized on 8 bits using the dictionary with the lowest bit rate. That dictionary corresponding to half the higher bit rate dictionary, the search is effected in the other half of the dictionary only if the distance between the three LSP and the chosen element in the dictionary exceeds a certain threshold; and
the third and final subvector made up of the last four LSP is quantized using a dictionary of size 512 (9 bits) and a dictionary of size 128 (7 bits).
On the other hand, as mentioned above in relation to the second variant (corresponding to multimode coding with advanced a posteriori decision) the choice is made to give preference to the high bit rate for functional unit 2 (LTP delay). In the NB-AMR coder, the open loop LTP delay search is effected twice per frame for the LTP delay of 24 bits and only once per frame for that of 20 bits. The aim is to give preference to the high bit rate for this functional unit. The open loop LTP delay calculation is therefore effected in the following manner:
Two open loop delays are calculated over the two supersubframes. If they are sufficiently close to allow differential coding, the open loop search is not effected over the entire frame. The results for the two supersubframes are used instead; and
If they are not sufficiently close, an open loop search is effected over the whole of the frame, focused around the two open loop delays found beforehand. A variant reducing complexity retains only the open loop delay of the first of them.
It is possible to make a partial selection to reduce the number of combinations to be explored after certain functional units. For example, after functional unit 1 (LPC), the combinations with 26 bits can be eliminated for this block if the performance of the 23 bits mode is sufficiently close or the 23 bits mode can be eliminated if its performance is too degraded compared to the 26 bits mode.
Thus the present invention can provide an effective solution to the problem of the complexity of multiple coding by mutualizing and accelerating the calculations executed by the various coders. The coding structures can therefore be represented by means of functional units describing the processing operations effected. The functional units of the different forms of coding used in multiple coding have strong relations that the present invention exploits. Those relations are particularly strong when different codings correspond to different modes of the same structure.
Note finally that from the point of view of complexity the present invention is flexible. It is in fact possible to decide a priori on the maximum multiple coding complexity and to adapt the number of coders explored as a function of that complexity.
Number | Date | Country | Kind |
---|---|---|---|
0314490 | Dec 2003 | FR | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/FR04/03009 | 11/24/2004 | WO | 6/8/2006 |