DOMAINE TECHNIQUE

Le domaine technique de l'invention est lié à l'observation d'un échantillon, en particulier un échantillon biologique, par imagerie sans lentille, en mettant en œuvre un algorithme de reconstruction holographique aux performances améliorées.

ART ANTERIEUR

L'observation d'échantillons, et en particulier des échantillons biologiques, par imagerie sans lentille connaît un développement important depuis ces dix dernières années. Cette technique permet d'observer un échantillon en le disposant entre une source de lumière et un photodétecteur matriciel, sans disposer de lentille de grossissement optique entre l'échantillon et le photodétecteur. Ainsi, le photodétecteur collecte une image de l'onde lumineuse transmise par l'échantillon.

Cette image est formée de figures d'interférence entre l'onde lumineuse émise par la source et transmise par l'échantillon, et des ondes de diffraction, résultant de la diffraction par l'échantillon de l'onde lumineuse émise par la source. Ces figures d'interférences sont parfois dénommées figures de diffraction, ou désignées par le terme anglais « diffraction pattern ».

Le document WO2008090330 décrit un dispositif permettant l'observation d'échantillons biologiques, en l'occurrence des cellules, par imagerie sans lentille. Le dispositif permet d'associer, à chaque cellule, une figure d'interférence dont la morphologie permet d'identifier le type de cellule. L'imagerie sans lentille apparaît alors comme une alternative simple, et peu onéreuse, à un microscope classique. De plus, son champ d'observation est nettement plus important que ne peut l'être celui d'un microscope. On comprend alors que les perspectives d'application liées à cette technologie sont importantes. Ce document mentionne également la possibilité d'illuminer l'échantillon à l'aide de sources de lumière de différentes longueurs d'onde.

Le document US2012/0218379, postérieur au document précédent, reprend l'essentiel des enseignements de WO2008090330, tout en mentionnant la possibilité d'utiliser un photodétecteur matriciel couleur, mais l'information couleur est ensuite traitée pour former une image monochrome.

D'une façon générale, l'image formée sur le photodétecteur matriciel, comportant des figures d'interférence, peut être traitée par un algorithme de propagation numérique, de manière à estimer des propriétés optiques de l'échantillon. De tels algorithmes sont bien connus dans le domaine de la reconstruction holographique. Pour cela, la distance entre l'échantillon et le photodétecteur étant connue, on applique un algorithme de propagation, prenant en compte cette distance, ainsi que la longueur d'onde. On peut alors reconstituer une image d'une propriété optique de l'échantillon. Un algorithme de reconstruction numérique est par exemple décrit dans US2012/0218379.

Il est également connu que de tels algorithmes peuvent générer une image affectée d'un bruit de fond important, qualifié par le terme anglais « Twin Image », signifiant « image jumelle ». Un tel bruit est dû au fait que l'image de départ, formée par le photodétecteur, ne contient qu'une information partielle quant à l'onde lumineuse collectée par le photodétecteur. En effet, une image ne donne qu'une information relative à la partie réelle de l'onde lumineuse, obtenue par l'intensité mesurée. Mais une telle image ne comporte pas d'information relative à la partie imaginaire de l'onde lumineuse à laquelle est exposé le détecteur, en particulier sa phase. L'algorithme de reconstruction utilise donc une information incomplète, ce qui se traduit par l'apparition d'un bruit sur l'image reconstruite.

Un tel bruit de fond peut compliquer l'interprétation des images reconstruites par propagation numérique ; il est important d'en diminuer l'influence par des algorithmes appropriés.

Pour ce faire, la publication « Lensless phase contrast microscopy based on multiwavelenght Fresnel diffraction », Optics Letters Vol. 39, N°2, 15 janvier 2014, décrit un algorithme permettant d'améliorer la qualité des images reconstruites. Cette publication décrit un dispositif d'imagerie sans lentille basé sur l'utilisation de trois sources de lumière de longueurs d'onde différentes (respectivement 685 nm, 785 nm, 940 nm). L'échantillon est illuminé successivement par ces trois sources de lumière. Le photodétecteur acquiert alors autant d'images que de sources de lumière, ces images étant formées dans un plan, dit plan du détecteur, selon lequel s'étend le capteur. A chaque image correspond une longueur d'onde.

Une première image, d'une première longueur d'onde, est rétropropagée, selon ladite première longueur d'onde, dans un plan selon lequel s'étend l'objet, dit plan de l'objet, de façon à obtenir, dans ce plan objet, un premier champ complexe. La phase de ce premier champ complexe, dans le plan objet, est multipliée par un ratio entre la première longueur d'onde et une deuxième longueur d'onde. Ce champ complexe est alors propagé, selon ladite deuxième longueur d'onde, du plan objet vers le plan du détecteur, suite à quoi son module est remplacé par le module de l'image acquise à ladite deuxième longueur d'onde. Il est alors rétropropagé dans le plan de l'objet, pour une deuxième itération. Le procédé itératif se poursuit jusqu'à une l'atteinte d'un critère de convergence. Le document WO2014035238, dont certains inventeurs sont les auteurs de la publication évoquée ci-dessus, comporte les mêmes enseignements.

La publication Bao P,« Optical surface profile measurement using phase retrieval by tuning the illumination wavelenght », Optics Communications 285, 5029-5036, 2012 décrit un algorithme itératif permettant une reconstruction d'une surface tridimensionnelle d'un objet transparent à partir d'une image acquise par un capteur d'image dans une configuration sans lentille. Cet algorithme vise à illuminer un objet à l'aide de deux longueurs d'onde. L'image acquise par le capteur d'image, à chaque longueur d'onde, est propagée dans un plan de l'objet. On estime ensuite une différence de phases, à chacune de ces longueurs d'onde, dans le plan objet. La surface tridimensionnelle est reconstruite à partir de la différence de phase ainsi estimée.

Le document WO2015/015023 décrit un algorithme itératif de reconstruction holographique basé sur un capteur d'image couleur dans une configuration sans lentille. Selon cet algorithme, l'échantillon est illuminé par une source de lumière polychromatique. L'image acquise est décomposée en différentes composantes spectrales. Chaque composante est rétropropagée dans un plan selon lequel s'étend l'échantillon, la propagation des composantes étant réalisée selon différentes distances les unes des autres. Ce document décrit ensuite un mélange des spectres de chaque hologramme dans le domaine de Fourier, permettant de générer une image unique, cette dernière étant utilisée dans une boucle d'itération suivante.

Les inventeurs proposent une méthode alternative à celles proposées dans les publications précédentes, permettant une estimation fiable des propriétés optiques d'un échantillon.

EXPOSE DE L'INVENTION

Un premier objet de l'invention est un procédé d'observation d'un échantillon, comportant les étapes suivantes :

i) illumination dudit échantillon à l'aide d'une source de lumière, apte à produire une onde lumineuse se propageant selon un axe de propagation.

