# Patent

## Title
Method and device for processing a sequence of video images

## Abstract

The invention discloses a new modulation scheme to prevent illegal copy in movie theaters. None of the existing modulation schemes make the right assumptions to modelize a camcorder device. In such classical schemes, camcorders are always considered as sampling devices without taking into account the effects of shutter speed. The shutter of the camcorder is classically compared to a low-pass filter. As a result, modulation effects can be easily removed by setting up shutter speed to low values (low cutoff frequency). The right mathematical modelization of the shutter effect shows that, in the spectral domain, the behavior of the camcorder is close to a cardinal sine function with a main lobe and sidelobes. Consequently, according to the invention, it is proposed to select a modulation frequency that not only generates visual artifacts once recorded by the camcorder under shutter-free conditions, but that can also generate artifacts despite low shutter speed settings by going through sidelobe of the shutter spectrum.

## Background
This application claims the benefit, under 35 U.S.C. §119 of EP Patent Application 06290155.8, filed Jan. 26, 2006.
FIELD OF THE INVENTION
The invention relates to a method for processing a sequence of video images represented by their video signals comprising a modulation step for modulating temporally the amplitude of the brightness of a set of pixels of each image of the sequence around a brightness value to be displayed for said image, this modulation being not visible to the human eye.
BACKGROUND OF THE INVENTION
The visual contents, whether these be fixed or moving images, are in general creations that benefit from guarantees of exclusivity associated with the creator's rights. Their reproduction is in general permitted only within a strictly defined framework that allows the creators and their beneficiaries to be remunerated.
To ensure that these legal rules are complied with correctly, many systems have been developed to prevent illegal copies or to make the quality of the copies sufficiently degraded to make them unusable.
Within this context, the patent application EP 1 237 369 aims to combat the copying of images by means of a camera while they are being displayed, for example using a camcorder in a movie theatre. In this document, it has been proposed to modulate temporally the amplitude of the brightness of selected pixels representing an anti-piracy message around the value to be displayed at a high rate that makes the message invisible to the human eye but generates artifacts in the sequence filmed by the camcorder. Such a solution requires a modulation at a rate higher than the flicker fusion frequency, which is of around 50 Hz, and therefore applies only to systems having a high image refresh rate, at least of around 100 Hz. Applied to systems with a lower display rate (50 Hz or 60 Hz for example), the modulation could be visible to the human eye and would visibly degrade the rendition of the displayed image.
The major drawback of this method is that the anti-piracy message can be easily removed by setting up the camcorder to low shutter speed settings.
SUMMARY OF THE INVENTION
It is an object of the present invention to elaborate a method which is both able to generate aliasing artifacts on camcorder copies of a digital motion picture, and resist to low shutter speed settings.
In the present invention, the mathematical modelization of the shutter effect shows that, in the spectral domain, the behavior of the camcorder is close to a cardinal sine function with a main lobe and sidelobes. Consequently, according to the invention, it is proposed to select a modulation frequency that not only generates visual artifacts once recorded by the camcorder under shutter-free conditions, but that can also generate artifacts despite low shutter speed settings by going through sidelobes of the shutter spectrum.
Thus, the object of the invention is solved by a method for processing a sequence of video images represented by their video signals comprising a modulation step for modulating temporally the amplitude of the brightness of a set of pixels of each image of the sequence around a brightness value to be displayed for said image, wherein said modulation is chosen to be not visible to the human eye while simultaneously generating artifacts due to aliasing when said video images are captured by a digital video capturing device. According to the invention, the modulation frequency is selected to be in a limited frequency range centered on a local extremum of the function sinc(πfτ), where the function sinc is the cardinal sine function and τ is the shutter integration time of the digital video capturing device.
In a standard case, the shutter integration time τ is adjusted to the inverse of the sampling frequency fs of the digital video capturing device which is for example a PAL camcorder or a NTSC camcorder.
According to the invention, in that standard case, to defeat PAL camcorders (sampling frequency equal to 50 Hz), the modulation frequency is thus selected to be in a limited frequency range centered on the frequency 71.5 Hz or 122.9 Hz. For a NTSC camcorders (sampling frequency equal to 60 Hz), the modulation frequency is selected to be in a limited frequency range centered on the frequency 85.8 Hz or 147.5 Hz.
The invention concerns also a device for processing a sequence of video images represented by their video signals comprising a modulator for modulating temporally the amplitude of the brightness of a set of pixels of each image of the sequence around a brightness value to be displayed for said image, wherein the modulation is chosen to be not visible to the human eye while simultaneously generating artifacts due to aliasing when said video images are captured by a digital video capturing device. According to the invention, the modulation frequency is selected to be in a limited frequency range centered on a local extremum of the function sinc(πfτ), where the function sinc is the cardinal sine function and τ is the shutter integration time of the digital video capturing device.
As mentioned before, when the shutter integration time τ is adjusted to the inverse of the sampling frequency fs of the digital video capturing device, the modulation frequency used by the modulator is selected to be in a limited frequency range centered on the frequency 71.5 Hz or 122.9 Hz to defeat PAL camcorders and selected to be in a limited frequency range centered on the frequency to 85.8 Hz or 147.5 Hz to defeat NTSC camcorders.
Furthermore, the images computed by the modulator are outputted at a frame rate that is equal to or higher than twice the modulation frequency Preferably, this frame rate is equal to a multiple of the frame rate at which the images are inputted into the modulator.



