WO2015041312A1 - 血管内血流動態の画像処理方法及びシステム - Google Patents

血管内血流動態の画像処理方法及びシステム Download PDF

Info

Publication number
WO2015041312A1
WO2015041312A1 PCT/JP2014/074801 JP2014074801W WO2015041312A1 WO 2015041312 A1 WO2015041312 A1 WO 2015041312A1 JP 2014074801 W JP2014074801 W JP 2014074801W WO 2015041312 A1 WO2015041312 A1 WO 2015041312A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
blood flow
analysis
blood
imaging
Prior art date
Application number
PCT/JP2014/074801
Other languages
English (en)
French (fr)
Inventor
恭輔 鎌田
英昭 林
Original Assignee
国立大学法人旭川医科大学
インフォコム株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=52688959&utm_source=***_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=WO2015041312(A1) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by 国立大学法人旭川医科大学, インフォコム株式会社 filed Critical 国立大学法人旭川医科大学
Priority to US15/023,015 priority Critical patent/US10898088B2/en
Priority to CN201480060640.XA priority patent/CN105705084B/zh
Priority to JP2015537973A priority patent/JP5953437B2/ja
Priority to EP14846348.2A priority patent/EP3047796B1/en
Publication of WO2015041312A1 publication Critical patent/WO2015041312A1/ja
Priority to US17/128,718 priority patent/US20210106237A1/en

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/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0071—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by measuring fluorescence emission
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0263—Measuring blood flow using NMR
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0275—Measuring blood flow using tracers, e.g. dye dilution
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4887—Locating particular structures in or on the body
    • A61B5/489—Blood vessels
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7278—Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B2018/00315—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body for treatment of particular body parts
    • A61B2018/00345—Vascular system
    • A61B2018/00404—Blood vessels other than those in or around the heart
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2505/00—Evaluating, monitoring or diagnosing in the context of a particular type of medical care
    • A61B2505/05—Surgical care
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00—Medical imaging apparatus involving image processing or analysis
    • A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0265—Measuring blood flow using electromagnetic means, e.g. electromagnetic flowmeter
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M5/00—Devices for bringing media into the body in a subcutaneous, intra-vascular or intramuscular way; Accessories therefor, e.g. filling or cleaning devices, arm-rests
    • A61M5/007—Devices for bringing media into the body in a subcutaneous, intra-vascular or intramuscular way; Accessories therefor, e.g. filling or cleaning devices, arm-rests for contrast media
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56366—Perfusion imaging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10016—Video; Image sequence
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10048—Infrared image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20212—Image combination
    • G06T2207/20221—Image fusion; Image merging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • G06T2207/30104—Vascular flow; Blood flow; Perfusion
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing

