[go: up one dir, main page]

WO2009081883A1 - レーザドップラー血流測定方法及び装置 - Google Patents

レーザドップラー血流測定方法及び装置 Download PDF

Info

Publication number
WO2009081883A1
WO2009081883A1 PCT/JP2008/073245 JP2008073245W WO2009081883A1 WO 2009081883 A1 WO2009081883 A1 WO 2009081883A1 JP 2008073245 W JP2008073245 W JP 2008073245W WO 2009081883 A1 WO2009081883 A1 WO 2009081883A1
Authority
WO
WIPO (PCT)
Prior art keywords
blood flow
light
laser
blood
optical system
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/JP2008/073245
Other languages
English (en)
French (fr)
Inventor
Tadashi Hachiga
Hiroki Ishida
Shunsuke Akiguchi
Hiroki Shirakawa
Tsugunobu Andoh
Yasushi Kuraishi
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of National Colleges of Technologies Japan
University of Toyama NUC
Original Assignee
Institute of National Colleges of Technologies Japan
University of Toyama NUC
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 Institute of National Colleges of Technologies Japan, University of Toyama NUC filed Critical Institute of National Colleges of Technologies Japan
Priority to US12/809,738 priority Critical patent/US9031640B2/en
Priority to JP2009547086A priority patent/JP5234470B2/ja
Publication of WO2009081883A1 publication Critical patent/WO2009081883A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/026Measuring blood flow
    • A61B5/0285Measuring or recording phase velocity of blood waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/41Detecting, measuring or recording for evaluating the immune or lymphatic systems
    • A61B5/411Detecting or monitoring allergy or intolerance reactions to an allergenic agent or substance
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/026Measuring blood flow
    • A61B5/0261Measuring blood flow using optical means, e.g. infrared light

