FR 2915604 A1 20081031 FR 0703039 A 20070426 PROCEDE ET DISPOSITIF DE LOCALISATION D'UN IRIS HUMAIN DANS UNE IMAGE La présente invention concerne le domaine du traitement d'images numériques, et plus particulièrement, les traitements appliqués à des images numériques de l'oeil humain dans des applications d'identification ou d'authentification. La reconnaissance d'iris constitue une technique d'identification biométrique éprouvée, pourvu que l'image sur laquelle sont appliqués les traitements d'analyse et d'identification soit exploitable. En particulier, la performance des algorithmes de reconnaissance dépend fortement de la netteté de l'image de l'iris à identifier. Or, dans la plupart des applications, et notamment dans les applications dites "embarquées" (par exemple, pour du contrôle d'accès à un téléphone ou ordinateur portable, pour une clé électronique, etc.), la caméra (capteur numérique et lentille) utilisée ne dispose pas de dispositif d'autofocus ou de mise au point automatique, ajustant la focale en fonction de la distance. De plus, pour obtenir une résolution suffisante de l'iris, les images sont prises à une distance relativement faible, généralement de l'ordre de 10 à 30 cm. Il en résulte une faible zone de netteté (plage de distance entre la caméra et l'oeil, dans laquelle l'image est nette). Cette faible profondeur de champ ajoutée au fait que l'oeil est sphérique peut engendrer des différences de netteté entre des zones d'une même image de l'oeil. L'utilisateur doit être très coopératif en se plaçant dans une zone de netteté relativement réduite. Généralement, l'utilisateur bouge pour se placer à une distance de netteté. En pratique, les images issues de systèmes d'acquisition numériques courants présentent les défauts suivants : ù apparition de tâches blanches provenant de réflexions spéculaires de sources lumineuses sur les obstacles existants dans l'environnement d'acquisition (lunettes, optique du capteur, ...), ù flou dû à de faibles mouvements de l'oeil (la pupille et l'iris étant toujours visibles), et à la qualité optique du capteur, ù contours de zones de l'image (pupille/iris, iris/cornée) moins nets, dus à un défaut de mise au point, et à la faible profondeur de champ du système optique, ù les cils couvrent une partie de l'iris. Lorsque les images sont de mauvaise qualité, les performances de certains procédés d'identification d'iris conventionnels peuvent diminuer, généralement en raison d'erreurs de localisation de l'iris. Un mode de réalisation de l'invention propose d'améliorer la localisation d'iris, notamment pour rendre les algorithmes d'identification moins sensibles à la qualité des images. Cet objectif est atteint par la prévision d'un procédé de localisation de l'iris dans une image d'un oeil, comprenant une étape de localisation de la pupille dans l'image. Selon un mode de réalisation, le procédé comprend des étapes de détection de positions de sauts d'intensités de pixels situés sur une ligne traversant la pupille et des zones de transition entre l'iris et la cornée, de part et d'autre de la pupille, et de détermination du centre et du rayon d'un cercle passant par les positions détectées des transitions. Selon un mode de réalisation, la détection des positions des transitions est effectuée en calculant l'énergie de composantes hautes fréquences d'un spectre de Fourier obtenu par une transformée de Fourier rapide sur une fenêtre mobile de type porte, appliquée aux intensités des pixels situés sur la ligne. Selon un mode de réalisation, le rayon de l'iris est déterminé en calculant la distance entre la position du centre de la pupille et l'une des positions sur la ligne, de transitions gauche et droite entre la cornée et l'iris. Selon un mode de réalisation, les positions sur la ligne des transitions gauche et droite entre la cornée et l'iris sont déterminées en recherchant deux positions où l'énergie haute fréquence est localement maximale. Selon un mode de réalisation, l'énergie E des composantes hautes fréquences est obtenue à l'aide de la formule suivante : E(m) = E feDHF Ax -2itifx II(x)IIm(x)e Ax x=1 2 dans laquelle I représente l'intensité d'un pixel à la position x sur la ligne, m représente la positon du centre de la fenêtre mobile de type porte IZm(x), et DHF représente le domaine des composantes hautes fréquences. Selon un mode de réalisation, la détection des positions des transitions est effectuée en appliquant un filtre gradient entre deux fenêtres mobiles successives de type porte, aux intensités des pixels situés sur la ligne. Selon un mode de réalisation, le rayon de l'iris est déterminé en calculant la distance entre la position du centre de la pupille et l'une des positions sur la ligne, des transitions gauche et droite entre la cornée et l'iris. Selon un mode de réalisation, les positions sur la ligne des transitions gauche et droite entre la cornée et l'iris sont déterminées en recherchant deux positions où la valeur du filtre gradient présente localement un extremum. Selon un mode de réalisation, le filtre gradient est calculé de la manière suivante : G(m)=ù I I(x) - - I I(x) XEI-Im XEIIm-1 dans laquelle I représente l'intensité d'un pixel à la position x, m et m-1 représentent la positon du centre des fenêtres mobiles de type porte fIm(x) et IIm_i(x), et Ax représente la largeur des fenêtres mobile. Selon un mode de réalisation, la largeur des fenêtres mobiles est égale à 16 pixels et le pas de déplacement des fenêtres mobiles est fixé à 4 pixels. Selon un mode de réalisation, la position du centre et la valeur du rayon de l'iris sont affinées par un opérateur intégro-différentiel. Selon un mode de réalisation, la localisation de la pupille est effectuée par une segmentation de l'image, pour délimiter le contour de la pupille, suivie d'une estimation du centre de la pupille par un calcul de centre géométrique du contour de la pupille. Selon un mode de réalisation, le rayon Rp de la 25 pupille est obtenu par l'expression suivante : Rp = max (w, h) (41)/4(1)2(4OY41 où w= ù et h= ù 01 01)1 7tl 01)2 avec (pi = E(xùX)2 et (1)2 = 1(yùy)2 ROIP ROIP x, y étant les coordonnées des pixels dans la zone de 30 recherche ROIP, et x, étant les coordonnées du centre de la pupille. Selon un mode de réalisation, la valeur du rayon de la pupille est affinée par un opérateur intégrodifférentiel. Selon un mode de réalisation, le procédé comprend une étape de prétraitement effectuant un filtrage des intensités des pixels de l'image pour éliminer les cils sur la cornée et l'iris, ainsi que les reflets sur la pupille. Selon un mode de réalisation, le filtrage effectué durant l'étape de prétraitement applique aux intensités des pixels de l'image des opérations de fermeture et d'ouverture morphologiques. Selon un mode de réalisation, l'étape de prétraitement comprend l'application d'un traitement de rehaussement de contraste aux intensités des pixels de l'image. Selon un mode de réalisation, le rehaussement de contraste est effectué uniquement sur une zone de l'image dans laquelle la pupille est supposée se trouver. L'invention concerne également un dispositif de localisation de l'iris dans une image d'un oeil, configuré pour localiser la pupille dans l'image. Selon un mode de réalisation, le dispositif est configuré pour détecter des positions de sauts d'intensités de pixels situés sur une ligne traversant la pupille et des zones de transition entre l'iris et la cornée, de part et d'autre de la pupille, et déterminer le centre et le rayon d'un cercle passant par les positions détectées des transitions. Selon un mode de réalisation, le dispositif est configuré pour détecter des positions des transitions en calculant l'énergie de composantes hautes fréquences d'un spectre de Fourier obtenu par une transformée de Fourier rapide sur une fenêtre mobile de type porte, appliquée aux intensités des pixels situés sur la ligne. Selon un mode de réalisation, le dispositif est configuré pour déterminer le rayon de l'iris en calculant la distance entre la position du centre de la pupille et l'une des positions sur la ligne, de transitions gauche et droite entre la cornée et l'iris, correspondant à l'énergie haute fréquence la plus élevée. Selon un mode de réalisation, le dispositif est configuré pour déterminer les positions sur la ligne des transitions gauche et droite entre la cornée et l'iris en recherchant deux positions où l'énergie haute fréquence est localement maximale. Selon un mode de réalisation, le dispositif est configuré pour calculer l'énergie E des composantes hautes fréquences à l'aide de la formule suivante : -2itifx 2 E(m) = E E I(x)II m (x)e Ax fEDHF x=1 dans laquelle I représente l'intensité d'un pixel à la position x sur la ligne, m représente la positon du centre de la fenêtre mobile de type porte IIm(x), et DHF représente le domaine des composantes hautes fréquences. Selon un mode de réalisation, le dispositif est configuré pour détecter les positions des transitions en appliquant un filtre gradient entre deux fenêtres mobiles successives de type porte, aux intensités des pixels situés sur la ligne. Selon un mode de réalisation, le dispositif est configuré pour déterminer le rayon de l'iris en calculant la distance entre la position du centre de la pupille sur la ligne et l'une des positions sur la ligne, des transitions gauche et droite entre la cornée et l'iris, correspondant à la valeur du filtre gradient la plus élevée en valeur absolue. Selon un mode de réalisation, le dispositif est configuré pour déterminer les positions sur la ligne des transitions gauche et droite entre la cornée et l'iris en recherchant deux positions où la valeur du filtre gradient présente localement un extremum. Selon un mode de réalisation, le dispositif est 5 configuré pour calculer le filtre gradient de la manière suivante : G(m) 1aI 11(x)- 1 11(x) / X2 XEIIm LxX XEIlm_1 dans laquelle I représente l'intensité d'un pixel à la position x, m et m-1 représentent la positon du centre 10 des fenêtres mobiles de type porte II11,(x) et 1-1ä1(x), et Ax représente la largeur des fenêtres mobile. Selon un mode de réalisation, la largeur des fenêtres mobiles est égale à 16 pixels et le pas de déplacement des fenêtres mobiles est fixé à 4 pixels. 15 Selon un mode de réalisation, le dispositif est configuré pour affiner la position du centre et la valeur du rayon de l'iris par un opérateur intégro-différentiel. Selon un mode de réalisation, le dispositif est configuré pour localiser la pupille par une segmentation 20 de l'image, pour délimiter le contour de la pupille, puis estimer le centre de la pupille par un calcul de centre géométrique du contour de la pupille. Selon un mode de réalisation, le dispositif est configuré pour calculer le rayon Rp de la pupille par 25 l'expression suivante : Rp = max (w, h ) '4PY'~2\ (4 ~/41 3 \) où w= ù et h= ù \u/ 01)1) (1)2) avec (1)1 = E(xùX)2 et (1)2 = E(yù@)2 ROIP ROIP x, y étant les coordonnées des pixels dans la zone de 30 recherche ROIP, et X, étant les coordonnées du centre de la pupille. Selon un mode de réalisation, le dispositif est configuré pour affiner la valeur du rayon de la pupille par un opérateur intégro-différentiel. Selon un mode de réalisation, le dispositif est configuré pour prétraiter l'image en effectuant un filtrage des intensités des pixels de l'image pour éliminer les cils sur la cornée et l'iris, ainsi que les reflets sur la pupille. Selon un mode de réalisation, le filtrage est effectué en appliquant aux intensités des pixels de l'image des opérations de fermeture et d'ouverture morphologiques. Selon un mode de réalisation, le prétraitement de l'image est effectué en appliquant un traitement de rehaussement de contraste aux intensités des pixels de l'image. Selon configuré uniquement pupille est supposée se trouver. Des exemples de réalisation de l'invention seront décrits dans ce qui suit, à titre non limitatif en relation avec les figures jointes parmi lesquelles : - la figure 1 représente sous forme de blocs un système de reconnaissance d'iris, - la figure 2 représente sous forme de blocs des modules du système de reconnaissance d'iris, selon un mode de réalisation de l'invention, - les figures 3a et 3b sont des images illustrant un prétraitement effectué par un module représenté sur la figure 2, - la figure 4 représente une courbe de rehaussement de contraste utilisée par le module de prétraitement, un mode de réalisation, le dispositif est pour dans rehausser le contraste de l'image une zone de l'image dans laquelle la - les figures 5a et 5b sont des images illustrant un procédé de localisation du contour de l'iris, selon l'invention - la figure 6 représente une courbe de variation 5 d'intensité de pixels de l'image le long d'un axe passant par la pupille, - la figure 7 représente une courbe de variation d'énergie de composantes hautes fréquences du spectre de Fourier de la courbe de la figure 6, et 10 - la figure 8 représente une courbe résultant de l'application d'un filtre gradient à la courbe de la figure 6. La figure 1 représente un système de reconnaissance d'iris IDS. Le système IDS est destiné à exploiter des 15 images d'oeil pour effectuer une identification ou authentification par reconnaissance iridienne. A cet effet, le système IDS comprend un capteur d'image numérique 1, un module de sélection d'image IMP, et un module de reconnaissance d'iris ID. Le capteur d'images 20 acquiert une séquence de trames vidéo ou d'images numériques SV d'un oeil 0 d'un utilisateur. Le module de sélection d'image IMP analyse en temps réel chaque trame ou image de la séquence d'images SV et sélectionne des images IT de la séquence SV ayant une qualité suffisante 25 pour effectuer une localisation de l'iris et une identification de l'utilisateur. Le module de reconnaissance d'iris ID utilise les images sélectionnées IT pour reconnaître les images d'iris. Le module ID comprend classiquement un module 30 d'extraction de signature d'iris ISE, et un module de comparaison MM connecté à une base de données de signatures d'iris IDB. Le module d'extraction ISE extrait une signature d'iris des images sélectionnées par le module IMP. Le module de comparaison MM recherche si la 35 signature d'iris fournie par le module ISE se trouve dans la base donnée IDB. Le module MM fournit un signal de détection DT indiquant si oui ou non l'iris de l'utilisateur a été reconnu. Par souci de clarté, seuls les modules ou étapes qui sont utiles à la compréhension de l'invention ont été représentés dans les figures et seront décrits par la suite. En particulier, les traitements de calcul de signature d'iris effectués par le module ISE et de reconnaissance d'iris effectués par le module MM n'ont pas été détaillés. D'une manière générale, le procédé d'extraction d'iris selon l'invention est compatible avec toute exploitation d'images d'iris. De même, le mode d'obtention des images à traiter par le procédé selon l'invention n'a pas été détaillé, l'invention étant là encore compatible avec n'importe quelle image numérique d'oeil, en niveau de gris ou en couleur. Par ailleurs, dans la description qui suit, on entendra par intensité d'un pixel le niveau de gris ou la luminance du pixel. La figure 2 représente un mode de réalisation selon l'invention de modules de localisation des contours de la pupille PDET et de l'iris IDET. Le module de localisation de la pupille PDET fait partie du module de sélection d'image IMP. Le module de localisation de l'iris IDET fournit au module d'extraction de signature d'iris ISE le centre et les rayons des contours interne (rayon de la pupille) et externe de l'iris dans les images sélectionnées IT de la séquence SV. Le module PDET comprend un module de prétraitement PRT, un module de délimitation d'une zone de recherche de la pupille SAD, un module de localisation de la pupille PLOC, et un module d'affinage PRD de la valeur du rayon Rp de la pupille. Le module PRT applique un prétraitement aux images 35 pour améliorer la détection de la pupille. Un tel prétraitement est par exemple décrit dans la demande de brevet EP 1 764 740. Ce prétraitement comprend tout d'abord une étape de filtrage permettant de supprimer des phénomènes indésirables produisant en général des erreurs de localisation. L'étape de filtrage comprend l'application d'opérations de fermeture et d'ouverture morphologiques à l'image. L'opération de fermeture morphologique (pl est effectuée par un élément structurant ayant la forme d'un disque dont le rayon en nombre de pixels correspond approximativement à la longueur d'un cil dans l'image. Dans une image de 320x240 pixels, ce rayon est par exemple de l'ordre de 10 pixels. Cette opération a pour effet d'éliminer les détails fins (de taille inférieure à celle de l'élément structurant) de l'image comme les cils, de lisser les contours des zones plus grandes que l'élément structurant, et de connecter certaines zones. L'opération d'ouverture morphologique Y1 est effectuée par un second élément structurant de la forme d'un disque dont la taille correspond sensiblement à celle des taches blanches des reflets de lumière sur l'oeil dans l'image. Ce disque présente un rayon compris entre deux et cinq fois celui du disque utilisé par l'opération de fermeture morphologique. Ainsi, dans une image de 320x240 pixels, ce rayon est par exemple de l'ordre de 20 pixels. L'opération d'ouverture morphologique a pour effet d'éliminer au moins en partie les taches blanches, de lisser les contours, mais préserve les zones sombres. Le filtrage M ainsi réalisé peut être modélisé par la formule suivante : M = yl • (pi (1) La figure 3a représente une image d'un oeil fournie par la caméra 1. La figure 3b représente la même image 35 après avoir subi le filtrage M. Il peut être constaté sur la figure 3b que les taches blanches présentes dans l'image de la figure 3a ont disparu. Il en est de même des cils qui ccuvraient en partie l'iris. I1 est à noter toutefois que ce traitement de filtrage a réduit le contraste entre les différentes régions de l'image regroupant des pixels contigus d'intensités voisines, c'est-à-dire notamment l'iris IR, la cornée CR et la pupille P. Ce traitement de filtrage a également diminué la circularité des régions de l'iris et de la pupille. La baisse de contraste peut être significative si l'intensité moyenne des images originales est faible. Un contraste insuffisant provoque des erreurs de localisation de la pupille et de l'iris. Une étape de rehaussement de contraste de l'image peut toutefois être réalisée. Pour limiter la durée de ce traitement, le rehaussement du contraste peut être effectué uniquement dans la région de l'iris IR et de la pupille P. A cet effet, le module PRT fait appel au une zone de recherche module SAD pour délimiter la pupille dans l'image. Le module SAD délimite laROIP de zone de recherche ROIP en comparant à un seuil l'intensité de chaque pixel de la région de faible intensité de l'image. Il calcule ensuite le centre géométrique et région. Pour vérifier la sombre) dans l'image, les déviations standard de cette présence de la pupille P (région le contraste de la zone de recherche ROIP est déterminé. Si le contraste est élevé, il est alors supposé que la pupille se trouve dans la zone de recherche, sinon, la pupille est considérée perdue et l'image ne sera pas sélectionnée pour effectuer une reconnaissance d'iris. Le contraste C d'une zone de l'image peut être obtenu par la formule suivante : Im ax ùlm in C= Im ax + lm in (2) dans laquelle Imax et Imin sont les intensités maximale et minimale dans la zone de recherche ROIP de l'image prétraitée. Le module SAD compare le contraste C à un seuil et s'il est supérieur au seuil, il considère que la pupille se trouve dans la zone de recherche ROIP. Dans le cas d'une séquence d'images, le module SAD recherche la pupille (région noire) dans une première image de la séquence SV, délimite la zone de recherche ROIP où la pupille P est supposée se trouver, et estime la position de la zone ROIP dans chaque image de la sequence en fonction d'une ou plusieurs images précédentes. Il existe plusieurs procédés de rehaussement de contraste d'une image. Toutefois, il s'agit toujours de diminuer les faibles intensités et d'augmenter les intensités les plus élevées, tout en conservant les formes initiales des régions de l'image. La figure 4 représente un exemple de courbe donnant l'intensité corrigée d'un pixel après rehaussement de contraste en fonction de l'intensité du pixel dans l'image à traiter. Entre 0 et une valeur de seuil sl d'intensité d'un pixel de l'image à traiter, l'intensité corrigée du pixel est fixée à une valeur minimum vmin, (par exemple égale à 20 pour 256 niveaux de gris). Entre la valeur sl et une valeur de seuil d'intensité s2, l'intensité corrigée du pixel varie entre la valeur minimum vmin et une valeur maximum d'intensité vmax, (par exemple fixée à 200 pour 256 niveaux de gris). La zone de la courbe entre les valeurs de seuil sl et s2 détermine une fonction d'étirement linéaire du contraste. Au-delà de la valeur s2, l'intensité corrigée du pixel est fixée à la valeur maximum vmax. Les valeurs de seuil sl et s2 peuvent être déterminées de manière adaptative en fonction de 35 l'intensité des pixels de l'image, de la façon suivante : sl = Imoy - 2a (3) (4) IN G= 1 (Ii ù Im oy)2 V-1 Ii étant l'intensité du pixel i dans la zone de recherche ROIP, a étant l'écart type des intensités des pixels dans la zone ROIP, et N étant le nombre de pixels dans la zone ROIP. La figure 5a représente une image fournie par le 10 module de prétraitement PRT. Les images prétraitées par le module de prétraitement PRT sont. fournies au module PLOC. Le module PLOC met en œuvre un procédé de localisation de la pupille. Le procédé de localisation de la pupille P peut 15 utiliser des opérateurs integro-différentiels permettant d'estimer le centre Cp(x,) de coordonnées estimées x,, et le rayon Rp de la pupille P. Pour plus d'informations sur ces opérateurs, on pourra se référer aux documents : "High confidence visual recognition of persons by a test 20 of statistical independence", J.G. Daugman, IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. 15, N 11, p. 1148-1161, 1993, et "High confidence personal identification by rapid video analysis of iris texture", J.Daugman, IEEE Conf. 25 Publication, n 408, European convention on security and detection, 16-18 May 1995. Les opérateurs integro-différentiels présentent de bonnes performances lorsque les images sont de bonne qualité (très peu de reflets lumineux dans l'image, fort 30 contraste entre la pupille et l'iris). Dans le cas contraire, le taux de bonne localisation diminue. En outre, la complexité algorithmique de ces opérateurs croît en fonction de la taille de l'image. Pour diminuer avec : N Im oy = E Ii (5) le temps de calcul de l'opérateur intégro-différentiel, le module PLOC peut utiliser la transformée de Hough décomposée suivant les gradients (GHT), combinée avec l'opérateur intégro-différentiel. La transformée GHT permet d'estimer rapidement la position du centre de la pupille, tandis que l'opérateur intégro-différentiel permet d'affiner la position du centre Cp (x, ) et le rayon Rp de la pupille au voisinage du centre estimé par la GHT. Pour plus d'informations sur ce procédé de localisation de la pupille, on pourra se référer au document : "Person identification technique using human iris recognition", C.L. TISSE, L. MARTIN, L. TORRES, M. ROBERT, The l5th International Conference on Vision Interface, S6.1, p. 294-300, May 27-29, 2002, Calgary, Canada 2002. Le procédé utilisant la transformée GHT pour localiser la pupille est plus rapide que le procédé basé sur un opérateur intégro-différentiel (complexité algorithmique plus faible) et présente de bonnes performances avec une image de bonne qualité (peu de reflets lumineux, faible présence de cils, région noire au voisinage de la pupille, fort contraste entre la pupille et l'iris, ...). Dans le cas contraire, ces performances diminuent. Afin de localiser efficacement la pupille P quelle que soit la qualité des images, avec une durée de traitement faible, le module PLOC selon un mode de réalisation, met en oeuvre un modèle de localisation basé sur la segmentation de la pupille P dans la zone de recherche ROIP. A cet effet, tous les pixels dans la zone ROIP de l'image prétraitée (filtrage + rehaussement de contraste) qui ont une intensité égale à si sont considérés comme appartenant à la pupille, et les autres pixels sont considérés comme n'appartenant pas à la pupille. L'image segmentée S peut ainsi être obtenue en affectant la valeur maximum 255 (blanc) aux pixels de la pupille et 0 (noir) aux autres pixels de la zone ROIP. La figure 5b est une image illustrant le résultat de la segmentation ainsi réalisée. Après avoir segmenté la pupille P dans la zone de recherche ROIP, le module PLOC estime la valeur du centre Cp(x,) de la pupille en calculant le centre géométrique des pixels (blancs) de la pupille dans l'image. Le rayon de la pupille Rp est ensuite obtenu à partir de l'expression suivante : Rp = max(w,h) (6) et h= (4\ (1)3 \n/ `1)2 avec 4) 1 = E(x-X)2 et (1)2 = E(y-3')2 ROIP ROIP (x,y) étant les coordonnées des pixels dans la zone de 15 recherche ROIP. Pour affiner la valeur du rayon Rp de la pupille P, le module d'affinage PRD peut appliquer l'opérateur intégro-différentiel à l'image prétraitée au voisinage des coordonnées (x,) du centre Cp et du rayon Rp de la 20 pupille, estimés précédemment. L'opérateur intégro- différentiel peut être défini par l'expression suivante : 0 I(x+rcos(8),+rsin(0))- 0=0 2 itr I(x + (r -1) cos(8), + (r -1) sin(8W (7) 27t(r -1) I(x,y) représentant l'intensité du pixel de coordonnées 25 (x, y) . Pour alléger les traitements d'image effectués par le module de sélection IMP, un traitement de sous- échantillonnage spatial peut être appliqué aux images de la séquence d'images SV avant leur traitement par le 30 module de prétraitement PRT. Le traitement de sous- où w= (Rp, x, ) = max rc[Rp-3,Rp+3] échantillonnage consiste à réduire la résolution des images d'un facteur par exemple égal à 2, pour passer d'images de 640x480 pixels à des images 320x240 pixels. Ce traitement de sous-échantillonnage affecte peu les performances du procédé de localisation selon l'invention. Une fois que la limite intérieure de l'iris IR a été localisée (par le centre Cp et rayon Rp de la pupille P), le module IDET localise la limite extérieure de l'iris, supposée circulaire, centrée au voisinage du centre Cp de la pupille. La localisation effectuée par le module IDET comprend donc la détermination du centre Ci et du rayon Ri du contour externe de l'iris. Pour déterminer le rayon externe Ri de l'iris IR, le module IDET peut, comme pour la pupille, utiliser l'opérateur intégro-différentiel. L'opérateur intégrodifférentiel est alors appliqué à une grille de petite taille (par exemple 5 x 5 pixels environ) centrée sur le centre Cp(x,) de la pupille P préalablement estimé (qui est supposé voisin du centre Ci de l'iris) et en supposant que le rayon de la pupille est situé dans une plage de valeurs comprise entre 2 et 4 fois le rayon de la pupille. Toutefois, cet opérateur ne permet pas d'estimer avec une grande précision le rayon Ri de l'iris et cette précision décroît lorsque le contraste entre l'iris IR et la cornée CR est faible. Par ailleurs, la complexité algorithmique de l'opérateur intégrodifférentiel croît en fonction de la largeur de la plage de valeurs estimée pour le rayon. Selon un mode de réalisation, ces deux problèmes sont surmontés par unprocédé d'estimation rapide du rayon Ri de l'iris IR. Ce procédé est basé sur la localisation des transitions entre l'iris IR et la cornée CR de part et d'autre de la pupille P. Ainsi, pour localiser le contour extérieur de l'iris, le module IDET se base sur l'utilisation d'un profil d'intensité des pixels situés sur un axe horizontal (Cp,x) passant par le centre Cp de la pupille P sur l'image prétraitée par le module PRT. La figure 5a représente l'axe (Cp,x), l'ordonnée y des pixels de cet axe étant fixée à l'ordonnée estimée @ du centre Cp de la pupille. D'une manière générale, l'axe (Cp,x) est choisi de manière à traverser les zones de transition entre l'iris IR et la cornée CR de part et d'autre de la pupille. La figure 6 représente un exemple de profil d'intensité I (x) sur l'axe (Cp, x) . Le profil d'intensité comporte des zones de pente élevée correspondant aux transitions gauches entre la cornée et l'iris TCI et entre l'iris et la pupille TIP, et aux transitions droites entre la pupille et l'iris TPI et entre l'iris et la cornée TIC. Selon un mode de réalisation, la localisation des transitions TOI, TIP, TPI, TIC est effectuée en appliquant au profil d'intensité I(x) la transformée de Fourier rapide (SFFT : Shifted Fast Fourier Transform) sur une fenêtre de type porte IIm(x) de centre m mobile, de largeur Ax, par exemple égale à 16 pixels, et dont le pas d'évolution est fixé par exemple à 4 pixels. Ensuite, l'énergie E(m) des composantes hautes fréquences du spectre de Fourier résultant est calculée à l'aide de la formule suivante : E(m) = E fEDHF ex -2itifx 2 1I(x)IIm (x)e Ax x=1 (8) dans laquelle DHF représente le domaine des composantes hautes fréquences [0, ..., Af/4], Af étant par exemple égal à 16. Dans cet exemple, le domaine des composantes hautes fréquences considéré correspond au quart supérieur du spectre de Fourier, c'est-à-dire aux quatre plus grandes valeurs. La figure 7 représente une courbe de variation de l'énergie des composantes hautes fréquences E(m) des pixels le long de l'axe (Cp,x). La transition gauche entre l'iris et la pupille TIP est détectée en une position ml sur l'axe (Cp, x) où l'énergie haute fréquence E(m) est maximale (pic central E(ml)). La transition gauche entre la cornée et l'iris TCI correspond à une position m2 où peut être observé un deuxième pic d'énergie E(m2) (autre maximum local). Le même raisonnement est appliqué pour localiser des positions m3 et m4 correspondant aux transitions droites entre la pupille et l'iris TIP et entre l'iris et la cornée TIC. Des pics d'énergie E(m3), E(m4) correspondant aux transitions droites TPI et TIC peuvent être observés aux positions m3 et m4. Les positions m2 et m4 sur l'axe (Cp,x) des transitons entre l'iris et la cornée TCI, TIC sont celles qui sont respectivement le plus à droite et le plus à gauche parmi les quatre positions ml-m4 sur l'axe (Cp,x), ou bien celles qui correspondent aux deux pics d'énergie les moins élevés parmi les quatre pics observés E(ml), E(m2), E(m3), E (m4) . Le rayon approximatif de l'iris Ri est obtenu à l'aide de l'expression suivante : Ril = l k - m4, Ri2 = I k - m2 l Si Ril Ri2 alors Si E(m4) > E(m2) alors Ri = Ril Sinon Ri = Ri2 Sinon Si E(m4) 20 TIC si les distances Rit et Ri2 sont proches. Dans le cas contraire, la valeur du rayon Ri est calculée avec la position m2 ou m4 correspondant au pic d'énergie E(m2), E(m4) le plus faible. Selon un autre mode de réalisation, la localisation des transitions TOI, TIC entre l'iris et la cornée le long de l'axe (Cp, x) passant par le centre de la pupille Cp est effectuée en appliquant un procédé dit de "filtrage gradient". Le procédé de filtrage gradient comprend l'application d'un filtre gradient G(m) entre deux fenêtres mobiles successives de type porte IIm_1(x), IIm(x) de centres m-1, m, espacées d'un pas d'évolution fixé par exemple à 4 pixels, chaque porte fIm(x) ayant une largeur Ax par exemple égale à 16. Le filtre gradient G(m) peut être exprimé de la manière suivante : G(m)=- E I(x) - 1 11(x) (10 ) ` x XEHm AxXEnm-1 La figure 8 représente un exemple de courbe de variation du filtre gradient G(m) appliqué aux pixels de l'axe (Cp,x). A la position ml correspondant à la transition gauche entre l'iris et la pupille TIP la courbe du filtre gradient G(m) présente une valeur minimum G(ml) (pic central négatif). A la position m2 correspondant à la transition gauche entre la cornée et l'iris TCI, la courbe présente un second pic négatif (autre valeur minimum locale) G(m2). A la transition droite entre la pupille et l'iris TPI, située à la position m3, la courbe G(m) présente une valeur maximum G(m3) (pic central positif). A la transition droite entre l'iris et la cornée TIC située à la position m4, un second pic positif G(m4) (autre valeur maximum locale) peut être observé sur la courbe du filtre gradient G(m). Ici encore, les positions m2 et m4 sur l'axe (Cp,x) des transitons entre l'iris et la cornée TCI, TIC sont celles qui sont respectivement le plus à droite et le plus à gauche parmi les quatre positions ml-m4, ou celles qui correspondent aux deux pics du filtre gradient les moins élevés en valeur absolue parmi les quatre pics observés G(ml), G(m2), G(m3), G(m4). Le rayon approximatif de l'iris est alors obtenu à l'aide de la formule suivante : Ril = I ?~ - m41, Ri2 = 1* - m2 l Si Ril Ri2 alors Si 1G(m4)1 > 1G(m2)I alors Ri = Ril Sinon Ri = Ri2 Sinon Si IG(m4)I Comme pour la détermination du rayon Rp de la pupille, la valeur du rayon Ri et la position du centre Ci de l'iris peuvent ensuite être affinées à l'aide de l'opérateur différentiel décrit ci-dessus, appliqué à l'image prétraitée. Le centre Ci(xi,i) de l'iris peut ainsi être obtenu avec une grille de taille 5 x 5 pixels centrée sur le centre de la pupille Cp(x,), et la valeur du rayon de l'iris est voisine de la valeur Ri précédemment calculée. Le centre Ci (xi,i) et le rayon Ri de l'iris peuvent être affinés à l'aide de la formule 35 suivante : (Ri, xi, i) = max rE Ri-4,Ri+8 XE X-2,X+2 yE @-2,+2 I(x+rcos(0),y+rsin(0)) ù 0=o 27tr I(x + (r ù 1) cos(8), y + (r -1) sin(0))l (12) 21t(r -1) Le procédé selon l'invention s'avère peu sensible aux occlusions de l'oeil par une paupière ou les cils, au 5 contraste entre l'iris et la cornée ou entre la cornée et la pupille, et à la présence de reflets lumineux parasites sur l'oeil. Il peut être mis en œuvre sans nécessiter d'importants moyens de calcul, tout en offrant un taux de bonne localisation proche de 100 %. Le procédé 10 de l'invention est donc adapté aux systèmes de reconnaissance d'iris à faible coût, c'est-à-dire ayant une caméra peu performante en ce qui concerne l'optique et le capteur d'image (sans dispositif d'autofocus, avec une faible profondeur de champ, une faible résolution 15 optique, .), et des moyens de calculs de faible capacité. Il apparaîtra clairement à l'homme de l'art que la..DTD: présente invention est susceptible de diverses variantes de réalisation et d'applications. En particulier, 20 l'invention n'est pas limitée à une localisation du contour extérieur de l'iris faisant intervenir un calcul d'énergie de composantes hautes fréquences ou d'un filtre gradient. D'autres procédés permettant de détecter des sauts d'intensités correspondant aux positions sur l'axe 25 (Cp, x) des transitions entre l'iris et la cornée peuvent être envisagés par l'homme de l'art. Il importe que cette détection soit effectuée parmi des pixels le long d'une ligne droite ou courbe, ou encore d'une ligne composée de segments de droite, traversant la pupille et des zones de 30 transition entre l'iris et la cornée, de part et d'autre de la pupille. Ainsi, il n'est pas nécessaire que cette ligne soit droite et passe exactement par le centre estimé de la pupille. Par ailleurs, d'autres procédés de localisation du contour de la pupille peuvent être envisagés, sans sortir du cadre de la présente invention. Il en est de même du procédé de prétraitement. Le procédé selon l'invention peut également être mis en œuvre sans prévoir de prétraitement préalable de l'image afin d'éliminer les cils et les reflets ou rehausser le contraste. L'absence de tout ou partie d'un tel prétraitement n'entraîne qu'une baisse de performances de la localisation de l'iris. REVENDICATIONS 1. Procédé de localisation de l'iris dans une image d'un oeil, comprenant une étape de localisation de la pupille dans l'image, caractérisé en ce qu'il comprend des étapes de détection de positions de sauts d'intensités de pixels situés sur une ligne (Cp,x) traversant la pupille et des zones de transition entre l'iris (IR) et la cornée (CR), de part et d'autre de la pupille (P), et de détermination du centre (Ci) et du rayon (Ri) d'un cercle passant par les positions détectées des transitions. 2. Procédé selon la revendication 1, dans lequel la détection des positions des transitions est effectuée en calculant l'énergie de composantes hautes fréquences d'un spectre de Fourier obtenu par une transformée de Fourier rapide sur une fenêtre mobile de type porte, appliquée aux intensités des pixels situés sur la ligne (Cp,x). 3. Procédé selon la revendication 2, dans lequel le rayon (Ri) de l'iris est déterminé en calculant la distance entre la position (x) du centre (Cp) de la pupille et l'une des positions (m2, m4) sur la ligne (Cp,x), de transitions gauche et droite entre la cornée et l'iris (TOI, TIC). 4. Procédé selon la revendication 3, dans lequel les positions (m2, m4) sur la ligne (Cp,x) des transitions gauche et droite entre la cornée et l'iris (TCI, TIC) sont déterminées en recherchant deux positions où l'énergie haute fréquence (E(m)) est localement maximale. 24 5. Procédé selon l'une des revendications 2 à 4, dans lequel l'énergie E(m) des composantes hautes fréquences est obtenue à l'aide de la formule suivante : -2nifx 2 E(m) = E E I(x)IIm (x)e x feDHF x=1 dans laquelle I(x) représente l'intensité d'un pixel à la position x sur la ligne (Cp,x), m représente la positon du centre de la fenêtre mobile de type porte IIm(x), et DHF représente le domaine des composantes hautes fréquences. 6. Procédé selon la revendication 1, dans lequel la détection des positions des transitions (TOI, TIP, TPI, TIC) est effectuée en appliquant un filtre gradient entre deux fenêtres mobiles successives de type porte, aux intensités des pixels situés sur la ligne (Cp,x). 7. Procédé selon la revendication 6, dans lequel le rayon (Ri) de l'iris (IR) est déterminé en calculant la distance entre la position (x) du centre (Cp) de la pupille et l'une des positions (m2, m4) sur la ligne (Cp,x), des transitions gauche et droite entre la cornée et l'iris (TCI, TIC). 8. Procédé selon l'une des revendications 6 et 7, dans lequel les positions (m2, m4) sur la ligne (Cp,x) des transitions gauche et droite entre la cornée et l'iris (TOI, TIC) sont déterminées en recherchant deux positions où la valeur du filtre gradient (G(m)) présente localement un extremum. 9. Procédé selon l'une des revendications 6 à 8, dans lequel le filtre gradient (G(m)) est calculé de la manière suivante .G(m)=ù /~j~ E I(x) ù 1 El(x) `-" XEIIm XEnm_1 dans laquelle I(x) représente l'intensité d'un pixel à la position x, m et m-1 représentent la positon du centre des fenêtres mobiles de type porte IZm(x) et llm_I(x), et Ax représente la largeur des fenêtres mobile. 10. Procédé selon l'une des revendications 2 à 9, dans lequel la largeur des fenêtres mobiles (IIm(x)) est égale à 16 pixels et le pas de déplacement des fenêtres 10 mobiles est fixé à 4 pixels. 11. Procédé selon l'une des revendications 1 à 10, dans lequel la position du centre (Ci) et la valeur du rayon (Ri) de l'iris sont affinées par un opérateur 15 intégro-différentiel. 12. Procédé selon l'une des revendications 1 à 11, dans lequel la localisation de la pupille (P) est effectuée par une segmentation de l'image, pour délimiter 20 le contour de la pupille, suivie d'une estimation du centre (Cp) de la pupille par un calcul de centre géométrique du contour de la pupille. 13. Procédé selon la revendication 12, dans lequel 25 le rayon Rp de la pupille est obtenu par l'expression suivante : Rp = max(w,h) où w=4 //2\/ et h=i4\//(i ~8 ^(1)1) \7ti 0t2) avec ~I = E(X-X)2 ROIPet 42 = E(YùY)2 ROIP 30 x, y étant les coordonnées des pixels dans la zone de recherche ROIP, et X, étant les coordonnées du centre (Cp) de la pupille. 14. Procédé selon la revendication 13, dans lequel la valeur du rayon (Rp) de la pupille est affinée par un opérateur intégro-différentiel. 15. Procédé selon l'une des revendications 1 à 14, comprenant une étape de prétraitement effectuant un filtrage des intensités des pixels de l'image pour éliminer les cils sur la cornée et l'iris, ainsi que les reflets sur la pupille. 16. Procédé selon la revendication 15, dans lequel le filtrage effectué durant l'étape de prétraitement applique aux intensités des pixels de l'image des opérations de fermeture et d'ouverture morphologiques. 17. Procédé selon la revendication 15 ou 16, dans lequel l'étape de prétraitement comprend l'application d'un traitement de rehaussement de contraste aux intensités des pixels de l'image. 18. Procédé selon la revendication 15 ou 17, dans lequel le rehaussement de contraste est effectué uniquement sur une zone de l'image (ROIP) dans laquelle la pupille est supposée se trouver. 19. Dispositif de localisation de l'iris dans une image d'un œil, configuré pour localiser la pupille dans l'image, caractérisé en ce qu'il est configuré pour détecter des positions de sauts d'intensités de pixels situés sur une ligne (Cp,x) traversant la pupille (P) et des zones de transition entre l'iris (IR) et la cornée (CR), de part et d'autre de la pupille (P), et déterminer le centre (Ci) et le rayon (Ri) d'un cercle passant par les positions détectées des transitions. 20. Dispositif selon la revendication 19, configuré pour détecter des positions des transitions en calculant l'énergie de composantes hautes fréquences d'un spectre de Fourier obtenu par une transformée de Fourier rapide sur une fenêtre mobile de type porte, appliquée aux intensités des pixels situés sur la ligne (Cp,x). 21. Dispositif selon la revendication 20, configuré pour déterminer le rayon (Ri) de l'iris en calculant la distance entre la position (x) du centre (Cp) de la pupille et l'une des positions (m2, m4) sur la ligne (Cp,x), de transitions gauche et droite entre la cornée et l'iris (TCI, TIC), correspondant à l'énergie haute fréquence (E(m3), E(m4)) la plus élevée. 22. Dispositif selon la revendication 21, configuré pour déterminer les positions (m2, m4) sur la ligne (Cp,x) des transitions gauche et droite entre la cornée et l'iris (TCI, TIC) en recherchant deux positions où l'énergie haute fréquence (E(m)) est localement maximale. 23. Dispositif selon l'une des revendications 20 à 22, configuré pour calculer l'énergie E(m) des composantes hautes fréquences à l'aide de la formule suivante : E(m) _ E feDHF Ox -2nifx E I(x)fl (x)e x=1 2 dans laquelle I(x) représente l'intensité d'un pixel à la position x sur la ligne (Cp,x), m représente la positon du centre de la fenêtre mobile de type porte Ilm(x), et DHF représente le domaine des composantes hautes fréquences. 24. Dispositif selon la revendication 19, configuré pour détecter les positions des transitions (TCI, TIP, TPI, TIC) en appliquant un filtre gradient entre deux fenêtres mobiles successives de type porte, aux intensités des pixels situés sur la ligne (Cp,x). 25. Dispositif selon la revendication 24, configuré pour déterminer le rayon (Ri) de l'iris en calculant la distance entre la position (x) du centre (Cp) de la pupille sur la ligne (Cp,x) et l'une des positions (m3, m4) sur la ligne (Cp,x), des transitions gauche et droite entre la cornée et l'iris (TCI, TIC), correspondant à la valeur du filtre gradient (G(m3), G(m4)) la plus élevée en valeur absolue. 26. Dispositif selon l'une des revendications 24 et 25, configuré pour déterminer les positions (m2, m4) sur la ligne (Cp,x) des transitions gauche et droite entre la cornée et l'iris (TCI, TIC) en recherchant deux positions où la valeur du filtre gradient (G(m)) présente localement un extremum. 27. Dispositif selon l'une des revendications 24 à 26, configuré pour calculer le filtre gradient (G(m)) de la manière suivante : G(m)=- /~}~ E I(x) - 1 E I(x) Ax xenm XEnm_1 dans laquelle I(x) représente l'intensité d'un pixel à la position x, m et m-1 représentent la positon du centre des fenêtres mobiles de type porte [lm(x) et IIm_1(x), et Ax représente la largeur des fenêtres mobile. 5 28. Dispositif selon l'une des revendications 20 à 27, dans lequel la largeur des fenêtres mobiles (IIm(x)) est égale à 16 pixels et le pas de déplacement des fenêtres mobiles est fixé à 4 pixels. 29. Dispositif selon l'une des revendications 19 à 28, configuré pour affiner la position du centre (Ci) et la valeur du rayon (Ri) de l'iris (IR) par un opérateur intégro-différentiel. 10 30. Dispositif selon l'une des revendications 19 à 29, configuré pour localiser la pupille (P) par une segmentation de l'image, pour délimiter le contour de la pupille (P), puis estimer le centre (Cp) de la pupille 15 par un calcul de centre géométrique du contour de la pupille. 31. Dispositif selon la revendication 30, configuré pour calculer le rayon Rp de la pupille par l'expression 20 suivante : Rp = max (w, h) 2 fi)" et h= \7t) \41 ~ avec 11)1 = E(xùX)2 ROIP x, y étant les coordonnées des pixels dans la zone. de 25 recherche ROIP, et x, étant les coordonnées du centre (Cp) de la pupille. 32. Dispositif selon la revendication 31, configuré pour affiner la valeur du rayon (Rp) de la pupille par un 30 opérateur intégro-différentiel. 33. Dispositif selon l'une des revendications 19 à 32, configuré pour prétraiter l'image en effectuant un où w= et (1)2 = ROIP filtrage des intensités des pixels de l'image pour éliminer les cils sur la cornée et l'iris, ainsi que les reflets sur la pupille. 34. Dispositif selon la revendication 33, dans lequel le filtrage est effectué en appliquant aux intensités des pixels de l'image des opérations de fermeture et d'ouverture morphologiques. 35. Dispositif selon l'une des revendications 33 et 34, dans lequel le prétraitement de l'image est effectué en appliquant un traitement de rehaussement de contraste aux intensités des pixels de l'image. 36. Dispositif selon la revendication 35, configuré pour rehausser le contraste de l'image uniquement dans une zone de l'image (ROIP) dans laquelle la pupille est supposée se trouver.20