Bit-accurate film grain simulation method based on pre-computed transformed coefficients

Information

  • Patent Grant
  • 9098916
  • Patent Number
    9,098,916
  • Date Filed
    Wednesday, October 26, 2005
    19 years ago
  • Date Issued
    Tuesday, August 4, 2015
    9 years ago
Abstract
Creation of a Bit-accurate film grain pattern for blending in an image block occurs by first establishing a set of bit-accurate transformed coefficients. The set of bit-accurate transformed coefficients undergo frequency filtering and a subsequent bit-accurate inverse transformation to yield the film grain pattern. The film grain pattern can then undergo blending with an image block to restore the look of film to the image.
Description
TECHNICAL FIELD

This invention relates to a technique for simulating film grain in an image.


BACKGROUND OF THE INVENTION

Motion picture films comprise silver-halide crystals dispersed in an emulsion, coated in thin layers on a film base. The exposure and development of these crystals form the photographic image consisting of discrete tiny particles of silver. In color negatives, the silver undergoes chemical removal after development and tiny blobs of dye occur on the sites where the silver crystals form. These small specks of dye are commonly called ‘grain’ in color film. Grain appears randomly distributed on the resulting image because of the random formation of silver crystals on the original emulsion. Within a uniformly exposed area, some crystals develop after exposure while others do not.


Grain varies in sizes and shapes. The faster the film, the larger the clumps of silver formed and blobs of dye generated, and the more they tend to group together in random patterns. The grain pattern is typically known as ‘granularity’. The naked eye cannot distinguish individual grains, which vary from 0.0002 mm to about 0.002 mm. Instead, the eye resolves groups of grains, referred to as blobs. A viewer identifies these groups of blobs as film grain. As the image resolution becomes larger, the perception of the film grain becomes higher. Film grain becomes clearly noticeable on cinema and high-definition images, whereas film grain progressively loses importance in SDTV and becomes imperceptible in smaller formats.


Motion picture film typically contains image-dependent noise resulting either from the physical process of exposure and development of the photographic film or from the subsequent editing of the images. The photographic film possesses a characteristic quasi-random pattern, or texture, resulting from physical granularity of the photographic emulsion. Alternatively, a similar pattern can be simulated over computed-generated images in order to blend them with photographic film. In both cases, this image-dependent noise is referred to as grain. Quite often, moderate grain texture presents a desirable feature in motion pictures. In some instances, the film grain provides visual cues that facilitate the correct perception of two-dimensional pictures. Film grain is often varied within a single film to provide various clues as to time reference, point of view, etc. Many other technical and artistic uses exist for controlling grain texture in the motion picture industry. Therefore, preserving the grainy appearance of images throughout image processing and delivery chain has become a requirement in the motion picture industry.


Several commercially available products have the capability of simulating film grain, often for blending a computer-generated object into a natural scene. Cineon® from Eastman Kodak Co, Rochester N.Y., one of the first digital film applications to implement grain simulation, produces very realistic results for many grain types. However, the Cineon® application does not yield good performance for many high-speed films because of the noticeable diagonal stripes the application produces for high grain size settings. Further, the Cineon® application fails to simulate grain with adequate fidelity when images are subject to previous processing, for example, such as when the images are copied or digitally processed.


Another commercial product that simulates film grain is Grain Surgery™ from Visual Infinity Inc., which is used as a plug-in of Adobe® After Effects®. The Grain Surgery™ product appears to generate synthetic grain by filtering a set of random numbers. This approach suffers from disadvantage of a high computational complexity.


None of these past schemes solves the problem of restoring film grain in compressed video. Film grain constitutes a high frequency quasi-random phenomenon that typically cannot undergo compression using conventional spatial and temporal methods that take advantage of redundancies in the video sequences. Attempts to process film-originated images using MPEG-2 or ITU-T Rec. H.264 ISO/IEC 14496-10 compression techniques usually either result in an unacceptably low degree of compression or complete loss of the grain texture.


Thus, there exists a need for a technique simulating film grain, especially a technique that affords relatively low complexity.


BRIEF SUMMARY OF THE INVENTION

Briefly, in accordance with the present principles, there is provided a method for simulating a film grain pattern. The method begins by obtaining a set of bit-accurate transformed coefficients. The set of bit-accurate transformed coefficients then undergoes filtering. Thereafter, the filtered set of bit-accurate transformed coefficients undergoes a bit-accurate inverse transform to yield a film grain pattern.





DETAILED DESCRIPTION OF THE DRAWINGS


FIG. 1 depicts a block schematic diagram of the combination of a transmitter and receiver in a film grain processing chain useful for practicing the technique of the present principles;



FIG. 2 depicts, in flow chart form, the steps of a method for creating a bit-accurate film grain pattern using a Gaussian random number generator;



FIG. 3 depicts, in flow chart form, the steps of a method for creating a bit-accurate film grain pattern using a Gaussian random number look up table;



FIG. 4 depicts, in flow chart form, the steps of a method for creating a bit-accurate film grain pattern creation using a single image of Gaussian noise;



FIG. 5 depicts, in flow chart form, the steps of a method for creating a bit-accurate film grain pattern creation using a set of pre-computed Discrete Cosine Transformed (DCT) coefficients of a single image of Gaussian noise; and



FIG. 6 depicts, in flow chart form, the steps of a method for creating a bit-accurate film grain pattern using pre-computed DCT coefficients of several images of Gaussian noise.





DETAILED DESCRIPTION

To understand the technique of the present principles for creating a bit-accurate film grain pattern, a brief overview of film grain simulation will prove helpful. FIG. 1 depicts a block schematic diagram of a transmitter 10, which receives an input video signal and, in turn, generates a compressed video stream at its output. In addition, the transmitter 10 also generates information indicative of the film grain (if any) present in the sample. In practice, the transmitter 10 could comprises part of a head-end array of a cable television system, or other such system that distributes compressed video to one or more downstream receivers 11, only one of which is shown in FIG. 1. The transmitter 10 could also take the form of encoder that presents media like DVDs. The receiver 11 decodes the coded video stream and simulates film grain in accordance with the film grain information and decoded video, both received from the transmitter 10 or directly from the media itself in the case of a DVD or the like, to yield an output video stream that has simulated film grain. The receiver 11 can take the form of a set-top box or other such mechanism that serves to decode compressed video and simulate film grain in that video.


