Noise reduction with integrated tonal noise reduction

Abstract

The system provides a technique for suppressing or eliminating tonal noise in and input signal. The system operates on the input signal at a plurality of frequency bins and uses information generated at a prior bin to assist in calculating values at subsequent bins. The system first identifies peaks in a signal and then determines if the peaks are from tonal effects. This can be done by comparing the estimated background noise of a current bin to the smoothed background noise of the same bin. The smoothed background noise can be calculated using an asymmetric IIR filter. When the ratio of the current background noise estimate to the currently calculated smoothed background noise is far greater than 1, tonal noise is assumed. When tonal noise is found, a number of suppression techniques can be applied to reduce the tonal noise, including gain suppression with fixed floor factor, an adaptive floor factor gain suppression technique, and a random phase technique.

Description

BACKGROUND OF THE SYSTEM

1. Technical Field

The system is directed to the field of sound processing. More particularly, this system provides a way to remove tonal noise without degrading speech or music.

2. Related Art

Speech enhancement often involves the removal of noise from a speech signal. It has been a challenging topic of research to enhance a speech signal by removing extraneous noise from the signal so that the speech may be recognized by a speech processor or by a listener. Various approaches have been developed over the past decades. Among them the spectral subtraction methods are the most widely used in real-time applications. In this method, an average noise spectrum is estimated and subtracted from the noisy signal spectrum, so that average signal-to-noise ratio (SNR) is improved.

However, prior art speech enhancement techniques do not always work when the noise is of a type referred to as “tonal” noise. Tonal noise can occur in homes, offices, cars, and other environments. An often quoted source of tonal noise in the home and office is the buzzing of fluorescent lights. Another is the hum of a computer or projector fan. In the car tonal noise can result from rumble strips, car engine, alternator whine, radio interference (“GSM buzz”), or a whistle from an open window. This tonal noise can negatively impact phone conversations and speech recognition, making speech a little more difficult to understand or recognize.

A speech processing system which examines an input signal for desired signal content may interpret the tonal noise as speech, may isolate a segment of the input signal with the tonal noise, and may attempt to process the tonal noise. The speech processing system consumes valuable computational resources not only to isolate the segment, but also to process the segment and take action based on the result of the processing. In a speech recognition system, the system may interpret the tonal noise as a voice command, execute the spurious command, and responsively take actions that were never intended.

Tonal noise appears as constant peaks in an acoustic frequency spectrum. By definition the peaks stand out from the broader band noise, often by 6 to 20 dB. Noise reduction typically attenuates all frequencies equally, so the remaining tonal noise is quieter, but is just as distinct after noise reduction as before. Therefore the existing noise removal approach does not really help reduce tonal noise relative to the broader background noise.

SUMMARY

The invention details an improvement to a noise removal system. Quasi-stationary tonal noise appears as peaks in a spectrum of normally broadband or diffuse noise. Noise reduction typically attenuates all frequencies equally, so tonal noise while quieter is just as distinct before noise reduction as after. The system identifies peaks, determines which peaks are likely to be tonal peaks, and applies an adaptive suppression to the tonal peaks. The system uses a technique of tonal noise reduction (TNR) that places greater attenuation at frequencies where tonal noise is found. The TNR system may do additional processing (phase randomization) to virtually eliminate any residual tonal sound. This system is not a simple passive series of notch filters and therefore does not remove speech or music that overlaps in frequencies. Moreover it is adaptive and does not do any additional filtering if tonal noise is not present.

BRIEF DESCRIPTION OF THE DRAWINGS

The invention can be better understood with reference to the following drawings and description. The components in the Figures are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the invention. Moreover, in the Figures, like reference numerals designate corresponding parts throughout the different views.

FIG. 1 is a PSD of normal car noise.

FIG. 2 is a PSD of tonal noise.

FIG. 3 the PSD of the tonal noise after prior art noise reduction.

FIG. 4 illustrates the PSD of the tonal noise, processed by the disclosed tonal noise reduction method.

FIG. 5 is a flow diagram illustrating the operation of the system in identifying and suppressing tonal noise.

FIG. 6 is a flow diagram illustrating the technique used by the system to estimate the smoothed background noise.

FIG. 7 is a flow diagram illustrating a technique for determining the presence of tonal peaks.

FIG. 8 is a flow diagram illustrating prior art technique for estimating a clean speech signal.

