Abstract
The question of eliminating the multiplicative noise has attracted much interest in many research studies. In this work, we are interested with the one that are follows the Gamma distribution. The rationale of this paper is to shed light on a brief comparative study of some local and nonlocal models, for denoising images contaminated with the noise of this type. The improved method of Split Bregman is used to implement those models. The totality of experiments indicates that the proposed nonlocal method gives better results than some other methods.
This work was supported by the PPR CNRST: Modèles Mathématiques appliqués à l’environement, à l’imagerie médicale et aux Biosystèmes.
You have full access to this open access chapter, Download conference paper PDF
Similar content being viewed by others
Keywords
1 Introduction
The main goal of an image denoising algorithm is to gain access though both noise suppression and the maintenance of important features such as contours and textures. In many real world applications, images are not usually contaminated with additive noise. This can be expressed in a word the presence of a noise called multiplicative, which appears in various image processing applications, for instance in Synthetic Aperture Radar (SAR) and Ultrasound imaging [10, 15]. For a mathematical description of such a problem we suppose that \(\varOmega \) is a connected open subset of . Let be the original image, from the degraded mechanism of multiplicative noise, the observation image f is obtained by
where \(\eta \) is multiplicative noise. In this work, we are interested on the multiplicative noise that follows the Gamma distribution. So far, a many variational models have been considered to manipulate the image restoration problem with the multiplicative noise, while a large amount of works on this subject are total variation (TV)-based. The first approach with the TV regularization devoted to multiplicative noise removal (suited for Gaussian multiplicative noise) also have proposed by Rudin et al. [11] (RLO-model). By using property of the Maximum a Posteriori (MAP) estimator and TV as edge preserver, Aubert and Aujol [2] have introduced a multiplicative denoising variational model (AA-model) based on the Bayes rule and Gamma distribution, whose energy functional can be written as follows
where the first term is the regularization term which imposes some prior constraints on the original image, and the second term is the fitting term which measures the violation of relation between u and the observation image f. Obviously, because of the non-convexity of (2), it does not have the global optimal solution. To make this problem treatable, many approaches [3, 4, 9, 13, 14] have been proposed in the literature. Among those approaches, Huang et al. [9] considered an exponential-transformation \(u\longrightarrow \exp (u)\) in the fitting term of (2) to obtain the following strictly convex model
where \(\lambda _{1}\) and \(\lambda _{2}\) are positive regularization parameters, and \(z=\log u\) is an auxiliary variable. In a recent work, Huang et al. [8] proposed a significant modification of the regular term in (2). In deed, they provided a non-convex Bayesian type model for multiplicative noise removal, which includes TV and the Weberized TV as regularizers [12]. This gives an effective improvement over the previous models in terms of visual appearance. This model can be described as follows
where the first two terms are the regularization terms, while the third one is the nonconvex data fitting term of (2). \(\alpha _{1},\alpha _{2}\) are regularization parameters. The main idea of the second regular term in (4) is to discover the universal influence of ambient intensity level u on human’s sensitivity to the local intensity increment \(|\nabla u|\), or so called Just Noticeable Difference (JND) [16], in the perception of both sound and light. This local fluctuation can be formulated by
according to Weber’s law [12, 16], when the mean intensity of the background is increasing with a higher value, then the intensity increments \(|\nabla u|\) also has higher value. In the other hand, Dong et al. [4] proposed to replace the classical TV-norm in [9, 13] with the discrete nonlocal norm one and they applied the split Bregman iteration to solve this two nonlocal models, by introducing the two constraints \(z=\log u\) and \(z_{0}=\log f\), this two models can be described as follows
Inspired by the idea of [8], we propose a nonlocal variational model for multiplicative noise removal
The rest of this paper proceeds as follows. In Sect. 2, first we describe the necessary definitions and tools of some nonlocal operators. After, we give a detailed implementation of the efficient computational method for obtaining the numerical solution of (7) using an improved split Bregman algorithm. Numerical experiments intended for the effectiveness of the proposed method are provided in Sect. 3. Finally, conclusions are made in Sect. 4.
2 Our Proposed Model
2.1 Nonlocal Differential Operators
First, we give a brief overview of the definitions of some nonlocal functionalities (for more details see [1, 5, 6]). Let , be a real function and \(\mu \) be a probability measure. Let a nonnegative symmetric function represents the weights between the current pixel and the pixels in \(\varOmega \). Given a pair of points \((x,y)\in \varOmega \times \varOmega \), then the nonlocal gradient is defined as
The inner product between two vectors and the associated norm of \(q_{1}\) at \(x\in \varOmega \) are defined as
Hence, the norm of \(\nabla _{NL}u\) at \(x\in \varOmega \) is referred to as
The nonlocal divergence operator can be defined by the standard adjoint relation with the nonlocal gradient, , as follows
which leads to the definition of nonlocal divergence of the vector q
The general nonlocal p-Laplacian operator for \(1\le p <\infty \), which is defined by
2.2 The Proposed Model
In this subsection, we describe a method to solve our model (7). We first take \(\psi (u)=\alpha _{1}+\alpha _{2}/u\), the minimization problem (7) can be rewritten as
Generally, the nonlocal TV norm in (14) is difficult to be computed straightforwardly. To avoid this drawback, we adopt the split Bregman iteration, which was initially introduced by Goldstein and Osher [7]. This method is originally designed for the \(L_{1}\) regularization \(+ L_{2}\) minimization problem. Note that the second term contains a log-term, which will not allow to apply directly this method. To overtake this difficulty, we proceed as in [4] and introduce an auxiliary variable z, such that \(z=u\) and go on to solve the following unconstrained optimization problem
where \(\lambda >0\) should be large enough so that z is sufficiently close to u in the sense of \(L_{2}\)-norm. Subsequently, the minimization problem (14) can be iteratively solved by minimizing the following two subproblems
The first subproblem (16) is a \(L_{1}\) regularization \(+L_{2}\) minimization problem, thus it can be efficiently solved by the classical split Bregman iteration. Let \(d=\nabla _{NL}u\), then (16) becomes
where \(\beta \) is a positive regularization parameter. To solve the problem (15) is equivalent to solve the following system
Taking in the system (18), then (18) can be written as
The minimizers of first and last subproblems from system (18) are characterized by the optimality condition given by the Euler-Lagrange formulation. The second subproblem from system (18) can be implemented via the generalized shrinkage formula
Which are outlined in the following system
Using the Gauss-Seidel iterative scheme, then the numerical algorithm summarizes all of these elements
Let \(\psi ^{'}(u_{j}^{k+1,n}) = -\frac{\alpha _{2}}{(u_{j}^{k+1,n})^{2}}\), \(|d^{k+1,n}| = \sqrt{\sum _{j}(d_{i,j}^{k+1,n})^{2}+(d_{j,i}^{k+1,n})^{2}}\) and
\(d^{k+1,n=0}=d^{k}\), \(b^{k+1,n=0}=b^{k}\) and \(z^{k+1,l=0}=z^{k}\). We choose the weight function as
where \(d(x,y)=\int _{\varOmega }G_{\varrho }(s)|u(x+s)-u(y+s)|^{2}ds\) is the distance between patches located at x and y, \(G_{\varrho }\) is a Gaussian kernel of standard deviation \(\varrho \) and \(h_{0}\) is a filtering parameter. In all experiments, we fixed the number of outer iterations in the Split Bregman to be four (i.e., \(k=4\)) and number of inner iterations to be two (i.e., \(n=2\)). We use patches of \(5\times 5\) to compute the weight and a search window of size \(11\times 11\).
3 Numerical Results
In this section, we present numerical results so as to illustrate the performance of our proposed model. These results are compared to those obtained by the AA-model proposed by Aubert et al. (2), the HNW-model proposed by Huang et al. (3) and the HXW-model proposed by Huang et al. (4). Furthermore, we compared our model with the second nonlocal DZK-model proposed by Dong et al. (6). It’s worth anticipating that all the numerical simulations (exclude DZK-model) are implemented by the approach described in the previous Sect. 2.2. In all tests, each pixel of an original images (Figs. 1(a)–3(a)) are degraded by a noise which follows a Gamma distribution. The noise level is controlled by the value of \(\sigma \) in the experiments, the noisy images with different levels (\(\sigma =0.08,0.1,0.2\)) are shown in Figs. 1(b)–3(b). For measuring the quality of the restored image, three tools are considered. The first one is the structural similarity index (SSIM), the second one is the peak signal to noise ratio (PSNR), and the last one is the signal to noise ratio (SNR), thus these measures are given by
where \(u^{*}\), u, and \(M\times N\) are respectively the restored image, the true image and the size of image. \(\bar{u}\) and \(\bar{u}^{*}\) represent the mean values of u and \(u^{*}\) in image domain \(\varOmega \). \(\sigma _{u}\) and \(\sigma _{u^{*}}\) denote the variances of u and \(u^{*}\). \(\sigma _{uu^{*}}\) is the covariance of u and \(u^{*}\). \(C_{1}\) and \(C_{2}\) are two variables to stabilize the division with a low denominator. Figures 1(c)–3(f) show the denoising results of the three noisy gray level images by different methods. In these experiments, it is clear that the restoration results obtained by the proposed method are visually better than those by the AA, HNW and HXW-models. Table 1 summarize the PSNRs and SSIMs results by different methods on the three noisy gray level images. From Table 1, we can see that PSNRs and SSIMs values of restored images using our method are wider than those restored by using the other three methods. In Table 2, we present the PSNRs and SSIMs results obtained by some specific cases of general model (7) on the three noisy gray level images. This can be included the models:
-
When \(\alpha _{1}=0\), this reduces to the nonlocal version of the model in [12] with the fitting term of (2).
-
When \(\alpha _{2}=0\), this reduces to the nonlocal version of the AA-model (2). We choose these tests to show that each model could be the best one among the three models. All the three nonlocal models obtain quite good results, than the other local models. From Table 2, It turns out that the results of the images restored, when noise level is small (\(\sigma =0.08,0.1\)) are better than those of other cases. But, when we increase \(\sigma \) (\(\sigma =0.2\)), we obtained almost the same results. Table 3 shows the SNRs results of our model compared with the nonlocal DZK-model on the two noisy gray level Barbara and Lena images under noise level \(\sigma =0.1\). From Table 3, we can see that our method gives better results than those obtained by DZK-model (6). The main reason for getting better results is that the nonlocal combination of TV and Weberized TV regularizers are dominated in our model. It can successfully remove the multiplicative noise following Gamma distribution, while preserving fine structures and textures as well.
4 Conclusion
In this work, we have proposed a general nonlocal model for multiplicative noise removal problem that following Gamma distribution under the combination of two nonlocal regularizers the weberized TV and the TV operator. By combining these two regularizers, our nonlocal model outperform the classical AA-, HNW-, HXW- and DZK-models. It can both preserve more fine structures and remove the noise of this type. We have tested all of those models under different noise levels and their performances are evaluated and compared both visually and quantitatively.
References
Andreu, F., Mazón, J.M., Rossi, J.D., Toledo, J.: A nonlocal \(p\)-Laplacian evolution equation with nonhomogeneous Dirichlet boundary conditions. SIAM J. Math. Anal. 40(5), 1815–1851 (2008)
Aubert, G., Aujol, J.F.: A variational approach to removing multiplicative noise. SIAM J. Appl. Math. 68(4), 925–946 (2008)
Bioucas-Dias, J.M., Figueiredo, M.A.T.: Total variation restoration of speckled images using a split-Bregman algorithm. In: Image Proceedings of IEEE (ICIP) (2009)
Dong, F., Zhang, H., Kong, D.X.: Nonlocal total variation models for multiplicative noise removal using split Bregman iteration. J. Math. Comput. Model. 55(3–4), 939–954 (2012)
Elmoataz, A., Lezoray, O., Bougleux, S.: Nonlocal discrete regularization on weighted graphs: a framework for image and manifold processing. IEEE Trans. Image Process. 17(7), 1047–1060 (2008)
Gilboa, G., Osher, S.: Nonlocal operators with application to image processing. Multiscale Model. Simul. 7(3), 1005–1028 (2008)
Goldstein, T., Osher, S.: The split Bregman method for \(L1\) regularized problems. SIAM J. Imaging Sci. 2(2), 323–343 (2009)
Huang, L.L., Xiao, L., Wei, Z.H.: Multiplicative noise removal via a novel variational model. J. Image Video Process. (2010). Article ID 250768
Huang, Y.M., Ng, M.K., Wen, Y.W.: A new total variation method for multiplicative noise removal. SIAM J. Imaging Sci. 2(1), 20–40 (2009)
Ramamoorthy, S., Siva Subramanian, R., Gandhi, D.: An efficient method for speckle reduction in ultrasound liver images for e-health applications. In: Natarajan, R. (ed.) ICDCIT 2014. LNCS, vol. 8337, pp. 311–321. Springer, Cham (2014). https://doi.org/10.1007/978-3-319-04483-5_32
Rudin, L., Lions, P.L., Osher, S.: Multiplicative denoising and deblurring: theory and algorithms. In: Osher, S., Paragios, N. (eds.) Geometric Level Set Methods in Imaging, Vision, and Graphics, pp. 103–119. Springer, New York (2003). https://doi.org/10.1007/0-387-21810-6_6
Shen, J.: On the foundations of vision modeling: I. Weber’s law and Weberized TV restoration. J. Phys. D: Nonlinear Phenom. 175(3–4), 241–251 (2003)
Shi, J., Osher, S.: A nonlinear inverse scale space method for a convex multiplicative noise model. SIAM J. Imaging Sci. 1(3), 294–321 (2008)
Steidl, G., Teuber, T.: Removing multiplicative noise by Douglas-Rachford splitting methods. J. Math Imaging Vis. 36(2), 168–184 (2010)
Sveinsson, J.R., Benediktsson, J.A.: Speckle reduction and enhancement of SAR images in the wavelet domain. In: Proceedings of the International Geoscience and Remote Sensing Symposium, vol. 1, pp. 63–66 (1996)
Weber, E.H.: De pulsu, resorptione, auditu et tactu, in Annotationes anatomicae et physio-logicae. Koehler, Leipzig (1834)
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2018 Springer International Publishing AG, part of Springer Nature
About this paper
Cite this paper
Karami, F., Meskine, D., Oubbih, O. (2018). A Nonlocal Model for Image Restoration with Gamma Distributed Multiplicative Noise. In: Mansouri, A., El Moataz, A., Nouboud, F., Mammass, D. (eds) Image and Signal Processing. ICISP 2018. Lecture Notes in Computer Science(), vol 10884. Springer, Cham. https://doi.org/10.1007/978-3-319-94211-7_42
Download citation
DOI: https://doi.org/10.1007/978-3-319-94211-7_42
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-319-94210-0
Online ISBN: 978-3-319-94211-7
eBook Packages: Computer ScienceComputer Science (R0)