FR 3040815 A1 20170310 FR 1558272 A 20150907 1 DOMAINE DE L'INVENTION L'invention concerne un procédé de caractérisation d'un matériau dans lequel une surface est formée, par l'analyse de tavelures générées sur la surface par un faisceau lumineux cohérent. ETAT DE LA TECHNIQUE L'analyse de tavelures générées sur la surface d'un matériau par un faisceau lumineux cohérent est déjà mise en oeuvre pour étudier des propriétés du matériau. Ainsi on connaît de la demande de brevet FR3007170 un procédé de détermination du caractère vivant ou non d'un tissu. Ce procédé comprend l'exposition de la surface à un faisceau lumineux cohérent pour générer sur ladite surface des tavelures. Les tavelures sont des figures d'interférences résultant des interférences du faisceau diffusé sur la surface. Une image des tavelures est acquise, puis cette image est analysée pour déterminer si le matériau constitutif de la surface est vivant ou non. Ce procédé a déjà rendu de grands services. Cependant, l'analyse qui est mise en oeuvre sur l'image des tavelures n'est pas suffisamment précise pour permettre de déterminer de quel matériau il s'agit. On connaît également du document W02006/085252 un procédé d'analyse de tavelures générées sur une surface, mais cette fois pour surveiller l'évolution temporelle de l'aspect des tavelures pour détecter l'existence d'un flux sanguin sous-jacent, et encore une fois décider si le tissu est vivant ou non. Ce procédé ne permet pas non plus de caractériser le matériau formant la surface. On connaît encore du document WO 2006/049396 un procédé d'analyse de tavelures pour déterminer si le matériau constitutif de la surface sur laquelle apparaissent les tavelures est un tissu biologique authentique ou une reproduction. Ce procédé se base sur l'analyse de nuances de gris dans une image de tavelures formées sur la surface. Ce procédé permet seulement de différencier un tissu biologique d'un tissu biologique et ne permet donc pas de différencier de nombreux matériaux, notamment des matériaux non biologiques. PRESENTATION DE L'INVENTION 3040815 2 L'invention a pour but de permettre de déterminer un matériau constitutif d'une surface par l'analyse de tavelures générées sur ladite surface. A cet égard, l'invention a pour objet un procédé de caractérisation d'un matériau 5 dans lequel est formée une surface, comprenant les étapes consistant à : - projeter sur la surface un faisceau lumineux cohérent, pour générer sur ladite surface des tavelures résultant d'interférences des rayons lumineux diffusés par ladite surface, - acquérir une image de ladite surface sur laquelle apparaissent les tavelures, 10 - traiter ladite image pour calculer au moins un critère de caractérisation du matériau, et - à partir des critères, déterminer le matériau constituant la surface, dans lequel l'étape de traitement de l'image comprend : - le calcul d'une fonction d'autocorrélation de l'intensité lumineuse sur l'image 15 et le calcul de la valeur d'au moins un critère établi à partir de ladite fonction, et le calcul de la valeur d'au moins un critère représentatif d'une distribution de déphasages des rayons lumineux diffusés par la surface. Avantageusement, mais facultativement, le procédé de caractérisation d'un 20 matériau selon l'invention peut en outre comprendre au moins l'une des caractéristiques suivantes : - le calcul de la valeur d'au moins un critère représentatif d'une distribution de déphasages des rayons lumineux diffusés peut comprendre le calcul d'une transformée en ondelettes de l'image, et le cas échéant le calcul de la valeur 25 d'au moins un critère parmi le groupe suivant : o moyenne des phases des coefficients d'ondelettes sur l'image, ou o écart-type des phases des coefficients d'ondelettes sur l'image. - le calcul de la valeur d'au moins un critère de caractérisation à partir de la fonction d'autocorrélation de l'intensité lumineuse peut comprendre : 30 o le calcul de la valeur d'au moins un critère établi à partir de la transformée de Fourier de la fonction d'autocorrélation, et o le calcul de la valeur d'au moins un critère établi à partir de la transformée en cosinus discrète de la fonction d'autocorrélation. 3040815 3 - un critère établi à partir de la transformée de Fourier de la fonction d'autocorrélation peut comprendre au moins l'un parmi le groupe suivant : o moyenne de la transformée de Fourier de la fonction d'autocorrélation sur l'image, ou 5 o écart-type de la transformée de Fourier de la fonction d'autocorrélation sur l'image. - le traitement de l'image peut en outre comprendre : o l'établissement d'une matrice d'autocorrélation de l'intensité lumineuse de l'image, et 10 o le calcul de la transformée en cosinus discrète de ladite matrice, et o le calcul de la valeur d'un critère de caractérisation comme la moyenne sur l'image de la trace de la matrice obtenue par la transformée en cosinus discrète de la matrice d'autocorrélation. - le traitement de l'image peut en outre comprendre l'établissement d'une 15 matrice d'autocorrélation de l'intensité lumineuse de l'image, et le calcul de la valeur d'un critère de caractérisation égal à la racine carrée de la somme des valeurs propres moyennes de ladite matrice sur l'image. - l'étape de détermination du matériau constituant la surface peut en outre comprendre le calcul d'un score global à partir des valeurs des critères, et la 20 classification du matériau par comparaison du score global avec une pluralité de plages prédéterminées correspondant à des matériaux respectifs. - Le procédé peut également comprendre une étape préliminaire d'élaboration des plages de scores, ladite étape comprenant : 25 o pour une pluralité de surfaces formées dans des matériaux déterminés, mettre en oeuvre des étapes : ^ de projection sur la surface d'un faisceau lumineux cohérent, pour générer sur ladite surface des tavelures résultant d'interférences des rayons lumineux diffusés par ladite 30 surface, ^ d'acquisition d'une image de ladite surface sur laquelle apparaissent les tavelures, ^ de traitement de ladite image pour calculer au moins un critère de caractérisation du matériau, ledit traitement 3040815 4 comprenant le calcul d'une fonction d'autocorrélation de l'intensité lumineuse sur l'image et le calcul de la valeur d'au moins un critère à partir de ladite fonction, le calcul de la valeur d'au moins un critère représentatif d'une distribution de 5 déphasages des rayons lumineux diffusés par la surface, et l'établissement d'un score global à partir desdites valeurs, o mettre en oeuvre un mécanisme d'apprentissage sur l'ensemble des scores globaux et l'ensemble des matériaux correspondants pour obtenir les plages correspondant à chaque matériau. 10 L'invention a également pour objet un procédé de contrôle biométrique, comprenant la mise en oeuvre du procédé la description qui précède sur une surface présentée par un individu lors du contrôle biométrique, le procédé comprenant en outre l'émission d'une alerte dans le cas où le matériau caractérisé 15 ne correspond pas à un tissu biologique. L'invention porte également sur un système de caractérisation d'un matériau comprenant : - une source de lumière cohérente, 20 - un capteur d'image, et - une unité de traitement, comprenant un processeur adapté pour : o traiter une image d'une surface sur laquelle apparaissent des tavelures, ledit traitement comprenant le calcul d'une fonction d'autocorrélation de l'intensité lumineuse sur l'image, 25 o calculer, à partir desdites images, la valeur d'au moins un critère de caractérisation du matériau formant la surface, comprenant au moins un critère calculé à partir de la fonction d'autocorrélation et au moins un critère représentatif d'une distribution de déphasages des rayons lumineux diffusés par la surface, et 30 o à partir des valeurs des critères, déterminer le matériau constituant la surface. L'invention a enfin pour objet un produit programme d'ordinateur, comprenant des instructions de code pour la mise en oeuvre d'un procédé 3040815 5 comprenant le traitement d'une image d'une surface pour calculer la valeur d'au moins un critère de caractérisation du matériau constitutif de la surface, ledit traitement comprenant : - le calcul d'une fonction d'autocorrélation de l'intensité lumineuse sur l'image 5 et le calcul de la valeur d'au moins un critère établi à partir de ladite fonction, et - le calcul de la valeur d'au moins un critère représentatif d'une distribution de déphasages des rayons lumineux diffusés par la surface, lorsqu'il est exécuté par un processeur. 10 L'invention proposée permet de caractériser le matériau formant une surface par l'analyse des tavelures générées sur la surface par un faisceau lumineux cohérent. En effet, des matériaux différents génèrent des tavelures de propriétés différentes. 15 L'utilisation d'au moins un critère basé sur l'autocorrélation de l'intensité lumineuse et d'au moins un critère basé sur les déphasages des rayons lumineux diffusés par la surface permet d'analyser des caractéristiques différentes de matériaux et ainsi de discriminer entre un nombre élevé de matériaux. En effet, l'autocorrélation de l'intensité lumineuse traduit un caractère diffusif 20 plus ou moins élevé du matériau (permettant de discriminer par exemple entre silicone et papier ou plâtre). La distribution de déphasages caractérise quant à elle le relief de la surface puisque le déphasage varie en fonction de la hauteur du point de la surface sur lequel le rayon a été diffusé. Elle permet donc de caractériser les matériaux selon le critère du relief (permettant de discriminer par exemple entre 25 papier et porcelaine). De plus, les critères basés sur l'autocorrélation de l'intensité lumineuse de l'image peuvent eux-mêmes refléter des caractéristiques différentes des tavelures de l'image et peuvent encore améliorer la précision de la caractérisation du procédé. Par exemple, des critères établis à partir de la densité spectrale de 30 puissance décrivent la distribution des tailles des figures de tavelures, tandis que d'autres critères peuvent décrire une distribution d'intensité des figures de tavelures. DESCRIPTION DES FIGURES 3040815 6 D'autres caractéristiques, buts et avantages de la présente invention apparaîtront à la lecture de la description détaillée qui va suivre, au regard des figures annexées, données à titre d'exemples non limitatifs et sur lesquelles : - La figure la représente schématiquement les principales étapes d'un 5 procédé de caractérisation de matériau selon un mode de réalisation de l'invention. - La figure 1 b représente schématiquement les principales étapes d'élaboration de plages de valeurs de score global pour caractériser un matériau, 10 - La figure 2 représente schématiquement un système pour la mise en oeuvre du procédé de caractérisation de matériau. - La figure 3 représente schématiquement un système d'identification ou d'authentification d'un individu. - La figure 4 représente un exemple de résultats obtenus en utilisant deux 15 critères basés sur la fonction d'autocorrélation dont un est obtenu à partir de la transformée en cosinus discrète de ladite fonction, et un critère représentatif de la distribution de déphasage. DESCRIPTION DETAILLEE D'AU MOINS UN MODE DE MISE EN OEUVRE DE 20 L'INVENTION En référence à la figure 1, on a représenté schématiquement les principales étapes d'un procédé de caractérisation d'un matériau constitutif d'une surface S. Ce procédé peut être mis en oeuvre lors d'un contrôle biométrique, au cours duquel une image d'un trait biométrique formé sur la surface est acquise. Typiquement, un trait 25 biométrique peut être une empreinte digitale ou la forme du visage. Le procédé a alors pour objectif de déterminer le matériau dans lequel la surface est formée, par exemple si elle est formée d'un tissu biologique, typiquement de l'épiderme, ou d'un matériau tel que le papier ou la silicone, pouvant être utilisé pour des fraudes. 30 Le procédé comprend une première étape 100 d'éclairage de la surface avec un faisceau lumineux cohérent, c'est-à-dire dont les rayons lumineux présentent une longueur de cohérence avantageusement supérieure à lcm, pour générer sur cette surface des tavelures. 3040815 7 Le procédé comprend en outre une deuxième étape 200, mise en oeuvre simultanément à la première, d'acquisition d'au moins une image de la surface S sur laquelle se trouvent les tavelures, de manière à ce que les tavelures apparaissent sur l'image. 5 On a représenté sur la figure 2 un système 1 adapté pour la mise en oeuvre du procédé de caractérisation du matériau. Pour la mise en oeuvre de l'étape 100 d'éclairage, ce système comprend une source 10 de lumière cohérente. Cette source est avantageusement un laser ou une diode laser monomode. Des sources lumineuses présentant une longueur de cohérence de l'ordre de 10 quelques dizaines de micromètres, telles que les LED, ou de l'ordre de quelques micromètres, comme la lumière blanche, sont proscrites. La source 10 lumineuse comprend avantageusement un polariseur adapté pour que le faisceau lumineux émis présente une polarisation rectiligne. L'éclairage de la surface S par le faisceau lumineux fait apparaître sur cette 15 surface des tavelures. Pour la mise en oeuvre de l'étape 200 d'acquisition d'image(s), le système 1 comprend avantageusement un capteur d'image 11 sensible à la même longueur d'onde que la lumière émise par la source 10, et comprenant un polariseur (non représenté) parallèle ou perpendiculaire à la direction de polarisation du faisceau 20 lumineux émis par la source 10. Le capteur d'image peut par exemple être une caméra numérique ou un appareil photographique numérique. Le système 1 comprend en outre de préférence une unité de traitement 12 comprenant typiquement au moins un processeur 120, et une mémoire 121, pour 25 réaliser le traitement 300 de l'image décrit ci-après. En référence à la figure 3, lorsque le procédé de caractérisation de matériau est mis en oeuvre lors d'un contrôle biométrique, le système 1 est avantageusement intégré à un système d'identification ou d'authentification 2. 30 Ce système 2 comporte un capteur de donnée biométrique 21, qui peut être distinct ou confondu avec le capteur d'image 11. Par exemple, lorsque le trait biométrique analysé est une empreinte digitale, le capteur d'image 11 peut être utilisé à la fois pour acquérir une image de l'empreinte et pour acquérir une image de tavelures générées sur le doigt où se trouve l'empreinte. 3040815 8 Le système d'identification ou d'authentification peut en outre comporter un lecteur 22 adapté pour obtenir, par lecture d'un document d'identité éventuellement numérique, une donnée biométrique de référence. Par exemple, la donnée biométrique peut être enregistrée dans une puce conservée dans le document 5 d'identité. Le lecteur 22 est dans ce cas adapté pour accéder au contenu mémorisé dans la puce. Alternativement, dans le cas par exemple ou le trait biométrique analysé est la forme du visage, une photographie de l'individu peut figurer dans le document d'identité et le lecteur 22 peut être un capteur d'image adapté pour acquérir une image de cette photographie. 10 Alternativement, le système d'identification ou d'authentification peut comporter une base de données 23 mémorisant un ensemble de données biométriques de référence. Ce système comporte en outre une unité de traitement 24, comprenant un processeur adapté pour comparer la donnée biométrique acquise sur l'individu avec 15 une donnée de référence obtenue à partir d'un document d'identité ou une donnée enregistrée dans la base de données. Avantageusement l'unité de traitement 24 est confondue avec l'unité 12, c'est-à-dire que la même unité réalise l'authentification ou l'identification et la détection de fraude par caractérisation du matériau formant la surface du trait biométrique analysé. 20 Dans le cas où l'identification ou l'authentification est mise en oeuvre pour contrôler l'accès physique de l'individu à une zone, le système peut comprendre un dispositif de restriction d'accès 25 pouvant être contrôlé par l'unité de traitement pour bloquer ou autoriser l'accès de l'individu à la zone. Le dispositif 25 est typiquement une barrière, une porte pouvant être verrouillée ou déverrouillée, etc. 25 Le système comprend enfin une alarme 26 pouvant être activée par l'unité de traitement 24 en cas de détection de fraude. L'alarme 26 peut être par exemple un signal sonore ou un signal visuel comme l'allumage d'une lumière ou l'affichage, sur un écran, d'un message d'alerte. 30 De retour à la figure 1, le procédé de caractérisation d'un matériau comprend, une fois une image acquise des tavelures réalisées sur la surface, une étape 300 de traitement de la ou des image(s) acquise(s) pour calculer des valeurs de critères statistiques de caractérisation du matériau constitutif de la surface. Cette 3040815 9 étape est avantageusement mise en oeuvre par l'unité de traitement 12 du système 1. Critères basés sur la fonction d'autocorrélation 5 L'étape 300 de traitement est décrite dans la suite pour une image. Dans le cas où plusieurs images de la surface sont acquises à l'étape 200, les images peuvent être traitées de manière identique, et les valeurs calculées lors du traitement 300 sont dans ce cas avantageusement moyennées sur l'ensemble des images traitées. 10 Le traitement d'une image comprend le calcul 310 de la fonction d'autocorrélation spatiale de l'intensité lumineuse sur l'image. Cette fonction est un paramètre statistique qui fournit une indication sur la taille des grains dans les tavelures et sur leur distribution. La fonction d'autocorrélation s'écrit comme suit : p(1, m) = ((I (x, y) - (I))I (x + 1, y + m) - (I))) o-2 Où a est l'écart-type de l'intensité lumineuse dans l'image, (I) est la valeur moyenne 15 dans l'image de l'intensité lumineuse I, (x,y) définit la position d'un pixel dans l'image, et I et m traduisent une distance dans les directions de x et y. A partir de la fonction d'autocorrélation, différents critères statistiques peuvent être utilisés. Avantageusement, les critères statistiques sont calculés uniquement sur les 20 pixels de l'image correspondant à une figure de tavelure. Pour identifier ces pixels, on calcule de préférence la surface s couverte par les tavelures, qui est définie comme l'intégrale double de la fonction d'autocorrélation : +00 +00 s= f f p(1, m) dldm -00 -00 Un premier critère de caractérisation de matériau peut être défini à partir d'une matrice A dite d'autocorrélation. Cette matrice A est obtenue en reformulant la 25 fonction d'autocorrélation comme suit, pour de petits décalages I et m : p(l,m) = [I (x , y) - I (x + x, y + 3,31)]2 (x,y)Ew Où w est l'ensemble des pixels de l'image. En appliquant l'expansion de Taylor à /(x + Ax, y + iy) on obtient : 3040815 10 a I a I 1(x+ Ax, y + Ay) = I (x, y) + -d x Ax + -dy Ay ai ai En posant /x = -ax et Iy = Ty, on a : 1(x + Ax, y + Ay) = 1(x, y) + + 1,4 Et en effectuant le changement de variable X=(x y), on obtient : 1(x, y) + + 1,4 = 1(X) + (V' I)AX Où w = (Ix,/y)T D'où : p(i, m) = [/xAx + /y3,37]2 = [/2x,Ax + 21x13,AxAy + 12 3A237] = (AX)T A (X)(AX) 5 Avec la matrice A d'autocorrélation définie en tout point (x,y) de l'image par : /2x lx1y1 A =w Ilxly 112 y Les valeurs propres notées À+ et À- de cette matrice indiquent le degré d'autocorrélation d'un pixel (x,y) par rapport aux pixels voisins. Les valeurs de À+ et À- sur l'image permettent donc d'obtenir la distribution d'intensité dans le domaine w, qui n'est pas le même pour différents types de matériaux. En particulier, un tissu 10 biologique présente une autocorrélation plus forte qu'un matériau utilisé pour une fraude. Par conséquent un premier critère de caractérisation Ci peut être défini à partir des valeurs propres de la matrice A. Avantageusement, ce critère est défini comme : C1= .\//1+ + A_ 15 Où A+ est la moyenne des valeurs propres À+ sur l'image, et A_est la moyenne des valeurs propres À- sur l'image. Ainsi le procédé comprend avantageusement une étape 311 de calcul de la valeur du critère Ci pour l'image acquise à l'étape 200. 20 Alternativement, ou en cumul avec la valeur du première critère de caractérisation Ci, la valeur d'un deuxième critère 02 peut être calculée à partir de la matrice A. 3040815 11 Pour calculer la valeur de ce deuxième critère, le procédé comprend une étape supplémentaire 312 de calcul de la transformée en cosinus discrète (ou DCT pour Discrete Cosine Transform) de la matrice d'autocorrélation A. Cette étape modifie la matrice A en une matrice A' de Toeplitz, c'est-à-dire 5 une matrice dont les coefficients situés sur chaque diagonale sont les mêmes. En l'espèce la matrice A' est donc de la forme : [a bl I-c a De ce fait, les valeurs propres de la matrice A' correspondent au minimum et maximum de puissance du spectre de l'intensité I au pixel (x,y) ; les valeurs propres de A' sont maximales quand l'autocorrélation dans l'image est forte. 10 On peut donc définir le deuxième critère 02 à partir de la trace de la matrice A', qui est une approximation de la somme des valeurs propres de A'. Ce critère peut par exemple être la moyenne sur l'image de la trace de A'. Plus cette moyenne est élevée, et plus l'autocorrélation dans l'image est élevée. Ainsi le procédé comprend avantageusement le calcul 313 de la valeur du 15 critère 02 à partir de l'image acquise à l'étape 200. Alternativement ou en cumul aux critères précédents, un critère de caractérisation établi à partir de la fonction d'autocorrélation peut être fourni par la transformée de Fourier de la fonction d'autocorrélation. Dans ce cas, l'étape de 20 traitement 300 comprend une étape 314 de calcul de la transformée de Fourier de la fonction d'autocorrélation de l'intensité lumineuse, notée F(p(1,m)). Or, d'après le théorème de Weiner-Khinchin, la transformée de Fourier de la fonction d'autocorrélation est égale à la densité spectrale de puissance de l'intensité : F(p(1,m)) = DSP (I) 25 La densité spectrale de puissance décrit la distribution des tailles des tâches dans la figure formée par les tavelures. Elle est donc exploitée pour caractériser les matériaux. On définit un troisième critère de corrélation C3 comme l'écart-type de la densité spectrale de puissance de l'intensité sur l'image, qui est calculé comme 30 l'écart-type de la transformée de Fourier de la fonction d'autocorrélation sur l'image. 3040815 12 On définit un quatrième critère de corrélation C4 comme la moyenne de la densité spectrale de puissance de l'intensité sur l'image, qui est calculée comme la moyenne de la transformée de Fourier de la fonction d'autocorrélation sur l'image. Le procédé comprend avantageusement une étape de calcul de la valeur 5 d'au moins un critère établi à partir de la transformée de Fourier, soit le critère C3 (étape 315), soit le critère C4 (étape 316). Le ou les critères décrits ci-avant établis à partir de la fonction d'autocorrélation sont avantageusement combinés à au moins un critère 10 supplémentaire représentatif de la distribution des déphasages des rayons lumineux diffusés par la surface. En d'autres termes, le procédé comprend le calcul de la valeur d'au moins un des critères décrits ci-avant, et de la valeur d'au moins un critère défini ci-après relatif au déphasage des rayons lumineux. 15 Critères de distribution de déphasages L'étape 300 de traitement de l'image comprend une étape 320 de la transformée en ondelettes de l'image afin de caractériser les déphasages des rayons lumineux diffusés par la surface. 20 En effet, les défauts de surface induisent sur les rayons lumineux des déphasages qui sont fonction des hauteurs des défauts. On peut donc exprimer le champ U en un point donné P comme une somme discrète des contributions des rayons lumineux provenant d'un ensemble de N défauts, les rayons lumineux provenant d'un même défaut présentant tous le 25 mêmes déphasage (pi, : U(P) = ane i(Pn n=1 Où a, désigne l'amplitude du champ des faisceaux lumineux diffusés par le défaut n. Lorsque le nombre de défauts N est important, le théorème de la limite centrale s'applique à l'équation précédente, et le champ en un point d'observation 30 U(P) suit une loi gaussienne ; les tavelures sont qualifiées de gaussiennes. 3040815 13 La transformée en ondelettes 320 de l'image permet d'estimer les valeurs du champ a, et du déphasage cp, en choisissant une échelle e adaptée pour que le nombre de défauts N considéré dans l'image soit suffisamment important pour que le champ suive une loi gaussienne. 5 Pour modéliser le champ gaussien, en adaptant les travaux de Mallat et al. dans la publication Characterization of signais from multiscale edges, IEEE Transaction on Pattern Analysis and Machine Intelligence, 14 : 710-732, 1992, on définit une fonction gaussienne e(x,y), et deux fonctions d'ondelettes 'pl et 1p2 comme étant respectivement les dérivées en x et en y de la fonction 0 : dB (x, y) Ipt oc, = dy 10 Pour une image présentant une intensité lumineuse I(x, y) E L2 (R2), la transformée en ondelettes à l'échelle e de l'intensité lumineuse I de l'image a deux composantes : I*(edee(x, y)) = ed(0, * I)(x, y) dx dx 1-) = I * tpe2 (x, y) = (edee(x, y)) = ed (Be * I)(x, y) dy dy Ainsi à une échelle e, l'amplitude de champs en (x,y) est définie par : a = ,/(1)D2 (e 2 Et la phase est définie par : (f) = tan- (Dj- 15 On définit alors un cinquième critère de corrélation C5 comme l'écart-type des phases des coefficients d'ondelettes sur l'image. On définit enfin un sixième critère de corrélation C6 comme la moyenne des phases des coefficients d'ondelettes sur l'image. Le procédé comprend donc avantageusement le calcul de la valeur d'au 20 moins l'un des critères C5 (étape 321), C6 (étape 322), ou des deux critères. Ainsi lors d'une étape 400, le processeur détermine la nature du matériau à partir des valeurs des critères calculées à l'étape précédente, et notamment à partir dx dB (x, y) lp 2 Oc, = 3040815 14 d'au moins l'un des quatre critères obtenus à partir de la fonction d'autocorrélation de l'intensité lumineuse, et d'au moins l'un des deux critères présentés obtenus à partir de la transformée en ondelettes de l'intensité lumineuse de l'image. De préférence, pour plus de précision dans la caractérisation du matériau, 5 les valeurs de tous les critères décrits ci-avant sont calculées et prises en compte dans l'étape 400. L'étape de caractérisation 400 du matériau comprend d'abord l'agrégation 410 par le processeur 120 des valeurs des critères obtenues pour former un score global. Le score global peut par exemple être un vecteur comprenant l'ensemble 10 des valeurs des critères. Alternativement, le score global peut être une combinaison linéaire des valeurs des critères ou, plus préférablement, un polynôme évalué en ces valeurs. Puis, le processeur 120 compare 411 le score global obtenu à un ensemble de plages prédéterminées correspondant à des matériaux respectifs. D'après la 15 plage dans laquelle se trouve le score global pour l'image analysée, on peut identifier le matériau correspondant. De préférence, les plages de scores correspondant aux matériaux sont élaborées préalablement lors d'une étape 500 représentée sur la figure 1 b. Pour ce faire, des étapes analogues à celles décrites ci-avant d'éclairage 100', d'acquisition 20 d'une image de tavelures 200', de traitement de l'image 300', et de calcul d'un score global 410', sont mises en oeuvre sur un grand nombre de surfaces, chaque surface étant formée en un matériau déterminé. Par exemple, des surfaces réalisées dans des matériaux tels que de l'épiderme, de la peau, de la cire d'abeille, de la colle, de la gélatine, du latex, de la pâte à modeler, de la porcelaine, de la silicone, etc., sont 25 analysées en vue d'établir des plages de score correspondantes. Les scores obtenues pour chaque surface sont ensuite agrégés au cours de l'étape 500 par un algorithme d'apprentissage, par exemple de type SVM (séparateurs à vaste marge, ou en anglais Support Vector Machine), analyse en composantes principales, ou autres, pour obtenir les plages de scores 30 correspondant à chaque matériau. Avantageusement, l'algorithme d'apprentissage élabore un polynôme permettant de séparer au mieux les matériaux testés. Ainsi l'évaluation ultérieure du même polynôme à partir des valeurs des critères calculées pour un matériau candidat permet de déterminer avec certitude la nature du matériau candidat. 3040815 15 En référence à la figure 4, on a représenté des exemples de scores globaux représentés pour un ensemble de matériaux, les scores étant établis à partir d'un critère calculé à partir de la transformée en ondelettes, un critère calculé à partir de la transformée en cosinus discrète de la fonction d'autocorrélation, et le critère C. 5 On constate que l'utilisation cumulée de plusieurs critères dont un critère calculé à partir de la transformée en ondelettes permet non seulement de distinguer l'épiderme de matériaux frauduleux, mais encore de distinguer les matériaux frauduleux entre eux. De retour à la figure 1, lorsque le procédé décrit ci-avant est mis en oeuvre 10 dans le cadre d'un contrôle biométrique, il comprend avantageusement une étape 600 de détection de fraude et d'émission d'une alerte, dans le cas où le matériau identifié n'est pas un tissu biologique, notamment pas de l'épiderme. Ceci implique que l'unité de traitement 24 commande l'activation de l'alarme 26, par exemple par émission d'un signal visuel (signal lumineux, affichage d'un 15 message d'alerte) ou sonore, et commande le verrouillage (ou le maintien à l'état fermé) du dispositif de restriction d'accès 25. Si au contraire le matériau est identifié comme un tissu biologique dont est constitué le trait biométrique sur lequel l'identification ou l'authentification est basée, le procédé peut continuer classiquement en procédé d'identification ou 20 d'authentification mis en oeuvre par le système 2 par comparaison d'une donnée biométrique acquise sur l'individu avec une donnée de référence. REVENDICATIONS1. Procédé de caractérisation d'un matériau dans lequel est formée une surface, comprenant les étapes consistant à : - projeter (100) sur la surface un faisceau lumineux cohérent, pour générer sur ladite surface des tavelures résultant d'interférences des rayons lumineux diffusés par ladite surface, - acquérir (200) une image de ladite surface sur laquelle apparaissent les tavelures, - traiter (300) ladite image pour calculer au moins un critère de caractérisation du matériau (C1, ..., C6), et - à partir des critères, déterminer (400) le matériau constituant la surface, dans lequel l'étape de traitement (300) de l'image comprend : - le calcul (310) d'une fonction d'autocorrélation de l'intensité lumineuse sur l'image et le calcul de la valeur d'au moins un critère (C1, ..., C4) établi à partir de ladite fonction, et - le calcul de la valeur d'au moins un critère représentatif (C5, C6) d'une distribution de déphasages des rayons lumineux diffusés par la surface. 2. Procédé de caractérisation selon la revendication 1, dans lequel le calcul (C5, C6) de la valeur d'au moins un critère représentatif d'une distribution de déphasages des rayons lumineux diffusés comprend le calcul (420) d'une transformée en ondelettes de l'image. 3. Procédé de caractérisation selon la revendication 2, comprenant le calcul de la valeur d'au moins un critère (C5, C6) parmi le groupe suivant : - moyenne (C6) des phases des coefficients d'ondelettes sur l'image, ou - écart-type (C5) des phases des coefficients d'ondelettes sur l'image. 4. Procédé de caractérisation selon l'une des revendications précédentes, dans lequel le calcul de la valeur d'au moins un critère de caractérisation à partir de la fonction d'autocorrélation de l'intensité lumineuse comprend : - le calcul de la valeur d'au moins un critère (C3, C4) établi à partir de la transformée de Fourier de la fonction d'autocorrélation, et 3040815 17 - le calcul de la valeur d'au moins un critère (C2) établi à partir de la transformée en cosinus discrète de la fonction d'autocorrélation. 5. Procédé de caractérisation selon la revendication 4, dans lequel le critère établi à 5 partir de la transformée de Fourier de la fonction d'autocorrélation comprend au moins l'un parmi le groupe suivant : - moyenne (04) de la transformée de Fourier de la fonction d'autocorrélation sur l'image, ou - écart-type (Ce) de la transformée de Fourier de la fonction d'autocorrélation 10 sur l'image. 6. Procédé de caractérisation selon l'une des revendications 4 ou 5, dans lequel le traitement (400) de l'image comprend : - l'établissement d'une matrice d'autocorrélation de l'intensité lumineuse de 15 l'image, et - le calcul de la transformée en cosinus discrète (312) de ladite matrice, et - le calcul de la valeur d'un critère (C2) de caractérisation comme la moyenne sur l'image de la trace de la matrice obtenue par la transformée en cosinus discrète de la matrice d'autocorrélation. 20 7. Procédé de caractérisation selon l'une des revendications précédentes, dans lequel le traitement de l'image comprend l'établissement d'une matrice d'autocorrélation de l'intensité lumineuse de l'image, et le calcul de la valeur d'un critère de caractérisation (Ci) égal à la racine carrée de la somme des valeurs 25 propres moyennes de ladite matrice sur l'image. 8. Procédé de caractérisation selon l'une des revendications précédentes, dans lequel l'étape de détermination (400) du matériau constituant la surface comprend le calcul (411) d'un score global à partir des valeurs des critères, et la classification du 30 matériau par comparaison (412) du score global avec une pluralité de plages prédéterminées correspondant à des matériaux respectifs. 3040815 18 9. Procédé de caractérisation selon la revendication 8, comprenant en outre une étape préliminaire (500) d'élaboration des plages de scores, ladite étape comprenant : - pour une pluralité de surfaces formées dans des matériaux déterminés, 5 mettre en oeuvre des étapes : o de projection sur la surface d'un faisceau lumineux cohérent, pour générer sur ladite surface des tavelures résultant d'interférences des rayons lumineux diffusés par ladite surface, o d'acquisition d'une image de ladite surface sur laquelle apparaissent 10 les tavelures, o de traitement de ladite image pour calculer au moins un critère de caractérisation du matériau, ledit traitement comprenant le calcul d'une fonction d'autocorrélation de l'intensité lumineuse sur l'image et le calcul de la valeur d'au moins un critère à partir de ladite fonction, 15 le calcul de la valeur d'au moins un critère représentatif d'une distribution de déphasages des rayons lumineux diffusés par la surface, et l'établissement d'un score global à partir desdites valeurs, - mettre en oeuvre un mécanisme d'apprentissage sur l'ensemble des scores globaux et l'ensemble des matériaux correspondants pour obtenir les plages 20 correspondant à chaque matériau. 10. Procédé de contrôle biométrique, comprenant la mise en oeuvre du procédé selon l'une des revendications précédentes sur une surface présentée par un individu lors du contrôle biométrique, le procédé comprenant en outre l'émission 25 d'une alerte dans le cas où le matériau caractérisé ne correspond pas à un tissu biologique. 11. Système (1) de caractérisation d'un matériau comprenant : - une source de lumière (10) cohérente, 30 - un capteur d'image (11), et - une unité de traitement (12), comprenant un processeur (120) adapté pour : o traiter une image d'une surface sur laquelle apparaissent des tavelures, ledit traitement comprenant le calcul d'une fonction d'autocorrélation de l'intensité lumineuse sur l'image, 3040815 19 o calculer, à partir desdites images, la valeur d'au moins un critère de caractérisation du matériau formant la surface, comprenant au moins un critère calculé à partir de la fonction d'autocorrélation et au moins un critère représentatif d'une distribution de déphasages des rayons 5 lumineux diffusés par la surface, et o à partir des valeurs des critères, déterminer le matériau constituant la surface. 12. Produit programme d'ordinateur, comprenant des instructions de code pour la 10 mise en oeuvre d'un procédé comprenant le traitement (300) d'une image d'une surface pour calculer la valeur d'au moins un critère de caractérisation du matériau constitutif de la surface, ledit traitement comprenant : - le calcul (310) d'une fonction d'autocorrélation de l'intensité lumineuse sur l'image et le calcul de la valeur d'au moins un critère (C1, ..., C4) établi à 15 partir de ladite fonction, et - le calcul de la valeur d'au moins un critère représentatif (C5, C6) d'une distribution de déphasages des rayons lumineux diffusés par la surface, lorsqu'il est exécuté par un processeur (120).