Difference between revisions of "Hageman:NAMICFluidMechDTITractography"

From NAMIC Wiki
Jump to: navigation, search
 
(16 intermediate revisions by the same user not shown)
Line 4: Line 4:
  
 
== Overview ==
 
== Overview ==
In this paper, we present a novel method modeled on the dynamics of a viscous fluid described by the second order non-linear Navier-Stokes equations, an extension of previous work by our group \cite{hageman-toga:2006}  The Navier-Stokes equations govern flow for a viscous Newtonian fluid, a fluid whose internal stress forces are linearly related to its strain forces \cite{versteeg}However, even though these equations are most commonly seen in the context of fluid mechanics, they have been shown to be successful in modeling a large number of diverse physical phenomena.  Our second order nonlinear-based approach is an extension of previous linear PDE methods.  We hypothesize that our method will have a better performance than competing tractography methods in the presence of noise and it will not require a white matter mask, commonly used by other PDE-based methods.  This is because our model contains a viscous force not present in previous methods, represented as an additional convection term in the PDE.  We model local viscosity of the fluid as a function of the local intervoxel and intravoxel anisotropy in the corresponding DTI image volume.  The incorporation of this convection term in our flow field calculation allows us to closely couple the magnitude of the fluid velocity to the magnitude of the underlying anisotropy of the DTI tensor field, providing a dampening force in background areas, such as gray matter and CSF.  This eliminates the need for the white matter mask used by other PDE-based methods to prevent the model from entering these areas.  To compute an estimate of the most likely connection path between two regions in the brain, we simulate the flow of an artificial fluid between those two points through a volume whose dimensions, pressure, and local viscosity are derived from the underlying DTI data.  We then numerically solve for the fluid velocity vector field.  The estimated connection path is then computed by finding the optimal path through the fluid velocity that simultaneously maximizes both the fluid velocity and its gradient.  We apply our method to DTI volumes from normal human subjects and to computer-generated DTI phantoms to validate our approach.    
+
 
 +
Currently, our project is focused on developing a novel method for diffusion tensor imaging (DTI) tractography modeled on the dynamics of a viscous fluid described by the second order non-linear Navier-Stokes equations.  Even though these equations are most commonly seen in the context of fluid mechanics, they have been shown to be successful in modeling a large number of diverse physical phenomena.  Our second order nonlinear-based approach is an extension of previous linear PDE methods, but our model contains a viscous force not present in previous methods, represented as an additional convection term in the PDE.  We model local viscosity of the fluid as a function of the local intervoxel and intravoxel anisotropy in the corresponding DTI image volume.  The incorporation of this convection term in our flow field calculation allows us to closely couple the magnitude of the fluid velocity to the magnitude of the underlying anisotropy of the DTI tensor field, providing a dampening force in background areas, such as gray matter and CSF.  This eliminates the need for the white matter mask used by other PDE-based methods to prevent the model from entering these areas.  To compute an estimate of the most likely connection path between two regions in the brain, we simulate the flow of an artificial fluid between those two points through a volume whose dimensions, pressure, and local viscosity are derived from the underlying DTI data.  We then numerically solve for the fluid velocity vector field.  The estimated connection path is then computed by finding the optimal path through the fluid velocity that simultaneously maximizes both the fluid velocity and its gradient.   
 +
 
 +
Computational fluid dynamics is a rich field and, in addition to this tractography method, we are investigating its application to the analysis of diffusion tensor imaging (DTI) datasets for registration and analysis of white matter pathology.  We are currently developing several useful and novel diffusion tensor imaging (DTI) analysis algorithms modeled on the principles of fluid mechanics for inclusion within the NA-MIC framework.  The goal of this project is to fully develop these methods, make them compatible with the NA-MIC ITK-based software infrastructure (i.e. Slicer), and promote their dissemination to the scientific community.
  
 
== Description ==
 
== Description ==
  
Descriptions and progress go hereHere is where you write the bulk of your project.  The details given in this Template Project page represent the minimal amount of information that should be included for each projectBy housing all of the Project information on a single page, Wiki maintanence is reducedNow we can link to this page from several other places in the Wiki, such as each site's Site Page, the NA-MIC Collaborations page, Project Week pages, and from the externally visible NA-MIC webpage.
+
Our algorithm simulates the flow of an artificial fluid through a volume whose dimensions, viscosity, and pressure tensor field are derived from a DTI volumeSpecific regions of interest are chosen as sources and/or sinks, and we simulate the flow of an artificial fluid governed by the Navier-Stokes equations.  The most likely connection path is then estimated using a generalized gradient vector flow (GGVF) based approach to compute the trajectory through the fluid velocity vector field that simultaneously maximizes the magnitude of the fluid velocity and its gradient along the path.  Our fluid model is valid only as a theoretical framework for generating a connectivity metric and does not try to model any aspect of the underlying diffusion process.
 +
 
 +
