Image restoration using \(L_1\) norm penalty function. (English) Zbl 1258.94013

Summary: The process of estimating an original image from a given blurred and noisy image is known as image restoration. It is an ill-posed inverse problem, since one of the ways of solving it requires finding a solution to a Fredholm integral equation of convolution type in two-dimensional space. The focus of the article is to achieve a quality edge preserving image restoration using a less expensive (fast) regularization technique with \(L_1\)-norm penalty function. \(L_1\)-norm-based approaches do not penalize edges or high frequency contents in the restored image compared to \(L_2\)-norm-based approaches. Total variation (TV) is an established \(L_1\)-norm regularization approach that performs edge preserving image restoration, but at a high computational cost. TV regularization requires linearization of a highly nonlinear penalty term, which increases the restoration time considerably for large scale images. In order to reduce the computational cost, we extend least absolute shrinkage and selection operator (LASSO), an \(L_1\)-norm minimization statistical modeling technique to image restoration. The penalty function of LASSO is an identity matrix so it is computationally fast. The metrics, like, residual error, peak signal to noise ratio (PSNR), restoration time, edge map of the restored image, and subjective visual evaluation are used to assess the performances of both methods. Based on our experimental results, we show that LASSO achieves similar quality of edge preserving restoration as TV regularization, and is approximately two times faster in computation compared to TV regularization on the same set of images. We also analyze the impact of the different degree of blurring caused by point spread functions (PSFs) corrupted by different signal to noise ratios (SNRs) on image restoration.


94A08 Image processing (compression, reconstruction, etc.) in information and communication theory
Full Text: DOI


[1] Gonzalez RC, Digital Image Processing,, 2. ed. (2002)
[2] DOI: 10.1109/79.581363
[3] Nashed MZ, IEEE Transactions on Antenna Propagation 29 pp 220– (1981)
[4] Tikhonov AN, Solution of Ill-posed Problems (1977)
[5] Hadamard J, Lectures on Cauchy’s Problem in Linear Partial Differential Equations (1923) · JFM 49.0725.04
[6] Engl HW, Surveys on Mathematics for Industry 3 pp 71– (1993)
[7] DOI: 10.1007/978-94-009-1740-8
[8] DOI: 10.1007/978-3-322-99202-4
[9] Hanke M, Surveys on Mathematics for Industry 3 pp 253– (1993)
[10] DOI: 10.1137/1.9780898719697 · Zbl 0890.65037
[11] Idier J, Proceedings of the 3rd International Conference on Inverse Problems in Engineering: Theory and Practice pp 23–
[12] DOI: 10.1016/0167-2789(92)90242-F · Zbl 0780.49028
[13] Koenker R, Journal of Royal Society of Statistics, Series B 66 (1) pp 145– (2004) · Zbl 1064.62038
[14] Vogel CR, Nonsmooth regularization, Inverse Problems in Geophysical Applications (1995)
[15] Tibshirani R, Journal of Royal Statistics Society, Series B 58 (1) pp 267– (1996)
[16] DOI: 10.1137/1.9780898717570 · Zbl 1008.65103
[17] DOI: 10.1088/0266-5611/10/6/003 · Zbl 0809.35151
[18] DOI: 10.1137/0917016 · Zbl 0847.65083
[19] Chan RH, Proceedings of the International Society of Photo-optical Instrumentation Engineers pp 314–
[20] Bioucas-Das J, Proceedings of IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP)
[21] Chan TF, SIAM Journal of Scientific Computing 20 (6) pp 1964– (1999) · Zbl 0929.68118
[22] DOI: 10.1109/TPAMI.2003.1227989 · Zbl 05112222
[23] Breiman L, Better subser selection using the non-negative garotte (1993)
[24] DOI: 10.1080/00401706.1993.10485033
[25] Grandvalet , Y . 1998.Least Absolute Shrinkage is Equivalent to Quadratic Penalization, Edited by: Niklasson , L , Boden , M and Ziemske , T . 201–206. Skovde, Sweden: Springer-Verlag. Perspectives in neural computing
[26] Grandvalet Y, Neural Information Processing Systems (NIPS) 11 pp 445– (1998)
[27] Guo W, IEEE Transactions on Signal Processing Letters 6 (7) pp 165– (1999)
[28] Molina R, IEEE Transactions on Image Processing 12 (12) pp 1642– (2003) · Zbl 05453161
[29] Yang Y, IEEE Transactions on Image Processing 4 (7) pp 896– (1995)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.