Apparatus and method for processing an audio signal using a harmonic post-filter

Information

  • Patent Grant
  • 11037580
  • Patent Number
    11,037,580
  • Date Filed
    Wednesday, February 27, 2019
    5 years ago
  • Date Issued
    Tuesday, June 15, 2021
    3 years ago
Abstract
An apparatus for processing an audio signal having associated therewith a pitch lag information and a gain information, includes a domain converter for converting a first domain representation of the audio signal into a second domain representation of the audio signal; and a harmonic post-filter for filtering the second domain representation of the audio signal, wherein the post-filter is based on a transfer function including a numerator and a denominator, wherein the numerator includes a gain value indicated by the gain information, and wherein the denominator includes an integer part of a pitch lag indicated by the pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag.
Description
BACKGROUND OF THE INVENTION

Transform-based audio codecs generally introduce inter-harmonic noise when processing harmonic audio signals, particularly at low bitrates.


This effect is further worsen when the transform-based audio codec operates at low delay, due to the worse frequency resolution and/or selectivity introduced by a shorter transform size and/or a worse window frequency response.


This inter-harmonic noise is generally perceived as a very annoying artifact, significantly reducing the performance of the transform-based audio codec when subjectively evaluated on highly tonal audio material.


Several solutions exist to improve the subjective quality of transform-based audio codecs on harmonics audio signals. All of them are based on prediction-based techniques, either in the transform-domain or in the time-domain.


Examples of transform-domain approaches are:

  • [1] H. Fuchs, “Improving MPEG Audio Coding by Backward Adaptive Linear Stereo Prediction”, 99th AES Convention, New York 1995, Preprint 4086.
  • [2] L. Yin, M. Suonio, M. Vaananen, “A New Backward Predictor for MPEG Audio Coding”, 103rd AES Convention, New York 1997, Preprint 4521
  • [3] Juha Ojanpera, Mauri Vaananen, Lin Yin, “Long Term Predictor for Transform Domain Perceptual Audio Coding”, 107th AES Convention, New York 1999, Preprint 5036.


Examples of time-domain approaches are:

  • [4] Philip J. Wilson, Harprit Chhatwal, “Adaptive transform coder having long term predictor”, U.S. Pat. No. 5,012,517, Apr. 30, 1991.
  • [5] Jeongook Song, Chang-Heon Lee, Hyen-O Oh, Hong-Goo Kang, “Harmonic Enhancement in Low Bitrate Audio Coding Using and Efficient Long-Term Predictor”, EURASIP Journal on Advances in Signal Processing 2010.
  • [6] Juin-Hwey Chen, “Pitch-based pre-filtering and post-filtering for compression of audio signals”, U.S. Pat. No. 8,738,385, May 27, 2014.


SUMMARY

According to an embodiment, an apparatus for processing an audio signal having associated therewith a pitch lag information and a gain information may have: a domain converter for converting a first domain representation of the audio signal into a second domain representation of the audio signal; and a harmonic post-filter for filtering the second domain representation of the audio signal, wherein the post-filter is based on a transfer function having a numerator and a denominator, wherein the numerator has a gain value indicated by the gain information, and wherein the denominator has an integer part of a pitch lag indicated by the pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag.


According to another embodiment, a method of processing an audio signal having associated therewith a pitch lag information and a gain information may have the steps of: converting a frequency representation of the audio signal into a time-domain representation of the audio signal; and filtering the time-domain representation of the audio signal by a harmonic post-filter, wherein the post-filter is based on a transfer function having a numerator and a denominator, wherein the numerator has a gain value indicated by the gain information, and wherein the denominator has an integer part of a pitch lag indicated by the pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag.


Another embodiment may have a system for processing an audio signal having an encoder for encoding an audio signal and a decoder having a processor, the processor having: a domain converter for converting a frequency representation of the audio signal into a time-domain representation of the audio signal; and a harmonic post-filter for filtering the time-domain representation of the audio signal, wherein the post-filter is based on a transfer function having a numerator and a denominator, wherein the numerator has a gain value indicated by a gain information, and wherein the denominator has an integer part of a pitch lag indicated by a pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag.


According to another embodiment, a method of processing an audio signal having a method of encoding an audio signal and a method of decoding may have the steps of: converting a frequency representation of the audio signal into a time-domain representation of the audio signal; and filtering the time-domain representation of the audio signal using a harmonic post-filter, wherein the post-filter is based on a transfer function having a numerator and a denominator, wherein the numerator has a gain value indicated by a gain information, and wherein the denominator has an integer part of a pitch lag indicated by a pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag.


Another embodiment may have a non-transitory digital storage medium having a computer program stored thereon to perform the method of processing an audio signal having associated therewith a pitch lag information and a gain information, having the steps of: converting a frequency representation of the audio signal into a time-domain representation of the audio signal; and filtering the time-domain representation of the audio signal by a harmonic post-filter, wherein the post-filter is based on a transfer function having a numerator and a denominator, wherein the numerator has a gain value indicated by the gain information, and wherein the denominator has an integer part of a pitch lag indicated by the pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag; when said computer program is run by a computer.


Another embodiment may have a non-transitory digital storage medium having a computer program stored thereon to perform the method of processing an audio signal having a method of encoding an audio signal and a method of decoding, having the steps of: converting a frequency representation of the audio signal into a time-domain representation of the audio signal; and filtering the time-domain representation of the audio signal using a harmonic post-filter, wherein the post-filter is based on a transfer function having a numerator and a denominator, wherein the numerator has a gain value indicated by a gain information, and wherein the denominator has an integer part of a pitch lag indicated by a pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag; when said computer program is run by a computer.


