# Patent

## Title
Methods and systems for noise reduction and image enhancement involving selection of noise-control parameter

## Abstract

Aspects of the present invention are related to systems and methods for image enhancement, wherein an input image may be separated into two images: a first image comprising image content, from the input image, which may be sharpened; and a second image comprising image content, from the input image, which may be attenuated. One aspect of the present invention relates to determining a value for a noise-control parameter that controls the amount of attenuation applied to the second image.

## Background
FIELD OF THE INVENTION
The present invention relates to image and video processing and, in particular, to methods and systems for noise reduction and image enhancement in an image or a video sequence.
BACKGROUND
The quality of a video sequence or of an image may vary widely depending on the source. For example, computer-generated content may be of the highest quality. Packaged media content, for example, Blu-ray media, may be of relatively high quality. Over-the-air high-definition (HD) broadcast content may be of relatively intermediate quality, while content distributed via cable and satellite, may be of a relatively lower quality. Internet protocol television (IPTV) and streamed content may be of relatively low quality. Methods and systems for video and image enhancement that automatically adapt to image- or video-content quality may be desirable.
SUMMARY
Some embodiments of the present invention comprise methods and systems for noise reduction and image enhancement in an input image or a video sequence.
According to a first aspect of the present invention, an edge-preserving filter, for generating, from an input image, a first image comprising image content to sharpen and a second image comprising image content to attenuate, may be controlled based on a noise estimate associated with the input image. In some of these embodiments, a parameter value associated with the edge-preserving filter may be proportional to the noise estimate. In some embodiments, the edge-preserving filter may comprise a bi-lateral filter.
In some embodiments of the present invention, the noise estimate may be based on a first noise estimate associated with, at least one of, the sensor, thermal and grain noise in the input image and a second noise estimate associated with the compression noise in the input image. In some embodiments, the compression noise in the input image may be determined from coding parameters associated with the codec (coder/decoder) by which the image was compressed. In alternative embodiments, the compression noise in the input image may be determined from the image data.
According to a second aspect of the present invention, an edge-preserving filter, for generating, from an input image, a first image comprising image content to sharpen and a second image comprising image content to attenuate, may be controlled based on a noise estimate associated with the input image and a level of sharpening. In some embodiments, the edge-preserving filter may comprise a bi-lateral filter.
In some embodiments of the present invention, the noise estimate may be adjusted to account for the level of sharpening. In some of these embodiments, the level of sharpening may be associated with an alpha parameter associated with an un-sharp masking filter.
According to a third aspect of the present invention, a noise-control-parameter value that controls the amount of attenuation applied to the second image may be automatically selected based on the amount of compression detected in the input image.
The foregoing and other objectives, features, and advantages of the invention will be more readily understood upon consideration of the following detailed description of the invention taken in conjunction with the accompanying drawings.



BRIEF DESCRIPTION OF THE SEVERAL DRAWINGS
FIG. 1 is a chart showing exemplary embodiments of the present invention comprising controlling the image separation effectuated by an edge-preserving filter using an estimate of the noise in the input image;
FIG. 2 is a chart showing exemplary embodiments of the present invention comprising estimation of compression noise by accessing a look-up table;
FIG. 3 is an exemplary, original, decoded image frame;
FIG. 4 is an enhanced image resulting from processing, according to embodiments of the present invention, of the exemplary, original, decoded image frame shown in FIG. 3, wherein the noise estimate used to control the bi-lateral filter parameter value was based only on sensor, thermal and grain noise;
FIG. 5 is an enhanced image resulting from processing, according to embodiments of the present invention, of the exemplary, original, decoded image frame shown in FIG. 3, wherein the noise estimate used to control the bi-lateral filter parameter value was based on sensor, thermal and grain noise and compression noise;
FIG. 6 is an exemplary, original, decoded image frame associated with a lightly compressed image;
FIG. 7 is an enhanced image resulting from processing, according to embodiments of the present invention, of the exemplary, original decoded image frame shown in FIG. 6, wherein the noise estimate used to control the bi-lateral filter parameter value was based only on compression noise;
FIG. 8 is an enhanced image resulting from processing, according to embodiments of the present invention, of the exemplary, original decoded image frame shown in FIG. 6, wherein the noise estimate used to control the bi-lateral filter parameter value was based on noise due to sensor, thermal and grain noise and compression noise;
FIG. 9 is a chart showing exemplary embodiments of the present invention comprising estimation of compression noise directly from image data;
FIG. 10 is a chart showing exemplary embodiments of the present invention comprising estimation of compression noise directly from image data for a still image or a single frame in a video sequence;
FIG. 11 is a picture illustrating exemplary pixel locations for determination of horizontal differences according to embodiments of the present invention;
FIG. 12 is a chart showing exemplary embodiments of the present invention comprising estimation of compression noise directly from image data for a frame in a video sequence;
FIG. 13 is chart showing exemplary embodiments of the present invention comprising un-sharp masking;
FIG. 14 is a chart showing exemplary embodiments of the present invention comprising edge-preserving filter parameter value determination based on the level of sharpening associated with the un-sharp masking;
FIG. 15 is a chart showing exemplary embodiments of the present invention comprising edge-preserving filter parameter value determination based on the level of sharpening associated with the un-sharp masking and based on temporal smoothing;
FIG. 16 is a plot showing two exemplary look-up tables for associating a sharpening factor with a multiplier;
FIG. 17 is a chart showing exemplary embodiments of the present invention comprising edge-preserving filter parameter value determination based on the level of sharpening associated with the un-sharp masking;
FIG. 18 is a chart showing an exemplary embodiment of the present invention comprising edge-preserving filter parameter value determination based on the level of sharpening associated with the un-sharp masking;
FIG. 19 is an exemplary, original noise image;
FIG. 20 shows a result of bi-lateral filtering, using a sigma range parameter value equal to the standard deviation of the noise, of the original noise image shown in FIG. 19;
FIG. 21 shows a result of bi-lateral filtering, according to embodiments of the present invention, of the original noise image shown in FIG. 19, where the bi-lateral filter parameter value was set based on the noise in the original noise image;
FIG. 22 shows a result of bi-lateral filtering, according to embodiments of the present invention, of the original noise image shown in FIG. 19, where the bi-lateral filter parameter value was set based on the noise in the original noise image and the level of sharpening;
FIG. 23 shows an exemplary, original image;
FIG. 24 shows a result of bi-lateral filtering, according to embodiments of the present invention, of the original image shown in FIG. 23, where the bi-lateral filter parameter value was set based on noise estimated from the original image;
FIG. 25 shows a result of bi-lateral filtering, according to embodiments of the present invention, of the original image shown in FIG. 23, where the bi-lateral filter parameter value was set based on noise estimated from the original image and the level of sharpening;
FIG. 26 is a picture depicting exemplary embodiments of the present invention comprising automatic selection of a noise-control parameter for controlling the amount of attenuation of an image channel;
FIG. 27 is a picture depicting exemplary embodiments of the present invention comprising automatic selection of a noise-control parameter for controlling the amount of attenuation of an image channel based on a noise estimate of the image noise in an input image and a noise estimate of the compression noise in the input image;
FIG. 28 is an exemplary mapping between the compression detected in an input image and a noise-control parameter according to embodiments of the present invention;
FIG. 29 is a picture depicting exemplary embodiments of the present invention comprising automatic selection of a noise-control parameter for controlling the amount of attenuation of an image channel and edge-preserving filter parameter value determination based on the level of sharpening associated with the un-sharp masking; and
FIG. 30 is a picture depicting exemplary embodiments of the present invention comprising automatic selection of a noise-control parameter for controlling the amount of attenuation of an image channel based on a noise estimate of the image noise in an input image and a noise estimate of the compression noise in the input image and edge-preserving filter parameter value determination based on the level of sharpening associated with the un-sharp masking and based on temporal smoothing.



DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS
Embodiments of the present invention will be best understood by reference to the drawings, wherein like parts are designated by like numerals throughout. The figures listed above are expressly incorporated as part of this detailed description.
It will be readily understood that the components of the present invention, as generally described and illustrated in the figures herein, could be arranged and designed in a wide variety of different configurations. Thus, the following more detailed description of the embodiments of the methods and systems of the present invention is not intended to limit the scope of the invention but it is merely representative of the presently preferred embodiments of the invention.
Elements of embodiments of the present invention may be embodied in hardware, firmware and/or software. While exemplary embodiments revealed herein may only describe one of these forms, it is to be understood that one skilled in the art would be able to effectuate these elements in any of these forms while resting within the scope of the present invention.
The quality of a video sequence or of an image may vary widely depending on the source. For example, computer-generated content may be of the highest quality, packaged media content, for example, Blu-ray media, may be of relatively high quality, over-the-air high-definition (HD) broadcast content may be of relatively intermediate quality, while content distributed via cable and satellite, may be of a relatively lower quality and Internet protocol television (IPTV) and streamed content may be of relatively low quality. Methods and systems for video and image enhancement that automatically adapt to image- or video-content quality may be desirable.
When there is little, or no, noise in an image, either a still image or a frame in a video sequence, the image may be sharpened aggressively, and small details in the image may not be attenuated. However, selective image sharpening may be required in the presence of noise. Stronger and larger-scale features, as defined relative to the degree of noise, may be identified and enhanced, while less significant and smaller-scale features may be unprocessed, or potentially attenuated to mitigate the noise present in the image.
U.S. patent application Ser. No. 12/228,774, entitled “Image Sharpening Technique,” filed Aug. 15, 2008, which is hereby incorporated by reference herein in its entirety, describes a technique for image sharpening wherein two image channels are generated from an image: a first channel which may include primarily texture information and a second channel which may include primarily edge and flat-region information. The first channel may be filtered to attenuate higher frequency content of the image in the first channel, and the second channel may be filtered to sharpen the image in the second channel. The filtered first channel and the filtered second channel may be combined to form an enhanced image associated with the input image.
In some embodiments of the present invention, an edge-preserving filter may be used to separate an input image into two channels: a first channel which may comprise image content to be sharpened; and a second channel which may comprise image content to attenuate. In some embodiments of the present invention, the edge-preserving filter may comprise a bi-lateral filter. In alternative embodiments, another edge-preserving filter may be used to perform the separation. In some embodiments, the assignment of the input-image content to the first channel and the second channel may be controlled by a parameter of the edge-preserving filter. For example, a range sigma parameter of a bi-lateral filter may control the assignment of input-image content to the first channel and the second channel in embodiments wherein the edge-preserving filter comprises a bi-lateral filter. In some embodiments, in order to maximize the amount of sharpening, the filter parameter value, for example, the range sigma parameter value in a bi-lateral filter, may be set based on noise statistics associated with the input image.
Some embodiments of the present invention may be understood in relation to FIG. 1. An input image 2 and a noise estimate 3 associated with the input image 2 may be made available to an edge-preserving filter 4. In some embodiments of the present invention, the edge-preserving filter may comprise a bi-lateral filter. In alternative embodiments, another edge-preserving filter may be used to perform the separation. The edge-preserving filter 4 may be used to divide the input image 2 into two channels 6, 10. A first channel 6 may correspond to significant features in the input image 2, and the second channel 10, formed by removing 8 the first channel 6 from the input image 2, may contain the residual difference between the first channel 6 and the input image 2. The residual difference may comprise noise and texture data. The first channel 6 may be sharpened 12, and the second channel 10 may be attenuated 14. The sharpened channel 16 and the attenuated channel 18 may be combined 20 to form an enhanced image 22. The value of the filter parameter of the edge-preserving filter 4 may control the assignment of input-image content to the first channel 6 and the second channel 10. The value of the filter parameter of the edge-preserving filter 4 may be based on the noise estimate 3. In some embodiments of the present invention, the value of the filter parameter may be proportional to the noise estimate 3. However, the value of the filter parameter may not be equal to zero. In some embodiments, when the noise estimate 3 indicates that the filter parameter should be set to zero, a small, pre-determined value may be assigned to the filter parameter. In alternative embodiments, when the noise estimate 3 indicates that the filter parameter should be set to zero, the edge-preserving filter 4 may pass the input image 2 directly through as the first channel 6, thereby effectuating only sharpened data. In some embodiments of the present invention (not shown), the enhanced image 22 may be up-sampled to a resolution greater than that of the input image 2.
Multiple noise processes may appear in an image or video sequence. At one extreme, an original image may be of very high quality, for example, computer generated imagery. However, other images, or video sequences, may be of lower quality, for example, a sequence or image may be corrupted, during the acquisition process by thermal, sensor, or other noise. In some situations, film grain, or other analog noise, may corrupt an image or video sequence. Furthermore, compression artifacts may corrupt an image or video sequence. The degree of compression artifacts may be related to the bit-rate associated with a video sequence or image. Exemplary bit-rates and quality levels may be high quality Blu-ray discs compressed at approximately 40 Mbps (Mega bits per second), lower quality over-the-air transmissions compressed at approximately 20 Mbps, further lower quality trans-coded cable transmissions compressed at approximately 12 Mbps and lowest quality satellite and IPTV services compressed at less than 10 Mbps.
In some embodiments of the present invention described in relation to FIG. 2, an input bitstream 30, associated with a still image or a video sequence, may be decoded by an appropriate decoder 32. The image data 34 may be made available, to an image-noise statistic estimator 38, from the decoder 32, and the decoded coding parameters 36 may be made available, to a codec (coder/decoder)-noise statistic estimator 40, from the decoder 32. In alternative embodiments, the coding parameters 36 may be made available to the codec-noise statistic estimator 40 through meta-data or other external means.
The image-noise statistic estimator 38 may analyze the image data 34 to determine the amount of thermal, sensor, grain, or other image noise present in the image data 34. Many existing methods are known in the art for estimating the image noise present in image data. In an exemplary embodiment of the present invention, the image-noise statistic estimator 38 may identify one, or more, substantially smooth regions in the image data 34 and may calculate the standard deviation of the pixel values within each region. The maximum standard deviation value may be associated with the image-noise statistic 44.
The codec-noise statistic estimator 40 may receive coding parameters 36 associated with the input bitstream 30. In some embodiments, the coding parameters 36 may comprise the quantization interval used for coding the current slice or frame in a video sequence. In alternative embodiments, the coding parameters 36 may comprise the number of bits used to represent the source data. The coding parameters 36 may be used in a look-up operation to determine the amount of noise due to compression. In some embodiments of the present invention, the codec type may be used to select a look-up table from multiple, pre-computed look-up tables 42. The coding parameters 36 may be used as an index into the selected table, and the output may be a measure of the image noise due to compression 46. The measurement provided by the look-up tables 42 may be normalized, either at the time the look-up tables 42 are created or in a post-processing operation (not shown) to have the same units of measurement as the noise estimate 44 provided by the image-noise statistic estimator 38.
In some embodiments of the present invention, a look-up table 42 may be constructed by selecting multiple images and video sequences that are representative of input data. The images and videos sequences may be compressed using a variety of codecs and codec settings, and each result may be subsequently decoded. The decoded results may be compared to the original data, and the noise due to the compression system may be computed, for example. The standard deviation of the error. This computation may be performed over all images and video sequences that are compressed using the same codec and parameter settings, and the result may be stored in the look-up table as the noise statistic for the combination of codec and codec parameter settings. In some embodiments, if additional values are needed, the test image and video sequence data may be compressed with the desired configuration settings. In alternative embodiments, noise statistics may be interpolated from values with a similar compression configuration.
The image-noise statistic 44 and the codec-noise statistic 46 may be combined by a combiner 48 to produce a final noise estimate 50, which may be made available to the edge-preserving filter to control the filter parameter. In one embodiment of the present invention, the maximum of the image-noise statistic 44 and the codec-noise statistic 46 may be assigned to the noise estimate 50. Alternative fusion methods for combining 48 the image-noise statistic 44 and the codec-noise statistic 46 may be used to produce a final noise estimate 50 which may be made available to the edge-preserving filter, and the value of the filter parameter of the edge-preserving filter may be based on the noise estimate. In some embodiments of the present invention, the value of the filter parameter may be proportional to the noise estimate.
The effectiveness of some embodiments of the present invention may be illustrated in FIGS. 3-8. FIG. 3 depicts an exemplary, original, decoded image 60 with two regions 62, 64 that exhibit coding artifacts shown inside two white circles superimposed on the image 60. One region 62 is along the road, and another region 64 is around a car. FIG. 4 shows an enhancement 70 of the decoded image 60 shown in FIG. 3. The enhanced image 70 was generated using a bi-lateral filter using only an estimate of the image noise to select the range sigma parameter value of the bi-lateral filter. The coding artifacts are visibly enhanced, for example, the regions 72, 74 that are shown in the two superimposed, white circles. However, the enhanced image 80 shown in FIG. 5, by contrast, was enhanced according to embodiments of the present invention wherein a noise estimate accounting for both the image noise, for example, sensor, thermal, grain and other image noise, and the compression noise in the image frame was used to select the range sigma parameter value of the bi-lateral filter. Using this noise estimate produces an enhanced image 80, but without amplifying the coding artifacts, for example, again examine the regions 82, 84 inside the superimposed, white circles.
A second exemplary, original, decoded image 90 is depicted in FIG. 6. Here, the image 90 is the result of light compression. Inspection of FIG. 6 shows a significant amount of grain noise, which may become much more visible if enhanced. Shown in FIG. 7 is an enhanced image 100 wherein the noise estimate controlling the parameter selection in the bi-lateral filter considered only the coding noise. As can be seen in FIG. 7, the enhanced image 100 contains amplified grain noise. However, the enhanced image 110 shown in FIG. 8 resulting from image enhancement according to embodiments of the present invention wherein the noise estimate used to select the range sigma parameter for the bi-lateral filter accounted for both the image noise, for example, sensor, thermal, grain and other noise, and the compression noise in an image frame. As can be seen from FIG. 8, the resulting image is enhanced but without amplifying grain noise artifacts.
Embodiments of the present invention described in relation to FIG. 2 comprise codec-noise statistic estimation from the coding parameters. In alternative embodiments described in relation to FIG. 9, the compression noise may be estimated from the decoded image data.
In some embodiments of the present invention, described in relation to FIG. 9, comprising a noise-estimation system 115, an input bitstream 120, associated with a still image or a video sequence, may be decoded by a decoder 122. The image data 124 may be made available, from the decoder 122, to an image-noise statistic estimator 126 and to a codec-noise statistic estimator 128. The image-noise statistic estimator 126 may estimate an image-noise statistic 130 associated with the amount of image noise present in the image data 124, and the codec-noise statistic estimator 128 may estimate a codec-noise statistic 132 associated with the compression noise. The image-noise statistic 130 and the codec-noise statistic 132 may be combined by a combiner 134 to produce a final noise estimate 136 associated with the image data. In one embodiment of the present invention, the maximum value of the image-noise statistic 130 and the codec-noise statistic 132 may be assigned to the noise estimate 136. Alternative fusion methods for combining 134 the image-noise statistic 130 and the codec-noise statistic 132 may be used to produce a final noise estimate 136.
The image-noise statistic estimator 126 may analyze the image data 124 to determine the amount of thermal, sensor, grain and other image noise present in the image data 124. Many existing methods are known in the art for estimating the image noise. In an exemplary embodiment of the present invention, the image-noise statistic estimator 126 may identify one, or more, substantially smooth regions in the image data 124 and may calculate the standard deviation of the pixel values within each region. The maximum standard-deviation value may be associated with the image-noise statistic 130.
In some embodiments of the present invention, the codec-noise statistic estimator 128 may calculate an estimate for the codec noise 132 according to FIG. 10. In these embodiments, a luminance image may be computed 140 from an input RGB (Red Green Blue), or other color space, image. The luminance may be denoted I(x, y), where x and y may denote the horizontal and vertical indices, respectively. A horizontal difference value at each point may be computed 142, thereby producing a plurality of horizontal difference values, and a vertical difference value at each point may be computed 144, thereby producing a plurality of vertical difference values, according to:

DH(x,y)=I(x+1,y)−I(x,y)

and

DV(x,y)=I(x,y+1)−I(x,y),

respectively, where DH(x, y) may denote a horizontal difference value at pixel location (x, y) and DV(x, y) may denote a vertical difference value at pixel location (x, y).

The standard deviation of the horizontal differences at a plurality of horizontal offsets may be computed 146, and the standard deviation of the vertical differences at a plurality of vertical offsets may be computed 148. In some embodiments of the present invention, a standard deviation value may be calculated for each offset within a coding block, thereby producing a plurality of horizontal standard deviation values and a plurality of vertical standard deviation values. The number of offsets may be determined by the structure of the image, or video, codec and any processing or scaling of the decoded data that may be performed prior to estimating the compression noise. In an exemplary embodiment comprising the use of eight offsets, the standard deviation values of the horizontal differences may be calculated 146 according to:



  
    
      
        
          Mean
          H
        
        ⁡
        
          [
          i
          ]
        
      
      =
      
        
          ∑
          
            x
            =
            0
          
          Height
        
        ⁢
        
          
            ∑
            
              y
              =
              0
            
            
              Width
              8
            
          
          ⁢
          
            
              
                D
                H
              
              ⁡
              
                (
                
                  x
                  ,
                  
                    
                      8
                      ·
                      y
                    
                    +
                    i
                  
                
                )
              
            
            
              Height
              ·
              
                (
                
                  Width
                  /
                  8
                
                )
              
            
          
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
  




  
    
      
        Mean
        ⁢
        
          
        
        ⁢
        
          
            2
            H
          
          ⁡
          
            [
            i
            ]
          
        
      
      =
      
        
          ∑
          
            x
            =
            0
          
          Height
        
        ⁢
        
          
            ∑
            
              y
              =
              0
            
            
              Width
              8
            
          
          ⁢
          
            
              
                
                  D
                  H
                
                ⁡
                
                  (
                  
                    x
                    ,
                    
                      
                        8
                        ·
                        y
                      
                      +
                      i
                    
                  
                  )
                
              
              2
            
            
              Height
              ·
              
                (
                
                  Width
                  /
                  8
                
                )
              
            
          
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
  




  
    
      
        
          STD
          H
        
        ⁡
        
          [
          i
          ]
        
      
      =
      
        sqrt
        ⁡
        
          (
          
            
              Mean
              ⁢
              
                
              
              ⁢
              
                
                  2
                  H
                
                ⁡
                
                  [
                  i
                  ]
                
              
            
            -
            
              
                
                  Mean
                  H
                
                ⁡
                
                  [
                  i
                  ]
                
              
              2
            
          
          )
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7.
      
    
  



