Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2006:(1071):242-250.

Framework for the Statistical Shape Analysis of Brain Structures using SPHARM-PDM

Framework for the Statistical Shape Analysis of Brain Structures using SPHARM-PDM

Martin Styner et al. Insight J. 2006.

Abstract

Shape analysis has become of increasing interest to the neuroimaging community due to its potential to precisely locate morphological changes between healthy and pathological structures. This manuscript presents a comprehensive set of tools for the computation of 3D structural statistical shape analysis. It has been applied in several studies on brain morphometry, but can potentially be employed in other 3D shape problems. Its main limitations is the necessity of spherical topology.The input of the proposed shape analysis is a set of binary segmentation of a single brain structure, such as the hippocampus or caudate. These segmentations are converted into a corresponding spherical harmonic description (SPHARM), which is then sampled into a triangulated surfaces (SPHARM-PDM). After alignment, differences between groups of surfaces are computed using the Hotelling T(2) two sample metric. Statistical p-values, both raw and corrected for multiple comparisons, result in significance maps. Additional visualization of the group tests are provided via mean difference magnitude and vector maps, as well as maps of the group covariance information.The correction for multiple comparisons is performed via two separate methods that each have a distinct view of the problem. The first one aims to control the family-wise error rate (FWER) or false-positives via the extrema histogram of non-parametric permutations. The second method controls the false discovery rate and results in a less conservative estimate of the false-negatives.

PubMed Disclaimer

Figures

Figure 1
Figure 1
Schematic view of our SPHARM-PDM based shape analysis pipeline.
Figure 2
Figure 2
The SPHARM shape description of a human lateral ventricle shown at 4 different degrees ( 1, 3, 6, 10 harmonics).
Figure 3
Figure 3
Left: Illustration of SPHARM correspondence via alignment of the spherical parametrization using the first order ellipsoid meridian and equator. Right: Set of 6 caudate structures with correspondence shown at selected locations via colored spheres.
Figure 4
Figure 4
Icosahedron subdivision for different levels of sampling. From left to right: Base icosahedron, subdivision factors 2, 4 and 6.
Figure 5
Figure 5
Left: Schematic view of testing procedure: A set of features per surface point (shown as colormaps) are analyzed using uni- or multi-variate statistics, resulting in a raw significance map (bottom) .Right: The raw significance maps are overly optimistic and need to be corrected to control for the multiple comparison problem. The corrected significance maps on the other hand are commonly pessimistic estimates.
Figure 6
Figure 6
Left: Scheme of P-value computation via permutation tests. The real group difference S0 (e.g. Hotelling T2) is compared to that the group differences Sj of random permutations of the group labels. The quantile in the Sj histogram associated with S0 is the p-value. Right: Example run on test data with 1000 random permutations. On the left the individual Sj are plotted and on the left the corresponding Sj histogram is shown. The black line indicates the value of S0.
Figure 7
Figure 7
Workflow of the shape analysis pipeline starting at the segmentation of the brain structures to the local final statistical analysis.
Figure 8
Figure 8
Quality Control visualization of the SPHARM correspondence with KWMeshVisu using the Phi-attribute file shown on 3 randomly selected cases. Same color represent the same ϕ parameter value of the spherical parametrization.
Figure 9
Figure 9
Descriptive statistics of right caudate analysis. The left visualization shows the mean difference as a colormap from green (0 mm difference) to red (1.6 mm difference) as well an vectorfield on the overall mean surface mesh. For two regions with large differences a zoomed view is provided. The right visualization shows the covariance ellipsoid colored with the magnitude of the largest ellipsoid axis (axis of largest covariance). The same two zoomed views are provided. The variance in the first zoom section is markedly reduced compared to the second zoom section.
Figure 10
Figure 10
Significance testing on right caudate study

References

    1. Thomson D. On Growth and Form. second edition. Cambridge University Press; 1942. 1.
    1. Davatzikos C, Vaillant M, Resnick S, Prince JL, Letovsky S, Bryan RN. A computerized method for morphological analysis of the corpus callosum. J. of Comp. Assisted Tomography. 1996 Jan./Feb;20:88–97. 1. - PubMed
    1. Joshi S, Miller M, Grenander U. On the geometry and shape of brain sub-manifolds. Pat. Rec. Art. Intel. 1997;11:1317–1343. 1.
    1. Csernansky JG, Joshi S, Wang LE, Haller J, Gado M, Miller JP, Grenander U, Miller MI. Hippocampal morphometry in schizophrenia via high dimensional brain mapping. Proc. Natl. Acad. Sci. USA. 1998 September;95:11406–11411. 1. - PMC - PubMed
    1. Csernansky JG, Joshi S, Wang LE, Haller J, Gado M, Miller JP, Grenander U, Miller MI. Hippocampal morphometry in schizophrenia via high dimensional brain mapping. Proc. Natl. Acad. Sci. USA. 1998 September;95:11406–11411. 1. - PMC - PubMed

LinkOut - more resources