BRIEF DESCRIPTION OF THE DRAWINGS
Exemplary embodiments of the invention are illustrated in the drawings and are explained in more detail in the following description. In the drawings:
FIG. 1 is a diagram illustrating the flicker sensitivity of the human eye;
FIG. 2 is a couple of diagrams illustrating the uniform sampling of a function x(t);
FIG. 3 is a couple of diagrams illustrating the sampling process inn the time domain and in the frequency domain;
FIG. 4 is a couple of diagrams illustrating a sampling operation with a low sampling frequency;
FIG. 5 is a couple of diagrams showing the frequency bandwidth of a base band signal and an amplitude-modulated signal;
FIG. 6 illustrates the amplitude of the function sinc(πfT) for the frequency band 0-50 Hz;
FIG. 7 illustrates the amplitude of the function sinc(πfT) for the full spectrum;
FIG. 8 represents four diagrams illustrating the method of the invention;
FIG. 9 shows an example of warning message; and
FIG. 10 shows a circuit implementation of the present invention.



DESCRIPTION OF PREFERRED EMBODIMENTS
The invention specifies a new modulation system to prevent illegal copy in movie theaters. Although other modulation schemes have already been proposed to defeat camcorders, none of those make the right assumptions to modelize a camcorder device. In such classical schemes, camcorders are always considered as sampling devices only while a good model should take into account the effects of shutter speed. The shutter of the camcorder is classically compared to a variable cutoff frequency low-pass filter located before the sampling process in a functional scheme. As a result, modulation effects can be easily removed by setting up shutter speed to low values (low cutoff frequency). In the present invention, the mathematical modelization of the shutter effect shows that, in the spectral domain, the behavior of the camcorder is close to a cardinal sine function with a main lobe and sidelobes. Consequently, according to the invention, it is proposed to select a modulation frequency that not only generates visual artifacts once recorded by the camcorder under shutter-free conditions, but that can also generate artifacts despite low shutter speed settings by going through sidelobes of the shutter spectrum.
Before explaining in detail the invention, some information about the human visual system, the Fourier transformation, the sampling process and the amplitude modulation are given.
According to several psychophysical studies, the Human Visual System or HVS is selectively sensitive to different ranges of both spatial and temporal frequencies. These studies usually consist in analyzing the reactions of a human observer to cyclic (sinusoidal) variations of brightness in either spatial and/or temporal domains. The term “spatial frequency” refers to the number of cycles per degree of visual angle, whereas the term “temporal frequency” refers to the number of cycles per second. In the present application, the term “frequency” will refer to temporal frequency.
The brightness flicker sensitivity of the HSV according to frequency is illustrated by FIG. 1. The brightness flicker becomes imperceptible for a frequency equal to or greater than 50 Hz.
Fourier transformation is a well-known tool to analyze a time-varying function in the frequency (or “spectral”) domain. Let x(t) be a continuous function. Its corresponding Fourier transform can be written as:




  
    
      
        
          F
          ⁡
          
            [
            
              x
              ⁡
              
                (
                t
                )
              
            
            ]
          
        
        =
        
          
            X
            ⁡
            
              (
              f
              )
            
          
          =
          
            
              ∫
              
                -
                ∞
              
              
                +
                ∞
              
            
            ⁢
            
              
                x
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                e
                
                  
                    -
                    j2π
                  
                  ⁢
                  
                    
                  
                  ⁢
                  ft
                
              
              ⁢
              
                d
                t
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




The corresponding inverse Fourier transform is:




  
    
      
        
          
            F
            
              -
              1
            
          
          ⁡
          
            [
            
              X
              ⁡
              
                (
                f
                )
              
            
            ]
          
        
        =
        
          
            x
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              ∫
              
                -
                ∞
              
              
                +
                ∞
              
            
            ⁢
            
              
                X
                ⁡
                
                  (
                  f
                  )
                
              
              ⁢
              
                e
                
                  j2π
                  ⁢
                  
                    
                  
                  ⁢
                  ft
                
              
              ⁢
              
                d
                f
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




A Fourier transform is defined by an amplitude and a phase angle although amplitude is the only interesting data to focus on. Examples of conversions from temporal to frequency domains of well-known functions and operators are given in the following table 1:















x(t) = F−1[X(f)]
X(F) = F[x(t)]



1
δ(f)


 







  e
  
    
      -
      
        
          j2π
          ⁢
          f
        
        0
      
    
    ⁢
    t
  




δ(f + f0)


 



a · x(t) + b · y(t)
a · X(f) + b · Y(f)



x(t) · y(t)
X(F) * Y(F)



x(t) * y(t)
X(F) · Y(f)


 







  
    e
    
      
        -
        
          
            j2π
            ⁢
            f
          
          0
        
      
      ⁢
      t
    
  
  ⁣
  
    ·
    
      x
      ⁡
      
        (
        t
        )
      
    
  




δ(f + f0) * X(f) = X(f + f0)


 







  
    cos
    ⁡
    
      (
      
        2
        ⁢
        
          πf
          0
        
        ⁢
        t
      
      )
    
  
  =
  
    
      
        e
        
          
            -
            
              
                j2π
                ⁢
                f
              
              0
            
          
          ⁢
          t
        
      
      +
      
        e
        
          
            
              j2π
              ⁢
              f
            
            0
          
          ⁢
          t
        
      
    
    2
  








  
    1
    2
  
  ⁢
  
    (
    
      
        δ
        ⁡
        
          (
          
            f
            +
            
              f
              0
            
          
          )
        
      
      +
      
        δ
        ⁡
        
          (
          
            f
            -
            
              f
              0
            
          
          )
        
      
    
    )
  






 



