Difference between revisions of "Projects:MultiscaleShapeAnalysis"

From NAMIC Wiki
Jump to: navigation, search
 
(22 intermediate revisions by 4 users not shown)
Line 1: Line 1:
Multiscale Shape Analysis applied to Caudate and Hippocampus
+
Back to [[NA-MIC_Internal_Collaborations:StructuralImageAnalysis|NA-MIC Collaborations]], [[Algorithm:Stony Brook|Stony Brook University Algorithms]], [[Algorithm:UNC|UNC Algorithms]], [[DBP1:Harvard|Harvard DBP1]]
  
Back to [[NA-MIC_Collaborations|NA-MIC_Collaborations]]
+
__NOTOC__
 
+
= Multiscale Shape Analysis applied to Caudate and Hippocampus =
'''Objective'''
 
  
 
We present a novel method of statistical surface-based morphometry based on the use of non-parametric permutation tests and a spherical wavelet (SWC) shape representation. As an application, we analyze two brain structures, the caudate nucleus and the hippocampus. We show that the results nicely complement the results obtained with shape analysis using a sampled point representation (SPHARM-PDM).
 
We present a novel method of statistical surface-based morphometry based on the use of non-parametric permutation tests and a spherical wavelet (SWC) shape representation. As an application, we analyze two brain structures, the caudate nucleus and the hippocampus. We show that the results nicely complement the results obtained with shape analysis using a sampled point representation (SPHARM-PDM).
  
''Overview''
+
= Description =
  
'Pre-processing'
+
''Pre-processing''
  
 
[[Image:UNCShape_OverviewAnalysis_MICCAI06.gif|thumb|right|300px|]]
 
[[Image:UNCShape_OverviewAnalysis_MICCAI06.gif|thumb|right|300px|]]
  
