Implementation of the N3 MRI bias field correction algorithm. More...
#include <itkN3MRIBiasFieldCorrectionImageFilter.h>
Implementation of the N3 MRI bias field correction algorithm.
The nonparametric nonuniform intensity normalization (N3) algorithm is a method for correcting nonuniformity associated with MR images. The algorithm assumes a simple parametric model (Gaussian) for the bias field but does not require tissue class segmentation. In addition, there are only a couple of parameters to tune with the default values performing quite well.
N3 has been publicly available as a set of perl scripts (http://www.bic.mni.mcgill.ca/software/N3/) but, with this class, has been reimplemented for the ITK library with only one minor variation involving the b-spline fitting routine. We replaced the original fitting approach with the itkBSplineScatteredDataPointSetToImageFilter which is not susceptible to ill-conditioned matrix calculation as is the original proposed fitting component.
Notes for the user: 0. Based on our experience with the filter, sometimes the scale of the bias field is too small particularly for large mesh element sizes. Therefore, we added an option (which is true by default) for finding the optimal scaling of the bias field which is based on minimizing the coefficient of variation of the corrected image. This cost function is given below in N3BiasFieldScaleCostFunction.
Contributed by Nicholas J. Tustison, James C. Gee in the Insight Journal paper:
J.G. Sled, A.P. Zijdenbos and A.C. Evans. "A Nonparametric Method for Automatic Correction of Intensity Nonuniformity in MRI Data" IEEE Transactions on Medical Imaging, Vol 17, No 1. Feb 1998.