The present invention is based on the finding that the subjective quality of an audio signal can be substantially improved by using a harmonic post-filter having a transfer function comprising a numerator and a denominator. The numerator of the transfer function comprises a gain value indicated by a transmitted gain information and the denominator comprises an integer part of a pitch lag indicated by a pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag.


Hence, it is possible to remove inter-harmonic noise introduced by a typical domain-changing audio decoder as an artifact. This harmonic post-filter is particularly useful in that it relies on transmitted information, i.e., the pitch gain and the pitch lag which are available anyway in a decoder, since this information is received from a corresponding encoder via a decoder input signal. Furthermore, the post-filtering is of specific accuracy due to the fact that not only the integer part of the pitch lag is accounted for, but, in addition, the fractional part of the pitch lag is accounted for. The fractional part of the pitch lag can be particularly introduced into the post-filter via a multi-tap filter which has filter coefficients actually depending on the fractional part of the pitch lag. This filter can be implemented as an FIR filter or can also be implemented as any other filter such as an IIR filter or a different filter implementation. Any domain change such as a time to frequency change or an LPC to time change or a time to LPC change or a frequency to time change can be advantageously improved by the post-filter concept of the invention. Advantageously, however, the domain change is a frequency to time domain change.


Hence, embodiments of the present invention reduces inter-harmonic noise introduced by a transform audio codec based on a long-term predictor working in the time domain. Contrary to [04]-[6], where both pre-filter before the transform coding and a post-filter after the transform decoding are used, the present invention may apply a post-filter only.


Furthermore, it has been noticed that the pre-filter employed in [04]-[6] has the tendency to introduce instabilities in the input signal given to the transform encoder. These instabilities are due to changes in gain and/or pitch lag from frame to frame. The transform coder has difficulties in encoding such instabilities, particularly at low bitrates, and one will sometimes introduce even more noise in the decoded signal compared to a situation without any pre- or post-filter.


Advantageously, the present invention does not employ any pre-filter at all and, therefore, completely avoids the problems involved with a pre-filter.


Furthermore, the present invention relies on a post-filter that is applied on the decoded signal after transform coding. This post-filter is based on a long-term prediction filter accounting for the integer part and the fractional part of the pitch lag that reduces the inter-harmonic noise introduced by the transform audio codec.


For better robustness, the post-filter parameters pitch lag and pitch gain are estimated at the encoder-side and transmitted in the bitstream. However, in other implementations, the pitch lag and pitch gain can also be estimated on the decoder-side based on the decoded audio signal obtained by an audio decoder comprising a frequency-time converter for converting a frequency-representation of the audio signal into a time-domain representation of the audio signal.


In an embodiment, the numerator additionally comprises a multi-tap filter for a zero fractional part of the pitch lag in order to compensate for a spectral tilt introduced by the multi-tap filter in the denominator, which depends on the fractional part of the pitch lag.


Advantageously, the post-filter is configured to suppress an amount of energy between harmonics in a frame, wherein the amount of energy suppressed is smaller than 20% of a total energy of the time-domain representation in the frame.


In a further embodiment, the denominator comprises a product between the multi-tap filter and the gain value.


In a further embodiment, the filter numerator further comprises a product of a first scalar value and a second scalar value, wherein the denominator only comprises the second scalar value rather than the first scalar value. These scalar values are set to predetermined values and have values greater than 0 and lower than 1; and, additionally, the second scalar value is lower than the first scalar value. Hence, it is possible in a very efficient way to set the energy removal characteristics which are typically unwanted and to additionally set the filter strength, i.e., how strong the filter attenuates inter-harmonic artifacts in a transform-domain decoder output signal.


The apparatus further comprises, in an embodiment, a filter controller for setting at least the second scalar value depending on a bitrate so that a higher value is set for a lower bitrate and vice versa.


Furthermore, the filter controller is configured for selecting, depending on the fractional part of the pitch lag, the corresponding multi-tap filter in a signal-dependent way in order to set the harmonic post-filter signal-adaptively, i.e., dependent on the actually provided fractional part value of the pitch lag.





BRIEF DESCRIPTION OF THE DRAWINGS

Embodiments of the present invention will be detailed subsequently referring to the appended drawings, in which:



FIG. 1 illustrates an embodiment of an inventive apparatus for processing an audio signal;



FIG. 2 illustrates an implementation of the harmonic post-filter represented as transfer functions in the z domain;



FIG. 3 illustrates a further embodiment for the harmonic post-filter represented by a transfer function in the z domain;



FIG. 4 illustrates an implementation of an encoder for generating an encoded signal to be decoded by a transform-domain audio decoder illustrated in FIG. 1;



FIG. 5 illustrates an implementation of the multi-tap filter as an FIR filter controlled by a filter controller;



FIG. 6 illustrates a cooperation between the filter controller and a memory having pre-stored tap weights depending on the fractional part;



FIG. 7a illustrates a frequency response of a filter having a zero α value.



FIG. 7b illustrates a frequency response of a harmonic post-filter having an α value equal to 1;



FIG. 7c illustrates a frequency response of a harmonic post-filter having an α value of 0.8;



FIG. 8a illustrates an embodiment of a harmonic post-filter having a β value equal to 0.4; and



FIG. 8b illustrates a frequency response of a harmonic post-filter having a β value of 0.2.





DETAILED DESCRIPTION OF THE INVENTION


