Noise-robust speech feature processing with empirical mode decomposition

In this article, a novel technique based on the empirical mode decomposition methodology for processing speech features is proposed and investigated. The empirical mode decomposition generalizes the Fourier analysis. It decomposes a signal as the sum of intrinsic mode functions. In this study, we implement an iterative algorithm to find the intrinsic mode functions for any given signal. We design a novel speech feature post-processing method based on the extracted intrinsic mode functions to achieve noise-robustness for automatic speech recognition. Evaluation results on the noisy-digit Aurora 2.0 database show that our method leads to significant performance improvement. The relative improvement over the baseline features increases from 24.0 to 41.1% when the proposed post-processing method is applied on mean-variance normalized speech features. The proposed method also improves over the performance achieved by a very noise-robust frontend when the test speech data are highly mismatched.


Introduction
State-of-the-art automatic speech recognition (ASR) systems can achieve satisfactory performance under wellmatched conditions.However, when there is a mismatch between the train data and the test data, the performance often degrades quite severely.The versatility of everyday environments requires ASR systems to function well in a wide range of unseen noisy conditions.Therefore, noise-robust speech processing technology for recognition is an important research topic.
Many techniques for noise-robustness have been proposed and put to tests.Speech enhancement methods, such as the well-known spectral subtraction [1] and Wiener filters [2], introduce pre-processing steps to remove the noise part or estimate the clean part given the noisy speech signal.Auditory frontend approaches incorporate knowledge of human auditory systems acquired from psychoacoustic experiments, such as critical bands and spectral/temporal masking effects [3,4], in the process of speech feature extraction.Noise-robust feature post-processing techniques, such as cepstral mean subtraction (CMS) [5], cepstral variance normalization (CVN) [6], and histogram equalization (HEQ) [7], aim to convert raw speech features to a form that is less vulnerable to the corruption of adverse environments.
In this article, we study a feature post-processing technique for noise-robust ASR based on the empirical mode decomposition (EMD) [8].Through EMD, a feature sequence (as a function of time) is decomposed into intrinsic mode functions (IMFs).The basic idea behind our proposed method is that the low-order IMFs contain high-frequency components and they are removed based on a threshold estimated from training data.Alternatively, the recombination weights can be decided using other algorithms [9].Since EMD is a temporal-domain technique, a comparison of EMD with common temporal processing techniques is in order.In the RASTA processing of speech [10], a filter combining temporal difference and integration effects is designed.It results in a bandpass filter, which discriminates speech and noise by their difference in temporal properties.The RASTA processing technique is generally considered very effective for both additive and convolution noises.However, a basic assumption underlying any filtering technique is that the signals being processed are approximately stationary, which may not be the case for speech or non-stationary noises.Furthermore, using linear filters implies a decomposition of signal into sinusoidal functions.In contrast, IMFs used in EMD are data driven, so they are theoretically more general than sinusoidal functions, and may lead to better signal-noise decomposition.A comparison between the results of using EMD and RASTA is given in Section 5.In modulation spectrogram approach [11], modulation patterns of temporal envelope signals of the critical-band channels are represented by the amplitudes at 4-Hz (via FFT) dynamically.This representation proves to be robust for syllable recognition under noise corruption.For a different application, critical parameters of central frequency may have to be tuned.In temporal modulation processing of speech signals [12], the DC part of the signal is denoised for better speech detection in noisy conditions, and to provide an SNR estimator via cross-correlation with low modulation-frequency (1.5-Hz) part of the signal.In contrast to the above reviewed methods of temporal processing, we note that the proposed EMD does not assume stationarity of signal, and there are no task-dependent parameters to be tuned when we extract IMFs.
The rest of article is organized as follows.In Section 2, we introduce the formulation of EMD and show that it is a generalization of the Fourier analysis.In Section 3, we design an iterative algorithm to extract IMFs for EMD.In Section 4, we describe the proposed EMDbased feature post-processing method and give a few illustrative examples.Experimental results are presented in Section 5. Finally, concluding remarks are summarized in Section 6.

