# Patent

## Title
Distance iris recognition

## Abstract

A system for one dimensional segmentation of an iris of an eye into a map of the iris and classification of the map into unaffected areas and affected areas. Also, the system may provide for regular shape fitting of the areas for normalization and identifying the unaffected areas as symmetric segments. Further, the system may assign weights to the unaffected areas and the affected areas of the map of the iris and an enrolled map of an iris and their corresponding bins for matching purposes.

## Background
This application is a Divisional Application of U.S. application Ser. No. 11/275,703, filed Jan. 25, 2006, which claims the benefit of U.S. Provisional Application No. 60/647,270, filed Jan. 26, 2005. This application is a Divisional of U.S. application Ser. No. 11/275,703, filed Jan. 25, 2006, which in turn is a continuation-in-part of U.S. application Ser. No. 11/043,366, filed Jan. 26, 2005.
U.S. application Ser. No. 11/275,703, filed Jan. 25, 2006, is hereby incorporated by reference. U.S. Provisional Application No. 60/647,270, filed Jan. 26, 2005, is hereby incorporated by reference. U.S. application Ser. No. 11/043,366, filed Jan. 26, 2005, is hereby incorporated by reference.


BACKGROUND
The present invention pertains to recognition systems and particularly to biometric recognition systems. More particularly, the invention pertains to iris recognition systems.
SUMMARY
The invention is a system that incorporates certain improvements which support biometrics technology for person recognition from afar.



BRIEF DESCRIPTION OF THE DRAWING
FIG. 1 shows a basis for iris center localization;
FIG. 2 shows iris references for image polar mapping;
FIG. 3 shows one dimensional polar segmentation;
FIG. 4 shows image segmentation omitting normalization for illustrative purposes;
FIG. 5 shows normalization and treating eyelid occlusions;
FIG. 6 shows segmentation and normalization conducted simultaneously;
FIG. 7a shows examples of side effects of closures, eyelashes, and slanted oriented irises;
FIG. 7b shows examples of the results of an enhanced feature extraction;
FIG. 8 shows a perspective plane orientation in image acquisitions;
FIG. 9 shows an outline of segmentation of the iris into various regions and classifying the regions as unaffected and affected;
FIG. 10 shows the symmetric portion of the iris representing the non-occluded areas where the borders of the iris and sclera are visible;
FIG. 11 shows a determination of the center of the ellipse with an intersection of circles and an elliptic fitting scheme;
FIG. 12 relates to an adjustment/calibration of an iris radial longitude;
FIG. 13 is associated with a program for regenerating circles form ellipses;
FIG. 14 may show a visualization of the mixture modeling approach as applied to iris segmentation;
FIGS. 15a, 15b and 15c show a mixture modeling based analysis;
FIGS. 16a, 16b and 16c show another mixture modeling based analysis; and
FIG. 17 illustrates an example computer system usable in conjunction with certain illustrative instances of the present system.



