Research Article  Open Access
A Class of FourthOrder TelegraphDiffusion Equations for Image Restoration
Abstract
A class of nonlinear fourthorder telegraphdiffusion equations (TDE) for image restoration are proposed based on fourthorder TDE and bilateral filtering. The proposed model enjoys the benefits of both fourthorder TDE and bilateral filtering, which is not only edge preserving and robust to noise but also avoids the staircase effects. The existence, uniqueness, and stability of the solution for our model are proved. Experiment results show the effectiveness of the proposed model and demonstrate its superiority to the existing models.
1. Introduction
The use of secondorder partial differential equations (PDE) has been studied as a useful tool for image restoration (noise removal). They include anisotropic diffusion equations [1–3] and total variation models [4] as well as curve evolution equations [5]. These secondorder techniques have been proved to be effective for removing noise without causing excessive smoothing of the edges. However, the images resulting from these techniques are often piecewise constant, and therefore, the processed image will look “blocky” [6–8]. To reduce the blocky effect, while preserving sharp jump discontinuities, many methods in the literature [9–21] have been proposed to solve the problem.
A rather detailed analysis of blocky effects associated with anisotropic diffusion which was carried out in [6]. Let denote the image intensity function, the time, the anisotropic diffusion equation as formulated by Perona and Malik [1] presented as (the PM model) where is the conductance coefficient, and denote the divergence and the gradient, respectively. Equation (1.1) is associated with the following energy functional: where is the image support, and is an increasing function associated with the diffusion coefficient as The PM model is then associated with an energydissipating process that seeks the minimum of the energy functional. From energy functional (1.2), it is obvious that level images are the global minima of the energy functional. Detailed analysis in [6] indicates that when there is no backward diffusion, a level image is the only minimum of the energy functional, so PM model will evolve toward the formation of a level image function. Since PM model is designed such that smooth areas are diffused faster than leee smooth ones, blocky effects will appear in the early stage of diffusion even though all the blocks will finally merge to form a level image. Similarly, when there is backward diffusion, however, any piecewise level image is a global minimum of the energy functional, so blocks will appear in the early stage of the diffusion and will.
In 2000, You and Kaveh [12] proposed the following fourthorder PDE (the YK model) where denotes Laplace operator. The YK model replaces the gradient operator in PM model with a Laplace operator. Due to the fact that the Laplace of an image at a pixel is zero only if the image is planar in its neighborhood, the YK fourthorder PDE attempts to remove noise and preserve edges by approximating an observed image with a piecewise planar image. It is well known that piecewise smooth images look more natural than the piecewise constant images. Therefore, the blocky effect will be reduced and the image will look more nature. However, fourthorder PDE has an inherent limitation. The proposed PDE tends to leave images with speckle artifacts.
Recently, Ratner and Zeevi [22] introduced the following telegraphdiffusion equation (TDE model) where is the damping coefficient and is the elasticity coefficient. Note that the TDE model is derived from (1.1) by adding second time derivative of the image. It is interesting to note that (1.5) converges to the diffusion equation (1.1) not only for large and , but also after very long time [5, 6]. While one may find when the noise is large, (1.5) will be unstable in presence of noise which is similar to that of the PM model [1]. In order to overcome this drawback, most recently, Cao et al. [23] proposed the following improved telegraphdiffusion equation (ITDE model) where is a Gaussian filtering. The TDE model (1.5), together with its improved version (1.6), has certain drawbacks. Since the model uses anisotropic diffusion, the filter inherits the staircase effects of anisotropic diffusion. The Gaussian filtering is not a perfect filtering method. Moreover, since the model uses anisotropic diffusion, the filter inherits the block effects of anisotropic diffusion.
Inspired by the ideas of [12, 22], it is natural to investigate a model inherits the advantages of the YK model and TDE model. Our proposed fourthorder telegraphdiffusion equation as follows: where denotes Laplace operator. To our best knowledge, secondorder derivative is sensitive to the noise. While one may find when the image is very noisy, (1.7) will be unstable which could not distinguish correctly the “true” edges and “false” edges. Considering that we can eliminate some noise before solving model (1.6), we reformulate model (1.6) as follows: where is a bilateral filtering [24, 25]; namely, with the normalization constant where will be a spatial Gaussian that decreases the influence of distant pixels, while will be a range Gaussian that decreases the influence of pixels with intensity values that are very different from those of ; for example, where parameters and dictate the amount of filtering applied in the domain and the range of the image, respectively.
The ability of edge preservation in the fourthorder TDEbased image restoration method strongly depends on the conductance coefficient . The desirable conductance coefficient should be diffused more in smooth areas and less around less intensity transitions, so that small variations in image intensity such as noise and unwanted texture are smoothed and edges are preserved. In the implementation of our proposed method, we use the following function: where is a threshold constant.
We will recall that the main purpose of the function is to provide adaptive smoothing. It should not only precisely locate the position of the main edges, but it should also inhibit diffusion at edges. This is exactly what bilateral filtering accomplishes. It is proposed as a tool to reduce noise and preserve edges, by means of exploiting all relevant neighborhoods. It combines gray levels not only based on their gray similarity but also their geometric closeness and prefers near values to distant values in both domain and range.
The remainder of this paper is organized as follows. In Section 2, we will show the existence and uniqueness of our proposed model. Section 3 presents a discredited numerical implementation of the proposed model. Numerical experiments are presented in Section 4, and the paper is concluded in Section 5.
2. Analysis of Our Proposed Model: Existence and Uniqueness of Weak Solutions
In this section, we establish the existence and uniqueness of the following problem: where is a bounded domain of with an appropriately smooth boundary, denotes the unit outer normal to , and . In this section, will represent a generic constant that may change from line to line even if in the same inequality.
The following standard notations are used throughout. We denote is a Hilbert space for the norm where is a positive integer and of order denotes the distributional derivative of .
We denote by the set of all functions such that for almost every in , belongs to . is a normed space for the norm where is a positive integer. is a normed space for the norm We denote by the dual of and introduce the following function space Obviously, is a Banach space with the norm .
Definition 2.1. Let be a fixed positive number. A function is a weak solution of the problem (2.1)–(2.3) provided (i), and ,(ii) for any and a. e. time , (iii).
2.1. FourthOrder Linear Equation: Existence and Uniqueness
Now, we consider the existence and uniqueness of weak solutions of the following linear TD problem where .
Definition 2.2. A function is called a weak solution of the problem (2.8)–(2.10) provided (i), and ,(ii)for any and a. e. , (iii).
Theorem 2.3. Suppose that and , then the problem (2.8)–(2.10) admits a unique weak solution .
In order to prove Theorem 2.3, we will prove the following two lemmas.
Lemma 2.4. Suppose that , then there exists a constant , depending only on , and , such that Moreover, there exists a constant , depending only on , , such that for each .
Proof. By the definition of in (1.8), we then calculate Moreover, since and are smooth, we have Therefore, we draw the conclusion.
Lemma 2.5. If there exists such that Then, one has the estimate
Proof. Multiplying (2.8) by and integrating over yields
for a.e. .
Note that
Since and satisfy (2.15), then , belong to and there exists a constant depending on and such that
for any , a.e. . Therefore, since is decreasing and is smooth, there exist two constants , such that
for any , a.e. .
Consequently, (2.18) yields
In terms of (2.17) and (2.21), it then follows that
for a.e. .
Now, write
Then, (2.22) reads
for a.e. . Applying the Gronwall inequality to (2.24) yields the estimate
Since
we obtain from (2.23)–(2.25) the estimate
Multiplying (2.8) by and integrating over , we deduce for a.e. that
where and .
Consequently,
Therefore,
This completes the proof of the lemma.
Remark 2.6. By Lemma 2.5, Theorem 2.3 can be proved by the Galerkin method (see [22]). Here, we omit the details of the proof of Theorem 2.3.
2.2. FourthOrder Nonlinear Equation: Existence and Uniqueness
Theorem 2.7 (see ([26], Schauder’s Fixed Point Theorem)). Suppose that denotes a real Banach space and is compact and convex, and assume also that is continuous. Then, has a fixed point in .
Theorem 2.8. Suppose that , then the problem (2.1)–(2.3) admits one and only one weak solution.
Proof. In the following, we prove the theorem in two parts.(1) Existence of a Solution
In this first section, we show the existence of a weak solution of (2.1)–(2.3) by Shauder’s fixed point theorem [22] and Theorem 2.3. We introduce the space
Obviously, is a Banach space with the norm
Recalling Lemma 2.5, we introduce the subspace of defined by
By construction, is a mapping from into . Since is compactly imbedded in , is a nonempty, convex, and weakly compact subset of .
In order to apply the Schauder fixed point theorem, we need to prove that is a continuous compact mapping from into . Let be a sequence in which converges weakly to some in and . We have to prove that converges weakly to . By using the classical theorem of compact inclusion in Sobolev spaces [27], the sequence contains a subsequence (still denoted by) and contains a subsequence (still denoted by) such that
Since is the weak solutions of (2.8)–(2.10) corresponding to , we have, for a.e.
Integrating (2.43) from 0 to yields
From (2.38)(2.39), we have
Considering the third term in (2.44)
According to Lemma 2.4 and (2.40), we deduce, as
and note also, as , .
Letting in (2.44) yields
Therefore, by the arbitrariness of , we obtain, for a.e.
which implies
By the uniqueness of the solution of linear problem (2.8)–(2.10) and (2.37), the sequence converges weakly in to . Hence, the mapping is weakly continuous from into . This in turn shows that mapping is compact. A similar argument shows that is a continuous mapping. Applying the Schauder fixed point theorem, we conclude that has a fixed point , which consequently solve (2.1)–(2.3). Using the classical theory of parabolic equations and the bootstrap argument [26], we can deduce that is a strong solution of (1.1)–(1.4) and .(2) Uniqueness of the Solution
Now, we turn to the proof of the uniqueness, following the idea in [22]. Let and be two solutions of (2.1)–(2.3), and we have, for almost every in and
in the distribution sense, where
It suffices to show that . To verify this, fix and set
for . Then, for each . Multiplying both sides of (2.52) by , and integrating over yields
Now, (), and so
As seen in the proof of Lemma 2.5, there exists positive constants and are positive constant depending on and such that
which imply from (2.58),
By using the Young inequality to (2.60), we obtain
Now, let us write
whereupon (2.61) can be rewritten as
Note that
Therefore, (2.63) implies
Choose so small such that
Then, if , we have
Consequently, the integral form of the Gronwall inequality implies on .
Finally, we apply the same argument on the intervals , , and so forth and eventually deduce that on . Thus, we obtain the uniqueness of weak solutions.
3. Discretised Numerical Scheme
In this section, we construct an explicit discrete scheme to numerically solve differential equation (2.1)–(2.3). Assume a space grid size of and a time step size of , we quantize the space and time coordinates as follows:
where is the size of the image. Let be the approximation to the value . We then use a fivestep approach to calculate (2.1)–(2.3).
(a) calculating the Laplace of the image intensity functions with symmetric boundary conditions
calculating the value of the following function
calculating the Laplace of as with symmetry boundary conditions
defining the discrete approximation: with
finally, the numerical approximation to the differential equation (2.1) is given as
4. Experimental Results
In this section, we present numerical results obtained by applying our proposed fourthorder TDE to image denosing. We test the proposed method on “Barbara” image with size (taken from USCSIPI image database) and “license plate” image with size . These two images are shown in Figure 1(a) and Figure 3(a), respectively. The value chosen for the time step size is 0.25. To quantify the achieved performance improvements, we adopt improvement in signal to noise ratio (ISNR), which is defined as where is the initial image (noised image) and is the denoised image. The value of ISNR is large, and the restored image is better.
(a)
(b)
(c)
(d)
(e)
(f)
We first study the effects of damping coefficient . Figure 1(a) shows the noisy “Barbara” image. Figures 1(c)–1(f) show the restored image using fourthorder TDE (, ) with , 40, 70 and 100, respectively. In Figure 2, we plot the ISNR with different damping coefficient . We can see that the INSR reaches a maximum at . The INSR value using fourthorder TDE is lower than INSR value using proposed domainbased fourthorder PDE at . When , the ISNR with fourthorder TDE is almost equal to the ISNR with domainbased fourthorder PDE.
(a)
(b)
(c)
(d)
(e)
(f)
Next, we test the proposed method for image restoration on 50 synthetic degraded images generated using random white Gaussian noise of variance . To verify the effectiveness of our proposed fourthorder TDE method for image restoration, it was evaluated in comparison with PM model [1], TDE model [22], and ITDE model [23]. Figure 3(a) shows the original “License plate” image with size . We then added random white Gaussian noise of variance to generate 50 degraded image. One such degraded image is shown in Figure 3(b). The results yield by PM secondorder PDE is shown in Figure 3(c). We observe that “PM secondorder PDE” can cause the processed image look block. The results yield by TDE and ITDE are shown in Figures 3(d) and 3(e), respectively. There two methods can leave much more sharp edges than PM secondorder PDE but inherits the blocky effects from PM secondorder PDE to some degree. Figure 3(f) is the restored image using our proposed fourthorder TDE. For the 50 random generated images, the mean of ISNR values for PM, TDE, ITDE, and our proposed fourthorder TDE are 5.7123, 5.8594, 5.9041, and 6.1688, respectively. Table 1 gives 10 of the 50 ISNR values. From Figure 3 and Table 1, we can conclude that our proposed method performs the best quality.

