Eye tracking system and method

Information

  • Patent Grant
  • 7653213
  • Patent Number
    7,653,213
  • Date Filed
    Friday, September 30, 2005
    19 years ago
  • Date Issued
    Tuesday, January 26, 2010
    15 years ago
Abstract
A method of tracking an expected location of a head in a computerized headtracking environment having a delayed processing requirement for locating a current head position, the method comprising the step of: utilizing previously tracked positions to estimate a likely future tracked position; outputting the likely future tracked position as the expected location of the head. Kalman filtering of the previously tracked positions can be utilized in estimating the likely future tracked position.
Description
FIELD OF THE INVENTION

The present invention relates to a system for accurate prediction of a current eye location and, in particular, discloses a system for head prediction suitable for utilisation in stereoscopic displays.


BACKGROUND OF THE INVENTION

Auto-stereoscopic displays give the observer the visual impression of depth, and are therefore specifically useful for applications in the CAD area, but also have applications in 3D gaming and motion picture entertainment. The impression of depth is achieved by providing the two eyes of the observer with different images which correspond to the view from the respective eye onto the virtual scene. For background information on Autostereoscopic Displays, reference is made to: “Autostereoscopic Displays and Computer Graphics”, by Halle in Computer Graphics, ACM SIGGRAPH, 31(2), May 1997, pp58-62.


Passive auto-stereoscopic displays require the observer to hold their head in a specified position, the sweet spot, where the eyes can observe the correct images. Such systems require the user to keep their head in this specified position during the whole experience and therefore have low market acceptance. When looked at from a position other than the sweet spot, the image looses the impression of depth and becomes inconsistent, resulting in eye strain as the brain attempts to make sense of the images it perceives. This eye strain can generate a feeling of discomfort very quickly which encumbers the market acceptance even more.


Active auto-stereoscopic displays in addition contain a device to track the position of the head and the eyes, typically a camera coupled with IR LED illumination, but other methods such as magnetic or capacitive methods are feasible. Once the position of the eyes relative to the display is known, the display is adjusted to project the two image streams to the respective eye locations. This adjustment can be achieved either by a mechanical device that operates a physical mask which is placed in front of the display or by a liquid crystal mask that blocks the view to the display from certain directions but allows the view from other directions, i.e. the current position of the eyes. Such displays allow the users head to be in a convenient volume in front of the auto-stereoscopic display while the impression of depth is maintained.


Although active auto-stereoscopic displays are much more practicable than passive displays, it has been found that such displays can suffer from the lag introduced by the head tracking system. When moving the head, the time between the actual head motion and the adjustment of the display to the new head position causes an offset sufficiently large to break the impression of depth and the consistency of the images with the previously described problems. This effect is particularly visible with mechanically adjusted displays.


Often applications for active auto-stereoscopic displays specifically use the head position of the observer not only to adjust the display to maintain the impression of depth but also to change the viewpoint of the scene. Such systems actively encourage the observer to move their head to get a view of the scene from different directions. In such applications visual consistency breakdowns during every head motion reduces the usability.


SUMMARY OF THE INVENTION

It is an object of the present invention to provide for a system for real time eye position prediction.


In accordance with a first aspect of the present invention, there is provided a method of tracking an expected location of a head in a computerized headtracking environment having a delayed processing requirement for locating a current head position, the method comprising the step of: utilizing previously tracked positions to estimate a likely future tracked position; outputting the likely future tracked position as the expected location of the head.


Preferably, Kalman filtering of the previously tracked positions can be utilized in estimating the likely future tracked position. The likely future tracked position can be utilized to control an auto-stereoscopic display for the display of images for eyes located at expected positions corresponding to the likely future tracked position.


In accordance with a further aspect of the present invention, there is provided a system for providing an expected location of a head the system comprising: video input means for providing at least one video signal of the head; first processing means for processing the video signal so as to output a substantially continuous series of current head location data; second processing means for processing predetermined one of the current head location data so as to output a predicted future expected location output of the head. The video input means preferably can include stereo video inputs. The second processing means can utilize a Kalman filtering of the current head location data. The system can be interconnected to an auto-stereoscopic display driven by the predicted expected location output of the head.


