This paper presents a method of obtaining geometric registrations between high-genus (g≥1) surfaces. Surface registration between simple surfaces, such as simply connected open surfaces, has been well studied. However, very few works have been carried out for the registration of high-genus surfaces. The high-genus topology of the surface poses a great challenge for surface registration. A possible approach is to partition surfaces into simply connected patches and registration can be done in a patch-by-patch manner. Consistent cuts are required, which are usually difficult to obtain and prone to error. In this work, we propose an effective way to obtain geometric registration between high-genus surfaces without introducing consistent cuts. The key idea is to conformally parameterize the surface into its universal covering space, which is either the Euclidean plane or the hyperbolic disk embedded in R2. Registration can then be done on the universal covering space by iteratively minimizing a shape mismatching energy measuring the geometric dissimilarity between the two surfaces. The Beltrami coefficient of the mapping is considered and adjusted in order to control the bijectivity of the mappings in each iteration. Our proposed algorithm effectively computes a smooth registration between high-genus surfaces that matches geometric information as much as possible. The algorithm can also be applied to find a smooth registration minimizing any general energy functionals. Numerical experiments on high-genus surface data show that our proposed method is effective for registering high-genus surfaces with geometric matching. We also applied the method to register anatomical structures for medical imaging, which demonstrates the usefulness of the proposed algorithm.

This work proposes a novel framework for optimization in the constrained diffeomorphism space for deformable surface registration. First the diffeomorphism space is modeled as a special complex functional space on the source surface, the Beltrami coefficient space. The physically plausible constraints, in terms of feature landmarks and deformation types, define subspaces in the Beltrami coefficient space. Then the harmonic energy of the registration is minimized in the constrained subspaces. The minimization is achieved by alternating two steps: 1) optimization - diffuse the Beltrami coefficient, and 2) projection - first deform the conformal structure by the current Beltrami coefficient and then compose with a harmonic map from the deformed conformal structure to the target. The registration result is diffeomorphic, satisfies the physical landmark and deformation constraints, and minimizes the conformality distortion. Experiments on human facial surfaces demonstrate the efficiency and efficacy of the proposed registration framework.

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.

This paper presents two algorithms, based on conformal geometry, for the multi-scale representations of geometric shapes and surface morphing. A multi-scale surface representation aims to describe a 3D shape at different levels of geometric detail, which allows analyzing or editing surfaces at the global or local scales effectively. Surface morphing refers to the process of interpolating between two geometric shapes, which has been widely applied to estimate or analyze deformations in computer graphics, computer vision and medical imaging. In this work, we propose two geometric models for surface morphing and multi-scale representation for 3D surfaces. The basic idea is to represent a 3D surface by its mean curvature function, H, and conformal factor function λ, which uniquely determine the geometry of the surface according to Riemann surface theory. Once we have the (λ,H) parameterization of the surface, post-processing of the surface can be done directly on the conformal parameter domain. In particular, the problem of multi-scale representations of shapes can be reduced to the signal filtering on the  and H parameters. On the other hand, the surface morphing problem can be transformed to an interpolation process of two sets of (λ,H) parameters. We test the proposed algorithms on 3D human face data and MRI-derived brain surfaces. Experimental results show that our proposed methods can effectively obtain multi-scale surface representations and give natural surface morphing results.

Detecting deformities on objects is a typical topic in shape analysis, and has much applications such as abnormalities detection in medical imaging (e.g. growth of tumor, spread of cancer). While many algorithms are already well-established in a 2-dimensional case when the object is indeed a surface, a model that still performs well in the general n-dimensional case is still missing. It is our goal in this paper to complete this missing piece, by introducing an indicator in order to effectively distinguish between normal and abnormal deformities. The proposed framework is closely related to the classic 2-dimensional conformal geometry and quasi-conformal geometry. In this work, we model abnormal deformations by anisotropic deformations. Given any two objects of the same dimension (with landmark constraints in between), we define the Anisotropic Indicator, a locally defined real-valued function on the original object, which demonstrates the abnormalities in the deformation between them. Both global and local features about the abnormalities between the two objects can be tracked by analyzing the indicator. We tested the algorithm by detecting deformations on synthetic data and real data, and results show that our algorithm can detect deformations of different types and degrees.

Registration, which aims to find an optimal one-to-one correspondence between different data, is an important problem in various fields. This problem is especially challenging when large deformations occur. In this paper, we present a novel algorithm to obtain diffeomorphic image or surface registrations with large deformations via quasi-conformal maps. The basic idea is to minimize an energy functional involving a Beltrami coefficient term, which measures the distortion of the quasi-conformal map. The Beltrami coefficient effectively controls the bijectivity and smoothness of the registration, even with very large deformations. Using the proposed algorithm, landmark-based registration between images or surfaces can be effectively computed. The obtained registration is guaranteed to be diffeomorphic (1-1 and onto), even with a large deformation or large number of landmark constraints. The proposed algorithm can also be combined with matching intensity (such as image intensity or surface curvature) to improve the accuracy of the registration. Experiments have been carried out on both synthetic and real data. Results demonstrate the efficacy of the proposed algorithm to obtain diffeomorphic registration between images or surfaces.

