Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2970526

Formats

Article sections

- Abstract
- 1 Introduction
- 2 A Variational Framework
- 3 Laplace-Beltrami Eigen-Features
- 4 Experimental Results
- 5 Conclusions and Future Work
- References

Authors

Related links

Inf Process Med Imaging. Author manuscript; available in PMC 2010 November 2.

Published in final edited form as:

Inf Process Med Imaging. 2009; 21: 467–478.

PMCID: PMC2970526

NIHMSID: NIHMS242513

See other articles in PMC that cite the published article.

We propose in this work a novel variational method for computing maps between surfaces by combining informative geometric features and regularizing forces including inverse consistency and harmonic energy. To tackle the ambiguity in defining homologous points on smooth surfaces, we design feature functions in the data term based on the Reeb graph of the Laplace-Beltrami eigenfunctions to quantitatively describe the global geometry of elongated anatomical structures. For inverse consistency and robustness, our method computes simultaneously the forward and backward map by iteratively solving partial differential equations (PDEs) on the surfaces. In our experiments, we successfully mapped 890 hippocampal surfaces and report statistically significant maps of atrophy rates between normal controls and patients with mild cognitive impairment (MCI) and Alzheimer’s disease (AD).

Surface mapping is an important problem in medical image analysis with applications in population studies[1] and the creation of 3D shape prior models[2]. While significant progresses have been made with the development of various techniques[3–9], most of them focus heavily on the *regularization* part of the problem. The lack of informative *data terms* makes it still a challenging problem to resolve the ambiguity in defining homologous points on smooth surfaces. To this end we propose in this paper a novel variational framework to compute maps directly between surfaces. By combining data terms from Laplace-Beltrami eigen-features and regularizing forces from inverse consistency and harmonic energy, our method can generate robust and high quality maps for a class of subcortical structures with clinical significance.

Many interesting techniques have been proposed in previous work to solve various surface mapping problems. A popular approach is to first map the surface onto the sphere and then solve the registration problem in this canonical domain [5–7, 9]. By viewing the surface as a subset of ^{3}, successful image registration techniques were applied to compute surface maps with possible landmark constraints [3, 4, 8]. Landmark curves detected by shape context features were used to guide the mapping of hippocampal surfaces[10]. The medial model provides a compact surface representation and is also very useful to construct maps between surfaces[11–13].

We propose in this work a variational approach to automatically compute both the forward and backward maps between two surfaces with inverse consistency. An important feature of our method is that we use the Reeb graph of the Laplace-Beltrami eigenfunctions to define intrinsic features for the data fidelity term in our energy function. These features provide quantitative descriptions of the salient geometry shared by elongated structures such as the hippocampus and significantly remove the ambiguity in finding correspondences for such surfaces. Another interesting feature of our method is that we incorporate both inverse consistency[14] and harmonic energy [15] as the regularization term to improve the quality of maps. The harmonic energy encourages the smoothness in the maps and the inverse consistency term helps reduce area distortions. In our experiments, we perform extensive tests to demonstrate the quality of the maps generated by our method. We also demonstrate the robustness of our method with its successful mapping of 890 hippocampal surfaces and report statistically significant results.

The rest of the paper is organized as follows. In section 2, we develop the variational framework that combines data terms of intrinsic features and regularization terms including both inverse consistency and harmonic energy. We propose the Laplace-Beltrami eigen-features for the modeling of elongated structures such as hippocampus in section 3. Experimental results are presented in section 4 to demonstrate the quality and robustness of our mapping algorithm. Finally conclusions and future work are discussed in section 5.

In this section, we propose a variational framework for the direct mapping of surfaces. Let (, *g*) and (, *h*) be two Riemann surfaces with their metric tensor *g* and *h*. To obtain inverse consistency, we compute two maps jointly in our method. We denote *u*_{1}: → as the map from to , and *u*_{2}: → the map from to . We also assume there are *L* feature functions defined on each surface and denote
${\xi}_{1}^{j}:{\mathcal{M}}_{1}\to \mathbb{R}$ and
${\xi}_{2}^{j}:{\mathcal{M}}_{2}\to \mathbb{R}(1\le j\le L)$ as the *j*-th feature function on and , respectively.

