# Patent

## Title
Method and apparatus for robust velocity prediction

## Abstract

One embodiment relates to a method of velocity prediction using a sensor array. Differential signals are obtained from the sensor array at multiple frame rates. A first velocity prediction is determined from the differential signals at a first frame rate, and a second velocity prediction is determined from the differential signals at a second frame rate. Another embodiment relates to an optical navigation apparatus which includes a sensor array, differential circuitry, driver circuitry, and a signal processing device. The driver circuitry is configured to drive the sensor array and differential circuitry so as to operate at multiple frame rates in an interleaved manner. The signal processing device is configured to compute a first velocity prediction from the differential signals at a first frame rate and to compute a second velocity prediction from the differential signals at a second frame rate. Other embodiments and features are also disclosed.

## Background
TECHNICAL FIELD
The present disclosure relates generally to velocity prediction. More particularly, the present disclosure relates to velocity prediction by optical navigation sensors.
BACKGROUND OF THE INVENTION
Velocity prediction may be performed by various device including two-dimensional (2D) optical navigation sensors and other apparatus. Two-dimensional optical navigation sensors include, for example, optical mouse devices for user interface applications. Such optical mouse devices are commonly employed with personal computers.
It is highly desirable to improve methods and apparatus for velocity prediction. It is particularly desirable to improve velocity prediction by two-dimensional optical navigation sensors, such as optical mouse devices.
SUMMARY
One embodiment relates to a method of velocity prediction using a sensor array. Differential signals are obtained from the sensor array at multiple frame rates. A first velocity prediction is determined from the differential signals at a first frame rate, and a second velocity prediction is determined from the differential signals at a second frame rate.
Another embodiment relates to a method of velocity prediction using a two-dimensional (2D) sensor array. Differential signals from the 2D sensor array at multiple frame rates. First and third velocity predictions are determined from the differential signals at a first frame rate, wherein the first and third velocity predictions measure velocity along two orthogonal axes. Second and fourth velocity predictions are determined from the differential signals at a second frame rate, wherein the second and fourth velocity predictions also measure velocity along the two orthogonal axes.
Another embodiment relates to an optical navigation apparatus which includes a sensor array, differential circuitry, driver circuitry, and a signal processing device. The driver circuitry is configured to drive the sensor array and differential circuitry so as to operate at multiple frame rates in an interleaved manner. The signal processing device is configured to compute a first velocity prediction from the differential signals at a first frame rate and to compute a second velocity prediction from the differential signals at a second frame rate.
Other embodiments and features are also disclosed.



BRIEF DESCRIPTION OF THE DRAWINGS
These and various other features and advantages of the present invention may be apparent upon reading of the following detailed description in conjunction with the accompanying drawings and the appended claims provided below.
FIG. 1 is a schematic diagram of a two-dimensional comb array for sensing 2D motion in accordance with an embodiment of the invention.
FIG. 2 is a plot showing an example of unwrapping phase angles from inverse tangent calculations so as to determine a velocity predictor.
FIG. 3 shows an example timing diagram which switches between a first state to collect signal data for the even frames and a second state to collect signal data for the odd frames in accordance with an embodiment of the invention.
FIG. 4 is a schematic diagram showing an apparatus for robust velocity prediction in accordance with an embodiment of the invention.



