# A Class of  Fourth-Order Telegraph-Diffusion Equations for Image Restoration

Authors: Weili Zeng; Xiaobo Lu; Xianghua Tan
Journal: Journal of Applied Mathematics
 (2011)

---

## Abstract

A class of nonlinear
fourth-order telegraph-diffusion equations (TDE)
for image restoration are proposed based on
fourth-order TDE and bilateral filtering. The
proposed model enjoys the benefits of both
fourth-order 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.


---

## Body

## 1. Introduction

The use of second-order 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 second-order 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 u denote the image intensity function, t the time, the anisotropic diffusion equation as formulated by Perona and Malik [1] presented as (the PM model) (1.1)ut-∇⋅(g(|∇u|)∇u)=0,
where g is the conductance coefficient, ∇· and ∇ denote the divergence and the gradient, respectively. Equation (1.1) is associated with the following energy functional:(1.2)E(u)=∫Ωf(|∇u|)dΩ,
where Ω is the image support, and f(·)≥0 is an increasing function associated with the diffusion coefficient as (1.3)g(s)=f′(s)s.
The PM model is then associated with an energy-dissipating 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 fourth-order PDE (the YK model)(1.4)ut+∇2(g(|∇2u|)∇2u)=0,
where ∇2 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 fourth-order 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, fourth-order PDE has an inherent limitation. The proposed PDE tends to leave images with speckle artifacts.Recently, Ratner and Zeevi [22] introduced the following telegraph-diffusion equation (TDE model) (1.5)utt+λut-∇⋅(g(|∇u|)∇u)=0,
where λ is the damping coefficient and g 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 g, 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 telegraph-diffusion equation (ITDE model)(1.6)utt+λut-∇⋅(g1(|∇Gσ*u|)∇u)=0,
where Gσ 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 fourth-order telegraph-diffusion equation as follows: (1.7)utt+λut+∇2(g2(|∇2u|)∇2u)=0,
where ∇2 denotes Laplace operator. To our best knowledge, second-order 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:(1.8)utt+λut+∇2(g3(|∇2Bσ(u)|)∇2u)=0,
where Bσ is a bilateral filtering [24, 25]; namely,(1.9)Bσ(u(x))=1W(x)∫ΩGσs(ξ,x)Gσr(u(ξ),u(x))u(ξ)dξ,
with the normalization constant                    (1.10)W(x)=∫ΩGσs(ξ,x)Gσr(u(ξ),u(x))dξ,
where Gσs will be a spatial Gaussian that decreases the influence of distant pixels, while Gσr will be a range Gaussian that decreases the influence of pixels ξ with intensity values that are very different from those of u(x); for example, (1.11)Gσs=exp⁡(-|ξ-x|22σs2),Gσr=exp⁡(-|u(ξ)-u(x)|22σr2),
where parameters σs and σrdictate the amount of filtering applied in the domain and the range of the image, respectively.The ability of edge preservation in the fourth-order TDE-based image restoration method strongly depends on the conductance coefficientg. 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: (1.12)g3(s)=11+(s/k)2,
where k is a threshold constant.We will recall that the main purpose of the functiong 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 Section2, 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:(2.1)utt+λut+∇2(g(|∇2Bσ(u)|)∇2u)=0,(x,t)∈ΩT=Ω×(0,T),(2.2)u(x,0)=u0(x),ut(x,0)=0,x∈Ω,(2.3)u=0,∂u∂n=0,(x,t)∈∂Ω×(0,T),
where Ω is a bounded domain of ℝN with an appropriately smooth boundary, n denotes the unit outer normal to Ω, and T>0. In this section, C 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 denoteHk(Ω) is a Hilbert space for the norm(2.4)‖u‖Hk(Ω):=(∑|s|≤k∫Ω|∂mu∂xm|2dx)1/2,
where k is a positive integer and ∂mu/∂xm of order |m|=∑j=1mmj≤k denotes the distributional derivative of u.We denote byLp(0,T;Hk(Ω)) the set of all functions u such that for almost every t in (0,T), u(t) belongs to Hk(Ω). Lp(0,T;Hk(Ω)) is a normed space for the norm     (2.5)‖u‖Lp(0,T,Hk(Ω)):={(∫0T‖u(t)‖Hk(Ω)p)1/p,(1≤p<∞),esssup⁡0≤t≤T⁡‖u(t)‖Hk(Ω),(p=∞),
where k is a positive integer. L∞(0,T;L2(Ω)) is a normed space for the norm(2.6)‖u‖L∞(0,T,L2(Ω)):=esssup⁡0≤t≤T⁡‖u(t)‖L2(Ω).
We denote by H2(Ω)′ the dual of H2(Ω) and introduce the following function space V(Ω)(2.7)V(Ω)={u∈H2(Ω):u=0,∂u∂n,x∈∂Ω}.
Obviously, V(Ω) is a Banach space with the norm ||·||V(Ω)=||·||H2(Ω).Definition 2.1.
LetT be a fixed positive number. A function u is a weak solution of the problem (2.1)–(2.3) provided (i)
u∈C([0,T];L2(Ω))∩L∞(0,T;V), ut∈L∞(0,T;L2(Ω)) and utt∈L∞(0,T;V′),(ii)
∫Ωuttφ+λ∫Ωutφ+∫Ωg(|ΔBσ*w|)ΔuΔφ=0 for any φ∈V(Ω) and a. e. time 0≤t≤T,(iii)
u(0)=u0,ut(0)=0.


### 2.1. Fourth-Order Linear Equation: Existence and Uniqueness

Now, we consider the existence and uniqueness of weak solutions of the following linear TD problem(2.8)utt+λut+∇2(g(|∇2Bσ(w)|)∇2u)=0,(x,t)∈ΩT=Ω×(0,T),(2.9)u(x,0)=u0(x),ut(x,0)=0,x∈Ω,(2.10)u=0,∂u∂n=0,(x,t)∈∂Ω×(0,T),
where w∈L∞([0,T];L2(Ω)).Definition 2.2.
A functionuw is called a weak solution of the problem (2.8)–(2.10) provided (i)
uw∈C([0,T];L2(Ω))∩L∞(0,T;V), (uw)t∈L∞(0,T;L2(Ω)) and (uw)tt∈L2(0,T;V′),(ii)
∫Ω(uw)ttφdx+λ∫Ω(uw)tφdx+∫Ωg(|∇2Bσ*w|)∇2uw∇2φdx=0for any φ∈V(Ω) and a. e. 0≤t≤T,(iii)
u(0)=u0,ut(0)=0.Theorem 2.3.
Suppose thatw∈L∞([0,T];L2(Ω)) and u0∈H2(Ω), then the problem (2.8)–(2.10) admits a unique weak solution uw.In order to prove Theorem2.3, we will prove the following two lemmas.Lemma 2.4.
Suppose thatu∈L2(Ω), then there exists a constant C1>0, depending only on σ, |Ω| and N, such that
(2.11)‖∇2Bσ(u)‖L∞(Ω)≤C1‖u‖L2(Ω).
Moreover, there exists a constant C2>0, depending only on σ, |Ω|, n such that
(2.12)‖g(|∇2Bσ(u)|)-g(|∇2Bσ(ν)|)‖L∞(Ω)≤C2‖u-ν‖L2(Ω),
for each u,ν∈L2(Ω).Proof.
By the definition ofBσ in (1.8), we then calculate
(2.13)‖∇2Bσ(u)‖L∞(Ω)=1W(x)∫Ω|exp⁡(-|ξ-x|2σs2)||exp⁡(-|u(ξ)-u(x)|2σr2)||u(ξ)|dξ≤C(|Ω|,σ,N)‖u‖L2(Ω)∫Ω|exp⁡(-|ξ-x|2σs2)||exp⁡(-|u(ξ)-u(x)|2σr2)|dξ≤C1‖u‖L2(Ω).
Moreover, since g(s) and Bσ are smooth, we have
(2.14)‖g(|∇2Bσ(u)|)-g(|∇2Bσ(ν)|)‖L∞(Ω)≤C‖|∇2Bσ(u)|2-|∇2Bσ(ν)|2‖L∞(Ω)≤C‖|∇2Bσ(u)|+|∇2Bσ(ν)|‖L∞(Ω)‖∇2Bσ(u-ν)‖L∞(Ω)=C(|Ω|,σ,N)(‖u‖L2(Ω)+‖ν‖L2(Ω))‖u-ν‖L2(Ω)≤C2‖(u-ν)‖L2(Ω).
Therefore, we draw the conclusion.Lemma 2.5.
If there existsM>0 such that
(2.15)‖w‖L∞(0,T;L2(Ω))+‖wt‖L∞(0,T;L2(Ω))≤M.
Then, one has the estimate
(2.16)‖u(t)‖L∞(0,T;H2(Ω))+‖ut(t)‖L∞(0,T;L2(Ω))+‖utt‖L2(0,T;V′)≤C‖u0‖H2(Ω).Proof.
Multiplying (2.8) by ut and integrating over Ω yields
(2.17)12ddt∫Ω(ut)2dx+λ∫Ω(ut)2dx+∫Ωg(|∇2Bσ(w)|)∇2u∇2utdx=0,
for a.e. 0≤t≤T.
Note that(2.18)∫Ωg(|∇2Bσ(w)|)ΔuΔutdx=ddt(12∫Ωg(|∇2Bσ(w)|)(∇2u)2dx)-12∫Ωgt(|∇2Bσ(w)|)(∇2u)2dx.
Since w and wt satisfy (2.15), then ∇2Bσ(w), ∇2Bσ(wt) belong to L∞(0,T;L∞(Ω)) and there exists a constant C depending on Bσ and Ω such that
(2.19)‖∇2Bσ(w)‖L∞(0,T;L∞(Ω)≤C‖u0‖H2(Ω),‖∇2Bσ(wt)‖L∞(0,T;L∞(Ω)≤C‖u0‖H2(Ω),
for any x∈V, a.e. 0≤t≤T. Therefore,since g is decreasing and t→g(t) is smooth, there exist two constants α, β>0 such that
(2.20)α≤g(|∇2Bσ(w)|)≤1,|gt(|∇2Bσ(w)|)|≤β‖u0‖H2(Ω),
for any x∈Ω, a.e. 0≤t≤T.
Consequently, (2.18) yields
(2.21)∫Ωg(|∇2Bσ*w|)ΔuΔutdx≥α2ddt(∫Ω(∇2u)2dx)-C‖u‖H2(Ω)2.
In terms of (2.17) and (2.21), it then follows that
(2.22)ddt(‖ut‖L2(Ω)2+‖u‖H2(Ω)2)≤C(‖ut‖L2(Ω)2+‖u‖H2(Ω)2),
for a.e. 0≤t≤T.
Now, write(2.23)η(t)=‖ut‖L2(Ω)2+‖u‖H2(Ω)2.
Then, (2.22) reads
(2.24)η′(t)≤Cη(t),
for a.e. 0≤t≤T. Applying the Gronwall inequality to (2.24) yields the estimate
(2.25)η(t)≤eCtη(0),0≤t≤T.
Since
(2.26)η(0)=‖ut(0)‖L2(Ω)2+‖u(0)‖H2(Ω)2=‖u(0)‖H2(Ω)2,
we obtain from (2.23)–(2.25) the estimate
(2.27)‖u(t)‖L∞(0,T;H2(Ω))+‖ut(t)‖L∞(0,T;L2(Ω))≤C‖u0‖H2(Ω)2.
Multiplying (2.8) by φ and integrating over Ω, we deduce for a.e. 0≤t≤T that
(2.28)∫Ωuttφdx+λ∫Ωutφdx+∫Ωg(|∇2Bσ*w|)∇2u∇2φdx=0,
where φ∈H02(Ω) and ∥φ∥H02(Ω)≤1.
Consequently,(2.29)‖utt‖V′≤C‖u‖H2(Ω).
Therefore,
(2.30)∫0T‖utt‖V′2dt≤C∫0T‖u‖H2(Ω)2dt≤C‖u0‖H2(Ω)2.
This completes the proof of the lemma.Remark 2.6.
By Lemma2.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. Fourth-Order Nonlinear Equation: Existence and Uniqueness

Theorem 2.7 (see ([26], Schauder’s Fixed Point Theorem)).
Suppose thatX denotes a real Banach space and K⊂X is compact and convex, and assume also that
(2.31)S:K⟶X
is continuous. Then, S has a fixed point in K.Theorem 2.8.
Suppose thatu0∈H2(Ω), 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
(2.32)W(0,T)={w∈L∞(0,T;H2(Ω)),wt∈L∞(0,T;L2(Ω)),wtt∈L∞(0,T;H2(Ω)′)}.
Obviously, W(0,T) is a Banach space with the norm
(2.33)‖w‖W(0,T)=‖w‖L∞(0,T;H2(Ω))+‖wt‖L∞(0,T;L2(Ω))+‖wtt‖L∞(0,T;H2(Ω)′).
Recalling Lemma 2.5, we introduce the subspace W0 of W(0,T) defined by
(2.34)W0={w∈W(0,T);w(0)=u0,wt(0)=0,‖w‖L∞(0,T;L2(Ω))+‖wt‖L∞(0,T;L2(Ω))+‖wtt‖L2(0,T;V′)≤C‖u0‖H2(Ω)}.
By construction, S:w→uw is a mapping from W0 into W0. Since W(0,T) is compactly imbedded in L2(0,T;L2(Ω)), W0 is a nonempty, convex, and weakly compact subset of L2(0,T;L2(Ω)).
In order to apply the Schauder fixed point theorem, we need to prove thatS is a continuous compact mapping from W0 into W0. Let {wk} be a sequence in W0 which converges weakly to some w in W0 and uk=S(wk). We have to prove that S(wk)=uk converges weakly to S(w)=uw. By using the classical theorem of compact inclusion in Sobolev spaces [27], the sequence {wk} contains a subsequence (still denoted by) {wk} and {uk} contains a subsequence (still denoted by) {uk} such that
(2.35)wk⟶winL2(0,T;L2(Ω))anda.e.on⁡Ω×(0,T),(2.36)g(|∇2Bσ*wk|)⟶g(|∇2Bσ*w|)inL2(0,T;L2(Ω))anda.e.on⁡ΩT,(2.37)uk⟶uinL2(0,T;H2(Ω)),(2.38)(uk)t⟶(u)tweaklyinL2(0,T;L2(Ω)),(2.39)(uk)tt⟶(u)ttweaklyinL2(0,T;V′),(2.40)∇2uk⟶∇2uweaklyinL2(0,T;L2(Ω)),(2.41)uk(x,0)⟶u(x,0)inL2(Ω),(2.42)(uk)t(x,0)⟶0inV′.
Since {uk} is the weak solutions of (2.8)–(2.10) corresponding to {wk}, we have, for a.e. 0≤t≤T(2.43)∫Ω(uk)ttφdx+λ∫Ω(uk)tφdx+∫Ωg(|∇2Bσ(wk)|)∇2uk∇2φdx=0.
Integrating (2.43) from 0 to T yields
(2.44)∫0T∫Ω(uk)ttφdxdt+λ∫0T∫Ω(uk)tφdxdt+∫0T∫Ωg(|∇2Bσ(wk)|)∇2uk∇2φdxdt=0.
From (2.38)-(2.39), we have
(2.45)∫0T∫Ω(uk)ttφdxdt⟶∫0T∫Ω(u)ttφdxdt,ask⟶∞,(2.46)λ∫0T∫Ω(uk)tφdxdt⟶λ∫0T∫Ω(u)tφdxdt.
Considering the third term in (2.44)
(2.47)|∫0T∫Ωg(|∇2Bσ(wk)|)∇2uk∇2φdxdt-∫0T∫Ωg(|∇2Bσ(w)|)ΔuΔφdxdt|≤∫0T∫Ω|g(|∇2Bσ(wk)|)-g(|∇2Bσ(w)|)|ΔukΔφdxdt+|∫0T∫Ωg(|∇2Bσ(wk)|)∇2(uk-u)∇2φdxdt|=A+B.
According to Lemma 2.4 and (2.40), we deduce, as k→∞(2.48)A≤C∫0T∫Ω‖wk-w‖L2(Ω)|Δuk||Δφ|≤C(φ)‖wk-w‖L2(Ω)‖Δuk‖L2(0,T;L2(Ω))⟶0,
and note also, as k→∞, B→∞.
Lettingk→∞ in (2.44) yields
(2.49)∫0T∫Ω(u)ttφdxdt+λ∫0T∫Ω(u)tφdxdt+∫0T∫Ωg(|∇2Bσ(w)|)∇2u∇2φdxdt=0.
Therefore, by the arbitrariness of φ∈V, we obtain, for a.e. t∈[0,T](2.50)∫Ω(u)ttφdx+λ∫Ω(u)tφdx+∫Ωg(|∇2Bσ(w)|)∇2u∇2φdx=0,
which implies
(2.51)u=uw=S(w).
By the uniqueness of the solution of linear problem (2.8)–(2.10) and (2.37), the sequence uk=S(wk) converges weakly in W0 to u=S(w). Hence, the mapping S is weakly continuous from W0 into W0. This in turn shows that mapping S is compact. A similar argument shows that S is a continuous mapping. Applying the Schauder fixed point theorem, we conclude that S has a fixed point u=S(u), which consequently solve (2.1)–(2.3). Using the classical theory of parabolic equations and the bootstrap argument [26], we can deduce that u is a strong solution of (1.1)–(1.4) and u∈C∞((0,T)×Ω).(2) Uniqueness of the Solution
Now, we turn to the proof of the uniqueness, following the idea in [22]. Let u1 and u2 be two solutions of (2.1)–(2.3), and we have, for almost every t in [0,T] and i=1,2(2.52)(u1-u2)tt+λ(u1-u2)t+∇2(β1∇2(u1-u2))=∇2((β1-β2)∇2u2),(x,t)∈ΩT,(2.53)ui(x,0)=u0(x),∂ui∂t(x,0)=0,x∈Ω,(2.54)ui=0,∂ui∂n=0,(x,t)∈∂Ω×(0,T),
in the distribution sense, where
(2.55)βi(t)=g(|∇2Bσ*ui|).
It suffices to show that u1-u2≡0. To verify this, fix 0<s<T and set
(2.56)νi(t)={∫tsui(τ)dτ,0<t≤s,0,s≤t<T,
for i=1,2. Then, νi(t)∈H02(Ω) for each 0≤t≤T. Multiplying both sides of (2.52) by ν1-ν2, and integrating over Ω×(0,s) yields
(2.57)∫0s∫Ω(-(u1-u2)t(ν1-ν2)t-λ(u1-u2)(ν1-ν2)t+β1∇2(u1-u2)∇2(ν1-ν2))dxdt=∫0s∫Ω(β1-β2)∇2u2∇2(ν1-ν2)dxdt.
Now, (νi)t=-u (0<t<), and so
(2.58)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+12∫Ωβ1(⋅,0)|∇2(ν1-ν2)(⋅,0)|2dx=∫0s(∫Ω(β1-β2)∇2u2∇2(ν1-ν2)dx+12∫Ω|∇2(ν1-ν2)|2∂β1∂tdx)dt.
As seen in the proof of Lemma 2.5, there exists positive constants C3 and C4 are positive constant depending on Bσ,Ω,Tand ||u0||H2(Ω) such that
(2.59)C3≤βi(⋅,0)≤1,C3≤βi≤1,|(β1)t|≤C4,x∈Ω,a.e.t∈(0,T),i=1,2,
which imply from (2.58),
(2.60)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+C32∫Ω|∇2(ν1-ν2)(⋅,0)|2dx≤∫0s(‖β1-β2‖L∞(Ω)(∫Ω|∇2u2|2dx)1/2(∫Ω|∇2(ν1-ν2)|2dx)1/2+C42∫Ω|∇2(ν1-ν2)|2dx)dt.
By using the Young inequality to (2.60), we obtain
(2.61)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+C32∫Ω|∇2(ν1-ν2)(⋅,0)|2dx≤C∫0s((∫Ω|u1-u2|2dx)1/2(∫Ω|∇2(ν1-ν2)|2dx)1/2+∫Ω|∇2(ν1-ν2)|2dx)dt≤C∫0s(∫Ω|u1-u2|2dx+∫Ω|∇2(ν1-ν2)|2dx)dt.
Now, let us write
(2.62)wi(⋅,t)=∫0tui(⋅,τ)dτ,0<t<T,
whereupon (2.61) can be rewritten as
(2.63)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+C32∫Ω|∇2(w1-w2)(⋅,s)|2dx≤C∫0s(∫Ω|u1-u2|2dx+∫Ω|∇2(w1-w2)(⋅,s)-∇2(w1-w2)(⋅,t)|2dx)dt.
Note that
(2.64)‖∇2(w1-w2)(⋅,s)-∇2(w1-w2)(⋅,t)‖L2(Ω)2≤2‖∇2(w1-w2)(⋅,s)‖L2(Ω)2+2‖∇2(w1-w2)(⋅,t)‖L2(Ω)2.
Therefore, (2.63) implies
(2.65)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+C32∫Ω|∇2(w1-w2)(⋅,s)|2dx≤C∫0s(∫Ω|u1-u2|2dx+2∫Ω|∇2(w1-w2)(⋅,t)|2dx)dt+2Cs‖∇2(w1-w2)(⋅,s)‖L2(Ω)2.
Choose T1 so small such that
(2.66)C32-2CT1≥C34.
Then, if 0<s≤T1, we have
(2.67)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+∫Ω|∇2(w1-w2)(⋅,s)|2dx≤C∫0s(∫Ω|u1-u2|2dx+2∫Ω|∇2(w1-w2)(⋅,t)|2dx)dt.
Consequently, the integral form of the Gronwall inequality implies u1-u2≡0 on (0,T1].
Finally, we apply the same argument on the intervals(T1,T2], (2T1,3T1], and so forth and eventually deduce that u1≡u2 on (0,T). Thus, we obtain the uniqueness of weak solutions.


## 2.1. Fourth-Order Linear Equation: Existence and Uniqueness

Now, we consider the existence and uniqueness of weak solutions of the following linear TD problem(2.8)utt+λut+∇2(g(|∇2Bσ(w)|)∇2u)=0,(x,t)∈ΩT=Ω×(0,T),(2.9)u(x,0)=u0(x),ut(x,0)=0,x∈Ω,(2.10)u=0,∂u∂n=0,(x,t)∈∂Ω×(0,T),
where w∈L∞([0,T];L2(Ω)).Definition 2.2.
A functionuw is called a weak solution of the problem (2.8)–(2.10) provided (i)
uw∈C([0,T];L2(Ω))∩L∞(0,T;V), (uw)t∈L∞(0,T;L2(Ω)) and (uw)tt∈L2(0,T;V′),(ii)
∫Ω(uw)ttφdx+λ∫Ω(uw)tφdx+∫Ωg(|∇2Bσ*w|)∇2uw∇2φdx=0for any φ∈V(Ω) and a. e. 0≤t≤T,(iii)
u(0)=u0,ut(0)=0.Theorem 2.3.
Suppose thatw∈L∞([0,T];L2(Ω)) and u0∈H2(Ω), then the problem (2.8)–(2.10) admits a unique weak solution uw.In order to prove Theorem2.3, we will prove the following two lemmas.Lemma 2.4.
Suppose thatu∈L2(Ω), then there exists a constant C1>0, depending only on σ, |Ω| and N, such that
(2.11)‖∇2Bσ(u)‖L∞(Ω)≤C1‖u‖L2(Ω).
Moreover, there exists a constant C2>0, depending only on σ, |Ω|, n such that
(2.12)‖g(|∇2Bσ(u)|)-g(|∇2Bσ(ν)|)‖L∞(Ω)≤C2‖u-ν‖L2(Ω),
for each u,ν∈L2(Ω).Proof.
By the definition ofBσ in (1.8), we then calculate
(2.13)‖∇2Bσ(u)‖L∞(Ω)=1W(x)∫Ω|exp⁡(-|ξ-x|2σs2)||exp⁡(-|u(ξ)-u(x)|2σr2)||u(ξ)|dξ≤C(|Ω|,σ,N)‖u‖L2(Ω)∫Ω|exp⁡(-|ξ-x|2σs2)||exp⁡(-|u(ξ)-u(x)|2σr2)|dξ≤C1‖u‖L2(Ω).
Moreover, since g(s) and Bσ are smooth, we have
(2.14)‖g(|∇2Bσ(u)|)-g(|∇2Bσ(ν)|)‖L∞(Ω)≤C‖|∇2Bσ(u)|2-|∇2Bσ(ν)|2‖L∞(Ω)≤C‖|∇2Bσ(u)|+|∇2Bσ(ν)|‖L∞(Ω)‖∇2Bσ(u-ν)‖L∞(Ω)=C(|Ω|,σ,N)(‖u‖L2(Ω)+‖ν‖L2(Ω))‖u-ν‖L2(Ω)≤C2‖(u-ν)‖L2(Ω).
Therefore, we draw the conclusion.Lemma 2.5.
If there existsM>0 such that
(2.15)‖w‖L∞(0,T;L2(Ω))+‖wt‖L∞(0,T;L2(Ω))≤M.
Then, one has the estimate
(2.16)‖u(t)‖L∞(0,T;H2(Ω))+‖ut(t)‖L∞(0,T;L2(Ω))+‖utt‖L2(0,T;V′)≤C‖u0‖H2(Ω).Proof.
Multiplying (2.8) by ut and integrating over Ω yields
(2.17)12ddt∫Ω(ut)2dx+λ∫Ω(ut)2dx+∫Ωg(|∇2Bσ(w)|)∇2u∇2utdx=0,
for a.e. 0≤t≤T.
Note that(2.18)∫Ωg(|∇2Bσ(w)|)ΔuΔutdx=ddt(12∫Ωg(|∇2Bσ(w)|)(∇2u)2dx)-12∫Ωgt(|∇2Bσ(w)|)(∇2u)2dx.
Since w and wt satisfy (2.15), then ∇2Bσ(w), ∇2Bσ(wt) belong to L∞(0,T;L∞(Ω)) and there exists a constant C depending on Bσ and Ω such that
(2.19)‖∇2Bσ(w)‖L∞(0,T;L∞(Ω)≤C‖u0‖H2(Ω),‖∇2Bσ(wt)‖L∞(0,T;L∞(Ω)≤C‖u0‖H2(Ω),
for any x∈V, a.e. 0≤t≤T. Therefore,since g is decreasing and t→g(t) is smooth, there exist two constants α, β>0 such that
(2.20)α≤g(|∇2Bσ(w)|)≤1,|gt(|∇2Bσ(w)|)|≤β‖u0‖H2(Ω),
for any x∈Ω, a.e. 0≤t≤T.
Consequently, (2.18) yields
(2.21)∫Ωg(|∇2Bσ*w|)ΔuΔutdx≥α2ddt(∫Ω(∇2u)2dx)-C‖u‖H2(Ω)2.
In terms of (2.17) and (2.21), it then follows that
(2.22)ddt(‖ut‖L2(Ω)2+‖u‖H2(Ω)2)≤C(‖ut‖L2(Ω)2+‖u‖H2(Ω)2),
for a.e. 0≤t≤T.
Now, write(2.23)η(t)=‖ut‖L2(Ω)2+‖u‖H2(Ω)2.
Then, (2.22) reads
(2.24)η′(t)≤Cη(t),
for a.e. 0≤t≤T. Applying the Gronwall inequality to (2.24) yields the estimate
(2.25)η(t)≤eCtη(0),0≤t≤T.
Since
(2.26)η(0)=‖ut(0)‖L2(Ω)2+‖u(0)‖H2(Ω)2=‖u(0)‖H2(Ω)2,
we obtain from (2.23)–(2.25) the estimate
(2.27)‖u(t)‖L∞(0,T;H2(Ω))+‖ut(t)‖L∞(0,T;L2(Ω))≤C‖u0‖H2(Ω)2.
Multiplying (2.8) by φ and integrating over Ω, we deduce for a.e. 0≤t≤T that
(2.28)∫Ωuttφdx+λ∫Ωutφdx+∫Ωg(|∇2Bσ*w|)∇2u∇2φdx=0,
where φ∈H02(Ω) and ∥φ∥H02(Ω)≤1.
Consequently,(2.29)‖utt‖V′≤C‖u‖H2(Ω).
Therefore,
(2.30)∫0T‖utt‖V′2dt≤C∫0T‖u‖H2(Ω)2dt≤C‖u0‖H2(Ω)2.
This completes the proof of the lemma.Remark 2.6.
By Lemma2.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. Fourth-Order Nonlinear Equation: Existence and Uniqueness

Theorem 2.7 (see ([26], Schauder’s Fixed Point Theorem)).
Suppose thatX denotes a real Banach space and K⊂X is compact and convex, and assume also that
(2.31)S:K⟶X
is continuous. Then, S has a fixed point in K.Theorem 2.8.
Suppose thatu0∈H2(Ω), 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
(2.32)W(0,T)={w∈L∞(0,T;H2(Ω)),wt∈L∞(0,T;L2(Ω)),wtt∈L∞(0,T;H2(Ω)′)}.
Obviously, W(0,T) is a Banach space with the norm
(2.33)‖w‖W(0,T)=‖w‖L∞(0,T;H2(Ω))+‖wt‖L∞(0,T;L2(Ω))+‖wtt‖L∞(0,T;H2(Ω)′).
Recalling Lemma 2.5, we introduce the subspace W0 of W(0,T) defined by
(2.34)W0={w∈W(0,T);w(0)=u0,wt(0)=0,‖w‖L∞(0,T;L2(Ω))+‖wt‖L∞(0,T;L2(Ω))+‖wtt‖L2(0,T;V′)≤C‖u0‖H2(Ω)}.
By construction, S:w→uw is a mapping from W0 into W0. Since W(0,T) is compactly imbedded in L2(0,T;L2(Ω)), W0 is a nonempty, convex, and weakly compact subset of L2(0,T;L2(Ω)).
In order to apply the Schauder fixed point theorem, we need to prove thatS is a continuous compact mapping from W0 into W0. Let {wk} be a sequence in W0 which converges weakly to some w in W0 and uk=S(wk). We have to prove that S(wk)=uk converges weakly to S(w)=uw. By using the classical theorem of compact inclusion in Sobolev spaces [27], the sequence {wk} contains a subsequence (still denoted by) {wk} and {uk} contains a subsequence (still denoted by) {uk} such that
(2.35)wk⟶winL2(0,T;L2(Ω))anda.e.on⁡Ω×(0,T),(2.36)g(|∇2Bσ*wk|)⟶g(|∇2Bσ*w|)inL2(0,T;L2(Ω))anda.e.on⁡ΩT,(2.37)uk⟶uinL2(0,T;H2(Ω)),(2.38)(uk)t⟶(u)tweaklyinL2(0,T;L2(Ω)),(2.39)(uk)tt⟶(u)ttweaklyinL2(0,T;V′),(2.40)∇2uk⟶∇2uweaklyinL2(0,T;L2(Ω)),(2.41)uk(x,0)⟶u(x,0)inL2(Ω),(2.42)(uk)t(x,0)⟶0inV′.
Since {uk} is the weak solutions of (2.8)–(2.10) corresponding to {wk}, we have, for a.e. 0≤t≤T(2.43)∫Ω(uk)ttφdx+λ∫Ω(uk)tφdx+∫Ωg(|∇2Bσ(wk)|)∇2uk∇2φdx=0.
Integrating (2.43) from 0 to T yields
(2.44)∫0T∫Ω(uk)ttφdxdt+λ∫0T∫Ω(uk)tφdxdt+∫0T∫Ωg(|∇2Bσ(wk)|)∇2uk∇2φdxdt=0.
From (2.38)-(2.39), we have
(2.45)∫0T∫Ω(uk)ttφdxdt⟶∫0T∫Ω(u)ttφdxdt,ask⟶∞,(2.46)λ∫0T∫Ω(uk)tφdxdt⟶λ∫0T∫Ω(u)tφdxdt.
Considering the third term in (2.44)
(2.47)|∫0T∫Ωg(|∇2Bσ(wk)|)∇2uk∇2φdxdt-∫0T∫Ωg(|∇2Bσ(w)|)ΔuΔφdxdt|≤∫0T∫Ω|g(|∇2Bσ(wk)|)-g(|∇2Bσ(w)|)|ΔukΔφdxdt+|∫0T∫Ωg(|∇2Bσ(wk)|)∇2(uk-u)∇2φdxdt|=A+B.
According to Lemma 2.4 and (2.40), we deduce, as k→∞(2.48)A≤C∫0T∫Ω‖wk-w‖L2(Ω)|Δuk||Δφ|≤C(φ)‖wk-w‖L2(Ω)‖Δuk‖L2(0,T;L2(Ω))⟶0,
and note also, as k→∞, B→∞.
Lettingk→∞ in (2.44) yields
(2.49)∫0T∫Ω(u)ttφdxdt+λ∫0T∫Ω(u)tφdxdt+∫0T∫Ωg(|∇2Bσ(w)|)∇2u∇2φdxdt=0.
Therefore, by the arbitrariness of φ∈V, we obtain, for a.e. t∈[0,T](2.50)∫Ω(u)ttφdx+λ∫Ω(u)tφdx+∫Ωg(|∇2Bσ(w)|)∇2u∇2φdx=0,
which implies
(2.51)u=uw=S(w).
By the uniqueness of the solution of linear problem (2.8)–(2.10) and (2.37), the sequence uk=S(wk) converges weakly in W0 to u=S(w). Hence, the mapping S is weakly continuous from W0 into W0. This in turn shows that mapping S is compact. A similar argument shows that S is a continuous mapping. Applying the Schauder fixed point theorem, we conclude that S has a fixed point u=S(u), which consequently solve (2.1)–(2.3). Using the classical theory of parabolic equations and the bootstrap argument [26], we can deduce that u is a strong solution of (1.1)–(1.4) and u∈C∞((0,T)×Ω).(2) Uniqueness of the Solution
Now, we turn to the proof of the uniqueness, following the idea in [22]. Let u1 and u2 be two solutions of (2.1)–(2.3), and we have, for almost every t in [0,T] and i=1,2(2.52)(u1-u2)tt+λ(u1-u2)t+∇2(β1∇2(u1-u2))=∇2((β1-β2)∇2u2),(x,t)∈ΩT,(2.53)ui(x,0)=u0(x),∂ui∂t(x,0)=0,x∈Ω,(2.54)ui=0,∂ui∂n=0,(x,t)∈∂Ω×(0,T),
in the distribution sense, where
(2.55)βi(t)=g(|∇2Bσ*ui|).
It suffices to show that u1-u2≡0. To verify this, fix 0<s<T and set
(2.56)νi(t)={∫tsui(τ)dτ,0<t≤s,0,s≤t<T,
for i=1,2. Then, νi(t)∈H02(Ω) for each 0≤t≤T. Multiplying both sides of (2.52) by ν1-ν2, and integrating over Ω×(0,s) yields
(2.57)∫0s∫Ω(-(u1-u2)t(ν1-ν2)t-λ(u1-u2)(ν1-ν2)t+β1∇2(u1-u2)∇2(ν1-ν2))dxdt=∫0s∫Ω(β1-β2)∇2u2∇2(ν1-ν2)dxdt.
Now, (νi)t=-u (0<t<), and so
(2.58)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+12∫Ωβ1(⋅,0)|∇2(ν1-ν2)(⋅,0)|2dx=∫0s(∫Ω(β1-β2)∇2u2∇2(ν1-ν2)dx+12∫Ω|∇2(ν1-ν2)|2∂β1∂tdx)dt.
As seen in the proof of Lemma 2.5, there exists positive constants C3 and C4 are positive constant depending on Bσ,Ω,Tand ||u0||H2(Ω) such that
(2.59)C3≤βi(⋅,0)≤1,C3≤βi≤1,|(β1)t|≤C4,x∈Ω,a.e.t∈(0,T),i=1,2,
which imply from (2.58),
(2.60)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+C32∫Ω|∇2(ν1-ν2)(⋅,0)|2dx≤∫0s(‖β1-β2‖L∞(Ω)(∫Ω|∇2u2|2dx)1/2(∫Ω|∇2(ν1-ν2)|2dx)1/2+C42∫Ω|∇2(ν1-ν2)|2dx)dt.
By using the Young inequality to (2.60), we obtain
(2.61)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+C32∫Ω|∇2(ν1-ν2)(⋅,0)|2dx≤C∫0s((∫Ω|u1-u2|2dx)1/2(∫Ω|∇2(ν1-ν2)|2dx)1/2+∫Ω|∇2(ν1-ν2)|2dx)dt≤C∫0s(∫Ω|u1-u2|2dx+∫Ω|∇2(ν1-ν2)|2dx)dt.
Now, let us write
(2.62)wi(⋅,t)=∫0tui(⋅,τ)dτ,0<t<T,
whereupon (2.61) can be rewritten as
(2.63)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+C32∫Ω|∇2(w1-w2)(⋅,s)|2dx≤C∫0s(∫Ω|u1-u2|2dx+∫Ω|∇2(w1-w2)(⋅,s)-∇2(w1-w2)(⋅,t)|2dx)dt.
Note that
(2.64)‖∇2(w1-w2)(⋅,s)-∇2(w1-w2)(⋅,t)‖L2(Ω)2≤2‖∇2(w1-w2)(⋅,s)‖L2(Ω)2+2‖∇2(w1-w2)(⋅,t)‖L2(Ω)2.
Therefore, (2.63) implies
(2.65)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+C32∫Ω|∇2(w1-w2)(⋅,s)|2dx≤C∫0s(∫Ω|u1-u2|2dx+2∫Ω|∇2(w1-w2)(⋅,t)|2dx)dt+2Cs‖∇2(w1-w2)(⋅,s)‖L2(Ω)2.
Choose T1 so small such that
(2.66)C32-2CT1≥C34.
Then, if 0<s≤T1, we have
(2.67)12∫Ω|(u1-u2)(⋅,s)|2dx+∫0s∫Ωλ|(u1-u2)|2dxdt+∫Ω|∇2(w1-w2)(⋅,s)|2dx≤C∫0s(∫Ω|u1-u2|2dx+2∫Ω|∇2(w1-w2)(⋅,t)|2dx)dt.
Consequently, the integral form of the Gronwall inequality implies u1-u2≡0 on (0,T1].
Finally, we apply the same argument on the intervals(T1,T2], (2T1,3T1], and so forth and eventually deduce that u1≡u2 on (0,T). 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 h and a time step size of τ, we quantize the space and time coordinates as follows: (3.1)t=n*τ,n=0,1,2,…,x=i*h,i=0,1,2,…,N,y=j*h,j=0,1,2,…,M,whereM×N is the size of the image. Let ui,jn be the approximation to the value u(ih,jh,nτ). We then use a five-step approach to calculate (2.1)–(2.3).(a) calculating the Laplace of the image intensity functions(3.2)∇2ui,jn=0.5[ui+1,jn+ui-1,jn+ui,j-1n+ui,j+1n-4ui,jnh2]+0.25[ui+1,j+1n+ui-1,j+1n+ui+1,j-1n+ui-1,j-1n-4ui,jnh2],
with symmetric boundary conditions (3.3)u-1,jn=u0,jn,uN+1,jn=uN,jn,j=0,1,2,…,N,ui,-1n=ui,0n,ui,M+1n=ui,Mn,i=0,1,2,…,M,(b) calculating the value of the following function(3.4)φi,jn=c2(|∇2ui,jn|)∇2ui,jn,(c) calculating the Laplace of φi,jn as(3.5)∇2φi,jn=0.5[φi+1,jn+φi-1,jn+φi,j-1n+φi,j+1n-4φi,jnh2]+0.25[φi+1,j+1n+φi-1,j+1n+φi+1,j-1n+φi-1,j-1n-4φi,jnh2],
with symmetry boundary conditions(3.6)φ-1,jn=φ0,jn,φN+1,jn=φN,jn,j=0,1,2,…,N,φi,-1n=φi,0n,φi,M+1n=φi,Mn,i=0,1,2,…,M,(d) defining the discrete approximation:(3.7)δui,jn=ui,jn-ui,jn-1τ,δ2ui,jn=δui,jn-δui,jn-1τ,
with (3.8)ui,j0=ui,j-1,i=0,1,2,…,M,j=0,1,2,…,N,(e) finally, the numerical approximation to the differential equation (2.1) is given as (3.9)αδ2ui,jn+βδui,jn=-∇2φi,jn-1.


## 4. Experimental Results

In this section, we present numerical results obtained by applying our proposed fourth-order TDE to image denosing. We test the proposed method on “Barbara” image with size225×225 (taken from USC-SIPI image database) and “license plate” image with size 240×306. 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(4.1)ISNR=10⋅log⁡10(∑i,j[u(i,j)-u0(i,j)]2∑i,j[u(i,j)-unew(i,j)]2),
where u0(·) is the initial image (noised image) and unew(·) is the denoised image. The value of ISNR is large, and the restored image is better.Comparison of different methods on “Barbara” image. (a) Original image. (b) Noised image. (c) Fourth-order TDE withλ=1. (d) Fourth-order TDE with λ=40. (e) Fourth-order TDE with λ=70. (f) Fourth-order TDE with λ=100.
(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 fourth-order TDE (k=0.5, a=4) with λ=1, 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 λ=40. The INSR value using fourth-order TDE is lower than INSR value using proposed domain-based fourth-order PDE at λ<9.5. When λ→100, the ISNR with fourth-order TDE is almost equal to the ISNR with domain-based fourth-order PDE.Figure 2
Graph of ISNR versus different damping coefficientλ for “Barbara”.Comparison of different methods on “License plate” image. (a) Original image. (b) Noisy image. (c) PM smodel withk=10. (d) TDE model with k=10, λ=40. (e) ITDE model with k=10, λ=40, σ1=0.1. (f) Proposed fourth-order TDE with k=1, a=4,λ=10.
(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σ=20. To verify the effectiveness of our proposed fourth-order 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 240×306. We then added random white Gaussian noise of variance σ=20 to generate 50 degraded image. One such degraded image is shown in Figure 3(b). The results yield by PM second-order PDE is shown in Figure 3(c). We observe that “PM second-order 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 second-order PDE but inherits the blocky effects from PM second-order PDE to some degree. Figure 3(f) is the restored image using our proposed fourth-order TDE. For the 50 random generated images, the mean of ISNR values for PM, TDE, ITDE, and our proposed fourth-order 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.Table 1
ISNR (in dB) values for “license plate” image using random white Gaussian noise of varianceσ=20.
method12345678910PM model5.72465.73095.74475.71275.75025.73545.74475.70495.72545.5188TDE model5.83355.85365.84435.87515.82715.67475.83275.85735.87675.8187ITDE model5.91475.91525.92175.91095.91985.91075.91225.90985.91505.9112Fourth-order TDE6.17196.13826.16416.12686.16266.14316.15436.17976.11626.1883Finally, we designed to further evaluate the good behavior of our proposed fourth-order 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 USC-SIPI 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.Table 2
ISNR (in dB) values for “license plate” image using different methods across five noise levels.
MethodWith Gaussian white noiseσ=10σ=15σ=20σ=25σ=30PM model6.99306.56475.72464.98594.3677TDE model7.10326.61435.83355.12914.5038ITDE model7.21656.83725.91475.30224.6195Fourth-order TDE7.36486.92966.17195.55364.8564Table 3
The mean of ISNR (in dB) values for eight standard test images using random white Gaussian noise of variance 20.
MethodsNoise imagesLenaBoatsCameramanPeppersHouseElainePM model4.32754.14874.03265.11204.27265.1737TDE model4.46624.27454.16895.20224.39965.2461ITDE model4.57434.39484.29175.28724.51275.3830Fourth-order TDE4.71254.52224.41345.44314.64755.5705


## 5. Conclusions

A class of nonlinear fourth-order telegraph-diffusion equation (TDE) for image restoration is presented in this paper. The proposed model first extends the second order TDE for image restoration to fourth-order TDE. Moreover, our proposed model combines nonlinear fourth-order 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 fourth-order TDE method.



---

*Source: 240370-2011-08-02.xml*