Hybridisation of the bi-dimensional empirical mode decomposition (BEMD) with denoising techniques has been proposed in the literature as an effective approach for image denoising. In this Letter, the Student's probability density function is introduced in the computation of the mean envelope of the data during the BEMD sifting process to make it robust to values that are far from the mean. The resulting BEMD is denoted tBEMD. In order to show the effectiveness of the tBEMD, several image denoising techniques in tBEMD domain are employed; namely, fourth order partial differential equation (PDE), linear complex diffusion process (LCDP), non-linear complex diffusion process (NLCDP), and the discrete wavelet transform (DWT). Two biomedical images and a standard digital image were considered for experiments. The original images were corrupted with additive Gaussian noise with three different levels. Based on peaksignal-to-noise ratio, the experimental results show that PDE, LCDP, NLCDP, and DWT all perform better in the tBEMD than in the classical BEMD domain. It is also found that tBEMD is faster than classical BEMD when the noise level is low. When it is high, the computational cost in terms of processing time is similar. The effectiveness of the presented approach makes it promising for clinical applications.
Introduction:The bi-dimensional empirical mode decomposition (BEMD) is the two-dimensional (2D) version of the original empirical mode decomposition of Huang et al. [1] which is an adaptive multi-resolution analysis technique that decomposes a given signal into a finite sum of components called intrinsic mode functions (IMFs), plus a residue. The empirical mode decomposition algorithm is based on a sifting process that successively estimates the mean envelope of an unknown signal as the average of the upper and lower envelopes, minus a residual component. In particular, low order IMFs represent fast oscillations (high frequency modes), and high order IMFs represent slow oscillations (low frequency modes). The IMFs are local and auto-adaptive to the input data. Hence, the merit of using the EMD is that it is driven by the input data. Unlike the wavelet transform, the EMD does not require a-priori basis function selection, since the IMFs are automatically derived from the input signal. In addition, the EMD is effective for the analysis of both stationary and non-stationary signals, and it is a non-linear decomposition technique.To address the problem of mean sensibility to outliers during the sifting process, an adjusted empirical mode decomposition method built on Student's probability density function (PDF) denoted by tEMD [2] where the data was transformed with the Student PDF prior to computation. Because Student's PDF has longer tails than the normal distribution, it is more robust to values that are far from the mean. In particular, it is less susceptible to signal variance fluctuations, which may lead to more representative mean envelope curves when applying the EMD algorithm. The tEMD was found to be mo...