Well-posed anisotropic diffusion for image denoising

A nonlinear iterative smoothing filter based on a second-order partial differential equation is introduced. It smooths out the image according to an anisotropic diffusion process. The approach is based on a smooth approximation of the total variation (TV) functional which overcomes the non-differentiability of the TV functional at the origin. In particular, the authors perform linear smoothing over smooth areas but selective smoothing over candidate edges. By relating the smoothing parameter to the time step, they arrive at a CFL condition which guarantees the causality of the discrete scheme. This allows the adoption of higher time discretisation steps, while ensuring the absence of artefacts deriving from the non-smooth behaviour of the TV functional at the origin. In particular, it is shown that the proposed approach avoids the typical staircase effects in smooth areas which occur in the standard time-marching TV scheme.

[1]  E. Giusti Minimal surfaces and functions of bounded variation , 1977 .

[2]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[3]  E. Giusti Functions of Bounded Variation and Caccioppoli Sets , 1984 .

[4]  Robert Hummel,et al.  Reconstructions from zero crossings in scale space , 1989, IEEE Trans. Acoust. Speech Signal Process..

[5]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[6]  R. LeVeque Numerical methods for conservation laws , 1990 .

[7]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[8]  P. Lions,et al.  Axioms and fundamental equations of image processing , 1993 .

[9]  Wenyuan Xu,et al.  Analysis and design of anisotropic diffusion for image processing , 1994, Proceedings of 1st International Conference on Image Processing.

[10]  C. Vogel,et al.  Analysis of bounded variation penalty methods for ill-posed problems , 1994 .

[11]  Raymond H. Chan,et al.  Continuation method for total variation denoising problems , 1995, Optics & Photonics.

[12]  Stan Z. Li,et al.  On Discontinuity-Adaptive Smoothness Priors in Computer Vision , 1995, IEEE Trans. Pattern Anal. Mach. Intell..

[13]  Yuying Li,et al.  A computational algorithm for minimizing total variation in image restoration , 1996, IEEE Trans. Image Process..

[14]  Wenyuan Xu,et al.  Behavioral analysis of anisotropic diffusion in image processing , 1996, IEEE Trans. Image Process..

[15]  K. E. Tait,et al.  Image recovery using the anisotropic diffusion equation , 1996, IEEE Trans. Image Process..

[16]  Satyanad Kichenassamy,et al.  The Perona-Malik Paradox , 1997, SIAM J. Appl. Math..

[17]  P. Lions,et al.  Image recovery via total variation minimization and related problems , 1997 .

[18]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[19]  J. Koenderink The structure of images , 2004, Biological Cybernetics.