(57)【特許請求の範囲】
【請求項１】
流速の情報が得られている血管内の注目領域内の各代表点において、体内の血管を流れる血流の流速情報を取得する取得手段と、
前記取得した各代表点での流速情報に基づいて、血管内の血液圧力を陽に含まず、粘性散逸項を含む演算式を用いて、各代表点における血流のエネルギー損失を演算する演算手段と、
前記演算した、前記各代表点におけるエネルギー損失の大きさを表す画像を生成する画像生成手段と、
を含む血流可視化診断装置。
【請求項２】
請求項１記載の血流可視化診断装置であって、
前記演算式が、
【数１】
0005438744.tif 000016
（ただしηは粘性）にて表される血流可視化診断装置。
【請求項３】
請求項１から２のいずれか一項に記載の血流可視化診断装置であって、
前記取得手段は、体表側から体内の血管を流れる血流の流速情報を検出する検出器から、当該検出した血流情報を取得し、
前記演算手段は、前記検出器との関係で定められた座標系で記述した前記演算式を用いて、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算し、
前記画像生成手段は、前記演算した、前記複数代表点におけるエネルギー損失の大きさを表す画像をデカルトの直交座標系に変換して生成する血流可視化診断装置。
【請求項４】
請求項３記載の血流可視化診断装置であって、
前記検出器は、超音波信号を体表側から体内の血管に向けて放射し、体内から反射された超音波信号を受信して、当該受信した超音波信号により血流の流速情報を検出する超音波検出器であり、
前記検出器との関係で定められた座標系は、超音波検出器の超音波信号の受信位置を原点とする極座標である血流可視化診断装置。
【請求項５】
コンピュータを、
流速情報が得られている血管内の注目領域内の各代表点において、体内の血管を流れる血流の前記流速情報を取得する取得手段と、
前記取得した各代表点での流速情報に基づいて、血管内の血液圧力を陽に含まず、粘性散逸項を含む演算式を用いて、各代表点における血流のエネルギー損失を演算する演算手段と、
前記演算した、前記各代表点におけるエネルギー損失の大きさを表す画像を生成する画像生成手段と、
として機能させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、体内の血管内を流れる血流に関する情報を可視化する血流可視化診断装置、及びプログラムに関する。
【背景技術】
【０００２】
体内の血管を流れる血流を可視化する技術としては、コンピュータシミュレーションを用いる方法、超音波診断装置を用いる方法、ＭＲＩ（Magnetic Resonance Imaging）を用いる方法等が知られている。
【０００３】
このうちコンピュータシミュレーションを用いる方法では、血流に関する種々の情報を得ることはできるが、現実に計測される患者の血管ではないのが一般的であり、診断には利用できない。また超音波診断装置やＭＲＩを用いる方法では、患者の実際の状況を検査することができ、血管内の血液の流速の情報は得られる。しかしながら、血流速を得て可視化するだけでは、血流の状況は知ることができるが、診断上必要な情報が必ずしも得られない。
【０００４】
例えば、手術等で取り付けられた人工弁の動作状態や、血管にできている瘤の状態等は、血流速の情報のみでは把握できず、血流のエネルギー効率といった観点が必要となる。
【０００５】
特許文献１には、超音波信号により実際に測定した血管形状と血流情報とにより、コンピュータシミュレーションを行い、血流や血圧の分布を演算して表示する技術が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第４２６９６２３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
このように従来の技術では、実際の患者の診断に必要な画像が必ずしも得られず、有用性に乏しいという問題点があった。
【０００８】
本発明は上記実情に鑑みて為されたもので、実際の患者において測定された血流のエネルギー効率に関係する量を可視化して表示し、有用性を向上できる血流可視化診断装置、及びプログラムを提供することを、その目的の一つとする。
【課題を解決するための手段】
【０００９】
上記従来例の問題点を解決する本発明は、血流可視化診断装置であって、体内の血管を流れる血流の流速情報を取得する取得手段と、前記取得した流速情報に基づいて、血管内の血液圧力を陽に含まない演算式により、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算する演算手段と、前記演算した、前記複数代表点におけるエネルギー損失の大きさを表す画像を生成する画像生成手段と、を含むこととしたものである。
【００１０】
ここで前記取得手段は、体表側から体内の血管を流れる血流の流速情報を検出する検出器から、当該検出した血流情報を取得し、前記演算手段は、前記検出器との関係で定められた座標系で記述した前記演算式を用いて、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算し、前記画像生成手段は、前記演算した、前記複数代表点におけるエネルギー損失の大きさを表す画像をデカルトの直交座標系に変換して生成するものであってもよい。
【００１１】
さらに前記検出器は、超音波信号を体表側から体内の血管に向けて放射し、体内から反射された超音波信号を受信して、当該受信した超音波信号により血流の流速情報を検出する超音波検出器であり、前記検出器との関係で定められた座標系は、超音波検出器の超音波信号の受信位置を原点とする極座標であることとしてもよい。
【００１２】
また本発明の一態様に係るプログラムは、コンピュータを、体内の血管を流れる血流の流速情報を取得する取得手段と、前記取得した流速情報に基づいて、血管内の血液圧力を陽に含まない演算式により、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算する演算手段と、前記演算した、前記複数代表点におけるエネルギー損失の大きさを表す画像を生成する画像生成手段と、として機能させることとしたものである。
【発明の効果】
【００１３】
本発明によると、診断に対する有用性を向上できる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施の形態に係る血流可視化診断装置の構成例及び接続例を表すブロック図である。
【図２】本発明の実施の形態に係る血流可視化診断装置が受け入れる血流の流速情報の例を表す説明図である。
【図３】本発明の実施の形態に係る血流可視化診断装置の例を表す機能ブロック図である。
【図４】本発明の実施の形態に係る血流可視化診断装置での演算内容を説明するための概略説明図である。
【図５】本発明の実施の形態に係る血流可視化診断装置による診断対象となる血管の例を表す説明図である。
【図６】本発明の実施の形態に係る血流可視化診断装置に接続される検出装置が流速を測定する座標系の例を表す説明図である。
【図７】本発明の実施の形態に係る血流可視化診断装置の動作例を表すフローチャート図である。
【発明を実施するための形態】
【００１５】
本発明の実施の形態について図面を参照しながら説明する。本発明の実施の形態に係る血流可視化診断装置１は、図１に例示するように、体内の血管を流れる血流の流速情報を検出する検出装置２に接続される。またこの血流可視化診断装置１は、制御部１１、記憶部１２、操作部１３、表示部１４及び入力インタフェース１５を含んで構成される。
【００１６】
制御部１１は、ＣＰＵ等のプログラム制御デバイスであり、記憶部１２に格納されたプログラムに従って動作する。本実施の形態では、この制御部１１は、検出装置２から体内の血管を流れる血流の流速情報を取得し、当該取得した流速情報に基づいて、血管内の血液圧力を陽に含まない演算式により、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算する。そして制御部１１は、当該演算した、複数代表点におけるエネルギー損失の大きさを表す画像を生成して、表示出力する。この制御部１１の具体的な動作については、後に詳しく述べる。
【００１７】
記憶部１２は、メモリデバイスやディスクデバイスを含んで構成される。この記憶部１２は、制御部１１によって実行されるプログラムが保持される。このプログラムは、ＤＶＤ−ＲＯＭ（Digital Versatile Disc Read Only Memory）等のコンピュータ可読な記録媒体に格納されて提供され、この記憶部１２に複写される。また、この記憶部１２は、制御部１１のワークメモリとしても動作する。
【００１８】
操作部１３は、マウスやキーボード等であり、利用者からの指示の入力を受け入れて、当該指示の内容を表す情報を、制御部１１に出力する。表示部１４は、液晶ディスプレイ等の表示デバイスであり、制御部１１から入力される指示に従って、情報を表示出力する。入力インタフェース１５は、検出装置２から、血流の流速に関する情報を受け入れて、制御部１１に出力する。
【００１９】
検出装置２は、例えば超音波診断装置であり、体表側から体内の血管に向けて超音波信号を放射し、体内から反射された超音波信号を受信する。具体的に超音波診断装置では、プローブ（探触子）が、予め定めた範囲（注目領域）内を走査しつつ周波数ｆ0の超音波信号を各方向に放射する。ある方向に放射された超音波信号は、そのプローブからその方向にある血管内の赤血球で反射してプローブへ戻ってくるが、赤血球の移動速度に応じたドップラーシフトが生じるので、その周波数はｆ0からずれて、ｆ0＋ｆdとなる。このｆd（ドップラー偏移周波数）によって赤血球の移動速度、すなわち血流の流速が検出できる。なお、超音波信号を二次元的に走査しつつ放射することで、二次元面内の各方向にある血流の流速の情報が得られる（図２）。この方法は、例えばカラードップラー法（ＣＤＩ）等として広く知られているので、ここでの詳しい説明を省略する。
【００２０】
検出装置２は、ほかには例えばＭＲＩでもよい。通常のＭＲＩが信号強度によって血管等を描出するのに対して（マグニチュード画像）、ＭＲＩ位相画像は傾斜磁場の方向により血流速に比例した情報を表しており、空間上の任意の点の血流速度について方向とその速度の大きさとを検出できる。各方向の傾斜磁場の位相画像のシリーズを重ね合わせて、強度画像の血管形状の上にさらに重ね合わせることにより、血管内の血流の分布をベクトルで可視化することもできる。この方法はphase velocity mappingとして広く知られており、ここでの詳しい説明は省略する。
【００２１】
次に本実施の形態における制御部１１の動作について説明する。本実施の形態の制御部１１は、図３に例示するように、機能的には、情報取得部２１と、演算部２２と、画像生成部２３と、表示制御部２４とを含んで構成される。ここで情報取得部２１は、検出装置２が出力する血流の流速情報を受け入れる。この情報は、既に例示したように、二次元面内の各方向（座標ｘ1，ｘ2で表される各位置）にある血流の流速の情報である。ここで流速の情報は、それぞれ、３次元的な各座標軸方向の成分としてベクトル（ｕ1，ｕ2，ｕ3）で表される。
【００２２】
演算部２２は、受け入れた血流の情報を参照し、流速の情報が得られている二次元面内の各位置をそれぞれ代表点として、各代表点の血流の流速の情報に基づき、血管内の血液圧力を陽に含まない演算式、
【数１】
0005438744.tif 000002
を用いて、当該複数の代表点における血流のエネルギー損失ΔＥを演算する。なおηは粘性を表している。
【００２３】
この（１）式は、所定の体積範囲（Ｖ）を積分しなければならないが、各代表点における血流のエネルギー損失ΔＥを演算する場合、被積分関数が略一定の局所的体積を積分したものとして、被積分関数そのものをエネルギー損失ΔＥを表すものとして扱ってよい。すなわち、この演算部２２は、（１）式の被積分関数（微分は差分で置換えて演算すればよい）を演算して、その値を出力する。この演算式がエネルギー損失を表すことについては後に説明する。
【００２４】
画像生成部２３は、演算部２２にて得られた二次元面内の各代表点の血流のエネルギー損失ΔＥに基づいて、二次元面内の対応する画素の画素値を定める。例えばエネルギー損失ΔＥが０に近いほど青色、大きくなるにつれて緑、黄色、橙、赤と変化するように画素値を定める。なお、エネルギー損失ΔＥそのものでなく、その対数であるlogΔＥに基づいて画素値を定めてもよい。
【００２５】
これにより画像生成部２３は、検出装置２としての超音波診断装置が生成した二次元の画像に対応する、エネルギー損失の分布を表す画像を生成する。表示制御部２４は、画像生成部２３が生成した画像を、表示部１４に出力して表示させる。
【００２６】
ここで、演算部２２が用いる（１）式により、エネルギー損失が演算できることについて説明する。血流は、非圧縮性の粘性流体であり、広く知られた流体力学の理論により、ナビエ・ストークスの方程式が成立する。
【００２７】
【数２】
0005438744.tif 000003
【数３】
0005438744.tif 000004
【００２８】
この（３）式で外力は０とした。またＰは静水圧、ρは密度、ηは粘性、ｔは時刻を表す。また、速度ｕは既に述べたように、３次元的な各座標軸方向の成分を有するベクトル（ｕ1，ｕ2，ｕ3）である。なお、血液では、密度ρは略1060ｋｇ／ｍ3であり、粘性ηは3.0から5.0ｋｇ／ｍ／ｓである。
【００２９】
この（３）式（ナビエ・ストークスの式）は、流体のエネルギーについては陽に表していない。流体のエネルギーは流体のもつ粘性による散逸によって減少する。この流体のエネルギーは、内部エネルギーをｅとして、次の（４）式で表すことができる。
【００３０】
【数４】
0005438744.tif 000005
【００３１】
ここで、Ｔは温度、Ｓはエントロピー、Ｖは体積、Ｐは圧力であり、これらにより（４）式のごとく内部エネルギーを表すことができることは、熱力学の基本法則による。さらにこの（４）式を微分することにより、
【数５】
0005438744.tif 000006
を得る。人体内の血流に関しては、断熱的かつ、非圧縮性の流れと合理的に仮定できるので、ｄＳ＝０、ｄＶ＝0となり、従って、（５）式は簡単に、
【数６】
0005438744.tif 000007
と表すことができる。
【００３２】
ある体積内のエネルギーＥの時間変化はしたがって、
【数７】
0005438744.tif 000008
となる。ところで（３）式から、
【数８】
0005438744.tif 000009
とできる。ここでさらに（２）式を用いて、
【数９】
0005438744.tif 000010
が得られる。なお、
【数１０】
0005438744.tif 000011
である。
【００３３】
これを血管等、注目領域全体の体積で積分すると、ガウスの定理と（２）式の連続の式とによって、
【数１１】
0005438744.tif 000012
が得られる。ここでｎは、注目領域の表面に立つ法線ベクトルであり、Ａは注目領域の表面面積を意味する。すなわち右辺側第１，第２項の積分は、面積分である。
【００３４】
血流においては、図４に示すように血管壁において流速は０であり、注目領域（Ｖ）に係る血流の流入路及び流出路（図４のＡ１，Ａ２）では生理学的な血管では例えば末梢血管などで口径が一定で直線化された部分では十分な層流が得られていると仮定することは合理的であり、この場合当該流入面に対して速度ｕは垂直に交わり、これらの面で∂ｕi／∂ｎi＝０と合理的に仮定できる。この仮定の下では、（１１）式の右辺第２項は０になるので、（１１）式は結局、
【数１２】
0005438744.tif 000013
と簡略化できる。
【００３５】
さらに、血流は心臓による周期Ｔの脈動に影響されるから、周期境界条件
【数１３】
0005438744.tif 000014
を適用でき、（１２）、（１３）式より、
【数１４】
0005438744.tif 000015
を得る。この（１４）式から、血流に関する合理的な仮定の下で、流体のエネルギー（左辺）が、右辺の式と等しくなることが示され、血流のエネルギー変化量ΔＥ、すなわち、エネルギー損失ΔＥが、（１）式で表されることが理解される。
【００３６】
この（１４）式からはさらに次のことが理解される。この（１４）式の左辺はさらに、血流の流入路と、流出路とにおける面積Ａに係る被積分関数の値の和で表すことができる。つまり、流入路における圧力Ｐと流出路における圧力Ｐとの差が、エネルギー損失ΔＥに影響することがわかる。いま、具体的に図５（ａ）から（ｇ）の血管形状を考える。図５（ａ）は健康な血管形状の例である。Ｙ字状に分岐した血管において、Ｖ字部の各血管から血液が流入するとすると、図５（ａ）では、合流点においてスムーズに血流が進む結果、流入路側における圧力Ｐや、流出路における圧力Ｐは略同じとなる。従って、エネルギー損失は比較的小さくなる。
【００３７】
ところが、図５（ｂ）のように一方の血管からの合流部において狭窄が生じていると、流入路側における圧力Ｐは上昇するが、流出路での圧力Ｐが減少するためエネルギー損失が大きくなる。また図５（ｃ）に示されるように、合流部において瘤ができていると、この瘤において血液の渦が生じて流入路側での圧力Ｐが減少するので、結果としてエネルギー損失が大きくなる。このことは局所的に見ても同じであり、狭窄や瘤の生じた位置でのエネルギー損失は大きくなる。
【００３８】
また図５（ｄ）は弁が開ききることのできる健康な血管を示す。ここでは（１４）式の右辺側に示される場所に応じた速度変化の項に相当する値は比較的０に近い値となる（どこでも流速がほとんど変わらない）。従って、エネルギー損失は比較的小さい値となる。
【００３９】
一方、図５（ｆ）は弁の開きが小さくなり、血管が狭くなった状態（血管の狭窄が生じている状態）を示している。ここでは狭い開口部から血流が勢いよく吹き出される（jet flow）。従ってこの場所での速度変化は大きくなり、この弁の場所でのエネルギー損失は比較的大きい値となる。
【００４０】
さらに図５（ｇ）は弁が閉じなくなり、血流が局所的に逆流する現象が生じた状態を表す。ここでは全体の流れに対して対流が生じることとなるので、速度変化が表れることとなり、当該部分でのエネルギー損失が比較的大きい値となる。
【００４１】
このように、血流の流速そのものではなく、エネルギー損失ΔＥを観測することによって、血管形状や弁の異常等によって起きる異常な血流が生じる場所を簡便に特定できることがわかる。すなわち、エネルギー損失が大きい場所では、血流を正常なものに近づけるために臓器への負担が強いられていることが理解され、患者への説明も容易になり、また医師も異常な血流のある箇所を容易に特定できる。
【００４２】
なお、ここまでの説明では、検出装置２は超音波診断装置であるとして説明したが、本実施の形態はこれに限られない。すなわち、注目領域内における各部の血流の流速の情報が得られれば、（１）式による演算を行って、各部でのエネルギー損失が演算できるからである。具体的には検出装置２はＭＲＩ（Magnetic Resonance Imager）であってもよい。
【００４３】
また検出装置２として、超音波診断装置のようにプローブから二次元的に広がりのある領域を走査する装置を用いる場合は、これらの検出装置２にて実際に得られている血流の流速情報は、図６に例示するように、プローブと観測対象の赤血球Ｂとを結ぶ線分を動径方向とした極座標（超音波信号の受信位置を原点とする極座標；ｒ，θ，φ）によって表される。
【００４４】
つまり情報取得部２１が受け入れる、検出装置２が出力する血流の流速情報は、二次元面内の各方向（座標ｘ1，ｘ2で表される各位置）にある血流の流速の情報であるが、流速の情報自体は走査した面の座標に関わらずその成分（ｕ1，ｕ2，ｕ3）が極座標系で表された値となる。
【００４５】
そこで演算部２２は、血管内の血液圧力を陽に含まない演算式である（１）式を極座標系で演算する。このように座標変換を行う方法については、広く知られているので、ここでの詳細な説明は省略する。すなわち演算部２２は、流速の情報が得られている二次元面内の各位置（デカルトの直交座標系で表される位置）にある点をそれぞれ代表点として、各代表点における、極座標で表された血流の流速の情報に基づき、極座標系における（１）式を用いて、当該複数の代表点における血流のエネルギー損失ΔＥ（スカラ量）を演算する。
【００４６】
画像生成部２３は、演算部２２にて得られた二次元面内の各代表点の血流のエネルギー損失ΔＥに基づいて、二次元面内の対応する画素の画素値を定める。これにより画像生成部２３は、演算された、複数代表点におけるエネルギー損失の大きさを表す画像をデカルトの直交座標系に変換して生成することとなる。表示制御部２４は、画像生成部２３が生成した画像を、表示部１４に出力して表示させる。
【００４７】
また検出装置２が他の装置である場合も、当該装置が各代表点で検出した血流の流速を表す速度ベクトルの座標系に変換した（１）式を用いてエネルギー損失ΔＥを演算することとする。
なお、ＭＲＩのphase velocity mappingを用いて血流の流速情報を得る場合には、位相画像の方向により、直交座標系となることもある。この場合は直交座標系での（１）式を用いてエネルギー損失ΔＥを演算する。
【００４８】
本実施の形態は、以上の構成を備えてなり、次のように動作する。血流可視化診断装置１は、図７に示すように、検出装置２が出力する血流の流速情報を受け入れる（Ｓ１）。この流速情報は、二次元の面（Ｘ，Ｙ座標系）内の各点Ｑi（座標（Ｘi，Ｙi））、（ｉ＝１，２，…Ｎ）における流速のベクトル値を関連付けたものである。
【００４９】
そして血流可視化診断装置１は、当該受け入れた血流の情報を参照し、流速の情報が得られている二次元面内の各位置Ｑiをそれぞれ代表点として、各代表点Ｑiの血流の流速の情報に基づき、血管内の血液圧力を陽に含まない（１）式を用いて、当該複数の代表点における血流のエネルギー損失ΔＥ（Ｑi）を演算する（Ｓ２）。ここでは既に述べたように、（１）式の被積分関数（微分は差分で置換えて演算すればよい）を演算して、その値を出力することとする。
【００５０】
血流可視化診断装置１は、処理Ｓ２にて得た二次元面内の各代表点Ｑiの血流のエネルギー損失ΔＥ（Ｑi）の対数に基づいて、二次元面内の対応する画素（Ｑiの座標（Ｘi，Ｙi）にある画素）の画素値を定める（Ｓ３）。例えばエネルギー損失ΔＥの対数 log ΔＥが小さいほど青色、大きくなるにつれて緑、黄色、橙、赤と変化するように画素値を定める。そして血流可視化診断装置１は、当該生成した画像を出力表示する（Ｓ４）。
【００５１】
本実施の形態によると、血流の異常に直接関係している、血管内の各部における血流のエネルギー損失の大きさを演算して表示するので、実際の患者の診断に必要な画像を得ることができ、有用性を向上できる。
【符号の説明】
【００５２】
１  血流可視化診断装置、２  検出装置、１１  制御部、１２  記憶部、１３  操作部、１４  表示部、１５  入力インタフェース、２１  情報取得部、２２  演算部、２３  画像生成部、２４  表示制御部。
【図１】
0005438744.tif 000017
【図２】
0005438744.tif 000018
【図３】
0005438744.tif 000019
【図４】
0005438744.tif 000020
【図５】
0005438744.tif 000021
【図６】
0005438744.tif 000022
【図７】
0005438744.tif 000023