FIG. 9 is a flow diagram illustrating the use of an adaptive factor to calculate a suppression gain value.

FIG. 10. is a flow diagram illustrating a suppression technique using random phases.

DETAILED DESCRIPTION OF THE SYSTEM

A typical frequency domain speech enhancement system usually consists of a spectral suppression gain calculation method, and a background noise power spectral density (PSD) estimation method. While spectral suppression is well understood, PSD noise estimation historically received less attention. However, it has been-found very important to the quality and intelligibility of the overall-system in recent years. Most spectral suppression methods can achieve good quality when background noise is stationary or semi-stationary over time and also smooth across frequencies. When tonal noise is present in the background a conventional spectral suppression method can suppress it, but cannot eliminate the tonal noise. The residual tonal noises are distinctive and can be annoying to the human ear. This system provides principles and techniques to remove the tonal noise completely without degrading speech quality.

Tonal noise reduction (TNR) of the system places greater-attenuation at the peak frequencies to the extent to which the peaks are greater than the diffuse noise. For example, if a peak is seen in a noise estimate that is 10 dB greater than the noise in the surrounding frequencies then an extra 10 dB of noise attenuation is done at that frequency. Thus, the spectral shape after TNR will be smooth across neighboring frequencies and tonal noise is significantly reduced.

At any given frequency the contribution of noise can be considered insignificant when the speech is greater than 12 dB above the noise. Therefore, when the signal is significantly higher than the noise, tonal or otherwise, NR, with or without TNR should not and does not have, any significant impact. Lower SNR signals will be attenuated more heavily around the tonal peaks, and those signals equal to the tonal noise peaks will be attenuated such that the resulting spectrum is flat around the peak frequency (its magnitude is equal to the magnitude of the noise in the neighboring frequencies).

Reducing the power of the tonal noise (while leaving its phase intact) may not completely remove the sound of the tones, because the phase at a given frequency still contributes to the perception of the tone. In one method, if the signal is close to the tonal noise, the phase at that frequency bin may be randomized. This has the benefit of completely removing the tone at that frequency. The system provides improved voice quality, reduced listener fatigue, and improved speech recognition.

Other systems, methods features and advantages of the invention will be, or will, become, apparent to one with skill in the art upon examination of the following figures and detailed description. It is intended that all such additional systems, methods, features and advantages be included within this description, be within the scope of the invention, and be protected by the following claims.

Methods to Detect Tonal Noise

Normal car noise is diffuse noise. Its power density smoothly decays when frequency increases. A spectrogram of normal car noise shows a relatively smooth and somewhat homogeneous distribution throughout the spectrogram. By contrast, tonal noise usually only covers certain frequencies and holds for a relative long period of time. A spectrogram of tonal noise shows a much uneven distribution.

A PSD of normal car noise is illustrated in FIG. 1. The graph shows how the power of a signal is distributed with frequency. As can be seen, normal road noise has more power at lower frequencies and has a substantially reduction in power with frequency so that at the higher frequencies, the power of the signal is relatively small. By contrast, the PSD of tonal noise, illustrated in FIG. 2, shows that the power is distributed in a number of peaks at varying frequencies. The PSD of the tonal noise signal of FIG. 2 is much more “peaky” than that of normal road noise.

Most conventional noise tracking algorithms with reasonable frequency resolution can track tonal noise in the background. Tonal noise usually shows in the noise spectrum as peaks standing much above their neighbors as illustrated at a number of frequencies in FIG. 2.

FIG. 5 is a flow diagram illustrating the operation of the system in identifying and suppressing tonal noise. At step 501 the system identifies the peaks of a background noise spectrum. At step 502 the tonal peaks that are to be suppressed are identified. At step 503, the tonal peaks are suppressed so that their impact on the signal is reduced.

Tonal Noise Peak Detection

It can be seen that to deal with tonal noise, one method is to first identify the peaks of tonal noise. FIG. 6 is a flow diagram illustrating the technique used by the system to identify peaks in an input signal. The system transforms the time domain signal into frequency domain. The frequency resolution may vary from systems to systems. In some embodiments of the system, the frequency resolution for this part of the system is 43 Hz per bin. The input signal is analyzed at each of the frequency bins. At step 601 the background noise estimate for a current bin under consideration is obtained. At step 602, the current background noise estimate is compared to the smoothed background noise for the prior bin (the bin analyzed just prior to the current bin). At decision block 603 it is determined if the current background noise estimate is greater than or equal to the smoothed background noise of the prior bin. If yes, a first algorithm is applied at step 604. If no, a second algorithm is applied at step 605.

