Multi-pose face tracking using multiple appearance models

<u style="single">Background</u>

0001Tracking faces in a video sequence is an important module for automated video surveillance. It is a precursor to many applications such as video based face recognition, tag and track of faces and multi-camera indexing. Face tracking in a video has been a long studied problem using many features, such as skin color and edge based face structure features. Tracking faces poses unique set of problems when compared with tracking other objects, such as people, cars, or other objects of interest.

0002Faces are approximately uniform in their color which makes the tracking solution possible using color as an appearance model. Many researchers have used features derived from the skin color of the face, such as color histogram, for face tracking. Using the face color as an appearance model for tracking provides invariance to different head pose variations. However, the problem of face tracking using color is challenging when the background is of similar color or in the presence of ambient illumination variations. Using edge information of faces as an appearance model for face tracking proved robust to illumination variations. However, the out-of-plane variations of the face pose worsens 2D edge model matching. A generalized tracking algorithm has also been used that models the appearance using a mixture of Gaussians. It may be used for tracking a face with pose changes, typically in-plane pose changes. Others may use an appearance model and embed the tracking in the particle filter framework.

0003Two main approaches have been used to track faces in videos. In one approach, local features of the face are detected (or manually marked) and the features are tracked over time. This is useful if the orientation of the face needs to be computed along with the face position (as used in Human Computer Interaction applications). Another approach utilizes global features of the face, such as a color histogram, that distinguishes the face from the background.

0004In surveillance videos, multiple faces need to be tracked with face sizes as small as 24x24 pixels, making it difficult to identify and track local features. The faces can undergo illumination changes (because of shadows and indoor lighting), can have partial occlusions, and have large pose changes. The background may also be cluttered, depending on the setup. These challenges need to be overcome for effective face tracking in a surveillance setup.

<u style="single">Brief Description of the Drawings</u>

0005<figref idref="f0001">FIGs. 1A and 1B</figref> are graphs illustrating difference of Gaussian (DOG) filter models according to an example embodiment.

0006<figref idref="f0002">FIGs. 2A, 2B and 2C</figref> illustrate the output of the filters of <figref idref="f0001">FIGs. 1A and 1B</figref> for a face according to an example embodiment.

0007<figref idref="f0002">FIG. 3</figref> is a pseudocode representation of predicting positions of a face in tracks according to an example embodiment.

0008<figref idref="f0003">FIG. 4</figref> is a flowchart with pseudocode illustrating a method of updating a particle filter according to an example embodiment.

0009<figref idref="f0004">FIG. 5</figref> is a flowchart representation of a method of tracking a face using multiple models according to an example embodiment.

0010<figref idref="f0004">FIG. 6</figref> is a block representation of a predicted position according to an example embodiment.

0011<figref idref="f0005">FIG. 7</figref> is a series of frames from a surveillance video illustrating face tracking between the frames according to an example embodiment.

0012<figref idref="f0006">FIG. 8</figref> is a series of frames from a surveillance video illustrating face tracking between the frames with illumination change according to an example embodiment.

0013<figref idref="f0007">FIG. 9</figref> is a block diagram of a typical computer system used to implement methods according to an example embodiment.

<u style="single">Detailed Description</u>

0014In the following description, reference is made to the accompanying drawings that form a part hereof, and in which is shown by way of illustration specific embodiments which may be practiced. These embodiments are described in sufficient detail to enable those skilled in the art to practice the invention, and it is to be understood that other embodiments may be utilized and that structural, logical and electrical changes may be made without departing from the scope of the present invention. The following description of example embodiments is, therefore, not to be taken in a limited sense, and the scope of the present invention is defined by the appended claims.

