FR 3031180 A1 20160701 FR 1463445 A 20141230 1 MÉTHODE ET DISPOSITF DE CONTRÔLE DE LA POSITION D'UNE PARTICULE TRANSPARENTE PAR PILE D'IMAGES HOLOGRAPHIQUES DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne de manière générale le domaine de l'analyse optique de particules. Elle trouve notamment à s'appliquer au diagnostic microbiologique, plus particulièrement à l'analyse de microorganismes par spectroscopie (notamment par spectroscopie Raman ou par spectroscopie de fluorescence) ou encore par diffraction ou par diffusion. ÉTAT DE LA TECHNIQUE ANTÉRIEURE L'analyse de particules de faible taille, telles que des particules biologiques, et plus particulièrement des microorganismes, peut être classiquement réalisée au moyen d'une sonde optique. Cette sonde optique illumine la particule à analyser au moyen d'un faisceau focalisé ou d'un faisceau laser et mesure une caractéristique physique du faisceau après interaction avec la particule, par exemple son spectre de fluorescence ou son spectre Raman. L'analyse d'une particule peut permettre de l'identifier ou de la caractériser. Le positionnement du faisceau d'analyse et plus précisément de la tache de focalisation ou du spot laser sur la particule est généralement effectué de manière visuelle sous microscopie. Le positionnement du faisceau d'analyse par rapport à la particule à analyser doit être particulièrement précis, la taille de la tache de focalisation ou du spot laser étant généralement du même ordre que celle de la particule à analyser. En outre, ce positionnement est rendu délicat voire impossible lorsque la particule est transparente dans le spectre visible. Dans la suite, nous désignerons par particule transparente une particule qui présente une très faible absorption de la lumière visible (c'est-à-dire un coefficient de transmission supérieur ou égal à 70% dans le spectre visible ou plus généralement dans une bande spectrale de 300 nm à 1000 nm) et un faible écart d'indice de réfraction par rapport à son milieu environnant. Certaines cellules (eucaryotes et 3031180 2 procaryotes non marquées), certaines particules en matériaux diélectriques peuvent ainsi être qualifiées de transparentes. Dans le cas d'une particule transparente, il est connu d'observer la particule dans des conditions de défocalisation, autrement de décaler l'objectif du microscope le long de 5 son axe optique, par rapport à la particule, de manière à ce que la particule en question ne soit plus située dans le plan focal. Dans ce cas-là, si la sonde optique utilise le même objectif, le faisceau d'analyse n'est plus focalisé sur la particule à analyser. Une fois la particule centrée sur l'axe optique, il est certes possible de se placer à nouveau dans des conditions de focalisation, mais la particule peut s'être déplacée, avoir migré ou s'être 10 décrochée de son support (cas par exemple de bactéries capturées par un support), de sorte qu'il n'est pas possible de garantir l'intégrité de l'analyse et encore moins d'effectuer un asservissement dynamique du faisceau d'analyse sur la position de la particule. Le but de la présente invention est par conséquent de permettre d'effectuer un 15 contrôle en temps réel de la position d'une particule transparente, en particulier pour garantir la fiabilité de son analyse au moyen d'une sonde optique. EXPOSÉ DE L'INVENTION La présente invention est définie par une méthode de contrôle de la position d'une 20 particule d'intérêt au sein d'un échantillon, l'échantillon étant éclairé au moyen d'une première source lumineuse pour former une image holographique dans le plan image d'un premier système optique, ladite image étant acquise au moyen d'un capteur d'images, la particule d'intérêt étant transparente à la longueur d'onde de la première source lumineuse: 25 - en un premier instant, on acquiert une première image de l'échantillon, le capteur d'images étant décalé le long de son axe optique par rapport au plan image, on déduit de ladite première image, par propagation numérique, une première pile d'images holographiques dudit échantillon, correspondant à différentes valeurs de décalage du capteur par rapport au plan image, et l'on sélectionne dans cette première pile l'image 3031180 3 présentant le plus fort contraste comme première image caractéristique de l'échantillon ; - en un second instant, on acquiert une seconde image de l'échantillon, le capteur d'images étant décalé le long de son axe optique par rapport au plan image, on déduit de ladite seconde image, par propagation numérique, une seconde pile d'images 5 holographiques dudit échantillon correspondant à différentes valeurs de décalage du capteur par rapport au plan image, et l'on sélectionne dans ladite seconde pile l'image présentant le plus fort contraste comme seconde image caractéristique de l'échantillon ; - on estime le déplacement de la particule d'intérêt entre le premier et le second instant à partir des première et seconde images caractéristiques. 10 L'estimation du déplacement de la particule d'intérêt peut être obtenue par corrélation des première et seconde images caractéristiques. Alternativement, l'estimation du déplacement de la particule d'intérêt peut être obtenue en minimisant un écart quadratique entre les première et seconde images caractéristiques. 15 Alternativement encore, l'estimation du déplacement de la particule d'intérêt est obtenue par minimisation de l'information mutuelle entre les première et seconde images caractéristiques. Lorsque l'estimation du déplacement est obtenue à partir de la valeur maximale de la corrélation entre les première et seconde images caractéristiques, l'estimation du 20 déplacement peut être validée si la valeur maximale de la corrélation est supérieure à une valeur de seuil prédéterminée. Préalablement au premier instant, on peut placer la particule d'intérêt dans une position de consigne prédéterminée et, au premier instant, on vérifie, dans la première image caractéristique, la présence d'une particule à la position de consigne. 25 Avantageusement, après le second instant, on déplace l'échantillon de manière à corriger le déplacement estimé et l'on vérifie, dans la seconde image caractéristique, la présence d'une particule à la position de consigne. La présence d'une particule à la position de consigne peut être vérifiée en comparant l'intensité de l'image en ce point par rapport à une valeur de seuil ith = i' +36 où ici, est la valeur moyenne de l'intensité des pixels dans l'image et 6 est l'écart-type de 3031180 4 l'intensité des pixels dans l'image. La sélection de la première image caractéristique (respectivement la seconde image caractéristique) dans la première pile (respectivement dans la seconde pile) peut être réalisée en interpolant les valeurs d'écart-type de l'intensité des pixels dans les 5 différentes images de la première pile (respectivement de la seconde pile) pour obtenir une courbe d'écart-type en fonction du décalage par rapport au plan image et en recherchant une valeur de décalage correspondant au maximum de la courbe ainsi obtenue. On peut avantageusement extraire de la première pile d'images (respectivement 10 de la seconde pile d'images), une première (respectivement une seconde) pile de blocs d'image centrés sur la particule d'intérêt et l'on recherche dans la première (respectivement la seconde) pile de blocs celui présentant le contraste maximal, l'image de la première (respectivement de la seconde) pile d'images correspondant au bloc de contraste maximal étant sélectionnée comme première (respectivement seconde) image 15 caractéristique. La présente invention concerne également un dispositif de contrôle de la position d'une particule d'intérêt au sein d'un échantillon, ledit dispositif comprenant une première source lumineuse, un premier système optique adapté à former une image holographique de l'échantillon dans un plan image, un capteur d'images disposé à une distance 20 prédéterminée du plan image et adapté à acquérir ladite image holographique formée par le premier système optique, le capteur d'images acquérant une première image holographique en un premier instant et une seconde image holographique en un second instant, ledit dispositif comprenant en outre : - des moyens de calcul pour calculer par propagation numérique une première 25 (respectivement une seconde) pile d'images holographiques dudit échantillon, correspondant à différentes valeurs de décalage du capteur par rapport au plan image, à partir de la première (respectivement de la seconde) image holographique ; - des moyens pour sélectionner au sein de la première (respectivement de la seconde) pile d'images holographiques l'image présentant le plus fort contraste comme 30 première (respectivement seconde) image caractéristique de l'échantillon ; 3031180 5 - des moyens d'estimation du déplacement de la particule d'intérêt entre les premier et second instants à partir des première et seconde images caractéristiques. Ledit dispositif de contrôle peut comprendre en outre une platine motorisée sur laquelle est placée l'échantillon et des moyens de commande adaptés à recevoir le 5 déplacement estimé par les moyens d'estimation et à corriger ce déplacement de manière à maintenir ladite particule d'intérêt en une position de consigne prédéterminée. La présente invention concerne enfin un système d'analyse d'une particule d'intérêt au sein d'un échantillon comprenant un dispositif de contrôle de position tel que défini plus haut, une seconde source lumineuse adaptée à émettre un faisceau sur ladite 10 position de consigne de l'échantillon, un second système optique pour recevoir le faisceau après avoir interagi avec l'échantillon et des moyens de mesure pour mesurer une caractéristique physique du faisceau reçu par le second système optique. Selon un premier type de montage du système, les premier et second systèmes optiques partagent en commun un objectif de microscope. 15 Le contrôle de position par le dispositif de contrôle et la mesure de la caractéristique physique par les moyens de mesure peuvent alors être multiplexés dans le temps. Alternativement, la seconde source de lumière est une source laser de longueur d'onde distincte de la longueur d'onde de la première source de lumière. Dans ce cas, le 20 contrôle de position par le dispositif de contrôle et la mesure de la caractéristique physique par les moyens de mesure sont avantageusement simultanés. BRÈVE DESCRIPTION DES DESSINS 25 D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de modes réalisation préférentiels de réalisation en référence aux figures jointes parmi lesquelles : La Fig. 1 représente de manière schématique un premier exemple d'un système d'analyse mettant en oeuvre la méthode de contrôle de position selon l'invention ; 30 La Fig. 2 représente de manière schématique un second exemple d'un système 3031180 6 d'analyse mettant en oeuvre la méthode de contrôle de position selon l'invention; Les Figs. 3A et 3B représentent le résultat d'une détection de particules dans une image ; La Fig. 4 représente de manière schématique un organigramme de la méthode de 5 contrôle de la position d'une particule d'intérêt selon un mode de réalisation de l'invention ; La Fig. 5 représente de manière schématique l'obtention d'une pile d'images holographiques à partir d'une image défocalisée ; La Fig. 6 représente une courbe donnant l'écart-type dans une image de la pile en 10 fonction de son décalage par rapport à la position de focalisation ; La Fig. 7 représente de manière schématique l'extraction d'une pile de blocs à partir d'une pile d'images ; La Fig. 8 illustre de manière schématique une courbe d'écart-type dans un bloc d'image en fonction du décalage par rapport à la position de focalisation. 15 EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS La méthode de contrôle de la position d'une particule transparente est avantageusement mise en oeuvre dans un système d'analyse de particules au moyen d'une sonde optique. Ce système d'analyse sera illustré à l'aide de deux exemples de réalisation. 20 La Fig. 1 représente un premier exemple de système d'analyse optique mettant en oeuvre la méthode de contrôle de la position de particule transparente selon l'invention. Ce système d'analyse, 100, comprend une première voie dite voie d'analyse et une seconde voie dite voie de contrôle. La voie de contrôle permet de contrôler la position d'une particule d'intérêt et, le 25 cas échéant d'asservir sa position relative par rapport au faisceau d'analyse. La voie d'analyse permet quant à elle de mesurer une caractéristique physique du faisceau d'analyse après que celui-ci a interagi avec la particule d'intérêt de position contrôlée. La voie de contrôle comprend une première source lumineuse, de préférence 3031180 7 spatialement cohérente, et éventuellement temporellement cohérente, 110. De préférence, la première source est de faible largeur spectrale, par exemple inférieure à 200 nm voire à 100 nm ou même 25 nm. Cette source lumineuse peut être notamment une diode laser ou une diode électroluminescente. Le faisceau lumineux émis par la source 5 est amené au moyen d'une fibre optique sous l'échantillon à analyser, 120. Le cas échéant, le faisceau émis par la fibre est collimaté au moyen d'un diaphragme, 125. L'échantillon peut se présenter sous la forme d'un milieu liquide tel que de l'eau, une solution tampon, un milieu de culture ou un milieu réactif, dans lequel se trouvent des particules. Les particules sont avantageusement fixées sur un support tel qu'une lame de 10 microscope par exemple. Alternativement, l'échantillon peut se présenter sous la forme d'un milieu solide tel qu'une gélose dans lequel se trouvent les particules en question. Les particules peuvent se situer à l'intérieur du milieu ou à la surface de l'échantillon. Alternativement enfin, l'échantillon peut se présenter sous la forme d'un milieu 15 gazeux, les particules se situant dans ce milieu ou de préférence fixées sur un support tel qu'une lame de microscope par exemple. Les particules à analyser, P, peuvent être des particules biologiques, des poussières (de polluants par exemple) ou encore des particules de matériaux diélectriques. En tout état de cause, ces particules sont transparentes au sens défini plus haut. 20 Dans le cas d'un échantillon liquide, l'échantillon peut être enfermé dans une chambre d'analyse, 130, délimitée verticalement par une lame inférieure 131, par exemple une lame de microscope conventionnelle, et une lame supérieure 132. La chambre d'analyse est alors délimitée latéralement par un joint étanche 135. Les lames inférieure et supérieure sont transparentes à la longueur d'onde de la première source lumineuse. 25 La voie de contrôle comprend en outre un premier système optique comportant un objectif de microscope 140. Ce premier système optique est optionnellement équipé d'un filtre optique (non représenté) pouvant être situé devant ou derrière l'objectif. Après diffraction par l'échantillon, le faisceau de la première source traverse 30 l'objectif et est réfléchi par la lame réfléchissante 145 vers un capteur d'images 160, CCD 3031180 8 ou CMOS, par exemple le capteur d'une caméra. Le premier système optique comprend en outre une lentille de tube 150 permettant de focaliser le faisceau sur le capteur d'images. Le plan image en sortie de l'objectif est rejeté à l'infini et la lentille de tube 150 le ramène dans le plan du capteur 160. 5 II conviendra de noter que la lentille de tube peut être placée en amont de la lame réfléchissante ou bien en aval de cette dernière comme représenté en Fig. 1. Dans les conditions de focalisation, le plan objet no contenant les particules et le plan image H, sont conjugués à travers le premier système optique. Le capteur d'images acquiert alors une image en transmission d'une portion du plan objet. 10 La position du capteur d'images peut être réglée le long de son axe optique A. Si l'on note z la coordonnée le long de l'axe optique en prenant pour origine le point O où celui-ci coupe H,, on comprendra que z représente l'écart de défocalisation du capteur d'images. L'image formée sur le capteur d'images est une image holographique dans la 15 mesure où elle résulte de l'interférence entre une onde diffractée par les particules et une onde de référence ayant traversé l'échantillon sans avoir interagi avec lui. Alternativement, il est possible de diviser le faisceau lumineux en deux composantes, par exemple au moyen d'une seconde lame semi-transparente (non représentée). La première composante sert alors d'onde de référence et la seconde 20 composante est diffractée par l'échantillon, les deux composantes interférant dans le plan image du système optique. Les images de l'échantillon ainsi acquises font ensuite l'objet d'un traitement par la méthode de contrôle de position selon l'invention. Le dispositif 100 comprend en outre une voie d'analyse (non représentée sur cette 25 figure). Celle-ci comporte une seconde source de lumière, par exemple un laser, permettant d'illuminer une particule d'intérêt et d'effectuer la mesure d'une caractéristique du faisceau ayant interagi avec elle. La voie d'analyse peut partager avec la voie de contrôle tout ou partie du premier système optique. Dans ce cas le faisceau laser est choisi pour être 30 orienté dans l'axe du premier système optique, D . Après interaction avec la particule, le 3031180 9 faisceau laser traverse alors à nouveau le premier système optique et est analysé au moyen de moyens de mesure, par exemple un spectromètre ou un détecteur d'image. La voie d'analyse consiste, par exemple, en une voie de spectrométrie Raman ou une voie de mesure de fluorescence. 5 Dans cet exemple de réalisation, la voie de contrôle et la voie d'analyse sont multiplexées temporellement. Autrement dit, pendant un premier intervalle de temps, on effectue un contrôle de la position de la particule d'intérêt au moyen de la première source de lumière, la seconde source de lumière étant alors éteinte. Ce contrôle permet par exemple d'asservir la position de la particule d'intérêt sur l'axe optique D du premier 10 système optique. Dans un second intervalle de temps, la première source de lumière est éteinte et la seconde source de lumière est allumée. Le faisceau laser ayant interagi avec la particule est analysée par le spectromètre. L'homme du métier comprendra que les voies de contrôle et d'analyse pourront être inversées. Plus précisément, le capteur d'images pourrait se trouver dans l'axe 15 optique D et le spectromètre dans l'axe optique A. La Fig. 2 représente un second exemple de système d'analyse optique dans lequel la méthode de contrôle de la position de particule transparente selon l'invention peut être avantageusement mise en oeuvre. 20 Ce système d'analyse, 200, comprend comme précédemment une voie de contrôle et une voie d'analyse. Sur cette figure les éléments portant les mêmes signes de référence qu'en Fig. 1 sont identiques à ceux précédemment décrits. La voie de contrôle est identique à celle précédemment décrite, à ceci près qu'elle 25 comprend une lame séparatrice (par exemple une première lame dichroïque) 147 en lieu et place de la lame semi-réfléchissante. La voie d'analyse comprend une seconde source de lumière, par exemple un laser, 170, un second système optique (non représenté) pour recueillir le faisceau ayant interagi avec la particule d'intérêt et des moyens de mesure (180) d'une caractéristique 30 physique de ce faisceau. 3031180 10 Le faisceau émis par le laser 170 est tout d'abord réfléchi par une seconde lame séparatrice (par exemple une seconde lame dichroïque) 175, vers l'échantillon. Le faisceau laser réfléchi sur la lame séparatrice est aligné selon l'axe optique D. Après interaction avec la particule d'intérêt, P*, le faisceau laser est réfléchi, traverse à nouveau le premier 5 système optique, les première et seconde lames séparatrices, le second système optique, pour être analysé par les moyens de mesure 180. Dans le système d'analyse 200, la voie de contrôle et la voie d'analyse sont multiplexées en longueur d'onde. Plus précisément, les première et seconde sources de lumière ont des spectres disjoints. Par exemple, la première source de lumière peut être 10 une diode laser GaInN avec un spectre centré autour de 450 nm et la seconde source de lumière peut être un laser YAG doublé émettant à 532 nm. La première lame séparatrice réfléchit la lumière à 450 nm vers le capteur d'images et laisse passer celle à 532 nm. Si les moyens de mesure sont constitués par un spectromètre Raman, la seconde lame séparatrice est choisie pour réfléchir la lumière à 532 nm et laisser passer les longueurs 15 d'onde supérieures. Le multiplexage en longueur d'onde présente l'avantage de pouvoir contrôler la position de la particule d'intérêt en temps réel pendant l'analyse et par conséquent de pouvoir effectuer un asservissement en position avec une constante de temps très faible. L'asservissement en position dans les dispositifs 100 et 200 peut être réalisé par exemple 20 en montant la chambre d'analyse 130 sur une platine motorisée en X,Y,Z, celle-ci recevant alors des commandes de la voie de contrôle pour maintenir la particule en une position prédéterminée, de préférence sur l'axe optique D. Quel que soit le dispositif d'analyse, on supposera dans la suite que l'on a 25 préalablement localisé des particules au sein d'une image de l'échantillon acquise par le capteur d'images. La localisation des particules peut être réalisée de manière classique au moyen d'un filtre morphologique lorsque les particules sont simplement translucides. Lorsque les particules sont transparentes, la localisation peut être par exemple réalisée au moyen de la méthode de localisation décrite dans la demande FR1460947. 3031180 11 Après localisation des particules, on sélectionne parmi elles une particule d'intérêt, P*, que l'on souhaite analyser. On suppose dans la suite que la particule d'intérêt a été placée, le cas échéant à l'aide de la platine motorisée, à une position de consigne prise conventionnellement sur l'axe optique. Comme indiqué dans la partie introductive, la 5 position de la particule d'intérêt peut dériver dans le temps. La méthode de contrôle de position permet de déterminer si la position de la particule a évolué entre un premier instant et un second instant et si nécessaire, de commander la platine motorisée de manière à ramener la particule d'intérêt à ladite position de consigne. 10 La Fig. 4 représente de manière schématique un organigramme de la méthode de contrôle de la position d'une particule d'intérêt selon un premier mode de réalisation de l'invention. A l'étape 410, on acquiert en un premier instant une image initiale 10 de l'échantillon au moyen de la voie de contrôle. Plus précisément cette image est acquise par 15 le capteur d'images, celui-ci étant décalé d'un écart z0 par rapport à sa position de focalisation. L'image 1 ainsi obtenue est dite défocalisée. A l'étape 420, on obtient une pluralité d'images holographiques /1,...,4 de l'échantillon. Ces images holographiques /1,...,/, sont calculées à partir de l'image initiale 10, au 20 moyen d'un modèle de propagation numérique comme expliqué ci-après. Les images supplémentaires calculées sont celles qui seraient observées pour différentes positions zi,...,z, du capteur d'images. On trouvera une description de la méthode de calcul d'image par propagation numérique dans l'article de Sang-Hyuk Lee et al. intitulé « Holographic microscopy of 25 holgraphically trapped three-dimensional structures » publié dans Optics Express, Vol. 15 ; N° 4, 19 Février 2007, pp. 1505-1512. Plus précisément, si l'on note hz(r) la fonction de propagation de Rayleigh- Sommerfeld, soit : 3031180 laelkR hz(r)=-- 271- az R où z est l'écart de défocalisation, autrement dit l'écart par rapport au plan de mise au point, r = (x, y) est la position dans le plan de l'image, R2 = r2 + z2 et k = 2rn/ À est le 5 nombre d'onde relatif au milieu de propagation, l'onde dans le plan de coordonnée z peut être exprimée sous la forme : a(r , z) = a V zlexp (iv(r , z)) (2-1) 1 2 +00 (2-2) 10 a(r,z)=- 4;r f B(q)H_,(q)exp(iqr)d2q où B(q) est la transformée de Fourier de b(r), intensité de l'onde diffractée dans le plan focal (l'intensité de l'onde de référence est ici supposée constante), H z(q) est la transformée de Fourier de h z(r) et q est la variable duale de r dans la transformée de 15 Fourier. On comprend ainsi que l'on obtient une pile d'images holographiques /1,...,4, pour des écarts de défocalisation le long de l'axe optique, A, chaque image holographique in étant définie à partir de l'amplitude complexe a(r,zn). Plus précisément, chaque image holographique In est convertie en intensité. Selon une première variante, l'image 20 holographique I est définie par l'intensité a(r,zi2 . Selon une seconde variante préférée, l'image holographique In est définie par l'intensité S (a (r , zn))2 où S(a(r,zn)) est la partie imaginaire de a(r , z') . 12 (1) La Fig. 5 représente un exemple de pile d'images holographiques /1,...,4, obtenue 25 à partir d'une image défocalisée 10. 3031180 13 A l'étape 430, on sélectionne dans la pile d'images holographiques, /1,...,4, une image de référence /'f. Cette image de référence est celle qui approche au mieux aux conditions de focalisation. Autrement dit, cette image est celle qui est la plus proche du 5 plan image du premier système optique lorsque les particules sont dans le plan objet de cet ensemble optique. Ce critère de contraste maximal pourra être par exemple un écart-type maximal ou une valeur moyenne de gradient maximale dans cette zone. On a représenté en Fig. 6 une courbe donnant l'écart-type dans une image de la pile 10 en fonction de décalage z par rapport au plan image. En fait, les valeurs d'écart-type sont obtenues pour des valeurs discrètes de décalage zn et l'on effectue une interpolation entre ces valeurs d'écart-type, par exemple à l'aide d'une fonction spline, pour obtenir une courbe d'écart-type C(z) en fonction du décalage. L'image de référence /'f est celle de la pile d'images présentant un écart-type le plus élevé. Selon une variante, notamment 15 lorsque le pas d'échantillonnage du décalage axial est relativement grossier, on pourra déterminer la valeur zfix. - arg max(C(z)) et effectuer un nouveau calcul de propagation numérique pour obtenir une image holographique à la position de focalisation. On comprendra que cette variante requiert moins de calcul en raison du faible nombre d'images dans la pile. 20 On peut optionnellement vérifier que l'image de référence /nef contient bien une particule en son centre. Pour ce faire, on vérifie que l'intensité du pixel au centre de l'image de référence est par exemple supérieure à une valeur de seuil : 25 im - i' +3o- (3) où icn, et a sont respectivement la valeur moyenne et l'écart-type de l'intensité dans l'image /-'f . 3031180 14 Si l'intensité du pixel n'est pas supérieure à ladite valeur de seuil, un signal d'erreur est généré et l'opérateur est invité à sélectionner une nouvelle particule. A l'étape 440, on détecte, dans l'image de référence, une pluralité de particules dans le voisinage de la particule d'intérêt, chaque particule étant associée à un couple de 5 coordonnées dans l'image. Cette détection peut être facilitée par un traitement préalable de l'image, par exemple une détection de contours et des opérations de filtrage morphologiques (ouverture, fermeture). Grâce à ce filtrage préalable, les objets ne correspondant pas à des particules sont exclus de l'image. A l'issue de l'étape de détection, on obtient une première image caractéristique de l'échantillon, centrée sur l'axe optique, 10 comprenant la particule d'intérêt et une pluralité de particules voisines. Cette image est stockée dans une mémoire temporaire. On a représenté en Fig. 3A un exemple d'une telle image caractéristique. Les positions des différentes particules sont illustrées en Fig. 3B. On remarque en particulier la position de la particule d'intérêt P* , ici sur l'axe optique. 15 A un instant ultérieur, on répète les étapes 410 à 440, notées sur la figure 410' à 440', on obtient alors, à l'issue de l'étape de détection, une seconde image caractéristique de l'échantillon. On estime alors, en 450, le déplacement de la particule d'intérêt entre le premier 20 et le second instants, à partir des première et seconde images caractéristiques. Cette estimation peut être obtenue au moyen d'une technique de recalage d'image, par exemple en effectuant une corrélation des deux images caractéristiques, ou bien en minimisant un écart quadratique moyen entre ces deux images, ou encore en minimisant l'information mutuelle entre ces deux images, de manière connue en soi. 25 Dans le cas d'une estimation de déplacement par corrélation, on pourra générer une alerte si le maximum de la fonction de corrélation est inférieur à une valeur de seuil prédéterminée. En effet, cela peut signifier que la distribution des particules dans l'échantillon, au voisinage de la particule d'intérêt, a été alors fortement perturbée. A l'étape 460, la seconde image caractéristique est ensuite translatée au moyen 30 du vecteur de déplacement déterminé à l'étape précédente. 3031180 15 On peut alors vérifier ensuite à nouveau en 470, dans la seconde image caractéristique translatée, si une particule est bien présente sur l'axe optique. Si ce n'est pas le cas, un signal d'erreur est généré et l'opérateur est invité à sélectionner une nouvelle particule d'intérêt. L'absence de la particule d'intérêt peut en effet provenir du fait que 5 celle-ci s'est décrochée de son support. Si la particule d'intérêt est bien présente sur l'axe optique, on valide alors le vecteur de déplacement. Le vecteur de déplacement ainsi validé permet de contrôler la position de la particule d'intérêt et notamment d'asservir sa position sur l'axe optique en commandant la platine motorisée. 10 Selon une variante non représentée du premier mode de réalisation, après que l'image de référence, notée dans la suite iile.)f a été sélectionnée en 430 ou 430', à partir de la pile d'images holographiques, on détermine une région d'intérêt (ROI) sous la forme d'un bloc d'image, B.(1e) , centré sur la position de la particule d'intérêt dans l'image de référence 15 sélectionnée. On notera que ce bloc peut avoir une taille plus grande ou plus faible que celle de la particule à analyser. Ce bloc d'image, B(1) , est extrait de l'image de référence ref sélectionnée, /r(l). On extrait alors des autres images de la pile les blocs d'image correspondant à la même position que le bloc extrait de l'image de référence. On obtient ainsi une pile de 20 blocs d'image, Bi,...,B, , c'est-à-dire une restriction de la pile d'images à la zone d'intérêt centrée sur la particule à analyser. On comprend que chaque bloc de la pile correspond à une position axiale différente et donc à des écarts de défocalisation différents. La Fig. 7 représente de manière schématique une pile d'images /1,...,4, et une pile 25 de blocs d'image Bi,...,B, centrés sur une particule d'intérêt P . On affine la sélection de l'image de référence, et donc du bloc image de référence pour la particule à analyser. En effet, lorsque toutes les particules ne sont pas dans un même plan orthogonal à l'axe optique, le plan de mise au point pourra différer d'une particule à l'autre. On recherche alors parmi les blocs de la pile, et avantageusement parmi 3031180 16 les blocs des images de part et d'autre de l'image de référence /, sélectionnée à l'étape 430,430' celui qui correspond au mieux aux conditions de mise au point. Cette sélection pourra être effectuée, au moyen d'un critère de contraste maximal, portant cette fois-ci sur les blocs de la pile. Le critère de contraste maximal pourra être par exemple un 5 écart-type maximal ou une valeur moyenne de gradient maximale sur le bloc. Le bloc ainsi sélectionné devient alors le nouveau bloc de référence, Br(e2), de la pile, correspondant à une image i(,2f? de la pile d'images holographiques. La Fig. 8 donne la courbe donnant l'écart-type dans un bloc d'image de la pile en 10 fonction du rang de l'image de la pile. On comprend que la première sélection Br(elf) qui était optimale sur l'image /(1) ne l'est plus lorsque la mesure de contraste est limitée à la région ref d'intérêt, l'image 1 (e2) étant plus près des conditions de focalisation pour la particule d'intérêt. Les étapes 440, 440',450-470, sont alors effectuées sur la base de l'image de 15 référence /(2) en lieu et place de l'image 1(1) ref ref " Selon un second mode de réalisation de l'invention, le contrôle de la position de la particule d'intérêt est réalisé en estimant explicitement la position de la particule d'intérêt à chaque instant au moyen d'une pile d'images holographiques relative à cet 20 instant. Si l'on dispose de la position de consigne (ou de la position en un premier instant) de la particule d'intérêt, on peut alors déterminer si cette particule s'est déplacée et, le cas échéant estimer son déplacement. Plus précisément, le second mode de réalisation concerne une méthode de contrôle de la position d'une particule d'intérêt au sein d'un échantillon, l'échantillon étant 25 éclairé au moyen d'une première source lumineuse pour former une image holographique dans le plan image d'un premier système optique, ladite image étant acquise au moyen d'un capteur d'images, la particule d'intérêt étant transparente à la longueur d'onde de la première source lumineuse, ladite méthode de contrôle étant avantageuse en ce que : 3031180 17 - on dispose d'une information quant à une première position de ladite particule d'intérêt en un premier instant ; - on acquiert, à un second instant, postérieur au premier instant, une image de l'échantillon, le capteur d'images étant décalé le long de son axe optique par rapport au 5 plan image, on déduit de ladite image, par propagation numérique, une pile d'images holographiques dudit échantillon correspondant à différentes valeurs de décalage du capteur par rapport au plan image, et l'on sélectionne dans ladite pile l'image présentant le plus fort contraste comme image caractéristique de l'échantillon ; - on déduit la position de la particule audit second instant en fonction de l'image 10 caractéristique ainsi sélectionnée ; - on estime le déplacement de la particule d'intérêt entre ladite première position et ladite deuxième position. La variante décrite en relation avec la Fig. 7 s'applique également au second mode de réalisation. 15 Enfin, d'autres modes de réalisation pourront être envisagés par l'homme du métier sans sortir pour autant du cadre de la présente invention. REVENDICATIONS1. Méthode de contrôle de la position d'une particule d'intérêt au sein d'un échantillon, l'échantillon étant éclairé au moyen d'une première source lumineuse pour former une image holographique dans le plan image d'un premier système optique, ladite image étant acquise au moyen d'un capteur d'images, la particule d'intérêt étant transparente à la longueur d'onde de la première source lumineuse, caractérisée en ce que : - en un premier instant, on acquiert (410) une première image de l'échantillon, le capteur d'images étant décalé le long de son axe optique par rapport au plan image, on déduit (420) de ladite première image, par propagation numérique, une première pile d'images holographiques dudit échantillon, correspondant à différentes valeurs de décalage du capteur par rapport au plan image, et l'on sélectionne (430) dans cette première pile l'image présentant le plus fort contraste comme première image caractéristique de l'échantillon ; - en un second instant, on acquiert (410') une seconde image de l'échantillon, le capteur d'images étant décalé le long de son axe optique par rapport au plan image, on déduit (420') de ladite seconde image, par propagation numérique, une seconde pile d'images holographiques dudit échantillon correspondant à différentes valeurs de décalage du capteur par rapport au plan image, et l'on sélectionne (430') dans ladite seconde pile l'image présentant le plus fort contraste comme seconde image caractéristique de l'échantillon ; - on estime (450) le déplacement de la particule d'intérêt entre le premier et le second instant à partir des première et seconde images caractéristiques. 25 2. Méthode de contrôle de la position d'une particule d'intérêt selon la revendication 1, caractérisée en ce que l'estimation du déplacement de la particule d'intérêt est obtenue par corrélation des première et seconde images caractéristiques. 30 3. Méthode de contrôle de la position d'une particule d'intérêt selon la 3031180 19 revendication 1, caractérisée en ce que l'estimation du déplacement de la particule d'intérêt est obtenue en minimisant un écart quadratique entre les première et seconde images caractéristiques. 5 4. Méthode de contrôle de la position d'une particule d'intérêt selon la revendication 1, caractérisée en ce que l'estimation du déplacement de la particule d'intérêt est obtenue par minimisation de l'information mutuelle entre les première et seconde images caractéristiques. 10 5. Méthode de contrôle de la position d'une particule d'intérêt selon la revendication 2, caractérisée en ce que l'estimation du déplacement est obtenue à partir de la valeur maximale de la corrélation entre les première et seconde images caractéristiques, l'estimation du déplacement étant validée si la valeur maximale de la corrélation est supérieure à une valeur de seuil prédéterminée. 15 6. Méthode de contrôle de la position d'une particule d'intérêt selon l'une des revendications précédentes, caractérisée en ce que, préalablement au premier instant, on place la particule d'intérêt dans une position de consigne prédéterminée et, qu'au premier instant, on vérifie, dans la première image caractéristique, la présence d'une particule à la 20 position de consigne. 7. Méthode de contrôle de la position d'une particule d'intérêt selon la revendication 6, caractérisée en ce qu'après le second instant, on déplace l'échantillon de manière à corriger le déplacement estimé et l'on vérifie, dans la seconde image caractéristique, la présence d'une particule à la position de consigne. 8. Méthode de contrôle de la position d'une particule d'intérêt selon la revendication 6 ou 7, caractérisée en ce que la présence d'une particule à la position de consigne est vérifiée en comparant l'intensité de l'image en ce point par rapport à une valeur de seuil ith - i' +3o- où i' est la valeur moyenne de l'intensité des pixels dans 3031180 20 l'image et a est l'écart-type de l'intensité des pixels dans l'image. 9. Méthode de contrôle de la position d'une particule d'intérêt selon l'une des revendications précédentes caractérisée en ce que la sélection de la première image 5 caractéristique (respectivement la seconde image caractéristique) dans la première pile (respectivement dans la seconde pile) est réalisée en interpolant les valeurs d'écart-type de l'intensité des pixels dans les différentes images de la première pile (respectivement de la seconde pile) pour obtenir une courbe d'écart-type en fonction du décalage par rapport au plan image et en recherchant une valeur de décalage correspondant au maximum de la 10 courbe ainsi obtenue. 10. Méthode de contrôle de la position d'une particule selon l'une des revendications précédentes, caractérisée en ce que l'on extrait de la première pile d'images (respectivement de la seconde pile d'images), une première (respectivement une seconde) 15 pile de blocs d'image centrés sur la particule d'intérêt et que l'on recherche dans la première (respectivement la seconde) pile de blocs celui présentant le contraste maximal, l'image de la première (respectivement de la seconde) pile d'images correspondant au bloc de contraste maximal étant sélectionnée comme première (respectivement seconde) image caractéristique. 20 11. Dispositif de contrôle de la position d'une particule d'intérêt au sein d'un échantillon, caractérisé en ce qu'il comprend une première source lumineuse, un premier système optique adapté à former une image holographique de l'échantillon dans un plan image, un capteur d'images disposé à une distance prédéterminée du plan image et adapté 25 à acquérir ladite image holographique formée par le premier système optique, le capteur d'images acquérant une première image holographique en un premier instant et une seconde image holographique en un second instant, ledit dispositif comprenant en outre : - des moyens de calcul pour calculer par propagation numérique une première (respectivement une seconde) pile d'images holographiques dudit échantillon, 30 correspondant à différentes valeurs de décalage du capteur par rapport au plan image, à 3031180 21 partir de la première (respectivement de la seconde) image holographique ; - des moyens pour sélectionner au sein de la première (respectivement de la seconde) pile d'images holographiques l'image présentant le plus fort contraste comme première (respectivement seconde) image caractéristique de l'échantillon ; 5 - des moyens d'estimation du déplacement de la particule d'intérêt entre les premier et second instants à partir des première et seconde images caractéristiques. 12. Dispositif de contrôle de la position d'une particule d'intérêt au sein d'un échantillon selon la revendication 11, caractérisé en ce qu'il comprend en outre une platine 10 motorisée sur laquelle est placé l'échantillon et des moyens de commande adaptés à recevoir le déplacement estimé par les moyens d'estimation et à corriger ce déplacement de manière à maintenir ladite particule d'intérêt en une position de consigne prédéterminée. 15 13. Système d'analyse d'une particule d'intérêt au sein d'un échantillon comprenant un dispositif de contrôle de position selon la revendication 9, une seconde source lumineuse adaptée à émettre un faisceau sur ladite position de consigne de l'échantillon, un second système optique pour recevoir le faisceau après avoir interagi avec l'échantillon et des moyens de mesure pour mesurer une caractéristique physique du 20 faisceau reçu par le second système optique. 14. Système d'analyse selon la revendication 10, caractérisé en ce que les premier et second systèmes optiques partagent en commun un objectif de microscope. 25 15. Système d'analyse selon la revendication 14, caractérisé en ce que le contrôle de position par le dispositif de contrôle et la mesure de la caractéristique physique par les moyens de mesure sont multiplexés dans le temps. 16. Système d'analyse selon la revendication 14, caractérisé en ce que la seconde 3 0 source de lumière est une source laser de longueur d'onde distincte de la longueur d'onde 3031180 22 de la première source de lumière, et que le contrôle de position par le dispositif de contrôle et la mesure de la caractéristique physique par les moyens de mesure sont simultanés.