This paper presents a novel algorithm to obtain landmark-based genus-1 surface registration via a special class of quasi-conformal maps called the Teichmüller maps. Registering shapes with important features is an important process in medical imaging. However, it is challenging to obtain a unique and bijective genus-1 surface matching that satisfies the prescribed landmark constraints. In addition, as suggested by [11], conformal transformation provides the most natural way to describe the deformation or growth of anatomical structures. This motivates us to look for the unique mapping between genus-1 surfaces that matches the features while minimizing the maximal conformality distortion. Existence and uniqueness of such optimal diffeomorphism is theoretically guaranteed and is called the Teichmüller extremal mapping. In this work, we propose an iterative algorithm, called the Quasi-conformal (QC) iteration, to find the Teichmüller extremal mapping between the covering spaces of genus-1 surfaces. By representing the set of diffeomorphisms using Beltrami coefficients (BCs), we look for an optimal BC which corresponds to our desired diffeomorphism that matches prescribed features and satisfies the periodic boundary condition on the covering space. Numerical experiments show that our proposed algorithm is efficient and stable for registering genus-1 surfaces even with large amount of landmarks. We have also applied the algorithm on registering vertebral bones with prescribed feature curves, which demonstrates the usefulness of the proposed algorithm.

No abstract uploaded!

No abstract uploaded!

As V. I. Arnold observed in the 1960s, the Euler equations of incompressible fluid flow correspond formally to geodesic equations in a group of volume-preserving diffeomorphisms. Working in an Eulerian framework, we study incompressible flows of shapes as critical paths for action (kinetic energy) along transport paths constrained to be shape densities (characteristic functions). The formal geodesic equations for this problem are Euler equations for incompressible, inviscid potential flow of fluid with zero pressure and surface tension on the free boundary. The problem of minimizing this action exhibits an instability associated with microdroplet formation, with the following outcomes: Any two shapes of equal volume can be approximately connected by an Euler spray---a countable superposition of ellipsoidal geodesics. The infimum of the action is the Wasserstein distance squared, and is almost never attained except in dimension 1. Every Wasserstein geodesic between bounded densities of compact support provides a solution of the (compressible) pressureless Euler system that is a weak limit of (incompressible) Euler sprays. Each such Wasserstein geodesic is also the unique minimizer of a relaxed least-action principle for a two-fluid mixture theory corresponding to incompressible fluid mixed with vacuum.

We deal with a stationary problem of a reaction–diffusion system with a conservation law under the Neu-mann boundary condition. It is shown that the stationary problem turns to be the Euler–Lagrange equation of an energy functional with a mass constraint. When the domain is the finite interval (0, 1), we investigate the asymptotic profile of a strictly monotone minimizer of the energy as d, the ratio of the diffusion coef-ficient of the system, tends to zero. In view of a logarithmic function in the leading term of the potential, we get to a scaling parameter κsatisfying the relation ε:=√d=√logκ/κ2. The main result shows that a sequence of minimizers converges to a Dirac mass multiplied by the total mass and that by a scaling with κthe asymptotic profile exhibits a parabola in the nonvanishing region. We also prove the existence of an unstable monotone solution when the mass is small.