Definitions

  • the present invention relates to a method and system for imaging blood flow dynamics in blood vessels.
  • CT Computer Tomography: computed tomography
  • MRI Magnetic Resonance Imaging: nuclear magnetic resonance imaging
  • a catheter The blood vessels are delineated by injecting a contrast agent from the patient, and the condition of the brain parenchyma, the movement of the blood vessels, and the blood flow are observed based on the images obtained by continuous imaging.
  • Perfusion Imaging the blood flow (perfusion) of tissue capillaries or a functional vasculature equivalent thereto, some mark (tracer) on the blood flow on the artery side, and the tracer goes into the blood flow.
  • An imaging method that observes passage through tissue which makes it possible to image tissue blood flow at the capillary level quantitatively or semi-quantitatively, and an image taken by such a method It is called a perfusion image.
  • Perfusion imaging is an analysis of tissues (capillaries).
  • the arrival time (AT: Arrival Time), time to maximum amplitude (TTP: Time To Peak), blood flow (BF: Blood Flow), blood volume (BV: Blood volume), mean transit time (MTT: Mean Transit Time), etc.
  • AT Arrival Time
  • TTP Time To Peak
  • BF Blood Flow
  • BV Blood volume
  • MTT Mean Transit Time
  • TDC Time Density Curve
  • CT Perfusion The method of imaging using a CT apparatus is called CT Perfusion
  • MRI Perfusion the method of imaging using an MRI apparatus.
  • Perfusion images using CT and MRI are rapidly injected with a contrast agent as a tracer to highlight and capture blood vessels, and the same cross-sectional image is continuously taken, and individual pixels of the obtained tomographic image are obtained. It can be obtained by calculating and analyzing TDC representing concentration change over time.
  • CT Perfusion and MRI Perfusion are both methods of analysis of tissue (capillaries).
  • CT Perfusion has X-ray exposure and is a pathological condition in which blood brain barrier (Blood Brain Barrier) is abnormal
  • blood brain barrier Blood Brain Barrier
  • MRI Perfusion has the disadvantage that quantitative data can not be obtained because patients who have been operated on for implanting magnetic bodies in the body are out of indication and linearity between contrast agent concentration and signal intensity can not be guaranteed.
  • CT Perfusion and MRI Perfusion can not be applied to patients who are allergic to contrast agents because they can not inject tracers intravenously.
  • Xe CT ⁇ How to use Xe CT>
  • Xe is diffused into brain tissue when Xe gas is inhaled and tomographic images are obtained sequentially with a CT apparatus, and tissue concentration (CT value) of CT image Take advantage of a slight rise.
  • CT value tissue concentration
  • TDC is calculated from the rise in CT value over time to image the blood flow state. Since Xe gas has excitatory action and anesthetic action and is difficult to use, it needs a closed circuit device to supply the gas, and consideration must also be given to X-ray exposure.
  • RI drug is injected or inhaled into the peripheral vein by a method using radioisotope-containing drug (hereinafter referred to as RI drug) as a tracer, and SPECT device (Single Photon Emission Computed) of the distribution of RI drug in the body over time
  • RI drug radioisotope-containing drug
  • SPECT device Single Photon Emission Computed
  • a catheter is guided to a target blood vessel from an artery such as the groin, elbow, wrist, etc., and a contrast agent is injected into the blood vessel and X-ray fluoroscopic imaging is performed to observe the travel of the blood vessel or the stenosis area.
  • the procedure can be used to treat simultaneously.
  • Images taken by digital subtraction angiography (DSA) can be observed with high contrast only for blood vessels injected with contrast agent.
  • a blood flow state can be imaged by setting an ROI to a blood vessel to be evaluated of a captured DSA image and analyzing TDC of a blood vessel (contrast agent) present in the set ROI.
  • Angiographic examination can easily determine the effect before and after endovascular surgery, but there is exposure to X-rays, and furthermore, patients who are allergic to contrast agents are not eligible for examination. Further, in order to directly insert a catheter into a blood vessel, the examination room and examination equipment are required to have a cleanliness level similar to that of the operating room. In addition, post-operative hemostasis requires absolute rest for several hours, and basically it is an examination that requires hospitalization.
  • Diagnosis in vascular diseases is performed not only by imaging blood flow conditions, but also by quantifying blood flow conditions and blood flow quantitatively and qualitatively.
  • the above information is calculated by analyzing the density change graph of the above-mentioned image with time, and using an ultrasonic Doppler wave using ultrasonic waves, and a blood flow measurement device by an electromagnetic blood flow meter (Electromagnetic Blood Flow Meter) There is a way.
  • ⁇ Method using ultrasonic Doppler> Using color Doppler imaging of an ultrasonic examination apparatus, it is possible to color blood actors in a living body in real time and obtain a blood circulation state by superimposing on a two-dimensional tomographic image B mode image it can.
  • the color Doppler method uses the fact that the frequency of the reflected sound wave changes due to the Doppler effect of ultrasonic waves, and determines whether the target object (blood) is approaching or away from the probe (probe). It is a technology to image Because the ultrasonic probe tip is large, it is impossible to evaluate small blood vessels, and blood flow measurement is impossible if the blood vessel diameter and the probe size are different, but it can be easily inspected anywhere and observation from multiple directions is possible. You can observe the results in real time.
  • the type of information that can be acquired by the device to be used also determines the availability of the inspection and the measurement method. For example, in an electromagnetic blood flow meter, instantaneous blood flow, average blood flow, stroke volume, etc. can be measured, but the output data is numerical data (graph), and in intraoperative fluorescence angiography, it is visually Although AT and TTP images can be estimated, it is impossible to create blood flow evaluation images to which BV, BF, MTT data and their information are added. Therefore, it is necessary for the operator to appropriately select an instrument capable of obtaining information necessary for the evaluation after the operation even from the measurable data contents to evaluate the operation.
  • Patent Document 1 discloses a method for evaluating the patency of a blood vessel that has undergone bypass graft surgery or the like using a fluorescence imaging agent that emits radiation of a specific wavelength.
  • BV Boood Volume: blood volume
  • BF Boood Flow: blood flow
  • MTT Mean Transit Time: average transit time
  • the image processing method for blood flow dynamics in blood vessels is an image output obtained by imaging a moving image with infrared light, with a part of a blood vessel into which a fixed amount of fluorescent contrast agent has been injected as an imaging target. It is characterized in that the shape of the luminance value time change curve is subjected to image analysis, and non-quantitative data of blood volume or blood flow is calculated based on the result of the image analysis.
  • the blood flow volume is measured by the electromagnetic blood flow meter with a part of the blood vessel as a measurement target, and the determination of the blood volume or the blood flow volume based on the measurement result and the non-quantitative data. Calculate the data.
  • an analysis image is generated from quantitative data of the blood volume or blood flow. Also, instead of the analysis image, an analysis moving image is generated.
  • the fluorescent contrast agent is indocyanine green or fluorescein.
  • a part of the blood vessel is imaged as a moving image by natural light, and an output image captured as a moving image by the natural light and an output image captured as a moving image as infrared light are fused.
  • the image processing system for blood flow dynamics in a blood vessel comprises an infrared imaging device for imaging a part of a blood vessel into which a predetermined amount of fluorescent contrast agent has been injected as a moving image by infrared light, and imaging using the imaging device And analyzing the shape of the luminance value time change curve of the image output, and calculating non-quantitative data of blood volume or blood flow based on the result of the image analysis.
  • the image analysis device includes an image analysis device that calculates quantitative data instead of non-quantitative data.
  • the apparatus further comprises an analysis image generation device that generates an analysis image from the quantitative data of the blood volume or blood flow volume.
  • a natural light imaging device for imaging a part of the blood vessel as a moving image by natural light, and an image for fusing an output image of the natural light imaging device and an analysis image generated from the analysis image generation device And a fusion device.
  • the image processing method and system for blood flow dynamics in blood vessels make it possible to estimate information such as BV, BF, MTT, blood vessel wall thickness and the like.
  • FIG. 5 is a flow diagram of a method that includes image fusion processing. It is a figure which shows the process sequence of an image fusion process.
  • the system according to the present embodiment converts an infrared light imaging device 100 that picks up a blood vessel as a moving image with infrared light, and a video signal output from the infrared light imaging device 100 to obtain a moving image.
  • a moving image conversion device 102 for creating a file an image analysis device 104 for analyzing a moving image file converted by the moving image conversion device 102, and calculating non-quantitative data of blood volume or blood flow based on the image analysis result,
  • an analysis image generation device for generating an analysis image from quantitative data of blood volume or blood flow obtained by the analysis device.
  • These devices may be provided, for example, inside one personal computer. Alternatively, they may be separate devices. Alternatively, some or all of the functions may be incorporated in the form of software in a personal computer, a measuring instrument, a display device, an analyzer, etc.
  • FIG. 2 shows an example of the configuration of a system according to this embodiment.
  • FIG. 2 an example of the configuration of the invention according to the present embodiment in micro-surgery (operation performed using a microscope) is described.
  • the implementation procedure of the method according to the present embodiment is as follows.
  • Indocyanine green is administered as a fluorescent contrast agent. Administration is by intravenous or transarterial routes. Example) Dilute 25 mg of indocyanine green to 10 ml and administer 2 ml once.
  • fluorescent contrast agent indocyanine green (ICG: Indocyanin green), fluorescein (FITC: Fluorescein isothiocyanate), etc.
  • ICG Indocyanin green
  • FITC Fluorescein isothiocyanate
  • the ICG polymerizes with .alpha.-lipopretein in blood to develop monochromatic fluorescence to infrared light. In the field of ophthalmology, this phenomenon is applied to fundus fluorescent angiography.
  • FITC is a fluorescent dye and emits green light when it is exposed to ultraviolet light.
  • any drug that emits light of a specific wavelength by irradiating light of a specific wavelength can be used appropriately.
  • Imaging is performed by the infrared imaging device 100.
  • the infrared imaging device 100 a general intraoperative microscope capable of infrared imaging can be used.
  • the analog (composite) video signal picked up by the infrared light imaging device 100 is analog-digital converted of the video signal using the video capture device of the video conversion device 102.
  • a digital video signal is captured as a moving image file on a local disk using a personal computer.
  • the moving image file may be a general purpose moving image file such as MPEG or MOV.
  • the image analysis device 104 analyzes the image output of the moving image pickup according to the following procedure.
  • fluorescent angiography agents are used during surgery in parallel to measuring quantitative BF values using an electromagnetic blood flow meter for one specific blood vessel that appears in the intraoperative field of view before vascular surgery.
  • the quantitative MTT data and the non-quantitative rBV and rBF data are calculated.
  • one point (one pixel) of a blood vessel in an image is monitored. This makes it possible to obtain information on the change in luminance value of one point (one pixel) of the blood vessel.
  • a target for obtaining information on change in luminance value may be one pixel, or information on an average value of luminance values of a plurality of pixels may be used.
  • the quantitative analysis shown in FIG. 3 is performed on all pixels of the moving image. The analysis results are calculated in time (X axis) units. AT, TTP, and MTT can be obtained as quantitative image data.
  • the analysis result is rBV and rBF, which are calculated from the integral value of the luminance value (Y axis), are obtained as non-quantitative image data.
  • ROI While viewing moving image data of a fluorescent contrast agent as a guide image, an ROI is created at the same blood vessel position as the specific blood vessel position measured by the electromagnetic blood flow meter (FIG. 4).
  • ROI statistical processing is performed on the rBF image.
  • An average non-quantitative blood flow value calculated by ROI statistical processing is defined as rBF (ROI).
  • a conversion factor value: Kbf that holds between eBF (ROI) and rBF (ROI) is calculated.
  • eBF (ROI) Kbf * rBF (ROI)
  • ROI creation is a region-of-interest creation process in which a region of interest on an image to be observed is surrounded by a two-dimensional closed curve.
  • ROI statistics is image processing which calculates the average image value per unit pixel inside ROI shape.
  • the output parameter image of the image analysis can be extracted from the shape of the luminance value time change curve obtained from the information of the change of the luminance value of one point of the blood vessel as shown in FIG.
  • the relationship is as shown in.
  • FIG. 1 An example of an analysis image before and after bypass surgery is shown in FIG.
  • the AT is shortened (blood flow has become faster) and the dark part has increased in the central part of the image.
  • the method according to the present embodiment it is possible to diagnose the condition of the blood vessel before or after the blood vessel surgery by measuring the quantitative BF value of the specific blood vessel using the electromagnetic blood flow meter, and in particular It can be provided in a form that is easy to visually recognize.
  • the area under the curve is defined as non-quantitative rBV, and non-quantitative rBV image for each pixel is easily You can ask for it.
  • the area centroid time when the AT time is defined as the origin zero [sec] is defined as the quantitative MTT.
  • a quantitative MTT image for each pixel can also be easily obtained.
  • Second Embodiment In the method of the second embodiment, the calculation method in the “quantitative analysis of moving image file” of the first embodiment is different, but other configuration of the system, other steps related to the method, and the like are the same as those of the first embodiment.
  • perfusion image method First Moment method for CT Perfusion and MRI Perfusion, which are analysis methods for tissues (capillary blood vessels), is used to analyze blood vessels using intraoperative fluorescence angiography (using a microscope camera) data.
  • the calculation method in the case of using in.
  • Non-quantitative rBF can be converted to a quantitative value into a BF [ml / min] image by measuring the quantitative BF [ml / min] of a specific blood vessel position using an electromagnetic blood flow meter.
  • BF [ml / min] Kbf * rBF: electromagnetic blood flow meter correction equation
  • BF [ml / min] is the “flow per unit tissue” quantitative blood flow (Blood Flow).
  • BF [ml / min] and BV [ml] can be directly calculated.
  • tissue capillary blood
  • BF [ml / min] and BV [ml] can be directly calculated.
  • blood vessel analysis it can be expected that the concept of "flowing around unit tissue" is not appropriate.
  • blood vessel analysis it is required to calculate BF [ml / min] and BV [ml] flowing through a single blood vessel.
  • Quantitative MTT is a calculated amount that can be defined by the same concept in both tissue analysis and blood vessel analysis. The area under the curve is defined as r ⁇ .
  • the rBF in FIG. 3 described above corresponds to r ⁇ .
  • Non-quantitative rMFV can be converted to a quantitative value into an MFV [cm / sec] image by measuring the quantitative MFV [cm / sec] of a specific blood vessel position using an electromagnetic blood flow meter.
  • MFV [cm / sec] Kbf * rMFV: electromagnetic blood flow meter correction formula
  • MFV [cm / sec] is an average blood flow velocity “per unit space”.
  • the space of the unit space is considered to indicate both blood vessels and tissues (capillaries).
  • BF [ml / min] a method of calculating quantitative BF is considered.
  • the following blood vessel diameter conversion formula holds between quantitative MFV and quantitative BF.
  • BF [ml / min] S [mm * mm] * MFV [cm / sec]: Vessel diameter conversion formula
  • S is defined as the inner blood vessel cross-sectional area at a certain blood vessel position.
  • This blood vessel diameter conversion formula means that if the inner blood vessel cross-sectional area S can be measured at the position of each blood vessel, the quantitative BF can be calculated. It is an important fact that the inner vessel cross-sectional area S is not a single value.
  • the thickness of the blood vessel differs for each blood vessel.
  • the thickness of the blood vessel differs depending on the position of even one blood vessel when strictly considered. These differences mean that the inner vessel cross-sectional area S is strictly given by the image.
  • the physical quantity obtained directly is non-quantitative rMFV.
  • the method of the third embodiment differs from that of the second embodiment in the calculation formula of BF, but the system configuration, other steps relating to the method, and the like are the same as those of the first embodiment.
  • the blood flow velocity (blood flow volume) at a certain blood vessel position is measured.
  • the blood vessel diameter at this one blood vessel position; R [mm] can be measured.
  • the inner blood vessel cross-sectional area: S is determined. In the first approximation, the inner vessel cross-sectional area: S approximates the same in the same single blood vessel measured by the electromagnetic blood flow meter.
  • the perfusion image (First Moment method) blood vessel analysis method can be approximately calculated.
  • Fourth Embodiment As a fourth embodiment, a method is described in which the same analysis as in the first to third embodiments can be performed by determining the operating conditions without performing the electromagnetic blood flow measurement.
  • the fourth embodiment is the same as the first to third embodiments except that the operation condition is determined without performing the electromagnetic blood flow measurement.
  • Quantitative BV and BF images can be approximately calculated by aligning the operating conditions without using an electromagnetic blood flow meter every operation.
  • the operating conditions include the sensitivity of the microscope camera, the magnification of the microscope camera, the working distance of the microscope camera, the angle of the microscope camera, the dose of the fluorescent contrast agent (ICG, FITC), etc. is there.
  • Working Distance is the distance from the tip of the objective lens of the microscope camera to the object in focus.
  • Angle is the angle at which the microscope camera is looking at the subject.
  • the dose of the fluorescent contrast agent refers to the amount of drug to be injected when injecting the fluorescent contrast agent (ICG, FITC) from the vein of the arm.
  • the fluorescent contrast agent is diluted to facilitate injection. It means equalizing the condition of dilution and injecting the same amount of injection after dilution.
  • the specific conditions for aligning the operating conditions may be determined for each institution such as a hospital, for example.
  • the quantitative BV and BF images can be obtained by using the electromagnetic blood flow meter to calculate the quantitative BV and BF images and aligning the operation conditions without using the electromagnetic blood flow meter. It can be made approximate calculation possible. For this reason, there are specific effects such as the blood flow volume which the bypass graft can supply to the tissue, and the prediction of the patency of the bypass over a long period.
  • the thickness information can provide information on the risk of surgical operations such as ease of rupture of the aneurysm and release of thrombus from the stenosis lesions.
  • the ICG can see 10 mm deep from its fluorescence frequency.
  • the blood vessel wall is about 0.1 to 10 mm, and the condition of the blood vessel wall is important at the time of surgery. Since the ICG brightness changes with the blood vessel wall thickness, the wall thickness can be estimated by keeping the dose, the microscope camera sensitivity, the distance to the object to be observed, and the magnification constant.
  • FITC can see 5 mm deep from its fluorescence frequency. Observation of vascular lesions having relatively thin blood vessel walls, patency of bypass, etc. is possible.
  • the purpose of this method is to know the blood vessel wall thickness from blood flow (BV) of blood vessel itself, blood volume (BF), mean transit time (MTT), and signal brightness, not brain tissue.
  • BV blood flow
  • BF blood volume
  • MTT mean transit time
  • signal brightness not brain tissue.
  • Image fusion processing With the output parameter image data of the above-mentioned fluorescence contrast agent analysis results alone, the relationship between anatomic positions such as arteries, veins, and cerebral sulcus is unclear.
  • the relationship of anatomic position can be grasped by performing image superposition processing (image fusion processing) of output parameter image data of a fluorescence contrast agent analysis result on anatomical image (moving image) data.
  • Image fusion processing is an analysis method of tissue (capillary blood vessels), but requires an analysis method of blood vessels.
  • the results of the above-described imaging of intravascular blood flow dynamics can be used.
  • a system relating to an image processing method of blood flow dynamics in blood vessels including image fusion processing is configured as shown in FIG.
  • the natural light imaging device 108 performs imaging as well as the imaging with the infrared light imaging device 100. Here, imaging with natural light is performed, and image or video data is recorded.
  • a natural light imaging device 108 for imaging a part of a blood vessel as a moving image by natural light, an output image of the natural light imaging device 108 and an analysis (moving image) image generated from an analysis image (moving image) generation device 106.
  • image (moving image) data of the following (A), (B) and (C) are used.
  • (gray scale image) means a value type image (used mainly in MRI and CT for medical images) in which the display color is applied by referring to the gray scale color table.
  • (Color scale image) is a meaning of a value type image (used mainly in nuclear medicine in medical images) to which a display color is applied by referring to a color table such as rainbow color.
  • (Color image) means an image whose color is fixed as shown in a photograph.
  • N pieces of moving image data mean an image consisting of N pieces of animation in the time axis direction.
  • Still image data is the meaning of an image such as a single picture.
  • the fluorescence contrast agent analysis according to the technique of the above embodiment 1 or 2 is a process of deriving (B) to (C).
  • the important point is that (A) is not used in the processing of fluorescence contrast agent analysis only, (B) and (C) are in the same positional relationship to create a still image from (B) moving image . Since the image fusion processes of (B) and (C) are at the same position, the positional relationship is always the same.
  • Whether (A) and (B) and (C) are in the same positional relationship depends on the device used. There are devices with the same position as the devices with different imaging positions of the natural light moving image and the infrared light moving image. Also, the number M of animations of natural light animations and the number N of animations of animations of fluorescent contrast agents are generally different. Further, it is different depending on the microscope device whether or not the timing at which the natural light moving image was captured and the timing at which the fluorescent contrast agent was captured were captured at the same time. It is not an essential condition of the image fusion process whether or not these are taken at the same time. In the above-described embodiment, the following three types of image fusion processing will be described.
  • Image fusion processing I multiplication type image fusion processing: Peak still image data (gray scale image) calculated from moving image data of fluorescent contrast agent and
  • C Output parameter still image data of image analysis (color scale image) Image fusion processing I to create image fusion still image data (color image) by image fusion processing. In this process, (one color image) is created from a combination of (one gray scale image) and (one color scale image).
  • Image fusion processing II multiplication type image fusion processing: (B) Moving image data of fluorescent contrast agent (gray scale image) and (C) Output parameter of image analysis Still image data (color scale image) image fusion Image fusion processing II that processes and creates image fusion moving image data (color image).
  • (N color images) are created from a combination of (N gray scale images) and (one color scale image).
  • Image fusion processing III Natural light moving image and image fusion processing: Image fusion moving image data (color image) created by image fusion processing and (A) moving image data of natural light (color image) or (A) still light of natural light Image fusion processing III for creating image fusion moving image data (color image) by image fusion processing of image data (color image).
  • (N color images) is created from a combination of (N color images) and (M or 1 color images).
  • FIG. 13 shows an example of a processing procedure of image superposition processing (image fusion processing). This processing procedure performs image fusion processing II and image fusion processing III, and generates (N color images) from a combination of (N color images) and (M or 1 color images). It is shown.
  • the image fusion processing I and the image fusion processing II perform the same multiplication type image fusion processing.
  • the difference between image fusion processing I and image fusion processing II is that (I) (grayscale image) is one sheet, while (II) is (grayscale image) N sheets.
  • (C) Output parameter of image analysis Still image data (color scale image) is one.
  • the image fusion processing I performs multiplication type image fusion processing only once.
  • Image fusion processing II is performed N times while changing the moving image data (gray scale image) of the fluorescent contrast agent.
  • the image fusion process I is a still image
  • the image fusion process II is a moving image.
  • the advantage of multiplication-type image fusion processing is that the colors of both grayscale and color images can be reliably reproduced.
  • the addition type image fusion processing can be calculated from the following basic equation.
  • Color of color scale image 1 (Red1, Green1, Blue1)
  • Color scale image 2 color (Red2, Green2, Blue2)
  • Image fusion processing result image color (FusionRed, FusionGreen, FusionBlue) ⁇ : synthesis ratio is a value represented by a number between 0.0 and 1.0.
  • FusionRed ⁇ * Red2 + (1.0- ⁇ ) * Red1
  • the addition type image fusion processing equation is named as the addition type image fusion processing equation because the colors of the two images are added.
  • the addition type image fusion processing is a method that can not reliably reproduce the colors (Red, Green, Blue) of the color scale image.
  • the multiplication-type image fusion processing expression is named as a multiplication-type image fusion processing expression in order to multiply the colors of two images.
  • the color (Red, Green, Blue) of the color scale image can be reliably reproduced by using the multiplication-type image fusion processing equation.
  • the gray scale image is black
  • the color scale image after fusion processing can always be black.
  • the position where the gray scale image is white makes the color scale image after the fusion process the same as the color scale image color (Red, Green, Blue) of the original image.
  • BF, BV, and MTT which could not be grasped conventionally by fluorescence contrast agent analysis, are stably grasped, and further, estimation of blood vessel wall thickness from signal luminance change, natural light and fusion (fusion) Provide anatomical orientation. Anatomical positional relationship can also be easily grasped by fusing the fluorescence contrast agent analysis result with the image captured under natural light.
  • BV derived from brain function response
  • BF change can also be detected.
  • monochrome moving picture editing may be applicable not only to fluorescence contrast agent analysis but also to gray scale moving pictures such as ordinary angiography. It is possible to determine the therapeutic effect immediately after blood vessel surgery by quantitatively comparing the state before blood vessel surgery and the state after blood vessel surgery using quantitative BV, BF, and MTT images.
  • the amount and timing of intravenous injection of the fluorescent contrast agent are important, and since the fluorescent contrast agent is metabolized by the liver, it also changes in the decay process of the brightness of the fluorescent contrast agent by the liver function.
  • the influence of the shift can be minimized even if the timing of injection of the fluorescent contrast agent is shifted. For this reason, it is an analysis method effective for therapeutic effect determination.
  • the method and system according to the present invention can be applied to neurosurgical surgery, orthopedic surgery, and vascular surgery in the field of ophthalmologic surgery where intraoperative microscopy is used, but is not particularly limited thereto, and various vascular surgery Can be applied to
  • Reference Signs List 100 infrared light imaging device 102 moving image conversion device 104 image analysis device 106 analysis image (moving image) generation device 108 natural light imaging device 110 image (moving image) fusion device

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Hematology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Vascular Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