FIG. 11 depicts an exemplary portion of an image 160 with pixel locations shown as squares. The pixel locations shown in cross hatch, for example 162, may be the locations used for computing the horizontal standard deviation associated with an offset of zero, STDV[0], and the pixel locations shown in white, for example 164, may be the locations used for computing the horizontal standard deviation associated with an offset of two, STDV[2]. The standard deviation values of the vertical differences may be similarly calculated 148 according to:



  
    
      
        
          Mean
          V
        
        ⁡
        
          [
          i
          ]
        
      
      =
      
        
          ∑
          
            x
            =
            0
          
          
            Height
            8
          
        
        ⁢
        
          
            ∑
            
              y
              =
              0
            
            Width
          
          ⁢
          
            
              
                D
                V
              
              ⁡
              
                (
                
                  
                    
                      8
                      ·
                      x
                    
                    +
                    i
                  
                  ,
                  y
                
                )
              
            
            
              
                (
                
                  Height
                  /
                  8
                
                )
              
              ·
              Width
            
          
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
  




  
    
      
        Mean
        ⁢
        
          
        
        ⁢
        
          
            2
            V
          
          ⁡
          
            [
            i
            ]
          
        
      
      =
      
        
          ∑
          
            x
            =
            0
          
          
            Height
            8
          
        
        ⁢
        
          
            ∑
            
              y
              =
              0
            
            Width
          
          ⁢
          
            
              
                
                  D
                  V
                
                ⁡
                
                  (
                  
                    
                      
                        8
                        ·
                        x
                      
                      +
                      i
                    
                    ,
                    y
                  
                  )
                
              
              2
            
            
              
                (
                
                  Height
                  /
                  8
                
                )
              
              ·
              Width
            
          
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
  




  
    
      
        
          STD
          V
        
        ⁡
        
          [
          i
          ]
        
      
      =
      
        sqrt
        ⁡
        
          (
          
            
              Mean
              ⁢
              
                
              
              ⁢
              
                
                  2
                  V
                
                ⁡
                
                  [
                  i
                  ]
                
              
            
            -
            
              
                
                  Mean
                  V
                
                ⁡
                
                  [
                  i
                  ]
                
              
              2
            
          
          )
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7.
      
    
  



Referring again to FIG. 10, statistics may be computed 150 for the horizontal standard deviation values, and statistics may be computed 152 for the vertical standard deviation values. The statistics computed 150 for the horizontal standard deviation values may be referred to as horizontal statistics and the statistics computed 152 for the vertical standard deviation values may be referred to as vertical statistics. In some embodiments of the present invention, the average of the horizontal standard deviation values and the maximum horizontal standard deviation value may be computed 150, and the average of the vertical standard deviation values and the maximum vertical standard deviation value may be computed 152. These values may be computed according to:



  
    
      STD
      
        H
        ⁢
        _
        ⁢
        MEAN
      
    
    =
    
      
        ∑
        
          i
          =
          0
        
        7
      
      ⁢
      
        
          
            STD
            H
          
          ⁡
          
            [
            i
            ]
          
        
        8
      
    
  




  
    
      
        STD
        
          H
          ⁢
          _
          ⁢
          MAX
        
      
      =
      
        max
        ⁡
        
          (
          
            
              STD
              H
            
            ⁡
            
              [
              i
              ]
            
          
          )
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
  




  
    
      STD
      
        V
        ⁢
        _
        ⁢
        MEAN
      
    
    =
    
      
        ∑
        
          i
          =
          0
        
        7
      
      ⁢
      
        
          
            STD
            V
          
          ⁡
          
            [
            i
            ]
          
        
        8
      
    
  




  
    
      
        STD
        
          V
          ⁢
          _
          ⁢
          MAX
        
      
      =
      
        max
        ⁡
        
          (
          
            
              STD
              V
            
            ⁡
            
              [
              i
              ]
            
          
          )
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
    ,
  



respectively. In alternative embodiments, other statistical values may be calculated, for example, the median, the skew, the kurtosis and other statistical measures.

The computed statistics associated with the horizontal standard deviation values may be combined 154 to form an estimate of the horizontal compression noise. The computed statistics associated with the vertical standard deviation values may be combined 156 to form an estimate of the vertical compression noise. In an exemplary embodiment of the present invention, an estimate of the horizontal compression noise may be calculated according to a weighted average given by:

NoiseCompression—H=4.64·STDH—MAX−4.26·STDH—MEAN+0.58,

and an estimate of the vertical compression noise may be calculated according to a weighted average given by:

NoiseCompression—V=4.64·STDV—MAX−4.26·STDV—MEAN+0.58,

where the values 4.64, 4.26 and 0.58 are exemplary weighting parameters. In alternative embodiments, other weighting values may be used.

The estimate for the horizontal compression noise and the estimate for the vertical compression noise may be combined 158 to form a single compression noise estimate. Any data fusion method known in the art may be used to combine the estimate for the horizontal compression noise and the estimate for the vertical compression noise. In some embodiments of the present invention, the compression-noise estimate may be determined according to:

NoiseCompression=max(NoiseCompression—H,NoiseCompression—V).

In alternative embodiments of the present invention, the codec-noise statistic estimator 128 may calculate an estimate for the codec noise 132 according to FIG. 12. In these embodiments, the input image may be a frame in a video sequence, and past frames may be considered in determining the noise estimate. In these embodiments, a luminance image may be computed 170 from an input RGB, or other color space, image, and the luminance may be denoted I(x, y), where x and y may denote the horizontal and vertical indices, respectively. A horizontal difference value at each point may be computed 172, thereby producing a plurality of horizontal difference values, and a vertical difference value at each point may be computed 174, thereby producing a plurality of vertical difference values, according to:

DH(x,y)=I(x+1,y)−I(x,y)

and

DV(x,y)=I(x,y+1)−I(x,y),

respectively, where DH(x, y) may denote a horizontal difference value at pixel location (x, y) and DV(x, y) may denote a vertical difference value at pixel location (x, y).

The standard deviation of the horizontal differences at a plurality of horizontal offsets may be computed 176, and the standard deviation of the vertical differences at a plurality of vertical offsets may be computed 178. In some embodiments of the present invention, a standard deviation value may be calculated for each offset within a coding block, thereby producing a plurality of horizontal standard deviation values and a plurality of vertical standard deviation values. The number of offsets may be determined by the structure of the video codec and any processing or scaling of the decoded data that may be performed prior to estimating the compression noise. In an exemplary embodiment comprising the use of eight offsets, the standard deviation values of the horizontal differences may be calculated 176 according to:



  
    
      
        
          Mean
          H
        
        ⁡
        
          [
          i
          ]
        
      
      =
      
        
          ∑
          
            x
            =
            0
          
          Height
        
        ⁢
        
          
            ∑
            
              y
              =
              0
            
            
              Width
              8
            
          
          ⁢
          
            
              
                D
                H
              
              ⁡
              
                (
                
                  x
                  ,
                  
                    
                      8
                      ·
                      y
                    
                    +
                    i
                  
                
                )
              
            
            
              Height
              ·
              
                (
                
                  Width
                  /
                  8
                
                )
              
            
          
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
  




  
    
      
        Mean
        ⁢
        
          
        
        ⁢
        
          
            2
            H
          
          ⁡
          
            [
            i
            ]
          
        
      
      =
      
        
          ∑
          
            x
            =
            0
          
          Height
        
        ⁢
        
          
            ∑
            
              y
              =
              0
            
            
              Width
              8
            
          
          ⁢
          
            
              
                
                  D
                  H
                
                ⁡
                
                  (
                  
                    x
                    ,
                    
                      
                        8
                        ·
                        y
                      
                      +
                      i
                    
                  
                  )
                
              
              2
            
            
              Height
              ·
              
                (
                
                  Width
                  /
                  8
                
                )
              
            
          
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
  




  
    
      
        
          STD
          H
        
        ⁡
        
          [
          i
          ]
        
      
      =
      
        sqrt
        ⁡
        
          (
          
            
              Mean
              ⁢
              
                
              
              ⁢
              
                
                  2
                  H
                
                ⁡
                
                  [
                  i
                  ]
                
              
            
            -
            
              
                
                  Mean
                  H
                
                ⁡
                
                  [
                  i
                  ]
                
              
              2
            
          
          )
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7.
      
    
  