DESCRIPTION
The invention may provide people identification and verification, using an “iris-at-a-distanceTM” or “Iris@aDistanceTM” (Honeywell International Inc.) system. Iris recognition has been recently recognized and has gained a lot of attention due to its high reliability in identifying humans. Its suitability as an exceptionally accurate biometric derives from its extremely data-rich physical structure, genetic independence (no two eyes are the same even for twins), stability over time, and non-contact means (a feature important for non-cooperative subjects). The present invention may help provide reliable calibration and an efficient segmentation (i.e., localization) of the iris-at-a-distance, resulting in better extraction of the iris features that are eventually converted into a numeric code. The iris codes may be compared with previously generated iris codes for verification and/or identification purposes.
As stated above, iris features are a reliable/accurate biometric measure. However, from a technical point of view, this accuracy lies heavily on how best the iris and pupil are segmented and extracted. When analyzing iris images of cooperative subjects, the segmentation approach is a relatively straightforward process of edge detection and circular fitting. However, this is often not the case for iris-at-a-distance systems, which often do not enjoy the cooperation of the subject. For most cases in an iris-at-a-distance system, only a portion of the iris is captured due to, for example, closure effect and/or eyelash and eyelid occlusions. Furthermore, given that the subject is typically not asked to cooperate, a tilted head or a rotated iris typically must also be considered. The present invention addresses these challenges, and in some cases, extracts accurate segments of the iris borders, among other things, in an iris-at-a-distance environment.
In one illustrative example, the process may include: 1) using a POSETM (i.e., polar segmentation) technique to move virtually immediately the analysis to a polar domain and execute a 1-D segmentation of the iris borders; 2) using one or more symmetry properties to detect one or more non-occluded areas of the iris—non-symmetric regions can correspond to areas partially covered by eyelashes, eyelids, and so forth (thus asymmetric)—and, in some cases, can limit the analysis to only those reliable segments where the iris and the sclera are detected (as noted in this description relative to symmetry); 3) once orientation is detected, identifying the nominal angles with the least likelihood of distortions (i.e., occluded or deformed due to orientation) by, for example, estimating the ellipse parameters from nominal angles, and computing a calibration factor; 4) using a rotated ellipse detection technique that uses overlapping variable circles to detect the iris borders modeled as elliptic or irregular shapes rather than circles, and/or using a least square fitting to estimate the elliptic parameters and orientation; and 5) using mixture modeling to handle variation in the iris textures. The sclera is the tough white fibrous outer envelope of tissue covering the entire eyeball except the cornea.
The present invention is well suited for high-security access control or “at-a-distance biometrics” applications where less control is exercised on subject positioning and/or orientations. Such operations may include, for example, subjects captured at variant ranges from the acquisition device, and/or may not have the subjects eye(s) directly aligned with the imaging equipment. Usually, for such applications, it is difficult to implement the level of control required by most of the existing art to enable reliable iris recognition. The present invention may help cope with asymmetry in acquired iris images, and may further help under uncontrolled environments as long as some of the iris annular is visible.
A leading concern in the development of biometric systems is how to avoid rejecting valid subjects or validating imposters. The iris recognition may provide a reliable solution by offering a much more discriminating biometric than other types of biometrics including face and fingerprint recognition techniques.
The following desirable properties make iris recognition technology as a potential reliable personal identification tool. For instance, this technology may have the following noted characteristics. First, it may have uniqueness and genetic independence. Eyes have an extraordinary structure that is unique to each human being. Unlike other well known biometric technologies, such as face-prints, fingerprints, irises are unique to a person and even among genetically twin individuals. Although the striking visual similarity of identical twins reveals the genetic penetrance of facial appearance, a comparison of genetically identical irises reveals just the opposite for iris patterns. The biomedical literature suggests that iris features are as distinct for each human as fingerprints or patterns of retinal blood vessels. Second, it has an extremely data-rich physical structure. There may be sufficient texture to make adequate discrimination between subjects. Third, there is no aging effect, that is, there is stability over life of the iris features. Fourth, it has non-invasiveness. The iris is considered an internal and unique organ, yet is externally visible and can be measured at a distance.
All these properties make iris recognition suitable for highly reliable personal identification then other notable biometrics such as facial recognition. It has been demonstrated that, unlike facial recognition, the phase structure and local features extracted from irises is purely epigenetic, so performance of iris recognition is not limited with application to identical twins or by the existence of partial genetic relationships.
Conversion of an iris annular image into a numeric code that can be easily manipulated may be essential to iris recognition. Computing iris features may use a good-quality segmentation process that focuses on the subject's iris and properly extracts its borders.
In support of an iris-at-a-distance application, one may use POSETM (Honeywell International Inc.) and a 1-D Polar based “segmentation approach”. POSE differs from the usual state-of-the art techniques in that it may conduct a one-dimensional segmentation process in the polar domain, replace the exhaustive search for geometric models (such as circles) and avoid the use of costly edge detections and curve fitting by simply executing a straightforward peak search on 1D signatures. The present approach may map immediately into the polar domain right from the start. To retrieve the iris annular region from the image, POSE may map the analysis at an earlier stage then previously done into the polar domain. By conducting the segmentation in polar domain, this may lead to a more efficient and faster process to execute not only the segmentation, but also calibration, and noise removal, all in one single step to generate a feature map for the encoding step.
In this framework, one may introduce a feature extraction technique in support of POSE to reduce the side effects of eye occlusions, eye closures, and deformation in captured iris images due to tilted heads and orientation perspectives projected in the image plane. The technique may be suited for the iris-at-a-distance applications, i.e., in cases where subjects are unaware that they are under surveillance, or in a crowded area, or even in cases where subject is aware of iris control but are non-cooperative. Such operations may include subjects captured at variant ranges from the acquisition device or may not have their eye directly aligned with the imaging equipment.
Here may be a designed and implemented an iris-at-a-distance biometric system specifically for these kinds of applications, which could function as an extremely reliable tool for people recognition and identification. Further, the present system may solve the asymmetry problem associated with image acquisition without the collaboration of the subjects and that it can operate under uncontrolled operations as long as some of the iris annular is visible.
Several issues may be researched further when developing an iris recognition solution. Relative to calibration issues, the iris inner and outer boundaries of iris may be approximated by ellipses than circles of irregular shapes using snake delineation. However, the two ellipses are usually not concentric. One may characterize the shape and texture of the structure of the iris having a large number of interlacing blocks such as freckles, coronas, furrows, crypts, and stripes.
A change in the camera-to-eye distance may result in variations in the size of the same iris. Preprocessing to reduce the side effects of non-uniform lighting or illumination reflections may be needed before one executes feature extractions procedures. Specular (mirror-like) reflections and the treatment of soft specular reflected dots may affect the segmentation analysis.
The outer boundaries of the iris may be captured with irregular edges due to presence of eyelids and eyelashes. Taken in tandem, these observations suggest that iris localization may be sensitive to a wide range of edge contrasts.
The orientation of head and eyes may result into different perspective of views of the iris circular shape. The captured shapes of the iris are usually far from circles or ellipses due to the orientation, tilt and slant angles.
Fundamental approaches to iris imaging matching may be built around Iridian (i.e., pertaining to an iris) solutions which inherit Daugman's algorithms for iris recognition. The Iridian technology algorithms may limit the extensibility of the iris recognition into real-time non-controlled environment. While certain operations are consistent with the preconditions, it is difficult to implement these existing technologies without the level of control required by the algorithms.
To extend to at-a-distance applications, some of the algorithm procedures that may require exhaustive processes which include Hough transforms. Some issues with a Hough method may include requiring threshold values to be chosen for edge detection, which may result into critical information (e.g., edge points) being removed/missed, and thus resulting into failure to detect the iris or pupil regions. The Haugh transform is computationally intensive due to its brute-force approach, and then may not be suitable for real-time applications. Furthermore, the method may fail where the image is subject to local noise in the eye image since it works on local spatial features.
A Daugman integro-differential operator may be seen as a variation of Haugh transform, since it appears to make use of derivatives of the image and perform a search to find geometric information that determines spatial parameters identifying the circles of the iris and pupil. The advantage the Daugman operator has over Haugh may be that it does not require threshold values as it is based on raw derivatives. However, it may fail when the image suffers from local spatial noise (e.g., specular reflections along the eye image, speckles due to digitization, and so forth).
As to eyelid occlusion and masking, one approach in dealing with the eyelid occlusions masks portions of the image may use linear fitting. However, the eyelid boundaries may be irregular due to the presence of eyelashes as well. Another approach in dealing with variable occlusions may be modeling the eyelids with parabolic curvatures and using the extracted configuration of model components to fine tune the image intensity derivative information. The alternative may cost a lot of computation given that is based on edge detection and non-linear curve fitting.
Another iris recognition technique may be based on a 1D process; but it is mostly tailored to improve the encoding scheme for better representation of the iris features rather than simplifying the segmentation process. Unlike the present approach which formulates the segmentation problem into a 1D framework, some methods construct a set of 1D intensity signals decomposed from the 2D constructed iris map. Gaussian moments applied on 1D representation of feature vectors may have been advocated by some as the best representation of local features that indirectly quantify the variations in textures due to coronas, stripes, furrows, and so forth. Such technique may still be based on 2D segmentation process to construct a 2D normalized polar map of the iris. It then may decompose the 2D map into a set of 1D intensity signals by simply averaging some segments of the rows of the map. Although, this effort may complement the present finding by executing the entire iris segmentation and encoding in a 1D domain. A 2D encoding scheme may appear to be much more reliable than a 1D process because of the inherited two dimensional relational spatial features in a 2D encoded signature.
A new 1D encoding scheme that may generate a 1D iris signature includes translation, rotation, illumination, and scale invariant. The merit of this method is that may allow users to enroll at lower level of iris image quality. This implies that the technique may not be as reliable as the 2D encoding scheme. However, that technique may support the concept of having a search mode before passing the limited potential subjects to a 2D encoding for final identifications. Once again, the present approach may implement a complete 1D segmentation and encoding iris technique.
Iris segmentation may be a factor to note. The first stage of iris recognition may be to isolate the actual iris region in a digital image of the eye. Often, for others, the iris region may be approximated by geometric models, i.e., two circles, to simplify the image processing segmentation of the iris and pupil.
Many others may make use of first derivatives of image intensity to signal the location of edges to segment the borders of these geometric models representing the boundaries of the iris. Unfortunately, most of existing algorithms may yield costly procedures for geometric model searches throughout the digital image.
Normalization may be a crucial step of the overall analysis given that that the eye image goes through various deformations. The intent may be to scale the extracted iris region cropped from the image to allow for a fair comparison with the database templates.
The dimensional inconsistencies among the captured iris images may be primarily due to many reasons such as a stretching of the iris caused by the pupil dilation from varying levels of illumination and lighting.
As to image range, the varying distance of image capture and imager orientation may be due to camera or head tilting and slanting. There may be local rotation of the eye within the eye socket. The subject or the subject's face might not be directly aligned with the acquisition device.
Image enhancement may be applied to minimize illumination artifacts, i.e., non-uniform brightness that illuminates areas more than others within the iris annular region. Reducing the illumination artifacts may improve subsequent encoding and feature extraction steps. The perspective orientation may be addressed before conducting feature extraction; however, this could add more computational burden on the system. The present segmentation algorithm does not appear to require these preprocessing steps to extract accurate features of the iris.
Encoding may be noted. In order to provide accurate recognition or identification of individual irises, one may need to extract the most discriminating information present in the polar presentation of the extracted iris. Just the significant features of the iris patterns may need to be encoded so that comparisons between two subjects can be made easy. The encoding scheme may be to generate a simpler template of a few bits that captures the essence of iris patterns. The extracted numeric code may then be used to compare it to multiple stored codes. Encoding the iris signature may include applying an encoding algorithm such as wavelet or Gabor filters or other techniques as listed below to extract textural information from images, such as the detailed patterns of the iris to produce a bitwise template containing a number of bits of information and excluding some of the corrupt areas using masking within the iris pattern. The choice of the encoding filters may be made on the basis of achieving the best recognition rate and preserving the iris patterns in the limited generated bitwise template.
In one's encoding process, regardless of the kind of filter chosen, one may generate more quantized phasor information represented by more than just two bits and prioritized with the most significant bits over the least significant bits when conducting matching. The merit of this scheme is that it may provide a quick way of matching subjects and also provide way to generate the most probable match instead of the best match when facing with poor quality iris images and iris patterns. For instance, one may conduct average weighting on the numeric code when conducting matching using any of the information divergence measure.
Encoding may include the actual encoding of the extracted features having different means of filtering and processing. The encoding mechanism may involve applying one or more selected filters to the segmented iris image(s). Some of the filters used by the state-of-the art techniques may include but are not limited to the following, such as wavelet/bank filters which may be also known as part of a multi-resolution technique. The wavelet approach may have an advantage over traditional Fourier transform in that the frequency data is localized. Gabor filters may also be capable of presenting a conjoint representation of the iris pattern in a spacial and frequency domain. Log Gabor filtering may be more reliable than Gabor filtering. Haar filters have been shown to possibly outperform Gabor filters. Laplacian filters may involve pyramid based decomposition to obtain a simplified version of the signal.
There may be matching which can be conducted using information theory measures, with a metric of the hamming distance. This process may execute matching between a query and encoded signatures. Several types of information measures, including but not limited to, may include a hamming code, a Euclidian code, a Jeffery code, a Kullback code, or any other standard information divergence measure which can be used to provide the weighted distance. The average weighted measure may be emphasized in the present approach. More weight may be distributed on the most significant bits versus lower bits. As a result, a phasor value may be represented by 8 or 16 segments (2^N) rather than just 2 bits in other codes. The weights may be distributed based upon the significance of the bit position.
POSE may perform iris recognition under suboptimal image acquisition conditions. The technique may be used for iris segmentation to detect all boundaries (inner, outer, eyelid and sclera and horizon) of the image iris. This technique may be well suited for high-security access control or iris-at-a-distance applications with less control exercised on subject positioning or orientations. Such operations may include subjects captured at variant ranges from the acquisition device or may not have their eye directly aligned with the imaging equipment. Usually, for such applications, it may be difficult to implement the level of controls required by much of the existing art to enable reliable iris recognition operations. The present approach of iris recognition may cope with asymmetry in acquired iris imaging and it may operate under any uncontrolled operations as long as some of the iris annular is visible.
There may be center pre-estimation. In the following analysis one does not approximate the iris boundaries as circles or ellipses which may take any other irregular shapes. The analysis still holds for any other deviated shapes due to distortions. One may localize the iris using a present simpler method and more general than done otherwise in the related art.
Assume that ρ(x,y) is a continuous density function on the blob region R as a function of the intensity function of the image. Then the center of mass may be defined by the following equation,




  
    
      C
      ⁡
      
        (
        
          
            
              x
              ^
            
            o
          
          ,
          
            
              y
              ^
            
            o
          
        
        )
      
    
    =
    
      (
      
        
          
            m
            x
          
          M
        
        ,
        
          
            m
            y
          
          M
        
      
      )
    
  
  ,




where M is the mass within the blob, and the moments of mass on both axes are,





  
    
      
        
          
            
              
                m
                x
              
              =
              
                ∫
                
                  
                    ∫
                    R
                  
                  ⁢
                  
                    
                      ρ
                      ⁡
                      
                        (
                        
                          x
                          ,
                          y
                        
                        )
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    y
                    ⁢
                    
                      d
                      x
                    
                    ⁢
                    
                      d
                      y
                    
                  
                
              
            
            ;
          
          ⁢
          
            
          
          &
        
        ⁢
        
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              m
              y
            
            =
            
              ∫
              
                
                  ∫
                  R
                
                ⁢
                
                  
                    ρ
                    ⁡
                    
                      (
                      
                        x
                        ,
                        y
                      
                      )
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  x
                  ⁢
                  
                    d
                    x
                  
                  ⁢
                  
                    d
                    y
                  
                
              
            
          
          ;
        
      
    
    
      
        (
        1
        )
      
    
  




