Behavioral analysis of anisotropic diffusion in image processing

©1996 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or distribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from the IEEE. This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.DOI: 10.1109/83.541424In this paper, we analyze the behavior of the anisotropic diffusion model of Perona and Malik (1990). The main idea is to express the anisotropic diffusion equation as coming from a certain optimization problem, so its behavior can be analyzed based on the shape of the corresponding energy surface. We show that anisotropic diffusion is the steepest descent method for solving an energy minimization problem. It is demonstrated that an anisotropic diffusion is well posed when there exists a unique global minimum for the energy functional and that the ill posedness of a certain anisotropic diffusion is caused by the fact that its energy functional has an infinite number of global minima that are dense in the image space. We give a sufficient condition for an anisotropic diffusion to be well posed and a sufficient and necessary condition for it to be ill posed due to the dense global minima. The mechanism of smoothing and edge enhancement of anisotropic diffusion is illustrated through a particular orthogonal decomposition of the diffusion operator into two parts: one that diffuses tangentially to the edges and therefore acts as an anisotropic smoothing operator, and the other that flows normally to the edges and thus acts as an enhancement operator

Modeling anisotropic diffusion using a departure from isotropy approach

There are a large number of finite volume solvers available for solution of isotropic diffusion equation. This article presents an approach of adapting these solvers to solve anisotropic diffusion equations. The formulation works by decomposing the diffusive flux into a component associated with isotropic diffusion and another component associated with departure from isotropic diffusion. This results in an isotropic diffusion equation with additional terms to account for the anisotropic effect. These additional terms are treated using a deferred correction approach and coupled via an iterative procedure. The presented approach is validated against various diffusion problems in anisotropic media with known analytical or numerical solutions. Although demonstrated for two-dimensional problems, extension of the present approach to three-dimensional problems is straight forward. Other than the finite volume method, this approach can be applied to any discretization method

Anisotropic enhanced backscattering induced by anisotropic diffusion

The enhanced backscattering cone displaying a strong anisotropy from a material with anisotropic diffusion is reported. The constructive interference of the wave is preserved in the helicity preserving polarization channel and completely lost in the nonpreserving one. The internal reflectivity at the interface modifies the width of the backscatter cone. The reflectivity coefficient is measured by angular-resolved transmission. This interface property is found to be isotropic, simplifying the backscatter cone analysis. The material used is a macroporous semiconductor, gallium phosphide, in which pores are etched in a disordered position but with a preferential direction

A network for multiscale image segmentation

Detecting edges of objects in their images is a basic problem in computational vision. The scale-space technique introduced by Witkin [11] provides means of using local and global reasoning in locating edges. This approach has a major drawback: it is difficult to obtain accurately the locations of the 'semantically meaningful' edges. We have refined the definition of scale-space, and introduced a class of algorithms for implementing it based on using anisotropic diffusion [9]. The algorithms involves simple, local operations replicated over the image making parallel hardware implementation feasible. In this paper we present the major ideas behind the use of scale space, and anisotropic diffusion for edge detection, we show that anisotropic diffusion can enhance edges, we suggest a network implementation of anisotropic diffusion, and provide design criteria for obtaining networks performing scale space, and edge detection. The results of a software implementation are shown

Enforcing the non-negativity constraint and maximum principles for diffusion with decay on general computational grids

In this paper, we consider anisotropic diffusion with decay, and the diffusivity coefficient to be a second-order symmetric and positive definite tensor. It is well-known that this particular equation is a second-order elliptic equation, and satisfies a maximum principle under certain regularity assumptions. However, the finite element implementation of the classical Galerkin formulation for both anisotropic and isotropic diffusion with decay does not respect the maximum principle. We first show that the numerical accuracy of the classical Galerkin formulation deteriorates dramatically with increase in the decay coefficient for isotropic medium and violates the discrete maximum principle. However, in the case of isotropic medium, the extent of violation decreases with mesh refinement. We then show that, in the case of anisotropic medium, the classical Galerkin formulation for anisotropic diffusion with decay violates the discrete maximum principle even at lower values of decay coefficient and does not vanish with mesh refinement. We then present a methodology for enforcing maximum principles under the classical Galerkin formulation for anisotropic diffusion with decay on general computational grids using optimization techniques. Representative numerical results (which take into account anisotropy and heterogeneity) are presented to illustrate the performance of the proposed formulation

Anisotropic Diffusion Limited Aggregation

Using stochastic conformal mappings we study the effects of anisotropic perturbations on diffusion limited aggregation (DLA) in two dimensions. The harmonic measure of the growth probability for DLA can be conformally mapped onto a constant measure on a unit circle. Here we map $m$ preferred directions for growth of angular width $\sigma$ to a distribution on the unit circle which is a periodic function with $m$ peaks in $[-\pi, \pi)$ such that the width $\sigma$ of each peak scales as $\sigma \sim 1/\sqrt{k}$, where $k$ defines the ``strength'' of anisotropy along any of the $m$ chosen directions. The two parameters $(m,k)$ map out a parameter space of perturbations that allows a continuous transition from DLA (for $m=0$ or $k=0$) to $m$ needle-like fingers as $k \to \infty$. We show that at fixed $m$ the effective fractal dimension of the clusters $D(m,k)$ obtained from mass-radius scaling decreases with increasing $k$ from $D_{DLA} \simeq 1.71$ to a value bounded from below by $D_{min} = 3/2$. Scaling arguments suggest a specific form for the dependence of the fractal dimension $D(m,k)$ on $k$ for large $k$, form which compares favorably with numerical results.Comment: 6 pages, 4 figures, submitted to Phys. Rev.

On the striated regularity for the 2D anisotropic Boussinesq system

In this paper, we investigate the global existence and uniqueness of strong solutions to 2D Boussinesq system with anisotropic thermal diffusion or anisotropic viscosity and with striated initial data. Using the key idea of Chemin to solve 2-D vortex patch of ideal fluid, namely the striated regularity can help to bound the gradient of the velocity, we can prove the global well-posedness of the Boussinesq system with anisotropic thermal diffusion with initial vorticity being discontinuous across some smooth interface. In the case of an anisotropic horizontal viscosity we can study the propagation of the striated regularity for the smooth temperature patches problem.Comment: 36 page