In accordance with a further aspect of the present invention, there is provided in a camera based face tracking system, a method of predicting the future position of a face, the method comprising the steps of: (a) providing a current prediction of the face position using facial features detected in a previous and current input image frame; (b) deriving a first covariance matrix from the current prediction and a previous covariance matrix; (c) utilizing the current prediction of the face position from step (a) and a Kalman filter to determine a corresponding projected point of the facial feature on the plane of at least one camera; (d) deriving a Jacobain of the projected points in the step (c); (e) deriving a residual covariance of the projected points in the step (c); (f) deriving a suitable filter gain for the Kalman filter; (g) deriving a suitable update coefficients for the first covariance matrix; (h) updating the Kalman filter utilizing the filter gain;


The face tracking system preferably can include multiple cameras observing a user and the steps (a) to (h) are preferably carried out for substantially each camera. Further, the method can also include the step of: (i) determining a corresponding expected eye position from the current state of the Kalman filter.


A noise component can be added to the first covariance matrix. The noise component preferably can include a translational noise component and a rotational noise component. Further, the residual covariance of the step (e) can be utilized to tune response of the Kalman filter.





BRIEF DESCRIPTION OF THE DRAWINGS

Preferred forms of the present invention will now be described by way of example only, with reference to the accompanying drawings in which:



FIG. 1 illustrates schematically a top view of a user using an autostereoscopic display in accordance with the preferred embodiment;



FIG. 2 illustrates schematically a side view of a user using an autostereoscopic display in accordance with the preferred embodiment;



FIG. 3 illustrates the processing chain of the system of the preferred embodiment; and



FIG. 4 illustrates the relationship between user space and image space.



FIG. 5 illustrates the head reference frame which is fixed relative to a head.





DESCRIPTION OF PREFERRED AND OTHER EMBODIMENTS

In the preferred embodiment, there is provided a method for reducing the adjustment lag of auto-stereoscopic displays to thereby improve their usability. Ideally, the method includes the utilisation of a prediction filter that is optimal for the requirements of auto-stereoscopic displays although other methods are possible.


Turning initially to FIG. 1 and FIG. 2, there is illustrated schematically an arrangement of a system for use with the preferred embodiment wherein a user 2 is located in front of an automatic stereoscopic display 3. Two cameras 4, 5 monitor the user and their video feeds are processed to derive a current facial position. The cameras 4, 5 are interconnected to a computer system implementing facial tracking techniques.



FIG. 3 illustrates schematically the subsequent operation of processing chain incorporating the preferred embodiment. The camera feeds e.g. 4,5, are fed to a visual head tracker 7 which tracks a current position of the user's head. The head tracker 7 can be one of many standard types available on the market. The system utilized in the preferred embodiment was that disclosed in International PCT patent application No. PCT/AU01/00249 entitled “Facial Image Processing System” assigned to the present applicant, the contents of which are incorporated herewith. The face tracking system 7 takes an input from the two cameras and derives a current face location 8. The face location is ideally derived in real time. Subsequently, a face location predictor 9 is implemented which takes the face location 8 and outputs a predicted face location 10. This is then fed to the Autostereoscopic display device 3 for use in outputting images to the user.


The functionality of the head location predictor 9 is to predict the position of the eyes of a person looking at the autostereoscopic display in a coordinate system fixed relative to the display. The eye position information is used in turn by the autostereoscopic display to produce different images when seen by the left and right eye of the user and thus create the illusion of depth.


Notation


The index i∈[0,n] is used for numerating the cameras. In an example embodiment the number of cameras n=1 or n=2


The index j∈[0,m] is used for numerating facial features of a user. In initial experiments a variable number of facial features was used with a typical value m=15.


Vectors are typically noted in bold letters, while scalars are usually noted in non-bold letters. 3D vectors are usually expressed in BOLD UPPERCASE while 2D vectors are usually noted in bold lowercase.


When writing a geometric vector, the reference frame (if any) is indicated to the top left of the vector, while the facial feature index (if any) is indicated on the bottom left of the vector. Thus jip represents the 2D projection of the facial feature j observed by camera i in its image plane referential. jiP represents the 3D position of the facial feature j in the reference frame of camera i.


Camera Projective Geometry