DETAILED DESCRIPTION
2D Comb Array for Tracking 2D Motion
FIG. 1 is a schematic diagram depicting an example 2D comb array 102 for detecting 2D motion in accordance with an embodiment of the invention. Such a 2D comb array may be utilized, for example, in an optical navigation sensor, such as an optical mouse device.
The example 2D array 102 is made up of 36 sub-arrays 104 organized in a 6-by-6 matrix. Of course, the array 102 may be of various matrix sizes, the 6×6 matrix size being shown here for purpose of discussion.
An expanded view of one such sub-array 104 is shown on the left side of the figure. In this example with M=4, each sub-array 104 comprises 16 detector elements or photo-detector cells organized in a 4-by-4 matrix. In other implementations, other values of M may be used to create the sub-arrays.
The 16 detector elements in each sub-array 104 are each identified as being a member of one of eight groups of elements. The group number associated with each detector element of each sub-array 104 is shown by the number (1, 2, 3, 4, 5, 6, 7, or 8) labeling the element in the expanded view. The signals from each group are electrically ganged or wired together for the entire array 102. The resultant group signals (numbered 1 through 8) are output from the array 102 (as shown on the right side of the figure).
Differential circuitry 106 is used to generate differential signals from pairs of the group signals. For this example with M=4, there are four differential outputs from the array (CC, CS, SC, and SS). A first differential signal CC is generated by the difference of signals 1 and 2. A second differential signal SC is generated by the difference of signals 3 and 4. A third differential signal CS is generated by the difference of signals 5 and 6. A fourth differential signal SS is generated by the difference of signals 7 and 8.
The spatial frequency distribution of the image captured on the detector array is approximately centered on the spatial frequency of the detector array. Through the processing of the four quasi-sinusoidal outputs (CC, CS, SC, and SS) representing separate in-phase and quadrature signals for motion along two orthogonal axes, the 2D movement of the surface relative to the detector array may be tracked.
Velocity Prediction Using 2D Comb Array
In accordance with an embodiment of the invention, at each sample frame, phase angle values φx and φy, and radius values Rx and Ry, may be computed according to the following equations.




  
    
      
        
          φ
          x
        
        =
        
          
            tan
            
              -
              1
            
          
          ⁡
          
            (
            
              
                CS
                +
                SC
              
              
                CC
                -
                SS
              
            
            )
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          1
        
        )
      
    
  
  
    
      
        
          R
          x
        
        =
        
          
            
              
                (
                
                  CC
                  -
                  SS
                
                )
              
              2
            
            +
            
              
                (
                
                  CS
                  +
                  SC
                
                )
              
              2
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          2
        
        )
      
    
  
  
    
      
        
          φ
          y
        
        =
        
          
            tan
            
              -
              1
            
          
          ⁡
          
            (
            
              
                CS
                -
                SC
              
              
                CC
                +
                SS
              
            
            )
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          3
        
        )
      
    
  
  
    
      
        
          R
          x
        
        =
        
          
            
              
                (
                
                  CC
                  +
                  SS
                
                )
              
              2
            
            +
            
              
                (
                
                  CS
                  -
                  SC
                
                )
              
              2
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          4
        
        )
      
    
  




Rx and Ry indicate the contrast of the detected quasi-sinusoidal signals. Changes in the phase angles φx and φy from a previous sample frame are indicative of motion along the two orthogonal axes (the x and y axes).

The changes in phase angles from a previous sample frame may be denoted by Δφx and Δφy and may be computed according to the following equations.




  
    
      
        
          Δφ
          x
        
        =
        
          
            
              φ
              
                x
                ,
                i1
              
            
            -
            
              
                φ
                
                  x
                  ,
                  
                    i
                    -
                    1
                  
                
              
              ⁢
              
                
              
              ⁢
              where
              ⁢
              
                
              
              ⁢
              
                φ
                
                  x
                  ,
                  i
                
              
            
          
          =
          
            
              tan
              
                -
                1
              
            
            ⁡
            
              (
              
                
                  
                    CS
                    i
                  
                  +
                  
                    SC
                    i
                  
                
                
                  
                    CC
                    i
                  
                  -
                  
                    SS
                    i
                  
                
              
              )
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          5
        
        )
      
    
  
  
    
      
        
          Δφ
          y
        
        =
        
          
            
              φ
              
                y
                ,
                i1
              
            
            -
            
              
                φ
                
                  y
                  ,
                  
                    i
                    -
                    1
                  
                
              
              ⁢
              
                
              
              ⁢
              where
              ⁢
              
                
              
              ⁢
              
                φ
                
                  y
                  ,
                  i
                
              
            
          
          =
          
            
              tan
              
                -
                1
              
            
            ⁡
            
              (
              
                
                  
                    CS
                    i
                  
                  -
                  
                    SC
                    i
                  
                
                
                  
                    CC
                    i
                  
                  +
                  
                    SS
                    i
                  
                
              
              )
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          6
        
        )
      
    
  




