[go: up one dir, main page]

JP2007135894A - Blood flow analysis device and simulation method based on human blood flow data - Google Patents

Blood flow analysis device and simulation method based on human blood flow data Download PDF

Info

Publication number
JP2007135894A
JP2007135894A JP2005334325A JP2005334325A JP2007135894A JP 2007135894 A JP2007135894 A JP 2007135894A JP 2005334325 A JP2005334325 A JP 2005334325A JP 2005334325 A JP2005334325 A JP 2005334325A JP 2007135894 A JP2007135894 A JP 2007135894A
Authority
JP
Japan
Prior art keywords
blood flow
information
component vector
data
blood vessel
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
Application number
JP2005334325A
Other languages
Japanese (ja)
Inventor
Takashi Kosugi
隆司 小杉
Haruo Isoda
治夫 礒田
Shuhei Yamashita
修平 山下
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
R Tech kk
Hamamatsu University School of Medicine NUC
Original Assignee
R Tech kk
Hamamatsu University School of Medicine NUC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by R Tech kk, Hamamatsu University School of Medicine NUC filed Critical R Tech kk
Priority to JP2005334325A priority Critical patent/JP2007135894A/en
Publication of JP2007135894A publication Critical patent/JP2007135894A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

【課題】インビボでの血流動態解析を行い血管壁剪断応力及び圧力を求めることができる血流解析装置、血流シミュレーション方法、及びプログラムを提供する。
【解決手段】空間3次元と時間から成る4次元の速度3成分ベクトル情報をもつヒト血流データが入力されるデータ入力部10と、入力されたデータから4次元の速度3成分ベクトル情報を取得し、当該4次元の速度3成分ベクトル情報に基づいて血管壁剪断応力及び圧力を算出するデータ解析部20と、データ解析部での解析により得られた血管壁剪断応力及び圧力などの解析結果を表示する表示部50を有する。
【選択図】図1
A blood flow analysis apparatus, a blood flow simulation method, and a program capable of performing blood flow dynamics analysis in vivo to obtain blood vessel wall shear stress and pressure are provided.
A data input unit to which human blood flow data having four-dimensional velocity three-component vector information composed of three-dimensional space and time is input, and four-dimensional velocity three-component vector information is obtained from the input data. The data analysis unit 20 calculates the vascular wall shear stress and pressure based on the four-dimensional velocity three-component vector information, and the analysis results such as the vascular wall shear stress and pressure obtained by the analysis in the data analysis unit. It has the display part 50 to display.
[Selection] Figure 1

Description

本発明は、ヒト血流データをもとにした血流解析を行う血流解析装置、血流シミュレーション方法、及びプログラムに関する。   The present invention relates to a blood flow analysis device, a blood flow simulation method, and a program for performing blood flow analysis based on human blood flow data.

粥状動脈硬化病変、脳動脈瘤、解離性大動脈瘤などの血管病変には好発部位があり、これらの発生や進行には血管壁に加わる剪断応力が大きな役割を果たすと考えられる。個々の患者等の血管壁剪断応力を正確かつ容易に求めることができれば、病変の予後推定や治療方針決定に役立てることができる可能性が高く、また、将来の血管病変発生を推定して、その予防に役立てることができる可能性もある。   Vascular lesions such as atherosclerotic lesions, cerebral aneurysms, and dissecting aortic aneurysms have recurrent sites, and it is considered that shear stress applied to the vascular wall plays a major role in the occurrence and progression of these lesions. If it is possible to accurately and easily determine the vascular wall shear stress of individual patients, etc., there is a high possibility that it will be useful for estimating the prognosis of the lesion and determining the treatment strategy. There is a possibility that it can be useful for prevention.

正確な血管壁剪断応力及び圧力を求めるには、血管壁近傍の血流動態を正確に解析することが必要である。近年のヒト血流解析方法には、計算流体力学(CFD:Computational Fluid Dynamics)のインシリコン血流解析、レーザードップラー流速測定法(LDV:Laser Doppler Velocimetry)や粒子画像流速測定法(PIV:Particle Image Velocimetry)によるインビトロ(in vitro)血流解析、及び位相コントラストシネ磁気共鳴画像法(PC cine MRI:Phase Contrast cine Magnetic Resonance Imaging)によるインビボ(in vivo)血流解析などがある。   In order to obtain accurate vascular wall shear stress and pressure, it is necessary to accurately analyze the blood flow dynamics in the vicinity of the vascular wall. Recent human blood flow analysis methods include computational fluid dynamics (CFD) in-silicon blood flow analysis, laser doppler velocimetry (LDV) and particle image velocimetry (PIV). In vitro blood flow analysis by Velocimetry, and in vivo blood flow analysis by phase contrast cine magnetic resonance imaging (PC) MRI (Phase Contrast Cine Magnetic Resonance Imaging).

特開2005−40299号公報JP 2005-40299 A

上述した血流解析方法において、CFDによる血流解析やインビトロ血流解析では、生体の血流動態を正確に表現することが難しく、計算条件や実験条件などに依存して結果が異なることが知られている。また、個々の患者毎に血流動態を表現することも困難である。
このような事情から、インビボで生体から血流動態を直接測定できることが望ましい。しかしながら、現在、インビボ血流解析では、血流速度ベクトルは比較的簡便に得られるようになったものの、血管壁剪断応力及び圧力を簡便かつ正確に得ることができない。
In the blood flow analysis method described above, it is difficult to accurately express the blood flow dynamics of a living body in CFD blood flow analysis or in vitro blood flow analysis, and the results differ depending on calculation conditions and experimental conditions. It has been. It is also difficult to express blood flow dynamics for each individual patient.
Under such circumstances, it is desirable that blood flow dynamics can be directly measured from a living body in vivo. However, at present, in vivo blood flow analysis, blood flow velocity vectors can be obtained relatively easily, but vascular wall shear stress and pressure cannot be obtained easily and accurately.

本発明は、インビボでの血流動態解析を行い血管壁剪断応力及び圧力を求めることができる血流解析装置、血流シミュレーション方法、及びプログラムを提供することを目的とする。   It is an object of the present invention to provide a blood flow analysis device, a blood flow simulation method, and a program that can perform blood flow dynamic analysis in vivo and obtain blood vessel wall shear stress and pressure.

本発明の血流解析装置は、空間3次元と時間から成る4次元の速度3成分ベクトル情報をもつヒト血流データが入力される入力手段と、上記入力手段より入力されたデータから4次元の速度3成分ベクトル情報を取得し、当該4次元の速度3成分ベクトル情報に基づいて血管壁に加わる剪断応力、又は剪断応力及び圧力を算出する解析手段と、上記解析手段にて得られた解析結果を表示する表示手段とを備えることを特徴とする。
本発明の血流シミュレーション方法は、空間3次元と時間から成る4次元の速度3成分ベクトル情報をもつヒト血流データが入力される入力手段より入力されたデータから4次元の速度3成分ベクトル情報を含む血流情報を取得し、当該血流情報に基づいて血管壁に加わる剪断応力及び圧力を算出し、上記血流情報及び血管壁に加わる剪断応力及び圧力に係る情報の提供を受け、当該情報を用いて計算流体力学による血流シミュレーションを行うことを特徴とする。
本発明のプログラムは、空間3次元と時間から成る4次元の速度3成分ベクトル情報をもつヒト血流データが入力される入力手段より入力されたデータから取得した4次元の速度3成分ベクトル情報に基づいて、血管形状内の所定点における4次元の速度3成分ベクトル情報を補間する補間処理ステップと、上記情報取得ステップ及び上記補間処理ステップにて得られた4次元の速度3成分ベクトル情報を用い、血管形状内の所定領域の血管壁に加わる剪断応力、又は剪断応力及び圧力を算出する情報算出ステップと、上記各ステップにて得られる任意の情報を表示部に表示させる表示ステップとをコンピュータに実行させることを特徴とする。
The blood flow analysis apparatus according to the present invention includes input means for inputting human blood flow data having four-dimensional velocity three-component vector information composed of three-dimensional space and time, and four-dimensional data from the data input from the input means. Analyzing means for acquiring velocity three-component vector information and calculating shear stress applied to the blood vessel wall, or shear stress and pressure based on the four-dimensional velocity three-component vector information, and an analysis result obtained by the analyzing means And display means for displaying.
According to the blood flow simulation method of the present invention, four-dimensional velocity three-component vector information is obtained from data input from an input means to which human blood flow data having four-dimensional velocity three-component vector information consisting of three-dimensional space and time is input. Blood flow information including the blood flow information, calculating the shear stress and pressure applied to the blood vessel wall based on the blood flow information, receiving the blood flow information and information relating to the shear stress and pressure applied to the blood vessel wall, It is characterized by performing blood flow simulation by computational fluid dynamics using information.
The program according to the present invention converts four-dimensional velocity three-component vector information acquired from data input from input means to which human blood flow data having four-dimensional velocity three-component vector information consisting of space three dimensions and time is input. Based on the interpolation processing step for interpolating the four-dimensional velocity three-component vector information at a predetermined point in the blood vessel shape, and using the information acquisition step and the four-dimensional velocity three-component vector information obtained in the interpolation processing step. An information calculation step for calculating a shear stress applied to a blood vessel wall in a predetermined region in the blood vessel shape, or a shear stress and a pressure, and a display step for causing the display unit to display arbitrary information obtained in the above steps are displayed on the computer. It is made to perform.

