# Image Enhancement under Data-Dependent Multiplicative Gamma Noise

Authors: Jidesh Pacheeripadikkal; Bini Anattu
Journal: Applied Computational Intelligence and Soft Computing
 (2014)

---

## Abstract

An edge enhancement filter is proposed for denoising
and enhancing images corrupted with data-dependent noise which is
observed to follow a Gamma distribution. The filter is equipped with
three terms designed to perform three different tasks. The first term is
an anisotropic diffusion term which is derived from a locally adaptivep-laplacian functional. The second term is an enhancement term or
a shock term which imparts a shock effect at the edge points making
them sharp. The third term is a reactive term which is derived based
on the maximum a posteriori (MAP) estimator and this term helps
the diffusive term to perform a Gamma distributive data-dependent
multiplicative noise removal from images. And moreover, this reactive
term ensures that deviation of the restored image from the original one
is minimum. This proposed filter is compared with the state-of-the-art
restoration models proposed for data-dependent multiplicative noise.


---

## Body

## 1. Introduction

Image restoration is an important activity in image processing. Many different kinds of degradation models have been discussed in the recent literature. The degradation due to an additive data independent noise was explored extensively in the image denoising literature; see [1–3] for details. Further the degradations due to the combined effect of random noise and linear shift invariant blur are also studied elaborately; see [4–7] for further details. Another category of degradation that is discussed in the recent literatures is the model with multiplicative data-dependent noise with linear shift-invariant blur. Nevertheless, a solution to data-dependent noise model should practically consider the noise distribution as well [8]. Different noise distributions considered under a multiplicative data-dependent noise set-up include Gamma, Poisson, and Gaussian distributions; refer to [8–10] for further details.


## 2. PDE Models for Multiplicative Noise

A common representation for a multiplicative denoising model, found in the image processing literature (assuming a shift-invariant nature of the linear blurring operator), is(1)
u
0
=
k
*
u
×
n
,
where k is a linear shift invariant blurring kernel, “*” denotes a linear convolution operator, and n is a multiplicative data-dependent noise. The distribution of the noise intensity varies with reference to the application under consideration.The first kind of model that was proposed for image restoration assuming a Gaussian multiplicative noise was RLO model by Lions et al. in [10]. The diffusion equation for the model is
(2)
u
t
=
div
⁡
∇
u
|
∇
u
+
β
|
-
λ
k
*
(
u
0
(
k
*
u
)
2
)
(
u
0
(
k
*
u
)
-
1
)
-
μ
k
*
(
u
0
(
k
*
u
)
2
)
,
where λ and μ are the parameters which control the regularization and the data fitting characteristics of the filter. These parameters are updated dynamically during the evolution process. We use the notation u
t to represent the partial derivative ∂
u
/
∂
t and div
⁡ denotes a divergence operator. We invite readers to refer to [10] for details on the update of these parameters and the theoretical analysis of the model. The solution to this model is defined in the space of bounded variations. The model is well suited for data independent Gaussian multiplicative noise with mean 1 and variance σ
2. However, most of the imaging modalities that are in place today do not produce images with Gaussian multiplicative noise. Many of these models produce speckled images whose intensities are Gamma distributed ultrasound and satellite images are good examples.Aubertt and Aujol [8] proposed a model (AA model) for restoring images corrupted by multiplicative Gamma noise. This model adopts the diffusion term from the total variation (TV) regularization model in [6] (ROF) and the reactive term is derived based on the Bayesian maximum a posteriori (MAP) estimator. The Euler Lagrange equation for this model is
(3)
u
t
=
-
div
⁡
(
∇
u
|
∇
u
|
)
+
λ
k
*
(
k
*
u
-
u
0
)
(
k
*
u
)
2
.
Here λ is a regularization parameter. And taking negative of the above equation gives the gradient descent solution:
(4)
u
t
=
div
⁡
(
∇
u
|
∇
u
|
)
-
λ
k
*
(
k
*
u
-
u
0
)
(
k
*
u
)
2
.
Since the diffusion term (the first term in the above equation) is adopted from the TV model, the diffusive behavior of this filter is quite similar to the ROF model; therefore, the process results in formation of staircases in restored images. The reactive term (second term in the above equation) is conditionally convex, and therefore uniqueness of the solution cannot be proven conclusively. Some modifications are suggested for this model in the recent literature. In [11] the authors propose a Weberized diffusive term for improving the intensity of the filtered image and thereby enhancing the visual quality of the output. We recall the Weber’s low; the intensity increment in the restored data is directly proportional to mean intensity of the background (see [11] for details). However, this model retains the reactive term in AA model and therefore the solution may not be unique. Another remarkable modification to AA model is the one proposed in [12]. In this model, the data fitting term is modified so that the functional remains convex under all conditions. Therefore, a unique solution is guaranteed under a gradient descent set-up. Nevertheless, the total variation based diffusive term may eventually form piecewise constant regions in the filtered output. There are few improvements suggested for the AA model to overcome the difficulties due to the presence of second-order TV diffusion term; refer to [13, 14] for details. Moreover, none of the models discussed above could really enhance the image features while denoising images. However, in many applications enhancement of images is vital.


