Holographic method for characterizing a particle in a sample

Patent metadata

Field	Value
Publication number	US-10845286-B2
Application number	US-201716086922-A
Country	US
Kind code	B2
Filing date	Mar 22, 2017
Priority date	Mar 23, 2016
Publication date	Nov 24, 2020
Grant date	Nov 24, 2020

Abstract

A method for holographic characterization of a particle contained in a sample, based on an image, or hologram, of the sample obtained by an image sensor when the sample is illuminated by a light source. The hologram is the subject of a holographic reconstruction, to obtain a reference complex image, representative of the light wave transmitted by the sample in a reconstruction plane. A holographic propagation operator is applied to the reference complex image, to obtain a plurality of secondary complex images, from which a profile is determined describing the change in an optical feature of the light wave transmuted by the sample along the axis of propagation of the light wave.

First claim

The invention claimed is:

1. A method for characterizing a particle contained in a sample, the method comprising:
a) illuminating the sample using a light source, the light source emitting an incident light wave propagating towards the sample along a propagation axis;
b) acquiring, using an image sensor, an image of the sample, formed in a detection plane, the sample being positioned between the light source and the image sensor, the image being representative of a light wave transmitted by the sample under effect of the illumination;
c) applying a propagation operator to the image acquired in b), to calculate a reference complex image, which is a complex image representative of the sample in a reconstruction plane;
d) selecting a radial position of the particle in a plane parallel to the detection plane;
e) from the reference complex image calculated in c), determining at least one characteristic quantity of the light wave transmitted by the sample, at a plurality of distances from the detection plane or from the reconstruction plane, wherein e) further includes :
ei) applying a propagation operator to the reference complex image, to calculate secondary complex images for a plurality of distances from the reconstruction plane or from the detection plane; and
eii) determining a characteristic quantity of the light wave transmitted by the sample at each of the plurality of distances, from the secondary complex images;

f) forming a profile representing a variation of the characteristic quantity determined in e) along an axis parallel to the propagation axis, the axis passing through the radial position selected in d); and
g) characterizing the particle using the profile formed in f),
wherein g) includes at least one of
classifying the particle according to a plurality of predetermined particle classes,
determining of a state of the particle, according to a plurality of predetermined states,
estimating a geometric parameter of the particle, and
estimating an optical parameter of the particle.

2. The method of claim 1 , wherein at least one characteristic quantity is determined from a modulus or an argument of a secondary complex image.

3. The method of claim 1 , wherein the characterization is performed by comparing the profile formed in f) with profiles determined in a learning phase.

4. The method of claim 1 , wherein, in d), the radial position of the particle is selected using the image acquired in b) or using the reference complex image calculated in c).

5. The method of claim 1 , wherein no magnifying optics are placed between the sample and the image sensor.

6. The method of claim 1 , wherein:
d) further includes selecting a plurality of radial coordinates that are representative of a same type of particle;
f) further includes forming a plurality of elementary profiles, at each selected radial coordinate, and combining the elementary profiles so as to form a profile representative of the particle.

7. The method of claim 1 , wherein, in c), the calculation of the reference complex image comprises:
ci) defining an initial image of the sample in the detection plane, from the image acquired by the image sensor;
cii) determining a complex image of the sample in a reconstruction plane by applying a propagation operator to the initial image of the sample defined in ci) or to the image of the sample in the detection plane resulting from the preceding iteration ;
ciii) calculating a noise indicator from the complex image determined in cii), the noise indicator depending on a reconstruction noise affecting the complex image;
civ) updating the image of the sample in the detection plane by adjusting phase values of the pixels of the image, the adjustment being carried out according to a variation of the indicator calculated in ciii) with the phase values;
cv) reiterating cii) to civ) until a convergence criterion is reached, so as to obtain a complex reference image of the sample in the detection plane, or in the reconstruction plane.

8. The method of claim 7 , wherein ciii) further includes:
for various pixels, calculating a quantity associated with each pixel, according to the value of the complex image determined in cii) for the pixel, or of a dimensional derivative of the complex image image for the pixel;
combining the quantities associated with various pixels, to obtain the noise indicator.

9. The method of claim 8 , wherein the noise indicator is a norm of order lower than or equal to 1, calculated from the quantities associated with each pixel.

10. The method of claim 8 , wherein, in ciii), the quantity associated with each pixel is calculated from the modulus of a dimensional derivative, for the pixel, of the complex image determined in cii).

11. The method of claim 8 , wherein:
in ci), the initial image of the sample is defined by a normalization of the image acquired by the image sensor by an image representative of the light wave emitted by the light source;
in ciii), the quantity associated with each pixel is calculated according to the value of the complex image determined in cii), for the pixel, subtracted from a strictly positive number, in particular the number 1.

12. The method of claim 7 , wherein, in civ), the value of the phase of each pixel is adjusted by forming a vector, referred to as the phase vector, each term of which corresponds to the value of the phase of a pixel of the image of the sample in the detection plane, this vector being updated, in each iteration, so as to minimize or to maximize the noise indicator calculated in ciii), on the basis of a gradient of the noise indicator according to each term of the phase vector.

13. The method of claim 1 , wherein the particle is chosen from: a cell or a microorganism or a microbead or an exosome or a droplet of an emulsion or a cell nucleus or a piece of cellular debris.

14. A device for characterizing a particle contained in a sample, the device including:
a light source configured to emit an incident light wave propagating towards the sample along a propagation axis;
an image sensor configured to acquire an image of the sample, formed in a detection plane, the sample being positioned between the light source and the image sensor, the image being representative of a light wave transmitted by the sample under effect of the illumination;
a holder, configured to hold the sample between the light source and the image sensor;
a processor, configured to receive the image of the sample acquired by the image sensor and
apply a propagation operator to the image acquired, to calculate a reference complex image, which is a complex image representative of the sample in a reconstruction plane;
select a radial position of the particle in a plane parallel to the detection plane:
from the reference complex image calculated, determine at least one characteristic quantity of the light wave transmitted by the sample, at a plurality of distances from the detection plane or from the reconstruction plane, by:
applying a propagation operator to the reference complex image, to calculate secondary complex images for a plurality of distances from the reconstruction plane or from the detection plane; and
determining a characteristic quantity of the light wave transmitted by the sample at each of the plurality of distances, from the secondary complex images;

forming a profile representing a variation of the characteristic quantity determined along an axis parallel to the propagation axis, the axis passing through the radial position selected; and
characterizing the particle using the profile formed,
wherein the characterizing includes includes at l

Assignees

• Commissariat Energie Atomique
• Horiba Abx Sas
• Iprasense Sas

Inventors

• Allier Cedric
• Ali Cherif Anais
• Blandin Pierre
• Esteban Geoffrey
• Herve Lionel
• Isebe Damien

Classifications

• G03H2001/0883

  Reconstruction aspect, e.g. numerical focusing · CPC title

• G03H2001/0816

  Iterative algorithms · CPC title

• G03H2001/0447

  In-line recording arrangement · CPC title

• G03H2001/0033

  in hologrammetry for measuring or analysing · CPC title

• G03H1/0866 Primary

  Digital holographic imaging, i.e. synthesizing holobjects from holograms · CPC title