本発明によれば、入力される4次元の速度ベクトルデータを収集可能な画像データに基づいて、3次元血流速度3成分ベクトル情報及び時間情報を有する4次元の速度3成分ベクトル情報を取得し、取得した4次元の速度3成分ベクトル情報を基に血管壁剪断応力及び圧力を算出するので、対象者を撮影して得られた画像データを用いたインビボでの血流解析により血管壁剪断応力及び圧力を正確かつ容易に求めることができる。   According to the present invention, based on image data that can collect input four-dimensional velocity vector data, three-dimensional blood velocity three-component vector information and four-dimensional velocity three-component vector information having time information are acquired. Since the vascular wall shear stress and pressure are calculated based on the acquired four-dimensional velocity three-component vector information, the vascular wall shear stress is obtained by in vivo blood flow analysis using image data obtained by photographing the subject. And the pressure can be determined accurately and easily.

以下、本発明の実施形態を図面に基づいて説明する。
図1は、本発明の一実施形態による血流解析装置の構成例を示すブロック図である。図1に示すように本実施形態による血流解析装置は、データ入力部10、データ解析部20、関心領域設定部30、結果出力部40、表示部50、及びデータ出力部60を有する。
Hereinafter, embodiments of the present invention will be described with reference to the drawings.
FIG. 1 is a block diagram illustrating a configuration example of a blood flow analysis device according to an embodiment of the present invention. As shown in FIG. 1, the blood flow analysis device according to the present embodiment includes a data input unit 10, a data analysis unit 20, a region of interest setting unit 30, a result output unit 40, a display unit 50, and a data output unit 60.