In the above equations, the current frame is denoted by i, such that phase angles for a current frame are denoted by the subscript i, and the phase angles for the previous frame are denoted by the subscript i−1.

While Δφx and Δφy are indicative of motion along the x and y axes, they do not completely reflect the actual two-dimensional motion. This is because the values of Δφx and Δφy are restricted to the range from −π to +π due to the inverse tangent function. In other words, the values of Δφx and Δφy are “wrapped” in the range [−π, +π].
Consider the functions ΔΦx and ΔΦy to be “unwrapped” versions of Δφx and ΔΦy, respectively. Hence, Δφx is a modulo function of ΔΦx, and Δφy is a modulo function of ΔΦy, where the values of Δφx and Δφy each “wraps” within the range [−π+ρ]. ΔΦx and ΔΦy are indicative of the actual (unwrapped) motion of the sensor relative to the surface.
Since Δφx and Δφy are computed from the differential signals output by the 2D comb array, they may be “unwrapped” to determine the functions ΔΦx and ΔΦy. Such unwrapping of an example function Δφ(t) to generate the corresponding function Φ(t) is illustrated in FIG. 2.
In FIG. 2, ΔΦ(t) is generated assuming an initial condition constraint of ΔΦ(0)=0. In other words, at t=0, it is assumed that there is no relative motion (i.e. the sensor is at rest relative to the surface).
With this initial assumption, the “actual velocities” ΔΦx and ΔΦy may be computed, for example, by tracking the average velocities over the past K frames (where K>2) and assuming that the next velocity with be within +/−π of the average velocity. This computation may be implemented according to the following equations.




  
    
      
        
          ΔΦ
          x
        
        =
        
          
            Δφ
            x
          
          -
          
            2
            ⁢
            π
            ×
            
              INTEGER
              ⁡
              
                (
                
                  
                    
                      Δφ
                      x
                    
                    -
                    
                      〈
                      
                        ΔΦ
                        x
                      
                      〉
                    
                    +
                    π
                  
                  
                    2
                    ⁢
                    π
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          7
        
        )
      
    
  
  
    
      
        
          ΔΦ
          y
        
        =
        
          
            Δφ
            y
          
          -
          
            2
            ⁢
            π
            ×
            
              INTEGER
              ⁡
              
                (
                
                  
                    
                      Δφ
                      y
                    
                    -
                    
                      〈
                      
                        ΔΦ
                        y
                      
                      〉
                    
                    +
                    π
                  
                  
                    2
                    ⁢
                    π
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          8
        
        )
      
    
  




In the above equations, the INTEGER function takes the largest integer value that is not greater than its argument. <ΔΦx> and <ΔΦy> are the average velocities over the past K frames (K>2) and may be computed according to the following equations.





  
    
      
        
          〈
          
            ΔΦ
            x
          
          〉
        
        =
        
          
            1
            K
          
          ⁢
          
            
              ∑
              
                j
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              ΔΦ
              
                x
                ,
                
                  (
                  
                    i
                    -
                    j
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          9
        
        )
      
    
  
  
    
      
        
          〈
          
            ΔΦ
            y
          
          〉
        
        =
        
          
            1
            K
          
          ⁢
          
            
              ∑
              
                j
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              ΔΦ
              
                y
                ,
                
                  (
                  
                    i
                    -
                    j
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          10
        
        )
      
    
  




<ΔΦx> and <ΔΦy> are utilized within the INTEGER functions so as to track the number of “2π” rotations that have occurred within a frame.

