International Journal of Antennas and Propagation

Volume 2015 (2015), Article ID 769478, 12 pages

http://dx.doi.org/10.1155/2015/769478

## Enhancing Image Denoising Performance of Bidimensional Empirical Mode Decomposition by Improving the Edge Effect

^{1}School of Computer and Communication Engineering, University of Science and Technology Beijing, Beijing 100083, China^{2}Department of Civil Engineering, North China Institute of Science and Technology, Yanjiao, Beijing 101601, China^{3}School of Government, Central University of Finance and Economics, Beijing 100081, China

Received 20 April 2015; Accepted 31 August 2015

Academic Editor: Atsushi Mase

Copyright © 2015 Feng-Ping An et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Bidimensional empirical mode decomposition (BEMD) algorithm, with high adaptive ability, provides a suitable tool for the noisy image processing, and, however, the edge effect involved in its operation gives rise to a problem—how to obtain reliable decomposition results to effectively remove noises from the image. Accordingly, we propose an approach to deal with the edge effect caused by BEMD in the decomposition of an image signal and then to enhance its denoising performance. This approach includes two steps, in which the first one is an extrapolation operation through the regression model constructed by the support vector machine (SVM) method with high generalization ability, based on the information of the original signal, and the second is an expansion by the closed-end mirror expansion technique with respect to the extrema nearest to and beyond the edge of the data resulting from the first operation. Applications to remove the Gaussian white noise, salt and pepper noise, and random noise from the noisy images show that the edge effect of the BEMD can be improved effectively by the proposed approach to meet requirement of the reliable decomposition results. They also illustrate a good denoising effect of the BEMD by improving the edge effect on the basis of the proposed approach. Additionally, the denoised image preserves information details sufficiently and also enlarges the peak signal-to-noise ratio.

#### 1. Introduction

In these years, digital image processing technology has been widely used in space exploration, remote sensing, biomedical and industrial inspection, and various other fields. Although it has made a great progress in many aspects, there are still some new challenges. One of them is closely related to image denoising, since noises can be introduced by almost all of the factors in the process of signal acquisition, conversion, and transfer, not to mention those involved in the original image from the instruments themselves, and environments. In order to obtain real information for a piece of image, denoising generally constitutes a necessary task before the implementation of the postanalysis (such as compression, split, and merge), aiming at improving the signal-to-noise ratio (SNR) and restoring true information of the image.

In the process of image denoising, it is necessary to eliminate the noise contaminating a given image as far as possible and also to preserve all important details of the image. To meet these requirements, a number of theoretical methods have been proposed and developed [1–5]. Relatively, techniques based on Fourier transform and wavelet transform have received more attention than the others in applications to image denoising [3, 6–8]. At the same time, it is also found that they have some limitations. The methods based on Fourier transform can run well for parts of the texture with some deterministic characteristics or of the change with gentle slopes, while they cannot give a satisfactory performance for those parts of the texture with sharp changes in the image. Although wavelet transform provides a useful description in time-frequency domain for some sudden changes of the image contents or slow ones of the texture, a complete description for both of them is still a hard work.

With the development of artificial intelligence algorithms, there have appeared a number of new attempts in image denoising. From a new view, the application of artificial neural network (ANN) fully shows advantages of parallel computing, nonlinear mapping, and adaptive capacity in image processing [9–13]. A good denoising effect especially can be presented by the method of pulse coupled ANN [9, 10, 12]. Another representative technique is the so-called support vector machine (SVM) method. Based on statistical learning theory and structural risk minimization principle, with a high degree of generalization ability, fascinating learning efficiency, and other characteristics [14], it has attracted attention of the denoising researches. Unfortunately, almost all of these methods based on artificial intelligence algorithms require a requisite of signal preprocessing so as to obtain optimal model assumptions and simultaneously to carry out the model training processes. Consequently, there are some difficulties to make full use of their adaptive ability, which limit their extensive applications in image denoising.

For the adaptive ability of image denoising, one of its most important senses addresses the idea that its denoised version can preserve the image features as far as possible, with little influence on the following analysis and processing of the image. However, such influence, even the image distortion, can be found after the actual denoising operation. It is just a view of the limitation of the techniques based on traditional Fourier transform, wavelet transform, and modern artificial intelligence algorithms. To be worthy of note, Huang et al. [15] proposed a feasible technical approach, referred to as empirical mode decomposition (EMD) method, able to solve this problem to some extent. The EMD is a technique suitable for nonstationary time-frequency analysis and has been applied in earthquake [16], structural diagnosis [17], biological [18], mechanical fault diagnosis [19, 20], marine [21], and other fields. Compared with traditional Fourier transform and wavelet transform, it abandons the requirement of choosing* a priori* basis function, completely relying on the data feature itself through data-driven implementation of the decomposition process. Thus, the local variation characteristics may be seen clearly in time-frequency domain from the signal by using this method.

