Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 3151303, 16 pages
http://dx.doi.org/10.1155/2016/3151303

Research Article

Variable Splitting Based Method for Image Restoration with Impulse Plus Gaussian Noise

Tingting Wu

College of Science, Nanjing University of Posts and Telecommunications, Nanjing 210023, China

Received 9 June 2016; Revised 7 October 2016; Accepted 23 October 2016

Academic Editor: Bogdan Dumitrescu

Copyright © 2016 Tingting Wu. 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.

Image denoising is a fundamental problem in realm of image processing. A large amount of literature is dedicated to restoring an image corrupted by a certain type of noise. However, little literature is concentrated on the scenario of mixed noise removal. In this paper, based on the model of two-phase method for image denoising proposed by Cai et al. (2008) and the idea of variable splitting, we are capable of decomposing the image denoising problem into subproblems with closed form. Numerical results illustrate the validity and robustness of the proposed algorithms, especially for restoring the images contaminated by impulse plus Gaussian noise.

1. Introduction

When recording a real-world scene by camera, we desire to acquire an ideal image which is a faithful representation of the scene. However, most observed images are more or less involved in blurry and noise. Hence, deblurring and denoising to the observed image are fundamental aspects in image processing. Let be an ideal image scaled in . Hereafter, we represent an image as () by stacking its columns as a vector. is matrix representation of spatially invariant point spread function (PSF) which characterizes the blurry effects on imaging system and is additive noise with some statistical distributions (e.g., Gaussian noise follows normal distribution, impulse noise follows binomial distribution or uniform distribution, and uniform white noise follows uniform distribution). Accordingly, the observed image involved in the blurring matrix and the additive noise can be boiled down to

The main task of image restoration is to recover the ideal image from the observed image . There are two difficulties in finishing this task. The first one is that the singular values of decay asymptotically to zero which results in a large condition number, and as a consequence (1) is an ill-posed inverse problem in mathematical sense. The second difficulty is that the inevitably additive noise in the observed image makes (1) difficult to handle. These two difficulties make the linear least squares solvers generally fail in deriving ideal image and regularization strategy is therefore introduced to improve the numerical performances, for example, Tikhonov regularization [1], Mumford-Shah regularization [2], and total variation (TV) regularization [3]. The main superiority of TV regularization is that it can regularize images without blurring their edges. The common models utilizing TV regularization for image restoration are the constrained model

and it is equivalent to (or rather equivalent under a proper , see [4]) the unconstrained model

where , are parameters depending on noise variance; is first-order derivative operator and is the TV-norm (details in Section 2); denotes the -norm in Euclidean space (if , herein should be comprehended as ). Commonly, is hinged on the statistical distribution of the additive noise ; for example, for impulse noise removal, for Gaussian noise removal, and for uniform white noise removal. It is therefore pivotal to figure out statistical information of before opting for mathematical model.

Recently, Huang et al. [5] proposed a modified TV regularization model for Gaussian noise removal

where , are positive parameters. Intuitively, a new fitting term is added to the unconstrained model (3). Essentially, by smoothing of the nonsmooth function via the Moreau approximation, the objective function of the following problem

can be replaced by

which is identical with (4). They employed the alternating minimization (AM) algorithm to handle model (4) in variables and , respectively. The numerical results therein indicated that the model was attractive. Lately, they extended model (4) to the mixed noise (impulse plus Gaussian noise) removal in [6]. The extended model reads

where contains the noise-free pixels (pixels are not contaminated by impulse noise); is the characteristic function of set . Specifically, for any , , can be defined as

which is convex and nonsmooth. Then, by using the incomplete data set in , TV minimization method is constructed for image restoration and suitable values can be filled in the detected noisy pixel locations.

As interposition, we expand in this paragraph about mixed noise removal. Gaussian noise [3] generally results from the analog-to-digital conversion of measured voltages and it corrupts an image by exerting subtle perturbations on the gray scale of all pixels. Impulse noise [7] typically arises from the malfunctioning arrays in camera sensors or faulty memory locations in hardware and it affects an image by seriously devastating the gray scales of some pixels (but keeping the other pixels noise-free). Salt-and-pepper noise [8] and random-valued noise [9] are two categories of impulse noise. Let be the dynamic interval of image . For given noise level (also the so-called noise intensity) , salt-and-pepper noise and random-valued noise contaminate image , respectively, by the following:

(i)Salt-and-pepper noise (see [10])

(ii)Random-valued noise (see [11])

 where represents the procedure of contaminating image with impulse noise and satisfies uniform distribution in interval .

Figure 1 shows distinctions of image contaminated by salt-and-pepper noise and random-valued noise. Visually, the salt-and-pepper noise corrupted image possesses black-white speckles while the random-valued noise corrupted image has speckles with various gray scales.

(a)

(b)

(c)

