(57)【特許請求の範囲】
【請求項１】
プロセッサにより実行される方法であって、
それぞれがボクセル値の３次元配列により描かれた、関心領域についての時間的に異なる複数のボリューム画像を生成する工程（４０４、４０８）と、
ボクセル値の複数の３次元配列を空間的および時間的に崩壊させる工程であって、空間的および時間的に崩壊した２次元画像を生成する工程（４１６、４２０）と、
を含み、
前記ボクセル値の複数の３次元配列を時間的に崩壊させることは、対応するボクセル位置についての最大ボクセル値および最小ボクセル値から選択された基準に従って実行されるものであり、
前記ボクセル値の複数の３次元配列を空間的に崩壊させることは、前記選択された基準に従って実行されるものである、
ことを特徴とする方法。
【請求項２】
前記画像が、磁気共鳴技術、ｘ線コンピュータ断層撮影技術および核医学撮像技術の中から選択された撮像方法により生成された医学診断画像である、請求項１に記載の方法。
【請求項３】
前記空間的および時間的に崩壊させる工程は、
第１に、ボクセル値の上記複数の３次元配列を時間的に崩壊させて、時間的に崩壊したボクセル値の３次元配列を生成し、第２に、前記時間的に崩壊したボクセル値の３次元配列を選択された方向に投影して、２次元画像を生成する（ａ）工程、および、
第１に、前記時間的に異なるボリューム画像のそれぞれを選択された方向に空間的に投影して、時間的に異なる複数の２次元画像を生成し、第２に、前記複数の２次元画像を時間的に崩壊させて、空間的および時間的に崩壊した２次元画像を生成する（ｂ）工程
のうち１つの工程を含む、請求項１または請求項２に記載の方法。
【請求項４】
前記ボリューム画像が、タイムオブフライト磁気共鳴アンジオグラフィック技術、および、位相コントラスト磁気共鳴アンジオグラフィック技術のうち１つの技術を用いて、生成される、請求項１から請求項３のいずれかに記載の方法。
【請求項５】
前記時間的に異なるボリューム画像が生成される間、監視された肉体的な動きに従って、前記関心領域からデータを取得することを制御する工程をさらに含む、請求項１から請求項４のいずれかに記載の方法。
【請求項６】
前記関心領域の１つ以上のボクセル位置について時間経過情報を記録する工程であって、前記時間経過情報が、血液流量率情報、脈管壁の変化情報、造影剤の伝搬情報、および造影剤のピーク到着時間情報のうち少なくとも１つの情報を含む工程をさらに含む、請求項１から請求項５のいずれかに記載の方法。
【請求項７】
それぞれがボクセル値の３次元配列により描かれた、関心領域についての時間的に異なる複数のボリューム画像（１３０−１３４）を生成する手段（１１４、１２２、１２４）と、
ボクセル値の複数の３次元配列を空間的および時間的に崩壊させる手段であって、空間的および時間的に崩壊した２次元画像を生成する手段（１４０、１４６、２４０、２４６）と、
を具備し、
前記崩壊させる手段が、時間的に異なる複数の画像を、時間的に崩壊した単一の画像に崩壊させる時間経過投影プロセッサ（１４０、２４０）を含み、
前記時間経過投影プロセッサが、対応するボクセル位置についての最大ボクセル値および最小ボクセル値から選択された基準に従って、前記ボクセル値の複数の３次元配列を時間的に崩壊させることを実行し、前記選択された基準に従って、前記ボクセル値の複数の３次元配列を空間的に崩壊させることを実行する、ように構成されたものである、
ことを特徴とする画像診断装置。
【請求項８】
前記崩壊させる手段が、
ボリューム画像から空間的に投影された２次元画像を生成する空間投影プロセッサ（１４６、２４６）を含み、
前記時間経過投影プロセッサおよび前記空間投影プロセッサのうち第１プロセッサは、再構築プロセッサに接続され、前記時間経過投影プロセッサおよび前記空間投影プロセッサのうち第２プロセッサは、時間的に崩壊されかつ空間的に投影された２次元画像を生成する前記第１プロセッサに接続されている、請求項７に記載の画像診断装置。
【請求項９】
前記生成する手段は、
試験領域（１０４）に磁場を生成する磁石（１０２）と、
前記試験領域にある選択された磁気双極子に磁気共鳴を励起させ操作する無線周波数トランスミッタ（１１０）と、
前記試験領域を横切る磁場勾配を選択的に形成する勾配磁場コイル（１０８）と、
前記無線周波数の磁気共鳴信号を受信および復調する受信器（１１４）と、
関心領域についての異なる時間に得られた複数のボリューム画像を生成する再構築プロセッサ（１２４）と、
を含む、請求項７または請求項８に記載の画像診断装置。
【発明の詳細な説明】
【０００１】
本発明は、画像診断法の分野に関する。画像診断法によって、磁気共鳴アンジオグラフィー（ＭＲＡ）技術および装置、ならびに、これらにより得られたイメージデータに対してディジタル処理を行う方法および装置に対して、特定のアプリケーションがもたらされている。ここでは主に磁気共鳴アンジオグラフィ（血管造影）に関連して本発明を説明し記載するが、本発明をその他の撮像物理療法および人体以外の対象物に合わせて修正できることも認識されるであろう。
【０００２】
一般にＭＲＩでは、ほぼ一様で時間的に一定な主磁場（Ｂ0）が、被撮像物または被試験物が置かれる試験領域に設定される。主磁場を用いて別の方法で調整された対象物における、選択された磁気双極子は、磁気共鳴無線周波数（ＲＦ）励起および操作を介して、磁気共鳴を励起するために傾けられる。この磁気共鳴は、一般的には、上記対象物における選択された領域から、検出可能な磁気共鳴エコーを誘導するために操作される。撮像時には、上記エコーは、主磁場に設定された磁気の勾配を介して空間的に符号化される。ＭＲＩスキャナーからのローデータは、一般的にｋ空間として知られている行列に集約される。逆フーリエ変換、２次元フーリエ変換、３次元フーリエ変換またはその他の変換を用いることにより、上記対象物のイメージ画像が、上記ｋ空間から再構築される。
【０００３】
脈管系を撮像するために、多くのＭＲアンジオグラフィー技術が開発されている。タイムオブフライト（ＴＯＦ）技術は、動いているスピンと静止しているスピンとを区別するために、スピンの横励起とこれによる磁気共鳴信号の取得との間の時間間隔を専ら利用している。この時間間隔の間、新規のスピンは、上記磁気共鳴信号が得られる領域に移動し、励起されたスピンは、上記領域の外へ移動する。対照的に、静止しているスピンは、ＲＦ励起とデータ取得との間の間隔において固定されたままとなり、この結果、静止したスピンにより生成された磁気共鳴信号の振幅は、移動しているスピンにより生成された磁気共鳴信号の振幅と実質的に異なる。このような磁気共鳴信号から画像が再構築されるときには、移動しているスピンに対応する画像ピクセルは、シーケンスによって、静止しているスピンに対応する画像ピクセルよりもずっと明るいかまたはずっと暗いかのいずれかとなる。この方法によれば、結果として得られる画像中で、移動する血液を運ぶ脈管系は、周囲の静止しているまたはゆっくり移動している組織よりも、明るくまたは暗く見える。
【０００４】
位相コントラスト技術は、移動しているスピンにより生成された磁気共鳴信号の位相が、静止しているまたはゆっくり移動しているスピンにより生成される磁気共鳴信号の位相と異なる、という事実を専ら利用している。位相コントラスト方法は、結果として得られる磁気共鳴信号の位相をスピン速度の関数に変調する、磁気共鳴パルスシークエンス間の磁場勾配を専ら用いる。したがって、再構築される画像におけるピクセルのコントラストまたは輝度を制御するために、上記磁気共鳴信号の位相を制御することができる。血液は比較的速く移動するので、脈管系は、結果として得られる画像においてより明るくまたはより暗く見える。
【０００５】
コントラストエンハンストＭＲアンジオグラフィーは、ＭＲアンジオグラフィーの診断能力を高めるために用いられている。コントラストエンハンストＭＲアンジオグラフィーでは、ガドリニウムのような造影剤が、スキャンに先行して患者に注入される。このような注入の時間は、ｋ空間の中心線（この中心線は画像のコントラストを支配する）がピーク動脈強化（peak arterial enhancement）の間、すなわち、ボルス状の造影剤が対象となる脈管構造を通って流れる時間に得られるように、注意深く計測される。
【０００６】
ＭＲデータの取得に用いられる方法に関係なく、ＭＲＩは、ボクセル明暗度の３次元（３Ｄ）配列を備えるボリューム画像データを生成するために、用いられうる。複数の隣接したスライスに適用される３次元パルスシーケンスまたは２次元（２Ｄ）パルスを用いて、このようなボリューム画像データを取得することができる。また、対象となる関心領域を網羅するために、３次元パルスシーケンスを複数のサブボリュームに適用することも可能である。
【０００７】
最大値投影法（ＭＩＰ）は、３次元ボリューム画像データセットをレンダリングするための一般的かつ有効的な方法であり、ＭＲアンジオグラフィー画像に対して特に有用である。投影画像は、脈管組織、ならびに、狭窄、テローム性動脈硬化および動脈瘤のような病気を検査する際に、特に有用である。臨床診断では、一般的に、投影画像を見ることは、それぞれのＭＲＩスライス画像を見ることよりも好ましい。
【０００８】
ＭＩＰ画像は、本質的には、３次元ボリュームを、指定されたビューイング点から見た２次元平面に、または、指定されたビューイング角度もしくは投影角度に沿った２次元平面に、投影したものである。ＭＩＰ画像を生成する際には、光線が、画像データ点の３次元配列を通って、２次元投影画像平面における各ピクセルに投影される。上記放射線に沿ったデータ点における最大値を有するデータ点の値が選択される。この結果、各放射線に対して選択された値は、結果として得られる２次元画像における対応するピクセルのグレースケールを制御するために用いられる。様々なビューイング角度またはビューイング位置から見たＭＩＰ画像を取得することができるので、放射線専門医は、患者を柔軟に診断することができる。同様に、最小値投影法（ＭｉｎＩＰ）は、黒血の血管造影に適用されている。ＭｉｎＩＰでは、流れる血からの共鳴反応は、結果として得られる画像において血液が黒く描かれるように、ＭＲシーケンスにより最小化される。各放射線に沿うボクセル明暗度のうちの最小値は、結果として得られる２次元投影画像における対応するピクセルに割り当てられる。
【０００９】
コントラストアンジオグラフィーでは、ボルス状の造影剤は、画像領域中を力学的に移動する。また、血液の流れは、心周期に従って周期的に速くなったり遅くなったりする。ＭＲ画像は、ちょうど良い時間にこの動きをスナップショットとして捉えられたものである。何度かの試験では、多数の画像が、ちょうど良い時間にまたは心周期の位相において得られている。この後、関心領域がシネマティックディスプレイに表示される。
【００１０】
現在のＭＩＰ技術には、例えば、心周期や造影物質の到着および出発等に影響される、血流作用における時変効果（time varying effect）が含まれていない。このＭＩＰ技術では、むしろ、得られる血流作用が静止しているものと仮定している。脈管の縁の鮮明度を向上させるために、再構築グリッド再位置決め（reconstruction grid repositioning）や様々な補間技術のような多くのアルゴリズムが開発されているが、このような技術では、上述した時変効果に起因する誤りを補正することができない。
【００１１】
本発明の第１の態様では、画像診断方法は、それぞれがボクセル値の３次元配列により描かれた、関心領域についての時間的に置換された複数のボリューム画像を生成する工程（４０４、４０８）と、選択された基準に従って、ボクセル値の複数の３次元配列を空間的および時間的に崩壊させる工程であって、空間的および時間的に崩壊した２次元画像を生成する工程（４１６、４２０）と、を具備する。
【００１２】
本発明の第２の態様では、それぞれがボクセル値の３次元配列により描かれた、関心領域についての時間的に置換された複数のボリューム画像（１３０−１３４）を生成する手段（１１４、１２２、１２４）と、選択された基準に従って、ボクセル値の複数の３次元配列を空間的および時間的に崩壊させる手段であって、空間的および時間的に崩壊した２次元画像を生成する手段（１４０、１４６、２４０、２４６）と、を具備する。
【００１３】
本発明の一利点は、本発明が、撮像される血管における血管管腔を描く際の精度を向上させることである。
【００１４】
本発明の別の利点は、本発明が、生体内での血流時間の変化を捉えることである。
【００１５】
本発明の別の利点は、本発明が、例えば変動する管腔の大きさのような血管壁の変動を補償することである。
【００１６】
本発明の別の利点は、本発明を無数の診断画像取得技術に迅速に適用できる点にある。
【００１７】
添付図面を参照した例を用いて、本発明を実施するための１つの方法を詳細に説明する。
【００１８】
図１を参照するに、磁気共鳴画像システム１００は、試験領域１０４を通るほぼ一様でありかつ時間的に一定な磁場Ｂ0を形成する超伝導永久磁石１０２を含む。図１には開磁石システム（open magnet system）が示されているが、本発明は、ボアタイプ磁石を用いたような、他の周知なタイプのＭＲＩスキャナにも、同様に適用される。本発明は、また、コンピュータ断層撮影、核医学およびディジタルＸ線等のような撮像物理療法にも適用される。勾配パルス増幅器１０６は、選択された数のまたは選択された対の全体勾配コイル１０８に対して、電流パルスを適用することにより、試験領域１０４のｘ、ｙおよびｚ軸に沿った磁場勾配を形成する。無線周波数トランスミッタ１１０は、任意的にはディジタルトランスミッタであり、全体無線周波数コイル１１２に接続され、試験領域１０４にＲＦパルスを発信する。
【００１９】
第２セットの、磁石コイル、勾配コイルおよび無線周波数コイルが、下部の磁極片として患者の下に設けられている。全体無線周波数コイル１１２またはこれに代えて（図示しない）表面コイルアレイが、結果として得られる共鳴信号を復調するための受信器１１４に接続されている。
【００２０】
シーケンス制御部１１６は、勾配パルス増幅器１０６およびトランスミッタ１１０を制御して、磁気共鳴アンジオグラフィースキャニングシーケンスを生成する。さらに具体的には、上記ディジタルトランスミッタは、人体における選択されたボリューム中で、磁気スピンを反転または励起させ、磁気共鳴を誘導し、磁気共鳴の重点を変え、磁気共鳴を操作し、空間的にまたは別の方法により磁気共鳴を符号化し、スピンを飽和させることを目的として、対象物に適用される一連のＲＦパルスおよび磁場勾配パルスを生成させられる。生成されるシーケンスは、３次元画像に積み重ねられる複数のパラレル２次元画像データを生成するための、従来の何らかの３次元または２次元の脈管系撮像シーケンスであろう。このＭＲＩ方法は、造影剤を用いるかまたは用いない、何らかの力学的方法（dynamic method）または何らかのゲート制御撮像方法（gated acquisition method）であろう。具体例としてのゲート制御撮像方法は、トリガとして例えばＱＲＳ複素数のような電気カルジオグラム信号を用いる。具体例としてのコントラストＭＲＡ技術は、Keyhole、TRICK、Spiral interleaved、および、centric encoding等を含む。上記シーケンスは、一般的には、例えばシーケンス制御部１１６のシーケンス制御メモリに格納された、ユーザ選択可能であって事前にプログラムされた複数の画像シーケンスのうちの１つであろう。
【００２１】
電気カルジオグラフのゲーティング（gating）が用いられた一実施形態では、選択的に用いられる電気カルジオグラフセンサ（ＥＣＧ）１１８は、患者１５０に取り付けられ、この患者の電気的なＥＣＧ信号を出力する。ＥＣＧゲートパルス生成器１２０は、シーケンス制御部１１６に対するＥＣＧゲーティング信号としてゲートパルスを出力する。このゲートパルスは、上記ＥＣＧ信号のＲ波ピーク値に同期した非常に小さいパルス幅を有する。このゲートパルスは、例えば開始時間およびその他のシーケンススキャニングパラメータを制御してＥＣＧゲートスキャンを生成するために、ＥＣＧゲーティング信号としてシーケンス制御部１１６により用いられる。
【００２２】
これらの信号は、受信器１１４により受信および変調され、ｋ空間データメモリ１２２における標本とされる。メモリ１２２内のデータは、従来技術を用いた再構築またはアレイプロセッサ１２４により、ボリュームデータとして再構築される。アレイプロセッサ１２４は、Ｎ個のボリューム画像メモリ、すなわち、メモリ１３０、１３２、．．．、１３４のそれぞれについて同じように選択されたボリュームに対して、ボリューム画像データを供給する。ここで、Ｎは、得られるボリューム画像データセットの数であり、少なくとも２つである。
【００２３】
一実施形態では、例えば、上記ボリューム画像データが２次元マルチスライス方法により生成されるのであれば、第１ボリューム画像メモリ１３０に供給される、各スライスについてのピクセル明暗度の２次元配列を生成するために、各データスライスは、再構築プロセッサ１２４により別々に２次元フーリエ変換がなされうる。マルチプルスライスについてのピクセル値を集めることにより、ボクセル明暗度の３次元配列が構築される。
【００２４】
３次元タイムオブフライト技術（3D time-of-flight techniques）、３次元位相コントラスト（3D phase contrast）または３次元位相感度技術（3D phase sensitive techniques）等を用いた、別の実施形態では、ｋ空間の３次元配列データは、データメモリ１２２に読み込まれる。この３次元ｋ空間配列は、プロセッサ１２４により３次元逆フーリエ変換がなされ、ボリューム画像メモリ１３０によりボクセル明暗度の３次元配列として格納される。
【００２５】
さらに別の実施形態では、逆フーリエ変換が各受信データラインに対して実行される。この各受信データラインは、変換データラインのすべてのセットが生成されるまでメモリに格納されている。この後、上記データに対して位相符号化方向への第２の逆フーリエ逆変換が実行され、メモリに格納されるべきスライスデータが生成されうる。撮像ボリューム全体についてのボクセル明暗度がボリュームメモリ１３０に格納されるまで、上記工程は、各スライスについて繰り返される。一実施形態では、３次元ボリュームイメージングは、関心領域を形成するためにこの後に結合される複数の隣接したサブボリュームに対して、用いられるであろう。
【００２６】
上述したスキャニング手順およびデータ取得手順とは無関係に、選択された領域についてのボリューム画像Ｍ1、Ｍ2、・・・、ＭNがＮ個の異なる時間に得られ、複数のボリューム画像データセットが、Ｎ個のボリュームメモリ、すなわち、メモリ１３０、１３２、．．．、１３４に格納される。
【００２７】
時間分解明暗度追跡（１４０）および空間明暗度投影（１４６）については、簡略化のために、本明細書では主に最大値追跡法および最大値投影法を参照して説明する。しかしながら、本発明では、血脈管が暗く見えるようなＭＲＡアプリケーション用の時間分解追跡および空間投影についての最小明暗値を用いることも考慮されており、本明細書での記載は、これらを用いることにも同様に適用されうる。メモリ１３０、１３２、．．．１３４内に格納されている複数のボリューム画像データセットは、時間的に置換されたメモリ１３０、１３２、．．．１３４内の各撮像ボリュームの対応するボクセルについての最大明暗度ボクセルを探し出す、時間経過投影プロセッサ１４０により、時間的に崩壊される。メモリ１３０、１３２、．．．、１３４に格納されたボクセル明暗度を比較し、かつ、時間的に崩壊した（時間崩壊）ボリューム画像データセットに対して、各空間位置ごとに見つけられた最大明暗値を割り当てることにより、選択されたボリュームについての時間的に崩壊した（時間崩壊）ボリュームデータセットが、生成されボリュームメモリ１４４に格納される。
【００２８】
一実施形態では、本来の関心領域内のすべての空間座標について、時間分解データが集められるとともに時間経過投影が実行される。しかしながら、本来の関心領域全体のすべての空間座標に対して、時間分解データを取得すること、およびまたは、時間経過投影を実行することは、必要ではない。別の実施形態では、時間分解データを取得すること、およびまたは、時間経過投影は、本来の関心領域より小さくかつこの関心領域内に含まれる、１つ以上の特定の関心サブボリュームに限定される。この結果、いくつかの実施形態では、３Ｄメモリ１４４は、関心サブボリューム内の選択された空間座標のみについての時間的に崩壊した（時間崩壊）データを受信する。時間非分解値は、上記関心サブボリュームの外側にある本来の関心領域内の他のすべての領域に対して、ロードされる。この方法によれば、データを取得するための時間およびまたは画像を再構築するための時間を削減することができる。
【００２９】
一実施形態では、時間経過投影プロセッサ１４０は、さらに、撮像されたボリュームについての時間的な情報および事象を記録するためのメモリまたはデータベース１４２を任意的に含む。このような時間的な情報を、医療診断で得られた血管造影図を補うために用いることが可能である。このような情報は、例えば、（コントラスト技術のための）ボルス伝搬の時間経過（time course of bolus propagation）およびボルス到着のピーク時間（peak time of bolus arrival）、（非コントラストおよびまたはゲーティング技術）新規かつ高速な血液の到着（fresh/fast blood arrival）、ならびに、動脈および静脈のような脈管間の時間差等を含む。このような情報自体を様々な肉体的および病理学上の条件の医療診断に用いることができ、これにより、必要な血管造影図に加えて補助的な診断情報を得ることができる。例えば、近年のおよびまたは緩慢に脈管が充満することによって、閉塞またはその他の循環系機能不全の可能性を認識することができる。時間経過情報のデータベース１４２を今後参照するときのために格納することもできる。一実施形態では、脈管特性的な血液の流量率情報のうちの患者特性的な情報を格納するために、データベース１４２を用いることもでき、さらに、一実施形態では、シーケンス制御部１１６により例えば次の取得のためのスキャンタイミングを調整するために、データベース１４２を用いることもできる。
【００３０】
投影プロセッサ１４６は、選択されたビューイング点からの、または、選択されたビューイング方向もしくはビューイング軸に沿った、時間的に崩壊（時間崩壊）したボリュームデータの最大値投影（ＭＩＰ）を形成し、２次元投影画像は、投影画像メモリ１４８に格納される。ビデオプロセッサ１５２は、投影画像メモリ１４８内の画像を、ビデオモニタ１５４のＭＩＰ画像として表示するのに適当なフォーマットに変換する。
【００３１】
本発明の時間分解明暗度追跡ＭＩＰ技術の一実施形態では、ＭＩＰプロセッサ１４６は、アンジオグラフィに適した付加的機能を供給する追加モジュールを有する１モジュールの画像プロセッサである。ここで、上記追加モジュールは、角度選択、リクロップ（recrop）、背景抑圧、補間、再構築グリッド再位置決め、およびこれらと同等な他のものを実現するためのプロセッサおよびプログラムモジュールを含む。
【００３２】
ビデオプロセッサ１５２は、また、３次元時間崩壊画像メモリに直接接続されており、このメモリから直接画像データを選択する。このビデオプロセッサは、様々なスライス画像またはこれと同等な他のものを選択することができる。このビデオプロセッサまたは組み合わさった複数のプロセッサ（図示しない）は、表面レンダリング（surface rendering）および他の従来方式の画像表示を処理することができる。
【００３３】
図１に示す実施形態では、時間経過投影１４０は、空間投影１４６に先行して実行されている。しかしながら、時間分解明暗度追跡および最大値投影を実行する順序を逆にすることが可能であることが認識されるであろう。図２は、図１に示した実施形態に類似した、本発明の別の実施形態を表しているが、この図２において、ボリュームメモリ１３０、１３２、．．．１３４に格納されている各ボリューム画像は、投影プロセッサ２４６を介して最大値投影法が実行され、２次元画像メモリ２３０、２３２、．．．２３４のうちの対応する番号のメモリに格納される、時間的に置換された複数の２次元投影画像が生成される。このように生成された２次元投影の集まりは、この後、時間経過投影プロセッサ２４０を用いて時間的に崩壊され、メモリ２４８に格納されるべき、時間的および空間的に投影された最終的な２次元画像とされる。なお、ビデオプロセッサ１５２は、利用者が投影画像の時間展開をシネマティック表示させることを選ぶことができるように、シネプロセッサ２６０により投影画像メモリ２３０、２３２、．．．２３４に接続されている。
【００３４】
図３は、時間経過投影プロセッサ１４０（図１）または時間経過投影プロセッサ２４０（図２）を適用した様子を示しており、試験状態にある血管の真腔についてのボクセル値の典型例を示している。動脈のような第１血管３００は、心臓のポンピング動作とともに膨張し、この第１血管の半径は、時間ｔ3において最大値となる。鰓動脈のような第２血管３０２の半径は、時間ｔ5経過後の短時間に心拍により最大値にまで膨張する。閉塞した動脈のような第３血管３０４は、心周期を通してほぼ一定の半径値を維持する。時間ｔ1で得られる単一の画像により、閉塞した動脈を特定することは効果的ではないであろう。時間ｔ5または時間ｔ5で得られる単一画像により、診断を行うのは困難であろう。最大明暗値画像または最小明暗値画像３０６において、各脈管の半径が最大となっているように見える。この技術は、また、その他の血流変化、造影剤の到着および出発、心周期の効果、管腔の大きさの変動、および、血管壁のその他の変化などのような様々な時変効果を含む。
【００３５】
ここで図４を参照するに、図１に示したＭＲＩシステムで実行することが可能な教示にかかる方法４００の一例が示されている。工程４０４では、ボリューム画像データが得られる。このボリュームデータは、好ましくは磁気共鳴データであり、より好ましくは磁気共鳴アンジオグラフィデータである。しかしながら、コンピュータ断層撮影法またはその他の容積測定画像療法のようなその他の物理療法から得られた体積データについても、用いることが考慮されている。ボリューム画像データが得られた後、工程４０８では、アレイプロセッサ１２４を用いてボリュームデータが再構築され、再構築された電子的画像は、メモリ１３０、１３２、．．．、１３４のうち１つのメモリに格納される。工程４０４および工程４０８は、それぞれが異なる時間に対応する、関心領域についての少なくとも２つのボリューム画像が生成されるまで繰り返される。例えば２、３、４、５等の所望値または所定値のデータセットが得られた際には、工程４１２で決定されているように、複数のボリュームセットが工程４１６で時間的に崩壊される。時間的崩壊工程では、時間経過投影プロセッサ１４０は、時間的に置換された複数のボリューム画像における各空間位置または対応するボクセルに対して、最大値明暗度を割り当てて、時間的に崩壊した３次元画像を生成する。ディスプレイ１５６にレンダリングするための最大値（または最小値）投影技術を用いて３次元画像を２次元平面に投影することにより、工程４２０では、ＭＲＡ空間投影画像が、上記時間的に崩壊した３次元画像から生成される。図５は、図４に示した実施形態に類似した、本発明の別の実施形態にかかる方法５００の一例を示しているが、この図５では、工程４１６および工程４２０の順序が逆となっている。
【図面の簡単な説明】
【図１】本発明の態様にかかる時間分解明暗値追跡を含む磁気共鳴撮像スキャナの概略を示す図
【図２】本発明の磁気共鳴撮像スキャナの別の実施形態を示す図
【図３】試験状態にある血管の真腔についてのボクセル値の典型例を提供する教示にかかる時間経過投影の適用を示す図
【図４】本発明の方法の一例を概説するフロー図
【図５】本発明の方法の別例を概説するフロー図
【図１】
0004856824.tif 000002
【図２】
0004856824.tif 000003
【図３】
0004856824.tif 000004
【図４】
0004856824.tif 000005
【図５】
0004856824.tif 000006