Code covered by the BSD License  

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

Highlights from
Joint Anisotropic Wiener filter for Diffusion Weighted MRI

5.0 | 1 rating Rate this file 10 Downloads (last 30 days) File Size: 11.5 MB File ID: #43992 Version: 1.1
image thumbnail

Joint Anisotropic Wiener filter for Diffusion Weighted MRI



18 Oct 2013 (Updated )

Denoises Diffusion Weighted MRI volumes corrupted with stationary Rician noise

| Watch this File

File Information

This module reduces stationary Rician noise on Diffusion Weighted MRI data (i.e. the data sets used to estimate diffusion tensor images; note it DOES NOT apply to anatomical MRI or any other kind of MRI data. It WILL NOT work with already reconstructed diffusion tensor volumes either). Use this m-file to denoise your DW-MRI data BEFORE estimating your diffusion tensor volume or your field of ODFs.

This software filters the image in the mean squared error sense using a Rician noise model. The N closest gradient directions to the direction being processed are filtered together to improve the results: the noise-free signal at each 3-D spatial location is seen as an n-dimensional vector which has to be estimated with the LMMSE method (Wiener filter) from a set of corrupted measurements. The covariance matrix of the noise-free vector and the cross covariance between this signal and the noise have to be estimated, which is done taking into account the image formation process.

All these estimations are performed as sample estimates in a 'shaped neighborhood' defined by the weights extracted from the structural similarity of the voxels following the same idea as in the Non-Local Means filter. The NLM distances are computed in a features space comparing only the local mean value and gradients of an RGB image created from the projections of the gradient images in a set of three independent directions.

This software works with very large 4-D data sets. This means that, depending on the actual size of your data, the filter configuration (i.e. input parameters), and the power of your computer, its execution time may range from few minutes to several tens of minutes. If you cannot afford such a long processing time, we advise you to use our C++/ITK versions of this software, that may be downloaded (the open-source code and also multi-platform, pre-compiled binaries) at:

The algorihtm is described in detail in the following two references:

[1] Antonio Tristan-Vega and Santiago Aja-Fernandez, 'DWI filtering using joint information for DTI and HARDI', Medical Image Analysis, Volume 14, Issue 2, Pages 205-218. 2010;

[2] Antonio Tristan-Vega, Veronique Brion, Gonzalo Vegas-Sanchez-Ferrero, and Santiago Aja-Fernandez, 'Merging squared-magnitude approaches to DWI denoising: An adaptive Wiener filter tuned to the anatomical contents of the image', In Proceedings of IEEE EMBC 2013, pp. 507-510. Osaka (Japan). 2013,

which we ask you to cite in case you use this tool for your research.

MATLAB release MATLAB 7.12 (R2011a)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (2)
16 Aug 2014 hariharan mahalingam

great job

Comment only
04 Nov 2013 pranjul


21 Oct 2013 1.1

Minor changes to the description of the tool, the m-files have not been changed at all

Contact us