Method and apparatus of creating a perceptual harmony map

Abstract

The present invention generally relates to a method and apparatus for creating an harmony map of an image. It also relates to a method and apparatus for assessing the quality of an image. The invention also relates to a method and apparatus for assessing the quality of an image such that an assessment score of the image is obtained from a perceptual map of the image created from the previous method. As an interesting tool for content creator and in the goal of maximizing the artistic effect, the proposed method provides perceptual harmony-guided quality map as well as a score of disharmony of a picture without reference image. The method links two distinct subjects: perceptual quality metric and harmonious hue templates in order to provide new tools to content creator whatever the expertise level.

Description

This application claims the benefit, under 35 U.S.C. §365 of International Application PCT/EP2014/056069, filed Mar. 26, 2014, which was published in accordance with PCT Article 21(2) on Oct. 2, 2014 in English and which claims the benefit of European patent application No. 13305394.2, filed Mar. 28, 2013.

1. FIELD OF INVENTION

The present invention generally relates to a method and apparatus for creating an harmony map of an image. It also relates to a method and apparatus for assessing the quality of an image.

2. TECHNICAL BACKGROUND

When manipulating, editing, improving image, the best quality as well as a certain artistic intent are usually the finality. Nevertheless, although the issues related to objective quality assessment have been widely studied in the context of low level artefacts (such as blur, blockiness, jitter, . . . ), the artistic intent is a problem somehow more tricky and subjective leading to high difficulties in modeling or generalization.

As an intermediate indicator, aesthetic quality metric based on features (colorfulness, line orientation, shape, . . . ) intuitively related to beauty and rules (composition, rules-of-third, skyline, . . . ) are showing up recently in the community.

Depending on the application context, some approaches take advantage of a reference source or do the best effort without any reference when providing an absolute quality measurement.

The problem solves by the invention is to define a method for assessing the color harmony of an image and for assessing what the quality of an image (picture) is.

3. SUMMARY OF THE INVENTION

The invention solves this problem by defining a method for creating a perceptual harmony map of an image which comprises:

• • calculating a multi-resolution decomposition of the image,
  • at each resolution level,
  • for each harmonious color template, computing an optimal rotation angle and an harmony distance for the hue value associated with each pixel of the image at the resolution level, and computing a harmony distance map by weighting the sum of said harmony distances,
  • calculating a contrast map and an entropy activity map,
  • obtaining a perceptual harmony map by integrating together the contrast map the entropy activity map and the harmony distance map, and
  • obtaining the perceptual map by accumulating the perceptual harmony maps.

The invention also relates to a method for assessing the quality of an image characterized in that an assessment score of the image is obtained from a perceptual map of the image created from the previous method.

As an interesting tool for content creator and in the goal of maximizing the artistic effect, the proposed method provides perceptual harmony-guided quality map as well as a score of disharmony of a picture without reference image.

The method links two distinct subjects: perceptual quality metric and harmonious hue templates in order to provide new tools to content creator whatever the expertise level.

The multi-resolution decomposition of the image is processed to calculate a perceptual harmony map at different resolutions of the image in order to mimic the human visual system.

Moreover, each perceptual harmony map is obtained by applying masking functions to the harmony distance map or, in other words, by integrating together the contrast map, the entropy activity map and the harmony distance map. Such masking of the harmony distance map provides information regarding local and global visibility.

The specific nature of the invention as well as other objects, advantages, features and uses of the invention will become evident from the following description of a preferred embodiment taken in conjunction with the accompanying drawings.

4. LIST OF FIGURES

The embodiments will be described with reference to the following figures:

FIG. 1 represents color templates;

FIG. 2 depicts a flowchart of the processing method for creating a perceptual harmony map of an image; and

FIG. 3 depicts a device according comprising means to implement the method for creating a perceptual harmony map of the image.

5. DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT OF THE INVENTION

Harmony is a subjective concept whose difficulties lie on its definition and its measurement. Fedorovskaya et al. (E. Fedorovskaya, C. Neustaedter, and W. Hao. Image harmony for consumer images. In Image Processing, 2008. ICIP 2008. 15^th IEEE International Conference on, pages 121-124, October) investigated through a series of experiments the identification of some image features involving harmony (dis)comfort. They found out that edge contrast, average lightness, range of lightness . . . may influence global harmony appreciation. Despite this previous reference, most of harmony investigations in image processing field relates to color combinations and complementaries. Starting with Itten contrast measurement where harmonious doublets and triplets of color are defined (J. Itten. The art of color: the subjective experience and objective rationale of color. Van Nostrand Reinhold, New York, 1973.), Matsuda (Y. Matsuda. Color design, Asakura Shoten, 1995) extended this work by defining HSV (Hue, Saturation, Value)-based templates that predict set of harmonious hues within an image.

