# Patent

## Title
Device and method for creating a saliency map of an image

## Abstract

The invention concerns a device and a method for creating a saliency map of an image. It comprises the steps of:

    
    
        Projection of said image according to the luminance component and if said image is a color image, according to the luminance component and according to the chrominance components,
        Perceptual sub-bands decomposition of said components according to the visibility threshold of a human eye,
        Extraction of the salient elements of the sub-bands related to the luminance component,
        Contour enhancement of said salient elements in each sub-band related to the luminance component,
        Calculation of a saliency map from the contour enhancement, for each sub-band related to the luminance component.
        Creation of the saliency map as a function of the saliency maps obtained for each sub-band.

## Background
This application claims the benefit, under 35 U.S.C. §365 of International Application PCT/EP04/053471, filed Dec. 14, 2004, which was published in accordance with PCT Article 21(2) on Jun. 30, 2005 in English and which claims the benefit of European patent application No. 03293216.2, filed Dec. 18, 2003.
The invention is related to a device and a method for creating a saliency map of an image.
The human information processing system is intrinsically a limited system and especially for the visual system. In spite of the limits of our cognitive resources, this system has to face up to a huge amount of information contained in our visual environment. Nevertheless and paradoxically, humans seem to succeed in solving this problem since we are able to understand our visual environment.
It is commonly assumed that certain visual features are so elementary to the visual system that they require no attentional resources to be perceived. These visual features are called pre-attentive features.
According to this tenet of vision research, human attentive behavior is shared between pre-attentive and attentive processing. As explained before, pre-attentive processing, so-called bottom-up processing, is linked to involuntary attention. Our attention is effortless drawn to salient parts of our view. When considering attentive processing, so-called top-down processing, it is proved that our attention is linked to a particular task that we have in mind. This second form of attention is thus a more deliberate and powerful one in the way that this form of attention requires effort to direct our gaze towards a particular direction.
The detection of the salient points in an image enable the improvement of further steps such as coding or image indexing, watermarking, video quality estimation.
The known approaches are more or less based on non-psycho visual features. In opposition with such methods, the proposed method relies on the fact that the model is fully based on the human visual system (HVS) such as the computation of early visual features.
In a first aspect, the invention proposes a method for creating a saliency map of an image comprising the steps of:

    
    
        projection of said image according to the luminance component and if said image is a color image, according to the luminance component and according to the chrominance components,
        perceptual sub-bands decomposition of said components according to the visibility threshold of a human eye,
        extraction of the salient elements of the sub-bands related to the luminance component,
        contour enhancement of said salient elements in each sub-band related to the luminance component,
        calculation of a saliency map from the contour enhancement, for each sub-band related to the luminance component.
        creation of the saliency map as a function of the saliency maps obtained for each sub-band.
    
    


In a second aspect, the invention proposes a device for creating a saliency map of an image characterized in that it comprises means for:

    
    
        Projecting said image according to the luminance component sand if said image is a color image, according to the luminance component and according to the chrominance components,
        Transposing into the frequential domains said luminance and chrominance signals,
        Decomposing into perceptual sub-bands said components of the frequential domain according to the visibility threshold of a human eye,
        Extracting the salient elements of the sub-bands related to the luminance component,
        Contour enhancing said salient elements in each sub-band related to the luminance component,
        Calculating a saliency map from the contour enhancement, for each sub-band related to the luminance component.
        Creating the saliency map as a function of the saliency maps obtained for each sub-band.
    
    





Other characteristics and advantages of the invention will appear through the description of a non-limiting embodiment of the invention, which will be illustrated, with the help of the enclosed drawings, wherein:
FIG. 1 represents a general flow-chart of a preferred embodiment of the method according to the invention applied to a black and white image,
FIG. 2 represents a general flow-chart of a preferred embodiment of the method according to the invention applied to a black and white image,
FIG. 3 represents the psycho visual spatial frequency partitioning for the achromatic component,
FIG. 4 represents the psycho visual spatial frequency partitioning for the chromatic components,
FIG. 5 represents the Daily Contrast Sensitivity Function,
FIGS. 6a and 6b represent respectively the visual masking and a non linear model of masking,
FIG. 7 represents the flow-chart of the normalisation step according to the preferred embodiment,
FIG. 8 represents the inhibition/excitation step,
FIG. 9 represents the profile of the filters to model facilitative interactions for θ=0,
FIG. 10 represents an illustration of the operator D(z),
FIG. 11 represents the chromatic reinforcement step,
FIG. 12 represents the non CRF exhibition caused by the adjacent areas of the CRF flanks,
FIG. 13 represents a profile example of the normalized weighting function for a particular orientation and radial frequency.