ii) acquisition, à l'aide d'un photodétecteur, d'une pluralité d'images de l'échantillon, formées dans un plan de détection, l'échantillon étant disposé entre la source de lumière et le photodétecteur, chaque image étant représentative d'une onde lumineuse, transmise par l'échantillon, sous l'effet de ladite illumination, chaque image étant acquise selon une bande spectrale différente l'une de l'autre,

le procédé étant caractérisé en ce qu'il comprend également les étapes suivantes : iii) détermination, à partir de chaque image acquise selon chaque bande spectrale, d'une amplitude complexe initiale de l'onde lumineuse transmise, selon ladite bande spectrale, dans ledit plan de détection,

iv) rétropropagation de chaque amplitude complexe établie dans le plan de détection, selon une bande spectrale, pour déterminer une amplitude complexe de l'onde transmise, selon chaque bande spectrale, dans un plan selon lequel s'étend l'échantillon,

v) combinaison d'une pluralité d'amplitudes complexes déterminées lors de l'étape iv), à différentes bandes spectrales, pour calculer une fonction de pondération dans le plan de l'échantillon,

vi) projection de ladite fonction de pondération dans le plan de détection, de façon à obtenir, pour chaque bande spectrale, une fonction de pondération dans ledit plan de détection,

vii) mise à jour de chaque amplitude complexe de l'onde lumineuse transmise, selon chaque bande spectrale, dans le plan de détection, à l'aide de ladite fonction de pondération déterminée, dans ladite bande spectrale, au cours de l'étape vi),

viii) répétition des étapes iv à vii jusqu'à l'atteinte d'un critère d'arrêt.

Ainsi chaque itération comprend une propagation, dite rétropropagation, du plan du détecteur vers le plan de l'échantillon, d'une pluralité d'amplitudes complexes, dans différentes bandes spectrales. Ces différentes amplitudes complexes sont combinées, dans le plan de l'échantillon, pour former une fonction de pondération. Cette combinaison, dans le plan de l'échantillon, d'une pluralité d'amplitudes complexes, correspondant à différentes plages spectrales, a pour effet de lisser le bruit affectant chacune d'entre elles, ce bruit étant la conséquence de la rétropropagation.

La fonction de pondération est ensuite propagée du plan de l'échantillon vers le plan de détection, où elle est utilisée pour former une nouvelle estimation de l'amplitude complexe de l'onde lumineuse à laquelle est exposé l'échantillon, et cela dans chaque bande spectrale considérée.

De façon alternative aux méthodes exposées dans l'art antérieur, la fonction de pondération, dans ledit plan de l'échantillon, peut être calculée en calculant une somme pondérée de différentes amplitudes complexes, ou de leurs logarithmes, dans le plan de l'échantillon, de l'onde lumineuse transmise, lesdites amplitudes complexes étant respectivement associées à différentes bandes spectrales. La fonction de pondération, dans ledit plan de l'échantillon, peut également être déterminée en calculant une somme pondérée du module et/ou de l'argument de différentes amplitudes complexes, dans le plan de l'échantillon, de l'onde lumineuse transmise, lesdites amplitudes complexes étant respectivement associées à différentes bandes spectrales.

Par onde lumineuse transmise, on entend l'onde lumineuse à laquelle est exposé le photodétecteur.

Le procédé peut comprendre l'une quelconque des caractéristiques suivantes, prises isolément ou selon les combinaisons techniquement réalisables :

lors de l'étape iii), le module de l'amplitude complexe de l'onde transmise dans une bande spectrale est déterminé en normalisant l'intensité de l'image mesurée par le photodétecteur 16, dans ladite bande spectrale, par une intensité de référence mesurée par ledit photodétecteur en l'absence d'échantillon.

Lors de l'étape iv), la détermination de l'amplitude complexe dans un plan de l'échantillon, dans une bande spectrale, est obtenue en appliquant un opérateur de propagation, selon un axe de propagation, à ladite amplitude complexe, définie dans la même bande spectrale, dans le plan de détection.

Lors de l'étape vi), ladite fonction de pondération, dans le plan de détection, est obtenue en appliquant un opérateur de propagation, selon l'axe de propagation, à la fonction de pondération, déterminée, dans le plan de l'échantillon, au cours de l'étape v).

Lors de l'étape vii), le module de l'amplitude complexe de l'onde lumineuse transmise, selon une bande spectrale, dans le plan de détection, est calculé en fonction du module de ladite amplitude complexe initiale, dans ladite bande spectrale.

Lors de l'étape vii), l'argument de l'amplitude complexe de l'onde lumineuse transmise, selon une bande spectrale, dans le plan de détection, est calculé en fonction de l'argument de la fonction de pondération déterminée, dans ledit plan de détection et à ladite bande spectrale, au cours de l'étape vi).

Lors de l'étape v), ladite fonction de pondération peut être commune à l'ensemble des bandes spectrales. Alternativement, cette étape v) peut comprendre la détermination d'une pluralité de fonctions de pondérations, chaque fonction de pondération étant associée à une bande spectrale.

Le procédé peut comporter, suite à l'étape viii), l'étape ix) : formation d'une image représentative du module ou de l'argument de l'amplitude complexe de l'onde, transmise par l'échantillon, dans le plan de l'échantillon ou dans le plan de détection, selon au moins une bande spectrale (λ,).

Un autre objet de l'invention est un dispositif d'observation d'un échantillon comportant : une source de lumière, apte à illuminer ledit échantillon,

un photodétecteur, l'échantillon étant disposé entre la source de lumière et le photodétecteur,

le photodétecteur étant apte à former une pluralité d'images, dans un plan de détection, d'une onde lumineuse transmise par l'échantillon sous l'effet de l'illumination par ladite source de lumière, chaque image étant obtenue selon une bande spectrale différente l'une de l'autre,

un processeur, apte à traiter ladite pluralité d'images en exécutant des instructions, programmées dans une mémoire, mettant en œuvre un procédé tel que précédemment décrit.

FIGURES

La figure 1 représente un premier exemple de dispositif pour la mise en œuvre de l'invention, l'échantillon analysé étant une lame d'anatomopathologie.

La figure 2 représente un premier exemple de dispositif pour la mise en œuvre de l'invention, l'échantillon analysé étant un liquide corporel comportant des particules.

La figure 3 représente le plan de détection, sur lequel une image est formée, ainsi que le plan de l'échantillon. Cette figure illustre également les liens entre les principales grandeurs mises en œuvre dans les différents modes de réalisation exposés.

La figure 4 représente un logigramme décrivant l'enchaînement des principales étapes d'un procédé itératif de reconstruction.

La figure 5 représente un deuxième exemple de dispositif pour la mise en œuvre de l'invention, l'échantillon analysé étant une lame d'anatomopathologie.

