1 February 2011 Adaptive mixed image denoising based on image decomposition
Author Affiliations +
Abstract
Image denoising while preserving image features is a key problem in image processing and computer vision. This letter proposes an adaptive mixed method for image restoration. First, this method decomposes a given image as the sum of two components: geometric structure and oscillating pattern according to Meyer's theory. Second, a coupled bidirectional diffusion equation is used to restore the structure part, and a nonlocal means filter is used to remove noise in the oscillating part. Experimental results show advantages of this method in feature-preserving denoising.
Liu, Fu, and Zhang: Adaptive mixed image denoising based on image decomposition

1.

Introduction

Image denoising while preserving image features such as edge, detail, and texture is a key problem in image processing and computer vision. In the past decade, variational approaches and partial differential equations have been widely used in image processing:1, 2 for example, the total variation (TV)3 and the anisotropic diffusion (AD)4 for image denoising, and the shock filter for image sharpening,5 etc. However, in the process of image denoising these algorithms based on variations and partial differential equations often blur fine details and particularly textures of image.

In 2005, Buades proposed a nonlocal means filter based on the information of a whole image, which can effectively preserve image details and textures.6 Due to relatively fewer pixels with similarity around edges, however, this method also blurs part of the image edges.

Observing the poor preservation of image texture by the TV denoising, Meyer proposed a TV-G image decomposition model,7 where the G space is used to model oscillating patterns instead of the Sobolev L2 space. However, the G norm is not easy to compute numerically. To overcome this difficulty, some algorithms are proposed to approximate Mayer's model.1, 8

In this letter, to better preserve image details and textures, we propose an adaptive mixed method for image restoration. First, according to Meyer's theory, this method decomposes a given image into two components: the geometric structure and the oscillating pattern. Second, a coupled bidirectional diffusion equation9 is used to restore the structure part, and a nonlocal means filter is used to remove noise in the oscillating part. The proposed method integrates advantages of image decomposition, the excellent capability of the shock filter in edge sharpening and that of the neighborhood filter in image texture denoising.

2.

Adaptive Mixed Image Restoration

2.1.

Structure-Texture Image Decomposition

Main information of an image is encoded in its edges, details, and textures. These components need to be represented in suitable spaces to capture their special characteristics in the image.

For a given image [TeX:] \documentclass[12pt]{minimal}\begin{document}$f(x): \Omega \subset \bf {R}^2 \rightarrow \bf {R}$\end{document} f(x):ΩR2R , we decompose f into the form f = u + v, where u is the structure part including edges, details, and flat areas in BV(Ω),1 and v is the oscillating part including textures and noise in G(Ω).7 In this letter, we adopt the following decomposition model:8

1

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} {\mathop{\rm inf}\limits_{(u,v)\in BV(\Omega )\times G_\mu (\Omega )}} \left[F(u,v)=\int _\Omega\! |Du|+\frac{1}{2\lambda }||f\!-\!u\!-\!v||_{L^2}^2\right]\!,\hspace*{-1pc}\nonumber\\ \end{eqnarray} \end{document} inf (u,v)BV(Ω)×Gμ(Ω)F(u,v)=Ω|Du|+12λ||fuv||L22,
where Gμ(Ω) = {vG(Ω)| ||v||G ⩽ μ}.

In a discrete framework Aujol proved the existence and uniqueness of a solution to Eq. 1, and the approximation of Eq. 1 to Meyer's model.8

2.2.

Other Related Work

In a report5 the special idea and technique developed in numerical solution of nonlinear hyperbolic equation were applied to feature-oriented image restoration for the first time. Then, Osher and Rudin introduced a novel edge sharpening technique called shock filter.5

In order to avoid unnatural artifacts of the shock-type processing,10 we introduced a soft edge sharpening algorithm in a coupled bidirectional diffusion (CBD) equation:9

2

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{equation} \frac{{\partial u}}{{\partial t}}=c_Nu_{NN}+c_Tu_{TT}-\omega _N{\rm{sign}}(G_\sigma \ast u_{NN})|\nabla u|, \end{equation} \end{document} ut=cNuNN+cTuTTωN sign (Gσ*uNN)|u|,
where sign is a sign function, (uNN, uTT) are the second directional derivatives along local normal and tangent directions, Gσ is a Gaussian kernel function with the deviation σ, and (cN, cT, ωN) are the forward and backward diffusion coefficients, respectively.

Another effective image denoising technique is the neighborhood filter,6, 11 which considers some similarity both in spatial location and in gray level of two pixels. Buades proposed the following nonlocal means (NL-means) filter:6

