(57)【要約】【課題】実際の心臓の動きを反映して、心臓の運動機能を正確に評価する手法を提供する。【解決手段】位置取得手段２０で、所定の時間間隔で被写体の心臓を撮影して得られた３次元ボリュームデータ群１１０のうちの異なる時に撮影された２以上の３次元ボリュームデータ１００内に存在する心臓の解剖学的に対応する点の位置を求め、機能情報取得手段３０で得られた２以上の前記３次元ボリュームデータ内に存在する心臓の前記各点の位置から３次元空間内の各点における心臓の機能を表す機能情報を取得する。表示手段４０で、機能情報に応じて前記各点における心臓の運動機能を前記被写体の心臓の形状に合わせて表示する。【選択図】図１
2008289799.tif 000002
【特許請求の範囲】
【請求項１】
所定の時間間隔で被写体の心臓を撮影して得られた３次元ボリュームデータ群を記憶するボリュームデータ記憶手段と、
前記３次元ボリュームデータ群のうちの異なる時に撮影された２以上の３次元ボリュームデータ内に存在する心臓の解剖学的に対応する点の位置を求める位置取得手段と、
該位置取得手段により得られた２以上の前記３次元ボリュームデータ内に存在する心臓の前記対応する各点の位置から、３次元空間内の各点における心臓の機能を表す機能情報を取得する機能情報取得手段と、
該機能情報取得手段により得られた前記機能情報に応じて前記各点における心臓の運動機能を前記被写体の心臓の形状に合わせて表示する表示手段とを備えたことを特徴とする心臓機能解析装置。
【請求項２】
前記機能情報が、前記各点における心筋の厚さから壁厚の変化量または増加率を表すものであることを特徴とする請求項１記載の心臓機能解析装置。
【請求項３】
前記機能情報が、前記各点における拡張期の心室の径と収縮期の心室の径の差を表すものであることを特徴とする請求項１記載の心臓機能解析装置。
【請求項４】
前記機能情報が、前記対応する各点の変位量を表すものであることを特徴とする請求項１記載の心臓機能解析装置。
【請求項５】
前記表示手段が、前記心臓の形状の上に前記機能情報の値に応じて異なる色を各点に配置して表示するものであることを特徴とする請求項１から４いずれか記載の心臓機能解析装置。
【請求項６】
コンピュータを、
所定の時間間隔で被写体の心臓を撮影して得られた３次元ボリュームデータ群のうちの異なる時に撮影された２以上の３次元ボリュームデータ内に存在する心臓の解剖学的に対応する各点の位置を求める位置取得手段と、
該位置取得手段により得られた２以上の前記３次元ボリュームデータ内に存在する心臓の前記対応する各点の位置から、３次元空間内の各点における心臓の機能を表す機能情報を取得する機能情報取得手段と、
該機能情報取得手段により得られた前記機能情報に応じて前記各点における心臓の運動機能を前記被写体の心臓の形状に合わせて表示する表示手段として機能させるプログラム。
【請求項７】
所定の時間間隔で被写体の心臓を撮影して得られた３次元ボリュームデータ群のうちの異なる時に撮影された２以上の３次元ボリュームデータ内に存在する心臓の解剖学的に対応する各点の位置を求める位置取得ステップと、
該位置取得ステップにより得られた２以上の前記３次元ボリュームデータ内に存在する心臓の前記対応する各点の位置から、３次元空間内の各点における心臓の機能を表す機能情報を取得する機能情報取得ステップと、
該機能情報取得ステップにより得られた前記機能情報に応じて前記各点における心臓の運動機能を前記被写体の心臓の形状に合わせて表示する表示ステップとを備えたことを特徴とする心臓機能解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、心臓の機能を表示する心臓機能解析装置、方法およびそのプログラムに関するものである。
【背景技術】
【０００２】
従来、複数時相を用いたCT画像を対象と心臓運動解析では、それぞれの時相から三次元画像の中から短軸方向の断面画像（短軸像）を切り出し、対応する位置の断層画像を対象に各時相から心筋壁の輪郭R1,R2,R3,R4を抽出し(図８ａの拡張期の図と図８ｂの収縮期の図を参照)、長軸Ｌの同じ位置の交線上の断面における心筋壁(内腔)の輪郭R1，R3を求め、その差分量を心筋壁運動量とする処理をおこなっている(図８ｃ参照)。この処理を心基部から心尖部までの長軸に直交する各短軸段面位置でおこない、まとめたものを三次元的な解析結果としている。
【０００３】
また、このような心臓運動解析の結果を表す表示方法としてブルズアイ（Bull’s Eye）表示法が用いられている。ブルズアイ表示法は、心臓を楕円体モデルで近似して考えたとき、楕円の長軸を横切る方向に一定の間隔で切ったスライス面を同心円状に配置して表示する手法である。この表示方法は、主に機能画像を対象として用いられ、心筋シンチグラフィ（SPECT）によって直接撮影されたものや、CT/MRIから得られた解析結果から得られた機能を表わすために用いられる（例えば、特許文献１など）。
【特許文献１】特開2005-27999公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、上述のCT/MRIから得られた解析結果は、短軸像から中心からの放射方向の動きに注目して求められているが、実際の心臓の動きは膨張や収縮だけでなく回転やねじれも加わるため解析結果は実際の心臓の動きに基づいた心臓の機能を反映したものではなかった。
【０００５】
そこで、本発明では、心臓の機能を正確に評価することが可能な心臓機能解析装置、方法およびそのプログラムを提供することを目的とするものである。
【課題を解決するための手段】
【０００６】
本発明の心臓機能解析装置は、所定の時間間隔で被写体の心臓を撮影して得られた３次元ボリュームデータ群を記憶するボリュームデータ記憶手段と、
前記３次元ボリュームデータ群のうちの異なる時に撮影された２以上の３次元ボリュームデータ内に存在する心臓の解剖学的に対応する点の位置を求める位置取得手段と、
該位置取得手段により得られた２以上の前記３次元ボリュームデータ内に存在する心臓の前記対応する各点の位置から、３次元空間内の各点における心臓の機能を表す機能情報を取得する機能情報取得手段と、
該機能情報取得手段により得られた前記機能情報に応じて前記各点における心臓の運動機能を前記被写体の心臓の形状に合わせて表示する表示手段とを備えたことを特徴とするものである。
【０００７】
また、本発明のプログラムは、コンピュータを、
所定の時間間隔で被写体の心臓を撮影して得られた３次元ボリュームデータ群のうちの異なる時に撮影された２以上の３次元ボリュームデータ内に存在する心臓の解剖学的に対応する各点の位置を求める位置取得手段と、
該位置取得手段により得られた２以上の前記３次元ボリュームデータ内に存在する心臓の前記対応する各点の位置から、３次元空間内の各点における心臓の機能を表す機能情報を取得する機能情報取得手段と、
該機能情報取得手段により得られた前記機能情報に応じて前記各点における心臓の運動機能を前記被写体の心臓の形状に合わせて表示する表示手段として機能させることを特徴とするものである。
【０００８】
また、本発明の心臓機能解析方法は、所定の時間間隔で被写体の心臓を撮影して得られた３次元ボリュームデータ群のうちの異なる時に撮影された２以上の３次元ボリュームデータ内に存在する心臓の解剖学的に対応する各点の位置を求める位置取得ステップと、
該位置取得ステップにより得られた２以上の前記３次元ボリュームデータ内に存在する心臓の前記対応する各点の位置から、３次元空間内の各点における心臓の機能を表す機能情報を取得する機能情報取得ステップと、
該機能情報取得ステップにより得られた前記機能情報に応じて前記各点における心臓の運動機能を前記被写体の心臓の形状に合わせて表示する表示ステップとを備えたことを特徴とするものである。
【０００９】
「心臓の機能を表す機能情報」とは、心臓の運動機能を評価する指標となるものであって、例えば、心臓の対応する点が移動した変位量や心臓の壁厚や心室の径などから得られるものである。
【００１０】
また、前記機能情報が、前記各点における心筋の厚さから壁厚の変化量または増加率を表すものであってもよい。
【００１１】
また、前記機能情報が、前記各点における拡張期の心室の径と収縮期の心室の径の差を表すものであってもよい。
【００１２】
さらに、前記機能情報が、前記対応する各点の変位量を表すものであってもよい。
【００１３】
さらにまた、前記表示手段が、前記心臓の形状の上に前記機能情報の値に応じて異なる色を各点に配置して表示するものが望ましい。
【発明の効果】
【００１４】
本発明によれば異なる時に撮影された２以上の３次元ボリュームデータ内に存在する心臓の解剖学的に対応する点の位置を求め、２以上の心臓の対応する各点での心臓の機能を求めて、この動きに応じた心臓の運動機能を被写体の心臓の形状に合わせて表示することにより、単に膨張や収縮だけでなく回転やねじれも加わる実際の心臓の機能を反映した心臓の運動機能を評価することができる。
【００１５】
また、機能を評価するものとして、心臓の対応する点の変位量、心筋の厚さから壁厚の変化量や増加率、拡張期と収縮期の心室の径の差を求めることにより、心臓の運動機能を定量的に表すことができる。
【００１６】
また、表示する際に心臓の形状の上に前記機能情報の値に応じて異なる色を各点に配置して表示することにより、心臓の機能の良いところと悪いところを確認することができる。
【発明を実施するための最良の形態】
【００１７】
以下、図面を参照して本発明の心機能表示装置の実施の形態を詳細に説明する。図１は本発明の心臓機能解析装置の概略構成図である。なお、図１のような心機能表示装置１の構成は、補助記憶装置に読み込まれた心臓機能解析処理プログラムをコンピュータ上で実行することにより実現される。このとき、この心臓機能解析処理プログラムは、ＣＤ−ＲＯＭ等の記憶媒体に記憶され、もしくはインターネット等のネットワークを介して配布され、コンピュータにインストールされる。
【００１８】
本発明の心臓機能解析装置１は、所定の時間間隔で被写体の心臓を撮影して得られた複数の３次元ボリュームデータ１００からなる３次元ボリュームデータ群１１０を記憶するボリュームデータ記憶手段１０と、２以上の３次元ボリュームデータ１００内に存在する心臓の解剖学的に対応する点の位置を求める位置取得手段２０と、２以上の３次元ボリュームデータ１００内に存在する心臓の解剖学的に対応する点の位置から、３次元空間内の各点の心臓の機能を表す機能情報を取得する機能情報取得手段３０と、機能情報に応じて各点における心臓の運動機能を被写体の心臓の形状に合わせて表示する表示手段４０とを備える。
【００１９】
３次元ボリュームデータ１００は、診断対象となる心臓などを断層面に垂直な方向に沿って順に得られる２次元の断層像データを積層することによって得られ、CT装置やMRI装置などで撮影された複数の断層画像を重ね合わせて生成される。例えば、CT装置を用いて得られたボリュームデータには、ボクセルごとにＸ線の吸収量を蓄えたデータが得られ、各ボクセルに対して１つの画素値（CT装置で撮影した場合は、Ｘ線の吸収量を示す値）が与えられたデータとなる。
【００２０】
３次元ボリュームデータ群１１０は、複数の３次元ボリュームデータ１００で構成され、例えば、被写体を一定の時間間隔Δｔで撮影を行ない、時間ｔ１、ｔ２、・・・、ｔｎで撮影した一連の３次元ボリュームデータ１００である。
【００２１】
記憶手段１０は、ハードディスクなどの大容量の記憶装置である。記憶手段１０には３次元ボリュームデータ群１１０が記憶される。
【００２２】
位置取得手段２０は、異なる時間に撮影された２つ以上の３次元ボリュームデータ１００間で、３次元ボリュームデータ１００内に存在する心臓の解剖学的に対応する点の位置を求める。
【００２３】
具体的には、時間ｔと時間ｔ＋Δｔの２つの時相の３次元ボリュームデータ１００間で３次元のボリュームデータの画素値が一致するように位置合わせをすることによって、心臓の解剖学的に対応する点の位置を求め、さらに、時間の経過に沿って、２つ以上の時相の３次元ボリュームデータ１００内の位置合わせを順次行うことで、２以上の時相の心臓の解剖学的に対応する点の位置を求めることができる。
【００２４】
位置合わせには、例えば、３次元空間を変形させることで３次元空間内の画像を変形させて位置合わせする手法を用いることができる。まず、図２ａに示すように３次元空間内に、一定の間隔で３次元の格子状にB-Splineの制御点（図２ａの四角）を配置し、この制御点を矢印の方向へ移動させることで、図２ｂのように３次元空間を変形させ、図３に示すような時間ｔにおける心臓の形状（実線）と時間ｔ＋Δｔにおける心臓の形状（破線）が一致するように３次元空間を変形させる（図４参照）。具体的には、２つの画像間の差の自乗が小さくなるような場所を求めるＳＩＤ（squared intensity difference）の手法や、相関係数が高くなるような場所を求めるＣＣ（correlation coefficient）の手法を用いて最も類似した場所を検索して位置合わせを行う。この位置合わせの結果に基づいて心臓の解剖学的に対応する点の位置を求めることができる。ここでは、B-Splineの制御点を移動させることで３次元空間を変形させる場合について説明したが、他の補間式を用いて空間を記述して３次元空間を変形するようにしてもよい。
【００２５】
機能情報取得手段３０は、２以上の３次元ボリュームデータ１００から、心臓の解剖学的に対応する各点から心臓の機能を表す機能情報を取得する。心臓の運動機能は、膨張や収縮を繰り返す心臓の動きによって評価されるが、心臓は中心からの放射方向へ動くだけでなく実施には心臓に回転やねじれが加わって動いている。心臓の機能は、心臓の対応する点が移動した変位量、心室の内径、心筋の厚さなどから評価することができるので、各時相の心臓の対応する点の変位量や心室の内径や心筋の厚さから機能情報を求める。
【００２６】
機能情報としては、具体的には、心壁の対応する各点の変位量、拡張期の心筋の厚さから収縮期の心筋の厚さの差を表した壁厚変化量、拡張期の心筋の厚さをA、収縮期の心筋の厚さをBとした時、B-A/Aで計算された値を表した壁厚増加率、拡張期の心室の径から収縮期の心室の径の差を表した壁運動量などを求める。
【００２７】
表示手段４０は、機能情報に応じて各点における心臓の機能を心臓の形状に合わせて表示する。具体的には、心壁の各点の変位を表す三次元のベクトルのx,y,z成分の大きさに対応して色成分r,g,b成分を割り当てて、実際の心臓の形状の上に、各点の機能情報の値に対応した色を配置して動きの方向を可視化する（図５参照。図5では、色の変わったところをハッチングを変えて表している。）。他の機能情報についても、心壁上の各点の機能情報の値に対応する色を配置して心臓の機能を可視化することができる。このように心臓の形状の上に機能に応じた色を配置することによって、心臓の各位置がどの程度機能しているかを即座に把握することができる。また、色に限らずテクスチャーを変えても表示してもよいし、機能情報を表す数値を心臓の各対応位置の上に表示してもよい。
【００２８】
次に本発明の心臓機能解析装置を用いて心臓の機能を表した画像を生成する処理の流れについて図６のフローチャートに従って説明する。
【００２９】
診断対象となる被写体をCT装置などにより時間間隔Δｔで撮影した３次元ボリュームデータ群１１０をハードディスク（記憶手段１０）に記憶する（S100）。
【００３０】
次に、位置取得手段２０で、これらの３次元ボリュームデータ群１１０内の各３次元ボリュームデータ内の心臓の解剖学的に対応する点の位置を求める（S101）。心臓の拡張期と収縮期によってねじれが加わるため、心臓の解剖学的に対応する各点は、図７に示すようにねじれた位置に存在する。
【００３１】
そこで、機能情報取得手段３０で、２以上の３次元ボリュームデータ内の心臓の解剖学的に対応する各点の位置、心筋の厚さ、心室の径から、機能情報として各対応する点の変位量・壁厚変化量・壁厚増加率・壁運動量などを求める（S102）。表示手段４０は、機能情報取得手段３０で得られた機能情報の値に応じた色を心臓の形状上の各点に配置した画像を表示する（S103）。
【００３２】
以上、詳細に説明したように、心臓の膨張や収縮だけでなく回転やねじれも加わる実際の心臓の動きを反映して評価された心臓の運動機能を表示することが可能になる。
【図面の簡単な説明】
【００３３】
【図１】本発明の心臓機能解析装置の概略構成図
【図２ａ】３次元空間を変形させる方法を説明するための図（その１）
【図２ｂ】３次元空間を変形させる方法を説明するための図（その２）
【図３】位置あわせを説明するための図（その１）
【図４】位置あわせを説明するための図（その２）
【図５】心臓の形状に合わせて機能情報に応じた色を配置した画像の一例
【図６】心臓機能解析装置の処理の流れを説明するための図
【図７】心臓の解剖学的に対応する点が動いた位置を表す図
【図８ａ】心臓の拡張期の図
【図８ｂ】心臓の収縮期の図
【図８ｃ】心臓の拡張期と収縮期の心筋壁の位置と心臓の機能を説明するための図
【符号の説明】
【００３４】
１    心臓機能解析装置
２    表示装置
１０  ボリュームデータ記憶手段
２０  位置取得手段
３０  機能情報取得手段
４０  表示手段
100  ３次元ボリュームデータ
110  ３次元ボリュームデータ群
【図１】
2008289799.tif 000003
【図２ａ】
2008289799.tif 000004
【図２ｂ】
2008289799.tif 000005
【図３】
2008289799.tif 000006
【図４】
2008289799.tif 000007
【図５】
2008289799.tif 000008
【図６】
2008289799.tif 000009
【図７】
2008289799.tif 000010
【図８ａ】
2008289799.tif 000011
【図８ｂ】
2008289799.tif 000012
【図８ｃ】
2008289799.tif 000013