In our variational framework, the maps are computed as the minimizer of the following energy function:

$$E={E}_{D}+{E}_{IC}+{E}_{H}.$$

(1)

There are three terms in the energy function and all of them are symmetric with respect to and . We call the first energy *E _{D}* as the data fidelity term and it is defined as:

$${E}_{D}=\sum _{j=1}^{L}{\alpha}_{D}^{j}\left[{\int}_{{\mathcal{M}}_{1}}{({\xi}_{1}^{j}-{\xi}_{2}^{j}\circ {u}_{1})}^{2}d{\mathcal{M}}_{1}+{\int}_{{\mathcal{M}}_{2}}{({\xi}_{2}^{j}-{\xi}_{1}^{j}\circ {u}_{2})}^{2}d{\mathcal{M}}_{2}\right].$$

(2)

This energy penalizes the mismatch of feature functions induced by the two maps, and the parameters ${\alpha}_{D}^{j}$ are used to assign proper weights for different feature functions.

The other two terms are for regularization, where *E _{IC}* encourages inverse consistency and

$${E}_{IC}={\alpha}_{IC}\left[{\int}_{{\mathcal{M}}_{1}}{(I-{u}_{2}\circ {u}_{1})}^{2}d{\mathcal{M}}_{1}+{\int}_{{\mathcal{M}}_{2}}{(I-{u}_{1}\circ {u}_{2})}^{2}d{\mathcal{M}}_{2}\right]$$

(3)

where *α _{IC}* is a regularization parameter. For inverse consistency, this energy penalizes the difference between the identity map

$${E}_{H}={\alpha}_{H}\left[{\int}_{{\mathcal{M}}_{1}}{\left|\right|{J}_{{u}_{1}}\left|\right|}^{2}d{\mathcal{M}}_{1}+{\int}_{{\mathcal{M}}_{2}}{\left|\right|{J}_{{u}_{2}}\left|\right|}^{2}d{\mathcal{M}}_{2}\right]$$

(4)

where *J*_{u1} and *J*_{u2} are the Jacobian of the two maps defined on the surfaces, and *α _{H}* is the regularization parameter for this term.

To minimize the energy function with respect to the maps *u*_{1} and *u*_{2}, we derive their gradient flows and compute them iteratively via the solution of these partial differential equations(PDEs) on the two surfaces and . For the energy *E _{D}*, the gradient flows of

$$\begin{array}{l}\frac{\partial {E}_{D}}{\partial {u}_{1}}=-2\sum _{j=1}^{L}{\alpha}_{D}^{j}({\xi}_{1}^{j}-{\xi}_{2}^{j}\circ {u}_{1}){\nabla}_{{\mathcal{M}}_{2}}{\xi}_{2}^{j}({u}_{1})\\ \frac{\partial {E}_{D}}{\partial {u}_{2}}=-2\sum _{j=1}^{L}{\alpha}_{D}^{j}({\xi}_{2}^{j}-{\xi}_{2}^{j}\circ {u}_{2}){\nabla}_{{\mathcal{M}}_{1}}{\xi}_{1}^{j}({u}_{2})\end{array}$$

(5)

where and denote the intrinsic gradient on the surfaces. For the energy *E _{IC}*, the gradient flows are:

$$\begin{array}{l}\frac{{E}_{IC}}{\partial {u}_{1}}=-2{\alpha}_{IC}\left[{J}_{{u}_{2}}^{T}(I-{u}_{2}\circ {u}_{1})+({u}_{2}^{-1}-{u}_{1})\mid {J}_{{u}_{2}^{-1}}\mid \right]\\ \frac{{E}_{IC}}{\partial {u}_{2}}=-2{\alpha}_{IC}\left[{J}_{{u}_{1}}^{T}(I-{u}_{1}\circ {u}_{2})+({u}_{1}^{-1}-{u}_{2})\mid {J}_{{u}_{1}^{-1}}\mid \right]\end{array}$$

(6)

where
${u}_{1}^{-1}$ and
${u}_{2}^{-1}$ are the inverse maps of *u*_{1} and *u*_{2}, and the determinant of their Jacobian are denoted as
$\mid {J}_{{u}_{2}^{-1}}\mid $ and
$\mid {J}_{{u}_{2}^{-1}}\mid $.

