Left-invariant evolutions of wavelet transforms on the Similitude Group

Enhancement of multiple-scale elongated structures in noisy image data is relevant for many biomedical applications but commonly used PDE-based enhancement techniques often fail at crossings in an image. To get an overview of how an image is composed of local multiple-scale elongated structures we construct a multiple scale orientation score, which is a continuous wavelet transform on the similitude group, SIM(2). Our unitary transform maps the space of images onto a reproducing kernel space defined on SIM(2), allowing us to robustly relate Euclidean (and scaling) invariant operators on images to left-invariant operators on the corresponding continuous wavelet transform. Rather than often used wavelet (soft-)thresholding techniques, we employ the group structure in the wavelet domain to arrive at left-invariant evolutions and flows (diffusion), for contextual crossing preserving enhancement of multiple scale elongated structures in noisy images. We present experiments that display benefits of our work compared to recent PDE techniques acting directly on the images and to our previous work on left-invariant diffusions on orientation scores defined on Euclidean motion group.Comment: 40 page

Scale spaces on Lie groups

In the standard scale space approach one obtains a scale space representation u:R of an image by means of an evolution equation on the additive group (R d ,¿+¿). However, it is common to apply a wavelet transform (constructed via a representation of a Lie-group G and admissible wavelet ¿) to an image which provides a detailed overview of the group structure in an image. The result of such a wavelet transform provides a function on a group G (rather than (R d ,¿+¿)), which we call a score. Since the wavelet transform is unitary we have stable reconstruction by its adjoint. This allows us to link operators on images to operators on scores in a robust way. To ensure -invariance of the corresponding operator on the image the operator on the wavelet transform must be left-invariant. Therefore we focus on left-invariant evolution equations (and their resolvents) on the Lie-group G generated by a quadratic form Q on left invariant vector fields. These evolution equations correspond to stochastic processes on G and their solution is given by a group convolution with the corresponding Green’s function, for which we present an explicit derivation in two particular image analysis applications. In this article we describe a general approach how the concept of scale space can be extended by replacing the additive group R d by a Lie-group with more structure. The Dutch Organization for Scientific Research is gratefully acknowledged for financial support This article provides the theory and general framework we applied in [9],[5],[8]

Invertible Orientation Scores as an Application of Generalized Wavelet Theory,”

Abstract -Inspired by the visual system of many mammals, we consider the construction of-and reconstruction from-an orientation score of an image, via a wavelet transform corresponding to the left-regular representation of the Euclidean motion group in ‫ތ‬ 2 ( ‫ޒ‬ 2 ) and oriented wavelet ψ ∈ ‫ތ‬ 2 ( ‫ޒ‬ 2 ). Because this representation is reducible, the general wavelet reconstruction theorem does not apply. By means of reproducing kernel theory, we formulate a new and more general wavelet theory, which is applied to our specific case. As a result we can quantify the well-posedness of the reconstruction given the wavelet ψ and deal with the question of which oriented wavelet ψ is practically desirable in the sense that it both allows a stable reconstruction and a proper detection of local elongated structures. This enables image enhancement by means of left-invariant operators on orientation scores

Numerical Approaches for Linear Left-invariant Diffusions on SE(2), their Comparison to Exact Solutions, and their Applications in Retinal Imaging

Left-invariant PDE-evolutions on the roto-translation group $SE(2)$ (and their resolvent equations) have been widely studied in the fields of cortical modeling and image analysis. They include hypo-elliptic diffusion (for contour enhancement) proposed by Citti & Sarti, and Petitot, and they include the direction process (for contour completion) proposed by Mumford. This paper presents a thorough study and comparison of the many numerical approaches, which, remarkably, is missing in the literature. Existing numerical approaches can be classified into 3 categories: Finite difference methods, Fourier based methods (equivalent to $SE(2)$-Fourier methods), and stochastic methods (Monte Carlo simulations). There are also 3 types of exact solutions to the PDE-evolutions that were derived explicitly (in the spatial Fourier domain) in previous works by Duits and van Almsick in 2005. Here we provide an overview of these 3 types of exact solutions and explain how they relate to each of the 3 numerical approaches. We compute relative errors of all numerical approaches to the exact solutions, and the Fourier based methods show us the best performance with smallest relative errors. We also provide an improvement of Mathematica algorithms for evaluating Mathieu-functions, crucial in implementations of the exact solutions. Furthermore, we include an asymptotical analysis of the singularities within the kernels and we propose a probabilistic extension of underlying stochastic processes that overcomes the singular behavior in the origin of time-integrated kernels. Finally, we show retinal imaging applications of combining left-invariant PDE-evolutions with invertible orientation scores.Comment: A final and corrected version of the manuscript is Published in Numerical Mathematics: Theory, Methods and Applications (NM-TMA), vol. (9), p.1-50, 201

Nonlinear Diffusion on the 2D Euclidean Motion Group

Linear and nonlinear diffusion equations are usually considered on an image, which is in fact a function on the translation group. In this paper we study diffusion on orientation scores, i.e. on functions on the Euclidean motion group SE(2). An orientation score is obtained from an image by a linear invertible transformation. The goal is to enhance elongated structures by applying nonlinear left-invariant diffusion on the orientation score of the image. For this purpose we describe how we can use Gaussian derivatives to obtain regularized left-invariant derivatives that obey the non-commutative structure of the Lie algebra of SE(2). The Hessian constructed with these derivatives is used to estimate local curvature and orientation strength and the diffusion is made nonlinearly dependent on these measures. We propose an explicit finite difference scheme to apply the nonlinear diffusion on orientation scores. The experiments show that preservation of crossing structures is the main advantage compared to approaches such as coherence enhancing diffusion

Left-invariant Stochastic Evolution Equations on SE(2) and its Applications to Contour Enhancement and Contour Completion via Invertible Orientation Scores

We provide the explicit solutions of linear, left-invariant, (convection)-diffusion equations and the corresponding resolvent equations on the 2D-Euclidean motion group SE(2). These diffusion equations are forward Kolmogorov equations for stochastic processes for contour enhancement and completion. The solutions are group-convolutions with the corresponding Green's function, which we derive in explicit form. We mainly focus on the Kolmogorov equations for contour enhancement processes which, in contrast to the Kolmogorov equations for contour completion, do not include convection. The Green's functions of these left-invariant partial differential equations coincide with the heat-kernels on SE(2), which we explicitly derive. Then we compute completion distributions on SE(2) which are the product of a forward and a backward resolvent evolved from resp. source and sink distribution on SE(2). On the one hand, the modes of Mumford's direction process for contour completion coincide with elastica curves minimizing $\int \kappa^{2} + \epsilon ds$, related to zero-crossings of 2 left-invariant derivatives of the completion distribution. On the other hand, the completion measure for the contour enhancement concentrates on geodesics minimizing $\int \sqrt{\kappa^{2} + \epsilon} ds$. This motivates a comparison between geodesics and elastica, which are quite similar. However, we derive more practical analytic solutions for the geodesics. The theory is motivated by medical image analysis applications where enhancement of elongated structures in noisy images is required. We use left-invariant (non)-linear evolution processes for automated contour enhancement on invertible orientation scores, obtained from an image by means of a special type of unitary wavelet transform