Anisotropic Diffusion Stencils: From Simple Derivations over Stability Estimates to ResNet Implementations

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Schrader, Karl, Weickert, Joachim, Krause, Michael
Natura: Preprint
Pubblicazione: 2023
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866914742791045120
author Schrader, Karl
Weickert, Joachim
Krause, Michael
author_facet Schrader, Karl
Weickert, Joachim
Krause, Michael
contents Anisotropic diffusion processes with a diffusion tensor are important in image analysis, physics, and engineering. However, their numerical approximation has a strong impact on dissipative artefacts and deviations from rotation invariance. In this work, we study a large family of finite difference discretisations on a 3 x 3 stencil. We derive it by splitting 2-D anisotropic diffusion into four 1-D diffusions. The resulting stencil class involves one free parameter and covers a wide range of existing discretisations. It comprises the full stencil family of Weickert et al. (2013) and shows that their two parameters contain redundancy. Furthermore, we establish a bound on the spectral norm of the matrix corresponding to the stencil. This gives time step size limits that guarantee stability of an explicit scheme in the Euclidean norm. Our directional splitting also allows a very natural translation of the explicit scheme into ResNet blocks. Employing neural network libraries enables simple and highly efficient parallel implementations on GPUs.
format Preprint
id arxiv_https___arxiv_org_abs_2309_05575
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Anisotropic Diffusion Stencils: From Simple Derivations over Stability Estimates to ResNet Implementations
Schrader, Karl
Weickert, Joachim
Krause, Michael
Numerical Analysis
Machine Learning
Image and Video Processing
Anisotropic diffusion processes with a diffusion tensor are important in image analysis, physics, and engineering. However, their numerical approximation has a strong impact on dissipative artefacts and deviations from rotation invariance. In this work, we study a large family of finite difference discretisations on a 3 x 3 stencil. We derive it by splitting 2-D anisotropic diffusion into four 1-D diffusions. The resulting stencil class involves one free parameter and covers a wide range of existing discretisations. It comprises the full stencil family of Weickert et al. (2013) and shows that their two parameters contain redundancy. Furthermore, we establish a bound on the spectral norm of the matrix corresponding to the stencil. This gives time step size limits that guarantee stability of an explicit scheme in the Euclidean norm. Our directional splitting also allows a very natural translation of the explicit scheme into ResNet blocks. Employing neural network libraries enables simple and highly efficient parallel implementations on GPUs.
title Anisotropic Diffusion Stencils: From Simple Derivations over Stability Estimates to ResNet Implementations
topic Numerical Analysis
Machine Learning
Image and Video Processing
url https://arxiv.org/abs/2309.05575