Insight Journal - N4ITK: Nick's N3 ITK Implementation For MRI Bias Field Correction

N4ITK: Nick's N3 ITK Implementation For MRI Bias Field Correction

Please use this identifier to cite or link to this publication: http://hdl.handle.net/10380/3053

Published in The Insight Journal - 2009 January - June.

Submitted by Nick Tustison on 04-30-2009.

Several algorithms exist for correcting the nonuniform intensity in magnetic resonance images caused by field inhomogeneities. These algorithms constitute important preprocessing steps for subsequent image analysis tasks. One such algorithm, known as parametric bias field correction (PABIC), has already been implemented in ITK. Another popular algorithm is the nonuniform intensity normalization (N3) approach. A particularly salient advantage of this algorithm is that it does not require a prior tissue model for its application. In addition, the source code for N3 is publicly available at the McConnell Brain Imaging Centre (Montreal Neurological Institute, McGill University) which includes source code and the coordinating set of perl scripts. This submission describes an implementation of the N3 algorithm for the Insight Toolkit given as a single class, viz. itk::N3MRIBiasFieldCorrectionImageFilter. We tried to maintain minimal difference between the publicly available MNI N3 implementation and our ITK im- plementation. The only intentional variation is the substitution of an earlier contribution, i.e. the class itk::BSplineScatteredDataPointSetToImageFilter, for the originally proposed least-squares approach for B-spline fitting used to model the bias field. In addition, we include a more extensive modification to the original N3 algorithm found in the class itk::N4MRIBiasFieldCorrectionImageFilter. The latter algorithm employs a multi-resolution approach, similar to FFD image registration strategies, and has a slightly modified iterative update scheme.

Code

Automatic Testing Results by Insight-Journal Dashboard on Fri Jan 22 22:26:54 2010 for revision #9

expertise: 5 sensitivity: 5

Automatic Testing Results by Insight-Journal Dashboard on Tue Oct 27 14:44:12 2009 for revision #6

expertise: 5 sensitivity: 5

Automatic Testing Results by Insight-Journal Dashboard on Tue Oct 20 11:18:04 2009 for revision #5

expertise: 5 sensitivity: 5

Automatic Testing Results by Insight-Journal Dashboard on Fri Jul 17 07:23:06 2009 for revision #4

expertise: 5 sensitivity: 5

Reviews

Good contribution by Dženan Zukić on 2011-05-30 09:14:31 for revision #11

expertise: 3 sensitivity: 5

Summary:

A valuable contribution, which spares us from implementing the N3 method. It has basic documentation and a very complex example, but it lacks a simplest example.

Open Science:

Source code provided, as well as some images.

Reproducibility:

It compiles after ITK was configured to use reivew code. The complex example runs noticeably faster than the code I used:

typedef itk::N4MRIBiasFieldCorrectionImageFilter<RealImageType, CharImageType> CorrecterType;

CorrecterType::Pointer correcter=CorrecterType::New();

correcter->SetInput(reader->GetOutput());

correcter->Update();

Open source Contributions:

Code is usable, but the examples provided are long, and a straightforward use of the filter takes much longer than the provided examples.

Quality of the data :

I did not check the included data, I used my own data.

Interest:

Everyone who uses MR images.

Comment by Itezziodiss Holland:

Comment by Lari7usa Vigorda:

Comment by Ldovandebrreina Holland:

Comment by Lari7usa Vigorda:

Comment by Ntlnoel Holland:

Comment by Ailiaalta Holland:

Comment by Lheeralalkris Olivier:

Comment by Tyamatriley Vigorda:

Comment by Lithrailesmered Olivier:

Comment by Elissierdeon Olivier:

Comment by Pshopgolf Olivier:

Comment by Lbbellquin Holland:

Comment by Lheeralalkris Olivier:

Comment by Tbootsbuy Vigorda:

Comment by Tsmccathedema Ken:

Comment by Rsalbarrantravi Olivier:

Comment by Bridesmaidbride Ken:

Comment by Pshopgolf Olivier:

Comment by Aollwaymyrt Vigorda:

Comment by Lheeralalkris Olivier:

Comment by Ingtroglinsterl Vigorda:

Comment by Rankeurensuan Ken:

Comment by Tbootsbuy Vigorda:

Comment by Ntveronesbrid Ken:

Comment by Mwangmajor Holland:

Comment by Uejesusrodg Holland:

Comment by Rsalbarrantravi Olivier:

Comment by Conecoyn Vigorda:

Comment by Lari7usa Vigorda:

Comment by Hlinksshopping Olivier:

Comment by :

Comment by Eveningdiscount Ken:

Comment by Eunicerosa Vigorda:

Comment by Sdubersterodo Olivier:

Comment by Aollwaymyrt Vigorda:

Comment by Dolanguiraorlan Holland:

Comment by Sdubersterodo Olivier:

Comment by Mwangmajor Holland:

Comment by Uejesusrodg Holland:

Comment by Tbootsbuy Vigorda:

Comment by Elissierdeon Olivier:

Very nice contribution for open source community by Andriy Fedorov on 2009-06-17 16:05:40 for revision #3

expertise: 2 sensitivity: 5

Summary:

The authors present an implementation of a published method for MRI bias correction. The code follows ITK style, is accompanied by guidelines on parameter initialization and references to the literature. The authors supply the code with the data they used to test the code.

Hypothesis:

n/a

Evidence:

n/a

Open Science:

The code uses ITK, the authors do provide images, and sufficient details to replicate the results.

Reproducibility:

The code is compiled with review enabled. We were able to reproduce the results authors present, and successfully apply the implementation to our data.

Use of Open Source Software:

The code promotes best practices of open source and ITK programming.

Open source Contributions:

The code is usable, the authors provide flexibility to work with either 2d or 3d images, the code includes not only filters, but also the test driver program to simplify usage.

Code Quality :

The coding style seems to be very good, the code is readable and follows ITK design. I did not observe any features that would make this code non-portable, but personally have compiled it only on linux.

Quality of the data :

The authors provide reasonable testing data sets.

Interest:

CORRECTION: the accompanying document does not stress the meaning of the shrink factor. My original experiments we using shrink factor 1, which results in no resampling. With the shrink factor 4, the code takes under 1 minute for the same data. At the same time, I cannot distinguish visually the bias field from the one produced with shrink factor 1. I withdraw my complaint about high execution time.

Free comment :

The code is rather slow on 3d data: 256x256x124 image takes about 50 min to process. The code is not multi-threaded.

Mask generation can be challenging due to the nature of input data -- field inhomogeneities may not (and do not, in our case) allow to generate the mask with simple Otsu threshold.

The algorithm requires setting of a relatively large number of parameters, and due to significant execution time, parameter tuning is not easy. We were able to achieve visually pleasing results by increasing the m_BiasFieldFullWidthAtHalfMaximum parameter from its default value.

A wrapper module for 3D Slicer (http://www.slicer.org) has been implemented by Nicolas Rannou (Brigham and Women's Hospital), and can be downloaded from NAMIC SandBox here: http://svn.na-mic.org/NAMICSandBox/trunk/MRIBiasCorrectionFilter/. This module has preview mode, which allows to apply correction to a single slice, making parameter selection easier. The use of Slicer allows to use additional filters to construct and edit mask that can be used in bias correction.

Cannot Compile the program by Sylvain Bouix on 2009-05-21 15:29:36 for revision #1

expertise: 2 sensitivity: 5

Add a new review

Quick Comments

Comment by Amir Khojaste

hello guys, I'm at my first steps to learn c++, I really appreciate your help if you could tell me how I can run this code? it has some errors, can i have the executable file of this project? thanks a lot. I would be so grateful if you also could email me at a.khojaste@queensu.ca
regards

Resources


Download All
Download Paper , View Paper
Download Source code
Github

Statistics more


Global rating:
Review rating:	[review]
Code rating:
Views:	46673
Downloads:	14768
Paper Quality:
1 comment

Information more


Paper Id:	640
Categories:	Filtering, IO
Keywords:	MRI, bias field, inhomogeneity, ninjas
Toolkits:	ITK
Revision:
Export citation:
Status:	Open for public review

Share

Linked Publications more


Diffeomorphic Demons Using ITK's Finite Difference Solver Hierarchy by Vercauteren T., Pennec X., Perchant A., Ayache N.
Reader/Writer for Analyze Object Maps for ITK by Hawley J., Johnson H.
A Label Geometry Image Filter for Multiple Object Measurement by Padfield D., Miller J.
A Skull-Stripping Filter for ITK by Bauer S., Fejes T., Reyes M.
Anisotropic Fast Marching in ITK by Mirebeau J.
Importing Contours from DICOM-RT Structure Sets with ITK4 by Dowling J.
An ITK Implementation of the Symmetric Log-Domain Diffeomorphic Demons Algorithm by Dru F., Fillard P., Vercauteren T.
Segmentation of Head and Neck CT Scans Using Atlas-based Level Set Method by Zhang X., Tian J., Wu Y., Zheng J., Deng K.
Incremental Delaunay Triangulation by Rigaud S., Gouaillard A.
An OpenCL implementation of the Gaussian pyramid and the resampler by Shamonin D., Staring M.

View license

Loading license...

Send a message to the author