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. Landmarkmatching 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 formulating the optimization problem to C and thereby linearizing the problem. Errors introduced near the pole are corrected using quasiconformal theories. Also, by adjusting the Beltrami differential of the mapping, a diffeomorphic (1-1, onto) spherical parameterization can be effectively obtained. The proposed algorithm has been tested on 38 human brain surfaces. Experimental results demonstrate that the computation of the landmark aligned spherical harmonic parameterization is significantly accelerated using the proposed algorithm.

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

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 using the linear Beltrami solver introduced in this paper. Using the Beltrami representation, the surface homeomorphism can be easily compressed by Fourier approximation, without distorting the bijectivity of the map. The storage requirement can be effectively reduced, which is useful for many practical problems in computer graphics and imaging. In this paper, we propose applying the algorithm to texture map compression and video compression. With our proposed algorithm, the storage requirement for the texture properties of a textured surface can be significantly reduced. Our algorithm can further be applied to compressing motion vector fields for video compression, which effectively improves the compression ratio.

We present a new approach to obtain diffeomorphic registrations with large deformations using landmark and intensity information via quasi-conformal maps. The basic idea is to minimize an energy functional involving a Beltrami coefficient term, which measures the distortion of the quasiconformal map. The Beltrami coefficient effectively controls the bijectivity and smoothness of the registration. In this paper, we first propose the quasi-conformal landmark registration (QCLR) algorithm to obtain diffeomorphic (1-1 and onto) registrations between images or surfaces. Using QCLR, landmark-aligned diffeomorphisms between images or surfaces can be obtained, even with a large geometric difference or a large number of landmark constraints. This algorithm is then extended to the quasi-conformal hybrid registration (QCHR) algorithm, which combines landmark and intensity (such as image intensity or surface curvature) information to achieve a more accurate registration result. Experiments have been carried out on both synthetic and real data. Results demonstrate the stability and efficacy of the proposed algorithm to obtain diffeomorphic registrations between images or surfaces.

No abstract

Abstract.To compare and integrate brain data, data from multiple subjects are typically mapped into a canonical space. One method to do this is to conformally map cortical surfaces to the sphere. It is well known that any genus zero Riemann surface can be conformally mapped to a sphere. Therefore, conformal mapping offers a convenient method to parameterize cortical surfaces without angular distortion, generating an orthogonal grid on the cortex that locally preserves the metric. To compare cortical surfaces more effectively, it is advantageous to adjust the conformal parameterizations to match consistent anatomical features across subjects. This matching of cortical patterns improves the alignment of data across subjects, although it is more challenging to create a consistent conformal (orthogonal) parameterization of anatomy across subjects when landmarks are constrained to lie at specific locations in the spherical parameter space. Here we propose a new method, based on a new energy functional, to optimize the conformal parameterization of cortical surfaces by using landmarks. Experimental results on a dataset of 40 brain hemispheres showed that the landmark mismatch energy can be greatly reduced while effectively preserving conformality. The key advantage of this conformal parameterization approach is that any local adjustments of the mapping to match landmarks do not affect the conformality of the mapping significantly. We also examined how the parameterization changes with different weighting factors. As expected, the landmark matching error can be reduced if it is more heavily penalized, but conformality is progressively reduced.

Surface parameterizations have been widely used in computer graphics and geometry processing. In particular, as simply-connected open surfaces are conformally equivalent to the unit disk, it is desirable to compute the disk conformal parameterizations of the surfaces. In this paper, we propose a novel algorithm for the conformal parameterization of a simply-connected open surface onto the unit disk, which significantly speeds up the computation, enhances the conformality and stability, and guarantees the bijectivity. The conformality distortions at the inner region and on the boundary are corrected by two steps, with the aid of an iterative scheme using quasi-conformal theories. Experimental results demonstrate the effectiveness of our proposed method.

In shape analysis, finding an optimal 1-1 correspondence between surfaces within a large class of admissible bijective mappings is of great importance. Such process is called surface registration. The difficulty lies in the fact that the space of all surface diffeomorphisms is a complicated functional space, making exhaustive search for the best mapping challenging. To tackle this problem, we propose a simple representation of bijective surface maps using Beltrami coefficients (BCs), 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 BCs. Hence, every bijective surface map can be represented by a unique BC. Conversely, given a BC, we can reconstruct the unique surface map associated to 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.