Variable metric techniques for forward-backward methods in imaging

Abstract Variable metric techniques are a crucial ingredient in many first order optimization algorithms. In practice, they consist in a rule for computing, at each iteration, a suitable symmetric, positive definite scaling matrix to be multiplied to the gradient vector. Besides quasi-Newton BFGS techniques, which represented the state-of-the-art since the 70’s, new approaches have been proposed in the last decade in the framework of imaging problems expressed in variational form. Such recent approaches are appealing since they can be applied to large scale problems without adding significant computational costs and they produce an impressive improvement in the practical performances of first order methods. These scaling strategies are strictly connected to the shape of the specific objective function and constraints of the optimization problem they are applied to; therefore, they are able to effectively capture the problem features. On the other side, this strict problem dependence makes difficult extending the existing techniques to more general problems. Moreover, in spite of the experimental evidence of their practical effectiveness, their theoretical properties are not well understood. The aim of this paper is to investigate these issues; in particular, we develop a unified framework for scaling techniques, multiplicative algorithms and the Majorization–Minimization approach. With this inspiration, we propose a scaling matrix rule for variable metric first order methods applied to nonnegatively constrained problems exploiting a suitable structure of the objective function. Finally, we evaluate the effectiveness of the proposed approach on some image restoration problems.

[1]  Saverio Salzo,et al.  The Variable Metric Forward-Backward Splitting Algorithm Under Mild Differentiability Assumptions , 2016, SIAM J. Optim..

[2]  H. Sebastian Seung,et al.  Learning the parts of objects by non-negative matrix factorization , 1999, Nature.

[3]  T. Chan,et al.  On the Convergence of the Lagged Diffusivity Fixed Point Method in Total Variation Image Restoration , 1999 .

[4]  Patrick L. Combettes,et al.  Proximal Splitting Methods in Signal Processing , 2009, Fixed-Point Algorithms for Inverse Problems in Science and Engineering.

[5]  I. Loris,et al.  On the convergence of variable metric line-search based proximal-gradient method under the Kurdyka-Lojasiewicz inequality , 2016 .

[6]  L. Lucy An iterative technique for the rectification of observed distributions , 1974 .

[7]  L. Zanni,et al.  New adaptive stepsize selections in gradient methods , 2008 .

[8]  L. Zanni,et al.  A scaled gradient projection method for constrained image deblurring , 2008 .

[9]  Jeffrey A. Fessler,et al.  Fast Monotonic Algorithms for Transmission Tomography , 1999, IEEE Trans. Medical Imaging.

[10]  Federica Porta,et al.  Variable Metric Inexact Line-Search-Based Methods for Nonsmooth Optimization , 2015, SIAM J. Optim..

[11]  Luca Zanni,et al.  A New Steplength Selection for Scaled Gradient Methods with Application to Image Deblurring , 2014, J. Sci. Comput..

[12]  J. M. Martínez,et al.  Gradient Method with Retards and Generalizations , 1998 .

[13]  Hedy Attouch,et al.  The Rate of Convergence of Nesterov's Accelerated Forward-Backward Method is Actually Faster Than 1/k2 , 2015, SIAM J. Optim..

[14]  Valeria Ruggiero,et al.  On the steplength selection in gradient methods for unconstrained optimization , 2018, Appl. Math. Comput..

[15]  W. Hager,et al.  The cyclic Barzilai-–Borwein method for unconstrained optimization , 2006 .

[16]  Alvaro R. De Pierro,et al.  On the convergence of an EM-type algorithm for penalized likelihood estimation in emission tomography , 1995, IEEE Trans. Medical Imaging.

[17]  M. Daube-Witherspoon,et al.  An Iterative Image Space Reconstruction Algorthm Suitable for Volume ECT , 1986, IEEE Transactions on Medical Imaging.

[18]  Juan Peypouquet,et al.  Splitting Methods with Variable Metric for Kurdyka–Łojasiewicz Functions and General Convergence Rates , 2015, J. Optim. Theory Appl..

[19]  H Lantéri,et al.  A general method to devise maximum-likelihood signal restoration multiplicative algorithms with non-negativity constraints , 2001, Signal Process..

[20]  H. Lantéri,et al.  Penalized maximum likelihood image restoration with positivity constraints:multiplicative algorithms , 2002 .