データ入力部10は、外部からデータを入力するためのものであり、例えば本実施形態では、3次元位相コントラストシネ磁気共鳴画像法(3D PC cine MRI)により得られた画像データが入力される。なお、この3次元位相コントラストシネ磁気共鳴画像法は、3次元空間におけるx、y、zの3軸すべてに速度エンコードを行い、時間軸を含めて4次元の速度ベクトルデータ(3次元空間速度3成分情報及び時間情報を有する速度ベクトルデータ)を収集可能な撮影法である。   The data input unit 10 is for inputting data from the outside. For example, in this embodiment, image data obtained by three-dimensional phase contrast cine magnetic resonance imaging (3D PC cine MRI) is input. This three-dimensional phase contrast cine magnetic resonance imaging method performs velocity encoding on all three axes x, y, and z in three-dimensional space, and includes four-dimensional velocity vector data (three-dimensional spatial velocity 3 This is an imaging method capable of collecting (speed vector data having component information and time information).

データ解析部20は、データ入力部10より供給されるデータの解析を行い、当該データに基づいて血流情報や血管壁剪断応力及び圧力などを算出する。データ解析部20は、血管情報作成部21、血流情報取得部22、補間処理部23、及び剪断応力及び圧力算出部24を有する。   The data analysis unit 20 analyzes data supplied from the data input unit 10 and calculates blood flow information, blood vessel wall shear stress, pressure, and the like based on the data. The data analysis unit 20 includes a blood vessel information creation unit 21, a blood flow information acquisition unit 22, an interpolation processing unit 23, and a shear stress and pressure calculation unit 24.

血管情報作成部21は、データ入力部10より供給されるデータに基づいて、血管形態(血管壁面形状、より詳しくは血管壁に相当するサーフェイスの形状)の情報を求めて血管形状ファイルを作成する。なお、この血管形状ファイルにおける血管形態の情報は、3次元空間での血管形状を経時的に得た情報、すなわち4次元的な血管形態の情報である。   Based on the data supplied from the data input unit 10, the blood vessel information creation unit 21 obtains information on the blood vessel form (blood vessel wall shape, more specifically, the shape of the surface corresponding to the blood vessel wall) and creates a blood vessel shape file. . The information on the blood vessel form in the blood vessel shape file is information obtained over time in a three-dimensional space, that is, information on a four-dimensional blood vessel form.

血流情報取得部22は、データ入力部10より供給されるデータから、経時的な3次元血流速度ベクトルデータ、すなわち4次元血流速度ベクトルデータ(4次元血流速度3成分情報)を含む血流情報(流速、流量等)を取得する。   The blood flow information acquisition unit 22 includes three-dimensional blood flow velocity vector data over time, that is, four-dimensional blood flow velocity vector data (four-dimensional blood flow velocity three-component information) from the data supplied from the data input unit 10. Acquire blood flow information (flow velocity, flow rate, etc.).

補間処理部23は、供給されるデータを基に作成された血管形状内における血流情報に係る補間処理を行う。具体的には、補間処理部23は、作成された血管形状内に所定の精度で格子点(計測補点)を生成し、その計測補点における血流情報を血流情報取得部22がデータより取得した血流情報を基に補間設定する。また、補間処理部23は、血管内の任意の指定位置における血流情報を、これら計測補点で補間設定された血流情報に基づいて算出する。なお、血管形状内における計測補点の生成精度は任意に指定可能である。   The interpolation processing unit 23 performs interpolation processing related to blood flow information in the blood vessel shape created based on the supplied data. Specifically, the interpolation processing unit 23 generates lattice points (measurement complements) with a predetermined accuracy in the created blood vessel shape, and the blood flow information acquisition unit 22 stores blood flow information at the measurement complements. Interpolation is set based on the acquired blood flow information. Further, the interpolation processing unit 23 calculates blood flow information at an arbitrary designated position in the blood vessel based on the blood flow information set by interpolation at these measurement complement points. In addition, the generation accuracy of the measurement complement point in the blood vessel shape can be arbitrarily designated.

剪断応力及び圧力算出部24は、血管情報作成部21により作成された血管形状ファイルと、血流情報取得部22及び補間処理部23により得られた4次元血流速度ベクトルデータとに基づいて、後述する関心領域設定部30等により設定された関心領域での血管壁剪断応力及び圧力を算出する。より詳細には、剪断応力及び圧力算出部24は、関心領域における4次元血流速度ベクトルデータから、所望の血管壁近傍における剪断速度(dv/dx)を計算する。ここで、dvは所望の血管壁に沿った血流の流速であり、dxは血管壁から流速dvの測定部位までの距離である。続いて、剪断応力及び圧力算出部24は、計算して得られた剪断速度と粘度とを乗算して血管壁剪断応力を算出する。   The shear stress and pressure calculation unit 24 is based on the blood vessel shape file created by the blood vessel information creation unit 21 and the four-dimensional blood flow velocity vector data obtained by the blood flow information acquisition unit 22 and the interpolation processing unit 23. The vascular wall shear stress and pressure in the region of interest set by the region-of-interest setting unit 30 described later are calculated. More specifically, the shear stress and pressure calculation unit 24 calculates a shear rate (dv / dx) in the vicinity of a desired blood vessel wall from the four-dimensional blood flow velocity vector data in the region of interest. Here, dv is the flow velocity of the blood flow along the desired blood vessel wall, and dx is the distance from the blood vessel wall to the measurement site of the flow velocity dv. Subsequently, the shear stress and pressure calculation unit 24 calculates the vascular wall shear stress by multiplying the calculated shear rate and viscosity.

関心領域設定部30は、血流情報や血管壁剪断応力及び圧力について算出やその結果出力を行う関心領域を設定する。関心領域設定部30は、3次元空間における複数の方向から任意の矩形がユーザーにより選択されることで、3次元空間上での関心領域を決定する。ユーザーによる一度の矩形選択では、3次元空間において視線方向へ無限遠方まで領域選択される。さらに、ユーザーにより任意の角度にオブジェクトを回転させることで関心領域の絞込みを行うことを可能にしている。なお、この関心領域設定部30により決定される領域が各種のデータ解析の対象範囲となるが、同時に存在可能な関心領域は1つに限定されず、複数存在しても良い。   The region-of-interest setting unit 30 sets a region of interest for calculation and output as a result of blood flow information, blood vessel wall shear stress, and pressure. The region-of-interest setting unit 30 determines a region of interest in the three-dimensional space by selecting an arbitrary rectangle from a plurality of directions in the three-dimensional space. When the user selects a rectangle once, the region is selected in the three-dimensional space up to infinity in the line-of-sight direction. Furthermore, the user can narrow down the region of interest by rotating the object at an arbitrary angle. The region determined by the region-of-interest setting unit 30 is a target range for various data analysis, but the number of regions of interest that can exist simultaneously is not limited to one, and a plurality of regions of interest may exist.

結果出力部40は、データ解析部20がデータの解析等を行って得られた血流情報や血管壁剪断応力及び圧力(以下、これらをまとめて「血流流れ場情報」とも称す。)を表示部50やデータ出力部60を介して出力する。表示部50は、例えばディスプレイなどの表示装置で構成され、結果出力部40からの指示に基づいて結果をユーザーが閲覧可能なように表示する。また、データ出力部60は、結果出力部40からの指示に基づいて結果を媒体(記録媒体や印刷媒体等)を介して出力したり、ファイル等のデータとして外部に出力したりする。   The result output unit 40 includes blood flow information, blood vessel wall shear stress and pressure (hereinafter, collectively referred to as “blood flow field information”) obtained by the data analysis unit 20 performing data analysis and the like. The data is output via the display unit 50 and the data output unit 60. The display unit 50 is configured by a display device such as a display, for example, and displays the result based on an instruction from the result output unit 40 so that the user can view the result. The data output unit 60 outputs a result via a medium (such as a recording medium or a print medium) based on an instruction from the result output unit 40, or outputs the result as data such as a file to the outside.

ここで、本実施形態において、表示部50を用いて表示される結果表示には、例えば以下の(a)〜(e)がある。
(a)入力されたデータより得られる血流速度ベクトルデータの経時的な3次元空間表示
(b)血流速度ベクトルデータに基づく、関心領域の血管断面における瞬時最高速度・最低速度・平均速度、経時的な最高速度・最低速度・平均速度の表示(そのグラフ化を含む。)
(c)血流情報に基づく、関心領域の血管断面における瞬時流量、経時的な流量、1心拍平均流量の表示(そのグラフ化表示を含む。)
(d)入力されたデータの解析結果として得られる関心領域の血管壁領域内における血管壁剪断応力及び圧力の経時的な3次元空間表示
(e)関心領域の血管壁領域内における血管壁剪断応力及び圧力の最大値・最小値・平均値・標準偏差値、及び経時的なそれらの表示(そのグラフ化表示を含む。)
Here, in the present embodiment, the result display displayed using the display unit 50 includes, for example, the following (a) to (e).
(A) Three-dimensional spatial display of blood flow velocity vector data obtained from input data over time (b) Instantaneous maximum velocity / minimum velocity / average velocity in the blood vessel cross section of the region of interest based on the blood flow velocity vector data Display of maximum speed, minimum speed, and average speed over time (including graphing)
(C) Display of instantaneous flow rate, temporal flow rate, and 1 heart rate average flow rate in the blood vessel cross section of the region of interest based on blood flow information (including its graphed display).
(D) Three-dimensional spatial display of vascular wall shear stress and pressure in the vascular wall region of the region of interest obtained as an analysis result of the input data over time (e) vascular wall shear stress in the vascular wall region of the region of interest And the maximum value, minimum value, average value, standard deviation value of pressure, and their display over time (including its graphed display)

次に、動作について説明する。
図2は、本実施形態による血流解析装置での解析処理動作例を示すフローチャートである。
Next, the operation will be described.
FIG. 2 is a flowchart showing an example of analysis processing operation in the blood flow analysis apparatus according to the present embodiment.

まず、医療用画像診断装置を用いて得られた対象者(患者等)に係る4次元的な画像データが、データ入力部10を介して読み込まれる(S11)。このとき、読み込まれた4次元的な画像データは、本実施形態における血流解析装置において処理可能なように必要に応じてデータ入力部10によりデータ変換処理が施される。画像データとしては、例えばDICOM(Digital Imaging and Communications in Medicine)形式の画像データが用いられる。   First, four-dimensional image data relating to a subject (patient or the like) obtained using the medical image diagnostic apparatus is read via the data input unit 10 (S11). At this time, the read four-dimensional image data is subjected to data conversion processing by the data input unit 10 as necessary so that it can be processed by the blood flow analysis apparatus according to the present embodiment. For example, image data in DICOM (Digital Imaging and Communications in Medicine) format is used as the image data.

本実施形態では、この4次元的な画像データは、上述したように3次元空間の3軸すべてに時間軸を含めた4次元の速度ベクトルデータを収集可能な3次元位相コントラストシネ磁気共鳴画像法(3D PC cine MRI)を用いて得られたMRI画像データとする。なお、本発明はこれに限定されるものではなく、4次元的な画像データは、CT(Computed Tomography)装置、MRI装置、回転血管撮影装置等の医療用画像診断装置を用いて得られる4次元の速度ベクトルデータが収集可能な任意の画像データを適用可能である。   In the present embodiment, the four-dimensional image data is obtained from the three-dimensional phase contrast cine magnetic resonance imaging method capable of collecting four-dimensional velocity vector data including the time axis in all three axes of the three-dimensional space as described above. It is set as the MRI image data obtained using (3D PC cine MRI). Note that the present invention is not limited to this, and the four-dimensional image data is obtained by using a medical image diagnostic apparatus such as a CT (Computed Tomography) apparatus, an MRI apparatus, and a rotating angiography apparatus. Any image data that can collect the velocity vector data can be applied.

次に、血管情報作成部21が、ステップS11において読み込まれたデータを用いて血管形状の解析を行い、3次元空間での経時的な血管形状を内容とする4次元的な血管形状ファイルを作成する(S12)。   Next, the blood vessel information creation unit 21 analyzes the blood vessel shape using the data read in step S11, and creates a four-dimensional blood vessel shape file containing the blood vessel shape over time in the three-dimensional space. (S12).

また、血管情報取得部22が、ステップS11において読み込まれたデータから4次元血流速度ベクトルデータを含む血流情報を取得する(S13)。続いて、補間処理部23が、ステップS12において作成された血管形状内に格子状の計算ポイント(計測補点)を生成して、各計算ポイントにおける血流情報をステップS13において取得した血流情報に基づいて算出することにより補間する(S14)。   In addition, the blood vessel information acquisition unit 22 acquires blood flow information including four-dimensional blood flow velocity vector data from the data read in step S11 (S13). Subsequently, the interpolation processing unit 23 generates grid-like calculation points (measurement complement points) in the blood vessel shape created in step S12, and the blood flow information obtained in step S13 is blood flow information at each calculation point. Is interpolated by calculation based on (S14).

次に、上述したステップS12〜S14により得られた血管形状及び血流情報(詳細には、それに含まれる血流速度ベクトルデータ)を参照し、剪断応力及び圧力算出部24が、予め関心領域設定部30により設定された関心領域での血管壁剪断応力及び圧力を、流速、その測定点の血管壁からの距離、及び流体である血液の粘性係数に基づいて算出する(S15)。
そして、上述したデータ解析部20の各機能部での処理により解析結果として得られた血流流れ場情報を、外部からの要求等に応じて表示部50に表示したり、データ出力部60を介して出力したりする(S16)。
Next, referring to the blood vessel shape and blood flow information (specifically, blood flow velocity vector data included therein) obtained in steps S12 to S14 described above, the shear stress and pressure calculation unit 24 sets the region of interest in advance. The vascular wall shear stress and pressure in the region of interest set by the unit 30 are calculated based on the flow velocity, the distance from the vascular wall at the measurement point, and the viscosity coefficient of blood as fluid (S15).
The blood flow field information obtained as an analysis result by the processing in each functional unit of the data analysis unit 20 described above is displayed on the display unit 50 according to a request from the outside, or the data output unit 60 is displayed. (S16).

図3は、上述した処理により得られる解析結果(流速分布)の表示例を示す図である。
本実施形態では、図3(a)、(b)に示すように、関心領域内の血管内部において指定された分解能で流速を算出した結果を、例えば色付のベクトル矢印で表示する。流速値は3次元ベクトル成分として流速の大小に応じて矢印が色付けされる。
FIG. 3 is a diagram illustrating a display example of an analysis result (flow velocity distribution) obtained by the above-described processing.
In this embodiment, as shown in FIGS. 3A and 3B, the result of calculating the flow velocity with the resolution specified inside the blood vessel in the region of interest is displayed, for example, with a colored vector arrow. The flow velocity value is colored as an arrow according to the magnitude of the flow velocity as a three-dimensional vector component.

なお、計測時間分のスライダを合わせて表示して、ユーザーが指定した任意時刻における情報を表示できるようにしても良いし、時刻を任意の速さで変化させたアニメーションとして情報を表示できるようにしても良い。また、3次元空間での各座標における最大流速が表示可能であり、血管形状がSTLファイルとして提供される場合には、ベクトル情報に重ねて表示する。また、すべての3次元オブジェクトは、ユーザーの指定に応じて回転や拡縮など描画変更可能なようにしても良い。   The slider for the measurement time may be displayed together so that the information at an arbitrary time specified by the user can be displayed, or the information can be displayed as an animation with the time changed at an arbitrary speed. May be. In addition, the maximum flow velocity at each coordinate in the three-dimensional space can be displayed, and when the blood vessel shape is provided as an STL file, it is displayed superimposed on the vector information. Also, all three-dimensional objects may be rendered changeable such as rotating or scaling according to the user's specification.

図4は、上述した処理により得られる解析結果(血流状態)の表示例を示す図である。
本実施形態では、図4(a)〜図4(d)に示すように、解析結果として得られる血流の状況を、流線図、流跡線図、流脈線図などにより表示することが可能である。図4(a)は、血流状態を流線図により示した例であり、流線はある瞬間の速度ベクトルをたどっていった線を示している。図4(b)は、血流状態を流跡線図により示した例であり、線は特定の流体の塊が経時的にたどった跡を示している。また、図4(c)、(d)は、血流状態を流脈線図により示した例であり、特定の一点を通過した流体のつながりを示している。なお、図4(c)は流脈を点を用いてポイント表示したものであり、図4(d)は流脈を線を用いてライン表示したものである。
FIG. 4 is a diagram illustrating a display example of an analysis result (blood flow state) obtained by the above-described processing.
In the present embodiment, as shown in FIGS. 4A to 4D, the state of blood flow obtained as an analysis result is displayed as a stream diagram, a trajectory diagram, a flow diagram, or the like. Is possible. FIG. 4A is an example showing a blood flow state by a flow diagram, and the stream line shows a line that follows a velocity vector at a certain moment. FIG. 4B is an example in which the blood flow state is shown by a trajectory diagram, and the line shows a trace of a specific fluid mass over time. FIGS. 4C and 4D are examples in which the blood flow state is shown by a flow diagram, and shows the connection of fluid that has passed a specific point. Note that FIG. 4C shows a flow point displayed using points, and FIG. 4D shows a flow line displayed using lines.

図5は、上述した処理により得られる解析結果(血管壁剪断応力)の表示例を示す図である。
本実施形態では、図5(a)、(b)に示すように、関心領域内の血管内部において指定された測定点における血管壁剪断応力を算出した結果を、例えばカラーマップ表示により表示する。
なお、図示しないが関心領域の血管断面における血流量の経時変化等もグラフ化などして表示可能であり、複数断面における血流量を1つに重ねて表示できるようにしても良い。
FIG. 5 is a diagram illustrating a display example of an analysis result (blood vessel wall shear stress) obtained by the above-described processing.
In the present embodiment, as shown in FIGS. 5A and 5B, the result of calculating the blood vessel wall shear stress at the measurement point designated inside the blood vessel in the region of interest is displayed by, for example, a color map display.
Although not shown, it is possible to display a change in blood flow over time in the blood vessel cross section of the region of interest in a graph or the like, and the blood flow volume in a plurality of cross sections may be displayed in an overlapping manner.

以上、説明したように本実施形態によれば、3次元位相コントラストシネ磁気共鳴画像法をはじめ、空間3次元と時間から成る4次元の速度3成分ベクトル情報をもつヒト血流データが入力されたデータなどにより得られた4次元の速度ベクトルデータを収集可能なデータに基づいて、関心領域における3次元血流速度3成分ベクトル情報及び時間情報を取得する。そして、得られた3次元血流速度3成分ベクトル情報及び時間情報を基に、所望の血管壁近傍における血管壁剪断応力及び圧力を算出し表示する。これにより、対象者(患者等)を実際に撮影して得られた画像データを用いて、インビボでの血流解析により血管壁剪断応力及び圧力を正確かつ容易に求めることができる。また、実際に対象者を撮影して得られた画像データを用いて血流解析を行うので、対象者毎の血流動態を容易に把握することができる。   As described above, according to the present embodiment, human blood flow data having four-dimensional velocity three-component vector information including three-dimensional space and time is input, including three-dimensional phase contrast cine magnetic resonance imaging. Three-dimensional blood flow velocity three-component vector information and time information in the region of interest are acquired based on data that can collect four-dimensional velocity vector data obtained from data or the like. Based on the obtained three-dimensional blood velocity three-component vector information and time information, the vascular wall shear stress and pressure near the desired vascular wall are calculated and displayed. Thereby, the blood vessel wall shear stress and pressure can be accurately and easily determined by in vivo blood flow analysis using image data obtained by actually photographing the subject (patient or the like). Moreover, since blood flow analysis is performed using image data obtained by actually photographing the subject, blood flow dynamics for each subject can be easily grasped.

なお、入力されたデータをデータ解析部20にて解析して得られる血流流れ場情報を解析結果として出力するだけではなく、当該血流流れ場情報を用いて計算流体力学(CFD)などによる血流シミュレーションを実行するようにしても良い。図6は、本発明の実施形態による血流解析装置の他の構成例を示すブロック図である。なお、この図6において、図1に示したブロック等と同一の機能を有するブロック等には同一の符号を付し、重複する説明は省略する。   In addition, the blood flow field information obtained by analyzing the input data by the data analysis unit 20 is not only output as an analysis result, but also by computational fluid dynamics (CFD) using the blood flow field information. A blood flow simulation may be executed. FIG. 6 is a block diagram illustrating another configuration example of the blood flow analysis device according to the embodiment of the present invention. In FIG. 6, blocks having the same functions as those shown in FIG. 1 are denoted by the same reference numerals, and redundant description is omitted.

図6に示す血流解析装置は、データ解析部20にて得られる血流流れ場情報を用いて血流シミュレーションを実行可能にした血流解析装置であり、データ入力部10、データ解析部20、関心領域設定部30、結果出力部40、表示部50、データ出力部60、シミュレーション実行部70、及び条件設定部80を有する。すなわち、図6に示す血流解析装置は、上述した図1に示した血流解析装置に、シミュレーション実行部70と条件設定部80とをさらに設けたものである。   The blood flow analysis device shown in FIG. 6 is a blood flow analysis device that enables blood flow simulation using the blood flow field information obtained by the data analysis unit 20, and includes a data input unit 10 and a data analysis unit 20. The region of interest setting unit 30, the result output unit 40, the display unit 50, the data output unit 60, the simulation execution unit 70, and the condition setting unit 80. That is, the blood flow analysis device shown in FIG. 6 is obtained by further providing a simulation execution unit 70 and a condition setting unit 80 in the blood flow analysis device shown in FIG.

シミュレーション実行部70は、入力されたデータをデータ解析部20にて解析して得られた血流流れ場情報が供給され、当該血流流れ場情報及び条件設定部80から入力されるシミュレーション条件を用いて血流シミュレーションを実行する。シミュレーション実行部70は、血管の時間的(経時的)なシミュレーションを行うものであり、例えば供給される血流流れ場情報を境界条件として用いて計算流体力学(CFD)による血流シミュレーションを行ったり、さらには構造力学(有限要素法:FEM)の処理を加えて計算流体力学(CFD)との連成問題を解き血流シミュレーションを行ったりする。
条件設定部80は、シミュレーション実行部70で実行するシミュレーションの条件を入力設定するものであり、境界条件や血管形状などを指定することが可能である。
The simulation execution unit 70 is supplied with blood flow field information obtained by analyzing the input data in the data analysis unit 20, and sets the simulation conditions input from the blood flow field information and the condition setting unit 80. To perform a blood flow simulation. The simulation execution unit 70 performs temporal (temporal) simulation of blood vessels, and performs, for example, blood flow simulation by computational fluid dynamics (CFD) using supplied blood flow field information as boundary conditions. In addition, structural mechanics (finite element method: FEM) processing is added to solve a coupled problem with computational fluid dynamics (CFD) to perform blood flow simulation.
The condition setting unit 80 is used to input and set conditions for the simulation executed by the simulation execution unit 70, and can specify boundary conditions, blood vessel shapes, and the like.

このように、対象者を実際に撮影して得られた画像データに基づいて得られた対象者毎の血流流れ場情報を用いて血流シミュレーションを行うことにより、生体内の血流動態を反映させた血流シミュレーションを行うことができる。   In this way, by performing blood flow simulation using blood flow field information for each subject obtained based on image data obtained by actually photographing the subject, blood flow dynamics in the living body can be determined. The reflected blood flow simulation can be performed.

また、入力されたデータを解析して得られた解析結果である血流流れ場情報を所定の記憶装置などに蓄積してデータベース化し、個々の血流流れ場情報を用いて予測を行ったり、血流シミュレーションを行ったりするようにしても良い。また、通信ネットワークを介して血流解析装置により血流解析を行ったり、解析結果を共有したりできるようにしても良い。   In addition, blood flow field information, which is the analysis result obtained by analyzing the input data, is stored in a predetermined storage device or the like to create a database, and prediction is performed using individual blood flow field information, A blood flow simulation may be performed. In addition, blood flow analysis may be performed by a blood flow analysis device or an analysis result may be shared via a communication network.

また、上述した実施形態では、データ解析部20にて、供給されたデータを基に得られた血管形状ファイルと、4次元血流速度ベクトルデータを含む血流情報とに基づいて、関心領域での血管壁剪断応力を算出するようにしているが、血管壁に加わる圧力を算出することも可能である。   In the above-described embodiment, the region of interest is determined based on the blood vessel shape file obtained based on the supplied data and the blood flow information including the four-dimensional blood flow velocity vector data by the data analysis unit 20. Although the blood vessel wall shear stress is calculated, the pressure applied to the blood vessel wall can also be calculated.

また、上述した実施形態では、4次元の速度ベクトルデータを収集可能なインビボの4次元的なデータを用いて血流解析等を行う場合について説明しているが、レーザードップラー流速測定法(LDV)や粒子画像流速測定法(PIV)により得られたインビトロの4次元的なデータを用いて同様に血流解析等を行うことも可能である。   In the above-described embodiment, a case where blood flow analysis or the like is performed using in vivo four-dimensional data capable of collecting four-dimensional velocity vector data has been described. However, laser Doppler velocimetry (LDV) Similarly, blood flow analysis and the like can be performed using in vitro four-dimensional data obtained by particle image velocimetry (PIV).

(本発明の他の実施形態)
なお、以上に説明した本実施形態の血流解析装置は、コンピュータのCPU又はMPU、RAM、ROMなどで構成されるものであり、RAMやROMに記憶されたプログラムが動作することによって実現できる。したがって、上記機能を果たすようにコンピュータを動作させるプログラムをコンピュータに読み込ませることによって実現できるものであり、そのプログラム及び当該プログラムをコンピュータが読み取り可能なように記録した、例えばCD−ROMのような記録媒体は、それぞれ本発明を構成する。上記プログラムを記録する記録媒体としては、CD−ROM以外に、フレキシブルディスク、ハードディスク、磁気テープ、光磁気ディスク、不揮発性メモリカード等を用いることができる。
また、コンピュータが供給されたプログラムを実行することにより上述の実施形態の機能が実現されるだけでなく、そのプログラムがコンピュータにおいて稼働しているOS(オペレーティングシステム)又は他のアプリケーションソフト等と共同して上述の実施形態の機能が実現される場合や、供給されたプログラムの処理のすべて又は一部がコンピュータの機能拡張ボードや機能拡張ユニットにより行われて上述の実施形態の機能が実現される場合も、かかるプログラムは本発明の実施形態に含まれる。
(Other embodiments of the present invention)
The blood flow analysis apparatus according to the present embodiment described above is configured by a CPU or MPU of a computer, a RAM, a ROM, and the like, and can be realized by operating a program stored in the RAM or ROM. Therefore, it can be realized by causing a computer to read a program for operating the computer so as to perform the above functions, and the program and the program are recorded so as to be readable by the computer, such as a CD-ROM. Each medium constitutes the present invention. As a recording medium for recording the program, a flexible disk, a hard disk, a magnetic tape, a magneto-optical disk, a nonvolatile memory card, and the like can be used in addition to the CD-ROM.
Moreover, not only the functions of the above-described embodiments are realized by executing a program supplied by a computer, but the program is also shared with an OS (operating system) or other application software running on the computer. When the functions of the above-described embodiment are realized, or when all or part of the processing of the supplied program is performed by a function expansion board or a function expansion unit of the computer, the functions of the above-described embodiment are realized. Such a program is included in the embodiment of the present invention.

例えば、本実施形態に示した血流解析装置は、図7に示すようなコンピュータ機能700を有し、そのCPU701により上述した実施形態での動作が実施される。
コンピュータ機能700は、図7に示すように、CPU701と、ROM702と、RAM703と、キーボード(KB)709のキーボードコントローラ(KBC)705と、表示部としてのCRTディスプレイ(CRT)710のCRTコントローラ(CRTC)706と、ハードディスク(HD)711及びフレキシブルディスク(FD)712のディスクコントローラ(DKC)707と、ネットワークインタフェースカード(NIC)708とが、システムバス704を介して互いに通信可能に接続された構成としている。
CPU701は、ROM702若しくはHD711に記憶されたソフトウェア、又はFD712より供給されるソフトウェアを実行することで、システムバス704に接続された各構成部を総括的に制御する。すなわち、CPU701は、上述したような動作を行うための処理プログラムを、ROM702、HD711、又はFD712から読み出して実行することで、上述した実施形態での動作を実現するための制御を行う。RAM703は、CPU701の主メモリ又はワークエリア等として機能する。
KBC705は、KB709や図示していないポインティングデバイス等からの指示入力を制御する。CRTC706は、CRT710の表示を制御する。DKC707は、ブートプログラム、種々のアプリケーション、ユーザファイル、ネットワーク管理プログラム、及び上述した実施形態における処理プログラム等を記憶するHD711及びFD712とのアクセスを制御する。NIC708はネットワーク713上の他の装置と双方向にデータをやりとりする。
For example, the blood flow analysis apparatus shown in the present embodiment has a computer function 700 as shown in FIG. 7, and the CPU 701 performs the operation in the above-described embodiment.
As shown in FIG. 7, the computer function 700 includes a CPU 701, a ROM 702, a RAM 703, a keyboard controller (KBC) 705 of a keyboard (KB) 709, and a CRT controller (CRTC) of a CRT display (CRT) 710 as a display unit. ) 706, a disk controller (DKC) 707 of a hard disk (HD) 711 and a flexible disk (FD) 712, and a network interface card (NIC) 708 are communicably connected to each other via a system bus 704. Yes.
The CPU 701 comprehensively controls each component connected to the system bus 704 by executing software stored in the ROM 702 or the HD 711 or software supplied from the FD 712. That is, the CPU 701 reads out a processing program for performing the above-described operation from the ROM 702, the HD 711, or the FD 712 and executes it, thereby performing control for realizing the operation in the above-described embodiment. The RAM 703 functions as a main memory or work area for the CPU 701.
The KBC 705 controls an instruction input from the KB 709 or a pointing device (not shown). The CRTC 706 controls display on the CRT 710. The DKC 707 controls access to the HD 711 and the FD 712 that store a boot program, various applications, user files, a network management program, a processing program in the above-described embodiment, and the like. The NIC 708 exchanges data with other devices on the network 713 in both directions.

