Abstract
The removal of Poisson or Poisson-Gaussian noise is often performed through the following three-step procedure. First, the noise variance is stabilized by applying either the Anscombe or the Generalized Anscombe root transformation (also called Anscombe transform) to the data, producing a signal in which the noise can be treated as additive Gaussian noise with unitary variance. Second, the noise is removed using a conventional denoising algorithm for additive white Gaussian noise. Third, an inverse transformation is applied to the denoised signal, obtaining the estimate of the signal of interest. The choice of the proper inverse transformation is crucial in order to minimize the bias error which arises when the nonlinear forward transformation is applied. We introduce the exact unbiased inverses of the Anscombe and Generalized Anscombe transformations and show that they play an integral part in ensuring accurate denoising results, particularly at the low-count regime, leading to state-of-the-art photon-limited imaging without any notable increase in the computational complexity compared to the other inverses. We also show that these inverses are optimal in the sense that they can be interpreted as maximum-likelihood inverses. Moreover, we thoroughly analyze the behaviour of the proposed inverses, which enables us to derive their closed-form approximations. Our recent work on iterative filtering of combinations of the noisy image with a progressively refined estimate leads to very efficient denoising and deblurring of even extremely low-count images (less than one count per pixel) using off-the-shelf Gaussian filters. |
SOFTWARE | PEOPLE | LINKS | REFERENCES |
Software |
Information and instructions |
Denoising software for Poisson and Poisson-Gaussian data for Matlab (ver. 7 or later) 5-Mbyte zip-file includes functions implementing the exact unbiased inverse in the iterative framework for Poisson denoising via VST v3.00, released March 16, 2016 | ||
Information and instructions |
Iterative Poisson image denoising software for Matlab (ver. 7 or later) 3.3-Mbyte zip-file Stand-alone package: includes all required invansc and BM3D components. v1.00, released March 16, 2016 | ||
Information and instructions |
Iterative Poisson image deblurring software for Matlab (ver. 7 or later) 1.5-Mbyte zip-file Stand-alone package: includes all required invansc and BM3D components. v1.00, released April 18, 2017 |
People |
Links |
References |
L. Azzari and A. Foi, “Variance Stabilization for Noisy+Estimate Combination in Iterative Poisson Denoising”, IEEE Signal Processing Letters, vol. 23, no. 8, pp. 1086-1090, August 2016. http://doi.org/10.1109/LSP.2016.2580600 Supplementary material
M. Mäkitalo and A. Foi, “Noise parameter mismatch in variance stabilization, with an application to Poisson-Gaussian noise estimation” IEEE Trans. Image Process., vol. 23, no. 12, pp. 5348-5359, December 2014. http://doi.org/10.1109/TIP.2014.2363735
M. Mäkitalo and A. Foi, “Optimal inversion of the generalized Anscombe transformation for Poisson-Gaussian noise”, IEEE Trans. Image Process., vol. 22, no. 1, pp. 91-103, January 2013. doi:10.1109/TIP.2012.2202675
M. Mäkitalo and A. Foi, “Optimal inversion of the Anscombe transformation in low-count Poisson image denoising”, IEEE Trans. Image Process., vol. 20, no. 1, pp. 99-109, January 2011. doi:10.1109/TIP.2010.2056693
M. Mäkitalo and A. Foi, “A closed-form approximation of the exact unbiased inverse of the Anscombe variance-stabilizing transformation”, IEEE Trans. Image Process., vol. 20, no. 9, pp. 2697-2698, September 2011. doi:10.1109/TIP.2011.2121085
M. Mäkitalo and A. Foi, “On the inversion of the Anscombe transformation in low-count Poisson image denoising”, Proc. Int. Workshop on Local and Non-Local Approx. in Image Process., LNLA 2009, Tuusula, Finland, pp. 26-32, August 2009. doi:10.1109/LNLA.2009.5278406