We use the [[Algorithm:UNC:Shape_Analysis#SPHARM-PDM|UNC Pipeline]] to pre-process the shapes. The input is a set of binary segmentation of a single brain structure. We use the following steps from the pipeline:
+
We use the [[Projects:ShapeAnalysisFrameworkUsingSPHARMPDM|UNC Pipeline]] to pre-process the shapes. The input is a set of binary segmentation of a single brain structure. We use the following steps from the pipeline:
  
 
# Morphological operations: fills any interior holes and applies a minimal smoothing operation
 
# Morphological operations: fills any interior holes and applies a minimal smoothing operation
Line 21: Line 20:
 
# Retriangulation: samples the parameterized and aligned surface meshes into triangulated surfaces (SPHARM-PDM) via icosahedron subdivision of the spherical parametrization. These SPHARM-PDM surfaces are all spatially aligned using rigid Procrustes alignment.
 
# Retriangulation: samples the parameterized and aligned surface meshes into triangulated surfaces (SPHARM-PDM) via icosahedron subdivision of the spherical parametrization. These SPHARM-PDM surfaces are all spatially aligned using rigid Procrustes alignment.
  
'Spherical Wavelet Features'
+
''Spherical Wavelet Features''
  
 
[[Image:Basis_membership.png|thumb|right|200px|Figure 2: Visualization of spherical wavelet basis functions and associated regions at three levels (columns). Top row : Values of single spherical Wavelet Basis Function shown on the sphere at scales 1 through 3.  Middle and Bottom row:  Regions of influence of the spherical wavelet basis functions shown on the sphere and on the original surface, each basis function region has a random color.]]
 
[[Image:Basis_membership.png|thumb|right|200px|Figure 2: Visualization of spherical wavelet basis functions and associated regions at three levels (columns). Top row : Values of single spherical Wavelet Basis Function shown on the sphere at scales 1 through 3.  Middle and Bottom row:  Regions of influence of the spherical wavelet basis functions shown on the sphere and on the original surface, each basis function region has a random color.]]
Line 27: Line 26:
 
For each triangulated SPHARM-PDM surface (and its corresponding spherical parameterization), a spherical wavelet description is computed. As a result, each shape is represented by a series of 3D spherical wavelet coefficients (SWC). Each 3D coefficient is associated with a basis function that describes a region of the surface. The size of that region depends on the scale of the coefficient. Therefore each 3D coefficient describes the shape at a specific scale and spatial region (See Figure 2).
 
For each triangulated SPHARM-PDM surface (and its corresponding spherical parameterization), a spherical wavelet description is computed. As a result, each shape is represented by a series of 3D spherical wavelet coefficients (SWC). Each 3D coefficient is associated with a basis function that describes a region of the surface. The size of that region depends on the scale of the coefficient. Therefore each 3D coefficient describes the shape at a specific scale and spatial region (See Figure 2).
  
'Statistics'
+
''Statistics''
  
We use the UNC [[Algorithm:UNC:Shape_Analysis#Statistics|statistical test toolbox]] that analyzes differences between two groups of surfaces described by a set of features. The group differences are computed locally for every feature using the standard robust Hotelling T^2 two sample metric. Statistical p-values, both raw and corrected for multiple comparisons are given as output. The toolbox outputs a global average P-value for all features, as well as a raw and corrected P-value for each feature. We use the 3D spherical wavelet coefficients as features.
+
We use the UNC [[Projects:LocalStatisticalAnalysisViaPermutationTests |statistical test toolbox]] that analyzes differences between two groups of surfaces described by a set of features. The group differences are computed locally for every feature using the standard robust Hotelling T^2 two sample metric. Statistical p-values, both raw and corrected for multiple comparisons are given as output. The toolbox outputs a global average P-value for all features, as well as a raw and corrected P-value for each feature. We use the 3D spherical wavelet coefficients as features.
  
'Significance Maps'
+
''Significance Maps''
  
 
If a feature (3D spherical wavelet coefficient) is found significant (i.e. its P-value is less than a pre-determined significance level, such as 0.05), we color all points that are in the support of that basis function at that scale with the corresponding P-value. This allows us to visualize both the raw and FDR corrected P-values as significance color maps on the surface of the mean shape of the structure under study. The color at each point is the P-value.
 
If a feature (3D spherical wavelet coefficient) is found significant (i.e. its P-value is less than a pre-determined significance level, such as 0.05), we color all points that are in the support of that basis function at that scale with the corresponding P-value. This allows us to visualize both the raw and FDR corrected P-values as significance color maps on the surface of the mean shape of the structure under study. The color at each point is the P-value.
  
'''Progress'''
+
''Progress''
  
 
We conducted statistical shape analysis of two brain structures, the caudate nucleus and hippocampus, using spherical wavelet coefficients (SWC) as features and compare the results obtained to shape analysis using a SPHARM-PDM representation.
 
We conducted statistical shape analysis of two brain structures, the caudate nucleus and hippocampus, using spherical wavelet coefficients (SWC) as features and compare the results obtained to shape analysis using a SPHARM-PDM representation.
  
''References''
+
= Key Investigators =
 
 
* D. Nain, M. Styner, M. Niethammer, J. J. Levitt, M E Shenton, G Gerig, A. Bobick, A. Tannenbaum. Statistical Shape Analysis of Brain Structures using Spherical Wavelets. Submitted to ISBI 2007.
 
* M. Styner, I. Oguz, S. Xu, C. Brechbuehler, D. Pantazis, J. Levitt, M. Shenton, G. Gerig. Framework for the Statistical Shape Analysis of Brain Structures using SPHARM-PDM. MICCAI 2006 Open Source Workshop. Insight Journal DSpace link: http://hdl.handle.net/1926/215
 
 
 
* Y. Gao, D. Nain, X. LeFaucheur, A. Tannenbaum. A Spherical Wavelet ITK Filter. In preparation for the Insight Journal, December 2006.
 
* Y. Gao, D. Nain, M. Styner, M. Niethammer, J. J. Levitt, M E Shenton, G Gerig, A. Bobick, A. Tannenbaum. Framework for the Statistical Shape Analysis of Brain Structures using Spherical Wavelets. In preparation for the Insight Journal, February/March 2007.
 
 
 
'''Key Investigators'''
 
  
* Core 1:
+
* Georgia Tech Algorithms: Delphine Nain, Yi Gao, Xavier Le Faucheur, Allen Tannenbaum
** Georgia Tech: Delphine Nain, Yi Gao, Xavier Le Faucheur, Allen Tannenbaum
+
* UNC Algorithms: Martin Styner
** UNC: Martin Styner
+
* Harvard DBP1: James Levitt, Marc Niethammer, Sylvain Bouix, Martha Shenton
* Core 3:
 
** Harvard PNL: James Levitt, Marc Niethammer, Sylvain Bouix, Martha Shenton
 
  
'''Links'''
+
= Publications =
 +
* [http://www.na-mic.org/publications/pages/display?search=MultiscaleShapeAnalysis NA-MIC Publications Database on Multiscale Shape Analysis]
  
* [[Algorithm:UNC:Shape_Analysis|UNC Shape Analysis]]
+
[[Category:Shape Analysis]] [[Category:Schizophrenia]]
* [[Algorithm:GATech|Georgia Tech Summary Page]]
 
* [[NA-MIC_Collaborations|NA-MIC_Collaborations]]
 

Latest revision as of 00:57, 16 November 2013

Home < Projects:MultiscaleShapeAnalysis
Back to NA-MIC Collaborations, Stony Brook University Algorithms, UNC Algorithms, Harvard DBP1


Multiscale Shape Analysis applied to Caudate and Hippocampus

We present a novel method of statistical surface-based morphometry based on the use of non-parametric permutation tests and a spherical wavelet (SWC) shape representation. As an application, we analyze two brain structures, the caudate nucleus and the hippocampus. We show that the results nicely complement the results obtained with shape analysis using a sampled point representation (SPHARM-PDM).

Description

Pre-processing

UNCShape OverviewAnalysis MICCAI06.gif

We use the UNC Pipeline to pre-process the shapes. The input is a set of binary segmentation of a single brain structure. We use the following steps from the pipeline:

  1. Morphological operations: fills any interior holes and applies a minimal smoothing operation
  2. Surface conversion: converts the processed binary segmentations to surface meshes
  3. Spherical Parameterization: computes a spherical parameterization for the surface meshes using a area-preserving, distortion minimizing spherical mapping
  4. Alignment: computes a spherical harmonic (SPARM) description of the parameterized surface meshes and aligns the spherical parameterizations using the first order ellipsoid from the spherical harmonic coefficients.
  5. Retriangulation: samples the parameterized and aligned surface meshes into triangulated surfaces (SPHARM-PDM) via icosahedron subdivision of the spherical parametrization. These SPHARM-PDM surfaces are all spatially aligned using rigid Procrustes alignment.

Spherical Wavelet Features

Figure 2: Visualization of spherical wavelet basis functions and associated regions at three levels (columns). Top row : Values of single spherical Wavelet Basis Function shown on the sphere at scales 1 through 3. Middle and Bottom row: Regions of influence of the spherical wavelet basis functions shown on the sphere and on the original surface, each basis function region has a random color.

For each triangulated SPHARM-PDM surface (and its corresponding spherical parameterization), a spherical wavelet description is computed. As a result, each shape is represented by a series of 3D spherical wavelet coefficients (SWC). Each 3D coefficient is associated with a basis function that describes a region of the surface. The size of that region depends on the scale of the coefficient. Therefore each 3D coefficient describes the shape at a specific scale and spatial region (See Figure 2).

Statistics

We use the UNC statistical test toolbox that analyzes differences between two groups of surfaces described by a set of features. The group differences are computed locally for every feature using the standard robust Hotelling T^2 two sample metric. Statistical p-values, both raw and corrected for multiple comparisons are given as output. The toolbox outputs a global average P-value for all features, as well as a raw and corrected P-value for each feature. We use the 3D spherical wavelet coefficients as features.

Significance Maps

If a feature (3D spherical wavelet coefficient) is found significant (i.e. its P-value is less than a pre-determined significance level, such as 0.05), we color all points that are in the support of that basis function at that scale with the corresponding P-value. This allows us to visualize both the raw and FDR corrected P-values as significance color maps on the surface of the mean shape of the structure under study. The color at each point is the P-value.

Progress

We conducted statistical shape analysis of two brain structures, the caudate nucleus and hippocampus, using spherical wavelet coefficients (SWC) as features and compare the results obtained to shape analysis using a SPHARM-PDM representation.

Key Investigators

  • Georgia Tech Algorithms: Delphine Nain, Yi Gao, Xavier Le Faucheur, Allen Tannenbaum
  • UNC Algorithms: Martin Styner
  • Harvard DBP1: James Levitt, Marc Niethammer, Sylvain Bouix, Martha Shenton

Publications