FR 2906387 A1 20080328 FR 0608284 A 20060921 1 PROCEDE ET DISPOSITIF DE SELECTION D'IMAGES DANS UNE SEQUENCE D'IMAGES D'IRIS REÇUE EN FLUX CONTINU 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 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 une image 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 (réelle ou simulée) 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 profondeur de champ (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. En pratique, l'utilisateur bouge pour se placer à une distance de netteté. 2906387 2 Un traitement préalable à la reconnaissance d'iris proprement dite est donc souvent mis en œuvre pour sélectionner une image suffisamment nette. Ce traitement préalable consiste à sélectionner parmi un certain nombre 5 d'images, une ou plusieurs images nettes qui sont soumises à un algorithme de reconnaissance d'iris. La présente invention concerne plus particulièrement le traitement de sélection appliqué à des images d'un même œil pour sélectionner des images 10 suffisamment nettes à envoyer au dispositif de reconnaissance. Certains dispositifs de sélection de l'art antérieur enregistrent une séquence d'images, puis calculent un indice de netteté à l'aide d'un opérateur de 15 netteté, et finalement sélectionnent les images présentant les meilleurs indices de netteté. Parmi les opérateurs de netteté les plus couramment utilisés, on peut citer Tenegrad, SML (Sum Modified Laplacian), SMD (Sum Modulus Difference), et FSWM (Frequency Selective 20 Weighted median Filter). FSWM s'avère le plus insensible au bruit et le plus performant notamment en termes de ressources de calcul nécessaires. Ce mode de sélection présente l'avantage d'être rapide, et de ne pas augmenter le coût du système. 25 Cependant, l'acquisition d'une image d'iris exploitable est difficile à obtenir en pratique, car l'utilisateur n'est pas très coopératif. La distance de netteté (environ 25 cm) n'est pas toujours respectée. Très souvent, l'utilisateur bouge afin de se placer à la bonne 30 distance de netteté, et parfois il ferme promptement les paupières. Les images fournies par la caméra sont souvent de mauvaise qualité (images défocalisées, images floues dues à des mouvements trop rapides de l'oeil, à la qualité du système optique, présence de cils devant l'iris, présence de taches blanches dans la pupille dues aux 2906387 3 réflexions d'une source lumineuse, etc.). Les images les plus nettes sélectionnées par le dispositif de sélection sont donc souvent de qualité insuffisante pour localiser précisément l'iris dans l'image ou pour reconnaître 5 celui-ci. La présente invention vise à améliorer les procédés de sélection existants, notamment en réduisant le taux d'images non exploitables qui sont fournies au traitement d'analyse et d'identification. 10 L'invention vise également à proposer un procédé de sélection qui soit compatible avec les systèmes embarqués, notamment en terme de ressources de calcul utilisées. Ces objectifs sont atteints par la prévision d'un 15 procédé de sélection d'images d'un ensemble d'images en fonction de critères de netteté et de contraste. Selon un mode de réalisation de l'invention, le procédé comprend des étapes consistant à : û présélectionner des images par une analyse simplifiée 20 de netteté et/ou de contraste de chaque image de l'ensemble d'images, et û sélectionner des images présélectionnées par une analyse plus fine de la netteté et/ou du contraste de chaque image présélectionnée. 25 Selon un mode de réalisation de l'invention, la présélection d'images comprend des étapes consistant à : û calculer un contraste de chaque image de l'ensemble d'images, û comparer le contraste obtenu à un seuil de contraste, 30 et û rejeter l'image si le contraste obtenu est inférieur au seuil de contraste. Selon un mode de réalisation de l'invention, le contraste CS est calculé à l'aide de la formule 35 suivante : 2906387 4 CS = Imax ùImin 1m ax + 1m in dans laquelle : Imax = max f 1 G r (x, y) Imin =min iN ,lGt (x, y) 5 Gt(x,y) représente le niveau de gris ou la luminance du pixel de coordonnées x et y dans une région de l'image t. Selon un mode de réalisation de l'invention, l'ensemble d'images est une séquence d'images reçue en 10 flux continu, la présélection d'images comprenant des étapes consistant à : ù localiser un motif dans chaque image de la séquence, ù estimer une vitesse de déplacement du motif entre deux images successives de la séquence, et 15 ù sélectionner les images pour lesquelles la vitesse estimée de déplacement du motif est inférieure à un seuil de vitesse. Selon un mode de réalisation de l'invention, la localisation du motif dans chaque image comprend des 20 étapes consistant à : ù estimer la position du motif dans l'image en fonction de la position du motif dans une image précédente de la séquence, ù définir une zone de recherche du motif centrée sur la 25 position estimée du motif dans l'image, et ù localiser le motif dans l'image en se limitant à la zone de recherche. Selon un mode de réalisation de l'invention, l'estimation de la position du motif dans l'image est 30 effectuée en mettant en oeuvre le modèle de Kalman standard. Selon un mode de réalisation de l'invention, le motif est l'image de la pupille d'un oeil. 2906387 5 Selon un mode de réalisation de l'invention, la pupille est localisée dans la zone de recherche en recherchant des pixels ayant un niveau de gris ou une luminance inférieure à un seuil d'intensité, et en 5 déterminant le centre et le rayon d'un cercle suivant approximativement le contour des pixels trouvés. Selon un mode de réalisation de l'invention, la sélection des images comprend une étape de mesure de la densité des pixels noirs dans une fenêtre de localisation 10 de la pupille. Selon un mode de réalisation de l'invention, la sélection des images comprend une étape de mesure de la netteté de la texture de l'image de l'iris de l'oeil. Selon un mode de réalisation de l'invention, la 15 mesure de la netteté de la texture de l'iris est appliquée aux pixels de l'iris situés sur un anneau sensiblement circulaire centré sensiblement sur le centre de la pupille et ayant au moins un pixel de large. Selon un mode de réalisation de l'invention, la 20 mesure de la netteté de la texture de l'iris est effectuée à partir de l'énergie des composantes moyennes fréquences des pixels situés sur l'anneau. Selon un mode de réalisation de l'invention, la mesure de la netteté de la texture de l'iris dans une 25 image t est effectuée en calculant un indice de netteté I à l'aide de la formule suivante : EMF I(t)_ ET ù E max dans laquelle : EMF représente l'énergie des composantes de 30 moyennes fréquences d'un profil d'intensité d'iris des pixels situés sur l'anneau et localisés sur l'anneau par un angle, ET est l'énergie totale du spectre de Fourier du profil d'intensité, et 2906387 6 Emax est l'énergie maximale du spectre de Fourier du profil d'intensité. Selon un mode de réalisation de l'invention, la sélection des images comprend une étape de calcul d'un 5 indice de qualité ) fonction d'une mesure de la densité des pixels noirs de la pupille, d'une mesure de la netteté de la texture de l'iris et d'une mesure de la vitesse de déplacement de la pupille par rapport à l'image précédente dans la séquence d'image, et une étape 10 de comparaison de l'indice de qualité à un seuil de qualité, les images sélectionnées ayant un indice de qualité supérieur au seuil de qualité. Selon un mode de réalisation de l'invention, l'indice de qualité Q pour une image t est obtenu à 15 l'aide de la formule suivante : Q(t) = ad(t) +(1- a)I(t) + V(t) + 1 E dans laquelle : a et E sont des coefficients déterminés expérimentalement, et égaux de préférence à 0,5 et 0,0001 20 respectivement, d est la densité des pixels noirs dans la fenêtre de localisation de la pupille dans l'image t, I est un indice de netteté de l'iris de l'image t, et 25 V est la vitesse estimée de la pupille dans l'image t par rapport à l'image précédente. Selon un mode de réalisation de l'invention, le seuil de qualité est ajusté par une machine d'apprentissage. 30 Selon un mode de réalisation de l'invention, le procédé comprend une étape préalable de sous- échantillonnage spatial des images de l'ensemble d'images. 1 2906387 7 Selon un mode de réalisation de l'invention, le procédé comprend des étapes de mémorisation d'un nombre prédéfini d'images sélectionnées, et de sélection parmi les images mémorisées d'un nombre prédéfini d'images 5 ayant une netteté et/ou un contraste le plus élevé. L'invention concerne également un dispositif de sélection d'images d'un ensemble d'images en fonction de critères de netteté et de contraste. Selon un mode de réalisation de l'invention, le 10 dispositif comprend : ù un module de présélection d'images pour présélectionner des images par une analyse simplifiée de netteté et/ou de contraste de chaque image de l'ensemble d'images, et ù un module de sélection d'image pour sélectionner des 15 images présélectionnées par une analyse plus fine de la netteté et/ou du contraste de chaque image présélectionnée. Selon un mode de réalisation de l'invention, le module de présélection est configuré pour : 20 ù calculer un contraste de chaque image de l'ensemble d'images, ù comparer le contraste obtenu à un seuil de contraste, et ù rejeter l'image si le contraste obtenu est inférieur au 25 seuil de contraste. Selon un mode de réalisation de l'invention, l'ensemble d'images est une séquence d'images reçue en flux continu, le module de présélection étant configuré pour : 30 ù localiser un motif dans chaque image de la séquence, ù estimer une vitesse de déplacement du motif entre deux images successives de la séquence, et ù sélectionner des images pour lesquelles la vitesse estimée de déplacement du motif est inférieure à un 35 seuil de vitesse. 2906387 8 Selon un mode de réalisation de l'invention, le module de présélection est configuré pour : ù estimer la position du motif dans l'image en fonction de la position du motif dans une image précédente de la 5 séquence, ù définir une zone de recherche du motif centrée sur la position estimée du motif dans l'image, et ù localiser le motif dans l'image en se limitant à la zone de recherche. 10 Selon un mode de réalisation de l'invention, le module de présélection est configuré pour estimer la position du motif dans l'image en mettant en oeuvre le modèle de Kalman standard. Selon un mode de réalisation de l'invention, le 15 motif est l'image de la pupille d'un oeil. Selon un mode de réalisation de l'invention, le module de présélection est configuré pour localiser la pupille dans la zone de recherche en recherchant des pixels ayant un niveau de gris ou une luminance 20 inférieure à un seuil d'intensité, et en déterminant le centre et le rayon d'un cercle suivant approximativement le contour des pixels trouvés. Selon un mode de réalisation de l'invention, le module de sélection est configuré pour mesurer la densité 25 des pixels noirs dans une fenêtre de localisation de la pupille. Selon un mode de réalisation de l'invention, le module de sélection est configuré pour mesurer la netteté de la texture de l'image de l'iris de l'oeil. 30 Selon un mode de réalisation de l'invention, le module de sélection est configuré pour mesurer la netteté de la texture de l'image de l'iris de l'oeil dans un anneau sensiblement circulaire centré sensiblement sur le centre de la pupille et ayant au moins un pixel de large. 2906387 9 Selon un mode de réalisation de l'invention, le module de sélection est configuré pour mesurer la netteté de la texture de l'iris à partir de l'énergie des composantes moyennes fréquences des pixels situés sur 5 l'anneau. Selon un mode de réalisation de l'invention, le module de sélection est configuré pour calculer un indice de netteté I(t) à l'aide de la formule suivante : EMF 1(0 = ET - E max 10 dans laquelle : EMF représente l'énergie des composantes de moyennes fréquences d'un profil d'intensité d'iris des pixels situés sur l'anneau et localisés sur l'anneau par un angle, 15 ET est l'énergie totale du spectre de Fourier du profil d'intensité, et Emax est l'énergie maximale du spectre de Fourier du profil d'intensité. Selon un mode de réalisation de l'invention, le 20 module de sélection est configuré pour calculer un indice de qualité ) fonction d'une mesure de la densité des pixels noir de la pupille, d'une mesure de la netteté de la texture de l'iris et d'une mesure de la vitesse de déplacement de la pupille par rapport à l'image 25 précédente dans la séquence d'image, et une étape de comparaison de l'indice de qualité à un seuil de qualité, les images sélectionnées ayant un indice de qualité supérieur au seuil de qualité. Selon un mode de réalisation de l'invention, le 30 seuil de qualité est ajusté par une machine d'apprentissage. Selon un mode de réalisation de l'invention, le dispositif comprend un module de sous-échantillonnage qui sous-échantillonne les images de l'ensemble d'images et 2906387 10 qui envoie des images sous-échantillonnées au module de présélection. Selon un mode de réalisation de l'invention, le dispositif comprend une zone mémoire pour mémoriser un 5 nombre prédéfini d'images sélectionnées par le module de sélection fine et un module de sélection pour sélectionner parmi les images mémorisées un nombre prédéfini d'images ayant une netteté et/ou un contraste les plus élevés. 10 La présente demande couvre également un procédé et un dispositif de sélection d'images d'iris sur la base de deux paramètres, à savoir la netteté de la texture de l'iris et la densité de pixels noirs dans la pupille. Plus précisément, le procédé de sélection d'images 15 d'iris selon l'invention comprend pour chaque image d'iris des étapes de : û détermination d'une zone de localisation de la pupille dans l'image, û mesure de la densité des pixels noirs dans la zone de 20 localisation de la pupille, û mesure de la netteté de la texture de l'iris dans l'image, û calcul d'un indice de qualité de l'image fonction de la mesure de la densité des pixels noirs de la pupille, et 25 de la mesure de la netteté de la texture de l'iris, et û sélection de l'image si l'indice de qualité de l'image est supérieur à un seuil de qualité. Selon un mode de réalisation de l'invention, la pupille est localisée dans une zone de recherche en 30 recherchant des pixels ayant un niveau de gris ou une luminance inférieure à un seuil d'intensité, et en déterminant le centre et le rayon d'un cercle suivant approximativement le contour des pixels trouvés. Selon un mode de réalisation de l'invention, la 35 mesure de la netteté de la texture de l'iris est 2906387 11 appliquée aux pixels de l'iris situés sur un anneau sensiblement circulaire centré sensiblement sur le centre de la pupille et ayant au moins un pixel de large. Selon un mode de réalisation de l'invention, la 5 mesure de la netteté de la texture de l'iris est effectuée à partir de l'énergie des composantes moyennes fréquences des pixels situés sur l'anneau. Selon un mode de réalisation de l'invention, la mesure de la netteté de la texture de l'iris dans l'image 10 t est effectuée en calculant un indice de netteté I (t) à l'aide de la formule suivante : EMF I(t) = ET - E max dans laquelle : EMF représente l'énergie des composantes de 15 moyennes fréquences d'un profil d'intensité d'iris des pixels situés sur l'anneau et localisés sur l'anneau par un angle, ET est l'énergie totale du spectre de Fourier du profil d'intensité, et 20 Emax est l'énergie maximale du spectre de Fourier du profil d'intensité. Selon un mode de réalisation de l'invention, le procédé comprend une étape de mesure de la vitesse de déplacement de la pupille par rapport à une image 25 précédente dans une séquence d'image, l'indice de qualité étant fonction de la vitesse mesurée de la pupille. Selon un mode de réalisation de l'invention, l'indice de qualité Q(t) pour une image t est obtenu à l'aide de la formule suivante : Q(t) = ad(t) + (1- a)I(t) + V(t) + 1 E 1 30 dans laquelle : 2906387 12 a et E sont des coefficients déterminés expérimentalement, et égaux de préférence à 0,5 et 0,0001 respectivement, d(t) est la densité des pixels noirs dans la 5 fenêtre de localisation de la pupille dans l'image t, I(t) est un indice de netteté de l'iris de l'image t, et V(t) est la vitesse estimée de la pupille dans l'image t par rapport à l'image précédente. 10 Selon un mode de réalisation de l'invention, le seuil de qualité est ajusté par une machine d'apprentissage. Le dispositif de sélection d'image selon l'invention met en oeuvre le procédé de sélection d'image 15 tel que défini ci-avant. La présente demande couvre également un procédé et un dispositif d'estimation de la netteté de l'iris dans des images d'iris, sur la base d'une analyse spectrale de pixels de l'iris situés sur un anneau centré sur la 20 pupille. Plus précisément, le procédé selon l'invention d'estimation de la netteté de l'iris dans des images d'iris comprend pour chaque image des étapes de : û détermination du centre et du rayon de la pupille dans 25 l'image, û détermination d'un anneau d'au moins un pixel de large situé sur l'iris et centré sur le centre de la pupille, et û mesure de la netteté de la texture de l'iris à partir 30 de l'analyse spectrale des pixels de l'image situés sur l'anneau. Selon un mode de réalisation de l'invention, la pupille est localisée dans une zone de recherche en recherchant des pixels ayant un niveau de gris ou une 35 luminance inférieure à un seuil d'intensité, et en 2906387 13 déterminant le centre et le rayon d'un cercle suivant approximativement le contour des pixels trouvés. Selon un mode de réalisation de l'invention, la mesure de la netteté de la texture de l'iris est 5 effectuée à partir de l'énergie des composantes moyennes fréquences des pixels situés sur l'anneau. Selon un mode de réalisation de l'invention, la mesure de la netteté de la texture de l'iris dans l'image t est effectuée en calculant un indice de netteté I(t) à 10 l'aide de la formule suivante : EMF I(t) _ dans laquelle : EMF représente l'énergie des composantes de moyennes fréquences d'un profil d'intensité d'iris des 15 pixels situés sur l'anneau et localisés sur l'anneau par un angle, ET est l'énergie totale du spectre de Fourier du profil d'intensité, et Emax est l'énergie maximale du spectre de Fourier 20 du profil d'intensité. Le dispositif d'estimation de la netteté de l'iris selon l'invention met en oeuvre le procédé de sélection d'image tel que défini ci-avant. 25 Ces objets, caractéristiques et avantages ainsi que d'autres de la présente invention seront exposés plus en détail dans la description suivante d'un mode de réalisation de l'invention, faite à titre non limitatif en relation avec les figures jointes parmi lesquelles : 30 - 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 un dispositif de sélection selon un mode de réalisation de l'invention, ET ù E max 2906387 14 - les figures 3A, 3B et 3C sont des images d'oeil illustrant un procédé de sélection d'images d'iris selon un mode de réalisation de l'invention, - les figures 4 et 5 sont des courbes illustrant le 5 procédé de sélection d'images d'iris selon un mode de réalisation de l'invention, - la figure 6 est un exemple d'image d'iris défocalisée susceptible d'être traitée par le procédé de sélection selon un mode de réalisation de l'invention, 10 - la figure 7 est un organigramme illustrant différentes étapes du procédé de sélection selon un mode de réalisation de l'invention, - les figures 8 à Il sont des courbes illustrant les performances du procédé de sélection d'images d'iris. 15 - la figure 12 représente sous forme de blocs un dispositif de sélection selon un autre mode de réalisation de l'invention. La figure 1 représente un système de reconnaissance d'iris IDS. Le système IDS est destiné à exploiter des 20 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 qui acquiert une séquence de trames vidéo ou d'images numériques SV d'un oeil O d'un utilisateur, un 25 dispositif de sélection d'image IMP qui reçoit la séquence d'images SV, et un dispositif de reconnaissance d'iris ID utilisant les images sélectionnées IT par le dispositif IMP pour reconnaître les images d'iris. Le dispositif ID comprend classiquement un module 30 d'extraction de signature d'iris ISE qui extrait une signature d'iris des images sélectionnées par le module IMP, et un module de comparaison MM connecté à une base de données de signatures d'iris IDB, qui recherche si la signature d'iris fournie par le module ISE se trouve dans 35 la base donnée IDB. Le module MM fournit un signal de 2906387 15 détection DT indiquant si oui ou non l'iris du sujet a été reconnu. Par souci de clarté, seules les modules ou étapes qui sont utiles à la compréhension de l'invention ont été 5 représentés dans les figures et seront décrits par la suite. En particulier, les modules d'extraction de signature d'iris ISE et de reconnaissance d'iris MM exploitant le procédé de sélection de l'invention n'ont pas été détaillées, l'invention étant compatible avec 10 toute exploitation aval des images sélectionnées. De même, le mode d'obtention des images à traiter par 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. 15 La figure 2 représente un mode de réalisation selon l'invention d'un dispositif de sélection IMP. Le dispositif IMP analyse en temps réel chaque trame ou image d'une séquence d'images SV, et sélectionne les images IT ayant une qualité suffisante pour permettre une 20 localisation de l'iris et une identification de l'utilisateur. Le dispositif IMP comprend un module de détection et de localisation de la pupille DLP, un module de présélection ou de sélection rapide d'images QSEL, un module de sélection ou de sélection fine d'images FSEL, 25 et un module de suivi de la pupille PPP. Le module DLP recherche la pupille (région noire) dans chaque image de la séquence dans une zone de recherche dite "région d'intérêt" ROI où la pupille est supposée se trouver. La détection de la pupille dans une 30 image est effectuée en mesurant le contraste dans la zone de recherche. La région ROI est avantageusement estimée en fonction d'une ou plusieurs images précédentes par le module PPP. Les figures 3A et 3B sont des images d'oeil sur 35 lesquelles la région d'intérêt ROI a été représentée. 2906387 16 Dans l'exemple des figures 3A et 3B, la région ROI présente une forme carrée couvrant sensiblement l'ensemble de l'iris. Si le rapport signal sur bruit de l'image est 5 suffisant, l'algorithme suivant peut être utilisé pour détecter la pupille. Tout d'abord, une valeur de contraste CS est calculée à l'aide de la formule suivante : CS= ImaxùImin (1) 1m ax + 1m in 10 dans laquelle : Imax =maxN1 (x, y) lm in = min N 1 G r (x, y) Gt(x,y) représente le niveau de gris ou la luminance du pixel i de la région ROI de l'image t ayant les 15 coordonnées (x, y) dans l'image. La mesure de Imax est obtenue à partir des niveaux de gris ayant un niveau de gris inférieur à un seuil TH, afin de ne pas prendre en compte les pixels saturés issus de réflexion spéculaire d'une source lumineuse émettant 20 dans le proche infrarouge. La valeur de contraste CS obtenue est ensuite comparée à un seuil de contraste CT. On considère en effet que la pupille peut être détectée lorsque le contraste dans la région ROI est suffisamment élevé, par exemple de l'ordre de 0,6. Si le contraste de 25 l'image est inférieur au seuil CT, l'image est rejetée. Ensuite, si le contraste CS est supérieur au seuil de contraste CT, la pupille est localisée dans l'image t par le module DLP en positionnant dans l'image un cercle P (représenté en trait blanc dans les figures 3A, 3B) de 30 centre C(t) et de rayon R(t) suivant approximativement le contour de la pupille. La position du centre C(t) est déterminé en recherchant les pixels ayant un niveau de gris inférieur à un seuil d'intensité IT. (2) (3) 2906387 17 Le seuil d'intensité IT est déterminé en fonction de la répartition des niveaux de gris Gr(x,y) des pixels dans la région ROI par l'équation suivante : IT = u - 6 (4) 5 dans laquelle p représente la moyenne et 6 l'écart type de l'histogramme des niveaux de gris ou de la luminance Gt(x,y) des pixels dans la région ROI. L'usage d'un tel seuil adaptatif rend la détection de la pupille insensible aux changements de luminosité des images. 10 Le centre C(t) et le rayon R(t) de la pupille sont ensuite déterminés. A cet effet, chaque point de la zone d'intérêt ROI est considéré comme le centre potentiel C(t) et les gradients des points situés sur des arcs de cercle centrés sur le centre potentiel considéré sont 15 mesurés. Les rayons de ces arcs de cercle varient dans une plage de rayons possibles d'une pupille, et le centre de la pupille correspond au point pour lequel la variation du gradient est la plus importante. Les variations du gradient sont déterminés à l'aide 20 d'opérateurs intégro-différentiels. Ainsi, le rayon R(t) de la pupille est déterminé à l'aide d'opérateurs intégro-différentiels par l'équation suivante : Gt{xi(r,0),y;(r,e)} ôc iEC R(t) =Argmax (5) 25 dans laquelle : C représente le cercle de rayon R, R variant entre des valeurs minimum Rmin et maximum Rmax estimées du rayon de la pupille, Argmax représente une fonction fournissant le rayon 30 R ayant un score maximum par rapport à une formule normée. Autrement dit, le rayon R(t) est obtenu prenant la valeur maximum de l'intégrale de la norme du gradient au 2906387 18 rayon R-1 moins intégrale de la norme du gradient au rayon R. Selon une autre méthode dite "par segmentation", le centre C(t) est déterminé en calculant le centre 5 géométrique de la région pupille segmentée dans la région ROI. La segmentation de la pupille dans la région ROI est obtenue par seuillage adaptatif d'histogramme. Ainsi, tous les pixels qui auront un niveau de gris inférieur au seuil IT prendront la valeur 255 comme niveau de gris, et 10 les autres prendront la valeur 0, d'où la segmentation de la pupille. Un autre procédé de détermination du centre C(t) et du rayon R(t) de la pupille est exposé dans le document Tisse, et al. "Person identification technique using 15 human iris recognition", Journal of System Research, vol.4, pp.67-75, 2003. Le module PPP effectue uneprédiction de la position P(t/t-1) du centre de la pupille dans l'image suivante t notamment à partir de la position du centre 20 C(t-l) de la pupille dans l'image courante, et des précédentes mesures du centre de la pupille. A cet effet, il met en oeuvre le modèle de Kalman standard SKF tel que décrit notamment dans le document "Video Pupil Tracking for Iris based Identification", W. Ketchantang, 25 S. Derrode, S. Bourennane, and L. Martin, ACIVS 2005, LNCS3708, pp 1-8, 2005. Le modèle de Kalman Standard prédit l'état d'un système à partir des précédentes mesures effectuées sur ce dernier de telle sorte à minimiser la matrice de 30 covariance de l'erreur de prédiction. Une fois que la position prédite P(t/t-1) est déterminée, le module PPP positionne également la région d'intérêt ROI(t) dans l'image suivante t. La région ROI(t) est centrée sur la position prédite P(t/t-1) et 35 présente des côtés de longueur choisie de l'ordre de 2906387 19 quatre fois le rayon de la pupille R(t-1) dans l'image courante t-1 (en fonction de la dilatation de la pupille). Le module de sélection rapide QSEL effectue une 5 estimation de la vitesse V(t) de déplacement de la pupille entre l'image t et l'image précédente t-1. Si la vitesse ainsi déterminée est trop élevée (supérieure à un seuil de vitesse VT), on considère que l'image de l'iris dans l'image courante t est floue, et donc non 10 exploitable pour effectuer une identification. Par conséquent, l'image t est rejetée. La vitesse estimée de déplacement V(t) de la pupille est obtenue à partir de la position prédite P(t/t-1) du centre de la pupille dans l'image t et de la 15 position mesurée C(t) du centre de la pupille dans l'image t : v() =iP(t / t -1)C(tI T dans laquelle AT est égal à l'intervalle de temps entre deux images consécutives de la séquence d'images SV (par 20 exemple 1/25 s), et P(t/tû1)C(t)II représente la distance en nombre de pixels entre les points P(t/t-1) et C(t). Tout d'abord, le module FSEL localise dans chaque image sélectionnée par le module QSEL un anneau sensiblement circulaire situé dans l'iris et centré sur 25 le centre de la pupille, comme illustré sur la figure 3C. Sur la figure 3C, l'anneau IR présente un rayon interne égal à environ 1,5 fois le rayon R(t) de la pupille dans l'image t et une largeur d'au moins un pixel. Tous les pixels de l'iris situés sur l'anneau IR sont extraits 30 pour former un profil d'intensité de l'iris. Le profil d'intensité donne l'intensité D(t,O) (luminance ou niveau de gris) de chaque pixel extrait de l'image t en fonction d'un angle 0 déterminant la position du pixel sur (6) 2906387 20 l'anneau IR. Le profil d'intensité de l'iris caractérise la qualité de la texture de l'iris. La prise en compte uniquement des pixels situés sur un anneau permet dans certains cas de s'affranchir d'un 5 masquage partiel de l'iris par les cils ou une paupière. La figure 4 représente un exemple de profil d'intensité d'iris D(t,O) en fonction de l'angle O. Le module FSEL calcule ensuite la transformée de Fourier (Transformée de Fourier Rapide FFT) à chaque profil 10 D(t,O) d'image. Le profil utilisé pour le calcul de la transformée de Fourier est constitué d'un nombre de points en puissance de 2 pour accélérer le calcul, par exemple 256 points. Dans l'exemple de la figure 4, ces 256 points représentent la totalité du cercle (256 points 15 correspondent à 360 ). La courbe de variation de l'amplitude de la transformée de Fourier TF du signal D(t,O) est représentée sur la figure 5. La transformée de Fourier du signal D(t,O) comprend un pic très marqué centré sur la fréquence centrale (FO/2 = 128). Compte 20 tenu de la symétrie du spectre de Fourier, on travaille uniquement avec la moitié gauche du spectre. Le module FSEL calcule ensuite un indice de netteté en utilisant avantageusement les énergies des composantes de moyennes fréquences du spectre du signal D(t,O). Le 25 domaine des moyennes fréquences est par exemple choisi égal à FO/2 +/- FO/4, FO étant la fréquence centrale du spectre. L'indice de netteté de l'iris est calculé à l'aide de la formule suivante : I(t) = EMF (7) ET ù E max dans laquelle : EMF est l'énergie des composantes de moyennes fréquences du signal D(t,O), ET est l'énergie totale du spectre de Fourier, et 30 2906387 21 Emax est l'énergie maximale du spectre de Fourier, qui correspond à l'énergie à la fréquence centrale FO du spectre. Dans le domaine des fréquences moyennes, il s'avère 5 que l'indice de netteté obtenu est moins sensible au bruit et à la présence de cils devant l'iris. Il est à noter que la mesure de netteté de l'image obtenue à l'aide de la formule précédente est indépendante du contraste. En effet, si l'image G' est 10 moins contrastée que l'image G, D' = X D avec 0 mémorisées (M REVENDICATIONS 1. Procédé de sélection d'images d'un ensemble d'images en fonction de critères de netteté et de contraste, caractérisé en ce qu'il comprend des étapes 5 consistant à : ù présélectionner des images (t) par une analyse simplifiée de netteté et/ou de contraste de chaque image de l'ensemble d'images, et ù sélectionner des images présélectionnées par une 10 analyse plus fine de la netteté et/ou du contraste de chaque image présélectionnée. 2. Procédé selon la revendication 1, dans lequel la présélection d'images comprend des étapes consistant à : 15 ù calculer un contraste (CS) de chaque image (t) de l'ensemble d'images, ù comparer le contraste obtenu à un seuil de contraste (CT), et ù rejeter l'image si le contraste obtenu est inférieur au 20 seuil de contraste. 3. Procédé selon la revendication 2, dans lequel le contraste CS de chaque image est calculé à l'aide de la formule suivante : 25 CS = Imax ù Im in 1m ax + 1m in dans laquelle : Imax = max N 1 ut (x, y) Imin = miniv1 Gt (x, y) Gt(x,y) représente le niveau de gris ou la luminance 30 du pixel de coordonnées x et y dans une région de l'image t. 29 2906387 30 4. Procédé selon l'une des revendications 1 à 3, dans lequel l'ensemble d'images est une séquence d'images (t) reçue en flux continu, la présélection d'images comprenant des étapes consistant à : localiser un motif (P) dans chaque image de la séquence, estimer une vitesse de déplacement du motif entre deux images successives de la séquence, et sélectionner les images pour lesquelles la vitesse estimée de déplacement du motif est inférieure à un seuil de vitesse (VT). 5. Procédé selon la revendication 4, dans lequel la localisation du motif (P) dans chaque image comprend des 15 étapes consistant à : estimer la position du motif dans l'image en fonction de la position du motif dans une image précédente de la séquence, définir une zone de recherche (ROI) du motif centrée 20 sur la position estimée du motif dans l'image (t), et localiser le motif dans l'image en se limitant à la zone de recherche. 6. Procédé selon la revendication 5, dans lequel 25 l'estimation de la position du motif dans l'image est effectuée en mettant en oeuvre le modèle de Kalman standard. 7. Procédé selon l'une des revendications 1 à 6, 30 dans lequel le motif est l'image de la pupille d'un oeil. 8. Procédé selon la revendication 7, dans lequel l'image de la pupille est localisée dans la zone de recherche (ROI) en recherchant des pixels ayant un niveau 35 de gris ou une luminance inférieure à un seuil 2906387 31 d'intensité (IT), et en déterminant le centre (C) et le rayon (R) d'un cercle suivant approximativement le contour des pixels trouvés. 5 9. Procédé selon l'une des revendications 7 et 8, dans lequel la sélection des images comprend une étape de mesure de la densité des pixels noirs dans une fenêtre (P) de localisation de la pupille dans chaque image. 10 10. Procédé selon l'une des revendications 7 à 9, dans lequel la sélection des images comprend une étape de mesure de la netteté de la texture de l'image de l'iris de l'oeil dans chaque image. 11. Procédé selon la revendication 10, dans lequel la mesure de la netteté de la texture de l'iris est appliquée aux pixels de l'image de l'iris situés sur un anneau sensiblement circulaire (IR) centré sensiblement sur le centre de l'image de la pupille et ayant au moins un pixel de large. 12. Procédé selon la revendication 11, dans lequel la mesure de la netteté de la texture de l'image de l'iris est effectuée à partir de l'énergie des composantes moyennes fréquences des pixels situés sur l'anneau (IR). 13. Procédé selon la revendication 11 ou 12, dans lequel la mesure de la netteté de la texture de l'image de l'iris dans une image t est effectuée en calculant un indice de netteté I(t) à l'aide de la formule suivante : EMF I(t) = ET ù E max dans laquelle : 2906387 32 EMF représente l'énergie des composantes de moyennes fréquences d'un profil d'intensité d'iris des pixels situés sur l'anneau (IR) et localisés sur l'anneau par un angle, 5 ET est l'énergie totale du spectre de Fourier du profil d'intensité, et Emax est l'énergie maximale du spectre de Fourier du profil d'intensité. 10 14. Procédé selon l'une des revendications 7 à 13, dans lequel la sélection des images comprend une étape de calcul d'un indice de qualité (Q(t)) pour chaque image fonction d'une mesure de la densité des pixels noirs de l'image de la pupille, d'une mesure de la netteté de la 15 texture de l'image de l'iris et d'une mesure de la vitesse de déplacement de la pupille dans l'image par rapport à l'image précédente dans la séquence d'image, et une étape de comparaison de l'indice de qualité de l'image à un seuil de qualité (QT), les images 20 sélectionnées ayant un indice de qualité supérieur au seuil de qualité. 15. Procédé selon la revendication 14, dans lequel l'indice de qualité Q(t) pour une image t est obtenu à 25 l'aide de la formule suivante : Q(t) = ad(t) + (1- a)I(t) + 1 V(t) + ù dans E dans laquelle : a et E sont des coefficients déterminés expérimentalement, et égaux de préférence à 0,5 et 0,0001 30 respectivement, d(t) est la densité des pixels noirs dans la fenêtre (P) de localisation de la pupille dans l'image t, 5 2906387 33 I(t) est un indice de netteté de l'iris de l'image t, et V(t) est la vitesse estimée de la pupille dans l'image t par rapport à l'image précédente. 16. Procédé selon la revendication 14 ou 15, dans lequel le seuil de qualité (QT) est ajusté par une machine d'apprentissage (APM). 10 17. Procédé selon l'une des revendications 1 à 16, comprenant une étape préalable de sous-échantillonnage spatial des images de l'ensemble d'images. 18. Procédé selon l'une des revendications 1 à 17, 15 comprenant des étapes de mémorisation d'un nombre prédéfini d'images sélectionnées, et de sélection parmi les images mémorisées d'un nombre prédéfini d'images ayant une netteté et/ou un contraste le plus élevé. 20 19. Dispositif de sélection d'images d'un ensemble d'images en fonction de critères de netteté et de contraste, caractérisé en ce qu'il comprend : û un module de présélection d'images (QSEL) pour 25 présélectionner des images par une analyse simplifiée de netteté et/ou de contraste de chaque image (t) de l'ensemble d'images, et û un module de sélection d'image (FSEL) pour sélectionner des images présélectionnées par une analyse plus fine 30 de la netteté et/ou du contraste de chaque image présélectionnée. 20. Dispositif selon la revendication 19, dans lequel le module de présélection (QSEL) est configuré 35 pour : 2906387 34 calculer un contraste (CS) de chaque image (t) de l'ensemble d'images, comparer le contraste obtenu à un seuil de contraste (CT), et 5 rejeter l'image si le contraste obtenu est inférieur au seuil de contraste. 21. Dispositif selon la revendication 20, dans lequel le contraste CS d'une image est calculé à l'aide 10 de la formule suivante : CS= ImaxùImin Imax + lm in dans laquelle : Imax = max N 1 G t (x, y) lm in = min N1 Gt (x, y) 15 GI(x,y) représente le niveau de gris ou la luminance du pixel de coordonnées x et y dans une région de l'image t. 22. Dispositif selon l'une des revendications 19 à 20 21, dans lequel l'ensemble d'images est une séquence d'images (t) reçue en flux continu, le module de présélection (QSEL) étant configuré pour : localiser un motif (P) dans chaque image de la séquence, 25 estimer une vitesse de déplacement du motif entre deux images successives de la séquence, et sélectionner des images pour lesquelles la vitesse estimée de déplacement du motif est inférieure à un seuil de vitesse (VT). 30 23. Dispositif selon la revendication 22, dans lequel le module de présélection (QSEL) est configuré pour : 2906387 ù estimer la position du motif dans l'image en fonction de la position du motif dans une image précédente de la séquence, ù définir une zone de recherche (ROI) du motif centrée 5 sur la position estimée du motif dans l'image (t), et ù localiser le motif dans l'image en se limitant à la zone de recherche. 24. Dispositif selon la revendication 23, dans 10 lequel le module de présélection (QSEL) est configuré pour estimer la position du motif dans l'image en mettant en œuvre le modèle de Kalman standard. 25. Dispositif selon l'une des revendications 19 à 15 24, dans lequel le motif est l'image de la pupille d'un oeil. 26. Dispositif selon la revendication 25, dans lequel le module de présélection (QSEL) est configuré 20 pour localiser la pupille dans la zone de recherche (ROI) en recherchant des pixels ayant un niveau de gris ou une luminance inférieure à un seuil d'intensité (IT), et en déterminant le centre (C) et le rayon (R) d'un cercle suivant approximativement le contour des pixels trouvés. 25 27. Dispositif selon l'une des revendications 25 et 26, dans lequel le module de sélection (FSEL) est configuré pour mesurer la densité des pixels noirs dans une fenêtre (P) de localisation de la pupille. 28. Dispositif selon l'une des revendications 25 à 27, dans lequel le module de sélection (FSEL) est configuré pour mesurer la netteté de la texture de l'image de l'iris de l'oeil. 30 35 2906387 36 29. Dispositif selon la revendication 28, dans lequel le module de sélection (FSEL) est configuré pour mesurer la netteté de la texture de l'image de l'iris de l'oeil dans un anneau (IR) sensiblement circulaire centré 5 sensiblement sur le centre de l'image de la pupille et ayant au moins un pixel de large. 30. Dispositif selon la revendication 29, dans lequel le module de sélection (FSEL) est configuré pour 10 mesurer la netteté de la texture de l'image de l'iris à partir de l'énergie des composantes moyennes fréquences des pixels situés sur l'anneau (IR). 31. Dispositif selon la revendication 29 ou 30, 15 dans lequel le module de sélection (FSEL) est configuré pour calculer un indice de netteté I(t) à l'aide de la formule suivante : EMF I(t) = ET ù E max dans laquelle : 20 EMF représente l'énergie des composantes de moyennes fréquences d'un profil d'intensité d'iris des pixels situés sur l'anneau (IR) et localisés sur l'anneau par un angle, ET est l'énergie totale du spectre de Fourier du 25 profil d'intensité, et Emax est l'énergie maximale du spectre de Fourier du profil d'intensité. 32. Dispositif selon l'une des revendications 25 à 30 31, dans lequel le module de sélection (FSEL) est configuré pour calculer un indice de qualité (Q(t)) fonction d'une mesure de la densité des pixels noir de la pupille, d'une mesure de la netteté de la texture de l'iris et d'une mesure de la vitesse de déplacement de la 2906387 37 pupille par rapport à l'image précédente dans la séquence d'image, et une étape de comparaison de l'indice de qualité à un seuil de qualité (QT), les images sélectionnées ayant un indice de qualité supérieur au 5 seuil de qualité. 33. Dispositif selon la revendication 32, dans lequel l'indice de qualité Q(t) pour une image t est obtenu à l'aide de la formule suivante : t = ad(t) + (1- a)I(t) + 1 QO 1 dans laquelle : a et g sont des coefficients déterminés expérimentalement, et égaux de préférence à 0,5 et 0,0001 respectivement, 15 d(t) est la densité des pixels noirs dans la fenêtre (P) de localisation de la pupille dans l'image t, I(t) est un indice de netteté de l'iris de l'image t, et V(t) est la vitesse estimée de la pupille dans 20 l'image t par rapport à l'image précédente. 31. Dispositif selon la revendication 29 ou 30, dans lequel le seuil de qualité (QT) est ajusté par une machine d'apprentissage (APM). 25 34. Dispositif selon l'une des revendications 19 à 33, comprenant un module de sous-échantillonnage (SES) qui sous-échantillonne les images de l'ensemble d'images et qui envoie des images sous-échantillonnées au module 30 de présélection (QSEL). 35. Dispositif selon l'une des revendications 19 à 34, comprenant une zone mémoire (IMEM) pour mémoriser un 10 V(t) ± 2906387 38 nombre prédéfini d'images sélectionnées par le module de sélection fine (FSEL) et un module de sélection (MSEL) pour sélectionner parmi les images mémorisées un nombre prédéfini d'images ayant une netteté et/ou un contraste 5 les plus élevés.