0015The functions or algorithms described herein may be implemented in software or a combination of software, hardware and human implemented procedures in one embodiment. The software may consist of computer executable instructions stored on computer readable media such as memory or other type of storage devices. The term "computer readable media" is also used to represent any means by which the computer readable instructions may be received by the computer, such as by different forms of wired or wireless transmissions. Further, such functions correspond to modules, which are software, hardware, firmware or any combination thereof. Multiple functions may be performed in one or more modules as desired, and the embodiments described are merely examples. The software may be executed on a digital signal processor, ASIC, microprocessor, or other type of processor operating on a computer system, such as a personal computer, server or other computer system.

0016A system and method are provided for tracking a face moving through multiple frames of a video. A predicted position of a face in a video frame is obtained. Localized search for both a color model and an edge model are performed to derive similarity values, represented as correlation values, for each position in the search area defined around the predicted position. The correlation values are then combined to determine a best position match to track a face in the video.

0017In the description, multiple independent appearance models used for face representation are first described, followed by a method of predicting the position of a face. Finally, methods of combining the results of the models to track the face are described.

0018In one embodiment, representations from an edge representation and a color representation for a face appearance model are combined to track faces. The color representation in one embodiment is a YCbCr color space of an input image. The edge representation in one embodiment is the output of a Difference of Gaussian (DOG) filter on a grayscale input image. The two models are independent in one embodiment, and may be used to represent the appearance of a face, and track it over multiple frames of a video. In one embodiment, each model provides correlation values that may be combined by multiplication to provide a final correlation matrix which may be used for tracking.

0019DOG filters as illustrated in <figref idref="f0001">FIGs. 1A and 1B. FIG. 1A</figref> illustrates a response of an on-center, off surround filter, and <figref idref="f0001">FIG. 1B</figref> illustrates the response of an off-center, on surround filter. A DOG filter has a response that is high around the edges of the image, i.e., where there are variations in intensity, and low (or zero) over uniform regions. In one embodiment, the DOG filter may be referred to as an edge based model or face sketch. Even though responsive to edges, DOG filters are un-oriented filters and hence are not sensitive to orientation. Cells with such properties have been found in the Lateral Geniculate Nucleus (LGN) layer of the human visual system, which is the first stage of visual processing. Such filters have contrast normalization properties which help in dealing with illumination changes. DOG filters have also been shown to be useful for face recognition. DOG filters provide good recognition rate and may be resistant to slight pose changes.

