Novel Schemes for Hyperbolic PDEs Using Osmosis Filters from Visual Computing

Recently a new class of generalised diffusion filters called osmosis filters has been proposed. Osmosis models are useful for a variety of tasks in visual computing. In this paper, we show that these filters are also beneficial outside image processing and computer graphics: We exploit their use for the construction of better numerical schemes for hyperbolic partial differential equations that model physical transport phenomena. Our novel osmosis-based algorithm is constructed as a two-step, predictor-corrector method. The predictor scheme is given by a Markov chain model of osmosis that captures the hyperbolic transport in its advection term. By design, it also incorporates a discrete diffusion process. The corresponding terms can easily be identified within the osmosis model. In the corrector step, we subtract a stabilised version of this discrete diffusion. We show that the resulting osmosis-based method gives correct, highly accurate resolutions of shock wave fronts in both linear and nonlinear test cases. Our work is an example for the usefulness of visual computing ideas in numerical analysis.

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

[2]  Michael Breuß,et al.  Highly Accurate Schemes for PDE-Based Morphology with General Convex Structuring Elements , 2011, International Journal of Computer Vision.

[3]  Michael Breuß,et al.  Adaptive Continuous-Scale Morphology for Matrix Fields , 2011, International Journal of Computer Vision.

[4]  Michael Breuß,et al.  A Shock-Capturing Algorithm for the Differential Equations of Dilation and Erosion , 2006, Journal of Mathematical Imaging and Vision.

[5]  Michael Breuß,et al.  Osmosis models for visual computing , 2011 .

[6]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Multidimensional Systems , 2002 .

[7]  Yehoshua Y. Zeevi,et al.  Forward-and-backward diffusion processes for adaptive image enhancement and denoising , 2002, IEEE Trans. Image Process..

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

[9]  E. Seneta.  Non-negative Matrices and Markov Chains , 2008 .

[10]  Michael Breuß,et al.  Can Variational Models for Correspondence Problems Benefit from Upwind Discretisations? , 2011, Journal of Mathematical Imaging and Vision.

[11]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[12]  Knut-Andreas Lie,et al.  Scale Space and Variational Methods in Computer Vision, Second International Conference, SSVM 2009, Voss, Norway, June 1-5, 2009. Proceedings , 2009, SSVM.

[13]  Thorsten Grahs,et al.  Image Processing for Numerical Approximations of Conservation Laws: Nonlinear Anisotropic Artificial Dissipation , 2001, SIAM J. Sci. Comput..

[14]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[15]  L. Rudin,et al.  Feature-oriented image enhancement using shock filters , 1990 .

[16]  L. Rudin.  Images, Numerical Analysis of Singularities and Shock Filters , 1987 .

[17]  Thomas Brox,et al.  Stabilised Nonlinear Inverse Diffusion for Approximating Hyperbolic PDEs , 2005, Scale-Space.

[18]  Michael Breuss,et al.  An analysis of the influence of data extrema on some first and second order central approximations of hyperbolic conservation laws , 2005 .

[19]  Joachim Weickert,et al.  Theoretical Foundations for Discrete Forward-and-Backward Diffusion Filtering , 2009, SSVM.

[20]  Alan S. Willsky,et al.  Image segmentation and edge enhancement with stabilized inverse diffusion equations , 2000, IEEE Trans. Image Process..

[21]  M. Breuß.  The correct use of the Lax–Friedrichs method , 2004 .

[22]  Stanley Osher,et al.  Shocks and other nonlinear filtering applied to image processing , 1991, Optics & Photonics.

[23]  S. Zalesak.  Introduction to “Flux-Corrected Transport. I. SHASTA, A Fluid Transport Algorithm That Works” , 1997 .

[24]  L. Álvarez,et al.  Signal and image restoration using shock filters and anisotropic diffusion , 1994 .

[25]  Joachim Weickert,et al.  Scale Space and PDE Methods in Computer Vision, 5th International Conference, Scale-Space 2005, Hofgeismar, Germany, April 7-9, 2005, Proceedings , 2005, Scale-Space.

[26]  Rachid Deriche,et al.  Image restoration via PDE , 1997, Defense + Security.

[27]  F C Billingsley,et al.  Applications of digital image processing. , 1970, Applied optics.

[28]  E. Seneta.  Non-negative Matrices and Markov Chains (Springer Series in Statistics) , 1981 .

[29]  J. Boris,et al.  Flux-Corrected Transport , 1997 .

[30]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[31]  M. Welk,et al.  Staircasing in semidiscrete stabilised inverse linear diffusion algorithms , 2007 .

[32]  G. Wei,et al.  Shock capturing by anisotropic diffusion oscillation reduction , 2000 .

[33]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[34]  Thorsten Grahs,et al.  Entropy-Controlled Artificial Anisotropic Diffusion for the Numerical Solution of Conservation Laws Based on Algorithms from Image Processing , 2002, J. Vis. Commun. Image Represent..