The standard deviation values of the vertical differences may be similarly calculated 178 according to:




  
    
      
        
          Mean
          V
        
        ⁡
        
          [
          i
          ]
        
      
      =
      
        
          ∑
          
            x
            =
            0
          
          
            Height
            8
          
        
        ⁢
        
          
            ∑
            
              y
              =
              0
            
            Width
          
          ⁢
          
            
              
                D
                V
              
              ⁡
              
                (
                
                  
                    
                      8
                      ·
                      x
                    
                    +
                    i
                  
                  ,
                  y
                
                )
              
            
            
              
                (
                
                  Height
                  /
                  8
                
                )
              
              ·
              Width
            
          
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
  




  
    
      
        Mean
        ⁢
        
          
        
        ⁢
        
          
            2
            V
          
          ⁡
          
            [
            i
            ]
          
        
      
      =
      
        
          ∑
          
            x
            =
            0
          
          
            Height
            8
          
        
        ⁢
        
          
            ∑
            
              y
              =
              0
            
            Width
          
          ⁢
          
            
              
                
                  D
                  V
                
                ⁡
                
                  (
                  
                    
                      
                        8
                        ·
                        x
                      
                      +
                      i
                    
                    ,
                    y
                  
                  )
                
              
              2
            
            
              
                (
                
                  Height
                  /
                  8
                
                )
              
              ·
              Width
            
          
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
  




  
    
      
        
          STD
          V
        
        ⁡
        
          [
          i
          ]
        
      
      =
      
        sqrt
        ⁡
        
          (
          
            
              Mean
              ⁢
              
                
              
              ⁢
              
                
                  2
                  V
                
                ⁡
                
                  [
                  i
                  ]
                
              
            
            -
            
              
                
                  Mean
                  V
                
                ⁡
                
                  [
                  i
                  ]
                
              
              2
            
          
          )
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7.
      
    
  



Statistics may be computed 180 for the horizontal standard deviation values, and statistics may be computed 182 for the vertical standard deviation values. In some embodiments of the present invention, the average of the horizontal standard deviation values and the maximum horizontal standard deviation value may be computed 850, and the average of the vertical standard deviation values and the maximum vertical standard deviation value may be computed 182. These values may be computed according to:



  
    
      STD
      
        H
        ⁢
        _
        ⁢
        MEAN
      
    
    =
    
      
        ∑
        
          i
          =
          0
        
        7
      
      ⁢
      
        
          
            STD
            H
          
          ⁡
          
            [
            i
            ]
          
        
        8
      
    
  




  
    
      
        STD
        
          H
          ⁢
          _
          ⁢
          MAX
        
      
      =
      
        max
        ⁡
        
          (
          
            
              STD
              H
            
            ⁡
            
              [
              i
              ]
            
          
          )
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
  




  
    
      STD
      
        V
        ⁢
        _
        ⁢
        MEAN
      
    
    =
    
      
        ∑
        
          i
          =
          0
        
        7
      
      ⁢
      
        
          
            STD
            V
          
          ⁡
          
            [
            i
            ]
          
        
        8
      
    
  




  
    
      
        STD
        
          V
          ⁢
          _
          ⁢
          MAX
        
      
      =
      
        max
        ⁡
        
          (
          
            
              STD
              V
            
            ⁡
            
              [
              i
              ]
            
          
          )
        
      
    
    ,
    
      i
      =
      
        0
        ⁢
        
          
        
        ⁢
        ...
        ⁢
        
          
        
        ⁢
        7
      
    
    ,
  



respectively. In alternative embodiments, other statistical values may be calculated, for example, the median, the skew, the kurtosis and other statistical measures.

The computed statistics associated with the horizontal standard deviation values may be combined 184 to form an estimate of the horizontal compression noise. The computed statistics associated with the vertical standard deviation values may be combined 186 to form an estimate of the vertical compression noise. In an exemplary embodiment of the present invention, an estimate of the horizontal compression noise may be calculated according to a weighted average given by:

NoiseCompression—H=4.64·STDH—MAX−4.26·STDH—MEAN+0.58,

and an estimate of the vertical compression noise may be calculated according to a weighted average given by:

NoiseCompression—V=4.64·STDV—MAX−4.26·STDV—MEAN+0.58,

where the values 4.64, 4.26 and 0.58 are exemplary weighting parameters. In alternative embodiments, other weighting values may be used.

The location of block boundaries in the current image frame may be estimated 188, 190 using the statistics previously computed. In an exemplary embodiment, a horizontal-block boundary location may be estimated 188 using the maximum horizontal standard deviation value according to:

BlockH[frame_num]=i, where STDH—MAX==STDH[i], 

where frame_num may denote a time index associated with the current frame. A vertical block-boundary location may be estimated 190 by the maximum vertical standard deviation value according to:

BlockV[frame_num]=i, where STDV—MAX==STDV[i], 

where frame_num may denote a time index associated with the current frame.

The number of unique horizontal block-boundary locations and the number of unique vertical block-boundary locations in a temporal block may be determined 192, 194 by counting the number of unique values for BlockH[i] and BlockV[i], respectively, where i is an index with values from frame_num to frame_num-N and N is a constant. If the number of unique values for the horizontal direction is above a threshold, then the estimate for the horizontal compression noise may be set equal to zero, and if the number of unique values for the vertical direction is above a threshold, then the estimate for the vertical compression noise may be set equal to zero.
The estimate for the horizontal compression noise and the estimate for the vertical compression noise may be combined 196 to form a single compression noise estimate. Any data fusion method known in the art may be used to combine the estimate for the horizontal compression noise and the estimate for the vertical compression noise. In some embodiments of the present invention, the compression-noise estimate may be determined according to:

NoiseCompress=max(NoiseCompression—H,NoiseCompression—V).

