The present invention relates to an audio/speech encoding apparatus, audio/speech decoding apparatus and audio/speech encoding and decoding methods using Huffman coding.
In signal compression, Huffman coding is widely used to encode an input signal utilizing a variable-length (VL) code table (Huffman table). Huffman coding is more efficient than fixed-length (FL) coding for the input signal which has a statistical distribution that is not uniform.
In Huffman coding, the Huffman table is derived in a particular way based on the estimated probability of occurrence for each possible value of the input signal. During encoding, each input signal value is mapped to a particular variable length code in the Huffman table.
By encoding signal values that are statistically more likely to occur using relatively short VL codes (using relatively few bits), and conversely encoding signal values that are statistically infrequently to occur using relatively long VL codes (using relatively more bits), the total number of bits used to encode the input signal can be reduced.
However, in some applications, such as audio signal encoding, the signal statistics may vary significantly from one set of audio signal to another set of audio signal. And even within the same set of audio signal.
If the statistics of the audio signal varies drastically from the statistics of the predefined Huffman table, the encoding of the signal can not be optimally done. And it happens that, to encode the audio signal which has different statistics, the bits consumption by Huffman coding is much more than the bits consumption by fixed length coding.
One possible solution is to include both the Huffman coding and fixed length coding in the encoding, and the encoding method which consumes fewer bits are selected. One flag signal is transmitted to decoder side to indicate which coding method is selected in encoder. This solution is utilized in a newly standardized ITU-T speech codec 0.719.
The solution solves the problem for some very extreme sequences in which the Huffman coding consumes more bits than the fixed length coding. But for other input signals which have different statistics from the Huffman table but still select the Huffman coding, it is still not optimal.
In ITU-T standardized speech codec G719, Huffman coding is used in encoding of the norm factors' quantization indices.
The structure of G.719 is illustrated in
At encoder side, the input signal sampled at 48 kHz is processed through a transient detector (101). Depending on the detection of a transient, a high frequency resolution or a low frequency resolution transform (102) is applied on the input signal frame. The obtained spectral coefficients are grouped into bands of unequal lengths. The noun of each band is estimated (103) and resulting spectral envelope consisting of the norms of all bands is quantized and encoded (104). The coefficients are then normalized by the quantized norms (105). The quantized norms are further adjusted (106) based on adaptive spectral weighting and used as input for bit allocation (107). The normalized spectral coefficients are lattice-vector quantized and encoded (108) based on the allocated bits for each frequency band. The level of the non-coded spectral coefficients is estimated, coded (109) and transmitted to the decoder. Huffman encoding is applied to quantization indices for both the coded spectral coefficients as well as the encoded norms.
At decoder side, the transient flag is first decoded which indicates the frame configuration, i.e., stationary or transient. The spectral envelope is decoded and the same, bit-exact, norm adjustments and bit-allocation algorithms are used at the decoder to recompute the bit-allocation which is essential for decoding quantization indices of the normalized transform coefficients. After de-quantization (112), low frequency non-coded spectral coefficients (allocated zero bits) are regenerated by using a spectral-fill codebook built from the received spectral coefficients (spectral coefficients with non-zero bit allocation) (113). Noise level adjustment index is used to adjust the level of the regenerated coefficients. High frequency non-coded spectral coefficients are regenerated using bandwidth extension. The decoded spectral coefficients and regenerated spectral coefficients are mixed and lead to normalized spectrum. The decoded spectral envelope is applied leading to the decoded full-band spectrum (114). Finally, the inverse transform (115) is applied to recover the time-domain decoded signal. This is performed by applying either the inverse modified discrete cosine transform for stationary modes, or the inverse of the higher temporal resolution transform for transient mode.
In encoder (104), the norm factors of the spectral sub bands are scalar quantized with a uniform logarithmic scalar quantizer with 40 steps of 3 dB. The codebook entries of the logarithmic quantizer are shown in
The encoding of quantization indices for norm factors is illustrated in
[1]
Diff_index(n)=Index(n)−Index(n−1)+15 for nε[1,43] (Equation 1)
And the differential indices are encoded by two possible methods, fixed length coding (305) and Huffman coding (306). The Huffman table for the differential indices is shown in
However, for an audio input signal, there is a physical phenomenon named as auditory masking. Auditory masking occurs when the perception of one sound is affected by the presence of another sound. As example, if there are two signals with similar frequencies existing at the same time: one powerful spike at 1 kHz and one lower-level tone at 1.1 kHz, the lower-level tone at 1.1 kHz will be masked (inaudible) due to existence of the powerful spike at 1 kHz.
The sound pressure level needed to make the sound perceptible in the presence of another sound (masker), is defined as masking threshold in audio encoding. The masking threshold depends upon the frequency, the sound pressure level of the masker. If the two sounds have similar frequency, the masking effect is large, and the masking threshold is also large. If the masker has large sound pressure level, it has strong masking effect on the other sound, and the masking threshold is also large.
According to the auditory masking theory above, if one sub band has very large energy, it would have large masking effect on other sub bands, especially on its neighboring sub bands. Then the masking threshold for other sub bands, especially the neighboring sub band, is large.
If the sound component in the neighboring sub band has small quantization errors (less than the masking threshold), the degradation on sound component in this sub band is not able to be perceived by the listeners.
It is not necessary to encode the normal factor with very high resolution for this sub band as long as the quantization errors below the masked threshold.
In this invention, apparatus and methods exploring audio signal properties for generating Huffman tables and for selecting Huffman tables from a set of predefined tables during audio signal encoding are provided.
Briefly, the auditory masking properties are explored to narrow down the range of the differential indices, so that a Huffman table which have fewer code words can be designed and used for encoding. As the Huffman table has fewer code words, it is possible to design the code codes with shorter length (consumes fewer bits). By doing this, the total bits consumption to encode the differential indices can be reduced.
By adopting Huffman codes which consume fewer bits, the total bits consumption to encode the differential indices can be reduced.
The main principle of the invention is described in this section with the aid of
In the encoder illustrated in
The differential indices for the modified indices are calculated according to the equation below:
[2]
Diff_index(n)=New_index(n)−New_index(n−1)+15 for nε[1,43] (Equation 2)
The range of the differential indices for Huffman coding is identified as shown in the equation below (504).
[3]
Range=[Min(Diff_index(n),Max(Diff_index(n))] (Equation 3)
According to the value of the range, the Huffman table which is designed for the specific range among a set of predefined Huffman table is selected (505) for encoding of the differential indices (506). As example, if among all the differential indices for the input frame, the minimum value is 12, and the maximum value is 18, then the Range=[12,18]. The Huffman table designed for [12,18] are selected as the Huffman table for encoding.
The set of predefined Huffman tables are designed (detail will be explained in later part) and arranged according to the range of the differential indices. The flag signal to indicate the selected Huffman table and the coded indices are transmitted to the decoder side.
Another method for selection of Huffman table is to calculate all the bits consumption using every Huffman table, then select the Huffman table which consumes fewest bits.
As example, a set of 4 predefined Huffman tables are shown in
Comparing the Huffman code length in
In the decoder illustrated in
[4]
Diff_index(n)=Index(n)+Index(n−1)−15 for nε[1,43] (Equation 4)
As shown in
The modification of the indices can be done as below (using sub band 2 as example). As shown in
For sub bands 1 and 3, because their energies are above the masking threshold, their indices are not changed. Then the differential indices are closer to the centre. Using sub band 1 as example:
[5]
Diff_index(1)=Index(1)−Index(0)+15 for nε[1,43] (Equation 5)
[6]
New_diff_index(1)=New_index(1)−New_index(0)+15 for nε[1,43] (Equation 6)
[7]
∵New_index(1)−New_index(0)<Index(1)−Index(0)
∴New_diff_index(1)−15<Diff_index(1)−15 (Equation 7)
In this invention, the design of the Huffman table can be done offline with a large input sequence database. The process is illustrated in
The energies of the sub bands processed by the psychoacoustic modelling (1001) to derive the masked threshold Mask(n). According to the derived Mask(n), the quantization indices of the norm factors for the sub bands whose quantization errors energy are below the masking threshold are modified (1002) so that the range of the differential indices can be smaller.
The differential indices for the modified indices are calculated (1003).
The range of the differential indices for Huffman coding is identified (1004). For each value of range, all the input signal which have the same range will be gathered and the probability distribution of each value of the differential index within the range is calculated.
For each value of range, one Huffman table is designed according to the probability. Some traditional Huffman table design methods can be used here to design the Huffman table.
In this embodiment, a method which can maintain the bits saving, but to restore the differential indices to a value closer to the original value is introduced.
As shown in
If they consume same number of hits in the selected Huffman table, the modified differential indices are restored to the original differential indices. If they don't consume same number of bits, the code words in the Huffman table which is closest to the original differential indices and consumes same number of bits are selected as the restored differential indices.
The merits of this embodiment are quantization error of the norm factor can be smaller while the bits consumption is the same as the embodiment 1.
In this embodiment, a method which avoids using of the psychoacoustic model but only use some energy ratio threshold is introduced.
As shown in
[8]
Energy(n)/Energy(n−1)<Threshold && Energy(n)/Energy(n+1)<Threshold (Equation 8)
The modification of the quantization index can be done as shown in the equation below:
[9]
where,
NFNEW
NFIndex(n) means the decoded norm factor for sub band nusing the original quantization index
Energy(n−1) means the energy for sub band n−1
Energy(n) means the energy for sub band n
Energy(n+1) means the energy for sub band n+1
The merit of this embodiment is the very complex and high complexity psychoacoustic modelling can be avoided.
In this embodiment, a method which narrows down the range of the differential indices while being able to perfectly reconstruct the differential indices is introduced.
As shown in
Diff_index(n)=Index(n)−Index(n−1)+15 (Equation 10)
where,
Diff_index(n) means differential index for sub band n
Index(n) means the quantization index for sub band n
Index(n−1) means the quantization index for sub band n−1
In order to reduce the range of the differential indices, a module is implemented to modify values of some differential indices (1302).
The modification is done according to the value of the differential index for the preceding sub band and a threshold.
One way to modify the differential index (when n≧1) can be done as shown in the equation below, the first differential index would not be modified so as to achieve perfect reconstruction in decoder side:
[11]
if Diff_index(n−1)>(15+Threshold),
Diff_index_new(n)=Duff_index(n)+Diff_index(n−1)−(15+Threshold); else if Diff_index(n−1)<(15−Threshold),
Diff_index_new(n)=Duff_index(n)+Diff_index(n−1)−(15−Threshold); otherwise
Diff_index_new(n)=Diff_index(n); (Equation 11)
where,
n≧1;
Diff_index(n) means differential index for sub band n;
Diff_index(n−1) means differential index for sub band n−1;
Diff_index_new(n) means the new differential index for sub band n;
Threshold means the value to examine whether to make the modification of the differential index:
The reason why this modification can reduce the range of the differential indices is explained as following: for audio/speech signal, it is true that the energy fluctuates from one frequency band to another frequency band. However, it is observed that, there is normally no abrupt change in energy from neighboring frequency bands. The energy gradually increases or decreases from one frequency band to another frequency band. The norm factors which represent the energy also gradually changes. The norm factor quantization indices would also gradually change, and then the differential indices would vary in a small range.
The abrupt energy change happens only when some main sound components which have large energy start to show effect in the frequency band or their effect start to diminish. The norm factors which represent the energy also have abrupt change from the preceding frequency band, the norm factor quantization indices would also suddenly increase or decrease by a large value. Then it resulted in a very large or very small differential index.
As an example, assume that there is one main sound component which has large energy in frequency sub band n. While in frequency sub band (n−1) and (n+1), there is no main sound component. Then according to the Huffman table in
[12]
∵Diff_index_new(n−1)<(15−Threshold)
∴Diff_index(n−1)−(15−Threshold)<0
∵Diff_index_new(n)=Diff_index(n)+Diff_index(n−1)−(15−Threshold);
∴Diff_index_new(n)<Diff_index(n) (Equation 12)
As shown in
The way to reconstruct the differential index(when n≧1), which is corresponding to the modification in encoder, can be done as shown in the equation below, the first differential index would be directly received as it is not modified at encoder side:
[13]
if Diff_index(n−1)>(15+Threshold),
Diff_index(n)=Diff_index_new(n)−Diff_index(n−1)+(15+Threshold); else if Diff_index(n−1)<(15−Threshold),
Diff_index(n)=Diff_index_new(n)−Duff_index(n−1)+(15−Threshold); otherwise
Diff_index(n)=Diff_index_new(n); (Equation 13)
where,
Diff_index(n) means differential index for sub band n;
Diff_index(n−1) means differential index for sub band n−1;
Diff_index_new(n) means the new differential index for sub band n:
Threshold means the value to examine whether to reconstruct the differential index:
As shown in the above Equation (11) and Equation (13), whether the modification of a differential index should be done and how much it should be modified is all dependent on the differential index for preceding frequency band. If the differential index for the preceding frequency band can be perfectly reconstructed, then the current differential index can also be perfectly reconstructed.
As shown in the above Equation (11) and Equation (13), the first differential index is not modified at encoder side, it is directly received and can be perfectly reconstructed, then the second differential index can be reconstructed according to the value of the first differential index; then the third differential index, the forth differential index, and so on, by following the same procedure, all the differential indices can be perfectly reconstructed.
The merit of this embodiment is that the range of the differential indices can be reduced, while the differential indices can still be perfectly reconstructed in decoder side. Therefore, the bits efficiency can be improved while retain the bit exactness of the quantization indices.
Further, although cases have been described with the embodiments above where the present invention is configured by hardware, the present invention may be implemented by software in combination with hardware.
Each function block employed in the description of the aforementioned embodiment may typically be implemented as an LSI constituted by an integrated circuit. These may be individual chips or partially or entirely contained on a single chip. “LSI” is adopted here but this may also be referred to as “IC,” “system LSI,” “super LSI” or “ultra LSI” depending on differing extents of integration.
Further, the method of circuit integration is not limited to LSI's, and implementation using dedicated circuitry or general purpose processors is also possible. After LSI manufacture, utilization of an FPGA (Field Programmable Gate Array) or a reconfigurable processor where connections and settings of circuit cells within an LSI can be reconfigured is also possible.
Further, if integrated circuit technology comes out to replace LSI's as a result of the advancement of semiconductor technology or a derivative other technology, it is naturally also possible to carry out function block integration using this technology. Application of biotechnology is also possible.
The disclosure of Japanese Patent Applications No. 2011-94295, filed on Apr. 20, 2011 and No. 2011-133432, filed on Jun. 15, 2011, including the specification, drawings and abstract is incorporated herein by reference in its entirety.
The encoding apparatus, decoding apparatus and encoding and decoding methods according to the present invention are applicable to a wireless communication terminal apparatus, base station apparatus in a mobile communication system, tele-conference terminal apparatus, video conference terminal apparatus and voice over internet protocol (VOIP) terminal apparatus.
Number | Date | Country | Kind |
---|---|---|---|
2011-094295 | Apr 2011 | JP | national |
2011-133432 | Jun 2011 | JP | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/JP12/01701 | 3/12/2012 | WO | 00 | 9/30/2013 |