Next Article in Journal
Word Spotting as a Service: An Unsupervised and Segmentation-Free Framework for Handwritten Documents
Next Article in Special Issue
Retinal Processing: Insights from Mathematical Modelling
Previous Article in Journal
A Pipelined Tracer-Aware Approach for Lesion Segmentation in Breast DCE-MRI
Previous Article in Special Issue
Multi-Frequency Image Completion via a Biologically-Inspired Sub-Riemannian Model with Frequency and Phase
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Liouville Integrability in a Four-Dimensional Model of the Visual Cortex

by
Ivan Galyaev
1,† and
Alexey Mashtakov
2,*,†,‡
1
V. A. Trapeznikov Institute of Control Sciences of RAS, 117997 Moscow, Russia
2
Ailamazyan Program Systems Institute of RAS, 152021 Veskovo, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Current address: Yaroslavl Region, Pereslavsky District, St. Petera I, 4a, 152021 Veskovo, Russia.
J. Imaging 2021, 7(12), 277; https://doi.org/10.3390/jimaging7120277
Submission received: 4 October 2021 / Revised: 28 November 2021 / Accepted: 13 December 2021 / Published: 17 December 2021

Abstract

:
We consider a natural extension of the Petitot–Citti–Sarti model of the primary visual cortex. In the extended model, the curvature of contours is taken into account. The occluded contours are completed via sub-Riemannian geodesics in the four-dimensional space M of positions, orientations, and curvatures. Here, M = R 2 × SO ( 2 ) × R models the configuration space of neurons of the visual cortex. We study the problem of sub-Riemannian geodesics on M via methods of geometric control theory. We prove complete controllability of the system and the existence of optimal controls. By application of the Pontryagin maximum principle, we derive a Hamiltonian system that describes the geodesics. We obtain the explicit parametrization of abnormal extremals. In the normal case, we provide three functionally independent first integrals. Numerical simulations indicate the existence of one more first integral that results in Liouville integrability of the system.

1. Introduction

