Mimetic finite difference methods in image processing

We introduce the use of mimetic methods to the imaging community, for the solution of the initial-value problems ubiquitous in the machine vision and image processing and analysis fields. PDE-based image processing and analysis techniques comprise a host of applications such as noise removal and restoration, deblurring and enhancement, segmentation, edge detection, inpainting, registration, motion analysis, etc. Because of their favorable stability and efficiency properties, semi-implicit finite difference and finite element schemes have been the methods of choice (in that order of preference). We propose a new approach for the numerical solution of these problems based on mimetic methods. The mimetic discretization scheme preserves the continuum properties of the mathematical operators often encountered in the image processing and analysis equations. This is the main contributing factor to the improved performance of the mimetic method approach, as compared to both of the aforementioned popular numerical solution techniques. To assess the performance of the proposed approach, we employ the Catte-Lions-Morel-Coll model to restore noisy images, by solving the PDE with the three numerical solution schemes. For all of the benchmark images employed in our experiments, and for every level of noise applied, we observe that the best image restored by using the mimetic method is closer to the noise-free image than the best images restored by the other two methods tested. These results motivate further studies of the application of the mimetic methods to other imaging problems. Mathematical subject classification: Primary: 68U10; Secondary: 65L12.

[1]  Fiorella Sgallari,et al.  Variational Numerical Methods for Solving Nonlinear Diffusion Equations Arising in Image Processing , 2002, J. Vis. Commun. Image Represent..

[2]  Mikhail Shashkov,et al.  The sensitivity and accuracy of fourth order finite-difference schemes on nonuniform grids in one dimension , 1995 .

[3]  K. Höllig,et al.  A Diffusion Equation with a Nonmonotone Constitutive Function , 1983 .

[4]  Joachim Weickert,et al.  PDE - Based Preprocessing of Medical Images , 2000, Künstliche Intell..

[5]  José E. Castillo,et al.  Linear Systems Arising for Second-Order Mimetic Divergence and Gradient Discretizations , 2005, J. Math. Model. Algorithms.

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

[7]  J. E. Castillo,et al.  A Matrix Analysis Approach to Higher-Order Approximations for Divergence and Gradients Satisfying a Global Conservation Law , 2003, SIAM J. Matrix Anal. Appl..

[8]  Joachim Weickert,et al.  Anisotropic diffusion in image processing , 1996 .

[9]  Mikhail Shashkov,et al.  Fourth- and sixth-order conservative finite difference approximations of the divergence and gradient , 2001 .

[10]  Fiorella Sgallari,et al.  Semi-implicit complementary volume scheme for solving level set like equations in image processing and curve evolution , 2003, Numerische Mathematik.

[11]  K. Mikula,et al.  A coarsening finite element strategy in image selective smoothing , 1997 .

[12]  Pavel B. Bochev,et al.  Principles of Mimetic Discretizations of Differential Operators , 2006 .

[13]  K. Höllig.  Existence of infinitely many solutions for a forward backward heat equation , 1983 .

[14]  Mark Nitzberg,et al.  Nonlinear Image Filtering with Edge and Corner Enhancement , 1992, IEEE Trans. Pattern Anal. Mach. Intell..

[15]  Xiaoping Li,et al.  Nonlinear diffusion with multiple edginess thresholds , 1994, Pattern Recognit..

[16]  James M. Hyman,et al.  The convergence of mimetic discretization for rough grids , 2004 .

[17]  Martin Rumpf,et al.  DIFFUSION MODELS AND THEIR ACCELERATED SOLUTION IN IMAGE AN D SURFACE PROCESSING , 2000 .

[18]  John F. Canny,et al.  A Computational Approach to Edge Detection , 1986, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[19]  M. Aurada,et al.  Convergence of adaptive BEM for some mixed boundary value problem , 2012, Applied numerical mathematics : transactions of IMACS.

[20]  Alessandro Sarti,et al.  Numerical solution of parabolic equations related to level set formulation of mean curvature flow , 1998 .

[21]  J. Weickert,et al.  Image Processing Using a Wavelet Algorithm for Nonlinear Diffusion , 1994 .

[22]  P. Lions,et al.  Image selective smoothing and edge detection by nonlinear diffusion. II , 1992 .

[23]  Eberhard Bänsch,et al.  Adaptivity in 3D image processing , 2001 .

[24]  Jitendra Malik,et al.  Anisotropic Diffusion , 1994, Geometry-Driven Diffusion in Computer Vision.

[25]  Karol Mikula,et al.  Semi-implicit finite volume scheme for solving nonlinear diffusion equations in image processing , 2001, Numerische Mathematik.

[26]  B. Aulbach,et al.  Linear Dynamic Processes with Inhomogeneous Time Scale , 1990, Nonlinear Dynamics and Quantum Dynamical Systems.

[27]  D Marr,et al.  Theory of edge detection , 1979, Proceedings of the Royal Society of London. Series B. Biological Sciences.

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

[29]  Max A. Viergever,et al.  Efficient and reliable schemes for nonlinear diffusion filtering , 1998, IEEE Trans. Image Process..

[30]  Martin Rumpf,et al.  An Adaptive Finite Element Method for Large Scale Image Processing , 1999, J. Vis. Commun. Image Represent..

[31]  Joachim Weickert,et al.  Efficient image segmentation using partial differential equations and morphology , 2001, Pattern Recognit..

[32]  M. Steinberg.  High-Order Mimetic Finite Difference Methods on Nonuniform Grids , 1995 .

[33]  Martin Bohner,et al.  Mimetic methods on measure chains , 2001 .

[34]  Joachim Weickert,et al.  Coherence-Enhancing Diffusion Filtering , 1999, International Journal of Computer Vision.

[35]  J. Kacur,et al.  Solution of nonlinear diffusion appearing in image smoothing and edge detection , 1995 .

[36]  G. I. Barenblatt,et al.  A degenerate pseudoparabolic regularization of a nonlinear forward-backward heat equation arising in the theory of heat and mass exchange in stably stratified turbulent shear flow , 1993 .

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