なお、上記実施形態は、何れも本発明を実施するにあたっての具体化のほんの一例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。   The above-described embodiments are merely examples of implementation in carrying out the present invention, and the technical scope of the present invention should not be construed as being limited thereto. That is, the present invention can be implemented in various forms without departing from the technical idea or the main features thereof.

本発明の一実施形態による血流解析装置の構成例を示す図である。It is a figure which shows the structural example of the blood-flow analyzer by one Embodiment of this invention. 本実施形態による血流解析装置での解析処理動作の一例を示すフローチャートである。It is a flowchart which shows an example of the analysis processing operation | movement with the blood-flow analyzer by this embodiment. 解析結果(流速分布)の出力例を示す図である。It is a figure which shows the example of an output of an analysis result (flow velocity distribution). 解析結果(血流状態)の出力例を示す図である。It is a figure which shows the example of an output of an analysis result (blood flow state). 解析結果(血管壁剪断応力)の出力例を示す図である。It is a figure which shows the example of an output of an analysis result (blood vessel wall shear stress). 本発明の第2の実施形態による血流解析装置の他の構成例を示す図である。It is a figure which shows the other structural example of the blood-flow analyzer by the 2nd Embodiment of this invention. 本実施形態における画像処理装置を実現可能なコンピュータ機能を示すブロック図である。It is a block diagram which shows the computer function which can implement | achieve the image processing apparatus in this embodiment.