The overall management of film grain requires the transmitter 10 (i.e., the encoder) provide information with respect to the film grain in the incoming video. In other words, the transmitter 10 “models” the film grain. Further the receiver 11 (i.e., decoder) simulates the film grain according to the film grain information received from the transmitter 10. The transmitter 10 enhances the quality of the compressed video by enabling the receiver 11 to simulate film grain in the video signal when difficulty exists in retaining the film grain during the video coding process.


In the illustrated embodiment of FIG. 1, the transmitter 10 includes a video encoder 12 which encodes the video stream using any of the well known video compression techniques such as the ITU-T Rec. H.264|ISO/IEC 14496-10 video compression standard. Optionally, a film grain remover 14, in the form of a filter or the like depicted in dashed lines in FIG. 1, could exist upstream of the encoder 12 to remove any film grain in the incoming video stream prior to encoding. To the extent that the incoming video contains no film grain, no need would exist for the film grain remover 14.


A film grain modeler 16 accepts the input video stream, as well as the output signal of the film grain remover 14 (when present). Using such input information, the film grain modeler 16 establishes the film grain in the incoming video signal. In its simplest form, the film grain modeler 16 could comprise a look up table containing film grain models for different film stocks. Information in the incoming video signal would specify the particular film stock originally used to record the image prior to conversion into a video signal, thus allowing the film grain modeler 16 to select the appropriate film grain model for such film stock. Alternatively, the film grain modeler 16 could comprise a processor or dedicated logic circuit that would execute one or more algorithms to sample the incoming video and determine the film grain pattern that is present.


The receiver 11 typically includes a video decoder 18 that serves to decode the compressed video stream received from the transmitter 10. The structure of the decoder 18 will depend on the type of compression performed by the encoder 12 within the transmitter 10. Thus, for example, the use within the transmitter 10 of an encoder 12 that employs the ITU-T Rec. H.264|ISO/IEC 14496-10 video compression standard to compress outgoing video will dictate the need for an H.264-compliant decoder 18. Within the receiver 11, a film grain simulator 20 receives the film grain information from the film grain model 16. The film grain simulator 20 can take the form of a programmed processor, or dedicated logic circuit having the capability of simulating film grain for combination via a combiner 22 with the decoded video stream.


Film grain simulation aims to synthesize film grain samples that simulate the look of the original film content. As described, film grain modeling occurs at the transmitter 10 of FIG. 1, whereas film grain simulation occurs at the receiver 11. In particular, film grain simulation occurs in the receiver 11 along with the decoding of the incoming video stream from the transmitter 10 upstream of the output of the decoded video stream. Note that the decoding process that occurs in the receiver 11 makes no use of images with added film grain. Rather, film grain simulation constitutes a post-processing method for synthesizing simulated film grain in the decoded images for display. For that reason, the ITU-T Rec. H.264|ISO/IEC 14496-10 video compression standard contains no specifications regarding the film grain simulation process. However, film grain simulation requires information concerning the grain pattern in the incoming video signal, which information typically undergoes transmission in a Supplemental Enhancement Information (SEI) message when using the ITU-T Rec. H.264|ISO/IEC 14496-10 video compression standard as specified by the Amendment 1 (Fidelity Range Extensions) of that compression standard.


In accordance with the present principles, disclosed herein are several implementations of a method for creating a bit-accurate film grain pattern, such as for storage in a database for film grain simulation in HD DVD systems. Such a film grain pattern database enables film grain simulation in accordance with the parameters in a film grain Supplemental Enhancement Information (SEI) message according to the H.264|MPEG-4 AVC standard as discussed above. More specifically, devices and systems that implement the existing or future HD DVD specifications can employ such a film grain pattern database to store samples of different film grain types, each characterized by different frequency model parameters.


As discussed in detail below, the method of the present principles for creating a bit-accurate pattern of film grain makes use of a basic strategy that enables implementation in several different ways within the receiver 11. The basic strategy for creating a bit-accurate film grain pattern begins by establishing a set of integer transformed coefficients, typically, although not necessarily, Discrete Cosine Transformed (DCT) coefficients, typically N×N in size where N is an integer. The term “image” will some times refer to such a set of integer transformed coefficients. The step of establishing a set or image of integer transformed coefficients can occur in several different ways. For example, the establishing step could entail the processor or discrete logic circuitry in the receiver 11 of FIG. 1 accessing a database of integer-transformed coefficients as described hereinafter. Alternatively, the processor or discrete logic circuitry in the receiver 11 could establish a set or image of Gaussian random numbers, by either reading a database of such numbers or generating them directly. Thereafter, an integer DCT would be performed on the Gaussian random number image to yield a set of integer DCT coefficients.


The next step in the basic process of creating a bit-accurate film grain pattern entails frequency filtering the integer-transformed coefficients according to a desired film grain size and shape. For example, the frequency filtering could employ a predefined set of cut frequencies fHL, fVL, fHH, and fVH that represent cut-off frequencies (in two dimensions) of a filter that characterizes the desired film grain pattern. Following the frequency filtering, the integer transformed coefficients undergo an inverse transform to create the bit-accurate film grain pattern for subsequent storage in a database. In some instances, scaling could occur following the inverse transform of the coefficients.


In practice, the film grain pattern creation method of the present principles makes use of integer DCT and inverse DCT operations, thus assuring bit accuracy which proves useful for verification and testing purposes. However, different implementations of the basic strategy as described below do not necessarily afford bit accuracy between each other. Choosing a particular implementation will depend on the desired tradeoff between memory requirements and computational cost.



FIG. 2 illustrates a first implementation of a method in accordance with the present principles for creating a bit-accurate pattern of film grain. The method of FIG. 2 commences upon execution of the start step 100 during which initialization occurs, although such initialization need not necessarily happen. Next step 101 occurs, initiating entry into a loop that repeats for all possible film grain sizes and shapes. Upon entry into the loop, step 102 occurs to achieve creation of an image of random Gaussian values. In practice, the Gaussian random image has a size of N×N values. In the method of FIG. 2, creation of the image of Gaussian random values occurs by generating such values via a Gaussian random number generator 103.


