Brevet FR2938639A1 — publié le 20100521

Le domaine de l'invention est celui de la calibration d'une matrice de détection pour la correction des variations de gain entre pixels. L'invention concerne plus précisément l'interférométrie et un procédé qui permet de réduire les effets des erreurs causées par les variations de gain entre pixels de la matrice de détection dans le plan focal d'un interféromètre, erreurs qui viennent dégrader la précision des mesures réalisées par l'interféromètre. En spectrométrie par Transformée de Fourier Statique, l'interféromètre élabore un interférogramme en créant un nombre fini de différences de marche optique réparties spatialement, et acquises par une matrice de détection. Ceci peut être réalisé à l'aide d'une multitude de miroirs en parallèle, appelés échelettes, et d'un système optique imageant ces échelettes sur une matrice. Ceci peut aussi être réalisé à l'aide d'un léger tilt entre les deux miroirs, la matrice effectuant alors l'acquisition des franges d'interférence créées par ce tilt entre les deux miroirs. Dans tous les cas, les gains de détection de la matrice de détection varient d'un pixel à l'autre, ce qui induit une distorsion lors de l'acquisition de chaque élément de l'interférogramme. On cherche donc à corriger l'écart de gain interpixels le plus précisément possible. La méthode utilisée en imagerie classique pour calibrer les réponses interpixels d'une matrice de détection consiste à envoyer un signal uniforme (ou connu) sur tous les pixels de la matrice de détection. Cela peut se faire par exemple en éclairant l'instrument avec une sphère intégrante, qui fournit alors un éclairage uniforme. La réponse des pixels du plan focal doit alors être uniforme. Si ce n'est pas le cas, on calcule pour chaque pixel le gain relatif nécessaire pour que chacun réponde de la même façon à cet éclairage uniforme. Mais dans un interféromètre, le signal uniforme en entrée ressort modifié par les interférences. Le signal n'étant plus spatialement uniforme, il devient alors difficile de calibrer la matrice de détection en sortie de l'interféromètre. La méthode classique d'éclairage uniforme n'est donc pas satisfaisante dans le cadre de l'interférométrie. Il existe donc un besoin pour une méthode de calibration des défauts de gain entre pixels d'une matrice de détection qui permette de corriger avec précision l'écart de gain interpixels dans le cadre d'une application interférométrique. L'invention a pour objectif de répondre à ce besoin, et propose à cet effet selon un premier aspect, un procédé de calibration d'une matrice de détection utilisée pour acquérir une image d'un interférogramme, comprenant les étapes selon lesquelles : - on élabore un premier interférogramme et au moins deux interférogrammes décalés, le premier interférogramme étant élaboré en créant un nombre fini de différences de marche, et les interférogrammes décalés étant élaborés, pour chacun d'entre eux, par décalage de chacune des différences de marche du premier interférogramme, de telle sorte qu'on dispose pour chaque différence de marche d'au moins trois échantillons positionnés sur une portion de sinusoïde ; - on détermine par interpolation pour chaque portion de sinusoïde une fonction coïncidant avec la portion de sinusoïde; - on détermine un coefficient de calibration pour une différence de marche à partir de la fonction coïncidant avec la portion de sinusoïde correspondant à cette différence de marche. Certains aspects préférés, mais non limitatifs, de ce procédé sont les suivants : ù le coefficient de calibration correspond à la moyenne de la fonction ; ù la fonction est sinusoïdale ; ù la fonction correspond au produit d'une fonction sinusoïdale et d'une fonction polynomiale ; ù les au moins trois échantillons correspondent à des différences de 3o marche régulièrement espacées ; les au moins trois échantillons correspondent à des différences de marche espacées de manière à présenter un décalage total supérieur au quart d'une longueur d'onde du spectre en entrée de l'instrument élaborant les interférogrammes ; les au moins trois échantillons correspondent à des différences de marche espacées de manière à présenter un décalage total d'une longueur d'onde du spectre en entrée de l'instrument ; on élabore trois interférogrammes décalés ; on déplace un élément mobile de l'instrument de mesure pour 1 o élaborer un interférogramme décalé ; l'élément mobile est monté sur un mécanisme piézoélectrique ; le spectre en entrée de l'instrument élaborant les interférogrammes est à bande étroite. Selon un second aspect, l'invention concerne un instrument de mesure 15 comprenant une matrice de détection pour acquérir une image d'un interférogramme calibrée conformément au procédé selon le premier aspect de l'invention. D'autres aspects, buts et avantages de la présente invention apparaîtront mieux à la lecture de la description détaillée suivante de formes 20 de réalisation préférées de celle-ci, donnée à titre d'exemple non limitatif et faite en référence à la figure unique annexée représentant l'acquisition d'interférogrammes décalés. Selon un premier aspect, l'invention concerne un procédé de calibration d'une matrice de détection utilisée pour acquérir une image d'un 25 interférogramme, l'interférogramme étant élaboré par un instrument de mesure, par exemple par un interféromètre statique à échelettes. Un exemple non limitatif d'application de l'invention est celui de l'interférométrie pour l'obtention d'un spectre dans l'infrarouge. Il s'agit par exemple de surveiller la pollution atmosphérique, notamment via l'obtention 30 des profils de concentration pour l'ozone et le monoxyde de carbone. Deux bandes spectrales étroites sont alors particulièrement visées : la bande B1 = 4 [1020cm-1, 1080cm-1] pour l'ozone et la bande B2 = [2132cm-', 2192cm-1] pour le monoxyde de carbone. D'une manière générale, et prenant l'exemple d'un interféromètre à échelettes, l'invention propose d'acquérir pour chaque échelette trois points (au moins) espacés à peu près régulièrement qui permettent de remonter à l'oscillation sinus, et donc de retrouver la valeur moyenne du signal. En effet, en faisant varier la différence de marche, on se déplace localement sur une portion de sinusoïde : le signal interférométrique varie effectivement de manière quasi-sinusoïdale lorsque l'on fait varier la différence de marche, et 1 o ce d'autant plus précisément que le spectre en entrée de l'instrument est à bande étroite. La période d'oscillation du signal interférométrique en fonction de la différence de marche est alors proche de la longueur d'onde centrale du spectre à bande étroite. En interpolant les points acquis par une fonction, on peut alors calculer un coefficient de calibration correspondant à l'échelette 15 et pouvant être appliqué à la matrice de détection pour corriger des écarts de gain interpixels. Ainsi le procédé selon le premier aspect de l'invention permet de calculer un coefficient de calibration pour chaque échelette, ce qui permet de déterminer le coefficient de calibration de chaque superpixel de la matrice de 20 détection. On décrit ci-après de manière plus précise une mise en oeuvre possible du procédé selon le premier aspect de l'invention. Selon une première étape, on élabore un premier interférogramme et au moins deux interférogrammes décalés. 25 Le premier interférogramme est élaboré en créant un nombre fini de différences de marche (chaque différence de marche correspondant à une échelette de l'interféromètre), tandis que les interférogrammes décalés sont élaborés successivement (dans le temps), pour chacun d'entre eux, par décalage de chacune des différences de marche du premier interférogramme 30 d'une distance connue. Considérant pour chaque différence de marche au moins trois points, notamment au moins trois points relativement régulièrement espacés le long d'une période de l'oscillation interférométrique, on dispose alors pour chaque différence de marche d'au moins trois échantillons placés sur une portion de sinusoïde (car placés le long de l'oscillation sinusoïdale du signal interférométrique). On retiendra qu'il n'y a pas de contrainte à espacer régulièrement les trois points ; ce qui importe étant que la distance entre les points soit connue. L'optimum consiste cependant en un espacement relativement régulier entre les points. Si deux points sont trop proches, ils risquent de ne plus apporter assez d'information pour résoudre le système d'équation dont il sera fait état par la suite. De manière avantageuse, les trois échantillons correspondent à des différences de marche distantes d'une longueur d'onde du spectre en entrée de l'instrument (4.63 pm au total dans le cas de la bande B2, soit lambda/2 entre chaque point). Un décalage total de l'ordre de grandeur d'une période (donc de lambda), entre le premier et le dernier échantillon est préféré car cela donne une meilleure estimation des paramètres recherchés. Mais un tel décalage n'est pas indispensable. En revanche, un décalage total trop faible (inférieur à lambda/4 par exemple) peut conduire à une mauvaise estimation des gains interpixels. Un décalage supérieur à lambda doit s'apprécier, pour chaque valeur, à modulo lambda. Selon une mise en oeuvre préférentielle de l'invention retenue dans la suite de la description, on élabore quatre interférogrammes décalés de manière à disposer de quatre échantillons d'interférogramme par échelette. Un interférogramme I est défini de la manière suivante : I(x)=LB(6).SF(6).(1+cos(27t .x)), OÙ û B est le spectre d'entrée de l'instrument (il correspond par exemple au spectre atmosphérique mesuré par l'instrument), en W/m2/sr/cm-1 ; û SF est le filtre de l'instrument, dans les simulations il correspond au filtre étroit. II définit le domaine spectral de mesure de l'instrument. C'est la présence de ce filtre étroit qui rend le signal sinusoïdal en fonction d'une variation de la différence de marche de l'ordre de grandeur de la longueur d'onde. ù x est la Différence de Marche (DDM), en cm ; ù a est le nombre d'onde en cm-1. On note Il, 12, 13 et 14 les 4 interférogrammes échantillonnés aux différences de marche suivantes conformément à la première étape du procédé selon le premier aspect de l'invention. Pour Il, x = DDM1 = [-0.13cm : 80pm : +8.06 cm] ; les différences de marche sont échantillonnées au pas moyen de 80pm entre -0.13 cm et +8.06 cm. Pour 12, x = DDM2 = DDM1+2.31 pm ; Pour 13, x = DDM3 = DDM1+4.05 pm ; Pour 14, x = DDM4= DDM1+6.36 pm. On relèvera que ces valeurs ne sont pas parfaitement régulières dans cet exemple, du fait notamment d'autres contraintes instrumentales indépendantes (deux bandes spectrales sont acquises par l'instrument). Pour chaque interférogramme, on dispose de 1024 valeurs. Notons i=1 1024. On relèvera que du fait des erreurs de réalisation des échelettes, l'échantillonnage n'est pas parfaitement régulier. On considère typiquement une erreur de réalisation de type bruit blanc d'amplitude +/-5 pm. A la suite de cette première étape, on dispose pour chaque différence de marche d'au moins trois échantillons placés sur une portion de sinusoïde. Reprenant l'exemple ici présenté, on considère ainsi pour chaque échelette i (correspondant à une différence de marche x;) quatre points régulièrement espacés permettant l'acquisition de quatre échantillons 11(x;), I2(x;),13(x;), et I4(xi) positionnés sur une portion de sinusoïde. On a représenté sur la figure 1 l'acquisition de ces échantillons respectivement pour une échelette i et une échelette i+k. Les échantillons correspondant à l'échelette i correspondent ainsi aux différences de marche DDM1(i), DDM2(i)= DDM1(i)+2.31 pm, DDM3(i)= DDM2(i)+1.74pm et DDM4(i)= DDM3(i)+2.31 pm. Les échantillons correspondant à l'échelette i+k correspondent quant à eux aux différences de marche DDM1(i+k), DDM2(i+k)= DDM1(i+k)+1.74pm, DDM3(i+k)= DDM2(i+k)+2.31 pm et DDM4(i+k)= DDM3(i+k)+2.31 pm, avec DDM1(i+k)= DDM1(i) +k.Te + alpha(k) où Te représente le pas moyen d'échantillonnage (80 pm dans l'exemple ici considéré) et alpha(k) l'erreur de régularité (connue) du pas d'échantillonnage. 1 o Le signal interférométrique variant de manière quasi-sinusoïdale lorsque l'on fait varier la différence de marche, les quatre échantillons acquis pour une échelette sont positionnées sur une portion de sinusoïde comme cela est représenté sur la figure 1. Comme on travaille sur une bande spectrale étroite centrée autour d'un 15 nombre a0, on peut alors écrire le spectre vu par l'instrument sous la forme du produit de convolution entre une fonction étroite (qui est le spectre étroit décalé en 0) centrée en 0 et un dirac centré en a0 : S(6)=B(6).SF(a) =B(6).*S0(6)Oa(6 -60) =C(6)Oa(6 -60) Or l'interférogramme est la transformée de Fourier cosinus du spectre 20 vu par l'instrument. I(x)=TFcos[C(6)Oa(6 -60)] =TFcos[C(6)]*TFcos(a(6 -60)) = TF cos[C(6 )] * cos(2i6 0 .x) L'interférogramme est donc localement le produit d'une fonction cosinus de période a0 et d'une fonction qui est l'inverse d'une fonction étroite, donc une fonction variant lentement en fonction de x. 25 On peut donc interpoler localement l'interférogramme par une fonction sinus. On considère alors que la fonction lentement variable est constante sur cet intervalle. Si on veut affiner l'évaluation (notamment dans le cas d'un spectre à bande plus large), on peut approximer la fonction lentement variable par un polynôme d'ordre adéquat (1 ou 2 en pratique suffit). Nous détaillons ci-après le cas où l'on considère la fonction lentement variable comme constante. Au cours d'une seconde étape, on détermine par interpolation les paramètres d'une fonction (typiquement une fonction sinusoïdale) coïncidant avec chaque portion de sinusoïde. On définit par exemple la fonction sinusoïdale suivante : F(xi)=Mi + C;.cos(27t6;'x;); où M; correspond au niveau moyen de la fonction ; C; correspond au contraste ; 6;' est le nombre d'onde en cm-' i correspond à l'échelette. x; est laDDM Dès lors, pour chaque échelette i, on cherche les paramètres de la fonction F (Mi, Ci et ai') tels que F(xi)=I(xi) et cela pour les quatre interférogrammes élaborés. Les paramètres sont par exemple obtenus par mise en oeuvre de l'algorithme des moindres carrés pour minimiser l'écart F(xi)ûI(xi) z On retiendra que dans le cadre de cette variante préférentielle, on dispose de quatre équations pour déterminer trois inconnues. Notons qu'il est possible de représenter la fonction lentement variable par un polynôme d'ordre 1 (une pente) pour améliorer le calcul, la fonction F est alors définie par : F(x;)= [ M; + C;.cos(27t6;'x;) ] * pi*(xi û xi). On dispose alors de quatre équations pour déterminer 4 inconnues. Au final, pour chaque échelette i, à partir des quatre valeurs d'interférogramme (I1(i), I2(i), I3(i), I4(i)), on obtient les inconnues Mi, Ci et ai'. Au cours d'une troisième étape, on détermine un coefficient de 3o calibration pour une différence de marche en utilisant les paramètres de la fonction coïncidant avec la portion de sinusoïde correspondant à cette différence de marche. Le coefficient de calibration correspond à une différence de marche est par exemple la moyenne de cette fonction. Revenant à l'exemple de l'interpolation par la fonction sinusoïdale F présentée ci-dessus, le coefficient de calibration pour la différence de marche i correspond ainsi au paramètre Mi. La Demanderesse a évalué l'écart entre l'interférogramme de départ Ij et l'interférogramme calculé par interpolation Icj, et cela pour chaque 1 o différence de marche x;. Ic1(x;)=F(x;=DDMI(i)) = M; + C;.cos(2n6;'x;) Ic2(x;)=F(x;=DDM2(i)) = M; + C;.cos(2n6;'x;) Ic3(x;)=F(x;=DDM3(i)) = M; + C;.cos(2n6;'x;) Ic4(x;)=F(x;=DDM4(i)) = M; + C;.cos(2n6;'x;) 15 La Demanderesse a pu constater que le rapport signal sur bruit (le bruit correspondant à l'écart type de la différence entre interférogrammes de départ et interférogrammes calculés) est supérieur à la spécification en vigueur dans le domaine d'application de l'invention. Des écarts de performance sont cependant notés pour les faibles 20 différences de marche (typiquement les échelettes 1 à 67) où le contraste interférométrique est important. Ainsi en prenant en compte ces premières échelettes, on diminue d'un facteur 8 le rapport signal sur bruit de l'interférogramme. La Demanderesse a ainsi pu évaluer que le bruit sur l'échelette 67 est environ 500 fois plus important que pour l'échelette 601. Il 25 est possible d'améliorer (lorsque cela est nécessaire) le résultat en prenant en compte une fonction d'ordre plus élevée pour la fonction TFcos[C(a)]. Les résultats sont en revanche particulièrement bons sur 95% de l'interférogramme (pour les grandes différences de marche, typiquement les échelettes 68 à 1024). Sur un interférogramme sans bruit, le procédé selon 30 le premier aspect de l'invention permet ainsi de connaître les variations de gain de la matrice de détection pour les échelettes 68 à 1024 avec une précision de 0.003%. Une fois la calibration de la matrice de détection réalisée (par calcul des coefficients de calibration pour l'ensemble des échelettes), lors de l'élaboration d'un nouvel interférogramme, on vient corriger l'image de l'interférogramme acquise par la matrice de détection à l'aide desdits coefficients de calibration (plus précisément on vient corriger chacun des superpixels de la matrice correspondant à une échelette du coefficient de calibration associé à cette échelette calculé lors de la calibration). 1 o Au cours de la première étape du procédé selon l'invention, on réalise l'acquisition d'un premier interférogramme et d'au moins deux interférogrammes décalés, chacun des interférogrammes décalés correspondant à des échantillons décalés d'une quantité donnée par rapport aux échantillons du premier interférogramme. Selon un mode de réalisation 15 possible de l'invention, on déplace un élément mobile de l'instrument de mesure pour élaborer un interférogramme décalé. Afin de générer une variation de signal connu pour chaque (super)pixel de la matrice de détection et élaborer ainsi un interférogramme décalé, on peut faire varier faiblement la différence de marche sur l'un des deux bras de 20 l'interféromètre, en positionnant l'un des miroirs de l'interféromètre sur un mécanisme piézo. Alternativement, on utilise un faible mouvement d'inclinaison (mouvement de tilt) ou un faible mouvement de translation d'une lame de verre (prismée dans le cas d'une translation), éventuellement déportée dans 25 un bras de l'interféromètre. Cette lame de verre peut être la lame compensatrice de l'interféromètre. Ceci permet d'obtenir un mouvement très pur mécaniquement, avec un dispositif de déplacement à base de piézo-électrique. On relèvera que l'invention n'est pas limitée à un procédé selon son 30 premier aspect mais s'étend également à un instrument de mesure comprenant une matrice de détection pour acquérir une image d'un 11 interférogramme calibrée conformément au procédé selon le premier aspect de l'invention. Les avantages de l'invention par rapport à l'existant sont notamment les suivants : - Calibration pouvant être réalisée dans un instrument embarqué en vol (par exemple un interféromètre embarqué dans un satellite) - Fiabilité accrue de la calibration - Très bonne précision, aucun élément externe (cache noir) n'est nécessaire, et toute la chaîne instrumentale est prise en compte. - Modulation de phase possible (pour les spectromètres acquérant une bande étroite). - Dans le cas d'un spectromètre à Transformée de Fourier statique à échelettes, la calibration peut être réalisée directement sur la scène observée. A ce titre, un dispositif piézoélectrique s'avère plus fiable qu'un dispositif à cache amovible qui viendrait s'intercaler alternativement sur les deux voies de l'interféromètre pour l'obturer. Un tel dispositif à cache amovible présente effectivement deux mécanismes, à priori volumineux. Par ailleurs, les caches intercalés peuvent avoir eux même des défauts s'ils ne sont pas parfaitement noirs. On notera en outre que dans l'infrarouge thermique, ces caches émettraient de la lumière, ce qui ne rendrait pas possible une calibration. Enfin, dans le cadre d'une utilisation opérationnelle, un tel dispositif à cache amovible présenterait une criticité importante et une perte de la mission en cas de panne en position obturée. Une calibration éléments par éléments présente quant à elle l'inconvénient de ne pouvoir être réalisée qu'à l'assemblage. Ensuite, un suivi de cette calibration par un dispositif partiel (éclairant uniquement la matrice par exemple) pourrait permettre de prendre en compte une partie des évolutions. Mais ce dispositif ne permet donc pas de prendre en compte 3o toutes les évolutions susceptibles d'arriver au cours de la vie de l'instrument (poussières se déposant en entrée de l'instrument par exemple). 12 Dans le cas d'une calibration sur un corps noir embarqué (ou sphère intégrante � blanche �), l'uniformité de l'éclairage est modifiée par les interférences, mais celles-ci proviennent d'une scène connue (le corps noir) et il est donc possible d'étalonner les défauts du détecteur à partir de cette scène connue à défaut d'être uniforme. Cependant, ce système est forcément présent sur une voie annexe de l'interféromètre. Il présente donc des défauts qui ne sont pas présents sur la voie de mesure de l'interféromètre (poussières...) ; à l'inverse, il ne passe pas par tous les éléments de l'interféromètre (miroir en entrée...) et ne calibre donc pas tous les défauts. Enfin, ce système nécessite de connaître parfaitement le profil spectral de la source d'éclairement, et son uniformité, qui ne peut être meilleure que quelques pourcents. Or on souhaite un étalonnage typiquement de l'ordre de 0.1%.15 REVENDICATIONS1. Procédé de calibration d'une matrice de détection utilisée pour acquérir une image d'un interférogramme, comprenant les étapes selon lesquelles : on élabore un premier interférogramme et au moins deux interférogrammes décalés, le premier interférogramme étant élaboré en créant un nombre fini de différences de marche, et les interférogrammes décalés étant élaborés, pour chacun d'entre eux, par décalage de chacune 1 o des différences de marche du premier interférogramme, de telle sorte qu'on dispose pour chaque différence de marche d'au moins trois échantillons positionnés sur une portion de sinusoïde ; - on détermine par interpolation pour chaque portion de sinusoïde une fonction coïncidant avec la portion de sinusoïde; 15 - on détermine un coefficient de calibration pour une différence de marche à partir de la fonction coïncidant avec la portion de sinusoïde correspondant à cette différence de marche. 2. Procédé selon la revendication 1, dans lequel le coefficient de calibration 20 correspond à la moyenne de la fonction. 3. Procédé selon l'une des revendications 1 à 2, dans lequel la fonction est sinusoïdale. 25 4. Procédé selon l'une des revendications 1 à 2, dans lequel la fonction correspond au produit d'une fonction sinusoïdale et d'une fonction polynomiale. 5. Procédé selon l'une des revendications 1 à 3, dans lequel les au moins 30 trois échantillons correspondent à des différences de marche régulièrement espacées. 6. Procédé selon l'une des revendications 1 à 5, dans lequel les au moins trois échantillons correspondent à des différences de marche espacées de manière à présenter un décalage total supérieur au quart d'une longueur d'onde du spectre en entrée de l'instrument élaborant les interférogrammes. 7. Procédé selon la revendication 6, dans lequel les au moins trois échantillons correspondent à des différences de marche espacées de manière à présenter un décalage total d'une longueur d'onde du spectre en entrée de l'instrument. 8. Procédé selon l'une des revendications 1 à 7, dans lequel on élabore trois interférogrammes décalés. 9. Procédé selon l'une des revendications 1 à 8, dans lequel on déplace un 15 élément mobile de l'instrument de mesure pour élaborer un interférogramme décalé. 10. Procédé selon la revendication 9, dans lequel l'élément mobile est monté sur un mécanisme piézoélectrique. 11. Procédé selon l'une des revendications 1 à 10, dans lequel le spectre en entrée de l'instrument élaborant les interférogrammes est à bande étroite. 12. Instrument de mesure comprenant une matrice de détection pour 25 acquérir une image d'un interférogramme calibrée conformément au procédé selon l'une des revendications 1 à 11. 10 20