JPS61198043A - Ｍｒ−ｃｔ装置の静磁界変動検出方法 - Google Patents
Ｍｒ−ｃｔ装置の静磁界変動検出方法Info
- Publication number
- JPS61198043A JPS61198043A JP60040099A JP4009985A JPS61198043A JP S61198043 A JPS61198043 A JP S61198043A JP 60040099 A JP60040099 A JP 60040099A JP 4009985 A JP4009985 A JP 4009985A JP S61198043 A JPS61198043 A JP S61198043A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- static magnetic
- image
- static
- pulse 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
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
- G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
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)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
産業上の利用分野
この発明は、ＭＲ−ＣＴ型装置静磁界変動検出方法に関
する。
する。
従来の技術
ＭＲ−ＣＴ型装置、強力な静磁界中に置かれた被検体（
人体）の核磁気共鳴（ＮＭＲ）現象を利用して内部組織
の種々のパラメータに対応してイメージを作り出すもの
である。共鳴周波数は磁界に比例するので、静磁界に対
して傾斜磁界を重畳し、ＮＭＲ信号の周波数スペクトル
を求めることによって空間分布情報が得られる。イメー
ジング技術の１つとして、たとえばスピンワーブ法が知
られている（特開昭５４−１５８９８８）。ＭＲ−ＣＴ
型装置静磁界の強度および均一度はこの空間分布情報の
精度に直接影響を与える。
人体）の核磁気共鳴（ＮＭＲ）現象を利用して内部組織
の種々のパラメータに対応してイメージを作り出すもの
である。共鳴周波数は磁界に比例するので、静磁界に対
して傾斜磁界を重畳し、ＮＭＲ信号の周波数スペクトル
を求めることによって空間分布情報が得られる。イメー
ジング技術の１つとして、たとえばスピンワーブ法が知
られている（特開昭５４−１５８９８８）。ＭＲ−ＣＴ
型装置静磁界の強度および均一度はこの空間分布情報の
精度に直接影響を与える。
発明が解決しようとする問題点
ところで、従来では、ＭＲ−ＣＴ型装置静磁界の強度お
よび均一度は据付時に調整されているが、静磁界強度は
マグネット温度の影響を受は変動しやすい。静磁界の変
動は画像の歪み、ぼけ等を生じ画質に重大な影響を与え
る。
よび均一度は据付時に調整されているが、静磁界強度は
マグネット温度の影響を受は変動しやすい。静磁界の変
動は画像の歪み、ぼけ等を生じ画質に重大な影響を与え
る。
この発明は、撮像前などにおいて、随時、簡単に静磁界
の変動を検出する方法を提供することを目的とする。
の変動を検出する方法を提供することを目的とする。
問題点を解決するための手段
この発明のＭＲ−ＣＴ装置の静磁界変動検出方法では、
静磁界に１軸方向の傾Ｍ磁界を重畳させながらデータ収
集する第１のパルスシーケンスと、静磁界に上記ｌ軸方
向に関して反対方向の傾斜磁界を重畳させながらデータ
収集する第２のパルスシーケンスとが実行され、これら
第１．第２のパルスシーケンスで得たデータの相関関係
より静磁界の基準値からの変動を検出する。
静磁界に１軸方向の傾Ｍ磁界を重畳させながらデータ収
集する第１のパルスシーケンスと、静磁界に上記ｌ軸方
向に関して反対方向の傾斜磁界を重畳させながらデータ
収集する第２のパルスシーケンスとが実行され、これら
第１．第２のパルスシーケンスで得たデータの相関関係
より静磁界の基準値からの変動を検出する。
作 用
静磁界が基準値からずれていると、第１のパルスシーケ
ンスで得た周波数情報と第２のパルスシーケンスで得た
周波数情報はともに、その静磁界のずれに対応する周波
数だけ同じようにずれることになる。そこで、いずれか
一方の周波数情報を基準値の静磁界に対応する中心周波
数に関して反転し他方と比較するなどの両者の相関をと
れば、直ちに静磁界のずれを知ることができる。
ンスで得た周波数情報と第２のパルスシーケンスで得た
周波数情報はともに、その静磁界のずれに対応する周波
数だけ同じようにずれることになる。そこで、いずれか
一方の周波数情報を基準値の静磁界に対応する中心周波
数に関して反転し他方と比較するなどの両者の相関をと
れば、直ちに静磁界のずれを知ることができる。
実施例
第３図はＭＲ−ＣＴ装置の概要を示すものである。この
図で静磁界コイル１、Ｇｚ（Ｚ方向傾斜磁界）コイル２
、ＧＶ（Ｙ方向傾斜磁界）コイル３、Ｇｘ（Ｘ方向傾斜
磁界）コイル４に電源１１〜１４よりそれぞれ励磁電流
が供給される。ｘ。
図で静磁界コイル１、Ｇｚ（Ｚ方向傾斜磁界）コイル２
、ＧＶ（Ｙ方向傾斜磁界）コイル３、Ｇｘ（Ｘ方向傾斜
磁界）コイル４に電源１１〜１４よりそれぞれ励磁電流
が供給される。ｘ。
Ｙ、Ｚの各方向は８４図に示す通り１人体の体軸方向を
Ｚとし、Ｘ、Ｙの各方向はＸ−Ｙ平面がＺ軸に直角な平
面をなす方向とする。被検体である人体はこれら各コイ
ル１〜４によって形成される磁界中に置かれ、人体の周
囲にアンテナコイル５が配置される。
Ｚとし、Ｘ、Ｙの各方向はＸ−Ｙ平面がＺ軸に直角な平
面をなす方向とする。被検体である人体はこれら各コイ
ル１〜４によって形成される磁界中に置かれ、人体の周
囲にアンテナコイル５が配置される。
制御コンピュータ２１は、後述の本発明方法で検出した
静磁界変動にもとづき静磁界電源１１を制御するととも
に、スピンワープ法などの所定のパルスシーケンスを実
行するため波形発生回路２２を制御する。波形発生回路
２２から発生する波形によって傾斜磁界電源１２〜１４
およびＲＦＦ幅器・高周波電源２３が所定のシーケンス
に沿って制御される。ＲＦ増増幅器嵩高周波電源２３ら
発生するＲＦパルスはスイッチ回路２４を経てアンテナ
コイル５に送られ、人体に１８０°パルスや９０°パル
スなどが与えられる。
静磁界変動にもとづき静磁界電源１１を制御するととも
に、スピンワープ法などの所定のパルスシーケンスを実
行するため波形発生回路２２を制御する。波形発生回路
２２から発生する波形によって傾斜磁界電源１２〜１４
およびＲＦＦ幅器・高周波電源２３が所定のシーケンス
に沿って制御される。ＲＦ増増幅器嵩高周波電源２３ら
発生するＲＦパルスはスイッチ回路２４を経てアンテナ
コイル５に送られ、人体に１８０°パルスや９０°パル
スなどが与えられる。
人体からのＦＩＤ （自由銹導減衰）信号やスピンエコ
ー信号はアンテナコイル５で受信され、これらＮＭＲ信
号は切り換えられたスイッチ回路２４によりＲＦＦ幅器
２５に送られ、さらに位相弁別検波回路２６およびイン
ターフェイス２７を経、この過程でＡ／Ｄサンプリング
されてディジタルデータにされ、主コンピユータ２８に
取り込まれる。主コンピユータ２８とイメージプロセッ
サ２９によって高速フーリエ変換を主体とするデータ処
理が行なわれて画像が作られ、ディスプレイ装置で表示
される。
ー信号はアンテナコイル５で受信され、これらＮＭＲ信
号は切り換えられたスイッチ回路２４によりＲＦＦ幅器
２５に送られ、さらに位相弁別検波回路２６およびイン
ターフェイス２７を経、この過程でＡ／Ｄサンプリング
されてディジタルデータにされ、主コンピユータ２８に
取り込まれる。主コンピユータ２８とイメージプロセッ
サ２９によって高速フーリエ変換を主体とするデータ処
理が行なわれて画像が作られ、ディスプレイ装置で表示
される。
スピンワープ法などの２次元フーリエ変換法（フーリエ
拳ン゛イマトグラフィ）やン゛イマトグラフィ（投影復
元法）により第４図のＸ−Ｙ平面に関する断層像を得る
場合、第１図のパルスシーケンス（１）のように、９０
°パルスの後発生するＦＩＤ信号をＡ／Ｄサンプリング
するときにＧｘ（Ｘ方向傾斜磁界）を与えて、Ｘ−Ｙ平
面上の対象のＸ軸への投影データを得る。 Ｇｙ、 Ｇ
ｚ（Ｙ、Ｚ方向傾斜磁界）については省略しているが、
Ｇｚは選択照射法によりｘ−Ｙ平面に対応するＺ方向の
狭い範囲を励起するためＲＦパルスとともに与えられ、
ＧＶはスピンワープ法（あるいは２次元フーリエ変換法
）においてＹ方向の位相エンコーディングを行なうため
に与えられる。
拳ン゛イマトグラフィ）やン゛イマトグラフィ（投影復
元法）により第４図のＸ−Ｙ平面に関する断層像を得る
場合、第１図のパルスシーケンス（１）のように、９０
°パルスの後発生するＦＩＤ信号をＡ／Ｄサンプリング
するときにＧｘ（Ｘ方向傾斜磁界）を与えて、Ｘ−Ｙ平
面上の対象のＸ軸への投影データを得る。 Ｇｙ、 Ｇ
ｚ（Ｙ、Ｚ方向傾斜磁界）については省略しているが、
Ｇｚは選択照射法によりｘ−Ｙ平面に対応するＺ方向の
狭い範囲を励起するためＲＦパルスとともに与えられ、
ＧＶはスピンワープ法（あるいは２次元フーリエ変換法
）においてＹ方向の位相エンコーディングを行なうため
に与えられる。
静磁界変動を測定するためには、このような第１図のパ
ルスシーケンス（りに加えて第１図のパルスシーケンス
（２）が実行される。このパルスシーケンス（２）はパ
ルスシーケンス（１）とＸ方向傾斜磁界Ｇｘの向きが正
反対になっているだけで他のパラメータは全く同じであ
る。
ルスシーケンス（りに加えて第１図のパルスシーケンス
（２）が実行される。このパルスシーケンス（２）はパ
ルスシーケンス（１）とＸ方向傾斜磁界Ｇｘの向きが正
反対になっているだけで他のパラメータは全く同じであ
る。
ここで、静磁界とＸ方向傾斜磁界Ｇｘとの合成の様子を
見ると第２図のＡのようになる。静磁界は第４図のよう
にＺ方向に向きその強さが予め設定された基準値Ｈａで
あるとする− Ｘｓ、Ｘ２は設定されている関心領域の
Ｘ方向の境界であり。
見ると第２図のＡのようになる。静磁界は第４図のよう
にＺ方向に向きその強さが予め設定された基準値Ｈａで
あるとする− Ｘｓ、Ｘ２は設定されている関心領域の
Ｘ方向の境界であり。
Ｘｏはその中点である。Ｘ方向傾斜磁界Ｇｘの勾配をｈ
ｏ ／ （Ｘ２ Ｘｌ）７！ｌ：すれば、パルスシー
ケンス（１）においてはＸ方向の各位置での合成磁界強
度は第２図Ａの実線のようになる。したがって収集され
たデータをフーリエ変換して周波数情報を得ればそれは
第２図Ｂの上のようになる。
ｏ ／ （Ｘ２ Ｘｌ）７！ｌ：すれば、パルスシー
ケンス（１）においてはＸ方向の各位置での合成磁界強
度は第２図Ａの実線のようになる。したがって収集され
たデータをフーリエ変換して周波数情報を得ればそれは
第２図Ｂの上のようになる。
ｆｏはＨａによって決まる周波数であり、ｆｔ−ｆｏ
＝ｆｏ−ｆ□はｈＯによって決まる。パルスシーケンス
（２）でデータを収集すると、そのときの合成磁界強度
は第２図Ａの破線のようになるから、周波数情報は第２
図Ｂの下のようになる。この場合、両パルスシーケンス
で得た周波数情報（第２図Ｂの上と下）の関係は、パル
スシーケンス（２）の情報（第２図Ｂの下）を周波数ｆ
Ｏで左右反転するとパルスシーケンス（１）の情報（第
２図Ｂの上）と全く同じになるという関係である。
＝ｆｏ−ｆ□はｈＯによって決まる。パルスシーケンス
（２）でデータを収集すると、そのときの合成磁界強度
は第２図Ａの破線のようになるから、周波数情報は第２
図Ｂの下のようになる。この場合、両パルスシーケンス
で得た周波数情報（第２図Ｂの上と下）の関係は、パル
スシーケンス（２）の情報（第２図Ｂの下）を周波数ｆ
Ｏで左右反転するとパルスシーケンス（１）の情報（第
２図Ｂの上）と全く同じになるという関係である。
ところが静磁界が変動してその強度が基準値ＨｏからＨ
ｏ’にずれたとすると、このような関係は成立しない。
ｏ’にずれたとすると、このような関係は成立しない。
この場合、パルスシーケンス（１）では、合成磁界強度
は第２図Ｃの実線のように。
は第２図Ｃの実線のように。
周波数情報は第２図りの上のようになり、パルスシーケ
ンス（２）では１合成磁界強度は第２図Ｃの破線のよう
に、周波数情報は第２図りの下のようになる。したがっ
て、もはや、第２図りの下のような周波数情報をｆＯを
中心に左右反転しても、第２図りの上のような周波数情
報とは一致しない。そこで、第２図りの下の情報を反転
したものと上の反転しない情報との相互相関関数を計算
し、その最大値となる周波数の値からｆｏ’を求めれば
、ｆｏ’−ｆｏからＨａ’−Ｈａ 、すなわち静磁界の
変動量が分る。
ンス（２）では１合成磁界強度は第２図Ｃの破線のよう
に、周波数情報は第２図りの下のようになる。したがっ
て、もはや、第２図りの下のような周波数情報をｆＯを
中心に左右反転しても、第２図りの上のような周波数情
報とは一致しない。そこで、第２図りの下の情報を反転
したものと上の反転しない情報との相互相関関数を計算
し、その最大値となる周波数の値からｆｏ’を求めれば
、ｆｏ’−ｆｏからＨａ’−Ｈａ 、すなわち静磁界の
変動量が分る。
したがって、この検出された変動量に応じて静磁界電源
１１を制御することにより静磁界を補正することができ
る。
１１を制御することにより静磁界を補正することができ
る。
なお、上記では静磁界の変動を補正するようにしたが、
１回のスキャンで１ラインデータを収集するときに、上
記のようにして検出された静磁界の変動量により、この
各ラインの位置を補正するようにしてもよい。
１回のスキャンで１ラインデータを収集するときに、上
記のようにして検出された静磁界の変動量により、この
各ラインの位置を補正するようにしてもよい。
また、スピンワープ法では、Ｙ方向の傾斜磁界Ｇｙの向
きが第１ラインと第ｎライン（ｎはライン（ビュー）数
）、第２ラインと第（ｎ −１）ライン、・・・、でそ
れぞれ正反対になるから、それらを基にスキャン中の静
磁界変動を求めることもできる。
きが第１ラインと第ｎライン（ｎはライン（ビュー）数
）、第２ラインと第（ｎ −１）ライン、・・・、でそ
れぞれ正反対になるから、それらを基にスキャン中の静
磁界変動を求めることもできる。
さらに、実際のスキャンの直前、直後に上記のようなパ
ルスシーケンス（１） 、 （２）による方法で静磁界
変動を求め、その変動がある範囲内ならばそのまま画像
を構成し、範囲外ならデータの並び換えを行なって（こ
の場合データ量は半減するカリ画像を作ることもできる
。
ルスシーケンス（１） 、 （２）による方法で静磁界
変動を求め、その変動がある範囲内ならばそのまま画像
を構成し、範囲外ならデータの並び換えを行なって（こ
の場合データ量は半減するカリ画像を作ることもできる
。
発明の効果
この発明によれば、簡単に静磁界の変動を検出すること
ができ、静磁界が基準値からずれたままで撮像を行なっ
た場合の画像のずれや歪みやぼけを防止できる。
ができ、静磁界が基準値からずれたままで撮像を行なっ
た場合の画像のずれや歪みやぼけを防止できる。
第１図はこの発明の一実施例を説明するためのタイムチ
ャート、第２図は合成磁界強度と周波数情報との関係を
説明するためのもので、第２図Ａ、ＣはＸ軸での合成磁
界強度分布のグラフ、第２図Ｂ、Ｄは周波数情報の周波
数域での位置を示す図、第３図はＭＲ−ＣＴ装置の概要
を示すブロック図、第４図は人体に対する各方向を説明
するための模式図である。 １・・・静磁界コイル ２・・・Ｚ方向傾斜磁界コイル ３・・・Ｙ方向傾斜磁界コイル ４・・・Ｘ方向傾斜磁界コイル ５・・・アンテナコイル
ャート、第２図は合成磁界強度と周波数情報との関係を
説明するためのもので、第２図Ａ、ＣはＸ軸での合成磁
界強度分布のグラフ、第２図Ｂ、Ｄは周波数情報の周波
数域での位置を示す図、第３図はＭＲ−ＣＴ装置の概要
を示すブロック図、第４図は人体に対する各方向を説明
するための模式図である。 １・・・静磁界コイル ２・・・Ｚ方向傾斜磁界コイル ３・・・Ｙ方向傾斜磁界コイル ４・・・Ｘ方向傾斜磁界コイル ５・・・アンテナコイル
Claims (1)
- （１）静磁界に１軸方向の傾斜磁界を重畳させながらデ
ータ収集する第１のパルスシーケンスと、静磁界に上記
１軸方向に関して反対方向の傾斜磁界を重畳させながら
データ収集する第２のパルスシーケンスとを有し、これ
ら第１、第２のパルスシーケンスで得たデータの相関関
係より静磁界の基準値からの変動を検出することを特徴
とするＭＲ−ＣＴ装置の静磁界変動検出方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60040099A JPS61198043A (ja) | 1985-02-28 | 1985-02-28 | Ｍｒ−ｃｔ装置の静磁界変動検出方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60040099A JPS61198043A (ja) | 1985-02-28 | 1985-02-28 | Ｍｒ−ｃｔ装置の静磁界変動検出方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS61198043A true JPS61198043A (ja) | 1986-09-02 |
Family
ID=12571415
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP60040099A Pending JPS61198043A (ja) | 1985-02-28 | 1985-02-28 | Ｍｒ−ｃｔ装置の静磁界変動検出方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS61198043A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006116216A (ja) * | 2004-10-25 | 2006-05-11 | Toshiba Corp | 磁気共鳴イメージング装置および画像補正評価方法 |
| CN105806928A (zh) * | 2016-03-04 | 2016-07-27 | 中国海洋石油总公司 | 一种静磁场核磁效应分析方法 |
-
1985
- 1985-02-28 JP JP60040099A patent/JPS61198043A/ja active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006116216A (ja) * | 2004-10-25 | 2006-05-11 | Toshiba Corp | 磁気共鳴イメージング装置および画像補正評価方法 |
| CN105806928A (zh) * | 2016-03-04 | 2016-07-27 | 中国海洋石油总公司 | 一种静磁场核磁效应分析方法 |
| CN105806928B (zh) * | 2016-03-04 | 2019-02-26 | 中国海洋石油集团有限公司 | 一种静磁场核磁效应分析方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3431203B2 (ja) | 検査空間における磁場のシミング方法 | |
| US6700376B2 (en) | Method and apparatus for correcting static magnetic field using a pair of magnetic fields which are the same or different from each other in intensity and direction | |
| US4672320A (en) | Imaging apparatus and method using nuclear magnetic resonance | |
| US20050218894A1 (en) | Magnetic resonance imaging device | |
| JPH0236260B2 (ja) | ||
| US5195524A (en) | Flow imaging method by means of an MRI apparatus and apparatus for realizing same | |
| EP0204320B1 (en) | Magnetic resonance imaging method | |
| US4791369A (en) | Method of measuring magnetic field error of an NMR imaging apparatus and of correcting distortion caused by the error | |
| EP0124108B1 (en) | Correction circuit for a static magnetic field of an nmr apparatus and nmr apparatus for utilizing the same | |
| JP4138258B2 (ja) | 磁気共鳴イメージング装置 | |
| JP2000279390A (ja) | 磁気共鳴イメージング装置 | |
| US4710715A (en) | Method of mapping magnetic field strength and tipping pulse accuracy of an NMR imager | |
| JPH0454448B2 (ja) | ||
| JPS61198043A (ja) | Ｍｒ−ｃｔ装置の静磁界変動検出方法 | |
| JP2585278B2 (ja) | 核磁気共鳴を用いた検査装置 | |
| JP3018076B2 (ja) | 核磁気共鳴を用いた検査装置 | |
| JPH01221153A (ja) | Ｍｒｉ装置 | |
| JP2607466B2 (ja) | 核磁気共鳴を用いた検査装置 | |
| JPS62246356A (ja) | 核磁気共鳴を用いる検査装置 | |
| JPH049414B2 (ja) | ||
| JPH08215167A (ja) | 空間磁場分布測定方法及び空間磁場分布測定装置 | |
| JPH0244219B2 (ja) | ||
| JP2602226B2 (ja) | 核磁気共鳴イメージング装置 | |
| JPS60149953A (ja) | 核磁気共鳴を用いた検査装置 | |
| JPH0337046A (ja) | 磁気共鳴撮像装置の磁場状態検出方法 |