The image of Gaussian random values created during step 102 then undergoes a bit-accurate transform, typically by way of an integer Discrete Cosine Transform (DCT) operation during step 104, although other bit-accurate transform techniques exist. The integer DCT coefficients established during step 102 then undergo frequency filtering during step 106. Typically, the frequency filter occurs using a predefined set of cut frequencies fHL, fVL, fHH, and fVH that represent cut-off frequencies (in two dimensions) of a filter that characterizes the desired film grain pattern. Next step 108 occurs during which the frequency filtered, transformed coefficients undergo an inverse bit-accurate transformation, typically, although not necessarily, an integer Inverse Discrete Cosine Transformation (IDCT), to yield a bit-accurate film grain pattern. Under some circumstances, scaling of the inversely transformed block of coefficient following step 108 will prove useful to vary the intensity of the bit-accurate film grain pattern.


Following step 108, step 110 occurs and the bit-accurate film grain that results from the integer IDCT transformation (and scaling if performed) undergoes storage in a database 111. Steps 102-110 undergo repeating during execution of the loop initiated during step 101. The loop continues for all possible film grain sizes and shapes, whereupon loop execution ends during step 112, after which the method ends during step 114.


To achieve bit accuracy, the method of FIG. 2 requires a bit-accurate Gaussian random number generator and bit-accurate DCT and Inverse DCT transform operations. Achieving bit accuracy among different implementations of the present method requires the use of a bit-accurate Gaussian random number generator. Which particular Gaussian random number generator is used remains unimportant as long as all implementations use the same generator. As example, one could use the integer approximation of a bit-accurate Gaussian generator described in the publication “Numerical Recipes in C: The Art of Scientific Computing” (ISBN 0-521-41108-5), which is based on the Box-Muller transformation.


In addition to a bit-accurate Gaussian random number generator, a bit-accurate implementation of a Discrete Cosine Transform remains necessary to provide bit accuracy for the whole system. The illustrated embodiments described herein make use of an integer approximation of the DCT. The use of any integer approximation remains possible if the result obtained by the integer transformation lies reasonably close to the result obtained by the floating point DCT. This requirement ensures that the filtering performed in the transformed domain constitutes frequency filtering.


In the illustrated embodiment, an integer version of the DCT occurs by scaling the floating point version of the DCT. The N×N DCT matrix can be computed as:

















for ( i = 0; i < N; j++)









for (j = 0; j < N; j++)









{









if(i == 0)k = sqrt(1/N);










else
 k = sqrt(2/N);









C(i,j) = k * cos( ((2*j+1)*i*PI) / (2*N));



}











where C(i,j) represents the coefficient of the transformation matrix at row j and column i. Then, the integer matrix is computed as:

















for ( i = 0; i < N; j++)









for (j = 0; j < N; j++)









{









Cint(i,j) = round( C(i,j) * scaling_factor)









}











where round(x) returns the nearest integer approximation of x, and the scaling_factor is a positive integer value. The value of the scaling factor determines the quality of the integer approximation of the DCT (the larger the scaling factor, the better the approximation) as well as the number of bits required to compute the transform (the smaller the scaling factor, the lesser bits are required).


Once the integer approximation of the DCT has been defined, the DCT transform is computed as follows:

B=((CintT*b*Cint)+2scalingfactor)>>2*scaling_factor

where CintT denotes the transposed version of the transformation matrix. Analogously, the inverse transform is computed as follows:

b=((Cint*B*CintT)+2scalingfactor)>>2*scaling_factor

In a particular embodiment where an 8×8 DCT is used, the integer approximation would be:







Cint
8
τ

=

(



6


6


6


6


6


6


6


6




8


7


4


2



-
2




-
4




-
7




-
8





7


3



-
3




-
7




-
7




-
3



3


7




7



-
2




-
8




-
4



4


8


2



-
7





6



-
6




-
6



6


6



-
6




-
6



6




4



-
8



2


7



-
7




-
2



8



-
4





3



-
7



7



-
3




-
3



7



-
7



3




2



-
4



7



-
8



8



-
7



4



-
2




)






with scale factor equal to 16.



FIG. 3 illustrates a second implementation of a method in accordance with the present principles for creating a bit-accurate pattern of film grain. The method of FIG. 3 commences upon execution of the start step 200 during which initialization occurs, although such initialization need not necessarily happen. Next step 201 occurs, initiating entry into a loop that repeats for all possible film grain sizes and shapes. Upon entry into the loop, step 202 executes to create a set or image of random Gaussian values. In practice, the Gaussian random image has a size of N×N values. Creation of the image of Gaussian random values occurs during execution of the method of FIG. 3 by reading a set of Gaussian random values from a Gaussian random number look-up table (LUT) 203.


The image of Gaussian random values obtained during step 202 then undergoes a bit-accurate transform, typically by way of an integer Discrete Cosine Transform (DCT) operation during step 204, although other bit-accurate transform techniques exist. The integer transformed coefficients established during step 202 then undergo frequency filtering during step 206. Typically, the frequency filter occurs using a predefined set of cut frequencies fHL, fVL, fHH, and fVH that represent cut-off frequencies (in two dimensions) of a filter that characterizes the desired film grain pattern. Next step 208 occurs during which the frequency filtered block of transformed coefficients undergoes an inverse bit-accurate transformation, typically, although not necessarily, an integer Inverse Discrete Cosine Transformation (IDCT), to yield a bit-accurate film grain pattern. Under some circumstances, scaling of the frequency filtered, inverse transformed block of coefficients will prove useful.


Thereafter, step 210 occurs and the film grain pattern resulting from the integer inverse transformation (and scaling if performed) undergoes storage in a database 211. The steps 202-210 within the loop initiated during step 201 undergo repeating for all possible film grain sizes and shapes, whereupon loop execution ends during step 212, after which the method ends during step 214. In this way, the database 211 stores a plurality of film grain patters for future use in simulating film grain in a video signal.


As compared to the method of FIG. 2, the film grain simulation method of FIG. 3 obviates the use of a Gaussian random number generator. Instead, the method of FIG. 3 makes use of the look-up table 203 that contains pre-computed Gaussian random values. The approach simplifies the required hardware needed for implementation.