A pinhole camera model is used for the projection of 3D Points onto the camera image plane. The model for projection is shown in FIG. 4. A point (20) iP=(iPx,iPy,iPz)T in the reference frame of camera i∈[0,n] projects onto the image plane (21) at a point ip=(ipx, ipy)T in the image plane reference frame of camera i∈[0,n] following the equations












i


p

=




i


o

+

(





f
x



i




P
x



i


P
z



i









f
y



i




P
y



i


P
z



i






)






Equation





1







where io=(iox,ioy)T is the principal point and if=(ifx, ify)T is the focal length of camera i∈[0,n]. In initial experiments, the image size is 640×480 pixels, the principal point is near the center of the image and the focal length is typically around 1800 pixels.


Reference Frames


System Reference Frame


The system reference frame S is fixed relative to the camera(s) and the autostereoscopic display and is shown in FIG. 2 and FIG. 4 with its x-axis 11 horizontal, the y-axis 13 pointing up and the z-axis 12 pointing toward a user.


A point iP expressed in the camera reference frame i is related to a point sP expressed in the system reference frame S with the equation

iP=SiRSP+iTS  Equation 2


Head Reference Frame


The head reference frame H is fixed relative to the head being tracked as shown in FIG. 5. The origin of the head reference frame is placed at the midpoint of the eyeball centers of the left and right eyes, with the x axis 15 aligned on the eyeball centers, the y axis 16 pointing up and the z axis 17 pointing toward the back of the head.


A point HP expressed in the head reference frame H is related to a point SP expressed in the system reference frame S with the equation:

SP=HSRHP+STH  Equation 3


Head Pose


The head pose is defined as the head translation and rotation expressed in the system reference frame and is described by the rotation matrix HSR and the translation vector STH


The head pose rotation HSR is stored using a vector of Euler angles e=(ex,ey,ez)T. If cx=cos(ex), sx=sin(ex), . . . then:















H
S


R

=


(



1


0


0




0



c
x




-

s
x






0



s
x




c
x




)



(




c
y



0



s
y





0


1


0





-

s
y




0



c
y




)



(




c
z




-

s
z




0





s
z




c
z



0




0


0


1



)








=

(





c
y



c
z






-

c
y




s
z





s
y








s
x



s
y



c
z


+


c
x



s
z








c
x



c
z


-


s
x



s
y



c
z







-

s
x




c
y








s
x



s

z
-




c
x



s
y



c
z







c
x



s
y



s
z


+


s
x



c
z







c
x



c
y





)








Equation





4







Head Model


The head model is a collection of 3D point jHM, j∈[0,m] expressed in the head reference frame. Each point represents a facial feature being tracked on the face.


Eye Position in the Head Model


The center of the eyeballs in the head model are noted jHE, j=0,1. The right eyeball center is noted 0HE and the left eyeball center is noted 1HE.


Head Pose Estimation Using Extended Kalman Filtering (EKF)


State of the Kalman Filter


Given








T
H



S

=


(


t
x

,

t
y

,

t
z


)

T


,











t




T
H



S


=


(



t
.

x

,


t
.

y

,


t
.

z


)

T







and
















t



e

=


(



e
.

x

,


e
.

y

,


e
.

z


)

T


,





the state of the Extended Kalman Filter is selected as the position (rotation and translation) and the corresponding velocity of the head expressed in the system reference frame.









x
=


(


e
x

,

e
y

,

e
z

,

t
x

,

t
y

,

t
z

,


e
.

x

,


e
.

y

,


e
.

z

,


t
.

x

,


t
.

y

,


t
.

z


)

T





Equation





5







Dynamics of the Head Motion


The position of the eyes in the system reference frame is predicted by modeling the motion of the head with a set of constant dynamics. In the example embodiment, a constant velocity model is used with the noise being modeled as a piecewise constant acceleration between each measurement. (Similar techniques are outlined in Yaakov Bar-Shalom, Xiao-Rong Li: Estimation and Tracking, Principles, Techniques, and Software, Artech House, 1993, ISBN 0-89006-643-4, at page 267).

xk+1|k=Fxk+Γνk  Equation 6


The transition matrix Fk+1 is:










F

k
+
1


=

(




I

6
×
6






I

6
×
6



T





0



I

6
×
6





)





Equation





7







where I6×6 is the 6×6 identity matrix and T is the sample time, typically 16.66 ms for a 60 Hz measurement frequency.