The above-described algorithm needs to assume that the initial velocity is known in order to work. Specifically, the algorithm assumes that the initial velocity is close to zero such that the sensor is at rest with respect to the surface.
Problems and Drawbacks with the Above-Described Technique
The above-described technique has substantial problems and drawbacks.
First, the above-described technique does not determine an absolute velocity, and so it relies on the assumption that the sensor begins at some known velocity (typically, at rest). However, this assumption may be violated. In the example of an optical mouse device, the mouse device may be brought close to the tracking surface at a non-zero initial velocity. In other words, the user's hand may be moving as he/she brings the mouse device to the tracking surface.
If the initial condition assumed is in error, then the velocity predictor technique described above would not detect the error, and furthermore may not readily recover from the error.
For simplicity of discussion, consider the motion in one of the two dimensions. For example, if the motion starts from rest (ΔΦ=0π/frame), then accelerates to a speed of ΔΦ=3π/frame, then decelerates back to rest, the velocity predictor technique described above would work as intended.
However, if the motion starts at 2π/frame, then accelerates to 3π/frame, and then decelerates to rest, then the velocity predictor technique described above would give erroneous results. The velocity prediction error occurs because the initial 2π/frame speed results in a same result from the inverse tangent function as an initial speed of zero (or 4π/frame for that matter). As such, the velocity predictor algorithm erroneously assumes that its starting condition is at rest (ΔΦ=0π/frame). When the motion accelerates to 3π/frame, the velocity predictor algorithm give an erroneous speed of only π/frame. When the motion decelerates to a resting state (zero velocity), the velocity predictor algorithm give an erroneous speed of −2π/frame. In other words, the velocity predictor algorithm continues to report a motion with a negative 2π/frame velocity offset from the actual motion.
Other problems may occur due to under-sampling of the motion with a frame rate which is too low compared with the motion. For example, the determination of the average velocities <ΔΦx> and <ΔΦy> may be erroneous due to such under-sampling.
Robust Velocity Prediction Using 2D Comb Array
The present disclosure provides methods and apparatus to achieve a robust velocity prediction using a 2D comb array.
In accordance with an embodiment of the invention, a robust velocity prediction may be achieved by utilizing a technique where the sensor alternates between two frame rates. The two frame rates are such that they have a non-integer ratio. In other words, neither frame rate is an integer multiple of the other.
This technique is now discussed in relation to a specific example. In this example, consider the first frame rate for even frames to be l/Δt0 and the second frame rate for odd frames to be 1/Δt1. Further, consider that Δt1/Δt0=r, where r is a non-integer number. The non-integer ratio may be selected depending on the specific application.
Consider the problematic example discussed above where the motion along one dimension starts at 2π/frame, then accelerates to 3π/frame, and then decelerates to rest or 0π/frame. As previously discussed, the single frame rate technique described above would give an erroneous velocity prediction.
Now, consider application of the two-frame-rate technique, where r=0.9, for example. When the motion in the one dimension is actually at rest or 0π/frame, the base sample rate would measure a velocity of −2π/frame. However, the slower sample rate which is 90% of the base rate would measure a velocity of −1.87π/frame. This difference between the two velocity measurements indicates that the velocity derived from the base sample rate is incorrect. Correcting the velocity prediction then becomes a matter of determining the nearest velocity that provides a ΔΦ=0π/frame measurement at both frame rates and predicted velocity accordingly. In this case, 0π/frame is the nearest velocity, followed by +/−20π/frame. Hence, the true velocity is determined to be 0π/frame. Note that, in this example, the technique would break down at +/−20π/frame for an initial condition. So, preferably such an initial condition would be beyond the conditions expected to be encountered for the particular application. The following is a more detailed discussion of this technique.
FIG. 3 shows an example timing diagram which switches between a first state (State 0) to collect signal data for the even frames and a second state (State 1) to collect signal data for the odd frames in accordance with an embodiment of the invention. As seen, each even frame is collected over a time period of Δt0, and each odd frame is collected over a time period of Δt1. In others, frames of width Δt0, and frames of width Δt1 are interleaved.
The differential signals obtained from the comb array during even frames (State 0) are utilized to compute a first velocity estimation V0, while the differential signals obtained from the comb array during odd frames (State 1) are utilized to compute a second velocity estimation V1. These computations may be implemented according to the following equations.




  
    
      
        
          V
          0
        
        =
        
          
            
              
                〈
                ΔΦ
                〉
              
              0
            
            
              Δ
              ⁢
              
                
              
              ⁢
              
                t
                0
              
            
          
          =
          
            
              2
              
                K
                ×
                Δ
                ⁢
                
                  
                
                ⁢
                
                  t
                  0
                
              
            
            ⁢
            
              
                ∑
                
                  j
                  =
                  1
                
                
                  K
                  /
                  2
                
              
              ⁢
              
                
              
              ⁢
              
                ΔΦ
                
                  (
                  
                    i
                    -
                    
                      2
                      ×
                      j
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          11
        
        )
      
    
  
  
    
      
        
          
            ΔΦ
            
              (
              
                i
                -
                
                  2
                  ×
                  j
                
              
              )
            
          
          ⁢
          
            Δφ
            
              (
              
                i
                -
                
                  2
                  ×
                  j
                
              
              )
            
          
        
        -
        
          2
          ⁢
          π
          ×
          
            INTEGER
            ⁡
            
              (
              
                
                  
                    Δφ
                    
                      (
                      
                        i
                        -
                        
                          2
                          ×
                          j
                        
                      
                      )
                    
                  
                  -
                  
                    
                      〈
                      ΔΦ
                      〉
                    
                    0
                  
                  +
                  π
                
                
                  2
                  ⁢
                  π
                
              
              )
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          12
        
        )
      
    
  
  
    
      
        
          V
          1
        
        =
        
          
            
              
                〈
                ΔΦ
                〉
              
              1
            
            
              Δ
              ⁢
              
                
              
              ⁢
              
                t
                1
              
            
          
          =
          
            
              2
              
                K
                ×
                Δ
                ⁢
                
                  
                
                ⁢
                
                  t
                  1
                
              
            
            ⁢
            
              
                ∑
                
                  j
                  =
                  1
                
                
                  K
                  /
                  2
                
              
              ⁢
              
                
              
              ⁢
              
                ΔΦ
                
                  (
                  
                    i
                    -
                    
                      2
                      ×
                      j
                    
                    +
                    1
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          13
        
        )
      
    
  
  
    
      
        
          ΔΦ
          
            
              (
              
                i
                -
                
                  2
                  ×
                  j
                
              
              )
            
            +
            1
          
        
        =
        
          
            Δφ
            
              (
              
                i
                -
                
                  2
                  ×
                  j
                
                +
                1
              
              )
            
          
          -
          
            2
            ⁢
            π
            ×
            
              INTEGER
              ⁡
              
                (
                
                  
                    
                      Δφ
                      
                        (
                        
                          i
                          -
                          
                            2
                            ×
                            j
                          
                          +
                          1
                        
                        )
                      
                    
                    -
                    
                      
                        〈
                        ΔΦ
                        〉
                      
                      1
                    
                    +
                    π
                  
                  
                    2
                    ⁢
                    π
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          14
        
        )
      
    
  




