JP4559501B2 - 心機能表示装置、心機能表示方法およびそのプログラム - Google Patents

Description

本発明は、心臓の機能を表示する心機能表示装置およびそのプログラムに関するものである。

従来より、心臓の機能を表す表示方法としてブルズアイ（Bull’s Eye）表示法が用いられている。ブルズアイ表示法は、図２に示すように、心臓を楕円体モデルで近似して考えたとき、楕円の長軸を横切る方向に一定の間隔で切ったスライス面P1, P2,・・・,Pi,・・・，Pn上の心臓の機能の評価を表す画像Q1,Q2,・・・,Qi，・・・,Qn（図１３参照）を同心円状に配置して表示する手法である。ブルズアイ表示法では、楕円体モデルの一方の長軸頂点に近いスライス面P1の機能を同心円の中心に近い場所Q1に配置し、もう一方の長軸頂点に近いスライス面Pnの機能を同心円の外側Qnに配置する。

この表示方法は、主に機能画像を対象として用いられる。機能画像には、心筋シンチグラフィ（SPECT）によって直接撮影されるものや、CT/MRIから得られた解析結果から得られた機能を表した画像がある。

一方、アンギオグラフィ装置などにより、心臓の血管に造影剤などを注入することにより血管の状態を幾何学的に解析して、血管の狭窄箇所を発見することが可能になってきた。しかし、血管の位置と機能画像の位置とが対応していないため、アンギオグラフィ装置によって血管の狭窄箇所を発見したとしても機能画像でどの程度の支障があるかを確認することができなかった。あるいは、機能画像上に表れている異常が血管のどの辺りであるのかを把握することができなかった。

そこで、三次元の画像データから冠状動脈を抽出し、抽出した冠状動脈のパスを座標変換することによりブルズアイ表示された機能画像上に血管の形状を重ね合わせて表示することにより、機能画像で観察される状態と血管の位置を対応させて観察することができるようにしたものが提案されている（例えば、特許文献１など）。
特開2005-27999公報

しかしながら、特許文献１のように冠状動脈のパスを座標変換することによりブルズアイ表示された機能画像上に線を描画する方法では、実際に心臓の周囲を走行している血管の形態が表現されているものではないため、いずれの冠状動脈を観察しているのか判別するのが困難である。

そこで、本発明では、心臓の血管の形態が明瞭に表現できるようにして各位置における心臓の機能を正確に評価することが可能な心機能表示装置およびそのプログラムを提供することを目的とするものである。

本発明の心機能表示装置は、被写体を撮影して得た三次元医用画像を構成する第１のボクセルデータを用いて心臓の各位置における機能を表した機能画像を生成する機能画像生成手段と、
前記被写体を撮影して得た三次元医用画像を構成する前記第１のボクセルデータと同一または異なる第２のボクセルデータのうち心臓の心筋外壁に沿う血管を含む領域のボクセルデータを用いて、前記血管の形態を描画した形態画像を生成する形態画像生成手段と、
前記機能画像の心臓の各位置と前記形態画像の心臓の各位置とを対応させて、前記機能画像と前記形態画像とを重ねて表示する重ね表示手段とを備えたことを特徴とするものである。

本発明のプログラムは、コンピュータを、
被写体を撮影して得た三次元医用画像を構成する第１のボクセルデータを用いて心臓の各位置における機能を表した機能画像を生成する機能画像生成手段と、
前記被写体を撮影して得た三次元医用画像を構成する前記第１のボクセルデータと同一または異なる第２のボクセルデータのうち心臓の心筋外壁に沿う血管を含む領域のボクセルデータを用いて、前記血管の形態を描画した形態画像を生成する形態画像生成手段と、
前記機能画像の心臓の各位置と前記形態画像の心臓の各位置とを対応させて、前記機能画像と前記形態画像とを重ねて表示する重ね表示手段として機能させることを特徴とするものである。

「心臓の各位置における機能」とは、心臓の各位置の働きをいい、具体的には、例えば、心筋の各位置における厚さ・拡張期や収縮期の厚さの変化から心臓の各位置における働きを推定することができるものであり、「機能画像」は、各心臓の各位置の機能を視認できるように表した画像である。具体的には、例えば、機能に応じて色を変えたり数値で各心臓の位置の機能をあらわすことができる。