Figure 1: Impulse noise. (a) white-black ideal image. (b) Contaminated by salt-and-pepper noise with . (c) Contaminated by random-valued noise with .

The observed image is contaminated by mixed noise via

where is additive Gaussian noise. It is tricky to retrieve the ideal image in (11) from its observation because the statistical distribution of noise in is superimposed.

Cai et al. [12] proposed a two-phase method for mixed noise removal. They first detect the noise-free set by median-type filter [10, 13] and then execute the procedure of image deblurring on set . The model therein is

where is the set of image edges and the option for is determined experimentally. Model (12) is attractive to render restoration with fine structures in image, such as neat edges and textures. However, the nonconvexity of objective function in model (12) makes it difficult to be minimized. A convex approximation of model (12) is considered therein. However, the computational effort on minimizing that convex approximation is time-consuming. Some literature is devoted to the case of mixed noise removal recently because of its realistic applications, for example, [6, 12, 14–16].

As aforementioned, Huang et al. [6] proposed a fast method for mixed noise removal by solving model (7). The two-phase pattern is exploited as in [12]. Concretely, one has the following.

(i)Phase  1: Noise Detection. Median-type filters are favorable for impulse noise detection. Adaptive median filter (AMF) [10] is adopted for salt-and-pepper noise detection and adaptive center-weighted median filter (ACWMF) [13] is chosen for random-valued noise detection. Assuming that is the observed image and is the restored image by median-type filter, the set of pixels corrupted by impulse noise, denoted by , is detected via

(a)salt-and-pepper noise: .(b)random-valued noise: .

(ii)Phase 2: Minimization on model (7). Model (7) is solved by AM algorithm. Specifically, the -related subproblem is solved by Chambolle’s projection algorithm [17] and the -related subproblem is solved by preconditioned conjugate gradient (PCG) method.

The numerical results in [6] illustrated that their algorithm is efficient and promising. It is about 8–10 times (resp., 15–20 times) faster than the algorithm proposed in [12] for impulse noise (resp., mixed noise removal). However, the properties of model (7) can be further exploited by variable splitting strategy and the resulting subproblems can be tractably tackled with closed-form solutions. This is also the main goal of our paper. Since model (4) is an exceptional case of model (7) with and , we concentrate on handling model (7) in this paper by variable splitting strategy.

The contributions of this paper are as follows. Firstly, we extend models (4) and (7) to deal with different kinds of image denoising, that is, Gaussian noise, salt-and-pepper noise plus Gaussian noise, and random-valued noise plus Gaussian noise in various ways of blur kernels. Secondly, the variable splitting strategies with efficient closed-form solution are investigated. Concretely, we propose the recursions of optimizations (4) and (7) by the classical alternating direction method (ADM). Furthermore, we contemplate coping with (7) with three-block ADM scheme. To guarantee the convergence of algorithm, the ADM with Gaussian-back-substitution is considered for handling model (7). Overall, compared with some standard methods, our proposed models and algorithms exhibit the attractive performance of restoring the images contaminated by impulse plus Gaussian noise.

The rest of the paper is organized as follows. In Section 2, we summarize some basic notations and properties that will be useful in the subsequent sections. In Sections 3 and 4, algorithms based on variable splitting are presented and followed by some numerical experiments in Section 5. The numerical comparisons with the algorithm for Gaussian noise removal in [5] and the algorithm for mixed noise removal in [6] are implemented therein. Finally, we end the paper with some conclusions in Section 6.

2. Preliminaries

In this section, we summarize some basic concepts and properties that will be used in subsequent analysis. Let be a vector and be an integer. By , we denote the -norm of and denotes the Euclidean norm for brevity. The inner product in Euclidean space is denoted . For any positive definite matrix , the matrix norm of is defined by . We denote by the identity operator. For any , we define a vector with the th component . With the definition of -norm, it follows that is the TV-norm of . Let be a function. The domain (resp., epigraph) of is denoted by (resp., ). If is nonempty, the function is said to be proper. The function is called lower semicontinuous (l.s.c.) if is a closed set in , and is a convex function if is convex set in . The subdifferential of , denoted by , is defined

For any , the proximity function of is the well-known soft-thresholding (also the so-called shrinkage). With special denotation , it is defined as

where all operations are calculated in componentwise and is taken as zero if . Let be nonempty set. The indicator function of a nonempty set is denoted by

The proximity operator of , denoted by , is defined as [18, 19]

Specially, the projection operator on a nonempty closed convex set , denoted by , is the proximity operator of .

In TV based image processing, the boundary conditions are indispensable to be taken into account for avoiding unwanted artifacts near the boundary. A common technique is assuming a certain behavior of the sharp image outside the boundary. The conventional boundary conditions are zero, periodic, and reflexive boundary conditions. We consider the reflexive boundary condition herein. By denoting as the first-order derivative operator, the horizontal direction first-order derivative and the vertical direction first-order derivative under reflexive boundary condition are defined as (for notational convenience, is used for this definition)

