FR 2980000 A1 20130315 FR 1158123 A 20110913 PROCEDE DE CORRECTION DES EFFETS DE PROXIMITE ELECTRONIQUE UTILISANT UNE DECONVOLUTION DU MOTIF A INSOLER PAR UNE METHODE PROBABILISTE La présente invention s'applique notamment au domaine de la lithographie électronique pour la gravure directe sur wafer ou pour la fabrication de masques. Aux dimensions critiques inférieures à 22 nm atteintes par les technologies en production, l'interaction du faisceau d'électrons avec la cible est notamment affectée par une diffusion des électrons autour de la trajectoire initiale (effet de forward scattering ou de diffusion vers l'avant) ainsi que par une rétrodiffusion (effet de back scattering ou de diffusion vers l'arrière). Ces effets, qualifiés d'effets de proximité, dépendent notamment des matériaux de la cible et de sa géométrie. Il est donc nécessaire de tenir compte des effets de proximité pour obtenir des motifs gravés fidèles au design défini. Pour cela, il est connu d'appliquer deux méthodes, en général de manière indépendante. D'une part, on modélise les effets de proximité pour moduler les doses de rayonnement électronique utilisées pour bombarder la cible. D'autre part, on modifie le design source en fonction de règles de dessin élaborées de manière plus ou moins empirique pour que le design cible soit conforme à l'objectif recherché. En ce qui concerne la modulation de dose, l'état de l'art de référence est constitué par la publication de Eisennmann (« Proxecco - Proximity Effect Correction by Convolution » - H. Eisennmann, T. Waas, H. Hartmann - J. Vac. Sci. Technol., vol B11, N°6, pp 2741-2745, Nov-Dec. 1993) et le logiciel commercialisé sous la marque ProxeccoTM. La modulation de dose est appliquée à un motif fracturé, chaque élément de base recevant une dose modulée qui est générée à partir d'un modèle des effets de proximité. Ledit modèle utilise généralement une fonction dite de diffusion ou d'étalement de point (ou PSF, point spread function). La PSF de Eisennmann est une paire de gaussiennes, l'une pour la diffusion vers l'avant et l'autre pour la diffusion vers l'arrière. Les doses locales sont déterminées par convolution de la PSF avec la géométrie de la cible. Celle-ci est le cas échéant corrigée selon un principe analogue à ceux de l'Optical Proximity Correction (OPC) utilisée en lithographie optique. La 35 correction géométrique à appliquer est en général déterminée par une optimisation automatique des dimensions des formes : une ligne sera par exemple amincie, épaissie, allongée ou rétrécie. Un dessinateur peut également proposer des géométries alternatives qui lui sont suggérées par son savoir-faire. Les corrections ainsi déterminées pour un type de motif permettent de générer des règles de correction qui sont ensuite appliquées automatiquement sur tous les points de la puce qui ont la même configuration. Pour un exposé détaillé des différentes modalités de correction géométrique on pourra se reporter notamment aux publications de Sewell et de Cook (« Control of Pattern Dimensions in Electron Lithography » - H. Sewell, J. Vac. Sci. Technol., 15, 927 (1978) ; « Dose, Shape and Hybrid Modifications for PYRAMID in Electron Beam Proximity Effect Correction » - Brian D. Cook, Soo Young Lee - IEEE transactions on semiconductor manufacturing, vol. 11, N°1, Feb. 1998). Il est également connu de décomposer l'exposition en plusieurs passes. L'étape d'exposition supplémentaire s'applique alors en général à des géométries plus fines afin notamment d'améliorer la résolution et la fenêtre de procédé. On parle alors d'étape de procédé ou de fonction eRIF (electron Resolution Improvement Feature). Des procédés de ce type ont été décrits notamment par Martin et Manakli (« New writing strategy in electron beam direct write lithography to improve critical dense fines patterning for sub45nm nodes » - L. Martin - EMLC 2009 ; demande internationale de brevet n° PCT/EP2011/055746 ayant pour titre « Procédé de lithographie électronique avec correction des arrondissements de coins », dont S. Manakli est inventeur et dont le titulaire est un des co-déposants de la présente demande de brevet. Le calcul séparé des paramètres de ces différents procédés conduit cependant à une solution sous-optimale à la fois en temps de conception et en temps d'exposition. Des solutions de calcul conjoint des paramètres de modulation de dose et de correction géométrique ont déjà été proposées. Une voie est celle proposée par la demande internationale de brevet n° PCT/EP2011/055863 ayant pour titre « Procédé de lithographie à optimisation combinée de l'énergie rayonnée et de la géométrie de dessin », dont S. Manakli est inventeur et dont le titulaire est un des co-déposants de la présente demande de brevet. Selon ce procédé, l'énergie rayonnée dans une dose appliquée à une zone et les dimensions du motif à graver sur cette zone sont calculées en combinaison par l'intermédiaire de la latitude en énergie du procédé. Ce procédé marche bien pour des motifs simples tels que des cellules à grande maille, des fins de lignes et des coins, notamment. Il est moins efficace pour des motifs plus complexes. Une autre voie d'optimisation conjointe des paramètres de dose rayonnée et de géométrie projetée consiste à minimiser l'écart entre le résultat de la convolution d'une dose rayonnée avec la PSF et le motif cible. Pour ce faire, il a été suggéré de déconvoluer le motif à graver par une méthode appropriée. L'utilisation de transformées de Fourier inverses en combinaison avec une PSF classique à deux gaussiennes telle que proposée par Eisennmann (référence citée ci-dessus) a été préconisée par Haslam ( Cependant, aux dimensions pertinentes pour les technologies actuellement en production ou en développement (dimension critique de une à deux dizaines de nanomètres), cette méthode n'est plus adaptée en raison de la coupure qu'elle réalise des hautes fréquences spatiales, inhérente à la l'inversion, qui empêche une prise en compte suffisamment précise des effets de diffusion vers l'avant qui dominent à ces distances. L'invention résout ce problème en proposant de combiner une PSF d'un type particulier avec une déconvolution du motif à graver utilisant une méthode itérative de type probabiliste. A cet effet, l'invention prévoit un procédé de lithographie par rayonnement d'au moins un motif à insoler sur un support résiné comprenant au moins une étape de génération d'une fonction choisie d'étalement de point dudit rayonnement et une étape d'application d'une méthode choisie de déconvolution audit motif à insoler par ladite fonction d'étalement de point , ledit procédé étant caractérisé en ce que ladite fonction d'étalement de point est choisie pour ne modéliser que les effets de diffusion dudit rayonnement vers l'avant et en ce que ladite méthode de déconvolution est choisie parmi les méthodes modélisant les interactions dudit rayonnement avec le support résiné par une distribution de probabilité conjointe des dites interactions. Avantageusement, ladite méthode de déconvolution comprend une étape de calcul d'un maximum de vraisemblance de chaque point source dudit rayonnement pour tous les points image du motif à insoler, chaque point image étant multiplié par la valeur de la fonction d'étalement de point entre le point source et le point image. Avantageusement, ladite méthode de déconvolution est une méthode de Lucy-Richardson. Avantageusement, le procédé de lithographie de l'invention comprend en outre une étape de tramage de l'image déconvoluée du motif à insoler. Avantageusement, le procédé de lithographie de l'invention comprend en 15 outre une étape de découpage de l'image déconvoluée en au moins deux niveaux distincts d'exposition. Avantageusement, les au moins deux niveaux distincts d'exposition sont choisis pour maximiser la la population en pixels de chacun des niveaux. Avantageusement, le procédé de lithographie de l'invention comprend en outre une étape de vectorisation de l'image déconvoluée du motif à insoler. Avantageusement, la fonction d'étalement de point est une gaussienne. Avantageusement, la largeur à mi-hauteur de la gaussienne est comprise entre 3nm et 15nm. Pour mettre en oeuvre le procédé, l'invention prévoit également un 30 programme d'ordinateur comprenant des instructions de code de programme permettant l'exécution du procédé de lithographie par rayonnement de l'invention lorsque le programme est exécuté sur un ordinateur, ledit programme comprenant une fonction de génération d'une fonction choisie d'étalement de point dudit rayonnement et un module d'application d'une 35 méthode choisie de déconvolution aux sorties éventuellement filtrées de 20 25 ladite fonction d'étalement des points dudit rayonnement appliquée audit motif à insoler, ledit programme étant caractérisé en ce que ladite fonction d'étalement de point est choisie pour ne modéliser que les effets de diffusion dudit rayonnement vers l'avant et en ce que ladite méthode de déconvolution est choisie parmi les méthodes modélisant les interactions des points dudit rayonnement par une distribution de probabilité conjointe des dites interactions. Avantageusement, le programme d'ordinateur de l'invention comprend en outre une fenêtre d'interface dans laquelle un opérateur peut programmer l'application à au moins un deuxième motif appartenant à une classe déterminée de règles de correction déterminées par déconvolution d'un premier motif appartenant à ladite classe. L'invention est particulièrement avantageuse car elle permet de parvenir de manière complètement automatique à une image rayonnée optimale. Les limitations du procédé ne dépendent plus que des limitations de la machine de lithographie et des arbitrages à effectuer pour favoriser le rendement global du procédé. Ce niveau de résultat n'aurait été accessible auparavant, pour les types de motifs complexes qui peuvent être ainsi traités, qu'au prix d'une modulation de dose de très grande précision et de corrections géométriques très sophistiquées nécessitant une longue mise au point. En outre, le procédé de l'invention génère une image déconvoluée en niveaux de gris qui pourra subir des traitements différents selon que l'on 25 utilise une machine d'exposition à matrice de points ou une machine de projection à faisceau formé (de type Variable Shaped Beam ou VSB). Dans ce deuxième cas, il est d'ailleurs possible de générer automatiquement des eRIF qui permettront d'améliorer la résolution et la vitesse d'écriture, ce qui constitue un avantage supplémentaire. 30 L'invention sera mieux comprise, ses différentes caractéristiques et avantages ressortiront de la description qui suit de plusieurs exemples de réalisation et de ses figures annexées dont : - La figure 1 illustre l'application de doses modulées selon l'art 35 antérieur; - La figure 2 représente un organigramme des traitements à appliquer pour mettre en oeuvre l'invention dans plusieurs de ses modes de réalisation ; - Les figures 3a et 3b représentent respectivement un motif à exposer et sa PSF associée, selon un mode de réalisation de l'invention ; - Les figures 4a et 4b représentent respectivement une image en niveaux de gris d'un motif à exposer et sa transformée en image binaire par tramage, selon un mode de réalisation de l'invention ; - Les figures 5a et 5b représentent respectivement une image en mode bitmap d'un motif à exposer et sa transformée en image vectorielle, selon un mode de réalisation de l'invention ; - Les figures 6a, 6b et 6c représentent trois niveaux d'exposition différents, selon un mode de réalisation de l'invention ; La figure 1 illustre l'application de doses modulées selon l'art antérieur. On constate que l'exposition est corrigée en fonction de la position des motifs. Les doses reçues sont plus élevées en bordure de plaquette qu'au centre. En effet, les motifs situés à l'intérieur reçoivent en proportion des électrons rétrodiffusés en provenance des zones contigües plus nombreux que les zones situées en bordure.. Selon Eisennmann (référence citée ci-dessus), ces modulations de doses sont calculées à partir de la PSF f(x, y) à double gaussienne suivante : / Ej Eji. 1 1 o f (0) E 3(110) mea ^ Red ) Avec les notations : o a est la largeur du rayonnement direct; o 6 est la largeur de rétrodiffusion ; o ri est le ratio des intensités des rayonnements direct et rétrodiffusé. o E est la position radiale d'un point (x, y). Les valeurs des paramètres a, 6 et ri peuvent être déterminées expérimentalement pour un procédé donné. Ces paramètres sont fonction de la tension d'accélération de la machine et de la cible. Selon Eisennmann, pour effectuer la correction illustrée par la figure 1, il suffit de traiter le phénomène à rayon d'action élevé (supérieur au micron) par rapport aux points d'impact, soit l'effet de diffusion vers l'arrière, et on peut ne pas tenir compte des effets de diffusion vers l'avant qui sont dominants à faible distance (rayon d'action d'un ordre de grandeur plus faible pour les basses tensions d'accélération (de l'ordre de 5kV) et deux ordres de grandeur pour les accélérations plus fortes (typiquement 50kV)) . Il est possible d'effectuer une déconvolution de l'effet de la PSF par une transformée de Fourier inverse: / p E Fm F(r) F(PSF)) Il est nécessaire de réaliser une troncature des fréquences spatiales élevées pour éviter l'annulation du dénominateur. Ces fréquences étant représentatives de la diffusion vers l'avant, cela n'est pas gênant dans l'application décrite par Eisennmann. L'image résultante r est reconstituée comme la somme d'un terme égal à cette déconvolution tronquée des hautes fréquences et d'un terme proportionnel à l'image source p : r E Fm _g con, E F (p)_E C E p Ainsi la partie du motif proche du centre du faisceau est-elle bien réintroduite 20 dans l'image résultante, mais sans correction. Cette méthode de déconvolution n'est pas applicable lorsque l'on s'intéresse précisément aux fréquences spatiales élevées, ce qui est le cas de la présente invention. L'invention propose une méthode permettant de surmonter cette limitation. 25 La figure 2 représente un organigramme des traitements à appliquer pour mettre en oeuvre l'invention dans plusieurs de ses modes de réalisation. Les premières étapes 210, 220 du procédé de l'invention sont indépendantes de l'outil d'exposition utilisé. Il s'agit d'abord de choisir le type de PSF utilisé. 30 Puisque l'on s'intéresse essentiellement aux phénomènes dominants à une échelle de l'ordre de la dizaine de nanomètres, il n'est pas nécessaire de modéliser dans la PSF l'effet de diffusion vers l'arrière. On choisira donc en règle générale au cours de l'étape 210 une PSF qui ne modélise que l'effet de diffusion vers l'avant, qui génère l'essentiel des effets de proximité à l'échelle à laquelle on se situe. Avantageusement, on choisira pour la PSF une gaussienne dont la largeur à mi-hauteur sera comprise entre 3 nm et 15nm. D'autres fonctions peuvent néanmoins convenir, notamment des fonctions symétriques de largeur à la base comprise typiquement entre 1 Onm et 50 nm. Une fois le type de PSF choisi, il s'agit d'appliquer au cours d'une étape 220 une méthode de déconvolution adaptée. Ladite méthode doit à la fois être adaptée à la fonction à déconvoluer, être convergente et réaliser cette convergence dans un temps admissible. Dans le cadre de la présente invention, on se propose d'utiliser des méthodes appartenant à la classe des méthodes probabilistes de déconvolution. Les méthodes de ce type utilisant des itérations sur une variable aléatoire conditionnelle semblent les plus adaptées. En effet, par l'effet de diffusion, la dose reçue en un point dépendra des doses rayonnées 15 en tous les points voisins. Il s'agit donc de calculer en tout point j la dose pià générer telle qu'en un point i la dose reçue effective soit 1-1, étant entendu que la dose pi rayonnera en d'autres points que i .selon une distribution donnée par la PSF que l'on peut discrétiser en un ensemble de valeurs ui,j . Il convient donc de résoudre le problème : 20 E p Une des méthodes testées, dont la convergence a été confirmée, est l'algorithme de Lucy-Richardson (« Bayesian-Based Iterative Method of Image Restoration » - W. Hadley-Richardson, Visibility Laboratory, University of California at San Diego - Journal of the Optical Society of America, Vol. 25 62, N°1, Jan. 1972). On suppose que pi suit une loi donnée (par exemple une loi de Poisson), et on recherche pour tout j la valeur la plus probable pour un ri connaissant par ailleurs ui 1 . Cela revient à résoudre l'équation récursive suivante : 30 p(toi) p(t) v Lu " Avec ci ^ lp,,u(jt) L'algorithme converge donc, bien que lentement. D'autres méthodes de la même famille, notamment celles dites « EM », pour Expectation - Maximization, basées sur une alternance d'étapes de calcul d'espérance mathématique et de maxima de vraisemblance, pourront être utilisées sans sortir du cadre de l'invention. L'image déconvoluée est obtenue en niveaux de gris, chaque niveau en chaque point codant l'intensité de la dose à rayonner. Une opération unique permet donc de réaliser à la fois la modulation de dose et la correction géométrique toutes deux optimales pour le motif à exposer. Cette opération unique s'intègre donc en lieu et place de ces deux opérations dans les procédés de l'art antérieur de préparation des données (dataprep) implantés dans les outils de lithographie électronique du commerce, sans modification autre que logicielle. Il peut cependant être également avantageux d'optimiser l'architecture de calcul de manière adaptée à la méthode de déconvolution choisie. Une architecture de calcul adaptée sera préférentiellement de nature hétérogène, c'est à dire contenir des processeurs standards et des unités de calculs particulières, par exemple des GPUs (Graphics Processing Unit, puces initialement dédiées au graphisme, mais qui ont aujourd'hui évolué pour développer une puissance de calcul vectoriel inédite) ou des FPGA (puces dont le circuit logique est programmable : il est donc possible de programmer une puces à l'efficacité optimisée pour un algorithme donné). Il est ensuite nécessaire de convertir cette image en niveaux de gris en format exploitable par le type d'outil d'exposition utilisé (étapes 230a et 230b de la figure 2). Deux types principaux d'outil sont généralement utilisés, des outils à matrice de points, comme celui commercialisé par la société Mapper et des outils à formation de faisceau, comme ceux commercialisé par la société Vistec (Vistec SB250 electron beam lithography system par exemple) Dans le premier cas, la conversion est une opération de tramage. Dans le deuxième cas, il s'agit d'une opération de vectorisation. Ces deux étapes seront décrites en commentaires respectivement aux figures 4a, 4b et 5a, 5b. Dans le cas de la formation de faisceaux, il est avantageux, pour des motifs à exposer complexes, d'interposer une opération 240 de découpage en niveaux permettant l'exposition d'eRIF. Cette opération est décrite en commentaire aux figures 6a, 6b, 6c. Le motif à exposer est ensuite fracturé en zones correspondants chacune à un faisceau (étape 250). Dans les deux procédés correspondant à l'outil à matrice de points et à l'outil 35 à formation de faisceau, l'étape 260 d'exposition peut être alors déclenchée. Les figures 3a et 3b représentent respectivement un motif à exposer et sa PSF associée, selon un mode de réalisation de l'invention. La figure 3a représente un motif à exposer 310a. Ledit motif comprend des 5 structures complexes, ayant chacune des dimensions critiques de l'ordre de la dizaine de nanomètres. La figure 3b représente la PSF choisie 310b. Comme indiqué, dans le cadre de l'invention, on choisit avantageusement une PSF constituée par une fonction gaussienne modélisant l'effet de diffusion vers l'avant, qui domine 10 aux dimensions d'intérêt pour les structures des types de motifs considérés. Les figures 4a et 4b représentent respectivement une image en niveaux de gris d'un motif à exposer et sa transformée en image binaire par tramage, selon un mode de réalisation de l'invention. 15 Un outil de type Mapper accepte en entrée une image binaire (pixel 0 ou 1) : le passage d'image en niveaux de gris à image binaire s'appelle le tramage (dithering). Il existe de nombreux algorithmes pour effectuer cette transformation, le plus connu étant Floyd-Steinberg. C'est un algorithme à diffusion d'erreur. Cet algorithme effectue un tramage par la diffusion de 20 l'erreur de quantification d'un pixel à ses voisins. Par exemple, un premier pourcentage de son erreur choisi en fonction du nombre de pixels est ajouté au pixel à sa droite, un deuxième pourcentage au pixel situé en bas à gauche, un troisième pourcentage au pixel situé en dessous et le complément au pixel en bas à droite. Les pourcentages sont choisis 25 notamment en fonction du nombre de pixels. Les figures 5a et 5b représentent respectivement une image en mode bitmap d'un motif à exposer et sa transformée en image vectorielle, selon un mode de réalisation de l'invention. 30 La vectorisation est une transformation de l'image comme ensemble de points (bitmap) en image dite vectorielle, qui représente des formes géométriques. Des zones sont alors délimitées correspondant à des doses différentes. Les niveaux de gris de la figure 5a sont restreints et leur étendue est transformée en polygones. Plusieurs algorithmes de vectorisation sont 35 utilisables, par exemple Potrace (http://potrace.sourceforganeti, http://potrace.sourceforge.net/potrace.pdf). Les figures 6a, 6b et 6c représentent trois niveaux d'exposition différents, selon un mode de réalisation de l'invention. II est souvent avantageux de réaliser plusieurs passes d'exposition. Ce découpage en niveaux permet de réaliser un moins grand nombre de shots et autorise en conséquence une vitesse d'écriture plus rapide. Pour décider du découpage en niveaux, on peut par exemple avantageusement réaliser un histogramme des valeurs prises par les pixels : plus un niveau de luminosité est peuplé (ie plus le nombre de pixels correspondant à ce niveau est élevé), plus il est judicieux de retenir ce niveau comme niveau d'exposition distinct. Chaque niveau extrait par l'étape ainsi décrite reçoit une dose donnée qui peut être exposée en une seule passe. Sur la figure 6b, les lignes fines sont 15 des eRIF générés automatiquement grâce à l'invention contrairement à l'art antérieur. Il est ensuite possible de réaliser l'étape de fracturation consistant à découper les images vectorielles en formes élémentaires (triangles rectangles et rectangles) de taille adaptée au procédé qui seront exposées 20 par l'outil VSB. Ce procédé de fracturation est bien connu de l'homme du métier de la lithographie électronique. Avantageusement, une aide spécifique à la conception peut être intégrée dans les outils dont dispose l'opérateur. 25 La correction géométrique suppose une action de l'opérateur qui, en fonction de son savoir faire et des résultats des expositions déjà faites a idée des modifications automatiques (raccourcissement/élongation de lignes, etc.) qui doivent êtres appliquées à tout ou partie du circuit. Pour aider cet opérateur lorsque les modifications ne sont pas triviales et 30 qu'une intervention manuelle reste nécessaire, un outil mettant en oeuvre le procédé décrit peut être mis en oeuvre dans son application de design: une fenêtre, ouverte à sa demande, affiche une solution optimale. L'opérateur peut alors choisir de généraliser à des motifs ayant les mêmes caractéristiques les règles de correction qu'il observe dans cette fenêtre, 35 créant ainsi une bibliothèque de règles de correction réutilisables pour des motifs similaires utilisant par exemple des algorithmes de reconnaissance automatique de forme. Les exemples décrits ci-dessus sont donnés à titre d'illustration de certains 5 modes de réalisation de l'invention. Ils ne limitent en aucune manière le champ de l'invention qui est défini par les revendications qui suivent. REVENDICATIONS1. Procédé de lithographie par rayonnement d'au moins un motif à insoler (310a) sur un support résiné comprenant au moins une étape de génération d'une fonction choisie d'étalement de point (310b) dudit rayonnement et une étape d'application (220) d'une méthode choisie de déconvolution audit motif à insoler par ladite fonction d'étalement de point , ledit procédé étant caractérisé en ce que ladite fonction d'étalement de point est choisie pour ne modéliser que les effets de diffusion dudit rayonnement vers l'avant et en ce que ladite méthode de déconvolution est choisie parmi les méthodes modélisant les interactions dudit rayonnement avec le support résiné par une distribution de probabilité conjointe des dites interactions. 2. Procédé de lithographie selon la revendication 1, caractérisé en ce que ladite méthode de déconvolution comprend une étape de calcul d'un maximum de vraisemblance de chaque point source dudit rayonnement pour tous les points image du motif à insoler, chaque point image étant multiplié par la valeur de la fonction d'étalement de point entre le point source et le point image. 3. Procédé de lithographie selon la revendication 2, caractérisé en ce que ladite méthode de déconvolution est une méthode de Lucy-25 Richardson. 4. Procédé de lithographie selon l'une des revendications 1 à 3, caractérisé en ce que qu'il comprend en outre une étape de tramage de l'image déconvoluée du motif à insoler. 30 5. Procédé de lithographie selon l'une des revendications 1 à 3, caractérisé en ce qu'il comprend en outre une étape de découpage de l'image déconvoluée en au moins deux niveaux distincts d'exposition. 6. Procédé de lithographie selon la revendication 5, caractérisé en ce que les au moins deux niveaux distincts d'exposition sont choisis pour maximiser la population en pixels de chacun des niveaux. 7. Procédé de lithographie selon l'une des revendications 5 à 6, caractérisé en ce qu'il comprend en outre une étape de vectorisation de l'image déconvoluée du motif à insoler. 8. Procédé de lithographie selon l'une des revendications 1 à 7, 10 caractérisé en ce que la fonction d'étalement de point est une gaussienne. 9. Procédé de lithographie selon la revendication 8, caractérisé en ce que la largeur à mi-hauteur de la gaussienne est comprise entre 3nm et 15nm. 15 10. Programme d'ordinateur comprenant des instructions de code de programme permettant l'exécution du procédé de lithographie par rayonnement selon l'une des revendications 1 à 9 lorsque le programme est exécuté sur un ordinateur, ledit programme comprenant une fonction de 20 génération d'une fonction choisie d'étalement de point dudit rayonnement et un module d'application d'une méthode choisie de déconvolution aux sorties éventuellement filtrées de ladite fonction d'étalement des points dudit rayonnement appliquée audit motif à insoler, ledit programme étant caractérisé en ce que ladite fonction d'étalement de point est choisie pour ne 25 modéliser que les effets de diffusion dudit rayonnement vers l'avant et en ce que ladite méthode de déconvolution est choisie parmi les méthodes modélisant les interactions des points dudit rayonnement par une distribution de probabilité conjointe des dites interactions. 30 11. Programme d'ordinateur selon la revendication 10, caractérisé en ce qu'il comprend en outre une fenêtre d'interface dans laquelle un opérateur peut programmer l'application à au moins un deuxième motif appartenant à une classe déterminée de règles de correction déterminées par déconvolution d'un premier motif appartenant à ladite classe. 35