Les figures 6A, 6B et 6C représentent des images reconstruites, dans un plan de l'échantillon, dans une première bande spectrale, selon un algorithme de reconstruction itératif, ces images étant respectivement obtenues suite à un nombre d'itérations respectivement égal à 1, 3 et 10. Les figures 7 A, 7B et 7C représentent des images reconstruites, dans un plan de l'échantillon, dans une deuxième bande spectrale, selon un algorithme de reconstruction itératif, ces images étant respectivement obtenues suite à un nombre d'itérations respectivement égal à 1, 3 et 10. Les figures 8A, 8B et 8C représentent des images reconstruites, dans un plan de l'échantillon, dans une troisième bande spectrale, selon un algorithme de reconstruction itératif, ces images étant respectivement obtenues suite à un nombre d'itérations respectivement égal à 1, 3 et 10. Les figures 9A, 9B et 9C représentent des images composites, obtenues par combinaison des images reconstruites dans la première, la deuxième et la troisième bande spectrales, ces images correspondant respectivement à 1, 3 et 10 itérations. Ces images sont des images en couleur, représentées ici en noir et blanc.

EXPOSE DE MODES DE REALISATION PARTICULIERS

La figure 1 représente un exemple de dispositif objet de l'invention. Une source de lumière 11 est apte à produire une onde lumineuse 12, dite onde lumineuse incidente, en direction d'un échantillon 10, selon un axe de propagation Z.

L'échantillon 10 peut être un échantillon biologique que l'on souhaite caractériser. Il peut par exemple s'agir d'une lame de tissu, ou lame d'anatomopathologie, comportant une fine épaisseur de tissu déposée sur une lame transparente 15. Par fine épaisseur, on entend une épaisseur de préférence inférieure à 100 μιη, et de préférence inférieure à 10 μιη, typiquement quelques micromètres. Un tel échantillon est représenté sur la figure 1. On observe que l'échantillon s'étend selon un plan P _¾ dit plan de l'échantillon, perpendiculaire à l'axe de propagation Z.

L'échantillon 10 peut également comporter un milieu 14, solide ou liquide, comportant des particules 1, 2, 3, 4, 5 à caractériser, un tel cas de figure étant représenté sur la figure 2. Il peut par exemple s'agir de particules biologiques dans un milieu de culture, ou dans un liquide corporel. Par particule biologique, on entend une cellule, une bactérie ou autre microorganisme, un champignon, une spore... Le terme particules peut également désigner des microbilles, par exemple des microbilles métalliques, des microbilles de verre ou des microbilles organiques, couramment mises en œuvre dans des protocoles biologiques. Il peut également s'agir de gouttelettes insolubles baignant dans un milieu liquide, par exemple des gouttelettes lipidiques dans une émulsion de type huile dans eau. Ainsi, le terme particule désigne à la fois des particules endogènes, initialement présentes dans l'échantillon examiné, et des particules exogènes, ajoutées à cet échantillon avant l'analyse. D'une façon générale, une particule a une taille avantageusement inférieure à 1 mm, voire inférieure à 500 μιη, et de préférence une taille comprise entre 0.5 μιη et 500 μιη.

La distance Δ entre la source de lumière et l'échantillon est de préférence supérieure à 1 cm. Elle est de préférence comprise entre 2 et 30 cm. De préférence, la source de lumière, vue par l'échantillon, est considérée comme ponctuelle. Cela signifie que son diamètre (ou sa diagonale) est préférentiellement inférieur au dixième, mieux au centième de la distance entre l'échantillon et la source de lumière. Ainsi, de préférence, la lumière parvient à l'échantillon sous la forme d'ondes planes, ou pouvant être considérées comme telles. La source de lumière 11 est apte à produire une pluralité d'ondes lumineuses incidentes

12i 12 _n , chaque i ^eme onde lumineuse 12, s' étendant selon une i ^eme bande spectrale λ,. Les bandes spectrales 12i...l2 _n sont différentes les unes des autres, et, de préférence, ne se recouvrent pas. Dans l'exemple de dispositif représenté sur les figures 1 et 2, la source de lumière comporte trois sources de lumières élémentaires de type diodes électroluminescentes (LED, acronyme du terme anglais Light Emitting Diode), lli, II2 et II3 émettant respectivement dans les bandes spectrales λι = 450nm - 465 nm ; λ2 = 520nm - 535 nm ; λ3 = 620nm - 630 nm. De préférence, il n'y a pas de recouvrement entre les différentes bandes spectrales ; un recouvrement négligeable, par exemple concernant moins de 25%, mieux moins de 10% de l'intensité lumineuse émise, est toutefois envisageable. Dans ces exemple, la source de lumière 11 comporte une diode LED multiple de marque Crée (marque déposée) et de référence XLamp (marque déposée) MC-E. Cette diode comporte quatre diodes électroluminescentes élémentaires adressables individuellement, dont trois seulement sont mises en œuvre dans le cadre de cette invention, la quatrième étant une LED blanche. Les sources de lumière élémentaires peuvent être des sources temporellement cohérentes, comme des diodes Laser. D'autres configurations de sources lumineuses 11 sont possibles et décrites par la suite.

La source de lumière 11 est de préférence ponctuelle. Elle peut notamment comprendre un diaphragme 18, ou filtre spatial. L'ouverture du diaphragme est typiquement comprise entre 5 μιη et 1mm, de préférence entre 50 μιη et 500 μιη, par exemple 150 μιη. Le diaphragme peut être remplacé par une fibre optique, dont une première extrémité est placée face à une source de lumière élémentaire lli, II2 ou II3, et dont une deuxième extrémité est placée face à l'échantillon.

La source de lumière 11 comporte de préférence un diffuseur 17, disposé entre chaque source de lumière élémentaire lli, II2 et II3 et le diaphragme 18. Les inventeurs ont constaté que l'usage d'un tel diffuseur permet de s'affranchir de contraintes de centrage de chaque source élémentaire par rapport à l'ouverture du diaphragme. Autrement dit, l'utilisation d'un tel diffuseur permet d'utiliser une source de lumière élémentaire 11,, avec 1 < i < 3, légèrement décentrée par rapport à l'ouverture du diaphragme 18. Dans cet exemple, le diaphragme est fourni par Thorlabs sous la référence P150S.

De préférence, chaque source de lumière élémentaire 11, est de faible largeur spectrale, par exemple inférieure à 100 nm, voire 20 nm. Le terme largeur spectrale désigne la largeur à mi- hauteur de la bande d'émission de la source de lumière considérée. Dans cet exemple, le diffuseur mis en œuvre est un diffuseur 40° (Référence Light Shaping Diffuser 40°, fabriqué par Luminit). La fonction d'un tel diffuseur est de répartir le faisceau lumineux, produit par une source de lumière élémentaire 11,, selon un cône d'angle a, a étant égal à 40° dans le cas présent. De préférence, l'angle de diffusion a varie entre 10° et 60°. L'échantillon 10 est disposé entre la source de lumière 11 et un photodétecteur matriciel 16. Ce dernier s'étend de préférence parallèlement, ou sensiblement parallèlement à la lame transparente 15 supportant l'échantillon.

