FR 2983317 A1 20130531 FR 1160915 A 20111129 B11352 - DD13119JBD 1 PROCÉDÉ DE FORMATION D'UN HOLOGRAMME SYNTHÉTIQUE DANS UNE IMAGE TRAMÉE Domaine de l'invention La présente invention concerne un procédé d'intégration d'un hologramme synthétique dans une image visible. Plus particulièrement, la présente invention concerne un procédé d'intégration d'un hologramme synthétique par codage d'ouvertures dans une image tramée dont le rendu visuel est en niveaux de gris. Exposé de l'art antérieur Dans de nombreux secteurs, notamment dans l'industrie du luxe (par exemple la parfumerie, la bijouterie ou la maroquinerie), ou dans le domaine de médicament, la lutte contre la copie de produits de marque est une préoccupation quotidienne. De nombreux procédés et dispositifs sont actuellement utilisés pour tenter de garantir l'authenticité des produits de marque. Parmi ces techniques, il a été proposé de reporter, sur les produits à identifier, des puces transparentes sur lesquelles sont formés des hologrammes. De tels hologrammes sont obtenus à partir d'une image visible, retrouvée à la lecture de l'hologramme. La présence de l'hologramme est difficilement détectable à l'oeil nu, et sa lecture directe sans dispositif de B11352 - DD13119JBD 2 lecture adapté ne permet pas d'identifier l'image visible source de l'hologramme. La figure 1 est un organigramme illustrant un procédé de formation d'un hologramme. La figure 2 illustre un exemple de rendu obtenu à l'aide du procédé de formation d'hologramme de la figure 1 dans le cas d'un hologramme synthétique par codage d'ouvertures. En figure 1, on part d'une image source de l'hologramme 10 (IMAGE H) dont on calcule, à une étape 12, la transformée de Fourier. Le calcul de cette transformée de Fourier permet d'obtenir une image d'amplitude de la transformée de Fourier 14 (AH) et une image de phase de la transformée de Fourier 16 (911). Les images d'amplitude 14 et de phase 16 sont ensuite combinées dans une nouvelle image qui correspond à l'hologramme fini (étape 18, HOLOGRAM). Différentes combinaisons des images d'amplitude 14 et de phase 16 sont possibles, l'idée étant que l'éclairement de l'hologramme avec un faisceau adapté, combiné à une optique réalisant une transformée de Fourier inverse, permet de lire l'image source de l'hologramme 10 (IMAGE H). L'hologramme 18 est constitué de nombreux pixels dont la forme est définie par les pixels correspondants des images obtenues par le calcul de la transformée de Fourier 12. L'hologramme 18 présente ainsi une résolution identique à celle de l'image initiale 10. La figure 2 illustre un exemple d'hologramme obtenu par le procédé de la figure 1, cet hologramme étant un hologramme synthétique par codage d'ouvertures, encore appelé hologramme synthétique à détour de phase. En figure 2 est illustrée une portion d'hologramme constituée de 5x5 pixels. Un repère (x, y) est utilisé pour identifier chaque pixel. Chaque pixel 22 comprend une zone sombre centrale 24, en forme d'ellipse, le fond du pixel étant clair. Le centre de la zone 24 est aligné, dans chaque pixel 22, sur le centre du pixel selon l'axe y. B11352 - DD13119JBD 3 La taille de la zone sombre 24 dans chaque pixel dépend du résultat du calcul de l'image d'amplitude 14 de la transformée de Fourier 12 de l'image source IMAGE H pour le pixel de mêmes coordonnées. Des zones 24 de tailles plus ou moins importantes sont donc prévues dans l'hologramme. En outre, le décalage du centre de chaque zone 24 par rapport au centre du pixel selon l'axe x dépend de la valeur du pixel de mêmes coordonnées de l'image de phase 16 (911). On notera que différents procédés de définition d'hologramme par codage d'ouvertures sont connus. Par exemple, la forme des zones 24 dans les pixels 22 peut différer de la forme elliptique proposée ici. Notamment, il a été proposé de former des zones 24 rectangulaires, de taille plus ou moins importante en fonction de l'image d'amplitude 14, et décalées plus ou moins du centre du pixel en fonction de l'image de phase 16 (selon l'axe x). La figure 3 illustre un dispositif de lecture d'un hologramme synthétique tel que celui de la figure 2, obtenu par un procédé tel que celui de la figure 1. En figure 3, le dispositif de lecture est un dispositif en réflexion. On notera qu'il existe des dispositifs de lecture similaires en transmission. Un faisceau lumineux 30 est transmis, par exemple par un cube séparateur de faisceau 32, en direction d'un hologramme 34 (non représenté en détail) formé sur une puce 36. Les pixels de l'hologramme 34 réfléchissent le faisceau lumineux. Le faisceau réfléchi par l'hologramme repasse dans le cube séparateur de faisceau 32 puis est combiné par une lentille 38 dont la fonction est de réaliser une transformée de Fourier inverse des rayons reçus de façon à afficher une image sur un dispositif d'affichage (non représenté) placé au point focal de la lentille. On notera que, si un codage par clef de phase de l'hologramme est prévu à la formation de l'hologramme, une lame 35 de phase similaire à la clef de phase doit être placée dans le B11352 - DD13119JBD 4 dispositif de lecture de la figure 3 pour réaliser la lecture de l'hologramme. Des techniques d'intégration d'hologrammes dans une image ont été proposées. La publication de Joseph Rosen et Bahram Javidi intitulée "Hidden images in halftone pictures", Appl. Opt. 40, No. 20, 3346-3353, présente l'une de ces techniques. Les figures 4A et 4B illustrent, respectivement, une première image 40 dont le rendu visuel est un rendu en niveaux de gris et une deuxième image 42, issue de la première image, dans laquelle un hologramme synthétique par codage d'ouvertures est intégré. La deuxième image est obtenue par le procédé décrit dans la publication ci-dessus. La figure 4A est une image obtenue par tramage d'une image en niveaux de gris. Le procédé de tramage, couramment utilisé par exemple dans les procédés d'impression, consiste à associer, à chaque teinte d'un pixel initial en niveaux de gris, un pixel comprenant une zone centrale blanche dans un fond noir (ou l'inverse), la zone centrale étant plus ou moins grande en fonction de la teinte du pixel équivalent de l'image en niveaux de gris. Visuellement, et si l'on s'éloigne de l'image, l'image tramée et l'image en niveaux de gris apparaissent identiques. On notera que les techniques de tramage sont généralement réalisées avec une résolution assez importante pour 25 que, à l'oeil nu, la trame n'apparaisse pas. L'image 42 de la figure 4B est obtenue à partir de l'image 40 de la figure 4A et intègre un hologramme synthétique. Dans l'image 42, l'hologramme est intégré en modifiant la taille et l'alignement des zones centrales des différents pixels, en 30 fonction de l'hologramme à intégrer. Une lecture de l'image de la figure 4B, par exemple à l'aide du dispositif de la figure 3, associé par exemple à un filtre, permet d'obtenir l'image source utilisée pour former l'hologramme. B11352 - DD13119JBD Des inconvénients de la technique mentionnée ci-dessus apparaissent en figure 4B : en effet, on voit dans cette figure que certaines zones claires décalées se superposent à d'autres zones claires de pixels voisins. Cette superposition implique 5 une perte d'information importante, qui n'est pas désirée. La lecture de l'hologramme codé se trouve donc particulièrement dégradée par cette intégration. En outre, la technique ci-dessus implique un mauvais contraste en niveaux de gris et des contraintes importantes au niveau du codage de phase. Résumé Un objet d'un mode de réalisation de la présente invention est de prévoir un procédé palliant tout ou partie des inconvénients des techniques d'intégration d'hologrammes dans 15 des images de l'art antérieur. Plus particulièrement, un objet d'un mode de réalisation de la présente invention est de prévoir un procédé permettant d'intégrer un hologramme dans une image en niveaux de gris dont la lecture directe est de bonne qualité. 20 Ainsi, un mode de réalisation de la présente invention prévoit un procédé d'intégration d'un hologramme synthétique dans une image directement observable d'une scène, comprenant les étapes suivantes : définir une première image de la scène, comprenant des pixels en niveaux de gris ; fabriquer des 25 première et deuxième matrices à partir de la première image, la première matrice comprenant des pixels d'une première et d'une deuxième teinte selon que le pixel correspondant de la première image a un niveau de gris inférieur ou supérieur à un seuil, chaque élément la deuxième matrice comprenant une valeur égale à 30 l'écart de niveau de gris entre le pixel correspondant de la première image et le pixel correspondant de la première matrice ; fabriquer une troisième matrice de pixels en calculant l'image de phase de la transformée de Fourier d'une image issue d'une deuxième image source de l'hologramme ; fabriquer une 35 quatrième matrice de pixels, chaque pixel de la quatrième B11352 - DD13119JBD 6 matrice comprenant une zone centrale dont la surface est déterminée à partir de l'élément correspondant de la deuxième matrice et dont un décalage par rapport au centre du pixel est issu du pixel correspondant de la troisième matrice ; et réaliser une lithographie d'une couche opaque en surface d'une plaquette selon le motif défini par la quatrième matrice de pixels, la lithographie étant prévue pour éliminer, dans un pixel, la couche opaque à l'intérieur de la zone centrale si le pixel équivalent de la première matrice est de la première teinte et pour éliminer la couche opaque à l'extérieur de la zone centrale si le pixel équivalent de la première matrice est de la deuxième teinte. Selon un mode de réalisation de la présente invention, la taille de la zone centrale de chacun des pixels de la quatrième matrice est calculée à partir de la valeur du pixel correspondant d'une cinquième matrice de pixels, chaque pixel de la cinquième matrice étant obtenu par calcul de l'amplitude diffractée d'une ouverture définie dans un pixel, le rapport entre la surface de l'ouverture et la surface du pixel étant égal à la valeur de l'élément correspondant de la deuxième matrice. Selon un mode de réalisation de la présente invention, l'amplitude Ah d'un pixel de la cinquième matrice (IMAGE GA) est liée à la valeur I de l'élément correspondant de la deuxième matrice par la relation suivante : IIAL = I'' .4 (Jignm ) J1 étant la première fonction de Bessel. Selon un mode de réalisation de la présente invention, la troisième matrice est obtenue en réalisant les étapes suivantes : (a) calculer une sixième matrice d'amplitude obtenue à partir de la deuxième image source de l'hologramme en réalisant les mêmes étapes que celles permettant d'obtenir la cinquième matrice à partir de la deuxième matrice ; (b) combiner les pixels de la sixième matrice d'amplitude à une distribution de phase aléatoire ou à une composante continue pour obtenir une B11352 - DD13119JBD 7 valeur complexe ; et (c) calculer l'image de phase de la transformée de Fourier de la valeur complexe. Selon un mode de réalisation de la présente invention, la troisième matrice est obtenue par un procédé d'optimisation 5 tirant profit de la cinquième matrice. Selon un mode de réalisation de la présente invention, le procédé d'optimisation consiste à : (a) calculer une sixième matrice d'amplitude obtenue à partir de la deuxième image source de l'hologramme en réalisant les mêmes étapes que celles 10 permettant d'obtenir la cinquième matrice à partir de la deuxième matrice ; (b) combiner les pixels de la sixième matrice d'amplitude à une distribution de phase aléatoire ou à une composante de phase continue pour obtenir une première valeur complexe ; (c) calculer la transformée de Fourier de la première 15 valeur complexe ; (d) combiner les pixels de la cinquième matrice de pixels et les pixels de l'image de phase issue de la transformée de Fourier de l'étape (c) pour obtenir une deuxième valeur complexe ; (e) calculer la transformée de Fourier inverse de la deuxième valeur complexe ; (f) combiner les pixels de la 20 sixième matrice d'amplitude aux pixels de l'image de phase issue du calcul de l'étape (e) pour obtenir une nouvelle première valeur complexe ; (g) réitérer l'étape (c) pour calculer la transformée de Fourier de la nouvelle première valeur complexe, la troisième matrice correspondant à l'image de phase obtenue de 25 cette dernière transformée de Fourier. Selon un mode de réalisation de la présente invention, le procédé d'optimisation comprend en outre la réitération des étapes (d) à (g), la troisième matrice correspondant à l'image de phase obtenue à l'étape (g) après plusieurs cycles de 30 réitération des étapes (d) à (g). Selon un mode de réalisation de la présente invention, au moins trois cycles sont réalisés. Selon un mode de réalisation de la présente invention, le décalage de la zone centrale par rapport au centre du pixel 35 de chaque pixel de la quatrième matrice est directement B11352 - DD13119JBD 8 proportionnel de la valeur du pixel correspondant de la troisième matrice. Selon un mode de réalisation de la présente invention, le seuil de niveau de gris est égal à 0,5, à 10 % près, sur une 5 échelle normalisée de niveaux de gris. Selon un mode de réalisation de la présente invention, la zone centrale de chaque pixel de la quatrième matrice est elliptique. Selon un mode de réalisation de la présente invention, 10 le rapport de forme de la zone centrale elliptique de chaque pixel de la quatrième matrice est égal à 1-1/2. Selon un mode de réalisation de la présente invention, la première teinte de la première matrice est une teinte minimum, la deuxième teinte de la première matrice est une 15 teinte maximum, la lithographie étant prévue pour éliminer la couche opaque à l'intérieur de la zone centrale d'un pixel si le pixel équivalent de la première matrice est de la première teinte, et pour éliminer la couche opaque à l'extérieur de la zone centrale d'un pixel si le pixel équivalent de la première 20 matrice est de la deuxième teinte. Selon un mode de réalisation de la présente invention, une taille minimum de la zone centrale des pixels de la quatrième matrice est imposée. Selon un mode de réalisation de la présente invention, 25 la taille minimum de la zone centrale des pixels de la quatrième matrice correspond à un niveau de gris de 0,02 pour les tons sombres des pixels et à 0,98 pour les tons clairs des pixels, sur une échelle normalisée de niveaux de gris. Brève description des dessins 30 Ces objets, caractéristiques et avantages, ainsi que d'autres seront exposés en détail dans la description suivante de modes de réalisation particuliers faite à titre non limitatif en relation avec les figures jointes parmi lesquelles : B11352 - DD13119JBD 9 la figure 1, précédemment décrite, est un organigramme illustrant un procédé de formation d'hologrammes synthétiques par codage d'ouvertures connu ; la figure 2 illustre un exemple de résultat obtenu à 5 l'aide du procédé de la figure 1 ; la figure 3 illustre un dispositif de lecture d'un hologramme tel que celui de la figure 2, obtenu par exemple par un procédé tel que celui de la figure 1 ; les figures 4A et 4B, précédemment décrites, 10 illustrent respectivement une image tramée et une image issue de l'image de la figure 4A dans laquelle un hologramme est intégré ; les figures aA à 5C illustrent une technique de formation de pixels par codage d'ouvertures équivalents en 15 termes de diffraction ; la figure 6 est un organigramme illustrant un procédé selon un mode de réalisation de la présente invention ; la figure 7 illustre des résultats d'étapes du procédé de la figure 6 ; 20 la figure 8 illustre une courbe de la valeur d'éléments d'une matrice issue d'une image initiale en niveaux de gris obtenue selon une étape d'un mode de réalisation de la présente invention ; la figure 9 est une vue en perspective d'un pixel d'un 25 hologramme selon un mode de réalisation de la présente invention ; la figure 10 est une courbe de l'amplitude de la trame holographique en fonction du niveau de gris d'un pixel tramé ; la figure 11 est un organigramme illustrant un procédé 30 selon une variante d'un mode de réalisation de la présente invention ; les figures 12A à 12D illustrent des résultats de lecture obtenus à partir d'hologrammes intégrés dans des images visibles formées par le procédé des figures 6 et 11, en fonction 35 du nombre d'itérations de ce procédé ; B11352 - DD13119JBD 10 la figure 13 illustre un résultat du procédé selon un mode de réalisation de la présente invention ; les figures 14A à 14D illustrent des résultats obtenus en réalisant des variantes d'un procédé selon un mode de 5 réalisation de la présente invention ; la figure 15 illustre une variante d'un procédé selon un mode de réalisation de la présente invention ; les figures 16A à 16C illustrent des résultats obtenus en réalisant des variantes d'un procédé selon un mode de 10 réalisation de la présente invention ; et les figures 17A à 17D, 18 et 19 illustrent des applications d'un procédé selon un mode de réalisation de la présente invention. Par souci de clarté, de mêmes éléments ont été 15 désignés par de mêmes références aux différentes figures. Description détaillée Les figures aA à 5C illustrent une propriété d'équivalence de pixels d'un hologramme synthétique par codage d'ouvertures dont on prévoit de tirer profit pour intégrer un 20 hologramme dans une image en niveaux de gris. Ces trois figures montrent, en trois dimensions, une portion d'une plaquette transparente 50 dont la surface définit un pixel d'un hologramme synthétique par codage d'ouvertures. Le pixel représenté présente une largeur "a". 25 Dans l'exemple de la figure aA, le pixel considéré comprend une zone opaque centrale 52, de forme elliptique, sur un fond transparent. En pratique, la zone opaque peut être formée sur la plaquette 50 en déposant, sur l'ensemble de cette plaquette 50, une couche opaque (absorbante ou réfléchissante), 30 par exemple en oxyde de platine (PtOx), que l'on vient graver par exemple par lithographie de façon à définir la zone centrale 52. Dans l'exemple représenté en figure aA, le centre de la zone opaque 52 est décalé du centre du pixel d'une distance b 35 vers la gauche (selon l'axe x). Comme nous l'avons vu B11352 - DD13119JBD 11 précédemment, le décalage b est lié au pixel correspondant de l'image de phase de la transformée de Fourier de l'image source de l'hologramme. La taille de la zone opaque 52, quant à elle, est image du pixel correspondant de l'image d'amplitude de la transformée de Fourier de l'image source de l'hologramme. La forme elliptique proposée ici présente l'intérêt d'être, lorsqu'elle est utilisée pour former des pixels d'une image obtenue par tramage d'une image en niveaux de gris, relativement douce à l'oeil. En particulier, cette forme est beaucoup plus agréable à l'oeil que les formes rectangulaires des figures 4A et 4B. La figure 5B illustre un pixel complémentaire du pixel de la figure aA. Dans ce pixel, les régions opaques et transparentes du pixel sont inversées par rapport au pixel de la figure aA, c'est-à-dire que la zone 54, de même surface que la zone 52 de la figure aA, est transparente et le reste de la surface du pixel est recouvert d'une portion opaque 56. Le pixel de la figure 5C est obtenu à partir du pixel de la figure 5B et comprend un fond opaque 56 dans lequel est définie une zone transparente 58. La zone transparente 58 est de même taille que la zone 54. La zone 58 est décalée par rapport à la zone transparente 54, selon l'axe x vers la droite, de telle façon que le déphasage équivalent entre les zones 54 et 58 soit un déphasage de n (en relation avec l'image de phase de l'hologramme). En pratique, un tel déphasage peut être obtenu en décalant la zone 58 par rapport à la zone 54 d'une distance égale à a/2. Conformément aux principes des écrans complémentaires de Babinet, le pixel de la figure aA est équivalent, en ce qui concerne le motif de diffraction qu'il implique lorsqu'il est éclairé, au pixel complémentaire et déphasé de la figure 5C. Compte tenu de cette propriété, l'idée est de proposer un procédé d'intégration d'un hologramme dans une image tramée présentant un rendu en niveaux de gris, dans lequel les cellules les plus foncées de l'image en niveaux de gris sont constituées B11352 - DD13119JBD 12 non pas d'une zone centrale opaque de taille importante sur un fond transparent, mais d'une zone centrale transparente de petite taille sur un fond opaque (équivalente). La zone centrale est ensuite déphasée pour intégrer l'hologramme, comme on le verra par la suite. La figure 6 est un organigramme d'un procédé proposé pour intégrer un hologramme dans une image visible. Des étapes initiales 60 et 62 du procédé consistent à choisir, d'une part, une image réelle en niveaux de gris d'une scène, IMAGE G, que l'on souhaite voir à l'oeil nu une fois l'hologramme intégré dans l'image, et d'autre part une image source de l'hologramme, IMAGE H, que l'on souhaite intégrer sous forme d'hologramme synthétique par codage d'ouvertures dans une image visible finale (IMAGE F). A partir de l'image en niveaux de gris 60 (IMAGE G), on détermine une première matrice image (IMAGE GBE,w) 64 en noir et blanc, cette image comprenant des pixels d'une première teinte, par exemple noire, au niveau des pixels de l'image 60 dont le niveau de gris est inférieur à un seuil a (pixels foncés) et des pixels d'une deuxième teinte, par exemple blanche, pour les pixels de l'image 60 présentant un niveau de gris supérieur au seuil a (pixels clairs). On notera que l'inverse est également possible, l'idée étant de distinguer les pixels de l'image 60 dont le niveau de gris est inférieur ou supérieur au seuil a. A partir de l'image 60, des moyens de traitement et de calcul sont mis en oeuvre pour déterminer une matrice IMAGE G', 66, présentant une taille identique à la résolution de l'image 60, IMAGE G, dans laquelle chaque élément de la matrice contient une valeur proportionnelle à l'écart entre le niveau de gris du pixel correspondant de l'image 60 et le niveau de gris du pixel correspondant de l'image 64 (IMAGE GB00. Le lien entre les images 60 et 64 et la matrice 66 sera décrit plus en détail ci-après en relation avec la figure 8. Les valeurs stockées dans la matrice 66 sont ainsi croissantes, jusqu'à un seuil, avec le B11352 - DD13119JBD 13 niveau de gris des pixels correspondants de l'image initiale 60 puis décroissantes lorsque le niveau de gris des pixels correspondants de l'image 60 augmente au-dessus du seuil (l'inverse est également possible). En pratique, le seuil de détection de la matrice 66 est le même seuil que le seuil a de détection de l'image 64. Il est choisi de préférence à un niveau de gris moyen de l'image 60, par exemple de 0,5 à 10 % près, de préférence égal à 0,5, sur une échelle normalisée de niveaux de gris. Une étape suivante 68 consiste à obtenir, à partir de la matrice 66 (IMAGE G'), une matrice image IMAGE GA 68 d'amplitude de la trame holographique. La détermination des pixels de l'image 68 sera décrite ci-après plus en détail, en relation avec les figures 9 et 10. En parallèle, une étape 70 consiste à calculer la transformée de Fourier d'une image issue de l'image source de l'hologramme 62 (IMAGE H). Plus précisément, on calcule la transformée de Fourier d'une image d'amplitude, IMAGE HA, obtenue à partir de l'image source IMAGE H de manière similaire à l'obtention de l'image 68, IMAGE GA, à partir de l'image 66, IMAGE G', comme cela sera décrit plus en détail ci-après. La transformée de Fourier calculée à l'étape 70 permet d'obtenir une image d'amplitude AH et une image de phase 72 ((pH). L'image d'amplitude AH n'est pas utilisée. Pour améliorer la qualité de l'hologramme, une distribution de phase aléatoire peut être ajoutée à l'image d'amplitude IMAGE HA avant le calcul de la transformée de Fourier. Cette méthode de brouillage de phase est connue de l'homme de l'art. Une étape 74, mise en oeuvre par des moyens de calcul et de traitement, consiste à combiner l'image de phase (pH issue de la transformée de Fourier 70, l'image en noir et blanc 64 (IMAGE GBE,w) et l'image d'amplitude 68 (IMAGE GA) pour obtenir une image finale 76 (IMAGE F) dont le rendu visuel est identique au rendu visuel de l'image 60 (IMAGE G) et qui intègre un B11352 - DD13119JBD 14 hologramme obtenu à partir de l'image source de l'hologramme 62 (IMAGE H). Pour cela, chaque pixel de l'image finale 76 comprend une zone centrale entourée d'un fond, défini de la façon 5 suivante : - si le pixel de mêmes coordonnées de l'image 64 (IMAGE GBE,w) est à une première teinte, par exemple noire, la zone centrale est transparente et le fond opaque et si le pixel de mêmes coordonnées de l'image 64 (IMAGE GBE,w) est à une 10 deuxième teinte, par exemple blanche, la zone centrale est opaque et le fond transparent ; - la taille de la zone centrale dans chaque pixel est calculée à partir du niveau de gris du pixel de mêmes coordonnées de l'image 68 (IMAGE GA) ; 15 - si le pixel de mêmes coordonnées de l'image de phase cppi de la transformée de Fourier de l'image 62 est foncé (valeurs de phase tendant vers -n), la zone centrale dans le pixel de mêmes coordonnées de l'image finale 76 est fortement décalée d'un côté du pixel, si le pixel de mêmes coordonnées de l'image 20 de phase de la transformée de Fourier 72 cppi de l'image 62 est clair (valeurs de phase tendant vers +n), la zone centrale est décalée dans le pixel de mêmes coordonnées de l'image 76 vers un côté opposé du pixel, et si le pixel de mêmes coordonnées de l'image de phase 72 présente une teinte moyenne, la zone 25 centrale du pixel de mêmes coordonnées de l'image finale 76 est centrée dans le pixel (le décalage de la zone centrale dans le pixel est proportionnel à la valeur du pixel équivalent de l'image ((pH)). Une fois l'image finale, IMAGE F, obtenue, on réalise 30 une lithographie d'une couche opaque définie en surface d'une plaquette selon le motif de cette image. La figure 7 illustre certaines étapes de l'organigramme de la figure 6 par un exemple d'images obtenues par ces étapes. Dans cette figure, les images obtenues aux 35 différentes étapes sont référencées par la référence de l'étape B11352 - DD13119JBD 15 correspondante, suivie d'un '. Dans cette figure, la trame des images obtenues est visible pour une meilleure compréhension. La figure 8 illustre le principe de formation des éléments de la matrice (IMAGE G') à partir de l'image 60 (IMAGE G). La matrice 66 est constituée d'éléments dont la valeur est définie en fonction de l'écart, en valeur absolue, entre le niveau de gris du pixel correspondant de l'image 60 et le niveau de gris du pixel correspondant de l'image 64 (IMAGE GB00. Plus particulièrement, cette figure illustre une courbe de la valeur G' d'un élément de la matrice 66 (IMAGE G'), en fonction du niveau de gris G du pixel équivalent de l'image 60 (IMAGE G). Dans cette courbe, on voit que pour un pixel de l'image IMAGE G présentant un niveau de gris G compris entre 0 15 et 0,5 (pixel noir à pixel moyennement gris), la valeur de G' d'un élément de la matrice 66 augmente jusqu'à une valeur 13. Si le niveau de gris G d'un pixel de l'image en niveaux de gris, IMAGE G, est compris entre 0,5 et 1 (pixel moyennement gris à pixel blanc), la valeur de l'élément 20 équivalent de la matrice IMAGE G' décroit entre le seuil 13 et une valeur nulle lorsque le pixel est très clair. Pour permettre un codage optimal des hologrammes, on pourra prévoir un seuil minimal de la valeur des éléments de la matrice IMAGE G', pour des pixels de l'image IMAGE G 60 25 complètement noirs ou complètement blancs (respectivement G = 0 ou G = 1). Cela permet d'optimiser le rendu de la lecture finale de l'hologramme, comme nous le verrons ci-après. Le seuil minimal sera typiquement de 2 % des niveaux de gris maximum ou minimum (la taille minimum de la zone centrale des pixels 30 correspondant à un niveau de gris de 0,02 pour les tons sombres des pixels et à 0,98 pour les tons clairs des pixels, sur une échelle normalisée de niveaux de gris). La matrice 66 (IMAGE G') est utilisée pour obtenir une image d'amplitude IMAGE _GA (68) de l'onde associée à B11352 - DD13119JBD 16 l'hologramme, qui est utilisée pour déterminer la taille des zones centrales des pixels de l'image finale (IMAGE F). On notera que, avec une forme de la zone centrale elliptique, on peut optimiser la forme de l'ellipse pour que celle-ci présente une taille maximale, dans la direction y, égale au côté du pixel avec un taux de remplissage du pixel égal à 50 %. Par des calculs géométriques simples, on obtient qu'une ellipse présentant un rapport entre grand côté et petit côté égal à n/2 permet d'obtenir un taux de remplissage du pixel égal à 50 % pour une dimension de l'ellipse, selon la direction y, égale à la largeur du pixel. Ce rapport permet en outre un rendu visuel très doux. Pour obtenir l'image d'amplitude 68, il faut considérer les équations de diffraction d'une onde plane par une ouverture. La figure 9 illustre, en trois dimensions, un pixel comprenant une zone centrale transparente entourée d'un fond opaque. La zone centrale présente des dimensions appelées wy dans la direction y et wx dans la direction x. Lorsqu'une onde plane est incidente sur l'ouverture de la figure 9, ou de façon plus générale un motif de diffraction de coordonnées (n, m) d'un hologramme synthétique par codage d'ouvertures tel que proposé ci-dessus, un motif de dimensions wxnm et wynm génère une onde diffractée sphérique d'amplitude elnm reliée aux paramètres de l'ouverture par la relation non-linéaire suivante : wY An - E, nm J1 (n.wxnm), (1) 2 J1 étant la première fonction de Bessel. On notera que pour un codage optimal de la phase, on considère une forme 30 elliptique de l'ouverture telle que wy > wx. Dans une image en niveaux de gris, la surface de la portion centrale foncée d'un pixel normalisé vis-à-vis de la surface du pixel fixe la tonalité perçue par l'observateur. Cette tonalité I est reliée à la forme du pixel par la 35 relation : B11352 - DD13119JBD H-wx.-WYn. I - 4 Pour obtenir l'image d'amplitude de l'image initiale 60 (IMAGE G), on applique directement l'équation (1) à des ouvertures dont la surface est définie de façon que le rapport 5 Ip entre la surface des ouvertures et la surface du pixel corresponde à la valeur de l'élément correspondant de la matrice 66 (IMAGE G'). Pour assurer une bonne définition de l'image d'amplitude IMAGE GA 68, le seuil 13 peut être fixé à 0,5 (la surface de l'ouverture représentant la moitié de la surface du 10 pixel). On obtient ainsi une relation permettant de calculer les dimensions des zones centrales des pixels de l'image IMAGE F. En pratique, on définit un rapport de forme de l'ouverture p = wy/wx et on limite wy à une valeur maximale égale à 1 (côté des pixels normalisé). Les deux équations ci- 15 dessus permettent d'obtenir une relation entre la valeur de l'élément correspondant Inm de la matrice IMAGE G' et l'amplitude de l'onde difractée par ce pixel (A4 ), en fonction de p, qui est la relation suivante : = o qn j (3) p 20 L'équation ci-dessus permet ainsi de lier les éléments de la matrice IMAGE G', et donc le niveau de gris du pixel correspondant de l'image initiale IMAGE G, à l'amplitude du pixel de mêmes coordonnées de l'image IMAGE GA, c'est-à-dire à la taille de la zone centrale des pixels de l'image finale 25 IMAGE F. La figure 10 donne une représentation de cette fonction de transfert pour différentes valeurs de p. On notera que le choix d'un motif elliptique limite la gamme de valeurs d'intensité à wx = 1 et wy = 1, soit une intensité maximale de 30 0,8 dans le pixel. Au-delà, le motif déborde de la cellule et le codage de la phase se dégrade. Dans le cas du procédé ci-dessus, l'intensité I est limitée au seuil a prédéfini par la courbe de la figure 8, c'est-à-dire à 0,5. Pour optimiser l'amplitude nm 17 (2) B11352 - DD13119JBD 18 de l'onde émise tout en gardant une forme de motif quasi-circulaire, on choisit préférentiellement p = 1-1/2, comme nous l'avons vu ci-dessus. Dans ce cas, l'équation (3) se simplifie en : Ahrim = 1,1 2 .J1(V8Inm) ( 4 ) On notera que, pour simplifier la programmation de cette fonction de transfert, la fonction de Bessel pourra être remplacée par son développement série avec une bonne précision, ce qui permet d'obtenir l'équation suivante : Anm =(Inm)- (Inm)2 + 3 (Imj3 (5) Ainsi, les pixels de l'image finale obtenus par le procédé de la figure 6 présentent une structure tramée dont la zone centrale a une amplitude de trame holographique définie par l'équation ci-dessus à partir de l'élément de mêmes coordonnées de la matrice IMAGE G'. Le procédé de la figure 6 permet avantageusement d'obtenir une image présentant un visuel identique à l'image en niveaux de gris IMAGE G, dissimulant un hologramme synthétique par codage d'ouvertures et dont la lecture est directe. Ce procédé n'est toutefois pas optimal du fait que l'image d'amplitude AH de la transformée de Fourier de l'image source de l'hologramme IMAGE H a été supprimée et remplacée par l'image d'amplitude IMAGE GA. La reconstruction de l'image IMAGE H est donc dégradée, comme on le verra à la suite. On notera que l'image d'amplitude IMAGE HA dont la transformée de Fourier est calculée à l'étape 70 de la figure 6 est obtenue à partir de l'image source IMAGE H 62 en suivant le même procédé que celui d'obtention de l'image 68 (IMAGE GA) à partir de la matrice 66 (IMAGE G'), c'est-à-dire en calculant 30 l'amplitude diffractée par des ouvertures dont la taille dépend du niveau de gris du pixel équivalent de l'image source IMAGE H. La figure 11 illustre un procédé d'optimisation permettant d'obtenir une image de phase 911 issue de la transformée de Fourier de l'image d'amplitude IMAGE HA et de B11352 - DD13119JBD 19 l'image d'amplitude IMAGE GA, permettant de pallier à la dégradation à la lecture du procédé décrit ci-dessus. Le procédé d'optimisation de l'image de phase imposée à l'image finale proposé ici est une méthode d'optimisation itérative décrite dans la publication de Gerchberg et Saxton intitulée "A pratical algorithm for the determination of the phase from image and diffraction plane pictures", Optik, Vol. 35, 237-246. Partant de l'image source de l'hologramme, IMAGE H (62), on détermine l'image d'amplitude de cette image, IMAGE HA (80), de façon similaire à la détermination de l'amplitude de la trame holographique IMAGE GA issue de la matrice IMAGE G' (66). A une étape suivante 82, l'image d'amplitude IMAGE HA 80 est combinée à une composante de phase d'initialisation (pi 83 15 issue soit d'une distribution de phase aléatoire du type brouillage de phase, soit d'une composante de phase continue, pour obtenir une valeur complexe (initialisation du processus). On calcule ensuite la transformée de Fourier de l'image complexe obtenue à l'étape 82, à une étape 84. 20 L'image d'amplitude A obtenue par la transformée de Fourier de l'étape 84 est laissée de côté (tout comme l'image d'amplitude AH en figure 6). A une étape 86, une valeur complexe est calculée par combinaison de l'image d'amplitude, IMAGE GA, 68, obtenue à partir de l'image en niveaux de gris 60 et l'image 25 de phase cp obtenue à l'étape 84. On calcule ensuite, à une étape 88, la transformée de Fourier inverse de la valeur complexe obtenue à l'étape 86. L'image d'amplitude A du résultat obtenu est laissée de côté. L'image de phase cp obtenue par cette transformée de Fourier 30 inverse, quant à elle, est introduite en entrée de l'étape 82 en lieu et place de la distribution de phase d'initialisation du processus (pi, en combinaison avec l'image d'amplitude IMAGE HA. L'étape 84 de calcul de la transformée de Fourier de la valeur complexe obtenue à l'étape 82 est réalisée de nouveau pour 35 obtenir une image de phase (pH optimisée (étape 90), du fait de B11352 - DD13119JBD 20 la prise en compte de l'image 68, IMAGE GA. Le cycle d'étapes 82 à 88 peut être réitéré plusieurs fois pour optimiser l'image de phase cppi obtenue à l'étape 90. L'image de phase utilisée directement dans remplacement de l'image de réalisation du cycle décrit intègre des informations l'hologramme, IMAGE H, mais obtenue à l'étape 90 peut être le procédé de la figure 6 en phase 72. En effet, du fait de la ci-dessus, cette image de phase (PH issues de l'image source de intègre également un correctif qui permet de contrecarrer le fait que la taille des zones centrales des pixels de l'image finale, IMAGE F, obtenue ne dépend pas de l'image source de l'hologramme mais des niveaux de gris de l'image visible en niveaux de gris, IMAGE G. La réalisation du cycle d'amélioration de l'image de phase cppi ci-dessus permet l'amélioration de la lecture de l'hologramme intégré dans l'image tramée. On notera que le cycle ci-dessus pourra être réitéré autant de fois que nécessaire pour que l'image de phase cppi finale de l'étape 90 comprenne un correctif efficace de la perte d'information liée à la non- utilisation de l'image d'amplitude de la transformée de Fourier de l'image source de l'hologramme. Au bout d'un certain nombre d'itérations, dont le verra ci-après, on obtient une est utilisée dans le procédé de nombre peut varier comme on le image de phase cppi optimisée qui la figure 6. La figure 12A illustre un exemple d'une image source de l'hologramme IMAGE H. Les figures 12B, 12C et 12D illustrent des résultats de simulation obtenus sur la base d'un dispositif de lecture tel que celui de la figure 3 sur une image tramée 30 obtenue par le procédé de la figure 6. Plus précisément, l'image 12B illustre la lecture sur une image obtenue par le procédé de la figure 6 sans réaliser de cycle d'ajustement de l'image de phase de la figure 11, l'image 12C illustre la lecture sur une image obtenue par le procédé de 35 la figure 6 en réalisant trois fois le cycle d'ajustement de la B11352 - DD13119JBD 21 figure 11 et la figure 12D illustre la lecture sur une image obtenue par le procédé de la figure 6 en réalisant 300 fois le cycle d'ajustement de la figure 11. De ces figures, on peut noter que la réalisation du procédé d'optimisation de la figure 11 permet d'obtenir, même avec un faible nombre de cycles, un correctif efficace de l'image de phase appliquée à l'image finale, prenant en compte la perte d'information initiale de l'amplitude de la transformée de Fourier de l'image source de l'hologramme. En effet, la lecture réalisée sur une image obtenue par le procédé de la figure 6 après trois réitérations seulement du cycle de la figure 11 (figure 12C) est de très bonne qualité, quasiment de même qualité que la lecture de la figure 12D. Dans le cadre du procédé ci-dessus, la qualité de la 15 visualisation de l'image en niveau de gris est un avantage recherché. La figure 13 illustre le résultat obtenu à partir d'une image représentée en insert, à une résolution d'écran qui laisse apparaître la trame (les zones centrales des pixels 20 présentant un rapport de forme p = 1,5). La qualité du rendu de gris est largement satisfaisante lorsque l'image est observée à une distance suffisante, l'éloignement de la figure introduisant un flou visuel qui gomme l'effet de la trame. Pour simuler cet effet de flou, on a cherché à rendre 25 artificiellement floue l'image avec un flou gaussien de largeur de 2 et de 4 pixels. Les figures 14A et 14B montrent les résultats obtenus, la figure 14D illustrant l'image initiale dont sont issues les images des figures 14A et 14B : on voit apparaître une granularité sur l'image liée à certaines 30 conjonctions de déplacement des motifs de trame. Comme on l'a présenté précédemment, la position des motifs de trame est définie par la valeur de la phase affectée au pixel. Si deux pixels successifs ont des valeurs de phase dont l'écart est important, la distance entre les motifs sera 35 également importante. Suivant que le fond de trame est clair ou B11352 - DD13119JBD 22 sombre cet écart se traduira par l'apparition d'un point blanc ou noir. La figure 15 montre un détail de l'image IMAGE G en niveaux de gris à partir duquel on calcule les images IMAGE GA 5 et cppl. Ces grandeurs permettent de calculer la répartition des motifs représentés sur l'image IMAGE F. On voit qu'aux endroits où la phase passe d'une valeur maximale (pixel blanc de (pH) à une valeur minimale (pixel noir de (pH) on obtient un décalage opposé des deux motifs qui conduit à l'apparition d'un "vide". 10 Après l'opération de flou gaussien, ce vide se traduit par l'apparition de pixels sombres sur l'image résultante (IMAGE F'). Ce phénomène est gênant si la détection de l'image reconstruite doit être rigoureusement identique à l'original. 15 C'est généralement le cas en préservation des données. Il est toutefois possible de réduire cet effet. Pour cela, on peut citer les méthodes suivantes : - après simulation de l'effet d'écriture et de détection, on modifie l'image source pour anticiper l'opaci- 20 fication et l'éclaircissement de certains pixels ; - on introduit directement dans la boucle de calcul de la phase une modification de l'amplitude pour prendre en compte la présence des sauts de phase ; - on sur-échantillonne l'image source IMAGE G pour 25 augmenter le flou sans perte de résolution image. Le résultat de cette dernière proposition est illustré en figure 14C, sur une image sur-échantillonnée par un facteur 2, en conservant un flou de 2 pixels image. On voit nettement l'amélioration sur le rendu de niveau de gris. 30 Les figures 16B et 16C illustrent le résultat d'un flou sur une image tramée de la façon proposée ci-dessus, et par les techniques antérieures, l'image initiale étant illustrée en figure 16A. Dans le cas de l'art antérieur, l'absence de cellules sombres et claires comme dans notre cas, provoque des 35 défauts d'une seule tonalité, claire dans notre exemple. Ces B11352 - DD13119JBD 23 défauts sont donc particulièrement visibles quand ils apparaissent sur le fond de tonalité opposée (figures 16B et 16C). Le procédé proposé ici résout ce problème. De nombreuses applications du procédé ci-dessus sont 5 envisageables. Une première application concerne le stockage analogique des images. La solution proposée permet une préservation pérenne des données graphiques grâce à un mode de sauvegarde visuel. La récupération des données se fait sans 10 décodage numérique et ne risque pas l'obsolescence du format de lecture. Cette solution a l'inconvénient d'un mode de récupération des données relativement long ; une solution serait donc de conserver la pérennité du stockage analogique tout en 15 utilisant la rapidité de récupération de données numériques. Prenons l'exemple d'une image de 1000x1000 pixels, codée sur 256 niveaux de gris (1 octet). Son poids numérique, au format bitmap sans compression, est d'environ 103 ko. L'idée est de coder une version numérique de l'image 20 source de l'hologramme IMAGE H. Pour faciliter la détection, une représentation binaire de l'image source de l'hologramme IMAGE H et choisie, en supposant que seulement 50 % de la surface de l'hologramme est exploitable. Cette restriction est liée à la difficulté qu'il y a à obtenir des images non bruitées en 25 optique cohérente. Cela conduit à une capacité de stockage en bits correspondant à 50 % du nombre de pixels. On obtient une capacité numérique d'environ 62 ko. Cette réduction de la taille du fichier peut être obtenue en stockant une image compressée de l'image d'origine IMAGE G, par exemple au format jpeg. 30 L'utilisateur peut donc avoir le choix entre une détection numérique rapide des images numériques au format compressé ou une détection analogique lente des fichiers sources au format natif. Il est à noter que la détection de l'image sur le 35 format compressé peut aussi être utilisée pour améliorer la B11352 - DD13119JBD 24 qualité de l'image source reconstruite, en aidant notamment à corriger l'effet de granularité. Les figures 17A à 17D illustrent deux exemples d'images IMAGE G suivant ce principe, en vue complète en ce qui concerne les images 17A et 17C, et en vue agrandie en ce qui concerne les figures 17B et 17D. Plus particulièrement, les figures 17A et 17B illustrent une image au format bimap présentant une taille de 1000 ko et une image au format jpeg présentant une taille de 62 ko. Dans le même domaine d'application, on peut également proposer une variante pour augmenter la capacité de surface du média de stockage analogique. En plus des données analogiques, le média de stockage doit pouvoir préserver les métadonnées associées aux documents visuels. Ces données sont contenues en arrière plan du fichier numérique et peuvent concerner toutes formes d'information (nom du fichier, date de création, format du fichier, commentaires, etc...). Ces informations peuvent être visualisées au format binaire (suite de caractères) et peuvent accepter une certaine dégradation (la forme passe derrière le fond). Ce type de données s'applique donc bien au stockage par holographie synthétique. La figure 18 présente cette application. Un faisceau laser traverse la représentation tramée de l'image, formée sur un support 100, dont le rendu visuel est en niveaux de gris. Derrière une lentille de Fourier 102, les métadonnées 104 associées sont visibles. Ces métadonnées sont accessible visuellement à l'aide d'un laser et d'une optique de conception simple. Avantageusement, l'accès à leur contenu ne dépend pas d'un format de décodage numérique. La renumérisation des données analogiques étant un processus long, la représentation tramée de l'image IMAGE G doit être imagée sur un capteur matriciel pour reconstituer la matrice des pixels de niveaux de gris d'origine. Ce processus se fait par une succession de prises d'images qui doivent ensuite B11352 - DD13119JBD 25 être raboutées. Si l'image présente peu de détails, l'opération de raboutage peut être compliquée. Une solution peut être de disposer au sein de l'image des hologrammes d'alignement comme présenté en figure 19. La représentation tramée 104 présente une série d'hologrammes 106. Chaque hologramme 106 forme dans le plan focal d'une lentille de Fourier 108 un point lumineux 110 à une coordonnée précise. Si le faisceau sonde est au croisement exact de quatre hologrammes, on détecte dans le plan de lecture quatre points lumineux 110 différents. Le positionnement de la plaquette sur laquelle sont formés les hologrammes sur l'axe optique peut donc être asservi à la détection de ces points lumineux. Le système peut balayer l'image en se positionnant en des endroits précis, pour lesquels l'intensité des quatre points est équilibrée, indépendamment du contenu visuel de l'image. Une autre application du procédé proposé ici, et qui sort du cadre du stockage analogique, consiste en la sécurisation de documents d'identité, et de façon plus générale de tout type de documents. On peut envisager de coder dans la photographie des cartes d'identité un hologramme permettant d'attester de l'authentification du document. Des modes de réalisation particuliers de la présente invention ont été décrits. Diverses variantes et modifications apparaîtront à l'homme de l'art. En particulier, d'autres applications de l'intégration d'un hologramme par codage d'ouvertures dans une image tramée dont le rendu est un rendu en niveaux de gris pourront être prévues. En outre, on notera que des procédés de brouillage de l'hologramme connus, impliquant l'impossibilité d'une lecture directe de l'hologramme, pourront être prévus en combinaison avec le procédé décrit ici. Dans ce cas, l'utilisation d'une clé de phase ou d'un masque de lecture pourra être prévu lors de la lecture de l'hologramme, par exemple avec un dispositif de lecture tel que celui de la figure 3. B11352 - DD13119JBD 26 Divers modes de réalisation avec diverses variantes ont été décrits ci-dessus. On notera que l'homme de l'art pourra combiner divers éléments de ces divers modes de réalisation et variantes sans faire preuve d'activité inventive. REVENDICATIONS1. Procédé d'intégration d'un hologramme synthétique dans une image directement observable d'une scène, comprenant les étapes suivantes : définir une première image (IMAGE G) de ladite scène, 5 comprenant des pixels en niveaux de gris ; fabriquer des première et deuxième matrices à partir de la première image, la première matrice (IMAGE GBE,w) comprenant des pixels d'une première et d'une deuxième teinte selon que le pixel correspondant de la première image (IMAGE G) 10 a un niveau de gris inférieur ou supérieur à un seuil (a), chaque élément la deuxième matrice (IMAGE G') comprenant une valeur égale à l'écart de niveau de gris entre le pixel correspondant de la première image (IMAGE G) et le pixel correspondant de la première matrice (IMAGE GBE,w) ; 15 fabriquer une troisième matrice de pixels ((pH) en calculant l'image de phase de la transformée de Fourier d'une image issue d'une deuxième image source de l'hologramme (IMAGE H) ; fabriquer une quatrième matrice de pixels (IMAGE F), 20 chaque pixel de la quatrième matrice comprenant une zone centrale dont la surface est déterminée à partir de l'élément correspondant de la deuxième matrice (IMAGE G') et dont un décalage par rapport au centre du pixel est issu du pixel correspondant de la troisième matrice ((pH) ; et 25 réaliser une lithographie d'une couche opaque en surface d'une plaquette selon le motif défini par la quatrième matrice de pixels (IMAGE F), la lithographie étant prévue pour éliminer, dans un pixel, la couche opaque à l'intérieur de la zone centrale si le pixel équivalent de la première matrice 30 (IMAGE GBE,w) est de la première teinte et pour éliminer la couche opaque à l'extérieur de la zone centrale si le pixel équivalent de la première matrice (IMAGE GBE,w) est de la deuxième teinte.B11352 - DD13119JBD 28 2. Procédé selon la revendication 1, dans lequel la taille de la zone centrale de chacun des pixels de la quatrième matrice (IMAGE F) est calculée à partir de la valeur du pixel correspondant d'une cinquième matrice de pixels (IMAGE GA), chaque pixel de la cinquième matrice (IMAGE GA) étant obtenu par calcul de l'amplitude diffractée d'une ouverture définie dans un pixel, le rapport entre la surface de l'ouverture et la surface du pixel étant égal à la valeur de l'élément correspondant de la deuxième matrice (IMAGE G'). 3. Procédé selon la revendication 2, dans lequel l'amplitude el d'un pixel de la cinquième matrice (IMAGE GA) est liée à la valeur I de l'élément correspondant de la deuxième matrice (IMAGE G') par la relation suivante : A = 1,1 . (V8I. ) J1 étant la première fonction de Bessel. 4. Procédé selon la revendication 2 ou 3, dans lequel la troisième matrice ((pH) est obtenue en réalisant les étapes suivantes : (a) calculer une sixième matrice d'amplitude (IMAGE HA) obtenue à partir de la deuxième image source de l'hologramme (IMAGE H) en réalisant les mêmes étapes que celles permettant d'obtenir la cinquième matrice (IMAGE GA) à partir de la deuxième matrice (IMAGE G') ; (b) combiner les pixels de la sixième matrice 25 d'amplitude (IMAGE HA) à une distribution de phase aléatoire ou à une composante continue pour obtenir une valeur complexe ; et (c) calculer l'image de phase de la transformée de Fourier de ladite valeur complexe. 5. Procédé selon la revendication 2 ou 3, dans lequel 30 la troisième matrice ((pH) est obtenue par un procédé d'optimisation tirant profit de la cinquième matrice (IMAGE GA). 6. Procédé selon la revendication 5, dans lequel le procédé d'optimisation consiste à : (a) calculer une sixième matrice d'amplitude 35 (IMAGE HA) obtenue à partir de la deuxième image source deB11352 - DD13119JBD 29 l'hologramme (IMAGE H) en réalisant les mêmes étapes que celles permettant d'obtenir la cinquième matrice (IMAGE GA) à partir de la deuxième matrice (IMAGE G') ; (b) combiner (82) les pixels de la sixième matrice 5 d'amplitude (IMAGE HA) à une distribution de phase aléatoire ou à une composante de phase continue pour obtenir une première valeur complexe ; (c) calculer (84) la transformée de Fourier de la première valeur complexe ; 10 (d) combiner (86) les pixels de la cinquième matrice de pixels (IMAGE GA) et les pixels de l'image de phase issue de la transformée de Fourier de l'étape (c) pour obtenir une deuxième valeur complexe ; (e) calculer (88) la transformée de Fourier inverse de 15 la deuxième valeur complexe ; (f) combiner (82) les pixels de la sixième matrice d'amplitude (IMAGE HA) aux pixels de l'image de phase issue du calcul (88) de l'étape (e) pour obtenir une nouvelle première valeur complexe ; 20 (g) réitérer l'étape (c) pour calculer (84) la transformée de Fourier de ladite nouvelle première valeur complexe, la troisième matrice ((pH) correspondant à l'image de phase obtenue de cette dernière transformée de Fourier. 7. Procédé selon la revendication 6, dans lequel le 25 procédé d'optimisation comprend en outre la réitération des étapes (d) à (g), la troisième matrice ((pH) correspondant à l'image de phase obtenue à l'étape (g) après plusieurs cycles de réitération des étapes (d) à (g). 8. Procédé selon la revendication 7, dans lequel au 30 moins trois cycles sont réalisés. 9. Procédé selon l'une quelconque des revendications 1 à 8, dans lequel le décalage de la zone centrale par rapport au centre du pixel de chaque pixel de la quatrième matrice (IMAGE F) est directement proportionnel de la valeur du pixel 35 correspondant de la troisième matrice ((pH).B11352 - DD13119JBD 30 10. Procédé selon l'une quelconque des revendications 1 à 9, dans lequel le seuil (a) de niveaux de gris est égal à 0,5, à 10 % près, sur une échelle normalisée de niveaux de gris. 11. Procédé selon l'une quelconque des revendications 1 5 à 10, dans lequel la zone centrale de chaque pixel de la quatrième matrice (IMAGE F) est elliptique. 12. Procédé selon la revendication 11, dans lequel le rapport de forme de la zone centrale elliptique de chaque pixel de la quatrième matrice (IMAGE F) est égal à 1-1/2. 10 13. Procédé selon l'une quelconque des revendications 1 à 12, dans lequel la première teinte de la première matrice (IMAGE GBE,w) est une teinte minimum, la deuxième teinte de la première matrice (IMAGE GBE,w) est une teinte maximum, la lithographie étant prévue pour éliminer la couche opaque à 15 l'intérieur de la zone centrale d'un pixel si le pixel équivalent de la première matrice (IMAGE GBE,w) est de la première teinte, et pour éliminer la couche opaque à l'extérieur de la zone centrale d'un pixel si le pixel équivalent de la première matrice (IMAGE GBE,w) est de la deuxième teinte. 20 14. Procédé selon l'une quelconque des revendications 1 à 13, dans lequel une taille minimum de la zone centrale des pixels de la quatrième matrice (IMAGE F) est imposée. 15. Procédé selon la revendication 14, dans lequel la taille minimum de la zone centrale des pixels de la quatrième 25 matrice (IMAGE F) correspond à un niveau de gris de 0,02 pour les tons sombres desdits pixels et à 0,98 pour les tons clairs desdits pixels, sur une échelle normalisée de niveaux de gris.