[21]  Émilie Chouzenoux,et al.  A Stochastic Majorize-Minimize Subspace Algorithm for Online Penalized Least Squares Estimation , 2015, IEEE Transactions on Signal Processing.

[22]  M. Bertero,et al.  Efficient gradient projection methods for edge-preserving removal of Poisson noise , 2009 .

[23]  P. L. Combettes,et al.  Variable metric quasi-Fejér monotonicity , 2012, 1206.5705.

[24]  J. Borwein,et al.  Two-Point Step Size Gradient Methods , 1988 .

[25]  Luca Zanni,et al.  Reconstruction of 3D X-ray CT images from reduced sampling by a scaled gradient projection algorithm , 2018, Comput. Optim. Appl..

[26]  Émilie Chouzenoux,et al.  A block coordinate variable metric forward–backward algorithm , 2016, Journal of Global Optimization.

[27]  D S Biggs,et al.  Acceleration of iterative image restoration algorithms. , 1997, Applied optics.

[28]  William H. Richardson,et al.  Bayesian-Based Iterative Method of Image Restoration , 1972 .

[29]  Michel Barlaud,et al.  Deterministic edge-preserving regularization in computed imaging , 1997, IEEE Trans. Image Process..

[30]  Valeria Ruggiero,et al.  Scaling Techniques for ε-Subgradient Methods , 2016, SIAM J. Optim..

[31]  Daniel Kolditz,et al.  Iterative reconstruction methods in X-ray CT. , 2012, Physica medica : PM : an international journal devoted to the applications of physics to medicine and biology : official journal of the Italian Association of Biomedical Physics.

[32]  S. Bonettini,et al.  New convergence results for the scaled gradient projection method , 2014, 1406.6601.

[33]  Valeria Ruggiero,et al.  Inertial Variable Metric Techniques for the Inexact Forward-Backward Algorithm , 2018, SIAM J. Sci. Comput..

[34]  Tieyong Zeng,et al.  Variational Approach for Restoring Blurred Images with Cauchy Noise , 2015, SIAM J. Imaging Sci..

[35]  Luca Zanni,et al.  Limited-memory scaled gradient projection methods for real-time image deconvolution in microscopy , 2014, Commun. Nonlinear Sci. Numer. Simul..

[36]  Zhi-Quan Luo,et al.  A Unified Convergence Analysis of Block Successive Minimization Methods for Nonsmooth Optimization , 2012, SIAM J. Optim..

[37]  Émilie Chouzenoux,et al.  Variable Metric Forward–Backward Algorithm for Minimizing the Sum of a Differentiable Function and a Convex Function , 2013, Journal of Optimization Theory and Applications.

[38]  Valeria Ruggiero,et al.  Steplength selection in gradient projection methods for box-constrained quadratic programs , 2019, Appl. Math. Comput..

[39]  William W. Hager,et al.  An efficient gradient method using the Yuan steplength , 2014, Comput. Optim. Appl..

[40]  Valeria Ruggiero,et al.  A Variable Metric Forward-Backward Method with Extrapolation , 2015, SIAM J. Sci. Comput..

[41]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[42]  José M. Bioucas-Dias,et al.  Adaptive total variation image deblurring: A majorization-minimization approach , 2009, Signal Process..

[43]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[44]  Bin Zhou,et al.  Gradient Methods with Adaptive Step-Sizes , 2006, Comput. Optim. Appl..

[45]  Erkki Oja,et al.  Unified Development of Multiplicative Algorithms for Linear and Quadratic Nonnegative Matrix Factorization , 2011, IEEE Transactions on Neural Networks.

[46]  Luca Zanni,et al.  Scaling techniques for gradient projection-type methods in astronomical image deblurring , 2013, Int. J. Comput. Math..

[47]  P. L. Combettes,et al.  Variable metric forward–backward splitting with applications to monotone inclusions in duality , 2012, 1206.6791.

[48]  Patrick L. Combettes,et al.  Signal Recovery by Proximal Forward-Backward Splitting , 2005, Multiscale Model. Simul..

[49]  Roger Fletcher,et al.  On the asymptotic behaviour of some new gradient methods , 2005, Math. Program..

[50]  Boris Polyak,et al.  Constrained minimization methods , 1966 .

[51]  A. Chambolle,et al.  On the Convergence of the Iterates of the “Fast Iterative Shrinkage/Thresholding Algorithm” , 2015, J. Optim. Theory Appl..