An harmonious color template is thus a set of HSV values (hue, saturation and value) that are considered as rendering/reflecting a global harmonious effect when present at the same time. Height harmonious color templates T_m (mε{i, I, L, T, V, X, Y, J}) are used as depicted on FIG. 1. The harmonious color templates O and N are not used. Each harmonious color template T_m is made of different portions/sectors to handle color complementary, color orthogonality and set of close hues. Each sector has a center α_m and a size w_m. A harmonious color template may be turned around the hue wheel. Thus an harmonious color template is defined by a template type and a sector angle α_m.

A complete implementation of the invention is depicted in FIG. 2. One can notice that the following method can be extended to video source by applying same process to consecutive frames.

At step 10, a multi-resolution decomposition of the image is calculated.

For example, a Discrete Wavelet Transform (DWT) which has proven to be efficient both in term of prediction and computation performances (A. Ninassi, O. LeMeur, P. Le Callet, and D. Barba. On the performance of human visual system based image qual-ity assessment metric using wavelet domain. Pages 680610-680610-12, 2008.). A CDF 9/7 (Cohen-Daubechies-Feauveau) kernel may be used.

Each decomposition level comprises three orientation subbands (horizontal, vertical, and oblique frequencies). The spatial frequency range of a decomposition level lε[0; L−1] is [2^−(l+1)·f_max; 2⁻¹·f_max] where f_max is the maximum spatial frequency of an input image Im.

At step 20, at each resolution level l, for each harmonious color template T_m, an optimal rotation angle α_m,l is computed and an harmony distance HD_m,l(h) is computed for the hue value h associated with each pixel (x,y) of the image Im at the resolution level l. A harmony distance map HDMap_l(x,y) is then calculated for each pixel of the image Im at the resolution level l by weighting the sum of said harmony distances HD_m,l(h).

According to an embodiment, the optimal rotation angle α_m,l which is computed for a harmonious color template T_m minimizes the energy E_αm,l defined by the Kullback-Liebler divergence between the original hue distribution M_h,l of the image Im at the resolution level l and the distribution P_h,l(m, α_m,l) of the harmonious color template T_m for the angle α_m,l such as described, for example by Baveye et al. (Y. Baveye, F. Urban, C. Chamaret, V. Demoulin, and P. Hellier. Saliency-guided consistent color harmonization. In S. Tominaga, R. Schettini, and A. Trmeau, editors, Computational Color Imaging, volume 7786 of Lecture Notes in Computer Science, pages 105-118. Springer Berlin Heidelberg, 2013).

Note that the distributions M_h,l and P_h,l(m, l, α_m,l) may be considered as being histograms with 360 bins.

The distributions M_h,l, i.e. the color histogram M_h,l, is computed in the HSV space as a normalized hue distribution weighted by saturation S(x,y) and value V(x,y):

$M_{h,l}=\frac{1}{{\Sigma }_{\left(x,y\right)}S\left[x,y\right]*V\left[x,y\right]}*\sum _{\left(x,y\right)\in \left\{\left(u,v\right)\backslash H\left[u,v\right]=h\right\}}S\left[x,y\right]*V\left[x,y\right]$

where (x,y) is a pixel belonging to the image Im at each resolution level l.

Mathematically speaking, the energy E_αm,l is given by:

$E_{{\alpha }_{m,l}}=\underset{\alpha }{\min }\sum _h{M}_{h,l}x\ln \left(\frac{M_{h,l}}{P_{h,l}\left(m,\alpha \right)}\right)$

with the distribution P_h,l(m, α) for a harmonious color template T_m with n sectors is given by:

$P_{h,l}\left(m,\alpha \right)=\sum _{k=1}^n{S}_k\left(m,{\alpha }_k\right)$

and a distribution of a sector S_k(m, α_k) given for each bin value i by:

