JPH01141656A - Ｎｍｒイメージング装置の磁場ドリフト補正方法 - Google Patents
Ｎｍｒイメージング装置の磁場ドリフト補正方法Info
- Publication number
- JPH01141656A JPH01141656A JP62301955A JP30195587A JPH01141656A JP H01141656 A JPH01141656 A JP H01141656A JP 62301955 A JP62301955 A JP 62301955A JP 30195587 A JP30195587 A JP 30195587A JP H01141656 A JPH01141656 A JP H01141656A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- phase
- drift
- main magnetic
- signal
- 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
Links
- 238000000034 method Methods 0.000 title claims description 21
- 238000003384 imaging method Methods 0.000 title claims description 18
- 238000012937 correction Methods 0.000 claims abstract description 10
- 238000001514 detection method Methods 0.000 claims abstract description 5
- 230000005415 magnetization Effects 0.000 claims description 11
- 239000013598 vector Substances 0.000 claims description 11
- 238000005070 sampling Methods 0.000 abstract description 5
- 230000006866 deterioration Effects 0.000 abstract description 4
- 238000005259 measurement Methods 0.000 description 14
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 2
- 241001203771 Eudonia echo Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000000498 cooling water Substances 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000002542 deteriorative effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
（産業上の利用分野）
本発明は、１ビュー期間中に、Ｒ影用のＮＭＲ信号の収
集と共に、主磁場ドリフト量測定の為の信号の収集を行
うＮＭＲイメージング装置の磁場ドリフト補正方法に関
し、更に詳しくは前記ドリフトｆｆｉ測定の為の信号の
処理方法の改善に関する。
集と共に、主磁場ドリフト量測定の為の信号の収集を行
うＮＭＲイメージング装置の磁場ドリフト補正方法に関
し、更に詳しくは前記ドリフトｆｆｉ測定の為の信号の
処理方法の改善に関する。
（従来の技術）
第４図は従来例のＮＭＲイメージング装置の概略ブロッ
ク図である。第４図において、１は被検体（図示せず）
を載置し、所定の位置に設置するテーブル、２は−様な
静磁場を発生する主磁場マグネットや、被検体にＲＦ
（Ｒａｄｉｏ−ｆｒｅｑｕｅｎｃｙ）パルスを印加し、
ＮＭＲ現象によって被検体に生じた信号を受信するＲＦ
コイルや、ＲＦ信号に３次元の位置情報を与える目的で
各方向に線形勾配磁場ａｘ、 ＧＶ、 Ｇ’ｌを印加す
る勾配磁場コイル等を配置したマグネットアセンブリで
ある。３は前記主磁場マグネットに電流を流す主磁場電
源、４は前記勾配磁場コイルにパルス電流を与える勾配
磁場駆動回路、５は所定の周波数のＲＦ基準信号を発振
するＲＦ発振回路、６はＲＦ基準信号をゲート・変調回
路７によって所定の幅（時間）の変調波パルス波形とし
た信号を、所定のゲインで増幅後ＲＦコイルに与えるＲ
Ｆパワーアンプ、８はＲＦコイルによって受信されたＮ
ＭＲ信号を増幅するプリアンプ、９は増幅されたＮＭＲ
信号を位相検波して復調する位相検波器、１０は検波さ
れた信号をＡ／Ｄ変換するＡ／Ｄ変換器、１１はＡ／Ｄ
変換された信号をスキャンシーケンスに応じた所定の画
像再構成処理により画像化し、デイスプレィ１２に表示
するコンビコータ、１３はＲＦパルスの送受信や勾配磁
場パルスの印加タイミング等のスキャンシーケンスをＩ
制御するシーケンスメモリ、１４はパルスシーケンスの
設定や画像処理を行う操作コンソールである。
ク図である。第４図において、１は被検体（図示せず）
を載置し、所定の位置に設置するテーブル、２は−様な
静磁場を発生する主磁場マグネットや、被検体にＲＦ
（Ｒａｄｉｏ−ｆｒｅｑｕｅｎｃｙ）パルスを印加し、
ＮＭＲ現象によって被検体に生じた信号を受信するＲＦ
コイルや、ＲＦ信号に３次元の位置情報を与える目的で
各方向に線形勾配磁場ａｘ、 ＧＶ、 Ｇ’ｌを印加す
る勾配磁場コイル等を配置したマグネットアセンブリで
ある。３は前記主磁場マグネットに電流を流す主磁場電
源、４は前記勾配磁場コイルにパルス電流を与える勾配
磁場駆動回路、５は所定の周波数のＲＦ基準信号を発振
するＲＦ発振回路、６はＲＦ基準信号をゲート・変調回
路７によって所定の幅（時間）の変調波パルス波形とし
た信号を、所定のゲインで増幅後ＲＦコイルに与えるＲ
Ｆパワーアンプ、８はＲＦコイルによって受信されたＮ
ＭＲ信号を増幅するプリアンプ、９は増幅されたＮＭＲ
信号を位相検波して復調する位相検波器、１０は検波さ
れた信号をＡ／Ｄ変換するＡ／Ｄ変換器、１１はＡ／Ｄ
変換された信号をスキャンシーケンスに応じた所定の画
像再構成処理により画像化し、デイスプレィ１２に表示
するコンビコータ、１３はＲＦパルスの送受信や勾配磁
場パルスの印加タイミング等のスキャンシーケンスをＩ
制御するシーケンスメモリ、１４はパルスシーケンスの
設定や画像処理を行う操作コンソールである。
このような構成のＮＭＲイメージング装置においては、
主磁場マグネットがコイルに低電流を流すことにより主
磁場を発生する常伝導型のコイルで構成されている場合
、コイル冷却水温度の変化、室温変化、′ＩＩｉ源の電
流検出部のシャント抵抗の温度変化等によるコイルＮ流
の変動により、主磁場が時間的に変動する磁場ドリフト
が生じる。又、主磁場マグネットが永久磁石の場合は、
室温変化に伴う磁石の磁性変化などにより磁場ドリフト
が生じる。この様な磁場ドリフトは、得られる画像にボ
ケ、ゴーストなどの画質劣化を引き起こすなどの悪Ｗｅ
Ｗがあるため、一般に、主磁場強度を測定し、その誤差
分くドリフト徂）に応じて、コイルに流す１！流値をフ
ィードバック調節したり、ＲＥ周波数を磁場に会わせて
シフトするいわゆるＮＭＲロック法が用いられる。しか
し、このＮＭＲロック法においては、ドリフトを抑えた
いスキャン中には、前記主磁場測定が困難なために、充
分な磁場ドリフトの低減ができなかった。
主磁場マグネットがコイルに低電流を流すことにより主
磁場を発生する常伝導型のコイルで構成されている場合
、コイル冷却水温度の変化、室温変化、′ＩＩｉ源の電
流検出部のシャント抵抗の温度変化等によるコイルＮ流
の変動により、主磁場が時間的に変動する磁場ドリフト
が生じる。又、主磁場マグネットが永久磁石の場合は、
室温変化に伴う磁石の磁性変化などにより磁場ドリフト
が生じる。この様な磁場ドリフトは、得られる画像にボ
ケ、ゴーストなどの画質劣化を引き起こすなどの悪Ｗｅ
Ｗがあるため、一般に、主磁場強度を測定し、その誤差
分くドリフト徂）に応じて、コイルに流す１！流値をフ
ィードバック調節したり、ＲＥ周波数を磁場に会わせて
シフトするいわゆるＮＭＲロック法が用いられる。しか
し、このＮＭＲロック法においては、ドリフトを抑えた
いスキャン中には、前記主磁場測定が困難なために、充
分な磁場ドリフトの低減ができなかった。
そこで、スキャン中の各ビュー毎に主磁場を測定して、
磁場ドリフトを低減する方法が提案されている（特開昭
θ０−２３６６３７号）。第５図は、この従来例のスキ
ャンシーケンスを表わすタイムチャートである。第５図
のように、始めに９０”のＲＦパルスと共にスライス勾
配Ｇｚを印加し選択励起を行い磁化ベクトルを回転座標
系上でＸ方向に倒し、投影のための勾配磁ＪＩＧＸ、
Ｇｙを印加し、イメージング（Ｉ彩）用のエコー信号１
００を得る。その後、主磁場強度を測定するために１８
０゛パルスによりて磁化ベクトルを反転させた後、再び
勾配磁場ＧＸ、 ＧＶを印加し、２度目のエコー信号が
最大になった時点で勾配磁場を０にして断層面全体から
のエコー信号、即ち磁場ドリフト邊胴定の為のエコー信
号２００の収集を行う。第６図はこのエコー信号２００
を周波数解析したフーリエスペクトラムである。このフ
ーリエスペクトラムの中心周波数ω１を求め、磁場強度
ＢＯとラーモア周波数ωの関係式、即ちω−γＢＯ（：
γ磁気回転比）より主磁場強度を求め、これに応じて画
像再構成データに補正を加えることにより、スキャン中
の磁場ドリフトの［９を低減することができる。
磁場ドリフトを低減する方法が提案されている（特開昭
θ０−２３６６３７号）。第５図は、この従来例のスキ
ャンシーケンスを表わすタイムチャートである。第５図
のように、始めに９０”のＲＦパルスと共にスライス勾
配Ｇｚを印加し選択励起を行い磁化ベクトルを回転座標
系上でＸ方向に倒し、投影のための勾配磁ＪＩＧＸ、
Ｇｙを印加し、イメージング（Ｉ彩）用のエコー信号１
００を得る。その後、主磁場強度を測定するために１８
０゛パルスによりて磁化ベクトルを反転させた後、再び
勾配磁場ＧＸ、 ＧＶを印加し、２度目のエコー信号が
最大になった時点で勾配磁場を０にして断層面全体から
のエコー信号、即ち磁場ドリフト邊胴定の為のエコー信
号２００の収集を行う。第６図はこのエコー信号２００
を周波数解析したフーリエスペクトラムである。このフ
ーリエスペクトラムの中心周波数ω１を求め、磁場強度
ＢＯとラーモア周波数ωの関係式、即ちω−γＢＯ（：
γ磁気回転比）より主磁場強度を求め、これに応じて画
像再構成データに補正を加えることにより、スキャン中
の磁場ドリフトの［９を低減することができる。
（発明が解決しようとする問題点）
しかし、上記の方法においては、磁場測定用の信号の収
集時間を長くしなければ、磁場測定の精度を上げること
ができない。例えば主磁場強度が０．３７（テスラ）の
場合、１ ｐＤＩ磁場が変動すると水素原子核のラーモ
ア周波数は１２．８Ｈｚ変動する。従って、０．ｌｐｐ
ｍの精度で磁場変動を測定するためには、１．２８Ｈｚ
の周波数を周波数解析して検出する必要があり、前記磁
場測定用の信号の収集期間を１／１．２８＝０．８秒と
しなければならない。この様な時間の磁場測定用シーケ
ンスを１ビュー朋問中に加えると、総スキャン時間が長
くなり、種々の弊害が生じる。従って、実際には前記サ
ンプリング期間を充分に長くとり、磁場変動を細かく検
出することは不可能であり、前記スキャン中の磁場ドリ
フトによる画像の劣化を充分に防ぐことができないとい
う問題があった。
集時間を長くしなければ、磁場測定の精度を上げること
ができない。例えば主磁場強度が０．３７（テスラ）の
場合、１ ｐＤＩ磁場が変動すると水素原子核のラーモ
ア周波数は１２．８Ｈｚ変動する。従って、０．ｌｐｐ
ｍの精度で磁場変動を測定するためには、１．２８Ｈｚ
の周波数を周波数解析して検出する必要があり、前記磁
場測定用の信号の収集期間を１／１．２８＝０．８秒と
しなければならない。この様な時間の磁場測定用シーケ
ンスを１ビュー朋問中に加えると、総スキャン時間が長
くなり、種々の弊害が生じる。従って、実際には前記サ
ンプリング期間を充分に長くとり、磁場変動を細かく検
出することは不可能であり、前記スキャン中の磁場ドリ
フトによる画像の劣化を充分に防ぐことができないとい
う問題があった。
従って本願発明の目的は、スキャン時間を長くすること
なく、精度の高い主磁場測定をスキャン中に行い、磁場
ドリフトによる画質の劣化を防ぎ、鮮明な画像を得るこ
とができるＮＭＲイメージング装置の磁場ドリフト補正
方法を実現することにある。
なく、精度の高い主磁場測定をスキャン中に行い、磁場
ドリフトによる画質の劣化を防ぎ、鮮明な画像を得るこ
とができるＮＭＲイメージング装置の磁場ドリフト補正
方法を実現することにある。
（問題点を解決するための手段）
上記の目的を達成するために、本発明のＮＭＲイメージ
ング装置の磁場ドリフト補正方法は以下のような構成を
している。
ング装置の磁場ドリフト補正方法は以下のような構成を
している。
即ち、１ビュー期間中に磁化ベクトルの位相をばらけさ
せたのと同量で極性が反対の勾配磁場を印加して主磁場
ドリフト１１１？１定の為の信号の収集を行ない、この
収集した信号を直交検波した後、■（同相）成分とＱ（
直交）成分の偏角から各ビュー毎の位相を求め、前記位
相の変動量に応じて前記主磁場ドリフトの補正を行うこ
とを特徴とする。
せたのと同量で極性が反対の勾配磁場を印加して主磁場
ドリフト１１１？１定の為の信号の収集を行ない、この
収集した信号を直交検波した後、■（同相）成分とＱ（
直交）成分の偏角から各ビュー毎の位相を求め、前記位
相の変動量に応じて前記主磁場ドリフトの補正を行うこ
とを特徴とする。
（実施例）
以下、図面を参照して本発明について詳細に説明する。
第１図は、本発明の一実施例の２ＤＦＴ （２次元フー
リエ変換）法によるパルスシーケンス（１ビュー分）を
表わすタイムチャートである。
リエ変換）法によるパルスシーケンス（１ビュー分）を
表わすタイムチャートである。
第１図のように、通常の２ＤＦＴ法と同様に、始めに９
０°のＲＦパルスと共にスライス勾配Ｇ５１ｉｃｅを印
加し、断層面の選択励起を行う。次に、ＲＦ倍信号２次
元の位と情報を周波数と位相情報として与える目的で、
読出し方向に勾配磁場Ｇｒｅａｄを印加し、同様に位相
エンコード方向に勾配磁場Ｇｏｈａｓｅ■を印加する。
０°のＲＦパルスと共にスライス勾配Ｇ５１ｉｃｅを印
加し、断層面の選択励起を行う。次に、ＲＦ倍信号２次
元の位と情報を周波数と位相情報として与える目的で、
読出し方向に勾配磁場Ｇｒｅａｄを印加し、同様に位相
エンコード方向に勾配磁場Ｇｏｈａｓｅ■を印加する。
その後、−度ばらけた磁化ベクトルを再び揃えるための
１８０゛パルスを印加し、読出し勾配Ｇｒｅａｄを印加
しながらイメージング用のスピンエコー信号３００を得
る。このスピンエコー信号は、それまでに期間ＴＥで印
加された読出し方向の勾配磁場Ｇｒｅａｄの総和がＯの
とき最大となり、その後次第に減衰する。ここまでは通
常のスキャンと同様であるが、この後、主磁場測定用の
Ｆ ■Ｑ （Ｆｒｅｅ Ｔｎｄｕｃｔｉｏｎ Ｄｅｃａ
ｙ）信号を得るために、磁化ベクトルの位相をばらけさ
せたのと同量で極性が反対の勾配Ｗｉ場を印加する。
１８０゛パルスを印加し、読出し勾配Ｇｒｅａｄを印加
しながらイメージング用のスピンエコー信号３００を得
る。このスピンエコー信号は、それまでに期間ＴＥで印
加された読出し方向の勾配磁場Ｇｒｅａｄの総和がＯの
とき最大となり、その後次第に減衰する。ここまでは通
常のスキャンと同様であるが、この後、主磁場測定用の
Ｆ ■Ｑ （Ｆｒｅｅ Ｔｎｄｕｃｔｉｏｎ Ｄｅｃａ
ｙ）信号を得るために、磁化ベクトルの位相をばらけさ
せたのと同量で極性が反対の勾配Ｗｉ場を印加する。
即ち、位相情報を与えるために印加した勾配磁場Ｇｐｈ
ａｓｅ■と面積（強度Ｘ時ＩＩ）が等しく符号が反対の
勾配磁場ＧｐｈａＳｅ■°を印加する。又、読出し方向
には、前記スピンエコー信号が最大の時点（ＴＥ ）以
降に印加した読出し勾配Ｇｒｅａｄ■と面積が等しく符
号が反対の勾配磁場Ｇｒｅａｄ■゛を印加する。ここで
、この期間の緩和ＴＩ 、Ｔ２を無視すれば、磁化ベク
トルに与えられる位相回転量θは、次式のようになる。
ａｓｅ■と面積（強度Ｘ時ＩＩ）が等しく符号が反対の
勾配磁場ＧｐｈａＳｅ■°を印加する。又、読出し方向
には、前記スピンエコー信号が最大の時点（ＴＥ ）以
降に印加した読出し勾配Ｇｒｅａｄ■と面積が等しく符
号が反対の勾配磁場Ｇｒｅａｄ■゛を印加する。ここで
、この期間の緩和ＴＩ 、Ｔ２を無視すれば、磁化ベク
トルに与えられる位相回転量θは、次式のようになる。
即ち、
θ−ｆ１″Ｑ−ｔｄｔ・・・（１）
ｉ
となり、１ビュー間で与えられた勾配磁場による磁化ベ
クトルの位相回転量は相殺され、Ｇｐｈａｓｅ■°及び
Ｇｒｅａｄ■゛を印加し終わった時点ｔａで断層面全体
の磁化ベクトルの位相が一致し、ＦＩＤ信号４００を誘
起する。但し、（１）式において、Ｇは読出し又は位相
エンコード方向の勾配磁場の強度であり、ｔ（口〜ｔ２
）は勾配磁場を印加した時間である。以下、このように
して１ビュー期間中に得られたＦＩＤ信号を主磁場測定
に用いて、主磁場ドリフト量を求める信号処理方法を説
明する。第２図は前記ＦＩＤ信号を直交検波したＩ （
Ｉｎｐｈａｓｅ）とＱ （Ｑｕａｄｒａｔｕｒｅ）成分
を表わす図である。第２図において、７’ｓは収集した
ＦＩＤ信号のデータ収集時間である。第２図のように、
始めに前記ＦＩＤ信号は、通常のイメージング用データ
と同様に、中心周波数ｆＯで直交検波され、■とＱの複
素数データ即ち、同相と直交成分に分離され、所定の時
＠間隔でＮ３点サンプリングされる。ここで、第３図は
サンプリングされた１１Ｑデータから磁場ドリフト補正
１を求めるフロー図である。まず、Ｋビュー目のＦＩＤ
信号の各サンプリング点（時刻） ｉでのデータ（ｉ
、Ｑｉの偏角から各サンプル点での位相θ１（ｉ−ｉ〜
ＮＳ）を次式のように弾出する。
クトルの位相回転量は相殺され、Ｇｐｈａｓｅ■°及び
Ｇｒｅａｄ■゛を印加し終わった時点ｔａで断層面全体
の磁化ベクトルの位相が一致し、ＦＩＤ信号４００を誘
起する。但し、（１）式において、Ｇは読出し又は位相
エンコード方向の勾配磁場の強度であり、ｔ（口〜ｔ２
）は勾配磁場を印加した時間である。以下、このように
して１ビュー期間中に得られたＦＩＤ信号を主磁場測定
に用いて、主磁場ドリフト量を求める信号処理方法を説
明する。第２図は前記ＦＩＤ信号を直交検波したＩ （
Ｉｎｐｈａｓｅ）とＱ （Ｑｕａｄｒａｔｕｒｅ）成分
を表わす図である。第２図において、７’ｓは収集した
ＦＩＤ信号のデータ収集時間である。第２図のように、
始めに前記ＦＩＤ信号は、通常のイメージング用データ
と同様に、中心周波数ｆＯで直交検波され、■とＱの複
素数データ即ち、同相と直交成分に分離され、所定の時
＠間隔でＮ３点サンプリングされる。ここで、第３図は
サンプリングされた１１Ｑデータから磁場ドリフト補正
１を求めるフロー図である。まず、Ｋビュー目のＦＩＤ
信号の各サンプリング点（時刻） ｉでのデータ（ｉ
、Ｑｉの偏角から各サンプル点での位相θ１（ｉ−ｉ〜
ＮＳ）を次式のように弾出する。
θｉ −ｔａｎ ”（Ｑｉ ／Ｉｉ ）・・・（２）次
に、（２）式で求められた各サンプリング点での位相の
平均を次式のように求める。
に、（２）式で求められた各サンプリング点での位相の
平均を次式のように求める。
θに一Σθｉ ／Ｎｓ ・・・（３）ｉ鼻１
ここで、このθにはにビュー目の位相を表わしている。
次に（３）式で求められたにビュー目の位相θにと１つ
前のに一１ビュー目の位相θに−１の位相差 Δ８に一θに−１−θｋ ・−（４）を求める。こ
の位相差Δθに即ち、ビュー間の位相の変動ｍは、前の
ビューの時点からの磁場の変動により生じたラーモア周
波数のずれを表わしている。従って、位相差Δθｋから
前のビューとのラーモア周波数差Δｆ、を次式から求め
る。
前のに一１ビュー目の位相θに−１の位相差 Δ８に一θに−１−θｋ ・−（４）を求める。こ
の位相差Δθに即ち、ビュー間の位相の変動ｍは、前の
ビューの時点からの磁場の変動により生じたラーモア周
波数のずれを表わしている。従って、位相差Δθｋから
前のビューとのラーモア周波数差Δｆ、を次式から求め
る。
Δｆ−Δθｋ ／ （２π−Ｔｖ ）−（５）但し、Ｔ
Ｖはに一１ビューとにビューとのＦＲ間間隔である。最
後に、この周波数のずれΔｆに応じてコイルに流す補正
電流値を計算してフィードバック調整し、主ｍｓの変動
を補正して終了する。このように本実施例の方法におい
ては、主１！１１８１定用のＦＩＤ信号を周波数解析せ
ずに、直接位相を求め、各ビュー毎の位相の変動ｍから
、主磁場ドリフト吊を求めているため、前記ドリフト１
の測定精度は、ＦＩＤ信号のサンプリングデータ長ＴＳ
に依存しない。従って、必要な主磁場ドリフト伍の測定
精度を得るためにＦＩＤ信号の収集時間を長くする必要
がない。従って、スキャン時間を長くすることなく、１
ビュー期間中に精度の高い磁場ドリフト測定を行うこと
ができる。
Ｖはに一１ビューとにビューとのＦＲ間間隔である。最
後に、この周波数のずれΔｆに応じてコイルに流す補正
電流値を計算してフィードバック調整し、主ｍｓの変動
を補正して終了する。このように本実施例の方法におい
ては、主１！１１８１定用のＦＩＤ信号を周波数解析せ
ずに、直接位相を求め、各ビュー毎の位相の変動ｍから
、主磁場ドリフト吊を求めているため、前記ドリフト１
の測定精度は、ＦＩＤ信号のサンプリングデータ長ＴＳ
に依存しない。従って、必要な主磁場ドリフト伍の測定
精度を得るためにＦＩＤ信号の収集時間を長くする必要
がない。従って、スキャン時間を長くすることなく、１
ビュー期間中に精度の高い磁場ドリフト測定を行うこと
ができる。
尚、本発明は、上記実施例に限定するものではなく、特
許請求の範囲内で種々の変形が可能である。例えば、前
記磁場ドリフトの補正は、主磁場マグネットに流す電流
値を調節する以外にも、ＲＦパルスの周波数を調節して
も良い。又、撮影用のＮＭＲ信号に次のような位相処理
を施しても良い。即ち、第にビュー目で得られたＮＭＲ
信号を・ 読出し方向に１次元フーリエ変換したデー
タの各々に、ｅｉΔθｋを掛けて、磁場ドリフトによっ
て生じる各ビュー間での位相エラーを補償しても良い。
許請求の範囲内で種々の変形が可能である。例えば、前
記磁場ドリフトの補正は、主磁場マグネットに流す電流
値を調節する以外にも、ＲＦパルスの周波数を調節して
も良い。又、撮影用のＮＭＲ信号に次のような位相処理
を施しても良い。即ち、第にビュー目で得られたＮＭＲ
信号を・ 読出し方向に１次元フーリエ変換したデー
タの各々に、ｅｉΔθｋを掛けて、磁場ドリフトによっ
て生じる各ビュー間での位相エラーを補償しても良い。
更に、１１ｉｌドリフト量測定のための信号はＦＩＤ信
号でなく、エコー信号でも良い。更に、上記の方法を適
宜組み合わせて用いても良い。更に、２次元フーリエ変
換法以外にも、マルチエコー法や３次元フーリエ変換法
やプロジェクション法に用いても良い。
号でなく、エコー信号でも良い。更に、上記の方法を適
宜組み合わせて用いても良い。更に、２次元フーリエ変
換法以外にも、マルチエコー法や３次元フーリエ変換法
やプロジェクション法に用いても良い。
（発明の効果）
以上説明の通り、本発明のＮＭＲイメージング装置の磁
場ドリフト補正方法によれば、以下の効果が得られる。
場ドリフト補正方法によれば、以下の効果が得られる。
即ち、撮影用のＮＭＲ信号を検出した後、１ビュー期間
中に磁化ベクトルの位相をばらけさせたのと同量で極性
が反対の勾配磁場を印加して前記ドリフトｌ測定の為の
信号の収集を行ない、この収集した信号を直交検波した
後、■（同相）成分とＱ（直交）成分の偏角から各ビュ
ー毎の位相を求め、前記位相の変動遣に応じて前記主磁
場ドリフトの補正を行うため、前記ドリフト量測定の為
の信号を周波数解析せずに位相量から前記ドリフト吊が
求められ、主磁場ドリフト伍の測定精度は、前記信号の
データ長（収集時間）に依存しない。従って、必要な主
磁場ドリフト量の測定粘度を得るために、この測定のた
めの信号を長く収集する必要がない。このため、スキャ
ン時間を長くすることなく、１ビュー期間中に精度の高
い１１場ドリフト測定を行うことができる。これにより
、スキャン中の磁場ドリフトによる画質の劣化を防ぎ、
鮮明な画像を得ることができる。
中に磁化ベクトルの位相をばらけさせたのと同量で極性
が反対の勾配磁場を印加して前記ドリフトｌ測定の為の
信号の収集を行ない、この収集した信号を直交検波した
後、■（同相）成分とＱ（直交）成分の偏角から各ビュ
ー毎の位相を求め、前記位相の変動遣に応じて前記主磁
場ドリフトの補正を行うため、前記ドリフト量測定の為
の信号を周波数解析せずに位相量から前記ドリフト吊が
求められ、主磁場ドリフト伍の測定精度は、前記信号の
データ長（収集時間）に依存しない。従って、必要な主
磁場ドリフト量の測定粘度を得るために、この測定のた
めの信号を長く収集する必要がない。このため、スキャ
ン時間を長くすることなく、１ビュー期間中に精度の高
い１１場ドリフト測定を行うことができる。これにより
、スキャン中の磁場ドリフトによる画質の劣化を防ぎ、
鮮明な画像を得ることができる。
第１図は、本発明の一実施例の２ＤＦＴ法によるパルス
シーケンスを表わすタイムチャート、第２図は、本発明
の一実施例の直交検波されたＩとＱ成分を表わす図、第
３図は本発明の一実施例の磁場ドリフト補正量を求める
フロー図、第４図は従来例のＮＭＲイメージング装置の
概略ブロック図、第５図は従来例のスキャンシーケンス
を表わすタイムチャート、第６図はエコー信＠２００を
周波数解析したフーリエスペクトラムである。 １・・・テーブル、２・・・マグネットアセンブリ、３
・・・１磁Ｓ電源、４・・・勾配磁場駆動回路、５・・
・ＲＦ発発註回路６・・・ＲＦパワーアンプ、７・・・
ゲート・変調回路、８・・・プリアンプ、９・・・位相
検波器、１ｏ・・・Ａ／Ｄ変換器、１１・・・コンピュ
ータ、１２・・・デイスプレィ、１３・・・シーケンス
メモリ、１４・・・操作コンソール、１００．３００・
・・イメージング用エコー信号２００・・・磁場測定用
エコー信号 ４００・・・磁場測定用ＦＩＤ信号
シーケンスを表わすタイムチャート、第２図は、本発明
の一実施例の直交検波されたＩとＱ成分を表わす図、第
３図は本発明の一実施例の磁場ドリフト補正量を求める
フロー図、第４図は従来例のＮＭＲイメージング装置の
概略ブロック図、第５図は従来例のスキャンシーケンス
を表わすタイムチャート、第６図はエコー信＠２００を
周波数解析したフーリエスペクトラムである。 １・・・テーブル、２・・・マグネットアセンブリ、３
・・・１磁Ｓ電源、４・・・勾配磁場駆動回路、５・・
・ＲＦ発発註回路６・・・ＲＦパワーアンプ、７・・・
ゲート・変調回路、８・・・プリアンプ、９・・・位相
検波器、１ｏ・・・Ａ／Ｄ変換器、１１・・・コンピュ
ータ、１２・・・デイスプレィ、１３・・・シーケンス
メモリ、１４・・・操作コンソール、１００．３００・
・・イメージング用エコー信号２００・・・磁場測定用
エコー信号 ４００・・・磁場測定用ＦＩＤ信号
Claims (1)
- 撮影面内の磁化ベクトルを励起し、所定の勾配磁場を印
加して撮影用のＮＭＲ信号の収集を行つた後、勾配磁場
を０にして、主磁場ドリフト量測定の為の信号の収集を
行うＮＭＲイメージング装置の磁場ドリフト補正方法に
おいて、１ビュー期間中に前記磁化ベクトルの位相をば
らけさせたのと同量で極性が反対の勾配磁場を印加して
前記ドリフト量測定の為の信号の収集を行ない、この収
集した信号を直交検波した後、Ｉ（同相）成分とＱ（直
交）成分の偏角から各ビュー毎の位相を求め、前記位相
の変動量に応じて前記主磁場ドリフトの補正を行うこと
を特徴とするＮＭＲイメージング装置の磁場ドリフト補
正方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP62301955A JPH01141656A (ja) | 1987-11-30 | 1987-11-30 | Ｎｍｒイメージング装置の磁場ドリフト補正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP62301955A JPH01141656A (ja) | 1987-11-30 | 1987-11-30 | Ｎｍｒイメージング装置の磁場ドリフト補正方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH01141656A true JPH01141656A (ja) | 1989-06-02 |
Family
ID=17903127
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP62301955A Pending JPH01141656A (ja) | 1987-11-30 | 1987-11-30 | Ｎｍｒイメージング装置の磁場ドリフト補正方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH01141656A (ja) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1102082A2 (en) * | 1999-11-15 | 2001-05-23 | General Electric Company | Method and apparatus for reducing image artifacts caused by magnet vibration in an MR imaging system |
| EP1102081A2 (en) * | 1999-11-15 | 2001-05-23 | General Electric Company | Method and apparatus to compensate for image artifacts caused by magnet vibration in an MR imaging system |
| JP2001190519A (ja) * | 1999-11-22 | 2001-07-17 | General Electric Co <Ge> | 磁気共鳴撮像時の分極用磁場の変動の補償 |
| WO2005003802A1 (en) * | 2003-07-07 | 2005-01-13 | Koninklijke Philips Electronics N.V. | A method of monitoring a magnetic field drift of a magnetic resonance imaging apparatus |
| JP2007260001A (ja) * | 2006-03-28 | 2007-10-11 | Univ Of Tsukuba | 磁気共鳴イメージング装置および撮像方法 |
| JP2010281812A (ja) * | 2009-04-27 | 2010-12-16 | Bruker Biospin Ag | 一定でないｒｆ位相を考慮しながらｎｍｒ送信周波数とｎｍｒ線の共鳴周波数を高精度で同期させる装置及びその同期システム |
| CN102866369A (zh) * | 2011-12-12 | 2013-01-09 | 中国科学院深圳先进技术研究院 | 磁共振的主磁场漂移矫正方法和*** |
-
1987
- 1987-11-30 JP JP62301955A patent/JPH01141656A/ja active Pending
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100747934B1 (ko) * | 1999-11-15 | 2007-08-08 | 제너럴 일렉트릭 캄파니 | 오차 정량화 방법 및 오차 정량화 시스템 |
| EP1102081A2 (en) * | 1999-11-15 | 2001-05-23 | General Electric Company | Method and apparatus to compensate for image artifacts caused by magnet vibration in an MR imaging system |
| EP1102081A3 (en) * | 1999-11-15 | 2003-05-28 | General Electric Company | Method and apparatus to compensate for image artifacts caused by magnet vibration in an MR imaging system |
| EP1102082A3 (en) * | 1999-11-15 | 2003-10-15 | General Electric Company | Method and apparatus for reducing image artifacts caused by magnet vibration in an MR imaging system |
| KR100758084B1 (ko) * | 1999-11-15 | 2007-09-11 | 제너럴 일렉트릭 캄파니 | 화상 아티팩트 감소 방법 및 화상 아티팩트 최소화 장치 |
| EP1102082A2 (en) * | 1999-11-15 | 2001-05-23 | General Electric Company | Method and apparatus for reducing image artifacts caused by magnet vibration in an MR imaging system |
| JP2001190519A (ja) * | 1999-11-22 | 2001-07-17 | General Electric Co <Ge> | 磁気共鳴撮像時の分極用磁場の変動の補償 |
| US7224164B2 (en) | 2003-07-07 | 2007-05-29 | Koninklijke Philips Electronics N.V. | Method of monitoring a magnetic field drift of a magnetic resonance imaging apparatus |
| WO2005003802A1 (en) * | 2003-07-07 | 2005-01-13 | Koninklijke Philips Electronics N.V. | A method of monitoring a magnetic field drift of a magnetic resonance imaging apparatus |
| JP2007260001A (ja) * | 2006-03-28 | 2007-10-11 | Univ Of Tsukuba | 磁気共鳴イメージング装置および撮像方法 |
| JP4576534B2 (ja) * | 2006-03-28 | 2010-11-10 | 国立大学法人 筑波大学 | 磁気共鳴イメージング装置および撮像方法 |
| JP2010281812A (ja) * | 2009-04-27 | 2010-12-16 | Bruker Biospin Ag | 一定でないｒｆ位相を考慮しながらｎｍｒ送信周波数とｎｍｒ線の共鳴周波数を高精度で同期させる装置及びその同期システム |
| CN102866369A (zh) * | 2011-12-12 | 2013-01-09 | 中国科学院深圳先进技术研究院 | 磁共振的主磁场漂移矫正方法和*** |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Morgan et al. | Correction of spatial distortion in EPI due to inhomogeneous static magnetic fields using the reversed gradient method | |
| US7218110B2 (en) | Method and apparatus for diffusion magnetic resonance imaging with the effects of eddy currents compensated | |
| US7689262B2 (en) | Method and magnetic resonance tomography apparatus for correcting changes in the basic magnetic field | |
| US20140232393A1 (en) | Mapping eddy current fields in mri system | |
| JP5051952B2 (ja) | 核スピントモグラフィ装置およびその作動方法 | |
| JPH07171122A (ja) | Ｅｐｉ及びｇｒａｓｅ ｍｒｉにおける読み出し傾斜磁界極性の補正方法 | |
| US6836113B2 (en) | Measurement and correction of gradient-induced cross-term magnetic fields in an EPI sequence | |
| US6906515B2 (en) | Magnetic resonance imaging device and method | |
| US5810727A (en) | Method of and device for the imaging of objects by means of magnetic resonance | |
| US4954779A (en) | Correction for eddy current caused phase degradation | |
| JPS63238856A (ja) | 核磁化分布決定方法及び装置 | |
| US20070170917A1 (en) | Magnetic resonance imaging with real-time magnetic filed mapping | |
| US6037775A (en) | Method and apparatus for magnetic field stabilization in a MRI system | |
| JPH0315455B2 (ja) | ||
| JPH01141656A (ja) | Ｎｍｒイメージング装置の磁場ドリフト補正方法 | |
| US4710715A (en) | Method of mapping magnetic field strength and tipping pulse accuracy of an NMR imager | |
| US7187171B2 (en) | Magnetic resonance imaging apparatus and its control method | |
| Butts et al. | Phase velocity mapping with a real time line scan technique | |
| US5905377A (en) | Method and apparatus for correcting gradient system and static magnetic field in magnetic resonance imaging | |
| JP3868754B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3112930B2 (ja) | 磁気共鳴装置 | |
| WO2012143540A1 (en) | Determining positions of a magnetic field probe in a magnetic resonance measurement | |
| JP2899821B2 (ja) | Ｍｒｉのｒｆ自動調整法 | |
| JPH0670911A (ja) | 磁気共鳴映像装置 | |
| JPH0244219B2 (ja) |