δ(t − t0)




  e
  
    -
    
      
        j2π
        ⁢
        f
        ⁢
        t
      
      0
    
  






 







  
    ∑
    
      n
      =
      
        -
        ∞
      
    
    
      n
      =
      
        +
        ∞
      
    
  
  ⁢
  
    δ
    ⁡
    
      (
      
        t
        -
        
          nT
          S
        
      
      )
    
  








  
    ∑
    
      n
      =
      
        -
        ∞
      
    
    
      n
      =
      
        +
        ∞
      
    
  
  ⁢
  
    
      1
      
        T
        S
      
    
    ⁢
    
      δ
      ⁡
      
        (
        
          f
          -
          
            n
            
              T
              S
            
          
        
        )
      
    
  













Fourier transformation turns out to be very useful when it comes to analyzing Analog-to-Digital and Digital-to-Analog Conversion schemes (ADC/DAC). The conversion of any function from analog to digital always starts with a time sampling process, which consists in switching from continuous to discrete representation spaces. This discretization process can be mathematically modelized by a multiplication between the function and an infinite series of Dirac's delta functions separated by Ts, where Ts is the sampling period and




  
    f
    S
  
  =
  
    1
    
      T
      S
    
  




is the sampling frequency:





  
    
      
        
          {
          
            x
            n
          
          }
        
        =
        
          
            
              ∑
              
                n
                =
                
                  -
                  ∞
                
              
              
                n
                =
                
                  +
                  ∞
                
              
            
            ⁢
            
              x
              ⁡
              
                (
                
                  nT
                  S
                
                )
              
            
          
          =
          
            
              x
              ⁡
              
                (
                t
                )
              
            
            ·
            
              
                ∑
                
                  n
                  =
                  
                    -
                    ∞
                  
                
                
                  n
                  =
                  
                    +
                    ∞
                  
                
              
              ⁢
              
                δ
                ⁡
                
                  (
                  
                    t
                    -
                    
                      nT
                      S
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




The sampling process of the function x(t) is illustrated by FIG. 2. Using transforms presented in Table 1, the Fourier transform of discrete function {xn} can be written as follows:




  
    
      
        
          F
          ⁡
          
            [
            
              {
              
                x
                n
              
              }
            
            ]
          
        
        =
        
          F
          ⁡
          
            [
            
              
                x
                ⁡
                
                  (
                  t
                  )
                
              
              ·
              
                
                  ∑
                  
                    n
                    =
                    
                      -
                      ∞
                    
                  
                  
                    n
                    =
                    
                      +
                      ∞
                    
                  
                
                ⁢
                
                  δ
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        nT
                        S
                      
                    
                    )
                  
                
              
            
            ]
          
        
      
    
  
  
    
      
        =
        
          
            F
            ⁡
            
              [
              
                x
                ⁡
                
                  (
                  t
                  )
                
              
              ]
            
          
          *
          
            f
            ⁡
            
              [
              
                
                  ∑
                  
                    n
                    =
                    
                      -
                      ∞
                    
                  
                  
                    n
                    =
                    
                      +
                      ∞
                    
                  
                
                ⁢
                
                  δ
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        nT
                        S
                      
                    
                    )
                  
                
              
              ]
            
          
        
      
    
  
  
    
      
        =
        
          
            X
            ⁡
            
              (
              f
              )
            
          
          *
          
            
              ∑
              
                n
                =
                
                  -
                  ∞
                
              
              
                n
                =
                
                  +
                  ∞
                
              
            
            ⁢
            
              
                1
                
                  T
                  S
                
              
              ⁢
              
                δ
                ⁡
                
                  (
                  
                    f
                    -
                    
                      n
                      
                        T
                        S
                      
                    
                  
                  )
                
              
            
          
        
      
    
  
  
    
      
        =
        
          
            f
            S
          
          ⁢
          
            
              ∑
              
                n
                =
                
                  -
                  ∞
                
              
              
                n
                =
                
                  +
                  ∞
                
              
            
            ⁢
            
              X
              ⁡
              
                (
                
                  f
                  -
                  
                    nf
                    S
                  
                
                )
              
            
          
        
      
    
  




In other words, sampling a function in the space/time domain introduces a periodization in the frequency domain as illustrated by FIG. 3. In this figure, the spectrum of the function x(t) is spread over the frequency band [0-F]. As shown by this figure, x(t) can be fully recovered from {xn} because F[x(t)] can be recovered from a low-pass filtering operation on F[{xn}]. However, in order to perform lossless Analog-to-Digital/Digital-to-Analog conversions, special care must be taken regarding chosen sampling times/frequencies. If the sampling frequency fs is lower than 2F, the periods of F[x(t)] tend to overlap, thus leading to the spectrum F[{xn}] illustrated by FIG. 4. This new spectrum does not allow to recover the original spectrum F[x(t)]. The overlapping phenomenon, called aliasing, can be prevented if the well-known Nyquist-Shannon sampling theorem is respected. This theorem states that when sampling a signal (e.g., converting from an analog signal to digital), the sampling frequency must be greater than twice the bandwidth of the input signal in order to be able to reconstruct the original perfectly from the sampled version: fs>2F.
Therefore, as any sampling device, a camcorder can generate aliasing artifacts. In the case of an anti-piracy or anti-camcorder system, it can be useful to make the camcorder generating such aliasing artifacts. Since it is not possible to act on the sampling frequency fs of the camcorder (50 Hz PAL/60 Hz NTSC for interleaved fields), the principle of the invention is to increase the bandwidth F of the video signal by using amplitude modulation with specific modulation frequencies that are not filtered by the camcorder.
Amplitude modulation (AM) is a form of modulation in which the amplitude of a carrier signal changes depending on the amplitude of a modulating signal. A basic AM operation consists in multiplying the modulating signal, for example x(t), with the carrier signal c(t) of frequency fm. The carrier signal c(t) is defined as:

c(t)=β cos(2πfmt)

The amplitude-modulated signal xAM(t) is:

xAM(t)=x(t)·(1+c(t))=x(t)+x(t)·β cos(2πfmt)

The Fourier transform of the signal xAM(t) is:




  
    
      
        
          
            X
            AM
          
          ⁡
          
            (
            f
            )
          
        
        =
        
          
            X
            ⁡
            
              (
              f
              )
            
          
          +
          
            
              X
              ⁡
              
                (
                f
                )
              
            
            *
            
              1
              2
            
            ⁢
            
              (
              
                
                  δ
                  ⁡
                  
                    (
                    
                      f
                      -
                      
                        f
                        m
                      
                    
                    )
                  
                
                +
                
                  δ
                  ⁡
                  
                    (
                    
                      f
                      +
                      
                        f
                        m
                      
                    
                    )
                  
                
              
              )
            
          
        
      
    
  
  
    
      
        =
        
          
            X
            ⁡
            
              (
              f
              )
            
          
          +
          
            
              1
              2
            
            ⁢
            
              (
              
                
                  X
                  ⁡
                  
                    (
                    
                      f
                      -
                      
                        f
                        m
                      
                    
                    )
                  
                
                +
                
                  X
                  ⁡
                  
                    (
                    
                      f
                      +
                      
                        f
                        m
                      
                    
                    )
                  
                
              
              )
            
          
        
      
    
  




A graphical representation of this Fourier transform is given at FIG. 5. It shows the amplitude spectrum X(f) of the modulating signal and the full amplitude spectrum XAM(f) of the amplitude-modulated signal.
In reference to FIG. 5, it can be easily understood how useful this amplitude modulation can be to generate aliasing artifacts on a recorded video sequence. If the modulation frequency fm is high enough to increase the bandwidth of the video signal over the limit fs/2 and high enough to make extra spectral content invisible to the human visual system, then this extra content will overlap during camcorder acquisition, creating artifacts a human observer is able to perceive. In order to create a warning message, modulation shall be applied to a selected set of pixels representing a warning message as disclosed in EP 1 237 369.
Until now, modulation-based anti-camcorder methods only considered camcorders as sampling devices, which is a wrong assumption as it will be seen later. In camcorder technology, the recording/sampling process of a visual signal first requires its exposure to CCD/CMOS sensors for a certain amount of time. Exposure is commonly denoted as “shutter speed”, although the shutter used in most video cameras is not mechanical and although it is a time period and not a speed. Chip camera “shutter” speed simply represents the inverse of the time during which the light-induced charge is allowed to electronically build in the chip before the cycle is repeated. In the frequency domain, shutter is commonly compared to a low pass filter, with exposure time being directly related to the cutoff frequency of the filter.
For a PAL camcorder device recording 25 frames per second, the actual sampling frequency fs is 50 Hz because of the use of interlaced frames. Based on the Nyquist-Shannon sampling theorem, to generate aliasing artifacts, the modulation frequency fm should be set to any value higher than fs/2, which leads us to the condition:

fm>25 Hz

But, to be sure the modulation effects will not be perceived by the legal audience, the modulation frequency fm should not be located in frequency bands where a human observer can perceive flicker (see FIG. 1), hence a second condition is:

fm>50 Hz

This new condition is a more constraining one. It is a well-known fact that cinema standards set refreshing rates to 24 images per second. In case of an attempt of piracy over a motion picture, setting up the shutter speed of the camcorder to 1/50, which is a standard (value for PAL camcorders), does not modify the motion since motion will never exceed 12 Hz. Furthermore, if a shutter speed of 1/50 equals a low pass filter with cutoff frequency of 50 Hz and, as the filtering process occurs before the sampling process in the camcorder, it would never be possible to generate a modulation effect that goes through the shuttering process with the condition fm>50 Hz. Then aliasing artifacts generation should be impossible.
But below is given a mathematical modelization of the shutter behavior showing that the shutter is not just a low pass filter and that an appropriate amplitude modulation can generate aliasing.
Exposure time is the interval during which the CCD/CMOS sensors are exposed to incident light by the shutter, while integration time is defined as the interval during which the clocks of the camcorder are set to trap and retain charge. For a CCD, the integration starts when the CCD is cleared. It is counted from the end of clearing, until the CCD starts to read out. Although integration time is delimited by the behaviour of the readout electronics of the camcorder, it can be considered that it equals exposure time. So the shutter effect can be modelized by the following equation:




  
    
      
        
          
            x
            shut
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            ∫
            
              t
              -
              T
            
            t
          
          ⁢
          
            
              x
              ⁡
              
                (
                τ
                )
              
            
            ⁢
            
              d
              τ
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




    
    
        where
        
            T is the shutter speed,
            x(t) is the input signal, and
            xshut(t) is the output signal.
        
        
    
    