Γ is the gain multiplying the process noise, with a value fixed at









Γ
=

(





I

6
×
1





T
2

2








I

6
×
1



T




)





Equation





8







where I6×1 is the 6×1 column vector fill with 1, and T is again the sample time.


Initialization of the Kalman Filter when Face is Found:


Upon detecting the face, the state is set to the estimated head pose obtained from an initial face searching algorithm. Many different example algorithms can be used for determining an initial position. In the preferred embodiment the techniques discussed in International PCT patent application No. PCT/AU01/00249 were used to provide an initial head pose estimate, with a null velocity.


The covariance matrix is empirically reset to










P

0



0



=


1
100



(





q
e



I

3
×
3






0

3
×
3





0

3
×
3





0

3
×
3







0

3
×
3






q
t



I

3
×
3






0

3
×
3





0

3
×
3







0

3
×
3





0

3
×
3






q
e



I

3
×
3






0

3
×
3







0

3
×
3





0

3
×
3





0

3
×
3






q
t



I

3
×
3






)






Equation





9







Iteration of the Kalman Filter during Tracking:


1. Prediction of the State


At the beginning of the each new image frame k+1, the state is predicted according to a constant velocity model:

xk+1|k=Fxk  Equation 10

2. Prediction of the Covariance Matrix


The covariance matrix is updated according to dynamics and process noise

Pk+1|k=FPkFT+Qk  Equation 11


Qk represents the process noise and is computed according to a piecewise constant white acceleration model (As for example set out in Yaakov Bar-Shalom, Xiao-Rong Li: Estimation and Tracking, Principles, Techniques, and Software, Artech House, 1993, ISBN 0-89006-643-4, at page 267).













Q
k

=

E


[

Γ








υ
k



(

Γ






υ
k


)


T


]








=

(





I
3




T
4

4



q
e




0




I
3




T
3

2



q
e




0




0




I
3




T
4

4



q
t




0




I
3




T
3

2



q
t








I
3




T
3

2



q
e




0




I
3



T
2



q
e




0




0




I
3




T
3

2



q
t




0




I
3



T
2



q
t





)








Equation





12







In the implementation, we have set the translation process noise was set: qe=0.01 m·s−2 and the rotation process noise qt=0.01 rad·s−2


3. Prediction of the Image Measurements


For a facial feature j observed from camera i, the expected projection jip(xk+1|k) is computed in the image plane according to the predicted state xk+1|k


The 3D point jHM, j∈[0,m]corresponding to the facial feature j is first rotated according to the state x of the Kalman filter into a point jSM

jSM(x)HSR(x)jHM+STH(x)  Equation 1


The point jSM is then expressed in the reference frame of camera i∈[0,n]

jiM(jSM)SiRjSM+iTS  Equation 2


The point jiM is then projected onto the image plane of camera i∈[0,n] into a point jip













j
i


p



(



j
i


M

)


=




i


o

+

(





f
x



i




M
x


j
i


M
z


j
i









f
y



i




M
y


j
i


M
z


j
i






)






Equation





3








4. Computation of the Jacobian of the Measurement Process


The projection can be summarised as:

jip(x)=jip(jiM(jSM(x)))  Equation 4


The Jacobian of the projection is thus

















j
i


p



(
x
)





x


=








j
i


p



(



j
i


M

)








j
i


M











j
i


M



(



j
S


M

)








j
S


M











j
S


M



(
x
)





x









with




Equation





5













j
i


p



(



j
i


M

)








j
i


M



=

(





f
x



i


M
z


j
i




0



-



f
x



i



M
x


j
i




(

M
z


j
i

)

2







0




f
y



i


M
z


j
i





-



f
y



i



M
y


j
i




(

M
z


j
i

)

2






)





Equation





6














j
i


M



(



j
S


M

)








j
S


M



=



S
i


R














j
S


M



(
x
)





x


=

(



0



a
0




a
1



1


0


0


0


0


0


0


0


0





a
2




a
3




a
4



0


1


0


0


0


0


0


0


0





a
5




a
6




a
7



0


0


1


0


0


0


0


0


0



)





Equation





7








a
0

=



-

s
y




c
z



M
x


j
H


+


s
y



s
z



M
y


j
H


+


c
y



M
z


j
H











a
1

=



-

c
y




s
z



M
x