Empirical mode decomposition
The EMD generalizes the Fourier series.Sinusoidal basis functions used in the Fourier analysis are generalized to data-dependent IMFs.Compared to a sinusoidal function, an IMF satisfies the generalized alternating property and the generalized zero-mean property, and relaxes the amplitude and frequency from being constant to being generally time-varying.

The Fourier series
A signal x(t) of finite duration, say T, can be represented by a Fourier series, which is a weighted sum of complex exponential functions with frequencies ω k = 2πk/T.That is, we can write (1) Defining we can re-write (1) as If x(t) is real, p k , q k in (2) are real.Equation 3 can be seen as a decomposition of x(t) in the vector space spanned by the following basis set: The following properties of about basis functions of the Fourier series are quite critical in the generalization to EMD.
• (alternating property) A basis function has alternating stationary points and zeros.That is, there is exactly one zero between two stationary points, and exactly one stationary point between two zeros.
• (zero-mean property) The maxima and minima of the basis functions are opposite in sign, and the average of the maxima and the minima is 0.

Empirical mode decomposition
In EMD, a real-valued signal x(t) is decomposed as Those c k (t)s in ( 5) are called IMFs.As generalization of sinusoidal function, an IMF is required to satisfy the following generalized properties.
• (generalized alternating property) The difference between the number of extrema (maxima and minima) and the number of zeros is either 0 or 1.
• (generalized zero-mean property) The average of the upper envelope (a smooth curve through the maxima) and the lower envelope (a smooth curve through the minima) is zero.
The amplitude and frequency of an IMF are defined as follows.Given a real-valued function c k (t), let d k (t) be the Hilbert transform of c k (t).A complex function f k (t) is formed by In (6), we identify a k (t) and ν k (t) as the time-dependent amplitude and the time-dependent frequency of f k (t).Note that the Fourier analysis is a special case of (6), since sinω k t is the Hilbert transform of cosω k t.While sinusoidal functions have constant amplitudes and frequencies, IMFs have time-varying amplitudes and frequencies.

Intrinsic mode functions
The core problem for EMD is to find IMFs given a signal.In the following subsections, we state the algorithm that we design for EMD and highlight properties of IMFs with an illuminating instance.

Extraction algorithm
The pseudocode of the extraction of IMFs is stated as follows. Require In this algorithm, there is an outer loop to control the number of IMFs and there is an inner loop to find the next IMF given the current remainder function.The spline interpolation is used to find envelopes (cf.Section 4.2).To guard against slow convergence, we enforce a criterion to terminate the iteration if the difference between the old and new candidates of h(t) is below a threshold.

An important property
In the extraction of IMFs, the remainder function r(t) is recursively replaced by the hypothetical function h(t), The envelopes u(t) and l(t) are smoother than r(t) as each envelope is the spline interpolation of a proper subset of points of r(t).Being the remainder after the subtraction of the envelope mean, h(t) approximates the time-varying local high-frequency part of r(t).Whenever h(t) is a valid IMF, it is set to c k (t) and subtracted, so the remaining part of signal is smoother.Thus, we expect c k (t) to be progressively smooth as k increases.
For an illustrative example, IMFs extracted from the log-energy sequence of an utterance in the Aurora 2.0 database with a signal-to-noise ratio (SNR) of 0 dB are shown in Figure 1.One can see clearly that the degree of oscillation decreases as k increases, as is predicted by our analysis.

EMD-based feature post-processing
The goal of speech feature post-processing is to reduce the mismatch between clean speech and noisy speech.In order to achieve this goal, we first look at the patterns introduced by the presence of noises of varying levels, then we propose a method to counter such patterns.
The patterns created by noises of several SNRs can be observed on the log-energy sequences of an underlying clean utterance in the Aurora 2.0 database, as shown at the top of Figure 2. We can see that the degree of oscillation of the speech feature sequence increases with the noise level.That is, the spurious spikes in the sequence basically stems from the noise signal, rather than from the speech signal.

