WO2009046985A2 - System and method for examining an object - Google Patents

System and method for examining an object Download PDF

Info

Publication number
WO2009046985A2
WO2009046985A2 PCT/EP2008/008547 EP2008008547W WO2009046985A2 WO 2009046985 A2 WO2009046985 A2 WO 2009046985A2 EP 2008008547 W EP2008008547 W EP 2008008547W WO 2009046985 A2 WO2009046985 A2 WO 2009046985A2
Authority
WO
WIPO (PCT)
Prior art keywords
images
sequence
light
detector
image
Prior art date
Application number
PCT/EP2008/008547
Other languages
French (fr)
Other versions
WO2009046985A3 (en
Inventor
Christoph Hauger
Werner Nahm
Theo Lasser
Marcel Leutenegger
Erica Martin-Williams
Antonio Lopez
Original Assignee
Carl Zeiss Surgical Gmbh
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 Carl Zeiss Surgical Gmbh filed Critical Carl Zeiss Surgical Gmbh
Publication of WO2009046985A2 publication Critical patent/WO2009046985A2/en
Publication of WO2009046985A3 publication Critical patent/WO2009046985A3/en
Priority to US12/757,000 priority Critical patent/US8929974B2/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/6408Fluorescence; Phosphorescence with measurement of decay time, time resolved fluorescence
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026Measuring blood flow
    • A61B5/0261Measuring blood flow using optical means, e.g. infrared light
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/14546Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring analytes not otherwise provided for, e.g. ions, cytochromes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47Scattering, i.e. diffuse reflection
    • G01N21/4795Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/6428Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/645Specially adapted constructive features of fluorimeters
    • G01N21/6456Spatial resolved fluorescence measurements; Imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B2018/00315Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body for treatment of particular body parts
    • A61B2018/00434Neural system
    • A61B2018/00446Brain
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • A61B2576/02Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/026Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0042Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
    • A61B5/14553Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases specially adapted for cerebral tissue
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/645Specially adapted constructive features of fluorimeters
    • G01N2021/6482Sample cells, cuvettes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30016Brain
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection

Definitions

  • the present invention relates to a system and a method for examining an object.
  • this invention relates to a system and a method for examining an object containing a flowing liquid.
  • this invention relates to a system and a method for examining a blood motion in a body part of a patient.
  • this invention relates to a system and a method for displaying functional regions in a brain of a patient.
  • measuring light is directed to a region of a body part of a patient to be examined.
  • the measuring light penetrates some millimetres into the region to be examined and interacts with matter in the region within a penetration depth of the measuring light.
  • the measuring light is reflected from particles moving relative to each other, wherein the measuring light experiences a change of a wavelength compared to the wavelength of the measuring light being incident, wherein the change depends on the motion of the reflecting particles, in particular on the motion of blood cells.
  • These changes of a wavelength of the reflected measuring light may be determined by detection using an appropriate detector.
  • US 5,361,769 describes a scanning method for measuring blood perfusion of a body part to be examined. Thereby, a laser beam scans the body part to be examined, wherein the reflected signal is detected by a detector.
  • US 2005/0187477 Al describes a laser-Doppler-perfusion imaging system to gain perfusion related data of a sample.
  • the body part to be examined is illuminated across an area using laser light and the illuminated region is imaged to a detector extending in an area. From a plurality of such recorded images of the body part to be examined, information about a concentration of moving particles and about a velocity of particles moving relative to each other or relative to fixed structures in a sample volume which is associated to each pixel of the detector can be gained.
  • the present invention has been accomplished taking the above problems into consideration.
  • An exposure time for each recording of an image of the sequence of first images is smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms . This is advantageous in order to prohibit that temporal oscillations of the received signals for each pixel of the image are averaged out over a time interval longer than 1 ms .
  • Temporal changes of image values in a high frequency range above 1 kHz, in particular above 10 kHz, are used to obtain a first analysis image value for each pixel. Changes of image values in a low frequency range below 100 Hz are used to obtain a second analysis image value for each pixel.
  • the first analysis image value may comprise perfusion data of the illuminated object.
  • the second analysis image value represents temporal changes of image values in the low frequency range below 100 Hz.
  • this second analysis image value is in particular sensitive to changes of detected intensities of light emanating from the object spanning a time interval of longer than 10 ms .
  • Determining the second analysis image value may thereby comprise temporally filtering the image values to eliminate portions of high frequencies above 100 Hz from the temporal changes of the image values.
  • the second analysis image value may also be obtained by selecting at least two images from the at least one sequence of first images and computationally combining the corresponding pixels. Further, averaging methods for averaging image values over certain time ranges may be applied.
  • temporal changes of image values are determined in two different time scales.
  • processes occurring in these two different time scales may be examined.
  • the first analysis image value may represent information about a motion of particles, in particular blood cells, in the illuminated object
  • the second analysis image value may represent other properties of the object, in particular properties of a liquid flowing through the object. This may for example comprise a composition of the liquid, a concentration of certain particles in the liquid, or a ratio of concentrations of two different particles within the liquid.
  • the measuring light illuminating the object comprises wavelengths in a measuring light wavelength range from one of 605 nm to 630 nm, 560 nm to 580 nm, and 794 nm to 814 nm and the recording the first images comprises detecting measuring light, wherein a ratio of an intensity of the detected measuring light having wavelengths outside the measuring light wavelength range from a total intensity of the detected measuring light amounts to less than 10 %.
  • the deoxyhaemoglobin absorbs at least three times stronger than oxyhaemoglobin.
  • Light of this wavelength range reflected from a body part filled with or flown through by blood may therefore indicate a ratio of deoxyhaemoglobin and oxyhaemoglobin in the blood.
  • the second analysis image value may indicated a change of a degree of oxygenation of haemoglobin in the blood. Changes in the degree of oxygenation may indicate activities of muscle cells, nerve cells or organs of an examined body part of a patient.
  • a further measured parameter is provided to examine blood flow and/or an oxygen transport in the blood.
  • this may be advantageous for displaying and identifying functional regions in the brain of a patient.
  • an increase of deoxyhaemoglobin may be detected in certain regions of the brain, further an increase of the tissue perfusion leading to washing out the deoxyhaemoglobin and thus leading to an increase of the oxyhaemoglobin.
  • a reaction of the organism in certain regions of the brain for example in response to a stimulation, may be temporally examined, in particular with respect to a temporal chronology or a mutual causality of the different processes.
  • Functional active regions may be for example a sensory region, a motor region, an acoustic region, a speech region, or a visual region, but also regions not yet known.
  • an isobestic point of the haemoglobin is located.
  • oxyhaemoglobin and deoxyhaemoglobin absorb incident light of this wavelength to a same extent.
  • detection in the aforementioned wavelength range allows to determine changes of a total haemoglobin content in the examined object.
  • a total haemoglobin content may locally vary depending on an activity of cells, such as nerve cells, etc. Therefore, detection in the aforementioned measuring light wavelength range is also suited for examination of functional regions in the brain of a patient. In order not to receive disturbing signals, measuring light having wavelengths outside the measuring light wavelength range is largely prohibited to contribute to the detected measuring light.
  • Indocyanin green absorbs light in a wavelength range from 790 nm to 820 nm, wherein a maximum of an absorption is at about 805 nm.
  • Indocyanin green (ICG) is a fluorescence dye that may be injected into a blood stream of a patient. Indocyanin green fluoresces in a wavelength range from 800 nm to 900 nm, in particular between 820 nm and 850 nm, further in particular between 830 nm and 840 nm, wherein a fluorescence maximum is at about 835 nm. According to an embodiment of the present invention therefore the object is illuminated with measuring light which is suited to excite fluorescence emission of ICG.
  • the excitation light is thereby used to obtain a first analysis image value according to another embodiment of the present invention which first analysis image value may represent perfusion data of the object containing a flowing liquid.
  • first analysis image value may represent perfusion data of the object containing a flowing liquid.
  • the recording the at least one sequence of first images by the first detector comprises detecting measuring light in a wavelength range from 790 nm to 820 nm.
  • perfusion data about the object containing a flowing liquid may be gained from the first analysis image value. It is advantageous to use measuring light in a first predetermined wavelength range from 560 nm to 580 nm or from 605 nm to 630 nm or from 794 nm to 814 nm. Thereby, beside the perfusion data, information about a degree of oxygenation of haemoglobin in the blood or about a total concentration of haemoglobin in the blood may be gained.
  • the method further comprises: recording at least one second image of the illuminated object by detecting measuring light in a second predetermined wavelength range from 560 nm to 580 nm or from 605 nm to 630 nm or from 794 nm to 814 nm, wherein a ratio of an intensity of detected measuring light having wavelengths outside the second predetermined wavelength range from a total intensity of detected measuring light is smaller than 10 %, wherein the first predetermined wavelength range is different from the second predetermined wavelength range.
  • measuring light is detected in a second predetermined wavelength range which second predetermined wavelength range is different from the first predetermined wavelength range.
  • the first analysis image value associated to pixels (i,j) in the sequence of first images, respectively is dependent on temporal changes of image values of the pixels in the at least one sequence of first images in a frequency range above 1 kHz, in particular above 10 kHz.
  • the recording the at least one sequence of first images comprises recording of temporally subsequent images having a temporal interval of smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.05 ms .
  • a temporal interval of smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.05 ms .
  • the recording the first images comprises detecting measuring light, wherein the ratio of an intensity of the detected measuring light having wavelengths larger than 850 nm and smaller than 500 nm from a total intensity of the detected measuring light amounts to less than 10 %.
  • the measuring light exhibits a coherence length corresponding to at least twice a penetration depth of the measuring light into the object.
  • a penetration depth of the measuring light into the object refers to a path within the object, wherein upon traversal of the measuring light through this path of the object an intensity of the measuring light drops to a fraction of 1/e (1/2.71828) of an intensity of light being incident onto the object.
  • a coherence length of the measuring light is an average length of a continuous wave train. Portions of such a wave train are in a fixed phase relation relative to each other and thus may form a defined interference pattern upon superposition.
  • measuring light on one hand being reflected at the surface of the object and on the other hand being reflected at a layer situated deeper by the penetration depth can interfere with each other.
  • reflected light capable of interference evolves from a depth within the object corresponding to the penetration depth of the measuring light.
  • the method further comprises illuminating the object with white light and imaging and detecting the object illuminated with white light to determine a white light image of the object.
  • the parameters determined from the first analysis image values and the second analysis image values and the second image can be put in relation to structures of the object, in particular organs, identifiable in the visible light.
  • regions of high or low perfusion, regions of high or low ICG concentration, regions of high or low oxygenation or regions of high or low haemoglobin concentration with blood vessels or other relevant structures of an examined body part.
  • These images may also be supplied to a microscope, in particular a surgical microscope, by known techniques.
  • ICG concentration data, data of a degree of oxygenation also preoperatively acquired data, such as nuclear spin data, functional nuclear spin data or x-ray data may be supplied for concurrent display.
  • the first analysis image value associated to a pixel represents a value comprising one of a perfusion, a concentration, an average velocity, and a measure for a velocity distribution, in particular a standard deviation of the velocity distribution, of particles moving relative to each other, in particular blood cells, and a combination of the former.
  • the detected temporal changes of image values may be caused by reflection of measuring light from a plurality of particles moving relative to each other and subsequent superposition.
  • a motion of a particle with a particular velocity involves that light reflected from this particle exhibits a wavelength different from the wavelength of the incident light.
  • This change in the wavelength involves a change in the frequency of the light reflected from the moving particle.
  • k is the wave vector (2 ⁇ / ⁇ ) of the light having wavelength ⁇ and v is the velocity of the particle with respect to a light source emitting the measuring light.
  • the fluctuations or temporal changes of the detected image values represent frequency shifts of reflected measuring light being reflected from moving particles.
  • a spectrum of frequency shifts may be obtained. This may comprise applying Fourier methods.
  • Determining a concentration of particles moving relative to each other may comprise integrating the spectrum of the frequency shifts.
  • Determining a perfusion may comprise computing a first moment of the spectrum of the frequency shifts.
  • Determining an average velocity may comprise forming a ratio of a perfusion and a concentration.
  • the method further comprises at least one of displaying, in particular in superposition, and storing at least one of the first images, the first analysis image values, the second analysis image values, and the white light image for a plurality of pixels.
  • the method further comprises repeating the recording the sequence of first images as temporal sequence and respectively evaluating the repeatedly recorded sequence of first images.
  • the object is a part of a brain of a patient, in particular an aneurysm in the head of the patient.
  • functional region in the brain of a patient may be effectively displayed or identified.
  • first analysis image value second analysis image value, second image
  • a method according to an embodiment of the present invention may be applied.
  • the method further comprises stimulating the patient.
  • the stimulating may thereby comprise requesting to perform a particular action, requesting to speak a word, requesting to move a hand or a finger or the like.
  • the stimulating may be an acoustic stimulating, visual stimulating or may be causing pain.
  • the stimulating may thereby be repeatedly performed, in equal or not equal intervals, wherein regular or irregular sections of a pause may be inserted.
  • the method further comprises: displaying a region of the object in dependence of a result of comparing one of the first analysis value, the second analysis value, and the image value of the at least one second image with preset threshold values for pixels to which the region of the object is imaged.
  • preset threshold values for pixels to which the region of the object is imaged.
  • functional regions in a brain of a patient may be identified by detecting that one of the first analysis value, the second analysis value and the image value of the at least one second image exceeds or falls below preset threshold values. This may indicate that a perfusion, a degree of oxygenation, a haemoglobin concentration or a ICG concentration or their temporal changes are above or below an expected value.
  • the method further comprises injecting indocyanin green into a blood stream of a patient.
  • the system further comprises a first filter for spectrally filtering light for recording images of the illuminated object by the first detector, wherein the first filter is configured such that a ratio of an intensity of light outside a measuring light wavelength range from 605 nm to 630 nm, 560 nm to 580 nm, or 794 nm to 814 nm emanating from the filter from a total intensity of light outside this wavelength range being incident onto the first filter amounts to less than 10 %.
  • Such a configured first filter allows to detect measuring light which is, besides the determination of a perfusion, suitable for determination of a degree of oxygenation of haemoglobin or for determination of a total concentration of haemoglobin. Disturbing light wavelengths are thus efficiently suppressed.
  • the system is suited to excite a fluorescence of indocyanin green (ICG).
  • ICG indocyanin green
  • the second detector is configured to detect fluorescence light emitted by ICG, wherein excitation light in the wavelength range between 790 nm and 820 nm is detected with ten times lower sensitivity than fluorescence light in a wavelength range from 820 nm to 850 nm.
  • a first filter is arranged in a beam path of the measuring light in a direction of travelling of the measuring light in front of the first detector, wherein a ratio of an intensity of light in a wavelength range between 820 nm and 850 nm emanating from the first filter from a total intensity of light outside the wavelength range being incident onto the first filter amounts to less than 10 %.
  • the second detector comprises a second filter arranged in a beam path of the measuring light in a direction of travelling of the measuring light in front of the second detector, wherein a ratio of an intensity of light in a wavelength range between 820 nm and 790 nm emanating from the second filter from a total intensity of light in the wavelength range being incident onto the second filter amounts to less than 10 %.
  • the system is suited to determine perfusion data as well as data about a degree of oxygenation or a total content of haemoglobin concurrently for the same object region.
  • the system further comprises a second detector for recording at least one second image of the illuminated object, wherein the second detector is adapted to detect measuring light in a second predetermined wavelength range from 560 nm to 580 nm or from 605 nm to 630 nm or from 794 nm to 814 nm and wherein the second detector is configured such that light in the second predetermined wavelength range having a further initial intensity emanating from the illuminated object along a beam path in a direction towards the second detector is detected as a further first intensity and such that light outside the second predetermined wavelength range having the further initial intensity emanating from the illuminated object along the beam path in the direction towards the second detector is detected by the second detector as a further second intensity, wherein the further first intensity is ten times, in particular one hundred times higher than the further second intensity and wherein the first predetermined wavelength range is different from the second predetermined wavelength range.
  • the evaluation system is further adapted to associate pixels (i,j) in the sequence of first images (1 ⁇ (1,J)) to which same locations (x,y) of the object are imaged with a second analysis image value, respectively, wherein the second analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1 ⁇ (1,J)) in a frequency range below 100 Hz.
  • temporal changes of a concentration of ICG, of a degree of oxygenation and a total content of haemoglobin may be determined.
  • the first analysis image value respectively associated with pixels in the sequence of first images is dependent of temporal changes of image values of the pixels of the at least one sequence of first images in a frequency range above 1 kHz, in particular above 10 kHz.
  • the first detector comprises a first filter arranged in a beam path of measuring light detected by the first detector, wherein a ratio of an intensity of light outside the first predetermined wavelength range emanating from the first filter from a intensity of light outside the first predetermined wavelength range being incident onto the filter amounts to less than 10 %.
  • the evaluation system comprises a low path filter to determine temporal changes of image values of the pixels in the at least one sequence of first images in a frequency range below 100 Hz.
  • the evaluation system is adapted to record images having a temporal interval of smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.05 ms .
  • a detector allowing to record images having a temporal interval of 0.05 ms, allows a detection of measuring light signals in a frequency range up to 10 kHz.
  • the detected image values represent a light signal sampled with a sampling frequency of 20 kHz.
  • a maximal velocity of a moving particle that may be determined using such a detector is 4 mm/s when utilizing measuring light having a wavelength of 800 nm. This velocity approximately corresponds to blood flow velocities in capillaries and micro-capillaries of a patient.
  • the light source is adapted to generate light having wavelengths in the region from 500 nm to 850 nm.
  • the light source is adapted to generate light having a coherence length, comprising a temporal coherence length, of greater than 2 mm.
  • the system further comprises an illumination apparatus for illuminating the object with white light and an imaging apparatus for imaging the object illuminated with white light to a white light image of the object and further a display apparatus for displaying images, in particular in superposition.
  • the system further comprises a storage system for storing at least one of first images of the illuminated object, first analysis image values, second analysis image values, the second image and the white light image.
  • the system further comprises a display apparatus for displaying at least one of the recorded first images, the first analysis image values, the second analysis image values, the second image, and the white light image.
  • Figure 1 illustrates an embodiment of a method for examining an object according to the invention
  • Figure 2 illustrates an evaluation method of a method for examining an object according to an embodiment of the present invention
  • Figure 3 illustrates a system for examining an object according to an embodiment of the present invention
  • Figure 4 illustrates a further system for examining an object according to an embodiment of the present invention
  • Figure 5 illustrates an absorption spectrum and a fluorescence emission spectrum of indocyanin green (ICG) ;
  • Figure 6 illustrates a system for examining an object according to an embodiment of the present invention
  • Figure 7a illustrates a white light image
  • Figure 7b illustrates an image of first analysis values (perfusion data) of an object, wherein the images were determined according to a method for examining an object according to an embodiment of the present invention
  • Figure 8a, Figure 8b, Figure 8c, Figure 8d, Figure 8e, Figure 8f show second images determined by detection of fluorescence of ICG, which were determined according to a method according to an embodiment of the present invention
  • Figure 9 illustrates a system for examining an object according to an embodiment of the present invention.
  • Figure 10a, Figure 10b, Figure 10c, Figure 1Od illustrate a method according to an embodiment of the present invention for examining an object.
  • Figure 1 shows a method for examining an object according to an embodiment of the present invention.
  • a sequence of images Bi, ... , Bjc, ... , B N is recorded from an object illuminated with measuring light.
  • a temporal interval t ⁇ between the recording of the first image B ⁇ and the recording of the second image B2 amounts to 0.05 ms.
  • the sequence of images B]_, ... , B ⁇ , ... , B N is illustrated in Figure Ia.
  • the image B ⁇ has an image value Ik ; (i,j)) at the location defined by the indices i and j.
  • the first index of the pixels of the image runs from 1 to L and the second index of the pixels of the image runs from 1 to M.
  • a temporal course of image values I(t,i,j) may be obtained.
  • An example of such a temporal course of image values I(t,i,j) is illustrated in Figure Ic. In this Figure, the image value I(t,i,j) of a pixel (i,j) is plotted versus the time t.
  • temporal course I(t,i,j) comprises changes or fluctuations which are temporarily spaced apart less than 1 ms (high frequency portion of the temporal changes of the image values) which are superimposed by temporal changes which are temporarily spaced apart in the order of 100 ms or greater (low frequency portion of the temporal changes of image values).
  • Temporal changes of image values being temporally spaced apart shorter than 0.5 ms are comprised in the image values but are not illustrated here.
  • Temporal changes of image values being spaced apart longer than I s or 10 s are comprised in the temporal changes of the image values but can not be illustrated due to the chosen section of a temporal development of image values.
  • Figure 2 schematically illustrates an evaluation method according to an embodiment of the present invention.
  • the left hand side of Figure 2a illustrates as a graph 11 temporal changes of images values in a frequency range above 1 kHz, in the following also denoted as high frequency range of temporal changes of the image values, or simply high frequency range.
  • this portion of temporal changes of image values carries information about a motion of particles within the illuminated object from which measuring light was reflected.
  • a further evaluation of this high frequency portion of temporal changes of image values comprises establishing a Fourier spectrum or power spectrum of temporal changes of image values.
  • the Fourier spectrum or the power spectrum may further be analysed with respect to its moments. Thereby, moments of this spectrum are for example an average value or an area.
  • first analysis values a concentration of particles moving relative to each other, a velocity of particles moving relative to each other and a perfusion of these particles may be obtained which parameter are denoted as first analysis values.
  • B ⁇ , ... , B ⁇ , ... , B N an image 21 of first analysis values is obtained as illustrated in the right hand side of Figure 2a.
  • the values of first analysis values within the image 21 of first analysis values are illustrated as contour lines 19 in the illustrated example.
  • Figure 2b shows a further evaluation step of the inventive method.
  • the temporal course of image values I(t,i,j) of the temporal course of image values illustrated in Figure Ic is illustrated in a frequency range below 100 Hz as curve 12.
  • the curve 12 does not comprise anymore the very fast temporal changes of image values occurring in time intervals below 1 ms, in particular 0.1 ms .
  • the temporal changes of image values may be analysed in different low frequency ranges, such as below 1 Hz, below 10 Hz, or below 100 Hz.
  • the curve 12 may be determined in different manners. For example it is possible, to apply a low pass filter to the image value I(t,i,j) to obtain the curve 12.
  • At least two images may be selected from the sequence of images B ⁇ , ... , B ⁇ , ... , B N to obtain the curve 12.
  • This may for example also comprise interpolation methods to determine values of image values at intermediate locations not measured.
  • an image 22 of second image values is obtained again showing contour lines 23. It is apparent that a shape of the contour lines 19 of the image 21 have a form different from a form of the contour lines 23 of the image 22 in Figure 2.
  • FIG. 3 illustrates a system 1 for examining an object according to an embodiment of the present invention.
  • the system 1 comprises an illumination apparatus 5, a beam splitter 24, an objective 26, a position acquisition apparatus 27, a first camera 14, a controller 42, and an evaluation system 40.
  • the illumination apparatus 5 comprises a laser 6 for generating light having a coherence length of at least 2 mm.
  • the laser 6 is tuneable to be able to generate laser light having different wavelengths.
  • the light 7 generated by the laser 6 is collimated by the collimator 8 to form measuring light 32. Measuring light 32 traverses the semi-transparent mirror 24 and the objective 26 to illuminate the object 28 in the object plane 30 with measuring light 32.
  • a position (x,y) of the object 28 is registered by the position acquisition apparatus 27 thereby allowing a conversion of coordinates in an object coordinate system into coordinates of a coordinate system of the system 1.
  • the light 34 emanating from a point (x,y) of the object 28 traverses the objective 26 and is reflected at the semi-transparent mirror 24 to enter the first camera 14.
  • the first camera 14 comprises a first filter 16 to filter out undesired wavelengths from the light 34 emanating from the object or to filter out stray light from the surrounding.
  • the filtered light 36 traverses the camera optics 18 imaging the filtered light 36 to a pixel (i,j) of the first detector 20. Different locations of the object 28 thus are imaged to different pixels of the first detector.
  • the first detector 20 is capable to record images having an exposure time of 0.05 ms . A temporal interval between recordings of two images amounts to 0.1 ms.
  • the first detector 20 and also the first filter 16 are connected with the camera controller 42.
  • the camera controller 42 is capable to send signals to the first filter 16 to change a filter characteristics of the first filter 16. This may for example be realized in that a wheel having different filters arranged in a circumferential direction is rotated to bring a desired filter into a beam path of the first camera 14.
  • the controller 42 may control the tuneable laser to choose a desired wavelength of the measuring light 32.
  • measuring light having different wavelengths may be advantageous for examining the object.
  • Wavelengths in a measuring light wavelength range from 605 nm to 630 ran, 560 nm to 580 ran or 794 nm to 814 nm may be utilized to measure a degree of oxygenation of haemoglobin or to measure a concentration of haemoglobin within the sample.
  • Images recorded by the first detector are acquired by the controller 42 and transferred to the evaluation system 40.
  • the evaluation system 40 comprises a processing system 44, a display device 46, a storage system 48 and an input system 50.
  • the processing system 44 comprises a computer capable of running image processing software.
  • the image processing software enables to read in a sequence of images generated by the first detector 20.
  • the software is adapted to evaluate the sequence of images per pixel to analyse temporal changes of image values I(t,i,j).
  • the image processing software is controlled via the input system 50 comprising a console and a cursor.
  • On the display device 46 comprising a screen the course of the analysis of the sequence of images may be monitored.
  • the evaluating the sequence of images per pixel comprises computing a spectrum of the temporal changes of image values and computing moments of this spectrum. Further, the temporal changes of image values are filtered with respect to their frequency range. Result images of the analysis may be displayed on the display device 46, in particular in superposition.
  • FIG. 4 shows a system 1 for examining an object according to an embodiment of the present invention.
  • the light 7 emitted by the laser 6 being collimated by the collimator 8 and leaving same as measuring light 32 does not travel along a common beam path together with the light 34 emanating from the object 28.
  • Other elements of the system 1 illustrated in Figure 4 are analogous to those of the system 1 illustrated in Figure 3 and thus a repeated description is avoided.
  • Figure 5 illustrates an excitation spectrum and a fluorescence emission spectrum of indocyanin green (ICG) .
  • Curve 60 shows an excitation spectrum of ICG, wherein a normalized absorption is plotted versus a wavelength. ICG absorbs light in a wavelength range from about 700 nm to 850 nm. Thereby, a maximum of absorption occurs at around 805 nm. By illuminating ICG with light between 790 nm and 820 nm a fluorescence emission of ICG is excited.
  • Curve 62 shows a fluorescence emission spectrum of ICG, wherein a normalized fluorescence is plotted against a wavelength. It is apparent that ICG is fluorescent in a wavelength range from about 800 nm to 900 nm.
  • a maximum of fluorescence emission is observed at around 835 nm.
  • FIG. 6 illustrates a system 1 for examining an object according to an embodiment of the present invention.
  • the system 1 comprises an illumination apparatus 5, a semi-transparent mirror 24, a semi-transparent mirror 24', an objective 26, a first camera 14, a second camera 64 and an evaluation system 74.
  • the illumination apparatus 5 comprises a laser 65 capable to generate light in a wavelength range from 790 nm to 820 nm.
  • the light 7 emitted by the laser 65 is collimated by the collimator 8 to form measuring light 32.
  • Measuring light 32 traverses the semi-transparent mirror 24 and the semi-transparent mirror 24' to traverse the objective 26 and illuminate the object 28.
  • the first camera 14 comprises a first filter 16 having a transmission characteristics to allow traversing light having wavelengths suited for excitation of fluorescence of ICG, i.e. in a wavelength range between 790 nm and 820 nm, as filtered light 36. Filtered light 36 traverses the camera optics 18 to be imaged onto a pixel (i,j) of the first detector.
  • the first detector 20 is capable to record images having an exposure time of 0.05 ms. A temporal interval between the recordings of two images amounts to 0.1 ms.
  • the light entering the second camera 64 having emanated from the location (x,y) of the object 28 as light 34 traverses a second filter 66 which is designed to attenuate light used for excitation of fluorescence emission of ICG, i.e. light having wavelengths between 790 nm and 820 nm, upon traversal through the filter 66 strongly in its intensity.
  • the filter 66 is designed such that only 10 % of an intensity of light in a wavelength range between 790 nm and 820 nm is allowed to traverse.
  • the filter 66 is transparent for fluorescence light of ICG, i.e. transparent in a wavelength range between about 820 nm and 850 nm.
  • the filtered light 86 emanating from the filter 66 traverses the camera optics 68 to be imaged onto a pixel (i,j) of the second detector 70.
  • the first detector 20 is capable to subsequently record images having a very short temporal interval (0.1 ms)
  • the second detector 70 records images having much larger temporal intervals (larger than 10 ms) .
  • a CCD-detector may be employed as a second detector 70.
  • the first detector 20 however may be a CMOS- detector.
  • the controller 72 controls the recording of images of the first detector 20 as well as images of the second detector 70 and acquires the thus obtained first and second images.
  • the first and second images are provided to the evaluation system 74, in particular to the processing system 76.
  • the processing system 76 comprises as the processing system 44 of the embodiment illustrated in Figure 3 an image processing software.
  • the image processing software is however also configured to process second images generated by the second detector 70.
  • the processing the images acquired by the first detector 20 is simplified in comparison to the processing the first images of the embodiment illustrated in Figure 3.
  • second analysis image values must be determined from a sequence of recorded first images, wherein the second analysis image values depend on temporal changes of image values in the sequence of first images in a frequency range below 100 Hz.
  • first analysis values are determined from the temporal changes of image values of the sequence of recorded first images, wherein the first analysis values represent perfusion data.
  • information about a motion of particles within the illuminated object are determined. Due to the choice of the measuring light being suited for excitation of a fluorescence of ICG, and the choice of the filter characteristics of the second filter 66, the second detector 70 substantially detects fluorescence light of ICG emanating from the object 28.
  • the embodiments illustrated in the Figures 3, 4 and 6 may further comprise a system for imaging the object as white light image.
  • a white light imaging system may comprise an illumination source such as a halogen lamp, an objective, and a further detector such as a CCD camera.
  • the white light imaging system may also be adapted as a microscope, in particular as a surgical microscope. Images generated by the white light imaging system may be provided to the evaluation system 40 and 74, respectively, and may be displayed together with the first and second images, the first and second analysis image values by the display device 46 or/and may be supplied for display in the surgical microscope.
  • Figure 7 shows images determined using systems according to embodiments of the present invention and using methods according to embodiments of the present invention, respectively.
  • Figure 7a shows a white light image of a portion of a brain of a patient.
  • a meshwork of capillaries and micro-capillaries is observed penetrating underlying brain tissue.
  • the bright red capillaries 90 and 92 are arteries and the dark red capillaries 93 represent veins.
  • Reference sign 91 denotes an aneurysm.
  • Figure 7b shows an image determined according to a method of the present invention showing a perfusion of blood cells in a region of the brain of the patient, wherein the region overlaps with the region of the brain illustrated in Figure 7a.
  • Values of the perfusion are represented by pseudo-colour values (or grey values, respectively) , comprising blue, yellow, red in increasing order.
  • pseudo-colour values or grey values, respectively
  • the capillaries 90, 92 and 93 exhibit a lower perfusion than adjacent brain tissue.
  • a possible explanation is that a motion of blood cells in the capillaries 90, 92 and 93 is that fast (such as greater than 10 mm/s) that such a motion cannot be resolved any more by the detector used here.
  • Figures 8a to 8f show images from a region of the brain overlapping with the regions illustrated in Figures 7a and 7b.
  • the images shown in Figure 8 are acquired by detection of fluorescence light of ICG by the second detector 70 illustrated in Figure 6.
  • the Figures 8a to 8f were recorded with a temporal interval of 0.5 s to I s and thus represent a temporal development of a fluorescence of ICG.
  • ICG was injected into the blood stream of the patient. After the injection of ICG initially ( Figures 8a to 8d) an increasing intensity of fluorescence of ICG in the arteries, such as arteries 90 and 92, is observed. That means that ICG is supplied to the brain via the arteries.
  • FIG. 8e shows a condition in which the brain is maximally filled with ICG. Subsequently outflow of ICG from the examined region via veins occurs. This is shown in Figure 8f.
  • vein 93 shows a high fluorescence intensity indicating that it guides blood.
  • aneurysm 91 also indicating that the aneurysm is flown through by blood.
  • Figure 9 shows a system 100 for examining an object according to an embodiment of the present invention.
  • the system comprises a laser 101 for generating laser light which illuminates the object 28 after collimating it by a collimator 105 and reflection at the semi-transparent mirror 131, after traversal through the semi-transparent mirror 130 and traversal through the objective 140.
  • the light of the laser 101 reflected from the object 28 is detected by the detector 109 after traversal through the semi-transparent mirror 130 through the two semi- transparent mirrors 131, 132 and reflection at the semi- transparent mirror 134 and traversal through camera optics 107.
  • Light source 111 emits light comprising wavelengths in the range from 790 nm to 820 nm and is thus suited for excitation of fluorescence of ICG.
  • the light generated by the light source 111 is collimated by the collimator 115, is reflected at the semi- transparent mirror 132, traverses through the semi-transparent mirror 130 and through the objective 140 to illuminate the object 28. There it excites fluorescence of ICG.
  • the fluorescence light of ICG emanating from the object 28 is detected by the detector 119 after traversal through several semi-transparent mirrors and traversal through the camera optics 117.
  • the detector 119 comprises a suitable filter to attenuate light having wavelengths outside a range from 820 nm to 850 nm from light detected by the detector 119 strongly in its intensity.
  • the system 100 further comprises a white light source 121 for illuminating the object 28 with white light.
  • the white light image of the object 28 is detected by the detector 129 after traversal through the objective and several semi- transparent mirrors and traversal through the camera optics 127. Sequences of images recorded by the detector 109 and images recorded by the detectors 119 and 129 are supplied to a monitor and evaluation unit 132.
  • the monitor and evaluation unit 132 is adapted to derive perfusion data from the sequence of images received by the detector 109 and to display these perfusion data together with the white light image and the ICG fluorescence image and to process these further.
  • Figure 10 shows methods according to embodiments of the present invention.
  • Figures 10a to 1Od show different illumination and detection scenarios according to embodiments of the present invention.
  • an illumination strength in dependence of a wavelength is illustrated, respectively, and in the Figures 10a to 1Od on the right hand side a detection sensitivity in dependence of the wavelength is illustrated, respectively.
  • the illumination in different wavelength ranges may be performed by one or more light sources and the detection in different wavelength ranges may be performed by one or more detectors.
  • a white light illumination wavelength range and a white light detection wavelength range is denoted, respectively.
  • an illumination wavelength range for illumination with light suitable for excitation of fluorescence of ICG and a ICG detection wavelength range suited for detection of fluorescence of ICG is denoted, respectively.
  • Reference sign LDI in these Figures denotes an illumination wavelength range and a detection wavelength range suited to detect laser Doppler signals from an illuminated sample which are suitable to obtain perfusion data about the object, respectively.
  • the scenario illustrated in Figure 10a may for example be realized using two light sources and three detectors.
  • the LDI illumination light is concurrently suited to excite fluorescence of ICG in the object.
  • the LDI and/or ICG illumination light may be in a wavelength range overlapping or not overlapping with the illumination wavelength range of the white light illumination.
  • An LDI detection wavelength range may in embodiments overlap with a detection wavelength range for white light detection or not.
  • the LDI detection signal may be separated from the white light detection signal by appropriate processing that may comprise an analysis of temporal changes of image values in a range above 1 kHz.
  • embodiments of the present invention provide a system and a method for examining an object containing a fluid liquid, wherein the object is illuminated with measuring light and images are temporarily shortly subsequently recorded.
  • the images are evaluated per pixel to determine perfusion data from a high frequency portion above 1 kHz and to determine further information about properties of the object from a low frequency portion below 100 Hz, such as a degree of oxygenation of haemoglobin, a concentration of haemoglobin or a concentration of ICG.
  • This information determined by evaluation is displayed in a form of an image in superposition with a white light image of the object.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • Optics & Photonics (AREA)
  • General Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Engineering & Computer Science (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • Cardiology (AREA)
  • Hematology (AREA)
  • Physiology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

A system and a method for examining an object containing a fluid liquid are provided, wherein the object is illuminated with measuring light and images are temporarily shortly subsequently recorded. The images are evaluated per pixel to determine perfusion data from a high frequency portion above 1 kHz and to determine further information about properties of the object from a low frequency portion below 100 Hz, such as a degree of oxygenation of haemoglobin, a concentration of haemoglobin or a concentration of ICG. This information determined by evaluation is displayed in a form of an image in superposition with a white light image of the object.

Description

System and Method for Examining an Object
Background of the invention
Field of the invention
The present invention relates to a system and a method for examining an object. In particular, this invention relates to a system and a method for examining an object containing a flowing liquid. Further, this invention relates to a system and a method for examining a blood motion in a body part of a patient. In particular, this invention relates to a system and a method for displaying functional regions in a brain of a patient.
Brief Description of Related Art
From the state of the art it is known to exploit the so-called "Doppler effect" to optically examine a motion of particles in a liquid. Thereby, measuring light is directed to a region of a body part of a patient to be examined. The measuring light penetrates some millimetres into the region to be examined and interacts with matter in the region within a penetration depth of the measuring light. Among other things, the measuring light is reflected from particles moving relative to each other, wherein the measuring light experiences a change of a wavelength compared to the wavelength of the measuring light being incident, wherein the change depends on the motion of the reflecting particles, in particular on the motion of blood cells. These changes of a wavelength of the reflected measuring light may be determined by detection using an appropriate detector. The basic principles for this so-called "laser Doppler method" are described in Stern, "In vivo evaluation of microcirculation by coherent light scattering", Nature, Vol. 254, pages 56 - 58, March 1975. Information about a motion of particles in a liquid are also referred to as perfusion data or perfusion in the following. US 4,109,647 describes how one may obtain information about the velocity of red blood cells at one point of a body part to be examined using this method.
US 5,361,769 describes a scanning method for measuring blood perfusion of a body part to be examined. Thereby, a laser beam scans the body part to be examined, wherein the reflected signal is detected by a detector.
US 2005/0187477 Al describes a laser-Doppler-perfusion imaging system to gain perfusion related data of a sample. The body part to be examined is illuminated across an area using laser light and the illuminated region is imaged to a detector extending in an area. From a plurality of such recorded images of the body part to be examined, information about a concentration of moving particles and about a velocity of particles moving relative to each other or relative to fixed structures in a sample volume which is associated to each pixel of the detector can be gained.
However, it became apparent that the previously described systems and methods often do not provide satisfying results.
For displaying and identifying functional regions in a brain of a patient so far a number of techniques have been used. In particular, it should be mentioned: functional magnet resonance tomography, electroencephalography, positron emission tomography and others. Due to the "blood oxygen level dependent (BOLD)- effect", it is possible in principle to detect a contrast in a body part to be examined due to the presence or non-presence of oxygenated and deoxygenated, respectively, haemoglobin using magnet resonance tomography. However, this method is very cost consuming and time consuming. Due to a long data recording time required by this method motions of the body part to be examined may occur during data recording. This in turn leads to difficulties in an evaluation of the data, in particular by the rotations and translations associated therewith, or also deformations of the data sets relative to each other which have to be determined to align the data relative to each other. Thus, there is a need to provide a system and a method diminishing the above-mentioned disadvantages.
Summary of the Invention
The present invention has been accomplished taking the above problems into consideration.
Therefore, it is an object of the present invention, to provide a system and a method for examining an object with which an examination of an object, in particular containing a flowing liquid, may be improved.
According to an embodiment of the present invention, a method for examining an object containing a flowing liquid is provided, wherein the method comprises: illuminating the object with measuring light and recording at least one sequence of first images of the illuminated object, wherein an exposure time for each recording of an image of the sequence of first images is smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms; evaluating the at least one sequence of first images (k = 1,...,N), by associating pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range above 1 kHz, in particular above 10 kHz; and evaluating the at least one sequence of first images (k = 1,...,N) by further associating pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a second analysis image value, respectively, wherein the second analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range below 100 Hz.
An exposure time for each recording of an image of the sequence of first images is smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms . This is advantageous in order to prohibit that temporal oscillations of the received signals for each pixel of the image are averaged out over a time interval longer than 1 ms . Temporal changes of image values in a high frequency range above 1 kHz, in particular above 10 kHz, are used to obtain a first analysis image value for each pixel. Changes of image values in a low frequency range below 100 Hz are used to obtain a second analysis image value for each pixel. The first analysis image value may comprise perfusion data of the illuminated object. Herein, in particular information about a concentration of particles moving relative to each other, about a velocity of particles moving relative to each other and about measures derived therefrom is comprised. Hereby, the second analysis image value represents temporal changes of image values in the low frequency range below 100 Hz. Thus, this second analysis image value is in particular sensitive to changes of detected intensities of light emanating from the object spanning a time interval of longer than 10 ms . Determining the second analysis image value may thereby comprise temporally filtering the image values to eliminate portions of high frequencies above 100 Hz from the temporal changes of the image values. The second analysis image value may also be obtained by selecting at least two images from the at least one sequence of first images and computationally combining the corresponding pixels. Further, averaging methods for averaging image values over certain time ranges may be applied.
According to an embodiment temporal changes of image values are determined in two different time scales. Thus, processes occurring in these two different time scales may be examined. While the first analysis image value may represent information about a motion of particles, in particular blood cells, in the illuminated object, the second analysis image value may represent other properties of the object, in particular properties of a liquid flowing through the object. This may for example comprise a composition of the liquid, a concentration of certain particles in the liquid, or a ratio of concentrations of two different particles within the liquid. According to an embodiment of the present invention, the measuring light illuminating the object comprises wavelengths in a measuring light wavelength range from one of 605 nm to 630 nm, 560 nm to 580 nm, and 794 nm to 814 nm and the recording the first images comprises detecting measuring light, wherein a ratio of an intensity of the detected measuring light having wavelengths outside the measuring light wavelength range from a total intensity of the detected measuring light amounts to less than 10 %.
In a wavelength range from 605 nm to 630 nm the deoxyhaemoglobin absorbs at least three times stronger than oxyhaemoglobin. Light of this wavelength range reflected from a body part filled with or flown through by blood may therefore indicate a ratio of deoxyhaemoglobin and oxyhaemoglobin in the blood. Thus, the second analysis image value may indicated a change of a degree of oxygenation of haemoglobin in the blood. Changes in the degree of oxygenation may indicate activities of muscle cells, nerve cells or organs of an examined body part of a patient. Thus, amongst the perfusion data derived from the first analysis image value, a further measured parameter is provided to examine blood flow and/or an oxygen transport in the blood. In particular, this may be advantageous for displaying and identifying functional regions in the brain of a patient. After stimulating or provoking the patient for example an increase of deoxyhaemoglobin may be detected in certain regions of the brain, further an increase of the tissue perfusion leading to washing out the deoxyhaemoglobin and thus leading to an increase of the oxyhaemoglobin. Thus, a reaction of the organism in certain regions of the brain, for example in response to a stimulation, may be temporally examined, in particular with respect to a temporal chronology or a mutual causality of the different processes. Functional active regions may be for example a sensory region, a motor region, an acoustic region, a speech region, or a visual region, but also regions not yet known. In a wavelength range from 560 nm to 580 nm, in particular at 570 nm and in a wavelength range between 794 nm to 814 nm, in particular at 804 nm, an isobestic point of the haemoglobin is located. At the isobestic point of the haemoglobin, oxyhaemoglobin and deoxyhaemoglobin absorb incident light of this wavelength to a same extent. Thus, detection in the aforementioned wavelength range allows to determine changes of a total haemoglobin content in the examined object. A total haemoglobin content may locally vary depending on an activity of cells, such as nerve cells, etc. Therefore, detection in the aforementioned measuring light wavelength range is also suited for examination of functional regions in the brain of a patient. In order not to receive disturbing signals, measuring light having wavelengths outside the measuring light wavelength range is largely prohibited to contribute to the detected measuring light.
According to an embodiment of the present invention, a method for examining an object containing a flowing liquid is provided, wherein the method comprises: illuminating the object with measuring light comprising wavelengths in a range from 790 nm to 820 nm; recording at least one sequence of first images of the illuminated object by a first detector, wherein an exposure time for each recording of an image of the sequence of first images is smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms; evaluating the at least one sequence of first images (k = 1,...,N), by associating pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value , depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)); and recording at least one second image of the illuminated object, wherein the second image is formed by detecting light having wavelengths in a range from 820 nm to 850 nm, in particular 830 nm to 840 nm, by a second detector and wherein a ratio of an intensity of light having wavelengths smaller than 820 nm detected by the second detector from a total intensity of light detected by the second detector amounts to less than 10 %.
Indocyanin green absorbs light in a wavelength range from 790 nm to 820 nm, wherein a maximum of an absorption is at about 805 nm. Indocyanin green (ICG) is a fluorescence dye that may be injected into a blood stream of a patient. Indocyanin green fluoresces in a wavelength range from 800 nm to 900 nm, in particular between 820 nm and 850 nm, further in particular between 830 nm and 840 nm, wherein a fluorescence maximum is at about 835 nm. According to an embodiment of the present invention therefore the object is illuminated with measuring light which is suited to excite fluorescence emission of ICG. The excitation light is thereby used to obtain a first analysis image value according to another embodiment of the present invention which first analysis image value may represent perfusion data of the object containing a flowing liquid. By detecting within the fluorescence emission wavelength range by a second detector information about a concentration of ICG within the object containing a flowing liquid may be obtained. Again, it is thus possible using the inventive method to simultaneously determine two parameters carrying independent information about a flowing liquid.
According to an embodiment, the recording the at least one sequence of first images by the first detector comprises detecting measuring light in a wavelength range from 790 nm to 820 nm.
According to an embodiment of the present invention, a method for examining an object containing a flowing liquid is provided, the method comprising: illuminating the object with measuring light and recording at least one sequence of first images of the illuminated object, wherein an exposure time for each recording of an image of the sequence of first images is smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms; evaluating the at least one sequence of first images (k = 1,...,N) by associating pixels (i,j) in the sequence of first images (Ik(i,j)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (I(i/j))f wherein the recording the at least one sequence of first images of the illuminated object comprises detecting measuring light in a first predetermined wavelength range from one of 560 nm to 580 nm, 605 nm to 630 nm, and 794 nm to 814 nm, wherein a ratio of an intensity of detected measuring light having wavelengths outside the first predetermined wavelength range from a total intensity of detected measuring light is smaller than 10 %.
Similar as in other embodiments of the present invention perfusion data about the object containing a flowing liquid may be gained from the first analysis image value. It is advantageous to use measuring light in a first predetermined wavelength range from 560 nm to 580 nm or from 605 nm to 630 nm or from 794 nm to 814 nm. Thereby, beside the perfusion data, information about a degree of oxygenation of haemoglobin in the blood or about a total concentration of haemoglobin in the blood may be gained.
According to an embodiment of the invention, the method further comprises: recording at least one second image of the illuminated object by detecting measuring light in a second predetermined wavelength range from 560 nm to 580 nm or from 605 nm to 630 nm or from 794 nm to 814 nm, wherein a ratio of an intensity of detected measuring light having wavelengths outside the second predetermined wavelength range from a total intensity of detected measuring light is smaller than 10 %, wherein the first predetermined wavelength range is different from the second predetermined wavelength range.
Thereby, measuring light is detected in a second predetermined wavelength range which second predetermined wavelength range is different from the first predetermined wavelength range. Thus, it is possible to simultaneously determine a degree of oxygenation of haemoglobin as well as a total concentration of haemoglobin in the object.
According to an embodiment of the present invention, the method further comprises evaluating the at least one sequence of first images (Jc = 1,...,N) by further associating pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a second analysis image value, respectively, wherein the second analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range below 100 Hz.
Using this method, it is thus possible to detect temporal changes of an ICG concentration, temporal changes of a degree of oxygenation and temporal changes of a total haemoglobin content in the object.
According to an embodiment of the present invention, the first analysis image value associated to pixels (i,j) in the sequence of first images, respectively, is dependent on temporal changes of image values of the pixels in the at least one sequence of first images in a frequency range above 1 kHz, in particular above 10 kHz.
According to an embodiment of the present invention, the recording the at least one sequence of first images comprises recording of temporally subsequent images having a temporal interval of smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.05 ms . Thus, it is possible to sample and detect a light signal emanating from the object that fast such that the detected image values comprise temporal changes of the light signals in a frequency range up to at least 10 kHz.
According to an embodiment of the present invention, the recording the first images comprises detecting measuring light, wherein the ratio of an intensity of the detected measuring light having wavelengths larger than 850 nm and smaller than 500 nm from a total intensity of the detected measuring light amounts to less than 10 %.
According to an embodiment of the present invention, the measuring light exhibits a coherence length corresponding to at least twice a penetration depth of the measuring light into the object. A penetration depth of the measuring light into the object refers to a path within the object, wherein upon traversal of the measuring light through this path of the object an intensity of the measuring light drops to a fraction of 1/e (1/2.71828) of an intensity of light being incident onto the object. A coherence length of the measuring light is an average length of a continuous wave train. Portions of such a wave train are in a fixed phase relation relative to each other and thus may form a defined interference pattern upon superposition. Thus, according to this embodiment, measuring light on one hand being reflected at the surface of the object and on the other hand being reflected at a layer situated deeper by the penetration depth can interfere with each other. Thus, reflected light capable of interference evolves from a depth within the object corresponding to the penetration depth of the measuring light.
According to an embodiment of the present invention, the method further comprises illuminating the object with white light and imaging and detecting the object illuminated with white light to determine a white light image of the object. By this provision, the parameters determined from the first analysis image values and the second analysis image values and the second image can be put in relation to structures of the object, in particular organs, identifiable in the visible light. Thereby, it is for example possible to correlate regions of high or low perfusion, regions of high or low ICG concentration, regions of high or low oxygenation or regions of high or low haemoglobin concentration with blood vessels or other relevant structures of an examined body part. These images may also be supplied to a microscope, in particular a surgical microscope, by known techniques. Besides the perfusion data, functional region determined therefrom, ICG concentration data, data of a degree of oxygenation, also preoperatively acquired data, such as nuclear spin data, functional nuclear spin data or x-ray data may be supplied for concurrent display.
According to an embodiment of the present invention, the first analysis image value associated to a pixel represents a value comprising one of a perfusion, a concentration, an average velocity, and a measure for a velocity distribution, in particular a standard deviation of the velocity distribution, of particles moving relative to each other, in particular blood cells, and a combination of the former.
The detected temporal changes of image values may be caused by reflection of measuring light from a plurality of particles moving relative to each other and subsequent superposition. A motion of a particle with a particular velocity involves that light reflected from this particle exhibits a wavelength different from the wavelength of the incident light. This change in the wavelength involves a change in the frequency of the light reflected from the moving particle. In particular, this frequency shift Aω (= 2π times the difference between the frequencies of the incident and reflected light) is
Aω = k -v ,
wherein k is the wave vector (2π/λ) of the light having wavelength λ and v is the velocity of the particle with respect to a light source emitting the measuring light.
The fluctuations or temporal changes of the detected image values represent frequency shifts of reflected measuring light being reflected from moving particles. From the image values a spectrum of frequency shifts may be obtained. This may comprise applying Fourier methods. Determining a concentration of particles moving relative to each other may comprise integrating the spectrum of the frequency shifts. Determining a perfusion may comprise computing a first moment of the spectrum of the frequency shifts. Determining an average velocity may comprise forming a ratio of a perfusion and a concentration. According to an embodiment of the present invention, the method further comprises at least one of displaying, in particular in superposition, and storing at least one of the first images, the first analysis image values, the second analysis image values, and the white light image for a plurality of pixels.
According to an embodiment of the present invention, the method further comprises repeating the recording the sequence of first images as temporal sequence and respectively evaluating the repeatedly recorded sequence of first images.
According to an embodiment of the present invention, the object is a part of a brain of a patient, in particular an aneurysm in the head of the patient. Using the inventive method, functional region in the brain of a patient may be effectively displayed or identified. Herein, it is in particular advantageous to be provided with several independent parameters (first analysis image value, second analysis image value, second image) , to perform a more accurate analysis. For examining an aneurysm in the head of a patient in particular a method according to an embodiment of the present invention may be applied.
According to an embodiment, the method further comprises stimulating the patient. The stimulating may thereby comprise requesting to perform a particular action, requesting to speak a word, requesting to move a hand or a finger or the like. The stimulating may be an acoustic stimulating, visual stimulating or may be causing pain. The stimulating may thereby be repeatedly performed, in equal or not equal intervals, wherein regular or irregular sections of a pause may be inserted.
According to an embodiment of the present invention, the method further comprises: displaying a region of the object in dependence of a result of comparing one of the first analysis value, the second analysis value, and the image value of the at least one second image with preset threshold values for pixels to which the region of the object is imaged. For example, functional regions in a brain of a patient may be identified by detecting that one of the first analysis value, the second analysis value and the image value of the at least one second image exceeds or falls below preset threshold values. This may indicate that a perfusion, a degree of oxygenation, a haemoglobin concentration or a ICG concentration or their temporal changes are above or below an expected value.
According to an embodiment of the present invention, the method further comprises injecting indocyanin green into a blood stream of a patient.
According to an embodiment of the present invention, a system for examining an object is provided, wherein the system comprises: a light source for illuminating the object with measuring light; a first detector for recording at least one sequence of first images of the illuminated object, wherein an exposure time for a recording of an image of the sequence of first images is adjustable to be smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms; an evaluation system for evaluating the at least one sequence of first images (k = 1,...,N), wherein the evaluation system is adapted to associate pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value depends on temporal changes of the image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (I]<(i,j)) in a frequency region above 1 kHz, in particular above 10 kHz, wherein the evaluation system is further adapted to associate pixels (i,j) in the sequence of first images (Ifc(i,j)) to which same locations (x,y) of the object are imaged with a second analysis image value, respectively, wherein the second analysis image value depends on temporal changes of the image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1, j)) in a frequency range below 100 Hz.
According to an embodiment of the present invention, the system further comprises a first filter for spectrally filtering light for recording images of the illuminated object by the first detector, wherein the first filter is configured such that a ratio of an intensity of light outside a measuring light wavelength range from 605 nm to 630 nm, 560 nm to 580 nm, or 794 nm to 814 nm emanating from the filter from a total intensity of light outside this wavelength range being incident onto the first filter amounts to less than 10 %.
Such a configured first filter allows to detect measuring light which is, besides the determination of a perfusion, suitable for determination of a degree of oxygenation of haemoglobin or for determination of a total concentration of haemoglobin. Disturbing light wavelengths are thus efficiently suppressed.
According to an embodiment of the present invention, a system for examining an object is provided, the system comprising: a light source for illuminating the object with measuring light comprising wavelengths in a range from 790 nm to 820 nm; a first detector for recording at least one sequence of first images of the illuminated object, wherein an exposure time for a recording of an image of the sequence of first images is adjustable to be smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms; an evaluation system for evaluating the at least one sequence of first images (k = 1,...,N), wherein the evaluation system is adapted to associate pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a first analysis image value, wherein the first analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (Ikfi,])), wherein the system further comprises: a second detector for detecting light having wavelengths in a range from 820 nm to 850 nm, in particular 830 nm to 840 nm, wherein the second detector is arranged for recording at least one second image of the illuminated object and is configured such that light in a first wavelength range between 820 nm and 850 nm having an initial intensity emanating from the illuminated object along a beam path in a direction towards the second detector is detected by the second detector as a first intensity and such that light in a second wavelength range between 790 nm and 820 nm having the initial intensity emanating from the illuminated object along the beam path in the direction towards the second detector is detected by the second detector as a second intensity, wherein the first intensity is ten times higher than the second intensity.
The system is suited to excite a fluorescence of indocyanin green (ICG). The second detector is configured to detect fluorescence light emitted by ICG, wherein excitation light in the wavelength range between 790 nm and 820 nm is detected with ten times lower sensitivity than fluorescence light in a wavelength range from 820 nm to 850 nm.
According to an embodiment of the present invention, a first filter is arranged in a beam path of the measuring light in a direction of travelling of the measuring light in front of the first detector, wherein a ratio of an intensity of light in a wavelength range between 820 nm and 850 nm emanating from the first filter from a total intensity of light outside the wavelength range being incident onto the first filter amounts to less than 10 %.
According to an embodiment of the present invention, the second detector comprises a second filter arranged in a beam path of the measuring light in a direction of travelling of the measuring light in front of the second detector, wherein a ratio of an intensity of light in a wavelength range between 820 nm and 790 nm emanating from the second filter from a total intensity of light in the wavelength range being incident onto the second filter amounts to less than 10 %.
According to an embodiment of the present invention, a system for examining an object is provided, the system comprising: a light source for illuminating the object with measuring light; a first detector for recording at least one sequence of first images of the illuminated object, wherein an exposure time for a recording of an image of the sequence of first images is adjustable to be smaller than 1 ms, in particular smaller than 0.01 ms, and further in particular 0.03 ms; an evaluation system for evaluating the at least one sequence of first images (k = 1,...,N), wherein the evaluation system is adapted to associate pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (Ik(i,j))f wherein the first detector is adapted to detect measuring light in a first predetermined wavelength range from 560 nm to 580 nm or from 605 nm to 630 nm or from 794 nm to 814 nm and wherein the first detector is configured such that light in the first predetermined wavelength range having an initial intensity emanating from the illuminated object along a beam path in a direction towards the first detector is detected by the first detector as a first intensity and such that light outside the first predetermined wavelength range having the initial intensity emanating from the illuminated object along the beam path in the direction towards the first detector is detected by the first detector as a second intensity, wherein the first intensity is ten times, in particular one hundred times higher than the second intensity.
The system is suited to determine perfusion data as well as data about a degree of oxygenation or a total content of haemoglobin concurrently for the same object region.
According to an embodiment of the present invention, the system further comprises a second detector for recording at least one second image of the illuminated object, wherein the second detector is adapted to detect measuring light in a second predetermined wavelength range from 560 nm to 580 nm or from 605 nm to 630 nm or from 794 nm to 814 nm and wherein the second detector is configured such that light in the second predetermined wavelength range having a further initial intensity emanating from the illuminated object along a beam path in a direction towards the second detector is detected as a further first intensity and such that light outside the second predetermined wavelength range having the further initial intensity emanating from the illuminated object along the beam path in the direction towards the second detector is detected by the second detector as a further second intensity, wherein the further first intensity is ten times, in particular one hundred times higher than the further second intensity and wherein the first predetermined wavelength range is different from the second predetermined wavelength range.
By this provision it is possible to determine perfusion data as well as also data about a degree of oxygenation and data about a total content of haemoglobin.
According to an embodiment of the present invention, the evaluation system is further adapted to associate pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a second analysis image value, respectively, wherein the second analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range below 100 Hz.
Thereby, temporal changes of a concentration of ICG, of a degree of oxygenation and a total content of haemoglobin may be determined.
According to an embodiment of the present invention, the first analysis image value respectively associated with pixels in the sequence of first images is dependent of temporal changes of image values of the pixels of the at least one sequence of first images in a frequency range above 1 kHz, in particular above 10 kHz.
According to an embodiment of the present invention, the first detector comprises a first filter arranged in a beam path of measuring light detected by the first detector, wherein a ratio of an intensity of light outside the first predetermined wavelength range emanating from the first filter from a intensity of light outside the first predetermined wavelength range being incident onto the filter amounts to less than 10 %. According to an embodiment of the present invention, the evaluation system comprises a low path filter to determine temporal changes of image values of the pixels in the at least one sequence of first images in a frequency range below 100 Hz.
According to an embodiment of the present invention, the evaluation system is adapted to record images having a temporal interval of smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.05 ms .
For example, providing a detector allowing to record images having a temporal interval of 0.05 ms, allows a detection of measuring light signals in a frequency range up to 10 kHz. Thus, the detected image values represent a light signal sampled with a sampling frequency of 20 kHz. A maximal velocity of a moving particle that may be determined using such a detector is 4 mm/s when utilizing measuring light having a wavelength of 800 nm. This velocity approximately corresponds to blood flow velocities in capillaries and micro-capillaries of a patient.
According to an embodiment of the present invention, the light source is adapted to generate light having wavelengths in the region from 500 nm to 850 nm.
According to an embodiment of the present invention, the light source is adapted to generate light having a coherence length, comprising a temporal coherence length, of greater than 2 mm.
According to an embodiment of the present invention, the system further comprises an illumination apparatus for illuminating the object with white light and an imaging apparatus for imaging the object illuminated with white light to a white light image of the object and further a display apparatus for displaying images, in particular in superposition.
According to an embodiment of the present invention, the system further comprises a storage system for storing at least one of first images of the illuminated object, first analysis image values, second analysis image values, the second image and the white light image.
According to an embodiment of the present invention, the system further comprises a display apparatus for displaying at least one of the recorded first images, the first analysis image values, the second analysis image values, the second image, and the white light image.
Brief Description of the Drawings
The forgoing as well as other advantageous features of the invention will be more apparent from the following detailed description of exemplary embodiments of the invention with reference to the accompanying drawings. It is noted that not all possible embodiments of the present invention necessarily exhibit each and every, or any, of the advantages identified herein.
Figure 1 illustrates an embodiment of a method for examining an object according to the invention;
Figure 2 illustrates an evaluation method of a method for examining an object according to an embodiment of the present invention;
Figure 3 illustrates a system for examining an object according to an embodiment of the present invention;
Figure 4 illustrates a further system for examining an object according to an embodiment of the present invention;
Figure 5 illustrates an absorption spectrum and a fluorescence emission spectrum of indocyanin green (ICG) ; Figure 6 illustrates a system for examining an object according to an embodiment of the present invention;
Figure 7a illustrates a white light image and
Figure 7b illustrates an image of first analysis values (perfusion data) of an object, wherein the images were determined according to a method for examining an object according to an embodiment of the present invention;
Figure 8a, Figure 8b, Figure 8c, Figure 8d, Figure 8e, Figure 8f show second images determined by detection of fluorescence of ICG, which were determined according to a method according to an embodiment of the present invention;
Figure 9 illustrates a system for examining an object according to an embodiment of the present invention; and
Figure 10a, Figure 10b, Figure 10c, Figure 1Od illustrate a method according to an embodiment of the present invention for examining an object.
Detailed Description of Exemplary Embodiments
In the exemplary embodiments described below, components that are alike in function and structure are designated as far as possible by alike reference numerals. Therefore, to understand the features of the individual components of a specific embodiment, the descriptions of other embodiments and of the summary of the invention should be referred to.
Figure 1 shows a method for examining an object according to an embodiment of the present invention. A sequence of images Bi, ... , Bjc, ... , BN is recorded from an object illuminated with measuring light. A temporal interval t^ between the recording of the first image B^ and the recording of the second image B2 amounts to 0.05 ms. The sequence of images B]_, ... , B^, ... , BN is illustrated in Figure Ia. The image B^ is illustrated in Figure Ib. Every image B^ (k = 1,...,N) is described by image values 1^(1,J)), wherein (i,j) denote indices of the pixels of the detector arranged in an area. The image B^ has an image value Ik;(i,j)) at the location defined by the indices i and j. In the illustrated example the first index of the pixels of the image runs from 1 to L and the second index of the pixels of the image runs from 1 to M. Using the temporal distances between the recordings of two subsequent images for every pixel (i,j) of the sequence of images a temporal course of image values I(t,i,j) may be obtained. An example of such a temporal course of image values I(t,i,j) is illustrated in Figure Ic. In this Figure, the image value I(t,i,j) of a pixel (i,j) is plotted versus the time t. It is apparent that the temporal course I(t,i,j) comprises changes or fluctuations which are temporarily spaced apart less than 1 ms (high frequency portion of the temporal changes of the image values) which are superimposed by temporal changes which are temporarily spaced apart in the order of 100 ms or greater (low frequency portion of the temporal changes of image values). Temporal changes of image values being temporally spaced apart shorter than 0.5 ms are comprised in the image values but are not illustrated here. Temporal changes of image values being spaced apart longer than I s or 10 s are comprised in the temporal changes of the image values but can not be illustrated due to the chosen section of a temporal development of image values.
Figure 2 schematically illustrates an evaluation method according to an embodiment of the present invention. The left hand side of Figure 2a illustrates as a graph 11 temporal changes of images values in a frequency range above 1 kHz, in the following also denoted as high frequency range of temporal changes of the image values, or simply high frequency range. In particular this portion of temporal changes of image values carries information about a motion of particles within the illuminated object from which measuring light was reflected. A further evaluation of this high frequency portion of temporal changes of image values comprises establishing a Fourier spectrum or power spectrum of temporal changes of image values. The Fourier spectrum or the power spectrum may further be analysed with respect to its moments. Thereby, moments of this spectrum are for example an average value or an area. From the moments of the Fourier spectrum or power spectrum a concentration of particles moving relative to each other, a velocity of particles moving relative to each other and a perfusion of these particles may be obtained which parameter are denoted as first analysis values. When the afore-explained analysis of temporal changes of image values in the high frequency range is performed for every pixel of the sequence of images, B^, ... , B^, ... , BN an image 21 of first analysis values is obtained as illustrated in the right hand side of Figure 2a. The values of first analysis values within the image 21 of first analysis values are illustrated as contour lines 19 in the illustrated example. Figure 2b shows a further evaluation step of the inventive method. On the left hand side in Figure 2b the temporal course of image values I(t,i,j) of the temporal course of image values illustrated in Figure Ic is illustrated in a frequency range below 100 Hz as curve 12. Thus, the curve 12 does not comprise anymore the very fast temporal changes of image values occurring in time intervals below 1 ms, in particular 0.1 ms . Depending on the application and on the object to be examined the temporal changes of image values may be analysed in different low frequency ranges, such as below 1 Hz, below 10 Hz, or below 100 Hz. Thereby, the curve 12 may be determined in different manners. For example it is possible, to apply a low pass filter to the image value I(t,i,j) to obtain the curve 12. Alternatively, at least two images may be selected from the sequence of images B^, ... , B^, ... , BN to obtain the curve 12. This may for example also comprise interpolation methods to determine values of image values at intermediate locations not measured. From the values or temporal changes of the curve 12 an image 22 of second image values is obtained again showing contour lines 23. It is apparent that a shape of the contour lines 19 of the image 21 have a form different from a form of the contour lines 23 of the image 22 in Figure 2.
Figure 3 illustrates a system 1 for examining an object according to an embodiment of the present invention. The system 1 comprises an illumination apparatus 5, a beam splitter 24, an objective 26, a position acquisition apparatus 27, a first camera 14, a controller 42, and an evaluation system 40. The illumination apparatus 5 comprises a laser 6 for generating light having a coherence length of at least 2 mm. Preferably, the laser 6 is tuneable to be able to generate laser light having different wavelengths. The light 7 generated by the laser 6 is collimated by the collimator 8 to form measuring light 32. Measuring light 32 traverses the semi-transparent mirror 24 and the objective 26 to illuminate the object 28 in the object plane 30 with measuring light 32. A position (x,y) of the object 28 is registered by the position acquisition apparatus 27 thereby allowing a conversion of coordinates in an object coordinate system into coordinates of a coordinate system of the system 1. Thereby, it is in particular possible to convert coordinates in images of the object generated from the first camera 14 into coordinate of a coordinate system of the object. The light 34 emanating from a point (x,y) of the object 28 traverses the objective 26 and is reflected at the semi-transparent mirror 24 to enter the first camera 14. The first camera 14 comprises a first filter 16 to filter out undesired wavelengths from the light 34 emanating from the object or to filter out stray light from the surrounding. The filtered light 36 traverses the camera optics 18 imaging the filtered light 36 to a pixel (i,j) of the first detector 20. Different locations of the object 28 thus are imaged to different pixels of the first detector. The first detector 20 is capable to record images having an exposure time of 0.05 ms . A temporal interval between recordings of two images amounts to 0.1 ms. The first detector 20 and also the first filter 16 are connected with the camera controller 42. The camera controller 42 is capable to send signals to the first filter 16 to change a filter characteristics of the first filter 16. This may for example be realized in that a wheel having different filters arranged in a circumferential direction is rotated to bring a desired filter into a beam path of the first camera 14. Correspondingly, the controller 42 may control the tuneable laser to choose a desired wavelength of the measuring light 32. Depending on the object 28 to be examined measuring light having different wavelengths may be advantageous for examining the object. Wavelengths in a measuring light wavelength range from 605 nm to 630 ran, 560 nm to 580 ran or 794 nm to 814 nm may be utilized to measure a degree of oxygenation of haemoglobin or to measure a concentration of haemoglobin within the sample. Images recorded by the first detector are acquired by the controller 42 and transferred to the evaluation system 40. The evaluation system 40 comprises a processing system 44, a display device 46, a storage system 48 and an input system 50. The processing system 44 comprises a computer capable of running image processing software. The image processing software enables to read in a sequence of images generated by the first detector 20. The software is adapted to evaluate the sequence of images per pixel to analyse temporal changes of image values I(t,i,j). The image processing software is controlled via the input system 50 comprising a console and a cursor. On the display device 46 comprising a screen the course of the analysis of the sequence of images may be monitored. The evaluating the sequence of images per pixel comprises computing a spectrum of the temporal changes of image values and computing moments of this spectrum. Further, the temporal changes of image values are filtered with respect to their frequency range. Result images of the analysis may be displayed on the display device 46, in particular in superposition. Thereby, it is for example possible, to display perfusion data of the examined object together with degrees of oxygenation or total concentrations of haemoglobin, wherein those parameters are concurrently determined by the inventive system 1. Thus, a more accurate analysis of different aspects of the object is concurrently possible. These images and other data may also be supplied for display into a surgical microscope. Figure 4 shows a system 1 for examining an object according to an embodiment of the present invention. In contrast to the embodiment illustrated in Figure 3, the light 7 emitted by the laser 6 being collimated by the collimator 8 and leaving same as measuring light 32 does not travel along a common beam path together with the light 34 emanating from the object 28. Other elements of the system 1 illustrated in Figure 4 are analogous to those of the system 1 illustrated in Figure 3 and thus a repeated description is avoided.
Figure 5 illustrates an excitation spectrum and a fluorescence emission spectrum of indocyanin green (ICG) . Curve 60 shows an excitation spectrum of ICG, wherein a normalized absorption is plotted versus a wavelength. ICG absorbs light in a wavelength range from about 700 nm to 850 nm. Thereby, a maximum of absorption occurs at around 805 nm. By illuminating ICG with light between 790 nm and 820 nm a fluorescence emission of ICG is excited. Curve 62 shows a fluorescence emission spectrum of ICG, wherein a normalized fluorescence is plotted against a wavelength. It is apparent that ICG is fluorescent in a wavelength range from about 800 nm to 900 nm. A maximum of fluorescence emission is observed at around 835 nm. To detect a content of ICG within a sample using fluorescence of ICG it is thus necessary to illuminate the sample with light in a wavelength range from about 790 nm to 820 nm and to determine a fluorescence intensity by detection of light having wavelengths between 820 nm and around 850 nm emanating from the sample.
Figure 6 illustrates a system 1 for examining an object according to an embodiment of the present invention. The system 1 comprises an illumination apparatus 5, a semi-transparent mirror 24, a semi-transparent mirror 24', an objective 26, a first camera 14, a second camera 64 and an evaluation system 74. The illumination apparatus 5 comprises a laser 65 capable to generate light in a wavelength range from 790 nm to 820 nm. The light 7 emitted by the laser 65 is collimated by the collimator 8 to form measuring light 32. Measuring light 32 traverses the semi-transparent mirror 24 and the semi-transparent mirror 24' to traverse the objective 26 and illuminate the object 28. Light 34 emanating from the illuminated object traverses the objective 26 and impinges onto the semi-transparent mirror 24' at which a portion of the impinging light is reflected to enter the second camera 64. Another portion of the light impinging onto the semi- transparent mirror traverses the same and is reflected at the semi-transparent mirror 24 to enter the first camera 14. The first camera 14 comprises a first filter 16 having a transmission characteristics to allow traversing light having wavelengths suited for excitation of fluorescence of ICG, i.e. in a wavelength range between 790 nm and 820 nm, as filtered light 36. Filtered light 36 traverses the camera optics 18 to be imaged onto a pixel (i,j) of the first detector. The first detector 20 is capable to record images having an exposure time of 0.05 ms. A temporal interval between the recordings of two images amounts to 0.1 ms. The light entering the second camera 64 having emanated from the location (x,y) of the object 28 as light 34 traverses a second filter 66 which is designed to attenuate light used for excitation of fluorescence emission of ICG, i.e. light having wavelengths between 790 nm and 820 nm, upon traversal through the filter 66 strongly in its intensity. In particular, the filter 66 is designed such that only 10 % of an intensity of light in a wavelength range between 790 nm and 820 nm is allowed to traverse. On the other hand, the filter 66 is transparent for fluorescence light of ICG, i.e. transparent in a wavelength range between about 820 nm and 850 nm. The filtered light 86 emanating from the filter 66 traverses the camera optics 68 to be imaged onto a pixel (i,j) of the second detector 70. While, as described above, the first detector 20 is capable to subsequently record images having a very short temporal interval (0.1 ms), the second detector 70 records images having much larger temporal intervals (larger than 10 ms) . Thus, a CCD-detector may be employed as a second detector 70. The first detector 20 however may be a CMOS- detector. The controller 72 controls the recording of images of the first detector 20 as well as images of the second detector 70 and acquires the thus obtained first and second images. The first and second images are provided to the evaluation system 74, in particular to the processing system 76. The processing system 76 comprises as the processing system 44 of the embodiment illustrated in Figure 3 an image processing software. The image processing software is however also configured to process second images generated by the second detector 70. The processing the images acquired by the first detector 20 is simplified in comparison to the processing the first images of the embodiment illustrated in Figure 3. In the embodiment shown in Figure 6 not necessarily second analysis image values must be determined from a sequence of recorded first images, wherein the second analysis image values depend on temporal changes of image values in the sequence of first images in a frequency range below 100 Hz. However, as in the embodiment illustrated in Figure 3, first analysis values are determined from the temporal changes of image values of the sequence of recorded first images, wherein the first analysis values represent perfusion data. Thus, also in this embodiment, information about a motion of particles within the illuminated object are determined. Due to the choice of the measuring light being suited for excitation of a fluorescence of ICG, and the choice of the filter characteristics of the second filter 66, the second detector 70 substantially detects fluorescence light of ICG emanating from the object 28.
The embodiments illustrated in the Figures 3, 4 and 6 may further comprise a system for imaging the object as white light image. Such a white light imaging system may comprise an illumination source such as a halogen lamp, an objective, and a further detector such as a CCD camera. The white light imaging system may also be adapted as a microscope, in particular as a surgical microscope. Images generated by the white light imaging system may be provided to the evaluation system 40 and 74, respectively, and may be displayed together with the first and second images, the first and second analysis image values by the display device 46 or/and may be supplied for display in the surgical microscope.
Figure 7 shows images determined using systems according to embodiments of the present invention and using methods according to embodiments of the present invention, respectively. Figure 7a shows a white light image of a portion of a brain of a patient. A meshwork of capillaries and micro-capillaries is observed penetrating underlying brain tissue. Thereby, the bright red capillaries 90 and 92 are arteries and the dark red capillaries 93 represent veins. Reference sign 91 denotes an aneurysm. Figure 7b shows an image determined according to a method of the present invention showing a perfusion of blood cells in a region of the brain of the patient, wherein the region overlaps with the region of the brain illustrated in Figure 7a. Values of the perfusion are represented by pseudo-colour values (or grey values, respectively) , comprising blue, yellow, red in increasing order. In particular, it is apparent that the capillaries 90, 92 and 93 exhibit a lower perfusion than adjacent brain tissue. A possible explanation is that a motion of blood cells in the capillaries 90, 92 and 93 is that fast (such as greater than 10 mm/s) that such a motion cannot be resolved any more by the detector used here.
Figures 8a to 8f show images from a region of the brain overlapping with the regions illustrated in Figures 7a and 7b. The images shown in Figure 8 are acquired by detection of fluorescence light of ICG by the second detector 70 illustrated in Figure 6. The Figures 8a to 8f were recorded with a temporal interval of 0.5 s to I s and thus represent a temporal development of a fluorescence of ICG. Before the image 8a was recorded ICG was injected into the blood stream of the patient. After the injection of ICG initially (Figures 8a to 8d) an increasing intensity of fluorescence of ICG in the arteries, such as arteries 90 and 92, is observed. That means that ICG is supplied to the brain via the arteries. In contrast, in a region of the vein 93, practically no fluorescence of ICG is observed. The same holds for the aneurysm 91. Figure 8e shows a condition in which the brain is maximally filled with ICG. Subsequently outflow of ICG from the examined region via veins occurs. This is shown in Figure 8f. Thereby, vein 93 shows a high fluorescence intensity indicating that it guides blood. The same holds for the aneurysm 91 also indicating that the aneurysm is flown through by blood. A combination and evaluation of the images of the brain shown in the Figures 7 and 8 thus allows a detailed examination of the brain.
Figure 9 shows a system 100 for examining an object according to an embodiment of the present invention. The system comprises a laser 101 for generating laser light which illuminates the object 28 after collimating it by a collimator 105 and reflection at the semi-transparent mirror 131, after traversal through the semi-transparent mirror 130 and traversal through the objective 140. The light of the laser 101 reflected from the object 28 is detected by the detector 109 after traversal through the semi-transparent mirror 130 through the two semi- transparent mirrors 131, 132 and reflection at the semi- transparent mirror 134 and traversal through camera optics 107. Light source 111 emits light comprising wavelengths in the range from 790 nm to 820 nm and is thus suited for excitation of fluorescence of ICG. Light generated by the light source 111 is collimated by the collimator 115, is reflected at the semi- transparent mirror 132, traverses through the semi-transparent mirror 130 and through the objective 140 to illuminate the object 28. There it excites fluorescence of ICG. The fluorescence light of ICG emanating from the object 28 is detected by the detector 119 after traversal through several semi-transparent mirrors and traversal through the camera optics 117. For this the detector 119 comprises a suitable filter to attenuate light having wavelengths outside a range from 820 nm to 850 nm from light detected by the detector 119 strongly in its intensity. The system 100 further comprises a white light source 121 for illuminating the object 28 with white light. The white light image of the object 28 is detected by the detector 129 after traversal through the objective and several semi- transparent mirrors and traversal through the camera optics 127. Sequences of images recorded by the detector 109 and images recorded by the detectors 119 and 129 are supplied to a monitor and evaluation unit 132. The monitor and evaluation unit 132 is adapted to derive perfusion data from the sequence of images received by the detector 109 and to display these perfusion data together with the white light image and the ICG fluorescence image and to process these further. Figure 10 shows methods according to embodiments of the present invention. Figures 10a to 1Od show different illumination and detection scenarios according to embodiments of the present invention. In the Figures 10a to 1Od on the left hand side an illumination strength in dependence of a wavelength is illustrated, respectively, and in the Figures 10a to 1Od on the right hand side a detection sensitivity in dependence of the wavelength is illustrated, respectively. Depending on the embodiment, the illumination in different wavelength ranges may be performed by one or more light sources and the detection in different wavelength ranges may be performed by one or more detectors. With reference sign W in Figures 10a to 1Od on the left hand side and the right hand side a white light illumination wavelength range and a white light detection wavelength range is denoted, respectively. With reference sign ICG in these Figures an illumination wavelength range for illumination with light suitable for excitation of fluorescence of ICG and a ICG detection wavelength range suited for detection of fluorescence of ICG is denoted, respectively. Reference sign LDI in these Figures denotes an illumination wavelength range and a detection wavelength range suited to detect laser Doppler signals from an illuminated sample which are suitable to obtain perfusion data about the object, respectively. The scenario illustrated in Figure 10a may for example be realized using two light sources and three detectors. Hereby, the LDI illumination light is concurrently suited to excite fluorescence of ICG in the object. In embodiments, the LDI and/or ICG illumination light may be in a wavelength range overlapping or not overlapping with the illumination wavelength range of the white light illumination. An LDI detection wavelength range may in embodiments overlap with a detection wavelength range for white light detection or not. In the case of an overlap, the LDI detection signal may be separated from the white light detection signal by appropriate processing that may comprise an analysis of temporal changes of image values in a range above 1 kHz.
Summarized, embodiments of the present invention provide a system and a method for examining an object containing a fluid liquid, wherein the object is illuminated with measuring light and images are temporarily shortly subsequently recorded. The images are evaluated per pixel to determine perfusion data from a high frequency portion above 1 kHz and to determine further information about properties of the object from a low frequency portion below 100 Hz, such as a degree of oxygenation of haemoglobin, a concentration of haemoglobin or a concentration of ICG. This information determined by evaluation is displayed in a form of an image in superposition with a white light image of the object.
While the invention has been described with respect to certain exemplary embodiments thereof, it is evident that many alternatives, modifications and variations will be apparent to those skilled in the art. Accordingly, the exemplary embodiments of the invention set forth herein are intended to be illustrative and not limiting in any way. Various changes may be made without departing from the spirit and scope of the present invention as defined in the following claims.

