FR 3028616 A1 20160520 FR 1460947 A 20141113 1 PROCEDE D'ANALYSE COMPRENANT LA DETERMINATION D'UNE POSITION D'UNE PARTICULE BIOLOGIQUE. DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne le domaine de l'analyse d'un échantillon recevant des particules biologiques, cette analyse consistant en particulier à déterminer une position d'au moins une desdites particules biologiques, selon un axe de la profondeur de l'échantillon. ÉTAT DE LA TECHNIQUE ANTÉRIEURE On connaît dans l'art antérieur de nombreux cas dans lesquels il est nécessaire de connaître la position d'une particule biologique, selon l'axe de la profondeur de l'échantillon recevant cette particule. Une solution couramment utilisée consiste à acquérir une série d'images de l'échantillon, à différentes profondeurs dans celui-ci. On recherche ensuite l'image sur laquelle on observe une image nette de la particule biologique. Un inconvénient de cette solution est qu'elle ne permet pas de connaître la position d'une particule biologique transparente dont l'image nette se confond avec l'image du milieu environnant. Un objectif de la présente invention est de proposer un procédé et un dispositif ne présentant pas cet inconvénient. En particulier, un objectif de la présente invention est de proposer un procédé et un dispositif permettant de déterminer la position de n'importe quelle particule biologique, selon au moins un axe de l'espace à trois dimensions. 3028616 2 EXPOSÉ DE L'INVENTION Cet objectif est atteint avec un procédé d'analyse d'un échantillon recevant des particules biologiques, parmi lesquelles une particule d'intérêt, le procédé 5 étant mis en oeuvre dans un dispositif comprenant l'échantillon disposé entre une première source lumineuse et un système optique, ledit procédé comprenant les étapes suivantes : définition d'un point de repère, situé sur une première interface de l'échantillon ou à une distance connue de celle-ci selon un axe parallèle à l'axe 10 optique du système optique ; - à l'aide de la première source lumineuse, illumination d'une région recevant la particule d'intérêt, dite région illuminée ; - à l'aide d'un capteur situé dans un plan image du système optique, acquisition d'une image de la région illuminée, dite image de référence, la particule 15 d'intérêt se trouvant hors du plan de mise au point du système optique, et la distance, selon un axe parallèle à l'axe optique du système optique, entre ledit plan de mise au point et le point de repère, étant une distance connue dite distance utile ; à partir de l'image de référence, construction numérique d'une série 20 d'images reconstruites, associées chacune à un décalage prédéterminé du plan de mise au point le long de l'axe optique dudit système optique ; à partir de la série d'images reconstruites, détermination de la distance selon un axe parallèle à l'axe optique du système optique, entre la particule d'intérêt et ledit plan de mise au point. 25 Avantageusement, le procédé comprend une détermination de la distance entre la particule d'intérêt et ledit point de repère, à partir de ladite distance entre la particule d'intérêt et ledit plan de mise au point ainsi que de ladite distance utile. 3028616 3 Il peut en outre comporter une détermination de la présence de ladite particule d'intérêt dans le milieu à partir de la distance entre ladite particule d'intérêt et ledit plan de mise au point. La détermination de la distance utile, comprend les sous-étapes suivantes : 5 illumination du point de repère par un faisceau laser, ledit faisceau laser étant focalisé par ledit système optique ; acquisition d'une image sur le capteur représentant la réflexion du faisceau laser sur une interface de l'échantillon ; ajustement de la distance utile selon un axe parallèle à l'axe optique du 10 système optique, entre l'échantillon et le système optique, en fonction de l'intensité de l'image formée sur le capteur. Typiquement, les particules biologiques adhèrent sur une deuxième interface de l'échantillon, les première et deuxième interfaces étant distinctes ou 15 confondues. La distance selon un axe parallèle à l'axe optique du système optique, entre la position du plan de mise au point associée à l'image de référence, et la projection du point de repère selon cet axe et sur la deuxième interface, est comprise par exemple entre +5 um et +2000 um ou entre -5 um et -2000 um. 20 La région illuminée peut comprendre une pluralité de particules biologiques. La première source lumineuse présente typiquement une largeur spectrale inférieure à 200 nm. Avantageusement, le procédé comprend en outre une détermination, à partir de ladite série d'images reconstruites, de la position de la particule d'intérêt 25 dans un plan orthogonal à l'axe optique du système optique. Chaque image reconstruite est formée par une partie réelle et une partie imaginaire, et parmi les parties réelles et imaginaires, seules les parties imaginaires des images reconstruites sont avantageusement utilisées pour déterminer la 3028616 4 distance selon un axe parallèle à l'axe optique du système optique, entre la particule d'intérêt et le point de repère. Chaque image reconstruite peut être associée à un décalage le long de l'axe optique du système optique et à une valeur d'un paramètre utile, de manière à 5 constituer une fonction décrivant l'évolution du paramètre utile selon ledit décalage, la détermination de la distance entre la particule et le plan de mise au point mettant en oeuvre une recherche d'une valeur remarquable, en particulier un extremum, un point d'inflexion ou un passage par zéro de ladite fonction. La détermination de la distance entre la particule d'intérêt et le plan de mise 10 au point comprend avantageusement les sous-étapes suivantes : à partir d'une première série d'images reconstruites associée à un premier pas des décalages du plan de mise au point, détermination d'une distance approchée entre la particule d'intérêt et ledit plan de mise au point ; - à partir d'une deuxième série d'images reconstruites associée à un 15 deuxième pas des décalages du plan de mise au point, le deuxième pas étant plus fin que le premier pas, détermination d'une distance précise entre la particule d'intérêt et ledit plan de mise au point. A partir de la distance entre la particule d'intérêt et ledit plan de mise au 20 point, le système optique est avantageusement déplacé par rapport à ladite particule d'intérêt de façon à focaliser le faisceau laser d'analyse (363) sur la particule d'intérêt. A partir de ladite distance entre la particule d'intérêt et ledit plan de mise au point, le système optique est avantageusement déplacé par rapport à ladite 25 particule d'intérêt de façon à ajuster la focalisation d'un photodétecteur situé dans le plan image. Le procédé peut comprendre en outre une étape de comptage du nombre de particules biologiques présentes dans l'échantillon. Ces particules biologiques 3028616 5 peuvent comprendre des bactéries, des spores, des cellules, des levures ou des micro-organismes. L'invention concerne également un dispositif d'analyse d'un échantillon recevant des particules biologiques, parmi lesquelles une particule d'intérêt, le 5 dispositif comprenant : une première source lumineuse ; un ensemble d'imagerie comprenant un système optique et un capteur, tels que le capteur se trouve dans un plan image du système optique, ledit plan image étant le conjugué, par le système optique, d'un plan de mise au point ; 10 un support adapté à recevoir l'échantillon, disposé entre la première source lumineuse et l'ensemble d'imagerie ; des moyens de translation, adaptés à déplacer le support relativement aux moyens d'imagerie, selon un axe parallèle à l'axe optique du système optique ; et 15 des moyens de calculs : reliés à une mémoire stockant une information relative à un point de repère situé sur une première interface de l'échantillon ou à une distance connue de celle-ci selon l'axe optique du système optique, 20 recevant en entrée une image acquise par le capteur, dite image de référence lorsque ledit plan de mise au point est décalé par rapport audit point de repère adaptés à réaliser, à partir de l'image de référence, une construction numérique d'une série d'images reconstruites, 25 associées chacune à un décalage prédéterminé du plan de mise au point le long de l'axe optique du système optique ; et fournissant en sortie la distance selon un axe parallèle à l'axe optique du système optique, entre la particule d'intérêt et le plan de mise au point. 3028616 6 BRÈVE DESCRIPTION DES DESSINS La présente invention sera mieux comprise à la lecture de la description d'exemples de réalisation donnés à titre purement indicatif et nullement limitatif, 5 en faisant référence aux dessins annexés sur lesquels : - la figure 1 illustre un premier mode de réalisation d'un dispositif d'analyse selon l'invention ; - la figure 2 illustre de manière schématique un premier mode de réalisation d'un procédé d'analyse selon l'invention ; 10 - la figure 3 illustre un deuxième mode de réalisation d'un dispositif d'analyse selon l'invention ; la figure 4 illustre de manière schématique un détail d'un deuxième mode de réalisation d'un procédé d'analyse selon l'invention ; les figures 5A à 5C illustrent un détail du procédé illustré en figure 4 ; 15 la figure 6 illustre une image de référence selon l'invention ; la figure 7 illustre de manière schématique un détail d'un troisième mode de réalisation d'un procédé d'analyse selon l'invention ; et - la figure 8 illustre un profil utilisé dans un quatrième mode de réalisation de procédé selon l'invention. 20 EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Le procédé d'analyse et le dispositif d'analyse selon l'invention seront décrits conjointement, le dispositif d'analyse selon l'invention étant adapté à la mise en oeuvre du procédé d'analyse selon l'invention. 25 La figure 1 illustre un premier mode de réalisation d'un dispositif d'analyse 100 selon l'invention. La figure 2 illustre un premier mode de réalisation d'un procédé d'analyse selon l'invention, mis en oeuvre dans le dispositif d'analyse 100. 3028616 7 Le dispositif d'analyse 100 comprend un support 110, adapté à recevoir un échantillon 111. Le support 110 est par exemple une pince, ou un plateau transparent, ou un plateau percé d'une ouverture pour laisser passer des rayons lumineux. 5 L'échantillon 111 est un milieu transparent ou translucide, c'est-à-dire présentant un coefficient de transmission supérieur ou égal à 70% dans le spectre visible ou plus généralement dans un spectre compris entre 300 nm et 1000 nm. L'échantillon 111 est délimité par des interfaces, ou frontières, ou limites entre l'échantillon et un milieu situé à son contact direct. 10 L'échantillon consiste en un liquide tel que de l'eau, une solution tampon, un liquide comprenant un réactif, un milieu de culture, et des particules biologiques 112 situées dans ce liquide. En variante, l'échantillon consiste en un milieu solide tel qu'une gélose, et des particules biologiques situées sur ce milieu solide. Selon une autre variante, l'échantillon consiste en un gaz dans lequel se trouvent des 15 particules biologiques. Ainsi, les particules biologiques peuvent être situées à l'intérieur de l'échantillon (par exemple baigner dans un liquide), ou affleurer en surface de l'échantillon (par exemple être situées sur une gélose). Les particules biologiques 112 désignent par exemple des bactéries, des spores, des cellules, des levures, ou tout type de micro-organisme. L'une de ces 20 particules biologiques est nommée particule d'intérêt 112A. Dans l'exemple représenté en figure 1, l'échantillon 111 comprend un milieu de culture enfermé dans une chambre fluidique délimitée par une lame inférieure 113, par exemple une lame de microscope standard, une lame supérieure 114, et un adhésif 115 entourant latéralement l'échantillon 111 et reliant ensemble les 25 lames supérieure et inférieure. Les lames supérieure et inférieure 113, 114 sont transparentes dans le visible et le cas échéant d'autres longueurs d'ondes utiles. La frontière entre l'échantillon et la lame inférieure 113 définit une interface inférieure 116 de l'échantillon 111. La frontière entre la lame supérieure 114 et l'échantillon définit une interface supérieure 117 de l'échantillon 111. 3028616 8 Une première source lumineuse 120 se trouve en amont de l'échantillon 111, dans le sens de propagation de la lumière, depuis la première source lumineuse vers l'échantillon 111. Dans la suite, les termes amont et aval se réfèrent au sens de propagation de la lumière, depuis la première source lumineuse 120 vers 5 l'échantillon 111. La première source lumineuse 120 est par exemple un laser, une diode électroluminescente, une lampe blanche notamment une lampe à vapeur de mercure filtrée ou non. La première source lumineuse peut comprendre une fibre optique pour amener la lumière sous la lame inférieure 113. 10 La première source lumineuse 120 est avantageusement temporellement cohérente. Elle présente une largeur spectrale avantageusement inférieure à 200 nm, voire inférieure à 100 nm ou même 25 nm. De préférence, la première source lumineuse 120 est spatialement cohérente. 15 On donnera dans la suite plus de détails sur la cohérence spatiale et temporelle de la première source lumineuse 120. La première source lumineuse illumine en transmission une région 119 de l'échantillon. La particule d'intérêt 112A se trouve dans la région 119. Un ensemble d'imagerie 130 est situé en-aval de la lame supérieure 114, 20 dans le sens de propagation de la lumière, depuis la première source lumineuse vers l'échantillon 111. Dans l'exemple décrit, l'ensemble d'imagerie 130 se trouve au-dessus de l'échantillon 111. L'ensemble d'imagerie 130 comprend un système optique 131 et un capteur 132. 25 Le système optique 131 consiste par exemple en un objectif, notamment un objectif de microscope. Il présente un plan objet 133, nommé également plan de mise au point, et un plan image 134. Le plan image 134 est le conjugué (ou image) du plan de mise au point par le système optique. En d'autres termes, un objet situé dans le plan de mise au point correspond à une image nette dans le plan image. 3028616 9 Le capteur 132 est par exemple un capteur matriciel de type CCD ou CMOS. Il est situé dans le plan image 134. Ainsi, le capteur 132 acquiert une image en transmission, d'une partie du plan de mise au point 133. Dans une configuration plan-infini, le plan de mise au point est le plan focal 5 objet du système optique, et le plan image est renvoyé à l'infini. Une optique de proximité, par exemple une lentille tube, permet de focaliser le plan image sur le capteur 132. Dans la suite du texte, le terme capteur 132 regroupe le capteur 132 et son optique de proximité. Le système optique 131 présente un axe optique 135. L'axe optique est 10 orthogonal au plan de mise au point et au plan image. L'axe optique 135 définit la profondeur de l'échantillon 111. Cet axe relie les interfaces inférieure et supérieure 116, 117 de l'échantillon 111. De préférence, l'axe optique 135 est sensiblement orthogonal aux interfaces inférieure et supérieure 116, 117. Le terme interface désigne une frontière de l'échantillon 111 étudié. 15 L'analyse au sens de l'invention comprend en particulier une détermination de la position de la particule d'intérêt 112A selon un axe parallèle à l'axe optique 135. Le dispositif 100 comprend en outre des moyens de translation 140, adaptés à déplacer l'ensemble d'imagerie 130 selon un axe parallèle à l'axe optique 135. 20 Dans une configuration plan-infini telle que définie ci-avant, les moyens de translation 140 peuvent ne déplacer que le système optique 131, le capteur 132 restant fixe. En variante, les moyens de translation 140 sont adaptés à déplacer le support 110 selon un axe parallèle à l'axe optique 135. Selon un mode de réalisation, les moyens de translation 140 réalisent en 25 outre des translations selon deux autres axes définissant un plan orthogonal à l'axe optique 135. Le capteur 132 est relié à des moyens de calcul 150, en particulier un processeur ou un microprocesseur. 3028616 10 Le dispositif 100 est adapté à la mise en oeuvre du procédé selon l'invention, illustré de façon schématique en figure 2. Le procédé comprend une première étape 21, de détermination d'un point de repère 118. Cette étape consiste à choisir, ou définir le point de repère 118. Ce 5 choix est généralement arbitraire. Le point de repère est situé sur une première interface de l'échantillon 111, ou à une distance connue de cette première interface, cette distance étant définie selon l'axe optique 135. En particulier, le point de repère est situé sur l'interface supérieure 117 ou 10 inférieure 116. Dans l'exemple représenté en figure 1, le point de repère est situé sur l'interface supérieure 117 de l'échantillon. Selon une variante non représentée, le point de repère ne se trouve pas directement sur une interface de l'échantillon, mais à une distance connue de celle-ci. Par exemple, le point de repère se trouve sur une face de la lame 15 supérieure 114 opposée à l'échantillon, l'épaisseur de la lame supérieure 114 étant connue. Dans le dispositif 100, les moyens de calcul 150 sont reliés à une mémoire 151 stockant une information relative au point de repère, en particulier la distance selon un axe parallèle à l'axe optique 135, entre ce point de repère et ladite 20 première interface et éventuellement ses coordonnées dans un plan orthogonal à l'axe optique 135. On illumine ensuite l'échantillon 111 à l'aide de la première source lumineuse 120, formant ainsi dans l'échantillon une région illuminée 119 telle que décrite ci-avant (étape 22). 25 Ensuite, on acquiert une image en transmission de la région illuminée 119, à l'aide de l'ensemble d'imagerie 130 (étape 23). Cette image est nommée image de référence. Lors de cette acquisition, le plan de mise au point 133 est situé à une distance connue D1 du point d'intérêt 118, cette distance étant définie selon un axe parallèle à l'axe optique 135. La particule d'intérêt 112A se trouve hors du plan 3028616 11 de mise au point 133. De préférence, lors de l'acquisition de l'image de référence, le plan de mise au point se trouve en dehors de l'échantillon 111. Par exemple, lorsque le point de repère 118 se trouve sur l'interface supérieure 117, le plan de mise au point est situé entre 5 um et 1500 um au-dessus ou en dessous (en amont 5 ou en aval) du point de repère, de préférence entre 5 um et 1000 um, voire entre 5 um et 800 um. Là encore, la distance est définie selon un axe parallèle à l'axe optique 135. D'une façon générale, lors de l'acquisition de l'image de référence, le plan de mise au point est situé à une distance inférieure à 2 um, et de préférence 10 supérieure à 5 um, de préférence dans la plage [5 um - 1500 um] voire [5um - 1000 um], voire [5um - 800um] d'une interface délimitant l'échantillon 111. En variante, on dispose d'une information a priori quant à la position d'une particule. Le plan de mise au point est alors décalé, par rapport à cet a priori, d'une valeur telle que définie ci-dessus. 15 Dans l'exemple représenté sur la figure 1, un exemple d'information a priori est l'observation de particules biologiques adhérentes à une surface, par exemple la lame supérieure 114 délimitant l'échantillon 111. La connaissance de la position de la lame permet d'établir un a priori sur la position des particules, ces dernières étant adhérentes à la face 117 de cette lame adjacente de l'échantillon. 20 On dit que l'image de référence est l'image d'une partie de la région illuminée, car elle est formée par des rayons lumineux provenant de cette région illuminée. L'image de référence est en particulier un hologramme, formé par l'interférence de rayons lumineux provenant de la première source lumineuse et diffractés par une particule biologique de l'échantillon, avec des rayons lumineux 25 provenant de la première source lumineuse 120 et ayant traversé l'échantillon sans être diffractés. L'image de référence comprend notamment un hologramme associé à la particule d'intérêt. Dans le dispositif 100, les moyens de calcul 150 sont reliés au capteur 132 pour recevoir l'image de référence. 3028616 12 Dans une étape 24, on met en oeuvre une construction numérique d'une série d'images reconstruites. Chaque image reconstruite correspond à un décalage prédéterminé, le long de l'axe optique 135, du plan de mise au point relativement à la position du plan de mise au point associée à l'image de référence. Il s'agit de 5 décalages virtuels, c'est-à-dire simulés par la construction numérique, et non pas mis en oeuvre physiquement. On parle généralement de propagation numérique. Ces décalages sont notés Di', Di2, ... Din, Di_i, ... Di_n. De préférence, les décalages s'étendent entre deux positions extrêmes vis-à-vis du plan de mise au point, ces deux positions extrêmes encadrant l'interface supérieure 117 et, de façon plus 10 générale, le point d'intérêt. Ces décalages peuvent être répartis selon un pas régulier, notamment un pas compris entre 0,10 um et 1 um, par exemple 0,25 um ou 0,20 um. Ainsi, on reconstruit numériquement les images qui seraient acquises par le capteur 132 si le plan de mise au point se trouvait successivement aux distances Di', Di2, ... Din, Di_i, ... Di_n de la position du plan de mise au point associée 15 à l'image de référence. L'image de référence est une image acquise expérimentalement, tandis que les images reconstruites forment des images construites par propagation numériques. La différence entre les deux positions extrêmes associées à ces décalages dépend de l'échantillon 111 observé, et notamment de son épaisseur. Lorsqu'on 20 ne dispose d'aucun a priori quant à la position des particules d'intérêt, les deux positions extrêmes sont déterminées de telle sorte qu'elles se situent de part et d'autre de l'échantillon 111. Lorsqu'on dispose d'un a priori sur la position d'une particule d'intérêt dans l'échantillon, les deux positions extrêmes sont établies pour être agencées de part et d'autre de cette position. 25 La reconstruction numérique utilise la transformée de Fourier de l'image de référence (ayant éventuellement subi un traitement préalable), à laquelle on applique un opérateur de propagation avant de la repasser dans l'espace réel. L'opérateur de propagation est une fonction du décalage associé à l'image reconstruite calculée. 3028616 13 Selon l'invention, l'opérateur de propagation est par exemple une intégrale basée sur l'équation de Rayleigh Sommerfeld. Dans l'article « 3D Localization of weak scatterers in digital holographic microscopy using Rayleigh-Sommerfeld back-propagation », 16 July 2012 / Vol. 20, 5 No 15 / OPTICS EXPRESS, 16735-16744, Wilson et al. décrit un exemple d'une propagation numérique mise en oeuvre à partir de l'image d'une microsphère. Un autre exemple de propagation numérique est décrit par Lee et al. dans l'article « Holographic microscopy of holographically trapped three-dimensinal structures », 19 February 2007 / Vol. 15, No. 4 / OPTICS EXPRESS 1505-1512. 10 Les images reconstruites sont des images complexes, c'est-à-dire avec une partie réelle et une partie imaginaire. Elles comprennent un ensemble de points, chaque point de l'image étant affecté d'une grandeur complexe. Les moyens de calcul 150 du dispositif 100 selon l'invention sont adaptés à réaliser, à partir de l'image de référence une construction numérique d'une série 15 d'images reconstruites, telle que décrite ci-dessus. On détermine ensuite, à partir de la série d'images reconstruites, la distance entre la particule d'intérêt 112A et le point de repère 118, cette distance étant définie selon un axe parallèle à l'axe 135 (étape 25). En variante, on se limite, à partir de la série d'images reconstruites, à la distance Dref entre la particule 20 d'intérêt 112A et le plan de mise au point 133. Pour cela, on calcule la distance Dref selon l'axe optique 135, entre la particule d'intérêt 112A, et le plan de mise au point associé à l'image de référence. Connaissant la distance D1 selon l'axe optique 135 entre ce plan de mise au point et le point de repère 118, on en déduit la distance DE entre la particule d'intérêt 25 112A et le point de repère 118 selon l'axe optique 135. Pour déterminer la distance Dref selon l'axe optique 135, entre la particule d'intérêt 112A, et le plan de mise au point associé à l'image de référence, on associe chaque image reconstruite à une valeur d'un paramètre utile, le paramètre utile étant fonction d'un paramètre de l'intensité complexe (ou amplitude 3028616 14 complexe) des points appartenant aux images reconstruites. Un paramètre de l'intensité complexe (ou amplitude complexe) d'une image reconstruite est par exemple la partie imaginaire, la partie réelle, le module ou la phase de l'image reconstruite. De préférence, ces valeurs du paramètre utile sont réunies sous la 5 forme d'un profil, représentant la valeur du paramètre utile en fonction du décalage de chaque image reconstruite, selon l'axe optique 135 par rapport au plan de mise au point 133. Ensuite, on recherche une valeur remarquable sur ce profil, par exemple un maximum ou un point d'inflexion ou un passage par zéro. Cette valeur 10 remarquable est associée au décalage selon l'axe optique 135, entre la particule d'intérêt et la position du plan de mise au point 133 associé à l'image de référence. De préférence, on fait en sorte que la valeur remarquable à rechercher soit un maximum et de façon générale un extremum, pour limiter l'effet de l'incertitude lors du passage par une dérivée seconde (détection d'un point d'inflexion), et pour 15 s'affranchir d'une correction des biais de mesure (détection d'un passage par zéro). Mais la valeur remarquable du paramètre utile peut également être un passage par zéro de la fonction décrivant l'évolution de ce paramètre en fonction du décalage entre l'image reconstruite et le plan de mise au point 133, ou un point d'inflexion de cette fonction, ou tout autre critère. 20 Ainsi, de façon générale, on détermine la distance Dref entre le plan de mise au point 133 de l'image de référence et une particule d'intérêt 112A, par la succession des étapes suivantes : - identification, sur l'image de référence, d'un hologramme associé à ladite particule d'intérêt, 25 reconstruction numérique d'images, en appliquant un opérateur de projection à ladite image de référence, ou à une région d'intérêt de l'image de référence comportant ledit hologramme, de façon à obtenir, pour une pluralité de décalages Di', Die, ... Din, Dig, ... DL, prédéterminés par rapport au plan de mise au 3028616 15 point 133, une image reconstruite, chaque image reconstruite correspondant à chaque décalage, extraction, sur chaque image reconstruite, d'un paramètre représentatif de la valeur complexe de chaque point de l'image, dit paramètre 5 utile, analyse de l'évolution de ce paramètre utile en fonction dudit décalage et identification d'une valeur remarquable de ce paramètre utile, détermination du décalage correspondant à ladite valeur remarquable de ce paramètre utile, la distance Dref entre la particule d'intérêt 112A et le plan de 10 mise au point 133 étant alors considérée comme égale à ce décalage. Connaissant la distance D1 entre le plan de mise au point 133 et le point de repère 118, il est alors possible de déduire la distance DE entre la particule d'intérêt 112A et le point de repère 118, par l'opération DE = Dref- D1 15 Par extraction, sur chaque image reconstruite, d'un paramètre utile, on entend le calcul d'un paramètre représentatif des différents points constituant l'image, par exemple en effectuant une somme ou une moyenne des grandeurs complexes associées à chacun de ces points. Le paramètre utile peut être le carré de la partie imaginaire de l'image, 20 déterminé à partir de la valeur de la partie imaginaire des différents points constituant l'image. Le paramètre utile correspond alors à la moyenne du carré de la partie imaginaire de l'ensemble des points constituant l'image. Les inventeurs ont constaté qu'un tel paramètre permettait d'obtenir une bonne précision de localisation. 25 Outre la partie imaginaire, le paramètre peut également comprendre la valeur de la partie réelle des points constituant l'image ou leur module, ou le carré de ces différentes grandeurs. Les inventeurs ont montré qu'il est possible de localiser la position, selon l'axe optique 135, d'une particule d'intérêt 112A, cette 3028616 16 position correspondant à une valeur remarquable du paramètre utile, et par exemple : un maximum : c'est le cas lorsqu'on observe l'évolution, selon l'axe optique 135, du carré de la partie imaginaire, ou du module, ou du module au carré 5 de l'image. un passage par zéro, par exemple lorsqu'on observe l'évolution, selon l'axe optique 135, de la partie réelle de l'image ou de son carré. De façon alternative ou complémentaire, l'identification d'une valeur 10 remarquable dudit paramètre d'intérêt permet de conclure à la présence d'une particule d'intérêt 112A dans le milieu, dans la plage de décalage considérée, la distance entre la particule d'intérêt 112A et le plan de mise au point Dref n'étant pas forcément mémorisée. A chaque particule d'intérêt 112A correspond un hologramme sur l'image de 15 référence. Il est alors possible de sélectionner une région d'intérêt se limitant à un hologramme, et d'appliquer les étapes précédemment décrites à chaque région d'intérêt associé à un hologramme. Les moyens de calcul 150 fournissent en sortie la distance DE, selon l'axe 135, entre la particule d'intérêt 112A et le point de repère 118. 20 Le recours à une propagation numérique permet de réduire le nombre d'images de l'échantillon à acquérir. En particulier, l'acquisition d'une unique image est suffisante. Le procédé selon l'invention est donc particulièrement rapide et automatisa ble. Dans certains cas, lorsque la particule d'intérêt se trouve dans le plan de mise 25 au point du système optique, son image ne peut pas être distinguée de l'image du milieu environnant (le milieu environnant correspondant à la portion d'échantillon 111 entourant la particule d'intérêt). C'est par exemple le cas lorsque la particule d'intérêt baigne dans un milieu de même coefficient de transmission ou d'indice de réfraction voisin (par exemple à 20% près), dans le spectre de la première source 3028616 17 lumineuse. Dans ce cas, le procédé selon l'art antérieur tel que décrit en introduction ne permet pas de déterminer la position de la particule d'intérêt. En revanche, on pourra identifier la particule d'intérêt sur une image reconstruite selon l'invention, car cette image reconstruite est une image complexe contenant 5 plus d'information qu'une image acquise. Ainsi, le procédé selon l'invention permet de déterminer la position d'une particule selon un axe parallèle à l'axe optique du système optique, même lorsque la particule ne peut être distinguée de son milieu environnant sur une image acquise lorsque la particule se trouve dans le plan de mise au point du système optique. 10 Enfin, le procédé selon l'invention permet d'associer une image et une position du plan de mise au point selon l'axe optique 135, cette position étant déterminée avec une grande précision (liée au pas des décalages mis en oeuvre à l'étape 23). Une telle précision ne serait pas facilement atteignable à l'aide de décalages physiques, et non virtuels. 15 L'image de référence peut comprendre l'hologramme d'une unique particule biologique, qui définit alors la particule d'intérêt selon l'invention. En variante, l'image de référence comprend l'hologramme de plusieurs particules biologiques, on choisit arbitrairement l'un de ces hologrammes pour définir la particule d'intérêt selon l'invention, et on utilise des portions des images 20 reconstruites, ou régions d'intérêt, centrées sur l'un de ces hologrammes. On peut mettre en oeuvre une étape préliminaire de normalisation de l'image de référence. Pour cela, on acquiert une image dite de fond, représentative des défauts du système optique (réflexions internes, poussières, etc.). L'image de fond correspond par exemple à l'image acquise en déplaçant l'échantillon dans un plan 25 parallèle au plan de mise au point, pendant la prise d'image. En variante, l'image de fond est une image moyenne formée à partir de plusieurs images statiques obtenues pour plusieurs positions de l'échantillon dans un plan parallèle au plan de mise au point. Ensuite, on divise, pixel par pixel, l'image de référence par l'image 3028616 18 de fond. On obtient ainsi une image de référence normalisée Une telle normalisation améliore les résultats de la méthode. Selon une première variante de l'invention, on réalise la propagation numérique à partir de l'image de référence utilisée directement, ou à partir de 5 l'image de référence normalisée telle que définie ci-dessus. L'image reconstruite associée à un décalage z est notée : Up (z) = TF-1{TF(U0) X H(z)} avec U0 l'image utilisée pour réaliser la propagation numérique (ici l'image de référence ou l'image de référence normalisée), TF l'opérateur transformée de 10 Fourier, TF' l'opérateur transformée de Fourier inverse, X la multiplication terme à terme de matrices, et H un opérateur de propagation basé sur l'intégrale de Rayleigh Sommerfeld. Dans l'exemple présenté ici, on a en particulier : H(u,v,z) = exp[-Izi* Im(p(u,v)) + i * z * Re (p(u,v))1 15 avec - Re l'opérateur partie réelle, lm l'opérateur partie imaginaire, j.2 = -1, - u et y les coordonnées dans l'espace de Fourier, associées aux coordonnées x, y un plan orthogonal à l'axe 135, dans l'espace réel, - p(u,y)=kx A x(I 1- u2 -y2 -1), et k=271. 20 Les coordonnées u et y sont définies de la façon suivante : 2 7 0 1 . 2 1 ï . . 2 7 Aix-1 2, 2 u=- --- 2, ,- -- ,... , 2) parpasde 1 . 2 2, 2 1_ i -2) Ap /^Tx ApAix 41, Nx (Ny-1 1 2 / 0 1 2 Nyx,Ap y =- --- parpasde A N ' Ai, N y Ai, N y N xA P Y Y P et avec - Op la taille du pas d'échantillonnage, défini dans le plan de mise au point (donc lié au pas de pixel du capteur 132 et au facteur de grandissement du système optique 131), 3028616 19 ^ N, et Ny les nombres de pixels selon x, respectivement y, dans l'image Uo, ^ la longueur d'onde centrale de la première source lumineuse. 5 On détermine ensuite la position de la particule d'intérêt selon l'axe optique du système optique, à partir de la partie imaginaire de chaque image reconstruite Up(z). En particulier, on recherche l'image reconstruite ou une portion d'image reconstruite, dont le carré de la partie imaginaire (c'est-à-dire la valeur moyenne du carré de la partie imaginaire des différents points composant l'image) est 10 maximum. Une portion d'image reconstruite est une partie d'une image reconstruite, centrée sur l'image de la particule d'intérêt. Si l'image de référence se rapporte à plusieurs particules biologiques, on peut utiliser la même image de référence pour déterminer la position de chacune de ces particules biologiques, en utilisant à chaque fois des images reconstruites centrées sur une particule 15 biologique différente. On va maintenant décrire, en référence aux figures 3 et 4, un deuxième mode de réalisation de procédé et de dispositif 300 selon l'invention. La figure 3 ne sera décrite que pour ses différences relativement à la figure 1. Les références numériques 111, 120, 112A, 116, 117, 118, 130, 131, 132, 135, 20 120, 150, 151 de la figure 1 correspondent respectivement aux références numériques 311, 320, 312A, 316, 317, 318, 330, 331, 332, 335, 320, 350, 351 de la figure 3. Chacune des caractéristiques distinguant le dispositif selon la figure 3 du dispositif selon la figure 1 pourra être isolée et combinée indépendamment au 25 dispositif de la figure 1 pour former de nombreuses autres variantes de l'invention. Selon un mode de réalisation avantageux tel que représenté en figure 3, les particules biologiques adhèrent directement sur une deuxième interface de l'échantillon. 3028616 20 La droite parallèle à l'axe 335 et passant par le point de repère 318 traverse cette deuxième interface. Cette deuxième interface est avantageusement confondue avec la première interface telle que définie ci-avant. Dans le cas contraire, on connaît de préférence 5 la position de la première interface relativement à la deuxième interface, selon un axe parallèle à l'axe optique 135. En particulier, cette position s'étendant selon un axe parallèle à l'axe optique 135 et passant par le point de repère 318. On dispose ainsi d'une connaissance a priori sur la position, selon l'axe optique 335, du point de repère 318 relativement aux particules biologiques. 10 La deuxième interface est sensiblement orthogonale à l'axe 335. Il s'agit en particulier de l'interface inférieure ou supérieure 316, 317 de l'échantillon. En pratique, on étudie par exemple : - une gélose, les particules biologiques s'étendant sur la face supérieure de celle-ci, ou 15 - des particules biologiques qui adhèrent naturellement sur une lame recouvrant l'échantillon, ladite lame étant située au-dessus de l'échantillon. Dans l'exemple représenté en figure 3, les particules biologiques adhèrent sur l'interface supérieure 317. 20 L'interface supérieure est légèrement inclinée relativement à un plan orthogonal à l'axe 335, typiquement d'un angle inférieur à 0,1 rad, voire même 0,05 rad ou 0,02 rad. Le procédé selon l'invention permet de s'affranchir de l'incertitude sur le positionnement des particules biologiques selon l'axe 335, lorsque la lame 25 supérieure, définissant l'interface supérieure de l'échantillon, est inclinée relativement à un plan orthogonal à l'axe 335. Plus généralement, le procédé selon l'invention permet de s'affranchir de l'incertitude sur le positionnement des particules biologiques selon l'axe 335, lorsque la lame supérieure, définissant l'interface supérieure de l'échantillon, est 3028616 21 déformée relativement à un plan orthogonal à l'axe 335. Cette déformation peut désigner toute déformation telle que les particules biologiques adhérant sur ladite deuxième interface soient situées à l'intérieur d'un cylindre délimité par deux surfaces planes perpendiculaires à l'axe optique 335 et distantes de moins de 5 100 um selon cet axe, voire moins de 50 um. De préférence, la distance selon un axe parallèle à l'axe optique 335, entre la position du plan de mise au point associée à l'image de référence, et la projection du point de repère 318 selon cet axe et sur la deuxième interface, est comprise entre +5 um et +1500 um ou entre -5 um et -1500 um. Avantageusement, cette 10 distance est comprise entre + 5 um et +1000 um ou - 5 um et -1000 um, voire entre + 5 um et +800 um ou - 5 um et -800 um, ou même entre + 5 um et +200 um ou 5 um et -200 um. On définit ainsi des plages de défocalisation optimale, offrant un calcul plus aisé de la distance entre la particule d'intérêt et le point de repère. Puisque les particules biologiques adhèrent sur ladite deuxième interface et 15 sont situées dans un cylindre de hauteur limitée, on s'assure facilement que lors de l'acquisition de l'image de référence, la particule d'intérêt se trouve à une distance du plan de mise au point située dans une plage de défocalisation optimale. La largeur des plages de défocalisation optimale peut dépendre des caractéristiques de la première source lumineuse 320, notamment sa cohérence 20 spatiale et temporelle. Ainsi, on pourra adapter la première source lumineuse en fonction d'une distance selon l'axe optique 335, entre plusieurs particules d'intérêt que l'on souhaite étudier, de sorte que ces particules d'intérêt se trouvent simultanément dans une plage de défocalisation optimale. 25 Par exemple, on peut utiliser une source de lumière blanche, pour une plage de défocalisation optimale allant + 5 um à +200 um, et de -5 um à -200 um. Plus on restreint la largeur spectrale de la première source lumineuse, plus la largeur de la plage de défocalisation optimale augmente. 3028616 22 En outre, plus le diamètre apparent de la première source lumineuse 320 est faible, plus la plage de défocalisation optimale est large. Dans l'exemple représenté en figure 3, le dispositif d'analyse 300 comprend des moyens pour déterminer la distance, selon un axe parallèle à l'axe optique 335, 5 entre le point de repère et la position du plan de mise au point du système optique associée à l'image de référence. Cette distance est nommée distance utile. Ces moyens comprennent en particulier un laser 360, adapté à illuminer le point de repère 318. Un obturateur 361 relié à une platine de translation 362 permet d'obturer la sortie du laser, lorsque l'on ne souhaite pas que le faisceau 10 laser 363 illumine le point de repère 318. Le faisceau laser 363 est avantageusement incident sur le point de repère, selon un axe parallèle à l'axe optique 335. Le faisceau laser 363 et le capteur d'image 332 sont optiquement couplés au même système optique 331, de telle sorte que le point de focalisation du faisceau 15 laser correspond au plan objet du capteur, ou plan de mise au point. Dans l'exemple représenté en figure 3, le point de repère se trouve sur l'interface supérieure 317, c'est pourquoi le faisceau laser 363 pénètre dans l'échantillon par cette interface supérieure 317. Le faisceau laser 363 traverse une première lame séparatrice 371, puis est 20 réfléchi sur une deuxième lame séparatrice 371, avant d'atteindre l'interface supérieure 317. Chaque lame séparatrice peut être une lame dichroïque. En variante, on utilise un cube ou un miroir semi-réfléchissant. La distance utile est déterminée en mettant en oeuvre les sous-étapes suivantes illustrées en figure 4. 25 Dans une sous-étape 41, on illumine le point de repère 318 à l'aide du faisceau laser. Ensuite, on ajuste la distance, notamment selon un axe parallèle à l'axe optique 335, entre l'échantillon, et le système optique 331, de sorte que le capteur 332 reçoive l'image d'une tache formée par la réflexion spéculaire du faisceau laser 3028616 23 363 sur le point de repère (sous-étape 42). Cet ajustement peut être réalisé grâce aux moyens de translation 340 décrits en référence à la figure 1. En pratique, la définition du point de repère 318 peut dépendre de la position arbitraire du faisceau laser 363, dans un plan orthogonal à l'axe optique 335. 5 Cet ajustement met en oeuvre un repérage de la réflexion spéculaire du faisceau laser 363 sur la surface recevant le point de repère. Un exemple d'un tel ajustement est illustré sur les figures 5A à 5C, le point de repère étant situé sur l'interface supérieure 317. L'image 5A correspond à une position du plan de mise au point 533 au-dessus 10 de l'interface supérieure 317. L'image obtenue sur le capteur est une tache large et peu lumineuse. L'image 5B correspond à une position du plan de mise au point sur l'interface supérieure 317. L'image obtenue sur le capteur est une tache étroite et très lumineuse. 15 L'image 5C correspond à une position du plan de mise au point au-dessous de l'interface supérieure 317. L'image obtenue sur le capteur est une tache large et peu lumineuse. En déplaçant le système optique selon l'axe 335, on observe successivement : une tache large et peu lumineuse ; 20 une première tache étroite et très intense correspondant à la réflexion du faisceau laser 363 focalisé par le système optique 331 sur la face 315 de la lame supérieure 314, opposée à l'échantillon, cette tache présentant une intensité maximale ; une deuxième tache étroite et intense correspondant à la réflexion du 25 faisceau laser 363 focalisé par le système optique 331 sur la face 317 de la lame supérieure 314, cette face étant adjacente de l'échantillon, et constituant alors une de ses interfaces, cette tache présentant un maximum secondaire de l'intensité, puis 3028616 24 une tache large et peu lumineuse, correspondant à la rétrodiffusion du faisceau laser dans l'échantillon 311. Ainsi, l'analyse de l'évolution de l'intensité lumineuse du signal de réflexion 5 spéculaire du faisceau laser, en fonction de l'écartement relatif entre le système optique 331 et l'échantillon 311 permet de déterminer la position du point d'intérêt 318, ce dernier correspondant, dans cet exemple, au point de focalisation du faisceau laser avec une surface apte à réfléchir le faisceau, en l'occurrence la face 317 ou la face 315. 10 Dans l'exemple précédemment décrit, l'image de la réflexion spéculaire la plus intense et l'image de la réflexion spéculaire présentant un maximum secondaire d'intensité, correspondent aux configurations dans lesquelles le plan objet passe par l'intersection du faisceau laser et, respectivement, la face de la lame supérieure 314 opposée à l'échantillon, et la face de la lame supérieure 314 15 adjacente de l'échantillon. On peut remarquer qu'il est particulièrement intéressant de repérer la réflexion du faisceau laser sur une face de la lame supérieure adjacente de l'échantillon, car elle sera visible même lorsque le système optique est formé par un objectif à immersion. 20 Ensuite, le procédé peut comprendre une sous-étape 43, consistant à déplacer d'une distance connue et selon l'axe optique 335, le support recevant l'échantillon relativement aux moyens d'imagerie. Là encore, on peut utiliser les moyens de translation tels que décrits en référence à la figure 1. Cette sous-étape peut être utile pour placer une particule d'intérêt dans une plage de défocalisation 25 optimale telle que définie ci-avant. On a représenté en figure 3 le faisceau laser 364 réfléchi sur l'interface supérieure 317, et imagé sur le capteur par le système optique 331. Le faisceau laser réfléchi 364 est réfléchi sur la deuxième lame séparatrice 372, puis sur la première lame séparatrice 371. 3028616 25 Les moyens de calcul 351 reçoivent en entrée une image acquise par le capteur 351, lorsque l'obturateur 361 est ouvert, et que la première source lumineuse 320 est éteinte. Ils pilotent les moyens de translation 340 pour mettre en oeuvre la sous-étape 42, et le cas échéant la sous-étape 43. 5 De préférence, le procédé selon l'invention comprend également une détermination, à partir de la série d'images reconstruites, de la position de la particule d'intérêt dans un plan orthogonal à l'axe optique du système optique. Un exemple d'une telle détermination sera décrit plus avant en référence à la figure 7. 10 Les moyens de translation 340 peuvent mettre en oeuvre une translation selon deux axes définissant ensemble un plan orthogonal à l'axe optique 335. On peut ainsi localiser successivement plusieurs particules d'intérêt d'un même échantillon. Après une unique étape de détermination d'un point de repère, on met en oeuvre plusieurs fois les séries d'étapes suivantes telles que décrites ci- 15 avant : illumination d'une région de l'échantillon ; acquisition d'une image de référence ; construction numérique d'une série d'images reconstruites ; et détermination de la distance entre une particule d'intérêt et le point 20 de repère, et/ou détection du nombre et de la position à partir de l'image de référence ou d'une ou plusieurs images reconstruites des particules d'intérêt. Entre deux séries de ces étapes, on translate l'échantillon dans un plan orthogonal à l'axe optique 335. 25 La figure 3 permet également d'illustrer une étape supplémentaire d'analyse mise en oeuvre dans un procédé et un dispositif selon l'invention. Le dispositif 300 comprend des moyens pour positionner le waist du faisceau laser 363 (ou col, c'est-à-dire emplacement du faisceau laser dont le diamètre est le plus étroit) sur la particule d'intérêt. 3028616 26 Le rayonnement émis par la particule d'intérêt, en réaction, est collecté par le système optique 331 et reçu par un spectromètre 380 qui analyse ce rayonnement. Le spectromètre 380 est par exemple un spectromètre analysant sur la diffusion Raman ou un spectromètre de fluorescence. Dans ce type d'analyse, il 5 est préférable que le faisceau laser d'excitation soit centré sur la particule examinée. Cela évite une perturbation du spectre de diffusion Raman par le voisinage de la particule. Le positionnement du waist du faisceau laser 363 est mis en oeuvre grâce à des moyens de translation tels qu'une platine de translation (non représentée), 10 adaptés à déplacer l'échantillon relativement aux moyens d'imagerie 331. Ces moyens de translation peuvent être formés par les moyens de translation 340. Les moyens de translation sont pilotés par les moyens de calcul 350. Les moyens de calcul 350 utilisent la position de la particule d'intérêt relativement au point de repère pour positionner le plan de mise au point du système optique 331 sur la 15 particule d'intérêt. On pourra envisager tout autre type d'analyse nécessitant un positionnement précis du faisceau laser sur la particule d'intérêt. On pourra utiliser une source laser distincte du laser 360 pour analyser la particule d'intérêt. De la même façon, un photodétecteur 332, optiquement couplé au système 20 optique 331, peut être utilisé, par exemple pour recueillir un signal de fluorescence émis par la particule d'intérêt 112A en réponse à un signal d'excitation. La connaissance d'une des distances précédemment évoquées (DE ou Dref) permet la focalisation du système optique sur la particule, optimisant la collecte du signal de fluorescence par le photodétecteur 332. 25 Dans la plupart des cas, un milieu comporte une pluralité de particules d'intérêt. Durant une première phase, on détermine et on mémorise la distance Dref ou DE correspondant à chaque particule d'intérêt 112A de l'échantillon. Durant une deuxième phase, le positionnement relatif du système optique 331 est ajusté de telle sorte que le plan de mise au point de ce système optique comprenne la 3028616 27 particule d'intérêt. Cet ajustement est réalisé, successivement, pour chaque particule d'intérêt, de manière à optimiser l'analyse. La figure 6 illustre un exemple d'une image de référence selon l'invention. Les axes des abscisses et des ordonnées sont gradués en pixels. On distingue 5 plusieurs figures de diffraction 61, associées chacune à une particule biologique. On va maintenant décrire, en référence à la figure 7, un exemple détaillé des étapes mises en oeuvre pour déterminer, à partir d'une première série d'images reconstruites, la distance entre la particule d'intérêt et le point de repère, selon un axe parallèle à l'axe optique du système optique. 10 Dans l'exemple représenté, on se place dans le mode de réalisation tel que décrit ci-avant, dans lequel on a : H (u, v, z) = exp[-lz I * Int(p(u, v)) + i * z * Re (p (u, 0)1 On met ensuite en oeuvre les sous-étapes suivantes : 15 - sous-étape 71 : On calcule, pour chaque image reconstruite de la première série d'images reconstruites, le carré de la partie imaginaire. On obtient une première série d'images reconstruites dites utiles. Le carré de la partie imaginaire d'une image reconstruite définit donc l'intensité de l'image reconstruite utile correspondante. 20 - sous-étape 72 : Chaque image reconstruite utile correspond à une même matrice de pixels. Selon l'image reconstruite utile, l'intensité associée à chaque pixel varie. Pour chaque pixel, on sélectionne l'intensité maximale parmi les intensités sur les différentes images reconstruites utiles. On forme ainsi l'image des maxima. 25 De façon similaire, on forme l'image des minima. On calcule ensuite une image des gradients 720, correspondant à la différence entre l'image des maxima et l'image des minima. On peut réaliser cette étape à l'aide d'une première matrice, dont une dimension correspond aux pixels, et l'autre dimension correspond au décalage du 3028616 28 plan de mise au point associé à l'image reconstruite utile. Ce décalage est mesuré selon l'axe optique du système optique, nommé axe (Oz) dans la suite. Par exemple, si les images reconstruites utiles présentent chacune X pixels en largeur et Y pixels en hauteur, et qu'on reconstruit N images reconstruites utiles, 5 la première matrice présente XY colonnes et N lignes. L'image des minima correspond alors à une deuxième matrice XY colonnes et 1 ligne. L'image des maxima correspond à une troisième matrice XY colonnes et 1 ligne. L'image des gradients correspond une quatrième matrice XY colonnes et 1 ligne, égale à la différence entre la troisième et la deuxième matrices. 10 - sous-étape 73 : On sélectionne les pixels présentant un fort gradient d'intensité selon l'axe (Oz). On définit un gradient moyen associé à ces pixels. La position des plus grandes valeurs de ce gradient moyen définit grossièrement la position des particules biologiques selon l'axe (Oz). 15 Pour ce faire, on peut sélectionner dans la quatrième matrice, les colonnes présentant les valeurs les plus élevées (par exemple premier percentile). On définit ainsi une série de M pixels. Ensuite, on sélectionne ces pixels dans la première matrice, et on combine (par exemple somme ou moyenne) les pixels sélectionnés associés à un même décalage selon l'axe (Oz). On forme ainsi une cinquième 20 matrice N lignes et 1 colonne. On recherche ensuite le maximum parmi les valeurs de cette cinquième matrice, et on relève le décalage selon l'axe (Oz), associé à ce maximum. La figure 730 illustre la cinquième matrice sous la forme d'un profil dans lequel l'axe des abscisses est un décalage selon l'axe (Oz), et l'axe des ordonnées 25 est la valeur correspondante dans la cinquième matrice. La sous-étape 73 définit une position moyenne selon l'axe (Oz), de plusieurs particules biologiques imagées sur l'image de référence. Cette position moyenne définit une position approchée d'une particule d'intérêt selon l'invention, en 3028616 29 supposant que les particules biologiques sont situées sensiblement dans un même plan orthogonal à l'axe (Oz). On exploite un gradient selon l'axe (Oz) associé à une sélection de pixels présentant chacun un fort gradient selon cet axe. Le maximum du gradient associé 5 à cette sélection est donc d'autant plus facile à repérer. - sous-étape 74 : On sélectionne parmi les images de la première série d'images reconstruites utiles, celle associée à la position selon l'axe (Oz), calculée à la sous-étape 73. Sur cette image, on réalise un seuillage, pour obtenir une image binaire 740. 10 Ensuite, on relève les coordonnées des objets binaires mis en évidence par le seuillage. On calcule ainsi une position approchée de chaque particule biologique, dans un plan orthogonal à l'axe (Oz). Cette étape peut comprendre une étape d'interpolation de chaque objet binaire par une ellipse, pour calculer en outre une forme géométrique approchée de chaque objet binaire, et donc de chaque 15 particule biologique. On choisit arbitrairement la particule d'intérêt en choisissant arbitrairement l'un des objets binaires. On a donc déterminé une position approchée de la particule d'intérêt, dans un plan orthogonal à l'axe (Oz). On définit sur l'image binaire 740 une région d'intérêt comprenant ce seul 20 objet binaire. La région d'intérêt est définie par des positions de pixels. La région d'intérêt est par exemple un carré de 16*16 pixels. - sous-étape 75 : On met en oeuvre cette étape à partir de la première série d'images reconstruites utiles, en sélectionnant dans chaque image reconstruite utile une 25 région correspondant à la région d'intérêt telle que définie ci-avant. En variante, on met en oeuvre cette étape à partir d'une deuxième série d'images reconstruites utiles. Chaque image de la deuxième série d'images reconstruites utiles est formée par le carré de la partie imaginaire d'une image d'une deuxième série d'images reconstruites. 3028616 La deuxième série d'images reconstruites est associée à un pas d'échantillonnage selon l'axe (Oz), inférieur au pas d'échantillonnage de la première série d'images reconstruites. Chaque image de la deuxième série d'images reconstruites correspond à la région d'intérêt telle que définie ci-avant. 5 L'amplitude de la plage d'échantillonnage selon l'axe (Oz) associée à la deuxième série d'images reconstruites, est inférieure à l'amplitude de la plage d'échantillonnage associée à la première série d'images reconstruites. La deuxième série d'images reconstruites peut être calculée uniquement pour la recherche d'une position précise de la particule d'intérêt. 10 On met en oeuvre les sous-étapes 72, et 73 sur la deuxième série d'images reconstruites utiles, pour calculer une position précise de la particule d'intérêt, selon l'axe (Oz). La figure 750 correspond à la figure 730, en remplaçant la première série d'images reconstruites utiles par la deuxième série d'images reconstruites utiles. 15 - sous-étape 76 : On calcule une position précise de la particule d'intérêt, dans un plan orthogonal à l'axe (Oz). Pour cela, on sélectionne parmi les images de la deuxième série d'images reconstruites utiles, celle associée à la position précise de la particule d'intérêt 20 selon l'axe (Oz). Sur cette image, on réalise un seuillage, pour obtenir une nouvelle image binaire. Ensuite, on relève les coordonnées du centre de l'objet binaire mis en évidence par le seuillage. De nombreuses variantes peuvent être envisagées. 25 En particulier, on peut se placer dans le mode de réalisation tel que décrit ci-avant, dans lequel on a : 3028616 31 - H(u, v, z) = exp(i *2 *if * z * w), avec w = T12 -2u _ v2, pour u2 +v2 REVENDICATIONS1. Procédé d'analyse d'un échantillon (111 ; 311) recevant des particules biologiques (112), parmi lesquelles une particule d'intérêt (112A), caractérisé en ce qu'il est mis en oeuvre dans un dispositif (100 ; 300) comprenant l'échantillon disposé entre une première source lumineuse (120 ; 320) et un système optique (131 ; 331), et en ce qu'il comprend les étapes suivantes : définition d'un point de repère (118 ; 318), situé sur une première interface (116, 117 ; 316, 317) de l'échantillon ou à une distance connue de celle-ci selon un axe parallèle à l'axe optique (135 ; 335) du système optique ; à l'aide de la première source lumineuse, illumination d'une région recevant la particule d'intérêt, dite région illuminée (119) ; à l'aide d'un capteur (132 ; 332) situé dans un plan image du système optique, acquisition d'une image de la région illuminée, dite image de référence, la particule d'intérêt (112A) se trouvant hors du plan de mise au point (133 ; 333) du système optique (131 ; 331), et la distance, selon un axe parallèle à l'axe optique du système optique, entre ledit plan de mise au point (133; 333) et le point de repère, étant une distance connue dite distance utile (D1) ; à partir de l'image de référence, construction numérique d'une série d'images reconstruites, associées chacune à un décalage prédéterminé (Di', Die, Din, Di_i, Di_n) du plan de mise au point le long de l'axe optique dudit système optique ; à partir de la série d'images reconstruites, détermination de la distance (Dref) selon un axe parallèle à l'axe optique du système optique, entre la particule d'intérêt et ledit plan de mise au point. 2. Procédé d'analyse selon la revendication 1, caractérisé en ce qu'il comporte une détermination de la distance (DE) entre la particule d'intérêt et ledit point de 3028616 36 repère (118 ; 318), à partir de ladite distance (Dref) entre la particule d'intérêt (112A) et ledit plan de mise au point (133 ; 333) ainsi que de ladite distance utile (D1) . 5 3. Procédé d'analyse selon la revendication 1 ou 2, caractérisé en ce ce qu'il comporte une détermination de la présence de ladite particule d'intérêt dans le milieu à partir de la distance entre ladite particule d'intérêt et ledit plan de mise au point (Dref). 10 4. Procédé d'analyse selon l'une des revendications 1 à 3, caractérisé en ce que la détermination de la distance utile (D1), comprend les sous-étapes suivantes : illumination du point de repère (118 ; 318) par un faisceau laser (363), ledit faisceau laser étant focalisé par ledit système optique (131; 331) ; 15 acquisition d'une image sur le capteur (132 ; 332) représentant la réflexion du faisceau laser sur une interface (116, 117 ; 316, 317) de l'échantillon (111 ; 311) ajustement de la distance utile (D1) selon un axe parallèle à l'axe optique du système optique, entre l'échantillon (111 ; 311) et le système optique 20 (131; 331), en fonction de l'intensité de l'image formée sur le capteur (132; 332). 5. Procédé d'analyse selon l'une des revendications 1 à 4, caractérisé en ce que les particules biologiques (112) adhèrent sur une deuxième interface (116, 25 117 ; 316, 317) de l'échantillon, la première et la deuxième interfaces étant distinctes ou confondues. 6. Procédé d'analyse selon l'une des revendications 1 à 5, caractérisé en ce que la distance (D1) selon un axe parallèle à l'axe optique du système optique, 3028616 37 entre la position du plan de mise au point associée à l'image de référence, et la projection du point de repère (118 ; 318) selon cet axe et sur la deuxième interface, est comprise entre +5 um et +2000 um ou entre -5 um et -2000 um. 5 7. Procédé d'analyse selon l'une quelconque des revendications 1 à 6, caractérisé en ce que la région illuminée (119) comprend une pluralité de particules biologiques (112). 8. Procédé d'analyse selon l'une quelconque des revendications 1 à 7, 10 caractérisé en ce que la première source lumineuse (120 ; 320) présente une largeur spectrale inférieure à 200 nm. 9. Procédé d'analyse selon l'une quelconque des revendications 1 à 8, caractérisé en ce qu'il comprend en outre une détermination, à partir de ladite 15 série d'images reconstruites, de la position de la particule d'intérêt (112A ; 312A) dans un plan orthogonal à l'axe optique du système optique. 10. Procédé d'analyse selon l'une quelconque des revendications 1 à 9, caractérisé en ce que chaque image reconstruite est formée par une partie réelle 20 et une partie imaginaire, et en ce que parmi les parties réelles et imaginaires, seules les parties imaginaires des images reconstruites sont utilisées pour déterminer la distance (DE) selon un axe parallèle à l'axe optique du système optique, entre la particule d'intérêt (112A; 312A) et le point de repère (118 ; 318). 25 11. Procédé d'analyse selon l'une des revendications 1 à 10, caractérisé en ce que chaque image reconstruite est associée à un décalage le long de l'axe optique (135 ; 335) du système optique et à une valeur d'un paramètre utile, de manière à constituer une fonction décrivant l'évolution du paramètre utile selon ledit décalage, et en ce que la détermination de la distance (Dref) entre la particule 3028616 38 (112A, 312A) et le plan de mise au point (133, 333) met en oeuvre une recherche d'une valeur remarquable, en particulier un extremum, un point d'inflexion ou un passage par zéro de ladite fonction. 5 12. Procédé d'analyse selon l'une quelconque des revendications 1 à 11, dans lequel que la détermination de la distance (Dref) entre la particule d'intérêt (112A; 312A) et le plan de mise au point (133 ; 333) comprend les sous-étapes suivantes : - à partir d'une première série d'images reconstruites associée à un premier pas des décalages du plan de mise au point (133 ; 333), détermination d'une 10 distance approchée entre la particule d'intérêt et ledit plan de mise au point ; à partir d'une deuxième série d'images reconstruites associée à un deuxième pas des décalages du plan de mise au point (133), le deuxième pas étant plus fin que le premier pas, détermination d'une distance précise 15 entre la particule d'intérêt et ledit plan de mise au point. 13. Procédé d'analyse selon l'une quelconque des revendications 1 à 12, caractérisé en ce qu'à partir de la distance (Dref) entre la particule d'intérêt (112A, 312) et ledit plan de mise au point (133, 333), le système optique (131, 331) est 20 déplacé par rapport à ladite particule d'intérêt (112A, 312A) de façon à focaliser le faisceau laser d'analyse (363) sur la particule d'intérêt (112A; 312A). 15. Procédé d'analyse selon l'une quelconque des revendications 1 à 14, caractérisé en ce qu'à partir de ladite distance (Dref) entre la particule d'intérêt 25 (112A, 312) et ledit plan de mise au point (133, 333), le système optique (131, 331) est déplacé par rapport à ladite particule d'intérêt (112A, 312A) de façon à ajuster la focalisation d'un photodétecteur (132 ; 332) situé dans le plan image dudit système optique . 3028616 39 16. Procédé d'analyse selon l'une quelconque des revendications 1 à 15, caractérisé en ce qu'il comprend en outre une étape de comptage du nombre de particules biologiques (112) présentes dans l'échantillon (111 ; 311). 5 17. Procédé d'analyse selon l'une quelconque des revendications 1 à 16, caractérisé en ce que les particules biologiques (112) consistent en des bactéries, des spores, des cellules, des levures ou des micro-organismes. 18. Dispositif (100 ; 300) d'analyse d'un échantillon (111 ; 311) recevant des 10 particules biologiques (112), parmi lesquelles une particule d'intérêt (112A; 312A), caractérisé en ce qu'il comprend : une première source lumineuse (120; 320) ; un ensemble d'imagerie (130; 330) comprenant un système optique (131; 331) et un capteur (132 ; 332), tels que le capteur se trouve dans un plan 15 image du système optique, ledit plan image étant le conjugué, par le système optique, d'un plan de mise au point (133, 333) ; un support (110) adapté à recevoir l'échantillon, disposé entre la première source lumineuse et l'ensemble d'imagerie ; des moyens de translation (140 ; 340), adaptés à déplacer le support (110) 20 relativement aux moyens d'imagerie (130 ; 330), selon un axe parallèle à l'axe optique (135 ; 335) du système optique ; et des moyens de calculs (150 ; 350) : reliés à une mémoire (151 ; 351) stockant une information relative à un point de repère (118 ; 318) situé sur une première 25 interface (116, 117; 316, 317) de l'échantillon ou à une distance connue de celle-ci selon l'axe optique du système optique ; recevant en entrée une image acquise par le capteur, dite image de référence lorsque ledit plan de mise au point est décalé par rapport audit point de repère ; 3028616 40 - adaptés à réaliser, à partir de l'image de référence, une construction numérique d'une série d'images reconstruites, associées chacune à un décalage prédéterminé (Di'; Di2; Din ; DL' ; Di_n) du plan de mise au point le long de l'axe optique du 5 système optique ; et fournissant en sortie la distance (Dref) selon un axe parallèle à l'axe optique (135 ; 335) du système optique, entre la particule d'intérêt (112A ; 312A) et le plan de mise au point (133 ; 333). 10