# Patent

## Title
Adaptive contrast enhancement method for video signals based on time-varying nonlinear transforms

## Abstract

An adaptive contrast enhancement (ACE) method and apparatus provide a natural enhancement in accordance with the time-varying characteristics of a video sequence. Characteristics of a time varying video sequence are specified and a nonlinear transform over the input video sequence is performed to enhance mainly the contrast of the input. A probability density function (PDF) of a time varying input video sequence is computed and then some predetermined video parameters relating to contrast is extracted from the PDF. Based upon the extracted video parameters, a nonlinear transform function is then constructed and updated as a look up table (LUT), which is synchronized with the associated video picture or field SYNC signal. The transform LUT is then applied to the input video to provide an enhanced video output signal.

## Background
FIELD OF THE INVENTION
The present invention relates generally to video processing, and more particularly to video signal enhancement.
BACKGROUND OF THE INVENTION
The development of modern digital video technology has brought significant enhancement in the video quality for consumers, such as in DVD players and in digital TVs (DTV) compared to the analog TV systems. However, such digital video systems only enhance the video quality in terms of signal to noise ratio (SNR) and resolution, without regard to other important issues relating to video enhancement. Such issues include contrast enhancement, brightness enhancement, and detail enhancement. Generally, video enhancement processes comprise a collection of techniques that seek to improve the visual appearance of video when displayed. This primarily includes gray level and contrast manipulation, noise reduction, edge crispening and sharpening. Compared to image restoration, video or image enhancement methods neither increase the inherent information content in the data nor require mathematical modeling. The basic principle of video enhancement is to manipulate a given sequence of images so that their appearance on display media can be improved. Because quantifying the criteria for enhancement is difficult, conventional video enhancement techniques are empirical and require interactive procedures to obtain satisfactory results.
Among the techniques for video enhancement, contrast enhancement is important because it plays a fundamental role in the overall appearance of an image to human being. A human being's perception is sensitive to contrast rather than the absolute values themselves. Hence, it is natural to enhance the contrast of an image in order to provide a good looking image to human beings.
Contrast enhancement involves considering the overall appearance of a given image rather than local appearances such as edge crispening or peaking. There are conventional models of contrast enhancement, and some examples include the root law and the logarithmic law. Image enhancement by contrast manipulation has been performed in various fields of medical image processing, astronomical image processing, satellite image processing, infrared image processing, etc. For example, histogram equalization is a useful method in X-ray image processing because it enhances the details of an X-ray image significantly to e.g. detect tumors easily.
Although several conventional methods for contrast enhancement exist, their primary application is limited to still images. Direct applications of such methods to moving images results in visual artifacts such as unnatural appearance or over-enhancement. Hence, such methods are not suitable for consumer products such as TV whose primary content is a sequence of images. Therefore, there is a need for a contrast enhancement method that is applicable to a video sequence which inherently has time-varying characteristics.
BRIEF SUMMARY OF THE INVENTION
The present invention addresses the above needs. It is an object of the present invention to provide an adaptive contrast enhancement (ACE) method and apparatus which provide a natural enhancement in accordance with the time-varying characteristics of a video sequence.
In one embodiment, an ACE method according to the present invention includes the steps of specifying the characteristics of a time varying video sequence and performing a nonlinear transform over the input video sequence to enhance mainly the contrast of the input. As such, first, a probability density function (PDF) of a time varying input video sequence is computed and then some predetermined video parameters relating to contrast is extracted from the PDF. Based upon the extracted video parameters, a nonlinear transform function is then constructed and updated as a look up table (LUT), which is synchronized with the associated video picture or field SYNC signal. The transform LUT is then applied to the input video to provide an enhanced video output signal.
An example method for adaptive contrast enhancement, according to the present invention includes the steps of: (i) obtaining a time varying video signal including a plurality of temporally ordered digital pictures, each one of the digital pictures represented by a set of samples, each one of the samples having a gradation level in a range from a lower limit C to an upper limit U; (ii) constructing a contrast enhancement transform including at least a first transform function and a second transform function by: for a first one of the digital pictures, selecting a first value for the samples with gradation values between the lower limit C and the upper limit U; dividing the set of the samples representing the first one of the digital pictures into a first portion having samples with a first mean value not greater than the first value and a second portion having samples with a second mean value not less than the first value; selecting the first transform function based on a distribution of the samples with the first mean value in the first portion, wherein the first transform function is based on a first enhancement function that is a varying function of gradation level, equaling zero at the lower limit C, equaling zero at the first value, and having only one local maxima in a range from the lower limit C to the first value; and selecting the second transform function in dependence on a distribution of the samples with the second mean value in the second portion, wherein the second transform function is based on a second enhancement function that is a varying function of gradation level, equaling zero at the first value, equaling zero at the upper limit U, and having only one local maxima in a range from the first value to the upper limit U; and (iii) enhancing contrast of a digital picture by applying the contrast enhancement transform to a set of samples representing the digital picture.
In another embodiment, the present invention provides an adaptive contrast enhancement device for enhancing a time varying video signal including a plurality of temporally ordered digital pictures, each one of the digital pictures represented by a set of samples, each one of the samples having a gradation level in a range from a lower limit C to an upper limit U. The adaptive contrast enhancement device comprises: (i) a distribution estimator that determines the number of samples in a picture, having gradation values; (ii) a mean estimator that: for a first one of the digital pictures, determines a first value for the samples with gradation values between the lower limit C and the upper limit U; for a first portion of the samples, determines a first mean value not greater than the first value, and for a second portion of the samples, determines a second mean value not less than the first value; and (ii) a contrast enhancer that enhances contrast of the digital picture by applying a contrast enhancement transform to a set of samples representing the digital picture, the contrast enhancer including a first transform function and a second transform function, wherein: the first transform function is based on a distribution of the samples with the first mean value in the first portion, and the second transform function in dependence on a distribution of the samples with the second mean value in the second portion.
As such, the present invention provides a contrast enhancement method and device that are applicable to a video sequence which inherently has time-varying characteristics.



