Min ZhangStony Brook UniversityRen GuoOregon State UniveristyWei ZengSchool of Computing and Information Sciences, Florida International UniversityFeng LuoRutgers UniversityShing Tung YauHarvard UniversityXianfeng GuStony Brook Univerisity
Computational GeometryDifferential GeometryGeometric Modeling and ProcessingConvex and Discrete Geometry mathscidoc:1612.01001
Graphical Models/Geometric Modeling and Processing 2014, 76, (5), 321-339, 2014.9
Ricci ﬂow deformsthe Riemannian metric proportionallyto the curvature, such that the curvatureevolves accordingto a heat diffusion process and eventually becomes constant everywhere. Ricci ﬂow has demonstrated its great potential by solving various problems in many ﬁelds, which can be hardly handled by alternative methods so far. This work introduces the uniﬁed theoretic framework for discrete Surface Ricci Flow, including all the common schemes: Tangential Circle Packing, Thurston’s Circle Packing, Inversive Distance Circle Packing and Discrete Yamabe Flow. Furthermore, this work also introduces a novel schemes, Virtual Radius Circle Packing and the Mixed Type schemes, under the uniﬁed framework. This work gives explicit geometric interpretation to the discrete Ricci energies for all the schemes with all back ground geometries, and the corresponding Hessian matrices. The uniﬁed frame work deepens our understanding to the the discrete surface Ricci ﬂow theory, and has inspired us to discover the new schemes, improved the ﬂexibility and robustness of the algorithms, greatly simpliﬁed the implementation and improved the efﬁciency. Experimental results show the uniﬁed surface Ricci ﬂow algorithms can handle general surfaces with different topologies, and is robust to meshes with different qualities, and is effective for solving real problems.
Morphing is the process of changing a geometric model or an image into another. The process generally involves rigid body motions and non-rigid deformations. It is well known that there exists a unique conformal mapping from a simply connected surface into a unit disk by the Riemann mapping theorem. On the other hand, a 3D surface deformable model can be built via various approaches such as mutual parameterization from direct interpolation or surface matching using landmarks. In this paper, a numerical methods of 3D surface morphing based on deformable model and conformal mapping is demonstrated.
We take the advantage of the unique representation of 3D surfaces by the mean curvatures $H$ and the conformal factors $\lambda$ associated with the Riemann mapping and build up the deformation model by consistently registering the landmarks on the conformal parametric domains. As a result, the correspondence of the $(H, \lambda)$ between two surfaces can be defined and a 3D deformation field can be reconstructed. Furthermore, by composition of the M\"obius transformation and the 3D deformation field, a smooth morphing sequence can be generated over a consistent mesh structure via the cubic spline homotopy. Several numerical experiments on the face morphing are presented to demonstrate the robustness of our approach.
Ronald Lok Ming LuiHarvard UniversityTsz Wai WongUniversity of California, Los AngelesWei ZengStony Brook UniversityXianfeng GuStony Brook UniversityPaul M. ThompsonUniversity of California, Los AngelesTony F. ChanThe Hong Kong University of Science and TechnologyShing-Tung YauHarvard University
Journal of Inverse Problem and Imaging, 4, (2), 311 - 333, 2010.5
We address the problem of detecting deformities on elastic surfaces. This is of great importance for shape analysis, with applications such as detecting abnormalities in biological shapes (e.g., brain structures). We propose an effective algorithm to detect abnormal deformations by generating quasi-conformal maps between the original and deformed surfaces. We firstly flatten the 3D surfaces conformally onto 2D rectangles using the discrete Yamabe flow and use them to compute a quasi-conformal map that matches internal features lying within the surfaces. The deformities on the elastic surface are formulated as non-conformal deformations, whereas normal deformations that preserve local geometry are formulated as conformal deformations. We then detect abnormalities by computing the Beltrami coefficient associated uniquely with the quasi-conformal map. The Beltrami coefficient is a complex-valued function defined on the surface. It describes the deviation of the deformation from conformality at each point. By considering the norm of the Beltrami coefficient, we can effectively segment the regions of abnormal changes, which are invariant under normal (non-rigid) deformations that preserve local geometry. Furthermore, by considering the argument of the Beltrami coefficient, we can capture abnormalities induced by local rotational changes. We tested the algorithm by detecting abnormalities on synthetic surfaces, 3D human face data and MRI-derived brain surfaces. Experimental results show that our algorithm can effectively detect abnormalities and capture local rotational alterations. Our method is also more effective than other existing methods, such as the isometric indicator, for locating abnormalities.
Chengfeng WenStony Brook UniversityDefeng WangThe Chinese University of Hong KongLin ShiThe Chinese University of Hong KongWinnie C.W. ChuThe Chinese University of Hong KongJack C.Y. ChengThe Chinese University of Hong KongRonald Lok Ming LuiThe Chinese University of Hong Kong
Computerized Medical Imaging and Graphics, 44, 1-12, 2015.9
The analysis of the vestibular system (VS) is an important research topic in medical image analysis. VS is a sensory structure in the inner ear for the perception of spatial orientation. It is believed several diseases, such as the Adolescent Idiopathic Scoliosis (AIS), are due to the impairment of the VS function. The morphology of the VS is thus of great research significance. A major challenge is that the VS is a genus-3 surface. The high-genus topology of the VS poses great challenges to find accurate pointwise correspondences between the surfaces and whereby perform accurate shape analysis. In this paper, we present a method to obtain the landmark constrained diffeomorphic registration between the VS surfaces based on the quasi-conformal theory. Given a set of corresponding landmarks on the VS surfaces, a diffeomorphism between the VS surfaces that matches the features consistently can be obtained. The basic idea is to iteratively search for an admissible Beltrami coefficient, which is associated to our desired landmark matching registration. With the obtained surface registrations, vertex-wise morphometric analysis can be carried out. Two types of geometric features are used for shape comparison. One is the collection of homotopic loops on each canals of the VS, which can be used to measure the local thickness of the canals. From the homotopic loops, centerlines can be extracted. By examining the deviations of the centerlines from the best fit planes, bendings of the canals can be detected. The second geometric feature is the minimal surface enclosed by the homotopic loop. From the minimal surfaces of each homotopic loops, cross-sectional area of the canals can be evaluated. To study the local shape difference more comprehensively, a complete shape index, which is defined using the Beltrami coefficients and surface curvatures, is used. We test proposed registration method on 15 VS of normal control subjects and 12 VS of patients suffering from AIS. Experimental results show the efficacy and accuracy of the proposed algorithm to compute the VS surface registration. Shape analysis has also been carried out using the proposed geometric features and shape index, which reveals shape differences in the posterior canal between normal and diseased AIS groups.
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.