The manipulation of surface homeomorphisms is an important aspect in three-dimensional modeling and surface processing. Every homeomorphic surface map can be considered as a quasi-conformal map, with its local nonconformal distortion given by its Beltrami differential. As a generalization of conformal maps, quasi-conformal maps are of great interest in mathematical study and real applications. Efficient and accurate computational construction of desirable quasi-conformal maps between general surfaces is crucial. However, in the literature we have reviewed, all existing computational works on construction of quasi-conformal maps to or from a compact domain require global parametrization onto the plane and are difficult to directly apply to maps between arbitrary surfaces. This work fills the gap by proposing to compute quasi-conformal homeomorphisms between arbitrary Riemann surfaces using discrete Beltrami flow, which is a vector field corresponding to the adjustment to the intrinsic Beltrami differential of the map. The vector field is defined by a partial differential equation in a local conformal coordinate. Based on this formulation and a composition formula, we can compute the Beltrami flow of any homeomorphism adjustment as a vector field on the target domain defined from the source domain, with appropriate boundary conditions and correspondences. Numerical tests show that our method provides a robust and efficient way of adjusting surface homeomorphisms. It is also insensitive to surface representation and has no limitation to the classes of surfaces that can be processed. Extensive numerical examples will be shown.

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.

In this paper, we propose a general framework for approximating differential operator directly on point clouds and use it for geometric understanding on them. The discrete approximation of differential operator on the underlying manifold represented by point clouds is based only on local approximation using nearest neighbors, which is simple, efficient and accurate. This allows us to extract the complete local geometry, solve partial differential equations and perform intrinsic calculations on surfaces. Since no mesh or parametrization is needed, our method can work with point clouds in any dimensions or co-dimensions or even with variable dimensions. The computation complexity scaled well with the number of points and the intrinsic dimensions (rather than the embedded dimensions). We use this method to define the Laplace-Beltrami (LB) operator on point clouds, which links local and global information together. With this operator, we propose a few key applications essential to geometric understanding for point clouds, including the computation of LB eigenvalues and eigenfunctions, the extraction of skeletons from point clouds, and the extraction of conformal structures from point clouds.

This paper proposes a novel approach for extracting two intrinsic feature curves on hippocampal (HC) surfaces. The hippocampus is a key target of study in medical imaging, as it degenerates in conditions such as epilepsy and Alzheimer's disease (AD), but its structure is complex. To facilitate HC morphometry, we generate two intrinsic feature curves that describe their global geometries. For example, the separation of them captures thickness changes in HC surfaces, which can be used to effectively measure HC atrophy found in patients with AD. They also separate HC surfaces into upper and lower surface patches where intrinsic shape analysis using conformal modules can be carried out. Based on these curves, we further propose a parameterization of HC surfaces called the eigen-harmonic parameterization (EHP). EHP maps each HC surface onto a parameter domain and imposes longitudinal and azimuthal coordinates on each surface, which follow the gradient and level sets of its first nontrivial Laplace--Beltrami eigenfunction, respectively. Each tubular domain is constructed according to the geometry of an individual HC surface. This gives a parameter domain with much less geometric distortion compared to spherical parameterization. With EHP, all HC surfaces are automatically registered with intrinsic feature curves preserved and geometric distortions minimized. This allows shape analysis on any number of HC surfaces to be performed consistently. We studied geometric changes over time in 138 HC surfaces of patients with AD and normal subjects scanned at two different times. We successfully located areas with significantly different shape changes over time between the two groups.

In this paper, we propose novel algorithms for inpainting and refinement of diffeomorphisms. We first represent a diffeomorphism by its Beltrami coefficient. Then it is possible to refine and inpaint the diffeomorphism by processing this Beltrami coefficient. With the inpainted/refined Beltrami coefficient, we construct a new diffeomorphism using the exact Beltrami holomorphic flow algorithm proposed in this paper. We apply our algorithms on several practical applications, which include the inpainting of a highly distorted diffeomorphism, the inpainting of image sequences of deforming shapes, the super-resolution of diffeomorphisms and the global parameterization of cortical surfaces by combining local parameterizations. Experiments show that our algorithm can solve these problems with natural and smooth results. We demonstrate how our proposed method can be widely applied in areas from texture mapping to video processing, and from computer graphics to medical imaging.

Surface registration is widely used in machine vision and medical imaging, where 1-1 correspondences between surfaces are computed to study their variations. Surface maps are usually stored as the 3D coordinates each vertex is mapped to, which often requires lots of storage memory. This causes inconvenience in data transmission and data storage, especially when a large set of surfaces are analyzed. To tackle this problem, we propose a novel representation of surface diffeomorphisms using Beltrami coefficients, which are complex-valued functions defined on surfaces with supreme norm less than 1. Fixing any 3 points on a pair of surfaces, there is a 1-1 correspondence between the set of surface diffeomorphisms between them and the set of Beltrami coefficients on the source domain. Hence, every bijective surface map can be represented by a unique Bel-trami coefficient. Conversely, given a Beltrami coefficient, we can reconstruct the unique surface map associated to it using the Beltrami Holomorphic flow (BHF) method introduced in this paper. Using this representation, 1/3 of the storage space is saved. We can further reduce the storage requirement by 90% by compressing the Beltrami coefficients using Fourier approximations. We test our algorithm on synthetic data, real human brain and hippocampal surfaces. Our results show high accuracy in the reconstructed data, while the amount of storage is greatly reduced. Our approach is compared with the Fourier compression of the coordinate functions using the same amount of data. The latter approach often shows jaggy results and cannot guarantee to preserve diffeomorphisms.