In the above equations, the current frame is denoted by i, and K is an even number and is greater than 2. <ΔΦ>0 is the average change in phase over the prior K/2 even frames, and <ΔΦ>1 is the average change in phase over the prior K/2 odd frames.

Ideally, V0 should equal to V1, such that the following equation is satisfied.




  
    
      
        
          
            
              〈
              ΔΦ
              〉
            
            1
          
          
            
              〈
              ΔΦ
              〉
            
            0
          
        
        =
        
          
            
              Δ
              ⁢
              
                
              
              ⁢
              
                t
                1
              
            
            
              Δ
              ⁢
              
                
              
              ⁢
              
                t
                0
              
            
          
          =
          r
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          15
        
        )
      
    
  




However, the above ideal equation (Eqn 15) is not necessarily satisfied because there may be extra or missing 2π “rotations” in <ΔΦ>0 as a result of the initial condition assumption and/or under-sampling or other reasons.

In accordance with a first embodiment, it is assumed that a faster frame rate yields a more accurate velocity estimation. In a first case, the second frame rate (for odd frames) is set to be faster than the first frame rate (for even frames). In other words, Δt1<Δt0, or 0≦r≦1. In this case, the velocity estimation based on the odd frames (i.e. <ΔΦ>) is assumed to be the more accurate velocity estimation. In a second case, the first frame rate (for even frames) is set to be faster than the second frame rate (for odd frames). In other words, Δt0<Δt1, or r>1. In this case, the velocity estimation based on the even frames (i.e. <ΔΦ>0) is assumed to be the more accurate velocity estimation.
The extra or missing 2π “rotations” in the velocity estimation based on the slower frame rate may be determined and accounted for as follows. For example, consider the case where the second frame rate (for odd frames) is set to be faster than the first frame rate (for even frames). A new or corrected value of <ΔΦ>0 may be denoted (<ΔΦ>0)new. (<ΔΦ>0)new satisfies the following equations.