The principles of biological visual systems are of great interest among researchers in many fields of science. An important problem is the development and study of realistic mathematical models describing a certain stage of visual signal processing. In this paper, we perform a study of a mathematical model that describes a mechanism of visual signal processing by the primary visual cortex of the human brain.
Human vision is a complex process that has not yet been fully understood. Perception of visual information starts in the eye. A ray of light hits the retina and causes the activation of light-sensitive receptors. The retina is a part of the brain hierarchically organized in several layers. The first layer consists of light-sensitive receptors. In the next layers of the retina, bipolar and ganglion cells are contained. Here, the primary processing of the visual signal takes place. Furthermore, the visual signal arrives through the optic nerve to the lateral geniculate nucleus (LGN) of the thalamus. Afterward, the signal is transmitted to the visual cortex of the brain.
The mechanism of visual signal processing by neurons of the retina and LGN have been extensively studied. The scale-space theory [1] proposed a mathematical model for this stage. This theory was inspired by the properties of the Gaussian kernel and its derivatives as regularized differential operators as well as solutions to the linear diffusion equation. The receptive fields of bipolar and ganglion cells as well as LGN neurons are well approximated by filter profiles based on the Gaussian kernel and Gaussian derivatives. Their work is modeled by the action of the filter to the input signal. Mathematically, this operation is defined as the convolution of two functions.
After LGN, the signal is transmitted to the visual cortex of the brain. Physiological research indicates that the visual cortex is composed of multiple layers. The research by Hubel and Wiesel [2] has made significant progress in understanding the principles of the primary visual cortex V1. In particular, they found that the receptive fields of V1 neurons are elongated rather than rounded. The cells of V1 are capable of detecting segments of contours with different orientations from the whole image. Mathematically, the work of V1 cells can be understood as a lift of a 2D input image to the extended space of positions and orientations SE ( 2 ) = R 2 × SO ( 2 ) .
The mathematical model of V1 as a sub-Riemannian structure on the Heisenberg group was proposed by J. Petitot [3]. Then, this model was refined by G. Citti and A. Sarti [4] as a sub-Riemannian structure on the Lie group SE ( 2 ) = R 2 × SO ( 2 ) . Here, SE ( 2 ) models the configuration space of neurons V1, which can be understood as the space of positions R 2 and orientations SO ( 2 ) . According to this model, the process of completion of occluded contours occurs by minimizing the excitation energy of neurons that perceive visual information in the areas of the observed scene that are occluded. This process can be interpreted as the action of the hypoelliptic diffusion operator studied in [5,6,7]. The resulting curves are sub-Riemannian length minimizers in SE ( 2 ) . Their exact parameterization was obtained in [8]. Such curves are used for image inpainting [9] and for the explanation of some visual illusions [10].
The principles of biological visual systems are actively used in computer vision. Based on these principles, effective methods of image processing are created: enhancement, segmentation, inpainting, and feature detection. For example, in [11], the authors describe an approach that is based on lifting an image into the extended space of positions and orientations. After such lifting, the sub-Riemannian length minimizers are used to detect the salient lines [12,13,14]. This approach is actively used in medical image processing.
The experimental data [15] (see also [16]) suggest that not only detectors of orientation but also detectors of curvature exist in V1. In this paper, we consider the mathematical model of the visual cortex ([17], p. 57) [18] that is obtained by extension of the classical Petitot–Citti–Sarti model. The extension is performed by taking into account the curvature of contours of the observed image. This leads to a sub-Riemannian structure in the four-manifold M = R 2 × SO ( 2 ) × R , where the fourth component means the curvature. In the four-dimensional model of the visual cortex, a partially occluded contour is completed via the planar projection of a sub-Riemannian length-minimizer in M that satisfies the boundary conditions (position, orientation, and curvature) obtained from the boundary of the occluded region. See Figure 1.
In [18], the authors show that the 4D model allows for better perceptual grouping and completion of complex images than the classical 3D model. They consider the completion of images via minimal surfaces in the 4D model. In this paper, we study the related problem of finding length-minimizers of the associated sub-Riemannian manifold. In analogy with [9], knowledge of length-minimizers provides a method for the completion of isophotes of corrupted images. Such a method is an alternative to the minimal surface method, proposed in [18], where the linearized Euler–Lagrange equation is used to compute the surface. Such a linearization provides an approximation for the geodesic flow. We expect that the usage of the precise length-minimizers will provide a more accurate method for image completion since it does not involve linearization.
In this article, we consider the problem of sub-Riemannian geodesics in the space M of positions, orientations, and curvatures. We formulate the problem as an optimal control problem in Section 2. In Section 3, we prove complete controllability of the system and the existence of optimal controls. Then, in Section 4, we apply a necessary optimality condition—the Pontryagin maximum principle (PMP)—and examine the Hamiltonian system of PMP. We obtain explicit parametrization of abnormal extremals and provide a numerical investigation resulting in the Liouville integrability of the normal Hamiltonian system of PMP. In conclusion, we summarize the main results.
Note that the question of integrability of the geodesic flow in M is important both for theory and for applications. Liouville integrability ensures the absence of undesirable chaotic behavior. It guarantees that the trajectories of the system remain close to each other under a small perturbation of the initial value. This property is highly important for the stability of numerical schemes of integration. Note that integrability is a rare phenomenon. In the common case, a system of ODE is not integrable. In [19], the authors show that, even in the simplest case of Carnot groups, the sub-Riemannian geodesic flow is not integrable already in dimension 8 (and consequently in higher dimensions).

2. Optimal Control Problem

