DBP2:UNC:Cortical Thickness Roadmap

Back to NA-MIC Internal Collaborations, UNC DBP 2

Brain development

Objective

We would like to create an end-to-end application within Slicer3 allowing individual and group analysis of regional and local cortical thickness. Such a workflow applied to the young brain (2-4 years old) is one goal of the UNC DBP. This page describes the technology roadmap for cortical thickness analysis in the NA-MIC Kit. The basic components necessary for this end-to-end application are:

• Tissue segmentation: Should be multi-modality, correcting for intensity inhomogeneity and work on non-skull-stripped data.
• Cortical thickness measurement: Local cortical thickness needs measurements at every location of the white-gray matter boundary, as well as at the gray-csf boundary. Regional analysis does not need such a dense measurement.
• Surface inflation applied on cortical data
• Cortical correspondence: Local analysis needs a full correspondence on both white-gray boundary and gray-csf boundary.
• Statistical analysis/Hypothesis testing: Measurements need to be compared and tested localy incorporating multiple-comparison correction, correlative analysis would be necessary too.

Roadmap

Starting with several MRI images (weighted-T1, weighted-T2...) we want to obtain cortical thickness maps for each subject, compute cortical correspondences between subjects, and analyze the cortical thickness at these corresponding locations. Ultimately, the NA-MIC Kit will provide a workflow for individual and group analysis of cortical thickness. It will be implemented as a set of Slicer3 modules that can be used interactively within the Slicer3 application as well as in batch on a computing cluster using BatchMake.

Next we discuss the main modules and details of current status and development work:

White matter/gray matter segmentation

Tissue segmentation

• UNC has an open source segmentation tool called itkEMS (binary download) implemented in an ITK framework for segmenting white matter and gray matter in the young brain. An initial Slicer3 command line module was implemented and works well. This initial Slicer3 command line module is based on the existing itkEMS xml I/O for parameters. A more extensive adaptation that would allow the direct modification of the itkEMS parameters within Slicer3 is in discussion.
  • Since this segmentation technique exists in an ITK framework, the integration into Slicer3 was low risk and has been completed.
• UNC also adapts the current Slicer3 EM Segment module to their young brain studies. Here, UNC adapted their UNC atlas of the 2 year old brain to provide priors for the EM Segment module. Initial tests have been successful, but more tests are needed.
  • This is a good test case for applying the Slicer3 EM Segment module to a slightly different application. UNC worked through the training material on the Slicer3 EM Segment module and discussed open issues regularly with Brad Davis (as well as Kilian Pohl).
• A tissue segmentation evaluation was performed to compare both algorithms: itkEMS and EMSegment

Local cortical thickness measurement

Cortical Thickness: white matter

Cortical Thickness: grey matter

• UNC has an algorithm to measure a local cortical thickness given a labeling of white matter and gray matter. This technique has been converted into a Slicer3 command line module.
  • This technique is non-symmetric and sparse (only computing distances where they can be computed reliably).
  • This module is well suited for regional analysis, but additional work to interpolate dense measurements along the boundaries from the sparse ones would be necessary for a local analysis.
  • A Slicer 3 external module has been implemented.
• Marc Niethammer developed a technique at a previous project week that is symmetric. We applied his tool in first tests successfully and this module seems to offer an alternative way to compute cortical thickness measurements for a local, symmetric analysis rather than regional analysis.

Local correspondence

• Regional as well as local subject comparisons are needed.
• Regional analysis will require precise deformable registration to a young brain atlas.
  • NA-MIC toolkit can be applied here, Slicer 3 has a b-spline based MI registration which needs to be tested.
    • The Slicer 3 registration modules have no possibility to save/load the computed transform (only the transformed image is saved), this is a necessity for this application!
    • This is an issue for both rigid, affine and b-spline based registration (all of these are needed for the full pipeline).
    • Testing will be performed to see whether registration quality is appropriate
    • An initial framework for the storage/retrieval and application of linear transforms has been implemented, and tests at UNC will soon take place.
  • Parcellation atlas for regional analysis is the current atlas used by the DBP.
• Local analysis requires techniques which are not currently in the NA-MIC Kit, specifically a point-to-point correspondence on the cortex (either white/gray or gray/csf) boundary needs to be established.
  • Freesurfer can be used for the local analysis (but it is not in the NA-MIC Kit).
  • Ipek Oguz is developing local analysis tools and should be available in Fall 2008. The tool necessitates the following steps (also performed within Freesurfer):
  1. Spherical topology of white matter: Marc Niethammer has Slicer module for this step, testing will be necessary to ensure appropriateness in the cortical case. Initial tests of that module were successful.
  2. Inflation of cortical surface: A Slicer module exists for this step, but likely not with the right properties. Tests and possible additional implementation will need to be done.
  3. Computation of correspondence on inflated surface using particle system: No module exists, a first command line prototype tool has been developed in collaboration with Utah (I Oguz, J Cates). Initial tests have shown to be very promising (ISBI 08 submission).