About ten years ago, Nunes et al. [22, 23] introduced the EMD into image processing and proposed the algorithm of bidimensional empirical mode decomposition (BEMD). Soon afterward, this algorithm received attention from some researchers and was put into image compression [24, 25], image texture classification [26], image denoising [27], image texture segmentation [28, 29], image scaling [30], texture synthesis [31], and image feature extraction [32]. Due to the noise characteristics of practical images, it is generally likely that the highest noise content of an image is extracted into the first mode component for the natural property of the BEMD process, and the remaining fewer noise contents sneak into the subsequent ones. On the basis of this fact, the image denoising algorithm by using the BEMD is devised [23, 27]. It is known that the BEMD has a close relation with correct interpolation operation. But, at the boundary region, incorrect interpolation may occur for lack of constraints from extrema points. This is the so-called edge effect and also an open problem in the BEMD. However, it has not been solved so well that further works are required urgently in the application of the BEMD.

In fact, for the decomposition of a one-dimensional signal by the EMD, the edge effects or end effects may be deleted more or less by some strategies such as mirror expansion [33], neural networks expansion [34, 35], ARMA model prediction of time series [36, 37], polynomial fitting extrapolation [38], and waveform feature match expansion [39]. Once these end effects are mitigated, more satisfactory decomposed results can be obtained, which shows that treatment of the end effect is beneficial for the application of the EMD. Accordingly, improvement of the edge effects for the image signals also should be significant to the development of the image denoising by using the BEMD.

To this end, we propose an approach of improving the edge effect of the BEMD to enhance its image denoising performance. This approach applies the SVM method [40, 41] to construct a regression model of the digital image. By using this model, an extrapolation operation is done to have a guess of the image beyond its boundaries. The extrema of the extrapolation signals neighboring the edges are taken as positions of the mirror, and then the closed-end mirroring expansion is introduced to further remove the edge effects of the BEMD. Some case studies of removing the Gaussian white noise, salt and pepper noise, and random noise from the noisy images illustrate a good denoising result of the BEMD coupled with the proposed approach, and the denoised version of the image can also maintain its original details to the full extent.

#### 2. Fundamentals of BEMD Algorithm

Similar to the EMD of the one-dimensional case, the BEMD uses the extrema that are found in the original image or obtained from the first derivative of the original or the higher-order derivative, to achieve the decomposition of the image signal. Distances between extrema may provide the information for characterizing the image on intrinsic length scales.

For a two-dimensional image denoted by , the basic procedure of the BEMD may be summarized as follows [15, 18, 22, 23]:(1)Initialize the image under consideration, , , , where and represent numbers of the rank on the discrete image plane.(2)Initialize the parameters, , and .(3)Extract extrema points of .(4)By cubic spline, interpolate between local maxima and between minima, respectively, to get two envelope surfaces and of .(5)Calculate the mean envelope surface in terms of these two envelope surfaces, given by(6)Update the original signal and designate a new one for iteration, given by(7)Calculate the standard deviation, SD, given by(8)Repeat steps to until the calculated standard deviation is less than a predetermined criterion (SD) (generally taken as 0.2–0.3) and stop the iteration. Right now, it is regarded that represents an intrinsic mode function (IMF); that is, .(9)Update this signal and obtain the residual signal(10)Repeat steps to over and over again, up to the th times, ; when the residual is a monotonic signal, stop the process of BEMD and finally obtain all IMF components.

Once the decomposition process has finished, the original image can be expressed as the sum of all IMFs and the residual , given by

Stopping criterion directly decides the end of the decomposition process and also determines the corresponding final residual term. Therefore, the choice of the value of SD has a close relation with the number of the IMFs obtained by the BEMD. It is necessary to properly choose the value of SD on the basis of the actual situation. In order to portray details of the image, it should be a best choice of a lower value of the SD, generally taken as 0.2–0.3.

In decomposing a one-dimensional signal through the EMD, the residual presents a trend varying with time, which can be ignored in the following analyses. However, the residual resulting from the decomposition of a two-dimensional signal through the BEMD often contains some information of the original image characteristics. It cannot be ignored in the following works, and its influence must be considered for some contributions to the original image.