(ΔΦ0)new=ΔΦ0+2π×N0  (Eqn 16)

(ΔΦ0)new×r=ΔΦ  (Eqn 17)

where N0 is the number of missing 2π rotations in <ΔΦ>0 in comparison to <ΔΦ>1. Substituting for N0 in the first equation above (Eqn 16) gives the following equation for computing (<ΔΦ>0)new.





  
    
      
        
          
            (
            
              
                〈
                ΔΦ
                〉
              
              0
            
            )
          
          new
        
        =
        
          
            
              〈
              ΔΦ
              〉
            
            0
          
          +
          
            2
            ⁢
            π
            ×
            
              INTEGER
              ⁡
              
                (
                
                  
                    
                      
                        〈
                        ΔΦ
                        〉
                      
                      1
                    
                    -
                    
                      
                        
                          〈
                          ΔΦ
                          〉
                        
                        0
                      
                      ×
                      r
                    
                  
                  
                    2
                    ⁢
                    π
                    ×
                    r
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          18
        
        )
      
    
  




In the above equation (Eqn 18), the INTEGER function takes the largest integer value that is not greater than its argument. As calculated according to the above equation (Eqn 18), the function (<ΔΦ>0)new becomes a velocity estimation which is corrected for extra or missing phase rotations.

In accordance with a second embodiment, the absolute values of the two velocity estimations (under the two different frame rates) are compared. This technique then assumes the one with the smaller absolute value (i.e. the one with the slower speed) is the accurate velocity estimation. This second approach may be particularly advantageous in overcoming the “run-away” problem in comb array based position tracking due to over-estimation of the velocity predictor with extra phase rotations.
If the velocity estimation based on odd frames (i.e. <ΔΦ>1) has the smaller absolute value, then <ΔΦ>1 is assumed to be accurate and the above equation (Eqn 18) may be utilized to correct the velocity estimation for the even frames. On the other hand, if the velocity estimation based on even frames (i.e. <ΔΦ>0) has the smaller absolute value, then <ΔΦ>0 is assumed to be accurate and the below equation (Eqn 19) may be utilized to correct the velocity estimation for the odd frames.




  
    
      
        
          
            (
            
              
                〈
                ΔΦ
                〉
              
              1
            
            )
          
          new
        
        =
        
          
            
              〈
              ΔΦ
              〉
            
            1
          
          +
          
            2
            ⁢
            π
            ×
            
              INTEGER
              ⁡
              
                (
                
                  
                    
                      
                        
                          〈
                          ΔΦ
                          〉
                        
                        0
                      
                      ×
                      r
                    
                    -
                    
                      
                        〈
                        ΔΦ
                        〉
                      
                      1
                    
                  
                  
                    2
                    ⁢
                    π
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          19
        
        )
      
    
  