Once integrated, the incoming signal x(t) enters the sampling process of the camcorder. The actual model of a recorded sequence xrec(t) can then be written as:




  
    
      
        
          
            x
            rec
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            ∑
            
              n
              =
              
                -
                ∞
              
            
            
              +
              ∞
            
          
          ⁢
          
            
              [
              
                
                  ∫
                  
                    t
                    -
                    T
                  
                  t
                
                ⁢
                
                  
                    x
                    ⁡
                    
                      (
                      τ
                      )
                    
                  
                  ⁢
                  
                    d
                    τ
                  
                
              
              ]
            
            ·
            
              
                ∑
                
                  n
                  =
                  
                    -
                    ∞
                  
                
                
                  n
                  =
                  
                    +
                    ∞
                  
                
              
              ⁢
              
                δ
                ⁡
                
                  (
                  
                    t
                    -
                    
                      nT
                      S
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




Using two changes of variable in Equation 3, it turns out xshut(t) can be written as:




  
    
      
        
          
            
              
                
                  x
                  shut
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
              ⁢
              
                
                  ∫
                  
                    t
                    -
                    τ
                  
                  t
                
                ⁢
                
                  
                    x
                    ⁡
                    
                      (
                      τ
                      )
                    
                  
                  ⁢
                  
                    d
                    τ
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    ∫
                    
                      -
                      T
                    
                    0
                  
                  ⁢
                  
                    
                      x
                      ⁡
                      
                        (
                        
                          φ
                          +
                          t
                        
                        )
                      
                    
                    ⁢
                    
                      d
                      φ
                    
                    ⁢
                    
                      
                    
                    ⁢
                    with
                    ⁢
                    
                      
                    
                    ⁢
                    φ
                  
                
                =
                
                  τ
                  -
                  t
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  ∫
                  
                    -
                    ∞
                  
                  
                    +
                    ∞
                  
                
                ⁢
                
                  
                    
                      x
                      ⁡
                      
                        (
                        
                          φ
                          +
                          t
                        
                        )
                      
                    
                    ·
                    
                      
                        rect
                        T
                      
                      ⁡
                      
                        (
                        
                          φ
                          +
                          
                            T
                            2
                          
                        
                        )
                      
                    
                  
                  ⁢
                  
                    d
                    φ
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    ∫
                    
                      -
                      ∞
                    
                    
                      +
                      ∞
                    
                  
                  ⁢
                  
                    
                      
                        x
                        ⁡
                        
                          (
                          χ
                          )
                        
                      
                      ·
                      
                        
                          rect
                          T
                        
                        ⁡
                        
                          (
                          
                            χ
                            -
                            t
                            +
                            
                              T
                              2
                            
                          
                          )
                        
                      
                    
                    ⁢
                    
                      d
                      χ
                    
                    ⁢
                    
                      
                    
                    ⁢
                    with
                    ⁢
                    
                      
                    
                    ⁢
                    χ
                  
                
                =
                
                  φ
                  +
                  t
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  ∫
                  
                    -
                    ∞
                  
                  
                    +
                    ∞
                  
                
                ⁢
                
                  
                    
                      x
                      ⁡
                      
                        (
                        χ
                        )
                      
                    
                    ·
                    re
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      ct
                      T
                    
                    ⁡
                    
                      (
                      
                        
                          -
                          χ
                        
                        +
                        t
                        -
                        
                          T
                          2
                        
                      
                      )
                    
                  
                  ⁢
                  
                    d
                    χ
                  
                
              
            
          
        
        
          
            
              
              ⁢
              
                (
                
                  
                    
                      rect
                      T
                    
                    ⁡
                    
                      (
                      .
                      )
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  is
                  ⁢
                  
                    
                  
                  ⁢
                  an
                  ⁢
                  
                    
                  
                  ⁢
                  even
                  ⁢
                  
                    
                  
                  ⁢
                  function
                
                )
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  [
                  
                    x
                    *
                    
                      rect
                      T
                    
                  
                  ]
                
                ⁢
                
                  (
                  
                    t
                    -
                    
                      T
                      2
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




where rectT(.) is the rectangular function, defined as:





  
    
      rect
      T
    
    ⁡
    
      (
      t
      )
    
  
  =
  
    {
    
      
        
          1
        
        
          if
        
        
          
            t
            ∈
            
              [
              
                
                  
                    -
                    T
                  
                  2
                
                ,
                
                  T
                  2
                
              
              ]
            
          
        
      
      
        
          0
        
        
          else
        
        
          
            
          
        
      
    
  




This result is very interesting, especially when looking at the corresponding Fourier transform of function xshut(t):




  
    
      
        
          
            
              
                F
                ⁡
                
                  [
                  
                    
                      x
                      shut
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                  ]
                
              
              =
              
                F
                ⁡
                
                  [
                  
                    
                      [
                      
                        x
                        *
                        
                          rect
                          T
                        
                      
                      ]
                    
                    ⁢
                    
                      (
                      
                        t
                        -
                        
                          T
                          2
                        
                      
                      )
                    
                  
                  ]
                
              
            
          
        
        
          
            
              =
              
                
                  F
                  ⁡
                  
                    [
                    
                      
                        [
                        
                          x
                          *
                          
                            rect
                            T
                          
                        
                        ]
                      
                      ⁢
                      
                        (
                        t
                        )
                      
                    
                    ]
                  
                
                ·
                
                  e
                  
                    
                      -
                      jπ
                    
                    ⁢
                    
                      
                    
                    ⁢
                    fT
                  
                
              
            
          
        
        
          
            
              =
              
                
                  F
                  ⁡
                  
                    [
                    
                      x
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                    ]
                  
                
                ·
                
                  F
                  ⁡
                  
                    [
                    
                      
                        rect
                        T
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                    ]
                  
                
                ·
                
                  e
                  
                    
                      -
                      jπ
                    
                    ⁢
                    
                      
                    
                    ⁢
                    fT
                  
                
              
            
          
        
        
          
            
              =
              
                
                  
                    X
                    ⁡
                    
                      (
                      f
                      )
                    
                  
                  .
                  T
                
                ⁢
                
                  
                
                ⁢
                sin
                ⁢
                
                  
                
                ⁢
                
                  
                    c
                    ⁡
                    
                      (
                      
                        π
                        ⁢
                        
                          
                        
                        ⁢
                        fT
                      
                      )
                    
                  
                  ·
                  
                    e
                    
                      
                        -
                        jπ
                      
                      ⁢
                      
                        
                      
                      ⁢
                      fT
                    
                  
                
              
            
          
        
        
          
            
              =
              
                
                  X
                  shut
                
                ⁡
                
                  (
                  f
                  )
                
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




While the term e−jπfT in the equation (6) has only an influence on the phase of Xshut(f), it turns out that the term sinc(πfT) defines the actual mathematical model for the low pass filtering operation that was induced by shutter adjustment. FIG. 6 is a graphical representation of the cardinal sine function sinc(πfT) for the frequency band 0-50 Hz. As it can be seen, when set to a speed of 1/50, the shutter behaves as a low pass filter with a cutoff frequency of 50 Hz or even lower if a −3 dB cutoff frequency is considered. Because of that, all attempts to generate modulation artifacts will either be visible to the legal audience (for 25 Hz<fm<50 Hz—see FIG. 1) or removed by the shutter (for fm>50 Hz).
Using a mathematical model to study the effect of the shutter is not useless though. FIG. 7 shows the graphical representation of function sinc(πfT) over the full frequency band, which turns out to be very interesting: although cutoff is detected at 50 Hz, a function sinc(.) presents severe sidelobe effects. This means it is still possible to let certain high frequencies go through the shutter. The function sinc(.) is known to present zero-crossing for points that are multiples of π (i.e. frequencies that are multiples of




  1
  T




for sinc(πfT)). In the present application, points which are worth focusing on are the ones that indicate the location of sidelobe extrema. Local extrema are found when a function has a derivative of zero:





  
    
      
        
          
            
              d
              
                (
                
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  
                    c
                    ⁡
                    
                      (
                      x
                      )
                    
                  
                
                )
              
            
            
              d
              x
            
          
          =
          0
        
        ⁢
        
          

        
        ⁢
        
          
            
              
                x
                ·
                
                  cos
                  ⁡
                  
                    (
                    x
                    )
                  
                
              
              -
              
                sin
                ⁡
                
                  (
                  x
                  )
                
              
            
            
              x
              2
            
          
          =
          0
        
        ⁢
        
          

        
        ⁢
        
          
            
              x
              ·
              
                cos
                ⁡
                
                  (
                  x
                  )
                
              
            
            -
            
              sin
              ⁡
              
                (
                x
                )
              
            
          
          =
          0
        
      
    
    
      
        (
        7
        )
      
    
  




It is not easy to find an analytic solution to Equation 7. But the function sinc(x) is polynomial:




  
    
      
        
          sin
          ⁢
          
            
          
          ⁢
          
            c
            ⁡
            
              (
              x
              )
            
          
        
        =
        
          
            ∑
            
              n
              =
              0
            
            
              +
              ∞
            
          
          ⁢
          
            
              
                
                  (
                  
                    -
                    1
                  
                  )
                
                n
              
              
                
                  (
                  
                    
                      2
                      ⁢
                      n
                    
                    +
                    1
                  
                  )
                
                !
              
            
            ⁢
            
              x
              
                2
                ⁢
                n
              
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




Otherwise, because of its close relationship with the sine function, the sidelobe peaks of the cardinal sine will be most likely located around odd multiples of




  ±
  
    π
    2
  




(except for the interval [−π, π] which only includes the main lobe). Since it is possible to give reliable roots approximations for the equation (7), it can be numerically solved using the Newton-Raphson iterative method for finding roots:





  
    
      
        
          
            
              
                x
                
                  n
                  +
                  1
                
              
              =
              
                
                  x
                  n
                
                -
                
                  
                    
                      
                        x
                        n
                      
                      .
                      
                        cos
                        ⁡
                        
                          (
                          
                            x
                            n
                          
                          )
                        
                      
                    
                    -
                    
                      sin
                      ⁡
                      
                        (
                        
                          x
                          n
                        
                        )
                      
                    
                  
                  
                    
                      d
                      
                        (
                        
                          
                            
                              x
                              n
                            
                            ·
                            
                              cos
                              ⁡
                              
                                (
                                
                                  x
                                  n
                                
                                )
                              
                            
                          
                          -
                          
                            sin
                            ⁡
                            
                              (
                              
                                x
                                n
                              
                              )
                            
                          
                        
                        )
                      
                    
                    /
                    
                      d
                      
                        x
                        n
                      
                    
                  
                
              
            
          
        
        
          
            
              =
              
                
                  x
                  n
                
                +
                
                  
                    
                      
                        x
                        n
                      
                      .
                      
                        cos
                        ⁡
                        
                          (
                          
                            x
                            n
                          
                          )
                        
                      
                    
                    -
                    
                      sin
                      ⁡
                      
                        (
                        
                          x
                          n
                        
                        )
                      
                    
                  
                  
                    
                      x
                      n
                    
                    ·
                    
                      sin
                      ⁡
                      
                        (
                        
                          x
                          n
                        
                        )
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




where nε[0,nmax] is the iterator, and

    
    
        x0 is the initialization root for every sidelobe of the function sinc(x).
    
    


The implementation scheme of the equation (9) gives the following points for the extrema of sinc(x):













Sidelobe
Peak location
Peak location in Hz
Peak value







Main
x = 0
f = 0
|sinc(x)| = 1


lobe





 


1st
x = ±4.4934095




  f
  =
  
    
      ±
      
        4.4934095
        πT
      
    
    =
    71.51
  




|sinc(x)| = 0.2172336


 


2nd
x = ±7.7252518




  f
  =
  
    
      ±
      
        7.7252518
        πT
      
    
    =
    122.95
  




|sinc(x)| = 0.1283746


 


3rd
x = ±10.904122




  f
  =
  
    
      ±
      
        10.904122
        πT
      
    
    =
    173.54
  




|sinc(x)| = 0.0913252









These points define frequencies that can be used as modulation frequencies fm to defeat shutter speed configurations of fs=1/T. As illustrated by a diagram of the FIG. 8, the sidelobes of the sampled modulated signal generates aliasing because the attenuation of the shutter spectrum gain is not strong enough around sidelobes to fully remove the modulated signal. A part of this signal will be able to go through the shutter and generates aliasing for the frequency band visible to the human eye.
FIG. 8 shows four diagrams in the frequency domain. A first diagram shows the spectrum of the base band signal X(f) representing the video signal. A second diagram shows the spectrum of the corresponding signal recorded by the camcorder. A third diagram shows the spectrum of the signal X(f) which is amplitude modulated on a carrier signal of frequency fm. The resulting signal is called XAM(f). A fourth diagram shows the spectrum of the sampled modulated signal when the modulation frequency is close to a peak location given in the previous table.
Consequently, according the invention, a set of pixels locations are selected in order to form an anti-copy pattern e.g a warning message. An amplitude modulation is applied to these selected pixels such that the amplitude of the brightness of the selected pixels is modulated at a rate substantially equal to one of the frequency values given for the sidelobes in the previous table. Consequently, the modulation frequency is chosen among the frequency values f different from zero for which the derivative of the function sinc(πfτ) is zero where τ is the shutter integration time of the camcorder. The chosen modulation frequency is thus a local extremum of the function sinc(πfτ).
For example, if the frequency of the first sidelobe




  (
  
    f
    =
    
      4.4934095
      
        π
        ⁢
        
          
        
        ⁢
        T
      
    
  
  )




is used as modulation frequency, the following modulation frequencies can be used:

    
    
        fm=71.51 Hz to overcome shutter settings of T= 1/50 in PAL camcorders;
        fm=85.82 Hz to overcome shutter settings of T= 1/60 in NTSC camcorders;
        fm=50.00 Hz in case shutter speed 1/T is set to higher values.

Of course, the modulation frequency can also be selected in a limited frequency range centered on the above-mentioned modulation frequencies fm, for example in the range [fm−Δ,fm+Δ] where Δ is comprised between 5 and 15 Hz. For example if Δ=10 Hz, fm is selected in the range [61.51 Hz, 81.51 Hz] to overcome shutter settings of T= 1/50 in PAL camcorders or [75.82 Hz, 95.82 Hz] to overcome shutter settings of T= 1/60 in NTSC camcorders

    
    


Other shutter configurations usually set exposure time T to lower values, and will let any of those two modulation frequencies go through.
Of course, the frequency of the second or third sidelobes can be used as modulation frequency.
An example of anti-copy pattern e.g a warning message is given in FIG. 10. Let M(x,y,t) be the brightness value (0 or 1) of the pixels of the warning message and C1(x,y,t), C2(x,y,t) and C3(x,y,t) the three color components of a pixel located on the coordinates (x,y) in space and the coordinate t in time. The triplet {C1,C2,C3} can be in a Red-Green-Blue color space (i.e. C1=Red, C2=Green and C3=Blue) or in a more HVS-friendly space, such as CIE XYZ or Yxy. The advantage of using non-RGB color spaces is to work in a linear theory-compliant domain.
For the modulation in the Red-Green-Blue color space, modulated pixel values {C1′,C2′,C3′} will be computed as follows:



  
    ∀
    
      i
      ∈
      
        [
        
          1
          ,
          3
        
        ]
      
    
  




  
    
      
        C
        i
        ′
      
      ⁡
      
        (
        
          x
          ,
          y
          ,
          t
        
        )
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    
                      C
                      i
                    
                    ⁡
                    
                      (
                      
                        x
                        ,
                        y
                        ,
                        t
                      
                      )
                    
                  
                  ·
                  
                    (
                    
                      1
                      +
                      
                        
                          β
                          i
                        
                        ⁢
                        
                          cos
                          ⁡
                          
                            (
                            
                              2
                              ⁢
                              π
                              ⁢
                              
                                
                              
                              ⁢
                              
                                f
                                m
                              
                              ⁢
                              t
                            
                            )
                          
                        
                      
                    
                    )
                  
                
                ⁢
                
                  
                
                ⁢
                if
                ⁢
                
                  
                
                ⁢
                
                  M
                  ⁡
                  
                    (
                    
                      x
                      ,
                      y
                      ,
                      t
                    
                    )
                  
                
              
              =
              1
            
          
        
        
          
            
              
                
                  C
                  i
                
                ⁡
                
                  (
                  
                    x
                    ,
                    y
                    ,
                    t
                  
                  )
                
              
              ⁢
              
                
              
              ⁢
              if
              ⁢
              
                
              
              ⁢
              not
            
          
        
      
    
  



where the modulation indexes βi are carefully chosen to be sure the pixel values do never reach saturation (for example 255 when the brightness values are coded in 8 bits).

FIG. 10 is a block diagram of a possible circuit implementation of the inventive method. The method is more particularly implemented in a modulator 10. It receives the pixel values of the images to be displayed. A signal P is also provided to the modulator. This signal indicates if the incoming pixel belongs to the anti-copy pattern or not.
Let I be the intensity of an incoming pixel (that is either R, G, B, or a combination of those) with a frame rate fr (e.g. fr=24 Hz), and I1, I2, I3 . . . the outcoming modulated pixels with a frame rate fr′. The amplitude of the luminance of the pixels of the anti-copy pattern is modulated with a modulation frequency fm. The frame rate fr′ is chosen to be equal or higher than twice the modulation frequency fm and preferably to be a multiple of fr (the conditions of the Nyquist-Shannon sampling theorem should be met by the video projector). For example, if fr′=3 fr, the modulator generates 3 pixel values I1, I2, I3 for the pixel value I. The pixel values I1, I2, I3 are defined such as:



  
    
      l
      n
    
    =
    
      l
      ⁡
      
        (
        
          1
          +
          
            β
            ⁢
            
              
            
            ⁢
            
              cos
              ⁡
              
                (
                
                  2
                  ⁢
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  
                    f
                    m
                  
                  ⁢
                  
                    n
                    /
                    
                      f
                      r
                      ′
                    
                  
                
                )
              
            
          
        
        )
      
    
  




  
    
      for
      ⁢
      
        
      
      ⁢
      n
    
    =
    
      1
      ⁢
      
        
      
      ⁢
      to
      ⁢
      
        
      
      ⁢
      3
      ⁢
      
        (
        
          =
          
            
              f
              r
              ′
            
            
              f
              r
            
          
        
        )
      
    
  



These modulated pixels are then displayed by a video projector 20.

## Claims

1. Method for processing a sequence of video images represented by their video signals comprising a modulation step for modulating temporally at a modulation frequency (fm) the amplitude of the brightness of a set of pixels of each image of the sequence around a brightness value to be displayed for said image, wherein said modulation is chosen to be not visible to the human eye while simultaneously generating artifacts due to aliasing when said video images are captured by a digital video capturing device,
wherein the modulation frequency (fm) is selected to be in a given frequency range centered on a local extremum of the function sinc(πfτ), where the function sinc is the cardinal sine function, τ is the shutter integration time of the digital video capturing device, and f is the frequency, wherein, if said local extremum corresponds to frequency F, said given frequency range extends from F−Δ to F+Δ, Δ being comprised between 5 and 15 Hz.

2. Method according to claim 1, wherein the shutter integration time τ is selected to be equal to the inverse of the sampling frequency fs of the digital video capturing device.

3. Method according to claim 2, wherein the modulation frequency is selected to be in said given frequency range centered on the frequency 71.5 Hz or 122.9 Hz when the inverse of the shutter integration time of the digital video capturing device is equal to 50 Hz.

4. Method according to claim 2, wherein the modulation frequency is selected to be in said given frequency range centered on the frequency 85.8 Hz or 147.5 Hz when the inverse of the shutter integration time of the digital video capturing device is equal to 60 Hz.

5. Device for processing a sequence of video images represented by their video signals comprising a modulator for modulating temporally at a modulation frequency (fm) the amplitude of the brightness of a set of pixels of each image of the sequence around a brightness value to be displayed for said image, wherein the modulation is chosen to be not visible to the human eye while simultaneously generating artifacts due to aliasing when said video images are captured by a digital video capturing device,
characterized in that the frequency (fm) used for the amplitude modulation is selected to be in given frequency range centered on a local extremum of the function sinc(πfτ) where the function sinc is the cardinal sine function, and τ is the shutter integration time of the digital video capturing device, and f is the frequency, wherein, if said local extremum corresponds to frequency F, said given frequency range extends from F−Δ to F+Δ, Δ being comprised between 5 and 15 Hz.

6. Device according to claim 5, wherein the shutter integration time τ is selected to be equal to the inverse of the sampling frequency fs of the digital video capturing device.

7. Device according to claim 6, wherein the modulation frequency is selected to be in said given frequency range centered on the frequency 71.5 Hz or 122.9 Hz when the inverse of the shutter integration time of the digital video capturing device is equal to 50 Hz.

8. Device according to claim 6, wherein the modulation frequency is selected to be in said given frequency range centered on the frequency 85.8 Hz or 147.5 Hz when the inverse of the shutter integration time of the digital video capturing device is equal to 60 Hz.

9. Device according to claim 5, wherein the images computed by the modulator are outputted at a frame rate that is equal to or higher than twice the modulation frequency.

10. Device according to claim 9, wherein the images computed by the modulator are outputted at a frame rate that is equal to a multiple of the frame rate at which the images are inputted into the modulator.

11. Device according to claim 5, wherein the digital video capturing device is a digital camcorder.

12. Method according to claim 1, wherein the digital video capturing device is a digital camcorder.

13. Method according to claim 2, wherein the digital video capturing device is a digital camcorder.

14. Method according to claim 3, wherein the digital video capturing device is a digital camcorder.

15. Method according to claim 4, wherein the digital video capturing device is a digital camcorder.