FR 3074949 A1 20190614 FR 1761926 A 20171211 L'invention concerne un procédé de traitement d'images, obtenues par détecteur à diffraction. Le domaine d'application de l'invention concerne l'analyse de matériaux cristallins ou polycristallins, et notamment l'analyse par diffraction rétrodiffusée en haute résolution (en anglais : HR-EBSD pour High Resolution Electron BackScatter Diffraction). Notamment, le détecteur peut être un détecteur à diffraction électronique produisant des figures de diffraction électronique, appelées figures de Kikuchi, dont l'analyse permet de calculer très précisément les variations (relatives) des paramètres de maille cristalline. On cherche notamment à comparer des paramètres de maille cristalline d'un même matériau se trouvant dans des états de contraintes différents, ce qui affecte ces paramètres. L'ensemble de l'information obtenue, point par point, lors du balayage de la surface d'un échantillon est habituellement présenté sous la forme d'images bidimensionnelles pour différentes composantes de la cristallographie, voire de la qualité de l'analyse, ou du degré de confiance qu'on lui accorde. On connaît dans l'état de l'art de la technique, la comparaison de figures de diffraction entre deux états différents à l'aide d'algorithmes basés sur une inter- corrélation, par exemple par le logiciel dénommé Crosscourt et par le logiciel dénommé Straincorrelator, également basé sur le logiciel Crosscourt. Cette inter-corrélation repose sur un découpage des images suivant des imagettes (également nommées zones d'intérêt), dont on suit le mouvement en vue d'obtenir un champ de déplacement discrétisé sur les centres de chaque imagette. Cette inter-corrélation d'imagettes possède de nombreux inconvénients. Chaque imagette fournit un déplacement moyen par inter-corrélation. Pour que le traitement fonctionne, la taille de l'imagette doit être grande pour contenir un nombre suffisant de bandes de Kikuchi. Afin d'obtenir un champ de déplacement assez riche, il est obligatoire de prendre un nombre important d'imagettes, en 30 pratique variant entre 25 et 100. Du fait de ce nombre élevé et de la taille également importante des imagettes, des chevauchements importants entre imagettes existent. D'une part, cela ralentit considérablement le traitement de corrélation car les pixels 2 d'une même imagette sont pris plusieurs fois en compte. D'autre part, ces chevauchements introduisent de fortes corrélations spatiales entre imagettes de sorte que l'augmentation du nombre d'imagettes, à partir d'un certain seuil, ne permet ni de diminuer l'incertitude de mesure, ni de réduire des biais systématiques. 5 En outre, le choix de positionnement des imagettes influence les résultats obtenus, ce qui est une manifestation de la non optimalité de la méthode utilisée. En cas de relativement grandes déformations, des rotations importantes des motifs peuvent être observées, supérieures à un degré. Le déplacement étant cherché par imagette au moyen d'une translation pure, 1 0 la base cinématique n'est alors qu'approximative pour reproduire la transformation observée sur les images. En conséquence, l'inter-corrélation provoque des erreurs importantes, par exemple, des contraintes évaluées dans les zones de fortes rotations. L'invention vise à obtenir un procédé et un dispositif de traitement d'images qui pallient les inconvénients de l'état de l'art de la technique en permettant de 15 calculer le champ de déplacement entre les images, avec une meilleure qualité et un coût de calcul plus bas. A cet effet, un premier objet de l'invention est un procédé de traitement d'images, obtenues par un détecteur à diffraction, d'un matériau cristallin ou polycristallin, dans lequel on mesure par le détecteur : 20 - une première image de diffraction du matériau se trouvant dans un état de référence, donnant des premières valeurs numériques de pixels en fonction de deux coordonnées de pixels, - au moins une deuxième image de diffraction du matériau se trouvant dans un état déformé par rapport à l'état de référence, la deuxième image donnant des 25 deuxièmes valeurs numériques de pixels en fonction des coordonnées de pixels (les valeurs numériques de pixels des images pouvant être des niveaux de gris ou des valeurs de luminance et étant notamment des valeurs de l'intensité du faisceau diffracté en un point du détecteur). Suivant l'invention, dans un calculateur, un champ de déplacement, pour 30 déplacer des pixels de la première image à des pixels d'une image déformée, en fonction : 3 - des deux coordonnées des pixels, - de coordonnées prédéterminées d'un centre, correspondant à la projection normale, dans le plan d'image du détecteur, d'un point source du faisceau diffracté dans le matériau, et 5 - de composantes d'un tenseur gradient de déformation élastique, est préenregistrée dans une mémoire. Au cours d'une première étape de calcul, on fait prendre au tenseur gradient de déformation élastique courant une valeur déterminée du tenseur gradient de déformation élastique. 10 Au cours d'une deuxième étape de calcul, on calcule le champ de déplacement courant à partir du tenseur gradient de déformation élastique courant et des coordonnées de pixels de la première image. Au cours d'une troisième étape de calcul, on calcule des troisièmes valeurs numériques de pixels d'une image déformée en corrigeant la deuxième image aux 15 coordonnées de pixels auxquelles a été additionné le champ de déplacement courant. Ainsi, l'image déformée est corrigée du champ de déplacement courant, et est appelé troisième image ou image déformée corrigée ayant des troisièmes valeurs de pixels, pour la distinguer des première et deuxième images. Au cours d'un algorithme itératif, on effectue des itérations des première, 20 deuxième et troisième étapes de calcul sur des valeurs déterminées de tenseur modifiées, jusqu'à remplir un critère de convergence sur la valeur déterminée du tenseur gradient de déformation élastique, pour calculer le champ de déplacement correspondant. On décrit ci-dessous des modes de réalisation des différentes étapes de ce 25 procédé. Ainsi, suivant un mode de réalisation, le détecteur à diffraction permet notamment de mesurer une configuration de référence de la géométrie de la diffraction, associant aux coordonnées de chaque pixel, une direction du faisceau électronique diffracté, et dont l'origine correspondant à la projection normale dans le plan d'image du détecteur, du point source du faisceau diffracté dans le matériau, 4 - une première figure de diffraction (de Kikuchi) du matériau se trouvant dans un état de référence, donnant une image dite de référence, soit une image en niveaux de gris dont la valeur numérique correspond en chaque pixel à une densité d'électrons diffractés dans la direction géométriquement reliée à la position du pixel, 5 - au moins une deuxième figure de diffraction se trouvant dans un état déformé par rapport à l'état de référence, cette deuxième image donnant des deuxièmes valeurs numériques de niveaux de gris en chaque pixel. Suivant un mode de réalisation, le champ de déplacement est défini en chaque pixel, pour apparier les pixels de la première image de diffraction aux pixels de 10 l'image de diffraction du cristal déformé, de manière à faire coïncider au mieux les niveaux de gris des pixels appariés, ce champ de déplacement étant ainsi associé à une direction du faisceau diffracté et prenant une forme algébrique dont l'expression est une fonction connue du tenseur gradient de déformation élastique. Suivant un mode de réalisation, au cours de la première étape de calcul, le 15 tenseur gradient de déformation élastique courant est initialisé à une valeur prédéterminée de tenseur et est actualisé au cours des itérations. Suivant un mode de réalisation, au cours de la troisième étape de calcul, on calcule les troisièmes valeurs numériques de pixel de la troisième image déformée, mentionnées ci-dessus, en interpolant la deuxième image aux coordonnées de pixels 20 auxquelles a été additionné le champ de déplacement courant. Dans le cadre de l'invention, on fournit les troisièmes valeurs de pixels d'une image déformée, qui sont une correction de la deuxième image, s'appuyant sur une majeure partie de l'image, en un temps plus court. Le calcul procuré par l'invention réduit en outre l'incertitude de mesure du champ de déplacement entre la première 25 image et la deuxième image. L'invention est ainsi spécialement adaptée, dans des modes réalisation, au calcul à partir de la troisième image, d'une déformation du matériau entre la première image et la deuxième image. e Suivant un mode de réalisation de l'invention, le tenseur r gradient de déformation élastique est égal à 5 Fie F2e F3e [Fe]= F4e F5e F6e F7e F8e 1 où file fr2e fr3e fr4e F5e fi6e fr7e fr8e sont les composantes du tenseur e F gradient de déformation élastique. Suivant un mode de réalisation de l'invention, le champ de déplacement, 5 uy, pour déplacer des pixels de la première image à des pixels d'une image déformée, en fonction : - des deux coordonnées x,y de pixels, - des coordonnées prédéterminées x*, y*, z* du centre, correspondant à la projection normale, dans le plan d'image du détecteur, du point source du faisceau 1 0 diffracté dans le matériau, et - des composantes Fie F2e ,F3e F4e F5 e F6e F7e 77e r8 du tenseur Êe gradient de déformation élastique, est égal à z* (Pie (x x*) - .Y*) fr3ez*) - x*) ux(x, = fr7e(x - x*) + fr8e(Y y*) z * z * CÈ4e(x - x*) + fr5e - y*) fr6ez*) (y _ y*) uy(x, y) = fr; (x - x*) fr8e(Y-y*) z * 15 Suivant un mode de réalisation de l'invention, l'algorithme itératif est effectué par une méthode de type Gauss-Newton. Suivant un mode de réalisation de l'invention, au cours de l'algorithme itératif, on calcule un vecteur (or ) de correction vérifiant l'équation [M]{6fre }= fY} 2 0 où [1%1] est une matrice hessienne de dimension 8 x 8, ayant comme coefficients 6 mi,(n -"= I(vf(x).svi (x, {fre}))(vf(x).(t. (x, {fre})) ROI où f représente les premières valeurs de pixels de la première image, x représente les deux coordonnées (x, y) de pixels, (Di (x; {fre }) est un champ de sensibilité du champ de déplacement par 5 rapport à la composante Fe du tenseur Ê e gradient de déformation élastique et est égal à la dérivée partielle du champ de déplacement par rapport à chaque composante Fe du tenseur Ê e gradient de déformation élastique, {y} est le vecteur de résidus ayant comme composantes y i(n) = (X) - (n) (X ,{fr e }))V f (X) - (1:. (X ,{fr e}) ROI 10 et (n)(x,{pels)\= (x) = g(x + u(x)) est la troisième valeur numérique de chaque pixel de l'image déformée et g représente les deuxièmes valeurs de pixels de la deuxième image, e le critère de convergence sur la valeur déterminée r de tenseur étant que la 15 norme du vecteur ore de correction soit inférieure à une borne positive prescrite Oc), non nulle, la valeur déterminée {F e } de tenseur étant incrémentée du vecteur { ( e de correction à chaque itération de la première étape de calcul. e Suivant un mode de réalisation de l'invention, le tenseur r gradient de 2 0 déformation élastique possède, comme composantes , huit composantes frle fr 2e fr 3e fr 4e fr e fr6e fr frge et une neuvième composante fixée à 1, à savoir selon l'équation suivante : 7 FFF3e 1 2 F4e F5e F6e Fe Fe 1 7 8 [Fe]= 5 10 le champ de sensibilité du champ u(x) de déplacement a comme composantes (te (x; {fre les composantes selon les équations suivantes : z * (x - x*) F.7(x - x*)+fr8e (y - y*) + z * z * - Y*) x2 =^ F7e(X X*) fr8e(y - y*) Z * (z*)2 (I)x3 =^ F7e (x - x*) + fr8e (y - y*) + z * (I) *(x X*)(file .X*) fie( y*) 3e 1*) (fr7e - .X*) fr8e - y*) + z *Y (I)x8 z * (Y Y*)(frie(x x*) fie( - Y*) fiez*) (fi;e(x - x*) + fr8e (y - y*) + z *Y 15 8 z * (x - x*) _ , F7e(X- x*) fr8e - y*) Z * * (.); Y*) F7e (X - .X*) F8e (y - y*) + z * (z*)2 = , F7e (x - x*) + Ê8e (y - y*) + z * z * (x _ x*)(fr4e (x x*) fr5e (y y*) fr6ez*) (fr7e (X - x*) fr8e (y - y*) + z *Y 5 (I) y8 * (y Y*)(F' 4e (X -X*) F' 5e (Y y*) F'6ez*) (fr7e (X - .X*) fr8e (y - y*) + z *Y Suivant un mode de réalisation de l'invention, le procédé de traitement d'images est exécuté pour la majeure partie ou pour l'ensemble des pixels de la première image et de la deuxième image. Suivant un mode de réalisation de l'invention, les première et deuxième 1 0 images sont obtenues après filtrage de valeurs de pixels surexposés en les remplaçant par une moyenne de pixels voisins de ceux-ci. On filtre ainsi des valeurs aberrantes (bruit classiquement défini comme « poivre et sel » qui intervient pour certains détecteurs). Suivant un mode de réalisation de l'invention, les première et deuxième 15 images sont obtenues après filtrage des valeurs de pixels par soustraction de tendances globales de niveau de gris, représentées par un polynôme d'ordre 2 ou 3 obtenu par un processus de régression. Suivant un mode de réalisation de l'invention, les première et deuxième images sont obtenues après filtrage des valeurs de pixels par un filtre de lissage 20 gaussien. Suivant un mode de réalisation de l'invention, le procédé de traitement d'images est exécuté pour plusieurs deuxièmes images. 9 Suivant un mode de réalisation de l'invention, l'on fournit sur une sortie au moins l'un parmi : - le tenseur gradient de déformation élastique, calculé comme remplissant le critère de convergence, 5 - un résidu, calculé comme étant la différence entre d'une part les troisièmes valeurs numériques de pixels de l'image déformée corrigée, ayant été calculées pour le champ de déplacement correspondant comme remplissant le critère de convergence sur la valeur déterminée du tenseur gradient de déformation élastique, et d'autre part les premières valeurs numériques de pixels de la première image, 1 0 - une moyenne quadratique des résidus calculés pour respectivement plusieurs deuxièmes images, - le champ de déplacement correspondant, calculé comme remplissant le critère de convergence sur la valeur déterminée du tenseur gradient de déformation élastique, 15 - une image déformée corrigée, correspondant aux troisièmes valeurs numériques de pixel ayant été calculées à partir du champ de déplacement correspondant, calculé comme remplissant le critère de convergence. Suivant un mode de réalisation de l'invention, les première et deuxième images sont obtenues après soustraction d'un fond d'image. 20 Un deuxième objet de l'invention est un dispositif de traitement d'images d'un matériau cristallin ou polycristallin, comportant un détecteur à diffraction permettant d'acquérir : - une première image de diffraction du matériau se trouvant dans un état de référence, donnant des premières valeurs numériques de pixels en fonction de deux 25 coordonnées de pixels, - au moins une deuxième image de diffraction du matériau se trouvant dans un état déformé par rapport à l'état de référence, la deuxième image donnant des deuxièmes valeurs numériques de pixels en fonction des coordonnées de pixels, caractérisé en ce que le dispositif comporte au moins un calculateur, 30 comportant au moins une mémoire, dans laquelle est enregistrée un champ de 10 déplacement, pour déplacer des pixels de la première image à des pixels d'une image déformée, en fonction : - des deux coordonnées de pixels, - de coordonnées prédéterminées d'un centre, correspondant à la projection 5 normale, dans le plan d'image du détecteur, d'un point source du faisceau diffracté dans le matériau, et - de composantes d'un tenseur gradient de déformation élastique, le calculateur étant configuré pour au cours d'une première étape de calcul, faire prendre au tenseur gradient de 10 déformation élastique courant une valeur déterminée, au cours d'une deuxième étape de calcul, calculer le champ de déplacement courant en à partir du tenseur gradient de déformation élastique courant et des coordonnées pour chacun des pixels de la première image, au cours d'une troisième étape de calcul, calculer des troisièmes valeurs 15 numériques de pixels d'une image déformée en appliquant la deuxième image aux coordonnées de pixels auxquelles a été additionné le champ de déplacement courant, au cours d'un algorithme itératif, effectuer des itérations des première, deuxième et troisième étapes de calcul sur des valeurs déterminées de tenseur réactualisées, jusqu'à remplir un critère de convergence sur la valeur déterminée du 20 tenseur gradient de déformation élastique, pour calculer le champ de déplacement correspondant. Un troisième objet de l'invention est un programme d'ordinateur, comportant des instructions de code pour la mise en oeuvre du procédé de traitement d'images d'un matériau cristallin ou polycristallin tel que décrit ci-dessus, lorsqu'il est exécuté 25 sur un calculateur. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif en référence aux dessins annexés, sur lesquels : - la figure 1 représente schématiquement un dispositif d'acquisition et de 30 traitement d'images suivant un mode de réalisation de l'invention, 11 - la figure 2 représente schématiquement en perspective la projection des faisceaux diffractés par un dispositif de traitement d'images suivant un mode de réalisation de l'invention ainsi que l'effet produit par le tenseur gradient de déformation élastique sur le faisceau diffracté, 5 - la figure 3 représente schématiquement en vue de côté la projection des faisceaux diffractés par un dispositif de traitement d'images suivant un mode de réalisation de l'invention, - la figure 4 représente schématiquement un organigramme du procédé de traitement d'images suivant un mode de réalisation de l'invention, 10 - La figure 5A montre une image d'une éprouvette, prise par un microscope électronique à balayage, - La figure 5B montre une image obtenue à partir de la figure 5A par un logiciel de l'état de l'art de la technique, - La figure 5C montre une échelle des couleurs de la figure 5B, 15 - les figures 6A et 6B montrent un exemple d'images de départ de diffraction d'une éprouvette, - la figure 6C représente la différence obtenue entre les images des figures 6A et 6B, - la figure 6D représente une valeur calculée pour un tenseur gradient de 2 0 déformation élastique apparaissant dans le procédé et le dispositif de traitement d'images suivant un mode de réalisation de l'invention, dans l'exemple des figures 6A et 6B, - les figures 6E, 6F, 6G, 6H et 61 représentent des images apparaissant dans différentes étapes du procédé et du dispositif de traitement d'images suivant des 25 modes de réalisation de l'invention, dans l'exemple des figures 6A et 6B, - les figures 7A, 7B, 7C, 7D et 7E représentent des composantes d'un tenseur de contrainte obtenues par le procédé et le dispositif de traitement d'images suivant des modes de réalisation de l'invention, à partir de la figure 5A, - La figure 8A représente une contrainte équivalente de von Mises obtenue 30 par le procédé et le dispositif de traitement d'images suivant des modes de réalisation 12 de l'invention, à partir de la figure 5A, et la figure 8B représente un histogramme de cette contrainte. A la figure 1, dans le procédé de traitement d'images suivant l'invention, au cours d'une étape préalable de mesure, on obtient par un détecteur 10 (ou capteur 5 10 ou caméra 10) à diffraction, des images numériques (ou clichés) d'un matériau CR cristallin ou polycristallin. Le détecteur 10 peut être un détecteur à diffraction rétrodiffusée en haute résolution d'images, dit EBSD ou HR-EBSD (en anglais : HR-EBSD pour High Resolution Electron Back-Scattered Diffraction). Bien entendu, d'autres types de détecteur 10 à diffraction peuvent être utilisés. 10 On mesure par le détecteur 10 une première image f du matériau CR se trouvant dans un état de référence, donnant des premières valeurs numériques f(x,y) de pixels (par exemple des niveaux de gris ou autres) en fonction des deux coordonnées x,y des pixels. On mesure également par le détecteur EBSD une ou plusieurs deuxième(s) image(s) g du matériau CR se trouvant dans un état déformé 15 par rapport à l'état de référence. La deuxième image g donne des deuxièmes valeurs numériques g(x,y) de pixel (par exemple des niveaux de gris) en fonction des deux coordonnées x,y des pixels. Les images f et g sont choisies lors de l'étape E2 de la figure 4. Aux figures 2 et 3, les 3 directions x, y et z sont orthonormales entre elles. Suivant un mode de réalisation de l'invention, l'état déformé et l'état de 2 0 référence sont obtenus en appliquant des sollicitations mécaniques différentes au matériau CR, par exemple en n'appliquant aucune sollicitation mécanique au matériau CR dans l'état de référence et en appliquant une sollicitation mécanique déterminée au matériau CR dans l'état déformé. Dans ce cas, il peut être prévu des moyens pour exercer et/ou contrôler une sollicitation mécanique sur le matériau CR. 25 Ainsi que représenté à la figure 1, le détecteur 10 fait partie d'un dispositif 1 de mesure comportant une source interne 2 d'émission d'un faisceau 3 de particules incidentes et le détecteur 10 à diffraction rétrodiffusée. La source 2 et le faisceau 3 de particules incidentes sont positionnés par rapport au matériau CR pour que le matériau CR émette, par interaction du matériau CR avec le faisceau 3 de particules 30 incidentes, un ou plusieurs faisceaux 4 de particules diffracté(s). Le faisceau 3 peut être ou comprendre un faisceau de rayons X, ou un faisceau d'électrons, ou autres. 13 Le faisceau 4 diffracté peut également correspondre à une diffraction de type Kossel, une diffraction de type Laue ou une diffraction dans un microscope électronique à transmission (notamment à diffraction TKD, qui est l'abréviation en anglais de « Transmission Kikuchi Diffraction », diffraction Kikuchi en transmission, réalisée 5 dans un microscope électronique à transmission). Le détecteur est positionné par rapport à la source 2, au faisceau 3 de particules incidentes et au matériau CR, pour recevoir ou intercepter le ou les faisceaux 4 de particules diffracté(s), issus du matériau CR en réponse au faisceau 3 de particules incidentes. Le détecteur 10 comporte par exemple un écran 11 de réception du ou des faisceau(x) 4 de particules 10 diffracté(s). L'écran 11 de réception possède une étendue déterminée, par exemple bidimensionnelle et plane, ou autre. L'écran 11 de réception est relié à une unité 12 d'enregistrement et de production des images f, g (pouvant être une unité 12 de calcul automatique) à partir du ou des faisceau(x) 4 de particules diffracté(s), reçu(s) ou intercepté(s) par l'écran 11 de réception. Le faisceau 3 de particules incidentes 15 peut être ou contenir un rayonnement monochromatique incident de longueur d'onde déterminée ou plusieurs rayonnements monochromatiques incidents de longueurs d'ondes déterminées différentes. Le détecteur 10 peut faire par exemple partie d'un microscope électronique, notamment un microscope électronique à transmission et/ou à balayage. L'écran 11 peut être par exemple un écran phosphorescent, ou 20 autre. Suivant des modes de réalisation de l'invention, le phénomène de diffraction du faisceau 3 de particules incidentes par le matériau CR donne lieu sur le détecteur 10 à un diagramme de Laue dans le cas où le faisceau 3 de particules incidentes est un faisceau de rayons X, et à un cliché de Kikuchi dans le cas où le faisceau 3 de 25 particules incidentes est un faisceau d'électrons. D'après la loi de Bragg, le rayonnement monochromatique incident 3 de longueur d'onde X. va diffracter sur les plans cristallins hkl (ou plans diffractants hkl) du matériau CR en respectant : = 2dhkr sin 0 (1) 30 où n est l'ordre de diffraction et dhki la distance inter-réticulaire des plans hkl. L'angle 0 est le demi-angle entre le faisceau incident 3 et le faisceau diffracté 4. Un 14 matériau CR cristallin dont la maille élémentaire est décrite par le repère (a, b, c) possède un réseau réciproque dont le repère est (a*, b*, c*). Par définition, a* est tel que a - a* = 1 , b - a* = c - a* = 0. Idem pour b* et c* par permutation circulaire. 5 Soit ki le vecteur d'onde du faisceau incident et kf le vecteur d'onde du faisceau diffracté. La diffraction résultant de processus élastiques (sans perte d'énergie), la norme des vecteurs est identique, kf = = X.-1, et donc kf décrit un cercle de rayon égal à la norme de k. Le vecteur de diffraction q, défini comme : q = kf - (2) 1 0 est donc un vecteur qui, si son origine est positionnée à l'extrémité du vecteur k, parcourt alors un cercle de même rayon 1/), centré sur l'origine du vecteur k. Dans ce même espace réciproque, la condition de Bragg peut s'exprimer par la condition de coïncidence entre le vecteur dila) de la famille de plans cristallins hkl, 15 et le vecteur de diffraction q. Seules quelques directions favorables donnent lieu à une diffraction cohérente. Suivant un mode de réalisation de l'invention, le faisceau d'électrons adopté en EBSD a une longueur d'onde X, très faible et 0 est généralement inférieur à 2°. Les électrons en condition de diffraction restent donc à proximité de la trace du plan 2 0 cristallin 40 sur lequel ils diffractent, ainsi que représenté à la figure 1. Suivant un mode de réalisation de l'invention, dans le cas où la source 2 est divergente, les faisceaux diffractés 4 se répartissent sur deux cônes de Kossel CO de manière symétrique par rapport à la trace du plan cristallin 40. Les cônes de Kossel CO captés par le détecteur 10 peuvent avoir une forme de deux hyperboles 25 extrêmement ouvertes, ainsi que représenté à la figure 1. Suivant un mode de réalisation de l'invention, la première image f de l'état de référence et la deuxième image g de l'état déformé sont obtenues lors d'une unique acquisition par balayage de la surface d'un échantillon du matériau CR. L'état de référence est alors pris au centre du grain (supposé le moins contraint) que l'on 30 compare alors avec des clichés de diffraction acquis en périphérie du grain. 15 Suivant un mode de réalisation de l'invention, aux figures 2 et 3, un centre 0 correspondant à la projection normale, dans le plan 110 d'image de l'écran 11 du détecteur 10, d'un point source S du faisceau diffracté 4 dans le matériau CR, a deux coordonnées prédéterminées (x*, y*) dans le plan 110 de l'écran 11. La projection 5 du ou des faisceau(x) diffracté(s) 4 sur le détecteur 10 est illustrée dans un exemple non limitatif aux figures 2 et 3. Le détecteur 10 est choisi comme référentiel dans toute la suite. A la figure 2, le coin B en bas à gauche de l'écran 11 plan est considéré comme origine, l'axe x est l'axe horizontal dans le plan de l'écran 110, l'axe y est l'axe vertical dans le plan de l'écran 110, et l'axe z est l'axe normal au 10 plan de l'écran 110. Le faisceau 3 d'électrons incidents est dirigé sur l'éprouvette avec une inclinaison de 70° par rapport à la normale et un point "central" S dans le volume d'interaction du matériau CR est considéré comme source "effective" des électrons diffractés 4. Le centre de projection 0 est la projection normale (par rapport au plan de l'écran 11) du point S sur le plan 110 de l'écran 11. Ses 15 coordonnées sont notées (x*, y*, 0). On note z* la distance entre le point source S et le centre 0, donc S a pour coordonnées (x*, y* , -z*). Suivant un mode de réalisation de l'invention, on note F un tenseur gradient de déformation selon l'équation ci-dessous, qui, appliqué à un point courant X dans la configuration initiale du matériau CR de référence, permet d'obtenir sa position x 2 0 dans la configuration déformée : dx [F]= Lw] = ax ax ax ax a Y az ay ay ay ax a Y az az az az _ax a A titre d'exemple, la figure 1 représente la diffraction du faisceau dans la chambre du microscope, se projetant sur l'écran de la caméra EBSD, représenté 25 pour deux états cristallins et le tenseur gradient de transformation associé pour 16 passer d'un état à l'autre. La figure 3 détaille la projection, et représente le point source S et le centre de projection. Suivant un mode de réalisation de l'invention, le tenseur gradient de déformation F est le produit de deux parties : la partie élastique Fe (tenseur gradient 5 de déformation élastique) et la partie plastique FP selon l'équation suivante : F = Fe FP Suivant un mode de réalisation de l'invention, la déformation plastique a pour effet de rendre les bandes de Kikuchi moins nettes, un effet difficile à quantifier et qui n'autorise pas la mesure de la déformation plastique par l'analyse d'images de 1 0 diffraction. En cas de déformation élastique hydrostatique, on observe une variation de la largeur des bandes de Kikuchi sur le détecteur. Au contraire, la déformation élastique déviatorique modifie la forme de la maille cristalline, c.-à-d. l'orientation relative des plans cristallins, donc les relations angulaires à l'intérieur du cristal. La projection sur un écran loin de la source amplifie l'écart angulaire entre deux 15 faisceaux et donne lieu à une variation mesurable entre les clichés. Suivant un mode de réalisation de l'invention, comme seules huit composantes de Fe sont mesurables, il faut définir une convention fixant ce degré de e liberté manquant. Suivant un mode de réalisation, on choisit le seul tenseur r ayant huit composantes F,fr2e fr3e fr4e fr5e fr6e fr7e fige et une neuvième 20 composante fixée à 1, à savoir selon l'équation ci-dessous : Fe Fe Fe 1 2 3 [F e = Fe Fe Fe 4 5 6 F7e F8e 1 _ _ Suivant un mode de réalisation de l'invention, le gradient de déformation élastique vrai est alors selon l'équation ci-dessous : Fe = Fe 17 Le facteur aziaz n'est pas mesurable mais peut être déterminé par des hypothèses supplémentaires (comme par exemple le choix d'un état de contrainte plane qui est souvent adopté). Suivant un mode de réalisation de l'invention, considérons un élément 4X 5 dans le réseau cristallin de référence du matériau CR. Supposons que cette direction particulière donne lieu à un faisceau diffracté 4 qui intersecte le détecteur 10 au point P, de coordonnées (x,y,0). Dans la suite, nous notons p le vecteur SP. Nous avons (p) = (4x, 4y, 4z)T = (x - x*,y - y*,z*)T = a4X, où a est l'échelle de projection. Si le volume d'interaction autour du point source S est soumis à une 10 déformation élastique décrite par le tenseur gradient de déformation élastique Fe, l'élément de référence 4X est transformé en 4x = Fe4X. Le rayon diffracté 4 initialement dans la direction p est redirigé vers le point P", tel que p" = aFe4X où comme précédemment le vecteur SP" est noté p". La droite SP" 15 intersecte le détecteur 10 au point P'. Le déplacement apparent u sur le détecteur 10 est u = p'- p (c'est-à-dire le vecteur PP') selon les équations suivantes : * u = p'-p"= z p" p = z e aFeAX - p pz a(F AX)z z* Fei) p = F- ep p (Fep)z (Fep)z Une fonction de déplacement, donnant le champ de déplacement (ux,uy) pour 2 0 déplacer (ou « advecter » ou apparier) les pixels de la première image f aux pixels d'une image déformée g, en fonction : - des deux coordonnées (x,y), - des deux coordonnées prédéterminées x*, y* du centre 0, correspondant à la projection normale, dans le plan 110 d'image du détecteur 10, du point source S 25 du faisceau diffracté dans le matériau, et - de composantes Fe , Fe , Fe , Fe e 4 F5 F6e F7 e , F8e du tenseur Fe gradient de déformation élastique, 18 est préenregistrée dans une mémoire de l'unité 12. Bien entendu, l'expression « chacun des pixels » ou « les pixels » peut être remplacée par « des pixels » pour désigner une partie des pixels de l'image. Suivant un mode de réalisation de l'invention, le champ (ux,uy) de 5 déplacement est égal à, selon les équations (9) ci-dessous : ux(x, y) = z * (fije(x - x*) + Ê2e (y - y*) + Ê3ez*) (x - x*) fr7e (x - x*) + Ê8e (y - y*) + z * (X,Y) = z * ( fi'4e(x - x*) + fr5e (Y Y*) fr6eI*) (y - y*) fr7e (x x*) fr8e (y - y*) Z * Le champ de déplacement (ux,uy) entre la première image f de diffraction, appelée image de référence, et la deuxième image g de diffraction, appelée image déformée, reflète la déformation élastique de la maille cristalline du matériau CR au 1 0 point étudié. On décrit ci-dessous un algorithme pour calculer le tenseur r gradient de déformation élastique. Cet algorithme peut être par exemple itératif sur les étapes E7, E8, E9, El0 et Ell décrites ci-dessous, ainsi que représenté à la figure 4. Le procédé de traitement d'images est exécuté par l'unité 12 de calcul 15 automatique, telle que par exemple un ou plusieurs calculateur(s) et/ou un ou plusieurs ordinateur(s), et/ou un ou plusieurs processeur(s) et/ou un ou plusieurs serveur(s) et/ou une ou plusieurs machine(s), pouvant être programmées à l'avance par un programme informatique préenregistré. Au cours de la première étape E6 de calcul, on fait prendre au tenseur e 2 0 r gradient de déformation élastique courant une valeur déterminée r du tenseur gradient de déformation élastique, comme par exemple une valeur initiale INIT avant la première itération. Suivant un mode de réalisation de l'invention, cette valeur initiale INIT avant la première itération peut être par exemple le tenseur identité. 25 Au cours d'une deuxième étape E7 de calcul, postérieure à l'étape E6, on calcule le champ (ux,uy) de déplacement courant en appliquant la fonction de 19 déplacement dépendant du tenseur Fe gradient de déformation élastique courant aux deux coordonnées x,y pour chacun des pixels de la première image! On calcule ainsi le champ de déplacement courant induit par le tenseur gradient de déformation élastique courant aux coordonnées de pixels de la première image. 5 Au cours d'une troisième étape E8 de calcul, postérieure à l'étape E7, on calcule des troisièmes valeurs numériques d'une image déformée g(x) = g(x ± 11(x)) en appliquant ou interpolant la deuxième image g aux deux coordonnées des pixels auxquelles a été additionné le champ de déplacement (ux,uy) courant, où x désigne les deux coordonnées x,y des pixels. Le champ de 10 déplacement courant (ux,uy) est ainsi calculé en fonction du tenseur r gradient de déformation élastique courant et des coordonnées pour chacun des pixels de la première image au cours de la deuxième étape E8 de calcul. Suivant un mode de réalisation de l'invention, on calcule une fonction coût A 15 (voir ci-dessous), à partir des différences quadratiques, pour la majeure partie ou l'ensemble des pixels, entre les troisièmes valeurs numériques de pixel de l'image déformée g(x) = g(x ± 11(x)) , ayant été calculées, et les premières valeurs numériques f(x) de pixels de la première imagef, on calcule le tenseur r gradient de déformation élastique courant 2 0 minimisant la fonction coût A, calculée à partir des différences, pour la majeure partie ou l'ensemble des pixels, entre les troisièmes valeurs numériques de pixel de l'image déformée g(x) = g(x ± 11(x)) , ayant été calculées, et les premières valeurs numériquesf(x) de pixels de la première image f Au cours de l'algorithme itératif, on effectue des itérations des première, 25 deuxième et troisième étapes E7, E8, E9, E10, El 1 de calcul en modifiant à chaque itération les valeurs déterminées de tenseur r (appelées valeurs déterminées de tenseur modifiées Fe), jusqu'à remplir un critère de convergence sur cette valeur 20 déterminée Fe du tenseur gradient de déformation élastique, autrement dit jusqu'à remplir un critère de convergence sur la correction à la valeur courante de Fe . On calcule lors de l'étape E7 par la fonction de déplacement le champ de déplacement e (71,,/,/y) qui correspond à ces valeurs déterminées de tenseur r . 5 Dans l'état de la technique, la corrélation d'images numériques globale consiste à corréler les deux images f et g. Suivant un mode de réalisation de l'invention, le procédé met en oeuvre une corrélation d'images numériques intégrée (en abrégé CINI). Dans le cas d'une diffraction électronique, le procédé suivant l'invention peut être dénommé en abrégé ADDICTED pour Alternative Dedicated 1 0 Digital Image Correlation Tailored to Electronic Diffraction, soit en français corrélation d'images numériques alternative, dédiée et adaptée à la diffraction électronique. Suivant un mode de réalisation de l'invention, la recherche du champ de déplacement (ux,uy) est réalisée en corrigeant l'image déformée g par le champ de 15 déplacement gu(x)= g(x±u(x)) lors de l'étape E8, de façon à s'approcher le ' plus possible de l'image de référence g'(x)= f(x) , ce qui revient à minimiser la norme quadratique du résidu r = (f(x)- g'(x)) sommé sur la région d'intérêt entière (ou fonction coût A), pouvant être les images entières f et g. Dans ce qui suit, l'image déformée g, ayant été corrigée par le champ de déplacement 2 0 gu (X) = g(x + u(x)) lors de l'étape E8 est appelée image déformée corrigée g. Le procédé de traitement d'images suivant l'invention peut donc être un procédé de correction de la ou des deuxième(s) image(s) g et le dispositif de traitement d'images suivant l'invention peut donc être un dispositif de correction de 25 la ou des deuxième(s) image(s) g. La fonction coût A à minimiser s'écrit selon l'équation suivante : 21 A= ROI ogu(x)= g(x + u(x)) est l'image gu déformée, corrigée du champ de déplacement (ux,uy) (au cours de l'algorithme de résolution numérique, il s'agira de la détermination courante) et la minimisation de cette fonction coût A conduira à 5 des corrections successives de cette détermination jusqu'à convergence. Suivant un mode de réalisation de l'invention, les champs (ux,uy) de déplacement (ou leurs corrections) sont recherchés comme des combinaisons linéaires (ou affines) de champs constituant une « base cinématique ». Ils sont susceptibles d'avoir un support quelconque et nécessitent alors un traitement global 10 de l'ensemble de la cinématique. Si de plus la dite base cinématique est issue d'un modèle physique, et donc restreinte à des mécanismes bien identifiés ayant une signature en champ de déplacement, alors elle sera dite "intégrée". La corrélation d'images numériques qui en résulte est alors également qualifiée d'intégrée. Dans le cas de la méthode ADDICTED, la CINI est l'outil choisi pour 15 analyser les figures de diffraction (Jet g sont les figures de diffraction respectives du matériau CR cristallin ou polycristallin de référence (libre de contrainte pour l'image f) et du cristal observé (pour l'image g) . Suivant un mode de réalisation de l'invention, on explicite la manière dont le champ de déplacement (ux,uy) observé dans la figure de diffraction dépend ici de la déformation élastique du matériau CR 20 diffractant Fe . Enfin, d'autres paramètres, P, peuvent influer sur le champ de déplacement (ux,uy) mesuré, par exemple des paramètres P de projection, comme l'inclinaison de l'éprouvette par rapport au capteur 10, la taille physique d'un pixel et le pas de balayage du microscope électronique à balayage, ces deux derniers étant résumés par le terme de décalage induit par le faisceau (en anglais : beam-induced 25 shift). Suivant un mode de réalisation de l'invention, on explicite le champ de déplacement (ux,uy) selon l'équation suivante : U = U(X, P) 22 Enfin, génériquement, ces champs peuvent avoir une dépendance non linéaire en Fe ou P. Suivant un mode de réalisation de l'invention, on développe ces dépendances autour de la détermination courante des paramètres pour accéder à une expression affine adaptée aux approches CINI. Ainsi, suivant un mode de réalisation 5 de l'invention, on écrit selon les équations suivantes : U = U(X;F^ e 6F^ P + 6P) = u(x; P) + (x; e,P)ke +If j (X; î'e,P)ÔPJ au '11 = ap, où (1), correspond au champ de sensibilité par rapport à la composante 1 0 (pour i allant par exemple de 1 à 8) du tenseur gradient r de déformation, et Ti est le champ de sensibilité par rapport au paramètre 131. Ainsi, par exemple, (I) est une matrice de taille (2Np, xel) x 8. Suivant un mode de réalisation de l'invention, le champ de sensibilité (te (X; (fr e }) du champ u(x) de déplacement a comme composantes les composantes F^e) et 1:Dyi (X.F^e) selon les équations suivantes : z* (X - X*) (1)xl fr7e (X - X*) fr8e (y - y*) + z * 1* (Y - Y*) = x2e (X - X*) fr8e (y - y*) + z * (z*)2 F - x*) + 'fr8e(y - y*) + z * =0 =0 15 20 23 0 * X*)( file .X*) fie( y*) -fr3ez*) (fr7e - .X*) fr8e - y*) + z *Y z * (y - Y*)(fr le (X - )(*) fie( - y*) fr 3e 1*) (I)x8 (fi;e(x - x*) + fr8e (y - y*) + z 5 z * (x - x*) - , F7e (x - x*) + Ê8e (y - y*) + z * z * - y*) 10 (1)y6 Fe (X - .X*) F8e (y - y*) + z * (z*)2 = , F7e(x- x*) + Ê8e (y - y*) + z * z * ()C X*)(fr4e (X .X*) fie( y*) F6ez*) (fi;e(x - x*) + fr8e (y - y*) + z *Y 1* (.32 Y*)(F'4e (x .X*) F' 5e (Y y*) F'6ez*) (fr7e (X - .X*) fr8e (y - y*) + z *Y Suivant un mode de réalisation de l'invention, l'algorithme itératif est de type Gauss-Newton. On cherche, par exemple, à minimiser la fonction coût A de façon 15 itérative par un algorithme de Gauss-Newton. Suivant un mode de réalisation de l'invention, au cours de l'étape E8 postérieure à l'étape E7, on calcule (1)y8 24 k"(n) , P)= k",(x) = g(x + u(x,Êe , P)) , qui est la troisième valeur numérique de pixel de l'image déformée corrigée. Suivant un mode de réalisation de l'invention, au cours de l'étape E8 (n) postérieure à l'étape E7, on calcule g. (x, F e= g , (x) = g (x + u(x, Êe)), 5 qui est la troisième valeur numérique de pixel de l'image déformée corrigée. Suivant un mode de réalisation de l'invention, au cours de l'algorithme itératif, au cours de l'étape El0 postérieure à l'étape E9, on calcule un vecteur {6fr e de correction vérifiant l'équation suivante : [1'1]{15 e = fYI 10 où [1%1] est une matrice hessienne de dimension 8 x 8. Suivant un mode de réalisation de l'invention, au cours de l'étape E9 postérieure à l'étape E8, on calcule les coefficients (n-1) ij de la matrice hessienne [1%1] selon l'équation suivante : M 1(.11 1 -1) = 1(V f (X) Se {ft e f (X) Se {ft e })) ROI où f est les valeurs de pixels de l'image de référence, x représente les deux 15 coordonnées x et y des pixels. (1), (x; (fre p est le champ de sensibilité du champ de déplacement (ux,uy) par rapport à la composante Fe du tenseur Ê e gradient de déformation élastique et est égale à la dérivée partielle du champ (ux,uy) de déplacement par rapport à chaque composante Fe du tenseur Ê e gradient de déformation élastique. Le champ de sensibilité (1), (x; (fre p est donc un vecteur à e 20 deux dimensions, qui a une composante en x, égale à (1) (x; F ) et une composante en y, égale à {y} est le second membre de la méthode de Newton-Raphson et est le vecteur de résidus ayant comme composantes 25 y(n) = (f(x) - k'(n)(x, {fi'e }))vf(x) (D, (x, {fi'e ROI \7f(x) est le gradient def(x), c'est-à-dire de f(x,y). Suivant un mode de réalisation de l'invention, au cours de l'étape E9 postérieure à l'étape E8, on calcule les coefficients y,(n) du second membre {y} 5 selon l'équation suivante : y i(n) =IV(X)-:e(n)(X,{fre}))vf(X)*(1)/(X,{fre}) ROI Suivant un mode de réalisation de l'invention, le critère de convergence sur la valeur déterminée r de tenseur étant que la norme du vecteur 16Fe de 1 0 correction ayant été calculé est inférieure à la borne positive prescrite 8,, non nulle. Suivant un mode de réalisation de l'invention, lors de l'étape El 1 postérieure à l'étape E10, on examine après chaque itération des étapes E7, E8, E9 et E10, si la norme du vecteur {or } de correction est inférieure à une borne positive prescrite non nulle. Par exemple, la borne Ô, peut être égale à 10-7 ou autres. Dans le cas 15 où la norme du vecteur {or } de correction ayant été calculé n'est pas inférieure à la borne positive prescrite Ô, (cas NON de la figure 4), on exécute l'étape E12 de mise à jour décrite ci-dessous. Suivant un mode de réalisation de l'invention, au cours de l'étape E12 de calcul, on fait prendre au tenseur r gradient de déformation élastique courant une 2 0 valeur déterminée {r17^ej de tenseur. Par exemple, lors de l'étape E12 de calcul de 117^ e mise à jour, la valeur déterminée (r de tenseur est incrémentée du vecteur Wej de correction, selon l'équation de mise à jour de la valeur déterminée {frej, suivante : 26 {Fe }= {Fe }+ {6fre } L'étape E12 est suivie de la première étape E7 de calcul de l'itération suivante. Suivant un mode de réalisation de l'invention, dans le cas où la norme du 5 vecteur {ore } de correction ayant été calculé est inférieure à la borne positive prescrite Ô, (cas OUI de la figure 4), on exécute l'étape E13 décrite ci-dessous. Suivant un mode de réalisation de l'invention, lors de l'étape E13, on retient e les valeurs déterminées du tenseur r gradient de déformation élastique calculées au cours de l'étape E12 de la dernière itération, qui est celle vérifiant le critère de 1 0 convergence, et le champ de déplacement (ux,uy) calculé au cours de la dernière étape E7, qui est calculé par le champ de déplacement (ux,uy) qui correspond à ces e valeurs déterminées de tenseur r , ayant été retenues. Suivant un mode de réalisation de l'invention, lors de l'étape E13, l'on fournit sur une sortie (pouvant être par exemple un affichage sur un écran, et/ou un 15 enregistrement dans une mémoire, et/ou l'envoi sur un port de sortie ou autres) un ou plusieurs parmi : - le champ de déplacement (ux,uy) correspondant, calculé comme remplissant e le critère de convergence sur la valeur déterminée F du tenseur gradient de déformation élastique, e 2 0 - le tenseur r gradient de déformation élastique, calculé comme remplissant le critère de convergence, - l'image corrigée ou image déformée corrigée gu , correspondant à l'image déformée g, ayant été corrigée par le champ de déplacement (ux,uy) correspondant, calculé lors de l'étape E8 comme remplissant le critère de convergence sur la valeur e 25 déterminée r du tenseur gradient de déformation élastique, c'est-à-dire les troisièmes valeurs numériques de pixel ayant été calculées à partir du champ de 27 déplacement (ux,uy) correspondant, calculé comme remplissant le critère de e convergence sur la valeur déterminée r du tenseur gradient de déformation élastique, - un champ de résidu r, calculé comme étant la différence entre d'une part les 5 troisièmes valeurs numériques de pixels de l'image déformée corrigée g' , ayant été calculées pour le champ de déplacement (ux,uy) correspondant comme remplissant le e critère de convergence sur la valeur déterminée r du tenseur gradient de déformation élastique, et d'autre part les premières valeurs numériques f(x,y) de pixels de l'image de référence f, 10 - une moyenne quadratique des résidus r calculés pour respectivement plusieurs deuxièmes images g. Suivant un mode de réalisation de l'invention, le procédé est exécuté pour la majeure partie (par exemple au moins 50 %, au moins 60%, au moins 70%, au moins 80% ou au moins 90% des pixels de la première image f et de la deuxième image g) 15 ou l'ensemble des pixels de la première image f et de la deuxième image g, et ce en une seule fois pour chaque image. Suivant un mode de réalisation de l'invention, le procédé est exécuté pour plusieurs deuxièmes images g. Le champ de résidu r collecte tous les artéfacts de la formation et 2 0 l'acquisition des images, et ainsi contient potentiellement des informations très riches des images analysées. Les signaux dans le résidu qui ne sont pas un bruit blanc indiquent souvent une exploitation incomplète de l'information ou un modèle de transformation cinématique entre images non adapté. Le résidu peut donc contenir un signal mais aussi des bruits blancs. Si l'on suppose que l'image de référence f est 25 formée à partir d'une image de référence parfaite fp sans bruit, superposée à un bruit bf supposé de distribution normale, sans corrélation spatiale, et bien connu sous le terme de « bruit blanc gaussien » : alors l'image f est égale à, selon l'équation suivante : f= fp+ bf 28 Suivant un mode de réalisation de l'invention, pour une série g d'images g de diffraction déformées, chaque image g' est égale à, selon l'équation suivante : g = g pl + big bi est le bruit pur. g 5 Après l'application du procédé suivant l'invention, on a, selon l'équation suivante : g g Le résidu est alors, lorsque l'image de référence et l'image déformée correspondent parfaitement après le recalage (dans ce cas f ) Ç p selon les 10 équations suivantes : f = fp +bf b f - g Suivant un mode de réalisation de l'invention, la même image de référence est utilisée au sein d'un même grain, ce qui correspond donc à une seule image f, mais une multitude d'images g. On observe alors que les résidus à convergence 15 obtenus pour l'ensemble des images g, peuvent être moyennés et, en notant (...) cette moyenne sur les différents points de mesure au sein d'un grain, la moyenne des résidus r est égale à, selon les équations suivantes : (ri ) bf (17)g ) bf Ainsi, suivant un mode de réalisation de l'invention, la moyenne des résidus r 20 fournit une estimation du bruit bf, et on soustrait de l'image de référencefp le bruit bf ainsi calculé. Ceci permet de diminuer les incertitudes de mesure par la diminution du résidu (variance divisée par deux). Suivant un mode de réalisation de l'invention, lors de l'étape E5 précédant l'étape E6, on filtre les valeurs de pixels f(x,y) de la première image f et les valeurs 25 de pixels g(x,y) de la deuxième image g par un filtre gaussien. Un lissage gaussien à très petite échelle de l'image de diffraction facilite beaucoup le calcul ultérieur, car l'image brute est corrompue par un bruit blanc très important. Ce lissage gaussien 29 consiste à convoluer l'image initiale f (et éventuellement g) par une fonction émolliente, G = f *G où G(x,y) est choisie comme une gaussienne, dépendant d'une longueur 5 interne et qui s'écrit : 1 G(x,y)= e 42 27r2 En pratique, la longueur choisie peut être de 1 à 2 pixels. Ce choix dépend du niveau de bruit sur l'image, lui-même fonction du temps d'acquisition, du nombre de pixels, des paramètres du faisceau, du matériau imagé. Le lissage gaussien élimine 10 de façon efficace le bruit haute fréquence de l'image de diffraction. Suivant un mode de réalisation de l'invention, lors de l'étape E4 précédant l'étape E5, on filtre les valeurs de pixels f(x,y) de la première image f et les valeurs de pixels g(x,y) de la deuxième image g par un filtrage de valeurs de pixels surexposés en les remplaçant par une moyenne de pixels voisins de ceux-ci. Dans le 15 cas de la non-uniformité de l'écran phosphorescent 110 utilisé pour capter les électrons, les images acquises f, g possèdent souvent en des points fixes des pixels surexposés. Ces pixels très brillants sont néfastes au calcul de corrélation d'images et il est recommandé de remplacer la valeur de niveau de gris de ces "mauvais pixels" par le niveau de gris moyen des voisins. 20 Suivant un mode de réalisation de l'invention, lors de l'étape E4 précédant l'étape E5, on filtre les valeurs de pixels f(x,y) de la première image f et les valeurs de pixels g(x,y) de la deuxième image g par un filtrage de ces valeurs de pixels par soustraction de tendances globales de niveaux de gris, représentées par un polynôme d'ordre 2 ou 3 obtenu par un processus de régression. En effet, les images f et g 25 peuvent présenter des variations globales de niveaux de gris liées à des fluctuations de l'énergie moyenne diffractée. Pour la corrélation d'images globale, il est préférable dans ce cas de soustraire les tendances globales de niveaux de gris. Suivant un mode de réalisation de l'invention, lors de l'étape E3 précédant l'étape E4, on soustrait aux valeurs de pixels f(x,y) de la première image f et aux X2 +y2 30 valeurs de pixels g(x,y) de la deuxième image g un fond d'image, ayant été calculé. En effet, les images acquises f et g peuvent être entachées d'un fond dû en particulier à la distribution d'énergie des électrons diffractés. Cette correction par soustraction du fond élimine les variations d'intensité moyennes de l'image, ce qui permet ensuite 5 d'augmenter le contraste et de rendre les bandes de Kikuchi plus nettes. Par exemple, afin d'y parvenir, on réalise une acquisition à faible grandissement du matériau CR cristallin ou polycristallin, afin de scanner un grand nombre de grains d'orientations différentes. Les images de diffraction ainsi acquises sont moyennées et une estimation du fond peut être obtenue. En principe, l'image de fond est unique 10 pour toute la zone étudiée d'une acquisition HR-EBSD. Bien entendu, dans d'autres modes de réalisation de l'invention, l'ordre des étapes E3, E4 et E5 pourrait être modifié ou l'une ou plusieurs de ces étapes E3, E4 et E5 pourrait être supprimée(s). L'étape E2 peut être précédée d'une étape El d'entrée pour entrer des images de Kikuchi, le centre 0 et les paramètres P. L'étape 15 E13 peut être suivie d'une étape de post-traitement, par exemple pour obtenir d'autres grandeurs calculées à partir des grandeurs obtenues à l'étape E13, comme par exemple une composante de déformation eyy suivant la direction y, ou autres. Suivant un mode de réalisation de l'invention, le dispositif 1 de traitement d'images d'un matériau cristallin ou polycristallin suivant l'invention, comporte des 20 moyens pour la mise en oeuvre du dispositif 1 de traitement d'image suivant l'invention. Le dispositif 1 de traitement d'images suivant l'invention comporte le détecteur 10 à diffraction permettant d'acquérir : - une première image f du matériau se trouvant dans un état de référence, donnant des premières valeurs numériques f(x,y) de pixels en fonction de deux 25 coordonnées (x,y) de pixels, - une ou plusieurs deuxième image g du matériau se trouvant dans un état déformé par rapport à l'état de référence, la deuxième image donnant des deuxièmes valeurs numériques g(x,y) de pixels en fonction des coordonnées (x,y) des pixels. Le dispositif 1 de traitement d'images comporte au moins un calculateur 12, 30 comportant au moins une mémoire, dans laquelle est enregistrée un champ (ux,uy) de 31 déplacement, pour faire passer les pixels de la première image aux pixels d'une image déformée, en fonction : - des deux coordonnées (x,y) des pixels, - de coordonnées prédéterminées (x*, y*) du centre 0, correspondant à la 5 projection normale, dans le plan d'image du détecteur, du point source S du faisceau diffracté dans le matériau, et - des composantes Fie ,F2e ,F3e ,F4e ,F5 e ,F6e ,F e e) du tenseur Êe gradient de déformation élastique. Le calculateur 12 est configure pour : 10 au cours d'une première étape de calcul E6 ou E12, faire prendre au tenseur 17^ e F gradient de déformation élastique courant une valeur déterminée I (r au cours d'une deuxième étape E7 de calcul, calculer le champ (ux,uy) de déplacement courant correspondant au tenseur {Fe gradient de déformation élastique courant et aux coordonnées x,y de pixels de la première imagef, 15 au cours d'une troisième étape E8 de calcul, calculer des troisièmes valeurs numériques de pixel d'une image déformée --g,i(x)= g(x±u(x)) en corrigeant la deuxième image g aux coordonnées de pixels du champ de déplacement (u,,ity) courant, au cours d'un algorithme itératif, effectuer des itérations des première, 20 deuxième et troisième étapes de calcul E12, E7, E8 sur des valeurs déterminées de tenseur modifiées r , jusqu'à remplir un critère de convergence sur la valeur déterminée { Fe , pour calculer le champ de déplacement (ux,uy) correspondant. L'invention concerne également un programme d'ordinateur, comportant des instructions de code pour la mise en oeuvre du procédé de traitement d'images d'un 25 matériau cristallin ou polycristallin suivant l'invention, lorsqu'il est exécuté sur un calculateur. Le programme d'ordinateur est enregistré dans une mémoire du calculateur 12. 32 Suivant un mode de réalisation de l'invention, les première et deuxième images f et g sont obtenues après correction d'une translation globale des images par interpolation des images, selon l'équation suivante des images corrigées f et k : k(x) = g(x + w) 5 j'(X) = f(x + w) où x représente les deux coordonnées x et y des pixels des images. Cette translation w peut être obtenue sur la base d'un étalonnage du détecteur avec une éprouvette standard et d'une direction connue de balayage par rapport à l'écran 110 et la taille physique d'un pixel de l'écran 110. Ainsi, avant de lancer le calcul de 1 0 corrélation intégrée, on peut corriger au préalable la dérive des images due au mouvement du point d'émission, c'est-à-dire le mouvement de balayage du faisceau d'électrons. Suivant un mode de réalisation de l'invention, les première et deuxième images f et g sont obtenues après correction par centrage au centre 0 de projection. 15 Par exemple, l'éprouvette (matériau CR) peut être inclinée à 70° en position d'acquisition HR-EBSD et l'écran 110 peut également être désorienté de quelques degrés aussi. Si l'acquisition est réalisée sur une grande zone, la distance z* entre le point S d'émission et l'écran peut varier. La valeur de z* a un effet sur le grandissement de l'image de diffraction. Cet effet peut être corrigé en ajustant la 2 0 taille de l'image de diffraction centrée au centre de projection 0 = (x*, y*)T avant de lancer la corrélation d'images, surtout lorsque la zone mesurée dépasse 100 pm, selon l'équation suivante des images corrigées "f et 'g : g(x- x*) g [(zg*/zf*)(x - x*)] "f(x- x*) [(zg*/zf*)(x - x*)] 25 où zg* est la distance entre le point d'émission d'image « déformée », et zf* est celle entre le point d'émission d'image de référence. En conséquence, l'effet des paramètres P de projection sur le champ Ti de sensibilité peut être traité et corrigé e au préalable et la résolution être effectué sur le tenseur r . 33 On décrit ci-dessous à titre d'exemple non limitatif un essai sur une éprouvette de traction polycristalline, en référence aux figures 5A, 5B, 5C, 6A à 61. On a réalisé la première image f à l'état de référence et la deuxième image g à l'état déformé d'une éprouvette. Pour réaliser l'état déformé, l'éprouvette d'acier 5 inoxydable de grade AISI 316L à gros grains a été polie et soumise à un chargement de traction à l'intérieur de la chambre d'un microscope électronique à balayage à l'aide d'une platine d'essai in-situ. La direction de chargement est horizontale dans les figures 5A et 5B. Lors d'une première phase élastique du chargement, une acquisition HR-EBSD a été faite en considérant une zone d'intérêt focalisée sur un 10 point triple de la microstructure, c'est-à-dire séparant trois grains Gl, G2 et G3. L'algorithme suivant l'invention a été utilisé pour traiter les données. La figure 5A montre dans une zone d'intérêt de l'éprouvette une image en électrons secondaires prise au microscope électronique à balayage. La figure 5B montre dans la zone d'intérêt de l'éprouvette de la figure 5A une figure de pôle 15 inverse obtenue par analyse EBSD standard. La figure 6A montre l'image f de référence obtenue en un point de la zone d'intérêt, en fonction des coordonnées dans le plan du détecteur, x en abscisse et y en ordonnées. La figure 6B montre la deuxième image g de l'état déformé, en fonction des coordonnées x en abscisse et des coordonnées y en ordonnées. La figure 6C montre l'écart initial ECI entre 20 l'image de la figure 6B et l'image de la figure 6A, en fonction des coordonnées x en e abscisse et des coordonnées y en ordonnées. La figure 6D montre le tenseur r gradient de déformation élastique, obtenu par la méthode ADDICTED suivant l'invention à partir des images f et g des figures 6A et 6B. La figure 6E montre selon l'échelle ECH de niveaux de gris représentée à droite la valeur numérique de la 25 composante u, du champ de déplacement ayant été calculé à partir des images des figures 6A et 6B par la méthode ADDICTED suivant l'invention, en fonction des coordonnées x en abscisse et des coordonnées y en ordonnées. La figure 6F montre selon l'échelle ECH de niveaux de gris (valeurs numériques de pixels) représentée à droite la valeur numérique de l'autre composante uy du champ de déplacement ayant 30 été calculé à partir des images des figures 6A et 6B par la méthode ADDICTED suivant l'invention, en fonction des coordonnées x en abscisse et des coordonnées y 34 en ordonnées. La figure 6G montre l'image déformée et corrigée gu , c'est-à-dire ayant été corrigée à partir de la figure 6B par le champ (ux,uy) de déplacement des figures 6E et 6F, ayant été calculé à partir des images des figures 6A et 6B par la méthode ADDICTED suivant l'invention, en fonction des coordonnées x en abscisse 5 et des coordonnées y en ordonnées. La figure 6H montre le champ de résidu r, égal à l'écart entre l'image déformée corrigée gu de la figure 6G, obtenue par la méthode ADDICTED suivant l'invention, et l'image f de la figure 6A. La figure 61 montre une image de diffraction du grain G3, où les signaux contenus dans l'image 6H sont visibles. 1 0 Il existe une rotation peu perceptible à l'oeil nu entre les deux images des figures 6A et 6B mais notable lorsque l'on calcule l'écart initial entre celles-ci, comme illustré à la figure 6C. On peut conclure des figures 6H et 6C que l'écart initial dominant illustré sur la figure 6C disparaît grâce à la méthode ADDICTED suivant l'invention à la figure 6H et qu'un champ de bandes "fantômes" devient 15 visible dans les résidus de la figure 6H. On observe que l'image de diffraction de la figure 61 présente une similitude des bandes Kikuchi avec les bandes "fantômes" présentes dans le résidu illustré sur la figure 6H. Ce phénomène indique que "l'image de fond" n'a pas été acquise de façon suffisamment précise et robuste et qu'en conséquence le signal de diffraction du 2 0 grain G3 est encore contenu dans cette image. En interprétant le champ de résidu fourni à la figure 6H par la méthode ADDICTED suivant l'invention, il est possible de révéler des informations auparavant cachées et qui pourraient être utiles pour construire une meilleure image de fond. Pour initialiser la toute première estimation {r17^e } lors de l'étape E6 (la 25 première image g de diffraction « déformée »), la matrice identité est un bon choix pour le tenseur Fe , du fait des faibles déformations attendues, et de l'absence d'informations supplémentaires. Pour initialiser lors de l'étape E6 le tenseur sur les autres pixels, il est généralement préférable de prendre la valeur de {Fe} ayant été calculée pour l'un de ses voisins, car le niveau de déformation est supposé être proche entre des points matériels voisins. La vitesse de convergence en fonction des conditions de calcul est résumée dans le Tableau 1 ci-dessous : 5 Tableau 1 : Nombre d'itérations pour différentes conditions de calcul Lissage gaussien Noyau de 2 pixels Noyau de 1 pixel Non lissé Initialisé par tenseur 17 65 700 -îl e r = matrice identité Initialisé par tenseur 14 50 300 F des pixels voisins On peut conclure que le lissage gaussien de l'image de diffraction réduit drastiquement le nombre d'itérations nécessaires pour converger et le lissage avec un noyau plus grand (2 pixels) rend la convergence plus rapide qu'avec un noyau plus 1 0 petit (1 pixel). L'initialisation avec l'estimation du tenseur F e des pixels voisins accélère également la convergence surtout pour les calculs sur l'image de diffraction non filtrée par le lissage. En termes de valeurs du tenseur Fe obtenues, différentes conditions de calcul ne conduisent pas à des valeurs du tenseur Fe identiques mais assez proches. L'écart type des valeurs du tenseur F e obtenues est inférieur à 10-4. 15 Les calculs présentés ci-après ont été effectués à partir de l'image g du matériau sous contrainte de la figure 5A, cette image est lissée avec un noyau e gaussien de 2 pixels et avec initialisation par la valeur du tenseur r du pixel e voisin. L'hypothèse de contrainte plane a été appliquée au tenseur F, qui a donc été transformé en (assimilé au) tenseur Fe (partie élastique Fe du tenseur gradient de 2 0 déformation F). La zone d'émission (interagissant avec le faisceau électronique) est de profondeur d'environ 25 nm, qui rend l'hypothèse de contraintes planes légitime. 36 Il faut noter que cette hypothèse de contrainte plane repose sur l'hypothèse que la normale à la surface soit bien connue, ce qui n'est en général pas le cas après déformation. Un relief sur la surface va dégrader la précision des mesures HREBSD. L'acquisition HR-EBSD a été réalisée à l'entrée en plasticité et il est alors 5 légitime de négliger l'effet du relief On suppose ici que la rotation peut être grande, mais que la déformation pure est faible (régime élastique). Il convient alors de garder un formalisme de grandes transformations sur Fe, et calculer la rotation R par une décomposition polaire, ce qui permet d'estimer le tenseur de déformation pure élastique de Green-Lagrange droit Ue lié aux contraintes selon l'équation suivante : 10 Fe=Rue Concernant les déformations pures, une hypothèse de petites perturbations est légitime, I - U e REVENDICATIONS1. Procédé de traitement d'images, obtenues par un détecteur à diffraction, d'un matériau cristallin ou polycristallin, dans lequel on mesure par le détecteur : - une première image (n de diffraction du matériau se trouvant dans un état de référence, donnant des premières valeurs numériques (4x,y)) de pixels en fonction de deux coordonnées (x,y) de pixels, - au moins une deuxième image (g) de diffraction du matériau se trouvant dans un état déformé par rapport à l'état de référence, la deuxième image donnant 10 des deuxièmes valeurs numériques (g(x,y)) de pixels en fonction des coordonnées (x,y) de pixels, caractérisé en ce que dans un calculateur, un champ (ux,uy) de déplacement, pour déplacer des pixels de la première image à des pixels d'une image déformée, en fonction : 15 - des deux coordonnées (x,y) des pixels, - de coordonnées prédéterminées (x*, y*) d'un centre (0), correspondant à la projection normale, dans le plan d'image du détecteur, d'un point source (S) du faisceau diffracté dans le matériau, et ( 2 e fr3e fr4e, fr5e .fr6e, Êle frse - de composantes fi ) d'un tenseur Ê e 20 gradient de déformation élastique, est préenregistré dans une mémoire, au cours d'une première étape (E6, E12) de calcul, on fait prendre au tenseur F gradient de déformation élastique une valeur déterminée ( e'el) du tenseur gradient de déformation élastique, 25 au cours d'une deuxième étape (E7) de calcul, on calcule le champ (ux,uy) de déplacement courant à partir du tenseur (Fe) gradient de déformation élastique courant et des coordonnées (x,y) de pixels de la première image (/), au cours d'une troisième étape (E8), on calcule des troisièmes valeurs numériques de pixel d'une image déformée Ciu(x) = g(x + u(x)) ) en 41 corrigeant la deuxième image (g) aux coordonnées de pixels auxquelles a été additionné le champ de déplacement (ux,uy) courant, au cours d'un algorithme itératif, on effectue des itérations des première, deuxième et troisième étapes (E12, E7, E8) de calcul sur des valeurs déterminées de 5 tenseur modifiées ( {ft' e} , jusqu'à remplir un critère de convergence sur la valeur déterminée ( e }) du tenseur gradient de déformation élastique, pour calculer le champ de déplacement (ux,uy) correspondant. 2. Procédé de traitement d'images suivant la revendication 1, caractérisé en ce que le tenseur Êe gradient de déformation élastique est égal à Fle F2e F3. fre4 fr5e fre 6 fr7e p8e 1 Fie Fe Fe Fe Fe F6. F7. F8e où sont les composantes du tenseur ee gradient de déformation élastique. 3. Procédé de traitement d'images suivant l'une quelconque des revendications précédentes, caractérisé en ce que le champ ux,uy de déplacement, 15 pour déplacer des pixels de la première image à des pixels d'une image déformée, en fonction : - des deux coordonnées x,y de pixels, - des coordonnées prédéterminées x*,y*, z* du centre (0), correspondant à la projection normale, dans le plan d'image du détecteur, du point source (S) du 20 faisceau diffracté dans le matériau, et ie F2e F4e Fse F6e F7e - des composantes F -e8 du tenseur gradient de déformation élastique, est égal à Êe 42 ux (x, y) z* (frie (X - X*)-1- fr2e (y - y*) + P3ez*) = (x - x*) P7e (x - x*) + fr8e (y - y*) + z * z (P4e (X - X*) + te (y - y*) + P (x, y)6ez*) = fr (y _ y*) e 4. Procédé de traitement d'images suivant l'une quelconque des revendications précédentes, caractérisé en ce que l'algorithme itératif est effectué par une méthode de type Gauss-Newton. 5 5. Procédé de traitement d'images suivant l'une quelconque des revendications précédentes, caractérisé en ce qu'au cours de l'algorithme itératif, on calcule un vecteur Wel de correction vérifiant l'équation [m]We 1= {7} où [NI] est une matrice hessienne de dimension 8 x 8, ayant comme 10 coefficients Ai-(n-1) = E (vf(x) (D, (x, e f (X) - 41 (X , e 1)) `r"- 41 ROI où f représente les premières valeurs de pixels de la première image, x représente les deux coordonnées x et y de pixels, Of (X; {Fe est un champ de sensibilité du champ (ux,uy) de déplacement par rapport à la composante frie du tenseur t'e gradient de déformation élastique et est égal à la dérivée partielle du champ (ux,uy) de déplacement par rapport à chaque composante frie du tenseur Êe gradient de déformation élastique, {y} est le résidu ayant comme coefficients E (f(x) - n (X , tfr e })» f (X) - erre iX e ROI 43 (n) (x, Ifre (x) = g (x + u(x)) est la troisième valeur numérique de chaque pixel de l'image déformée et g représente les deuxièmes valeurs de pixels de la deuxième image, le critère de convergence sur la valeur déterminée Fe de tenseur étant que la 5 nonne du vecteur tee de correction soit inférieure à une borne positive prescrite (8), non nulle, la valeur déterminée IF } de tenseur étant incrémentée du vecteur 'tee de correction à chaque itération de la première étape (E12) de calcul. 6. Procédé de traitement d'images suivant la revendication 5, caractérisé en 10 ce que le tenseur Êe gradient de déformation élastique possède, comme fie, je F 3e , F F 5e , F 6e , F; Fse et une composantes Fie , huit composantes neuvième composante fixée à 1, à savoir selon l'équation suivante : Fie fr2e F3e Wel= File És se p6e F; fr8e 1 le champ de sensibilité du champ u(x) de déplacement a comme 15 composantes (Di e )) les composantes selon les équations suivantes : z * (x - x*) cle x2 = fr7e (X - X*) ftie - y*) + z * (z*)2 ^Fie (X x*) P8e (Y Y*)+ z * ) x4 = 0 CI) xl = ^ F7e (X - P8e (y - y*) + z * z * (Y - Y*) (1) x3 = 44 z * (x - x*)(Pe (x - x*) + .Ê2e (y - y*) + Ê3ez*) (Dx7 (fr; (X - X*) fr8e (y - y*) + z *)2 5 (1) x8 - z * (Y Y*)(-frie (x x*) fie (y - Y*) kz*) (P.; (x - x*) + fr; (y - y*) + z *)2 (I) =0 Yi ei y2 = 0 CD y3 =0 'o (13 y8 = (fre(X_X*)+fre(y_y*)+Z*)2 7. Procédé de traitement d'images suivant l'une quelconque des 15 revendications précédentes, caractérisé en ce qu'il est exécuté pour la majeure partie ou l'ensemble des pixels de la première image (f) et de la deuxième image (g). z * (x - x*) (byel - fr7e (X - Xee) Pse (Y - Y*) z (Dy' p7e (x - x*2.; (y +/-8- .Y*) e (y- y*) + z * (z*)2 (13 y6 - fr.; - X*) ftie (y - y*) + z * z * (x - x*)(fr: (x - x*) + P5e (y - y*) + Ê6ez*) (1) Y7 (fr; (X - X*) fr8e - y*) z *)2 Z * (y - y*)(fr: (x - x*) + te (y - y*) + frjz*) 45 8. Procédé de traitement d'images suivant l'une quelconque des revendications précédentes, caractérisé en ce que les première et deuxième images sont obtenues après filtrage (E4) de valeurs de pixels surexposés en les remplaçant par une moyenne de pixels voisins de ceux-ci. 5 9. Procédé de traitement d'images suivant l'une quelconque des revendications précédentes, caractérisé en ce que les première et deuxième images sont obtenues après filtrage (E4) des valeurs de pixels par soustraction de tendances globales de niveau de gris, représentées par un polynôme d'ordre 2 ou 3 obtenu par un processus de régression. 10 10. Procédé suivant l'une quelconque des revendications précédentes, caractérisé en ce que les première et deuxième images sont obtenues après filtrage (E5) des valeurs de pixels par un filtre de lissage gaussien. 11. Procédé de traitement d'images suivant l'une quelconque des revendications précédentes, caractérisé en ce qu'il est exécuté pour plusieurs 15 deuxièmes images (g). 12. Procédé de traitement d'images suivant l'une quelconque des revendications précédentes, caractérisé en ce que l'on fournit (E13) sur une sortie au moins l'un parmi: - le tenseur (Fe) gradient de déformation élastique, calculé comme remplissant le critère de convergence, - un résidu (r), calculé comme étant la différence entre d'une part les troisièmes valeurs numériques de pixels de l'image déformée corrigée ( go ), ayant été calculées pour le champ de déplacement (u',uy) correspondant comme remplissant le critère de convergence sur la valeur déterminée ( fre} ) du tenseur gradient de déformation élastique, et d'autre part les premières valeurs numériques (gx,y)) de pixels de la première image (f), - une moyenne quadratique des résidus (r ) calculés pour respectivement plusieurs deuxièmes images (g), 46 - le champ de déplacement (u.'uy) correspondant, calculé comme remplissant le critère de convergence sur la valeur déterminée ( \fre} ) du tenseur gradient de déformation élastique, - une image déformée corrigée ( gu ), correspondant aux troisièmes valeurs 5 numériques de pixels ayant été calculées à partir du champ de déplacement (ux,uy) correspondant, calculé comme remplissant le critère de convergence. 13. Procédé de traitement d'images suivant l'une quelconque des revendications précédentes, caractérisé en ce que les première et deuxième images sont obtenues (E3) après soustraction d'un fond d'image. 10 14. Dispositif de traitement d'images d'un matériau cristallin ou polycristallin, comportant un détecteur (10) à diffraction permettant d'acquérir : - une première image (f) de diffraction du matériau se trouvant dans un état de référence, donnant des premières valeurs numériques (ftx,y)) de pixels en fonction de deux coordonnées (x,y) de pixels, 15 - au moins une deuxième image (g) de diffraction du matériau se trouvant dans un état déformé par rapport à l'état de référence, la deuxième image donnant des deuxièmes valeurs numériques (g(x,y)) de pixels en fonction des coordonnées (x,y) de pixels, caractérisé en ce que le dispositif comporte au moins un calculateur (12), 20 comportant au moins une mémoire, dans laquelle est enregistrée un champ (ux,uy) de déplacement, pour déplacer des pixels de la première image à des pixels d'une image déformée, en fonction : - des deux coordonnées (x,y) de pixels, - de coordonnées prédéterminées (x*,y*,z*) d'un centre (0), correspondant à 25 la projection normale, dans le plan d'image du détecteur, d'un point source (S) du faisceau diffracté dans le matériau, et (Fie ,F2e ,F3 e ,F4e.,F5e ,F6e ,F7 e - de composantes ) d'un tenseur ( t''e ) gradient de déformation élastique, le calculateur (12) étant configuré pour : 47 au cours d'une première étape (E6, E12) de calcul, faire prendre au tenseur (F ) gradient de déformation élastique courant une valeur déterminée ( fre5), au cours d'une deuxième étape (E7) de calcul, calculer le champ (ux,uy) de déplacement courant à partir du tenseur (Fe) gradient de déformation élastique 5 courant et des coordonnées (x,y) pour chacun des pixels de la première image (f), au cours d'une troisième étape (E8) de calcul, calculer des troisièmes valeurs numériques de pixels d'une image déformée (g1,(x) = g(x + u(x)) ) en appliquant la deuxième image (g) aux coordonnées de pixels auxquelles a été additionné le champ de déplacement (u.'uy) courant, 10 au cours d'un algorithme itératif, effectuer des itérations des première, deuxième et troisième étapes (E12, E7, E8) de calcul sur des valeurs déterminées de tenseur réactualisées ( , jusqu'à remplir un critère de convergence sur la valeur déterminée (5Fe) du tenseur gradient de déformation élastique, pour calculer le champ de déplacement (ux,uy) correspondant. 15 15. Programme d'ordinateur, comportant des instructions de code pour la mise en oeuvre du procédé de traitement d'images d'un matériau cristallin ou polycristallin suivant l'une quelconque des revendications 1 à 13, lorsqu'il est exécuté sur un calculateur.