MRI Denoising Software

During last years, we proposed several denoising methods dedicated to 3D MRI. In order to make easier the use of these denoising filters, we decided to develop a Matlab-based package with GUI containing all the proposed methods. In addition, we included SPM routines for NIFTI reading and writing to propose an efficient format compatibility with standard software. Finally, we enable the processing of several files in 2-clicks to make easier the processing of large database. The different filters cover a large range of clinical situations and you should be able to find a way to denoise your 3D MRI images with this package (MRI, fMRI, ASL...). A package dedicated to 4D images (DWI, r-fMRI...) can be found here.

Matlab package:

The Matlab package MRIdenoisingPackage_r01 contains five denoising filters for 3D MRI compiled for Linux (64 bits), Windows (64 bits) and MAC (64bits) .

All the methods are copyrighted and are the property of Pierrick Coupé and Jose Manjon.

Installation instructions:

    1. Download the package

    2. Extract the files from .zip

    3. Open Matlab

    4. Select the directory of the MRIdenoisingPackage_r01 as current directory

    5. Run the script MainMRIDenoising.m (by opening it and pressing ‘F5′ or by writing ‘MainMRIDenoising’ in the prompt).

    6. Follow the instructions of the GUI

      1. Select one or several 3D NIFTI image(s)

      2. (you can try the image provided as example in the directory Data in order to verify your architecture compatibility)

      1. Select the desired filter and tune the parameters according to your application

      1. Press run

    1. Wait and check the result on the figure. The denoising image(s) should be saved as .nii in the same directory than the noisy image(s) with the selected suffix (e.g., _denoised). One figure will be display by selected file.


Filters description and usage:

In this section, I tried to coarsely define the underlying assumptions of each method and the possible applications for fully-automatic processing. It might occur that your problem is non listed in the proposed applications. In such situation, the smoothing parameter can be tune to be adjusted for your application.

ONLM filter

Description: The first denoising filter proposed in the package is the Optimized Nonlocal Means (ONLM) [1]. This filter is currently the state-of-the-art for 3D MRI denoising and has been well-validated. The impact of this filter on segmentation or cortical surface extraction has been investigated (see here for details). A Rician-adapted version (ORNLM) has been proposed in [3] to handle the intensity bias introduced by Rician noise. This filter requires the estimation of the noise variance to be fully automatic. In the package, we use our RMAD method [5] in order to estimate the noise variance. The ONLM filter is included in several software such as MedINRIA and Minctool.

Utilization: This filter is theoretically dedicated to MRI acquired with

    • 1 coil (stationary noise)

    • without partial Fourier (non correlated noise)

In this case, the Rician option should be activated. This filter obtained good results in other clinical setups. For ASL or fMRI (separated 3D images), the option Rician should be non-activated since the difference of two Rician distributions tends to a Gaussian one. In case of zero-padding in k-space or partial Fourier acquisition, the noise variance case be under-estimated, we suggest to use increase smoothing parameter or to use AONLM which may be more robust to correlated noise. In case of SENSE or GRAPPA acquisition, we suggest to use AONLM.


AONLM filter

Description: The second denoising filter proposed in the package is the Adaptive Optimized Nonlocal Means (AONLM) [4]. This filter has been designed for spatially varying noise typically presents in parallel imaging. By using integrated noise estimation, this filter is fully automatic and quiet robust. This is obtained at the expense of longer computational time. The AONLM is included in VBM8 (SPM toolbox).

Utilization: This filter is theoretically dedicated to MRI acquired with

    • multi-coils (non stationary noise)

    • SENSE reconstruction (SENSE results in Rican noise, GRAPPA results in non-central Chi noise)

    • without partial Fourier (non correlated noise)

In this case, the Rician option should be activated. This filter obtained good results in many clinical setups especially on SENSE without partial Fourier. In practice, AONLM will be more robust than ONLM face of partial Fourier. As ONLM, AONLM can be use for ASL or fMRI (3D images separated) using Gaussian model. Finally, for GRAPPA reconstruction no adapted methods are proposed in this package. However, by using Gaussian noise model (i.e., non-activation of Rician option) good results can be obtained. Only the intensity bias correction will not be achieved.


Multi-resolution ONLM filter

Description: The third denoising filter proposed in the package is the Multi-resolution Optimized Nonlocal Means (ONLM) [1]. Based on a similar approach than ONLM, this filter is more efficient on very low SNR image compared to ONLM. A Rician-adapted version has been proposed in order to handle the intensity bias introduced by Rician noise. This filter requires the estimation of the noise variance to be fully automatic. In the package, we use our RMAD method [5] in order to estimate the noise variance.

Utilization: This filter is theoretically dedicated to MRI acquired with

    • 1 coil (stationary noise)

    • without partial Fourier (non correlated noise).