0020For the edge based model, the output of <i>Difference of Gaussians</i> (DOG) on the grayscale face image is used. There are two types of filters, on-center off-surround and off center, on-surround. Though, the output of one filter is sufficient to obtain the output of another (basically the output of off-center is the negative of the on-center output), two filters are used so that the representation of the appearance model is positive. In one embodiment, the positive kernel is defined and the negative output may be obtained by negative of the ON output before thresholding. The DOG is defined as: <maths id="math0001"><math display="block"><mi mathvariant="italic">ON</mi><mfenced><mi>x</mi><mo>⁢</mo><mi>y</mi></mfenced><mo>=</mo><mfrac><mn>1</mn><msqrt><mn>2</mn><mo>⁢</mo><mi>π</mi><mo>⁢</mo><msub><mi>σ</mi><mi>c</mi></msub></msqrt></mfrac><mo>⁢</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mfenced><msup><mi>x</mi><mn>2</mn></msup><mo>+</mo><msup><mi>y</mi><mn>2</mn></msup></mfenced><mrow><mn>2</mn><mo>⁢</mo><msubsup><mi>σ</mi><mi>c</mi><mn>2</mn></msubsup></mrow></mfrac></mrow></msup><mo>-</mo><mfrac><mn>1</mn><msqrt><mn>2</mn><mo>⁢</mo><mi>π</mi><mo>⁢</mo><msub><mi>σ</mi><mi>s</mi></msub></msqrt></mfrac><mo>⁢</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mfenced><msup><mi>x</mi><mn>2</mn></msup><mo>+</mo><msup><mi>y</mi><mn>2</mn></msup></mfenced><mrow><mn>2</mn><mo>⁢</mo><msubsup><mi>σ</mi><mi>s</mi><mn>2</mn></msubsup></mrow></mfrac></mrow></msup></math><img file="EP2017769A2_D0001.tif" /></maths> where, σ<i>_C</i> (3) is the standard deviation for the center Gaussian and σ<i>_S</i> = K σ<i>_C</i> is for the surround Gaussian, where K is a constant. The DOGs are normalized such that the sum of the positive values equals 1 and the sum of negative values equals 1. This results in zero response for uniform regions. Using the center surround filters we obtain two set of features: <maths id="math0002"><math display="block"><mtable columnalign="left"><mtr><mtd><mi mathvariant="italic">ONOUT</mi><mo>=</mo><mi mathvariant="italic">ON</mi><mo>*</mo><mi>I</mi></mtd></mtr><mtr><mtd><mi mathvariant="italic">OFFOUT</mi><mo>=</mo><mo>-</mo><mi mathvariant="italic">ONOUT</mi></mtd></mtr><mtr><mtd><mi mathvariant="italic">FONOUT</mi><mo>=</mo><msup><mfenced open="[" close="]"><mi mathvariant="italic">ONOUT</mi></mfenced><mo>+</mo></msup></mtd></mtr><mtr><mtd><mi mathvariant="italic">FOFFOUT</mi><mo>=</mo><msup><mfenced open="[" close="]"><mi mathvariant="italic">OFFOUT</mi></mfenced><mo>+</mo></msup></mtd></mtr><mtr><mtd><mi>where</mi><mspace width="2em" /><msup><mfenced open="[" close="]"><mi>x</mi></mfenced><mo>+</mo></msup><mo>=</mo><mi>x</mi><mo>,</mo><mi>x</mi><mo>></mo><mn>0</mn><mo>,</mo><mn>0</mn><mspace width="1em" /><mi>otherwise</mi></mtd></mtr></mtable></math><img file="EP2017769A2_D0002.tif" /></maths>

0021<figref idref="f0002">FIGs. 2A, 2B and 2C</figref> illustrate the output of the DOG filter. <figref idref="f0002">FIG. 2A</figref> is the input image. <figref idref="f0002">FIG. 2B</figref> is the ON filter output and <figref idref="f0002">FIG. 2C</figref> is the OFF filter output.

0022Color template appearance models are basically snapshots of the objects in an appropriate color space. In one embodiment, the YCbCr color space is used for representing the face. YCbCr color space may be more discriminative than other color spaces between skin and non-skin pixels. The YCbCr space may be created by converting RGB space to YCbCr space using various conversion functions. Representing the face in such a space may provide maximal separation between face and non-face pixels and hence help to track the face better. In one embodiment, illumination and pose changes are handled by combining results from the DOG and color template appearance models. Using both models and combining them facilitates tracking faces reliably in surveillance videos under head pose and ambient light illumination variations.

0023An overall computer implemented method 300 for tracking faces is illustrated in a pseudocode in <figref idref="f0002">FIG. 3</figref>. The face tracking method 300 uses detections from a face detector and tracks the face over time. An input for method 300 is a set of detections from each model for a particular frame of video.

0024For each track of video, as indicated at 305, a predicted position of the face is obtained at 310. If the predicted position is out of frame, the track may be deleted at 315. At 320, the best match around the predicted position is found.

0025At 325, detections with the tracks are associated based on an overlap between the position of the best match and the position of the detection. If more than one track is associated with a detection at 330, associations other than the one that has the best overlap are removed.

0026At 335, for each detection that is not associated with a track, a new track is started on the detection. If a track is not associated with detections and the track has lived for less than two frames, it may be deleted at 340.