Claims

What is claimed is :
1. A method for examining an object containing a flowing liquid, the method comprising:
illuminating the object with measuring light and recording at least one sequence of first images of the illuminated object, wherein an exposure time for each recording of an image of the sequence of first images is smaller than 1 ms, in particular smaller than 0.1 ms and further in particular smaller than 0.03 ms;
evaluating the at least one sequence of first images (k = 1,...,N) by associating pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range above 1 kHz, in particular above 10 kHz; and
evaluating the at least one sequence of first images
(k = 1,...,N), by further associating pixels (i,j) in the sequence of first images (Ijr(i,j)) to which same locations
(x,y) of the object are imaged with a second analysis image value, respectively, wherein the second analysis image value depends on temporal changes of image values
(I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range below
100 Hz.
2. The method according to claim 1, wherein the measuring light illuminating the object comprises wavelengths in a measuring light wavelength range from one of 605 nm to 630 nm, 560 nm to 580 nm, and 794 nm to 814 nm, and wherein the recording the first images comprises detecting measuring light, wherein a ratio of an intensity of the detected measuring light having wavelengths outside the measuring light wavelength range from a total intensity of the detected measuring light amounts to less than 10 %.
3. A method for examining an object containing a flowing liquid, the method comprising:
illuminating the object with measuring light comprising wavelengths in a range from 790 nm to 820 nm;
recording at least one sequence of first images of the illuminated object by a first detector, wherein an exposure time for each recording of an image of the sequence of first images is smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms;
evaluating the at least one sequence of first images (k = 1,...,N) by associating pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)); and
recording at least one second image of the illuminated object, wherein the second image is formed by detecting light having wavelengths in a range from 820 nm to 850 nm, in particular 830 nm to 840 nm, by a second detector and wherein a ratio of an intensity of light detected by the second detector having wavelengths smaller than 820 nm from a total intensity of light detected by the second detector amounts to less than 10 %.
4. The method according to claim 3, wherein the recording the at least one sequence of first images by the first detector comprises detecting measuring light in a wavelength range from 790 nm to 820 nm.
5. A method for examining an object containing a flowing liquid, the method comprising:
illuminating the object with measuring light and recording at least one sequence of first images of the illuminated object, wherein an exposure time for each recording of an image of the sequence of first images is smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms; and
evaluating the at least one sequence of first images (k = 1,...,N) by associating pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (Ik(i, j ) ) ;
wherein the recording the at least one sequence of first images of the illuminated object comprises detecting measuring light in a first predetermined wavelength range from one of 560 nm to 580 nm, 605 nm to 630 nm, and 794 nm to 814 nm, wherein a ratio of an intensity of detected measuring light having wavelengths outside the first predetermined wavelength range from a total intensity of detected measuring light amounts to less than 10 %.
6. The method according to claim 5, further comprising:
recording at least one second image of the illuminated object by detecting measuring light in a second predetermined wavelength range from one of 560 nm to 580 nm, 605 nm to 630 nm, and 794 nm to 814 nm, wherein a ratio of an intensity of detected measuring light having wavelengths outside the second predetermined wavelength range from a total intensity of detected measuring light amounts to less than 10 %, wherein the first predetermined wavelength range is different from the second predetermined wavelength range.
7. The method according to one of claims 3 to 6, further comprising:
evaluating the at least one sequence of first images
(k = 1,...,N) by further associating pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a second analysis image value, respectively, wherein the second analysis image value depends on temporal changes of image values
(I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range below 100 Hz.
8. The method according to one of claims 3 to 7, wherein the first analysis image value respectively associated with pixels (i,j) in the sequence of first images (1^(1,J)) depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range above 1 kHz, in particular above 10 kHz.
9. The method according to one of claims 1 to 8, wherein the recording the at least one sequence of first images comprises recording temporally subsequent images with a temporal interval of smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.05 ms.
10. The method according to one of claims 1 to 9, wherein the recording the first images comprises detecting measuring light, wherein a ratio of an intensity of the detected measuring light having wavelengths greater than 850 nm and smaller than 500 nm from a total intensity of the detected measuring light amounts to less than 10 %.
11. The method according to one of claims 1 to 10, wherein the measuring light has a coherence length corresponding at least to twice a penetration depth of the measuring light into the object.
12. The method according to one of claims 1 to 11, further comprising:
illuminating the object with white light; and
imaging and detecting the object illuminated with white light to determine a white light image of the object.
13. The method according to one of claims 1 to 12, wherein the first analysis image value associated with a pixel represents a value comprises one of a perfusion, a concentration, an average velocity, a measure for a velocity distribution, in particular a standard deviation of the velocity distribution, of particles moving relative to each other, in particular blood cells, and a combination of the same.
14. The method according to one of claims 1 to 13, further comprising at least one of displaying, in particular in superposition, and storing at least one of the first images, the first analysis image values, the second analysis image values, and the white light image for a plurality of pixels.
15. The method according to one of claims 1 to 14, further comprising repeating the recording the sequence of first images as temporal sequence and respective evaluating the repeatedly recorded sequence of first images.
16. The method according to one of claims 1 to 15, wherein the object is part of a brain of a patient, in particular an aneurysm in the head of the patient.
17. The method according to one of claims 1 to 16, further comprising:
stimulating the patient.
18. The method according to one of claims 1 to 17, further comprising:
displaying a region of the object in dependence of a result of comparing of at least one of the first analysis value, the second analysis value, and the image value of the at least one second image to preset threshold values for pixels to which the region of the object is imaged.
19. A system for examining an object, the system comprising:
a light source for illuminating the object with measuring light;
a first detector for recording at least one sequence of first images of the illuminated object, wherein an exposure time for a recording an image of the sequence of first images is adjustable to be smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms;
an evaluation system for evaluating the at least one sequence of first images ( k = 1,...,N), wherein the evaluation system is adapted,
to associate pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range above 1 kHz, in particular above 10 kHz, and to associate pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a second analysis image value, respectively, wherein the second analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixel (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range below 100 Hz.
20. The system according to claim 19, further comprising a first filter for spectrally filtering light for recording images of the illuminated object by the first detector, wherein the first filter is configured such that a ratio of an intensity of light emanating from the filter outside a measuring light wavelength range from one of 605 nm to 630 nm, 560 nm to 580 nm, and 794 nm to 814 nm from a total intensity of light being incident onto the first filter outside of this wavelength range amounts to less than 10 %.
21. A system for examining an object, the system comprising:
a light source for illuminating the object with measuring light comprising wavelengths in a range from 790 nm to 820 nm;
a first detector for recording at least one sequence of first images of the illuminated object, wherein an exposure time for recording an image of the sequence of first images is adjustable to be smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms;
an evaluation system for evaluating the at least one sequence of first images (k = 1,...,N), wherein the evaluation system is adapted,
to associate pixels (i,j) in the sequence of first images
(I]<(i,j)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value depends on temporal changes of the image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J));
wherein the system further comprises:
a second detector for detecting light having wavelengths in a range from 820 nm to 850 nm, in particular from 830 nm to 840 nm, wherein the second detector is arranged for recording at least one second image of the illuminated object and is configured such that
light in a first wavelength range between 820 nm and 850 nm having an initial intensity emanating from the illuminated object along a beam path in a direction towards the second detector is detected by the second detector as a first intensity, and
light in a second wavelength range between 790 nm and 820 nm having the initial intensity emanating from the illuminated object along the beam path in the direction towards the second detector is detected by the second detector as a second intensity, wherein the first intensity is ten times higher than the second intensity.
22. The system according to claim 21, wherein the first detector detects light in a wavelength range from 790 nm to 820 nm, and wherein a first filter is arranged in a beam path of the measuring light in a direction of travelling of the measuring light in front of the first detector, wherein a ratio of an intensity of light having wavelengths between 820 nm and 850 nm emanating from the first filter from a total intensity of light outside the wavelength range being incident onto the first filter amounts to less than 10 %.
23. The system according to claim 21 or 22, wherein the second detector comprises a second filter arranged in a beam path of the measuring light in a direction of travelling of the measuring light in front of the second detector, wherein a ratio of an intensity of light in a wavelength range between 820 nm and 790 nm emanating from the second filter from a total intensity of light in the wavelength range being incident onto the second filter amounts to less than 10 %.
24. A system for examining an object, the system comprising:
a light source for illuminating the object with measuring light;
a first detector for recording at least one sequence of first images of the illuminated object, wherein an exposure time for recording an image of the sequence of first images is adjustable to be smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.03 ms;
an evaluation system for evaluating the at least one sequence of first images (k = 1,...,N), wherein the evaluation system is adapted,
to associate pixels (i,j) in the sequence of first images (1^(1,J)) to which same locations (x,y) of the object are imaged with a first analysis image value, respectively, wherein the first analysis image value depends on temporal changes of the image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J));
wherein the first detector is adapted,
to detect measuring light in a first predetermined wavelength range from one of 560 nm to 580 nm, 605 nm to 630 nm, and 794 nm to 814 nm,
and the first detector is configured such that
light in the first predetermined wavelength range having an initial intensity emanating from the illuminated object along a beam path in a direction towards the first detector is detected by the first detector as a first intensity and
light outside the first predetermined wavelength range having the initial intensity emanating from the illuminated object along the beam path in the direction towards the first detector is detected by the first detector as a second intensity, wherein the first intensity is ten times higher, in particular one hundred times, than the second intensity.
25. The system according to claim 24, further comprising:
a second detector for recording at least one second image of the illuminated object, wherein the second detector is adapted to detect measuring light in a second predetermined wavelength range from one of 560 nm to 580 nm, 605 nm to 630 nm, and 794 nm to 814 nm and wherein the second detector is configured such that
light in the second predetermined wavelength range having a further initial intensity emanating from the illuminated object along a beam path in a direction towards the second detector is detected by the second detector as a further first intensity and
light outside the second predetermined wavelength range having the further initial intensity emanating from the illuminated object along the beam path in the direction towards the second detector is detected by the second detector as a further second intensity, wherein the further first intensity is ten times, in particular one hundred times, higher than the further second intensity and wherein the first predetermined wavelength range is different from the second predetermined wavelength range.
26. The system according to one of claims 21 to 25, wherein the evaluation system is further adapted, to associate pixels (i,j) in the sequence of first images
(Ij4 (i,j)) to which same locations (x,y) of the object are imaged with a second analysis image value, respectively, wherein the second analysis image value depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (1^(1,J)) in a frequency range below 100 Hz.
27. The system according to one of claims 21 to 26, wherein the first analysis image value associated to respective pixels
(i,j) in the sequence of first images (I]4 (i,j)) depends on temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (I]4(IfJ)) in a frequency region above 1 kHz, in particular above 10 kHz.
28. The system according to one of claims 24 to 27, wherein the first detector comprises a first filter arranged in a beam path of the measuring light detected by the first detector, wherein a ratio of an intensity of light outside the first predetermined wavelength range emanating from the first filter from an intensity of light outside the first predetermined wavelength range being incident onto the first filter amounts to less than 10 %.
29. The system according to claim 19 or 20 or, in combination with claim 26, according to one of claims 21 to 28, wherein the evaluation system comprises a low path filter to determine temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (I]4(^j)) in a frequency range below 100 Hz.
30. The system according to claim 19 or 20 or, in combination with claim 26, according to one of claims 21 to 29, wherein the evaluation system is adapted to determine temporal changes of image values (I(t,i,j)) of the pixels (i,j) in the at least one sequence of first images (I]4(^j)) in a frequency range below 100 Hz by selecting at least two images from the at least one sequence of first images.
31. The system according to one of claims 19 to 30, wherein the first detector is adapted to record images having a temporal interval of smaller than 1 ms, in particular smaller than 0.1 ms, and further in particular smaller than 0.05 ms.
32. The system according to one of claims 19 to 31, wherein the light source is adapted to generate light comprising wavelengths in a range from 500 nm to 850 nm.
33. The system according to one of claims 19 to 32, wherein the light source is adapted to generate light having a coherence length of greater than 2 mm.
34. The system according to one of claims 19 to 33, further comprising an illumination apparatus for illuminating the object with white light, and an imaging apparatus for imaging the object illuminated with white light to a white light image of the object.
35. The system according to one of claims 19 to 34, further comprising a storage system for storing at least one of first images of the illuminated object, first analysis image values, second analysis image values, the second image, and the white light image.
36. The system according to one of claims 19 to 35, further comprising a display apparatus for displaying at least one of the recorded first images, the first analysis image values, the second analysis image values, the second image, and the white light image, in particular in superposition.
PCT/EP2008/008547 2007-10-09 2008-10-09 System and method for examining an object WO2009046985A2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US12/757,000 US8929974B2 (en) 2007-10-09 2010-04-08 System and method for examining an illuminated object

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102007048362.9 2007-10-09
DE102007048362.9A DE102007048362B9 (en) 2007-10-09 2007-10-09 System and method for examining an object

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US12/757,000 Continuation-In-Part US8929974B2 (en) 2007-10-09 2010-04-08 System and method for examining an illuminated object