We consider the following control system:
x ˙ = u 1 cos θ , y ˙ = u 1 sin θ , θ ˙ = u 1 k , k ˙ = u 2 , ( x , y , θ , k ) = q SE ( 2 ) × R = M , ( u 1 , u 2 ) R 2 .
For the trajectory corresponding to the control ( u 1 ( t ) , u 2 ( t ) ) , t [ 0 , T ] , T > 0 , we define a cost functional—the sub-Riemannian length of this trajectory:
l = 0 T u 1 2 ( t ) + u 2 2 ( t ) d t .
We study the problem of finding a Lipschitzian curve γ : [ 0 , T ] M —that is a trajectory of system (1), satisfying the given boundary conditions
γ ( 0 ) = q 0 , γ ( T ) = q 1 , q i M ,
and having the minimal sub-Riemannian length l ( γ ) min .
Remark 1.
It is easy to check that system (1) is invariant under parallel translations and rotations in the plane ( x , y ) . Due to this, without loss of generality, we can reduce the problem for an arbitrary q 0 = ( x 0 , y 0 , θ 0 , k 0 ) to the case q 0 = ( 0 , 0 , 0 , k 0 ) .

3. Complete Controllability of the System

The first question that arises when studying problems (1)–(3) is the existence of an admissible trajectory connecting boundary conditions (3). If for any q 0 , q 1 M the answer is positive, then the control system is called completely controllable. Let us investigate the controllability of system (1) using the technique of geometric control theory [20].
System (1) has the following form:
γ ˙ = u 1 X 1 + u 2 X 2 ,
where the vector fields near the controls are given by
X 1 = cos θ sin θ k 0 , X 2 = 0 0 0 1 .
We investigate the controllability using Chow–Rashevskii theorem. In our case, it suffices to check that the rank condition is satisfied. To do this, we calculate the following Lie brackets of the fields X i :
X 3 = [ X 1 , X 2 ] = 0 0 1 0 , X 4 = [ X 1 , X 3 ] = sin θ cos θ 0 0 .
For the matrix composed of the vector fields X 1 , , X 4 we have
det cos θ 0 0 sin θ sin θ 0 0 cos θ k 0 1 0 0 1 0 0 1 .
We conclude that the rank of the matrix is four, so the vector fields X i are linearly independent. Therefore, they define a basis of the tangent space T q M at every point q. Thus, we see that all of the conditions of the Chow–Rashevskii theorem are satisfied and we obtain the following result.
Theorem 1.
The control system (1) is completely controllable.
Remark 2.
Since the values u i are unbounded, condition (4) is equivalent to γ ˙ Δ γ = span ( X 1 ( γ ) , X 2 ( γ ) ) , where the family of planes Δ is called the distribution. Due to condition (5), the growth vector of the distribution Δ equals ( 2 , 3 , 4 ) . Such systems are called structures of Engel type.
Furthermore, a question of the existence of optimal trajectories arises: does there always exist an admissible trajectory satisfying conditions (3), on which the functional (2) reaches its minimum value? For our problems (1)–(3), the answer is positive. The existence of optimal trajectories is guaranteed by the Filippov theorem [20,21].

4. Pontryagin Maximum Principle

Before proceeding to the examination of extremal trajectories, let us reduce the problem under consideration to a simpler one. By virtue of Cauchy–Schwarz inequality, the original problem is equivalent to the problem of minimizing the action functional
J = 0 T u 1 2 ( t ) + u 2 2 ( t ) 2 d t min .
Apply to problems (1), (3), (6) a necessary condition of optimality—Pontryagin maximum principle (PMP)—we introduce the Pontryagin function:
h u ν = p , i = 1 2 u i X i + ν / 2 i = 1 2 u i 2 , p T * M , ν 0 .
PMP states that, if ( u ( t ) , q ( t ) ) , t [ 0 , T ] , is an optimal process, then the following conditions hold:
  • Hamiltonian system p ˙ = h u ν q , q ˙ = h u ν p ;
  • Maximum condition h u ( t ) ν ( p ( t ) , q ( t ) ) = max u R 2 h u ν ( p ( t ) , q ( t ) ) ;
  • Nontriviality condition ( p ( t ) , ν ) ( 0 , 0 ) t [ 0 , T ] .
Denote h i = p , X i . The Pontryagin function takes the form
h u ν = u 1 h 1 + u 2 h 2 + ν / 2 u 1 2 + u 2 2 .
In the formulation of PMP, without loss of generality, it suffices to consider two cases: ν = 0 , the abnormal case, and ν = 1 , the normal case. Next, we consider both cases in detail.

4.1. Abnormal Case ν = 0

