TY - JOUR

T1 - Covariance-preconditioned iterative methods for nonnegatively constrained astronomical imaging

AU - Bardsley, Johnathan M.

AU - Nagy, James G.

PY - 2006

Y1 - 2006

N2 - We consider the problem of solving ill-conditioned linear systems Ax = b subject to the nonnegativity constraint x ≥ 0, and in which the vector b is a realization of a random vector b, i.e., b is noisy. We explore what the statistical literature tells us about solving noisy linear systems; we discuss the effect that a substantial black background in the astronomical object being viewed has on the underlying mathematical and statistical models; and, finally, we present several covariance-based preconditioned iterative methods that incorporate this information. Each of the methods presented can be viewed as an implementation of a preconditioned modified residual-norm steepest descent algorithm with a specific preconditioner, and we show that, in fact, the well-known and often used Richardson-Lucy algorithm is one such method. Ill-conditioning can inhibit the ability to take advantage of a priori statistical knowledge, in which case a more traditional preconditioning approach may be appropriate. We briefly discuss this traditional approach as well. Examples from astronomical imaging are used to illustrate concepts and to test and compare algorithms.

AB - We consider the problem of solving ill-conditioned linear systems Ax = b subject to the nonnegativity constraint x ≥ 0, and in which the vector b is a realization of a random vector b, i.e., b is noisy. We explore what the statistical literature tells us about solving noisy linear systems; we discuss the effect that a substantial black background in the astronomical object being viewed has on the underlying mathematical and statistical models; and, finally, we present several covariance-based preconditioned iterative methods that incorporate this information. Each of the methods presented can be viewed as an implementation of a preconditioned modified residual-norm steepest descent algorithm with a specific preconditioner, and we show that, in fact, the well-known and often used Richardson-Lucy algorithm is one such method. Ill-conditioning can inhibit the ability to take advantage of a priori statistical knowledge, in which case a more traditional preconditioning approach may be appropriate. We briefly discuss this traditional approach as well. Examples from astronomical imaging are used to illustrate concepts and to test and compare algorithms.

KW - Image restoration

KW - Linear models

KW - Preconditioning

KW - Statistical methods

KW - Weighted least squares

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

U2 - 10.1137/040615043

DO - 10.1137/040615043

M3 - Article

AN - SCOPUS:33750149637

SN - 0895-4798

VL - 27

SP - 1184

EP - 1197

JO - SIAM Journal on Matrix Analysis and Applications

JF - SIAM Journal on Matrix Analysis and Applications

IS - 4

ER -