EP1740974A1 - Procede et systeme de mesure de vitesse du flux sanguin - Google Patents

Procede et systeme de mesure de vitesse du flux sanguin

Info

Publication number
EP1740974A1
EP1740974A1 EP05744612A EP05744612A EP1740974A1 EP 1740974 A1 EP1740974 A1 EP 1740974A1 EP 05744612 A EP05744612 A EP 05744612A EP 05744612 A EP05744612 A EP 05744612A EP 1740974 A1 EP1740974 A1 EP 1740974A1
Authority
EP
European Patent Office
Prior art keywords
streak
speed
image
sin
streaks
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
EP05744612A
Other languages
German (de)
English (en)
Other versions
EP1740974B1 (fr
Inventor
François LACOMBE
Georges Le Goualher
Aymeric Perchant
Nicholas Ayache
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.)
Mauna Kea Technologies SA
Original Assignee
Mauna Kea Technologies SA
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 Mauna Kea Technologies SA filed Critical Mauna Kea Technologies SA
Publication of EP1740974A1 publication Critical patent/EP1740974A1/fr
Application granted granted Critical
Publication of EP1740974B1 publication Critical patent/EP1740974B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
    • A61B3/1241—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes specially adapted for observation of ocular blood flow, e.g. by fluorescein angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/1025—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for confocal scanning
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/704—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
    • G01F1/708—Measuring the time taken to traverse a fixed distance
    • G01F1/7086—Measuring the time taken to traverse a fixed distance using optical detecting arrangements
    • 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
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1468—Optical investigation techniques, e.g. flow cytometry with spatial resolution of the texture or inner structure of the particle
    • G01N15/147—Optical investigation techniques, e.g. flow cytometry with spatial resolution of the texture or inner structure of the particle the analysis being performed on a sample stream
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00—Microscopes
    • G02B21/0004—Microscopes specially adapted for specific applications
    • G02B21/002—Scanning microscopes
    • G02B21/0024—Confocal scanning microscopes (CSOMs) or confocal "macroscopes"; Accessories which are not restricted to use with CSOMs, e.g. sample holders
    • G02B21/0028—Confocal scanning microscopes (CSOMs) or confocal "macroscopes"; Accessories which are not restricted to use with CSOMs, e.g. sample holders specially adapted for specific applications, e.g. for endoscopes, ophthalmoscopes, attachments to conventional microscopes
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00—Microscopes
    • G02B21/0004—Microscopes specially adapted for specific applications
    • G02B21/002—Scanning microscopes
    • G02B21/0024—Confocal scanning microscopes (CSOMs) or confocal "macroscopes"; Accessories which are not restricted to use with CSOMs, e.g. sample holders
    • G02B21/008—Details of detection or image processing, including general computer control
    • G02B21/0084—Details of detection or image processing, including general computer control time-scale detection, e.g. strobed, ultra-fast, heterodyne detection
    • 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/1027—Determining speed or velocity of a particle
    • 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/1447—Spatial selection
    • G01N2015/145—Spatial selection by pattern of light, e.g. fringe pattern