BRIEF DESCRIPTION OF THE DRAWINGS
These and other features, aspects and advantages of the present invention will become understood with reference to the following description, appended claims and accompanying figures where:
FIG. 1A is a block diagram of an embodiment of a device for performing the adaptive contrast enhancement method according to the present invention;
FIG. 1B is a block diagram of another embodiment of a device for performing the adaptive contrast enhancement method according to the present invention;
FIG. 2 shows an example of a transform that can be used to increase the dynamic range of samples having gradation levels within a certain range;
FIG. 3 shows an example of a transform that can be used to increase the dynamic range of samples having gradation levels within another range;
FIGS. 4A–C show examples of enhancement functions;
FIG. 5A shows an example of the enhancement functions ƒl(x) and ƒu(x);
FIG. 5B shows a plot of ƒl(x)·gl(ml) and ƒu(x)·gu(mu);
FIG. 5C shows a plot of then transform function that is obtained using the terms ƒl(x)·gl(ml) and ƒu(x)·gu(mu);
FIG. 6 shows an example flowchart of an embodiment of a adaptive contrast enhancement method according to the present invention; and
FIG. 7 shows another block diagram of an embodiment of a device for performing the adaptive contrast enhancement method according to the present invention.



Like reference characters refer to like elements.
DETAILED DESCRIPTION OF THE INVENTION
While this invention is susceptible of embodiments in many different forms, there are shown in the drawings and will herein be described in detail, preferred embodiments of the invention with the understanding that the present disclosure is to be considered as an exemplification of the principles of the invention and is not intended to limit the broad aspects of the invention to the embodiments illustrated.
As noted above, an Adaptive Contrast Enhancement (ACE) method according to the present invention specifies the characteristics of a time varying video sequence and performs a nonlinear transform over the input video sequence to enhance mainly the contrast of the input. Referring to FIG. 1A, an example architecture of an ACE apparatus (device) 10 implementing an example ACE method according to the present invention, is shown. In a functional block 14, a probability density function (PDF) of a time varying input video sequence is computed and then some predetermined video parameters relating to contrast is extracted from the PDF. Based upon the extracted video parameters, a nonlinear transform function is then constructed and updated as a look up table (LUT) in the block 14, which is synchronized with the associated video picture or field SYNC signal. The transform LUT is then applied to the input video to enhance the input signal in a functional block 16. Therefore, the primary steps of the ACE method are the extraction of the video parameters relating to the contrast of the input video signal and the construction of the transform function, described further below.
In the example of FIG. 1A, a memory device 12 is used to delay the input video for one frame/field period to apply the constructed transform to the video frame/field used for the transform construction. The incoming picture is stored in the memory 10 while the transform LUT is constructed using parameters obtained from the picture. The memory 10 is provided to delay the input video for one frame or field period so that the transform LUT can be applied to the picture that was used to construct the transform LUT, as shown in the functional block 14. Alternatively, as shown in another example ACE apparatus 20 of FIG. 1B, the memory device can be removed from the architecture because a video sequence typically has a high correlation in temporal direction, and therefore, in most applications, the LUT transform that is constructed from one picture can be applied to the subsequent picture in the video sequence.
The example ACE method according to the present invention is now described in more detail. In this description, In(·) denotes a picture (fame, or, field) of an incoming video sequence at time instant n, wherein the picture In(·) comprises samples to be enhanced whose values are, in general, from the gradation levels {C,C+1, . . . , U}, where C is associated with the black gradation level and U is associated with the white gradation level in a video system. The determination of the values of C and U is based on each particular application (e.g., C=0, U=255, etc.). The range represented by {C,C+1, . . . , U} can be narrower than the real dynamic range of the input picture system.
The average brightness, or, the mean of the input video picture In(·) is denoted as m, and h(x) denotes the PDF of In(·), where h(x) represents the number of the samples in In(·) whose gradation level equals to x. It is preferable to use the mean m, since good results are obtained with this parameter, however, it should be understood that another value deviating from the mean m could be used instead. Conceptually any one of a number of values between C and U could be used, however, it is preferable to use the mean or a value very close to the mean so that the average brightness will not be changed, or at least will not be changed significantly.
A value




  N
  =
  
    
      ∑
      
        x
        =
        C
      
      U
    
    ⁢
    
      h
      ⁡
      
        (
        x
        )
      
    
  




represents the total number of samples in the input video picture In(·) to be enhanced. Then the mean of those samples can be computed as





  m
  =
  
    
      ∑
      
        x
        =
        C
      
      U
    
    ⁢
    
      x
      ·
      
        
          h
          ⁡
          
            (
            x
            )
          
        
        /
        
          N
          .
        
      
    
  




