EP3218769A1 - Procede d'analyse comprenant la determination holographique d'une position d'une particule biologique - Google Patents

Procede d'analyse comprenant la determination holographique d'une position d'une particule biologique

Info

Publication number
EP3218769A1
EP3218769A1 EP15801133.8A EP15801133A EP3218769A1 EP 3218769 A1 EP3218769 A1 EP 3218769A1 EP 15801133 A EP15801133 A EP 15801133A EP 3218769 A1 EP3218769 A1 EP 3218769A1
Authority
EP
European Patent Office
Prior art keywords
particle
interest
optical system
image
distance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
EP15801133.8A
Other languages
German (de)
English (en)
Other versions
EP3218769B1 (fr
Inventor
François PERRAUT
Pierre Joly
Quentin JOSSO
Meike KLOSTER-LANDSBERG
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Biomerieux SA
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Original Assignee
Biomerieux SA
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Biomerieux SA, Commissariat a lEnergie Atomique et aux Energies Alternatives CEA filed Critical Biomerieux SA
Publication of EP3218769A1 publication Critical patent/EP3218769A1/fr
Application granted granted Critical
Publication of EP3218769B1 publication Critical patent/EP3218769B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/50—Depth or shape recovery
    • G06T7/55—Depth or shape recovery from multiple images
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1434—Optical arrangements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
    • G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
    • G01N21/453—Holographic interferometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G01N21/4788—Diffraction
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/0005—Adaptation of holography to specific applications
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0402—Recording geometries or arrangements
    • G03H1/041—Optical element in the object space affecting the object beam, not otherwise provided for
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • G06T15/10—Geometric effects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/97—Determining parameters from multiple pictures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/01—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials specially adapted for biological cells, e.g. blood cells
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1429—Signal processing
    • G01N15/1433—Signal processing using image recognition
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N2015/1006—Investigating individual particles for cytology
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1434—Optical arrangements
    • G01N2015/144—Imaging characterised by its optical setup
    • G01N2015/1445—Three-dimensional imaging, imaging in different image planes, e.g. under different angles or at different depths, e.g. by a relative motion of sample and detector, for instance by tomography
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1434—Optical arrangements
    • G01N2015/1452—Adjustment of focus; Alignment
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1434—Optical arrangements
    • G01N2015/1454—Optical arrangements using phase shift or interference, e.g. for improving contrast
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/0005—Adaptation of holography to specific applications
    • G03H2001/0033—Adaptation of holography to specific applications in hologrammetry for measuring or analysing
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/0005—Adaptation of holography to specific applications
    • G03H2001/005—Adaptation of holography to specific applications in microscopy, e.g. digital holographic microscope [DHM]
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0402—Recording geometries or arrangements
    • G03H2001/0428—Image holography, i.e. an image of the object or holobject is recorded
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
    • G03H2001/0447—In-line recording arrangement
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
    • G03H2001/0452—Digital holography, i.e. recording holograms with digital recording means arranged to record an image of the object
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0465—Particular recording light; Beam shape or geometry
    • G03H2001/0467—Gated recording using pulsed or low coherence light source, e.g. light in flight, first arriving light
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • G03H2001/0883—Reconstruction aspect, e.g. numerical focusing

Definitions

  • the present invention relates to the field of the analysis of a sample receiving biological particles, this analysis consisting in particular of determining a position of at least one of said biological particles, along an axis of the depth of the sample.
  • a commonly used solution is to acquire a series of images of the sample at different depths in the sample. We then look for the image on which we observe a clear image of the biological particle.
  • a disadvantage of this solution is that it does not allow to know the position of a transparent biological particle whose clear image merges with the image of the surrounding medium.
  • An object of the present invention is to provide a method and a device not having this disadvantage.
  • an object of the present invention is to provide a method and a device for determining the position of any biological particle, along at least one axis of the three-dimensional space.
  • This objective is achieved with a method of analyzing a sample receiving biological particles, among which a particle of interest, the method being implemented in a device comprising the sample disposed between a first light source and an optical system said method comprising the steps of:
  • the invention therefore advantageously comprises a step of detecting the reference point, by placing the optical system so that its focus plane is placed on the reference point, on said first interface or at the known distance thereof, and arrangement of the optical system such that the particle of interest is out of the focus plane of the optical system.
  • the reference image is acquired for this arrangement of the optical system, such that the particle of interest is outside the focus plane of the optical system.
  • This arrangement of the optical system can correspond directly to the arrangement of the optical system implemented to detect the reference point.
  • the reference image, or defocused image is acquired directly after detection of the reference point.
  • the method comprises a step of shifting the focusing plane of the optical system relative to said cue point, by a translation of a support receiving the sample relative to the optical system.
  • the reference image, or defocused image is acquired after this offset step.
  • the invention makes it possible to precisely know a distance between the plane of focus and said reference point. It is thus possible to precisely know a distance between the particle of interest and the focusing plane, which then makes it possible to better position an analysis device to study the particle of interest.
  • the method comprises a determination of the distance between the particle of interest and said reference point, from said distance between the particle of interest and said focusing plane and said useful distance.
  • It may further comprise a determination of the presence of said particle of interest in the medium from the distance between said particle of interest and said focusing plane.
  • the determination of the useful distance comprises the following sub-steps: illumination of the reference point by a laser beam, said laser beam being focused by said optical system; acquiring an image on the sensor representing the reflection of the laser beam on an interface of the sample;
  • the biological particles adhere to a second interface of the sample, the first and second interfaces being distinct or merged.
  • the distance along an axis parallel to the optical axis of the optical system, between the position of the focusing plane associated with the reference image, and the projection of the reference point along this axis and on the second interface, is included for example between +5 ⁇ and +2000 ⁇ or between -5 ⁇ and -2000 ⁇ .
  • the illuminated region may comprise a plurality of biological particles.
  • the first light source typically has a spectral width of less than 200 nm.
  • the method further comprises determining, from said series of reconstructed images, the position of the particle of interest in a plane orthogonal to the optical axis of the optical system.
  • Each reconstructed image is formed by a real part and an imaginary part, and among the real and imaginary parts, only the imaginary parts of the reconstructed images are advantageously used to determine the distance along an axis parallel to the optical axis of the optical system, between the particle of interest and the landmark.
  • Each reconstructed image may be associated with an offset along the optical axis of the optical system and with a value of a useful parameter, so as to constitute a function describing the evolution of the useful parameter according to said offset, the determination of the distance between the particle and the plan of point implementing a search of a remarkable value, in particular an extremum, a point of inflection or a zero crossing of said function.
  • the optical system is advantageously displaced with respect to said particle of interest so as to focus the analysis laser beam (363) on the particle of interest. 'interest.
  • the optical system is advantageously displaced with respect to said particle of interest so as to adjust the focus of a photodetector located in the image plane.
  • the method may further comprise a step of counting the number of biological particles present in the sample.
  • biological particles may comprise bacteria, spores, cells, yeasts or microorganisms.
  • the invention also relates to a device for analyzing a sample receiving biological particles, among which a particle of interest, the device comprising:
  • an imaging assembly comprising an optical system and a sensor, such that the sensor is in an image plane of the optical system, said image plane being the conjugate, by the optical system, of a focusing plane;
  • a support adapted to receive the sample, arranged between the first light source and the imaging assembly;
  • translational means adapted to move the support relative to the imaging means, along an axis parallel to the optical axis of the optical system
  • the device may also include a laser, adapted to provide a laser beam aligned with the optical axis of the optical system, and coupled to the optical system so that the focusing point of the laser beam corresponds to the focusing plane.
  • FIG. 1 illustrates a first embodiment of an analysis device according to the invention
  • FIG. 2 schematically illustrates a first embodiment of an analysis method according to the invention
  • FIG. 3 illustrates a second embodiment of an analysis device according to the invention
  • FIG. 4 schematically illustrates a detail of a second embodiment of an analysis method according to the invention
  • FIGS. 5A to 5C illustrate a detail of the method illustrated in FIG. 4;
  • FIG. 6 illustrates a reference image according to the invention;
  • FIG. 7 schematically illustrates a detail of a third embodiment of an analysis method according to the invention.
  • FIG. 8 illustrates a profile used in a fourth embodiment of the method according to the invention. DETAILED PRESENTATION OF PARTICULAR EMBODIMENTS
  • analysis method and the analysis device according to the invention will be described together, the analysis device according to the invention being adapted to the implementation of the analysis method according to the invention.
  • FIG. 1 illustrates a first embodiment of an analysis device 100 according to the invention.
  • FIG. 2 illustrates a first embodiment of an analysis method according to the invention, implemented in the analysis device 100.
  • the analysis device 100 comprises a support 110, adapted to receive a sample 111.
  • the support 110 is for example a clamp, or a transparent plate, or a plate pierced with an opening for passing light rays.
  • the sample 111 is a transparent or translucent medium, that is to say having a transmission coefficient greater than or equal to 70% in the visible spectrum or more generally in a spectrum between 300 nm and 1000 nm.
  • Sample 111 is delimited by interfaces, or boundaries, or boundaries between the sample and a medium in direct contact with it.
  • the sample consists of a liquid such as water, a buffer solution, a liquid comprising a reagent, a culture medium, and biological particles 112 located in this liquid.
  • the sample consists of a solid medium such as agar, and biological particles located on this solid medium.
  • the sample consists of a gas in which there are biological particles.
  • the biological particles may be located within the sample (for example bathing in a liquid), or flush on the surface of the sample (for example being located on an agar).
  • the biological particles 112 denote for example bacteria, spores, cells, yeasts, or any type of microorganism.
  • One of these biological particles is named particle of interest 112A.
  • the sample 111 comprises a culture medium enclosed in a fluidic chamber delimited by a lower blade 113, for example a standard microscope slide, an upper blade 114, and an adhesive 115 surrounding the liner laterally. Sample 111 and connecting together the upper and lower blades.
  • the upper and lower blades 113, 114 are transparent in the visible and, where appropriate, other useful wavelengths.
  • the boundary between the sample and the lower blade 113 defines a lower interface 116 of the sample 111.
  • the boundary between the upper blade 114 and the sample defines an upper interface 117 of the sample 111.
  • a first light source 120 is located upstream of the sample 111, in the direction of propagation of the light, from the first light source 111.
  • the terms upstream and downstream refer to the direction of propagation of the light, from the first light source 120 to the sample 111.
  • the first light source 120 is for example a laser, a light-emitting diode, a white lamp including a mercury vapor lamp filtered or not.
  • the first light source may include an optical fiber for bringing light under the lower blade 113.
  • the first light source 120 is advantageously temporally coherent. It has a spectral width advantageously less than 200 nm, or even less than 100 nm or even 25 nm.
  • the first light source 120 is spatially coherent.
  • the first light source illuminates in transmission a region 119 of the sample.
  • the particle of interest 112A is in region 119.
  • An imaging assembly 130 is located downstream of the upper blade 114, in the propagation direction of the light, from the first light source to the sample 111. In the example described, the imaging assembly 130 is above sample 111.
  • the imaging assembly 130 includes an optical system 131 and a sensor
  • the optical system 131 consists for example of an objective, in particular a microscope objective. It has an object plane 133, also called a focus plane, and an image plane 134.
  • the image plane 134 is the conjugate (or image) of the focusing plane by the optical system. In other words, an object in the focus plane is a sharp image in the image plane.
  • the sensor 132 is for example a matrix sensor of the CCD or CMOS type. I l is located in the image plane 134. Thus, the sensor 132 acquires a transmission image of a part of the focusing plane 133.
  • the focus plane is the object focal plane of the optical system, and the image plane is returned to infinity.
  • a proximity optics for example a tube lens, makes it possible to focus the image plane on the sensor 132.
  • sensor 132 groups together the sensor 132 and its proximity optics.
  • the optical system 131 has an optical axis 135.
  • the optical axis is orthogonal to the plane of focus and to the image plane.
  • the optical axis 135 defines the depth of the sample 111. This axis connects the lower and upper interfaces 116, 117 of the sample 111.
  • the optical axis 135 is substantially orthogonal to the lower and upper interfaces 116, 117
  • the term interface refers to a boundary of the sample 111 studied.
  • the analysis within the meaning of the invention comprises in particular a determination of the position of the particle of interest 112A along an axis parallel to the optical axis 135.
  • the device 100 further comprises translation means 140, adapted to move the imaging assembly 130 along an axis parallel to the optical axis 135.
  • the translation means 140 may move only the optical system 131, the sensor 132 remaining fixed.
  • the translation means 140 are adapted to move the support 110 along an axis parallel to the optical axis 135.
  • the translation means 140 furthermore perform translations along two other axes defining a plane orthogonal to the optical axis 135.
  • the sensor 132 is connected to computing means 150, in particular a processor or a microprocessor.
  • the device 100 is adapted to the implementation of the method according to the invention, illustrated schematically in FIG.
  • the method includes a first step 21, determining a landmark 118.
  • This step consists of choosing, or defining the landmark 118. This choice is generally arbitrary.
  • the reference point is located on a first interface of the sample 111, or at a known distance from this first interface, this distance being defined along the optical axis 135.
  • the cue point is located on the upper interface 117 or lower 116.
  • the cue point is located on the upper interface 117 of the sample.
  • the reference point is not located directly on an interface of the sample, but at a known distance from it.
  • the landmark is on one side of the upper blade 114 opposite the sample, the thickness of the upper blade 114 being known.
  • the calculation means 150 are connected to a memory 151 storing information relating to the reference point, in particular the distance along an axis parallel to the optical axis 135, between this reference point and said first interface and possibly its coordinates in a plane orthogonal to the optical axis 135.
  • the sample 111 is then illuminated with the first light source 120, thereby forming an illuminated region 119 in the sample as described above (step 22).
  • a transmission image of the illuminated region 119 is acquired using the imaging assembly 130 (step 23).
  • This image is named reference image.
  • the focusing plane 133 is located at a known distance D1 from the point of interest 118, this distance being defined along an axis parallel to the optical axis 135.
  • the particle of interest 112A is outside plan 133.
  • the focusing plane is outside the sample 111.
  • the focusing plane is located between 5 ⁇ and 1500 ⁇ above or below (upstream or downstream) of the reference point, preferably between 5 ⁇ and 1000 ⁇ , or even between 5 ⁇ and 800 ⁇ .
  • the distance is defined along an axis parallel to the optical axis 135.
  • the focusing plane when acquiring the reference image, is located at a distance greater than 2 ⁇ , and preferably greater than 5 ⁇ , preferably in the range [5 ⁇ - 1500 ⁇ ] or [5 ⁇ - 1000 ⁇ ], or [5 ⁇ - 800 ⁇ ] of an interface delimiting the sample 111.
  • it has a priori information qua nt the position of a particle. The focus plane is then shifted, relative to this a priori, a value as defined above.
  • an example of a priori information is the observation of biological particles adhering to a surface, for example the upper blade 114 delimiting the sample 111.
  • the knowledge of the position of the blade allows to establish a priori on the position of the particles, the latter being adherent to the face 117 of this adjacent blade of the sample.
  • the reference image is said to be the image of a portion of the illuminated region, since it is formed by light rays from this illuminated region.
  • the reference image is in particular a hologram, formed by the interference of light rays originating from the first light source and diffracted by a biological particle of the sample, with light rays coming from the first light source 120 and having passed through. the sample without being diffracted.
  • the reference image includes in particular a hologram associated with the particle of interest.
  • the calculation means 150 are connected to the sensor 132 to receive the reference image.
  • a digital construction of a series of reconstructed images is implemented.
  • Each reconstructed image corresponds to a predetermined offset, along the optical axis 135, of the focusing plane relative to the position of the focus plane associated with the reference image.
  • each reconstructed image is computed into a reconstruction plane, associated with a predetermined offset from the focus plane. More precisely, each reconstruction plane corresponds to the image, by the optical system, of an object plane offset from the focusing plane of this predetermined offset.
  • the offsets extend between two extreme positions vis-à-vis the plane of focus, these two extreme positions flanking the upper interface 117 and, more generally, the point of interest.
  • These offsets can be distributed according to a regular step, in particular a step between 0.10 ⁇ and 1 ⁇ , for example 0.25 ⁇ or 0.20 ⁇ .
  • the images acquired by the sensor 132 would be digitally reconstructed if the focusing plane was successively at the distances Di 1 , Di 2 ,... Di n , DU, ... Di n of the position of the plane in focus associated with the reference image.
  • the reference image is an experimentally acquired image, while the reconstructed images form digital propagation images.
  • the difference between the two extreme positions associated with these offsets depends on the sample 111 observed, and in particular on its thickness.
  • the two extreme positions are determined so that they are located on either side of the sample 111. from a priori on the position of a particle of interest in the sample, the two extreme positions are established to be arranged on either side of this position.
  • the digital reconstruction uses the Fourier transform of the reference image (possibly having undergone a preliminary treatment), to which a propagation operator is applied before returning it to real space.
  • the propagation operator is a function of the offset associated with the calculated reconstructed image.
  • the propagation operator is for example an integral based on the Rayleigh Sommerfeld equation.
  • Wilson et al. describes an example of a digital propagation implemented from the image of a microsphere. Another example of digital propagation is described by Lee et al. in the article “Microscopy Holography of Holographically Trapped Three-dimensional Structures,” 19 February 2007 / Vol. 15, No. 4 / OPTICS EXPRESS 1505-1512.
  • the reconstructed images are complex images, that is to say with a real part and an imaginary part. They comprise a set of points, each point of the image being assigned a complex magnitude.
  • the computing means 150 of the device 100 according to the invention are adapted to produce, from the reference image, a digital construction of a series of reconstructed images, as described above.
  • the distance between the particle of interest 112A and the reference point 118 is then determined from the series of reconstructed images, this distance being defined along an axis parallel to the axis 135 (step 25). As a variant, it is limited, from the series of reconstructed images, to the distance D ref between the particle of interest 112A and the focusing plane 133.
  • the distance D re f along the optical axis 135 is calculated between the particle of interest 112A and the focusing plane associated with the reference image. Knowing the distance Di along the optical axis 135 between this plane of focus and the reference point 118, the distance DF between the particle of interest 112A and the reference point 118 is deduced along the optical axis 135.
  • each reconstructed image is associated with a value of a useful parameter.
  • the useful parameter being a function of a parameter of the complex intensity (or complex amplitude) of the points belonging to the reconstructed images.
  • a parameter of the complex intensity (or complex amplitude) of a reconstructed image is for example the imaginary part, the real part, the module or the phase of the reconstructed image.
  • these values of the useful parameter are combined in the form of a profile, representing the value of the useful parameter as a function of the offset of each reconstructed image, along the optical axis 135 with respect to the focusing plane 133.
  • a remarkable value is sought on this profile, for example a maximum or a point of inflection or a passage through zero.
  • This remarkable value is associated with the shift along the optical axis 135 between the particle of interest and the position of the focusing plane 133 associated with the reference image.
  • the remarkable value to be searched for is a maximum and generally an extremum, to limit the effect of the uncertainty when passing through a second derivative (detection of a point of inflection), and to overcome a correction of measurement bias (detection of a zero crossing).
  • the remarkable value of the useful parameter can also be a zero crossing of the function describing the evolution of this parameter as a function of the offset between the reconstructed image and the focusing plane 133, or a point of inflection of this parameter. function, or any other criteria.
  • the distance D ref between the focusing plane 133 of the reference image and a particle of interest 112A is determined by the following succession of steps: - identification, on the reference image, of a hologram associated with said particle of interest,
  • the distance D ref between the particle of interest 112A and the focusing plane 133 then being considered equal to this offset.
  • the useful parameter may be the square of the imaginary part of the image, determined from the value of the imaginary part of the different points constituting the image.
  • the useful parameter then corresponds to the average of the square of the imaginary part of the set of points constituting the image. The inventors found that such a parameter made it possible to obtain good location accuracy.
  • the parameter can also include the value of the real part of the points constituting the image or their module, or the square of these different quantities.
  • the inventors have shown that it is possible to locate the position, along the optical axis 135, of a particle of interest 112A, this position corresponding to a remarkable value of the useful parameter, and for example:
  • a zero crossing for example when we observe the evolution, along the optical axis 135, of the real part of the image or its square.
  • the identification of a remarkable value of said parameter of interest makes it possible to conclude that a particle of interest 112A is present in the medium, in the offset range under consideration, the distance between the particle of interest. interest 112A and the focus plane D ref is not necessarily stored.
  • Each particle of interest 112A corresponds to a hologram on the reference image. It is then possible to select a region of interest limited to a hologram, and to apply the previously described steps to each region of interest associated with a hologram.
  • the calculation means 150 output the distance DF, along the axis 135, between the particle of interest 112A and the reference point 118.
  • the use of digital propagation makes it possible to reduce the number of images of the sample to be acquired. In particular, the acquisition of a single image is sufficient.
  • the method according to the invention is therefore particularly fast and automatable.
  • the particle of interest when the particle of interest is in the plane of focus of the optical system, its image can not be distinguished from the image of the surrounding medium (the surrounding medium corresponding to the sample portion 111 surrounding the particle of interest). This is for example the case when the particle of interest is immersed in a medium of the same transmission coefficient or refractive index neighbor (for example to 20%), in the spectrum of the first light source.
  • the method according to the prior art as described in the introduction does not make it possible to determine the position of the particle of interest.
  • the method according to the invention makes it possible to determine the position of a particle along an axis parallel to the optical axis of the optical system, even when the particle can not be distinguished from its surroundings on an image acquired when the particle is found in the development plan of the optical system.
  • the method according to the invention makes it possible to associate an image and a position of the focusing plane along the optical axis 135, this position being determined with great precision (related to the pitch of the offsets implemented at the step 23). Such precision would not be easily achievable by physical, not virtual, offsets.
  • the reference image may comprise the hologram of a single biological particle, which then defines the particle of interest according to the invention.
  • the reference image comprises the hologram of several biological particles, one of these holograms is arbitrarily chosen to define the particle of interest according to the invention, and portions of the reconstructed images, or regions, are used. of interest, centered on one of these holograms.
  • the background image corresponds for example to the image acquired by moving the sample in a plane parallel to the plane of focus, during the image taking.
  • the background image is an average image formed from several static images obtained for several positions of the sample in a plane parallel to the focusing plane.
  • the reference image is divided by the background image. This results in a standardized reference image.
  • Such normalization improves the results of the method.
  • digital propagation is carried out from the reference image used directly, or from the standard reference image as defined above.
  • t p (z) TF- 1 ⁇ TF (U 0 ) XH (z) ⁇
  • U Q the image used to carry out the digital propagation (here the reference image or the standardized reference image), TF the transformed Fourier operator, TF '1 the inverse Fourier transformed operator, X the multiplication term futures of matrices, and H a propagation operator based on the Rayleigh Sommerfeld integral.
  • H (u, v, z) exp [- ⁇ z ⁇ * Im (jp (u, v)) + i * z * Re (p (u, v))] with
  • ⁇ p ⁇ u, v) kx A p x (V l - u 2 - v 2 - l), and k - T
  • N x and N y the number of pixels along x, y, respectively, in the image
  • ⁇ 1 the central wavelength of the first light source.
  • the position of the particle of interest along the optical axis of the optical system is then determined from the imaginary part of each reconstructed image U p (z).
  • the reconstructed image or a reconstructed image portion whose square of the imaginary part (that is to say the average value of the square of the imaginary part of the different points composing the image) is maximum.
  • a reconstructed image portion is a part of a reconstructed image, centered on the image of the particle of interest. If the reference image relates to several biological particles, the same reference image can be used to determine the position of each of these biological particles, each time using reconstructed images centered on a different biological particle.
  • FIGS. 3 and 4 A second method and device embodiment 300 according to the invention will now be described with reference to FIGS. 3 and 4.
  • Figure 3 will only be described for its differences with respect to Figure 1.
  • 120, 150, 151 of FIG. 1 respectively correspond to the reference numerals 311, 320, 312A, 316, 317, 318, 330, 331, 332, 335, 320, 350, 351 of FIG.
  • each of the characteristics distinguishing the device according to FIG. 3 from the device according to FIG. 1 may be isolated and combined independently with the device of FIG. 1 to form many other variants of the invention.
  • the biological particles adhere directly to a second interface of the sample.
  • This second interface is advantageously confused with the first interface as defined above.
  • the position of the first interface relative to the second interface is preferably known, along an axis parallel to the optical axis 135. In particular, this position extending along an axis parallel to the optical axis 135 and passing through the landmark 318.
  • the second interface is substantially orthogonal to the axis 335. This is in particular the lower or upper interface 316, 317 of the sample.
  • the biological particles adhere to the upper interface 317.
  • the upper interface is slightly inclined relative to a plane orthogonal to the axis 335, typically an angle less than 0.1 rad, or even 0.05 rad or 0.02 rad.
  • the method according to the invention makes it possible to overcome the uncertainty in the positioning of the biological particles along the axis 335, when the upper blade, defining the upper interface of the sample, is inclined relative to a plane orthogonal to the axis 335. More generally, the method according to the invention makes it possible to overcome the uncertainty in the positioning of the biological particles along the axis 335, when the upper blade, defining the upper interface of the sample, is deformed relative to a orthogonal plane to the axis 335.
  • This deformation may designate any deformation such that the biological particles adhering to said second interface are located inside a cylinder delimited by two plane surfaces perpendicular to the optical axis 335 and distant from less than 100 ⁇ along this axis, or even less than 50 ⁇ .
  • the distance along an axis parallel to the optical axis 335, between the position of the focusing plane associated with the reference image, and the projection of the reference point 318 along this axis and on the second interface is between +5 ⁇ and +1500 ⁇ or between -5 ⁇ and -1500 ⁇ .
  • this distance is between + 5 ⁇ and +1000 ⁇ ⁇ - 5 ⁇ and -1000 ⁇ , or even between + 5 ⁇ and +800 ⁇ ⁇ - 5 ⁇ and -800 ⁇ , or even between + 5 ⁇ and +200 ⁇ or - 5 ⁇ and -200 ⁇ .
  • Optimal defocusing ranges are thus defined, providing an easier calculation of the distance between the particle of interest and the reference point.
  • the biological particles adhere to said second interface and are located in a cylinder of limited height, it is easy to ensure that when acquiring the reference image, the particle of interest is at a distance from the plane of focus in an optimal defocus range.
  • the width of the optimal defocusing ranges may depend on the characteristics of the first light source 320, in particular its spatial and temporal coherence.
  • the first light source as a function of a distance along the optical axis 335, between several particles of interest that it is desired to study, so that these particles of interest are simultaneously in a range of optimal defocus.
  • a white light source can be used, for an optimum defocusing range of + 5 ⁇ ⁇ +200 ⁇ , and from -5 ⁇ to -200 ⁇ . More the spectral width of the first light source is limited, the width of the optimal defocusing range increases.
  • the analysis device 300 comprises means for determining the distance, along an axis parallel to the optical axis 335, between the reference point and the position of the system debugging plane. optical associated with the reference image. This distance is named useful distance.
  • These means include in particular a laser 360, adapted to illuminate the reference point 318.
  • a shutter 361 connected to a translation plate 362 makes it possible to close off the laser output, when it is not desired for the laser beam 363 to illuminate. landmark 318.
  • the laser beam 363 is advantageously incident on the reference point, along an axis parallel to the optical axis 335.
  • the laser beam 363 and the image sensor 332 are optically coupled to the same optical system 331, so that the focusing point of the laser beam corresponds to the object plane of the sensor, or focusing plane.
  • the reference point is on the upper interface 317, which is why the laser beam 363 enters the sample via this upper interface 317.
  • the laser beam 363 passes through a first separator blade 371, then is reflected on a second separator blade 371, before reaching the upper interface 317.
  • Each separator blade may be a dichroic blade. Alternatively, a cube or a semi-reflective mirror is used.
  • the useful distance is determined by implementing the following substeps illustrated in FIG. 4.
  • a sub-step 41 the reference point 318 is illuminated with the aid of the laser beam. Then, the distance is adjusted, in particular along an axis parallel to the optical axis 335, between the sample and the optical system 331, so that the sensor 332 receives the image of a spot formed by the specular reflection of the image. laser beam 363 at the cue point (substep 42). This adjustment can be achieved thanks to the translation means 340 described with reference to FIG.
  • the definition of the reference point 318 may depend on the arbitrary position of the laser beam 363, in a plane orthogonal to the optical axis 335.
  • This adjustment implements a registration of the specular reflection of the laser beam 363 on the surface receiving the reference point.
  • An example of such an adjustment is shown in FIGS. 5A-5C, the reference point being located on the upper interface 317.
  • the image 5A corresponds to a position of the focusing plane 533 above the upper interface 317.
  • the image obtained on the sensor is a wide and dim spot.
  • the image 5B corresponds to a position of the focusing plane on the upper interface 317.
  • the image obtained on the sensor is a narrow and very bright spot.
  • the image 5C corresponds to a position of the plane of focus below the upper interface 317.
  • the image obtained on the sensor is a wide spot and poor light.
  • the analysis of the evolution of the light intensity of the specular reflection signal of the laser beam, as a function of the relative spacing between the optical system 331 and the sample 311 makes it possible to determine the position of the point of interest 318, the latter corresponding, in this example, to the focal point of the laser beam with a surface capable of reflecting the beam, in this case the face 317 or the face 315.
  • the image of the most intense specular reflection and the image of the specular reflection having a secondary maximum of intensity correspond to the configurations in which the object plane passes through the intersection of the laser beam and respectively, the face of the upper blade 314 opposite the sample, and the face of the upper blade 314 adjacent to the sample.
  • the method may comprise a sub-step 43, consisting in moving the support receiving the sample relative to the imaging means by a known distance and along the optical axis 335.
  • This substep may be useful for placing a particle of interest in an optimal defocusing range as defined above.
  • FIG. 3 shows the laser beam 364 reflected on the upper interface 317, and imaged on the sensor by the optical system 331. The reflected laser beam 364 is reflected on the second separator plate 372, then on the first separator plate 371 .
  • the computing means 351 receive as input an image acquired by the sensor 351, when the shutter 361 is open, and the first light source 320 is off. They control the translation means 340 to implement the substep 42, and if necessary the substep 43.
  • the method according to the invention also comprises a determination, from the series of reconstructed images, of the position of the particle of interest in a plane orthogonal to the optical axis of the optical system.
  • a determination from the series of reconstructed images, of the position of the particle of interest in a plane orthogonal to the optical axis of the optical system.
  • the translation means 340 may implement a translation along two axes together defining a plane orthogonal to the optical axis 335. It is thus possible to locate successively several particles of interest of the same sample. After a single step of determining a landmark, the following series of steps are implemented several times as described above:
  • FIG. 3 also makes it possible to illustrate an additional analysis step implemented in a method and a device according to the invention.
  • the device 300 includes means for positioning the waist of the laser beam 363 (or neck, i.e., location of the narrowest diameter laser beam) on the particle of interest.
  • the radiation emitted by the particle of interest is collected by the optical system 331 and received by a spectrometer 380 which analyzes this radiation.
  • the spectrometer 380 is for example a spectrometer analyzing the Raman scattering or a fluorescence spectrometer. In this type of analysis, it is preferable that the excitation laser beam is centered on the examined particle. This avoids disturbance of the Raman scattering spectrum by the vicinity of the particle.
  • the positioning of the waist of the laser beam 363 is implemented by means of translation means such as a translation plate (not shown), adapted to move the sample relative to the imaging means 331.
  • translation means such as a translation plate (not shown), adapted to move the sample relative to the imaging means 331.
  • These translation means can be formed by the translation means 340.
  • the translation means are controlled by the calculation means 350.
  • the calculation means 350 use the position of the particle of interest relative to the reference point to position the system debugging plane. 331 optics on the particle of interest.
  • the imaging means are moved relative to the particle of interest to focus an analysis laser beam in a plane orthogonal to the optical axis of the imaging means, this plane receiving said particle of interest, and the waist of the analysis laser beam is moved in this plane so as to place it exactly on the particle of interest.
  • a laser source separate from the laser 360 can be used to analyze the particle of interest.
  • a photodetector 332, optically coupled to the optical system 331, may be used, for example to collect a fluorescence signal emitted by the particle of interest 112A in response to an excitation signal.
  • the knowledge of one of the previously mentioned distances allows the focusing of the optical system on the particle, optimizing the collection of the fluorescence signal by the photodetector 332.
  • a medium has a plurality of particles of interest.
  • the distance Dref or DF corresponding to each particle of interest 112A of the sample is determined and stored.
  • the relative positioning of the optical system 331 is adjusted so that the focusing plane of this optical system comprises the particle of interest. This adjustment is made, successively, for each particle of interest, so as to optimize the analysis.
  • FIG. 6 illustrates an example of a reference image according to the invention.
  • the abscissa and ordinate axes are graduated in pixels.
  • FIG. 7 a detailed example of the steps used to determine, from a first series of reconstructed images, the distance between the particle of interest and the reference point, according to FIG. an axis parallel to the optical axis of the optical system.
  • H (u, v, z) exp [- ⁇ z ⁇ * Im (j> (u, v)) + i * z * Re (p (u, v))]
  • Each useful reconstructed image corresponds to the same matrix of pixels. Depending on the useful reconstructed image, the intensity associated with each pixel varies.
  • An image of the gradients 720 corresponding to the difference between the image of the maxima and the image of the minima is then calculated.
  • This step can be carried out using a first matrix, one dimension of which corresponds to the pixels, and the other dimension corresponds to the shift of the focusing plane associated with the useful reconstructed image.
  • This offset is measured along the optical axis of the optical system, named axis (Oz) in the following.
  • the first matrix has XY columns and N rows.
  • the image of the minima then corresponds to a second matrix XY columns and 1 line.
  • the image of maxima corresponds to a third matrix XY columns and 1 line.
  • the image of the gradients corresponds to a fourth matrix XY columns and 1 line, equal to the difference between the third and the second matrices.
  • Pixels with a strong intensity gradient along the axis (Oz) are selected.
  • An average gradient associated with these pixels is defined.
  • the position of the highest values of this average gradient roughly defines the position of the biological particles along the axis (Oz).
  • the columns with the highest values can be selected in the fourth matrix.
  • Figure 730 illustrates the fifth matrix in the form of a profile in which the x-axis is an offset along the axis (Oz), and the y-axis is the corresponding value in the fifth matrix.
  • Sub-step 73 defines an average position along the axis (Oz) of several biological particles imaged on the reference image. This average position defines an approximate position of a particle of interest according to the invention, assuming that the biological particles are located substantially in the same plane orthogonal to the axis (Oz).
  • a gradient is used along the axis (Oz) associated with a selection of pixels each having a strong gradient along this axis.
  • the maximum of the gradient associated with this selection is therefore all the easier to spot.
  • One of the images of the first series of useful reconstructed images is selected, that associated with the position along the axis (Oz), calculated in the substep 73.
  • a thresholding is performed to obtain a binary image 740.
  • the coordinates of the binary objects highlighted by the thresholding are recorded.
  • An approximate position of each biological particle is thus calculated in a plane orthogonal to the axis (Oz).
  • This step may comprise a step of interpolating each binary object by an ellipse, in order to further calculate an approximate geometrical shape of each binary object, and therefore of each biological particle.
  • the particle of interest is arbitrarily chosen by arbitrarily choosing one of the binary objects.
  • An approximate position of the particle of interest has therefore been determined in a plane orthogonal to the axis (Oz).
  • a region of interest comprising this single binary object is defined on the binary image 740.
  • the region of interest is defined by pixel positions.
  • the region of interest is for example a square of 16 * 16 pixels.
  • This step is implemented from the first series of useful reconstructed images, by selecting in each useful reconstructed image a region corresponding to the region of interest as defined above.
  • this step is implemented from a second series of useful reconstructed images.
  • Each image of the second series of useful reconstructed images is formed by the square of the imaginary part of an image of a second series of reconstructed images.
  • the second series of reconstructed images is associated with a sampling step along the axis (Oz), less than the sampling interval of the first series of reconstructed images.
  • Each image of the second series of reconstructed images corresponds to the region of interest as defined above.
  • the amplitude of the sampling range along the axis (Oz) associated with the second series of reconstructed images is less than the amplitude of the sampling range associated with the first series of reconstructed images.
  • the second set of reconstructed images can be calculated only for finding a precise position of the particle of interest.
  • Sub-steps 72 and 73 are implemented on the second set of reconstructed useful images to calculate a precise position of the particle of interest along the axis (Oz).
  • Figure 750 corresponds to Figure 730, replacing the first set of reconstructed useful images with the second set of useful reconstructed images.
  • a precise position of the particle of interest is calculated in a plane orthogonal to the axis (Oz). For this, we select from the images of the second series of useful reconstructed images, that associated with the precise position of the particle of interest along the axis (Oz).
  • a thresholding is performed to obtain a new binary image. Then, we find the coordinates of the center of the binary object highlighted by the thresholding.
  • a useful reconstructed image is then defined as being the squared module of the corresponding reconstructed image.
  • Each useful reconstructed image is then defined as the squared module of the corresponding reconstructed image.
  • An image of the gradients is then formed as defined above, and the approximate positions of the biological particles are deduced, by thresholding, in a plane orthogonal to the axis (Oz).
  • the image obtained here by thresholding is named the first thresholded image.
  • the approximate position of the particle of interest is defined in a plane orthogonal to the axis (Oz).
  • first binary object associated with a biological particle.
  • a first region of interest receiving this single binary object is defined.
  • a region corresponding to this first region of interest is selected in each useful reconstructed image.
  • the position is sought along the axis (Oz) of the image having a maximum value of standard deviation. This position defines an approximate position of the particle of interest along the axis (Oz), assuming that all the biological particles are located substantially in the same plane orthogonal to the axis (Oz).
  • the particle of interest is defined by arbitrarily choosing a second binary object. By searching for the center of this second binary object, we deduce the precise position of the particle of interest, in a plane orthogonal to the axis (Oz).
  • the distance along this axis is calculated between the particle of interest and the focusing plane associated with the axis. reference image.
  • the calculation makes it possible to determine the absolute value of the distance, along an axis parallel to the axis 135, between the particle of interest and the focusing plane associated with the reference image. Then we choose one or the other sign knowing the experimental setup.
  • a series of secondary images is determined, each secondary image being formed by the real part of a corresponding reconstructed image; in the same way as in the sub-step 72 illustrated in FIG. 7, the series of secondary images (three-dimensional matrix) is transformed into a two-dimensional matrix, one dimension corresponding to the pixels, and the other dimension corresponds to an offset along the axis (Oz).
  • the pixels having a strong intensity gradient along the axis (Oz) are selected.
  • An average gradient associated with these pixels is thus defined. It is this gradient which is illustrated in FIG. 8.
  • the abscissa axis is an offset along the axis (Oz), in particular a shift of the focusing plane relative to the position of the focusing plane associated with the reference image.
  • the y-axis refers to the real part of the secondary images.
  • the sign of the distance along an axis parallel to the axis 135 is determined between the particle of interest and the focusing plane associated with the reference image:
  • the reference image is below the focusing plane (in the direction of propagation of the light).
  • the digital construction step of a series of reconstructed images may implement other propagation operators than those exemplified.
  • the series of reconstructed images may also be exploited in different ways, to deduce the distance between the focus plane associated with the reference image, and the particle of interest.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Pathology (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Dispersion Chemistry (AREA)
  • Computing Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • Computer Graphics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

Procédé d'analyse holographique d'un échantillon (111) incluant une particule d'intérêt (112A) disposé entre une source lumineuse (120) et un système optique (131) et comprenant les étapes suivantes : - définition d'un point de repère (118) situé sur une interface (116, 117) de l'échantillon, ou à une distance connue de celle-ci; - positionnement du plan objet (133) du système optique à une distance contrôlée D1 du point de repère, la particule d'intérêt (112A) se trouvant hors dudit plan objet (133); - acquisition d'une image holographique en transmission de l'échantillon; - reconstruction numérique d'images, associées chacune à un décalage virtuel prédéterminé (Dii) du plan objet le long de l'axe optique; et - déduction de la distance (D F) entre la particule d'intérêt et le point de repère.

Description

PROCEDE D'ANALYSE COMPRENANT LA DETERMINATION
HOLOGRAPHIQUE 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. 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é étant mis en œuvre 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 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 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 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.
L'invention comprend donc avantageusement une étape de détection du point de repère, par un placement du système optique de façon à ce que son plan de mise au point soit placé sur le point de repère, sur ladite première interface ou à la distance connue celle-ci, et disposition du système optique de telle sorte que la particule d'intérêt se trouve hors du plan de mise au point du système optique.
L'image de référence est acquise pour cette disposition du système optique, telle que la particule d'intérêt se trouve hors du plan de mise au point du système optique.
Cette disposition du système optique peut correspondre directement à la disposition du système optique mise en œuvre pour détecter le point de repère. Dans ce cas, l'image de référence, ou image défocalisée, est acquise directement après détection du point de repère.
En variante, après détection du point de repère, le procédé comprend une étape de décalage du plan de mise au point du système optique relativement audit point de repère, par une translation d'un support recevant l'échantillon relativement au système optique. L'image de référence, ou image défocalisée, est acquise après cette étape de décalage. L'invention permet de connaître avec précision une distance entre le plan de mise au point et ledit point de repère. On peut ainsi connaître avec précision une distance entre la particule d'intérêt et le plan de mise au point, ce qui permet ensuite de positionner au mieux un dispositif d'analyse pour étudier la particule d'intérêt. 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.
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 : - 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 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 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 μιη et +2000 μιη ou entre -5 μιη et -2000 μιη.
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 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 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 à 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 œuvre 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 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 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 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 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 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 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 ;
- 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
- 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,
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, 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.
Le dispositif peut comprendre également un laser, adapté à fournir un faisceau laser aligné avec l'axe optique du système optique, et couplé au système optique de sorte que le point de focalisation du faisceau laser correspond au plan de mise au point. 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, 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 ;
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 ; 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. 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 œuvre du procédé d'analyse selon l'invention.
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 œuvre dans le dispositif d'analyse 100.
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.
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.
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 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 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 ensem ble les 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.
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 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.
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.
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, 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.
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. Le capteur 132 est par exemple un capteur matriciel de type CCD ou CMOS. I l 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 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 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é.
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. 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 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. Le dispositif 100 est adapté à la mise en œuvre 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 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 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 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 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).
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 Dl 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 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 μιη et 1500 μιη au-dessus ou en dessous (en amont ou en aval) du point de repère, de préférence entre 5 μιη et 1000 μιη, voire entre 5 μιη et 800 μιη. 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 supérieure à 2 μιη, et de préférence supérieure à 5 μιη, de préférence dans la plage [5 μιη - 1500 μιη] voire [5μιη - 1000 μιη], voire [5μιη - 800μιη] d'une interface délimitant l'échantillon 111. En variante, on dispose d'une information a priori qua nt à 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.
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.
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 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. Dans une étape 24, on met en œuvre 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 décalages virtuels, c'est-à-dire simulés par la construction numérique, et non pas mis en œuvre physiquement. En d'autres termes, chaque image reconstruite est calculée en un plan de reconstruction, associé à un décalage prédéterminé par rapport au plan de mise au point. Plus précisément, chaque plan de reconstruction correspond à l'image, par le système optique, d'un plan objet décalé du plan de mise au point de ce décalage prédéterminé. On parle généralement de propagation numérique. Ces décalages sont notés Dii, Di2, ... Din, DU, ... 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 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 μιη et 1 μιη, par exemple 0,25 μιη ou 0,20 μιη. 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 Dii, Di2, ... Din, DU, ... Di-n de la position du plan de mise au point associée à 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 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. 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.
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 holographie microscopy using Rayleigh-Sommerfeld back-propagation », 16 July 2012 / Vol. 20, No 15 / OPTICS EXPRESS, 16735-16744, Wilson et al. décrit un exemple d'une propagation numérique mise en œuvre à 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 « Holographie microscopy of holographically trapped three-dimensinal structures », 19 February 2007 / Vol. 15, No. 4 / OPTICS EXPRESS 1505-1512.
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 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 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 Di 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 DF entre la particule d'intérêt 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 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 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 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 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.
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,
- 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 Dii, Di2, ... Din, DU, ... Di-n prédéterminés par rapport au plan de mise au 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 utile,
- analyse de l'évolution de ce paramètre utile en fonction dudit décalage et identification d'une valeur remarquable de ce pa ramè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 mise au point 133 étant alors considérée comme égale à ce décalage.
Connaissant la distance Di entre le plan de mise au point 133 et le point de repère 118, il est alors possible de déduire la distance DF entre la particule d'intérêt 112A et le point de repère 118, par l'opération DF = Dref - Di
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, 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.
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 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é 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 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 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 DF, selon l'axe 135, entre la particule d'intérêt 112A et le point de repère 118.
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 automatisable. Dans certains cas, lorsque la particule d'intérêt se trouve dans le plan de mise 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 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 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.
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 œuvre à l'étape 23). Une telle précision ne serait pas facilement atteignable à l'aide de décalages physiques, et non virtuels.
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'hologram me 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 reconstruites, ou régions d'intérêt, centrées sur l'un de ces hologrammes.
On peut mettre en œuvre 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 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 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 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 :
t p(z) = TF-1{TF(U0) X H(z)}
avec UQ 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 Fourier, TF'1 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[—\z\ * Im(jp(u, v)) + i * z * Re(p(u, v))] avec
■ Re l'opérateur partie réelle, Ira l'opérateur partie imaginaire, i2 =— 1,
■ u et v 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,
2π
■ p{u, v) = k x Ap x(V l - u2 - v2 - l), et k -- T
Les coordonnées u et v sont définies de la façon suivante
λ -2- . _ _L . . Nx - 1 par pas de
A p ' A p A, { Nx
À ( 0 À r i A { Ny ~ l _ par pas de
{"> 2 2 Ny x Ap
J '' P 2 ;- ; Ap
J { "y J et avec
■ Δρ 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 a u facteur de grandissement du système optique 131),
■ Nx et Ny les nombres de pixels selon x, respectivement y, dans l'image
U0,
■ 1 la longueur d'onde centrale de la première source lumineuse.
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 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 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,
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 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.
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 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.
La deuxième interface est sensiblement orthogonale à l'axe 335. I l 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
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.
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 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 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 100 μιη selon cet axe, voire moins de 50 μιη.
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 μιη et +1500 μιη ou entre -5 μιη et -1500 μιη. Avantageusement, cette distance est comprise entre + 5 μιη et +1000 ίη ου - 5 μιη et -1000 μιη, voire entre + 5 μιη et +800 μίη ου - 5 μιη et -800 μιη, ou même entre + 5 μιη et +200 μιη ou - 5 μιη et -200 μιη. 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 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 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.
Par exemple, on peut utiliser une source de lumière blanche, pour une plage de défocalisation optimale allant + 5 μπΊ θ +200 μιη, et de -5 μιη à -200 μιη. Plus on restreint la largeur spectrale de la première source lumineuse, plus la largeur de la plage de défocalisation optimale augmente.
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, 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 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 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 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 œuvre les sous-étapes suivantes illustrées en figure 4.
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 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.
Cet ajustement met en œuvre 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 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.
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 ;
une première tâche é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 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
une tâche 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 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.
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 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.
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 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.
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. I ls pilotent les moyens de translation 340 pour mettre en œuvre la sous-étape 42, et le cas échéant la sous-étape 43.
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.
Les moyens de translation 340 peuvent mettre en œuvre 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 œuvre plusieurs fois les séries d'étapes suivantes telles que décrites ci- 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 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. La figure 3 permet également d'illustrer une étape supplémentaire d'analyse mise en œuvre 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.
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 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 œuvre grâce à des moyens de translation tels qu'une platine de translation (non représentée), 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 particule d'intérêt.
En particulier, les moyens d'imagerie sont déplacés par rapport à la particule d'intérêt pour focaliser un faisceau laser d'analyse dans un plan orthogonal à l'axe optique des moyens d'imagerie, ce plan recevant ladite particule d'intérêt, et le waist du faisceau laser d'analyse est déplacé dans ce plan de façon à le placer exactement sur la 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 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 (DF OU Dref) permet la focalisation du système optique sur la particule, optimisant la collecte du signal de fluorescence par le photodétecteur 332.
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 DF 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 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 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 œuvre 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.
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[—\z\ * Im(j>(u, v)) + i * z * Re(p(u, v))]
On met ensuite en œuvre les sous-étapes suivantes :
• 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.
• 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.
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 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, 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.
• 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).
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 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 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 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é à 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. 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 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 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 œuvre cette étape à partir de la première série d'images reconstruites utiles, en sélectionnant dans chaque image reconstruite utile une région correspondant à la région d'intérêt telle que définie ci-avant.
En variante, on met en œuvre 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.
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. 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.
On met en œuvre 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.
• 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 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.
En particulier, on peut se placer dans le mode de réalisation tel que décrit ci-avant, dans lequel on a :
H(u, v, z) = exp(i * 2 * π * z * w), avec w = pour
u2 + v2 < T
X2
H(u, v, z) = 0 sinon. On définit alors une image reconstruite utile comme étant le module au carré de l'image reconstruite correspondante.
Ensuite, on met en œuvre des étapes similaires à celles décrites ci-avant, à la différence près que la recherche d'une position selon l'axe (Oz) met en œuvre la recherche d'une image reconstruite utile présentant une valeur maximale d'écart- type.
Selon une autre variante, on se place dans le mode de réalisation tel que décrit ci-avant, dans lequel on a :
H(u, v, z) = exp(i * 2 * π * z * w), avec w = \—— u2— v2, pour u2 + v2 <
λ2
H(u, v, z) = 0 sinon. On définit alors chaque image reconstruite utile comme étant le module au carré de l'image reconstruite correspondante.
On forme ensuite une image des gradients telle que définie ci-avant, et on en déduit, par seuillage, les positions approchées des particules biologiques, dans un plan orthogonal à l'axe (Oz). L'image obtenue ici par seuillage est nommée première image seuillée. On définit ainsi notamment la position approchée de la particule d'intérêt, dans un plan orthogonal à l'axe (Oz).
Sur la première image seuillée, on choisit arbitrairement un premier objet binaire associé à une particule biologique. On définit une première région d'intérêt recevant ce seul objet binaire.
Ensuite, on sélectionne une région correspondant à cette première région d'intérêt, dans chaque image reconstruite utile. A partir de la série d'images de la première région d'intérêt ainsi réalisée, on recherche la position selon l'axe (Oz) de l'image présentant une valeur maximale d'écart-type. Cette position définit une position approchée de la particule d'intérêt, selon l'axe (Oz), en supposant que toutes les particules biologiques sont situées sensiblement dans un même plan orthogonal à l'axe (Oz).
Ensuite, on revient à l'image des gradients et on forme une deuxième image seuillée, avec un seuil supérieur au seuil utilisé pour former la première image seuillée (par exemple 1,5 fois supérieur). Sur la deuxième image seuillée, on définit la particule d'intérêt en choisissant arbitrairement un deuxième objet binaire. Par recherche du centre de ce deuxième objet binaire, on en déduit la position précise de la particule d'intérêt, dans un plan orthogonal à l'axe (Oz).
A partir de cette deuxième image seuillée, on définit une deuxième région d'intérêt recevant ce seul deuxième objet binaire. Ensuite, on sélectionne une région correspondant à cette deuxième région d'intérêt, dans chaque image reconstruite utile. A partir de la série d'images ainsi réalisée, on recherche la position selon l'axe (Oz) de l'image présentant une valeur maximale d'écart-type. Cette position définit une position précise de la particule d'intérêt, selon l'axe (Oz). La figure 8 est utilisée pour illustrer un cas particulier de l'invention. On se réfère aux références numériques de la figure 1.
Pour déterminer la distance entre la particule d'intérêt et le point de repère, selon un axe parallèle à l'axe 135, on calcule la distance selon cet axe, entre la particule d'intérêt et le plan de mise au point associé à l'image de référence.
En pratique, le calcul permet de déterminer la valeur absolue de la distance, selon un axe parallèle à l'axe 135, entre la particule d'intérêt et le plan de mise au point associé à l'image de référence. On choisit ensuite l'un ou l'autre signe connaissant le montage expérimental.
Par exemple,
si le point de repère est situé au-dessus de l'interface supérieure 117, et
si le plan de mise au point associé à l'image de référence est situé au- dessus de l'interface supérieure 117,
alors on sait que la particule d'intérêt est sous le plan de mise au point associé à l'image de référence.
Lorsqu'il n'est pas possible de déterminer avec certitude le signe de cette distance, on peut utiliser un profil selon l'axe (Oz) relatif à la partie réelle des images reconstruites. On se place alors dans le mode de réalisation tel que décrit ci-avant, dans lequel on a :
H(u, v, z) = exp[—\z\ * Im(j>(u, v)) + i * z * Re(p(u, v))] En particulier, on peut mettre en œuvre les étapes suivantes :
à partir de la série d'images reconstruites, on détermine une série d'images secondaires, chaque image secondaire étant formée par la partie réelle d'une image reconstruite correspondante ; de la même façon qu'à la sous-étape 72 illustrée en figure 7, on transforme la série d'images secondaires (matrice à trois dimensions), en une matrice à deux dimensions, dont une dimension correspond aux pixels, et l'autre dimension correspond à un décalage selon l'axe (Oz).
dans ladite matrice à deux dimensions, on sélectionne les pixels présentant un fort gradient d'intensité selon l'axe (Oz). On définit ainsi un gradient moyen associé à ces pixels. C'est ce gradient qui est illustré en figure 8. L'axe des abscisses est un décalage selon l'axe (Oz), en particulier un décalage du plan de mise au point relativement à la position du plan de mise au point associée à l'image de référence. L'axe des ordonnées se rapporte à la partie réelle des images secondaires.
en fonction de la forme du profil, on détermine le signe de la distance, selon un axe parallèle à l'axe 135, entre la particule d'intérêt et le plan de mise au point associé à l'image de référence :
si on passe, pour une abscisse croissante, d'un maximum positif à un minimum négatif (cas de la figure 8), l'image de référence est au-dessus du plan de mise au point (dans le sens de propagation de la lumière),
si on passe, pour une abscisse croissante, d'un minimum négatif à un maximum positif (cas contraire de la figure 8), l'image de référence est au-dessous du plan de mise au point (dans le sens de propagation de la lumière).
L'invention n'est pas limitée aux exemples précédemment cités, et de nombreuses variantes pourront être envisagées sans sortir du cadre de la présente invention. Par exemple, l'étape de construction numérique d'une série d'images reconstruites pourra mettre en œuvre d'autres opérateurs de propagation que ceux cités en exemple. La série d'images reconstruites pourra également être exploitée de différentes façons, pour en déduire la distance entre le plan de mise au point associé à l'image de référence, et la particule d'intérêt.

Claims

REVENDICATIONS
1. 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 œuvre 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 ;
- détection de ce point de repère, par un placement du système optique (131 ; 331) de façon à ce que son plan de mise au point (133 ; 333) soit placé sur ladite première interface ou à la distance connue de celle-ci, et disposition du système optique de telle sorte que la particule d'intérêt se trouve hors du plan de mise au point 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 holographique de la région illuminée, dite image de référence ou image défocalisée, 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 (Dl) ;
- à partir de l'image de référence, construction numérique d'une série d'images reconstruites, associées chacune à une simulation d'un décalage prédéterminé (Dii, Di2, Din, DU, 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 distan (Dref) selon un axe parallèle à l'axe optique du système optique, entre 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 étape de décalage du plan de mise au point (133 ; 333) du système optique relativement audit point de repère, par une translation d'un support (110) recevant l'échantillon relativement au système optique (131 ; 331).
3. Procédé d'analyse selon la revendication 1 ou 2, caractérisé en ce qu'il comporte une détermination de la distance (DF) entre la particule d'intérêt et ledit point de 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 (Di) .
4. Procédé d'analyse selon l'une quelconque des revendications 1 à 3, caractérisé en 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).
5. Procédé d'analyse selon l'une des revendications 1 à 4, caractérisé en ce que la détermination de la distance utile (Dl), 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) ;
- 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 (Dl) selon un axe parallèle à l'axe optique du système optique, entre l'échantillon (111 ; 311) et le système optique (131 ; 331), en fonction de l'intensité de l'image formée sur le capteur (132 ; 332).
6. Procédé d'analyse selon l'une des revendications 1 à 5, caractérisé en ce que les particules biologiques (112) adhèrent sur une deuxième interface (116, 117 ; 316, 317) de l'échantillon, constituée d'une interface inférieure (116 ; 316) ou supérieure (117 ; 317) de l'échantillon, la première et la deuxième interfaces étant distinctes ou confondues.
7. Procédé d'analyse selon l'une des revendications 1 à 6, caractérisé en ce que la distance (Dl) 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 (118 ; 318) selon cet axe et sur la deuxième interface, est comprise entre +5 μιη et +2000 μιη ou entre -5 μιη et -2000 μιη.
8. Procédé d'analyse selon l'une quelconque des revendications 1 à 7, caractérisé en ce que la région illuminée (119) com prend une pluralité de particules biologiques (112).
9. Procédé d'analyse selon l'une quelconque des revendications 1 à 8, caractérisé en ce que la première source lumineuse (120 ; 320) présente une largeur spectrale inférieure à 200 nm.
10. Procédé d'analyse selon l'une quelconque des revendications 1 à 9, caractérisé en ce qu'il comprend en outre une détermination, à partir de ladite 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.
11. Procédé d'analyse selon l'une quelconque des revendications 1 à 10, caractérisé en ce que chaque image reconstruite est formée par une partie réelle 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 une distance (DF) 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).
12. Procédé d'analyse selon l'une des revendications 1 à 11, 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 (112A, 312A) et le plan de mise au point (133, 333) met en œuvre une recherche d'une valeur remarquable, en particulier un extremum, un point d'inflexion ou un passage par zéro de ladite fonction.
13. Procédé d'analyse selon l'une quelconque des revendications 1 à 12, 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 simulés du plan de mise au point (133 ; 333), 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 deuxième pas des décalages simulés 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 entre la particule d'intérêt et ledit plan de mise au point.
14. Procédé d'analyse selon l'une quelconque des revendications 1 à 13, 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 déplacé par rapport à ladite particule d'intérêt (112A, 312A) de façon à focaliser un 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 (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 .
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).
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 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 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) 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 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 holographique 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, associées chacune à une simulation d'un décalage prédéterminé (Dii ; Di2 ; Din ; DU ; Di-n) du plan de mise au point le long de l'axe optique du 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).
EP15801133.8A 2014-11-13 2015-11-13 Procede d'analyse comprenant la determination holographique d'une position d'une particule biologique. Active EP3218769B1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1460947A FR3028616A1 (fr) 2014-11-13 2014-11-13 Procede d'analyse comprenant la determination d'une position d'une particule biologique.
PCT/EP2015/076530 WO2016075279A1 (fr) 2014-11-13 2015-11-13 Procede d'analyse comprenant la determination holographique d'une position d'une particule biologique