In this case, the Rician option should be activated. This filter obtained very good results in other clinical setups. For ASL or fMRI (separated 3D images), the option Rician should be non-activated since the difference of two Rician distributions tends to a Gaussian one. In case of zero-padding in k-space or partial Fourier acquisition, the noise variance case be under-estimated, we suggest to use increase smoothing parameter or to use AONLM which may be more robust to correlated noise. In case of SENSE or GRAPPA acquisition, we suggest to use AONLM.


ODCT filter

Description: The fourth denoising filter proposed in the package is the Oracle-based DCT filter [6]. This filter obtained really good results and is the fastest proposed methods. A Rician-adapted version has been proposed in order to handle the intensity bias introduced by Rician noise. This filter requires the estimation of the noise variance to be fully automatic. In the package, we use our RMAD method [5] in order to estimate the noise variance.

Utilization: This filter is theoretically dedicated to MRI acquired with

    • 1 coil (stationary noise)

    • without partial Fourier (non correlated noise)

In this case, the Rician option should be activated. This filter obtained very good results in other clinical setups. For ASL or fMRI (separated 3D images), the option Rician should be non-activated since the difference of two Rician distributions tends to a Gaussian one. In case of zero-padding in k-space or partial Fourier acquisition, the noise variance case be under-estimated, we suggest to use increase smoothing parameter or to use AONLM which may be more robust to correlated noise. In case of SENSE or GRAPPA acquisition, we suggest to use AONLM.


PRINLM filter

Description: The last denoising filter proposed in the package is the Prefiltered Rotationally Invariant NonLocal Means filter [6]. This approach used ODCT results as preprocessing step. This filter obtained the best results over the proposed methods and is the second fastest. A Rician-adapted version has been proposed in order to handle the intensity bias introduced by Rician noise. This filter requires the estimation of the noise variance to be fully automatic. In the package, we use our RMAD method [5] in order to estimate the noise variance.

Utilization: This filter is theoretically dedicated to MRI acquired with

    • 1 coil (stationary noise)

    • without partial Fourier (non correlated noise)

In this case, the Rician option should be activated. This filter obtained very good results in other clinical setups. For ASL or fMRI (separated 3D images), the option Rician should be non-activated since the difference of two Rician distributions tends to a Gaussian one. In case of zero-padding in k-space or partial Fourier acquisition, the noise variance case be under-estimated, we suggest to use increase smoothing parameter or to use AONLM which may be more robust to correlated noise. In case of SENSE or GRAPPA acquisition, we suggest to use AONLM.

Acknowledgments:

We thank the Pr. D. Louis Collins for providing the example MRI data included in this package. This MRI dataset remains the property of the Image Processing Lab. This package benefits from the use of nifti routines from the well-known SPM software (http://www.fil.ion.ucl.ac.uk/spm/). All these routines remain the property of the Wellcome Trust Centre for Neuroimaging. In addition, this package benefits from the use of 3D wavelet Matlab toolbox of the Polytechnic University of Brooklyn (http://eeweb.poly.edu/iselesni/WaveletSoftware/index.html). This Matlab toolbox remains the property of authors.


References:

    1. P. Coupé, P. Yger, S. Prima, P. Hellier, C. Kervrann, C. Barillot. An Optimized Blockwise NonLocal Means Denoising Filter for 3-D Magnetic Resonance Images. IEEE Transactions on Medical Imaging, 27(4):425–441, 2008.

    2. P. Coupé, P. Hellier, S. Prima, C. Kervrann, C. Barillot. 3D Wavelet Subbands Mixing for Image Denoising. International Journal of Biomedical Imaging, 2008.

    3. N. Wiest-Daesslé, S. Prima, P. Coupé, S.P. Morrissey, C. Barillot. Rician noise removal by non-Local Means filtering for low signal-to-noise ratio MRI: applications to DT-MRI. In MICCAI’08, pages 171–179, 2008.

    4. J. V. Manjon, P. Coupé, L. Martí-Bonmatí, D. L. Collins, M. Robles. Adaptive Non-Local Means Denoising of MR Images with Spatially Varying Noise Levels. Journal of Magnetic Resonance Imaging, 31(1):192–203, 2010.

    5. P. Coupé, J. V. Manjon, E. Gedamu, D. Arnold, M. Robles, D. L. Collins. Robust Rician Noise Estimation for MR Images. Medical Image Analysis, 14(4) : 483–493, 2010.

    6. J. V. Manjon, P. Coupé, A. Buades, D. L. Collins, M. Robles. New Methods for MRI Denoising based on Sparseness and Self-Similarity. Medical Image Analysis, 16(1): 18-27, 2012.

    7. P. Coupé, J. V Manjon, M. Robles, D. L. Collins. Adaptive multiresolution non-local means filter for three-dimensional magnetic resonance image denoising. IET Image Processing, 6(5): 558-568, 2012.