Variable exponent functionals in image restoration
Introduction
Image denoising is one of the fundamental problems in image processing with numerous applications. The aim of image denoising is to design methods which can selectively smooth a noisy image without losing significant features such as edges.
Variational denoising methods are widely studied numerically and theoretically in recent years. In variational framework, the denoising problem can be expressed as follows: given an original image f, it is assumed that it has been corrupted by some additive noise n. Then the problem is to recover the true image u fromLet us consider the following representative minimization problemwhere is a constant and is a scalar parameter. The first term in the energy functional of (1.1) is a regularization term and the second term is a fidelity term. As , it is the widely used Rudin–Osher–Fatemi (ROF) model proposed in 1992 [12]. The considerable advantage of the ROF model is that it can well preserve edge sharpness and location while smooth out noise. Mathematically, it is reasonable since its solution belongs to bounded variation (BV) space which allows discontinuities in functions. However, the ROF model favors solutions that are piecewise constant which often causes the staircasing effect [11], [14], [15]. The staircasing effect creates false edges which are misleading and not satisfactory in visual effects.
Choosing in (1.1) results in isotropic diffusion which solves the staircasing effect problem but it oversmoothes images such that the edges are blurred and dislocated. A fixed value of results in anisotropic diffusion between the ROF model and the isotropic smoothing. However, there is a trade-off between piecewise smooth regions reconstruction and edge preservation.
Since different values of p should have different advantages, it encourages one to combine their benefits with a variable exponent. Blomgren et al. proposed the following minimization problem in [1]where , and p is a monotonically decreasing function. This model is a variable exponent model. It chooses diffusion speed through exponent and then can reduce the staircasing effect. Since p depends on , it is hard to establish the lower semi-continuity of the energy functional. Bollt et al. proved that this problem with an or norm fidelity term has a minimizer in [2], however, nothing about the associated heat equations was discussed.
Later, Chen et al. proposed the following model in [3]where is the Gaussian kernel, are fixed parameter, and is a user-defined threshold. Mathematically, the energy minimization problem and the associated heat flow were discussed.
Inspired by the above models, we propose the modelwhere and .
Clearly in the regions with edges, since the image gradient is large, model (1.4) approximates the ROF model, so the edges will be preserved; In relatively smooth regions since image gradient is small, model (1.4) approximates isotropic smoothing, so they will be processed into piecewise smooth regions. In other regions, the diffusion is properly adjusted by the function .
The proposed model (1.4) is simpler than (1.3) in the formulation. Meanwhile, model (1.4) is more automatic than (1.3) since no user-defined threshold is needed in (1.4). Chen et al. studied problem (1.3) in BV framework [3], however, in this paper we will study problem (1.4) in the variable exponent space .
The paper is organized as follows: in Section 2 we give some important lemmas and then prove the existence and uniqueness of the solution of the minimization problem (1.4). In Section 3 we prove the existence, uniqueness and stability of the solution of the heat flow problem and discuss the long-time behavior. In Section 4 we provide our numerical algorithm and experimental results to illustrate the effectiveness of our model in image restoration. Finally, we conclude the paper in Section 5.
Section snippets
The minimization problem
Let be a bounded open set with Lipschitz boundary, . By the definition of and Gaussian convolution, we obtain . Then there exists a constant , such that . Therefore, and . Meanwhile, since , we get in the proposed model (1.4).
Variable exponent spaces. Let be a measurable function, called variable exponent on . By we denote the family of all measurable functions on . Let
The associated heat flow to problem (1.4)
Using the steepest descent method, the associated heat flow to problem (1.4) is given byFirstly, we derive another definition of weak solution of problem (3.1), (3.2), (3.3). DenoteThen (3.1) is equivalent to , where denotes the Gateaux derivative of F about u.
Suppose u be a classical solution of (3.1), (3.2), (3.3). For each ,
Numerical results
We consider dimension . Suppose the image size is . Set be the time step and be the space step. Let . Definewhere . Then the finite difference scheme of the heat flow (3.1), (3.2), (3.3) is given by
Conclusion
In this paper, we have studied a variational exponent functional to recover images based on the models (1.2), (1.3). The significant difference between our model and (1.3) is that in our model (1.4) can approximate 1 (but larger than 1) while in (1.3) will be equal to 1 in regions with large gradient. However, theoretically, the two models are discussed in different spaces. (1.3) is studied in BV space while (1.4) is studied in variable exponent Sobolev space .
The
Acknowledgements
This work is partially supported by the National Science Foundation of Shanghai (10ZR1410200), the National Science Foundation of China (10901104, 10871126) and the Research Fund for the Doctoral Program of Higher Education (200802691037).
References (16)
- et al.
Existence of solutions for p(x)-Laplacian Dirichlet Problem
Nonlinear Anal.
(2003) - et al.
Nonlinear total variation based noise removal algorithms
Physica D
(1992) - P. Blomgren, T.F. Chan, P. Mulet, C. Wong, Total variation image restoration: numerical methods and extensions, in:...
- et al.
Graduated adaptive image denoising: local compromise between total variation and isotropic diffusion
Adv. Comput. Math.
(2009) - et al.
Variable exponent, linear growth functionals in image restoration
SIAM J. Appl. Math.
(2006) - L.C. Evans, Partial Differential Equations, Graduate Studies in Mathematics, vol. 19, A.M.S., Providence, RI,...
- et al.
On the generalized Orlicz–Sobolev space
J. Gansu Edu. College
(1998) On the existence of minimizers of the variable exponent Dirichlet energy integral
Commun. Pure Appl. Anal.
(2006)
Cited by (128)
Single image noise level estimation by artificial noise
2023, Signal ProcessingAsymptotic behavior of solutions for a new general class of parabolic Kirchhoff type equation with variable exponent sources
2023, Journal of Mathematical Analysis and ApplicationsA truncated generalized Huber prior for image smoothing
2023, Applied Mathematical ModellingOn Global Solution for a Class of p(x)-Laplacian Equations with Logarithmic Nonlinearity
2024, Mediterranean Journal of Mathematics