FIG. 4 illustrates a third implementation of a method in accordance with the present principles for creating a bit-accurate pattern of film grain. The method of FIG. 4 commences upon execution of the start step 300 during which initialization occurs, although such initialization need not necessarily happen. Next, step 302 occurs to create a set or image of random Gaussian values. In practice, the Gaussian random image has a size of N×N values. Creation of the image of Gaussian random values occurs in the method of FIG. 4, by reading a set of Gaussian random values from a Gaussian random number look up table (LUT) 303. The method used to create the Gaussian random number LUT 303 does not require bit accuracy. As long as all implementations of the film grain pattern database creation use the same LUT, the result assures bit accuracy. However, when using a bit-accurate method to create such a LUT, then LUT generation can occur at the beginning of the film grain pattern database creation process, avoiding the need for permanent storage of the LUT. Given that a single transform (step 304) needs to be computed for the generation of the entire film grain pattern database, the randomness of the film grain patterns stored in the database is reduced when compared to the result obtained with previously illustrated implementations (FIGS. 2 and 4). However, the third implementation illustrated in FIG. 4 also has lower computational needs, which proves useful for hardware and/or real-time implementations. The image of Gaussian random values created during step 302 then undergoes a bit-accurate transform, typically by way of an integer Discrete Cosine Transform (DCT) operation during step 304; although other bit-accurate transform techniques exist. Following step 304, step 305 occurs which initiates entry into a loop that repeats for all possible film grain sizes and shapes. Step 306, the first step within the loop, initiates frequency filtering of the integer DCT coefficients established during step 304. Typically, the frequency filter occurs using a predefined set of cut frequencies fHL, fVL, fHH, and fVH that represent cut-off frequencies (in two dimensions) of a filter that characterizes the desired film grain pattern. Following step 306, the frequency filtered integer DCT coefficients undergoes an inverse bit-accurate transformation, typically, although not necessarily, an integer Inverse Discrete Cosine Transformation (IDCT), during step 308 to yield a bit-accurate film grain pattern. Under some circumstances, scaling of film grain pattern created from the inversely transformed coefficients following step 308 will prove useful. Thereafter, step 310 occurs and the film grain pattern undergoes storage in a database 311.


Steps 306-310 within the loop initiated during step 305 undergo for all possible film grain sizes and shapes, whereupon loop execution ends during step 312, after which the method ends during step 314. In this way, the database 311 stores film grain patterns for all sizes and shapes of film grain.



FIG. 5 illustrates a fourth implementation of a method in accordance with the present principles for creating a bit-accurate pattern of film grain. The method of FIG. 5 commences upon execution of the start step 400 during which initialization occurs, although such initialization need not necessarily happen. Next, step 402 occurs, during which an image or set of DCT coefficients is read from a look-up table 403 created from an image of Gaussian random values. In practice, the image of DCT coefficients read from the LUT 403 has a size of N×N. Thus, unlike the previously described implementations of FIGS. 2-4, the film grain simulation method of FIG. 5 obviates the need to perform a separate integer DCT operation since the values read from the LUT 403 have already undergone such a transformation prior to loading into the LUT.


Following step 402, step 405 occurs which initiates entry into a loop that repeats for all possible film grain sizes and shapes. Step 406, the first step within the loop, initiates frequency filtering of the integer DCT coefficients obtained during step 402. Typically, the frequency filter occurs using a predefined set of cut frequencies fHL, fVL, fHH, and fVH that represent cut-off frequencies (in two dimensions) of a filter that characterizes the desired film grain pattern. Following step 406, the frequency filtered integer DCT coefficients undergo an inverse bit-accurate transformation, typically, although not necessarily, an integer Inverse Discrete Cosine Transformation (IDCT), during step 408 to yield a pattern of film grain. Under some circumstances, scaling of the inversely transformed block of coefficients following step 408 will prove useful. Thereafter, step 410 occurs and the film grain pattern resulting from the integer IDCT transformation (and scaling if performed) undergoes storage in a database 411.


Steps 406-410 within the loop initiated during step 405 undergo repeating for all possible film grain sizes and shapes, whereupon loop execution ends during step 412, after which the method ends during step 414. In this way, the database 411 stores film grain patterns for all sizes and shapes of film grain.


The method of FIG. 5 makes use of a single pre-computed block of transformed coefficients to generate the film grain patterns that populate the database 411. Compared to the previous implantations described in connection with FIGS. 2-4, implementation of FIG. 5 eliminates the step of performing an integer DCT transform, thus reducing computational cost.



FIG. 6 illustrates a fifth implementation of a method in accordance with the present principles for creating a bit-accurate pattern of film grain. The method of FIG. 6 commences upon execution of the start step 500 during which initialization occurs, although such initialization need not necessarily happen. Next, step 501 occurs, initiating entry into a loop that repeats for all possible film grain sizes and shapes. Step 502, the first step in the loop initiates reading of an image of set of integer DCT coefficients from a look up table (LUT) 503. In practice, the image of DCT coefficients has a size of N×N. Like the implementation of FIG. 5, the implementation of FIG. 6 establishes an image of integer transformed coefficients image by obtaining a set of integer DCT coefficients from the LUT 503. In practice, the values within the LUT 503 originate from an image of Gaussian random values that subsequently undergo an integer DCT. Thus, unlike the previously described implementations of FIGS. 2-4, the implementation of FIG. 6, like that of FIG. 5, obviates the need to perform a separate integer DCT since the values read from the look-up table 503 have already undergone such a transformation prior to loading in the table.


Following step 502, step 506 occurs during which the integer DCT coefficients undergo frequency filtering. Typically, the frequency filter occurs using a predefined set of cut frequencies fHL, fVL, fHH, and fVH that represent cut-off frequencies (in two dimensions) of a filter that characterizes the desired film grain pattern. Following step 506, the frequency filtered integer DCT coefficients undergo an inverse bit-accurate transformation, typically, although not necessarily, an integer Inverse Discrete Cosine Transformation (IDCT), to yield a block of film grain during step 508. Under some circumstances, scaling of the inversely transformed block of coefficient following step 508 will prove useful. Thereafter, step 510 occurs and the block of film grain that results from the integer IDCT transformation (and scaling if performed) undergoes storage in a database 511.


Steps 502-510 within the loop initiated during step 501 undergoes repeating for all possible film grain sizes and shapes, whereupon loop execution ends during step 512, after which the method ends during step 514. In this way, the database 511 stores film grain patterns for all sizes and shapes of film grain.


The implementation of film grain pattern creation described with respect to FIG. 6, like the implementation of FIGS. 2 and 3, makes use of different sets of integer DCT coefficients. Using different sets of integer DCT coefficients allows for creation of a richer database of film blocks patterns as compared to the implementation of FIGS. 1, 4 and 5, which utilizes a single block of integer transformed coefficients. While the implementation described with respect to FIGS. 1, 4 and 5 affords reduced storage and computational requirements, such implementations lower afford performance because the same noise pattern serves as the basis for generating all film grain patterns.