「形態」とは、表れている形やありさまをいい、「血管の形態」とは、血管の形や血管の走行の状態などを表したものをいう。

また、前記機能画像は、前記心臓の心基部から心尖部に向けて伸びる長軸を横切る方向に切った複数のスライス面における心臓の機能を表す画像を同心円状に配置した画像が望ましい。

また、前記機能画像は、前記心臓の心基部に近いスライス面の機能を表す画像を前記同心円上に配置した画像の中心に近い位置に配置し、前記スライス面の位置が前記心臓の心尖部に近づくほど該スライス面の機能を表す画像を前記中心より遠い位置に配置した画像が望ましい。

さらに、前記形態画像は、前記心臓の心筋外壁に沿う血管を含む該心筋外壁から所定の距離内に存在するボクセルデータに対してＭＩＰ処理を施して、該距離内に存在する血管を投影した画像であってもよい。

さらにまた、前記形態画像上に描画された血管の前記心臓上の範囲が、前記機能画像に表された機能の前記心臓上の範囲より広い範囲を表すものであってもよい。

本発明によれば、心臓の各位置における機能を表した機能画像と、血管の形態を描画した形態画像とを、機能画像の心臓の各位置と形態画像の心臓の各位置とが対応するように重ねて表示することにより、各心臓の機能と心臓の解剖学的な位置を対応させて観察することができる。

また、心臓の心基部から心尖部に向けて伸びる長軸を横切る方向に切った複数のスライス面における心臓の機能を表す画像を同心円状に配置し、心基部の画像を同心円の中心付近に配置し、心臓の心尖部に近づくほど同心円の外側に配置することにより、重ねて表示した形態画像の冠状動脈が中心から外側に向かって分岐していくので、心臓の解剖学的な位置が把握しやすい画像を生成することができる。

さらに、心臓の心筋外壁に沿う血管を含む該心筋外壁から所定の距離内に存在するボクセルデータに対してＭＩＰ処理（最大値投影法）を施して形態画像を生成することにより、血管の形態を詳細に表した画像を生成することができ、心臓の解剖学的な位置がさらに把握しやすくなる。

さらにまた、形態画像に描かれた心臓の血管の範囲を機能画像に表された機能の範囲より広くすることにより、機能画像と形態画像を重ねて表示しても機能画像の表示範囲に限定されることなく心臓の解剖学的な位置を確認できる部分があるので、心臓の位置と機能との対応が把握しやすい。

以下、図面を参照して本発明の心機能表示装置の実施の形態を詳細に説明する。図１は本発明の心機能表示装置の概略構成図である。なお、図１のような心機能表示装置１の構成は、補助記憶装置に読み込まれた心機能表示処理プログラムをコンピュータ上で実行することにより実現される。このとき、この心機能表示処理プログラムは、ＣＤ−ＲＯＭ等の記憶媒体に記憶され、もしくはインターネット等のネットワークを介して配布され、コンピュータにインストールされる。

本発明の心機能表示装置１は、被写体を撮影して得た三次元医用画像を構成する第１のボクセルデータ１００を記憶する第１の記憶手段１０と、第１のボクセルデータ１００を用いて心臓の各位置における機能を表した機能画像２００を生成する機能画像生成手段２０と、前記被写体を撮影して得た三次元医用画像を構成する第２のボクセルデータ１１０を記憶する第２の記憶手段３０と、被写体を撮影して得た三次元医用画像を構成する第２のボクセルデータ１１０のうち心臓の心筋外壁に沿う血管を含む領域のボクセルデータを用いて、前記血管の形態を描画した形態画像を生成する形態画像生成手段４０と、第１のボクセルデータ１００と第２のボクセルデータ１１０を位置合せする位置合せ手段５０と、形態画像２１０と機能画像２００とを重ねて表示装置２上に表示する重ね表示手段６０とを備える。

第１、第２のボクセルデータ１００、１１０は、ＣＴ（Computed Tomography）装置やＭＲＩ（Magnetic Resonance Imaging）装置などにより撮影されて得られた３次元画像である。第１のボクセルデータ１００は、機能画像２００を生成するために用いられ、第２のボクセルデータ１１０は、形態画像２１０を生成するために用いられる。

