Nonlinear diffusion filtering is usually performed with explicit schemes. They are only stable for very small time steps, which leads to poor efficiency and limits their practical use. Based on a recent discrete nonlinear diffusion scale-space framework we present semi-implicit schemes which are stable for all time steps. These novel schemes use an additive operator splitting (AOS) which guarantees equal treatment of all coordinate axes. They can be implemented easily in arbitrary dimensions, have good rotational invariance and reveal a computational complexity and memory requirement which is linear in the number of pixels. Examples demonstrate that, under typical accuracy requirements, AOS schemes are at least ten times more efficient than the widely-used explicit schemes.

Abstract—Diffusions are useful for image processing and computer vision because they provide a convenient way of smoothing noisy data, analyzing images at multiple scales, and enhancing discontinuities. A number of diffusions of image brightness have been defined and studied so far; they may be applied to scalar and vector-valued quantities that are naturally associated with intervals of either the real line, or other flat manifolds. Some quantities of interest in computer vision, and other areas of engineering that deal with images, are defined on curved manifolds; typical examples are orientation and hue that are defined on the circle. Generalizing brightness diffusions to orientation is not straightforward, especially in the case where a discrete implementation is sought. An example of what may go wrong is presented. A method is proposed to define diffusions of orientation-like quantities. First a definition in the continuum is discussed, then a discrete orientation diffusion is proposed. The behavior of such diffusions is explored both analytically and experimentally. It is shown how such orientation diffusions contain a nonlinearity that is reminiscent of edge-process and anisotropic diffusion. A number of open questions are proposed at the end. Index Terms—Orientation analysis, texture analysis, diffusions, scale-space.

We introduce a fast, multiscale algorithm for image segmentation. Our algorithm uses modern numeric techniques to nd an approximate solution to normalized cut measures in time that is linear in the size of the image with only a few dozen operations per pixel. In just one pass the algorithm provides a complete hierarchical decomposition of the image into segments. The algorithm detects the segments by applying a process of recursive coarsening in which the same minimization problem is represented with fewer and fewer variables producing an irregular pyramid. During this coarsening process we may compute additional internal statistics of the emerging segments and use these statistics to facilitate the segmentation process. Once the pyramid is completed it is scanned from the top down to associate pixels close to the boundaries of segments with the appropriate segment. The algorithm is inspired by algebraic multigrid (AMG) solvers of minimization problems of heat or electric networks. We demonstrate the algorithm by applying it to real images.

. This paper gives an overview of scale-space and image enhancement techniques which are based on parabolic partial differential equations in divergence form. In the nonlinear setting this filter class allows to integrate a-priori knowledge into the evolution. We sketch basic ideas behind the different filter models, discuss their theoretical foundations and scale-space properties, discrete aspects, suitable algorithms, generalizations, and applications. 1 Introduction During the last decade nonlinear diffusion filters have become a powerful and well-founded tool in multiscale image analysis. These models allow to include a-priori knowledge into the scale-space evolution, and they lead to an image simplification which simultaneously preserves or even enhances semantically important information such as edges, lines, or flow-like structures. Many papers have appeared proposing different models, investigating their theoretical foundations, and describing interesting applications. For a n...

The completion of interrupted lines or the enhancement of flow-like structures is a challenging task in computer vision, human vision, and image processing. We address this problem by presenting a multiscale method in which a nonlinear diffusion filter is steered by the so-called interest operator (second-moment matrix, structure tensor). An m-dimensional formulation of this method is analysed with respect to its well-posedness and scale-space properties. An efficient scheme is presented which uses a stabilization by a semi-implicit additive operator splitting (AOS), and the scale-space behaviour of this method is illustrated by applying it to both 2-D and 3-D images.