The foregoing describes several different implementations of a technique for creating at least one, and preferably a plurality of film grain patterns for simulating film grain in a video signal.

Claims
  • 1. A method of creating a pattern of bit-accurate film grain, comprising the steps of: (a) establishing a set of bit-accurate transformed coefficients by scaling a floating point version of a Discrete Cosine Transformed (DCT) matrix of coefficients;(b) frequency filtering the matrix of bit-accurate transformed coefficients; and(c) performing a bit-accurate inverse transformation on the frequency filtered transformed coefficients to yield a bit-accurate film grain pattern for blending in an image; and(d) automatically repeating steps (a) - (c) for different sizes and shapes of film grain.
  • 2. The method according to claim 1 further comprising step of (e) storing the film grain pattern in a database.
  • 3. The method according to claim 1 wherein the matrix of DCT coefficients is obtained from a database of stored values.
  • 4. The method according to claim 1 wherein the matrix of DCT coefficients is obtained by performing an integer Discrete Cosine Transformation on an image of Gaussian random numbers.
  • 5. The method according to claim 1 wherein the frequency filtering step makes use of a predefined set of cut frequencies fHL, fVL, and fHH, and fVH that represent cut-off frequencies, in two dimensions of a filter that characterizes the desired film grain pattern.
  • 6. The method according to claim 1 wherein the step of performing a bit-accurate inverse transformation on the frequency filtered transformed coefficients comprises the step of performing an integer inverse Discrete Cosine Transform.
  • 7. Apparatus comprising one of a processor and dedicated logic circuit for creating a pattern of bit-accurate film grain by (a) establishing a set of bit-accurate transformed coefficients by scaling a floating point version of a Discrete Cosine Transformed (DCT) matrix of coefficients; (b) frequency filtering the matrix of bit-accurate transformed coefficients; (c) performing a bit-accurate transformation on the frequency filtered transformed coefficients to yield a bit-accurate film grain pattern for blending in an image, and (d) automatically and repeating steps (a)-(c) for different possible sizes and shapes of film grain.
  • 8. The apparatus according to claim 7 further comprising a first memory for storing the film grain pattern.
  • 9. The apparatus according to claim 8 further comprising a second memory for storing at least one set of the integer Discrete Cosine Transformed coefficients.
  • 10. The apparatus according to claim 7 further comprising a second memory for storing a plurality of sets of the integer Discrete Cosine Transformed coefficients.
  • 11. The apparatus according to claim 7 further comprising: a Gaussian random generator for generating a Gaussian random number image;means for performing an integer Discrete Cosine transformation on the Gaussian random number image to yield the set of set of Discrete Cosine bit-accurate transformed coefficients.
  • 12. Apparatus for creating a pattern of bit-accurate film grain, comprising: means for establishing a set of bit-accurate transformed coefficients by scaling a floating point version of a Discrete Cosine Transformed (DCT) matrix of coefficients;means for frequency filtering the matrix of bit-accurate transformed coefficients; andmeans for automatically performing a bit-accurate transformation on the frequency filtered transformed coefficients to yield a bit-accurate film grain pattern for blending in an image for different sizes and shapes of film grain.
CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims the benefit, under 35 U.S.C. §365 of International Application PCT/US2005/038802 filed Oct. 26, 2005 which was published in accordance with PCT Article 21(2) on May 26, 2006 in English and which claims the benefit of U.S. provisional patent application No.'s 60/628,837 filed Nov. 17, 2004.

