Method for counting and characterization of particles in a fluid in movement

Patent metadata

Field	Value
Publication number	US-10467764-B2
Application number	US-201715856679-A
Country	US
Kind code	B2
Filing date	Dec 28, 2017
Priority date	Dec 28, 2016
Publication date	Nov 5, 2019
Grant date	Nov 5, 2019

Abstract

A method allowing particles to be tracked in a moving fluid, via an optical method. The particles are in motion in a fluidic chamber. An image of the fluidic chamber is acquired, so as to obtain three-dimensional positions of particles in the fluidic chamber at a first time. Three-dimensional positions of 10 particles at a second time are also obtained, the second time being subsequent to the first time. On the basis of the obtained three-dimensional positions, potential movements of particles, between said times, are established. On the basis of a model of movement of the particles, potential movements are validated. The validated movements allow the particles in the fluid to be counted. In addition, if 15 the particles are of different nature, the movement model may comprise a component of movement of the particles with respect to the fluid that is characteristic of this difference. Determining this component then allows the particles to be characterized.

First claim

The invention claimed is:

1. A method for counting particles moving in a fluid flowing through a fluidic chamber, the method comprising:
a) placing the fluidic chamber between a light source and an image sensor, the image sensor lying in a detection plane;
b) illuminating the fluidic chamber with the light source, the light source emitting an incident light wave that propagates along a propagation axis, and acquiring, with the image sensor, a first image representative of an exposure wave to which the image sensor is exposed, the image sensor including various pixels, each pixel being associated with a radial coordinate in the detection plane;
c) on the basis of the acquired first image, obtaining three-dimensional coordinates of particles, in the fluidic chamber, at a first time;
d) obtaining three-dimensional coordinates of particles in the fluidic chamber at a second time, subsequent to the first time;
e) on the basis of the coordinates of the particles obtained at the first time and at the second time, determining potential movements of the particles between said times;
f) acquiring a model of the movement of the fluid in the fluidic chamber;
g) on the basis of the model of the movement of the fluid acquired in step f), validating movements among the potential movements calculated in step e); and
h) on the basis of the movements validated in step g), determining a number of particles and/or coordinates of the particles at the first time and/or at the second time,
wherein step c) comprises
ci) obtaining a first image of interest from the first image acquired in step b), and applying a digital propagation operator to the first image of interest in order to propagate the first image of interest, b a plurality of reconstruction distances, along the propagation axis, so as to obtain a first stack of images, including as many reconstructed complex images as there are reconstruction distances, each reconstructed complex image being representative of the exposure wave to which the image sensor is exposed;
cii) for at least one radial coordinate defined in the first image of interest, determining a reconstruction distance that maximizes variation in a component of each reconstructed complex image forming the first stack of images, along an axis parallel to the propagation axis and passing through said radial coordinate, the determined reconstruction distance forming a transverse coordinate associated with said radial coordinate, a value of the component calculated at the reconstruction distance being a maximum value associated with said radial coordinate, substep cii) being carried out for all or some radial coordinates associated with pixels of the first image of interest;
ciii) establishing a list of three-dimensional positions, each three-dimensional position including a radial coordinate and the associated transverse coordinate determined in substep cii), each three-dimensional position being associated with the maximum value obtained in substep cii); and
civ) selecting three-dimensional positions depending on the associated maximum values.

2. The method of claim 1 , wherein step c) comprises:
on the basis of each reconstructed complex image, obtaining radial coordinates of particles in the fluidic chamber at the first time.

3. The method of claim 1 , wherein the first image of interest is:
the first image acquired in step b);
or the first image acquired in step b), from which an image of the fluidic chamber is subtracted, the image of the fluidic chamber is acquired by the image sensor, prior or subsequently to the acquisition of the first image, the subtraction being weighted by a weighting term;
or the first image acquired in step b), from which an average of images acquired prior and subsequently to the acquisition of the first image is subtracted.

4. The method of claim 1 , wherein, in substep cii), the component includes a real part, or an imaginary part, or a modulus, or a phase of each reconstructed complex image forming the first stack of images.

5. The method of claim 1 , wherein substep civ) comprises:
forming a first maxima image, each pixel of the first maxima image is associated with a three-dimensional position determined in substep ciii) and is assigned the maximum value determined, in substep ciii), for said three-dimensional position;
selecting, in the first maxima image, pixels having values that are maximum in a neighbouring zone defined around each pixel; and
calculating, for each selected pixel, a signal-to-noise ratio depending on the maximum value and the values of pixels of the first maxima image located in a calculation zone lying around said selected pixel;
such that each three-dimensional position is selected depending on the signal-to-noise ratio calculated for the selected pixels of the first maxima image.

6. The method of claim 1 , wherein step d) includes acquiring, with the image sensor, a second image, each pixel of the second image is associated with a radial coordinate in the detection plane.

7. The method of claim 6 , wherein step d) comprises:
di) obtaining a second image of interest from the acquired second image and applying a digital propagation operator to the second image of interest in order to propagate the second image of interest, by a plurality of reconstruction distances, along the propagation axis, so as to obtain a second stack of images, including as many reconstructed complex images as there are reconstruction distances, each reconstructed complex image being representative of an exposure wave to which the image sensor is exposed at the second time;
dii) for at least one radial coordinate defined in the second image of interest, determining a reconstruction distance that maximizes variation in a component of each reconstructed complex image forming the second stack of images; along an axis parallel to the propagation axis and passing through said radial coordinate, the determined reconstruction distance forming a transverse coordinate associated with said radial coordinate, a value of the component calculated at said reconstruction distance being a maximum value associated with the radial coordinate, substep dii) being carried out for all or some radial coordinates associated with pixels of the second image of interest;
diii) establishing a list of three-dimensional positions, each three-dimensional position including a radial coordinate and the associated transverse coordinate determined in substep dii), each three-dimensional position being associated with the maximum value obtained in substep dii); and
div) selecting three-dimensional positions depending on the associated maximum values.

8. The method of claim 7 , wherein, in substep di), the second image of interest is:
the acquired second image;
or the acquired second image, from which an image of the fluidic chamber is subtracted, the image of the fluidic chamber acquired by the image sensor, prior or subsequently to the acquisition of the second image, the subtraction being weighted by a weighting term between 0 and 1;
or the acquired second image; from which an average of images acquired prior and subsequently to the acquisition of the second image is subtracted.

9. The method of claim 4 ; wherein, in substep dii), the component includes a real part, or an imaginary part, or a modulus, or a phase of each reconstructed complex image forming the second stack of images.

10. The method of claim 7 , wherein substep div) includes:
forming a second maxima image, each pixel of the second maxima image is associated with a three-dimensional position determined in substep diii) and is assigned the maximum value determined, in substep diii), for the three-dimensional position;
selecting

Assignees

• Commissariat Energie Atomique

Inventors

• Joly Pierre
• Rousier Rodrigue
• Elvira David

Classifications

• H04N23/56

  provided with illuminating means · CPC title

• G01N15/1429

  Signal processing · CPC title

• G01N2015/1486

  Counting the particles · CPC title

• G06T2207/20224

  Image subtraction · CPC title

• G03H2001/0447

  In-line recording arrangement · CPC title