[[Image:NAMICFMech_VisFig.jpg|Viscosity Map]]
 +
 
 +
Viscosity maps derived from 2D slices of DTI data from human control subjects.  Viscosity values were calculated from the corresponding diffusion tensor image and are color-coded according to the legend bar seen on the right side of the figureA. Axial slice taken at the level of the internal capsule.  The corpus callosum, marked with a star is a highly organized white matter tract and is therefore characterized by low viscosity.  Conversely, the lateral ventricle, marked with a delta contains CSF and is highly viscous.  BA mid-sagittal slice.  As in A, the corpus callosum is marked with a star and is characterized by low viscosity.  In contrast, the lateral ventricular space, marked with a delta contains CSF and therefore has no architecture.  Consequently, it is highly viscous.
  
''A Special Note: Dashed Boxes''
+
[[Image:NAMICFMech_PressureFig.jpg|Pressure Map]]
  
If a blank space is the first character for a Wiki line, the following text will be listed in a dashed box, such as the one appearing at the top of each Site and Project Page containing the links to go back to previous pagesNote, that if you entered the edit mode via any of the subheading edit links, upon a "Save page" the blank space at the top of the Wiki page will be removedTherefore, you will need to manually add the blank space to maintain the dashed text box at the top of the pageDashed boxes are also used at the bottom of the Project Pages to link to applicable Project Week pages.
+
Representations of the pressure tensor derived from 2D slices from human control subjects.  At each voxel in the image, the pressure force is represented by a tensor glyph, an ellipsoid whose axis is obtained from a diagonalization of the corresponding pressure tensor.  The color of the ellipsoid represents the dominant diffusion direction, according to the color coded axes in the figure with the superior-inferior z-axis (blue) coming out of the page, the anterior-posterior y-axis (green) vertical, and the left-right x-axis (red) horizontal.  A. A 2D axial slice taken at the level of the internal capsuleThe white box marks the enlarged area shown in B.  B.  Enlarged view from A.  The posterior limb of the corpus callosum, marked with a star is a highly organized white matter tract, and the pressure force acts on the artificial fluid co-linear with the fiber tractIn contrast, the lateral ventricular space, marked with a delta contains little structureConsequently, the pressure force is isotropic in that region.
  
''A Special Note about Publications''
+
[[Image:NAMICflowfieldlesion.JPG|400 px|Fluid Velocity Field Solution]]
  
"In Press" publications will be listed in two places on the Wiki (in the "In Press" sections on the Publications page and on the appropriate Project page(s))"In Press" publications will not be added to the Publications DB.
+
We solve our modified Navier-Stokes fluid mechanics model with these variables to get a fluid velocity field which is then used as a metric of regional connectivityTracts are generated using a modified method based on a generalized gradient vector field approach.
  
"In Print" publications will be added to the Publications DB with the appropriate Grant number (linking it to the search on the Publications DB) and with the appropriate Project page URL added as a keyword (linking it into the Project page).  Once the publication is added to the DB, it will be removed from it's listing in the two "In Press" sections mentioned above.
+
[[Image:NAMICFMech_CspFig.jpg|Human Corticospinal Tracts]]
  
Write access for the publications database will be given out to a few select individualsIf you do not currently have write access to the publications database, please send the following information to John Melonakos at (jmelonak at ece.gatech.edu), for "In Print" publications only: 1) the final printed version of the PDF, 2) the bibtex entry for the publication, 3) the NA-MIC Project URLs against which this publication should be linked, and 4) any other files that you would like to add to the publication database for this publication.
+
Segmentation results of the corticospinal tracts using our method in human control DTI images.  ROIs were placed within the brainstem below the level of the crossing pontine fibers and within the corona radiata above the level of the corpus callosumCross-sections of the approximate location of these ROIs are shown Figure 8D, drawn in white, superimposed on the corresponding axial DEC slices. Figure 8A shows the estimated connection paths between these ROIs generated by our method. An axial slice of the tensor glyphs at the level of the mid-brain is shown for spatial reference. The tracts show a prominent lateral course at the level of the mid-pons (Figure 8A: white arrow).  This corresponds to a strong lateral diffusion component at that point as seen in the directionally encoded color (DEC) image of the axial slice at the mid-pontine level (Figure 8B, 8C: white arrow).  The tensor glyphs and DEC image are color-coded by the axes shown in the figure.
 +
 
 +
