# Patent

## Title
Movement detection method and device with a thresholding algorithm per envelope

## Abstract

A movement detection method that includes calculating a first mean of a signal designed to be supplied by at least one pixel of a pixel matrix which corresponds to a captured image. The movement detection method further includes calculating a second signal, second mean, third mean and fourth mean wherein a movement is detected from the result of at least one comparison made of the signal and the second, third and fourth mean.

## Background
TECHNICAL FIELD
This document relates to the field of movement detection and more particularly that of image sensors, such as CMOS imaging devices used in the visible or infrared range, wherein a movement detection method is used.
STATE OF THE PRIOR ART
Movement detection involves detecting the movement of moving elements with respect to fixed elements in a field of captured images. These elements may be for example vehicles or even people. Such movement detection consists of isolating, among the signals supplied by an image sensor, those related to the moving elements, for example by detecting the significant variations on the mean or variance of a signal of a pixel or a group of pixels indicating a change in the nature of the element captured in this pixel or group of pixels, with respect to those related to the fixed elements of which the mean or variance remains for example substantially constant in time. For this purpose, detection methods or algorithms are used.
A first approach consists of using a “recursive average” algorithm for such movement detection. This algorithm is based on an estimated background calculation, which is to say of the fixed elements found in all of the captured images. Which is to say Xn is the background, or the mean, corresponding to an image n, Sn is the signal corresponding to the acquired image n and 1/N is a weighting coefficient, therefore:




  
    X
    n
  
  =
  
    
      X
      
        n
        -
        1
      
    
    -
    
      
        1
        N
      
      ⁢
      
        X
        
          n
          -
          1
        
      
    
    +
    
      
        1
        N
      
      ⁢
      
        S
        n
      
    
  




A comparison is then made between a chosen threshold value Th and |Sn−Xn|. If the value obtained is positive, this means that a movement has been detected. Xn and Sn may be variables obtained from a signal supplied by a pixel or by considering several signals supplied by several pixels, for example a group of pixels located next to one another that form a macropixel, like a single signal, taking for example the mean of these signals.
Such an algorithm has especially as disadvantage a lack of robustness in the detection as the detection threshold Th is determined à priori, prior to the algorithm being used, and is global for all of the pixels of the matrix. This disadvantage results in low precision of the location of the movements in the captured images. The low pass type filtering carried out by this algorithm induces dephasing, and consequently a delay in the response with respect to the signal. This delay results in a “drag effect” when the moving elements detected are shown, which is to say that the signals supplied by the pixels which no longer capture part of a moving element are still considered as signals which represent part of the moving element.
A second approach consists of using a “sigma-delta” algorithm for the movement detection. This algorithm permits significant variations of the signal to be detected by calculating two variables that can be assimilated to the mean value and the variance of the signal. FIG. 1 is a diagrammatical representation of a movement detection device using a sigma-delta algorithm.
Firstly, the sigma-delta mean M1n is calculated, with a constant incrementation and decrementation value, for example 1, of the signal Sn corresponding to the acquired image n. For this purpose, the signal Sn is sent as an input of the first means of calculating the sigma-delta mean 2. These means 2 first carry out an initialisation M10=S0. For the following images, which is to say for n>0, these means 2 compare M1n-1 and Sn. If M1n-1<Sn, then the value of M1n-1 is incremented such that M1n=M1n-1+1. If M1n-1>Sn, then the value of M1n-1 is decremented such that M1n=M1n-1−1. The value of the signal Δn=|M1n−Sn| is calculated by a subtractor 4 and absolute value calculation means 6. The calculation of N×Δn is then made by the multiplier 8, wherein N is a constant corresponding to the adaptive threshold of the algorithm whose value is chosen in function of the complexity of the scene. The calculation of a sigma-delta mean M2n of N·Δn is then made and sent as an input of second means of calculating the sigma-delta mean 10. Next, the initialisation M20=0 is carried out first. For the following images, which is to say for n>0, a comparison is made of M2n-1 and N·Δn. If M2n-1<N·Δn, then the value of M2n-1 is incremented such that M2n=M2n-1+1. If M2n-1>N·Δn, then the value of M2n-1 is decremented such that M2n=M2n-1−1. Finally, a comparison is made by a comparator 12 of the signal Δn and M2n. If M2n<Δn, this means that a movement has been detected.
As for the recursive average algorithm, the variables Sn, M1n, Δn and M2n may be obtained from a signal supplied by a pixel or by considering several signals supplied by several pixels, for example a group of pixels next to one another, like a single signal, taking for example the mean of these signals.
However, such an algorithm especially has the disadvantage of not filtering high frequency parasite movements which are considered as movements to be detected (for example, a movement of the leaves of a tree or snow falling). Furthermore, the constant N used must be determined à priori, which reduces the adaptability of the detection carried out by this algorithm.
DESCRIPTION OF THE INVENTION
Thus there is a need to propose a method of movement detection which permits the detection of high frequency parasite movements to be reduced or eliminated and which offers more efficient detection, for example in terms of precision of locating the movements, with respect to the methods of the prior art, and which reduces or eliminates the “drag” effect obtained by the methods of the prior art.
Also, there is a need to propose a method of movement detection which requires few calculation and memory hardware resources to be used, and that can be installed analogically in a very low consumption imaging device (with for example a mean consumption equal to approximately several hundred μW).
For this purpose, one embodiment proposes a method of movement detection comprising at least the following steps:

    
    
        the calculation of a first mean M1n of a signal Sn designed to be supplied by at least one pixel of a pixel matrix which corresponds to an n-th captured image, in function of the value of the signal Sn and/or a previous value M1n-1;
        the calculation of a signal Δn=|M1n−Sn|;
        the calculation of a second mean M2n of the signal Δn in function of the value of a previous value M2n-1 and/or the value of the signal Δn;
        the calculation of a third mean M3n of the signal (Sn−k1·M2n) in function of a previous value M3n-1 and/or the value of the signal (Sn−k1·M2n);
        the calculation of a fourth mean M4n of the signal (Sn+k2·M2n) in function of a previous value M4n-1 and/or the value of the signal (Sn+k2·M2n);
    
    


