Patent Grant
10930292
The present invention relates to an audio processor and a method for processing an audio signal, a decoder and a method for decoding an audio signal, and an encoder and a method for encoding an audio signal. Furthermore, a calculator and a method for determining phase correction data, an audio signal, and a computer program for performing one of the previously mentioned methods are described. In other words, the present invention shows a phase derivative correction and bandwidth extension (BWE) for perceptual audio codecs or correcting the phase spectrum of bandwidth-extended signals in QMF domain based on perceptual importance.
Perceptual Audio Coding
The perceptual audio coding seen to date follows several common themes, including the use of time/frequency-domain processing, redundancy reduction (entropy coding), and irrelevancy removal through the pronounced exploitation of perceptual effects [1]. Typically, the input signal is analyzed by an analysis filter bank that converts the time domain signal into a spectral (time/frequency) representation. The conversion into spectral coefficients allows for selectively processing signal components depending on their frequency content (e.g. different instruments with their individual overtone structures).
In parallel, the input signal is analyzed with respect to its perceptual properties, i.e. specifically the time- and frequency-dependent masking threshold is computed. The time/frequency dependent masking threshold is delivered to the quantization unit through a target coding threshold in the form of an absolute energy value or a Mask-to-Signal-Ratio (MSR) for each frequency band and coding time frame.
The spectral coefficients delivered by the analysis filter bank are quantized to reduce the data rate needed for representing the signal. This step implies a loss of information and introduces a coding distortion (error, noise) into the signal. In order to minimize the audible impact of this coding noise, the quantizer step sizes are controlled according to the target coding thresholds for each frequency band and frame. Ideally, the coding noise injected into each frequency band is lower than the coding (masking) threshold and thus no degradation in subjective audio is perceptible (removal of irrelevancy). This control of the quantization noise over frequency and time according to psychoacoustic requirements leads to a sophisticated noise shaping effect and is what makes a the coder a perceptual audio coder.
Subsequently, modern audio coders perform entropy coding (e.g. Huffman coding, arithmetic coding) on the quantized spectral data. Entropy coding is a lossless coding step, which further saves on bit rate.
Finally, all coded spectral data and relevant additional parameters (side information, like e.g. the quantizer settings for each frequency band) are packed together into a bitstream, which is the final coded representation intended for file storage or transmission.
Bandwidth Extension
In perceptual audio coding based on filter banks, the main part of the consumed bit rate is usually spent on the quantized spectral coefficients. Thus, at very low bit rates, not enough bits may be available to represent all coefficients in the precision that may be used for achieving perceptually unimpaired reproduction. Thereby, low bit rate requirements effectively set a limit to the audio bandwidth that can be obtained by perceptual audio coding. Bandwidth extension [2] removes this longstanding fundamental limitation. The central idea of bandwidth extension is to complement a band-limited perceptual codec by an additional high-frequency processor that transmits and restores the missing high-frequency content in a compact parametric form. The high frequency content can be generated based on single sideband modulation of the baseband signal, on copy-up techniques like used in Spectral Band Replication (SBR) [3] or on the application of pitch shifting techniques like e.g. the vocoder [4].
Digital Audio Effects
Time-stretching or pitch shifting effects are usually obtained by applying time domain techniques like synchronized overlap-add (SOLA) or frequency domain techniques (vocoder). Also, hybrid systems have been proposed which apply a SOLA processing in subbands. Vocoders and hybrid systems usually suffer from an artifact called phasiness [8] which can be attributed to the loss of vertical phase coherence. Some publications relate improvements on the sound quality of time stretching algorithms by preserving vertical phase coherence where it is important [6][7].
State-of-the-art audio coders [1] usually compromise the perceptual quality of audio signals by neglecting important phase properties of the signal to be coded. A general proposal of correcting phase coherence in perceptual audio coders is addressed in [9].
However, not all kinds of phase coherence errors can be corrected at the same time and not all phase coherence errors are perceptually important. For example, in audio bandwidth extension it is not clear from the state-of-the-art, which phase coherence related errors should be corrected with highest priority and which errors can remain only partly corrected or, with respect to their insignificant perceptual impact, be totally neglected.
Especially due to the application of audio bandwidth extension [2][3][4], the phase coherence over frequency and over time is often impaired. The result is a dull sound that exhibits auditory roughness and may contain additionally perceived tones that disintegrate from auditory objects in the original signal and hence being perceived as an auditory object on its own additionally to the original signal. Moreover, the sound may also appear to come from a far distance, being less “buzzy”, and thus evoking little listener engagement [5]
Therefore, there is a need for an improved approach.
According to an embodiment, an audio processor for processing an audio signal may have: an audio signal phase measure calculator configured for calculating a phase measure of an audio signal for a time frame; a target phase measure determiner for determining a target phase measure for said time frame; a phase corrector configured for correcting phases of the audio signal for the time frame using the calculated phase measure and the target phase measure to achieve a processed audio signal.
According to another embodiment, a decoder for decoding an audio signal may have: an audio processor according to claim 1; a core decoder configured for core decoding an audio signal in a time frame with a reduced number of subbands with respect to the audio signal; a patcher configured for patching a set of subbands of the core decoded audio signal with the reduced number of subbands, wherein the set of subbands forms a first patch, to further subbands in the time frame, adjacent to the reduced number of subbands, to achieve an audio signal with a regular number of subbands; wherein the audio processor is configured for correcting the phases within the subbands of the first patch according to a target function.
According to another embodiment, an encoder for encoding an audio signal may have: a core encoder configured for core encoding the audio signal to achieve a core encoded audio signal having a reduced number of subbands with respect to the audio signal; a fundamental frequency analyzer for analyzing the audio signal or a low-pass filtered version of the audio signal for achieving a fundamental frequency estimate of the audio signal; a parameter extractor configured for extracting parameters of subbands of the audio signal not included in the core encoded audio signal; an output signal former configured for forming an output signal having the core encoded audio signal, the parameters, and the fundamental frequency estimate.
According to another embodiment, a method for processing an audio signal may have the steps of: calculating a phase measure of an audio signal for a time frame with an audio signal phase measure calculator; determining a target phase measure for said time frame with a target phase measure determiner; correcting phases of the audio signal for the time frame with a phase corrector using the calculated phase measure and the target phase measure to achieve a processed audio signal.
According to another embodiment, a method for decoding an audio signal may have the steps of: decoding an audio signal in a time frame with a reduced number of subbands with respect to the audio signal; patching a set of subbands of the decoded audio signal with the reduced number of subbands, wherein the set of subbands forms a first patch, to further subbands in the time frame, adjacent to the reduced number of subbands, to achieve an audio signal with a regular number of subbands; correcting the phases within the subbands of the first patch according to a target function with the audio processor.
According to another embodiment, a method for encoding an audio signal may have the steps of: core encoding the audio signal with a core encoder to achieve a core encoded audio signal having a reduced number of subbands with respect to the audio signal; analyzing the audio signal or a low-pass filtered version of the audio signal with a fundamental frequency analyzer for achieving a fundamental frequency estimate of the audio signal; extracting parameters of subbands of the audio signal not included in the core encoded audio signal with a parameter extractor; forming an output signal having the core encoded audio signal, the parameters, and the fundamental frequency estimate with an output signal former.
According to another embodiment, a non-transitory digital storage medium may have a computer program stored thereon to perform any of the inventive methods.
According to another embodiment, an audio signal may have: a core encoded audio signal having a reduced number of subbands with respect to an original audio signal; a parameter representing subbands of the audio signal not included in the core encoded audio signal; a fundamental frequency estimate of the audio signal or the original audio signal.
The present invention is based on the finding that the phase of an audio signal can be corrected according to a target phase calculated by an audio processor or a decoder. The target phase can be seen as a representation of a phase of an unprocessed audio signal. Therefore, the phase of the processed audio signal is adjusted to better fit the phase of the unprocessed audio signal. Having a, e.g. time frequency representation of the audio signal, the phase of the audio signal may be adjusted for subsequent time frames in a subband, or the phase can be adjusted in a time frame for subsequent frequency subbands. Therefore, a calculator was found to automatically detect and choose the most suitable correction method. The described findings may be implemented in different embodiments or jointly implemented in a decoder and/or encoder.
Embodiments show an audio processor for processing an audio signal comprising an audio signal phase measure calculator configured for calculating a phase measure of an audio signal for a time frame. Furthermore, the audio signal comprises a target phase measure determiner for determining a target phase measure for said time frame and a phase corrector configured for correcting phases of the audio signal for the time frame using the calculated phase measure and the target phase measure to obtain a processed audio signal.
According to further embodiments, the audio signal may comprise a plurality of subband signals for the time frame. The target phase measure determiner is configured for determining a first target phase measure for a first subband signal and a second target phase measure for a second subband signal. Furthermore, the audio signal phase measure calculator determines a first phase measure for the first subband signal and a second phase measure for the second subband signal. The phase corrector is configured for correcting the first phase of the first subband signal using the first phase measure of the audio signal and the first target phase measure and for correcting a second phase of the second subband signal using the second phase measure of the audio signal and the second target phase measure. Therefore, the audio processor may comprise an audio signal synthesizer for synthesizing a corrected audio signal using the corrected first subband signal and the corrected second subband signal.
In accordance with the present invention, the audio processor is configured for correcting the phase of the audio signal in horizontal direction, i.e. a correction over time. Therefore, the audio signal may be subdivided into a set of time frames, wherein the phase of each time frame can be adjusted according to the target phase. The target phase may be a representation of an original audio signal, wherein the audio processor may be part of a decoder for decoding the audio signal which is an encoded representation of the original audio signal. Optionally, the horizontal phase correction can be applied separately for a number of subbands of the audio signal, if the audio signal is available in a time-frequency representation. The correction of the phase of the audio signal may be performed by subtracting a deviation of a phase derivative over time of the target phase and the phase of the audio signal from the phase of the audio signal.
Therefore, since the phase derivative over time is a frequency
with φ being a phase), the described phase correction performs a frequency adjustment for each subband of the audio signal. In other words, the difference of each subband of the audio signal to a target frequency can be reduced to obtain a better quality for the audio signal.
To determine the target phase, the target phase determiner is configured for obtaining a fundamental frequency estimate for a current time frame and for calculating a frequency estimate for each subband of the plurality of subbands of the time frame using the fundamental frequency estimate for the time frame. The frequency estimate can be converted into a phase derivative over time using a total number of subbands and a sampling frequency of the audio signal. In a further embodiment, the audio processor comprises a target phase measure determiner for determining a target phase measure for the audio signal in a time frame, a phase error calculator for calculating a phase error using a phase of the audio signal and the time frame of the target phase measure, and a phase corrector configured for correcting the phase of the audio signal and the time frame using the phase error.
According to further embodiments, the audio signal is available in a time frequency representation, wherein the audio signal comprises a plurality of subbands for the time frame. The target phase measure determiner determines a first target phase measure for a first subband signal and a second target phase measure for a second subband signal. Furthermore, the phase error calculator forms a vector of phase errors, wherein a first element of the vector refers to a first deviation of the phase of the first subband signal and the first target phase measure and wherein a second element of the vector refers to a second deviation of the phase of the second subband signal and the second target phase measure. Additionally, the audio processor of this embodiment comprises an audio signal synthesizer for synthesizing a corrected audio signal using the corrected first subband signal and the corrected second subband signal. This phase correction produces corrected phase values on average.
Additionally or alternatively, the plurality of subbands is grouped into a baseband and a set of frequency patches, wherein the baseband comprises one subband of the audio signal and the set of frequency patches comprises the at least one subband of the baseband at a frequency higher than the frequency of the at least one subband in the baseband.
Further embodiments show the phase error calculator configured for calculating a mean of elements of a vector of phase errors referring to a first patch of the second number of frequency patches to obtain an average phase error. The phase corrector is configured for correcting a phase of the subband signal in the first and subsequent frequency patches of the set of frequency patches of the patch signal using a weighted average phase error, wherein the average phase error is divided according to an index of the frequency patch to obtain a modified patch signal. This phase correction provides good quality at the crossover frequencies, which are the border frequencies between two subsequent frequency patches.
According to a further embodiment, the two previously described embodiments may be combined to obtain a corrected audio signal comprising phase corrected values which are good on average and at the crossover frequencies. Therefore, the audio signal phase derivative calculator is configured for calculating a mean of phase derivatives over frequency for a baseband. The phase corrector calculates a further modified patch signal with an optimized first frequency patch by adding the mean of the phase derivatives over frequency weighted by a current subband index to the phase of the subband signal with the highest subband index in a baseband of the audio signal. Furthermore, the phase corrector may be configured for calculating a weighted mean of the modified patch signal and the further modified patch signal to obtain a combined modified patch signal and for recursively updating, based on the frequency patches, the combined modified patch signal by adding the mean of the phase derivatives over frequency, weighted by the subband index of the current subband, to the phase of the subband signal with the highest subband index in the previous frequency patch of the combined modified patch signal.
To determine the target phase, the target phase measure determiner may comprise a data stream extractor configured for extracting a peak position and a fundamental frequency of peak positions in a current time frame of the audio signal from a data stream. Alternatively, the target phase measure determiner may comprise an audio signal analyzer configured for analyzing the current time frame to calculate a peak position and a fundamental frequency of peak positions in the current time frame. Furthermore, the target phase measure determiner comprises a target spectrum generator for estimating further peak positions in the current time frame using the peak position and the fundamental frequency of peak positions. In detail, the target spectrum generator may comprise a peak detector for generating a pulse train of a time, a signal former to adjust a frequency of the pulse train according to the fundamental frequency of peak positions, a pulse positioner to adjust the phase of the pulse train according to the position, and a spectrum analyzer to generate a phase spectrum of the adjusted pulse train, wherein the phase spectrum of the time domain signal is the target phase measure. The described embodiment of the target phase measure determiner is advantageous for generating a target spectrum for an audio signal having a waveform with peaks.
The embodiments of the second audio processor describe a vertical phase correction. The vertical phase correction adjusts the phase of the audio signal in one time frame over all subbands. The adjustment of the phase of the audio signal, applied independently for each subband, results, after synthesizing the subbands of the audio signal, in a waveform of the audio signal different from the uncorrected audio signal. Therefore, it is e.g. possible to reshape a smeared peak or a transient.
According to a further embodiment, a calculator is shown for determining phase correction data for an audio signal with a variation determiner for determining a variation of the phase of the audio signal in a first and a second variation mode, a variation comparator for comparing a first variation determined using the phase variation mode and a second variation determined using the second variation mode, and a correction data calculator for calculating the phase correction in accordance with the first variation mode or the second variation mode based on a result of the comparing.
A further embodiment shows the variation determiner for determining a standard deviation measure of a phase derivative over time (PDT) for a plurality of time frames of the audio signal as the variation of the phase in the first variation mode or a standard deviation measure of a phase derivative over frequency (PDF) for a plurality of subbands as the variation of the phase in the second variation mode. The variation comparator compares the measure of the phase derivative over time as the first variation mode and the measure of the phase derivative over frequency as the second variation mode for time frames of the audio signal. According to a further embodiment, the variation determiner is configured for determining a variation of the phase of the audio signal in a third variation mode, wherein the third variation mode is a transient detection mode. Therefore, the variation comparator compares the three variation modes and the correction data calculator calculates the phase correction in accordance with the first variation mode, the second variation, or the third variation mode based on a result of the comparing.
The decision rules of the correction data calculator can be described as follows. If a transient is detected, the phase is corrected according to the phase correction for transients to restore the shape of the transient. Otherwise, if the first variation is smaller or equal than the second variation, the phase correction of the first variation mode is applied or, if the second variation is larger than the first variation, the phase correction in accordance with the second variation mode is applied. If the absence of a transient is detected and if both the first and the second variation exceed a threshold value, none of the phase correction modes are applied.
The calculator may be configured for analyzing the audio signal, e.g. in an audio encoding stage, to determine the best phase correction mode and to calculate the relevant parameters for the determined phase correction mode. In a decoding stage, the parameters can be used to obtain a decoded audio signal which has a better quality compared to audio signals decoded using state of the art codecs. It has to be noted that the calculator autonomously detects the right correction mode for each time frame of the audio signal.
Embodiments show a decoder for decoding an audio signal with a first target spectrum generator for generating a target spectrum for a first time frame of a second signal of the audio signal using first correction data and a first phase corrector for correcting a phase of the subband signal in the first time frame of the audio signal determined with a phase correction algorithm, wherein the correction is performed by reducing a difference between a measure of the subband signal in the first time frame of the audio signal and the target spectrum. Additionally, the decoder comprises an audio subband signal calculator for calculating the audio subband signal for the first time frame using a corrected phase for the time frame and for calculating audio subband signal for a second time frame different from the first time frame using the measure of the subband signal in the second time frame or using a corrected phase calculation in accordance with a further phase correction algorithm different from the phase correction algorithm.
According to further embodiments, the decoder comprises a second and a third target spectrum generator equivalent to the first target spectrum generating and a second and a third phase corrector equivalent to the first phase corrector. Therefore, the first phase corrector can perform a horizontal phase correction, the second phase corrector may perform a vertical phase correction, and the third phase corrector can perform phase correction transients. According to a further embodiment the decoder comprises a core decoder configured for decoding the audio signal in a time frame with a reduced number of subbands with respect to the audio signal. Furthermore, the decoder may comprise a patcher for patching a set of subbands of the core decoded audio signal with a reduced number of subbands, wherein the set of subbands forms a first patch, to further subbands in the time frame, adjacent to the reduced number of subbands, to obtain an audio signal with a regular number of subbands. Furthermore, the decoder can comprise a magnitude processor for processing magnitude values of the audio subband signal in the time frame and an audio signal synthesizer for synthesizing audio subband signals or a magnitude of processed audio subband signals to obtain a synthesized decoded audio signal. This embodiment can establish a decoder for bandwidth extension comprising a phase correction of the decoded audio signal.
Accordingly, an encoder for encoding an audio signal comprising a phase determiner for determining a phase of the audio signal, a calculator for determining phase correction data for an audio signal based on the determined phase of the audio signal, a core encoder configured for core encoding the audio signal to obtain a core encoded audio signal having a reduced number of subbands with respect to the audio signal, and a parameter extractor configured for extracting parameters of the audio signal for obtaining a low resolution parameter representation for a second set of subbands not included in the core encoded audio signal, and an audio signal former for forming an output signal comprising the parameters, the core encoded audio signal, and the phase correction data can form an encoder for bandwidth extension.
All of the previously described embodiments may be seen in total or in combination, for example in an encoder and/or a decoder for bandwidth extension with a phase correction of the decoded audio signal. Alternatively, it is also possible to view all of the described embodiments independently without respect to each other.
Embodiments of the present invention will be detailed subsequently referring to the appended drawings, in which:
In the following, embodiments of the invention will be described in further detail. Elements shown in the respective figures having the same or a similar functionality will have associated therewith the same reference signs.
Embodiments of the present invention will be described with regard to a specific signal processing. Therefore,
1 Introduction
Perceptual audio coding has proliferated as mainstream enabling digital technology for all types of applications that provide audio and multimedia to consumers using transmission or storage channels with limited capacity. Modern perceptual audio codecs are expected to deliver satisfactory audio quality at increasingly low bit rates. In turn, one has to put up with certain coding artifacts that are most tolerable by the majority of listeners. Audio Bandwidth Extension (BWE) is a technique to artificially extend the frequency range of an audio coder by spectral translation or transposition of transmitted lowband signal parts into the highband at the price of introducing certain artifacts.
The finding is that some of these artifacts are related to the change of the phase derivative within the artificially extended highband. One of these artifacts is the alteration of phase derivative over frequency (see also “vertical” phase coherence) [8]. Preservation of said phase derivative is perceptually important for tonal signals having a pulse-train like time domain waveform and a rather low fundamental frequency. Artifacts related to a change of the vertical phase derivative correspond to a local dispersion of energy in time and are often found in audio signals which have been processed by BWE techniques. Another artifact is the alteration of the phase derivative over time (see also “horizontal” phase coherence) which is perceptually important for overtone-rich tonal signals of any fundamental frequency. Artifacts related to an alteration of the horizontal phase derivative correspond to a local frequency offset in pitch and are often found in audio signals which have been processed by BWE techniques.
The present invention presents means for readjusting either the vertical or horizontal phase derivative of such signals when this property has been compromised by application of so-called audio bandwidth extension (BWE). Further means are provided to decide if a restoration of the phase derivative is perceptually beneficial and whether adjusting the vertical or horizontal phase derivative is perceptually advantageous.
Bandwidth-extension methods, such as spectral band replication (SBR) [9], are often used in low-bit-rate codecs. They allow transmitting only a relatively narrow low-frequency region alongside with parametric information about the higher bands. Since the bit rate of the parametric information is small, significant improvement in the coding efficiency can be obtained.
Typically the signal for the higher bands is obtained by simply copying it from the transmitted low-frequency region. The processing is usually performed in the complex-modulated quadrature-mirror-filter-bank (QMF) [10] domain, which is assumed also in the following. The copied-up signal is processed by multiplying the magnitude spectrum of it with suitable gains based on the transmitted parameters. The aim is to obtain a similar magnitude spectrum as that of the original signal. On the contrary, the phase spectrum of the copied-up signal is typically not processed at all, but, instead, the copied-up phase spectrum is directly used.
The perceptual consequences of using directly the copied-up phase spectrum is investigated in the following. Based on the observed effects, two metrics for detecting the perceptually most significant effects are suggested. Moreover, methods how to correct the phase spectrum based on them are suggested. Finally, strategies for minimizing the amount of transmitted parameter values for performing the correction are suggested.
The present invention is related to the finding that preservation or restoration of the phase derivative is able to remedy prominent artifacts induced by audio bandwidth extension (BWE) techniques. For instance, typical signals, where the preservation of the phase derivative is important, are tones with rich harmonic overtone content, such as voiced speech, brass instruments or bowed strings.
The present invention further provides means to decide if—for a given signal frame—a restoration of the phase derivative is perceptually beneficial and whether adjusting the vertical or horizontal phase derivative is perceptually advantageous.
The invention teaches an apparatus and a method for phase derivative correction in audio codecs using BWE techniques with the following aspects:
A time-domain signal x(m), where m is discrete time, can be presented in the time-frequency domain, e.g. using a complex-modulated Quadrature Mirror Filter bank (QMF). The resulting signal is X(k,n), where k is the frequency band index and n the temporal frame index. The QMF of 64 bands and the sampling frequency fs of 48 kHz are assumed for visualizations and embodiments. Thus, the bandwidth fBW of each frequency band is 375 Hz and the temporal hop size thop (17 in
The resulting signal is X(k,n), where k is the frequency band index and n the temporal frame index. X(k,n) is a complex signal. Thus, it can also be presented using the magnitude Xmag(k,n) and the phase components Xpha(k,n) with j being the complex number
X(k,n)=Xmag(k,n)ejX
The audio signals are presented mostly using Xmag(k,n) and Xpha(k,n) (see
3 Audio Data
The audio data used to show an effect of a described audio processing are named ‘trombone’ for an audio signal of a trombone, ‘violin’ for an audio signal of a violin, and ‘violin+clap’ for the violin signal with a hand clap added in the middle.
4 Basic Operation of SBR
Assuming a signal X(k,n), the bandwidth-extended (BWE) signal Z(k,n) is obtained from the input signal X(k,n) by copying up certain parts of the transmitted low-frequency frequency band. An SBR algorithm starts by selecting a frequency region to be transmitted. In this example, the bands from 1 to 7 are selected:
∀1≤k≤7:Xtrans(k,n)=X(k,n). (2)
The amount of frequency bands to be transmitted depends on the desired bit rate. The figures and the equations are produced using 7 bands, and from 5 to 11 bands are used for the corresponding audio data. Thus, the cross-over frequencies between the transmitted frequency region and the higher bands are from 1875 to 4125 Hz, respectively. The frequency bands above this region are not transmitted at all, but instead, parametric metadata is created for describing them. Xtrans (k,n) is coded and transmitted. For the sake of simplicity, it is assumed that the coding does not modify the signal in any way, even though it has to be seen that the further processing is not limited to the assumed case.
In the receiving end, the transmitted frequency region is directly used for the corresponding frequencies.
For the higher bands, the signal may be created somehow using the transmitted signal. One approach is simply to copy the transmitted signal to higher frequencies. A slightly modified version is used here. First, a baseband signal is selected. It could be the whole transmitted signal, but in this embodiment the first frequency band is omitted. The reason for this is that the phase spectrum was noticed to be irregular for the first band in many cases. Thus, the baseband to be copied up is defined as
∀1≤k≤6:Xbase(k,n)=Xtrans(k+1,n). (3)
Other bandwidths can also be used for the transmitted and the baseband signals. Using the baseband signal, raw signals for the higher frequencies are created
Yraw(k,n,i)=Xbase(k,n), (4)
where Yraw(k,n,i) is the complex QMF signal for the frequency patch i. The raw frequency-patch signals are manipulated according to the transmitted metadata by multiplying them with gains g(k,n,i)
Y(k,n,i)=Yraw(k,n,i)g(k,n,i). (5)
It should be noted that the gains are real valued, and thus, only the magnitude spectrum is affected and thereby adapted to a desired target value. Known approaches show how the gains are obtained. The target phase remains non-corrected in said known approaches.
The final signal to be reproduced is obtained by concatenating the transmitted and the patch signals for seamlessly extending the bandwidth to obtain a BWE signal of the desired bandwidth. In this embodiment, i=7 is assumed.
Z(k,n)=Xtrans(k,n),
Z(k+6i+1,n)=Y(k,n,i). (6)
It is assumed that the parametric representation of the higher bands is perfect, i.e., the magnitude spectrum of the reconstructed signal is identical to that of the original signal
Zmag(k,n)=Xmag(k,n). (7)
However, it should be noted that the phase spectrum is not corrected in any way by the algorithm, so it is not correct even if the algorithm worked perfectly. Therefore, embodiments show how to additionally adapt and correct the phase spectrum of Z(k,n) to a target value such that an improvement of the perceptual quality is obtained. In embodiments, the correction can be performed using three different processing modes, “horizontal”, “vertical” and “transient”. These modes are separately discussed in the following.
Zmag(k,n) and Zpha(k,n) are depicted in
5 Meaning of the Phase Spectrum in the QMF Domain
Often it is thought that the index of the frequency band defines the frequency of a single tonal component, the magnitude defines the level of it, and the phase defines the ‘timing’ of it. However, the bandwidth of a QMF band is relatively large, and the data is oversampled. Thus, the interaction between the time-frequency tiles (i.e., QMF bins) actually defines all of these properties.
A time-domain presentation of a single QMF bin with three different phase values, i.e., Xmag(3,1)=1 and Xpha(3,1)=0, π/2, or π is depicted in
Considering a case where only one frequency band is non-zero for all temporal frames, i.e.,
∀n∃:Xmag(3,n)=1. (8)
By changing the phase between the temporal frames with a fixed value α, i.e.,
Xpha(k,n)=Xpha(k,n−1)+α, (9)
a sinusoid is created. The resulting signal (i.e., the time-domain signal after inverse QMF transform) is presented in
Correspondingly, if the phase is selected randomly, the result is narrow-band noise (see
Considering a case where only one temporal frame is non-zero for all frequency bands, i.e.,
∀k∃:Xmag(k,3)=1. (10)
By changing the phase between the frequency bands with a fixed value α, i.e.,
Xpha(k,n)=Xpha(k−1,n)+α, (11)
a transient is created. The resulting signal (i.e., the time-domain signal after inverse QMF transform) is presented in
Correspondingly, if the phase is selected randomly, the result is a short noise burst (see
6 Measures for Describing Perceptually Relevant Properties of the Phase Spectrum
As discussed in Section 4, the phase spectrum in itself looks quite messy, and it is difficult to see directly what its effect on perception is. Section 5 presented two effects that can be caused by manipulating the phase spectrum in the QMF domain: (a) constant phase change over time produces a sinusoid and the amount of phase change controls the frequency of the sinusoid, and (b) constant phase change over frequency produces a transient and the amount of phase change controls the temporal position of the transient.
The frequency and the temporal position of a partial are obviously significant to human perception, so detecting these properties is potentially useful. They can be estimated by computing the phase derivative over time (PDT)
Xpdt(k,n)=Xpha(k,n+1)−Xpha(k,n) (12)
and by computing the phase derivative over frequency (PDF)
Xpdf(k,n)=Xpha(k+1,n)−Xpha(k,n). (13)
Xpdt(k,n) is related to the frequency and Xpdf(k,n) to the temporal position of a partial. Due to the properties of the QMF analysis (how the phases of the modulators of the adjacent temporal frames match at the position of a transient), π is added to the even temporal frames of Xpdf(k,n) in the figures for visualization purposes in order to produce smooth curves.
Next it is inspected how these measures look like for our example signals.
For the trombone, Xpdt is relatively noisy. On the contrary, the Xpdf appears to have about the same value at all frequencies. In practice, this means that all the harmonic components are aligned in time producing a transient-like signal. The temporal locations of the transients are determined by the Xpdf values.
The same derivatives can also be computed for the SBR-processed signals Z(k,n) (see
Correspondingly, PDF of the frequency patches is otherwise identical to that of the baseband, but at the cross-over frequencies the PDF is, in practice, random. At the cross-over, the PDF is actually computed between the last and the first phase value of the frequency patch, i.e.,
Zpdt(7,n)=Zpha(8,n)−Zpha(7,n)=Ypha(1,n,i)−Ypha(6,n,i) (14)
These values depend on the actual PDF and the cross-over frequency, and they do not match with the values of the original signal.
For the trombone, the PDF values of the copied-up signal are correct apart from the cross-over frequencies. Thus, the temporal locations of the most of the harmonics are in the correct places, but the harmonics at the cross-over frequencies are practically at random locations. The perceptual effect of this is discussed in Section 7.
7 Human Perception of Phase Errors
Sounds can roughly be divided into two categories: harmonic and noise-like signals. The noise-like signals have, already by definition, noisy phase properties. Thus, the phase errors caused by SBR are assumed not to be perceptually significant with them. Instead, it is concentrated on harmonic signals. Most of the musical instruments, and also speech, produce harmonic structure to the signal, i.e., the tone contains strong sinusoidal components spaced in frequency by the fundamental frequency.
Human hearing is often assumed to behave as if it contained a bank of overlapping band-pass filters, referred to as the auditory filters. Thus, the hearing can be assumed to handle complex sounds so that the partial sounds inside the auditory filter are analyzed as one entity. The width of these filters can be approximated to follow the equivalent rectangular bandwidth (ERB) [11], which can be determined according to
ERB=24.7(4.37fc+1), (15)
where fc is the center frequency of the band (in kHz). As discussed in Section 4, the cross-over frequency between the baseband and the SBR patches is around 3 kHz. At these frequencies the ERB is about 350 Hz. The bandwidth of a QMF frequency band is actually relatively close to this, 375 Hz. Hence, the bandwidth of the QMF frequency bands can be assumed to follow ERB at the frequencies of interest.
Two properties of a sound that can go wrong due to erroneous phase spectrum were observed in Section 6: the frequency and the timing of a partial component. Concentrate on the frequency, the question is, can human hearing perceive the frequencies of individual harmonics? If it can, then the frequency offset caused by SBR should be corrected, and if not, then correction is not required.
The concept of resolved and unresolved harmonics [12] can be used to clarify this topic. If there is only one harmonic inside the ERB, the harmonic is called resolved. It is typically assumed that the human hearing processes resolved harmonics individually and, thus, is sensitive to the frequency of them. In practice, changing the frequency of resolved harmonics is perceived to cause inharmonicity.
Correspondingly, if there are multiple harmonics inside the ERB, the harmonics are called unresolved. The human hearing is assumed not to process these harmonics individually, but instead, their joint effect is seen by the auditory system. The result is a periodic signal and the length of the period is determined by the spacing of the harmonics. The pitch perception is related to the length of the period, so human hearing is assumed to be sensitive to it. Nevertheless, if all harmonics inside the frequency patch in SBR are shifted by the same amount, the spacing between the harmonics, and thus the perceived pitch, remains the same. Hence, in the case of unresolved harmonics, human hearing does not perceive frequency offsets as inharmonicity.
Timing-related errors caused by SBR are considered next. By timing the temporal position, or the phase, of a harmonic component is meant. This should not be confused with the phase of a QMF bin. The perception of timing-related errors was studied in detail in [13]. It was observed that for the most of the signals human hearing is not sensitive to the timing, or the phase, of the harmonic components. However, there are certain signals with which the human hearing is very sensitive to the timing of the partials. The signals include, for example, trombone and trumpet sounds and speech. With these signals, a certain phase angle takes place at the same time instant with all harmonics. Neural firing rate of different auditory bands were simulated in [13]. It was found out that with these phase-sensitive signals the produced neural firing rate is peaky at all auditory bands and that the peaks are aligned in time. Changing the phase of even a single harmonic can change the peakedness of the neural firing rate with these signals. According to the results of the formal listening test, human hearing is sensitive to this [13]. The produced effects are the perception of an added sinusoidal component or a narrowband noise at the frequencies where the phase was modified.
In addition, it was found out that the sensitivity to the timing-related effects depends on the fundamental frequency of the harmonic tone [13]. The lower the fundamental frequency, the larger are the perceived effects. If the fundamental frequency is above about 800 Hz, the auditory system is not sensitive at all to the timing-related effects.
Thus, if the fundamental frequency is low and if the phase of the harmonics is aligned over frequency (which means that the temporal positions of the harmonics are aligned), changes in the timing, or in other words the phase, of the harmonics can be perceived by the human hearing. If the fundamental frequency is high and/or the phase of the harmonics is not aligned over frequency, the human hearing is not sensitive to changes in the timing of the harmonics.
8 Correction Methods
In Section 7, it was noted that humans are sensitive to errors in the frequencies of resolved harmonics. In addition, humans are sensitive to errors in the temporal positions of the harmonics if the fundamental frequency is low and if the harmonics are aligned over frequency. SBR can cause both of these errors, as discussed in Section 6, so the perceived quality can be improved by correcting them. Methods for doing so are suggested in this section.
8.1 Correcting Frequency Errors—Horizontal Phase Derivative Correction
As discussed in Section 7, humans can perceive an error in the frequency of a harmonic mostly when there is only one harmonic inside one ERB. Furthermore, the bandwidth of a QMF frequency band can be used to estimate ERB at the first cross over. Hence, the frequency has to be corrected only when there is one harmonic inside one frequency band. This is very convenient, since Section 5 showed that, if there is one harmonic per band, the produced PDT values are stable, or slowly changing over time, and can potentially be corrected using low bit rate.
Embodiments show the phase corrector 70 being configured for correcting subband signals 95 of different subbands of the audio signal 55 within the time frame 75, so that frequencies of corrected subband signals 95 have frequency values being harmonically allocated to a fundamental frequency of the audio signal 55. The fundamental frequency is the lowest frequency occurring in the audio signal 55, or in other words, the first harmonics of the audio signal 55.
Furthermore, the phase corrector 70 is configured for smoothing the deviation 105 for each subband 95 of the plurality of subbands over a previous time frame, the current time frame, and a future time frame 75a to 75c and is configured for reducing rapid changes of the deviation 105 within a subband 95. According to further embodiments, the smoothing is a weighted mean, wherein the phase corrector 70 is configured for calculating the weighted mean over the previous, the current and the future time frames 75a to 75c, weighted by a magnitude of the audio signal 55 in the previous, the current and the future time frame 75a to 75c.
Embodiments show the previously described processing steps vector based. Therefore, the phase corrector 70 is configured for forming a vector of deviations 105, wherein a first element of the vector refers to a first deviation 105a for the first subband 95a of the plurality of subbands and a second element of the vector refers to a second deviation 105b for the second subband 95b of the plurality of subbands from a previous time frame 75a to a current time frame 75b. Furthermore, the phase corrector 70 can apply the vector of deviations 105 to the phases 45 of the audio signal 55, wherein the first element of the vector is applied to a phase 45a of the audio signal 55 in a first subband 95a of a plurality of subbands of the audio signal 55 and the second element of the vector is applied to a phase 45b of the audio signal 55 in a second subband 95b of the plurality of subbands of the audio signal 55.
From another point of view, it can be stated that the whole processing in the audio processor 50 is vector-based, wherein each vector represents a time frame 75, wherein each subband 95 of the plurality of subband comprises an element of the vector. Further embodiments focus on the target phase measure determiner which is configured for obtaining a fundamental frequency estimate 85b for a current time frame 75b, wherein the target phase measure determiner 65 is configured for calculating a frequency estimate 85 for each subband of the plurality of subbands for the time frame 75 using the fundamental frequency estimate 85 for the time frame 75. Furthermore, the target phase measure determiner 65 may convert the frequency estimates 85 for each subband 95 of the plurality of subbands into a phase derivative over time using a total number of subbands 95 and a sampling frequency of the audio signal 55. For clarification it has to be noted that the output 85 of the target phase measure determiner 65 may be either the frequency estimate or the phase derivative over time, depending on the embodiment. Therefore, in one embodiment the frequency estimate already comprises the right format for further processing in the phase corrector 70, wherein in another embodiment the frequency estimate has to be converted into a suitable format, which may be a phase derivative over time.
Accordingly, the target phase measure determiner 65 may be seen as vector based as well. Therefore, the target phase measure determiner 65 can form a vector of frequency estimates 85 for each subband 95 of the plurality of subbands, wherein the first element of the vector refers to a frequency estimate 85a for a first subband 95a and a second element of the vector refers to a frequency estimate 85b for a second subband 95b. Additionally, the target phase measure determiner 65 can calculate the frequency estimate 85 using multiples of the fundamental frequency, wherein the frequency estimate 85 of the current subband 95 is that multiple of the fundamental frequency which is closest to the center of the subband 95, or wherein the frequency estimate 85 of the current subband is a border frequency of the current subband 95 if none of the multiples of the fundamental frequency are within the current subband 95.
In other words, the suggested algorithm for correcting the errors in the frequencies of the harmonics using the audio processor 50 functions as follows. First, the PDT is computed and the SBR processed signal Zpdt. Zpdt(k,n)=Zpha(k,n+1)−Zpha(k,n). The difference between it and a target PDT for the horizontal correction is computed next:
Dpdt(k,n)=Zpdt(k,n)+Zthpdt(k,n). (16a)
At this point the target PDT can be assumed to be equal to the PDT of the input of the input signal
Zthpdt(k,n)=Xpdt(k,n). (16b)
Later it will be presented how the target PDT can be obtained with a low bit rate.
This value (i.e. the error value 105) is smoothened over time using a Hann window W(l). Suitable length is, for example, 41 samples in the QMF domain (corresponding to an interval of 55 ms). The smoothing is weighted by the magnitude of the corresponding time-frequency tiles
Dsmpdt(k,n)=circmean{Dpdt(k,n+l),W(l)Zmag(k,n+l)},−20≤l≤20, (17)
where circmean {a, b} denotes computing the circular mean for angular values a weighted by values b. The smoothened error in the PDT Dsmpdt(k,n) is depicted in
Next, a modulator matrix is created for modifying the phase spectrum in order to obtain the desired PDT
Qpha(k,n+1)=Qpha(k,n)−Dsmpdt(k,n). (18)
The phase spectrum is processed using this matrix
Zchpha(k,n)=Zpha(k,n)+Qpha(k,n). (19)
Using Xpdt(k,n) as a target PDT, it is likely to transmit the PDT-error values Dsmpdt(k,n) for each time-frequency tile. A further approach calculating the target PDT such that the bandwidth for transmission is reduced is shown in section 9.
In further embodiments, the audio processor 50 may be part of a decoder 110. Therefore, the decoder 110 for decoding an audio signal 55 may comprise the audio processor 50, a core decoder 115, and a patcher 120. The core decoder 115 is configured for core decoding an audio signal 25 in a time frame 75 with a reduced number of subbands with respect to the audio signal 55. The patcher patches a set of subbands 95 of the core decoded audio signal 25 with a reduced number of subbands, wherein the set of subbands forms a first patch 30a, to further subbands in the time frame 75, adjacent to the reduced number of subbands, to obtain an audio signal 55 with a regular number of subbands. Additionally, the audio processor 50 is configured for correcting the phases 45 within the subbands of the first patch 30a according to a target function 85. The audio processor 50 and the audio signal 55 have been described with respect to
According to further embodiments, the patcher 120 is configured for patching a set of subbands 95 of the audio signal 25, wherein the set of subbands forms a second patch, to further subbands of the time frame, adjacent to the first patch and wherein the audio processor 50 is configured for correcting the phase 45 within the subbands of the second patch. Alternatively, the patcher 120 is configured for patching the corrected first patch to further subbands of the time frame, adjacent to the first patch.
In other words, in the first option the patcher builds an audio signal with a regular number of subbands from the transmitted part of the audio signal and thereafter the phases of each patch of the audio signal are corrected. The second option first corrects the phases of the first patch with respect to the transmitted part of the audio signal and thereafter builds the audio signal with the regular number of subbands with the already corrected first patch.
Further embodiments show the decoder 110 comprising a data stream extractor 130 configured for extracting a fundamental frequency 114 of the current time frame 75 of the audio signal 55 from a data stream 135, wherein the data stream further comprises the encoded audio signal 145 with a reduced number of subbands. Alternatively, the decoder may comprise a fundamental frequency analyzer 150 configured for analyzing the core decoded audio signal 25 in order to calculate the fundamental frequency 140. In other words, options for deriving the fundamental frequency 140 are for example an analysis of the audio signal in the decoder or in the encoder, wherein in the latter case the fundamental frequency may be more accurate at the cost of a higher data rate, since the value has to be transmitted from the encoder to the decoder.
In an alternative embodiment an intelligent gap filling encoder may be used for encoding the audio signal 55. Therefore, the core encoder encodes a full bandwidth audio signal, wherein at least one subband of the audio signal is left out. Therefore, the parameter extractor 165 extracts parameters for reconstructing the subbands being left out from the encoding process of the core encoder 160.
The described methods 2300, 2400 and 2500 may be implemented in a program code of a computer program for performing the methods when the computer program runs on a computer.
8.2 Correcting Temporal Errors—Vertical Phase Derivative Correction
As discussed previously, humans can perceive an error in the temporal position of a harmonic if the harmonics are synced over frequency and if the fundamental frequency is low. In Section 5 it was shown that the harmonics are synced if the phase derivative over frequency is constant in the QMF domain. Therefore, it is advantageous to have at least one harmonic in each frequency band. Otherwise the ‘empty’ frequency bands would have random phases and would disturb this measure. Luckily, humans are sensitive to the temporal location of the harmonics only when the fundamental frequency is low (see Section 7). Thus, the phase derivate over frequency can be used as a measure for determining perceptually significant effects due to temporal movements of the harmonics.
Regarding further embodiments, the plurality of subbands 95 is grouped into a baseband 30 and a set of frequency patches 40, the baseband 30 comprising one subband 95 of the audio signal 55 and the set of frequency patches 40 comprises the at least one subband 95 of the baseband 30 at a frequency higher than the frequency of the at least one subband in the baseband. It has to be noted that the patching of the audio signal has already been described with respect to
A further embodiment is depicted at the bottom of
In a further embodiment, the audio signal phase derivative calculator 210 is configured for calculating a mean of phase derivatives over frequency 215 for a plurality of subband signals comprising higher frequencies than the baseband signal 30 to detect transients in the subband signal 95. It has to be noted that the transient correction is similar to the vertical phase correction of the audio processor 50′ with the difference that the frequencies in the baseband 30 do not reflect the higher frequencies of a transient. Therefore, these frequencies have to be taken into consideration for the phase correction of a transient.
After the initialization step, the phase correct 70′ is configured for recursively updating, based on the frequency patches 40, the further modified patch signal 40″ by adding the mean of the phase derivatives over frequency 215, weighted by the subband index of the current subband 95, to the phase of the subband signal with the highest subband index in the previous frequency patch. The advantageous embodiment is a combination of the previously described embodiments, where the phase corrector 70′ calculates a weighted mean of the modified patch signal 40′ and the further modified patch signal 40″ to obtain a combined modified patch signal 40′″. Therefore, the phase corrector 70′ recursively updates, based on the frequency patches 40, a combined modified patch signal 40′″ by adding the mean of the phase derivatives over frequency 215, weighted by the subband index of the current subband 95 to the phase of the subband signal with the highest subband index in the previous frequency patch of the combined modified patch signal 40′″. To obtain the combined modified patches 40a′″, 40b′″, etc., the switch 220b is shifted to the next position after each recursion, starting at the combined modified 48″′ for the initialization step, switching to the combined modified patch 40b′″ after the first recursion and so on.
Furthermore, the phase corrector 70′ may calculate a weighted mean of a patch signal 40′ and the modified patch signal 40″ using a circular mean of the patch signal 40′ in the current frequency patch weighted with a first specific weighting function and the modified patch signal 40″ in the current frequency patch weighted with a second specific weighting function.
In order to provide an interoperability between the audio processor 50 and the audio processor 50′, the phase corrector 70′ may form a vector of phase deviations, wherein the phase deviations are calculated using a combined modified patch signal 40′″ and the audio signal 55.
The second correction mode is therefore applied on the combined modified patch signal 40′″ to obtain the modified patch signal 40″ for the second time frame 75b. Additionally, the first correction mode is applied on the patches of the audio signal 55 in the second time frame 75b to obtain the patch signal 40′. Again, a combination of the patch signal 40′ and the modified patch signal 40″ results in the combined modified patch signal 40″. The processing scheme described for the second time frame is applied to the third time frame 75c and any further time frame of the audio signal 55 accordingly.
According to a further embodiment, the patcher 120 is configured for patching the set of subbands 95 of the audio signal 25, wherein the set of subbands forms a further patch, to further subbands of the time frame, adjacent to the patch, and wherein the audio processor 50′ is configured for correcting the phases within the subbands of the further patch. Alternatively, the patcher 120 is configured for patching the corrected patch to further subbands of the time frame adjacent to the patch.
A further embodiment is related to a decoder for decoding an audio signal comprising a transient, wherein the audio processor 50′ is configured to correct the phase of the transient. The transient handling is described in other word in section 8.4. Therefore, the decoder 110 comprises a further audio processor 50′ for receiving a further phase derivative of a frequency and to correct transients in the audio signal 32 using the received phase derivative or frequency. Furthermore, it has to be noted that the decoder 110′ of
In other words, the suggested algorithm for correcting the errors in the temporal positions of the harmonics functions as follows. First, a difference between the phase spectra of the target signal and the SBR-processed signal (Ztvpha(k,n) and Zpha) is computed
Dpha(k,n)=Zpha(k,n)−Ztvpha(k,n), (20a)
which is depicted in
Ztvpha(k,n)=Xpha(k,n) (20b)
Later it will be presented how the target phase spectrum can be obtained with a low bit rate.
The vertical phase derivative correction is performed using two methods, and the final corrected phase spectrum is obtained as a mix of them.
First, it can be seen that the error is relatively constant inside the frequency patch, and the error jumps to a new value when entering a new frequency patch. This makes sense, since the phase is changing with a constant value over frequency at all frequencies in the original signal. The error is formed at the cross-over and the error remains constant inside the patch. Thus, a single value is enough for correcting the phase error for the whole frequency patch. Furthermore, the phase error of the higher frequency patches can be corrected using this same error value after multiplication with the index number of the frequency patch.
Therefore, circular mean of the phase error is computed for the first frequency patch
Davgpha(n)=circmean{Dpha(k,n)},8≤k≤13. (21)
The phase spectrum can be corrected using it
Ycv1pha(k,n,i)=Ypha(k,n,i)−i·Davgpha(n). (22)
This raw correction produces an accurate result if the target PDF, e.g. the phase derivative over frequency Xpdf(k,n), is exactly constant at all frequencies. However, as can be seen in
The other correction method begins by computing a mean of the PDF in the baseband
Xavgpdf(n)=circmean{Xbasepdf(k,n)}. (23)
The phase spectrum can be corrected using this measure by assuming that the phase is changing with this average value, i.e.,
Ycv2pha(k,n,1)=Xbasepha(6,n)+k·Xavgpdf(n),
Ycv2pha(k,n,i)=Ycvpha(6,n,i−1)+k·Xavgpdf(n), (24)
wherein Ycvpha is the combined patch signal of the two correction methods.
This correction provides good quality at the cross-overs, but can cause a drift in the PDF towards higher frequencies. In order to avoid this, the two correction methods are combined by computing a weighted circular mean of them
Ycvpha(k,n,i)=circmean{Ycv12pha(k,n,i,c),Wfc(k,c)}, (25)
where c denotes the correction method (Ycv1pha or Ycv2pha and Wfc(k,c) is the weighting function
Wfc(k,1)=[0.2,0.45,0.7,1,1,1],
Wfc(k,2)=[0.8,0.55,0.3,0,0,0]. (26a)
The resulting phase spectrum Ycvpha(k,n,i) suffers neither from discontinuities nor drifting. The error compared to the original spectrum and the PDF of the corrected phase spectrum are depicted in
The corrected phase spectrum Zcvpha(k,n) is obtained by concatenating the corrected frequency patches Ycvpha(k,n,i). To be compatible with the horizontal-correction mode, the vertical phase correction can be presented also using a modulator matrix (see Eq. 18)
Qpha(k,n)=Zcvpha(k,n)−Zpha(k,n). (26b)
8.3 Switching Between Different Phase-Correction Methods
Sections 8.1 and 8.2 showed that SBR-induced phase errors can be corrected by applying PDT correction to the violin and PDF correction to the trombone. However, it was not considered how to know which one of the corrections should be applied to an unknown signal, or if any of them should be applied. This section proposes a method for automatically selecting the correction direction. The correction direction (horizontal/vertical) is decided based on the variation of the phase derivatives of the input signal.
Therefore, in
Furthermore, the variation determiner 275 may be configured for determining a standard deviation measure of a phase derivative over time (PDT) for a plurality of time frames of the audio signal 55 as the variation 290a of the phase in the first variation mode and for determining a standard deviation measure of a phase derivative over frequency (PDF) for a plurality of subbands of the audio signal 55 as the variation 290b of the phase in the second variation mode. Therefore, the variation comparator 280 compares the measure of the phase derivative over time as the first variation 290a and the measure of the phase derivative over frequency as a second variation 290b for time frames of the audio signal.
Embodiments show the variation determiner 275 for determining a circular standard deviation of a phase derivative over time of a current and a plurality of previous frames of the audio signal 55 as the standard deviation measure and for determining a circular standard deviation of a phase derivative over time of a current and a plurality of future frames of the audio signal 55 for a current time frame as the standard deviation measure. Furthermore, the variation determiner 275 calculates, when determining the first variation 290a, a minimum of both circular standard deviations. In a further embodiment, the variation determiner 275 calculates the variation 290a in the first variation mode as a combination of a standard deviation measure for a plurality of subbands 95 in a time frame 75 to form an averaged standard deviation measure of a frequency. The variation comparator 280 is configured for performing the combination of the standard deviation measures by calculating an energy-weighted mean of the standard deviation measures of the plurality of subbands using magnitude values of the subband signal 95 in the current time frame 75 as an energy measure.
In an advantageous embodiment, the variation determiner 275 smoothens the averaged standard deviation measure, when determining the first variation 290a, over the current, a plurality of previous and a plurality of future time frames. The smoothing as weighted according to an energy calculated using corresponding time frames and a windowing function. Furthermore, the variation determiner 275 is configured for smoothing the standard deviation measure, when determining the second variation 290b over the current, a plurality of previous, and a plurality of future time frames 75, wherein the smoothing is weighted according to the energy calculated using corresponding time frames 75 and a windowing function. Therefore, the variation comparator 280 compares the smoothened average standard deviation measure as the first variation 290a determined using the first variation mode and compares the smoothened standard deviation measure as the second variation 290b determined using the second variation mode.
An advantageous embodiment is depicted in
The second processing path comprises a PDF calculator 300b for calculating a phase derivative over frequency 305b from the audio signal 55 or a phase of the audio signal. A circular standard deviation calculator 310b forms a standard deviation measures 335b of the phase derivative over frequency 305. The standard deviation measure 305 is smoothened by a smoother 340b to form a smooth standard deviation measure 345b. The smoothened average standard deviation measures 345a and the smoothened standard deviation measure 345b are the first and the second variation, respectively. The variation comparator 280 compares the first and the second variation and the correction data calculator 285 calculates the phase correction data 295 based on the comparing of the first and the second variation.
Further embodiments show the calculator 270 handling three different phase correction modes. A figurative block diagram is shown in
The variation comparator 280 has to determine a suitable correction mode based on three variations. Based on this decision, the correction data calculator 285 calculates the phase correction data 295 in accordance with a third variation mode if a transient is detected. Furthermore, the correction data calculator 85 calculates the phase correction data 295 in accordance with a first variation mode, if an absence of a transient is detected and if the first variation 290a, determined in the first variation mode, is smaller or equal than the second variation 290b, determined in the second variation mode. Accordingly, the phase correction data 295 is calculated in accordance with the second variation mode, if an absence of a transient is detected and if the second variation 290b, determined in the second variation mode, is smaller than the first variation 290a, determined in the first variation mode.
The correction data calculator is further configured for calculating the phase correction data 295 for the third variation 290c for a current, one or more previous and one or more future time frames. Accordingly, the correction data calculator 285 is configured for calculating the phase correction data 295 for the second variation mode 290b for a current, one or more previous and one or more future time frames. Furthermore, the correction data calculator 285 is configured for calculating correction data 295 for a horizontal phase correction and the first variation mode, calculating correction data 295 for a vertical phase correction in the second variation mode, and calculating correction data 295 for a transient correction in the third variation mode.
In other words, the PDT of the violin is smooth over time whereas the PDF of the trombone is smooth over frequency. Hence, the standard deviation (STD) of these measures as a measure of the variation can be used to select the appropriate correction method. The STD of the phase derivative over time can be computed as
Xstdt1(k,n)=circstd{Xpdt(k,n+l)},−23≤l≤0,
Xstdt2(k,n)=circstd{Xpdt(k,n+l)},0≤l≤23,
Xstdt(k n)=min{Xstdt1(k,n),Xstdt2(k,n)}, (27)
and the STD of the phase derivative over frequency as
Xstdf(n)=circstd{Xpdf(k,n)},2≤k≤13, (28)
where circstd{ } denotes computing circular STD (the angle values could potentially be weighted by energy in order to avoid high STD due to noisy low-energy bins, or the STD computation could be restricted to bins with sufficient energy). The STDs for the violin and the trombone are shown in
The used correction method for each temporal frame is selected based on which of the STDs is lower. For that, Xstdt(k,n) values have to be combined over frequency. The merging is performed by computing an energy-weighted mean for a predefined frequency range
The deviation estimates are smoothened over time in order to have smooth switching, and thus to avoid potential artifacts. The smoothing is performed using a Hann window and it is weighted by the energy of the temporal frame
where W(l) is the window function and
is the sum of Xmag(k,n) over frequency. A corresponding equation is used for smoothing Xstdf(n).
The phase-correction method is determined by comparing Xsmstdt(n) and Xsmstdf(n). The default method is PDT (horizontal) correction, and if Xsmstdf(n)<Xsmstdt(n), PDF (vertical) correction is applied for the interval [n−5, n+5]. If both of the deviations are large, e.g. larger than a predefined threshold value, neither of the correction methods is applied, and bit-rate savings could be made.
8.4 Transient Handling—Phase Derivative Correction for Transients
The violin signal with a hand clap added in the middle is presented
The solution to the problem is straightforward. First, the transients are detected using a simple energy-based method. The instant energy of mid/high frequencies is compared to a smoothened energy estimate. The instant energy of mid/high frequencies is computed as
The smoothing is performed using a first-order IIR filter
Xsmmagmh(n)=0.1·Xmagmh(n)+0.9·Xsmmagmh(n−1). (32)
If Xmagmh(n)/Xsmmagmh(n)>θ, a transient has been detected. The threshold θ can be fine-tuned to detect the desired amount of transients. For example, θ=2 can be used. The detected frame is not directly selected to be the transient frame. Instead, the local energy maximum is searched from the surrounding of it. In the current implementation the selected interval is [n−2, n+7]. The temporal frame with the maximum energy inside this interval is selected to be the transient.
In theory, the vertical correction mode could also be applied for transients. However, in the case of transients, the phase spectrum of the baseband often does not reflect the high frequencies. This can lead to pre- and post-echoes in the processed signal. Thus, slightly modified processing is suggested for the transients.
The average PDF of the transient at high frequencies is computed
Xavghipdf(n)=circmean{Xpdf(k,n)},−11≤k≤36. (33)
The phase spectrum for the transient frame is synthesized using this constant phase change as in Eq. 24, but Xavgpdf(n) is replaced by Xavghipdf(n). The same correction is applied to the temporal frames within the interval [n−2, n+2] (π is added to the PDF of the frames n−1 and n+1 due to the properties of the QMF, see Section 6). This correction already produces a transient to a suitable position, but the shape of the transient is not necessarily as desired, and significant side lobes (i.e., additional transients) can be present due to the considerable temporal overlap of the QMF frames. Hence, the absolute phase angle has to be correct, too. The absolute angle is corrected by computing the mean error between the synthesized and the original phase spectrum. The correction is performed separately for each temporal frame of the transient.
The result of the transient correction is presented in
9 Compression of the Correction Data
Section 8 showed that the phase errors can be corrected, but the adequate bit rate for the correction was not considered at all. This section suggests methods how to represent the correction data with low bit rate.
9.1 Compression of the PDT Correction Data—Creating the Target Spectrum for the Horizontal Correction
There are many possible parameters that could be transmitted to enable the PDT correction. However, since Dsmpdt(k,n) is smoothened over time, it is a potential candidate for low-bit-rate transmission.
First, an adequate update rate for the parameters is discussed. The value was updated only for every N frames and linearly interpolated in between. The update interval for good quality is about 40 ms. For certain signals a bit less is advantageous and for others a bit more. Formal listening tests would be useful for assessing an optimal update rate. Nevertheless, a relatively long update interval appear to be acceptable.
An adequate angular accuracy for Dsmpdt(k,n) was also studied. 6 bits (64 possible angle values) is enough for perceptually good quality. Furthermore, transmitting only the change in the value was tested. Often the values appear to change only a little, so uneven quantization can be applied to have more accuracy for small changes. Using this approach, 4 bits (16 possible angle values) was found to provide good quality.
The last thing to consider is an adequate spectral accuracy. As can be seen in
9.1.1 Using Frequency Estimation for Compressing PDT Correction Data
As discussed in Section 5, the phase derivative over time basically means the frequency of the produced sinusoid. The PDTs of the applied 64-band complex QMF can be transformed to frequencies using the following equation
The produced frequencies are inside the interval finter(k)=[fc(k)−fBW, fc(k)+fBW], where fc(k) is the center frequency of the frequency band k and fBW is 375 Hz. The result is shown in
The same plot can be applied to the direct copy-up Zfreq(k,n) and the corrected Zchfreq(k,n) SBR (see
Since the frequencies of Xfreq(k,n) are spaced by the same amount, the frequencies of all frequency bands can be approximated if the spacing between the frequencies is estimated and transmitted. In the case of harmonic signals, the spacing should be equal to the fundamental frequency of the tone. Thus, only a single value has to be transmitted for representing all frequency bands. In the case of more irregular signals, more values are needed for describing the harmonic behavior. For example, the spacing of the harmonics slightly increases in the case of a piano tone [14]. For simplicity, it is assumed in the following that the harmonics are spaced by the same amount. Nonetheless, this does not limit the generality of the described audio processing.
Thus, the fundamental frequency of the tone is estimated for estimating the frequencies of the harmonics. The estimation of fundamental frequency is a widely studied topic (e.g., see [14]). Therefore, a simple estimation method was implemented to generate data used for further processing steps. The method basically computes the spacings of the harmonics, and combines the result according to some heuristics (how much energy, how stable is the value over frequency and time, etc.). In any case, the result is a fundamental-frequency estimate for each temporal frame Xf
Here, the fundamental frequency Xf
Alternatively, the fundamental frequency could be estimated in the decoding stage, and no information has to be transmitted. However, better estimates can be expected if the estimation is performed with the original signal in the encoding stage.
The decoder processing begins by obtaining a fundamental-frequency estimate Xf
The frequencies of the harmonics can be obtained by multiplying it with an index vector
∀κ∃:Xharm(κ,n)=κ·Xf
The result is depicted in
The transmitted parameter of the algorithm is the fundamental frequency Xf
The next step of the algorithm is to find a suitable value for each frequency band. This is performed by selecting the value of Xharm(κ,n) which is closest to the center frequency of each band fc(k) to reflect that band. If the closest value is outside the possible values of the frequency band (finter(k)), the border value of the band is used. The resulting matrix Xehfreq(k,n) contains a frequency for each time-frequency tile.
The final step of the correction-data compression algorithm is to convert the frequency data back to the PDT data
where mod( ) denotes the modulo operator. The actual correction algorithm works as presented in Section 8.1. Zthpdt(k,n) in Eq. 16a is replaced by Xehpdt(k,n) as the target PDT, and Eqs. 17-19 are used as in Section 8.1. The result of the correction algorithm with compressed correction data is shown in
Embodiments use more accuracy for low frequencies and less for high frequencies, using the total of 12 bits for each value. The resulting bit rate is about 0.5 kbps (without any compression, such as entropy coding). This accuracy produces equal perceived quality as no quantization. However, significantly lower bit rate can probably be used in many cases producing good enough perceived quality.
One option for low-bit-rate schemes is to estimate the fundamental frequency in the decoding phase using the transmitted signal. In this case no values have to be transmitted. Another option is to estimate the fundamental frequency using the transmitted signal, compare it to the estimate obtained using the broadband signal, and to transmit only the difference. It can be assumed that this difference could be represented using very low bit rate.
9.2 Compression of the PDF Correction Data
As discussed in Section 8.2, the adequate data for the PDF correction is the average phase error of the first frequency patch Davgpha(n). The correction can be performed for all frequency patches with the knowledge of this value, so the transmission of only one value for each temporal frame may be used. However, transmitting even a single value for each temporal frame can yield too high a bit rate.
Inspecting
Hence, the PDF (or the location of a transient) can be transmitted only sparsely in time, and the PDF behavior in between these time instants could be estimated using the knowledge of the fundamental frequency. The PDF correction can be performed using this information. This idea is actually dual to the PDT correction, where the frequencies of the harmonics are assumed to be equally spaced. Here, the same idea is used, but instead, the temporal locations of the transients are assumed to be equally spaced. A method is suggested in the following that is based on detecting the positions of the peaks in the waveform, and using this information, a reference spectrum is created for phase correction.
9.2.1 Using Peak Detection for Compressing PDF Correction Data—Creating the Target Spectrum for the Vertical Correction
The positions of the peaks have to be estimated for performing successful PDF correction. One solution would be to compute the positions of the peaks using the PDF value, similarly as in Eq. 34, and to estimate the positions of the peaks in between using the estimated fundamental frequency. However, this approach would involve a relatively stable fundamental-frequency estimation. Embodiments show a simple, fast to implement, alternative method, which shows that the suggested compression approach is possible.
A time-domain representation of the trombone signal is shown in
Using the transmitted metadata, a time-domain signal is created, which consists of impulses in the positions of the estimated peaks (see
The waveform of signals having vertical phase coherence is typically peaky and reminiscent of a pulse train. Thus, it is suggested that the target phase spectrum for the vertical correction can be estimated by modeling it as the phase spectrum of a pulse train that has peaks at corresponding positions and a corresponding fundamental frequency.
The position closest to the center of the temporal frame is transmitted for, e.g., every 20th temporal frame (corresponding to an interval of −27 ms). The estimated fundamental frequency, which is transmitted with equal rate, is used to interpolate the peak positions in between the transmitted positions.
Alternatively, the fundamental frequency and the peak positions could be estimated in the decoding stage, and no information has to be transmitted. However, better estimates can be expected if the estimation is performed with the original signal in the encoding stage.
The decoder processing begins by obtaining a fundamental-frequency estimate Xf
Ztvpha(k,n)=Xevpha(k,n). (37)
The suggested method uses the encoding stage to transmit only the estimated peak positions and the fundamental frequencies with the update rate of, e.g., 27 ms. In addition, it should be noted that errors in the vertical phase derivate are perceivable only when the fundamental frequency is relatively low. Thus, the fundamental frequency can be transmitted with a relatively low bit rate.
The result of the correction algorithm with compressed correction data is shown in
9.3 Compression of the Transient Handling Data
As transients can be assumed to be relatively sparse, it can be assumed that this data could be directly transmitted. Embodiments show transmitting six values per transient: one value for the average PDF, and five values for the errors in the absolute phase angle (one value for each temporal frame inside the interval [n−2, n+2]). An alternative is to transmit the position of the transient (i.e. one value) and to estimate the target phase spectrum Xetpha(k,n) as in the case of the vertical correction.
If the bit rate needed to be compressed for the transients, similar approach could be used as for the PDF correction (see Section 9.2). Simply the position of the transient could be transmitted, i.e., a single value. The target phase spectrum and the target PDF could be obtained using this location value as in Section 9.2.
Alternatively, the transient position could be estimated in the decoding stage and no information has to be transmitted. However, better estimates can be expected if the estimation is performed with the original signal in the encoding stage.
All of the previously described embodiments may be seen separately from the other embodiments or in a combination of embodiments. Therefore,
Accordingly, the decoder 110″ comprises a third target spectrum generator 65c, wherein the third target spectrum generator 65c generates a target spectrum for a third time frame of the subband of the audio signal 32 using third correction data 295c. Furthermore, the decoder 110″ comprises a third phase corrector 70c for correcting a phase 45 of the subband signal and the time frame of the audio signal 32 determined with a third phase correction algorithm, wherein the correction is performed by reducing a difference between a measure of the time frame of the subband of the audio signal and the target spectrum 85c. The audio subband signal calculator 350 can calculate the audio subband signal for a third time frame different from the first and the second time frames using the phase correction of the third phase corrector.
According to an embodiment, the first phase corrector 70a is configured for storing a phase corrected subband signal 91a of a previous time frame of the audio signal or for receiving a phase corrected subband signal of the previous time frame 375 of the audio signal from a second phase corrector 70b of the third phase corrector 70c. Furthermore, the first phase corrector 70a corrects the phase 45 of the audio signal 32 in a current time frame of the audio subband signal based on the stored or the received phase corrected subband signal of the previous time frame 91a, 375.
Further embodiments show the first phase corrector 70a performing a horizontal phase correction, the second phase corrector 70b performing a vertical phase correction, and the third phase corrector 70c performing a phase correction for transients.
From another point of view,
A second demultiplexer 130 (DEMUX) first divides the received metadata 135 into activation data 365 and correction data 295a-c for the different correction modes. Based on the activation data, the computation of the target spectrum is activated for the right correction mode (others can be idle). Using the target spectrum, the phase correction is performed to the received BWE signal using the desired correction mode. It should be noted that as the horizontal correction 70a is performed recursively (in other words: dependent on previous signal frames), it receives the previous correction matrices also from other correction modes 70b, c. Finally, the corrected signal, or the unprocessed one, is set to the output based on the activation data.
After having corrected the phase data, the underlying BWE synthesis further downstream is continued, in the case of the current example the SBR synthesis. Variations might exist where exactly the phase correction is inserted into the BWE synthesis signal flow. Advantageously, the phase-derivative correction is done as an initial adjustment on the raw spectral patches having phases Zpha(k,n) and all additional BWE processing or adjustment steps (in SBR this can be noise addition, inverse filtering, missing sinusoids, etc.) are executed further downstream on the corrected phases Zcpha(k,n).
Many other embodiments can be thought of where the signal processor blocks are switched. For example, the magnitude processor 125′ and the block A may be swapped. Therefore, the block A works on the reconstructed audio signal 35, where the magnitude values of the patches have already been corrected. Alternatively, the audio subband signal calculator 350 may be located after the magnitude processor 125′ in order to form the corrected audio signal 355 from the phase corrected and the magnitude corrected part of the audio signal.
Furthermore, the decoder 110″ comprises a synthesizer 100 for synthesizing the phase and magnitude corrected audio signal to obtain the frequency combined processed audio signal 90. Optionally, since neither the magnitude nor the phase correction is applied on the core decoded audio signal 25, said audio signal may be transmitted directly to the synthesizer 100. Any optional processing block applied in one of the previously described decoders 110 or 110′ may be applied in the decoder 110″ as well.
According to embodiments, the calculator 270 comprises a set of correction data calculators 285a-c for correcting the phase correction in accordance with a first variation mode, a second variation mode, or a third variation mode. Furthermore, the calculator 270 determines activation data 365 for activating one correction data calculator of the set of correction data calculators 285a-c. The output signal former 170 forms the output signal comprising the activation data, the parameters, the core encoded audio signal, and the phase correction data.
Embodiments show the calculator 270 comprising a metadata former 390, which forms a metadata stream 295′ comprising the calculated correction data 295a, 295b, or 295c and the activation data 365. The activation data 365 may be transmitted to the decoder if the correction data itself does not comprise sufficient information of the current correction mode. Sufficient information may be for example a number of bits used to represent the correction data, which is different for the correction data 295a, the correction data 295b, and the correction data 295c. Furthermore, the output signal former 170 may additionally use the activation data 365, such that the metadata former 390 can be neglected.
From another point of view, the block diagram of
The correction-mode-computation block first computes the correction mode that is applied for each temporal frame. Based on the activation data 365, correction-data 295a-c computation is activated in the right correction mode (others can be idle). Finally, multiplexer (MUX) combines the activation data and the correction data from the different correction modes.
A further multiplexer (not depicted) merges the phase-derivative correction data into the bit stream of the BWE and the perceptual encoder that is being enhanced by the inventive correction.
The methods 5800 and 5900 as well as the previously described methods 2300, 2400, 2500, 3400, 3500, 3600 and 4200, may be implemented in a computer program to be performed on a computer.
It has to be noted that the audio signal 55 is used as a general term for an audio signal, especially for the original i.e. unprocessed audio signal, the transmitted part of the audio signal Xtrans(k,n) 25, the baseband signal Xbase(k,n) 30, the processed audio signal comprising higher frequencies 32 when compared to the original audio signal, the reconstructed audio signal 35, the magnitude corrected frequency patch Y(k,n,i) 40, the phase 45 of the audio signal, or the magnitude 47 of the audio signal. Therefore, the different audio signals may be mutually exchanged due to the context of the embodiment.
Alternative embodiments relate to different filter bank or transform domains used for the inventive time-frequency processing, for example the short time Fourier transform (STFT) a Complex Modified Discrete Cosine Transform (CMDCT), or a Discrete Fourier Transform (DFT) domain. Therefore, specific phase properties related to the transform may be taken into consideration. In detail, if e.g. copy-up coefficients are copied from an even number to an odd number or vice versa, i.e. the second subband of the original audio signal is copied to the ninth subband instead of the eighth subband as described in the embodiments, the conjugate complex of the patch may be used for the processing. The same applies to a mirroring of the patches instead of using e.g. the copy-up algorithm, to overcome the reversed order of the phase angles within a patch.
Other embodiments might resign side information from the encoder and estimate some or all useful correction parameters on decoder site. Further embodiments might have other underlying BWE patching schemes that for example use different baseband portions, a different number or size of patches or different transposition techniques, for example spectral mirroring or single side band modulation (SSB). Variations might also exist where exactly the phase correction is concerted into the BWE synthesis signal flow. Furthermore, the smoothing is performed using a sliding Hann window, which may be replaced for better computational efficiency by, e.g. a first-order IIR.
The use of state of the art perceptual audio codecs often impairs the phase coherence of the spectral components of an audio signal, especially at low bit rates, where parametric coding techniques like bandwidth extension are applied. This leads to an alteration of the phase derivative of the audio signal. However, in certain signal types the preservation of the phase derivative is important. As a result, the perceptual quality of such sounds is impaired. The present invention readjusts the phase derivative either over frequency (“vertical”) or over time (“horizontal”) of such signals if a restoration of the phase derivative is perceptually beneficial. Further, a decision is made whether adjusting the vertical or horizontal phase derivative is perceptually advantageous. The transmission of only very compact side information is needed to control the phase derivative correction processing. Therefore, the invention improves sound quality of perceptual audio coders at moderate side information costs.
In other words, spectral band replication (SBR) can cause errors in the phase spectrum. The human perception of these errors was studied revealing two perceptually significant effects: differences in the frequencies and the temporal positions of the harmonics. The frequency errors appear to be perceivable only when the fundamental frequency is high enough that there is only one harmonic inside an ERB band. Correspondingly, the temporal-position errors appear to be perceivable only if the fundamental frequency is low and if the phases of the harmonics are aligned over frequency.
The frequency errors can be detected by computing the phase derivative over time (PDT). If the PDT values are stable over time, differences in them between the SBR-processed and the original signals should be corrected. This effectively corrects the frequencies of the harmonics, and thus, the perception of inharmonicity is avoided.
The temporal-position errors can be detected by computing the phase derivative over frequency (PDF). If the PDF values are stable over frequency, differences in them between the SBR-processed and the original signals should be corrected. This effectively corrects the temporal positions of the harmonics, and thus, the perception of modulating noises at the cross-over frequencies is avoided.
Although the present invention has been described in the context of block diagrams where the blocks represent actual or logical hardware components, the present invention can also be implemented by a computer-implemented method. In the latter case, the blocks represent corresponding method steps where these steps stand for the functionalities performed by corresponding logical or physical hardware blocks.
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 are advantageously 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.
| Number | Date | Country | Kind |
|---|---|---|---|
| 14175202 | Jul 2014 | EP | regional |
| 15151478 | Jan 2015 | EP | regional |
This application is a continuation of copending U.S. patent application Ser. No. 15/392,776, filed Dec. 28, 2016, which is a continuation of copending International Application No. PCT/EP2015/064443, filed Jun. 25, 2015, which is incorporated herein in its entirety by this reference thereto, which claims priority from European Applications Nos. EP 14 175 202.2, filed Jul. 1, 2014, and EP 15 151 478.3, filed Jan. 16, 2015, which are each incorporated herein in its entirety by this reference thereto.
| Number | Name | Date | Kind |
|---|---|---|---|
| 4802225 | Patterson | Jan 1989 | A |
| 5001758 | Galand et al. | Mar 1991 | A |
| 5142584 | Ozawa | Aug 1992 | A |
| 5794186 | Bergstrom et al. | Aug 1998 | A |
| 5809459 | Bergstrom et al. | Sep 1998 | A |
| 6226661 | Savell | May 2001 | B1 |
| 6292775 | Holmes | Sep 2001 | B1 |
| 6701297 | Main | Mar 2004 | B2 |
| 6745155 | Andringa et al. | Jun 2004 | B1 |
| 6980933 | Cheng et al. | Dec 2005 | B2 |
| 7146503 | Johnston et al. | Dec 2006 | B1 |
| 7676336 | Herre et al. | Mar 2010 | B2 |
| 8036891 | Fulop | Oct 2011 | B2 |
| 8255228 | Hilpert et al. | Aug 2012 | B2 |
| 8886499 | Matsumoto | Nov 2014 | B2 |
| 9240196 | Disch et al. | Jan 2016 | B2 |
| 9424847 | Ishikawa et al. | Aug 2016 | B2 |
| 9881624 | Choo et al. | Jan 2018 | B2 |
| 10140997 | Disch et al. | Nov 2018 | B2 |
| 10192561 | Disch et al. | Jan 2019 | B2 |
| 10529346 | Disch et al. | Jan 2020 | B2 |
| 20030004697 | Ferris | Jan 2003 | A1 |
| 20040225505 | Andersen et al. | Nov 2004 | A1 |
| 20050165587 | Cheng et al. | Jul 2005 | A1 |
| 20050171785 | Nomura et al. | Aug 2005 | A1 |
| 20050177360 | Schuijers et al. | Aug 2005 | A1 |
| 20060014299 | Troup | Jan 2006 | A1 |
| 20060143000 | Setoguchi | Jun 2006 | A1 |
| 20060147056 | Klayman | Jul 2006 | A1 |
| 20070019746 | Lin et al. | Jan 2007 | A1 |
| 20070027678 | Hotho et al. | Feb 2007 | A1 |
| 20070067162 | Villemoes et al. | Mar 2007 | A1 |
| 20070094009 | Ryu et al. | Apr 2007 | A1 |
| 20070100610 | Disch et al. | May 2007 | A1 |
| 20070233466 | Tammi | Oct 2007 | A1 |
| 20070238415 | Sinha et al. | Oct 2007 | A1 |
| 20080052068 | Aguilar et al. | Feb 2008 | A1 |
| 20080154615 | Oomen | Jun 2008 | A1 |
| 20080170711 | Breebaart et al. | Jul 2008 | A1 |
| 20080235034 | Lee et al. | Sep 2008 | A1 |
| 20080243493 | Rault et al. | Oct 2008 | A1 |
| 20080255828 | Chesnutt et al. | Oct 2008 | A1 |
| 20090006103 | Koishida et al. | Jan 2009 | A1 |
| 20090012782 | Geiser et al. | Jan 2009 | A1 |
| 20090299756 | Davis et al. | Dec 2009 | A1 |
| 20090304198 | Herre et al. | Dec 2009 | A1 |
| 20090326942 | Fulop | Dec 2009 | A1 |
| 20100010808 | Sugiyama et al. | Jan 2010 | A1 |
| 20100063811 | Gao | Mar 2010 | A1 |
| 20100198588 | Sudo et al. | Aug 2010 | A1 |
| 20100286805 | Gao et al. | Nov 2010 | A1 |
| 20110004478 | Jax et al. | Jan 2011 | A1 |
| 20110046964 | Moon et al. | Feb 2011 | A1 |
| 20110132179 | Saino | Jun 2011 | A1 |
| 20110173005 | Hilpert et al. | Jul 2011 | A1 |
| 20110173006 | Nagel et al. | Jul 2011 | A1 |
| 20110206209 | Ojala | Aug 2011 | A1 |
| 20110206223 | Ojala | Aug 2011 | A1 |
| 20110216918 | Nagel et al. | Sep 2011 | A1 |
| 20110246205 | Lin et al. | Oct 2011 | A1 |
| 20110280421 | Macours | Nov 2011 | A1 |
| 20110288873 | Nagel et al. | Nov 2011 | A1 |
| 20110305352 | Villemoes et al. | Dec 2011 | A1 |
| 20120010880 | Nagel et al. | Jan 2012 | A1 |
| 20120020499 | Neusinger et al. | Jan 2012 | A1 |
| 20120239388 | Sverrisson | Sep 2012 | A1 |
| 20120278088 | Villemoes | Nov 2012 | A1 |
| 20120303362 | Duni et al. | Nov 2012 | A1 |
| 20130010985 | Disch et al. | Jan 2013 | A1 |
| 20130058498 | Disch et al. | Mar 2013 | A1 |
| 20130089215 | Kon | Apr 2013 | A1 |
| 20130114817 | Wu et al. | May 2013 | A1 |
| 20130117029 | Liu et al. | May 2013 | A1 |
| 20130166286 | Matsumoto | Jun 2013 | A1 |
| 20130173273 | Kuntz et al. | Jul 2013 | A1 |
| 20130262096 | Wilhelms-tricarico et al. | Oct 2013 | A1 |
| 20130304480 | Kuntz et al. | Nov 2013 | A1 |
| 20140088978 | Mundt et al. | Mar 2014 | A1 |
| 20140214413 | Atti et al. | Jul 2014 | A1 |
| 20140297295 | Villemoes et al. | Oct 2014 | A1 |
| 20140372131 | Disch et al. | Dec 2014 | A1 |
| 20150078571 | Kurylo et al. | Mar 2015 | A1 |
| 20150149156 | Atti et al. | May 2015 | A1 |
| 20150230041 | Fejzo et al. | Aug 2015 | A1 |
| 20150340045 | Hardwick et al. | Nov 2015 | A1 |
| 20160006453 | Jalali et al. | Jan 2016 | A1 |
| 20160118056 | Choo et al. | Apr 2016 | A1 |
| 20160134985 | Hashimoto et al. | May 2016 | A1 |
| 20160291056 | Pickerd et al. | Oct 2016 | A1 |
| 20170110133 | Disch et al. | Apr 2017 | A1 |
| 20170110134 | Disch et al. | Apr 2017 | A1 |
| 20170110135 | Disch et al. | Apr 2017 | A1 |
| 20170270937 | Nagel et al. | Sep 2017 | A1 |
| Number | Date | Country |
|---|---|---|
| 1217111 | May 1999 | CN |
| 1527995 | Sep 2004 | CN |
| 1647155 | Jul 2005 | CN |
| 1669073 | Sep 2005 | CN |
| 1754205 | Mar 2006 | CN |
| 1781141 | May 2006 | CN |
| 1875402 | Dec 2006 | CN |
| 1917495 | Feb 2007 | CN |
| 1934618 | Mar 2007 | CN |
| 1950815 | Apr 2007 | CN |
| 1969487 | May 2007 | CN |
| 101046964 | Oct 2007 | CN |
| 101051456 | Oct 2007 | CN |
| 101086845 | Dec 2007 | CN |
| 101091209 | Dec 2007 | CN |
| 101373594 | Feb 2009 | CN |
| 101443843 | May 2009 | CN |
| 101641734 | Feb 2010 | CN |
| 101960515 | Jan 2011 | CN |
| 102089807 | Jun 2011 | CN |
| 102194457 | Sep 2011 | CN |
| 102334158 | Jan 2012 | CN |
| 102483921 | May 2012 | CN |
| 102027537 | Oct 2012 | CN |
| 102741921 | Oct 2012 | CN |
| 102985970 | Mar 2013 | CN |
| 103037290 | Apr 2013 | CN |
| 103258539 | Aug 2013 | CN |
| 103490678 | Jan 2014 | CN |
| 103548077 | Jan 2014 | CN |
| 103621110 | Mar 2014 | CN |
| 103632678 | Mar 2014 | CN |
| 1259955 | Jan 2006 | EP |
| 1903560 | Mar 2008 | EP |
| 1081684 | Apr 2010 | EP |
| 2019391 | Jan 2013 | EP |
| 2631906 | Aug 2013 | EP |
| 2720222 | Apr 2014 | EP |
| 2411976 | May 2014 | EP |
| 2545553 | Jul 2014 | EP |
| 2004-53895 | Feb 2004 | JP |
| 2012515362 | Jul 2012 | JP |
| 2013521536 | Jun 2013 | JP |
| 2013135433 | Jul 2013 | JP |
| 2015-508911 | Mar 2015 | JP |
| 2325046 | May 2008 | RU |
| 2407072 | Dec 2010 | RU |
| 2418324 | May 2011 | RU |
| 2452044 | May 2012 | RU |
| 201248619 | Dec 2012 | TW |
| 9732413 | Sep 1997 | WO |
| 1998057436 | Dec 1998 | WO |
| 03042979 | May 2003 | WO |
| 2005073960 | Aug 2005 | WO |
| 2007068861 | Jun 2007 | WO |
| 2009008068 | Jan 2009 | WO |
| 2010108895 | Sep 2010 | WO |
| 2012025283 | Mar 2012 | WO |
| 2012131438 | Oct 2012 | WO |
| 2013124445 | Aug 2013 | WO |
| 2013127801 | Sep 2013 | WO |
| Entry |
|---|
| Daudet, Laurent et al., “MDCT Analysis of Sinusoids: Exact Results and Applications to Coding Artifacts Reduction”, IEEE Transactions on Speech and Audio Processing, vol. 12, No. 3, May 2004, 302-312. |
| Dietz, Martin et al., “Spectral Band Replication, a Novel Approach in Audio Coding”, Audio Engineering Society Convention Paper 5553 Presented at the 112th Convention. XP009020921, May 2002, 1-8. |
| Dorran, David et al., “Time-Scale Modification of Music Using a Synchronized Subband/Time-Domain Approach”, IEEE International Conference on Acoustics, Speech, and Signal Processing, May 2004, IV225-IV228. |
| Ekstrand, Per , “Bandwidth Extension of Audio Signals by Spectral Band Replication”, Proc 1st IEEE Benelux Workshop on Model Based Processing and Coding of Audio (MPCA-2002) XP-000962047, Nov. 15, 2002, 53-58. |
| Fitz, Kelly R. et al., “A Unified Theory of Time-Frequency Reassignment”, arXiv preprint arXiv:0903.3080, Mar. 2009, 1-38. |
| Fulop, Sean A. et al., “Algorithms for computing the time-corrected instantaneous frequency (reassigned) spectrogram, with applications”, The Journal of the Acoustical Society of America 119.1, 01.06, pp. 360-371. |
| Gerzon, Michael A. , “The Presentation of Statistical Information on a Circle”, Journal of the Audio Engineering Society 22.7, 1974, p. 536. |
| Griesinger, David , “The Relationship Between Audience Engagement and the Ability to Perceive Pitch, Timbre, Azimuth and Envelopment of Multiple Sources”, 26th Tonmeistertagung-VDT International Convention, Nov. 2010, 456-480. |
| Jaillet, F. et al., “On the structure of the phase around the zeros of the short-time Fourier transform”, Proceedings of NAG/DAGA International Conference on Acoustics, Rotterdam, Netherlands., Mar. 2009, pp. 1584-1587. |
| Kim, Junghoe et al., “Enhanced Stereo Coding with Phase Parameters for MPEG Unified Speech and Audio Coding”, Audio Engineering Society Convention 127. Oct. 12, 2009 http://www.aes.org/e-lib/browse.cfm?elib=15070, Oct. 1, 2009. |
| Kim, Kijun et al., “Improvement in Parametric High-Band Audio Coding by Controlling Temporal Envelope With Phase Parameter”, Audio Engineering Society Convention Paper 8982 Presented at the 135th Convention, New York, Oct. 2013, 1-7. |
| Klapuri, Anssi P. , “Multiple Fundamental Frequency Estimation Based on Harmonicity and Spectral Smoothness”, IEEE Transactions on Speech and Audio Processing, vol. 11, No. 6, Nov. 2003, 804-816. |
| Laitinen, Mikko-Ville et al., “Sensitivity of Human Hearing to Changes in Phase Spectrum”, J. Audio Eng. Soc., vol. 61, No. 11 XP40633294A, Nov. 2013, 860-877. |
| Laroche, Jean , “Frequency-Domain Techniques for High-Quality Voice Modification”, Proc. of the 6th Int. Conference on Digital Audio Effects (DAFx-03), London, UK, Sep. 2003, 322-328. |
| Laroche, Jean et al., “Phase-Vocoder: About This Phasiness Business”, Applications of Signal Processing to Audio and Acoustics. IEEE ASSP Workshop, 1997, 19-22. |
| Larsen, Erik et al., “Audio Bandwidth Extension Application of Psychoacoustics, Signal Processing and Loudspeaker Design”, John Wiley and Sons, Ltd. Chapters 5 and 6, 2004, C1-91. |
| Moore, Brian C. et al., “Suggested Formulae for Calculating Audiotry-Filter Bandwidths and Excitation Patterns”, The Journal of the Acoustical Society of America 74, 750 (1983); doi: 10.1121/1.389861, Sep. 1983, 750-753. |
| Nagel, Frederik et al., “A Phase Vocoder Driven Bandwidth Extension Method With Novel Transient Handling for Audio Codecs”, Audio Engineering Society Convention Paper Presented at the 126th Convention, Munich, Germany, May 2008, 1-8. |
| Nelson, Douglas J. , “Cross-spectral Methods for Processing Speech”, The Journal of the Acoustical Society of America 110.5, Nov. 2001, 2575-2592. |
| Painter, Ted et al., “Perceptual Coding of Digital Audio”, Proceedings of the IEEE, vol. 88, No. 4, Apr. 2000, 451-513. |
| Robel, Axel , “A New Approach to Transient Processing in the Phase Vocoder”, 6th International Conference on Digital Audio Effects (DAFx), Sep. 2003, 344-349. |
| Shackleton, Trevor M. et al., “The Role of Resolved and Unresolved Harmonics in Pitch Perception and Frequency Modulation Dicrimination”, J. Acoust. Soc. Am. 95, Jun. 1994, 3529-3540. |
| Stylianou, Yannis, “Applying the Harmonic Plus Noise Model in Concatenative Speech Synthesis”, IEEE Transactions on Speech and Audio Processing, vol. 9, No. 1, Jan. 2001, pp. 21-29, Jan. 2001, pp. 21-29. |
| Sean Fulop, “The Reassigned Spectrogram”, Speech Spectrum Analysis, Signals and Communication Technology, DOI:10,1007/978-3-642-17478-0_6, May 2011, pp. 127-163. |
| Philippe Guillemain, et al. “Characterization of Acoustic Signals Through Continuous Linear Time-Frequency Representations.” Proceedings of the IEEE vol. 84, No. 4, Apr. 1996, pp. 561-585. |
| Cheng Liu, et al., “Instantaneous frequency estimation by the reassigned Stockwell spectrogram”, 11th International Conference on Information Science, Signal Processing and their Applications (ISSPA). IEEE, Jul. 2012, pp. 1235-1240. |
| Axel Roebel, “Transient detection and preservation in the phase vocoder”, International Computer Music Conference (ICMG), Oct. 2003, Singapore, pp. 247-250. |
| “Group delay and phase delay”, Wikipedia, 6 pages, downloaded May 4, 2020. (Year: 2020), 6 pp. |
| Enderlein, J. , et al., “Two channel SAW spread spectrum transmission system for audio signals”, IEEE Ultrasonics Symposium, 1996, 1996. |
| Number | Date | Country | |
|---|---|---|---|
| 20190156842 A1 | May 2019 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | 15392776 | Dec 2016 | US |
| Child | 16258604 | US | |
| Parent | PCT/EP2015/064443 | Jun 2015 | US |
| Child | 15392776 | US |