[[Image:Hageman_FullBrainSlicerTractography.jpg|Full brain tracts segmented using our fluid mechanics based tractography method.]]
  
 
== Key Investigators ==
 
== Key Investigators ==
Line 31: Line 40:
 
''In Print''
 
''In Print''
 
* [http://www.na-mic.org/Special:Publications?text=Projects:MultiscaleShapeSegmentation&submit=Search&keywords=checked NA-MIC Publications Database].
 
* [http://www.na-mic.org/Special:Publications?text=Projects:MultiscaleShapeSegmentation&submit=Search&keywords=checked NA-MIC Publications Database].
 +
* Hageman NS, Toga AW, Narr K, Shattuck DW (2009).  A diffusion tensor imaging tractography algorithm based on Navier-Stokes fluid mechanics. IEEE Trans. in Medicial Imaging, 28(3): 348-360.
  
 
<!--''In Press''-->
 
<!--''In Press''-->
  
 
''In Submission''
 
''In Submission''
* Hageman NS, Toga AW, Narr K, Shattuck DW (2008).  A diffusion tensor imaging tractography algorithm based on Navier-Stokes fluid mechanics. IEEE Trans. in Medicial Imaging, In Submission.
+
* Hamilton L, Nuechterlein K, Hageman NS, Woods R, Asarnow R, Alger J, Gaser C, Toga AW, Narr K (2008).  Mean Diffusivity and Fractional Anisotropy as Indicators of Schizophrenia and Genetic Vulnerability, ''Neuroimage'', In Submission.
* Hamilton L, Nuechterlein K, Hageman NS, Woods R, Asarnow R, Alger J, Gaser C, Toga AW, Narr K (2008).  Mean Diffusivity and Fractional Anisotropy as Indicators of Schizophrenia and Genetic Vulnerability, Neuroimage, In Submission.
 
 
 
  
 
== Links ==
 
== Links ==
Line 44: Line 52:
 
* [http://www.loni.ucla.edu LONI Website]
 
* [http://www.loni.ucla.edu LONI Website]
  
  Project Week Results:  [[2008_Winter_Project_Week:Fluid_Mechanics_Tractography|2008 Winter]]
+
  Project Week Results:  [[2008_Winter_Project_Week:Fluid_Mechanics_Tractography|2008 Winter]], [[2008_Summer_Project_Week|2008 Summer]], [[2009_Winter_Project_Week|2009 Winter]], [[2009_Summer_Project_Week|2009 Summer]], [[2010_Summer_Project_Week|2010 Summer]]

Latest revision as of 18:11, 22 June 2010

Home < Hageman:NAMICFluidMechDTITractography


Fluid Mechanics Based DTI Tractography

Overview

Currently, our project is focused on developing a novel method for diffusion tensor imaging (DTI) tractography modeled on the dynamics of a viscous fluid described by the second order non-linear Navier-Stokes equations. Even though these equations are most commonly seen in the context of fluid mechanics, they have been shown to be successful in modeling a large number of diverse physical phenomena. Our second order nonlinear-based approach is an extension of previous linear PDE methods, but our model contains a viscous force not present in previous methods, represented as an additional convection term in the PDE. We model local viscosity of the fluid as a function of the local intervoxel and intravoxel anisotropy in the corresponding DTI image volume. The incorporation of this convection term in our flow field calculation allows us to closely couple the magnitude of the fluid velocity to the magnitude of the underlying anisotropy of the DTI tensor field, providing a dampening force in background areas, such as gray matter and CSF. This eliminates the need for the white matter mask used by other PDE-based methods to prevent the model from entering these areas. To compute an estimate of the most likely connection path between two regions in the brain, we simulate the flow of an artificial fluid between those two points through a volume whose dimensions, pressure, and local viscosity are derived from the underlying DTI data. We then numerically solve for the fluid velocity vector field. The estimated connection path is then computed by finding the optimal path through the fluid velocity that simultaneously maximizes both the fluid velocity and its gradient.

Computational fluid dynamics is a rich field and, in addition to this tractography method, we are investigating its application to the analysis of diffusion tensor imaging (DTI) datasets for registration and analysis of white matter pathology. We are currently developing several useful and novel diffusion tensor imaging (DTI) analysis algorithms modeled on the principles of fluid mechanics for inclusion within the NA-MIC framework. The goal of this project is to fully develop these methods, make them compatible with the NA-MIC ITK-based software infrastructure (i.e. Slicer), and promote their dissemination to the scientific community.

Description

Our algorithm simulates the flow of an artificial fluid through a volume whose dimensions, viscosity, and pressure tensor field are derived from a DTI volume. Specific regions of interest are chosen as sources and/or sinks, and we simulate the flow of an artificial fluid governed by the Navier-Stokes equations. The most likely connection path is then estimated using a generalized gradient vector flow (GGVF) based approach to compute the trajectory through the fluid velocity vector field that simultaneously maximizes the magnitude of the fluid velocity and its gradient along the path. Our fluid model is valid only as a theoretical framework for generating a connectivity metric and does not try to model any aspect of the underlying diffusion process.

Viscosity Map

Viscosity maps derived from 2D slices of DTI data from human control subjects. Viscosity values were calculated from the corresponding diffusion tensor image and are color-coded according to the legend bar seen on the right side of the figure. A. Axial slice taken at the level of the internal capsule. The corpus callosum, marked with a star is a highly organized white matter tract and is therefore characterized by low viscosity. Conversely, the lateral ventricle, marked with a delta contains CSF and is highly viscous. B. A mid-sagittal slice. As in A, the corpus callosum is marked with a star and is characterized by low viscosity. In contrast, the lateral ventricular space, marked with a delta contains CSF and therefore has no architecture. Consequently, it is highly viscous.

Pressure Map

Representations of the pressure tensor derived from 2D slices from human control subjects. At each voxel in the image, the pressure force is represented by a tensor glyph, an ellipsoid whose axis is obtained from a diagonalization of the corresponding pressure tensor. The color of the ellipsoid represents the dominant diffusion direction, according to the color coded axes in the figure with the superior-inferior z-axis (blue) coming out of the page, the anterior-posterior y-axis (green) vertical, and the left-right x-axis (red) horizontal. A. A 2D axial slice taken at the level of the internal capsule. The white box marks the enlarged area shown in B. B. Enlarged view from A. The posterior limb of the corpus callosum, marked with a star is a highly organized white matter tract, and the pressure force acts on the artificial fluid co-linear with the fiber tract. In contrast, the lateral ventricular space, marked with a delta contains little structure. Consequently, the pressure force is isotropic in that region.

Fluid Velocity Field Solution

We solve our modified Navier-Stokes fluid mechanics model with these variables to get a fluid velocity field which is then used as a metric of regional connectivity. Tracts are generated using a modified method based on a generalized gradient vector field approach.

Human Corticospinal Tracts

Segmentation results of the corticospinal tracts using our method in human control DTI images. ROIs were placed within the brainstem below the level of the crossing pontine fibers and within the corona radiata above the level of the corpus callosum. Cross-sections of the approximate location of these ROIs are shown Figure 8D, drawn in white, superimposed on the corresponding axial DEC slices. Figure 8A shows the estimated connection paths between these ROIs generated by our method. An axial slice of the tensor glyphs at the level of the mid-brain is shown for spatial reference. The tracts show a prominent lateral course at the level of the mid-pons (Figure 8A: white arrow). This corresponds to a strong lateral diffusion component at that point as seen in the directionally encoded color (DEC) image of the axial slice at the mid-pontine level (Figure 8B, 8C: white arrow). The tensor glyphs and DEC image are color-coded by the axes shown in the figure.

Full brain tracts segmented using our fluid mechanics based tractography method.

Key Investigators

  • UCLA: Nathan Hageman
  • UCLA: Arthur Toga, Ph.D

Publications

In Print

  • NA-MIC Publications Database.
  • Hageman NS, Toga AW, Narr K, Shattuck DW (2009). A diffusion tensor imaging tractography algorithm based on Navier-Stokes fluid mechanics. IEEE Trans. in Medicial Imaging, 28(3): 348-360.


In Submission

  • Hamilton L, Nuechterlein K, Hageman NS, Woods R, Asarnow R, Alger J, Gaser C, Toga AW, Narr K (2008). Mean Diffusivity and Fractional Anisotropy as Indicators of Schizophrenia and Genetic Vulnerability, Neuroimage, In Submission.

Links

  • Hageman NS, Shattuck DW, Narr K, Toga AW (2006). A diffusion tensor imaging tractography method based on Navier-Stokes fluid mechanics. Proceedings of the 2006 IEEE International Symposium on Biomedical Imaging: From Nano to Macro (ISBI 2006), Arlington, VA, USA, 6-9 April 2006. p. 798-801 (PDF)
  • LONI Website
Project Week Results:  2008 Winter, 2008 Summer, 2009 Winter, 2009 Summer, 2010 Summer