第１、第２の記憶手段１０、３０は、大容量の記憶装置であり、ハードディスクや画像サーバなどである。被写体をＣＴ装置やＭＲＩ装置などを用いて撮影した第１、第２のボクセルデータ１００、１１０が記憶される。

機能画像２００は、心臓の動きや心室の内径や心筋の厚さなどに応じて、心臓の機能の評価を表す評価値の分布を表す画像である。具体的には、あるフェーズにおける心室の径を表した心室内径画像、拡張期における心室の径を表した拡張末期心室内径の画像、収縮期における心室の径を表した収縮末期心室内径の画像、区画された領域毎の駆出率を表した局所駆出率の画像、あるフェーズでの心筋の厚みを表した壁厚の画像、拡張期における心筋の厚さを表した拡張末期壁厚の画像、収縮期における心筋の厚さを表した収縮末期壁厚の画像、拡張期の心筋の厚さから収縮期の心筋の厚さの差を表した壁厚変化量の画像、拡張期の心筋の厚さをA、収縮期の心筋の厚さをBとした時、B-A/Aで計算された値を表した壁厚増加率の画像、拡張期の心室の径から収縮期の心室の径の差を表した壁運動量の画像、心筋シンチグラフィの画像などがある。

心筋シンチグラフィを表す機能画像２００は、心筋に集まる薬を被検者の腕から注射して、その薬の分布を体の外から測定したデータを画像として表したもので、薬を変えることで、心筋の血流状態、心筋組織の代謝、神経の働きなどを表すことができる。

また、機能画像２００で心臓の動きの評価を表す場合には、動いている心臓を複数フェーズ撮影して画像間の違いから心臓の機能を評価する値を得て機能画像２００を作成する。このように複数フェーズ撮影する場合には、被曝することがないＭＲＩ装置で複数フェーズ撮影した画像から機能画像２００を生成するのが望ましい。

一方、形態画像２１０は心臓の心筋外壁に沿う血管などの形態を表す画像であるが、各臓器の構造が明瞭に表れるＣＴ装置で撮影した三次元画像から生成するのが好ましい。

また、形態画像２１０の作成に用いる第２のボクセルデータには、形態を詳細に表すために、ＣＴ装置やＭＲＩ装置などで細かいピッチ（例えば、１〜３ｍｍ間隔）で撮影した断層画像から生成したボクセルデータを用いるのが好ましい。一方、機能画像２００の作成するために心臓の動きに応じて複数のフェーズで撮影する必要がある場合は、細かいピッチで撮影を行ったのでは、複数フェーズを撮影することができない場合があるので、第１のボクセルデータは、ＭＲＩ装置などで広いピッチ（例えば、５〜１０ｍｍ間隔）で撮影した断層画像から生成したボクセルデータを用いることが多い。

機能画像生成手段２０は、第１の記憶手段１０に記憶されている第１のボクセルデータ１００より機能画像２００を生成する。ここでは、心筋の壁厚を表した心機能画像２００をブルズアイ表示法で生成する場合について具体的に説明する。

まず、第１のボクセルデータ１００から、図２に示すような楕円形状Ｅで表面モデルを近似して、心臓の長軸Ａ_１、短軸Ａ_２を決定する。この長軸Ａ_１は、心臓の心基部から心尖部に向けて伸び、さらに、心室領域の中心部を貫くように決定する。また、短軸Ａ_２は長軸Ａ_１に直交する。

この長軸Ａ_１を横切る方向（つまり短軸Ａ_２の方向）に一定の間隔で心臓を切ったスライス面P1,P2,…,Pi,…,Pn-1,Pn上の心臓の断面画像を第１のボクセルデータ１００から生成して心筋の壁厚を求める。例えば、図４に示すようなスライス面Pi上の断面から心内膜の輪郭と心外膜の輪郭を抽出し、各スライス面Pi上の心室領域内の長軸が通る点Cから放射線状に広がる線l1,l2,l3,…,lj,…,lｍ上の心内膜の輪郭と心外膜の輪郭の間の距離を壁厚ｄとして求める。