One method for implementing the technique of FIG. 6 is the application of an asymmetric IIR (infinite impulse response) filter to detect the location as well as magnitude of tonal noise peaks.

As noted at step 601, the background noise estimate B_n(k) at n th frame and k th frequency bin is estimated. The smoothed background noise B_n(k) for this kth bin can be calculated by an asymmetric IIR filter. The background noise estimate B_n(k) of the present bin is compared to the smoothed background noise B_n(k−1) of the prior bin (step 602). Depending on the results of the comparison, different-branches of the asymmetrical IIR filter are applied.

when B_n(k)≧ B_n(k−1) (step 603 is true) the following is applied.B_n(k)=β₁*B_n(k)+(1−β₁)*B_n(k−1)(step 604)

when B_n(k)< B_n(k−1) (step 603 is false) then apply:B_n(k)=β₂*B_n(k)+(1−β₂)*B_n(k−1)(step 605)

Here β₁ and β₂ are two parameters in the range from 0 to 1. They are used to adjust the rise and fall adaptation speed. By choosing β₂ to be greater than or equal to β₁, the smoothed background noise follows closely to the noise estimation except at the places where there are tonal peaks. The smoothed background can then be used to remove tonal noise in the next step. Note that the same filter can be run through the noise spectrum in forward or reverse direction, and also for multiple passes as desired.

Identifying Tonal Noise Peaks

FIG. 7 is a flow diagram illustrating a ratio technique for determining the presence of tonal peaks. At step 701 the smoothed background noise for the current bin is calculated. (This can be done as described in FIG. 6). At step 702 the smoothed background noise of the current bin is compared to the background noise estimate of the current bin. At decision block 703 it is determined if the ratio is much greater than 1. If so, it is presumed that the peak at that bin is a tonal peak at step 704. If not, the peak at that bin is presumed to be normal noise at step 705.

One method for implementing the technique of FIG. 7 is described here. The ratio between non-smoothed (B_n(k)) and smoothed ( B_n(k)) (step 701) background noise is given by:ξ_n(k)=B_n(k)/B_n(k)(Step 702).

The value of ξ_n(k) is normally around 1 (step 703 is false) meaning the non-smoothed background noise is approximately equal to the smoothed background noise and is thus normal noise (step 705). However when there is tonal noise in the background, large values of ξ_n(k) are found (step 703 is true) at different frequencies. Therefore a large ξ_n(k) is used as an indicator of tonal noise (step 704).

The system tracks which bins have noise due to tonal effects and which bins have noise considered to bet normal noise.

Methods to Remove Tonal Noise

Non-Adaptive

Once the peaks that require processing have been determined, corrective action can be taken. FIG. 8 is a flow diagram illustrating a non-adaptive technique for estimating a clean speech signal. At step 801 the spectral magnitude of the noisy speech signal at the current bin is determined. At step 802 a suppression gain value is applied to the spectral magnitude. At step 803 an estimate of clean speech spectral magnitude is generated.

The system of FIG. 8 can be implemented as follows. In a classical additive noise model, noisy speech is given byy(t)=x(t)+d(t)

Where x(t) and d(t) denote the speech and the noise signal, respectively.

Let |Y_n,k|, |X_n,k|, and |D_n,k| designate the short-time spectral magnitude of noisy speech, speech and noise, respectively, at n th frame and k th frequency bin. The noisy speech spectral magnitude can be known (step 801), but the actual values of the noise and clean speech are not known. To obtain a cleaned up speech signal requires manipulation of the noisy speech spectral magnitude. The noise reduction process consists in the application (step 802) of a spectral gain value G_n,k to each short-time spectrum value. An estimate of the clean speech spectral-magnitude can be obtained (step 803) as:|{circumflex over (X)}_n,k|=G_n,k·|Y_n,k|

Where G_n,k is the spectral suppression gain. Various methods have been introduced in the literatures on how to calculate this gain. Examples include the decision-directed approach proposed in Ephraim, Y.; Malah, D.; Speech Enhancement Using A Minimum-Mean Square Error Short-Time Spectral Amplitude Estimator, IEEE Trans. on Acoustics, Speech, and Signal Processing Volume 32, Issue 6, December 1984 Pages: 1109-1121.