Measurement of the center of the inner boundary may be of interest. The segmentation analysis does not necessarily rely on exact measurements of centers of both boundaries. Additional steps to extract the exact center of the inner iris ellipse may be noted.

Iris center localization may be shown in FIG. 1. The blob analysis likely will not necessarily lead to an accurate center of the pupil 11. Hence, further processing may be required to come up with a more accurate measure to locate the pupil center. There may be a center 14, C({circumflex over (x)}o,ŷo). Assuming that the pre-estimated center 12, C(xo,yo), is within the pupil 11 circle, then one may show easily that the actual center of the pupil is given by

r(θ)=Radius w.r.t. C({circumflex over (x)}o,ŷo)  (2)

Thus,

xo={circumflex over (x)}o+Δr(θo)cos(θo)

and

yo=ŷo+Δr(θo)sin(θo)

where,





  
    Δ
    ⁢
    
      
    
    ⁢
    
      r
      ⁡
      
        (
        
          θ
          o
        
        )
      
    
  
  =
  
    arg
    ⁢
    
      
    
    ⁢
    
      
        max
        θ
      
      ⁢
      
        (
        
          
            
              r
              ⁡
              
                (
                θ
                )
              
            
            -
            
              r
              ⁡
              
                (
                
                  θ
                  +
                  π
                
                )
              
            
          
          2
        
        )
      
    
  




In particular, if one has non-deformed iris 13 shape (i.e., virtually perfect circles), the equations above may be simplified to the following:




  
    
      x
      o
    
    =
    
      
        (
        
          
            
              r
              0
            
            ⁡
            
              (
              
                x
                ,
                0
              
              )
            
          
          +
          
            
              r
              π
            
            ⁡
            
              (
              
                x
                ,
                0
              
              )
            
          
        
        )
      
      2
    
  




  
    
      y
      o
    
    =
    
      
        (
        
          
            
              r
              
                
                  -
                  π
                
                /
                2
              
            
            ⁡
            
              (
              
                0
                ,
                y
              
              )
            
          
          +
          
            
              r
              
                π
                /
                2
              
            
            ⁡
            
              (
              
                0
                ,
                y
              
              )
            
          
        
        )
      
      2
    
  



The segmentation process may include the 1D POSE technique, although other approaches may be used. To simplify the segmentation process, one may propose mapping the analysis at an early stage to polar domain to execute the segmentation using 1D analysis procedure. In addition, the analysis may be conducted only on a subsection of the image surrounding the iris 13. In many other techniques, extraction of polar representation occurs near the end of the analysis.
A rapid polar conversion from an approximate center may permit a fast 1D segmentation in polar domain. The POlar SEgmentation (POSE) may yield rapid extraction of the apparent pupil and iris boundaries using one dimension signal processing. The analysis may detect all kind of boundaries; including non-elliptic boundaries (i.e., geometrically and biologically distorted images). The approach may handle line-of-sight boundaries at the far side of a significantly rotated head, eye, and/or iris.
Detection of a starting point within the pupil may be sufficient to initiate the mapping to a polar domain. Further adjustment of the pupil center may be considered as part of the POSE technique in the polar domain. The central point estimate may be based on any blob analysis, thresholding (presuming that the pupil resides at the darkest contrast of the eye image), or other approaches. Unlike known related art, the pupil center may be used as a basis of the analysis. The iris region may be normalized further to be centric/centered with respect to the same pupil center.
One may start the process by a rapid computation of a pupil-centered polar representation of the region-of-interest (ROI) surrounding the iris 13. The ROI may be mapped to the polar domain with respect to the estimated pupil center 12, C(xo,yo):




  
    
      I
      R
    
    ⁡
    
      (
      
        x
        ,
        y
      
      )
    
  
  ⁢
  
    
  
  ⁢
  
    ⟶
    
      C
      ⁡
      
        (
        
          
            x
            o
          
          ,
          
            y
            o
          
        
        )
      
    
  
  ⁢
  
    
      I
      R
    
    ⁡
    
      (
      
        θ
        ,
        r
      
      )
    
  




Based upon the predefined/estimated center, one may proceed by estimating an approximate iris 13 width. Then one may use the predefined center and iris width to execute the mapping immediately before any actual segmentation is executed.
FIG. 2 shows the referenced parameters that are used to map the image domain to the polar domain, that is, the Figure shows iris references for image polar mapping One may then construct a 1D function in the polar domain, which then may be convolved to locate the boundaries of the pupil 11 and the iris 13 of an eye 10. The analytical 1D function may be defined at each angle as a function of the radius variable, r:




  
    
      
        
          
            f
            θ
          
          ⁡
          
            (
            r
            )
          
        
        =
        
          ∫
          
            
              ∫
              
                
                  Δ
                  ⁢
                  r
                
                ,
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  θ
                
              
            
            ⁢
            
              
                
                  I
                  R
                
                ⁡
                
                  (
                  
                    θ
                    ,
                    r
                  
                  )
                
              
              ⁢
              
                
              
              ⁢
              
                d
                θ
              
              ⁢
              
                d
                r
              
            
          
        
      
    
    
      
        (
        
          3.
          ⁢
          
            
          
          ⁢
          a
        
        )
      
    
  




The estimated location of the iris 13 with respect to the center 12, C(xo,yo) may be computed as follows:

riris,rpupil=maxrk{ƒ0(r){circle around (×)}u(r)}  (3.b)

where u(r) may be a step function used to detect the contrast changes in the function f( ).

FIG. 3 shows 1D polar segmentation. The present approach may be further simplified to 1D match filter to detect the exact edges of boundaries of the iris 13. The segmentation may then be conducted at all desired predefined angles. The peaks 15 and 16 from left to right represent the radii of the pupil 11 and the iris 13, respectively.
FIG. 4 shows image segmentation omitting normalization for illustrative purposes (note that normalization is normally conducted on 1D signals extracted during the segmentation process). As shown in this Figure, eyelids 17 and 18, and/or eyelashes 19 may obscure some of the iris 13 annular region. Related art techniques may treat these obscurations as noise and tend to isolate the obscured region by first detecting the edges of the eyelids/eyelashes by fitting multiple linear/nonlinear curvatures/lines to the upper and lower eyelids using, for example, a Hough transformation, in many cases, thus adding more computational burden to already complicated procedures. Other alternative techniques make use of the Gabor filters and thresholding on a variance of intensities and connectivity criterion so that estimated curvatures lead to connected points of eyelash representation. Unlike these techniques, the present isolation of affected areas of the iris 13 may be done by comparing the expected radius segment 20 (median value of estimated segments 20, the expected length may be limited to only measurements taken at the nominal angles with least likelihood of distortions or noise) with the measured width. If there is a significant reduction in measured width, then the data may be treated as noisy and treated separately.
FIG. 5 shows normalization and treating eyelid 17, 18 occlusions 21. There may be two present alternatives that may handle the eyelid occlusions 21. First, one may mask all the corresponding pixels 22 at the affected angles. Although, this approach is simpler than the following approach, it tends to exclude some pixels that may end up being crucial for discriminations. Second, one may mask only the affected region 23 of the iris 13. The region captured in FIG. 4, under the eyelid-occluding the top portion 21, may be accounted for by normalizing the region based on the expected width value (i.e., median value of the measurements). Thus, just the pixels 23 that correspond to the iris 13 occluded underneath the eyelid 18 are masked out, as shown in FIG. 5.
Iris 13 normalization may be noted. For the purpose of achieving more accurate recognition results, it may be necessary to correct for the shape deformation and bring uniformity into the shape of the iris 13 before passing the polar data into the encoding scheme. Typical techniques conduct normalization towards the end. With the present technique, normalization may be conducted during the segmentation process.
It is expected that at least four major artifacts may result into deformation of the iris 13 circular shape, thus making a normalization scheme necessary while mapping the iris pixels into the rubber polar presentation. First, there may be range differences. Irises of different people may be captured at different ranges of the camera even under constrained conditions. Second, there may be iris dilation. It may be natural to expect that the size of iris change because of illumination changes. Usually, this type of change may be uniform across the annular iris shape. Third, there may be a natural iris deformation. There may be cases where the iris is asymmetric by nature. Pupil 11 region is not necessarily concentric within iris 13 region, and may be slightly nasal. Although, this may be used as an additional discriminating feature for individuals, one may bring all textures into a single referenced shape for simplicity. The present normalization scheme may be able to preserve the discriminating features of deformed iris 13 into the new presentation. Fourth, there may be perspective orientations which may include any line-of-sight boundaries at the far side of a significantly rotated head, eye 10, and/or iris 13.
Unlike the related art, the present normalization process may be conducted as the iris 13 segment is extracted from POSE technique. For instance, let




  
    
      
        
          
            S
            θ
          
          ⁡
          
            (
            r
            )
          
        
        =
        
          
            ∫
            Δθ
          
          ⁢
          
            
              
                I
                R
              
              ⁡
              
                (
                
                  θ
                  ,
                  r
                
                )
              
            
            ⁢
            
              
            
            ⁢
            
              d
              θ
            
          
        
      
    
    
      
        (
        
          4.
          ⁢
          a
        
        )
      
    
  