To derive the gradients of the harmonic energy, we denote (*x*^{1}, *x*^{2}) and (*y*^{1}, *y*^{2}) as the local coordinates of and . In the local coordinates, the maps can be represented as
${u}_{1}=({u}_{1}^{1},{u}_{1}^{2})$ and
${u}_{2}=({u}_{2}^{1},{u}_{2}^{2})$. The gradient flow of *u*_{1} can then be expressed in the form of Einstein’s summation as [15]:

$$\frac{\partial {E}_{H}}{\partial {u}_{1}^{r}}=-{\alpha}_{H}\phantom{\rule{0.16667em}{0ex}}\left[{\mathrm{\Delta}}_{{\mathcal{M}}_{1}}{u}_{1}^{r}+{g}^{\alpha \beta}{\mathrm{\Gamma}}_{2,pq}^{r}\frac{\partial {u}_{1}^{p}}{\partial {x}^{\alpha}}\frac{\partial {u}_{1}^{q}}{\partial {x}^{\beta}}\right]\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}(r=1,2)$$

(7)

where is the Laplace-Beltrami operator on , *g ^{αβ}* = (

$$\frac{\partial {E}_{H}}{\partial {u}_{2}^{s}}=-{\alpha}_{H}\left[{\mathrm{\Delta}}_{{\mathcal{M}}_{2}}{u}_{2}^{s}+{h}^{pq}{\mathrm{\Gamma}}_{1,\alpha \beta}^{s}\frac{\partial {u}_{2}^{\alpha}}{\partial {y}^{p}}\frac{\partial {u}_{2}^{\beta}}{\partial {y}^{q}}\right]\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}(s=1,2)$$

(8)