Definitions

  • the present invention relates to a method for measuring the speed of a microscopic object in motion within a flow, such as a blood flow, from a light scanning microscope. It applies in particular but not exclusively to the study of microcirculation in which the problem is to detect particles such as red blood cells or leukocytes in movement, to estimate the direction of movement of these particles as well as their speeds.
  • the present invention can however be applied to other fields such as for example microfluidics.
  • the speed of movement of an object from an imaging system is estimated conventionally from a sequence of images representing this moving object. The basic assumption being that the acquisition rate of the imaging system is such that the object makes small displacements from one image to another.
  • the measurement is made on a video sequence, using two slots placed on the screen broadcasting the sequence.
  • the two slits are parallel and oriented perpendicular to the vessel on which the measurement is made.
  • a photodiode is present in front of each of the two slots.
  • the device used offers two measurement modes. In a first mode, the two slits are spaced (in equivalent distance on the fabric) from 45 ⁇ m to 70 ⁇ m.
  • the speed measurement is carried out by a calculation of 1 inter-correlation of the two signals from the photodiodes.
  • the two slits are very close together: 7.4 ⁇ m in equivalent tissue distance, a distance slightly less than the average diameter of a red blood cell.
  • the SLO for "Scanning Laser Ophtal oscope" in English is a device whose principle is based on non-fiber confocal microscopy. Images are captured at a rate of 50 interlaced images per second. The speed measurement method used by this device is based on tracking cell.
  • the same image is made up of interlaced fields corresponding to two instants spaced 20 ms apart, a globule in motion therefore appears in two different places on the image, once on the even lines and another on the odd lines.
  • the speed measurement is immediate.
  • the large field of the device up to 1200 ⁇ m makes it possible to measure speeds of several cm / s).
  • the limitations of this device are mainly due to the monitoring of the globules: marked globules of large size (12 ⁇ m for leukocytes) and in small numbers are required.
  • Doppler effect describes the frequency shift experienced by a wave reflected by a moving object relative to the observer.
  • moving objects are red blood cells.
  • a monofrequential wave for example a laser
  • the speed measurement by Doppler effect has the advantage of being very fast and precise.
  • the maximum measurable speeds vary from 1 mm / s to several cm / s.
  • a drawback of the Doppler comes from the difficulty of precisely identifying the area over which the speed is measured, especially at depth. The reflected light can come from different vessels where blood flows at different speeds.
  • Spatio-Temporal Analysis Such a method for measuring the speed of leukocytes is described in the document "Measuring microcirculation using spa tiotemporal image analysis.” By Yoshinobu Sato et al., CVRMed, pages 302-308, 1995.
  • This analysis comprises three main stages: i) First, the vessel is extracted from the image by a segmentation made on the histogram of the temporal variances of the pixels of the frames of the sequence, ii) A spatiotemporal image is then constructed from successive images of the vessel. This image can be three-dimensional, or two-dimensional if each image of the sequence is projected on the axis or the contours of the vessel, iii) The traces left in the space-time image by the movement of leukocytes are then reinforced by application of a bank of filters with selective orientation (for example a Gabor filter bank). The traces are extracted by thresholding the best responses. The last step is to connect the traces together to reconstruct the entire trajectories of the leukocytes.
  • the present invention aims to remedy the aforementioned drawbacks by proposing a method capable of measuring the speed of red blood cells in particular.
  • the object of the present invention is to measure the speed of rapidly moving particles by means of a scanning imaging system. By rapid movement is meant a speed beyond about 2mm / s.
  • this method comprises the following steps: acquisition of an image by light scanning in x and y of a plane containing said object, this plane is also called subsurface image field located a few ⁇ m in a sample or tissue organic; - Detection in the plane (x, y) of a streak generated by the displacement of said object during the acquisition of said image; determining the slope of said streak in the plane (x, y ); estimation of the speed V g of said object from the slope thus determined.
  • a device of the SLO type or any other scanning image acquisition device (x, y) for which the scanning speeds are adapted to the implementation of the method according to the invention.
  • a microscope, confocal or not, fiber laser, in particular single fiber with distal scanning. This distal scanning can be carried out by a micromirror, by displacement of lenses or optics, by spectral scanning ...
  • a single image acquired by a two-dimensional plane scan can be used.
  • a morphological image is used.
  • the invention includes a streak detection step in which the following steps are carried out: accentuation of a set of streaks in the image by applying a filter; application of a threshold so as to keep the most important streaks; - adjustment of a line or an ellipse on each of these streaks; and identifying said streak.
  • a streak detection step in which the following steps are carried out: accentuation of a set of streaks in the image by applying a filter; application of a threshold so as to keep the most important streaks; - adjustment of a line or an ellipse on each of these streaks; and identifying said streak.
  • a first simplification consists in considering the trajectory of the spot as horizontal on the acquisition window. This simplification is justified by the ratio of one thousand existing between the horizontal speed and the vertical speed of the spot, the calculation confirms that the vertical position of the spot varies by less than 0.1 ⁇ m on the acquisition window.
  • a second approximation is to consider the time necessary for the spot to cover the acquisition window horizontally as negligible, that is to say to consider the objects as stationary during the course of a scanning line.
  • the speed of the objects observed in the vessels considered, the red blood cells have a speed less than 20 mm / sec) compared to the horizontal speed of the spot (> 1 m / s) justifies this approximation.
  • the scanning is modeled by instantaneous horizontal scanning lines spatially spaced by a distance V y // X , with f x scanning frequency "x". Once the scan has been modeled, it remains to model the moving objects, for example in the following red blood cells.
  • red blood cells from the simple rod to a realistic three-dimensional model.
  • the simplest model that can account for the observed deformations consists in representing a red blood cell by a vertical rod.
  • V y tan (1) g cos6 "with V y the vertical speed of the light spot used for scanning, V g the speed of the sought globule , ⁇ the angle between the vector V g and the fast scan x axis.
  • Tan ( ⁇ ) is the slope of the streaks. If the trajectory of the globules is supposed to be collinear with the edges of the vessel carrying the object then it is possible to know ⁇ by detecting the edges of the vessel and therefore to access the speed V g of the globule. In the general case where ⁇ is unknown but where the length L of the streak is known, two possibilities are possible: 1) The size D of the globule is known. In this case, we consider t v the visibility time of the globule, defined as the time between the first intersection and the last intersection of the trajectories. The globule has a vertical speed V g sin ( ⁇ ) and a vertical extent D.
  • the object such as a red blood cell
  • a vertical rod of length D the object, such as a red blood cell
  • the red blood cell can be represented in the form of a solid sphere or of a biconcave shape having a symmetry of revolution.
  • the orthogonal projection in the observation plane is assimilated to a disk of radius R. On the acquired image, this moving disk appears as an ellipse whose angle ⁇ between the streak and the axis " x "is given by:
  • a confocal microscopy system with light scanning in fiber mode used to measure the speed of the microscopic object in motion within a flow, such as a blood flow
  • this system comprising: means for acquiring an image by light scanning in x and y of a plane containing said object; - Means for detecting a streak generated by the displacement of said object during the acquisition of said image; means for determining the slope of said streak; and means for estimating the speed V g of said object from the slope thus determined.
  • FIG. 2 is a very simplified diagram illustrating the scanning mode of the imaging system of Figure 1;
  • Figure 3 is a schematic image representing streaks, this image being from a simulated acquisition; and
  • FIG. 4 is a schematic view of a step of selecting streaks.
  • the arterioles have a diameter which varies between 50 ⁇ m and 100 ⁇ m; the capillaries are much thinner with a diameter of 3 ⁇ m to 8 ⁇ m; finally, the venules have a diameter of 30 ⁇ m to 50 ⁇ m; - the speed of the red blood cells in these vessels is in a range from less than 1 mm / s for the smallest vessels to a few tens of mm / s for the arterioles; - red blood cells are cells with an average diameter of approximately 7 ⁇ m, to be compared with diameters of 10 ⁇ m to 15 ⁇ m for white blood cells.
  • the method described in this document WO2004 / 008952A1 provides for the focusing of the beam at the fiber outlet in order to excite a subsurface plane and produce a confocal image.
  • the method provides for producing a divergent beam at the outlet of the fiber capable of exciting a microvolume of the tissue from the surface.
  • the excitation signal is deflected at a speed corresponding to the acquisition of a sufficient number of images per second for use in real time and the fluorescence signal is detected at a detection frequency corresponding to a minimum frequency fiber sampling one by one.
  • the injection system comprises several optical elements 5 preceded by a scanning of fibers 6, such as a deflector, making it possible to scan the fibers one by one at very high speed. Each fiber is used in turn to convey the illumination beam and also the corresponding return beam from the object observed.
  • the spatial resolution is obtained by focusing the laser beam at a point and by the confocal character residing in the spatial filtering of the object observed by the same fibers as those used for the illumination. This makes it possible to receive, by means of a photodetector 9, exclusively the signal coming from the object observed and to produce an image point by point.
  • the image guide 1 consists of a very large number of flexible optical fibers, for example 30,000 fibers 2 ⁇ m in diameter and spaced 3.3 ⁇ m apart. In practice, it is possible to use either all of the fibers of the image guide, or a selected subset of these fibers, for example centered.
  • the electronic and computer means 7 for controlling, analyzing and digitally processing the detected signal and for viewing include in particular the following cards: - a synchronization card 8 which has the following functions: - to synchronously control the scanning; to know at all times the position of the laser spot thus scanned; and to manage all the other cards via a microcontroller itself which can be controlled; a detector card 9 which includes an analog circuit which in particular performs an impedance matching, an analog digital converter then a programmable logic component (for example an FPGA circuit) which formats the signal; a digital acquisition card 10 which makes it possible to process a stream of digital data at variable frequency and to display it on a screen 11; - a graphics card 12. Alternatively, one can use a single card combining the functionality of these different cards.
  • FIG. 2 is a very simplified diagram illustrating the scanning mode of the imaging system of FIG. 1.
  • the scanning laser spot is symbolized by the dotted lines 13 which describe a conventional scanning trajectory in a square scanning window 14.
  • the trajectory of the laser spot 13 is a "Z" from top to bottom.
  • the horizontal speed V x along the horizontal axis A x is assumed to be much higher compared to the speed V y along the axis A y . This hypothesis amounts to neglecting the time for the globule to catch up with the spot between two horizontal scans.
  • the speed V y can be 3 mm / s, while that of V x can be 5 m / s.
  • the image guide 1 is also shown in a cross-sectional view.
  • the optical fibers are represented in the form of ordered circles.
  • the imaged area 15 only corresponds to a number limited optical fibers located inside a rectangle.
  • the laser beam is successively injected into each of the optical fibers.
  • the image shown in Figure 3 simulates an acquisition, that is to say that each fiber has been injected only once. This image shows oblique streaks corresponding to the interaction between the scanning system and the moving particles.
  • each streak is then framed by an ellipse to define a slope.
  • the electronic and computer means 7 determine the slope of each streak so as to calculate the speed of each red blood cell.
  • the present invention therefore makes it possible to determine the speed from a single image. This makes it possible in particular to avoid problems of shaky image matching. It advantageously makes it possible to apprehend significant speeds with regard to the field of view and the acquisition frequency. For example, an image acquisition system at 11 Hz makes it possible to apprehend speeds of the order of 5 to 25 mm / sec, i.e. speeds that are impossible to estimate with most of the techniques of the art prior.
  • the globule speed is measured, the vertical movement of which does not go against the spot.
  • the angle of the trajectory must be between the horizontal and the critical angle for which the vertical speed of the globule becomes equal to that of the spot.
  • the invention is not limited to the examples which have just been described and numerous modifications can be made to these examples without departing from the scope of the invention.

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Ophthalmology & Optometry (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Dispersion Chemistry (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Optics & Photonics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Fluid Mechanics (AREA)
  • Vascular Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Hematology (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Microscoopes, Condenser (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
  • Measuring Volume Flow (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

L'invention concerne un procédé pour mesurer la vitesse d'une particule telle qu'un globule rouge en mouvement au sein d'un écoulement, tel qu'un flux sanguin, à partir d'un microscope à balayage lumineux, ce procédé comprenant les étapes suivantes: acquisition d'une image par balayage lumineux en x et y d'un plan contenant ledit objet; détection dans le plan (x, y) d'une strie engendrée par le déplacement dudit objet lors de l'acquisition de ladite image ; détermination de la pente de ladite strie dans le plan (x, y) ; estimation de la vitesse vg dudit objet à partir de la pente ainsi déterminée.

Description

" Procédé et système de mesure de vitesse du flux sanguin "
La présente invention se rapporte à un procédé pour mesurer la vitesse d'un objet microscopique en mouvement au sein d'un écoulement, tel qu'un flux sanguin, à partir d'un microscope à balayage lumineux. Elle s'applique notamment mais non exclusivement à l'étude de la microcirculation dans lequel la problématique est de détecter des particules telles que des globules rouges ou des leucocytes en déplacement, d'estimer le sens de déplacement de ces particules ainsi que leurs vitesses. La présente invention peut toutefois s'appliquer à d'autres domaines tels que par exemple la microfluidique. L'estimation de la vitesse de déplacement d'un objet à partir d'un système d'imagerie se fait classiquement à partir d'une séquence d'images représentant cet objet en mouvement. L'hypothèse de base étant que la cadence d'acquisition du système d'imagerie est telle que l'objet effectue de petits déplacements d'une image à l'autre. La présence de l'objet, à différentes positions, sur la série temporelle d'images permet alors de remonter, moyennant un étalonnage du système d'acquisition, à la vitesse de ce dernier. On connaît le document "Erythrocyte velocity measurement in microvessels by a two-slit photometric method" de H. ayland et RC. Johnson, publié dans J. Appl . Physiol. 22(2) :333-337, 1967, et décrivant une méthode photométrique à deux fentes . La mesure par méthode photométrique à deux fentes est probablement la technique la plus ancienne de mesure automatisée de vitesse du flux sanguin. Cette méthode mesure la vitesse des globules rouges et s'applique de préférence aux capillaires et aux petites veinules dans lesquelles les globules circulent en file, isolés ou par petits agrégats. La mesure se fait sur une séquence vidéo, à l'aide de deux fentes placées sur l'écran diffusant la séquence. Les deux, fentes sont parallèles et orientées perpendiculairement au vaisseau sur lequel est réalisée la mesure. Une photodiode est présente en face de chacune des deux fentes. L'appareil utilisé propose deux modes de mesures. Dans un premier mode, les deux fentes sont espacées (en distance équivalente sur le tissu) de 45 μm à 70 μm. La mesure de vitesse est réalisée par un calcul de 1' inter-corrélation des deux signaux issus des photodiodes. Dans le second mode de mesure, les deux fentes sont très rapprochées : 7.4 μm en distance équivalente tissu, soit une distance légèrement inférieure au diamètre moyen d'un globule rouge. En conséquence, deux signaux consécutifs générés respectivement par la diode amont et la diode aval sont causés par un même globule, et il est donc possible de calculer la vitesse de ce globule. Cependant, les limitations de la méthode sont la nécessité de faire la mesure sur des vaisseaux très fins afin de limiter l'observation à un globule rouge, et les vitesses mesurées, qui en raison de la fréquence image de 30 images par seconde, ne peuvent dépasser 2 mm/s. On connaît également le méthode de la projection Spatio- Temporelle ("Line shift Diagram"). Cette méthode est une extension de la précédente. Au lieu d'échantillonner le signal en deux points, l'utilisateur sélectionne une zone d'intérêt, c'est à dire un rectangle inscrit dans un vaisseau. A chaque image, une moyenne des niveaux de gris est calculée sur la largeur du vaisseau en chaque point de l'axe du vaisseau, le signal de la zone d'intérêt est projeté en une dimension sur l'axe du vaisseau. Ensuite les signaux unidimensionnels obtenus pour chaque image sont alignés verticalement pour donner une image spatiotemporelle qui laisse apparaître les traces des globules. La vitesse est évaluée en corrélant les signaux adjacents. Cette méthode est utilisée par les logiciels Caplmage® et Capiscope®, de pair avec l'appareil d'acquisition Cytoscan®. Une telle méthode est notamment décrite dans le document "Orthogonal Polarisation Spectral Imagina : A new method for study of microcircula tion . " de . Groner, I. . inkelman, A. G. Harris, G. Inde, G.I. Bouma, K. Messmer, et R.G. Nadeau, publié dans Nature Médecine, 5 : 1209-1213, 1999. La principale limitation est la gamme de vitesse mesurée, qui ne peut dépasser 2 mm/s en raison de la fréquence des images (25 images/secondes ou 50 images/secondes si les mesures sont réalisées alternativement sur les deux champs entrelacés du flux vidéo) . Le SLO pour "Scanning Laser Ophtal oscope" en langue anglaise, est un appareil dont le principe repose sur la microscopie confocale non fibrée . Les images sont capturées à une cadence de 50 images entrelacées par seconde. La méthode de mesure de vitesse utilisée par cet appareil repose sur le suivi de cellule. Une même image est constituée de champs entrelacés correspondants à deux instants espacés de 20 ms, un globule en mouvement apparaît donc à deux endroits différents sur l'image, une fois sur les lignes paires et une autre sur les lignes impaires. Une fois les deux images du globule localisées, la mesure de la vitesse est immédiate. Le grand champ de l'appareil (jusqu'à 1200 μm permet de mesurer des vitesses de plusieurs cm/s) . Les limitations de cet appareil sont principalement dues au suivi des globules : il faut des globules marqués de grosse taille (12 μm pour les leucocytes) et en faible nombre. Une telle méthode semble inapplicable à la mesure de vitesse de globules rouges, dont les dimensions sont plus petites (7 - 8 μm) , de concentration beaucoup plus forte (1 000 fois plus importante que celle des globules blancs) et plus difficiles à marquer. Une autre méthode de mesure qui est couplée couramment avec le Cytoscan® ou le SLO est la mesure de vitesse de globules rouges par effet Doppler. L'effet Doppler décrit le décalage fréquentiel que subit une onde réfléchie par un objet en mouvement par rapport a l'observateur. Dans le cas du flux sanguin, les objets en mouvements sont les globules rouges. Une onde monofrequentielle (par exemple un laser) de longueur d'onde donnée est envoyée sur un vaisseau sanguin. La mesure de la vitesse par effet Doppler a l'avantage d'être très rapide et précise. Suivant le matériel et les logiciels d'analyse utilisés les vitesse maximales mesurables varient de 1 mm/s à plusieurs cm/s. Un inconvénient du Doppler vient de la difficulté à identifier précisément la zone sur laquelle la vitesse est mesurée, surtout en profondeur. La lumière réfléchie peut provenir de différents vaisseaux où le sang circule à des vitesses différentes. On connaît encore la méthode par Analyse Spatio-Temporelle. Une telle méthode pour la mesure de la vitesse de leucocytes est décrite dans le document "Measuring microcircula tion using spa tiotemporal image analysis . " par Yoshinobu Sato et al., CVRMed, pages 302-308, 1995. Cette analyse comporte trois étapes principales : i) Dans un premier temps, le vaisseau est extrait de l'image par une segmentation faite sur l'histogramme des variances temporelles des pixels des trames de la séquence, ii) Une image spatiotemporelle est alors construite à partir des images successives du vaisseau. Cette image peut être tridimensionnelle, ou bidimensionnelle si chaque image de la séquence est projetée sur l'axe ou les contours du vaisseau, iii) Les traces laissées dans l'image spatiotemporelle par le mouvement des leucocytes sont ensuite renforcées par application d'un banc de filtres à orientation sélective (par exemple un banc de filtres de Gabor) . Les traces sont extraites par seuillage des meilleures réponses. La dernière étape consiste à connecter entre elles les traces pour reconstituer les trajectoires entières des leucocytes. L'obtention de ces traces permet ensuite par un calcul de leurs tangentes d'avoir une estimation de la vitesse des leucocytes sur leurs trajectoires . Egalement dans le domaine de l'analyse spatio-temporelle, on connaît les documents : "Two-photon imaging of neocortical microcirculation . " D. Kleinfeld and W. Denk; In Imaging Neurons: A Laboratory Manual (R. Yuste, F. Lanni, and A. Konnerth, editors) , 1999, Cold Spring Harbor Laboratory Press, NY, pp. 23.1-23.15; accessible sur Internet à l'adresse suivante : http: //physics .ucsd.edu/neurophysics/publications/kleinfeld_denk_c shl_2003.pdf; et "Two-photon imaging of capillary blood flow in olfactory bulb glomeruli " E. Chaigneau et al., PNAS, 28 Octobre, 2003; 100(22): 13081 - 13086; accessible sur Internet à l'adresse suivante : http://www.pnas.org/cgi/content/full/100/22/13081. Ces documents décrivent une méthode de mesure de vitesse en réalisant, à l'aide d'un microscope non fibre, plusieurs balayages successifs sur l'axe du vaisseau sanguin. On obtient donc plusieurs images "unidimensionnelles" d'un même segment à des instants différents successifs. Ces images sont mises bout à bout de façon à former une visualisation globale dont l'ordonnée est un repère temporel. Le mouvement des particules dans le vaisseau sanguin se traduit sous forme de bandes obliques. On détermine la vitesse instantanée en calculant la pente de chaque bande sur la visualisation globale. Toutefois, l'inconvénient de cette méthode réside sur le fait qu'il est indispensable de positionner le dispositif d'acquisition parallèlement à l'axe du vaisseau sanguin. Par ailleurs, la constitution de la visualisation globale nécessite l'acquisition de plusieurs images d'un même segment. Dans la plupart des techniques qui viennent d'être présentées, la mesure de la vitesse repose sur une analyse entre au moins deux images successives d'une acquisition; ce qui par ailleurs soulève des problèmes de bougé et de mise en correspondance entre images. Dès lors, la gamme de vitesse mesurable par ces techniques dépend du champ de vue ainsi que de la fréquence image. Or, pour des systèmes d'imagerie à balayage présentant un petit champ de vue, par exemple de l'ordre de 166 μm x 118 μm, et des cadences d'acquisition, par exemple de l'ordre de 12 images par seconde, une particule animée d'une vitesse supérieure à 1.8 mm/s traversera le champ d'observation entre deux images successives, ce qui rend impossible l'utilisation de la plupart des techniques citées précédemment. La présente invention vise à remédier aux inconvénients précités en proposant un procédé apte à mesurer la vitesse des globules rouges notamment. La présente invention a pour but la mesure de vitesse de particules en mouvement rapide au moyen d'un système d'imagerie à balayage. Par mouvement rapide, on entend une vitesse au delà d'environ 2mm/s. On atteint le but recherché avec un procédé pour mesurer la vitesse d'un objet microscopique, tel que par exemple des globules rouges et des leucocytes en mouvement au sein d'un écoulement, tel qu'un flux sanguin, à partir d'un microscope à balayage lumineux. Selon l'invention, ce procédé comprend les étapes suivantes : acquisition d'une image par balayage lumineux en x et y d'un plan contenant ledit objet, ce plan est aussi appelé champ imagé subsurfacique situé à quelques μm dans un échantillon ou un tissu biologique; - détection dans le plan (x, y) d'une strie engendrée par le déplacement dudit objet lors de l'acquisition de ladite image; détermination de la pente de ladite strie dans le plan (x, y); estimation de la vitesse Vg dudit objet à partir de la pente ainsi déterminée. On peut utiliser un microscope, confocal ou non, à balayage lumineux, notamment laser, en mode fibre ou non fibre. Pour mettre en œuvre le mode de réalisation non fibre, on peut utiliser, en adaptant l'électronique de traitement, un appareil de type SLO ou tout autre appareil d'acquisition d'image à balayage (x, y) pour lequel les vitesses de balayage sont adaptées à la mise en œuvre du procédé selon l'invention. En complément de ce qui précède, on peut utiliser de façon non limitative un microscope, confocal ou non, laser fibre, notamment mono-fibre, à balayage distal. Ce balayage distal peut être réalisé par un micromiroir, par déplacement de lentilles ou d'optiques, par balayage spectral... On peut aussi utiliser un système mono-fibré à balayage proximal où le balayage lumineux est obtenu par déviation d'une extrémité de la fibre optique au sein d'une tête optique proche de l'objet observé. Contrairement aux méthodes d'analyse spatiotemporelle de
Kleinfeld et de Chaigneau où l'on acquiert une succession d'images "unidimensionnelles" d'un même segment dans l'axe du vaisseau sanguin, dans la présente invention, on peut utiliser une seule image acquise par un balayage plan bidimensionnel . Contrairement à l'art antérieur où l'on mesure des stries sur une image spatio-temporelle (t=f(x)), ici les stries sont complètement différentes puisqu'elles sont issues d'une image plan (y=f(x)), la notion temporelle étant induite dans le système de balayage point à point en "Z". On se sert avantageusement d'une image morphologique. Par exemple, dans le document de Yoshinobu Sato selon l'art antérieur, pour obtenir des stries, on réalise plusieurs images dans lesquelles on extrait une seule ligne dans chaque image, et on élabore ainsi l'image finale t=f(x). L'utilisation d'un système d'imagerie à balayage, système dans lequel l'objet en déplacement va être observé pendant le temps de génération d'une image (temps de parcours trame), conduit à l'apparition de stries permettant de remonter à une estimation de la vitesse d'objets en déplacement à partir d'une seule image. Ces stries proviennent donc de l'interaction entre l'objet en mouvement et le système de balayage. L'invention est notamment remarquable par le fait que ces stries sont généralement considérées par l'homme du métier comme des parasites à éliminer dans un système d'imagerie à balayage. On tire donc profit des éléments considérés comme des aberrations dans les microscopes à balayage laser, notamment en mode fibre. Pour ce faire, l'invention comporte une étape de détection de stries dans laquelle on réalise les étapes suivantes : accentuation d'un ensemble de stries de l'image par application d'un filtre; application d'un seuil de façon à conserver les stries les plus importantes; - ajustement d'une droite ou d'une ellipse sur chacune de ces stries; et identification de ladite strie. Selon l'invention, compte tenu des ordres de grandeurs des vitesses de balayage et de celles des objets observés, il est possible de simplifier le modèle de la trajectoire du spot lumineux. Une première simplification consiste à considérer la trajectoire du spot comme horizontale sur la fenêtre d'acquisition. Cette simplification est justifiée par le rapport mille existant entre la vitesse horizontale et la vitesse verticale du spot, le calcul confirme que la position verticale du spot varie de moins de O.lμm sur la fenêtre d'acquisition. Une seconde approximation est de considérer le temps nécessaire au spot pour couvrir horizontalement la fenêtre d'acquisition comme négligeable, c'est à dire de considérer les objets comme immobiles pendant le parcours d'une ligne de balayage. La vitesse des objets observés (dans les vaisseaux considérés, les globules rouges ont une vitesse inférieure à 20 mm/sec) par rapport à la vitesse horizontale du spot (> 1 m/s) justifie cette approximation. Au final le balayage est modélisé par des lignes de balayage horizontal instantané espacées spatialement d'une distance Vy//X, avec fx fréquence du balayage "x". Une fois le balayage modélisé, il reste à modéliser les objets en mouvement, par exemple dans la suite les globules rouges. Plusieurs modèles ont été envisagés pour décrire les globules rouges : du simple bâtonnet à un modèle tridimensionnel réaliste. Le modèle le plus simple pouvant rendre compte des déformations observées consiste à représenter un globule rouge par un bâtonnet vertical. Dans le cas de l'utilisation de ce modèle, on obtient une relation simple liant l'angle α des stries observées (par rapport à l'axe "x") à la vitesse verticale du balayage Vy et à la vitesse de déplacement horizontal Vgcos (θ) du globule : Vy tan = (1) gcos6» avec Vy la vitesse verticale du spot lumineux utilisé pour le balayage, Vg la vitesse du globule recherché, θ l'angle entre le vecteur Vg et l'axe x de balayage rapide. Tan(α) est la pente des stries . Si la trajectoire des globules est supposée colinéaire au bords du vaisseau véhiculant l'objet alors il est possible de connaître θ en détectant les bords du vaisseau et donc d'accéder à la vitesse Vg du globule. Dans le cas général où θ est inconnu mais où la longueur L de la strie est connue, deux possibilités sont envisageables : 1) La taille D du globule est connue. Dans ce cas, on considère tv le temps de visibilité du globule, défini comme le temps écoulé entre la première intersection et la dernière intersection des trajectoires. Le globule a une vitesse verticale Vgsin(θ) et une étendue verticale D. On assimile donc l'objet à une barre perpendiculaire à la direction d'écoulement, de longueur finie suffisamment grande pour s'étaler sur plusieurs lignes de balayage horizontal. Le spot est ponctuel et a une vitesse verticale Vy, tv dépend donc uniquement de la différence de vitesse verticale du spot et du globule, et de D selon la relation :
D (2) -Fgsin0
Les étendues verticale L|sin(α)| et horizontale L|cos(α)| de la strie observée sont directement reliées à tv et aux vitesses verticale du spot et horizontale du globule par :
Z,|sinα| = tvJ | ( 3 ) Les relations 2 et 3 permettent d'écrire
Ce qui, combiné à la relation (1), donne un système de deux équations à deux inconnues θ et Vy. La résolution doit se faire selon le signe de Vy-Vgsin(θ), c'est à dire selon que le spot se déplace plus ou moins vite verticalement que le globule. Dans le cas où les deux vitesses sont égales, le temps de visibilité devient infini et la strie devient infiniment longue. On a alors Vgsin(θ) = Vy , en remplaçant dans la relation (1), il vient tan(α) = tan(θ). Par conséquent, la strie est exactement la trajectoire du globule et la vitesse Vg est donnée par :
[5)
Si le spot lumineux se déplace verticalement plus vite que le globule (Vy>Vgsin (θ) ) , alors la résolution du système (3) donne :
D tan# = tanα (6) Lcosa
Si le spot lumineux se déplace verticalement moins vite que le globule (Vy<Vgsin (θ) ) , alors la résolution du système (3) donne
D tan< = tanα +• (7) Lcosa
2) Lorsqu'un balayage inversé est disponible Si l'on dispose d'un autre balayage, inversé par rapport au précédent, c'est à dire de vitesse verticale -Vy, une strie d'angle -α et de longueur L' sera observée. L' est donné par : (8)
Le rapport des relations (4) et (8) permet d'écrire :
deux cas sont alors à distinguer : | Vy| > | Vgsin (θ) | et |vy|<|vgsin(θ) | (si | Vy | = | Vgsin (θ) | , alors la relation (5) s 'applique) . Si | Vy I > | Vgsin (θ) | alors la résolution donne : tan θ = tanα L ~ L' L + L' (10)
VB Si I Vv|< | Vgsin(θ) | alors la résolution donne
L + L' tan# = tanα L -L' (11)
Dans tout ce qui précède, l'objet, tel qu'un globule rouge, est assimilé à un bâtonnet vertical de longueur D. L'homme du métier comprendra aisément qu'on peut remplacer cette hypothèse par celle, plus réaliste, d'une sphère (ou disque), voire d'une ellipse ou d'une structure plus compliquée encore et évaluer l'écart entre l'état observé dans ce cas et celui observé dans le cas simple. Selon la précision recherchée, il conviendra de prendre ou non en compte cet écart. A titre d'exemple, le globule rouge peut être représenté sous la forme d'une sphère solide ou d'une forme biconcave possédant une symétrie de révolution. Dans les deux cas, on assimile la projection orthogonale dans le plan d'observation à un disque de rayon R. Sur l'image acquise, ce disque en mouvement apparaît comme une ellipse dont l'angle α entre la strie et l'axe "x" est donné par :
La longueur de l'axe principal est donnée par V avec Vr = — et R le rayon du disque considéré. y
Des simulations numériques effectuées avec des modèles plus complexes de représentation des globules rouges ( surfaces toroidales) montrent que le modèle du disque est suffisant pour représenter les globules rouges. Si l'orientation θ de la trajectoire du globule est inconnue, mais la longueur L des stries est disponible, il est possible de remonter partiellement à une information de vitesse. Les relations (12) et (13) forment un système liant le couple (α, L) qui sont les paramètres observables lors des acquisitions au couple (θ, Vg) qui sont les paramètres recherchés. Selon un autre aspect de l'invention, il est proposé un système de microscopie confocale à balayage lumineux en mode fibre, utilisé pour mesurer la vitesse de l'objet microscopique en mouvement au sein d'un écoulement, tel qu'un flux sanguin, ce système comprenant : des moyens pour acquérir une image par balayage lumineux en x et y d'un plan contenant ledit objet; - des moyens pour détecter une strie engendrée par le déplacement dudit objet lors de l'acquisition de ladite image; des moyens pour déterminer la pente de ladite strie; et des moyens pour estimer la vitesse Vg dudit objet à partir de la pente ainsi déterminée. D'autres avantages et caractéristiques de l'invention apparaîtront à l'examen de la description détaillée d'un mode de mise en œuvre nullement limitatif, et des dessins annexés, sur lesquels : la figure 1 est un schéma général d'un exemple de système d'imagerie confocale fibrée mettant en œuvre le procédé selon
1 ' invention; - la figure 2 est un schéma très simplifié illustrant le mode de balayage du système d'imagerie de la figure 1; la figure 3 est une image schématique représentant des stries, cette image étant issue d'une acquisition simulée; et la figure 4 est une vue schématique d'une étape de sélection de stries. Bien que l'invention n'y soit pas limitée, on va maintenant décrire le procédé selon l'invention mis en œuvre dans un microscope confocal à balayage laser en mode fibre, ce procédé s 'appliquant au domaine de la microcirculation, dont les ordres de grandeurs sont les suivants : - les artérioles ont un diamètre qui varie entre 50 μm et 100 μm; les capillaires sont beaucoup plus fins avec un diamètre de 3 μm à 8 μm; enfin les veinules ont un diamètre de 30 μm à 50 μm; - la vitesse des globules rouges dans ces vaisseaux est comprise dans une gamme allant de moins de 1 mm/s pour les vaisseaux les plus petits à quelques dizaines de mm/s pour les artérioles; - les globules rouges sont des cellules dont le diamètre moyen est d'environ 7 um, à comparer avec des diamètres de 10 μm à 15 μm pour les globules blancs. D'une manière générale, pour mettre en œuvre la présente invention, on peut se baser sur le système décrit dans le document W0 2004/008952A1, « Procédé et appareillage d'imagerie de fluorescence haute résolution par fibre optique et notamment d'imagerie confocale », Mauna Kea Technologies, dans lequel on utilise un guide d' image fait de plusieurs milliers de fibres optiques, un signal d'excitation étant émis par une source, dévié et injecté tour à tour dans l'une des fibres dudit guide, chaque point d'excitation du tissu en sortie de fibre émettant en retour un signal de fluorescence collecté par ladite fibre, puis détecté et numérisé pour former un élément d'image. Selon un premier aspect, le procédé décrit dans ce document W02004/008952A1 prévoit la focalisation du faisceau en sortie de fibre pour exciter un plan subsurfacique et réaliser une image confocale. Selon un second aspect, le procédé prévoit de produire un faisceau divergent en sortie de fibre susceptible d'exciter un microvolume du tissu depuis la surface. Le signal d'excitation est dévié à une vitesse correspondant à l'acquisition d'un nombre d'images par seconde suffisant pour une utilisation en temps réel et l'on détecte le signal de fluorescence à une fréquence de détection correspondant à une fréquence minimale d'échantillonnage des fibres une à une. Sur la figure 1 on voit un faisceau ordonné de fibres optiques souples formant un guide d'image 1 avec, sur son extrémité proximale, une source lumineuse 2 et un système d'injection de fibres permettant d'illuminer les fibres une à une et, sur son extrémité distale, une tête optique 3 permettant de focaliser le faisceau sortant de la fibre illuminée en un point situé à une profondeur donnée de l'objet observé 4. Le système d'injection comprend plusieurs éléments optiques 5 précédés d'un système de balayage de fibres 6, tel qu'un déviateur, permettant de balayer les fibres une à une à très grande vitesse. Chaque fibre est utilisée tour à tour pour véhiculer le faisceau d'illumination et également le faisceau de retour correspondant provenant de l'objet observé. La résolution spatiale est obtenue par focalisation du faisceau laser en un point et par le caractère confocal résidant dans le filtrage spatial de l'objet observé par les mêmes fibres que celles ayant servi à l'illumination. Cela permet de réceptionner, au moyen d'un photodétecteur 9, exclusivement le signal provenant de l'objet observé et de réaliser une image point par point. Le guide d'image 1 est constitué d'un très grand nombre de fibres optiques souples, par exemple 30 000 fibres de 2 μm de diamètre et espacées de 3,3 μm. En pratique, on peut utiliser soit l'ensemble des fibres du guide d'image, soit un sous-ensemble choisi de ces fibres, par exemple centré. Les moyens électroniques et informatiques 7 de commande, d'analyse et de traitement numérique du signal détecté et de visualisation comprennent notamment les cartes suivantes: - une carte de synchronisation 8 qui a pour fonctions : - de commander de manière synchronisée le balayage; de connaître à tout instant la position du spot laser ainsi balayé ; et de gérer toutes les autres cartes par l'intermédiaire d'un microcontrôleur lui-même pouvant être piloté ; - une carte détecteur 9 qui comprend un circuit analogique qui réalise notamment une adaptation d'impédance, un convertisseur analogique numérique puis un composant logique programmable (par exemple un circuit FPGA) qui met en forme le signal ; - une carte d'acquisition numérique 10 qui permet de traiter un flot de données numériques à fréquence variable et de l'afficher sur un écran 11; - une carte graphique 12. En variante, on peut utiliser une seule carte regroupant les fonctionnalités de ces différentes cartes. Ces moyens électroniques et informatiques 7, aptes à réaliser les étapes du procédé selon l'invention, peuvent se présenter sous forme d'un micro ordinateur doté de moyen de traitement nécessaire pour calculer la vitesse des globules rouges. La figure 2 est un schéma très simplifié illustrant le mode de balayage du système d'imagerie de la figure 1. Le spot laser de balayage est symbolisé par les pointillés 13 qui décrivent un trajectoire de balayage classique dans une fenêtre carrée de balayage 14. La trajectoire du spot laser 13 est un "Z" de haut en bas. La vitesse horizontale Vx selon l'axe horizontale Ax est supposée très supérieure devant la vitesse Vy selon l'axe Ay. Cette hypothèse revient à négliger le temps de rattrapage du globule par le spot entre deux balayage horizontaux. C'est aussi confondre l'effet observé sur les lignes paires avec celui observé sur les lignes impaires. A titre d'exemple, la vitesse Vy peut être de 3 mm/s, alors que celle de Vx peut être de 5 m/s. Sur la figure 2, à l'intérieur de la fenêtre de balayage 14, est également représenté le guide d'image 1 selon une vue en coupe transversale. Les fibres optiques sont représentées sous forme de cercles ordonnés. La zone imagée 15 ne correspond qu'à un nombre limité de fibres optiques situées à l'intérieur d'un rectangle. Le faisceau laser est successivement injecté dans chacune des fibres optiques. L'image représentée sur la figure 3, simule une acquisition, c'est à dire que chaque fibre n'a été injectée qu'une seule fois. Cette image fait apparaître des stries obliques correspondant à l'interaction entre le système de balayage et les particules en mouvement. L'apparition de ces stries s'explique par l'interaction entre l'image des globules rouges et le mécanisme de formation de l'image. Le spot laser effectue le balayage selon une trajectoire en forme de Z, une mesure est effectuée en un ensemble de positions du laser, par exemple 896 mesures par lignes sur 640 lignes. Un globule rouge en mouvement va être intersecté à une position donnée sur une ligne du balayage. A la ligne suivante, ce globule est encore intersecté, il s'est cependant déplacé par rapport à la ligne précédente. Ce phénomène continue tant que l'intersection entre la ligne de balayage et le globule existe. Ce phénomène créé alors une strie dont l'inclinaison est fonction de la vitesse du globule rouge. Sur la figure 4, après acquisition d'une image, les stries sont mises en évidence. Un seuillage permet de ne conserver que les stries les plus importantes. Chaque strie est ensuite encadrée par une ellipse permettant de définir une pente. Puis, les moyens électroniques et informatiques 7 déterminent la pente de chaque strie de façon à calculer la vitesse de chaque globule rouge. La présente invention permet donc de déterminer la vitesse à partir d'une seule image. Ceci permet d'éviter notamment des problèmes de bougés de mise en correspondance d'images. Elle permet avantageusement d'appréhender des vitesses importantes en regard du champ de vue et de la fréquence d'acquisition. A titre d'exemple, un système d'acquisition d'image à 11 Hz permet d'appréhender des vitesses de l'ordre de 5 à 25 mm/sec, soit des vitesses impossibles à estimer avec la plupart des techniques de l'art antérieur. De préférence, on mesure la vitesse de globule dont le mouvement vertical ne va pas à l' encontre du spot. De manière général l'angle de la trajectoire doit se situer entre l'horizontale et l'angle critique pour lequel la vitesse verticale du globule devient égale à celle du spot. Bien sûr, l'invention n'est pas limitée aux exemples qui viennent d'être décrits et de nombreux aménagements peuvent être apportés à ces exemples sans sortir du cadre de l'invention. On peut notamment envisager des mesures faites sur plusieurs stries, voire sur plusieurs images afin d'améliorer la précision.

Claims

REVENDICATIONS
1. Procédé pour mesurer la vitesse d'un objet microscopique en mouvement au sein d'un écoulement, tel qu'un flux sanguin, à partir d'un microscope à balayage lumineux, ce procédé comprenant les étapes suivantes : acquisition d'une image par balayage lumineux en x et y d'un plan contenant ledit objet; détection dans le plan (x, y) d'une strie engendrée par le déplacement dudit objet lors de l'acquisition de ladite image; détermination de la pente de ladite strie dans le plan (x, y) ; et estimation de la vitesse Vg dudit objet à partir de la pente ainsi déterminée.
2. Procédé selon la revendication 1, caractérisé en ce que l'étape de la détection de la strie comprend les étapes suivantes accentuation d'un ensemble de stries de l'image par application d'un filtre; application d'un seuil de façon à conserver les stries les plus importantes; ajustement d'une droite ou d'une ellipse sur chacune de ces stries; et - identification de ladite strie.
3. Procédé selon la revendication 1 ou 2, caractérisé en ce que la vitesse Vg dudit objet est donnée par la relation suivante Vg*cos (θ) = Vy/tan (α) avec Vy la vitesse verticale du spot lumineux utilisé pour le balayage, et "α" l'angle entre l'axe horizontal "x" et la strie.
4. Procédé selon la revendication 3, caractérisé en ce que l'angle θ est obtenu par détection des bords de vaisseaux sanguins véhiculant l'objet.
5. Procédé selon la revendication 3, caractérisé en ce que pour obtenir l'angle θ, on assimile ledit objet à un bâtonnet vertical de hauteur finie D et on calcule l'angle θ à partir de la relation suivante :
où L est la longueur de la strie ,
6. Procédé selon la revendication 5, caractérisé en ce que dans le cas où Vg*sin (θ) < Vy :
D tan# = tanα -- Lcosa
7 . Procédé selon la revendication 5, caractérisé en ce que dans le cas où Vg*sin (θ) > Vy :
D tan# = tanα +- Lcosa v.-W Lsinαry + tan α
8. Procédé selon la revendication 3, caractérisé en ce qu'on acquiert une seconde image du même plan mais dans un balayage inversé, et on utilise la relation suivante : où L' est la longueur de la strie lors du balayage inversé.
9. Procédé selon la revendication 8 et 5, caractérisé en ce que lorsque |vg*sin(θ)| < | Vy | : L - V tan θ - tanα L + L' où L est la longueur de la strie dans un premier sens de balayage, et L' la longueur de la strie lors du balayage inversé.
10. Procédé selon la revendication 5 et caractérisé en ce que lorsque |vg*sin(θ)| > | Vy | : L + L' tan# = tanα- L -L'
11. Procédé selon la revendication 3, caractérisé en ce que lorsque Vg*sin(θ) = Vy , pour déterminer Vg et θ on utilise en outre la relation suivante :
12. Procédé selon la revendication 3, caractérisé en ce que ledit objet étant un globule rouge, on assimile sa forme sur l'image acquise à une ellipse de rayon R dont l'angle α entre la est donné par : 2cos(0) tan(2α) = -2sin(0) et la longueur de l'axe principale est donnée par ΛJVT 2 -2Vr sin(θ) + 2 -|Vr|Λ/vr 2 - 4Vr sin(θ)+ avec v Vτ = — V* vv
13. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on utilise un microscope confocal .
14. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on utilise un microscope à balayage lumineux en mode fibre.
15. Procédé selon l'une quelconque des revendications 1 à 13, caractérisé en ce qu'on utilise un microscope à balayage lumineux non fibre.
16. Système de microscopie à balayage lumineux, utilisé pour mesurer la vitesse d'un objet microscopique en mouvement au sein d'un écoulement, tel qu'un flux sanguin, ce système mettant en œuvre un procédé selon l'une quelconque des revendications précédentes; ce système comprenant : - des moyens pour acquérir une image par balayage lumineux en x et y d'un plan contenant ledit objet; des moyens pour détecter dans le plan (x, y) une strie engendrée par le déplacement dudit objet lors de l'acquisition de ladite image; - des moyens pour déterminer dans le plan (x, y) la pente de ladite strie; et des moyens pour estimer la vitesse Vg dudit objet à partir de la pente ainsi déterminée.
17. Système selon la revendication 16, caractérisé en ce qu'au cours de la détection de la strie, le système comprend : des moyens pour accentuer un ensemble de stries de l'image par application d'un filtre; des moyens pour appliquer un seuil de façon à conserver les stries les plus importantes; des moyens pour ajuster une ellipse sur chacune de ces stries; et des moyens pour identifier ladite strie.
18. Système selon la revendications 16 ou 17, caractérisé en ce qu'on utilise un microscope confocal.
19. Système selon l'une quelconque des revendications 16 à
18, caractérisé en ce qu'on utilise un microscope à balayage lumineux en mode fibre.
20. Système selon l'une quelconque des revendications 16 à 18, caractérisé en ce qu'on utilise un microscope à balayage lumineux non fibre.
EP05744612A 2004-04-02 2005-03-25 Procede et systeme de mesure de vitesse du flux sanguin Active EP1740974B1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0403519A FR2868279B1 (fr) 2004-04-02 2004-04-02 Procede et systeme de mesure de vitesse du flux sanguin
PCT/FR2005/000718 WO2005098474A1 (fr) 2004-04-02 2005-03-25 Procede et systeme de mesure de vitesse du flux sanguin

Publications (2)

Publication Number Publication Date
EP1740974A1 true EP1740974A1 (fr) 2007-01-10
EP1740974B1 EP1740974B1 (fr) 2010-09-01

Family

ID=34944528

Family Applications (1)

Application Number Title Priority Date Filing Date
EP05744612A Active EP1740974B1 (fr) 2004-04-02 2005-03-25 Procede et systeme de mesure de vitesse du flux sanguin

Country Status (12)

Country Link
US (1) US7911590B2 (fr)
EP (1) EP1740974B1 (fr)
JP (1) JP4758418B2 (fr)
CN (1) CN100592107C (fr)
AT (1) ATE479909T1 (fr)
AU (1) AU2005230544B2 (fr)
CA (1) CA2561770C (fr)
DE (1) DE602005023292D1 (fr)
ES (1) ES2352559T3 (fr)
FR (1) FR2868279B1 (fr)
IL (1) IL178347A (fr)
WO (1) WO2005098474A1 (fr)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7414729B2 (en) * 2005-10-13 2008-08-19 President And Fellows Of Harvard College System and method for coherent anti-Stokes Raman scattering endoscopy
CN101520888B (zh) * 2008-02-27 2012-06-27 中国科学院自动化研究所 基于方向场的视网膜图像血管的增强方法
US8144958B2 (en) 2008-09-11 2012-03-27 Carl Zeiss Meditec Ag Medical systems and methods
US8600134B2 (en) 2009-01-30 2013-12-03 Mauna Kea Technologies Method and system for processing images acquired in real time through a medical device
JP5582520B2 (ja) * 2009-02-18 2014-09-03 国立大学法人名古屋大学 血流観測装置
JP5721412B2 (ja) * 2010-12-02 2015-05-20 キヤノン株式会社 眼科装置、血流速算出方法およびプログラム
JP5721411B2 (ja) * 2010-12-02 2015-05-20 キヤノン株式会社 眼科装置、血流速算出方法およびプログラム
JP5792967B2 (ja) * 2011-02-25 2015-10-14 キヤノン株式会社 画像処理装置及び画像処理システム
US9986909B2 (en) 2011-02-25 2018-06-05 Canon Kabushiki Kaisha Image processing apparatus and image processing system for displaying information about ocular blood flow
US9292936B2 (en) * 2013-01-09 2016-03-22 Omiimii Ltd. Method and apparatus for determining location
US9844320B2 (en) * 2014-01-29 2017-12-19 University Of Rochester System and method for observing an object in a blood vessel
DE102017201252A1 (de) * 2017-01-26 2018-07-26 Universität Ulm Verfahren und Vorrichtung zur Untersuchung von Zellen
GB201707239D0 (en) 2017-05-05 2017-06-21 Univ Edinburgh Optical system and method
GB201711699D0 (en) 2017-07-20 2017-09-06 Univ Bristol Microfluidics analysis system
CN108120851B (zh) * 2017-12-20 2024-03-08 中国医学科学院生物医学工程研究所 一种血液仿体的流速测量方法及仪器
JPWO2021079701A1 (fr) * 2019-10-25 2021-04-29
CN115128299B (zh) * 2022-08-31 2022-12-27 之江实验室 一种测量非透明流场的光声粒子图像测速系统和方法

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS612002A (ja) * 1984-06-14 1986-01-08 Hamamatsu Photonics Kk テレビ画像中の粒子像等の移動を測定する装置
JPH0664074B2 (ja) * 1986-01-08 1994-08-22 株式会社日立製作所 移動物体の速度測定装置
JPH0830711B2 (ja) * 1990-03-28 1996-03-27 満雄 上和野 高速画像処理装置
JPH05264433A (ja) * 1992-03-16 1993-10-12 Mitsuo Kamiwano 気泡測定装置
US5333044A (en) * 1992-11-24 1994-07-26 The United States Of America As Represented By The Department Of Energy Fluorescent image tracking velocimeter
DE4322043C2 (de) * 1993-07-02 1995-07-20 Heidelberg Engineering Optisch Verfahren und Gerät zur Messung der Fließgeschwindigkeit, insbesondere des Blutes
JP3366398B2 (ja) * 1993-10-25 2003-01-14 株式会社日立メディコ Ｍｒｉ装置
JP3313841B2 (ja) * 1993-09-24 2002-08-12 興和株式会社 血流測定装置
US6540981B2 (en) * 1997-12-04 2003-04-01 Amersham Health As Light imaging contrast agents
EP0957750A1 (fr) * 1995-10-23 1999-11-24 Cytometrics, Inc. Procede et appareil d'analyse par imagerie reflechie
US6061176A (en) * 1998-09-14 2000-05-09 Shih; Song Hsin Microscope system for observation and display of microcirculation at multiple body areas
US6473698B1 (en) * 1999-04-22 2002-10-29 University Of Louisville Research Foundation, Inc. Method and apparatus for automated rolling leukocyte velocity measurement in vivo
US20020160470A1 (en) * 2000-11-13 2002-10-31 Genoptix Methods and apparatus for generating and utilizing linear moving optical gradients
FR2842407B1 (fr) * 2002-07-18 2005-05-06 Mauna Kea Technologies "procede et appareillage d'imagerie de fluorescence confocale fibree"
WO2004051310A1 (fr) * 2002-12-02 2004-06-17 Koninklijke Philips Electronics N.V. Outil de segmentation permettant d'identifier des regions de flux dans un systeme d'imagerie
WO2004073501A2 (fr) * 2003-02-20 2004-09-02 Gutin Mikhail Tomographie par coherence optique avec balayage par coherence 3d
US7238150B2 (en) * 2004-06-28 2007-07-03 Houwen Otto H Method for calculating the turbulence factor for a decanting centrifuge

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO2005098474A1 *

Also Published As

Publication number Publication date
US20080045848A1 (en) 2008-02-21
CA2561770C (fr) 2015-06-23
JP4758418B2 (ja) 2011-08-31
ATE479909T1 (de) 2010-09-15
DE602005023292D1 (de) 2010-10-14
AU2005230544B2 (en) 2011-01-06
US7911590B2 (en) 2011-03-22
ES2352559T3 (es) 2011-02-21
IL178347A (en) 2013-12-31
WO2005098474A1 (fr) 2005-10-20
CA2561770A1 (fr) 2005-10-20
CN100592107C (zh) 2010-02-24
JP2007530197A (ja) 2007-11-01
CN1957266A (zh) 2007-05-02
AU2005230544A1 (en) 2005-10-20
FR2868279A1 (fr) 2005-10-07
IL178347A0 (en) 2007-02-11
EP1740974B1 (fr) 2010-09-01
FR2868279B1 (fr) 2006-06-09

Similar Documents

Publication Publication Date Title
CA2561770C (fr) Procede et systeme de mesure de vitesse du flux sanguin
CA2491748C (fr) Procede et appareillage d&#39;imagerie de fluorescence haute resolution par fibre optique et notamment d&#39;imagerie confocale
EP3132237B9 (fr) Appareil et procédé de microscopie à balayage de faisceau optique
EP3343201B1 (fr) Procédé et dispositif de comptage et de caractérisation de particules dans un fluide en mouvement
FR2725532A1 (fr) Microscope autofocus
FR2636429A1 (fr) Appareil de mesure de particules
FR2893132A1 (fr) Dispositif d&#39;analyse a balayage d&#39;echantillons biologiques par fluorescence
EP2795390B1 (fr) Microscopie optique non-linéaire quantitative utilisant un faisceau mis en forme
US20200142174A1 (en) Microscopy method and apparatus for optical tracking of emitter objects
FR2488691A1 (fr) Procede et dispositif pour la detection et la quantification d&#39;agglutinats en temps reel
EP2356429B1 (fr) Dispositif d&#39;analyse d&#39;un melange polyphasique via un faisceau de lumiere retrodiffusee par celui-ci
EP2649431B1 (fr) Systeme et procede d&#39;imagerie multitechniques pour l&#39;analyse chimique, biologique ou biochiimique d&#39;un echantillon.
WO2022017784A1 (fr) Systèmes et procédés d&#39;analyse microscopique d&#39;un échantillon
FR2707018A1 (fr)
JP7270907B2 (ja) 細胞観察システムおよび細胞観察方法
WO2023275664A1 (fr) Procédé d&#39;imagerie, endoscope et produit programme d&#39;ordinateur
WO2021219629A1 (fr) Procede d&#39;obtention du profil d&#39;une surface en mouvement par rapport au systeme
JP2023138481A (ja) 単一分子局在化顕微鏡法のための方法
WO2017064440A1 (fr) Dispositif optique de détection
EP3170045A1 (fr) Systeme de trajectographie automatique pour microscope optique base sur une camera rapide intelligente

Legal Events

Date Code Title Description
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

17P Request for examination filed

Effective date: 20061102

AK Designated contracting states

Kind code of ref document: A1

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

AX Request for extension of the european patent

Extension state: AL BA HR LV MK YU

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

GRAC Information related to communication of intention to grant a patent modified

Free format text: ORIGINAL CODE: EPIDOSCIGR1

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

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

AX Request for extension of the european patent

Extension state: AL BA HR LV MK YU

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: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: FRENCH

REF Corresponds to:

Ref document number: 602005023292

Country of ref document: DE

Date of ref document: 20101014

Kind code of ref document: P

REG Reference to a national code

Ref country code: NL

Ref legal event code: VDEP

Effective date: 20100901

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

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: 20100901

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: 20100901

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: 20100901

REG Reference to a national code

Ref country code: ES

Ref legal event code: FG2A

Effective date: 20110209

LTIE Lt: invalidation of european patent or patent extension

Effective date: 20100901

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: 20100901

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: 20100901

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: 20100901

REG Reference to a national code

Ref country code: IE

Ref legal event code: FD4D

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: 20100901

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: 20101202

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: 20100901

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 FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20100901

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

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: 20100901

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: 20110101

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: 20100901

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: 20100901

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: 20110103

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: 20100901

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

26N No opposition filed

Effective date: 20110606

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: 20100901

BERE Be: lapsed

Owner name: MAUNA KEA TECHNOLOGIES

Effective date: 20110331

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602005023292

Country of ref document: DE

Effective date: 20110606

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 NON-PAYMENT OF DUE FEES

Effective date: 20110331

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: BE

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

Effective date: 20110331

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

Ref country code: LI

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

Effective date: 20110331

Ref country code: CH

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

Effective date: 20110331

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

Ref country code: LU

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

Effective date: 20110325

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: 20100901

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: 20101201

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

Effective date: 20100901

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

Ref country code: ES

Payment date: 20150326

Year of fee payment: 11

Ref country code: IT

Payment date: 20150324

Year of fee payment: 11

REG Reference to a national code

Ref country code: DE

Ref legal event code: R082

Ref document number: 602005023292

Country of ref document: DE

Representative=s name: GRAMM, LINS & PARTNER PATENT- UND RECHTSANWAEL, DE

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 12

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

Ref country code: IT

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

Effective date: 20160325

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 13

REG Reference to a national code

Ref country code: ES

Ref legal event code: FD2A

Effective date: 20170427

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

Ref country code: ES

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

Effective date: 20160326

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 14

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

Ref country code: DE

Payment date: 20230330

Year of fee payment: 19

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

Ref country code: GB

Payment date: 20240328

Year of fee payment: 20

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

Ref country code: FR

Payment date: 20240328

Year of fee payment: 20