be the 1D representation of the polar map at angle Θ; for simplicity, one may assume that uniform integration is used along the Θ axis. Then, the normalized signal may be based on the interpolation/decimation of the signal. The normalized output is stated as follows:

Pθ(r)=L↑[h↓Sθ(r)]  (4.b)

L may be the desired dimension (i.e., number of rows) for the polar representation of the iris 13 texture. The variable h may vary based upon the deformation type and amount of degradations. FIG. 6 shows segmentation and normalization of the iris 13 conducted simultaneously. A feature extraction technique may be noted.

The present approach addresses a number of challenges and may makes headway towards the commercialization of a reliable iris system without any constrains on the subject of interest. Interest may be in particular interest in recognizing/identifying subjects in large open spaces like airport gates, plazas, potentially crossroads, and sensitive checkpoints. These applications may constitute iris-at-a-distance with no constrains imposed in image acquisitions, and which may be referred to hereon as “iris-at-a-distance” (IAD).
A comprehensive iris-at-a-distance system may depend primarily on two different technologies—optical design and a computer vision solution. The computer vision solution may be stressed in this section.
The present system may be the only one that features a reliable segmentation stage for irregular acquired iris 13 images. The practical aspects of fielding an iris-at-a-distance system may be noted. There are a number of requirements that a surveillance system may need to fulfill to function properly and be commercially viable. The occlusion issues of eyelashes 19 and eyelids 17, 18 may be some of the challenges. The POSE segmentation technique may be used to solve some of the challenges. POSE may have been proven to be a more reliable mechanism than such things as spatial thresholding or Hough space based techniques.
Spatial thresholding may do a questionable job of classifying all the relevant iris pixels of the eye image. Edge detection and Hough transform may require an intensive computational load to execute all necessary steps to extract the iris 13 segmentation. In addition, fitting the iris into predefined shapes such as circles, ellipses or predefined curves (representing eyelashes) may be closest to reality only under constrained conditions where the subject iris is placed through guidelines to picture a perfect perspective of the iris 13. This scenario appears to be far from real in the iris-at-a-distance approach.
For the iris-at-a-distance, it may be challenging to implement the level of controls required by much of the related art to enable reliable iris recognition operations. Three major challenges arise when dealing with such application to uncontrolled operations are closure of the captured eye 10, eyelash 19 and eyelid 17, 18 occlusions 21, and artifacts due to subjects not having their eye directly aligned with the imaging equipment. As shown in FIG. 7, the effects of eyelashes 19 or eye-closure or any line-of-sight boundaries at the far side of significantly rotated head, eye 10, or iris 13 may result in some deformation due to orientation perspectives.
FIG. 7a shows an example of side effects of closures, eyelashes 19, and slanted oriented irises 13. Many of the related art algorithms fail when faced with closure, occlusions, and deformation in image acquisition.
The present system may use a feature extraction technique based on the POSE segmentation approach to handle such scenarios. The present approach may cope well with asymmetry in acquired images and it may operate under uncontrolled operations. FIG. 7b shows the outcome of the present approach with enhanced feature extraction 24 results.
Perspective plane transformation may be noted. A present preprocessing approach may be based upon perspective plane orientation that addresses some of the deformations in image acquisitions. Such operations may include subjects captured at variant angles from the acquisition device or the subjects may have their eye 10 or iris 13 not directly looking into the imaging equipment.
The present preprocessing approach for estimating the orientation in space of the face surface from an imaging angle may be based upon the fact that some of these orientation angles are provided by the face recognition tool. One may derive an accurate model of the captured image and its respective projection in the imaging plane. The perspective projection may have a dominant and fundamental role in the present preprocessing operations to detect iris 13 positions with respect to the face, and not with respect to the imaging plane. Segmentation of the iris 13 may be considerably simplified if the effects of the perspective projections are eliminated. The POSE technique may be easily applied once on the recovered face image with aligned perspective. The orientation estimation may be essential if the camera is situated at an angle of the normal axis of the eye gaze. For instance, the segmentation and feature extraction 24 procedures may be considerably simplified if the effects of the perspective projection are eliminated first, thus reducing the asymmetry in the captured frames and producing accurate segmentations.
FIG. 8 shows a perspective plane orientation in image acquisitions. To solve the perspective transformation, the analysis may involve a high-order partial-differential equation. In the following, one may assume a pinhole perspective projection model to provide an estimate of the geometric view of the actual face from the camera 25 perspective. One may assign a camera plane or imaging plane 27 coordinate system, {right arrow over (xc)}=(xc,yc) (imaging plane coordinates) and the actual face plane 26 coordinate {right arrow over (xo)}=(xo,yo). One may use the elevation angle of the face normal for representing the orientation of the face. The elevation angle, α, may be the angle between the face normal and the camera system axis 28. One may use the tilt angle, β, to denote the angle between x-axis and the projection of the face normal onto the imaging plane 27. One may denote θ to be the angle between the upper eyelid 18 center and the back projection of the image tilt vector (cos β, sin β) on the face plane 26. The tilt vector may be associated with the elevation angle and indicate how much the surface of the face is tilted from the perspective of the imaging plane 27. Thus, the coordinate transformation from the actual face to the imaging plane 27 may be given by




  
    
      
        
          
            
              x
              →
            
            c
          
          =
          
                        ⁢
            
              T
              α
            
            ⁢
                        ⁢
            
              
                x
                →
              
              o
            
          
        
        ⁢
        
          

        
        ⁢
        where
        ⁢
        
          

        
        ⁢
        
          
                        =
            
              [
              
                
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      β
                    
                  
                  
                    
                      
                        -
                        sin
                      
                      ⁢
                      
                        
                      
                      ⁢
                      β
                    
                  
                
                
                  
                    
                      sin
                      ⁢
                      
                        
                      
                      ⁢
                      β
                    
                  
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      β
                    
                  
                
              
              ]
            
          
          ;
        
        ⁢
        
          

        
        ⁢
        
          
            
              T
              α
            
            =
            
              [
              
                
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      α
                    
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    1
                  
                
              
              ]
            
          
          ;
        
        ⁢
        
          

        
        ⁢
        and
        ⁢
        
          

        
        ⁢
        
                    =
          
            
              [
              
                
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      θ
                    
                  
                  
                    
                      
                        -
                        sin
                      
                      ⁢
                      
                        
                      
                      ⁢
                      θ
                    
                  
                
                
                  
                    
                      sin
                      ⁢
                      
                        
                      
                      ⁢
                      θ
                    
                  
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      θ
                    
                  
                
              
              ]
            
            .
          
        
      
    
    
      
        (
        1
        )
      
    
  




To extract the actual coordinates of a point P(x,y) captured in the imaging plane 27, one may express the original coordinates as an inverse transformation of a point on the image plane to a point on the face surface plane 26 due to the perspective back projection.

One may use an enhanced feature extraction 24 approach. The purpose of this effort may be to explore that which will enhance current capabilities of the POSE technique in iris 13 recognition as applied to more unconstrained iris capture in more diverse environments where occlusions 21 and irregularities are expected. The present approach may operate under any of these uncontrolled operations as long as some of the iris 13 annular is visible in the acquired image. An overall view of the present approach is shown in FIG. 9.
The approach may be initiated by a step 31 segmenting the iris 13 region using the POSE technique; then one may do a step 32 to classify the iris regions into unaffected versus affected regions. In a step 33, one may process the unaffected regions to fit regular or irregular, i.e., elliptic, iris shapes. This may involve an elliptic fitting step 38 for normalization. The present curvature fitting approach may be general enough to cover virtually all possible cases for irregular shapes using snake delineation via POSE. In a step 34, one may then process the affected regions to cover any possible regions that can be calibrated using the parameters of the extracted shape, i.e., calibration step 37 to reduce occlusion effects. In a step 35, one may cluster out the areas where iris 13 is completely covered with eyelashes 18 or eyelids 17 and 18 using the mixture modeling technique step 36. There may be inputs from steps 36, 37 and 38 for an encoding step 39, where the map of the iris 13 may be converted into a numeric bitwise code.
The first stage of iris recognition is to isolate the actual iris region in a digital eye image. The POSE technique may be successfully applied to extract the iris region with least computation possible, i.e., a 1D based analysis. The success of the present POSE segmentation may depend on the imaging quality of eye images. For instance, one may presume that images are passed through preprocessing stages to eliminate artifacts due to specula reflections or any other types of background noise. In addition, POSE may localize the analysis around the pupil center and immediately map it to the polar domain, without a need to adjust for non-concentric to the iris center which is a process that may be required by the related art.
FIG. 9 shows a feature extraction 24 approach using POSE 31. Iris classification step 32 may use symmetry in POSE edges. POSE may or the related art may segment the iris 31 as long as it is visible. POSE may require only a portion not necessarily the entire iris be visible. For instance, the eyelids 17 and 18 and eyelashes 19 may normally occlude the upper and lower parts of the iris 13 region. The related art may make use of a Hough transform to solve the extraction of eyelashes 19 or eyelids 17 and 18 approximating the upper 18 and lower 17 eyelids with parabolic arcs. In performing the related-art edge detections, the derivative of horizontal direction may be adjusted for detecting the eyelids 17 and 18. This may work when the eye 10 or head is not tilted from the horizontal axis. In addition, there may be a number of problems with the Hough transform method. First, it may require thresholding to be chosen for edge detection, and this may result in critical edge points being removed, resulting in a failure to detect arcs (or even circles for the iris borders). Moreover, the Hough method may be computationally intensive due to its brute-force approach and thus may not be practical for the iris-at-a-distance application.
An approach may be introduced for detecting the orientation of the eye 10 without having to add extra processes. This approach may be based on detecting the symmetry of the POSE resulting edge points. As shown in FIG. 10, the symmetric portion 41 of the iris 13 may represent the non-occluded areas where the borders of the iris 13 and sclera are clearly visible. One may indicate that the asymmetric regions represent the eyelids 17 and 18 and eyelashes 19 (obviously the asymmetry is due to occlusion of the sclera). There may be an iris image classification using asymmetry in detected POSE edges. This technique may provide a way to identify best visible iris areas from affected areas. The affected area may be either masked out or subjected to further analysis to extract more texture patterns if deemed necessary.
The elliptic-fitting step 38 based normalization may be used, as noted in FIG. 9 and introduced in FIG. 11. Firstly, the edge map may be generated by radially POSE segmentation over all angles. The symmetric points 43 of the edge map may then be selected for the normalization process. From these symmetric arc points 43, votes may be cast in within POSE space (i.e., ranges are predefined by POSE edge map) for the parameters of variable circles 42 passing through each edge point 43 of the symmetric arcs. The ellipse 44 detection scheme may then be based on the overlapping of these variable circles 42. These edge points 43 may be used as the center coordinates of these variable circles 42 which are able to define any circle according to the equation,