3

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} &&\hspace*{-.5pc}NL_hu(x)\nonumber\\ &&\hspace*{.3pc}=\frac{1}{c(x)}\int _\Omega {\rm exp}\left[-\frac{\{G_\sigma \ast [u(x+\cdot )-u(y+\cdot )]^2\}(0)}{h^2}\right]\nonumber\\ &&\hspace*{.8pc}\times\, u(y)dy, \nonumber\\ &&\hspace*{-.5pc}\{G_\sigma \ast [u(x+\cdot )-u(y+\cdot )]^2\}(0)\nonumber\\ &&\hspace*{.3pc}=\int _\Omega G_\sigma (z)[u(x+z)-u(y+z)]^2dz, \end{eqnarray} \end{document} NLhu(x)=1c(x)Ω exp {Gσ*[u(x+·)u(y+·)]2}(0)h2×u(y)dy,{Gσ*[u(x+·)u(y+·)]2}(0)=ΩGσ(z)[u(x+z)u(y+z)]2dz,
where [TeX:] \documentclass[12pt]{minimal}\begin{document}$c(x)=\int _\Omega {\rm exp}\,[-\,\{G_\sigma \ast [u(x+\cdot )-u(y+\cdot )]^2\}\,(0)/\break h^2]\,dy$\end{document} c(x)=Ω exp [{Gσ*[u(x+·)u(y+·)]2}(0)/h2]dy is a normalization factor, and h is a filtering parameter related to noise level. The nonlocal means filter gives better results in denoising textured patterns.

2.3.

Proposed Mixed Model

In order to better preserve image details and textures, we propose an adaptive mixed method based on image decomposition (AM-ID) for image restoration, which fuses different advantages of previously mentioned methods, such as the excellent capability of the shock filter in edge sharpening and that of the neighborhood filter in image texture denoising. This method decomposes a given image f into two components: geometric structure u such as edge and detail, and oscillating pattern v including texture and noise according to Eq. 1. Then, the coupled bidirectional diffusion [Eq. 2] is used to restore the structure part, and the nonlocal means filter [Eq. 3] is used to remove noise in the oscillating part. For the sake of clarity, the AM-ID method is illustrated in Fig. 1.

Fig. 1

The proposed mixed method.

020502_1_1.jpg

Therefore, for a given degraded image f, we solve the following problem to reconstruct its restored version g:

4

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} \hspace*{-.3pc}\left\lbrace\! \begin{array}{l}{\mathop {\rm inf}\limits_{(u,v)\in BV(\Omega )\times G_\mu (\Omega )}}\!\left(F(u,v)\,{=}\int _\Omega |Du|+\frac{1}{2\lambda }||f-u-v||_{L^2}^2\right)\\[12pt] \frac{{\partial u}}{{\partial t}}=c_Nu_{NN}+c_Tu_{TT}-\omega _N\hbox{sign}(G_\sigma \ast u_{NN})|\nabla u|\\[5pt] NL_hv(x)=\frac{1}{c(x)}\int _\Omega {\rm exp}\left[-\frac{\{G_\sigma \ast [v(x+\cdot )-v(y+\cdot )]^2\}(0)}{h^2}\right]v(y)dy\\[5pt] g=\tilde{u}+\tilde{v} \end{array}\right.\! ,\hspace*{-2pc}\nonumber\\ \end{eqnarray} \end{document} inf (u,v)BV(Ω)×Gμ(Ω)F(u,v)=Ω|Du|+12λ||fuv||L22ut=cNuNN+cTuTTωNsign(Gσ*uNN)|u|NLhv(x)=1c(x)Ω exp {Gσ*[v(x+·)v(y+·)]2}(0)h2v(y)dyg=ũ+ṽ,
where [TeX:] \documentclass[12pt]{minimal}\begin{document}$\tilde{u}$\end{document} ũ and [TeX:] \documentclass[12pt]{minimal}\begin{document}$\tilde{v}$\end{document} ṽ are solutions of the CBD and the NL-means processings, respectively.

3.

Numerical Implementation and Experimental Results

3.1.

Numerical Implementation

In order to solve the model 4, we need to numerically compute Eqs. 1, 2, 3. For Eq. 1, we use the alternate projection algorithm.8 For Eq. 2, we use the numerically coupled scheme.9 Finally, for Eq. 3, a simple accelerating trick is adopted.12 On denoising oscillating patterns by the NL-means method, a smaller parameter h should be chosen to better preserve image textures.

3.2.

Experimental Results

In this section we give some experimental results to show the performance of the proposed method. Although many test images are used, we only show experiments on two typical images (Barbara and Lena) because of the imposed page limit. The test images contain rich details, textures, and large scale features. The corrupted versions are obtained by adding a Gaussian noise.

In Fig. 2, we compare the proposed method (AM-ID) with the anisotropic diffusion (AD) on a noisy Barbara (256 × 256 pixels) image. It is clear that, the given image has been successfully decomposd into the structure u and oscillating patterns v (texture and noise). At the same time, one can observe the edge sharpening of the CBD method in [TeX:] \documentclass[12pt]{minimal}\begin{document}$\tilde{u}$\end{document} ũ , and the powerful denoising of the NL-means method in [TeX:] \documentclass[12pt]{minimal}\begin{document}$\tilde{v}$\end{document} ṽ . These processings have resulted in a better image restoration, where the noise is removed effectively, edges are sharpened reasonably, and most textures and details are preserved in [TeX:] \documentclass[12pt]{minimal}\begin{document}$\tilde{u}+\tilde{v}$\end{document} ũ+ṽ . As for the AD method, although it can effectively remove the noise, it has smoothed out most of image textures. Without distinguishing image texture from noise, it elminates both noise and textures.

