Source author record

Gonzalo Sanguinetti

Gonzalo Sanguinetti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Neurons and Cognition Computer Vision math-ph math.AP math.GR math.MP math.RT physics.data-an

Catalog footprint

What is connected

5works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

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.

preprint2015arXiv

Sub-Riemannian Fast Marching in SE(2)

We propose a Fast Marching based implementation for computing sub-Riemanninan (SR) geodesics in the roto-translation group SE(2), with a metric depending on a cost induced by the image data. The key ingredient is a Riemannian approximation of the SR-metric. Then, a state of the art Fast Marching solver that is able to deal with extreme anisotropies is used to compute a SR-distance map as the solution of a corresponding eikonal equation. Subsequent backtracking on the distance map gives the geodesics. To validate the method, we consider the uniform cost case in which exact formulas for SR-geodesics are known and we show remarkable accuracy of the numerically computed SR-spheres. We also show a dramatic decrease in computational time with respect to a previous PDE-based iterative approach. Regarding image analysis applications, we show the potential of considering these data adaptive geodesics for a fully automated retinal vessel tree segmentation.

preprint2015arXiv

Sub-Riemannian mean curvature flow for image processing

In this paper we reconsider the sub-Riemannian cortical model of image completion introduced by Sarti and Citti in 2003. This model combines two mechanisms, the sub-Riemannian diffusion and the concentration, giving rise to a diffusion driven motion by curvature. In this paper we give a formal proof of the existence of viscosity solutions of the sub-Riemannian motion by curvature. Results of completion and enhancement on a number of natural images are shown and compared with other models.

preprint2011arXiv

Coherent states of the Euclidean group and activation regions of primary visual cortex

The uncertainty principle of SE(2) allows to construct a coherent states transform that is strictly related to the Bargmann transform for the second Heisenberg group H2. The corresponding target space is characterized constructively and related to the almost complex structure of SE(2) as a contact manifold. Such a coherent state transform provides a model for neural activity maps in the primary visual cortex, that are then described in terms of minimal uncertainty states. The results of the model are compared with the experimental measurements.

preprint2010arXiv

An uncertainty principle underlying the pinwheel structure in the primary visual cortex

The visual information in V1 is processed by an array of modules called orientation preference columns. In some species including humans, orientation columns are radially arranged around singular points like the spokes of a wheel, that are called pinwheels. The pinwheel structure has been observed first with optical imaging techniques and more recently by in vivo two-photon imaging proving their organization with single cell precision. In this research we provide evidence that pinwheels are de facto optimal distributions for coding at the best angular position and momentum. In the last years many authors have recognized that the functional architecture of V1 is locally invariant with respect to the symmetry group of rotations and translations SE(2). In the present study we show that the orientation cortical maps used to construct pinwheels can be modeled as coherent states, i.e. the configurations best localized both in angular position and angular momentum. The theory we adopt is based on the well known uncertainty principle, introduced by Heisenberg in quantum mechanics and later extended to many other groups of invariance. Here we state a corresponding principle in the cortical geometry with SE(2) symmetry, and by computing its minimizers we obtain a model of orientation activity maps in the cortex. As it is well known the pinwheels configuration is directly constructed from these activity maps, and we will be able to formally reproduce their structure starting from the group symmetries of the functional architecture of the visual cortex. The primary visual cortex is then modeled as an integrated system in which the set of simple cells implements the SE(2) group, the horizontal connectivity implements its Lie algebra and the pinwheels implement its minimal uncertainty states.