where the new or corrected value of <ΔΦ>1 is denoted (<ΔΦ>1)new.

In accordance with a third embodiment, a single velocity predictor value is estimated by mixing the delta-phases from both frame rates into a single average velocity V. This technique assumes that some of the initial condition constraint errors or other estimation errors may be removed or mitigated by this mixed averaging. The average velocity V may be determined using the equations below.




  
    
      
        V
        =
        
          
            1
            K
          
          ⁢
          
            
              ∑
              
                j
                =
                1
              
              
                K
                /
                2
              
            
            ⁢
            
              
            
            ⁢
            
              [
              
                
                  
                    ΔΦ
                    
                      (
                      
                        i
                        -
                        
                          2
                          ×
                          j
                        
                      
                      )
                    
                  
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      t
                      0
                    
                  
                
                +
                
                  
                    ΔΦ
                    
                      (
                      
                        j
                        -
                        
                          2
                          ×
                          j
                        
                        +
                        1
                      
                      )
                    
                  
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      t
                      1
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        
          Eqn
          ⁢
          
            
          
          ⁢
          20
        
        )
      
    
  




The averaging of this approach may be implemented with less logic than the prior two approaches.

Note that for 2D position tracking, one of the above three approaches would be performed along each of the two orthogonal directions. Also, note that while the above descriptions discusses velocities in terms of changes in phase (i.e. phase velocities), these phase velocities correspond to linear velocities, where the conversion between the two depend upon the dimensions and geometries of the specific apparatus.
In the above three approaches, the ratio r of the two frame rates is preferably selected to allow the entire range of anticipated velocities to be discriminated without requiring resolution beyond the capabilities of the measurement system.
While the technique described above uses two frame rates, other embodiments may utilize three or more frame rates. For example, under the third approach, the delta-phases from the three or more frame rates may be averaged to compute the single average velocity V.
The above-described technique for tracking 2D motion using a 2D comb detector array may be advantageously implemented so as to require less signal processing power than existing techniques. Such lower power consumption may be particularly suitable for power-sensitive applications, such as, for example, wireless mouse devices.
FIG. 4 is a schematic diagram showing an apparatus for robust velocity prediction in accordance with an embodiment of the invention. The apparatus includes a 2D comb detector array 102 and differential circuitry 106 as discussed above in FIG. 1. The apparatus further includes multiple-frame-rate driver circuitry 402 and a signal processing device 404 configured to process the multiple-frame-rate data.
The driver circuitry 402 is configured to provide a timing signals so as to drive the 2D comb detector array 102 and differential circuitry 106 at two (or more) frame rates, where frames at one rate are interleaved with frames from the other rate.
The signal processing device 404 is configured to receive the differential signals at the multiple frame rates from the differential circuitry 106 and to process those signals to generate a robust 2D velocity prediction. The robust 2D velocity prediction may be implemented using the multiple frame rate technique discussed above. In one implementation, the signal processing device 404 may include two (or more) buffers 410, each buffer being configured to hold data associated with a particular frame rate of image acquisition.
The signal processing device 404 may be implemented with hardware signal processing circuitry and/or a combination of hardware and software. In the case of the hardware/software combination, the software comprises processor-executable instructions to perform at least a portion of the technique.
The foregoing description of specific embodiments and examples of the invention have been presented for the purpose of illustration and description, and although the invention has been described and illustrated by certain of the preceding examples, it is not to be construed as being limited thereby. They are not intended to be exhaustive or to limit the invention to the precise forms disclosed, and many modifications, improvements and variations within the scope of the invention are possible in light of the above teaching. It is intended that the scope of the invention encompass the generic area as herein disclosed, and by the claims appended hereto and their equivalents.

## Claims