Fig. 2

Denoising the noisy Barbara image (from top-left to bottom-right): Original u, v + 100 by the method 1, result by the AD method, f, [TeX:] \documentclass[12pt]{minimal}\begin{document}$\tilde{u}$\end{document} ũ by the CBD processing, [TeX:] \documentclass[12pt]{minimal}\begin{document}$\tilde{v}+100$\end{document} ṽ+100 by the NL-means processing, result by the AM-ID method, respectively.

020502_1_2.jpg

In Fig. 3, we also test related methods on a noisy Lena (256 × 256 pixels) image. Once more, a better restoration of image textures and details is obtained by the proposed method. For example, see the feathers on the hat and Lena's eyes.

Fig. 3

Denoising the noisy Lena image (from left to right): original image, its noisy version, results by the AD and the AM-ID methods, respectively.

020502_1_3.jpg

Finally, we report the peak signal to noise ratio (dB) results and the computational time (seconds) in the parentheses, respectively: the noisy Barbara image (19.01), the AD method (20.42, 0.60), the AM-ID method (22.08, 6.56); the noisy Lena image (21.09), the AD method (22.59, 0.45), and the AM-ID method (23.18, 6.64). All methods are implemented using the MATLAB program on a notebook computer with 1.86 GHz CPU and 3.00 GB Memory. Because of the complexity of the proposed method, a relatively longer computational time is needed.

4.

Conclusions

Image denoising while preserving image features is an important task in image restoration, particularly for noisy textured images. This letter proposes an adaptive mixed method for image restoration, which fuses different advantages of several algorithms. Results in numerical experiments show that the proposed method produces better image restoration with effective noise removal and nice feature preserving of textures and details.

Acknowledgments

The research has been supported by the National Natural Science Foundation of China (Nos. 60933008, 61070094, 61020106001), and the China Postdoctoral Science Foundation (No. 20090460089).

References

1.  G. Aubert and P. Kornprobst, Mathematical Problems in Image Processing: Partial Differential Equations and the Calculus of Variations, 2nd. ed. Springer-Verlag, New York (2006). Google Scholar

2.  Y. Wang, R. Niu, X. Yu, L. Zhang, and H. Shen, “Image restoration and enhancement based on tunable forward-and-backward diffusion,” Opt. Eng. 49, 057004 (2010). 10.1117/1.3431657 Google Scholar

3.  L. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D 60, 259–268 (1992). 10.1016/0167-2789(92)90242-F Google Scholar

4.  P. Perona and J. Malik, “Scale-space and edge detection using anisotropic diffusion,” IEEE Trans. Pattern Anal. Mach Intell 12, 629–639 (1990). 10.1109/34.56205 Google Scholar

5.  S. Osher and L. I. Rudin, “Feature-oriented image enhancement using shock filters,” SIAM (Soc. Ind. Appl. Math.) J. Numer. Anal. 27, 919–940 (1990). Google Scholar

6.  A. Buades, B. Coll, and J. M. Morel, “A review of image denoising algorithms, with a new one,” Multiscale Model. Simul. 4, 490–530 (2005). 10.1137/040616024 Google Scholar

7.  Y. Meyer, Oscillating Patterns in Image Processing and Nonlinear Evolution Equations, American Mathematical Society, Boston, (2001). Google Scholar

8.  J. F. Aujol, G. Aubert, L. Blanc-Féraud, and A. Chambolle, “Image decomposition into a bounded variation component and an oscillating component,” J. Math. Imaging Vision 22, 71–88 (2005). 10.1007/s10851-005-4783-8 Google Scholar

9.  S. Fu, Q. Ruan, C. Mu, and W. Wang, “Geometry-driven nonlinear equation with an accelerating coupled scheme for image enhancement,” Lecture Notes Comput. Sci. 4487, 490–496 (2007). 10.1007/978-3-540-72584-8 Google Scholar

10.  S. Fu, Q. Ruan, W. Wang, F. Gao, and H. D. Cheng, “A feature-dependent fuzzy bidirectional flow for adaptive image sharpening,” Neurocomputing 70, 883–895 (2007). Google Scholar

11.  D. Q. Chen, H. Zhang, and L. Z. Cheng, “Nonlocal variational model and filter algorithm to remove multiplicative noise,” Opt. Eng. 49, 077002 (2010). 10.1117/1.3463019 Google Scholar

12.  L. Condat, “A simple trick to speed up the non-local means,” Tech. Rep., GREYC Laboratory, Centre National de la Recherche Scientifique (CNRS), France (2010). Google Scholar

© (2011) Society of Photo-Optical Instrumentation Engineers (SPIE)
Run Liu, Run Liu, Shujun Fu, Shujun Fu, Caiming Zhang, Caiming Zhang, } "Adaptive mixed image denoising based on image decomposition," Optical Engineering 50(2), 020502 (1 February 2011). https://doi.org/10.1117/1.3542041 . Submission:
JOURNAL ARTICLE
3 PAGES


SHARE
Back to Top