符号の説明Explanation of symbols

10 データ入力部
20 データ解析部
21 血管情報作成部
22 血流情報取得部
23 補間処理部
24 剪断応力及び圧力算出部
30 関心領域設定部
40 結果出力部
50 表示部
60 データ出力部
70 シミュレーション実行部
80 条件設定部
DESCRIPTION OF SYMBOLS 10 Data input part 20 Data analysis part 21 Blood vessel information creation part 22 Blood flow information acquisition part 23 Interpolation processing part 24 Shear stress and pressure calculation part 30 Area of interest setting part 40 Result output part 50 Display part 60 Data output part 70 Simulation execution part 80 Condition setting section

Claims (6)

空間3次元と時間から成る4次元の速度3成分ベクトル情報をもつヒト血流データが入力される入力手段と、
上記入力手段より入力されたデータから4次元の速度3成分ベクトル情報を取得し、当該4次元の速度3成分ベクトル情報に基づいて血管壁に加わる剪断応力、又は剪断応力及び圧力を算出する解析手段と、
上記解析手段にて得られた解析結果を表示する表示手段とを備えることを特徴とする血流解析装置。
Input means for inputting human blood flow data having four-dimensional velocity three-component vector information composed of three-dimensional space and time;
Analysis means for acquiring four-dimensional velocity three-component vector information from data input from the input means and calculating shear stress applied to the blood vessel wall or shear stress and pressure based on the four-dimensional velocity three-component vector information When,
A blood flow analysis apparatus comprising: display means for displaying an analysis result obtained by the analysis means.
上記解析手段は、
上記入力手段より入力されたデータから4次元の速度3成分ベクトル情報を取得する情報取得手段と、
上記情報取得手段にて取得した4次元の速度3成分ベクトル情報に基づいて、血管形状内の所定点における4次元の速度3成分ベクトル情報を補間する補間処理手段と、
上記情報取得手段及び上記補間処理手段により得られた4次元の速度3成分ベクトル情報を用い、血管形状内の所定領域の血管壁に加わる剪断応力、又は剪断応力及び圧力を算出する情報算出手段とを備えることを特徴とする請求項1記載の血流解析装置。
The analysis means is
Information acquisition means for acquiring four-dimensional velocity three-component vector information from data input from the input means;
Interpolation processing means for interpolating the four-dimensional velocity three-component vector information at a predetermined point in the blood vessel shape based on the four-dimensional velocity three-component vector information acquired by the information acquisition means;
Information calculation means for calculating shear stress, or shear stress and pressure applied to a blood vessel wall in a predetermined region in the blood vessel shape, using the four-dimensional velocity three-component vector information obtained by the information acquisition means and the interpolation processing means; The blood flow analysis apparatus according to claim 1, comprising:
上記解析手段は、上記入力手段より入力されたデータに基づいて血管形状の情報を作成する血管情報作成手段をさらに備えることを特徴とする請求項2記載の血流解析装置。   The blood flow analysis apparatus according to claim 2, wherein the analysis means further comprises blood vessel information creation means for creating blood vessel shape information based on the data input from the input means. 上記解析手段にて得られた解析結果を用いて計算流体力学による血流シミュレーションを行うシミュレーション実行手段をさらに備えることを特徴とする請求項1記載の血流解析装置。   The blood flow analysis apparatus according to claim 1, further comprising simulation execution means for performing a blood flow simulation by computational fluid dynamics using the analysis result obtained by the analysis means. 空間3次元と時間から成る4次元の速度3成分ベクトル情報をもつヒト血流データが入力される入力手段より入力されたデータから4次元の速度3成分ベクトル情報を含む血流情報を取得し、当該血流情報に基づいて血管壁に加わる剪断応力及び圧力を算出し、
上記血流情報及び血管壁に加わる剪断応力及び圧力に係る情報の提供を受け、当該情報を用いて計算流体力学による血流シミュレーションを行うことを特徴とする血流シミュレーション方法。
Obtaining blood flow information including four-dimensional velocity three-component vector information from data input from input means to which human blood flow data having four-dimensional velocity three-component vector information consisting of space three dimensions and time is input; Calculate the shear stress and pressure applied to the blood vessel wall based on the blood flow information,
A blood flow simulation method characterized by receiving the blood flow information and information related to shear stress and pressure applied to a blood vessel wall and performing blood flow simulation by computational fluid dynamics using the information.
空間3次元と時間から成る4次元の速度3成分ベクトル情報をもつヒト血流データが入力される入力手段より入力されたデータから4次元の速度3成分ベクトル情報を取得する情報取得ステップと、
上記情報取得ステップにて取得した4次元の速度3成分ベクトル情報に基づいて、血管形状内の所定点における4次元の速度3成分ベクトル情報を補間する補間処理ステップと、
上記情報取得ステップ及び上記補間処理ステップにて得られた4次元の速度3成分ベクトル情報を用い、血管形状内の所定領域の血管壁に加わる剪断応力、又は剪断応力及び圧力を算出する情報算出ステップと、
上記各ステップにて得られる任意の情報を表示部に表示させる表示ステップとをコンピュータに実行させるためのプログラム。
An information acquisition step of acquiring four-dimensional velocity three-component vector information from data input from input means for inputting human blood flow data having four-dimensional velocity three-component vector information consisting of space three dimensions and time;
An interpolation processing step for interpolating the four-dimensional velocity three-component vector information at a predetermined point in the blood vessel shape based on the four-dimensional velocity three-component vector information acquired in the information acquisition step;
An information calculation step for calculating shear stress or shear stress and pressure applied to the blood vessel wall in a predetermined region in the blood vessel shape using the four-dimensional velocity three-component vector information obtained in the information acquisition step and the interpolation processing step. When,
A program for causing a computer to execute a display step of displaying arbitrary information obtained in the above steps on a display unit.
JP2005334325A 2005-11-18 2005-11-18 Blood flow analysis device and simulation method based on human blood flow data Pending JP2007135894A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005334325A JP2007135894A (en) 2005-11-18 2005-11-18 Blood flow analysis device and simulation method based on human blood flow data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005334325A JP2007135894A (en) 2005-11-18 2005-11-18 Blood flow analysis device and simulation method based on human blood flow data