FIG. 1 illustrates an apparatus for processing an audio signal having associated therewith a pitch lag information and a gain information. This gain information can be transmitted to a decoder 100 via a decoder input 102 receiving an encoded signal or, alternatively, this information can be calculated in the decoder itself, when this information is not available. However, for a more robust operation, it is advantageous to calculate the pitch lag information and the pitch gain information on the encoder-side.


The decoder 100 comprises e.g. a frequency-time converter for converting a frequency-time representation of the audio signal into a time-domain representation of the audio signal. Thus, the decoder is not a pure time-domain speech codec, but comprises a pure transform domain decoder or a mixed transform domain decoder or any other coder operating in a domain different from a time domain. Furthermore, it is advantageous that the second domain is the time domain.


The apparatus furthermore comprises a harmonic post-filter 104 for filtering the time-domain representation of the audio signal, and this harmonic post-filter is based on a transfer function comprising a numerator and a denominator. Particularly, the numerator comprises a gain value indicated by the gain information and the denominator comprises an integer part of a pitch lag indicated by the pitch lag information and, importantly, further comprises a multi-tap filter depending on a fractional part of the pitch lag.


An implementation of this harmonic post filter with a transfer function H(z) is illustrated in FIG. 2. This filter receives the decoder output signal 106 and subjects this decoded output signal to a post-filtering operation to obtain a post-filtered output signal 108. This post-filtered output signal can be output as the processed signal or can be further processed by any procedure for removing any discontinuities introduced by the post-filtering operation which, of course, is signal-dependent, i.e., can vary from frame to frame. This discontinuity removal operation can be any of the well-known discontinuity removal operation such as cross-fading, which means that an earlier frame is faded out and, at the same time, a new frame is faded in and, advantageously, the fading characteristic is so that the fading factors add up to one throughout the cross-fading operation. However, other discontinuity removal such as low-pass filtering or LPC filtering can be applied as well.


The apparatus for processing an audio signal illustrated in FIG. 1 furthermore comprises a multi-tap filter information storage 112 and a filter controller 114. Particularly, the filter controller 114 receives side information 116 from the decoder 100, and this side information can, for example, be the pitch gain information g and the pitch lag information, i.e., information on the integer part Tint of the pitch lag and the fractional part Tfr of the pitch lag. This information is useful for setting the harmonic post-filter from frame to frame and, additionally, for selecting a multi-tap filter information B(z,Tfr). Furthermore, additional information such as the bitrate applied by the decoder or the sampling rate underlying the decoded signal can also be used by the filter control 114 in order to particularly set the scalar values α, β for a certain encoder and/or decoder setting with respect to bitrate and sampling rate.



FIG. 2 illustrates a pole/zero representation of a filter transfer function H(z) in the z domain as known in the art. Naturally, there are numerous other representations of the harmonic post-filter, which are all filter representations, which can be converted to the kind of pole/zero representation in the z domain. Hence, the present invention is applicable for each filter which is describable in any way by such a transfer function as illustrated in the specification.



FIG. 3 illustrates an embodiment of the harmonic post-filter again described a as a transfer function in the pole/zero notation in the z domain.


The filter can be described as follows:







H


(
z
)


=


1
-

αβ






gB


(

z
,
0

)





1
-

β






gB


(

z
,

T
fr


)




z

-

T
int











with g the decoded gain, Tint and Tfr the integer and fractional part of the decoded pitch lag, α and β two scalars that weight the gain, and B(z,Tfr) a low-pass FIR filter whose coefficients depends on the fractional part of the decoded pitch lag.


Note that B(z,0) in the numerator of H(z) is used to compensate for the tilt introduced by B(z,Tfr).


β is used to control the strength of the post-filter. A β equals to 1 produces full effects, suppressing the maximum possible amount of energy between the harmonics. A β equals to 0 disables the post-filter. Generally, a quite low value is used to not suppress too much energy between the harmonics. The value can also depend on the bitrate with a higher value at a lower bitrate, e.g. 0.4 at low bitrate and 0.2 at a high bitrate.


α is used to add a slight tilt to the frequency response of H(z), in order to compensate for the slight loss in energy in the low frequencies. The value of α is generally chosen close to 1, e.g.


0.8.


An example of B(z,TFr) is given in FIG. 6. The order and the coefficients of B(z,Tfr) can also depend on the bitrate and the output sampling rate. A different frequency response can be designed and tuned for each combination of bitrate and output sampling rate.


Particularly, it has been found out that even values for a between 0.6 and lower than 1.0 are useful and that, additionally, values for β between 0.1 and 0.5 have been proved to be useful as well.


Furthermore, the multi-tap filter can have a variable number of taps. It has been found that for certain implementations, four taps are sufficient, where one tap is z+1. However, smaller filters with only two taps or even larger filters with more than four taps are useful for certain implementations.



FIG. 6 illustrates an implementation of filters B(z) for different fractional values of the pitch lag and, particularly, for a pitch lag resolution of ¼. For this implementation, four different filter descriptions for the multi-tap filter in the denominator of the transfer function of the harmonic post-filter are illustrated. However, it has been found that the filter coefficients do not necessarily have to indicate exactly the illustrated values in FIG. 6, but certain variations of +/−0.05 can be useful in other implementations as well.