Or, simply the mean can be computed by summing up the gradation levels of the samples in the input video picture to be enhanced and dividing by N.

Based on the mean, m, two parameters ml and mu are computed which are associated with the mean of samples lower than or equal to the mean m and the mean of samples greater than or equal to the mean m, respectively, wherein:




  
    
      
        
          
            m
            l
          
          =
          
            
              
                ∑
                
                  x
                  =
                  C
                
                m
              
              ⁢
              
                
                  h
                  ⁡
                  
                    (
                    x
                    )
                  
                
                ⁢
                x
              
            
            
              
                ∑
                
                  x
                  =
                  C
                
                m
              
              ⁢
              x
            
          
        
        ⁢
        
          

        
        ⁢
        and
      
    
    
      
        (
        1
        )
      
    
  
  
    
      
        
          m
          u
        
        =
        
          
            
              ∑
              
                x
                =
                m
              
              U
            
            ⁢
            
              
                h
                ⁡
                
                  (
                  x
                  )
                
              
              ⁢
              x
            
          
          
            
              ∑
              
                x
                =
                m
              
              U
            
            ⁢
            x
          
        
      
    
    
      
        (
        2
        )
      
    
  




such that C≦ml≦m and m≦mu≦U. The parameters ml and mu roughly indicate how the corresponding samples are distributed in the regions (C,m) and (m,U), respectively. In one example, the values mu and ml represent the mean brightness of the sub-images (the first is the sub-image which can virtually comprise of the samples less than or equal to the mean, and the second is the sub-image which can virtually comprise of the samples greater than or equal to the mean, respectively). For example, ml≈m implies that the samples in the region (C,m) are mostly distributed near to the mean m, whereas ml≈C implies that the samples in the region (C,m) are mostly distributed near to C. Similarly, mu≈m implies that the samples in the region (m,U) are mostly distributed near to the mean m, whereas mu≈U implies that the samples in the region (m,U) are mostly distributed near to U. The value of m can be selected as desired depending on application (e.g., m=128).
FIG. 2 shows an example representation of a transform, wherein x denotes the input gradation level of the input sample and y denotes the transformed output. In this example, the gradation values in the region (p,k) are mapped to the values in the region (p,k′). Hence, the dynamic range, Din, for the samples whose gradation levels are in the region (p,k) is increased to Dout as a result of the transform. The overall contrast of the image is enhanced if the input picture has more samples in the region (p,k) than in the region (k,q) because the dynamic range for the samples in the region (p,k) has been increased.
Therefore, the example transform illustrated in FIG. 2 is suitable to enhance the contrast when the samples are distributed more in the region (p,k) region than in the region (k,q). However, if more samples are distributed in the region (k,q) than in the region (p,k), then another example transform function shown in FIG. 3 can increase the overall contrast of the picture.
Further, ACE can be performed according to the present invention by combining the example transforms of FIGS. 2–3 in accordance with the sample distribution of the input image. The transform function is adjusted adaptively in accordance with the distribution of the gradation levels in the input picture in order to increase the overall contrast of the picture. An example transform ψ(x) for the ACE according to the present invention is represented by:




  
    
      
        
          ψ
          ⁡
          
            (
            x
            )
          
        
        =
        
          {
          
            
              
                
                  
                  ⁢
                  
                    
                      x
                      +
                      
                        
                          
                            g
                            l
                          
                          ⁡
                          
                            (
                            
                              m
                              l
                            
                            )
                          
                        
                        ·
                        
                          
                            f
                            l
                          
                          ⁡
                          
                            (
                            x
                            )
                          
                        
                      
                    
                    ,
                  
                
              
              
                
                  
                  ⁢
                  
                    C
                    ≤
                    x
                    ≤
                    m
                  
                
              
            
            
              
                
                  
                    x
                    +
                    
                      
                        
                          g
                          u
                        
                        ⁡
                        
                          (
                          
                            m
                            u
                          
                          )
                        
                      
                      ·
                      
                        
                          f
                          u
                        
                        ⁡
                        
                          (
                          x
                          )
                        
                      
                    
                  
                  ,
                
              
              
                
                  m
                  ≤
                  x
                  ≤
                  U
                
              
            
            
              
                
                  
                  ⁢
                  
                    x
                    ,
                  
                
              
              
                
                  
                  ⁢
                  otherwise
                
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




wherein ƒl(x) and ƒu(x) are enhancement basis functions, and gl(ml) and gu(mu) are adaptive gain adjusting functions. The enhancement basis functions ƒl(x) and ƒu(x) determine the general characteristics of the enhancement, whereas the adaptive gain adjusting functions gl(ml) and gu(mu) control the manner and degree of enhancement associated with the input picture.
The following constraints apply to the enhancement basis functions:
ƒl(x) is positive and only defined in region (C,m) and ƒu(x) is positive and only defined in region (m,U),

ƒl(C)=ƒl(m)=0, and ƒu(m)=ƒu(U)=0,  (5)

ƒl(x) has a local maxima in (C,m) and ƒu(x) has a local maxima in (m,U),  (6)

Arbitrary enhancement basis functions can be used so long as they meet the constraints given in relations (4)–(6) above. For example, FIGS. 4A–C show three example plots of the lower enhancement basis function, ƒl(x), satisfying the constraints given in relations (4), (5) and (6) above.
The followings are constraints apply to the adaptive gain adjusting functions:

gl(C)=1, and gl(m)=−1,  (7)

gu(m)=1, and gu(U)=−1,  (8)

gl(ml) is a monotonically decreasing function with respect to mu,  (9)

gu(mu) is a monotonically decreasing function with respect to mu,  (10)

wherein the value of gl(ml) changes from 1 to −1 as ml varies from C to m and the value of gu(mu) changes from 1 to −1 as mu varies from m to U. The gain parameters adjust the enhancement functions in accordance with the sample distribution as discussed in relation to FIGS. 2–3 above. Further, the negative sign of gl(ml) and gu(mu) revert the curvature direction of the enhancement basis functions ƒl(x) and ƒu(x).
If the enhancement basis functions ƒl(x) and ƒu(x) are negative in the regions (C,m) and (m,U), respectively, the conditions given in relations (7)–(10) are changed, respectively, to:

gl(C)=−1, and gl(m)=1,  (11)

gu(m)=−1, and gu(U)=1,  (12)

gl(ml) is a monotonically increasing function with respect to ml,  (13)

gu(mu) is a monotonically increasing function with respect to mu.  (14)

Referring back to relation (3) above, the input gradation levels in (C,U) are changed as gl(ml)·ƒl(x) if x∈(C,m) and as gu(mu)·ƒu(x) if x∈(m,U). The transform given in relation (9) maps m to m because it is required that ƒl(m)=ƒu(m)=0, keeping the mean brightness of the transformed picture without significant change compared to the mean brightness of the input picture.
The characteristics of the transform in relation (3) change from picture to picture as the characteristics or the gradation level distribution of the video signal changes from picture to picture. Specifically, the values of the parameters m, ml and mu vary from picture to picture, and so does the final transform ψ(x) in relation (3). This provides the adaptive nature of an example contrast enhancement method according to the present invention.
Referring to FIG. 5A an example of the enhancement function ƒu(x) which satisfies the conditions given in relations (4)–(6), is illustrated. Further, FIG. 5B shows an example plot of the term ƒu(x)·gu(mu) of relation (3), wherein e.g. gu(mu)=−0.25. And, FIG. 5C is an example of a complete plot of the transform ψ(x), for m≦x≦u, in relation (3) to be applied to the input video.
Referring to FIG. 6 an example flowchart of the steps of the above embodiment of the contrast enhancement method of the present invention is shown. Specific functional forms of ƒl(x), ƒu(x), gl(ml) and gu(mu) that satisfy the respective conditions in relation (4) through (10) predetermined. Generally, the contrast enhancement method based on the transform given in relation (3) further includes the steps of:
Computing the PDF, h(x), of the incoming picture In(·) (step 22);
Computing the mean, m, and the values ml and mu (step 24);
Computing the gain functions gl(ml) and gu(mu) (step 26);
Using the transform in relation (3) to construct the transform LUT (step 27); and
Applying the LUT to the incoming video input signal to generate an enhanced video signal (step 28).
As noted above, the LUT update can be synchronized with a picture SYNC signal, wherein the LUT is applied to transform the input picture that was used to construct the transform and stored in a picture memory (e.g., FIG. 1A). Or, the LUT can be applied to the next input picture if no memory is incorporated (e.g., FIG. 1B). This process is a pixel by pixel operation which outputs ψ(x) for the input pixel gradation level x.
In the above, the example ACE method is based upon the enhancement functions combined with the adaptive gain adjustment, depending on the sample distribution of the incoming picture. Any arbitrary function for ƒl(x), ƒu(x), gl(ml) and gu(mu) can be used as long as they are subject to the conditions given in relations (4)–(10) above.
In another example ACE method according to the present invention, a second order polynomial is used as one specific embodiment of the enhancement functions used in relation (3), wherein:

ƒl(x)=Kl(x−C)(m−x),  (15)

ƒu(x)=Ku(m−x)(x−U),  (16)

where Kl and Ku are pre-determined constants. These enhancement functions satisfy the constraints in relations (4)–(6). Accordingly, the transform in relation (3) an be represented as:




  
    
      
        
          ψ
          ⁡
          
            (
            x
            )
          
        
        =
        
          {
          
            
              
                
                  
                  ⁢
                  
                    
                      x
                      +
                      
                        
                          
                            
                              g
                              l
                            
                            ⁡
                            
                              (
                              
                                m
                                l
                              
                              )
                            
                          
                          ·
                          
                            
                              K
                              l
                            
                            ⁡
                            
                              (
                              
                                x
                                -
                                C
                              
                              )
                            
                          
                        
                        ⁢
                        
                          (
                          
                            m
                            -
                            x
                          
                          )
                        
                      
                    
                    ,
                  
                
              
              
                
                  
                  ⁢
                  
                    C
                    ≤
                    x
                    ≤
                    m
                  
                
              
            
            
              
                
                  
                    x
                    +
                    
                      
                        
                          
                            g
                            u
                          
                          ⁡
                          
                            (
                            
                              m
                              u
                            
                            )
                          
                        
                        ·
                        
                          
                            K
                            u
                          
                          ⁡
                          
                            (
                            
                              m
                              -
                              x
                            
                            )
                          
                        
                      
                      ⁢
                      
                        (
                        
                          x
                          -
                          U
                        
                        )
                      
                    
                  
                  ,
                
              
              
                
                  m
                  ≤
                  x
                  ≤
                  U
                
              
            
            
              
                
                  
                  ⁢
                  
                    x
                    ,
                  
                
              
              
                
                  
                  ⁢
                  otherwise
                
              
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




The choices of the gain functions, gl(ml) and gu(mu), are versatile, satisfying the constraints given in relations (11)–(14). Different gain functions result in different characteristics of ψ(x) and, hence, said gain function choices can be varied depending on specific application. A simple example choice of gl(ml) and gu(mu) can be:

gl(ml)=2(ml−C)/(m−C)−1 and gu(mu)=2(mu−m)/(U−m)−1.  (18)

Direct applications (using constants for Kl and Ku) of the transform given in relation (17) may result in the gray inversion problem, depending on the value of m where the gray inversion implies, wherein ψ(x1)>ψ(x2) for some x1<x2.
Such gray inversion problem can be prevented with an algebraic manipulation such as:




  
    
      ψ
      ⁡
      
        (
        
          x
          1
        
        )
      
    
    <
    
      ψ
      ⁡
      
        (
        
          x
          2
        
        )
      
    
  
  ⁢
  
    ↔
    iff
  
  ⁢
  
    
      x
      1
    
    <
    
      x
      2
    
  




wherein the above values Kl and Ku are bounded such as:




  
    K
    l
  
  ≤
  
    
      1
      
        m
        -
        C
      
    
    ⁢
    
      
    
    ⁢
    and
    ⁢
    
      
    
    ⁢
    
      K
      u
    
  
  ≤
  
    
      1
      
        U
        -
        m
      
    
    .
  




As such, for the boundary values of Kl and Ku, the transform given in relation (17) becomes:




  
    
      
        
          ψ
          ⁡
          
            (
            x
            )
          
        
        =
        
          {
          
            
              
                
                  
                  ⁢
                  
                    
                      x
                      +
                      
                        
                          
                            g
                            l
                          
                          ⁡
                          
                            (
                            
                              m
                              l
                            
                            )
                          
                        
                        ·
                        
                          
                            
                              (
                              
                                x
                                -
                                C
                              
                              )
                            
                            ⁢
                            
                              (
                              
                                m
                                -
                                x
                              
                              )
                            
                          
                          
                            m
                            -
                            C
                          
                        
                      
                    
                    ,
                  
                
              
              
                
                  
                  ⁢
                  
                    C
                    ≤
                    x
                    ≤
                    m
                  
                
              
            
            
              
                
                  
                    x
                    +
                    
                      
                        
                          g
                          u
                        
                        ⁡
                        
                          (
                          
                            m
                            u
                          
                          )
                        
                      
                      ·
                      
                        
                          
                            (
                            
                              m
                              -
                              x
                            
                            )
                          
                          ⁢
                          
                            (
                            
                              x
                              -
                              U
                            
                            )
                          
                        
                        
                          U
                          -
                          m
                        
                      
                    
                  
                  ,
                
              
              
                
                  m
                  ≤
                  x
                  ≤
                  U
                
              
            
            
              
                
                  
                  ⁢
                  
                    x
                    ,
                  
                
              
              
                
                  
                  ⁢
                  otherwise
                
              
            
          
        
      
    
    
      
        (
        19
        )
      
    
  




Furthermore, depending on applications, the transform can be represented by:




  
    
      
        
          ψ
          ⁡
          
            (
            x
            )
          
        
        =
        
          {
          
            
              
                
                  
                  ⁢
                  
                    
                      x
                      +
                      
                        α
                        ·
                        
                          
                            g
                            l
                          
                          ⁡
                          
                            (
                            
                              m
                              l
                            
                            )
                          
                        
                        ·
                        
                          
                            
                              (
                              
                                x
                                -
                                C
                              
                              )
                            
                            ⁢
                            
                              (
                              
                                m
                                -
                                x
                              
                              )
                            
                          
                          
                            m
                            -
                            C
                          
                        
                      
                    
                    ,
                  
                
              
              
                
                  
                  ⁢
                  
                    C
                    ≤
                    x
                    ≤
                    m
                  
                
              
            
            
              
                
                  
                    x
                    +
                    
                      α
                      ·
                      
                        
                          g
                          u
                        
                        ⁡
                        
                          (
                          
                            m
                            u
                          
                          )
                        
                      
                      ·
                      
                        
                          
                            (
                            
                              m
                              -
                              x
                            
                            )
                          
                          ⁢
                          
                            (
                            
                              x
                              -
                              U
                            
                            )
                          
                        
                        
                          U
                          -
                          m
                        
                      
                    
                  
                  ,
                
              
              
                
                  m
                  ≤
                  x
                  ≤
                  U
                
              
            
            
              
                
                  
                  ⁢
                  
                    x
                    ,
                  
                
              
              
                
                  
                  ⁢
                  otherwise
                
              
            
          
        
      
    
    
      
        (
        20
        )
      
    
  




where α is a pre-determined gain to adjust the overall degree of enhancement. Note that ψ(x)=x when α=0 in relation (20), meaning no change.
FIG. 7 shows a block diagram of another example ACE apparatus 30, implementing an example ACE method, according to the present invention. Such an ACE apparatus is for a contrast enhancement method based on the transform in relation (3), wherein specific functional forms of ƒl(x), ƒu(x), gl(ml) and gu(mu) satisfying the respective conditions disclosed in relations (4) through (10) above, are pre-determined. In this example, the ACE apparatus comprises: A memory device 12; A PDF estimator block 32 that computes the PDF, h(x), for an input picture In(·) video signal; A mean compute block 34 that computes the mean, m, and ml and mu values based on the PDF; A transform compute block 36 that computes gl(ml) and gu(mu) and ψ(x); and a LUT transform block 38 that constructs said transform LUT and applies the LUT to the incoming video input signal to generate enhanced video output signal, as described. The update can be synchronized with a picture SYNC signal as shown.
The incoming picture is stored in the memory 10 while the transform LUT is constructed using parameters obtained from the picture. As noted above, the memory 12 is provided to delay the input video for one frame or field period so that the transform LUT can be applied to the picture that was used to construct the transform LUT, as shown in functional block 16. A video sequence typically has a high correlation in the temporal direction, and therefore, in most applications, the LUT transform that is constructed from one picture can be applied to the subsequent picture in the video sequence. As shown in FIG. 1B above, the incoming picture is not stored in a memory while the transform LUT is constructed using the parameters that are obtained from the incoming picture. The transform that had been constructed from the previous picture in the video sequence is applied to this incoming picture. Similarly, the transform that is being constructed from this incoming picture will be applied to the subsequent picture in the video sequence. Applying the transform LUT to the input picture is a pixel by pixel operation that outputs ψ(x) for the input pixel gradation level x. In both embodiments, the LUT transform can be updated in a manner that is synchronized with a picture SYNC signal.
The various components of the arrangements in FIGS. 1A–B and 7 can be implemented in many ways known to those skilled in the art, such as for example, as program instructions for execution by a processor, as logic circuits such as ASIC, etc. The present invention has been described in considerable detail with reference to certain preferred versions thereof; however, other versions are possible. Therefore, the spirit and scope of the appended claims should not be limited to the description of the preferred versions contained herein.

## Claims

1. A method for adaptive contrast enhancement, comprising the steps of:
obtaining a time varying video signal including a plurality of temporally ordered digital pictures, each one of the digital pictures represented by a set of samples, each one of the samples having a gradation level in a range from a lower limit C to an upper limit U;
constructing a contrast enhancement transform including at least a first transform function and a second transform function by performing steps of:
(a) for a first one of the digital pictures, selecting a first value for the samples with gradation values between the lower limit C and the upper limit U,
(b) dividing the set of the samples representing the first one of the digital pictures into a first portion having samples with a first distribution value not greater than the first value and a second portion having samples with a second distribution value not less than the first value,
(c) selecting the first transform function based on a distribution of the samples with the first distribution value in the first portion, and
(d) selecting the second transform function in dependence on a distribution of the samples with the second distribution value in the second portion; and

enhancing contrast of a digital picture by applying the contrast enhancement transform to a set of samples representing the digital picture.

2. The method of claim 1, wherein step (a) further includes the steps of:
calculating a mean value m from the gradation level of all of the samples representing the first one of the digital pictures, and selecting the first value to be the mean value m.

3. The method of claim 1, wherein step (b) further includes the steps of:
calculating a first mean value ml from the gradation level of the samples in the first portion and selecting the first distribution value to be the first mean value ml, and calculating a second mean value mu from the gradation level of the samples in the second portion and selecting the second distribution value to be the second mean value mu.

4. The method of claim 3, wherein:
the mean value ml represents the mean brightness of the samples in the first portion, and the mean value mu represents the mean brightness of the samples in the second portion.

5. The method of claim 3, further comprising the steps of:
selecting the digital picture, which is enhanced when performing the step of enhancing the contrast, from a set of digital pictures including the first one of the digital pictures and one of the digital pictures that is temporally subsequent with respect to the first one of the digital pictures.

6. The method of claim 3, wherein the digital picture that is enhanced when performing the step of enhancing the contrast is an immediately temporally subsequent picture with respect to the first one of the digital pictures.

7. The method of claim 3, further comprising the steps of:
using the contrast enhancement transform to construct a look-up table for receiving input gradation levels in a range from the lower limit C to the upper limit U and for providing corresponding output gradation levels having an increased dynamic range in at least two regions as compared to corresponding regions of the input gradation levels; and
performing the step of enhancing the contrast by applying the look-up table to the set of the samples and thereby inherently applying the contrast enhancement transform to the set of the samples.

8. The method of claim 7, wherein the step of constructing the contrast enhancement transform includes the steps of:
constructing the first transform function based on a first enhancement function that is a varying function of gradation level, equaling zero at the lower limit C, equaling zero at the first value, and having only one local maxima in a range from the lower limit C to the first value; and
constructing the second transform function based on a second enhancement function that is a varying function of gradation level, equaling zero at the first value, equaling zero at the upper limit U, and having only one local maxima in a range from the first value to the upper limit U.

9. A method for adaptive contrast enhancement, comprising the steps of:
obtaining a time varying video signal including a plurality of temporally ordered digital pictures, each one of the digital pictures represented by a set of samples, each one of the samples having a gradation level in a range from a lower limit C to an upper limit U;
constructing a contrast enhancement transform including at least a first transform function and a second transform function by performing steps of:
(a) for a first one of the digital pictures, selecting a first value for the samples with gradation values between the lower limit C and the upper limit U,
(b) dividing the set of the samples representing the first one of the digital pictures into a first portion having samples with a first mean value not greater than the first value and a second portion having samples with a second mean value not less than the first value,
(c) selecting the first transform function based on a distribution of the samples with the first mean value in the first portion, wherein the first transform function is based on a first enhancement function that is a varying function of gradation level, equaling zero at the lower limit C, equaling zero at the first value, and having only one local maxima in a range from the lower limit C to the first value, and
(d) selecting the second transform function in dependence on a distribution of the samples with the second mean value in the second portion, wherein the second transform function is based on a second enhancement function that is a varying function of gradation level, equaling zero at the first value, equaling zero at the upper limit U, and having only one local maxima in a range from the first value to the upper limit U; and

enhancing contrast of a digital picture by applying the contrast enhancement transform to a set of samples representing the digital picture.

10. The method of claim 9, wherein the step of constructing the contrast enhancement transform further includes the steps of:
choosing the first enhancement function to be non-negative in a region from the lower limit C to the first value;
selecting a first gain function that monotonically decreases as a function of the first mean value, the first gain function equaling one at the lower limit C and equaling minus one at the first value;
choosing the second enhancement function to be non-negative in a region from the first value to the upper limit U; and
selecting a second gain function that monotonically decreases as a function of the second mean value, the second gain function equaling one at the first value and equaling minus one at the upper limit U;
constructing the first transform function by multiplying the first enhancement function by the first gain function; and
constructing the second transform function by multiplying the second enhancement function by the second gain function.

11. The method of claim 10, further comprising the steps of:
calculating a mean value m from the gradation level of all of the samples representing the first one of the digital pictures, and selecting the first value to be the mean value m;
calculating a mean value ml from the gradation level of the samples in the first portion and selecting the first mean value to be the mean value ml; and
calculating a mean value mu from the gradation level of the samples in the second portion and selecting the second mean value to be the mean value mu.

12. The method of claim 11, wherein:
the mean value ml represents the mean brightness of the samples in the first portion, and the mean value mu represents the mean brightness of the samples in the second portion.

13. The method of claim 9, wherein the step of constructing the contrast enhancement transform further includes the steps of:
choosing the first enhancement function to be non-positive in a region from the lower limit C to the first value;
selecting a first gain function that monotonically increases as a function of the first mean value, the first gain function equaling minus one at the lower limit C and equaling one at the first value;
choosing the second enhancement function to be non-positive in a region from the first value to the upper limit U; and
selecting a second gain function that monotonically increases as a function of the second mean value, the second gain function equaling minus one at the first value and equaling one at the upper limit U;
constructing the first transform function by multiplying the first enhancement function by the first gain function; and
constructing the second transform function by multiplying the second enhancement function by the second gain function.

14. The method of claim 13, further comprising the steps of:
calculating a mean value m from the gradation level of all of the samples representing the first one of the digital pictures, and selecting the first value to be the mean value m;
calculating a mean value ml from the gradation level of the samples in the first portion and selecting the first mean value to be the mean value ml; and calculating a mean value mu from the gradation level of the samples in the second portion and selecting the second mean value to be the mean value mu.

15. The method of claim 14, wherein:
the mean value ml represents the mean brightness of the samples in the first portion, and the mean value mu represents the mean brightness of the samples in the second portion.

16. The method of claim 9, wherein the step of constructing the contrast enhancement transform includes the steps of:
when the first enhancement function is chosen to be non-negative in a region from the lower limit C to the first value, then selecting a first gain function that monotonically decreases as a function of the first mean value;
when the second enhancement function is chosen to be non-negative in a region from the first value to the upper limit U, then selecting a second gain function that monotonically decreases as a function of the second mean value;
when the first enhancement function is chosen to be non-positive in the region from the lower limit C to the first value, then selecting a first gain function that monotonically increases as a function of the first mean value;
when the second enhancement function is chosen to be non-positive in the region from the first value to the upper limit U, then selecting a second gain function that monotonically increases as a function of the second mean value;
constructing the first transform function by multiplying the first enhancement function by the first gain function; and
constructing the second transform function by multiplying the second enhancement function by the second gain function.

17. The method of claim 16, further comprising the steps of:
calculating a mean value m from the gradation level of all of the samples representing the first one of the digital pictures, and selecting the first value to be the mean value m;
calculating a mean value ml from the gradation level of the samples in the first portion and selecting the first mean value to be the mean value ml; and
calculating a mean value mu from the gradation level of the samples in the second portion and selecting the second mean value to be the mean value mu.

18. The method of claim 17, wherein:
the mean value ml represents the mean brightness of the samples in the first portion, and the mean value mu represents the mean brightness of the samples in the second portion.

19. The method of claim 17, further comprising the steps of:
selecting the first enhancement function to be Kl(x−C)(m−x), where Kl is a constant and x is an input sample; and
selecting the second enhancement function to be Ku(m−x)(x−U), where Ku is a constant.

20. The method of claim 19, further comprising the steps of:
choosing Kl to be not greater than 1/(m−C); and
choosing Ku to be not greater than 1/(U−m).

21. The method of claim 19, further including the steps of:
selecting the first gain function to be:


2(ml−C)/(m−C)−1; and


selecting the second gain function to be:


2(mu−m)/(u−m)−1.

22. The method of claim 21, further including the steps of:
choosing Kl to be 1/(m−C); and
choosing Ku to be 1/(U−m).

23. The method of claim 16, further comprising the steps of:
constructing the first transform function by multiplying the first enhancement function by the first gain function and by a first pre-determined gain factor to adjust the overall degree of enhancement; and
constructing the second transform function by multiplying the second enhancement function by the second gain function and by a second pre-determined gain factor to adjust the overall degree of enhancement.

24. An adaptive contrast enhancement device for enhancing a time varying video signal including a plurality of temporally ordered digital pictures, each one of the digital pictures represented by a set of samples, each one of the samples having a gradation level in a range from a lower limit C to an upper limit U, comprising:
a distribution estimator that determines the number of samples in a picture, having gradation values
a mean estimator that:
for a first one of the digital pictures, determines a first value for the samples with gradation values between the lower limit C and the upper limit U;
for a first portion of the samples, determines a first mean value not greater than the first value, and
for a second portion of the samples, determines a second mean value not less than the first value,

a contrast enhancer that enhances contrast of the digital picture by applying a contrast enhancement transform to a set of samples representing the digital picture, the contrast enhancer including a first transform function and a second transform function, wherein:
the first transform function is based on a distribution of the samples with the first mean value in the first portion, and
the second transform function in dependence on a distribution of the samples with the second mean value in the second portion.

25. The adaptive contrast enhancement device of claim 24, wherein:
the first transform function is based on a first enhancement function that is a varying function of gradation level, equaling zero at the lower limit C, equaling zero at the first value, and having only one local maxima in a range from the lower limit C to the first value, and
the second transform function is based on a second enhancement function that is a varying function of gradation level, equaling zero at the first value, equaling zero at the upper limit U, and having only one local maxima in a range from the first value to the upper limit U.

26. The adaptive contrast enhancement device of claim 25, wherein the mean estimator:
calculates a mean value m from the gradation level of all of the samples representing the first one of the digital pictures, and selects the first value to be the mean value m;
calculates a mean value ml from the gradation level of the samples in the first portion and selects the first mean value to be the mean value ml; and
calculates a mean value mu from the gradation level of the samples in the second portion and selects the second mean value to be the mean value mu.

27. The adaptive contrast enhancement device of claim 26, wherein:
the mean value ml represents the mean brightness of the samples in the first portion, and the mean value mu represents the mean brightness of the samples in the second portion.

28. The adaptive contrast enhancement device of claim 25, wherein:
the contrast enhancer further includes a first gain function and a second gain function, such that:
when the first enhancement function is non-negative in a region from the lower limit C to the first value, the first gain function monotonically decreases as a function of the first mean value;
when the second enhancement function is non-negative in a region from the first value to the upper limit U, then the second gain function monotonically decreases as a function of the second mean value;
when the first enhancement function is non-positive in the region from the lower limit C to the first value, then the first gain function monotonically increases as a function of the first mean value;
when the second enhancement function is non-positive in the region from the first value to the upper limit U, then the second gain function monotonically increases as a function of the second mean value;
the first transform function is based on a mathematical product of the first enhancement function and the first gain function; and
the first transform function is based on a mathematical product of the second enhancement function and the second gain function.

29. The adaptive contrast enhancement device of claim 28, the mean estimator:
calculates a mean value m from the gradation level of all of the samples representing the first one of the digital pictures, and selects the first value to be the mean value m;
calculates a mean value ml from the gradation level of the samples in the first portion and selects the first mean value to be the mean value ml; and
calculates a mean value mu from the gradation level of the samples in the second portion and selects the second mean value to be the mean value mu.

30. The adaptive contrast enhancement device of claim 29, wherein:
the first enhancement function is selected to be Kl(x−C)(m−x), where Kl is a constant and x is an input sample; and
the second enhancement function is selected to be Ku(m−x)(x−U), where Ku is a constant.

31. The adaptive contrast enhancement device of claim 30, wherein:
Kl is selected to be not greater than 1/(m−C); and
Ku is selected to be not greater than 1/(U−m).

32. The adaptive contrast enhancement device of claim 30, wherein:
the first gain function is selected to be:


2(ml−C)/(m−C)−1; and


the second gain function is selected to be:


2(mu−m)/(u−m)−1.

33. The adaptive contrast enhancement device claim 32,wherein:
Kl is selected to be 1/(m−C); and
Ku is selected to be 1/(U−m).

34. The adaptive contrast enhancement device of claim 28, wherein:
the first transform function is based on a mathematical product of: (i) the first enhancement function, (ii) the first gain function and (iii) a first pre-determined gain factor to adjust the overall degree of enhancement; and
the first transform function is based on a mathematical product of: (i) the second enhancement function, (ii) the second gain function and (iii) a second pre-determined gain factor to adjust the overall degree of enhancement.