The present invention relates to head animation, and more particularly, to generating an animated three-dimensional video head.
Virtual spaces filled with avatars are an attractive way to allow for the experience of a shared environment. However, animation of a photo-realistic avatar generally requires intensive graphic processes, particularly for rendering facial features.
Accordingly, there exists a significant need for improved rendering of facial features. The present invention satisfies this need.
The present invention provides a technique for for translating facial animation values to head mesh positions for rendering facial features of an animated avatar. In the method, an animation vector a of dimension Na is provided. Na is the number of facial animation values in the animation vector. A mapping algorithm F is applied to the animation vector to generate a target mix vector g of dimension M. M is the number of targets associated with the head mesh positions. The head mesh positions are deformed based on the target mix vector g.
In more detailed features of the invention, the animation vector may be determined by sensing facial features, associated with independent targets, of a human subject. Alternatively, the animation vector is obtained by applying facial feature position sensing to annotated vertices on a model. The targets may comprise independent targets and dependent targets. The dependent targets may be generated by linear superposition of vertex coordinates of the independent targets.
A calibration vector c may be applied to the target mix vector g to generate a calibrated target mix vector gc. The head mesh positions are then deformed based on the calibrated target mix vector.
Groups that associate sets of the animation values with sets of targets may be defined. The mapping algorithm F may be independently applied to each grouped set of animation values to generate corresponding target mix group-vectors. The target mix group-vectors may be combined to generate the target mix vector g of dimension M, where M is the number of targets associated with the head mesh positions; and
Other features and advantages of the present invention should be apparent from the following description of the preferred embodiments taken in conjunction with the accompanying drawings, which illustrate, by way of example, the principles of the invention.
The present invention provides a technique for translating an animation vector to a target mix vector. With reference to
The animation vector is an abstracted sensing result and is a most compact representation of the facial expression as determined by audio-visual sensing. By definition, the animation vector is zero for the neutral facial expression, shown by the first morph target. The target mix vector describes the contribution of each individual target to the current expression.
Different mapping algorithms may be used for mapping the animation vectors to the target mix vectors. The common goal of the mapping algorithms is to provide a reasonable interpolation between the points in animation space associated with the targets. Each mapping algorithm is exactly defined by a set of parameters. The parameters vary from algorithm to algorithm.
Calibration may be performed by multiplying the target mix vector with a diagonal matrix. Since the matrix is diagonal, it is henceforth referred to as the calibration vector. The calibration vector may be eliminated if the head model driven by the target mix vector is derived from the same human actor used to obtain the animation vector.
An overview of the translation from animation vector to target mix vector follows. Be a the animation vector of dimension Na (number of animation values), g the target mix vector of dimension M (number of independent targets) and p1, . . . , pL the parameters of mapping algorithm F( ), then
g=F(a, p1, . . . , pL) (1)
The calibrated target mix vector is obtained by multiplying with the diagonal matrix defined by the calibration vector c:
Further, be t1, . . . , tM the parameterization of the targets associated with the components of the target mix vector and t0 the parameterization of the neutral model, then the parameterization of the current expression can be obtained by a simple matrix multiplication:
t=(t1−t0 t2−t0 . . . tM−t0)·gc+t0 (3)
The matrix (t1−t0 t2−t0 . . . tM−t0) is referred to as the target matrix T.
A description of the mapping algorithm follows. Every target ti is associated with an animation vector a1. The target and animation vectors are connected by this deliberate association only. This will become apparent in the formulas, where targets enter only as the i-th unity vector ei representing the target mix vector that results in exactly that target. The parameterization of the target t1 is not relevant for computation of the mapping parameters. (This means that is does not matter if the target is defined by vertex positions, morph link positions or muscle tensions.)
The animation vector can be set manually or it can be derived from a reference model with targets, if the model is equipped with ground truth anchors that enable the application of the sensing algorithm to the model and it's deformations AND if the reference model implements all needed targets. The reference model must have a human geometry, since the purpose of the ground truth anchors is to simulate tracking on the model. Manual editing is necessary if the animation vector contains elements that cannot be derived from visual sensing, such as lip synch animation values.
The mapping is basically a multidimensional interpolation in a between the target points. The mapping parameters are determined by minimizing the error in reproducing the target points. Depending on the mapping algorithm, perfect reproduction of the target points may not be possible.
The parameters p1, . . . , pL of the mapping are determined by solving the set of equations
where ei is a unity vector representing exact reproduction of the target points.
Targets can be divided into independent groups of targets, such as eye region targets and mouth region targets. Different mapping algorithms can be applied to the different groups to achieve more flexibility.
A description of types of mapping algorithms follow. The simplest mapping algorithm is the linear mapping:
F(a,P)=P·a (5)
The parameter matrix is determined by solving the equation
P·(a1 a2 . . . aM)=11 (6)
using singular value decomposition (SVD). If Na<M, the equation (6) is overdetermined and SVD will return the “solution” that satisfies eq. (4). If Na>M, the equation is underdetermined and the SVD solution will be the vector with the smallest norm that satisfies equation (6). Representative SVD techniques are described in great detail in “Numerical Recipes in C, the Art of Scientific Computing”, William H. Press et al., Cambridge University Press, second edition 1993.
A more general mapping is achieved by using a set of basis functions as input. Obviously, the linear method is a special case of this more general method.
F(a,P)=P·B(a) (7)
The solution is analog to the solution of the linear problem. Since the number of basis functions is independent of the number of animation values Na, it is always possible to choose exactly M functions, so that the system is neither over- or underdetermined:
P=(B(a1) B(a2) . . . B(aM))−1 (8)
The basis functions can be chosen manually by carefully analyzing the animation vectors of the participating targets. It is very tedious to design a basis function manually such that it only responds when the associated target is acted and not responds when any other target is acted. Off-diagonal elements of the P matrix lead to corrections and decouple the targets such that this desired behavior is achieved.
The target matrix, calibration matrix and mapping parameter matrix can be combined into one matrix by simple multiplication, which can be done ahead of time:
The decoder matrix D, offset vector (neutral target) t0 and definition of basis functions B(a) are the most compact representation of the algorithm parameters.
A description of radial/angular basis function mapping follows. The basis functions can also be a set of automatically determined functions such as radial basis functions. Certain properties of the mapping algorithm can be designed into the basis functions. A desirable property of the algorithm is scale linearity:
F(λa)=λF(a) (10)
This means that the “amplitude” of an expression is translated linearly to the model. It can be obtained by using basis functions of the form
Because of their linear or locally linear behavior, the following functions are useful for {tilde over (b)} (∂):
{tilde over (b)}(∂)=∂ (preferred) (12)
Equations (13) and (14) can be locally approximated by Equation (12) and have interesting saturation characteristics. The parameter α determines the localization of these basis functions.
All mapping algorithms are somewhat arbitrary because of their nature as interpolation algorithms. The way interpolation is done between targets is determined by the choice of basis functions. It seems to be reasonable to shoot for the algorithm that delivers the most linear interpolation possible while still reproducing all targets.
If the basis functions are determined automatically, it is easily possible to add K dependent targets that are created by linear superposition of independent targets. This enables one to have more control over the interpolation process by providing additional support points. Eq. (4) is then generalized to:
Each dependent target is defined by its animation vector ai and target mix vector gi with i>M, which defines the superposition of independent targets. Eq. (8) has to be modified to yield the solution to this more general problem:
P=(g1 g2 . . . gM+K)·(B(a1) B(a2) . . . B(aM+K))31 1 (16)
A human anatomy reference model with targets ti is created. The vertices that correspond to tracking (motion capture) nodes are annotated. The displacement of these annotated vertices relative to the neutral model are used to compute the animation vector ai for each target. Alternatively, the animation vectors can be obtained by applying tracking to a human actor who performs expressions that are equivalents to all of the targets. This will remove the need to annotate certain vertices.
Advantageously, 8 to 22 (or more) facial tracking nodes may be used to define and animate mouth, eyes, eyebrows, nose and head angle. The components of the animation vector may be displacement values relative to neutral face values, either of the tracking node coordinates directly or of combinations of tracking nodes, such as the distance between nodes or linear combinations of node coordinates. Representative facial components of the animation vector for the mouth may include vertical mouth position, horizontal mouth position, mouth width, lip distance, and mouth corner positions (left and right).
Optionally, groups are defined which consist of a selection of certain components of the animation vector and certain targets. (Example: All components of the animation vector describing the mouth and all targets containing mouth expressions) The algorithm can then be applied independently to each group.
Having obtained the associations (ai, ti) for all targets, the algorithm parameters are determined by the process described above. Our preferred implementation uses the mapping according to Eq. (7) with basis functions according to Eqs. (11) and (12).
Morphing of a texture map on a deformed three-dimensional head mesh is described in U.S. Pat. No. 6,272,231, titled WAVELET-BASED FACIAL MOTION CAPTURE FOR AVATAR ANIMATION. Imaging systems for acquiring images and image mapping are described in U.S. patent application Ser. No. 09/724,320, titled METHOD AND APPARATUS FOR RELIEF TEXTURE MAP FLIPPING. The entire disclosures of U.S. Pat. No. 6,272,231 and U.S. patent application Ser. No. 09/724,320 are incorporated herein by reference.
Although the foregoing discloses the preferred embodiments of the present invention, it is understood that those skilled in the art may make various changes to the preferred embodiments without departing from the scope of the invention. The invention is defined only by the following claims.
Priority is claimed as a continuation-in-part of U.S. application Ser. No. 09/929,823 filed Aug. 13, 2001, which is incorporated herein by this reference.
Number | Name | Date | Kind |
---|---|---|---|
4725824 | Yoshioka | Feb 1988 | A |
4805224 | Koezuka et al. | Feb 1989 | A |
4827413 | Baldwin et al. | May 1989 | A |
5159647 | Burt | Oct 1992 | A |
5168529 | Peregrim et al. | Dec 1992 | A |
5187574 | Kosemura et al. | Feb 1993 | A |
5220441 | Gerstenberger | Jun 1993 | A |
5280530 | Trew et al. | Jan 1994 | A |
5333165 | Sun | Jul 1994 | A |
5383013 | Cox | Jan 1995 | A |
5430809 | Tomitaka | Jul 1995 | A |
5432712 | Chan | Jul 1995 | A |
5511153 | Azarbayejani et al. | Apr 1996 | A |
5533177 | Wirtz et al. | Jul 1996 | A |
5550928 | Lu et al. | Aug 1996 | A |
5581625 | Connell | Dec 1996 | A |
5588033 | Yeung | Dec 1996 | A |
5680487 | Markandey | Oct 1997 | A |
5699449 | Javidi | Dec 1997 | A |
5714997 | Anderson | Feb 1998 | A |
5715325 | Bang et al. | Feb 1998 | A |
5719954 | Onda | Feb 1998 | A |
5736982 | Suzuki et al. | Apr 1998 | A |
5764803 | Jacquin et al. | Jun 1998 | A |
5774591 | Black et al. | Jun 1998 | A |
5802220 | Black et al. | Sep 1998 | A |
5809171 | Neff et al. | Sep 1998 | A |
5828769 | Burns | Oct 1998 | A |
5917937 | Szeliski et al. | Jun 1999 | A |
5982853 | Liebermann | Nov 1999 | A |
5995119 | Cosatto et al. | Nov 1999 | A |
6011562 | Gagne et al. | Jan 2000 | A |
6044168 | Tuceryan et al. | Mar 2000 | A |
6052123 | Lection et al. | Apr 2000 | A |
6307576 | Rosenfield | Oct 2001 | B1 |
6351269 | Georgiev | Feb 2002 | B1 |
6504546 | Cosatto et al. | Jan 2003 | B1 |
6600491 | Szeliski et al. | Jul 2003 | B1 |
6611268 | Szeliski et al. | Aug 2003 | B1 |
6654018 | Cosatto et al. | Nov 2003 | B1 |
Number | Date | Country |
---|---|---|
4406020 | Jun 1995 | DE |
0807902 | Nov 1997 | EP |
WO9953443 | Oct 1999 | WO |
Number | Date | Country | |
---|---|---|---|
20030043153 A1 | Mar 2003 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 09929823 | Aug 2001 | US |
Child | 10216662 | US |