FR 2774761 A1 19990813 FR 9801693 A 19980212 PROCEDE ET DISPOSITIF DE DETERMINATION DE LA REPARTITION D'INTENSITE ET DE PHASE DANS DIFFERENTS PLANS DE COUPE D'UN FAISCEAU LASER Domaine technique et art antérieur L'invention concerne un procédé et un dispositif pour la détermination de la répartition d'intensité et/ou de phase dans différents plans de coupe d'un faisceau spatialement cohérent, notamment du faisceau lumineux issu d'un laser. Les dispositifs connus, pour réaliser une telle détermination, se répartissent en trois catégories 1) les dispositifs interférométriques, 2) les dispositifs de type Hartmann-Shack, 3) les dispositifs de prélèvement de profils spatiaux. Les dispositifs et méthodes interférométriques utilisent la combinaison du front d'onde à mesurer avec une onde ayant une relation de phase avec celle-ci. Les dispositifs et méthodes de type Hartmann Shack utilisent un masque de phase (matrices de lentilles, trous...) pour évaluer la pente de la surface d'onde en chaque point. La reconstruction permet alors de calculer la surface d'onde. Le problème de l'obtention de la phase à partir de simples profils d'intensités a déjà été étudié dans la littérature sous le nom anglais de "phase retrieval problem". L'étude de différentes solutions existantes est proposée dans les article "Phase retrieval algorithms : a comparison", Applied Optics, vol. 21, p. 2758, Août 1982 et "The phase retrieval problem", IEEE transactions on Antennas and Propagation AP29, p.386, Mars 1981. On connaît en particulier l'algorithme proposé par Gerchberg et Saxton en 1971 (revue Optik 34, 275 (1971)). Ces auteurs utilisent un plan image, et un plan de diffraction en champ lointain, au foyer d'une lentille. Selon cet algorithme, on mesure un profil d'intensité dans un premier plan de coupe, en attribuant un profil de phase arbitraire en tout point de ce plan ; puis, on calcule, à l'aide des équations de propagation du faisceau, le profil d'intensité et de phase dans un second plan de coupe : on remplace le profil d'intensité calculé dans ce second plan par le profil mesuré, et on calcule à nouveau les profils d'intensité et de phase dans le premier plan. On itère jusqu'à convergence. D'après l'article de D.L. Misell ("A method for the solution of the phase problem in electron microscopy") paru dans J. Phys. D. Appl. Phys., vol. 6, 1973, p. L6-L9 et le commentaire sur cet article de Gerchberg et Saxton (mêmes références, p. L31) les deux plans peuvent être quelconques mais proches l'un de 1 'autre. Dans la littérature récente, de nombreuses études ont été faites sur ces méthodes, appliquées à la microscopie électronique et à l'émission des antennes. Ces études proposent des améliorations algorithmiques. On peut citer par exemple l'article "Radiation Pattern Evaluation from near-field Intensities on planes", IEEE Transactions on Antennas and Propagation, vol. 44, nO 5, Mai 1996, qui propose de ne pas utiliser nécessairement le plan objet et le plan de Fourrier. Des essais effectués avec des images de faisceaux laser et ces méthodes n'ont pas permis d'aboutir à des résultats suffisamment bons pour que l'on soit en mesure de se passer de systèmes de type Hartmann-Shack. Il se pose donc le problème de trouver un procédé et un dispositif permettant d'améliorer la précision de la reconstruction de phase, et de concurrencer des systèmes commerciaux coûteux, par exemple de type Hartmann-Shack ou Zygo. Exposé de l'invention L'invention a tout d'abord pour objet un procédé itératif de détermination de la répartition d'intensité et de phase d'un faisceau lumineux cohérent dans un plan, comportant - la mesure de l'intensité Ii, i=1... N, du faisceau, dans N plans, N > 3, dont le plan dans lequel on veut déterminer la répartition d'intensité et de phase, - le choix, pour le plan i=l, d'une matrice de phase pl de départ, et le calcul d'une matrice d'amplitude complexe par produit terme à terme de la matrice de phase J(Pl avec la matrice d'amplitude A1 correspondante, - pour chaque plan j > 1 : * détermination d'une matrice complexe propagée B'j, à partir de la matrice mesurée d'intensité I l du plan j-l et de la matrice de phase p du plan j-l, * extraction, de B'j, de la matrice de phase du plan i, - itération du procédé jusqu'à convergence (et j=l lorsque j-1=N). Le procédé selon l'invention met donc en oeuvre N plans, avec N > 3. L'un de ces plans est celui dont on veut déterminer la répartition d'intensité et de phase. Le procédé débute avec un profil de phase (nul, par exemple) sur un premier plan. On effectue ensuite, par calcul, des propagations de faisceau d'un plan à un autre en conservant, à l'issue de chaque itération, le profil de phase calculé pour un plan, et en introduisant le profil d'intensité mesuré dans le même plan de coupe. On passe ainsi du plan n" 1 au plan n" 2, ... au plan n" N, puis on revient dans le plan n" 1. On constate la convergence du profil de phase vers le profil recherché : en effet, les profils d'intensité et de phase calculés dans l'un des N plans permettent de retrouver fidèlement les profils d'intensité mesurés dans n'importe quel autre plan. Pour chaque plan j(j > l), on calcule un profil de phase à partir de la matrice complexe traduisant l'état du faisceau dans ce plan. Cette matrice complexe résulte elle-même d'un calcul, qui traduit les équations de propagation, mettant en oeuvre la matrice (mesuré) d'intensité et la matrice (calculée) de phase dans le plan j-i. Autrement dit, pour chaque plan j, on calcule, par propagation, le profil d'intensité et de phase dans ce plan, puis on remplace le profil d'intensité calculé de ce plan par le profil d'intensité mesuré. On peut alors, à partir de ce profil d'intensité mesuré et du profil de phase calculé, effectuer un calcul de propagation pour le plan j+l. On peut utiliser un cheminement quelconque entre les N plans, selon le besoin. Par conséquent, selon l'invention, au lieu d'évoluer entre deux plans de coupe en aller et retour, on met en oeuvre au moins trois plans de coupe, et on passe, par exemple, du premier plan au second, puis au troisième, .... puis au Nième, et on revient au premier plan, et ainsi de suite. La précision de la reconstruction de phase sur les plans intermédiaires est ainsi notablement améliorée, ce qui permet de concurrencer des systèmes commerciaux coûteux par exemple de type Hartmann-Shack ou Zygo. Selon un perfectionnement, la phase prise en compte dans le premier plan de coupe à chaque itération est une moyenne entre les profils de phase déterminés à partir de tous les couples de plans par un allerretour. Cette méthode permet d'introduire un coefficient de pondération dans le calcul de la moyenne, ce coefficient pouvant dépendre d'une fonction de mérite du calcul en cours. On optimise ainsi la convergence sur les N profils. Deux versions du procédé selon l'invention sont donc proposées selon le cheminement du calcul entre les plans de mesure. La première méthode parcourt les plans de manière séquentielle, alors que la seconde effectue des aller-retour entre des couples de plans, et permet de pondérer le poids relatif des différents plans mesurés. L'invention concerne également un dispositif de mise en oeuvre de l'invention. Ainsi, l'invention concerne un dispositif, pour déterminer l'intensité et la phase d'un faisceau lumineux dans un plan comportant - des moyens pour mesurer l'intensité Ii (i=l, . . N) du faisceau dans N plans, N > 3, - des moyens de calcul programmés pour, à partir d'une matrice de phase pl de départ correspondant à un plan * calculer une matrice d'amplitude complexe (B1) par produit de la matrice de phase .icpl avec la matrice d'amplitude (Al) correspondante, * pour chaque plan j > l : calculer une matrice complexe propagée B'j à partir de la matrice mesurée d'amplitude Aj-1 du plan j-l et de la matrice de phase ç du plan j-i, extraire une matrice de phase du plan j de la matrice B'j, itérer le calcul ci-dessus en faisant j=l lorsque j-l=N. Ce dispositif utilise autant de caméras que de plans de coupe ; de préférence les images mesurées sont numérisées avant traitement. Il peut comporter un système de mise en forme du faisceau, dont il limite l'étendue pour l'adaptation aux caméras et à la longueur du banc ; une étape d'étalonnage est alors incorporée dans les étapes du procédé principal. L'adjonction d'une lentille ou de tout autre objet de phase, permettant d'adapter le faisceau au dispositif de mesure (maîtrise de la divergence...) permet d'augmenter le champ d'application du dispositif. L'étape d'étalonnage permet ensuite de s'affranchir, par calcul, de la présence du composant ajouté, et d'accéder ainsi aux caractéristiques du faisceau. Un analyseur de phase utilisant un procédé selon l'invention permet d'obtenir la répartition de l'amplitude complexe du champ, représentant à la fois l'intensité et la phase du faisceau dans un plan donné (en général un plan image). Cette mesure permet alors de calculer, directement et de façon fiable, l'évolution du faisceau dans d'autres plans. Brève description des figures De toute façon, les caractéristiques et avantages de l'invention apparaîtront mieux à la lumière de la description qui va suivre. Cette description porte sur les exemples de réalisation, donnés à titre explicatif et non limitatif, en se référant à des dessins annexés sur lesquels - La figure 1 représente un dispositif pour la mise en oeuvre de l'invention. - La figure 2 représente des étapes d'un traitement préalable de matrices d'intensité. - La figure 3 représente les étapes d'un premier mode de réalisation d'un procédé selon l'invention. - La figure 4 représente les étapes d'un deuxième mode de réalisation d'un procédé selon l'invention. - La figure 5 représente schématiquement un traitement dit "de propagation". Description détaillée de modes de réalisation de l'invention L'invention va d'abord être décrite à l'aide du mode de réalisation schématisé sur la figure 3. Un dispositif de mise en oeuvre est représenté schématiquement sur la figure 1. L'invention permet de calculer le champ complexe d'un faisceau laser 4 dans la pupille de départ définie par convention par le plan P'1 représenté en figure 1 à partir de mesures d'intensité effectuées dans N plans Pl, P2, P3 (N > 3) (Pi et P'i étant optiquement conjugués). On décrit ici le cas où N=3. Soit, Ai .icP1 l'amplitude complexe dans le plan i, l On mesure Ii=lAil2 intensité dans le plan i, à l'aide de caméras CCD 1, 2, 3. Les caméras sont reliées à un calculateur 6. Le procédé et la méthode de calcul sont implémentés par programmation de ce calculateur. Les instructions appropriées de programmation sont mémorisées sur des disques magnétiques ou dans des mémoires conventionnelles RAM ou ROM. Pour chacun des plans de mesure Pl, la mesure Il est fournie à une mémoire de calculateur 6. Sont également mémorisées - la valeur de la longueur d'onde du faisceau analysé, notée X, - les échelles des mesures exprimées en m/pixel (distance entre deux points successifs de chaque matrice, après interpolation) notées echx et echy, - les distances relatives entre les N plans de mesure d (Pj, P1) =d1,1. A partir des mesures d'intensité Ii, la méthode selon l'invention permet de retrouver la phase pl dans le plan numéro 1. On pourra alors calculer la phase ç du faisceau dans n'importe quel plan. La mesure de Ii fournie par chaque caméra se présente, après numérisation, sous la forme d'une matrice Ii(n',p'), comportant n' lignes et p' colonnes de nombres réels. Chaque matrice est rééchantillonnée en nxp où n et p sont des puissances de 2, proches de n' et p', pour pouvoir utiliser des transformées de Fourier rapides. Un traitement est préalablement effectué sur les matrices, pour s'affranchir de certains problèmes liés à la mesure : ce traitement est représenté schématiquement sur la figure 2. Une étape d'interpolation (étape 2-2) permet de passer d'un profil à n'xp' points au profil ayant nxp points, avec n et p puissances de 2 les plus proches de n' et p'. On définit donc n et p puis on interpole les données pour obtenir la matrice I', correspondant à I (étape 2-3). On centre le profil dans la fenêtre de manière à avoir le centre de gravité de la répartition d'intensité ainsi considérée au centre de la matrice (étape 2-4). On va utiliser, pour les calculs, l'amplitude dans les plans de mesure. On prend donc la racine, terme à terme, de la matrice en intensité (étape 2-5) pour obtenir enfin A' (étape 2-6). On effectue ensuite une normalisation des profils pour s'assurer de la conservation de l'énergie. On divise pour cela la matrice A' par la racine de la somme des carrés de ses termes (étape 2-7) pour obtenir l'amplitude recherchée A (étape 2-8). La figure 3 donne ensuite les étapes principales du calcul selon l'invention, dans le cas de trois plans de mesure. On réalise l'acquisition des trois intensités li, I2, I3 (étape 3-1) . On traite ensuite ces trois images numérisées comme décrit ci-dessus en liaison avec la figure 2 (étape 3-2). On obtient les amplitudes A1, A2 et A3 de dimensions nxp (étape 3-3). On choisit (étape 3-4) une matrice de phase de départ ei#1 de dimensions nxp (avec par exemple (Pi=0), que l'on multiplie avec Al (étape 3-5) pour obtenir la matrice B1 complexe nxp B1=A1 e(Pl On donne plus loin plus de détails sur cette multiplication. On applique ensuite à cette matrice B1 un traitement dit "de propagation" (étape 3-6), qui sera décrit plus loin, et qui permet d'obtenir une matrice B'2 à partir de la matrice B1, la distance de propagation étant la distance d12 entre les plans i=1 et i=2. On récupère l'information #2 de phase contenue dans B'2 (comme décrit ci-dessous) . Cette information est multipliée par l'amplitude A2 (étape 3-7). On obtient alors la matrice B2 à laquelle on applique la propagation sur une distance d23 (distance entre les plans i=2 et i=3) . On obtient la matrice B'3 (étape 3-8). De même, on extrait p3 de B'3 et on calcule B3, produit de #3 et de A3 (étape 3-5). En appliquant la fonction de propagation à B3 sur une distance -d13 (pour revenir au plan 1), on trouve le profil reconstruit B'1 (étape 3-10), contenant un nouveau terme cpi, qui sert a la boucle suivante (avec k=k+l) On peut aussi propager B3 vers un quatrième plan si on a choisi un nombre de plans supérieur à 3. Le procédé s'arrête, par exemple lorsque k atteint une valeur limite kmaXt ou lorsque le résultat est jugé satisfaisant. La valeur de B1 calculée correspond alors au profil d'amplitude complexe dans le premier plan (i=1). Une version améliorée de ce procédé, dite version "pondérée", est représentée schématiquement sur la figure 4. On réalise l'acquisition des trois intensités I1, 12, I3 (étape 4-1) . On traite ensuite ces trois images numérisées comme décrit ci-dessus (étape 4-2). On obtient donc les matrices réelles d'amplitude Al, A2, A3, de dimensions nxp (étape 4-3). On choisit une matrice de phase de départ et(P1 de dimensions nxp (étape 4-4, avec par exemple #1=0). On multiplie J(Pl avec A1 pour obtenir la matrice complexe B1 (étape 4-5) B1=A1 j(P1 Là encore, plus de détails sur cette multiplication seront donnés plus loin. On applique à cette matrice la fonction de propagation, associée aux distances d12 et d13 pour obtenir les matrices B'2 et B'3 (étapes 4-6 et 4-7). Ces dernières contiennent les informations de phase 92 dans le plan i=2, et V3 dans le plan i=3. On multiplie #2 et A2 pour obtenir B2 (étape 4-8), 3 et A3 pour obtenir B3 (étape 4-9) . On applique ensuite à chacune des matrices B2 et B3 le traitement de propagation, respectivement avec une distance -d12 et -d13 (étapes 4-10 et 4-11). Des deux matrices B'11 et B'12 ainsi obtenues, on extrait les deux phases cela et (P12. On obtient (étape 412) un nouveau terme pl en effectuant une moyenne pondérée de ces deux phases (le type de pondération est défini selon les besoins). Ce terme sert à la boucle suivante avec k=k+l. On peut également travailler avec un plus grand nombre N de plans (N > 3), la moyenne pondérée au début de la boucle k+l se faisant alors sur les N-l phases déterminées lors de la boucle k. On peut faire aussi une moyenne sur N-p phases, avec p > 1, mais cela revient à ignorer p-l plans, ou à pondérer p-l plan par 0 tout se passe alors comme si les mesures étaient faites sur N- (p-i) plans. On s'arrête lorsque k atteint une valeur limite kma. ou lorsque le résultat est jugé satisfaisant. La phase pl ainsi trouvée, combinée avec l'amplitude A., donne l'amplitude complexe dans le premier plan. Le traitement de propagation, correspondant aux étapes 3-6, 3-8, 3-10 et 4-6, 4-7, 4-10 et 4-ll des figures 3 et 4, est illustré sur la figure 5. A partir d'une matrice B(x,y) de départ, de dimensions nxp, on calcule une "matrice propagée" Bp (xp,yp), également de dimensions nxp. La matrice complexe B décrit l'amplitude et la phase du faisceau dans un certain plan. La matrice complexe Bp décrit l'amplitude et la phase du faisceau dans un autre plan. On peut passer de B à Bp, c'est-àdire d'un plan à un autre, à l'aide des équations de propagation du faisceau. On peut par exemple passer de B(x,y) à Bp(xp,yp) par convolution de B par la matrice de propagation dans le vide de l'optique de Fourrier. Une première étape 5-1, dite de "zero padding", consiste à centrer la matrice B de dimension nxp dans une matrice 2nx2p, composée uniquement de zéros, afin de faire une convolution linéaire. On obtient donc une matrice B(x,y) de taille 2nx2p. Par une transformation de Fourier rapide (étape 5-2), on obtient une matrice B'( > ,v), de dimensions 2nx2p. Une multiplication par un facteur (étape 5-3) permet obtenir une matrice B'p(up,vp) de dimensions 2nx2p. Par transformée de Fourier inverse (étape 5-4), on obtient une matrice propagée Bp(xp,yp) de dimensions 2nx2p. L'extraction (étape 5-5) est l'étape inverse qui consiste à retirer les termes (le plus souvent nuls) périphériques de la matrice Bp pour revenir à la dimension de départ. La FFT est l'algorithme de transformée de Fourier rapide qui existe dans la plupart des bibliothèques de fonctions. Tout autre algorithme de propagation pourrait être utilisé avec succès. En ce qui concerne le traitement des phases et les multiplications, on notera d'abord que l'ensemble des calculs de propagation s'effectue avec des matrices d'amplitudes complexes regroupant amplitude et phase. C'est l'information de phase qui est conservée à chaque étape alors qu ' on reprend les profils d'amplitude acquis. Lorsqu'on manipule une matrice B d'amplitude complexe (B=IBI.ei(P), on cherche à extraire la phase p. Il existe deux possibilités de gestion de la phase - selon une première possibilité, le calcul de la phase peut se faire en calculant l'argument de l'amplitude complexe. Lors de la multiplication de la phase par un profil enregistré Ai, on multiplie terme à terme la matrice Ai avec la matrice des e1. On obtient alors la nouvelle amplitude complexe B ; - selon une deuxième possibilité pour alléger la charge de calcul, on peut ne calculer la phase # d'une matrice B que lorsque l'on veut son expression en fin de reconstruction. Pour les calculs intermédiaires, on travaille directement avec la matrice ei(P L'obtention de ce terme se fait directement en divisant terme à terme la matrice complexe B par sa matrice module Si. . Lorsque l'on veut ensuite multiplier l'amplitude enregistrée par la phase on multiplie terme à terme la matrice A par la matrice ei# ainsi calculée. La méthode selon l'invention, grâce à un nombre de plans de référence supérieur ou égal à 3, permet en quelque sorte de "forcer" le profil calculé en chacun des plans de mesure. Il en résulte un profil intensité/phase de départ beaucoup plus stable et plus proche de la réalité mesurée. Un excellent accord a pu être constaté entre les profils mesurés et les profils calculés à partir du plan origine. Ceci prouve la véracité du profil de phase reconstitué. Cet accord est nettement meilleur que ce que l'on peut espérer obtenir avec un analyseur de type Hartmann-Shack ou Zygo, et ceci pour un coût moindre. REVENDICATIONS 1. Procédé itératif de détermination de la répartition d'intensité et de phase d'un faisceau lumineux cohérent (4), dans un plan, comportant - la mesure de l'intensité Ii, i=i. . . N, du faisceau, dans N plans, N > 3, dont le plan dans lequel on veut déterminer la répartition d'intensité et de phase, - le choix, pour le plan i=l, d'une matrice de phase (Pi de départ et le calcul d'une matrice d'amplitude complexe, par produit terme à terme de la matrice de phase ex91 avec la matrice d'amplitude A1 correspondante, - pour chaque plan j > 1 * détermination d'une matrice complexe propagée B'j à partir de la matrice mesurée d'amplitude Ajl du plan j-l et de la matrice de phase ç du plan j-l, # extraction, de B'j, d'une matrice de phase du plan , - itération du procédé jusqu'à convergence, avec j=l lorsque j-l=N. 2. Procédé selon la revendication 1, chaque matrice Ai d'intensité comportant n lignes et p colonnes, n et p étant des puissances de 2. 3. Procédé selon la revendication 1 ou 2, la phase prise en compte dans le plan j=1 à chaque itération étant une moyenne pondérée entre les N-l profils de phase déterminés antérieurement. 4. Procédé selon l'une des revendications 1 à 3, la matrice complexe propagée B'j étant obtenue - par calcul du produit de la matrice mesurée d'amplitude Ajl du plan j-l et de la matrice de phase eiç du plan j-l, - par convolution de ce produit par la matrice de propagation dans le vide de l'optique de Fourrier. 5. Procédé selon l'une des revendications 1 à 4, la matrice de phase ç du plan j-1 étant obtenue par calcul de l'argument de la matrice complexe propagée B'j-1 associée au plan j-i. 6. Procédé selon l'une des revendications 1 à 4, la matrice de phase ç du plan j-1 étant obtenue par division, terme à terme, de la matrice complexe propagée B'j-1 par sa matrice module IB'i-îI 7. Procédé selon l'une des revendications 1 à 6, le faisceau étant préalablement mis en forme par des moyens qui en limitent son étendue, et le procédé comportant en outre une étape d'étalonnage. 8. Dispositif pour déterminer l'intensité et la phase d'un faisceau lumineux (4) dans un plan comportant - des moyens (1, 2, 3) pour mesurer l'intensité A (i=l,...N) du faisceau dans N plans, N > 3, - des moyens de calcul (6) programmés pour, à partir d'une matrice de phase (Pi de départ, correspondant à un plan i=1 * calculer une matrice d'amplitude complexe (B1) par produit de la matrice de phase ei#1 avec la matrice d'intensité (A1) correspondante, * pour chaque plan j > 1 calculer une matrice complexe propagée B'j à partir de la matrice mesurée d'amplitude A du plan j-1 et de la matrice de phase p du plan j-i, extraire une matrice de phase du plan j de la matrice B'j. * itérer le calcul ci-dessus, en faisant j=l lorsque j-1=N. 9. Dispositif selon la revendication 8, les moyens (6) de calcul étant en outre programmés pour transformer la matrice d'intensité mesurée en matrice à n lignes et p colonnes, n et p étant des puissances de 2. 10. Dispositif selon l'une des revendications 8 ou 9, les moyens (6) de calcul étant en outre programmés pour produire à la fin de chaque itération, une moyenne pondérée des N-l profils de phases déterminés au cours de cette itération. 11. Dispositif selon l'une des revendications 8 à 10, les moyens (6) de calcul étant programmés pour calculer une matrice complexe propagée B'j - par produit de la matrice mesurée d'intensité Ajl du plan j-l et de la matrice de phase eiç du plan j-l, - puis par convolution de ce produit par la matrice de propagation dans le vide de l'optique de Fourrier. 12. Dispositif selon l'une des revendications 8 à 11, les moyens (6) de calcul étant programmés pour fournir la matrice de phase ç du plan j-l par calcul de l'argument de la matrice complexe B'j-1 associée au plan j-i. 13. Dispositif selon l'une des revendications 8 à 11, les moyens (6) de calcul étant programmés pour fournir la matrice de phase ç du plan j-l par division, terme à terme, de la matrice complexe propagée B'j par sa matrice module lB'-i 14. Dispositif selon l'une des revendications 8 à 12, comportant en outre des moyens pour mettre en forme le faisceau et en limiter l'étendue afin de l'adapter aux caméras, les moyens (6) de calcul étant en outre programmés pour que le calcul soit affranchi de la présence desdits moyens de mise en forme du faisceau.