Basic idea
Since the spikes introduced by the noise are manifest in the low-order IMFs, we propose to subtract these IMFs to alleviate mismatch.That is, for a noisy speech signal x(t) with EMD we simply subtract a small number, say N, of IMFs from x(t), i.e., At the bottom of Figure 2, EMD post-processed sequences of the same instances are shown.Comparing them to the original sequences at the top, we can see that the mismatch between clean and noisy speech is significantly reduced.

Implementation details
The spline interpolation is employed to find upper and lower envelopes during the process of IMF extraction.For upper envelopes (and similarly for lower envelopes), we use the local maximum points and the end points as the interpolation points.These interpolation points divide the entire time span into segments, and each segment, say segment i, is interpolated by a cubic function, where the parameters a i , b i , g i , δ i can be decided by requiring the overall interpolation function to be continuous up to the second-order derivatives [13].
In the extraction algorithm, we also guard against perpetual changes in the extraction process of IMF via a threshold on the standard deviation (SD), which is defined as follows: where T is the total number of points in the sequence, h o (t) and h n (t) are the old and new candidates for IMF.
In our experiments, we set s = 0.25 [8].The number of iterations needed to find the first IMF varies with the input signal.The histogram (statistics) of this iteration scheme applied on a data set is given in Figure 3.

Experiments
The proposed EMD-based approach to noise-robustness is evaluated on the Aurora 2.0 database [14].After the baseline results are reproduced, we first apply the commonly used per-utterance mean-variance normalized (MVN) on the speech features to boost the performance, then we apply the proposed EMD-based postprocessing to achieve further improvement.After seeing significant performance gain over the baseline, we apply the proposed method to ETSI advanced frontend (AFE) speech features [15] to see if further improvement can be achieved on speech features that are already very noise-robust to begin with.We also compare EMD with the RASTA processing method.

Aurora database
The Aurora 2.0 noisy digit database is widely used for the evaluation of noise-robust frontends [14].Eight types of additive noises are artificially added to clean speech data with SNR levels ranging from 20 to -5 dB.The data may be further convolved with two types of convolution noises.The multi-train recognizer is trained by a data set (called the multi-train set) consisting of clean and multicondition noisy speech samples.The clean-train recognizer is trained by a data set (called the clean-train set) consisting of clean speech samples only.Test data in Set A are matched to the multi-condition train data, test data in Set B are not matched to the multi-condition train data, and test data in Set C are further mismatched due to convolution.Note that the proportion of the data amounts of Set A, Set B, and Set C is 2 : 2 : 1.

Frontend and backend
The baseline speech feature vector consists of the static features of 12 mel-frequency cepstral coefficients (MFCC) C 1 , . .., C 12 and the log energy.Dynamic features of velocity (delta) and acceleration (delta-delta) are also derived, resulting in a 39-dimension vector per frame.
The standard backend recognizer of Aurora evaluation [14] is adopted.That is, we use 16-state whole-word models for digits, a 3-state silence model, and a 1-state short-pause model.The state of the short-pause model is tied to the middle state of the silence model.The state-emitting probability density is a 3-component Gaussian mixture for a word state, and a 6-component Gaussian mixture for a silence/short-pause state.

Results
Three sets of experiments have been carried out in this research.In the first set of experiments, noisy feature sequences are replaced by the corresponding clean feature sequences.This is possible in Aurora 2.0 because clean and noisy speech data are "parallel", i.e. each noisy utterance has a corresponding clean utterance.The results are compared to case where each sequence is post-processed by EMD.In the second set of experiments, various aspects of EMD are investigated.In the final set of experiments, the proposed EMD method is compared to the well-known temporal filtering method of RASTA.

