Chi-square unbiased risk estimate for denoising magnitude MR images

In this article we develop Stein-type results for unbiased estimation of the risk associated with parametric estimators of the noncentrality parameter of chi-squared random variables on two degrees of freedom. These results allow for estimator adaptivity, and thus can be used to optimize the parameters of a broad class of typical denoising functions, subject only to weak smoothness assumptions. We show how to apply these results to the problem of enhancing magnitude magnetic resonance images, which are known to be corrupted by Rician noise. As an example, we propose a transform-domain point-wise estimator based on linear expansion of thresholds. Finally, we apply this estimator to synthetic and real image data in conjunction with the undecimated Haar wavelet transform, and conclude that it is able to outperform previous wavelet-based techniques and compares favorably with a more recent approach based on non-local means.