By denoting as the discrete cosine transform (DCT) and as its inverse, (also the so-called Laplacian operator) under reflective boundary condition can be diagonalized by DCT (see [20]) as

where is a diagonal matrix. Similarly, the blurring matrix (under the circumstance that its PSF is double symmetry) can be diagonalized by DCT (see [21, Chapter 4]) as

where is a diagonal matrix.

3. Alternating Direction Method for Model (7)

Consider the following optimization problem with special structure:

where and are l.s.c. proper convex functions; and are given matrices; is a known vector; and are nonempty closed convex sets (we denote and instead of and in (20) as being merely for convenience of analysis in subsequent sections). Suppose that is an optimal solution of (20). For some , it follows from the first-order optimality condition in optimization that

where is the Lagrange multiplier of (20) [4, 22]. Owing to the special structure of (20), that is, both objective function and linear constraint are separable in variables and , the methods based on variable splitting are recommended, for example, the classical ADM which was proposed originally in [23] and was studied intensively in literature (e.g., [24–28]). Numerous practical applications in diversified realms have illustrated its robustness and efficiency (see [29] for an overview of ADM). By introducing the augmented Lagrangian function of (20) on as

where is the penalty parameter of the linear constraint, the recursions of the well-known ADM can be summarized as Algorithm 1.

Choose arbitrary , and

for     do

               (1)

               (2)

               (3)

Algorithm 1: Alternating direction method.

We now concentrate on dealing with model (7) by ADM, that is, the second phase in [6]. The case in model (7) is considered in this section. By introducing auxiliary variable , we rewrite model (7) as

By denoting ,

the linear constraint in (23) can be reformulated as . Furthermore, by setting , , and

the optimization (23) complies with the pattern of optimization (20) and hence can be handled by ADM.

We elaborate on the recursions for optimization (23) by utilizing ADM as follows. Firstly, the augmented Lagrangian function for (23) can be specified as

(i) The -related subproblem for (23) exploiting ADM amounts to

  which can be approximated efficiently by Chambolle’s projection algorithm [17].(ii) The -related subproblem (involved in the duplets ) for (23) utilizing ADM is

According to the first-order optimality condition and the fact that , it follows that satisfies

Consequently, and can be solved orderly by

We hence only need to concentrate on tackling the first linear system in (30). It can be solved by preconditioned conjugate gradient (PCG) method as suggested in [6] because its coefficient matrix is positive definite.

Remark 1. Specially, if , the model (7) with reduces to the model (4) proposed in [5]. Because of under this circumstance, the -related subproblem is identical to (27) while the -related subproblem decays to

Likewise, PCG method is a better option for the first linear system in (31). Specially, if the PSF of blurring matrix is doubly symmetric, it can be diagonalized by DCT (see [21]). Hence, linear system (31) can be easily solved.

Note that when is adopted in model (7), the resulted recursion for -related subproblem is identical to the case of in (27). However, the -related subproblem is quite different with the case of . The -related subproblem of model (7) with is

which is tricky to be solved because of the nondifferential term (although it can be solved by iterative methods, the computational effort is intensive). It hence inspires our further investigation in following section.

4. Alternating Direction Method with Three Variables for Model (7)

Although the algorithm of ADM in Section 3 is valid for model (7) with and it is favorable to handle model (4) for Gaussian noise removal (i.e., model (7) with and ), we have deadlock over the model (7) with . Consequently, we contemplate coping with model (7) by ADM with three separable variables.

A general extension of optimization (20) is

where , , are lower semicontinuous proper convex functions; , , and are given matrices; is a known vector; , , and are nonempty closed convex sets. Suppose that is an optimal solution of optimization (33). Accordingly, based on the first-order optimality condition in optimization and for some , it follows that

where is the Lagrange multiplier of (33). By defining the augmented Lagrangian function of (33) as

with being the penalty parameter of the linear constraint, the recursions of ADM with three separable variables for (33) are summarized as in Algorithm 2.

Choose arbitrary , , and

for    do

            (1)

            (2)

            (3)

            (4)

Algorithm 2: Alternating direction method with three separable variables.

We now consider utilizing Algorithm 2 to solve the model (7). It will be shown that all the resulting subproblems for model (7) exploiting Algorithm 2 possess closed-form solutions.

By introducing auxiliary variables , , and , the model (7) can be reformulated as

By analogical manipulations as in Section 3, we define ,

Thereby, the linear constraint in (36) can be realigned as in the form of . Furthermore, by denoting , , and ,

the objective function in (36) can be split into three-block functions in variables , , and , respectively. Overall, optimization (36) favors the pattern of (33). Similarly, the augmented Lagrangian function of (36) can be specified as