Le terme sensiblement parallèlement signifie que les deux éléments peuvent ne pas être rigoureusement parallèles, une tolérance angulaire de quelques degrés, inférieure à 10° étant admise.

Le photodétecteur 16 est un imageur, apte à former une image selon un plan de détection P. Dans l'exemple représenté, il s'agit d'un photodétecteur matriciel comportant une matrice de pixels, de type CCD ou un CMOS. Les CMOS sont les photodétecteurs préférés, car la taille des pixels est plus faible, ce qui permet d'acquérir des images dont la résolution spatiale est plus favorable. Dans cet exemple, le détecteur est un capteur CMOS commercialisé par Omnivision sous la référence OV5647. Il s'agit d'un capteur CMOS RGB, comprenant 2592 x 1944 pixels, avec un pas inter-pixels est de 1.4 μιη. La surface utile du photodétecteur est de 3.6 x 2.7 mm ² . Le plan de détection P s'étend de préférence perpendiculairement à l'axe de propagation Z de l'onde lumineuse incidente 12. De préférence, le photodétecteur comprend une matrice de pixels, au-dessus de laquelle est disposée une fenêtre de protection transparente. La distance entre la matrice de pixels et la fenêtre de protection est généralement comprise entre quelques dizaines de μιη à 150 ou 200 μιη. Les photodétecteurs dont le pas inter pixel est inférieur à 3 μιιι sont préférés, afin d'améliorer la résolution spatiale de l'image. Le photodétecteur peut comprendre un système de renvoi d'image de type miroir vers une matrice de pixels, auquel cas le plan de détection correspond au plan selon lequel s'étend le système de renvoi d'image. D'une façon générale, le plan de détection P correspond au plan dans lequel une image est formée.

La distance d entre l'échantillon 10 et la matrice de pixels du photodétecteur 16 est, dans cet exemple, égale à 300 μιη. D'une manière générale, et cela quel que soit le mode de réalisation, la distance d entre l'échantillon et les pixels du photodétecteur est préférentiellement comprise entre 50 μιη et 2 cm, de préférence comprise entre 100 μιη et 2 mm. On remarque l'absence d'optique de grossissement entre le photodétecteur 16 et l'échantillon 10. Cela n'empêche pas la présence éventuelle de microlentilles de focalisation au niveau de chaque pixel du photodétecteur 16, ces dernières n'ayant pas de fonction de grandissement de l'image. La figure 3 représente un échantillon 10, comportant des objets diffractants 32 disposés autour de zones non diffractantes ou peu diffractantes 31, qualifiées par la suite de zones pauvres. L'échantillon peut être solide, par exemple dans le cas d'un tissu déposé sur une lame d'anatomopathologie. Il peut également être liquide, par exemple dans le cas d'un liquide corporel ou d'un milieu de culture cellulaire. Le photodétecteur 16 est apte à réaliser une image I, d'une onde lumineuse 22, transmise par l'échantillon 10 lorsque ce dernier est éclairé par une onde incidente 12,, dans la i ^eme bande spectrale λ,. La bande spectrale de l'onde lumineuse transmise 22, comporte tout ou partie de la bande spectrale de l'onde incidente 12,. L'onde lumineuse 22,, transmise par l'échantillon dans la bande spectrale λ,, résulte de l'interaction de l'échantillon 10 avec l'onde lumineuse incidente 12i, produite par la source de lumière élémentaire 11,.

Sous l'effet de l'onde lumineuse incidente 12,, l'échantillon 10 peut engendrer une onde diffractée, susceptible de produire, au niveau du plan de détection P, des interférences, en particulier avec une partie de l'onde lumineuse incidente 12, transmise par l'échantillon. Ces interférences donnent lieu, sur l'image acquise par le photodétecteur, à une pluralité de figures de diffraction élémentaires, chaque figure de diffraction élémentaire 36 comportant une zone centrale et plusieurs anneaux de diffraction concentriques. Chaque figure de diffraction élémentaire 36 est due à un objet diffractant 32 dans l'échantillon.

Par ailleurs, l'échantillon peut absorber une partie de l'onde lumineuse incidente 12,. Ainsi, l'onde lumineuse 22,, dans une bande spectrale λ,, transmise par l'échantillon, et à laquelle est exposé le photodétecteur matriciel 16, peut comprendre :

une composante résultant de la diffraction, précédemment décrite, cette composante de diffraction pouvant notamment se traduire par la présence de figures de diffraction élémentaires sur le photodétecteur 16, chaque figure de diffraction élémentaire pouvant être associée à un élément diffractant 32 de l'échantillon. Un tel élément diffractant peut être une cellule, ou une particule ou tout autre objet diffractant 32 présent dans l'échantillon 10.

une composante résultant de l'absorption de l'onde lumineuse incidente 12, dans l'échantillon. Un processeur 20, par exemple un microprocesseur, est apte à traiter chaque image générée par le photodétecteur matriciel 16. En particulier, le processeur est un microprocesseur relié à une mémoire programmable 23 dans laquelle est stockée une séquence d'instructions pour effectuer les opérations de traitement d'images et de calculs décrites dans cette description. Il peut également être relié à un écran d'affichage 24.

On décrit ci-après les étapes d'un procédé itératif d'obtention d'une image de l'échantillon 10, en lien avec les figures 3 et 4. l ^ere étape : initialisation

Lors d'une première étape 100 d'acquisition d'images, chaque source de lumière élémentaire 11, de la source de lumière 11 est successivement activée, chaque source de lumière produisant une onde lumineuse incidente (12i, .. 12 _N ), dans une bande spectrale (À _I ,...À _N ), selon un axe de propagation Z, en direction de l'échantillon 10.

A chaque acquisition, le photodétecteur matriciel recueille une image I, correspondant à une bande spectrale λ,, l'indice i, relatif à la bande spectrale, étant compris entre 1 et N, N étant le nombre de bandes spectrales considérées. Dans l'exemple représenté sur les figures 1 et 2, la source de lumière 11 comporte trois sources de lumière élémentaires lli, 11 ₂ et 11 ₃ . Le photodétecteur recueille trois images li, l ₂ , b, correspondant respectivement aux bandes spectrales λι, λ ₂ et λ3.

L'échantillon est placé à une coordonnée axiale z = 0, selon l'axe de propagation Z. On désigne par r une coordonnée radiale, c'est-à-dire une coordonnée dans un plan perpendiculaire à l'axe de propagation Z. Le plan z = d correspond au plan de détection, tandis que le plan z = 0 correspond à un plan passant par l'échantillon, dit plan de l'échantillon et noté Po.

Si 7 ^=d (r)= if (r) désigne la valeur de l'intensité collectée, dans la bande spectrale λ,, par le pixel du détecteur de coordonnée radiale r dans le plan de détection, il est possible d'établir, à l'aide de l'image I, une amplitude com (r) de l'onde 22, au niveau dudit pixel de coordonnée r, dont le module peu l'expression :

