"体内の血管を流れる血流の流速情報を取得する取得手段と、前記取得した流速情報に基づいて、血管内の血液圧力を陽に含まない演算式により、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算する演算手段と、前記演算した、前記複数代表点におけるエネルギー損失の大きさを表す画像を生成する画像生成手段と、を含む血流可視化診断装置。請求項1記載の血流可視化診断装置であって、前記取得手段は、体表側から体内の血管を流れる血流の流速情報を検出する検出器から、当該検出した血流情報を取得し、前記演算手段は、前記検出器との関係で定められた座標系で記述した前記演算式を用いて、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算し、前記画像生成手段は、前記演算した、前記複数代表点におけるエネルギー損失の大きさを表す画像をデカルトの直交座標系に変換して生成する血流可視化診断装置。請求項2記載の血流可視化診断装置であって、前記検出器は、超音波信号を体表側から体内の血管に向けて放射し、体内から反射された超音波信号を受信して、当該受信した超音波信号により血流の流速情報を検出する超音波検出器であり、前記検出器との関係で定められた座標系は、超音波検出器の超音波信号の受信位置を原点とする極座標である血流可視化診断装置。コンピュータを、体内の血管を流れる血流の流速情報を取得する取得手段と、前記取得した流速情報に基づいて、血管内の血液圧力を陽に含まない演算式により、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算する演算手段と、前記演算した、前記複数代表点におけるエネルギー損失の大きさを表す画像を生成する画像生成手段と、として機能させるプログラム。このように従来の技術では、実際の患者の診断に必要な画像が必ずしも得られず、有用性に乏しいという問題点があった。本発明は上記実情に鑑みて為されたもので、実際の患者において測定された血流のエネルギー効率に関係する量を可視化して表示し、有用性を向上できる血流可視化診断装置、及びプログラムを提供することを、その目的の一つとする。上記従来例の問題点を解決する本発明は、血流可視化診断装置であって、体内の血管を流れる血流の流速情報を取得する取得手段と、前記取得した流速情報に基づいて、血管内の血液圧力を陽に含まない演算式により、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算する演算手段と、前記演算した、前記複数代表点におけるエネルギー損失の大きさを表す画像を生成する画像生成手段と、を含むこととしたものである。ここで前記取得手段は、体表側から体内の血管を流れる血流の流速情報を検出する検出器から、当該検出した血流情報を取得し、前記演算手段は、前記検出器との関係で定められた座標系で記述した前記演算式を用いて、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算し、前記画像生成手段は、前記演算した、前記複数代表点におけるエネルギー損失の大きさを表す画像をデカルトの直交座標系に変換して生成するものであってもよい。さらに前記検出器は、超音波信号を体表側から体内の血管に向けて放射し、体内から反射された超音波信号を受信して、当該受信した超音波信号により血流の流速情報を検出する超音波検出器であり、前記検出器との関係で定められた座標系は、超音波検出器の超音波信号の受信位置を原点とする極座標であることとしてもよい。また本発明の一態様に係るプログラムは、コンピュータを、体内の血管を流れる血流の流速情報を取得する取得手段と、前記取得した流速情報に基づいて、血管内の血液圧力を陽に含まない演算式により、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算する演算手段と、前記演算した、前記複数代表点におけるエネルギー損失の大きさを表す画像を生成する画像生成手段と、として機能させることとしたものである。本発明によると、診断に対する有用性を向上できる。本発明の実施の形態について図面を参照しながら説明する。本発明の実施の形態に係る血流可視化診断装置1は、図1に例示するように、体内の血管を流れる血流の流速情報を検出する検出装置2に接続される。またこの血流可視化診断装置1は、制御部11、記憶部12、操作部13、表示部14及び入力インタフェース15を含んで構成される。制御部11は、CPU等のプログラム制御デバイスであり、記憶部12に格納されたプログラムに従って動作する。本実施の形態では、この制御部11は、検出装置2から体内の血管を流れる血流の流速情報を取得し、当該取得した流速情報に基づいて、血管内の血液圧力を陽に含まない演算式により、血管内の注目領域内の複数代表点における血流のエネルギー損失を演算する。そして制御部11は、当該演算した、複数代表点におけるエネルギー損失の大きさを表す画像を生成して、表示出力する。この制御部11の具体的な動作については、後に詳しく述べる。記憶部12は、メモリデバイスやディスクデバイスを含んで構成される。この記憶部12は、制御部11によって実行されるプログラムが保持される。このプログラムは、DVD-ROMDigitalVersatileDiscReadOnlyMemory等のコンピュータ可読な記録媒体に格納されて提供され、この記憶部12に複写される。また、この記憶部12は、制御部11のワークメモリとしても動作する。操作部13は、マウスやキーボード等であり、利用者からの指示の入力を受け入れて、当該指示の内容を表す情報を、制御部11に出力する。表示部14は、液晶ディスプレイ等の表示デバイスであり、制御部11から入力される指示に従って、情報を表示出力する。入力インタフェース15は、検出装置2から、血流の流速に関する情報を受け入れて、制御部11に出力する。検出装置2は、例えば超音波診断装置であり、体表側から体内の血管に向けて超音波信号を放射し、体内から反射された超音波信号を受信する。具体的に超音波診断装置では、プローブ探触子が、予め定めた範囲注目領域内を走査しつつ周波数f0の超音波信号を各方向に放射する。ある方向に放射された超音波信号は、そのプローブからその方向にある血管内の赤血球で反射してプローブへ戻ってくるが、赤血球の移動速度に応じたドップラーシフトが生じるので、その周波数はf0からずれて、f0fdとなる。このfdドップラー偏移周波数によって赤血球の移動速度、すなわち血流の流速が検出できる。なお、超音波信号を二次元的に走査しつつ放射することで、二次元面内の各方向にある血流の流速の情報が得られる図2。この方法は、例えばカラードップラー法CDI等として広く知られているので、ここでの詳しい説明を省略する。検出装置2は、ほかには例えばMRIでもよい。通常のMRIが信号強度によって血管等を描出するのに対してマグニチュード画像、MRI位相画像は傾斜磁場の方向により血流速に比例した情報を表しており、空間上の任意の点の血流速度について方向とその速度の大きさとを検出できる。各方向の傾斜磁場の位相画像のシリーズを重ね合わせて、強度画像の血管形状の上にさらに重ね合わせることにより、血管内の血流の分布をベクトルで可視化することもできる。この方法はphasevelocitymappingとして広く知られており、ここでの詳しい説明は省略する。次に本実施の形態における制御部11の動作について説明する。本実施の形態の制御部11は、図3に例示するように、機能的には、情報取得部21と、演算部22と、画像生成部23と、表示制御部24とを含んで構成される。ここで情報取得部21は、検出装置2が出力する血流の流速情報を受け入れる。この情報は、既に例示したように、二次元面内の各方向座標x1,x2で表される各位置にある血流の流速の情報である。ここで流速の情報は、それぞれ、3次元的な各座標軸方向の成分としてベクトルu1,u2,u3で表される。演算部22は、受け入れた血流の情報を参照し、流速の情報が得られている二次元面内の各位置をそれぞれ代表点として、各代表点の血流の流速の情報に基づき、血管内の血液圧力を陽に含まない演算式、を用いて、当該複数の代表点における血流のエネルギー損失ΔEを演算する。なおηは粘性を表している。この1式は、所定の体積範囲Vを積分しなければならないが、各代表点における血流のエネルギー損失ΔEを演算する場合、被積分関数が略一定の局所的体積を積分したものとして、被積分関数そのものをエネルギー損失ΔEを表すものとして扱ってよい。すなわち、この演算部22は、1式の被積分関数微分は差分で置換えて演算すればよいを演算して、その値を出力する。この演算式がエネルギー損失を表すことについては後に説明する。画像生成部23は、演算部22にて得られた二次元面内の各代表点の血流のエネルギー損失ΔEに基づいて、二次元面内の対応する画素の画素値を定める。例えばエネルギー損失ΔEが0に近いほど青色、大きくなるにつれて緑、黄色、橙、赤と変化するように画素値を定める。なお、エネルギー損失ΔEそのものでなく、その対数であるlogΔEに基づいて画素値を定めてもよい。これにより画像生成部23は、検出装置2としての超音波診断装置が生成した二次元の画像に対応する、エネルギー損失の分布を表す画像を生成する。表示制御部24は、画像生成部23が生成した画像を、表示部14に出力して表示させる。ここで、演算部22が用いる1式により、エネルギー損失が演算できることについて説明する。血流は、非圧縮性の粘性流体であり、広く知られた流体力学の理論により、ナビエストークスの方程式が成立する。この3式で外力は0とした。またPは静水圧、ρは密度、ηは粘性、tは時刻を表す。また、速度uは既に述べたように、3次元的な各座標軸方向の成分を有するベクトルu1,u2,u3である。なお、血液では、密度ρは略1060kg/m3であり、粘性ηは30から50kg/m/sである。この3式ナビエストークスの式は、流体のエネルギーについては陽に表していない。流体のエネルギーは流体のもつ粘性による散逸によって減少する。この流体のエネルギーは、内部エネルギーをeとして、次の4式で表すことができる。ここで、Tは温度、Sはエントロピー、Vは体積、Pは圧力であり、これらにより4式のごとく内部エネルギーを表すことができることは、熱力学の基本法則による。さらにこの4式を微分することにより、を得る。人体内の血流に関しては、断熱的かつ、非圧縮性の流れと合理的に仮定できるので、dS=0、dV=0となり、従って、5式は簡単に、と表すことができる。ある体積内のエネルギーEの時間変化はしたがって、となる。ところで3式から、とできる。ここでさらに2式を用いて、が得られる。なお、である。これを血管等、注目領域全体の体積で積分すると、ガウスの定理と2式の連続の式とによって、が得られる。ここでnは、注目領域の表面に立つ法線ベクトルであり、Aは注目領域の表面面積を意味する。すなわち右辺側第1,第2項の積分は、面積分である。血流においては、図4に示すように血管壁において流速は0であり、注目領域Vに係る血流の流入路及び流出路図4のA1,A2では生理学的な血管では例えば末梢血管などで口径が一定で直線化された部分では十分な層流が得られていると仮定することは合理的であり、この場合当該流入面に対して速度uは垂直に交わり、これらの面で∂ui/∂ni=0と合理的に仮定できる。この仮定の下では、11式の右辺第2項は0になるので、11式は結局、と簡略化できる。さらに、血流は心臓による周期Tの脈動に影響されるから、周期境界条件を適用でき、12、13式より、を得る。この14式から、血流に関する合理的な仮定の下で、流体のエネルギー左辺が、右辺の式と等しくなることが示され、血流のエネルギー変化量ΔE、すなわち、エネルギー損失ΔEが、1式で表されることが理解される。この14式からはさらに次のことが理解される。この14式の左辺はさらに、血流の流入路と、流出路とにおける面積Aに係る被積分関数の値の和で表すことができる。つまり、流入路における圧力Pと流出路における圧力Pとの差が、エネルギー損失ΔEに影響することがわかる。いま、具体的に図5aからgの血管形状を考える。図5aは健康な血管形状の例である。Y字状に分岐した血管において、V字部の各血管から血液が流入するとすると、図5aでは、合流点においてスムーズに血流が進む結果、流入路側における圧力Pや、流出路における圧力Pは略同じとなる。従って、エネルギー損失は比較的小さくなる。ところが、図5bのように一方の血管からの合流部において狭窄が生じていると、流入路側における圧力Pは上昇するが、流出路での圧力Pが減少するためエネルギー損失が大きくなる。また図5cに示されるように、合流部において瘤ができていると、この瘤において血液の渦が生じて流入路側での圧力Pが減少するので、結果としてエネルギー損失が大きくなる。このことは局所的に見ても同じであり、狭窄や瘤の生じた位置でのエネルギー損失は大きくなる。また図5dは弁が開ききることのできる健康な血管を示す。ここでは14式の右辺側に示される場所に応じた速度変化の項に相当する値は比較的0に近い値となるどこでも流速がほとんど変わらない。従って、エネルギー損失は比較的小さい値となる。一方、図5fは弁の開きが小さくなり、血管が狭くなった状態血管の狭窄が生じている状態を示している。ここでは狭い開口部から血流が勢いよく吹き出されるjetflow。従ってこの場所での速度変化は大きくなり、この弁の場所でのエネルギー損失は比較的大きい値となる。さらに図5gは弁が閉じなくなり、血流が局所的に逆流する現象が生じた状態を表す。ここでは全体の流れに対して対流が生じることとなるので、速度変化が表れることとなり、当該部分でのエネルギー損失が比較的大きい値となる。このように、血流の流速そのものではなく、エネルギー損失ΔEを観測することによって、血管形状や弁の異常等によって起きる異常な血流が生じる場所を簡便に特定できることがわかる。すなわち、エネルギー損失が大きい場所では、血流を正常なものに近づけるために臓器への負担が強いられていることが理解され、患者への説明も容易になり、また医師も異常な血流のある箇所を容易に特定できる。なお、ここまでの説明では、検出装置2は超音波診断装置であるとして説明したが、本実施の形態はこれに限られない。すなわち、注目領域内における各部の血流の流速の情報が得られれば、1式による演算を行って、各部でのエネルギー損失が演算できるからである。具体的には検出装置2はMRIMagneticResonanceImagerであってもよい。また検出装置2として、超音波診断装置のようにプローブから二次元的に広がりのある領域を走査する装置を用いる場合は、これらの検出装置2にて実際に得られている血流の流速情報は、図6に例示するように、プローブと観測対象の赤血球Bとを結ぶ線分を動径方向とした極座標超音波信号の受信位置を原点とする極座標;r,θ,φによって表される。つまり情報取得部21が受け入れる、検出装置2が出力する血流の流速情報は、二次元面内の各方向座標x1,x2で表される各位置にある血流の流速の情報であるが、流速の情報自体は走査した面の座標に関わらずその成分u1,u2,u3が極座標系で表された値となる。そこで演算部22は、血管内の血液圧力を陽に含まない演算式である1式を極座標系で演算する。このように座標変換を行う方法については、広く知られているので、ここでの詳細な説明は省略する。すなわち演算部22は、流速の情報が得られている二次元面内の各位置デカルトの直交座標系で表される位置にある点をそれぞれ代表点として、各代表点における、極座標で表された血流の流速の情報に基づき、極座標系における1式を用いて、当該複数の代表点における血流のエネルギー損失ΔEスカラ量を演算する。画像生成部23は、演算部22にて得られた二次元面内の各代表点の血流のエネルギー損失ΔEに基づいて、二次元面内の対応する画素の画素値を定める。これにより画像生成部23は、演算された、複数代表点におけるエネルギー損失の大きさを表す画像をデカルトの直交座標系に変換して生成することとなる。表示制御部24は、画像生成部23が生成した画像を、表示部14に出力して表示させる。また検出装置2が他の装置である場合も、当該装置が各代表点で検出した血流の流速を表す速度ベクトルの座標系に変換した1式を用いてエネルギー損失ΔEを演算することとする。なお、MRIのphasevelocitymappingを用いて血流の流速情報を得る場合には、位相画像の方向により、直交座標系となることもある。この場合は直交座標系での1式を用いてエネルギー損失ΔEを演算する。本実施の形態は、以上の構成を備えてなり、次のように動作する。血流可視化診断装置1は、図7に示すように、検出装置2が出力する血流の流速情報を受け入れるS1。この流速情報は、二次元の面X,Y座標系内の各点Qi座標Xi,Yi、i=1,2,...Nにおける流速のベクトル値を関連付けたものである。そして血流可視化診断装置1は、当該受け入れた血流の情報を参照し、流速の情報が得られている二次元面内の各位置Qiをそれぞれ代表点として、各代表点Qiの血流の流速の情報に基づき、血管内の血液圧力を陽に含まない1式を用いて、当該複数の代表点における血流のエネルギー損失ΔEQiを演算するS2。ここでは既に述べたように、1式の被積分関数微分は差分で置換えて演算すればよいを演算して、その値を出力することとする。血流可視化診断装置1は、処理S2にて得た二次元面内の各代表点Qiの血流のエネルギー損失ΔEQiの対数に基づいて、二次元面内の対応する画素Qiの座標Xi,Yiにある画素の画素値を定めるS3。例えばエネルギー損失ΔEの対数logΔEが小さいほど青色、大きくなるにつれて緑、黄色、橙、赤と変化するように画素値を定める。そして血流可視化診断装置1は、当該生成した画像を出力表示するS4。本実施の形態によると、血流の異常に直接関係している、血管内の各部における血流のエネルギー損失の大きさを演算して表示するので、実際の患者の診断に必要な画像を得ることができ、有用性を向上できる。"