Abstract:
A method and system in which facial image representations stored in a database are defined by facial coordinates in a plane common to other images in the database in order to facilitate comparison or likeness of the facial images by comparing the common plane facial coordinates, the common plane being determined by the locations of the eyes and mouth corners; at least one input operatively connected to the at least one processor and configured to input the corners of the eyes and mouth coordinates; the at least one processor configured to convert inputted coordinates for the corners of the eyes and mouth into estimated common plane coordinates by minimizing the error between the inputted corners of the eyes and mouth coordinates and the estimated coordinates corners of the eyes and mouth obtained from the least square estimation model of the common plane coordinates of the corners of eyes and mouth.

Description:
GOVERNMENT INTEREST 
     The invention described herein may be manufactured, used, and/or licensed by or for the United States Government without the payment of royalties 
    
    
     BACKGROUND OF THE INVENTION 
     Face recognition includes the process of recognizing an individual by comparing the captured face image against one or more stored face images to identify a match. The stored images are usually called the gallery or watch list databases. The captured image or video is usually called the probe. 
     In national security or military applications, there is a need for nighttime personnel target identification. At night or in darkness (without illumination), a visible face image (produced using visible light) is of limited use for identification. Thermal face images are found to be more useful at night since the thermal face images can be acquired by thermal sensors without external illumination. The thermal camera measures the heat (temperature) that is emitted (radiated) from the human face. However, most watch list databases contain visible imagery. For a given picture, visible sensors measure the light reflected from the facial surface at a single observation angle. The thermal face images and the visible face images are very different in appearance; but there is a need for cross modality face recognition, i.e., thermal-to-visible face recognition. Thermal-to-visible face recognition would include the identification of a person in a thermal image by comparing the person&#39;s thermal facial image to many visible facial images in a database or watch list. 
     Another challenge in facial recognition occurs during adverse conditions, such as during night time surveillance when frequently the images of the subject are not in a still or frontal position, but instead consist of multiple 3-D head rotations (pose angles). Different pose angles and different illumination conditions can result in different images. For example, the face image of a person in a side view is very different from the image of the same person in a frontal view. It is beneficial to consider the pose changes in face recognition. 
     The research on thermal face recognition is summarized in (1) C. S. Martin et. al., “Recent advances on face recognition using thermal infrared images,” Reviews, Refinements and New Ideas in Face Recognition, InTech Open Access Publisher, edited by Peter M. Corcoran, Chapter 5, pp. 95-112, July 2011; (2) M. K. Bhowmik, et. al., “Thermal infrared face recognition—a biometric identification technique for robust security system,” Reviews, Refinements and New Ideas in Face Recognition, InTech Open Acess Publisher, edited by Peter M. Corcoran, Chapter 6, pp. 114-138, July 2011 ; (3) L. B. Wolff , D. A. Socolinsky, and C. K. Eveland, “Face recognition in the thermal infrared,” Computer Vision Beyond the Visible Spectrum Book, pp. 167-191, Springer London, 2006). Thermal imagery has less texture information. However, thermal imagery is less sensitive to illumination variation and expression changes. It has certain applications for nighttime personnel identification that is crucial to national security and military operations. U.S. Pat. No. 7,406,184, to L. B. Wolf, entitled “Method and apparatus for using thermal infrared for face recognition,” Jul. 29, 2008, discloses a method of incorporating the use of thermal face imagery into an end-to-end face recognition system to mitigate the influence of varying ambient illumination on systems using visible imagery. In this method, both visible face images and thermal face images are used to create a face representation template that is matched or compared with a stored database or gallery of face templates. The system requires both visible imagery and thermal imagery of an individual (wanted individual) in the probe. For a nighttime operation especially without an external illumination, this condition is not easy to meet. 
     Cross-modality face recognition research, in which the gallery contains images in one modality (e.g., visible) and the probe contains images in another modality (e.g., thermal), generally work using the existing techniques that were designed for visible face recognition (such as using a variation of edge information). U.S. Pat. No. 7,512,255, entitled “Multi-modal face recognition,” discloses a method to identify an individual using visible and infrared images using a sequence of multimodal data (3D from multiple visible light cameras, 2D infrared) and employing an Annotated Face Model (AFM), which is a 3D segmented facial model. 
     Generally, there are three common approaches to the problem of pose changes in facial recognition using two dimensional data. First, the training set contains many face images with different poses, which requires multiple datasets. An extension of this type of approach is to train each pose-specific face classifier with multiple examples of that pose. Then, the output of these classifiers is fused to give a face recognizer that can process a wider range of facial poses (for example, see U.S. Pat. No. 7,542,592 to Huang, Pose-invariant Face Recognition System and Process”). Second, for each input image, the estimate of the pose is calculated and then the input image is normalized to a virtual frontal view pose before it is matched against the gallery (see for example, U.S. Pub. Pat. App. No. 2010/0284577, “Pose-variant face recognition using multiscale local descriptors,” Nov. 11, 2010). 
     A third approach is to use multiple 2D images to estimate a 3D face model. The 3D geometry information can be used for pose correction (see U.S. Pub. Pat, Application No. 2010/0149177, “Generation of normalized 2D imagery and IS systems via 2D to 3D lifting of multifeatured objects,” Jun. 17, 2010, M. I. Miller). 
     When using 3D sensor data, a 3D face model can be established, which requires a high-resolution 3D sensor (see U.S. Pub. Pat. Application No. 2006/0120571 A1, “System and method for passive face recognition,” Jun. 8, 2006, P. H. Tu, et. al). 
     Visible sensor measures the reflectance of light from the object. Infrared sensor measures the heat transmitted from the object. Face images that are acquired from visible sensor and infrared sensor represent different face phenomena. Pose rotations of the human head can produce significant changes and distortions in the facial appearance of a person in the face image. Such effects cause degradation of face recognition performance, resulting in possible performance degradation. 
     Sensors of different modalities measure different properties of the object. Visible light and thermal infrared are in different spectrum. Therefore, visible sensors and thermal sensors capture different physical properties of the human face. Because of these physical variations, they are translated into different features. 
     One example is that edges are not aligned to each other in thermal and visible images. Human visual system and man-made algorithms do use edge as the primary information for recognition purposes. However, if e edge information is used to match thermal and visible face images, they do not quite match. 
     Other examples are that when imagery of a person is acquired from different sensors at different time, they do not have the same 3D orientation angles. This is a very common phenomenon. In a practical situation, the camera might not be exactly looking straightly at the face when the image is acquired. If one wants to match face images with different poses, one needs to address the 3D transformation among them first. 
     There exists a need to identify a personnel target when the probe image and gallery image are in different modalities. There exists a need to utilize common structures. There also exists a need to (biometric landmarks) characterize both images from gallery and probe, e.g., visible and thermal images. 
     It is known to use biometric landmarks such as the eyes and mouth to capture the predominate identifying features of an image. For example, Intel® Perceptual Computing SDK—How to use the Face Detaction Software, https://software.intel.com/en-us, Oct. 30, 2012, discloses The Intel® Perceptual Computing SDK, which is a library of pattern detection and recognition algorithm implementations exposed through standardized interfaces. The SDK provides a suite of face analysis algorithms including face location detection, landmark detection, face recognition and face attribute detection. The Intel face detection algorithm locates the rectangle position of a face or multiple faces from an image or a video sequence in real-time capture or playback mode. The detection algorithm locates the 6 point or 7 point landmarks namely, the outer and inner corners of the eyes, the tip of the noise, and the outer corners of the mouth. 
     Since the face images do not have the same 3D orientation angles when they are acquired from different sensors and/or at different times, there is a need to develop a new 3D registration method via a single frame for face images having 3D pose angles. 
     It is difficult for a machine to search facial photos for a match in part because when a person&#39;s photo appears in a non-standard position, the orientation of the face is different; i.e., the faces are generally pointed at an angle to the vertical plane of the camera. Because the orientation of the face is different, the coordinates of the facial features do not match those of a pose position in which the subject&#39;s head is upright and aligned with the camera (commonly referred to as a mug shot). 
     SUMMARY OF THE INVENTION 
     A preferred embodiment and method of the present invention provides a method and system for face recognition when face images in the gallery and probe are in different modalities and/or have different 3D pose angles due to the fact that they are acquired by different sensors and/or at different times. A preferred embodiment converts the coordinates of the eyes and mouth of a facial image appearing in a randomly oriented photograph to virtual coordinates corresponding to an estimate of how the head and facial coordinates would appear if the subject&#39;s head were to be turned such that the centers of the eyes and corner of the mouth were oriented in a vertical plane; i.e., the roll, pitch, and yaw are zero and the scale is one. The estimate is obtained using models based upon a non-linear Gaussian Least Square Differential Correction. The matching is performed based on matching the virtual coordinates of facial images in the probe against ones in the gallery. 
     The present invention provides a new system to mitigate degradation sources due to cross modality. It utilizes common biometric landmarks characterizing both visible and thermal images. Since the edge information and texture information do not correspond to each other in visible and thermal imagery, only common features remain for visible and thermal images of the same person. They are biometric landmarks, such as eye locations, mouth location, etc. The present invention also provides a new 3D registration method via a single frame for face images having 3D pose angles. 
     A preferred embodiment of the present invention comprises a system for facial recognition using an image of a human in which the eyes and mouth corners are not in a vertical plane when the image was taken to determine estimated virtual plane coordinates corresponding to estimated location of eyes and mouth corners of the human when the eyes and mouth corners are in a vertical plane; the system comprising: 
     at least one processor configured to determine the virtual plane coordinates; 
     at least one input operatively connected to the at least one processor and configured to input the first corners of the eyes and mouth coordinates F 1 , F 2 , F 3 , F 4 , F 5 , F 6 , F 7 , F 8 , F 9 , F 10 , F 11 , F 12  from the image into the at least one processor comprising data points in a vector
 
 F =( F   1   , . . . ,F   12 ) t  
 
