Abstract
Several image restoration applications require the solution of nonnegatively constrained minimization problems whose objective function is typically constituted by the sum of a data fit function and a regularization function. Newton projection methods are very attractive because of their fast convergence, but they need an efficient implementation to avoid time consuming iterations. In this paper we present NPTool, a set of Matlab functions implementing Newton projection methods for image denoising and deblurring applications. They are specifically thought for two different data fit functions, the Least Squares function and the Kullback–Leibler divergence, and two regularization functions, Tikhonov and Total Variation, giving the opportunity of solving a large variety of restoration problems. The package is easily extensible to other linear or nonlinear data fit and regularization functions. Some examples of its use are included in the package and shown in this paper.
Similar content being viewed by others
References
Bardsley, J., Luttman, A.: Total variation-penalized Poisson likelihood estimation for ill-posed problems. Adv. Comput. Math. 31, 35–39 (2009)
Bardsley, J., Nagy, J.: Covariance-preconditioned iterative methods for nonnegatively constrained astronomical imaging. SIAM J. Matrix Anal. Appl. 27, 1184–1198 (2006)
Beck, A., Teboulle, M.: Fast gradient-based algorithms for constrained total variation image denoising and deblurring problems. IEEE Trans. Image Process. 18, 2419–2434 (2009)
Bertero, M., Boccacci, P.: Introduction to Inverse Problems in Imaging. Institute of Physics Publishing, Bristol (1998)
Bertsekas, D.: Projected Newton methods for optimization problem with simple constraints. SIAM J. Control Optim. 20(2), 221–245 (1982)
Bertsekas, D.: Nonlinear Programming, 2nd edn. Athena Scientific (1999)
Dai, Y.H., Fletcher, R.: Projected Barzilai-Borwein methods for large-scale box-constrained quadratic programming. Numer. Math. 100, 51–47 (2005)
Csiszár, I.: Why least squares and maximum entropy? An axiomatic approach to inference for linear inverse problems. Ann. Stat. 19, 2032–2066 (1991)
Dobson, D., Santosa, F.: Recovery of blocky images from noisy and blurred data. SIAM J. Appl. Math. 56, 1181–1198 (1996)
Gafni, E.M., Bertsekas, D.P.: Two-metric projection methods for constrained optimization. SIAM J. Control. Optim. 22, 936–964 (1984)
Hansen, P.C.: Rank–Deficient and Discrete Ill-Posed Problems. SIAM, Philadelphia (1997)
Hansen, P.C., Nagy, J., O’Leary, D.P.: Deblurring Images. Matrices, Spectra and Filtering. SIAM, Philadelphia (2006)
Krishnan, D., Lin, P., Yip, A.M.: A primal-dual active-set method for non-negativity constrained total variation deblurring problems. IEEE Trans. Image Process. 16, 2766–2777 (2007)
Landi, G., Loli Piccolomini, E.: A fast projected quasi-Newton method for nonnegative tikhonov regularization. Int. J. Math. Comput. Sci. 3(3), 199–213 (2008)
Landi, G., Loli Piccolomini, E.: A projected Newton-CG method for nonnegative astronomical image deblurring. Numer. Algorithms 48, 279–300 (2008)
Landi, G., Loli Piccolomini, E.: An improved Newton projection method for nonnegative deblurring of Poisson-corrupted images with Tikhonov regularization. Numer. Algorithms 60, 169–188 (2012)
Landi, G., Loli Piccolomini, E.: An efficient method for nonnegatively constrained total variation-based denoising of medical images corrupted by Poisson noise. Comput. Med. Imaging Graph. 36, 38–46 (2012)
Lucy, L.B.: An iterative technique for the rectification of observed images. Astron. J. 79, 745–754 (1974)
Nagy, J., Palmer, K., Perrone, L.: Iterative methods for image deblurring: a Matlab object-oriented approach. Numer. Algorithms 36, 73–93 (2004)
Nagy, J., Strakos, Z.: Enforcing nonnegativity in image reconstruction algorithms. Math. Model. Est. Imag. 4121, 182–190 (2000) (Wilson, D.C., et al. eds.)
Nocedal, J., Wright, S.J.: Numerical Optimization. Springer, New York (1999)
Piana, M., Bertero, M.: Projected Landweber method and preconditioning. Inverse Probl. 13, 441–463 (1997)
Rudin, L.I., Osher, S., Fatemi, E.: Nonlinear total variation based noise removal algorithms. Phys. D 60, 259–268 (1992)
Vogel, C.R.: Computational Methods for Inverse Problems. SIAM, Philadelphia (2002)
Vogel, C.R., Oman, M.E.: Fast, robust total variation–based reconstruction of noisy, blurred images. IEEE Trans. Image Process. 7, 813–824 (1998)
Wang, Z., Bovik, A.C.: Mean squared error: love it or leave it? IEEE Signal Process. Mag. 98, 98–117 (2009)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Landi, G., Piccolomini, E.L. NPTool: a Matlab software for nonnegative image restoration with Newton projection methods. Numer Algor 62, 487–504 (2013). https://doi.org/10.1007/s11075-012-9602-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-012-9602-x