wherein a movement is considered as detected from the result of at least one comparison made on the signals Sn, (M3n−M2n) and (M4n+M2n);
where n: natural whole number;
and k1 and k2: non-null positive real numbers.
Consequently, the significant variations of the signal emitted by a pixel or a macropixel may be detected by “thresholding” the variations of the signal by the envelope of the signal Sn itself, which is to say by defining adaptive upper and lower threshold values for the movement detection depending on the activity of the captured signal Sn. Moreover, this thresholding can keep the sign of the emitted signal derivative, which it is not possible with a thresholding realized on the absolute value of the difference between the emitted signal and an estimated background.
The adaptability of the detection is also improved by eliminating certain constants that had to determined à priori in the methods of the prior art. The sensitivity of the detection is also adapted locally to the activity of the pixels, which is to say individually for each pixel.
The coefficients k1 and k2 allow a direct adaptation of the upper and lower threshold values depending on the disturbance level in the analyzed images. These coefficients allow an expansion or a reduction of the envelope defined by the detection threshold around Sn. For example, it is possible to impose “harder” threshold detection conditions with an expansion of the envelope, thus enables to escape some noise disturbances during movement detection. The values of k1 and k2 may be chosen in function of the environment in which the movement detection is realized, that is in function of the ambient “agitation” in the analyzed images. In an urban environment, the values of k1 and k2 may be equals to 1. Preferably, the values of k1 and/or k2 may be comprised between 0 and 3. Moreover, the values of k1 and k2 may be similar or not, depending on whether a symmetrical or an unsymmetrical envelope is needed, in function of disturbances that are not to be detected as movements.
The first mean M1n may be obtained from the following equation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        1
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          1
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            1
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          1
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            1
          
        
        ⁢
        
          S
          n
        
      
    
  
  ,




where M1−1=0, and
1/N1: non-null positive real number.
By choosing an appropriate value of 1/N1, the time constant of the calculation of this first mean may be chosen, which may be slow to estimate the background of the images captured.
In one variant, the first mean M1n may be obtained at least by the following calculation steps:

    
    
        M10=S0;
    
    


and for n>0:

    
    
        M1n=M1n-1+c1 when M1n-1<Sn;
        M1n=M1n-1−c1 when M1n-1>Sn;
    
    


Where c1: non-null positive real number.
The second mean M2n may be obtained from the following equation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        2
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          2
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            2
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          2
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            2
          
        
        ⁢
        
          Δ
          n
        
      
    
  
  ,




where M2−1=0, et
1/N2: non-null positive real number.
The second mean M2n may be also obtained at least by the following calculation steps:

    
    
        M20=Δ0;
    
    


and for n>0:

    
    
        M2n=M2n-1+c2 when M2n-1<Δn;
        M2n=M2n-1−c2 when M2n-1>Δn;
    
    


where c2: non-null positive real number.
The third mean M3n may be obtained from the following equation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        3
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          3
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            3
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          3
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            3
          
        
        ⁢
        
          (
          
            
              S
              n
            
            -
            
              
                k
                1
              
              ⁢
              M
              ⁢
              
                
              
              ⁢
              
                2
                n
              
            
          
          )
        
      
    
  
  ,




Where M3−1=0, and
1/N3: non-null positive real number.
In this case, the value of 1/N3 may be adjusted in function of the value |M3n−Sn|.
The third mean M3n may be obtained at least by the following calculation steps:

    
    
        M30=(S0−k1·M20);
    
    


and for n>0:

    
    
        M3n=M3n-1+c3 when M3n-1<(Sn−k1·M2n);
        M3n=M3n-1−c3 when M3n-1>(Sn−k1·M2n);
    
    


where c3: non-null positive real number.
The fourth mean M4n may be obtained from the following equation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        4
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          4
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            4
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          4
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            4
          
        
        ⁢
        
          (
          
            
              S
              n
            
            +
            
              
                k
                2
              
              ⁢
              M
              ⁢
              
                
              
              ⁢
              
                2
                n
              
            
          
          )
        
      
    
  
  ,




where M4−1=0, and
1/N4: non-null positive real number.
In this case, the value of 1/N4 may be adjusted in function of the value |M4n−Sn|.
An asymmetry of the thresholding envelope may be realised by the choice of the values of time constants 1/N3 and 1/N4 which may be dynamically adapted during the method.
For example, when M4n<Sn, N4 may be small, which corresponds to a short time constant τ, for example lower than about 0.5 second (for example equal to 380 ms). When M4n>Sn, N4 may be high, which corresponds to a long time constant τ, for example higher than about 1 second (for example equal to 1.18 s), at a sampling period equal to 25 Hz. The time constant τ of a mean may be equal to the ratio: sampling period of the image capture/ln(1−1/N)−1, the sampling period corresponding to the period of the signal Sn period (noted period (Sn)). Thus, when M4n<Sn, the value of 1/N4 may verify the relation:
period (Sn)/ln(1−1/N4)−1<0.5 s;
and when M4n>Sn, the value of 1/N4 may verify the relation:
period (Sn)/ln(1−1/N4)−1>1 s.
Moreover, when M3n<Sn, N3 may be high, which corresponds to a long time constant τ, for example higher than about 1 second, and when M3n>Sn, N3 may be small, which corresponds to a short time constant τ, for example smaller than about 0.5 second. Thus, when M3n<Sn, the value of 1/N3 may verify the relation:
period (Sn)/ln(1−1/N3)−1>1 s;
and when M3n>Sn, the value of 1/N3 may verify the relation:
period (Sn)/ln(1−1/N3)−1<0.5 s.
A short time constant may be equal to about 380 ms, and a long time constant may be equal to about 1.18 s. The time constants are here considered at a sampling period of 25 Hz.
Thus, an unsymmetrical thresholding envelope allows adapting the thresholding of positive variations of the signal (comparison of Sn and M4n) and the thresholding of negative variations of the signal (comparison of Sn and M3n). For example, it is possible to realize a filtering of a camera blooming with repeated flashes by imposing “harder” threshold detection conditions on the variations (for example positive in the case of a grey-levels scene and white flashes corresponding to a maximum value of the output signal) induced by flashes et keeping a suitable sensitivity on negatives variations induces by the other darker objects.
The fourth mean M4n may be obtained at least by the following calculation steps:

    
    
        M40=(S0+k2·M20);
    
    


and for n>0:

    
    
        M4n=M4n-1+c4 when M4n-1<(Sn+k2·M2n);
        M4n=M4n-1−c4 when M4n-1>(Sn+k2·M2n);
    
    