where is Lagrange multiplier.

We now elucidate the procedures of handling optimization (36) by Algorithm 2, and it will be shown that all the resulting subproblems possess closed-form solutions (either linear least squares problem or problem involved in soft-thresholding).

(i) The -related subproblem for (36) utilizing Algorithm 2 can be formulated as

which essentially is a linear least squares problem. By first-order optimality condition in optimization, it follows that satisfies the linear system

As aforementioned, when the first-order partial derivatives operator is defined as in (17), the above linear system is tractable to be solved by DCT as in (18).

(ii) The -related subproblem for (36) exploiting Algorithm 2 can be rewritten as

which is also a linear least squares problem. Likewise, by the first-order optimality condition, it follows that

which is homogeneous to the first linear system in (31). Accordingly, it can be handled as promoted in Remark 1.

(iii) The -related subproblem for (36) applying Algorithm 2 is involved in the triplets as follows:

It is obvious that the unconstrained optimization (44) is separable in variables , , and . Hence, can be acquired simultaneously (in parallel manner if possible) as follows:

(a)The -related subproblem is

 and its closed-form solution is obtained by the soft-thresholding defined in (14). Concretely, we have .(b) The -related subproblem is

 which can be facilely solved by .(c) The -related subproblem is

 Combining the definition of , the following follows. If ,

 If ,

Remark 2. Literally, we take a fixed penalty for the foregoing inferences. In practice, it is attractive to choose dynamically, for example, the continuation strategy (e.g., [30]) and self-adaptive strategy (e.g., [28]). Moreover, it is reasonable to choose differently in the last three terms of (39) so as to penalize those three different linear constraints in (36); that is, the augmented Lagrangian function in (39) can be rewritten as

where , , are penalties corresponding to different linear constraints. Hence, the recursions for -, -, and -related subproblems can be regained accordingly. Generally, those skills on penalties can numerically render impressive performance (see, e.g., [7, 31, 32] and references therein).

Recently, A series variant algorithms for the ADM with more than two variables [33–35] are proposed, which are proved to be convergent and illustrate attractive numerical performances in tackling optimization problems with separable structure. Generally, the three-block ADM may not converge, which was demonstrated by Chen et al. [36] using counterexamples. Nevertheless, if all the objective functions are strongly convex, Han and Yuan [37] proved the global convergence of the three-block ADM scheme. Then, this condition was relaxed and only two or more functions in the objective are required to be strongly convex to ensure the convergence. More recently, it was studied in [38, 39] that the global convergence of the three-block ADM scheme can be guaranteed only when one of the objective functions is strongly convex. However, neither of the three-block functions in the objective function in (36) is strongly convex and then the convergence of Algorithm 2 is still ambiguous. Under some special circumstances (e.g., the matrix involved in the linear constraint is the identity matrix or it is sparse), the performances of some variant algorithms can be matchable with that of Algorithm 2, even fairly close to that of Algorithm 2 [36]. We recommend herein the ADM with Gaussian-back-substitution algorithm proposed in [40] for handling model (7).

By denoting variable and defining the sparse matrix as

where is blurring matrix and is defined in (37), ADM with Gaussian-back-substitution algorithm for the concrete optimization (36) is summarized in Algorithm 3.

Choose arbitrary , , ,

and

for    do

            (1)

            (2)

            (3)

            (4)

            (5)

Algorithm 3: Alternating direction method with Gaussian-back-substitution for (36).

The augmented Lagrangian function in Algorithm 3 is as defined in (39). The detailed recursions for solving (36) utilizing Algorithm 3 can be inferred analogically as the scenario of using Algorithm 2. The convergence analysis of Algorithm 3 was expanded in [40]. The option on parameter is recommended for practical numerical experiments.

5. Numerical Experiments

We test the aforementioned algorithms on the fundamental image processing problem, image denoising (both Gaussian noise removal and mixed noise removal). We first concentrate on the scenario of Gaussian noise removal in Section 5.1, with numerical comparisons to the algorithm in [5] (abbr, HNW08) and then turn to the case of mixed noise removal in Section 5.2, with numerical comparisons to the algorithm in [6] (abbr, HNW09). We abbreviate hereafter Algorithms 1, 2, and 3 as “ADM2,” “ADM3,” and “ADMGB,” respectively. All the codes are written in MATLAB 7.9 and are run on a desktop computer with Intel(R) Core (TM) 2.66 HZ and 2 G memory. All the testing images for numerical experiments in this section are illustrated in Figure 2. To measure the quality of restored image, we define the signal-to-noise ratio (SNR) as

where is the restoration of the ideal image . The stopping criterion for testing algorithms is set as

