Brevet FR3087009A1 — publié le 20200410

Procédé de détermination de paramètres d'une particule Description DOMAINE TECHNIQUE Le domaine technique de l'invention est la caractérisation de particules, par exemple des particules biologiques, à partir d'une image. ART ANTERIEUR L'imagerie holographique connaît, depuis plusieurs années, des développements importants, en particulier dans le domaine de l'analyse d'échantillons. Un grand nombre d'applications ont été décrites, pour caractériser des échantillons, en particulier des échantillons destinés à des applications en biologie ou pour du diagnostic. De nombreuses publications concernent par exemple des échantillons comportant des particules, l'objectif étant de caractériser les particules. De façon usuelle, on forme une image d'un échantillon, désignée par le terme hologramme, et un opérateur de propagation holographique est appliqué à l'image. L'image de l'échantillon peut être acquise selon une configuration dite défocalisée, comme décrit dans le document W02016/097092. Dans ce document, on décrit un procédé d'imagerie selon une configuration défocalisée, pour identifier un microorganisme, à partir d'une image d'un échantillon acquise par un capteur d'image décalé par rapport à un plan de mise au point d'un système optique. Un opérateur de propagation holographique est alors appliqué à l'image défocalisée. Selon une autre configuration, l'image de l'échantillon peut être acquise selon une configuration d'imagerie sans lentille, aucune optique de formation d'image n'étant disposée entre l'échantillon et le capteur d'image. Le document W02016/151249 décrit un procédé d'analyse de cellules, disposées dans un milieu de culture, sans marquage. Le document W02016/151248 décrit un procédé pour identifier des particules, par exemple des particules sanguines. Les procédés décrits dans ces documents sont basés sur l'établissement d'un profil, représentant une évolution, le long d'un axe optique, d'une grandeur caractéristique d'une onde lumineuse d'exposition, à laquelle est exposé le capteur d'image. Généralement, à chaque particule de l'échantillon est associé un profil. La particule est caractérisée en comparant le profil à des profils étalons, ces derniers étant obtenus de façon expérimentale, en mettant en oeuvre des échantillons dont les particules ont des caractéristiques connues. 2 Les documents précédemment évoqués sont basés sur une analogie entre un profil établi en appliquant un opérateur de reconstruction holographique à une image d'échantillon, et une série de profils étalons, obtenus expérimentalement. Afin d'obtenir une reconstruction holographique de bonne qualité, l'image de l'échantillon peut faire l'objet d'un algorithme de 5 reconstruction itératif. Des algorithmes de reconstruction itératifs sont par exemple décrits dans W02016189257 ou dans W02017162985. Les procédés préalablement décrits supposent que les profils mesurés à partir de l'échantillon et les profils étalons soient préférentiellement établis en mettant en oeuvre le même algorithme de reconstruction holographique. De plus, il est préférable, sans pour autant que cela soit nécessaire, que les profils étalons soient obtenus dans 10 des conditions expérimentales les plus proches possibles des conditions dans lesquelles l'image de l'échantillon est formée. Les conditions expérimentales sont par exemple le type de source de lumière, le type de capteur, la taille des pixels, la distance entre le capteur et l'échantillon. Les inventeurs ont cherché à perfectionner les procédés décrits ci-avant, de façon à être moins liés aux conditions expérimentales et à l'algorithme mis en oeuvre pour établir les profils étalons. 15 EXPOSE DE L'INVENTION Un premier objet de l'invention est un procédé de détermination d'un paramètre d'une particule présente dans un échantillon, l'échantillon s'étendant entre un capteur d'image et une source de lumière, le capteur d'image s'étendant selon un plan de détection, le procédé comportant les étapes suivantes : 20 a) illumination de l'échantillon par la source de lumière, la source de lumière émettant une onde lumineuse incidente se propageant selon un axe de propagation; b) acquisition d'une image de l'échantillon par le capteur d'image, le capteur d'image étant exposé à une onde lumineuse d'exposition; c) détermination d'une position de la particule dans le plan de détection ou dans un plan 25 parallèle au plan de détection ; d) à partir de l'image acquise, application d'un opérateur de propagation, selon une pluralité de distances du plan de détection, de façon à estimer, à chaque distance, une amplitude complexe de l'onde lumineuse d'exposition ; e) à partir de l'amplitude complexe estimée, en différentes distances, obtention d'un profil 30 représentant une évolution d'une propriété de l'onde lumineuse d'exposition selon un axe parallèle à l'axe de propagation et passant par la position de la particule. La particule peut être associée à un ensemble de paramètres, comportant au moins une taille de la particule et un indice de réfraction de la particule. 3 Le procédé peut comporter : f) une modélisation d'une particule, prenant en compte une valeur de chaque paramètre de l'ensemble de paramètres, et une modélisation d'une onde lumineuse d'exposition, se propageant vers le capteur d'image, sous l'effet d'une illumination, par la source de 5 lumière de la particule modélisée ; g) à partir de l'onde lumineuse d'exposition modélisée, formation d'un profil représentant une évolution de l'amplitude complexe de l'onde lumineuse d'exposition modélisée, selon un axe parallèle à l'axe de propagation; h) comparaison entre le profil obtenu lors de l'étape e) et le profil modélisé résultant de 10 l'étape g), de façon à déterminer la valeur d'au moins un paramètre de la particule. Selon un mode de réalisation, l'indice de réfraction est une grandeur complexe : il comporte une partie réelle et une partie imaginaire. Selon un mode de réalisation, l'ensemble de paramètres comporte également une position de la particule selon l'axe de propagation. 15 Selon un mode de réalisation, l'étape f) comporte des modélisations de particules présentant respectivement différentes valeurs d'au moins un paramètre, de façon à obtenir, suite à l'étape g), une base de données de profils modélisés, à chaque profil modélisé étant associé un ensemble de paramètres. L'étape h) peut comporter une minimisation d'un écart entre le profil résultant de l'étape e) et les profils de la base de données, les valeurs respectives des paramètres 20 de la particule étant celles minimisant l'écart. L'écart peut être un écart quadratique. Selon un mode de réalisation, les étapes f) à h) sont effectuées de façon itérative, la valeur d'un paramètre de la particule résultant d'une itération étant utilisée pour initialiser une itération suivante. Le paramètre peut être l'indice de réfraction, et/ou la taille et/ou une distance par rapport au plan de détection. 25 Selon un mode de réalisation, l'étape f) comporte des modélisations de particules de différentes tailles et/ou de différents indices de réfraction, de façon à obtenir, suite à l'étape g), une base de données de profils modélisés, à chaque profil modélisé étant associé une taille et/ou un indice de réfraction. L'étape h) peut comporter une détermination de la taille de la particule et/ou de l'indice de réfraction de la particule. 30 Selon un mode de réalisation, l'étape f) comporte des modélisations de particules situées à différentes distances du plan de détection, de façon à obtenir, suite à l'étape g) une base de données de profils modélisés, chaque profil modélisé étant associé à une distance par rapport au plan de détection. L'étape h) peut comporter une détermination de la distance entre la particule et le plan de détection. Selon un mode de réalisation, les étapes f) à h) sont mises en oeuvre de façon itérative, de telle sorte qu'au cours de chaque itération, le profil modélisé lors de l'étape g) se rapproche 5 progressivement du profil déterminé lors de l'étape e). Selon ce mode de réalisation, l'étape h) d'une itération peut comporter une détermination d'un écart entre le profil modélisé lors de l'étape g) de la même itération et le profil résultant de l'étape e), ainsi qu'une détermination d'un gradient de l'écart en fonction d'au moins un paramètre de l'ensemble de paramètres, de façon à déterminer les valeurs des paramètres de la particule modélisée lors de l'étape f) de 10 l'itération suivante. Quel que soit le mode de réalisation, le procédé peut comporter : une application d'un opérateur de propagation à l'image acquise, en considérant plusieurs distances par rapport au plan de détection, de façon à obtenir une pile d'images reconstruites, à chaque image reconstruite étant associée une distance par 15 rapport au plan de détection ; à partir d'une image reconstruite, une estimation d'une taille de la particule ; de telle sorte que lors de l'étape h) comporte une prise en compte de la taille de la particule estimée à partir de l'image reconstruite. L'image reconstruite, sur la base de laquelle l'estimation de la taille de la particule est effectuée, est de préférence une image reconstruite dans un plan 20 de reconstruction passant par la particule. Une telle image forme une image d'observation de la particule. La position du plan de reconstruction correspondant à l'image d'observation peut être déterminée par la mise en oeuvre d'un algorithme de focalisation numérique. Quel que soit le mode réalisation, l'étape h) peut comprendre une détermination d'un profil modélisé le plus proche du profil résultant de l'étape e). Les valeurs des paramètres de la 25 particule correspondent alors respectivement aux valeurs des paramètres associés au profil modélisé le plus proche, c'est-à-dire aux valeurs des paramètres de la particule modélisée ayant permis d'obtenir le profil modélisé le plus proche. La particule peut notamment être une cellule, le procédé comportant une détermination d'un état de vie ou de mort de la cellule en fonction de l'indice de réfraction. De préférence, l'indice 30 de réfraction comporte une partie réelle et une partie imaginaire, l'état de vie ou de mort de la cellule étant défini en fonction d'une comparaison entre la partie réelle et la partie imaginaire. 5 Selon un mode de réalisation, il n'y a pas d'optique de formation d'image entre l'échantillon et le capteur d'image. Selon un autre mode de réalisation, un système optique s'étend entre l'échantillon et le capteur d'image, le système optique définissant un plan image et un plan objet, et dans lequel lors de l'étape b), l'image est acquise selon une configuration défocalisée, 5 le plan de détection étant décalé par rapport au plan image, et/ou un plan de l'échantillon, selon lequel s'étend l'échantillon, étant décalé par rapport au plan objet. Un autre objet de l'invention est un dispositif d'observation d'un échantillon, l'échantillon comportant des particules, le dispositif comportant : - une source de lumière, configurée pour illuminer l'échantillon ; 10 - un support d'échantillon, configuré pour recevoir l'échantillon ; - un capteur d'image, disposé de telle sorte que lorsque l'échantillon est disposé sur le support, il s'étend entre le capteur d'image et la source de lumière, le capteur d'image étant configuré pour acquérir une image de l'échantillon ; - un processeur, configuré pour mettre en oeuvre les étapes d) à h) d'un procédé selon le 15 premier objet de l'invention, à partir d'une image acquise par le capteur d'image. Selon un mode de réalisation, aucune optique de formation d'image n'est disposée entre le capteur d'image et le processeur. Selon un autre mode de réalisation, le dispositif comporte un système optique s'étendant entre l'échantillon et le capteur d'image, le système optique définissant un plan image et un plan objet, le dispositif comportant un moyen d'ajustement du 20 système optique, ou de l'échantillon, ou du capteur d'image, de telle sorte que : - le plan de détection soit décalé du plan image ; - et/ou le plan objet soit décalé d'un plan d'échantillon, selon lequel s'étend l'échantillon. D'autres avantages et caractéristiques ressortiront plus clairement de la description qui va suivre de modes particuliers de réalisation de l'invention, donnés à titre d'exemples non limitatifs, et 25 représentés sur les figures listées ci-dessous. FIGURES La figure 1 montre un dispositif permettant une mise en oeuvre de l'invention. Les figures 2A à 2F sont des résultats de modélisation d'une figure de diffraction ou de profils correspondant à des cellules dont les paramètres sont connus. Il en est de même des figures 3A 30 à 3F. Les figures 4A à 4D schématisent les principales étapes de différents modes de réalisation. Les figures 5A à 5C illustrent une détermination de paramètres d'une cellule. 6 Les figures 6A et 6B montrent l'évolution d'une erreur de quantification affectant respectivement les estimations de la partie réelle et de la partie imaginaire de l'indice de réfraction en fonction de la taille de particules. La figure 7 montre un exemple d'application de l'invention à la classification de cellules vivantes 5 et de cellules mortes. La figure 8 montre un autre dispositif permettant une mise en oeuvre de l'invention. EXPOSE DE MODES DE REALISATION PARTICULIERS La figure 1 représente un exemple de dispositif selon l'invention. Une source de lumière 11 est 10 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 d'illumination A2L. L'échantillon 10 est un échantillon que l'on souhaite caractériser. Il comprend notamment un milieu 10m dans lequel baignent des particules 10t. Le milieu 10m peut être un milieu liquide. Il 15 peut comprendre un liquide corporel, obtenu par exemple à partir de sang ou d'urine ou de lymphe ou de liquide céphalorachidien. Il peut également s'agir d'un milieu de culture, comportant des nutriments permettant le développement de microorganismes ou de cellules. Par particule, on entend notamment, et de façon non exhaustive : - une cellule, qu'il s'agisse d'une cellule de culture ou une cellule corporelle, par exemple 20 une cellule sanguine ; - un microorganisme, par exemple une bactérie ou une levure ou une microalgue ; - une particule solide, par exemple une microbille, la microbille pouvant être fonctionnalisée, de façon à favoriser un greffage avec un analyte ; - une particule formant une émulsion dans le milieu 10m, en particulier une particule 25 insoluble dans le milieu 10m, un exemple étant une gouttelette lipidique dans un milieu aqueux. Une particule 10i peut être solide ou liquide. L'échantillon 10 est, dans cet exemple, contenu dans une chambre fluidique 15. La chambre fluidique 15 est par exemple une chambre fluidique de type Gene Frame® d'épaisseur e=250 30 p.m. L'épaisseur e de l'échantillon 10, selon l'axe de propagation, varie typiquement entre 10 p.m et 1 cm, et est de préférence comprise entre 20 p.m et 500 wn. L'échantillon s'étend selon un plan P10, dit plan de l'échantillon, perpendiculaire à l'axe de propagation Z. Le plan de 7 l'échantillon est défini par les axes X et Y représentés sur la figure 1. L'échantillon est maintenu sur un support 10s à une distance d d'un capteur d'image 16. La concentration de particules peut varier entre 50 par microlitre et 500 000 par microlitre. Elle peut être par exemple égale à 100000 par microlitre. 5 La distance D entre la source de lumière 11 et la chambre fluidique 15 est de préférence supérieure à 1 cm. Elle est de préférence comprise entre 2 et 30 cm. Avantageusement, la source de lumière 11, 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 la chambre fluidique 15 et la source de lumière. Sur la figure 1, la source de 10 lumière est une diode électroluminescente. Elle est généralement associée à un diaphragme 18, ou filtre spatial. L'ouverture du diaphragme est typiquement comprise entre 5 pm et 1 mm, de préférence entre 50 pm et 500 pm. Dans cet exemple, le diaphragme est fourni par Thorlabs sous la référence P1505 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 15 une deuxième extrémité est placée en regard de l'échantillon 10. Le dispositif représenté sur la figure 1 comporte également 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 11 20 selon un cône d'angle a. De préférence, l'angle de diffusion a varie entre 100 et 800 . Alternativement, la source de lumière peut être une source laser, telle une diode laser. Dans ce cas, il n'est pas utile de lui associer un filtre spatial ou un diffuseur. De préférence, la bande spectrale d'émission 4X 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 25 ladite bande spectrale. L'échantillon 10 est disposé entre la source de lumière 11 et le capteur d'image 16 précédemment évoqué. Le capteur d'image 16 définit un plan de détection P0, s'étendant 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 30 pas être rigoureusement parallèles, une tolérance angulaire de quelques degrés, inférieure à 20 ou 10° étant admise 8 Le capteur d'image 16 est apte à former une image /0 de l'échantillon 10 selon le plan de détection P0. Dans l'exemple représenté, il s'agit d'un capteur d'image 16 comportant une matrice de pixels, de type CCD ou un CMOS. Le plan de détection P0 s'étend de préférence perpendiculairement à l'axe de propagation Z de l'onde lumineuse incidente 12. La distance d 5 entre l'échantillon 10 et la matrice de pixels du capteur d'image 16 est préférentiellement comprise entre 50 (lm et 2 cm, de préférence comprise entre 100 (Lm et 2 mm. On remarque, dans ce mode de réalisation, l'absence d'optique de grossissement ou de formation d'image entre le capteur d'image 16 et l'échantillon 10. Cela n'empêche pas la présence éventuelle de microlentilles de focalisation au niveau de chaque pixel du capteur 10 d'image 16, ces dernières n'ayant pas de fonction de grandissement de l'image acquise par le capteur d'image, leur fonction étant d'optimiser l'efficacité de détection. La source de lumière 11 peut comporter des sources de lumières élémentaires, émettant selon différentes bandes spectrales. Le capteur d'image est alors configuré pour acquérir, simultanément ou successivement, une image /0 dans chaque bande spectrale. Ainsi, le terme 15 image /0 acquise par le capteur d'image peut désigner un ensemble d'images acquises dans différentes bandes spectrales, suite à l'illumination de l'échantillon dans les différentes bandes spectrales. Comme évoqué dans les demandes citées dans l'art antérieur, sous l'effet de l'onde lumineuse incidente 12, les particules 10i présentes dans l'échantillon peuvent engendrer une onde 20 diffractée 13, susceptible de produire, au niveau du plan de détection P0, des interférences, en particulier avec une partie 12' 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 16, désignée par le terme "onde lumineuse d'exposition", peut comprendre : 25 une composante 13 résultant de la diffraction de l'onde lumineuse incidente 12 par chaque particule de l'échantillon ; une composante 12' résultant de la transmission de l'onde lumineuse incidente 12 par l'échantillon, une partie de cette dernière pouvant être absorbée dans l'échantillon. Ces composantes forment des interférences dans le plan de détection. Aussi, chaque image 30 acquise par le capteur d'image comporte des figures d'interférences (ou figures de diffraction), chaque figure d'interférence pouvant être associée à une particule 10i de l'échantillon. 9 Un processeur 20, par exemple un microprocesseur, est apte à traiter chaque image Io acquise par le capteur d'image 16. En particulier, le processeur est un microprocesseur relié à une mémoire programmable 22 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 5 processeur peut être couplé à un écran 24 permettant l'affichage d'images acquises par le capteur d'image 16 ou calculées par le processeur 20. Une image /0 acquise par le capteur d'image 16, également appelée hologramme, peut faire l'objet d'une reconstruction, dite reconstruction holographique. Comme décrit en lien avec l'art antérieur, on peut appliquer, à l'image acquise par le capteur d'image, un opérateur de 10 propagation holographique h, de façon à calculer une amplitude complexe A (x, y, z) représentative de l'onde lumineuse d'exposition 14, et cela en tout point de coordonnées (x, y, z) de l'espace, et plus particulièrement entre le capteur d'image et l'échantillon. Les coordonnées (x, y) désignent des coordonnées, dites coordonnées radiales, parallèlement au plan de détection Po. La coordonnée z est une coordonnée selon l'axe de propagation Z, 15 traduisant une distance entre l'échantillon 10 et le capteur d'image 16. L'amplitude complexe peut être obtenue par une des expressions suivantes : A(x, y, z) = lo(x, y, z) * h , * désignant l'opérateur produit de convolution, ou, et de préférence, A (x, y, z) = 1-10(x, y, z) * h, ou encore : 20 A (x, y, z) =10(x.Y.z) r n, étant une moyenne de l'image acquise. /7) L'opérateur de propagation h a pour fonction de décrire la propagation de la lumière entre le capteur d'image 16 et un point de coordonnées (x, y, z), situé à une distance I z I du capteur d'image. Il est alors possible de déterminer une propriété de l'onde lumineuse d'exposition 14, par 25 exemple le module M(x, y, z) et/ou la phase rp (x, y, z), à la distance I z I avec : - M (x, y, z) = abs [A (x, y, z)] ; - y(x, y, z) = ar g [A (x, y, z)] ; Les opérateurs abs et ar g désignent respectivement le module et l'argument. La distance I z I est une distance de reconstruction. 30 L'opérateur de propagation est par exemple la fonction de Fresnel-Helmholtz, telle que : 7 h(x, y, z) = - e 127 exp(jm). /az 10 L'expression complexe A(x, y, z) de l'onde lumineuse 14, en tout point de coordonnées (x, y, z) de l'espace, est telle que : A (x, y, z) = M(x, y, z)ei9Kx'Y'z). L'expression complexe A est une grandeur complexe dont l'argument et le module sont respectivement représentatifs de la phase et de l'intensité de l'onde lumineuse d'exposition 14 5 détectée par le capteur d'image 16 pour former l'image lo. En déterminant l'amplitude complexe, pour une même position radiale (x, y), le long de l'axe Z, selon une pluralité de coordonnées z, il est possible de former un profil représentatif de l'onde lumineuse d'exposition. Il peut s'agir d'un profil de la phase ou du module de l'onde lumineuse d'exposition. De façon générale, il s'agit d'un profil d'une propriété optique de l'onde lumineuse 10 d'exposition, le terme propriété optique désignant une propriété obtenue à l'aide de l'amplitude complexe A(x, y, z), et représentative de cette dernière. Il peut s'agir du module, de la phase, de la partie réelle, de la partie imaginaire, ou de leurs combinaisons. Selon un mode de réalisation, on applique un produit de convolution de l'image lo par l'opérateur de propagation h. Cela permet d'obtenir une image complexe A, représentant une 15 distribution spatiale de l'expression complexe A dans un plan de reconstruction Pz, s'étendant à une distance IzI du plan de détection Po. Dans cet exemple, le plan de détection P0 a pour équation z = O. L'image complexe A, correspond à une image complexe de l'échantillon dans le plan de reconstruction Pz. Elle représente également une distribution spatiale bidimensionnelle des propriétés optiques de l'onde lumineuse d'exposition 14. Un tel procédé, 20 désigné par le terme reconstruction holographique, permet notamment de reconstruire une image du module ou de la phase de l'onde lumineuse d'exposition 14 dans le plan de reconstruction Il est possible de former des images M, et cpz représentant respectivement le module ou la phase d'une image complexe complexe Az dans un plan Pz situé à une distance I z I du plan de 15 détection P0, avec M, = mod (Ai) et cpz = arg(Az). Lorsque le plan de reconstruction correspond à un plan selon lequel s'étend l'échantillon, les images M, et (pz permettent une observation de l'échantillon avec une résolution spatiale correcte. Lorsque des images complexes Az sont formées selon différentes distances de reconstruction, on obtient une pile d'images complexes Az, ... Azn, chaque image complexe étant définie selon 30 des coordonnées radiales (x, y), respectivement à des distances de reconstruction z' par rapport au plan de détection. Chaque image complexe A, de la pile d'images correspond à une distance de reconstruction z. A partir de la valeur de chaque image complexe L1 en une coordonnée radiale (x1, y3, correspondant à une particule 10i, on peut obtenir la valeur d'une propriété optique selon les différentes distances de reconstruction, dont l'interpolation permet de former un profil Fxbyi. Lorsque la propriété optique considérée est le module, le profil Fxiyi est noté M. Lorsque la propriété optique considérée est la phase, le profil Frnyi est noté 5 qi. Le profil correspond à une évolution de la propriété optique selon un axe parallèle à l'axe de propagation Z, et passant par la coordonnée radiale (xi, yi) de la particule 10i. Les images complexes A x Az peuvent être formées en utilisant un algorithme de reconstruction holographique itératif, tel que décrit dans W02016189257 ou dans W02017162985. Selon de tels algorithmes, on ajuste progressivement la phase de l'onde 10 lumineuse d'exposition) dans le plan de détection. Dans W02016189257, l'ajustement de la phase est effectué par itérations en moyennant, au cours de chaque itération, la phase d'ondes lumineuses reconstruites dans le plan de l'échantillon, selon différentes bandes spectrales. Dans W02017162985, l'ajustement de la phase est effectué par itération de façon à minimiser, lors de chaque itération, le bruit de reconstruction d'une image complexe reconstruite dans le plan 15 de l'échantillon. Selon une première possibilité, décrite dans W02017050672, une pile d'images complexes Azn est obtenue en mettant en oeuvre un algorithme itératif de reconstruction à partir d'une image Io acquise par le capteur d'image. L'algorithme itératif est appliqué de façon à reconstruire une image complexe, dite image complexe de référence A10, dans le plan de 20 l'échantillon P10. On suppose que l'image complexe de référence A10 forme un bon descripteur de l'onde lumineuse d'exposition 14. Elle comporte en particulier des estimations fiables du module et à la phase de l'onde lumineuse d'exposition 14 dans le plan de l'échantillon P10. Les autres images complexes formant la pile d'images Az, Azn sont obtenues par une simple convolution de l'image complexe de référence A10 par un opérateur de propagation h tel que 25 précédemment décrit. Selon une autre possibilité, chaque image complexe Az de la pile d'images Az, ... Azn est obtenue par application successive d'un algorithme itératif de reconstruction à partir de l'image /0 acquise par le capteur d'image, en considérant les différentes distances de reconstruction 30 Le recours à une pile d'images complexes n'est pas nécessaire pour établir un profil Fr Un profil F peut être obtenu en estimant l'amplitude complexe A(x1,y1,z) de l'onde lumineuse d'exposition en différentes coordonnées z, pour des mêmes coordonnées radiales (xi,yi), et 12 cela à partir d'une image acquise /0 ou de plusieurs images acquises dans des bandes spectrales différentes. Un aspect important de l'invention est l'utilisation de profils modélisés établis à partir de particules modélisées 10(par). Pour cela, on prend en compte un ensemble de paramètres par 5 d'une particule. Puis on modélise, par une modélisation numérique, l'amplitude complexe d'une onde lumineuse d'exposition 14',,,,d, se propageant entre la particule modélisée et le capteur d'image 16, et résultant d'une illumination de la particule modélisée 10(par) par la source de lumière 11. La modélisation peut notamment se basée sur la théorie de Mie. La théorie de Mie est un 10 modèle de diffusion élastique permettant d'obtenir une solution aux équations de Maxwell décrivant une onde lumineuse diffractée par une particule sphérique illuminée par une onde lumineuse incidente monochromatique de longueur d'onde À. Outre sa forme sphérique, la particule est caractérisée par un indice de réfraction n, ce dernier pouvant notamment être un indice de réfraction complexe n = Re(n) + jIm(n), avec /2 = -1. Re et lm sont 15 respectivement les opérateurs partie réelle et partie imaginaire. Une particule modélisée peut également être caractérisée par sa dimension, par exemple son diamètre ou son rayon. L'application du modèle de diffusion de Mie permet de simuler une figure de diffraction I mod(par) à différentes distances de l'échantillon. Cela permet notamment de former une figure de diffraction I Qmod(par) dans le plan de détection Po, comme 20 représenté sur les figures 2A à 2C et 3A à 3C. La notation I OEntod(par) désigne le fait que la figure de diffraction est modélisée dans le plan de détection Pu, en fonction d'un ensemble de paramètres par. L'ensemble de paramètres par peut former un vecteur de paramètres, au sens ou l'ensemble de paramètres peut être représenté sous une forme vectorielle. Les figures 2A, 2B et 2C sont des simulations d'une figure de diffraction omtod (par), dans le 25 plan de détection, résultant de la diffraction d'une onde lumineuse de longueur d'onde 450 nm, par une particule modélisée 10(par) dont le diamètre est respectivement égal à 20 pm, 15 i_xm et 10 p.rn se propageant dans un milieu d'indice égal à 1.33. L'indice de réfraction n de la particule est de 1.37, soit un saut d'indice an égal à 0.04 par rapport à l'indice du milieu. Sur ces figures, la partie imaginaire Im(n) de l'indice de réfraction est considérée comme nulle. 30 Afin de prendre en compte les imperfections du capteur d'image 16, les figures de diffraction modélisées ont été échantillonnées spatialement en considérant un pas de pixels, ce dernier étant dans cet exemple égal à 1.67 (Lm. Les figures ont ensuite été floutées par application d'un 13 filtrage gaussien pour tenir compte du bruit du capteur d'image. La distance entre le capteur d'image 16 et la particule modélisée 10(par) a été considérée égale à 1000 p.m. Chacune de ces figures a été obtenue en modélisant, dans le plan de détection Po, l'amplitude complexe A(x,y,z) de l'onde lumineuse d'exposition 14rod se propageant vers le plan de 5 détection, l'onde lumineuse d'exposition modélisée résultant de l'illumination de la particule modélisée. On a ensuite extrait le module de l'amplitude complexe A(x, y, z), dans le plan de détection Po, de façon à pouvoir simuler la figure de diffraction formée sur l'image acquise par le capteur d'image. A partir de la figure de diffraction modélisée 107'd(par), on peut former un profil, dit profil 10 modélisé F (par), représentant une évolution, parallèlement à l'axe de propagation Z, de l'amplitude complexe A(x, y, z) de l'onde lumineuse d'exposition modélisée 14rnad. Le profil modélisé F(par) est de préférence établi, à partir de la figure de diffraction modélisée /0,',',t, de la même façon que le profil Fzy a été établi en partant de l'image acquise I. Lorsque le profil Fzy a été établi en formant une pile d'images complexes à partir de l'image acquise /0, le profil 15 modélisé F (par) est établi en formant une pile d'images complexes Azi,mod ...Azwinm à partir de la figure de diffraction modélisée Le profil modélisé F (par) représente une modélisation de l'évolution de l'amplitude complexe Amoa(x, y, z) de l'onde lumineuse d'exposition modélisée 14,'od, et passant par une position radiale déterminée. La position radiale peut notamment correspondre au centre de la figure de diffraction, ce qui correspond 20 également aux coordonnées du centre de la particule modélisée dans un plan parallèle au plan de détection. Les figures 2D, 2E et 2F sont des profils F (par) du module de l'expression complexe modélisée Amod(x, y, z), passant par une coordonnée radiale (x, y) située au centre d'une particule 10(par) modélisée, dont le diamètre est respectivement égal à 20 [lm, 15 prn et 10 pm. Sur 25 chaque figure, on a considéré trois distances entre le capteur d'image et la particule modélisée, les distances s'élevant respectivement à 1000 pm, 1800 pm et 2600 p.m. Sur les graphes, l'abscisse 0 correspond à la position de la particule. L'unité de l'axe des abscisses est le pm. Les figures 3A, 3B et 3C sont des simulations d'une figure de diffraction résultant de la diffraction d'une onde lumineuse de longueur d'onde 450 nm, par une particule dont le diamètre est 30 respectivement égal à 20 pm, 15 pm et 10 Inn se propageant dans un milieu d'indice égal à 1.33. L'indice de réfraction n de chaque particule est de 1.37, soit un saut d'indice An égal à 0.01 par rapport à l'indice du milieu. Les figures 3D, 3E et 3F sont des profils F (par) du module de 14 l'expression complexe modélisée, passant par une coordonnée radiale (x, y) située au centre d'une particule modélisée 10(par), dont le diamètre est respectivement égal à 20 pm, 15 prn et 10 pm. Sur chaque figure, on a considéré trois distances entre le capteur d'image et la particule, les distances s'élevant respectivement à 1000 gm, 1800 pm et 2600 pm. 5 La comparaison des figures 2D/3D, 2E/3E et 2F/3F montre que pour une même taille de particule, une faible variation d'indice de réfraction a une influence importante sur le profil. Les figures 2A à 2F et 3A à 3F montrent qu'on peut modéliser l'amplitude complexe de l'onde lumineuse d'exposition 14 à laquelle est exposé le capteur d'image 16. Cela permet en particulier de modéliser numériquement des profils F (par), par exemple des profils du module ou de la 10 phase, correspondant à des particules dont les paramètres (taille, indice de réfraction n) sont connus. Plus généralement, la modélisation permet d'établir un profil F (par), représentatif de l'amplitude complexe de l'onde lumineuse d'exposition, en considérant une particule dont le paramètre, ou les paramètres par sont connus. Par profil représentatif de l'amplitude complexe de l'onde lumineuse d'exposition, on entend un profil d'une grandeur établie à partir de 15 l'amplitude complexe : la grandeur peut être établie à partir du module, ou de la phase, ou de la partie réelle, ou de la partie imaginaire. Les paramètres de la particule forment un ensemble par pouvant comprendre : - un indice de réfraction n de la particule, cet indice pouvant être une grandeur complexe, ou une variation de l'indice de réfraction n de la particule par rapport au milieu ambiant 20 10m; - et/ou une dimension de la particule, de type rayon r ou diamètre d; - et/ou une distance z de la particule par rapport au plan de détection. Chaque profil modélisé dépend également de la longueur d'onde X. Il s'agit d'une différence notable par rapport à l'art antérieur, dans lequel les profils étalons sont 25 obtenus expérimentalement, en utilisant des échantillons connus. L'invention permet d'éviter le recours à des échantillons connus pour former des profils étalons. De plus, le fait de disposer de profils modélisés permet d'obtenir des valeurs quantitatives précises des paramètres d'une particule. L'invention permet notamment d'obtenir une valeur quantitative de l'indice de réfraction d'une particule. 30 Dans les paragraphes qui précèdent, en lien avec les figures 2A à 2F et 3A à 3F, on a décrit l'obtention de profils modélisés F(par) selon deux étapes : - la modélisation d'une figure de diffraction /0,mod dans le plan de détection ; 15 à partir de la figure de diffraction modélisée, l'application d'un algorithme de reconstruction numérique pour établir un profil F(par), parallèlement à l'axe de propagation Z, d'une amplitude complexe A,'d(x, y, z) d'une onde lumineuse d'exposition modélisée 14',,,d. 5 De façon alternative, les profils modélisés peuvent être obtenus sans nécessairement modéliser une figure de diffraction /0,mod dans le plan de détection. On modélise alors l'amplitude complexe de l'onde lumineuse d'exposition 14,,,,d le long de l'axe de propagation Z. Cependant, les inventeurs estiment qu'il est préférable de procéder à une modélisation de la figure de diffraction 10,,,,,d(par), puis à la formation du profil F(par) correspondant à la particule 10 modélisée de la même façon qu'on forme un profil F , passant par une coordonnée (xi, yi), à partir de l'image /0 acquise par le capteur d'image. Quelle que soit la façon de les obtenir, les profils F(par) ainsi modélisés peuvent être mémorisés dans une base de données, de façon à être mis à profit dans un procédé dont les principales étapes sont décrites ci-après, en lien avec les figures 4A à 4D. A chaque profil est 15 associé un ensemble de paramètres (par). Etape 100 : Illumination de l'échantillon 10 à l'aide de la source lumineuse 11. Etape 110: Acquisition d'une image /0 de l'échantillon 10 par le capteur d'image 16, cette image formant un hologramme. Un des intérêts de la configuration sans lentille, représentée sur la figure 1, est le large champ observé, permettant d'adresser simultanément un volume 20 d'échantillon élevé. Cela permet d'observer simultanément plusieurs particules, et ainsi d'obtenir une caractérisation rapide de l'échantillon. Le champ observé dépend de la taille du capteur d'image, en étant légèrement inférieur à la surface de détection de ce dernier, du fait de l'espacement entre les pixels du capteur et l'échantillon. Le champ observé est généralement supérieur à 10 mm2, et est typiquement compris entre 10 mm2 et 50 mm2, ce qui est 25 significativement plus élevé qu'avec un microscope. Etape 120 : Détection de particules dans l'échantillon. L'image acquise /0 comporte généralement un grand nombre de figures d'interférence. Du fait du recouvrement entre les différentes figures d'interférence, l'image acquise n'est généralement pas aisément exploitable pour localiser les particules présentes dans le champ observé. Ces dernières sont plus facilement 30 identifiables à partir d'une image complexe reconstruite par application d'un opérateur de propagation holographique h à l'image acquise J. 16 Aussi, l'étape 120 comporte une reconstruction d'au moins une image, dite image d'observation I', de l'échantillon. On applique un opérateur de reconstruction holographique à l'image acquise, selon une distance de reconstruction, de façon à obtenir une image complexe représentant l'amplitude complexe de l'onde lumineuse d'exposition 14 dans un plan de 5 reconstruction parallèle au plan de détection et situé à la distance de reconstruction de ce dernier. L'image d'observation I' peut être l'image du module ou de la phase de l'image complexe ainsi reconstruite. Le plan de reconstruction dans lequel l'image d'observation est définie est de préférence un plan P10 selon lequel s'étend l'échantillon 10. Sa position peut être déterminée a priori, ou par un algorithme de focalisation numérique, ce type d'algorithme étant 10 connu de l'homme du métier. Sur l'image d'observation I', les particules 10i apparaissent suffisamment contrastées pour être aisément discernables du milieu ambiant 10m. Etape 130 : détermination d'une position radiale (xi, y() de chaque particule 10i. Il s'agit d'obtenir une position radiale représentative de chaque particule discernable sur l'image 15 d'observation ou sur l'image acquise /0, lorsque cette dernière est exploitable. Un algorithme de segmentation peut être appliqué à l'image d'observation, de façon à extraire des régions d'intérêt RO( i correspondant respectivement à chaque particule 10i. Etape 140 : formation d'un profil Frisi pour chaque particule 10i à laquelle correspond une position radiale (xi,yi). Il s'agit d'établir un profil Fxby i représentatif de l'évolution de 20 l'amplitude complexe de l'onde lumineuse d'exposition parallèlement à l'axe de propagation Z, le profil Fx i passant par la position radiale (xi, yi) associée à chaque particule 10. Le profil Fxisi peut être un profil du module de l'onde lumineuse d'exposition : Exisi = Mxiyi. Le profil F.whyi peut être un profil du module de l'onde lumineuse d'exposition :F,. = Il peut également s'agir d'un profil représentant une évolution de la partie réelle ou de la partie 25 imaginaire de la particule. Etape 150 : confrontation du profil Fxii,i formé pour chaque particule avec des profils modélisés F (par), chaque profil modélisé étant paramétré par un ensemble de paramètres par comme précédemment évoqué. L'ensemble de paramètres par comporte par exemple un indice de réfraction n, ou un saut d'indice de réfraction an du à la particule. L'étape 150 30 comporte une détermination du profil modélisé F (par) le plus proche du profil F formé pour la particule. Les paramètres part de la particule 10i sont les paramètres associés au profil modélisé le plus proche Ft (par) 17 Sur le procédé représenté sur la figure 4A, on compare le profil à chaque profil modélisé F (par), de façon à déterminer les paramètres du profil minimisant la comparaison. La comparaison peut par exemple comporter un écart quadratique, auquel cas : par, = argmin[F, y F (par)] ( 1) par L' 5 Lorsque le profil est un profil représentant une évolution du module de l'amplitude complexe de l'onde lumineuse d'exposition, pari = argmin[Mx y, - M(par)]2 (2), par OÙ M (par) correspond à des profils modélisés du module d'une amplitude complexe d'une onde lumineuse d'exposition modélisée 14mod en prenant en compte différentes valeurs du vecteur de paramètres par. 10 Lorsque le profil Fi est un profil représentant une évolution de la phase de l'amplitude complexe 2 de l'onde lumineuse d'exposition, pari = argmin[cpx y - cp(par).1 (3) par où co(par) correspond à des profils modélisés de la phase d'une amplitude complexe d'une onde lumineuse d'exposition modélisée l4mod en prenant en compte différentes valeurs du vecteur de paramètres par. 15 Selon un mode de réalisation, la détermination des paramètres d'une particule peut combiner différents profils, par exemple sous la forme d'une minimisation d'une somme pondérée, de type : pari = argmin (ki[Mxi-Yi- M(par)]2 ± k2Pz1svi cp - (par)2] ) (4) où k1 et k2 sont des par scalaires, formant des termes de pondération. 20 De préférence, le vecteur de paramètres par d'une particule comprend au moins l'indice de réfraction. Lorsque l'indice de réfraction est exprimé sous la forme d'une grandeur complexe, les paramètres comprennent la partie réelle de l'indice de réfraction Re(n) et sa partie imaginaire Im(n). Comme précédemment indiqué, les paramètres peuvent comprendre une dimension (diamètre ou rayon) ou une distance z de la particule par rapport au plan de 25 détection, selon l'axe de propagation Z. L'étape 150 requiert le recours à des profils modélisés F (par), pour différentes valeurs du vecteur de paramètres par. Comme décrit en lien avec les figures 2D à 2F, ou 3D à 3F, de tels profils peuvent être obtenus par modélisation numérique, en particulier selon le modèle de diffusion de Mie. 18 Lorsque l'algorithme se base sur une base de données de profils, la formation de cette dernière. fait l'objet des étapes 90 et 95. Au cours d'une étape 90, on prend en compte les paramètres que l'on souhaite déterminer, ces paramètres formant un ensemble de paramètres, pouvant prendre la forme d'un vecteur de paramètres. Au cours d'une étape 95, on effectue différentes 5 modélisations pour différentes valeurs du vecteur de paramètres par de façon à obtenir, pour chaque valeur du vecteur de paramètres par, un profil modélisé F(par). L'étape 95 peut également comporter une interpolation entre des profils modélisés F(par), F(par'), où par' est un vecteur de paramètres dont les valeurs sont proches du vecteur par. L'interpolation permet d'obtenir des profils correspondant à des paramètres compris entre par et par'. 10 Selon un autre mode de réalisation, illustré sur la figure 4B, on peut appliquer, lors de l'étape 150, un algorithme itératif permettant d'estimer un profil F(par) le plus proche du profil mesuré F. Selon ce mode de réalisation, on ajuste progressivement le vecteur de paramètres par de la particule modélisée de telle sorte que le profil F(par) correspondant au vecteur de paramètres par tende vers le profil mesuré Fxbyi. Au cours de chaque itération q, on détermine 15 un écart Eq entre le profil mesuré F et un profil F(par). La notation para correspond au vecteur de paramètres pris en compte lors de l'itération q. L'écart Eq peut être un écart quadratique moyen ou un écart moyen entre le profil mesuré E et le profil F(par). Il s'agit de préférence d'un scalaire. Le vecteur de paramètres parq+1 de l'itération suivante q + 1 est déterminé de façon à minimiser l'écart 8'1+1 entre le profil mesuré F et le profil F(par") 20 établi à l'aide du vecteur de paramètres para+1. Ainsi, chaque itération q vise à estimer un vecteur de paramètres q+1 permettant de minimiser l'écart. Le procédé se réitère jusqu'à un arrêt des itérations, cet arrêt correspondant à l'atteinte d'un critère de convergence ou d'un nombre d'itérations prédéterminé. Le vecteur de paramètres para" considéré lors de l'itération suivante peut être estimé par un 25 algorithme de type descente de gradient, au cours duquel à chaque itération q, on détermine un gradient de l'écart Veq, ce dernier correspondant à une variation de l'écart eq en fonction d'un ou plusieurs paramètres, et de préférence chaque paramètre, du vecteur de paramètres para. Le vecteur de paramètres para+1 pris en compte lors de l'itération suivante est déterminé en fonction de Vs", de façon à minimiser l'écart eq÷1. 30 Sur la figure 4B, on a représenté les sous-étapes 151, 152, et 153, correspondant respectivement à la formation du profil modélisé F(par), le calcul de l'écart Eq et de son gradient V e , de façon à définir les paramètres para +1 à prendre en compte lors de l'itération suivante. 19 Lors de la première itération (q=1), l'algorithme itératif est initialisé par un vecteur de paramètres initial par'. Le vecteur de paramètres initial peut être déterminé a priori. Selon un mode de réalisation, les deux modes de réalisation précédemment décrits sont combinés : on utilise une base de données de profils et on détermine le vecteur de paramètres 5 pari minimisant la comparaison entre les profils de la base de données F(par) et le profil mesuré Fxtyi. Puis on ajuste progressivement le vecteur de paramètres correspondant à la particule observée en mettant en oeuvre un algorithme d'ajustement itératif, par exemple de type descente de gradient. L'algorithme d'ajustement itératif est alors initialisé par le vecteur de paramètres par,. Autrement dit, selon ce mode de réalisation, on met en oeuvre les étapes 10 100 à 150 représentées sur la figure 4A. Le vecteur de paramètres pari, résultant de l'étape 150 de la figure 4A, forme le vecteur initial para-1- permettant une mise en oeuvre de l'étape 150 de la figure 4B. Selon une variante, dont les étapes sont représentées sur la figure 4C, l'étape 150 met en oeuvre un algorithme de minimisation, tel que décrit en lien avec la figure 4A. L'étape 150 est initialisée 15 sur la base d'un a priori basé sur une estimation d'un des paramètres recherchés. Dans l'exemple représenté sur la figure 4C, on prend en compte une dimension ri d'une particule 10i, la dimension ayant été estimée par l'image d'observation /' résultant de l'étape 120. Le fait d'initialiser l'algorithme de minimisation en utilisant un a priori, basé sur une estimation d'un paramètre, permet de réduire le temps de calcul et d'améliorer la précision de la détermination 20 du vecteur des paramètres. Les figures SA à 5C illustrent cet aspect. La figure .5A montre une application de l'invention sur une cellule flottante 101 de type CHO, baignant dans un milieu de culture. La cellule 10i est illuminée par une source de lumière 11 de longueur d'onde 450 nm. La source de lumière est une diode électroluminescente émettant dans la bande spectrale 450nm - 465 nm. Elle est disposée à une distance de 1400 ttm d'un 25 capteur d'image monochrome (référence IDS 3840x2748 pixels). Les paramètres recherchés sont par = [ri, Ani]. On a acquis une image 1, de la cellule, à partir de laquelle une image d'observation l' a été obtenue, par reconstruction holographique. Le rayon de la cellule ri a été estimé par l'image d'observation, représentant le module de l'amplitude complexe dans le plan de reconstruction. L'image d'observation /' a également permis de définir les coordonnées 30 y1) du centre de la cellule 10i dans le plan de détection. A partir de l'image acquise, on a déterminé un profil , du module de l'onde lumineuse d'exposition 14, le long d'un axe 20 parallèle à l'axe de propagation Z et passant par les coordonnées radiales de la cellule. Le profil F est appelé profil mesuré. La courbe a montre un profil F (init) utilisé pour l'initialisation de l'algorithme, le profil utilisant le rayon estimé sur une image d'observation de l'échantillon I'. Le rayon est estimé égal à 8.83 5 Ftm. Ce profil correspond aux paramètres [ri = 8.15 Fun; Re(Ani) = 0.025]. La courbe b correspond au profil mesuré F. yi . Les courbes c et ci correspondent au profil modélisé le plus proche Ft (par), respectivement sans considérer et en considérant la partie imaginaire Im(n) de l'indice de réfraction n. Les paramètres pari associés aux courbes c et cl sont respectivement : [ri = 8.15 ptm,Re(ani) = 0.027] et [ri = 8.15 11m, Re(Ani) = 0.025, /m(Ani) = 0.01]. Ces 10 paramètres sont considérés comme correspondant à ceux de la cellule examinée. La figure 5B montre un autre exemple d'application sur une autre cellule du même type. Dans cet exemple, on a estimé, à partir d'une image d'observation l', le rayon ri de la cellule à 6.74 et L'algorithme de minimisation a été initialisé en considérant le profil modélisé représenté sur la courbe a, les paramètres associés à ce profil étant [n = 6.74 pn, Re(Ani) = 0.01, 15 /m(Ani) = 0.00]. La courbe b correspond au profil mesuré La courbe c correspond au profil modélisé le plus proche r(par), les paramètres associés à ce profil étant [ri = 6.09 Ryn, Re(Ani) = 0.014, im(Ant) = 0.007]. Ces paramètres sont considérés comme correspondant à ceux de la cellule examinée. La figure 5C complète la figure 5B. Sur la figure 5C, on a représenté le profil mesuréF (courbe 20 b), ce dernier étant également représenté sur la figure 5B. Sur la figure 5C, le profil a correspond au profil modélisé avec lequel l'algorithme de minimisation est initialisé, ce dernier correspondant aux paramètres [ri = 6.68 pm, Re (Ani) = 0.025, tin(ûni) = 0.0]. La courbe c correspond au profil modélisé le plus proche, les paramètres associés à ce profil étant [ri 6.25 pm, Re(An() = 0.012, im(ani) = 0.008]. Ces paramètres sont considérés comme 25 correspondant à ceux de la cellule examinée. Bien que l'algorithme de minimisation soit initialisé en considérant des valeurs d'indice différentes, il permet d'obtenir des profils les plus proches F* (par) (profils c des figures 5B et 5C) dont les paramètres respectifs sont cohérents : [ri = 6.09 pm, Re(Ani) = 0.014, im(Ani) = 0.007] pour la figure 5B et [ri = 6.25 p,m, Re(Ani) = 0.012, /m(Ani) = 0.008] sur la figure 5C. 30 Selon un mode de réalisation, dont les étapes sont illustrées sur la figure 4D, l'étape 150 de l'algorithme est réitérée. De préférence, l'étape 150 prend en compte une estimation de la taille de la particule résultant d'une image d'observation. Au cours d'une première itération, l'étape 21 150 résulte en une première estimation de l'indice de réfraction. Cette première estimation est utilisée pour définir le profil d'initialisation F (init) de la deuxième itération. Le nombre d'itérations peut être supérieur à 2, mais les inventeurs ont constaté que des résultats satisfaisants étaient obtenus au bout de deux itérations. Les figures 6A et 6B représentent 5 l'erreur de quantification affectant respectivement la partie réelle et la partie imaginaire d'un indice de réfraction sans réitération de l'étape 150 (courbes ci ) et avec réitération de l'étape 150 (courbes b). Les erreurs sont représentées en fonction du rayon des cellules. Ces courbes ont été obtenues sur la base de modélisations. La figure 7 montre un exemple permettant de classifier des cellules en fonction de la partie 10 imaginaire et de la partie réelle de l'indice de réfraction. Sur cette figure, on a représenté la partie réelle et la partie imaginaire de l'indice de réfraction de différentes cellules, les valeurs représentées étant obtenues en mettant en oeuvre l'invention. On constate que la comparaison entre la partie réelle et la partie imaginaire est un indicateur pertinent relatif à l'état de la cellule, parmi un état vivant ou un état mort. On a tracé une ligne séparant deux populations de cellules 15 en fonction de leurs états respectifs. La figure 8 schématise un dispositif permettant une mise en oeuvre de l'invention. Contrairement au dispositif représenté sur la figure 1, le dispositif de la figure 8 comporte un système optique 19 de formation d'image. Le système optique 19 définit un plan image et un plan objet. Le système optique peut être une lentille ou un objectif. Lors de l'acquisition de 20 l'image de l'échantillon, le capteur d'image est disposé selon une configuration défocalisée. 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 [Lm -500 Rrn. L'invention pourra être mise en oeuvre pour la caractérisation de particules dans le domaine de 25 la biologie ou de la santé. D'autres applications peuvent être envisagées, par exemple dans le contrôle de l'environnement, ou de procédés industriels, ou dans le domaine de l'agroalimentaire. REVENDICATIONS1. Procédé de détermination de paramètres d'une particule (10;) présente dans un échantillon (10), l'échantillon s'étendant entre un capteur d'image (16) et une source de lumière (11), le capteur d'image s'étendant selon un plan de détection (Po), le procédé comportant les étapes suivantes : a) illumination de l'échantillon par la source de lumière, la source de lumière émettant une onde lumineuse incidente (12) se propageant selon un axe de propagation (Z) ; b) acquisition d'une image de l'échantillon (/0) par le capteur d'image (16), le capteur d'image étant exposé à une onde lumineuse d'exposition (14); c) détermination d'une position (xi, yi) de la particule dans le plan de détection ou dans un plan parallèle au plan de détection ; d) à partir de l'image acquise (4), application d'un opérateur de propagation (h), selon une pluralité de distances (z) du plan de détection, de façon à estimer, à chaque distance, une amplitude complexe (A(x,, yi, z)) de l'onde lumineuse d'exposition (14); e) à partir de l'amplitude complexe estimée, en différentes distances (z), obtention d'un profil (F) représentant une évolution de l'amplitude complexe de l'onde lumineuse d'exposition (A(x,, yi, z)) selon un axe parallèle à l'axe de propagation (2) et passant par la position (xi, yi) de la particule ; le procédé étant caractérisé en ce que la particule est associée à un ensemble de paramètres (par), comportant au moins une taille de la particule et un indice de réfraction (n, an) de la particule, le procédé comportant également : f) une modélisation d'une particule (10(par)), prenant en compte une valeur de chaque paramètre de l'ensemble de paramètres (par), et une modélisation d'une onde lumineuse d'exposition (14mod), se propageant vers le capteur d'image, sous l'effet d'une illumination, par la source de lumière (11), de la particule modélisée ; g) à partir de l'onde lumineuse d'exposition modélisée (14'od), formation d'un profil (F (par)) représentant une évolution de l'amplitude complexe de l'onde lumineuse d'exposition modélisée, selon un axe parallèle à l'axe de propagation (Z) ; h) comparaison entre le profil (F) obtenu lors de l'étape e) et le profil (F (par)) modélisé résultant de l'étape g), de façon à déterminer la valeur d'au moins un paramètre (pari) de la particule. 23 Z. Procédé selon la revendication 1, dans lequel l'indice de réfraction (n) comporte une partie réelle (Re(n)) et une partie imaginaire (Im(n)). 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'ensemble de paramètres (par) comporte également une distance (z) de la particule par rapport au plan de 5 détection (P0). 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel : l'étape f) comporte des modélisations de particules présentant respectivement différentes valeurs d'au moins un paramètre, de façon à obtenir, suite à l'étape g), une base de données de profils modélisés, à chaque profil modélisé étant associé un 10 ensemble de paramètres ; l'étape h) comporte une minimisation d'un écart entre le profil (F,i) résultant de l'étape e) et les profils (F(par)) de la base de données, les valeurs respectives des paramètres de la particule étant celles minimisant l'écart. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel les étapes f) à 15 h) sont effectuée de façon itérative, la valeur d'un paramètre de la particule (n, An) résultant d'une itération étant utilisée pour initialiser une itération suivante. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel : - l'étape f) comporte des modélisations de particules de différentes tailles et/ou de différents indices de réfraction, de façon à obtenir, suite à l'étape g), une base de 20 données de profils modélisés, à chaque profil modélisé étant associé une taille et/ou un indice de réfraction; - l'étape h) comporte une détermination de la taille de la particule et/ou de l'indice de réfraction de la particule. 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel : 25 - l'étape f) des modélisations de particules situées à différentes distances du plan de détection, de façon à obtenir, suite à l'étape g) une base de données de profils modélisés, chaque profil modélisé étant associé à une distance par rapport au plan de détection ; - l'étape h) comporte une détermination de la distance entre la particule et le plan de 30 détection. 24 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel les étapes f) à h) sont mises en oeuvre de façon itérative, de telle sorte qu'au cours de chaque itération, le profil modélisé lors de l'étape g) se rapproche progressivement du profil déterminé lors de l'étape e). 9. Procédé selon la revendication 8, dans lequel l'étape h) d'une itération (q) comporte : 5 - une détermination d'un écart (eq) entre le profil modélisé F (parq) lors de l'étape g) de la même itération et le profil (Fx,),,) résultant de l'étape e) ; - une détermination d'un gradient de l'écart (eq) en fonction d'au moins un paramètre de l'ensemble de paramètres, de façon à déterminer les valeurs des paramètres de la particule modélisée lors de l'étape f) de l'itération suivante (q + 1). 10 10. Procédé selon l'une quelconque des revendications précédentes, comportant également : - une application d'un opérateur de propagation (h) à l'image acquise, en considérant plusieurs distances (z) par rapport au plan de détection, de façon à obtenir une pile d'images reconstruites, à chaque image reconstruite étant associée une distance par rapport au plan de détection ; 15 - à partir d'une image reconstruite (/'), estimation d'une taille (ri, di) de la particule ; de telle sorte que lors de l'étape h) comporte une prise en compte de la taille de la particule estimée à partir de l'image reconstruite. 11. Procédé selon l'une quelconque des revendications précédentes, dans lequel la particule est une cellule, le procédé comportant une détermination d'un état de vie ou de mort de la cellule 20 en fonction de l'indice de réfraction. 12. Procédé selon la revendication 11, dans lequel l'indice de réfraction comporte une partie réelle et une partie imaginaire, l'état de vie ou de mort de la cellule étant défini en fonction d'une comparaison entre la partie réelle (Re(n)) et la partie imaginaire (Im(n)). 13. Procédé selon l'une quelconque des revendications précédentes, dans lequel il n'y a pas 25 d'optique de formation d'image entre l'échantillon et le capteur d'image. 14. Procédé selon l'une quelconque des revendications 1 à 12, dans lequel un système optique (19) s'étend entre l'échantillon et le capteur d'image, le système optique définissant un plan image et un plan objet, et dans lequel lors de l'étape b), l'image est acquise selon une configuration défocalisée, le plan de détection étant décalé par rapport au plan image, et/ou un 30 plan de l'échantillon, selon lequel s'étend l'échantillon, étant décalé par rapport au plan objet. 25 15. Dispositif d'observation d'un échantillon (10), l'échantillon comportant des particules, le dispositif comportant : - une source de lumière (11), configurée pour illuminer l'échantillon ; - un support d'échantillon (10s), configuré pour recevoir l'échantillon ; 5 - un capteur d'image (16), disposé de telle sorte que lorsque l'échantillon est disposé sur le support, il s'étend entre le capteur d'image et la source de lumière, le capteur d'image étant configuré pour acquérir une image de l'échantillon ; - un processeur (20), configuré pour mettre en oeuvre les étapes d) à h) de l'une quelconque des revendications précédentes, à partir d'une image acquise par le capteur 10 d'image. 16. Dispositif selon la revendication 15, dans lequel aucune optique de formation d'image n'est disposée entre le capteur d'image et le processeur. 17. Dispositif selon la revendication 15, comportant un système optique (19) s'étendant entre l'échantillon et le capteur d'image, le système optique définissant un plan image et un plan 15 objet, le dispositif comportant un moyen d'ajustement du système optique, ou de l'échantillon, ou du capteur d'image, de telle sorte que : - le plan de détection soit décalé du plan image ; - et/ou le plan objet soit décalé d'un plan d'échantillon, selon lequel s'étend l'échantillon.