Publications (1)

Publication Number Publication Date
JP2007135894A true JP2007135894A (en) 2007-06-07

Family

ID=38199456

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005334325A Pending JP2007135894A (en) 2005-11-18 2005-11-18 Blood flow analysis device and simulation method based on human blood flow data

Country Status (1)

Country Link
JP (1) JP2007135894A (en)

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010207531A (en) * 2009-03-12 2010-09-24 Toshiba Corp Mri apparatus and data processor
WO2011030415A1 (en) * 2009-09-09 2011-03-17 株式会社ユネクス Device for examining vascular function
JP2011062358A (en) * 2009-09-17 2011-03-31 Tohoku Univ Ultrasonic diagnostic apparatus, blood flow visualization apparatus, and control program
JP2011062357A (en) * 2009-09-17 2011-03-31 Ge Medical Systems Global Technology Co Llc Ultrasonic diagnostic apparatus, blood flow visualization apparatus, and control program
JP2013504394A (en) * 2009-09-16 2013-02-07 モナシュ ユニバーシティ Imaging method
JP2014046109A (en) * 2012-09-03 2014-03-17 Toshiba Corp Medical image processor
JP2014236990A (en) * 2009-04-01 2014-12-18 コヴィディエン リミテッド パートナーシップ Microwave ablation system with user-controlled ablation size and method of use
JP2015013206A (en) * 2014-09-19 2015-01-22 株式会社東芝 Mri apparatus and data processing apparatus
JP2015044037A (en) * 2010-08-12 2015-03-12 ハートフロー, インコーポレイテッド Method and system for patient-specific blood flow modeling
JP2015526724A (en) * 2012-08-15 2015-09-10 アスペクト・エーアイ・リミテッド Measurement of fluid properties using MRI
JP2016140683A (en) * 2015-02-05 2016-08-08 東芝メディカルシステムズ株式会社 Medical image processing apparatus and magnetic resonance imaging apparatus
JP2016146995A (en) * 2015-02-13 2016-08-18 東芝メディカルシステムズ株式会社 Medical image processing device, and magnetic resonance imaging device
JP2017225807A (en) * 2016-05-02 2017-12-28 パイ メディカル イメージング ビー ヴイPie Medical Imaging B.V. Method and apparatus for automatically determining relationships between data
JP2018180707A (en) * 2017-04-06 2018-11-15 富士通株式会社 Flow line visualization device, flow line visualization method, and flow line visualization program
WO2019058657A1 (en) * 2017-09-19 2019-03-28 富士フイルム株式会社 Fluid analysis device, method for operating fluid analysis device, and fluid analysis program
US10354050B2 (en) 2009-03-17 2019-07-16 The Board Of Trustees Of Leland Stanford Junior University Image processing method for determining patient-specific cardiovascular information
US10537282B2 (en) 2016-04-12 2020-01-21 Canon Medical Systems Corporation Medical image processing apparatus, medical image diagnosis apparatus and recording medium
US11107587B2 (en) 2008-07-21 2021-08-31 The Board Of Trustees Of The Leland Stanford Junior University Method for tuning patient-specific cardiovascular simulations
JP2023027392A (en) * 2009-09-23 2023-03-01 ライトラボ・イメージング・インコーポレーテッド Lumen morphology and vascular resistance measurements data collection systems, apparatus and methods
US12082912B2 (en) 2009-09-23 2024-09-10 Lightlab Imaging, Inc. Lumen morphology and vascular resistance measurements data collection systems apparatus and methods

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07116143A (en) * 1993-10-25 1995-05-09 Hitachi Medical Corp Indicating method of measuring fluid in apparatus for mri
JPH09187436A (en) * 1996-01-05 1997-07-22 Hitachi Medical Corp Three-dimensional blood flow velocity image projection method in magnetic resonance diagnosing apparatus

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07116143A (en) * 1993-10-25 1995-05-09 Hitachi Medical Corp Indicating method of measuring fluid in apparatus for mri
JPH09187436A (en) * 1996-01-05 1997-07-22 Hitachi Medical Corp Three-dimensional blood flow velocity image projection method in magnetic resonance diagnosing apparatus