Feature replacement experiments
The first set of experiments is designed to determine which speech feature sequence is to be applied the EMD-based post-processing.For each of the 13 static features, we replace noisy feature sequences with clean feature sequences (RwC: replace with clean).Based on the results summarized in Table 1, it is clear that replacing noisy log-energy sequence leads to the most significant improvement.The performance level decreases as we move down the table from C 1 to C 12 .Thus, unless otherwise stated, in the remaining investigation, we focus on using log-energy sequences as the targets to be processed by the proposed EMD.In addition, we apply the proposed EMD to noisy feature sequences and the results are also shown in Table 1.It is interesting to see that EMD even leads to better performance than clean feature replacement in the cases from C 2 to C 12 .Furthermore, applying EMD to all features does not yield better performance than EMD on log energy alone, although the performance levels are quite close.Higher-order cepstral features provide information for the more delicate structures in the speech signal.It is more difficult to recover such information lost in the presence of noise through EMD.In contrast, the loss of information conveyed by log energy due to noise is relatively easy to recover.

Effectiveness of EMD
The recognition accuracy rates of clean-train tasks averaged over 0-20 dB noisy test data with different degrees of feature post-processing are listed in Table 2.The row of "baseline" shows the results of using the raw speech features extracted by the ETSI standard frontend.The row of "MVN" shows the results after the application of the mean-variance normalization (MVN).MVN achieves 24.0% relative improvement.
The proposed EMD-based method is applied to the log-energy feature sequences, by subtracting the first IMF for each utterance.Applying EMD on the MVN feature sequence, the relative improvement improves from 24.0 to 41.1%.The results show that the EMDbased post-processing of subtracting IMFs from the speech feature sequences significantly reduce the mismatch between clean and noisy feature sequences.
It is very encouraging to see that the case of most significant improvement by EMD is with Set C (66.4-75.3%).We note that Set C contains arguably the most mismatched data because that convolution noises are applied to the utterances in addition to additive noises.With only MVN, the accuracy level on Set C is significantly below Set A or Set B. After EMD, the accuracy levels of the three sets become very close.Thus, EMD does increase the noise-robustness of the ASR system.
Detailed comparison between the word accuracy rates of the MVN method and the proposed EMD-based post-processing method are broken down in Table 3.In addition, we present a scatter plot of the word accuracy rates in Figure 4.It can be clearly seen that the recognition accuracy is improved by EMD.
In addition to ETSI basic frontend feature sequence, we also apply the proposed EMD-based method on ETSI AFE feature sequence.It is important for us to point out that AFE is a strongly noise-robust frontend, which combines modules for voice activity detection (VAD), Wiener-filter noise reduction, and blind  equalization.From Table 3, we can see that while AFE already achieves a relative improvement of 67.1% over the baseline, the application of EMD further improves the performance, achieving further improvements in Sets A and C. The improvement on the most mismatched test data set (Set C) is the most significant (from 85.6 to 86.1%).We also compare subtracting different numbers of IMFs.Essentially, the more IMFs are subtracted, the smoother the resultant sequence becomes.Recognition accuracies when subtracting 1 IMF (MVN+EMD1) and 2 IMFs (MVN+EMD2) are listed in Table 4. From the results, we can see that for the noisier 0 and -5 dB data, MVN+EMD2 yields better accuracy.The results confirm that we should subtract fewer IMFs in higher SNRs, because the interference of noise is not as severe as in lower-SNR cases.
Based on the arguments given in Section 4, it is clear that the noise level and the number of IMFs to be subtracted from the signal to reduce mismatch are closely related.Therefore, we use a scheme that allows the  The scatter plot of word accuracy before and after EMD process.The x-axis is the word accuracy rate before the proposed EMD processing and the y-axis is the word accuracy rate after the EMD processing.A point in the plot corresponds to a test data subset in the Aurora 2.0 corpus, and there are 70 points.The line is x = y, so we can see that the EMD processing technique improves the recognition accuracy.number of IMFs to be subtracted from speech feature sequences to vary from utterance to utterance.We calculate the average oscillation frequency of the log-energy feature sequences from the clean-train data and use it as a threshold.If the oscillation frequency of the remainder is lower than the threshold, we stop finding and subtracting the next IMF.The results of recognition experiments are listed in Table 5.We can see that this scheme, denoted by MVN+EMDd, does outperform the schemes of subtracting a fixed number (1 or 2) of IMFs.We also inspect the number of IMFs, N in (9), subtracted in the dynamic scheme of EMD. Figure 5 shows the average of N on the test set as a function of SNR, for the MVN feature and the AFE feature.As expected, it increases as SNR decreases, i.e., as the noise level increases.