L'exposant d exprime le fait que l'amplitude complexe est déterminée selon le plan P de l'échantillon, d'équation z = d. L'amplitude complexe a (r) comporte un module et un argument, de telle sorte que :

af {r) = Mf où :

Mf (r) désigne le module de l'amplitude complexe de l'onde lumineuse détectée par le photodétecteur, dans la i ^eme bande spectrale λ,, à une coordonnée radiale r dans le plan de détection ;

- <pf (r) désigne la phase de l'amplitude complexe de l'onde lumineuse détectée par le photodétecteur, dans la i ^eme bande spectrale λ, et à ladite coordonnée radiale r dans plan de détection.

Or, le photodétecteur matriciel ne fournit aucune information relative à la phase de l'onde lumineuse. Aussi, lors de l'étape 100, on considère que e ⁷ ^^ est égal à une valeur initiale arbitraire, par exemple égale à 1.

L'amplitude complexe af (r) peut être exprimée de façon normalisée, selon l'expression :

Af (r) =

¹ 1 jmean

où :

jmean désigne l'intensité moyenne l'onde lumineuse 12, émise par la source de lumière 11 dans la i ^eme bande spectrale λ, ; cette intensité moyenne peut être déterminée expérimentalement, en disposant le photodétecteur 16 face à la source de lumière 11, sans échantillon interposé entre l'un et l'autre, et en calculant la moyenne des pixels de l'image acquise par le photodétecteur 16.

Af (r) désigne l'amplitude complexe normalisée de l'onde lumineuse 22, détectée par le photodétecteur matriciel 16 dans la i ^ème bande spectrale λ,.

La normalisation peut également être réalisée en divisant l'amplitude complexe a (r) par I™ ^ean (r), ce terme représentant l'intensité lumineuse, à la coordonnée radiale r, mesurée en l'absence d'échantillon.

L'amplitude complexe normalisée Af (r) comporte un module et un argument, de telle sorte que :

où :

mf (r) désigne le module de l'amplitude complexe normalisée Af (r) ;

<pf (r) désigne la phase de l'amplitude complexe normalisée, qui est également la phase de l'amplitude complexe œf (r). La première étape 100 permet d'attribuer à chaque amplitude complexe a (r) ou à chaque amplitude complexe normalisée Af (r) une valeur initiale à partir de l'image I, détectée par le photodétecteur dans la i ^eme bande spectrale λ,, de telle sorte que :

ou

L'indice p correspond au rang de l'itération du procédé itératif décrit par la suite. L'étape 100 étant une étape d'initialisation, on attribue la valeur 1 à cet indice.

En adressant tout ou partie des pixels r du photodétecteur 16, on obtient une image complexe, ou champ complexe, de l'onde lumineuse 22, au niveau du plan du détecteur, cette image rassemblant les amplitudes complexes af (r) ou les amplitudes complexes normalisées Af (r).

Dans la suite de la description, on ne considère que l'amplitude complexe normalisée Af (r), sachant que le raisonnement s'applique également à l'amplitude complexe a (r).

Cette première étape est répétée pour chaque bande spectrale (λι...λ _Ν ) détectée par le photodétecteur.

2' ^eme étape : rétropropagation dans le plan Pp de l'échantillon

Durant une deuxième étape 200, on estime l'amplitude complexe normalisée A° _{i p} (r ^' ) de l'onde 22, à laquelle est exposé le détecteur, dans le plan de l'échantillon P ₀ . Cette estimation est réalisée par rétropropagation de l'amplitude complexe normalisée Af _p (r), déterminée dans le plan de détection P, du plan de détection P vers le plan de l'échantillon P ₀ .

L'indice p désigne le rang de l'itération. Lors de la première itération (p = 1), on utilise l'amplitude complexe normalisée initiale (r) = Af (r) obtenue suite à la première étape 100. Lors des itérations suivantes (p > 1), on utilise l'amplitude complexe résultant de l'itération précédente, comme cela sera détaillé par la suite. Selon les principes bien connus de la reconstruction holographique numérique, en effectuant un produit de convolution entre l'amplitude complexe de l'onde lumineuse 22,, relative à la bande spectrale λ,, déterminée dans le plan de détection z = d, et un opérateur de propagation h(r,z), il est possible de reconstruire une amplitude complexe de cette même onde lumineuse en tout point de coordonnées (r,z) de l'espace, et en particulier dans le plan Po de l'échantillon. Autrement dit, l'amplitude complexe normalisée A ^z _{i p} j) de l'onde lumineuse 22, peut s'obtenir, à un point de coordonnées (r,z), à partir de Af _p ^d (r), selon l'opération :

A , _p (r) = A ; ^d (r) * h _i (r, z - d),

Où ή ί désigne l'opérateur de propagation dans la bande spectrale λ,.

Lorsque la reconstruction s'effectue selon le sens la propagation de la lumière, par exemple de l'échantillon vers le photodétecteur, on parle de propagation. Lorsque la reconstruction s'effectue selon le sens inverse de la propagation de la lumière, par exemple du photodétecteur vers l'échantillon, on parle de rétropropagation.

L'opérateur de propagation peut notamment reposer sur le modèle de diffraction de Fresnel. Dans cet exemple, l'opérateur de propagation est la fonction de Fresnel-Helmholtz, telle que :

1 f ^*2

h(r, z) = — e ^J exp(jn—).

où λ désigne la longueur d'onde.

Ainsi,