Hypothesis testing

Regional testing can be done with standard statistical tools as there are a limited, relatively small number of regions.

For local analysis, such standard statistical tools are not applicable. Currently there is no local hypothesis testing framework with multiple comparison and generalized linear model computation in Slicer 3. Freesurfer has such a framework readily available. The UNC shape analysis testing framework and extensions of it should be applicable here.

• Freesurfer can compute group difference and correlative analysis on local cortical thickness.
• UNC shape analysis testing framework allows local hypothesis testing with multiple comparison correction, but no General Linear Model. An extension in this regard is currently in development (M Styner). A Slicer 3 module allowing of the UNC shape analysis framework is now available. A more generic statistical testing module for Slicer 3 is in development.
• An alternative possible testing framework is also being developed at UNC (outside of NAMIC) as an open source tool.

Performance characterization and validation

• Characterize response based on signal noise, patient motion, etc.
• Comparison to other tools (FreeSurfer, itkEMS, UNC cortical thickness).

Schedule

• Tissue segmentation

• 09/11/2007 - Slicer registration modules input/output transforms (J Miller) -- Affine and rigid on TransformIO branch --
• 09/11/2007 - Slicer module to resample data through a transformation (J Miller)
• 09/11/2007 - B-spline transform nodes (A Yarmarkovich)
• 10/15/2007 - White matter/gray matter segmentation of the young brain using itkEMS as a Slicer3 module (C Vachet) -- 99% DONE --
• 01/01/2008 - White matter/gray matter segmentation of the young brain using the Slicer3 EM Segment module (C Vachet) -- DONE --
• 01/02/2008 - Separate use case scenario wiki page: how to adjust/adapt the atlas definition in the EM Segment module
• 01/07/2008 - Registration to atlas in EM Segment (B Davis) -- DONE --
• 02/15/2008 - Comparison study itkEMS vs EM Segment in Slicer 3 (C Vachet) -- DONE --

• Cortical thickness

• 12/23/2007 - UNC Cortical thickness measurement as a Slicer3 module (C Vachet) -- DONE --
• 01/15/2008 - Niethammer's Laplacian Cortical thickness measurement module code working at UNC as a Slicer3 module (C Vachet) - at AHM
• 03/01/2008 - Regional comparison between UNC and Laplacian cortical thickness (C Vachet)

• Cortical correspondence

• 03/15/2008 - Slicer 3 Deformable registration of young brain regional atlas (C Vachet) - after AHM
• 05/01/2008 - Regional analysis of cortical thickness as a Slicer3 module (C Vachet)
• 11/01/2007 - Niethammer Spherical topology tool tested on cortical data (I Oguz) -- DONE --
• 02/01/2007 - Inflation/Unfolding Slicer module tested on cortical data (I Oguz)
• 06/01/2008 - Slicer 3 module for cortical correspondence on inflated surface (I Oguz)
• 08/01/2008 - Fully evaluation of cortical computation on autism DBP datasets (C Vachet)

• Statistical analysis and Hypothesis testing

• 02/01/2008 - Slicer 3 Statistical Shape Analysis module (C Vachet, N Augier)
• 04/01/2008 - Slicer 3 Statistical Cortical Thickness Analysis module (C Vachet)
• 06/01/2008 - GLM model for UNC statistical shape analysis (M Styner, N Augier)
• xx/xx/xxxx - Nonparametric hypothesis testing in ITK (J Miller)

• Workflow/cohesive tools

• 06/01/2008 - Groupwise regional analysis of cortical thickness as a NA-MIC Workflow (C Vachet)
• 02/01/2009 - Groupwise local analysis of cortical thickness as a NA-MIC Workflow (C Vachet)

Team and Institute

• Co-PI: Heather Cody Hazlett, PhD, (heather_cody at med.unc.edu, Ph: 919-966-4099)
• Co-PI: Joseph Piven, MD
• NA-MIC Engineering Contact: Jim Miller, GE Research
• NA-MIC Algorithms Contact: Martin Styner, UNC
• Other team members
  • UNC DBP: Clement Vachet
  • UNC Algorithm: Ipek Oguz, Nicolas Augier, Marc Niethammer