The medial axis is an important shape representation that finds a wide range of applications in shape analysis. For large-scale shapes of high resolution, a progressive medial axis representation that starts with the lowest resolution and gradually adds more details is desired. In this paper, we propose a fast and robust geometric algorithm that computes progressive medial axes of a large-scale planar shape. The key ingredient of our method is a novel structural analysis of merging medial axes of two planar shapes along a shared boundary. Our method is robust by separating the analysis of topological structure from numerical computation. Our method is also fast and we show that the time complexity of merging two medial axes is O(n log nv), where n is the number of total boundary generators, nv is strictly smaller than n and behaves as a small constant in all our experiments. Experiments on large-scale polygonal data and comparison with state-of-the-art methods show the efficiency and effectiveness of the proposed method.

Here we propose a novel method to compute Teichmller shape space based shape index to study brain morphometry. Such a shape index is intrinsic, and invariant under conformal transformations, rigid motions and scaling. We conformally map a genus-zero open boundary surface to the Poincar disk with the Yamabe flow method. The shape indices that we compute are the lengths of a special set of geodesics under hyperbolic metric. Tests on longitudinal brain imaging data were used to demonstrate the stability of the derived feature vectors. In leave-one-out validation tests, we achieved 100% accurate classification (versus only 68% accuracy for volume measures) in distinguishing 11 HIV/AIDS individuals from 8 healthy control subjects, based on Teichmller coordinates for lateral ventricular surfaces extracted from their 3D MRI scans.

Fusion of images with same or different modalities has been conquering medical imaging field more rapidly due to the presence of highly accessible patients’ information in recent years. For example, cross platform non-rigid registration of CT with MRI images has found a significant role in different clinical application. In some instances labelling of anatomical features by medical experts are also involved to further improve the accuracy and authenticity of the registration. Being motivated by these, we propose a new algorithm to compute diffeomorphic hybrid multi-modality registration with large deformations. Our iterative scheme consists of mainly two steps. First, we obtain the optimal Beltrami coefficient corresponding to the diffeomorphic mapping that exactly superimposes the feature points. The second step detects the intensity difference in the framework of mutual information. A non-rigid deformation which minimizes the intensity difference is then obtained. Experiments have been carried out on both synthetic and real data. Results demonstrate the stability and efficacy of the proposed algorithm to obtain diffeomorphic image registration.

We propose a novel method to apply Teichmller space theory to study the signature of a family of nonintersecting closed 3D curves on a general genus zero closed surface. Our algorithm provides an efficient method to encode both global surface and local contour shape information. The signatureTeichmller shape descriptoris computed by surface Ricci flow method, which is equivalent to solving an elliptic partial differential equation on surfaces and is numerically stable. We propose to apply the new signature to analyze abnormalities in brain cortical morphometry. Experimental results with 3D MRI data from Alzheimers disease neuroimaging initiative (ADNI) dataset [152 healthy control subjects versus 169 Alzheimers disease (AD) patients] demonstrate the effectiveness of our method and illustrate its potential as a novel surface-based cortical morphometry measurement in AD research.

This paper presents a method to compute the quasi-conformal parameterization (QCMC) for a multiply-connected 2D domain or surface. QCMC computes a quasi-conformal map from a multiply-connected domain S onto a punctured disk D_S associated with a given Beltrami differential. The Beltrami differential, which measures the conformality distortion, is a complex-valued function μ: S→ℂ with supremum norm strictly less than 1. Every Beltrami differential gives a conformal structure of S. Hence, the conformal module of D_S, which are the radii and centers of the inner circles, can be fully determined by μ, up to a Möbius transformation. In this paper, we propose an iterative algorithm to simultaneously search for the conformal module and the optimal quasi-conformal parameterization. The key idea is to minimize the Beltrami energy subject to the boundary constraints. The optimal solution is our desired quasi-conformal parameterization onto a punctured disk. The parameterization of the multiply-connected domain simplifies numerical computations and has important applications in various fields, such as in computer graphics and vision. Experiments have been carried out on synthetic data together with real multiply-connected Riemann surfaces. Results show that our proposed method can efficiently compute quasi-conformal parameterizations of multiply-connected domains and outperforms other state-of-the-art algorithms. Applications of the proposed parameterization technique have also been explored.