Brevet FR2957146A1 — publié le 20110909

DOMAINE DE L'INVENTION L'invention concerne le domaine de l'interférométrie à transformée de Fourier statique. Plus précisément, elle concerne un procédé de métrologie pour interféromètre à transformée de Fourier statique. ETAT DE L'ART En interférométrie à transformée de Fourier statique, on élabore un interférogramme en créant un nombre fini de différences de chemin optique réparties spatialement, et acquises par une matrice de détection. Ceci peut être réalisé à l'aide d'une multitude de surfaces réfléchissantes en parallèle, appelées échelettes, et d'un système optique imageant ces échelettes sur un détecteur. Ce concept d'interféromètre statique ou à échelettes est largement décrit dans le brevet FR2787186. Les instruments de type interféromètre statique connaissent divers applications, comme par exemple la surveillance de la pollution atmosphérique par spectrométrie infrarouge, qui est réalisée en embarquant l'instrument sur un satellite en orbite. De manière générale, l'utilisation de l'interféromètre statique implique des spécifications très sévères sur l'instrument. Un point critique des spécifications est la connaissance très précise des différences de chemin optique causées par les miroirs comportant une multitude de surfaces réfléchissantes en parallèle, ou échelettes. La précision requise est en général de l'ordre du nanomètre. De plus, un suivi régulier des différences de chemin optique est 25 nécessaire car leurs valeurs fluctuent en raison d'effets thermiques ou mécaniques. Il existe divers procédés de métrologie permettant d'obtenir une mesure des différences de chemin optique causées par les échelettes. Un de ces procédés consiste à illuminer chaque échelette avec un laser 30 accordable, et à balayer une bande de longueurs d'onde donnée. Ceci permet de mesurer les différences de chemin optique causées par la géométrie des échelettes. Un inconvénient de ce procédé est qu'il n'est pas applicable à toutes les longueurs d'onde, car il est difficile de trouver un laser accordable à certaines longueurs d'ondes, comme par exemple à 10pm. Un autre inconvénient encore est que ce procédé est long, et 5 nécessite de très nombreuses images laser. Par conséquent, il rend l'instrument de mesure indisponible pendant une longue période, ce qui est préjudiciable pour de nombreuses applications, comme par exemple dans le cas de mesures effectuées à partir d'un satellite en orbite. 10 Enfin ce type de procédé nécessite un dispositif complexe et lourd à mettre en oeuvre. PRESENTATION DE L'INVENTION L'invention propose de pallier au moins un de ces inconvénients. 15 A cet effet, l'invention propose un procédé de métrologie d'un interféromètre statique comprenant une lame séparatrice apte à séparer un faisceau incident vers un premier et un deuxième miroirs, le premier miroir présentant N surfaces réfléchissantes sensiblement planes et parallèles entre elles, lesdites surfaces étant décalées entre elles et réalisant N 20 premières hauteurs de marche, le deuxième miroir présentant M surfaces réfléchissantes sensiblement planes et parallèles entre elles, lesdites surfaces étant décalées entre elles et réalisant M deuxièmes hauteurs de marche, les hauteurs de marche correspondant à la distance entre chaque surface et la lame, de sorte que les deux miroirs fournissent une pluralité de 25 faisceaux réfléchis présentant entre eux des différences de chemin optique égales aux N*M sommes croisées des premières hauteurs de marche et des deuxièmes hauteurs de marche, ledit procédé étant caractérisé en ce qu'il comprend au moins les étapes consistant à : - acquérir un interférogramme d'intensité (I(x,y)) à partir d'un faisceau 30 incident monochromatique, - estimer lesdites N+M premières et deuxièmes hauteurs de marche par résolution d'un système d'équations liant l'intensité I(x,y) mesurée aux N+M premières et deuxièmes hauteurs de marche, ledit système exploitant les corrélations et dépendances géométriques des miroirs, et - déduire les N*M différences de chemin optique à partir des N+M premières et deuxièmes hauteurs de marche. L'invention est avantageusement complétée par les caractéristiques suivantes, prises seules ou en une quelconque de leur combinaison techniquement possible : - le système d'équations liant l'intensité I(x,y) aux premières et deuxièmes hauteurs de marche hMm, hM2 s'écrit I(x,y)=cos(2rr.6.(hMl+hM2)), où a est le nombre d'onde du faisceau incident monochromatique ; - les premières et deuxièmes hauteurs de marche sont chacune prises comme égales à la somme d'une composante connue (pMl ;pM2) et d'une composante inconnue (ami ;aM2) et en ce que l'on estime les N+M composantes inconnues par résolution dudit système liant l'intensité aux hauteurs de marche; - l'étape d'estimation des N+M composantes inconnues est itérée en intégrant à la composante connue utilisée pour l'itération courante la composante inconnue déterminée à l'itération précédente ; - la composante inconnue est la somme d'une composante axiale selon la direction du faisceau incident qui tombe sur ce miroir, et d'une composante transversale supérieure ou égale à zéro représentant les inclinaisons par rapport à au moins une des directions de l'espace des surfaces réfléchissantes du premier miroir et/ou des surfaces réfléchissantes du deuxième miroir, l'étape d'estimation des hauteurs de marche comprenant l'estimation desdites composantes axiales et transversales ; - la résolution du système liant l'intensité aux hauteurs de marche comprend une étape d'estimation des sommes croisées (aM,+aM2) des composantes inconnues à partir de l'équation I(x,y) = cos(2rr.6(pM1+pM2+aM1+aM2)), pour en déduire une estimation desdites N+M premières et deuxièmes hauteurs de marche ; - l'estimation des sommes (aM1+aM2) des composantes inconnues est mise en oeuvre en appliquant la fonction cosinus inverse à l'équation I(x,y)=cos(2Tr.a.(pM1+pM2+aM1+aM2)) ; - on estime les composantes inconnues aM1,aM2 représentées par un vecteur A, à partir de l'estimation des sommes (aM1+aM2) desdites composantes inconnues, représentées par un vecteur Y, en inversant l'équation matricielle Y = W * A, où W est la matrice de passage entre Y et A; - l'estimation des composantes inconnues ami, aM2 est mis en oeuvre en effectuant un développement limité de l'équation I(x,y)=cos(2Tr.6.(pM1+pM2+aM1+aM2)) autour de l'infiniment petit 2Tr.6.(aM1+aM2) ; - le procédé comprend une étape d'estimation de l'amplitude du faisceau incident monochromatique ; - l'étape d'estimation de l'amplitude a du faisceau incident monochromatique de nombre d'onde a est mise en oeuvre en minimisant selon a la fonction f(a) = I(x,y) û a.cos(2Tr.6.6(x,y)). L'invention présente de nombreux avantages. Un avantage de l'invention est de permettre une mesure des 20 différences de chemin optique en utilisant simplement une source monochromatique de type laser. Un autre avantage est de permettre une mesure rapide des différences de chemin optique, à partir d'un unique interférogramme. Un autre avantage encore est que l'invention permet d'obtenir une 25 mesure des différences de chemin optique avec une précision élevée. Un autre avantage encore est que l'invention nécessite des équipements moins encombrants que ceux utilisés dans l'art antérieur. PRESENTATION DES FIGURES 30 D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, et qui doit être lue en regard des dessins annexés sur lesquels : - la figure 1 est un schéma optique général d'un mode de réalisation d'un interféromètre statique ; - la figure 2 est un exemple d'image spatiale d'un interférogramme obtenue avec un interféromètre statique ; - la figure 3 est un schéma bloc d'une réalisation possible du procédé selon l'invention. DESCRIPTION DETAILLEE Présentation de l'interféromètre En figure 1, on a représenté un schéma optique de principe d'un interféromètre 1 statique. L'interféromètre 1 statique comprend une lame séparatrice 3, apte à séparer un faisceau 2 incident en plusieurs faisceaux dirigés sur deux miroirs constitués d'un premier miroir M1 et d'un deuxième miroir M2. Le premier miroir M1 présente un ensemble de N surfaces (SI, S2, S3... ,SN) réfléchissantes sensiblement planes et parallèles entre elles. Les surfaces (SI, S2, S3... ,SN) réfléchissantes sont décalées l'une par rapport à l'autre et réalisent N premières hauteurs de marche hMl(x). De même, le miroir M2 présente un ensemble de surfaces (S'1, S'2, S'3... ,S'M) réfléchissantes sensiblement planes et parallèles entre elles. Les surfaces (S'1, S'2, S'3... ,S'M) réfléchissantes sont décalées l'une par rapport à l'autre et réalisent M deuxièmes hauteurs de marche hM2(y). Les premières et deuxièmes hauteurs de marche hMl(x), hM2(y) correspondent à la distance entre chaque surface (SI, S2, S3... ,SN, S'1, S'2, S'3... ,S'M) et la lame 3 séparatrice, selon la direction (x,y) du faisceau issu de la lame 3 séparatrice et arrivant sur la surface réfléchissante. Chaque hauteur de marche est avantageusement définie comme la distance entre le centre de la surface réfléchissante et la lame 3 séparatrice selon la direction du faisceau issu de la lame 3 et dirigé au centre de ladite surface réfléchissante. Comme cela sera explicité par la suite, ces distances sont exprimées algébriquement (signe positif ou négatif selon la convention d'orientation des axes (x,y)). Pour le miroir M1, le décalage se fait par exemple selon la direction x du faisceau 2 arrivant sur ce miroir et issu de la lame 3 séparatrice. Pour le miroir M2, le décalage se fait par exemple selon la direction y du faisceau 2 arrivant sur ce miroir et issu de la lame 3 séparatrice. Ceci n'est pas obligatoire et la direction de décalage des surfaces peut être selon toute direction de l'espace, et dépend notamment de l'orientation angulaire de la lame séparatrice 3. En figure 1, la lame séparatrice 3 est orientée à 45°. Si les hauteurs de marche sont ordonnées, le décalage des surfaces réfléchissantes donne au miroir un aspect semblable à des marches d'escalier. Il est possible d'avoir des hauteurs de marche non ordonnées donnant au miroir l'aspect de marches d'escaliers ascendantes puis descendantes, ou le contraire. Il est clair qu'il est possible d'utiliser deux miroirs M1 et M2 ayant le même nombre de surfaces réfléchissantes, ce qui implique que N=M. Un nombre typique de surfaces réfléchissantes pour un miroir selon l'invention est « N=32 », donné à titre d'exemple non limitatif. Les surfaces (SI, S2, S3... ,SN et/ou S'1, S'2, S'3... ,S'M) réfléchissantes sont sensiblement planes et parallèles, puisqu'il est difficile d'obtenir des surfaces (SI, S2, S3... ,SN; S'1, S'2, S'3... ,S'M) réfléchissantes parfaitement planes et parallèles. En effet, la fabrication des miroirs (Ml, M2) impose un usinage très précis, et l'homme du métier est conscient qu'un tel usinage implique des erreurs de réalisation. De manière avantageuse, l'un des deux miroirs (par exemple Ml) possède des surfaces (SI, S2, S3... ,SN) réfléchissantes présentant des décalages très fins (de l'ordre de quelques dizaines à une centaine de microns) selon l'axe x du faisceau 2 incident à ce miroir, tandis que l'autre miroir M2 possède des surfaces (S'1, S'2, S'3... ,S'M) réfléchissantes présentant des décalages plus grands selon l'axe y du faisceau 2 incident à ce miroir. Ceci permet d'obtenir des différences de chemin optique plus finement échantillonnées. Par exemple, l'un des miroirs (par exemple M1) a un décalage selon l'axe x entre la surface SI réfléchissante, qui est la surface la moins décalée selon l'axe x, et la surface SN réfléchissante, qui est la surface la plus décalée selon l'axe x, égal au décalage selon l'axe y entre deux surfaces réfléchissantes consécutives (Si et S'i+1) de l'autre miroir M2. On voit en figure 1 que ces surfaces (SI, S2, S3... ,SN; S'1, S'2, S'3... ,S'M) réfléchissantes ont par exemple une forme rectangulaire. Par ailleurs, les premiers et deuxièmes miroirs M1 et M2 sont conjugués par rapport à la lame séparatrice 3, ce qui signifie qu'ils sont images l'un de l'autre par rapport à la dite lame. Il est à noter que pour obtenir cette conjugaison, il n'est pas obligatoire que l'axe médian 11 du miroir M1 et l'axe médian 12 du miroir M2 soient perpendiculaires, comme cela est représenté en figure 1. L'interféromètre 1 statique comporte en outre une matrice 10 de détection qui détecte l'intensité des faisceaux lumineux renvoyés par les miroirs M1 et M2. Enfin, l'interféromètre 1 statique comprend une optique 13 d'entrée par laquelle entrent des faisceaux lumineux émis par une source 4 lumineuse, et une optique 14 de sortie intercalée entre la lame séparatrice 3 et la matrice 10 de détection. L'optique 13 d'entrée par laquelle entrent les faisceaux lumineux peut tout aussi bien se situer entre les miroirs et le détecteur, ce qui fait pénétrer les faisceaux lumineux à travers l'autre face de la lame séparatrice 3. Mise en oeuvre de l'interférométrie Dans le procédé selon l'invention, un faisceau 2 incident monochromatique est émis en direction des miroirs M1, M2. Ce type de faisceau est en général émis par une source 4 monochromatique. La source 4 monochromatique est typiquement un laser. Ce type de source est aisé à utiliser et est peu encombrant. En général, l'interféromètre étant utilisé pour analyser un rayonnement compris dans une bande spectrale donnée (par exemple, la bande spectrale du rayonnement infrarouge), l'on choisit alors dans le procédé de métrologie selon l'invention un faisceau 2 incident de longueur d'onde comprise dans ladite bande spectrale. Le faisceau 2 incident émis par la source 4 monochromatique traverse l'optique 13 d'entrée, pour être dirigé vers la lame séparatrice 3. La lame séparatrice 3 sépare le faisceau 2 incident par la source 4 monochromatique de telle sorte que chaque miroir (M1,M2) reçoit alors un faisceau 2 incident sur une surface de chacun desdits miroirs (M1,M2). Etant donné que les surfaces (SI, S2, S3... ,SN; S'1, S'2, S'3... ,S'M) sont décalées entre elles, par exemple selon l'axe x et/ou y des faisceaux incidents à ces surfaces, les faisceaux réfléchis par les surfaces suivent un chemin de longueur différente et présentent alors entre eux un nombre fini N*M de différences de chemin optique 6(x,y), égales aux N*M sommes croisées des N premières hauteurs de marche (hMl(x)) des surfaces (SI, S2, S3... ,SN) réfléchissantes du miroir M1 et des M deuxièmes hauteurs de marche (hM2(y)) des surfaces (S'1, S'2, S'3... ,S'M) réfléchissantes du miroir M2. La figure 2 illustre le trajet d'un faisceau séparé par la lame 3 séparatrice et rencontrant les surfaces S1 et S'1 dont les hauteurs de marche sont différentes, ce qui induit une différence de chemin optique. Comme évoqué plus haut, les hauteurs de marche sont définies de manière algébrique, selon une convention d'orientation, ce qui explique que les différences de chemin optique 6(x,y), sont égales aux N*M sommes croisées des N premières hauteurs de marche (hMl(x)) et des M deuxièmes hauteurs de marche (hM2(y)). Dans l'exemple de la figure 2, les premières hauteurs de marche (hMl(x)) sont définies comme positives, au vu de l'orientation de l'axe x, tandis que les deuxièmes hauteurs de marche (hM2(y)) sont définies comme négatives, au vu de l'orientation de l'axe y. Il est clair qu'il s'agit simplement d'une convention de signe et qu'il est possible de procéder à une définition inverse, en inversant les sens des axes x et y. Par ailleurs, les différences de chemin optique 6(x,y) peuvent également être dues au fait que les surfaces réfléchissantes (Si, S2, S3,... et/ou Si', S2', S3',...) ne sont pas parfaitement planes et parallèles entre elles. Il existe d'autres causes provoquant une variation des hauteurs de marche et donc des différences de chemin optique 6(x,y), comme par exemple des effets thermiques, des effets de contrainte mécanique liée au support, des effets liés à la gravité (notamment pour les applications spatiales), l'ensemble desdites causes étant bien connu de l'homme du métier. Comme les faisceaux réfléchis par les premiers et deuxièmes miroirs M1 et M2 présentent entre eux des différences de chemin optique 6(x,y), la matrice 10 de détection acquiert un interférogramme, dont l'intensité I(x,y) dépend des différences de chemin optique 6(x,y) des faisceaux réfléchis par les miroirs M1 et M2. L'intensité I(x,y) dépend également du nombre d'onde a du faisceau 2 incident monochromatique. On rappelle que le nombre d'onde a est défini comme l'inverse de la longueur d'onde À, c'est- 1 à-dire6 = ù. X On a représenté en figure 2 un exemple d'interférogramme d'intensité I(x,y) obtenu avec un faisceau 2 incident monochromatique. En pratique, l'interférogramme d'intensité I(x,y) est obtenu en illuminant simultanément avec un faisceau 2 incident la surface S'l réfléchissante et la surface SI réfléchissante, la surface S'i réfléchissante et la surface S2 réfléchissante et l'ensemble des combinaisons surfaces Si et S'i réfléchissantes, pour i allant de 1 à N et j allant de 1 à M. On obtient donc un interférogramme dont l'intensité I(x,y) est enregistrée sur une matrice de N*M facettes, en général de formes carrées. Chaque facette carrée de position (x,y) sur la matrice de taille N*M correspond à une différence de chemin optique 6(x,y), ce qui donne en tout N*M mesures d'intensité. Estimation des différences de chemin optique Etant donné que l'interféromètre 1 statique possède une configuration de miroirs particulière, il existe des corrélations entre les différentes intensités obtenues dans l'interférogramme d'intensité I(x,y). Ainsi, toutes les différences de chemin optique 6(x,y) liées à une même surface du premier miroir M1 seront corrélées. De même, toutes les différences de chemin optique 6(x,y) liées à une même surface du deuxième miroir M2 seront corrélées. Dans l'interférogramme d'intensité I(x,y) de la figure 2, cela signifie que dans une ligne x donnée, les valeurs des différences de chemin optique 6(x,y) sont corrélées. Il est en de même pour une colonne y donnée. Le procédé selon l'invention tire avantageusement parti du fait que l'interférogramme d'intensité I(x,y) est obtenu en croisant les hauteurs de marche de deux miroirs, ce qui entraîne une corrélation entre les intensités obtenues. Pour mesurer les différences de chemin optique en exploitant ces propriétés de corrélation, le procédé selon l'invention comprend notamment les étapes illustrées en figure 3. La première étape, déjà décrite, consiste en l'acquisition 20 d'un interférogramme d'intensité I(x,y), via une illumination des surfaces réfléchissantes par un faisceau 2 incident monochromatique. A ce stade, l'on dispose de N*M mesures I(x,y), acquises par la matrice 10 de détection, mais également de N*M inconnues, correspondant 25 aux N*M différences de chemin optique 6(x,y) à estimer. Estimer les N*M inconnues 6(x,y) directement et uniquement à partir des N*M mesures I(x,y) donnerait une estimation peu précise, étant donné que la relation liant les mesures aux inconnues n'est pas linéaire et que les mesures sont nécessairement bruitées. 30 En effet, la relation existant entre l'intensité I(x,y) mesurée et les différences de chemin optique 6(x,y) s'écrit en général selon le système d'équations I(x,y) = cos(2rr.6.(6(x,y))), où a est le nombre d'onde du faisceau 2 incident monochromatique. Le procédé selon l'invention exploite la corrélation qui existe entre les différences de chemin optique b(x,y) mesurées. Cette corrélation se traduit mathématiquement sur les lignes et colonnes de l'interférogramme. Par exemple, si les surfaces réfléchissantes de chaque miroir sont parfaitement parallèles, pour une ligne x1 donnée de l'interférogramme, les différences de chemin optique vont s'écrire ; b(x1,Y1) = hM1(xl) + hM2(Y1) b(xl,Y2) = h M1(xl) + h M2(Y2) b(xl,Y3) = h M1(xl) + h M2(Y3) (...) b(xl,YM) = hM1(xl) + hM2(YM) Comme on peut le constater, le terme hM1(xl), correspondant à la hauteur de marche de la surface S1 du premier miroir M1, est présent pour toutes les différences de chemin optique de la ligne x1. Ces propriétés sont vérifiées pour toutes les lignes et colonnes (x,y) de l'interférogramme. Ainsi, plutôt que de chercher à estimer de manière abrupte les différences de chemin optique, le procédé selon l'invention se base sur l'estimation des N+M premières et deuxièmes hauteurs de marche hM1,hM2 des surfaces réfléchissantes des premier et deuxième miroirs M1 et M2. Ainsi, le procédé selon l'invention met en oeuvre une étape 22 d'estimation desdites N+M premières et deuxièmes hauteurs de marche (hM1(x);hM2(y)) de chaque surface réfléchissante (Si, S2, S3... ,SN; S'1, S'2, S'3... ,S'M) par résolution d'un système d'équations liant l'intensité I(x,y) mesurée aux N+M premières et deuxièmes hauteurs de marche hM1(x) et hM2(y) des surfaces réfléchissante (S1, S2, S3... ,SN; S'1, S'2, S'3... ,S'M). Ce système d'équations exploite donc les corrélations et dépendances géométriques des miroirs de l'instrument. Le système d'équations liant l'intensité I(x,y) aux premières et deuxièmes hauteurs de marche est préalablement modélisé lors d'une étape 21, et pourra typiquement consister en une modélisation s'écrivant I(x,y)=cos(2rr.a.(hM1+hM2)). Le nombre d'inconnues décrivant les différences de chemin optique est ainsi grandement réduit. La résolution peut faire intervenir tout moyen de linéarisation, d'inversion, ou de résolution numérique connu. Cette résolution peut notamment faire intervenir dans une étape intermédiaire une estimation des sommes croisées des premières et deuxièmes hauteurs de marche, pour en déduire par la suite les hauteurs de marche elles-mêmes. Il est possible de ne pas utiliser un partie des mesures I(x,y) jugées comme trop bruitées ou indésirables, et donc d'utiliser uniquement un sous-ensemble des N*M mesures I(x,y) pour l'estimation des N+M inconnues. Enfin, la dernière étape 23 consiste à déduire les N*M différences de chemin optique 6(x,y) des N+M hauteurs de marche. Cette déduction est immédiate puisque les différences de chemin optique (6(x,y)) sont égales aux N*M sommes croisées des premières hauteurs de marche hMl(x) et des deuxièmes hauteurs de marche hM2(y), selon la relation 6(x,y) = hMl(x) + hM2(y). Il est à noter que le procédé selon l'invention présente de nombreux avantages. Comme on peut le voir en figure 1, un avantage de l'invention est de permettre une mesure des différences de chemin optique en utilisant une source monochromatique de type laser, qui est beaucoup plus simple et moins onéreuse qu'un laser accordable. De plus, une seule acquisition d'interférogramme est nécessaire, ce qui permet une mesure rapide des différences de chemin optique. Le temps nécessaire pour mettre en oeuvre le procédé de métrologie selon l'invention est donc bien inférieur à celui des procédés de l'art antérieur. Enfin, le procédé selon l'invention permet une mesure des différences de chemin optique avec une bonne précision. Par ailleurs, on modélise avantageusement dans le procédé selon l'invention les premières et deuxièmes hauteurs de marche (hMl(x) ;hM2(y)) comme égales à la somme d'une composante connue (pMl;pM2) et d'une composante inconnue (aMi;aM2). Cela revient à poser hMl(x) = pM1 + ami et hM2(y) = PM2 + aM2. Dans ce cas, l'étape 22 d'estimation des N+M premières et deuxièmes hauteurs de marche (hM1,hM2) nécessite uniquement l'estimation des N+M composantes inconnues (aM1;aM2). Il est très fréquent que l'on dispose d'une connaissance préalable 5 des hauteurs de marche, qui peut être une connaissance a priori ou provenant d'une précédente mesure. Cette connaissance peut provenir notamment de données du constructeur, donnant une estimation de la position des marches ce qui permet d'avoir une connaissance des hauteurs de marche avec une 10 certaine marge d'erreur. Cette connaissance peut alternativement provenir d'une mesure selon un procédé de l'art antérieur. Enfin, cette connaissance peut provenir d'une itération précédente du procédé de métrologie selon l'invention, comme cela sera explicité par la suite. La composante inconnue (ami ;aM2) correspond à l'erreur de 15 connaissance, qui peut correspondre par exemple à une évolution temporelle des gradients de température dans l'instrument (et donc des indices de réfraction et des différences de chemin optique), et que l'on cherche à connaître le plus précisément possible. Avantageusement, on itère l'étape 22 d'estimation des N+M 20 composantes inconnues (ami ;aM2) en intégrant à la composante connue (pM1 ;pM2) utilisée pour l'itération courante la composante inconnue (ami ;aM2) déterminée à l'itération précédente. Ceci signifie que la composante connue (pM1(n);pM2(n)) à l'itération numéro « n » est prise comme égale à la somme de la composante connue 25 (pM1(n-1);pM2(n-1)) à l'itération numéro « n-1 » et de la composante inconnue (aM1(n 1);aM2(n-1)) déterminée à l'itération numéro « n-1 ». Ainsi, l'on affine à chaque itération l'estimation des hauteurs de marche, en prenant en compte la mesure effectuée à l'itération précédente. La mesure est donc rendue plus précise. 30 Dans un mode de réalisation particulier, on prend également en compte l'éventuelle inclinaison des surfaces réfléchissantes dans le sens de leur longueur. Pour ce faire, il est possible de modéliser la composante inconnue (ami ;aM2) comme la somme d'une composante axiale (a'M, ;a'M2) selon la direction (x ; y) du faisceau 2 incident qui tombe sur ce miroir, et d'une composante transversale (iMl ;iM2) supérieure ou égale à zéro représentant les inclinaisons par rapport à au moins une des directions de l'espace (x,y) des surfaces (SI, S2, S3... ,SM) réfléchissantes du premier miroir (Ml) et/ou des surfaces (S'1, S'2, S'3... ,S'M) réfléchissantes du deuxième miroir (M2). En général, chaque composante transversale (iMl ;iM2) est modélisée selon : /ml (x, y) = bM (x).y et im2 (x, y) = bM2 (Y)•x , où (0M, , 0M2) sont des coefficients d'inclinaisons supérieurs ou égaux à zéro. Cette modélisation correspondant à une pente d'inclinaison constante selon la longueur de la surface réfléchissante. Dans ce cas, le nombre d'inconnues augmente, mais reste tout de même strictement inférieur aux N*M mesures. L'étape d'estimation 22 comprend alors l'estimation des composantes axiales (a'M, ;a'M2) et des coefficients d'inclinaison (b'Ml ;b'M2), ce qui permet d'en déduire les composantes inconnues aM,,aM2, et donc les hauteurs de marche hMl,hM2. On décrira à présent des exemples de mise en oeuvre de l'estimation des différences de chemin optique. Exemples de réalisation de l'estimation des différences de chemin optique Comme décrit précédemment, il est possible d'exprimer à l'étape 21 l'intensité I(x,y) en fonction des premières et deuxièmes hauteurs de marche selon le système d'équations : I(x,y)=cos(2rr.a.(hMl+hM2)) Cette modélisation de l'interférogramme d'intensité I(x,y) est obtenue 30 à partir d'un faisceau 2 incident monochromatique de nombre d'onde a. Il est à noter que d'autres modélisations sont possibles, notamment si l'on prend en compte une fonction d'auto-apodisation (diminution du contraste quand la différence de chemin optique augmente), ou des effets de gains inter-pixels sur la matrice 10 de détection. On suppose ici que les effets décrits ci-dessus ont été corrigés et que la ligne de base, c'est à dire la partie non modulée du signal de l'interférogramme, a été retirée. De plus, le contraste de l'interférogramme est supposé avoir été ramené à une valeur unitaire, par un étalonnage préalable des mesures. En se plaçant dans le cas où l'on dispose d'une connaissance a priori des différences de chemin optique, il est possible de réécrire ce système d'équations selon : I(x,y)=cos(2rr.a.(pM1+pM2+aM1+aM2)), où pM1, pM2 sont les composantes connues des hauteurs de marche, et ami, aM2 sont les composantes inconnues à estimer. On dispose donc à l'issue de l'étape 21 d'un système d'équations liant l'intensité I(x,y) mesurée aux composantes inconnues ami, aM2 des 15 hauteurs de marche. Dans un premier mode de réalisation, l'on met en oeuvre l'étape 22 en résolvant les équations I(x,y)=cos(2rr.6.(pM1+pM2+aM1+aM2)) par développement du cosinus et développement limité autour du terme 2rr.6.(aM1+aM2), qui est un infiniment petit. 20 Ainsi, l'équation I(x,y)=cos(2rr.6.(pM1+pM2+aM1+aM2)) devient: I(x,y)=cos(2TT.a.(pM1+pM2+aM1+aM2)) I(x,y)=cos(2rr.6.(pM1+pM2))*cos(2rr6(aM1+aM2)) - sin(2rr.6.(pM1+pM2)) sin(2rr6(aM1+aM2)) 25 I(x,y)=cos(2rr.6.(pM1+pM2))*[cos(2rr6(aM1))cos(2rr6(aM2)) û sin(2rr6(aM1))sin(2rr6(aM2))]- sin(2rr.6.(pM1+pM2))*2ua(aM1+aM2) I(x,y)=cos(2rr.6.(pM1+pM2))*[1 û 2rr6(aM1))2rr6(aM2)]-sin(2rr.6.(pM1+pM2))*2rr6(aM1+aM2) 30 On résout alors ce système d'équations en déterminant les N+M composantes inconnues ami et aM2 à partir des N*M mesures I(x,y), ou d'un sous-ensemble de ces N*M mesures. Ce sous-ensemble est choisi par exemple en éliminant les mesures trop bruitées. Une fois les composantes inconnues aMi et aM2 estimées, l'on peut déterminer immédiatement les premières et deuxièmes hauteurs de marche (via les relations hMl(x) = pM1 + aMi et hM2(y) = PM2 + aM2) et par la suite les différences de chemin optique à l'étape 23, via la relation 6(x,y) = hMl(x) + hM2(y). Comme énoncé précédemment, il est avantageux d'itérer l'étape 22 d'estimation en intégrant à la composante connue (pMi ;pM2) utilisée pour l'itération courante la composante inconnue (ami ;aM2) déterminée à l'itération précédente. Dans un deuxième mode de réalisation de l'invention, l'étape 22 d'estimation des premières et deuxièmes hauteurs de marche hMi,hM2, impliquant la résolution du système d'équations liant les premières et deuxièmes hauteurs de marche à l'intensité I(x,y), est mise en oeuvre en estimant dans une étape intermédiaire les sommes croisées aM,+aM2 des composantes inconnues, pour en déduire par la suite les dites composantes inconnues aMi, aM2, et par là les premières et deuxièmes hauteurs de marche hMi,hM2. Ainsi, les hauteurs de marche ne sont pas estimées directement mais déduites des sommes croisées de leurs composantes inconnues. Pour ce faire, il est avantageux d'appliquer la fonction cosinus 20 inverse à l'équation : I(x, y) = cos(27c6 (PM1 (x) + aMi (x) + PM2 (y) + aM2 (x))) , Dans un premier temps, on cherche donc à estimer les sommes croisées aMi(x)+aM2(y) des composantes inconnues. Pour estimer ces sommes croisées, l'équation précédente est alors 25 modifiée selon : I(x, y) = cos(2ir N +pp + 2ir (aMi (x) + aM2 (y))) où l'on a posé : N = Int(6 (pMi(x) + pM 2 (Y)) + 1 / 2)) (p = 2it6 (p11 (x) + pM 2 (Y)) ù 2Nit 30 La fonction « Int » est ici la fonction partie entière, qui associe à un nombre réel l'entier qui lui est immédiatement inférieur ou égal. L'entier N et le réel cp sont connus. On pose alors : m = (p + 2ir (aMl (x, y) + aM2 (x, y)) L'équation à résoudre s'écrit donc : I (x, y) = cos(2ir N + m) Le résultat de la fonction cosinus inverse, que l'on va appliquer, dépend de la valeur de : m=qp+27ro(ami (x)+aM2(y)) La valeur de m peut être comprise dans quatre intervalles : - Si m est compris entre 0 et rr, alors cos-1(I(x,y))=m ; - Si m est compris entre û rr et 0, alors cos-1(I(x,y))=-m ; - Si m est supérieur à rr et inférieur à 27, alors cos-1(I(x,y))=2rr-m ; - Si m est supérieur à -27 et inférieur à -rr, on obtient alors la relation cos-1(I(x,y))=2rr+m. Néanmoins, m n'est pas connu, seul cp est connu. On calcule alors les quatre résultats possibles de la fonction cosinus inverse, et l'on choisit la valeur dont l'écart avec cp est minimal. Dans la grande majorité des cas, cela permet de détermineraMi(x)+aM2(y) . On peut avantageusement choisir de filtrer les résultats incertains, c'est-à-dire ne pas considérer les points pour lesquels la valeur du premier minimum (écart minimum entre le résultat du cosinus inverse et la valeur de cp) est proche de celle du second minimum (second écart minimum entre le résultat du cosinus inverse et la valeur de cp). Dans le cas où cp est proche de 0 ou de rr ou de -rr, la valeur choisie n'est en générale pas la bonne. Les indéterminations en 0 ne peuvent pas être évitées. Les indéterminations en rr ou û rr sont uniquement causées par la différence entre cp et m, qui est d'autant plus faible que l'erreur de connaissance, c'est-à-dire la somme aM1(x)+aM2(y)des composantes inconnues, est faible. Il est donc possible d'éviter les indéterminations en rr ou û rr en minimisant 30 l'erreur de connaissance. Dans un deuxième temps, on estime les N+M composantes inconnues aMi,aM2 à partir des sommes croisées (aMi+aM2) déterminées par inversion du cosinus, pour en déduire les N+M premières et deuxièmes hauteurs de marche hM1, hM2. Il s'agit en fait d'une reconstruction mathématique de variables inconnues (ami, aM2) à partir de la connaissance de leurs sommes croisées aM1+aM2. Divers algorithmes peuvent être utilisés, comme par exemple des algorithmes de type réseaux de neurones, ou des procédés d'inversion matricielle comme une inversion par moindre carrés. On détaille ici une inversion par moindre carrés. Soit Y le vecteur des 10 sommes des composantes inconnues aMi(x)+aM2(y) . Soit A le vecteur des composantes inconnues aM1,aM2 à estimer. Dans le cas où l'on cherche en outre à estimer l'éventuelle inclinaison des surfaces réfléchissantes des miroirs, l'on décompose comme décrit précédemment les composantes inconnues ami, aM2 en la 15 somme de composantes axiales (a1M1 ;a1M2) et transversales (iM1 ;iM2), avec iMl (x, y) = bM (x).y et iM2 (x, y) = bM2 (y).x , les coefficients b1M1 et b1M2 étant des coefficients d'inclinaison. Dans ce cas, le vecteur A comprend les inconnues a1M1, a1M2, b1M1 et b1M2• Le problème à résoudre s'écrit donc sous forme matricielle : 20 Y=W*A On désigne par W la matrice de passage entre Y et A. Cette matrice est triviale, puisqu'il s'agit d'une matrice de passage entre un vecteur d'inconnues et un vecteur des sommes croisées de ces inconnues. Dans le cas où les composantes transversales (iM1 ;iM2) sont prises 25 comme égales à zéro, la matrice W est une matrice de 1 et de 0 qui permet de passer du vecteur des composantes inconnues (ami ;aM2) à leurs sommes croisées (ami + aM2). Si les composantes transversales (IM1 ;iM2) sont prises comme non nulles, la matrice W est une matrice comprenant des valeurs de x et y, ce qui permet de passer du vecteur des composantes 30 inconnues (a1M1, b1M1, ;a1M2, b1M2,,...) à leur sommes croisées (a1M1+b1Mi*y ; a1M2+b1M2*x...). Il est clair que le nombre d'inconnues augmente. Etant donné que nous disposons d'une estimation des sommes aMi(x)+aM2(y)des composantes inconnues, c'est à dire du vecteur Y, obtenue par l'inversion du cosinus précédemment décrite, il est possible d'estimer A. On peut par exemple estimer A en effectuant une inversion par 5 moindre carrés : A=(WTW +X,I)-l*WT *y, où I est la matrice identité, À un coefficient réel, et «T» la fonction de transposition d'une matrice. Le coefficient réel À est un coefficient de régularisation, car la matrice 10 WTW est en général mal conditionnée. II est déterminé empiriquement et est typiquement de l'ordre de 10-7. Dans le cas simple où l'on suppose les marches parfaitement parallèles (composantes transversales iM1,iM2 nulles), alors la matrice pseudo-inverse de W s'écrit de façon suivante: pour chacune de ses lignes (la ligne i correspondant au miroir Si) on a des 15 valeurs nulles sauf pour les M colonnes correspondant aux différences de chemin optique formées par la combinaison du miroir Si et d'un miroir S'i. Les valeurs non nulles sont toutes égales à 1/N. A l'issue de cette inversion par moindre carrées, on dispose d'une estimation des composantes inconnues aM1,aM2. On en déduit 20 immédiatement les hauteurs de marche hM1, hM2, par les relations hM1= pM1+aM1 et hM2=pM2+aM2, ce qui était l'objectif de l'étape d'estimation 22. A l'étape 23, on en déduit par conséquent les différences de chemin optique, par la relation 6(x,y) = hM1(x) + hM2(y). Dans ce deuxième mode de réalisation, il est très avantageux d'itérer 25 l'étape d'estimation des N+M composantes inconnues en intégrant à la composante connue (pM1 ;pM2) utilisée pour l'itération courante la composante inconnue (ami ;aM2) déterminée à l'itération précédente de ce mode de réalisation. En effet, les résultats issus de l'inversion du cosinus sont d'autant plus précis que l'erreur de connaissance est faible. 30 Des simulations ont montré dans le cas d'un faisceau 2 incident monochromatique de longueur d'onde À=10pm, et en supposant l'interférogramme d'intensité I(x,y) non bruité, que deux itérations permettent de déterminer les différences de chemin optique 6(x,y) avec une erreur de l'ordre du nanomètre, voire inférieure. Dans ce deuxième mode de réalisation, on met en oeuvre au cours de l'étape 22 une estimation des sommes croisées aM,+aM2 des composantes inconnues par l'inversion du cosinus dans l'expression I(x,y)=cos(2rr.6.(pMl+pM2+aM1+aM2)). Cette estimation peut également se faire sous certaines conditions par linéarisation de l'équation I(x,y)=cos(2rr.6.(pMl+pM2+aM,+aM2)). Il est alors nécessaire que le terme 2rr.6.(aMl +aM2) soit petit devant le chiffre 1. Par ailleurs, il est parfois nécessaire d'estimer l'amplitude du faisceau 2 incident monochromatique, que l'on ne connaît pas toujours parfaitement. Cette estimation est requise dès lors que l'amplitude n'est pas connue à mieux que 0,1%. Le procédé selon l'invention comprend alors une étape 24 d'estimation de l'amplitude, nommée a, du faisceau 2 incident monochromatique. En général, cette étape d'estimation de l'amplitude a du faisceau est mise en oeuvre en minimisant l'écart entre l'intensité I(x,y) de l'interférogramme et l'intensité modélisée en utilisant les différences de chemin optique 6(x,y) connues a priori ou estimées par le procédé selon l'invention. L'étape 24 d'estimation de l'amplitude a du faisceau 2 monochromatique peut en particulier être effectuée en minimisant la fonction f(a) = I(x,y) ù a.cos(2rr.6.6(x,y)), où a est l'amplitude du faisceau que l'on cherche à estimer, a le nombre d'onde du faisceau et 6(x,y) la différence de chemin optique. Avantageusement, on utilise la différence de chemin optique 6(x,y) estimée par le procédé à l'issue de l'étape 23. En effet, c'est à l'issue cette étape que l'on dispose de la meilleure connaissance des différences de chemin optique 6(x,y). Il a été constaté que le procédé selon l'invention permet d'obtenir une précision de mesure améliorée des différences de chemin optique 6(x,y) même dans le cas d'un interférogramme d'intensité I(x,y) très bruité (rapport entre le niveau de la ligne de base et le bruit de l'interférogramme égal à 100), et d'un faisceau 2 incident monochromatique à basse longueur d'onde. Dans ce cas, la précision reste de l'ordre de quelques nanomètres, y compris en présence d'un tilt des échelettes. Le procédé de métrologie selon l'invention trouve de nombreuses applications et notamment dans le domaine spatial, et permet alors d'obtenir une connaissance des différences de chemin optique précise et mise à jour en fonction de l'évolution des conditions de vol de l'aéronef embarquant l'interféromètre. Cette connaissance des différences de chemin optique permet d'effectuer diverses observations grâce à l'interféromètre : caractérisation de la présence de gaz dans l'atmosphère (profil de présence de gaz de type CO2, ozone, CO, CH4...) en spectroscopie infrarouge, mesure de pics de pollutions.... Ces mesures sont obtenues après traitement des interférogrammes mesurés par l'interféromètre, ledit traitement nécessitant la mesure des différences de chemin optique en vol et/ou au sol par le procédé de métrologie selon l'invention. Ainsi, l'utilisation d'un interféromètre mettant en oeuvre le procédé de métrologie selon l'invention, notamment à bord d'un satellite, permet une mise à jour et une connaissance fine des différences de chemin optique, et donc des mesures de spectroscopie précises. Le procédé de métrologie selon l'invention est rapide, précis, et nécessite une instrumentation beaucoup moins lourde et complexe que les procédés de l'art antérieur, ce qui facilite son utilisation en vol.25 REVENDICATIONS1. Procédé de métrologie d'un interféromètre (1) statique comprenant une lame (3) séparatrice apte à séparer un faisceau incident (2) vers 5 un premier (Ml) et un deuxième miroirs (M2), le premier miroir (Ml) présentant N surfaces (SI, S2, S3... ,SN) réfléchissantes sensiblement planes et parallèles entre elles, lesdites surfaces étant décalées entre elles et réalisant N premières hauteurs de marche (hMl(x)), 10 le deuxième miroir (M2) présentant M surfaces (S'1, S'2, S'3... ,S'M) réfléchissantes sensiblement planes et parallèles entre elles, lesdites surfaces étant décalées entre elles et réalisant M deuxièmes hauteurs de marche (hM2(y)), les hauteurs de marche (hMl(x),hM2(y)) correspondant à la distance 15 entre chaque surface (SI, S2, S3... ,SN, S'1, S'2, S'3... ,S'M) et la lame (3), de sorte que les deux miroirs (M1,M2) fournissent une pluralité de faisceaux réfléchis présentant entre eux des différences de chemin optique (b(x,y)) égales aux N*M sommes croisées des premières hauteurs de marche (hMl(x)) et des deuxièmes hauteurs de marche (hM2(y)), 20 ledit procédé étant caractérisé en ce qu'il comprend au moins les étapes consistant à : - acquérir un interférogramme d'intensité (I(x,y)) à partir d'un faisceau (2) incident monochromatique, - estimer (22) lesdites N+M premières et deuxièmes hauteurs de 25 marche (hMl(x);hM2(y)) par résolution d'un système d'équations liant l'intensité I(x,y) mesurée aux N+M premières et deuxièmes hauteurs de marche, ledit système exploitant les corrélations et dépendances géométriques des miroirs, et - déduire (23) les N*M différences de chemin optique (b(x,y)) à partir 30 des N+M premières et deuxièmes hauteurs de marche. 2. Procédé selon la revendication 1, dans lequel le système d'équations liant l'intensité I(x,y) aux premières et deuxièmes hauteurs de marche hMm, 22hM2 s'écrit I(x,y)=cos(2rr.a.(hM1+hM2)), où a est le nombre d'onde du faisceau (2) incident monochromatique. 3. Procédé selon l'une des revendications 1 ou 2, dans lequel les premières et deuxièmes hauteurs de marche (hM1(x);hM2(y)) sont chacune prises comme égales à la somme d'une composante connue (pM1;pM2) et d'une composante inconnue (aM1;aM2), et en ce que l'on estime les N+M composantes inconnues (aM1;aM2) par résolution dudit système. 4. Procédé selon la revendication 3, dans lequel on itère l'étape d'estimation des N+M composantes inconnues en intégrant à la composante connue (pM1 ;pM2) utilisée pour l'itération courante la composante inconnue (ami ;aM2) déterminée à l'itération précédente. 5. Procédé selon l'une des revendications 3 à 4, dans lequel la composante inconnue (ami ;aM2) est la somme d'une composante axiale (a1M1 ;a1M2) selon la direction (x ; y) du faisceau (2) incident qui tombe sur ce miroir, et d'une composante transversale (iM1 ;iM2) supérieure ou égale à zéro représentant les inclinaisons par rapport à au moins une des directions de l'espace (x,y) des surfaces (Si, S2, S3... ,SM) réfléchissantes du premier miroir (M1) et/ou des surfaces (S'1, S'2, S'3... ,S'M) réfléchissantes du deuxième miroir (M2), l'étape d'estimation (22) des premières et deuxièmes hauteurs de 25 marche comprenant l'estimation desdites composantes axiales et transversales. 6. Procédé selon la revendication 5, dans lequel chaque composante transversale (iM1 ;iM2) est modélisée selon : 30 iMl (x, y) = bM (x).y et iM2 (x, y) = bM2 (Y)•x , où (bM1 , bM2) sont des coefficients d'inclinaisons supérieurs ou égaux à zéro. 7. Procédé selon l'une des revendications 3 à 6, dans lequel on met en oeuvre pour résoudre ledit système l'étape consistant à estimer les sommes croisées (aM,+aM2) des composantes inconnues 5 à partir de l'équation I(x,y) = cos(2rr.6(pMl+pM2+aM,+aM2)), pour en déduire une estimation (22) desdites N+M premières et deuxièmes hauteurs de marche. 8. Procédé selon la revendication 7, dans lequel l'estimation des sommes 10 (aM,+aM2) des composantes inconnues est mise en oeuvre en appliquant la fonction cosinus inverse à l'équation I(x,y)=cos(2rr.6.(pMl+pM2+aM,+aM2)). 9. Procédé selon l'une des revendications 7 à 8, dans lequel on estime les composantes inconnues aMi,aM2 représentées par un vecteur A, à partir de 15 l'estimation des sommes (aMi+aM2) desdites composantes inconnues, représentées par un vecteur Y, en inversant l'équation matricielle Y = W A, où W est la matrice de passage entre Y et A. 10. Procédé selon l'une des revendications 3 à 6, dans lequel l'estimation 20 des composantes inconnues ami, aM2 est mis en oeuvre en effectuant un développement limité de l'équation I(x,y)=cos(2rr.6.(pMi+pM2+aM,+aM2)) autour de l'infiniment petit 2rr.6.(aMl+aM2). 11. Procédé selon l'une des revendications 1 à 10, comprenant en outre 25 une étape (24) d'estimation de l'amplitude a du faisceau (2) incident monochromatique. 12. Procédé selon la revendication 11, dans lequel l'étape (24) d'estimation de l'amplitude a du faisceau (2) incident monochromatique de nombre 30 d'onde a est mise en oeuvre en minimisant selon a la fonction f(a) = I(x,y) ù a.cos(2rr.6.6(x,y)). 13. Procédé selon l'une des revendications 1 à 12, dans lequell'interféromètre étant utilisé pour analyser un rayonnement compris dans une bande spectrale, le faisceau (2) incident est émis par une source (4) monochromatique de longueur d'onde comprise dans ladite bande spectrale.5