【課題】　本発明では、術中顕微鏡で撮影される蛍光造影剤の動画データに関する解析技術が提供され、BV、BF、MTTなどの情報を推定可能なパーフュージョン（Perfusion）解析方法を蛍光造影剤解析へ応用する事で、蛍光造影剤解析でも、BV、BF、MTTなどの情報、血管壁厚などを推定することを可能とする方法およびシステムを提供することを課題とする。 【解決手段】 本発明に係る血管内血流動態の画像処理方法は、一定量の蛍光造影剤が注入された血管の一部を撮像対象として、赤外光により動画撮像をし、該動画撮像した画像出力の輝度値時間変化曲線の形状を画像解析し、該画像解析の結果に基づいて、血液量又は血流量の非定量データを算出することを特徴とする。

Description

血管内血流動態の画像処理方法及びシステム
　本発明は、血管内血流動態の画像処理方法及びシステムに関する。
　脳動脈瘤やもやもや病などの脳血管疾患において、その疾患の診断や治療結果の判定をするために、様々な検査が行われている。例えばCT（Computed Tomography：コンピュータ断層撮影）検査やMRI（Magnetic Resonance Imaging：核磁気共鳴画像法）検査では、断層画像やその画像を再構成した画像を利用し、血管造影検査(Angiography)では、カテーテルから造影剤を注入することで血管を描出し、連続で撮影することで得られる画像を元に、脳実質の状態や血管の走行、血液の流れを観察する。
　血管外科手術の診断において、臓器や血管の形態を画像として観察し、虚血部位や範囲、発症からの経過時間などを判定するだけでなく、血管に流れる血流の状態を画像化、数値化することで、疾患の診断や治療結果の効果を定量的、定性的に判定することが可能である。
　パーフュージョンイメージング(Perfusion Imaging)法は、組織の毛細血管あるいはそれに準ずる機能血管系の血流(灌流)に、動脈側で血流に何らかの目印(トレーサ)をつけ、トレーサが血流にのって組織を通過してゆく様子を観察する撮影方法で、毛細血管レベルの組織血流を定量的、あるいは半定量的に画像化することが可能となっており、このような方法で撮影した画像をパーフュージョンイメージと言う。パーフュージョンイメージング法は組織(毛細血管)の解析法である。
　また、血流の判定を行うための指標としては、到達時間（AT:Arrival Time）、最大振幅までの時間（TTP:Time To Peak）、血流量(BF：Blood Flow)、血液量（BV：Blood Volume）、平均通過時間（MTT：Mean Transit Time）などが利用される。例えば、診断画像を用いて血流量を定量評価する場合、トレーサ物質による画像の経時的濃度変化を示すTime Density Curve（TDC）を算出し、それを解析することで上記の評価値を算出する。さらに、これらの評価値及び、TDCの解析に基づき、BF、BV、MTTの各定性画像を作成し、臨床評価に利用している。
　血流状態を画像化・数値化することが可能な検査には、以下のようなものがある。
＜CT装置、若しくはMRI装置を利用する方法＞
　CT装置を用いて撮影する方法をCT Perfusion、MRI装置を用いて撮影する方法をMRI Perfusionと呼ぶ。CT及びMRIを用いたパーフュージョンイメージは、血管を強調して撮影するために造影剤をトレーサとして急速静注し、同一断面画像を連続で撮影を行い、得られた断層画像の個々のピクセルの経時的な濃度変化を表すTDCを算出し、解析することで得られる。CT Perfusion及びMRI Perfusionは共に組織(毛細血管)の解析法である。
　CT PerfusionはX線の被曝があり、血液脳関門(Blood Brain Barrier)に異常をみる病態で定量性を求める場合には特別の配慮が必要となるが、標準的なCT装置及び解析ソフトウェアがあれば簡単に測定することが可能である。一方、MRI Perfusionは体内に磁性体を埋め込む手術を行った患者は適応外で、造影剤濃度と信号強度の間の線形性が保証されないため、定量的データは得にくいという欠点がある。また、CT Perfusion、MRI Perfusionは共に造影剤にアレルギーを持った患者には、トレーサを静注することができないため、適応できない。
＜Xe CTを利用する方法＞
　非放射性Xe(Xenon)ガスをトレーサとして利用する方法で、Xeガスを吸入しCT装置にて断層像を継時的に得るとXeが脳組織に拡散され、CT画像の組織濃度(CT値)が軽度上昇することを利用する。この継時的なCT値の上昇からTDCを算出し血流状態を画像化する。
　Xeガスには興奮作用、麻酔作用があり利用しにくく、ガスを供給する閉鎖回路装置が必要で、X線被曝にも配慮が必要である。
＜SPECT装置、若しくはPET装置を利用する方法＞
　放射性同位元素を含んだ薬品(以後RI薬品と呼ぶ)をトレーサとして利用する方法で、RI薬品を抹消静脈に注射若しくは吸入し、RI薬品の体内分布を継時的にSPECT装置(Single Photon Emission Computed Tomography)やPET装置(Positron Emission Tomography)を利用して体外から計測することでTDCを算出し、血流状態を画像化する撮影方法である。
　PETを利用した検査は、現在利用可能なパーフュージョンイメージング法で最も定量性に優れており、酸素摂取率やエネルギー代謝率などを同時に測定できる利点がある。しかし、SPECT、PET撮影は共にRI薬品を扱うための核医学設備が必要で、検査による被曝も多い。
＜血管造影検査を利用する方法＞
　鼠径部、肘、手首などの動脈からカテーテルを目的血管に誘導し、造影剤を血管内に注入してX線透視撮影することで、血管の走行や狭窄部位などを観察する方法であり、この手法を用いて治療を同時に行うことができる。デジタル・サブトラクション・アンギオグラフィー法(DSA)で撮影され、造影剤が注入された血管のみを高いコントラストで観察可能となる。撮影されたDSA画像の評価したい血管にROIを設定し、設定されたROI内に存在する血管(造影剤)のTDCを解析することで、血流状態を画像化することができる。
　血管造影検査は、血管内手術前後の効果判定を簡単に行うことができるが、X線による被曝があり、さらに、造影剤にアレルギーを持った患者は、検査適応外となる。また、血管内にカテーテルを直接挿入するため、検査室及び検査機器には手術室に準ずる清潔度が要求される。また、術後は止血のため数時間は絶対安静が必要で、基本的には入院が必要な検査である。
　血管疾患における診断は、血流状態を画像化するだけでなく、血流状態、血流量を定量的、定性的に数値化することでも行われる。これらの情報は、前記の画像の経時的濃度変化グラフを解析して算出する方法と、超音波を利用した超音波ドップラー、電磁血流計(Electromagnetic Blood Flow Meter)による血流計測機器を利用する方法がある。
＜超音波ドップラーを利用した方法＞
　超音波検査装置のカラードップラー法(Color Doppler Imaging)を利用すると、リアルタイムに生体内の血行動体に色を付け、2次元断層画像であるBモード画像上に重ね合わせて血行状態を得ることができる。カラードップラー法は、超音波のドップラー効果により、反射した音波の周波数が変化することを利用し、目的の物体(血液)が探触子(プローブ)に近づいているのか、遠ざかっているのかを判定し画像化する技術である。
　超音波プローブチップは大きいため小さな血管の評価はできず、血流計測も血管径とプローブサイズが異なると不可能であるが、場所を選ばず簡単に検査ができ、多方向からの観察が可能で、リアルタイムに結果を観察することができる。
＜電磁血流計を利用した方法＞
　磁界中で導体が動くと起電力が発生するフレミングの法則に基づくもので、血流を電流と見なし、血流と磁界の両方の直角方向に起電力が発生することを利用して、瞬時血流量、平均血流量、一回拍出量などを測定する。
　電磁血流計を利用した方法は、手術中の血流量計測に用いられ、リアルタイムな計測結果を測定することができる。しかし電磁血流計は、プローブを直接血管に装着する必要があり、目的血管を露出させないと計測ができない。
　さらに、血流の評価は、血液の流れ自体を観察することでも行われる。
＜術中蛍光血管撮影法を利用した方法＞
　トレーサとして静脈内投与された蛍光血管造影剤から励起される近赤外領域の蛍光を、手術用顕微鏡の近赤外領域光下でビデオ撮影することで血流を評価する方法である。
　血管内の蛍光血管造影剤の状態(血流の有無)を把握することはできるが、血流状態、血流量など詳細に信号変化を解析することは不可能である。撮影には目的血管を手術用顕微鏡装置の近赤外光下でビデオ撮影する必要があるが、術中に血流状態をリアルタイムに確認することが可能である。
　上記のように、血管疾患における血流の評価には様々な種類の検査や測定法が存在するが、臨床現場においては、患者の病態や用途、診断タイミングに応じて、最も適切な検査、測定法を検討し、採用している。例えば、動脈瘤クリッピング後の血流状態の評価を行う場合、手術中にリアルタイムで血流の評価を行うことが要求される。この時、例えばCTやMRI、血管造影検査など、手術室に備え付けられていない専門の装置を必要とする検査は適応できず、電磁血流計や術中蛍光血管造影法を利用した血流の評価が選択される。
　さらに、使用する装置で取得可能な情報の種類によっても、その検査、測定法利用の可否が決まる。例えば、電磁血流計では瞬時血流量、平均血流量、一回拍出量などが測定可能であるが、出力されるデータは数値データ(グラフ)であり、術中蛍光血管撮影法では、目視でAT、TTP画像は推定可能であるが、BV、BF、MTTデータ及びそれらの情報を付加した血流評価画像の作成は不可能である。そのため、術者は術後の評価に必要な情報を得ることができる計器を、その測定可能なデータ内容からも適宜選択し、手術の評価を行う必要がある。つまり、例えばBV、BF、MTTの値、及び血流評価画像を得たい場合、例えば、患者の病状の落ち着いた手術数日後にCT Perfusion Imageを撮影し、さらに撮影した画像を解析する必要がある。つまり、手術中リアルタイムにBV、BF、MTTの値、及び血流評価画像を得るのは困難である。
　例えば、特許文献１には、特定の波長の放射線を放射する蛍光イメージング剤を用いて、バイパスグラフト手術などを受けた血管の開存性を評価する方法が開示されている。
特開２０１１－１４７７９７号公報
　しかしながら、従来の技術においては、術中蛍光血管撮影法からBV（Blood Volume:血液量）、BF（Blood Flow:血流量）、MTT（Mean Transit Time:平均通過時間）などの所望の情報を推定することは非常に困難であった。
　本発明では、術中顕微鏡で撮影される蛍光造影剤の動画データに関する解析技術が提供され、BV、BF、MTTなどの情報を推定可能なパーフュージョン（Perfusion）解析方法を蛍光造影剤解析へ応用する事で、蛍光造影剤解析でも、BV、BF、MTTなどの情報、血管壁厚などを推定することを可能とする方法およびシステムを提供することを課題とする。
 本発明に係る血管内血流動態の画像処理方法は、一定量の蛍光造影剤が注入された血管の一部を撮像対象として、赤外光により動画撮像をし、該動画撮像した画像出力の輝度値時間変化曲線の形状を画像解析し、該画像解析の結果に基づいて、血液量又は血流量の非定量データを算出することを特徴とする。
 本発明の一態様によれば、前記非定量データに代えて、前記血液量又は血流量の定量データを算出する。
 本発明の一態様によれば、前記血管の一部を測定対象として、電磁血流量計により血流量を測定し、この測定結果と前記非定量データとに基づいて、血液量又は血流量の定量データを算出する。