■¾ » = A{ _v (r) = Af;"(r) . h _M (r, -d)

où

r' désigne les coordonnées radiales dans le plan du photodétecteur (z = d), r désigne les coordonnées radiales dans le plan de reconstruction (z = 0),

λ, désigne la longueur d'onde centrale de la bande spectrale considérée.

A° _p (r) est donc obtenu par rétropropagation de Af _p (r) selon la distance d séparant le plan de détection P du plan de l'échantillon Po.

Cette deuxième étape est répétée pour chaque bande spectrale (ÀI...ÀN) émise par la source de lumière 11 ou, de façon plus générale, pour chaque bande spectrale (ÀI...ÀN) respectivement associée à chaque image (II... IN) détectée par le photodétecteur 16. II est possible, à ce stade, d'éta blir une image du module ou de la phase de l'amplitude A° _p (r) complexe de chaque onde lumineuse 22,, dans le plan de l'échantillon Po, l'amplitude complexe étant normalisée ou non, en calculant la valeur de A° _p (r) aux différentes coordonnées r dans le plan de l'échantillon.

Chaque image du module de l'amplitude complexe A° _p (r) est représentative de l'intensité de l'onde lumineuse au niveau de l'échantillon, tandis que chaque image de l'argument de l'amplitude complexe A° _p (r) est repréntative de la phase de l'intensité de l'onde lumineuse au niveau de l'échantillon.

Lorsque, comme dans le cas présent, on met en œuvre trois bandes spectrales centrées respectivement sur des longueurs d'onde dans le bleu, le vert et le rouge, l'information contenue dans les trois images permet d'obtenir une image en couleur de l'échantillon.

On note que l'amplitude complexe normalisée A° _p (r) est équivalente à une fonction de transmission de l'onde incidente 12, par l'échantillon 10 à la coordonnée radiale r.

3 ^ième étape : détermination de la fonction de pondération

Au cours de l'étape 300, on constitue, dans le plan de l'échantillon, une fonction de pondération, notée F _p (r), de l'amplitude complexe de l'onde lumineuse transmise par l'échantillon dans les différentes bandes spectrales λ, considérées.

Selon cet exemple, la fonction de pondération F _p °(r), dans le plan de l'échantillon, peut être commune à chaque bande spectrale. Elle est obtenue en combinant les amplitudes complexes normalisées A° _p (r) de l'onde lumineuse transmise par l'échantillon, dans le plan P ₀ de l'échantillon et dans les différentes bandes spectrales λ.

Selon un exemple, la fonction de pondération est obtenue par une somme pondérée de chaque amplitude complexe déterminée lors de l'étape 200, dans le plan Po de l'échantillon, selon l'expression :

F _p ° (r) = -L _M ° _p (r)

i

où ki désigne un facteur de pondération, positif, associé à la i ^ème bande spectrale λ,.

Les facteurs de pondération peuvent être égaux les uns aux autres, par exemple égaux à 1/3. D'autres façons de déterminer la fonction de pondération, dans le plan de l'échantillon, sont détaillées par la suite.

4' ^eme étape : propagation de la fonction de pondération dans le plan du détecteur

L'étape 400 vise à propager, du plan de l'échantillon Po vers le plan du détecteur P, la fonction de pondération F _p (r) déterminée, lors de l'étape précédente, dans le plan de l'échantillon Po.

L'opérateur de propagation étant dépendant de la longueur d'onde, cette propagation est effectuée pour chaque bande spectrale λ, considérée.

Ainsi, pour chaque bande spectrale λ,, F^ _p (r) = F _p (r) * h i (r, z = d). Lorsque l'opérateur de propagation est un opérateur de Fresnel-Helmholtz tel que précédemment défini,

L'opérateur de propagation dépendant de la longueur d'onde, on détermine autant de fonctions de pondération, dans le plan de détection, que de bandes spectrales considérées.

r' désigne les coordonnées radiales dans le plan de l'échantillon (z = 0),

r désigne les coordonnées radiales dans le plan de reconstruction, c'est-à-dire le plan du détecteur (z = d),

λ, désigne la longueur d'onde centrale de la bande spectrale considérée. 5' ^eme étape : mise à jour de l'amplitude complexe dans le plan du détecteur

Dans l'étape 500, on utilise la valeur de la fonction de pondération, dans le plan de détection z = d, pour mettre à jour l'estimation de l'amplitude complexe normalisée de l'onde lumineuse 22,, à laquelle est exposé le photodétecteur 16, dans la bande spectrale λ,.

La formule de mise à jour est :

où :

|F^ _p (r) | désigne le module de F^ _p (r) ;

mf (r) désigne le module de l'amplitude complexe initiale normalisée déterminé, à partir de l'image I,, lors de la première étape 100. Ce terme a une fonction d'attache aux données mesurées ;

<p _p désigne une estimation de la phase de l'amplitude complexe de l'onde 22, dans la i ^eme bande spectrale λ, ;

Afp Çr) désigne l'amplitude complexe de l'onde lumineuse 22i transmise par l'échantillon, dans le plan du photodétecteur 16, cette amplitude complexe constituant la base de l'itération suiva nte.

A la suite de cette éta pe, une nouvelle itération peut commencer, la donnée d'entrée de cette nouvelle itération p + 1 étant Af _{>p+ 1} r) = A _p r), cette nouvelle itération débutant par la rétropropagation de chaque amplitude complexe normalisée Af _p+1 r), pour les différentes bandes spectrales considérées, dans le plan Po de l'échantillon, selon l'étape 200. Les étapes 200 à 500 sont réalisées de manière itérative, soit selon un nom bre d'itérations _ma x prédéterminé, soit jusqu'à l'atteinte d'un critère de convergence, ce dernier pouvant être, par exemple, exprimé sous la forme d'un écart entre l'estimation de deux mêmes grandeurs entre deux itérations successives. Lorsque cet écart est inférieur à un seuil ε donné, le critère de convergence est atteint. Par exemple, le processus est arrêté lorsque l'une de ces conditions est atteinte :

Arg (Al _p {r) - A° _p+1 {r)) < ε ;

Cette liste n'est pas limitative.

A l'issue du procédé, on dispose d'une estimation de l'amplitude complexe de l'onde lumineuse 22i, transmise par l'échantillon, et à laquelle est exposé le photodétecteur, dans le plan du détecteur P, d'équation z = d et/ou dans le plan de l'échantillon Po, d'équation z = 0, et cela pour chaque bande spectrale considérée. En utilisant les différentes amplitudes complexes reconstruites A _p r) dans le plan de l'échantillon, on obtient une représentation précise de ce dernier, dans chacune des bandes spectrales considérées, en particulier en formant des images à partir du module ou de la phase desdites amplitudes complexes.

Comme précédemment évoqué, lorsque les bandes spectrales sont distribuées dans le spectre visible, les images du module ou de la phase peuvent être combinées, par exemple superposées de façon à obtenir des représentations en couleur. II est rappelé que cet algorithme, bien que décrit en relation avec une amplitude complexe normalisée A _t , s'applique également à l'amplitude complexe non normalisée a _t .

Apport de la fonction de pondération

Un des points importants de cet algorithme itératif est la constitution de la fonction de pondération F°(r) dans le plan de l'échantillon. En effet, d'une façon générale, la détermination de l'amplitude complexe d'une onde lumineuse à partir d'une image acquise par un photodétecteur est insuffisante, car l'information quant à la phase de l'onde n'est pas enregistrée par le photodétecteur, ce dernier n'étant sensible qu'à l'intensité, correspondant au module de l'amplitude complexe de l'onde.

Ainsi, comme indiqué lors de la description de l'étape 100, l'amplitude complexe f(r) ou l'amplitude complexe normalisée Af(r) déterminées au cours de cette étape ne comportent pas d'information quant à la phase de l'onde lumineuse qu'elles représentent. Ce manque d'information se traduit, lors de la rétropropagation du plan du détecteur P vers le plan de l'échantillon P ₀ , faisant l'objet de l'étape 200, par la formation d'artéfacts désignés par le terme « twin image ».