where t represents the transpose, F 1  represents the horizontal coordinate of the left eye outer corner, F 2  represents the vertical coordinate of the left eye outer corner, F 3  represents the horizontal coordinate of the left eye inner corner, F 4  represents the vertical coordinate of the left eye inner corner, F 5  represents the horizontal coordinate of the right eye outer corner, F 6  represents the vertical coordinate of the right eye outer corner, F 7  represents the horizontal coordinate of the right eye inner corner, F 8  represents the vertical coordinate of the right eye inner corner, F 9  represents the horizontal coordinate of the left mouth corner, F 10  represents the vertical coordinate of the left mouth corner, F 11  represents the horizontal coordinate of the right mouth corner, and F 12  represents the vertical coordinate of the right mouth corner;
 
     the at least one processor configured to convert the first coordinates F 1 , F 2 , F 3 , F 4 , F 5 , F 6 , F 7 , F 8 , F 9 , F 10 , F 11 , F 12  into second coordinates for the corners of the eye and mouth in a virtual vertical plane comprising P 1 =(−x e1 ,y e1 ,0), P 2 =(−x e2 ,y e2 ,0), P 3 =(x e1 ,y e1 ,0), P 4 =(x e2 ,y e2 ,0), P 5 =(−x m ,y m ,0), and P 6 =(x m ,y m ,0), where P 1  is the estimated left eye outer corner coordinates, P 2  is, the estimated left eye inner corner coordinates, P 3  is the estimated right eye outer corner coordinates, P 4  is the estimated right eye inner corner coordinates, P 5  is the estimated left mouth corner coordinates and P 6  is the estimated right mouth corner coordinates and x and y represent horizontal and vertical distances from a facial reference point, and to determine the head orientation of the human subject using roll, yaw and pitch relative to the virtual vertical plane where, θ represents the yaw, φ represents the pitch, and Ψ represents the roll; 
     the at least one processor configured to solve for the parameter vector V p  comprising 9 parameters x e1 , y e1 , x e2 , y e2 , x m , y m , θ, φ, ψ using the following equation:
 
 V   p =( x   e1   ,y   e1   ,x   e2   ,y   e2   ,x   m   ,y   m ,θ,φ,ψ) t  
 