PCT Information
Filing Document Filing Date Country Kind 371c Date
PCT/US2005/038802 10/26/2005 WO 00 5/11/2007
Publishing Document Publishing Date Country Kind
WO2006/055208 5/26/2006 WO A
US Referenced Citations (117)
Number Name Date Kind
4897775 Kaassens Jan 1990 A
4935816 Faber Jun 1990 A
4998167 Jaqua Mar 1991 A
5028280 Ihara et al. Jul 1991 A
5140414 Mowry Aug 1992 A
5216556 Steinberg et al. Jun 1993 A
5262248 Ihara et al. Nov 1993 A
5285402 Keith Feb 1994 A
5285482 Sehier et al. Feb 1994 A
5335013 Faber Aug 1994 A
5374954 Mowry Dec 1994 A
5406326 Mowry Apr 1995 A
5453840 Parker et al. Sep 1995 A
5457491 Mowry Oct 1995 A
5461596 Barrett Oct 1995 A
5471572 Buchner et al. Nov 1995 A
5475425 Przyborski et al. Dec 1995 A
5526446 Adelson et al. Jun 1996 A
5550815 Cloonan et al. Aug 1996 A
5629769 Cookingham et al. May 1997 A
5641596 Gray et al. Jun 1997 A
5659382 Rybczynski Aug 1997 A
5687011 Mowry Nov 1997 A
5706361 Kent et al. Jan 1998 A
5709972 Cookingham et al. Jan 1998 A
5715008 Sekiguchi et al. Feb 1998 A
5742892 Chaddha Apr 1998 A
5751398 Beard May 1998 A
5767860 Zimmer et al. Jun 1998 A
5805782 Foran Sep 1998 A
5817447 Yip Oct 1998 A
5831673 Przyborski et al. Nov 1998 A
5845017 Keyes Dec 1998 A
5917609 Breeuwer et al. Jun 1999 A
6067125 May May 2000 A
6134200 Timmermans Oct 2000 A
6216838 Bacher Apr 2001 B1
6219838 Cherichetti et al. Apr 2001 B1
6233647 Bentz et al. May 2001 B1
6269180 Sevigny Jul 2001 B1
6285711 Ratakonda et al. Sep 2001 B1
6327304 Miller et al. Dec 2001 B1
6327391 Ohnishi et al. Dec 2001 B1
6370192 Pearlstein et al. Apr 2002 B1
6373992 Nagao Apr 2002 B1
6441918 Hori Aug 2002 B1
6459699 Kimura et al. Oct 2002 B1
6496221 Wolf et al. Dec 2002 B1
6559849 Anderson et al. May 2003 B1
6587509 Suzuki et al. Jul 2003 B1
6650327 Airey et al. Nov 2003 B1
6665369 Ukita Dec 2003 B2
6667815 Nagao Dec 2003 B1
6724942 Arai Apr 2004 B1
6728317 Demos Apr 2004 B1
6744928 Juri et al. Jun 2004 B1
6839152 Fan et al. Jan 2005 B2
6868190 Morton Mar 2005 B1
6963668 Engeldrum et al. Nov 2005 B2
6990251 Edgar Jan 2006 B2
6990252 Shekter Jan 2006 B2
6995793 Albadawi et al. Feb 2006 B1
7065255 Chen et al. Jun 2006 B2
7092016 Morton et al. Aug 2006 B2
7106907 Lin et al. Sep 2006 B2
7245783 Fielding Jul 2007 B2
7286565 Carr Oct 2007 B1
7362911 Frank Apr 2008 B1
7596239 Winger et al. Sep 2009 B2
7630005 Takei Dec 2009 B2
7653132 Dang Jan 2010 B2
7664337 Balram et al. Feb 2010 B2
7738722 Gomila et al. Jun 2010 B2
7742655 Gomila et al. Jun 2010 B2
7856055 Zhou et al. Dec 2010 B2
7958532 Paul et al. Jun 2011 B2
20010056568 Hirotsu et al. Dec 2001 A1
20020003903 Engeldrum et al. Jan 2002 A1
20020016103 Behnke Feb 2002 A1
20020034337 Shekter Mar 2002 A1
20020133764 Wang Sep 2002 A1
20020154140 Tazaki Oct 2002 A1
20020163657 Bogdanowicz Nov 2002 A1
20020171649 Fogg Nov 2002 A1
20030011615 Tidwell Jan 2003 A1
20030031128 Kim et al. Feb 2003 A1
20030043922 Kalker Mar 2003 A1
20030063778 Rowe et al. Apr 2003 A1
20030068097 Wilson et al. Apr 2003 A1
20030086623 Berkner et al. May 2003 A1
20030101453 Matsuyama May 2003 A1
20030206231 Chen et al. Nov 2003 A1
20030206662 Avinash Nov 2003 A1
20030218610 Mech et al. Nov 2003 A1
20040013308 Jeon et al. Jan 2004 A1
20040073926 Nakamura et al. Apr 2004 A1
20040101059 Joch et al. May 2004 A1
20060082649 Gomila et al. Apr 2006 A1
20060083316 Cooper et al. Apr 2006 A1
20060083426 Cooper et al. Apr 2006 A1
20060083429 Joly Apr 2006 A1
20060104366 Huang et al. May 2006 A1
20060182183 Winger Aug 2006 A1
20060183275 Schoner Aug 2006 A1
20060256853 Schlockermann et al. Nov 2006 A1
20060291557 Tourapis Dec 2006 A1
20070002947 Lu et al. Jan 2007 A1
20070030996 Winger et al. Feb 2007 A1
20070036452 Llach et al. Feb 2007 A1
20070047658 Tourapis et al. Mar 2007 A1
20070058866 Boyce et al. Mar 2007 A1
20070058878 Gomilla et al. Mar 2007 A1
20070070241 Boyce et al. Mar 2007 A1
20070104380 Gomila et al. May 2007 A1
20070117291 Cooper et al. May 2007 A1
20070297515 Gomila et al. Dec 2007 A1
20080252781 DeWaele et al. Oct 2008 A1
Foreign Referenced Citations (62)
Number Date Country
1530252 Sep 2004 CN
0364285 Apr 1990 EP
0622000 Oct 1992 EP
0575006 Dec 1993 EP
0622000 Mar 2000 EP
1175091 Jan 2002 EP
1215624 Jun 2002 EP
1511320 Mar 2005 EP
2312124 Oct 1997 GB
1156069 Jun 1989 JP
3187661 Aug 1991 JP
4097681 Mar 1992 JP
5252459 Sep 1993 JP
8079765 Mar 1996 JP
9062718 Mar 1997 JP
9139940 May 1997 JP
9247681 Sep 1997 JP
10509297 Sep 1998 JP
11 250246 Sep 1999 JP
200041242 Aug 2000 JP
2001357090 Dec 2001 JP
2001357090 Dec 2001 JP
2001357095 Dec 2001 JP
2002057719 Feb 2002 JP
2002374541 Dec 2002 JP
2003024326 Jan 2003 JP
200324326 Jan 2003 JP
2003163853 Jun 2003 JP
2003179923 Jun 2003 JP
2004120057 Apr 2004 JP
2004135169 Apr 2004 JP
2005080301 Mar 2005 JP
2007507172 Mar 2007 JP
2007521573 Aug 2007 JP
2007529945 Oct 2007 JP
2073913 Sep 1991 RU
2088962 Aug 1997 RU
2139637 Oct 1999 RU
2139637 Jan 2001 RU
WO9314591 Jul 1993 WO
WO 9404960 Mar 1994 WO
WO9520292 Jul 1995 WO
WO9710676 Mar 1997 WO
WO9722204 Jun 1997 WO
WO9841026 Sep 1998 WO
WO0018109 Mar 2000 WO
WO0146992 Jun 2001 WO
WO0174064 Oct 2001 WO
WO0177871 Oct 2001 WO
WO0233958 Apr 2002 WO
WO0251160 Jun 2002 WO
WO 03005731 Jan 2003 WO
WO2004077348 Sep 2004 WO
WO2004095829 Nov 2004 WO
WO2004095829 Nov 2004 WO
WO2004104931 Dec 2004 WO
WO2005027045 Mar 2005 WO
W02005032143 Apr 2005 WO
WO2005032143 Apr 2005 WO
WO2005039188 Apr 2005 WO
WO2006022705 Mar 2006 WO
WO2006057703 Jun 2006 WO
Non-Patent Literature Citations (94)
Entry
Oktem et al., “Transform Domain Algorithm for Reducing Effect of Film-Grain Noise in Image Compression” Electrnonic Letters, Oct. 14, 1999, vol. 35, No. 21.
Auto FX Software: Dreamsuite Series Two, Film Grain, http://web.archive.org/web/20040805085520/www.autofx.com/dreamsuite2/effect—pages/filmgrain.html, Aug. 5, 2004.
“SEI message for film grain encoding: syntax and results”, Christina Gomila, Sep. 2-5, 2003.
Office Action from U.S. Appl. No. 10/556,834 mailed Aug. 19, 2008.
Office Action from U.S. Appl. No. 10/556,834 mailed Feb. 20, 2009.
Notice of Non-Compliant Amendment from U.S. Appl. No. 10/556,834 mailed Jun. 16, 2009.
Office Action from U.S. Appl. No. 10/569,318 mailed Aug. 15, 2008.
Office Action from U.S. Appl. No. 10/569,318 mailed Feb. 13, 2009.
Office Action from U.S. Appl. No. 10/569,318 mailed Jul. 31, 2009.
Office Action from U.S. Appl. No. 10/572,820 mailed Jun. 1, 2009.
Office Action from U.S. Appl. No. 10/571,148 mailed Apr. 15, 2009.
Office Action from U.S. Appl. No. 10/575,676 mailed Dec. 15, 2008.
Office Action from U.S. Appl. No. 10/575,676 mailed Mar. 13, 2009.
Office Action from U.S. Appl. No. 11/246,848 mailed Jun. 26, 2009.
Office Action from U.S. Appl. No. 11/252,177 mailed May 29, 2008.
Office Action from U.S. Appl. No. 11/252,177 mailed Nov. 5, 2008.
Office Action from U.S. Appl. No. 11/252,177 mailed Apr. 17, 2009.
Office Action from U.S. Appl. No. 11/252,177 mailed Sep. 18, 2009.
Al-Shaykh et al, “Lossy Compression of Images Corrupted by Film Grain Noise,” School of Electrical and Computer Engineering, 1996 IEEE.
Al-Shaykh et al. “Restoration of Lossy Compressed Noisy Images,” IEEE Transactions on Image Processing, vol. 8, No. 10, Oct. 1999.
Al-Shaykh, “Lossy Compression of Noisy Images,” IEEE Transactions on Image Processing, vol. 7, No. 12, Dec. 1998.
Brightwell et al., “Automated Correction of Film Unsteadiness, Dirt and Grain,” International Broadcasting Convention, Sep. 16-20, 1994, Conference Publication No. 397, IEE, 1994.
Campisi et al, “Signal-Dependent Film Grain Noise Generation Using Homomorphic Adaptive Filtering,” IEE Proceedings, Image Signal Process, vol. 147, No. 3, Jun. 2000.
Chao et al, “Constant Quality Rate Control for Streaming MPEG-4 Fgs. Video,” Integrated Media Systems Center and Department of Electrical Engineering, University of Southern California, 2000 IEEE.
Chavel et al., “Film Grain Noise in Partially Coherent Imaging,” Optical Engineering, vol. 19, No. 3, May-Jun. 1980.
Fischer et al, “Image Sharpening Using Permutation Weighted Medians,” Department of Electrical Engineering, University of Delaware, fischer@ee.udel.edu, paredesj@ee.udel.edu, arce@ee.udel.edu.
Gomila, “SEI Message for Film Grain Encoding,” XP-002308742, Joint Video Team (JVT) of ISO/IEC MPEG & ITU-T VCEG, Document: JVT-H022, 8th Meeting, May 23-27, 2003, Geneva, CH.
Illingworth et al, “Vision, Image and Signal Processing,” The Institution of Electrical Engineers, IEE Proceedings, Jun. 2000, vol. 147, No. 3.
McLean et al, “Telecine Noise Reduction,” XP-002285972, 2001 The Institute of Electrical Engineers.
Oktem et al, “Transform Domain Algorithm for Reducing Effect of Film-Grain Noise in Image Compression,” Electronics Letters, Oct. 14, 1999, vol. 35, No. 21.
Peng et al, “Adaptive Frequency Weighting for Fine-Granularity-Scalability,” Visual Communications and Image Processing 2002, Proceedings of SPIE, vol. 4671, 2002 SPIE 0277-786X/02.
Prades-Nebot et al, “Rate Control for Fully Fine-Grained Scalable Video Coders,” Visual Communications and Image Processing 2002, Proceedings of SPIE, vol. 4671 (2002), SPIE 0277-786X/02.
Schaar et al, “Fine-Granularity-Scalability for Wireless Video and Scalable Storage,” Visual Communications and Image Processing 2002, Proceedings of SPIE, vol. 4671 (2002) SPIE 0277-786X/02.
Shahnaz et al, “Image Compression in Signal-Dependent Noise,” Applied Optics, vol. 38, No. 26, Sep. 10, 1999.
Yan et al, “Efficient Video Coding with Hybrid Spatial and Fine-Grain SNR Scalabilities,” Department of Electronic Engineering, Beijing Institute of Technology, China.
Yan et al, “Film Grain Noise Removal and Generation for Color Images,” Department Electrical and Computer Engineers, University of Toronto, dimitris@comm.toronto.edu.
Yan et al, “Signal-Dependent Film Grain Noise Removal and Generation Based on Higher-Order Statistics,” University of Toronto, Department of Electrical and Computer Engineering, dimitris@comm.toronto.edu, 1997 IEEE.
Yoshida, “Go with the Grain, Film R&D Chief Urges, for Art's Sake,” EE Times, Feb. 7, 2005.
Zhang et al, “Constant Quality Constrained Rate Allocation for FGS Video Coded Bitstreams,” Visual Communications and Image Processing 2002, Proceedings of SPIE, vol. 4671 (2002) SPIE 0277-786X/02.
International Search Report Dated Feb. 22, 2006.
Office Action for U.S. Appl. No. 10/552,179 mailed Sep. 1, 2010.
Office Action from U.S. Appl. No. 10/556,833 mailed May 10, 2010.
Final Office Action from U.S. Appl. No. 10/556,833 mailed Oct. 20, 2010.
Office Action ADV ACT from U.S. Appl. No. 10/556,833 mailed Jan. 7, 2011.
Office Action from U.S. Appl. No. 10/556,833 mailed Feb. 15, 2011.
Office Action from U.S. Appl. No. 10/556,833 mailed May 10, 2011.
Office Action from U.S. Appl. No. 10/556,833 mailed Sep. 30, 2011.
Final Office Action from U.S. Appl. No. 10/571,148 mailed May 12, 2010.
Office Action from U.S. Appl. No. 10/571,148 mailed Sep. 10, 2010.
Office Action ADV ACTION from U.S. Appl. No. 10/572,820 mailed Mar. 19, 2010.
Office Action from U.S. Appl. No. 10/572,820 mailed May 11, 2010.
Final Office Action from U.S. Appl. No. 10/572,820 mailed Sep. 30, 2010.
Office Action ADV ACTION from U.S. Appl. No. 10/575,676 mailed Jun. 7, 2010.
Office Action from U.S. Appl. No. 10/575,676 mailed Jul. 19, 2011.
Finala Office Action from U.S. Appl. No. 10/575,676 mailed Oct. 26, 2011.
Office Action from U.S. Appl. No. 10/581,151 mailed Jan. 20, 2010.
Office Action from U.S. Appl. No. 11/246,474 mailed Sep. 1, 2010.
Office Action from U.S. Appl. No. 11/246,848 mailed Jul. 6, 2010.
Final Office Action from U.S. Appl. No. 11/246,848 mailed Oct. 13, 2010.
Office Action from U.S. Appl. No. 11/268,070 mailed Jul. 21, 2010.
Final Office Action from U.S. Appl. No. 11/268,070 mailed Nov. 12, 2010.
Office Action ADV ACTION from U.S. Appl. No. 11/268,070 mailed Jan. 11, 2011.
Office Action from U.S. Appl. No. 11/268,070 mailed Jun. 22, 2011.
Office Action from U.S. Appl. No. 11/273,067 mailed Aug. 4, 2010.
Office Action from U.S. Appl. No. 11/667,581 mailed Jun. 8, 2011.
Final Office Action from U.S. Appl. No. 11/667,581 mailed Oct. 14, 2011.
Office Action from U.S. Appl. No. 11/667,816 mailed Sep. 15, 2011.
Final Office Action from U.S. Appl. No. 11/667,816 mailed Dec. 5, 2011.
Office Action from U.S. Appl. No. 11/667,846 mailed Oct. 4, 2011.
Office Action from U.S. Appl. No. 11/284,378 mailed Dec. 22, 2010.
Final Office Action from U.S. Appl. No. 11/285,540 mailed Nov. 23, 2010.
Office Action from U.S. Appl. No. 11/285,540 mailed Mar. 30, 2011.
Office Action from U.S. Appl. No. 12/589,217 mailed Jul. 13, 2010.
Final Office Action from U.S. Appl. No. 12/589,217 mailed Oct. 28, 2010.
Office Action ADV ACTION from U.S. Appl. No. 12/589,217 mailed Dec. 2, 2010.
Byun et al: Power Efficient MPEG-4 Decoder Featuring Low-Complexity Error Resilience, ASIC, 2002 Proceedings IEEE Asia Pacific Conference Aug. 6-8, 2002 Piscataway, NJ USA IEEE.
Conklin et al: “Dithering 5-Tap Filter for Inloop Deblocking,” JVT of ISO/IEC MPEG & ITU-T VCEG, 3rd Meeting, Fairfax, VA May 6-10, 2002, pp. 1016.
Gomila et al: “Film Grain Modeling vs. Encoding”, JVT of ISO/IEC MPEG & ITU-T VCEG 11th Meeting, Munich, DE, Mar. 15-19, 2004.
Naderi et al: “Estimation of Images Degraded by Film Grain Noise”, Applied Optics, vol. 17, Issue 8, pp. 1223-1237, Jan. 1, 1978.
Pirsch et al: “VLSI Architures for Video Compression—A Survey”, Proceedings of IEEE, New York, USA, vol. 83, No. 2, Feb. 1, 1995, pp. 220-246.
Puri et al: “Video Coding Using the H.264/MPEG-4 AVC Compression Standard”, Signal Processing Image Communication, Elsevier Science Publishers, Amsterdam, NL, vol. 19, No. 9. Oct. 1, 2004, pp. 7993-7894.
Takashi et al: “a 60-Mw MPEG-4 Video Codec Using Clustered Voltage Scaling w/Variable Supply Voltage Scheme”, IEEE Journal of Solid State Circuits, Piscataway, NJ, Nov. 1998,v.33 #11.
Naderi et al., “Estimation of Images Degraded by Film-Grain Noise”, 1978 Optical Society of America, 1978, pp. 1228-1237.
Byun et al., “Power Efficient MPEG-4 Decoder Architecture Featuring Low-Complexity Error Resilience”, 2000 IEEE.
Pirsch et al., “VLSI Architectures for Video Compression—A Survey”, Proceedings ofthe IEEE, vol. 83, No. 2, Feb. 1995.
Conklin et al., “Dithering 5-Tap Filter for Inloop Deblocking”, JVT of ISO/IEC PEG & ITU-T VCEG, Document: JVT-C056, 3rd eeting: Fairfax, VA, May 6-10, 2002.
Bjontegaard—“Addition—of—comfort—noise—as—post—processing”: ITU Telecommunications Standarization Sector, Sunriver, Oregon, Sep. 8, 1997. pp. 102.
Lerner—“Fixed—vs.—floating—point: A—surprisingly—hard—choice”. eetimes, Feb. 6, 2007. pp. 1-4.
McMahon—EtAl—“High—Quality—SD—and—HD—AVC—Test—Results”. Joint Video Team (JVT of ISO/IEC MPEG & ITU-T VCEG) Geneva, Switzerland. Oct. 9, 2002.
Schlockerman—EtAl: “Film Grain—coding—in—H.264/AVC”; JVT (JVT of ISO/IEC MPEG & ITU-T VCEG); 9th Meeting, San Diego, California Sep. 2, 2003.
Wiegand—EtAl—“Overview—of—the—H.264/AVC—Video—Coding—Standard” IEEE Transactions on Circuits and Systems for Video Technology, vol. 13, No. 7. Jul. 2003.
Sullivan—EtAl—“The H.264/AVC advanced—video—coding—standard:—Overview—and—Introduction—to—the—Fidelity—Range—Extensions”. Proceedings of SPIE 5558 App of Digital Image Processing XXVII, Nov. 2, 2004.
Schlockerman—EtAl: “Film—Grain—coding—in—H.264/AVC”; JVT (JVT on of ISO/IEC MPEG & ITU-T VCEG); 9th Meeting, San Diego, California Sep. 2, 2003.
Sullivan—EtAl—“The—H.264/AVC—advanced—video—coding—standard:—Overview—and—Introduction—to—the—Fidelity—Range—Extensions”. Proceedings of SPIE 5558 App of Digital Image Processing XXVII, Nov. 2, 2004.
Related Publications (1)
Number Date Country
20070269125 A1 Nov 2007 US
Provisional Applications (1)
Number Date Country
60628837 Nov 2004 US