1. A method of velocity prediction using a sensor array, the method comprising:
obtaining differential signals from the sensor array at multiple frame rates;
determining a first velocity prediction from the differential signals at a first frame rate; and
determining a second velocity prediction from the differential signals at a second frame rate; and
processing the first and second velocity predictions so as to derive a correction for extra or missing phase rotations, wherein a velocity prediction from the faster of the first and second frame rates is utilized to correct a velocity prediction from the slower of the first and second frame rates.

2. The method of claim 1, wherein the sensor array comprises an array of detector cells configured as a comb array.

3. The method of claim 1, wherein the first frame rate is 1/Δt0, the second frame rate is 1/Δt1, and the sensor array obtains the differential signals at the two frame rates by interleaving frames of width Δt0 and frames of width Δt1.

4. The method of claim 3, wherein the ratio of the first frame rate to the second frame rate is a non-integer number.

5. A method of velocity prediction using a sensor array, the method comprising:
obtaining differential signals from the sensor array at multiple frame rates;
determining a first velocity prediction from the differential signals at a first frame rate; and
determining a second velocity prediction from the differential signals at a second frame rate; and
processing the first and second velocity predictions so as to derive a correction for extra or missing phase rotations, wherein a slower velocity prediction of the first and second velocity predictions is utilized to correct a faster velocity prediction of the first and second velocity predictions.

6. The method of claim 1, further comprising:
averaging velocity predictions from the multiple frame rates so as to determine an average velocity prediction.

7. The method of claim 1, further comprising:
a third velocity prediction from the differential signals at the first frame rate, wherein the first and third velocity predictions measure velocity along two orthogonal axes; and
determining a fourth velocity prediction from the differential signals at the second frame rate, wherein the second and fourth velocity predictions also measure velocity along the two orthogonal axes.

8. The method of claim 7, wherein the 2D sensor array comprises an array of detector cells configured as a 2D comb array.

9. The method of claim 7, wherein the first frame rate is 1/Δt0, the second frame rate is 1/Δt1, and the sensor array obtains the differential signals at the two frame rates by interleaving frames of width Δt0 and frames of width Δt1.

10. The method of claim 9, wherein the ratio of the first frame rate to the second frame rate is a non-integer number.

11. The method of claim 7, wherein said processing of the first and second velocity predictions corrects extra or missing phase rotations along a first axis, and further comprising:
processing the third and fourth velocity predictions so as to derive a correction for extra or missing phase rotations along a second axis.

12. The method of claim 7, further comprising:
averaging the first and second velocity predictions so as to determine a first average velocity prediction along a first axis; and
averaging the third and fourth velocity predictions so as to determine a second average velocity prediction along a second axis.

13. An optical navigation apparatus, the apparatus comprising:
a sensor array configured to generate multiple output signals based on optical signals detected by groups of detector elements in the array;
differential circuitry configured to receive the multiple output signals and to generate differential signals;
driver circuitry configured to drive the sensor array and differential circuitry so as to operate at multiple frame rates in an interleaved manner; and
a signal processing device configured to compute a first velocity prediction from the differential signals at a first frame rate and to compute a second velocity prediction from the differential signals at a second frame rate,
wherein the sensor array comprises a two-dimensional comb array,
wherein the signal processing device is further configured to compute a third velocity prediction from the differential signals at a first frame rate and to compute a fourth velocity prediction from the differential signals at a second frame rate, and
wherein the first and second velocity predictions measure motion alone a first axis and the third and fourth velocity predictions measure motion alone a second axis which is orthogonal to the first axis.

14. The apparatus of claim 13, wherein the first frame rate is 1/Δt0, the second frame rate is 1/Δt1, and the sensor array obtains the differential signals at the two frame rates by interleaving frames of width Δt0 and frames of width Δt1.

15. The apparatus of claim 14, wherein the ratio of the first frame rate to the second frame rate is a non-integer number.

16. The apparatus of claim 13, wherein the signal processing device is further configured to process the first and second velocity predictions so as to derive a correction for extra or missing phase rotations.

17. The apparatus of claim 13, wherein the signal processing device is further configured to average velocity predictions from the multiple frame rates so as to determine an average velocity prediction.