Abstract:
A method ( 50 ) and apparatus ( 40 ) for reducing noise in a video signal calculates ( 53 ) a discrimination value of a filter based on a magnitude difference between a current frame and a prior frame of a sequence of frames and modifies ( 54 ) the current frame using the modified discrimination value. The output of the filter is calculated as follows: 
         g ( i, j, k )= f ( i, j, k )±β 
     wherein:
   f(i, j, k) is the current frame input to the filter;   g(i, j, k) is the current frame output from the filter; and   β is the discrimination value of the filter, which is calculated as follows:   
 
     
       
         
           
             β 
             = 
             
               { 
               
                 
                   
                     
                       
                         
                           κδ 
                           ′ 
                         
                         + 
                         χ 
                       
                       , 
                     
                   
                   
                     
                       
                         
                           δ 
                           ′ 
                         
                         = 
                         
                           max 
                            
                           
                             { 
                             
                               δ 
                               , 
                               
                                 δ 
                                 0 
                               
                             
                             } 
                           
                         
                       
                       , 
                     
                   
                   
                     
                       δ 
                       &lt; 
                       
                         δ 
                         1 
                       
                     
                   
                 
                 
                   
                     
                       0 
                       , 
                     
                   
                   
                     
                         
                     
                   
                   
                     
                       δ 
                       ≥ 
                       
                         δ 
                         1 
                       
                     
                   
                 
               
             
           
         
       
       
         
           
             in which: 
             κ and X are constants less than one; 
             δ is a magnitude difference between the current frame and the prior frame; and 
             δ 0  and δ 1  are constants that satisfy δ 0 ≦δ 1 .

Description:
FIELD OF THE INVENTION 
       [0001]    The present invention relates generally to methods and apparatuses for processing digital information, and more particularly to a method and apparatus for processing digital information prior to transmission over a network. 
       BACKGROUND 
       [0002]    Motion compensated temporal filters (MCTF) are used in pixel processing of video images. MCTF is an adaptive temporal filter along estimated motion trajectories. The temporal correlation is exploited through temporal filtering. MCTF is most effective in processing sequences of images in which signal is highly temporal correlated along the motion trajectories but the noise, by nature, is totally uncorrelated. In other words, the image does not change from frame to frame but the noise does. The noise probability distribution is also exploited. MCTF is primarily used to reduce the zero mean stationary white Gaussian noise. Such noise has a bell shaped distribution, and it is of primarily small magnitude. Large magnitude noise rarely happens. Based on these two observations, the magnitude differences between the current frame and the previous filtered frame are first calculated along the estimated motion trajectories. Small differences are considered primarily of noise. Large differences are considered primarily of signal content. 
         [0003]    However, existing MCTFs do not remove all of the noise inherent in a video signal, thereby wasting valuable bandwidth resources. 
         [0004]    The present invention is therefore directed to the problem of developing a method and apparatus for performing filtering of noise in a video encoding process that removes as much of the noise as possible without disrupting the video image. 
       SUMMARY OF THE INVENTION 
       [0005]    The present invention solves these and other problems by providing a method and apparatus for filtering noise in a video encoding process, in which an adaptation scheme based on differences between frames is employed, thereby allowing the signal content to pass through while still effectively suppressing the random noise. 
         [0006]    According to one aspect of the present invention, a method for filtering noise in a video encoding process determines a discrimination value of a filter based on a magnitude difference between a current frame and a prior frame of a sequence of frames and modifies the current frame using the determined discrimination value. 
         [0007]    According to another aspect of the present invention, the current frame may be modified as follows: 
         [0000]        g ( i, j, k )= f ( i, j, k )±β 
         [0000]    wherein f(i, j, k) is the current frame input to the filter; g(i, j, k) is the current frame output from the filter; and ft is the discrimination value of the filter. 
         [0008]    According to still another aspect of the present invention, the discrimination value may be determined as follows: 
         [0000]    
       
         
           
             β 
             = 
             
               { 
               
                 
                   
                     
                       
                         δ 
                         - 
                         
                           
                             ( 
                             
                               1 
                               / 
                               Δ 
                             
                             ) 
                           
                            
                           
                             δ 
                             2 
                           
                         
                       
                       , 
                     
                   
                   
                     
                       δ 
                       &lt; 
                       Δ 
                     
                   
                 
                 
                   
                     
                       0 
                       , 
                     
                   
                   
                     
                       δ 
                       ≥ 
                       Δ 
                     
                   
                 
               
             
           
         
       
     
         [0000]    in which δ is the magnitude difference between the current frame and the prior filtered frame along the motion trajectories, and Δ is a predetermined filter strength. 
         [0009]    According to yet another aspect of the present invention, the discrimination value may be determined as follows: 
         [0000]    
       
         
           
             β 
             = 
             
               { 
               
                 
                   
                     
                       
                         
                           κδ 
                           ′ 
                         
                         + 
                         χ 
                       
                       , 
                     
                   
                   
                     
                       
                         
                           δ 
                           ′ 
                         
                         = 
                         
                           max 
                            
                           
                             { 
                             
                               δ 
                               , 
                               
                                 δ 
                                 0 
                               
                             
                             } 
                           
                         
                       
                       , 
                     
                   
                   
                     
                       δ 
                       &lt; 
                       
                         δ 
                         1 
                       
                     
                   
                 
                 
                   
                     
                       0 
                       , 
                     
                   
                   
                     
                         
                     
                   
                   
                     
                       δ 
                       ≥ 
                       
                         δ 
                         1 
                       
                     
                   
                 
               
             
           
         
       
     
         [0000]    in which κ and X are constants less than one; δ is a magnitude difference between the current frame and the prior filtered frame along the motion trajectories; δ 0  and δ 1  are constants that satisfy δ 0 ≦δ 1 ; and κ, X δ 0 , and δ 1  can be modified to adjust the filter strength. 
         [0010]    According to still another aspect of the present invention, an apparatus for encoding a video signal includes two encoders and a noise reduction preprocessor. The first encoder receives the video signal. The second encoder receives a plurality of parameters from the first encoder to use in encoding the video signal. The noise reduction preprocessor receives pixel data from the first encoder, filters the pixel data by determining a discrimination value of a filter based on a magnitude difference between a current frame and a prior frame of a sequence of frames, and by modifying the current frame using the modified discrimination value, and then provides preprocessed video to the second encoder. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS  
         [0011]      FIG. 1  depicts a plot of a linear adaptive filter used in a Motion Compensation Temporal Filter (MCTF) and the optimal filter according to one aspect of the present invention. 
           [0012]      FIG. 2  depicts a mapping of the linear adaptive filter in  FIG. 1  to a different domain and a plot of an exemplary embodiment of a MCTF according to another aspect of the present invention. 
           [0013]      FIG. 3  depicts the exemplary embodiment of  FIG. 2  in the original domain along with the linear adaptive filter and the optimal filter. 
           [0014]      FIG. 4  depicts a block diagram of an exemplary embodiment of a pixel processing system, which employs the MCTF of the present invention according to still another aspect of the present invention. 
           [0015]      FIG. 5  depicts a flow chart of an exemplary embodiment of a method for filtering images in a video signal according to yet another aspect of the present invention. 
           [0016]      FIG. 6  depicts a block diagram of an exemplary embodiment of a MCTF for filtering images in a video signal according to yet another aspect of the present invention. 
           [0017]      FIG. 7  depicts a plot of a processing order used in a bi-directional MCTF according to yet another aspect of the present invention. 
       
    
    
     DETAILED DESCRIPTION  
       [0018]    It is worthy to note that any reference herein to “one embodiment” or “an embodiment” means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the invention. The appearances of the phrase “in one embodiment” in various places in the specification are not necessarily all referring to the same embodiment. 
         [0019]    MCTF is used in a dual-pass video encoder. MCTF is employed between the first pass and second pass encoders. The first pass encoder provides motion vectors, film mode status, scene change decisions, and complexities along with the video data to MCTF. MCTF filters the video data using the motion vectors, film mode status, scene change decisions, and complexities from first pass encoder. MCTF then provides the filtered video data to second pass encoder, which also receives information from first pass encoder. 
         [0020]      FIG. 4  shows the block diagram of the overall encoding process for performing MPEG2 encoding. The architecture implements certain video compression features to improve the video quality over current video encoder products, including motion compensated temporal filtering that tracks motion of objects to efficiently remove video noise. There are five major components in the system—a first pass encoder  41 , a noise reduction preprocessor  42 , a second pass encoder  43 , a master compression controller  44  and a packet processor  45 . There is also a Video FIFO and Packet Creator  46 . 
         [0021]    The first pass encoder  41  is a 110 MHz C-Cube DvXpertII (E4) MPEG2 video compression chip. It performs the following functions:
       Captures the input CCIR656 (D1) format video;   Captures VBI data;   Horizontal filters and decimates the pixel data in ½, ⅔, or ¾ resolution mode;   Detelecines (inverse 3:2 pulldown);   Detects special events (scene changes);   Trial encodes the video frames to generate complexity measure for the frames; and   Determines P or B picture coding type of a frame.       
 
         [0029]    The noise reduction preprocessor  42  is a 300 MHz Equator MAPCA media processor chip. It performs the following functions:
       Lookahead frame delay;   Temporal filtering; and.   Spatial filtering.       
 
         [0033]    The motion compensated temporal filter is an IIR filter. The motion vectors generated by first pass encoding are used to trace motion of objects. 
         [0034]      FIG. 6  shows the block diagram of the MCTF. The forward temporal filter  61  processes each frame (P and B) in the video sequence using the forward motion vectors from the first pass encoder  41 . After a certain delay, the backward temporal filter  63  processes the B frames in the video sequence using the backward motion vectors from the first pass encoder  41 . If there are two B frames between two consecutive P frames, minimum two-frame delay is required between the forward  61  and backward  63  temporal filters. The backward temporal filter  63  does not process the P frames in the video sequence. 
         [0035]    For P frames, unidirectional temporal filtering is employed, the forward temporal filter generates the following output gp(i, j, k) for every pixel at indexes (i, j) on (x, y) coordinate and index k on time t axis. 
         [0000]        g   p ( i, j, k )= f   p ( i, j, k )±β 
         [0000]    where
       f p (i, j, k) is the input pixel of the current P frame at input time k.   β is the discrimination value of the filter and it is a function of the absolute difference between f p (i, j, k) and g p (i+d (x, j) ,j+d (y, f) , k−M−1).   (d (x, f) , d (y, f) ) is the forward motion vector at (i, j) with half-pel precision. g p (i+d (x, f) , j+d (y, f) , k−M−1) is the prior filtered P frame along the motion trajectories at input time k−M−1.   M is the number of B frames between the two consecutive P frames.       
 
         [0040]    For B frames, bi-directional temporal filtering is employed, the forward temporal filter generates the following output g′ B (i, j, k) for every pixel at indexes (i, j) on (x, y) coordinate and index k on time t axis. 
         [0000]        g′   B ( i, j, k )= f   B ( i, j, k )±β 1   
         [0041]    where
       f B (i, j, k) is the input pixel of the current B frame at input time k.       
 
         [0043]    For the first B frame between the two consecutive P frames, β 1  is the discrimination value of the filter and it is a function of the absolute difference between f B (i, j, k) and g p (i+d (x, f) , j+d (y, f) , k−1). 
         [0044]    For the second B frame between the two consecutive P frames, β 1  is the discrimination value of the filter and it is a function of the absolute difference between f B (i, j, k) and g p (i+d (x, f _, j+d (y, f) , k−2). 
         [0045]    (d (x, f) , d (y, f) ) is the forward motion vector at (i, j) with half-pel precision. g p (i+d (x, f) , j+d (y, f) , k−1) or g p (i+d (x, f) , j+d (y, f) , k−2) is the prior forward temporal filtered P frame along the motion trajectories at input time k−1 or k−2 respectively. 
         [0046]    For B frames, the backward temporal filter generates the final output g B (i, j, k) for every pixel at indexes (i, j) on (x, y) coordinate and index k on time t axis. 
         [0000]        g   B ( i, j, k )= g′   B ( i, j, k )±β 2   
         [0047]    where
       g′ B (i, j, k) is the input pixel to the backward temporal filter of the current B frame at input time k. g′ B (i, j, k) has been previously filtered by the forward temporal filter.       
 
         [0049]    For the first B frame, β 2  is the discrimination value of the filter and it is a function of the absolute difference between g′ B (i, j, k) and g p (i+d (x, b) , j+d (y, b) , k+2). 
         [0050]    For the second B frame, β 2  is the discrimination value of the filter and it is a function of the absolute difference between g′ B (i, j, k) and g p (i+d (x, b) , j+d (y, b) , k+1). 
         [0051]    (d (x, b) , d (y, b) ) is the backward field motion vector at (i, j) with half-pel precision. 
         [0052]    g p (i+d (x, b) , j+d (y, b) , k+1) or g p (i+d (x, b) , j+d (y, b) , k+2) is the prior forward temporal filtered P frame along the motion trajectories at input time k+1 or k+2 respectively. 
         [0053]    It is important to insert a minimum two-frame delay between the forward and backward temporal filters. The two-frame delay allows the forward temporal filter  61  process the P frame at input time k+1 or k+2 prior to the backward temporal filter  63  process the B frame at input time k. 
         [0054]      FIG. 7  shows the order of processing for the bi-directional filtering. The forward temporal filter processes each frame (both P and B) in the display order as shown by the arrows  1 ,  2 ,  3 , and  4 . The backward temporal filter processes the B frames in the display order as shown by the arrows  5  and  6 . It does not process the P frames. 
         [0055]    The second pass encoder  43  is a 110 MHz C-Cube DvXpertII (E4) MPEG2 video compression chip. It performs the following functions:
       Refines the motion vectors generated by the first pass encoding;   Calculates need parameters for statmux;   Performs MPEG2 Main Profile at Main Level video encoding;   Rate control; and   Determine if a frame should be coded as I-frame.       
 
         [0061]    The Master Compression Controller (MCC)  44  is a 50 MHz PPC860T PowerPC controller. It performs the following functions.
       Communicates with the Controller/Element Manager;   Functions as the host processor for the first pass encoder, second pass encoder and noise reduction preprocessor to perform configuration and status inquiry tasks;   Decodes the VBI data captured on the first pass encoder, delays the VBI decoded VBI data and passes it to the second pass encoder;   Communicates with the QLP to relay need parameter and bit rate messages to the second pass encode; and   Maintains the queue of bit rate values received from the QLP for buffer modeling.       
 
         [0067]    The Packet Processor is an 80 MHz DSP56309. It performs the following functions.
       Generates PTS and DTS;   Pulls bits out from the video FIFO to create MPEG2 packets at a fixed rate, or at a bit rate specified by the transmission bit rate parameter received from the QLP;   Models the decoder buffer and computes the minimum and maximum transmission bit rate constraints to send to the QLP in statmux; and   Performs other MPEG2 system multiplex functions.       
 
         [0072]    According to one aspect of the present invention, half-pel motion compensation is employed in the MCTF to improve the accuracy of the estimated motion trajectories. This increases the signal correlation to allow more noise reduction without trailing or blocking artifacts. 
         [0073]    According to another aspect of the present invention, the adaptation scheme of the MCTF is optimized. The optimal adaptation scheme is first determined by experiment. By adjusting the filter continuously at various levels using a series of images with variety of contents, the optimal adaptation area is identified. The “forbidden area” shown in  FIG. 1  represents the area in which filtering would degrade the image. The area in white represents the area in which potential noise that may be filtered without adversely affecting the image. The optimal adaptation allows more noise reduction without adversely affecting the signal content. In one version of MCTF, a linear adaptation scheme is employed. Another MCTF is developed, which employs a non-linear adaptation scheme to closely approximate the optimal scheme, and which also reduces the computational complexity if parameters are selected appropriately. 
       The Optimal Adaptation Scheme 
       [0074]    The MCTF is described by the following equation: 
         [0000]        g ( i, j, k )=α f ( i, j, k )+(1−α) g ( i−d   x   , j−d   y   , k− 1)  (1) 
         [0000]    where g(i−d x , j−d y , k−1) is the previous filtered frame; f(i, j, k) and g(i, j, k) are the current input and output frames, respectively; d x  and d y  are the horizontal and vertical motion displacements, respectively; the filter coefficient α adapts to the magnitude difference δ between f(i, j, k) and g(i−d x , j−d y , k−1) . The optimal adaptation scheme is established experimentally to allow maximum noise suppression without adversely affecting the video content. The optimal adaptation scheme is shown in  FIG. 1  in the α−δ domain along with the maximum strength linear adaptation scheme. From  FIG. 1 , we can see the optimal adaptation scheme is a non-linear function which depends on δ γ , where γ&lt;1. Obviously, the implementation of the optimal adaptation scheme is computationally intensive in the α−δ domain because of the need to calculate δ γ . Also, any adaptation scheme cannot travel through the shaded area. Otherwise, it will adversely affect the video content. 
         [0075]    While computationally intensive, the non-linear adaptation scheme will result in an improvement in noise reduction. Nevertheless, the following exemplary embodiment seeks to obtain this improvement while concomitantly decreasing the computation complexity. 
       Exemplary Embodiment 
       [0076]    An exemplary embodiment of a Motion Compensated Temporal Filter (MCTF) is developed in this section. Let δ represent the magnitude difference between f(i, j, k) and g(i−d x , j−d y ,k−1). 
         [0000]      δ=| f ( i, j, k )− g ( i−d   x   , j−d   y   , k− 1)| or  g ( i−d   y   j−d   y   , k− 1)= f ( i, j, k )±δ  (2) 
         [0077]    Substituting equation (2) into equation (1), we obtain: 
         [0000]        g ( i, j, k )= f ( i, j, k )±(1−α)δ= f ( i, j, k )±β  (3) 
         [0000]    where β represents the discrimination between the filtered and original frames. 
         [0000]      β=(1−α)δ  (4) 
         [0078]    Equation (3) represents a new implementation of a MCTF, where the discrimination β is first calculated. Adding or subtracting β from the input, we can then obtain the filter output. If f(i, j, k) is greater then g(i−d x , j−d y , k−1), subtraction is performed. Otherwise addition is performed. Comparing equation (3) to equation (1), we can see equation (1) consists of two multiplications while equation (3) consists of only one. Because a is a fixed-point number and the pixel values are 8-bit precision, the multiplications have to be performed with 16-bit precision. Although equation (3) has only one multiplication, the additional computation involved in the decision of whether adding or subtracting β from the input cancels out the savings. 
         [0079]    Equation (4) translates an adaptation scheme in the original α-δ domain into the new, β-δ domain. As will be shown, a non-linear adaptation in the α-δ domain is translated into a linear adaptation in the β-δ domain. This allows a non-linear adaptation in the original domain to be implemented as a linear adaptation in the new domain. This is the fundamental advantage of the exemplary embodiment. 
         [0080]    To get more insight into the exemplary embodiment, let us first consider the linear adaptation scheme depicted in  FIG. 1 . In the linear adaptation scheme, a is a linear function of δ. 
         [0000]    
       
         
           
             
               
                 
                   α 
                   = 
                   
                     { 
                     
                       
                         
                           
                             
                               
                                 ( 
                                 
                                   1 
                                   / 
                                   Δ 
                                 
                                 ) 
                               
                                
                               δ 
                             
                             , 
                           
                         
                         
                           
                             δ 
                             &lt; 
                             Δ 
                           
                         
                       
                       
                         
                           
                             0 
                             , 
                           
                         
                         
                           
                             δ 
                             ≥ 
                             Δ 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
           
         
       
     
         [0081]    Clearly, the bigger is Δ, the stronger is the filter. Substitute equation (5) into equation (4), we obtain: 
         [0000]    
       
         
           
             
               
                 
                   β 
                   = 
                   
                     { 
                     
                       
                         
                           
                             
                               δ 
                               - 
                               
                                 
                                   ( 
                                   
                                     1 
                                     / 
                                     Δ 
                                   
                                   ) 
                                 
                                  
                                 
                                   δ 
                                   2 
                                 
                               
                             
                             , 
                           
                         
                         
                           
                             δ 
                             &lt; 
                             Δ 
                           
                         
                       
                       
                         
                           
                             0 
                             , 
                           
                         
                         
                           
                             δ 
                             ≥ 
                             Δ 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0082]    Equation (6) describes the adaptation scheme in the new β-δ domain. Clearly, it is a non-linear function. In another word, a linear adaptation scheme in the α-δ domain (original implementation) is translated into a non-linear adaptation scheme in the β-δ domain (exemplary embodiment). This non-linear adaptation scheme in the β-δ domain is shown in  FIG. 2  along with a linear adaptation scheme. From equation (6), we can see that β attains its maximum value at 
         [0000]      δ m =(Δ/2). 
         [0083]    The maximum value is β m =(Δ/4). Clearly, to avoid blocking artifacts, we must choose small Δ. This limits the capability of the noise suppression. 
         [0084]    In accordance with another aspect of the present invention, a linear adaptation in the β-δ domain can be written as: 
         [0000]    
       
         
           
             
               
                 
                   β 
                   = 
                   
                     { 
                     
                       
                         
                           
                             
                               
                                 κδ 
                                 ′ 
                               
                               + 
                               χ 
                             
                             , 
                           
                         
                         
                           
                             
                               
                                 δ 
                                 ′ 
                               
                               = 
                               
                                 max 
                                  
                                 
                                   { 
                                   
                                     δ 
                                     , 
                                     
                                       δ 
                                       0 
                                     
                                   
                                   } 
                                 
                               
                             
                             , 
                           
                         
                         
                           
                             δ 
                             &lt; 
                             
                               δ 
                               1 
                             
                           
                         
                       
                       
                         
                           
                             0 
                             , 
                           
                         
                         
                           
                               
                           
                         
                         
                           
                             δ 
                             ≥ 
                             
                               δ 
                               1 
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   7 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where κ and X are constants less than one. Thus, equation (7) represents a linear calculation of the discrimination value for the filter, which discrimination value is then added or subtracted to/from the input frame to obtain the output frame. 
         [0085]    In accordance with another aspect of the present invention, β should be truncated down to an integer to avoid trailing or blocking artifacts. In general, because the constants κ and X are fixed numbers, equation (7) has to be implemented with 16-bit precision. 
         [0086]    From equation (4), we can see a linear adaptation in the β-δ domain is translated into a non-linear adaptation in the α-δ domain. The non-linear adaptation in the α-δ domain is shown in  FIG. 3  along with the optimal adaptation. From  FIG. 3 , we can see the non-linear adaptation closely approximates the optimal adaptation. Also we can see, the implementation of the non-linear adaptation in the α-δ domain is rather complex as compared to the implementation of the same adaptation in the β-δ domain, in which the adaptation becomes linear. 
         [0087]    In accordance with still another aspect of the present invention, from equation (7), filters with different strength can be implemented by selecting different slopes. This is another advantage over the non-linear adaptation of equation (6). From equation (6), the slope, at δ=0, is actually fixed regardless of the strength of the filter. A stronger filter can be implemented either by increasing the slope or by increasing the maximum discrimination β. Equation (7) allows us to increase one without changing the other. This allows the adaptation scheme be more manageable. 
         [0088]    In accordance with yet another aspect of the present invention, from equation (7), κ and X are used primarily to adjust the filter strength where δ is small, while δ 0  and δ 1  are used to adjust the filter strength where δ is relatively large. δ 0  is also used to adjust the maximum discrimination value to allow as much noise reduction as possible without disrupting the video image. Equation (7) allows us to adjust one aspect of the filter without adversely affecting other aspects of the filter. This allows the adaptation scheme be more manageable. 
         [0089]    An example of linear adaptation in the β-δ domain is given as: 
         [0000]    
       
         
           
             
               
                 
                   β 
                   = 
                   
                     { 
                     
                       
                         
                           
                             
                               
                                 ( 
                                 
                                   
                                     δ 
                                     ′ 
                                   
                                   + 
                                   3 
                                 
                                 ) 
                               
                               / 
                               4 
                             
                             , 
                           
                         
                         
                           
                             
                               
                                 δ 
                                 ′ 
                               
                               = 
                               
                                 max 
                                  
                                 
                                   { 
                                   
                                     δ 
                                     , 
                                     5 
                                   
                                   } 
                                 
                               
                             
                             , 
                           
                         
                         
                           
                             δ 
                             &lt; 
                             7 
                           
                         
                       
                       
                         
                           
                             0 
                             , 
                           
                         
                         
                           
                               
                           
                         
                         
                           
                             δ 
                             ≥ 
                             7 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   8 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where we choose κ=¼, X=¾, δ 0 =5, and δ 1 =7. 
         [0090]    From equation (8), we can see, the implementation of the linear adaptation, which maps to a non-linear adaptation in the a-d domain, can be performed with 8-bit precision. In other words, the exemplary embodiment of the non-linear adaptation in the α-δ domain not only closely approximates the optimal adaptation, but also reduces the computational complexity (note that the division can be implemented by right shifting). 
         [0091]    The optimal adaptation scheme is established experimentally to allow maximum noise suppression without adversely affecting the video content. A new implementation of MCTF is developed, which employs non-linear adaptation scheme to closely approximate the optimal scheme. The new implementation also reduces the computational complexity if parameters are selected appropriately. 
         [0092]    Turning to  FIG. 5 , shown therein is an exemplary embodiment of a method for filtering a video signal according to still another aspect of the present invention. 
         [0093]    In element  51 , the next frame in the sequence of frames is received by the processor. We shall call this frame the current frame. 
         [0094]    In element  52 , the magnitude difference between the current frame and a prior filtered frame along the motion trajectories in the sequence of frames is determined. 
         [0095]    In element  53 , a discrimination value is determined based on the magnitude difference. This discrimination value can be obtained from either equation (6) or equation (7), although equation (7) is considerably easier to implement and manage and it also allows more noise reduction without adversely affecting the video content. 
         [0096]    In element  54 , the current frame is modified using the discrimination value calculated above. The discrimination value is added to or subtracted from the current frame to obtain the filtered version of the current frame. 
         [0097]    As the discrimination value is modified based on the magnitude difference between the current frame and the immediately preceding frame before each frame is “filtered,” the discrimination value thereby adapts to the magnitude difference between frames. 
         [0098]    In element  55 , the strength of the filter is adjusted by modifying certain variables in the determination of the discrimination value. Various techniques can be employed for adjusting the strength of the filter, such as: (1) by modifying the slope of the equation calculating the discrimination value fi of the filter; (2) by modifying the maximum discrimination β; (3) in cases in which δ is small, by modifying the offset of the equation X; (4) in cases where δ is small by modifying the slope of the equation κ; (5) in cases where δ is relatively large by modifying the cutoff value of the equation δ 1 ; and (6) in cases where δ is relatively large by modifying the saturation value of the equation δ 0 . These are only some examples of how to adjust the strength of the filter. 
         [0099]    In step  56 , the maximum discrimination value is regulated by modifying the saturation value of the equation (e.g., δ 0  in some embodiments) to avoid trailing or blocking artifacts. 
         [0100]    Although various embodiments are specifically illustrated and described herein, it will be appreciated that modifications and variations of the invention are covered by the above teachings and are within the purview of the appended claims without departing from the spirit and intended scope of the invention. For example, certain filter values are depicted, however, others are possible depending upon the application without departing from the scope of the present invention. Furthermore, these examples should not be interpreted to limit the modifications and variations of the invention covered by the claims but are merely illustrative of possible variations.