WO2024262621A1 - ラマン散乱光を用いたシステムおよび方法 - Google Patents
ラマン散乱光を用いたシステムおよび方法 Download PDFInfo
- Publication number
- WO2024262621A1 WO2024262621A1 PCT/JP2024/022648 JP2024022648W WO2024262621A1 WO 2024262621 A1 WO2024262621 A1 WO 2024262621A1 JP 2024022648 W JP2024022648 W JP 2024022648W WO 2024262621 A1 WO2024262621 A1 WO 2024262621A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- interference
- scattered light
- light
- interference pattern
- raman scattered
- 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.)
- Pending
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration or pH-value ; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid or cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration or pH-value ; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid or cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/65—Raman scattering
Definitions
- the present invention relates to a system and method using Raman scattered light.
- JP Patent Publication 2013-517516 discloses an apparatus and method for reducing laser speckle by using stimulated Raman scattering in optical fibers.
- the fiber core diameter and length are selected to obtain the desired output light, and tunable despeckle is formed by combining two optical fibers in parallel and adjusting the amount of light in each path with the aid of a rotatable waveplate and polarizer as a beam splitter.
- JP Patent Publication No. 2010-532699 discloses a laser speckle imaging system and method.
- This document provides an apparatus and method for measuring perfusion in tissue, the method including the steps of recording images of the tissue under laser light, calculating a plurality of contrast images from the plurality of images of the tissue, determining a power spectrum of the scattered light from the plurality of contrast images, and determining the perfusion from the power spectrum.
- Body fluids from a living organism such as blood, contain blood cells and plasma as their main components, with blood cells including red blood cells, white blood cells, and platelets. Furthermore, in order to check the health of a living organism, it is necessary to measure glucose, hemoglobin A1c, creatinine, albumin, and other components in the blood in addition to the main components, and plasma glucose level is required as the measurement standard for blood glucose level. Therefore, for example, even if information on glucose concentration is obtained from a living organism, the value or evaluation of the information as information may differ depending on whether it is possible to determine the origin of the information.
- Raman spectroscopy is a method of evaluating materials using Raman scattered light.
- the light interacts with the material to emit Raman scattered light with a different wavelength from the incident light, and a spectrum is obtained that shows the intensity as a function of frequency or wavelength.
- Chemical bond information can be obtained from the peak positions in the spectrum, and molecular structure information can be obtained from the waveform of the entire spectrum.
- it is not easy to simultaneously obtain dynamic or spatial information about the source of Raman scattered light, and when the source (target) moves or fluctuates over time, it is not easy to identify the source or perform an analysis that takes into account the displacement of the source.
- One aspect of the present invention is a system, for example a monitoring system, having a first information acquisition device that acquires first interference information including temporal fluctuations in an interference pattern of Raman scattered light from an object, and a first analysis device that analyzes the first interference information.
- the interference pattern is typically two-dimensional, but may be three-dimensional.
- the first analysis device may include a behavior analysis device that analyzes the behavior of a specific component or an object related thereto in the object based on the first interference information.
- behavior refers to a phenomenon including fluctuations or displacements such as movement, shading, and pulsation of a specific component or an object related thereto. For example, the movement of an object flowing in a flow channel and the accompanying change in concentration of a component (fluctuations in content (content)).
- the first analysis device includes a first pattern analysis device configured to analyze the temporal behavior of a specific component that causes Raman scattered light in the object or an object related thereto (source) based on the asymmetric intensity distribution of the interference pattern.
- first interference information is obtained, which includes two-dimensional or three-dimensional temporal fluctuations of an interference pattern, typically speckles, which were required to be removed or reduced as noise in conventional Raman spectroscopy.
- the interference pattern is generated by the interference and diffraction properties of Raman scattered light, and includes information on the order of wavelength from the source of the Raman scattered light to a sensor (detector, for example, an image element or an image element array) that acquires the interference pattern.
- the Raman scattered light includes information on a specific component that is the cause of the Raman scattered light.
- the inventors of the present application found that the interference pattern has an asymmetric intensity distribution, and that the asymmetric intensity distribution suggests various behaviors of the specific component or an object related to the specific component (for example, an object containing more of the specific component, an object containing less of the specific component, an object in which the specific component fluctuates in association with the behavior, etc.). Therefore, in this system, the behavior of the specific component or an object related to the specific component in the target can be analyzed by analyzing the fluctuations in the asymmetric intensity distribution of the interference pattern of the Raman scattered light over time using the first pattern analysis device.
- the first information acquisition device may be configured to acquire first interference information including an interference pattern of one or more monochromatic Raman scattered lights. Coherent broadband Raman scattered light including one or more peaks reflecting the components contained in the object is generated, and the interference pattern can be acquired. On the other hand, by irradiating the object with laser light of a specific wavelength so as to obtain narrowband Raman scattered light of a specific wavelength of a specific component, it is possible to acquire an interference pattern of monochromatic Raman scattered light of a specific wavelength, and by analyzing this, the behavior of a specific component or an object containing it can be analyzed with greater precision.
- the first pattern analysis device may be configured to analyze the behavior of a specific component or an object related to the component in the target object based on the intensity difference between multiple symmetrical regions (e.g., regions to the left and right, above and below the center of the interference pattern) included in the interference pattern.
- multiple symmetrical regions e.g., regions to the left and right, above and below the center of the interference pattern
- the first analysis device may include a second pattern analysis device configured to determine the concentration of a predetermined component or an object associated with the behavior in the object based on the intensity of the center of the interference pattern and a central region around it. This second pattern analysis device may be configured to determine the concentration of a predetermined component or an object associated with the behavior in the object based on the intensity difference in the central region of the interference pattern of the multiple monochromatic Raman scattered lights included in the first interference information.
- the first analysis device may include a third pattern analysis device configured to obtain physical property values (e.g., pulse rate, flow rate, hematocrit value, etc.) associated with the behavior of the predetermined component or an object associated with the behavior in the object.
- the first analysis device can analyze the change in intensity over time in a specific area of the interference pattern to obtain information about the source of Raman scattered light and its surroundings over time. For example, information such as the movement of the source or an object that includes the source, and changes in concentration of a specific component contained in the source can be obtained.
- the first information acquisition device may be configured to acquire first interference information including an interference pattern of Raman scattered light by at least two laser lights. For example, by irradiating Stokes light and pump light, Raman scattered light can be induced, and the first information including the interference pattern can be acquired efficiently. Furthermore, the location in the depth direction where the spot is formed by at least two laser lights can be controlled. Therefore, an interference pattern that reflects the state of the depth where the spot is formed can be acquired, rather than an interference pattern that mainly reflects the state of the surface of the object.
- the temporal behavior can be analyzed from the first interference information including the interference pattern to identify the source of the interference pattern, and the specific component can be qualitatively analyzed because it is Raman scattered light from that source.
- the first information acquisition device may acquire first interference information including an interference pattern of scattered light of any of stimulated Raman scattering (SRS), Raman amplification, inverse Raman, coherent anti-Stokes Raman scattering (CARS), coherent Stokes Raman scattering (CSRS), and hyper-Raman scattering.
- SRS stimulated Raman scattering
- Raman amplification Raman amplification
- inverse Raman coherent anti-Stokes Raman scattering
- CARS coherent anti-Stokes Raman scattering
- CSRS coherent Stokes Raman scattering
- hyper-Raman scattering hyper-Raman scattering
- the system may have a second information acquisition device configured to acquire second interference information including temporal fluctuations in an interference pattern from the object of a monochromatic laser light that is irradiated onto the object to generate Raman scattered light, and a second analysis device configured to analyze the behavior of the object based on the second interference information.
- the interference pattern of the laser light may be two-dimensional or three-dimensional. An interference pattern of the laser light (laser speckle) can be obtained together with the interference pattern of the Raman scattered light (Raman speckle), and interference information including a wider range of information including the surface of the object can be obtained regardless of the source of the Raman scattered light.
- the system may have a third information acquisition device configured to acquire an image of the surface of the object irradiated with the laser light to generate Raman scattered light, and an image analysis device configured to analyze the irradiation position of the laser light from the image.
- the system may have a generation device configured to generate first interference information including an interference pattern by irradiating the object with the laser light to generate Raman scattered light and receiving the Raman scattered light from the object.
- An example of a generation device is one that irradiates at least two laser lights onto blood in a blood vessel of a living body as the object, and generates first interference information including an interference pattern of the Raman scattered light.
- the first information acquisition device may acquire first interference information including an interference pattern of Raman scattered light in which at least two laser lights are irradiated onto blood in a blood vessel of a living body as the object.
- the predetermined component that is the source of the Raman scattered light may include at least one of glucose, hemoglobin A0, hemoglobin A1c, glycoalbumin, albumin, anhydroglucitol, fructosamine, insulin, glucagon, creatinine, and albumin.
- the first analytical device may analyze the temporal behavior of plasma components and/or blood cell components in blood.
- the system may also have an injection device that injects a drug into the living body based on the state of the living body obtained by the first analytical device.
- the monitoring system includes an information acquisition device that acquires information about Raman scattered light from an object, and an analysis device that analyzes first interference information.
- the method includes the information acquisition device acquiring first interference information including temporal fluctuations in an interference pattern of the Raman scattered light, and the analysis device analyzing the temporal behavior of a predetermined component that causes Raman scattered light in the object or an object related thereto, based on the asymmetric intensity distribution of the interference pattern.
- Acquiring the first interference information may include acquiring first interference information including an interference pattern of at least one monochromatic Raman scattered light.
- Acquiring the first interference information may include acquiring first interference information including an interference pattern of Raman scattered light caused by at least two laser lights.
- the analyzing may include analyzing the behavior of a predetermined component or an object related thereto in the object over time based on the intensity difference of a predetermined region included in the interference pattern included in the first interference information.
- the analyzing may include analyzing the behavior of a predetermined component or an object containing the predetermined component in the object over time based on the intensity difference of a plurality of regions included in the interference pattern included in the first interference information.
- the analyzing may include analyzing the behavior of a predetermined component or an object related thereto in the object based on the intensity difference of a plurality of symmetrical regions included in the interference pattern.
- the analyzing may determine the concentration of the predetermined component or an object related thereto in the object associated with the behavior based on the intensity of the center and its periphery of the interference pattern.
- the determining the concentration may include determining the concentration of the predetermined component or an object related thereto in the object based on the intensity difference of the central region included in the interference pattern of a plurality of monochromatic Raman scattered lights included in the first interference information.
- the method may include an information acquisition device acquiring second interference information including temporal fluctuations in an interference pattern from the object of monochromatic laser light irradiated to the object to generate Raman scattered light, and an analysis device analyzing the behavior of the object based on the second interference information.
- Acquiring the first interference information may include acquiring first interference information including an interference pattern of Raman scattered light when blood in a blood vessel of a living organism is irradiated with at least two laser lights as the object.
- the method may further include injecting a drug into the living organism based on the state of the living organism obtained by the analysis.
- a system having a first information acquisition device that acquires first interference information including an interference pattern and a first analysis device that analyzes the first interference information may be implemented in a computer, and a program (program product) including instructions for causing a computer to function as the above-mentioned system or the first analysis device is included in the present invention, and the program may be provided recorded on an appropriate computer-readable recording medium.
- FIG. 11 shows an example in which the estimated blood glucose value obtained in this example is compared with the actual measured value using a consensus error grid.
- FIG. 12 shows a different example of an image element device.
- FIG. 13 shows an overview of different examples of systems for acquiring CARS speckle.
- FIG. 14 shows an example of a monitoring system.
- FIG. 15 shows an example of measurement using a monitoring system.
- FIG. 16 is a flowchart showing a process of measurement using the monitoring system.
- This system 1 shows an overview of a system (biological monitoring system, monitoring system) 1 that acquires Raman scattered light and analyzes information contained in the Raman scattered light in relation to the present invention.
- This system 1 includes a generating device (detecting device) 10 that irradiates a target with laser light to generate Raman scattered light, receives the Raman scattered light, and generates and provides information (first interference information) 81 that includes an interference pattern (scattering pattern, Raman speckle) 81p of the Raman scattered light, and a management device (monitor) 50 that acquires and analyzes the first interference information 81 and monitors and/or manages the living body 5.
- a generating device (detecting device) 10 that irradiates a target with laser light to generate Raman scattered light, receives the Raman scattered light, and generates and provides information (first interference information) 81 that includes an interference pattern (scattering pattern, Raman speckle) 81p of the Raman scattered light
- a management device (monitor) 50 that
- This system 1 has blood 7 flowing through a blood vessel 6 of a sample (living body) 5 as an observation target (target, monitoring target), and acquires an interference pattern 81p of Raman scattered light 31 derived from blood.
- the generating device (optical device, inspection device) 10 in this example acquires CARS (Coherent Anti-Stokes Raman Scattering) light 31 as Raman scattered light.
- the biological monitoring system 1 may also include a medication system (injection device) 90 that injects medicine to maintain the health of the living body 5.
- An example of a biomonitoring system 1 is a wearable mobile terminal with built-in communication functions and a user interface, such as a smart watch.
- An example of a biomonitoring system 1 may be a desktop or laptop type terminal with built-in communication functions and a user interface.
- An example of a medication system 90 is a system that injects a drug through the skin 5a of a living body 5, and may include an injector 91, a supply device (supply unit) 92 that supplies a specific drug to the injector 91, and a controller 93 that controls them.
- a biological monitoring system (biological management system) 1 is a measurement system (blood glucose measurement device) that measures blood glucose levels.
- the biological monitoring system 1 includes a generating device 10 that includes a Raman optical device (optical system) 19 that irradiates blood 7 flowing through a blood vessel 6 with laser light to generate CARS light 31.
- the generating device 10 shown in FIG. 1 is an example, and shows a generating device that observes a part of a finger of a living body (patient, user) 5, for example, blood 7 flowing through a blood vessel 6 in a nail bed as shown in image 83 in FIG. 2.
- the generating device 10 is not limited to this example as long as it can non-invasively irradiate blood 7 flowing through a blood vessel 6 through the skin of the living body 5 with laser light to generate information including Raman scattered light (CARS light) 31.
- the generating device 10 may be equipped with a minimally invasive method such as controlling the light path with an implant embedded in the living body or forming a blood flow just under the skin by embedding an artificial blood vessel (bioport) in the living body.
- the example generating device 10 includes a laser source 11 for generating (generating) pump light (first laser light) 35 and Stokes light (second laser light) 36.
- An example of the pump light 35 is a monochromatic laser (narrowband laser) with a wavelength of 1030 nm.
- An example of the Stokes light 36 is one or more monochromatic lasers (narrowband lasers) with a predetermined wavelength within the wavelength range of 1100-1300 nm.
- the Stokes light 36 may be generated by a tunable laser, or a laser with a desired wavelength may be selected from broadband lasers with wavelengths of about 1100-1300 nm using a bandpass filter (BP).
- BP bandpass filter
- the wavelength bands of the pump light 35 and Stokes light 36 irradiated to the target to generate the Raman scattered light 31 are not limited to the above.
- a probe light may be irradiated to the target, for example, a blood vessel 6 of a living body 5, to generate the Raman scattered light 31.
- FIG 3 shows an example of a generating device 10 that generates CARS light 31 by irradiating a target with pump light 35 and Stokes light 36 from a laser source 11.
- monochromatic pump light 35 of 1030 nm and two types of monochromatic Stokes light 36 with wavelengths of 1155 nm and 1159 nm selected from broadband laser light 36b are supplied.
- These pump lights 35 and Stokes lights 36 are irradiated onto blood 7 in a blood vessel 6, which is the target, to generate CARS light 31 whose source is blood 7, which is detected (acquired) by the generating device 10.
- CARS light 31 with a wavelength of 928.7 nm corresponding to the Stokes light 36 with a wavelength of 1155 nm and CARS light 31 with a wavelength of 925.8 nm corresponding to the Stokes light 36 with a wavelength of 1159 nm are generated and are received by the two-dimensional image sensor array (EMCCD) 13 in a time-division manner.
- the method of acquiring a plurality of monochromatic CARS lights 31 is not limited to time-division, and as described later, the image sensor array 13 may be divided to simultaneously acquire CARS lights 31 with different wavelengths (in parallel, in parallel).
- the image sensor array 13 of the generating device 10 generates the first interference information 81 including the interference pattern 81p of the CARS light 31 of each wavelength.
- the interference pattern 81p includes the laser speckle of the CARS light 31 of each wavelength.
- the interference pattern 81p including the laser speckle (Raman speckle, CARS speckle) of the CARS light 31 is generated by the interference and diffraction of the CARS light 31, which is Raman scattered light, and includes information on the order of wavelength from the source of the CARS light 31 to the image sensor array 13, which is a sensor that acquires the interference pattern 81p.
- an example of the source of the CARS light 31 is a component contained in the blood 7 in the blood vessel 6.
- the source (location of generation) of the CARS light 31 is inside the living body. Therefore, the interference pattern 81p obtained from the CARS light 31 can be used as one of those called biospeckles.
- the first interference information 81 including the interference pattern 81p of the CARS light 31 will be referred to as CARS speckle 81.
- the CARS speckle 81 may include information on the interference pattern 81p as a speckle image, or may be information converted or compressed into data such as the average brightness (average intensity) of a predetermined area (e.g., the central area) of the speckle image (interference pattern) 81p, or the brightness difference (intensity difference) of multiple predetermined areas (e.g., multiple areas symmetrical around the center, left and right, top and bottom) of the speckle image 81p.
- a predetermined area e.g., the central area
- the speckle image (interference pattern) 81p the brightness difference (intensity difference) of multiple predetermined areas (e.g., multiple areas symmetrical around the center, left and right, top and bottom) of the speckle image 81p.
- FIG 4 shows an example of a spectrum (CARS spectrum) 84 obtained by dispersing CARS scattered light 31 obtained using broadband Stokes light.
- the CARS spectrum shown in Figure 4(a) is an example of a CARS spectrum 84a of a glucose solution (20% solution) and a CARS spectrum 84b of water.
- the CARS spectrum 84a of the glucose solution has a peak Pg specific to glucose that has a top near a wavelength of 928.7 nm and a bottom near a wavelength of 925.8 nm.
- the wavelength of the Stokes light 36 for generating the CARS light 31a with a wavelength of 928.7 nm is 1155 nm
- the wavelength of the Stokes light 36 for generating the CARS light 31b with a wavelength of 925.8 nm is 1159 nm.
- the generating device 10 of the monitoring system 1 shown in FIG. 1 irradiates the pump light 35 and Stokes light 36 supplied from the laser source 11 onto the object, in this example, the skin 5a of a living body (finger) 5 in which blood vessels 6 exist, and focuses the pump light 35 and Stokes light 36 on the blood vessels 6 through the skin 5a.
- the generating device 10 receives the CARS light (Raman scattered light) 31 generated by the blood 7 in the blood vessels 6 with the image sensor 13, and is configured to generate CARS speckles 81.
- the generating device 10 includes a supply optical system (laser light supply optical path, fifth optical path) 25 that guides the pump light 35 and Stokes light 36 supplied from the laser source 11 to a lens (objective lens) 12 facing the object, and a first optical path (first interference pattern acquisition optical system) 21 that guides the CARS light 31 generated by the object and input through the objective lens 12 to a first image sensor array (EMCCD) 13.
- a supply optical system laser light supply optical path, fifth optical path
- first optical path first interference pattern acquisition optical system
- An example of the objective lens 12 is an immersion objective lens, such as a silicone immersion objective lens (silicone oil immersion objective lens), which allows for a larger numerical aperture and a finer resolution to be obtained.
- the objective lens 12 may also be an air immersion lens (air objective lens), which allows for a higher intensity signal to be obtained.
- the generating device 10 further includes a second optical path (second interference pattern acquisition optical system) 22 that guides reflected light (scattered light, second scattered light) 32 of the pump light 35 reflected from the skin 5a (at or near the surface of the skin 5a) via the objective lens 12 to a second image sensor array (CCD) 14, and a third optical path (image acquisition optical system) 23 that guides visible light 33 reflected from the surface of the skin 5a where the object is present via the objective lens 12 to a camera 15.
- the generating device 10 may also include a fourth optical path (spectrum acquisition optical system) 24 that guides the CARS light 34 to a spectroscope (spectrometer) 16 for spectroscopic analysis.
- Each of the optical paths 21 to 25 may include an appropriate optical element corresponding to the purpose, such as a prism, a mirror, or a filter.
- a dichroic mirror 25a that reflects light with a wavelength of 980 nm or less may be provided upstream of the objective lens 12 to separate the short-wavelength CARS light 31 and the visible light image 33 input from the objective lens 12 into the first optical path 21 and the third optical path 23.
- a dichroic mirror 21a that reflects light with a wavelength of 805 nm or less may be provided at the branch point between the first optical path 21 and the third optical path 23 to separate the visible light image 33 into the third optical path 23.
- a filter group 21b that allows light with a wavelength of 900 to 960 nm to pass may be provided upstream of the image sensor array 13 that captures the CARS speckles 81.
- the optical path 25 that guides the pump light 35 and the Stokes light 36 to the objective lens 12 may be provided with a beam splitter 25b for guiding the scattered light 32 of the pump light 35 input from the objective lens 12 to the second optical path 22.
- the generating device 10 may be equipped with a scanning device (scanner) 17 that can adjust the spot positions (irradiation positions) of the pump light 35 and the Stokes light 36 in three dimensions by manipulating the objective lens 12.
- the scanner 17 may include a function of automatic control based on the visible light image 33 obtained by the camera 15. It may also be controlled based on the analysis results of the management device 50 described below or an automatic sequence (protocol) for measurement.
- the scanner 17 may also be equipped with a function of scanning in the depth direction (Z direction) so that the spot positions of the pump light 35 and the Stokes light 36 can be adjusted to the position of the blood vessel 6 under the skin.
- an interference pattern of CARS light (CARS speckle) 81p is generated in the first image sensor array 13 of the first optical path 21 to which the CARS light 31 is guided.
- the first interference information 81 including the interference pattern 81p generated by the generating device 10 is sent to the management device (management device, monitor device) 50.
- the scattered light 32 of the pump light 35 from the skin 5a passes through the dichroic mirror 25a and is guided to the second optical path 22, and an interference pattern of the pump light (laser speckle, pump laser speckle) 82p is generated in the second image sensor array 14.
- the second interference information 82 including the interference pattern 82p generated by the generating device 10 is sent to the management device 50.
- a filter 22a for separating light of a wavelength other than the scattered light 32 of the pump light, for example, the Stokes light 36, may be provided upstream of the second image sensor array 14.
- the image (image information) 83 obtained by the visible light camera 15 and the spectral information 84 obtained by the spectrometer 16 generated by the generation device 10 may also be sent to the management device 50.
- FIG. 5 shows examples of some information obtained by the generating device 10.
- FIG. 5(a) shows an example of a visible light image 83, in which a blood vessel 6 and pump light 35 and Stokes light 36 irradiated in accordance with the position of the blood vessel 6 are shown.
- the visible light image 83 it can be confirmed that the pump light 35 and Stokes light 36 are irradiated to form a spot on the blood vessel 6 or in its vicinity.
- FIG. 5(b) shows an example of second interference information (hereinafter, pump laser speckle) 82 including an interference pattern 82p of the scattered light 32 of the pump light 35.
- pump laser speckle second interference information
- the scattered light 32 of the pump light 35 includes information on the surface structure of the skin 5a irradiated with the pump light 35 and the movement of the skin 5a, i.e., information on the relative minute three-dimensional movement (displacement) between the living body (finger) 5 and the pump light 35, and such information is reflected in the pump laser speckle 82. Therefore, by referring to the temporal displacement of the pump laser speckle 82, the relative behavior (movement, displacement) of the source of the CARS light 31, i.e., the spots (irradiation positions) of the pump light 35 and the Stokes light 36, and the finger 5 can be grasped in three dimensions, and this can be one factor in determining whether the spots are being maintained in the desired location.
- Figure 5(c) shows an example of an interference pattern 85 of CARS light obtained by broadband Stokes light.
- Figure 5(d) shows an example of an interference pattern 81p of CARS speckle 81 (81a) obtained by narrowband Stokes light 36 with a wavelength of 1155 nm
- Figure 5(e) shows an example of an interference pattern 81p of CARS speckle 81 (81b) obtained by narrowband Stokes light 36 with a wavelength of 1159 nm.
- the interference pattern 81p in this example was acquired by an EMCCD (1024 x 1024 pixels) at a frame rate of 4.9 ms (202.55 Hz), but the frame rate etc. are not limited to this example.
- the bandwidth of CARS light 31 obtained by narrowband Stokes light 36 is narrow, the interference or diffraction pattern of the CARS speckle 81 appears clearly. Therefore, by referring to the CARS speckle 81 of the CARS light 31 generated by the narrowband Stokes light (monochromatic Stokes light) 36, the behavior of the source of the CARS light 31, in this example the behavior of the composition components in the blood 7 in the blood vessel 6, can be grasped with greater accuracy.
- the CARS speckle 81a shown in FIG. 5(d) corresponds to an interference pattern 81p of CARS light 31 at 928.7 nm, which is the top of the glucose peak in the CARS spectrum
- the CARS speckle 81b shown in FIG. 5(e) corresponds to an interference pattern 81p of CARS light 31 at 925.8 nm, which is the bottom of the glucose peak in the CARS spectrum. Therefore, by obtaining the intensity difference between these CARS speckles 81a and 81b, for example, the difference in the average brightness of a predetermined region (the center and the central region around it) 88 including the center 89 of the interference pattern 81p as shown in FIG. 5(f), information corresponding to the glucose concentration (content rate, amount) can be obtained.
- the inventors of the present application found that the interference pattern 81p shown in this example has an asymmetric intensity distribution, and that the asymmetric intensity distribution suggests various behaviors of a specific component or an object related to that specific component (e.g., an object containing more of the specific component, an object containing less of the specific component, an object whose specific component fluctuates with its behavior, etc.). For example, even if the intensity of the interference pattern 81p is the same, the intensity difference between the left and right regions 87a and 87b may be reversed, which allows information to be obtained regarding the situation of the movement of the source.
- One example is thought to include information regarding the entry and exit of constituent components such as plasma or blood cells in blood vessels into and out of the spot (irradiation position, focus) of the laser that generates Raman scattering light. That is, the movement of blood cells (red blood cells, white blood cells, etc.) and plasma in blood vessels, and their behavior is reflected in the interference pattern 81p, and the behavior of a specific component or an object related to it in the target object can be analyzed by analyzing the fluctuations in the asymmetric intensity distribution of the interference pattern over time.
- constituent components such as plasma or blood cells in blood vessels into and out of the spot (irradiation position, focus) of the laser that generates Raman scattering light. That is, the movement of blood cells (red blood cells, white blood cells, etc.) and plasma in blood vessels, and their behavior is reflected in the interference pattern 81p, and the behavior of a specific component or an object related to it in the target object can be analyzed by analyzing the fluctuations in the asymmetric intensity distribution of the interference pattern over time.
- the system may be configured to analyze the behavior of a specific component in an object or an object related thereto based on the intensity difference between multiple symmetrical regions 87a and 87b (in this example, regions to the left and right of the center 89 of the interference pattern 81p) contained in the interference pattern 81p.
- the multiple symmetrical regions may be upper and lower regions, and the shape of the region may be a square, or may be another polygon, semicircle, sector, etc.
- the size and shape of the region for which the intensity (brightness) of the interference pattern 81p is to be determined may be determined according to the application.
- the management device 50 of the biological monitoring system 1 shown in FIG. 1 may include a laser control device 51 that controls the laser source 11 of the generating device 10, an information acquisition device (input device, input interface) 52 that acquires information such as CARS speckles 81 from the generating device 10, an analysis device 53 that analyzes the acquired information, a database (library) 54 that stores the acquired information, information for analysis, and further programs 54p, a medication control device 95 that controls the medication system 90, and an input/output interface 60 that controls a user interface 61 such as a display or touch panel.
- a laser control device 51 that controls the laser source 11 of the generating device 10
- an information acquisition device (input device, input interface) 52 that acquires information such as CARS speckles 81 from the generating device 10
- an analysis device 53 that analyzes the acquired information
- a database (library) 54 that stores the acquired information, information for analysis, and further programs 54p
- a medication control device 95 that controls the medication system 90
- an input/output interface 60 that controls
- the input/output interface 60 may include a function 60a for communicating with an external device wirelessly or via a wire, and an alarm function 60b for reporting the results of the analysis device 53 and/or an alarm resulting therefrom to the user interface 61 or a predetermined institution (such as a hospital or security service provider).
- An example of the management device 50 is an information processing device such as a computer or server equipped with computer resources such as a CPU and memory, and performs a predetermined function by downloading a program 54p recorded in an appropriate computer-readable recording medium such as a database 54.
- the information acquisition device (input device) 52 includes a first information acquisition device (first information acquisition function, information acquisition unit) 52a that acquires first interference information (CARS speckle) 81 including temporal fluctuations in a two-dimensional or three-dimensional interference pattern 81p of Raman scattered light from an object, a second information acquisition device (second information acquisition function, information acquisition unit) 52b that acquires second interference information (pump laser speckle) 82 including temporal fluctuations in a two-dimensional or three-dimensional interference pattern 82p from an object of monochromatic laser light (pump light in this example) 35 that is irradiated to the object to generate Raman scattered light, a third information acquisition device (third information acquisition function, information acquisition unit) 52c that acquires image information 83 of the surface of the object to which the laser light is irradiated to generate Raman scattered light, and a fourth information acquisition device (fourth information acquisition function, information acquisition unit) 52d that acquires a Raman spectrum 84 obtained by dispersing the Raman scattered light.
- first information acquisition device
- the first interference information 81 in this example may include interference patterns 81a and 81b of two monochromatic Raman scattered lights 31 generated by two monochromatic Stokes lights 36, respectively, as described above.
- the first interference information 81 may include only an interference pattern 81p of one monochromatic Raman scattered light 31.
- it is possible to estimate the glucose concentration by only acquiring an interference pattern 81p of Raman scattered light 31 of a wavelength corresponding to the top of the glucose peak.
- an interference pattern 81p of Raman scattered light 31 of multiple wavelengths, such as the top and bottom of the peak may be acquired.
- the first interference information 81 in this example is a CARS speckle, and includes an interference pattern 81p of Raman scattered light (CARS light) 31 caused by at least two laser lights, the pump light 35 and the Stokes light 36.
- the CARS light 31 may be generated by irradiating a probe light having a wavelength of, for example, about 780 nm in addition to the pump light 35 and the Stokes light 36, or may be TD-CARS light (time-dependent CARS light) generated by a probe light delayed with respect to the pump light 35 and the Stokes light 36.
- TD-CARS light time-dependent CARS light
- CARS speckles 81 are acquired, but the method is not limited to inverse Raman and coherent anti-Stokes Raman scattering (CARS), and may be an interference pattern of Raman scattered light obtained by stimulated Raman scattering (SRS), Raman amplification, coherent Stokes Raman scattering (CSRS), hyper-Raman scattering, etc.
- the interference pattern is generated by a two-dimensional image sensor array such as a CCD, but since the scattered light is not completely collimated light, a three-dimensional interference pattern may be obtained, for example, the scattered light may be split by a beam splitter or the like to acquire multiple interference patterns at different depths. More information related to the interference and/or diffraction contained in the scattered light can be acquired with even greater precision.
- the analysis device 53 includes a first analysis device (CARS speckle analysis device) 56 that analyzes the behavior of a specific component or an object containing the same in an object such as blood based on CARS speckle 81, which is the first interference information; a second analysis device (pump laser speckle analysis device) 57 that analyzes the behavior (relative behavior) of the object based on pump laser speckle 82, which is the second interference information; a third analysis device (image analysis device) 58 that analyzes image information 83; a fourth analysis device (spectrum analysis device) 59 that analyzes CARS spectrum 84; and an automatic analysis device 55 for automatically performing these analyses by the analysis device 53.
- CARS speckle analysis device CARS speckle analysis device
- the first analysis device 56 includes a first pattern analysis device (first behavior analysis device, first pattern analysis function, analysis unit) 56a configured to analyze the temporal behavior of a specific component that causes Raman scattered light in the object or an object related to the specific component, based on the asymmetric intensity distribution of the interference pattern 81p included in the first interference information (CARS speckle) 81.
- the first pattern analysis device 56a may be configured to analyze the behavior of a specific component in the object or an object related to the specific component, based on the intensity difference between multiple symmetrical regions (e.g., left and right regions 87a and 87b) included in the interference pattern 81p.
- the first analysis device 56 may further include a second pattern analysis device 56b (second behavior analysis device, analysis function, analysis unit) configured to determine the concentration of a predetermined component in the object associated with the behavior based on the intensity of the center of the interference pattern 81p and a central region 88 around it.
- the second pattern analysis device 56b may be configured to determine the concentration of a predetermined component in the object or an object associated therewith based on the intensity difference of the central region 88 of the interference patterns 81a and 81b of the multiple monochromatic Raman scattering lights included in the first interference information 81.
- the first analysis device 56 may further include a third pattern analysis device (third behavior analysis device) 56c configured to obtain physical property values related to the behavior of a predetermined component in the object or an object associated therewith.
- FIGS. 6 and 7 show an example of a two-dimensional CARS speckle 81 obtained from a living body using pump light 35 and narrowband Stokes light 36 with a wavelength of 1155 nm.
- An example of a CARS speckle 81 is a bundle of frames each including an interference pattern 81p along the time at which the frames were acquired.
- the frame rate may be in the range of several tens of microseconds to several hundreds of milliseconds, or several hundred microseconds to several tens of milliseconds.
- the frame rate needs to be sufficiently shorter than the movement (displacement) speed of the object whose behavior is to be obtained, and needs to be such that an interference pattern 81p of a measurable intensity can be obtained.
- Figure 6 shows the time change (vertical direction) of intensity (brightness) in the horizontal direction (X-axis, center line) passing through the center of a CARS speckle 81 containing multiple frames.
- Figure 7 shows an enlarged portion of the CARS speckle 81.
- the source the component that appears as a characteristic
- the CARS speckles 81 reflect the temporal behavior and concentration of glucose as a component in the living body and objects related to it, such as tissues, cells, blood cells and/or plasma in blood, and lymphocytes and/or plasma in lymph.
- 6(a) and 7(a) show CARS speckles 81 obtained from subcutaneous tissue, and almost no change in intensity over time is observed.
- 6(b) and 7(b) show CARS speckles 81 obtained from within a blood vessel 6, and it can be seen that the intensity (intensity and brightness) repeatedly changes at regular intervals. Therefore, by obtaining the intensity difference in frame units (over time) of a specified area including the center of the CARS speckle 81 using the first analysis device 56, it is possible to grasp physical properties associated with the behavior (flow) of the blood 7 in the blood vessel 6, which is the subject of this example, such as pulse rate and blood flow velocity.
- the first analysis device 56 which is a CARS speckle analysis device, may include a fourth pattern analysis device 56d that determines whether or not the CARS speckle 81 is a signal from a blood vessel 6 by performing a comprehensive analysis of the CARS speckle 81.
- Figure 8 shows the intensity on the X-axis of CARS speckles 81 derived from glucose in blood.
- Figure 8(a) shows the intensity change on the X-axis of high-intensity frames included in the course of CARS speckles 81 (dynamic speckle image, speckle video, speckle flow), i.e., the intensity change on the X-axis of interference pattern 81p.
- Figure 8(b) shows the intensity change on the X-axis of low-intensity frames included in the course of CARS speckles 81, i.e., the intensity change on the X-axis of interference pattern 81p.
- FIG. 9 shows an example of the results of analyzing CARS speckles 81 derived from glucose in blood in the time direction.
- the vertical axis (Y axis) shows the results of calculating the intensity difference (difference in brightness) by dividing a predetermined region from the center of the CARS speckle 81 into multiple regions, in this example, the left and right regions 87a and 87b shown in FIG. 5(f)
- the horizontal axis (X axis) shows the average intensity of a predetermined region corresponding to almost the entire CARS speckle 81, for example the intensity (brightness) of the central region 88
- the Z axis shows the time elapse (number of frames), showing the results calculated for each frame plotted.
- the information contained in the CARS speckle 81 fluctuates in a spiral manner over time by processing it as described above. This is assumed to indicate that blood 7 flows through the blood vessel 6 forming the spot (focus) of the pump light 35 and the Stokes light 36, and red blood cells enter and leave the spot.
- the first pattern analyzer 56a can therefore analyze the behavior of red blood cells in blood 7 based on the intensity difference between multiple regions of the CARS speckle 81.
- the second pattern analyzer 56b can extract CARS speckle 81 of the plasma portion suitable for measuring blood glucose levels from the blood-derived CARS speckle 81 and obtain the intensity of its central region 88 to obtain the glucose concentration mainly reflected in the CARS speckle 81.
- the third pattern analyzer 56c can extract frames to which red blood cells are the main contributor and frames to which plasma is the main contributor from among the many frames of the CARS speckle 81, making it possible to obtain, for example, the hematocrit value, which is the ratio of red blood cells to plasma, as a physical property value related to the behavior of red blood cells.
- FIG. 10 shows an example in which CARS speckles 81a and 81b are obtained in sequence, for example, in 3-second intervals (3-second intervals, 3-second time divisions), from narrowband Stokes light 36 with a wavelength of 1155 nm and narrowband Stokes light 36 with a wavelength of 1159 nm.
- the second pattern analyzer 56b extracts a frame reflecting an interference pattern 81p of a plurality of monochromatic Raman scattering lights contained in the CARS speckle 81, in this example, an interference pattern 81p of two-color CARS speckles 81a and 81b at 928.7 nm, which is the top of the glucose peak in the CARS spectrum, and 925.8 nm, which is the bottom of the peak, due to the narrowband Stokes light 36 of each of the above wavelengths, and obtains the intensity difference of each predetermined region, for example, a region 88 including the center.
- This intensity difference corresponds to the glucose concentration in the plasma portion of the blood 7, so that it is possible to obtain quantitative and highly accurate information on the blood glucose level, including changes over time.
- the CARS speckle of the frame corresponding to the plasma portion of the 3D profile shown in FIG. 9 may be selected.
- Figure 11 shows the estimated blood glucose level obtained by the above process using the system 1 of this example, compared with the actual measured value (reference) obtained by a glucose meter (Glucometer (SMBG) Nipro Freestyle Freedom Lite) using a consensus error grid.
- Glucometer Glucometer
- the system 1 of this example can provide a function as a sensor that constantly measures blood glucose levels in real time, a function as a blood glucose monitor, and a function to manage the human body based on blood glucose levels.
- the system 1 can also provide a function to constantly monitor other physical properties (parameters) in the blood, not just blood glucose levels, in real time.
- Figure 12 shows different examples of an image sensor array (image generating device, image sensor) 13 for generating an interference pattern.
- image sensor array image generating device, image sensor
- the image sensor array 13 shown in Figure 12 has a Bayer pattern (Bayer pattern, Bayer arrangement, Bayer filter) as shown in Figure 12 (a), and as shown in Figure 12 (b), three types of pixels or filters 13a to 13c sensitive to three types of light with different wavelengths are arranged in a checkerboard pattern.
- this image sensor array 13 multiple, for example, three types of monochromatic interference patterns can be acquired in parallel.
- the original interference pattern obtained by the image sensor array 13 is a mosaic image.
- an image (interference pattern) for each color can be obtained as shown in Figure 12 (c).
- FIG. 13 shows another example of a generating device 10 that generates CARS light 31 by irradiating a target with pump light 35 and Stokes light 36 from a laser source 11.
- this generating device 10 monochromatic pump light 35 and monochromatic (narrowband) Stokes light 36 of four different wavelengths selected from broadband laser light 36b are irradiated to blood 7 in a blood vessel 6, which is a target, to generate CARS light 31 of four different wavelengths originating from blood 7.
- These CARS lights 31 of different wavelengths are spatially separated by a beam splitter 18 and guided to one or more imaging element arrays (image sensors) 13, and interference patterns of CARS light 31 of each wavelength are generated in parallel (simultaneously).
- the beam splitter 18 may be adjusted so that interference patterns of CARS light 31 of different wavelengths are generated in different areas of one image sensor 13, or may be adjusted so that interference patterns of CARS light 31 of different wavelengths are generated in different image sensors 13.
- the narrowband Stokes light 36 may be generated from the broadband Stokes light 36b using a bandpass filter, a diffraction grating, or the like.
- the narrowband 6 may be generated using a tunable laser.
- FIGS. 14 and 15 show an example of a biomonitoring system (monitoring system, monitor) 1 that measures the blood glucose level in the blood flowing through the blood vessels in the human wrist (more specifically, the inside of the wrist (palm side)) 5.
- this monitoring system 1 is a portable desktop-type device that is approximately 250 mm in length and width and 120 mm in height, and includes a housing 66 that incorporates the generating device 10 and the management device 50, and a cover 65 that also serves as a user interface 61.
- a display 61 with a touch panel function appears as a user interface
- a mechanism 64 for setting the user's wrist 5 appears on the top surface of the housing 66.
- the top surface of the housing 66 includes an arm mounting section 67 that is recessed in the shape of a hand so that the hand (arm) can be placed on it and the wrist can be easily aligned to a predetermined position, a palm rest 68 for holding the wrist in front of the objective lens 12 of the generating device 10, and a wrist strap (wristband) 69 for holding the wrist.
- the user opens the cover 65 of the biomedical management system 1 and places the inside of the wrist 5 on the top surface of the housing 66 facing the objective lens 12 of the generating device 10, thereby enabling the blood glucose level to be measured non-invasively and in real time and displayed on the display 61.
- the biomedical management system 1 may include a medication device 90 in addition to the management device 50 in the housing 66, and may be configured to automatically receive processing according to the measurement results.
- step 101 the management system 1 is installed at a predetermined location on the human body and measurement is started.
- the inside of the wrist 5 is set to face the objective lens 12 of the generation device 10, and then the start button of the management system 1 is clicked.
- step 102 the blood vessel 6 to be measured is automatically selected. This process may be performed automatically (autonomously) by the scanning device 17 of the generation device 10 using the image information 83 obtained by the generation device 10, or may be performed by the image analysis device 58 of the management device 50.
- An example of a blood vessel 6 selected based on the image information 83 is a capillary vessel with a blood vessel diameter of 6 ⁇ m or less and a blood flow velocity of 300 ⁇ m/sec or less that can be observed from a flat portion of the wrist 5.
- step 103 detailed alignment of the laser (laser spot, irradiation position) to be irradiated onto the selected blood vessel 6 is performed.
- This process may also be performed automatically (autonomously) by the scanning device 17 of the generating device 10 using the image information 83, or may be performed by the image analysis device 58 of the management device 50.
- Image information 83 of the pump light 35 and Stokes light 36 actually irradiated onto the skin 5a of the wrist 5 may be obtained, and the laser spot may be aligned with the blood vessel 6.
- step 104 the generating device 10 irradiates the skin 5a of the wrist 5 with pump light 35 and Stokes light 36, generates CARS light 31, and acquires interference information (CARS speckle) 81 including an interference pattern 81p.
- the generating device 10 generates pump laser speckle 82 and image information 83 together with the CARS speckle 81, and supplies them to the managing device 50.
- step 105 the managing device 50 analyzes this information (data) acquired via the acquiring device 52 with the analyzing device 53 to verify the accuracy of the information. Items to be verified include (1) whether the signal (data) is derived from blood, (2) whether any signals other than blood are mixed in, and (3) whether there is positional stability.
- Whether or not a signal is derived from blood may be determined, for example, by the fourth pattern analyzer 56d of the first analyzer 56 comparing the intensity (average intensity) of the central region 88 of the interference pattern 81p contained in the CARS speckle 81 with that observed when tissue is observed, and determining whether or not there is an oscillatory increase or decrease in brightness intensity on a time scale of 10-100 ms.
- Whether or not a signal other than blood is mixed in may be determined, for example, by the first pattern analyzer 56a of the first analyzer 56, by determining whether or not a stable spiral structure is observed in the 3D plot image shown in FIG. 9 after a predetermined time (predetermined number of frames) has elapsed.
- the criterion may be whether or not a stable spiral structure is observed for 1 second or more, at about 10 times/second.
- the positional stability may be determined based on the behavior of the object (skin 5a or human body (wrist)) 5 obtained by analyzing the second interference information (pump laser speckle) 82, which includes temporal fluctuations in the interference pattern 82p from the object (skin 5a) of the monochromatic laser light (pump light 35 in this example) irradiated to the object to generate Raman scattered light, acquired by the information acquisition device 52, using the second analysis device 57. For example, if a fluctuation of 10% or more is observed in the brightness (intensity) of the central region of the interference pattern 82, it may be determined that the object's position is unstable.
- the image information 83 may be analyzed by the image analysis device 58, and it may be determined that the object's position is unstable because vibrations or shifts are observed in the relative positions of the spots of the pump light 35 and/or Stokes light 36 and the blood vessel 6.
- the fourth analyzer (spectral analyzer) 59 of the first analyzer 53 may intermittently or periodically analyze the CARS spectrum 84 to determine whether a predetermined CARS signal, for example a peak or trough corresponding to glucose, has been obtained.
- a predetermined CARS signal for example a peak or trough corresponding to glucose
- the analysis device 53 performs data processing.
- the first analysis device 56 analyzes the first interference information (CARS speckle) 81 including the temporal fluctuation of the interference pattern 81p of the Raman scattered light acquired by the information acquisition device 52.
- the first pattern analysis device 56a analyzes the temporal behavior of objects (blood cells (red blood cells) and plasma in this example) related to a specific component (glucose in this example) that causes Raman scattered light in the blood vessel 6, which is the target, based on the asymmetric intensity distribution of the interference pattern 81p. More specifically, the first pattern analysis device 56a analyzes the behavior of the blood cells and plasma in the blood vessel based on the intensity difference between multiple symmetrical regions (left and right regions) 87a and 87b included in the interference pattern 81p.
- the second pattern analyzer 56b determines the concentration of a predetermined component or an object related to the behavior in the object based on the intensity of the center of the interference pattern 81p and the surrounding region 88.
- the second pattern analyzer 56b determines the value of glucose in the plasma (plasma blood glucose level), which is one of the objects whose behavior has been confirmed in the blood 7, which is the object. More specifically, a frame corresponding to the bottom (the part where the intensity difference between the left and right regions 87a and 87b is small) of the spiral pattern of the 3D plot image shown in FIG. 9 generated by the first pattern analyzer 56a is selected from the CARS speckle 81, and the plasma blood glucose level is determined based on the intensity of the central region 88.
- the second pattern analyzer 56b can determine the plasma blood glucose level based on the intensity difference of the central region 88 included in the interference pattern 81p of multiple monochromatic Raman scattering lights (in this example, the interference pattern of one peak wavelength of glucose and the interference pattern of the bottom wavelength) of CARS speckle 81.
- the second pattern analysis device 56b is not limited to obtaining plasma glucose levels, but can also obtain blood glucose levels of blood cells (red blood cells), and may select multiple frames that are most suitable for calculating the glucose level from the CARS speckles 81 obtained from the generation device (detection device) 10, and obtain the desired glucose level based on those frames.
- the third pattern analyzer 56c of the first analyzer 56 may obtain a physical property value related to the behavior of a predetermined component or an object related thereto in the object, for example, a hematocrit value indicating the ratio of red blood cells to plasma.
- the first analyzer 56 may obtain information on components other than glucose that cause CARS scattered light, that is, components that have a large effect on CARS scattering or substances related thereto, using the interference pattern 81p of multiple monochromatic CARS scattered lights.
- hemoglobin A0, hemoglobin A1c, glycoalbumin, albumin, anhydroglucitol, fructosamine, insulin, and glucagon can be listed as targets (components) for which abnormalities in the body can be detected in relation to blood glucose levels.
- targets (components) for which abnormalities in renal function can be detected include creatinine and albumin.
- Other blood components may be measured and monitored using this management system 1. Additionally, the management system 1 may acquire other information regarding blood vessels or blood, such as blood flow and pulsation (beats), and these may also be subject to monitoring.
- the input/output interface 60 displays the various measured results on the user interface 61.
- the input/output interface 60 may provide the measured results to a specified external system.
- the medication control device 95 determines whether or not medication is required based on the measured results. If the medication control device 95 determines that medication, for example, injection of insulin, is required, in step 109, it causes the medication system 90 to perform a procedure to inject a specified drug into the user.
- the system 1 of this example has a generating device (detecting device) 10 that irradiates at least two laser beams 35 and 36 onto blood 7 in a blood vessel 6 of a living body 5 as an object, and generates first interference information (CARS speckle) 81 including an interference pattern of the Raman scattered light.
- the first analyzing device 56 can also analyze the temporal behavior of plasma components and/or blood cell components in the blood, and in the above example, calculates the blood glucose level from the glucose concentration contained in the plasma components.
- the first analyzing device 56 can analyze the quantitative behavior of components other than glucose, for example, components including at least one of hemoglobin A1c, creatinine, and albumin.
- components contained in other bodily fluids can also be measured in a similar manner.
- the main components are not limited to plasma components and red blood cells, but may also include other blood cell components, for example, white blood cells and/or platelets.
- the target component whose concentration is to be measured may include any of the components that are the subject of testing for bodily fluids such as blood.
- the subject of measurement may be any chemical, molecule, compound, composition, microorganism, or aggregate, including, but not limited to, blood cells, amino acids, peptides, polypeptides, proteins, glycoproteins, lipoproteins, nucleosides, nucleotides, oligonucleotides, nucleic acids, sugars, carbohydrates, oligosaccharides, polysaccharides, fatty acids, lipids, hormones, metabolites, cytokines, chemokines, receptors, neurotransmitters, antigens, allergens, antibodies, substrates, cofactors, inhibitors, drugs, pharmaceuticals, nutrients, prions, toxins, poisons, explosives, pesticides, chemical warfare agents, biological hazards, radioisotopes, vitamins, heterocyclic aromatic compounds, carcinogens, mutagens, narcotics, amphetamines, barbiturates, hallucinogens, waste products, and/or pollutants.
- Microorganisms include, but are not limited
- the Raman scattered light is not limited to CARS light, but may be scattered light obtained using other known methods such as stimulated Raman scattering (SAS) and surface enhanced Raman scattering (SERS), and it is possible to obtain information including their interference patterns.
- SAS stimulated Raman scattering
- SERS surface enhanced Raman scattering
- the system 1 of this example has a first information acquisition device that acquires first interference information including temporal fluctuations in a two-dimensional or three-dimensional interference pattern of Raman scattered light from an object, and a first analysis device that analyzes the behavior of a specific component in the object or an object containing the same based on the first interference information.
- the first analysis device may include a behavior analysis device that analyzes the temporal behavior of a specific component in the object or an object containing the same based on the intensity difference in a specific two-dimensional or three-dimensional region included in the interference pattern included in the first interference information.
- the first analysis device may include a behavior analysis device that analyzes the temporal behavior of a specific component in the target object or an object containing the same, based on the intensity differences between multiple two-dimensional or three-dimensional regions included in the interference pattern contained in the first interference information. For example, even if the intensity of the interference pattern is the same, the intensity difference between the left and right regions may be reversed, thereby making it possible to obtain information on the status of the movement of the source.
- One example is the discovery that information can be obtained on the entry and exit of constituent components such as plasma or blood cells in blood vessels into and out of the laser spot that generates Raman scattered light.
- the first analysis device may include a behavior analysis device that analyzes the quantitative behavior of a specific component in the target or an object containing the same based on the intensity difference in a specific region included in the interference pattern of multiple monochromatic Raman scattered lights included in the first interference information.
- the concentration of the specific component and its fluctuation can be obtained quantitatively with greater accuracy.
- One example is obtaining an interference pattern of Raman scattered light at the wavelengths of the top and bottom of the peak of the specific component included in the Raman spectrum.
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Optics & Photonics (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
監視システム(1)は、対象物(6)からラマン散乱光の干渉パターン(81p)の時間的な変動を含む第1の干渉情報(81)を取得するように構成された第1の情報取得装置(52a)と、第1の干渉情報を分析する第1の分析装置(56)とを有し、第1の分析装置は、干渉パターンの非対称な強度分布に基づき、対象物内の前記ラマン散乱光の発生要因となる所定の成分またはそれに関連する物体の時間的な挙動を分析するように構成された第1のパターン解析装置(56a)を含む。
Description
本発明は、ラマン散乱光を用いたシステムおよび方法に関するものである。
国際公開WO2011/071921(日本国特表2013-517516)号公報には、装置および方法は、光ファイバにおける誘導ラマン散乱を使用することでレーザースペックルを低減することが開示されている。ファイバコア径および長さは所望の出力光が得られるように選択され、2つの光ファイバを並列に合成することによって、および、ビームスプリッタとしての回転可能波長板および偏光器の助けにより各経路の光量を調整することによって、調整可能デスペックラが形成される。
国際公開WO2009/008745(日本国特表2010-532699)号公報には、レーザースペックル画像化システムおよび方法が開示されている。この文献においては、組織内の灌流を計測するための装置及び方法を提供し、その方法は、レーザー光下の組織の画像を記録し、組織の複数の画像から複数のコントラスト画像を計算し、複数のコントラスト画像から散乱光のパワースペクトルを決定し、パワースペクトルから灌流を決定するステップを含む。
生体の体液、例えば、血液には主な組成成分として血球と血漿とが含まれ、血球としては、赤血球、白血球、血小板などが含まれる。また、生体の健康状態をチェックするためには、主な組成成分とは別に、血液中のグルコース、ヘモグロビンA1c、クレアチニン、アルブミンなどの測定が求められ、血糖値の測定基準としては血漿血糖値が要求される。したがって、例えば、生体からグルコースの濃度に関する情報が得られたとしても、その情報が何に由来するものかが判断可能か否かにより情報としての価値または評価が異なることがある。
ラマン分光法は、ラマン散乱光を用いて物質の評価を行う分光法(スペクトロスコピー)であり、光を物質に照射すると、光が物質と相互作用することで入射光と異なる波長を持つラマン散乱光が出射され、その強度を周波数または波長の関数として示すスペクトルを取得する。そのスペクトルに含まれるピーク位置からは化学結合の情報が得られ、スペクトル全体の波形からは分子構造の情報などを得ることができる。しかしながら、ラマン散乱光の発生源の動的な情報あるいは空間的な情報を同時に得ることは容易ではなく、発生源(ターゲット)が動いたり、時間的に変動したりした場合に、発生源の特定や、発生源の変位などを加味した分析は容易ではない。例えば、血管中を流れる血液に含まれる特定の組成成分に含まれる所定の成分の濃度を精度よく測定することは難しい。生体などから得られるラマン散乱光を用いて、対象物、またはそれに含まれる成分の挙動をさらに精度よく分析できるシステムおよび方法が求められている。
本発明の一態様は、対象物からラマン散乱光の干渉パターンの時間的な変動を含む第1の干渉情報を取得する第1の情報取得装置と、第1の干渉情報を分析する第1の分析装置とを有するシステム、例えば、モニタリグシステムである。干渉パターンは典型的には2次元であるが、3次元であってもよい。第1の分析装置は、第1の干渉情報に基づき対象物内の所定の成分またはそれに関連する物体の挙動を分析する挙動解析装置を含んでもよい。本明細書において挙動は所定の成分またはそれに関連する物体の移動、濃淡、脈動などの変動または変位を含む現象を示す。例えば、流路内を流れる物体の移動や、それに伴う成分の濃度変化(含有量(含有率)の変動)である。挙動の一例は、血球(赤血球、白血球など)および血漿の血管内の移動である。第1の分析装置は、干渉パターンの非対称な強度分布に基づき、対象物内のラマン散乱光の発生要因となる所定の成分またはそれに関連する物体(発生源)の時間的な挙動を分析するように構成された第1のパターン解析装置を含む。
このシステムにおいては、従来のラマン分光法において、ノイズとして除去または低減することが求められていた干渉パターン、典型的にはスペックルの2次元または3次元の時間的な変動を含む第1の干渉情報を取得する。干渉パターンは、ラマン散乱光の干渉・回折性により生じ、ラマン散乱光の発生源から干渉パターンを取得するセンサー(ディテクタ、例えば画像素子、画像素子アレイ)に至るまでの波長オーダーの情報を含む。そして、ラマン散乱光は、ラマン散乱光の発生要因となる所定の成分の情報を含む。さらに、本願の発明者らは、その干渉パターンに非対称な強度分布があることを見出し、その非対称な強度分布が、所定の成分またはその所定の成分に関連する物体(例えば、所定の成分をより多く含む物体、所定の成分をより少なく含む物体、挙動に付随して所定の成分が変動する物体など)の様々な挙動を示唆するものであることを見出した。したがって、本システムにおいては、第1のパターン解析装置により、ラマン散乱光の干渉パターンの非対称な強度分布の変動を時間に即して分析することにより、対象物内の所定の成分またはそれに関連する物体の挙動を分析できる。
第1の情報取得装置は、1つまたは複数の単色のラマン散乱光の干渉パターンを含む第1の干渉情報を取得するように構成されていてもよい。コヒーレントなブロードバンドにより対象物に含まれる成分を反映した1つまたは複数のピークを含むラマン散乱光が生成され、その干渉パターンを取得できる。一方、特定の成分の特定の波長のナローバンドのラマン散乱光が得られるように特定の波長のレーザー光で対象物を照射することにより、特定の波長の単色のラマン散乱光の干渉パターンを取得することが可能であり、それを解析することにより、特定の成分のまたはそれを含む物体の挙動をより精度よく分析できる。
第1のパターン解析装置は、干渉パターンに含まれる対称な複数の領域(例えば、干渉パターンの中心に対して左右、上下などの領域)の強度差に基づき、対象物内の所定の成分またはそれに関連する物体の挙動を解析するように構成されていてもよい。
第1の分析装置は、干渉パターンの中央およびその周辺の中央領域の強度に基づき、挙動に関連付けられた対象物内の所定の成分またはそれに関連する物体の濃度を求めるように構成された第2のパターン解析装置を含んでもよい。この第2のパターン解析装置は、第1の干渉情報に含まれる複数の単色のラマン散乱光の干渉パターンの中央領域の強度差分に基づき、対象物内の所定の成分またはそれに関連する物体の濃度を求めるように構成されていてもよい。第1の分析装置は、対象物内の所定の成分またはそれに関連する物体の挙動に関連する物性値(例えば、脈拍、流速、ヘマトクリット値など)を得るように構成された第3のパターン解析装置を含んでもよい。
このシステムにおいては、第1の分析装置により、干渉パターンの所定の領域の強度の時間変化を分析(解析)することによりラマン散乱光の発光源およびその周囲の時間的な変動に関する情報を得ることができる。例えば、発生源または発生源を含む物体の動き(移動)、発生源に含まれる所定の成分の濃度変化などの情報を得ることができる。
第1の情報取得装置は、少なくとも2つのレーザー光によるラマン散乱光の干渉パターンを含む第1の干渉情報を取得するように構成されていてもよい。例えば、ストークス光と、ポンプ光とを照射することによりラマン散乱光を誘導することができ、効率よく干渉パターンを含む第1の情報を取得できる。さらに、少なくとも2つのレーザー光によりスポットを形成する深度方向の場所を制御できる。したがって、対象物の表面の状態を主に反映した干渉パターンではなく、スポットが形成される深度の状態を反映した干渉パターンを取得できる。例えば、皮膚下の血管を流れる血液にスポットを形成するように少なくとも2つのレーザー光を照射することにより、血液中の特定の成分の有無あるいは濃度のみならず、血液に含まれる特定の成分または、その特定の成分を含む血漿および血球に関する時間的な挙動を分析できる。あるいは、干渉パターンを含む第1の干渉情報から時間的な挙動を分析して、干渉パターンの発生源を特定し、その発生源からのラマン散乱光であることから、特定の成分を定性的に解析できる。第1の情報取得装置は、誘導ラマン散乱(SRS)、ラマン増幅、逆ラマン、コヒーレント・反ストークスラマン散乱(CARS)、コヒーレント・ストークスラマン散乱(CSRS)、ハイパーラマン散乱のいずれかの散乱光の干渉パターンを含む第1の干渉情報を取得してもよい。
このシステムは、ラマン散乱光を発生させるために対象物に照射される単色のレーザー光の対象物からの干渉パターンの時間的な変動を含む第2の干渉情報を取得するように構成された第2の情報取得装置と、第2の干渉情報に基づき、対象物の挙動を分析するように構成された第2の分析装置とを有していてもよい。レーザー光の干渉パターンは2次元または3次元であってもよい。ラマン散乱光の干渉パターン(ラマンスペックル)とともに、レーザー光の干渉パターン(レーザースペックル)を得ることができ、ラマン散乱光の発生源に関わらず、対象物の表面を含めたより広い範囲の情報を含む干渉情報を得ることができる。
このシステムは、ラマン散乱光を発生させるためにレーザー光が照射される対象物の表面の画像を取得するように構成された第3の情報取得装置と、画像からレーザー光の照射位置を解析するように構成された画像解析装置とを有していてもよい。このシステムは、対象物にラマン散乱光を発生させるためのレーザー光を照射して対象物からのラマン散乱光を受光することにより干渉パターンを含む第1の干渉情報を生成するように構成された生成装置を有してもよい。生成装置の一例は、対象物として、生体の血管中の血液に少なくとも2つのレーザー光を照射し、そのラマン散乱光の干渉パターンを含む第1の干渉情報を生成するものである。第1の情報取得装置は、対象物として、生体の血管中の血液に少なくとも2つのレーザー光が照射されたラマン散乱光の干渉パターンを含む第1の干渉情報を取得してもよい。ラマン散乱光の発生源となる所定の成分は、グルコース、ヘモグロビンA0、ヘモグロビンA1c、グリコアルブミン、アルブミン、アンヒドログルシトール、フルクトサミン、インシュリン、グルカゴン、クレアチニンおよびアルブミンの少なくともいずれかを含んでもよい。第1の分析装置は、血液中の、血漿成分および/または血球成分の時間的な挙動を分析してもよい。このシステムは、第1の分析装置により得られる生体の状態に基づいて、生体に薬剤を注入する注入装置を有していてもよい。
本発明の他の態様の1つは、モニタリングシステムの制御方法である。モニタリングシステムは、対象物からラマン散乱光の情報を取得する情報取得装置と、第1の干渉情報を分析する分析装置とを有する。この方法は、情報取得装置が、ラマン散乱光の干渉パターンの時間的な変動を含む第1の干渉情報を取得することと、分析装置が、干渉パターンの非対称な強度分布に基づき、対象物内のラマン散乱光の発生要因となる所定の成分またはそれに関連する物体の時間的な挙動を分析することを有する。第1の干渉情報を取得することは、少なくとも1つの単色のラマン散乱光の干渉パターンを含む第1の干渉情報を取得することを含んでもよい。第1の干渉情報を取得することは、少なくとも2つのレーザー光によるラマン散乱光の干渉パターンを含む第1の干渉情報を取得することを含んでもよい。
分析することは、第1の干渉情報に含まれる干渉パターンに含まれる所定の領域の強度差に基づき、対象物内の所定の成分またはそれに関連する物体の時間的な挙動を解析することを含んでもよい。分析することは、第1の干渉情報に含まれる干渉パターンに含まれる複数の領域の強度差分に基づき、対象物内の所定の成分またはそれを含む物体の時間的な挙動を解析することを含んでもよい。分析することは、干渉パターンに含まれる対称な複数の領域の強度差に基づき、対象物内の所定の成分またはそれに関連する物体の挙動を解析することを含んでもよい。分析することは、干渉パターンの中央およびその周辺の強度に基づき、挙動に関連付けられた対象物内の所定の成分またはそれに関連する物体の濃度を求めてもよい。濃度を求めることは、第1の干渉情報に含まれる複数の単色のラマン散乱光の干渉パターンに含まれる、中央領域の強度差分に基づき、対象物内の所定の成分またはそれに関連する物体の濃度を求めることを含んでいてもよい。
当該方法は、情報取得装置が、ラマン散乱光を発生させるために対象物に照射される単色のレーザー光の対象物からの干渉パターンの時間的な変動を含む第2の干渉情報を取得することと、分析装置が、第2の干渉情報に基づき、対象物の挙動を分析することとを含んでいてもよい。第1の干渉情報を取得することは、対象物として、生体の血管中の血液に少なくとも2つのレーザー光を照射されたラマン散乱光の干渉パターンを含む第1の干渉情報を取得することを含んでいてもよい。当該方法は、さらに、分析することにより得られる生体の状態に基づいて、生体に薬剤を注入することを有していてもよい。
干渉パターンを含む第1の干渉情報を取得する第1の情報取得装置と、その第1の干渉情報を分析する第1の分析装置とを有するシステム、特に、第1の分析装置はコンピュータに実装されていてもよく、コンピュータを上記のシステムまたは第1の分析装置として機能させるための命令を含むプログラム(プログラム製品)は本発明に含まれ、プログラムは、コンピュータに読み取り可能な適当な記録媒体に記録して提供されてもよい。
図1に本発明に関連して、ラマン散乱光を取得し、そのラマン散乱光に含まれる情報を解析するシステム(生体監視システム、モニタリングシステム)1の概要を示している。このシステム1は、対象物にレーザー光を照射してラマン散乱光を生成し、そのラマン散乱光を受光してラマン散乱光の干渉パターン(散乱パターン、ラマンスペックル)81pを含む情報(第1の干渉情報)81を生成および提供する生成装置(検出装置)10と、第1の干渉情報81を取得して分析し、生体5を監視および/または管理する管理装置(モニター)50とを含む。このシステム1は、サンプル(生体)5の血管6を流れる血液7を観察対象(対象物、モニタリングの対象)とし、血液由来のラマン散乱光31の干渉パターン81pを取得する。本例の生成装置(光学装置、検査装置)10は、ラマン散乱光としてCARS(コヒーレント反ストークスラマン散乱、Coherent Anti-Stokes Raman Scattering)光31を取得する。生体監視システム1は、生体5の健康を維持するための薬剤を注入する投薬システム(注入装置)90を備えていてもよい。
生体監視システム1の一例はスマートウォッチなどの通信機能およびユーザーインターフェイスを内蔵したウェアラブルな携帯端末である。生体監視システム1の一例は、通信機能およびユーザーインターフェイスを内蔵したデスクトップまたはラップトップ型の端末であってもよい。投薬システム90の一例は、生体5の皮膚5aを介して薬剤を注入するシステムであり、インジェクタ91と、インジェクタ91に所定の薬剤を供給する供給装置(供給ユニット)92と、それらを制御するコントローラ93とを含んでいてもよい。
生体監視システム(生体管理システム)1の一例は血糖値を測定する測定システム(血糖値測定装置)である。生体監視システム1は、血管6を流れる血液7にレーザー光を照射してCARS光31を生成するラマン光学装置(光学系)19を含む生成装置10を含む。図1に示した生成装置10は一例であり、生体(患者、ユーザー)5の指の一部、例えば、図2に画像83で示すような、爪床の血管6を流れる血液7を観察対象とした生成装置を示している。生成装置10は、非侵襲で生体5の皮膚を介して血管6を流れる血液7にレーザー光を照射してラマン散乱光(CARS光)31を含む情報を生成できればよく、この例に限定されるものではない。また、生成装置10は、生体内に埋め込まれたインプラントで光路を制御したり、生体内に人口血管(バイオポート)を埋め込んで血液の流れを皮膚直下に形成したりするなどの低侵襲的な手法を備えたものであってもよい。
一例の生成装置10は、ポンプ光(第1のレーザー光)35およびストークス光(第2のレーザー光)36を発生(生成)するためのレーザー源11を含む。ポンプ光35の一例は、波長1030nmの単色レーザー(ナローバンドレーザー)である。ストークス光36の一例は、波長1100-1300nm内の所定の波長の1つまたは複数の単色レーザー(ナローバンドレーザー)である。ストークス光36はチューナブルレーザーにより生成されてもよく、波長1100-1300nm程度のブロードバンドレーザーからバンドパスフィルター(BP)により所望の波長のレーザーを選択してもよい。ラマン散乱光31を生成するためにターゲットに照射されるポンプ光35およびストークス光36の波長帯は上記に限定されない。また、ラマン散乱光31を生成するために、ポンプ光35およびストークス光36に加えてプローブ光をターゲット、例えば、生体5の血管6に向けて照射してもよい。
図3に、レーザー源11からポンプ光35およびストークス光36をターゲットに照射してCARS光31を生成する生成装置10の一例を示している。この例では、1030nmの単色のポンプ光35と、ブロードバンドのレーザー光36bから選択された波長1155nmおよび波長1159nmの2種類の単色のストークス光36とが供給される。これらのポンプ光35およびストークス光36は、ターゲットである血管6中の血液7に照射され、血液7を発生源とするCARS光31が生成され、生成装置10はそれを検出(取得)する。波長1155nmのストークス光36および波長1159nmのストークス光36を順番に照射することにより、波長1155nmのストークス光36に対応した波長928.7nmを中心とするCARS光31と、波長1159nmのストークス光36に対応した波長925.8nmのCARS光31が生成され、時分割で2次元の撮像素子アレイ(EMCCD)13により受光される。複数の単色のCARS光31を取得する方法は時分割に限らず、後述するように、撮像素子アレイ13を分割して波長の異なるCARS光31を同時に(平行して、並列に)取得するようにしてもよい。
本例においては、生成装置10の撮像素子アレイ13により、それぞれの波長のCARS光31の干渉パターン81pを含む第1の干渉情報81が生成される。干渉パターン81pは、それぞれの波長のCARS光31のレーザースペックルを含む。CARS光31のレーザースペックル(ラマンスペックル、CARSスペックル)を含む干渉パターン81pはラマン散乱光であるCARS光31の干渉・回折性により生じ、CARS光31の発生源から干渉パターン81pを取得するセンサーである撮像素子アレイ13に至るまでの波長オーダーの情報を含む。本例において、CARS光31の発生源の一例は血管6中の血液7に含まれる成分である。したがって、CARS光31の発生源(発生個所)は生体内部となる。このため、CARS光31から得られる干渉パターン81pはバイオスペックルと称されるものの1つとして用いることができる。なお、以降において、CARS光31の干渉パターン81pを含む第1の干渉情報81をCARSスペックル81として参照する。CARSスペックル81は、干渉パターン81pの情報をスペックル画像として備えていてもよく、スペックル画像(干渉パターン)81pの所定のエリア(例えば、中央の領域)の平均輝度(平均強度)、スペックル画像81pの複数の所定のエリア(例えば、中心の周りに対称な複数の領域、左右、上下)の輝度の差分(強度差分)などのデータに変換または圧縮された情報であってもよい。
図4に、ブロードバンドのストークス光を用いて得られたCARS散乱光31を分光したスペクトル(CARSスペクトル)84の一例を示している。図4(a)に示したCARSスペクトルはグルコース溶液(20%溶液)のCARSスペクトル84aと、水のCARSスペクトル84bとの一例である。グルコース溶液のCARSスペクトル84aは、波長928.7nm付近にトップがあり、波長925.8nm付近にボトムがあるグルコース特有のピークPgを備えている。波長928.7nmのCARS光31aを生成するためのストークス光36の波長は1155nm、波長925.8nmのCARS光31bを生成するためのストークス光36の波長は1159nmである。図4(b)に示すように、これらの波長の単色のストークス光36をポンプ光35とともにターゲットに照射することにより、グルコースのピークPgのトップとを示すCARS光31aと、ボトムを示すCARS光31bを取得でき、それぞれのCARSスペックル81を取得できる。
図1に示す監視システム1の生成装置10は、レーザー源11から供給されるポンプ光35およびストークス光36を対象物、本例では血管6が存在する生体(指)5の皮膚5aに照射し、皮膚5aを介して血管6にポンプ光35およびストークス光36を集光する。その結果として、血管6の血液7により生成されたCARS光(ラマン散乱光)31を生成装置10はイメージセンサー13で受光し、CARSスペックル81を生成するように構成されている。生成装置10は、レーザー源11から供給されるポンプ光35およびストークス光36を対象物に面したレンズ(対物レンズ)12に導く供給光学系(レーザー光供給光路、第5の光路)25と、対象物で発生したCARS光31であって、対物レンズ12を介して入力されたCARS光31を第1の撮像素子アレイ(EMCCD)13に導く第1の光路(第1の干渉パターン取得光学系)21とを含む。
対物レンズ12の一例はシリコン浸対物レンズ(シリコーンオイル浸対物レンズ)などの液浸対物レンズであり、開口数を大きくし、より小さい分解能を得ることができる。対物レンズ12は、空浸レンズ(エアー対物レンズ)であってもよく、より高い強度の信号を得ることができる。
生成装置10は、さらに、対物レンズ12を介して皮膚5aから(皮膚5aの表面またはその近傍で)反射されたポンプ光35の反射光(散乱光、第2の散乱光)32を第2の撮像素子アレイ(CCD)14に導く第2の光路(第2の干渉パターン取得光学系)22と、対物レンズ12を介して対象物が存在する皮膚5aの表面で反射された可視光33をカメラ15に導く第3の光路(画像取得光学系)23とを含む。生成装置10は、分光分析用にCARS光34を分光器(スペクトロメーター)16へ導く第4の光路(スペクトル取得光学系)24を含んでいてもよい。
それぞれの光路21~25は、プリズム、ミラー、フィルターなどの目的に対応した適当な光学素子を含んでいてもよい。例えば、対物レンズ12の上流に波長980nm以下の光を反射するダイクロイックミラー25aを設け、対物レンズ12から入力された波長が短いCARS光31および可視光画像33を第1の光路21および第3の光路23に分離してもよい。また、第1の光路21と第3の光路23との分岐に波長805nm以下の光を反射するダイクロイックミラー21aを設けて可視光画像33を第3の光路23に分離するようにしてもよい。CARSスペックル81を捉える撮像素子アレイ13の上流には、900-960nmの光が通過するようなフィルター群21bを設けてもよい。ポンプ光35およびストークス光36を対物レンズ12に導く光路25には、対物レンズ12から入力されたポンプ光35の散乱光32を第2の光路22に導くためのビームスプリッタ25bを設けてもよい。
生成装置10は、対物レンズ12を操作してポンプ光35およびストークス光36のスポット位置(照射位置)を3次元で調整できるスキャン装置(スキャナー)17を備えていてもよい。スキャナー17は、カメラ15により得られる可視光画像33に基づいて自動的に制御する機能を含んでいてもよい。また、後述する管理装置50の解析結果または測定のための自動的なシーケンス(プロトコル)に基づいて制御されてもよい。スキャナー17は、皮膚5aに沿った2次元方向(X-Y方向)にスキャンできる機能に加えて、皮膚下の血管6の位置にポンプ光35およびストークス光36のスポット位置を調整できるように深度方向(Z方向)にスキャンできる機能を備えていてもよい。
生成装置10においては、CARS光31が導かれる第1の光路21の第1の撮像素子アレイ13にCARS光の干渉パターン(CARSスペックル)81pが生成される。生成装置10により生成された干渉パターン81pを含む第1の干渉情報81は管理装置(マネージメント装置、モニター装置)50に送られる。ポンプ光35の皮膚5aからの散乱光32は、ダイクロイックミラー25aを透過して第2の光路22に導かれ、第2の撮像素子アレイ14にポンプ光の干渉パターン(レーザースペックル、ポンプレーザースペックル)82pが生成される。生成装置10により生成された干渉パターン82pを含む第2の干渉情報82は管理装置50に送られる。第2の撮像素子アレイ14の上流には、ポンプ光の散乱光32以外の波長の光、例えば、ストークス光36を分離するためのフィルター22aなどを設けておいてもよい。生成装置10により生成される可視光カメラ15により得られる画像(画像情報)83および分光器16で得られるスペクトル情報84も管理装置50に送られてもよい。
図5に、生成装置10により得られるいくつかの情報の例を示している。図5(a)は、可視光画像83の一例を示しており、血管6と、血管6の位置に合わせて照射されたポンプ光35およびストークス光36が映し出されている。可視光画像83を参照することにより、ポンプ光35およびストークス光36が血管6またはその付近にスポットを形成するように照射されていることを確認できる。図5(b)は、ポンプ光35の散乱光32の干渉パターン82pを含む第2の干渉情報(以降においては、ポンプレーザースペックル)82の一例を示す。ポンプ光35の散乱光32には、ポンプ光35が照射されている皮膚5aの表面の構造および皮膚5aの動き、すなわち、生体(指)5とポンプ光35との間の相対的な微小な3次元方向の動き(変位)の情報が含まれており、それらの情報はポンプレーザースペックル82に反映される。したがって、ポンプレーザースペックル82の時間的な変位を参照することによりCARS光31の発生源、すなわち、ポンプ光35とストークス光36のスポット(照射位置)と指5との相対的な挙動(移動、変位)を3次元で把握でき、所望の場所にスポットが維持できているかを判定するための1つの要素となり得る。
図5(c)にブロードバンドストークス光により得られたCARS光の干渉パターン85の一例を示している。図5(d)に波長1155nmのナローバンドストークス光36により得られたCARSスペックル81(81a)の干渉パターン81pの一例を示し、図5(e)に波長1159nmのナローバンドストークス光36により得られたCARSスペックル81(81b)の干渉パターン81pの一例を示している。本例の干渉パターン81pは、EMCCD(1024×1024ピクセル)によりフレームレートが4.9m秒(202.55Hz)で取得されたものであるが、フレームレートなどは本例に限定されない。ナローバンドストークス光36により得られるCARS光31のバンド幅は狭いので、そのCARSスペックル81は干渉または回折パターンが鮮明に現れる。このため、ナローバンドストークス光(単色のストークス光)36によるCARS光31のCARSスペックル81を参照することにより、CARS光31の発生源の挙動、本例においては血管6の血液7中の組成成分の挙動をより精度よく把握できる。
図5(d)に示したCARSスペックル81aは、CARSスペクトルのグルコースのピークのトップである928.7nmのCARS光31の干渉パターン81pに対応しており、図5(e)に示したCARSスペックル81bは、CARSスペクトルのグルコースのピークのボトムである925.8nmのCARS光31の干渉パターン81pに対応している。したがって、これらのCARSスペックル81aおよび81bの強度差、図5(f)に示すように、例えば、干渉パターン81pの中心89を含む所定の領域(中央およびその周辺の中央領域)88の輝度の平均値の差分を取得することにより、グルコースの濃度(含有率、含有量)に対応する情報を取得することができる。
また、本例に示す干渉パターン81pについて、本願の発明者らは、非対称な強度分布があることを見出し、その非対称な強度分布が、所定の成分またはその所定の成分に関連する物体(例えば、所定の成分をより多く含む物体、所定の成分をより少なく含む物体、挙動に付随して所定の成分が変動する物体など)の様々な挙動を示唆するものであることを見出した。たとえば、干渉パターン81pの強度が同一であっても、左右の領域87aおよび87bの強度差が反転することがあり、それにより発生源の移動の状況に関わる情報を得ることができる。一例は、血管中を血漿または血球などの組成成分がラマン散乱光を発生させるレーザーのスポット(照射位置、フォーカス)に入ったり、出たりすることに関する情報が含まれていると考えられる。すなわち、血球(赤血球、白血球など)および血漿の血管内の移動であり、その挙動が干渉パターン81pに反映され、干渉パターンの非対称な強度分布の変動を時間に即して分析することにより、対象物内の所定の成分またはそれに関連する物体の挙動を分析できる。
図5(f)に示すように、例えば、干渉パターン81pに含まれる対称な複数の領域(本例では、干渉パターン81pの中心89に対して左右の領域)87aおよび87bの強度差に基づき、対象物内の所定の成分またはそれに関連する物体の挙動を解析するように構成されていてもよい。対称な複数の領域は上下の領域であってもよく、領域の形状は方形であってもよく、その他の多角形、半円、扇状などであってもよい。干渉パターン81pの強度(輝度))を求める領域のサイズ、形状は、アプリケーションに即して決定してもよい。
図1に示した生体監視システム1の管理装置50は、生成装置10のレーザー源11を制御するレーザー制御装置51と、生成装置10からCARSスペックル81などの情報を取得する情報取得装置(入力装置、入力インターフェイス)52と、取得された情報を分析する分析装置53と、取得された情報や分析用の情報、さらに、プログラム54pなどを格納するデータベース(ライブラリ)54と、投薬システム90を制御する投薬制御装置95と、ディスプレイやタッチパネルなどのユーザーインターフェイス61を制御する入出力インターフェイス60を含んでもよい。入出力インターフェイス60は、無線または有線により外部の機器と通信する機能60aと、分析装置53の結果および/またはそれによるアラームをユーザーインターフェイス61や所定の機関(病院、セキュリティーサービスプロバイダーなど)に通報するアラーム機能60bとを備えていてもよい。管理装置50の一例はCPUおよびメモリなどのコンピュータ資源を備えたコンピュータまたはサーバーなどの情報処理装置であり、データベース54などのコンピュータに読み取り可能な適当な記録媒体に記録されたプログラム54pをダウンロードして所定の機能を発揮する。
情報取得装置(入力装置)52は、対象物からラマン散乱光の2次元または3次元の干渉パターン81pの時間的な変動を含む第1の干渉情報(CARSスペックル)81を取得する第1の情報取得装置(第1の情報取得機能、情報取得部)52aと、ラマン散乱光を発生させるために対象物に照射される単色のレーザー光(本例ではポンプ光)35の対象物からの2次元または3次元の干渉パターン82pの時間的な変動を含む第2の干渉情報(ポンプレーザースペックル)82を取得する第2の情報取得装置(第2の情報取得機能、情報取得部)52bと、ラマン散乱光を発生させるためにレーザー光が照射される対象物の表面の画像情報83を取得する第3の情報取得装置(第3の情報取得機能、情報取得部)52cと、ラマン散乱光を分光したラマンスペクトル84を取得する第4の情報取得装置(第4の情報取得機能、情報取得部)52dとを含む。本例の第1の干渉情報81は、上述したように2つの単色のストークス光36によりそれぞれ生成される2つの単色のラマン散乱光31の干渉パターン81aおよび81bを含んでもよい。第1の干渉情報81は、1つの単色のラマン散乱光31の干渉パターン81pのみを含んでいてもよい。例えば、グルコースのピークのトップに相当する波長のラマン散乱光31の干渉パターン81pのみを取得しても、グルコースの濃度を推定することは可能である。しかしながら、グルコースなどの監視する目標となる成分の定量的な情報を精度よく取得するためにはピークのトップとボトムと、など複数の波長のラマン散乱光31の干渉パターン81pを取得してもよい。
本例の第1の干渉情報81は、CARSスペックルであり、ポンプ光35およびストークス光36の少なくとも2つのレーザー光によるラマン散乱光(CARS光)31の干渉パターン81pを含む。CARS光31は、ポンプ光35、ストークス光36に加えて、例えば波長が780nm前後のプローブ光を照射することにより生成されたものであってもよく、ポンプ光35およびストークス光36に対して遅延したプローブ光により生成されたTD-CARS光(時間依存型のCARS光)であってもよい。複数のレーザー光35および36(さらにプローブ光)により散乱光を生成することにより散乱光を誘導生成することが可能となり、共鳴成分が得られるなど、ラマン散乱光を取得することに関するメリットは多い。また、2つ以上のレーザー光を用いることにより、散乱光を生成する位置を前後(深さ)方向も含めた3次元で制御することが可能となる。したがって、皮膚下の血管6中の血液7や他の組織を対象物(ターゲット)として、ターゲットの情報がより強く反映された散乱光を取得できる。
本例ではCARSスペックル81を取得しているが、逆ラマン、コヒーレント・反ストークスラマン散乱(CARS)に限定されず、誘導ラマン散乱(SRS)、ラマン増幅、コヒーレント・ストークスラマン散乱(CSRS)、ハイパーラマン散乱などにより得られるラマン散乱光の干渉パターンを取得するものであってもよい。また、本例においては、干渉パターンをCCDなどの2次元の撮像素子アレイにより生成しているが、散乱光は完全にコリメートされた光ではないので、3次元の干渉パターンを得てもよく、例えば、散乱光をビームスプリッタなどで分割して異なる深さにおける複数の干渉パターンを取得するようにしてもよい。散乱光に含まれる干渉および/または回折に関わる、より多くの情報をさらに精度よく取得できる。
分析装置53は、第1の干渉情報であるCARSスペックル81に基づき、血液などの対象物内の所定の成分またはそれを含む物体の挙動を分析する第1の分析装置(CARSスペックル分析装置)56と、第2の干渉情報であるポンプレーザースペックル82に基づき、対象物の挙動(相対的な挙動)を分析する第2の分析装置(ポンプレーザースペックル分析装置)57と、画像情報83を分析する第3の分析装置(画像分析装置)58と、CARSスペクトル84を分析する第4の分析装置(スペクトル分析装置)59と、分析装置53によるこれらの分析を自動的に行うための自動分析装置55とを含む。
第1の分析装置56は、第1の干渉情報(CARSスペックル)81に含まれる干渉パターン81pの非対称な強度分布に基づき、対象物内のラマン散乱光の発生要因となる所定の成分またはそれを含むそれに関連する物体の時間的な挙動を分析するように構成された第1のパターン解析装置(第1の挙動解析装置、第1のパターン解析機能、解析ユニット)56aを含む。第1のパターン解析装置56aは、干渉パターン81pに含まれる対称な複数の領域(例えば、左右の領域87aおよび87b)の強度差に基づき、対象物内の所定の成分またはそれに関連する物体の挙動を解析するように構成されていてもよい。
第1の分析装置56は、さらに、干渉パターン81pの中央およびその周辺の中央領域88の強度に基づき、挙動に関連付けられた対象物内の所定の成分またはそれに関連する物体の濃度を求めるように構成された第2のパターン解析装置56b(第2の挙動解析装置、解析機能、解析ユニット)を含んでもよい。第2のパターン解析装置56bは、第1の干渉情報81に含まれる複数の単色のラマン散乱光の干渉パターン81aおよび81bの中央領域88の強度差分に基づき、対象物内の所定の成分またはそれに関連する物体の濃度を求めるように構成されていてもよい。第1の分析装置56は、さらに、対象物内の所定の成分またはそれに関連する物体の挙動に関連する物性値を得るように構成された第3のパターン解析装置(第3の挙動解析装置)56cを含んでいてもよい。
図6および図7に、ポンプ光35と波長1155nmのナローバンドストークス光36とにより生体中から得られた2次元のCARSスペックル81の一例を示している。CARSスペックル81の一例は、干渉パターン81pをそれぞれ含むフレームを、そのフレームが取得された時間に沿って束ねたものである。フレームレートは数10μ秒から数100m秒の範囲であってもよく、数100μ秒から数10m秒の範囲であってもよい。フレームレートは挙動を得ようとしている対象物の移動(変位)速度よりも十分に短ければよく、測定できる程度の強度の干渉パターン81pが得られる程度であればよい。
図6は、多数のフレームを含むCARSスペックル81の中心を通る横方向(X軸、中心線)上の強度(輝度)の時間変化(縦方向)を示している。図7にCARSスペックル81の一部を拡大して示している。これらのCARSスペックル81は、生成する際のストークス光36の波長より、発生源(特徴として現れる成分)は生体中のグルコースであり、生体中の成分としてのグルコースおよびそれに関連する物体、例えば、組織、細胞、血液中の血球および/または血漿、リンパ液中のリンパ球および/または血漿などの時間的な挙動や濃度などを反映していると想定される。
図6(a)および図7(a)は、皮下組織から得られたCARSスペックル81を示しており、時間的な強度変化はほとんど見られない。図6(b)および図7(b)は、血管6内から得られたCARSスペックル81を示しており、一定の間隔で強弱(強度が変化、輝度が変化)を繰り返していることが分かる。したがって、第1の分析装置56により、CARSスペックル81の中心を含む所定の領域のフレーム単位(時間経過)の強度差を得ることにより、本例の対象物である血管6内の血液7の挙動(流れ)に伴う物性値、例えば、脈拍、血流速度などを把握できる。また、CARSスペックル81の強度が時間経過により変動していることにより、血液7がCARS光31の発生源であることを推測できる。したがって、CARSスペックルの分析装置である第1の分析装置56は、CARSスペックル81の全体的な解析により、血管6からの信号であるか否かを判断する第4のパターン解析装置56dを含んでいてもよい。
図8に、血液中のグルコース由来のCARSスペックル81のX軸上の強度を示している。図8(a)は、CARSスペックル81の経過(動的スペックル画像、スペックルビデオ、スペックル流)に含まれる強度の高いフレームのX軸上の強度変化、すなわち、干渉パターン81pのX軸上の強度変化を抜き出して示している。図8(b)は、CARSスペックル81の経過に含まれる強度の低いフレームのX軸上の強度変化、すなわち、干渉パターン81pのX軸上の強度変化を抜き出して示している。これらを観察すると、左右の強度に規則的な非対称性があることが認められる。
図9に、血液中のグルコース由来のCARSスペックル81を時間方向に解析した結果の一例を示している。図9は、縦軸(Y軸)に、CARSスペックル81の中心から所定の領域を複数、本例では、図5(f)に示す左右の領域87aおよび87bに分けて強度差(輝度の差分)を算出した結果を示し、横軸(X軸)に、CARSスペックル81のほぼ全体に対応する所定の領域の平均強度、例えば中央領域88の強度(輝度)を示し、Z軸に時間経過(フレーム数)を示し、フレームごとに算出した結果をプロットした様子を示している。本図(3Dプロファイル)に示されているように、CARSスペックル81に含まれる情報は、上記のように処理することにより、時間的に螺旋を描くように変動していることが分かる。これは、そのポンプ光35およびストークス光36のスポット(フォーカス)を形成している血管6中を血液7が流れ、そのスポットに赤血球が出入りしていることを示していると想定される。
このため、第1のパターン解析装置56aにより、CARSスペックル81の複数の領域の強度差分に基づき血液7の中の赤血球の挙動を解析できる。この結果、例えば、CARSスペックル81の多数のフレームの中から赤血球ではなく血漿が主に寄与しているフレームを抽出することができる。したがって、第2のパターン解析装置56bにより、血液由来のCARSスペックル81の中から血糖値の測定に適した血漿部分のCARSスペックル81を抽出し、その中央領域88の強度を求めることにより、CARSスペックル81に主に反映されているグルコースの濃度を求めることができる。また、第3のパターン解析装置56cにより、CARSスペックル81の多数のフレームの中から赤血球が主に寄与しているフレームと、血漿が主に寄与しているフレームとを抽出することが可能となり、赤血球の挙動に関連する物性値として、例えば、赤血球と血漿との比であるヘマトクリット値を求めることが可能となる。
図10に、波長1155nmのナローバンドストークス光36によるCARSスペックル81aと、波長1159nmのナローバンドストークス光36によるCARSスペックル81bとを順番に、例えば3秒単位(3秒間隔、3秒の時間分割)で取得した例を示している。図10(a)に示すCARSスペックルの左右の領域の強度の差分から、それぞれのCARSスペックル81aおよび81bの血漿部分を特定することが可能である。したがって、図10(b)に示すCARSスペックル81aおよび81bの平均強度から、図10(a)の情報に基づき血漿部分に相当するフレームを抽出して、それらの強度差を取得することにより、血漿血糖値の値を導くことができる。
したがって、第2のパターン解析装置56bにより、CARSスペックル81に含まれる複数の単色のラマン散乱光の干渉パターン81p、本例においては、上記のそれぞれの波長のナローバンドストークス光36による、CARSスペクトルのグルコースのピークのトップである928.7nmと、ピークのボトムである925.8nmの2色のCARSスペックル81aおよび81bの干渉パターン81pが反映されているフレームを抽出し、それぞれの所定の領域、例えば、中心を含む領域88の強度差分を取得できる。この強度差分は、血液7の血漿部分のグルコース濃度に対応するため、血糖値の定量的で精度の高い情報を、時間変化を含めて取得することが可能となる。血糖値の算出には、図9に示した3Dプロファイルの血漿部分に当たるフレームのCARSスペックルを選択してもよい。
図11に、上記の過程により本例のシステム1を用いて得られた血糖値の推定値と、グルコースメータ(Glucometer(SMBG)ニプロ社製Freestyle Freedom Lite)により得られた実測値(リファレンス)とをコンセンサスエラーグリッドを用いて対比して示している。本システム1による推定値は、実測値(リファレンス)に対しゾーンA(95%以内)に収まり、ISO15197によるグルコース測定器認証の条件の1つをほぼ満たしていることが分かる。したがって、本例のシステム1により、血糖値を常時、リアルタイムで測定するセンサーとしての機能、血糖値モニターとしての機能、さらに、血糖値に基づき人体を管理する機能とを提供できる。さらに、システム1により、血糖値に限らず、血液中の他の物性値(パラメータ)を常時、リアルタイムでモニターする機能もあわせて提供できる。
図12に、干渉パターンを生成するための撮像素子アレイ(画像生成デバイス、イメージセンサー)13の異なる例を示している。図3に示した生成装置10の例では、波長の異なる複数のCARS光31により生成される複数の干渉パターンを1つの撮像素子アレイ13により時分割で取得している。図12に示した撮像素子アレイ13は、図12(a)に示すようなベイヤーパターン(Bayer Pattern、ベイヤー配列、ベイヤーフィルター)を備えており、図12(b)に示すように、波長の異なる3種類の光に敏感な3種類の画素またはフィルター13a~13cが市松模様状に配列されている。この撮像素子アレイ13においては、複数、例えば3種類の単色の干渉パターンを並列に取得できる。しかしながら、撮像素子アレイ13で得られたオリジナルの干渉パターンはモザイク画像となる。このため、各色の画像ごとに、画素間のデータを補完するなどの方法によりデモザイク処理を行うことにより、図12(c)に示すように、各色の画像(干渉パターン)を得ることができる。
図13に、レーザー源11からポンプ光35およびストークス光36をターゲットに照射してCARS光31を生成する生成装置10の他の一例を示している。この生成装置10では、単色のポンプ光35と、ブロードバンドのレーザー光36bから選択された4つの異なる波長の単色(ナローバンド)のストークス光36とが、ターゲットである血管6中の血液7に照射され、血液7を発生源とする4つの異なる波長のCARS光31が生成される。これらの波長の異なるCARS光31はビームスプリッタ18により空間的に分離され、1つまたは複数の撮像素子アレイ(イメージセンサー)13に導かれ、それぞれの波長のCARS光31による干渉パターンが並列に(同時に)生成される。ビームスプリッタ18は、波長の異なるCARS光31を1つのイメージセンサー13の異なる領域に干渉パターンが生成されるように調整されていてもよく、異なるイメージセンサー13に異なる波長のCARS光31の干渉パターンがそれぞれ生成されるように調整されていてもよい。ナローバンドのストークス光36は、ブロードバンドのストークス光36bからバンドパスフィルターや回折格子(グレーチング)などを用いて生成されてもよい。ナローバンドの6はチューナブルレーザーを用いて生成されてもよい。
図14および図15に、人体の手首(さらに具体的には手首の内側(手のひら側))5の血管を流れる血液中の血糖値を測定する生体監視システム(モニタリングシステム、モニター)1の一例を示している。図14(a)に示すように、このモニタリングシステム1は、縦横が250mm、高さが120mm程度の可搬式(ポータブル)なデスクトップタイプの装置であり、生成装置10および管理装置50を内蔵したハウジング66と、ユーザーインターフェイス61を兼ねたカバー65とを含む。
図14(b)に示すように、カバー65をオープンするとユーザーインターフェイスとしてタッチパネル機能が付いたディスプレイ61が現れ、ハウジング66の上面にユーザーの手首5をセットする機構64が現れる。ハウジング66の上面には手(アーム)を載せて手首を所定の位置に合わせやすいように手の形状に凹んだアーム搭載部67と、手首を生成装置10の対物レンズ12の前面に保持するためのパームレスト68と、手首を保持するためのリストストラップ(リストバンド)69とを含む。
図15に示すように、ユーザーは、この生体管理システム1のカバー65を開けて、ハウジング66の上面に、手首5の内側を生成装置10の対物レンズ12に対峙するようにセットすることで、非侵襲で、リアルタイムに血糖値を測定でき、ディスプレイ61に表示できる。さらに、この生体管理システム1は、ハウジング66に、管理装置50に加えて投薬装置90を含でいてもよく、測定結果に応じた処理が自動的に受けられるように構成されていてもよい。
図16に管理システム1の自動分析装置55により自動的に測定を行う制御プロセス100の概要を示している。ステップ101において、人体の所定の場所に管理システム1を設置して測定を開始する。図15に示す管理システム1であれば、手首5の内側を生成装置10の対物レンズ12に面するようにセットした後、管理システム1のスタートボタンをクリックする。ステップ102において、測定対象となる血管6を自動的に選択する。この処理は、生成装置10により得られる画像情報83を用いて生成装置10のスキャン装置17が自動的(自律的)に行ってもよく、管理装置50の画像解析装置58が行ってもよい。画像情報83に基づき選択される血管6の一例は、血管径が6μm以下で、血流速度が300μm/秒以下の毛細血管で、手首5の平面状の部分から観察できるものである。
次に、ステップ103において、選択された血管6に対して照射するレーザー(レーザースポット、照射位置)の詳細な位置合わせ(アライメント)を行う。この処理も、画像情報83を用いて生成装置10のスキャン装置17が自動的(自律的)に行ってもよく、管理装置50の画像解析装置58が行ってもよい。実際にポンプ光35とストークス光36とを手首5の皮膚5aに照射した画像情報83を取得して、レーザースポットと血管6との位置合わせを行ってもよい。
ステップ104において、生成装置10でポンプ光35とストークス光36とを手首5の皮膚5aに照射し、CARS光31を生成し、干渉パターン81pを含む干渉情報(CARSスペックル)81を取得する。生成装置10は、CARSスペックル81とともに、ポンプレーザースペックル82および画像情報83を生成して管理装置50に供給する。ステップ105において、管理装置50は、取得装置52を介して取得したこれらの情報(データ)を分析装置53により分析して情報の精度を検証する。検証するアイテムとしては、(1)血液由来の信号(データ)であるか、(2)血液以外の信号が混入していないか、(3)位置安定性があるか、などを挙げることができる。
血液由来の信号であるか否かの判断は、例えば、第1の分析装置56の第4のパターン解析装置56dにより、CARSスペックル81に含まれる干渉パターン81pの中央領域88の強度(平均強度)が組織を観察したときと比較し、10-100m秒の時間スケールの輝度強弱の振動的な増減の有無で判断してもよい。血液以外の信号が混入していなかの判断については、例えば、第1の分析装置56の第1のパターン解析装置56aにより、図9に示した3Dプロット画像において、所定の時間(所定のフレーム数)が経過したときに、安定した、らせん構造が観察されるか否かで判断してもよい。例えば、1秒以上で、10回/秒程度の安定した、らせん構造が観察されるか否かを基準としてもよい。
位置安定性の判断については、情報取得装置52より取得された、ラマン散乱光を発生させるために対象物に照射される単色のレーザー光(本例においてはポンプ光35)の対象物(皮膚5a)からの干渉パターン82pの時間的な変動を含む第2の干渉情報(ポンプレーザースペックル)82を、第2の分析装置57で分析して得られる対象物(皮膚5aまたは人体(手首))5の挙動で判断してもよい。例えば、干渉パターン82の中央領域の輝度(強度)に10%以上の変動が観察される場合は、対象物の位置が安定していないと判断してもよい。画像情報83を画像分析装置58により分析して、ポンプ光35および/またはストークス光36のスポットと血管6との相対的な位置に振動またはずれが観察されることから、対象物の位置が安定していないと判断してもよい。
第1の分析装置53の第4の分析装置(スペクトル分析装置)59が、断続的に、または定期的にCARSスペクトル84を解析して、所定のCARS信号、例えば、グルコースに対応するピークまたはボトムが得られているか否かを確認する処理をおこなってもよい。
情報の精度が検証されると、ステップ106において、分析装置53がデータ処理を行う。情報取得装置52が取得したラマン散乱光の干渉パターン81pの時間的な変動を含む第1の干渉情報(CARSスペックル)81を第1の分析装置56が解析する。具体的には、ステップ106aにおいて、第1のパターン解析装置56aが、干渉パターン81pの非対称な強度分布に基づき対象物である血管6内のラマン散乱光の発生要因となる所定の成分(本例ではグルコース)に関連する物体(本例では血球(赤血球)および血漿)の時間的な挙動を分析する。さらに具体的には、第1のパターン解析装置56aが、干渉パターン81pに含まれる対称な複数の領域(左右の領域)87aおよび87bの強度差に基づき、血管中の血球および血漿の挙動を解析する。
ステップ106bにおいて、第2のパターン解析装置56bが、干渉パターン81pの中央およびその周辺の領域88の強度に基づき、挙動に関連付けられた対象物内の所定の成分またはそれに関連する物体の濃度を求める。本例では、第2のパターン解析装置56bは、対象物である血液7の中で挙動が確認された物体の1つである、血漿中のグルコースの値(血漿血糖値)を求める。さらに具体的には、第1のパターン解析装置56aにより生成される図9に示した3Dプロット画像のらせん状のパターンのボトム(左右の領域87aおよび87bの強度差が小さい部分)に対応するフレームをCARSスペックル81から選択して、その中央領域88の強度に基づき血漿血糖値を求める。第2のパターン解析装置56bは、CARSスペックル81に含まれる複数の単色のラマン散乱光の干渉パターン(本例では、グルコースの1つのピークの波長の干渉パターンおよび、ボトムの波長の干渉パターン)81pに含まれる中央領域88の強度差分に基づき血漿血糖値の値を求めることができる。第2のパターン解析装置56bは、血漿血糖値に限らず、血球(赤血球)の血糖値を得ることも可能であり、生成装置(検出装置)10から得られたCARSスペックル81の中で、血糖値を算出するために最も適した複数のフレームを選択し、それらに基づき所望の血糖値を取得してもよい。
さらに、ステップ106cにおいて、第1の分析装置56の第3のパターン解析装置56cは、対象物内の所定の成分またはそれに関連する物体の挙動に関連する物性値、例えば、赤血球と血漿との比率を示すヘマトクリット値を得てもよい。第1の分析装置56は、複数の単色のCARS散乱光の干渉パターン81pを用いて、グルコース以外の、CARS散乱光の発生要因となる成分、すなわち、CARS散乱に大きな影響を与える成分またはそれに関連する物質の情報を取得してもよい。例えば、血糖値に関連して体内の異常を検出できる対象(成分)としては、グルコースの他に、ヘモグロビンA0、ヘモグロビンA1c、グリコアルブミン、アルブミン、アンヒドログルシトール、フルクトサミン、インシュリン、およびグルカゴンを挙げることができる。腎機能の異常の検出できる対象(成分)としては、クレアチニンおよびアルブミンなどを挙げることができる。この管理システム1を用いて、その他の血中成分を測定してモニターしてもよい。また、管理システム1は、血管または血液に関し、血流、脈動(拍動)などの他の情報を取得してもよく、これらを含めてモニタリングの対象としてもよい。
ステップ107において、入出力インターフェイス60は、測定された様々な結果をユーザーインターフェイス61に表示する。入出力インターフェイス60は、外部の所定のシステムに対して測定された結果を提供してもよい。ステップ108において、投薬制御装置95は測定された結果より投薬の要否を判断する。投薬制御装置95が、投薬、例えばインスリンの注入が必要であると判断すると、ステップ109において投薬システム90にユーザーに対して所定の薬剤を注入する処置を行う。
このように、本例のシステム1は、対象物として、生体5の血管中6の血液7に少なくとも2つのレーザー光35および36を照射し、そのラマン散乱光の干渉パターンを含む第1の干渉情報(CARSスペックル)81を生成する生成装置(検出装置)10を有する。また、第1の分析装置56は、血液中の、血漿成分および/または血球成分の時間的な挙動を分析でき、上記においては血漿成分に含まれるグルコース濃度から血糖値を算出している。同様に、第1の分析装置56は、グルコース以外の成分、たとえば、ヘモグロビンA1c、クレアチニンおよびアルブミンの少なくともいずれかを含む成分の定量的な挙動を分析することができる。
また、上記においては、体液の典型的な例として血管を流れる血液を例に説明しているが、体液はリンパ管を流れるリンパ液などの他の体液に含まれる成分であっても同様に測定できる。主たる構成成分は、血漿成分、赤血球に限らず、他の血球成分、例えば、白血球および/または血小板を含んでいてもよい。また、濃度を測定するターゲットとなる成分は血液などの体液の検査対象となるいずれかの成分を含んでいてもよい。測定対象は、任意の化学物質、分子、化合物、組成物、微生物または凝集物であってよく、例えば、血球、アミノ酸、ペプチド、ポリペプチド、タンパク質、糖タンパク質、リポタンパク質、ヌクレオシド、ヌクレオチド、オリゴヌクレオチド、核酸、糖、炭水化物、オリゴ糖、多糖、脂肪酸、脂質、ホルモン、代謝産物、サイトカイン、ケモカイン、受容体、神経伝達物質、抗原、アレルゲン、抗体、基質、補助因子、阻害物質、薬剤、製剤、栄養分、プリオン、トキシン、毒物、爆発物、農薬、化学兵器物質、生物学的有害物質、放射線同位体、ビタミン、複素環式芳香族化合物、発癌物質、変異誘発物質、麻薬、アンフェタミン、バルビツール酸塩、幻覚発現物質、廃棄物、および/または汚染物質が挙げられるが、これらに限定されるものではない。微生物としては、ウィルス、細菌、細胞などが挙げられるが、これらに限定されるものではない。
ラマン散乱光は、CARS光に限らず、誘導ラマン散乱(SAS、Stimulated Raman Scattering)、表面増強ラマン散乱(SERS、Surface Enhanced Raman Scattering)などの他の公知の方法を用いて取得された散乱光であってもよく、それらの干渉パターンを含む情報を取得することが可能である。
本例のシステム1は、対象物からラマン散乱光の2次元または3次元の干渉パターンの時間的な変動を含む第1の干渉情報を取得する第1の情報取得装置と、第1の干渉情報に基づき、対象物内の所定の成分またはそれを含む物体の挙動を分析する第1の分析装置とを有する。第1の分析装置は、第1の干渉情報に含まれる干渉パターンに含まれる2次元または3次元の所定の領域の強度差に基づき、対象物内の所定の成分またはそれを含む物体の時間的な挙動を解析する挙動解析装置を含んでもよい。
上述したように、本願の発明者らは、干渉パターンに含まれる複数の領域の強度差分を得ることにより、発生源の時間的な変動に関するより詳しい情報が得られることを見出した。したがって、第1の分析装置は、第1の干渉情報に含まれる干渉パターンに含まれる2次元または3次元の複数の領域の強度差分に基づき、対象物内の所定の成分またはそれを含む物体の時間的な挙動を解析する挙動解析装置を含んでもよい。例えば、干渉パターンの強度が同一であっても、左右の領域の強度差が反転することがあり、それにより発生源の移動の状況に関わる情報を得ることができる。一例は、血管中を血漿または血球などの組成成分がラマン散乱光を発生させるレーザーのスポットに入ったり、出たりすることに関する情報が得られることを見出したことである。
第1の分析装置は、第1の干渉情報に含まれる複数の単色のラマン散乱光の干渉パターンに含まれる、所定の領域の強度差分に基づき、対象物内の所定の成分またはそれを含む物体の定量的な挙動を解析する挙動解析装置を含んでもよい。所定の成分の濃度およびその変動をより精度よく定量的に得ることができる。一例は、ラマンスペクトルに含まれる所定の成分のピークのトップとボトムとの波長のラマン散乱光の干渉パターンを取得することである。それらの干渉パターンを解析することにより、所定の成分またはそれが含まれる物体の時間的な挙動と、定量的な挙動とを得ることができ、発生源を特定し、その発生源における所定の成分の濃度を精度よく取得できる。
また、上記においては、本発明の特定の実施形態を説明したが、様々な他の実施形態および変形例は本発明の範囲および精神から逸脱することなく当業者が想到し得ることであり、そのような他の実施形態および変形は以下の請求の範囲の対象となり、本発明は以下の請求の範囲により規定されるものである。
Claims (24)
- 対象物からラマン散乱光の干渉パターンの時間的な変動を含む第1の干渉情報を取得するように構成された第1の情報取得装置と、
前記第1の干渉情報を分析する第1の分析装置とを有し、
前記第1の分析装置は、前記干渉パターンの非対称な強度分布に基づき、前記対象物内の前記ラマン散乱光の発生要因となる所定の成分またはそれに関連する物体の時間的な挙動を分析するように構成された第1のパターン解析装置を含む、システム。 - 請求項1において、
前記第1の情報取得装置は、少なくとも1つの単色の前記ラマン散乱光の前記干渉パターンを含む前記第1の干渉情報を取得するように構成されている、システム。 - 請求項1または2において、
前記第1のパターン解析装置は、前記干渉パターンに含まれる対称な複数の領域の強度差に基づき、前記対象物内の前記所定の成分またはそれに関連する物体の前記挙動を解析するように構成されている、システム。 - 請求項1ないし3のいずれかにおいて、
前記第1の分析装置は、前記干渉パターンの中央およびその周辺の中央領域の強度に基づき、前記挙動に関連付けられた前記対象物内の前記所定の成分またはそれに関連する物体の濃度を求めるように構成された第2のパターン解析装置を含む、システム。 - 請求項4において、
前記第2のパターン解析装置は、前記第1の干渉情報に含まれる複数の単色の前記ラマン散乱光の前記干渉パターンの前記中央領域の強度差分に基づき、前記対象物内の所定の成分またはそれに関連する物体の濃度を求めるように構成されている、システム。 - 請求項1ないし5のいずれかにおいて、
前記第1の分析装置は、前記対象物内の前記所定の成分またはそれに関連する物体の前記挙動に関連する物性値を得るように構成された第3のパターン解析装置を含む、システム。 - 請求項1ないし6のいずれかにおいて、
前記第1の情報取得装置は、少なくとも2つのレーザー光によるラマン散乱光の前記干渉パターンを含む前記第1の干渉情報を取得するように構成されている、システム。 - 請求項1ないし7のいずれかにおいて、
前記第1の情報取得装置は、誘導ラマン散乱(SRS)、ラマン増幅、逆ラマン、コヒーレント・反ストークスラマン散乱(CARS)、コヒーレント・ストークスラマン散乱(CSRS)、ハイパーラマン散乱のいずれかの散乱光の前記干渉パターンを含む前記第1の干渉情報を取得するように構成されている、システム。 - 請求項1ないし8のいずれかにおいて、
前記ラマン散乱光を発生させるために前記対象物に照射される単色のレーザー光の前記対象物からの前記干渉パターンの時間的な変動を含む第2の干渉情報を取得するように構成された第2の情報取得装置と、
前記第2の干渉情報に基づき、前記対象物の挙動を分析するように構成された第2の分析装置とを有する、システム。 - 請求項1ないし9のいずれかにおいて、
前記ラマン散乱光を発生させるためにレーザー光が照射される前記対象物の表面の画像を取得するように構成された第3の情報取得装置と、
前記画像から前記レーザー光の照射位置を解析するように構成された画像解析装置とを有する、システム。 - 請求項1ないし10のいずれかにおいて、
前記対象物に前記ラマン散乱光を発生させるためのレーザー光を照射して前記対象物からのラマン散乱光を受光することにより前記干渉パターンを含む前記第1の干渉情報を生成するように構成された生成装置を有する、システム。 - 請求項1ないし10のいずれかにおいて、
前記対象物として、生体の血管中の血液に少なくとも2つのレーザー光を照射し、そのラマン散乱光の前記干渉パターンを含む前記第1の干渉情報を生成するように構成された生成装置を有し、
前記所定の成分は、グルコース、ヘモグロビンA0、ヘモグロビンA1c、グリコアルブミン、アルブミン、アンヒドログルシトール、フルクトサミン、インシュリン、グルカゴン、クレアチニンおよびアルブミンの少なくともいずれかを含む、システム。 - 請求項12において、
前記第1の分析装置は、前記血液中の、血漿成分および/または血球成分の前記挙動を分析する、システム。 - 請求項11ないし13のいずれかにおいて、
前記第1の分析装置により得られる前記生体の状態に基づいて、前記生体に薬剤を注入する注入装置を有する、システム。 - モニタリングシステムの制御方法であって、
前記モニタリングシステムは、対象物からラマン散乱光の情報を取得する情報取得装置と、前記第1の干渉情報を分析する分析装置とを有し、
当該方法は、前記情報取得装置が、前記ラマン散乱光の干渉パターンの時間的な変動を含む第1の干渉情報を取得することと、
前記分析装置が、前記干渉パターンの非対称な強度分布に基づき、前記対象物内の前記ラマン散乱光の発生要因となる所定の成分またはそれに関連する物体の時間的な挙動を分析することを有する、方法。 - 請求項15おいて、
前記第1の干渉情報を取得することは、少なくとも1つの単色の前記ラマン散乱光の前記干渉パターンを含む前記第1の干渉情報を取得することを含む、方法。 - 請求項15または16において、
前記第1の干渉情報を取得することは、少なくとも2つのレーザー光によるラマン散乱光の前記干渉パターンを含む前記第1の干渉情報を取得することを含む、方法。 - 請求項15ないし17のいずれかにおいて、
前記分析することは、前記干渉パターンに含まれる対称な複数の領域の強度差に基づき、前記対象物内の前記所定の成分またはそれに関連する物体の前記挙動を解析することを含む、方法。 - 請求項15ないし18のいずれかにおいて、
前記分析することは、前記干渉パターンの中央およびその周辺の中央領域の強度に基づき、前記挙動に関連付けられた前記対象物内の前記所定の成分またはそれに関連する物体の濃度を求めることを含む、方法。 - 請求項19において、
前記濃度を求めることは、前記第1の干渉情報に含まれる複数の単色の前記ラマン散乱光の前記干渉パターンに含まれる、前記中央領域の強度差分に基づき、前記対象物内の所定の成分またはそれに関連する物体の濃度を求めることを含む、方法。 - 請求項15ないし20のいずれかにおいて、
前記情報取得装置が、前記ラマン散乱光を発生させるために前記対象物に照射される単色のレーザー光の前記対象物からの前記干渉パターンの時間的な変動を含む第2の干渉情報を取得することと、
前記分析装置が、前記第2の干渉情報に基づき、前記対象物の挙動を分析することとを有する、方法。 - 請求項15ないし21のいずれかにおいて、
前記第1の干渉情報を取得することは、前記対象物として、生体の血管中の血液に少なくとも2つのレーザー光を照射されたラマン散乱光の前記干渉パターンを含む前記第1の干渉情報を取得することを含む、方法。 - 請求項22において、
投薬装置により、前記分析することにより得られる前記生体の状態に基づいて、前記生体に薬剤を注入することを有する、方法。 - 請求項1ないし14のいずれかに記載のシステムとしてコンピュータを稼働させるための命令を有するプログラム。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2023-102381 | 2023-06-22 | ||
| JP2023102381 | 2023-06-22 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2024262621A1 true WO2024262621A1 (ja) | 2024-12-26 |
Family
ID=93935864
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/JP2024/022648 Pending WO2024262621A1 (ja) | 2023-06-22 | 2024-06-21 | ラマン散乱光を用いたシステムおよび方法 |
Country Status (1)
| Country | Link |
|---|---|
| WO (1) | WO2024262621A1 (ja) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070121119A1 (en) * | 2005-10-21 | 2007-05-31 | Southwest Research Institute | Spatial Heterodyne Wide-Field Coherent Anti-Stokes Raman Spectromicroscopy |
| JP2009204613A (ja) * | 2008-02-27 | 2009-09-10 | Gwangju Inst Of Science & Technology | コヒーレントアンチストークスラマン散乱を用いた映像獲得装置および方法 |
| JP2018536856A (ja) * | 2015-11-17 | 2018-12-13 | 韓國科學技術院Korea Advanced Institute Of Science And Technology | 混沌波センサを利用した試料特性探知装置 |
| US20190385281A1 (en) * | 2017-02-14 | 2019-12-19 | ContinUse Biometrics Ltd. | System and method for monitoring a sample with at least two wavelengths |
| WO2022084539A2 (en) * | 2020-10-23 | 2022-04-28 | Rsp Systems A/S | A raman probe and apparatus and method for non-invasive in vivo measurement of analyte presence or concentration |
| KR20220067222A (ko) * | 2020-11-17 | 2022-05-24 | 한국전기연구원 | 광열 모듈레이션 레이저 스펙클 영상을 이용한 포도당 농도 측정 시스템, 방법, 및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한 프로그램을 기록한 기록 매체 |
-
2024
- 2024-06-21 WO PCT/JP2024/022648 patent/WO2024262621A1/ja active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070121119A1 (en) * | 2005-10-21 | 2007-05-31 | Southwest Research Institute | Spatial Heterodyne Wide-Field Coherent Anti-Stokes Raman Spectromicroscopy |
| JP2009204613A (ja) * | 2008-02-27 | 2009-09-10 | Gwangju Inst Of Science & Technology | コヒーレントアンチストークスラマン散乱を用いた映像獲得装置および方法 |
| JP2018536856A (ja) * | 2015-11-17 | 2018-12-13 | 韓國科學技術院Korea Advanced Institute Of Science And Technology | 混沌波センサを利用した試料特性探知装置 |
| US20190385281A1 (en) * | 2017-02-14 | 2019-12-19 | ContinUse Biometrics Ltd. | System and method for monitoring a sample with at least two wavelengths |
| WO2022084539A2 (en) * | 2020-10-23 | 2022-04-28 | Rsp Systems A/S | A raman probe and apparatus and method for non-invasive in vivo measurement of analyte presence or concentration |
| KR20220067222A (ko) * | 2020-11-17 | 2022-05-24 | 한국전기연구원 | 광열 모듈레이션 레이저 스펙클 영상을 이용한 포도당 농도 측정 시스템, 방법, 및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한 프로그램을 기록한 기록 매체 |
Non-Patent Citations (2)
| Title |
|---|
| FANTUZZI ERIC M., HEUKE SANDRO, LABOUESSE SIMON, GUDAVIČIUS DOMINYKAS, BARTELS RANDY, SENTENAC ANNE, RIGNEAULT HERVÉ: "Random illumination coherent anti-Stokes Raman scattering microscopy (RIM-CARS)", OPTICA IMAGING CONGRESS 2024 (3D, AOMS, COSI, ISA, PCAOP), OPTICA PUBLISHING GROUP, 15 March 2023 (2023-03-15), pages JF2A.15, XP093251176, ISBN: 978-1-957171-37-1, DOI: 10.1364/3D.2024.JF2A.15 * |
| FANTUZZI ERIC M., HEUKE SANDRO, LABOUESSE SIMON, GUDAVIČIUS DOMINYKAS, BARTELS RANDY, SENTENAC ANNE, RIGNEAULT HERVÉ: "Wide-field coherent anti-Stokes Raman scattering microscopy using random illuminations", NATURE PHOTONICS, NATURE PUBLISHING GROUP UK, LONDON, vol. 17, no. 12, 1 December 2023 (2023-12-01), London, pages 1097 - 1104, XP093251194, ISSN: 1749-4885, DOI: 10.1038/s41566-023-01294-x * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109124649B (zh) | 用于通过拉曼光谱进行无创伤性活体测量的设备 | |
| Wang et al. | Single-cell label-free photoacoustic flowoxigraphy in vivo | |
| CN106604677B (zh) | 用于通过拉曼光谱进行透皮体内测量的方法和装置 | |
| JP5047788B2 (ja) | 分光システムのためのオートフォーカス機構 | |
| US8054463B2 (en) | Method and system for measuring sub-surface composition of a sample | |
| CN109171762A (zh) | 对生物体进行监视的监视器和系统 | |
| JP2009536317A (ja) | ラマン分析 | |
| US11125692B2 (en) | Determination method, determination apparatus, and recording medium | |
| CN1981189B (zh) | 具有多探针的光谱系统 | |
| Burmeister et al. | Infrared Spectroscopy Spectroscopic Considerations for Noninvasive Blood Glucose Measurements with Near Infrared Spectroscopy | |
| WO2024262621A1 (ja) | ラマン散乱光を用いたシステムおよび方法 | |
| Pain et al. | Visualizing odor representation in the brain: a review of imaging techniques for the mapping of sensory activity in the olfactory glomeruli | |
| US7544503B2 (en) | Apparatus for the pH determination of blood and method therefor | |
| CN108603841A (zh) | 分析方法和分析装置 | |
| JP7500117B2 (ja) | 体液に含まれる成分の濃度を測定するシステムおよび方法 | |
| Karakostas et al. | Portable optical blood scattering sensor | |
| JP7573912B2 (ja) | インプラント、生体監視システムおよび生体管理システム | |
| US20230285136A1 (en) | Implant, monitoring system for organism, and management system for organism | |
| GB2632096A (en) | A method and a device for testing a model |
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: 24826022 Country of ref document: EP Kind code of ref document: A1 |
|
| ENP | Entry into the national phase |
Ref document number: 2025528143 Country of ref document: JP Kind code of ref document: A |