In some embodiments of the present invention, the parameters of an edge-preserving filter may be adjusted based on a sharpening value used in a sharpening filter. In some exemplary embodiments, the edge-preserving filter may comprise a bi-lateral filter. In exemplary embodiments of the present invention shown in FIG. 13 and FIG. 14, an edge-preserving filter 204 may be used to separate an input image 200 into two channels: a first channel 206 which may comprise image content to be sharpened by un-sharp masking (USM) 212; and a second channel 210 which may comprise image content to attenuate. The level of sharpening associated with the un-sharp masking 212 may be controlled by the value of the USM alpha 213 set for the USM filter 212. In these embodiments, a filter parameter value 203 of the edge-preserving filter 204 may control the assignment of input-image content to the first channel 206 and the second channel 210. In some embodiments of the present invention wherein the edge-preserving filter comprises a bi-lateral filter, the filter parameter value may be the range sigma parameter value associated with the bi-lateral filter.
In some embodiments, in order to maximize the amount of sharpening, the filter parameter 203 may be set based on the noise statistics associated with the input image 200, and additionally, in order to mitigate noise amplification due to sharpening, the filter parameter 203 may be set based on the sharpening value 213 used in the sharpening filter 212.
An input image 200 and a filter parameter value 203 associated with the input image may be made available to an edge-preserving filter 204. The edge-preserving filter 204 may be used to divide the input image 200 into two channels 206, 210. The first channel 206 may correspond to significant features in the input image 200, and the second channel 210, formed by removing 208 the first channel 206 from the input image 200, may contain the residual difference between the first channel 206 and the input image 200. The residual difference 210 may comprise noise and texture data. The first channel 206 may be sharpened using un-sharp masking 212, and the second channel 210 may be attenuated 214, wherein the amount of attenuation may be controlled by a noise control parameter 215. The sharpened channel 216 and the attenuated channel 218 may be combined 220 to form an enhanced image 222. The filter parameter value 203 of the edge-preserving filter 204 may control the assignment of input-image content to the first channel 206 and the second channel 210. In some embodiments of the present invention (not shown), the enhanced image 222 may be up-sampled to a resolution greater than that of the input image 200.
Referring to FIG. 14, the input image 200 may be received by a noise estimator 230 which may generate a noise estimate 231 based on an estimate of the image noise and the codec noise present in the input image 200. The noise estimate 231 may be determined according to the previously described embodiments of the present invention. The noise estimate 231 may be adjusted by a USM sharpening compensator 232 to account for the sharpening value, and the resulting sharpening-adjusted value 233 may be temporally filtered 234 to reduce noise in the control process. The output of the temporal filter 234 may be the filter parameter value 203 provided to the edge-preserving filter 204.
In some embodiments of the present invention described in relation to FIG. 14 and FIG. 15, the noise estimator 230 may comprise an image-noise statistic estimator 250 and a codec-noise statistic estimator 252 for estimating image noise and codec noise associated with an input image 200, respectively. The noise estimator 230 may comprise a combiner 254 for combining the image noise estimated by the image-noise statistic estimator 250 and the codec noise codec-noise statistic estimator 252. The noise estimate 231 may be provided to the USM sharpening compensator 232. In some embodiments, the noise estimate 231 may be converted to a standard deviation value in the noise estimator 230. In alternative embodiments, the noise estimate may be converted to a standard deviation value in the USM sharpening compensator 232.
The USM sharpening compensator 232 may adjust the noise estimate 231 by a value associated with the level of sharpening 213. In some embodiments of the present invention, a sharpening compensator controller 256 may use the USM alpha value 213 as an index into a look-up table 258 to determine a multiplier 257 which may be used to multiplicatively adjust 260 the standard-deviation-value noise estimate 231. FIG. 16 depicts two plots 280, 282 each associated with an exemplary look-up table derived to maintain a constant noise level with increasing sharpening factor.
In addition to adjusting the noise estimate 231 by the value 257 associated with the level of sharpening 213, the noise estimate 231 may be multiplicatively 260 adjusted by a value associated with control parameter, which may be denoted M, that controls how quickly the estimated parameter values may change as a function of time. In some embodiments the multiplicative adjustment 260 may be 1−M. In some embodiments of the present invention, a smaller value for M may correspond to a faster change, while a larger value for M may correspond to a slower change. In some embodiments of the present invention, M may be set to 0.5. A filter parameter value associated with a previous frame may be retrieved from storage 262 and multiplied 264 by the control parameter M. The results of the two multiplications 260, 264 may be added 266 to form the edge-preserving filter parameter 203 for the next frame.
In alternative embodiments described in relation to FIG. 17, an edge-preserving filter parameter may be determined for filtering of a still image, or a video frame. An input image 300 may be received by a noise estimator 302 which may generate a noise estimate 303 based on an estimate of the image noise and the codec noise present in the input image 300. The noise estimate 303 may be determined according to the previously described embodiments of the present invention. The noise estimate 303 may be adjusted by a USM sharpening compensator 304 to account for the sharpening value 305, and the resulting sharpening-adjusted value 306 may be the filter parameter value 306 provided to the edge-preserving filter 204.
In some embodiments of the present invention described in relation to FIG. 17 and FIG. 18, the noise estimator 302 may comprise an image-noise statistic estimator 310 and a codec-noise statistic estimator 312 for estimating image noise and codec noise associated with an input image 300, respectively. The noise estimator 302 may comprise a combiner 314 for combining the image noise estimated by the image-noise statistic estimator 310 and the codec noise codec-noise statistic estimator 312. The noise estimate 303 may be provided to the USM sharpening compensator 304. In some embodiments, the noise estimate 303 may be converted to a standard deviation value in the noise estimator 302. In alternative embodiments, the noise estimate 303 may be converted to a standard deviation value in the USM sharpening compensator 304.
The USM sharpening compensator 304 may adjust the noise estimate 303 by a value associated with the level of sharpening 305. In some embodiments of the present invention, a sharpening compensator controller 316 may use the USM alpha value 305 as an index into a look-up table 318 to determine a multiplier 319 which may be used to multiplicatively adjust 320 the standard-deviation-value noise estimate 303 to produce a sharpening-adjusted edge-preserving filter parameter value 306. FIG. 16 depicts two plots 280, 282 each associated with an exemplary look-up table derived to maintain a constant noise level with increasing sharpening factor.
The effectiveness of embodiments of the present invention may be illustrated in FIGS. 19-25. FIG. 19 depicts an exemplary, original, synthetic noise target 350. FIG. 20 shows the result 352 of bi-lateral filtering, using a sigma range parameter value equal to the standard deviation of the noise, of the original image 250 shown in FIG. 19. An enhanced image 354 is shown in FIG. 21. This enhanced image 354 exhibits significant amplification of the image noise. FIG. 22 depicts an enhanced image 356 enhanced according to embodiments of the present invention in which the bi-lateral filter sigma range parameter is set according to the noise statistics of the image in addition to accounting for the sharpening strength of the sharpening filter. The noise in FIG. 22 is similar to the noise in FIG. 20, but the degree of enhancement is different.
A second exemplary, original image 360 is depicted in FIG. 23. A superimposed, white circle 361 is shown on the image 360 in a region wherein the noise in the image 360 is visible. FIG. 24 depicts an image 362 which results from bi-lateral filtering according to embodiments of the present invention with a range sigma parameter that has been set accounting for the noise, but not the sharpening process. Inspection of the region inside the overlaid white circle 363 shows that the noise is amplified relative to the input image 360. FIG. 25 depicts an image 364 which results from bi-lateral filtering according to embodiments of the present invention with a range sigma parameter that has been set accounting for both noise and the sharpening process. Inspection of the region inside the overlaid white circle 365 shows that the noise is not amplified relative to the original image 360, but the desired increase in sharpening of image content is visible.
In some embodiments of the present invention, a noise-control parameter that controls the amount of attenuation of the noise and texture data may be automatically selected based on the amount of compression detected in the input image.
In exemplary embodiments of the present invention shown in FIG. 26, an edge-preserving filter 400 may be used to separate an input image 402 into two channels: a first channel 404 which may comprise image content to be sharpened by un-sharp masking (USM) 406; and a second channel 408 which may comprise image content to attenuate 410.
Then input image 402 may be made available to the edge-preserving filter 400. The edge-preserving filter 400 may be used to divide the input image 402 into two channels 404, 408. The first channel 404 may correspond to significant features in the input image 402, and the second channel 408, formed by removing 410 the first channel 404 from the input image 402, may contain the residual difference between the first channel 404 and the input image 402. The residual difference 408 may comprise noise and texture data. The first channel 404 may be sharpened using un-sharp masking 406, and the second channel 408 may be attenuated 416, wherein the amount of attenuation may be controlled by a noise-control parameter 412. The sharpened channel 414 and the attenuated channel 416 may be combined 418 to form an enhanced image 420.
The noise-control parameter 412 may be selected by a noise-control-parameter selector 422 based on the amount of compression detected in the input image 402. The input image 402 may be made available to a noise estimator 424 which may provide to the noise-control-parameter selector 422 an estimate 426 of the image noise, also considered the additive noise, and the compression noise in the input image 402.
In some embodiments of the present invention described in relation to FIG. 27, the noise estimator 424 may comprise an image-noise statistic estimator 450 and a codec-noise statistic estimator 452 for estimating image noise 451, which may be denoted NA, and codec noise 453, which may be denoted NC=Noisecompression, respectively. In some embodiments, the noise estimator 424 may comprise a combiner 454 for combining the image noise estimate 451 and the codec-noise estimate 453. In alternative embodiments (not shown), the noise estimator may not comprise a combiner.
The image-noise estimate 451 and the compression-noise estimate 453 may be provided to the noise-control-parameter selector 422. The noise-control-parameter selector 422 may comprise a noise-control-parameter-value calculator 456, which may receive the image-noise estimate 451 and the compression-noise estimate 453. The noise-control-parameter-value calculator 456 may determine a ratio between the image-noise estimate 451 and the total noise estimate and calculate a noise-control-parameter value 457, which may be denoted β, according to




  
    β
    =
    
      
        
          N
          A
        
        
          
            N
            A
          
          +
          
            N
            C
          
        
      
      ⁢
      
        E
        1
      
    
  
  ,




