Nonlinear total variation based noise removal algorithms. (English) Zbl 0780.49028

Summary: A constrained optimization type of numerical algorithm for removing noise from images is presented. The total variation of the image is minimized subject to constraints involving the statistics of the noise. The constraints are imposed using Lagrange multipliers. The solution is obtained using the gradient-projection method. This amounts to solving a time dependent partial differential equation on a manifold determined by the constraints. As \(t\to\infty\) the solution converges to a steady state which is the denoised image. The numerical algorithm is simple and relatively fast. The results appear to be state-of-the-art for very noisy images. The method is noninvasive, yielding sharp edges in the image. The technique could be interpreted as a first step of moving each level set of the image normal to itself with velocity equal to the curvature of the level set divided by the magnitude of the gradient of the image, and a second step which projects the image back onto the constraint set.


49N70 Differential games and control
49N75 Pursuit and evasion games
94A08 Image processing (compression, reconstruction, etc.) in information and communication theory
Full Text: DOI


[1] Frieden, B. R., Restoring with maximum likelihood and maximum entropy, J. Opt. Soc. Am., 62, 511 (1972)
[2] Phillips, D. L., A technique for the numerical solution of certain integral equations of the first kind, J. ACM, 9, 84 (1962) · Zbl 0108.29902
[3] Twomey, S., On the numerical solution of Fredholm integral equations of the first kind by the inversion of the linear system produced by quadrature, J. ACM, 10, 97 (1963) · Zbl 0125.36102
[4] Twomey, S., The application of numerical filtering to the solution of integral equations encountered in indirect sensing measurements, J. Franklin Inst., 297, 95 (1965)
[5] Hunt, B. R., The application of constrained least squares estimation to image restoration by digital computer, IEEE Trans. Comput., 22, 805 (1973)
[6] Rudin, L., Images, numerical analysis of singularities and shock filters, Caltech., C.S. Dept. Report #TR:5250:87 (1987)
[7] Osher, S.; Rudin, L. I., Feature oriented image enhancement using shock filters, SIAM J. Num. Anal., 27, 919 (1990) · Zbl 0714.65096
[8] Alvarez, L.; Lions, P. L.; Morel, J. M., Image selective smoothing and edge detection by nonlinear diffusion, SIAM J. Num. Anal., 29, 845 (1992) · Zbl 0766.65117
[9] Osher, S.; Sethian, J., Fronts propagating with curvature dependent speed: Algorithms based on a Hamilton-Jacobi formulation, J. Comput. Phys., 79, 12 (1985)
[10] Geman, D.; Reynolds, G., Constrained restoration and the recovery of discontinuities (1990), preprint
[13] Dodge, Y., Statistical data analysis based on the \(L_1\) norm and related methods (1987), North-Holland: North-Holland Amsterdam
[14] Rosen, J. G., The gradient projection method for nonlinear programming, Part II, nonlinear constraints, J. Soc. Indust. Appl. Math., 9, 514 (1961) · Zbl 0231.90048
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.