ブルズアイ表示する際、心筋の壁厚ｄに応じて色を変えて表示し、心臓の心基部に近いスライス面P1,P2,・・・の壁厚ｄを中心の近くの同心円上に表示し、心臓の心尖部に近いスライス面Pn,Pn-1,・・・の壁厚ｄを中心より遠い外周側の同心円上に表示する。また、スライス面上の各線l1,l2,l3,…,lj,…,lｍ上の壁厚ｄを基準線Ｌと各線l1,l2,l3,…,lj,…,lｍ間の角度θに対応させて、スライス面Piの壁厚ｄを同心円Qi上に表示した短軸画像を生成する。図３に示すように、このようにして作成された同心円Q1,Q2,…,Qi,…,Qn-1,Qnに表示された短軸画像を重ねてブルズアイ表示した機能画像２００を生成する。図５に心筋の壁厚ｄをブルズアイ表示した一例を示す。図５の各領域は壁厚ｄに応じて分けられた状態を示し、各領域は壁厚ｄに応じて色を変えるようにしたものが好ましい。あるいは、領域内にマウスなどを持っていくと数値で厚さが示されるようにしてもよい。

形態画像生成手段４０は、第２の記憶手段３０に記憶されている第２のボクセルデータ１１０から形態画像２１０を生成する。形態画２１０では、心筋外壁に沿って走行する冠状動脈の情報が重要であることから、冠状動脈を含むような心臓の表面モデルを作成する。まず、図６(a)のような形状をした心臓より、心筋外壁に沿って走行する冠状動脈を抽出し（図６(b)参照）、その芯線間を補間したりスプラインなどの関数当てはめることによって連続曲面を推定して表面モデルを生成する。図６(c)に生成した表面モデルの一例を示す。さらに、この表面モデルを基準に、機能画像生成手段２０と同様に、図２に示すような楕円形状Ｅで表面モデルを近似して、心臓の長軸Ａ_１、短軸Ａ_２の方向を決定する。

図７(a)に示すように、長軸Ａ１を横切るスライス面Piを生成する。血管の画素値は大きい値を持つので、血管を含むボクセルデータに対してMIP処理を施すことで血管を投影することができる。そこで、図７(a) ,(b)（図７(b)は図７(a)の一部を拡大した図）に示すように、各スライス画像Piの中心Ｃ（スライス面が長軸Ａ_１と交差する点）から各方向に放射状に伸びた線が表面モデルとの交差する一定の距離Ｄ内に存在するボクセルデータを用いてMIP処理を行い、MIP処理で各線上l1,l2,l3,…,lj,…,lｍの距離ｄ内を探索して得られた最大となる画素値（以下、MIP値という）を投影してブルズアイ表示をする。この距離Ｄは血管構造が得られるように、抽出時に得られた血管Ｂの血管径をもとに決定する。スライス面P1,P2,…,Pi,…,Pn-1,Pnの間隔はなるべく細かくしてMIP処理を行なうようにすることで、心筋外壁に沿う冠状動脈などの血管の走行の状態が詳細に観察できる。また、心基部に近いスライス面P1,P2,・・・のMIP値を中心の近くの同心円上に表示し、心臓の心尖部に近いスライス面Pn,Pn-1,・・・のMIP値を中心より遠い外周側の同心円上に表示すると、図８に示すように、心筋外壁に沿う冠状動脈の分岐が同心円の中心から外側に向かって分岐して行くので、心筋外壁に沿って心基部から心尖部に向かって走行する管状動脈が観察しやすくなる。また、MIP処理を用いることで、血管の形態を詳細に表すことができる。

位置合せ手段５０は、第１のボクセルデータ１００の心臓と第２のボクセルデータ１１０の心臓を位置合わせする。例えば、第１のボクセルデータ１００がＭＲＩ装置で心臓を撮影した画像であり、第２のボクセルデータ１１０がＣＴ装置で心臓を撮影した画像である場合には、同じ被写体を撮影した場合であっても、全く同じ位置になるように撮影することは難しい。そこで、機能画像２００の対応する心臓の位置に心臓の周囲に存在する血管の形態を投影するためには、第１のボクセルデータ１００の心臓の形状と第２のボクセルデータ１１０の心臓の形状とを位置合わせを行なう。

この位置合わせはレジストレーションと呼ばれる種々の手法を用いることができる。ここで、第１のボクセルデータ１００の画像Ａの点ｘ_Ａを第２のボクセルデータ１１０の画像Rの点ｘ_Ｒに変換する場合について具体的に考える。

変換関数をTとすると、

と表すことができる。

