TY - JOUR

T1 - Tikhonov regularized Poisson likelihood estimation

T2 - Theoretical justification and a computational method

AU - Bardsley, Johnathan M.

AU - Laobeul, N'Djekornom

PY - 2008/1

Y1 - 2008/1

N2 - The noise contained in images collected by a charge coupled device camera is predominantly of Poisson type. This motivates the use of the negative logarithm of the Poisson likelihood in place of the ubiquitous least squares fit-to-data. However if the underlying mathematical model is assumed to have the form z = Au, where A is a linear, compact operator, Poisson likelihood estimation is ill-posed, and hence some form of regularization is required. In Bardsley, J.M. and Vogel, C.R., 2004, A nonnegatively constrained convex programming method for image reconstruction. SIAM Journal on Scientific Computing, 25(4), pp. 1326-1343, a numerical method is presented and analyzed for Tikhonov regularized Poisson likelihood estimation, but no theoretical justification of the approach is given. Our primary objective in this article is to provide such a theoretical justification. We then briefly present the computational method of Bardsley, J.M. and Vogel, C.R., 2004, A nonnegatively constrained convex programming method for image reconstruction. SIAM Journal on Scientific Computing, 25(4), pp. 1326-1343, which is very effective and computationally efficient for this problem. The practical validity of the approach is then demonstrated on a synthetic example from astronomical imaging.

AB - The noise contained in images collected by a charge coupled device camera is predominantly of Poisson type. This motivates the use of the negative logarithm of the Poisson likelihood in place of the ubiquitous least squares fit-to-data. However if the underlying mathematical model is assumed to have the form z = Au, where A is a linear, compact operator, Poisson likelihood estimation is ill-posed, and hence some form of regularization is required. In Bardsley, J.M. and Vogel, C.R., 2004, A nonnegatively constrained convex programming method for image reconstruction. SIAM Journal on Scientific Computing, 25(4), pp. 1326-1343, a numerical method is presented and analyzed for Tikhonov regularized Poisson likelihood estimation, but no theoretical justification of the approach is given. Our primary objective in this article is to provide such a theoretical justification. We then briefly present the computational method of Bardsley, J.M. and Vogel, C.R., 2004, A nonnegatively constrained convex programming method for image reconstruction. SIAM Journal on Scientific Computing, 25(4), pp. 1326-1343, which is very effective and computationally efficient for this problem. The practical validity of the approach is then demonstrated on a synthetic example from astronomical imaging.

KW - Ill-posed problems

KW - Image reconstruction

KW - Maximum likelihood estimation

KW - Nonnegatively constrained minimization

KW - Regularization

UR - http://www.scopus.com/inward/record.url?scp=40849092848&partnerID=8YFLogxK

U2 - 10.1080/17415970701404235

DO - 10.1080/17415970701404235

M3 - Article

AN - SCOPUS:40849092848

SN - 1741-5977

VL - 16

SP - 199

EP - 215

JO - Inverse Problems in Science and Engineering

JF - Inverse Problems in Science and Engineering

IS - 2

ER -