## 3. Image Enhancement Models

Though the models discussed above are capable of restoring images from their blurred and noisy observations, they lack the capability to enhance the images. Image enhancement is an essential step in many image processing applications. Shock filters were introduced in the literature to enhance images from their blurred observations. A classical shock filter introduced by Osher and Rudin in [15] follows the Euler-Lagrange equation:
(5)
u
t
=
-
sign
⁡
(
u
η
η
)
|
∇
u
|
,
where sign
⁡ denotes the mathematical sign
⁡ function, which takes the values from the set {
-
1,0
,
1
}. This filter enhances the high frequency components present in the images; however, the noise components also get enhanced, consequently. Therefore, a preprocessing on the input image becomes inevitable.In another work, Alvarez and Mazorra [16] proposed a diffusion coupled shock filter. This model enhances and denoises simultaneously, so that the noise features do not blow up in the restored version of the image. The model is formulated as
(6)
u
t
=
-
sign
⁡
(
u
η
η
)
|
∇
u
|
+
λ
u
ξ
ξ
,
were ξ denotes the direction normal to the gradient or parallel to the level line and λ is a regularization parameter. The term u
ξ
ξ is mean curvature of the level curve. The first term signifies a shock filter and the second one is a diffusion term. This model is quite effective in enhancing images in case of noise intervention. There are quite a few modifications suggested for this model, as well. The modifications are suggested for both diffusion and the shock terms. We invite readers to refer to [17–19] for details.All these facts discussed above motivated us to propose a filter which enhances the edge features while restoring the images from their blurred and noisy observations, especially when the noise is multiplicative and Gamma distributed. The detailed mathematical formulation and experimental results are given in succeeding sections.


## 4. The Proposed Model

The energy functional for the new model is(7)
min
⁡
u
E
(
u
)
=
∫
Ω
‍
|
∇
u
|
p
(
s
)
d
Ω
+
λ
∫
Ω
‍
(
log
⁡
⁡
u
+
u
0
u
)
d
Ω
,
where u
0 is the observed image, Ω denotes the image domain, and the function p
(
s
) (where s denotes |
∇
u
|) is defined as
(8)
p
(
s
)
=
1
+
1
1
+
s
2
.
This functional is a modified form of the functional proposed by Aubertt and Aujol in [8] and Blomgren et al. in [20]. The reactive term of the proposed model is borrowed from AA model, which is in fact derived based on the Bayesian MAP estimator. The details of the p
-laplacian diffusion term in the above equation can be found in [20]. The proposed functional is not purely convex and the details are provided in the Appendix to follow. The Euler-Lagrange equation can be derived as
(9)
u
t
=
-
∇
·
(
∇
u
|
∇
u
|
2
-
p
(
s
)
)
+
λ
2
(
u
-
u
0
)
u
2
.
The solution to this equation is derived at the steady state (when u
t
=
0). Here we note that this evolution PDE does not have any enhancing capability on its own. Nevertheless, enhancement comes handy in many imaging problems because of the common device artifacts resulting in blurred versions of the captured images. Therefore, we modify the PDE defined in (9) to incorporate the enhancing capabilities as well. The new proposed PDE takes the following form:
(10)
u
t
=
-
∇
·
(
∇
u
|
∇
u
|
2
-
p
(
s
)
)
-
λ
1
2
π
arctan
⁡
(
G
σ
*
u
η
η
k
(
t
)
)
|
∇
u
|
+
λ
2
(
u
-
u
0
)
u
2
,
where λ
1 and λ
2 are the parameters that are evaluated experimentally. Here λ
2 is evaluated as λ
2
=
0.5
/
σ
2, where σ is the noise variance and λ
1 is assigned a value 0.02 experimentally to get the desired enhancement. The term G
σ
u stands for a Gaussian convolved version of the image function u and σ here denotes the standard deviation of the Gaussian function used for smoothing the image. The diffusion term causes a denoising effect, whereas the shock term (the second term in the above equation) is responsible for the enhancement property of the filter (see [18] for details) and the reactive term ensures the deviation from the actual solution to be minimal (so that restoration is more accurate). The parameters λ
1 and λ
2 control the enhancement and the denoising behavior of the filter. The function p
(
·
) as already defined in (8) depends on the gradient of the image function. In the regions where gradient value is zero or in the constant intensity regions the function takes a value 2 and in high gradient regions the functions assume a value 1 and in all other regions it assumes a value between 1 and 2. Since the diffusion process is driven by the value of the function p
(
·
), this diffusion term is commonly known as adaptive p-laplacian term. The p-laplace term reduces the staircase effect due to the prompt switching of the filter between norms in the range [
1
-
2
]. In order to avoid the shock components being active during the initial phases of evolution, a time-dependent function k
(
t
) is used in the shock term. This function helps diffusion process to dominate the evolution during the initial stages of the evolution process; refer to [18, 19] for details of this function. Here we assume that the noise variance and the blurring function are known in advance and the noise distribution is assumed to be Gamma.