Cited By (89)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11107587B2 (en) 2008-07-21 2021-08-31 The Board Of Trustees Of The Leland Stanford Junior University Method for tuning patient-specific cardiovascular simulations
JP2010207531A (en) * 2009-03-12 2010-09-24 Toshiba Corp Mri apparatus and data processor
US12176094B2 (en) 2009-03-17 2024-12-24 The Board Of Trustees Of Leland Stanford Junior University Systems and methods for processing electronic images
US10354050B2 (en) 2009-03-17 2019-07-16 The Board Of Trustees Of Leland Stanford Junior University Image processing method for determining patient-specific cardiovascular information
JP2016127996A (en) * 2009-04-01 2016-07-14 コヴィディエン リミテッド パートナーシップ Microwave ablation system with user-controlled ablation size, and method of use
US9277969B2 (en) 2009-04-01 2016-03-08 Covidien Lp Microwave ablation system with user-controlled ablation size and method of use
US9867670B2 (en) 2009-04-01 2018-01-16 Covidien Lp Microwave ablation system and user-controlled ablation size and method of use
US10111718B2 (en) 2009-04-01 2018-10-30 Covidien Lp Microwave ablation system with user-controlled ablation size and method of use
JP2014236990A (en) * 2009-04-01 2014-12-18 コヴィディエン リミテッド パートナーシップ Microwave ablation system with user-controlled ablation size and method of use
US10499998B2 (en) 2009-04-01 2019-12-10 Covidien Lp Microwave ablation system with user-controlled ablation size and method of use
US8845542B2 (en) 2009-09-09 2014-09-30 Unex Corporation Blood vessel function inspecting apparatus
US9545241B2 (en) 2009-09-09 2017-01-17 Unex Corporation Blood vessel function inspecting apparatus
WO2011030415A1 (en) * 2009-09-09 2011-03-17 株式会社ユネクス Device for examining vascular function
JP2014193409A (en) * 2009-09-16 2014-10-09 Monash Univ Imaging method
JP2013504394A (en) * 2009-09-16 2013-02-07 モナシュ ユニバーシティ Imaging method
JP2011062357A (en) * 2009-09-17 2011-03-31 Ge Medical Systems Global Technology Co Llc Ultrasonic diagnostic apparatus, blood flow visualization apparatus, and control program
JP2011062358A (en) * 2009-09-17 2011-03-31 Tohoku Univ Ultrasonic diagnostic apparatus, blood flow visualization apparatus, and control program
US12082912B2 (en) 2009-09-23 2024-09-10 Lightlab Imaging, Inc. Lumen morphology and vascular resistance measurements data collection systems apparatus and methods
US12274532B2 (en) 2009-09-23 2025-04-15 Lightlab Imaging, Inc. Lumen morphology and vascular resistance measurements data collection systems apparatus and methods
US12138017B2 (en) 2009-09-23 2024-11-12 Lightlab Imaging, Inc. Lumen morphology and vascular resistance measurements data collection systems apparatus and methods
JP2023027392A (en) * 2009-09-23 2023-03-01 ライトラボ・イメージング・インコーポレーテッド Lumen morphology and vascular resistance measurements data collection systems, apparatus and methods
JP7494281B2 (en) 2009-09-23 2024-06-03 ライトラボ・イメージング・インコーポレーテッド Systems, devices and methods for lumen morphology and vascular resistance measurement data collection - Patents.com
US12121325B2 (en) 2009-09-23 2024-10-22 Lightlab Imaging, Inc. Lumen morphology and vascular resistance measurements data collection systems apparatus and methods
US12285240B2 (en) 2009-09-23 2025-04-29 Lightlab Imaging, Inc. Lumen morphology and vascular resistance measurements data collection systems apparatus and methods
US10080614B2 (en) 2010-08-12 2018-09-25 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US10321958B2 (en) 2010-08-12 2019-06-18 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US9585723B2 (en) 2010-08-12 2017-03-07 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US9697330B2 (en) 2010-08-12 2017-07-04 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US9706925B2 (en) 2010-08-12 2017-07-18 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US9743835B2 (en) 2010-08-12 2017-08-29 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US9801689B2 (en) 2010-08-12 2017-10-31 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US9839484B2 (en) 2010-08-12 2017-12-12 Heartflow, Inc. Method and system for image processing and patient-specific modeling of blood flow
US10080613B2 (en) 2010-08-12 2018-09-25 Heartflow, Inc. Systems and methods for determining and visualizing perfusion of myocardial muscle
US9855105B2 (en) 2010-08-12 2018-01-02 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US9861284B2 (en) 2010-08-12 2018-01-09 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US9271657B2 (en) 2010-08-12 2016-03-01 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US9888971B2 (en) 2010-08-12 2018-02-13 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US9167974B2 (en) 2010-08-12 2015-10-27 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US12016635B2 (en) 2010-08-12 2024-06-25 Heartflow, Inc. Method and system for image processing to determine blood flow
US12357387B2 (en) 2010-08-12 2025-07-15 Heartflow, Inc. Method and system for image processing to determine blood flow
US10052158B2 (en) 2010-08-12 2018-08-21 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US10092360B2 (en) 2010-08-12 2018-10-09 Heartflow, Inc. Method and system for image processing and patient-specific modeling of blood flow
US9268902B2 (en) 2010-08-12 2016-02-23 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US9235679B2 (en) 2010-08-12 2016-01-12 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
JP7590499B2 (en) 2010-08-12 2024-11-26 ハートフロー, インコーポレイテッド Method and system for patient-specific blood flow modeling - Patents.com
US10149723B2 (en) 2010-08-12 2018-12-11 Heartflow, Inc. Method and system for image processing and patient-specific modeling of blood flow
US10154883B2 (en) 2010-08-12 2018-12-18 Heartflow, Inc. Method and system for image processing and patient-specific modeling of blood flow
US10159529B2 (en) 2010-08-12 2018-12-25 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US10166077B2 (en) 2010-08-12 2019-01-01 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US10179030B2 (en) 2010-08-12 2019-01-15 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US9226672B2 (en) 2010-08-12 2016-01-05 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US9449147B2 (en) 2010-08-12 2016-09-20 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US10327847B2 (en) 2010-08-12 2019-06-25 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US9152757B2 (en) 2010-08-12 2015-10-06 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US10376317B2 (en) 2010-08-12 2019-08-13 Heartflow, Inc. Method and system for image processing and patient-specific modeling of blood flow
US10441361B2 (en) 2010-08-12 2019-10-15 Heartflow, Inc. Method and system for image processing and patient-specific modeling of blood flow
US10478252B2 (en) 2010-08-12 2019-11-19 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US10492866B2 (en) 2010-08-12 2019-12-03 Heartflow, Inc. Method and system for image processing to determine blood flow
US9149197B2 (en) 2010-08-12 2015-10-06 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US10531923B2 (en) 2010-08-12 2020-01-14 Heartflow, Inc. Method and system for image processing to determine blood flow
US11793575B2 (en) 2010-08-12 2023-10-24 Heartflow, Inc. Method and system for image processing to determine blood flow
US10682180B2 (en) 2010-08-12 2020-06-16 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US10702340B2 (en) 2010-08-12 2020-07-07 Heartflow, Inc. Image processing and patient-specific modeling of blood flow
US10702339B2 (en) 2010-08-12 2020-07-07 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
JP2015044036A (en) * 2010-08-12 2015-03-12 ハートフロー, インコーポレイテッド Method and system for patient-specific modeling of blood flow
US11033332B2 (en) 2010-08-12 2021-06-15 Heartflow, Inc. Method and system for image processing to determine blood flow
US11083524B2 (en) 2010-08-12 2021-08-10 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US11090118B2 (en) 2010-08-12 2021-08-17 Heartflow, Inc. Method and system for image processing and patient-specific modeling of blood flow
JP2015044037A (en) * 2010-08-12 2015-03-12 ハートフロー, インコーポレイテッド Method and system for patient-specific blood flow modeling
US11116575B2 (en) 2010-08-12 2021-09-14 Heartflow, Inc. Method and system for image processing to determine blood flow
US11135012B2 (en) 2010-08-12 2021-10-05 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US11154361B2 (en) 2010-08-12 2021-10-26 Heartflow, Inc. Method and system for image processing to determine blood flow
US11298187B2 (en) 2010-08-12 2022-04-12 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
US11583340B2 (en) 2010-08-12 2023-02-21 Heartflow, Inc. Method and system for image processing to determine blood flow
US12029494B2 (en) 2010-08-12 2024-07-09 Heartflow, Inc. Method and system for image processing to determine blood flow
JP2023112190A (en) * 2010-08-12 2023-08-10 ハートフロー, インコーポレイテッド Method and system for modeling patient-specific blood flow
JP2015526724A (en) * 2012-08-15 2015-09-10 アスペクト・エーアイ・リミテッド Measurement of fluid properties using MRI
JP2014046109A (en) * 2012-09-03 2014-03-17 Toshiba Corp Medical image processor
JP2015013206A (en) * 2014-09-19 2015-01-22 株式会社東芝 Mri apparatus and data processing apparatus
JP2016140683A (en) * 2015-02-05 2016-08-08 東芝メディカルシステムズ株式会社 Medical image processing apparatus and magnetic resonance imaging apparatus
US10120052B2 (en) 2015-02-05 2018-11-06 Toshiba Medical Systems Corporation Medical image processing apparatus and medical image processing method
US9962106B2 (en) 2015-02-13 2018-05-08 Toshiba Medical Systems Corporation Medical image processing apparatus and medical image processing method
JP2016146995A (en) * 2015-02-13 2016-08-18 東芝メディカルシステムズ株式会社 Medical image processing device, and magnetic resonance imaging device
US10537282B2 (en) 2016-04-12 2020-01-21 Canon Medical Systems Corporation Medical image processing apparatus, medical image diagnosis apparatus and recording medium
JP2017225807A (en) * 2016-05-02 2017-12-28 パイ メディカル イメージング ビー ヴイPie Medical Imaging B.V. Method and apparatus for automatically determining relationships between data
JP2018180707A (en) * 2017-04-06 2018-11-15 富士通株式会社 Flow line visualization device, flow line visualization method, and flow line visualization program
US11730384B2 (en) 2017-09-19 2023-08-22 Fujifilm Corporation Fluid analysis apparatus, method for operating fluid analysis apparatus, and fluid analysis program
JPWO2019058657A1 (en) * 2017-09-19 2020-09-24 富士フイルム株式会社 Fluid analysis device and operation method of fluid analysis device and fluid analysis program
WO2019058657A1 (en) * 2017-09-19 2019-03-28 富士フイルム株式会社 Fluid analysis device, method for operating fluid analysis device, and fluid analysis program

