JPH057571A - 磁気共鳴影像方法及び装置 - Google Patents

磁気共鳴影像方法及び装置

Info

Publication number
JPH057571A
JPH057571A JP3335331A JP33533191A JPH057571A JP H057571 A JPH057571 A JP H057571A JP 3335331 A JP3335331 A JP 3335331A JP 33533191 A JP33533191 A JP 33533191A JP H057571 A JPH057571 A JP H057571A
Authority
JP
Japan
Prior art keywords
gradient
pulse
data
magnetic resonance
sequence
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.)
Pending
Application number
JP3335331A
Other languages
English (en)
Inventor
Johannes J M Cuppen
ヨセフス マリア クツペン ヨハネス
Miha Fuderer
フデラー ミハ
Antoon F Mehlkope
フラン メルコツプ アントン
Michael J Duijvestijn
ヨゼフ デユイブステイン ミカエル
Yperen Gerrit H Van
ヘンドリツク ヴアン イピーレン ゲリツト
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Philips Gloeilampenfabrieken NV
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 Philips Gloeilampenfabrieken NV filed Critical Philips Gloeilampenfabrieken NV
Publication of JPH057571A publication Critical patent/JPH057571A/ja
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (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)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Medical Treatment And Welfare Office Work (AREA)

Abstract