where NA and NC denote the image-noise estimate 451 and the compression-noise estimate 453, respectively, and E1 denotes the beta-preference parameter which controls the amount of attenuation when there is no compression. In alternative embodiments, the noise-control-parameter value 457 may be determined using a look-up-table, mapping or other function, that takes the values of NA and NC as input and provides a β value as output. An exemplary mapping 500 is illustrated in FIG. 28.

In some embodiments of the present invention, the calculated noise-control-parameter value 457 may be temporally smoothed by adjusting the calculated noise-control-parameter value 457 according to a control parameter, which may be denoted Mnc, that controls how quickly the noise-control-parameter values may change as a function of time. In some embodiments of the present invention, a smaller value for Mnc may correspond to a faster change, while a larger value for Mnc may correspond to a slower change. In some embodiments of the present invention, Mnc may be set to 0.5. A noise-control-parameter value associated with a previous frame 458 may be retrieved from storage 460 and multiplied 462 by the control parameter Mnc. The noise-control-parameter value 457 calculated by the noise-control-parameter-value calculator 456 for the current frame may be multiplied 464 by 1−Mnc. The results of the two multiplications may be added 468 to form the noise-control-parameter value for the next frame 412.
In some embodiments of the present invention described in relation to FIG. 29, the parameter values of an edge-preserving filter may be adjusted based on a sharpening value used in a sharpening filter, and a noise-control parameter value controlling an attenuator may be selected based on the amount of compression detected in an input image. In some exemplary embodiments, the edge-preserving filter may comprise a bi-lateral filter.
In exemplary embodiments of the present invention shown in FIG. 29, an edge-preserving filter 550 may be used to separate an input image 552 into two channels: a first channel 554 which may comprise image content to be sharpened by un-sharp masking (USM) 556; and a second channel 558 which may comprise image content to attenuate 560. The level of sharpening associated with the un-sharp masking 556 may be controlled by the value of the USM alpha 562 set for the USM filter 556. In these embodiments, a filter parameter value 564 of the edge-preserving filter 550 may control the assignment of input-image content to the first channel 554 and the second channel 558. In some embodiments of the present invention wherein the edge-preserving filter comprises a bi-lateral filter, the filter parameter value may be the range sigma parameter value associated with the bi-lateral filter.
In some embodiments, in order to maximize the amount of sharpening, the filter parameter 564 may be set based on the noise statistics associated with the input image 552, and additionally, in order to mitigate noise amplification due to sharpening, the filter parameter 564 may be set based on the sharpening value 562 used in the sharpening filter 556.
The input image 552 and the filter parameter value 564 associated with the input image 552 may be made available to the edge-preserving filter 550. The edge-preserving filter 550 may be used to divide the input image 552 into two channels 554, 558. The first channel 554 may correspond to significant features in the input image 552, and the second channel 558, formed by removing 566 the first channel 554 from the input image 552, may contain the residual difference between the first channel 554 and the input image 552. The residual difference 558 may comprise noise and texture data. The first channel 554 may be sharpened using un-sharp masking 556, and the second channel 558 may be attenuated 560, wherein the amount of attenuation may be controlled by a noise control parameter 568. The sharpened channel 570 and the attenuated channel 572 may be combined 574 to form an enhanced image 576. The filter parameter value 564 of the edge-preserving filter 550 may control the assignment of input-image content to the first channel 554 and the second channel 558. In some embodiments of the present invention (not shown), the enhanced image 576 may be up-sampled to a resolution greater than that of the input image 552.
The input image 552 may be received by a noise estimator 578 which may generate a noise estimate 580 based on an estimate of the image noise and the codec noise present in the input image 552. The noise estimate 580 may be determined according to the previously described embodiments of the present invention and may comprise an image-noise estimate, which may be denoted NA, and a compression-noise estimate, which may be denoted NC. The noise estimate 580 may be adjusted by a USM sharpening compensator 582 to account for the sharpening value 562, and the resulting sharpening-adjusted value 584 may be temporally filtered 586 to reduce noise in the control process. The output of the temporal filter 586 may be the filter parameter value 564 provided to the edge-preserving filter 550.
The noise-control parameter 568 may be selected by a noise-control-parameter selector 588 based on the amount of compression detected in the input image 552. The noise estimator 578 may make the noise estimate 580 available to the noise-control-parameter selector 588.
Some embodiments of the present invention described in relation to FIG. 29 may be further understood in relation to FIG. 30.
The noise estimator 578 may comprise an image-noise statistic estimator 600 and a codec-noise statistic estimator 602 for estimating image noise and codec noise associated with an input image 552, respectively. The noise estimator 578 may comprise a combiner 604 for combining the image noise estimate 606 estimated by the image-noise statistic estimator 600 and the codec noise estimate 608 estimated by the codec-noise statistic estimator 602. The combined noise estimate 610 may be provided to the USM sharpening compensator 582. In some embodiments, the combined noise estimate 610 may be converted to a standard deviation value in the noise estimator 578. In alternative embodiments, the combined noise estimate 610 may be converted to a standard deviation value in the USM sharpening compensator 582.
The USM sharpening compensator 582 may adjust the combined noise estimate 610 by a value associated with the level of sharpening 562. In some embodiments of the present invention, a sharpening compensator controller 612 may use the USM alpha value 562 as an index into a look-up table 614 to determine a multiplier 616 which may be used to multiplicatively adjust 618 the standard-deviation-value combined noise estimate 610. FIG. 16 depicts two plots 280, 282 each associated with an exemplary look-up table derived to maintain a constant noise level with increasing sharpening factor.
In addition to adjusting the combined noise estimate 610 by the value 616 associated with the level of sharpening 562, the combined noise estimate 610 may be multiplicatively 618 adjusted by a value associated with control parameter, which may be denoted M, that controls how quickly the estimated parameter values may change as a function of time. In some embodiments the multiplicative adjustment 618 may be 1−M. In some embodiments of the present invention, a smaller value for M may correspond to a faster change, while a larger value for M may correspond to a slower change. In some embodiments of the present invention, M may be set to 0.5. A filter parameter value 620 associated with a previous frame may be retrieved from storage 622 and multiplied 624 by the control parameter M. The results of the two multiplications may be added 626 to form the edge-preserving filter parameter 564 for the next frame.
The image-noise estimate 606 and the compression-noise estimate 608 may be provided to the noise-control-parameter selector 588. The noise-control-parameter selector 588 may comprise a noise-control-parameter-value calculator 630, which may receive the image-noise estimate 606 and the compression-noise estimate 608. The noise-control-parameter-value calculator 630 may determine a ratio between the image-noise estimate 606 and the total noise estimate and calculate a noise-control-parameter value 632, which may be denoted β, according to




  
    β
    =
    
      
        
          N
          A
        
        
          
            N
            A
          
          +
          
            N
            C
          
        
      
      ⁢
      
        E
        1
      
    
  
  ,




