Geodesic
A nonnegatively constrained trust region algorithm for the restoration of images with an unknown blur
Electronic transactions on numerical analysis, Tome 20 (2005), pp. 139-153
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 CG 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.
Classification : 65K10, 65F22
Keywords: constrained optimization, phase diversity, astronomical imaging 
@article{ETNA_2005__20__a7,
     author = {Bardsley,  Johnathan M.},
     title = {A nonnegatively constrained trust region algorithm for the restoration of images with an unknown blur},
     journal = {Electronic transactions on numerical analysis},
     pages = {139--153},
     year = {2005},
     volume = {20},
     zbl = {1121.65325},
     language = {en},
     url = {http://geodesic.mathdoc.fr/item/ETNA_2005__20__a7/}
}
TY  - JOUR
AU  - Bardsley,  Johnathan M.
TI  - A nonnegatively constrained trust region algorithm for the restoration of images with an unknown blur
JO  - Electronic transactions on numerical analysis
PY  - 2005
SP  - 139
EP  - 153
VL  - 20
UR  - http://geodesic.mathdoc.fr/item/ETNA_2005__20__a7/
LA  - en
ID  - ETNA_2005__20__a7
ER  - 
%0 Journal Article
%A Bardsley,  Johnathan M.
%T A nonnegatively constrained trust region algorithm for the restoration of images with an unknown blur
%J Electronic transactions on numerical analysis
%D 2005
%P 139-153
%V 20
%U http://geodesic.mathdoc.fr/item/ETNA_2005__20__a7/
%G en
%F ETNA_2005__20__a7
Bardsley,  Johnathan M. A nonnegatively constrained trust region algorithm for the restoration of images with an unknown blur. Electronic transactions on numerical analysis, Tome 20 (2005), pp. 139-153. http://geodesic.mathdoc.fr/item/ETNA_2005__20__a7/