JP2009000538A - Nmrイメージング操作中に誘導される渦電流の測定及び補償のための装置 - Google Patents

Nmrイメージング操作中に誘導される渦電流の測定及び補償のための装置 Download PDF

Info

Publication number
JP2009000538A
JP2009000538A JP2008186366A JP2008186366A JP2009000538A JP 2009000538 A JP2009000538 A JP 2009000538A JP 2008186366 A JP2008186366 A JP 2008186366A JP 2008186366 A JP2008186366 A JP 2008186366A JP 2009000538 A JP2009000538 A JP 2009000538A
Authority
JP
Japan
Prior art keywords
magnetic field
gradient
signal
field gradient
spin
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
JP2008186366A
Other languages
English (en)
Inventor
Weiguo Zhang
ウェイグオ・ツァン
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.)
Toshiba America MRI Inc
Original Assignee
Toshiba America MRI Inc
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 Toshiba America MRI Inc filed Critical Toshiba America MRI Inc
Publication of JP2009000538A publication Critical patent/JP2009000538A/ja
Pending legal-status Critical Current

Links

Images

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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56518Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field

Abstract

【課題】EC測定又はMRIシーケンス中に生ずる大きなバックグランド静磁場勾配を評価し、補償するための装置を提供する。
【解決手段】Bo磁場の不均一性によりEC測定域内にあるバックグランド磁場勾配、及びMRI装置内の他の静磁場勾配の定量値をも決定する手段と、バックグランド磁場勾配及び他の静磁場勾配に対応する方向に沿って同じマグニチュードの小さい反対の磁場勾配を加える手段とを含む装置により補償される。
【選択図】図2

Description