　本発明の一態様によれば、前記血液量又は血流量の定量データから解析画像を生成する。また、前記解析画像に代えて、解析動画を生成する。
 本発明の一態様によれば、前記蛍光造影剤がインドシアニングリーン、あるいはフルオレセインである。
 本発明の一態様によれば、前記血管の一部を自然光により動画撮像し、該自然光により動画撮像した出力画像と前記赤外光により動画撮像した出力画像とを融合させることを特徴とする。
 本発明に係る血管内血流動態の画像処理システムは、一定量の蛍光造影剤が注入された血管の一部を赤外光により動画として撮像する赤外光撮像装置と、該撮像装置により撮像された画像出力の輝度値時間変化曲線の形状を画像解析し、該画像解析の結果に基づいて血液量又は血流量の非定量データを算出する画像解析装置と、を備えることを特徴とする。
 本発明の一態様によれば、前記画像解析装置において、非定量データに代えて、定量データを算出する画像解析装置を備える。
 本発明の一態様によれば、前記血液量又は血流量の定量データから解析画像を生成する解析画像生成装置をさらに備える。
 本発明の一態様によれば、前記血管の一部を自然光により動画として撮像する自然光撮像装置と、該自然光撮像装置の出力画像と前記解析画像生成装置から生成された解析画像とを融合させる画像融合装置と、を備える。