## 5. Numerical Implementation

Explicit central difference schemes are employed for representing the finite differences of the terms in the proposed filter in (10) except the shock term, for which upwind [21] scheme is applied. Discretization of diffusion equation is
(11)
u
x
x
(
δ
+
u
y
2
)
-
2
u
x
×
u
y
×
u
x
y
+
u
y
y
(
δ
+
u
x
2
)
(
δ
+
u
x
2
+
u
y
2
)
4
-
(
p
(
s
)
)
/
2
,
where δ is a small positive value used to avoid “division by zero” in homogeneous intensity regions. The term |
∇
u
| in the shock term is discretized using the upwind scheme proposed in [21] to make the process stable
(12)
∥
∇
u
∥
=
D
x
2
+
D
y
2
,
where D
x
=
minmod
(
u
x
+
(
x
,
y
)
,
u
x
-
(
x
,
y
)
),  D
y
=
minmod
(
u
y
+
(
x
,
y
)
,
u
y
-
(
x
,
y
)
), and the minmod operator is defined as
(13)
minmod
(
x
,
y
)
=
{
min
⁡
⁡
(
∥
x
∥
,
∥
y
∥
)
if
x
y
>
0
,
0
Otherwise
.The first-order forward, backward, and central difference operators are as follows:(14)
u
x
+
(
x
,
y
)
=
u
(
x
+
1
,
y
)
-
u
(
x
,
y
)
,
u
x
-
(
x
,
y
)
=
u
(
x
,
y
)
-
u
(
x
-
1
,
y
)
,
u
x
(
x
,
y
)
=
(
u
x
+
+
u
x
-
)
2
,
u
y
+
(
x
,
y
)
=
u
(
x
,
y
+
1
)
-
u
(
x
,
y
)
,
u
y
-
(
x
,
y
)
=
u
(
x
,
y
)
-
u
(
x
,
y
-
1
)
,
u
y
(
x
,
y
)
=
(
u
y
+
+
u
y
-
)
2
.
Here the super scripts +/− denote the forward and backward finite differences, respectively, and nonsuperscripted symbols stand for the finite central differences. Finally the diffusion equation is discretized using the above-mentioned finite difference schemes.


## 6. Experimental Results

We have tested various images with different characteristics; however, the visual results are provided only for two test images “peppers” and “Lena” to maintain the brevity in explanation. However, we confirm that the model works for other images as well. These test images are chosen in such a way that they belong to two different classes. The image “pepper” is a low textured cartoon type image with homogeneous intensity regions and image “Lena” is a partially textured one. The images are normalized to the range[
0
-
1
] in our experiments. The blurring artifact is simulated using a spatially invariant Gaussian kernel with σ value 3. We transform the image to the Fourier domain and multiply it with the frequency response of the Gaussian kernel and transform the image back to the spatial domain to get the blurred version of the input image. Further, the multiplicative Gamma noise (with unit mean) is added with different noise variances to the images under testing. We have generated the Gamma noise (with the specified variance) using a function written in MATLAB. We also note that all the methods considered in this work are simulated (programmed) using MATLAB 7.9 (R2009b) on Linux platform. The process is iterated till the error in consecutive iterations becomes less than a threshold value. The overall algorithm for the model is summarized in Algorithm 1.Algorithm 1:Image enhancement filter.
(1) Input the noisy imageu
0 corrupted by out-of-focus blur and data-dependent
multiplicative Gamma noise.
(2) Calculate the regularization weightλ
2 as 0.5
/
σ
n
2 where σ
n is the noise
standard deviation andλ
1
=
0.02 (this was set empirically).
(3) Specify the stopping threshold valuet.
(4) Select a time step ofΔ
t
=
0.1 for the PDE evolution.
(5) Initialize the current restored imageu
r
n
=
0 and u
n
+
1
r
=
u
0
(6)while  (∥
u
r
n
+
1
-
u
r
n
∥/∥
u
r
n
∥) >
t  do,
(7)   Assignu
r
n
=
u
r
n
+
1
(8)   Evaluateu
r
n
+
1 using the PDE defined in (10).
(9)end whileThe visual results in Figures1 and 2 for the test images “peppers” and “Lena” show the performances of various restoration methods under a multiplicative Gamma noise set-up. Here we consider the multiplicative noise restoration models AA [8] and RLO [10], for demonstrating and comparing the restoration capacities of various filters under consideration and the one proposed in this paper. As evident from these figures, the RLO model is not that effective for the multiplicative Gamma noise, as it is designed mainly for handling Gaussian distributed multiplicative noise. AA model performs better than RLO model; nevertheless, it forms staircases on the sharp edges when it eventually reaches to the steady state, as evident from Figures 1(c), 1(g), 2(c) and 2(g). The proposed model is capable of enhancing and restoring images from their noisy observations without causing any staircase effect; see Figures 1(d), 1(h), 2(d), and 2(h), for a comparative analysis. Moreover, the deblurring/deconvolution capability of the proposed model is evident in these figures. We have tested the models for various noise variance values (the results for two different noise variances and two different blurs are provided in Figures 1 and 2) and different quantity of blurs. In all these different degradation scenarios, the proposed model is observed to perform at a better scale in terms of visual representation compared to the other contemporary and relevant methods discussed in the literature.Images restored using various methods, images of “peppers”: (a) and (e): noisy images with mean 1, noise variances 0.1 and 0.2, and Gaussian blurs with kernel spread of 3 & 4, respectively; (b) and (f) restored using RLO method [10]; (c) and (g) restored using AA filter [8]; (d) and (h) restored using the proposed method.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)Images restored using various methods, images of “Lena”: (a) and (e): noisy images with mean 1, noise variances 0.1 and 0.2, and Gaussian blurs with kernel spread of 3 & 4, respectively; (b) and (f) restored using RLO method [10]; (c) and (g) restored using AA filter [8]; (d) and (h) restored using the proposed method.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)We perform a quantitative study of various restoration models including the one proposed in this work, using standard statistical measures like Pratt's figure of merit (PFOM) [22] and structural similarity index (SSIM) [23]. Table 1 is dedicated for comparing the methods under consideration in terms of these two statistical measures. From this table we may infer that the proposed model has a better performance rate in terms of these statistical measures compared to other methods under consideration. The PFOM value signifies the edge preservation capability of the filter under consideration. The values are in the range [
0-1
]; the value 1 stands for a perfect edge preservation. Finally, SSIM denotes the structure similarity between the restored and the original images. This value also lies in the range [
0-1
]. And the value approaches 1 as the structure preservation approaches the perfection level. The numerical values tabulated for the two different statistical measures prove beyond doubt that the proposed model outperforms other models considered here, in terms of edge and structure preservation. The figures shown are highly in favor of the noise removal and edge and structure preservation capability of the proposed model.Table 1
Statistical measures {PFOM, SSIM} plotted for different images for different models.
Images
RLO
AA
The proposed model
Lena
{0.73, 0.71}
{0.79, 0.76}
{0.88, 0.86}
Cameraman
{0.71, 0.68}
{0.79, 0.74}
{0.91, 0.88}
Peppers
{0.74, 0.71}
{0.81, 0.77}
{0.93, 0.91}
Woman
{0.73, 0.69}
{0.80, 0.79}
{0.90, 0.87}
Baboon
{0.68, 0.67}
{0.73, 0.71}
{0.82,0.78}


## 7. Conclusion

In this work a filter is proposed for enhancing the images while restoring them from their blurred and noisy observations. The images are assumed to be corrupted with multiplicative Gamma noise and shift invariant Gaussian blur. Further the presence of a switchingp-laplacian diffusion term helps in reducing the staircase effect, which otherwise appears in all second-order nonlinear PDE models. The Experiments are performed both in terms of visual and quantitative statistical measures and all the results are in favor of the claim that the proposed model outperforms the other models considered in this work.



---

*Source: 981932-2014-06-01.xml*