0027A particle filter in one embodiment is used to provide the prediction of the position of the track over time. Other filters may be used in further embodiments. Four states of a face may be tracked in one embodiment. They are top-left position (x, y) and velocity (vx, vy). A weight for each particle is also maintained. For each track, the following variables are maintained: <ol id="ol0001"><li>1. <maths id="math0003"><math display="inline"><msubsup><mi>p</mi><mi>i</mi><mi>x</mi></msubsup><mfenced><mi>t</mi></mfenced></math><img file="EP2017769A2_D0003.tif" /></maths> top-left <i>x</i> position of the particle <i>i</i> at time <i>t</i></li><li>2. <maths id="math0004"><math display="inline"><msubsup><mi>p</mi><mi>i</mi><mi>y</mi></msubsup><mfenced><mi>t</mi></mfenced></math><img file="EP2017769A2_D0004.tif" /></maths> top-left <i>y</i> position of the particle i at time t</li><li>3. <maths id="math0005"><math display="inline"><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">xv</mi></msubsup><mfenced><mi>t</mi></mfenced></math><img file="EP2017769A2_D0005.tif" /></maths><i>velocity</i> in x direction of the particle i at time t</li><li>4. <maths id="math0006"><math display="inline"><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">yv</mi></msubsup><mfenced><mi>t</mi></mfenced></math><img file="EP2017769A2_D0006.tif" /></maths><i>velocity</i> in $y$ direction of the particle i at time t</li><li>5. <i>wt_i</i>(<i>t</i>) - weight for particle i at time t</li></ol>

0028When a face is detected by the detector, a new track is created and the corresponding particles are initialized. Each track is initialized with a set of 1000 particles with an initial constant velocity of 1 and the position of the particles distributed around the top-left position of the detected face. In one embodiment, 1000 particles are provided with an initial position around the top-left corner of the detected face, and the initial weight set to a constant <maths id="math0007"><math display="block"><msubsup><mi>p</mi><mi>i</mi><mi>x</mi></msubsup><mfenced><mn>1</mn></mfenced><mo>=</mo><msup><mi mathvariant="italic">obs</mi><mi>x</mi></msup><mo>+</mo><mi>σ</mi><mo>*</mo><mi>G</mi><mfenced><mi>i</mi></mfenced></math><img file="EP2017769A2_D0007.tif" /></maths><maths id="math0008"><math display="block"><msubsup><mi>p</mi><mi>i</mi><mi>y</mi></msubsup><mfenced><mn>1</mn></mfenced><mo>=</mo><msup><mi mathvariant="italic">obs</mi><mi>y</mi></msup><mo>+</mo><mi>σ</mi><mo>*</mo><mi>G</mi><mfenced><mi>i</mi></mfenced></math><img file="EP2017769A2_D0008.tif" /></maths><maths id="math0009"><math display="block"><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">xv</mi></msubsup><mfenced><mn>1</mn></mfenced><mo>=</mo><msub><mi>X</mi><mi mathvariant="italic">VEL</mi></msub><mo>+</mo><mi>σ</mi><mo>*</mo><mi>G</mi><mfenced><mi>i</mi></mfenced></math><img file="EP2017769A2_D0009.tif" /></maths><maths id="math0010"><math display="block"><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">xy</mi></msubsup><mfenced><mn>1</mn></mfenced><mo>=</mo><msub><mi>Y</mi><mi mathvariant="italic">VEL</mi></msub><mo>+</mo><mi>σ</mi><mo>*</mo><mi>G</mi><mfenced><mi>i</mi></mfenced></math><img file="EP2017769A2_D0010.tif" /></maths><maths id="math0011"><math display="block"><mi>w</mi><mo>⁢</mo><msub><mi>t</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac></math><img file="EP2017769A2_D0011.tif" /></maths> where <i>obs^x</i>(<i>t</i>), <i>obs^y</i>(<i>t</i>) is the observed (x, y) position (from the face detector) of the track at time $t$, <i>X_VEL</i>, <i>Y_VEL</i> is the initial velocity, <i>G</i>(<i>i</i>) is a Gaussian random number, σ is the variance, and N is the number of particles. The particles would be distributed around the detected position depending on σ. σ is set to 2 in one embodiment of the tracker.

