Implementation of an Optimal First-Order Method for Strongly Convex Total Variation Regularization
Abstract
We present a practical implementation of an optimal first-order method, due to Nesterov, for large-scale total variation regularization in tomographic reconstruction, image deblurring, etc. The algorithm applies to -strongly convex objective functions with -Lipschitz continuous gradient. In the framework of Nesterov both and are assumed known -- an assumption that is seldom satisfied in practice. We propose to incorporate mechanisms to estimate locally sufficient and during the iterations. The mechanisms also allow for the application to non-strongly convex functions. We discuss the iteration complexity of several first-order methods, including the proposed algorithm, and we use a 3D tomography problem to compare the performance of these methods. The results show that for ill-conditioned problems solved to high accuracy, the proposed method significantly outperforms state-of-the-art first-order methods, as also suggested by theoretical results.
Cite
@article{arxiv.1105.3723,
title = {Implementation of an Optimal First-Order Method for Strongly Convex Total Variation Regularization},
author = {Tobias Lindstrøm Jensen and Jakob Heide Jørgensen and Per Christian Hansen and Søren Holdt Jensen},
journal= {arXiv preprint arXiv:1105.3723},
year = {2011}
}
Comments
23 pages, 4 figures