Particularly, as illustrated in FIG. 1, the tap weights illustrated in FIG. 6 are stored within the memory 112 for the multi-tap filter information. The filter controller 114 receives the fractional part Tfr from line 116 of FIG. 1 and, in response to this value, addresses the memory 112 in order to retrieve, via a retrieval line 200 the specific filter information for the specific fractional part of the pitch lag. This information is then forwarded via an output line 202 to the harmonic post-filter 104 so that the harmonic post-filter is correctly set. A certain implementation of the multi-tap FIR filter is illustrated in FIG. 5. The weight indication w1 to w4 corresponds to the notation in FIG. 6 and the filter controller 114 applies, in response to the actual fractional part of the pitch lag the corresponding weights for a certain audio frame. The other portions such as delay portions 501, 502, 503 and the combiner 505 can be implemented as illustrated. In this context, it is emphasized that the delay value 501 is, in the z notation a negative delay value, since it has been found out that an FIR filter representation having a negative delay value in addition to a positive delay value such as 503 and 504 is particularly useful.


Subsequently, an encoder implementation having certain functional blocks and operating without any pre-filter is illustrated in FIG. 4. The filter portion illustrated in FIG. 4 comprises a pitch estimator 402, a pitch refiner 404, a fractional part estimator 406, a transient detector 408, a gain estimator 410 and a gain quantizer 412. The information provided by the gain quantizer 412, the fractional part estimator 406, the pitch refiner 404 and the decision bit generated by the transient detector 408 are input into an encoded signal former 414. The encoded signal former provides an encoded signal 102, which is then input into the decoder 100 illustrated in FIG. 1. The encoded signal 102 will comprise additional signal information not illustrated in FIG. 4.


Subsequently, the functionality of the pitch estimator 402 is described.


One pitch lag (integer part+fractional part) per frame is estimated (frame size e.g. 20 ms). This is done in 3 steps to reduce complexity and improves estimation accuracy.


A pitch analysis algorithm that produces a smooth pitch evolution contour is used (e.g. Open-loop pitch analysis described in Rec. ITU-T G.718, sec. 6.6). This analysis is generally done on a subframe basis (subframe size e.g. 10 ms), and produces one pitch lag estimate per subframe. Note that these pitch lag estimates do not have any fractional part and are generally estimated on a downsampled signal (sampling rate e.g. 6400 Hz). The signal used can be any audio signal, e.g. a LPC weighted audio signal as described in Rec. ITU-T G.718, sec. 6.5.


The pitch refiner operates as follows:


The final integer part of the pitch lag is estimated on an audio signal x[n] running at the core encoder sampling rate, which is generally higher than the sampling rate of the downsampled signal used in a. (e.g. 12.8 kHz, 16 kHz, 32 kHz . . . ). The signal x[n] can be any audio signal e.g. an LPC weighted audio signal.


The integer part of the pitch lag is then the lag dm that maximizes the autocorrelation function







C


(
d
)


=




n
=
0

N




x


[
n
]




x


[

n
-
d

]








with d around a pitch lag T estimated in step 1.a.

T−δ1≤d≤T+δ2


The fractional part estimator 406 operates as follows:


The fractional part is found by interpolating the autocorrelation function C(d) computed in step 2.b. and selecting the fractional pitch lag which maximizes the interpolated autocorrelation function. The interpolation can be performed using a low-pass FIR filter as described in e.g. Rec. ITU-T G.718, sec. 6.6.7.


The transient detector 408 illustrated in FIG. 4 is configured for generating a decision bit.


If the input audio signal does not contain any harmonic content, then no parameters are encoded in the bitstream. Only 1 bit is sent such that the decoder knows whether he has to decode the post-filter parameters or not. The decision is made based on several parameters:


a. Normalized correlation at the integer pitch lag estimated in step 1.b.







norm
.
corr
.

=





n
=
0

N




x


[
n
]




x


[

n
-

d
m


]










n
=
0

N




x


[
n
]




x


[
n
]
















n
=
0

N




x


[

n
-

d
m


]




x


[

n
-

d
m


]
















The normalized correlation is 1 if the input signal is perfectly predictable by the integer pitch lag, and 0 if it is not predictable at all. A high value (close to 1) would then indicate a harmonic signal. For a more robust decision, the normalized correlation of the past frame can also be used in the decision., e.g.:


If (norm.corr(curr.)*norm.corr.(prev.))>0.25, then the current frame contains some harmonic content (bit=1)


b. Features computed by a transient detector (e.g. Temporal flatness measure, Maximal energy change), to avoid activating the post-filter on a signal containing a transient. e.g. If (tempFlatness>3.5 or maxEnergychange>3.5) then set bit=0 and do not send any parameters


Furthermore, the gain estimator 410 calculates a gain to be input into the gain quantizer 412


The gain is generally estimated on the input audio signal at the core encoder sampling rate, but it can also be any audio signal like the LPC weighted audio signal. This signal is noted y[n] and can be the same or different than x[n].


The prediction yP[n] of y[n] is first found by filtering y[n] with the following filter

P(z)=B(z,Tfr)z−Tint

with Tint the integer part of the pitch lag (estimated in 1.b.) and B(z,Tfr) a low-pass FIR filter whose coefficients depend on the fractional part of the pitch lag Tfr (estimated in 1.c.).


One example of B(z) when the pitch lag resolution is ¼:







T
fr

=



0
4



B


(
z
)



=


0.0000


z

-
2



+

0.2325


z

-
1



+

0.5349


z
0


+

0.2325


z
1











T
fr

=



1
4



B


(
z
)



=


0.0152


z

-
2



+

0.3400


z

-
1



+

0.5094


z
0


+

0.1353


z
1











T
fr

=



2
4



B


(
z
)



=


0.0609


z

-
2



+

0.4391


z

-
1



+

0.4391


z
0


+

0.0609


z
1











T
fr

=



3
4



B


(
z
)



=


0.1353


z

-
2



+

0.5094


z

-
1