$\hspace{1em}\{\begin{array}{cc}\exp \left(-\frac{1}{1-\left(\frac{2xi-{\alpha }_k}{w_k}\right)}\right)& \mathrm{sih}\in \left[{\alpha }_k-\frac{w_k}{2},{\alpha }_k+\frac{w_k}{2}\right]\\ 0& \mathrm{otherwise}\end{array}$

The invention is not limited by the way the distribution is defined.

According to an embodiment, the harmony distance HD_m,l(h) is an arc-length distance on the hue wheel (measured in degrees) given by:

${\mathrm{HD}}_{m,l}\left(h\right)=\{\begin{array}{cc}0& \mathrm{if}h-{\alpha }_k<\frac{w_k}{2}\\ h-{\alpha }_k-\frac{w_k}{2}& \mathrm{otherwise}\end{array}$

where ∥•∥ is the arc-length distance, α_k is the closest sector center, w_k its size in degrees and h is the hue value of a pixel of the image Im at a resolution level l.

According to an embodiment, the harmony distance map HDMap_l(x,y) is the sum of said harmony distances HD_m,l(h) weighted by the relative energy calculated over the harmonious color templates.

Mathematically speaking, the harmony distance map HDMap_l(x,y) is given by:

${\mathrm{HDMap}}_l\left(x,y\right)=\sum _{T_m}\frac{E_{\mathrm{Total}}-E_{{\alpha }_{m,l}}}{E_{\mathrm{Total}}}{\mathrm{xHD}}_{m,l}\left(h\right)$

At step 30, a contrast map (M_l,o(m,n) is calculated for each resolution level l.

The contrast map models the visibility change of the signal due to contrast values.

According to an embodiment, when a DWT is used at step 10, the contrast masking is defined as the wavelet transformed values at site (m,n), weighted by the CSF function that describes, at a resolution level l the variations in visual sensitivity to the spatial frequency and orientation o.

Mathematically speaking, given a CSF value N_l,o as the mean value of a 2D CSF from Daly (S. Daly. Digital images and human vision. Chapter The visible differences predictor: an algorithm for the assessment of mage) over the spatial frequencies covered by a subband S_l,o, the contrast masking map is given by:CM_l,o(m,n)=s_l,o(m,n)·N_l,o

At step 40, an entropy activity map SML_l(x,y) is calculated for each resolution level l.

The entropy activity map reflects the uncertainty of the image such as in texture area.

According to an embodiment, the entropy activity map SML_l is evaluated by the computation of entropy on a n-by-n neighborhood.

Preferably, the entropy activity map SML_l is evaluated by computing as the minimum of the directional gradients (horizontal, vertical, diagonal) for each pixel (x,y) of the image Im. Such a semi-local activity map SML_l is then given by

${\mathrm{SML}}_l\left(x,y\right)=\min \left(\frac{\partial I}{\partial x},\frac{\partial I}{\partial y},\frac{1}{2}·\left(\frac{\partial I}{\partial x}+\frac{\partial I}{\partial y}\right)\right)$

where I is the input image.

The harmony distance maps are computed at differents resolutions level to match the HVS model.

At step 50, for each resolution level l, the contrast map CM_l,o (m,n), the entropy activity map SML_l(x,y) and the harmony distance map HDMap_l(x,y) are integrated together to get a perceptual harmony map PHDMap_l(x,y).

According to an example, said perceptual harmony distance map PHDMap_l(x,y) is given by:

${\mathrm{PHDMap}}_l\left(m,n\right)=\frac{{\mathrm{HDMap}}_l\left(m,n\right)}{1+\frac{1}{2}\left({\Sigma }_o{\mathrm{CM}}_{l,o}\left(m,n\right)+{\mathrm{SML}}_i\left(m,n\right)\right)}$

where (m,n) are the coordinates corresponding to the pixel (x,y) of the image Im at resolution level l.

At step 60, the perceptual harmony maps PHDMap_l(m,n) are accumulated to built a perceptual map PHMap (x,y).

Mathematically speaking, the perceptual map PHMap (x,y) is given by:

$\mathrm{PHMap}\left(x,y\right)=\sum _l{\mathrm{PHDMap}}_l\left(m,n\right)$

Note the perceptual map PHMap (x,y) is visually close to a harmony distance map computed by accumulating the harmony distance map HDMap_l(x,y), but integrates masking effects by decreasing impact of small group of pixels.

According to an aspect of the invention, a method for assessing the quality of an image is defined by obtaining a rating or assessment score R of the image from the perceptual map PHMap (x,y) (step 70 in FIG. 2).

According to an embodiment, the rating R is computed from by means of a Minkowski summation given by:

$R={\left(\frac{1}{W·H}\sum _y^H\sum _x^W{\left(\mathrm{PHMap}\left(x,y\right)\right)}^{\beta }\right)}^{\frac{1}{\beta }}$

where W and H are respectively the width and height of the image Im and β is a parameter set for example to 2.

The color information of an image at a resolution level is carried by the lowest subband.

Consequently, according to a variant of the method, the contrast and entropy activity map are integrated with the harmony distance map at a given resolution level only for the lowest subband of this resolution level.

FIG. 3 represents an exemplary architecture of a processing device 300 according to a specific and non limiting embodiment. The processing device can be for example a tablet, a PDA or a cell phone. Processing device 300 comprises following elements that are linked together by a data and address bus 30:

• • a microprocessor 31 (or CPU), which is, for example, a DSP (or Digital Signal Processor);
  • a ROM (or Read Only Memory) 32;
  • a RAM (or Random Access Memory) 33;
  • one or several Input/Output interface(s) 35, for example a keyboard, a mouse; and
  • a battery 36.

Each of these elements of FIG. 3 are well known by those skilled in the art and won't be disclosed further. The processing device 300 may comprise display means such as a screen for displaying the processed images. Algorithm of the processing method according to the invention are stored in the ROM 32. RAM 33 comprises in a register, the program executed by the CPU 31 and uploaded after switch on of the processing device 300. When switched on, the CPU 21 uploads the program in the RAM and executes the corresponding instructions. The images to be processed are received on one of the Input/Output interfaces 35. One of the Input/Output interface 35 is adapted to transmit the images processed according to the invention.

According to variants, processing devices 300 compatible with the invention are implemented according to a purely hardware realisation, for example in the form of a dedicated component (for example in an ASIC (Application Specific Integrated Circuit) or FPGA (Field-Programmable Gate Array) or VLSI (Very Large Scale Integration) or of several electronic components integrated into a device or even in a form of a mix of hardware elements and software elements.

While not explicitly described, the present embodiments and variants may be employed in any combination or sub-combination.

Claims

1. A method performed by a device for creating a perceptual harmony map of an image from at least one harmonious color template defined by a center and at least a sector wherein the method comprises: calculating, by a processor, a multi-resolution decomposition of the image,at each resolution level, for each harmonious color template, computing, by said processor, an optimal rotation angle of said harmonious color template around its center and an harmony distance for the hue value associated with each pixel of the image at the resolution level, said harmony distance quantifying the distance between said hue value and a sector of said harmonious color template, and computing a harmony distance map by weighting the sum of harmony distances calculated for the harmonious color templatescalculating, by said processor, a contrast map and an entropy activity map reflecting the uncertainty of the image,obtaining, by said processor, a perceptual harmony map by integrating together the contrast map the entropy activity map and the harmony distance map, andobtaining, by said processor, the perceptual map by accumulating the perceptual harmony maps calculated over the different resolution levels.

2. The method of claim 1, wherein an assessment score of the image is obtained, by said processor, from the perceptual map.

3. A processor readable non-transitory medium having stored therein instructions for causing a processor to perform at least the steps of the method according to claim 1.

4. A processor readable non-transitory medium having stored therein instructions for causing a processor to perform at least the steps of the method according to claim 2.

5. Non-transitory storage medium carrying instructions of program code for executing steps of the method according to claim 1, when said program is executed on a computing device.

6. Device for creating a perceptual harmony map of an image from at least one harmonious color template defined by a center and at least a sector wherein it comprises a processor configured to: calculate a multi-resolution decomposition of the image,at each resolution level, for each harmonious color template, compute an optimal rotation angle of said harmonious color template around its center and an harmony distance for the hue value associated with each pixel of the image at the resolution level, said harmony distance quantifying the distance between said hue value and a sector of said harmonious color template, and compute a harmony distance map by weighting the sum of harmony distances calculated for the harmonious color templatescalculate a contrast map and an entropy activity map reflecting the uncertainty of the image,obtain a perceptual harmony map by integrating together the contrast map the entropy activity map and the harmony distance map, andobtain the perceptual map by accumulating the perceptual harmony maps calculated over the different resolution levels.

7. A device for assessing the quality of an image, wherein it comprises a processor configured to obtain an assessment score of the image from a perceptual map of the image, the device wherein the processor is further configured to: calculate a multi-resolution decomposition of the image,at each resolution level, for each harmonious color template, compute an optimal rotation angle of said harmonious color template around its center and an harmony distance for the hue value associated with each pixel of the image at the resolution level, said harmony distance quantifying the distance between said hue value and a sector of said harmonious color template, and compute a harmony distance map by weighting the sum of harmony distances calculated for the harmonious color templatescalculate a contrast map and an entropy activity map reflecting the uncertainty of the image,obtain a perceptual harmony map by integrating together the contrast map the entropy activity map and the harmony distance map, andobtain the perceptual map by accumulating the perceptual harmony maps calculated over the different resolution levels.