(57)【要約】 （修正有） 【目的】 あるクラスの影像エラー発生源、例えば信号
取得中外部磁界により安定磁界の影響する信号取得中の
対象の剛体運動と対象の変形・計器ドリフト現象による
磁気共鳴撮像中の影像エラー減少方法を提供する。 【構成】 第１のスピンエコーパルス・勾配シーケンス
ｐｓ１は対象に印加される。パルス・勾配シーケンスｐ
ｓ１は対象の核スピンを時点ｔ１で送信手段により送信
／受信コイルを介して対象に印加される励起パルスｅｐ
１より成り、またディフューズされた核磁化を再集束す
る時点ｔ２で印加される再集束パルスｒｐ１からなる。
プログラム手段を用いて核磁化分布はフーリエ変換によ
りデータセットから再構成され、分布は磁気共鳴影像の
形式で表示される。時点ｔ１′での励起パルスｅｐ２と
ｔ２′での再集束パルスｐｐ２と、選択勾配Ｇｚと勾配
Ｇｘ・Ｇｙとからなる第２の勾配シーケンスｐｓ２が印
加される。磁気共鳴影像が積分変換により補正される。

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明は、安定磁界に配置された
対象からの磁気共鳴信号から形成された磁気共鳴影像中
の影像エラーを減少させる磁気共鳴影像方法に係る。
【０００２】本発明は又、安定磁界に配置された対象か
らの磁気共鳴信号から形成された磁気共鳴影像の影像エ
ラーを減少させ、安定磁界を発生する手段と、ＲＦ電磁
パルスを発生する手段と、安定磁界に重畳された勾配磁
界用手段とからなる磁気共鳴影像装置に係る。
【０００３】
【従来の技術】前記の種類の磁気共鳴影像方法及び装置
は欧州特許出願第２１２５２６号で公知である。引用の
出願は、エラーが対象の周期的運動、例えば呼吸により
生じる磁気共鳴影像で影像エラーを減少させる方法を開
示する。磁気共鳴信号が、公知の方法、例えば欧州特許
第２１２５２６号に開示の如くいわゆるスピンワープシ
ーケンスによって、得られる場合、対象の周期的運動中
にゴーストエラー及びぼけのような影像エラーは共鳴信
号の信号サンプルから生じた磁気共鳴影像中に生じる。
欧州特許第２１２５２６号は、如何にかかる影像エラー
がパルス及び勾配シーケンスで、所定のダイヤグラムに
応じてパルス及び勾配シーケンスでの位相符号化勾配の
振幅を選択することにより減少されうるかを示す。従っ
て、欧州特許第２１２５２６号は受信した磁気共鳴信号
の周期的変化により本質的に生じる影像エラーを減少さ
せることを目的とする。しかし、影像エラーの他の発生
源は、例えば安定磁界を妨害する可変外部磁界に存在す
る。これらの外部磁界は、ＭＲＩ装置の近辺のエレベー
タを動かしたり、通る電車又は列車等が通り過ぎること
により生じうる。これは、安定磁界が変動することを意
味する。最後に、影像再構成の後、磁界妨害は磁気共鳴
影像中に影像エラーを生じる。呼吸中に生じる如く対象
の変形をしない対象の剛体運動は又影像エラーを生じ
る。更に、安定磁界が永久磁石により発生される場合、
ドリフト現象は温度変動により安定磁界中に生じる。送
信／受信チェーンの増幅器又は勾配増幅器の遅いドラフ
ト現象は同様の問題を生じうる。
【０００４】
【発明が解決しようとする課題】本発明の目的は、ある
クラスの影像エラー発生源、例えば信号取得中外部磁界
により安定磁界の影響する信号取得中の対象の剛体運動
と対象の変形及び計器ドリフト現象により生じる磁気共
鳴影像中の影像エラーを減少させる方法を提供すること
である。
【０００５】
【課題を解決するための手段】これを達成する為、本発
明による磁気共鳴影像方法は、磁気共鳴信号は少なくと
も第１及び第２のパルス及び勾配シーケンスにより得ら
れ、その信号サンプルは少なくとも第１及び第２のデー
タセットに蓄積され、各データセットのデータラインは
周波数領域で互いに交差し、交差データラインのデータ
の組合せから影像エラー源用推定が決定され、その後、
決められた推定に基づいて少なくとも１つのデータセッ
トのデータラインのデータが補正され、その後磁気共鳴
影像が積分変換により補正されたデータセットから決め
られるか、又はフーリエ変換を用いて、第３のデータセ
ットがデータセットの１つから決められ、第４のデータ
が推定から決められ、その後磁気共鳴影像が第３のデー
タセットの第４のデータセットとのコンボリューション
により決められることを特徴とする。
【０００６】本発明は、異なる影像エラー発生源が周波
数又はフーリエ領域でのデータラインのデータの取得中
少なくとも略一定であり、又は影像エラー発生源のモデ
ルが形成されうるという事実を用いる。実際、データラ
インの信号サンプリングは略１５ｍｓ又はそれ以下を要
する。これは、影像エラー発生源の結果が１又は少しの
スカラーパラメータにより表わされうることを意味す
る。本発明による方法は、これらのパラメータの推定及
びそれに基づいた磁気共鳴影像での影像エラーの補正に
基づいている。第１及び第２のデータセットの交差ライ
ンは影像補正がなされるのに基いたデータラインエラー
に関して情報を提供する。補正は望ましくは積分変換に
先立って実行されるが、変換がフーリエ変換のような直
線的変換である場合、適用可能であるコンボリューショ
ン定理に基いた変換の後に実行されうる。後者の場合に
おいて、第１の２つの多次元フーリエ変換は、非補正デ
ータ及びエラー発生源の推定で実行され、第３及び第４
のデータセットはコンボリューションによりそれから補
正された影像を得るよう形成される。例えば呼吸運動用
センサ又は外部磁界のような他の影像エラー発生源の測
定用センサのような生理系センサの使用を必要としない
で種々の影像エラー発生源による影像エラーの補正を可
能にする方法が提供される。
【０００７】日本国特開平２−４５０３５号は、ＭＲＩ
装置外部の磁界の変動を測定し、ＭＲＩ装置の別なコイ
ルによりそれに基いた適切な補償磁界を発生するために
磁界を測定するセンサを用いる方法を示していることが
注目される。従って、この方法は、標準ＭＲＩ装置に比
べて追加ハードウェアを必要とする。
【０００８】ニューヨーク、１９９０年８月１８−２４
日の医学での磁気共鳴の学会、アブストラクトの本５６
１頁の論文「変換運動アーティファクト補正のアルゴリ
ズム」に対象の剛体運動の場合の影像を補正する方法が
説明され、その方法は単一影像に印加された後処理に基
いていることが注目される。対象の運動がデータライン
のデータの取得中充分遅いとすると、対象シフトにより
生じる位相エラーは、対象が有限の次元を有する前の仮
定下で信号の大きさから推定される。
【０００９】本発明による方法の実施例は、第１のパル
ス及び勾配シーケンスは、励起パルスと勾配方向が一致
しないパルス符号化勾配及び読出し勾配とからなり、第
２のパルス及び勾配シーケンスは少なくとも第１のパル
ス及び勾配シーケンスからずれる励起パルスと、位相符
号化勾配及び読出し勾配とからなり、位相符号化勾配及
び読出し勾配が変換されたことを特徴とする剛体及び非
剛体運動で生じる磁気共鳴影像中の影像エラーを減少す
るのに特に適している。相互に交差するデータラインの
剛体は周波数領域、いわゆるＫ平面で得られ、複素（振
幅、位相）磁気共鳴信号の信号サンプルはマップされ
る。データセットの１つに蓄積されたデータに印加され
る多次元積分変換、例えばフーリエ変換は、例えばグレ
ースケール影像として表示スクリーン上に表示されうる
対象の磁気共鳴分布となる。本発明により、データライ
ンは、第１のデータセットのデータラインｊのサンプル
ｉが第２のデータセットのデータラインｉのサンプルｊ
に対応するよう位相符号化勾配の選択で符号化される。
影像エラーを含む、未知数を含む式の結果の数から、推
定に影像エラーに対して決められ、データセットの１つ
のデータは該推定により補正される。
【００１０】剛体運動で生じる影像エラーを減少させる
のに特に適している本発明による方法の更なる実施例
は、データラインの交差の点のデータの位相差がデータ
の組み合せとして用いられ、データラインと関連する対
象シフトが影像エラー発生源の推定として該組み合せか
ら決められることを特徴とする。この点で、対象が１つ
のデータラインのデータの取得中動かないとする。運動
なしに、Ｋ空間の交差でのデータは等しいが、対応する
格子点のデータがもはや等しくないように運動はデータ
ラインのシフトを生じる。推定を用いると、データは、
データセットの対応するデータが等しくなるように補正
される。補正されたデータセットの１つから又はその組
み合せからの影像の再構成が実質的により少ない影像エ
ラーを含む影像を発生する。
【００１１】非剛体運動、例えば呼吸による運動で生じ
た影像エラーを減少させるのに特に適している本発明に
よる方法の更なる実施例は、データラインの交差の点の
データの位相差がデータの組み合せとして用いられ、対
象のストレッチ係数が影像エラー発生源の推定とに該組
み合せから得られることを特徴とする。この実施例にお
いて、運動は主に腹部−臀部であり、運動は対象の直線
的伸張又は圧縮として作れる。即ち、運動の影響下で組
織が局部的（例えばｙ方向の）距離ｃ，ｙに亘って変位
するとする。ここでｃはデータライン毎に異なる伸張係
数である。補正は又２つの「互いに垂直な」データセッ
トを基にして再び実行される。
【００１２】該実施例において、水又は脂肪の予飽和が
公知の方法で実現されうる。その場合に、更に少ない影
像エラーを含む影像が得られる。
【００１３】安定磁界の磁界不均一性をより少なく受け
やすい本発明による方法の実施例は、第１のパルス及び
勾配シーケンスは、励起パルスと、勾配方向が一致しな
いパルス符号化勾配及び読出し勾配とからなり、第２の
パルス及び勾配シーケンスは少なくとも第１のパルス及
び勾配シーケンスからずれる励起パルスと、位相符号化
勾配及び読出し勾配とからなり、第２のパルス及び勾配
シーケンスは同時に読出し勾配を印加される別な位相符
号化勾配からなることを特徴とする。水又は脂肪のよう
な異なる成分又は他の成分を含む対象の安定磁界の磁界
不均一及び化学シフトが第１及び第２のデータセットの
共鳴信号の共鳴の共鳴に差を生じるなら、差は対応する
サンプルのサンプリング時間で生じる。これらの時間差
は、データセットの全データラインに対する全対応サン
プルに対して励起とサンプリング間の時間経過を保つこ
とにより除去される。
【００１４】該時間差を推定する本発明による方法の他
の実施例は、第１のパルス及び勾配シーケンスは、励起
パルスと、勾配方向が一致しないパルス符号化勾配及び
読出し勾配とからなり、第２のパルス及び勾配シーケン
スは少なくとも第１のパルス及び勾配シーケンスからず
れる励起パルスと、位相符号化勾配及び読出し勾配とか
らなり、別な勾配は第１のパルス及び勾配シーケンスの
位相符号化パルスに重畳され、また読出し勾配中に印加
され、第２のパルス及び勾配シーケンスの位相符号化勾
配に読出し勾配中に印加される逆の別勾配が重畳される
ことを特徴とする。
【００１５】推定は該実施例で繰返して決定されうる。
【００１６】外部磁界により影像エラーを減少するのに
特に適している本発明による方法の実施例は、第１のパ
ルス及び勾配シーケンスは、励起パルスと、勾配方向が
一致しないパルス符号化勾配及び読出し勾配とからな
り、第２のパルス及び勾配シーケンスは又少なくとも第
１のパルス及び勾配シーケンスからずれる励起パルス
と、位相符号化勾配及び読出し勾配とからなり、読出し
勾配と一致する更なる勾配は、読出し勾配の勾配からず
れた勾配方向を印加され、３次元影像の場合に、更なる
位相符号化勾配が第１及び第２のパルス及び勾配シーケ
ンスに印加され、第３のパルス及び勾配シーケンスは推
定及び補正用の第３のデータセットを得るため印加さ
れ、励起パルスと磁気共鳴信号から読出す間アクティブ
である３つの勾配とからなり、その勾配方向は一致して
いないことを特徴とする。安定磁界が１つのデータライ
ンのデータの取得中ほどよく一定であるとすると、影像
エラーはデータラインに亘って直線的に増える位相エラ
ーとして作られうる。磁界妨害が影像に比較的に強い影
響を有するいわゆる勾配エコー法での如く、データライ
ンの取得時間がエコー時間に比べて短い場合、影像エラ
ーはデータラインに亘って一定の位相シフトとして更に
作られうる。データセットの影像エラー減少は、２次元
の場合１つ又は数データラインからなるデータセットを
基にして実行されうる、３次元の場合、補正さるべきデ
ータセットより実質的に少ないデータラインからなるデ
ータセットを基にして実行されうる。
【００１７】外部磁界による影像エラーを減少させるの
に適し、該位相推定を用いる本発明による方法の実施例
は、非常に単純である。２次元の場合において、補正用
に用いるべきデータセットは単に１つのデータラインか
らなる。雑音効果に関して改善を行なう為に、補正デー
タセットは数データラインからなり、位相エラーを推定
するように平均値がとられうる。
【００１８】全実施例において、推定は影像品質を改善
するために雑音に関して重み付けされうる。
【００１９】
【実施例】本発明を以下図面を参照して詳細に説明す
る。
【００２０】図１は、ＲＦ電磁パルスを送信／受信コイ
ル４を介して対象５に送信する送信手段２と、均一安定
磁界に位置する対象５にＲＦ電磁パルスで発生される各
磁気共鳴信号を受信する受信手段３とからなる本発明に
よる磁気共鳴影像装置１を概略的に示す。送信／受信コ
イルは単一コイルであるが、代りに別体の送信コイル及
び別体の受信コイルを設けてもよい。装置１は安定磁界
を発生する手段６からなる。手段６は磁石コイル７から
なり、抵抗磁石又は超伝導磁石の場合には、直流電源８
からなる。永久磁石の場合には、直流電源８はない。装
置１の動作中対象は磁石コイル７内に配置され、僅か過
剰な該スピン（磁気モーメントを有する核）は平衡の状
態で安定磁界の方向に向けられる。巨視的な見地から、
これは平衡磁化である磁化Ｍ0 として考えられる。装置
１は又送信手段２と受信手段３に結合される処理手段９
と、処理手段９及び送信手段に結合される処理コンピュ
ータ１０と、共鳴信号から決定され、受信手段３で受信
され、復調され、その信号サンプリング（共鳴信号の検
出）の後核磁化分布を表示する表示手段１１とからな
る。特に、送信手段２は、搬送波信号を発生するＲＦ発
振器１３と、搬送波信号の振幅及び／又は位相及び周波
数変調用変調器１４と、電力増幅器１５と、送信／受信
コイル４に結合される方向性結合器１６とからなる。Ｒ
Ｆ発振器１３は処理手段９に結合され、変調器１４は処
理コンピュータ１０に結合される。励起パルスが、プロ
グラム手段１２の制御下で送信手段２を介して、例えば
陽子のいわゆるラーモア周波数の近辺の周波数内容で、
対象５に印加される時、磁気共鳴信号が発生し、それか
ら陽子核スピン分布は例えばフーリエ変換によりプログ
ラムされた手段１２で決定されうる。共鳴信号を受信す
る受信手段３は方向性結合器１６と、受信及び復調ユニ
ット１７とからなる。ユニット１７は、出力信号が第１
及び第２のアナログディジタル変換器１８，１９により
夫々サンプルされる二重位相感知検出器でよい。Ａ／Ｄ
変換器１８及び１９は処理手段９に結合される。或い
は、送信手段２及び受信手段３はいわゆる位相コーヒー
レンスディジタル送信器／受信器により形成されうる。
装置１は又安定均一磁界に重畳される勾配磁界を発生す
る手段２０からなる。手段２０は勾配磁界Ｇｘ，Ｇｙ及
びＧｚを各々発生する。勾配磁石コイル２１，２２及び
２３と、処理コンピュータ１０で活動され：別々に活動
可能である勾配磁石コイル２１，２２及び２３に通電す
るよう働く電源２４とからなる。図示の実施例では、勾
配磁石コイルの空間の装置は、勾配磁界の磁界方向が均
一安定磁界の方向に一致し、勾配方向が図１の３つの相
互に垂直な軸Ｘ，Ｙ及びＺにより示される如く、互いに
垂直に延在するようにされている。フーリエ変換により
再構成されうる核ヒスピン分布からの磁気共鳴信号はい
わゆるパルス及び勾配シーケンスにより得られる。
【００２１】図２の（Ａ）は本発明による方法の実施例
のパルス及び勾配シーケンスを示す。この方法によれ
ば、第１のいわゆるスピンエコーパルス及び勾配シーケ
ンスｐｓ１は時間ｔの関数として対象５に印加される。
パルス及び勾配シーケンスｐｓ１は、対象５の核スピン
を励起するため時点ｔ＝ｔ１で送信手段２により送信／
受信コイル４を介して対象５に印加される励起パルスｅ
ｐ１よりなり、また、ディフェーズされた核磁化を再集
束するため時点ｔ＝ｔ２で印加される再集束パルスｒｐ
１からなる。時点ｔ＝ｔ３で、受信手段３は、対象５か
ら生じ、送信／受信コイル４で検出される磁気共鳴信号
を受信する。ＲＦ電磁励起パルスの適用中勾配磁界Ｇｚ
は対象のスライスの選択に対しアクティブであり、即ち
磁気共鳴信号は対象５のスライスから選択的に生じる。
更に、パルスｅｐ１とｒｐ１との間に、核スピンの位相
符号化する位相符号勾配Ｇｙが印加され、核スピンの周
波数符号化するための磁気共鳴信号の読出し中読出勾配
Ｇｘが印加される。勾配Ｇｘ，Ｇｙ及びＧｚは勾配磁石
コイル２１，２２及び２３を介して安定磁界に重畳され
る。シーケンスｐｓ１は、位相符号化勾配Ｇｙの異なる
値に対して多数回例えば２５６回繰り返される。Ｇｙの
各値に対して、磁気共鳴信号が受信及び復調ユニット１
７で受信され、復調され、サンプルの形式で、第１のデ
ータセットの蓄積の為Ａ／Ｄ変換器１８及び１９を介し
て処理手段９に印加される。各Ｇｙに対して、Ｋ空間の
データラインが満たされる。
【００２２】プログラムされた手段１２を用いて、核磁
化分布はフーリエ変換によりデータセットから再構成さ
れ、分布は表示手段１１により磁気共鳴影像の形式で表
示されうる。本発明により、時点ｔ＝ｔ１’での励起パ
ルスｅｐ２と、時点ｔ＝ｔ２’での再集束パルスｒｐ２
と、第１のシーケンスｐｓ１に関して交換された選択勾
配Ｇｚ及び更なる勾配Ｇｘ及びＧｙとからなる第２のパ
ルス及び勾配シーケンスｐｓ２が印加される。これは、
第１のシーケンスにより得られたデータラインに垂直に
延在するＫ空間でのデータラインとなる。第２のシーケ
ンスにより得られた磁気共鳴信号のサンプルは第２のデ
ータセットに蓄積される。
【００２３】図２の（Ｂ）は図２の（Ａ）に示すシーケ
ンスに関連したＫ空間での交差ラインを示す。２次元影
像に関する本実施例では、核スピン密度関数Ｐ0 （ｘ，
ｙ）のフーリエ変換ＦはＫ空間の軸Ｋｘ及びＫｙに沿っ
てプロットされる。本発明はいわゆるＫ軌道記述の言葉
で下記に説明される。Ｋ軌道記述の本質は、時間マップ
での勾配が磁気共鳴信号Ｓ（ｔ）で直接に対象の空間・
周波数領域内容を変えることにある。
【００２４】
【数１】
【００２５】Ｋ（ｔ）は勾配の整数であり、又空間・周
波数座標のベクトル、即ち
【００２６】
【数２】
【００２７】であり、ここで、Ｆは対象の磁気共鳴分布
であるｆ（ｘ）のフーリエ変換である。更に、上式
【００２８】
【数３】
【００２９】において、ｘはベクトルの形式での空間座
標の２次元又は２次元セットであり、γはいわゆる磁気
回転比であり、Ｔ2 は横緩和時間であり、Ｂは磁界の強
度である。Ｓ（ｔ）は対象に対応する空間周波数分布Ｆ
〔Ｋ（ｔ）〕の重み付け（及び雑音による妨害）観測と
して考えられるべきものである。磁気共鳴影像におい
て、核スピンの符号化処理は時間的に変わる勾配で決定
された一組の軌道に沿ったサンプリングからなり、復号
化処理は元の空間分布の推定である離散的影像にＫ空間
の該サンプルを変換することからなる。パルス及び勾配
シーケンスがいわゆるフーリエツォイクマトグラフィ
（ｚｅｕｇｍａｔｏｇｒａｐｈｙ）方法と関連する場
合、一連の異なる勾配は、勾配Ｇの瞬時値で、式ｄＫ／
ｄｔ＝Ｇ（一定測定勾配に対して、Ｋｘ＝γ・Ｇｘ・
ｔ）によって決められる速度でＫｘＫｙ平面の「偏り」
として考えられ、Ｓ（ｔ）の瞬時値はＫｘＫｙ平面に到
る点でＦ値を発生する。Ｋ空間のサンプルのフーリエ変
換は対象の核スピン磁化分布を発生する。Ｋ空間の直交
格子でプロットされるサンプルは離散的フーリエ領域影
像のデータ値として考えられるべきであって、離散的空
間領域影像はフーリエ領域影像を２Ｄ−ＦＴ又は３Ｄ−
ＦＴにすることによって得られうる。サンプルが（いわ
ゆるエコー平面パルス又は勾配シーケンスに対する場合
である如く）Ｋ空間の格子点と一致しない場合、ある種
の補間方法は、Ｆ（Ｋ）が離散的ＦＴの実行前格子点に
対して推定されることを確実にするよう用いられるべき
である。
【００３０】図６は、図２の（Ａ）に示す如き第１実施
例に対する式を示し、ここでは、記載シーケンスｐｓ１
及びｐｓ２は図２の（Ｂ）に示す如くＫ空間にアップさ
れた第１及び第２のデータセットを作った。この実施例
は剛体運動で生じた影像エラーを減少させるのに適して
いる。かかる剛体運動は、例えば対象が更に押さえつけ
手段を用いるにも拘らず動きを行なう人間対象の膝の影
像の場合、生きている対象の手足のＭＲ影像に生じう
る。そのような場合には、各データラインはその元の位
置の対象から生じるデータを含まない。データラインは
その元の位置に関して（δｘ，δｙ）画素によりシフト
された対象からのデータも一般的に含む。シフトは、勿
論、従来公知でなく、データライン毎に変化しうる。図
示の実施例では、水成分又は脂肪成分のいずれかが公知
の方法で予め飽和され、安定磁界の異質性が無視出来る
ほど小さいとする。対象の位置は、第１のデータセット
ではδｘ及びδｙがＫｘに依存せず、専らＫｙに依存す
るように、１つのデータラインのサンプリング中固定さ
れると考えられるとした。プログラムされた手段がわず
かにより複雑になるとしても、これらの仮定が有効でな
い場合、それらは原証的に上記の方法に対応する。同様
に、エラーがＫｘだけを依存することが第２のデータセ
ットに対して成り立つ。各シフトは、第１のデータセッ
トＡ用
【００３１】
【数４】
【００３２】として示される。本発明による方法の普遍
性を損なうことなく、たとえば、パルス及び勾配シーケ
ンスｐｓ１の位相符号化ステップｎa の数及びデータラ
インの数がパルス及び勾配シーケンスｐｓ２での位相符
号化ステップｎ6 の数に等しく、サンプル又は周波数符
号化ステップｎS の数がデータセット当たりの位相符号
化ステップの数に等しいとする。更に、シーケンスｐｓ
１及びｐｓ２での位相符号化ステップは、第１のデータ
セットＡのデータラインｊのサンプルｉがＫ空間でのデ
ータラインの交差の点での対応データを得るためにデー
タセットＢのデータラインｉのサンプルｊに対応するよ
うに選択されるべきである。本実施例では、
【００３３】
【数５】
【００３４】を有する式（ｎa ＋ｎb ）・ｎS が未知で
ある場合がある。（即ち、
【００３５】
【数６】
【００３６】が実際の影像データ及びデータセットのデ
ータラインの（ｎa＋ｎb ）の未知のエラーで未知であ
る）本発明による方法において、先ずデータラインのエ
ラーが決められ、それを基にしてデータラインが補正さ
れる。原則的に二倍の全取得時間が用いられるが、本発
明により得られた影像が２つの等しいデータセットの単
純平均の場合より更によいことが注目される。対象５が
完全に停止していた場合、第１及び第２のデータセット
Ａ及びＢの対応データは「理想的」データＳ0 （Ｋｘ，
Ｋｙ）と等しくなる。しかし、剛体運動の場合には、各
実際のデータＳA （Ｋｘ，Ｋｙ）及びＳB （Ｋｘ，Ｋ
ｙ）は図６の式〔１〕及び〔２〕により与えられる。本
発明により、実際のデータＳA 及びＳB は、それらが互
いに及びＳ0 に略等しくなるようにデータラインシフト
の推定に基づいて補正される。先ず、複素データは、＊
が複素共役を示す式〔３〕で表わされる如くＫ空間の格
子点で互いに分割される。１つの画素の大きさのオーダ
での小さいシフトに対し、そこでのｅべき数は、近似の
虚数部をとることで式〔４〕が得られるようにそのべき
級数展開の第１の２つの項で近似される。それを基にし
て、量Ｈｘ（Ｋｘ，Ｋｙ）及びＨｙ（Ｋｘ，Ｋｙ）が各
式〔５〕及び〔６〕で定義され、その量は式〔４〕から
得られる如く決め、ｘ方向およびｙ方向でのデータセッ
トＡ及びＢのデータラインと関連した対象５の差シフト
は各差シフトの１つが毎回０となる。式〔５〕及び
〔６〕において、文字εはこの方法で０による分割を防
ぐために小さい数、例えばε＝０．５、を示す。これら
の所定の量及び式〔４〕を基にして、式〔７〕から〔１
０〕に示す如く、データセットＡ及びＢのデータライン
と関連した対象５の各シフト８に対する所定の推定があ
る。ここで〜は推定又は近似を示す。定義Ｈｘ及び式
〔４〕による近似での式〔１１〕及び〔１２〕を介した
【００３７】
【数７】
【００３８】に対する近似として、式〔１３〕による表
現が用いられ、式〔１１〕に続く式〔１２〕は、
【００３９】
【数８】
【００４０】とし、Ｋｘは値の対称的範囲を有し、一方
式〔１２〕に続く式〔１３〕は式〔１２〕の第２の項が
絶対値で他の項に対して大きいとする。第１の近似にお
いて、式〔１２〕の第１及び第３の項がインコーヒーレ
ントランダム変数を含み、第２項に対して
【００４１】
【数９】
【００４２】からｎS として関連するので、後者のアプ
ローチは正当化される。式〔１２〕の最終項は、Ｋｘ用
値の対称範囲の仮定によって零である。他のシフトδに
対する推定が同様な方法で決定される時、データ値ＳA
及びＳB は式〔１４〕及び〔１５〕により補正されう
る。式〔１４〕及び〔１５〕による近似が更にＳ0 から
あまりにも離間している場合、近似処理が繰返されう
る。多数の繰返し、典型的には１から４回の繰返しの
後、磁気共鳴影像は補正データ
【００４３】
【数１０】
【００４４】から再構成され、ここでｎは繰返の数を示
す。上記の方法は雑音によりもたらされた不確かさを考
慮するよう適合されうる。式〔１６〕から明らかな如
く、推定は次に重み付け係数Ｗ（Ｋｘ，Ｋｙ）により重
み付けされる。式〔１７〕は、雑音による推定のエラー
が最小化されるので、望まれるべき重み付け係数を示
す。式〔１７〕において、Ｎ2 はデータ中の雑音の分散
である。
【００４５】図７は、上記シーケンスｐｓ１及びｐｓ２
が、図２の（Ｂ）に示す如くＫ空間にマップされた第１
及び第２のデータセットを得るのに用いられる図２の
（Ａ）の第２の実施例に対する式を示す。この実施例で
は、非剛体運動、例えば生きている対象５の呼吸運動で
生じる影像エラーを減少するのに適している。例えば、
横影像は生きている対象５の胸部領域からなりうる。胸
部領域において、運動は主に腹部−臀部であり、運動は
直線的伸張又は圧縮として作られる、即ち、組織が局部
的にｙ方向の距離ｃ，ｙに亘って変位され、ここでｃは
データライン毎に異なるストレッチ係数であるとする。
図６による第１の実施例における様なデータセットの大
きさに関して仮定すると、以下のことがより成り立つ。
図７の式〔１８〕において、「理想的」信号Ｓ0 （Ｋ
ｘ，Ｋｙ）が与えられ、式〔１９〕において、対象の伸
張により歪まされた信号即ち第１のデータセットＡの歪
みデータが与えられる。小さいゆがみ対して，、式〔１
９〕の第１のｅべき数はそのべき級数展開の最初の２つ
の項で近似され、従って式〔２０〕を得る。コンボリュ
ーション演算子
【００４６】
【数１１】
【００４７】の項で書かれた式〔２０〕は、Ｇ（Ｋｙ）
が関数２πｙ／ｎSのフーリエ変換である式〔２１〕を
作り、該フーリエ変換は式〔２２〕で示される。同様
に、式〔２３〕は第２のデータセットＢの歪みデータの
近似を生じる。第２の実施例において、歪んだＳA 及び
ＳB は、それらが出来るだけ「理想的」信号Ｓ0 に近似
するよう補正されるべきである。第１の実施例でのよう
に、量Ｈ（Ｋｘ，Ｋｙ）は式〔２１〕及び〔２３〕に基
づいた式〔２４〕で定義され、ここでＳ0 （Ｋｘ，Ｋ
ｙ）は、第１の例では（ＳA ＋ＳB ）／２で近似され
る。Ｈ（Ｋｘ，Ｋｙ）は伸張係数ＣA 及びＣB で表わさ
れる式〔２５〕に示す形式に略等しいことが単に示され
うる。伸張が、呼吸による運動で仮定されうる如く、デ
ータライン間で多少ランダムであると仮定されうる場
合、伸張係数は式〔２６〕及び〔２７〕に示す如く推定
され、データは分かった推定により補正されうる。雑音
及び他の不確かさを考慮する時、本発明による方法は、
Ｗ（Ｋｘ，Ｋｙ）が第１の例での如く重み付け係数であ
る式〔２８〕から〔３４〕に示す如くわずかに適合され
うる。コンボリューションの期間を制限するために、Ｇ
（Ｋｙ）の最初の４から６の最も大きい素子だけを用い
ることが許される。方法がわずかだけより複雑になるに
しても、エラーが運動のより複雑なパターン、例えば対
象の患者らの頭部の横影像の頭部の回転によるが、呼吸
により生じる、腹のサジタル影像の長手方向（尾−頭方
向）の非直線的伸張により生じる場合、わずか影像エラ
ーを含む影像はそれにも拘らず得られる。
【００４８】図３の（Ａ）は、励起パルスｅｐ３と、再
集束パルスｒｐ３と、スライス選択勾配Ｇｚと、位相符
号化勾配Ｇｙと、測定勾配Ｇｘとからなる第１のパルス
及び勾配シーケンスｐｓ３からなる本発明による方法の
パルス及び勾配シーケンスの更なる例を示す。パルスｅ
ｐ３及びｒｐ３は時点ｔ１及びｔ２で夫々印加され、磁
気共鳴信号は略時点ｔ３で検出される。シーケンスｐｓ
３は位相符号化勾配の多数のステップ、例えば２５６ス
テップに対し繰り返される。更に、時点ｔ１’での励起
パルスｅｐ４と時点ｔ２’での再集束パルスｒｐ４とか
らなるパルス及び勾配シーケンスｐｓ４が示される。シ
ーケンスｐｓ４は、位相符号化ステップの数の半分がＫ
空間の一方の側に位置し、固定、追加勾配Ｇｙが各シー
ケンスｐｓ４の磁気共鳴信号から読出す間に印加される
シーケンスｐｓ３からずれる。図３の（Ｂ）に示す如
く、走行のこの方法は、斜データラインがシーケンスｐ
ｓ４と関連する間、水平データラインがシーケンスｐｓ
３と関連するＫ空間での交差データラインのパターンに
なる。
【００４９】図３の（Ｃ）は本発明による方法の他の例
と関連するＫ空間での交差データラインを示す。関連し
たデータセットを有するデータラインのこのパターン
は、通例のシーケンスのように周波数符号化勾配を含む
が、Ｇｙと追加勾配の和と通例スピンエコー又はスピン
ラップシーケンスのような通例シーケンスでの如き波形
を有するＧｙ及び追加勾配間の差とで各形成される位相
符号化勾配が印加される２つのパルス及び勾配シーケン
ス（図示せず）によって得られる。追加勾配は又、読出
勾配中正及び負の感覚で印加される。従って、データラ
インを相互に、垂直に交差するパターンが、データライ
ンがＫ空間の軸に関して＋４５°及び−４５°の各角度
で被うＫ空間で得られる。図３の（Ａ）から図３の
（Ｃ）に示すシーケンスは、安定磁界の異質又は化学シ
フトによる問題が起こる図２に示す如き対応サンプル間
の時間差による問題を起こさない。図３の（Ａ）に示す
シーケンスにおいて、サンプルｉの励起とサンプリング
の間の時間経過は両データセットの全データラインに対
して一定に保たれる。推定及び補正式は類似した方法で
決定されうる。
【００５０】図４の（Ａ）は例えばＭＲＩ装置近くの動
くエレベータ、通る電車等からの外部磁界による影像エ
ラーを減少させるに適した本発明によるパルス及び勾配
シーケンスの例を示す。この図は、励起パスｅｐ５と公
知のスピンワープ又は磁界エコーシーケンスでのような
勾配とからなり、位相符号化勾配の代わりに、シーケン
スｐｓ５でのように、勾配Ｇｘと等しく、それと同時に
印加される単一勾配が印加されるが、その勾配方向が異
なるシーケンスｅｐ５からずれた磁界エコーシーケンス
ｐｓ６からなるパルス及び勾配シーケンスｐｓ５を示
す。これは図４の（Ｂ）に示す如くＫ空間の交差ライン
のパターンになる。補正さるべき第１のデータセットＡ
は図４の（Ｂ）の水平データラインのデータを含み、補
正用に使われる第２のデータセットＢは図４の（Ｂ）の
斜データラインのデータを含む。磁界妨害は１つのデー
タライン用取得時間に関して遅いが、完全磁気共鳴影像
用データを得るための全体取得時間に関して速いとす
る。図示の例は、磁界変化で生じた効果が最も大きい磁
界エコーシーケンスに使われるが、なかんずくスピンエ
コーシーケンス（９０°−１８０°エコー）のような他
のシーケンスに対しても用いられうる。磁界エコーシー
ケンスにおいて、安定磁界の磁界妨害はデータラインの
サンプルされたデータの周波数オフセットを生じる。磁
界がデータラインの取得中一定であるとするので、エラ
ーはデータラインに亘って位相エラーを直線的に増すと
してモデル化されうる。取得時間が励起及び取得間の時
間経過であるエコー時間と比較して短い場合、エラーは
データラインに亘って一定の位相シフトで近似されう
る。本実施例において、データセットは夫々がＮS サン
プルを有するＮA データラインを含み、但しＮA＝
ＮS 、データセットＢはＮS サンプルを有する１つのデ
ータラインを含む。セットＢは、「理想的」場合におい
て、ＳB （ＫB ）がＳA （ＫB ，ＫB ）と同じ信号を含
むようにセットアップされる。前記の理由によって、励
起パルスｅｐ５とＳA （Ｋｘ，Ｋｙ）のサンプリング間
の時間経過が励起パルスｅｐ６とＳB （Ｋｘ）のサンプ
リング間の時間経過に等しいことが確実とされる。
【００５１】図４の（Ｃ）は信号Ｓ0 （Ｋｘ，Ｋｙ）が
理想的に生じる仮定に基づいて図４の（Ａ）の実施例に
対する式を示す。式〔３５〕及び〔３６〕は該位相エラ
ーにより検出される実際のデータを示す。ここで、Ｎｘ
はサンプルの数であり、ｔac q はデータラインの取得時
間である図示のアプローチはｔacq の小さい値に対して
成り立つ。式〔３７〕から分かる如く、φB は、中央サ
ンプリングが現実でなければならないという知識を基に
補正され、その後、全体データラインＢは式〔３８〕に
応じて補正される。次に、データセットＡの全データラ
インの位相エラーは式〔３９〕により補正され、該位相
エラー推定はデータセットＡの位相補正に対して用いら
れる。磁気共鳴影像は補正データセットＡから再構成さ
れる。本発明による方法は更に改良されうる。例えば、
推定補正係数は雑音の効果を減少させる為に重み付けさ
れうる。データセットＢがいくつかのデータラインから
なる場合、雑音の効果は位相エラーの推定を平均化する
ことで更に減少されうる。データセットＢがいくつかの
データラインからなる場合にＫ空間のデータライン間に
ある距離が存在する場合、この事実は直線的に増加する
位相エラー成分を推定するのに用いられうる。物理的不
完全により、データセットＢのデータラインは、Ｓ
B （Ｋｘ）が正確にＳA （Ｋｘ，Ｋｙ）に対応せず、む
しろＳA （Ｋｘ，Ｋｙ＋Δ）に対応するようにＫ空間で
わずかにシフトされた。一般的に、Δは１より少ない非
整数である。その場合において、隣る値ＳB （Ｋｘ）と
ＳB （Ｋｘ−１）からの値ＳB ’（Ｋｘ−Δ）及びＳA
（Ｋｘ，Ｋｙ）とＳA （Ｋｘ−１，Ｋｙ）からの値
ＳA ’（Ｋｘ−Δ、Ｋｙ）を補間することが必要であ
る。次に位相エラーφA （Ｋｙ）はＳB ’（Ｋｘ−Δ）
とＳA ’（Ｋｘ−Δ，Ｋｙ）間の位相差として計算され
る。
【００５２】図５は、外部磁界による影像エラーが減少
し、図４の（Ａ）に対応するパルス及び勾配が対応する
方法で表わされる３次元影像用パルス及び勾配シーケン
スを示す。図４の（Ａ）と比較して、励起パルスｅｐ７
を有する更なるパルス及び勾配シーケンスｐｓ７が設け
られる。シーケンスｐｓ５はシーケンスｐｓ６のような
更なる位相符号化勾配Ｇｚからなる。シーケンスｐｓ７
は、読出し勾配Ｇｘと、更にこれと一致した勾配Ｇｙ及
びＧｚとからなる。シーケンスｐｓ７の各勾配の勾配方
向が異なる。シーケンスｐｓ５，ｐｓ６及びｓｐ７はＫ
空間のデータラインを交差する３つのデータセットＡ，
Ｂ及びＣを発生する。データセットＡは、補正の後そこ
から磁気共鳴影像を再構成するのに用いられ、「単一
斜」データセット及び「二重斜」データセットＣは補正
目的の為に用いられる。データが立方であるとし、励起
とサンプリング間の時間経過が対応するサンプルに対し
同じであることが確実とされる。
【００５３】図９は図８に示す式と類似の図５の３次元
実施例に対する式を示す。式〔４１〕から〔４３〕は近
似データからなり、式〔４４〕から〔４９〕は順次にＳ
A 用推定ステップを提供する。
【図面の簡単な説明】
【図１】本発明による磁気共鳴影像装置の概略を示す図
である。
【図２】（Ａ）は本発明による方法の実施例のパルス及
び勾配シーケンスを示す図、（Ｂ）は（Ａ）に示すシー
ケンスと関連したＫ空間の交差データラインを示す図で
ある。
【図３】（Ａ）は本発明による方法のパルス及び勾配シ
ーケンスの更なる実施例を示し、（Ｂ）は図３の（Ａ）
に示すシーケンスに関連したＫ空間の交差データライン
を示し、（Ｃ）は本発明による方法の他の実施例に関連
したＫ空間の交差データラインを示す図である。
【図４】（Ａ）は外部磁界により影像エラーを減少させ
るのに適している本発明によるパルス及び勾配シーケン
スの実施例を示し、（Ｂ）は図４の（Ａ）に示すシーケ
ンスに関連したＫ空間の交差データラインを示す図であ
る。
【図５】外部磁界により影像エラーを減少させるための
３次元影像用パルス及び勾配シーケンスを示す図であ
る。
【図６】図２の（Ａ）に示す如く第１の実施例用の式を
示す図である。
【図７】図２の（Ａ）に示す如く第２の実施例用の式を
示す図である。
【図８】図４の（Ａ）による実施用の式を示す図であ
る。
【図９】図５による実施例用の式を示す図である。
【符号の説明】
１ 磁気共鳴影像装置 ２ 送信手段 ３ 受信手段 ４ 送信／受信コイル ５ 対象 ６ 磁界発生手段 ７ 磁石コイル ８，２４ 電力供給源 ９ 処理手段 １０ 処理コンピュータ １１ 表示手段 １２ プログラム手段 １３ ＲＦ発振器 １４ 変調器 １５ 電力増幅器 １６ 方向性結合器 １７ 復調ユニット １８，１９ アナログディジタル変換器 ２０ 磁界勾配発生手段 ２１，２２，２３ 勾配磁石コイル
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.5 識別記号 庁内整理番号 ＦＩ 技術表示箇所 9118−2Ｊ Ｇ０１Ｎ 24/08 Ｙ (72)発明者 ミハ フデラー オランダ国 5621 ビーエー アインドー フエン フルーネヴアウツウエツハ １番 地 (72)発明者 アントン フラン メルコツプ オランダ国 5621 ビーエー アインドー フエン フルーネヴアウツウエツハ １番 地 (72)発明者 ミカエル ヨゼフ デユイブステイン オランダ国 5621 ビーエー アインドー フエン フルーネヴアウツウエツハ １番 地 (72)発明者 ゲリツト ヘンドリツク ヴアン イピー レン オランダ国 5621 ビーエー アインドー フエン フルーネヴアウツウエツハ １番 地

