Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA class of fast geodesic shooting algorithms for template matching and its applications via the $N$-particle system of the Euler-Poincare equations
66
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The Euler-Poincare (EP) equations describe the geodesic motion on the diffeomorphism group. For template matching (template deformation), the Euler-Lagrangian equation, arising from minimizing an energy function, falls into the Euler-Poincare theory and can be recast into the EP equations. By casting the EP equations in the Lagrangian (or characteristics) form, we formulate the equations as a finite dimensional particle system. The evolution of this particle system describes the geodesic motion of landmark points on a Riemann manifold. In this paper we present a class of novel algorithms that take advantage of the structure of the particle system to achieve a fast matching process between the reference and the target templates. The strong suit of the proposed algorithms includes (1) the efficient feedback control iteration, which allows one to find the initial velocity field for driving the deformation from the reference template to the target one, (2) the use of the conical kernel in the particle system, which limits the interaction between particles and thus accelerates the convergence, and (3) the availability of the implementation of fast-multipole method for solving the particle system, which could reduce the computational cost from $O(N^2)$ to $O(Nlog N)$, where $N$ is the number of particles. The convergence properties of the proposed algorithms are analyzed. Finally, we present several examples for both exact and inexact matchings, and numerically analyze the iterative process to illustrate the efficiency and the robustness of the proposed algorithms.
rate research
Read More
We study a class of partial differential equations (PDEs) in the family of the so-called Euler-Poincare differential systems, with the aim of developing a foundation for numerical algorithms of their solutions. This requires particular attention to the mathematical properties of this system when the associated class of elliptic operators possesses non-smooth kernels. By casting the system in its Lagrangian (or characteristics) form, we first formulate a particles system algorithm in free space with homogeneous Dirichlet boundary conditions for the evolving fields. We next examine the deformation of the system when non-homogeneous constant stream boundary conditions are assumed. We show how this simple change at the boundary deeply affects the nature of the evolution, from hyperbolic-like to dispersive with a non-trivial dispersion relation, and examine the potentially regularizing properties of singular kernels offered by this deformation. From the particle algorithm viewpoint, kernel singularities affect the existence and uniqueness of solutions to the corresponding ordinary differential equations systems. We illustrate this with the case when the operator kernel assumes a conical shape over the spatial variables, and examine in detail two-particle dynamics under the resulting lack of Lipschitz-continuity. Curiously, we find that for the conically-shaped kernels the motion of the related two-dimensional waves can become completely integrable under appropriate initial data. This reduction projects the two-dimensional system to the one-dimensional completely integrable Shallow-Water equation [Camassa, R. and Holm, D. D., Phys. Rev. Lett., 71, 1961-1964, 1993], while retaining the full dependence on two spatial dimensions for the single channel solutions.
In this paper, we propose a numerical method to solve the classic $L^2$-optimal transport problem. Our algorithm is based on use of multiple shooting, in combination with a continuation procedure, to solve the boundary value problem associated to the transport problem. We exploit the viewpoint of Wasserstein Hamiltonian flow with initial and target densities, and our method is designed to retain the underlying Hamiltonian structure. Several numerical examples are presented to illustrate the performance of the method.
We study Euler-Poincare systems (i.e., the Lagrangian analogue of Lie-Poisson Hamiltonian systems) defined on semidirect product Lie algebras. We first give a derivation of the Euler-Poincare equations for a parameter dependent Lagrangian by using a variational principle of Lagrange dAlembert type. Then we derive an abstract Kelvin-Noether theorem for these equations. We also explore their relation with the theory of Lie-Poisson Hamiltonian systems defined on the dual of a semidirect product Lie algebra. The Legendre transformation in such cases is often not invertible; so it does not produce a corresponding Euler-Poincare system on that Lie algebra. We avoid this potential difficulty by developing the theory of Euler-Poincare systems entirely within the Lagrangian framework. We apply the general theory to a number of known examples, including the heavy top, ideal compressible fluids and MHD. We also use this framework to derive higher dimensional Camassa-Holm equations, which have many potentially interesting analytical properties. These equations are Euler-Poincare equations for geodesics on diffeomorphism groups (in the sense of the Arnold program) but where the metric is H^1 rather than L^2.
We find necessary and sufficient conditions for the foliation defined by level sets of a function f(x_{1},...,x_{n}) to be totally geodesic in a torsion-free connection and apply them to find the conditions for d-webs of hypersurfaces to be geodesic, and in the case of flat connections, for d-webs (d > n) of hypersurfaces to be hyperplanar webs. These conditions are systems of generalized Euler equations, and for flat connections we give an explicit construction of their solutions.
The Green Nagdhi equations are frequently used as a model of the wave-like behaviour of the free surface of a fluid, or the interface between two homogeneous fluids of differing densities. Here we show that their multilayer extension arises naturally from a framework based on the Euler Poincare theory under an ansatz of columnar motion. The framework also extends to the travelling wave solutions of the equations. We present numerical solutions of the travelling wave problem in a number of flow regimes. We find that the free surface and multilayer waves can exhibit intriguing differences compared to the results of single layer or rigid lid models.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