+

0.3400


z
0


+

0.0152


z
1








The gain g is then computed as follows:






g
=





n
=
0


N
-
1





y


[
n
]





y
p



[
n
]








n
=
0


N
-
1






y
p



[
n
]





y
p



[
n
]









and limited between 0 and 1.


Finally, the gain is quantized e.g. on 2 bits, using e.g. uniform quantization.


If the gain is quantized to 0, then no parameters are encoded in the bitstream, only the one decision bit (bit=0).


As outlined before, the post-filter is applied on the output audio signal after the transform decoder. It processes the signal on the frame-by-frame basis, with the same frame size as used it the encoder-side such as 20 ms. As illustrated, it is based on a long-term prediction filter H(z) whose parameters are determined from the parameters estimated at the encoder-side and decoded from the bitstream. This information comprises the decision bit, the pitch lag and the gain. If the decision bit is 0, then the pitch lag and the gain are not decoded and are assumed to be 0 not written at all into the bitstream.


As discussed, if the filter parameters are different from one frame to the next frame, a discontinuity can be introduced at the border between the two frames. To avoid discontinuity, a discontinuity remover is applied such as a cross-fader or any other implementation for that purpose.


Furthermore, several different ways to set the harmonic post-filter are illustrated in FIGS. 7a to 8b. The plots illustrate the frequency domain transfer function. The horizontal axis is related to the normalized frequency 1 and the vertical axis is the magnitude of the filter response in dB. It is emphasized that in all illustrations but FIG. 7b, the filter introduces an amplification for low frequencies, i.e., a certain positive dB magnitude value.


Particularly, FIG. 7a illustrates a transfer function, implementing the filter in FIG. 3, with the certain parameter values as indicated above. Furthermore, the a value, i.e., the first scalar value is set to 0. FIG. 7b illustrates a similar situation, but now with an α value equal to 1. The other parameters are identical to FIG. 7a.



FIG. 7c illustrates a further implementation where a is equal to 0.8 which has a slight tilt and a boosting of the lower frequencies. Again, FIG. 7 has the same other parameters as indicated in FIG. 7a. It becomes clear that α equal to 1 removes the tilt and all harmonic frequencies have a gain of 1. The drawback of this setting is a loss of energy at the frequencies between the harmonics. Therefore, a value of a equal to 0.8 as in FIG. 7c is advantageous. This value adds a slight tilt compared to the a equal to 1 situation of FIG. 7b. In order to compensate the loss of energies at the frequencies between the harmonics, this slight tilt may be used.


Furthermore, FIGS. 8a and 8b illustrate filter settings for a value of a equal to 0.8 and different β-values, i.e., a β-value of 0.4 in FIG. 8a and a β-value of 0.2 in FIG. 8b. It becomes clear that a β-value of 0.4 has a stronger post-filtering effect compared to a β-value of 0.2 and, therefore, a β-value of 0.4 is used at lower bitrates in order to remove inter-harmonic noise introduced by such a low bitrate.


On the other hand, β equal 0.2 has a less strong effect for suppressing energy between the harmonics and, therefore, this β-value is advantageous for high bitrates due to the fact that at such higher bitrates, not so much inter-harmonic noise exists.


Although some aspects have been described in the context of an apparatus, it is clear that these aspects also represent a description of the corresponding method, where a block or device corresponds to a method step or a feature of a method step. Analogously, aspects described in the context of a method step also represent a description of a corresponding block or item or feature of a corresponding apparatus. Some or all of the method steps may be executed by (or using) a hardware apparatus, like for example, a microprocessor, a programmable computer or an electronic circuit. In some embodiments, some one or more of the most important method steps may be executed by such an apparatus.


The inventive transmitted or encoded signal can be stored on a digital storage medium or can be transmitted on a transmission medium such as a wireless transmission medium or a wired transmission medium such as the Internet.


Depending on certain implementation requirements, embodiments of the invention can be implemented in hardware or in software. The implementation can be performed using a digital storage medium, for example a floppy disc, a DVD, a Blu-Ray, a CD, a ROM, a PROM, and EPROM, an EEPROM or a FLASH memory, having electronically readable control signals stored thereon, which cooperate (or are capable of cooperating) with a programmable computer system such that the respective method is performed. Therefore, the digital storage medium may be computer readable.


Some embodiments according to the invention comprise a data carrier having electronically readable control signals, which are capable of cooperating with a programmable computer system, such that one of the methods described herein is performed.


Generally, embodiments of the present invention can be implemented as a computer program product with a program code, the program code being operative for performing one of the methods when the computer program product runs on a computer. The program code may, for example, be stored on a machine readable carrier.


Other embodiments comprise the computer program for performing one of the methods described herein, stored on a machine readable carrier.


In other words, an embodiment of the inventive method is, therefore, a computer program having a program code for performing one of the methods described herein, when the computer program runs on a computer.


A further embodiment of the inventive method is, therefore, a data carrier (or a non-transitory storage medium such as a digital storage medium, or a computer-readable medium) comprising, recorded thereon, the computer program for performing one of the methods described herein. The data carrier, the digital storage medium or the recorded medium are typically tangible and/or non-transitory.


A further embodiment of the invention method is, therefore, a data stream or a sequence of signals representing the computer program for performing one of the methods described herein. The data stream or the sequence of signals may, for example, be configured to be transferred via a data communication connection, for example, via the internet.


A further embodiment comprises a processing means, for example, a computer or a programmable logic device, configured to, or adapted to, perform one of the methods described herein.


A further embodiment comprises a computer having installed thereon the computer program for performing one of the methods described herein.