The Pontryagin function is h u 0 = u 1 h 1 + u 2 h 2 . This is a linear function, unbounded when h 1 2 + h 2 2 0 . Thus, the maximum condition is satisfied if and only if h 1 = h 2 0 . The maximized Hamiltonian in this case is H = max u R 2 h u 0 = 0 .
Let a , b , c , d R be the components of the covector p in canonical Darboux coordinates. They change by the law p i ˙ = h u 0 q i .
The Hamiltonian system of PMP has the form
x ˙ = u 1 cos θ , y ˙ = u 1 sin θ , θ ˙ = u 1 k , k ˙ = u 2 , a ˙ = 0 , b ˙ = 0 , c ˙ = u 1 ( a sin θ + b cos θ ) , d ˙ = u 1 c .
The condition h 1 = h 2 0 implies
a cos θ + b sin θ + c k 0 , d 0 .
By virtue of system (7), the second identity implies d ˙ = u 1 c 0 , while differentiation of the first identity implies u 2 c 0 . Thus, we have
( u 1 2 + u 2 2 ) c 2 0 .
Since every admissible curve of positive length is a Lipschitz reparameterization of an arclength parameterized admissible one (see ([22], Lemma 3.16)), one can chose the natural parameterization u 1 2 + u 2 2 = 1 for every trajectory of system (1) that is not a fixed point. Thus, identity (9) is equivalent to c 0 on the intervals of time where the trajectory is not a fixed point.
The identity c 0 implies c ˙ 0 . From the third equation of system (7), it follows that u 1 ( a sin θ + b cos θ ) 0 . On the other hand, the first identity of system (8) takes the form a cos θ + b sin θ 0 . Thus, we have u 1 2 ( a 2 + b 2 ) 0 . Note that, if a 2 + b 2 = 0 , then a = b = c = d = 0 , which contradicts the nontriviality condition in PMP. Thus, the case u 1 0 remains to be considered.
In this case, the covector p is constant and nonzero. We obtain that the abnormal extremals have the form γ ( t ) = ( 0 , 0 , 0 , k 0 + U ( t ) ) , where U ( t ) = 0 t u 2 ( τ ) d τ , and u 2 ( t ) is any real-valued L ( 0 , T ) function. It is easy to see that the trajectory is not optimal if u 2 changes its sign. This holds since, when u 2 changes its sign, the trajectory is followed in opposite directions. Choosing the natural parameterization u 1 2 + u 2 2 = 1 on optimal trajectories, we obtain u 2 = ± 1 . This results in the optimal trajectory γ ( t ) = ( 0 , 0 , 0 , k 0 ± t ) . Thus, we obtain the following result.
Theorem 2.
The abnormal extremal trajectories in problems (1)–(3) have the form γ ( t ) = ( 0 , 0 , 0 , k 0 + 0 t u 2 ( τ ) d τ ) , where u 2 ( t ) is any real-valued L ( 0 , T ) function. Naturally parameterized abnormal optimal trajectories have the form γ ( t ) = ( 0 , 0 , 0 , k 0 ± t ) .

4.2. Normal Case ν = 1

