Brevet FR3097639A1 — publié le 20201225

Description Titre de l'invention : Procédé de reconstruction holographique Domaine technique [0001] Le domaine technique de l'invention est la reconstruction d'images holographiques, en particulier pour caractériser un échantillon, par exemple un échantillon biologique. Technique antérieure [0002] 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 (l'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. [0003] 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 ». [0004] 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. [0005] Afin d'obtenir une observation satisfaisante de l'échantillon, des algorithmes itératifs de reconstruction d'image ont été développés, tels que ceux décrits dans W02016189257 ou dans W02017162985. Ces algorithmes comportent une application itérative d'un opérateur de propagation holographique, de façon à propager l'hologramme, formé dans le plan de détection, dans un plan de reconstruction, ce dernier correspondant généralement un plan de l'échantillon, selon lequel s'étend l'échantillon. Le plan de l'échantillon est généralement parallèle au plan de détection. Les algorithmes décrits dans l'art antérieur comportent des successions de propagation / rétropropagation, d'images entre le plan de détection et le plan de l'échantillon. En effet, l'image acquise par le capteur d'image ne comporte pas d'information relative à la phase de l'onde lumineuse d'exposition. L'objectif de ces algorithmes est d'estimer, de 2 façon itérative, la phase de l'onde lumineuse d'exposition clans le plan de détection. Cela permet de former une image correcte de l'échantillon, dans le plan de reconstruction. Ainsi, ces algorithmes permettent d'obtenir des propriétés optiques de Fonde lumineuse d'exposition. 11 peut s'agir par exemple du module ou de la phase. [0006] Les inventeurs proposent un procédé d'observation d'un échantillon par une méthode d'imagerie holographique, le procédé comprenant une étape de reconstruction d'une image complexe de l'échantillon, sur la base de laquelle on peut obtenir une représentation spatiale de paramètres de l'échantillon. Exposé de l'invention [0007] Un premier objet de l'invention est un procédé d'observation d'un échantillon, l'échantillon s'étendant selon un plan d'échantillon définissant des positions radiales, des paramètres de l'échantillon étant définis à chaque position radiale, le procédé comportant : a. illumination de l'échantillon à l'aide d'une source de lumière, apte à émettre une onde lumineuse incidente se propageant vers l'échantillon ; b. acquisition, à l'aide d'un capteur d'image, d'une image de l'échantillon, formée dans un plan de détection, l'échantillon étant disposé entre la source de lumière et le capteur d'image; c. à partir de l'image acquise par le capteur d'imagc, obtention d'une image de l'échantillon correspondant à une distribution d'au moins un paramètre de l'échantillon dans le plan d'échantillon ; [0008] le procédé étant caractérisé en ce que l'étape c) comporte les étapes suivantes: prise en compte une image d'initial isation; application d'un opérateur de propagation holographique à l'image d'initialisation, ou à une image complexe de l'échantillon résultant de l'étape (v) d'une itération précédente, de façon à estimer une image dans le plan de détection, par exemple l'image acquise par le capteur d'image dans le plan de détection ; comparaison de l'image acquise lors de l'étape h) et de l'image estimée lors de l'étape ii); iv. en fonction de la comparaison mise à jour des paramètres décrivant l'échantillon; v. tant qu'un critère de convergence ou un nombre prédéterminé d'itérations n'a pas été atteint, mise à jour d'une image complexe de l'échantillon, dans le plan de l'échantillon, à partir des paramètres mis à jour lors de l'étape iv), puis réitération des étapes (ii) à (y) ; vi. utilisation de paramètres de l'échantillon, résultant de la dernière étape (iv) 3 effectuée (c'est-à-dire la dernière étape (iv) avant l'étape (vi)), pour alimenter un algorithme d'intelligence artificielle à apprentissage supervisé, de façon à actualiser des paramètres de l'échantillon; vii. utilisation des paramètres résultant de l'étape (vi) pour mettre à jour l'image d'initialisation; viii. réitération, au moins une fois, des étapes (ii) à (v), en utilisant, lors de la première réitération, l'image d'initial isation mise à jour lors de l'étape (vii); ix. obtention d'une image de l'échantillon à partir des paramètres mis à jour lors l'étape iv) de la dernière itération des étapes (ii) à (y). [0009] L'étape iv) peut comporter un calcul d'un indicateur de validité, de telle sorte que les paramètres de l'échantillon sont mis à jour de façon à faire tendre l'indicateur de validité vers une valeur prédéterminée. Lors de l'étape iv), les paramètres de l'échantillon sont alors mis à jour de façon à minimiser l'indicateur de validité. [0010] Par image complexe de l'échantillon, on entend une image complexe d'une onde lumineuse d'exposition, dans le plan de l'échantillon, l'onde lumineuse d'exposition se propageant jusqu'au capteur d'image. [0011] L'algorithme d'intelligence artificielle à apprentissage supervisé peut par exemple mettre en oeuvre un réseau de neurones. Le réseau de neurones peut notamment être un réseau de neurones convolutif. [0012] L'étape iv) peut comporter une détermination d'un gradient de l'indicateur de validité en fonction d'au moins un paramètre, de telle sorte que les paramètres sont mis à jour pour réduire l'indicateur de validité de l'itération suivante. L'étape iv) peut notamment mettre en oeuvre un algorithme de type descente de gradient. [0013] Selon un mode de réalisation : l'étape iv) résulte en la détermination d'une image de chaque paramètre dans le plan de l'échantillon ; et/ou lors de l'étape vi), au moins une image d'un paramètre forme une couche d'entrée de l'algorithme d'intelligence artificielle à apprentissage supervisé ; et/ou lors de l'étape viii), l'algorithme d'intelligence artificielle à apprentissage supervisé délivre une image de sortie, correspondant à une image d'un paramètre de l'échantillon actualisé par l'algorithme. [0014] Par image d'un paramètre, on entend une distribution spatiale du paramètre dans le plan de l'échantillon. [0015] Selon un mode de réalisation, lors de l'étape (iii), le calcul de l'image de l'échantillon dans le plan de détection comporte l'application d'un produit de convolution par un noyau de convolution, le noyau de convolution représentant une étendue spatiale de la source de lumière. [0016] Selon un mode de réalisation, les paramètres décrivant l'échantillon comportent : un premier paramètre, représentant une absorbance de l'échantillon ; un deuxième paramètre, représentant une différence de trajet optique, selon Faxe de propagation de Fonde lumineuse incidente. [0017] Selon un mode de réalisation, lors de l'étape (vi), l'algorithme d'intelligence arti- ficielle à apprentissage supervisé est alimenté par : une image du premier paramètre, correspondant à une distribution spatiale du premier paramètre, mise à jour lors de la dernière itération des étapes (ii) à (v) précédent l'étape (vi); une image du deuxième paramètre, correspondant à une distribution spatiale du deuxième paramètre, mise à jour lors de la dernière itération des étapes (ii) à (y) précédent l'étape (vi) ; [0018] L'algorithme d'intelligence artificielle à apprentissage supervisé permet alors d'obtenir une image du deuxième paramètre actualisé. [0019] Selon un mode de réalisation, suite à une étape (viii), les étapes (vi) à (viii) sont re- nouvelées au moins une fois, voire au moins deux fois. Ainsi, suite à chaque étape (viii), une série d'itérations des étapes (ii) à (y) est effectuée. Après chaque série d'itérations : soit les étapes (vi) à (viii) sont renouvelées, préalablement à une autre série d'itérations des étapes (ii) à (y) ; soit l'étape (ix) est mise en oeuvre. [0020] Selon un mode de réalisation, aucune optique de formation d'image n'est disposée entre l'échantillon et le capteur d'image. [0021] Selon un mode de réalisation, un système optique, de type lentille ou objectif, est disposé entre l'échantillon et le capteur d'image, le système optique définissant un plan image et un plan objet, le procédé étant tel que lors de l'étape b) : le plan objet est décalé du plan d'échantillon selon une distance de défoca- lisation objet ; et/ou le plan image est décalé du plan de détection selon une distance de défocalisation image. [0022] Le procédé peut comporter : [0023] d) caractérisation de l'échantillon à partir de l'image de l'échantillon résultant de l'étape c), ou à partir de chaque image de l'échantillon résultant de l'étape c). [0024] Un deuxième objet de l'invention est un dispositif d'observation d'un échantillon comportant : une source de lumière, configurée pour émettre une onde lumineuse incidente pour illuminer l'échantillon ; un support d'échantillon, configuré pour recevoir l'échantillon ; un capteur d'image, configuré pour acquérir une image de l'échantillon lorsque l'échantillon est disposé sur le support d'échantillon ; [0025] un processeur, programmé pour exécuter des instructions permettant de mettre en oeuvre les étapes i) à viii) ou i) à ix) d'un procédé selon le premier objet de l'invention à partir d'une image acquise par le capteur d'image. [0026] L'invention sera mieux comprise à la lecture de l'exposé des exemples de réalisation présentés, dans la suite de la description, en lien avec les figures listées ci-dessous. Brève description des dessins [0027] [fie. 1A] La figure lA représente un dispositif permettant une mise en oeuvre de l'invention, selon une configuration en imagerie sans lentille. [0028] [fie.1B] La figure 1B représente un dispositif permettant une mise en oeuvre de l'invention, selon une configuration défocalisée. [0029] [fig.2] La figure 2 montre l'évolution d'une onde lumineuse et illustre le retard induit par une différence de trajet optique dans l'échantillon. Le retard génère un déphasage de l'onde lumineuse se propageant entre un échantillon et le capteur d'image. [0030] [fig.3A] La figure 3A représente les principales étapes d'un mode de réalisation de l'invention. [0031] [fig.313] La figure 3B illustre une indétermination d'une propriété optique de I 'échantillon. [0032] [fig.3C] La figure 3C schématise un réseau de neurones convolutif mis en oeuvre dans l'invention. [0033] [fig.4A] [0034] [fig.413] [0035] [fig.4C1 [0036] [fig.4D] [0037] [fig.4E] [0038] [fig.4F] [0039] [fig.4G] Les figures 4A à 4G représentent des images illustrant la phase d'apprentissage du réseau de neurones à l'aide d'un échantillon modélisé. [0040] [fig.5A] [0041] [fig.5131 [0042] [fig.5C] [0043] [fig.5D] [0044] [fig.5E] [0045] [fig.5F] Les figures 5A à 5F montrent un exemple d'application de l'invention sur un échantillon comportant des cellules flottantes. [0046] [fig.6A] [0047] [fig.6B] 6 100481 Ifie.6C1 [0049] [fig.6D] Les figures 6A à 6D montrent un autre exemple d'application de l'invention sur un échantillon comportant des cellules flottantes [0050] [fig.7A] [0051] [fig.713] [0052] [fig.7C] Les figures 7A à 7C montrent un exemple d'application de l'invention sur un échantillon comportant des cellules adhérentes. EXPOSE DE MODES DE REALISATION PARTICULIERS [0053] La figure 1A représente un exemple de dispositif 1 selon l'invention. Une source de lumière 11 est apte à émettre une onde lumineuse 12, dite onde lumineuse incidente, se propageant en direction d'un échantillon 10, selon un axe de propagation Z. L'onde lumineuse est émise selon une bande spectrale Ak, comportant une longueur d'onde X. La longueur d'onde X peut être une longueur d'onde centrale de ladite bande spectrale. [0054] L'échantillon 10 est un échantillon que l'on souhaite caractériser. Il peut notamment s'agir d'un milieu 10,' comportant des particules 10p. Les particules 10p peuvent être des particules sanguines, par exemple des globules rouges. Il peut également s'agir de cellules, de microorganismes, par exemple des bactéries ou des levures, des microalgues, des microbillcs, ou des gouttelettes insolubles dans le milieu liquide, par exemple des nanoparticules lipidiques. De préférence, les particules 10' ont un diamètre, ou sont inscrites dans un diamètre, inférieur à 1 mm, et de préférence inférieure à 100 prn. Il s'agit de microparticules (diamètre inférieur à 1 mm) ou de nanoparticules (diamètre inférieur à 1 in). Le milieu 10', dans lequel baignent les particules, peut être un milieu liquide, par exemple une phase liquide d'un liquide corporel, d'un milieu de culture ou d'un liquide prélevé dans l'environnement ou dans un procédé industriel. Il peut également s'agir d'un milieu solide ou ayant la consistance d'un gel, par exemple un substrat de type gélose, propice à la croissance de colonies bactériennes. [0055] L'échantillon peut également être un échantillon solide, par exemple une fine lame de tissu biologique, de type lame d'anatomopathologie, ou un extrait sec d'un fluide, par exemple d'un fluide biologique. [0056] L'échantillon est de préférence transparent ou suffisamment translucide pour pouvoir peimettre la formation d'une image par le capteur d'image. [0057] L'échantillon 10 peut être contenu dans une chambre fluidique 16, par exemple une micro-cuvette, d'utilisation courante dans les dispositifs de type point of care, dans laquelle l'échantillon 10 pénètre, par exemple par capillarité. L'épaisseur e de l'échantillon 10, selon l'axe de propagation varie typiquement entre 20 im et 1 cm, et est de préférence comprise entre 50 Rm et 500 Rm, par exemple 150 !lm. [0058] L'échantillon s'étend selon un plan P io, dit plan de l'échantillon, perpendiculaire à 7 l'axe de propagation. Il est maintenu sur un support 10s. Le plan de l'échantillon est repéré par deux axes orthogonaux X et Y, définissant respectivement des coordonnées X et Y. Chaque couple de coordonnée (X, Y) correspond à une position radiale 1.. Les positions radiales sont définies dans le plan de l'échantillon et dans un plan de détection décrit ci-après. [0059] La distance D entre la source de lumière 11 et l'échantillon 10 est de préférence su- pé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. [0060] La source de lumière 11 peut être une diode électroluminescente ou une diode laser. Elle peut être associée à diaphragme 18, ou filtre spatial. L'ouverture du diaphragme est typiquement comprise entre 5 pin et 1 mm, de préférence entre 10 prn et 200 !lm ou 500 pm. Dans cet exemple, le diaphragme est fourni par Thorlabs sous la référence Pi 50S et son diamètre est de 150 pm. 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. [0061] Le dispositif peut comporter 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 la source de lumière élémentaire 11 selon un cône d'angle a, a étant égal à 30° dans le cas présent. De préférence, l'angle de diffusion a varie entre 100 et 80°. [0062] De préférence, la bande spectrale d'émission A?, de l'onde lumineuse incidente 12 a une largeur inférieure à 100 nm. Par largeur de bande spectrale, on entend une largeur à mi-hauteur de ladite bande spectrale. Dans la suite du texte, la bande spectrale est désignée par une longueur d'onde X représentative de la bande spectrale, et correspondant par exemple à la longueur d'onde centrale. [0063] L'échantillon 10 est disposé entre la source de lumière 11 et un capteur d'image 20. Le capteur d'image 20 définit un plan de détection Po. Ce dernier s'étend de préférence parallèlement, ou sensiblement parallèlement au plan P10 selon lequel s'étend 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 100 étant admise. [0064] Le capteur d'image 20 est apte à former une image selon le plan de détection P0. Dans l'exemple représenté, il s'agit d'un capteur d'image comportant une matrice de 8 pixels, de type 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. Le plan de détection Po s'étend de préférence perpendiculairement à l'axe de propagation Z de l'onde lumineuse incidente 12. Ainsi, le plan de détection Po est parallèle au plan de l'échantillon P10. Le capteur d'image comporte des pixels, à chaque pixel étant associée une position radiale 1-, dans le plan de détection Po. [0065] La distance d entre l'échantillon 10 et la matrice de pixels du capteur d'image 20 est préférentiellement comprise entre 50 prn et 2 cm, de préférence comprise entre 100 prn et 2 mm. [0066] Sur le dispositif représenté sur la figure 1A, on remarque l'absence d'optique de grossissement ou de formation d'image 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. [0067] 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 P0, 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 14, transmise par l'échantillon, et à laquelle est exposé le capteur d'image 20, comporte est formée suite à l'absorption et la diffraction de l'onde lumineuse incidente 12 par l'échantillon. Ainsi, l'échantillon entraîne une absorption d'une partie de l'onde lumineuse incidente, ainsi qu'un déphasage de cette dernière. Le déphasage est dû à une variation d'indice de réfraction (ou indice optique) lorsque la lumière se propage à travers l'échantillon. [0068] L'onde lumineuse 14 peut également être désignée par le terme onde lumineuse d'exposition. 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 programmable 31 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 32 permettant l'affichage d'images acquises par le capteur d'image 20 ou calculées par le processeur 30. [0069] L'image acquise par le capteur d'image forme un hologramme. Elle ne permet géné- ralement pas une représentation visuelle satisfaisante de l'échantillon, en particulier lorsque l'échantillon comporte des éléments diffractants rapprochés les uns des autres. C'est notamment le cas lorsque l'échantillon comporte des particules rapprochées les unes des autres, ou lorsque l'échantillon est une fine lame de tissu biologique. 9 [0070] La figure 1B schématise un autre dispositif 11 permettant une mise en oeuvre de l'invention. Contrairement au dispositif représenté sur la figure 1A, le dispositif de la figure 1B comporte un système optique 15 de formation d'image. Le système optique 15 définit un plan image P et un plan objet Pu. Le système optique peut être une lentille ou un objectif. Lors de l'acquisition de l'image de l'échantillon, le capteur d'image est disposé selon une configuration défocal iséc. Le plan de détection est décalé par rapport au plan image et/ou le plan selon lequel s'étend l'échantillon est décalé par rapport au plan objet. Le décalage à est généralement faible, en étant de préférence inférieur à 1 mm, et typiquement dans une plage 50 ni - 500 !lm. [0071] Quel que soit le dispositif utilisé, l'échantillon être décrit par des paramètres. chaque position radiale correspond un ou plusieurs paramètres. Les paramètres correspondant à une même position peuvent former un vecteur F(r), chaque vecteur étant défini à une position radiale r dans le plan de l'échantillon. Chaque terme de chaque vecteur correspond à une propriété optique de l'échantillon. Chaque position radiale correspond à un ou plusieurs pixels du capteur d'image. [0072] Chaque vecteur de paramètres F(r) est de dimension W. W est un entier strictement positif. W correspond au nombre de paramètres considérés à chaque position radiale. Chaque vecteur F0) comporte W termes Fw(r), de telle sorte que [0073] F1(r) F(r) = F,,(r) Fu(r) [0074] L'ensemble des vecteur F(r). pour les différentes positions radiales considérées, forme un jeu de paramètres Jg rassemblant les paramètres de l'échantillon. [0075] On se base sur l'exemple décrit dans la demande FR1859618, dans lequel l'échantillon peut être décrit par une absorbance a(r) (premier terme du vecteur) ainsi qu'une différence de trajet optique L ( r ) (deuxième terme du vecteur), ces propriétés étant susceptibles de varier en fonction de la bande spectrale d'illumination Ainsi, chaque position radiale r, l'échantillon peut être décrit par W = 2 paramètres différents 100761 F', _ 1( r) = F i(r)- (r) 100771 et 100781 F',_ 2(r ) = F2(r) = 1,(r). [0079] L'absorbance a (r) correspond à une capacité de l'échantillon à absorber tout ou partie de l'oncle lumineuse d'illumination. Lorsqu'une particule est considérée comme transparente, a (r)- O. [0080] Dans d'autres modélisations, les paramètres peuvent être l'indice optique, c'est-à-dire l'indice de réfraction, de l'échantillon, sachant qu'il peut s'agir d'une grandeur complexe. Ainsi, les paramètres peuvent comprendre la partie réelle de l'indice optique, et/ou la partie imaginaire de l'indice optique. [0081] La différence de trajet optique L ( r ) dépend de l'épaisseur e (r) de l'échantillon, parallèlement à l'axe de propagation de la lumière, ainsi que de la différence d'indice induite par l'échantillon. Par exemple, lorsque l'échantillon comporte des particules 10, baignant dans un milieu 10', chaque particule 10p induit une différence de trajet optique L(r) telle que : [0082] L ( = (n p - n m) X e (r ) 100831 ot e(r) est l'épaisseur de la particule à la position radiale r; 100841 et n r, et el en correspondent respectivement aux indices de réfraction de la particule 10p et du milieu 100851 Sur la figure 2, on a représenté, par des pointillés, le front d'onde de l'onde lumineuse 12, incidente à l'échantillon, ainsi que de l'onde lumineuse d'exposition 14, à laquelle est exposé le capteur d'image 20, l'onde lumineuse d'exposition se propageant de l'échantillon vers le capteur d'image. Dans cet exemple, p > m . Le front d'onde est plan avant d'atteindre l'échantillon. En aval de l'échantillon, le front d'onde est déformé du fait de l'apparition de différences de chemin optique induites par la particule. [0086] Soit A in. l'image de l'amplitude complexe de l'onde lumineuse d'exposition 14 dans le plan de l'échantillon P 10. Cette image, qui est une image complexe, peut également être considérée comme une image complexe de l'échantillon. A chaque position radiale , l'amplitude complexe A 10 ( r peut être définie à partir du vecteur de paramètres F (r ) correspondant à la position radiale r Lorsque le vecteur de paramètres F ( r ) comporte les termes a(r) et L ( r ) , l'amplitude complexe A 10 ( r ) peut être exprimée selon l'expression suivante : [0087] A 10( r ) b (r)exp (21 n ( r ) + ot (r)) (1) [0088] Le terme b (r) est une amplitude représentative de l'onde lumineuse incidente 12 at- teignant l'échantillon. Cette amplitude peut être mesurée par le capteur d'image 20 en l'absence d'échantillon 10 sur le support 10s. La source de lumière 11 éclaire alors directement le capteur d'image. A partir de l'image Jo, b(r) acquise par le capteur d'image 20 dans la bande spectrale d'émission a, À , on obtient par exemple l'amplitude b (r) selon l'expression : 100891 b (r) 0, b(r). (1') 100901 L'expression (I) permet de déterminer une expression d'une image complexe à partir des vecteurs de paramètres déterminés dans les différentes positions radiales. 11 [0091] Le terme 2in-L (r ) détermine la phase de l'amplitude complexe A 10 À apparaît que ce terme est X périodique. Cela signifie que différentes particules dont la différence de trajet optique est égale à L (r ) ci À , q étant un entier, permettent d'obtention d'une même amplitude complexe Alo(r ) entre l'échantillon et le capteur d'image. Autrement dit, il existe potentiellement une infinité d'objets susceptibles de générer une même amplitude complexe entre l'échantillon et le capteur d'image. [0092] Un des objectifs des algorithmes de reconstruction holographique est de reconstruire les propriétés optiques d'un objet à partir d'une image acquise par un capteur d'image, formant un hologramme. Or, l'hologramme ne comporte qu'une information partielle relative à l'onde lumineuse d'exposition, à laquelle est exposé le capteur d'image. En particulier, l'image acquise par le capteur d'image, c'est-à-dire l'hologramme, ne comporte pas d'information relative à la phase de l'onde lumineuse d'exposition. La mise en oeuvre d'algorithmes de reconstruction holographique itératifs permet d'estimer, de façon itérative, une information relative à la phase de l'onde lumineuse d'exposition, qui correspond au déphasage de l'onde lumineuse incidente induit par l'objet, Lorsque l'onde lumineuse d'exposition est exprimée selon l'expression (1), le terme relatif à la phase correspond au terme . 2 L ( r ) Un problème de 1 Tr l'estimation d'un terme de phase à l'aide des algorithmes actuels, usuellement désigné par le terme "phase unwrapping", est d'identifier la valeur correcte du terme de phase parmi les infinités de valeurs. X périodiques, qui résultent en une même expression de l'amplitude complexe de l'onde lumineuse d'exposition. [0093] Le procédé décrit ci-dessous, dont les principales étapes sont décrites sur la figure 3A, permet d'adresser cette question. [0094] Etapc 100 : Illumination de l'échantillon et acquisition d'une image /0 dans chaque bande spectrale d'illumination A À . [0095] Etape 110 : Initialisation, Au cours de cette étape, on prend en compte une image complexe d'initialisation A° dans le plan de l'échantillon P10, L'indice 10 du symbole A ° désigne le fait que l'image est acquise dans le plan de l'échantillon P10 L'exposant 0 du symbole A 0 désigne le fait qu'il s'agit d'une image d'initialisation. [0096] En considérant une onde lumineuse d'exposition telle que définie par l'expression (1), l'initialisation revient à considérer, pour chaque position radiale r, des paramètres de l'échantillon formant un vecteur initial de paramètres, tel que : [0097][a°(r) Fa ( r ) = (r ) J 100981 Les termes composant le vecteur initial peuvent être définis de façon arbitraire, ou en 12 fonction d'un a priori sur l'échantillon. Par exemple, on peut attribuer à chaque terme une même valeur, par exemple une valeur nulle. Les vecteurs de paramètres F0 r) définis lors de cette étape forment un ensemble d'initialisation g° décrivant l'échantillon 10. Ils forment également une image complexe initiale A 10C r) dans le plan de l'échantillon P10. [0099] Les étapes 120 à 160 décrites ci-après sont mises en oeuvre de façon itérative, selon un rang d'itération n. n est un entier compris entre 1 (première itération) et N, N correspondant au nombre d'itérations. Au cours de l'étape 110, n = 0. Dans les notations utilisées, le rang d'itération est présenté sous la forme d'un exposant. [0100] A chaque itération est associé un ensemble gn de vecteurs de paramètres Fn(r). Chaque vecteur de paramètres Fn(r), associé à une position radiale r , est, tel que : [0101]ia7(r)] Fn (r) = Ln (r) 101021 Chaque itération revient à mettre à jour les vecteurs de paramètres Fn( r), c'est-à-dire les termes aï7(r) et Ln ( r), et cela pour les différentes positions radiales 1". considérées. [0103] Etape 120 : Estimation d'une image dans le plan de détection. [0104] Pour chaque position radiale r , à partir des paramètres Fn -1 ( r) résultant de l'étape 110, ou de l'étape 150 d'une itération précédente, détermination d'une amplitude complexe -1(r) de l'onde lumineuse d'exposition 14, dans le plan de 10 l'échantillon P10. Lors de la première itération, " = I. Dans cet exemple, l'amplitude complexe est déterminée à partir de l'expression (1), à partir de l'ensemble initial de vecteurs de paramètres g° ou de l'ensemble de vecteurs de paramètres an -1 résultant d'une itération précédente. Les amplitudes complexes A nio-(r) définies pour les différentes positions radiales r considérées forment une image complexe A - I de 10 l'onde lumineuse d'exposition 14, dans le plan de l'échantillon P10. L'image complexe A A n - 1 est également appelée image complexe de l'échantillon. 10 [0105] Ainsi, l'image complexe- - 1 prise en compte dans l'étape 120 est élaborée à [0 partir d'un jeu de paramètres an - 1 initial (lorsque n = 1) ou résultant d'une itération précédente des étapes 120 à 150). [0106] A partir de l'image complexe An - on applique un opérateur de propagation holo- 10 graphique hP' de façon à obtenir une image complexe A n de l'onde o 0 lumineuse d'exposition 14, dans le plan de détection, selon l'expression : [0107] A- " - Ar] *h 0 lo -1'10 Po (2) 13 [0108] - est un opérateur de propagation holographique, permettant une pro- pagation du plan de l'échantillon P Io vers le plan de détection P0. Cet opérateur dépend de la longueur d'onde À . Il peut s'agir d'un opérateur de Fresnel, par exemple 101091 1 -127TI2 ± '2 h (x, y, z ) - iA z e exp TU À ) (3) [0110] avec [0111] r = Y) [0112] D'une façon générale, l'opérateur de propagation holographique h Po modélise - un transport de l'onde lumineuse d'exposition 14 entre au moins deux points distants l'un de l'autre. Dans l'application décrite, le produit de convolution décrit en lien avec l'équation (2) modélise un transport de l'onde lumineuse d'exposition 14 entre le plan de l'échantillon P 0 et le plan de détection Po. [0113] En considérant le carré du module de l'onde lumineuse d'exposition 14, on obtient une estimation An de l'image /0 acquise par le capteur d'image. Ainsi, 10, 101141 An n A " IO - O ri ( 04) [0115] A7) est l'image complexe conjuguée de l'image complexe A in). [0116] L'expression (4) revient à adopter un modèle de mesure simple, dans laquelle l'intensité de l'estimation de l'image acquise par le capteur correspond au carré du module de l'image complexe A n 0 101171 Selon une alternative, il est possible de prendre en compte la cohérence spatiale de la source de lumière, en considérant un noyau de convolution K, ainsi que des inhomogénéités d'éclairage, de telle sorte que : [0118] A n =B[ (ArgA(7)1') K I (4') [0119] Le noyau de convolution K traduit une surface de la source de lumière parallèlement au plan de détection. [0120] B peut être obtenu par calibration, par exemple par une acquisition en l'absence d'objet entre la source et le capteur d'image [0121] D'une façon générale, au cours de cette étape, l'estimation A" peut être décrite selon /0 l'expression : [0122] [0123] où A- I correspond à l'image complexe de l'onde lumineuse d'exposition 14 dans 10 le plan de l'échantillon, initiale (lorsque " = 1) ou résultant d'une itération précédente (lorsque n > 1), et ni est une fonction prenant en compte l'expression (4) (ou (4')) ainsi que l'expression (2). 14 [0124] L'image complexe de l'onde lumineuse d'exposition Io -1 dépend des paramètres contenus dans les vecteurs Fil - 1 ( r ) décrivant l'échantillon, en l'occurrence l'absorbance a n -1(r) et la différence de chemin optique n - -1(r), selon l'expression (1). [0125] Ainsi, on peut écrire : [0126] An 101271 Où a 1'1 - I et L n sont respectivement des images de l'absorbance et de la différence de chemin optique résultant de finitialisation ou d'une itération précédente. r est une fonction prenant en compte les expressions (1), (2) ainsi que (4) ou (4'). 101281 Etape 130 : comparaison de l'image A" estimée lors de l'étape 120 avec l'image -1-0 /0 acquise par le capteur d'image 20 lors de l'étape 100. La comparaison peut être exprimée sous forme d'une différence ou d'un ratio, ou d'un écart quadratique. [0129] Etape 140 : Calcul d'un indicateur de validité à partir de la comparaison effectuée, lors de l'étape 130, dans chaque bande spectrale. L'indicateur de validité re- présente la pertinence de l'ensemble an des vecteurs F(r) décrivant l'échantillon. L'indice tivn signifie que l'indicateur de validité est établi sachant l'ensemble an des vecteurs Fn(r). Dans cet exemple, l'indicateur de validité est d'autant plus faible que l'échantillon est correctement décrit par l'ensemble g. 101301 L'indicateur de validité er comporte un critère d'erreur ce dernier quantifiant une erreur globale de l'image estimée V' par rapport à l'image mesurée /a. Par erreur o globale, on entend une erreur pour chaque position radiale. [0131] Le critère d'erreur E_C yr: est établi à partir de la comparaison des images An et /0. Par /0 exemple, [0132] [0133] où : AT, est le nombre de positions radiales considérées ; est l'opérateur écart-type, moyennant une prise en compte d'un modèle de bruit. [0134] L'indice QI T-n; attribué au critère d'erreur EgiFir représente le fait que cet indicateur est calculé dans le plan de détection Po, sachant l'ensemble des vecteurs gn pris en compte lors de l'itération. [0135] Le critère d'erreur ,.iayet est un critère d'attache aux données, dans le sens ou sa mini- misation permet de se rapprocher des données mesurées, en l'occurrence de l'image 0. Aussi, lorsque tend vers /0, c'est-à-dire lorsque l'ensemble an des vecteurs 0 décrivent correctement l'échantillon 10, ÉE;Fss. tend vers 1. Au cours de l'étape 150, on peut appliquer un algorithme de minimisation. de type descente de gradient, de façon à s'approcher progressivement, à chaque itération, de l'ensemble gn permettant une minimisation satisfaisante de l'indicateur de validité er,n. Aussi, l'objectif de cette étape est d'établir un ensemble an de vecteurs F n (r ) visant à obtenir, suite à une ré- itération des étapes 110 à 140, un indicateur de validité plus faible que l'indicateur de validité ert: de l'intération courante [0136] Cette étape permet de mettre à jour au moins un terme F ) de chaque vecteur Fn ). [0137] Pour cela, on définit, pour chaque position radiale, un gradient Gn (r) de ^ l'indicateur de validité er,, à l'égard du paramètre optique correspondant au terme [0139] Un algorithme de descente de gradient définit alors une direction d n ainsi qu'un pas d'avancement Cnw. Le terme F Jr) de chaque vecteur de paramètre est mis à jour selon l'expression : [0140] F I ( r ) = F ( r ) +d, ci , (12) 101411 L'indicateur de validité ..E;rr, est une variable scalaire. Cependant, il dépend de l'ensemble gn des vecteurs paramètres sur la base desquels il a été établi, par l'intermédiaire de l'image Ain estimée lors de l'étape 120. 0 [0142] Le gradient Gn( r) peut être défini pour chaque terme 171)1, ( r ) des vecteurs J F n ( r ) considérés lors de l'itération de rang " [0143] Selon un premier mode de réalisation, l'indicateur de validité ne prend en compte que le critère d'erreur : E [0144] Dans une variante, détaillée par la suite, l'indicateur de validité comporte également un critère morphologique, correspondant à la prise en compte de contraintes géométriques ou optiques de l'échantillon ou des particules formant l'échantillon. [0145] Etape 150 : Mise à jour des paramètres de l'échantillon, formant les vecteurs (r ) par minimisation de l'indicateur de validité ett. Les paramètres sont mis à jour en ap- de telle sorte que : z.> [0138] 11) 16 pliquant l'expression (12). [0146] Etape 160: nouvelle itération des étapes 120 à 150, en prenant en compte, lors de l'étape 120 de l'itération suivante (n +1), l'ensemble an mis à jour lors de l'étape 150 de l'itération dernièrement effectuée. [0147] Les étapes 120 à 160 sont itérées jusqu'à cc que la valeur de l'indicateur de validité e(F., soit considérée comme représentative d'une bonne description de l'échantillon par l'ensemble g" de vecteurs En (r). N désigne le rang de la dernière itération. [0148] En prenant en compte un indicateur tel que défini dans les équations (10) et (13), les itérations cessent lorsque la valeur de l'indicateur de validité Ery's est suffisamment faible, ou lorsqu'un nombre prédéterminé d'itérations a été atteint, ou lorsque l'indicateur de validité n'évolue plus significativement entre deux itérations successives. [0149] Suite à la dernière itération, on dispose. de paramètres FA (r ) de l'échantillon. Cependant, lorsque des paramètres conditionnent la phase de l'onde lumineuse d'exposition 14, tel que la différence de trajet optique, une indétermination peut subsister, du fait de la périodicité précédemment décrite. 10150] La figure 3B schématise une telle indétermination. Sur cette figure, on suppose que l'on ne considère qu'un seul paramètre L, correspondant à une différence de trajet optique. L'axe des ordonnées représente un indicateur s;L. en fonction de valeurs du paramètre L (axe des abscisses), cette valeur correspondant à une position radiale. Le point A correspond à une valeur du paramètre L estimée lors d'une première itération, suite à une prise en compte d'une image initiale. La mise en oeuvre d'un algorithme de type descente de gradient conduit vers une minimisation de l'indicateur Ely- en direction du point B, de dernier correspondant à un minimum local. Or, le graphe montre que la valeur minimale de l'indicateur £;,,, qui correspond à la reconstruction optimale, est obtenue sur le point C. En partant du point A, le point C ne peut être atteint par un algorithme de minimisation de type descente de gradient. Il résulte de cc qui précède que la qualité de la reconstruction dépend de l'image initiale considérée. Lors de l'étape 110, on ne dispose d'aucun a priori sur l'échantillon. Aussi, même si les itérations des étapes 120 à 160 permettent d'obtenir une reconstruction de bonne qualité, l'atteinte d'une reconstruction optimale peut nécessiter une nouvelle initialisation, tenant compte d'un a priori sur l'échantillon. Une nouvelle initial isation peut permettre de démarrer l'algorithme de descente de gradient par un point D. Il devient alors possible d'atteindre le point C par descente de gradient. C'est l'objectif des étapes 170 à 180 décrites ci-après : [0151] Etape 170 : Actualisation de tout ou partie des paramètres de l'échantillon par un al- 17 gorithme d'intelligence artificielle à apprentissage supervisé. [0152] A l'issue de l'étape 160, on dispose de vecteurs de paramètres FN(r) définies pour chaque position radiale r , et formant un jeu de paramètres el. Certains termes de ces vecteurs, voire tous les termes de ces vecteurs, peuvent être utilisés en tant que données d'entrée d'un algorithme d'intelligence artificielle à apprentissage supervisé. Le jeu de paramètres gA' comporte R vecteurs FN r ) , chaque vecteur comportant W termes F.N. R désigne le nombre de positions radiales considérées. A chaque w position radiale r , certains de ces termes, voire tous les termes, peuvent former les données d'entrée de l'algorithme, comme décrit par la suite. [0153] L'algorithme d'intelligence artificielle peut être un réseau de neurones, par exemple un réseau de neurones convolutif, usuellement décrit par l'acronyme CNN (Convolutivc Neuron Network). [0154] Dans cet exemple, le réseau de neurones comporte deux couches d'entrée IN. Chaque couche d'entrée représente une distribution spatiale (ou image) d'un paramètre FAL décrivant l'échantillon, telle que mise à jour lors de la dernière itération des étapes 120 à 160 précédent l'étape 170. Dans cet exemple, la première couche d'entrée représente une distribution du premier paramètre Fx r r , en l'occurence l'absorbance a dans le plan de l'échantillon tandis que la deuxième couche d'entrée représente une distribution du deuxième paramètre F'2, r ), en l'occurence la différence de trajet optique TPV (r), dans le plan de l'échantillon. [0155] D'une façon générale, l'algorithme est mis en oeuvre à partir d'au moins une couche d'entrée IN, correspondant à une distribution spatiale d'un paramètre FN de rang -14' dans le plan de l'échantillon, résultant de la dernière itération n = N des étapes 120 à 160. Dans l'exemple considéré, on utilise, en guise de couches d'entrée, les deux distributions spatiales des paramètres F--\ et Fly résultant de la dernière itération N des étapes 120à 160. [0156] Entre les couches d'entrée IN et la couche de sortie OUT, le réseau de neurones comporte 20 couches Ll, L2.. .L20, dont les rangs sont compris entre 1 (couche adjacente de la couche IN) à 20 (couche adjacente de la couche OUT). Chaque couche comporte 20 plans. Une couche est obtenue par une convolution des 20 plans de la couche de rang précédent par un noyau de convolution de taille 3 x 3. La couche IN est considérée comme la couche de rang 0. Le réseau de neurones peut comporter une ou plusieurs couches de sortie OUT. Chaque couche de sortie représente une image d'un paramètre, dans le plan de l'échantillon. Dans cet exemple, le réseau de neurones ne comporte qu'une seule couche de sortie, correspondant à une image, dite image de 18 sortie, du deuxième paramètre F2. c'est-à-dire la différence de trajet optique L, dans le plan de l'échantillon. L'image de sortie comporte des paramètres, dits paramètres de sortie, actualisés par l'algorithme. [0157] La figure 3C schématise une architecture d'un tel réseau. [0158] De façon alternative, il est possible de recourir à d'autres architectures de réseaux de neurones, ou encore à des algorithmes d'intelligence artificielle à apprentissage supervisé, par exemple un algorithme de type SVM (Support Vector Machine). [0159] Suite à l'étape 170, on dispose, pour chaque position radiale : de paramètres F w(r) non actualisés par l'algorithme : il s'agit, dans cet exemple, des paramètres de rang w = 1, en l'occurrence l'absorbance a (r): ils correspondent aux paramètres résultant de paramètres F. ,(r) actualisés par l'algorithme CNN: il s'agit, dans cet exemple, des paramètres de rang w =2, en l'occurrence la différence de trajet optique L (r) 101601 L'indice ° représente le fait que les paramètres ont été actualisés par l'algorithme CNN. 101611 Selon un mode de réalisation, l'algorithme peimet une actualisation de l'ensemble des paramètres. 101621 Dans cet exemple, l'étape 170 permet de constituer des vecteurs F r ) formant un ensemble de paramètres g actualisés par l'algorithme. [0163] Le réseau de neurones convolutif CNN a préalablement fait l'objet d'un ap- prentissage, sur la base de données d'apprentissage maîtrisées, comme décrit par la suite (cf. étape 90). [0164] Etape 180 : Réitération des étapes 120 à 160 en prenant en compte les paramètres de sortie résultant de l'étape 170 pour former une image complexe d'initialisation. Au cours de cette étape, à partir des paramètres de l'échantillon résultant de l'étape 170, on met à jour l'image d'initialisation A 0, - Selon cet exemple, l'image d'initialisation rI 10 A 0 °, est établie, à partir des valeurs de a (r) et de L r ) , résultant de l'étape 170, selon l'expression (1). 101651 Les itérations des étapes 120 à 160 sont alors réitérées. Durant la première étape 140 de la réitération, on utilise l'image d'initialisation A 0. résultant de l'étape 180. Suite 2-1 aux itérations des étapes 120 à 160 : soit on considère qu'un recours à un réseau de neurones est nécessaire, auquel cas les étapes 170 et 180 sont renouvelées, et une nouvelle réitération des étapes 120 à 160 est effectuée, en utilisant l'image résultant de l'algorithme en tant qu'image d'initialisation ; 19 soit on estime que les paramètres g 7"; résultant de la dernière étape 150 correspondent à une bonne représentation de l'échantillon, auquel cas on sort de l'algorithme, ce qui correspond à l'étape de sortie d'algorithme 190. Les paramètres considérés comme représentatifs de l'échantillon sont alors ceux correspondant à la dernière itération N des étapes 120 à 160 : il s'agit de l'ensemble g \' des vecteurs EN résultant de la dernière itération. 101661 Etape 190: sortie d'algorithme. 101671 Au cours de cette étape, on forme une image d'un paramètre, ce qui correspond à une distribution spatiale du paramètre considéré. Cela permet d'avoir une représentation de l'échantillon. Au cours de cette étape, on peut former respectivement différentes images de différents paramètres. Dans l'exemple décrit, on peut former une image de de la différence de trajet optique et/ou une image de l'absorbance. Ces images peuvent permettre de caractériser l'échantillon. 101681 Ainsi, d'une façon générale, le procédé comporte : - une première série d'itérations des étapes 120 à 160, de façon à obtenir un ensemble g?'" de vecteurs F-f\! comportant les paramètres caractérisant l'échantillon ; - suite aux premières réitérations des étapes 120 à 160, un recours à l'algorithme d'intelligence artificielle à apprentissage supervisé, et par exemple au réseau de neurones convolutif précédemment décrit, en utilisant tout ou partie des paramètres caractérisant l'échantillon (étape 170). - une utilisation des paramètres de l'échantillon actualisés par l'algorithme d'intelligence artificielle à apprentissage supervisé, ces derniers étant utilisés en tant que paramètres d'initialisation pour une deuxième série d'itérations des étapes 120 à 160 (étape 180). - suite à la deuxième série d'itérations des étapes 120 à 160, l'algorithme d'intelligence artificielle à apprentissage supervisé peut être à nouveau mise en oeuvre, préalablement à une troisième série itérations des étapes 120 à 160 etc... [0169] Ainsi, on peut attribuer un rang x à chaque série d'itérations des étapes 120 à 160. Lors de la première série d'itérations, X = 1. La série d'itérations de rang x permet d'obtenir un ensemble de paramètres f . Tout ou partie des paramètres f ' X peuvent être actualisés par l'algorithme d'intelligence artificielle, de façon à obtenir des paramètres actualisés X, o. Ces derniers forment des paramètres d'initialisation utilisés dans une série d'itérations de rang X + 1. [0170] La dernière série d'itérations, de rang X, des étapes 120 à 160 permet d'obtenir des paramètres JN, X, ces derniers permettant d'obtenir au moins une image d'un paramètre de l'échantillon (étape 190). L'image du paramètre de l'échantillon, ou chaque image du paramètre de l'échantillon, peut permettre une caractérisation de l'échantillon. [0171] Le nombre de séries d'itérations des étapes 120 à 160 à effectuer peut être défini à priori, sur la base de calibrations à l'aide d'échantillons de calibration considérés comme comparables à l'échantillon analysé. Le nombre de séries d'itérations des étapes 120 à 160 peut aussi être établi au cas par cas, par exemple en comparant les paramètres résultant de deux séries différentes. Les étapes 170 et 180 sont répétées autant de fois que le nombre de séries d'itérations des étapes 120 à 160, moins 1. [0172] Un aspect intéressant de l'invention est que les paramètres résultant du réseau de neurones sont utilisés non pas pour obtenir une représentation finale d'échantillon, niais pour effectuer une initialisation d'un algorithme de reconstruction holographique. Cela permet d'utiliser les performances liées au recours à l'intelligence artificielle, tout en prenant en compte les données mesurées. [0173] L'algorithme de réseau de neurones convolutif mis en oeuvre dans l'étape 170 suppose une phase préalable d'apprentissage à l'aide d'échantillons connus, par exemple des échantillons numériques obtenus par simulations. Les échantillons d'apprentissage doivent de préférence être représentatifs des échantillons analysés par la suite. L'apprentissage fait l'objet de l'étape 90. [0174] Etape 90 : apprentissage. [0175] La figure 4A montre une partie d'un exemple d'échantillon numérique, obtenu en simulant des sphères. Les inventeurs ont simulé 4000 sphères s'étendant dans le champ d'observation du capteur d'image, c'est-à-dire selon environ 3 mm2. Ces sphères sont considérées comme représentatives de cellules. Les rayons des sphères sont choisis aléatoirement entre 4 'lm et 10 pm. Leur différence d'indice optique L (r), par rapport au milieu dans lequel baignent les sphères, est choisi aléatoirement entre 0.01 et 0.05. La répartition des sphères est aléatoire. La simulation permet de simuler des sphères regroupées. L'apprentissage est également effectué en prenant en compte une distance, entre l'échantillon et le capteur d'image, correspondant de préférence à la distance considérée lors des essais sur des échantillons réels. [0176] Les inventeurs ont simulé 1000 échantillons : ils ont ainsi établi des images d'absorbance et de différence de trajet optique sur les 1000 échantillons simulés. Chaque image s'étend sur 1000 x 1000 pixels. La figure 4A est un extrait d'une image de différence de trajet optique. L'échelle en niveaux de gris correspond à la valeur de la différence de trajet optique pour chaque pixel. [0177] A partir de chaque image, représentant une distribution spatiale des paramètres de l'échantillon dans le plan de l'échantillon, on a simulé des images acquises par le capteur d'image, en mettant en oeuvre un modèle, tel que décrit en lien avec 21 l'expression (5). La figure 4B correspond à une simulation d'image acquise par le capteur d'image (hologramme) à partir de la figure 4A. L'application du modèle suppose la prise en compte : d'une distance entre l'échantillon et le capteur d'image, en l'occurrence 1270 itm ; d'une bande spectrale en l'occurrence centrée sur 450 nm et de largeur spectrale égale à 15 nm ; d'un diamètre de la source de lumière, en l'occurrence 50 pm ; d'une distance entre l'échantillon et le capteur d'image, en l'occurrence 50 mm. [0178] A partir de chaque simulation d'image acquise par le capteur d'image, on a mis en oeuvre un algorithme itératif de reconstruction, tel que décrit en lien avec les étapes 120 à 160, de façon à obtenir des images complexes d'échantillons d'apprentissage, dans le plan d'échantillon. A partir des images complexes reconstruites, on a obtenu une image de chaque paramètre dans le plan de l'échantillon. Les figures 4C et 4D sont respectivement des images de l'absorbance or(r) et de la différence de trajet optique L(r) résultant d'une mise en oeuvre de l'algorithme itératif de reconstruction à partir de la figure 4B. [0179] En considérant les 1000 échantillons numériques simulés et, pour chaque échantillon simulé; une image reconstruite de l'absorbance et une image reconstruite de la différence de trajet optique, on a obtenu un total de 2000 images reconstruites, telles que celles représentées sur les figures 4C et 4D. A partir de ces images, on a extrait de manière aléatoire: 10000 vignettes de dimensions 121x121 représentant l'absorbance, dont un exemple est donné sur la figure 4E : il s'agit d'absorbances estimées par les itérations des étapes 120 à 160. 10000 vignettes de dimensions 121x121 représentant le trajet optique, dont un exemple est donné sur la figure 4F. il s'agit de différences de trajet optique estimées par les itérations des étapes 120 à 160. [0180] Ces vignettes ont été utilisées en tant que données d'entrée d'apprentissage du réseau de neurones convolutif. [0181] A partir des paramètres simulés, tels que représentés sur la figure 4A, on a extrait 10000 vignettes, de dimension 121x121, correspondant aux vraies différences de trajet optiques, ces vignettes étant utilisées en tant que données de sortie d'apprentissage du réseau de neurones. Un exemple est donné sur la figure 4G. [0182] Les figures 4E, 4F et 4G correspondent à une même partie d'un échantillon numérique. Elles forment un jeu d'apprentissage (2 entrées, 1 sortie). 10000 jeux d'apprentissage, tels que représentés sur les figures 4E à 4G, ont été formés. Les jeux d'apprentissage ont été utilisés pour paramétrer le réseau de neurones convolutif tel que décrit en lien avec la figure 3C. [0183] Variante [0184] Selon une variante, l'indicateur de validité comporte également un critère mor- phologique, correspondant à la prise en compte de contraintes géométriques ou optiques de l'échantillon ou des particules formant l'échantillon. l'indicateur de validité Er. prend également en compte un critère morphologique ere;ri. A la différence du critère d'erreur ebn, qui est défini à partir de données mesurées ou estimées dans le plan de détection Po, le critère morphologique eiltri est défini dans le plan de l'échantillon P10. 101851 D'une façon générale, le critère morphologique ef, Li-2n dépend de la valeur des termes des vecteurs de paramètres déterminés lors de l'étape 110 ou lors de l'étape 150 d'une itération précédente, ou de leurs dérivées spatiales. Il est représentatif de la morphologie de l'échantillon, telle que déterminée à partir des vecteurs de paramètres. Autrement dit, le critère morphologique représente un critère d'attache à des données morphologiques de l'échantillon, ces dernières pouvant être définies par des hypothèses. [0186] Le critère morphologique Ezen peut prendre en compte une dérivée spatiale de la différence de trajet optique, de façon à prendre en compte une forme prédéfinie d'une particule. Par exemple, lorsque l'échantillon comporte des cellules adhérentes, la forme prédéfinic peut être un hémisphère, un tel cas de figure étant représenté sur la figure 2. Lorsque l'échantillon comporte des cellules flottantes, la forme prédéfinic peut être une sphère lorsque les cellules sont sphériques. [0187] Par exemple, si l'amplitude complexe de ronde lumineuse d'exposition 14 est définie selon l'expression (1), chaque vecteur de paramètre comporte un terme L "(r) un exemple de critère morphologique est : [0188]1:r) \,2 (-0; )2 (1 017n = 2, 1 ay [0189] Ce critère tend à diminuer lorsque la grandeur L "(r) ne présente un minium d'oscillations, ce qui est par exemple le cas lorsque les particules ont une morphologie sphérique ou hémisphérique de particules. Les valeurs de L n(r) pour lesquelles le critère est minimal correspondent donc à des particules, par exemple sphériques ou hémisphériques. isolées les unes des autres, avec un minimum d'oscillation de L n(r) entre les particules ou sur ces dernières. [0190] Le critère morphologique Erc-Jyti est minimal lorsque les vecteurs de paramètres F"(r) formant l'ensemble gn décrivent des objets répondant à des hypothèses morphologiques établies à priori. 23 [0191] Lorsque l'indicateur dc validité e&a inclut une prise en compte du critère mor- phologique efol,pn, il peut être défini sous la forme d'une somme pondérée du critère d'erreur et du critère morphologique Erft T. L'expression de l'indicateur de validité peut alors être, par exemple : 101921 _ n = Yeroir-n [0193] où Y est un scalaire positif. [0194] Mise en oeuvre du procédé sur des données réelles. [0195] Le procédé décrit en lien avec les étapes 100 à 190 a été mis en oeuvre en utilisant un échantillon comportant cellules flottantes de type cellules ovariennes de hamster (CHO, acronyme de Chinese Hamster Ovary). Les cellules baignaient dans un milieu de culture, dans une chambre fluidique (l'épaisseur 20 pm. [0196] La source de lumière était une diode électroluminescente dont la bande spectrale d'émission était centrée sur 450 nm, et filtrée par un diaphragme de 50 pm. Le capteur d'image était un capteur monochrome comportant 3240 x 2748 pixels de côté 1.67 pm. La figure SA montre un hologramme obtenu, c'est-à-dire une image acquise par le capteur d'image. La figure 5B montre une image du deuxième paramètre L (r) résultant dc plusieurs itérations des étapes 120 à 160. Ces itérations permettent de définir une image des paramètres a (r) et L (r ). [0197] La figure SC montre une image obtenue en appliquant un réseau de neurones à partir des images des paramètres a (r) et L (r). La figure 5C montre une image de sortie du réseau de neurones, correspondant à une image du paramètre L (r). La com- paraison entre les figures 5B et 5C montre que le réseau de neurones effectue des ajustements très significatifs, certains d'entre eux étant entourés par un cercle en pointillés. On a constaté que les ajustements les plus significatifs sont réalisés suite aux premières itérations des étapes 120 à 160. On constate par exemple que des pixels sombres de la figure 5B apparaissent sous la forme de pixels clairs sur la figure 5C. [0198] L'image L (r) de la figure 5C ainsi que l'image du paramètre a (r) ont été utilisées pour établir une image complexe d'initialisation A selon l'expression (1). f10 En se basant sur cette image d'initialisation, les étapes de reconstruction itérative, décrite en lien avec les étapes 120 à 160, ont été réitérées. [0199] La figure 5D montre une image du paramètre L (r) à la suite d'une réitération des étapes 120 à 160. On observe des modifications par rapport à l'image du paramètre ' L (r) représentée sur la figure 5C, issue de la mise en oeuvre du réseau de neurones : la reconstruction permet d'identifier des objets non détectés sur les images précédentes (tirets en pointillés). Elle permet également de retirer certains objets détectés par le 24 réseau de neurones, ces derniers étant repérés (flèches blanches). Le fait d'appliquer une deuxième reconstruction itérative, sur la base de l'image de la figure 5C, permet d'affiner la précision avec laquelle l'image du paramètre I. (r) est établie. [0200] La réitération des étapes 120 à 160 a permis de définir des deuxièmes images des pa- ramètres a (r) et L (r). Ces images ont été utilisées en tant qu'images d'entrée du réseau de neurones convolutif. La figure 5E montre l'image de sortie du réseau de neurones I (r). représentant une actualisation de l'image du paramètre L (t). Les nouveaux objets, issus de la reconstruction itérative, sont confirmés sur cette image (traits en pointillés). L'image du paramètre is °(r) de la figure 5E ainsi que l'image du paramètre a (r) ont été utilisées pour établir une image complexe d'initialisation A 0, selon l'expression (1). En se basant sur cette image d'initialisation, les étapes de n'0 reconstruction itérative, décrite en lien avec les étapes 120 à 160, ont été réitérées une nouvelle fois. La figure 5F montre une image du paramètre L(r) à la suite de cette nouvelle réitération des étapes 120 à 160. [0201] Les figures GA à 6D représentent une deuxième série d'essais sur des cellules flottantes de type CHO (Chinese Hamster Ovary), selon une densité différente que lors de l'essai décrit en lien avec les étapes 5A à 5F. La figure 6A montre une partie d'une image acquise par le capteur d'image La figure 6B représente une image de la différence de trajet optique L(r) estimée après des premières itérations des étapes 120 à 160. On observe des repliements de phases, qui correspondent à des niveaux de gris sombres au niveau de certaines cellules. Les images de l'absorbance et de la différence de chemin optique, résultant de la dernière itération des étapes 120 à 160, ont été utilisées en tant qu'images d'entrée du réseau de neurone convolutif (étape 170). L'image de la figure 6C correspond à une image de la différence de trajet optique résultant du réseau de neurones : en comparant les figures 6B et 6C, on peut observer que le phénomène de repliement de phase a été corrigé par le réseau de neurones. L'image de l'absorbance, résultant des étapes 120 à 160, ainsi que l'image de la différence de trajet optique représentée sur la figure GC ont été utilisées pour former une image complexe d'initialisation. Cette dernière a été utilisée dans une nouvelle mise en oeuvre des étapes 120 à 160. La figure 6D montre une image de la différence de trajet optique résultant de cette nouvelle itération. 102021 Au cours d'une troisième série d'essais, on a utilisé un échantillon comportant des cellules adhérentes de type PC12 dans un milieu de culture. Un autre réseau de neurones a été utilisé, dont la phase d'apprentissage était basée sur des modélisations d'échantillons comportant des cellules adhérentes. 102031 La source utilisée était une LED émettant dans une bande spectrale centrée sur 605 nm, et de largeur spectrale égale à 20 nm. L'échantillon était disposé à une distance de 3500 p.m. La figure 7A montre une image de la différence de chemin optique L (r) obtenue suite à des premières itérations d'étapes 120 à 160. Les premiers paramètre a (") et les deuxièmes paramètres .L (r) résultant obtenus suite aux premières itérations des étapes 120 à 160 ont été utilisés en tant que données d'entrée du réseau de neurones convolutif. La figure 7B représente l'image de sortie du deuxième paramètre L (r), telle que fournie par le réseau de neurones. La comparaison des images 7A et 7B permet d'apprécier l'effet du réseau de neurones convolutif. Les deuxièmes paramètres actualisés ont été utilisés pour former une image complexe d'initialisation, cette dernière ayant été utilisée pour une deuxième mise en oeuvre des étapes 120 à 160. La figure 7C montre l'image résultant de la reconstruction. Cette image montre une évolution spatiale plus continue du paramètre L (r) sans les discontinuités, apparaissant sur la figure 7A, représentatives du phénomène de repliement de phase. 102041 Le procédé décrit ci-avant peut permettre de caractériser l'échantillon, à partir des pa- ramètres déterminés suite à l'étape 180. Par caractériser, il est notamment entendu, à titre non exhaustif: l'obtention d'une image de l'échantillon permettant une observation de ce dernier ; une détermination d'un nombre de particules composant l'échantillon ; une classification de particules en classes de particules ; un dénombrement de particules ; une identification de particules par rapport à des paramètres d'identification ; un dénombrement de particules identifiées ; une estimation de caractéristiques géométriques de l'échantillon ; une estimation de caractéristiques géométriques de particules composant l'échantillon; une détermination de la viabilité de l'échantillon ; une estimation d'une teneur en masse sèche dans l'échantillon, et par exemple d'une ou de plusieurs particules de l'échantillon. [0205] L'invention pourra être mise en oeuvre, sur des échantillons biologiques, dans le domaine de la santé, par exemple dans l'aide au diagnostic, ou dans l'étude de processus biologiques. Elle peut également être mise en oeuvre par des échantillons prélevés dans l'environnnement, ou dans des installations industrielles, par exemple dans le domaine de l'agroalimentaire. 26 [Revendication 1] REVENDICATIONSProcédé d'observation d'un échantillon (10), l'échantillon s'étendant selon un plan d'échantillon (P10) définissant des positions radiales, des paramètres ((L ( r ) , a ( r ), F1v( r ) ) de l'échantillon étant définis à chaque position radiale (r ), le procédé comportant : a. illumination de l'échantillon à l'aide d'une source de lumière (11), apte à émettre une onde lumineuse incidente (12) se propageant vers l'échantillon ; h. acquisition, à l'aide d'un capteur d'image (20), d'une image ( /0) de l'échantillon (10), formée dans un plan de détection (P0 ), l'échantillon étant disposé entre la source de lumière (11) et le capteur d'image (20) ; c. à partir de l'image acquise par le capteur d'image, obtention d'une image de l'échantillon, [image de l'échantillon correspondant à une distribution d'au moins un paramètre de l'échantillon dans le plan d'échantillon ; le procédé étant caractérisé en ce que l'étape c) comporte les étapes suivantes: i. prise en compte une image d'initialisation ii. application d'un opérateur de propagation holographique (h) à l'image d'initialisation, ou à une image complexe de [échantillon résultant de l'étape (y) d'une itération précédente, de façon à estimer l'image acquise par le capteur d'image ("n ) /0 dans le plan de détection ; comparaison de l'image (10) acquise lors de l'étape b) et de l'image estimée lors de l'étape ii) ); /0 - iv. en fonction de la comparaison, mise à jour des paramètres décrivant l'échantillon (F (r), ( r ), L n (r)); v. tant qu'un critère de convergence ou un nombre prédéterminé d'itérations n'a pas été atteint, mise à jour dune image complexe de l'échantillon (n n ), dans le plan de l'échantillon (P10), à partir des paramètres mis à jour lors de l'étape iv), puis réitération des étapes (ii) à (y) ; vi. utilisation de paramètres de l'échantillon, résultant de la 27 dernière étape (i v) effectuée, pour alimenter un algorithme d'intelligence artificielle à apprentissage supervisé, de façon à actualiser des paramètres de l'échantillon (F ( r ) , L ( r ) ); vii. utilisation des paramètres résultant de l'étape (vi) pour mettre à jour l'image d'initial isation (A ); viii. réitération, au moins une fois, des étapes (ii) à (v), en utilisant, lors de la première réitération, l'image d'initialisation (A o ) 10 mise à jour lors de l'étape (vii); ix. obtention d'une image de l'échantillon à partir des paramètres mis à jour lors l'étape iv) de la dernière itération des étapes (fi) à (v). [Revendication 2] Procédé selon la revendication 1, dans lequel l'étape iv) comporte un calcul d'un indicateur de validité (f.C,n), de telle sorte que les paramètres de l'échantillon sont mis à jour de façon à faire tendre [indicateur de validité vers une valeur prédéterminée. [Revendication 3] Procédé selon la revendication 2, dans lequel lors de l'étape iv), les pa- ramètres de l'échantillon sont mis à jour de façon à minimiser l'indicateur de validité. [Revendication 4] Procédé selon l'une quelconque des revendications précédentes, dans lequel l'algorithme d'intelligence artificielle à apprentissage supervisé met en oeuvre un réseau de neurones. [Revendication 5] Procédé selon la revendication 4, dans lequel le réseau de neurones est un réseau de neurones convolutif. [Revendication 6] Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape iv) comporte une détermination d'un gradient (G "v ( r ) ) } de l'indicateur de validité (j..n) en fonction d'au moins un paramètre ( r ) ), de telle sorte que les paramètres (Fn » sont mis à jour pour réduire l'indicateur de validité de l'itération suivante (E [Revendication 7] Procédé selon la revendication 6, dans lequel l'étape iv) met en oeuvre un algorithme de type descente de gradient. [Revendication 81 Procédé selon l'une quelconque des revendications précédentes, dans lequel : l'étape iv) résulte en la détermination d'une image de chaque 28 paramètre dans le plan de l'échantillon ; lors de l'étape vi), au moins une image d'un paramètre forme une couche d'entrée de l'algorithme d'intelligence artificielle à apprentissage supervisé ; lors de l'étape viii), l'algorithme d'intelligence artificielle à apprentissage supervisé délivre une image de sortie, correspondant à une image d'un paramètre de l'échantillon actualisé par l'algorithme. [Revendication 9] Procédé selon l'une quelconque des revendications précédentes, dans [Revendication 10] lequel lors de l'étape (iii), le calcul de l'image de l'échantillon ( A17) dans Io le plan de détection (Po) comporte l'application d'un produit de convolution par un noyau de convolution (K), le noyau de convolution représentant une étendue spatiale de la source de lumière. Procédé selon l'une quelconque des revendications précédentes, dans lequel les paramètres décrivant l'échantillon comportent : - un premier paramètre (L ( r ) ) représentant une absorbance de l'échantillon ; - un deuxième paramètre (a ( r ) ), représentant une différence de trajet optique, selon l'axe de propagation (Z) de l'onde lumineuse incidente. [Revendication 11] Procédé selon la revendication 10, dans lequel lors de l'étape (vi), l'algorithme d'intelligence artificielle à apprentissage supervisé est alimenté par : une image du premier paramètre, correspondant à une distribution spatiale du premier paramètre, mise à jour lors de la dernière itération des étapes (ii) à (v) précédent l'étape (vi); une image du deuxième paramètre, correspondant à une distribution spatiale du deuxième paramètre, mise à jour lors de la dernière itération des étapes (ii) à (v) précédent l'étape (vi) ; et dans lequel l'algorithme d'intelligence artificielle à apprentissage supervisé permet d'obtenir une image du deuxième paramètre actualisé. [Revendication 12] Procédé selon l'une quelconque des revendications précédentes, dans 29 [Revendication 13] lequel suite à une étape (viii), les étapes (vi) à (viii) sont renouvelées au moins une fois. [Revendication 14] Procédé selon l'une quelconque des revendications précédentes, dans lequel aucune optique de formation d'image n'est disposée entre l'échantillon et le capteur d'image. Procédé selon l'une quelconque des revendications 1 à 12 dans lequel un système optique, de type lentille ou objectif, est disposé entre l'échantillon et le capteur d'image, le système optique définissant un plan image (P,) et un plan objet (Pu), le procédé étant tel que lors de l'étape h) : le plan objet (P0) est décalé du plan d'échantillon (P10) selon une distance de défocalisation objet ; et/ou le plan image (P1) est décalé du plan de détection (P0) selon une distance de défocalisation image. [Revendication 15] Dispositif d'observation d'un échantillon (10) comportant : une source de lumière (11), configurée pour émettre une onde lumineuse incidente pour illuminer l'échantillon ; un support d'échantillon (10s), configure pour recevoir l'échantillon ; un capteur d'image (20), configure pour acquérir une image dc l'échantillon lorsque l'échantillon est disposé sur le support d'échantillon (10s) ; un processeur (30), programmé pour exécuter des instructions permettant de mettre en oeuvre les étapes i) à viii) d'un procédé selon l'une quelconque des revendications précédentes à partir d'une image acquise par le capteur d'image.