An adaptive field rule for non-parametric MRI intensity inhomogeneity estimation algorithm

A widely accepted magnetic resonance imaging (MRI) model states that the observed voxel intensity is a piecewise constant signal intensity function corresponding to the tissue spatial distribution, corrupted with multiplicative and additive noise. The multiplicative noise is assumed to be a smooth bias field, it is called intensity inhomogeneity (IIH) field. Our approach to IIH correction is based on the definition of an energy function that incorporates some smoothness constraints into the conventional classification error function of the IIH corrected image. The IIH field estimation algorithm is a gradient descent of this energy function relative to the IIH field. We call it adaptive field rule (AFR). We comment on the likeness of our approach to the self-organizing map (SOM) learning rule, on the basis of the neighboring function that controls the influence of the neighborhood on each voxel's IIH estimation. We discuss the convergence properties of the algorithm. Experimental results show that AFR compares well with state of the art algorithms. Moreover, the mean signal intensity corresponding to each class of tissue can be estimated from the image data applying the gradient descent of the proposed energy function relative to the intensity class means. We test several variations of this gradient descent approach, which embody diverse assumptions about available a priori information.