Finally, we designed to further evaluate the good behavior of our proposed fourthorder TDE with white Gaussian noise across 5 noise levels. We added Gaussian white noise across five different variances to the original image. The ISNR values are given in Table 2. Our method obtains higher ISNR than the original method. Furthermore, a ISNR analysis conducted on the standard test image taken from USCSIPI image database (http://sipi.usc.edu/database/) is listed in Table 3. We also note that our method obtains higher mean of ISNR than the original methods.


5. Conclusions
A class of nonlinear fourthorder telegraphdiffusion equation (TDE) for image restoration is presented in this paper. The proposed model first extends the second order TDE for image restoration to fourthorder TDE. Moreover, our proposed model combines nonlinear fourthorder TDE with bilateral filtering, which is not only edge preserving and robust to noise but also avoids the staircase effects. Finally, we study the existence, uniqueness, and stability of the proposed model. A set of numerical experiments is presented to show the good performance of our proposed model. Numerical results indicate that the proposed model recovers well edges and reduces noise. In [28, 29], the authors point out that the main disadvantage of higher order methods is the complexity of computation. In the future, we will explore a fast and efficient algorithm for our proposed fourthorder TDE method.
Acknowledgments
This work was supported by National Natural Science Foundation of China under Grant no. 60972001 and National Key Technologies R & D Program of China under Grant no. 2009BAG13A06. The authors would like to acknowledge Professor Qilin Liu and Yuxiang Li from Department of Mathematics for many fruitful discussions. The authors also thank the anonymous reviewer for his or her constructive and valuable comments, which helped in improving the presentation of our work.
References
 P. Perona and J. Malik, “Scalespace and edge detection using anisotropic diffusion,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 12, no. 7, pp. 629–639, 1990. View at: Publisher Site  Google Scholar
 F. Catté, P.L. Lions, J.M. Morel, and T. Coll, “Image selective smoothing and edge detection by nonlinear diffusion,” SIAM Journal on Numerical Analysis, vol. 29, no. 1, pp. 182–193, 1992. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 G. W. Wei, “Generalized PeronaMalik equation for image restoration,” IEEE Signal Processing Letters, vol. 6, no. 7, pp. 165–167, 1999. View at: Publisher Site  Google Scholar
 L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D, vol. 60, no. 1–4, pp. 259–268, 1992. View at: Google Scholar
 B. B. Kimia, A. Tannenbaum, and S. W. Zucker, “On the evolution of curves via a function of curvature. I. The classical case,” Journal of Mathematical Analysis and Applications, vol. 163, no. 2, pp. 438–458, 1992. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. L. You, W. Xu, A. Tannenbaum, and M. Kaveh, “Behavioral analysis of anisotropic diffusion in image processing,” IEEE Transactions on Image Processing, vol. 5, no. 11, pp. 1539–1553, 1996. View at: Google Scholar
 K. Joo and S. Kim, “PDEbased image restoration, I: antistaircasing and antidiffusion,” Tech. Rep. 200307, Department of Mathematics, University of Kentucky, 2003. View at: Google Scholar
 J. Savage and K. Chen, “On multigrids for solving a class of improved total variation based staircasing reduction models,” in Image Processing Based on Partial Differential Equations: Mathematics and Visualization, Math. Vis., pp. 69–94, Springer, Berlin, Germany, 2007. View at: Publisher Site  Google Scholar
 E. Bae, J. Shi, and X.C. Tai, “Graph cuts for curvature based image denoising,” IEEE Transactions on Image Processing, vol. 20, no. 5, pp. 1199–1210, 2011. View at: Publisher Site  Google Scholar
 Q. Chen, P. Montesinos, Q. S. Sun, and D. Shen Xia, “Ramp preserving PeronaMalik model,” Signal Processing, vol. 90, no. 6, pp. 1963–1975, 2010. View at: Publisher Site  Google Scholar
 Y. Shih, C. Rei, and H. Wang, “A novel PDE based image restoration: convectiondiffusion equation for image denoising,” Journal of Computational and Applied Mathematics, vol. 231, no. 2, pp. 771–779, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y.L. You and M. Kaveh, “Fourthorder partial differential equations for noise removal,” IEEE Transactions on Image Processing, vol. 9, no. 10, pp. 1723–1730, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Lysaker, A. Lundervold, and X. C. Tai, “Noise removal using fourthorder partial differential equation with applications to medical magnetic resonance images in space and time,” IEEE Transactions on Image Processing, vol. 12, no. 12, pp. 1579–1589, 2003. View at: Publisher Site  Google Scholar
 Q. Liu, Z. Yao, and Y. Ke, “Entropy solutions for a fourthorder nonlinear degenerate problem for noise removal,” Nonlinear Analysis: Theory, Methods & Applications, vol. 67, no. 6, pp. 1908–1918, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. B. Greer and A. L. Bertozzi, “H^{1} solutions of a class of fourth order nonlinear equations for image processing,” Discrete and Continuous Dynamical Systems. Series A, vol. 10, no. 12, pp. 349–366, 2004. View at: Google Scholar  Zentralblatt MATH
 J. B. Greer and A. L. Bertozzi, “Traveling wave solutions of fourth order PDEs for image processing,” SIAM Journal on Mathematical Analysis, vol. 36, no. 1, pp. 38–68, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Q. Liu, Z. Yao, and Y. Ke, “Solutions of fourthorder partial differential equations in a noise removal model,” Electronic Journal of Differential Equations, no. 120, pp. 1–11, 2007. View at: Google Scholar  Zentralblatt MATH
 S. Didas, J. Weickert, and B. Burgeth, “Properties of higher order nonlinear diffusion filtering,” Journal of Mathematical Imaging and Vision, vol. 35, no. 3, pp. 208–226, 2009. View at: Publisher Site  Google Scholar
 M. Lysaker and X. C. Tai, “Iterative image restoration combining total variation minimization and a secondorder functional,” International Journal of Computer Vision, vol. 66, no. 1, pp. 5–18, 2006. View at: Publisher Site  Google Scholar
 F. Li, C. Shen, J. Fan, and C. Shen, “Image restoration combining a total variational filter and a fourthorder filter,” Journal of Visual Communication and Image Representation, vol. 18, no. 4, pp. 322–330, 2007. View at: Publisher Site  Google Scholar
 D. Yi and S. Lee, “Fourthorder partial differential equations for image enhancement,” Applied Mathematics and Computation, vol. 175, no. 1, pp. 430–440, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 V. Ratner and Y. Y. Zeevi, “Image enhancement using elastic manifolds,” in the 14th Edition of the International Conference on Image Analysis and Processing (ICIAP '07), pp. 769–774, September 2007. View at: Publisher Site  Google Scholar
 Y. Cao, J. Yin, Q. Liu, and M. Li, “A class of nonlinear parabolichyperbolic equations applied to image restoration,” Nonlinear Analysis: Real World Applications, vol. 11, no. 1, pp. 253–261, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 C. Tomasi and R. Manduchi, “Bilateral filtering for gray and color images,” in Proceedings of IEEE International Conference on Computer Vision, pp. 839–846, Bombay, India, January 1998. View at: Google Scholar
 M. Zhang and B. K. Gunturk, “Multiresolution bilateral filtering for image denoising,” IEEE Transactions on Image Processing, vol. 17, no. 12, pp. 2324–2333, 2008. View at: Publisher Site  Google Scholar
 L. C. Evans, Partial Differential Equations, vol. 19 of Graduate Studies in Mathematics, American Mathematical Society, Providence, RI, USA, 1998.
 R. A. Adams, Sobolev Spaces, vol. 6 of Pure and Applied Mathematics, Academic Press, London, UK, 1975.
 N. Komodakis and N. Paragios, “Beyond pairwise energies: efficient optimization for higherorder mrfs,” in IEEE Computer Society Conference on Computer Vision and Pattern Recognition Workshops (CVPR '09), pp. 2985–2992, June 2009. View at: Publisher Site  Google Scholar
 X. C. Tai, J. Hahn, and G. J. Chung, “A fast algorithm for Euler's elastic model using augmented Lagrangian method,” SIAM Journal on Imaging Sciences, vol. 4, pp. 313–344, 2011. View at: Google Scholar
Copyright
Copyright © 2011 Weili Zeng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.