　本発明に係る血管内血流動態の画像処理方法及びシステムでは、BV、BF、MTTなどの情報、血管壁厚などを推定することが可能となる。
第１実施形態に係る方法のフロー図である。 システム構成の一例を示す図である。 第１実施形態の画像解析における出力パラメータ画像の定義を説明する図である。 ROIの作成位置を説明する図である。 第１実施形態の画像解析における出力パラメータ画像の定義を説明する図である。 第２実施形態の画像解析における出力パラメータ画像の定義を説明する図である。 BFの算出方法を説明するフロー図である。 BFと血管断面積との関係を説明する図である。 第２実施形態のBFの算出方法を説明する図である。 第３実施形態のBFの算出方法を説明する図である。 バイパス手術の前後における解析画像の変化を示す図である。 画像フュージョン処理を含む方法のフロー図である。 画像フュージョン処理の処理手順を示す図である。
　以下に、本発明の血管内血流動態の画像処理方法及びシステムを実施するための形態について説明する。本実施形態では、脳血管のバイパス手術において、手術室中でBVおよびBF定量画像を算出して診断する方法を例として説明する。
　なお、本発明は以下の実施形態に限定されるものではないことは自明である。
　以下の実施形態において、用語と略語(解析出力画像名)は以下のように対応している。
　　　AT:　Arrival Time　到達時間　[sec]
　　　TTP:　Time To Peak　最大振幅までの時間　[sec]
　　　MTT:　Mean Transit Time　平均通過時間　[sec]
　　　MFV:　Mean Flow Velocity　定量血流速度　[cm/sec]
　　　BV:　Blood Volume　定量血液量　[ml]
　　　BF:　Blood　Flow　定量血流量　[ml/min]
　　　rMFV:　relative Mean Flow Velocity　非定量血流速度
　　　rBV:　relative Blood Volume　非定量血液量
　　　rBF:　relative Blood Flow　非定量血流量
　　　eBF:　電磁血流計を用いて計測した定量血流量　[ml/min]
　　　Peak:　最大振幅画像
　　　Kbf:　定量化換算係数
　　　S:　血管断面積(内側)　[mm*mm]
　　　Fusion：　画像重ね合わせ
　　　ROI:　Region of Interest　関心領域
　　　ICG:　Indocyanin green　インドシアニングリーン
　　　FITC:　Fluorescein isothiocyanate　フルオレセイン
　なお本実施形態において、画像出力とは、画像、動画、あるいは動画のコマから抽出された画像による出力を意味する。
＜第１実施形態＞
　まず、第１実施形態に係る血管内血流動態の画像処理のシステムの構成について説明する。本実施形態に係るシステムは、図１に示すように、血管を赤外光により動画として撮像する赤外光撮像装置１００と、赤外光撮像装置１００から出力されるビデオ信号を変換して動画ファイルを作成する動画変換装置１０２と、動画変換装置１０２によって変換された動画ファイルを解析し、該画像解析結果に基づいて血液量又は血流量の非定量データを算出する画像解析装置１０４と、画像解析装置１０４で得られた血液量又は血流量の定量データから解析画像を生成する解析画像生成装置１０６とを備える。
　これらの装置は、たとえば、一つのパーソナル・コンピュータ内部に設けられていてもよい。あるいは、個別の装置としてもよい。あるいは、ソフトウェアの形態でパーソナル・コンピュータや測定機器、表示装置、分析装置などに、その一部または全部の機能が組み込まれていてもよい。
　図２には、本実施形態に係るシステムの構成の一例が示されている。図２では、マイクロ・サージャリー（顕微鏡を使って行う手術）における、本実施形態に係る発明の構成の一例が記載されている。
　本実施形態に係る方法の実施手順は以下の通りである。
［蛍光造影剤の投与］
　蛍光造影剤としてインドシアニングリーンを投与する。経静脈、または経動脈経路で投与を行う。例）インドシアニングリーン25mgを10mlに希釈して一回2ml投与する。
　蛍光造影剤としては、インドシアニングリーン（ICG: Indocyanin green）、フルオレセイン（FITC: Fluorescein isothiocyanate）などを用いることができる。ICGは血液内でα－リポプレテインと重合することで赤外線光に対して、モノクロの蛍光を発色する。本現象を利用して眼科分野では眼底蛍光血管撮影に応用されている。FITCは蛍光色素であり、紫外線を当てると緑色の光を出す。
　上記の蛍光造影剤以外にも、たとえば、特定の波長の光を照射することにより、特定の波長の光を放出するような薬剤であれば適宜用いることができる。
［電磁血流計を用いた定量BF値の計測］
　電磁血流計(Electromagnetic Blood Flow Meter)を用いてある１つの特定の血管位置の定量血流値を計測する。この定量計測値をeBF(ROI)と定義する。実際の電磁血流計を用いた計測操作はバイパスグラフトを含む血管構造を電磁血流計プローブで挟んで行う。
［動画撮像］
　赤外光撮像装置１００で撮像を行う。赤外光撮像装置１００としては、赤外光撮影が可能な一般的な手術中顕微鏡を用いることができる。
［動画ファイル変換］
　赤外光撮像装置１００で撮像したアナログ(コンポジット)ビデオ信号を動画変換装置１０２のビデオキャプチャー装置を用いてビデオ信号のアナログ-デジタル変換する。デジタルビデオ信号はパソコンを用いてローカルディスクに動画ファイルとして取り込む。動画ファイルは、MPEGやMOVなどの汎用動画ファイルであってよい。
［動画ファイルの定量解析］
　画像解析装置１０４においては、以下の手順により、動画撮像した画像出力の解析を行う。
　ここでは、血管手術前に術中視野内に見える、ある１つの特定の血管に対して電磁血流計を用いた定量BF値を計測することと並行して、手術中に蛍光血管撮影薬を用いたPerfusion解析のFirst Moment法計算処理を行う事で、定量MTTデータと非定量rBV、rBFデータとを算出する。
　まず、時間経過とともに、画像の中の血管の1点(１ピクセル)をモニターする。これにより、当該血管の1点(１ピクセル)の輝度値の変化の情報を得ることができる。ここで、輝度値の変化の情報を得る対象を１ピクセルとしてもよいし、複数ピクセルの各輝度値の平均値の情報を用いてもよい。
　輝度値時間変化曲線を解析する事で、ICG薬剤が塊で血管を通過する状態を定量分析可能である。図３に示す定量解析を動画画像の全ピクセルで実行する。解析結果が時間(X軸)単位で算出されるAT、TTP、MTTは定量画像データとして求まる。解析結果が輝度値(Y軸)の積分値から算出されるrBV、rBFは非定量画像データとして求まる。
　蛍光造影剤の動画データをガイド画像として見ながら電磁血流計で計測した特定の血管位置と同一の血管位置にROIを作成する（図４）。次に、rBF画像に対するROI統計処理を行う。ROI統計処理により算出された平均非定量血流値をrBF(ROI)と定義する。eBF(ROI)とrBF(ROI)の間に成り立つ換算係数値：Kbfを算出する。
　　eBF(ROI) = Kbf * rBF(ROI)
　ここでROIとは、関心領域(ROI:Region of Interest)である。ROI作成とは、観測する画像上の関心領域を2次元閉曲線で囲む関心領域作成処理である。ROI統計とは、ROI形状内部の単位ピクセルあたりの平均画像値を算出する画像処理である。
　特定の血管位置(ROI)で成り立つ定量化換算式
　　eBF(ROI) = Kbf * rBF(ROI)
は、画像空間全***置に対しても成り立つ定量化換算式であると仮定する。
　すなわち、下記式のように、Perfusion解析のFirst Moment法計算処理から算出した非定量rBF画像全体に乗じる事で定量BFへ変換することができる。
　　BF = Kbf * rBF
　Perfusion解析のFirst Moment法計算で成り立つ関係式
　　rBF = rBF / MTT、BF = BV / MTT
と、非定量rBF画像を定量BF画像へ変換する定量化換算式
　　BF = Kbf * rBF
が成立する事から、下記のような非定量rBV画像を定量BV画像へ変換する定量化換算式が成り立つ。
　　BV = Kbf * rBV
　なお上記のパラメータは、以下に示す関係にある。
　　rBF = rBV / MTT　：　非定量関係式
　　BF = BV / MTT　：　定量関係式
　　eBF(ROI) = Kbf * rBF(ROI)　:　特定の血管位置(ROI)で成立する定量化換算式
　　BF = Kbf * rBF　：　画像全体で成立するrBFとBFの間の定量化換算式
　　BV = Kbf * rBV　：　画像全体で成立するrBVとBVの間の定量化換算式
　画像解析の出力パラメータ画像は、図５に示すような血管の1点の輝度値の変化の情報から得られた輝度値時間変化曲線の形状より抽出することができ、各パラメータの関係は図５に示すような関係にある。
　バイパス手術の前後における解析画像の一例を図１１に示す。この例では、バイパス手術後に、ATが短縮し（血流が早くなった）、画像中央部で色の濃い部分が増えた事がわかる。
　定量的な術中蛍光血管撮影解析を行うためには、蛍光造影剤の静脈注射をするタイミングにより、安定した結果を得ることが困難である。輝度値時間変化曲線から積分値を求め面積重心を求めることで、注射タイミングに依存しない安定したBF、BV、MTTを求めることができる。また、電磁血流計によるグラフト血流を参考にすることで、グラフト血管径、蛍光造影剤の投与量なども考慮して、より正確なBFを算出することもできる。必要に応じては、頚動脈に蛍光造影剤を注入することで一定のタイミングで検査を繰り返すことができる。
　上記のように、本実施形態に係る方法では、電磁血流計を用いて特定血管の定量BF値を計測することで、血管手術前あるいは後の血管の状態を診断することが可能となり、特に視覚的に認識しやすい形態で提供することができる。
　また、組織(毛細血管)の解析法として用いられる CT Perfusion及びMRI Perfusionのパーフュージョン解析のFirst Moment法では、カーブ下面積が非定量rBVと定義され、各ピクセル毎の非定量rBV画像は容易に求める事ができる。また、AT時間を原点ゼロ[sec]と定義した場合の面積重心時間が定量MTTと定義される。各ピクセル毎の定量MTT画像も容易に求める事ができる。
　従来のPerfusion解析方法では非定量rBVと非定量rBFの間には、rBF = rBV / MTT という計算式から成り立つ事が解っている。従来のPerfusion解析のFirst Moment法だけでは、定量BV、BF画像は推定不可能であるが、本実施形態に係る方法によれば、rBF = rBV / MTT という計算式を用いる事でピクセル毎の非定量rBV画像と定量MTT画像から非定量rBF画像を容易に求める事ができる。
＜第２実施形態＞
　第２実施形態の方法では、第１実施形態の「動画ファイルの定量解析」における計算手法が異なるが、システムの構成、方法に係る他のステップなどは第１実施形態と同様である。
　以下では、組織(毛細血管)の解析法であるCT Perfusion及びMRI Perfusion用のパーフュージョンイメージ法(First Moment法)を、術中蛍光血管撮影(顕微鏡カメラを用いた)データを用いて血管の解析法で使用した場合の計算手法について説明する。
　パーフュージョンイメージ法(First Moment法)を用いた組織(毛細血管)解析では、AT時間を原点ゼロ[sec]と定義した場合の面積重心時間が定量MTTと定義され、カーブ下面積が非定量rBVと定義される。非定量rBVと非定量rBFと定量MTTの間には、次の関係式が成立している。
　　rBF = rBV / MTT
　組織(毛細血管)解析法の電磁血流計補正式を考える。非定量rBFは電磁血流計を用いて特定血管位置の定量BF[ml/min]を計測する事で、BF[ml/min]画像へ定量値換算できる。
　　BF[ml/min] = Kbf * rBF　：　電磁血流計補正式
ここでBF[ml/min]は「単位組織あたりを流れる」定量血流量(Blood Flow)である。
　組織(毛細血管)解析からは、BF[ml/min]及びBV[ml]が直接算出できる。血管解析では「単位組織あたりを流れる」との概念は適切では無い事が予想できる。血管解析ではある1本の血管を流れるBF[ml/min]及びBV[ml]を算出する事が求められている。
　パーフュージョンイメージ法(First Moment法)を用いた組織(毛細血管)解析を血管解析に適応した場合の厳密解について考察する。血管解析ではAT時間を原点ゼロ[sec]と定義した場合の面積重心時間が定量MTTと定義される。この関係式は組織(毛細血管)解析法と変わらない。
　定量MTTは組織解析、血管解析共に同じ概念で定義できる算出量である。カーブ下面積をrβと定義する。非定量rβと定量MTTそれとrMFV：非定量平均血流速度(Mean Flow Velocity)の間には、
　　rMFV = rβ / MTT
関係式が成立している。この実施形態では、前述した図３のrBFがrβに相当する。
　上記の関係式の説明を図６に示す。この関係式は血管径を変えて行ったファントム実験のデータ解析から導き出す事が出来た。ファントム実験の結果rβ/MTTはrBFには比例していなかった。実験からは、「rβ/MTT」が「rBF/S」に比例する結論を示していた。
　血管解析法の電磁血流計補正式を考える。非定量rMFVは電磁血流計を用いて特定血管位置の定量MFV[cm/sec]を計測する事でMFV[cm/sec]画像へ定量値換算できる。
　　MFV[cm/sec] = Kbf * rMFV　：　電磁血流計補正式
　MFV[cm/sec]は「単位空間あたりを流れる」平均血流速度である。ここで単位空間の空間とは血管と組織(毛細血管)の両方を示していると考えられる。
　血管解析法の血管径換算式を考える。最終的に算出したい定量値はBF[ml/min]であるため定量BFを算出する方法を考える。定量MFVと定量BFの間には次の血管径換算式が成立する。
　 BF[ml/min] = S[mm*mm] * MFV[cm/sec]　：　血管径換算式
　BFの算出過程の概要を図７のフロー図に示す。
　また、図８に示すように、Sは、ある血管位置での内側血管断面積と定義する。
　この血管径換算式は、各血管の位置で内側血管断面積Sを測定できれば、定量BFを算出可能である事を意味している。内側血管断面積Sは1つの値で無い事が重要な事実である。血管の太さ(血管断面積)は血管毎に異なる。また厳密に考えると1本の血管でも位置が異なるとで血管の太さ(血管断面積)は異なる。これらの違いは内側血管断面積Sが厳密には画像で与えられる事を意味している。
　パーフュージョンイメージ(First Moment法)血管解析法を使用すると、直接得られる物理量は非定量rMFVである。この非定量rMFVから、電磁血流計補正式(MFV = Kbf * rMFV)と血管径換算式(BF= S * MFV)を用いる事で定量BFを算出することができる（図９）。
　第２実施形態の方法によって、より厳密なBFの値を算出することが可能となる。
＜第３実施形態＞
　第３実施形態の方法では、第２実施形態とはBFの算出式が異なるが、システムの構成、方法に係る他のステップなどは第１実施形態と同様である。
　BFを厳密に算出するためには、内側血管断面積:S[mm*mm]を厳密に計測する必要がある。しかし、現実には内側血管断面積:S[mm*mm]を厳密に計測する事は難しい。
　電磁血流計補正を行う段階では、ある１ヶ所の血管位置の血流速度(血流量)を測定している。この１ヶ所の血管位置での血管径;R[mm]は計測する事は可能である。計測した血管径;R[mm]から内側血管断面積:Sを求める。初めの近似は電磁血流計で計測した同じ１本の血管では内側血管断面積:Sは同一と近似する。
　電磁血流計では計測してない別の血管の内側血管断面積:Sについても同一と近似する。この近似により、厳密計測が難しかった内側血管断面積:S画像は、定数として取り扱う事が可能になる。
　第３実施形態の方法によって、パーフュージョンイメージ(First Moment法)血管解析法が近似的に算出可能となる。
＜第４実施形態＞
　第４実施形態として、電磁血流計測定を行わずとも、手術条件を定めることにより、第１から第３実施形態と同様の解析を行うことができる方法について説明する。第４実施形態においては、電磁血流計測定を行わずに手術条件を定めること以外は、上記第１から第３実施形態と同様である。
　電磁血流計測定を毎回行う事は簡単では無い。電磁血流計測定の計測操作は血管構造を損傷するリスクを含むため、できるかぎり回避したいのが現状である。手術毎に、電磁血流計を用いなくても手術条件をそろえる事で、定量BV、BF画像を近似的に算出できる。ここで、手術条件とは、顕微鏡カメラの感度、顕微鏡カメラの倍率、顕微鏡カメラのWorking Distance(ワーキング・ディスタンス)、顕微鏡カメラのAngle(アングル)、蛍光造影剤(ICG、FITC)の投与量などである。
　Working Distance(ワーキング・ディスタンス)とは顕微鏡カメラの対物レンズの先端から焦点があっている被写体までの距離の事である。Angle(アングル)とは、顕微鏡カメラが被写体を見ている角度の事である。
　蛍光造影剤の投与量とは、蛍光造影剤(ICG、FITC)を腕の静脈から注射するときの、注射する薬の量の事である。一般的には、蛍光造影剤を注射しやすいように希釈する。希釈の状態を同じくする事と希釈後の注射薬を同じ量だけ注射する事を意味する。なお手術条件をそろえるための具体的な条件は、例えば病院などの機関ごとに決定されても良い。
　第４実施形態に係る方法では、電磁血流計を用いる事で定量BV、BF画像を算出可能にする技術と電磁血流計を用いなくとも手術条件をそろえる事で、定量BV、BF画像を近似算出可能にすることができる。このため、バイパスグラフトが組織に供給できる血流量、長期にわたるバイパスの開存度の予測といった具体的な効果がある。
［血管壁厚の計測］
　脳動脈瘤、頚動脈狭窄病変においては、その厚さ情報は動脈瘤の破れやすさ、狭窄病変からの血栓遊離など手術操作に関わるリスクへの情報を提供できる。
　ICGはその蛍光周波数から10mmの深部を透見することができる。血管壁は0.1-10mm程度であり、手術の際にはその血管壁の性状が重要である。ICG輝度は血管壁厚により変化するため、投与量、顕微鏡カメラ感度、観察する対象物までの距離、倍率を一定に保つことで壁厚を推定することができる。
　FITCはその蛍光周波数から5mmの深部を透見することができる。比較的薄い血管壁を有している血管病変、バイパスの開存などの観察が可能である。
　本方法では脳の組織ではなく、血管自体の血流量(BV)、血液量(BF)、平均通過時間(MTT)、さらに信号輝度から血管壁厚を知ることを目的としている。
［画像フュージョン処理］　　　
　上記の蛍光造影剤解析結果の出力パラメータ画像データのみでは、動脈・静脈、脳溝などの解剖学的位置の関係は不明瞭である。解剖学的画像（動画）データの上に蛍光造影剤解析結果の出力パラメータ画像データを画像重ね合わせ処理(画像フュージョン処理)することで、解剖学的位置の関係を把握することできる。
　画像フュージョン処理は、組織(毛細血管)の解析法であるが、血管の解析法を必要としている。上記の血管内血流動態の画像処理の結果を利用することができる。画像フュージョン処理を含む血管内血流動態の画像処理方法に係るシステムは、図１２のような構成となっている。
　画像フュージョン処理を含む血管内血流動態の画像処理においては、赤外光撮像装置１００での撮像とともに、自然光撮像装置１０８でも撮像を行う。ここでは、自然光による撮像が行われ、画像または動画データが録画される。
　画像融合装置１１０では、血管の一部を自然光により動画として撮像する自然光撮像装置１０８と、該自然光撮像装置１０８の出力画像と解析画像（動画）生成装置１０６から生成された解析（動画）画像とを融合させる。
　画像フュージョン処理では、次の(A)、(B)および(C)の画像(動画)データを利用する。
(A)自然光のM枚動画データ(カラー画像)
(B)蛍光造影剤のN枚動画データ(グレースケール画像)
(C)解析結果の静止画データ(カラースケール画像)[BF、BV、MTTなどの複数のデータを含むものがある]。
　ここで(グレースケール画像)とは、グレースケールのカラーテーブルを参照する事で表示色を当てはめるバリュータイプ画像(主に医療画像ではMRI、CTで利用)の意味である。
(カラースケール画像)とは、レインボーカラー等のカラーテーブルを参照する事で表示色を当てはめるバリュータイプ画像(主に医療画像では核医学で利用)の意味である。
(カラー画像)とは写真のようにカラー色が固定されている画像の意味である。
　ここでN枚動画データとは、時間軸方向にN枚のアニメーションに成る画像の意味である。静止画データとは、１枚の写真のような画像の意味である。
　上記実施形態１または２の技術にかかる蛍光造影剤解析は、(B)から(C)を導き出す処理である。重要な点は、蛍光造影剤解析だけの処理では(A)は使わない点、(B)動画から(C)静止画を作り上げるため(B)と(C)は同じ位置関係にある点である。(B)と(C)の画像フュージョン処理は同一位置なので位置関係は必ず一致する。
　(A)と(B)、(C)が同じ位置関係にあるかは、使用する機器によって異なる。自然光動画と赤外光動画の撮像位置が異なる機器と位置が同じ機器が存在する。また、自然光動画のアニメーション枚数Mと蛍光造影剤の動画のアニメーション枚数Nは、一般的には異なる。また、自然光動画を撮影したタイミングと蛍光造影剤の動画を撮影したタイミングが同一時間に撮影したものか否かも顕微鏡機器によって異なる。これらを同一時間に撮影したか否かは、当該画像フュージョン処理の必須の条件ではない。上記の
　本実施形態では、次の３種類の画像フュージョン処理を説明する。
　(I)画像フュージョン処理I(掛け算型の画像フュージョン処理):蛍光造影剤の動画データから計算したPeak静止画データ(グレースケール画像)と(C)画像解析の出力パラメータ静止画データ(カラースケール画像)を画像フュージョン処理して画像フュージョン静止画データ(カラー画像)を作成する画像フュージョン処理I。
　この処理では、(1枚のグレースケール画像)と(1枚のカラースケール画像)の組み合わせから(1枚のカラー画像)を作成する。
　(II)画像フュージョン処理II(掛け算型の画像フュージョン処理):(B)蛍光造影剤の動画データ(グレースケール画像)と(C)画像解析の出力パラメータ静止画データ(カラースケール画像)を画像フュージョン処理して画像フュージョン動画データ(カラー画像)を作成する画像フュージョン処理II。
　この処理では、(N枚のグレースケール画像)と(1枚のカラースケール画像)の組み合わせから(N枚のカラー画像)を作成する。
　(III)画像フュージョン処理III(自然光動画と画像フュージョン処理):画像フュージョン処理で作成できた画像フュージョン動画データ(カラー画像)と(A)自然光の動画データ(カラー画像)または(A)自然光の静止画データ(カラー画像)を画像フュージョン処理して画像フュージョン動画データ(カラー画像)を作成する画像フュージョン処理III。
　この処理では、(N枚のカラー画像)と(M枚または1枚のカラー画像)の組み合わせから(N枚のカラー画像)を作成する。
　図１３には、画像重ね合わせ処理(画像フュージョン処理)の処理手順の一例が示されている。この処理手順は、画像フュージョン処理II及び画像フュージョン処理IIIを行い、(N枚のカラー画像)と(M枚または1枚のカラー画像)の組み合わせから(N枚のカラー画像)を作成する手順を示すものである。
　画像フュージョン処理Iと画像フュージョン処理IIは同じ掛け算型の画像フュージョン処理を行っている。画像フュージョン処理Iと画像フュージョン処理IIの違いは(I)は、(グレースケール画像)が1枚であるのに対して、(II)は(グレースケール画像)がN枚であることである。(C)画像解析の出力パラメータ静止画データ(カラースケール画像)は共に１枚である。
　画像フュージョン処理Iは掛け算型の画像フュージョン処理を１回だけ行っている。画像フュージョン処理IIは蛍光造影剤の動画データ(グレースケール画像)を変えながらN回分処理する。画像フュージョン処理の結果画像では画像フュージョン処理Iは静止画となり、画像フュージョン処理IIは動画となる。
　ここで、画像フュージョン処理Iと画像フュージョン処理IIで共通に使う掛け算型の画像フュージョン処理の詳細を説明する。掛け算型の画像フュージョン処理の利点は、グレースケール画像とカラー画像の双方の色を確実に再現できる点にある。
　以下に、掛け算型の画像フュージョン処理の具体的な処理手順をグレースケール画像1枚とカラースケール画像1枚の場合で説明する。はじめに、足し算型の画像フュージョン処理は最適でない事を説明する。
　足し算型の画像フュージョン処理は、次のような基本的な式から算出できる。
カラースケール画像1の色＝(Red1、Green1、Blue1)
カラースケール画像2の色＝(Red2、Green2、Blue2)
画像フュージョン処理結果画像の色=(FusionRed、FusionGreen、FusionBlue)
α:合成比率は、0.0～1.0の間の数字で表される値である。
　　FusionRed = α*Red2 + (1.0-α)*Red1
　　FusionGreen = α*Green2 + (1.0-α)*Green1
　　FusionBlue = α*Blue2 + (1.0-α)*Blue1
　足し算型の画像フュージョン処理式は、2つの画像の色を足しているため、足し算型の画像フュージョン処理式と名付けた。
　足し算型の画像フュージョン処理式で、グレースケール画像1枚とカラースケール画像1枚の場合は、
カラースケール画像1の色＝(Red1、Green1、Blue1)を
グレースケール画像の色＝(Gray、Gray、Gray)へ置き換える。
カラースケール画像2の色＝(Red2、Green2、Blue2)を
カラースケール画像の色＝(Red、Green、Blue)へ置き換える。
　　FusionRed = α*Red + (1.0-α)*Gray
　　FusionGreen = α*Green + (1.0-α)*Gray
　　FusionBlue = α*Blue + (1.0-α)*Gray
　足し算型の画像フュージョン処理式では、カラースケール画像の色(Red、Green、Blue)を再現出来ていない事が解る。
　この様に、足し算型の画像フュージョン処理は、カラースケール画像の色(Red、Green、Blue)を確実には再現出来ない手法である。これに対して、掛け算型の画像フュージョン処理式は、２つの画像の色を掛け算するため、掛け算型の画像フュージョン処理式と名付けた。
　掛け算型の画像フュージョン処理式を用いれば、カラースケール画像の色(Red、Green、Blue)を確実に再現することができる。グレースケール画像が黒い位置は、フュージョン処理後のカラースケール画像を必ず黒くできる。グレースケール画像が白い位置は、フュージョン処理後のカラースケール画像を元画像のカラースケール画像の色(Red、Green、Blue)と同じくできる。
　自然光のM枚動画データの中から、特定の時間の自然光の1枚静止画だけを使用するケースがある。動画にすると複雑な場合など、自然光の1枚静止画だけを用いる事で簡単化できる。
　以上に本発明の実施形態を説明した。第１実施形態の方法の一部を、第４実施形態に置き換えることが可能である。また、第１実施形態の方法の一部を、第２実施形態又は第３実施形態に置き換えることが可能である。すなわち、上記の実施形態を組み合わせることが可能である。
　以上説明したように、血管外科手術においては、術中に血流状態を把握できることが極めて重要である。脳血流遮断が完璧か、吻合血管の提供する血液量が十分かなどの血管手術の治療効果判定が重要である。上記に説明した実施形態に係る発明では、蛍光造影剤解析では従来把握できなかったBF、BV、MTTを安定して把握、さらに信号輝度変化より血管壁厚の推定、自然光と融合(fusion)による解剖学的オリエンテーションを提供する。蛍光造影剤解析結果を自然光の下で撮像された画像に融合(fusion)することで、解剖学的な位置関係も容易に把握することができる。
　特に、グラフトの長期開存度の予測、グラフトからの血流評価を含む過灌流に陥るリスクの予測が可能となり、脳表の血流変化状態を把握することで、脳機能反応に由来するBV、BF変化の検出もできる。さらにモノクロ動画編集は蛍光造影剤解析のみならず、通常の血管撮影のようなグレースケール動画も同様な解析に応用できる可能性がある。血管手術前の状態と血管手術後の状態を定量BV、BF、MTT画像を用いて定量比較する事で血管手術の直後に治療効果判定を行う事が可能である。
　従来の手法においては、蛍光造影剤の静脈注射の量とタイミングが重要であり、蛍光造影剤は肝臓で代謝されるため、肝機能によりその蛍光造影剤の輝度の減衰経過にも変化する。上記実施形態に係る発明から算出される定量BV、BF、MTT画像は蛍光造影剤を注射するタイミングがずれてもずれの影響は最小限にすることができる。このため、治療効果判定に有効な解析方法に成っている。
　また、バイパス手術後に血流の過潅流による症状が出る場合がある。このようなケースでは、バイパス手術後のBV、BF定量画像を観測する事で予測診断が可能である。過灌流は出血、痙攣などのリスクがあるため、術中からこの症候を予測できことは極めて有用である。従来の技術で提供された、AT、TTPのみの情報では、この症候の予測診断が困難であったが、上記実施形態に係る発明ではこのような予測にも対応することができる。
　また、フュージョン画像を利用することで、その領域について過灌流が起こりやすいかをも予測することができる。
　なお本発明に係る方法およびシステムは、術中顕微鏡が使用される脳神経外科手術、整形外科手術、眼科手術分野の血管外科手術に適応できるが、特にこれらに限定されるものではなく、多様な血管手術に応用することができる。
１００　赤外光撮像装置
１０２　動画変換装置
１０４　画像解析装置
１０６　解析画像（動画）生成装置
１０８　自然光撮像装置
１１０　画像（動画）融合装置