Finding surface mappings with least distortion arises from many applications in various fields. Extremal Teichmüller maps are surface mappings with least conformality distortion. The existence and uniqueness of the extremal Teichmüller map between Riemann surfaces of finite type are theoretically guaranteed (see Fletcher and Markovic, Quasiconformal maps and Teichmüller theory, Oxford Graduate Texts in Math., vol. 11, Oxford University Press, Oxford, 2007). Recently, a simple iterative algorithm for computing the Teichmüller maps between connected Riemann surfaces with given boundary value was proposed by Lui, Lam, Yau, and Gu in Teichmüller extremal mapping and its applications to landmark matching registration, arXiv:1211.2569. Numerical results were reported in the paper to show the effectiveness of the algorithm. The method was successfully applied to landmarkmatching registration. The purpose of this paper is to prove the iterative algorithm proposed in loc. cit., indeed converges.

A fresh framework for mesh optimization, the Filtered Hooke's Optimization, is proposed. With the notion of the elasticity theory, the Hooke's Optimization is developed by modifying the Hooke's law, in which an elastic force is simulated on the edges of a mesh so that adjacent vertices are either attracted to each other or repelled from each other, so as to regularize the mesh in terms of triangulation. A normal torque force is acted on vertices to guarantee smoothness of the surface. In addition, a filtering scheme, called the Newtonian Filtering, is proposed as a supplementary tool for the proposed Hooke's Optimization to preserve the geometry of the mesh. Numerical simulations on meshes with different geometry indicate an impressive performance of our proposed framework to significantly improve the mesh triangulation without noteworthy distortions of the mesh geometry.

Finding a meaningful 1–1 correspondence between different data, such as images or surface data, has important applications in various fields. It involves the optimization of certain energy functionals over the space of all diffeomorphisms. This type of optimization problems (called the diffeomorphism optimization problems, DOPs) is especially challenging, since the bijectivity of the mapping has to be ensured. Recently, a method, called the Beltrami holomorphic flow (BHF), has been proposed to solve the DOP using quasi-conformal theories (Lui et al. in J Sci Comput 50(3):557–585, 2012). The optimization problem is formulated over the space of Beltrami coefficients (BCs), instead of the space of all diffeomorphisms. BHF iteratively finds a sequence of BCs associated with a sequence of diffeomorphisms, using the gradient descent method, to minimize the energy functional. The use of BCs effectively controls the smoothness and bijectivity of the mapping, and hence makes it easier to handle the constrained optimization problem. However, the algorithm is computationally expensive. In this paper, we propose an efficient splitting algorithm, based on the classical alternating direction method of multiplier (ADMM), to solve the DOP. The basic idea is to split the energy functional into two energy terms: one involves the BC whereas the other involves the quasi-conformal map. Alternating minimization scheme can then be applied to minimize the energy functional. The proposed method significantly speeds up the previous BHF approach. It also extends the previous BHF algorithm to Riemann surfaces of arbitrary topologies, such as multiply-connected shapes. Experiments have been carried out on synthetic together with real surface data, which demonstrate the efficiency and efficacy of the proposed algorithm to solve the DOP.

Surface registration between cortical surfaces is crucial in medical imaging for performing systematic comparisons between brains. Landmark-matching registration that matches anatomical features, called the sulcal landmarks, is often required to obtain a meaningful 1-1 correspondence between brain surfaces. This is commonly done by parameterizing the surface onto a simple parameter domain, such as the unit sphere, in which the sulcal landmarks are consistently aligned. Landmark matching surface registration can then be obtained from the landmark aligned parameterizations. For genus-0 closed brain surfaces, the optimized spherical harmonic parameterization, which aligns landmarks to consistent locations on the sphere, has been widely used. This approach is limited by the loss of bijectivity under large deformations and the slow computation. In this paper, we propose FLASH, a fast algorithm to compute the optimized spherical harmonic parameterization with consistent landmark alignment. This is achieved by formulating the optimization problem to ℂ and thereby linearizing the problem. Errors introduced near the pole are corrected using quasiconformal theories. Also, by adjusting the Beltrami differential of the mapping, a diffeomorphic (1-1, onto) spherical parameterization can be effectively obtained. The proposed algorithm has been tested on 38 human brain surfaces. Experimental results demonstrate that the computation of the landmark aligned spherical harmonic parameterization is significantly accelerated using the proposed algorithm.