Similar Documents

Publication Publication Date Title
JP2007135894A (en) Blood flow analysis device and simulation method based on human blood flow data
JP7462816B2 (en) System and method for automated detection and visualization of disturbed blood flow using vector flow data - Patents.com
JP7232195B2 (en) Systems and methods for simultaneous visualization and quantification of intravascular wall shear stress
RU2491637C2 (en) Thickness gauge for measuring image objects
EP3600058B1 (en) System and method for concurrent visualization and quantification of blood flow using ultrasound vector flow imaging
JP6653673B2 (en) Blood flow analyzer, method, and program
JP2013192781A (en) Medical image display controller, method and program
KR102140383B1 (en) Method for microscope based augmented reality navigation, device and computer readable medium for performing the method
WO2018159708A1 (en) Blood flow analyzing device and method, and program
JP6934948B2 (en) Operation method of fluid analyzer and fluid analyzer and fluid analysis program
US20220268612A1 (en) Fluid analysis apparatus, fluid analysis method, and fluid analysis program
EP3607527B1 (en) Quantitative evaluation of time-varying data
JP7368500B2 (en) Fluid analysis equipment, methods and programs
KR101501516B1 (en) The method and apparatus for measuring a captured object using brightness information and magnified image of a captured image
JP6811872B2 (en) Fluid analysis device and operation method of fluid analysis device and fluid analysis program
JP5196714B2 (en) Image display device and image display program
JP6746639B2 (en) Image processing apparatus and method, and diagnostic support system
WO2020075571A1 (en) Fluid analysis device, method and program
KR101517753B1 (en) The method and apparatus for simultaneously displaying a plurality of parameters related to an ultrasound image

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20081107

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110906

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20120214