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.

Planar polynomial curves have rational offset curves, if they are either Pythagorean-hodograph (PH) or indirect Pythagorean-hodograph (iPH) curves. In this paper, we derive an algebraic and two geometric characterizations for planar quartic iPH curves. The characterizations are given in terms of quantities related to the Bézier control polygon of the curve, and naturally extend to quartic and cubic PH and quadratic iPH curves.

In the past decades, many methods for computing conformal mesh parameterizations have been developed in response to demand of numerous applications in the field of geometry processing. Spectral conformal parameterization (SCP) (Mullen et al. in Proceedings of the symposium on geometry processing, SGP ’08. Eurographics Association, Aire-la-Ville, Switzerland, pp 1487–1494, 2008) is one of these methods used to compute a quality conformal parameterization based on the spectral techniques. SCP focuses on a generalized eigenvalue problem (GEP) $L_{C} \bf{f} = \lambda B \bf{f}$ whose eigenvector(s) associated with the smallest positive eigenvalue(s) provide the conformal parameterization result. This paper is devoted to studying a novel eigensolver for this GEP. Based on structures of the matrix pair $(L_{C}, B)$, we show that this GEP can be transformed into a small-scale compressed and deflated standard eigenvalue problem with a symmetric positive definite skew-Hamiltonian operator. We then propose a symmetric skew-Hamiltonian isotropic Lanczos algorithm (SHILA) to solve the reduced problem. Numerical experiments show that our compressed deflating technique can exclude the impact of convergence from the kernel of $L_{C}$ and transform the original problem to a more robust system. The novel SHILA method can effectively avoid the disturbance of duplicate eigenvalues. As a result, based on the spectral model of SCP, our numerical eigensolver can compute the conformal parameterization accurately and efficiently.

Shape analysis is a central problem in the field of computer vision. In 2D shape analysis, classification and recognition of objects from their observed silhouettes are extremely crucial but difficult. It usually involves an efficient representation of 2D shape space with a metric, so that its mathematical structure can be used for further analysis. Although the study of 2D simply-connected shapes has been subject to a corpus of literatures, the analysis of multiply-connected shapes is comparatively less studied. In this work, we propose a representation for general 2D multiply-connected domains with arbitrary topologies using conformal welding. A metric can be defined on the proposed representation space, which gives a metric to measure dissimilarities between objects. The main idea is to map the exterior and interior of the domain conformally to unit disks and circle domains (unit disk with several inner disks removed), using holomorphic 1-forms. A set of diffeomorphisms of the unit circle S1 can be obtained, which together with the conformal modules are used to define the shape signature. A shape distance between shape signatures can be defined to measure dissimilarities between shapes. We prove theoretically that the proposed shape signature uniquely determines the multiply-connected objects under suitable normalization. We also introduce a reconstruction algorithm to obtain shapes from their signatures. This completes our framework and allows us to move back and forth between shapes and signatures. With that, a morphing algorithm between shapes can be developed through the interpolation of the Beltrami coefficients associated with the signatures. Experiments have been carried out on shapes extracted from real images. Results demonstrate the efficacy of our proposed algorithm as a stable shape representation scheme.

Surface parameterizations have been widely applied to computer graphics and digital geometry processing. In this paper, we propose a novel stretch energy minimization (SEM) algorithm for the computation of equiareal parameterizations of simply connected open surfaces with very small area distortions and highly improved computational efficiencies. In addition, the existence of nontrivial limit points of the SEM algorithm is guaranteed under some mild assumptions of the mesh quality. Numerical experiments indicate that the accuracy, effectiveness, and robustness of the proposed SEM algorithm outperform the other state-of-the-art algorithms. Applications of the SEM on surface remeshing, registration and morphing for simply connected open surfaces are demonstrated thereafter. Thanks to the SEM algorithm, the computation for these applications can be carried out efficiently and reliably.