r2=(x−xo)2+(y−yo)2 

The circle radius may vary based upon how distant the center from the extreme edge point as shown in FIG. 11. The corresponding radius variation may be computed as





  
    r
    ⁢
    
      
        d
        r
      
      
        d
        y
      
    
  
  =
  
    
      (
      
        y
        -
        
          y
          o
        
      
      )
    
    ⁢
    
      (
      
        1
        -
        
          
            (
            
              a
              /
              b
            
            )
          
          2
        
      
      )
    
  




and the resulting radius for each iteration may be defined as





  
    r
    k
  
  =
  
    
      
        a
        k
        2
      
      -
      
        
          (
          
            1
            -
            
              
                a
                k
                2
              
              /
              
                b
                2
              
            
          
          )
        
        ⁢
        
          
            (
            
              b
              -
              dy
            
            )
          
          2
        
      
    
  




A maximum point in the POSE space may correspond to the intersection of majority of constructed circles and thus to the corresponding radii and center 45 of the ellipse 44 as illustrated in FIG. 11.

Calibration may be used to reduce occlusion effects. An automatic POSE segmentation model may prove to be successful in detecting virtually all iris regions including area portions occluded by the eyelid 17 or 18. However to conduct an accurate normalization to bring the analysis into uniformity and make significantly accurate matching against templates in the database, one may need to resample the points detected along the radial axis based on the actual radius of the iris 13 and not on the detected one as it does not represent the complete radius. Thus, measurements may be adjusted and rescaled accordingly based upon the best elliptic fitting to the edge map points detected in the nominal areas (i.e., symmetric edges). The approach may be shown in the following steps.
First, one may estimate the ellipse parameters from nominal map edges with the following,




  
    
      
        (
        
          
            x
            -
            
              x
              o
            
          
          a
        
        )
      
      2
    
    +
    
      
        (
        
          
            y
            -
            
              y
              o
            
          
          b
        
        )
      
      2
    
  
  =
  1




Second, for all angles θ, one may compute calibration factor, as follows,




  
    
      λ
      =
      
        
          r
          ~
        
        r
      
    
    ;
    
      
    
    ⁢
    
      given
      ⁢
      
        
      
      ⁢
      that
    
  




  
    r
    =
    
      
        
          
            (
            
              
                a
                
                  -
                  2
                
              
              +
              
                (
                
                  
                    
                      tan
                      2
                    
                    ⁢
                    θ
                  
                  
                    b
                    2
                  
                
                )
              
            
            )
          
          
            -
            1
          
        
        ⁢
        
          (
          
            1
            +
            
              
                tan
                2
              
              ⁢
              θ
            
          
          )
        
      
    
  



Third, one may rescale the map based upon calibration factors.

One may have an adjustment/calibration of the iris radial longitude. FIG. 12 and its notation may be observed.

(y−yo)=(x−xo)tan θ

r2=(x−xo)2+(y−yo)2 

r2=(x−xo)2(1+tan2θ)

((x−xo)/a)2+((y−yo)/b)2=1

((x−xo)/a)+tan2 θ/b2(x−xo)=1

(1/a2+tan2 θ/b2)(x−xo)2=1

(x−xo)=(1/a2+tan2 θ/b2)−1

r2=(1/a2+tan2 θ/b2)−1(1+tan2 θ)

The present technique is well suited for iris-at-a-distance. Such operations may include subjects captured at variant ranges from the acquisition device or may not have their eye directly aligned with the imaging equipment. Usually for such applications, it is difficult to implement the level of controls required by most of the related art to enable reliable iris recognition operations. The present concept here may cope with asymmetry in acquired iris imaging and it may operate under much any uncontrolled operations as long as some of the iris annular is visible. A purpose of the present system is to have capabilities in iris recognition as applied to more unconstrained iris capture in rather diverse environments. The present analyses may detect many sorts of boundaries, including irregular and non-elliptic boundaries, i.e., it may detect geometrically and biologically distorted images. The present approach may handle most any line-of-sight boundaries at the far side of a significantly rotated head level iris. The pupil 11 region is not necessarily concentric with the iris 13 region and may usually be slightly nasal.
Reliable iris recognition technology may be a reality. However, computing iris features requires a good-quality segmentation process that focuses on the iris boundaries extraction.
In the related art, it may usually be difficult to implement the level of controls required in that art to enable reliable recognition operations. The major challenges when dealing with such application of uncontrolled operations may include closure, occlusion, and artifacts due to subjects not having their eye 10 aligned with the imaging equipment.
The present system may improve the capabilities of POSE as to application to quite unconstrained iris 13 capture in more diverse environments where occlusions and irregularities are expected. The preprocessing approach here may be based on perspective orientation that addresses these deformations in image acquisition. One may assume a perspective projection model to provide an estimate of the geometric view of the actual face from the camera perspective.
POSE may map the analysis at an early stage into the polar domain. By conducting segmentation in the polar domain, one may be led to a more efficient and faster process to execute not only segmentation plus calibration and noise removal in one single step to generate the feature amp for the encoding scheme.
A least squares solution may used for the iris fitting step. The general quadratic curve equation,

ax2+2bxy+cy2+2dx+2fy+g=0,

is an ellipse, upon defining



  
    
      Δ
      =
      
        
        
          
            
              a
            
            
              b
            
            
              c
            
          
          
            
              b
            
            
              c
            
            
              f
            
          
          
            
              d
            
            
              f
            
            
              g
            
          
        
        
      
    
    ,
    
      

    
    ⁢
    
      J
      =
      
        
        
          
            
              a
            
            
              b
            
          
          
            
              b
            
            
              c
            
          
        
        
      
    
    ,
    and
  




  
    l
    =
    
      a
      +
      
        c
        .
      
    
  



For simplicity, one may rescale the equation by dividing the constants by g, assuming g≠0, thus:




  
    
      V
      →
    
    =
    
      
        [
        
          
            a
            g
          
          ⁢
          
            
          
          ⁢
          
            b
            g
          
          ⁢
          
            
          
          ⁢
          
            c
            g
          
          ⁢
          
            
          
          ⁢
          
            d
            g
          
          ⁢
          
            
          
          ⁢
          
            f
            g
          
        
        ]
      
      T
    
  
  ,




let g=1, and
measurement k may be formulated as follows:

mk=[xk2 2xkyk yk2 2xk 2yk]{right arrow over (V)}

∀k=1→N measurements, for each measurement (xk, yk), one may define {right arrow over (m)}k=[xk2 2xk yk yk2 2xk 2yk].




  
    
      Let
      ⁢
      
        
      
      ⁢
      M
    
    =
    
      [
      
        
          
            
              m
              1
            
          
        
        
          
            
              m
              2
            
          
        
        
          
            ⋮
          
        
        
          
            
              m
              N
            
          
        
      
      ]
    
  
  ,
  
    N
    ×
    5
    ⁢
    
      
    
    ⁢
    
      dim
      .
    
  




One may solve for the ellipse 44 with
{right arrow over (V)}=(MTM)−1MT{right arrow over (1)}, where




  
    
      1
      →
    
    =
    
      [
      
        
          
            
              -
              1
            
          
        
        
          
            
              -
              1
            
          
        
        
          
            
              -
              1
            
          
        
        
          
            ⋮
          
        
        
          
            
              -
              1
            
          
        
      
      ]
    
  
  ,




having an N×1 dim.

Assuming Δ≠0, J>0, and Δ/I<0. In that case, the center 45 for the ellipse 44 is (xo, yo), where

xo=(cd−bf)/(b2−ac),

yo=(af−bd)/(b2−ac),

the angle of rotation is

Ø=1⁄2cot−1((c−a)/2b)

and the semi-axis lengths are