Claims (14)

    【特許請求の範囲】
  1. 【請求項１】 安定磁界に配置された対象からの磁気共
    鳴信号から形成された磁気共鳴影像中の影像エラーを減
    少させる磁気共鳴影像方法であって、磁気共鳴信号は少
    なくとも第１及び第２のパルス及び勾配シーケンスによ
    り得られ、その信号サンプルは少なくとも第１及び第２
    のデータセットに蓄積され、各データセットのデータラ
    インは周波数領域で互いに交差し、交差データラインの
    データの組合せから影像エラー源用推定が決定され、そ
    の後、決められた推定に基いて少なくとも１つのデータ
    セットのデータラインのデータが補正され、その後磁気
    共鳴影像が積分変換により補正されたデータセットから
    決められるか、又はフーリエ変換を用いて、第３のデー
    タセットがデータセットの１つから決められ、第４のデ
    ータが推定から決められ、その後磁気共鳴影像が第３の
    データセットの第４のデータセットとのコンボリューシ
    ョンにより決められることを特徴とする磁気共鳴影像方
    法。
  2. 【請求項２】 第１のパルス及び勾配シーケンスは、励
    起パルスと、勾配方向が一致しないパルス符号化勾配及
    び読出し勾配とからなり、第２のパルス及び勾配シーケ
    ンスは少なくとも第１のパルス及び勾配シーケンスから
    ずれる励起パルスと、位相符号化勾配及び読出し勾配と
    からなり、位相符号化勾配及び読出し勾配が交換された
    ことを特徴とする請求項１記載の磁気共鳴影像方法。
  3. 【請求項３】 データラインの交差の点のデータの位相
    差がデータの組み合せとして用いられ、データラインに
    関連した対象シフトが影像エラー源の推定として該組み
    合せから決められることを特徴とする請求項２記載の磁
    気共鳴影像方法。
  4. 【請求項４】 データラインの交差の点のデータの位相
    差がデータの組み合せとして用いられ、対象のストレッ
    チ係数が影像エラー源の推定として該組み合せから生じ
    ることを特徴とする請求項２記載の磁気共鳴影像方法。
  5. 【請求項５】 対象の水又は脂肪成分は、磁気共鳴信号
    の取得以前に予め飽和されることを特徴とする請求項
    ２，３又は４記載の磁気共鳴影像方法。
  6. 【請求項６】 第１のパルス及び勾配シーケンスは、励
    起パルスと、勾配方向が一致しないパルス符号化勾配及
    び読出し勾配とからなり、第２のパルス及び勾配シーケ
    ンスは少なくとも第１のパルス及び勾配シーケンスから
    ずれる励起パルスと、位相符号化勾配及び読出し勾配と
    からなり、第２のパルス及び勾配シーケンスは同時に読
    出し勾配を印加される別な位相符号化勾配からなること
    を特徴とする請求項１記載の磁気共鳴影像方法。
  7. 【請求項７】 第１のパルス及び勾配シーケンスは、励
    起パルスと、勾配方向が一致しないパルス符号化勾配及
    び読出し勾配とからなり、第２のパルス及び勾配シーケ
    ンスは少なくとも第１のパルス及び勾配シーケンスから
    ずれる励起パルスと、位相符号化勾配及び読出し勾配と
    からなり、別な勾配は第１のパルス及び勾配シーケンス
    の位相符号化パルスに重畳され、また読出し勾配中に印
    加され、第２のパルス及び勾配シーケンスの位相符号化
    勾配に読出し勾配中に印加される逆の別勾配が重畳され
    ることを特徴とする請求項１記載の磁気共鳴影像方法。
  8. 【請求項８】 第１のパルス及び勾配シーケンスは、励
    起パルスと、勾配方向が一致しないパルス符号化勾配及
    び読出し勾配とからなり、第２のパルス及び勾配シーケ
    ンスは又少なくとも第１のパルス及び勾配シーケンスか
    らずれる励起パルスと、位相符号化勾配及び読出し勾配
    とからなり、読出し勾配と一致する更なる勾配は読出し
    勾配の勾配からずれた勾配方向を印加され、３次元影像
    の場合に、更なる位相符号化勾配が第１及び第２のパル
    ス及び勾配シーケンスに印加され、第３のパルス及び勾
    配シーケンスは推定及び補正用の第３のデータセットを
    得るため印加され、励起パルスと磁気共鳴信号から読出
    す間アクティブである３つの勾配とからなり、その勾配
    方向は一致していないことを特徴とする請求項１記載の
    磁気共鳴影像方法。
  9. 【請求項９】 データラインの交差の点のデータの位相
    はデータの組み合せとしてとられ、位相は第１及び第２
    のデータセットの影像エラー源の推定として該組み合せ
    から順次に推定され、データラインのデータは推定位相
    に基づいて補正され、３次元影像の場合、第３のデータ
    セットの位相が推定され、データラインのデータが推定
    位相に基づいて補正されることを特徴とする請求項８記
    載の磁気共鳴影像方法。
  10. 【請求項１０】 ２次元影像の場合、第２のデータセッ
    トは１つのデータラインからなり、３次元影像の場合、
    第３のデータセットは１つのデータラインからなること
    を特徴とする請求項８又は９記載の磁気共鳴影像方法。
  11. 【請求項１１】 ２次元影像の場合、第２のデータセッ
    トはいくつかのデータラインからなり、推定はその位相
    を平均化することで決められ、３次元影像の場合、第３
    のデータセットは数データラインからなり、推定はその
    位相を平均化することにより決められることを特徴とす
    る請求項８又は９記載の磁気共鳴影像方法。
  12. 【請求項１２】 推定は繰返して決められることを特徴
    とする請求項１乃至７のうちいずれか一項記載の磁気共
    鳴影像方法。
  13. 【請求項１３】 推定は雑音に関して重み付けされるこ
    とを特徴とする請求項１乃至１２のうちいずれか一項記
    載の磁気共鳴影像方法。
  14. 【請求項１４】 安定磁界に配置された対象からの磁気
    共鳴信号から形成された磁気共鳴影像の影像エラーを減
    少させ、安定磁界を発生する手段と、ＲＦ電磁パルスを
    発生する手段と、安定磁界に重畳された勾配磁界用手段
    とからなる磁気共鳴影像装置であって、対象から磁気共
    鳴信号を少なくとも第１及び第２のパルス及び勾配シー
    ケンスにより得るためのプログラムされた手段からな
    り、該信号の信号サンプルは少なくとも第１及び第２の
    データセットに蓄積され、各データセットのデータライ
    ンは周波数領域で互いに交差し、交差データラインのデ
    ータの組合せから影像エラー源の推定を決定し、決めら
    れた推定に基いて少なくとも１つのデータセットのデー
    タラインのデータを順次補正し、その後磁気共鳴影像が
    積分変換により補正されたデータセットから決められる
    のに適しているプログラム手段を特徴とする磁気共鳴影
    像装置。
