Brevet FR3021127A1 — publié le 20151120

B13227 - DD14758ST 1 SYSTÈME ET PROCÉDÉ D'ACQUISITION D'IMAGES HYPERSPECTRALES Domaine La présente demande concerne un système et un procédé d'acquisition d'images hyperspectrales ou multispectrales, et la calibration d'un tel système ou procédé. Exposé de l'art antérieur On appelle généralement image hyperspectrale ou multispectrale d'une scène une série de plusieurs images élémentaires en deux dimensions de la scène, chaque image représentant une composante de la scène dans une bande restreinte de longueurs d'ondes. Chaque image élémentaire correspond généralement à l'intégration de l'intensité lumineuse dans une bande spectrale spécifique. Dans la suite, l'expression "image hyperspectrale" est considérée comme étant équivalente à l'expression "image multispectrale", et désigne une série d'au moins deux images élémentaires d'une même scène dans des bandes de longueurs d'ondes distinctes. Pour acquérir une image hyperspectrale, une méthode consiste à acquérir successivement les différentes images de la série à l'aide d'un capteur d'images, en disposant à chaque acquisition, entre la scène et le capteur, un filtre optique passe bande de bande passante étroite, par exemple de quelques nanomètres à quelques dizaines de nanomètres, correspondant à B13227 - DD14758ST 2 l'une des composantes de l'image hyperspectrale. A chaque nouvelle acquisition, un nouveau filtre est disposé entre la scène et le capteur, de façon à sélectionner successivement les différentes bandes de longueurs d'ondes de l'image hyperspectrale. Un inconvénient majeur de cette méthode réside dans la nécessité de faire défiler successivement les différents filtres optiques entre la scène et le capteur. Il en résulte des systèmes d'acquisition relativement encombrants et des temps d'acquisition relativement longs. Une fois l'acquisition terminée, l'infor- mation contenue dans chaque image élémentaire correspond à l'intégration de l'intensité lumineuse dans la bande passante du filtre correspondant. Une autre méthode consiste à utiliser, à la place du filtre, un prisme, pour disperser spatialement les différentes longueurs d'ondes de la scène. Un système matriciel d'obturateurs programmables, du type couramment désigné dans la technique par l'acronyme SLM (de l'anglais "Spatial Light Modulator" - modulateur spatial de lumière) peut être placé entre la scène et le prisme, en amont du capteur d'images. Le SLM est par exemple constitué par une matrice de micro-miroirs ou par une matrice de cellules à cristaux liquides. En fonction des résolutions spatiale et spectrale recherchées, une acquisition ou un nombre restreint (c'est-à-dire inférieur au nombre total de bandes spectrales de l'image hyperspectrale) d'acquisitions peut suffire pour reconstruire l'image hyperspectrale. Ce type de méthode est généralement désignée dans la technique par l'acronyme CASSI (de l'anglais "Coded Aperture Snapshot Spectral Imaging" - imagerie spectrale instantanée par ouverture codée). Un inconvénient de cette méthode réside toutefois dans la complexité du système optique nécessaire pour réaliser l'acquisition, et dans la résolution, relativement faible, des images élémentaires qui peuvent être acquises simultanément par le capteur. D'autres méthodes ont été proposées, par exemple dans 35 l'article intitulé "Compressive sensing spectrometry based on B13227 - DD14758ST 3 liquid crystal devices" de Yitzhak August et al., utilisant les théories de l'acquisition compressive pour reconstruire les différentes composantes d'un spectre étendu à partir d'un nombre réduit d'acquisitions. Dans l'article susmentionné, un filtre à cristaux liquides commandable électriquement est disposé entre la scène et un capteur d'images. Lors d'une acquisition, une tension de commande est appliquée au filtre, de façon que ce dernier sélectionne non pas une bande passante étroite, comme dans les méthodes traditionnelles, mais une bande passante large de forme irrégulière, correspondant à une juxtaposition des longueurs d'ondes des différentes composantes de l'image hyperspectrale auxquelles sont affectés des coefficients de pondération. Plusieurs acquisitions sont successivement effectuées, en modifiant à chaque fois la commande du filtre à cristaux liquides, ce qui revient à modifier la réponse du filtre et donc les coefficients de pondération affectés aux composantes de l'image hyperspectrale. Les auteurs ont montré qu'en utilisant les modèles mathématiques de l'acquisition compressive, il est possible de reconstruire une série de composantes élémentaires du spectre étendu à partir d'un nombre d'acquisition inférieur au nombre total de bandes spectrales de la série reconstruite. Un inconvénient de cette méthode réside toutefois dans la complexité de la calibration des filtres à cristaux liquides, et dans les possibilités réduites en termes de compression ou de 25 super-résolution. Il existe un besoin pour un système et un procédé d'acquisition d'images hyperspectrales, palliant tout ou partie des inconvénients des solutions connues. Résumé 30 Ainsi, un mode de réalisation prévoit un système d'acquisition d'une image hyperspectrale, comportant : un capteur d'images en niveaux de gris ; et un élément diffuseur et dispersif placé sur le trajet optique entre le capteur et une scène, cet élément comprenant une matrice de cellules à cristaux liquides 35 commandables individuellement, dans lequel chaque cellule peut B13227 - DD14758ST 4 recevoir une tension de commande choisie parmi une série d'au moins trois tensions de commande distinctes. Selon un mode de réalisation, le système comporte en outre un module de commande adapté à commander le capteur et l'élément pour acquérir successivement un nombre entier M supérieur à 1 d'images élémentaires en niveaux de gris de la scène, en appliquant, lors de chaque acquisition, une commande spécifique à l'élément. Selon un mode de réalisation, pour chaque commande 10 appliquée à l'élément, l'élément présente une fonction d'étalement de point qui dépend de la longueur d'onde des rayons qu'il reçoit. Selon un mode de réalisation, les M commandes appliquées à l'élément diffuseur et dispersif sont telles que les MxK fonctions d'étalement de point de l'élément, correspondant aux M 15 commandes appliquées et à K bandes spectrales de l'image hyperspectrale, avec K entier supérieur à 1, soient irrégulières et soient toutes différentes les unes des autres. Selon un mode de réalisation, le système comporte en outre un module de traitement adapté à reconstruire l'image 20 hyperspectrale à partir des M images élémentaires acquises par le capteur. Un autre mode de réalisation prévoit un procédé de commande du système susmentionné, comportant une phase d'acquisition au cours de laquelle le capteur et l'élément sont 25 commandés pour acquérir successivement un nombre entier M supérieur à 1 d'images élémentaires en niveaux de gris de la scène, en appliquant, lors de chaque acquisition, une commande spécifique à l'élément. Selon un mode de réalisation, le procédé de commande 30 comporte en outre une phase de reconstruction de l'image hyperspectrale à partir des M images élémentaires acquises pendant la phase d'acquisition. Selon un mode de réalisation, le procédé de commande comporte en outre une phase préalable de calibration au cours de 35 laquelle les MxK fonctions d'étalement de point de l'élément, B13227 - DD14758ST correspondant aux M commandes appliquées et à K bandes spectrales de l'image hyperspectrale, avec K entier supérieur à 1, sont déterminées. Selon un mode de réalisation, la phase de calibration 5 comprend l'acquisition successivement par le capteur de MxK images de taches résultant de la diffusion, par l'élément, pour les M commandes de l'élément et pour les K bandes spectrales lumineuses de l'image hyperspectrale, d'une source lumineuse ponctuelle réglable en longueur d'onde. Selon un mode de réalisation, les MxK images acquises lors de la phase de calibration sont mises en correspondance avec un modèle comportemental théorique de l'élément. Brève description des dessins Ces caractéristiques et avantages, ainsi que d'autres, 15 seront exposés en détail dans la description suivante de modes de réalisation particuliers faite à titre non limitatif en relation avec les figures jointes parmi lesquelles : la figure 1 représente de façon schématique et partielle un exemple d'un mode de réalisation d'un système d'acquisition 20 d'images hyperspectrales ; la figure 2 est un schéma illustrant le principe de fonctionnement du système d'acquisition d'images hyperspectrales de la figure 1 ; et la figure 3 représente de façon schématique et partielle 25 le système d'acquisition d'images hyperspectrales de la figure 1 pendant une phase de calibration. Description détaillée Par souci de clarté, de mêmes éléments ont été désignés par de mêmes références aux différentes figures. 30 La figure 1 représente de façon schématique et partielle un exemple d'un mode de réalisation d'un système d'acquisition d'images hyperspectrales. Le système de la figure 1 comprend un capteur d'images 101, par exemple une matrice de photodétecteurs, adapté à fournir 35 des images en niveau de gris. Le capteur d'images 101 est un B13227 - DD14758ST 6 capteur large bande, c'est-à-dire qu'il est sensible à toutes les longueurs d'ondes des images hyperspectrales que l'on souhaite acquérir. Le système de la figure 1 comprend en outre, sur le trajet optique entre le capteur 101 et un objet ou une scène 103 (OBJ) dont on souhaite acquérir une image hyperspectrale, un élément diffuseur et dispersif programmable 105, comprenant une matrice de cellules à cristaux liquides commandables individuellement. Les cristaux liquides de l'élément 105 sont des cristaux liquides présentant un indice de réfraction dépendant de la longueur d'onde, et de préférence des cristaux liquides dont la dispersion en longueur d'onde varie en fonction de la tension de polarisation appliquée, par exemple des cristaux liquides du type commercialisé par la société MERCK sous la référence TL216 ou sous la référence E44. A titre d'exemple non limitatif, chaque cellule élémentaire peut présenter des dimensions de l'ordre de 1 à 10 gm pour des applications dans le domaine du visible. En pratique, l'élément 105 peut comprendre une couche continue de cristaux liquides dont une première face est revêtue par une électrode commune continue et dont une deuxième face opposée à la première face est recouverte par un réseau matriciel d'électrodes élémentaires discrètes commandables individuellement. Dans ce cas, une cellule élémentaire de l'élément 105 est constituée par la portion de la couche de cristaux liquides disposée entre une électrode élémentaire et l'électrode commune. Dans cet exemple, l'objet 103 dont on souhaite acquérir une image hyperspectrale est supposé être situé à l'infini. Le cas échéant, un système optique complémentaire, non représenté, peut être disposé entre l'objet et l'élément 105 de façon que l'objet 103 apparaisse à l'infini du point de vue de l'élément 105. Le système de la figure 1 peut en outre comprendre une ou plusieurs mémoires (non représentées) adaptées à stocker une ou plusieurs images numériques acquises par le capteur 101, un 35 module de calcul (non représenté), par exemple un microprocesseur, B13227 - DD14758ST 7 adapté à traiter des images acquises par le capteur 101, et/ou un module de commande du capteur 101 et de l'élément diffuseur et dispersif 105. Le fonctionnement du système de la figure 1 va 5 maintenant être décrit. Lors d'une phase d'acquisition d'une image hyperspectrale de l'objet 103, le capteur 101 acquière successivement M images ym de l'objet 103 vu à travers l'élément diffuseur et dispersif 105, où M est un nombre entier supérieur à 1 et m est 10 un nombre entier compris dans la plage allant de 1 à M. Lors de l'acquisition d'une image ym, l'élément diffuseur et dispersif 105 est modulé par une commande Cm, c'est-à-dire un jeu de signaux de commande des différentes cellules à cristaux liquides de l'élément 105. Chaque jeu de commande permet 15 de générer une carte spécifique d'indices de réfraction sur la matrice de cristaux liquides. Dans les modes de réalisation décrits, les cellules ne sont pas commandées en obturation / ouverture comme dans un SLM d'un système de type CASSI, mais chaque cellule peut recevoir une tension de commande choisie parmi 20 une série d'au moins trois tensions distinctes, et de préférence parmi une série d'au moins dix tensions distinctes. De préférence, l'élément 105 est tel que, pour une commande Cm donnée de l'élément 105, il n'existe pas de discontinuité d'indice de réfraction à la frontière entre deux cellules voisines de la matrice, de façon 25 que la couche de cristaux liquides de l'élément 105 soit perçue comme un élément non discret (c'est-à-dire continu ou analogique). Les M commandes Cm appliquées à l'élément 105 lors des M acquisitions d'images par le capteur 105 sont toutes différentes les unes des autres. A chaque commande Cm de l'élément 105, 30 l'élément 105 présente une fonction d'étalement du point ou PSF (de l'anglais "Point Spread Function") qui varie en fonction de la longueur d'onde des rayons reçus par l'élément 105 (caractère dispersif de l'élément 105). Pour une commande Cm donnée de l'élément 105 et pour une longueur d'onde donnée, la PSF de 35 l'élément 105 forme une tache de diffusion, de préférence B13227 - DD14758ST 8 irrégulière, pouvant s'étaler sur une pluralité de photodétecteurs du capteur 101 et de préférence sur la majeure partie de la surface du capteur 101 (caractère diffusant de l'élément 105). Si on définit la PSF de l'élément 105, à une longueur d'onde donnée, par l'image de la tâche, dans le plan du capteur 101, résultant de la diffusion par l'élément 105 d'un faisceau provenant d'une source lumineuse ponctuelle à cette longueur d'onde, les commandes Cm sont de préférence choisies de façon que les MxK PSF de l'élément 105 (où K est un nombre entier supérieur à 1 correspondant au nombre de bandes spectrales de l'image hyperspectrale que l'on souhaite acquérir) soient suffisamment décorrélées les unes des autres. A titre d'exemple, un groupe de M commandes de l'élément 105 peut être sélectionné parmi un groupe de commandes générées sur la base de polynômes de Zernike de façon à minimiser le maximum d'intercorrélation et le maximum d'autocorrélation des PSF. A titre d'exemple, un groupe de commandes adapté pour une longueur d'onde donnée peut être sélectionné, puis un biais adapté à minimiser la corrélation des PSF dans l'axe spectral peut être déterminé pour chaque commande du groupe. De préférence, pour une commande donnée, la déviation spectrale de la PSF entre deux bandes spectrales successives de l'image hyperspectrale à acquérir est inférieure ou égale à un pixel du capteur, de façon à éviter certains artefacts spectraux. Du point de vue de l'acquisition compressive, on peut en outre chercher à minimiser la cohérence mutuelle entre la matrice d'acquisition (relative aux PSF) et la base dans laquelle les images hyperspectrales que l'on souhaite acquérir sont parcimonieuses. Par ailleurs, pour définir des commandes permettant de générer des PSF orthogonales pour une longueur d'onde donnée, on peut effectuer une analyse en composantes principales non négatives d'un groupe de PSF générées à partir de commandes issues de polynômes de Zernike. Il est alors possible de revenir aux commandes associées aux PSF orthogonales à l'aide d'un algorithme permettant de retrouver les fronts d'onde correspondant aux PSF générées. B13227 - DD14758ST 9 Les commandes Cm sont de préférence choisies de façon que, pour une commande Cm donnée et pour une longueur d'onde donnée, la PSF de l'élément 105 soit sensiblement la même en tout point de l'élément 105. Si ce n'est pas le cas, une correction de la distorsion peut être incorporée à l'algorithme de reconstruction de l'image hyperspectrale. A chaque commande Cm de l'élément 105, différentes cellules à cristaux liquides de l'élément 105 sont commandées pour présenter des indices de réfraction différents (à longueur d'onde constante), de façon à créer des décalages de phase du front d'onde en sortie de l'élément 105, permettant de générer des taches de diffusion ou PSF irrégulières. A titre d'exemple les commandes Cm correspondent à des commandes permettant de générer des PSF connues, par exemple des PSF issues d'aberrations du front d'onde correspondant aux polynômes de Zernike ou à une combinaison de ces polynômes. On notera en particulier, à titre d'exemple non limitatif, le cas d'aberrations d'astigmatisme ou d'aberrations sphériques couplées à un biais ou tilt. Comme cela sera expliqué plus en détail ci-après, les inventeurs ont déterminé qu'à partir des M images floues en niveau de gris ym acquises par le capteur 101 lors de la phase d'acquisition, il est possible de reconstituer une image hyperspectrale x comportant K images élémentaires xxk correspondant aux composantes de l'image x dans les différentes bandes spectrales Xk de l'image x (avec k entier allant de 1 à K). Les calculs permettant de reconstruire l'image hyperspectrale x à partir des M images ym vont être détaillés ci-après. Ces calculs peuvent être mis en oeuvre soit à l'aide d'un module de calcul intégré au système d'acquisition, soit par une unité de traitement externe. On désigne ici par hm,xk le noyau de convolution correspondant à l'effet de la PSF de l'élément 105 sur la lumière de la bande spectrale Xk, lorsque l'élément 105 reçoit la commande Cm. Si on considère que les images élémentaires xxk recherchées ont une résolution PxQ supérieure d'un facteur s à la résolution B13227 - DD14758ST 10 PsxQs du capteur 101, avec P, Q, s, Ps et Qs entiers, Ps=P/s et 45=Q/s, et si on considère que les MxK PSF de l'élément 105 peuvent être définies avec un niveau de granularité ou pas de discrétisation identique à celui des images xxk, chaque image ym acquise par le capteur 101 peut être exprimée analytiquement par la formule suivante : yin - Ds* X2,k), k=1 où * représente l'opérateur de convolution (non cyclique), et où Ds désigne un opérateur de décimation par le facteur s, par exemple un opérateur de moyennisation, correspondant à la décimation spatiale, par le capteur 101, de l'image hyperspectrale x. La figure 2 représente une modélisation schématique de cette équation. Dans un premier temps, on considère que l'on cherche à obtenir une image hyperspectrale x à la même résolution spatiale PsxQs que le capteur 101, et que les MxK PSF de l'élément 105 peuvent être définies avec un niveau de granularité correspondant à la résolution Ps x Qs du capteur. L'opérateur Ds peut donc être omis (cas s = 1). Les opérations de convolution dans le domaine spatial correspondent à des multiplications point à point ou produits de Hadamard dans le domaine de Fourier. L'équation susmentionnée peut alors être approximée comme suit : Yin - k=1 où '.' représente le produit de Hadamard, Ym est la transformée de Fourier de l'image ym, Hm,xk est la transformée de Fourier du noyau de convolution hm,xk (éventuellement complété par des zéro si les dimensions du noyau hm,xk sont inférieures aux dimensions de l'image xxk), et Xxk est la transformée de Fourier de l'image élémentaire xxk. On notera que les effets de bords peuvent, le cas échéant, être compensés en utilisant des méthodes de type B13227 - DD14758ST 11 zero-padding, c'est-à-dire en augmentant la taille des matrice et en complétant par des zéros. Pour chaque pixel de coordonnées spatiales p,q de l'image hyperspectrale x, où p est un entier allant de 1 à P et q 5 est un entier allant de 1 à Q, cette équation peut s'exprimer sous la forme d'une simple multiplication matricielle comme suit : Y (p, g) = G (p, q)X(p, g), où Y(p,q) est un vecteur de dimension M, contenant les M valeurs du pixel de coordonnées p,q dans les M images en transformée de 10 Fourier Ym, X(p,q) est un vecteur de dimension K contenant les K valeurs du pixel de coordonnées p,q dans les K images en transformée de Fourier Xxk, et G(p,q) est une matrice de dimensions MxK contenant les MxK valeurs des pixels de coordonnées p,q dans les MxK transformées de Fourrier Hm,xk des noyaux de convolution 15 hm,u. On considère le cas où le nombre M d'acquisitions effectuées par le capteur 101 est égal au nombre K d'images élémentaires de l'image hyperspectrale x recherchée. Les matrices G(p,q) sont alors des matrices carrées. Si les MxK PSF de l'élément 20 105 sont suffisamment décorrélées les unes des autres, les matrices carrées G(p,q) peuvent être inversées ou pseudoinversées. Les M commandes Cm de l'élément 105 sont de préférence choisies de façon que les matrices carrées G(p,q) soient inversables. Chaque vecteur X(p,q) de l'image hyperspectrale en 25 transformée de Fourier peut directement être calculé par une simple multiplication matricielle selon la formule suivante : X (p, g) = G(p,q)-1-Y (p, g), où G(p,q)-1 désigne la matrice G(p,q) inversée ou pseudo-inversée. Les différentes images élémentaires xxk de l'image hyperspectrale 30 x peuvent alors être reconstruites par de simples opérations de transformée de Fourier inverse. Un avantage du système d'acquisition d'images hyperspectrales décrit ci-dessus est qu'il est particulièrement B13227 - DD14758ST 12 simple et peu encombrant. En particulier, il n'est pas nécessaire de prévoir un mécanisme permettant de changer mécaniquement le filtre optique situé entre le capteur et la scène entre les différentes acquisitions. De plus, la réalisation de l'élément 105 est relativement simple, puisque l'élément 105 ne doit pas réaliser une fonction optique traditionnelle précise de filtrage, mais seulement introduire un flou dans l'image, la seule contrainte étant que les MxK PSF de l'élément 105 présentent un bon niveau de décorrélation les unes par rapport aux autres. On notera que les matrices G(p,q)-1 peuvent être prédéterminées, et stockées dans une mémoire d'un module de reconstruction de l'image hyperspectrale. Les matrices G(p,q)-1 peuvent en outre optionnellement être préconditionnées pour éviter certaines erreurs liées aux effets de bord. Plus généralement, si le nombre de valeurs de pixels effectivement acquises est supérieur au nombre de valeurs de pixels à reconstruire, par exemple si le nombre M d'acquisitions effectuées par le capteur 101 est inférieur au nombre K d'images élémentaires de l'image hyperspectrale x que l'on cherche à reconstruire, et/ou si la résolution spatiale PxQ de l'image hyperspectrale x est supérieure à la résolution pxq du capteur 101 (cas s > 1), l'image x peut être reconstruite à l'aide d'une méthode de régularisation itérative du type utilisé dans le domaine de l'acquisition compressive, correspondant à la résolution d'un problème de minimisation formulé comme suit : M K 2\ argminx 1(x) + y Fshirt,2,1,* x2A, - Yrn m=1 k=1 2/ où J(x) correspond à un opérateur de régularisation permettant d'exacerber la structure interne de l'image hyperspectrale et y est un paramètre scalaire de régularisation. A titre d'exemple non limitatif, l'opérateur J(x) est basé sur cinq opérateurs de gradients V1, V2, V3, V4 et VS, correspondant à des contraintes qui s'appliquent au cube de dimensions PxQxK de l'image hyperspectrale x, respectivement dans les directions spatiales B13227 - DD14758ST 13 verticale (P), horizontale (H) et diagonales et dans la direction spectrale (K), s'exprimant comme suit : Xp+1,q,k Xp,q,k ,p 1, 4X)p,q,k = ,p=Puq=1, IXpa,k+1-Xp,q,k ,k B13227 - DD14758ST 14 Un avantage du système et du procédé d'acquisition proposés est qu'ils permettent d'acquérir des images hyperspectrales présentant une résolution spatiale supérieure à celle du capteur et/ou présentant un nombre de bandes spectrales 5 distinctes supérieur au nombre d'acquisitions effectuées par le capteur, de façon particulièrement simple. En particulier, l'acquisition d'une image hyperspectrale présentant une résolution supérieure à celle du capteur peut être effectuée par un système ne comportant pas de prisme ni de dispositif d'obturation 10 matriciel à ouverture codée. Un autre avantage majeur réside dans le fait que, lors de l'acquisition, aucun filtrage spectral n'est effectué, ce qui implique que le choix des bandes reconstruites ne dépend pas de l'acquisition mais résulte uniquement d'un choix effectué à la 15 reconstruction. Il est ainsi possible, à partir d'une même série de M images ym d'une scène acquises par le capteur 101, de choisir le nombre et la position des bandes reconstruites. Un compromis est bien entendu souvent à faire entre la taille des bandes spectrales, leur nombre, la résolution spatiale de l'image 20 reconstruite, et la qualité/fiabilité de l'image reconstruite. Quelle que soit la méthode choisie pour reconstruire l'image hyperspectrale x à partir des M images ym acquises par le capteur 101, il est nécessaire de connaître les différentes PSF hm,u de l'élément 105. Les PSF hm,u peuvent par exemple être 25 déterminées lors d'une phase de calibration du système, ou être déterminées par simulation à partir de modèles de réponse théorique de l'élément 105, ou encore être déterminées par une méthode combinant calibration et simulation. La figure 3 illustre de façon schématique et partielle 30 le système d'acquisition d'images hyperspectrales de la figure 1 pendant une phase de calibration permettant de déterminer les différentes PSF hm,u de l'élément 105. Dans cet exemple, pour la calibration, l'objet 103 de la figure 1 est remplacé par une source lumineuse ponctuelle 303 35 réglable en longueur d'onde, permettant de balayer successivement B13227 - DD14758ST 15 les K bandes spectrales Xk des images hyperspectrales que l'on souhaite acquérir. La source 303 est par exemple une source laser réglable en longueur d'onde. Alternativement, la source 303 peut comprendre une source large bande couplée à une pluralité de filtres optiques passe-bande activables individuellement correspondant aux différentes bandes spectrales des images hyperspectrales que l'on souhaite acquérir. Lors d'une phase de calibration, on peut prévoir, dans un premier temps, d'appliquer à l'élément 105 la commande C1 et de commander la source 303 pour émettre uniquement dans la bande spectrale Xi. Une image peut alors être acquise par le capteur 101, définissant la PSF hi,x1 de l'élément 105 pour la commande C1 et pour la longueur d'onde Xi. L'opération peut être répétée K fois sans modifier la commande de l'élément 105 mais en modifiant à chaque fois la longueur d'onde d'émission de la source, de façon à balayer les K bandes spectrales des images hyperspectrales que l'on souhaite acquérir. L'ensemble de l'opération peut alors être répété M fois en modifiant à chaque fois la commande appliquée à l'élément 105 de façon à balayer les M commandes Cm susceptibles d'être appliquées à l'élément 105 et obtenir ainsi les M x K PSF hm,xk de l'élément 105. Cette méthode permet d'obtenir des PSF hm,xk présentant un niveau de granularité ou pas de discrétisation identique à celui du capteur 101, et présentant un niveau de bruit correspondant à celui du capteur. Toutefois, il peut être avantageux de définir les PSF hm,xk avec un pas de discrétisation inférieur à celui du capteur et/ou avec un niveau de bruit inférieur à celui de capteur (notamment si on souhaite reconstruire une image hyperspectrale présentant une résolution spatiale supérieure à celle du capteur ou présentant un nombre d'images élémentaires xxk supérieur au nombre d'images ym acquises par le capteur). Pour cela, les PSF peuvent être déterminées par simulation, à l'aide d'un modèle comportemental théorique de l'élément 105. B13227 - DD14758ST 16 Une autre méthode particulièrement avantageuse de détermination des PSF hm,u, combinant calibration et simulation, va maintenant être décrite. Dans un premier temps, MxK images hsm,xk correspondant 5 aux MxK PSF de l'élément 105 à la résolution du capteur 101, peuvent être acquises en utilisant une source réglable en longueur d'onde selon la méthode de calibration susmentionnée. On prévoit ensuite de mutualiser les informations collectées lors de cette phase d'acquisition pour réduire le bruit et éventuellement 10 augmenter la résolution des PSF acquises. Pour chaque bande de longueurs d'ondes Xk des images hyperspectrales que l'on souhaite acquérir, un modèle comportemental théorique hu(1-2) peut être défini, où S2 désigne un jeu d'un nombre restreint de paramètres de l'élément 105. Le 15 modèle théorique hxk(1-2) peut par exemple présenter un pas de discrétisation inférieur à celui du capteur 101. Pour chaque commande Cm de l'élément 105, on fait correspondre les mesures effectuées hsm,xk avec le modèle comportemental théorique hxk de façon à déterminer un jeu de 20 paramètres S.2m définissant le comportement réel de l'élément 105 pour la commande Cm. Le jeu S.2m peut être déterminé à l'aide d'un algorithme de minimisation d'erreur, par exemple un algorithme itératif, résolvant le problème de minimisation suivant (dans le cas où l'on souhaite diminuer l'effet d'un bruit de type 25 gaussien): ç)nl = ar gminn(C (,(1) + YII12Sm,k - D shAk(S1)1122), k =1 où Ds est un opérateur de décimation correspondant au moyennage effectué par le capteur 101 lors de l'acquisition, par rapport à une image présentant la même résolution que le modèle théorique 30 hxk, où y est un paramètre scalaire de régularisation, et où C(S2) est un opérateur de contraintes sur un jeu de paramètres S). Chacune des MxK PSF hm,xk de l'élément 105 peut alors être déterminée par la formule suivante: B13227 - DD14758ST 17 hrn,,Ak - hAk(nrii) - Ceci permet de prendre en compte d'éventuelles imperfections du système d'acquisition non modélisables à l'aide d'a priori théoriques. Avantageusement, une fonction I d'interpolation et de lissage peut être appliquée. Les PSF hm,xk de l'élément 105 peuvent alors être déterminées par la formule suivante: = hAk(nm) - (hsm,Ak - D shAk (flm)) - Des modes de réalisation particuliers ont été décrits. 10 Diverses variantes et modifications apparaîtront à l'homme de l'art. En particulier, une nouvelle solution permettant d'acquérir des images hyperspectrales a été proposée. A partir de l'enseignement décrit ci-dessus, l'homme de l'art saura adapter 15 certains aspects de la méthode et du système décrits en fonction de l'application visée. Par exemple, bien que l'on ait représenté sur les figures uniquement des exemples de réalisation dans lesquels l'élément diffuseur et dispersif 105 fonctionne en transmission, les modes de réalisation décrits ne se limitent pas 20 à ce cas particulier. Ainsi, on pourra prévoir un système d'acquisition du type décrit en relation avec les figures 1 à 3, dans lequel l'élément 105 fonctionne en réflexion. REVENDICATIONS1. Système d'acquisition d'une image hyperspectrale (x), comportant : un capteur (101) d'images en niveaux de gris ; et un élément (105) diffuseur et dispersif placé sur le 5 trajet optique entre le capteur (101) et une scène (103), cet élément (105) comprenant une matrice de cellules à cristaux liquides commandables individuellement, dans lequel chaque cellule peut recevoir une tension de commande choisie parmi une série d'au moins trois tensions de commande distinctes. 10 2. Système selon la revendication 1, comportant en outre un module de commande adapté à commander le capteur (101) et l'élément (105) pour acquérir successivement un nombre entier M supérieur à 1 d'images élémentaires (ym) en niveaux de gris de la scène (103), en appliquant, lors de chaque acquisition, une 15 commande (Cm) spécifique à l'élément (105). 3. Système selon la revendication 2, dans lequel, pour chaque commande (Cm) appliquée à l'élément (105), l'élément (105) présente une fonction d'étalement de point qui dépend de la longueur d'onde des rayons qu'il reçoit. 20 4. Système selon la revendication 3, dans lequel les M commandes appliquées à l'élément (105) diffuseur et dispersif sont telles que les MxK fonctions d'étalement de point (hm,xk) de l'élément (105), correspondant aux M commandes appliquées et à K bandes spectrales de l'image hyperspectrale (x), avec K entier 25 supérieur à 1, soient irrégulières et soient toutes différentes les unes des autres. 5. Système selon l'une quelconque des revendications 2 à 3, comportant en outre un module de traitement adapté à reconstruire l'image hyperspectrale (x) à partir des M images 30 élémentaires (ym) acquises par le capteur (101). 6. Procédé de commande d'un système selon l'une quelconque des revendications 1 à 5, comportant une phase d'acquisition au cours de laquelle le capteur (101) et l'élément (105) sont commandés pour acquérir successivement un nombre entierB13227 - DD14758ST 19 M supérieur à 1 d'images élémentaires (ym) en niveaux de gris de la scène (103), en appliquant, lors de chaque acquisition, une commande (Cm) spécifique à l'élément (105). 7. Procédé selon la revendication 6, comportant en 5 outre une phase de reconstruction de l'image hyperspectrale (x) à partir des M images élémentaires (ym) acquises pendant la phase d'acquisition. 8. Procédé selon la revendication 6 ou 7, comportant en outre une phase préalable de calibration au cours de laquelle les 10 MxK fonctions d'étalement de point (hm,u) de l'élément (105), correspondant aux M commandes appliquées et à K bandes spectrales de l'image hyperspectrale (x), avec K entier supérieur à 1, sont déterminées. 9. Procédé selon la revendication 8, dans lequel la 15 phase de calibration comprend l'acquisition successivement par le capteur (101) de MxK images de taches résultant de la diffusion, par l'élément (105), pour les M commandes (Cm) de l'élément (105) et pour les K bandes spectrales lumineuses de l'image hyperspectrale, d'une source lumineuse (303) ponctuelle réglable 20 en longueur d'onde. 10. Procédé selon la revendication 9, dans lequel les MxK images acquises lors de la phase de calibration sont mises en correspondance avec un modèle comportemental théorique de l'élément (105).