As mentioned in [12], it is an old and open problem for choosing the optimal parameters in algorithms; we therefore tone those parameters (e.g., , , and ) by trial and error. As suggested in Remark 2, we opt for the penalty differently in Algorithms 2 and 3 for various linear constraints.

Figure 2: Ideal images: shape.jpg, cameraman.png, Barbara.png, man.png, boat.png, and hill.png.

5.1. Gaussian Noise Removal

In this subsection, we consider the case of Gaussian noise removal, that is, the model (4). The ideal images are degraded by various blurs and are further contaminated by additive Gaussian noise with noise variance . Those procedures can be realized via fspecial and imnoise in MATLAB Image Processing Toolbox (IPT). The concrete parameters for both commands in experiments are summarized in Table 1.

Table 1: Problems for Gaussian noise removal.

Blur
Uniform
Gaussian
Out-of-focus

Shape
Hsize = 7
hsize = 7, sigma = 2
radius = 3

Cameraman
Hsize = 9
hsize = 9, sigma = 2
radius = 7

Barbara
Hsize = 11
hsize = 11, sigma = 5
radius = 9

We test HWN08 (codes are available at http://www.math.hkbu.edu.hk/~mng/) and ADM2 in this numerical experiments. We choose the weights and in model (4) for both testing algorithms. The penalty is adopted for ADM2. Both HNW08 and ADM2 require Chambolle’s projection algorithm [17] to approximate the solution of -related subproblem. We take 10 steps for this approximation at each iteration. The initial point is chosen as the observed degraded image and and . The stopping criterion in (53) is set as . The numerical results are reported in Table 2. The datum therein, such as “,” represents the iteration number, computing time in seconds and SNR, respectively. We can see that ADM is faster than AM algorithm. The reason is that the dual variable is exploited and is updated at each iteration. Figure 3 illustrated some devastated images and their restorations by ADM2.

Table 2: Computational results on Gaussian noise removal by solving model (4).

Blur

Algorithm
Shape
Cameraman
Barbara

Uniform
0.005
HNW08
62/1.94/23.91
61/7.33/21.77
51/31.63/18.23

ADM2
24/0.75/24.00
26/3.59/21.68
26/17.38/18.19

0.01
HNW08
53/1.67/20.87
51/6.72/19.91
52/32.61/17.70

ADM2
30/0.92/21.69
29/4.05/20.12
31/20.41/17.68

Gaussian
0.005
HNW08
82/2.56/20.89
80/9.70/20.28
58/31.22/18.02

ADM2
35/1.05/21.22
34/4.70/20.39
30/19.88/18.01

0.01
HNW08
58/1.80/17.93
54/6.48/19.53
57/35.38/17.64

ADM2
35/1.08/19.12
36/4.92/19.74
34/22.55/17.61

Out-of-focus
0.005
HNW08
58/1.83/23.37
66/7.91/20.43
61/32.91/18.03

ADM2
20/0.77/23.43
29/3.98/20.43
33/22.11/18.01

0.01
HNW08
48/1.23/19.34
55/6.59/18.81
57/31.56/17.39

ADM2
23/0.80/20.38
34/4.61/19.09
34/22.25/17.46

(a)

(b)

Figure 3: (a) Image devastated by out-of-focus blur and Gaussian noise (). (b) Restored images by ADM2.

5.2. Mixed Noise Removal

We now test the case of mixed noise removal, that is, the model (7), and report numerical comparisons between HNW09, ADM2, ADM3, and ADMGB. The ideal image of cameraman is degraded by Gaussian blur with hsize = 7 andsigma = 2 and out-of-focus blur with radius = 3, respectively. The blurred images are further corrupted by additive Gaussian noise with noise variance as and impulse noise with various noise level . The median-type filters are attractive in removing impulse noise. We choose AMF for salt-and-pepper noise detection while ACWMF for random-valued noise detection as in [6].

We first experiment on the scenario of salt-and-pepper noise plus Gaussian noise removal. The parameters in model (7) are taken as , , and for all testing algorithms. Moreover, we choose for ADM2; , , and for ADM3 and ADMGB. The window size for median-type filter AMF is set as 19 as in [12]. Because the - and -related subproblems in both HNW09 and ADM2 should be approximated by Chambolle’s projection method [17] and PCG method, respectively, we set 5-step Chambolle’s approximation for -related subproblem while as tolerance for pcg in MATLAB for -related subproblem. The initial points in experiments are taken , , and for all testing algorithms and the stopping criterion in (53) is set as . The numerical results are listed in Table 3. We plot the evolutions of SNR with respect to iteration numbers and computing time in Figure 4. It manifests the efficiency of ADM3 and ADMGB as the resulting subproblems possess closed-form solutions. Figure 5 illustrates some restored images by ADM3 with different impulse noise level. We also test algorithms on other images for salt-and-pepper noise plus Gaussian noise removal. The ideal images are degraded by out-of-focus blur with radius = 7 and contaminated by salt-and-pepper noise with noise level plus Gaussian noise with noise variance . Figure 6 shows the restored images by ADM3.

Table 3: Computational results on salt-and-pepper noise plus Gaussian noise removal by solving model (7).

Blur

HNW09
ADM2
ADM3
ADMGB

Gaussian
30%
23/39.08/20.65
14/30.72/20.68
45/6.53/20.95
46/10.08/20.92

50%
29/49.66/20.41
15/32.61/20.42
50/7.66/20.57
51/11.28/20.59

70%
36/61.31/20.09
17/36.58/20.10
66/9.95/20.12
67/13.55/20.13

90%
56/89.73/18.74
30/64.38/18.71
143/21.48/18.77
154/28.22/18.62

Out-of-focus
30%
19/32.78/22.66
13/27.95/22.30
46/6.67/22.59
47/10.13/22.48

50%
24/40.27/21.95
14/30.08/21.72
53/7.91/21.96
54/11.20/21.92

70%
35/59.88/20.89
17/37.61/20.74
85/11.89/20.82
94/19.78/20.89

90%
57/97.95/18.43
33/71.48/18.57
144/21.27/18.58
146/30.77/18.42

Figure 4: Evolutions of SNR with respect to iteration numbers and computing time (cameraman image degraded by Gaussian blur and contaminated by salt-and-pepper noise () plus Gaussian noise ()).

(a)

(b)

(c)

Figure 5: (a) Images degraded by Gaussian blur and contaminated by salt-and-pepper noise (, 50%, 70%, and 90%) plus Gaussian noise (). (b) Filtered images by AMF. (c) Restored images by ADM3.

(a)

(b)

Figure 6: (a) Devastated images by out-of-focus blur (radius = 7) and salt-and-pepper noise () plus Gaussian noise (). (b) Restored images by ADM3. From left to right: Barbara (96/80.39/17.77), man (95/80.13/20.98), boat (95/79.80/21.04), and hill (93/78.55/21.26).

We now experiment on the scenario of random-valued noise plus Gaussian noise removal. Generally, Gaussian noise mixed with random-valued noise is more intractable to be removed than that of Gaussian noise mixed with salt-and-pepper noise. The reason is that it is difficult to distinguish pixels contaminated by Gaussian noise with those contaminated by small value random-valued noise. We tackle model (7) by choosing , , and for all testing algorithms. We take for ADM2 and , , and for ADM3 and ADMGB. The median-type filter ACWMF is run 4 times in phase of noise detection as in [12]. The - and -related subproblems in HNW09 and ADM2 are still approximated by Chambolle’s projection method and PCG method as the scenario of random-valued noise plus Gaussian noise removal. Likewise, the initial points are also taken as , , and and the stopping criterion in (53) is set as . The numerical results are listed in Table 4. Figure 7 plots the evolutions of SNR with respect to iteration numbers and computing time. Some restored images by ADM3 are illustrated in Figure 8. Those images are also tested for random-valued noise plus Gaussian noise removal. All images are degraded by out-of-focus blur with radius = 7 and contaminated by random-valued noise with noise level and Gaussian noise with noise variance . Figure 9 shows the restored images by ADM3.

Table 4: Computational results on random-valued noise plus Gaussian noise removal by solving model (7).

Blur

HNW09
ADM2
ADM3
ADMGB

Gaussian
20%
30/40.95/19.08
24/43.78/19.29
42/6.28/19.28
43/8.77/19.25

30%
31/41.80/18.91
24/41.80/19.01
43/6.28/19.14
45/9.64/19.12

40%
31/39.66/18.50
25/42.42/18.61
43/6.33/18.54
45/9.17/18.77

50%
33/41.25/17.40
27/47.44/17.22
47/6.89/17.28
49/10.03/17.26

Out-of-focus
20%
27/35.69/19.46
21/38.78/19.77
38/5.78/19.73
40/8.33/19.73

30%
27/36.95/19.20
22/40.64/19.45
40/5.91/19.40
42/8.80/19.50

40%
30/40.58/18.48
23/42.22/18.95
41/6.19/18.93
45/9.30/18.66

50%
39/48.56/17.05
25/44.56/17.15
45/6.63/17.67
49/10.19/17.20

Figure 7: Evolutions of SNR with respect to iteration numbers and computing time (cameraman image degraded by out-of-focus blur and contaminated by random-valued noise () plus Gaussian noise ()).

(a)

(b)

(c)

Figure 8: (a) Devastated images by out-of-focus blur (radius = 7) and random-valued noise (, 30%, 40%, and 50%) plus Gaussian noise (). (c) Restored images by ADM3.

(a)

(b)

Figure 9: (a) Devastated images by out-of-focus blur (radius = 7) and random-valued noise () plus Gaussian noise (). (b) Restored images by ADM3. From left to right: Barbara (42/35.17/17.17), man (42/34.89/19.19), boat (41/33.92/18.90), and hill (42/35.13/19.66).

6. Conclusions

The variable splitting is a hot-investigated strategy in both variational inequalities and optimization. Numerous practical applications illustrate the attractive and efficient performance of this strategy. We illustrate in this paper that the variable splitting can be efficient and promising in image restoration. Note that the numerical effects for mixed noise removal utilizing model (7) with are similar as that of (see also [6]). As stated in this paper, we only consider the image denoising with blurry. However, the mixed noise may be occurring in image inpainting, superresolution, and decomposition, even in the area of video processing, which may be difficult to be handled. They are thus some of our future research topics.

Competing Interests

The author declares that there are no competing interests.

Acknowledgments

The author was supported by the Natural Science Foundation of China (Grants nos. 11501301 and 11401319), Jiangsu Planned Projects for Postdoctoral Research Funds (Grant no. 1501071B), and the Foundation of Jiangsu Key Lab for NSLSCS (Grant no. 201601).

References

A. Tikhonov and V. Arsenin, Solution of Ill-Posed Problems, John Wiley & Sons, New York, NY, USA, 1977.

D. Mumford and J. Shah, “Optimal approximations by piecewise smooth functions and associated variational problems,” Communications on Pure and Applied Mathematics, vol. 42, no. 5, pp. 577–685, 1989.

L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D, vol. 60, no. 1–4, pp. 259–268, 1992.

R. T. Rockafellar, Convex Analysis, Princeton University Press, Princeton, NJ, USA, 1970.

Y. M. Huang, M. K. Ng, and Y.-W. Wen, “A fast total variation minimization method for image restoration,” Multiscale Modeling & Simulation, vol. 7, no. 2, pp. 774–795, 2008.

Y.-M. Huang, M. K. Ng, and Y.-W. Wen, “Fast image restoration methods for impulse and Gaussian noises removal,” IEEE Signal Processing Letters, vol. 16, no. 6, pp. 457–460, 2009.

D. R. Han, X. M. Yuan, and W. X. Zhang, “An augmented Lagrangian based parallel splitting method for separable convex minimization with applications to image processing,” Mathematics of Computation, vol. 83, no. 289, pp. 2263–2291, 2014.

R. H. Chan, C.-W. Ho, and M. Nikolova, “Salt-and-pepper noise removal by median-type noise detectors and detail-preserving regularization,” IEEE Transactions on Image Processing, vol. 14, no. 10, pp. 1479–1485, 2005.

Y. Q. Dong, R. H. Chan, and S. F. Xu, “A detection statistic for random-valued impulse noise,” IEEE Transactions on Image Processing, vol. 16, no. 4, pp. 1112–1120, 2007.

H. Hwang and R. A. Haddad, “Adaptive median filters: new algorithms and results,” IEEE Transactions on Image Processing, vol. 4, no. 4, pp. 499–502, 1995.

T. Chen and H. R. Wu, “Space variant median filters for the restoration of impulse noise corrupted images,” IEEE Transactions on Circuits and Systems II: Analog and Digital Signal Processing, vol. 48, no. 8, pp. 784–789, 2001.

J.-F. Cai, R. H. Chan, and M. Nikolova, “Two-phase approach for deblurring images corrupted by impulse plus Gaussian noise,” Inverse Problems and Imaging, vol. 2, no. 2, pp. 187–204, 2008.

S.-J. Ko and Y. H. Lee, “Center weighted median filters and their applications to image enhancement,” IEEE Transactions on Circuits and Systems, vol. 38, no. 9, pp. 984–993, 1991.

C. Cocianu and A. Stan, “Neural architectures for correlated noise removal in image processing,” Mathematical Problems in Engineering, vol. 2016, Article ID 6153749, 14 pages, 2016.

S. Morillas, V. Gregori, and A. Hervás, “Fuzzy peer groups for reducing mixed Gaussian-impulse noise from color images,” IEEE Transactions on Image Processing, vol. 18, no. 7, pp. 1452–1466, 2009.

M.-G. Shama, T.-Z. Huang, J. Liu, and S. Wang, “A convex total generalized variation regularized model for multiplicative noise and blur removal,” Applied Mathematics and Computation, vol. 276, pp. 109–121, 2016.

A. Chambolle, “An algorithm for total variation minimization and applications,” Journal of Mathematical Imaging and Vision, vol. 20, no. 1-2, pp. 89–97, 2004.

P. L. Combettes and V. R. Wajs, “Signal recovery by proximal forward-backward splitting,” Multiscale Modeling & Simulation, vol. 4, no. 4, pp. 1168–1200, 2005.

J.-J. Moreau, “Fonctions convexes duales et points proximaux dans un espace hilbertien,” Comptes Rendus de l'Académie des Sciences, vol. 255, pp. 2897–2899, 1962.

M. J. Buckley, “Fast computation of a discretized thin-plate smoothing spline for image data,” Biometrika, vol. 81, no. 2, pp. 247–258, 1994.

P. C. Hansen, J. G. Nagy, and D. P. O'Leary, Deblurring Images: Matrices, Spectra, and Filtering, vol. 3 of Fundamentals of Algorithms, SIAM, Philadelphia, Pa, USA, 2006.

D. P. Bertsekas, Constrained Optimization and Lagrange Multiplier Methods, Academic Press, New York, NY, USA, 1982.

D. Gabay and B. Mercier, “A dual algorithm for the solution of nonlinear variational problems via finite element approximation,” Computers & Mathematics with Applications, vol. 2, no. 1, pp. 17–40, 1976.

M. Fukushima, “Application of the alternating direction method of multipliers to separable convex programming problems,” Computational Optimization and Applications, vol. 1, no. 1, pp. 93–111, 1992.

R. Glowinski and P. Le Tallec, Augmented Lagrangian and Operator-Splitting Methods in Nonlinear Mechanics, vol. 9 of SIAM Studies in Applied Mathematics, SIAM, Philadelphia, Pa, USA, 1989.

D. R. Han, H. J. He, H. Yang, and X. M. Yuan, “A customized Douglas-Rachford splitting algorithm for separable convex minimization with linear constraints,” Numerische Mathematik, vol. 127, no. 1, pp. 167–200, 2014.

D. R. Han and X. M. Yuan, “Local linear convergence of the alternating direction method of multipliers for quadratic programs,” SIAM Journal on Numerical Analysis, vol. 51, no. 6, pp. 3446–3457, 2013.

B. S. He, L.-Z. Liao, D. R. Han, and H. Yang, “A new inexact alternating directions method for monotone variational inequalities,” Mathematical Programming, vol. 92, no. 1, pp. 103–118, 2002.

S. Boyd, N. Parikh, E. Chu, B. Peleato, and J. Eckstein, “Distributed optimization and statistical learning via the alternating direction method of multipliers,” Foundations and Trends® in Machine Learning, vol. 3, no. 1, pp. 1–122, 2011.

E. T. Hale, W. Yin, and Y. Zhang, “Fixed-point continuation for -minimization: methodology and convergence,” SIAM Journal on Optimization, vol. 19, no. 3, pp. 1107–1130, 2008.

M. K. Ng, P. Weiss, and X. Yuan, “Solving constrained total-variation image restoration and reconstruction problems via alternating direction methods,” SIAM Journal on Scientific Computing, vol. 32, no. 5, pp. 2710–2736, 2010.

Z. N. Zhu, G. Ch. Cai, and Y.-W. Wen, “Adaptive box-constrained total variation image restoration using iterative regularization parameter adjustment method,” International Journal of Pattern Recognition and Artificial Intelligence, vol. 29, no. 7, Article ID 1554003, 20 pages, 2015.

C. Chen, R. H. Chan, S. Ma, and J. Yang, “Inertial proximal ADMM for linearly constrained separable convex optimization,” SIAM Journal on Imaging Sciences, vol. 8, no. 4, pp. 2239–2267, 2015.

C. H. Chen, M. Li, X. Liu, and Y. Y. Ye, “On the convergence of multi-block alternating direction method of multipliers and block coordinate descent method,” https://arxiv.org/abs/1508.00193

W. H. Yang and D. R. Han, “Linear convergence of the alternating direction method of multipliers for a class of convex optimization problems,” SIAM Journal on Numerical Analysis, vol. 54, no. 2, pp. 625–640, 2016.

C. H. Chen, B. S. He, Y. Y. Ye, and X. M. Yuan, “The direct extension of ADMM for multi-block convex minimization problems is not necessarily convergent,” Mathematical Programming, vol. 155, no. 1-2, pp. 57–79, 2016.

D. R. Han and X. M. Yuan, “A note on the alternating direction method of multipliers,” Journal of Optimization Theory and Applications, vol. 155, no. 1, pp. 227–238, 2012.

X. J. Cai, D. R. Han, and X. M. Yuan, “On the convergence of the direct extension of ADMM for three-block separable convex minimization models with one strongly convex function,” Computational Optimization and Applications, 2016.

M. Li, D. F. Sun, and K.-C. Toh, “A convergent 3-block semi-proximal ADMM for convex minimization problems with one strongly convex block,” Asia-Pacific Journal of Operational Research, vol. 32, no. 4, Article ID 1550024, 19 pages, 2015.

B. He, M. Tao, and X. Yuan, “Alternating direction method with Gaussian back substitution for separable convex programming,” SIAM Journal on Optimization, vol. 22, no. 2, pp. 313–340, 2012.