where c4: non-null positive real number.
The value of the first mean M1n and/or the second mean M2n and/or the third mean M3n and/or the fourth mean M4n may be respectively greater than a first non-null minimum threshold value SM1n and/or a second non-null minimum threshold value SM2n and/or a third non-null minimum threshold value SM3n and/or a fourth non-null minimum threshold value SM4n.
These threshold values may be comprised between approximately 1/250×the dynamic of the signal Sn and 1/25×the dynamic of the signal Sn, and for example equal to approximately 1/50×the dynamic of the signal Sn, the dynamic of the signal Sn corresponding to the possible maximum value of |Sn| (for example equal to 256 for a 8 bits coded signal).
The comparison may be made between the signals Sn and (M3n−M2n) and between the signals Sn and (M4n+M2n), wherein a movement is considered as detected when Sn>(M4n+M2n) or Sn<(M3n−M2n).
In one variant, the method may further comprise, between the calculation step of the fourth mean M4n of the signal (Sn+k2·M2n) and the comparison step, a calculation step of a fifth mean M5n of the signal Sn in function of a previous value M5n-1 and/or the value of the signal Sn, wherein said comparison may be made between the signals M5n and (M3n−M2n) and between the signals M5n and (M4n+M2n), wherein a movement is considered as detected when M5n>(M4n+M2n) or M5n<(M3n−M2n).
In this case, band pass filtering may be generated to eliminate the high frequency parasite movements, as well as the “drag” effect observed in the movement detection methods of the prior art.
The fifth mean M5n may be obtained from the following equation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        5
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          5
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            5
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          5
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            5
          
        
        ⁢
        
          S
          n
        
      
    
  
  ,




where M5−1=0, and
1/N5: non-null positive real number.
In this case, the high frequency parasites of the signal are filtered even more efficiently.
The values of the weighting coefficients 1/N, and the values of the incrementation and decrementation coefficients c in the case of a calculation of a sigma-delta type mean, may be calculated in function of the speed of acquisition of the images, and may be easily determined.
The values of the weighting coefficients 1/N which may be used for the calculation of means M1n and/or M2n may be chosen such that the time constants τ of these means are long, for example higher than about 1 second.
The values of the incrementation and decrementation coefficients c which may be used for the calculation of means M1n and/or M2n and/or M3n and/or M4n may be chosen and adapted during the method in order to satisfy the relation c<|dSn/dn|. In one variant, when the incrementation and decrementation coefficients have fixed values, for example equal to 1, it is possible to adapt a refresh period Tn of the method, that is the period of which the steps of the method are realized, in order to satisfy the relation c<|ΔSn/Tn|.
Another embodiment also relates to a movement detection device comprising at least:

    
    
        means of calculating, or a calculator of, a first mean M1n of a signal Sn designed to be supplied by at least one pixel of a pixel matrix and corresponding to a n-th captured image, in function of the value of the signal Sn and/or a previous value M1n-1;
        means of calculating, or a calculator of, a signal Δn=|M1n−Sn|;
        means of calculating, or a calculator of, a second mean M2n of the signal Δn in function of a previous value M2n-1 and/or the value of the signal Δn;
        means of calculating, or a calculator of, a third mean M3n of the signal (Sn−k1·M2n) in function of a previous value M3n-1 and/or the value of the signal (Sn−k1·M2n);
        means of calculating, or a calculator of, a fourth mean M4n of the signal (Sn+k2·M2n) in function of a previous value M4n-1 and/or the value of the signal (Sn+k2·M2n);
    
    


a movement is considered as detected from the result at the output of the means of comparison, or a comparator, making at least one comparison of the signals Sn, (M3n−M2n) and (M4n+M2n);
where n: natural whole number;
and k1 and k2: non-null positive real numbers.
The means of calculating, or the calculator of, the first mean M1n may carry out at least the following operation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        1
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          1
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            1
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          1
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            1
          
        
        ⁢
        
          S
          n
        
      
    
  
  ,




where M1−1=0, and
1/N1: non-null positive real number.
The means of calculating, or the calculator of, the first mean M1n may at least comprise:

    
    
        means of initialising the value of M10 to the value of S0;
        means of comparing, or a comparator of, the value of the signal M1n-1 and the value of the signal Sn;
        means of incrementing and decrementing the value of M1n by a constant c1, where c1: non-null positive real number.
    
    


The means of calculating, or the calculator of, the second mean M2n may carry out at least the following operation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        2
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          2
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            2
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          2
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            2
          
        
        ⁢
        
          Δ
          n
        
      
    
  
  ,




where M2−1=0, and
1/N2: non-null positive real number.
The means of calculating, or the calculator of, the second mean M2n may at least comprise:

    
    
        means of initialising the value of M20 to the value of Δ0;
        means of comparing, or a comparator of, the value of the signal M2n-1 and the value of the signal Δn;
        means of incrementing and decrementing the value of M2n by a constant c2, where c2: non-null positive real number.
    
    


The means of calculating, or the calculator of, the third mean M3n may carry out at least the following operation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        3
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          3
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            3
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          3
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            3
          
        
        ⁢
        
          (
          
            
              S
              n
            
            -
            
              
                k
                1
              
              ⁢
              M
              ⁢
              
                
              
              ⁢
              
                2
                n
              
            
          
          )
        
      
    
  
  ,




where M3−1=0, and
1/N3: non-null positive real number.
In this case, the means of calculating, or the calculator of, the third mean M3n may further comprise means of adjusting the value of 1/N3 in function of the value |M3n−Sn| such that:

    
    
        when M3n<Sn, the value of 1/N3 verify the relation: period (Sn)/ln(1−1/N3)−1>1 s;
    
    


and when M3n>Sn, the value of 1/N3 verify the relation: period (Sn)/ln(1−1/N3)−1<0.5 s.
The means of calculating, or the calculator of, the third mean M3n at least comprise:

    
    
        means of initialising the value of M30 to the value of (S0−k1·M20);
        means of comparing, or a comparator of, the value of the signal M3n-1 and the value of the signal (Sn−k1·M2n);
        means of incrementing and decrementing the value of M3n by a constant c3, where c3: non-null positive real number.
    
    


The means of calculating, or the calculator of, the fourth mean M4n may carry out at least the following operation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        4
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          4
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            4
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          4
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            4
          
        
        ⁢
        
          (
          
            
              S
              n
            
            +
            
              
                
                  k
                  2
                
                ·
                M
              
              ⁢
              
                
              
              ⁢
              
                2
                n
              
            
          
          )
        
      
    
  
  ,