Musical Tone Noise

One problem with the spectral suppression methods is the possible presence of musical tone noise. In order to eliminate or mask the music noise, the suppression gain should be floored:G_n,k=max(σ,G_n,k)

Here σ is a constant which has the value between 0 and 1.

Noise reduction methods based on the above spectral gain have good performance for normal car noise. However when there is tonal noise at the background, these methods can only suppress the tonal noise but can not eliminate it. Referring now to FIG. 3, the PSD of a signal after prior art noise reduction is shown. The signal still has peaks at the frequencies where tonal noise is present. Thus, the overall signal is suppressed, but the tonal noise remains.

Adaptive Method

In order to remove tonal noise, instead of using a constant floor σ, the system uses a variable floor that is specified at each frequency bin. FIG. 9 is a flow-diagram illustrating the use of an adaptive factor to calculate a suppression gain value. At step 901 the smoothed background noise and the background noise estimate values are determined for a current frequency bin.

At step 902 the smoothed background value and background noise estimate value are used to generate a ratio. This ratio is used at step 903 to calculate the value for the adaptive factor to be used for the current bin. At step 904 the adaptive factor is used to generate the suppression gain value for the current bin. In this manner each frequency bin has a changing suppression gain floor that is dependent on the values of the ratio at that bin. The operation of the system of Figure is described as follows:

At a frequency bin estimate the background noise B_n(k) and calculate the smoothed background noise B_n(k) (step 901). The techniques above may be used to generate the values. At step 902 calculate the ratio ξ_n(k) as described above. This can then be used at step 903 to generate an adaptive factor σ that is related to the current frequency bin. The adaptive factor is defined by:σ_n,k=σ·ξ_n(k)

The tonal noise suppression gain to be applied to the signal (step 904) is then given by:Ĝ_n,k=max(σ_n,k,G_n,k)

Random Technique

Applying the above adaptive suppression gain to the spectral magnitude can achieve improved tonal noise removal. However, when there are severe tonal noises in the background, using the original noisy phase may make the tonal sound still audible in the processed signal. For, further smoothing, an alternate technique is to replace the original phases by random phases in the frequency bins whenever the adaptive suppression gain applied to the original noisy signal is less than the smoothed background noise.

FIG. 10 is a flow diagram illustrating a suppression technique using random phases. At step 1001 apply the adaptive gain suppression technique of FIG. 9. At step 1002 compare the result (multiplied by the noisy signal) to the smoothed background noise value for the current frequency bin. At decision block 1003 determine if the result is less than the smoothed background noise value. If no, the generated result can be used. If the result is less than the smoothed background noise then at step 1005 replace the original phase with a random phase. Step 1002 and 1003 can be implemented as follows:If Ĝ_n,k·|Y_n,k|< B_n(k)

The estimate of the clean speech spectral magnitude can be obtained (step 1001) as:|{circumflex over (X)}_n,k|=Ĝ_n,k·|Y_n,k|

The estimate of the complex clean speech is given by:{circumflex over (X)}_n,k=|{circumflex over (X)}_n,k|·(R_n,k+I_n,k·j)

Here R_n,k, I_n,k are two Gaussian random numbers with zero mean and unit variance.

FIG. 4 illustrates the PSD of the tonal noise processed by the disclosed tonal noise reduction method. As can be seen, the resulting waveform has fewer peaks and a more smooth profile.

The illustrations have been discussed with reference to functional blocks identified as modules and components that are, not intended to represent discrete structures and may be combined or further sub-divided. In addition, while various embodiments of the invention have been described, it will be apparent to those of ordinary skill in the art that other embodiments and implementations are possible that are within the scope of this invention. Accordingly, the invention is not restricted except in light of the attached claims and their equivalents.

Claims

1. A method of identifying tonal noise comprising: transforming an input signal into a plurality of frequency bins;at each bin calculating a smoothed background noise and a background noise estimate;at each bin comparing the smoothed background noise to the background noise estimate;calculating a ratio of the background noise estimate to the smoothed background noise for a bin;comparing the ratio to a predetermined threshold value;identifying whether a peak in the bin is a tonal peak or a non-tonal noise peak based on the comparison between the ratio and the predetermined threshold value;identifying the bin as having the tonal peak in response to a determination that the ratio of the background noise estimate to the smoothed background noise is greater than the predetermined threshold value; andattenuating at least a portion of the tonal peak of the input signal to generate an output signal with reduced tonal noise.