Publications (2)

Publication Number Publication Date
EP3218769A1 true EP3218769A1 (fr) 2017-09-20
EP3218769B1 EP3218769B1 (fr) 2019-03-20

Family

ID=52423893

Family Applications (1)

Application Number Title Priority Date Filing Date
EP15801133.8A Active EP3218769B1 (fr) 2014-11-13 2015-11-13 Procede d'analyse comprenant la determination holographique d'une position d'une particule biologique.

Country Status (6)

Country Link
US (1) US10255688B2 (fr)
EP (1) EP3218769B1 (fr)
JP (1) JP6649379B2 (fr)
CN (1) CN107003229B (fr)
FR (1) FR3028616A1 (fr)
WO (1) WO2016075279A1 (fr)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10401793B2 (en) * 2010-06-17 2019-09-03 Purdue Research Foundation Digital holographic method of measuring cellular activity and measuring apparatus with improved stability
FR3044415B1 (fr) * 2015-11-27 2017-12-01 Biomerieux Sa Procede de determination de la reaction d'un microorganisme a son exposition a un antibiotique
WO2017130613A1 (fr) * 2016-01-27 2017-08-03 株式会社 日立ハイテクノロジーズ Dispositif d'observation
FR3049348B1 (fr) * 2016-03-23 2023-08-11 Commissariat Energie Atomique Procede de caracterisation d’une particule dans un echantillon
EP3252455A1 (fr) * 2016-05-30 2017-12-06 Biomérieux Dispositif et procede d'acquisition d'une particule presente dans un echantillon
FR3075370B1 (fr) 2017-12-15 2021-03-05 Biomerieux Sa Procede de calibration d'un dispositif d'analyse et dispositif associe
JP7054343B2 (ja) * 2017-12-26 2022-04-13 川崎重工業株式会社 分注装置及び分注方法
US10895843B2 (en) * 2018-03-02 2021-01-19 Regents Of The University Of Minnesota System, devices, and methods for digital holography
EP3575773A1 (fr) * 2018-05-28 2019-12-04 Universität für Bodenkultur Wien Procédé permettant de déterminer une distribution de particules tridimensionnelles dans un milieu
CN110044898A (zh) * 2019-03-29 2019-07-23 江苏大学 一种基于反射增强衍射的全息病害孢子检测装置与方法
CN113495470A (zh) * 2020-03-18 2021-10-12 中央大学 使用全像光学元件的无透镜数字全像取像系统
EP3922979B1 (fr) * 2020-06-08 2023-06-07 Universität für Bodenkultur Wien Procédé de détermination d'une longueur focale d'une particule dans un milieu

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8331019B2 (en) * 2007-01-26 2012-12-11 New York University Holographic microscopy of holographically trapped three-dimensional nanorod structures
EP2329324B1 (fr) * 2008-10-03 2014-01-29 Universite Libre De Bruxelles Procédé pour étudier des objets nanométriques utilisant la microscopie holographique
FR2958298B1 (fr) * 2010-04-06 2014-10-17 Commissariat Energie Atomique Procede de detection d'amas de particules biologiques
US10921255B2 (en) * 2014-12-09 2021-02-16 Bioaxial Sas Optical measuring device and process
FR3030749B1 (fr) * 2014-12-19 2020-01-03 Commissariat A L'energie Atomique Et Aux Energies Alternatives Methode d'identification de particules biologiques par piles d'images holographiques defocalisees
FR3031180A1 (fr) * 2014-12-30 2016-07-01 Commissariat Energie Atomique Methode et dispositf de controle de la position d'une particule transparente par pile d'images holographiques
US10690898B2 (en) * 2016-09-15 2020-06-23 Molecular Devices (Austria) GmbH Light-field microscope with selective-plane illumination