心臓の形状を剛体として考えた場合には、3次元空間内のX軸方向、Y軸方向、Z軸方向の各方向への移動（ｔ_ｘ、ｔ_ｙ、ｔ_ｚ）と、各軸の回転（α、β、γ）により変換関数が表され、下式のようになる。

心臓の形状は変化するため非剛体として考えた場合には、自由度が増えるためＴが多項式やスプライン関数などを含む複雑な関数で表される。

変換関数は、具体的には、心臓の特徴となる解剖学的な点をいくつか（例えば、３点）検出して、この点から求めるようにしてもよい。あるいは、一つの画像の輪郭を決定し、もう一つの画像の表面に相当する点列の距離が最も小さくなるように、フィティングを繰り返して変換関数を求めることもできる。あるいは、画像中の全てのボクセルの画素値を用いて類似度を調べることによって重ね合わせるようにして求めてもよい。このように画像中の全ての画素値を用いるようにすると、画素上の雑音成分が打ち消しあうので比較的安定した結果を得ることができる。（詳細は、日本放射線技術学会雑誌、第５３巻、第１号、P60〜P65、２００３年１月などを参照）。

位置合わせは、上述のいずれのレジストレーションの手法を用いてもよいが、例えば、心臓を剛体として考えて、第１、第２のボクセルデータ１００，１１０の画素値を用いて相関係数を用いて評価し、形態画像の座標系から機能画像の座標系へ変換するときの回転量、移動量などの変位量を求める。

重ね表示手段６０は、図５に示す機能画像２００上に図８の形態画像２１０を重ね合わせて表示する。位置合せ手段５０より求めた形態画像の座標系から機能画像の座標系へ変換するときの変位量を用いて、機能画像２００上の心臓の位置と形態画像２１０の心臓の位置が一致するように重ね合わせる。図９に図５に示す機能画像２００と図８の形態画像２１０を重ね合わせて表示した一例を示す。図９に示すように、重ね合わせて表示することによって、血管の位置から心臓のどのあたりの機能が正常であるか異常であるかなどを観察することができる。また、心基部に近い画像を同心円の中心付近に置き、心尖部に近い画像を同心円の外側に置くことにより、冠状動脈は心基部から心尖部に向かって分岐するので走行の状態が認識しやすくなり、解剖学的な位置と各位置での心臓の機能がどのような状態であるかを把握しやすくなる。

次に、心機能表示装置１を用いて、機能画像２００上に形態画像２１０を重ねた画像を表示するときの流れについて図１０のフローチャートに従って説明する。

まず、被写体の心臓をＣＴ装置で撮影し、撮影した画像を第１のボクセルデータ１００として第１の記憶手段１０に記憶する（S100）。さらに、同じ被写体をＭＲＩ装置で撮影した画像を第２のボクセルデータ１１０として第２の記憶手段２０に記憶する（S101）。

機能画像生成手段２０で、第１のボクセルデータ１００より心臓の機能をブルズアイ表示した機能画像２００を生成する（S102）。

次に、第２のボクセルデータ１１０を用いて、形態画像生成手段４０で形態画像２１０を生成する。各スライス面で長軸が交差する点から放射状に伸びた各線l1,l2,l3,…,lm上で血管を含む一定の距離内Dのボクセルデータを用いてMIP計算を行い、MIP計算で得られたMIP値を投影して、ブルズアイ表示した形態画像２１０を生成する（S103）。

位置合せ手段５０で、第１のボクセルデータ１００の心臓と第２のボクセルデータ１１０の心臓の形状を位置合わせする。具体的には、例えば、心臓を剛体として考えて、ボクセルデータ１００，１１０の画素値を用いて相関係数を用いて評価し、機能画像の座標系から形態画像の座標系へ変換するときの回転量、移動量などの変位量を求めることができる（S104）。

さらに、重ね表示手段６０で、位置合せ手段５０で得られた変位量を用いて、機能画像２００上の心臓の位置と形態画像２１０の心臓の位置が一致するように機能画像２００の座標系を形態画像２１０の座標系に変換して、機能画像２２０と形態画像２１０とを重ね合わせて表示する（S105）。

上述では、第１のボクセルデータと第２のボクセルデータが異なる場合について説明したが、第１と第２のボクセルデータを同じボクセルデータを用いてもよい。同じボクセルデータを用いる場合には、上述の位置合わせを行なう必要がない。