Les inventeurs ont constaté que ces artéfacts affectent principalement les zones pauvres 31 disposées au voisinage d'éléments diffractants 32, c'est-à-dire les zones situées entre deux éléments diffractants adjacents 32. En outre, ils ont constaté que ces artéfacts sont susceptibles de fluctuer en fonction de la longueur d'onde. Ainsi, en combinant les amplitudes complexes dans différentes longueurs d'onde, rétropropagées dans le plan de l'échantillon, les artéfacts, au niveau des zones pauvres 31, se moyennent statistiquement. Ce lissage statistique augmente alors le rapport signal sur bruit de l'image complexe rétropropagée dans le plan de l'échantillon. D'une façon générale, le procédé revient à :

obtenir une estimation initiale de l'amplitude complexe de l'onde 22, transmise par l'échantillon, dans le plan du détecteur, et cela dans plusieurs bandes spectrales (étape 100) ;

- rétropropager chacune de ces amplitudes complexes dans le plan de l'échantillon, pour obtenir, à chaque bande spectrale, une amplitude complexe dans le plan de l'échantillon A° _p (r) (étape 200)

constituer une fonction de pondération de chaque amplitude complexe F _p (r) dans le plan de l'échantillon (étape 300), de manière à réduire l'influence des interfacts de twin image ;

propager ladite fonction de pondération dans le plan du détecteur, pour au moins une bande spectrale, (étape 400) ;

mettre à jour l'estimation de l'amplitude complexe Af _p (r) de l'onde 22, transmise par l'échantillon, dans le plan du détecteur, et cela dans plusieurs bandes spectrales, à l'aide de la fonction de pondération F^ _p (r) propagée dans le plan du détecteur (étape 500).

La formule de mise à jour de l'étape 500 montre qu'à chaque itération, le module mf (r) (respectivement Ai (r)) de l'amplitude complexe normalisée A _p (r) (respectivement de l'amplitude complexe af _p (r)), dans le plan de détection, correspond à celui déterminé, au cours de l'étape 100, par chaque image I, formée par le photodétecteur 16 dans la bande spectrale λ,. Autrement dit, au cours des différentes itérations, le module, dans le plan de détection, de l'amplitude complexe af (r) ou de l'amplitude complexe normalisée A _p (r) n'évolue pas et correspond à celui dérivé de l'intensité mesurée par le photodétecteur. En revanche, l'algorithme tend à faire évoluer, à chaque mise à jour, l'argument des expressions complexe Af _p (r) ou af _p (r), et en particulier l'estimation de la phase <pf _p (r), cette dernière étant considérée comme égale à la phase de la fonction de pondération F _p (r) propagée dans le plan du détecteur, à chaque longueur d'onde λ,. Aussi, selon cet algorithme, chaque itération comprend :

une mise à jour de l'amplitude complexe A° _p (r) de chaque onde lumineuse dans le plan de l'échantillon P ₀ (étape 200) ;

une mise à jour de l'argument de chaque amplitude complexe, Af _p (r) et en particulier de sa phase, dans le plan de détection (étape 500).

Elaboration de la fonction de pondération.

Une première façon de constituer la fonction de pondération est d'effectuer une pondération équitable entre les différentes bandes spectrales λ, considérées.

Par exemple, la fonction de pondération peut prendre la forme F _p (r) = —∑i kiA° _p (r) , ou k , désigne le facteur de pondération, ou poids, attribué à la i ^eme bande spectrale λ,, comme précédemment décrit en liaison avec l'étape 300. Chaque facteur de pondération k _t est positif et peut avoir la même valeur, par exemple 1/3.

Selon une variante, et cela s'applique notamment dans le cas où l'échantillon analysé est coloré, selon une plage spectrale λ ₀ , les modules des amplitudes complexes de premières ondes lumineuses 22, dont les bandes spectrales λ, sont proches de la plage spectrale λο ont une valeur plus élevée que les modules des amplitudes complexes de deuxièmes ondes lumineuses dont les bandes spectrales sont plus éloignées de la longueur d'onde λο. Dans un tel cas de figure, il est préférable de sous-pondérer les amplitudes complexes des premières ondes lumineuses, et de surpondérer les amplitudes complexes des deuxièmes ondes lumineuses.

Par exemple, si l'échantillon est coloré à l'aide d'un colorant bleu, qui correspond dans notre exemple à la première bande spectrale λι, le facteur de pondération ki est plus faible que les facteurs de pondération k ₂ et k ₃ respectivement associés aux bandes spectrales λ ₂ (verte) et λ ₃ (rouge).

Selon une autre variante, le module et l'argument de chaque amplitude complexe sont pondérés par des facteurs de pondérations indépendants, de telle sorte que

ki et k'i étant des facteurs de pondération respectivement associés au module et à l'argument de l'amplitude complexe de l'onde lumineuse 22,, dans le plan de l'échantillon, dans la bande spectrale λ,.

Selon une autre variante, la combinaison des amplitudes complexes A° _p (r) prend la forme d'une somme de logarithmes, selon l'expression : ln(F _p °(r)) = ——∑; /q ln[j4°„(r)l.

Selon une autre variante, on détermine, dans le plan de l'échantillon, non pas une fonction de pondération F _p (r), mais une pluralité de fonctions de pondération F _p j), chaque fonction étant associée à une bande spectrale λ,.

Chaque fonction de pondération F° _p (r), associée à une i ^eme longueur d'onde, est obtenue en combinant plusieurs amplitudes complexes A° _p (r), respectivement associées à différentes bandes spectrales.

Selon un premier exemple, en considérant trois bandes spectrales :

Ainsi, selon ce mode de réalisation, la fonction de pondération prend la forme d'un vecteur F _p (r), de dimension N, N étant le nombre de bandes spectrales considérées, dont chaque terme F? _p (r) est une fonction de pondération associée à une bande spectrale λ,. Cette fonction de pondération peut être obtenue par le produit matriciel suivant :

r) = K Â

Où K est une matrice de pondération, chaque terme ky de la matrice de pondération représentant le poids associé à l'amplitude complexe Α° _ρ (τ) associée à la bande spectrale λ} pour le calcul de la fonction de pondération associée à la bande spectrale λ,.

La matrice K est une matrice carrée de dimension N x N, N étant le nombre de bandes spectrales considérées.

La fonction de pondération est de préférence normalisée, de telle sorte que chaque terme

F _j ⁰ _p s'exprime sous la forme : le terme ^{■ ■} étant un terme de normalisation.

Selon un deuxième exemple de ce mode de réalisation, toujours en considérant trois bandes spectrales,

Ainsi, selon ce mode de réalisation, la fonction de pondération prend la forme d'un vecteur Fp r), de dimension 2N, N étant le nombre de bandes spectrales considérées, dont chaque terme représente soit le module soit l'argument d'une fonction de pondération F° _p r) associée à une bande spectrale λ,. Cette fonction de pondération peut être obtenue par le produit matriciel suivant :

F°(r-) = K °

Où K est une matrice de pondération, de dimension 2N x2N, chaque terme kg de la matrice de pondération représentant le poids associé soit à l'argument, soit à la phase, de l'amplitude complexe Aj _p (r) associée à la bande spectrale λ} .