The Pontryagin function h u 1 takes the form h u 1 = u 1 h 1 + u 2 h 2 u 1 2 + u 2 2 / 2 . The maximum condition gives the expression for the extremal controls:
h u 1 u i = h i u i = 0 u i = h i .
The maximized Hamiltonian H = max u R 2 h u 1 takes the form
H = h 1 2 + h 2 2 2 .
By definition h i = p , X i , i = 1 , , 4 we have
h 1 = a cos θ + b sin θ + c k , h 2 = d , h 3 = c , h 4 = a sin θ + b cos θ . a = ( h 1 + k h 3 ) cos θ h 4 sin θ , b = h 4 cos θ + ( h 1 + k h 3 ) sin θ , c = h 3 , d = h 2 .
The Hamiltonian system in canonical coordinates has the form
x ˙ = ( a cos θ + b sin θ + c k ) cos θ , y ˙ = ( a cos θ + b sin θ + c k ) sin θ , θ ˙ = k ( a cos θ + b sin θ + c k ) , k ˙ = d , a ˙ = 0 , b ˙ = 0 , c ˙ = ( b cos θ a sin θ ) ( k c + a cos θ + b sin θ ) , d ˙ = c ( c k + a cos θ + b sin θ ) .
We rewrite this system in the coordinates h i :
x ˙ = h 1 cos θ , y ˙ = h 1 sin θ , θ ˙ = h 1 k , k ˙ = h 2 , h ˙ 1 = h 2 h 3 , h ˙ 2 = h 1 h 3 , h ˙ 3 = h 1 h 4 , h ˙ 4 = k h 1 ( k h 3 + h 1 ) .
By choosing the natural parameterization u 1 2 + u 2 2 = 1 on the extremal trajectories, one fixes the level surface of the Hamiltonian H = ( h 1 2 + h 2 2 ) / 2 = 1 / 2 . A polar angle α S 1 is introduce into the plane ( h 1 , h 2 ) :
h 1 = cos α , h 2 = sin α .
By rewriting Hamiltonian system (10) we obtain the following result.
Theorem 3.
The naturally parameterized normal extremal trajectories in problems (1)–(3) are solutions to the system
x ˙ = cos α cos θ , y ˙ = cos α sin θ , θ ˙ = k cos α , k ˙ = sin α , α ˙ = h 3 , h ˙ 3 = h 4 cos α , h ˙ 4 = k cos α ( k h 3 + cos α ) .
The question of the existence of an analytic expression for extremal trajectories arises: is the Hamiltonian system integrable? To prove the Liouville integrability of system (10), it suffices to find four functionally independent first integrals in involution. One such integral is the Hamiltonian H, and two more first integrals a and b follow directly from the representation of the Hamiltonian system in canonical coordinates. These three integrals are functionally independent and are in involution. The existence of the remaining first integral is under investigation.
Numerical experiments have been carried out, indicating the presence of the fourth integral. To this end, we consider the four-dimensional system, which is decoupled from the rest of the variables:
α ˙ = h 3 , h ˙ 3 = h 4 cos α , h ˙ 4 = k cos α ( k h 3 + cos α ) , k ˙ = sin α .

Numerical Simulations

The existence of the first integral of system (12) was studied using the Poincare map. The method consists of three steps. In the first step, one needs to find a periodic trajectory of the system. The second step consists of the construction of a manifold transversal to the flow of the system in a neighborhood of a point of the periodic trajectory. In the third step, small perturbations of the initial point are considered and the corresponding trajectories are computed numerically until they intersect the transversal submanifold several times (N times). Such points are called orbits of the Poincare map. If the dynamics exhibit chaotic behavior, the trajectories strongly diverge under a small perturbation of the initial point and the orbits form a set of points chaotically distributed in the transversal submanifold. In the case of an integrable system, an orbit forms a set of points lying in a submanifold of a smaller dimension in the transversal manifold. Next, we perform this method for system (12).
  • Step 1.
There exists a one-parameter family of periodic trajectories
α ( t ) = π 2 + t h 3 ( 0 ) , h 3 ( t ) = h 3 ( 0 ) , h 4 ( t ) = 0 , k ( t ) = sin ( t h 3 ( 0 ) ) h 3 ( 0 ) ,
obtained by changing the initial value h 3 ( 0 ) . In this case, the period is T = | 2 π h 3 ( 0 ) | . With such initial values, the point ( h 3 , h 4 ) = ( h 3 ( 0 ) , 0 ) is fixed. For the next steps, we choose a periodic trajectory passing through the initial point α ( 0 ) = π 2 , h 3 ( 0 ) = 1 , h 4 ( 0 ) = 0 , k ( 0 ) = 0 .
  • Step 2.
It can be checked that the hyperspace k = 0 is transversal to the periodic trajectory at the initial point. Indeed, the tangent vector to the periodic trajectory at the initial point is ( 0 , 0 , 0 , 1 ) , which is orthogonal to the hyperspace k = 0 .
  • Step 3.