j
H


-


c
y



c
z



M
y


j
H











a
2

=



(



c
x



s
y



c
z


-


s
x



s
z



)



M
x


j
H


-


(



c
x



s
y



s
z


+


s
x



s
z



)



M
y


j
H


-


c
x



c
y



M
z


j
H











a
3

=



s
x



c
y



c
z



M
x


j
H


-


s
x



c
y



s
z



M
y


j
H


+


s
x



s
y



M
z


j
H











a
4

=



(



c
x



c
z


-


s
x



s
y



s
z



)



M
x


j
H


-


(



s
x



s
y



c
z


+


c
x



s
z



)



M
y


j
H











a
5

=



(



s
x



s
y



c
z


+


c
x



s
z



)



M
x


j
H


+


(



c
x



c
z


-


s
x



s
y



s
z



)



M
y


j
H


-


s
x



c
y



M
z


j
H











a
6

=



-

c
x




c
y



c
z



M
x


j
H


+


c
x



c
y



c
z



M
y


j
H


-


c
x



s
y



M
z


j
H











a
7

=



(



c
x



s
y



s
z


+


s
x



c
z



)



M
x


j
H


+


(



c
x



s
y



c
z


-


s
x



s
z



)



M
y


j
H








Equation





8








5. Computation of Residual Covariance


Once the Jacobian of the projection has been computed, the residual covariance can be calculated for the measurement of feature j from camera i










S

k
+
1



j
i

=







j
i


p







(

x


k
+
1

|
k


)





x


k
+
1

|
k






P






k
+
1

|
k









j
i


p








(

x


k
+
1

|
k


)

T





x


k
+
1

|
k





+

V

k
+
1



j
i






Equation





9







where jiVk+1 represents the measurement noise of the facial feature j observed by camera i at frame k+1. We have set empirically










V

k
+
1



j
i

=

α







(










dx
2










dxdy











dxdy






dy
2





)


-
1







Equation





10







where dx and dy represent the image gradient in the image patch used to locate the facial feature j observed by camera i. The coefficient α can be used to tune the responsiveness of the filter (control the balance between the measurements and the process dynamics).


6. Computation of the Filter Gain


The filter gain can then be computed:










W

k
+
1



j
i

=


P


k
+
1

|
k








j
i


p








(

x


k
+
1

|
k


)

T





x


k
+
1

|
k






S

k
+
1


-
1


j
i






Equation





11








7. Update of the State Covariance

Pk+1|k+1=Pk+1|kjiWk+1Sk+1jiWk+1T  Equation 12

8. Update of the State of the Kalman Filter


If we note jizk+1 the measurement of the projection of the facial feature j from camera i taken at frame k+1 (obtained from zero-mean normalized cross-correlation), then the state of the filter can be updated with the equation

xk+1|k+1=xk+1|k+jiWk+1(jizk+1jip(xk+1|k))  Equation 13

9. Loop until all the Measurements Have Been Entered


Assign the new predicted state xk+1|k=xk+1|k+1. Return to step 0 until all the measurements from all the cameras are entered in the Kalman filter. The current implementation first loop on all the features measured by camera 0 (camera A), then loops on all the features measured by camera 1 (camera B).


When all the measurements have been entered, the last state represents the head position and velocity for the image frame k+1.


Prediction of the Eye Position


Given the current state of the Kalman filter x(T0) at time T0, a prediction of the head translation and rotation at time T1=T0+ΔT can be derived from Equation 10 as

x(T1)=FT)x(T0)  Equation 14


with the transition matrix










F






(

Δ





T

)


=

(




I

6
×
6






I

6
×
6



Δ





T





0



I

6
×
6





)





Equation





15







The predicted position of the eyeball center in the system reference frame at time T1 can then be computed from Equation 1

jSE(T1)=jSE(x(T1))=HSR(x(T1))jHE+STH(x(T1))  Equation 16


This completes the cycle of computation for one measurement frame and the predicted eye position jSE(T1) are forwarded to the autostereoscopic display to coincide the emitted image with the actual position of the eyes.


The foregoing describes only preferred forms of the present invention. Modifications, obvious to those skilled in the art can be made thereto without departing from the scope of the invention.

