Automatic computation of surface correspondence via harmonic map is an active research field in computer vision, computer graphics and computational geometry. It may help document and understand physical and biological phenomena and also has broad applications in biometrics, medical imaging and motion capture industries. Although numerous studies have been devoted to harmonic map research, limited progress has been made to compute a diffeomorphic harmonic map on general topology surfaces with landmark constraints. This work conquers this problem by changing the Riemannian metric on the target surface to a hyperbolic metric so that the harmonic mapping is guaranteed to be a diffeomorphism under landmark constraints. The computational algorithms are based on Ricci flow and nonlinear heat diffusion methods. The approach is general and robust. We employ our algorithm to study the

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.

We propose a new method to obtain landmark-matching transformations between n-dimensional Euclidean spaces with large deformations. Given a set of feature correspondences, our algorithm searches for an optimal folding-free mapping that satisfies the prescribed landmark constraints. The standard conformality distortion defined for mappings between 2-dimensional spaces is first generalized to the n-dimensional conformality distortion K(f) for a mapping f between n-dimensional Euclidean spaces (n≥3). We then propose a variational model involving K(f) to tackle the landmark-matching problem in higher dimensional spaces. The generalized conformality term K(f) enforces the bijectivity of the optimized mapping and minimizes its local geometric distortions even with large deformations. Another challenge is the high computational cost of the proposed model. To tackle this, we have also proposed a numerical method to solve the optimization problem more efficiently. Alternating direction method with multiplier (ADMM) is applied to split the optimization problem into two subproblems. Preconditioned conjugate gradient method with multi-grid preconditioner is applied to solve one of the sub-problems, while a fixed-point iteration is proposed to solve another subproblem. Experiments have been carried out on both synthetic examples and lung CT images to compute the diffeomorphic landmark-matching transformation with different landmark constraints. Results show the efficacy of our proposed model to obtain a folding-free landmark-matching transformation between n-dimensional spaces with large deformations.

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.

Surface registration, which transforms different sets of surface data into one common reference space, is an important process which allows us to compare or integrate the surface data effectively. If a nonrigid transformation is required, surface registration is commonly done by parameterizing the surfaces onto a simple parameter domain, such as the unit square or sphere. In this work, we are interested in looking for meaningful registrations between surfaces through parameterizations, using prior features in the form of landmark curves on the surfaces. In particular, we generate optimized conformal parameterizations which match landmark curves exactly with shape-based correspondences between them. We propose a variational method to minimize a compound energy functional that measures the harmonic energy of the parameterization maps and the shape dissimilarity between mapped points on the landmark curves. The novelty is that the computed maps are guaranteed to align the landmark features consistently and give a shape-based diffeomorphism between the landmark curves. We achieve this by intrinsically modeling our search space of maps as flows of smooth vector fields that do not flow across the landmark curves. By using the local surface geometry on the curves to define a shape measure, we compute registrations that ensure consistent correspondences between anatomical features. We test our algorithm on synthetic surface data. An application of our model to medical imaging research is shown, using experiments on brain cortical surfaces, with anatomical (sulcal) landmarks delineated, which show that our computed maps give a shape-based alignment of the sulcal curves without significantly impairing conformality. This ensures correct averaging and comparison of data across subjects.