New initial points were chosen from a small neighborhood of the initial point of the periodic trajectory. The trajectories departing from these points were computed. Points of intersections of every such trajectory with the transversal hyperspace were computed (for the first time, the second time, , the N-th time).
In Figure 2, the orbits of the Poincare map in the space ( α , h 3 , h 4 ) are shown. The red dot corresponds to the orbit of a periodic trajectory with initial value α ( 0 ) = π 2 , h 3 ( 0 ) = 1 , h 4 ( 0 ) = 0 . For the rest of the trajectories, the following parameters were used. For the orange trajectory, α ( 0 ) = 1.56 , h 3 ( 0 ) = 0.94 , h 4 ( 0 ) = 0.02 . For the green trajectory, α ( 0 ) = 1.55 , h 3 ( 0 ) = 1.06 , h 4 ( 0 ) = 0.05 . For the black trajectory, α ( 0 ) = 1.6 , h 3 ( 0 ) = 1.14 , h 4 ( 0 ) = 0.02 . For the blue trajectory, α ( 0 ) = 1.58 , h 3 ( 0 ) = 1.24 , h 4 ( 0 ) = 0.01 . The number of iterations of the Poincare map for all trajectories was chosen as N = 1000 .
The trajectories were computed by numerical integration (NDSolve in Wolfram Mathematica). The instances when the trajectory intersects the transversal hyperspace were determined using the numerical solution (FindRoot) of the equation k ( t ) = 0 along the given trajectory with initial approximation t = | 2 π h 3 ( 0 ) | at each iteration i = 1 , , N .
It is remarkable that the points of the Poincare map fill continuous closed curves. This indicates that, for a small perturbation of the initial conditions, the trajectories remain close over a long time interval. This situation arises when the system is integrable (at least in some domain). Thus, the presented numerical experiments indicate the presence of a fourth independent first integral, which results in Liouville integrability of the Hamiltonian system (11).

5. Conclusions

In this paper, we consider the problem of sub-Riemannian geodesics in the four-dimensional manifold M = R 2 × SO ( 2 ) × R . This problem arises when modelling the mechanism of occluded contours completion in the extended model of the visual cortex by J. Petitot, G. Citti, and A. Sarti. The extension of the classical 3D model is performed by taking the curvature of the contours into account.
The following main results are obtained in the article. Complete controllability of the system and existence of optimal controls are proven. The Hamiltonian system of PMP is derived. The explicit parametrization of abnormal trajectories is found. In the normal case, three functionally independent first integrals are found. A hypothesis is formulated, confirmed by numerical experiments: the normal Hamiltonian system of PMP is Liouville integrable.
Integrability is an important property of the model. This property indicates the adequateness of the model (the visual signal is processed in a deterministic way). It ensures the stability of numerical methods of integration, which appear in brain-inspired image processing algorithms (such as image enchantment, inpainting, and salient curve extraction).

Author Contributions

Conceptualization, A.M.; methodology, A.M.; software, I.G. and A.M.; validation, I.G. and A.M.; formal analysis, I.G. and A.M.; investigation, I.G. and A.M.; resources, I.G. and A.M.; data curation, A.M.; writing—original draft preparation, I.G. and A.M.; writing—review and editing, I.G. and A.M.; visualization, I.G. and A.M.; supervision, A.M.; project administration, A.M.; funding acquisition, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

The work by Alexey Mashtakov was supported by the Russian Science Foundation under grant 17-11-01387-P and performed in Ailamazyan Program Systems Institute of Russian Academy of Sciences. The work of Ivan Galyaev was supported by the Russian Science Foundation under grant 21-71-20034 and performed in V.A. Trapeznikov Institute of Control Sciences of Russian Academy of Sciences. Section 1, Section 2, Section 4.2 and Section 5 of the paper were written by A. Mashtakov. Section 3 and Section 4.1 of the paper were written by I. Galyaev.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

https://notebookarchive.org//2021-10-1uutypp (accessed on 15 December 2021).

Acknowledgments

The authors are grateful to Yu.L. Sachkov for many fruitful suggestions and valuable comments on the work.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
PMPPontryagin maximum principle