This work presents two mechanisms for processing fingerprint images; shape-adapted smoothing based on second moment descriptors and automatic scale selection based on normalized derivatives. The shape adaptation procedure adapts the smoothing operation to the local ridge structures, which allows interrupted ridges to be joined without destroying essential singularities such as branching points and enforces continuity of their directional fields. The scale selection procedure estimates local ridge width and adapts the amount of smoothing to the local amount of noise. In addition, a ridgeness measure is defined, which reflects how well the local image structure agrees with a qualitative ridge model, and is used for spreading the results of shape adaptation into noisy areas. The combined approach makes it possible to resolve fine scale structures in clear areas while reducing the risk of enhancing noise in blurred or fragmented areas. The result is a reliable and adaptively detailed estimate of the ridge orientation field and ridge width, as well as a smoothed grey-level version of the input image. We propose that these general techniques should be of interest to developers of automatic fingerprint identification

We introduce a family of first-order multidimensional ordinary differential equations (ODEs) with discontinuous right-hand sides and demonstrate their applicability in image processing. An equation belonging to this family is an inverse diffusion everywhere except at local extrema, where some stabilization is introduced. For this reason, we call these equations "stabilized inverse diffusion equations" ("SIDEs"). Existence and uniqueness of solutions, as well as stability, are proven for SIDEs. A SIDE in one spatial dimension may be interpreted as a limiting case of a semi-discretized Perona-Malik equation [14], [15]. In an experimental section, SIDEs are shown to suppress noise while sharpening edges present in the input signal. Their application to image segmentation is also demonstrated.

Abstract—In this paper, we use partial-differential-equationbased filtering as a preprocessing and post processing strategy for computer-aided cytology. We wish to accurately extract and classify the shapes of nuclei from confocal microscopy images, which is a prerequisite to an accurate quantitative intranuclear (genotypic and phenotypic) and internuclear (tissue structure) analysis of tissue and cultured specimens. First, we study the use of a geometry-driven edge-preserving image smoothing mechanism before nuclear segmentation. We show how this filter outperforms other widely-used filters in that it provides higher edge fidelity. Then we apply the same filter, with a different initial condition, to smooth nuclear surfaces and obtain sub-pixel accuracy. Finally we use another instance of the geometrical filter to correct for misinterpretations of the nuclear surface by the segmentation algorithm. Our prefiltering and post filtering nicely complements our initial segmentation strategy, in that it provides substantial and measurable improvement in the definition of the nuclear surfaces. Index Terms—Cytology, differential geometry, dynamic surfaces, image processing, level sets, Riemannian geometry, segmentation, surface evolution. I.

. Poor efficiency is a typical problem of nonlinear diffusion filtering, when the simple and popular explicit (Euler-forward) scheme is used: for stability reasons very small time step sizes are necessary. In order to overcome this shortcoming, a novel type of semi-implicit schemes is studied, so-called additive operator splitting (AOS) methods. They share the advantages of explicit and (semi-)implicit schemes by combining simplicity with absolute stability. They are reliable, since they satisfy recently established criteria for discrete nonlinear diffusion scale-spaces. Their efficiency is due to the fact that they can be separated into one-dimensional processes, for which a fast recursive algorithm with linear complexity is available. AOS schemes reveal good rotational invariance and they are symmetric with respect to all axes. Examples demonstrate that, under typical accuracy requirements, they are at least ten times more efficient than explicit schemes. 1 Introduction Although non...

The goal of this paper is to present segmentation algorithms which combine regularization by nonlinear partial differential equations (PDEs) with a watershed transformation with region merging. We develop efficient algorithms for two wellfounded PDE methods. They use an additive operator splitting (AOS) leading to recursive and separable filters. Further speed-up can be obtained by embedding AOS schemes into a pyramid framework. Examples are presented which demonstrate that the preprocessing by these PDE techniques eases and stabilizes the segmentation. The typical CPU time for segmenting a 256 2 image on a workstation is less than 2 seconds. Key Words: Nonlinear diffusion, additive operator splitting, Gaussian pyramid, watershed segmentation, region merging CR Subject Classification: I.4.6, I.4.3, I.4.4. 1 Introduction Segmentation is one of the bottlenecks of many image analysis and computer vision tasks ranging from medical image processing to robot navigation. Ideally it sho...