tion applications have no information about the con-
tent of the original (perfectly focused) image. For that
reason, in this work we have selected four state-of-
the-art non-referenced sharpness measures (Ferzli and
Karam, 2009) based on different features, such as im-
age statistics and frequency content. In this selection
we have considered those approaches that can be ap-
plied at a local spatial scale and, hence, can be consid-
ered as information measures in scale-space (Sporring
and Weickert, 1999), since our intention is not to eval-
uate the sharpness/quality of a whole image, but the
sharpness of an image region.
The organization of the rest of this paper fol-
lows. Section 2 presents the four state-of-the-art non-
reference sharpness metrics. In Section 3, we exper-
imentally study the behavior and sensitivity of each
measure to forward and reverse diffusion. A reduced
set of measures showing a good behavior are further
analyzed. Particularly in Section 3.1, for one of the
previous techniques we propose an strategy to detect
the time where the reverse diffusion should be stop
to prevent that the image blows up. Experiments on
real and test images are shown in Section 4. Finally,
conclusions are outlines in Section 5.
2 NON-REFERENCED LOCAL
SHARPNESS MEASURES
This section presents four non-referenced state-of-
the-art techniques that can be applied to measure the
local sharpness of an image region.
2.1 Variance
One of the simplest local sharpness measures is the
variance of an image window. That is,
f
var
(x
0
,y
0
) =
1
N
∑
(x,y)∈Ω
x
0
,y
0
[u(x,y) −u]
2
(1)
where Ω
x
0
,y
0
is a support window around pixel x
0
,y
0
and N is the total number of pixels in the support win-
dow; u(x, y) is the image; and u is the mean value of
the image in the window. The variance is a very sim-
ple and efficient metric that is quite robust to noise
and, intuitively, its value increases as the image gets
sharper since the intensity variation will be higher
than blurred images (Batten, 2000).
2.2 Sum of Modified Laplacians (SML)
Another family of metrics are based on the computa-
tion of second-order derivatives of the image, particu-
larly, the Laplacian. This type of metrics act as a high
pass filter in the frequency domain. Thus, they are
characterized by a good degree of accuracy but they
are very sensitive to noise (since they are based on
the direct computation of image derivatives). Partic-
ularly, here we analyze the Sum of Modified Lapla-
cians (SML), that is given by the following formula
(Aydin and Akgul, 2008):
f
SML
(x
0
,y
0
) =
∑
(x,y)∈Ω
x
0
,y
0
ML(I(x, y)) (2)
where Ω
x
0
,y
0
is a support window around pixel x
0
,y
0
,
and ML(I(x,y)) is the modified Laplacian measure
given by:
ML(I(x, y)) = |−I(x + s, y) + 2I(x,y) −I(x −s, y)|
+ |−I(x, y + s) + 2I(x,y) −I(x, y −s)|
The reason to compute the absolute value of the
Laplacian is to avoid that the horizontal and verti-
cal derivatives may cancel each other. Here, s is
a step variable that is used to set the distance be-
tween the central pixel and the pixels used to com-
pute the second order derivative. Using different s
values may be useful to cope with different sizes of
texture elements. In our experiments, a value of s = 1
was found to produce the best results. The modified
Laplacian measure can be interpreted as an approxi-
mation of the Frobenius matrix norm (Horn and John-
son, 1990) of the Hessian matrix, which introduces
some connections with the detection of salient image
points provided by Lindeberg’s blob detector (Linde-
berg, 1993).
2.3 Frequency Metric
This image sharpness metric has been introduced by
(Shaked and Tastl, 2005). The sharpness is measured
by means of a localized frequency analysis. If u(x,y)
is an image, and |U(ξ
x
,ξ
y
)| is its magnitude spec-
trum, the fractal image model proposed in (Shaked
and Tastl, 2005) states that natural images follow a
fractal behavior given by
|U(ξ
x
,ξ
y
)| =
α
||(ξ
x
,ξ
y
)||
2H+2
2
(3)
where α is a constant, H is the Hurst parameter (Man-
delbrot and Wallis, 1969), and ||·||
2
refers to the Eu-
clidean norm.
According to this model, it is theoretically pos-
sible to reconstruct the original image if the Hurst
parameter H is known. While this is not realistic,
it comes to show that the frequency distribution can
certainly be useful to estimate an image degradation
(for instance, due to blurring). Finally. their proposed
VISAPP2013-InternationalConferenceonComputerVisionTheoryandApplications
70