0029For any time t the prediction of the position from the particle filter is obtained as follows: <maths id="math0012"><math display="block"><mover><mi>x</mi><mo>^</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mstyle displaystyle="false"><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msubsup><mi>p</mi><mi>i</mi><mi>x</mi></msubsup></mstyle></math><img file="EP2017769A2_D0012.tif" /></maths><maths id="math0013"><math display="block"><mover><mi>y</mi><mo>^</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mstyle displaystyle="false"><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msubsup><mi>p</mi><mi>i</mi><mi>y</mi></msubsup></mstyle></math><img file="EP2017769A2_D0013.tif" /></maths>

0030The particle filter is also updated at various times as illustrated at 400 in <figref idref="f0003">FIG. 4</figref>. Let the observed position (either from detection or the obtained by search) be <i>obs^x</i>(<i>t</i>),<i>obx</i>^y(<i>t</i>). The velocity estimate may be obtained at 405 by <maths id="math0014"><math display="block"><msubsup><mi mathvariant="italic">obs</mi><mi>v</mi><mi>x</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>=</mo><msup><mi mathvariant="italic">obs</mi><mi>x</mi></msup><mfenced><mi>t</mi></mfenced><mo>-</mo><msup><mi mathvariant="italic">obs</mi><mi>x</mi></msup><mo>⁢</mo><mfenced><mi>t</mi><mo>-</mo><mn>1</mn></mfenced><mspace width="1em" /><mi>and</mi><mspace width="1em" /><msubsup><mi mathvariant="italic">obs</mi><mi>v</mi><mi>y</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>=</mo><msup><mi mathvariant="italic">obs</mi><mi>y</mi></msup><mfenced><mi>t</mi></mfenced><mo>-</mo><msup><mi mathvariant="italic">obs</mi><mi>y</mi></msup><mo>⁢</mo><mfenced><mi>t</mi><mo>-</mo><mn>1</mn></mfenced></math><img file="EP2017769A2_D0014.tif" /></maths>

0031The following steps are used to update the particle filter: Recompute the weights of the particles at 410 based on the observation: <maths id="math0015"><math display="block"><mi>w</mi><mo>⁢</mo><msub><mi>t</mi><mi>i</mi></msub><mo>=</mo><mi mathvariant="italic">exp</mi><mfenced><mo>-</mo><mn>0.5</mn><mo>*</mo><mfenced><mfrac><mrow><msup><mfenced><msubsup><mi>p</mi><mi>i</mi><mi>x</mi></msubsup><mo>-</mo><msup><mi mathvariant="italic">obs</mi><mi>x</mi></msup></mfenced><mn>2</mn></msup><mo>+</mo><msup><mfenced><msubsup><mi>p</mi><mi>i</mi><mi>y</mi></msubsup><mo>-</mo><msup><mi mathvariant="italic">obs</mi><mi>y</mi></msup></mfenced><mn>2</mn></msup></mrow><msub><mi>σ</mi><mi>d</mi></msub></mfrac><mo>+</mo><mfrac><mrow><msup><mfenced><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">xv</mi></msubsup><mo>-</mo><msup><mi mathvariant="italic">obs</mi><mi mathvariant="italic">xc</mi></msup></mfenced><mn>2</mn></msup><mo>+</mo><msup><mfenced><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">xv</mi></msubsup><mo>-</mo><msup><mi mathvariant="italic">obs</mi><mi mathvariant="italic">xc</mi></msup></mfenced><mn>2</mn></msup></mrow><msub><mi>σ</mi><mi mathvariant="italic">dv</mi></msub></mfrac></mfenced></mfenced></math><img file="EP2017769A2_D0015.tif" /></maths> where σ<i>_d</i> and σ<i>_dv</i> are for position and velocity respectively. This determines the distribution of the weight values across the particles.

