A nonnegatively constrained trust region algorithm for the restoration of images with an unknown blur

Research output: Contribution to journalArticlepeer-review

Abstract

We consider a large-scale optimization problem with nonnegativity constraints that arises in an application of phase diversity to astronomical imaging. We develop a cost function that incorporates information about the statistics of atmospheric turbulence, and we use Tikhonov regularization to induce stability. We introduce an efficient and easily implementable algorithm that intersperses gradient projection iterations with iterations from a well-known, unconstrained Newton/trust region method. Due to the large size of our problem and to the fact that our cost function is not convex, we approximately solve the trust region subproblem via the Steihaug-Toint truncated CO iteration. Iterations from the trust region algorithm are restricted to the inactive variables. We also present a highly effective preconditioner that dramatically speeds up the convergence of our algorithm. A numerical comparison using real data between our method and another standard large-scale, bound constrained optimization algorithm is presented.

Original languageEnglish
Pages (from-to)139-153
Number of pages15
JournalElectronic Transactions on Numerical Analysis
Volume20
StatePublished - 2005

Keywords

  • Astronomical imaging
  • Constrained optimization
  • Phase diversity

Fingerprint

Dive into the research topics of 'A nonnegatively constrained trust region algorithm for the restoration of images with an unknown blur'. Together they form a unique fingerprint.

Cite this