where is the Laplace-Beltrami operator on , *h ^{pq}* = (

By combining the above results, we have the gradient descent flows of *u*_{1} and *u*_{2} to minimize the energy:

$$\{\begin{array}{l}\frac{\partial {u}_{1}}{\partial t}=-\frac{\partial {E}_{D}}{\partial {u}_{1}}-\frac{\partial {E}_{IC}}{\partial {u}_{1}}-\frac{\partial {E}_{H}}{\partial {u}_{1}}\hfill \\ \frac{\partial {u}_{2}}{\partial t}=-\frac{\partial {E}_{D}}{\partial {u}_{2}}-\frac{\partial {E}_{IC}}{\partial {u}_{2}}-\frac{\partial {E}_{H}}{\partial {u}_{2}}\hfill \end{array}$$

(9)

To numerically solve these two equations on and , we use the approach of solving PDEs on implicit surfaces[17, 18, 16] and represent and as a signed distance function and *ψ*, respectively. With the implicit representation, we have a very simple formulation for the gradient flow of the harmonic energy. Using the signed distance function, we can express gradient operators on surfaces in terms of conventional gradient operators in ^{3}. For example, the intrinsic gradient *f* of a function *f* on can be expressed as *Π*_{}_{}*f*, where *f* is the gradient of *f* in ^{3} and *Π*_{}* _{}* =

$$\{\begin{array}{l}\frac{\partial {u}_{1}}{\partial t}=2\sum _{j=1}^{L}{\alpha}_{D}^{j}({\xi}_{1}^{j}-{\xi}_{2}^{j}\circ {u}_{1}){\mathrm{\Pi}}_{\nabla \psi}\nabla {\xi}_{2}^{j}({u}_{1})+{\alpha}_{IC}{\mathrm{\Pi}}_{\nabla \psi}\left[{J}_{{u}_{2}}^{T}(I-{u}_{2}\circ {u}_{1})+({u}_{2}^{-1}-{u}_{1})\mid {\mathrm{\Pi}}_{\nabla \phi}{J}_{{u}_{2}^{-1}}\mid \right]+{\alpha}_{H}{\mathrm{\Pi}}_{\nabla \psi}\nabla \xb7({\mathrm{\Pi}}_{\nabla \phi}{J}_{{u}_{1}}^{T})\hfill \\ \frac{\partial {u}_{2}}{\partial t}=2\sum _{j=1}^{L}{\alpha}_{D}^{j}({\xi}_{2}^{j}-{\xi}_{1}^{j}\circ {u}_{2}){\mathrm{\Pi}}_{\nabla \phi}\nabla {\xi}_{1}^{j}({u}_{2})+{\alpha}_{IC}{\mathrm{\Pi}}_{\nabla \phi}\left[{J}_{{u}_{1}}^{T}(I-{u}_{1}\circ {u}_{2})+({u}_{1}^{-1}-{u}_{2})\mid {\mathrm{\Pi}}_{\nabla \psi}{J}_{{u}_{1}^{-1}}\mid \right]+{\alpha}_{H}{\mathrm{\Pi}}_{\nabla \phi}\nabla \xb7({\mathrm{\Pi}}_{\nabla \psi}{J}_{{u}_{2}}^{T})\hfill \end{array}$$

(10)

For numerical efficiency, all computations are performed on narrow bands surrounding the surfaces. More details of the computational schemes can be found in [16].

To use the above variational framework, it is important to design proper features that can capture the common geometry across surfaces. This is generally a difficult problem and the solution is application dependent. In this section, we propose two feature functions to characterize the global geometry of elongated structures with neuroanatomical significance such as hippocampus and putamen. Both features are invariant to scale differences and natural pose variations. For a surface , the first feature function *ξ*_{1}: → characterizes the tail-to-head trend of elongated structures, and the second feature *ξ*_{2}: → describes the lateral profile of the surface. We next develop the algorithm to compute both features from the Laplace-Beltrami eigenfunction of .

For a surface , the eigenfunctions of its Laplace-Beltrami operator is de-fined as [19, 20]:

$${\mathrm{\Delta}}_{\mathcal{M}}f=-\lambda f$$

(11)

The spectrum of is discrete and the eigenvalues can be ordered as 0 = *λ*_{0} ≤ *λ*_{1} ≤ *λ*_{2} ≤ ···. For *λ _{i}*, the corresponding eigenfunction of

$${f}_{1}=arg\underset{f\perp {f}_{0},\left|\right|f\left|\right|\phantom{\rule{0.16667em}{0ex}}=1}{min}{\int}_{\mathcal{M}}{\mid {\nabla}_{\mathcal{M}}f\mid}^{2}d\mathcal{M}.$$

(12)

with

$${\lambda}_{1}={\int}_{\mathcal{M}}{\mid {\nabla}_{\mathcal{M}}f\mid}^{2}d\mathcal{M}.$$

(13)

Thus this function *f*_{1} can be viewed as the smoothest, non-constant map from to in the space orthogonal to *f*_{0}. To numerically compute the eigenfunction, we represent as a triangular mesh = (, ), where is the set of vertices and is the set of triangles. By using the weak form of (11) and the finite element method, we can compute the eigenfunction by solving a generalized matrix eigenvalue problem:

$$Qf=\lambda Uf$$

(14)

where *Q* and *U* are matrices derived with the finite element method.

To show how the eigenfunction *f*_{1} can model the global shape of elongated structures, we construct its Reeb graph to obtain an explicit representation [22]. For a function *f*_{1} defined on a manifold , its Reeb graph is defined as the quotient space with its topology defined via the equivalent relation *x* *y* if *f*_{1}(*x*) = *f*_{1}(*y*) for *x*, *y* . To numerically construct the Reeb graph, we trace a set of level contours of the eigenfunction *f*_{1} on the triangular mesh representation of . To ensure the level contours distribute evenly over the entire surface, we use an adaptive sampling scheme developed in []. These level contours are used as the nodes of the Reeb graph and the connectivity of these nodes are established according to the neighboring relations of level contours. As an example, we show in Fig. 1(a) the Reeb graph of a hippocampus. By representing each node of the Reeb graph as the centroid of the level contour, we can see this graph has a chain structure and provides a compact model of the essentially one dimensional, tail-to-head trend of the hippocampus.

The Reeb graph of a hippocampus. (a) The Reeb graph of the *f*_{1}. (b) The surface patches interpolating the level contours.

Based on the Reeb graph of *f*_{1}, we define the first feature function *ξ*_{1}. Because the eigenfunction is generally a Moss function [23], its Reeb graph has a tree structure for genus zero surfaces. For hippocampus, the Reeb graph of *f*_{1} typically has a chain structure as shown in Fig. 1. In the case there are branches in the Reeb graph, we prune the smaller branches according to the size of the associated level contour to ensure the pruned graph has a chain structure. We order the nodes on this chain with the increase of the function *f*_{1} and denote the level contour at each node as *C _{i}*(

$$C=A\mathcal{V}$$

(15)

where *C* = [*C*_{1}, *C*_{2}, ···, *C _{N}*]

$${\left|\right|{\xi}_{1}(C)-A{\xi}_{1}(\mathcal{V})\left|\right|}^{2}+\beta {\xi}_{1}{(\mathcal{V})}^{T}Q{\xi}_{1}(\mathcal{V})$$

(16)

where *ξ*_{1}(*C*) and *ξ*_{1}() are vectors of the values of *ξ*_{1} on the level contours and the vertices of the mesh, respectively. The matrix *Q* is the same as in (14) and the term *ξ*_{1}()^{T} Qξ_{1}() encourages smoothness of the function *ξ*_{1}. For this least square problem, we have the solution for *ξ*_{1}() as

$${\xi}_{1}(\mathcal{V})={({A}^{T}A+\beta Q)}^{-1}{A}^{T}{\xi}_{1}(C).$$

(17)

We define the second feature function *ξ*_{2} also based on the eigenfunction of the Laplace-Beltrami operator to characterize the lateral profile of elongated structures. For each level contour *C _{i}*, we generate a surface patch that approximates the minimal surface of

$${\xi}_{2}(\mathcal{V})={({A}^{T}A+\beta Q)}^{-1}{A}^{T}{\xi}_{2}(C)$$

(18)

where *ξ*_{2}(*C*) and *ξ*_{2}() are the vectors of values of *ξ*_{2} on the level contours and the vertices, respectively.

As an illustration, we show the feature functions of a hippocampus in Fig: 2 (a) and (b) with the parameter *β* = 10, *N* = 100, *K* = 100, which clearly illustrates the power of the eigen-features in characterizing the relative locations of points on the surfaces.

In this section we present experimental results to demonstrate our method. We will first illustrate our algorithm on the mapping of two types of surfaces: hippocampus and putamen. After that, we apply our method to a data set of 890 hippocampal surfaces and present statistically significant results.

In the first experiment, we tested our algorithm on nine hippocampal surfaces. We chose one of them as the atlas and computed maps between this surface and the other eight surfaces. The parameters in the energy function are
${\alpha}_{D}^{1}=100,{\alpha}_{D}^{2}=50$, *α _{IC}* = 6,

A visualization of the hippocampal mapping results. The texture pattern on the atlas surface in the center is projected onto the other eight surfaces.

To illustrate the importance of the Laplace-Beltrami eigen-features and inverse consistency, we turned off these terms separately and measured quantitatively their impact on the quality of maps. In our experiments, we computed the standard deviation of |*J*_{u1}| and |*J*_{u2}| as an indicator of the quality of maps because it measures the area distortion resulting from the maps. As shown in Fig. 4(a) and (b), smaller area distortions have been achieved for all the surfaces by incorporating both the eigen-features and inverse consistency.

In the second experiment, we tested our method on a group of 30 putamen surfaces. By choosing one of the surfaces as the atlas, we computed the maps from the other 29 surfaces to the atlas with the same parameters as in the first experiment. Using the correspondences across the 30 putamen surfaces that have been established by the maps, we can build a shape prior model for the putamen by performing a principal component analysis (PCA) [2]. The mean shape of the putamen is shown at the top of Fig. 5(a), where a texture pattern is also generated for visualization. The eigenvalues computed from the PCA process are plotted at the bottom of Fig. 5(a) where we can see most of the variations are captured by the first few components. As an illustration of the shape prior model, we plotted the shapes generated with the first three principle components. The texture pattern on the mean shape is also plotted on each synthesized surface in Fig. 5(b) to visualize correspondences across the shapes. From the results we can see the shape prior model is able to generate valid shapes with a quite large range (±5*σ _{i}*) of coefficients for the principal components.

In the third experiment, we apply our mapping algorithm to a data set of 890 left hippocampal surfaces. The hippocampi were segmented automatically with the algorithm in [25] from the screening and 12-month follow-up scans of 445 subjects from the Alzheimer’s Disease Neuroimaging Initiative(ADNI)[26]. Among the subjects, there are 136 normal controls (NC), 228 patients with mild cognitive impairment(MCI), and 81 patients with Alzheimer’s disease(AD). The goal of this experiment is to test the differences of normal aging, MCI, and AD in terms of the atrophy rates in the hippocampus.

As a first step in this experiment, we computed maps from all 890 surfaces to the atlas surface used in the first experiment as shown in the center of Fig. 3, and the same parameters as in the first experiment were used for all surfaces. In order to compute such a large amount of surface maps, we implemented our method on a grid of more than 1000 CPUs and this enables us to compute all the maps in parallel.

To perform statistical analysis, we first generated a regular mesh representation of the atlas surface with 2000 vertices. Using the maps computed above, we projected the triangular mesh of the atlas onto the 890 surfaces. As a result, all surfaces were represented with the same triangulation and their vertices have one-to-one correspondences. To quantify the atrophy rate of the hippocampus locally, we define a thickness measure at each vertex of the mapped surfaces as the distance from the vertex to the Reeb graph of the first Laplace-Beltrami eigenfunction as shown in Fig. 1(a). By viewing this Reeb graph as a medial core of the hippocampus, our thickness measure is very similar to the definition in [1] except that our definition is completely intrinsic and invariant to pose variations such as translation, rotation and reflection. For the hippocampal surfaces of a subject, the atrophy rate at each corresponding vertex then equals the change of the thickness over the 12-month period as a percentage of the thickness at the screening scan.

Given the atrophy rates at corresponding vertices of all subjects, we performed two group analyses. In the first group study, we compared the NC group with the MCI group. At each vertex, a one-tailed t-test is applied to test the hypothesis that the MCI group has a higher atrophy rates than the NC group. The resulting p-value map is plotted onto the mean shape of the NC group in Fig. 6(a) and (b) in the top and bottom view. This map shows the significance level across different regions of the hippocampus. To correct for multiple comparisons, we performed 10,000 permutation tests and an overall statistically significance p-value 0.0063 was obtained. In the second group analysis, we compared the NC group with the AD group. Similarly, a one-tailed t-test was applied to test the hypothesis that the AD group has a higher atrophy rate than the NC group. The p-value map of this group analysis is plotted in Fig. 6(c) and (d). We also applied 10,000 permutation tests to correct for multiple comparisons and the overall p-value is 0.001, which shows the map is statistically significant. By comparing the maps in Fig. 6(c)(d) and Fig. 6(a)(b), we can clearly see the expansion of the regions with higher atrophy rates from the MCI to the AD group.

We proposed a general framework for surface mapping with both eigen-feature-based data terms and regularizing terms for inverse consistency and smoothness. We successfully demonstrated its application in computing maps between elongated structures such as hippocampus and putamen. For future work, we will design new eigen-features and apply the current framework to other anatomical structures. We will also perform validations to compare our method with previous techniques.

1. Thompson PM, Hayashi KM, de Zubicaray GI, Janke AL, Rose SE, Semple J, Hong MS, Herman DH, Gravano D, Doddrell DM, Toga AW. Mapping hippocampal and ventricular change in Alzheimer disease. NeuroImage. 2004;22(4):1754–1766. [PubMed]

2. Cootes T, Taylor C, Cooper D, Graham J. Active shape models-their training and application. Computer Vision and Image Understanding. 1995;61(1):38–59.

3. Christensen GE, Rabbitt RD, Miller MI. Deformable templates using large deformation kinematics. IEEE Trans Imag Process. 1996;5(10):1435–1447. [PubMed]

4. Joshi S, Miller MI. Landmark matching via large deformation diffeomorphisms. IEEE Trans Imag Process. 2000;9(8):1357–1370. [PubMed]

5. Gerig G, Styner M, Jones D, Weinberger D, Lieberman J. Shape analysis of brain ventricles using SPHARM. Proc Workshop on Mathematical Methods in Biomedical Image Analysis. 2001:171–178.

6. Davies RH, Twining CJ, Allen PD, Cootes TF, Taylor CJ. Shape discrimination in the hippocampus using an MDL model. Proc IPMI. 2003:38–50. [PubMed]

7. Gu X, Wang Y, Chan TF, Thompson PM, Yau ST. Genus zero surface conformal mapping and its application to brain surface mapping. IEEE Trans Med Imag. 2004;23(8):949–958. [PubMed]

8. Wang L, Miller JP, Gado MH, McKeel DW, Rothermich M, Miller MI, Morris JC, Csernansky JG. Abnormalities of hippocampal surface structure in very mild dementia of the alzheimer type. NeuroImage. 2006;30(1):52–60. [PMC free article] [PubMed]

9. Yeo B, Sabuncu M, Vercauteren T, Ayache N, Fischl B, Golland P. Spherical demons: Fast surface registration. Proc MICCAI. 2008:745–753. [PMC free article] [PubMed]

10. Shi Y, Thompson PM, de Zubicaray G, Rose SE, Tu Z, Dinov I, Toga AW. Direct mapping of hippocampal surfaces with intrinsic shape context. NeuroImage. 2007;37(3):792–807. [PMC free article] [PubMed]

11. Pizer SM, Fritsch DS, Yushkevich PA, Johnson VE, Chaney EL. Segmentation, registration, and measurement of shape variation via image object shape. IEEE Trans Med Imag. 1999;18(10):851–865. [PubMed]

12. Styner M, Gerig G, Lieberman J, Jones D, Weinberger D. Statistical shape analysis of neuroanatomical structures based on medial models. Med Image Anal. 2003;7(3):207–220. [PubMed]

13. Yushkevich PA, Zhang H, Gee JC. Continuous medial representation for anatomical structures. IEEE Trans Med Imag. 2006;25(12):1547–1564. [PubMed]

14. Christensen G, Johnson H. Consistent image registration. IEEE Trans Med Imag. 2001;20(7):568–582. [PubMed]

15. Eells J, Lemaire L. Two reports on harmonic maps. World Scientific. 1995

16. Mémoli F, Sapiro G, Osher S. Solving variational problems and partial differential equations mapping into general target manifolds. Journal of Computational Physics. 2004;195(1):263–292.

17. Osher S, Sethian J. Fronts propagation with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations. Journal of computational physics. 1988;79(1):12–49.

18. Bertalmío M, Cheng L, Osher S, Sapiro G. Variational problems and partial differential equations on implicit surfaces. Journal of Computational Physics. 2001;174(2):759–780.

19. Reuter M, Wolter F, Peinecke N. Laplace-Beltrami spectra as Shape-DNA of surfaces and solids. Computer-Aided Design. 2006;38:342–366.

20. Qiu A, Bitouk D, Miller MI. Smooth functional and structural maps on the neocortex via orthonormal bases of the Laplace-Beltrami operator. IEEE Trans Med Imag. 2006;25(10):1296–1306. [PubMed]

21. Belkin M, Niyogi P. Laplacian eigenmaps for dimensionality reduction and data representation. Neural Computation. 2003;15(6):1373–1396.

22. Reeb G. Sur les points singuliers d’une forme de Pfaff completement integrable ou d’une fonction nemérique. Comptes Rendus Acad Sciences. 1946;222:847–849.

23. Uhlenbeck K. Generic properties of eigenfunctions. Amer J of Math. 1976;98(4):1059–1078.

24. Shewchuk JR. Delaunay refinement algorithms for triangular mesh generation. Comput Geom Theory & Applications. 2002;22(1–3):21–74.

25. Morra J, Tu Z, Apostolova L, Green A, Avedissian C, Madsen S, Parikshak N, Hua N, Toga A, Jack C, Weiner M, Weiner M, Thompson P. the ADNI. Validation of a fully automated 3D hippocampal segmentation method using subjects with Alzheimer’s disease, mild cognitive impairment, and elderly controls. NeuroImage. in press. [PMC free article] [PubMed]

26. Mueller S, Weiner M, Thal L, Petersen R, Jack C, Jagust W, Trojanowski J, Toga A, Beckett L. The Alzheimer’s disease neuroimaging initiative. Clin North Am. 2005;15:869–877. xi–xii. [PMC free article] [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |