Surface conformal maps between genus-0 surfaces play important roles in applied mathematics and engineering, with applications in medical image analysis and computer graphics. Previous work (Gu and Yau in Commun Inf Syst 2(2):121–146, 2002) introduces a variational approach, where global conformal parameterization of genus-0 surfaces was addressed through minimizing the harmonic energy, with two weaknesses: its gradient descent iteration is slow, and its solutions contain undesired parameterization foldings when the underlying surface has long sharp features. In this paper, we propose an algorithm that significantly accelerates the harmonic energy minimization and a method that iteratively removes foldings by taking advantages of the weighted Laplace–Beltrami eigen-projection. Experimental results show that the proposed approaches compute genus-0 surface harmonic maps much faster than the existing algorithm in Gu and Yau (Commun Inf Syst 2(2):121–146, 2002) and the new results contain no foldings.

Surface parameterizations and registrations are important in computer graphics and imaging, where 1-1 correspondences between meshes are computed. In practice, surface maps are usually represented and stored as three-dimensional coordinates each vertex is mapped to, which often requires lots of memory. This causes inconvenience in data transmission and data storage. To tackle this problem, we propose an effective algorithm for compressing surface homeomorphisms using Fourier approximation of the Beltrami representation. The Beltrami representation is a complex-valued function defined on triangular faces of the surface mesh with supreme norm strictly less than 1. Under suitable normalization, there is a 1-1 correspondence between the set of surface homeomorphisms and the set of Beltrami representations. Hence, every bijective surface map is associated with a unique Beltrami representation. Conversely, given a Beltrami representation, the corresponding bijective surface map can be exactly reconstructed using the linear Beltrami solver introduced in this paper. Using the Beltrami representation, the surface homeomorphism can be easily compressed by Fourier approximation, without distorting the bijectivity of the map. The storage requirement can be effectively reduced, which is useful for many practical problems in computer graphics and imaging. In this paper, we propose applying the algorithm to texture map compression and video compression. With our proposed algorithm, the storage requirement for the texture properties of a textured surface can be significantly reduced. Our algorithm can further be applied to compressing motion vector fields for video compression, which effectively improves the compression ratio.

A numerical method for computing the extremal Teichmüller map between multiply-connected domains is presented. Given two multiply-connected domains, there exists a unique Teichmüller map (T-Map) between them minimizing the conformality distortion. The extremal T-Map can be considered as the ‘most conformal’ map between multiply-connected domains. In this paper, we propose an iterative algorithm to compute the extremal T-Map using the Beltrami holomorphic flow (BHF). The BHF procedure iteratively adjusts the initial map based on a sequence of Beltrami coefficients, which are complex-valued functions defined on the source domain. It produces a sequence of quasi-conformal maps, which converges to the T-Map minimizing the conformality distortion. We test our method on synthetic data together with real human face data. Results show that our algorithm computes the extremal T-Map between two multiply-connected domains of the same topology accurately and efficiently.

Registration, which aims to find an optimal 1-1 correspondence between shapes, is an important process in different research areas. Landmark-based surface registration has been widely studied to obtain a mapping between shapes that matches important features. Obtaining a unique and bijective surface registration that matches features consistently is generally challenging, especially when a large number of landmark constraints are enforced. This motivates us to search for a unique landmark matching surface diffeomorphism, which minimizes the local geometric distortion. For this purpose, we propose a special class of diffeomorphisms called the Teichmüller mappings (T-Maps). Under suitable conditions on the landmark constraints, a unique T-Map between two surfaces can be obtained, which minimizes the maximal conformality distortion. The conformality distortion measures how far the mapping deviates from a conformal mapping, and hence it measures the local geometric distortion. In this paper, we propose an efficient iterative algorithm, called the quasi-conformal (QC) iteration, to compute the T-Map. The basic idea is to represent the set of diffeomorphisms using Beltrami coefficients (BCs) and look for an optimal BC associated to the desired T-Map. The associated diffeomorphism can be efficiently reconstructed from the optimal BC using the linear Beltrami solver (LBS). Using BCs to represent diffeomorphisms guarantees the diffeomorphic property of the registration, even with very large deformation. Using our proposed method, the T-Map can be accurately and efficiently computed. The obtained registration is guaranteed to be bijective. The proposed algorithm can also be extended to compute T-Map with soft landmark constraints. We applied the proposed algorithm to real applications, such as brain landmark matching registration, constrained texture mapping, and human face registration. Experimental results shows that our method is both effective and efficient in computing a nonoverlap landmark matching registration with the least amount of conformality distortion.

Shape registration has a wide range of applications in geometric modeling, medical imaging, and computer vision. This paper focuses on the registration of the genus-3 vestibular systems and studies the geometric differences between the normal and Adolescent Idiopathic Scoliosis (AIS) groups. The non-trivial topology of the VS poses great technical challenges to the geometric analysis. To tackle these challenges, we present an effective and practical solution to register the vestibular systems. We first extract six geodesic landmarks for the VS, which are stable, intrinsic, and insensitive to the VS’s resolution and tessellation. Moreover, they are highly consistent regardless of the AIS and normal groups. The detected geodesic landmarks partition the VS into three patches, a topological annulus and two topological disks. For each pair of patches of the AIS subject and the control, we compute a bijective map using the holomorphic 1-form and harmonic map techniques. With a carefully designed boundary condition, the three individual maps can be glued in a seamless manner so that the resulting registration is a homeomorphism with exact landmark matching. Our method is robust, automatic and efficient. It takes only a few seconds on a low-end PC, which significantly outperforms the non-rigid ICP algorithm. We conducted a student’s t-test on the test data. Computational results show that using the mean curvature measure E_H, our method can distinguish the AIS subjects and the normal subjects.