また、第１のボクセルデータと第２のボクセルデータがＣＴ装置によって撮影されたものとＭＲＩ装置によって撮影されたものである場合や、共にＭＲＩ装置あるいはＣＴ装置によって撮影されたものであってもスライス間隔が異なる場合には、上述のレジストレーションの手法を用いて位置合わせが行なわれるが、心筋シンチグラフィを対象とする場合は、画像の撮影目的が異なり形態画像と機能画像との関連が小さいため、画素値の相関を用いたレジストレーションでは上手く位置合わせできないことがある。そこでこのような場合には、それぞれの画像から対応する特徴点を選択し、その点間でレジストレーションをおこない座標変換のための回転量、移動量を算出するようにしてもよい。

また、上述では機能画像２００の座標系を形態画像２１０の座標系に変換する場合について説明したが、心臓上の各位置が一致するように変換するのであれば、形態画像２１０の座標系を機能画像２００の座標系に変換しても、双方を変形するようにしてもよい。

あるいは、機能画像２００と形態画像２１０のいずれか一方の画像を生成し、位置合わせの情報を用いて他方の画像を先に生成した画像の位置に対応させて描画するようにしてもよい。

また、上述では、長軸を横切るスライス面を等間隔の同心円状に形態画像と機能画像を表示する場合について説明したが、図１１(a)示すように、心臓の表面に沿って放射線状に伸びる線Ｓ_１、Ｓ_２、Ｓ_３、Ｓ_４の距離に対応させて、図１１(b)に示すように、機能画像と形態画像の値を配置するようにして、表面形状上の距離に依存したブルズアイ表示が可能になる。このように、表面形状の距離に依存したブルズアイ表示を行うと、図１２に示すような画像が生成される。

上述では、心基部に違い画像を同心円の中心付近に置き、心尖部に近い画像を同心円の外側に置いたブルズアイ表示について説明したが、従来から行なわれているように、心尖部に近い画像を同心円の中心付近に置き、心基部に近い画像を同心円の外側に置いたブルズアイ表示をしてもよい。

上述の図９では、形態画像２１０と機能画像２００が表す心臓上の範囲が略同じ場合について説明したが、図１４に示すように、形態画像２１０上に描かれた血管の範囲が、機能画像２００に表された機能の心臓上の範囲より広い範囲を表すようにしてもよい。このように形態画像２１０が表す心臓の範囲が機能画像２００より広ければ、形態画像２１０と機能画像２００とを重ねて表示しても機能画像２００の表示範囲に限定されることなく心臓の解剖学的な位置を確認できる部分があるので、心臓の位置と機能との対応が把握しやすい。

図１４では、同心円の外側の形態画像２１０が機能画像２００より広い場合について示しているが、同心円の内側の形態画像２１０が機能画像２００より広くなっていてもよい。さらに、同心円の内側と外側の両方の形態画像２１０が機能画像２００より広くなっていてもよい。

以上、詳細に説明したように、ブルズアイ表示した心臓の機能画像上に冠状動脈などの血管の形態を表した形態画像を重ね合わせることにより、心臓の位置と機能を対比して観察することが可能になる。

本発明の心機能表示装置の概略構成図 心臓とスライス面の関係を説明するための図 ブルズアイ表示を説明するための図 心筋の厚さブルズアイ表示の関係を説明するための図 ブルズアイ表示の一例 血管を含む表面モデルの生成方法を説明するための図 形態画像の生成方法を説明するための図 形態画像の一例 形態画像と機能画像を重ね合わせた画像の一例（その１） 本発明の心機能表示装置の動作を説明するためのフローチャート 心臓の表面に沿って放射線状に伸びる線の距離に対応させてブルズアイ表示する方法を説明するための図 放射線状に伸びる線の距離に対応させてブルズアイ表示の一例 従来のブルズアイ表示を説明するための図 形態画像と機能画像を重ね合わせた画像の一例（その２）

符号の説明

１ 心機能表示装置
２ 表示装置
１０ 第１の記憶手段
２０ 機能画像生成手段
３０ 第２の記憶手段
４０ 形態画像生成手段
５０ 位置合せ手段
６０ 重ね表示手段
100 第１のボクセルデータ
110 第２のボクセルデータ
200 機能画像
210 形態画像

Claims (5)