#### 3. Treatment of the Edge Effect

It is an open problem that there is a certain edge effect in the process of the BEMD. Generally, an actual image signal is not so long or even very short, and serious edge effects can be seen sometimes for its decomposition. Such decomposition is a process screening one after another for many times. With the progression of this process, the edge effect will be more serious, leading to distortions of the IMFs of the image signal decomposed by the BEMD. To this end, some treatments should be taken to mitigate or eliminate edge effects occurring in the screening process of the BEMD operation. Combining advantages of the SVM method with those of the mirror expansion, a new approach to deal with the edge effect caused by BEMD in the decomposition of an image signal is proposed as follows.

##### 3.1. SVM Regression Model and Its Extrapolation

###### 3.1.1. Regression Model

Given training sampling data , where is an input vector of columns, , and is the output value corresponding to , , may be transformed into a high dimensional feature space by the SVM regression through a nonlinear mapping . In this feature space, one can operate a linear regression to get its regression function, expressed asIn this equation, represents the inner product, is the complexity of the function described by the original data and transforming them into a high-dimensional space , and is a constant, as the bias term, . Equation (6) requires constraint conditions:The objective function iswhere and are the slack variables, denoting the upper and lower errors of the training under the constraint of the error and is a constant to control the penalty degree for the sampling data out of , . Vapnik proposed that the risk is measured by the -insensitive loss function as follows [40]:When the absolute value of the difference between the actual and predicted values is less than , the difference for both of them will be negligible, or otherwise it should be equal to .

The optimization problem presented by (8) may be solved by means of Lagrangian penalty method, in which the Lagrangian function is expressed as [40]where , , , and are Lagrange multipliers and all of them are greater than zero. According to the KKT conditions, one has [40]At last, the nonlinear regression problem is transformed into a dual problem of (7), given bywhere () is the kernel function of the SVM. After solving this optimization problem, one can determine parameters and . After that, substituting them into (11) yields the other parameters of the SVM and finally gets the regression model as follows:

Since the regression model is driven by and based on characteristics of the training data, an obvious adaptive ability can be effectively guaranteed in the regression process if appropriate input nodes and related parameters are chosen. Therefore, the SVM regression modeling may provide a good foundation for the expansion of the image data in the following operation.

###### 3.1.2. Parameter Optimization

In order for support vector machine (SVM) algorithm to have a better effect, this paper uses the particle swarm algorithm to improve support vector machine (SVM) parameter, such as , , and . The basic principle is as follows.

The algorithm divides particles into two parts according to the fitness value; the fitness value is the standard of measuring the particle’s quality. It distinguishes the particle’s quality of the wheat from the chaff through comparing a single particle fitness value and the group average fitness value and then adjusts reasonably the particle’s search ability and learning abilities , . It is advantageous to the particle swarm fast convergence and improves the accuracy of particle swarm search. is randomness; it can make , avoid convergence of all particles in a certain direction and reduce the likelihood of being trapped in local optimum. In the end, it gets the optimal parameters of the support vector machine (SVM) model. The calculation method of , , and is as follows:

In the above formula, , , , , , and are constant numbers from 0 to 1, , is a random number from 0 to 1, is the th particle fitness value of the* iter* generation, and is the average fitness of all particles of the* iter* generation.

When , the particle’s quality is poor and (14) will reduce the self-learning ability by enhancing the global search ability and the social learning ability . It will make the particles “fly” quickly toward the direction of the optimal solution and improve the quality of particles. When , the particle’s quality is better and (15) will reduce social learning ability by enhancing the particle’s local search ability and self-learning ability . The algorithm is conducive to finding the global optimal solution of support vector machine (SVM) algorithm parameters and reducing the likelihood of being trapped in local optimum.

###### 3.1.3. Extrapolation Expansion

