FR 3047077 A1 20170728 FR 1650575 A 20160125 1 Procédé d'observation d'un échantillon par imagerie sans lentille Description DOMAINE TECHNIQUE Le domaine technique de l'invention est lié à l'observation d'un échantillon, en particulier un échantillon biologique, par imagerie sans lentille. L'échantillon peut notamment être un tissu biologique, l'observation étant réalisée à des fins d'analyse histologique. 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 capteur d'image, sans disposer de lentille de grossissement optique entre l'échantillon et le capteur d'image. Ainsi, le capteur d'image 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 de lumière 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 de lumière. 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 W02008090330 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. D'une façon générale, l'image formée sur le capteur d'image, comportant des figures d'interférence, peut être traitée par un algorithme de reconstruction holographique, 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 3047077 2 capteur d'image étant connue, on applique un algorithme de propagation, prenant en compte cette distance, ainsi que la longueur d'onde de l'onde lumineuse émise par la source de lumière. 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. 5 On a récemment démontré l'intérêt de l'imagerie sans lentille pour l'analyse histologique de prélèvements tissulaires et tumoraux prenant la forme d'une tranche fine, ou lamelle, déposée sur un support transparent. Le large champ d'observation conféré par l'imagerie sans lentille, couplé à son faible coût constituent des avantages décisifs par rapport à l'usage d'un microscope classique. De plus, le recours à des algorithmes de reconstruction holographique permet de 10 former une image d'absorption ou une image de phase de l'échantillon observé, d'une qualité suffisante pour effectuer une analyse histologique fiable. La publication Greenbaum A, « wide-field computational imaging of pathology slides using lensfree on-chip microscopy », Sci. Transi. Med 6, 267ra175 (2014), décrit par exemple l'application de l'imagerie sans lentille à des lames d'anatomo-pathologie formées à en prélevant différents 15 types de tissus biologiques. Il en est de même quant à la publication Luo W. « Synthetic aperture- based on-chip microscopy », Light : Science & Applications (2015) 4, e261. Dans cette publication, une huile de couplage optique est appliquée entre le capteur d'image et l'échantillon de façon à améliorer le couplage optique. Cependant, le temps nécessaire à l'obtention d'une image d'un échantillon de 1 mm2 dépasse 45 minutes, ce qui est encore trop 20 important, et peu compatible avec une utilisation en routine. Les inventeurs proposent un procédé d'observation d'un échantillon, et en particulier un prélèvement tissulaire, par imagerie sans lentille, mettant en oeuvre des procédés de reconstruction holographiques, permettant la formation d'une image plus nette. Ce procédé est particulièrement adapté à l'observation de tissus biologiques. 25 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 à émettre une onde lumineuse incidente se propageant vers l'échantillon ; 30 ii) acquisition, à l'aide d'un capteur d'image s'étendant selon un plan de détection, d'au moins une image de l'échantillon, l'échantillon étant disposé entre la source de lumière 3047077 3 et le capteur d'image, chaque image étant représentative d'une onde lumineuse transmise par l'échantillon sous l'effet de ladite illumination ; III application d'un opérateur de propagation à l'image acquise, lors de l'étape ii), de façon à former une image représentative de ladite onde lumineuse transmise par l'échantillon, 5 dans un plan passant par ledit échantillon ; le procédé étant caractérisé en ce qu'il comprend également, préalablement à l'étape ii), l'application d'un liquide d'imprégnation sur l'échantillon (10), de telle sorte que l'échantillon s'imprègne dudit liquide d'imprégnation, ledit liquide d'imprégnation ayant un indice de réfraction strictement supérieur à 1. 10 L'onde lumineuse transmise par l'échantillon correspond à une onde lumineuse à laquelle est exposé le capteur d'image. L'imprégnation de l'échantillon peut être réalisée avant l'étape i) ou entre les étapes i) et ii). De préférence, l'image formée lors de l'étape iii), représentative de l'onde lumineuse transmise par l'échantillon, c'est-à-dire de l'onde lumineuse à laquelle est exposé le capteur d'image, est 15 établie à partir de la phase, dans le plan de l'échantillon, de ladite onde lumineuse transmise par l'échantillon. Ainsi, l'étape iii) permet de déterminer la phase de cette onde lumineuse dans un plan passant par l'échantillon. Le plan de passant par l'échantillon peut être un plan selon lequel s'étend l'échantillon. Il est de préférence sensiblement parallèle au plan de détection. Le procédé peut comprendre l'application d'un filtre passe-haut à l'image formée lors de l'étape 20 iii). De préférence, le liquide d'imprégnation a un indice de réfraction compris entre 1.2 et 1.8 ou entre 1.3 et 1.7. Selon un mode de réalisation, ledit liquide d'imprégnation s'étend de l'échantillon jusqu'au capteur d'image. Le capteur d'image comprend une pluralité de pixels, et le liquide 25 d'imprégnation s'étend de l'échantillon jusqu'à ces pixels. Chaque pixel peut être couplé à une microlentille de focalisation, auquel cas le liquide d'imprégnation peut s'étendre de l'échantillon jusqu'à chaque microlentille. Selon ce mode de réalisation, le liquide d'imprégnation peut être appliqué sur l'échantillon ou sur le capteur d'image, après quoi l'échantillon est rapproché du capteur d'image, jusqu'à ce que le liquide d'imprégnation s'étende de l'échantillon jusqu'au 30 capteur d'image. Selon un mode de réalisation préféré, une lame transparente est disposée entre l'échantillon et le capteur d'image, le liquide d'imprégnation s'étendant entre l'échantillon et ladite lame transparente. 3047077 4 L'échantillon peut notamment être une lamelle de tissu biologique, en particulier apte à la réalisation d'analyses histologiques. Il peut en particulier être déposé sur une lame support transparente. De préférence, l'échantillon est rendu transparent ou translucide après l'application du liquide 5 d'imprégnation. L'application du liquide d'imprégnation a alors pour fonction d'améliorer la transparence de l'échantillon. Selon un mode de réalisation, l'étape ii) comporte la formation d'une pluralité d'images, chaque image étant acquise dans une bande spectrale différente l'une de l'autre l'étape iii) comportant les sous-étapes suivantes : 10 a) détermination, à partir de chaque image acquise selon une bande spectrale, d'une amplitude complexe initiale de l'onde lumineuse transmise, selon ladite bande spectrale, dans ledit plan de détection ; b) à partir de chaque amplitude complexe établie dans ledit plan de détection, selon une bande spectrale, détermination d'une amplitude complexe de l'onde transmise, selon 15 chaque bande spectrale, dans le plan selon lequel s'étend l'échantillon ; c) combinaison d'une pluralité d'amplitudes complexes déterminées lors de la sous-étape b), à différentes bandes spectrales, pour calculer une fonction de pondération, dans le plan de l'échantillon ; d) projection de ladite fonction de pondération dans le plan de détection de façon à 20 obtenir, pour chaque bande spectrale, une fonction de pondération dans ledit plan de détection ; e) mise à jour de chaque amplitude complexe de l'onde lumineuse transmise, selon chaque bande spectrale, dans le plan de détection, à l'aide de de ladite fonction de pondération obtenue, dans ladite bande spectrale, au cours de la sous-étape d) ; 25 f) répétition éventuelle des sous-étapes b) à e) jusqu'à l'atteinte d'un critère d'arrêt. Le procédé peut alors comporter une sous-étape g) de formation d'une image représentative de l'argument, dans le plan de l'échantillon, de l'amplitude complexe de l'onde transmise par l'échantillon selon au moins une bande spectrale (X,). Cette image peut être représentative d'une combinaison, par exemple une somme ou une moyenne, de l'argument de l'amplitude 30 complexe, dans le plan de l'échantillon, de ladite onde transmise, selon différentes bandes spectrales. De préférence, aucune optique de grossissement n'est disposée entre l'échantillon et le capteur d'image. 3047077 5 FIGURES La figure 1A représente un premier exemple de dispositif pour la mise en oeuvre de l'invention, l'échantillon analysé étant une lame d'anatomo-pathologie. Les figures 1B et 1C représentent différentes configurations de l'échantillon. 5 La figure 2A représente le plan de détection, sur lequel une image est formée, ainsi que le plan selon lequel s'étend l'échantillon. Cette figure illustre également les liens entre les principales grandeurs mises en oeuvre dans un algorithme de reconstruction avantageux. La figure 2B représente un logigramme décrivant l'enchaînement des principales étapes de l'algorithme de reconstruction avantageux illustré sur la figure 2A. 10 La figure 3A représente une image d'un échantillon dit de référence, acquise par le capteur d'image. L'échantillon de référence est un échantillon configuré selon l'art antérieur. La figure 3B représente une image d'un échantillon dit de test, acquise par le capteur d'image. L'échantillon de référence et l'échantillon de test sont formés d'une lame de tissu prélevé dans un intestin. L'échantillon de test est configuré conformément à un mode de réalisation de 15 l'invention. Les figures 4A et 4B représentent respectivement une image du module et une image dite de phase d'une onde lumineuse transmise par l'échantillon de référence. Les figures 4C et 4D représentent respectivement un détail des figures 4A et 4B. Les figures 5A et 5B représentent respectivement une image du module et une image dite de 20 phase d'une onde lumineuse transmise par l'échantillon de test. Les figures 5C et 5D représentent respectivement un détail des figures 5A et 5B. Les figures 5E et 5F sont des vues de l'échantillon de test obtenues selon une modalité de microscopie classique, le champ d'observation correspondant respectivement aux détails illustrés sur les figures 5A et 5B. La figure 6A représente une image, dite de phase, de l'échantillon de test. La figure 6B est une 25 image obtenue de ce même échantillon selon une modalité de microscopie classique. 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 à émettre une onde lumineuse 12, dite onde lumineuse incidente, se propageant en 30 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 notamment s'agir d'une lame de tissu destinée à une analyse histologique, ou lame d'anatomopathologie, comportant une fine épaisseur de tissu déposée sur une lame transparente 15. Par 3047077 6 fine épaisseur, on entend une épaisseur de préférence inférieure à 100 um, et de préférence inférieure à 10 um, typiquement quelques micromètres. Un tel échantillon est représenté sur la figure 1. L'échantillon s'étend selon un plan Po, dit plan de l'échantillon, perpendiculaire à l'axe de propagation Z. 5 Cette lame de tissu est obtenue selon des procédés de préparation connus, à partir d'un échantillon tissulaire prélevé par biopsie ou frottis, puis préparé de façon à se présenter sous la forme d'une fine épaisseur déposée sur une lame transparente, cette dernière servant de support. De tels procédés sont connus dans le domaine de l'histologie. Ils comportent par exemple une découpe d'un tissu congelé, ou une inclusion d'un tissu prélevé dans une matrice 10 de paraffine. La distance 4 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 15 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. Dans l'exemple de la figure 1A, la source de lumière 11 comporte trois diodes électroluminescentes (LED, acronyme du terme anglais Light Emitting Diode) 111, 112 et 113 émettant respectivement dans une première bande spectrale Xi = 450nm - 465 nm, une 20 deuxième bande spectrale X2 = 520nm - 535 nm et une troisième bande spectrale X3 = 620nm - 630 nm. Ces trois diodes électroluminescentes constituent des sources élémentaires de la source de lumière 11.Dans cet exemple, la source de lumière est une diode électroluminescente (LED, acronyme de Light Emitting Diode) fournie par CREE sous la référence XLamp MCE. Les trois diodes électroluminescentes élémentaires 111, 112 et 113 la composant sont activées 25 simultanément. De façon alternative, ces diodes électroluminescentes peuvent être activées successivement. 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 um et 1mm, de préférence entre 50 um et 500 um. Dans cet exemple, le diaphragme est fourni 30 par Thorlabs sous la référence P150S et son diamètre est de 150 um. Le diaphragme peut être remplacé par une fibre optique, dont une première extrémité est placée face à la source de lumière 11 et dont une deuxième extrémité est placée face à l'échantillon 10. 3047077 7 Le dispositif comporte de préférence un diffuseur 17, disposé entre la source de lumière 11 et le diaphragme 18. L'usage d'un tel diffuseur permet de s'affranchir de contraintes de centrage de la source de lumière 11 par rapport à l'ouverture du diaphragme 18. La fonction d'un tel diffuseur est de répartir le faisceau lumineux, produit par une source de lumière élémentaire 5 11,, (13) 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 80°. Dans cet exemple, le diffuseur est fourni par Luminit sous la référence Luminit LSD 40°. L'échantillon 10 est disposé entre la source de lumière 11 et un capteur d'image 20. Ce dernier s'étend de préférence parallèlement, ou sensiblement parallèlement à la lame transparente 15 10 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 à 20° ou 10° étant admise. Le capteur d'image 20 est apte à former une image selon un plan de détection P20. Dans l'exemple représenté, il s'agit d'un capteur d'image comportant une matrice de pixels, de type 15 CCD ou un CMOS. Les CMOS sont les capteurs 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 capteur d'image 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 carrés de taille 1.4 um. La surface utile du capteur est de 3.6 x 2.7 mm2. Les capteurs dont la taille des 20 pixels est inférieure à 3 um sont préférés, afin d'améliorer la résolution spatiale de l'image. Le plan de détection P20 s'étend de préférence perpendiculairement à l'axe de propagation Z de l'onde lumineuse incidente 12. Le capteur d'image comporte un filtre de Bayer, de telle sorte que chaque pixel est sensible à une bande spectrale choisie parmi le bleu, le rouge ou le vert. Ainsi, lorsque l'échantillon est 25 exposé à la source de lumière 11, le capteur d'image 20 acquiert une image I pouvant être décomposée en : une première image I1 dans la première bande spectrale Xi d'émission de la première diode électroluminescente 111, cette image étant formée à partir des pixels exposés à une longueur d'onde transmise par le filtre bleu du filtre de Bayer; 30 une deuxième image /2 dans la deuxième bande spectrale X2 d'émission de la deuxième diode électroluminescente 112, cette image étant formée à partir des pixels exposés à une longueur d'onde transmise par le filtre vert du filtre de Bayer; 3047077 8 une troisième image /3 dans la troisième bande spectrale X3 d'émission de la troisième diode électroluminescente 113, cette image étant formée à partir des pixels exposés à une longueur d'onde transmise par le filtre rouge du filtre de Bayer. D'une façon générale, selon ce mode de réalisation, le capteur d'image 20 permet l'acquisition 5 d'images IL de l'échantillon 10 dans différentes bandes spectrales X. Chaque image IL est représentative d'une onde lumineuse 22, transmise par l'échantillon 10, et à laquelle est exposée le capteur d'image, dans chaque bande spectrale X. 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 10 toutefois envisageable. D'autres configurations sont possibles, par exemple l'utilisation d'un capteur d'image monochrome, acquérant une image IL de l'échantillon lorsque ce dernier est successivement illuminé par une onde incidente 12' et cela dans différentes bandes spectrales X. Chaque onde incidente 12, peut être émise par une source de lumière 11, émettant dans une desdites bandes 15 spectrales, ou par une source de lumière blanche filtrée par un filtre optique dont la bande passante correspond à ladite bande spectrale X. La distance d entre l'échantillon 10 et la matrice de pixels du capteur d'image 20 est, dans cet exemple, égale à 300 um. 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 capteur d'image est préférentiellement comprise 20 entre 50 um et 2 cm, de préférence comprise entre 100 um et 2 mm. On remarque l'absence d'optique de grossissement entre le capteur d'image 20 et l'échantillon 10. Cela n'empêche pas la présence éventuelle de microlentilles de focalisation au niveau de chaque pixel du capteur d'image 20, ces dernières n'ayant pas de fonction de grandissement de l'image acquise par le capteur d'image. 25 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 P20, des interférences, en particulier avec une partie de l'onde lumineuse incidente 12, transmise par 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 X, transmise par l'échantillon, et à laquelle est exposé 30 le capteur d'image 20, peut comprendre : une composante résultant de la diffraction de l'onde lumineuse incidente 12, par l'échantillon ; 3047077 9 une composante résultant de l'absorption de l'onde lumineuse incidente 12, dans l'échantillon. Un processeur 30, par exemple un microprocesseur, est apte à traiter chaque image acquise par le capteur d'image 20. En particulier, le processeur est un microprocesseur relié à une mémoire 5 programmable 32 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. Le processeur peut être couplé à un écran 34 permettant l'affichage d'images acquises par le capteur d'image 20 ou calculées par le processeur 30. Le dispositif comporte également un actionneur 40 permettant un déplacement relatif de 10 l'échantillon 10 par rapport au capteur d'image 20, ce déplacement étant de préférence réalisé selon une direction perpendiculaire à l'axe de propagation Z, c'est-à-dire parallèlement au plan de détection P20. Un tel déplacement permet d'effectuer un balayage d'un échantillon de grande surface, typiquement de quelques cm2. Le champ d'observation du capteur d'image est voisin de la taille du capteur CMOS, soit 3.6 x 2.7 mm soit 9.7 mm2, du fait de la faible distance d entre 15 l'échantillon et le capteur. Il est donc nécessaire de déplacer le capteur d'image 20 relativement à l'échantillon 10 pour effectuer une observation complète de ce dernier. Un aspect important de l'invention porte sur la configuration de l'échantillon 10. Préalablement à l'acquisition d'images par le capteur d'image, un liquide d'imprégnation 25 est déposé sur cet échantillon. Ce liquide d'imprégnation 25 est un liquide ou un gel dont l'indice de réfraction est 20 strictement supérieur à 1, et de préférence supérieur à 1.2 Avantageusement, cet indice de réfraction est proche de l'indice de réfraction de l'échantillon, ou compris entre l'indice de réfraction de l'échantillon et l'indice de réfraction du capteur d'image, en particulier des microlentilles de focalisation associées à chaque pixel. Par indice de réfraction proche de celui de l'échantillon, on entend compris entre +50% et - 50% ou +30% et -30% de l'indice de 25 réfraction de l'échantillon. L'indice de réfraction du liquide d'imprégnation est par exemple compris entre 1 et 2. Avantageusement, il est compris entre 1.2 et 1.8. Lorsque l'échantillon 10 comporte un tissu congelé, l'indice de réfraction est de préférence compris entre 1.3 et 1.5. Lorsque l'échantillon 10 est un tissu inclus dans une matrice de paraffine, l'indice de réfraction est de préférence voisin de 1.5, par exemple compris entre 1.4 et 1.7. 30 Le liquide d'imprégnation 25 est placé au contact de l'échantillon 10 et ce dernier, du fait de sa porosité, s'en imprègne. L'échantillon devient alors transparent ou translucide sous l'effet de cette imprégnation. Les inventeurs ont constaté qu'une telle configuration permettait d'obtenir 3047077 10 une image de bonne qualité de l'échantillon, comme décrit ultérieurement, en lien avec les essais expérimentaux réalisés. Le liquide d'imprégnation peut être choisi parmi un agent connu pour éclaircir les tissus. Il peut par exemple s'agir de BABB, acronyme de « Benzyl alcohol / Benzyl Benzoate », agent éclaircicant connu de l'homme du métier. Il peut également s'agir de 5 glycérol, d'éthanol, ou de DMSO, acronyme de Diméthylsulfoxyde. D'autres agents connus pour éclaircir les tissus, dénommés par le terme anglais « clearing agents », sont mentionnés dans la publication Zhu D « Recent Progress in Tissue Optical Clearing » Laser&Photonics Reviews. 2013 7(5) :732-757. Comme illustré sur la figure 1A, le liquide d'imprégnation 25 peut remplir un espace comprenant 10 l'échantillon 10, cet espace étant délimité par la lame 15 précédemment décrite, dite lame support, et une autre lame, dite lame de confinement 16, s'étendant de référence parallèlement, ou sensiblement parallèlement, à la lame support 15 et au plan de détection P20. Cette configuration permet de former des interfaces planes entre les différents milieux traversés par l'onde lumineuse transmise par l'échantillon jusqu'au capteur d'image 20. 15 La lame de confinement 16 est optionnelle et peut être omise, comme représenté sur la figure 1B. Dans cette configuration, le liquide d'imprégnation 25 s'étend entre l'échantillon 10 et le capteur d'images 20. Dans cet exemple de réalisation, le liquide d'imprégnation s'étend jusqu'aux microlentilles de focalisation disposées face à chaque pixel. Cela permet d'assurer une adaptation de l'indice de réfraction depuis l'échantillon jusqu'à chaque pixel. 20 Selon le mode de réalisation représenté sur la figure 1C, le liquide d'imprégnation est disposé au contact de l'échantillon et ne s'étend pas jusqu'au capteur d'image. Cette configuration présente cependant l'inconvénient d'une interface air/liquide d'imprégnation non plane et mal maîtrisée, par rapport aux configurations illustrées sur les figures 1A et 1B. Le recours à un liquide d'imprégnation entre un échantillon et un capteur d'images est connu 25 de l'art antérieur, et en particulier dans la publication de W. Luo précédemment citée. Cependant, un élément important de l'invention est que l'échantillon 10 est placé au contact du liquide d'imprégnation 25 et s'imprègne de ce dernier, contrairement à l'art antérieur, selon lequel le liquide d'imprégnation s'étend jusqu'à une lame transparente délimitant l'échantillon, et n'est pas directement en contact avec l'échantillon. L'imprégnation du liquide d'imprégnation 30 25 dans l'échantillon 10 modifie son aspect visuel en augmentant sa transparence, en réduisant la diffusion de l'onde lumineuse incidente 12, dans l'échantillon. L'invention tire profit que le liquide d'imprégnation 25 pénètre dans l'échantillon 10, ce dernier étant poreux. 3047077 11 Les images acquises sur le capteur d'image 20 ne permettent pas d'obtenir une représentation suffisamment précise de l'échantillon observé, notamment lorsque ce dernier est une lame de tissu. Comme décrit en lien avec l'art antérieur, on peut appliquer, à chaque image I acquise par le capteur d'image, dans une bande spectrale X, un opérateur de propagation h, de façon à 5 calculer une grandeur représentative de l'onde lumineuse 22 transmise par l'échantillon 10, et à laquelle est exposée le capteur d'image 20. Un tel procédé, désigné par le terme reconstruction holographique, permet notamment de reconstruire une image du module ou de la phase de cette onde lumineuse dans un plan parallèle au plan de détection P20, et notamment dans le plan Po selon lequel s'étend l'échantillon. Pour cela, on effectue un produit de 10 convolution entre une image I acquise par le capteur d'image 20, désignée par le terme hologramme, et un opérateur de propagation h. Il est alors possible de reconstruire une amplitude complexe a de l'onde lumineuse 22 en tout point de coordonnées (x, y, z) de l'espace, et en particulier dans un plan situé à une distance I z I du capteur d'image 20, ce plan pouvant être le plan de l'échantillon Po. L'amplitude complexe a est une grandeur complexe 15 dont l'argument et le module sont respectivement représentatifs de la phase et de l'intensité de l'onde lumineuse à laquelle est exposée le capteur d'image. L'opérateur de propagation h(x, y, z) a pour fonction de décrire la propagation de la lumière entre le capteur d'image 20 et un point de coordonnées (x, y, z) situé à une distance I z I du capteur d'image 20. Il est alors possible de déterminer le module M (x, y, z) et/ou la phase cf:1 20 (x,y,z) l'onde lumineuse 22, à cette distance I z I, dite distance de reconstruction, avec : - M(x, y, z) = abs [a (x, y, z)] - ço(x, y, z) = arg [a (x, y, z)] Les opérateurs abs et arg désignent respectivement le module et l'argument. L'expression complexe de l'onde lumineuse 22, à une coordonnée (x,y,z) est donnée par 25 a(x, y, z) = II (x, y) * h(x, y, z) le symbole * désignant un produit de convolution. .\//(x, y) désigne la racine carrée de l'intensité des pixels de coordonnées (x, y) de l'image I acquise dans le plan de détection P20, z désignant ici une distance de reconstruction par rapport au capteur d'image. L'opérateur de propagation peut être par exemple la fonction de Fresnel-Helmholtz, telle que : 30 h(x, y, z) = e j2nAexp Az x ). jAz où X désigne la longueur d'onde. 3047077 12 Lorsque la distance de reconstruction z correspond à la distance d entre le capteur d'image 20 et le plan de l'échantillon Po, on peut obtenir une image I° du module M° et/ou de la phase (p° de l'onde lumineuse 22 transmise par l'échantillon dans le plan de l'échantillon Po, chacune de ces images, dites images reconstruites, constituant une représentation de cet échantillon 10. 5 La suite de la description, relative à des essais expérimentaux montre que l'image représentant la phase de l'onde lumineuse 22 dans le plan de l'échantillon Po comporte une information riche quant à l'échantillon, et peut être utilisée à des fins d'analyse histologique. Les inventeurs ont développé un procédé de reconstruction holographique optimisé, mettant à profit l'acquisition d'images IL de l'échantillon dans plusieurs bandes spectrales X. Pour cela, 10 l'échantillon est éclairé simultanément en activant les trois sources de lumière élémentaires précédemment décrites 111, 112 et 113. Le capteur d'image 20 acquiert une image /, à partir de laquelle trois images I1, 12, 13 sont formées, respectivement selon les première, deuxième et troisième bande spectrale Xi, X2, X3 respectivement associées auxdites sources de lumière élémentaires, comme précédemment décrit. Le procédé mis en oeuvre est un procédé itératif, 15 comportant les étapes décrites ci-dessous, en lien avec les figures 2A et 2B. ère 1 é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 simultanément activée. Le capteur d'image acquiert une image /, à partir de laquelle le processeur forme des image IL correspondant respectivement à chaque 20 bande spectrale À1, l'indice i, relatif à la bande spectrale, étant compris entre 1 et N, N étant le nombre de bandes spectrales considérées. Dans cet exemple, N = 3. La formation de chaque image IL peut comprendre une étape dite de démosaïquation, comportant une interpolation de l'intensité entre les pixels exposés à une même bande spectrale. L'échantillon est placé à une coordonnée axiale z = 0, selon l'axe de propagation Z. On désigne 25 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 P20, tandis que le plan z = 0 correspond au plan de l'échantillon Po. Si IL -d(r)=1"1(r) désigne la valeur de l'intensité collectée, dans la bande spectrale par le pixel du capteur d'image de coordonnée radiale r dans le plan de détection P20, il est possible 30 d'établir, à l'aide de l'image IL une amplitude complexe al-d(r) = ad (r) (r) de l'onde 22, au niveau dudit pixel de coordonnée r, dont le module peut être exprimé selon l'expression : 3047077 13 Md (r) = ,j1 "I (r) L'exposant d exprime le fait que l'amplitude complexe est déterminée selon le plan de détection P20, d'équation z = d. Comme précédemment évoqué, l'amplitude complexe ad (r) (r) comporte un module et un argument, de telle sorte que : 5 dît (r) = (r)e-RPi (r) où : Md (r) désigne le module de l'amplitude complexe de l'onde lumineuse 22, détectée par le capteur d'image, dans la Me bande spectrale .11, à une coordonnée radiale r dans le plan de détection ; 10d (f) (r) désigne la phase de l'amplitude complexe de l'onde lumineuse 22, détectée par le capteur d'image, dans la Me bande spectrale À1 et à ladite coordonnée radiale r dans plan de détection. Or, le capteur d'images ne fournit aucune information relative à la phase de l'onde lumineuse 22,. Aussi, lors de l'étape 100, on considère que e -19f l(r) est égal à une valeur initiale arbitraire, 15 par exemple égale à 1. L'amplitude complexe dît (r) peut être exprimée de façon normalisée, selon l'expression : (r) 11`-1, (r) - 'imean où : Irne an désigne l'intensité moyenne d'une image formée dans la Me bande spectale X ; 20 cette intensité moyenne peut être déterminée expérimentalement, en disposant le le capteur d'image 20 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 capteur d'image 20 ; (r) désigne l'amplitude complexe normalisée de l'onde lumineuse 22, détectée par 25 le capteur d'image 20 dans la Me bande spectrale X. La normalisation peut également être réalisée en divisant l'amplitude complexe ad (r) (r) par irnean - ), ou sa racine carrée, ce terme représentant l'intensité lumineuse, à la coordonnée radiale r, mesurée en l'absence d'échantillon. 3047077 14 L'amplitude complexe normalisée Act (r) comporte un module et un argument, de telle sorte que : Ad (r) = md (r)eW (r) où : mci (r) désigne le module de l'amplitude complexe normalisée Act (r) dans le plan de 5 détection ; d (r) désigne la phase de l'amplitude complexe normalisée, qui est également la phase de l'amplitude complexe ad (r), (r), dans le plan de détection. La première étape 100 permet d'attribuer à chaque amplitude complexe ad (r) (r) ou à chaque amplitude complexe normalisée Act (r) une valeur initiale à partir de l'image IL acquise par le 10 capteur d'image 20 dans la Me bande spectrale X, de telle sorte que : a i,p=1(r) = M (r) = \il(r) OU (r) Adi,p=1(r) = mid (r) = mean L L'indice p correspond le rang de l'itération du procédé itératif décrit par la suite. L'étape 100 15 étant une étape d'initialisation, on attribue la valeur 1 à cet indice. En adressant tout ou partie des pixels du capteur d'image 20, on obtient une image complexe, ou champ complexe, de l'onde lumineuse 22, au niveau du plan de détection P20, cette image rassemblant les amplitudes complexes ad (r) ou les amplitudes complexes normalisées Act (r). Dans la suite de la description, on ne considère que l'amplitude complexe normalisée Act (r), 20 sachant que le raisonnement s'applique également à l'amplitude complexe ad (r). Cette première étape est répétée pour chaque bande spectrale (X1...XN) détectée par le capteur d'image 20. 2ime étape : rétro-propagation dans le plan P10 de l'échantillon Durant une deuxième étape 200, on estime l'amplitude complexe normalisée 11°p(r) de l'onde 25 22, à laquelle est exposé le capteur d'image, dans le plan de l'échantillon Po. Cette estimation est réalisée par rétro-propagation de l'amplitude complexe normalisée Acilp(r) déterminée dans le plan de détection P20, cette rétropropagation étant effectuée du plan de détection P20 vers le plan de l'échantillon Po. 3047077 15 Lors de la première itération (p = 1), on utilise l'amplitude complexe normalisée Acilp=i(r) = Ad(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. Comme précédemment décrit, en effectuant un produit de convolution entre l'amplitude 5 complexe de l'onde lumineuse 22,, relative à la bande spectrale À1, déterminée dans le plan de détection z = d, et un opérateur de propagation h, 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 AL p (r) de l'onde lumineuse 22, peut s'obtenir, à un point de coordonnées (r, z), à 10 partir de Aîp=d(r), selon l'opération : g,,p(r) = AU,d (r) * z - d), Où hi, désigne l'opérateur de propagation h dans la bande spectrale Ài. 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 j2n-- 2 15 h(r, , z) = e Âexp -). Riz où X désigne la longueur d'onde. Ainsi, Aîe (r) = 119,p(r) = AU,d (r) * hAi(r , -d) 1 - j27c± (r - r')2 = - j Àide if Aq,,p(r) exp) dr' À1d 20 où r' désigne les coordonnées radiales dans le plan P20 du capteur d'image, d'équation z = d, r désigne les coordonnées radiales dans le plan de l'échantillon Po, qui est ici le plan de reconstruction, d'équation z = 0, 25 X désigne la longueur d'onde centrale de la bande spectrale considérée. Lorsque la reconstruction s'effectue selon le sens la propagation de la lumière, par exemple de l'échantillon 10 vers le capteur d'image 20, on parle de propagation. Lorsque la reconstruction s'effectue selon le sens inverse de la propagation de la lumière, par exemple du capteur d'image 20 vers l'échantillon 10, on parle de rétro-propagation. Aci"p(r) est donc obtenu par rétro- 3047077 16 propagation de il'ilp(r) selon la distance d séparant le plan de détection P20 du plan de l'échantillon Po. Cette deuxième étape est répétée pour chaque bande spectrale (Xi...XN) émise par la source de lumière 11 ou, de façon plus générale, pour chaque bande spectrale (Xi...XN) respectivement 5 associée à chaque image (11...IN) détectée par le capteur d'image 20. Il est possible, à ce stade, d'établir une image du module ou de la phase de l'amplitude complexe Aci"p(r) 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 11°p(r) aux différentes coordonnées r dans le plan de l'échantillon Po. 10 Chaque image du module de l'amplitude complexe 11°p(r) est représentative de l'intensité de l'onde lumineuse dans le plan de l'échantillon Po, tandis que chaque image de l'argument de l'amplitude complexe 11°p(r) est représentative de la phase (p° de l'onde lumineuse 22, dans le plan Po de l'échantillon. L'amplitude complexe normalisée 11°p(r) correspond à une fonction de transmission de l'onde incidente 12, par l'échantillon 10 à la coordonnée radiale r. 15 3ime étape : détermination de la fonction de pondération Au cours de l'étape 300, on constitue, dans le plan de l'échantillon Po, une fonction de pondération, notée Fpc)(r), de l'amplitude complexe de l'onde lumineuse transmise par l'échantillon dans les différentes bandes spectrales X, considérées. Selon cet exemple, la fonction de pondération Fpc)(r), dans le plan de l'échantillon est commune à chaque bande 20 spectrale. Elle est obtenue en combinant les amplitudes complexes normalisées 11°p(r) de l'onde lumineuse transmise par l'échantillon, dans le plan Po de l'échantillon et dans les différentes bandes spectrales X. Dans cet exemple, elle 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 : 25 Fp° (r) = 1 killci"p (r) où ki désigne un facteur de pondération associé à la Me bande spectrale Xi. Les facteurs de pondération peuvent être égaux les uns aux autres, par exemple égaux à 1/3. 4ime étape : propagation de la fonction de pondération dans le plan du capteur d'image 3047077 17 L'étape 400 vise à propager, du plan de l'échantillon Po vers le plan du capteur d'image P20, la fonction de pondération F°(r), déterminée lors de l'étape précédente. L'opérateur de propagation h étant dépendant de la longueur d'onde, cette propagation est effectuée pour chaque bande spectrale À1 considérée. 5 Ainsi, pour chaque bande spectrale .11, F idp (r) = Fp° (r) * z = d). Lorsque l'opérateur de propagation est un opérateur de Fresnel-Helmholtz tel que précédemment défini, j2n-c Fidp(r) = I1id e if Fp(r') exp(jn- r')2)dr' L'opérateur de propagation dépendant de la longueur d'onde, on détermine autant de fonctions 10 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 P20 du capteur d'image, d'équation z = d, r désigne les coordonnées radiales dans le plan de l'échantillon Po, qui est ici le plan de reconstruction, d'équation z = 0, 15 X désigne la longueur d'onde centrale de la bande spectrale considérée. 5'me étape : Mise à jour de l'amplitude complexe dans le plan du capteur d'image Dans l'étape 500, on utilise la valeur de la fonction de pondération Figp(r), dans le plan de détection z = d, pour mettre à jour, dans le même plan, l'estimation de l'amplitude complexe normalisée Acil,p(r) de l'onde lumineuse 22' à laquelle est exposé le capteur d'image 20, dans la 20 bande spectrale X. La formule de mise à jour est : Fcl (r) (r) = (r)d ) x (- ml d (r) x e-19',P(r) t,p r où : Figp(r)1 désigne le module de Figp(r) ; 25 (r) désigne le module de l'amplitude complexe initiale normalisée Agp=1 déterminé, Ad p_1 à partir de chaque image 1,, lors de la première étape 100. Ce terme a une fonction d'attache aux données mesurées ; Cep désigne une estimation de la phase de l'amplitude complexe de l'onde 22, dans la Me bande spectrale .11 ; 3047077 18 Agp(r) désigne l'amplitude complexe de l'onde lumineuse 22, transmise par l'échantillon 10, dans le plan du le capteur d'image 20, cette amplitude complexe constituant la base de l'itération suivante, de rang p + 1. A la suite de cette étape, une nouvelle itération peut commencer, la donnée d'entrée de cette 5 nouvelle itération de rang p + 1 étant il'ilp+i(r) = il'ilp(r), cette nouvelle itération débutant par la rétro-propagation de chaque amplitude complexe normalisée il'ilp+i(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 nombre d'itérations pmax prédéterminé, soit jusqu'à l'atteinte d'un critère de convergence, ce dernier pouvant être, par 10 exemple, exprimé sous la forme d'un écart entre l'estimation de deux mêmes grandeurs entre deux itérations successives. A l'issue du procédé, on dispose d'une estimation de l'amplitude complexe de l'onde lumineuse 22,, transmise par l'échantillon, et à laquelle est exposé le capteur d'image 20, dans le plan de détection P20, d'équation z = d et/ou dans le plan de l'échantillon Po, d'équation z = 0, et 15 cela pour chaque bande spectrale considérée. En utilisant les différentes amplitudes complexes reconstruites Aci"p(r) dans le plan de l'échantillon Po, 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. Cet algorithme, bien que décrit en relation avec une amplitude complexe normalisée Al , 20 s'applique également à l'amplitude complexe non normalisée ai . A partir de l'amplitude complexe obtenue dans le plan de l'échantillon dans les différentes bandes spectrales À1, il est possible d'établir une image Ip représentant : - soit l'intensité, dans le plan Po de l'échantillon, de l'onde 22, transmise par l'échantillon dans une bande spectrale Ài , 25 soit la phase (p? , dans le plan Po de l'échantillon, de l'onde 22, transmise par l'échantillon dans une bande spectrale À. Il est également possible de combiner les modules ou les arguments de l'amplitude complexe 4, dans différentes bandes spectrales À1, par exemple sous la forme d'une moyenne, de façon à former une image /° représentant l'intensité ou la phase (p° de l'onde lumineuse 22 30 transmise par l'échantillon 10. Essais réalisés. 3047077 19 Des essais ont été réalisés en utilisant un échantillon obtenu à partir d'une biopsie d'un colon humain. Suite à la biopsie, le tissu a été immergé dans du formaldéhyde, ayant une fonction d'agent fixateur, puis déshydraté en le plongeant dans des solutions comportant une concentration croissante en alcool. Il a ensuite été inclus dans une matrice de paraffine puis 5 découpé selon de fines lamelles d'épaisseur comprise entre 3 et 4 um. Lors de ces essais, le liquide d'imprégnation 25 utilisé est une huile à immersion Olympus, dont l'indice de réfraction s'élève à 1.516 à 23°C. La lamelle de paraffine, formant l'échantillon test 10test, a ensuite été déposée sur une lame support 15, et une goutte de liquide d'imprégnation 25 a été appliquée sur la lamelle. L'échantillon 10 a ensuite été rapproché du capteur d'image 20 précédemment 10 décrit, de telle sorte que le liquide d'imprégnation s'étende entre ce capteur et l'échantillon, selon la configuration illustrée sur la figure 1B. Un autre lamelle de tissu n'a pas été mouillée par le liquide d'imprégnation et a également été observée, de façon à constituer un échantillon de référence 10ref. Les figures 3A et 3B représentent une image acquise par le capteur d'image 20, désignée par le 15 terme hologramme, en considérant respectivement l'échantillon de référence 10ref et l'échantillon test 10test. L'image acquise est ici une image dans la première bande spectrale Xi (en l'occurrence 450 nm -465 nm), représentée ici en niveaux de gris. On observe que l'image de l'échantillon de test (figure 3B) est nettement moins contrastée que l'image de l'échantillon de référence (figure 3A), l'échantillon de test ayant été rendu transparent par le liquide 20 d'imprégnation 25. Les figures 4A et 4B représentent les résultats de reconstructions holographiques obtenues à partir de l'hologramme de l'échantillon de référence 10ref illustré sur la figure 3A, en mettant en oeuvre l'algorithme de reconstruction décrit en lien avec les figures 2A et 2B. La figure 4A représente le module m° 1(r) de l'amplitude complexe normalisée Aci)=1(r), dans le plan de 25 l'échantillon Po, tandis que la figure 4B représente la phase e_1(r) de l'amplitude complexe normalisée 11?=1(r) dans le plan de l'échantillon Po. Aucune de ces figures n'est vraiment exploitable. Les figures 4C et 4D correspondent respectivement à des détails des images 4A et 4B, dans les régions d'intérêt délimitées par un cadre en pointillés. 30 Les figures 5A et 5B représentent les résultats de reconstructions holographiques obtenues à partir de l'hologramme de l'échantillon de test 10test illustré sur la figure 3B. La figure 5A représente le module m° 1(r) de l'amplitude complexe normalisée Aci)=1(r), dans le plan de 3047077 20 l'échantillon Po, tandis que la figure 5B représente la phase (fel(r) de l'amplitude complexe normalisée Aci)=1(r) dans le plan de l'échantillon Po. La figure 5A ne comporte pas d'information exploitable, ce qui est attendu vu que l'échantillon de test 10test a été rendu transparent par le liquide d'imprégnation 25. En revanche, la figure 5B est particulièrement riche en information 5 et constitue une représentation exploitable de l'échantillon observé. Les figures 5C et 5D correspondent à des détails de l'image 5B, dans les régions d'intérêt délimitées, sur la figure 5B, par un cadre en pointillés. Le procédé de reconstruction mis en oeuvre est rapide: la durée de l'acquisition s'élève à 0.5 seconde, tandis que la durée de la reconstruction est de l'ordre de la seconde, soit 1.5 seconde pour obtenir une image exploitable, le champ d'observation étant de 10 9.7 mm2. L'échantillon de test 10test a été observé par microscopie après coloration HES, acronyme de Hématéine - Eosine - Safran, coloration courante dans le domaine de l'histologie. Les figures 5E et 5F représentent des images, correspondant sensiblement aux régions d'intérêt respectivement observées sur les figures 5C et 5D. La cohérence entre l'image de phase 15 reconstruite et l'image obtenue au microscope est manifeste, ce qui atteste de la fiabilité et de l'intérêt de l'invention. La position relative de l'échantillon de test 10test et du capteur d'image 20 a été modifiée, grâce à la platine de translation 40. On a obtenu, par balayage, une représentation de l'ensemble de l'échantillon de test. La figure 6A représente une image de la phase e_1(r) de l'amplitude 20 complexe normalisée Aci)=1(r), dans la première bande spectrale Xi, reconstruite dans le plan de l'échantillon Po, après application d'un filtre passe-haut, de type filtre de variance. Un tel filtre permet de calculer la valeur de la variance dans le voisinage de chaque pixel, la valeur du pixel étant remplacée par la valeur de la variance ainsi calculée. Le voisinage de chaque pixel peut être par exemple une zone de 5 par 5 pixels. D'autres filtres, de type passe haut, par exemple 25 des filtres utisés pour la détection de contour, seraient utilisables. La figure 6B représente une observation de l'échantillon de test au microscope classique, ce qui constitue une méthode de référence. La cohérence entre ces deux images est élevée. D'une manière générale, une image de la phase de l'onde transmise par l'échantillon, c'est-à-dire de l'onde à laquelle est exposé le capteur d'image, représente une distribution spatiale de 30 l'indice de réfraction de l'échantillon, et cela dans la mesure où l'épaisseur de l'échantillon est considérée comme constante, ce qui est le cas lors de l'observation de lamelles de tissus biologiques. 3047077 21 Le procédé décrit ci-dessus permet d'obtenir une représentation d'un échantillon exploitable à pour l'histologie et l'aide au diagnostic, et cela en apportant les avantages suivants : par rapport à un microscope, le champ d'observation est plus élevé, ce qui permet l'observation d'un échantillon de grande surface, bien plus rapidement : 5 le procédé ne comporte pas de coloration de l'échantillon ; par rapport à d'autres méthodes holographiques, l'imprégnation d'un échantillon par un liquide d'imprégnation permet d'obtenir des hologrammes dont la qualité est suffisante pour obtenir des images reconstruites permettant une caractérisation correcte, et rapide de l'échantillon. 10 L'invention pourra être appliquée pour l'observation d'un échantillon tissulaire en histologie mais également, en dehors de la biologie, d'un échantillon poreux et dont la transparence est augmentée sous l'effet d'une imprégnation par le liquide d'imprégnation. Les applications visées peuvent être dans les domaines de l'agroalimentaire, l'analyse d'échantillons industriels, par exemple des filtres, ou dans le domaine du contrôle de l'environnement. 15 REVENDICATIONS1. Procédé d'observation d'un échantillon (10), comportant les étapes suivantes : i) illumination dudit échantillon à l'aide d'une source de lumière (11), apte à émettre une onde lumineuse incidente (12) se propageant vers l'échantillon (10); acquisition, à l'aide d'un capteur d'image (20) s'étendant selon un plan de détection (P20), d'au moins une image (I, I,) de l'échantillon, l'échantillon étant disposé entre la source de lumière (11) et le capteur d'image (20), chaque image étant représentative d'une onde lumineuse (22, 22i) transmise par l'échantillon sous l'effet de ladite illumination ; iii) application d'un opérateur de propagation à l'image acquise (I, Ii), lors de l'étape ii), de façon à former une image (I°, Ip) représentative de ladite onde lumineuse (22, 22,) transmise par l'échantillon, dans un plan (Po) passant par ledit échantillon ; le procédé étant caractérisé en ce qu'il comprend également, préalablement à l'étape ii), l'application d'un fluide dit d'imprégnation (25) sur l'échantillon (10), de telle sorte que l'échantillon s'imprègne dudit liquide d'imprégnation, ledit liquide d'imprégnation ayant un indice de réfraction strictement supérieur à 1. 2. Procédé selon la revendication 1, dans lequel l'image (I°, formée lors de l'étape iii), représentative de l'onde lumineuse (22, 22,) transmise par l'échantillon, est établie à partir de la phase ((p° , (p?), dans le plan (Po) passant par l'échantillon, de ladite onde lumineuse (22, 22,) transmise par l'échantillon. 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel le liquide d'imprégnation (25) a un indice de réfraction compris entre 1.2 et 1.8 ou entre 1.3 et 1.7. 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel ledit liquide d'imprégnation (25) s'étend de l'échantillon (10) jusqu'au capteur d'image (20). 5. Procédé selon la revendication 4, dans lequel le capteur d'image (20) comprenant une pluralité de pixels, ledit liquide d'imprégnation s'étend de l'échantillon (10) jusqu'auxdits pixels. 6. Procédé selon l'une quelconque des revendications 4 ou 5, dans lequel le liquide d'imprégnation (25) est appliqué sur l'échantillon (10) ou sur le capteur d'image (20), après quoi 3047077 23 l'échantillon est rapproché du capteur d'image (20), jusqu'à ce que le liquide d'imprégnation s'étende de l'échantillon jusqu'au capteur d'image (20). 7. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel une lame transparente (16) est disposée entre l'échantillon (10) et le capteur d'image (20), le liquide d'imprégnation 5 (25) s'étendant entre l'échantillon et ladite lame transparente (16). 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'échantillon (10) est une lamelle de tissu biologique. 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'échantillon (10) est rendu transparent ou translucide après l'application du liquide d'imprégnation (25). 10 10. Procédé selon l'une quelconque des revendications précédentes dans lequel l'étape ii) comporte la formation d'une pluralité d'images (Ii), chaque image étant acquise dans une bande spectrale (X,) différente l'une de l'autre, l'étape iii) comportant alors les sous étapes suivantes : a) détermination, à partir de chaque image acquise (Ii) selon une bande spectrale (XI d'une amplitude complexe initiale (11`iip=1(r)) de l'onde lumineuse transmise (22,), selon 15 ladite bande spectrale, dans ledit plan de détection (P20) ; b) à partir de chaque amplitude complexe ( il'ilp=1(r),Acilp(r)) établie dans ledit plan de détection (P20), selon une bande spectrale (XI détermination d'une amplitude complexe (11?p(r)) de l'onde transmise (22,), selon chaque bande spectrale, dans le plan (Po) selon lequel s'étend l'échantillon ; 20 c) combinaison d'une pluralité d'amplitudes complexes (11?p(r)) déterminées lors de la sous-étape b), à différentes bandes spectrales, pour calculer une fonction de pondération (Fpc)(r)) ,dans le plan de l'échantillon (Po) ; d) projection de ladite fonction de pondération (Fpc)(r)) dans le plan de détection (P20) de façon à obtenir, pour chaque bande spectrale (XI une fonction de pondération (Fpdi(r)) 25 dans ledit plan de détection (P20) ; e) mise à jour de chaque amplitude complexe (11`iip(r)) de l'onde lumineuse (22,) transmise, selon chaque bande spectrale (XI dans le plan de détection (P20), à l'aide de 3047077 24 de ladite fonction de pondération (Fid,p(r)) obtenue, dans ladite bande spectrale (Xi), au cours de la sous-étape d) ; f) répétition des sous-étapes b) à e) jusqu'à l'atteinte d'un critère d'arrêt. 11. Procédé selon la revendication 10, comportant une sous étape g) de formation d'une image 5 (/, I.) représentative de l'argument, dans le plan de l'échantillon (Po), de l'amplitude complexe de l'onde (22;) transmise par l'échantillon (10) selon au moins une bande spectrale (Xi). 12. Procédé selon l'une quelconque des revendications précédentes, dans lequel aucune optique de grossissement n'est disposée entre l'échantillon (10) et le capteur d'image (20).