JPH10512482A - 画像アーティファクトを抑制したデジタル消去磁気共鳴血管造影法 - Google Patents

画像アーティファクトを抑制したデジタル消去磁気共鳴血管造影法

Info

Publication number
JPH10512482A
JPH10512482A JP10508762A JP50876298A JPH10512482A JP H10512482 A JPH10512482 A JP H10512482A JP 10508762 A JP10508762 A JP 10508762A JP 50876298 A JP50876298 A JP 50876298A JP H10512482 A JPH10512482 A JP H10512482A
Authority
JP
Japan
Prior art keywords
dimensional
image
nmr
contrast
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP10508762A
Other languages
English (en)
Other versions
JP3130320B2 (ja
Inventor
エイ. ミストレッタ,チャールス,
グリスト,トーマス,エム.
フライン,リチャード
コロセック,フランク
Original Assignee
ウイスコンシン アルムニ リサーチ ファウンデーション
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ウイスコンシン アルムニ リサーチ ファウンデーション filed Critical ウイスコンシン アルムニ リサーチ ファウンデーション
Publication of JPH10512482A publication Critical patent/JPH10512482A/ja
Application granted granted Critical
Publication of JP3130320B2 publication Critical patent/JP3130320B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/281Means for the use of in vitro contrast agents
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5619Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences by temporal sharing of data, e.g. keyhole, block regional interpolation scheme for k-Space [BRISK]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56308Characterization of motion or flow; Dynamic imaging

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Vascular Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

(57)【要約】 対象物の動的なMRAの検討は、三次元の高速傾斜呼び戻しエコーパルスシーケンスを利用して実行される。結果として生じる一続きの再構築された各画像のフレーム率は、k空間の中央領域をk空間の周辺領域よりも高率でサンプリングすることで増大する。コントラスト仲介物が関心のある領域に入り込むときに、信号強度の変動により生じたアーティファクトは、収集されたデータの再正規化により減少される。どの三次元データセットが診断画像を再構築するのに利用されるべきかを、オペレータが選択できるように、k空間の中心近傍を通過する各平面を利用して、二次元の画像フレームが再構築される。

Description

【発明の詳細な説明】 画像アーティファクトを抑制したデジタル消去磁気共鳴血管造影法 [発明の背景] 本発明の背景は、磁気共鳴血管造影法(MRA)にあり、特に、NMR信号を 強調するコントラスト仲介物を用いた人体血管の動的な検討にある。 人体血管の診断検討は、多くの医学的な応用物を有する。デジタル消去血管造 影法(DSA)のようなX線画像方法は、心臓および一体化した血管を含む心臓 血管系における視覚化の広範囲な利用を見出した。腎臓の動脈や静脈,および首 部や頭部における頚動脈や静脈内の循環をあらわす画像は、計り知れない診断の 有効性をもつ。しかし不幸なことに、これらのX線画像方法は、患者が潜在的に 有害な電離放射線を被ることになるとともに、画像化されるべき血管内にコント ラスト仲介物を注入するために、侵入するカテーテルを利用する必要がある。 これらのX線技術による利点の一つに、画像データが高率(すなわち、時間的 に高い解像度)で得られるので、コントラスト仲介物を注入している間に、連続 した画像を得ることができることにある。こうした「動的な検討」は、コントラ スト仲介物の塊が関心のある血管を通過する画像の選択を可能にする。連続する 画像の初期のものは、推測される血管に十分なコントラストを有していない。ま た、後期の画像では、コントラスト仲介物が血管に達して、周囲組織に放散する ので、分析を行なうのが難しくなる。「リアルタイムデジタルX線消去イメージ ング」と表題された米国特許第4,204,225号公報に開示されるような消去方法は 、こうした画像における診断の有効性をかなり強調するのに利用されることがあ る。 磁気共鳴血管造影法(MRA)は、人体血管の画像を生成するのに、核磁気共 鳴(NMR)現象を利用する。人体組織のような個体が一様な磁場(極性場 B0 )に晒されると、組織内における個々のスピンによる磁気モーメントが、この 極性場に整列しようとするが、各々は特定のラーモア周波数でランダムな状態に 極性場の周りを歳差運動する。仮に、個体すなわち組織が、X−Y平面内にあっ て、ラーモア周波数に近い磁場(励磁場 B1)に晒されると、正味の直交磁気 モーメントMtを生成するのに、正味の整列したモーメントMzが、X−Y平面内 において回転すなわち「tipped(傾斜)」されるようになる。励起されたスピン によって信号は放出され、励起信号B1が終結した後で、この信号は受信され、 画像の形態に生成されることになる。 画像を生成するのにこれらの信号を利用する際に、傾斜磁場(Gx,Gyおよび Gz)が利用される。一般に、画像化されるべき領域は、独特な局部的方法を用い ることに従って、傾斜が変化する一続きの測定周期によりスキャンされる。各々 の測定は「view(ビュー)」のような技術を参照するとともに、ビュー数が画像 の解像度を決定する。非常によく知られている再構築技術技術の一つを用いて画 像を再構築するのに、受信したNMR信号の結果の様子すなわちビューが、デジ タル化され生成される。全体のスキャン時間は、画像に対して得られる測定周期 すなわちビューの数によってある程度決定されることから、画像の解像度を犠牲 にして、得られた画像の数を減らすことにより、スキャン時間を減らすことが可 能となる。 医学的画像を生成するのに一般的に利用される大部分のNMRスキャンは、臨 床上有用な画像のための必要なデータを得るために、多くの時間を必要とする。 このスキャン時間を分よりむしろ秒に減少することは、MRI方法を利用する臨 床上の動的検討を行なうのに大きな障害になる。最も一般的なMRI方法は通常 非トリガー用に利用されており、時間決定されたイメージングは、ピーター マ ンスフィールド氏(J.Phys.C.10: 55-58行,1977年)により初めて述べられ たような、エコー−プレーナ イメージング(EPI)パルスシーケンスを用い るためにある。EPIスキャンでは、200〜300ミリ秒オーダーで測定された期間 に現れる動的プロセスのイメージングが概ね可能となる。しかし、時間決定され たEPIは、RF励起中の長い時間間隔(例えば、100ms)によって、血管と周 囲組織との間で低いコントラストを表わすことになるため、強調されたMRAの コントラストに対し他方では相応しくないものとなる。さらに、EPIは流れに 関連したアーティファクトの変化に対して、感度を強調するとともに、k空間の T2 *変調によりEPI画像がぼやけることがある。 MRAに適用できるパルスシーケンスを用いたMRIスキャンの時間的な解像 度を増加させるのに、相当数の方法が考え出されている。「MR透視診断」なる 技術で知られており、米国特許第4,830,012号に記述される方法では、完全な画 像のために必要なN相のエンコードしたビューを、連続的にかつ繰り返して得る ことにより、対象物がスキャンされる。完全に新規な一組のNビューを待つとい うよりも、むしろ次の画像を再構築する前ではあるが、最も新らしいNビューを 利用することで、画像が遥かに高率で再構築される。言い換えれば、動的検討に おいては、以前の画像を再構築するのに利用したビューのみならず、新たに得ら れたビューから、画像が再構築される。MR透視診断では非常に高い時間率が達 成される一方、全体の画像コントラストを支配するk空間における中央のビュー が、遥かに遅い固有のスキャン率で依然として新しくされるので、MRAに対す る画像コントラストは不満足なものとなる。 MRI画像における時間的な解像度を上げる他の方法は、「keyhole(キーホ ール)」イメージングのような技術に引用される。例えば、アール.エイ.ジョ ーンズ氏などによる、SMR第11年ミーティング 1992年 要約 1138の「動的 でコントラストが強調された、k空間代用品を用いたラットの局部的大脳虚血の NMR環流イメージング」に記述されるように、一続きの画像は、コントラスト 仲介物が対象物に注入される動的な検討中に得られる。連続する画像のなかの最 初の画像は、位相エンコードのビュー全て(例えば、128ビュー)が得られる参 照ビューである。これに続く各画像は、中央のビュー(例えば、中央の32ビュー )を得るだけではあるが、生成される。これらのキーホールスキャンは、明らか に完全なスキャンに比べて遥かに速く得られるとともに、時間率は比例して増加 する。最も新しいk空間のビューを、参照スキャンからの外側で周辺にあるk空 間のビューと組み合せて用いることで、キーホール画像が再構築される。不幸な ことに、空間周波数が低い状況では、再構築された画像の変化が動的な検討の進 展を捕らえなくなり、k空間のイメージングが相応しくないものとなる。これは 、小領域におけるコントラストの変化が検討される場合の問題であり、こうした 検討中は、時間的な解像度のゲインが失われた点に対し、得られた中央ビューの 数を増やさなければならない。 k空間のキーホールイメージング方法に関連して、制限された撮像領域(FO V)の動的イメージングとして、技術的に知られた方法がある。例えば、ヒュー 氏およびパリッシュ氏により刊行された「医学における磁気共鳴,Vol.31,691 -694頁,1994年」や、フレデリクソン氏およびペルク氏による 第3SMR,1 ,197,1995年に記述されるように、この方法は、画像の変化する部分がFOV 全体の僅か片側半分を占有する動的検討に対し適用される。画像の静止部分を表 わす参照画像は、検討を開始する際に生成されるとともに、画像の動的な中央部 分だけを取り囲む一続きの画像が、位相エンコードのビュー数の半分を利用して 生成される。これらの動的なビューは、ビュー数の半分だけ(奇数或いは偶数ビ ューのいずれか一方)が必要であるため、高い時間率で得ることが可能になる。 画像の動的および静的な位置は、対応する一続きのFOV画像全体を生成するこ とで組み合わされる。勿論、画像の静的部分に各変化が現れたら、これらの変化 領域から入手した情報は、小さなFOVにアライアズ(aliased)されるアーティ ファクトを正確に取り去らなくなるであろう。 MR血管造影法は、調査するアクティブな領域が存在している。2つの基本的 技術が計画され、評価されている。第1クラスのタイムオブフライト(TOF) 技術は、周囲組織に関係する血液の動きを利用する方法からなる。大部分の共通 する検討法は、流れている血液と静止組織との間に存在する信号飽和度の違いを 利用することにある。これは、フローに関連した強調法として知られているが、 実際には血液組織におけるコントラストの改良は、多くの励起パルスを得て飽和 されるようになる静止した組織によるものなので、この効果は名前とは違ったも のになる。血液流は、励起された部分を通じて動いているが、これはより少ない 励起パルスを得てスピンにより絶えずリフレッシュされることから、飽和するの は少ない。その結果は、高信号の血液と低信号の静止組織との間で、所望の画像 コントラストになる。 MR法はまた、米国特許Re.32,701号公報に開示されるように、収集された 信号位相に対するエンコードの動きを発展させている。これらはMRA技術の第 2クラスを形成し、位相コントラスト(PC)法として知られている。一般に、 大部分のPC MRA技術は2つの画像を収集しており、各画像は、同一の速度 構成要素に対し、異なる感度を有している。血管造影法的な画像は、一対の速度 エンコード画像間の位相または複雑な違いのいずれかを配列することで、その後 得られる。位相コントラストMRA法は、三つの直交方向全ての速度構成要素に 感度を有するものに拡張されている。 近年において著しい進歩が成されているものの、多くの臨床現場では、MRA が今でも調査用の工具であると見なされているし、臨床上の患者に日常的に用い られていない。病変を隠したり、場合によっては、病変を真似ることのある多様 な有害画像アーティファクトの存在によって、TOF若しくはPC技術の一方を より広範囲に適用することが妨げられている。こうしたアーティファクトは、一 般に妥協した感度と同じような低い特質性を結果的にもたらす。 MRAの診断能力を高めるために、ガドリニウムのようなコントラスト仲介物 が、MRAスキャンの前に患者に注入されることがある。米国特許第5,417,213 号公報に記述されるように、その妙技は、コントラスト仲介物の塊が関心のある 血管を流れている瞬間に、中央のk空間ビューを得ようとするものである。これ は、日常的な臨床手続の一部分として成し遂げるには簡単なタイミングではない し、多くの解決法が提案されている。 係属中の米国特許出願08/622,718号「三次元デジタル消去磁気共鳴血管造影法 」では、動的な検討中に連続する三次元MRAデータのセット(組)を素早く得 るための方法が記述されている。この方法が、臨床的に適切な画像解像度で、所 望の時間的な解像度を提供している間、2つの問題が存在する。第1に、血管内 で失われた信号の形態を成す画像アーティファクトが生成されることである。こ うしたアーティファクトは、誤った診断を導くことがある。第2に、莫大な量の MRAデータが、動的な検討中に生成されることにあり、画像を再構築するのに 必要な生成時間が超過する。 [発明の概要] 本発明は、信号強度の変化により生じる画像アーティファクトが抑圧され、臨 床的な画像が実質的に生成する時間量を減らして再構築される、MRAの動的な 検討に対する改良を行うことにある。より特別には、本発明は、コントラスト仲 介物が対象物に注入される動的な検討の間における、連続する各三次元NMRデ ータセットを収集するための方法であり、この方法は、動的な検討の間に、コン トラスト仲介物がどのようにして収集されたNMRデータを重み付けするのかを 表わす信号C(t)を生成するとともに、画像を再構築するのに利用されるNM Rデータが等しく重み付けされるように、収集したNMRデータを正規化するこ とを含んでいる。加えて、一続きの三次元NMRデータセットの各々から三次元 画像を再構築するよりもむしろ、各々の三次元NMRデータセットの部分体を利 用して、二次元画像が生成される。オペレータは、1乃至2の三次元NMRデー タを選択するために、これらの各二次元画像を利用してもよく、この三次元NM Rデータから、最終の臨床的な画像が再構築されることになり、また、重み付け 信号C(t)を生成するのに、これらの各二次元画像における選択された関心の ある領域のデータを利用してもよい。 本発明に共通する目的は、動的な検討の間に生成された画像アーティファクト を抑制することにある。画像アーティファクトは、コントラスト仲介物が関心の ある領域を通過する間に、結果として生じる信号強度の差により生成されること がわかっている。より特別には、あらゆる特殊な画像に対しビューを収集する間 に起こる変動は、各血管内の信号ロスを生み出す。この問題の解決法は、コント ラスト仲介物が通過することによる信号強度の変化を表わす重み付け信号C(t )を生成するとともに、この重み付け信号C(t)を利用して、収集された画像 データを正規化することによって、前記変動を測定することにある。 本発明の他の共通する目的は、臨床的および診断的な画像を生成するのに必要 な処理時間を減少することにある。各々の三次元NMRデータセットを再構築す るというよりもむしろ、各々の収集された三次元NMRデータセットを利用して 、二次元画像が再構築される。各二次元画像は素早く再構築されるとともに、オ ペレータはどの画像フレームすなわち各フレームを完全に再構築すべきかを、各 二次元画像から選択できる。これらの同一の二次元画像は、重み付け信号C(t )を生成するのにも利用できる。 本発明における前述および他の目的と利点は、後に続く記述から明確になるで あろう。この記述中において、その部分を構成し、好ましい実施態様を説明する ために示された添付図面に対し参照が行なわれる。こうした実施態様は、必ずし も本発明の全範囲を示してはいないが、それ故に、本発明の範囲を解釈するため の特許請求の範囲に対し、参照が行なわれるものである。 [図面の説明] 図1は、本発明を利用したMRI装置のブロック図である。 図2は、図1におけるMRI装置の一部を構成する変換器の電気的ブロック図 である。 図3は、本発明の好ましい実施態様に利用されるパルスシーケンスのグラフ図 である。 図4は、本発明の好ましい実施態様を実行したときのデータがサンプルされる 三次元k空間のグラフ図である。 図5は、図4における三次元k空間がサンプルされたなかの整列を示すグラフ 図である。 図6は、図4における三次元k空間のサンプリングを示したグラフ図であり、 動的な検討における各画像フレームが本発明の一実施態様に従って再構築された 時間を示している。 図7は、信号強度を示すグラフ図であり、画像アーティファクトを生成する強 度の窪みを示している。 図8は、動的な検討中における各画像フレームに対するデータセットの概略図 であって、いかにしてMRA画像を生成するのに、それらのデータセットが組み 合わされるのかを示している。 図9は、動的な検討中に生成される一般的な重み信号C(t)のグラフ図である 。 図10は、図9における情報が、対応する画像データを再度重み付けすなわち 正規化するのに利用された後における、結果として生じたk空間の重み付けのグ ラフ図である。 [好ましい実施態様の説明] 最初に図1を参照すると、本発明を具体化する好ましいMRI装置の全体構成 が示されている。装置の動作は、キーボードと、制御パネル102と、表示器104と を含むオペレータコンソール100から制御される。オペレータコンソール100は、 スクリーン104上での各画像の生成と表示をオペレータが制御できるように、別 個のコンピュータ装置107とリンク116を通じて通信を行なう。コンピュータ装置 107は、お互いが背面を通して通信を行なう複数のモジュールを含んでいる。こ れらのモジュールは、画像プロセッサモジュール106と、CPUモジュール108と 、メモリモジュール113を含んでおり、画像データ配列を記憶するフレームバッ ファのような技術として知られている。コンピュータ装置107は、画像データや プログラムを記憶するディスクメモリ111やテープドライブ112に連結するととも に、高速のシリアルリンク115を通して、別個の装置制御部122と通信を行なう。 装置制御部122は、背面で共に接続された一組のモジュールを含んでいる。こ れらのモジュールは、CPUモジュール119と、シリアルリンク125を通してオペ レータコンソール100に接続するパルス発生器モジュール121とを含んでいる。こ のシリアルリンク125を通して、装置制御部122は実行されるべきスキャンシーケ ンスを指摘するオペレータからの命令を受け取る。パルス発生器モジュール121 は、所望のスキャンシーケンスを実行するのに、装置の各構成を動作させる。そ れにより、生成されるべきRFパルスのタイミング,強さおよび形を表わすとと もに、データ収集窓のタイミングおよび長さを表わすデータを生成する。パルス 発生器モジュール121は、スキャン中に生成されるべき傾斜パルスのタイミング と形を表わすのに、一組の傾斜アンプ127に接続している。また、パルス発生器 モジュール121は、電極からのECG(心電図)信号や、肺からの呼吸信号のよ うな、患者に接続された複数の異なるセンサからの信号を受け取る生理的収集コ ントローラ129からの患者データを受信する。さらにまた、パルス発生器モジュ ール121は、患者およびマグネット装置の状態に連動する様々なセンサからの信 号を受け取るスキャンルームインターフェース回路133に接続している。また、 スキャンルームインターフェース回路133を通して、望ましいスキャン位置に患 者を動かす命令を、患者位置決め装置134が受け取る。 パルス発生器モジュール121によって生成された傾斜波形が、Gx,Gyおよび Gzの各アンプから構成される傾斜アンプ装置127に与えられる。各々の傾斜アン プは、位置エンコードの収集信号のために用いられる傾斜磁場を生成するために 、一般に139で示された組立体の対応する傾斜コイルを励起する。傾斜コイル組 立体139は、極性マグネット140と全身RFコイル152を含むマグネット組立体141 の部分を形成する。装置制御部122内の変換器モジュール150は、RFアンプ151 により増幅され、送信/受信スイッチ154によりRFコイル152に連結されるパル スを生成する。患者内で励起された原子核により放射された結果の信号が、同じ RFコイル152によって感知され、送信/受信スイッチ154を通してプリアンプ15 3に連結されてもよい。増幅されたNMR信号は、変換器150の受信部により復調 され、濾波され、デジタル化される。送信モード中ではRFアンプ151をコイル1 52に電気的に接続し、受信モード中ではプリアンプ153を接続するために、送信 /受信スイッチ154はパルス発生器モジュール121からの信号により制御される。 また、送信/受信スイッチ154により、別個のRFコイル(例えば、ヘッドコイ ルや表面コイル)を、送信モード或いは受信モードのいずれか一方に用いること が可能である。 RFコイル152により取り上げられたNMR信号は、変換器モジュール150によ りデジタル化され、装置制御部122のメモリモジュール160に転送される。スキャ ンが完全で、完全なデータ配列メモリモジュール160に収集されているときに、 アレイプロセッサ161は画像データ配列に対してのデータをフーリエ変換処理す る。この画像データは、シリアルリンク115を通してコンピュータ装置107に伝達 され、そこでディスクメモリ111に記憶される。オペレータコンソール100から受 け取った命令に応答して、この画像データをテープドライブ112に保管させても よく、さもなければ、画像プロセッサ106によって画像データを更に生成し、オ ペレータコンソール100に伝達して、表示器104上に表示してもよい。 特に、図1および図2を参照すると、変換器150はパワーアンプ151を通してコ イル152AでRF励起磁場B1を生成し、その結果としてコイル152Bに誘起され た信号を受け取る。先に指摘したように、各コイル152A,152Bは、図2に示す ように別個としてもよく、さもなければ、図1に示すような単一の全身コイルと してもよい。RF励起磁場の基底すなわちキャリア周波数は、CPUモジュール 119およびパルス発生器モジュール121からの一組のデジタル信号を受信する周波 数シンセサイザ200の制御下で生成される。これらのデジタル信号は、出力201に てRFキャリア信号の周波数と位相を表わす。命令されたRFキャリアは、変調 器およびアップコンバータ202に供給され、そこでRFキャリアの振幅が、パル ス発生器121からも受け取られる信号R(t)に応答して変調される。この信号R( t)は、生成されるべきRF励起パルスの包絡線を定めるものであって、一続き の記憶された信号値を連続的に読み出すことによって、モジュール121に生成さ れる。これらの記憶されたデジタル値は、すなわち、オペレータコンソール100 から変換され、あらゆる所望のRFパルス包絡線を生成することが可能になる。 出力205で生成されたRF励起パルスの振幅は、背板118からデジタル命令を受 け取る励起減衰器206により減衰される。減衰されたRF励起パルスは、RFコ イル152Aを駆動するパワーアンプ151に供給される。変換器122のこの箇所にお けるさらに詳細な記述は、参考としてその中に具体化されている米国特許第4,95 2,877号が参照となる。 図1および図2をさらに参照すると、対象物により生成された信号は、受信コ イル152Bで取り上げられ、プリアンプ153を通して受信減衰器207の入力に供給 される。受信減衰器207はさらに、背板118から受け取ったデジタル減衰信号によ り、決められた量で信号を増幅する。 受け取った信号は、ラーモア周波数、もしくはラーモア周波数の周辺にあり、 この高周波信号は、第1にNMR信号をライン201上のキャリア信号と混合し、 その後で結果として生じる差信号を、ライン204上の2.5MHzの基準信号と混合 するダウンコンバータ208により、2段階の過程にて低い周波数に変換される。 低い周波数に変換されたNMR信号は、アナログ−デジタル(A/D)コンバー タ209の入力に供給されるとともに、このアナログ−デジタル(A/D)コンバ ータ209は、アナログ信号をサンプルし、デジタル化するとともに、16ビットの 同相(I)値と、受信信号に対応する16ビットの直角位相(Q)値とを生成するデジ タル検出器および信号プロセッサ210に信号を供給する。その結果として生じる 受信信号のデジタル化されたI値およびQ値の流れは、背板118を通してメモリ モジュール160に出力され、そこで各値が画像を再構築するのに利用される。 250KHzのサンプリング信号や、5,10および60MHzの基準信号と同じよ うに、2.5MHzの基準信号は、基準周波数発生器203により共通の20MHzの原 クロック信号から生成される。受信器のより詳細な記述は、米国特許第4,992,73 6号公報が参照となる。 本発明は、複数の様々なパルスシーケンスで利用されることがあるものの、本 発明の好ましい実施態様では、図3に表わされる3D(三次元)傾斜呼び戻しエ コーパルスシーケンスを利用する。ゼネラル エレクトリック社の1.5テスラM Rスキャナを利用したパルスシーケンス“3dfgre”は、再バージョンレベル5.5 の装置を用いて、商標名“SIGNA”の下で販売された。それは、以下に述べ るようなk空間のサンプリングパターンを行なうことができるように、様々な容 積からのデータを補正するのに部分修正された。 図3を特に参照すると、米国特許第4,431,968号に教示されるように、関心の ある3D容積中に横方向磁化を生成するために、60°のフリップ角を有するRF 励起パルスが、スラブ選択傾斜パルス222の存在するところで生成される。Z軸 に沿って向けられた位相エンコード傾斜パルス224と、Y軸に沿って向けられた 位相エンコード傾斜パルス226が、これに続いて起こる。X軸に沿って向けられ た読出し位相エンコード傾斜パルス228がこの後に続くと、部分的なエコー(60 %)NMR信号230が収集され、この信号が上述したようにデジタル化される。 この収集の後で、米国特許第4,665,365号に教示されるように、パルスシーケン スが繰り返される前に、巻き戻りの傾斜パルス232,234が磁化を繰り返し段階的 に行なう。 公知技術のように、図4で表わされた3Dのk空間をサンプル化するために、 パルスシーケンスは繰り返されるとともに、連続する値を通して位相エンコード パルス224,226が段階的に行なわれる。好ましい実施態様によれば、16の位相エ ンコードがZ軸に沿って用いられ、128の位相エンコードがY軸に沿って用いら れる。したがって、特にY軸の各位相エンコードに対して、16の異なるZ軸位相 エンコードを有する16の収集が、kz軸に沿ってサンプル化を完全に行なうのに 実行される。これは、ky軸に沿ってデジタル化を完全に行なうのに、128の異な るZ軸の位相エンコードで128回繰り返される。以下の議論から明らかになるよ うに、このサンプリングが行なわれる順序は重要である。 kx軸に沿ったサンプリングは、各々のパルスシーケンス間の読出し傾斜パル ス228が存在するところで、エコー信号230をサンプリングすることで行なわれる 。kx軸に沿った部分的なサンプリングだけが行なわれて、ホモダイン再構築を 利用するか、若しくはゼロフィルを行なうことにより、失ったデータが算出され ることは、これにより、パルスシーケンスのエコータイム(TE)を、1.8から2.0 ms未満に短縮でき、パルスの反復率(TR)を10.0ms未満に短縮できる。 特に図4を参照すると、好ましい実施態様に従って動的な検討を行なうのに、 サンプルされるべきk空間が、“A”から“D”で示される4つの領域に分割さ れる。これらの領域“A”から“D”の境界線は、ky軸に沿って配置されると ともに、ky=0に対して対称である。中央の領域である“A”は、ky=−16か らky=+15に配列している中央のk空間の領域を占有し、公知技術のように、 これらの各「中央」サンプルが、再構築される画像の全体的なコントラストを決 定する大部分の情報を包含する。ここで記述されるように、動的な検討中におけ る各々のフレーム画像の基礎を形成し、結果として生じる時間的なフレーム率を 決定するのが、この中央のk空間領域である。 残っている3つの「周辺」にあるk空間領域B〜Dは、中央領域Aに対向する 各側に分割配置される。これらの領域は、次の範囲でk空間を占有する。 領域Bは、ky=−17からky=−32と、ky=+16からky=+31。 領域Cは、ky=−33からky=−48と、ky=+32からky=+47。 領域Dは、ky=−49からky=−64と、ky=+48からky=+63。 動的な検討中において、k空間の中央領域は周辺領域よりも、高率でサンプル 化される。好ましい実施態様では、中央領域Aと連続する周辺領域B〜Dの一つ を交互にサンプリングすることで、これが達成される。従って、次の各サンプリ ングシーケンスのどちらかが、動的な検討中に行なわれる。 AB AC AD AB AC AD AB AC AD AD AC AB AD AC AB AD AC AB 後者のサンプリングシーケンスは、図5にてグラフ的に示されているが、ここ では、水平軸が動的検討中における実時間を示しており、垂直軸がky軸に沿っ たk空間内のサンプルされるべき領域となっている。各々のk空間領域A〜Dが サンプル化される間の時間周期が標識化され、下つき文字を表示した数の領域が 、動的な検討中にサンプルされている。ここでは、中央のk空間領域Aが周辺の k空間領域B〜Dよりも高い時間率でサンプルされているのが容易にわかる。 本実施態様では、以下に議論される他の再構築計画を実行するために、領域A 〜Dの全てが、動的な検討の始めと再度最後にスキャンされる。上述した交互の シーケンスが、240で表わされた動的検討の重大な期間中にその後行なわれるこ とは、見分けることができよう。特別な検討中に現れた関心のあるコントラスト 変化を取り囲む必要がある限り、この交互のシーケンスを拡張することができる 。 本発明を実行するのに、他の方法でk空間を分割できることは、当業者により 見分けることができよう。例えば、領域の数は変更することができ、また、各領 域の境界線がスライス選択軸kzに沿って配置されるように、各領域を向けるこ とが可能である。さらに、k空間を円形の中央領域と、これを取り囲む環状の周 囲領域とに分割することもできる。 動的な検討の間に収集されたデータは、この動的な検討中に現れるコントラス トの変化を克明に表現した連続するフレーム画像F1〜Fnを再構築するのに、複 数の方法で利用されることがある。図6に示した一つの実施態様では、F1から F7で示された画像フレームは、各々の中央k空間領域で収集したもの(A1〜A7 )からのデータを用いて再構築される。これは、周りを取り囲む周辺のk空間 領域B〜Dからの時間的に近接したデータと組み合わせた特別な中央のk空間領 域データを用いて、フレーム画像を再構築するのに十分なデータセットがあれば 、このデータセットを形成することで完成される。各々の画像フレームデータの セットは、動的な検討中の特別な時間で、対象物を克明に表現する。 こうした画像フレームデータセットを形成する一つの方法は、中央のk空間領 域Aを収集する時間に最も近い、周辺領域から収集したデータを利用することで ある。従って、F2からF6の各フレーム画像に対し、図6にて収集され克明に表 現されたデータは、次のような各データセットに形成されてもよい。 F2 → A2+B2+C2+D2 3 → A3+B2+C2+(D2またはD3) F4 → A4+B2+(C2またはC3)+D3 5 → A5+(B2またはB3)+C3+D3 6 → A6+B3+C3+(D3またはD4) 画像フレームのデータに最も近い時間のデータを選択する方法は、“nearest neighbor(ニアレスト ネイバー)”法として、ここに引用される。k空間の周 辺領域に対し最も近いデータが、時にはフレーム時間に近く、また他の場合には 、フレーム時間が2つのサンプル期間の間の途中にある。 各々のフレームF2からF6で、データセットを形成する他の方法は、各々の周 辺領域に対して収集されたデータの近接する2セット間で、改変を行うことであ る。図6にて収集され克明に表現されたデータから、フレーム画像F2からF6を 形成する改変方法の例は、次のようになる。 F2 → A2+(B1+B2)/2+(4C2+C1)/5+(5D2+D3)/6 F3 → A3+(5B2+B1)/6+(5C2+C3)/6+(D2+D3)/2 F4 → A4+(5B2+B3)/6+(C2+C3)/2+(5D3+D2)/6 F5 → A5+(B2+B3)/2+(5C3+C2)/6+(5D3+D4)/6 F6 → A6+(5B3+B2)/6+(4C3+C4)/5+(D3+D4)/2 上述の係属中の米国特許出願第622,718号に開示されるように、領域B,C, Dと同様に各領域A間を改変することにより、中間のフレーム画像(例えば、F2.5 ,F3.5,F4.5およびF5.5)を生成することも可能である。 上で表したように、各々のフレームでのデータセットを形成するのに用いられ る各ビューは、関心のある領域(ROI)にコントラスト仲介物が流れ込む間の 時間中に亘って収集される。その結果、コントラスト仲介物の量、つまり信号の 強さは、各々のデータセットを形成するのに利用されるビューの収集の間に変更 されるかもしれない。こうした信号強度の変動が、結果として正確な診断を阻害 する再構築された画像による画像アーティファクトを生み出すことがわかってい る。このようなアーティファクトの一つが、腎臓動脈のちょうど下にある大動脈 の再構築画像全体の断面図として、図7のグラフに示されている。280で示す強 度の落ち込みは、動脈中央の陰影領域として現れたアーティファクトである。本 発明の一つの見地は、こうした画像アーティファクトを抑制して、収集されたN MRデータを補正することにある。 三次元フーリエ変換法を用いて、三次元フレーム画像の対応するセットを再構 築するのに、画像フレームデータが用いられてもよい。こうした6つの二次元フ レーム画像が、フレーム画像データ250〜255として図8に示されているが、実際 には30個のこうしたデータセットが、通常の動的な検討中に生成される。こうし たスキャンから三次元画像を再構築するには膨大な処理時間を必要とする。本発 明の他の見地は、この処理時間を減らすことにある。 図8を参照すると、これは時間フレーム250〜255の各々における三次元k空間 データの部分体から、二次元画像を生成することにより達成される。例えば、各 々の三次元k空間データセット250〜255と、対応する各画像260〜265を生成する のにそこで行われる二次元フーリエ変換とから、一般には中央のスライスである 一つのk空間データの二次元スライスを選択してもよい。結果として生じるプレ ビューの画像が平均化され、画質を改良することができるように、他のk空間を 用いてこの処理を任意に繰り返してもよい。こうした二次元の再構築は、画像サ イズと使用されるコンピュータの各々に依存して、約0.5秒から5秒を必要とし 、オペレータは、データが収集されるのとほぼ同じ速さの、(通常は)30の時間 的フレームすべてを表した2次元画像を、持つことができる。好ましくは、最初 の二次元画像フレーム260は、デフォルトの参照フレームとして用いられ、後に 続く二次元画像フレーム261〜265が、連続する二次元の消去画像を形成するのに 生成されるときに、二次元の画像フレーム261〜265から差し引かれる。他のもの から二次元消去画像のどれかを差し引くことにより、どちらか一つの参照フレー ムが単に選ばれてもよい。例えば、二次元消去フレーム画像(261〜260)が、二 次元消去画像(263〜260)から差し引かれると、その結果は二次元画像フレーム 261が参照となる二次元消去画像(例えば、263〜261)となる。これらの二次元 画像260〜265あるいは関連する各消去画像は、以下に述べるような2つの目的に 対し利用されることがある。 第1に、二次元画像260〜265の各々を再検討するとともに、関心のある血管を 最も克明に表わした1乃至2の二次元画像を選択してもよい。例えば図8を参照 すると、仮にコントラスト仲介物が到達した後の、コントラスト仲介物による強 化の頂点で、二次元画像261および二次元画像264が血管を克明に表現しているな らば、対応する三次元のフレームデータセット251,254が、完全なる処理のため に選択される。k空間差分データセットが形成され、三次元差分画像データセッ トを得るのに、このk空間差分データセットが三次元フーリエ変換されている間 に、好ましい実施態様では、k空間のデータセット251,254の各々が、個々にフ ーリエ変換され、三次元差分画像データセット270を形成するのに、画像空間内 に差し引かれる。この2つの最善のデータセットによる処理は、30のデータセッ ト全てを再構築するのに比べて15の要因だけでよく、全体的な処理時間の減少を 生み出す。普通の価格のコンピュータを利用すると、30のデータセットを再構築 するのに必要な時間は、画像サイズや使用する受信コイルの数に依存して、5分 から1時間である。上述の提案では、再構築の時間が1分から4分未満に減少す る。仮に再構築のために追加のデータセットが選択されるならば、これらの時間 見積りは増加する。 全体的な血管の構造および状態を評価するために、この三次元画像データの三 次元配列270を、単一の二次元投影画像272に投影するのが、通常はさらに役に立 つ。これを行うために最も一般的に利用される技術は、投影画像272内の各ピク セル(画素)からの光線を三次元配列270の各画像データ点に通過投影するとと もに、その中から最大値を有するデータ点を選択することである。各々の光線に 対し選択された最大値は、投影画像272に対応するピクセルの輝度を制御するの に利用される。この方法は「最大値ピクセル技術」として以下に参照されるが、 実行するのが非常に簡単で、感覚的に満足する画像が与えられる。目下これは、 好ましい方法である。 投影画像272を形成するのに利用され、利用できる情報のより多くを保有する 別の技術は、「インテグレーション法」のような以下のものを参照することであ る。この投影法は、「適応性のあるNMR血管造影法における投影法」と題する 米国特許第5,204,627号公報に記述されている。この方法によれば、各々の投影 画像ピクセルの輝度が、全てのデータ点の合計により投影光線に沿って決定され る。 投影画像272を生成するのに利用されるさらに別の技術は、三次元領域成長法 を用いることである。成長されるべき三次元の画像データセット内における領域 の起点を、オペレータが決定する。成長した領域はその後ぼんやりとし、血管端 部のちょうど外側のボクセルを含むマスクを作るのに閾値化されるが、このボク セルは、領域が成長する過程において削られることになる。この方法は、血管の 端部が保持されており、しかも、描出過程に含まれる視覚的な指示を利用するこ とにより、血管の重腹部分を導き出すことが可能な、非常にスムーズな血管の表 現を提供する。 上で表したように、各々のフレームでのデータセットを形成するのに用いられ る各ビューは、関心のある領域(ROI)にコントラスト仲介物が流れ込む間の 時間中に亘って収集される。その結果、コントラスト仲介物の量、つまり信号の 強さは、各々のデータセットを形成するのに利用されるビューの収集の間に変更 されるかもしれない。こうした信号の重み付けの変動は、結果的に再構築された 画像の画像アーティファクトに生じることが議論されている。これらのアーティ ファクトは、正確な診断を阻害する。こうした一つのアーティファクトが図7の グラフに表わされているが、これは、人造動脈のちょうど下にある大動脈の再構 築画像全体の断面図である。280で示す強度の落ち込みは、動脈中央の陰影領域 として現れたアーティファクトである。この問題の解決は、NMR信号のコント ラスト増大に依存した時間に関する情報を収集すること、および、動的な検討中 に収集されるNMR信号の再重み付けすなわち正規化するのに、その情報を利用 することである。 特に図9を参照すると、関心のある領域にコントラスト仲介物が入り込んだと きのNMR信号の重み変化を表わす一般的なコントラスト曲線は、コントラスト 仲介物が到達する前に測定された、あるバックグラウンド値(B)を有する。N MRのk空間データ(Ik)は、動脈血流(Iart)による要素と、バックグラウ ンド(Ibkg)による要素とを含んでいる。動脈血流による要素(Iart)は、基 線値Bより上方の曲線(C(t))で、重みの増加を得る。バックグラウンド要 素(Ibkg)が、コントラスト仲介物が最初に通過する間で、一定の重み付けを 有するものと仮定すると、 Ik(t)=B*Ibkg+[B+C(t)]*Iart (1) 本発明に基づく再重み付けが、バックグラウンドのゴースティング(ghostin g)に起因するものでないことを保証するために、主要なコントラストk空間か ら、コントラスト仲介物が到達する前に入手したk空間の基線データを差し引く ことで形成された消去k空間データΔIkで、全ての再重み付けが行なわれてい る。これは、次のようにして与えられる。 ΔIk=Ik(t)−Ik(0)=C(t)*Iart (2) 我々にとって、変化する重み付けを差し引くためにC(t)で割ることは、零 で割る可能性があることから実際にはできない。一つの選択としては、例えば、 0.1などの比をfとして、C(t)+f*Bで割ることがある。この割られたも のは、次のようにして与えられる新規な連続するk空間差分データΔIk’を生 成する。 ΔIk’=ΔIk/[f*B+C(t)] (3) 小さなfによって、これは図10に示すように、コントラスト曲線の開始部を 鋭くし、また、等しい重み付けを有するk空間データの高い部分を生成する。再 重み付けされたビューが平坦な高い部分284に収集されるならば、この曲線282の 鋭さは、ゴースティングのない一様な画像を導き出す。 全てのアーティファクトを排除するために、コントラスト仲介物が関心のある 領域に到達し、重み付け要因の高い部分が確立された後にのみ、k空間の各ビュ ーが収集されるべきである。一つの解決法は、コントラスト仲介物を注入するす ぐ前あるいは注入したすぐ後に、ABCDのシーケンスを行ない、その後で、コ ントラスト仲介物が到達するまで、Aでの収集(ここでいうAは、上述した中央 のk空間の収集を引用している)を継続し、その地点で、ABACADのシーケ ンスが開始されるようにすることである。その後、ニアレスト ネイバー法は、 動的な検討における最初の三次元フレームデータセットに対し、後に続くB,C およびDの各k空間で収集したものを見つけるために、最初のAの収集したもの から時間的に後を探す。その後で、各Aで収集したものを時間的に平均化し、外 側にあるk空間位置の最も近くにある隣のものを選択することが、上述したよう に行える。A9を収集する間にコントラスト仲介物が到達したと仮定すると、シ ーケンスはこのようになる。 A1,B2,C3,D4,A5,A6,A7,A8,A9,B10,A11,C12 , A13,D14,A15,B16,A17,C18,A19,C20... 基線は、k空間のビューA1+B2+C3+D4から形成されるとともに、次 のk空間から差し引かれる。 A9+B10+C12+D14 B10+(A9+A11)/2+C12+D14 A11+B10+C12+D14 このK空間データは、その後、重み付け信号C(t)を利用した等式(3)に 基づき補正される。一つの実施態様では、画像データが収集されるときに重み付 け信号C(t)が生成され、この画像データの正規化をすぐに行うことができる 。しかし、次の議論で指摘されるように、他の実施態様では、k空間データが収 集されているときに、重み付け信号C(t)がわからなくてもよい。従って、k 空間データが収集されているときのデータを、時間的に印付けすることが必要で ある。これは、k空間データが収集された瞬間の重み付け信号C(t)の値が正 確にわかっているときに、等式(3)による補正が後で行なわれることを可能に する。 重み付け信号C(t)に対する各サンプル値Ciを生成するのに、多くの方法 を利用してよい。例えば、上述したゼネラル エレクトリック社のMRI装置で 利用できる「MR スマートプレップ“MR Smartprep”」を用いて、各サンプル 値Ciを生成してもよい。この方法は、上述した画像データの収集と交互に行な われるパルスシーケンスを実行することによって、関心のある領域からのNMR 信号強度を監視する。仮に、このパルスシーケンスが動的な検討の間に頻繁に十 分実行されるならば、画像データが収集されるときに、重み付け信号C(t)の 最新値を維持できる。その後、すぐに正規化を実行できる。重み付け信号C(t )を集めるのに利用されるパルスシーケンスは、例えば、コントラスト仲介物が 到達したことを示すために、腹部の大動脈あるいは他の動脈に対するスライス垂 線を定めるスライス選択RFパルスを有する標準の傾斜エコーシーケンスとする ことがある。収集された信号のフーリエ変換を行なうというよりも、むしろ重み 付け信号C(t)を表わすのに、統合された信号を用いてよい。特別な血管内に 現れるコントラスト信号をより正確に定めるために、関心のある血管を通過 する円筒励起を定めるRFパルスをさらに用いてもよい。 代わりに、収集された画像データから、重み付け信号C(t)の値を引き出す ようにしてもよい。例えば、動的な検討中における周期的な値Ciを生成するの に、上述したような再構築される二次元画像を利用してもよい。この場合、オペ レータは最初の二次元画像の関心のある領域を確認するとともに、基線値Bを形 成するのに、その中で強度の値が統合される。同じ強度の値が、後に続く二次元 画像に統合される。その後、図9に示すように、コントラスト増大値Ciを生成 するのに、各々のフレームにおける関心のある領域の各信号から、基線値Bが差 し引かれる。 この方法の利点は、放射線学者に関心のある血管内の正確な位置で、画像コン トラストを監視することを可能にし、それによって関心のある領域に、最善のア ーティファクト抑制が現れてくることにある。この方法の不利な点は、画像デー タが収集されるときに、C(t)の値がわからないことと、正規化の段階が後ろ 向きに行なわれなければならないことにある。 コントラスト増大値Ciを生成するのに利用される方法に関わらず、これらの 増大値は、動的な検討中に分離した瞬間で、単にコントラストの増大を表わすだ けである。コントラスト増大値Ciを生成するのに二次元の各画像が利用される ならば、例えば、各々の収集された画像フレームに対し、一つの値Ciが生成さ れる。より明確には、中央のk空間ビュー(すなわち、上述のシーケンスA)が 、各々の画像フレームに対し収集される瞬間で、各々の値Ciがコントラストの 増大を表わす。他方で「MR スマートプレップ」法が利用されるならば、一つ の画像フレームにつき1回よりも頻繁に、そのパルスシーケンスが交互に重なり 合って実行されるが、スキャン時間を増加することなく、しかも、動的な検討の 時間的解決を減らすことなくして、それ以上に余分な頻度でパルスシーケンスを 実行することはできない。従って、コントラスト増大曲線C(t)を生成するの に、補間法を利用したコントラスト増大値Ciにスムーズな曲線が適合すること が、本発明の教示となっている。 C(t)=Interp(Ci) ここで、“Interp”は補間を表わしている。 各々の収集されたNMR信号と結合した「消印」は、その後、対応する時間( t)でコントラスト増大値(C)を検索するのに利用されてもよい。
───────────────────────────────────────────────────── フロントページの続き (81)指定国 EP(AT,BE,CH,DE, DK,ES,FI,FR,GB,GR,IE,IT,L U,MC,NL,PT,SE),OA(BF,BJ,CF ,CG,CI,CM,GA,GN,ML,MR,NE, SN,TD,TG),AP(KE,LS,MW,SD,S Z,UG),UA(AM,AZ,BY,KG,KZ,MD ,RU,TJ,TM),AL,AM,AT,AU,AZ ,BB,BG,BR,BY,CA,CH,CN,CZ, DE,DK,EE,ES,FI,GB,GE,HU,I S,JP,KE,KG,KP,KR,KZ,LK,LR ,LS,LT,LU,LV,MD,MG,MK,MN, MW,MX,NO,NZ,PL,PT,RO,RU,S D,SE,SG,SI,SK,TJ,TM,TR,TT ,UA,UG,UZ,VN (72)発明者 グリスト,トーマス,エム. アメリカ合衆国,ウイスコンシン州 53711,マディソン,カウンシル クレス ト 3805番地 (72)発明者 フライン,リチャード アメリカ合衆国,ウイスコンシン州 53705,マディソン,ヒルクレスト ドラ イヴ 3712番地 (72)発明者 コロセック,フランク アメリカ合衆国,ウイスコンシン州 53711,マディソン,リチャードソン ス トリート 2759番地

Claims (1)

  1. 【特許請求の範囲】 1.a)対象物にコントラスト仲介物を注入し、 b)NMRパルスシーケンスを実行し、動的な検討の時間区分の間に、選択 された三次元k空間からの各NMRサンプル信号を収集して、三次元データセッ トを形成し、 c)動的な検討の間の連続する各時間区分に前記b)の手順を繰り返し行な い、対応する複数の追加の三次元データセットを形成し、 d)各三次元データセットの部分体を選択し、 e)この選択されたデータの部分体を利用して、各々の三次元データセット から二次元画像を再構築し、動的な検討の間の対象物を表現する連続する二次元 画像を形成し、 f)前記各三次元データセットの一乃至それ以上を選択し、 g)各々の選択された三次元データセットから、三次元フレーム画像を再構 築する各手順からなる、 動的な検討の間に収集されたNMRデータから、対象物のコントラスト増大 磁気共鳴画像を生成する方法。 2.h)三次元フレーム画像を投影することにより、二次元画像を生成する手順 をさらに含む請求項1記載の方法。 3.前記三次元データセットの一つが第2の三次元データセットとして選択され 、この第2の三次元データセットから第2の三次元フレーム画像が再構築される とともに、前記2つの三次元フレーム画像から三次元差分フレーム画像が生成さ れる請求項1記載の方法。 4.h)三次元差分フレーム画像を投影することにより、二次元画像を生成する 手順をさらに含む請求項3記載の方法。 5.手順d)で選択されたデータの部分体は、前記三次元空間内の単一の二次元 平面からのサンプル信号を含む請求項1記載の方法。 6.前記データの部分体は、前記三次元k空間の中央に近くにあるサンプル信号 にさらに限定される請求項5記載の方法。 7.h)再構築された二次元画像の関心のある領域を選択し、前記関心のある領 域内の画像データからコントラストの増大を表わす信号を生成する手順をさらに 含む請求項1記載の方法。 8.a)対象物にコントラスト仲介物を注入し、 b)NMRパルスシーケンスを実行し、動的な検討の時間区分の間に、選択 された三次元k空間からの各NMRサンプル信号を収集して、三次元データセッ トを形成し、 c)動的な検討の間の連続する各時間区分に前記b)の手順を繰り返し行な い、対応する複数の追加の三次元データセットを形成し、 d)動的な検討の間の連続する時間間隔で、コントラスト仲介物による各N MRサンプル信号の増大した振幅を表わすコントラスト増大値を生成し、 e)各NMRサンプル信号が収集された時間に対応して、この時間からのコ ントラスト増大値を利用して、各三次元データセットの一つのNMRサンプル信 号を正規化し、 f)正規化されたサンプル信号を利用して画像を再構築する各手順からなる 、動的な検討の間に収集されたNMRデータから、対象物のコントラスト増大磁 気共鳴画像を生成する方法。 9.NMRパルスシーケンスを実行し、対象物の関心のある領域にコントラスト 仲介物が入り込む前に、前記選択された三次元k空間からサンプル信号を収集し て、基線のデータセットを形成するとともに、前記正規化の手順e)は、三次元 データセットの一つの対応するNMRサンプル信号から、基線のデータセットを 差し引いて、差分データセットを形成することを含む請求項8記載の方法。 10.前記正規化の手順e)は、前記差分データセットを前記コントラスト増大値 で割る手順をさらに含む請求項9記載の方法。 11.対象物の関心のある領域からのNMR信号を収集し、このNMR信号を利用 して、前記コントラスト増大値を生成するために、各コントラスト増大値は、動 的な検討の間にNMRパルスシーケンスを実行することにより、手順d)で生成 される請求項8記載の方法。 12.各コントラスト増大値は、前記三次元データセットの各々からなる部分体を 処理することで、手順d)で生成される請求項8記載の方法。 13.前記部分体は、前記三次元k空間内の二次元平面からの各NMRサンプル値 であるとともに、前記部分体の処理は、該部分体からの二次元画像を再構築する ことを含む請求項12記載の方法。 14.前記二次元画像の関心のある選択された領域における各値から、各コントラ スト増大値が計算される請求項13記載の方法。 15.連続する時間間隔で各コントラスト増大値が生成される間を補間することに よりNMRサンプル信号が収集される各時間に対応して、この各時間に各コント ラスト増大値が計算される請求項8記載の方法。