A further embodiment according to the invention comprises an apparatus or a system configured to transfer (for example, electronically or optically) a computer program for performing one of the methods described herein to a receiver. The receiver may, for example, be a computer, a mobile device, a memory device or the like. The apparatus or system may, for example, comprise a file server for transferring the computer program to the receiver.


In some embodiments, a programmable logic device (for example, a field programmable gate array) may be used to perform some or all of the functionalities of the methods described herein. In some embodiments, a field programmable gate array may cooperate with a microprocessor in order to perform one of the methods described herein. Generally, the methods may be performed by any hardware apparatus.


While this invention has been described in terms of several embodiments, there are alterations, permutations, and equivalents which fall within the scope of this invention. It should also be noted that there are many alternative ways of implementing the methods and compositions of the present invention. It is therefore intended that the following appended claims be interpreted as including all such alterations, permutations and equivalents as fall within the true spirit and scope of the present invention.

Claims
  • 1. An apparatus for processing an audio signal having associated therewith a pitch lag information and a gain information, comprising: a domain converter for converting a first domain representation of the audio signal into a second domain representation of the audio signal, the second domain representation being a time domain representation; anda harmonic post-filter for filtering the second domain representation of the audio signal, wherein the harmonic post-filter is based on a along-term prediction filter working in the time-domain,wherein the harmonic post-filter is based on a transfer function comprising a numerator and a denominator, wherein the numerator comprises a gain value indicated by the gain information, and wherein the denominator comprises an integer part of a pitch lag indicated by the pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag indicated by the pitch lag information, orwherein the harmonic post-filter comprises a multi-tap filter being a finite impulse response (FIR) filter and comprising at least three taps, wherein the multi-tap filter in a denominator comprises, for a zero fractional part, four taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.2 and 0.3, the third tap is between 0.5 and 0.6, and the fourth tap is between 0.2 and 0.3, wherein the multi-tap filter comprises, for a first fractional part, four filter taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.3 and 0.4, the third tap is between 0.45 and 0.55, and the fourth tap is between 0.1 and 0.2, wherein the multi-tap filter comprises, for a second fractional part, four filter taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.35 and 0.45, the third tap is between 0.35 and 0.45, and the fourth tap is between 0.0 and 0.1, wherein the multi-tap filter comprises, for a third fractional part, four filter taps, wherein the first tap is between 0.1 and 0.2, the second tap is between 0.45 and 0.55, the third tap is between 0.3 and 0.4, and the fourth tap is between 0.0 and 0.1, wherein the third fractional part is greater than the second fractional part, and wherein the second fractional part is greater than the first fractional part.
  • 2. The apparatus of claim 1, wherein the transfer function of the harmonic post-filter comprises, in the numerator, a further multi-tap FIR filter for a zero fractional part of the pitch lag.
  • 3. The apparatus of claim 1, wherein the denominator comprises a product between the multi-tap filter and the gain value.
  • 4. The apparatus of claim 1, wherein the numerator furthermore comprises a product of a first scalar value and a second scalar value, wherein the denominator comprises the second scalar value and not the first scalar value, wherein the first scalar value and the second scalar value are predetermined and comprise values greater than 0 and wherein the second scalar value is lower than the first scalar value.
  • 5. The apparatus of claim 4, further comprising: a filter controller configured for setting the second scalar value depending on a bitrate, by which the frequency-time converter is operated, wherein the second scalar value is set to a first value, when the bitrate comprises a first value, wherein the second scalar value is set to a second value, when the bitrate comprises a second value, wherein the second value of the bitrate is lower than the first value of the bitrate, and wherein the second value of the second scalar value is greater than the first value of the second scalar value.
  • 6. The apparatus of claim 4, wherein the first scalar value is set between 0.6 and 1.0 and wherein the second scalar value is set between 0.1 and 0.5.
  • 7. Apparatus of claim 6, wherein a cross-fading characteristic of the fading out and the fading in is so that fading factors add up to one throughout the cross-fading operation.
  • 8. The apparatus of claim 1, wherein the harmonic post-filter comprises the transfer function H(z) in a pole-zero representation based on the following equation:
  • 9. The apparatus of claim 1, wherein the harmonic post-filter is configured to comprise a negative spectral tilt for compensating a loss in energy at frequencies between harmonics, orwherein the harmonic post-filter is configured to suppress an amount of energy between harmonics in a frame, wherein the amount of energy suppressed is smaller than 20% of a total energy of the time-domain representation in the frame.
  • 10. The apparatus of claim 1, wherein the domain converter is a frequency-time converter, wherein the first domain is a frequency domain, orwherein the domain converter is an LPC residual-time converter, wherein the first domain is an LPC residual domain.
  • 11. Apparatus of claim 1, wherein the long-term prediction filter, on which the harmonic post-filter is based, is configured to account for an integer part of a pitch lag indicated by the pitch lag information and a fractional part of the pitch lag indicated by the pitch lag information.
  • 12. Apparatus of claim 1, wherein the long-term prediction filter, on which the harmonic post-filter is based, comprises parameters, wherein the parameters are determined from parameters decoded a bitstream comprising the audio signal and the pitch lag information and the gain information.
  • 13. Apparatus of claim 12, wherein the bitstream further comprises a decision bit, and wherein the apparatus is configured to not decode any pitch lag or gain, or to assume the pitch lag and the gain as not written into the bitstream, or to assume the pitch lag and the gain as a zero value, when the decision bit is equal to zero.
  • 14. Apparatus of claim 1, wherein the harmonic post-filter comprises filter parameters derived from the pitch lag information and the gain information, wherein the harmonic post-filter is configured to have different parameters from a frame to a next frame, and wherein the apparatus further comprises a discontinuity remover for reducing a discontinuity at a border between the frame and the next frame.
  • 15. Apparatus of claim 14, wherein the discontinuity remover comprises at least one of a cross-fader, a low-pass filter, or an LPC filter.
  • 16. Apparatus of claim 14, wherein the discontinuity remover is configured to fade out a post filtered audio signal of the frame and, at the same time, to fade in a post filtered audio signal of the next frame.
  • 17. A method of processing an audio signal comprising a method of encoding an audio signal and a method of decoding comprising: converting a frequency representation of the audio signal into a time-domain representation of the audio signal; andfiltering the time-domain representation of the audio signal using a harmonic post-filter, wherein the harmonic post-filter is based on a long-term prediction filter working in the time-domain,wherein the harmonic post-filter is based on a transfer function comprising a numerator and a denominator, wherein the numerator comprises a gain value indicated by the gain information, and wherein the denominator comprises an integer part of a pitch lap indicated by the pitch lag information and a multi-tap filter depending on a fractional part of the pitch lap indicated by the pitch lag information, orwherein the harmonic post-filter comprises a multi-tap filter being a finite impulse response (FIR) filter and comprising at least three taps, wherein the multi-tap filter in a denominator comprises, for a zero fractional part, four taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.2 and 0.3, the third tap is between 0.5 and 0.6, and the fourth tap is between 0.2 and 0.3, wherein the multi-tap filter comprises, for a first fractional part, four filter taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.3 and 0.4, the third tap is between 0.45 and 0.55, and the fourth tap is between 0.1 and 0.2, wherein the multi-tap filter comprises, for a second fractional part, four filter taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.35 and 0.45, the third tap is between 0.35 and 0.45, and the fourth tap is between 0.0 and 0.1, wherein the multi-tap filter comprises, for a third fractional part, four filter taps, wherein the first tap is between 0.1 and 0.2, the second tap is between 0.45 and 0.55, the third tap is between 0.3 and 0.4, and the fourth tap is between 0.0 and 0.1, wherein the third fractional part is greater than the second fractional part, and wherein the second fractional part is greater than the first fractional part.
  • 18. A non-transitory digital storage medium having a computer program stored thereon to perform the method of processing an audio signal having associated therewith a pitch lag information and a gain information, comprising: converting a frequency representation of the audio signal into a time-domain representation of the audio signal; andfiltering the time-domain representation of the audio signal by a harmonic post-filter, wherein the harmonic post-filter is based on a long-term prediction filter working in the time-domain, wherein the harmonic post-filter is based on a transfer function comprising a numerator and a denominator, wherein the numerator comprises a gain value indicated by the gain information, and wherein the denominator comprises an integer part of a pitch lag indicated by the pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag indicated by the pitch lag information, orwherein the harmonic post-filter comprises a multi-tap filter being a finite impulse response (FIR) filter and comprising at least three taps, wherein the multi-tap filter in a denominator comprises, for a zero fractional part, four taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.2 and 0.3, the third tap is between 0.5 and 0.6, and the fourth tap is between 0.2 and 0.3, wherein the multi-tap filter comprises, for a first fractional part, four filter taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.3 and 0.4, the third tap is between 0.45 and 0.55, and the fourth tap is between 0.1 and 0.2, wherein the multi-tap filter comprises, for a second fractional part, four filter taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.35 and 0.45, the third tap is between 0.35 and 0.45, and the fourth tap is between 0.0 and 0.1, wherein the multi-tap filter comprises, for a third fractional part, four filter taps, wherein the first tap is between 0.1 and 0.2, the second tap is between 0.45 and 0.55, the third tap is between 0.3 and 0.4, and the fourth tap is between 0.0 and 0.1, wherein the third fractional part is greater than the second fractional part, and wherein the second fractional part is greater than the first fractional part;when said computer program is run by a computer.
  • 19. A non-transitory digital storage medium having a computer program stored thereon to perform the method of processing an audio signal comprising a method of encoding an audio signal and a method of decoding comprising: converting a frequency representation of the audio signal into a time-domain representation of the audio signal; andfiltering the time-domain representation of the audio signal using a harmonic post-filter, wherein the harmonic post-filter is based on a long-term prediction filter working in the time-domain, wherein the harmonic post-filter is based on a transfer function comprising a numerator and a denominator, wherein the numerator comprises a gain value indicated by the gain information, and wherein the denominator comprises an integer part of a pitch lag indicated by the pitch lag information and a multi-tap filter depending on a fractional part of the pitch lag indicated by the pitch lag information, orwherein the harmonic post-filter comprises a multi-tap filter being a finite impulse response (FIR) filter and comprising at least three taps, wherein the multi-tap filter in a denominator comprises, for a zero fractional part, four taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.2 and 0.3, the third tap is between 0.5 and 0.6, and the fourth tap is between 0.2 and 0.3, wherein the multi-tap filter comprises, for a first fractional part, four filter taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.3 and 0.4, the third tap is between 0.45 and 0.55, and the fourth tap is between 0.1 and 0.2, wherein the multi-tap filter comprises, for a second fractional part, four filter taps, wherein the first tap is between 0.0 and 0.1, the second tap is between 0.35 and 0.45, the third tap is between 0.35 and 0.45, and the fourth tap is between 0.0 and 0.1, wherein the multi-tap filter comprises, for a third fractional part, four filter taps, wherein the first tap is between 0.1 and 0.2, the second tap is between 0.45 and 0.55, the third tap is between 0.3 and 0.4, and the fourth tap is between 0.0 and 0.1, wherein the third fractional part is greater than the second fractional part, and wherein the second fractional part is greater than the first fractional part;when said computer program is run by a computer.