1. 被写体を撮影して得た三次元医用画像を構成する第１のボクセルデータを用いて心臓の各位置における機能を表した機能画像を生成する機能画像生成手段と、
   前記被写体を撮影して得た三次元医用画像を構成する前記第１のボクセルデータと同一または異なる第２のボクセルデータのうち心臓の心筋外壁に沿う血管を含む領域のボクセルデータを用いて、前記血管の太さを含む該血管の形態および該血管以外の前記心筋外壁に沿う組織の形態を描画した形態画像を生成する形態画像生成手段と、
   前記機能画像の心臓の各位置と前記形態画像の心臓の各位置とを対応させて、前記機能画像と前記形態画像とを重ねて表示する重ね表示手段とを備え、
   前記機能画像が、前記心臓の心基部から心尖部に向けて伸びる長軸を横切る方向に切った複数のスライス面における心臓の機能を表す画像を同心円状に配置した画像であり、
   前記心臓の心基部に近いスライス面の心臓の機能を表す画像を前記同心円状に配置した画像の中心に近い位置に配置し、前記スライス面の位置が前記心臓の心尖部に近づくほど該スライス面の心臓の機能を表す画像を前記中心より遠い位置に配置した画像であることを特徴とする心機能表示装置。
2. 前記形態画像が、前記心臓の心筋外壁に沿う血管を含む該心筋外壁から所定の距離内に存在するボクセルデータに対してＭＩＰ処理を施して、該距離内に存在する血管を投影した画像であることを特徴とする請求項１記載の心機能表示装置。
3. 前記形態画像上に描画された血管の前記心臓上の範囲が、前記機能画像に表された機能の前記心臓上の範囲より広い範囲を表すものであることを特徴とする請求項１または２記載の心機能表示装置。
4. 被写体を撮影して得た三次元医用画像を構成する第１のボクセルデータを用いて心臓の各位置における機能を表した機能画像を生成する機能画像生成ステップと、
   前記被写体を撮影して得た三次元医用画像を構成する前記第１のボクセルデータと同一または異なる第２のボクセルデータのうち心臓の心筋外壁に沿う血管を含む領域のボクセルデータを用いて、前記血管の太さを含む該血管の形態および該血管以外の前記心筋外壁に沿う組織の形態を描画した形態画像を生成する形態画像生成ステップと、
   前記機能画像の心臓の各位置と前記形態画像の心臓の各位置とを対応させて、前記機能画像と前記形態画像とを重ねて表示する重ね表示ステップとを備え、
   前記機能画像が、前記心臓の心基部から心尖部に向けて伸びる長軸を横切る方向に切った複数のスライス面における心臓の機能を表す画像を同心円状に配置した画像であり、
   前記心臓の心基部に近いスライス面の心臓の機能を表す画像を前記同心円状に配置した画像の中心に近い位置に配置し、前記スライス面の位置が前記心臓の心尖部に近づくほど該スライス面の心臓の機能を表す画像を前記中心より遠い位置に配置した画像であることを特徴とする心機能表示方法。
5. コンピュータを、
   被写体を撮影して得た三次元医用画像を構成する第１のボクセルデータを用いて心臓の各位置における機能を表した機能画像を生成する機能画像生成手段と、
   前記被写体を撮影して得た三次元医用画像を構成する前記第１のボクセルデータと同一または異なる第２のボクセルデータのうち心臓の心筋外壁に沿う血管を含む領域のボクセルデータを用いて、前記血管の太さを含む該血管の形態および該血管以外の前記心筋外壁に沿う組織の形態を描画した形態画像を生成する形態画像生成手段と、
   前記機能画像の心臓の各位置と前記形態画像の心臓の各位置とを対応させて、前記機能画像と前記形態画像とを重ねて表示する重ね表示手段として機能させるプログラムであり、
   前記機能画像が、前記心臓の心基部から心尖部に向けて伸びる長軸を横切る方向に切った複数のスライス面における心臓の機能を表す画像を同心円状に配置した画像であり、
   前記心臓の心基部に近いスライス面の心臓の機能を表す画像を前記同心円状に配置した画像の中心に近い位置に配置し、前記スライス面の位置が前記心臓の心尖部に近づくほど該スライス面の心臓の機能を表す画像を前記中心より遠い位置に配置した画像であることを特徴とするプログラム。