where NA and NC denote the image-noise estimate 606 and the compression-noise estimate 608, respectively, and E1 denotes the beta-preference parameter which controls the amount of attenuation when there is no compression. In alternative embodiments, the noise-control-parameter value 632 may be determined using a look-up-table, mapping or other function, that takes the values of NA and NC as input and provides a β value as output. An exemplary mapping 500 is illustrated in FIG. 28.

In some embodiments of the present invention, the calculated noise-control-parameter value 632 may be temporally smoothed by adjusting the calculated noise-control-parameter value 632 according to a control parameter, which may be denoted Mnc, that controls how quickly the noise-control-parameter values may change as a function of time. In some embodiments of the present invention, a smaller value for Mnc may correspond to a faster change, while a larger value for Mnc may correspond to a slower change. In some embodiments of the present invention, Mnc may be set to 0.5. A noise-control-parameter value associated with a previous frame 634 may be retrieved from storage 636 and multiplied 638 by the control parameter Mnc. The noise-control-parameter value 632 calculated by the noise-control-parameter-value calculator 630 for the current frame may be multiplied 640 by 1−Mnn. The results of the two multiplications may be added 642 to form the noise-control-parameter value for the next frame 568.
Some embodiments of the present invention described herein comprise an edge-preserving filter. In some embodiments, the edge-preserving filter may comprise a bi-lateral filter. A person of ordinary skill in the art will recognize the existence of many edge-preserving filters and many forms of bi-lateral filters.
Although the charts and diagrams in the figures may show a specific order of execution, it is understood that the order of execution may differ from that which is depicted. For example, the order of execution of the blocks may be changed relative to the shown order. Also, as a further example, two or more blocks shown in succession in the figure may be executed concurrently, or with partial concurrence. It is understood by those with ordinary skill in the art that software, hardware and/or firmware may be created by one of ordinary skill in the art to carry out the various logical functions described herein.
Some embodiments of the present invention may comprise a computer program product comprising a computer-readable storage medium having instructions stored thereon/in which may be used to program a computing system to perform any of the features and methods described herein. Exemplary computer-readable storage media may include, but are not limited to, flash memory devices, disk storage media, for example, floppy disks, optical disks, magneto-optical disks, Digital Versatile Discs (DVDs), Compact Discs (CDs), micro-drives and other disk storage media, Read-Only Memory (ROMs), Programmable Read-Only Memory (PROMs), Erasable Programmable Read-Only Memory (EPROMS), Electrically Erasable Programmable Read-Only Memory (EEPROMs), Random-Access Memory (RAMS), Video Random-Access Memory (VRAMs), Dynamic Random-Access Memory (DRAMs) and any type of media or device suitable for storing instructions and/or data.
The terms and expressions which have been employed in the foregoing specification are used therein as terms of description and not of limitation, and there is no intention in the use of such terms and expressions of excluding equivalence of the features shown and described or portions thereof, it being recognized that the scope of the invention is defined and limited only by the claims which follow.

## Claims

1. A system for enhancing an input image, said system comprising:
an edge-preserving filter for generating, from an input image, a first image comprising image content, from said input image, to sharpen and a second image comprising image content, from said input image, to attenuate; and
a noise-control-parameter selector for selecting a value for a noise-control parameter that controls the amount of attenuation applied to said second image, wherein said noise-control-parameter selector selects said value based on a first noise estimate of an additive noise in said input image and a second noise estimate of a compression noise in said input image, and wherein said value is a scalar value.

2. The system as described in claim 1, wherein said edge-preserving filter comprises a bi-lateral filter.

3. The system as described in claim 1, wherein said value is proportional, by a proportionality constant, to the ratio of the first noise estimate and the sum of the first noise estimate and the second noise estimate.

4. The system as described in claim 3, wherein said value is determined by accessing a look-up-table using said first noise estimate and said second noise estimate as indices into said look-up-table.

5. The system as described in claim 3, wherein said proportionality constant is a beta preference parameter associated with no compression in said input image.

6. The system as described in claim 1 further comprising a noise estimator for estimating said first noise estimate and said second noise estimate.

7. The system as described in claim 6, wherein said noise estimator comprises:
an image-noise statistic estimator for estimating at least one of grain noise, sensor noise and thermal noise in said input image, thereby producing an image-noise statistic for generating said first noise estimate; and
a codec-noise statistic estimator for estimating compression noise in said input image, thereby producing a codec-noise statistic for generating said second noise estimate.

8. The system as described in claim 7, wherein said codec-noise statistic estimator estimates said codec noise statistic from said input image.

9. The system as described in claim 7, wherein said codec-noise statistic estimator estimates said codec-noise statistic from coding parameters associated with said input image.

10. The system as described in claim 9, wherein said codec-noise statistic estimator comprises a plurality of look-up tables, wherein each look-up table in said plurality of look-up tables is associated with a codec and at least one codec coding parameter.

11. The system as described in claim 1 further comprising a decoder for decoding a bitstream associated with said input image.

12. The system as described in claim 1 further comprising:
a sharpener for sharpening said first image;
an attenuator for attenuating said second image; and
an image combiner for combining said sharpened first image and said attenuated second image.

13. The system as described in claim 12 further comprising an up-sampler for up-sampling said combined image.

14. The system as described in claim 12, wherein said sharpener comprises an un-sharp masking filter.

15. The system as described in claim 12, wherein said noise-control parameter controls said attenuator.

16. The system as described in claim 1, wherein said noise-control-parameter selector comprises a noise-control-parameter-value calculator.

17. The system as described in claim 16, wherein said noise-control-parameter-value calculator calculates a first value for said noise-control parameter based on a first noise estimate of an additive noise in said input image and a second noise estimate of a compression noise in said input image.

18. The system as described in claim 17, wherein said first value for said noise-control parameter calculated by said noise-control-parameter-value calculator is combined with a previously stored value for said noise-control parameter.

19. A method for calculating a value for a noise-control parameter for controlling the amount of attenuation of a first image, said method comprising:
receiving a first noise estimate of the additive noise present in a first image;
receiving a second noise estimate of the compression noise present in said first image;
calculating a first value for said noise-control parameter based on said first noise estimate and said second noise estimate, wherein said first value is a scalar value; and
combining said first value for said noise-control parameter with a previously stored value for said noise-control parameter.