FIG. 1 represents the general flow-chart of a preferred embodiment of the method according to the invention applied to a black and white image.
The algorithm is divided in three main parts.
The first one named visibility is based on the fact that the human visual system (HVS) has a limited sensitivity. For example, the HVS is not able to perceive with a good precision all signals in your real environment and is insensible to small stimuli. The goal of this first step has to take into account these intrinsic limitations by using perceptual decomposition, contrast sensitivity functions (CSF) and masking functions.
The second part is dedicated to the perception concept. The perception is a process that produces from images of the external world a description that is useful to the viewer and not cluttered with irrelevant information. To select relevant information, a center surround mechanism is notably used in accordance with biological evidences.
The last step concerns some aspects of the perceptual grouping domain. The perceptual grouping refers to the human visual ability to extract significant images relations from lower level primitive image features without any knowledge of the image content and group them to obtain meaningful higher-level structure. The proposed method just focuses on contour integration and edge linking.
Steps E3, E4 are executed on the signal in the frequential domain.
Steps E1, E6 and E9 are done in the spatial domain.
Steps E7 and E8 are done in the frequential or spatial domain. If they are done in the frequential domain, a Fourier transformation has to be carried on before step E7 and an inverse Fourier transformation has to be carried out before step E9.
In step E1, the luminance component is extracted from the considered image.
In step E2, the luminance component is transposed into the frequency domain by using known transformations such as the Fourier transformation in order to be able to apply in step E3, the perceptual sub-band decomposition on the image.
In step E3, a perceptual decomposition is applied on the luminance component. This decomposition is inspired from the cortex transform and based on the decomposition proposed in the document “The computation of visual bandwidths and their impact in image decomposition and coding”, International Conference and Signal Processing Applications and Technology, Santa-Clara, Calif., pp. 776-770, 1993. This decomposition is done according to the visibility threshold of a human eye.
The decomposition, based on different psychophysics experiments, is obtained by carving up the frequency domain both in spatial radial frequency and orientation. The perceptual decomposition of the component A leads to 17 psycho visual sub-bands distributed on 4 crowns as shown on FIG. 3.
The shaded region on the FIG. 3 indicates the spectral support of the sub-band belonging to the third crown and having an angular selectivity of 30 degrees, from 15 to 45 degrees.
Four domains (crowns) of spatial frequency are labeled from I to IV:
I: spatial frequencies from 0 to 1.5 cycles per degree;
II: spatial frequencies from 1.5 to 5.7 cycles per degree;
III: spatial frequencies from 5.7 to 14.2 cycles per degree;
IV: spatial frequencies from 14.2 to 28.2 cycles per degree.
The angular selectivity depends on the considered frequency domain. For low frequencies, there is no angular selectivity.
The main properties of these decompositions and the main differences from the cortex transform are a non-dyadic radial selectivity and an orientation selectivity that increases with the radial frequency.
Each resulting sub-band may be regarded as the neural image corresponding to a population of visual cells tuned to a range of spatial frequency and a particular orientation. In fact, those cells belong to the primary visual cortex (also called striate cortex or V1 for visual area 1). It consists of about 200 million neurons in total and receives its input from the lateral geniculate nucleus. About 80 percent of the cells are selective for orientation and spatial frequency of the visual stimulus.
On the image spatial spectrum, a well-known property of the HVS is applied, which is known as the contrast sensitivity function (CSF). The CSF applied is a multivariate function mainly depending on the spatial frequency, the orientation and the viewing distance.
Biological evidences have shown that visual cells response to stimuli above a certain contrast. The contrast value for which a visual cell response is called the visibility threshold (above this threshold, the stimuli is visible). This threshold varies with numerous parameters such as the spatial frequency of the stimuli, the orientation of the stimuli, the viewing distance, . . . This variability leads us to the concept of the CSF which expresses the sensitivity of the human eyes (the sensitivity is equal to the inverse of the contrast threshold) as a multivariate function. Consequently, the CSF permits to assess the sensitivity of the human eyes for a given stimuli.
In step E4, a 2 D anisotropic CSF designed by Dally is applied. Such a CSF is described in document “the visible different predictor: an algorithm for the assessment of image fidelity”, in proceedings of SPIE Human vision, visual processing and digital display III, volume 1666, pages 2-15, 1992.
The CSF enables the modelisation of an important property of the eyes, as the SVH cells are very sensitive to the spatial frequencies.
On FIG. 5, the Dally CSF is illustrated.
Once the Dally function has been applied, an inverse Fourrier Transformation is applied on the signal in step E5, in order to be able to apply the next step E6.
For natural pictures, the sensitivity can be modulated (increased or decreased the visibility threshold) by the presence of another stimulus. This modulation of the sensitivity of the human eyes is called the visual masking, as done in step E6.
An illustration of masking effect is shown on the FIGS. 6a and 6b. Two cues are considered, a target and a masker where CT and CM are the contrast threshold of the target in the presence of the masker and the contrast of the masker respectively. Moreover, CT0 is the contrast threshold measured by a CSF (without masking effect).
When CM varies, three regions can be defined:

    
    
        At low values of CM, the detection threshold remains constant. The visibility of the target is not modified by the masker.
        When CM tends toward CT0, the masker eases the detection of the target by decreasing the visibility threshold. This phenomenon is called facilitative or pedestal effect.
        When CM increases, the target is masked by the masker. His contrast threshold increases.
    
    


The visual masking method is based on the detection of a simple signal as sinusoidal patterns.
There are several other methods to achieve the visual masking modeling based on psychophysics experiments: for instance, a best method refers to the detection of quantization noise.
It is obvious that the preferred method is a strong simplification as regard the intrinsic complexity of natural pictures. Nevertheless, numerous applications (watermarking, video quality assessment) are built around such principle with interesting results compared to the complexity.
In the context of sub-band decomposition, masking has been intensively studied leading to define three kinds of masking: intra-channel masking, inter-channel masking and inter-component masking.
The intra-channel masking occurs between signals having the same features (frequency and orientation) and consequently belonging to the same channel. It is the most important masking effect.
The inter-channel masking occurs between signals belonging to different channels of the same component.
The inter-component masking occurs between channels of different components (the component A and one chromatic component for example). These two last visual masking are put together and are just called inter-masking in the following.
For the achromatic component, we used the masking function designed by Dally in document entitled “A visual model for Optimizing the Design of Image Processing Algorithms”, in IEEE international conferences on image processing, pages 16-20, 1994, in spite of the fact that this model doesn't take into account the pedestal effect. The strength of this model lies in the fact that it has been optimized with a huge amount of experimental results.
The variation of the visibility threshold is given by:

Ti,j,Aintra(m,n)=(1+(k1(k2|Ri,j(m,n)|)s)b1/b 

where Ri,j is a psycho visual channel stemming from the perceptual channel decomposition (For example, the shaded region on the FIG. 2.1 leads to the channel RIII,2). The values k1, k2, s, b are given below:

k1=0.0153
k2=392.5
The below table gives the values of s and b according to the considered sub-band:












Sub-band
s
b

















I
0.75
4


II
1
4


III
0.85
4


IV
0.85
4









We get the signal R′i,j(x,j) at the output of the masking step.

R′i,j(x,y)=Ri,j(x,y)/Ti,j(x,y)

Then, in step E7, the step of normalization enables to extract the main important information from the sub-band. Step E7 is detailed on FIG. 7.
In step S1, a first sub-band R′i,j(x,y) is selected. The steps S2 to S4 and S8 are carried on for each sub-band R′i,j(x,y) of the 17 sub-bands.
The steps S5 to S7 are done for the second crown (II).
I represents the spatial radial frequency band, I belongs to {I, II, III, IV}.
J represents the orientation, j belongs to {1, 2, 3, 4, 5, 6},
(x,y) represent the spatial coordinates.
In other embodiments, the different steps can be carried out on all the sub-bands.
Steps S2 and S3 aim to modelize the behavior of the classical receptive field (CRF).
The concept of CRF permits to establish a link between a retinal image and the global percept of the scene. The CRF is defined as a particular region of visual field within which an appropriate stimulation (with preferred orientation and frequency) provokes a relevant response stemming from visual cell. Consequently, by definition, a stimulus in the outer region (called surround) cannot activate the cell directly.
The inhibition and excitation in steps S2 and S3 are obtained by a Gabor filter, which is sensible as for the orientation and the frequency.
The Gabor filter can be represented as following:

gabor(x,y,σx,σy,f,θ)=Gσx,σy(xθ,yθ)cos(2πfxθ)

f being the spatial frequency of the cosinus modulation in cycles per degree (cy/°).
(xθ, yθ) are obtained by a translation of the original coordinates (x0, y0) and by a rotation of θ,



  
    
      [
      
        
          
            
              x
              θ
            
          
        
        
          
            
              y
              θ
            
          
        
      
      ]
    
    =
    
      
        [
        
          
            
              
                cos
                ⁢
                
                  
                
                ⁢
                θ
              
            
            
              
                sin
                ⁢
                
                  
                
                ⁢
                θ
              
            
          
          
            
              
                
                  -
                  sin
                
                ⁢
                
                  
                
                ⁢
                θ
              
            
            
              
                cos
                ⁢
                
                  
                
                ⁢
                θ
              
            
          
        
        ]
      
      ⁡
      
        [
        
          
            
              
                x
                -
                
                  x
                  0
                
              
            
          
          
            
              
                y
                -
                
                  y
                  0
                
              
            
          
        
        ]
      
    
  




  
    
      
        G
        
          
            σ
            x
          
          ,
          
            σ
            y
          
        
      
      ⁡
      
        (
        
          x
          ,
          y
        
        )
      
    
    =
    
      A
      ⁢
      
        
      
      ⁢
      exp
      ⁢
      
        {
        
          
            -
            
              
                (
                
                  x
                  
                    
                      2
                    
                    ⁢
                    
                      σ
                      x
                    
                  
                
                )
              
              2
            
          
          -
          
            
              (
              
                y
                
                  
                    2
                  
                  ⁢
                  
                    σ
                    y
                  
                
              
              )
            
            2
          
        
        }
      
    
  



    
    
        A representing the amplitude,
    
    


σx et σy representing the width of the gaussian envelop along the x and y axis respectively.




  
    excitation
    ⁡
    
      (
      
        x
        ,
        y
        ,
        
          σ
          x
        
        ,
        
          σ
          y
        
        ,
        f
        ,
        θ
      
      )
    
  
  =
  
    {
    
      
        
          
            gabor
            ⁡
            
              (
              
                x
                ,
                y
                ,
                
                  σ
                  x
                
                ,
                
                  σ
                  y
                
                ,
                f
                ,
                θ
              
              )
            
          
        
        
          
            
              if
              -
              
                1
                /
                
                  (
                  
                    4
                    ⁢
                    f
                  
                  )
                
              
            
            ≤
            
              x
              θ
            
            ≤
            
              1
              /
              
                (
                
                  4
                  ⁢
                  f
                
                )
              
            
          
        
      
      
        
          0
        
        
          otherwise
        
      
    
  




In order to obtain elliptic shapes, we take different variances σx<σy.
Finally, we get the output:

Ri,jExc(x,y)=R′i,j(x,y)*excitation(x,y,σx,σy,f,θ)

In step S3, the inhibition is calculated by the following formula:




  
    inhibition
    ⁡
    
      (
      
        x
        ,
        y
        ,
        
          σ
          x
        
        ,
        
          σ
          y
        
        ,
        f
        ,
        θ
      
      )
    
  
  =
  
    {
    
      
        
          0
        
        
          
            
              si
              -
              
                1
                /
                
                  (
                  
                    4
                    ⁢
                    f
                  
                  )
                
              
            
            ≤
            
              x
              θ
            
            ≤
            
              1
              /
              
                (
                
                  4
                  ⁢
                  f
                
                )
              
            
          
        
      
      
        
          
            
            
              gabor
              ⁡
              
                (
                
                  x
                  ,
                  y
                  ,
                  
                    σ
                    x
                  
                  ,
                  
                    σ
                    y
                  
                  ,
                  f
                  ,
                  θ
                
                )
              
            
            
          
        
        
          
            sinon
            .
          
        
      
    
  




And finally:

Ri,jInh(x,y)=R′i,j(x,y)*inhibition(x,y,σx,σy,fθ)

In step S4, the difference between the excitation and the inhibition is done. The positive components are kept, the negative components are set to “0”. This is the following operation,

R′′i,j(x,y)=|Ri,jExc(x,y)−Ri,jInh(x,y)|>0 

In step S5, for each orientation, for each sub-band of the second domain, two convolution products are calculated:

Li,j0(x,y)=R′′i,j(x,y)*Bi,j0(x,y)

Li,j1(x,y)=R′′i,j(x,y)*Bi,j1(x,y)

Bi,j0(x,y) and Bi,j1(x,y) are 2 half-butterfly filters. The profile of these filters allow the modelling of facilitative interactions for θ=0 given on FIG. 9. These filters are defined by using a bipole/butterfly filter.
It consists of a directional term Dθ(x,y) and a proximity term generated by a circle Cr blurred by a gaussian filter Gσx,σy(x,y).



  
    
      
        B
        
          
            θ
            
              i
              ,
              j
            
          
          ,
          α
          ,
          r
          ,
          σ
        
      
      ⁡
      
        (
        
          x
          ,
          y
        
        )
      
    
    =
    
      
        
          
            D
            
              θ
              
                i
                ,
                j
              
            
          
          ⁡
          
            (
            
              x
              ,
              y
            
            )
          
        
        ·
        
          C
          r
        
      
      *
      
        
          G
          
            
              σ
              x
            
            ,
            
              σ
              y
            
          
        
        ⁡
        
          (
          
            x
            ,
            y
          
          )
        
      
    
  




  with




  
    
      
        D
        
          θ
          
            i
            ,
            j
          
        
      
      ⁡
      
        (
        
          x
          ,
          y
        
        )
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    
                      cos
                      ⁡
                      
                        (
                        
                          
                            
                              π
                              /
                              2
                            
                            α
                          
                          ⁢
                          φ
                        
                        )
                      
                    
                    ⁢
                    si
                  
                
                
                  
                    φ
                    <
                    α
                  
                
              
              
                
                  0
                
                
                  
                    sinon
                    .
                  
                
              
            
            ⁢
            
              

            
            ⁢
            and
            ⁢
            
              

            
            ⁢
            φ
          
          =
          
            arctan
            ⁡
            
              (
              
                
                  y
                  ′
                
                /
                
                  x
                  ′
                
              
              )
            
          
        
        ,
      
    
  



where (x′,y′)T is the vector (x,y)T rotated by θi,j. The parameter α defines the opening angle 2α of the bipole filter. It depends on the angular selectivity γ of the considered sub-band. We take α=0.4×γ. The size of the bipole filter is about twice the size of the CRF of a visual cell.

In step S6, we compute the facilitative coefficient:




  
    
      f
      
        i
        ,
        j
      
      iso
    
    ⁡
    
      (
      
        x
        ,
        y
      
      )
    
  
  =
  
    D
    ⁡
    
      (
      
        
          
            
              L
              
                i
                ,
                j
              
              1
            
            ⁡
            
              (
              
                x
                ,
                y
              
              )
            
          
          +
          
            
              L
              
                i
                ,
                j
              
              0
            
            ⁡
            
              (
              
                x
                ,
                y
              
              )
            
          
        
        
          max
          ⁡
          
            (
            
              β
              ,
              
                
                
                  
                    
                      L
                      
                        i
                        ,
                        j
                      
                      1
                    
                    ⁡
                    
                      (
                      
                        x
                        ,
                        y
                      
                      )
                    
                  
                  -
                  
                    
                      L
                      
                        i
                        ,
                        j
                      
                      0
                    
                    ⁡
                    
                      (
                      
                        x
                        ,
                        y
                      
                      )
                    
                  
                
                
              
            
            )
          
        
      
      )
    
  




with,

    
    
        β a constant,
    
    






  
    D
    ⁡
    
      (
      z
      )
    
  
  =
  
    {
    
      
        
          
            
              
                
                  0
                  ,
                  
                    z
                    ≤
                    
                      s
                      1
                    
                  
                  ,
                
              
            
            
              
                
                  
                    α
                    1
                  
                  ,
                  
                    z
                    ≤
                    
                      s
                      2
                    
                  
                  ,
                
              
            
            
              
                ...
              
            
            
              
                
                  
                    α
                    
                      N
                      -
                      1
                    
                  
                  ,
                  
                    z
                    ≤
                    
                      s
                      
                        N
                        -
                        1
                      
                    
                  
                
              
            
          
          ⁢
          where
          ⁢
          
            
          
          ⁢
          
            α
            i
          
        
        ≤
        1
      
      ,
      
        i
        ∈
        
          [
          
            
              0
              ⁢
              
                
              
              ⁢
              ...
              ⁢
              
                
              
              ⁢
              N
            
            -
            1
          
          ]
        
      
    
  




An illustration of the operator D(z) is given on FIG. 9.
To ease the application of the facilitative coefficient, the operator D(z) ensures that the facilitative coefficient is constant by piece as shown on FIG. 9.
In step S7, the facilitative coefficient is applied to the normalized result obtained in step S4.

R′′′i,j(x,y)=R′′i,j(x,y)×(1+fi,jiso(x,y))

Going back to step E8 of FIG. 1, after step S7 of FIG. 7, the four saliency maps obtained for the domain II are combined to get the whole saliency map according to the following formula:

fixation(x,y)=α×R′′′II,0(x,y)+β×R′′′II,1(x,y)+χ×R′′′II,2(x,y)+δ×R′′′II,3(x,y)

α, β, χ, δ represent weighting coefficients which depend on the application (watermarking, coding . . . ).
In other embodiments, the saliency map can be obtained by a calculation using the whole 17 sub-bands and not only the sub-bands of domain II.
FIG. 2 represents the general flow-chart of a preferred embodiment of the method according to the invention applied to a colour image.
Steps T1, T4, T′4, T′′4, T5 and T8 are done in the spatial domain.
Steps T2, T′2, T′′2, T3, T′3, T′′3 are done in the frequencial domain.
A Fourier transformation is applied on three components between step T1 and steps T2, T′2, T′′2.
An inverse Fourier transformation is applied between respectively T3, T′3, T′′3 and T4, T′4 and T′′4.
Steps T6 and T7 can be done in the frequential or spatial domain. If they are done in the frequential domain, a Fourier transformation is done on the signal between steps T5 and T6 and an inverse Fourier transformation is done between steps T7 and T8.
Step T1 consists in converting the RGB luminances into the Krauskopf's opponent-colors space composed by the cardinal directions A, Cr1 and Cr2.
This transformation to the opponent-colors space is a way to decorrelate color information. In fact, it's believed that the brain uses 3 different pathways to encode information: the first conveys the luminance signal (A), the second the red and green components (Cr1) and the third the blue and yellow components (Cr2).
These cardinal directions are in closely correspondence with signals stemming from the three types of cones (L,M,S).
Each of the three components RGB firstly undergoes a power-law non linearity (called gamma law) of the form xγ with γ ̃2.4. This step is necessary in order to take into account the transfer function of the display system. The CIE (French acronym for “commission internationale de I'éclairage”) XYZ tristimulus value which form the basis for the conversion to an HVS color space is then computed by the following equation:




  
    (
    
      
        
          X
        
      
      
        
          Y
        
      
      
        
          Z
        
      
    
    )
  
  =
  
    
      (
      
        
          
            0.412
          
          
            0.358
          
          
            0.18
          
        
        
          
            0.213
          
          
            0.715
          
          
            0.072
          
        
        
          
            0.019
          
          
            0.119
          
          
            0.95
          
        
      
      )
    
    ⁢
    
      (
      
        
          
            R
          
        
        
          
            G
          
        
        
          
            B
          
        
      
      )
    
  




The response of the (L,M,S) cones are computed as follows:




  
    (
    
      
        
          L
        
      
      
        
          M
        
      
      
        
          S
        
      
    
    )
  
  =
  
    
      (
      
        
          
            0.240
          
          
            0.854
          
          
            
              -
              0.0448
            
          
        
        
          
            
              -
              0.389
            
          
          
            1.160
          
          
            0.085
          
        
        
          
            
              -
              0.001
            
          
          
            0.002
          
          
            0.573
          
        
      
      )
    
    ⁢
    
      (
      
        
          
            X
          
        
        
          
            Y
          
        
        
          
            Z
          
        
      
      )
    
  




From the LMS space, one has to obtain an opponent color space. There is a variety of opponent color space, which differ in the way to combine the different cones responses. From experimental experiments, the color space designed by Krauskopf have been validated and it is given by the following transformation:




  
    (
    
      
        
          A
        
      
      
        
          
            Cr
            ⁢
            
              
            
            ⁢
            1
          
        
      
      
        
          
            Cr
            ⁢
            
              
            
            ⁢
            2
          
        
      
    
    )
  
  =
  
    
      (
      
        
          
            1
          
          
            1
          
          
            0
          
        
        
          
            1
          
          
            
              -
              1
            
          
          
            0
          
        
        
          
            
              -
              0.5
            
          
          
            
              -
              0.5
            
          
          
            1
          
        
      
      )
    
    ⁢
    
      (
      
        
          
            L
          
        
        
          
            M
          
        
        
          
            S
          
        
      
      )
    
  




Then, in step T2, a perceptual decomposition is applied to the luminance component. Preliminary to step T2 and further to step T1, the luminance component is transposed into the frequency domain by using known transformations such as the Fourier transformation in order to be able to apply in step T2, the perceptual sub-band decomposition on the image.
The perceptual sub-band decomposition of step T2 is the same as the step E3 of FIG. 1, and thus will not be described here, as described earlier.
Concerning the decomposition of chromatic components Cr2 and Cr1, of steps T′2 and T′′2 as shown on FIG. 4, the decomposition leads to 5 psychovisual sub-bands for each of these components distributed on 2 crowns. Preliminary to steps T′2, T′′2 and further to step T1, the chrominance components are transposed into the frequency domain by using known transformations such as the Fourier transformation in order to be able to apply in step T′2 and T′′2, the perceptual sub-band decomposition on the image.
Two domains of spatial frequency are labeled from I to II:
I: spatial frequencies from 0 to 1.5 cycles per degree,
II: spatial frequencies from 1.5 to 5.7 cycles per degree.
In steps T3, T′3 and T′′3, a contract sensitivity function (CSF) is applied.
In step T3, the same contrast sensitivity as in step E4 of FIG. 1 is performed on the luminance component and thus will not be described here.
In step T′3 and T′′3, the same CSF is applied on the two chromatic components Cr1 and Cr2. On the two chromatic components, a two-dimensional anisotropic CSF designed by Le Callet is applied. It is described in document <<critères objectifs avec références de qualité visuelle des images couleurs>> of Mr Le Callet, university of Nantes, 2001.
This CSF uses two low-pass filters with a cut-off frequency of about 5.5 cyclse per degree and 4.1 cycles per degree respectively for Cr1 and Cr2 components.
In order to permit the direct comparison between early visual features stemming from different visual modalities (achromatic and chromatic components), the sub-bands related to the visibility are weighted. The visibility threshold being defined as the stimulus's contrast at a particular point for which the stimulus just becomes visible.
An inverse Fourier transformation is then applied on the different components (not shown on FIG. 2) in order to be able to apply the masking in the spatial domain.
Then, an intra masking is applied on the different sub-bands for the chromatic components Cr1 and Cr2 during step T′4 and T′′4 and for the achromatic component in step T4. This last step has already been explained in the description of FIG. 1, step E6. Thus, it will not be described again here.
Intra channel masking is incorporated as a weighing of the outputs of the CSF function. Masking is a very important phenomenon in perception as it describes interactions between stimuli. In fact, the visibility threshold of a stimulus can be affected by the presence of another stimulus.
Masking is strongest between stimuli located in the same perceptual channel or in the same sub-band. We apply the intra masking function designed by Dally on the achromatic component as described on FIG. 1, step E6 and, on the color component, the intra masking function described in document of P. Le Callet and D. Barba, “Frequency and spatial pooling of visual differences for still image quality assessment”, in Proc. SPIE Human Vision and Electronic Imaging Conference, San Jose, Calif., Vol. 3959, January 2000.
These masking functions consist of non linear transducer as expressed in document of Legge and Foley, “Contrast Masking in Human Vision”, Journal of the Optical Society of America, Vol. 70, N° 12, pp. 1458-1471, December 1980.
Visual masking is strongest between stimuli located in the same perceptual channel (intra-channel masking). Nevertheless, as shown in numerous studies, there are several interactions called inter-component masking providing a masking or a pedestal effect. From psychophysics experiments, significant inter-components interactions involving the chromatic components have been elected. Consequently, the sensitivity of achromatic component could be increased or decreased by the Cr1 component. The influence of the Cr2 on the achromatic component is considered insignificant. Finally, the Cr1 can also modulate the sensitivity of Cr2 component (and vice versa).
Then in step T5, a chromatic reinforcement is done.
The colour is one of the strongest attractor of the attention and the invention wants to take advantage of this attraction strength by putting forward the following property: the existence of regions showing a sharp colour and fully surrounded of areas having quite other colours implies a particular attraction to the borders of this region.
To avoid the difficult issue of aggregating measures stemming from achromatic and chromatic components, the colour facilitation consists in enhancing the saliency of achromatic structure by using a facilitative coefficient computed on the low frequencies of the chromatic components.
In the preferred embodiment, only a sub-set of the set of achromatic channels is reinforced. This subset contains 4 channels having an angular selectivity equal to π/4 and a spatial radial frequency (expressed in cyc/deg) belonging to [1.5, 5.7]. One notes these channels Ri,j where i represents the spatial radial frequency and j pertains to the orientation. In our example, j is equal to {0, π/4, π/2, 3π/4}. In order to compute a facilitative coefficient, one determines for each pixel of the low frequency of Cr1 and Cr2 the contrast value related to the content of adjacent areas and to the current orientation of the reinforced achromatic channel as illustrated on the FIG. 11. On FIG. 11, the contrast value is obtained by computing the absolute difference between the average value of the set A and the average value of the set B. The sets A and B belong to the low frequency of Cr1 or Cr2 and are oriented in the preferred orientation of the considered achromatic channel.
The chromatic reinforcement is achieved by the following equation, on an achromatic (luminance) channel Rij(x,y).

Ri,j(1)(x,y)=Ri,j(x,y)×(1+|A−B|Cr1+|A−B|Cr2)|i=II 

where,

Ai,j(1)(x,y) represents the reinforced achromatic sub-band,
Rj,j(x,y) represents an achromatic sub-band.
|A−B|k represents the contrast value computed around the current point on the chromatic component k in the preferred orientation of the sub-band Rij(x,y), as shown on FIG. 7. In the embodiment, the sets A and B belong to the sub-band of the first crown (low frequency sub-band) of the chromatic component k with an orientation equal to π/4.
In other embodiments, all the sub-bands can be considered.
In step T6, a center/surround suppressive interaction is carried on.
This operation consists first in a step of inhibition/excitation.
A two-dimensional difference-of-Gaussians (DoG) is used to model the non-CRF inhibition behavior of cells. The DoGσxex,σyex,σxinh,σyinh(x,y) is given by the following equation:



  
    
      
        DoG
        
          
            σ
            x
            ex
          
          ,
          
            σ
            y
            ex
          
          ,
          
            σ
            x
            inh
          
          ,
          
            σ
            y
            inh
          
        
      
      ⁡
      
        (
        
          x
          ,
          y
        
        )
      
    
    =
    
      
        
          G
          
            
              σ
              x
              inh
            
            ,
            
              σ
              y
              inh
            
          
        
        ⁡
        
          (
          
            x
            ,
            y
          
          )
        
      
      -
      
        
          G
          
            
              σ
              x
              ex
            
            ,
            
              σ
              y
              ex
            
          
        
        ⁡
        
          (
          
            x
            ,
            y
          
          )
        
      
    
  




  with




  
    
      
        G
        
          
            σ
            x
          
          ,
          
            σ
            y
          
        
      
      ⁡
      
        (
        
          x
          ,
          y
        
        )
      
    
    =
    
      
        1
        
          2
          ⁢
          
            
              π
              ⁡
              
                (
                
                  
                    σ
                    x
                  
                  ⁢
                  
                    σ
                    y
                  
                
                )
              
            
            2
          
        
      
      ⁢
      
        exp
        ⁡
        
          (
          
            
              -
              
                
                  x
                  2
                
                
                  2
                  ⁢
                  
                    σ
                    x
                    2
                  
                
              
            
            -
            
              
                y
                2
              
              
                2
                ⁢
                
                  σ
                  y
                  2
                
              
            
          
          )
        
      
    
  



a two-dimensional gaussian.

Parameters (σxex,σyex) and (σxinh,σyinh) correspond to spatial extends of the Gaussian envelope along the x and y axis of the central Gaussian (the CRF center) and of the inhibitory Gaussian (the surround) respectively. These parameters have been experimentally determined in accordance with the radial frequency of the second crown (the radial frequency ƒε[1.5, 5.7] is expressed in cycles/degree). Finally, the non-classical surround inhibition can be modeled by the normalized weighting function wσxex,σyex,σxinh,σyinh(x,y)given by the following equation:



  
    
      
        w
        
          
            σ
            x
            ex
          
          ,
          
            σ
            y
            ex
          
          ,
          
            σ
            x
            inh
          
          ,
          
            σ
            y
            inh
          
        
      
      ⁡
      
        (
        
          x
          ,
          y
        
        )
      
    
    =
    
      
        1
        
          
            
            
              H
              ⁡
              
                (
                
                  DoG
                  
                    
                      σ
                      x
                      ex
                    
                    ,
                    
                      σ
                      y
                      ex
                    
                    ,
                    
                      σ
                      x
                      inh
                    
                    ,
                    
                      σ
                      y
                      inh
                    
                  
                
                )
              
            
            
          
          1
        
      
      ⁢
      
        H
        ⁡
        
          (
          
            
              DoG
              
                
                  σ
                  x
                  ex
                
                ,
                
                  σ
                  y
                  ex
                
                ,
                
                  σ
                  x
                  inh
                
                ,
                
                  σ
                  y
                  inh
                
              
            
            ⁡
            
              (
              
                
                  x
                  ′
                
                ⁢
                
                  y
                  ′
                
              
              )
            
          
          )
        
      
    
  




  
    with
    ,
    
      

    
    ⁢
    
      
        H
        ⁡
        
          (
          z
          )
        
      
      =
      
        {
        
          
            
              
                
                  0
                  ,
                
              
              
                
                  z
                  <
                  0
                
              
            
            
              
                
                  z
                  ,
                
              
              
                
                  z
                  ≥
                  0
                
              
            
          
          ,
        
      
    
  



(x′,y′) is obtained by translating the original coordinate system by (x0,y0) and rotating it by θi,j expressed in radian,





  
    
      [
      
        
          
            
              x
              ′
            
          
        
        
          
            
              y
              ′
            
          
        
      
      ]
    
    =
    
      
        [
        
          
            
              
                cos
                ⁢
                
                  
                
                ⁢
                
                  θ
                  
                    i
                    ,
                    j
                  
                
              
            
            
              
                sin
                ⁢
                
                  
                
                ⁢
                
                  θ
                  
                    i
                    ,
                    j
                  
                
              
            
          
          
            
              
                
                  -
                  sin
                
                ⁢
                
                  
                
                ⁢
                
                  θ
                  
                    i
                    ,
                    j
                  
                
              
            
            
              
                cos
                ⁢
                
                  
                
                ⁢
                
                  θ
                  
                    i
                    ,
                    j
                  
                
              
            
          
        
        ]
      
      ⁡
      
        [
        
          
            
              
                x
                -
                
                  x
                  0
                
              
            
          
          
            
              
                y
                -
                
                  y
                  0
                
              
            
          
        
        ]
      
    
  
  ,




∥·∥ denotes the L1 norm, i.e the absolute value.
The FIG. 12 shows the structure of non-CRF inhibition.
The FIG. 13 shows a profile example of the normalized weighting function wσxex,σyex,σxinh,σyinh(x,y).
The response Ri,j(2)(x,y) of cortical cells to a particular sub-band Ri,j(1)(x,y) is computed by the convolution of the sub-band Ri,j(1)(x,y) with the weighting function wσxex,σyex,σxinh,σyinh(x,y):

Ri,j(2)(x,y)=H(Ri,j(1)(x,y)−Ri,j(1)(x,y)*wσxex,σyex,σxinh,σyinh(x,y))|i=II

with H(z) defined as has been described above.
In step T7, a facilitative interaction is carried on.
This facilitative interaction is usually termed contour enhancement or contour integration.
Facilitative interactions appear outside the CRF along the preferred orientation axis. These kinds of interactions are maximal when center and surround stimuli are iso-oriented and co-aligned. In other words, as shown by several physiologically observations, the activity of cell is enhanced when the stimuli within the CRF and a stimuli within the surround are linked to form a contour.
Contour integration in early visual preprocessing is simulated using two half butterfly filter Bi,j0 and Bi,j1. The profiles of these filters are shown on the 9 and they are defined by using a bipole/butterfly filter. It consists of a directional term Dθ(x,y) and a proximity term generated by a circle Cr blurred by a gaussian filter Gσx,σy(x,y).



  
    
      
        B
        
          
            θ
            
              i
              ,
              j
            
          
          ,
          α
          ,
          r
          ,
          σ
        
      
      ⁡
      
        (
        
          x
          ,
          y
        
        )
      
    
    =
    
      
        
          
            D
            
              θ
              
                i
                ,
                j
              
            
          
          ⁡
          
            (
            
              x
              ,
              y
            
            )
          
        
        ·
        
          C
          r
        
      
      *
      
        
          G
          
            
              σ
              x
            
            ,
            
              σ
              y
            
          
        
        ⁡
        
          (
          
            x
            ,
            y
          
          )
        
      
    
  




  
    
      with
      ⁢
      
        
      
      ⁢
      
        
          D
          
            θ
            
              i
              ,
              j
            
          
        
        ⁡
        
          (
          
            x
            ,
            y
          
          )
        
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    
                      
                        cos
                        ⁡
                        
                          (
                          
                            
                              
                                π
                                /
                                2
                              
                              α
                            
                            ⁢
                            φ
                          
                          )
                        
                      
                      ⁢
                      si
                      ⁢
                      
                        
                      
                      ⁢
                      φ
                    
                    <
                    α
                  
                
              
              
                
                  
                    0
                    ⁢
                    
                      
                    
                    ⁢
                    
                      sinon
                      .
                    
                  
                
              
            
            ⁢
            
              

            
            ⁢
            and
            ⁢
            
              
            
            ⁢
            φ
          
          =
          
            arc
            ⁢
            
              
            
            ⁢
            
              tan
              ⁡
              
                (
                
                  
                    y
                    ′
                  
                  /
                  
                    x
                    ′
                  
                
                )
              
            
          
        
        ,
      
    
  



where (x′,y′)T is the vector (x,y)T rotated by θi,j. The parameter α defines the opening angle 2α of the bipole filter. It depends on the angular selectivity γ of the considered sub-band. One takes α=0.4×γ. The size of the bipole filter is about twice the size of the CRF of a visual cell.

The two half butterfly filter Bi,j0 and Bi,j1 are after deduced from the butterfly filter by using appropriate windows.
For each orientation, sub-band and location, one computes the facilitative coefficient:




  
    
      f
      
        i
        ,
        j
      
      iso
    
    ⁡
    
      (
      
        x
        ,
        y
      
      )
    
  
  =
  
    D
    ⁡
    
      (
      
        
          
            
              L
              
                i
                ,
                j
              
              1
            
            ⁡
            
              (
              
                x
                ,
                y
              
              )
            
          
          +
          
            
              L
              
                i
                ,
                j
              
              0
            
            ⁡
            
              (
              
                x
                ,
                y
              
              )
            
          
        
        
          max
          ⁡
          
            (
            
              β
              ,
              
                
                
                  
                    
                      L
                      
                        i
                        ,
                        j
                      
                      1
                    
                    ⁡
                    
                      (
                      
                        x
                        ,
                        y
                      
                      )
                    
                  
                  -
                  
                    
                      L
                      
                        i
                        ,
                        j
                      
                      0
                    
                    ⁡
                    
                      (
                      
                        x
                        ,
                        y
                      
                      )
                    
                  
                
                
              
            
            )
          
        
      
      )
    
  




with,

β a constant,

Li,j0(x,y)=Ri,j(2)(x,y)*Bi,j0(x,y),

Li,j(x,y)=Ri,j(2)(x,y)*Bi,j1(x,y),





  
    D
    ⁡
    
      (
      z
      )
    
  
  =
  
    {
    
      
        
          
            
              0
              ,
              
                z
                ≤
                
                  s
                  1
                
              
              ,
            
          
        
        
          
            
              
                α
                1
              
              ,
              
                z
                ≤
                
                  s
                  2
                
              
              ,
            
          
        
        
          
            ⋯
          
        
        
          
            
              
                α
                
                  N
                  -
                  1
                
              
              ,
              
                z
                ≤
                
                  s
                  
                    N
                    -
                    1
                  
                
              
            
          
        
      
      ⁢
      where
    
  




An illustration of the operator D(z) is given on FIG. 9.
The sub-band Ri,j(3) resulting from the facilitative interaction is finally obtained by weighting the sub-band Ri,j(2) by a factor depending on the ratio of the local maximum of the facilitative coefficient fi,jiso(x,y) and the global maximum of the facilitative coefficient computed on all sub-bands belonging to the same range of spatial frequency:




  
    
      
        
          R
          
            i
            ,
            j
          
          
            (
            3
            )
          
        
        ⁡
        
          (
          
            x
            ,
            y
          
          )
        
      
      =
      
        
          
            R
            
              i
              ,
              j
            
            
              (
              2
              )
            
          
          ⁡
          
            (
            
              x
              ,
              y
            
            )
          
        
        ×
        
          (
          
            1
            +
            
              
                η
                iso
              
              ×
              
                
                  
                    max
                    
                      (
                      
                        x
                        ,
                        y
                      
                      )
                    
                  
                  ⁢
                  
                    (
                    
                      
                        f
                        
                          i
                          ,
                          j
                        
                        iso
                      
                      ⁡
                      
                        (
                        
                          x
                          ,
                          y
                        
                        )
                      
                    
                    )
                  
                
                
                  
                    max
                    j
                  
                  ⁢
                  
                    (
                    
                      
                        max
                        
                          (
                          
                            x
                            ,
                            y
                          
                          )
                        
                      
                      ⁢
                      
                        (
                        
                          
                            f
                            
                              i
                              ,
                              j
                            
                            iso
                          
                          ⁡
                          
                            (
                            
                              x
                              ,
                              y
                            
                            )
                          
                        
                        )
                      
                    
                    )
                  
                
              
              ⁢
              
                
                  f
                  
                    i
                    ,
                    j
                  
                  iso
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                  
                  )
                
              
            
          
          )
        
      
    
    
  
  
    i
    =
    II
  




From a standard butterfly shape, this facilitative factor permits to improve the saliency of isolated straight lines. ηiso permits to control the strength of this facilitative interaction.
In step E8, the saliency map is obtained by summing all the resulting sub-bands obtained in step E7.




  
    S
    ⁡
    
      (
      
        x
        ,
        y
      
      )
    
  
  =
  
    
      ∑
      
        
          i
          =
          II
        
        ,
        j
      
    
    ⁢
    
      
        R
        
          i
          ,
          j
        
        
          (
          3
          )
        
      
      ⁡
      
        (
        
          x
          ,
          y
        
        )
      
    
  




In other embodiments of the invention, one can use all the sub-bands and not only the sub-bands of the second crown.
Although cortical cells tuned to horizontal and vertical orientations are almost as numerous as cells tuned to other orientations, we don't introduce any weighting. This feature of the HVS is implicitly mimic by the application of 2D anisotropic CSF.

## Claims

1. In a computer device, a method for creating a saliency map of an image wherein it comprises the steps of:
Projection of said image according to the luminance component and if said image is a color image, according to the luminance component and according to the chrominance components,
Perceptual sub-bands decomposition, by the computer device, of said components according to the visibility threshold of a human eye,
Extraction of the salient elements of the sub-bands related to the luminance component,
Contour enhancement of said salient elements in each sub-band related to the luminance component,
Calculation, by the computer device, of a saliency map from the contour enhancement, for each sub-band related to the luminance component,
Creation of the saliency map as a function of the saliency maps obtained for each sub-band.

2. Method according to claim 1 wherein it comprises, further to the perceptual sub-bands decomposition by the computer device,
a step of achromatic contrast sensitivity function for the luminance component and if said image is a color image,
a step of chromatic sensitivity function for the chromatic components.

3. Method according to claim 2 wherein it comprises a step of visual masking, further to the step of contrast sensitivity function, for each sub-band of the luminance component and of the chromatic components.

4. Method according to claim 3 wherein, when said image is a color image, it comprises a step of chromatic reinforcement of the luminance sub-bands.

5. Method according to claim 1 wherein the perceptual sub-bands decomposition, by the computer, is obtained by carving-up the frequency domain both in spatial radial frequency and orientation.

6. Method according to claim 5 wherein the perceptual decomposition of the luminance component leads to 17 psycho visual sub-bands distributed on four crowns.

7. Method according to claim 5 wherein the perceptual decomposition of the chromatic components leads to 5 psycho visual sub-bands distributed on two crowns for each chromatic component.

8. Method according to claim 4 wherein the chromatic reinforcement of the luminance component is done on the sub-bands of the second crown and based on the sub-bands of the first crown of the chromatic components.

9. A computer readable medium encoded with a computer program, that when executed by a computer devices comprises the steps of:
Projection of said image according to the luminance component and if said image is a color image, according to the luminance component and according to the chrominance components,
Perceptual sub-bands decomposition, by the computer device of said components according to the visibility threshold of a human eye,
Extraction of the salient elements of the sub-bands related to the luminance component,
Contour enhancement of said salient elements in each sub-band related to the luminance component,
Calculation, by the computer device, of a saliency map from the contour enhancement, for each sub-band related to the luminance component,
Creation of the saliency map as a function of the saliency maps obtained for each sub-band.