JP3335331A 1990-12-21 1991-12-18 磁気共鳴影像方法及び装置 Pending JPH057571A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
NL9002842 1990-12-21
NL9002842A NL9002842A (nl) 1990-12-21 1990-12-21 Magnetische resonantie werkwijze en inrichting ter reductie van beeldfouten in een magnetisch resonantie beeld.

Publications (1)

Publication Number Publication Date
JPH057571A true JPH057571A (ja) 1993-01-19

Family

ID=19858192

Family Applications (1)

Application Number Title Priority Date Filing Date
JP3335331A Pending JPH057571A (ja) 1990-12-21 1991-12-18 磁気共鳴影像方法及び装置

Country Status (5)

Country Link
US (1) US5262725A (ja)
EP (1) EP0492706B1 (ja)
JP (1) JPH057571A (ja)
DE (1) DE69125123T2 (ja)
NL (1) NL9002842A (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20020033550A (ko) * 2000-10-30 2002-05-07 추후제출 Ｍｒ 촬상 방법, 위상 에러 측정 방법 및 ｍｒｉ 장치
KR101036046B1 (ko) * 2008-07-22 2011-05-20 가천의과학대학교 산학협력단 초고자장 ｍｒｉ 장비에서 ｅｐｉ 영상 왜곡 보정 방법

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5517122A (en) * 1993-11-04 1996-05-14 The Regents Of University Of California T2 restoration and noise suppression of hybrid MR images using Wiener and linear prediction techniques
DE4343023C1 (de) * 1993-12-16 1995-03-09 Bruker Medizintech Magnetresonanz-Bildgebung mit kombiniertem Rückprojektions- und Fouriertransformationsverfahren
DE4416363C2 (de) * 1994-05-09 1996-03-21 Siemens Ag Rekonstruktion von Bildern aus MR-Signalen in inhomogenen Magnetfeldern
DE4427496C2 (de) * 1994-08-03 1999-02-18 Siemens Ag Verfahren zum Betrieb eines Kernspintomographiegerätes zur Gewinnung mindestens zweier unterschiedlich gewichteter Bilder
US5785042A (en) * 1995-02-14 1998-07-28 Duke University Magnetic resonance imaging method providing for correction of striation artifacts
JP3505294B2 (ja) * 1995-03-28 2004-03-08 ジーイー横河メディカルシステム株式会社 Ｍｒｉ装置
US5672969A (en) * 1996-04-25 1997-09-30 General Electric Company Reduction of Nyquist ghost artifacts in oblique echo planar imaging
JP3384944B2 (ja) * 1996-07-11 2003-03-10 ジーイー横河メディカルシステム株式会社 Ｍｒｉ装置
JP4191884B2 (ja) * 2000-08-18 2008-12-03 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 画像処理方法、画像処理装置および画像撮影装置
JP3720752B2 (ja) * 2001-10-26 2005-11-30 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー ０次位相検出方法およびｍｒｉ装置
US7511495B2 (en) * 2005-04-25 2009-03-31 University Of Utah Systems and methods for image reconstruction of sensitivity encoded MRI data
US8085041B2 (en) * 2008-04-10 2011-12-27 General Electric Company Three-point method and system for fast and robust field mapping for EPI geometric distortion correction
JP5638393B2 (ja) * 2008-10-20 2014-12-10 株式会社日立メディコ 磁気共鳴イメージング装置及び方法
US8405392B2 (en) * 2008-12-19 2013-03-26 Sloan Kettering Institute For Cancer Research Corrected nuclear magnetic resonance imaging using magnetization transfer
US8427147B2 (en) * 2009-04-10 2013-04-23 Wisconsin Alumni Research Foundation Magnetic resonance imaging with fat suppression by combining phase rotating data with phase shifted data in K-space
FR2951835B1 (fr) * 2009-10-26 2013-10-18 Bruker Biospin Dispositif de correction de signaux de consigne et systeme de generation de gradients comportant un tel dispositif
US8653816B2 (en) * 2009-11-04 2014-02-18 International Business Machines Corporation Physical motion information capturing of a subject during magnetic resonce imaging automatically motion corrected by the magnetic resonance system
JP5611882B2 (ja) * 2010-05-31 2014-10-22 株式会社東芝 磁気共鳴イメージング装置
JP5925529B2 (ja) * 2011-03-31 2016-05-25 株式会社東芝 磁気共鳴イメージング装置
US9475930B2 (en) 2012-08-17 2016-10-25 Metabolix, Inc. Biobased rubber modifiers for polymer blends
EP3004225A1 (en) 2013-05-30 2016-04-13 Metabolix, Inc. Recyclate blends
DE112016001879T5 (de) * 2015-04-24 2018-01-04 Koninklijke Philips N.V. Ein Mehrkanaliges Sende-/Empfangs Hochfrequenz (HF) System
EP3422037A1 (en) * 2017-06-27 2019-01-02 Koninklijke Philips N.V. Method and device for determining a motion field from k-space data
CN112156383B (zh) * 2020-08-31 2022-10-11 上海沈德医疗器械科技有限公司 一种基于k空间能谱分析的磁共振温度修正方法及系统

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6117054A (ja) * 1984-07-03 1986-01-25 Yokogawa Medical Syst Ltd 核磁気共鳴断層撮影装置
DE3434161A1 (de) * 1984-09-18 1986-03-27 Bruker Medizintechnik Gmbh, 7512 Rheinstetten Verfahren zum messen der magnetischen kernresonanz
US4706026A (en) * 1985-08-16 1987-11-10 General Electric Company Method for reducing image artifacts due to periodic signal variations in NMR imaging
JPS62231647A (ja) * 1986-03-31 1987-10-12 株式会社東芝 磁気共鳴映像装置
NL8602019A (nl) * 1986-08-07 1988-03-01 Philips Nv Magnetische resonantiewerkwijze en -inrichting voor het elimineren van fasefouten in beeldelementen van een complexe afbeelding van een kernmagnetisatieverdeling.
US4833407A (en) * 1987-06-24 1989-05-23 Picker International, Inc. Scan time reduction using conjugate symmetry and recalled echo
DE3823398A1 (de) * 1988-07-09 1990-01-11 Spectrospin Ag Verfahren zur erzeugung einer folge von spinechosignalen, die verwendung dieses verfahrens bei der kernspintomographie und vorrichtung zum durchfuehren bzw. zur verwendung dieses verfahrens
EP0351446B1 (de) * 1988-07-22 1992-04-01 Siemens Aktiengesellschaft Pulssequenz für ein Kernspin-Tomographiegerät
JPH0245037A (ja) * 1988-08-08 1990-02-15 Toshiba Corp 磁気共鳴イメージング装置
DE4003547C2 (de) * 1989-02-24 1999-10-07 Siemens Ag Abtastung von Kernresonanzsignalen bei allgemeinen Gradientenformen

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20020033550A (ko) * 2000-10-30 2002-05-07 추후제출 Ｍｒ 촬상 방법, 위상 에러 측정 방법 및 ｍｒｉ 장치
KR101036046B1 (ko) * 2008-07-22 2011-05-20 가천의과학대학교 산학협력단 초고자장 ｍｒｉ 장비에서 ｅｐｉ 영상 왜곡 보정 방법

Also Published As

Publication number Publication date
NL9002842A (nl) 1992-07-16
DE69125123D1 (de) 1997-04-17
EP0492706B1 (en) 1997-03-12
EP0492706A1 (en) 1992-07-01
US5262725A (en) 1993-11-16
DE69125123T2 (de) 1997-08-21

Similar Documents

Publication Publication Date Title
JPH057571A (ja) 磁気共鳴影像方法及び装置
CN105308469B (zh) 一种mr成像方法、mr设备以及相关数据载体
US8148984B2 (en) Method for magnitude constrained phase contrast magnetic resonance imaging
EP2626718A1 (en) MRI with motion correction using navigators acquired using a Dixon technique
EP3635425B1 (en) Mr imaging using a stack-of-stars acquisition with variable contrast
JP3970371B2 (ja) Ｍｒｉシーケンス
US20170307716A1 (en) Propeller mr imaging with artefact suppression
US10996301B2 (en) Dual-echo dixon-type water/fat separation MR imaging
US10996305B2 (en) Propeller MR imaging with relaxation weightings varying between k-space blades and k-space data lines
US10859652B2 (en) MR imaging with dixon-type water/fat separation
US10746832B2 (en) DTI with correction of motion-induced diffusion gradient inconsistency
US11543482B2 (en) Magnetic resonance imaging using motion-compensated image reconstruction
US11226385B2 (en) Dixon type water/fat separation MR imaging with improved fat shift correction
JP7128972B2 (ja) 歪み補正を伴うｅｐｉ ｍｒ画像
EP2392935A1 (en) EPI distortion correction using non-phase encoded reference echoes
JP2024506575A (ja) Ｄｉｘｏｎタイプ水／脂肪分離ＭＲ画像