For a two-dimensional image , it can be described by a certain type of digital signal and forms a data matrix. In the theoretical frame of the SVM modeling, the data can be used to construct the training sampling data and to get the corresponding regression model. After that, the regression model is extrapolated towards outside regions of the image to realize its extrapolation expansion whose procedure may be summarized as follows.(1)For all of the sampling data in the matrix, indexes of the column are regarded as discrete variable in the SVM regression model, and elements of the matrix are output values of for the model.(2)Choosing these data as a sequence of the sampling, one can calculate the related parameters in the SVM modeling and then construct the corresponding SVM regression model, where the radial basis kernel (RBF) is used as the kernel function and the penalty coefficient takes a value of 10.(3)The regression model trained in step extrapolates data values toward the outside regions on the left and right sides of the image, respectively, where represents the index of the expanded column and is the forecasting value of extrapolation at the position .(4)Further, indexes of the row are regarded as the variable of SVM regression model, and elements of the matrix are values of corresponding to different .(5)The sampling data may be given and, similar to step , applied to construct SVM regression model.(6)The obtained regression model trained from step extrapolates data values toward the outside regions of the image up and down, respectively, where represents the row number and is the forecasting value at the position .(7)All the forecasting values resulting from steps and are added into the original matrix according to the corresponding indexes of the expanded rows and columns, which gives the extrapolation version of edge expansion of the given two-dimensional image.

##### 3.2. Closed-End Mirror Expansion of the Extrapolation Version

After the extrapolation expansion of the SVM regression model is carried out, the closed-end mirror expansion method will be applied for the two-dimensional image. As a further expansion, it makes full use of the known information involved in the original data for the extrapolation has provided larger reliable expansion regions. So it may have a better performance for mitigating the edge effect caused by the decomposition process of BEMD.

In operating the closed-end mirror expansion of the extrapolation version resulting from the SVM regression model, the procedure can be briefly summarized as follows:(1)At first, identify if the forecasting values in the expanded columns are local extrema or not. If it is true, the extrapolation operation driven by the SVM regression model is terminated, or otherwise, the operation continues to get data values until a local extremum is captured.(2)Similarly, identify the forecasting values in the expanded rows. If there exist some local extrema, one or more than one, the extrapolation is terminated, or otherwise, it continues to get data values until a local extremum is captured.(3)Finally, put “mirrors” on those positions of local extrema close to the edges determined in steps and and expand the extrapolation version of the image by traditional closed-end mirror expansion technique. In this way, a close sequence is formed, which will be used to realize the operation of the BEMD.

#### 4. Image Denoising by Using the BEMD

##### 4.1. Image Denoising

In the process of image acquisition, compression, conversion, and transmission, some noises may be mixed up with pure image information. Therefore, the real signal of the image may be generally regarded as a compound of the noise and pure image information, which can be expressed aswhere is the real signal of the image with noises, is the image signal, and is the noise.

If the BEMD is operated to these two types of signals, assuming that has the same numbers of IMFs as , one should havewhere and are the IMF components of and and and are their residuals. Image denoising means that the signal content is removed from . The residual resulting from the BEMD of can be neglected for it is generally very small. Most importantly, after removing all , the basic characteristic information still remains in the obtained image.

The noisy image can be decomposed into some IMFs by the BEMD. These IMFs contain high frequency components in space involved in the original pure image and also include contributions from the IMFs of the noise signal. Accordingly, it is a key problem for the image denoising to remove the latter from the decomposed results of the image. It is known that the BEMD algorithm can adaptively break the original image signal into finite components on different intrinsic scales. That is to say, it should be of benefit to the extraction of information details of the image, to lay a good foundation for the following denoising processing. In general, the noise has major contributions to some IMFs, but quite few to the others. Thus, only by removing those IMFs full of the noise contents can we reconstruct the real image from the original noisy one by synthesizing the remaining IMFs.

##### 4.2. Image Denoising Procedure by the Decomposition of BEMD

In order to obtain reliable decomposition results of the BEMD, the treatment of edge effects is introduced. To some extent, it provides a modified version of the BEMD with edge treatment for enhancing the image denoising performance. The flowchart of its image denoising may be schemed as shown in Figure 1. Details of the procedure are given as follows:(1)Firstly, input the image signal waiting for processing, and set .(2)Expand each edge of the image by means of the proposed approach to get the expansion signal. After that, decompose this signal by the BEMD, which is referred to as the first decomposition. In this decomposition, the noisy image signal is just broken down into one component of IMF and a residual . This IMF is regarded as the component full of noise and edge characteristics involved in the original image. At this time, the parameter is replaced with .(3)Further decompose the obtained IMF by the modified version of the BEMD, which is referred to as the second decomposition. In this decomposition, one can extract the residual, designated as .(4)Synthesize two residual components resulting from above two decomposition operations and obtain the reconstruction image . This image is the denoised result after the th denoising. If ( usually takes 2), take as a new original image and return to step , or otherwise terminate the denoising calculation process.