FR 3054354 A1 20180126 FR 1657064 A 20160722 1 La présente invention concerne un procédé de production d'image et plus particulièrement une amélioration d'un tel procédé dans le cas où ce procédé de production d'image produit une image présentant une épaisseur telle que se produise un effet d'ombrage. Certains procédés de production d'image produisent une image présentant une épaisseur. Quel que soit son mode de production, lorsqu'une image est observée, ceci s'effectue, le plus souvent, sous illumination multidirectionnelle. Dans le cas d'une image définie en épaisseur ceci a pour effet préjudiciable que chaque élément d'image projette une ombre, fonction de l'opacité de l'élément, sur les éléments d'image voisins. L'étendue du voisinage ainsi affecté est fonction de l'épaisseur des éléments d'image. Un tel effet d'ombrage modifie l'image effectivement 20 réalisée relativement à une image objective souhaitée, en ce que l'image effectivement réalisée est plus sombre et plus floue que l'image objective. Il peut être noté que ce phénomène d'ombrage préjudiciable n'existe pas pour une image imprimée de manière 25 classique. En effet la profondeur/épaisseur de l'encre d'impression est suffisamment négligeable par rapport aux dimensions transverses d'un élément d'image pour qu'un ombrage puisse se produire, ou du moins soit perceptible. La présente invention remédie à cet inconvénient. Le 30 principe de l'invention consiste à améliorer la consigne de production fournie au procédé de production d'image de telle manière à ce que l'image produite sur la base de cette consigne de production améliorée soit la plus proche possible de l'image objective que l'on souhaite et que les 35 conséquences de l'effet d'ombrage soient les plus réduites possibles ou encore pour corriger les effets de la dynamique, des distorsions liées à l'utilisation du laser ou des effets de la diffusion du matériau constituant le substrat recevant 3054354 2 l'impression. Pour cela l'invention modélise le dispositif de production d'image et son effet d'ombrage par une fonction de transfert fournissant une image prédite lorsque le dispositif 5 de production d'image est appliqué à une consigne de production. Cette fonction est ensuite « inversée » pour déterminer une consigne de production améliorée optimale telle que le dispositif de production d'image appliqué à cette consigne de production améliorée produise une image 10 réalisée la plus proche possible de l'image objective. L'invention a pour objet un procédé de production d'une image objective au moyen d'un dispositif de production d'image apte à produire une juxtaposition d'éléments d'image dans un substrat présentant une épaisseur et produisant une 15 image en épaisseur telle que se produise un effet d'ombrage de chaque élément d'image sur ses éléments d'image voisins, comprenant les étapes suivantes : détermination d'une consigne de production à partir de l'image objective, production d'une image réalisée au moyen du dispositif de production d'image utilisant ladite consigne de production, où l'étape de détermination remplace la consigne de production par une consigne de production améliorée déterminée par les étapes suivantes : modélisation du procédé initial de production d'image par une fonction de transfert transformant une consigne de production en une image prédite, détermination d'une consigne de production améliorée dont l'image prédite par la fonction, est la plus proche possible de l'image objective de telle manière à réduire les conséquences de l'effet d'ombrage dans l'image réalisée. Selon une autre caractéristique, l'effet d'ombrage se produit lorsque l'épaisseur du substrat de l'image est au moins égale au double de la dimension transverse d'un élément d'image. Selon une autre caractéristique le dispositif de 35 production d'image est un dispositif de gravure laser apte à noircir un substrat transparent dans son épaisseur et une consigne de production est une carte de tir laser. Selon une autre caractéristique le dispositif de 3054354 3 production d'image est un dispositif de micro perforation apte à réaliser des micro trous dans l'épaisseur d'un substrat transparent et une consigne de production est une carte de perforation. 5 Selon une autre caractéristique la modélisation du dispositif de production d'image détermine une transmittance de la consigne de production comme résultat intermédiaire. Selon une autre caractéristique la fonction de transfert combine une première fonction de transmittance transformant 10 une consigne de production en une transmittance, une deuxième fonction de réflectance transformant une transmittance en une réflectance et une troisième fonction d'image transformant une réflectance en une image prédite. Selon une autre caractéristique la première fonction est 15 définie par une forme paramétrique, et les paramètres sont déterminés par calibrage. Selon une autre caractéristique la première fonction est déterminée empiriquement, et tabulée. Selon une autre caractéristique la deuxième fonction est 20 donnée par la formule suivante : R = F 2(T) - b.lc1*T2 + (1 - b).1(1.(T * M).T 1- krT2 1- k2.(T2 *M) avec R une matrice de réflectance, beq0,1] un paramètre scalaire, k1,k2E[0,1] deux constantes optiques, T une matrice de transmittance, * l'opérateur de convolution et M un noyau de 25 convolution. Selon une autre caractéristique la deuxième fonction inclut de plus un bruit K, selon la formule : R = F 2(T) -(b.lerT22 + (1- b),Ici.(T * M).T *K avec 1- krT 1- k2.(T2 *M) * l'opérateur de convolution et K un noyau de convolution. 30 Selon une autre caractéristique K est un noyau bidimensionnel gaussien de déviation standard a. Selon une autre caractéristique 4-41-/J.U1-ri) et k2=rgfi, avec rs indice de réflexion spéculaire de l'interface air image, ri indice de réflectance interne de Fresnel de 35 l'interface air image et rg indice de réflectance intrinsèque 3054354 4 d'un substrat réfléchissant lambertien. Selon une autre caractéristique le noyau de convolution M est une matrice carrée unité de dimensions dxd, normalisée, M = [1/d2]cixd- 5 Selon une autre caractéristique la troisième fonction est la fonction identité. Selon une autre caractéristique la fonction de transfert est calibrée en identifiant les paramètres b, d et éventuellement la déviation standard o, par résolution de 10 l'équation suivante : b,d,o- = a rg min A(IRp,IPp ) = a rg min 3' 3.(IRp,F(Sp , avec ô1] bEloa] d>0 d>0 ô-k0 erk0 S une consigne de production, IR une image obtenue au moyen du dispositif de production utilisant la consigne de production S, IP une image prédite par la fonction de 15 transfert F à partir de la consigne de production S, soit IP = F(S), p un indice parcourant les pixels, F la fonction de transfert et 8 une distance colorimétrique entre deux pixels. Selon une autre caractéristique la résolution est alimentée par des paires de calibrage comprenant une consigne 20 de production et une image réalisée au moyen du dispositif de production utilisant la consigne de production. Selon une autre caractéristique l'étape de détermination d'une consigne de production améliorée SA comprend un calcul itératif automatique minimisant une distance entre l'image 25 objectif IO et une image prédite IP par la fonction de transfert F appliquée à la consigne de production améliorée SA. Selon une autre caractéristique l'étape de détermination d'une consigne de production améliorée SA comprend la 30 résolution de l'équation suivante : SA = a rg min (3(/0p,IPp))2 = a rg min (A(I0P F(SP )))2 , avec S 7 g 7ei S une consigne de production, IO une image objective, IP une image prédite par la fonction de transfert à partir de la consigne de production S, soit IP = F(S), p un indice parcourant les pixels, F la fonction de transfert et .8 une distance colorimétrique entre deux pixels. 3054354 5 L'invention concerne encore une image produite par un tel procédé. L'invention concerne encore un dispositif de production d'image apte à produire une juxtaposition d'éléments d'image 5 dans un substrat présentant une épaisseur et produisant une image en épaisseur telle que se produise un effet d'ombrage de chaque élément d'image sur ses éléments d'image voisins, comprenant : un moyen d'acquisition d'une image objectif IO, un moyen de détermination d'une consigne de production S à 10 partir de l'image objective IO, un dispositif de production d'une image réalisée IR à partir de ladite consigne de production S, où le moyen de détermination comprend encore des moyens de mise en oeuvre d'un tel procédé de sorte à déterminer une consigne de production améliorée SA pour 15 remplacer la consigne de production S. D'autres caractéristiques, détails et avantages de l'invention ressortiront plus clairement de la description détaillée donnée ci-après à titre indicatif en relation avec 20 des dessins sur lesquels : - la figure 1 illustre en vue de dessus un support comprenant une image réalisée par gravure laser, - la figure 2 illustre en vue de profil le support de la figure 1, 25 - la figure 3 illustre l'action d'un éclairant multidirectionnel sur un élément d'une image gravée laser et le problème d'ombrage, - la figure 4 présente un exemple de fonction t(s). 30 Les figures 1 et 2 illustrent un mode de réalisation d'un dispositif de production d'image. Tel qu'illustré à la figure 1, ce dispositif crée une image composée d'une juxtaposition d'éléments d'image 2 ou pixels, dans un support ou substrat 1. La particularité des procédés de production d'image 35 envisagés dans la présente est que les éléments d'image 2 présentent une épaisseur importante, comme cela est plus particulièrement visible en vue de profil à la figure 2. Les éléments d'images 2 comprennent une teinte, une opacité et/ou 3054354 6 une hauteur qui peuvent être variables. Il en résulte, tel qu'illustré à la figure 3, lorsque le support 1 est soumis à un éclairage 3, typiquement multidirectionnel, ce qui est généralement le cas lorsque 5 l'image est observée, que chaque élément d'image 2 projette sur ses voisins une ombre 4. Ceci a pour effet préjudiciable une augmentation du flou et un assombrissement de l'image produite. Il est considéré qu'un tel effet d'ombrage se produit, au 10 moins de manière à être détecté, lorsque l'épaisseur de l'image, soit plus précisément l'épaisseur du substrat supportant l'image 1, est au moins égale à deux fois la dimension transverse d'un élément d'image 2. L'épaisseur du substrat 1 détermine la hauteur/épaisseur maximale d'un 15 élément d'image 2. La dimension transverse d'un élément d'image 2 est une de ses dimensions dans le plan image, soit selon les cas sa largeur, sa longueur, ou son diamètre. Le rapport épaisseur/dimension transverse est au moins égal à deux. Il ne présente pas de limite supérieure. Pour des 20 réalisations par gravure laser, des points laser de diamètre typique 40pm peuvent être réalisés dans une couche sensible au laser d'épaisseur pouvant atteindre 400pm, soit un rapport 10. Dans les modes de réalisation classiquement utilisés, ce rapport est d'environ 5, l'épaisseur du substrat 1 étant 25 comprise entre 150p et 250pm. Les dispositifs de production d'image créant une image en épaisseur telle que se produise un effet d'ombrage comprennent les dispositifs par gravure laser. Le principe de production d'une image par gravure laser 30 utilise un support 1 dont au moins une partie est sensible au laser en ce qu'elle noircit au point où un tir laser est réalisé et rencontre ladite partie. L'épaisseur d'un tel élément d'image 2 et donc son niveau d'opacité, sont fonction de l'intensité et/ou durée du tir laser. Un laser permet de 35 tirer avec précision sur un point donné et avec une intensité finement modulée. Une série de tir laser permet ainsi de créer, élément d'image 2 par élément d'image 2, une image composée d'éléments d'image 2 juxtaposés d'opacité variable 3054354 7 et choisie. L'utilisation directe de ce principe permet de produire une image monochrome. La résolution des niveaux de gris est importante en liaison avec la résolution d'intensité du 5 laser. Il est encore connu de réaliser une image polychrome en appliquant le principe ci-dessus à un support 1 comprenant, superposé avec la couche sensible au laser, une matrice de points de couleur. Un élément noirci au laser permet de 10 proportionnellement masquer un point de couleur. Avec une matrice comprenant des couleurs choisies dans une base colorimétrique, des points laser masquent proportionnellement et sélectivement les points de couleur de manière à exprimer une teinte par composition des couleurs de base non masquées. 15 Les dispositifs de production d'image créant une image en épaisseur telle que se produise un effet d'ombrage comprennent encore les dispositifs par micro perforation. Le principe de production d'une image par micro perforation utilise un support dans lequel sont réalisés des 20 micro trous. En faisant varier la taille, le nombre et/ou la profondeur des micro trous, il est possible de créer un élément d'image d'opacité variable déterminée. Une série de micro trous permet ainsi de créer, élément d'image 2 après élément d'image 2, une image composée d'éléments d'image 2 25 juxtaposés d'opacité variable et choisie. L'utilisation directe de ce principe permet de produire une image monochrome. Tous ces dispositifs ont en commun de produire une image en épaisseur. Cette épaisseur entraîne un effet préjudiciable 30 d'ombrage qui floute et assombrit l'image réalisée par rapport à l'image objective. L'invention est applicable à tous ces dispositifs. Afin de réduire les conséquences de l'effet d'ombrage l'invention propose de réutiliser le dispositif de production 35 d'image dans un procédé de production d'image, après amélioration de la consigne de production. Le dispositif de production d'image produit une image I, ou image réalisée IR, à partir d'une consigne de production S. Ce procédé de 3054354 8 production d'image peut être quelconque. Dans le cas d'un dispositif de micro perforation, la consigne de production S est une carte de perforation définissant où un élément d'image doit être réalisé et 5 combien de trous, de quelles tailles et de quelles profondeurs comporte cet élément d'image. Dans le cas d'un dispositif de gravure laser, la consigne de production S est une carte de tir. Une carte de tir S est une matrice bidimensionnelle comprenant autant d'éléments que 10 d'éléments d'image 2, chaque élément étant spatialement associé à un élément d'image let comprenant une intensité de tir du laser. Cette intensité est typiquement indiquée par une valeur comprise entre 0 et 255. 0 indiquant, par exemple, une absence de tir laser et donc un point « blanc » et 255 15 indiquant une intensité maximale et donc un point « noir ». Alternativement une carte de tir S peut être normalisée, chacun de ses éléments étant compris entre 0 et 1. Selon l'art antérieur une consigne de production S est déterminée par analyse de l'image objective IO que l'on veut 20 produire. Cette consigne de production S est fournie au dispositif de production d'image qui réalise la production (série de tirs laser ou de micro trous) selon la consigne de production S pour produire une image réalisée IR qui est alors perturbée par l'effet d'ombrage, et donc différente de 25 l'image objective IO. Selon l'invention, le procédé réutilise le dispositif de production d'image, mais au préalable remplace la consigne de production S de l'art antérieur par une consigne de production améliorée SA. Cette consigne de production 30 améliorée SA est fournie au dispositif de production d'image qui réalise la production d'une image de manière identique, mais en utilisant la consigne de production améliorée SA, pour produire une nouvelle image réalisée IR améliorée. Cette image est améliorée en ce qu'elle est moins sujette aux 35 conséquences de l'effet d'ombrage. Ceci est réalisé en prenant en compte l'effet d'ombrage de manière à en inverser les effets, dans la mesure du possible, sous forme de corrections intégrées dans la 3054354 9 nouvelle consigne de production SA. Pour cela le procédé selon l'invention procède en deux étapes permettant de calculer une consigne de production SA améliorée. Une première étape consiste à modéliser le dispositif de 5 production d'image sous forme d'une fonction de transfert F transformant une consigne de production S en une image prédite IP. Cette image issue de la fonction de transfert est nommée image prédite IP. Au cours d'une deuxième étape, le procédé réalise une 10 inversion, au moins locale, de cette fonction de transfert F afin de déterminer une consigne de production améliorée SA, telle que l'image IP prédite, transformée de cette consigne de production SA par la fonction de transfert F, IP = F(SA), soit la plus proche possible de l'image objective IO, soit 15 F(SA) = IP IO. La fonction de transfert F qui transforme une consigne de production S en une image prédite IP peut être obtenue selon toute méthode et/ou modélisation. Selon une caractéristique avantageuse, il peut être 20 utilisé pour représenter le dispositif de production d'image IP = F(S), une modélisation déterminant la transmittance T de la consigne de production S comme résultat/variable intermédiaire. Ceci revient à décomposer la fonction de transfert F en deux fonctions G1, G2 telles que IP = F(S) = 25 G2(G1(S)), où G1 est une fonction qui fournit la transmittance T en fonction de la consigne de production S, T = G1(S) et G2 est une fonction qui fournit l'image prédite IP en fonction de la transmittance T. Pour mémoire, la transmittance T est une matrice homogène 30 à une image I,IR,IP, soit une matrice de mêmes dimensions qu'une image I,IR,IP et dont chaque élément est une transmittance élémentaire t du pixel de l'image I,IR,IP associé. Une transmittance élémentaire t comprend autant de composantes qu'un pixel : pour une image couleur un pixel 35 comprend typiquement trois composantes dans un système colorimétrique donné (XYZ, RGB, Lab, etc.), la transmittance élémentaire t associée comprend alors trois composantes associées correspondantes ; pour une image monochrome un 3054354 10 pixel comprend une unique composante scalaire, la transmittance élémentaire t associée comprend alors une unique composante scalaire associée. La transmittance élémentaire t représente la fraction de flux lumineux 5 traversant le pixel associé. La transmittance est encore appelée transparence. Un tel découpage faisant apparaître la grandeur transmittance T est avantageux en ce qu'il donne un sens physique à la modélisation et surtout en ce qu'il permet 10 d'utiliser une modélisation existante, faisant intervenir la transmittance T. Ainsi, il est possible d'utiliser, par exemple, un des nombreux modèles prédictifs existants déterminant une réflectance R à partir de la transmittance T. Pour mémoire, la réflectance R est une matrice homogène à 15 une image I,IR,IP, soit une matrice de mêmes dimensions qu'une image I,IR,IP et dont chaque élément est une réflectance élémentaire r du pixel de l'image I,IR,IP associé. Une réflectance élémentaire r comprend autant de composantes qu'un pixel : pour une image couleur un pixel 20 comprend typiquement trois composantes, la réflectance élémentaire r associée comprend alors trois composantes associées correspondantes ; pour une image monochrome un pixel comprend une unique composante scalaire, la réflectance élémentaire r associée comprend alors une unique composante 25 scalaire associée. La réflectance élémentaire r représente la fraction de flux lumineux réfléchie par le pixel associé. On sait d'après les lois de l'optique qu'un flux lumineux rencontrant un matériau transparent est en partie transmis, en partie réfléchi et en partie absorbé. De la conservation 30 des énergies, il découle que la somme de ces trois phénomènes est constante. La transmittance t, la réflectance r et l'absorbance a sont liées par la relation t + r + a = 1. L'utilisation d'un modèle prédictif, déterminant la réflectance R en fonction de la transmittance T, permet 35 d'écrire la fonction de transfert IP = F(S) sous la forme d'une combinaison de trois fonctions, soit IP = F3(F2(F1(S))). Ceci permet de décomposer la fonction de transfert F en trois sous fonctions/modèles : une première 3054354 11 fonction Fl de transmittance transformant une consigne de production S en une transmittance T, T = Fl(S), une deuxième fonction F2 de réflectance transformant une transmittance T en une réflectance R, R = F2(T) et une troisième fonction F3 5 d'image transformant une réflectance R en une image prédite IP, IP = F3(R). La première fonction Fl de transmittance est la plus spécifique. Elle dépend directement du dispositif de production d'image. 10 Selon une mode de réalisation elle peut être est définie par une forme paramétrique. Les paramètres sont ensuite déterminés par apprentissage/calibrage au moyen d'un jeu de paires de calibrage, chaque paire comprenant une consigne de production S et la transmittance T résultante obtenue au 15 moyen du dispositif de production d'image. L'homme du métier sait déterminer une forme paramétrique adaptée au phénomène à modéliser. Un exemple illustratif va maintenant être donné dans le cas particulier simplifié d'un dispositif de gravure laser 20 appliqué à la production d'une image monochrome. Dans ce cas, la transmittance T est avantageusement homogène à la consigne de production S. T est une matrice de mêmes dimensions que S, dont chaque élément est un scalaire t. Ces hypothèses simplificatrices permettent de travailler pixel par pixel et 25 de déterminer une fonction t=t(s), où t est la transmittance élémentaire d'un pixel et s est l'intensité laser du point/pixel homologue associé de la consigne de production S, encore noté Sp. p est un indice parcourant les points/pixels d'une matrice : I, IP, IR, IO, T, R, ou encore ici S. La 30 fonction t=t(s) est une fonction scalaire. Elle peut être étendue à l'ensemble des pixels d'une matrice. La fonction t=t(s) peut être déterminée au moyen d'images réalisées IR par le dispositif de production d'image sur la base de consignes de production S étalons. Ces étalons 35 comportent, par exemple, des intensités laser équi réparties dans tout l'intervalle d'intensité, soit 0-255. La transmittance t est ensuite déterminée sur l'image réalisée IR, par exemple par mesure photométrique. Il est ainsi 3054354 12 possible de déterminer une valeur de transmittance t associée à chaque valeur s d'intensité et donc de tabuler la fonction t=t(s). Cette fonction est nécessairement injective et peut aisément être inversée. 5 La figure 4 illustre un exemple d'une telle fonction t(s). Pour une intensité s nulle, le support 1 reste transparent et une transmittance t maximale est obtenue. Pour une intensité s maximale, le substrat est opacifié et la transmittance t est minimale. Entre ces deux points la 10 fonction est strictement monotone. De nombreux modèles prédictifs, déterminant la réflectance R en fonction de la transmittance T, existent et sont décrits dans la littérature : réflexion diffuse de la lumière, réflexion surfacique de Fresnel, diffusion latérale 15 de la lumière, gain ponctuel physique et optique, modèle de Yule-Nielsen, modèle spectral de Neugebauer, modèle de Ruckdeshel et Hauser, modèle de Clapper-Yule, pour ne citer que les plus connus. Ils peuvent tous être utilisés pour réaliser la deuxième fonction F2. 20 Un mode de réalisation préféré est proposé et décrit sur la base d'un modèle de Clapper-Yule « faible diffraction » modifié. Selon ce modèle, la deuxième fonction F2 est donnée par la formule : b.k1.T2 (1- b).1er(T *M).T R = F2(T)= + où : 1- ii-J2 1- k2.(T2*M) 25 R est la matrice de réflectance, bE[0,1] est un paramètre scalaire, k1,k2E[0,1] sont deux constantes optiques, T est la matrice de transmittance, * est l'opérateur de convolution et 30 M est un noyau de convolution. Ce modèle permet de calculer une réflectance R en fonction d'une transmittance T. Ce modèle est modifié, relativement au modèle de Clapper-Yule, notamment par l'adjonction de la convolution et d'un noyau de convolution 35 M. Ce noyau M permet avantageusement d'introduire explicitement une modélisation du phénomène d'ombrage caractéristique du problème traité. 3054354 13 L'opération de convolution * est bien connue dans le domaine du traitement d'image. Pour mémoire, la convolution A*N d'une matrice A avec un noyau N, matrice carré de dimension impaire, produit une matrice B de mêmes dimensions 5 que A. Un élément bi,j de B est défini par = (A * N);,1 = k . Selon une autre caractéristique, la deuxième fonction F2 comprend encore un opérateur K. Cet opérateur K est destiné à modéliser la fonction de transfert d'un scanner utilisé lors 10 de l'acquisition de l'image réalisée IR, et que l'on souhaite comparer avec une image prédite IP. Cet opérateur peut comprendre toute modification de l'image telle que : mise à l'échelle, interpolation, linéarisation, etc. 15 Selon un mode de réalisation, la deuxième fonction F2 est définie par la formule : R = F 2(T) - b.k1.T22 + (1- (T * M) * K où : 1- k2.T 1- k2.(T2 *M) R, b, kl, k2, T, M désignent les mêmes éléments que précédemment, 20 * est l'opérateur de convolution et K est un noyau de convolution. L'opérateur K est ici exprimé sous la forme d'une convolution avec un noyau de convolution K. Selon un mode de réalisation possible, l'opérateur K 25 modélise le scanner par un simple bruit K. Selon un mode de réalisation ce bruit est, par exemple, un bruit bidimensionnel gaussien de déviation standard o. L'homme du métier sait déterminer un tel noyau K. Selon un mode de réalisation possible le1=a-rj.rd(1-ri) et 30 k2=rg.ri, avec rs indice de réflexion spéculaire de l'interface air image, ri indice de réflectance interne de Fresnel de l'interface air image et rg indice de réflectance intrinsèque d'un substrat réfléchissant lambertien. Ces trois indices sont connus ou déterminable par l'homme du métier. Les deux 35 premiers sont connus : rs =0,04, rf =0,6 . rg peut être déterminé 3054354 14 par la formule rg - rf.1?',+(1-rs).(1-rf)r où Rw est la réflectance du support d'image 1 blanc, soit le support d'image 1 avant toute opération de production d'image. rg dépend du matériau. Rw peut être obtenue par mesure de la réflectance sur un 5 support d'image 1 vierge. Il a été vu que le noyau de convolution M peut avantageusement être utilisé pour modéliser l'effet d'ombrage. La convolution est particulièrement adaptée en ce qu'elle permet d'introduire une influence d'un point sur son 10 voisinage. Or, le phénomène d'ombrage est exactement cela : chaque point projette une ombre sur ses voisins. L'homme du métier sait définir ou modifier le contenu du noyau M afin de modéliser au mieux le phénomène physique d'ombrage. 15 Selon un possible mode de réalisation simple, l'ombrage peut être considéré introduire un flou. Un noyau de convolution M possible pour réaliser un flou est, de manière connue, une matrice carrée unité [1], soit une matrice carrée ne contenant que des 1. Afin que la convolution soit 20 unitaire, et conserve la luminance, la norme du noyau M est avantageusement normalisée à 1. Soit M = [1/d2]cixd, avec d la dimension du noyau M. Un tel noyau unité normalisé réalise, par l'effet de la convolution, un moyennage spatial d'un point avec son voisinage. Un noyau de convolution est une 25 matrice carrée de dimension dxd, avec préférentiellement d impair. La valeur de d détermine ici la taille du voisinage impacté par l'effet d'ombrage. Selon une caractéristique, la dimension d est un paramètre, pris comme inconnue, et est déterminée lors du 30 calibrage de la fonction de transfert F. La troisième fonction F3 d'image permet de transformer une réflectance R en une image prédite IP, IP = F3(R). La réflectance R est la quantité de lumière réfléchie lorsque le support 1 est éclairé par un éclairage multidirectionnel 3. 35 Aussi la matrice R est, sous certaines hypothèses, proportionnelle à l'image. Selon une caractéristique, en considérant que les facteurs de proportionnalité sont égaux à 3054354 15 1, la fonction F3 est la fonction identité. Soit IP = F3(R) = R. Il est ainsi possible de construire une fonction de transfert F. Cette fonction de transfert F calcule une image 5 prédite IP en fonction d'une consigne de production S. Selon les hypothèses retenues pour F, certains paramètres restent à déterminer. Ainsi, selon les hypothèses illustratives indiquées précédemment, la fonction de transfert F est déterminée par 10 la connaissance des paramètres b et d. Si un opérateur K est de plus considéré, il faut aussi connaitre les paramètres définissant cet opérateur K, soit dans le cas illustratif d'un bruit gaussien, sa déviation standard o. Avec une autre modélisation de F, des paramètres 15 différents pourraient devoir être déterminés. La méthode de détermination par calibrage/apprentissage reste cependant identique. Selon une caractéristique, la détermination des paramètres définissant la fonction de transfert F est 20 réalisée par calibrage, ou apprentissage. Ceci permet de résoudre, les paramètres, par exemple b, d et éventuellement la déviation standard o, étant les inconnues, une équation de la forme : paramètres = a rg min A(IRp,IPp) = a rg min A(IRp, F (S p)) 25 Cette équation est une fonction objectif. La résolution consiste, de manière connue, à déterminer les valeurs optimales des paramètres mentionnés à gauche afin de minimiser l'expression à droite. Cette résolution, de par le nombre d'inconnues important, est le plus souvent réalisée 30 par un procédé itératif, mis en oeuvre automatique au moyen d'un solveur informatique. Le principe est de faire varier les paramètres recherchés, indiqués à gauche, dans l'expression de droite, jusqu'à obtenir les valeurs optimales des paramètres qui minimisent l'expression. 35 Dans le cas présent l'expression à minimiser compare une image réalisée IR, réalisée physiquement, obtenue par le dispositif de production d'image en utilisant la consigne de 3054354 16 production S, avec une image prédite IP, obtenue au moyen de la fonction de transfert F à partir de cette même consigne de production S, IP = F(S). Ces deux images IR et IP sont comparées pixel à pixel. IRp est un pixel de l'image réalisée 5 IR, telle qu'observée ou scannée. IPp ou F(Sp) est un pixel de l'image prédite IP par la fonction de transfert F à partir de la consigne de production S. p est un indice parcourant tous les pixels de chacune des deux images réalisée IR ou prédite IP. Chaque pixel d'une image est comparé avec le 10 pixel homologue de l'autre image au moyen d'un opérateur de différence/distance A fournissant une valeur positive ou nulle indicative de la différence. La fonction objective à minimiser somme toutes ces différences, pour l'ensemble des pixels. 15 Selon un mode de réalisation, l'opérateur A est une distance colorimétrique. L'homme du métier connait de nombreuses formulations de telles distances colorimétriques. Dans le cas d'images monochromes, chaque pixel étant scalaire, l'opérateur A peut se résumer à une valeur absolue 20 de la différence. Ceci permet de déterminer les paramètres et ainsi de définir complètement la fonction de transfert F. Pour un tel mode de réalisation il convient d'alimenter la résolution par des paires de calibrage. Chaque paire de 25 calibrage comprend une consigne de production S et une image réalisée IR correspondante obtenue par le procédé initial de production au moyen de la consigne de production S. Ceci nécessite de réaliser physiquement une série d'images IR à partir d'une série de consignes de production S. Chaque image 30 réalisée IR est scannée afin de pouvoir être comparée avec l'image prédite IP = F(S) correspondante, issue de la même consigne de production S. Le nombre de telles paires de calibrage dépend du nombre de paramètres inconnus à déterminer. 35 Dans l'exemple illustratif, où F dépend de b, d et e, un jeu de quelques dizaines de paires de calibrage est suffisant pour calibrer la fonction de transfert F et déterminer les valeurs optimales de chacun des paramètres b, d et o. Ce jeu 3054354 17 de calibrage comprend par exemple : 20 ou 30 consignes de production uniformes comprenant chacune une unique intensité laser, ces intensités étant variables de manière équitablement répartie entre l'intensité nulle et l'intensité 5 maximale, d'autres consignes de production comprenant des motifs périodiques de fréquences variables en ligne et/ou en colonne, ou encore des consignes de production comprenant des motifs indicatifs de l'application. Une fois la fonction de transfert F déterminée, une image 10 prédite IP = F(S) obtenue à partir d'une consigne de production S donnée peut être calculée. En comparant cette image prédite IP avec l'image objectif IO et en tentant de minimiser leur distance, il est possible d'optimiser une nouvelle fonction objectif et de déterminer une consigne de 15 production améliorée SA. Ceci peut être réalisé en résolvant une équation de la forme : SA = a rg min y (A(Iop,ipp))2 = a rg min (A(I0 p, F (S p)))2 , où S S est une consigne de production, 20 IO est une image objective, IP est une prédite obtenue à partir de S, soit IP = F(S), p est un indice parcourant les pixels, F est la fonction de transfert et A est une distance colorimétrique entre deux pixels. 25 Un exemple de distance colorimétrique, pouvant être utilisé ici et/ou dans les résolutions précédentes, est la distance colorimétrique AE94 tel que définie par la commission internationale de l'éclairage (CIE). 30 REVENDICATIONS1. Procédé de production d'une image objective IO au moyen d'un dispositif de production d'image apte à produire une juxtaposition d'éléments d'image (2) dans un substrat (1) présentant une épaisseur et produisant une image en épaisseur telle que se produise un effet d'ombrage de chaque élément d'image (2) sur ses éléments d'image voisins, comprenant les étapes suivantes : - détermination d'une consigne de production S à partir de l'image objective IO, - production d'une image réalisée IR au moyen du dispositif de production d'image utilisant ladite consigne de production S, caractérisé en ce que l'étape de détermination remplace la consigne de production S par une consigne de production améliorée SA déterminée par les étapes suivantes : - modélisation du dispositif de production d'image par une fonction de transfert F transformant une consigne de 20 production S en une image prédite IP, - détermination d'une consigne de production améliorée SA dont l'image prédite IP par la fonction de transfert F, IP = F(SA), est la plus proche possible de l'image objective IO, de telle manière à réduire les conséquences de l'effet 25 d'ombrage dans l'image réalisée IR. 2. Procédé selon la revendication 1, où l'effet d'ombrage se produit lorsque l'épaisseur du substrat (1) de l'image est au moins égale au double de la dimension transverse d'un 30 élément d'image (2). 3. Procédé selon l'une quelconque des revendications 1 ou 2, où le dispositif de production d'image est un dispositif de gravure laser apte à noircir un substrat (1) transparent 35 dans son épaisseur et où une consigne de production S est une carte de tir laser. 4. Procédé selon l'une quelconque des revendications 1 ou 3054354 19 2, où le dispositif de production d'image est un dispositif de micro perforation apte à réaliser des micro trous dans l'épaisseur d'un substrat transparent et où une consigne de production S est une carte de perforation. 5 5. Procédé selon l'une quelconque des revendications 1 à 4, où la modélisation du dispositif de production d'image détermine une transmittance T de la consigne de production S comme résultat intermédiaire. 10 6. Procédé selon l'une quelconque des revendications 1 à 5, où la fonction de transfert F combine une première fonction Fl de transmittance transformant une consigne de production S en une transmittance T, T = Fl(S), une deuxième fonction F2 15 de réflectance transformant une transmittance T en une réflectance R, R = F2(T) et une troisième fonction F3 d'image transformant une réflectance R en une image prédite IP, IP = F3(R), soit IP = F(S) = F3(F2(F1(S))). 7. Procédé selon la revendication 6, où la première fonction Fl est définie par une forme paramétrique, et les paramètres sont déterminés par calibrage. 8. Procédé selon l'une quelconque des revendications 6 ou 25 7, où la première fonction Fl est déterminée empiriquement, et tabulée. 9. Procédé selon l'une quelconque des revendications 6 à 8, où la deuxième fonction F2 est donnée par la formule 30 suivante : 1/412 a_ bg * AI ).T R = F 2(T) - + avec 1- k2/2 1- k2.(T2*M) R une matrice de réflectance, bE[0,1] un paramètre scalaire, k1,lc2Eh1] deux constantes optiques, 35 T une matrice de transmittance, * l'opérateur de convolution et M un noyau de convolution. 3054354 20 10. Procédé selon la revendication 9, où la deuxième fonction F2 inclut de plus un bruit K, selon la formule : R=F2(T)-(b.k1.T2 +(1-b).44r*AlT)*K. , avec 1-k2.T2 1-k2.(r*A1-) 5 * l'opérateur de convolution et K un noyau de convolution. 11. Procédé selon la revendication 10, où K est un noyau bidimensionnel gaussien de déviation standard o. 10 12. Procédé selon l'une quelconque des revendications 9 à 11, où k1-(11-ri) et k2=rg.ri, avec 7-, indice de réflexion spéculaire de l'interface air image, ri indice de réflectance interne de Fresnel de l'interface air 15 image et rg indice de réflectance intrinsèque d'un substrat réfléchissant lambertien. 13. Procédé selon l'une quelconque des revendications 9 à 12, où le noyau de convolution M est une matrice carrée unité 20 de dimensions dxd, normalisée, M = [1/d2]cixd. 14. Procédé selon l'une quelconque des revendications 6 à 13, où la troisième fonction F3 est la fonction identité. 15. Procédé selon l'une quelconque des revendications 6 à 14, où la fonction de transfert F est calibrée en identifiant les paramètres b, d et éventuellement la déviation standard o, par résolution de l'équation suivante : b, d, a = a rg min = a rg min A(IR p ,F (S P I avec 0 11 `74 bERO p d> d>0 &a0 S une consigne de production, IR une image réalisée obtenue au moyen du dispositif de production utilisant la consigne de production S, IP une image prédite par la fonction de transfert F à partir de la consigne de production S, soit IP = F(S), p un indice parcourant les pixels, 3054354 21 F la fonction de transfert et A une distance colorimétrique entre deux pixels. 16. Procédé selon la revendication 15, où la résolution est 5 alimentée par des paires de calibrage comprenant une consigne de production S et une image réalisée IR au moyen du dispositif de production utilisant la consigne de production S. 10 17. Procédé selon l'une quelconque des revendications 1 à 16, où l'étape de détermination d'une consigne de production améliorée SA comprend un calcul itératif automatique minimisant une distance entre l'image objectif IO et une image prédite IP par la fonction de transfert F appliquée à 15 la consigne de production améliorée SA. 18. Procédé selon l'une quelconque des revendications 1 à 17, où l'étape de détermination d'une consigne de production améliorée SA comprend la résolution de l'équation suivante : 20 SA= a rg min y (A(Io,ipp))2 a rg min (0(I0 p , F (S p)))2 , avec S 71;i g S une consigne de production, IO une image objective, IP une image prédite par la fonction de transfert à partir de la consigne de production S, soit IP = F(S), 25 p un indice parcourant les pixels, F la fonction de transfert et 8 une distance colorimétrique entre deux pixels. 19. Image produite par un procédé selon l'une quelconque des 30 revendications précédentes. 20. Dispositif de production d'image apte à produire une juxtaposition d'éléments d'image (2) dans un substrat (1) présentant une épaisseur et produisant une image en épaisseur 35 telle que se produise un effet d'ombrage de chaque élément d'image (2) sur ses éléments d'image voisins, comprenant : - un moyen d'acquisition d'une image objectif IO, 3054354 22 - un moyen de détermination d'une consigne de production S à partir de l'image objective IO, - un dispositif de production d'une image réalisée IR à partir de ladite consigne de production S, 5 caractérisé en ce que le moyen de détermination comprend encore des moyens de mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 18 de sorte à déterminer une consigne de production améliorée SA pour remplacer la consigne de production S.