Definitions

  • the present invention relates to a blood flow measurement method for measuring a blood flow state of a living body by a laser Doppler method and a blood flow measurement device used for the measurement.
  • Data on human skin blood flow relate to the state of metabolic activity of skin tissue, the state of metabolism, the degree of aging, skin cancer, and various other skin conditions. It is important in judging the state of
  • the first cause of death in Japanese is malignant neoplasm
  • the second is heart disease
  • the third is cerebrovascular disease
  • the death from vascular lesions is almost the same rate as malignant neoplasm.
  • three-dimensional mapping information of surrounding new blood vessels and blood flow information in the new blood vessels are important.
  • vascular lesions have frequent occurrence sites, and spatiotemporal information such as blood flow velocity distribution in blood vessels is indispensable.
  • Recent blood flow measurement methods include laser Doppler flow velocity measurement method, particle image flow velocity measurement method (PIV), micro PIV method, ultrasonic blood flow meter method, MRI and CT method, and a method by fusing with numerical simulation. is there.
  • Experiments on living organisms include experiments in the body (In vivo) and experiments in the body (In vitro).
  • PIV measurement for obtaining two-dimensional flow information is in vitro measurement.
  • current MRI and CT performance has insufficient spatial and temporal resolution, it is difficult to obtain detailed hemodynamic information by in vivo measurement.
  • the laser Doppler measurement method is widely used as a blood flow measurement method because it can continuously measure instantaneous values of blood flow velocity and / or blood flow without invading the skin.
  • the blood flow measurement method using the laser Doppler method uses the monochromaticity and coherence of laser light. When laser light is irradiated to blood vessels in skin tissue, the light scattered by blood cells moving in the blood vessels causess a frequency shift due to the Doppler effect, and the blood flow rate is obtained based on the frequency shift amount.
  • a sensor element in which a pair of an optical fiber for transmitting laser light incident on the skin surface from a light source and an optical fiber for receiving light scattered from the skin are provided. Some of them are embedded in a probe support member facing the skin.
  • a laser irradiation method of this laser Doppler method a single beam method in which one beam is incident on a subject and a single beam are divided into two parallel lights and intersected by a lens or the like to measure an intersection (measurement). There is a two-beam method in which a point) is incident on the subject.
  • a differential type method and a reference beam method as a two-beam type measurement method.
  • the reference light method one laser beam is divided into two strong and weak beams, and interference between light shifted in frequency by the Doppler effect and light scattered by a stationary tissue is used. Since both of these scattered lights have a frequency difference of about several hundred Hz to several tens of Hz, interference between the two is detected as an electrical beat signal by a light receiving element such as a photodiode, and the detected frequency
  • the intensity is a signal having a value corresponding to the speed and number of blood cells. Further, the values of these signals are integrated and converted into signals relating to the blood flow velocity and flow rate of the blood vessel to obtain the blood flow rate of the blood vessel.
  • one laser beam is divided into two beams, the two beams are condensed and intersected, and interference of scattered light occurs depending on the direction of laser light irradiation at the intersecting position. .
  • the interval between the interference fringes varies depending on the Doppler frequency shift amount of the scattered light by the blood cells as the observation particles, and this is observed to obtain the blood flow velocity.
  • the differential method has an advantage that the scattered light can be collected by a wide light receiving aperture and a relatively strong signal intensity can be obtained.
  • a sensor unit that includes a plurality of sensor elements that detect the surface blood flow of a subject and generate an output signal, and a predetermined output signal from the sensor element.
  • a surface blood flow measuring device including a signal conversion unit that converts the measurement signal into a measurement signal and a display unit that displays the measurement signal as a measurement value.
  • an output signal is emitted from each sensor element according to the surface blood flow of the subject, and is converted into a predetermined measurement signal by an output signal conversion unit.
  • a plurality of sensor elements are provided in the sensor unit, a plurality of locations can be measured simultaneously or selectively at one measurement site. Furthermore, by averaging the fluctuations of each measurement value, calibration becomes possible and quantitative measurement can be performed.
  • each light receiving unit receives scattered light from a subject and measures the blood flow in each light receiving unit, thereby simultaneously recording blood flows at different measurement depths.
  • the laser light is converted into sheet light by a cylindrical lens, is incident on a mirror disposed at the focal position of the lens via the lens, and is reflected in the deflection direction of the mirror.
  • a blood flow distribution measuring device that irradiates and measures a blood flow two-dimensionally.
  • Non-Patent Documents 1 to 5 a laser Doppler capable of simultaneous multipoint measurement using a laser beam sheet light as a device for measuring the flow velocity of fluid in a non-living pipe or the like An anemometer was proposed.
  • Japanese Patent Laid-Open No. 5-15501 JP-A-8-182658 Japanese Patent Application Laid-Open No. Hei 7-100119 Transactions of the Society of Instrument and Control Engineers Voi.39, No.3 (2003) 218-224 Transactions of the Japan Society of Mechanical Engineers (Part B), Vol. 69, No. 677 (2003-1) 25-30 Experiments in Fluids24 (1998) 70-76 Experiments in Fluids 36 (2004) 274-281 SICE 2002 Aug. 5-7, 2002, Osaka, 2199-2204
  • a conventional single-flux blood flow meter uses a single beam (about 1 mm in diameter) for incident laser light, and detects scattered light with a set of sensors.
  • Some blood flowmeters may detect scattered light using a plurality of optical fibers, or may include a plurality of optical fibers for incident laser light and optical signal detection. It can be said that it is multipoint measurement by the combination of measurement. Therefore, the blood flow volume in a wide area cannot be obtained efficiently, and the resolution is low, and only a limited blood flow volume can be obtained.
  • the blood flow imaging apparatus measures a blood flow distribution by scanning a laser beam.
  • a blood flow imaging device using a two-beam blood flow meter measures one point with two light beams, and obtains a blood flow distribution by scanning the laser beam.
  • the apparatus using the two-beam method can perform point measurement with high spatiotemporal resolution, it is difficult to measure multiple points at the same time. The error in speed is large, and in this case as well, a wide range of blood flow cannot be measured.
  • the laser power in each sensor element is lowered, so that only the blood flow rate at a relatively shallow part of the skin surface can be measured. Furthermore, it is difficult to measure multiple points at the same time, and when the position of the blood vessel is shifted, the error in the measured average blood flow velocity increases.
  • the measurement depth can be individually adjusted by changing the laser light output and the distance between the light transmitting optical fiber and the light receiving optical fiber on the skin contact surface of the sensor element, that is, the distance between the light transmitting part and the light receiving part. Although it can, it does not measure blood flow at various measurement depths simultaneously. Therefore, there is a problem that the blood flow distribution in the depth direction cannot be known, and the state analysis of the skin tissue cannot be performed accurately.
  • the blood flow measuring devices disclosed in Patent Documents 1 and 2 also use optical fibers, which are measurement methods based on single-point measurement and have the same problems as described above.
  • the blood flow measuring device disclosed in Patent Document 3 is also not capable of detecting a three-dimensional blood vessel arrangement.
  • Non-Patent Documents 1 to 5 know the position of the conduit and the direction of fluid flow in advance, and cannot detect the arrangement of capillaries in invisible places.
  • the present invention has been made in view of the above-described conventional background art, and a laser Doppler blood flow measurement method capable of performing multidimensional measurement efficiently with high measurement accuracy and a wide range by a simple optical system and apparatus, and An object is to provide an apparatus.
  • the present invention divides a laser beam from a laser light source, makes each branched laser beam a sheet-like thin sheet light, crosses each other at a predetermined position in a living body, and the linear irradiation site where the sheet light intersects
  • the Doppler frequency is formed by the blood flow of the blood vessel in the living body among the scattered light that is incident on a plurality of light receiving elements arranged linearly at the imaging position.
  • the shifted scattered light is converted into an electric signal for each of the photoelectric conversion elements, and the laser light of the plurality of light receiving elements is shifted by a Doppler frequency and incident on the blood vessel in the living body depending on the position of the light receiving element.
  • the position information is detected and the blood flow velocity in the blood vessel is calculated for each photoelectric conversion element by the Doppler frequency shift of the laser light due to the blood flow.
  • the intersecting position of the sheet-like laser light is scanned in a direction intersecting with the longitudinal direction of the linear scattered light, and the blood vessel is detected for each photoelectric conversion element from the scanned linear scattered light.
  • the blood flow inside is calculated
  • the crossing position of the sheet-like laser light is scanned in a direction crossing the longitudinal direction of the linear scattered light, and the intravascular blood vessel is scanned for each photoelectric conversion element from the scanned linear scattered light.
  • Blood flow velocity in each blood vessel or blood flow volume in a predetermined region in the living body is calculated based on the blood flow velocity at each of the obtained linear irradiation sites. is there.
  • the scanning is performed in two directions orthogonal to the longitudinal direction of the linear irradiation site, and the linear irradiation site is sequentially shifted in a direction parallel to the longitudinal direction, so that a three-dimensional blood vessel in a predetermined region in the living body is obtained. It seeks the placement and makes it displayable. In addition, a three-dimensional blood flow state of the living tissue is calculated and displayed by the scanning.
  • the present invention also includes a laser light source, a beam splitting device that splits the laser light from the laser light source, and a cylindrical optical system that converts the laser light split by the beam splitting device into a planar thin sheet light.
  • a condensing optical system that intersects each sheet light with each other at a predetermined position in the living body, an imaging optical system that linearly forms an image of scattered light at a linear irradiation site where the sheet light intersects, and the imaging
  • a plurality of light receiving elements arranged at an imaging position of the optical system, a photoelectric conversion element that converts the scattered light incident on the light receiving element into an electrical signal for each photoelectric conversion element, and blood flow in the living tissue.
  • Blood flow calculation means for calculating the blood flow velocity at the linear irradiation site where the laser beam is shifted by the Doppler frequency for each photoelectric conversion element based on the electrical signal, and the intersection position of the sheet-like laser beam
  • the linear A laser beam scanning unit that scans in a direction intersecting with the longitudinal direction of the projecting site, and calculates the blood vessel arrangement in the living body in the region scanned by the laser beam scanning unit based on the calculation by the blood flow calculation unit This is a laser Doppler blood flow measurement device.
  • each blood flow rate and each blood vessel arrangement for each photoelectric conversion element obtained by the blood flow calculation unit is used to calculate each in-vivo body.
  • a blood flow rate calculating means for calculating the blood flow rate of the blood vessel is provided.
  • the imaging optical system is located on the same side as the condensing optical system with respect to the living tissue to be measured, and the condensing optical system is located between the branched sheet lights. Alternatively, the imaging optical system may be located on the opposite side of the condensing optical system across the living tissue to be measured.
  • the light receiving element is a solid-state image pickup element in which the light receiving elements are integrally arranged on a plane.
  • the laser Doppler blood flow measuring method and apparatus of the present invention it is possible to directly measure the blood flow velocity in the human body with high-order spatial resolution. Further, in a relatively thick blood vessel, it is possible to measure the blood flow velocity distribution and its time variation and to measure the blood flow velocity in the capillary. By scanning with laser light, it is possible to measure the effective inner diameter of a thick blood vessel, blood flow velocity distribution and its time fluctuation, and the average blood flow velocity of capillaries in a predetermined measurement range. The measurement result can be three-dimensionally mapped within a predetermined measurement range. In addition, there is no effect on blood flow due to contact, and it is suitable for measurement of inflammatory site, ear capillaries, aorta close to skin, blood flow of veins and blood flow distribution in blood vessels, and highly reproducible measurement results. can get.
  • FIG. 1 It is the schematic which shows arrangement
  • mapping figure which shows the blood vessel arrangement
  • a laser Doppler blood flow measuring device 10 is a laser light source as shown in FIG.
  • a semiconductor laser 12 such as a laser diode, and acoustooptic elements 14 and 15 for branching a continuous wave laser beam from the semiconductor laser 12 into two beams having slightly different frequencies are provided.
  • the wavelength of the laser light is preferably about 750 to 1500 nm.
  • a near infrared laser having a wavelength of 785 nm is used. This is because near infrared light has high transmittance into living tissue.
  • the light emitting side of the acoustooptic device 15 includes a reflecting prism 16 that reflects the branched laser light in two directions, and reflecting mirrors 18 and 19 that reflect the laser light reflected by the reflecting prism 16 at a right angle. ing. Furthermore, rod lenses 20 and 21 provided at positions parallel to each other where the two laser beams reflected by the reflecting mirrors 18 and 19 enter are arranged. Each rod lens 20, 21 spreads the incident laser light and makes it incident on the flat portion 22 a of the first cylindrical lens 22.
  • the first cylindrical lens 22 is positioned so that the flat portion 22 a faces the rod lenses 20 and 21, and the two incident laser beams are emitted from the cylindrical surface-shaped curved surface portion 22 b of the first cylindrical lens 22.
  • the first cylindrical lens 22 is formed into two sheet-like sheet lights Ls that are parallel to the central axis of the cylinder and orthogonal to the plane portion 22a.
  • the sheet light Ls further enters the flat portion 24a of the second cylindrical lens 24 that forms the condensing optical system.
  • the flat surface portion 24 a of the second cylindrical lens 24 is positioned parallel to the same side as the flat surface portion 22 a of the first cylindrical lens 22, and is a cylindrical curved surface portion that is the other surface of the second cylindrical lens 24.
  • 24b is arranged so that the central axes of the cylinders are orthogonal to each other in a twisted positional relationship with respect to the first cylindrical lens 22.
  • the emitted light travels in the same optical path regardless of which of the plane portions 22a and 24a and the curved surface portions 22b and 24b is the incident surface.
  • the surface of the optical system component used in this embodiment is coated with an antireflection film for this wavelength in order to suppress the attenuation of the wavelength of the laser beam used.
  • the two sheet lights Ls incident on the second cylindrical lens 24 are refracted by the curved surface portion 24b and deflected so as to cross each other.
  • the intersecting portion of the sheet light Ls is formed as a portion having a high light intensity in one line, and becomes a linear irradiation site Lx of the living tissue.
  • Scattered light at the linear irradiation site Lx is converged by a lens system 30 that forms an imaging optical system as shown in FIG.
  • the lens system 30 forms an image of the scattered light from the linear irradiation site Lx in a line shape along the linear irradiation site Lx on the end surface portion 32a of the optical fiber array 32 that is a light receiving element.
  • each optical receptacle 40 is provided with the number of avalanche photodiodes 42 as photoelectric conversion elements, and converts the light guided by the optical fibers 34 into electric signals.
  • the output of each avalanche photodiode 42 is input to each frequency filter 46 through each amplifier 44 provided for the number of optical fibers 34, and a signal of a predetermined frequency is A / D converted, so that the number of optical fibers 34 is obtained.
  • a data recorder 48 provided only for each.
  • the digital signal recorded in each data recorder 48 is input to the computer 50 via the USB connector, and a predetermined calculation analysis process is performed.
  • the laser Doppler blood flow measurement method of this embodiment will be described below.
  • the laser light emitted from the semiconductor laser 12 is modulated by the acoustooptic elements 14 and 15 and branched, and as shown in FIGS. While making it thin sheet light Ls, it cross
  • FIG. In this state, the blood flow rate in the blood vessel 52 and, if necessary, the blood flow rate can be obtained from the blood flow rate with a spatial resolution determined by the magnification of the lens system 30 as an imaging optical system with the optical fiber 34 as a unit. it can.
  • the magnification of the optical system is 1 using the optical fiber 34 having a diameter of 250 ⁇ m
  • the spatial resolution is 250 ⁇ m.
  • the intersecting portion of the sheet light Ls is formed as a single line having a high light intensity and becomes a linear irradiation site Lx of the living tissue.
  • This linear portion Lx has interference fringes due to the phase difference of the incident wavefront of the laser light, and appears as interference fringes as shown in FIG.
  • the Doppler frequency of the scattered light changes.
  • the Doppler frequency of the scattered light changes due to the Doppler effect of the laser light scattered by the blood cells in the bloodstream irradiated with the sheet light Ls due to the flow velocity of the bloodstream. There is no change in the Doppler frequency of scattered light in a stationary tissue.
  • Each scattered light at the linear irradiation site Lx is imaged at a point corresponding to each on the optical fiber array 32 by the lens system 30.
  • the scattered light having the changed frequency is input to the avalanche photodiode 42 through the optical fiber 34 and detected as an electric beat signal.
  • the detected frequency change and intensity are values corresponding to the velocity and number of blood cells. . This is calculated
  • the input signal obtained by the avalanche photodiode 42 is data as shown in the graph (a) of FIG. 8, for example, and this data is A / D converted to a computer.
  • the input digital signal is subjected to fast Fourier transform and further subjected to noise removal processing and the like to calculate a Doppler frequency.
  • the Doppler frequency is obtained by the appearance of a peak at the Doppler frequency in the frequency spectrum of the detected signal.
  • the blood cell velocity is calculated from the frequency, and the blood flow velocity is calculated because the blood cell velocity is the blood flow velocity.
  • the modulation is applied by the acoustooptic device 14, the direction of blood flow can also be identified. However, when it is not necessary to identify the direction of flow, measurement is possible without frequency modulation by the acoustooptic device 14.
  • the linear irradiation site Lx is scanned in a direction orthogonal to the longitudinal direction thereof, and a flow rate is obtained for each optical fiber 34 at a predetermined sampling period, and a blood flow mapping within a predetermined plane is formed. The Further, this is scanned in another direction orthogonal to the longitudinal direction of the linear irradiation site Lx, and three-dimensional blood flow velocity or blood flow volume can be mapped.
  • the linear irradiation site Lx can be scanned by placing the living body on a table movable in the XYZ axis directions orthogonal to each other and moving in the three axis directions. It is.
  • scanning can be performed by driving an optical system such as the cylindrical lenses 22 and 24. At this time, the crossing position of the sheet light Ls and the scanning direction thereof can be easily obtained by calculation for obtaining the intersection of the straight lines.
  • a linear irradiation site Lx is formed in the tissue of the living body B, and the linear irradiation site Lx is orthogonal to the longitudinal direction.
  • the linear irradiation site Lx is further translated by the irradiation length in the Y-axis direction, and scanning is similarly repeated.
  • FIG. 9B a small distance is moved in the orthogonal Z-axis direction, scanning is performed in the same manner as described above, and this movement and scanning are repeated within a predetermined range.
  • the three-dimensional blood vessel arrangement in the predetermined range of the living body B and the blood flow velocity or blood flow volume of the blood vessel can be obtained.
  • the linear irradiation site Lx is parallel to the X-axis direction, and as shown in FIG. 9C, scanning is performed in the Y-axis direction and shifted in the X-axis direction. Then, it further moves in the Z-axis direction and repeats scanning.
  • the inner diameter of the blood vessel 52 and the blood flow state of the capillary can be measured and three-dimensionally mapped.
  • the spatial resolution in this case can be about several hundred microns. That is, since the inner diameter of the optical fiber 34 of the optical fiber lay 32 that receives the scattered light at the linear irradiation site Lx of the intersecting line of the sheet light Ls becomes spatial resolution, the diameter of the optical fiber 34 can be selected to be 0.1 mm. For example, if the resolution is 0.1 mm and the diameter of the optical fiber 34 is 0.25 mm, the resolution is 0.25 mm. Furthermore, higher resolution can be achieved in proportion to the magnification of the optical system.
  • a three-dimensional mapping stereoscopic image of the blood flow of the capillary is obtained.
  • the blood flow velocity distribution in an arbitrary blood vessel for example, the carotid artery
  • spatiotemporal information of blood flow velocity distribution can be obtained.
  • mapping information and blood flow velocity information of capillaries and new blood vessels that are three-dimensionally arranged can be obtained, and it is also possible to obtain information on the progress of cancer and the effectiveness of anticancer agents.
  • the effective inner diameter of a thick blood vessel can be estimated.
  • the same components as those of the above-described embodiment are denoted by the same reference numerals and description thereof is omitted.
  • the laser light is branched by the beam splitter 54 to form the sheet light Ls.
  • the sheet light Ls is transmitted through the thin living body B, and an image of scattered light transmitted through the linear irradiation site Lx is transmitted to the optical fiber array 32 via the lens system 30. It is incident.
  • the degree of freedom of arrangement of the optical system is high, and an image can be easily obtained.
  • a third embodiment of the laser Doppler blood flow measuring method and apparatus of the present invention will be described with reference to FIG.
  • the same components as those of the above-described embodiment are denoted by the same reference numerals and description thereof is omitted.
  • a cylindrical lens 22 that forms sheet light is disposed on the semiconductor laser 12 side with respect to the beam splitter 54.
  • a cylindrical lens 56 is disposed in the imaging optical system, and an optical fiber ray 32 is disposed between the sheet lights Ls.
  • the scattered light reflected by the linear irradiation site Lx is guided to the optical fiber ray 32, the imaging optical system is located between the sheet lights Ls, the arrangement of the optical system becomes compact, and the apparatus is It can be downsized.
  • the laser Doppler blood flow measuring apparatus of this embodiment is an optical system similar to that shown in FIG. 11, and the cylindrical lens 22 that forms sheet light is placed on the side of the living body B to be measured with respect to the beam splitter 54. It is arranged. Further, a concave lens 58 that expands the laser light is disposed in front of the cylindrical lens 22, and the lens system 30 of the optical system for increasing the number of lenses is controlled by increasing the number of convex lenses.
  • the laser Doppler blood flow measuring method and apparatus is not limited to the above embodiment, and the light receiving element and the photoelectric conversion element may be a solid-state imaging element such as a CCD, and the solid-state imaging element may be a line sensor or An image sensor can be used.
  • the optical fiber can be appropriately selected from glass fiber, plastic fiber, and the like.
  • the sheet light Ls by the laser beam is scanned, the linear irradiation site Lx of the sheet light Ls is transmitted through the thin living body B, the image is detected, and the in vivo tissue
  • the blood flow rate was calculated and mapped.
  • the Z-axis direction which is the depth direction of the living body B
  • the mapping results are shown in FIGS. 13 (b) and 14 (b).
  • the blood vessel position in the mapping region and the blood flow velocity are shown in each figure (b) for the blood flow of the blood vessels shown in each photograph (a) in FIGS.
  • the sheet light Ls allows the blood flow velocity to pass through a thin living tissue, and the arrangement of blood vessels can also be confirmed.
  • the arrangement of blood vessels can be visually recognized in three dimensions, and the blood vessel growth of cancer cell tissue can be reliably grasped. .
  • the position of the blood vessel in each mapping region and the blood flow velocity are shown in FIGS. This shows that the three-dimensional blood vessel arrangement can be surely grasped also by the reflection type.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Hematology (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Immunology (AREA)
  • Vascular Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

半導体レーザ12からのレーザ光を分岐してシリンドリカルレンズ22によりシート光Lsを形成し、所定位置で互いに交差させる。シート光Lsが交差した線状照射部位Lxでの散乱光を直線状に結像させるレンズ系30を有する。レンズ系30の結像位置に、複数の光ファイバ34の光ファイバアレイ32を備える。血流によるドップラー周波数シフトした散乱光を、光ファイバ34毎に電気信号に変換するアバランシェフォトダイオード42を備える。レーザ光がドップラー周波数シフトした散乱光より、線状照射部位Lxでの生体組織内の血流速を、各光ファイバ34毎に演算するとともに走査して、生体内の所定領域での血管の血流速を算出する、測定精度が高く、広い範囲を効率的に多次元計測が可能なレーザドップラー血流測定方法及び装置。 

Description

レーザドップラー血流測定方法及び装置
 この発明は、生体の血流状態をレーザドップラー方式によって測定する血流測定方法及びその測定に使用する血流測定装置に関する。
 ヒトの皮膚の血流に関するデータは、皮膚組織の代謝活性の状態、新陳代謝の状態、老化度、皮膚ガン、その他種々の皮膚の状態に関係するものであり、皮膚の血流の測定は、皮膚の状態を判断する上で重要である。
 また、日本人の死亡原因第1位は悪性新生物、第2位は心疾患、第3位は脳血管疾患であり、血管系病変による死亡は悪性新生物とほぼ同じ割合である。しかも悪性新生物の発生や進行を認識するには、周囲の新生血管の3次元マッピング情報や、新生血管内の血流情報が重要である。さらに、血管病変には好発部位があり、血管内の血流速分布などの時空間情報が不可欠である。
 近年の血流測定方法には、レーザドップラー流速測定法、粒子画像流速測定法(PIV),マイクロPIV法、超音波血流計法、MRIやCT法、及び数値シミュレーションとの融合による方法などがある。生体に関する実験は体内(In vivo)での実験、あるいは体外(In vitro)での実験がある。一般に二次元流れの情報が得られるPIV計測はインビトロ計測である。また、現在のMRIやCTの性能では空間および時間分解能が十分でないため、インビボ計測により詳細な血行動態の情報を得るのは困難である。
 そこで、レーザドップラー方式による測定方法は、血流測定方法として皮膚に侵襲を与えず、血流速及び/または血流量の瞬間値を連続的に測定できることから広く利用されている。レーザドップラー方式による血流測定方法は、レーザ光の単色性、可干渉性を利用したもので、レーザ光が皮膚組織内の血管に照射されると、血管内で動いている血球で散乱した光がドップラー効果により周波数シフトを生じるので、この周波数シフト量を基にして血流速を求めるものである。
 レーザドップラー方式の血流測定装置としては、光源から皮膚表面にレーザ光を入射させる送光用光ファイバと、皮膚からの散乱光を受光する受光用光ファイバとが一対になったセンサ素子を設け、これらを皮膚に対面するプローブ支持部材に埋め込んだものがある。このレーザドップラー法のレーザ照射方式としては、1本のビームを被検体内に入射する1光束方式と、1本のビームを2本の平行光に分け、レンズ等で交差させて、交差点(測定点)を被検体内に入射させる2光束式とがある。
 さらに、2光束式の測定方法には差動型法と参照光法がある。参照光法は、1本のレーザビームを強弱の2本のビームに分割し、ドップラー効果により周波数シフトした光と、静止した組織により散乱した光との干渉を利用している。この両者の散乱光には数百Hz~数十Hz程度の周波数の差が生じているので、両者の光の干渉をフォトダイオード等の受光素子により電気的なビート信号として検出し、検出した周波数と強度が、血球の速度及び個数に対応した値の信号となる。さらに、これらの信号の値を積分して、その血管の血流の流速と流量に関する信号に変換し、血管の血流量を求めるものである。
 また、差動型法では、1本のレーザビームを2本のビームに分割し、この2本のビームを集光して交差させ、その交差位置でレーザ光照射方向により散乱光の干渉が生じる。この干渉縞の間隔は、観測粒子である血球による散乱光のドップラー周波数シフト量により差が生じるので、これを観測して血流速を求めるものである。差動型法は、広い受光開口で散乱光を集光することができ、比較的強い信号強度が得られる利点がある。
 また、特許文献1に開示されているように、被検体の表面血流を検出して出力信号を発する複数のセンサ素子を備え被検体に当接するセンサ部と、センサ素子からの出力信号を所定の測定信号に変換する信号変換部と、測定信号を測定値として表示する表示部とを備えた表面血流測定装置も提案されている。この表面血流測定装置は、センサ部を被検体に当てると、被検体の表面血流に応じて各センサ素子から出力信号が発せられ、出力信号変換部にて所定の測定信号に変換される。センサ部には、センサ素子が複数設けられているので、一つの測定部位で、同時に又は選択的に複数の箇所を測定できる。さらに、各測定値の変動を平均化することにより、キャリブレーションが可能となり、定量的に測定することができるものである。
 また、特許文献2に開示されているように、レーザドップラー方式の血流測定方法において、被検体にレーザ光を入射させる少なくとも一つの送光部に対し、該送光部からの距離が異なる複数の位置を受光部とし、各受光部で被検体からの散乱光を受光してそれぞれの受光部における血流量を測定し、それにより測定深度の異なる血流量を同時に記録する方法も提案されている。この場合、受光部は、送光部からの距離ごとにそれぞれ複数設けることが好ましい。また、被検体に入射させるレーザ光のパワーを変化させることが好ましい。これにより、皮膚組織における血流量の変化を解析でき、皮膚表層の浅い部位から深い部位に至る広い深度範囲の血流量をそれぞれの深度ごとに同時に記録することができる。
 さらに、特許文献3に開示されているように、レーザ光をシリンドリカルレンズによりシート光に変換し、レンズを介してレンズの焦点位置に配置されたミラーに入射され、ミラーの偏向方向に反射させて照射し、血流量を2次元測定する血流分布測定装置もある。
 その他、非特許文献1~5に開示されているように、生体ではない管路等において、流体の流速を測定する装置として、レーザ光のシート光を用いた多点同時測定が可能なレーザドップラー流速計が提案されていた。
特開平5-15501号公報 特開平8-182658号公報 特開平7-100119号公報 計測自動制御学会論文集Voi.39,No.3(2003)218-224 日本機械学会論文集(B編),69巻677号(2003-1)25-30 Experiments in Fluids24(1998)70-76 Experiments in Fluids 36(2004)274-281 SICE 2002 Aug. 5-7, 2002, Osaka, 2199-2204
 従来の1光束方式の血流計は、入射レーザ光に1本のビーム(直径1mm程度)を用い、散乱光を1組のセンサで検出している。また、一部の血流計は、散乱光を複数の光ファイバ等で検出するものや、入射レーザ光用と光信号検出用の光ファイバをそれぞれ複数備えたものもあるが、これらは1点計測の組み合わせによる多点計測と言える。従って、広い領域の血流量を効率的に求めることができず、しかも分解能が低く、限られた部位の血流量しか得られないものである。
 一方、2光束式の場合は2本のレーザビームの交差点が測定点になるので、1光束式の場合より、測定体積が小さく取れ、空間分解能(数百μm程度)が高くなる。また、血流画像化装置は、レーザビームをスキャンすることにより、血流分布を測定する。2光束式の血流計による血流画像化装置は、2光束で1点の計測を行い、そのレーザ光をスキャンすることにより血流分布を得ている。しかし、2光束方式での装置は、高時空間分解能で点計測を行うことができるが、多点を同時に計測することは困難であり、血管の位置ずれが生じると、計測される平均血流速度の誤差は大きくなるものであり、この場合も広い範囲の血流量を測定できるものではない。
 また、2光束方式の場合、光ファイバを分岐させて複数のセンサ素子を設けると、個々のセンサ素子におけるレーザパワーが低くなるので、皮膚表面の比較的浅い部分の血流量しか測定できない。さらに、多点を同時に計測することは困難であり、血管の位置ずれが生じると、計測される平均血流速度の誤差は大きくなる。また、レーザ光の出力や、センサ素子の皮膚接触面における送光用光ファイバと受光用光ファイバとの間隔、即ち送光部と受光部の間隔を変えることにより、測定深度を個別的に調整することはできても、種々の測定深度の血流量を同時に測定することはしていない。そのため、深さ方向の血流量分布を知ることができず、皮膚組織の状態分析を正確に行うことができないという問題がある。
 また、特許文献1,2に開示された血流測定装置も光ファイバを用いたもので、1点計測を基にした測定法であり、上記と同様の問題がある。特許文献3に開示された血流測定装置も、立体的な血管の配置を検知することが出来るものではない。
 その他、非特許文献1~5は、予め管路の位置や流体の流れの方向が分かっているものであり、目に見えない箇所の毛細血管の配置を検知することが出来るものではない。
 この発明は、上記従来の背景技術に鑑みて成されたもので、簡単な光学系及び装置により、測定精度が高く、広い範囲を効率的に多次元計測が可能なレーザドップラー血流測定方法及び装置を提供することを目的とする。
 この発明は、レーザ光源からのレーザ光を分岐して、分岐された各レーザ光を面状の薄いシート光にして生体内の所定位置で互いに交差させ、前記シート光が交差した線状照射部位での散乱光を線状に結像させて、その結像位置で線状に配置された複数の受光素子に入射した前記各散乱光のうちの、前記生体内の血管の血流によりドップラー周波数シフトした散乱光を、前記光電変換素子毎に電気信号に変換し、前記複数の受光素子のうちの前記レーザ光がドップラー周波数シフトして入射した前記受光素子の位置により、前記生体内の血管の位置情報を検知するとともに、前記血流による前記レーザ光のドップラー周波数シフトにより、前記血管内の血流速を前記光電変換素子毎に演算するレーザドップラー血流測定方法である。
 さらに、前記シート状のレーザ光の交差位置を、前記線状の散乱光の長手方向と交差する方向に走査して、走査された前記線状の散乱光から、前記光電変換素子毎に前記血管内の血流を求め、前記生体内の所定領域での血管の太さと配置を算出するものである。
 さらに、前記シート状のレーザ光の交差位置を前記線状の散乱光の長手方向と交差する方向に走査して、走査された前記線状の散乱光から、前記光電変換素子毎に前記血管内の血流速を求め、求めた個々の前記線状照射部位での前記血流速を基にして、前記生体内の所定領域での各血管内の血流速または血流量を算出するものである。
 前記走査は、前記線状照射部位の長手方向と直交する2方向に行い、前記線状照射部位を順次前記長手方向と平行な方向にずらして、前記生体内の所定領域での3次元的血管配置を求め、表示可能にするものである。また、前記走査により、前記生体組織の3次元的血流状態を演算し、表示可能にするものである。
 またこの発明は、レーザ光源と、前記レーザ光源からのレーザ光を分岐するビーム分割装置と、前記ビーム分割装置により分岐されたレーザ光を面状の薄いシート光にする円筒光学系と、分岐した前記各シート光を生体内の所定位置で互いに交差させる集光光学系と、前記シート光が交差した線状照射部位での散乱光を直線状に結像させる結像光学系と、前記結像光学系の結像位置に配置された複数の受光素子と、前記受光素子に入射した前記散乱光を前記光電変換素子毎に電気信号に変換する光電変換素子と、前記生体組織内の血流により前記レーザ光がドップラー周波数シフトした前記線状照射部位での血流速を、前記電気信号を基にして前記光電変換素子毎に演算する血流演算手段と、前記シート状のレーザ光の交差位置を前記線状照射部位の長手方向と交差する方向に走査するレーザ光走査手段とを備え、前記レーザ光走査手段により走査された領域の前記生体内の血管配置を、前記血流演算手段による演算を基に算出するレーザドップラー血流測定装置である。
 さらに、前記レーザ光走査手段により走査された所定の各位置での散乱光について、前記血流演算手段により求めた前記光電変換素子毎の各血流速と前記血管配置から、前記生体内の各血管の血流量を算出する血流量演算手段を設けたものである。
 また、前記結像光学系は、測定対象の生体組織に対して前記集光光学系と同じ側に位置し、分岐した前記各シート光の間に前記集光光学系が位置したものである。または、前記結像光学系は、測定対象の生体組織を挟んで前記集光光学系とは反対側に位置したものでも良い。
 また、前記受光素子は、平面上に受光素子が一体的に配置された固体撮像素子から成るものである。
 この発明のレーザドップラー血流測定方法及び装置によれば、ヒト体内での血流速の高次空間分解能での直接測定が可能となる。さらに、比較的太い血管では血流速分布とその時間変動の計測及び毛細血管の血流速の測定が可能である。そして、レーザ光を走査することにより、所定の測定範囲での、太い血管での有効内径、血流速分布とその時間変動、及び毛細血管の平均血流速の測定が可能であり、それらの測定結果を所定の測定範囲で3次元マッピングすることが可能となる。また、接触による血流への影響がなく、炎症部位、耳の毛細血管、皮膚に近い大動脈、静脈の血流と血管内の血流分布の測定に適しており、再現性の高い測定結果が得られる。
 また、予防医療にも効果があり、皮膚ガンなどの新生血管の血流速分布計測とその時空間情報、脳血管疾患と密接に関係する頸動脈血流速分布の時空間情報の計測、その他、基礎的な応用としてアレルギー性皮膚炎症の発生機序の解明にも極めて有効である。
この発明の第一実施形態のレーザドップラー血流測定装置の送光部の光学系を示す概略図である。 この発明の第一実施形態のレーザドップラー血流測定装置の受光部の結像光学系を示す概略図である。 この発明の第一実施形態のレーザドップラー血流測定装置のファイバアレイと光ファイバの配置等を示す概略図である。 この発明の第一実施形態のレーザドップラー血流測定装置の光電変換素子からコンピュータまでの機能を示すブロック図である。 この発明の第一実施形態のレーザドップラー血流測定方法を示す概念図である。 この発明の第一実施形態のレーザドップラー血流測定方法の光学系の配置等を示す概略図である。 この発明の第一実施形態のレーザドップラー血流測定方法によるレーザ光の干渉を示す写真である。 この発明の第一実施形態のレーザドップラー血流測定方法の処理を示すフローチャートである。 この発明の第一実施形態のレーザドップラー血流測定方法によるレーザ光の交差部の走査を示す概念図である。 この発明の第二実施形態のレーザドップラー血流測定装置の光学系の配置を示す概略図である。 この発明の第三実施形態のレーザドップラー血流測定装置の光学系の配置を示す概略図である。 この発明の第四実施形態のレーザドップラー血流測定装置の光学系の配置を示す概略図である。 この発明の実施例のレーザドップラー血流測定方法による血流を測定したマウスの耳を表側から見た部位の写真(a)と、レーザ光を所定深さで走査して、その測定結果のマッピングを行った図(b)である。 この発明の実施例のレーザドップラー血流測定方法による血流を測定したマウスの耳を裏側から見た部位の写真(a)と、レーザ光を所定深さで走査して、その測定結果のマッピングを行った図(b)である。 この発明の実施例のレーザドップラー血流測定方法による血流を測定したマウスの耳の血管配置を示すマッピング図である。 この発明の実施例のレーザドップラー血流測定方法による血流を測定したマウスの耳の3次元マッピング図である。 この発明の実施例のレーザドップラー血流測定方法による反射型の装置で血流を測定した写真(a)マッピングを行った図(b)(c)である。
符号の説明
10 レーザドップラー血流測定装置
12 半導体レーザ
14,15 音響光学素子
16 反射プリズム
18,19 反射鏡
20,21 ロッドレンズ
22,24 シリンドリカルレンズ
30 レンズ系
32 ファイバアレイ
34 光ファイバ
42 アバランシェフォトダイオード
Ls シート光
Lx 線状照射部位
 以下、この発明の実施形態について図面に基づいて説明する。図1~図9は、この発明の第一実施形態のレーザドップラー血流測定装置を示すもので、この実施形態のレーザドップラー血流測定装置10は、図1に示すように、レーザ光源であるレーザダイオード等の半導体レーザ12と、この半導体レーザ12からの連続発振レーザ光を周波数の若干異なる2本のビームに分岐させる音響光学素子14,15を備えている。レーザ光の波長は、波長750~1500nm程度が好ましく、例えば近赤外の波長785nmのレーザを用いる。これは、近赤外光は生体組織内への透過率が高いからである。
 さらに、音響光学素子15の光出射側には、分岐されたレーザ光を2方向に反射する反射プリズム16と、反射プリズム16により反射されたレーザ光を直角に反射する反射鏡18,19を備えている。さらに、反射鏡18,19により反射した2本のレーザ光が入射する互いに平行な位置に設けられたロッドレンズ20,21が配置されている。各ロッドレンズ20,21は、入射したレーザ光を広げて第一のシリンドリカルレンズ22の平面部22aに入射させる。
 第一のシリンドリカルレンズ22は、平面部22aがロッドレンズ20,21側に面して位置し、入射した2本のレーザ光は、第一のシリンドリカルレンズ22の円筒面状の曲面部22bから出射し、第一のシリンドリカルレンズ22の円筒の中心軸と平行であって平面部22aと直交するシート状の2本のシート光Lsに形成される。このシート光Lsは、さらに、集光光学系を形成した第二のシリンドリカルレンズ24の平面部24aに入射する。
 第二のシリンドリカルレンズ24の平面部24aは第一のシリンドリカルレンズ22の平面部22aと同じ側を向いて平行に位置し、第二のシリンドリカルレンズ24の他方の面である円筒面状の曲面部24bは、第一のシリンドリカルレンズ22に対して、円筒の中心軸が互いにねじれの位置関係で直交するように配置されている。なお、シリンドリカルレンズ22,24は、平面部22a,24aと曲面部22b,24bのいずれの面が入射面となっても、出射光は同様の光路へ進む。
 この実施形態に用いた光学系部品の表面には、用いたレーザ光の波長の減衰を抑えるために、この波長に対する反射防止膜がコーティングされている。
 第二のシリンドリカルレンズ24に入射した2本のシート光Lsは、曲面部24bで屈折して互いに交差するように偏向する。シート光Lsの交差部は、1本の線状に光強度の高い部分として形成され、生体組織の線状照射部位Lxとなる。
 線状照射部位Lxでの散乱光は、図2に示すような結像光学系を形成したレンズ系30により収束される。レンズ系30は、線状照射部位Lxの散乱光を、受光素子である光ファイバアレイ32の端面部32aに線状照射部位Lxに沿って線状に結像させるものである。
 光ファイバレイ32は、図3に示すように、ライトガイド束36に束ねられた光ファイバ34を経て、他方の端部に導かれ、光ファイバ34毎に光プラグ38に繋がっている。光プラグ38は、各々光レセプタクル40に接続可能に設けられている。各光レセプタクル40には、図4に示すように、光電変換素子であるアバランシェフォトダイオード42が光ファイバ34の数だけ設けられ、光ファイバ34により導かれた光を各々電気信号に変換する。各アバランシェフォトダイオード42の出力は、光ファイバ34の本数だけ設けられた各アンプ44を介して各々周波数フィルタ46に入力し、所定の周波数の信号がA/D変換されて、光ファイバ34の数だけ設けられたデータレコーダ48に各々記録される。各データレコーダ48に記録されたデジタル信号は、USBコネクタを介してコンピュータ50に入力され、所定の演算解析処理が行われる。
 次に、この実施形態のレーザドップラー血流測定方法について以下に説明する。このレーザドップラー血流測定方法は、半導体レーザ12から出射されたレーザ光を、音響光学素子14,15により変調して分岐させ、図5、図6に示すように、シリンドリカルレンズ22により面状の薄いシート光Lsにするとともに、シリンドリカルレンズ24により互いに交差させ、血管52のある生体B内の所定位置に線状照射部位Lxを形成する。この状態で、血管52内の血流速及びその血流速から必要に応じて血流量を、光ファイバ34を単位として結像光学系であるレンズ系30の倍率により決まる空間分解能で求めることができる。ここでは、250μmの径の光ファイバ34を用いて光学系の倍率を1とすると、空間分解能は250μmとなる。
 シート光Lsの交差部は、1本の線状に光強度の高い部分として形成され、生体組織の線状照射部位Lxとなる。この線状部位Lxは、レーザ光の入射波面の位相差によって干渉縞ができており、図7に示すような干渉縞となって現れる。この線状の領域を赤血球がある速度で通過すると、散乱光のドップラー周波数が変化する。この線状照射部位Lxでは、シート光Lsが照射された血流中の血球により散乱したレーザ光が、血流の流速によるドップラー効果により、散乱光のドップラー周波数が変化する。なお、静止した組織では散乱光のドップラー周波数の変化はない。線状照射部位Lxでの各散乱光は、それぞれレンズ系30により光ファイバレイ32上のそれぞれに対応する点に結像される。
 この周波数変化した散乱光を、光ファイバ34を介してアバランシェフォトダイオード42に入力し、電気的なビート信号として検出し、検出した周波数変化と強度が、血球の速度及び個数に対応した値である。これを各受光部である光ファイバ34毎に求めて、対応する部位の血流の流速を算出し、その算出結果を線状照射部位Lx上で積分し、その線上の血流量を求める。ここでの処理は、図8に示すように、アバランシェフォトダイオード42により得られた入力信号は、例えば図8のグラフ(a)のようなデータであり、このデータがA/D変換されてコンピュータ入力され、そのデジタル信号を高速フーリエ変換し、さらにノイズ除去処理等を施して、ドップラー周波数を算出する。ドップラー周波数は、図8(b)に示すように、検出した信号の周波数スペクトルにおいて、ドップラー周波数でのピークが表れることにより求まる。ドップラー周波数が求まると、その周波数から血球の速度が算出され、血球の速度は血流の流速であることから血流速が算出される。
 また、この実施形態では、音響光学素子14により変調をかけているので、血液の流れの方向も識別することができる。ただし、流れの方向を識別する必要がない場合は、音響光学素子14による周波数変調をしなくても計測が可能である。
 さらに、この線状照射部位Lxをその長手方向に対して直交する方向に走査して、所定のサンプリング周期で各光ファイバ34毎に流量を求め、所定面内での血流量のマッピングが形成される。またこれを、線状照射部位Lxの長手方向に対して直交する他の方向に走査して、3次元的な血流速または血流量のマッピングが可能となる。
 ここで、線状照射部位Lxの走査は、測定対象の生体が小さい場合は、生体自体を互いに直交するXYZ軸方向に移動可能な台に載せて、この3軸方向に移動させることにより走査可能である。また、ヒト等の大きな生体の場合は、シリンドリカルレンズ22,24等の光学系を駆動して走査することができる。この時のシート光Lsの交差位置とその走査方向は、直線の交点を求める演算により容易に求めることができる。
 この走査は、X軸方向の流速を求める場合、図9(a)に示す通り、生体Bの組織内に線状照射部位Lxを形成し、その線状照射部位Lxを、その長手方向と直交するX軸方向に相対的に所定距離移動させる。この後、さらにY軸方向に線状照射部位Lxを、その照射長さ分だけ平行移動させ、同様に走査を繰り返す。次に、図9(b)に示す通り、直交するZ軸方向に微小距離移動してさらに上記と同様に走査し、所定範囲でこの移動と走査を繰り返す。これにより、生体Bの所定範囲の3次元的血管配置とその血管の血流速または血流量を求めることが出来る。また、Y軸方向の流速を求める場合は、線状照射部位LxをX軸方向と平行にして、図9(c)に示す通り、Y軸方向に走査するとともにX軸方向にシフトして走査し、さらにZ軸方向に移動し、走査を繰り返す。
 この実施形態のレーザドップラー血流測定方法及び装置によれば、血管52の内径や毛細血管の血流状態を計測し、三次元マッピングすることができる。この場合の空間分解能は、数百ミクロン程度にすることができる。即ち、シート光Lsの交差線の線状照射部位Lxでの散乱光を受光する光ファイバレイ32の光ファイバ34の内径が、空間分解能になるので、光ファイバ34の径を0.1mmに選べば、分解能は0.1mm、光ファイバ34の径を0.25mmに選べば、分解能は0.25mmになる。さらに、光学系の倍率に比例して、より高い分解能とすることができる。
 以上の走査を繰り返すことにより、毛細血管の血流の3次元マッピングの立体像が得られる。これにより、任意の血管内(例えば頸動脈)の血流速分布及びその時間変化を測定することができる。また、血流速分布の時空間情報が得られる。さらに、立体的に配置している毛細血管、新生血管のマッピング情報と血流速度情報が得られ、ガンの進行状態や抗ガン剤の効き目に関する情報を得ることも可能となる。また、太い血管での有効内径を見積もることも出来る。
 次に、この発明のレーザドップラー血流測定方法及び装置の第二実施形態について、図10を基にして説明する。ここで、上述の実施形態と同様の構成は、同一の符号を付して説明を省略する。この実施形態のレーザドップラー血流測定方装置は、ビームスプリッタ54によりレーザ光を分岐してシート光Lsを形成したものである。
 この実施形態のレーザドップラー血流測定装置は、シート光Lsを薄い生体Bに透過させて、その線状照射部位Lxで透過した散乱光の像を、レンズ系30を介して光ファイバレイ32に入射させたものである。この場合、光学系の配置の自由度が高く、像を得やすい。
 次に、この発明のレーザドップラー血流測定方法及び装置の第三実施形態について、図11を基にして説明する。ここで、上述の実施形態と同様の構成は、同一の符号を付して説明を省略する。この実施形態のレーザドップラー血流測定装置は、シート光を形成するシリンドリカルレンズ22を、ビームスプリッタ54に対して半導体レーザ12側に配置したものである。さらに、結像光学系にシリンドリカルレンズ56を配置するとともに、光ファイバレイ32をシート光Lsの間に配置したものである。
 この場合、線状照射部位Lxで反射した散乱光を、光ファイバレイ32に導くもので、結像光学系がシート光Lsの間に位置して、光学系の配置がコンパクトになり、装置を小型化することができる。
 次に、この発明のレーザドップラー血流測定方法及び装置の第四実施形態について、図12を基にして説明する。ここで、上述の実施形態と同様の構成は、同一の符号を付して説明を省略する。この実施形態のレーザドップラー血流測定方装置は、図11に示すものと同様の光学系であって、シート光を形成するシリンドリカルレンズ22を、ビームスプリッタ54に対して測定対象の生体B側に配置したものである。さらに、レーザ光を拡張する凹レンズ58をシリンドリカルレンズ22の前に配置し、結増光学系のレンズ系30を凸レンズの枚数を増やして収差を抑えたものである。
 この実施形態の場合、上述の実施形態と同様の効果に加えて、シート光Lsの幅が数mmの場合に有効であり、収差が少ないものである。
 尚、この発明のレーザドップラー血流測定方法及び装置は、上記実施形態に限定されるものではなく、受光素子及び光電変換素子は、CCD等の固体撮像素子でも良く、固体撮像素子はラインセンサ又はイメージセンサを用いることができる。また、光ファイバは、ガラスファイバやプラスチックファイバ等、適宜選択可能である。
 次に、この発明のレーザドップラー血流測定方法と装置を用いた血流測定の実施例を示す。この実施例では、図9に示すように、レーザ光によるシート光Lsを走査し、シート光Lsの線状照射部位Lxを、薄い生体Bを透過させて、その像を検出し、生体内組織の血流速を算出してマッピングした。生体Bの深さ方向であるZ軸方向について、表面側からZ=0.25mmとZ=0.50mmの異なる深さでXY軸方向の走査を行った。そのマッピング結果を図13(b),図14(b)に示す。図13、図14の各写真(a)に表れた血管の血流について、各々マッピング領域の血管の位置と、血流速を各図(b)に示す。
 この実施例によれば、シート光Lsにより、薄い生体組織を透過して血流の速度を測定することができ、血管の配置も確認することができた。
 図15、図16は、上記実施例により、さらに深さ方向にシート光Lsによる走査を行い、Z=0.75mm、Z=1.00mmの深さでの血管の配置を計測し、それを立体的に表示したものである。
 この発明のレーザドップラー血流測定方法と装置によれば、図16に示すように、血管の配置を立体的に視認することが出来、ガン細胞組織の血管成長等を確実に把握することが出来る。
 また、図17に示す実施例は、反射型の血流測定装置による測定例を示すもので、図17(a)に示す写真の破線で囲んだ領域について、表面側からZ=0.25mmとZ=0.75mmの異なる深さで、XY軸方向の走査を行った。各々マッピング領域の血管の位置と、血流速を図(b)(c)に示す。反射型によっても、確実に立体的な血管配置を把握することが出来たことを示す。
 

Claims (10)

  1.  レーザ光源からのレーザ光を分岐して、分岐された各レーザ光を面状の薄いシート光にして生体内の所定位置で互いに交差させ、
     前記シート光が交差した線状照射部位での散乱光を線状に結像させて、
     その結像位置で線状に配置された複数の受光素子に入射した前記各散乱光のうちの、前記生体内の血管の血流によりドップラー周波数シフトした散乱光を、複数の光電変素子により電気信号に変換し、
     前記複数の受光素子のうちの、前記レーザ光がドップラー周波数シフトして入射した前記受光素子の位置により、前記生体内の血管の位置情報を検知するとともに、
     前記血流による前記レーザ光のドップラー周波数シフトにより、前記血管内の血流速を前記光電変換素子毎に演算することを特徴とするレーザドップラー血流測定方法。
  2.  前記シート状のレーザ光の交差位置を、前記線状の散乱光の長手方向と交差する方向に走査して、走査された前記線状の散乱光から、前記光電変換素子毎に前記血管内の血流を求め、前記生体内の所定領域での血管の太さと配置を算出することを特徴とする請求項1記載のレーザドップラー血流測定方法。
  3.  前記シート状のレーザ光の交差位置を前記線状の散乱光の長手方向と交差する方向に走査して、走査された前記線状の散乱光から、前記光電変換素子毎に前記血管内の血流速を求め、求めた個々の前記線状照射部位での前記血流速を基にして、前記生体内の所定領域での血管毎の血流速または血流量を算出することを特徴とする請求項2記載のレーザドップラー血流測定方法。
  4.  前記走査は、前記線状照射部位の長手方向と直交する2方向に行い、前記線状照射部位を順次前記長手方向と平行な方向にずらして、前記生体内の所定領域での3次元的血管配置を求め、前記血管配置を表示可能にすることを特徴とする請求項3記載のレーザドップラー血流測定方法。
  5.  前記走査により、前記生体組織内の血管の3次元的血流状態を演算し、表示可能にすることを特徴とする請求項4記載のレーザドップラー血流測定方法。
  6.  レーザ光源と、前記レーザ光源からのレーザ光を分岐するビーム分割装置と、前記ビーム分割装置により分岐されたレーザ光を面状の薄いシート光にする円筒光学系と、分岐した前記各シート光を生体内の所定位置で互いに交差させる集光光学系と、前記シート光が交差した線状照射部位での散乱光を直線状に結像させる結像光学系と、前記結像光学系の結像位置に配置された複数の受光素子と、前記受光素子に入射した前記散乱光を電気信号に変換する光電変換素子と、前記生体組織内の血流により前記レーザ光がドップラー周波数シフトした前記線状照射部位での血流速を、前記電気信号を基にして前記光電変換素子毎に演算する血流演算手段と、前記シート状のレーザ光の交差位置を前記線状照射部位の長手方向と交差する方向に走査するレーザ光走査手段とを備え、前記レーザ光走査手段により走査された領域の前記生体内の血管配置を、前記血流演算手段による演算を基に算出することを特徴とするレーザドップラー血流測定装置。
  7.  前記レーザ光走査手段により走査された所定の各位置での散乱光について、前記血流演算手段により求めた前記光電変換素子毎の各血流速と前記血管配置から、前記生体内の各血管の血流量を算出する血流量演算手段を設けたことを特徴とする請求項6記載のレーザドップラー血流測定装置。
  8.  前記結像光学系は、測定対象の生体組織に対して前記集光光学系と同じ側に位置し、分岐した前記各シート光の間に前記集光光学系が位置したことを特徴とする請求項6記載のレーザドップラー血流測定装置。
  9.  前記結像光学系は、測定対象の生体組織を挟んで前記集光光学系とは反対側に位置したことを特徴とする請求項6記載のレーザドップラー血流測定装置。
  10.  前記受光素子は、平面上に光電変換素子が一体的に配置された固体撮像素子から成ることを特徴とする請求項6記載のレーザドップラー血流測定装置。
     
PCT/JP2008/073245 2007-12-21 2008-12-19 レーザドップラー血流測定方法及び装置 Ceased WO2009081883A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US12/809,738 US9031640B2 (en) 2007-12-21 2008-12-19 Laser doppler blood flow measuring method and device
JP2009547086A JP5234470B2 (ja) 2007-12-21 2008-12-19 レーザドップラー血流測定方法及び装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2007330112 2007-12-21
JP2007-330112 2007-12-21

Publications (1)

Publication Number Publication Date
WO2009081883A1 true WO2009081883A1 (ja) 2009-07-02

Family

ID=40801184

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2008/073245 Ceased WO2009081883A1 (ja) 2007-12-21 2008-12-19 レーザドップラー血流測定方法及び装置

Country Status (3)

Country Link
US (1) US9031640B2 (ja)
JP (1) JP5234470B2 (ja)
WO (1) WO2009081883A1 (ja)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012176095A (ja) * 2011-02-25 2012-09-13 Canon Inc 画像処理装置及び画像処理システム
WO2016117106A1 (ja) * 2015-01-23 2016-07-28 オリンパス株式会社 外科処置装置
WO2016170897A1 (ja) * 2015-04-21 2016-10-27 オリンパス株式会社 医療装置、医療装置の作動方法
WO2016171274A1 (ja) * 2015-04-23 2016-10-27 オリンパス株式会社 内視鏡装置
JP2019097574A (ja) * 2017-11-28 2019-06-24 愛知時計電機株式会社 測定装置
JP2020188916A (ja) * 2019-05-21 2020-11-26 愛知時計電機株式会社 測定装置
US10856751B2 (en) 2015-03-25 2020-12-08 Olympus Corporation Blood-vessel recognition blood-flow measurement method
WO2020255840A1 (ja) * 2019-06-20 2020-12-24 ソニー株式会社 情報処理装置、情報処理方法、及びプログラム
JP2021025823A (ja) * 2019-08-01 2021-02-22 株式会社ジェイ・エム・エス 流路内の流動解析装置
JPWO2021060364A1 (ja) * 2019-09-27 2021-04-01
JP2021108791A (ja) * 2020-01-07 2021-08-02 国立大学法人富山大学 悪性黒色腫診断装置
CN113423335A (zh) * 2019-02-11 2021-09-21 株式会社高迎科技 血流测量装置及血流测量方法

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102010055772B4 (de) * 2010-12-23 2014-06-05 Carl Zeiss Meditec Ag Anordnung und Verfahren zur quantitativen Bestimmung des Blutflusses innerhalb von Blutgefäßen
CN103932662A (zh) * 2011-12-30 2014-07-23 广州宝胆医疗器械科技有限公司 多普勒激光气管镜系统
RU2511606C2 (ru) * 2012-08-22 2014-04-10 Федеральное государственное унитарное предприятие "Всероссийский научно-исследовательский институт автоматики им. Н.Л. Духова" (ФГУП "ВНИИА") Устройство доплеровского измерителя скорости на основе интерферометра фабри-перо с волоконным вводом излучения
US9121861B2 (en) * 2013-03-15 2015-09-01 National Security Technologies, Llc Photonic Doppler velocimetry lens array probe incorporating stereo imaging
US9844320B2 (en) * 2014-01-29 2017-12-19 University Of Rochester System and method for observing an object in a blood vessel
US10799129B2 (en) * 2016-01-07 2020-10-13 Panasonic Intellectual Property Management Co., Ltd. Biological information measuring device including light source, light detector, and control circuit
JP6376412B2 (ja) * 2016-02-16 2018-08-22 ウシオ電機株式会社 レーザシート光源装置
US12144655B2 (en) 2021-08-23 2024-11-19 Samsung Electronics Co., Ltd. Method and device for liveness detection
CN114947795A (zh) * 2022-04-19 2022-08-30 中国科学院苏州纳米技术与纳米仿生研究所 一种血流流速的检测方法和检测装置
WO2024091180A1 (en) * 2022-10-26 2024-05-02 National University Of Singapore Light sheet imaging apparatus and method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07100119A (ja) * 1993-10-04 1995-04-18 Kowa Co 血流分布測定装置
JPH1085195A (ja) * 1996-09-19 1998-04-07 Advance Co Ltd レーザー血流計

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3941477A (en) * 1974-10-17 1976-03-02 Deutsche Forschungs-Und Versuchsanstalt Fur Luft-Und Raumfahrt E.V. Measuring device for the measurement of fluid flow rates
US4154529A (en) * 1977-03-21 1979-05-15 Andrew Corporation System for detecting reflected laser beams
JP2641220B2 (ja) * 1987-10-28 1997-08-13 興和株式会社 速度分布測定装置
JP3091269B2 (ja) 1990-08-31 2000-09-25 花王株式会社 表面血流測定装置
US6485413B1 (en) * 1991-04-29 2002-11-26 The General Hospital Corporation Methods and apparatus for forward-directed optical scanning instruments
JP3313841B2 (ja) * 1993-09-24 2002-08-12 興和株式会社 血流測定装置
JP3109066B2 (ja) 1994-12-28 2000-11-13 花王株式会社 血流測定方法及び血流測定装置
US6540981B2 (en) * 1997-12-04 2003-04-01 Amersham Health As Light imaging contrast agents
US5999836A (en) * 1995-06-06 1999-12-07 Nelson; Robert S. Enhanced high resolution breast imaging device and method utilizing non-ionizing radiation of narrow spectral bandwidth
US6258046B1 (en) * 1995-07-06 2001-07-10 Institute Of Critical Care Medicine Method and device for assessing perfusion failure in a patient by measurement of blood flow
JP3679547B2 (ja) * 1997-04-02 2005-08-03 キヤノン株式会社 眼底血管観察装置
US6006128A (en) * 1997-06-02 1999-12-21 Izatt; Joseph A. Doppler flow imaging using optical coherence tomography
US5865828A (en) * 1997-08-08 1999-02-02 Jeng; James C. Coaxial dual laser
US6549801B1 (en) * 1998-06-11 2003-04-15 The Regents Of The University Of California Phase-resolved optical coherence tomography and optical doppler tomography for imaging fluid flow in tissue with fast scanning speed and high velocity sensitivity
US6569104B2 (en) * 1998-07-16 2003-05-27 Canon Kabushiki Kaisha Blood vessel detecting apparatus
US6332683B1 (en) * 1999-10-15 2001-12-25 Canon Kabushiki Kaisha Fundus examination apparatus
GB0010123D0 (en) * 2000-04-27 2000-06-14 Univ Nottingham Planar light sheet anemometers
JP4708543B2 (ja) * 2000-06-14 2011-06-22 キヤノン株式会社 眼血流計
JP2002034921A (ja) * 2000-07-19 2002-02-05 Canon Inc 眼底検査装置
WO2006111909A1 (en) * 2005-04-20 2006-10-26 Cvl Cosmetics S.A. Instrument and method for high-speed perfusion imaging
JP6316298B2 (ja) * 2012-09-10 2018-04-25 オレゴン ヘルス アンド サイエンス ユニバーシティ Oct血管造影による局所循環の定量化
US8843186B2 (en) * 2012-11-21 2014-09-23 Folim G. Halaka Non-invasive reagentless glucose determination
US9025159B2 (en) * 2012-12-10 2015-05-05 The Johns Hopkins University Real-time 3D and 4D fourier domain doppler optical coherence tomography system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07100119A (ja) * 1993-10-04 1995-04-18 Kowa Co 血流分布測定装置
JPH1085195A (ja) * 1996-09-19 1998-04-07 Advance Co Ltd レーザー血流計

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
NORIYUKI FURUICHI ET AL.: "Taten Doji Sokutei LDV to CCD Area Image Sensor o Heiyo shita 2 Seibun Ryusoku Bunpu Doji Sokuteiho no Kaihatsu", TRANSACTIONS OF THE SOCIETY OF INSTRUMENT AND CONTROL ENGINEERS, vol. 39, no. 3, 31 March 2003 (2003-03-31), pages 218 - 224 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012176095A (ja) * 2011-02-25 2012-09-13 Canon Inc 画像処理装置及び画像処理システム
WO2016117106A1 (ja) * 2015-01-23 2016-07-28 オリンパス株式会社 外科処置装置
US10856751B2 (en) 2015-03-25 2020-12-08 Olympus Corporation Blood-vessel recognition blood-flow measurement method
WO2016170897A1 (ja) * 2015-04-21 2016-10-27 オリンパス株式会社 医療装置、医療装置の作動方法
JP6095863B1 (ja) * 2015-04-21 2017-03-15 オリンパス株式会社 医療装置
WO2016171274A1 (ja) * 2015-04-23 2016-10-27 オリンパス株式会社 内視鏡装置
JP2019097574A (ja) * 2017-11-28 2019-06-24 愛知時計電機株式会社 測定装置
CN113423335A (zh) * 2019-02-11 2021-09-21 株式会社高迎科技 血流测量装置及血流测量方法
US12193802B2 (en) 2019-02-11 2025-01-14 Koh Young Technology Inc. Blood flow measurement device and blood flow measurement method
JP7244351B2 (ja) 2019-05-21 2023-03-22 愛知時計電機株式会社 測定装置
JP2020188916A (ja) * 2019-05-21 2020-11-26 愛知時計電機株式会社 測定装置
JPWO2020255840A1 (ja) * 2019-06-20 2020-12-24
WO2020255840A1 (ja) * 2019-06-20 2020-12-24 ソニー株式会社 情報処理装置、情報処理方法、及びプログラム
JP7537431B2 (ja) 2019-06-20 2024-08-21 ソニーグループ株式会社 情報処理装置、情報処理方法、及びプログラム
JP2021025823A (ja) * 2019-08-01 2021-02-22 株式会社ジェイ・エム・エス 流路内の流動解析装置
JPWO2021060364A1 (ja) * 2019-09-27 2021-04-01
WO2021060364A1 (ja) * 2019-09-27 2021-04-01 富士フイルム株式会社 流量測定装置
JP7355836B2 (ja) 2019-09-27 2023-10-03 富士フイルム株式会社 流量測定装置
JP2021108791A (ja) * 2020-01-07 2021-08-02 国立大学法人富山大学 悪性黒色腫診断装置
JP7385828B2 (ja) 2020-01-07 2023-11-24 学校法人金城学院 悪性黒色腫診断装置

Also Published As

Publication number Publication date
US9031640B2 (en) 2015-05-12
US20100280398A1 (en) 2010-11-04
JP5234470B2 (ja) 2013-07-10
JPWO2009081883A1 (ja) 2011-05-06

Similar Documents

Publication Publication Date Title
JP5234470B2 (ja) レーザドップラー血流測定方法及び装置
EP0904011B1 (en) Apparatus for imaging microvascular blood flow
JP5844510B2 (ja) 生体環境中での構造および流れの撮像
US20090219544A1 (en) Systems, methods and computer-accessible medium for providing spectral-domain optical coherence phase microscopy for cell and deep tissue imaging
EP1362252A4 (en) METHOD AND APPARATUS FOR DIFFERENTIAL PHASE OPTICAL COHERENCE TOMOGRAPHY
US5778878A (en) Laser doppler microscopy methods and instruments
JP2008542758A (ja) スペクトルコード化ヘテロダイン干渉法を画像化に使用可能なシステム、方法、及び装置
WO2018015517A1 (en) An integrated lens free imaging device
CN105996999B (zh) 基于oct测量样品深度分辨衰减系数的方法和系统
WO2020167870A1 (en) Transparent ultrasound transducers for photoacoustic imaging
CN115517628B (zh) 一种用于检测自由移动小鼠循环肿瘤细胞的光声系统
Bauer-Marschallinger et al. Fiber-optic annular detector array for large depth of field photoacoustic macroscopy
CN105865613A (zh) 海洋立体监测中的水下光学检测与成像传感方法及系统
JP6265364B2 (ja) 細胞識別装置及び細胞識別方法、並びに、細胞識別方法のプログラム及びそのプログラムを記録した記録媒体
CN106510644B (zh) 基于光纤束的医用光学相干层析成像二维前向扫描探头
Shabairou et al. All-optical, an ultra-thin endoscopic photoacoustic sensor using multi-mode fiber
CN203828901U (zh) 用于频域oct系统的光谱仪
US12270646B2 (en) Broadband interferometric confocal microscope
Akiguchi et al. Measurement of blood flow velocity in a model of stenosis in vitro and in mesenteric vessels in vivo using non-invasive micro multipoint laser Doppler velocimetry
CN110575142A (zh) 一种单光谱仪多光束光学相干层析成像仪
Wurzinger et al. Combined photoacoustic, pulse-echo laser ultrasound, and speed-of-sound imaging using integrating optical detection
JP7029788B2 (ja) 流量測定装置、流量測定方法および流量測定プログラム
CN212300424U (zh) 结合oct成像以及压力传感的光纤探头
Fedosov et al. Bioflow measuring: laser Doppler and speckle techniques
CN211213098U (zh) 一种单光谱仪多光束光学相干层析成像仪

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

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2009547086

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 12809738

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 08863983

Country of ref document: EP

Kind code of ref document: A1