EMD and RASTA
Since EMD is essentially a technique that alters feature sequences in the temporal domain, it is of interest to compare its effectiveness with common temporaldomain techniques.The proposed EMD method is compared to the RASTA processing since both are temporal processing techniques.The results are summarized in It is important to point out that EMD processing is an utterance-level method, so the latency is generally longer than using frame-level methods such as the  RASTA filter or the advanced ETSI frontend.There is a trade-off between complexity, latency, and accuracy.In certain scenarios where low latency is critical, fast online/sequential methods without significant sacrifice in performance may be preferred to batch techniques.

Conclusion
In this article, we propose a feature post-processing scheme for noise-robust speech recognition frontend based on EMD.We introduce EMD as generalization of the Fourier analysis.Our motivation is that speech signals are non-stationary and non-linear, so EMD is theoretically superior to Fourier analysis for signal decomposition.We implement an algorithm to find IMFs.Based on properties of the extracted IMFs, we propose to subtract low-order IMFs to reduce the mismatch between clean and noisy data.Evaluation results on the Aurora 2.0 database show that the proposed method can effectively improve recognition accuracy.Furthermore, with the ETSI AFE speech features, which are very noise-robust by design, the application of EMD method further improves recognition accuracy, which is very remarkable.

Figure 1 Figure 2
Figure1The intrinsic mode functions extracted from the log-energy sequence of the utterance MKG_677884ZA, which is corrupted by the subway noise with the signal-to-noise ratio of 0 dB.

Figure 3
Figure3The histogram of the number of iterations needed to find the first IMF c 1 (t) for the 8440 utterances of clean-train dataset of Aurora 2.0.The actual counts are 171, 3382, 3866, 878, 126, and 17.

Figure 4
Figure 4  The scatter plot of word accuracy before and after EMD process.The x-axis is the word accuracy rate before the proposed EMD processing and the y-axis is the word accuracy rate after the EMD processing.A point in the plot corresponds to a test data subset in the Aurora 2.0 corpus, and there are 70 points.The line is x = y, so we can see that the EMD processing technique improves the recognition accuracy.

Table 1
The word accuracy rates of clean-train tasksThe noisy feature sequences are replaced with the clean feature sequences or they are processed by the proposed EMD-based method.Each number in the table is the average word accuracy over 10 test subsets, 4 subsets from Set A, 4 subsets from Set B, and 2 subsets from Set C for each SNR.RwC: Replaced with Clean; L.E.: the log-energy sequence; C i : the ith MFCC sequence; all: the entire feature vector; none: no replacement or post-processing (baseline).

Table 2
Word accuracy rates of the Aurora 2.0 clean-train tasks for the 0-20 dB SNR test data, using the proposed method Baseline: raw features; MVN: mean variance normalization; MVN+EMD(e): EMD applied on the log-energy sequence; AFE: advanced frontend; AFE+EMD: EMD applied on the combined C 0 /log-energy sequence of AFE.

Table 3
The word accuracy rates of MVN and the proposed EMD method for every noise condition and every test subset (70 subsets total) of the Aurora 2.0 clean-train tasks

Table 6 ,
and it is clearly seen that EMD outperforms the RASTA in this evaluation.The results support our analysis in Section 1 from the theoretical perspective that EMD is potentially more effective on non-stationary signals than conventional techniques based on temporal filtering.Decomposition with IMFs is more general than decomposition with sinusoidal functions, in allowing time-varying amplitudes and frequencies for input signals.

Table 5
Word accuracy rates of Aurora 2.0 clean-train tasks for 0-20 dB SNR test data The average of the number of IMFs extracted for the MVN and AFE features as a function of SNR.For each utterance, the extraction of IMFs stops when the oscillation of r(t) is below a threshold determined by the train data set.It is clear that the average number increases with the level of noise.

Table 6
The comparison of RASTA and the proposed EMD method for every noise condition and every test subset (70 subsets total) of the Aurora 2.0 clean-train tasks