where M4−1=0, and
1/N4: non-null positive real number.
In this case, the means of calculating, or the calculator of, the fourth mean M4n may further comprise means of adjusting the value of 1/N4 in function of the value |M4n−Sn| such that:

    
    
        when M4n<Sn| the value of 1/N4 verify the relation: period (Sn)/ln(1−1/N4)−1<0.5 s;
    
    


and when M4n>Sn, the value of 1/N4 verify the relation: period (Sn)/ln(1−1/N4)−1>1 s.
The means of calculating, or the calculator of, the fourth mean M4n may at least comprise:

    
    
        means of initialising the value of M40 to the value of (S0+k2·M20);
        means of comparing, or a comparator of, the value of the signal M4n-1 and the value of the signal (Sn+k2·M2n);
        means of incrementing and decrementing the value of M4n by a constant c4, where c4: non-null positive real number.
    
    


The means of comparison, or the comparator, may be means of comparing, or a comparator of, the signals Sn and (M3n−M2n) and means of comparing, or a comparator of, the signals Sn and (M4n+M2n), wherein a movement may be considered as detected when Sn>(M4n+M2n) or Sn<(M3n−M2n).
The device may further comprise means of calculating, or a calculator of, a fifth mean M5n of the signal Sn in function of a previous value M5n-1 and/or the value of the signal Sn, wherein the means of comparison, or the comparator, may be means of comparing, or a comparator of, the signals M5n and (M3n−M2n) and means of comparing, or a comparator of, the signals M5n and (M4n+M2n), wherein a movement is considered as detected when M5n>(M4n+M2n) or M5n<(M3n−M2n).
In this case, the means of calculating, or the calculator of, the fifth mean M5n may carry out at least the following operation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        5
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          5
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            5
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          5
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            5
          
        
        ⁢
        
          S
          n
        
      
    
  
  ,




where M5−1=0, and
1/N5: non-null positive real number.
The means of comparison, or the comparator, may at least comprise an operational amplifier, or transconductance amplifier.
Finally, this invention also relates to an image capture device comprising at least one pixel matrix and one movement detection device as previously described.



BRIEF DESCRIPTION OF THE DRAWINGS
This invention will be more clearly understood upon reading the following description of embodiments provided purely by way of illustration and in no way restrictively, in reference to the appended drawings in which:
FIG. 1 shows diagrammatically a sigma-delta algorithm movement detection device,
FIG. 2 shows diagrammatically a movement detection device with a thresholding algorithm per envelope,
FIG. 3 shows signals obtained during the implantation of a movement detection method with a thresholding algorithm per envelope,
FIG. 4 shows part of an image capture device comprising an example of a movement detection device with a thresholding algorithm per envelope.



Identical, similar or equivalent parts of the various figures described below have the same numerical references so as to facilitate the passage from one figure to another.
The different parts shown in the figures are not necessarily to a uniform scale, to make the figures easier to read.
The different possibilities (variants and embodiments) should be understood as not being mutually exclusive and may be combined with one another.
DETAILED DESCRIPTION OF SPECIFIC EMBODIMENTS
A movement detection device 100 and method with a thresholding algorithm per envelope, according to one specific embodiment will now be described in relation to FIG. 2.
A first step of this method is to estimate the background of the images in which a movement is to be detected, wherein the background is formed by the fixed elements in the captured images. For this purpose, firstly, using means 102, a first recursive average M1n is calculated wherein an important weighting coefficient 1/N1 of a signal Sn corresponding to the signal supplied by a pixel, or a group of pixels also called macropixel, of an n-th acquired image. Wherein:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        1
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          1
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            1
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          1
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            1
          
        
        ⁢
        
          S
          n
        
      
    
  
  ,