発明の分野
本発明は、一般的には、核磁気共鳴(NMR)イメージング技術に関する。更に詳細には、本発明は、磁気共鳴イメージング(MRI)システムにおける磁場勾配スイッチングの結果として誘導された渦電流の定量的(quantitatively)測定及び補償装置に関する。
発明の背景
磁気共鳴分光学(MRS)にその基礎を有している磁気共鳴イメージング(MRI)は、核磁気共鳴(NMR)現象に敏感である相当数の原子核を有している物体(人間の身体のような)の内部構造を表わすディジタル化した目に見える画像を得るために広く受け入れられ、且つ商業的に利用可能な技術となった。MRIでは、映像化されるべき体内の核は、強力な主磁場Hoを核にかけることによって、極性を与えられる。選択された核は、特定のNMR周波数の無線周波数(RF)信号を加えることによって励起される。局部に集中した磁場を空間的に分散し、それから、核から生じたRF応答を適切に分析することによって、核の位置の関数として相対的にNMR応答のマップ又は画像が決定できる。フーリエ分析につづいて、空間のNMR応答を表わすデータがCRT上に表示できる。
図1に示したように、NMRイメージングシステムは典型的に、静磁場を加える磁石10と、3つの直角座標に沿って空間的に分布した磁場を加えるための勾配磁場コイル14と、選択した核に及び核からRF信号を送受信するためのRFコイル15及び16とを含む。コイル16によって受信されたNMR信号は、データをディスプレイ24上に表示される画像に処理するコンピュータ19に送られる。磁気共鳴画像は“ピクセル”と呼ばれる画素で構成されている。ピクセルの強さは、画像に表わされる物体の対応するボリューム要素又は“ボクセル”の内容のNMR信号強さに比例する。コンピュータ19はまた、それぞれRF増幅器/検出器21と22及び勾配増幅器20を介してRFコイル15と16及び勾配磁場コイル14の作動を制御する。
奇数のプロトン(陽子)及び/又は中性子を有している核のみが磁気モーメントを有しており、従ってNMR現象を受け易い。MRIでは、強力な静磁場核を整合するのに使用されて、平衡状態で主磁場に平行に整合した全磁化ベクトルを発生する。単一のRFパルスとして第1の磁場を横切って加えられる第2の磁場が、エネルギーを核に送り、そしてそれが全磁化ベクトルを例えば90°だけフリップさせる。ある例では、このRFパルスはまた、第2の強力RFパルスに続いていて、全磁化ベクトルを更に180°フリップさせる(時折、“フロップ”と呼ばれる)。励起後、核は歳差運動をし(precess)そして、次第に静磁場との整合に戻り、緩和する。核が歳差運動し緩和するにつれて、核は自由誘導減衰(FID)として知られている、弱いが検出可能な電気エネルギーを周りのコイルに誘導する。集合的にここではMR信号と呼ぶ、FID信号(及び/又はその磁気勾配−再焦点磁場エコー)は、空間に核の画像を作るため、コンピュータによって分析される。
MR信号を生じ、且つMR信号を得るため、種々のコイルがRF励起パルス及び勾配磁場を作り出す作用は、取得パルスシーケンス(acquisition pulse sequence)と呼ばれている。簡単な取得シーケンスのグラフ図が図2に示されている。この基本的実例では、勾配磁場Tgradは核のスライスを特定のRF周波数に敏感にするため、最初に主磁場に沿って付加される。その後、概略的に言うと、180°RFパルスに続く“エコー”のように現われるMR信号を発生するため、特定の核のスライスが、90°励起パルスによって選択され、180°RFパルスが続いて印加される。180°RFパルス中心と検出したMR信号中心との間の時間が、90°RFパルス中心と180°RFパルス中心との間の時間と同じであれば、作られたMR信号は、通常対称な“スピン−エコー”信号と呼ばれ、そして図2に示された加えたパルス及び磁場の特定のタイミングは、スピン−エコーシーケンス(spin−echo sequence)として知られている。更に、MRエコー信号中心は、加えた磁場勾配を調整することによりシフトできるので、非対称スピン−エコー(ASE)がその代わり作り出される。
NMR周波数及び主Bo磁場は、ラーモア関係によって関連づけられる。下記の式によって表わされたこの関係は、核の歳差運動の角周波数ωが、磁場Boと所謂、磁気回転比γとの積であり、各核種に対する基礎的物理定数であることを示している。
Figure 2009000538
磁気共鳴イメージング(MRI)及びボリューム−局所的(volume−localized)磁気共鳴分光学(MRS)では、磁場勾配のスイッチングが、殆んどすべてのパルスシーケンスの最も肝要な部分である。そのような作用が、磁石内の周りの伝導性材料内に渦電流(EC)を誘導し、そして、望ましくないEC磁場勾配及びBo(主磁場)振動を生じ、そしてそれ等が、最終的に画像のアーチファクトを生ずることはよく知られた事実である。通常、いくつかのMRI及びMRSシステムは、勾配磁場に誘導された渦電流の効果を最小にしようとして、システムの整合/較正手順中に調整される“プレエンファシス”ネットワークを使用する。残念乍ら、磁場に誘導された渦電流の精確な測定なしでは、勾配磁場に誘導された渦電流を効果的に最小にするためプレエンファシスネットワークを調整するのは難かしい。
渦電流を直接測定するため、小さいピックアップコイルを使用できることが一般的に知られているが、渦電流測定のためのピックアップコイルの使用は、典型的に、追加のハードウエアを必要とし、且つその使用に専門的知識を必要とする。例えば、J. Phys. E18, 1014 (1985), J.チャンキジイ, J.レフェブレ及びA.ブリグエットによる“電流コイルによって発生した磁場のマッピング方法”参照。例えば、磁場勾配パルスを切った後、小さいRFコイルを使用して、小さい試料の自由誘導減衰(FID)信号を測定することによるような、渦電流を測定するためのNMR測定技術に依存している他の多少、直接的でない方法もまた知られている。例えば、J. Phys. E19, 708 (1986), E.山本, H.河野による“NMRによる勾配時間−形状測定”;Magn. Reson. Med 29, 119−121 (1993), R.E.ワイソング及びI.J.ローウエによる“補償回路設定のため、渦電流に誘導された勾配の簡単な測定方法”;Magn. Reson. Med. 31, 73−76 (1994), Q.リュー, D.G.ヒューズP.S.アレンによる“40cmボア(穴径)の超電導性磁石内の渦電流磁場の定量特性”を参照。いくつかのNMR渦電流測定方法もまた、全試料を使用する。例えば、Magn. Reson. Med. 20, 268−284 (1991), Ch.ボッシュ, R.グリュッター及びE.マーチンによる“渦電流により発生した磁場の時間的及び空間的分析”参照。残念乍ら、すべての上記の公知の方法は、実施が非常に困難であるか、及び/又はEC磁場勾配とBo(主磁場)振動との間を区別する能力に欠けている−双方とも渦電流に関連しているが、画質にかなり異なる効果を有している。
故障防止システム実施問題のため、又は渦電流の悪影響を最小にするためのプレエンファシス装置の調節を助けるための特殊な装置を必要とせずに、磁場勾配スイッチングにより誘導された渦電流を、都合よく測定する能力を有することは、MRIシステムにとって、極めて有用である。従って、大きい及び小さい磁場不均一性が存在するとき、EC誘導磁場勾配及びBo振動の明確な定量(quantative)測定値を得るため装置の必要性がある。
発明の概要
本発明は、NMRイメージング操作中に誘導される渦電流の効果を測定し、補償するための装置を提供する。立方体又は円筒状の試料が、主磁場に関して中央にある位置に、且つ所望の測定方向に平行なその長手方向の軸線に向けて、MRIシステムのイメージングボリューム内に置かれる。磁場勾配パルスが、渦電流誘導のため、並びに、スライス選択スピン−エコー信号の発生のために加えられる。スピン−エコー信号は、勾配パルスを誘導する各渦電流の終結直後に得られる。2つのスライスが、主磁場の中心から同じ対称距離の所望の測定方向に沿って選択される。2つのスピン−エコー信号は、各スライスに対し;2つのエコー信号間で極性を反転した渦電流誘導勾配パルスにより得られる。磁場勾配及びBo振動を誘導する渦電流の定量値(Quantitative value)はそれから、得られたNMR信号の歳差運動周波数に基づいて決定される。
画質は、適切なプレエンファシスネットワークを設定するため、磁場勾配及びBo振動の定量値を加えることによって、渦電流効果を補償することにより改善される。更に、MRI/MRS取得シーケンス(acquisition sequence)内の勾配パルスは、更に勾配−誘導渦電流を補償するようにバランスしてもよい。開示された発明の他の見地は、勾配スイッチングの時間経過の測定、長い時定数を有している渦電流を測定するため、信号取得シーケンスの変更、プレエンファシス調整を助けるため、渦電流の繰返し測定、及び適度に大きいバックグランド磁場(主Bo磁場)不均一性の存在するとき、EC−誘導磁場勾配及びBo振動の測定を含む。
詳細な説明
本発明によって得られる目的及び利点は、特に添付図面を参照して、現在好ましい実施例の下記の詳細な説明を注意深く検討することにより理解されるであろう。
図1はMRIシステムを示している。そのようなシステムの1つの実例が、東芝OPARTTMMRIシステムである。そのようなMRIシステムは、患者イメージングボリューム11内に、実質的に均等、均一な分極磁場Boを発生する大型分極磁石構造体10を含む。適切なキャリッジ12が患者13の体の所望の部分をイメージボリューム11内に挿入する。磁場勾配が、電磁勾配磁場コイル14によって選択的に作り出される。RF核章動(nutation)パルスが、RFコイル15によってイメージボリューム内で患者の組織に伝えられる。MR信号を構成するRF応答が、適切なRF検出コイル構造体16を経て患者の組織から受けとられる。
MRIデータを得るためMRIシステムは、プログラム可能なコンピュータ/プロセッサ19の制御の下で、MRIパルスシーケンス、コントローラ17及び18を経て磁場勾配及びRF章動パルスを発生する。更に、プロセッサ19が勾配パルス増幅器20及びRF電源及び増幅器回路21及び22を制御する。MR信号(RF検出器)回路は、シールドしたMRIシステムガントリ内に置かれたMR信号RFコイルと適切にインタフェースされる。受信したMR応答は、ディジタル化装置によってディジタル化され、そしてイメージ処置のためのプロセッサアレイ等及び適切なコンピュータプログラム記憶媒体(図示せず)を含むプロセッサ19に送られる、その場合に、プログラムは記憶され、そしてMR信号データの取得及び処理を制御し、且つ制御端末24のCRT上にイメージ表示を作り出すために選択的に利用される。MRIシステムは、イメージングシーケンスコントローラ17及び18を通してオペレータ制御を働かせるための適切なキーボードスイッチ等を含む制御端末24を備えている。イメージ(画像)はまた、プリンティング装置25により、直接フィルム又は他の適切な媒体上に記録される。
システムコンピュータ/プロセッサ19に関連して、オペレータは典型的に、MRIシーケンス及びデータ処理技術のための選択メニューを提供される。本発明の実施例では、オペレータに利用できるこれ等の選択の1つは、MRI装置内の磁場勾配のスイッチングによって生ずる渦電流を測定し、分析するためのNMR信号データを得るため、本発明によってMRIパルスシーケンスを提供するプログラムである。MRI装置を制御し、且つ本発明の上記の方法を行なうための適した詳細なコンピュータプログラムの生成は、本開示の完全さから見て、十分当業者の能力範囲内にあると信じられる。
EC磁場勾配及びBo振動の定量(Quantitation)
本発明の第1の実施例では、スライス選択勾配を有している簡単なスピン−エコーシーケンスが、短時間−一定(持続)渦電流測定に使用される。この実例では、スライス選択に使用される磁場勾配パルスもまた、EC誘導勾配パルスでもあり、そのために渦電流が測定される。他の実施例は、異なるEC−誘導勾配パルスを利用している。図2は、短時間−一定渦電流を測定するのに使用されるスピン−エコーシーケンスの実例を示している。好ましくは、短時間−一定渦電流を測定するとき、スピン−エコー信号は、持続時間TgradのEC−誘導勾配パルスを切った直後にサンプルされ、そして、EC磁場勾配及びBo振動の定量が要求される期間Tacq中に得られる。
本発明によれば、立方体又は長方形寸法を有しているMR活性被験物体(試料)が、測定方向に沿って向けられた長手方向の軸線を有している主Bo磁場の中心に位置づけされる。好ましくは、合計4つのスピン−エコー信号が、測定方向に沿って、主磁場の中心から距離△L離れた位置に対称に置かれた2つのスライスから得られる。勾配の極性は、得られる第1の2つのスピン−エコー信号と得られる第2の2つの信号との間に反転される。磁場勾配の大きさ(マグニチュード)及び2つのスライスの分離/厚さはオペレータ制御を経て調整可能であるか、又は一定の適切な値に設定できる。
本発明の実施例の下記の論述を簡易化する目的のために、勾配スイッチングによって発生した渦電流の一次効果が、誘導勾配パルスの同じ方向に沿ってEC磁場勾配を誘起するだけでなく、ある量のBo振動を誘導すると仮定する。更に、EC磁場勾配及びBo振動の極性が、勾配パルスを誘導する渦電流の極性が反転したとき、反転すると仮定する。本発明による測定は、試料空間の三次元のどこに対しても行なうことができるが、本論述は、単一X, Y, Z座標方向に沿う場合に、−例えば、主磁場の“Z”方向に限定される。当業者はこのようにして得られたいかなる結果も、残りの2つの空間次元に容易に拡大できることは理解されるであろう。
従って、上記のようにして得られた4つのNMR信号は数学的に下記のように表わすことができる:
Figure 2009000538
この場合、Sn(t)は、得られたMR信号;△Zは、磁石の中心からのスライスのオフセット;ρ(X,Y,△Z)及びρ(X,Y,Z)は、それぞれスライス△Z及び−△Zのスピン分布;△B(X,Y,△Z)及び△B(X,Y,−△Z)は、それぞれスライス△Z及び−△Zのバックグラウンド磁場異質性であり;△B(t’)は時間変化(time−varying)EC誘導Bo振動であり;そして△G(t’)は時間変化EC磁場勾配である。
下式であるとすれば、
Figure 2009000538
当業者は、4つのスピン−エコー信号の時間変化スピン歳差運動周波数が、下記の4つの式により、時間変化EC誘導磁場勾配及びBo振動に関係があることは理解されるであろう。
Figure 2009000538
この場合fn(t)は、スピン−エコー信号の時間変化歳差運動周波数であり;△G(t)及び△B(t)は、それぞれ、誘導渦電流によって作り出された磁場勾配及びBo振動であり;γは、磁気回転比であり;そしてφは、選択したスライス内の磁場不均一性によって誘導された位相角である。
従って、時間変化EC磁場勾配、△G(t)及びBo振動、△B(t)は、下記の2式によって定量化(quantitate)できる。
Figure 2009000538
更に、上記の時間変化“歳差運動周波数”は、下記の関係によって4つのNMR信号から定量的に得ることができる:
Figure 2009000538
この場合、n=1,2,3,4のときのfn(i)は、第1及び第2の対のスピン−エコー信号に対応している4つのNMR信号の歳差運動周波数を表わしており;i=0,1,2 −−− N−2のときのiは、Nが、スピン−エコーNMR信号から得た所定数のデータ試料であるときの特定のデータ試料を表わしている;Su n及びSv nは、得られた複合(complex)NMR信号データの実部及び虚部であり;U(i)及びV(i)(u(i)ドット及びv(i)ドットを意味する)は、それぞれu(i)及びv(i)に対応する時間導関数であり;そして△tは、得られたNMR信号データの“ティック(tick)”又は、サンプリング率(秒)である。
従って、EC磁場勾配及びBo振動に対する定量的値は、[数8]を用いて4つの歳差運動周波数に対して得られた定量的値から決定される。
定量の実例として、図3-6及び7-10は、本発明の装置によって得られた4つのNMR信号の歳差運動周波数を示している。この実例では、通常の市販の“ベイビイオイル(baby−oil)から成っている10cm立方体形状の試料が、上記のように本発明の装置によるEC測定値を得るためのMR活性(active)被験体として使用される。2mmの厚さのスライスが使用され、そしてスライスは、主磁石(即ち、主Bo磁場)の中心から4cmの距離で得られた、従って2つのスライス間に8cmの距離が与えられる。振幅0.95ガウス/cm及び持続時間100msを有している勾配パルスTgradが使用される。4つのスピン−エコー信号は、取得時間(acquisition period) 20.48msに亘り、512の複合(complex)サンプリング点を用いて得られた。従って、連続繰返し時間(周期)TR=0.5秒及び信号エコー時間TE=56msのとき4つの平均値が得られ、従って全セットのデータに対する合計取得時間は4秒となる。第1のセットの4つのデータ信号は、調整したプレエンファシスネットワークを有しているMRIシステムを用いて得られた。第2のセットの4つのデータ信号は、プレエンファンシスネットワーク“調整誤り(mis−adjusted)”のとき、同じシステムを用いて得られる。
プレエンファシスネットワークは、渦電流誘導勾配磁場を最小にするため、MRIシステム調整中に調整される各物理的チャンネルに対する3つの時定数を利用する。“調整誤り”は、簡単にプレエンファシスネットワークをスイッチオフすることによって達成された。従って、図3−6は、予め調整されたプレエンファシスネットワークで得られた4つのNMR信号の歳差運動周波数を示しており、そして図7−10は、故意に調整を間違えたプレエンファシスにより得られた対応する信号データを示している。
図11及び12は、[数8]を用いて決定したEC磁場勾配の時間経過及び図3−6と図7−10に示された歳差運動周波数データを例示している。従って図11は、正しく調整されたプレエンファシスのときのEC磁場勾配の時間経過を例示し、そして図12はプレエンファシスを切ったときのEC磁場勾配の時間経過を例示している。図11及び12では、EC磁場勾配は、1000の誘導勾配マグニチュードにより表わされている。
スイッチングの勾配遷移特性
本発明の他の実施例によれば、EC勾配を定量するための上記の装置もまた、NMRイメージング装置における勾配スイッチング中(即ち、勾配マグニチュードの急激変化)の時間経過を特徴づけるのに使用された。本発明により作り出されるように、勾配スイッチングの時間経過の勾配強さ対時間特性グラフは、MRIハードウエアの機能評価にも、MRI/MRSパルスシーケンスにおける勾配パルスの自動バランシングにも有用であることが証明できる。磁場対時間特性データは、上述のようにスピン−エコー信号シーケンスを使用し、且つ、勾配磁場がスイッチされる直前又は直後にスピン−エコー信号データのサンプリングを開始し、そして、スピン−エコー信号200サンプリング時間Tacqによって、図2に例示したように、勾配スイッチングの全遷移時間に亘ってサンプリングを続け、且つ勾配パルスTgradのエッジ201を下降することによって得られる。
図13及び14は、本発明の装置を用いて得られたNMRイメージング装置の勾配スイッチングの時間経過グラフの実例を示している。この実例では、X座標方向(即ち、主磁場の従来のZ座標方向に垂直)の磁場勾配は、0.95ガウス/cmから−0.95ガウス/cmに“スイッチ”される。図13に示したように、勾配遷移の時間経過は、勾配磁場がスイッチされた直後に始まり、2,048msの期間に亘りMRスピン−エコー信号をサンプリングすることによって得られる。所望により、勾配スイッチング測定からの時間経過データは、多少より好都合な勾配スイッチング“モデル”に“フィット”する。例えば、1つの好都合なモデルは、遅れ時間及び、スルーレートによって簡単に特性を表わされるモデルである。そのようなモデルによれば、勾配スイッチングは、図14の実線で例示されているように、遷移の最終磁場マグニチュードに対する一定のスリューレートにおけるマグニチュードの直線的変化につづく短時間の間持続する一定値の磁場マグニチュードにより表わされる。
図15は、本発明の装置を用いて測定したMRI装置における勾配スイッチングの時間経過グラフの第2の実例を示している。この実例では、信号エコー時間、TE,20msを有しているスピン−エコーシーケンスは、1ガウス/cmから0ガウス/cmに“スイッチされた”125msの長い磁場勾配パルスに使用される。MR信号データは、512データサンプリング点を用いて勾配遷移中に、1試料当り4μsのサンプリング速さで得られる。
長い時定数を有している渦電流の測定
勾配スイッチングの発生後、測定値を得るための時間のウィンドウ(window of time)は、“T2”又は、試料の見かけのスピン−スピン緩和時間として知られているものによって制限されるので、図2のスピン−エコーシーケンスは、長い時定数を有している渦電流の測定には有用ではない。従って、本発明の他の実施例によって、長い時定数を有している渦電流は、図16に示されたような異なるセットのNMRパルスシーケンスを用いて測定される。図2に示されたように、短い時定数の渦電流を評価するのに使用される特定のスピン−エコーシーケンスと異なり、長い時定数の過電流を測定するのに好ましいパルスシーケンスのタイプは、図16に示された3つの実例によって例示されているように、EC誘導勾配パルスをスイッチオフするとき起る勾配パルス遷移周期の終結後発生されるNMR信号を作り出さなければならない。NMR信号データ試料取得“ウィンドウ”は、従って長い時定数渦電流を測定するのに必要な、長い信号データ取得周期をカバーするため、EC勾配パルスの終結後、所定時間生ずるようにフレキシブルに配置できる。図16は、選択した実例として、それぞれ(a),(b)及び(c)で、スピンエコー、磁場エコー及び刺激されたエコー信号取得パルスシーケンスを例示しているが、他のスライス選択取得システムもまた、長い時定数を有している渦電流を測定するのに利用されてもよい。
NMR信号データは、上述の2スライス、4−信号サンプリング手順を用いることによって得られ、そしてEC磁場勾配△G(t)及びBo振動△B(t)に対する定量値は、[数8]を用いて決定される。必要があれば、選択したパルスシーケンスは、多数のサンプリングウィンドウを経て定量を達成するため、勾配スイッチング後、何回繰返されてもよい。スライス−選択勾配は、EC勾配が各スライスの2つのエコー信号間で極性を反転される間同じであるから、スライス選択勾配パルスによって誘導されたいかなる渦電流も、[数8]を用いて定量されるように、EC磁場勾配△G(t)及びBo振動△B(t)から除去される。
図17及び18は、本発明により磁場−エコーパルスシーケンスを用いて測定した長い時定数を有しているEC勾配の実例の時間経過グラフである。図17は、プレエンファシスオフのときに作られたEC勾配を示しており、そして図18は、プレエンファシスオンのとき作られ、そして正しく調整されたEC勾配を示している。この実例では、一連の512スライス選択磁場−エコー信号が0.75ガウス/cmの1秒の長い勾配パルスの後に得られる。磁場エコー信号は、40μsサンプリングピッチで、32データ点のとき20ms毎に得られる。従って、この実例では、長い時定数のEC勾配測定のためのデータ取得時間ウィンドウは20msの時間分解能(time resolution)のとき10.24秒をカバーする。各信号の32データ点は、各信号に対する歳差運動周波数を得るのに使用される。それから、歳差運動周波数が各々の得られた信号を表わす単一の周波数に平均され、そしてEC勾配及びBo振動が、それから上述のように、4つのスライス選択測定値からの平均歳差運動周波数から決定される。
プレエンファシス調整のための渦電流の繰返し測定及び表示
本発明の上述の装置による渦電流の測定は繰返して行なわれ、そして測定したEC勾配及びBo振動の時間経過が、MRIシステムに関連したプレエンファシス調整に応答して特定のMRIシステムを評価するため表示される。例えば、信号を平均せずに、0.2秒のパルスシーケンス繰返し周期(TR)を用いて、0.4秒の時間分解能が単一のTR周期に2つのスライスを得ることにより、上記の短い時定数EC測定装置を用いて達成できる。EC勾配及びBo振動は、別々に定量されるので、プレ−エンファシスネットワークのBo及び勾配成分の調整は別々に行なわれ、監視される。
大きいバックグランド磁場の不均一性の定量及び補償、
不均一なバックグランド磁場が、下記の関係によって数学的に表わされる支配的に大きな静直線勾配を含んでいると仮定すれば:
Figure 2009000538
そのとき、渦電流を測定し、そしてEC勾配及びBo振動を決定するための本発明の上述の4信号、2スライス法はまた、例えば下式によりZ方向に対するバックグランド静勾配を評価するのに使用できる。
Figure 2009000538
渦電流を測定するための本発明の前述の実施例は、主に、他の磁場勾配の存在しないときに得られるスライス選択エコー信号の取得に依存している。特に、測定方向に垂直な方向に、大きなバックグランド磁場勾配又はバックグランド磁場の不均一性が存在するときは、得られた信号振幅が振動し、そしてEC測定値に誤差を生ずる。従って、本発明の他の実施例により、そのような誤差は、最初に静バックグランド勾配の決定し、それから同じマグニチュードであるが、対応する方向に沿って静勾配と反対の小さい磁場勾配の使用によって補償することにより修正される。
図19は、本発明の更に他の実施例による大きな直線バックグランド勾配を定量するのに使用されるパルスシーケンスの実例を示している。図19に示されたシーケンスは、基本的には測定方向に直角な方向に沿ったライン選択のときの、ライン選択非対称スピン−エコー(ASE)シーケンスである。エコー信号は、測定方向に沿って方向づけられた周波数−符号化(encoding)読み出し勾配で得られる。それ等の信号は2回測定される。1回は第1の(正の)方向の読み出し勾配で、そして2回目は、反対(負)の方向に向けられた読み出し勾配で測定される。
得られたMR信号の一次元(ID)フーリエ変換は、数学的に下記のように表わされる。
Figure 2009000538
この場合、ρ(Z)は、選択したラインのスピン分布である;Gz(Gzバーを意味する)は、測定又は読み出し方向に沿った静バックグランド磁場勾配である(この場合には、“Z”軸と呼ばれる);そしてφ(Z)は、読み出し勾配の調整ミスによる位相誤差である。
従って、信号位相は下記のように読み出しラインに沿った空間位置の関数として表わされる。
Figure 2009000538
従って、バックグランド磁場勾配は、下記の関係により上記の信号位相から決定される。
Figure 2009000538
この場合、ФおよびФ(φドットおよびφドットを意味する)は、選択したラインに沿った信号位相の空間的導関数である;そして△Tはスピン−エコー位置からのASE信号のタイムシフトである。
バックグランド磁場勾配の効果は、それから、同じマグニチュードを有しているが、バックグランド磁場勾配の極性とは反対の極性を有している追加の勾配によって、いくつかのイメージングシーケンスに補償される。そのような反対の勾配は、静磁場勾配の補償が必要な、取得シーケンスの任意の時間周期にイメージング勾配の頂部に適用さられる。図20−23は、短い時定数の渦電流測定のための前述のシーケンスを用いて得られたスピン−エコー信号を示している。図20及び21は、この実例では、“Y”座標方向に沿って主に、1.2ppm/cmの静勾配より成っている静勾配を補償せず得られたNMR信号を示している。図22及び23は、同じ静勾配を補償して得られたNMR信号を示している。図24及び25は、それぞれ、図21及び23のデータから、本発明によって定量された渦電流の時間経過を示している。
本発明は現在、最も実用的であり且つ好ましいと考えられているものに関連して記述されたが、本発明は、開示された実施例に限定されるべきではなく、添付の特許請求の範囲及び精神に含まれる種々の変形及び同等の装置を含む意図を有していると理解されるべきである。
渦電流の悪影響を最小にするためのプレエンファシス装置の調節を助けるための特殊な装置を必要とせずに、磁場勾配スイッチングにより誘導された渦電流を、都合よく測定することができる。従って、大きい及び小さい磁場不均一性が存在するとき、EC誘導磁場勾配及びBo振動の明確な定量測定値を得ることができる。
MRIシステムの概略図である。 勾配及びBo振動を誘導した渦電流を定量するための本発明の1実施例によるスピンーエコーシーケンスの実例のRF及び勾配波形タイミングダイアグラムを例示している。 正しく調整されたプレエンファシスで得られたNMR信号の歳差運動周波数を例示しているグラフである。 正しく調整されたプレエンファシスで得られたNMR信号の歳差運動周波数を例示しているグラフである。 正しく調整されたプレエンファシスで得られたNMR信号の歳差運動周波数を例示しているグラフである。 正しく調整されたプレエンファシスで得られたNMR信号の歳差運動周波数を例示しているグラフである。 故意に調整を誤ったプレエンファシスで得られたNMR信号の歳差運動を例示しているグラフである。 故意に調整を誤ったプレエンファシスで得られたNMR信号の歳差運動を例示しているグラフである。 故意に調整を誤ったプレエンファシスで得られたNMR信号の歳差運動を例示しているグラフである。 故意に調整を誤ったプレエンファシスで得られたNMR信号の歳差運動を例示しているグラフである。 図3−6に示されたデータから決定されたEC磁場勾配の時間経過のグラフである。 図7−10に示されたデータから決定されたEC磁場勾配の時間経過のグラフである。 本発明により測定したMRIシステムにおける勾配スイッチングの時間経過のグラフである。 本発明により測定したMRIシステムにおける勾配スイッチングの時間経過のグラフである。 本発明により測定したMRIシステムにおける第2の実例の勾配スイッチングの時間経過のグラフである。 長い時定数を有している渦電流測定のための他のMRIパルスシーケンスを例示している概略図である。 本発明の他の実施例により測定した、長い時定数を有しているEC勾配の実例の時間経過グラフである。 本発明の他の実施例により測定した、長い時定数を有しているEC勾配の実例の時間経過グラフである。 更に他の本発明の実施例による静バックグランド磁場勾配を測定するためのライン−選択非対称スピン−エコー(ASE)シーケンスの実例のダイアグラムである。 スピン−エコー信号の振幅及び位相のグラフである。 スピン−エコー信号の振幅及び位相のグラフである。 スピン−エコー信号の振幅及び位相のグラフである。 スピン−エコー信号の振幅及び位相のグラフである。 渦電流の時間経過のグラフである。 渦電流の時間経過のグラフである。
符号の説明
10…磁石 14…勾配磁場コイル 15…RFコイル 16…RFコイル 17…勾配パネルシーケンス 18…パルスシーケンス 19…制御コンピューター/プロセッサー 20…勾配増幅器 21…RF増幅器 22…RF増幅/検出器 23…デジタル化装置 25…プリンター