Claims (12)

  1.  一定量の蛍光造影剤が注入された血管の一部を撮像対象として、赤外光により動画撮像をし、
     該動画撮像した画像出力の輝度値時間変化曲線の形状を画像解析し、
     該画像解析の結果に基づいて、血液量又は血流量の非定量データを算出することを特徴とする血管内血流動態の画像処理方法。
  2.  前記非定量データに代えて、前記血液量又は血流量の定量データを算出することを特徴とする請求項１に記載の血管内血流動態の画像処理方法。
  3.  前記血管の一部を測定対象として、電磁血流量計により血流量を測定し、
     この測定結果と前記非定量データとに基づいて、血液量又は血流量の定量データを算出することを特徴とする請求項１に記載の血管内血流動態の画像処理方法。
  4.  前記血液量又は血流量の定量データから解析画像を生成することを特徴とする請求項１ないし３のいずれかに記載の血管内血流動態の画像処理方法。
  5.  前記解析画像に代えて、解析動画を生成することを特徴とする請求項４に記載の血管内血流動態の画像処理方法。
  6.  前記血管の一部を自然光により動画撮像し、該自然光により動画撮像した出力画像と前記解析画像とを融合させることを特徴とする請求項１ないし４のいずれかに記載の血管内血流動態の画像処理方法。
  7.  前記蛍光造影剤がインドシアニングリーンである請求項１ないし６のいずれかに記載の血管内血流動態の画像処理方法。
  8.  前記蛍光造影剤がフルオレセインである請求項１ないし７のいずれかに記載の血管内血流動態の画像処理方法。
  9.  一定量の蛍光造影剤が注入された血管の一部を赤外光により動画として撮像する赤外光撮像装置と、
     該撮像装置により撮像された画像出力の輝度値時間変化曲線の形状を画像解析し、該画像解析の結果に基づいて血液量又は血流量の非定量データを算出する画像解析装置と、を備えることを特徴とする血管内血流動態の画像処理システム。
  10.  前記画像解析装置において、非定量データに代えて、定量データを算出する画像解析装置を備えることを特徴とする請求項９に記載の血管内血流動態の画像処理システム。
  11.  前記血液量又は血流量の定量データから解析画像を生成する解析画像生成装置をさらに備えることを特徴とする請求項１０に記載の血管内血流動態の画像処理システム。
  12.  前記血管の一部を自然光により動画として撮像する自然光撮像装置と、
     該自然光撮像装置の出力画像と前記解析画像生成装置から生成された解析画像とを融合させる画像融合装置と、を備えることを特徴とする請求項１１に記載の血管内血流動態の画像処理システム。