0032The particles are resampled at 415 based on the weights. Around (50+P)% of the new particles are generated by resampling the original set of particles based on their weights, where P ranges from 0 to 50. A sampling with replacement algorithm may be used to resample the particles.

0033The remaining (50-P)% of the particles may be populated using the observation at 445. These particles would have their states distributed around the observation, just like in the initialization. The state of the particles is then updated at 450. In one embodiment, the state of the particles are updated using: <maths id="math0016"><math display="block"><msubsup><mi>p</mi><mi>i</mi><mi>x</mi></msubsup><mo>⁢</mo><mfenced><mi>t</mi><mo>+</mo><mn>1</mn></mfenced><mo>=</mo><msubsup><mi>p</mi><mi>i</mi><mi>x</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>+</mo><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">xv</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>+</mo><mi>σ</mi><mo>*</mo><mi>G</mi><mfenced><mi>i</mi></mfenced></math><img file="EP2017769A2_D0016.tif" /></maths><maths id="math0017"><math display="block"><msubsup><mi>p</mi><mi>i</mi><mi>y</mi></msubsup><mo>⁢</mo><mfenced><mi>t</mi><mo>+</mo><mn>1</mn></mfenced><mo>=</mo><msubsup><mi>p</mi><mi>i</mi><mi>y</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>+</mo><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">yv</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>+</mo><mi>σ</mi><mo>*</mo><mi>G</mi><mfenced><mi>i</mi></mfenced></math><img file="EP2017769A2_D0017.tif" /></maths><maths id="math0018"><math display="block"><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">xv</mi></msubsup><mo>⁢</mo><mfenced><mi>t</mi><mo>+</mo><mn>1</mn></mfenced><mo>=</mo><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">xv</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>+</mo><mi>σ</mi><mo>*</mo><mi>G</mi><mfenced><mi>i</mi></mfenced></math><img file="EP2017769A2_D0018.tif" /></maths><maths id="math0019"><math display="block"><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">yx</mi></msubsup><mo>⁢</mo><mfenced><mi>t</mi><mo>+</mo><mn>1</mn></mfenced><mo>=</mo><msubsup><mi>p</mi><mi>i</mi><mi mathvariant="italic">yv</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>+</mo><mi>σ</mi><mo>*</mo><mi>G</mi><mfenced><mi>i</mi></mfenced></math><img file="EP2017769A2_D0019.tif" /></maths><maths id="math0020"><math display="block"><mi>w</mi><mo>⁢</mo><msub><mi>t</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac></math><img file="EP2017769A2_D0020.tif" /></maths> where σ is a constant.

0034A search may be performed to find the best match as illustrated in <figref idref="f0004">FIG. 5</figref> at 500. In order to find the location of the track in the current frame, a search using the appearance model around the predicted position is done in one embodiment to find the best match which is then used to update the appearance model.

0035At 510, the predicted position is obtained from the particle filter (x, y). At 515, the area of search is fixed to be (x₁, y₁, x₂y₂) where, <maths id="math0021"><math display="block"><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mi>x</mi><mo>-</mo><mi mathvariant="italic">appModelWidth</mi></math><img file="EP2017769A2_D0021.tif" /></maths><maths id="math0022"><math display="block"><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mi>y</mi><mo>-</mo><mi mathvariant="italic">appModelWidth</mi></math><img file="EP2017769A2_D0022.tif" /></maths><maths id="math0023"><math display="block"><msub><mi>x</mi><mn>2</mn></msub><mo>=</mo><mi mathvariant="italic">x</mi><mo mathvariant="italic">+</mo><mi mathvariant="italic">appModelHeight</mi></math><img file="EP2017769A2_D0023.tif" /></maths><maths id="math0024"><math display="block"><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><mi mathvariant="italic">y</mi><mo mathvariant="italic">+</mo><mi mathvariant="italic">appModelHeight</mi></math><img file="EP2017769A2_D0024.tif" /></maths>