where x e1 , y e1 , x e2 , y e2 , x m , y m , represent the virtual plane coordinates of the corners of the eyes and mouth in the vertical virtual plane, and wherein the error to be minimized is the error between the inputted corners of the eyes and mouth coordinates F 1 , F 2 , F 3 , F 4 , F 5 , F 6 , F 7 , F 8 , F 9 , F 10 , F 11 , F 12  and the least square estimation model function ƒ(V p ) coordinates of the estimated inputted corners of the eyes and mouth coordinate values comprising horizontal coordinates f 1 (V p ), f 3 (V p ), f 5 (V p ), f 7 (V p ) f 9 (V p ), f 11 (V p ), and vertical coordinates ƒ 2 (V p ), ƒ 4 (V p ), ƒ 6 (V p ), ƒ 8 (V p ), ƒ 10 (V p ), ƒ 12 (V p ) and where the least square model function ƒ(V p ) is computed using
 
     f 1 (V p )=[1,0,0]T P 1    
     f 2 (V p )=[0,1,0]T P 1    
     f 3 (V p )=[1,0,0]T P 2    
     f 4 (V p )=[0,1,0]T P 2    
     f 5 (V p )=[1,0,0]T P 3    
     f 6 V p )=[0,1,0]T P 3    
     f 7 (V p )=[1,0,0]T P 4    
     f 8 (V p )=[0,1,0]T P 4    
     f 9 (V p )=[1,0,0]T P 5    
     f 10 (V p )=[0,1,0]T P 5    
     f 11 (V p )=[1,0,0]T P 6    
     f 12 (V p )=[0,1,0]T P 6    
     where T correlates to the head orientation θ, φ, ψ and the matrix
 
 T=T   θ   ·T   φ   ·T   ψ 
 
     where 
               T   θ     =     [           cos   ⁢           ⁢   θ           sin   ⁢           ⁢   θ         0               -   sin     ⁢           ⁢   θ           cos   ⁢           ⁢   θ         0           0       0       1         ]                   T   ϕ     =     [           cos   ⁢           ⁢   ϕ         0         sin   ⁢           ⁢   θ             0       1       0               -   sin     ⁢           ⁢   ϕ         0         cos   ⁢           ⁢   ϕ           ]                   T   ψ     =       [         1       0       0           0         cos   ⁢           ⁢   ψ           sin   ⁢           ⁢   ψ             0           -   sin     ⁢           ⁢   ψ           cos   ⁢           ⁢   ψ           ]     .           