In this calculation of the first mean M1n, the choice of the value of the weighting coefficient 1/N1 is important so that only the fixed elements in the acquired images are conserved. For an image capture device operating at 25 Hz, for example a weighting coefficient such as 1/N1=1/(27 is chosen. The value of the weighting coefficient 1/N1 is chosen in function of the amplitude and the frequency of the variations of the signal Sn.
Consequently, the value M1n calculated for a pixel or a macropixel corresponds to the luminous intensity emitted by the fixed element captured by this pixel or this macropixel, even if one or several moving elements temporarily pass before this pixel or this macropixel during the movement detection.
In one variant of this embodiment, the first mean M1n may not be a recursive average, but a sigma-delta mean, for example with a low level of incrementation and decrementation c1 (for example equal to 1), also permitting the estimation of the background to be calculated. A similar effect may be obtained by choosing a low refresh period. This low level of incrementation and decrementation corresponds to an important time constant in a recursive mean calculation. For this purpose, firstly the initialisation M10=S0 is carried out. For the following images, which is to say for n>0, M1n-1 and Sn are compared. If M1n-1<Sn, then the value of M1n-1 is incremented such that M1n=M1n-1+c1. If M1n-1>Sn, then the value of M1n-1 is decremented such that M1n=M1n-1−c1.
Then the signal Δn=|M1n−Sn| is calculated using a subtractor 104 and absolute value calculation means 106.
Then a second sigma-delta mean M2n, with a constant level of incrementation and decrementation c2 (for example equal to 1) of the signal Δn is calculated using calculating means 108.
To calculate this second mean M2n, firstly the initialisation M20=Δ0 is carried out. For the following images, which is to say for n>0, M2n-1 and Δn are compared. If M2n-1<Δn, then the value of M2n-1 is incremented such that M2n=M2n-1+c2. If M2n-1>Δn, then the value of M2n-1 is decremented such that M2n=M2n-1−c2.
It is also possible that this second mean M2n is not a sigma-delta type mean, but a recursive type mean. In this case, the calculation means 108 calculate the second recursive mean M2n with a high time constant 1/N2, of the signal Δn. This second recursive mean M2n is obtained with the following calculation:




  
    M
    ⁢
    
      
    
    ⁢
    
      2
      n
    
  
  =
  
    
      M
      ⁢
      
        
      
      ⁢
      
        2
        
          n
          -
          1
        
      
    
    -
    
      
        1
        
          N
          2
        
      
      ⁢
      M
      ⁢
      
        
      
      ⁢
      
        2
        
          n
          -
          1
        
      
    
    +
    
      
        1
        
          N
          2
        
      
      ⁢
      
        Δ
        n
      
    
  




The weighting coefficient 1/N2 is chosen in function of the amplitude and the frequency of the variations of the envelope signal Δn. For example N2=26.
Then the signal Sn−k1·M2n is calculated using a subtracter 110 and a multiplier 111. In this embodiment, the value of k1 is equal to 1.
The device 100 also comprises means 114 of calculating a third recursive mean M3n. This third recursive mean M3n is obtained using the following calculation:




  
    M
    ⁢
    
      
    
    ⁢
    
      3
      n
    
  
  =
  
    
      M
      ⁢
      
        
      
      ⁢
      
        3
        
          n
          -
          1
        
      
    
    -
    
      
        1
        
          N
          3
        
      
      ⁢
      M
      ⁢
      
        
      
      ⁢
      
        3
        
          n
          -
          1
        
      
    
    +
    
      
        1
        
          N
          3
        
      
      ⁢
      
        (
        
          
            S
            n
          
          -
          
            
              k
              1
            
            ·
            
              2
              n
            
          
        
        )
      
    
  




The weighting coefficient 1/N3 is chosen in function of the amplitude and the frequency of the variations of the signal (Sn−k1·M2n). For example N3=26.
This third mean M3n may also be obtained by calculating the delta-sigma mean of the signal (Sn−k1·M2n) with an incrementation constant c3.
Then the signal Sn+k2·M2 is calculated using an adder 112 and a multiplier 113. In this embodiment, the value of k1 is equal to 1.
Then, calculating means 116 permit a fourth recursive mean M4n to be calculated. This fourth recursive mean M4n is obtained by the following calculation:




  
    M
    ⁢
    
      
    
    ⁢
    
      4
      n
    
  
  =
  
    
      M
      ⁢
      
        
      
      ⁢
      
        4
        
          n
          -
          1
        
      
    
    -
    
      
        1
        
          N
          4
        
      
      ⁢
      M
      ⁢
      
        
      
      ⁢
      
        4
        
          n
          -
          1
        
      
    
    +
    
      
        1
        
          N
          4
        
      
      ⁢
      
        (
        
          
            S
            n
          
          +
          
            
              
                k
                2
              
              ·
              M
            
            ⁢
            
              
            
            ⁢
            
              2
              n
            
          
        
        )
      
    
  




The time constant 1/N4 is chosen in function of the amplitude and the frequency of the variations of the signal (Sn+k2·M2n). For example N4=26.
It is possible that the values of 1/N3 and/or 1/N4 are constant throughout the movement detection method. It is also possible that one or both of these two values vary during the movement detection method in function of the value of |M3n−Sn| for M3n and |M4n−Sn| for M4n. Consequently, a time constant 1/N3 may be chosen that is more rapid when M3n>Sn, or even a time constant 1/N4 that is slower when Sn>M4n.
In one variant, the third mean and/or the fourth mean may also be of the sigma-delta type, with incrementation and decrementation constants c3 and c4.
Preferably, a third mean M3n and a fourth mean M4n of the recursive types will be chosen.
Finally, to determine the presence or absence of movements in the captured image n, the device 100 comprises two comparators 122 and 124 which make a comparison of the signals (M3n−M2n), obtained by a subtractor 118, and Sn, and of the signals (M4n+M2n), obtained by an adder 120, and Sn, wherein a movement is considered as detected when Sn>(M4n+M2n) or Sn<(M3n−M2n).
In one variant, it is also possible that the presence or absence of movements in the captured image n is determined by comparing the signals (M4n+M2n) and (M3n−M2n) not with Sn but with a fifth recursive type mean M5n of this signal Sn. This fifth recursive type mean M5n is obtained by the following calculation:




  
    M
    ⁢
    
      
    
    ⁢
    
      5
      n
    
  
  =
  
    
      M
      ⁢
      
        
      
      ⁢
      
        5
        
          n
          -
          1
        
      
    
    -
    
      
        1
        
          N
          5
        
      
      ⁢
      M
      ⁢
      
        
      
      ⁢
      
        5
        
          n
          -
          1
        
      
    
    +
    
      
        1
        
          N
          5
        
      
      ⁢
      
        S
        n
      
    
  




The weighting coefficient 1/N5 used for this calculation is preferably rapid (for example N5=26). The weighting coefficient 1/N5 is chosen in function of the amplitude and the frequency of the variations of the signal Sn. In this case, the device 100 comprises means of calculating this fifth mean, not shown, of which the output is sent to the input of the comparators 122 and 124 instead of the signal Sn.
This method is based on the generation of an envelope whose upper and lower values form threshold values which permit the presence or absence of movement to be determined. The thresholding of the variations detected is therefore adaptive.
In general, the incrementation and decrementation values c1 and/or c2 and/or the time constants 1/N1, 1/N2, 1/N3, 1/N4 and 1/N5 used in the means calculations are adapted in function of the pixel resolution of the captured images, which is to say the number of levels of grey onto which the processed signal is encoded, as well as the operating frequency of the movement detection device 100.
FIG. 3 shows an example of a signal Sn from several images captured by a pixel matrix as well as the different signals calculated during a movement detection method with a thresholding algorithm per envelope as previously described.
In this FIG. 3, the x axis shows the evolution of the signals, graduated in the number of captured images, and the y axis shows the value of these signals graduated in the levels of grey. The curve 130 shows the signal Sn obtained at the output of a pixel or a macropixel. It is this signal that is sent to the input of the device 100. The curve 132 shows the first mean M1n, in this case recursive, obtained at the output of the calculation means 102. This first mean M1n shows the background of the image captured, which is to say the fixed element(s) captured by the pixel or the group of pixels. In FIG. 3, it can be seen that this first mean M1n varies very little, which effectively corresponds to the value of the fixed elements captured. The curve 134 shows the signal Δn. It may be seen in FIG. 3 that this signal Δn follows the most significant variations of the output signal Sn according to the curve 130, thus creating an estimation of the variations, which is to say of the moving elements in the captured images. Finally, the curves 136, 138 and 140 respectively show the signals M2n, M3n and M4n. The curves 138 and 140 show the upper and lower boundaries of the thresholding envelope calculated around the signal Sn. Consequently, by correctly choosing the value of the weighting coefficient for the calculation of the second mean M2n, the high frequency parasite movements corresponding to the low variations recorded by the image capture device are not considered as detected movements.
One example of an image capture device 200, or optical imaging device, permitting a movement detection method with a thresholding algorithm per envelope previously described to be implemented is shown in FIG. 4.
The image capture device 200 comprises pixel matrix 202 and a movement detection device 204. Each pixel of the matrix 202 is here formed by a photodiode and addressing and reading transistors. The movement detection method with a thresholding algorithm per envelope previously described may be applied to the signal Sn supplied by a pixel, by connecting a movement detection device similar to the device 204 to each column of pixels of the matrix 202. It is also possible that the movement detection method is applied to a signal Sn corresponding to the signals supplied by several pixels, for example the mean of these signals. Consequently, by considering the macropixels, it is possible to operate the image capture device 200 in low resolution zones, and only to detail these zones by using a movement detection method with a thresholding algorithm per envelope for each pixel of a macropixel when a movement is detected on this macropixel. It is therefore possible to reduce the number of movement detection devices 204 used in the image capture device 200 by only using a single movement detection device 204 per column of macropixels. A macropixel may for example be a square of 12×12 pixels, or any other value, for example 4×4 pixels as in FIG. 4. The values of each macropixel of the matrix 202 are read line by line.
The movement detection devices 204 shown in FIG. 4 comprises a comparator 206, for example an operational amplifier (or transconductance amplifier), a plurality of switched capacities 208, analogue memory registers 210 (four memory registers 210 are shown in FIG. 4 but the movement detection device 204 may comprise a different number of memory registers adapted to the number of values to be stored while the movement detection method is in use), an address demultiplexer 212a, an address multiplexer 212b, a multiplexer 214 for controlling the writing in the switched capacities 208, a multiplexer 216 supplying the values to be written in the switched capacities 208, a SRAM memory 218, and a capacitor 222.
The operation of the movement detection device 204 will now be described in relation to the implementation of the movement detection method previously described.
In the case of a first mean M1n of the recursive type, the value of the signal Mn-1 stored in one of the memory registers 210 is supplied to the input of the multiplexer 216, by means of the address multiplexer 212b. When n=0, M−1=0. The first recursive mean M1n is then calculated. For this purpose, the signal Mn-1 is sent to the terminals of a first switched capacity 208 whose value is equal to (N1-1)/N1. The value of the signal Mn-1 is stored at the terminals of this capacity.
The signal Sn supplied by the first macropixel of the matrix 202 is then supplied to the input of the multiplexer 216. The value of the signal Sn is stored at the terminals of a second switched capacity whose value C4 is equal to 1/N1.
By connecting in parallel the two switched capacities, the signal




  
    M
    ⁢
    
      
    
    ⁢
    
      1
      n
    
  
  =
  
    
      M
      ⁢
      
        
      
      ⁢
      
        1
        
          n
          -
          1
        
      
    
    -
    
      
        1
        
          N
          1
        
      
      ⁢
      M
      ⁢
      
        
      
      ⁢
      
        1
        
          n
          -
          1
        
      
    
    +
    
      
        1
        
          N
          1
        
      
      ⁢
      
        S
        n
      
    
  




is obtained at the terminals of these capacities.

It is also possible that the mean M1n is a sigma-delta mean. For this purpose, the initialisation M10=S0 is carried out by storing the value of the signal S0 in one of the memory registers 210. For the following images, M1n-1 and Sn are compared by applying one of the two signals to the positive input of the comparator 206 by means of the capacitor 222 in which this signal is stored, and by applying the other of the two signals to the negative input of the comparator 206. The result obtained at the output of the comparator 206 permits one of the values +c1 or −c1 applied to the input of the multiplexer 216 (+c and −c1 in FIG. 4) to be selected. If M1n-1<Sn, then the value of M1n-1 is incremented such that M1n=M1n-1+c1. If M1n-1>Sn, then the value of M1n-1 is decremented such that M1n=M1n-1−c1. In this embodiment, c1=1. The addition operation of +/−c1 to M1n-1 is carried out by means of the capacities 208 by storing at the terminals two of these capacities the +/−c1 and M1n-1 values, and by connecting these capacities in series so as to carry out an addition.
Whether the mean M1n is of the recursive or sigma-delta type, its value is stored in one of the memory registers 210. In this embodiment, the movement detection device 204 comprises four memory registers 210 designed to store four different means M1n, M2n, M3n and M4n.
The signal Δn=|M1n−Sn| is then calculated. This operation may be carried out by the switched capacities 208.
In a similar manner to the calculation of the first mean M1n when the latter is of the recursive type, the calculation of the second recursive type or sigma-delta mean M2n is carried out. When M2n is of the recursive type, then two other switched capacities 208 are used respectively of values (N2−1)/N2 and 1/N2. When the second mean M2n is of the sigma-delta type, this second mean is calculated in a similar manner to the calculation of the first mean M1n, using two values +c2 and −c2 applied to the input of the multiplexer 216 instead of the values +c1 and −c1.
Then the calculation (Sn−k1·M2n) is carried out by the switched capacities 208, then that of the third mean M3n of this signal. Finally, the signal (Sn+k2·M2n) and that of the fourth mean M4n of this signal are calculated.
A comparison is then made using the amplifier 206 of the signal Sn and the signal (M3n−M2n). The value obtained at the output of the amplifier 206 is representative of a detection or non-detection of a movement on the macropixel considered. If Sn<(M3n−M2n), then it is considered that a movement has been detected on the macropixel considered.
Next, a comparison is made using the amplifier 206 of the signal Sn and the signal (M4n+M2n). The value obtained at the output of the amplifier 206 is representative of a detection or non detection of a movement on the macropixel considered. If Sn>(M4n+M2n), then it is considered that a movement has been detected on the macropixel considered.
The operation is then repeated for the following macropixels, line after line.
In one variant, the comparison of the signal Sn and the signal (M3n−M2n) may be replaced by a comparison of a fifth recursive type mean M5n of this signal Sn where the signal (M3n−M2n). The value obtained at the output of the amplifier 206 is representative of a detection or non-detection of a movement on the macropixel considered. If M5n<(M3n−M2n), then it is considered that a movement has been detected on the macropixel considered. This fifth recursive mean M5n is obtained by the following calculation:




  
    M
    ⁢
    
      
    
    ⁢
    
      5
      n
    
  
  =
  
    
      M
      ⁢
      
        
      
      ⁢
      
        5
        
          n
          -
          1
        
      
    
    -
    
      
        1
        
          N
          5
        
      
      ⁢
      M
      ⁢
      
        
      
      ⁢
      
        5
        
          n
          -
          1
        
      
    
    +
    
      
        1
        
          N
          5
        
      
      ⁢
      
        S
        n
      
    
  




This calculation may be carried out by the switched capacities of value (N5−1)/N5 and 1/N5.
In the same way, the comparison of Sn and the signal (M4n+M2n) may be replaced by a comparison of M5n and (M4n+M2n). If M5n>(M4n+M2n), then it is considered that a movement has been detected on the macropixel considered.
For the mean M2n, a non-null minimum threshold value may be imposed, below which this mean is not allowed to descend. This minimum threshold value SM2n may be applied to the input of the multiplexer 216. When the value of this mean is below this threshold value, this value is then replaced by the threshold value. For example, in this embodiment, the value of SM2n may be equal to about 1/50×the dynamic of Sn.
When a movement is detected on the macropixel, it is possible, for an image n for which the movement detection has been carried out on macropixels, to store the values of the macropixels in the memory 218, then, on the macropixel(s) where movements are recorded, to implement the movement detection methods previously described for each of the pixels of the macropixel. Consequently, the location of the movements detected may be defined precisely, without processing all of the pixels of the captured images.
It may be seen that the method may be implemented using few hardware calculation (an operational amplifier, several switched capacities with a clock frequency for the instructions of several tens of kHz and several multiplexers/demultiplexers) and memory resources (several analogue registers per pixel and a SRAM memory for example).
The device shown in FIG. 4 permits an analogue implementation of the movement detection method previously described. However, it is also possible to use a digital implementation of the movement detection methods by connecting the pixel matrix 202 to signal digital processing means, for example a circuit of the DSP or FPGA type or a microprocessor, wherein the movement detection method is programmed.
The signals obtained at the output of the movement detection devices may be used to display an image on which the background captured forms a black background onto which the moving elements detected are shown in white, for example as shown on the image on the right in FIG. 4.

## Claims

1. A movement detection method comprising at least the following steps:
the calculation of a first mean M1n of a signal Sn designed to be supplied by at least one pixel of a pixel matrix which corresponds to an n-th captured image, in function of the value of the signal Sn and/or a previous value M1n-1;
the calculation of a signal Δn=|M1n−Sn|;
the calculation of a second mean M2n of the signal Δn in function of a previous value M2n-1 and/or the value of the signal Δn;
the calculation of a third mean M3n of the signal (Sn−k1·M2n) in function of a previous value M3n-1 and/or the value of the signal (Sn−k1·M2n);
the calculation of a fourth mean M4n of the signal (Sn+k2·M2n) in function of a previous value M4n-1 and/or the value of the signal (Sn+k2·M2n);
detecting a movement from the result of at least one comparison made of the signals Sn, (M3n−M2n) and (M4n+M2n);
where n: natural whole number;
and k1 and k2: non-null positive real numbers.

2. The method according to claim 1, wherein the first mean M1n is obtained from the following equation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        1
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          1
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            1
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          1
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            1
          
        
        ⁢
        
          S
          n
        
      
    
  
  ,



where M1−1=0, and
1/N1: non-null positive real number.

3. The method according to claim 1, wherein the first mean M1n is obtained at least by the following calculation steps:
M10=S0;
and for n>0:
M1n=M1n-1+c1 when M1n-1<Sn;
M1n=M1n-1−c1 when M1n-1>Sn;
where c1: non-null positive real number.

4. The method according to claim 1, wherein the second mean M2n is obtained from the following equation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        2
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          2
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            2
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          2
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            2
          
        
        ⁢
        
          Δ
          n
        
      
    
  
  ,



where M2−1=0, and
1/N2: non-null positive real number.

5. The method according to claim 1, wherein the second mean M2n is obtained at least by the following calculation steps:
M20=Δ0;
and for n>0:
M2n=M2n-1+c2 when M2n-1<Δn;
M2n=M2n-1−c2 when M2n-1>Δn;
where c2: non-null positive real number.

6. The method according to claim 1, wherein the third mean M3n is obtained from the following equation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        3
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          3
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            3
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          3
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            3
          
        
        ⁢
        
          (
          
            
              S
              n
            
            -
            
              
                
                  k
                  1
                
                ·
                M
              
              ⁢
              
                
              
              ⁢
              
                2
                n
              
            
          
          )
        
      
    
  
  ,



where M3−1=0, and
1/N3: non-null positive real number.

7. The method according to claim 6, wherein the value of 1/N3 is adjusted in function of the value |M3n−Sn|:
when M3n<Sn, the value of 1/N3 verify the relation: period (Sn)/ln(1−1/N3)−1>1 s;
when M3n>Sn, the value of 1/N3 verify the relation: period (Sn)/ln(1−1/N3)−1<0.5 s.

8. The method according to claim 1, wherein the third mean M3n is obtained at least by the following calculation steps:
M30=(S0−k1·M20)
and for n>0:
M3n=M3n-1+c3 when M3n-1<(Sn−k1·M2n);
M3n=M3n-1−c3 when M3n-1>(Sn−k1·M2n)
where c3: non-null positive real number.

9. The method according to claim 1, wherein the fourth mean M4n is obtained from the following equation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        4
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          4
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            4
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          4
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            4
          
        
        ⁢
        
          (
          
            
              S
              n
            
            +
            
              
                
                  k
                  2
                
                ·
                M
              
              ⁢
              
                
              
              ⁢
              
                2
                n
              
            
          
          )
        
      
    
  
  ,



where M4−1=0, and
1/N4: non-null positive real number.

10. The method according to claim 9, wherein the value of 1/N4 is adjusted in function of the value |M4n−Sn|:
when M4n<Sn, the value of 1/N4 verify the relation: period (Sn)/ln(1−1/N4)−1<0.5 s;
when M4n>Sn, the value of 1/N4 verify the relation: period (Sn)/ln(1−1/N4)−1>1 s.

11. The method according to claim 1, wherein the fourth mean M4n is obtained at least by the following calculation steps:
M40=(S0+k2·M20)
and for n>0:
M4n=M4n-1+c4 when M4n-1<(Sn+k2·M2n);
M4n=M4n-1−c4 when M4n-1>(Sn+k2·M2n);
where c4: non-null positive real number.

12. The method according to claim 1, wherein the value of the second mean M2n is greater than a second non-null minimum threshold value SM2n.

13. The method according to claim 12, the value of SM2n is comprised between approximately 1/250×the dynamic of the signal Sn and 1/25×the dynamic of the signal Sn.

14. The method according to claim 1, wherein said comparison is made of the signals Sn and (M3n−M2n) and of the signals Sn and (M4n+M2n), wherein a movement is considered as detected when Sn>(M4n+M2n) or Sn<(M3n−M2n).

15. The method according to claim 1, further comprising, between the calculation step of the fourth mean M4n of the signal (Sn+M2n) and the comparison step, a calculation step of a fifth mean M5n of the signal Sn in function of a previous value M5n-1 and/or the value of the signal Sn, wherein said comparison is made on the signals M5n and (M3n−M2n) and on the signals M5n and (M4n+M2n), wherein a movement is considered as detected when M5n>(M4n+M2n) or M5n<(M3n−M2n).

16. The method according to claim 15, wherein the fifth mean M5n is obtained from the following equation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        5
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          5
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            5
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          5
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            5
          
        
        ⁢
        
          S
          n
        
      
    
  
  ,



where M5−1=0, and
1/N5: non-null positive real number.

17. A movement detection device at least comprising:
means of calculating a first mean M1n of a signal Sn designed to be supplied by at least one pixel of a pixel matrix and corresponding to a n-th captured image, in function of the value of the signal Sn and/or a previous value M1n-1;
means of calculating a signal Δn=|M1n−Sn|;
means of calculating a second mean M2n of the signal Δn in function of a previous value M2n-1 and/or the value of the signal Δn;
means of calculating a third mean M3n of the signal (Sn−k1·M2n) in function of a previous value M3n-1 and/or the value of the signal (Sn−k1·M2n);
means of calculating a fourth mean M4n of the signal (Sn+k2·M2n) in function of a previous value M4n-1 and/or the value of the signal (Sn+k2·M2n);
wherein a movement is detected from the result obtained at the output of a means of comparison which make at least one comparison of the signals Sn, (M3n−M2n) and (M4n+M2n)
where n: natural whole number;
and k1 and k2: non-null positive real numbers.

18. The device according to claim 17, wherein the means of calculating the first mean M1n carry out at least the following operation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        1
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          1
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            1
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          1
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            1
          
        
        ⁢
        
          S
          n
        
      
    
  
  ,



where M1−1=0, and
1/N1: non-null positive real number.

19. The device according to claim 17, wherein the means of calculating the first mean M1n at least comprise:
means of initialising the value of M10 to the value of S0;
means of comparing the value of the signal M1n-1 and the value of the signal Sn;
means of incrementing and decrementing the value of M1n by a constant c1, where c1: non-null positive real number.

20. The device according to claim 17, wherein the means of calculating the second mean M2n carry out at least the following operation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        2
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          2
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            2
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          2
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            2
          
        
        ⁢
        
          Δ
          n
        
      
    
  
  ,



where M2−1=0, and
1/N2: non-null positive real number.

21. The device according to claim 17, wherein the means of calculating the second mean M2n at least comprise:
means of initialising the value of M20 to the value of Δ0;
means of comparing the value of the signal M2n-1 and the value of the signal Δn;
means of incrementing and decrementing the value of M2n by a constant c2, where c2: non-null positive real number.

22. The device according to claim 17, wherein the means of calculating the third mean M3n carry out at least the following operation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        3
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          3
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            3
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          3
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            3
          
        
        ⁢
        
          (
          
            
              S
              n
            
            -
            
              
                
                  k
                  1
                
                ·
                M
              
              ⁢
              
                
              
              ⁢
              
                2
                n
              
            
          
          )
        
      
    
  
  ,



where M3−1=0, and
1/N3: positive real number.

23. The device according to claim 22, wherein the means of calculating the third mean M3n further comprise means of adjusting the value of 1/N3 in function of the value |M3n−Sn| such that:
when M3n<Sn, the value of 1/N3 verify the relation: period (Sn)/ln(1−1/N3)−1>1 s;
when M3n>Sn, the value of 1/N3 verify the relation: period (Sn)/ln(1−1/N3)−1<0.5 s.

24. The device according to claim 17, wherein the means of calculating the third mean M3n at least comprise:
means of initialising the value of M30 to the value of (S0−k1·M20);
means of comparing the value of the signal M3n-1 and the value of the signal (Sn−k1·M2n);
means of incrementing and decrementing the value of M3n by a constant c3, where c3: non-null positive real number.

25. The device according to claim 17, wherein the means of calculating the fourth mean M4n carry out at least the following operation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        4
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          4
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            4
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          4
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            4
          
        
        ⁢
        
          (
          
            
              S
              n
            
            +
            
              
                
                  k
                  2
                
                ·
                M
              
              ⁢
              
                
              
              ⁢
              
                2
                n
              
            
          
          )
        
      
    
  
  ,



where M4−1=0, and
1/N4: non-null positive real number.

26. The device according to claim 25, wherein the means of calculating the fourth mean M4n further comprise means of adjusting the value of 1/N4 in function of the value |M4n−Sn| such that:
when M4n<Sn, the value of 1/N4 verify the relation: period (Sn)/ln(1−1/N4)−1<0.5 s;
when M4n>Sn, the value of 1/N4 verify the relation: period (Sn)/ln(1−1/N4)−1>1 s.

27. The device according to claim 17, wherein the means of calculating the fourth mean M4n at least comprise:
means of initialising the value of M40 to the value of (S0+k2·M20);
means of comparing the value of the signal M4n-1 and the value of the signal (Sn+k2·M2n);
means of incrementing and decrementing the value of M4n by a constant c4, where c4: non-null positive real number.

28. The device according to claim 17, wherein the means of comparison are means of comparing the signals Sn and (M3n−M2n) and means of comparing the signals Sn and (M4n+M2n), wherein a movement is considered as detected when Sn>(M4n+M2n) or Sn<(M3n−M2n).

29. The device according to claim 17, further comprising means of calculating a fifth mean M5n of the signal Sn in function of a previous value M5n-1 and/or the value of the signal Sn, wherein the means of comparison are means of comparing the signals M5n and (M3n−M2n) and means of comparing the signals M5n and (M4n+M2n), wherein a movement is considered as detected when M5n>(M4n+M2n) or M5n<(M3n−M2n).

30. The device according to claim 29, wherein the means of calculating the fifth mean M5n carry out at least the following operation:




  
    
      M
      ⁢
      
        
      
      ⁢
      
        5
        n
      
    
    =
    
      
        M
        ⁢
        
          
        
        ⁢
        
          5
          
            n
            -
            1
          
        
      
      -
      
        
          1
          
            N
            5
          
        
        ⁢
        M
        ⁢
        
          
        
        ⁢
        
          5
          
            n
            -
            1
          
        
      
      +
      
        
          1
          
            N
            5
          
        
        ⁢
        
          S
          n
        
      
    
  
  ,



where M5−1=0, and
1/N5: non-null positive real number.

31. An image capture device comprising at least one pixel matrix and one movement detection device according to claim 17.