0036That is, the area of search is twice the size of the appearance model. The size of the search area is determined by the face motion dynamics and also by tracking accuracy. If the search area is large, then the face is likely to be within the search area and hence the search mechanism is more likely to locate the face correctly. On the other hand, if the search area is large, then the time taken to perform the search becomes more, leading to reduction in computational speed. Hence, the search area is determined based on the scene and face dynamics. For scale = 1-p,1,1+q, where p and q are constants that are optimized based on face motion, the appearance model may be resized to match the scale at 520. At 525, <i>Ĉ</i> = correlation values are obtained by template matching within the search area for color template appearance model. <i>D̂</i> = correlation values are obtained at 530 by template matching within the search area for DOG appearance model. The correlation values are then multiplied at 540, <i>F̂</i> = <i>Ĉ</i> * <i>D̂</i> where <i>F̂</i> is the final correlation value matrix. Also at 540, the best position for scale s is <maths id="math0025"><math display="inline"><mfenced><msubsup><mi>x</mi><mi>b</mi><mi>s</mi></msubsup><mo>⁢</mo><msubsup><mi>y</mi><mi>b</mi><mi>s</mi></msubsup></mfenced><mo>=</mo><msub><mi>argmax</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></msub><mfenced><mi>F</mi></mfenced></math><img file="EP2017769A2_D0025.tif" /></maths>with correlation value <i>V_s</i>. Out of three positions the best one is chosen based on the maximum correlation value. <maths id="math0026"><math display="block"><mfenced><msub><mi>x</mi><mi>b</mi></msub><mo>⁢</mo><msub><mi>y</mi><mi>b</mi></msub></mfenced><mo>=</mo><msub><mi>argmax</mi><mi>s</mi></msub><mo>⁢</mo><msub><mi>V</mi><mi>s</mi></msub></math><img file="EP2017769A2_D0026.tif" /></maths>

0037The best match is searched for at 545. <figref idref="f0004">FIG. 6</figref> is a graphical representation of the search for the best match at 600. x , y at 610 indicates the predicted position. Solid line 615 indicates the frame, dashed line 620 indicates the search area and dotted line 625 indicates the size of the appearance model.

0038The correlation values obtained by search using the two appearance models are multiplied to obtain the final correlation matrix as indicated above at 540. The final correlation matrix is then used to find the best match position at 545. Multiplication ensures that the final best position is the one that has high correlation value for both the appearance models. One interesting aspect of this approach is that the range of the correlation values is different for each appearance model. For example, the range for color template appearance model in one embodiment is from 0.96 to 1 while the range for DOG appearance model is from 0.3 to 1. Because of the difference in range the final position estimate may be influenced heavily by the DOG appearance model. In one embodiment, the correlation values are normalized prior to multiplication. In further embodiments, the correlation values are weighted such that correlation values of models that match detected video conditions are more heavily weighted.

0039The predicted position from the particle filter gives an approximate location of the track. The dimension of the face (width, height) is obtained from the appearance model. When the face is detected, the width and height provided by the face detector may be used to construct the appearance model. After that, the face is searched across scale space to determine the appropriate width and height. From the predicted position, a search area is constructed

0040If the tracker loses a face then it may be difficult for the tracker to reacquire the face since the estimated position heavily influences the predicted position (since the latest observation is incorporated into the framework). Thus the particle filter framework in one embodiment, maintains multiple distributions of particles. In one embodiment, a boosted particle filter uses AdaBoost (adaptive boost - modifying a classifier in favor of instances misclassified by previous classifiers) detections to form part of the proposal distribution. In particular, based on a parameter, some particles may be sampled from proposal distributions based on AdaBoost detections and some are sampled based on state estimates. This helps the tracker to rely both on state estimates and AdaBoost detections, thus complementing each other. In particular, even though some particles are distributed based on the wrong estimate, some others are still distributed based on the detection (which is assumed correct) and hence the tracker would be able to reacquire the track.