whereby the virtual plane coordinates are used to search a database of a plurality of images of people represented by either virtual plane or vertical plane coordinates in order to identify or recognize an individual.
 
     These and other aspects of the present invention will be described in more detail below in conjunction with the following drawings. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
       The present invention can best be understood when reading the following specification with reference to the accompanying drawings, which are incorporated in and form a part of the specification, illustrate alternate embodiments of the present invention, and together with the description, serve to explain the principles of the invention. In the drawings: 
         FIG. 1  is an illustration of a flow chart providing a partial overview of the present invention including a description of a vector comprising ten parameters which include, inter alia, yaw, pitch and roll of the subject&#39;s pose. 
         FIG. 2  is an illustration of a comparison between prior art coordinate systems and the coordinate system in accordance with the principles of the present invention and part of the conversion procedure. 
         FIG. 3  is an illustration which continues from  FIG. 2  and illustrates, inter alia, the conversion and/or transposition of coordinates which take into account, inter alia, yaw, pitch and roll. 
         FIG. 4  is an illustration showing a flow chart illustrating the solving for unknown parameters using the Gaussian Least Squares Differential Correction (GLSDC). 
         FIG. 5  is an illustration showing a flow chart illustrating the solving for unknown parameters using the Gaussian Least Squares Differential Correction (GLSDC). 
         FIG. 6  is schematic representation of modules I, II and II. 
         FIG. 7  is an overall schematic representation of a preferred embodiment of the present invention for identifying a personnel target when the probe image and gallery image are in different modalities and face images having 3D pose angles due to the fact that they are acquired from different sensors and/or at different times. 
         FIG. 8  is a diagram for solving a transformation parameter vector. 
         FIG. 9  illustrates a graphic view showing the stereotactic registration. 
         FIG. 10A  is an illustration showing the pose of a person and the rotation directions for yaw, roll and tilt. 
         FIG. 10B  is an illustration showing an image onto which the facial plane is marked in order to show the roll of the face from the center of the facial plane 
         FIG. 11  is a schematic representation of a mask labelled with the fiducial points or facial features of the eye corners and mouth cornets. Lines L1 and L2 define an origin (at the intersection). 
         FIG. 12  illustrates a graphic view showing the stereotactic-based matching that is formulated as a multiple hypotheses testing. 
         FIG. 13  illustrates a graphic view showing the procedure of simulating varying 3D posing angles. 
         FIG. 14  shows an example where facial points of three real faces are displayed in the UCS. 
         FIG. 15  shows face images with various pose views. 
         FIG. 16  shows the results of simulating 3D pose angles between −10 to 10 degrees. 
         FIG. 17  shows the results of simulating 3D pose angles between −20 to 20 degrees. 
         FIG. 18  is an illustration schematically illustrating that matching can be implemented by exploring multiple frames through temporal information. 
         FIG. 19  shows when 10 frames are used, the matching error is reduced to 0.74 from the error range between 5 and 10. 
         FIG. 20  illustrates matching result of one thermal probe face to two visible gallery faces. 
         FIG. 21A  is a flow chart of the methodology of a preferred embodiment of the present invention. 
         FIG. 21B  is a continuing flow chart of the methodology of a preferred embodiment of the present invention. 
         FIG. 22  depicts a high level block schematic diagram of a general purpose computer. 
     
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
     The present invention now will be described more fully hereinafter with reference to the accompanying drawings, in which embodiments of the invention are shown. However, this invention should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art. In the drawings, the thickness of layers and regions may be exaggerated for clarity. Like numbers refer to like elements throughout. As used herein the term “and/or” includes any and all combinations of one or more of the associated listed items. 
     The terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the full scope of the invention. As used herein, the singular forms “a”, “an” and “the” are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms “comprises” and/or “comprising,” when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. 
     Unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the relevant art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein. 
     A flow chart of a preferred embodiment of the present invention is shown in  FIG. 6 , which contains 3 modules. The thermal test image and visible reference image are input to the Module I, which estimates coordinates of features for both thermal test image and visible reference image. These are biometric landmarks of the human face, such as eyes, nose, mouth, ears, etc. Many available algorithms have been researched (see (1) U.S. Pub. Pat. App. No. 2012/0308124, “Method and system for localizing parts of an object in an image for computer vision applications,” Dec. 6, 2012, P. N. Belhumeur; (2) European Pat. No. 1 811 456 B1, “Face feature point detector and feature point detector,” Nov. 12, 2004, K. Kinoshita; (3) P. Viola and M. Jones, “Robust real-time object detection,” Second International Workshop On Statistical And Computational Theories Of Vision—Modeling, Learning, Computing, And Sampling, pp. 1-25, Vancouver, Canada, Jul. 13, 2001). Visible and thermal images are acquired from different sensors. They need to be registered first. Module II ( 45 ) performs the stereotactic registration based on those features, or sub-set of those features. For example, when the ear features are available, the algorithm uses features including ears. When the ear features are not available, the algorithm uses the sub-set of features that exclude ears. The algorithm takes available and reliable features and performs the constraint optimization using those features. The registration transformation includes shifts, 3D rotations, and scale. Module III ( 46 ) performs the matching based on multiple hypotheses testing theory, which uses the features other than edges to correlate the visible and thermal face images. 
     Pictures of subjects are taken with their faces not pointed directly at the camera so that the face images may contain yaw (pose rotation left and right), pitch (head tile rotation (up and down), and roll (in plane rotation). Also, pictures of a subject may be taken using a different modality (such as a thermal camera instead of a conventional camera). Accordingly, facial recognition in pictures are taken with the head of a subject not pointed directly at the camera may be effected by yaw (pose rotation left and right) pitch (head tilt up and down), and roll (in plane rotation). In accordance with a preferred embodiment of the present invention, yaw, pitch and roll are incorporated into a Unified Coordinate system (UCS). The Unified Coordinate System includes yaw, pitch and roll of the subject&#39;s head and defines by measures the location of the inner and outer points of the two eyes and inner and outer points of the mouth, Optionally, a preferred embodiment comprises the development of a parameter vector V p  which comprises the six coordinates of the eyes and mouth, the yaw, pitch and roll of the subject&#39;s head and the scale (size) of the face. The parameter vector is defined as V p =(x e1 ,y e1 ,x e2 ,y e2 ,x m ,y m ,θ,φ,ψ,α) t  where x e1 , y e1 , x e2 , y e2 , x m , &amp; y m  are the coordinates used to define the previously mentioned six points of defining the eyes and mouth and θ defines the yaw, φ defines the pitch and Ψ defines the roll and α defines the scale. 
     As illustrated in  FIG. 1 , pictures of subjects are taken with their faces not pointed directly at the camera so that face images may contain yaw (pose rotation left and right), pitch (head tile rotation (up and down), and roll (in plane rotation). Also, pictures of a subject may be taken using a different modality (such as a thermal camera instead of a conventional camera). Accordingly, facial recognition in pictures are taken with the head of a subject not pointed directly at the camera may be effected by yaw (pose rotation left and right) pitch (head tilt up and down), and roll (in plane rotation). In accordance with a preferred embodiment of the present invention, yaw, pitch and roll are incorporated into a Unified Coordinate system (UCS). The Unified Coordinate System includes yaw, pitch and roll of the subject&#39;s head and defines by measures the location of the inner and outer points of the two eyes and inner and outer points of the mouth, Optionally, a preferred embodiment comprises the development of a parameter vector V p  which comprises the six coordinates of the eyes and mouth, the yaw, pitch and roll of the subject&#39;s head and the scale (size) of the face. The parameter vector is defined as V p =(x e1 ,y e1 ,x e2 ,y e2 ,x m ,y m ,θ,φ,ψ,α) T  where x e1 , y e1 , x e2 , y e2 , x m , &amp; y m  are the coordinates used to define the previously mentioned six points of defining the eyes and mouth and θ defines the yaw, φ defines the pitch and Ψ defines the roll and α defines the scale. 
     Referring now to  FIG. 2 , shown are two facial images  10 A and  10 B. Facial image  10 A overlays a grid using coordinates u, v in the world coordinate system and Face image  10 B overlays a grid using the Unified Coordinate system of the present invention. For the image  10 A, the biometric features of the outer and inner corners of the left eye,  16 A,  17 A, outer and inner of the right eye,  18 A,  19 A, and the corners of the mouth  20 A,  21 A are shown. Similarly, for the image  10 B, the biometric features of the outer and inner corners of the left eye,  16 B,  17 B, outer and inner of the right eye,  18 B,  19 B, and the corners of the mouth  20 B,  21 B are shown. Image  10 A is a conventional image taken with unknown values relating to yaw, pitch or roll. The image  10 B is taken with the face in a pose such that the eyes and the mouth are in the same vertical plane. For the purposes of the Unified Coordinate System (UCS) the origin point 0 (0,0,0) of the face plane is defined as the intersection of two lines L 1  and L 2 , as shown in  FIG. 2 . The line L 1  connects the center of two corners of left eye and the right corner of mouth. The line L 2  connects the center of two corners of right eye and the left corner of mouth. In order to have this plane that passes through these four points being vertical, the human head is tilted from the normal frontal pose. These lines form the central or face plane of the image  10 B from which the yaw, pitch and roll are to be measured. By definition, the yaw, pitch and roll is zero when the center of the eyes and corners of the lips are in the same vertical plane. The intersection of the lines further defines the origin (0,0,0) of the Unified Coordinate System developed in accordance with the principles of the present invention. The Unified Coordinate System defines the locations of the eyes and mouth on the basis of the distance from the middle of the face in the “x” direction and the distance from the origin. 
     Since the face is generally symmetrical in the lateral or “x” direction, measurements for the six points P1 through P6. The left eye outer corner is represented by P 1 =(−x e1 ,y e1 ,0). The left eye inner corner is represented by P 2 =(−x e2 ,y e2 ,0). The right eye outer corner is represented by P 3 =(x e1 ,y e1 ,0). The right eye inner corner is represented by P 4 =(x e2 ,y e2 ,0). The left mouth corner is represented by P 5 =(−x m ,y m ,0) and the right mouth corner is P 6 =(x m ,y m ,0). Due to symmetry, the 6 points may be represented by six parameters X e1 , X e2 , Y e1 , Y e2 , X m  and Y m , where x e1  and y e1  represent the outer corners of two eyes, x e2  and y e2  represent the inner corners of two eyes, x m  and y m  represent the corners of mouth, which form the first six parameters for the parameter vector V p . These 6 parameters define the face plane. The remaining parameters represent the yaw θ, pitch Φ, roll Ψ, and scale. 
     Referring now to the conventional image  10 A of  FIG. 2 , the locations for the biometric features of the outer and inner corners of the left eye,  16 A,  17 A, outer and inner of the right eye,  18 A,  19 A, and the corners of the mouth  20 A,  21 A may be determined by measuring the distance from the U and V axes. However, the yaw θ, pitch Φ, roll Ψ are unknown. 
     Facial images points  10 A are measurements in the conventional system with the origin located to the left of and below the facial image  10 A as depicted in  FIG. 2 . To convert coordinates  10 A into coordinates of the Unified Coordinate, the first step, depicted in Box  32  in  FIG. 4  is to determine the average coordinate positions in the horizontal and vertical positions as shown in Box  32 . The center of mass equations are utilized as U M =ΣU i /6 where i varies from 1 to 6 and V M =ΣV i /6 (see Box  32 ). The center of mass U M  and V M  is then subtracted from the U and V coordinates to yield the corrected Ū as shown in Boxes  33  and  34 . 
     Box  22  of  FIG. 2  is directed to the conversion of the points  16 A,  17 A,  18 A,  19 A,  20 A, and  21 A (defined as measurements using the U, V, coordinates into the coordinates P1 through P6 of the Unified Coordinate System (UCS) used with the preferred embodiment of the present invention. Note that the pointes  16 A through  21 A do not take into account the yaw, pitch and roll, which is done using the following methodology. 
     In order to convert the measurement of the Image  10 A (points  16 A- 21 A) into the Unified Coordinate System, the yaw, pitch and roll must be taken into account. For each landmark point (corners of eyes and mouth) on a face, the transformation equation is written as,
 
 P   Mn   =α·T·P   n   (1)
 
where P Mn =[u,v,r] t  represents the coordinates of a landmark of an acquired face image in the measurement (image) domain, P n =[x,y,z] t  represents the coordinates of a landmark in the UCS domain, αis the scale factor, and t is the transformation matrix. This transformation matrix is written as,
 
 T=T   si   ·T   phi   ·T   te   (2)
 
where T te  represents the in-lane rotation transformation in the (x,y) plane, that is,
 
                       T   te     =     [           cos   ⁢           ⁢   θ           sin   ⁢           ⁢   θ         0               -   sin     ⁢           ⁢   θ           cos   ⁢           ⁢   θ         0           0       0       1         ]       ,           (   3   )               
T phi  represents the pose-rotation transformation in the (x,z) plane, that is,
 
                       T   phi     =     [           cos   ⁢           ⁢   ϕ         0         sin   ⁢           ⁢   θ             0       1       0               -   sin     ⁢           ⁢   ϕ         0         cos   ⁢           ⁢   ϕ           ]       ,           (   4   )               
T si  represents the tilt-rotation transformation in the (y,z) plane, that is,
 
                       T   si     =     [         1       0       0           0         cos   ⁢           ⁢   ψ           sin   ⁢           ⁢   ψ             0           -   sin     ⁢           ⁢   ψ           cos   ⁢           ⁢   ψ           ]       ,           (   5   )               
where θ, φ, and ψ are the angles of in-plane, pose, and tilt rotations, respectively, as shown in  FIGS. 2, 10A, 10B .
 
     For a 2-D face image, the coordinates of each landmark result in z=0 for the UCS domain. For the coordinates of each landmark in the measurement domain, the r is not available. The measureable coordinates of a landmark in the measurement domain are denoted as P Mn,2 =[u,v] t , t represents the transpose operator. 
     From equation (1), in order to estimate the transformation that transforms a face image into the UCS, ten parameters are estimated in the UCS domain, which are defined as a parameter vector:
 
 V   p =( x   e1   ,y   e1   ,x   e2   ,y   e2   ,x   m   ,y   m ,θ,φ,ψ,α) t .  (6)
 
     The coordinates of each landmark in a face image can be represented by this parameter vector, V p , as shown in (1). The coordinates of each face landmark can also be obtained by measurements from the face image, which is denoted as  P   Mn,2 =[ū, v ] t , n=1,2, . . . ,6. 
     For a face image, 6 landmarks are used in the measurement domain which results 6 pairs of (u,v) values. Therefore, 12 equations can be formed from equation (1) to estimate 10 parameters in (6). 
     The aim is to make the difference between the measurement coordinates and the estimated coordinates of face landmarks minimal. That is,
 
min|   P     Mn,2   −P   Mn,2 |.  (7)
 
     In the following, minimization problem is formulated into a cost function via a procedure that we developed in S. S Young, “Optimized target localization in stereotactic radiosurgery using real-time digital portal images,” Phys. Med. Biol. Vol. 41, pp. 1621-1632, 1996. In this way, the estimate of ten parameters minimizes the cost function. 
     Define V p  as the parameter vector,
 
 V   p =( x   e1   ,y   e1   ,x   e2   ,y   e2   ,x   m   ,y   m ,θ,φ,ψ,α) t .  (8)
 
     The function F, which is the function of the parameter vector V p , can be defined as
 
 F= ƒ( V   p )  (9)
 
where
 
 F= ( F   1   , . . . ,F   12 ) t   (10)
 
ƒ( V   p )=(ƒ 1 ( V   p ), . . . ,ƒ 12 ( V   p )) t   (11)
 
                     F   k     =     {               F       2   ⁢   n     -   1       =           P   _       Mn   ,   2       ⁡     (   1   )       =     u   _         ,           k   =     odd   ⁢           ⁢   value                     F     2   ⁢   n       =           P   _       Mn   ,   2       ⁡     (   2   )       =     v   _         ,           k   =     even   ⁢           ⁢   value                       (   12   )                   f   k     ⁡     (     V   p     )       =     {                 f       2   ⁢   n     -   1       ⁡     (     V   p     )       =         P     Mn   ,   2       ⁡     (   1   )       =   u       ,           k   =     odd   ⁢           ⁢   value                       f     2   ⁢   n       ⁡     (     V   p     )       =         P     Mn   ,   2       ⁡     (   2   )       =   v       ,           k   =     even   ⁢           ⁢   value                       (   13   )               
where P Mn,2 (1)=u and P Mn,2 (2)=v are calculated from equation (1) using the parameter vector V p , and n=1, . . . ,6 represent 6 landmarks and k=1, . . . ,12 represent 12 equations that are resulted from 6 landmark measurements. The measurements  P   Mn,2 (1)=ū and  P   Mn,2 (1)=ū may be obtained from the inquiring face image by a user mouse interface or by using any of the facial landmark detection algorithms.
 
The goal is to find an estimate of V p , denoted as {tilde over (V)} p .This estimate minimizes
 
| F− ƒ( {tilde over (V)}   p )|.  (14)
 
To estimate ten parameters V p =(x e1 ,y e1 ,x e2 ,y e2 ,x m ,y m ,θ,φ,ψ,α) t , a non-linear least square algorithm, such as the Gaussian least square differential correction (GLSDC) can be applied. The function F can be approximated by
 
                     Δ   ⁢           ⁢   F     ≈       (       ∂   f       ∂     V   p         )     ⁢   Δ   ⁢           ⁢       V   p     .               (   15   )               
If ΔF is given, then ΔV p  can be found as follows:
 
                     Δ   ⁢           ⁢     V   p       =         [         (       ∂   f       ∂     V   p         )     t     ⁢     (       ∂   f       ∂     V   p         )       ]       -   1       ⁢       (       ∂   f       ∂     V   p         )     t     ⁢   Δ   ⁢           ⁢     F   .               (   16   )               
The algorithm can be summarized in the following.  FIG. 8  illustrates a schematic representation of this algorithm, which may comprise:
 
Step (1) Input measurements F (Box  200 ).
 
Step (2) Initialize the parameter vector (Box  210 )
 
 V   p =( x   e1   ,y   e1   ,x   e2   ,y   e2   ,x   m   ,y   m ,θ,φ,ψ,α) t .  (17)
 
Step (3) Calculate the estimated coordinate P Mn,2 =ƒ(V p ) using (13). (Box  220 )
 
Step (4) Calculate ΔF=F−ƒ(V p ) (Box  230 ).
 
Step (5) Calculate the estimate correction vector ΔV p  using (16) (Box  240 ).
 
Step (6) Update the parameter vector V p =V p +ΔV p  (Box  250 ).
 
Step (7) Stop iteration if ΔF does not change significantly from one iteration to another (Box  260 ), otherwise go to step (4).
 
Step ( 8 ) Output the estimated parameter vector V p  (Box  270 ).
 
     Referring now to  FIG. 6 , a flow chart of the overall preferred embodiment system is shown, which contains 3 modules ( 44 A,  44 B),  45  and  46 . The thermal test image and visible reference image are input to the Module I (Boxes  45 A,  45 B). Module I estimates coordinates of features for both thermal test image and visible reference image. These are biometric landmarks of the human face, such as eyes, nose, mouth, ears, etc. Visible and thermal images are acquired from different sensors, and are first registered. Module II (Box  45 ) performs registration based on those features, or sub-set of those features. For example, when the ear features are available, the algorithm uses features including ears. When the ear features are not available, the algorithm uses the sub-set of features that exclude ears. The algorithm takes available and reliable features and performs the constraint optimization using those features. The registration transformation includes shifts, 3D rotations, and scale. Module III (Box  46 ) performs the matching based on multiple hypotheses testing theory, which uses the features other than edges to correlate the visible and thermal face images 
       FIG. 7  shows another preferred embodiment of the invention. Assume the gallery contains multiple 2D face images. Each of them is input to algorithm of the estimate coordinates of face features (Box  150 ). After this, the face features corresponding to each face are input to the stereotactic registration (Box  152 ). In the step of stereotactic registration, the transformation parameters of features of each face are estimated via the generalized 3D registration algorithm (Box  154 ). Then the features of each face are transformed to the common coordinate system, the unified coordinates system (UCS) (Box  156 ), which is defined in details in the section of Stereotactic Registration. The 2D face image in the probe is also input to the algorithm of the estimate coordinates of face features (Box  150 ). Then, the features of the probe face are passed into the stereotactic registration step (Box  152 ). The probe face features are also transformed into the UCS domain (Box  156 ). Now, both transformed gallery face features and probe face features are input to the multiple hypotheses testing matching (Box  160 ). The matching algorithm generates the matching scores of the probe face with each face in the gallery. The minimum matching score corresponds to the correct match. 
     Many algorithms are available for finding available coordinates of face features. In this invention, the stereotactic registration (Box  45 ) and the stereotactic-based matching (Box  46 ) are described in the following. 
     Stereotactic Registration 
     When a face image is acquired under an un-cooperative condition, the face image can contain three rotations, such as, in-plane, pose, and tile rotations as shown in  FIG. 10A . The common terms of roll, jaw, and pitch can be used to describe these three rotations as illustrated in  FIG. 10A . 
     In order to match face images that are taken by different sensors and at different times, the face images from gallery and probe need to be transformed into a common coordinate system. This common coordinate system should be invariant to the sensor position with respect to the human subject when the face image is taken. In another word, this common coordinate system should be invariant of relative pose of the human head. 
     The following describes, a method of a generalized yaw, pitch, roll, scale, and shift transformation to a unified coordinate system (UCS) via a single frame. 
     Unified Coordinate System (UCS) 
     Before describing the UCS, a human face is described in three dimensions. Looking at a side view of a human face, there is a depth from eyes to mouth. This depth is unique for each person. Therefore, for a normal posed human head, which is normally termed as a frontal pose, the centers of corners of each eye and the corners of mouth are on two different vertical planes. 
     Using a unified coordinate system (UCS), four points on a human face are labelled. 
     This face plane is different from the Face Plane that was used in U.S. Pat. No. 7,221,809, “Face recognition system and method,” May 22, 2007, Z. J. Geng. The face plane that is used in the conjunction with a preferred embodiment of the present invention is the plane that passes through centers of two eyes and outer corners of mouth and is vertical. The Face Plane that was used in U.S. Pat. No. 7,221,809 was the plane that passed through centers of two eyes and outer corners of mouth but was not defined as vertical. And it also used the location of the nose tip to derive the pan (title) angle. Usually, it is difficult to locate the nose tip accurately. In accordance with the current invention, use of the nose tip is avoided. 
     This Unified Coordinate System nomenclature is unique for each person. That is, it is invariant in coordinates of the sensor with respect to the person, and relative pose of the human head. A face image is transformed from any other position into the UCS. Then the face recognition or face matching is performed in this UCS. The matching could be between any desired features. 
     Since face images at any other positions with possible three rotation angles are transformed into the Unified Coordinate System, the face recognition problem can be addressed with un-cooperative or uncontrolled conditions where face images have 3D rotational angles; allowing the address un-cooperative face recognition based on one single 2D frame via performing the generalized yaw, pitch, roll, scale, and shift transformation to a unified coordinate system (UCS). 
     Shift is related to the relative position of the object to the center of the image plane. Scale is related to the object range and sensor focal plane properties (number of pixels, element spacing, etc.) Yaw, pitch, and roll are related to the pose rotation, head tilt rotation and in plane rotation. 
     Generalized Yaw, Pitch, Roll, Scale, and Shift Transformation to a Unified Coordinate System (UCS) via a Single Frame 
     The face plane under the UCS is illustrated in  FIG. 10B . Assume a face is symmetric with respect to a center line, L c , as shown in  FIG. 10B . 
     In the face plane of the Unified Coordinate, the eye corners and mouth corners are in a same vertical plane where z=0. There are 6 primary points, or landmarks, to describe a face; i.e., the outer corners of two eyes, inner corners of two eyes, and corners of mouth. 
     The origin point O (0,0,0) of the face plane is defined as the intersection of two lines L 1  and L 2 , as shown in  FIG. 11 . The line L 1  connects the center of two corners of left eye and the right corner of mouth. The line L 2  connects the center of two corners of right eye and the left corner of mouth. 
     Since these 6 primary points are symmetric with respect to the center line L c , there are 6 parameters to define the face plane. They are (x e1 ,y e1 ,x e2 ,y e2 ,x m ,y m ), where x e1  and y e1  represent the outer corners of two eyes, x e2  and y e2  represent the inner corners of two eyes, x m  and y m  represent the corners of mouth, as shown in  FIG. 11 . 
     Therefore, these 6 primary points are presented in a three dimensional space as follows. 
     The left eye outer corner is represented by P 1 =(−x e1 ,y e1 ,0). The left eye inner corner is represented by P 2 =(−x e2 ,y e2 ,0). The right eye outer corner is represented by P 3 =(x e1 ,y e1 ,0). The right eye inner corner is represented by P 4 =(x e2 ,y e2 ,0). The left mouth corner is represented by P 5 =(−x m ,y m ,0) and the right mouth corner is P 6 =(x m ,y m ,0). These parameters are obtained via a minimization procedure by solving the equation (7). 
     Stereotactic-Based Matching 
     After the facial features are transformed into the UCS, matching can be performed. In the approach of stereotactic-based matching, the matching problem is formulated as a multiple hypotheses testing as illustrated in  FIG. 12 . 
     Multiple Hypotheses Testing 
     In this multiple hypotheses testing formation, we will present the gallery model, the probe model, and the decision rule. 
     Gallery Model 
     As shown in  FIG. 12 , the gallery contains Mface images that belong to Msubjects. In the gallery model, the gallery is denoted by a set of Mhypotheses,
 
 H   m   , m= 1,2 , . . . ,M   (18)
 
with a prior probability
 
 P   m =Prob(Gallery generating  H   m )  (19)
 
     Probe Model: 
     The probe contains the test face image. The probe observation vector, R, contains N landmark coordinates (r xi ,r yi ), that is,
 
 {right arrow over (R)}=[r   x1   r   y1   , . . . r   xN   r   yN ]  (20)
 
The coordinates of landmarks of the gallery are denoted as (x mi ,y mi ) as shown in  FIG. 12 . The coordinates from the probe image and the gallery image are related by some noise, this noise is assumed under the added white Gaussian noise model, called AWGN model. Therefore, it follows that
 
 r   xi   =x   mi   +n   xi   , r   yi   =y   mi   +n   yi   , i= 1,2 , . . . N   (21)
 
where n xi &#39;s and n yi &#39;s are independent identically distributed (iid) normal or Gaussian random variables with zero mean and variance σ n   2 , respectively, that is,
 
 E ( n   xi )=0 , E ( n   xi   2 )=σ n   2   , E ( n   yi )=0 , E ( n   yi   2 )=σ n   2 .  (22)
 
Then, the conditional probability in which the probability of the observation vector given the gallery generating H m  is represented as follows,
 
                   DECISION   ⁢           ⁢   RULE   ⁢       :     ⁢     
     ⁢           [       P   (       R   ⇀     /     H   m       ]     =       1       (         2   ⁢           ⁢   π       ⁢     σ   n       )       2   ⁢   N         ⁢       exp   [         -       ∑     i   =   1     N     ⁢       (       r   xi     -     x     m   ⁢           ⁢   i         )     2         +       (       r   yi     -     y     m   ⁢           ⁢   i         )     2         2   ⁢     σ   n   2         ]     .                   (   23   )               
Under the above model, the multiple hypotheses testing theory allows the development of the decision rule as shown in the following:
 
                       p   ⁡     (       R   ⇀     /     H   k       )         &lt;     not   ⁢           ⁢     H   k           &gt;     not   ⁢           ⁢     H   m           ⁢     p   ⁡     (       R   ⇀     /     H   m       )               (   24   )               
or
 
                     max   m     ⁢       [     p   ⁡     (       R   ⇀     /     H   m       )       ]     .             (   25   )               
By manipulating the equation in (23), it follows that
 
                       min   m     ⁢       ∑     i   =   1     N     ⁢       (       r   xi     -     x     m   ⁢           ⁢   i         )     2         +         (       r   yi     -     y     m   ⁢           ⁢   i         )     2     .             (   26   )               
In summary, the decision rule is to calculate the minimum value according (26) among Msubjects in the gallery. The subject m that results the minimum value in equation (26) is claimed as the right match to the test subject. The theory of multiple hypotheses testing shows that the cost function of this decision rule in equation (26) is optimal to minimize probability of error under the AWGN measurement model.
 
     EXAMPLE 
       FIG. 13  illustrates a simulation of varying 3D posing angles. The procedure is shown in this figure. The facial landmarks are obtained from several faces. These facial landmark points are input to a transformation with roll, pitch, yaw and scale. The outputs are our measured facial points. Now, these measured points pass into our generalized 3D estimation algorithm to transform to the unified coordinate system (UCS).  FIG. 14  shows an example where facial points of three real faces are displayed in the unified coordinate system (UCS). 
       FIG. 15  shows face images with various pose views. It illustrates that two corners of eyes are not visible beyond 30° pose view.  FIGS. 16 and 17  show the results of simulating 3D pose angles that are below 30 degrees. 
       FIG. 16  shows the results of simulating 3D pose angles between −10 to 10 degrees. 20 random angles between −10 to 10 degrees were chosen to obtain the face with 20 different poses. The first 20 points (from 1 to 20) on the horizontal axis represent 20 cases of face 1. The next 20 points (from 21 to 40) represent 20 cases of face 2. From 41 to 60 are 20 cases of face 3. The vertical axis is the matching error score. The small score represents the correct matching. From  FIG. 16 , it can be seen that, even with different pose angles, the scores of matching the face with its own face are much smaller than ones of the face matching with other faces. 
       FIG. 17  shows the results of simulating 3D pose angles between −20 to 20 degrees. 100 cases were attempted for each face. On the horizontal axis, 1 to 100 are for face 1; 101 to 200 are for face 2; 201 to 300 are for face 3. Again, the results in  FIG. 17  show that the face matches its own face and not matches other faces. 
     In the present invention, the stereotactic matching can be implemented by exploring multiple frames through temporal information as shown in  FIG. 18 . 
     Since many sensors can capture a video of a person, each of these frames can be transformed into the UCS domain of that person. The data can be averaged to obtain a better estimate of the coordinates of the landmarks in the UCS domain. This results in a better matching score of that person since the output error of a correct match goes down by 1/M where M is the number of temporal images used. Meanwhile, for an incorrect match, the averaging would not decrease the error since there is still a mismatch. 
       FIG. 19  shows that when 10 frames are used, the matching error is reduced to 0.74 from the error range between 5 and 10. 
       FIG. 20  illustrates a matching result of matching one thermal probe face image to two visible gallery images. The matching scores are displayed for the probe face with each face in the gallery. The smallest score corresponds to the correct matching. 
       FIG. 22  depicts a high level block diagram of a general purpose computer suitable for use in performing the functions described herein, including the steps shown in the block diagrams, schematic representations, and/or flowcharts. As depicted in  FIG. 22 , the system  500  includes a processor element  502  (e.g., a CPU) for controlling the overall function of the system  500 . Processor  502  operates in accordance with stored computer program code, which is stored in memory  504 . Memory  504  represents any type of computer readable medium and may include, for example, RAM, ROM, optical disk, magnetic disk, or a combination of these media. The processor  502  executes the computer program code in memory  504  in order to control the functioning of the system  500 . Processor  502  is also connected to network interface  505 , which transmits and receives network data packets. Also included are various input/output devices  506  (e.g., storage devices, including but not limited to, a tape drive, a floppy drive, a hard disk drive or compact disk drive, a receiver, a transmitter, a speaker, a display, a speech synthesizer, an output port, and a user input device (such as a keyboard, a keypad, a mouse and the like). 
     Although various preferred embodiments of the present invention have been described herein in detail to provide for complete and clear disclosure, it will be appreciated by those skilled in the art, that variations may be made thereto without departing from the spirit of the invention. 
     As used herein the terminology yaw means movement of a subject to the left or right. 
     As used herein the terminology “pitch” means the tilting of the subject&#39;s head on the y axis. 
     As used herein the terminology “roll” means the in plane rotation on the z-axis. 
     As used herein, virtual neutral pose position or virtual vertical plane position is defined as when the head is facing the camera or image producer and tilting such that the midpoints of the eyes and corners of the mouth are on a vertical plane and the corners of the eyes and corners of the mouth are equally distant from a face center line (as shown in  FIG. 11 ). 
     It should be emphasized that the above-described embodiments are merely possible examples of implementations. Many variations and modifications may be made to the above-described embodiments. All such modifications and variations are intended to be included herein within the scope of the disclosure and protected by the following claims. The invention has been described in detail with particular reference to certain preferred embodiments thereof, but it will be understood that variations and modifications can be effected within the spirit and scope of the invention.