References

  1. Ter Haar Romeny, B.M. Front-End Vision and Multi-Scale Image Analysis. In Multi-Scale Computer Vision Theory and Applications, Written in Mathematics; Computational Imaging and Vision; Springer: Dordrecht, The Netherlands, 2003; Volume 27. [Google Scholar]
  2. Hubel, D.H.; Wiesel, T.N. Receptive fields, binocular interaction and functional architecture in the cat’s visual cortex. J. Physiol. 1962, 160, 106–154. [Google Scholar] [CrossRef] [PubMed]
  3. Petitot, J. The neurogeometry of pinwheels as a sub-Riemannian contact structure. J. Physiol. Paris 2003, 97, 265–309. [Google Scholar] [CrossRef]
  4. Citti, G.; Sarti, A. A cortical based model of perceptual completion in the roto-translation space. J. Math. Imaging Vis. 2006, 24, 307–326. [Google Scholar] [CrossRef]
  5. Boscain, U.; Gauthier, J.-P.; Chertovskih, R.; Remizov, A. Hypoelliptic diffusion and human vision: A semidiscrete new twist. Siam J. Imaging Sci. 2014, 7, 669–695. [Google Scholar] [CrossRef]
  6. Boscain, U.V.; Chertovskih, R.; Gauthier, J.-P.; Prandi, D.; Remizov, A. Highly corrupted image inpainting through hypoelliptic diffusion. J. Math. Imaging Vis. 2018, 60, 1231–1245. [Google Scholar] [CrossRef] [Green Version]
  7. Duits, R.; Franken, E. Left-invariant parabolic evolutions on SE(2) and contour enhancement via invertible orientation scores Part I: Linear left-invariant diffusion equations on SE(2). Q. Appl. Math. 2010, 68, 255–292. [Google Scholar] [CrossRef] [Green Version]
  8. Sachkov, Y.L. Cut locus and optimal synthesis in the sub-Riemannian problem on the group of motions of a plane. ESAIM Control. Optim. Calc. Var. 2010, 17, 293–321. [Google Scholar] [CrossRef] [Green Version]
  9. Mashtakov, A.P.; Ardentov, A.A.; Sachkov, Y.L. Parallel Algorithm and Software for Image Inpainting via Sub-Riemannian Minimizers on the Group of Rototranslationse. Numer. Math. Theory Methods Appl. 2013, 65, 55–77. [Google Scholar]
  10. Franceschiello, B.; Mashtakov, A.; Citti, G.; Sarti, A. Geometrical optical illusion via sub-Riemannian geodesics in the roto-translation group. Differ. Geom. Its Appl. 2019, 65, 55–77. [Google Scholar] [CrossRef] [Green Version]
  11. Duits, R.; Felsberg, M.; Granlund, G.; Romeny, B. Image Analysis and Reconstruction using a Wavelet Transform Constructed from a Reducible Representation of the Euclidean Motion Group. Int. J. Comput. Vis. 2007, 72, 79–102. [Google Scholar] [CrossRef] [Green Version]
  12. Bekkers, E.J.; Duits, R.; Mashtakov, A.; Sanguinetti, G.R. A PDE Approach to Data-driven Sub-Riemannian Geodesics in SE(2). Siam J. Imaging Sci. 2015, 8, 2740–2770. [Google Scholar] [CrossRef] [Green Version]
  13. Mashtakov, A.; Duits, R.; Sachkov, Y.; Bekkers, E.J.; Beschastnyi, I. Tracking of Lines in Spherical Images via Sub-Riemannian Geodesics in SO(3). J. Math. Imaging Vis. 2017, 58, 239–264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Duits, R.; Ghosh, A.; Dela Haije, T.; Mashtakov, A. On sub-Riemannian geodesics in SE(3) whose spatial projections do not have cusps. J. Dyn. Control Syst. 2016, 22, 771–805. [Google Scholar] [CrossRef] [Green Version]
  15. Yue, X.; Robert, S.; Ungerleider, L.G. Curvature processing in human visual cortical areas. NeuroImage 2020, 222, 117295. [Google Scholar] [CrossRef] [PubMed]
  16. Zucker, S.W. The computational connection in vision: Early orientation selection. Behav. Res. Methods Instrum. Comput. 1986, 18, 608–617. [Google Scholar] [CrossRef] [Green Version]
  17. Petitot, J. Landmarks for Neurogeometry. In Neuromathematics of Vision; Citti, G., Sarti, A., Eds.; Spring: Berlin/Heidelberg, Geramny, 2014. [Google Scholar]
  18. Citti, G.; Giovannardi, G.; Ritore, M.; Sarti, A. Submanifolds of Fixed Degree in Graded Manifolds for Perceptual Completion. In Geometric Science of Information, Proceedings of the 5th International Conference, GSI 2021, Paris, France, 21–23 July 2021; Springer Nature Switzerland AG: Cham, Switzerland, 2021; Volume 12829. [Google Scholar]
  19. Lokutsievskiy, L.V.; Sachkov, Y.L. Liouville integrability of sub-Riemannian problems on Carnot groups of step 4 or greater. Sb. Math. 2018, 209, 672–713. [Google Scholar] [CrossRef]
  20. Agrachev, A.A.; Sachkov, Y.L. Control Theory from the Geometric Viewpoint; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  21. Zelikin, M.I. Optimal Control and Variational Calculus; URSS: Moscow, Russia, 2004. [Google Scholar]
  22. Agrachev, A.; Barilari, D.; Boscain, U. A Comprehensive Introduction to Sub-Riemannian Geometry; Cambridge Studies in Advanced Mathematics; Cambridge University Press: Cambridge, UK, 2019. [Google Scholar]
