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.

Adolescent Idiopathic Scoliosis (AIS) characterized by the 3D spine deformity affects about 4% schoolchildren worldwide. One of the prominent theories of the etiopathogenesis of AIS was proposed to be the poor postural balance control due to the impaired vestibular function. Thus, the morphometry of the vestibular system (VS) is of great importance for studying AIS. The VS is a genus-3 structure situated in the inner ear and consists of three semicircular canals lying perpendicular to each other. The high-genus topology of the surface poses great challenge for shape analysis. In this work, we propose an effective method to analyze shapes of high-genus surfaces by considering their geodesic spectra. The key is to compute the canonical hyperbolic geodesic loops of the surface, using the Ricci flow method. The Fuchsian group generators are then computed which can be used to determine the geodesic spectra. The geodesic spectra effectively measure shape differences between high-genus surfaces up to the hyperbolic isometry. We applied the proposed algorithm to the VS of 12 normal and 15 AIS subjects. Experimental results show the effectiveness of our algorithm and reveal statistical shape difference in the VS between right-thoracic AIS and normal subjects.

In this paper, we consider the problem of approximating the longest path in a grid graph that possesses a square-free 2-factor. By (1) analyzing several characteristics of four types of cross cells and (2) presenting three cycle-merging operations and one extension operation, we propose an algorithm that finds in an n-vertex grid graph G, a long path of length at least 5/6 n + 2. Our approximation algorithm runs in quadratic time. In addition to its theoretical value, our work has also an interesting application in image-guided maze generation.

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.

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.