Claims (2)

  1. EC測定又はMRIシーケンス中に生ずる大きなバックグランド静磁場勾配を評価し、そして補償するための装置であって、
    Bo磁場の不均一性によりEC測定域内にあるバックグランド磁場勾配、及びMRI装置内の他の静磁場勾配の定量値をも決定する手段と;
    バックグランド磁場勾配及び他の静磁場勾配に対応する方向に沿って同じマグニチュードの小さい反対の磁場勾配を加える手段とを含む装置。
  2. 非対称MRIパルスシーケンスがNMR信号を得るために使用され、そしてバックグランド静磁場勾配に対する定量値が、下記の2つの関係式によって決定され、
    Figure 2009000538
    ここに、ρ(Z)が、選択したラインのスピン分布であり;Gz(Gzバーを意味する)が、測定又は読み出し方向(この場合には“Z”軸線と呼ばれる)に沿った静バックグランド磁場勾配であり;φ0(Z)は、読み出し勾配の整調誤りによる誤差からの位相コントリビューションであり、そして、信号位相が下記の2つの関係式により読み出し方向に沿った空間位置の関数であり;
    Figure 2009000538
    そしてバックグランド磁場勾配が、下記の関係により信号位相から決定され;
    Figure 2009000538
    ここに、ФおよびФ(φドットおよびφドットを意味する)は、選択した読み出し方向に沿った信号位相の空間的導関数であり;そして△Tは、スピン−エコー位置からのASE信号の時間シフト量である、請求項1に記載の装置。
