Lok Ming Lui

Brain Surface Conformal Parameterization Using Riemann Surface Structure

In medical imaging, parameterized 3-D surface models are useful for anatomical modeling and visualization, statistical comparisons of anatomy, and surface-based registration and signal processing. Here we introduce a parameterization method based on Riemann surface structure, which uses a special curvilinear net structure (conformal net) to partition the surface into a set of patches that can each be conformally mapped to a parallelogram. The resulting surface subdivision and the parameterizations of the components are intrinsic and stable (their solutions tend to be smooth functions and the boundary conditions of the Dirichlet problem can be enforced). Conformal parameterization also helps transform partial differential equations (PDEs) that may be defined on 3-D brain surface manifolds to modified PDEs on a two-dimensional parameter domain. Since the Jacobian matrix of a conformal parameterization is diagonal, the modified PDE on the parameter domain is readily solved. To illustrate our techniques, we computed parameterizations for several types of anatomical surfaces in 3-D magnetic resonance imaging scans of the brain, including the cerebral cortex, hippocampi, and lateral ventricles. For surfaces that are topologically homeomorphic to each other and have similar geometrical structures, we show that the parameterization results are consistent and the subdivided surfaces can be matched to each other. Finally, we present an automatic sulcal landmark location algorithm by solving PDEs on cortical surfaces. The landmark detection results are used as constraints for building conformal maps between surfaces that also match explicitly defined landmarks.

Optimized Conformal Surface Registration with Shape-based Landmark Matching

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.

Teichmuller Mapping (T-Map) and Its Applications to Landmark Matching Registration

Registration, which aims to find an optimal 1-1 correspondence between shapes, is an important process in different research areas. Landmark-based surface registration has been widely studied to obtain a mapping between shapes that matches important features. Obtaining a unique and bijective surface registration that matches features consistently is generally challenging, especially when a large number of landmark constraints are enforced. This motivates us to search for a unique landmark matching surface diffeomorphism, which minimizes the local geometric distortion. For this purpose, we propose a special class of diffeomorphisms called the Teichmuller mappings (T-Maps). Under suitable conditions on the landmark constraints, a unique T-Map between two surfaces can be obtained, which minimizes the maximal conformality distortion. The conformality distortion measures how far the mapping deviates from a conformal mapping, and hence it measures the local geometric distortion. In this paper, we propose an eff...

Texture Map and Video Compression Using Beltrami Representation

Surface parameterizations and registrations are important in computer graphics and imaging, where 1-1 correspondences between meshes are computed. In practice, surface maps are usually represented and stored as three-dimensional coordinates each vertex is mapped to, which often requires lots of memory. This causes inconvenience in data transmission and data storage. To tackle this problem, we propose an effective algorithm for compressing surface homeomorphisms using Fourier approximation of the Beltrami representation. The Beltrami representation is a complex-valued function defined on triangular faces of the surface mesh with supreme norm strictly less than 1. Under suitable normalization, there is a 1-1 correspondence between the set of surface homeomorphisms and the set of Beltrami representations. Hence, every bijective surface map is associated with a unique Beltrami representation. Conversely, given a Beltrami representation, the corresponding bijective surface map can be exactly reconstructed usin...

FLASH: Fast Landmark Aligned Spherical Harmonic Parameterization for Genus-0 Closed Brain Surfaces

Surface registration between cortical surfaces is crucial in medical imaging for performing systematic comparisons between brains. Landmark-matching registration that matches anatomical features, called the sulcal landmarks, is often required to obtain a meaningful 1-1 correspondence between brain surfaces. This is commonly done by parameterizing the surface onto a simple parameter domain, such as the unit sphere, in which the sulcal landmarks are consistently aligned. Landmark-matching surface registration can then be obtained from the landmark aligned parameterizations. For genus-0 closed brain surfaces, the optimized spherical harmonic parameterization, which aligns landmarks to consistent locations on the sphere, has been widely used. This approach is limited by the loss of bijectivity under large deformations and the slow computation. In this paper, we propose FLASH, a fast algorithm to compute the optimized spherical harmonic parameterization with consistent landmark alignment. This is achieved by f...

Landmark- and Intensity-Based Registration with Large Deformations via Quasi-conformal Maps

Registration, which aims to find an optimal one-to-one correspondence between different data, is an important problem in various fields. This problem is especially challenging when large deformations occur. In this paper, we present a novel algorithm to obtain diffeomorphic image or surface registrations with large deformations via quasi-conformal maps. The basic idea is to minimize an energy functional involving a Beltrami coefficient term, which measures the distortion of the quasi-conformal map. The Beltrami coefficient effectively controls the bijectivity and smoothness of the registration, even with very large deformations. Using the proposed algorithm, landmark-based registration between images or surfaces can be effectively computed. The obtained registration is guaranteed to be diffeomorphic (1-1 and onto), even with a large deformation or large number of landmark constraints. The proposed algorithm can also be combined with matching intensity (such as image intensity or surface curvature) to improve the accuracy of the registration. Experiments have been carried out on both synthetic and real data. Results demonstrate the efficacy of the proposed algorithm to obtain diffeomorphic registration between images or surfaces.

Shape-based diffeomorphic registration on hippocampal surfaces using beltrami holomorphic flow

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 (Alzheimers 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.

Solving PDEs on manifolds with global conformal parametriazation

In this paper, we propose a method to solve PDEs on surfaces with arbitrary topologies by using the global conformal parametrization. The main idea of this method is to map the surface conformally to 2D rectangular areas and then transform the PDE on the 3D surface into a modified PDE on the 2D parameter domain. Consequently, we can solve the PDE on the parameter domain by using some well-known numerical schemes on ℝ2. To do this, we have to define a new set of differential operators on the manifold such that they are coordinates invariant. Since the Jacobian of the conformal mapping is simply a multiplication of the conformal factor, the modified PDE on the parameter domain will be very simple and easy to solve. In our experiments, we demonstrated our idea by solving the Navier-Stoke’s equation on the surface. We also applied our method to some image processing problems such as segmentation, image denoising and image inpainting on the surfaces.

Shape analysis of vestibular systems in adolescent idiopathic scoliosis using geodesic spectra

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.

Folding-Free Global Conformal Mapping for Genus-0 Surfaces by Harmonic Energy Minimization

Surface conformal maps between genus-0 surfaces play important roles in applied mathematics and engineering, with applications in medical image analysis and computer graphics. Previous work (Gu and Yau in Commun Inf Syst 2(2):121–146, 2002) introduces a variational approach, where global conformal parameterization of genus-0 surfaces was addressed through minimizing the harmonic energy, with two weaknesses: its gradient descent iteration is slow, and its solutions contain undesired parameterization foldings when the underlying surface has long sharp features. In this paper, we propose an algorithm that significantly accelerates the harmonic energy minimization and a method that iteratively removes foldings by taking advantages of the weighted Laplace–Beltrami eigen-projection. Experimental results show that the proposed approaches compute genus-0 surface harmonic maps much faster than the existing algorithm in Gu and Yau (Commun Inf Syst 2(2):121–146, 2002) and the new results contain no foldings.