0041In a typical surveillance video, the face may not be detected when it first appears in the scene. This may be due to the quality of the face image such as the size of the face (smaller than the minimum size that can be detected) or the illumination conditions. Hence, the face may be detected only after some frames in the scene. In order to establish the full track of the face, the tracker traces back the face starting from the frame where it was detected and going backwards in time up to the point where it can be reliably tracked. Even though the track is traced only after the face has moved past the frames, the information is used later for better description of where the person was.

0042Backtracking is achieved in a similar manner as normal tracking (which is described before) except that the frames are run backwards. Thus, the additional task can be achieved within the existing framework. Even though the track trace cannot be viewed in real time, it is used for forensic analysis.

0043An example with a realistic surveillance video with multiple people walking and complex backgrounds is illustrated in several frames of a video in <figref idref="f0005">FIG. 7</figref>. Initially, faces in the video are tagged (detected) by means of the face detector, and then the tagged faces are tracked. The output of the tracker is shown for several frames in <figref idref="f0005">FIG. 7</figref>. Two faces are tracked in the sequence of video images. The faces are identified with squares drawn around them. Other attributes may be used to identify the faces, such as highlighting. A person in a dark shirt in Frame 13 is so indicated, along with a person immediately to the left of him in Frame 20. As seen in subsequent frames, the tracker was able to track both people consistently despite various pose changes, including rapid changes over time as they reached the bottom of the stairs and turned, as the person in the dark shirt turned his head in Frame 100. Note that the scale also changed as they walked toward the camera.

0044<figref idref="f0006">FIG. 8</figref> show several frames which illustrate the handling of illumination changes. The illumination changes drastically between Frame 15 and Frame 17, from light to dark. A color template appearance model used alone would lose tracking with such an illumination change. However, the combined models of color template plus DOG appearance model performs better. The person is tagged and tracked over Frames 5 to 18. As the light changes, the track is still maintained.

0045A block diagram of a computer system that executes programming for performing the above algorithms and methods is shown in <figref idref="f0007">FIG. 9</figref>. A general computing device in the form of a computer 910, may include a processing unit 902, memory 904, removable storage 912, and non-removable storage 914. Memory 904 may include volatile memory 906 and non-volatile memory 908. Computer 910 may include - or have access to a computing environment that includes - a variety of computer-readable media, such as volatile memory 906 and non-volatile memory 908, removable storage 912 and non-removable storage 914. Computer storage includes random access memory (RAM), read only memory (ROM), erasable programmable read-only memory (EPROM) &amp; electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, compact disc read-only memory (CD ROM), Digital Versatile Disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium capable of storing computer-readable instructions, as well as data, including video frames. Computer 910 may include or have access to a computing environment that includes input 916, output 918, and a communication connection 920. The computer may operate in a networked environment using a communication connection to connect to one or more remote computers. The remote computer may include a personal computer (PC), server, router, network PC, a peer device or other common network node, or the like. The communication connection may include a Local Area Network (LAN), a Wide Area Network (WAN) or other networks.

0046Computer-readable instructions stored on a computer-readable medium are executable by the processing unit 902 of the computer 910. A hard drive, CD-ROM, and RAM are some examples of articles including a computer-readable medium.

0047A robust face tracking system tracks faces under varying illumination and head pose changes is presented. The system may handle illumination changes by means of combining multiple appearance models. It can also handle both in-plane and out-of-plane pose changes reliably. DOG filters are used to generate appearance models for the face. The appearance models are combined with a color appearance model to obtain robust tracking of a face through illumination, pose and scale changes

0048The Abstract is provided to comply with 37 C.F.R. §1.72(b) to allow the reader to quickly ascertain the nature and gist of the technical disclosure. The Abstract is submitted with the understanding that it will not be used to interpret or limit the scope or meaning of the claims.