Difference between revisions of "Documentation/4.1/Modules/BRAINSMush"

From Slicer Wiki
Jump to: navigation, search
(Prepend documentation/versioncheck template. See http://na-mic.org/Mantis/view.php?id=2887)
 
(57 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 +
<noinclude>{{documentation/versioncheck}}</noinclude>
 
[[Documentation-3.6|Return to Slicer 3.6 Documentation]]
 
[[Documentation-3.6|Return to Slicer 3.6 Documentation]]
  
Line 9: Line 10:
 
{|
 
{|
 
|[[Image:screenshotBlankNotOptional.png|thumb|280px|User Interface]]
 
|[[Image:screenshotBlankNotOptional.png|thumb|280px|User Interface]]
|[[Image:screenshotBlank.png|thumb|280px|Output]]
+
|[[Image:T1.png|thumb|280px|T1 Image]]
|[[Image:screenshotBlank.png|thumb|280px|Caption]]
+
|[[Image:T2.png|thumb|280px|T2 Image]]
 +
|}
 +
 
 +
{|
 +
|[[Image:Mush.png|thumb|280px|Mush Image]]
 +
|[[Image:MushMask.png|thumb|280px|Mush Brain Extraction Mask(Coronal)]]
 +
|[[Image:MushMask2.png|thumb|280px|Mush Brain Extraction Mask(Sagittal)]]
 
|}
 
|}
  
Line 26: Line 33:
 
* Steven A. Dunn: University of Iowa
 
* Steven A. Dunn: University of Iowa
 
* Vincent Magnotta, PhD: University of Iowa
 
* Vincent Magnotta, PhD: University of Iowa
* Contact: Vincent Magnotta, vincent-magnotta@uiowa.edu
+
* Contacts: Vincent Magnotta, [mailto:vincent-magnotta@uiowa.edu vincent-magnotta@uiowa.edu]; Steven Dunn, [mailto:steven-dunn@uiowa.edu steven-dunn@uiowa.edu]
  
 
===Module Description===
 
===Module Description===
Line 36: Line 43:
  
  
<center><math>MUSH = {\color{red}a}T1 + {\color{blue}b}T2</math></center>
+
<center><math>Mush = {\color{red}a}T1 + {\color{blue}b}T2</math></center>
  
  
Line 43: Line 50:
 
the mixture model corresponding to the weighted sum image of coefficients a and b. This is very fast because each step
 
the mixture model corresponding to the weighted sum image of coefficients a and b. This is very fast because each step
 
only involves the calculation of two jointly weighted statistics. The result of optimization is the pair of linear coefficients
 
only involves the calculation of two jointly weighted statistics. The result of optimization is the pair of linear coefficients
that minimizes the sum of squares error
+
that minimizes the sum of squares error:
  
This method was applied to a sample of 20 MR brain scans and its results were compared to those obtained by 3dSkullStrip, 3dIntracranial, BET and BET2. The average Jaccard metrics for the twenty subjects was 0.66(BET), 0.61(BET2), 0.88(3dIntracranial), 0.91(3dSkullStrip) and 0.94(MUSH).
 
  
== Usage ==
+
<center><math>Error = (mean - 1000.0)^2 + (variance - 0.0)^2</math></center>
  
===Use Cases, Examples===
 
  
This module is especially appropriate for these use cases:
+
The image is thresholded by the mean signal intensity plus-or-minus five standard deviations.
  
* Use Case 1:
+
This method was applied to a sample of 20 MR brain scans and its results were compared to those obtained by 3dSkullStrip, 3dIntracranial, BET and BET2. The average Jaccard metrics for the twenty subjects was 0.66(BET), 0.61(BET2), 0.88(3dIntracranial), 0.91(3dSkullStrip) and 0.94(MUSH).
* Use Case 2:
 
  
Examples of the module in use:
+
== Usage ==
  
* Example 1
+
===Use Case===
* Example 2
+
The processing involved in the module is not designed for interactive use. It may take up to 15 minutes( or longer, depending on processor speed) to process a dataset.
  
===Tutorials===
+
===Quick Tour of Features and Use===
  
Links to tutorials explaining how to use this module:
+
A list panels in the interface, their features, what they mean, and how to use them. For instance:
  
* Tutorial 1
+
{|
** Data Set 1
+
|-
 +
|* '''Input panel:'''
 +
Specify the two input images(typically T1 and T2) as well as an optional ROI mask. This mask can be specified to aid the creation of the MUSH image, but is not necessary. Providing an ROI mask occasionally produces greater contrast in the output MUSH image(by constraining the pixels over which the mean and variance are calculated), so it is provided as an option for those who want it.
 +
|[[Image:MushInputPanel.png|thumb|280px|Input Panel]]
  
===Quick Tour of Features and Use===
+
|-
 +
|* '''Output panel:'''
 +
Specify(optionally) the weights file, which stores the final values for the MUSH equation. Specify also the filename of the MUSH image as well as the brain volume mask.
 +
|[[Image:MushOutputPanel.png|thumb|280px|Output Panel]]
  
A list panels in the interface, their features, what they mean, and how to use them. For instance:
+
|-
 +
|* '''Seed point panel:'''
 +
Specify the seed point for mask generation. This defaults to the center of the brain in a standard MRI image. Normally this shouldn't need to be changed.
 +
|[[Image:MushSeedPointPanel.png|thumb|280px|Seed Point Panel]]
  
{|
+
|-
|
+
|* '''Target statistic parameters panel:'''
* '''Input panel:'''
+
Specify a variety of advanced parameters. In most cases these will not change. Desired mean and desired variance allow the user to modify the values used in the MUSH image generation equation as listed above. The lower and upper threshold factors are used to threshold the brain mask from the MUSH image. This is done in two passes; the pre-factors specify the initial thresholding, while the remaining two factors specify the second (and primary) thresholding. While the non-pre factors could likely change, the pre-factors most likely will not. However, they have been parameterized in case it is ever needed. The bounding box size specifies an initial cubic brain mask to be used in the event an ROI mask is not provided. The bounding box start is the XYZ point-coordinate of its start.
** '''First input'''
+
|[[Image:MushTargetStatisticParametersPanel.png|thumb|280px|Target Statistic Parameters Panel]]
** '''Second input'''
 
* '''Parameters panel:'''
 
** '''First parameter'''
 
** '''Second parameter'''
 
* '''Output panel:'''
 
** '''First output'''
 
** '''Second output'''
 
* '''Viewing panel:'''
 
|[[Image:screenshotBlankNotOptional.png|thumb|280px|User Interface]]
 
 
|}
 
|}
  
 
== Development ==
 
== Development ==
 
===Notes from the Developer(s)===
 
 
Algorithms used, library classes depended upon, use cases, etc.
 
  
 
===Dependencies===
 
===Dependencies===
  
Other modules or packages that are required for this module's use.
+
[http://www.itk.org ITK]
 
 
===Tests===
 
 
 
On the [http://www.cdash.org/CDash/index.php?project=Slicer3 Dashboard], these tests verify that the module is working on various platforms:
 
 
 
* MyModuleTest1 [http://viewvc.slicer.org/viewcvs.cgi/trunk MyModuleTest1.cxx]
 
* MyModuleTest2 [http://viewvc.slicer.org/viewcvs.cgi/trunk MyModuleTest2.cxx]
 
  
 
===Known bugs===
 
===Known bugs===
  
Links to known bugs in the Slicer3 bug tracker
+
Follow this [http://www.nitrc.org/tracker/?group_id=183 link] to the bug tracker at NITRC.
  
* [http://www.na-mic.org/Bug/view.php?id=000 Bug 000: description]
 
 
 
===Usability issues===
 
 
Follow this [http://na-mic.org/Mantis/main_page.php link] to the Slicer3 bug tracker. Please select the '''usability issue category''' when browsing or contributing.
 
  
 
===Source code & documentation===
 
===Source code & documentation===
  
Links to the module's source code:
+
Available at [http://www.nitrc.org/projects/brainsmush NITRC]
 
 
Source code:
 
*[http://viewvc.slicer.org/viewcvs.cgi/trunk file.cxx ]
 
*[http://viewvc.slicer.org/viewcvs.cgi/trunk file.h ]
 
 
Doxygen documentation:
 
*[http://www.na-mic.org/Slicer/Documentation/Slicer3-doc/html/classes.html class1]
 
  
 
== More Information ==  
 
== More Information ==  
  
 
===Acknowledgment===
 
===Acknowledgment===
Include funding and other support here.
+
This work was developed by the University of Iowa Departments of Radiology and Psychiatry. This software was supported in part of NIH/NINDS award NS050568.
  
 
===References===
 
===References===
Publications related to this module go here. Links to pdfs would be useful.
 

Latest revision as of 07:31, 14 June 2013

Home < Documentation < 4.1 < Modules < BRAINSMush


For the latest Slicer documentation, visit the read-the-docs.


Return to Slicer 3.6 Documentation

Gallery of New Features


BRAINSMush

BRAINSMush

User Interface
T1 Image
T2 Image
Mush Image
Mush Brain Extraction Mask(Coronal)
Mush Brain Extraction Mask(Sagittal)

General Information

Module Type & Category

Type: CLI

Category: Segmentation

Authors, Collaborators & Contact

  • Ronald Pierson, University of Iowa
  • Gregory Harris, University of Iowa
  • Hans Johnson, PhD: University of Iowa
  • Steven A. Dunn: University of Iowa
  • Vincent Magnotta, PhD: University of Iowa
  • Contacts: Vincent Magnotta, vincent-magnotta@uiowa.edu; Steven Dunn, steven-dunn@uiowa.edu

Module Description

BRAINSMush uses the Maximize Uniformity Summation Heuristic(MUSH) optimizer, as developed at the University of Iowa, to extract the brain and surface CSF from a multi-modal imaging study. It forms a linear combination of multi-modal MR imaging data to make the signal intensity within the brain as uniform as possible. This resulting image is then thresholded to obtain the brain and surface CSF region.

MUSH uses a T1-weighted and T2-weighted image as its inputs, and their mean and variance are calculated. A linear combination is then found that approaches the desired mean and variance( by default 1000.0 and 0.0, respectively) by only varying the coefficients a and b in the following equation:


[math]Mush = {\color{red}a}T1 + {\color{blue}b}T2[/math]


Within the region of interest, the MUSH optimizer finds the number of voxels, the sum voxels and the sum of squares for both images separately. Then a 2 by 2 Levenberg-Marquardt optimizer repeatedly reconstructs the mean and variance of the mixture model corresponding to the weighted sum image of coefficients a and b. This is very fast because each step only involves the calculation of two jointly weighted statistics. The result of optimization is the pair of linear coefficients that minimizes the sum of squares error:


[math]Error = (mean - 1000.0)^2 + (variance - 0.0)^2[/math]


The image is thresholded by the mean signal intensity plus-or-minus five standard deviations.

This method was applied to a sample of 20 MR brain scans and its results were compared to those obtained by 3dSkullStrip, 3dIntracranial, BET and BET2. The average Jaccard metrics for the twenty subjects was 0.66(BET), 0.61(BET2), 0.88(3dIntracranial), 0.91(3dSkullStrip) and 0.94(MUSH).

Usage

Use Case

The processing involved in the module is not designed for interactive use. It may take up to 15 minutes( or longer, depending on processor speed) to process a dataset.

Quick Tour of Features and Use

A list panels in the interface, their features, what they mean, and how to use them. For instance:

* Input panel:

Specify the two input images(typically T1 and T2) as well as an optional ROI mask. This mask can be specified to aid the creation of the MUSH image, but is not necessary. Providing an ROI mask occasionally produces greater contrast in the output MUSH image(by constraining the pixels over which the mean and variance are calculated), so it is provided as an option for those who want it.

Input Panel
* Output panel:

Specify(optionally) the weights file, which stores the final values for the MUSH equation. Specify also the filename of the MUSH image as well as the brain volume mask.

Output Panel
* Seed point panel:

Specify the seed point for mask generation. This defaults to the center of the brain in a standard MRI image. Normally this shouldn't need to be changed.

Seed Point Panel
* Target statistic parameters panel:

Specify a variety of advanced parameters. In most cases these will not change. Desired mean and desired variance allow the user to modify the values used in the MUSH image generation equation as listed above. The lower and upper threshold factors are used to threshold the brain mask from the MUSH image. This is done in two passes; the pre-factors specify the initial thresholding, while the remaining two factors specify the second (and primary) thresholding. While the non-pre factors could likely change, the pre-factors most likely will not. However, they have been parameterized in case it is ever needed. The bounding box size specifies an initial cubic brain mask to be used in the event an ROI mask is not provided. The bounding box start is the XYZ point-coordinate of its start.

Target Statistic Parameters Panel

Development

Dependencies

ITK

Known bugs

Follow this link to the bug tracker at NITRC.


Source code & documentation

Available at NITRC

More Information

Acknowledgment

This work was developed by the University of Iowa Departments of Radiology and Psychiatry. This software was supported in part of NIH/NINDS award NS050568.

References