Priority Claims (1)
Number Date Country Kind
14178820 Jul 2014 EP regional
CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation of co-pending U.S. patent application Ser. No. 15/417,231 filed Jan. 27, 2017, which is a continuation of International Application No. PCT/EP2015/066998, filed Jul. 24, 2015, which is incorporated herein by reference in its entirety, and additionally claims priority from European Application No. 14178820.8-1910, filed Jul. 28, 2014, which is incorporated herein by reference in its entirety. The present invention is related to audio processing and, particularly, to audio processing using a harmonic post filter.

US Referenced Citations (20)
Number Name Date Kind
5012517 Wilson et al. Apr 1991 A
5265167 Akamine Nov 1993 A
5293449 Tzeng Mar 1994 A
5359696 Gerson et al. Oct 1994 A
5568688 Andrews et al. Oct 1996 A
5752222 Nishiguchi et al. May 1998 A
5752223 Aoyagi et al. May 1998 A
5774835 Ozawa Jun 1998 A
6058350 Ihara May 2000 A
6058360 Bergstrom May 2000 A
7707034 Sun et al. Apr 2010 B2
8738385 Chen May 2014 B2
10083706 Markovic et al. Sep 2018 B2
10242688 Ravelli Mar 2019 B2
20080319740 Su Dec 2008 A1
20120101824 Chen Apr 2012 A1
20130096912 Resch et al. Apr 2013 A1
20130332151 Fuchs et al. Dec 2013 A1
20150051905 Gao Feb 2015 A1
20160225384 Resch Aug 2016 A1
Foreign Referenced Citations (9)
Number Date Country
698877 Feb 1996 EP
2757560 Jul 2014 EP
H10214100 Aug 1998 JP
2004302257 Oct 2004 JP
2010520505 Jun 2010 JP
2013120225 Jun 2013 JP
2013533983 Aug 2013 JP
2014510301 Apr 2014 JP
2121173 Oct 1998 RU
Non-Patent Literature Citations (12)
Entry
Mustapha et al.; An adaptive post-filtering technique based on a least square approach; 1999 IEEE Workshop on Speech Coding Proceedings. Model, Coders, and Error Criteria (Cat. No. 99EX351); pp. 156-158 (Year: 1999).
Analysis-by-Synthesis Principles; lEEEXplore, Wiley-IEEE EBooks; Jan. 1, 2001, Wiley-IEEE Press 2001, pp. 65-89 (Year: 2001).
Backward-Adaptive Code Excited Linear Prediction; IEEE Xplore; Wiley IEEE EBooks, Jan. 1, 2001, Wiley-IEEE Press 2001 (Year: 2001).
Li et al.; An Improved 1.2 kb/s speech coder based on MELP; IEEE conferences Jan. 1, 2004, Proceedings 7th International Conference on Signal Processing 2004, Proceedings ICSP '04, 2004, pp. 590-593. (Year: 2004).
“EVS Codec Error Concealment of Lost Packets (3GPP TS 26.447 version 12.0.0 Release 12)”, Universal Mobile Telecommunications System (UMTS), LTE, ETSI TS 126 447 V12.0.0, Oct. 2014, Oct. 2014, 1-82.
“Information technology—MPEG audio technologies—Part 3: Unified speech and audio coding”, ISO/IEC JTC 1/SC 29, ISO/IEC FDIS 23003-3:2011(E), ISO/IEC JTC 1/SC 29/WG 11, Sep. 20, 2011, Sep. 20, 2011, 1-291.
“ITU-T G.718”, Series G: Transmission Systems and Media, Digital Systems and Networks Digital Terminal Equipments—Coding of Voice and Audio Signals. Frame Error Robust Narrow-Band and Wideband Embedded Variable Bit-Rate Coding of Speech and Audio From 8-32 kbit/s, Jun. 2008, 00-249.
Chen, Juin-Hwey et al., “Adpative Postfiltering for Quality Enhancement of Coded Speech”, IEEE Transactions on Speech and Audio Processing, vol. 3, No. 1, Jan. 1995, 59-71.
Fuchs, Hendrik , “Improving MPEG Audio Coding by Backward Adaptive Linear Stereo Prediction”, Presented at the 99th Convention AES an Audio Engineering Society Preprint, New York, Oct. 1995, i-28.
Ojanperae, Juha et al., “Long Term Predictor for Transform Domain Perceptual Audio Coding”, Presented at the 107th Convention AES an Audio Engineering Society Preprint, New York, Sep. 1999, i-25.
Song, Jeongook et al., “Harmonic Enhancement in Low Bitrate Audio Coding Using an Efficient Long-Term Predictor”, Hindawi Publishing Corporation EURASIP Journal on Advances in Signal Processing vol. 2010, Article ID 939542, doi: 10.1155/2010/939542, 2010, 1-9.
Yin, Lin et al., “A New Backward Predictor for MPEG Audio Coding”, Presented at the 103rd Convention AES an Audio Engineering Society Preprint, New York, Sep. 1997, i-12.
Related Publications (1)
Number Date Country
20190198034 A1 Jun 2019 US
Continuations (2)
Number Date Country
Parent 15417231 Jan 2017 US
Child 16288018 US
Parent PCT/EP2015/066998 Jul 2015 US
Child 15417231 US