A First-Order Framework for Inverse Imaging Problems

Christopher Anand(anandcmcmaster.ca)

Abstract: We argue that some inverse problems arising in imaging can be efficiently treated using only single-precision (or other reduced-precision) arithmetic, using a combination of old ideas (first-order methods, polynomial preconditioners), and new ones (bilateral filtering, total variation). Using single precision, and having structures which parallelize in the ways needed to take advantage of low-cost/high-performance multi-core/SIMD architectures, this framework is especially suited to embedded image reconstruction applications like medical imaging. We show with a simulated magnetic resonance imaging problem that this method can be numerically effective. Since the convergence/error analysis is particularly simple for pure quadratic objectives, this approach can also be used in embedded environments with fixed computation budgets, or certification requirements. Simple analysis for the quadratic case also serves as a basis for the analysis of nonlinear problems solved via a sequence of quadratic approximations. We include one example of a nonlinear, nonquadratic penalty function.

Keywords: inverse problem, image reconstruction, single-precision

Category 1: Applications -- Science and Engineering (Biomedical Applications )

Category 2: Applications -- Science and Engineering (Other )

Citation: AdvOL2007-001, McMaster University, January/2006

Download: [PDF]

Entry Submitted: 01/25/2007
Entry Accepted: 01/25/2007
Entry Last Modified: 01/25/2007

Modify/Update this entry