JP2008186366A 2000-02-24 2008-07-17 Nmrイメージング操作中に誘導される渦電流の測定及び補償のための装置 Pending JP2009000538A (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US09/512,323 US6448773B1 (en) 2000-02-24 2000-02-24 Method and system for measuring and compensating for eddy currents induced during NMR imaging operations

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2001049792A Division JP4201993B2 (ja) 2000-02-24 2001-02-26 Nmrイメージング操作中に誘導される渦電流の測定および補償のための装置及びシステム

Publications (1)

Publication Number Publication Date
JP2009000538A true JP2009000538A (ja) 2009-01-08

Family

ID=24038616

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2001049792A Expired - Lifetime JP4201993B2 (ja) 2000-02-24 2001-02-26 Nmrイメージング操作中に誘導される渦電流の測定および補償のための装置及びシステム
JP2008186366A Pending JP2009000538A (ja) 2000-02-24 2008-07-17 Nmrイメージング操作中に誘導される渦電流の測定及び補償のための装置

Family Applications Before (1)

Application Number Title Priority Date Filing Date
JP2001049792A Expired - Lifetime JP4201993B2 (ja) 2000-02-24 2001-02-26 Nmrイメージング操作中に誘導される渦電流の測定および補償のための装置及びシステム

Country Status (2)

Country Link
US (1) US6448773B1 (ja)
JP (2) JP4201993B2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10588686B2 (en) 2012-06-26 2020-03-17 Covidien Lp Surgical instruments with structures to provide access for cleaning

Families Citing this family (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10104365C1 (de) * 2001-02-01 2002-08-22 Bruker Biospin Gmbh Supraleitendes Magnetsystem und magnetisches Resonanzspektrometer sowie Verfahre zu dessen Betrieb
CN100354643C (zh) * 2002-03-21 2007-12-12 深圳安科高技术股份有限公司 一种梯度预加强电路
JP3907182B2 (ja) * 2002-05-07 2007-04-18 株式会社東芝 磁気共鳴映像装置
ITSV20020030A1 (it) * 2002-07-01 2004-01-02 Esaote Spa Metodo di compensazione di correnti parassite causate da gradienti nelle macchine per il rilevametno d'immagini in risonanza magnetica nucle
JP3930439B2 (ja) 2003-02-06 2007-06-13 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 渦電流補正方法および磁気共鳴撮影装置
US6850066B2 (en) 2003-05-15 2005-02-01 Ge Medical Systems Global Technology Company Llc Systems and methods for gradient compensation in magnetic resonance imaging
US7112964B2 (en) * 2004-08-02 2006-09-26 General Electric Company Eddy current measurement and correction in magnetic resonance imaging systems with a static phantom
US7375522B2 (en) * 2006-08-28 2008-05-20 Wisconsin Alumni Research Foundation Method for aligning multiple MR images acquired with alternating readout gradient
JP5366484B2 (ja) * 2007-09-28 2013-12-11 株式会社東芝 磁気共鳴イメージング装置およびこの磁気共鳴イメージング装置における脂肪抑制効果の分析方法
US9297876B2 (en) * 2009-06-10 2016-03-29 Hitachi Medical Corporation Magnetic resonance imaging apparatus and eddy current compensation method
JP2012183233A (ja) * 2011-03-07 2012-09-27 Toshiba Corp 磁気共鳴イメージング装置
US9322892B2 (en) 2011-12-20 2016-04-26 General Electric Company System for magnetic field distortion compensation and method of making same
US9279871B2 (en) 2011-12-20 2016-03-08 General Electric Company System and apparatus for compensating for magnetic field distortion in an MRI system
US9274188B2 (en) 2012-11-30 2016-03-01 General Electric Company System and apparatus for compensating for magnetic field distortion in an MRI system
DE102013202217B4 (de) 2013-02-12 2015-05-28 Siemens Aktiengesellschaft MR-Anlage mit gepulsten Ausgleichsmagnetfeldgradienten
US9268003B2 (en) * 2013-02-13 2016-02-23 Sunnybrook Health Sciences Centre System and method for measuring induced radio frequency current using phase contrast magnetic resonance imaging
US9709653B2 (en) 2013-02-19 2017-07-18 Toshiba Medical Systems Corporation Mapping eddy current fields in MRI system
US9784715B2 (en) * 2013-04-19 2017-10-10 Zetec, Inc. Eddy current inspection probe based on magnetoresistive sensors
DE102015214925B4 (de) * 2014-09-11 2019-06-06 Siemens Healthcare Gmbh Verfahren zum Betrieb einer Magnetresonanzeinrichtung und Magnetresonanzeinrichtung
CA3115673A1 (en) * 2014-11-11 2016-05-19 Hyperfine Research, Inc. Pulse sequences for low field magnetic resonance
WO2016182407A1 (ko) 2015-05-14 2016-11-17 아탈라에르긴 자기 공명 영상 스캐너
WO2016195281A1 (ko) 2015-05-21 2016-12-08 아탈라에르긴 경사자장을 발생시키기 위해 복수의 코일을 이용하는 경사자장 발생 모듈
US10571537B2 (en) * 2015-05-21 2020-02-25 Bilkent University Multi-purpose gradient array for magnetic resonance imaging
US10641858B2 (en) 2017-04-06 2020-05-05 Bilkent University Spatiotemporal magnetic field monitoring with hall effect sensors during the MRI scan
TW202012951A (zh) 2018-07-31 2020-04-01 美商超精細研究股份有限公司 低場漫射加權成像
CN111856360A (zh) * 2019-04-24 2020-10-30 通用电气精准医疗有限责任公司 用于获取磁共振成像数据的方法和磁共振成像系统
US11510588B2 (en) 2019-11-27 2022-11-29 Hyperfine Operations, Inc. Techniques for noise suppression in an environment of a magnetic resonance imaging system

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5099208A (en) * 1989-10-05 1992-03-24 Vanderbilt University Method for magnetic resonance imaging and related apparatus
JPH0584227A (ja) * 1991-09-30 1993-04-06 Toshiba Corp 磁気共鳴診断装置
JPH0686764A (ja) * 1992-09-09 1994-03-29 Yokogawa Medical Syst Ltd Mri装置
JPH06133945A (ja) * 1992-10-30 1994-05-17 Toshiba Corp Mri及びmrsにおける静磁場均一性の調整方法
JPH10272120A (ja) * 1996-12-30 1998-10-13 General Electric Co <Ge> 磁気共鳴システムの磁場を補償する方法
JPH1147110A (ja) * 1997-07-30 1999-02-23 Shimadzu Corp 磁気共鳴イメージング装置

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4698591A (en) 1986-01-03 1987-10-06 General Electric Company Method for magnetic field gradient eddy current compensation
US4924186A (en) * 1986-09-29 1990-05-08 Kabushiki Kaisha Toshiba Magnetic resonance imaging system with auxiliary compensation coil
US4950994A (en) 1988-03-07 1990-08-21 General Electric Company Gradient and polarizing field compensation
US4885542A (en) 1988-04-14 1989-12-05 The Regents Of The University Of California MRI compensated for spurious NMR frequency/phase shifts caused by spurious changes in magnetic fields during NMR data measurement processes
US4970457A (en) 1989-04-05 1990-11-13 The Regents Of The University Of California MRI compensated for spurious rapid variations in static magnetic field during a single MRI sequence
US4978919A (en) * 1989-04-27 1990-12-18 Picker International, Inc. Measurement and calibration of eddy currents for magnetic resonance imagers
US4965521A (en) 1989-08-11 1990-10-23 Spectroscopy Imaging Systems Method and apparatus for compensating eddy current effects in a magnetic resonance device having pulsed magnetic field gradients
US5151656A (en) 1990-12-11 1992-09-29 General Electric Company Correction of nmr data acquired by an echo-planar technique
DE69211806T2 (de) 1991-10-25 1996-10-31 Univ Queensland Korrektur von Signalverzerrungen in einem magnetischen Kernresonanzapparat
DE4203036C2 (de) 1992-02-04 1994-04-14 Bruker Medizintech Verfahren zur Erzeugung eines seitenbandfreien NMR-Spektrums
US5493224A (en) 1992-03-03 1996-02-20 Hitachi, Ltd. Ultra high-speed magnetic resonance imaging method and apparatus
DE69326202T2 (de) 1992-06-29 2000-04-13 Koninkl Philips Electronics Nv Verfahren und Gerät zur Kernresonanzabbildung
DE69329873T2 (de) 1992-10-26 2001-09-13 Koninkl Philips Electronics Nv Kompensierung von Wirbelströmen in der Bilderzeugung mittels magnetischer Resonanz
DE69418404T2 (de) 1993-09-16 1999-11-11 Koninkl Philips Electronics Nv Korrektur der Polarität des Auslesegradienten in Bilderzeugung durch EPI und GRASE magnetische Resonanz
US5672969A (en) 1996-04-25 1997-09-30 General Electric Company Reduction of Nyquist ghost artifacts in oblique echo planar imaging
US5864233A (en) 1997-02-22 1999-01-26 General Electric Company Method to reduce eddy current effects in diffusion-weighted echo planar imaging
US6127826A (en) * 1999-01-27 2000-10-03 Picker International, Inc. EPI image based long term eddy current pre-emphasis calibration

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5099208A (en) * 1989-10-05 1992-03-24 Vanderbilt University Method for magnetic resonance imaging and related apparatus
JPH0584227A (ja) * 1991-09-30 1993-04-06 Toshiba Corp 磁気共鳴診断装置
JPH0686764A (ja) * 1992-09-09 1994-03-29 Yokogawa Medical Syst Ltd Mri装置
JPH06133945A (ja) * 1992-10-30 1994-05-17 Toshiba Corp Mri及びmrsにおける静磁場均一性の調整方法
JPH10272120A (ja) * 1996-12-30 1998-10-13 General Electric Co <Ge> 磁気共鳴システムの磁場を補償する方法
JPH1147110A (ja) * 1997-07-30 1999-02-23 Shimadzu Corp 磁気共鳴イメージング装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10588686B2 (en) 2012-06-26 2020-03-17 Covidien Lp Surgical instruments with structures to provide access for cleaning

Also Published As

Publication number Publication date
JP4201993B2 (ja) 2008-12-24
JP2001258865A (ja) 2001-09-25
US6448773B1 (en) 2002-09-10

Similar Documents

Publication Publication Date Title
JP4201993B2 (ja) Nmrイメージング操作中に誘導される渦電流の測定および補償のための装置及びシステム
KR100481740B1 (ko) 와전류들에의해유도된공간적및시간적으로변화하는자계들을측정하고보상하는방법
Boesch et al. Temporal and spatial analysis of fields generated by eddy currents in superconducting magnets: optimization of corrections and quantitative characterization of magnet/gradient systems
Barmet et al. A transmit/receive system for magnetic field monitoring of in vivo MRI
US6891373B2 (en) Method to determine the ADC coefficients in diffusion-weighted magnetic resonance imaging given use of steady-state sequences
US5450010A (en) Magnetic resonance imaging method and apparatus employing eddy current compensation by modification of gradient size
JP3402916B2 (ja) 核スピントモグラフ装置のマグネット装置のシム調整方法及び方法を実施する装置
US20040113615A1 (en) Characterization of the effect of spatial gradient field distortions in diffusion-weighted imaging
US6507190B1 (en) Method and apparatus for compensating polarizing fields in magnetic resonance imaging
US7689262B2 (en) Method and magnetic resonance tomography apparatus for correcting changes in the basic magnetic field
EP0395248B1 (en) Magnetic field and eddy current measuring method
JP2716889B2 (ja) 磁石の高速補正法
EP0100183B1 (en) Nuclear magnetic resonance method and apparatus
EP3060116B1 (en) Mr imaging with temperature mapping
US8717022B2 (en) Magnetic field gradient monitor apparatus and method
US5451877A (en) Method for the compensation of eddy currents caused by gradients in a nuclear magnetic resonance apparatus
Bartusek et al. Determination of pre-emphasis constants for eddy current reduction
US6127826A (en) EPI image based long term eddy current pre-emphasis calibration
US5856744A (en) Method and apparatus for B0 eddy current compensation by B1 phase modulation
Vashaee et al. B1 mapping with a pure phase encode approach: quantitative density profiling
JP2002529175A (ja) Mriの勾配電磁界歪みの自動測定
EP0439119A2 (en) Magnetic resonance imaging system
US20010054898A1 (en) Magnetic resonance imaging compensated for very rapid variations in static magnetic field
US6229309B1 (en) MR method
JPH0418856B2 (ja)

Legal Events

Date Code Title Description
RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20100811

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110426

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110714

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110830

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20111130

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20111205

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20120418

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20120522