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.

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.

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.

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.

Finding surface mappings with least distortion arises from many applications in various fields. Extremal Teichmüller maps are surface mappings with least conformality distortion. The existence and uniqueness of the extremal Teichmüller map between Riemann surfaces of finite type are theoretically guaranteed (see Fletcher and Markovic, Quasiconformal maps and Teichmüller theory, Oxford Graduate Texts in Math., vol. 11, Oxford University Press, Oxford, 2007). Recently, a simple iterative algorithm for computing the Teichmüller maps between connected Riemann surfaces with given boundary value was proposed by Lui, Lam, Yau, and Gu in Teichmüller extremal mapping and its applications to landmark matching registration, arXiv:1211.2569. Numerical results were reported in the paper to show the effectiveness of the algorithm. The method was successfully applied to landmarkmatching registration. The purpose of this paper is to prove the iterative algorithm proposed in loc. cit., indeed converges.