Figure 1. In the four-dimensional model of the visual cortex, an occluded contour is completed via the planar projection of a sub-Riemannian length-minimizer in the space R 2 × SO ( 2 ) × R = M ( x , y , θ , k ) of positions, orientations, and curvatures. In the left column, we show an example of the image with partially occluded contours and the complete image. In the right column, we show a trajectory that satisfies the given boundary conditions. The curvature is visualized as its reciprocal—the radius of the osculating circle.
Figure 1. In the four-dimensional model of the visual cortex, an occluded contour is completed via the planar projection of a sub-Riemannian length-minimizer in the space R 2 × SO ( 2 ) × R = M ( x , y , θ , k ) of positions, orientations, and curvatures. In the left column, we show an example of the image with partially occluded contours and the complete image. In the right column, we show a trajectory that satisfies the given boundary conditions. The curvature is visualized as its reciprocal—the radius of the osculating circle.
Jimaging 07 00277 g001
Figure 2. Orbits of the Poincare map in the space ( α , h 3 , h 4 ) are formed by intersection points of the transversal hyperspace k = 0 with trajectories close to the periodic one (red dot). Different orbits are depicted in different colors. Starting points are indicated for each trajectory.
Figure 2. Orbits of the Poincare map in the space ( α , h 3 , h 4 ) are formed by intersection points of the transversal hyperspace k = 0 with trajectories close to the periodic one (red dot). Different orbits are depicted in different colors. Starting points are indicated for each trajectory.
Jimaging 07 00277 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Galyaev, I.; Mashtakov, A. Liouville Integrability in a Four-Dimensional Model of the Visual Cortex. J. Imaging 2021, 7, 277. https://doi.org/10.3390/jimaging7120277

AMA Style

Galyaev I, Mashtakov A. Liouville Integrability in a Four-Dimensional Model of the Visual Cortex. Journal of Imaging. 2021; 7(12):277. https://doi.org/10.3390/jimaging7120277

Chicago/Turabian Style

Galyaev, Ivan, and Alexey Mashtakov. 2021. "Liouville Integrability in a Four-Dimensional Model of the Visual Cortex" Journal of Imaging 7, no. 12: 277. https://doi.org/10.3390/jimaging7120277

APA Style

Galyaev, I., & Mashtakov, A. (2021). Liouville Integrability in a Four-Dimensional Model of the Visual Cortex. Journal of Imaging, 7(12), 277. https://doi.org/10.3390/jimaging7120277

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop