Surface registration, which transforms different sets of surface data into one common reference space, is an important process which allows us to compare or integrate the surface data effectively. If a nonrigid transformation is required, surface registration is commonly done by parameterizing the surfaces onto a simple parameter domain, such as the unit square or sphere. In this work, we are interested in looking for meaningful registrations between surfaces through parameterizations, using prior features in the form of landmark curves on the surfaces. In particular, we generate optimized conformal parameterizations which match landmark curves exactly with shape-based correspondences between them. We propose a variational method to minimize a compound energy functional that measures the harmonic energy of the parameterization maps and the shape dissimilarity between mapped points on the landmark curves. The novelty is that the computed maps are guaranteed to align the landmark features consistently and give a shape-based diffeomorphism between the landmark curves. We achieve this by intrinsically modeling our search space of maps as flows of smooth vector fields that do not flow across the landmark curves. By using the local surface geometry on the curves to define a shape measure, we compute registrations that ensure consistent correspondences between anatomical features. We test our algorithm on synthetic surface data. An application of our model to medical imaging research is shown, using experiments on brain cortical surfaces, with anatomical (sulcal) landmarks delineated, which show that our computed maps give a shape-based alignment of the sulcal curves without significantly impairing conformality. This ensures correct averaging and comparison of data across subjects.

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.

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 Beltrami 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.

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.

We develop a new algorithm to automatically register hippocampal (HP) surfaces with complete geometric matching, avoiding the need to manually label landmark features. A good registration depends on a reasonable choice of shape energy that measures the dissimilarity between surfaces. In our work, we first propose a complete shape index using the Beltrami coefficient and curvatures, which measures subtle local differences. The proposed shape energy is zero if and only if two shapes are identical up to a rigid motion. We then seek the best surface registration by minimizing the shape energy. We propose a simple representation of surface diffeomorphisms using Beltrami coefficients, which simplifies the optimization process. We then iteratively minimize the shape energy using the proposed Beltrami Holomorphic flow (BHF) method. Experimental results on 212 HP of normal and diseased (Alzheimer's disease) subjects show our proposed algorithm is effective in registering HP surfaces with complete geometric matching. The proposed shape energy can also capture local shape differences between HP for disease analysis.

The study of 2D shapes is a central problem in the field of computer vision. In 2D shape analysis, classification and recognition of objects from their observed silhouettes are extremely crucial and yet difficult. It usually involves an efficient representation of 2D shape space with natural metric, so that its mathematical structure can be used for further analysis. Although significant progress has been made for the study of 2D simply-connected shapes, very few works have been done on the study of 2D objects with arbitrary topologies. In this work, we propose a representation of general 2D domains with arbitrary topologies using conformal geometry. A natural metric can be defined on the proposed representation space, which gives a metric to measure dissimilarities between objects. The main idea is to map the exterior and interior of the domain conformally to unit disks and circle domains, using holomorphic 1-forms. A set of diffeomorphisms from the unit circle S^1 to itself can be obtained, which together with the conformal modules are used to define the shape signature. We prove mathematically that our proposed signature uniquely represents shapes with arbitrary topologies. We also introduce a reconstruction algorithm to obtain shapes from their signatures. This completes our framework and allows us to move back and forth between shapes and signatures. Experiments show the efficacy of our proposed algorithm as a stable shape representation scheme.

There is an unmet medical need to identify neuroimaging biomarkers that is able to accurately diagnose and monitor Alzheimer's disease (AD) at very early stages and assess the response to AD-modifying therapies. To a certain extent, volumetric and functional magnetic resonance imaging (fMRI) studies can detect changes in structure, cerebral blood flow and blood oxygenation that are able to distinguish AD and mild cognitive impairment (MCI) subjects from normal controls. However, it has been challenging to use fully automated MRI analytic methods to identify potential AD neuroimaging biomarkers. We have thus proposed a method based on independent component analysis (ICA), for studying potential AD-related MR image features, coupled with the use of support vector machine (SVM) for classifying scans into categories of AD, MCI, and normal control (NC) subjects. The MRI data were selected from Open Access Series of Imaging Studies (OASIS) and the Alzheimer's Disease Neuroimaging Initiative (ADNI) databases. The experimental results showed that our ICA-based method can differentiate AD and MCI subjects from normal controls, although further methodological improvement in the analytic method and inclusion of additional variables may be required for optimal classification.

This paper proposes a novel algorithm to extract feature landmarks on the vestibular system (VS), for the analysis of Adolescent Idiopathic Scoliosis (AIS) disease. AIS is a 3-D spinal deformity commonly occurred in adolescent girls with unclear etiology. One popular hypothesis was suggested to be the structural changes in the VS that induce the disturbed balance perception, and further cause the spinal deformity. The morphometry of VS to study the geometric differences between the healthy and AIS groups is of utmost importance. However, the VS is a genus-3 structure situated in the inner ear. The high-genus topology of the surface poses great challenge for shape analysis. In this work, we present a new method to compute exact geodesic loops on the VS. The resultant geodesic loops are in Euclidean metric, thus characterizing the intrinsic geometric properties of the VS based on the real background geometry. This leads to more accurate results than existing methods, such as the hyperbolic Ricci flow method. Furthermore, our method is fully automatic and highly efficient, e.g., one order of magnitude faster than. We applied our algorithm to the VS of normal and AIS subjects. The promising experimental results demonstrate the efficacy of our method and reveal more statistically significant shape difference in the VS between right-thoracic AIS and normal subjects.

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.

In shape analysis, finding an optimal 1-1 correspondence between 3D surfaces within a large class of admissible bijective mappings is of great importance. Such a process is called surface registration. The difficulty lies in the fact that the space of all surface diffeomorphisms is a complicated functional space, making it challenging to exhaustively search for the best mapping. To tackle this problem, we propose a simple representation of bijective surface maps using Beltrami coefficients (BCs)—complex-valued functions defined on surfaces with supremum 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 BCs. Hence, every bijective surface map may be represented by a unique BC. Conversely, given a BC, we can reconstruct the unique surface map associated with it using the Beltrami Holomorphic flow (BHF) method. Using BCs to represent surface maps is advantageous because it is a much simpler functional space, which captures many essential features of a surface map. By adjusting BCs, we equivalently adjust surface diffeomorphisms to obtain the optimal map with desired properties. More specifically, BHF gives us the variation of the associated map under the variation of BC. Using this, a variational problem over the space of surface diffeomorphisms can be easily reformulated into a variational problem over the space of BCs. This makes the minimization procedure much easier. More importantly, the diffeomorphic property is always preserved. We test our method on synthetic examples and real medical applications. Experimental results demonstrate the effectiveness of our proposed algorithm for surface registration.