JP10508762A 1996-07-26 1997-02-03 対象物のコントラスト増大磁気共鳴画像を生成するmri装置 Expired - Lifetime JP3130320B2 (ja)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US08/686,795 US5881728A (en) 1996-07-26 1996-07-26 Digital subtraction magnetic resonance angiography with image artifact suppression
US686,795 1996-07-26
US08/686,795 1996-07-26
PCT/US1997/001781 WO1998004928A1 (en) 1996-07-26 1997-02-03 Digital subtraction magnetic resonance angiography with image artifact suppression

Publications (2)

Publication Number Publication Date
JPH10512482A true JPH10512482A (ja) 1998-12-02
JP3130320B2 JP3130320B2 (ja) 2001-01-31

Family

ID=24757798

Family Applications (1)

Application Number Title Priority Date Filing Date
JP10508762A Expired - Lifetime JP3130320B2 (ja) 1996-07-26 1997-02-03 対象物のコントラスト増大磁気共鳴画像を生成するmri装置

Country Status (6)

Country Link
US (1) US5881728A (ja)
EP (1) EP0852731B1 (ja)
JP (1) JP3130320B2 (ja)
AU (1) AU2257497A (ja)
DE (1) DE69735748T2 (ja)
WO (1) WO1998004928A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002045584A1 (fr) * 2000-12-04 2002-06-13 Hitachi Medical Corporation Procede de mesure dans un dispositif d'imagerie par resonance magnetique et dispositif d'imagerie par resonance magnetique
JP2006122301A (ja) * 2004-10-28 2006-05-18 Ge Medical Systems Global Technology Co Llc Mri装置
JP2006314471A (ja) * 2005-05-11 2006-11-24 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2015526227A (ja) * 2012-08-27 2015-09-10 コーニンクレッカ フィリップス エヌ ヴェ ダイナミック造影mri用の適応するキーホール圧縮

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6043654A (en) * 1997-11-14 2000-03-28 Picker International, Inc. Multi-volume slicing and interleaved phase-encoding acquisition for 3 D fast spin echo (FSE)
US6195579B1 (en) * 1998-12-18 2001-02-27 Wisconsin Alumni Research Foundation Contrast detection and guided reconstruction in contrast-enhanced magnetic resonance angiography
US6259940B1 (en) * 1999-04-28 2001-07-10 Mayo Foundation For Medical Education And Research Method of performing magnetic resonance angiography using two-dimensional imaging and de-rated gradients
WO2000065995A1 (fr) * 1999-04-30 2000-11-09 Hitachi Medical Corporation Procede d'imagerie par resonance magnetique et dispositif connexe
DE60132687T2 (de) 2000-03-27 2009-01-29 Koninklijke Philips Electronics N.V. Kernspin-verfahren zur bilderzeugung eines zeitabhängigen kontrastes
EP1204879A1 (en) * 2000-07-31 2002-05-15 Koninklijke Philips Electronics N.V. Magnetic resonance method and apparatus
US6725077B1 (en) * 2000-12-29 2004-04-20 Ge Medical Systems Global Technology Company, Llc Apparatus and method for just-in-time localization image acquisition
US20030060731A1 (en) * 2001-01-26 2003-03-27 Fleischhacker Mark G. Non-metallic guide wire
DE10119660B4 (de) * 2001-04-20 2006-01-05 Siemens Ag Verfahren zur schnellen Gewinnung eines Magnetresonanzbildes
US20030065490A1 (en) * 2001-09-27 2003-04-03 Lowry David F. System and technique for interrogation of an object
DE10147919A1 (de) * 2001-09-28 2003-04-30 Siemens Ag Verfahren zur Erzeugung von Bilddaten mittels magnetischer Resonanz
US6728569B2 (en) 2001-10-25 2004-04-27 Evanston Northwestern Healthcare Corp. Scoutless whole-body imaging with fast positioning
DE10161925B4 (de) * 2001-12-17 2005-08-18 Siemens Ag Verfahren zum Aufnehmen von Magnetresonanzsignalen eines Objekts mittels einer Magnetresonanzanlage und Rekonstruieren eines Bildes anhand der aufgenommenen Magnetresonanzsignale, sowie Computerprogrammprodukt und Magnetresonanzanlage
US6903551B2 (en) * 2002-05-01 2005-06-07 Brigham & Women's Hospital, Inc. Variable-density parallel magnetic resonance imaging
EP1601987A1 (en) * 2003-03-03 2005-12-07 Koninklijke Philips Electronics N.V. Magnetic resonance imaging method
DE10353342B4 (de) * 2003-11-14 2008-07-17 Siemens Ag Verbesserte MRT-Bildgebung auf Basis konventioneller PPA-Rekonstruktionsverfahren
JP4434753B2 (ja) * 2004-01-13 2010-03-17 株式会社東芝 磁気共鳴イメージング装置および磁気共鳴イメージング収集データの処理方法
DE102004043809B4 (de) * 2004-09-08 2008-01-24 Charité-Universitätsmedizin Berlin Verfahren zum Kalibrieren einer kontrastmittelgestützten Perfusionsbildgebung
US20070055138A1 (en) * 2005-08-22 2007-03-08 Edelman Robert R Accelerated whole body imaging with spatially non-selective radio frequency pulses
JP4964450B2 (ja) * 2005-10-18 2012-06-27 譲二 中川原 画像処理方法、画像処理プログラム、及び画像処理装置
US20080150942A1 (en) * 2006-10-13 2008-06-26 Calgary Scientific Inc. Visualization of volumetric medical imaging data
US8214013B2 (en) * 2007-01-26 2012-07-03 Mayo Foundation For Medical Education And Research Accelerated shells trajectory MRI acquisition
CN101903911A (zh) * 2007-12-20 2010-12-01 皇家飞利浦电子股份有限公司 身体和身体轮廓的3d重建
US20100145185A1 (en) * 2008-12-04 2010-06-10 Xiaole Hong Extending the resolution of mri data by combining subsets from plural image acquisitions
US8653817B2 (en) 2010-04-02 2014-02-18 General Electric Company Accelerated pseudo-random data magnetic resonance imaging system and method
US8738113B2 (en) 2010-09-30 2014-05-27 University Of Utah Research Foundation Retrospectively correlated turbo spin echo imaging
DE102012107926A1 (de) * 2012-08-28 2014-05-28 Digital Medics Gmbh Verfahren und Vorrichtung zur Ermittlung einer Kontrastmittelanreicherung
US10670678B2 (en) 2015-06-15 2020-06-02 Koninklijke Philips N.V. MR imaging using stack-of stars acquisition
KR101744424B1 (ko) * 2015-12-24 2017-06-08 (의료)길의료재단 Mra 영상을 이용하여 관류 특성 확인용 추가 영상을 얻기 위한 방법 및 시스템
DE102016202663A1 (de) * 2016-02-22 2017-08-24 Siemens Healthcare Gmbh MRT-Vorschau
US9786069B2 (en) * 2016-03-07 2017-10-10 Siemens Healthcare Gmbh Refined reconstruction of time-varying data
CN107993271A (zh) * 2017-12-26 2018-05-04 上海交通大学 一种磁共振动态成像采样方法和图像重建方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4204226A (en) * 1978-05-16 1980-05-20 Wisconsin Alumni Research Foundation Real-time digital X-ray time interval difference imaging
US4204225A (en) * 1978-05-16 1980-05-20 Wisconsin Alumni Research Foundation Real-time digital X-ray subtraction imaging
DE3603579A1 (de) * 1986-02-06 1987-08-13 Henkel Kgaa Verwendung ethoxylierter fettamine als loesungsvermittler
JPH0634785B2 (ja) * 1988-02-05 1994-05-11 株式会社日立メディコ 核磁気共鳴イメージング装置
JPH02140145A (ja) * 1988-11-21 1990-05-29 Toshiba Corp 磁気共鳴イメージング装置
US5190744A (en) * 1990-03-09 1993-03-02 Salutar Methods for detecting blood perfusion variations by magnetic resonance imaging
US5348011A (en) * 1991-11-14 1994-09-20 Picker International, Inc. Shared excitation phase encode grouping for improved throughput cardiac gated MRI cine imaging
US5522390A (en) * 1991-11-21 1996-06-04 U.S. Philips Corporation Magnetic resonance imaging method
US5417213A (en) * 1993-06-07 1995-05-23 Prince; Martin R. Magnetic resonance arteriography with dynamic intravenous contrast agents
US5377680A (en) * 1993-08-04 1995-01-03 General Electric Company MRI cardiac image produced by temporal data sharing
US5451876A (en) * 1993-10-18 1995-09-19 General Electric Company MRI system with dynamic receiver gain
US5685305A (en) * 1994-08-05 1997-11-11 The United States Of America As Represented By The Department Of Health And Human Services Method and system for MRI detection of abnormal blood flow

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002045584A1 (fr) * 2000-12-04 2002-06-13 Hitachi Medical Corporation Procede de mesure dans un dispositif d'imagerie par resonance magnetique et dispositif d'imagerie par resonance magnetique
US7167740B2 (en) 2000-12-04 2007-01-23 Hitachi Medical Corporation Measuring method in magnetic resonance imaging device and magnetic resonance imaging device
JP2006122301A (ja) * 2004-10-28 2006-05-18 Ge Medical Systems Global Technology Co Llc Mri装置
JP4699739B2 (ja) * 2004-10-28 2011-06-15 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Mri装置
JP2006314471A (ja) * 2005-05-11 2006-11-24 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2015526227A (ja) * 2012-08-27 2015-09-10 コーニンクレッカ フィリップス エヌ ヴェ ダイナミック造影mri用の適応するキーホール圧縮