PCT/JP2014/074801 2013-09-20 2014-09-19 血管内血流動態の画像処理方法及びシステム WO2015041312A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
US15/023,015 US10898088B2 (en) 2013-09-20 2014-09-19 Method and system for image processing of intravascular hemodynamics
CN201480060640.XA CN105705084B (zh) 2013-09-20 2014-09-19 血管内血流动态的图像处理方法及***
JP2015537973A JP5953437B2 (ja) 2013-09-20 2014-09-19 血管内血流動態の画像処理方法及びシステム
EP14846348.2A EP3047796B1 (en) 2013-09-20 2014-09-19 Method and system for image processing of intravascular hemodynamics
US17/128,718 US20210106237A1 (en) 2013-09-20 2020-12-21 Method and system for image processing of intravascular hemodynamics

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2013-194898 2013-09-20
JP2013194898 2013-09-20

Related Child Applications (2)

Application Number Title Priority Date Filing Date
US15/023,015 A-371-Of-International US10898088B2 (en) 2013-09-20 2014-09-19 Method and system for image processing of intravascular hemodynamics
US17/128,718 Continuation US20210106237A1 (en) 2013-09-20 2020-12-21 Method and system for image processing of intravascular hemodynamics

Publications (1)

Publication Number Publication Date
WO2015041312A1 true WO2015041312A1 (ja) 2015-03-26

Family

ID=52688959

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/074801 WO2015041312A1 (ja) 2013-09-20 2014-09-19 血管内血流動態の画像処理方法及びシステム

Country Status (5)

