Difference between revisions of "Modules:UnbiasedNonLocalMeans-Documentation-3.6"

From Slicer Wiki
Jump to: navigation, search
(Created page with 'Return to Slicer 3.6 Documentation Gallery of New Features __NOTOC__ ===Module Name=== UnbiasedNonLocalMeans == Gen…')
 
 
(11 intermediate revisions by 3 users not shown)
Line 1: Line 1:
 
[[Documentation-3.6|Return to Slicer 3.6 Documentation]]
 
[[Documentation-3.6|Return to Slicer 3.6 Documentation]]
  
[[Announcements:Slicer3.4#Highlights|Gallery of New Features]]
+
[[Announcements:Slicer3.6#Highlights|Gallery of New Features]]
 
__NOTOC__
 
__NOTOC__
 
===Module Name===
 
===Module Name===
Line 18: Line 18:
  
 
===Module Description===
 
===Module Description===
Filters a set of diffusion weighted images using Unbiased Non Local Means for Rician noise. It exploits not only the spatial redundancy, but the redundancy in similar gradient directions as well; it takes into account the N closest gradient directions to the direction being processed (a maximum of 5 gradient directions is allowed to keep a reasonable computational load, since we do not use neither similarity maps nor block-wise implementation). The noise parameter is automatically estimated. A complete description of the algorithm may be found in "DWI filtering using joint information for DTI and HARDI", by Antonio Tristan Vega and Santiago Aja-Fernandez (under review).
+
This module reduces noise (or unwanted detail) on a set of diffusion weighted images. For this, it filters the images using a Unbiased Non Local Means for Rician noise algorithm. It exploits not only the spatial redundancy, but the redundancy in similar gradient directions as well; it takes into account the N closest gradient directions to the direction being processed (a maximum of 5 gradient directions is allowed to keep a reasonable computational load, since we do not use neither similarity maps nor block-wise implementation).
 +
 
 +
The noise parameter is automatically estimated in the same way as in the jointLMMSE module.
 +
 
 +
A complete description of the algorithm may be found in:
 +
 
 +
Antonio Tristán-Vega and Santiago Aja-Fernández, "DWI filtering using joint information for DTI and HARDI", Medical Image Analysis, Volume 14, Issue 2, Pages 205-218. 2010.
 +
 
 +
Please, note that the execution of this filter is extremely slow, son only very conservative parameters (block size and search size as small as possible) should be used. Even so, its execution may take several hours. The advantage of this filter over joint LMMSE is its better preservation of edges and fine structures.
  
 
== Usage ==
 
== Usage ==
  
 
===Examples, Use Cases & Tutorials===
 
===Examples, Use Cases & Tutorials===
 +
 +
To use this filter, just load a DWI volume (gradient directions must be known) and apply the filter, as illustrated below. A brief summary of the parameters used is given below.
 +
 +
Before filtering:
 +
{|
 +
|[[Image:orig.png|thumb|600px|Original unfiltered image]]
 +
|}
 +
 +
 +
After filtering:
 +
{|
 +
|[[Image:UNLM.png|thumb|600px|Result after filtering with the unbiased non-local means method. This method performs relatively well, but has very high computational cost.]]
 +
|}
  
 
===Quick Tour of Features and Use===
 
===Quick Tour of Features and Use===
It is very easy to use it. Just select a DWI, set the parameters (if you really need it), and you're ready to go.
+
Just select a DWI, set the parameters (if you really need it), and you're ready to go.
  
 
* '''Input DWI Volume:''' set the DWI volume
 
* '''Input DWI Volume:''' set the DWI volume
Line 47: Line 68:
 
===Usability issues===
 
===Usability issues===
  
This filter is painfully slow. May take hours.
+
This filter is very slow. May take hours.
  
 +
===Source code & documentation===
 +
Source Code: Follow this [http://www.na-mic.org/ViewVC/index.cgi/trunk/Applications/CLI/DiffusionApplications/dwiUNLM link]
  
===Source code & documentation===
+
Doxygen documentation:
 +
*[http://www.na-mic.org/Slicer/Documentation/Slicer3-doc/html/classitk_1_1ComputeRestrictedHistogram.html itkComputeRestrictedHistogram.h]
 +
*[http://www.na-mic.org/Slicer/Documentation/Slicer3-doc/html/classitk_1_1OtsuStatistics.html itkOtsuStatistics.h]
 +
*[http://www.na-mic.org/Slicer/Documentation/Slicer3-doc/html/classitk_1_1OtsuThreshold.html itkOtsuThreshold.h]
 +
*[http://www.na-mic.org/Slicer/Documentation/Slicer3-doc/html/classitk_1_1UNLMFilter.html itkUNLMFilter.h]
  
 
== More Information ==  
 
== More Information ==  
  
 
===Acknowledgment===
 
===Acknowledgment===
Antonio Tristan Vega, Santiago Aja Fernandez. University of Valladolid (SPAIN). Partially founded by grant number TEC2007-67073/TCM from the Comision Interministerial de Ciencia y Tecnologia (Spain).
+
Antonio Tristan Vega, Santiago Aja Fernandez. University of Valladolid (SPAIN). Partially funded by grant number TEC2007-67073/TCM from the Ministerio de Ciencia e Innovación (Spain) and "FEDER" European Regional Development Fund.
  
 
===References===
 
===References===
 +
 +
 +
* Antonio Tristán-Vega and Santiago Aja-Fernández, "DWI filtering using joint information for DTI and HARDI", Medical Image Analysis, Volume 14, Issue 2, Pages 205-218 (2010).
 +
* Santiago Aja-Fernández, Antonio Tristán-Vega, and Carlos Alberola-López, "Noise estimation in single- and multiple-coil magnetic resonance data based on statistical models". Magnetic Resonance Imaging, Volume 27, Pages 1397-1409 (2009).

Latest revision as of 08:30, 17 May 2010

Home < Modules:UnbiasedNonLocalMeans-Documentation-3.6

Return to Slicer 3.6 Documentation

Gallery of New Features

Module Name

UnbiasedNonLocalMeans

General Information

Module Type & Category

Type: Interactive

Category: CLI/DiffusionApplications

Authors, Collaborators & Contact

  • Author: Antonio Tristán Vega , Santiago Aja Fernández
  • Contact: atriveg@bwh.harvard.edu

Module Description

This module reduces noise (or unwanted detail) on a set of diffusion weighted images. For this, it filters the images using a Unbiased Non Local Means for Rician noise algorithm. It exploits not only the spatial redundancy, but the redundancy in similar gradient directions as well; it takes into account the N closest gradient directions to the direction being processed (a maximum of 5 gradient directions is allowed to keep a reasonable computational load, since we do not use neither similarity maps nor block-wise implementation).

The noise parameter is automatically estimated in the same way as in the jointLMMSE module.

A complete description of the algorithm may be found in:

Antonio Tristán-Vega and Santiago Aja-Fernández, "DWI filtering using joint information for DTI and HARDI", Medical Image Analysis, Volume 14, Issue 2, Pages 205-218. 2010.

Please, note that the execution of this filter is extremely slow, son only very conservative parameters (block size and search size as small as possible) should be used. Even so, its execution may take several hours. The advantage of this filter over joint LMMSE is its better preservation of edges and fine structures.

Usage

Examples, Use Cases & Tutorials

To use this filter, just load a DWI volume (gradient directions must be known) and apply the filter, as illustrated below. A brief summary of the parameters used is given below.

Before filtering:

Original unfiltered image


After filtering:

Result after filtering with the unbiased non-local means method. This method performs relatively well, but has very high computational cost.

Quick Tour of Features and Use

Just select a DWI, set the parameters (if you really need it), and you're ready to go.

  • Input DWI Volume: set the DWI volume
  • Output DWI Volume: the filtered DWI volume
  • Estimation radius: This is the 3D radius of the neighborhood used for noise estimation. Noise power is estimated as the mode of the histogram of local variances
  • Search radius: This filter works by computing the weighted average of pixels in a large neighborhood. The search radius is the 3D radius of that neighborhood.
  • Comparison radius: The weights of the average are computed as the negative exponential of the (normalized) distance between two neighborhoods which have a 3D radius of this size.
  • h: The normalization constant for the distance between neighborhoods. This parameter is related to the noise power, and should be in the range 0.8->1.2. Higher values produce more blurring (and more noise reduction), while lower values better preserve edges.
  • Number of neighborhood gradients: This filter works gathering joint information from the N closest gradient directions to the one under study. This parameter is N. It is strongly recommended to keep this parameter less than 5 to keep a reasonable computational load. Even so, this filter is extremely slow.

Development

Dependencies

Volumes. Needed to load DWI volumes

Known bugs

Usability issues

This filter is very slow. May take hours.

Source code & documentation

Source Code: Follow this link

Doxygen documentation:

More Information

Acknowledgment

Antonio Tristan Vega, Santiago Aja Fernandez. University of Valladolid (SPAIN). Partially funded by grant number TEC2007-67073/TCM from the Ministerio de Ciencia e Innovación (Spain) and "FEDER" European Regional Development Fund.

References

  • Antonio Tristán-Vega and Santiago Aja-Fernández, "DWI filtering using joint information for DTI and HARDI", Medical Image Analysis, Volume 14, Issue 2, Pages 205-218 (2010).
  • Santiago Aja-Fernández, Antonio Tristán-Vega, and Carlos Alberola-López, "Noise estimation in single- and multiple-coil magnetic resonance data based on statistical models". Magnetic Resonance Imaging, Volume 27, Pages 1397-1409 (2009).