Claims
  • 1. A face tracking system for predicting the future position of a face comprising: multiple cameras for observing a user, each camera carrying out the steps of:(a) providing a current prediction of the face position using facial features detected in a previous and current input image frame;(b) deriving a first covariance matrix from the current prediction and a previous covariance matrix;(c) utilizing said current prediction of the face position from step (a) and a Kalman filter to determine a corresponding projected point of the facial feature on the plane of at least one camera;(d) deriving a Jacobian of the projected points in said step (c);(e) deriving a residual covariance of the projected points in said step (c);(f) deriving a suitable filter gain for said Kalman filter;(g) deriving a suitable update coefficients for said first covariance matrix;(h) updating said Kalman filter utilizing said filter gain.
  • 2. The face tracking system of claim 1 further comprising each camera carrying out the step of: (i) determining a corresponding expected eye position from the current state of the Kalman filter.
  • 3. The face tracking system of claim 1 wherein a noise component is added to said first covariance matrix.
  • 4. The face tracking system of claim 3 wherein said noise component includes a translational noise component and a rotational noise component.
  • 5. The face tracking system of claim 4 wherein said residual covariance of said step (e) is utilized to tune response of the Kalman filter.
  • 6. The face tracking system of claim 1, wherein said face tracking system further comprises an auto-stereoscopic display driven by said prediction of future position of the face.
  • 7. A system for tracking an expected location of a head in a computerized headtracking environment having a delayed processing requirement for locating a head position, the system comprising: a device, such that the system utilizes previously tracked positions to estimate a likely future tracked position; andoutputs the likely future tracked position as the expected location of the head to said device, wherein said likely future tracked position is utilized to control an auto-stereoscopic display for the display of images for eyes located at expected positions corresponding to said likely future tracked position, and wherein Kalman filtering of the previously tracked positions is utilized in estimating said likely future tracked position.
  • 8. The face tracking system of claim 7 wherein utilizing previously tracked positions to estimate a likely future tracked position further includes the steps of: (a) providing the expected location of the head using facial features detected in a previous and current tracked positions;(b) deriving a first covariance matrix from the current prediction and a previous covariance matrix;(c) utilising the current prediction from step (a) and a Kalman filter to determine a corresponding projected point of the facial feature on the plane of at least one camera;(d) deriving a Jacobain of the projected points from step (c);(e) deriving a residual covariance of the projected points from step (c);(f) deriving a suitable filter gain for the Kalman filter;(g) deriving a suitable update coefficients for the first covariance matrix; and(h) updating the Kalman filter utilising the filter gain.
  • 9. A system for providing an expected location of a head the system comprising: video input means for providing at least one video signal of the head;first processing means for processing the video signal so as to output a substantially continuous series of current head location data;second processing means for processing predetermined one of the current head location data so as to output a predicted future expected head position as the expected location of the head, thereby substantially overcoming a delay in processing the head position, wherein said second processing means utilizes a Kalman filtering of the current head location data; andan auto-stereoscopic display driven by said predicted expected location output of said head.
  • 10. The system of claim 9 wherein said video input means includes stereo video inputs.
  • 11. The system of claim 9 further comprising an auto-stereoscopic display driven by said predicted expected location output of said head.
Priority Claims (1)
Number Date Country Kind
2003901528 Mar 2003 AU national
Parent Case Info

This application is a continuation of pending International Patent Application No. PCT/AU2004/000413 filed on Mar. 31, 2004 which designates the United States and claims priority of Australian Patent Application No. 2003901528 filed on Mar. 31, 2003.

US Referenced Citations (7)
Number Name Date Kind
5771066 Barnea Jun 1998 A
5912721 Yamaguchi et al. Jun 1999 A
5978143 Spruck Nov 1999 A
6404900 Qian et al. Jun 2002 B1
6459446 Harman Oct 2002 B1
6674877 Jojic et al. Jan 2004 B1
7127081 Erdem Oct 2006 B1
Foreign Referenced Citations (5)
Number Date Country
1 154 655 Nov 2001 EP
2 324 428 Oct 1998 GB
9720244 Jun 1997 WO
0241128 May 2002 WO
02054132 Jul 2002 WO
Related Publications (1)
Number Date Country
20060146046 A1 Jul 2006 US
Continuations (1)
Number Date Country
Parent PCT/AU2004/000413 Mar 2004 US
Child 11241669 US