Also Published As

Publication number Publication date
JP6649379B2 (ja) 2020-02-19
US10255688B2 (en) 2019-04-09
WO2016075279A1 (fr) 2016-05-19
EP3218769B1 (fr) 2019-03-20
CN107003229B (zh) 2020-02-18
JP2018502283A (ja) 2018-01-25
CN107003229A (zh) 2017-08-01
FR3028616A1 (fr) 2016-05-20
US20170309036A1 (en) 2017-10-26

Similar Documents

Publication Publication Date Title
EP3218769B1 (fr) Procede d&#39;analyse comprenant la determination holographique d&#39;une position d&#39;une particule biologique.
EP3274689B1 (fr) Procédé et dispositif d&#39;analyse de particules
EP3234550B1 (fr) Méthode d&#39;identification de particules biologiques par piles d&#39;images holographiques défocalisées
EP3274694B1 (fr) Procédé de détermination de l&#39;état d&#39;une cellule
EP3241014B1 (fr) Système d&#39;analyse d&#39;un échantillon transparent avec contrôle de position, et procédé associé
EP3465153B1 (fr) Dispositif et procede d&#39;acquisition d&#39;une particule presente dans un echantillon
FR3049348A1 (fr) Procede de caracterisation d’une particule dans un echantillon
EP3519899B1 (fr) Dispositif d&#39;observation d&#39;un échantillon et procédé d&#39;observation d&#39;un échantillon
EP3559631B1 (fr) Procédé de numération de particules dans un échantillon par imagerie sans lentille
EP3637194B1 (fr) Procédé de détermination de paramètres d&#39;une particule
EP3584560A1 (fr) Procédé d&#39;observation d&#39;un échantillon par imagerie sans lentille, avec prise en compte d&#39;une dispersion spatiale dans l&#39;échantillon
EP3397944B1 (fr) Dispositif et procédé d&#39;observation bimodale d&#39;un objet
WO2018060589A1 (fr) Procédé de numération de leucocytes dans un échantillon
WO2021260321A1 (fr) Procédé de mise au point pour système d&#39;imagerie holographique
EP3724725B1 (fr) Procede de calibration d&#39;un dispositif d&#39;analyse et dispositif associe
EP3462381A1 (fr) Procédé de détection de microorganismes dans un échantillon
FR3081552A1 (fr) Dispositif et procede d&#39;observation d&#39;un echantillon fluorescent par imagerie defocalisee
EP3545362B1 (fr) Procédé de formation d&#39;une image de haute résolution par imagerie sans lentille
WO2019243725A1 (fr) Procede et dispositif de comptage de thrombocytes dans un echantillon
EP4307051A1 (fr) Procede et systeme de caracterisation de microorganismes par microscopie holographique numerique
FR3085474A1 (fr) Methode de mesure de la rugosite de ligne d&#39;un reseau de lignes par diffraction de rayons x

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20170517

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
REG Reference to a national code

Ref country code: DE

Ref legal event code: R079

Ref document number: 602015026856

Country of ref document: DE

Free format text: PREVIOUS MAIN CLASS: G03H0001080000

Ipc: G03H0001040000

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

RIC1 Information provided on ipc code assigned before grant

Ipc: G03H 1/04 20060101AFI20180911BHEP

INTG Intention to grant announced

Effective date: 20181010

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602015026856

Country of ref document: DE

REG Reference to a national code

Ref country code: AT

Ref legal event code: REF

Ref document number: 1111155

Country of ref document: AT

Kind code of ref document: T

Effective date: 20190415

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: FRENCH

REG Reference to a national code

Ref country code: NL

Ref legal event code: MP

Effective date: 20190320

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: NO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190620

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG4D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190621

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190620

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 1111155

Country of ref document: AT

Kind code of ref document: T

Effective date: 20190320

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190720

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190720

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602015026856

Country of ref document: DE

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

26N No opposition filed

Effective date: 20200102

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20191113

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20191130

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20191130

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20191130

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20191113

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20191130

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20151113

Ref country code: MT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190320

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20231123

Year of fee payment: 9

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20231120

Year of fee payment: 9

Ref country code: DE

Payment date: 20231120

Year of fee payment: 9