EP2673778A1 - Post-processing including median filtering of noise suppression gains - Google Patents
Post-processing including median filtering of noise suppression gainsInfo
- Publication number
- EP2673778A1 EP2673778A1 EP12707413.6A EP12707413A EP2673778A1 EP 2673778 A1 EP2673778 A1 EP 2673778A1 EP 12707413 A EP12707413 A EP 12707413A EP 2673778 A1 EP2673778 A1 EP 2673778A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- post
- raw
- gains
- processing
- banded
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012805 post-processing Methods 0.000 title claims abstract description 83
- 238000001914 filtration Methods 0.000 title claims abstract description 67
- 230000001629 suppression Effects 0.000 title description 82
- 230000005236 sound signal Effects 0.000 claims abstract description 122
- 238000012545 processing Methods 0.000 claims abstract description 120
- 238000000034 method Methods 0.000 claims abstract description 101
- 238000009499 grossing Methods 0.000 claims description 43
- 230000003595 spectral effect Effects 0.000 claims description 35
- 230000000694 effects Effects 0.000 claims description 32
- 238000002592 echocardiography Methods 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 9
- 230000004907 flux Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 description 40
- 238000001228 spectrum Methods 0.000 description 27
- 230000006870 function Effects 0.000 description 16
- 230000014509 gene expression Effects 0.000 description 14
- 230000015654 memory Effects 0.000 description 14
- 230000001419 dependent effect Effects 0.000 description 11
- 238000001514 detection method Methods 0.000 description 9
- 230000005534 acoustic noise Effects 0.000 description 7
- 238000013459 approach Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 230000035945 sensitivity Effects 0.000 description 6
- 238000003786 synthesis reaction Methods 0.000 description 6
- 238000004590 computer program Methods 0.000 description 5
- 230000009471 action Effects 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 230000008447 perception Effects 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 230000001131 transforming effect Effects 0.000 description 4
- 230000006978 adaptation Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000010348 incorporation Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000003638 chemical reducing agent Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 238000013479 data entry Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 239000012925 reference material Substances 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
- G10L21/0216—Noise filtering characterised by the method used for estimating noise
- G10L2021/02161—Number of inputs available containing the signal or the noise to be suppressed
- G10L2021/02166—Microphone arrays; Beamforming
Definitions
- the present disclosure relates generally to signal processing, in particular of audio signals.
- Acoustic signal processing is applicable today to improve the quality of sound signals such as from microphones.
- many devices such as handsets operate in the presence of sources of echoes, e.g., loudspeakers.
- signals from microphones may occur in a noisy environment, e.g., in a car or in the presence of other noise.
- there may be sounds from interfering locations e.g., out-of-location conversation by others, or out-of-location interference, wind, etc. Acoustic signal processing is therefore an important area for invention.
- An acoustic noise reduction system typically includes a noise estimator and a gain calculation module to determine suppression probability indicators, e.g., as a set of noise reduction gains that are determined, for example, on a set of frequency bands, and applied to the (noisy) input audio signal after transformation to the frequency domain and banding to the set of frequency bands to attenuate noise components.
- the acoustic noise reduction system may include one microphone input, or a plurality of microphone inputs and downmixing, e.g. , beamforming to generate one input audio signal.
- the acoustic noise reduction system may further include echo reduction, and may further include out-of-location signal reduction.
- the gains determined for each band for leveling and/or dynamic equalization include statistical outliers, e.g., isolated values, and such outliers might cause artifacts such as musical noise.
- Median filtering the gains, e.g., noise reduction gains, or leveling and/or dynamic equalization gains across frequency bands can reduce musical noise artifacts.
- Gain values may vary significantly across frequencies, and in such a situation, running a relatively wide median filter along frequency bands has the risk of disrupting the continuity of temporal envelope, which is the inherent property for many signals and is crucial to perception as well. Whilst offering greater immunity to the outliers, a longer median filter can reduce the spectral selectivity of the processing, and potentially introduce greater discontinuities or jumps in the gain values across frequency and time.
- FIG. 1 shows one example of processing of a set of one or more input audio
- microphone signals e.g., microphone signals 101 from differently located microphones, including an embodiment of the present invention.
- FIG. 2 shows diagrammatically sets of raw banded gains and the frequency
- FIG. 3A shows a simplified block diagram of a post-processor that includes a median filter according to an embodiment of the present invention.
- FIG. 3B shows a simplified flowchart of a method of post-processing that includes median filtering according to an embodiment of the present invention.
- FIG. 4 shows one example of an apparatus embodiment configured to determine a set of post-processed gains for suppression of noise, and in some versions, simultaneous echo suppression, and in some versions, simultaneous suppression of out-of-location signals.
- FIG. 5 shows one example of an apparatus embodiment in more detail.
- FIG. 6 shows an example embodiment of a gain calculation element that includes a spatially sensitive voice activity detector and a wind activity detector.
- FIG. 7 shows a flowchart of an embodiment of a method of operating a processing apparatus to suppress noise and out-of-location signals and, in some embodiments, echoes.
- FIG. 8 shows a simplified block diagram of a processing apparatus embodiment for processing one or more audio inputs to determine a set of raw gains, to post-process the raw gains including median filtering the determined raw gains, and to generate audio output that has been modified by application of the post-processed gains.
- Embodiments of the present invention include a method, an apparatus, and logic encoded in one or more computer-readable tangible media to carry out the method.
- One embodiment includes a method of applying post-processing to raw banded gains to improve the raw banded gains for applying to one or more input audio signals.
- the raw banded gains at a plurality of frequency bands comprising one or more frequency bins are determined by input processing the one or more input audio signals.
- the raw banded gains are to carry out one or more of reducing noise, reducing out-of-location signals, reducing echoes, perceptual domain-based leveling, perceptual domain-based dynamic range control, and perceptual domain-based dynamic equalization.
- the method comprises applying post-processing to the raw banded gains to generate banded post- processed gains.
- Generating of a particular post-processed gain for a particular frequency band includes at least median filtering using gain values for frequency bands adjacent to the particular frequency band.
- the post-processing is according to one or more properties, including an end condition and a width for the median filtering. At least one property of the post-processing depends on signal classification of the one or more input audio signals.
- One embodiment includes a method of processing one or more input audio signals.
- the method includes input processing one or more input audio signals to determine raw banded gains for applying to an audio signal, the raw banded gains being at a plurality of frequency bands comprising one or more frequency bins.
- the raw banded gains are to carry out one or more of reducing noise, reducing out-of-location signals, reducing echoes, perceptual domain-based leveling, perceptual domain-based dynamic range control, and perceptual domain-based dynamic equalization.
- the method further comprises applying post-processing to the raw banded gains to generate banded post- processed gains.
- Generating of a particular post-processed gain for a particular frequency band includes at least median filtering using gain values for frequency bands adjacent to the particular frequency band.
- the post-processing is according to one or more properties, including an end condition and a width for the median filtering. At least one property of the post-processing depends on signal classification of the one or more input audio signals.
- One embodiment includes an apparatus to post-process raw banded gains for
- the raw banded gains at a plurality of frequency bands comprising one or more frequency bins are determined by input processing the one or more input audio signals.
- the raw banded gains are to carry out one or more of reducing noise, reducing out-of-location signals, reducing echoes, perceptual domain-based leveling, perceptual domain-based dynamic range control, and perceptual domain-based dynamic equalization.
- the apparatus comprises a post-processor accepting the raw banded gains and applying post-processing to the raw banded gains to generate banded post-processed gains to apply to the one or more input signals.
- the post-processor includes a median filter to carry out median filtering of the raw banded gains.
- Generating by the post-processor of a particular post-processed gain for a particular frequency band includes the median filtering using gain values for frequency bands adjacent to the particular frequency band.
- the post-processing is according to one or more properties, including an end condition and a width for the median filtering. At least one property of the post-processing depends on signal classification of the one or more input audio signals.
- One embodiment includes an apparatus to process one or more input audio signals.
- the apparatus comprises an input processor accepting the one or more input audio signals and input processing the one or more input audio signals to generate raw banded gains at a plurality of frequency bands comprising one or more frequency bins.
- the raw banded gains are to carry out one or more of reducing noise, reducing out-of-location signals, reducing echoes, perceptual domain-based leveling, perceptual domain-based dynamic range control, and perceptual domain-based dynamic equalization.
- the apparatus further comprises a post-processor accepting the raw banded gains and applying post-processing to the raw banded gains to generate banded post-processed gains to apply to the one or more input signals.
- the post-processor includes a median filter to carry out median filtering of the raw banded gains.
- Generating by the post-processor of a particular post- processed gain for a particular frequency band includes the median filtering using gain values for frequency bands adjacent to the particular frequency band.
- the post-processing is according to one or more properties, including an end condition and a width for the median filtering. At least one property of the post-processing depends on signal classification of the one or more input audio signals.
- One embodiment includes a system for post-processing raw banded gains to
- the system comprises means for post-processing raw banded gains to generate banded post-processed gains, the raw banded gains determined by a means for input processing one or more input audio signals to generate the raw banded gains at a plurality of frequency bands comprising one or more frequency bins.
- the raw banded gains are to carry out one or more of reducing noise, reducing out-of-location signals, reducing echoes, perceptual domain-based leveling, perceptual domain-based dynamic range control, and perceptual domain-based dynamic equalization.
- the banded post-processed gains are for applying to the one or more audio signals.
- Generating a particular post-processed gain for a particular frequency band includes at least median filtering using gain values for frequency bands adjacent to the particular frequency band.
- the post-processing is according to one or more properties, including an end condition and a width for the median filtering. At least one property of the post-processing depends on signal classification of the one or more input audio signals.
- One embodiment includes a system for processing one or more input audio
- the system comprises means for input processing the one or more input audio signals to generate raw banded gains at a plurality of frequency bands comprising one or more frequency bins.
- the raw banded gains are to carry out one or more of reducing noise, reducing out-of-location signals, reducing echoes, perceptual domain-based leveling, perceptual domain-based dynamic range control, and perceptual domain-based dynamic equalization.
- the system further comprises means for post-processing raw banded gains to generate banded post-processed gains to apply to the one or more input audio signals to carry out the one or more of reducing noise, reducing out-of-location signals, reducing echoes, perceptual domain-based leveling, perceptual domain-based dynamic range control, and perceptual domain-based dynamic equalization.
- the generating by the means for post-processing of a particular post-processed gain for a particular frequency band includes at least median filtering using gain values for frequency bands adjacent to the particular frequency band.
- the post-processing applied by the means for post-processing is according to one or more properties, including an end condition and a width for the median filtering. At least one property of the post-processing depends on signal classification of the one or more input audio signals.
- the post-processing includes at least one of frequency- band-to-frequency-band smoothing and smoothing across time.
- one of or both the width and the end conditions of the median filtering depend on signal classification of the one or more input audio signals.
- the classification includes whether the input audio signals are likely or not to be voice. In some embodiments, the classification includes whether the input audio signals are likely or not to be wind.
- the frequency bands are on a perceptual or logarithmic scale.
- the raw banded gains determined from one or more input audio signals are for reducing noise.
- the raw banded gains are determined from more than one input audio signal and are for reducing noise and out-of- location signals.
- the raw banded gains are determined from one or more input audio signals and one or more reference signals, and are for reducing noise and echoes.
- One embodiment includes a tangible computer-readable storage medium
- One embodiment includes program logic that when executed by at least one
- processor causes carrying out a method as described herein.
- Particular embodiments may provide all, some, or none of these aspects, features, or advantages. Particular embodiments may provide one or more other aspects, features, or advantages, one or more of which may be readily apparent to a person skilled in the art from the figures, descriptions, and claims herein.
- One aspect of the invention is processing of one or more input audio signals, including input processing to generate raw gains for noise reduction, or for other forms of input signal improvement.
- the processing includes applying post-processing to the raw gains, including median filtering of raw gains for gain smoothing.
- Another aspect of the invention is the post-processing that includes the median filtering of raw gains determined by input processing, e.g., for noise reduction or for other input processing.
- a median filter replaces a particular raw gain value with the median of a predefined number of raw gain values, e.g., by the median of the particular raw gain value and a predefined set of neighboring raw gain values.
- a median filter has one or more properties, e.g., the number of valued of which the median is determined, and the end conditions.
- FIG. 1 shows one example of processing of a set of one or more input audio signals, e.g., microphone signals 101 from differently located microphones, including an embodiment of the present invention.
- the processing is by time frames of a number, e.g., M samples.
- there is only one input e.g., one microphone
- there is a plurality, denoted P of inputs e.g., microphone signals 101.
- An input processor 105 accepts sampled input audio signal(s) 101 and forms a banded instantaneous frequency domain amplitude metric 119 of the input audio signal(s) 101 for a plurality B of frequency bands. In some embodiments in which there is more than one input audio signal, the metric 119 is mixed-down from the input audio signal.
- the amplitude metric represents the spectral content. In many of the embodiments described herein, the spectral content is in terms of the power spectrum. However, the invention is not limited to processing power spectral values. Rather, any spectral amplitude dependent metric can be used. For example, if the amplitude spectrum is used directly, such spectral content is sometimes referred to as spectral envelope. Thus, the phrase "power (or other amplitude metric) spectrum" is sometimes used in the description.
- the input processor 105 determines a set of raw banded gains 111 to apply to the instantaneous amplitude metric 119.
- the input processing further includes determining a signal classification of the input audio signal(s), e.g., an indication of whether the input audio signal(s) is/are likely to be voice or not as determined by a voice activity detector (VAD), and/or an indication of whether the input audio signal(s) is/are likely to be wind or not as determined by a wind activity detector (WAD), and/or an indication that the signal energy is rapidly changing as indicated, e.g., by the spectral flux exceeding a threshold.
- VAD voice activity detector
- WAD wind activity detector
- a feature of embodiments of the present invention includes applying post- processing to the raw gains to improve the quality of the output.
- the post-processing includes median filtering of the raw gains determined by the input processing.
- a median filter considers a set of raw gains and outputs the gain that is the median of the set of raw gains.
- a set of B raw gains is determined every frame, so that there is a time sequence of sets of B raw gains over B frequency bands.
- the median filter extends across frequency.
- FIG. 2 shows diagrammatically sets of raw banded gains, one set for each of the present time, one frame back, two frames back, three frames back, etc., and further shows the coverage of an example median filter that includes five raw gain values centered around a frequency band b c in the present frame.
- filter width we mean the width of the filter in the frequency band domain.
- the post-processing produces a set of post-processed gains 125 that are applied to the instantaneous power (or other amplitude metric) 119 to produce output, e.g., as a plurality of processed frequency bins 133.
- An output synthesis filterbank 135 (or for subsequent coding, a transformer/remapper) converts these frequency bins to desired output 137.
- Input processing element 105 includes an input analysis filterbank, and a raw gain calculator.
- the input analysis filterbank for the case of one input audio signal 101, includes a transformer to transform the samples of a frame into frequency bins, and a banding element to form frequency bands, most of which include a plurality of frequency bins.
- the input analysis filterbank for the case of a plurality of input audio signals 101, includes a transformer to transform the samples of a frame of each of the input audio signals into frequency bins, a downmixer, e.g., a beamformer to downmix the plurality into a single signal, and a banding element to form frequency bands, most of which include a plurality of frequency bins.
- the transformer implements short time Fourier transform
- the transformer uses a discrete finite length Fourier transform (DFT) implemented by a fast Fourier transform (FFT).
- DFT discrete finite length Fourier transform
- FFT fast Fourier transform
- Other embodiments use different transforms.
- the B bands are at frequencies whose spacing is
- the bands are monotonically non-decreasing.
- a reasonable number, e.g., 90% of the frequency bands include contribution from more than one frequency bin, and in particular embodiments, each frequency band includes contribution from two or more frequency bins.
- the bands are monotonically increasing in a log-like manner.
- the bands are on a psycho-acoustic scale, that is, the frequency bands are spaced with a scaling related to psycho-acoustic critical spacing, such banding called "perceptually-banding" herein.
- the band spacing is around 1 ERB or 0.5 Bark, or equivalent bands with frequency separation at around 10% of the center frequency.
- a reasonable range of frequency spacing is from 5-20% or
- the input processing includes noise reduction
- the input processing a plurality of input audio signals are accepted by the input processor, and the input processing includes reducing out-of location signals.
- An example of input processing that includes reducing out-of location signals is described in concurrently filed
- the input processing also includes echo reduction.
- One example of input processing that includes echo reduction is described in concurrently filed International Application No.
- the input processing includes echo reduction
- one or more reference signals also are included and used to obtain an estimate of some property of the echo, e.g., of the power (or other amplitude metric) spectrum of the echo.
- the resulting raw banded gains achieve simultaneous echo reduction and noise reduction.
- the post- processed gains are accepted by an element 123 that modifies the gains to include additional echo suppression.
- the result is a set of post-processed gains 125 that are applied to the input signal or signals, e.g., that are used to process, in the frequency domain, .as frequency bins, the input audio signal if there is one input, or a downmix of the input audio signals if there is a plurality of input audio signals, e.g., from differently located microphones.
- Gain application module 131 accepts the banded post-processed gains 125 and applies such gains to the input audio signal or signals.
- the processed data 133 may then be converted back to the sample domain by an output synthesis filterbank 135 to produce a frame of M signal samples 137.
- the signal 133 is subject to transformation or remapping, e.g., to a form ready for coding according to some coding method.
- the invention is not limited to the input processing and gain calculation described in International Application No. PCT/US , U.S . 61/441 ,611 , or even to noise reduction.
- the input processing is to reduce noise (and possibly echo and out of location signals)
- the input processing may be, additionally or primarily, to carry out one or more of perceptual domain-based leveling, perceptual domain-based dynamic range control, and perceptual domain-based dynamic equalization that take into account the variation in the perception of audio depending on the reproduction level of the audio signal, as described, for example, in commonly owned WO 2004111994.
- the raw banded gains calculated per WO 2004111994 are post- processed, including median filtering, to determine post-processed gains 125 to apply to the (transformed) input.
- FIG. 3 A shows a simplified block diagram of a post-processor 121 that includes a median filter 305 according to an embodiment of the present invention.
- the postprocessor 121 accepts raw gains 111 and in embodiments in which the post-processing changes according to signal classification, one or more signal classification indicators 115, e.g., the outputs of one or more of a VAD and a WAD. While not included in all embodiments, some embodiments of the post-processor include a minimum gain processor 303 to ensure that the gains do not fall below a predefined, possibly frequency- dependent value.
- some embodiments of the post-processor include a smoothing filter 307 that processes the gains after median filtering to smooth frequency-band-to-frequency-band variations, and/or to smooth time variations.
- FIG. 3B shows a simplified flowchart of a method of post-processing 310 that includes in 311 accepting raw gains, and in embodiments in which the post-processing changes according to signal classification, one or more signal classification indicators 115.
- the post-processing includes median filtering 315 according to embodiments of the present invention. The inventors have found that median filtering is a powerful nonlinear smoothing technique, which works well for eliminating undesired outliers when compared with only using a smoothing method.
- Some embodiments include in step 313 ensuring that the gains do not fall below a predefined minimum, which may be frequency band dependent. Some embodiments further include, in step 317, band-to-band and/or time smoothing, e.g., linear smoothing using, e.g., a weighted moving average.
- band-to-band and/or time smoothing e.g., linear smoothing using, e.g., a weighted moving average.
- a median filter 315 of raw banded gain values is characterized by: 1) the number of raw banded gains to include to determine the median value, 2) the frequency band positions of the raw banded gains that are included; 3) the edge conditions, i.e., the conditions used to extend the raw banded gains to allow calculation of the median at the edges of time and frequency band; and 4) how the characterization of the median filter is affected by the signal classification, e.g., one or more of the presence of voice and the presence of wind.
- Some embodiments include a mechanism to control one or more of the median filtering characteristics over frequency and/or time based on signal classification. For example, in one embodiment that includes voice activity detection, one or more of the median filtering characteristics vary in accordance to whether the input is ascertained by a VAD to be voice or not. In one embodiment that includes wind activity detection, one or more of the median filtering characteristics vary in accordance to whether the input is ascertained by a WAD to be wind or not.
- Examples of different edge conditions include (a) extrapolating of interior values for the edges; (b) using the minimum gain value to extend the raw banded gains at the edges, (c) using a zero gain value to extend the raw banded gains at the edges (d) duplicating the central filter position value to extend the raw banded gains at the edges, and (e) using a maximum gain value to extend the raw banded gains at the edges. Additional post-processing
- the post-processor 121 includes a minimum gain processor 303 that carries out step 313 to ensure the gains do not fall below a predefined minimum gain value.
- the minimum gain processor ensures minimum values in a frequency-band dependent manner.
- the manner of prevention minimum is dependent on the activity classification 115, e.g., whether voice or not.
- Gain b ' RAW Gain b ' MIN + (l - Gain b ' MIN ) ⁇ Gain b ' s
- Gain' b RAW Gain' b M[N +Gain' b s
- Gain ⁇ Gainb' m Gain ⁇ ' ⁇ Gain l,MiN
- the range of the maximum suppression depth or minimum gain may range from -80dB to
- Gain ⁇ MjN is increased, e.g., in a frequency-band dependent way (or in another embodiment, by the same amount for each band b). In one embodiment, the amount of increase in the minimum is larger in the mid- frequency bands, e.g., bands between 500 Hz to 2 kHz.
- the post-processor 121 includes a smoothing filter 307, e.g., a linear smoothing filter that carries out one or both of frequency band-to-band smoothing and time smoothing. In some embodiments, such smoothing is varied according to signal classification 115.
- a smoothing filter 307 e.g., a linear smoothing filter that carries out one or both of frequency band-to-band smoothing and time smoothing. In some embodiments, such smoothing is varied according to signal classification 115.
- smoothing 317 uses a weighted moving average with a fixed kernel.
- One example uses a binomial approximation of a Gaussian weighting kernel for the weighted moving average.
- a 5 -point binomial smoother has a kernel
- a 3-point binomial smoother has a kernel - ⁇ - [l 2 l] .
- Many other weighted moving average filters are known, and any such filter can suitably be modified to be used for the band-to-band smoothing of the gain.
- the band-to-band smoothing is controlled by the signal
- a VAD e.g., a spatially-selective VAD is included, and if the VAD determines there is voice, the degree of smoothing is increased when noise is detected.
- 5 -point band-to-band weighted average smoothing is carried out in the case the VAD indicates voice is detected, else, when the VAD determines there is no voice, no smoothing is carried out.
- time smoothing of the gains also is included. In some embodiments, time smoothing of the gains also is included.
- the gain of each of the B bands is smoothed by a first order smoothing filter:
- Gain b ⁇ Smoothed a b Gain b + (l - a b )Gain b Smoothed ⁇ gv
- Gain b is the current time-frame gain
- G ain b, Smoothed ⁇ ev is Gain b,Smoothed from the previous M-sample frame.
- a b is a time constant which may be frequency band dependent and is typically in the range of 20 to 500ms. In one embodiment a value of 50ms was used.
- the amount of time smoothing is controlled by the signal classification of the current frame.
- the signal classification of the current frame is used to control the values of first order time constants used to filter the gains over time in each band.
- the parameters of post-processing are controlled by the immediate signal classifier (VAD, WAD) value that has low latency and is able to achieve a rapid transition of the post-processing from noise into voice (or other desired signal) mode.
- VAD immediate signal classifier
- WAD voice-to-live
- the speed with which more aggressive post-processing is reinstated after detection of voice, i.e., at the trail out, has been found to be less important, as it affects intelligibility of voice to a lesser extent.
- the band-to-band median filtering is controlled by the signal classification.
- a VAD is included, and if the VAD determines it is likely that there is no voice, a 7 point T-shaped median filter with 5-point band-to-band and 3 -point time median filtering is carried out, with edge processing including extending minimum gain values or a zero value at the edges to compute the median value.
- a 5-point T-shaped time-frequency median filtering is carried out with three frequency bands in the current time frame, and using two previous time frames, and in a second embodiment, a three point memory less frequency-band only median filter, with the edge values extrapolated at the edges to calculate the median, is used.
- the median value is the median value, such that the median filter is a median filter.
- the post-processing e.g., the median filtering depends on the classification of the signal, and one such classification, in some embodiments, is whether there is wind or not.
- a WAD is included, and if the WAD determines there is no wind, and a VAD indicates there is no voice, fewer raw gain values are included in the median filter.
- the median filtering should be shorter, e.g., by using 3-point band-to-band median filter, with extrapolating the edge values applied at the edges.
- more median filtering can be used, e.g., 5-point band-to-band median filtering is carried out, with edge processing including extending minimum gain values or a zero value at the edges to compute the median value. If the WAD indicated wind is likely, and the VAD indicates voice is unlikely, even more median filtering can be used, e.g., a 7-point band-to-band median filtering can be carried out, with edge processing including extending minimum gain values or a zero value at the edges to compute the median value.
- An acoustic noise reduction system typically includes a noise estimator and a raw gain calculation module to determine a set of noise reduction gains that are determined, for example, on a set of frequency bands, and applied to the (noisy) input audio signal after transformation to the frequency domain and banding to the set of frequency bands to attenuate noise components.
- the acoustic noise reduction system may include one microphone, or a plurality of inputs from differently located microphones and
- the acoustic noise reduction system may further include echo reduction, and may further include out-of- location signal reduction.
- FIG. 4 shows one example of an apparatus configured to determine a set of post- processed gains for suppression of noise, and in some versions, simultaneous echo suppression, and in some versions, simultaneous suppression of out-of-location signals.
- the inputs include a set of one or more input audio signals 101 , e.g., signals from differently located microphones, each in sets of M samples per frame.
- input audio signals 101 e.g., signals from differently located microphones, each in sets of M samples per frame.
- there are two or more input audio signals e.g., signals from spatially separated microphones.
- one or more reference signals 103 are also accepted, e.g., in frames of M samples.
- a first input processing stage 403 determines a banded signal power (or other amplitude metric) spectrum 413 denoted P , and a banded measure of the instantaneous power 417 denoted Y .
- P a banded signal power
- Y a banded measure of the instantaneous power 417
- each of the spectrum 413 and instantaneous banded measure 417 is of the inputs after being mixed down by a downmixer, e.g., a beamformer.
- the first input processing stage 403 When echo suppression is included, the first input processing stage 403 also determines a banded power spectrum estimate of the echo 415, denoted E , the determining being from a previously calculated power spectrum estimates of the echo using a filter with a set of adaptively determined filter coefficients. In those versions that include out-of-location signal suppression, the first input processing stage 403 also determines spatial features 419 in the form of banded location probability indicators 419 that are usable to spatially separate a signal into the components originating from the desired location and those not from the desired direction.
- Embodiments of the second stage 405 include a noise power (or other amplitude metric) spectrum calculator 421 to determine a measure of the noise power (or other amplitude metric) spectrum, denoted E , and a signal classifier 423 to determine a signal classification 115, e.g., one or more of a voice activity detector (VAD), a wind activity detector, and a power flux calculator.
- VAD voice activity detector
- FIG. 4 shows the signal classifier 423 including a VAD.
- FIG. 5 shows one embodiment 500 of the elements of FIG. 4 in more detail, and includes, for the example embodiment of noise, echo, and out-of-location noise suppression, the suppressor 131 that applied the post-processed gains 125 and the output synthesizer (or transformer or remapper) 135 to generate the output signal 137.
- the first stage processor 403 of FIG. 4 includes
- the input(s) frame(s) 101 are transformed by inputs transformer(s) 503 to determine transformed input signal bins, the number of frequency bins denoted by N.
- the frequency domain signals from the input transformers 503 are accepted by a banded spatial feature calculator to determine banded location probability indictors, each between 0 and 1.
- the signals are combines by combiner 511, in one embodiment a summer, to produce a combined reference input.
- Y' j is used as a good-enough approximation to P .
- the L B filter coefficients for filter 517 are determined by an adaptive filter
- the updating is triggered by a voice activity signal denoted S as determined by a voice activity detector (VAD) 525 using P (or N'3 ⁇ 4, and E .
- S a voice activity detector
- the VAD derived in the echo update voice- activity detector 525 and filter updater 527 serves the specific purpose of controlling the adaptation of the echo prediction.
- a VAD or detector with this purpose is often referred to as a double talk detector.
- the echo filter coefficient updating of updater 527 is gated, with updating occurring when the expected echo is significant compared to the expected noise and current input power, as determined by the VAD 525 and indicated by a low value of local signal activity S.
- the input transformers 503, 511 determine the short time Fourier transform (STFT).
- the following transform and inverse pair is used for the forward transform in elements 503 and 511 , and in output synthesis element 135.
- W Center' o
- n and v n are appropriate window functions
- x n represents the last 2N input samples with x N-l representing the most recent sample
- X n represents the N complex- valued frequency bins in increasing frequency order.
- the inverse transform or synthesis is represented in the last two equation lines.
- y n represents the 2N output samples that result from the individual inverse transform prior to overlapping, adding and discarding as appropriate for the designed windows. It should be noted, that this transform has an efficient implementation as a block multiply and FFT. Note that the use of x n and X n in the above expressions of transform is for convenience.
- the window functions u n and v n for the above transform in one embodiment is the sinusoidal window family, of which one suggested embodiment is
- windows also known as prototype filters, can be of length greater or smaller than the examples given herein.
- the downmixer is a beamformer 507 designed to achieve some spatial selectivity towards the desired position.
- the beamformer 507 is a linear time invariant process, i.e., a passive beamformer defined in general by a set of complex-valued frequency-dependent gains for each input channel. For the example of a two-microphone array, with the desired sound source located broad side to the array, i.e., at the
- one embodiment uses for beamformer 507 a passive beamformer 107 that determines the simple sum of the two input channels.
- beamformer 507 weights the sets of inputs (as frequency bins) by a set of complex valued weights.
- the beamforming weights of beamformer 107 are determined according to maximum-ratio combining (MRC).
- MRC maximum-ratio combining
- the beamformer 507 uses weights determined using zero-forcing. Such methods are well known in the art.
- n 0 where 3 ⁇ 4is the banded instantaneous power of the mixed-down, e.g., beamformed signal, W b is the normalization gain and w b n are elements from a banding matrix.
- W b is the normalization gain
- w b n are elements from a banding matrix.
- 3 ⁇ 4' PREV is a previously, e.g., the most recently determined signal power (or other frequency domain amplitude metric) estimate
- ap 3 ⁇ 4 is a time signal estimate time constant
- Y ⁇ n is an offset.
- a suitable range for the signal estimate time constant ap h was found to be between 20 to 200 ms.
- the offset Y ⁇ n is added to avoid a zero level power spectrum (or other amplitude metric spectrum) estimate, ⁇ min can b e measured, or can be selected based on a priori knowledge.
- Y ⁇ n for example, can be related to the threshold of hearing or the device noise threshold.
- the adaptive filter 517 includes determining the instantaneous echo power spectrum (or other amplitude metric spectrum), denoted T b for band b by using an L tap adaptive filter described by L-l
- n ⁇ F b X b ' ,
- One embodiment includes time smoothing of the instantaneous echo from echo prediction filter 517 to determine the echo spectral estimate E b ' .
- a first order time smoothing filter is used as follows
- E b ' T b for T b ' > E b ' pi ev
- Eb' a E ,b T b + I 1 - ⁇ *E,b tepr ev for T b ⁇ E b' Pr ev
- ctg b is a first order smoothing time constant
- the noise power spectrum calculator 523 uses a minimum follower with exponential growth:
- b is a parameter that specifies the rate over time at which the minimum follower can increase to track any increase in the noise.
- the parameter b is best expressed in terms of the rate over time at which minimum follower will track. That rate can be expressed in dB/sec, which then provides a mechanism for determining the value of b .
- the range is 1 to 30dB/sec. In one embodiment, a value of 20dB/sec is used.
- Examples of such different approaches include but are not limited to alternate methods of determining a minimum over a window of signal observation, e.g., a window of 1 and 10 seconds. In addition or alternate to the minimum, such different approaches might also determine the mean and variance of the signal during times that it is classified as likely to be noise or that voice is unlikely.
- the one or more leak rate parameters of the minimum follower are controlled by the probability of voice being present as determined by voice activity detecting (VAD).
- VAD element 525 determines an overall signal activity level denoted S as
- ⁇ ⁇ , ⁇ ⁇ > ⁇ are margins for noise end echo, respectively and Y s ' ens is a settable sensitivity offset. These parameters may in general vary across the bands. In one embodiment, the values of ⁇ ⁇ , ⁇ ⁇ are between 1 and 4. In a particular embodiment, ⁇ ⁇ , ⁇ ⁇ are each 2. Y' sens is set to be around expected microphone and system noise level, obtained by experiments on typical components. Alternatively, one can use the threshold of hearing to determine a value for Y sens . [0097] In one embodiment, the echo filter coefficient updating of updater 527 is gated, as follows.
- the adaptive filter coefficients are updated as: mnoo, , , (max(0, 3 ⁇ 4 - y N N b ' )-n)x b '
- F b F b + ⁇ 2 if 5 ⁇ S thresh ,
- ⁇ ⁇ is a tuning parameter tuned to ensure stability between the noise and echo estimate.
- a typical value for ⁇ ⁇ is 1.4 (+3dB).
- a range of values 1 to 4 can be used, ⁇ is a tuning parameter that affects the rate of convergence and stability of the echo estimate. Values between 0 and 1 might be useful in different embodiments.
- ⁇ is a tuning parameter tuned to ensure stability between the noise and echo estimate.
- X s ' ens is set to avoid unstable adaptation for small reference signals. In one embodiment X s ' ens is related to the threshold of hearing.
- the choice of value for S thresh depends on the number of bands. S thresh is between 1 and B, and for one embodiment having 24 bands to 8kHz, a suitable range was found to be between 2 and 8, with a particular embodiment using a value of 4.
- Embodiments of the present invention use spatial information in the form of one or more measures determined from one or more spatial features in a band b that are monotonic with the probability that the particular band b has such energy incident from a spatial region of interest. Such quantities are called spatial probability indicators.
- the one or more spatial probability indicators are functions of one or more banded weighted covariance matrices of the input audio signals.
- R' b ⁇ W b Xl,n ⁇ ⁇ ⁇ 3 ⁇ 4 r l3 ⁇ 4 ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇
- the 1 ⁇ 43 ⁇ 4 n provide an indication of how each bin is weighted for contribution to the bands.
- the one or more covariance matrices are smoothed over time.
- the banding matrix includes time dependent weighting for a weighted moving average, denoted as W3 ⁇ 4 / with elements 1 ⁇ 43 ⁇ 4 n ; , where / represents the time frame, so that, over L time frames, N-lL-l
- R' b ⁇ w b ⁇ [x ... Xp H [X ... Xp .
- each band covariance matrix R'3 ⁇ 4 is a 2x2 Hermetian positive definite matrix
- Rb' 2i Rb' ⁇ 2 > where the overbar is used to indicate the complex conjugate.
- ratio a quantity that is monotonic with the ratio of
- ⁇ is a small offset added to avoid singularities
- ⁇ can be thought of as the smallest expected value for R ⁇ j .
- it is the determined, or estimated (a priori) value of the noise power (or other frequency domain amplitude metric) in band b for the microphone and related electronics. That is, the minimum sensitivity of any preprocessing used.
- the coherence feature is
- One feature of some embodiments of the noise, echo and out-of-location signal suppression is that, based on the a priori expected or current estimate of the desired signal features—the target values, e.g., representing spatial location, gathered from statistical data— each spatial feature in each band can be used to create a probability indicator for the feature for the band b.
- desired location are modeled as Gaussian distributions that present a robust way of capturing the region of interest for probability indicators derived from each spatial feature and band.
- the function f Rb (ARatio') is a smooth function.
- the ratio probability indicator function is
- Width Rati 0 b is a width tuning parameter expressed in log units, e.g., dB.
- the Width Ratio £ is related to but does not need to be determined from actual data. It is set to cover the expected variation of the spatial feature in normal and noisy conditions, but also needs only be as narrow as is required in the context of the overall system to achieve the desired suppression.
- the function fp b ⁇ APhase') is a smooth function.
- f Rh ⁇ APhase b )
- Widthp ⁇ b is a width tuning parameter expressed in units of phase.
- Width ph ase 3 ⁇ 4 is related to but does not need to be determined from actual data.
- no target is used, and in one embodiment,
- CF actor ⁇ is a tuning parameter that may be a constant value in the range of 0.1 to 10; in one embodiment, a value of 0.25 was found to be effective.
- FIG. 6 shows one example of the calculation in element 529 of the raw gains, and includes a spatially sensitive voice activity detector (VAD) 621, and a wind activity detector (WAD) 623.
- VAD voice activity detector
- WAD wind activity detector
- Alternate versions of noise reduction may not include the WAD, or the spatially sensitive VAD, and further may not include echo suppression or other reduction.
- the embodiment shown in FIG. 6 includes additional echo suppression, which may not be included in simpler versions.
- the spatial probability indicators are used to determine what is referred to as the beam gain, a statistical quantity denoted BeamGain 3 ⁇ 4 that can be used to estimate the in-beam and out-of-beam power from the total power, e.g., using an out-of- beam spectrum calculator 603, and further, can be used to determine the out-of-beam suppression gain by a spatial suppression gain calculator 611.
- BeamGain 3 ⁇ 4 a statistical quantity denoted BeamGain 3 ⁇ 4 that can be used to estimate the in-beam and out-of-beam power from the total power, e.g., using an out-of- beam spectrum calculator 603, and further, can be used to determine the out-of-beam suppression gain by a spatial suppression gain calculator 611.
- the probability indicators are scaled such that the beam gain has a maximum value of 1.
- the beam gain is
- BeamGain BeamGain min + ( 1- BeamGain min )RPI -PP1 ' "CPl' j ,.
- Some embodiments use BeamGain ⁇ of 0.01 to 0.3 (-40dB to -lOdB).
- One embodiment uses a BeamGain m ⁇ utilizat of 0.1.
- Power' b j n B eam an ⁇ ⁇ Power' b,OutOfBeam are statistical measures used for suppression.
- One version of gain calculation uses a spatially-selective noise power spectrum calculator 605 that determines an estimate of the noise power (or other metric of the amplitude) spectrum.
- One embodiment of the invention uses a leaky minimum follower, with a tracking rate determined by at least one leak rate parameter.
- the leak rate parameter need not be the same as for the non-spatially-selective noise estimation used in the echo coefficient updating.
- N' b $ the spatially-selective noise spectrum estimate.
- N b ' t S , (l + a b )Nfc, Spr ev )
- N b sp T ev i the already determined, i.e., previous value of N' b ⁇ .
- the leak rate parameter a b is expressed in dB/s such that for a frame time denoted T, (l + a b )VT is between 1.2 and 4 if the probability of voice is low, and 1 if the probability of voice is high.
- echo gating is used, i.e.,
- N b ' t S , (l + a b )Nfc, Spr EV ) if N b ' t Spi ey > 2E b , else
- N b' ,S N b' ,S Pr ev ⁇ [00139] That is, the noise estimate is updated only if the previous noise estimate suggests the noise level is greater, e.g., greater than twice the current echo prediction. Otherwise the echo would bias the noise estimate.
- FIGS. 4, 5 and 6 One feature of the noise reducer shown in FIGS. 4, 5 and 6 includes
- the gain calculator 529 includes an element 613 to calculates a probability indicator, expressed as a gain for the intermediate signal, e.g., the frequency bins Y n based on the spatially- selective estimates of the noise power (or other frequency domain amplitude metric) spectrum, and further on the instantaneous banded input power 3 ⁇ 4' in a particular band.
- this probability indicator is referred to as a gain, denoted Gain N . It should be noted however that this gain Gain N is not directly applied, but rather combined with additional gains, i.e., additional probability indicators in a gain combiner 615 to achieve a single gain to apply to achieve a single suppressive action.
- the element 613 is shown with echo suppression, and in some versions does not include echo suppression.
- 3 ⁇ 4' is the instantaneous banded power (or other frequency domain amplitude metric)
- N3 ⁇ 4 s is the banded spatially- selective (out-of-beam) noise estimate
- Some embodiments of input processing for noise reduction include not only noise suppression, but also simultaneous suppression of echo.
- element 613 includes echo suppression and in gain calculator 529, the probability indicator for suppressing echoes is expressed as a gain denoted Gain ' b ⁇ +E .
- the above noise suppression gain expression in the case of also including echo suppression, becomes f max(0,3 ⁇ 4 - p N N b ' iS - fi E ' E b )) GainEx Pb
- Gain b ' N+E (“Gain ⁇ )
- Y b is again the instantaneous banded power
- N b s , E b ' are the banded spatially- selective noise and banded echo estimates
- ⁇ ⁇ ' , ⁇ ⁇ ' are scaling parameters in the range of 1 to 4, to allow for error in the noise and echo estimates and to offset the gain curve accordingly.
- Gain N+E Several of the expressions for Gain N+E described herein have the instantaneous banded input power (or other frequency domain amplitude metric) Y b in both the numerator and denominator. This works well when the banding is properly designed as described herein, with log-like or perceptually spaced frequency bands.
- the denominator uses the estimated banded power spectrum (or other amplitude metric spectrum) P b , so that the above expression for
- the suppression gain expressions above can be generalized as functions on the domain of the ratio of the instantaneous input power to the expected undesirable signal power, sometimes called "noise" for simplicity.
- the undesirable signal power is the sum of the estimated (location-sensitive) noise power and predicted or estimated echo power. Combining the noise and echo together in this way provides a single probability indicator in the form of a suppressive gain that causes simultaneous attenuation of both undesirable noise and of undesirable echo.
- the echo can achieve a level substantially higher than the level of the noise, such suppression may not lead to sufficient echo attenuation.
- an additional scaling of the probability indicator or gain is used, such additional scaling based on the ratio of input audio signal to echo power alone.
- f A ( ⁇ ) , f B ( ⁇ ) a pair of suppression gain functions, each having desired properties for suppression gains, e.g., as described above, including, for example being smooth.
- each of f A ( ⁇ ) , f B ( ⁇ ) has sigmoid function characteristics.
- signal power in yet another embodiment, can be applied for both noise and echo suppression, and o suppression.
- the suppression probability indicator for in-beam signals expressed as a beam gain 612, called the spatial suppression gain, and denoted Gaini' s is determined by a spatial suppression gain calculator 611 in element 529 (FIG. 5) as
- the spatial suppression gain 612 is combined with other suppression gains in gain combiner 615 to form an overall probability indicator expressed as a suppression gain.
- the overall probability indicator for simultaneous suppression of noise, echo, and out-of - beam signals expressed as a gain Gaini' RAW > i s m one embodiment the product of the gains:
- Gainb' Gaini', $ ⁇ Gain ⁇ N+E .
- Gaini', 0 A + 0.9 Gaini', s - Gaini', N+E .
- this noise and echo suppression gain is combined with the spatial feature probability indicator or gain for forming a raw combined gain, and then post-processed by a post-processor 625 and by the post processing step to ensure stability and other desired behavior.
- the gain function specific to the echo is the gain function specific to the echo
- gain calculator 529 includes a determiner of the additional echo suppression gain and a combiner 627 of the additional echo suppression gain with the post-processed gain to result in the overall B gains to apply. The inventors discovered that such an embodiment can provide a more specific and deeper attenuation of echo, since the echo probability indicator or gain is not subject to the smoothing and continuity
- FIG. 7 shows a flowchart of a method 700 of operating a processing apparatus 100 to suppress noise and out-of-location signals and in some embodiments echo in a number P>1 of signal inputs 101, e.g., from differently located microphones.
- method 700 includes processing a Q> ⁇ reference inputs 102, e.g., ⁇ 2 inputs to be rendered on Q loudspeakers, or signals obtained from Q loudspeakers.
- method 700 comprises: accepting 701 in the processing
- apparatus a plurality of sampled input audio signals 101, and forming 703, 707, 709 a mixed-down banded instantaneous frequency domain amplitude metric 417 of the input audio signals 101 for a plurality of frequency bands, the forming including
- the forming includes in 703 transforming the input audio signals to frequency bins, downmixing, e.g., beamforming 707 the frequency data, and in 709 banding.
- the method includes calculating the power (or other amplitude metric) spectrum of the signal.
- the downmixing can be before transforming, so that a single mixed-down signal is transformed.
- the system may make use of an estimate of the banded echo reference, or a similar representation of the frequency domain spectrum of the echo reference provided by another processing component or source within the realized system.
- the method includes determining in 705 banded spatial features, e.g., location probability indicators 419 from the plurality of sampled input audio signals.
- the method includes accepting 713 one or more reference signals and forming in 715 and 717 a banded frequency domain amplitude metric representation of the one or more reference signals.
- the representation in one embodiment is the sum.
- the method includes predicting in 721 a banded frequency domain amplitude metric representation of the echo 415 using adaptively determined echo filter coefficients.
- the predicting in one embodiment further includes voice-activity detecting— VAD— using the estimate of the banded spectral amplitude metric of the mixed-down signal 413, the estimate of banded spectral amplitude metric of noise, and the previously predicted echo spectral content 415.
- the coefficients are updated or not according to the results of voice-activity detecting. Updating uses an estimate of the banded spectral amplitude metric of the noise, previously predicted echo spectral content 415, and an estimate of the banded spectral amplitude metric of the mixed-down signal 413.
- the estimate of the banded spectral amplitude metric of the mixed-down signal is in one embodiment the mixed-down banded instantaneous frequency domain amplitude metric 417 of the input audio signals, while in other embodiments, signal spectral estimation is used.
- the method 700 includes: a) calculating in 723 raw
- suppression gains including an out-of-location signal gain determined using two or more of the spatial features 419, and a noise suppression gain determined using spatially- selective noise spectral content; and b) combining the raw suppression gains to a first combined gain for each band.
- the noise suppression gain in some embodiments includes suppression of echoes, and its calculating 723 also uses the predicted echo spectral content 415.
- the method 700 further includes in 725 carrying out
- spatially-selective voice activity detection determined using two or more of the spatial features 419 to generate a signal classification, e.g., whether voice or not.
- wind detection is used such that the signal classification further includes whether the signal is wind or not.
- the method 700 further includes carrying out post-processing on the first
- the post-processing includes ensuring minimum gain, e.g., in a band dependent manner.
- the post- processing includes carrying out median filtering of the combined gains, e.g., to ensure there are no outlier gains.
- Some embodiments of post-processing include ensuring smoothness by carrying out time and/or band-to-band smoothing.
- the post-processing 725 is according to the signal
- classification e.g., whether voice or not, or whether wind or not, and in some
- the characteristics of the median filtering vary according to the signal classification, e.g., whether voice or not, or whether wind or not.
- the method includes calculating in 726 an additional echo suppression gain.
- the additional echo suppression gain is included in the first combined gain which is used as a final gain for each band, and in another embodiment, the additional echo suppression gain is combined with the results of applying post-processing to the first combined gain to generate a final gain for each band.
- the method includes applying in 727 the final gain, including interpolating the gain for bin data to carry out suppression on the bin data of the mixed-down signal to form suppressed signal data 133, and applying in 729 one or both of a) output synthesis and transforming to generate output samples, and b) output remapping to generate output frequency bins.
- noise reduction is only one example of input processing that determines gains that can be post-processed by the post-processing method that includes median filtering described in embodiments of the present invention.
- a processing system-based apparatus A processing system-based apparatus
- FIG. 8 shows a simplified block diagram of one processing apparatus
- embodiment 800 for processing one or more of audio inputs 101 e.g., from microphones (not shown).
- the processing apparatus 800 is to determine a set of gains, to post-process the gains including median filtering the determined gains, and to generate audio output 137 that has been modified by application of the gains.
- One version achieves one or more of perceptual domain-based leveling, perceptual domain-based dynamic range control, and perceptual domain-based dynamic equalization that takes into account the variation in the perception of audio depending on the reproduction level of the audio signal.
- Another version achieved noise reduction.
- One noise reduction version includes echo reduction, and in such a version, the processing apparatus also accepts one or more reference signals 103, e.g., from one or more loudspeakers (not shown) or from the feed(s) to such loudspeaker(s).
- the processing apparatus 800 is to generate audio output 137 that has been modified by suppressing, in one embodiment noise and out-of-location signals, and in another embodiment also echoes as specified in accordance to one or more features of the present invention.
- the apparatus for example, can implement the system shown in FTG. 6, and any alternates thereof, and can carry out, when operating, the method of FTG. 7 including any variations of the method described herein.
- Such an apparatus may be included, for example, in a headphone set such as a Bluetooth headset.
- the audio inputs 101 , the reference input(s) 103 and the audio output 137 are assumed to be in the form of frames of M samples of sampled data.
- a digitizer including an analog-to-digital converter and quantizer would be present.
- a de- quantizer and a digital-to-analog converter would be present.
- FIG. 8 includes a processing system 803 that is
- the processing system 803 includes at least one processor 805, which can be the processing unit(s) of a digital signal processing device, or a CPU of a more general purpose processing device.
- the processing system 803 also includes a storage subsystem 807 typically including one or more memory elements.
- the elements of the processing system are coupled, e.g., by a bus subsystem or some other interconnection mechanism not shown in FIG. 8. Some of the elements of processing system 803 may be integrated into a single circuit, using techniques commonly known to one skilled in the art.
- the storage subsystem 807 includes instructions 811 that when executed by the processor(s) 805, cause carrying out of the methods described herein.
- the storage subsystem 807 is configured to store one or more tuning parameters 813 that can be used to vary some of the processing steps carried out by the processing system 803.
- the system shown in FIG. 8 can be incorporated in a specialized device such as a headset, e.g., a wireless Bluetooth headset.
- the system also can be part of a general purpose computer, e.g., a personal computer configured to process audio signals.
- the post-processing e.g., the median filtering is controlled by signal classification as determined by a VAD.
- the invention is not limited to any particular type of VAD, and many VADs are known in the art.
- the inventors When applied to suppression, the inventors have discovered that suppression works best when different parts of the suppression system are controlled by different VADs, each such VAD custom designed for the functions of the suppressor in which it is used in, rather than having an "optimal" VAD for all uses. Therefore, in some versions of the input processing for noise reduction, a plurality of VADs, each controlled by a small set of tuning parameters that separately control sensitivity and selectivity, including spatial selectivity, such parameters tuned according to the suppression elements in which the VAD is used.
- Each of the plurality of the VADs is an instantiation of a universal VAD that determines indications of voice activity from Y .
- the universal VAD is controlled by a set of parameters and uses an estimate of noise spectral content, the banded frequency domain amplitude metric representation of the echo, and the banded spatial features.
- the set of parameters includes whether the estimate of noise spectral content is spatially selective or not.
- the type of indication of voice activity that a particular instantiation determines is controlled by a selection of the parameters.
- BeamGain BeamGain m ⁇ n + (1- BeamGain m ⁇ n )RPI ⁇ ,' ⁇ yCPI 3 ⁇ 4,
- BeamGainExp is a parameter that for larger values increases the aggressiveness of the spatial selectivity of the VAD, and is 0 for a non-spatially-selective VAD
- N v N s denotes either the total noise power (or other frequency domain amplitude metric) estimate N , or the spatially-selective noise estimate N s determined using the out-of- beam power (or other frequency domain amplitude metric)
- ⁇ ⁇ , ⁇ ⁇ > 1 are margins for noise end echo, respectively
- Y s ' ens is a settable sensitivity offset.
- the values of ⁇ ⁇ , ⁇ ⁇ are between 1 and 4.
- BeamGainExp is between 0.5 to 2.0 when spatial selectivity is desired, and is 1.5 for one embodiment of a spatially- selective VAD, e.g., used to control post-processing in some embodiments of the invention.
- three spatial probability indicators namely the ratio probability indicator, the phase probability indicator, and the coherence probability indicator.
- the above expression also controls the operation of the universal voice activity detecting method.
- a binary decision or classifier can be obtained by considering the test S > 3 ⁇ 4 rei3 ⁇ 4 as indicating the presence of voice. It should also be apparent that the value S can be used as a continuous indicator of the instantaneous voice level.
- an improved useful universal VAD for operations such as transmission control or controlling the post processing could be obtained using a suitable "hang over" or period of continued indication of voice after a detected event. Such a hang over period may vary from 0 to 500ms, and in one embodiment a value of 200ms was used.
- the noise in the above expression is N s determined using an out-of-beam estimate of power (or other frequency domain amplitude metric).
- Y sens is set to be around expected microphone and system noise level, obtained by experiments on typical components.
- processor may refer to any device or portion of a device that processes electronic data, e.g., from registers and/or memory to transform that electronic data into other electronic data that, e.g., may be stored in registers and/or memory.
- a "computer” or a “computing machine” or a “computing platform” may include one or more processors.
- the methodologies described herein are, in some embodiments, performable by one or more processors that accept logic, instructions encoded on one or more computer- readable media. When executed by one or more of the processors, the instructions cause carrying out at least one of the methods described herein. Any processor capable of executing a set of instructions (sequential or otherwise) that specify actions to be taken is included. Thus, one example is a typical processing system that includes one or more processors. Each processor may include one or more of a CPU or similar element, a graphics processing unit (GPU), field-programmable gate array, application-specific integrated circuit, and/or a programmable DSP unit.
- GPU graphics processing unit
- DSP programmable digital signal processor
- the processing system further includes a storage subsystem with at least one storage medium, which may include memory embedded in a semiconductor device, or a separate memory subsystem including main RAM and/or a static RAM, and/or ROM, and also cache memory.
- the storage subsystem may further include one or more other storage devices, such as magnetic and/or optical and/or further solid state storage devices.
- a bus subsystem may be included for communicating between the components.
- the processing system further may be a distributed processing system with processors coupled by a network, e.g., via network interface devices or wireless network interface devices.
- the processing system requires a display, such a display may be included, e.g., a liquid crystal display (LCD), organic light emitting display (OLED), or a cathode ray tube (CRT) display.
- a display e.g., a liquid crystal display (LCD), organic light emitting display (OLED), or a cathode ray tube (CRT) display.
- the processing system also includes an input device such as one or more of an alphanumeric input unit such as a keyboard, a pointing control device such as a mouse, and so forth.
- the term storage device, storage subsystem, or memory unit as used herein, if clear from the context and unless explicitly stated otherwise, also encompasses a storage system such as a disk drive unit.
- the processing system in some configurations may include a sound output device, and a network interface device.
- a non-transitory computer-readable medium is configured with, e.g., encoded with instructions, e.g., logic that when executed by one or more processors of a processing system such as a digital signal processing device or subsystem that includes at least one processor element and a storage subsystem, cause carrying out a method as described herein. Some embodiments are in the form of the logic itself.
- a non- transitory computer-readable medium is any computer-readable medium that is not specifically a transitory propagated signal or a transitory carrier wave or some other transitory transmission medium. The term "non-transitory computer-readable medium" thus covers any tangible computer-readable storage medium.
- Non-transitory computer- readable media include any tangible computer-readable storage media and may take many forms including non- volatile storage media and volatile storage media.
- Non- volatile storage media include, for example, static RAM, optical disks, magnetic disks, and magneto-optical disks.
- Volatile storage media includes dynamic memory, such as main memory in a processing system, and hardware registers in a processing system.
- the storage subsystem thus a computer- readable storage medium that is configured with, e.g., encoded with instructions, e.g., logic, e.g., software that when executed by one or more processors, causes carrying out one or more of the method steps described herein.
- the software may reside in the hard disk, or may also reside, completely or at least partially, within the memory, e.g., RAM and/or within the processor registers during execution thereof by the computer system.
- the memory and the processor registers also constitute a non- transitory computer- readable medium on which can be encoded instructions to cause, when executed, carrying out method steps.
- the computer-readable medium is shown in an example embodiment to be a single medium, the term "medium” should be taken to include a single medium or multiple media (e.g., several memories, a centralized or distributed database, and/or associated caches and servers) that store the one or more sets of instructions.
- a non-transitory computer-readable medium e.g., a computer- readable storage medium may form a computer program product, or be included in a computer program product.
- the one or more processors operate as a standalone device or may be connected, e.g., networked to other processor(s), in a networked deployment, or the one or more processors may operate in the capacity of a server or a client machine in server-client network environment, or as a peer machine in a peer-to- peer or distributed network environment.
- the term processing system encompasses all such possibilities, unless explicitly excluded herein.
- the one or more processors may form a personal computer (PC), a media playback device, a headset device, a hands-free communication device, a tablet PC, a set-top box (STB), a personal digital assistant (PDA), a game machine, a cellular telephone, a Web appliance, a network router, switch or bridge, or any machine capable of executing a set of instructions (sequential or otherwise) that specify actions to be taken by that machine.
- PC personal computer
- PDA personal digital assistant
- game machine e.g., a single processor and a single storage subsystem, e.g., a single memory that stores the logic including instructions
- embodiments of the present invention may be embodied as a method, an apparatus such as a special purpose apparatus, an apparatus such as a data processing system, logic, e.g., embodied in a non- transitory computer-readable medium, or a computer-readable medium that is encoded with instructions, e.g., a computer-readable storage medium configured as a computer program product.
- the computer-readable medium is configured with a set of instructions that when executed by one or more processors cause carrying out method steps.
- aspects of the present invention may take the form of a method, an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects.
- the present invention may take the form of program logic, e.g., a computer program on a computer-readable storage medium, or the computer-readable storage medium configured with computer-readable program code, e.g., a computer program product.
- the short time Fourier transform is used to obtain the frequency bands
- the invention is not limited to the STFT.
- Transforms such as the STFT are often referred to as circulant transforms.
- Most general forms of circulant transforms can be represented by buffering, a window, a twist (real value to complex value transformation) and a DFT, e.g., FFT.
- a complex twist after the DFT can be used to adjust the frequency domain representation to match specific transform definitions.
- the invention may be implemented by any of this class of transforms, including the modified DFT (MDFT), the short time Fourier transform (STFT), and with a longer window and wrapping, a conjugate quadrature mirror filter (CQMF).
- MDFT modified DFT
- STFT short time Fourier transform
- CQMF conjugate quadrature mirror filter
- MDCT Modified discrete cosine transform
- MDST modified discrete sine transform
- any one of the terms comprising, comprised of or which comprises is an open term that means including at least the elements/features that follow, but not excluding others.
- the term comprising, when used in the claims should not be interpreted as being limitative to the means or elements or steps listed thereafter.
- the scope of the expression a device comprising A and B should not be limited to devices consisting of only elements A and B.
- Any one of the terms including or which includes or that includes as used herein is also an open term that also means including at least the elements/features that follow the term, but not excluding others. Thus, including is synonymous with and means comprising.
- Coupled when used in the claims, should not be interpreted as being limitative to direct connections only.
- the terms “coupled” and “connected,” along with their derivatives, may be used. It should be understood that these terms are not intended as synonyms for each other, but may be.
- the scope of the expression “a device A coupled to a device B” should not be limited to devices or systems wherein an input or output of device A is directly connected to an output or input of device B. It means that there exists a path between device A and device B which may be a path including other devices or means in between.
- “coupled to” does not imply direction.
- a device A is coupled to a device B
- a device B is coupled to a device A
- Coupled may mean that two or more elements are either in direct physical or electrical contact, or that two or more elements are not in direct contact with each other but yet still co-operate or interact with each other.
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Circuit For Audible Band Transducer (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
Description
Claims
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201161441611P | 2011-02-10 | 2011-02-10 | |
PCT/US2012/024372 WO2012109385A1 (en) | 2011-02-10 | 2012-02-08 | Post-processing including median filtering of noise suppression gains |
Publications (2)
Publication Number | Publication Date |
---|---|
EP2673778A1 true EP2673778A1 (en) | 2013-12-18 |
EP2673778B1 EP2673778B1 (en) | 2018-10-10 |
Family
ID=45809594
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP12707413.6A Active EP2673778B1 (en) | 2011-02-10 | 2012-02-08 | Post-processing including median filtering of noise suppression gains |
EP12707412.8A Active EP2673777B1 (en) | 2011-02-10 | 2012-02-08 | Combined suppression of noise and out-of-location signals |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP12707412.8A Active EP2673777B1 (en) | 2011-02-10 | 2012-02-08 | Combined suppression of noise and out-of-location signals |
Country Status (4)
Country | Link |
---|---|
EP (2) | EP2673778B1 (en) |
JP (1) | JP6002690B2 (en) |
CN (2) | CN103354937B (en) |
WO (2) | WO2012109385A1 (en) |
Families Citing this family (42)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012107561A1 (en) | 2011-02-10 | 2012-08-16 | Dolby International Ab | Spatial adaptation in multi-microphone sound capture |
JP5744236B2 (en) | 2011-02-10 | 2015-07-08 | ドルビー ラボラトリーズ ライセンシング コーポレイション | System and method for wind detection and suppression |
US8804977B2 (en) | 2011-03-18 | 2014-08-12 | Dolby Laboratories Licensing Corporation | Nonlinear reference signal processing for echo suppression |
CN103325380B (en) | 2012-03-23 | 2017-09-12 | 杜比实验室特许公司 | Gain for signal enhancing is post-processed |
CN104050969A (en) | 2013-03-14 | 2014-09-17 | 杜比实验室特许公司 | Space comfortable noise |
GB201309777D0 (en) * | 2013-05-31 | 2013-07-17 | Microsoft Corp | Echo suppression |
GB201309773D0 (en) | 2013-05-31 | 2013-07-17 | Microsoft Corp | Echo removal |
GB201309779D0 (en) | 2013-05-31 | 2013-07-17 | Microsoft Corp | Echo removal |
GB201309771D0 (en) | 2013-05-31 | 2013-07-17 | Microsoft Corp | Echo removal |
JP6087762B2 (en) * | 2013-08-13 | 2017-03-01 | 日本電信電話株式会社 | Reverberation suppression apparatus and method, program, and recording medium |
TWI498884B (en) * | 2013-09-09 | 2015-09-01 | Pegatron Corp | Electric device with environment sound filtering function and method for filtering environment sound |
EP2854133A1 (en) | 2013-09-27 | 2015-04-01 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Generation of a downmix signal |
CN104681034A (en) | 2013-11-27 | 2015-06-03 | 杜比实验室特许公司 | Audio signal processing method |
EP3080975B1 (en) * | 2013-12-12 | 2017-07-12 | Koninklijke Philips N.V. | Echo cancellation |
US10079941B2 (en) | 2014-07-07 | 2018-09-18 | Dolby Laboratories Licensing Corporation | Audio capture and render device having a visual display and user interface for use for audio conferencing |
CN104200817B (en) * | 2014-07-31 | 2017-07-28 | 广东美的制冷设备有限公司 | Sound control method and system |
JP6195548B2 (en) * | 2014-08-19 | 2017-09-13 | 日本電信電話株式会社 | Signal analysis apparatus, method, and program |
US9953661B2 (en) * | 2014-09-26 | 2018-04-24 | Cirrus Logic Inc. | Neural network voice activity detection employing running range normalization |
DK3245798T3 (en) * | 2015-01-14 | 2018-08-13 | Widex As | PROCEDURE TO OPERATE A HEARING SYSTEM AND HEARING SYSTEM |
WO2017001611A1 (en) | 2015-06-30 | 2017-01-05 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Method and device for the allocation of sounds and for analysis |
US9401158B1 (en) * | 2015-09-14 | 2016-07-26 | Knowles Electronics, Llc | Microphone signal fusion |
CN107493247B (en) * | 2016-06-13 | 2021-10-22 | 中兴通讯股份有限公司 | Self-adaptive equalization method and device and equalizer |
US9906859B1 (en) * | 2016-09-30 | 2018-02-27 | Bose Corporation | Noise estimation for dynamic sound adjustment |
CN106782591B (en) * | 2016-12-26 | 2021-02-19 | 惠州Tcl移动通信有限公司 | Device and method for improving speech recognition rate under background noise |
US11513205B2 (en) | 2017-10-30 | 2022-11-29 | The Research Foundation For The State University Of New York | System and method associated with user authentication based on an acoustic-based echo-signature |
CN107886965B (en) * | 2017-11-28 | 2021-04-20 | 游密科技(深圳)有限公司 | Echo cancellation method for game background sound |
CN108447500B (en) * | 2018-04-27 | 2020-08-18 | 深圳市沃特沃德股份有限公司 | Method and device for speech enhancement |
US10847178B2 (en) * | 2018-05-18 | 2020-11-24 | Sonos, Inc. | Linear filtering for noise-suppressed speech detection |
CN112437957B (en) * | 2018-07-27 | 2024-09-27 | 杜比实验室特许公司 | Forced gap insertion for full listening |
CN111145770B (en) * | 2018-11-02 | 2022-11-22 | 北京微播视界科技有限公司 | Audio processing method and device |
KR102517975B1 (en) * | 2019-01-29 | 2023-04-04 | 삼성전자주식회사 | Residual echo estimator to estimate residual echo based on time correlation, non-transitory computer-readable medium storing program code to estimate residual echo, and application processor |
CN109949820B (en) * | 2019-03-07 | 2020-05-08 | 出门问问信息科技有限公司 | Voice signal processing method, device and system |
CN110097884B (en) * | 2019-06-11 | 2022-05-17 | 大众问问(北京)信息科技有限公司 | Voice interaction method and device |
CN110136739B (en) * | 2019-07-03 | 2021-04-09 | 四川大学 | Method for removing impulse noise in voice signal |
CN112420073B (en) * | 2020-10-12 | 2024-04-16 | 北京百度网讯科技有限公司 | Voice signal processing method, device, electronic equipment and storage medium |
CN112235693B (en) * | 2020-11-04 | 2021-12-21 | 北京声智科技有限公司 | Microphone signal processing method, device, equipment and computer readable storage medium |
US20240171907A1 (en) * | 2021-02-04 | 2024-05-23 | Neatframe Limited | Audio processing |
GB202101561D0 (en) * | 2021-02-04 | 2021-03-24 | Neatframe Ltd | Audio processing |
AU2022233253A1 (en) * | 2021-03-11 | 2023-09-21 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Audio decorrelator, processing system and method for decorrelating an audio signal |
US11849291B2 (en) * | 2021-05-17 | 2023-12-19 | Apple Inc. | Spatially informed acoustic echo cancelation |
CN113689878A (en) * | 2021-07-26 | 2021-11-23 | 浙江大华技术股份有限公司 | Echo cancellation method, echo cancellation device, and computer-readable storage medium |
EP4198976B1 (en) * | 2021-12-17 | 2023-10-25 | GN Audio A/S | Wind noise suppression system |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2913105B2 (en) * | 1989-03-10 | 1999-06-28 | 日本電信電話株式会社 | Sound signal detection method |
DE4405723A1 (en) * | 1994-02-23 | 1995-08-24 | Daimler Benz Ag | Method for noise reduction of a disturbed speech signal |
JP3484801B2 (en) * | 1995-02-17 | 2004-01-06 | ソニー株式会社 | Method and apparatus for reducing noise of audio signal |
US5587998A (en) * | 1995-03-03 | 1996-12-24 | At&T | Method and apparatus for reducing residual far-end echo in voice communication networks |
US6351731B1 (en) * | 1998-08-21 | 2002-02-26 | Polycom, Inc. | Adaptive filter featuring spectral gain smoothing and variable noise multiplier for noise reduction, and method therefor |
JP3435357B2 (en) * | 1998-09-07 | 2003-08-11 | 日本電信電話株式会社 | Sound collection method, device thereof, and program recording medium |
KR100467020B1 (en) | 2002-07-26 | 2005-01-24 | 삼성전자주식회사 | Semiconductor Device With Self-Aligned Junction Contact Hole And Method Of Fabricating The Same |
JP4247037B2 (en) * | 2003-01-29 | 2009-04-02 | 株式会社東芝 | Audio signal processing method, apparatus and program |
AU2004248544B2 (en) | 2003-05-28 | 2010-02-18 | Dolby Laboratories Licensing Corporation | Method, apparatus and computer program for calculating and adjusting the perceived loudness of an audio signal |
DE602004017603D1 (en) * | 2004-09-03 | 2008-12-18 | Harman Becker Automotive Sys | Speech signal processing for the joint adaptive reduction of noise and acoustic echoes |
US8280730B2 (en) * | 2005-05-25 | 2012-10-02 | Motorola Mobility Llc | Method and apparatus of increasing speech intelligibility in noisy environments |
EP1931169A4 (en) * | 2005-09-02 | 2009-12-16 | Japan Adv Inst Science & Tech | Post filter for microphone array |
US8811627B2 (en) * | 2005-10-26 | 2014-08-19 | Nec Corporation | Echo suppressing method and apparatus |
US8068619B2 (en) * | 2006-05-09 | 2011-11-29 | Fortemedia, Inc. | Method and apparatus for noise suppression in a small array microphone system |
US8005237B2 (en) * | 2007-05-17 | 2011-08-23 | Microsoft Corp. | Sensor array beamformer post-processor |
JP4850191B2 (en) * | 2008-01-16 | 2012-01-11 | 富士通株式会社 | Automatic volume control device and voice communication device using the same |
CN101510426B (en) * | 2009-03-23 | 2013-03-27 | 北京中星微电子有限公司 | Method and system for eliminating noise |
-
2012
- 2012-02-08 WO PCT/US2012/024372 patent/WO2012109385A1/en active Application Filing
- 2012-02-08 JP JP2013553528A patent/JP6002690B2/en active Active
- 2012-02-08 EP EP12707413.6A patent/EP2673778B1/en active Active
- 2012-02-08 EP EP12707412.8A patent/EP2673777B1/en active Active
- 2012-02-08 CN CN201280008279.7A patent/CN103354937B/en active Active
- 2012-02-08 CN CN201280008266.XA patent/CN103348408B/en active Active
- 2012-02-08 WO PCT/US2012/024370 patent/WO2012109384A1/en active Application Filing
Non-Patent Citations (1)
Title |
---|
See references of WO2012109385A1 * |
Also Published As
Publication number | Publication date |
---|---|
WO2012109385A1 (en) | 2012-08-16 |
CN103354937A (en) | 2013-10-16 |
EP2673777A1 (en) | 2013-12-18 |
JP2014510452A (en) | 2014-04-24 |
EP2673778B1 (en) | 2018-10-10 |
CN103348408B (en) | 2015-11-25 |
WO2012109384A1 (en) | 2012-08-16 |
JP6002690B2 (en) | 2016-10-05 |
CN103354937B (en) | 2015-07-29 |
EP2673777B1 (en) | 2018-12-26 |
CN103348408A (en) | 2013-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2673778B1 (en) | Post-processing including median filtering of noise suppression gains | |
EP2880655B1 (en) | Percentile filtering of noise reduction gains | |
US8712076B2 (en) | Post-processing including median filtering of noise suppression gains | |
US9173025B2 (en) | Combined suppression of noise, echo, and out-of-location signals | |
US12112768B2 (en) | Post-processing gains for signal enhancement | |
US10482896B2 (en) | Multi-band noise reduction system and methodology for digital audio signals | |
US8861745B2 (en) | Wind noise mitigation | |
WO2002005262A2 (en) | Sub-band exponential smoothing noise canceling system | |
US10262673B2 (en) | Soft-talk audio capture for mobile devices | |
US20200286501A1 (en) | Apparatus and a method for signal enhancement | |
Ahmed et al. | Comparison of LMS and NLMS algorithm with the using of 4 Linear Microphone Array for Speech Enhancement. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
17P | Request for examination filed |
Effective date: 20130910 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
DAX | Request for extension of the european patent (deleted) | ||
17Q | First examination report despatched |
Effective date: 20150318 |
|
RAP1 | Party data changed (applicant data changed or rights of an application transferred) |
Owner name: DOLBY LABORATORIES LICENSING CORPORATION |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R079 Ref document number: 602012052000 Country of ref document: DE Free format text: PREVIOUS MAIN CLASS: G10L0021020000 Ipc: G10L0021020800 |
|
RIC1 | Information provided on ipc code assigned before grant |
Ipc: G10L 21/0208 20130101AFI20180323BHEP |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: GRANT OF PATENT IS INTENDED |
|
RIN1 | Information on inventor provided before grant (corrected) |
Inventor name: DICKINS, GLENN N. |
|
INTG | Intention to grant announced |
Effective date: 20180504 |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE PATENT HAS BEEN GRANTED |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: EP Ref country code: AT Ref legal event code: REF Ref document number: 1052179 Country of ref document: AT Kind code of ref document: T Effective date: 20181015 |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R096 Ref document number: 602012052000 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: NL Ref legal event code: MP Effective date: 20181010 |
|
REG | Reference to a national code |
Ref country code: LT Ref legal event code: MG4D |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: MK05 Ref document number: 1052179 Country of ref document: AT Kind code of ref document: T Effective date: 20181010 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: BG Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190110 Ref country code: FI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: LT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: NO Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190110 Ref country code: IS Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190210 Ref country code: LV Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: HR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: PL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: ES Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: AT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: PT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190210 Ref country code: GR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190111 Ref country code: RS Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: SE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: AL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R097 Ref document number: 602012052000 Country of ref document: DE |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: DK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: IT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: CZ Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 |
|
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: SK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: RO Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: SM Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: EE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 |
|
26N | No opposition filed |
Effective date: 20190711 |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: PL |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: SI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 Ref country code: LU Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190208 Ref country code: MC Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 |
|
REG | Reference to a national code |
Ref country code: BE Ref legal event code: MM Effective date: 20190228 |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: MM4A |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: CH Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190228 Ref country code: LI Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190228 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190208 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: BE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190228 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: TR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MT Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190208 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: CY Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: HU Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO Effective date: 20120208 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20181010 |
|
P01 | Opt-out of the competence of the unified patent court (upc) registered |
Effective date: 20230512 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: DE Payment date: 20240123 Year of fee payment: 13 Ref country code: GB Payment date: 20240123 Year of fee payment: 13 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: FR Payment date: 20240123 Year of fee payment: 13 |