Also Published As

Publication number Publication date
DE69735748T2 (de) 2006-12-28
JP3130320B2 (ja) 2001-01-31
AU2257497A (en) 1998-02-20
WO1998004928A1 (en) 1998-02-05
EP0852731B1 (en) 2006-04-26
DE69735748D1 (de) 2006-06-01
US5881728A (en) 1999-03-16
EP0852731A1 (en) 1998-07-15

Similar Documents

Publication Publication Date Title
JPH10512482A (ja) 画像アーティファクトを抑制したデジタル消去磁気共鳴血管造影法
JP4011585B2 (ja) Mriシステム
EP1010016B1 (en) Gated time-resolved contrast-enhanced 3d mr angiography
US5873825A (en) Three dimensional digital subtraction magnetic resonance angiography with limited k-space mask
JP3742114B2 (ja) ゲート化され時間解像されたコントラスト増幅3次元mr血管造影法
JP3878176B2 (ja) インターリービングされた投影−再構成データを使用する三次元位相コントラスト磁気共鳴イメージング
JP4889903B2 (ja) 患者の磁気共鳴血管造影図を作成するためのmriシステムの作動方法
EP1145029B1 (en) Contrast guided reconstruction in contrast-enhanced magnetic resonance angiography
JP4484371B2 (ja) 血管セグメンテーションを伴う造影剤増幅型磁気共鳴血管造影法
JP3964202B2 (ja) Mrイメージングにおいて造影強調後の定常自由歳差運動を使用する方法及び装置
US6505064B1 (en) Diagnostic imaging systems and methods employing temporally resolved intensity tracing
JP4477176B2 (ja) Mrアンジオグラフィのための選択的な動脈画像及び静脈画像の取得の方法及びシステム
US6556856B1 (en) Dual resolution acquisition of magnetic resonance angiography data with vessel segmentation
JP2004524079A (ja) アンダーサンプリングされた3d投影イメージング法を用いる磁気共鳴血管造影法
JP2002534180A (ja) 角度付けインターリービングされた投影データを用いる位相コントラスト磁気共鳴血流イメージング
US6044290A (en) Time-resolved digital subtraction magnetic resonance angiography using echo-planar imaging
WO2009047690A2 (en) Segmented multi-shot mri involving magnetization preparation
US20020133070A1 (en) Method for performing magnetic resonance angiography with subtraction of projection images
WO2001075469A1 (en) Magnetic resonance angiography with automated vessel segmentation
WO2008041208A2 (en) Segmented multi-shot mri involving magnetization preparation

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20071117

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081117

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091117

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091117

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101117

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111117

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111117

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121117

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121117

Year of fee payment: 12

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131117

Year of fee payment: 13

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term