=Paper= {{Paper |id=Vol-2688/invited1 |storemode=property |title=On Measures of Visual Contrast and Their Use in Image Processing (invited paper) |pdfUrl=https://ceur-ws.org/Vol-2688/invited1.pdf |volume=Vol-2688 |authors=Azeddine Beghdadi |dblpUrl=https://dblp.org/rec/conf/cvcs/Beghdadi20 }} ==On Measures of Visual Contrast and Their Use in Image Processing (invited paper)== https://ceur-ws.org/Vol-2688/invited1.pdf
 On Measures of Visual Contrast and Their Use
             in Image Processing

                                  Azeddine Beghdadi

            L2TI, Institut Galilée, Université Sorbonne Paris Nord, France
                        azeddine.beghdadi@univ-paris13.fr



        Abstract. Visual contrast is one of the most studied notions in the
        field of Visual Neuroscience and Psychophysics. It is a measure asso-
        ciated with a psycho-physical sensation that is not easy to define in
        an objective and unique way. Indeed, the contrast as defined through
        Weber’s famous experiment is associated with the subjective notion of
        just noticeable difference between a stimulus observed against an uni-
        form background. In this paper a critical review of different representa-
        tive measures of visual contrast and their uses in various applications is
        presented. This study provides also some insights on how to define the
        contrast and how to chose the most appropriate measure for developing
        contrast based methods for visual information processing and analysis.
        Perspectives and challenges that remain to be addressed are also dis-
        cussed in light of new trends in visual information processing. Through
        this study, it becomes clear that the concept of contrast and its use
        are highly application-dependent and that there is no universal contrast
        measure. It is also shown that, given the large number of psycho-physical
        parameters involved, it is not easy to define a contrast measure that is
        easy to use in the various methods of image processing and analysis.
        Simplifying contrast models without neglecting the most fundamental
        aspects seems to be the most pragmatic and practicable solution.

        Keywords: Contrast Measures · Just Noticeable Difference (JND) · Vi-
        sual Contrast · Image Processing and Analysis · Perceptual Image com-
        pression · Perceptual watermarking.


1     Introduction

Visual contrast is one of the most studied notions in the field of Psychophysics
and Visual Neuroscience [13, 20]. It is a measure associated with a psycho-
physical sensation that is not easy to define objectively and in a unique way.
Indeed, contrast as defined through Weber’s famous experiment is associated
with the subjective notion of Just Noticeable Difference (JND) between a stim-
ulus observed against an uniformed background. From the study conducted by
    Copyright © 2020 for this paper by its authors. Use permitted under Creative
    Commons License Attribution 4.0 International (CC BY 4.0). Colour and Visual
    Computing Symposium 2020, Gjøvik, Norway, September 16-17, 2020.
2         A. Beghdadi

Hecht [25] it seems that Bouguer (1760) was the first to study the differential
sensibility of the Human Visual System (HVS). The idea of measuring the just
noticeable increment needed to discern one stimulus from another was later stud-
ied by Weber and Fechner [13]. Indeed, the idea of analysing and measuring the
minimum increment of intensity between two stimuli in order to discern them
is at the basis of the very definition of what is commonly known as the Weber-
Fechner contrast [17]. In fact it was Fechner, a few years later, who formalized
Weber’s results in an experimental and theoretical framework. Weber-Fechner
law has led to many experiments and sophisticated models for understanding
the notion of perceptual contrast [25,33]. Since these pioneer works, many visual
contrast measures have been proposed in the literature [16, 24, 31, 35, 39, 41, 42].
    The taking into account of psycho-physical factors such as perceptual con-
trast in a visual information processing and transmission system is essentially
linked to the fact that the observer is the key element in any chain of acquisition,
processing and transmission of visual information. Indeed, the human observer
is often the supreme judge in the evaluation of the different stages of the process-
ing and transmission chain. It is therefore quite natural to think of developing
methods inspired by the mechanisms of the human visual system to incorpo-
rate perceptual criteria that meet the requirements of the observer. It is worth
noticing that among the perceptual aspects of the HVS, visual contrast is one
of the widely investigated psycho-visual aspects in vision research [13,24,25,48].
According to this principle and reasoning several methods of image processing
and analysis based on contrast measure have been developed [3, 48]. Indeed,
contrast plays a prominent role in important applications such as medical imag-
ing [19], image quality enhancement (IQE) [6] and image quality assessment
(IQA) [14, 15], image fusion [9, 44], and other applications such as image quan-
tization and compression [12, 28]. In this article we will review some of these
applications and discuss the most suitable contrast measures in each case. This
is not an exhaustive study of all known methods, but we will limit ourselves
to a few representative studies. The main objectives of this contribution are as
follows:

    – to discuss the fundamental criteria and factors that define visual contrast
      and present a critical review of the most representative contrast measures
      and associated models,
    – provide a brief description and discussion on some applications in the field
      of visual information processing based on perceptual contrast measure
    – provide some insights on how to use and chose the appropriate contrast mea-
      sure in some selected visual information processing and analysis applications,

    The paper is organized as follows: first some historical contrast models are
discussed in Section 2 followed by a classification discussion of some represen-
tative contrast measures in Sections 3, 4, and 5. Section 6 is dedicated to some
selected contrast based applications. Finally the paper ends by concluding re-
marks, challenges and some future directions of research in Section 7.
        On Measures of Visual Contrast and Their Use in Image Processing             3

2    Perceptual Contrast: History and Basic Notions

This section presents a brief historical review of the research on contrast mea-
sures and associated models developed by the scientific research community in
psycho-physics, optics, neuroscience and digital visual information processing.
The notion of visual contrast has been introduced in a clear and well defined
theoretical and experimental framework for the first time by Fechner [17]. Since
this pioneering work based on Weber experiments and model, several studies
have been carried out to enrich existing models with advances on both the theo-
retical and experimental levels [24, 25, 33, 35]. The introduction of the frequency
sensitivity aspect of contrast has been established through psycho-physical ex-
periments [13]. The directional selectivity of HVS has been clearly demonstrated
by Hubel and Wiesel [26], the two Nobel Prize laureates. Other aspects related to
the distance or viewing angle parameter have been introduced explicitly accord-
ing to optical models [33] or implicitly through a multi-resolution representation
of visual contrast [14, 44]. However, the colour aspect was neglected for a long
time in the first experiments. This is due to the fact that the notion of contrast
is much more related to the detection of details and more particularly the con-
tours of objects which is traditionally regarded as an achromatic process [13].
However, it has been shown that the colour aspect also plays a major role in
contrast sensitivity [22]. Other spatio-temporal and 3D aspects could be incor-
porated in the definition of the spatio-temporal contrast. It is worth noticing
that to the best of the author’s knowledge there is no measure of contrast in
an analytical form that integrates all these psycho-physical and geometrical pa-
rameters. Rather simplified and mathematically tractable expressions are often
used to define an objective measure of visual contrast. The essential criteria and
HVS properties that could be taken into account in the visual contrast measure
are given below.

 – sensitivity to the relative change of luminance
 – luminance adaptation phenomena
 – frequency selectivity
 – directional sensitivity
 – multiscale/multiresolution aspects
 – color aspects
 – viewing distance (or viewing angle)
 – temporal aspects (in the case of spatio-temporal visual signals)

Note that it is not easy to define a contrast measure that integrates all these
properties and criteria. Simplifications are often used, keeping only a few prop-
erties and aspects to establish a measure of visual contrast. There are two ways
to define contrast, depending on whether we associate a single value to the entire
image or a value for each pixel or group of pixels. In the first case it is global con-
trast, while in the second case it is local contrast. The contrast measure could be.
computed in the spatial domain, frequency domain, and even multi-resolution
or multi-scale representations.
4      A. Beghdadi

    It should be noted that given the different forms of representation of the
visual signal, the large number of contrast measures proposed in the literature
and the different contexts and fields of application, it is not easy to classify
all the contrast measures. In the following we classify contrast measures into
three categories that is psycho-physic and neuroscience based contrasts, local
structure based contrasts and statistical information based contrasts. Here we
limit ourselves to a few representative contrast measures from each category.


3   Psycho-physic and neuroscience based contrasts
In this first category are some representative contrast measures based on psycho-
physical experiments or models from theoretical and experimental studies in
neuroscience. Figure 1 illustrates the simultaneous contrast phenomena. It also
illustrates one of the most relevant parameters that should be taken into account
in the contrast definition, that is the influence of the surround and background
luminance in the visual perception of stimuli. This figure represents the foveal
image model used in several psycho-physical experiments such as the Moon and
Spencer experiment [33].




              Fig. 1. An example of imultaneous contrast experiment.


    What should be observed in Figure 1 is the effect of the background on
the appearance of the small dark disc and the ring around it. Indeed, as can
be seen, despite having the same luminance and the same gradient in the five
configurations the disc and the ring appear different according to the position
in the non-uniform background. It can be concluded from this example that the
gradient alone cannot account for the visual appearance of the visual signals. The
background against which the object is observed is important. This observation
leads us to question the many contrast measures based solely on the luminance
difference between stimuli, i.e. gradient. This is one of the reasons why it is
important to define contrast as a measure of relative variation in luminance, i.e.
a relative ratio. This is the case with the contrast measure proposed by Weber
and Fechnner [17] described below.
        On Measures of Visual Contrast and Their Use in Image Processing        5

3.1   Weber-Fechner contrast

Weber-Fechner contrast is one of the simplest contrast measures [17]. However,
it only applies to simple scenarios in which a uniform luminance background
L contains an object with an incremental luminance ∆L. The aim of Weber-
Fechner (W-F) contrast is to determine the value of ∆L, referred to as the
JND, which makes the object (target) just visible. The W-F contrast measure
is defined by:
                                          ∆L
                                    CW =                                      (1)
                                            L
One of the most important results of the experiments conducted by Weber and
Fechner is that this ratio remains constant over a fairly wide range of luminance
values. This value is called the Just Noticeable Contrast (JNC) and is of the
order of .02. It should be noted that many sophisticated contrast measures, used
in visual information processing and coding, are in one form or another based
on the Weber-Fechner definition [2, 14, 35, 42, 44, 50].


3.2   Michelson contrast

The Michelson contrast was first introduced in a purely physical context and
concerns the measurement of the visibility of interference fringes produced by
thin films [32]. However, it has been widely used in psycho-physical experiments
and more particularly the study of the frequency sensitivity of the HVS to per-
ceptual contrast [13, 45]. The Michelson contrast is defined as follows:

                                     Lmax − Lmin
                              CM =               ,                            (2)
                                     Lmax + Lmin
where Lmin and Lmax correspond to the minimum and maximum luminance
values in the optical image, respectively. Although this contrast has been used
extensively by the scientific community of vision research, it has several limita-
tions. Indeed, its use in the case of natural images can lead to over- or under-
estimation of contrast and in particular in the case of images contaminated by
impulse-type noise. It is also the case of images containing some singularities,
or isolated points even if they are perceptually invisible. Moreover, it does not
take into account the influence of the background and in particular the lumn-
inance adaptation phenomenon. It also does not consider the frequency aspect
although the stimulus signal shape is designed from a sinsusoidal function. It is
nevertheless surprising that such a simple contrast model has been widely used
by the vision research community for almost a century..


3.3   Moon-Spencer contrast

The main idea behind Moon and Spencer’s [33] model is to apply Holladay’s prin-
ciple [33] that any non-uniform background may be replaced by another uniform
luminance that produces the same perceptual effect. This leads to the definition
6      A. Beghdadi

of luminance adaptation that can be calculated according to a simplified model.
Based on this principle, Moon and Spencer proposed a simple enough model to
express the adaptive luminance, which is given below.

                             LA = αS LS + αB LB ,                             (3)

Where LS and LB are the luminance of the surround (immediate neighborhood)
and that of the background (or far surrounding), respectively. The two weight-
ing parameters, αS and αB are set experimentally to the values .923 and .077,
respectively. Moon and Spencer define the minimum perceptible contrast as:
                                    √ 
                          CW A + LA 2          if LA ≥ LS
                           LS
                  Cmin = C          q
                                          2
                                            2              .              (4)
                          W A + LS             if LA < L S
                           LS           LA


Where CW corresponds to the Weber-Fechner JNC contrast, and A is a constant
determined experimentally from psycho-visual tests and Hecht’s law [25] which
is equal to 0.808. Note that this contrast is more interesting in the sense that
it corresponds more or less to realistic configurations. It has been successfully
used and adapted to digital images in various applications [6, 27, 28, 34].


3.4   Lillesaeter contrast

Noting the asymmetry in Weber’s definition of contrast, Lillesaeter proposed two
measures of contrast [30]. In the first one only luminance is taken into account
while in the second definition he proposes to include the geometry and shapes
of the objects observed in the image. Indeed, it is observed that the negative
and positive Weber-Fechner contrasts with the same absolute increment are not
perceived equally. The Lillesaeter contrast is then defined as

                                            LO
                                  C = log                                     (5)
                                            LB

where LO and LB correspond to the average luminance of the object and the
background, respectively. Note that when LO and LB are very close to each
other the Lillesaeter contrast is equivalent to that of the Weber-Fechner contrast
measure
                     LO − LB
                C=              ≈ log LO − log LB (if |C| << 1).               (6)
                        LB
The second definition of Lillesaeter contrast incorporates the object contour
geometry as perceived by human. The idea of taking into account the geometry
of perceived objects is relevant but impractical in the evaluation of contrast in
digital images. Indeed, it leads to the computation of curvilinear integral which
requires the exact knowledge of the object contours. This leads inevitably to an
ill-posed problem which is image segmentation.
        On Measures of Visual Contrast and Their Use in Image Processing             7

3.5   DOG based contrast

Based on models describing the retinal ganglion cells and Lateral Geniculate
Nucleus (LGN) responses [29] to visual stimuli on the other, Tadmor and Tol-
hurst [42] proposed three measures of local contrast. The principle is based on
the use of linear filtering by two isotropic Gaussian Impulse Responses (IR) of
different sized kernels. The bandpass behaviour of the HVS is then modelled
through the Difference Of Gaussian (DOG) model [29]. The two Gaussian IRs
are associated with two regions, a center zone Ωc and a surround zone Ωs , to
mimic the receptive fields of the ON and OFF cells [20, 29].
    The two responses are given by:

                              Iσc (i, j) = (I ∗ hσc )(i, j),                       (7)

and
                              Iσs (i, j) = (I ∗ hσs )(i, j)                        (8)

where I is the input image signal, hσc and hσc are the two Gaussian IR associated
with the center and surround zones .
   The convolutions are performed in the sliding windows Ωc (i, j) and Ωs (i, j) of
odd size [−3σc , +3σc ]x[−3σc , +3σc ] and [−3σs , +3σs ]x[−3σs , +3σs ], respectively.
   Three local contrasts are then defined as follows:
                                          Iσc (i, j) − Iσs (i, j)
                         C1DOG (i, j) =                           ,                (9)
                                                 Iσc (i, j)

                                          Iσc (i, j) − Iσs (i, j)
                         C2DOG (i, j) =                           ,               (10)
                                                 Iσs (i, j)

                                          Iσc (i, j) − Iσs (i, j)
                         C3DOG (i, j) =                           ,               (11)
                                          Iσc (i, j) + Iσs (i, j)
The global contrast is derived by averaging the local contrasts. It can be noticed
that these contrasts do not take into account the directional and frequency se-
lectivity nor the colorfulness aspects. Note that this contrast expressed as a ratio
between a differential signal component and a low-pass component is somehow
inspired by Weber-Fechner’s simple model and corresponds well to the notion of
visual contrast.


4     Local structure based contrast measures

In this section, we introduce and discuss some representative contrast measures
that could be used in digital image processing and analysis applications. These
are essentially measures that explicitly or implicitly incorporate some features
of the SVH.
8       A. Beghdadi

4.1   Edginess based contrast measure

It is well established that one of the primitives of the image signal that is most
related to contrast and therefore to the visibility of detail is contour informa-
tion. Indeed, the most representative contours of the image signal correspond to
the spatial frequencies where the CSF reaches its maximum values [13, 29, 45].
Inspired by the contrast defined by Gordon and Rangayan in their contrast en-
hancement method [19], Beghdadi and Le Negrate introduced a new measure of
contrast incorporating edginess information [4]. This local contrast measure is
computed using a sliding window of odd size . For each Wij window, the mean
edge gray-level at the center pixel (i, j), is computed as
                                      P
                                              Φ(∆kl )f (k, l)
                                    (k,l)∈Wij
                        E(i, j) =           P                      .           (12)
                                                    Φ(∆kl )
                                        (k,l)∈Wij


In (12), f (k, l) corresponds to the gray level at the pixel (k, l) and Φ(∆kl ) is
an increasing monotonic function of the gradient operator at (k, l). A simple
function would be ∆nkl , with n > 0. The local contrast is expressed as:

                                       |E(i, j) − f (i, j)|
                           C(i, j) =                        .                  (13)
                                        E(i, j) + f (i, j)

Note that this contrast measure does not take into account the frequency selec-
tivity nor the directional selectivity aspects. Furthermore, when used in CE
it may introduce halo effects around the edges which is the result of over-
enhancement [5, 38]


4.2   Multiresolution contrasts

Toet was the first to propose a contrast measures taking into account the multi-
resolution aspect for image fusion [44]. It is based on the Burt and Adelson
pyramid decomposition scheme [10]. The contrast is expressed as

                                       gk (i, j) − gk−1 (i, j)
                         Ck (i, j) =                           ,               (14)
                                               gk (i, j)

Where the components gk (i, j) and gk−1 (i, j) are the gray-level of pixel (i, j) in
the Gaussian pyramid at the k th and (k − 1)th levels, respectively. Note that this
expression of local contrast is also inspired by Weber-Frechner’s intuitive defini-
tion. Indeed, the numerator is nothing more than a differential signal (intensity
increment) and the denominator is the signal against which the increment is
measured (reference signal). Here is also,the directional sensitivity and the col-
orfulness aspects are note taken into account. This contrast measure has inspired
several works and led to various interesting applications such as image fusion [44],
contrast enhancement [40] and image distortion prediction [14].
        On Measures of Visual Contrast and Their Use in Image Processing           9

4.3   Bandlimited contrast

By exploiting the results of psycho-physical experiments on the frequency sen-
sitivity of the HVS to contrast and in particular its behaviour as a bandpass
filter, Peli introduced the notion of band-limited contrast [35]. The image signal
is analyzed using a bank of cos-log type isotropic band-pass filters to extract
the different components describing the signal at different frequency bands. The
image component captured by the k th channel is given by

                              gk (i, j) = (f ∗ hk )(i, j).                      (15)

Where hk , is the impulse response corresponding to the k th band-pass filter and
gk the associated filtered component. For each pixel (i, j) in the k th component,
the contrast is expressed as:

                                               gk (i, j)
                                 Ck (i, j) =                                    (16)
                                               bk (i, j)

where the bk is given by:
                                            k−1
                                            X
                              bk (i, j) =         gm (i, j).                    (17)
                                            m=0

    Here, too, it can be assumed that this measure is somehow inspired by the
intuitive idea of Weber and Fechner’s model. Indeed, the numerator is the differ-
ential signal , i.e. pass-band signal, and the denominator contains the sum of all
the lower frequency components, i.e. baseband signal. This ratio does measure a
relative change in signal amplitude just as in the Weber-Fechner contrast model.
    Like Toet’s contrast, Peli’s contrast can be used for complex natural images,
unlike other contrasts such as Weber or Michelson. However, the lack of direc-
tional selectivity and colorfulness aspects of Peli’s contrast limits its use in some
real-world applications where these aspects play prominent roles.


4.4   Daly contrast

The contrast model proposed by Daly [14] is essentially based on the cortex trans-
form introduced by Watson [1]. It should be noted that in this contrast model
both frequency selectivity and directional selectivity are taken into account by
means of two families of linear filters. The input signal f is first analysed by
means of two cascades of isotropic band-pass linear filters called ”dom” and
”fan” corresponding to the frequency selectivity for the former and the direc-
tional selectivity for the latter [14]. The filtered versions of the input signal f
are given by:
                              gkl (i, j) = (hkl ∗ f )(i, j)                    (18)
where k and l are the dom and fan filter indices as defined in [1]. Daly contrast
is then defined by
10      A. Beghdadi


                                          gkl (i, j) − g kl (i, j)
                           Ckl (i, j) =                            .              (19)
                                                 g kl (i, j)
    Where g kl is the mean of the band (k, l). Note that this measure is unstable
because of the denominator, which tends to zero. Daly proposes two solutions
to overcome this problem. He then introduces two contrast measures where the
denominator is replaced in one case by the baseband signal mean and in the other
by the baseband signal calculated at each pixel. These two modified contrast
measures are given by:

                                          gkl (i, j) − g kl (i, j)
                           Ckl (i, j) =                            .              (20)
                                                    gK

and
                                          gkl (i, j) − g kl (i, j)
                           Ckl (i, j) =                            .              (21)
                                                 g K (i, j)
Note that Daly contrast model does not incorporate the colorfulness aspect. This
contrast measure has been used successfully in the design of image distortion
prediction models [14].


4.5   Isotropic contrast

The consideration of the multi-scale aspect in the visual contrast measure is
to some extent related to the characteristics of the HVS [26]. One of the first
contrast measures based on wavelet analysis was introduced by Winkler and
Vandergheynst [50]. It is surprising that most of the contrast measures proposed
so far did not explicitly take into account the multiscale aspect. However, it
should be noted that Toet’s proposed measure is somehow quite close in that it
introduces the multi-resolution aspect. The main idea of the measure proposed
by Winkler and Vandergheynst was to overcome the limitations of the contrast
proposed by Peli. They proposed a contrast measure using a directional wavelet
decomposition based on a translation invariant multiresolution representation
using 2-D analytical filters. By combining the different analytic oriented filter
responses they derived the isotropic contrast expressed as follows.
                                      r P
                                        2 |gkl (i, j)|2
                                                 l
                            Ck (i, j) =                          ,                (22)
                                                 gk (i, j)
     where gkl (i, j) is the gray-level at pixel (i, j) in the band-limited directional
filtered image obtained by filtering the input signal f by the directional wavelet
at resolution k and direction l. Similarly to Peli’s contrast, the denominator cor-
responds to the baseband signal, i.e. the filtered signal with the scaling function
at scale k. It has been demonstrated that in contrast to the Peli’s model, this
new contrast gives a flat response to sinusoidal patterns [50]. However, it is im-
portant to note that although in its design this contrast uses directional filters,
        On Measures of Visual Contrast and Their Use in Image Processing       11

it provides an isotropic contrast measure. This could be beneficial for certain
applications where directionality is not important, such as in the case of digital
watermarking [46].

4.6   Directional bandlimited contrast
Based on the work of Peli and Daly, Dauphin et al. [16] proposed a contrast
where directional selectivity is taken into account in the final contrast measure.
Indeed, other contrasts such as those of Daly and Winkler-Vandergheynst inte-
grate directional selectivity in the analysis of signal components but the final
contrast is rather isotropic. Whereas, in the model defined in [16] the final con-
trast measure is anisotropic in the sense that the final response emphasizes the
most directional salient signal components in the signal. A non-linear operation
of type max is thus used for the calculation of the local contrast. The image sig-
nal is analysed using a multichannel Gabor decomposition. The local directional
bandlimited contrast is computed as
                                         max (|gkl (i, j)|)
                                           l
                           Ck (i, j) =                                        (23)
                                               gk (i, j)
Where gk (i, j) is the gray-level associated to the frequency sub-band ρk and to
one of the four directions (0, π/4, π/2, 3π/4) represented by l. The normaliza-
tion term gk (i, j) represents the total energy of the background below the band
(k) which is obtained by filtering the original image by a Gaussian filter with a
standard deviation
                                    σk = 0.75ρk                              (24)
This local bandlimited directional contrast does not incorporate the colorfulness
aspect. It has been compared to Peli’s contrast and has been proven more efficient
and less complex than the wavelet-based contrast proposed in [50].

4.7   Multiscale color contrast
The RAMMG contrast measure proposed by Rizzi et al. [39] is one of the few
measures of contrast that incorporates both multi-resolution and colour aspects.
The image signal is decomposed using a pyramidal scheme in the CIELAB colour
space. At each level of resolution each pixel is associated with a local contrast
defined as the response of a pseudo-Laplacian computed by convolving the input
signal I by the mask S given by:
                                         √      √ 
                                            2      2
                                               1
                                    1     2      2 
                            S=        √  1 0 1 .                           (25)
                                 4 + 2 2 √2 √2
                                           2   1  2

The local contrast at the k th level of the pyramid is computed as follows.

                            Ck (i, j) = (Dk ∗ S)(i, j),                       (26)
12      A. Beghdadi

where Dk (i, j) is the absolute difference of the luminance between the current
pixel (i, j) and the central pixel at the k th level of resolution. The RAMMG global
contrast is obtained by averaging all local contrasts across all the different levels
of the pyramidal decomposition.

                                           X Xl XK l  W   H
                                    1
                 C RAM M G =                         Ck (i, j).                 (27)
                               Wl × Hl × K   i=1 j=1
                                                k=1


where K is the total number of decomposition levels and Wk and Hk represent
the width and height of the image at the k th level respectively. A very similar
global contrast measure called Gobal Contrast Factor (GCF) has been proposed
in [31]. GCF contrast has been found somehow consistent with subjective ranking
of a relatively wide range of natural images with varying contrast. However, it
should be noted that both contrast RAMMG et GFC are not relative measures of
the variation in the energy of the image signal and therefore cannot be included
in the family of contrasts conforming to the notion of contrast as defined by
Weber-Fechner.


5     Statistical features based contrasts

Very few contrasts measures based on statistical information have been intro-
duced in the literature. These measures are often related to the pixel values
distribution, such as the grey-level histogram or the 2D distribution computed
from the grey-level cooccurrence matrix (GLCM). Here we limit ourselves to
three contrast metrics based on some simple statistical features of pixel values.


5.1    Texture contrast measure

Haralick was the first who introduced the idea of using some statistical invariant
features for texture analysis. The set of these spatial descriptors introduced
in [23] are based on the GLCM computed from the digital image. Among the
Haralick’s texture descriptors a global contrast is defined. It is computed as
follows.

                                   K−1
                                   X K−1
                                       X
                            CH =              )(i − j)2 pij                     (28)
                                    i=0 j=0

Where i and j are the grey-levels of adjacent pixels in a defined neighbourhood
and pij is the joint mass probability function computed from the GLCM. Al-
though this measure has always been identified as a contrast, it does not meet
the basic criteria to be truly considered as a measure of contrast in the usual
sense and in line with psycho-physical experiences and the notion of contrast
that has been well established since the 19th century.
        On Measures of Visual Contrast and Their Use in Image Processing           13

5.2   Mutual information based contrast measure
Another way to exploit inter-pixel correlation, directly related to image contrast,
is to consider the measure of mutual information extracted from the GLCM. A
new global contrast measure based on mutual information has thus been intro-
duced for the first time to quantify the side effects such as saturation or halo
effect that could result from contrast enhancement [38]. This contrast measure
is defined by:
                           K−1
                                                                 !
                           X K−1X                    pxy (i, j)
                   CM I =           pxy (i, j) log2               ,           (29)
                            i=0 j=0
                                                    px (i)py (j)
where pxy is the joint probability mass function of the gray-level, whereas px and
py represent the marginal probabilities computed from the GLCM. While this
contrast is simple to compute, however, it does not provide information directly
related to visual contrast as it is purely based on statistical analysis of the signal
values distribution.

5.3   Root Mean Square contrast
The Root Mean Square (RMS) of the luminance in natural images has been
considered by Bex and Makhous [7] as a potential contrast measure in their
study on human observer sensitivity to contrast. As noticed by these authors,
this measure when divided by the average mean luminance of the image is a
good predictor of the relative contrast.
    Another version of this RMS based contrast has been proposed by Frasor
and Geisler [18] to make it more suitable to natural images . Local contrast
is measured in different randomly selected patches in the image. The contrast
associated with a given patch is calculated as follows:
                                 v
                                 u       N
                                 u 1 X         (Li − L)2
                       CRM S = t            wi           .               (30)
                                    wN i=1         L2

                                           N
                                           X
                                    wN =         wi                              (31)
                                           i=1

where N is the total number of pixels in the patch, Li is the luminance of the ith
pixel, L the patch luminance and wi is a windowed isotropic weighting function
given by:
                                               π
                               wi = 0.5 × cos ri                                (32)
                                               p
             p
where, ri = (xi − xc )2 + (yi − yc )2 , p is the radius of the patches , (xi , yi ) is
the position of the ith pixel within the patch, and (xc , yc ) is the center of the
patch. The patch luminance is given by
                                          N
                                       1 X
                                 L=          w i Li                              (33)
                                      wN i=1
14      A. Beghdadi

    It could be noticed that none of these statistical information based contrast
measures take into account the spatial frequency and directional content of the
visual signal. Nor do they incorporate other important aspects such as the lumi-
nance influence of near and far surrounds, the viewing distance and chromatic
aspect. Table 1 summarizes the key features of these representative local and
global contrast measures.

Concluding remarks From this brief review, we conclude that it is uneasy
to find a simple contrast measure that incorporates all the relevant perceptual
aspects related to the contrast notion. Therefore, most often simple measures
derived from psycho-physical experiments are used and adapted to digital im-
ages to solve a number of real problems such as image and video processing,
analysis and compression. Despite the enormous amount of work dedicated to
visual contrast, it is still not easy to evaluate and compare the different contrast
measures objectively. However, a few studies limited to subjective evaluation
have been carried out in the literature [41].
    Moreover, through this critical analysis on the existing contrast measures, it
becomes important to answer some relevant questions. One important question
is: what are the most relevant visual signal characteristics that the contrast mea-
sure should capture? To the best of author’s knowledge none of the published
work addressed this issue properly. However, some interesting studies have been
dedicated to address this critical question [36]. Peli conducted a thorough ex-
perimental study in 1997 on how to define the contrast measure [36]. To be
consistent with the findings on the human visual perception the study provides
guidelines confirming that computational contrast metrics should take into ac-
count multiscale aspects [26].
    It is also worth noticing that there is no clear and objective criteria nor
ground truth data to exploit in comparing the proposed contrast measures. Al-
though, some attempts have been made in the study of different measures of
contrast in digital images [41].


6     A Brief Overview of the use of Contrast in VIPA
An important issue to consider is how the contrast measure is defined and used
in developing various methods for Visual Information Processing and Analysis
(VIPA). We limit ourselves here to a few applications where the notion of con-
trast plays a predominant role. It should be noted that the choice of contrast
measure depends on the application. It is not always easy to make an appropriate
choice and pragmatic solutions are often used, based mainly on feedback.

6.1   Image Quality Assessment and Enhancement
Most objective HVS-based image quality methods explicitly or implicitly incor-
porate contrast measure [11]. Peli’s contrast and its variants have been success-
fully incorporated in the design of IQA measures [15] and distortion predictor
        On Measures of Visual Contrast and Their Use in Image Processing             15

            Table 1. Summary of some representative contrast measures.


Contrast                    Year     Key features
Weber-Fechner contrast [17] 1860     Defined on simple images (an object of uniform
                                     intensity in a uniform background). Could not be
                                     applied to natural and complex images
Michelson contrast [32]      1927    used for sinusoidal signals, not adapted for natural
                                     images. Does not integrate any HVS properties.
Moon-Spencer contrast [33] 1944      considers object embedded in a non-uniform lu-
                                     minance background. Can be applied to digital
                                     images of natural scenes with a few adjustments
                                     and adaptations.
Haralick contrast [23]       1973    capture the average local variations and spatial
                                     dependence of the pixels computed from gray-level
                                     co-occurrence matrix
Edginess based contrast [4] 1986     uses local edginess information in an image and
                                     quantifies local sharpness of the contour by mea-
                                     suring the visibility of the salient features by using
                                     a sliding window.
Multiresolution contrast [44] 1989   a spatial multi-resolution contrast based on Gaus-
                                     sian pyramid decomposition.
Bandlimited contrast [35]    1990    An image is decomposed into several channels
                                     using a bank of cosine-log bandpass filters, well
                                     suited for complex images and often used in the
                                     computation of the quality of encoded images.
Daly contrast [14]           1992    It is band-limited contrast based on the cortex
                                     transform and used by Daly in the Visible Differ-
                                     ence Predictor (VDP) model [48].
Lillesaeter contrast [30]    1993    not easy to use in practice since it requires be-
                                     forehand knowledge of the object contours and
                                     computation of the curvilinear integral along the
                                     boundaries of the objects contained in the ob-
                                     served image.
Winkler-Vandergheynst        1999    uses non-seperable directional filters. It has the
contrast [50]                        advantage of giving a flat response instead of an
                                     oscillating response to sinusoid gratings.
DOG based contrast [42]      2000    Based on Lateral Geniculate Nucleus responses
                                     model [20]. It uses Difference Of Gaussian (DOG)
                                     to model the bandlimited responses.
Directional contrast [16]    2003    It is local directional bandlimited contrast com-
                                     puted from multichannel Gabor decomposition
                                     and a nonlinear operation.
RAMMG contrast [39]          2004    Multilevel analysis (pyramid representation), 8-
                                     neighborhood pyramid subsampling of the image
                                     to various levels in the CIELAB.
RMS Contrast [18]            2006    a based on the RMS contrast introduced by
                                     Bex and Makous [7] using a windowed isotropic
                                     weighting function. It is a single scale contrast
                                     based on only pixel values.
MI based contrast [38]       2015    based on the mutual information computed from
                                     the gray-level co-occurrence matrix.
16     A. Beghdadi

models [14]. Another application where the concept of contrast is important con-
cerns the improvement of image quality and in particular Contrast Enhancement
(CE). The CE methods could be roughly classified into two categories, namely
direct methods and indirect methods [4]. The main idea of the direct methods is
to estimate the local contrast and to amplify it by means of a monotonic trans-
formation and to deduce the intensity of the pixel corresponding to this new
contrast. The contrast defined by Beghdadi and Le Négrate [4] and its variants
have been successfully used in many CE methods [3]. The edginess based contrast
measure has been also extended to stereo images by incorporating the depth in-
formation [21]. Another direct method based on the bandlimited contrast [35]
proposed by Tan et al. [43] which operates in the discrete cosine transform do-
main. Unfortunately most CE methods suffer from some side effects and there
is no unified framework to control these effects. Indeed, this pre-processing that
aims at amplifying the visibility of details by increasing gradient and sharpness
may introduce some artefacts such noise amplification, saturation and halo effect.
It is therefore useful to quantify these undesirable side effects. Many objective
measures for Contrast Enhancement Evaluation (CEE) based on the local or
global contrast have been studied in [37]. A critical study of CEE metrics has
revealed that the mutual information contrast measure is the most promising in
terms of simplicity and efficacy [5]. However, this metric has limitations and it
is in our interest to develop other measures that integrate the multi-scale and
multi-directional aspects. From this point of view the three contrast measures
defined in [8, 16, 50] are good candidates.


6.2   Visual Data Protection and Compression

The protection and coding of visual data are two classic problems of very active
research. Here, it is more precisely about image watermarking using perceptual
approaches [3]. The aim is to insert a watermark that is both invisible and re-
sistant to various attacks. This is then a very difficult problem where one tries
to achieve the best compromise between two antagonistic criteria that are ro-
bustness and transparency. Indeed, robustness requires putting more energy into
the watermark, which inevitably makes it visible and therefore breaks the trans-
parency. The visibility of the watermark is very much linked to the notion of JND
defined in the contrast measure. In this type of application where one seeks to
insert in a robust and transparent way the watermark the multi-scale or pyrami-
dal approach through JND measures related to the contrasts defined in [44,50] is
the most promising solution as shown by the studies in [34,46]. The other appli-
cation where contrast plays an important role is image compression with quality
control. This involves using contrast as a measure of the visibility of distortion
and artifacts inherent in lossy compression methods. The JND measure, again
very much related to contrast, and the visual masking phenomenon are the most
important parameters to be considered in the quantization and coefficient selec-
tion scheme in the transformed domain. The idea of exploiting contrast in the
quantifization of the image signal dates back to the work of Kretz [28] based on
        On Measures of Visual Contrast and Their Use in Image Processing        17

the contrast of Moon and Spencer [33]. Since this pioneering work other more so-
phisticated contrast measurements have been successfully introduced into image
compression and video coding models [12, 47]. Another interesting perceptual
coding scheme proposed in [51] based on the Watson-Solomon Contrast Gain
Control (CGC) [49] model has been proposed for High Efficiency Video Coding
(HEVC). It is worth noticing that most of visually lossless coding methods in the
literature exploit the contrast measure in an explicit or implicit manner [3, 48].
The only contrast model for perceptual coding and quantization that seems to
be complete and efficient is the one introduced recently in [51]. Therefore, we
recommend to use this model for perceptual coding.


6.3   Image fusion

Image fusion is becoming an active field of research and especially with the reviv-
ing of artificial intelligence based approaches. The use of perceptual information,
and particularly contrast, in visual data fusion schemes seems to be the most
promising approach in various applications [8, 40]. There are several ways to
merge information, depending on the application and the data analysis method
used. A first approach is to exploit the multi-scale representation of visual in-
formation in the development of the fusion scheme. Perceptual contrast is one
of the signal features that could be used in the design of the fusion scheme. The
idea behind the use of perceptual contrast is to exploit the most relevant per-
ceptual information that is contained in the contrast map. The strategy consists
then in using the contrast measure in the weighting function used in the fusion
scheme. One of the attractive image fusion approach based on this idea is to
use the directional wavelet based contrast as done in [8]. This strategy could be
used not only in multimodal medical imaging but also in various applications
such multi-modal video-surveillance, multi-focus based computing photography
and hyperspectral imaging, to name a few. It is also the case of the contrast
enhancement method based on a perceptual fusion scheme proposed in [40].


7     Conclusions and challenges

Through this panoramic and chronological study on visual contrast, and the
underlying models and experimental studies, it becomes clear that it is not easy
to express a contrast measure where all the relevant psycho-physical factors and
parameters related to the notion of visual contrast are taken into account. Note
also that the notion of contrast and its use are very application-dependent and
the best way to exploit a visual contrast model to solve real problems is to
simplify it while keeping the most fundamental aspects.
    Furthermore, it is difficult to classify the existing definitions and measures
related to the notion of contrast. This is mainly due to the various forms and
representations of visual and optical signals. Indeed, development of imaging
technologies has led to various image modalities. Therefore, it is now necessary
18      A. Beghdadi

to rethink and define the concept of contrast according to the modality of the
visual signal under consideration.
    It should also be noted that the absence of a universal definition of visual
contrast has opened the door to formulations and extensions of this concept
to other contexts where one seeks only to quantify the difference between two
stimuli or elements of the signal. This has led, for example, to the consideration
of signal gradient measure as contrast in several studies. A unifying framework
with clear criteria for defining the contrast will certainly help in avoiding such
confusions and mistakes and allows to progress properly in exploiting vision
research results in developing efficient contrast based VIPA methods.
    We can also see through this study that the measure of perceptual contrast
in the case of colour images is not sufficiently studied. According to the author’s
knowledge at the present time there is no well established definition or measure
of chromatic contrast that is recognized by the scientific community in the field
of vision research or digital image processing.
    The temporal aspect is also important to introduce in the contrast measure.
Another aspect that could be taken into account is the inter-channel interactions
in the definition of multi-scale contrast. With the renewed interest in artificial
intelligence approaches for solving complex problems and in particular feature-
based learning approaches, the contrast may play a key role in the design of
perceptual loss functions in convolutional neural network architectures.


References
 1. Andrew B, W.: The cortex transform: Rapid computation of simulated neural
    images. Computer Vision Graphics and Image Processing 39(3), 311–327 (1987)
 2. Beghdadi, A., Dauphin, G., Bouzerdoum, A.: Image analysis using local band direc-
    tional contrast. In: Proc of the International Symposium on Intelligent Multimedia,
    Video and Speech Processing, ISIMP’04 (2004)
 3. Beghdadi, A., Larabi, M.C., Bouzerdoum, A., Iftekharuddin, K.M.: A survey of
    perceptual image processing methods. Signal Processing: Image Communication
    28(8), 811–831 (2013)
 4. Beghdadi, A., Le Negrate, A.: Contrast enhancement technique based on sc detec-
    tion of edges. Computer Vision, Graphics, and Image Processing 46(2), 162–174
    (1989)
 5. Beghdadi, A., Qureshi, M.A., Deriche, M.: A critical look to some contrast en-
    hancement evaluation measures. In: 2015 Colour and Visual Computing Sympo-
    sium (CVCS). pp. 1–6. IEEE (2015)
 6. Belkacem-Boussaid, K., Beghdadi, A.: A new image smoothing method
    based on a simple model of spatial processing in the early stages of hu-
    man vision. IEEE Transactions on Image Processing 9(2), 220–226 (2000).
    https://doi.org/10.1109/83.821735
 7. Bex, P.J., Makous, W.: Spatial frequency, phase, and the contrast of natural im-
    ages. JOSA A 19(6), 1096–1106 (2002)
 8. Bhatnagar, G., Raman, B.: A new image fusion technique based on directive con-
    trast. ELCVIA: electronic letters on computer vision and image analysis 8(2),
    18–38 (2009)
         On Measures of Visual Contrast and Their Use in Image Processing            19

 9. Bhatnagar, G., Wu, Q.J., Liu, Z.: Directive contrast based multimodal medical
    image fusion in nsct domain. IEEE transactions on multimedia 15(5), 1014–1024
    (2013)
10. Burt, P., Adelson, E.: The laplacian pyramid as a compact image code. IEEE
    Transactions on communications 31(4), 532–540 (1983)
11. Chandler, D.M.: Seven challenges in image quality assessment: past, present, and
    future research. ISRN Signal Processing 2013 (2013)
12. Chandler, D.M., Hemami, S.S.: Dynamic contrast-based quantization for lossy
    wavelet image compression. IEEE Transactions on Image Processing 14(4), 397–
    410 (2005)
13. Cornsweet, T.: Visual Perception. Academic Press (1970)
14. Daly, S.J.: Visible differences predictor: an algorithm for the assessment of image
    fidelity. In: Human Vision, Visual Processing, and Digital Display III. vol. 1666,
    pp. 2–16. International Society for Optics and Photonics (1992)
15. Damera-Venkata, N., Kite, T.D., Geisler, W.S., Evans, B.L., Bovik, A.C.: Image
    quality assessment based on a degradation model. IEEE transactions on image
    processing 9(4), 636–650 (2000)
16. Dauphin, G., Beghdadi, A., de Lesegno, P.V.: A local directional bandlimited con-
    trast. In: Seventh International Symposium on Signal Processing and Its Applica-
    tions. vol. 2, pp. 197–200. IEEE (2003)
17. Fechner, G.: Elemente der psychophysik (leipzig: Breitkopf & hartel). In: English
    translation of Vol. 1 by HE Adler 1966. Holt, Rinehart, and Winston New York
    (1860)
18. Frazor, R.A., Geisler, W.S.: Local luminance and contrast in natural images. Vision
    research 46(10), 1585–1598 (2006)
19. Gordon, R., Rangayyan, R.M.: Feature enhancement of film mammograms using
    fixed and adaptive neighborhoods. Applied optics 23(4), 560–564 (1984)
20. Grossberg, S., Mingolla, E., Todorovic, D.: A neural network architecture for preat-
    tentive vision. IEEE Transactions on Biomedical Engineering 36(1), 65–84 (1989)
21. Hachicha, W., Beghdadi, A., Cheikh, F.A.: Combining depth information and local
    edge detection for stereo image enhancement. In: 2012 Proceedings of the 20th
    European Signal Processing Conference (EUSIPCO). pp. 250–254. IEEE (2012)
22. Hansen, T., Gegenfurtner, K.R.: Color contributes to object-contour perception in
    natural scenes. Journal of Vision 17(3), 14–14 (2017)
23. Haralick, R.M., Shanmugam, K., Dinstein, I.H.: Textural features for image classi-
    fication. IEEE Transactions on systems, man, and cybernetics (6), 610–621 (1973)
24. Haun, A.M., Peli, E.: Measuring the perceived contrast of natural images. In: SID
    Symposium Digest of Technical Papers. vol. 42, pp. 302–304. Wiley Online Library
    (2011)
25. Hecht, S.: The visual discrimination of intensity and the weber-fechner law. The
    Journal of general physiology 7(2), 235–267 (1924)
26. Hubel, D.H., Wiesel, T.N.: Brain mechanisms of vision. WH Freeman (1979)
27. Iordache, R., Beghdadi, A., de Lesegno, P.V.: Pyramidal perceptual filtering using
    moon and spencer contrast. In: Image Processing, 2001. Proceedings. 2001 Inter-
    national Conference on. vol. 3, pp. 146–149. IEEE (2001)
28. Kretz, F.: Subjectively optimal quantization of pictures. IEEE Transactions on
    Communications 23(11), 1288–1292 (1975)
29. Levine, M.D.: Vision in man and machine. McGraw-Hill College (1985)
30. Lillesaeter, O.: Complex contrast, a definition for structured targets and back-
    grounds. JOSA A 10(12), 2453–2457 (1993)
20      A. Beghdadi

31. Matkovic, K., Neumann, L., Neumann, A., Psik, T., Purgathofer, W.: Global con-
    trast factor-a new approach to image contrast. Computational Aesthetics 2005,
    159–168 (2005)
32. Michelson, A.: Studies in Optics. The Univ. of Chicago Science Series, University
    Press (1927), https://books.google.fr/books?id=FXazQgAACAAJ
33. Moon, P., Spencer, D.E.: Visual data applied to lighting design. JOSA 34(10),
    605–617 (1944)
34. Nguyen, P.B., Beghdadi, A., Luong, M.: Perceptual watermarking using a new
    just-noticeable-difference model. Signal Processing: Image Communication 28(10),
    1506–1525 (2013)
35. Peli, E.: Contrast in complex images. JOSA A 7(10), 2032–2040 (1990)
36. Peli, E.: In search of a contrast metric: Matching the perceived contrast of gabor
    patches at different phases and bandwidths. Vision Research 37(23), 3217–3224
    (1997)
37. Qureshi, M.A., Beghdadi, A., Deriche, M.: Towards the design of a consistent image
    contrast enhancement evaluation measure. Signal Processing: Image Communica-
    tion 58, 212–227 (2017)
38. Qureshi, M.A., Deriche, M.A., Beghdadi, A., Mohandes, M.: An information based
    framework for performance evaluation of image enhancement methods. In: 2015
    International Conference on Image Processing Theory, Tools and Applications,
    IPTA 2015, Orleans, France, November 10-13, 2015. pp. 519–523 (2015)
39. Rizzi, A., Algeri, T., Medeghini, G., Marini, D.: A proposal for contrast measure
    in digital images. In: Second European Conference on Color in Graphics, Imaging,
    and Vision (CGIV). pp. 187–192. Society for Imaging Science and Technology,
    Aachen (2004)
40. Saleem, A., Beghdadi, A., Boashash, B.: A distortion-free contrast enhancement
    technique based on a perceptual fusion scheme. Neurocomputing 226, 161–167
    (2017)
41. Simone, G., Pedersen, M., Hardeberg, J.Y.: Measuring perceptual contrast in dig-
    ital images. Journal of Visual Communication and Image Representation 23(3),
    491–506 (2012)
42. Tadmor, Y., Tolhurst, D.: Calculating the contrasts that retinal ganglion cells and
    lgn neurones encounter in natural scenes. Vision research 40(22), 3145–3157 (2000)
43. Tang, J., Peli, E., Acton, S.: Image enhancement using a contrast measure in the
    compressed domain. IEEE signal processing LETTERS 10(10), 289–292 (2003)
44. Toet, A.: Image fusion by a ration of low-pass pyramid. Pattern Recognition Letters
    9(4), 245–253 (1989)
45. Triantaphillidou, S., Jarvis, J., Psarrou, A., Gupta, G.: Contrast sensitivity in im-
    ages of natural scenes. Signal Processing: Image Communication 75, 64–75 (2019)
46. Vandergheynst, P., Kutter, M., Winkler, S.: Wavelet-based contrast computation
    and application to digital image watermarking. In: Wavelet Applications in Signal
    and Image Processing VIII. vol. 4119, pp. 82–93. International Society for Optics
    and Photonics (2000)
47. Watson, A.B.: Dct quantization matrices visually optimized for individual images.
    In: Human vision, visual processing, and digital display IV. vol. 1913, pp. 202–216.
    International Society for Optics and Photonics (1993)
48. Watson, A.B. (ed.): Digital Images and Human Vision. MIT Press, Cambridge,
    MA, USA (1993)
49. Watson, A.B., Solomon, J.A.: Model of visual contrast gain control and pattern
    masking. J. Opt. Soc. Am. A 14(9), 2379–2391 (Sep 1997)
        On Measures of Visual Contrast and Their Use in Image Processing         21

50. Winkler, S., Vandergheynst, P.: Computing isotropic local contrast from oriented
    pyramid decompositions. In: Image Processing, 1999. ICIP 99. Proceedings. 1999
    International Conference on. vol. 4, pp. 420–424. IEEE (1999)
51. Zhang, Y., Alam, M.M., Chandler, D.M.: Visually lossless perceptual image coding
    based on natural-scene masking models. Recent Advances in Image and Video
    Coding p. 1 (2016)