Publications (2)

Publication Number Publication Date
WO2009046985A2 true WO2009046985A2 (en) 2009-04-16
WO2009046985A3 WO2009046985A3 (en) 2009-07-09

Family

ID=40151643

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2008/008547 WO2009046985A2 (en) 2007-10-09 2008-10-09 System and method for examining an object

Country Status (3)

Country Link
US (1) US8929974B2 (en)
DE (1) DE102007048362B9 (en)
WO (1) WO2009046985A2 (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2606817A1 (en) * 2011-12-23 2013-06-26 Commissariat à l'Énergie Atomique et aux Énergies Alternatives System for reconstructing optical properties of a scattering medium, including a source of pulsed radiation and at least two detectors of two different types, and associated reconstruction method
EP2863801A4 (en) * 2012-06-21 2015-11-11 Novadaq Technologies Inc Quantification and analysis of angiography and perfusion
US9610021B2 (en) 2008-01-25 2017-04-04 Novadaq Technologies Inc. Method for evaluating blush in myocardial tissue
US9816930B2 (en) 2014-09-29 2017-11-14 Novadaq Technologies Inc. Imaging a target fluorophore in a biological material in the presence of autofluorescence
US10041042B2 (en) 2008-05-02 2018-08-07 Novadaq Technologies ULC Methods for production and use of substance-loaded erythrocytes (S-IEs) for observation and treatment of microvascular hemodynamics
US10219742B2 (en) 2008-04-14 2019-03-05 Novadaq Technologies ULC Locating and analyzing perforator flaps for plastic and reconstructive surgery
US10265419B2 (en) 2005-09-02 2019-04-23 Novadaq Technologies ULC Intraoperative determination of nerve location
US10434190B2 (en) 2006-09-07 2019-10-08 Novadaq Technologies ULC Pre-and-intra-operative localization of penile sentinel nodes
US10492671B2 (en) 2009-05-08 2019-12-03 Novadaq Technologies ULC Near infra red fluorescence imaging for visualization of blood vessels during endoscopic harvest
US10631746B2 (en) 2014-10-09 2020-04-28 Novadaq Technologies ULC Quantification of absolute blood flow in tissue using fluorescence-mediated photoplethysmography
US10992848B2 (en) 2017-02-10 2021-04-27 Novadaq Technologies ULC Open-field handheld fluorescence imaging systems and methods

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8224021B2 (en) * 2008-03-14 2012-07-17 Millivision Technologies, Inc. Method and system for automatic detection of a class of objects
US9066686B2 (en) 2008-07-10 2015-06-30 Novadaq Technologies Inc. Functional optical coherent imaging
DE102008040802B4 (en) * 2008-07-28 2014-09-18 Carl Zeiss Meditec Ag Method for the quantitative representation of blood flow as well as surgical microscope and analysis system
DE102011003606A1 (en) * 2011-02-03 2012-08-09 Siemens Aktiengesellschaft Method for removing unwanted perfusion-related signals from MR images, magnetic resonance system, computer program product and electronically readable data carriers
WO2013160861A1 (en) * 2012-04-27 2013-10-31 Aïmago S.A. Optical coherent imaging medical device
JP6023883B2 (en) * 2012-07-10 2016-11-09 アイマゴ ソシエテ アノニムAimago S.A. Perfusion assessment multi-modality optical medical device
DE102012019469A1 (en) * 2012-09-28 2014-04-03 Carl Zeiss Meditec Ag Method for the realization of OCT and other image recordings of an eye
DE102014103044A1 (en) * 2014-03-07 2015-09-10 Carl Zeiss Meditec Ag Surgical assistance system
CN107735015B (en) 2015-05-07 2021-09-10 史赛克欧洲运营有限公司 Method and system for laser speckle imaging of tissue using a color image sensor
CA3102138A1 (en) * 2018-06-08 2019-12-12 East Carolina University Determining peripheral oxygen saturation (spo2) and hemoglobin concentration using multi-spectral laser imaging (msli) methods and systems
CA3110871A1 (en) * 2018-09-05 2020-04-23 East Carolina University Systems for detecting vascular and arterial disease in asymptomatic patients and related methods
US11647889B2 (en) 2019-03-26 2023-05-16 East Carolina University Near-infrared fluorescence imaging for blood flow and perfusion visualization and related systems and computer program products
US10878568B1 (en) * 2019-08-08 2020-12-29 Neuroptica, Llc Systems and methods for imaging disease biomarkers

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050187477A1 (en) * 2002-02-01 2005-08-25 Serov Alexander N. Laser doppler perfusion imaging with a plurality of beams
US20050254008A1 (en) * 2002-06-14 2005-11-17 Ferguson R D Monitoring blood flow in the retina using a line-scanning laser ophthalmoscope
WO2006111836A1 (en) * 2005-04-20 2006-10-26 Ecole Polytechnique Federale De Lausanne (Epfl) Instrument and method for high-speed perfusion imaging

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4109647A (en) 1977-03-16 1978-08-29 The United States Of America As Represented By The Secretary Of The Department Of Health, Education And Welfare Method of and apparatus for measurement of blood flow using coherent light
SE468925B (en) 1991-08-22 1993-04-19 Gert Nilsson A METHOD AND APPARATUS SHOULD REDUCE THE DISTANCE-BASED FACTOR IN Saturation of STRAIGHT MOVEMENTS WITH AN IMAGING LASER-DOUBLE TECHNIQUE, SPECIFICALLY IN SEATING BLOOD PERFUSION THROUGH
GB9610700D0 (en) 1996-05-22 1996-07-31 Moor Instr Ltd Apparatus for imaging microvascular blood flow
AT409451B (en) 1999-12-14 2002-08-26 Hoffmann La Roche DEVICE FOR DETERMINING THE LOCAL DISTRIBUTION OF A MEASURED VALUE
US7580185B2 (en) * 2002-08-28 2009-08-25 Carl Zeiss Surgical Gmbh Microscopy system, microscopy method and a method of treating an aneurysm
US20130296715A1 (en) * 2005-04-20 2013-11-07 Ecole Polytechnique Federale De Lausanne (Epfl) Instrument and method for high-speed perfusion imaging
US7450243B2 (en) * 2006-07-10 2008-11-11 The Board Of Trustees Of The University Of Illinois Volumetric endoscopic coherence microscopy using a coherent fiber bundle
DE102007018641B8 (en) * 2007-04-19 2009-10-08 Carl Zeiss Surgical Gmbh Navigation system for brain surgery
CN101784227B (en) * 2007-07-06 2013-12-04 工业研究有限公司 Laser speckle imaging systems and methods

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050187477A1 (en) * 2002-02-01 2005-08-25 Serov Alexander N. Laser doppler perfusion imaging with a plurality of beams
US20050254008A1 (en) * 2002-06-14 2005-11-17 Ferguson R D Monitoring blood flow in the retina using a line-scanning laser ophthalmoscope
WO2006111836A1 (en) * 2005-04-20 2006-10-26 Ecole Polytechnique Federale De Lausanne (Epfl) Instrument and method for high-speed perfusion imaging

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DRAIJER M J; HONDEBRINK E; STEENBERGEN W; VAN LEEUWEN T G: "Laser Doppler perfusion imaging with a high-speed CMOS-camera" PROCEEDINGS OF THE SPIE - THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING, vol. 6631, 17 June 2007 (2007-06-17), - 19 June 2007 (2007-06-19) pages 66310N-1-66310N-7, XP002516360 Novel Optical Instrumentation for Biomedical Applications III, Munich, Germany *
SEROV A ET AL: "LASER DOPPLER PERFUSION IMAGING WITH A COMPLIMENTARY METAL OXIDE SEMICONDUCTOR IMAGE SENSOR" OPTICS LETTERS, OSA, OPTICAL SOCIETY OF AMERICA, WASHINGTON, DC, US, vol. 27, no. 5, 1 March 2002 (2002-03-01), pages 300-302, XP001077767 ISSN: 0146-9592 *
SEROV A N ET AL: "Quasiparallel laser-Doppler perfusion imaging using a CMOS image sensor" PROCEEDINGS OF THE SPIE - THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING SPIE-INT. SOC. OPT. ENG USA, vol. 5067, no. 1, 2003, pages 73-84, XP002516361 ISSN: 0277-786X *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10265419B2 (en) 2005-09-02 2019-04-23 Novadaq Technologies ULC Intraoperative determination of nerve location
US10434190B2 (en) 2006-09-07 2019-10-08 Novadaq Technologies ULC Pre-and-intra-operative localization of penile sentinel nodes
US11564583B2 (en) 2008-01-25 2023-01-31 Stryker European Operations Limited Method for evaluating blush in myocardial tissue
US9610021B2 (en) 2008-01-25 2017-04-04 Novadaq Technologies Inc. Method for evaluating blush in myocardial tissue
US10835138B2 (en) 2008-01-25 2020-11-17 Stryker European Operations Limited Method for evaluating blush in myocardial tissue
US9936887B2 (en) 2008-01-25 2018-04-10 Novadaq Technologies ULC Method for evaluating blush in myocardial tissue
US10219742B2 (en) 2008-04-14 2019-03-05 Novadaq Technologies ULC Locating and analyzing perforator flaps for plastic and reconstructive surgery
US10041042B2 (en) 2008-05-02 2018-08-07 Novadaq Technologies ULC Methods for production and use of substance-loaded erythrocytes (S-IEs) for observation and treatment of microvascular hemodynamics
US10492671B2 (en) 2009-05-08 2019-12-03 Novadaq Technologies ULC Near infra red fluorescence imaging for visualization of blood vessels during endoscopic harvest
FR2985023A1 (en) * 2011-12-23 2013-06-28 Commissariat Energie Atomique SYSTEM FOR RECONSTRUCTING OPTICAL PROPERTIES OF A DIFFUSING MEDIUM, COMPRISING A PULSE RADIATION SOURCE AND AT LEAST TWO DETECTORS OF TWO DIFFERENT TYPES, AND ASSOCIATED RECONSTRUCTION METHOD
US9407796B2 (en) 2011-12-23 2016-08-02 Commissariat á l'énergie atomique et aux énergies alternatives System for reconstructing optical properties in a diffusing medium, comprising a pulsed radiation source and at least two detectors of two different types, and associated reconstruction method
EP2606817A1 (en) * 2011-12-23 2013-06-26 Commissariat à l'Énergie Atomique et aux Énergies Alternatives System for reconstructing optical properties of a scattering medium, including a source of pulsed radiation and at least two detectors of two different types, and associated reconstruction method
US10278585B2 (en) 2012-06-21 2019-05-07 Novadaq Technologies ULC Quantification and analysis of angiography and perfusion
US11284801B2 (en) 2012-06-21 2022-03-29 Stryker European Operations Limited Quantification and analysis of angiography and perfusion
EP2863801A4 (en) * 2012-06-21 2015-11-11 Novadaq Technologies Inc Quantification and analysis of angiography and perfusion
US10488340B2 (en) 2014-09-29 2019-11-26 Novadaq Technologies ULC Imaging a target fluorophore in a biological material in the presence of autofluorescence
US9816930B2 (en) 2014-09-29 2017-11-14 Novadaq Technologies Inc. Imaging a target fluorophore in a biological material in the presence of autofluorescence
US10631746B2 (en) 2014-10-09 2020-04-28 Novadaq Technologies ULC Quantification of absolute blood flow in tissue using fluorescence-mediated photoplethysmography
US10992848B2 (en) 2017-02-10 2021-04-27 Novadaq Technologies ULC Open-field handheld fluorescence imaging systems and methods
US11140305B2 (en) 2017-02-10 2021-10-05 Stryker European Operations Limited Open-field handheld fluorescence imaging systems and methods
US12028600B2 (en) 2017-02-10 2024-07-02 Stryker Corporation Open-field handheld fluorescence imaging systems and methods

Also Published As

Publication number Publication date
DE102007048362A8 (en) 2009-08-13
DE102007048362B9 (en) 2017-08-17
DE102007048362A1 (en) 2009-04-23
US8929974B2 (en) 2015-01-06
DE102007048362B4 (en) 2017-04-27
US20110090325A1 (en) 2011-04-21
WO2009046985A3 (en) 2009-07-09

Similar Documents

Publication Publication Date Title
US8929974B2 (en) System and method for examining an illuminated object
Srinivasan et al. OCT methods for capillary velocimetry
US8406847B2 (en) Biological observation apparatus and method
CA2824134C (en) Methods, systems and computer program products for noninvasive determination of blood flow distribution using speckle imaging techniques and hemodynamic modeling
Sakadžić et al. Large arteriolar component of oxygen delivery implies a safe margin of oxygen supply to cerebral tissue
EP2967380B1 (en) Methods, systems and computer program products for non-invasive determination of blood flow distribution using speckle imaging techniques and hemodynamic modeling
CN107595250B (en) Blood flow imaging method and system based on motion and graph mixed contrast
US20090054788A1 (en) Method and apparatus for displaying a field of a brain of a patient and navigation system for brain surgery
JP5075116B2 (en) Spectroscopic determination of analyte concentration
US20060147897A1 (en) Characterization of arteriosclerosis by optical imaging
CN107773217B (en) Living tissue microcirculation metabolism dynamic measuring device and method
Golub et al. Erythrocyte-associated transients in PO2 revealed in capillaries of rat mesentery
JP2006042955A (en) In vivo material optometric device
Liu et al. Quantitative assessment of optical coherence tomography angiography algorithms for neuroimaging
US11937951B2 (en) Method and apparatus for medical imaging using differencing of multiple fluorophores
EP4262545B1 (en) Device, method and systems for providing imaging of one or more aspects of blood perfusion
JPH07120384A (en) Method and apparatus for optical measurement
JP3825459B2 (en) Biological light measurement device
Kuo et al. The use of intensity‐based Doppler variance method for single vessel response to functional neurovascular activation
JPWO2015037055A1 (en) Fluorescence image acquisition device
CN112022085B (en) Method for calculating blood vessel flow in retina
US10682053B2 (en) Ophthalmic photothermal optical coherence tomography apparatus
Du et al. Optical Techniques to Study Drug-Induced Neurovascular and Cellular Changes In Vivo
Zherebtsov et al. Enhancing Transcranial Blood Flow Visualization with Dynamic Light Scattering Technologies: Advances in Quantitative Analysis
JP2013192834A (en) Biological light measuring apparatus

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 08838041

Country of ref document: EP

Kind code of ref document: A2

NENP Non-entry into the national phase in:

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 08838041

Country of ref document: EP

Kind code of ref document: A2