a′={[2(af2+cd2+gb2−2bdf−acg)]/(b2−ac)[(c−a)(1+(4b2/(a−c)2)1/2−(c+a)]}1/2,

and

b′={[2(af2+cd2+gb2−2bdf−acg)]/(b2−ac)[(a−c)(1+(4b2/(a−c)2)1/2−(c+a)]}1/2.

Below is sketch of a program which may be applicable and scaled with a ratio a/b for regenerating circles from ellipses, in conjunction with FIG. 13.












function ellipse = find_me_ellipse(edges, rmin, rmax, num)


%***************************************************************************


%ellipse = find_me_ellipse(edges, rmin, rmax) -- returns the coordinates








%
of ellipse in an image using a simplified version of Hough



transform to detect overlapping circles at the center.


%
The image may be rescaled to transfer elliptic shapes into



circular shapes so that existing technology of circular



fitting can be applied to locate elliptic center and



eventually elliptic parameters as illustrated herein.


%







%INPUTS:









%
edges
- the image in which to find ellipses


%
rmin
- lower radius to search for


%
rmax
- upper radius to search for


%
num
- number of points to search for the ellipse shape







%OUTPUTS:









%
ellipse.x,.y.r
- ellipse coordinates







****************************************************************************


[rows, cols] = size(edges);


if ~exist(‘rmin’),  rmin = 5; end;


if ~exist(‘rmax’),  rmax = round((min(cols,rows))/2); end;








amax
= rmax(1);


ellipse.x
= 1;


ellipse.y
= 1;


ellipse.a
= 0;


ellipse.b
= 0;


ratioMax
= amax/bmax;


ratioMin
= amin/bmin;


ratio
= ratioMin:(ratioMax−ratioMin)/(num−1):ratioMax;


nradii
= (bmax−bmin+1);







% perform the circular Hough transform


hough = hough_circle(edges, bmin, bmax, ratio);


maxtotal = 0;


% find the maximum in the Hough space,


% and hence the parameters of the circle


for k=1:num,









scale = hough(:,:,n, k);



for n+1:nradii









layer = hough(:,:,n, k);



maxlayer = max2(layer);



if maxlayer > maxtotal



maxtotal = maxlayer



radius  = bmin+n−1;



[row,col] = find(layer == maxlayer);



ellipse.y = row(1);



ellipse.x = round(col(1)*scale);



ellipse.b = radius;



ellipse.a = round(radius*scale);



end % if;









end; % for n







end; % for scale.


%===========================================================================


function hough = hough_circle(edges, rmin, rmax, ratio)


%===========================================================================


% h = hough_circle(edge_img, rmin, rmax)








%
- takes an edge map image, and performs the Hough transform


%


% INPUTS:









%
edges:
the edge map image


%
rmin, rmax
limit search to this minmax range.


%







% OUTPUTS:









%
hough:
Hough transform







****************************************************************************








[rows,cols]
= size(edges);







[rowList, colsList] = find(edges ~=0);








len
= length(rowsList);


nradii
= (rmax−rmin+1);


num
= length(ratio);


hough
= zeros(rows,cols,bradii, num);







% for each scale,


for k = 1:num









scale = ratio(k);



XList = round(colsList/scale);



%for each edge point, draw circles of different radii



for index=1:len,









circle.y = rowsList(index);  % def circle coordinates



circle.x = XList(index);



for n=1:nradii









circle.r = n+rmin−1;



hough(:, :, n, scale) = add_circle(hough(:, :, n, scale), circle);









end % for n









end % for index







end; %scale


%









Texture segmentation may be accomplished using a mixture modeling (MM) technique. Segmentation may have originated in the related art based on a mixture of normals representation at the pixel level. This method may feature a far better adaptability to scene understanding and may handle even bimodal backgrounds (i.e., means to discarding moving objects of non-interest in a dynamic scene). In an attempt to mediate the computational burden required by the just noted approach, an adaptive mixture modeling approach may be implemented at the block level instead of pixel level analysis and simplified divergence measures may be introduced between the normals of incoming pixels/blocks and existing model distributions. Both alternatives may feature a far better adaptability to changes in image content and handle complex scenes.
Here, one may introduce a modification of mixture modeling approach of these two above-mentioned alternatives to handle the variation in the iris textures to take advantage of the statistical MM capabilities. The present approach may use a multi-normal representation at the pixel level. However, in the present approach, one does not necessarily model the pixels/blocks over-time; instead, one may model a block of pixels spatially with initialization executed on a predefined region of interest part of the iris area.
One may use the present histogram based approach for statistical distribution initialization instead of the related art approaches involving approximation, or expectation-maximization (EM) algorithm. In contrast to both EM and the approximation, the present histogram algorithm may provide more reliable initial statistical support that facilitates fast convergence and stable performance of the distribution classification.
One may also introduce an apparently reliable divergence measure as a modification to the spectral information divergence (SID) measure and relative entropy (RE) measure. The present divergence measure may be a simplified version of the SID measure used to match criterions between normals of scanned pixels and existing model normals. The measure may be a far superior measure to the fixed value (2-standard deviations), and much more reliable than the predefined STD. The measure may be for dual use, as described below, for the segmentation of the iris texture as well as for measuring similarities between iris signatures.
In the iris texture analysis, for every scanned pixel, one may check whether there is a match or non-match using the new divergence information measure. When a match is found, the model update may be performed using the first-in-first-out (FIFO) method while updating the weights of the model mix. When a match is not found, the update may be performed in a way that guarantees the classification of the scanned distribution in the outer set, associated with the eyelashes 19 or eyelid 17, 18 regions (non-iris region).
The present mixture modeling based concept may allow one to identify the eyelashes 19 and eyelids 17 and 18 pixels in each iris image while updating the description of each class of the iris texture's mixture model. The matched pixels may then be assembled into a continuous texture using a generic connected component algorithm.
Initialization 51 may be performed. A goal of the initialization phase is to set initial statistical parameters for the pixel clusters corresponding the iris 13 and non-iris 50 regions with reasonable values that represent their classes' distributions. These initial values may be used as a starting point for the dynamic changes in the pixel values across the iris 13 region and non-iris 50 regions. One may extract a predefined region 52 of the iris 13 based upon the symmetry based classifier with a sufficient number of pixels and then process them before expanding to the rest of the image. Each pixel 54 of the selected sub-region X may be considered as a mixture of three spatial-varying normal distributions:




  X
  ⁢
  
    ∼
  
  ⁢
  
    
      ∑
      i
    
    ⁢
    
      
        ω
        i
      
      ⁢
      
        N
        (
        
          
            μ
            i
          
          ,
          
            ∑
            i
          
        
        )
      
    
  




where 0≦ωi≦1;





  
    
      ∑
      i
    
    ⁢
    
      ω
      i
    
  
  =
  1




are the mixing proportions weights of the X variations within the iris region perspective. The notation N(μi,Σi) may represent a normal distribution with mean μ and covariance matrix Σ. For simplicity, one may consider in the present analysis that neighboring pixels are not correlated with each other and with identical standard deviation, thus the covariance matrix may be simplified to a scalar multiple of the identity matrix, i.e., N(μi,σi).

The initialization 51 in FIG. 14 may reveal iris distributions N({right arrow over (μ)}1,σ12), N({right arrow over (μ)}2,σ22), and N({right arrow over (μ)}3,σ32) as designated by 61, 62 and 63, respectively. Incoming evidence 67 from scanned pixel 54 may contribute ω1(I), ω2(I) and ω3(I) to distributions 61, 62 and 63, respectively. There may non-iris distributions N({right arrow over (μ)}4,σ42), N({right arrow over (μ)}5,σ52) and N({right arrow over (μ)}6,τ62), as designated by 64, 65 and 66, respectively. Incoming evidence 68 from scanned pixel 54 may contribute ν1(t), ν2(t) and ν3(t) to distributions 64, 65 and 66, respectively.
Related mixture modeling methods may initialize the pixel values either with random numbers, an approximation algorithm, or EM algorithm. The related art initialization methods appeared to be applied to time-varying distribution and even then did not appear to prove practical in most scenarios as they all appear to result in slow learning (approximation) and crude estimates (EM algorithm—a very computationally intensive technique) which usually may result into very biased results unless a large number of pixels is considered to bring the initialization phase under control. In contrast to the related art, and to avoid computational burden, one may presently make use of a very simple yet powerful method of clustering the accumulated pixels based upon their histogram and approximate their distribution based upon the variation in the histogram contrast. One may obtain the estimates of the three distributions by clustering the collected pixels into three classes that center within each class with largest dissimilarity among the three clusters. The weights may be linearly scaled based upon the fraction of pixels of each individual class with respect the overall accumulated pixels. The result may be a mixture model for three normal distributions per iris pixel with initial weights. These normal distributions may represent three potential different states for each iris pixel. One may intuitively claim that at least 3 classes are expected to categorize the depicted gray-scale iris pixels due to eyelash covering some of the iris regions as well as the texture within the iris may be categorized into at least 3 levels or classes of contrast. In experiments, the present histogram based clustering technique has appeared to prove a superior initialization method that caters to fast learning and exceptional stability to the subsequent main stage of iris segmentation. This appears to be especially true when initialization 51 happens in a not well predefined iris initial region.
It may be important to note that in both during initialization 51 and subsequent processing, one may chose to represent each iris pixel component of a pixel as a mixture of a predefined number (e.g., three may be used throughout the experiments) of normal distributions. One may arrive at this choice (i.e., K=3) after experimenting with clustering the resulted pixel contrast using the present optical and acquisition system. However, other system arrangements may require larger number of classes to cover the dynamic changes within the iris texture contrast variations. From various experiments, it appears that three normals per iris may be a sufficiently rich representation scheme to capture natural texture variation. Adding more normals appears to simply increase the computational load without improving the quality of segmentation.
One may have mixture model based iris segmentation. The initial mixture model may be updated dynamically thereafter while scanning the entire image. The update mechanism may be based on new evidence (new scanned pixels). None of the existing iris distributions may be dropped or replaced. However, distributions' parameters may change as it is dynamically updated. As for the non-iris distributions 64, 65 and 66, one of the existing distributions may be dropped and replaced with a new distribution if there is no match. While scanning throughout the image, at every point in space the distribution with the strongest evidence may be considered to represent the pixel's most probable iris state.
FIG. 14 may present a visualization of the mixture of normals model. The Figure may also reveal the mixture modeling approach as applied to iris segmentation. The present procedure may be used to update the mixing proportions weights of the mixture models with no thresholds or learning parameters. When a match is found, the model update may be performed and weights may be updated using adaptive weighting. When a match is not found, the update may be performed in a way that guarantees the inclusion of the incoming distributions in the non-iris distribution set. Once the minimum number of occurrence is reached, the least weighted distribution of non-iris may be replaced with the new most frequent distribution.
The algorithm may work as follows. First, the existing distributions of the predefined region of the iris may be initialized and weights may be determined based upon the fraction of the population set. Second, the algorithm may select the first Ko distributions that account for eyelash 19, eyelids 17 and 18, limbic and other possible items. One may typically set Ko=3 associated with three possible non-iris regions (i.e., eyelash 19, eyelid 17, 18, limbic and so forth). Third, the algorithm may then check for any scanned pixel value to be ascribed to any of the existing normal distributions. The matching criterion that one may use is an information divergence measure which may be a key contribution of the present approach from other similar methods. A divergence measure may be described herein.
Fourth, the algorithm may update the mixture of distributions and their parameters. The nature of the update may depend on the outcome of the matching operation. If a match is found within the iris or non-iris distributions, then the matched distribution may be updated using the method of moments. If a match is not found and the minimum number is reached, then the weakest distribution of the non-iris may be replaced with a new distribution. If the minimum required number is not reached, a counter may be incremented without replacing any distributions. The update performed in this case may preserve the nature of the iris distribution mixture while guaranteeing the inclusion of the new distribution in the non-iris set, which may be an aspect of the present method. The matching and model updates operations may be quite involved and are described in detail in the following. There may be the information divergence measure. The SID information measure between two distributions

ƒi ̃N[{right arrow over (μ)}i,σi2]

and

g ̃N[{right arrow over (μ)}g,σg2]

may be defined as:





  
    J
    ⁡
    
      (
      
        
          f
          i
        
        ,
        g
      
      )
    
  
  =
  
    ∫
    
      
        [
        
          
            
              f
              i
            
            ⁡
            
              (
              x
              )
            
          
          -
          
            g
            ⁡
            
              (
              x
              )
            
          
        
        ]
      
      ⁢
      
        log
        ⁡
        
          (
          
            
              
                f
                i
              
              ⁡
              
                (
                x
                )
              
            
            
              g
              ⁡
              
                (
                x
                )
              
            
          
          )
        
      
      ⁢
      
        d
        x
      
    
  




The SID measure may be the symmetric measure of the relative entropy measure. The SID measure may be combined with the spectral angle mapper (SAM). For the purpose of the present algorithm, one may need to use the SID measure over the other measure because it is conjugative (i.e., symmetric around its arguments). The relative entropy which is usually a logical choice in most of divergence measurements does not appear to satisfy this property and thus is not necessarily used in the present approach. With little manipulation, one may show that




  
    J
    ⁡
    
      (
      
        
          f
          i
        
        ,
        g
      
      )
    
  
  =
  
    
      
        3
        2
      
      ⁢
      
        
          (
          
            
              
                σ
                i
              
              
                σ
                g
              
            
            -
            
              
                σ
                g
              
              
                σ
                
                  i
                  ⁢
                  
                    
                  
                
              
            
          
          )
        
        2
      
    
    +
    
      
        1
        2
      
      ⁢
      
        
          (
          
            
              1
              
                σ
                g
                2
              
            
            +
            
              1
              
                σ
                i
                2
              
            
          
          )
        
        2
      
      ⁢
      
        
          (
          
            
              
                μ
                →
              
              g
            
            -
            
              
                μ
                →
              
              i
            
          
          )
        
        T
      
      ⁢
      
        (
        
          
            
              μ
              →
            
            g
          
          -
          
            
              μ
              →
            
            i
          
        
        )
      
    
  




In the related art, new scanned pixels data points may be modeled with fixed predefined distributions regardless of the application and operation conditions. The distribution may have had been assumed constant where these constants were based on some experimental observations. One may argue that the variance cannot be predefined as it varies with the operating environments; rather, one may provide better estimates of the incoming distribution that is dependent upon the current operations. One may formulate this approximation to simplify the J(ƒi,g) divergence measure and introduce a H(ƒi,g) divergence measure as more compact that results into similar results than the SID measure. Present estimates may be based on the fact that the incoming distribution is modeled as follows:

μg=xt and σg2=αiσi2 

where xt is the scanned data point. The choice of variance may be expressed as a scalar. Plugging this equation into the measure equation may result into:





  
    J
    ⁡
    
      (
      
        
          f
          i
        
        ,
        g
      
      )
    
  
  =
  
    
      
        3
        2
      
      ⁢
      
        
          (
          
            
              α
              i
            
            -
            
              1
              
                α
                i
              
            
          
          )
        
        2
      
    
    +
    
      
        1
        
          2
          ⁢
          
            
          
          ⁢
          
            σ
            i
            2
          
        
      
      ⁢
      
        
          (
          
            1
            +
            
              1
              
                α
                i
                2
              
            
          
          )
        
        2
      
      ⁢
      
        
          (
          
            
              
                μ
                →
              
              g
            
            -
            
              
                μ
                →
              
              i
            
          
          )
        
        T
      
      ⁢
      
        (
        
          
            
              μ
              →
            
            g
          
          -
          
            
              μ
              →
            
            i
          
        
        )
      
    
  




Based on the above assumptions, one may deduce the SID measure into a more compact information divergence measure. For instance, one may assume that αi=1; thus, the first term of the above measure is null, dropping the scalar 2 factor of second term as it would not affect the present comparison among all distributions when evaluated. This may result into a simple information measure (IM), given as follows:




  
    H
    ⁡
    
      (
      
        
          f
          i
        
        ,
        g
      
      )
    
  
  =
  
    
      1
      
        σ
        i
        2
      
    
    ⁢
    
      
        (
        
          
            
              μ
              →
            
            g
          
          -
          
            
              μ
              →
            
            i
          
        
        )
      
      T
    
    ⁢
    
      (
      
        
          
            μ
            →
          
          g
        
        -
        
          
            μ
            →
          
          i
        
      
      )
    
  




Once all of the divergence measures have been calculated one may find the distribution for which:





  
    H
    ⁢
    
      (
      
        
          f
          o
        
        ,
        g
      
      )
    
  
  =
  
    
      min
      i
    
    ⁢
    
      {
      
        H
        ⁡
        
          (
          
            
              f
              i
            
            ,
            g
          
          )
        
      
      }
    
  




and one may have a match between ƒo and g if and only if H(ƒi,g)≦κ, where κ is predefined cutoff value. This may apply for just the iris region matching. No thresholds are used if the match is in the non-iris region.

There may be a distribution update when there is a match. If the incoming distribution matches to one of the existing distributions, one may pool them together to a new normal distribution. This new normal distribution may be considered to represent the current state of the iris region or the non-iris region based upon the matched distribution and its associated category. In what follows, there may be described an effective approach to update the weights and the distribution parameters once there is a match ƒo. The weights may be updated as follows:




   
  
    {
    
      
        
          
            
              
                ω
                o
              
              ⁡
              
                (
                t
                )
              
            
            =
            
              
                
                  100
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      ω
                      o
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        1
                      
                      )
                    
                  
                
                +
                1
              
              101
            
          
        
      
      
        
          
            
              
                
                  ω
                  i
                
                ⁡
                
                  (
                  t
                  )
                
              
              
                i
                ≠
                o
              
            
            =
            
              
                100
                101
              
              ⁢
              
                
                  ω
                  o
                
                ⁡
                
                  (
                  
                    t
                    -
                    1
                  
                  )
                
              
            
          
        
      
    
  




One may make sure that the new weights still sum to 1,





  
    
      
        
          
            ∑
            
              
                ω
                i
              
              ⁡
              
                (
                t
                )
              
            
          
          =
          
          ⁢
          
            
              
                
                  ω
                  i
                
                ⁡
                
                  (
                  t
                  )
                
              
              +
              
                
                  ∑
                  
                    i
                    ≠
                    o
                  
                
                ⁢
                
                  
                    ω
                    i
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
              
            
            =
            
              
                100
                101
              
              ⁢
              
                
                  ∑
                  
                    i
                    ≠
                    o
                  
                
                ⁢
                
                  
                    ω
                    i
                  
                  ⁡
                  
                    (
                    
                      t
                      -
                      1
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        
          
          ⁢
          
            
              
                100
                101
              
              ⁢
              
                
                  ω
                  o
                
                ⁡
                
                  (
                  
                    t
                    -
                    1
                  
                  )
                
              
            
            +
            
              1
              101
            
          
        
      
    
    
      
        
          =
          
          ⁢
          
            
              
                100
                101
              
              ⁢
              
                (
                
                  
                    
                      ∑
                      
                        i
                        ≠
                        o
                      
                    
                    ⁢
                    
                      
                        ω
                        i
                      
                      ⁡
                      
                        (
                        
                          t
                          -
                          1
                        
                        )
                      
                    
                  
                  +
                  
                    
                      ω
                      o
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        1
                      
                      )
                    
                  
                
                )
              
            
            +
            
              1
              101
            
          
        
      
    
    
      
        
          
            =
            
            ⁢
            
              
                
                  100
                  101
                
                +
                
                  1
                  101
                
              
              =
              1
            
          
          ;
        
      
    
  
  +




To update the distribution parameters, one may introduce a learning parameter ε which is depicted as




  ɛ
  =
  
    1
    -
    
      
        H
        ⁡
        
          (
          
            
              f
              o
            
            ,
            g
          
          )
        
      
      
        K
        opt
      
    
  




Although it has been stated intuitively that Kopt may be inversely driving the selection of the learning parameter without formulating the relationship between the two variables, one may argue here that the actual driver of the learning parameter is the resulting divergence measure and how similar the incoming distribution is to the matched class. The constant Kopt should not be varying so as not to bias the overall measurement. One may here introduce a relationship between the two quantities and establish a guideline on how to set the learning parameter. This may also be a key contribution of our overall algorithm that works much more efficiently then many other approaches. Using the method of moments, one may get:

μi(t)=εμi(t−1)+(1−ε)μg 

σi2(t)≅εσi2(t−1)+(1−ε)σg2+ε(1−ε)(x(t)−μi(t−1))T(x(t)−μi(t−1))

while the unmatched distributions keep the same parameters as were previously set.

There may be a distribution update when there is a no-match. In the case where a match is not found, that is, where all divergence measures were found to be significantly bigger with respect to the Kopt value, no updates are conducted and a scanned pixel is considered to represent another cluster of pixels other than the iris region.
Relative to experiments and their results, the following may be noted. FIGS. 15a, 15b and 15c show a mixture modeling based analysis. FIG. 15a may here be regarded as a resulting rubber sheet map with respect to the pupil 11. FIG. 15b shows a mixture modeling based segmentation. For illustrative purposes, the iris 13 may be presented in red, the eyelid 17, 18 in green, and the eyelash 19 in blue. One may notice that the eyelash 19 covering the iris 13 may then be presented by the purple color as the mix of both red and blue. FIG. 15b may be regarded as revealing a segmentation of iris 13 versus eyelashes 19 and eyelids 17 and 18. FIG. 15c shows a deduced (binary) mask 71 representing the eyelashes 19 and eyelids 17 and 18. FIGS. 16a, 16b and 16c also show a mixture modeling based analysis. FIG. 16a shows an iris 13 normalized rubber sheet map. FIG. 16b shows a segmentation of iris 13 versus eyelashes 19 and eyelids 17 and 18. FIG. 16c shows a binary mask 72 representing the eyelashes 19 and eyelids 17 and 18.
In summary, there may an acquiring an image of an eye, an estimating a center of the pupil, an identifying a border of the pupil, an adjusting the center of the pupil, and a segmenting an iris from the border of the pupil to at least an outside border of the iris. One may further estimate a shape of the iris. In some examples, there may be an adjustment an eye image orientation perspective based on availability of face orientation angles. It may be noted that the terms “region” and “area” may be used interchangeably. These terms may be used interchangeably with “segment”.
The border of the pupil may be an inside border of the iris. The shape of the iris may be estimated as a regular shape while extracting the irregular actual shape of the inside border and the outside border of the iris. There may be two symmetric segments from a curvature perspective opposite of each other encompassing the pupil border and the outside border of the iris. There may be an identifying the border of the pupil with a snake delineation of the border of the pupil and the outside border of the iris. The snake delineation may be performed with one dimensional segmentation. The one dimensional segmentation may be performed in a polar domain originating at about the center of the pupil. The border of the pupil and the outside border of the iris may be irregular shapes. The irregular shapes may be normalized in a polar map.
The shape of the iris may be classified into unaffected regions and affected regions which can be normalized. The scaling for the normalizing of the affected regions may vary based on an obscuration by extraneous objects of the affected regions. One may identify symmetric segments of the unaffected regions. The affected regions may have obscured affected regions and non-obscured affected regions. There may be clustering to divide the affected regions into obscured affected regions and non-obscured affected regions. The clustering may be texture clustering using mixture modeling. The obscured affected regions may be masked.
The non-obscured affected regions may be a part of a fitting mechanism. An irregular actual shape may be estimated from unaffected regions and non-obscured affected regions as a regular shape. The irregular actual shape may be estimated from unaffected regions and non-obscured affected regions as a regular shape with a least square fitting approach applied to an elliptic quadratic curve. The least square fitting approach may include scaling the regular shape from ellipses to circles, performing a transform to determine at least one parameter, and scaling the circles back to ellipses.
An irregular actual shape may be estimated from unaffected regions and non-obscured affected regions as a regular shape with elliptic fitting using overlapping variable circles.
An irregular actual shape may be estimated from unaffected regions and non-obscured affected regions as a regular shape with a scaling mechanism to transfer an elliptic fitting approach into a circular fitting approach.
One may assign weighting to the unaffected regions and the affected regions. Hamming distance may be one of several approaches that may be used. The affected regions and the unaffected regions may be encoded into symmetric bins and asymmetric bins, respectively. There may be matching of the symmetric bins and asymmetric bins with symmetric bins and asymmetric bins of other irises. A weighting accorded a match of the symmetric bins may be significantly greater than a weighting accorded a match of the asymmetric bins. A weighting accorded a match of the symmetric bins with non-symmetric bins or vice versa, may be less than a weighting accorded a match of symmetric bins with symmetric bins or vice versa, and may be greater than a weighting accorded a match of the non-symmetric bins with non-symmetric bins.
FIG. 17 illustrates an example computer system usable in conjunction with certain illustrative instances of the invention. Computer system 100 may have processor(s) 102. The computer system 100 may also include a memory unit 130, processor bus 122, and input/output controller hub (ICH) 124. The processor(s) 102, memory unit 130, and ICH 124 may be coupled to the processor bus 122. The processor(s) 102 may have a suitable processor architecture. The computer system 100 may have one, two, three, or more processors, any of which may execute a set of instructions in accordance with illustrative examples of the present invention.
The memory unit 130 may include an operating system 140, which includes an I/O scheduling policy manager 132 and I/O schedulers 134. The memory unit 130 may store data and/or instructions, and may comprise any suitable memory, such as a dynamic random access memory (DRAM), for example. The computer system 100 may also include IDE drive(s) 108 and/or other suitable storage devices. A graphics controller 104 may control the display of information on a display device 106, according to the illustrative examples of the invention.
The input/output controller hub (ICH) 124 may provide an interface to I/O devices or peripheral components for the computer system 100. The ICH 124 may comprise any suitable interface controller to provide for any suitable communication link to the processor(s) 102, memory unit 130 and/or to any suitable device or component in communication with the ICH 124. For one instance of the invention, the ICH 124 may provide suitable arbitration and buffering for each interface.
For one example of the invention, the ICH 124 may provide an interface to one or more suitable integrated drive electronics (IDE) drives 108, such as a hard disk drive (HDD) or compact disc read-only memory (CD ROM) drive, or to suitable universal serial bus (USB) devices through one or more USB ports 110. For one instance, the ICH 124 may also provide an interface to a keyboard 112, a mouse 114, a CD-ROM drive 118, and one or more suitable devices through one or more firewire ports 116. The ICH 124 may also provide a network interface 120 though which the computer system 100 can communicate with other computers and/or devices.
In one example, the computer system 100 may include a machine-readable medium that stores a set of instructions (e.g., software) embodying any one, or all, of the methodologies for dynamically loading object modules described herein. Furthermore, software may reside, completely or at least partially, within memory unit 130 and/or within the processor(s) 102.
In the present specification, some of the matter may be of a hypothetical or prophetic nature although stated in another manner or tense.
Although the invention has been described with respect to at least one illustrative example, many variations and modifications will become apparent to those skilled in the art upon reading the present specification. It is therefore the intention that the appended claims be interpreted as broadly as possible in view of the prior art to include all such variations and modifications.

## Claims

1. A method of iris capture, comprising:
one dimensional segmentation of an iris into a map of the iris; and
classification of the map into unaffected areas and affected areas;
assigning weights to the unaffected areas and the affected areas of the map of the iris and an enrolled map of an iris; and
wherein:
higher weights are attributed to unaffected areas of the map of the iris and the enrolled map of an iris;
lower weights are attributed to affected areas of the map of the iris and the enrolled map of an iris; and
weights between the higher weights and the lower weights are attributed to unaffected and affected areas of the map of the iris and the enrolled map of an iris.

2. The method of claim 1, further comprising regular shape fitting of the areas for normalization.

3. The method of claim 2, further comprising identifying the unaffected areas as symmetric areas.

4. The method of claim 3, further comprising calibrating the affected areas.

5. The method of claim 4, further comprising clustering the affected areas.

6. The method of claim 5, further comprising assigning a weighting to the unaffected areas and the affected areas.

7. The method of claim 6, wherein the affected areas are affected by one or more eyelids and/or one or more eyelashes.