Selon ce mode de réalisation, chaque coordonnée du vecteur A _p représente soit le module, soit l'argument, d'une amplitude complexe, Aj _p (r) dans une bande spectrale j.

De même que dans l'exemple précédent, la fonction de pondération est de préférence normalisée, de telle sorte que chaque terme F° _p s'exprime sous la forme :

et

6

Arg (¾ (r)) = ¹ Y _j arg( _p (r))

Quel que soit le cas de figure, les coefficients d'une matrice de pondération peuvent être préalablement déterminés, soit de façon arbitraire, soit en se basant sur des essais expérimentaux.

Par exemple, il est possible d'établir un coefficient de régression linéaire oty entre deux composantes i et j du vecteur A _p (r), en considérant une pluralité de positions axiales (r) dans le plan de l'échantillon, de façon à obtenir une statistique suffisante. Le coefficient ky de la matrice de pondération peut alors être déterminé en fonction de ce coefficient de régression linéaire oty, éventuellement affecté d'un terme tenant compte de la dispersion autour du modèle de régression linéaire. Dans un tel cas de figure, la diagonale de la matrice de pondération peut être constituée de coefficients k _û égaux à 1.

Cela permet d'établir une fonction de pondération F° _p , associé à la longueur d'onde λ,, prenant en compte la corrélation entre les différents termes du vecteur Ap (r).

Variantes concernant la source de lumière ou le photodétecteur.

Dans les exemples donnés en lien avec les figures 1 et 2, la source de lumière 11, apte à émettre une onde lumineuse 12 dans différentes bandes spectrales, comporte trois sources de lumière élémentaires lli, II ₂ , II ₃ , prenant la forme de diodes électroluminescentes, émettant respectivement dans une première bande spectrale λι, une deuxième bande spectrale λ ₂ et une troisième bande spectrale λ3, les bandes spectrales étant différentes les unes des autres, et, de préférence, ne se recouvrant pas. La source de lumière 11 peut également comporter une source de lumière blanche ll _w , placée en amont d'un dispositif de filtrage 19, par exemple une roue à filtres, apte à interposer un filtre de bande passante λ, entre la source de lumière blanche et l'échantillon, comme représenté sur la figure 5, de telle sorte que l'image I, formée par le photodétecteur 16 soit représentative de ladite bande passante λ,. Plusieurs filtres, présentant des bandes passantes différentes les unes des autres, sont alors successivement interposés entre la source de lumière ll _w et l'échantillon 10.

Selon une variante, le dispositif de filtrage 19 peut également être un filtre tribande, définissant une pluralité de bandes spectrales. Un exemple de filtre convenant à l'application est le filtre Tri-Band 458, 530 & 628 nm - Edmund Optics, définissant des bandes spectrales respectivement centrées sur les longueurs d'onde 457 nm, 530 nm et 628 nm. Cela permet une illumination de l'échantillon simultanément à l'aide de 3 longueurs d'onde.

L'utilisation d'un diffuseur 17, tel que précédemment décrit, entre la source de lumière et le diaphragme 18 est préférable, et cela quel que soit le mode de réalisation.

Le photodétecteur 16 peut, comme précédemment décrit, être un photodétecteur matriciel RGB, ce qui permet d'acquérir successivement ou simultanément différentes des images ΙΙ..Ι,...ΙΝ différentes bandes spectrales λι..λί...λΝ. Dans ce cas, la source de lumière peut être une source de lumière blanche ll _w , auquel cas les différentes images peuvent être acquises simultanément. Il peut également s'agir d'un photodétecteur 16 monochrome, auquel cas la source de lumière 11 est apte à générer, successivement, une onde lumineuse dans différentes bandes spectrales λι..λί...λΝ. Dans une telle configuration, la source lumineuse comporte soit plusieurs sources de lumière élémentaires lli, II2, II ₃ , soit un dispositif de filtrage 19, comme précédemment décrit. Dans un tel cas de figure, l'échantillon est successivement exposé à des ondes lumineuses incidentes 12Ι..12,...12Ν, N étant le nom bre de bandes spectrales considérées. Une image I, (l≤i≤N), représentative de l'onde lumineuse 22, transmise par l'échantillon est alors acquise à chaque exposition.

Essais réalisés.

Des essais ont été réalisés selon la configuration représentée sur la figure 1 et décrite ci-dessus. L'échantillon est une lame d'anatomopathologie, comportant une coupe de colon coloré par Hematoxyline Eosine Safran. La source de lumière est disposée à une distance Δ égale à 5 cm de l'échantillon, cette distance séparant le diaphragme 18 de l'échantillon 10.

Les figures 6A, 6B et 6C représentent une image d u module |^4° _p (r) | de l'amplitude complexe ^ι _,ρ ") de l'onde 22i transmise par l'échantillon, dans le plan P ₀ de l'échantillon, dans la première bande spectrale λι s'étendant entre 450 et 465 nm, ces images éta nt obtenues après un nom bre d'itérations p respectivement égal à 1, 3 et 10.

Les figures 7 A, 7B et 7C représentent une image du module |^2 _,p ( ^r ) | de l'amplitude complexe ^2 _,p ( ^r ) de l'onde 222 transmise par l'échantillon, dans le plan Po de l'échantillon, dans la deuxième bande spectrale λ2 s'étendant entre 520 et 535 nm, ces images étant obtenues après un nom bre d'itérations p respectivement égal à 1, 3 et 10.

Les figures 8A, 8B et 8C représentent une image du module |^3 _,p (r) | de l'amplitude complexe ° _p (r) de l'onde 223 transmise par l'échantillon, dans le plan Po de l'échantillon, dans la troisième ba nde spectrale λ3 s'étendant entre 620 et 630 nm, ces images étant obtenues après un nom bre d'itérations p respectivement égal à 1, 3 et 10. On remarque que le niveau de gris moyen de ces images est plus important que le niveau de gris des images des figures 6A, 6B, 6C, 7A, 7B et 7C. Cela est dû a u fait que l'échantillon à la coloration rouge-violet de l'échantillon.

Les figures 9A, 9B et 9C représentent respectivement la com binaison des images 6A-7A-8A, 6B- 7B-8B; 6C-7C-8C. Ces figures permettent d'avoir une représentation en couleur de l'échantillon, en prenant en compte simultanément les trois bandes spectrales λι, λ2 et λ3. Sur chaque série d'images, on observe une augmentation du contraste en fonction du nombre d'itérations. On peut également constater la formation d'images dont la résolution spatiale est satisfaisante lorsque le nombre d'itérations est inférieur ou égal à 10, ce qui limite le temps de calcul à quelques secondes. Le procédé est donc apte l'observation d'échantillons, sur champ large, et à cadences élevées. Il permet d'obtenir des images selon une ou plusieurs bandes spectrales, le rendant compatible avec les procédés de coloration couramment pratiqués dans le domaine de l'anatomocytopathologie.