FR 3138691 A1 20240209 FR 2208107 A 20220804 L’invention concerne un procédé d’analyse d’une surface réfléchissante (S) par déflectométrie, avec un système comprenant un afficheur (Disp) et une optique d’imagerie (IO), le procédé comprenant les étapes consistant à : A afficher successivement sur l’afficheur (Disp) n images lumineuses Ili, une image lumineuse Ili présentant une intensité lumineuse binarisée obtenue par tramage du motif M(x-xi), B détecter sur le détecteur matriciel (Det) n images Idi obtenues par réflexion des n images lumineuses Ili sur la surface réfléchissante (S), C déterminer une fonction dite fonction de phase absolue objet (fap S ) à partir des n images détectées Idi, la fonction g vérifiant en outre la relation : g(k.M) = g(M) avec k réel quelconque, D comparer la fonction de phase absolue objet (fap S ) à une fonction de phase absolue (fap Mref ) dite de référence, et en déduire des informations sur la forme de ladite surface. Figure à publier : figure 3 Procédé de déflectométrie amélioré et système associé DOMAINE DE L’INVENTION La présente invention concerne le domaine de la déflectométrie pour l’analyse de la forme d’une surface. ETAT DE LA TECHNIQUE On connait le principe de la déflectométrie à décalage de phase ou PSD (pour « Phase Shifted Deflectometry » en anglais) pour analyser une surface S d’un objet Obj. Un montage classique d’un dispositif PSD 5 est illustré . Il comprend un afficheur Disp0 configuré pour afficher successivement 4 images Il0i(x,y), i indice variant de 1 à 4, les 4 images étant un motif sinusoïdal en intensité MS(x) selon un axe X respectivement décalé de 0, π/2, π, 3π/2. La illustre en A le motif sinusoïdal MS(x) et en B les quatre profils en intensité selon x des images/motifs MSi(x) avec MS1(x) correspondant à une phase de 0, MS2(x) correspondant à une phase de π/2, MS3(x) correspondant à une phase de π et MS4(x) correspondant à une phase de 3π/2. Sur l’exemple de la l’amplitude de 1 de la sinusoïde est codée en B entre 0 et 255 (8 bits). On définit le repère (X,Y) dans le plan de l’image affichée. En considérant que Il0i(x,y) représente l'intensité entre 0 et 1 (amplitude) d'une image et (x,y) les positions horizontale et verticale en pixel sur l'écran, on peut dire que par exemple MSi(x) = Il0i(x,0) car le motif est variable uniquement selon x. Ainsi, pour tout y, MSi(x) = Il0i(x,y) = 0.5.[cos((2π/Λ).x + φi) +1 ] où Λ est la période spatiale en pixel sur l'écran du sinus (ou cosinus) et φi est la phase de la ième image et prend pour valeur par exemple 0, π/2, π, 3π/2. L’afficheur Disp0 est disposé de sorte que les images lumineuses affichées Il0i se réfléchissent sur la surface S à caractériser (présentant localement une normale N) de l’objet Obj. On s’intéresse ici à la réflexion spéculaire et la réflexion sur la surface est oblique. Le dispositif comprend également une optique d’imagerie (non représentée) et un détecteur matriciel Det0, l’optique d’imagerie étant configurée pour imager la surface S sur le détecteur Det0. L’ensemble optique d’imagerie et détecteur est par exemple une caméra, qui prend séquentiellement 4 images Id0i correspondant à la partie du motif sinusoïdal MSi (affiché sur Disp0) réfléchie par S, imagée, et détectée par Det0. Les 4 images Id0i sont ensuite traitées par une unité de traitement (non représentée) qui détermine, à partir de ces 4 images détectées, les propriétés de la surface de l’objet Obj, et plus particulièrement son relief, afin typiquement d'identifier les défauts de surface. Dans une configuration de déflectométrie typique telle que celle représentée sur la , la caméra doit être maintenue inclinée par rapport à la partie de la surface à inspecter pour pouvoir imager la lumière réfléchie provenant de l’afficheur. L’afficheur et l’axe optique de la caméra doivent ainsi être positionnés au-delà d’un angle minimum par rapport à la partie de la surface à inspecter. En effet sans ajouter d'angle la camera bloquerait la lumière de l'écran ou l'inverse. Pour récupérer toute la lumière il suffit d'atteindre l'angle minimum à partir duquel la visibilité est totale. L’angle du rayon réfléchi par la surface S dépend de la pente locale de la surface au point d’impact du rayon. On dénomme (x,y) les coordonnées des images affichées Il0i(x,y) (repère X,Y), et (x’,y’) les coordonnées des images détectées Id0i(x’,y’), il existe typiquement un facteur de grandissement entre les deux systèmes. L’idée est que lorsque la surface S présente localement en un point P (x p ,y p ) une pente locale α P (α xP , α yP ) par rapport à la surface de référence plane, une frange réfléchie par cette surface S va être décalée sur le détecteur d’un décalage dx’ P dans la direction x et dy’ P dans la direction y, ces deux quantités dépendant directement et respectivement des pentes locales (α xP , α yP ). Ainsi en connaissant les paramètres géométriques du système, à partir du décalage dP(dx’ P , dy’ P ) on peut déduire la pente local α P (α xP , α yP ) sur la surface S en P. A partir de ces pentes α P on obtient par dérivation la courbure de S en chaque point et par intégration l’altitude de S en chaque point (distribution des hauteurs). Ainsi, dans le traitement on compare les 4 images détectées Id0i avec 4 images Id0 ref i obtenues par réflexion des images Il0i sur une surface de référence plane Sref positionnée à la place de l’objet Obj. Les images Id0 ref i sont des motifs sinusoïdaux non déformés semblables aux images affichées, au facteur de grandissement près. On définit un ensemble dit objet constitué des 4 images Id0i et un ensemble dit référence constitué des 4 images de références Id0 ref i. Pour la comparaison, on calcule une fonction de phase dite absolue pour chacun des deux ensembles d’images objet et référence, respectivement dénommées fonction de phase absolue objet fap0 S et fonction de phase absolue référence fap0 Ref , en appliquant une fonction mathématique g0 sur les 4 images correspondant à 4 variables que nous dénommerons pour la formule f1, f2, f3 et f4, égale à : g0 = arctg((f4-f2)/(f1-f3)) (1) Ainsi : fap0 S (x',y') = g0(Id01(x',y'), Id02(x',y'), Id03(x',y'), Id04(x',y')) = arctan[(Id04(x',y') –Id02(x',y')) / (Id01(x',y')-Id03(x',y'))]. fap0 Ref (x',y') = g0(Id0 ref 1(x',y'), Id0 ref 2(x',y'), Id0 ref 3(x',y'), Id0 ref 4(x',y') ) = arctan[(Id0 ref 4(x',y') –Id0 ref 2(x',y')) / (Id0 ref 1(x',y')-Id0 ref 3(x',y'))]. A titre illustratif, la fonction de phase absolue référence calculée avec les 4 motifs sinusoïdaux décalés de la est également illustrée . La fonction de phase absolue de référence calculée à partir des 4 motifs sinusoïdaux décalés est illustrée (courbe 20). On constate que cette fonction est une rampe sur des intervalles de 2π et qu’elle est donc bijective sur un intervalle de 2π: à une valeur de x (ou x’) dans cet intervalle correspond une valeur unique de y (ou y’): l’utilisation de 4 images avec un shift de phase différent permet de créer en chaque point une fonction de phase bijective. On calcule également de la même manière la fonction de phase absolue objet fap0 S à partir des 4 images déformées détectées Id0i(x’,y’). On obtient ainsi un « mapping » de la surface en phase absolue objet et en phase absolue référence. A partir de fap0 S et fap0 Ref , on reconstitue les trajets des différents rayons via un algorithme. Le caractère bijectif des fonctions de phase absolue référence et objet permet de relier une frange droite des images de références avec la frange déformée correspondante des images détectées, et ainsi de retrouver le décalage dP (mesuré en nombre de pixels dans le repère capteur puis transformé en distance) entre ces deux franges en un point P de la surface S. Le caractère bijectif est nécessaire pour créer un « mapping » parfait pour bien identifier les franges entres elles. A partir de ce décalage, on lie ce déplacement/décalage dP déterminé dans le repère du capteur x’ et y’ à des angles vus par le système d'imagerie optique couplé au capteur. Pour ça on utilise des formules de conjugaison et/ou d'optique géométrique. Une fois qu'on sait de quelle direction provient chaque phase absolue, on sait que chaque rayon frappe la surface S en un point, est réfléchi par les lois de Snell-Descartes jusqu'à un point de l'écran. On détermine ainsi, via les lois de Snell-Descartes, la pente de la surface α P au point P de la surface S. La méthode est ensuite réitérée dans la direction Y (le motif des quatre images affichées est MSi(y)) pour avoir les informations des pentes dans les deux directions X et Y. A partir de l’information sur les pentes selon X et Y, on en déduit la courbure et/ou l’altitude comme expliqué ci-dessus. La méthode PSD traite donc en tout au moins 8 images, un ensemble objet de 4 images détectées via un décalage du motif selon X et un ensemble objet de 4 images détectées via un décalage du motif selon Y. Il a été démontré que 4 phases telles que décrites précédemment par direction sont nécessaires et suffisantes pour l’obtention de la fonction de phase absolue. La méthode du PSD classique telle que décrite ci-dessus présente des inconvénients. Pour utiliser au mieux la dynamique de la caméra, on utilise plusieurs franges du motif affiché. Cela induit dans la fonction de phase absolue des sauts de phase correspondant à des discontinuités de la phase sur plusieurs périodes (visibles sur la rampe 20 de la ). Pour les traiter il convient de « dérouler la phase » c’est-à-dire de mettre en œuvre des algorithmes complexes en deux dimensions qui recherchent les sauts de phase réels par rapport aux erreurs de calculs. Dans le cas de plusieurs périodes, la fonction artan de g0 ne donne pas la phase absolue, il faut dérouler la phase pour l'obtenir. La mise en œuvre de ce « déroulement de la phase » allonge le temps d’analyse et la complexité des algorithmes associés, et génère des erreurs aux discontinuités. Comme la mise au point de la caméra est réglée sur la surface à caractériser, et non pas sur l’image virtuelle de l’image affichée, cela induit des problèmes de contraste (flou de défocalisation) des franges détectées. Pour palier cela dans les PSD classiques on ferme l’optique d’imagerie (ce qui augmente la dimension de la tâche d’Airy) et on augmente les temps d’intégration du détecteur, ce qui limite la vitesse d’analyse. Classiquement, les systèmes PSD présentent une configuration en incidence oblique. Il est alors difficile d’être au point sur l’ensemble de la surface. Lorsque l’on cherche à analyser des surfaces très plates (type « wafer » par exemple) ne pas être au point sur l’ensemble de la surface est une limitation. En effet si la profondeur de champs diminue le fait que l’on soit en incidence oblique implique qu'une portion décroissante de la surface à mesurer se retrouve dans la zone de netteté, et on perd alors des détails sur les bords. En outre en incidence oblique on doit fermer le système optique pour gagner en profondeur de champ. En plus d'augmenter la profondeur de champ, fermer le système optique diminue le flou de profondeur et cela permet de regagner du contraste lorsqu'il y a plusieurs franges à distinguer à l'écran. Le coût est un temps d'exposition plus long puisque le système est fermé et que moins de lumière est collectée. Enfin, les pixels des écrans utilisés pour la mise en œuvre de la PSD n’émettent pas une intensité transmise de façon linéaire par rapport à l’information initiale. Du fait de cette non linéarité l‘intensité émise par l’écran ne présente pas exactement une variation sinusoïdale : ce n’est pas un vrai sinus qui est généré. Cela induit au moment du calcul de la phase absolue l’apparition d’un bruit périodique dit de non linéarité que l’on n’arrive pas à supprimer facilement. Pour tenter de l’atténuer il convient : - de réaliser un calibrage sophistiqué du système ou - de caractériser finement la non linéarité de l’écran pour envoyer au pixel un signal prédéformé pour tenter de corriger la non linéarité et d’émettre un « vrai » sinus. Une telle démarche est par exemple décrite dans la publication «Non linearity response correction in phase shifting deflectometry » de Nguyen et al (2018 Meas. Sci. Technol. 29,045012). Un but de la présente invention est de remédier aux inconvénients précités en proposant un procédé et un système d’analyse d’une surface réfléchissante par déflectométrie amélioré, qui ne présente pas certains des inconvénients précités et présente en outre de meilleures performances qu’un procédé de déflectométrie à décalage de phase classique. DESCRIPTION DE L’INVENTION La présente invention a pour objet un procédé d’analyse d’une surface réfléchissante par déflectométrie, avec un système comprenant un afficheur disposé de sorte que des images lumineuses affichées sur l’afficheur se réfléchissent sur ladite surface réfléchissante et une optique d’imagerie configurée pour imager ladite surface sur un détecteur matriciel, le procédé comprenant les étapes consistant à : A. afficher successivement sur l’afficheur n images lumineuses Ili indicée i avec n≥2, les n images lumineuses étant obtenues par n décalages spatiaux xi selon une direction X d’un même motif, une image lumineuse Ili présentant une intensité lumineuse binarisée obtenue par tramage du motif M(x-xi), B. détecter sur le détecteur matriciel n images Idi dénommées images détectées, obtenues par réflexion des n images lumineuses Ili sur la surface réfléchissante, C. déterminer une fonction dite fonction de phase absolue objet à partir des n images détectées Idi, en appliquant une fonction mathématique g sur lesdites n images détectées, la fonction g et le nombre d’images n étant déterminés de sorte que la fonction de phase absolue soit bijective, la fonction g vérifiant en outre la relation : g(k.M) = g(M) avec k réel quelconque, D. comparer la fonction de phase absolue objet à une fonction de phase absolue dite de référence, déterminée en remplaçant la surface à caractériser par un miroir plan définissant un plan dit de référence, et en déduire des informations sur la forme de ladite surface. Selon une variante, le motif est une sinusoïde, n=4, les quatre décalages spatiaux correspondent à des décalages de phase respectivement égal à 0, π/2, π, 3π/2. Selon un mode de réalisation la fonction g est définie par : g= arctg[(fc4-fc2)/(fc1-fc3)] avec fci correspondant à n=4 variables. Selon un mode de réalisation, chaque image lumineuse n’affiche qu’une période du motif sinusoïdal. Selon une autre variante, le motif est une gaussienne, n=2, un premier décalage spatial correspondant à une gaussienne centrée sur un bord de l’image lumineuse Il1, et un deuxième décalage spatial correspondant à une gaussienne centrée sur un autre bord de l’image Il2, et dans lequel la fonction g est définie par : g=Ln(fc1) – Ln(fc2) avec fci correspondant à n=2 variables. Selon un mode de réalisation, le motif est invariant par transformée de Fourrier. Selon un mode de réalisation, le tramage est réalisé par l’algorithme de Floyd-Steinberg. Selon un mode de réalisation, l’optique d’imagerie présente un axe optique confondu avec une normale au plan de référence, dans lequel lesdites images lumineuses éclairent ladite surface perpendiculairement audit plan de référence, dans lequel on définit : - un centre optique O de l’optique d’imagerie, - une distance Ddp positive entre l’afficheur et le plan de référence Pref et une distance Dpo positive entre le plan de référence et le centre optique O, et dans lequel l’étape D comprend les sous étapes consistant à : - déterminer pour des points de la surface S, un angle β déterminé à partir du décalage entre la fonction phase absolue objet et la fonction phase absolue référence audit point de la surface, - déterminer ladite pente α avec la formule suivante : Selon un autre aspect, l’invention concerne un système d’analyse d’une surface réfléchissante par déflectométrie comprenant : - un afficheur disposé de sorte que des images lumineuses affichées sur l’afficheur se réfléchissent sur ladite surface réfléchissante, et configuré pour afficher successivement n images lumineuses Ili indicée i avec n≥2, les n images lumineuses étant obtenues par n décalages spatiaux xi selon une direction X d’un même motif M(x), une image lumineuse Ili(x,y) présentant une intensité lumineuse binarisée obtenue par tramage du motif M(x-xi), - une optique d’imagerie et un détecteur matriciel, l’optique d’imagerie étant configurée pour imager ladite surface sur ledit détecteur matriciel, le détecteur matriciel étant configuré pour détecter n images Idi dénommées images détectées, obtenues par réflexion des n images lumineuses Ili sur la surface réfléchissante, - une unité de traitement configurée pour : - déterminer une fonction dite fonction de phase absolue objet, à partir des n images détectées Idi, en appliquant une fonction mathématique g sur lesdites n images détectées, la fonction g et le nombre d’images n étant déterminés de sorte que la fonction de phase absolue soit bijective, la fonction g vérifiant en outre la relation : g(k.M) = g(M) avec k réel quelconque, - comparer la fonction de phase absolue objet à une fonction de phase absolue dite de référence, déterminée en remplaçant la surface à caractériser par un miroir plan définissant un plan dit de référence, et en déduire des informations sur la forme de ladite surface. Selon un mode de réalisation, l’optique d’imagerie présente un axe optique disposé perpendiculairement audit plan de référence, et dans lequel le système optique comprend en outre une lame séparatrice configurée pour envoyer lesdites images lumineuse sur ladite surface perpendiculairement audit plan de référence. Selon un autre aspect, l’invention concerne un programme d'ordinateur comportant des instructions qui conduisent le système selon l’invention à exécuter les étapes du procédé selon l’invention. La description suivante présente plusieurs exemples de réalisation du dispositif de l’invention : ces exemples sont non limitatifs de la portée de l’invention. Ces exemples de réalisation présentent à la fois les caractéristiques essentielles de l’invention ainsi que des caractéristiques additionnelles liées aux modes de réalisation considérés. L’invention sera mieux comprise et d’autres caractéristiques, buts et avantages de celle-ci apparaîtront au cours de la description détaillée qui va suivre et en regard des dessins annexés donnés à titre d’exemples non limitatifs et sur lesquels : La déjà citée illustre un montage classique d’un dispositif de déflectométrie à décalage de phase. La déjà citée illustre en A le motif sinusoïdal MS(x) et en B les quatre motifs MSi(x) décalés et la fonction de phase absolue associée. La illustre le procédé selon l’invention. La illustre le système selon l’invention. La illustre un exemple d’image binarisée par tramage d’un motif sinusoïdal affichée sur l’afficheur dans le procédé/système selon l’invention. La illustre en A la répartition d’intensité de l’image affichée dans trois plans successifs lors de la propagation de la lumière dans l’espace : I : sur l’afficheur (d0=0); II : à une distance d1 de Disp pour laquelle le moyennage n’est pas encore suffisant ; III : à une distance d2 suffisamment éloignée de l’écran Disp pour que le moyennage soit abouti. En B est illustré, pour les trois situations, un exemple de la variation de l’intensité I de l’image selon l’axe X. La schématise un exemple de la MTF de l’image détectée: en A dans la situation I et en B dans la situation III. La illustre un mode de réalisation de l’invention dans lequel chaque image lumineuse n’affiche qu’une période du motif sinusoïdal et la fonction de phase absolue référence associée. La illustre un exemple de mesure réalisée avec un système selon l’invention affichant une seule période du motif sinusoïdale binarisé pour la mesure d’une surface S concave. La concavité a été déduite pour faire apparaître les défauts de la surface. La illustre une mesure de la même surface concave (concavité déduite) effectuée avec un système PSD classique La illustre la mesure d’un wafer 2 pouces en silicium effectuée avec le système selon l’invention. La illustre la mesure de ce même wafer avec un appareil du commerce de type interféromètre laser (Zygo). La illustre une deuxième variante de l’invention dans laquelle le motif est une gaussienne et la fonction de phase absolue référence associée. On a n=2, l’image Il1 est constituée du motif M1(x) fonction gaussienne centrée sur un bord de l’image, et l’image Il2 est constituée du motif M2(x) fonction gaussienne centrée sur le bord opposé de l’image. La illustre le décalage dP correspondant à une pente α en un point C de la surface à caractériser. La est un schéma optique illustrant différentes grandeurs utiles à la compréhension de l’invention DESCRIPTION DETAILLEE DE L’INVENTION L’objet de la présente invention est un procédé 100 d’analyse d’une surface réfléchissante S par déflectométrie, illustré , qui constitue une amélioration du procédé PSD classique tel que décrit précédemment. Selon un autre aspect l’invention concerne également un système 10 illustré qui met en œuvre le procédé 100. Le système 10 comprend un afficheur Disp disposé de sorte que des images lumineuses affichées sur l’afficheur se réfléchissent sur la surface réfléchissante S, et une optique d’imagerie IO configurée pour imager la surface sur un détecteur matriciel Det. L’ensemble optique d’imagerie/détecteur forme par exemple une caméra CA. Le procédé selon l’invention est une généralisation du PSD classique dans le sens qu’il ne fonctionne pas uniquement avec un motif M(x) sinusoïdal mais également avec d’autres motifs, et permet de résoudre certains problèmes inhérents au PSD, ce qui en améliore les performances. Le procédé comprend une première étape A dans laquelle on affiche successivement sur l’afficheur Disp n images lumineuses Ili indicée i avec n≥2, les n images lumineuses étant obtenues par n décalages spatiaux xi selon une direction X d’un même motif M(x). On définit X et Y les deux directions du repère de l’image Ili(x,y). A la différence d’une image lumineuse utilisée dans un PSD classique, les images lumineuses Ili présentent une intensité lumineuse binarisée obtenue par tramage du motif M(x-xi) également dénommé Mi(x). Après de nombreuses expériences et simulations les inventeurs ont établis qu’un tel tramage (« dithetering » en anglais) était de manière surprenante compatible d’une caractérisation de la surface par déflectométrie à décalage de phase, en améliorait très nettement les performances tout en permettant l’utilisation d’autres motifs qu’une sinusoïde (voir plus loin). On entend par « image binarisée obtenue par tramage » la transformation d’une image à niveaux de gris, codés typiquement entre 1 et 255 (8 bits), en une image ne comprenant que des points présentant une intensité non nulle, préférentiellement d’intensité maximale (codée 1 ou 255) et des points présentant une intensité nulle (codés 0). Le tramage génère ainsi une image dont la densité spatiale de points varie selon le motif M(x) le long de la direction X. Il existe de nombreux algorithmes de tramage et les inventeurs en ont testés plusieurs qui donnent de bons résultats. On cite à titre non limitatif : - l’algorithme consistant à générer une image aléatoire dans laquelle chaque pixel a une probabilité de la forme (sin(x)+1)/2 d’être à 1 (ou 0). - l’algorithme de Floyd-Steinberg, D’autres algorithmes de « dithering » sont accessibles aisément sur internet. La illustre un exemple d’image binarisée par tramage d’un motif sinusoïdal de période T. Dans une deuxième étape B on détecte sur le détecteur matriciel Det n images Idi(x’,y’) dénommées images détectées, obtenues par réflexion des n images lumineuses Ili sur la surface réfléchissante S. Les coordonnées (x’, y’) de l’image détectées sont reliées aux coordonnées (x, y) de l’image affichée par un facteur de grandissement. On dénomme toujours (X, Y) les directions du repère de l’image détectée. Cette étape B est identique à ce qui se pratique classiquement en PSD, sauf que classiquement on a toujours n=4 et M(x) sinusoïde. Dans une étape C, on détermine une fonction dite fonction de phase absolue objet fap S à partir des n images détectées Idi, en appliquant une fonction mathématique g sur les n images détectées, la fonction g et le nombre d’images n étant déterminés de sorte que la fonction de phase absolue soit bijective, c’est-à-dire: y’ = fap S (x’) avec y’ unique pour chaque x’ (2) Lorsque le motif est sinusoïdal et n=4 la fonction mathématique g est égale à g0 donné par la formule (1). La fonction g vérifie en outre la relation : g(k.M) = g(M) (3) avec k réel quelconque. Cette condition (3) provient du fait que l’opérateur g doit être insensible à la réflectivité de la surface, l’intensité réfléchie par la surface S et détectée par le détecteur Det étant proportionnelle à cette réflectivité. La fonction g appliquée aux images ne doit donc pas garder cette réflectivité qui doit s’annuler. Un exemple est une annulation de la réflectivité par le biais d’un ratio entre les images (voir formule 1). Lorsque qu’un motif autre qu’une sinusoïde est utilisé, la fonction g est différente de la formule (1). Dans une étape D, on compare la fonction de phase absolue objet fap S à une fonction de phase absolue fap Mref dite de référence, déterminée en remplaçant la surface S à caractériser par un miroir plan Mref définissant un plan Pref dit de référence illustré et on en déduit des informations sur la forme de ladite surface, dans le même esprit que ce qui est réalisé en PSD classique. La fonction de phase absolue référence est obtenue à partir des 4 images détectées de références Idrefi. A noter que l’optique d’imagerie IO est configurée pour imager la surface S et non pas l’image virtuelle de l’afficheur (de même que pour une PSD classique). L’image affichée sur Disp n’est donc pas au point sur le détecteur Det, elle est floutée. Comme dans le PSD classique les étapes A à D sont réitérées selon l’axe Y pour obtenir une information de pente selon la direction Y. Lorsque l’on ne change pas les paramètres géométriques du système on peut utiliser une même fonction de phase absolue référence pour plusieurs mesures de différentes surfaces. Si l’on change les paramètres, par exemple la distance afficheur-surface, et/ou la distance surface-optique d’imagerie, et/ou la focale de l’optique d’imagerie, il convient de procéder à la détection des 4 images Id ref i (avec le miroir Mref à la place de S) et à la détermination de la fonction de phase absolue référence à partir de ces 4 images dans la nouvelle configuration. Le système 10 selon l’invention illustré comprend donc : - un afficheur Disp disposé de sorte que des images lumineuses affichées sur l’afficheur se réfléchissent sur la surface réfléchissante S, et configuré pour afficher successivement les n images lumineuses Ili binarisées, - une optique d’imagerie IO et un détecteur matriciel Det, l’optique d’imagerie étant configurée pour imager la surface S sur le détecteur Det, le détecteur matriciel étant configuré pour détecter les n images Idi obtenues par réflexion des n images lumineuses Ili sur la surface réfléchissante S, - une unité de traitement (PU) configurée pour mettre en œuvre les étapes C de détermination de la phase absolue objet et D de comparaison des fonctions de phase absolues objet et référence et de détermination d’informations sur la forme de S. De manière non limitative, la illustre une configuration en incidence normale qui est préférée pour la mise en œuvre du procédé selon l’invention. On entend par « incidence normale » une configuration dans laquelle l’axe optique OA de l’optique d’imagerie IO est perpendiculaire au plan de référence Pref et dans laquelle le système comprend une lame séparatrice BS configurée pour envoyer les images lumineuses sur la surface S perpendiculairement au plan de référence Pref. Cette configuration en incidence normale est peu utilisée en PSD classique. Comme on cherche à déterminer des angles de pente α très petits, l’ensemble de la surface S est au point sur le détecteur. Les inventeurs ont établi que la binarisation de l’image affichée Il, conjuguée à la défocalisation de l’image de Il réalisée par IO, permet d’obtenir sur la surface S une lumière issue de l’afficheur correspondant à « vrai » sinus, plutôt qu’un sinus dégradé comme c’est le cas en PSD classique. Cela est dû au fait que l’image binarisée se moyenne en se propageant dans l’espace et donne, si l’afficheur est disposé à une distance suffisante de la surface S, une fonction sinusoïdale proche d’un sinus parfait. La illustre en A la répartition d’intensité de l’image affichée Il dans trois plans successifs lors de la propagation de la lumière dans l’espace : I : sur l’afficheur (d0=0); II : à une distance d1 de Disp pour laquelle le moyennage n’est pas encore suffisant ; III : à une distance d2 suffisamment éloignée de l’écran Disp pour que le moyennage soit abouti. En B est illustré, pour les trois situations, un exemple de la variation de l’intensité I de l’image selon l’axe X (qui oscille entre 1 (255) et 0) pour une coupe transversale selon un axe 50, faisant apparaître le moyennage au fil de la propagation. Le flou gaussien de l’image binarisée correspond à la propagation de la lumière en espace libre. En d’autres termes, l’image en I correspond à ce que le détecteur détecte lorsque la caméra est au point sur l’image virtuelle Ivir (l’image en I est ici une image réelle de l’image virtuelle de l’écran Ivir créée par la surface S), et l’image en III correspond à ce que le détecteur détecte lorsque la caméra est au point sur la surface S, soit lorsque la mesure est effectuée. La schématise un exemple de la MTF (pour « Modulation Transfert Function » en anglais ou fonction de transfert de modulation, fonction de la fréquence spatiale fs en cycle/mm) de l’image détectée par Det : en A dans la situation I (focus sur l’image virtuelle de l’afficheur) pour laquelle l’image présente des hautes fréquences spatiales, et en B dans la situation III (focus sur la surface S) pour laquelle les hautes fréquences spatiale ont été éliminées : les fréquences spatiales élevées présentes dans l’image binarisée sont éliminées lors de la propagation et un « vrai » sinus est créé sur la surface S. Ainsi, la binarisation de l’image affichée permet de diminuer considérablement voir de supprimer le bruit de phase périodique de non linéarité présent dans le PSD classique, qui n’existe plus ici puisque l’image affichée est binarisée et non plus analogique. Cela est vrai car les écrans LCD sont bien homogènes pour une intensité donnée, mais sont non linéaires pour des variations d’intensité. Cette suppression du bruit évite des calibrations complexes du système et en simplifie grandement la mise en œuvre. Pour la mesure le réglage de la distance afficheur-surface doit être suffisante pour que le moyennage par propagation soit effectif, c’est à dire pour que s’opère le filtrage des hautes fréquences de l’image binarisée. Ce réglage s’effectue par exemple en vérifiant visuellement que le motif détecté présente une variation considérée comme satisfaisante ou en caractérisant la courbe de variation de l’intensité détectée selon un axe parallèle à X (voir courbes de la partie B pour II et IIII). Pour le réglage de l'ouverture du système, on ouvre au maximum pour diminuer la profondeur de champ et maximiser le moyennage (pas trop cependant si l'optique utilisée est moins performante aux ouvertures maximales). En résumé, les variables d'ajustement pour bien moyenner le « dither » sont : - prendre un écran avec une plus grande densité de pixel, - ouvrir le système optique, - éloigner l'écran, - éventuellement rapprocher l'objectif de la cible, - de manière optionnelle coller un dépoli sur l'écran afin d'exacerber le moyennage en plus de la propagation. Le système 10 selon l’invention avec une configuration en incidence normale permet de rapprocher l’optique d’imagerie IO de la surface à caractériser (voir plus loin). Cela permet d’augmenter la valeur relative de l’écart entre frange droite et frange courbe précité, directement relié à la valeur de la pente α recherchée, et donc de mesurer des valeurs de α plus petites qu’avec le PSD classique. La configuration en incidence normale permet également une mise au point sur toute la surface, ce qui est très utile pour la caractérisation de surface très plane, soit avec des pentes α très faibles. A une distance suffisante de l’afficheur la propagation de l’onde émise s’apparente à une transformée de Fourier. Préférentiellement on choisit un motif M(x) invariant par transformée de Fourier, pour obtenir sur la surface S un motif de même nature que celui affiché sur Disp. Comme expliqué plus haut, selon une première variante, le motif M(x) est un motif sinusoïdal (binarisé). Dans ce cas on a n=4 et les quatre décalages spatiaux xi correspondent à des décalages de phase respectivement égal à 0, π/2, π, 3π/2. Préférentiellement, la fonction mathématique g appliquée aux 4 images détectées Idi (objet) et aux 4 images référence Id ref i, correspondant à 4 variables dénommées fci, est : g= arctg[(fc4-fc2)/(fc1-fc3)]. avec fci = Idi (pour le calcul de fap S ) et fci = Id ref i (pour le calcul de fap ref ) Cette formule est identique à la formule (1) utilisée pour un procédé PSD classique. Selon un mode de réalisation du procédé 100 et du système 10, selon l’invention, chaque image lumineuse n’affiche qu’une période du motif, tel qu’illustré pour le motif sinusoïdal : on diminue la fréquence spatiale du sinus du motif. En effet avec l’invention on ouvre l’objectif, ce qui optimise les fréquences spatiales près de l’échantillon à caractériser, et diminue le contraste des franges affichées sur l’écran. Pour récupérer un bon contraste sur l’afficheur il convient donc de baisser la fréquence spatiale des franges affichées. La fonction de phase absolue objet calculée avec la formule (1) est égale à la courbe 80 qui ne présente plus de discontinuités (sauts de phase). N’utiliser qu’une seule période permet de ne pas avoir à « dérouler » la phase (voir état de la technique) ce qui rend l’algorithme pour remonter au dy et aux pentes beaucoup moins complexe et plus rapide. En outre avec une seule frange on a plus de signal car on ne perd pas de contraste sur les franges détectées. En effet si on utilise plusieurs franges, i.e une plus grande fréquence spatiale sur l'afficheur, comme dit plus haut le contraste des hautes fréquences spatiale diminuent. Une diminution du contraste implique une perte sur la dynamique de la camera. Avec une seule frange, en regardant la MTF on voit que les basses fréquences sont bien mieux transmises que les hautes fréquences, peut important la mise au point. On garde donc un bon contraste, et on a plus de signal car max-min tend vers 1 : on exploite alors au maximum la dynamique de la camera. La illustre un exemple de mesure réalisée avec un système selon l’invention affichant une seule période du motif sinusoïdale binarisé pour la mesure d’une surface S concave. La concavité a été déduite pour faire apparaître les défauts de la surface. On voit sur l’échelle de droite que la sensibilité de mesure sur les défauts résiduels de la surface est d’environ +/- 30 nm. On voit clairement les défauts de la surface. Le système présente les caractéristiques suivantes : - Distance centre optique O de l’optique d’imagerie- échantillon : 47 cm - Distance écran-Pref : 71 cm - Focale de l’optique d’imagerie : 50 mm - Résolution du détecteur : 4512x4512 pixels, taille d’un pixel 2,74 µm - Ouverture du dispositif d’imagerie : N=2 - Ecran utilisé : écran LCD d’ordinateur classique A titre de comparaison, la illustre une mesure de la même surface concave (concavité déduite) effectuée avec un système PSD classique affichant une période d’un motif analogique, sans précorrection de la non linéarité de l’afficheur et avec le même système que décrit ci-dessus. On voit sur l’échelle de droite que la sensibilité de mesure sur les défauts résiduels de la surface est d’environ +/- 150 nm. Cela est dû au bruit de phase qui vient ajouter une modulation parasite. Cette modulation correspond aux zones sombres et claires dans l’image. Du fait de ce signal parasite la précision de mesure est inférieure et les défauts sont moins visibles. La illustre une mesure effectuée avec le système selon l’invention d’une surface de type wafer 2 pouces en silicium (courbure globale déduite) Les sillons de polissage d'épaisseur nanométrique sont apparents et l’image de la surface est complète. L’amplitude de variation d’épaisseur mesuré est de +/- 100 nm soit 200 nm (voir échelle de droite). La illustre une mesure de ce même wafer (courbure globale déduite) effectuée avec un appareil du commerce de type interféromètre laser Zygo. On voit une zone 11 sans mesures correspondant à un artefact de mesure. L’objet est en effet trop courbé pour que la mesure soit possible avec le Zygo, celui-ci présentant une dynamique de mesure de +/-80 nm soit 160 nm maximum. Les zones correspondant à une variation d’épaisseur allant au-delà ne sont pas traitées et apparaissent en blanc. Le système selon l’invention est ainsi plus performant que le Zygo pour cette mesure : la dynamique est bien plus grande et les détails à haute fréquences spatiales (sillons) sont mieux définis. Les inventeurs ont montré que le procédé selon l’invention était compatible avec d’autres motifs que le motif sinusoïdal. Selon une deuxième variante, le motif est une gaussienne et n=2. L’image Il1 est constituée du motif M1(x) fonction gaussienne centrée sur un bord de l’image, et l’image Il2 est constituée du motif M2(x) fonction gaussienne centrée sur le bord opposé de l’image, tel qu’illustré . Le maximum de la gaussienne est codé à 1. Les deux gaussiennes ont les mêmes paramètres (amplitude, largeur). Les inventeurs ont montré qu’une fonction g définie par la formule (4) ci-dessous est bijective et vérifie la condition (3) : g=Ln(fc1) – Ln(fc2) (4) La fonction de phase absolue référence calculée avec la fonction g de la formule (4) est également illustrée (courbe 12). Elle ne présente pas de discontinuité. A noter que la fonction g de la formule (4) vérifie bien la condition g(kM) = g(M), et qu’une fonction gaussienne est invariante par transformée de Fourrier. Il convient de s’assurer que le maximum de la gaussienne ne sature pas la camera et que son minimum de la gaussienne n'est pas en dessous du niveau de bruit de la camera. L’avantage de ce motif gaussien est que 2x2 images sont nécessaires et suffisantes pour la mise en œuvre du procédé selon l’invention, à comparer aux 2x4 images avec l’utilisation d’un motif sinusoïdal. Cela permet une acquisition, un traitement, et donc une mesure plus rapide. Le système selon l’invention présente une grande souplesse car on peut modifier les paramètres géométriques en fonction de la précision recherchée, soit l’ordre de grandeur de l’angle α de la pente locale sur la surface S. La illustre le décalage dP correspondant à un angle α en un point C de S. Ce décalage dP mesuré dans le repère (x’, y’) du détecteur est la distance entre la phase absolue d'une frange lorsque α = 0 (référence) et la phase absolue de la même frange lorsque la surface C présente un angle α localement. Ce décalage se mesure séparément selon x’ (motif affiché selon x) et selon y’ (motif affiché selon y). On définit le centre optique O de l’optique d’imagerie et on dénomme β l’angle correspondant à ce décalage dP mesuré par rapport au centre O (O est le sommet de β). Ainsi, l’angle β est déterminé à partir du décalage dP entre la fonction phase absolue objet et la fonction phase absolue référence au point C de la surface. Les inventeurs ont déterminé une relation entre α et β en fonction des paramètres du système pour la configuration en incidence normale préférée. Cette relation ne fait intervenir que deux grandeurs : une distance Ddp positive entre l’afficheur et le plan de référence Pref, et une distance Dpo positive entre le plan de référence Pref et le centre optique O. Le schéma de la illustre ces différentes grandeurs dans un schéma optique pour un exemple numérique avec Ddp = 11 cm et Dpo = 31 cm. P Disp est le plan dans lequel est disposé l’afficheur Disp et P Det est le plan dans lequel est disposé le détecteur Det. Le point A’ est un point de l’afficheur. Le point A est l’image virtuelle de A’ lorsque le miroir de référence Mref est disposé en Pref (perpendiculairement à OA). Un rayon issu de A’ est réfléchit au point C (rayon 3) est focalisé au point F dans le plan P Det . Le point A’’ est l’image virtuelle de A’ lorsque c’est la surface S avec en C une pente locale α qui réfléchit le rayon issu de A’. Ce rayon réfléchit par C (rayon 4) est focalisé au point G dans le plan P Det . Cela est valable pour des petits angles. Le rayon est réfléchi par un point très proche du point C. Pour des petits angles, et pour établir la formule ci-dessous, on a considéré que ce point est C. La distance FG correspond au décalage dP. Ainsi, l’angle β est l'angle formé entre les points F, O (centre optique-sommet) et G où F est la position de la phase absolue sur le détecteur d'une frange de l'écran (A') après réflexion sur la référence et G est la position de la phase absolue sur le détecteur de la même frange après réflexion sur l'échantillon présentant un angle α. La relation établie par les inventeurs est : (5) Cette formule est bien entendu à appliquer selon X et Y indépendamment. La phase absolue d'une frange de référence tombe en F, la phase absolue de la même frange après réflexion sur l'échantillon tombe en G. Lorsque la phase varie dans la direction X, alors F dans le repère (x’, y’) est en x’ F et G en x’ G , et on a : dx' = x’ G -x’ F , distance entre F et G Pour le cas simplifié d’une lentille mince on a: βx = arctan(dx’/f) f focale de l’optique d’imagerie Pour un objectif plus complexe, on tient compte de la position des plans principaux. α X est alors déterminée avec la formule (5). Ainsi, selon un mode de réalisation du procédé 100, selon l’invention dans l’étape D, on détermine d’abord l’angle β puis l’angle α à l’aide de la formule (5). On voit grâce à cette formule (5) que la relation entre α et β peut être ajustée en fonction des valeurs des distances Dpo et Ddp choisies. Il s’agit d’un avantage important du système/procédé selon l’invention. Dans un système PSD classique en incidence oblique classiquement ces deux distances sont égales et on a donc α = β. Voir par exemple la publication « Review of phase measuring deflectometry » de Huang et al (Optics and Lasers Engineering, volume 107, august 2018, pages 247-257). Lorsqu’on souhaite mesurer des pentes α très faibles, on cherche à augmenter le facteur de proportionnalité entre α et β. Une étude de la formule (5) montre qu’augmenter ce facteur de proportionnalité revient à diminuer/minimiser le coefficient Dpo/Ddp. A focale constante cela revient à rapprocher la caméra (diminuer Dpo) et éloigner l'écran (augmenter Ddp). Une mesure des pentes α très faibles s’obtient également en augmentant la focale de l’optique d’imagerie. Le document «Phase measuring deflectometry for obtaining 3D shape of specular surface: a review of the state of the art » de Zhang et al (Optical Engineering, vol 60(2), 020903-1; 2021) décrit une méthode PSD classique avec ces deux distances Dpo et Ddp différentes. Mais il est alors nécessaire d’effectuer une calibration dans l’espace de la caméra très complexe pour obtenir in fine les bonnes valeurs (voir §4 Error source analysis). Procédé d’analyse d’une surface réfléchissante (S) par déflectométrie, avec un système comprenant un afficheur (Disp) disposé de sorte que des images lumineuses affichées sur l’afficheur se réfléchissent sur ladite surface réfléchissante (S) et une optique d’imagerie (IO) configurée pour imager ladite surface sur un détecteur matriciel (Det), le procédé comprenant les étapes consistant à : A afficher successivement sur l’afficheur (Disp) n images lumineuses Ili indicée i avec n≥2, les n images lumineuses étant obtenues par n décalages spatiaux xi selon une direction X d’un même motif M(x), une image lumineuse Ili présentant une intensité lumineuse binarisée obtenue par tramage du motif M(x-xi), B détecter sur le détecteur matriciel (Det) n images Idi dénommées images détectées, obtenues par réflexion des n images lumineuses Ili sur la surface réfléchissante (S), C déterminer une fonction dite fonction de phase absolue objet (fap S ) à partir des n images détectées Idi, en appliquant une fonction mathématique g sur lesdites n images détectées, la fonction g et le nombre d’images n étant déterminés de sorte que la fonction de phase absolue soit bijective, la fonction g vérifiant en outre la relation : g(k.M) = g(M) avec k réel quelconque, D comparer la fonction de phase absolue objet (fap S ) à une fonction de phase absolue (fap Mref ) dite de référence, déterminée en remplaçant la surface (S) à caractériser par un miroir plan (Mref) définissant un plan (Pref) dit de référence, et en déduire des informations sur la forme de ladite surface. Procédé selon la revendication précédente dans lequel le motif est une sinusoïde, n=4, les quatre décalages spatiaux correspondent à des décalages de phase respectivement égal à 0, π/2, π, 3π/2. Procédé selon la revendication précédente dans lequel la fonction g est définie par : g= arctg[(fc4-fc2)/(fc1-fc3)] avec fci correspondant à n=4 variables, fci correspondant aux images détectées avec ladite surface réfléchissante pour le calcul de la fonction de phase absolue objet, et aux images détectées avec ledit miroir de référence pour le calcul de la fonction de phase absolue référence. Procédé selon l’une des revendications 2 ou 3 dans lequel chaque image lumineuse n’affiche qu’une période du motif sinusoïdal. Procédé selon la revendication 1 dans lequel le motif est une gaussienne, n=2, un premier décalage spatial correspondant à une gaussienne centrée sur un bord de l’image lumineuse Il1 et un deuxième décalage spatial correspondant à une gaussienne centrée sur un autre bord de l’image Il2, et dans lequel la fonction g est définie par : g=Ln(fc1) – Ln(fc2) avec fci correspondant à n=2 variables, fci correspondant aux images détectées avec ladite surface réfléchissante pour le calcul de la fonction de phase absolue objet, et aux images détectées avec ledit miroir de référence pour le calcul de la fonction de phase absolue référence. Procédé selon l’une des revendications précédentes dans lequel le motif est invariant par transformée de Fourrier. Procédé selon l’une des revendications précédentes dans lequel le tramage est réalisé par l’algorithme de Floyd-Steinberg. Procédé selon l’une des revendications précédentes dans lequel l’optique d’imagerie présente un axe optique confondu avec une normale au plan de référence, dans lequel lesdites images lumineuses éclairent ladite surface perpendiculairement audit plan de référence, dans lequel on définit : un centre optique O de l’optique d’imagerie, une distance Ddp positive entre l’afficheur et le plan de référence Pref et une distance Dpo positive entre le plan de référence et le centre optique O, et dans lequel l’étape D comprend les sous étapes consistant à : déterminer pour des points de la surface S, un angle β déterminé à partir du décalage entre la fonction phase absolue objet et la fonction phase absolue référence audit point de la surface, déterminer ladite pente α avec la formule suivante : Système (10) d’analyse d’une surface réfléchissante (S) par déflectométrie comprenant : un afficheur (Disp) disposé de sorte que des images lumineuses affichées sur l’afficheur se réfléchissent sur ladite surface réfléchissante (S), et configuré pour afficher successivement n images lumineuses Ili indicée i avec n≥2, les n images lumineuses étant obtenues par n décalages spatiaux xi selon une direction X d’un même motif M(x), une image lumineuse Ili(x,y) présentant une intensité lumineuse binarisée obtenue par tramage du motif M(x-xi), une optique d’imagerie (IO) et un détecteur matriciel (Det), l’optique d’imagerie étant configurée pour imager ladite surface (S) sur ledit détecteur matriciel, le détecteur matriciel étant configuré pour détecter n images Idi dénommées images détectées, obtenues par réflexion des n images lumineuses Ili sur la surface réfléchissante (S), une unité de traitement (PU) configurée pour : déterminer une fonction dite fonction de phase absolue objet (fap S ), à partir des n images détectées Idi, en appliquant une fonction mathématique g sur lesdites n images détectées, la fonction g et le nombre d’images n étant déterminés de sorte que la fonction de phase absolue soit bijective, la fonction g vérifiant en outre la relation : g(k.M) = g(M) avec k réel quelconque, comparer la fonction de phase absolue objet (fap S ) à une fonction de phase absolue (fap Mref ) dite de référence, déterminée en remplaçant la surface (S) à caractériser par un miroir plan (Mref) définissant un plan (Pref) dit de référence, et en déduire des informations sur la forme de ladite surface. Système selon la revendication précédente dans lequel l’optique d’imagerie présente un axe optique (OA) disposé perpendiculairement audit plan de référence, et dans lequel le système optique comprend en outre une lame séparatrice (BS) configurée pour envoyer lesdites images lumineuse sur ladite surface perpendiculairement audit plan de référence. Programme d'ordinateur comportant des instructions qui conduisent le système de la revendication 9 à exécuter les étapes du procédé selon l’une des revendications 1 à 8.