Statistical noise analysis in GRAPPA using a parametrized noncentral Chi approximation model

The characterization of the distribution of noise in the magnitude MR image is a very important problem within image processing algorithms. The Rician noise assumed in single-coil acquisitions has been the keystone for signal-to-noise ratio estimation, image filtering, or diffusion tensor estimation for years. With the advent of parallel protocols such as sensitivity encoding or Generalized Autocalibrated Partially Parallel Acquisitions that allow accelerated acquisitions, this noise model no longer holds. Since Generalized Autocalibrated Partially Parallel Acquisitions reconstructions yield the combination of the squared signals recovered at each receiving coil, noncentral Chi statistics have been previously proposed to model the distribution of noise. However, we prove in this article that this is a weak model due to several artifacts in the acquisition scheme, mainly the correlation existing between the signals obtained at each coil. Alternatively, we propose to model such correlations with a reduction in the number of degrees of freedom of the signal, which translates in an equivalent nonaccelerated system with a minor number of independent receiving coils and, consequently, a lower signal-to-noise ratio. With this model, a noncentral Chi distribution can be assumed for all pixels in the image, whose effective number of coils and effective variance of noise can be explicitly computed in a closed form from the Generalized Autocalibrated Partially Parallel Acquisitions interpolation coefficients. Extensive experiments over both synthetic and in vivo data sets have been performed to show the goodness of fit of out model.