Country Link
US (2) US10898088B2 (ja)
EP (1) EP3047796B1 (ja)
JP (1) JP5953437B2 (ja)
CN (1) CN105705084B (ja)
WO (1) WO2015041312A1 (ja)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105559810A (zh) * 2015-12-10 2016-05-11 上海交通大学 血管单位时间血流量与血流速度的计算方法
JP2018534965A (ja) * 2015-09-23 2018-11-29 ノバダック テクノロジーズ ユーエルシー 組織の治癒を評価するための方法およびシステム
EP3416547A4 (en) * 2016-02-16 2019-08-14 Novadaq Technologies ULC FACILITATION OF EVALUATION OF BLOOD FLOW AND TISSUE PERFUSION BY FLUORESCENCE-MEDIATED PHOTOPLETHYSMOGRAPHY
KR20210114167A (ko) * 2020-03-10 2021-09-23 부산대학교 산학협력단 실시간 영상을 이용한 대장 혈류속도 분석방법
US20220280049A1 (en) * 2019-08-01 2022-09-08 Osaka University Blood vessel wall thickness estimation method, blood vessel wall thickness estimation device, and blood vessel wall thickness estimation system
US11439306B2 (en) * 2017-12-28 2022-09-13 Leica Instruments (Singapore) Pte. Ltd. Apparatus and method for measuring blood flow direction using a fluorophore
JP2022544052A (ja) * 2019-07-31 2022-10-17 ▲蘇▼州▲潤▼▲邁▼▲徳▼医▲療▼科技有限公司 造影画像間隔時間に基づいて血流速度を補正する方法及び装置
WO2023063324A1 (ja) * 2021-10-12 2023-04-20 国立大学法人東京工業大学 血液ポンプシステム及び血液循環システム
JP2023516882A (ja) * 2020-02-03 2023-04-21 カール ツアイス メディテック アクチエンゲゼルシャフト 術野内の血管の部分を通る血液の容積流量を特定するコンピュータ実施方法、コンピュータプログラム、及び手術システム
US11696695B2 (en) 2013-06-14 2023-07-11 Stryker European Operations Limited Quantification of absolute blood flow in tissue using fluorescence mediated photoplethysmography

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105705084B (zh) 2013-09-20 2019-07-12 国立大学法人旭川医科大学 血管内血流动态的图像处理方法及***
US10586314B2 (en) * 2016-04-12 2020-03-10 Shenzhen Everbest Machinery Industry Co., Ltd Image fusion method, apparatus, and infrared thermal imaging device
KR102571187B1 (ko) 2016-12-09 2023-08-28 퍼퓨전 테크 에이피에스 해부학적 구조에서 관류를 평가하기 위한 시스템 및 방법
CN110177505B (zh) * 2017-01-11 2023-03-14 株式会社岛津制作所 荧光成像装置和荧光成像***
JP6708143B2 (ja) * 2017-02-07 2020-06-10 株式会社島津製作所 時間強度曲線測定装置
JP7481264B2 (ja) * 2018-05-17 2024-05-10 ロンドン ヘルス サイエンシーズ センター リサーチ インコーポレイテッド 動的血管造影画像化
KR20210044767A (ko) 2018-06-14 2021-04-23 퍼퓨전 테크 에이피에스 자동 관류 측정 시스템 및 방법
US11911139B2 (en) 2018-06-14 2024-02-27 Perfusion Tech Aps System and method for automatic perfusion measurement
US11253165B2 (en) 2019-04-24 2022-02-22 Mohammad Nasser Intravascular MRI probe assembly
CN111429457B (zh) * 2020-06-03 2020-09-11 中国人民解放军总医院 图像局部区域亮度智能评价方法、装置、设备及介质
JP7274803B2 (ja) * 2020-11-04 2023-05-17 株式会社レイマック 非接触血管解析装置
CN113616226B (zh) * 2021-09-14 2023-06-23 上海联影医疗科技股份有限公司 一种血管分析方法、***、设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008220926A (ja) * 2007-03-09 2008-09-25 Korea Advanced Inst Of Science & Technology 下肢潅流量測定装置
JP2010017396A (ja) * 2008-07-11 2010-01-28 Olympus Corp 生体観察方法および装置
JP2011147797A (ja) 1999-09-24 2011-08-04 Natl Research Council Of Canada 術中血管造影を行なうための方法および装置
JP2013003495A (ja) * 2011-06-21 2013-01-07 Mitaka Koki Co Ltd 顕微鏡システム
JP2013039223A (ja) * 2011-08-16 2013-02-28 Fujifilm Corp 画像表示方法および装置

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5150292A (en) * 1989-10-27 1992-09-22 Arch Development Corporation Method and system for determination of instantaneous and average blood flow rates from digital angiograms
DE4214402C2 (de) * 1992-04-30 1997-04-24 Pulsion Verwaltungs Gmbh & Co Vorrichtung zum Bestimmen des Füllungszustandes eines Blutkreislaufs
US5865757A (en) 1994-12-01 1999-02-02 Hoeft; Andreas Process and a device for determining the cerebral blood flow CBF, and the intracerebral blood volume ICBV
JP3799173B2 (ja) 1998-09-10 2006-07-19 キヤノン株式会社 眼底血流計測装置
US6844465B2 (en) * 1998-10-30 2005-01-18 Ajinomoto Co., Inc. Method for preparing highly stable crystals of aspartame derivative
US20050182434A1 (en) 2000-08-11 2005-08-18 National Research Council Of Canada Method and apparatus for performing intra-operative angiography
US6915154B1 (en) 1999-09-24 2005-07-05 National Research Council Of Canada Method and apparatus for performing intra-operative angiography
DE10120980B4 (de) 2001-05-01 2009-12-03 Pulsion Medical Systems Ag Verfahren, Vorrichtung und Computerprogramm zur Bestimmung des Blutflusses in einer Gewebe- oder Organregion
US8620410B2 (en) * 2002-03-12 2013-12-31 Beth Israel Deaconess Medical Center Multi-channel medical imaging system
JP3920140B2 (ja) * 2002-05-13 2007-05-30 株式会社東芝 Ｍｒｉ装置及びフロー定量化装置
JP2007525250A (ja) * 2003-06-02 2007-09-06 ザ・ジェネラル・ホスピタル・コーポレイション 組織血流の遅延−補償型算出方法
WO2007004026A1 (en) * 2005-06-30 2007-01-11 Ge Healthcare Limited Evaluating cardiac function with dynamic imaging techniques and contrast media
DE102005059520B4 (de) 2005-12-13 2019-06-13 Iprm Intellectual Property Rights Management Ag Dilutionsvorrichtung und Computerprogramm
EP2173425B1 (en) * 2007-07-18 2012-11-21 Silk Road Medical, Inc. Systems for establishing retrograde carotid arterial blood flow
DE102008040804B4 (de) * 2008-07-28 2021-07-29 Carl Zeiss Meditec Ag Verfahren, Operationsmikroskop und Analysesystem zur quantitativen Darstellung des Blutflusses
DE102008040803A1 (de) * 2008-07-28 2010-02-04 Carl Zeiss Surgical Gmbh Verfahren zur quantitativen Darstellung des Blutflusses
DE102009022497A1 (de) * 2009-05-25 2010-12-02 Olympus Winter & Ibe Gmbh Operationsleuchte, Bildgebungssystem und Verwendung einer Operationsleuchte
US20110028850A1 (en) 2009-07-28 2011-02-03 Thomas Schuhrke Process for quantitative display of blood flow
JP5794597B2 (ja) 2011-07-09 2015-10-14 ガウス サージカルＧａｕｓｓ Ｓｕｒｇｉｃａｌ 体外血液量の推定及び外科的摘出サンプルの計数システム及び方法
CN105705084B (zh) 2013-09-20 2019-07-12 国立大学法人旭川医科大学 血管内血流动态的图像处理方法及***

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011147797A (ja) 1999-09-24 2011-08-04 Natl Research Council Of Canada 術中血管造影を行なうための方法および装置
JP2008220926A (ja) * 2007-03-09 2008-09-25 Korea Advanced Inst Of Science & Technology 下肢潅流量測定装置
JP2010017396A (ja) * 2008-07-11 2010-01-28 Olympus Corp 生体観察方法および装置
JP2013003495A (ja) * 2011-06-21 2013-01-07 Mitaka Koki Co Ltd 顕微鏡システム
JP2013039223A (ja) * 2011-08-16 2013-02-28 Fujifilm Corp 画像表示方法および装置

Non-Patent Citations (1)

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

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11696695B2 (en) 2013-06-14 2023-07-11 Stryker European Operations Limited Quantification of absolute blood flow in tissue using fluorescence mediated photoplethysmography
JP2018534965A (ja) * 2015-09-23 2018-11-29 ノバダック テクノロジーズ ユーエルシー 組織の治癒を評価するための方法およびシステム
JP2020171727A (ja) * 2015-09-23 2020-10-22 ノバダック テクノロジーズ ユーエルシー 組織の治癒を評価するための方法およびシステム
CN105559810A (zh) * 2015-12-10 2016-05-11 上海交通大学 血管单位时间血流量与血流速度的计算方法
EP3416547A4 (en) * 2016-02-16 2019-08-14 Novadaq Technologies ULC FACILITATION OF EVALUATION OF BLOOD FLOW AND TISSUE PERFUSION BY FLUORESCENCE-MEDIATED PHOTOPLETHYSMOGRAPHY
US10646128B2 (en) 2016-02-16 2020-05-12 Novadaq Technologies ULC Facilitating assessment of blood flow and tissue perfusion using fluorescence-mediated photoplethysmography
US11701016B2 (en) 2016-02-16 2023-07-18 Stryker European Operations Limited Facilitating assessment of blood flow and tissue perfusion using fluorescence-mediated photoplethysmography
US11439306B2 (en) * 2017-12-28 2022-09-13 Leica Instruments (Singapore) Pte. Ltd. Apparatus and method for measuring blood flow direction using a fluorophore
JP2022544052A (ja) * 2019-07-31 2022-10-17 ▲蘇▼州▲潤▼▲邁▼▲徳▼医▲療▼科技有限公司 造影画像間隔時間に基づいて血流速度を補正する方法及び装置
JP7260218B2 (ja) 2019-07-31 2023-04-18 ▲蘇▼州▲潤▼▲邁▼▲徳▼医▲療▼科技有限公司 造影画像間隔時間に基づいて血流速度を補正する方法及び装置
US20220280049A1 (en) * 2019-08-01 2022-09-08 Osaka University Blood vessel wall thickness estimation method, blood vessel wall thickness estimation device, and blood vessel wall thickness estimation system
JP2023516882A (ja) * 2020-02-03 2023-04-21 カール ツアイス メディテック アクチエンゲゼルシャフト 術野内の血管の部分を通る血液の容積流量を特定するコンピュータ実施方法、コンピュータプログラム、及び手術システム
JP7454686B2 (ja) 2020-02-03 2024-03-22 カール ツアイス メディテック アクチエンゲゼルシャフト 術野内の血管の部分を通る血液の容積流量を特定するコンピュータ実施方法、コンピュータプログラム、及び手術システム
KR102386872B1 (ko) * 2020-03-10 2022-04-13 부산대학교 산학협력단 실시간 영상을 이용한 대장 혈류속도 분석방법
KR20210114167A (ko) * 2020-03-10 2021-09-23 부산대학교 산학협력단 실시간 영상을 이용한 대장 혈류속도 분석방법
WO2023063324A1 (ja) * 2021-10-12 2023-04-20 国立大学法人東京工業大学 血液ポンプシステム及び血液循環システム

Also Published As

Publication number Publication date
US20210106237A1 (en) 2021-04-15
JPWO2015041312A1 (ja) 2017-03-02
US10898088B2 (en) 2021-01-26
EP3047796A4 (en) 2017-05-31
JP5953437B2 (ja) 2016-07-20
EP3047796A1 (en) 2016-07-27
EP3047796B1 (en) 2020-11-04
CN105705084A (zh) 2016-06-22
US20160262638A1 (en) 2016-09-15
CN105705084B (zh) 2019-07-12

Similar Documents

Publication Publication Date Title
US20210106237A1 (en) Method and system for image processing of intravascular hemodynamics
US11694339B2 (en) Method and apparatus for determining blood velocity in X-ray angiography images
Schroeder et al. Noninvasive detection and evaluation of atherosclerotic coronary plaques with multislice computed tomography
Schoepf et al. CT of coronary artery disease
Bauer et al. Noncalcified atherosclerotic plaque burden at coronary CT angiography: a better predictor of ischemia at stress myocardial perfusion imaging than calcium score and stenosis severity
Ohno et al. Oxygen-enhanced MR imaging: correlation with postsurgical lung function in patients with lung cancer
WO2007113705A1 (en) Determining tissue surrounding an object being inserted into a patient
Gremse et al. Virtual elastic sphere processing enables reproducible quantification of vessel stenosis at CT and MR angiography
JP2004320771A (ja) ディジタルサブトラクション血管造影法を実施するための方法
Temov et al. Coronary computed tomography angiography investigation of the association between left main coronary artery bifurcation angle and risk factors of coronary artery disease
Nishida et al. Model-based iterative reconstruction for multi–detector row CT assessment of the Adamkiewicz artery
JP2021522941A (ja) 動的血管造影画像化
Li et al. Evaluation of the early enhancement of coronary atherosclerotic plaque by contrast-enhanced MR angiography
Counseller et al. Recent technologies in cardiac imaging
Hetterich et al. X-ray phase-contrast computed tomography of human coronary arteries
Pereira et al. Non-invasive imaging techniques and assessment of carotid vasa vasorum neovascularization: Promises and pitfalls
JP6981807B2 (ja) 医用情報処理装置、ｘ線ｃｔ装置、医用情報処理プログラム、医用情報処理方法及び医用情報処理システム
Knollmann et al. Quantification of atherosclerotic coronary plaque components by submillimeter computed tomography
Hansch et al. Quantitative evaluation of MR perfusion imaging using blood pool contrast agent in subjects without pulmonary diseases and in patients with pulmonary embolism
TWI770235B (zh) 以動態影像判斷血流量變化及血管阻塞區域之方法
Zhang et al. 3D Cardiac microvessels embolization imaging based on X-ray phase contrast imaging
Mišič et al. Morphological Characteristics and Frequency of Accessory Right Hepatic Veins–Evaluation with Computed Tomography
Zheng et al. Accuracy of MDCT for detection and identification of carotid atherosclerotic plaque in a rabbit model
Ruzsics et al. Automated multidetector computed tomography evaluation of subacutely infarcted myocardium
Liu et al. Introduction to Imaging Methods

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

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2015537973

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

REEP Request for entry into the european phase

Ref document number: 2014846348

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2014846348

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 15023015

Country of ref document: US