2. The method of claim 1 where the step of comparing comprises comparing the smoothed background noise to the background noise estimate of the same bin.

3. The method of claim 1 where the threshold value is greater than 1.

4. The method of claim 2 wherein the step of determining a smoothed background noise for a current frame n is accomplished by Bn(k)=β1*Bn(k)+(1−β1)*Bn(k−1)when Bn(k)≧ Bn(k−1)where Bn(k) is the background noise estimate of the present frame n at frequency bin k and Bn(k−1) is the smoothed background noise of the prior bin k−1.

5. The method of claim 4 wherein the step of determining a smoothed background noise is given by Bn(k)=β2*Bn(k)+(1−β2)*Bn(k−1)when Bn(k)< Bn(k−1).

6. The method of claim 1 wherein the ratio ξn(k) is given by ξn(k)=Bn(k)/Bn(k).

7. A method of removing tonal noise from a signal comprising: determining a short-time spectral magnitude |Yn,k| of a noisy speech signal at an nth frame and kth frequency bin;calculating a background noise estimate of the noisy speech signal at the kth frequency bin;calculating a smoothed background noise of the noisy speech signal at the kth frequency bin;calculating a ratio of the background noise estimate and the smoothed background noise;calculating an adaptive suppression gain value Ĝn,k based on the ratio of the background noise estimate and the smoothed background noise; andattenuating at least a portion of a tonal noise in the noisy speech signal to generate an estimated clean speech signal |{circumflex over (X)}n,k| by |{circumflex over (X)}n,k|=Ĝn,k|Yn,k|.

8. The method of claim 7 wherein Ĝn,k is generated by Ĝn,k=max(σn,k,Gn,k)where σn,k is an adaptive gain factor related to a current frequency bin.

9. The method of claim 8 where σn,k is generated by σn,k=σ·ξn(k)where σ is a constant factor andξn(k) is the ratio between the background noise estimate and the smoothed background noise at bin k.

10. The method of claim 9 where ξn(k)=Bn(k)/Bn(k)where Bn(k) is the background noise estimate of the current frame n at frequency k and Bn(k) is the smoothed background noise of the same bin.

11. The method of claim 10 further including the step of comparing Ĝn(k)·|Yn(k)| to Bn(k).

12. The-method of claim 11 further including the step of accepting |{circumflex over (X)}n,k| when Ĝn(k)·|Yn(k)≧ Bn(k).

13. The method of claim 11 further including the step of replacing the original phase with a random phase when Ĝn,k·|Yn,k|< Bn(k).

14. The method of claim 1 where the input signal comprises an audio signal with speech content and tonal noise content.

15. The method of claim 1 where the input signal comprises an audio signal with tonal noise content and diffuse noise content, and where the step of attenuating comprises attenuating the tonal peak associated with the tonal noise content by a greater amount than the diffuse noise content.

16. The method of claim 1 where the step of calculating the smoothed background noise comprises calculating the smoothed background noise by an asymmetric infinite impulse response filter.

17. The method of claim 5 where βhd 1 and β2 are two parameters in a range from 0 to 1, and where β2 is greater than β1.

18. A method of attenuating tonal noise comprising: determining a short-time spectral magnitude |Yn,k| of an audio input signal;transforming the input signal into a plurality of frequency bins;calculating a background noise estimate of the input signal at a first bin of the plurality of frequency bins;calculating a smoothed background noise of the input signal at the first bin;calculating a ratio of the background noise estimate and the smoothed background noise;comparing the ratio to a predetermined threshold value;identifying whether a peak in the first bin is a tonal noise peak or a non-tonal noise peak in response to the comparison between the ratio and the predetermined threshold value;identifying the first bin as having the tonal noise peak in response to a determination that the comparison meets a predetermined condition;calculating an adaptive suppression gain value Ĝn,k based on the ratio; andattenuating at least a portion of the tonal noise peak of the input signal to generate an audio output signal |{circumflex over (X)}n,k| with reduced tonal noise by |{circumflex over (X)}n,k|=Ĝn,k|Yn,k|.

19. The method of claim 7 wherein the step of calculating the adaptive suppression gain value Ĝn,k comprises changing a suppression gain floor associated with the adaptive suppression gain value Ĝn,k that is dependent on the ratio of the background noise estimate and the smoothed background noise.