JP3701540B2 - 磁場測定方法、勾配コイル製造方法、勾配コイルおよび磁気共鳴撮影装置 - Google Patents
磁場測定方法、勾配コイル製造方法、勾配コイルおよび磁気共鳴撮影装置 Download PDFInfo
- Publication number
- JP3701540B2 JP3701540B2 JP2000117880A JP2000117880A JP3701540B2 JP 3701540 B2 JP3701540 B2 JP 3701540B2 JP 2000117880 A JP2000117880 A JP 2000117880A JP 2000117880 A JP2000117880 A JP 2000117880A JP 3701540 B2 JP3701540 B2 JP 3701540B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- gradient
- measurement
- coil
- space
- 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.)
- Expired - Fee Related
Links
- 238000002595 magnetic resonance imaging Methods 0.000 title claims description 14
- 238000004519 manufacturing process Methods 0.000 title claims description 9
- 238000000691 measurement method Methods 0.000 title claims description 8
- 238000005259 measurement Methods 0.000 claims description 99
- 238000000034 method Methods 0.000 claims description 37
- 230000005284 excitation Effects 0.000 claims description 24
- 230000003068 static effect Effects 0.000 claims description 18
- 239000000523 sample Substances 0.000 description 27
- 230000006870 function Effects 0.000 description 23
- 238000012545 processing Methods 0.000 description 17
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 15
- 238000010586 diagram Methods 0.000 description 13
- 238000013480 data collection Methods 0.000 description 9
- 238000003384 imaging method Methods 0.000 description 9
- 230000005415 magnetization Effects 0.000 description 7
- 238000013461 design Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 230000002093 peripheral effect Effects 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- ARUVKPQLZAKDPS-UHFFFAOYSA-L copper(II) sulfate Chemical compound [Cu+2].[O-][S+2]([O-])([O-])[O-] ARUVKPQLZAKDPS-UHFFFAOYSA-L 0.000 description 2
- 239000004615 ingredient Substances 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 239000007864 aqueous solution Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- GGVOVPORYPQPCE-UHFFFAOYSA-M chloronickel Chemical compound [Ni]Cl GGVOVPORYPQPCE-UHFFFAOYSA-M 0.000 description 1
- 229910000365 copper sulfate Inorganic materials 0.000 description 1
- 229910000366 copper(II) sulfate Inorganic materials 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 239000000696 magnetic material Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- 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/56518—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Signal Processing (AREA)
- General Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、磁場測定方法、勾配コイル（ｃｏｉｌ）製造方法、勾配コイルおよび磁気共鳴撮影装置に関し、特に、勾配磁場の印加停止後に残留する磁場を測定する方法、残留磁場の高次項が小さい勾配コイルとその製造方法、並びに、そのような勾配コイルを有する磁気共鳴撮影装置に関する。
【０００２】
【従来の技術】
磁気共鳴撮影（ＭＲＩ：Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）装置では、マグネットシステム（ｍａｇｎｅｔ ｓｙｓｔｅｍ）の内部空間すなわち静磁場を形成した空間に撮影する対象を搬入し、勾配磁場および高周波磁場を印加して対象内に磁気共鳴信号を発生させ、その受信信号に基づいて断層像を生成（再構成）する。
【０００３】
静磁場を発生するのに永久磁石を用いるマグネットシステムでは、互いに対向する１対の永久磁石の先端に静磁場空間の磁束分布を均一化するための磁極片をそれぞれ設け、勾配磁場発生用の勾配コイルはこの磁極片の磁極面に沿って設けられる。
【０００４】
【発明が解決しようとする課題】
上記のようなマグネットシステムでは、勾配コイルが磁極片に近接しているので勾配磁場により磁極片が磁化され、その残留磁化が形成する磁場により、スピン（ｓｐｉｎ）の位相が、あたかも時定数が極めて長い渦電流が存在するかのような影響を受ける。このため、精密な位相制御が要求される例えばファースト・スピンエコー（ＦＳＥ：Ｆａｓｔ Ｓｐｉｎ Ｅｃｈｏ）法等による撮影に支障をきたす。
【０００５】
これまでは、残留磁化による磁場を測定する適切な方法がなかったので、その影響を除去する方策を立てることができなかった。
そこで、本発明の課題は、勾配磁場の印加停止後に残留する磁場を測定する方法、測定によって判明した残留磁場の特性を反映させて勾配コイルを製造する方法、その方法によって製造された勾配コイル、および、そのような勾配コイルを有する磁気共鳴撮影装置を実現することである。
【０００６】
【課題を解決するための手段】
（１）上記の課題を解決するための１つの観点での発明は、測定空間に想定した球の表面の測定点に磁気共鳴可能な試料を配置し、前記測定空間に勾配磁場を印加し、前記勾配磁場の印加停止後にＲＦ励起を行って前記試料が生じるＦＩＤ信号を測定し、前記ＦＩＤ信号の位相の微分値に基づいて前記測定点における磁場強度を求め、前記磁場強度を前記測定空間における磁場強度分布を表す球面調和関数にフィッティングすることを特徴とする磁場測定方法である。
【０００７】
この観点での発明では、残留磁化が生じる磁場を、球面上の測定点に配置した試料のＦＩＤ信号を利用して測定し、測定値を球面調和関数にフィッティングして、残留磁場の特性を表す球面調和関数を特定する。
【０００８】
（２）上記の課題を解決するための他の観点での発明は、測定空間に想定した球の表面の測定点に磁気共鳴可能な試料を配置し、前記球の中心に磁気共鳴可能な試料を配置し、前記測定空間に勾配磁場を印加し、前記勾配磁場の印加停止後にＲＦ励起を行って前記試料が生じるＦＩＤ信号を測定し、前記測定点および前記球の中心の試料がそれぞれ生じるＦＩＤ信号の位相の微分値の差に基づいて前記測定点における磁場強度を求め、前記磁場強度を前記測定空間における磁場強度分布を表す球面調和関数にフィッティングすることを特徴とする磁場測定方法である。
【０００９】
この観点での発明では、球の中心に配置した試料のＦＩＤ信号を測定してこれをリファレンスとする。
（３）上記の課題を解決するための他の観点での発明は、測定空間に想定した球の表面の測定点に磁気共鳴可能な試料を配置し、１つの極性における最大勾配から反対極性における最大勾配まで順次に勾配を変化させ次いで前記反対極性における最大勾配から前記１つの極性における最大勾配まで順次に勾配を変化させながら間欠的に勾配磁場を印加し、前記勾配磁場の間欠期間中にＲＦ励起を行って前記試料が生じるＦＩＤ信号を測定し、前記ＦＩＤ信号の位相の微分値を求め、前記勾配の一巡に伴う前記微分値のヒステリシスに基づいて前記測定点における磁場強度を求め、前記磁場強度を前記測定空間における磁場強度分布を表す球面調和関数にフィッティングすることを特徴とする磁場測定方法である。
【００１０】
この観点での発明では、勾配磁場を１つの極性における最大勾配と反対極性における最大勾配の間で往復変化させる。これにより、残留磁化の履歴を一定にして磁場測定を行うことができる。
【００１１】
（４）上記の課題を解決するための他の観点での発明は、測定空間に想定した球の表面の測定点に磁気共鳴可能な試料を配置し、前記球の中心に磁気共鳴可能な試料を配置し、１つの極性における最大勾配から反対極性における最大勾配まで順次に勾配を変化させ次いで前記反対極性における最大勾配から前記１つの極性における最大勾配まで順次に勾配を変化させながら間欠的に勾配磁場を印加し、前記勾配磁場の間欠期間中にＲＦ励起を行って前記試料が生じるＦＩＤ信号を測定し、前記測定点および前記球の中心の試料がそれぞれ生じるＦＩＤ信号の位相の微分値の差を求め、前記勾配の一巡に伴う前記微分値の差のヒステリシスに基づいて前記測定点における磁場強度を求め、前記磁場強度を前記測定空間における磁場強度分布を表す球面調和関数にフィッティングすることを特徴とする磁場測定方法である。
【００１２】
この観点での発明では、球の中心に配置した試料のＦＩＤ信号を測定してこれをリファレンスとする。
（５）上記の課題を解決するための他の観点での発明は、前記測定を複数の測定点について逐次に行うことを特徴とする（１）ないし（４）のうちのいずれか１つに記載の磁場測定方法である。
【００１３】
この観点での発明では、複数の測定点について逐次に測定を行うので、同一試料を使い回しすることができる。
（６）上記の課題を解決するための他の観点での発明は、勾配コイルを製造するにあたり、空間に想定した球面上に前記勾配コイルが発生する磁界を計算し、前記磁界を球面調和関数にフィッティングし、前記フィッティングした球面調和関数における高次項のうち、（１）ないし（５）のうちのいずれか１つに記載の方法で測定した磁場を表す球面調和関数における最大の高次項に相当する高次項が小さくなるように、コイルの電流経路を決定することを特徴とする勾配コイル製造方法である。
【００１４】
この観点での発明では、測定によって判明した残留磁場の最大の高次項に相当する高次項が小さくなるように、勾配コイルの電流経路を決めるので、影響が小さい残留磁化を生じる勾配コイルを製造することができる。
【００１５】
（７）上記の課題を解決するための他の観点での発明は、電流経路を流れる電流によって勾配磁場を発生する勾配コイルであって、前記電流経路が下記の手順で決定されたことを特徴とする勾配コイルである。
【００１６】
記
（イ）空間に想定した球面上に前記勾配コイルが発生する磁界を計算する。
（ロ）前記磁界を球面調和関数にフィッティングする。
【００１７】
（ハ）前記フィッティングした球面調和関数における高次項のうち、（１）ないし（５）のうちのいずれか１つに記載の方法で測定した磁場を表す球面調和関数における最大の高次項に相当する高次項が小さくなるように、コイルの電流経路を決定する。
【００１８】
この観点での発明では、測定によって判明した残留磁場の最大の高次項に相当する高次項が小さくなるように、勾配コイルの電流経路が決められているので、影響が小さい残留磁化を生じる勾配コイルとなる。
【００１９】
（８）上記の課題を解決するための他の観点での発明は、静磁場、勾配磁場および高周波磁場を用いて獲得した磁気共鳴信号に基づいて画像を構成する磁気共鳴撮影装置であって、前記勾配磁場を発生する勾配コイルとして、下記の手順で決定された電流経路を有する勾配コイルを具備することを特徴とする磁気共鳴撮影装置である。
【００２０】
記
（イ）空間に想定した球面上に前記勾配コイルが発生する磁界を計算する。
（ロ）前記磁界を球面調和関数にフィッティングする。
【００２１】
（ハ）前記フィッティングした球面調和関数における高次項のうち、（１）ないし（５）のうちのいずれか１つに記載の方法で測定した磁場を表す球面調和関数における最大の高次項に相当する高次項が小さくなるように、コイルの電流経路を決定する。
【００２２】
この観点での発明では、測定によって判明した残留磁場の最大の高次項に相当する高次項が小さくなるように勾配コイルの電流経路が決められているので、影響が小さい残留磁化を生じる勾配コイルを用いて適正な磁気共鳴撮影をすることができる。
【００２３】
【発明の実施の形態】
以下、図面を参照して本発明の実施の形態を詳細に説明する。なお、本発明は実施の形態に限定されるものではない。図１に磁気共鳴撮影装置のブロック（ｂｌｏｃｋ）図を示す。本装置は本発明の実施の形態の一例である。本装置の構成によって、本発明の装置に関する実施の形態の一例が示される。本装置の動作によって、本発明の方法に関する実施の形態の一例が示される。
【００２４】
図１に示すように、本装置はマグネットシステム１００を有する。マグネットシステム１００は、主磁場マグネット部１０２、勾配コイル部１０６およびＲＦ（ｒａｄｉｏ ｆｒｅｑｕｅｎｃｙ）コイル部１０８を有する。これら主磁場マグネット部１０２および各コイル部は、いずれも空間を挟んで互いに対向する１対のものからなる。また、いずれも概ね円板状の形状を有し中心軸を共有して配置されている。マグネットシステム１００の内部空間（ボア：ｂｏｒｅ）に、対象３００がクレードル（ｃｒａｄｌｅ）５００に搭載されて図示しない搬送手段により搬入および搬出される。
【００２５】
主磁場マグネット部１０２はマグネットシステム１００の内部空間に静磁場を形成する。静磁場の方向は概ね対象３００の体軸方向と直交する。すなわちいわゆる垂直磁場を形成する。主磁場マグネット部１０２は例えば永久磁石等を用いて構成される。なお、永久磁石に限らず超伝導電磁石あるいは常伝導電磁石等を用いて構成しても良いのはもちろんである。
【００２６】
勾配コイル部１０６は静磁場強度に勾配を持たせるための勾配磁場を生じる。発生する勾配磁場は、スライス（ｓｌｉｃｅ）勾配磁場、リードアウト（ｒｅａｄ ｏｕｔ）勾配磁場およびフェーズエンコード（ｐｈａｓｅ ｅｎｃｏｄｅ）勾配磁場の３種であり、これら３種類の勾配磁場に対応して勾配コイル部１０６は図示しない３系統の勾配コイルを有する。
【００２７】
３系統の勾配コイルは、互いに直交する３方向において静磁場にそれぞれ勾配を付与するための３つの勾配磁場をそれぞれ発生する。３方向のうちの１つは静磁場の方向（垂直方向）であり、通常これをｚ方向とする。他の１つは水平方向であり、通常これをｙ方向とする。残りの１つはｚ，ｙ方向に垂直な方向であり、通常これをｘ方向とする。x方向は垂直面内でｚ方向に垂直であり、水平面内でｙ方向に垂直である。以下、ｘ，ｙ，ｚを勾配軸ともいう。
【００２８】
ｘ，ｙ，ｚはいずれもスライス勾配の軸とすることができる。いずれか１つをスライス勾配の軸としたとき、残り２つのうちの１つをフェーズエンコード勾配の軸とし、他をリードアウト勾配の軸とする。３系統の勾配コイルについては後にあらためて説明する。
【００２９】
ＲＦコイル部１０８は静磁場空間に対象３００の体内のスピンを励起するためのＲＦ励起信号を送信する。ＲＦコイル部１０８は、また、励起されたスピンが生じる磁気共鳴信号を受信する。ＲＦコイル部１０８は図示しない送信用のコイルおよび受信用のコイルを有する。送信用のコイルおよび受信用のコイルは、同じコイルを兼用するかあるいはそれぞれ専用のコイルを用いる。
【００３０】
勾配コイル部１０６には勾配駆動部１３０が接続されている。勾配駆動部１３０は勾配コイル部１０６に駆動信号を与えて勾配磁場を発生させる。勾配駆動部１３０は、勾配コイル部１０６における３系統の勾配コイルに対応して、図示しない３系統の駆動回路を有する。
【００３１】
ＲＦコイル部１０８にはＲＦ駆動部１４０が接続されている。ＲＦ駆動部１４０はＲＦコイル部１０８に駆動信号を与えてＲＦ励起信号を送信し、対象３００の体内のスピンを励起する。
【００３２】
ＲＦコイル部１０８には、また、データ（ｄａｔａ）収集部１５０が接続されている。データ収集部１５０はＲＦコイル部１０８が受信した受信信号を取り込み、それをビューデータ（ｖｉｅｗ ｄａｔａ）として収集する。
【００３３】
勾配駆動部１３０、ＲＦ駆動部１４０およびデータ収集部１５０には制御部１６０が接続されている。制御部１６０は、勾配駆動部１３０ないしデータ収集部１５０をそれぞれ制御して撮影を遂行する。
【００３４】
データ収集部１５０の出力側はデータ処理部１７０に接続されている。データ処理部１７０は、例えばコンピュータ（ｃｏｍｐｕｔｅｒ）等を用いて構成される。データ処理部１７０は図示しないメモリ（ｍｅｍｏｒｙ）を有する。メモリはデータ処理部１７０用のプログラムおよび各種のデータを記憶している。本装置の機能は、データ処理部１７０がメモリに記憶されたプログラムを実行することによりを実現される。
【００３５】
データ処理部１７０は、データ収集部１５０から取り込んだデータをメモリに記憶する。メモリ内にはデータ空間が形成される。データ空間は２次元フ−リエ（Ｆｏｕｒｉｅｒ）空間を構成する。データ処理部１７０は、これら２次元フ−リエ空間のデータを２次元逆フ−リエ変換して対象３００の画像を生成（再構成）する。２次元フ−リエ空間をｋスペース（ｋ−ｓｐａｃｅ）ともいう。
【００３６】
データ処理部１７０は制御部１６０に接続されている。データ処理部１７０は制御部１６０の上位にあってそれを統括する。データ処理部１７０には、また、表示部１８０および操作部１９０が接続されている。表示部１８０は、グラフィックディスプレー（ｇｒａｐｈｉｃ ｄｉｓｐｌａｙ）等で構成される。操作部１９０はポインティングデバイス（ｐｏｉｎｔｉｎｇ ｄｅｖｉｃｅ）を備えたキーボード（ｋｅｙｂｏａｒｄ）等で構成される。
【００３７】
表示部１８０は、データ処理部１７０から出力される再構成画像および各種の情報を表示する。操作部１９０は、操作者によって操作され、各種の指令や情報等をデータ処理部１７０に入力する。操作者は表示部１８０および操作部１９０を通じてインタラクティブ（ｉｎｔｅｒａｃｔｉｖｅ）に本装置を操作する。
【００３８】
図２に、本装置で撮影を行うときのパルスシーケンス（ｐｕｌｓｅ ｓｅｑｕｅｎｃｅ）の一例を示す。このパルスシーケンスは、グラディエントエコー（ＧＲＥ：Ｇｒａｄｉｅｎｔ Ｅｃｈｏ）法のパルスシーケンスである。
【００３９】
すなわち、（１）はＧＲＥ法におけるＲＦ励起用のα°パルスのシーケンスであり、（２）、（３）、（４）および（５）は、同じくそれぞれ、スライス勾配Ｇｓ、リードアウト勾配Ｇｒ、フェーズエンコード勾配ＧｐおよびグラディエントエコーＭＲのシーケンスである。なお、α°パルスは中心信号で代表する。パルスシーケンスは時間軸ｔに沿って左から右に進行する。
【００４０】
同図に示すように、α°パルスによりスピンのα°励起が行われる。フリップアングル（ｆｌｉｐ ａｎｇｌｅ）α°は９０°以下である。このときスライス勾配Ｇｓが印加され所定のスライスについての選択励起が行われる。
【００４１】
α°励起後、フェーズエンコード勾配Ｇｐによりスピンのフェーズエンコードが行われる。次に、リードアウト勾配Ｇｒにより先ずスピンをディフェーズ（ｄｅｐｈａｓｅ）し、次いでスピンをリフェーズ（ｒｅｐｈａｓｅ）して、グラディエントエコーＭＲを発生させる。グラディエントエコーＭＲの信号強度は、α°励起からエコータイム（ｅｃｈｏ ｔｉｍｅ）ＴＥ後の時点で最大となる。グラディエントエコーＭＲはデータ収集部１５０によりビューデータとして収集される。
【００４２】
このようなパルスシーケンスが周期ＴＲ（ｒｅｐｅｔｉｔｉｏｎ ｔｉｍｅ）で６４〜５１２回繰り返される。繰り返しのたびにフェーズエンコード勾配Ｇｐを変更し、毎回異なるフェーズエンコードを行う。これによって、ｋスペースを埋める６４〜５１２ビューのビューデータが得られる。
【００４３】
磁気共鳴撮影用パルスシーケンスの他の例を図３に示す。このパルスシーケンスは、スピンエコー（ＳＥ：Ｓｐｉｎ Ｅｃｈｏ）法のパルスシーケンスである。
【００４４】
すなわち、（１）はＳＥ法におけるＲＦ励起用の９０°パルスおよび１８０°パルスのシーケンスであり、（２）、（３）、（４）および（５）は、同じくそれぞれ、スライス勾配Ｇｓ、リードアウト勾配Ｇｒ、フェーズエンコード勾配ＧｐおよびスピンエコーＭＲのシーケンスである。なお、９０°パルスおよび１８０°パルスはそれぞれ中心信号で代表する。パルスシーケンスは時間軸ｔに沿って左から右に進行する。
【００４５】
同図に示すように、９０°パルスによりスピンの９０°励起が行われる。このときスライス勾配Ｇｓが印加され所定のスライスについての選択励起が行われる。９０°励起から所定の時間後に、１８０°パルスによる１８０°励起すなわちスピン反転が行われる。このときもスライス勾配Ｇｓが印加され、同じスライスについての選択的反転が行われる。
【００４６】
９０°励起とスピン反転の間の期間に、リードアウト勾配Ｇｒおよびフェーズエンコード勾配Ｇｐが印加される。リードアウト勾配Ｇｒによりスピンのディフェーズが行われる。フェーズエンコード勾配Ｇｐによりスピンのフェーズエンコードが行われる。
【００４７】
スピン反転後、リードアウト勾配ＧｒでスピンをリフェーズしてスピンエコーＭＲを発生させる。スピンエコーＭＲの信号強度は、９０°励起からＴＥ後の時点で最大となる。スピンエコーＭＲはデータ収集部１５０によりビューデータとして収集される。このようなパルスシーケンスが周期ＴＲで６４〜５１２回繰り返される。繰り返しのたびにフェーズエンコード勾配Ｇｐを変更し、毎回異なるフェーズエンコードを行う。これによって、ｋスペースを埋める６４〜５１２ビューのビューデータが得られる。
【００４８】
なお、撮影に用いるパルスシーケンスはＧＲＥ法またはＳＥ法に限るものではなく、例えば、ＦＳＥ（Ｆａｓｔ Ｓｐｉｎ Ｅｃｈｏ）法、ファーストリカバリＦＳＥ（Ｆａｓｔ Ｒｅｃｏｖｅｒｙ Ｆａｓｔ Ｓｐｉｎ Ｅｃｈｏ）法、エコープラナー・イメージング（ＥＰＩ：Ｅｃｈｏ Ｐｌａｎａｒ Ｉｍａｇｉｎｇ）等、他の適宜の技法のものであって良い。
【００４９】
データ処理部１７０は、ｋスペースのビューデータを２次元逆フ−リエ変換して対象３００の断層像を再構成する。再構成した画像はメモリに記憶し、また、表示部１８０で表示する。
【００５０】
図４に、勾配コイル部１０６付近のマグネットシステム１００の構造を断面図により模式的に示す。同図において、Ｏは静磁場の中心すなわちマグネットセンター（ｍａｇｎｅｔ ｃｅｎｔｅｒ）であり、ｘ，ｙ，ｚは前述した３方向である。
【００５１】
マグネットセンターＯを中心とする半径Ｒの球形領域ＳＶ（ｓｐｈｅｒｉｃ ｖｏｌｕｍｅ）が撮影領域であり、マグネットシステム１００はこのＳＶにおいて静磁場および勾配磁場が所定の精度を持つように構成される。
【００５２】
１対の主磁場マグネット部１０２は互いに対向する１対の磁極片（ポールピース：ｐｏｌｅ ｐｉｅｃｅ）２０２を有する。磁極片２０２は例えば軟鉄等の高透磁率の磁性材料で構成され、静磁場空間における磁束分布を均一化する働きをする。
【００５３】
磁極片２０２は概ね円板形状を成すが、周縁部が板面に垂直な方向（ｚ方向）、すなわち、磁極片２０２同士が互いに対向する方向に突出している。突出部は磁極片２０２の周縁部における磁束密度の低下を補う働きをする。
【００５４】
突出した周縁部の内側に形成される磁極片２０２の凹部に、勾配コイル部１０６が設けられる。勾配コイル部１０６は、Ｘコイル２０４、Ｙコイル２０６およびＺ２０８コイルを有する。
【００５５】
Ｘコイル２０４、Ｙコイル２０６およびＺコイル２０８は、それぞれ、本発明の勾配コイルの実施の形態の一例である。各コイルはいずれも概ね円板形状を成し、磁極片２０２の磁極面に、図示しない適宜の取り付け手段により、順次に層を成すように取り付けられている。
【００５６】
図５に、Ｘコイル２０４、Ｙコイル２０６およびＺ２０８コイルの電流経路のパターン（ｐａｔｔｅｒｎ）を略図によって示す。同図に示すように、Ｘコイル２０４は、円の中心に近い部分ではｙ方向に平行な直線状の複数の主電流経路（メインパス：ｍａｉｎ ｐａｓｓ）を有する。これらメインパスは円の中心を通るｙ軸に関して対称である。メインパスの帰路（リターンパス：ｒｅｔｕｒｎ ｐａｓｓ）は円周に沿って形成される。最外側のリターンパスの半径すなわちＸコイル２０４の外径はｒ００である。
【００５７】
Ｙコイル２０６は、円の中心に近い部分ではｘ方向に平行な直線状の複数のメインパスを有する。これらメインパスは円の中心を通るｘ軸に関して対称である。メインパスのリターンパスは円周に沿って形成される。最外側のリターンパスの半径すなわちＹコイル２０６の外径はｒ００である。
【００５８】
Ｚコイル２０８は同心円を成す複数の電流経路を有する。これら電流経路は全てがメインパスである。各メインパスの半径は内側から順にｒ１，ｒ２，・・・，ｒＭである。ｒＭはＺコイル２０８の外径となる。
【００５９】
このようなマグネットシステム１００について、撮影領域すなわち球形領域ＳＶにおける残留磁場を測定する方法を説明する。先ず、球形領域ＳＶの表面に複数の測定点を設定する。測定点Ｐを、図６に示すように、３次元の極座標で表すとした場合、複数の測定点は、半径ｒを一定とし、角度θ，φの組み合わせを複数通りに設定することにより行われる。角度θはいわば緯度に相当し、角度φはいわば経度に相当する。
【００６０】
角度θは、図７（ａ）に示すように、例えば２２．５度刻みで設定する。これによって７つの緯度が設定される。角度φは、同図の（ｂ）に示すように、例えば３０度刻みで設定する。これによって１２の経度が設定される。
【００６１】
このような緯度と経度の交点を測定点Ｐとする。したがって、球面上の測定点数は８４となる。これらの測定点に加えて、マグネットセンターＯにも測定点Ｐ０を設定する。測定点Ｐ０はリファレンス（ｒｅｆｅｒｅｎｃｅ）用である。リファレンスが不要な場合は測定点Ｐ０は省略可能である。
【００６２】
測定点には、磁気共鳴信号を生じる試料に信号検出用のＲＦコイルを取り付けたプローブ（ｐｒｏｂｅ）が設置される。試料としては、例えば、硫酸銅（ＣｕＳＯ４）や塩化ニッケル（ＮｉＣｌ）等の水溶液を微小な容器に封入したものが用いられる。試料の体積は、撮影空間の１ボクセル（ｖｏｘｅｌ）に相当する大きさとする。
【００６３】
そのようなプローブを適宜の支持手段で支持して上記の測定点に順次に位置決めし、球面上の各測定点およびマグネットセンターＯにおける磁場を測定する。プローブは少なくとも２つ用い、球面上の測定点とマグネットセンターＯの磁場を同時に測定する。リファレンスを必要としないときは、球面上の測定点だけで測定する。
【００６４】
図８に、磁場測定用のパルスシーケンスを示す。同図において（１）は勾配磁場のシーケンス、（２）はＲＦ励起のシーケンス、（３）は磁気共鳴信号すなわちＦＩＤ（Ｆｒｅｅ Ｉｎｄｕｃｔｉｏｎ Ｄｅｃａｙ）信号のシーケンスである。
【００６５】
同図に示すように、１つの軸に勾配磁場Ｇを印加し、それが終了した後にＲＦ励起を行う。ＲＦ励起に基づきプローブ中の試料がＦＩＤ信号を生じる。このＦＩＤ信号を、プローブに備わるＲＦコイルで受信し、受信データをデータ収集部１５０で収集する。
【００６６】
勾配磁場Ｇが例えば図９に実線で示すようなものであるとしたとき、それによって磁極片２０２が磁化されることにより、勾配磁場の印加をやめた後にも、例えば一点鎖線で示すような残留磁場が存在する。このため、図８に斜線で示すように、あたかも勾配磁場が残留するかのような状態になる。
【００６７】
このような残留磁場の下で測定されるＦＩＤ信号は次式で表される。
【００６８】
【数１】
【００６９】
ここで、
ρ：スピン密度
γ：磁界回転比
Ｇｒｍ：磁場強度
（１）式に示すように、磁場強度ＧｒｍはＦＩＤ信号の位相に影響を与える。そこで、複素数で表現したＦＩＤ信号
【００７０】
【数２】
【００７１】
から、
【００７２】
【数３】
【００７３】
により位相を求め、その時間微分を用いて
【００７４】
【数４】
【００７５】
により磁場強度を求める。これによって、測定点Ｐにおける磁場の測定値が得られる。
実際には、測定空間に静磁場が存在するので、マグネットセンターＯで測定したＦＩＤ信号の位相φ０（ｔ）の微分値を求め、
【００７６】
【数５】
【００７７】
により磁場強度を求める。これにより、静磁場の強度変動による測定誤差を除去することができる。以上の計算は、データ処理部１７０によって行う。
勾配磁場の印加の経過を反映して、磁極片２０２では磁化の履歴が生じている。磁化の履歴はマイナー・ヒステリシスループ（ｍｉｎｏｒ ｈｙｓｔｅｒｅｓｉｓ ｌｏｏｐ）を描く。このような履歴を持つものについて常に同一の条件で測定を行うために、磁場測定は次のようにして行われる。
【００７８】
図１０に、測定時の１つの測定点あたりの勾配磁場印加のシーケンスを示す。同図に示すように、勾配磁場は、その強度を順次変化させながら間欠的に印加する。すなわち、例えば、シーケンスの前半は、本装置で発生可能な最大の負の勾配−Ｇｍａｘから最大の正の勾配＋Ｇｍａｘまで順次に変化させながら印加し、後半は、最大の正の勾配＋Ｇｍａｘから最大の負の勾配−Ｇｍａｘまで順次に変化させながら印加する。このように勾配の強度を一巡させることにより、シーケンスを開始する前と後における磁極片２０２の磁化状態を、同一にすることができる。
【００７９】
勾配印加の間欠期間中は斜線で概念的に示すように勾配磁場が残留する。この期間中に、前述のようなＲＦ励起とＦＩＤ測定を毎回行い、各回の測定信号に基づいて、上記のようにしてそれぞれ磁場強度Ｇｒｍを求める。これによって、１つの測定点あたり例えば１５個の磁場測定値が得られる。
【００８０】
このような磁場測定値について、印加した勾配との対応を求めると、例えば図１１に示すようなヒステリシスループが得られる。そこで、このヒステリシスループのＧｒｍ軸方向の幅、すなわち、ヒステリシスループが横切るＧｒｍ軸上の２箇所の値の差を測定点Ｐの磁場強度Ｂｒ（ｋ）とする。
【００８１】
このようにして求めた８４個の測定点Ｐの磁場測定値を球面調和関数にフィッティング（ｆｉｔｔｉｎｇ）して、残留磁場の空間成分を求める。フィッティングはデータ処理部１７０により行う。以下、磁場測定値を単に測定値ともいう。
【００８２】
球面調和関数展開式は次式で与えられる。
【００８３】
【数６】
【００８４】
ここで、
【００８５】
【数７】
【００８６】
はルジャンドル（Ｌｅｇｅｎｄｒｅ）の倍関数、
【００８７】
【数８】
【００８８】
は定数である。
フィッティングは、例えば次のようにして行われる。
１． ｍ＝０の場合
この場合、（６）式は次のようになる。
【００８９】
【数９】
【００９０】
ｚ座標が異なる測定点での測定値を
【００９１】
【数１０】
【００９２】
とすると次式が成り立つ。
【００９３】
【数１１】
【００９４】
この方程式を解くことにより、
【００９５】
【数１２】
【００９６】
成分が得られる。これらは通常
【００９７】
【数１３】
【００９８】
成分と呼ばれる。
２． ｍ≠０の場合
例えばｍ＝１のとき（６）式は次式のようになる。
【００９９】
【数１４】
【０１００】
球面上のあるｘｙ平面内の測定値を
【０１０１】
【数１５】
【０１０２】
とすると、円周上の測定点における測定値は、
【０１０３】
【数１６】
【０１０４】
このとき、
【０１０５】
【数１７】
【０１０６】
【数１８】
【０１０７】
各ｘｙ平面について
【０１０８】
【数１９】
【０１０９】
を求めれば次式が成り立つ。
【０１１０】
【数２０】
【０１１１】
この方程式を解けば、
【０１１２】
【数２１】
【０１１３】
が求まる。これらは通常
【０１１４】
【数２２】
【０１１５】
と呼ばれる。
３．ｍ＝２の場合
この場合は（１０）式において
【０１１６】
【数２３】
【０１１７】
とすれば、同様に
【０１１８】
【数２４】
【０１１９】
すなわち
【０１２０】
【数２５】
【０１２１】
を求めることができる。
図１２に、以上のような複数の測定点における磁場測定とそれに基づく球面調和関数のフィッティングのフロー（ｆｌｏｗ）図を示す。
【０１２２】
同図に示すように、ステップ（ｓｔｅｐ）７０２で、球面上の測定点を設定する。測定点は前述の要領で設定される。測定点はＬ個である。
次に、ステップ７０４で、球の中心すなわちマグネットセンターＯと球面上の１つの測定点にそれぞれプローブをセット（ｓｅｔ）する。
【０１２３】
次に、ステップ７０６で、図１０に示したような測定シーケンスを実行して、２つのプローブのＦＩＤ信号を測定し、それらの位相の微分値の差に基づいて磁場測定値Ｇｒｍを求める。
【０１２４】
次に、ステップ７０８で、図１１に示したようなＧｒｍのヒステリシスに基づいて残留磁場Ｂｒ（ｋ）を求める。
次に、ステップ７１０で、全ての測定点について残留磁場Ｂｒ（ｋ）を求めたかどうかを判定し、否の場合はステップ７１２で測定点を変更し、ステップ７０６，７０８の動作を繰り返す。
【０１２５】
全ての測定点について残留磁場Ｂｒ（ｋ）を求めた終わったときは、ステップ７１４で、Ｂｒ（ｋ）を球面調和関数にフィッティングする。この結果、ステップ７１６で、残留磁場の空間成分
【０１２６】
【数２６】
【０１２７】
が求まる。
これら残留磁場の空間成分は、球形領域ＳＶすなわち撮影空間における残留磁場の特性を表す項目であるから、以上の測定および関数フィッティングにより、マグネットシステム１００の残留磁場特性が解析（測定）できたことになる。
【０１２８】
そこで、このような解析結果を勾配コイルの設計に反映させれば、マグネットシステム１００の残留磁場特性を改善することができる。
図１３に、残留磁場解析結果を反映させた勾配コイルの設計手順の一例をフロー図によって示す。ここではＺコイルの設計手順を示す。
【０１２９】
同図に示すように、ステップ９０２で、Ｚコイルの電流経路の半径ｒ１，ｒ２，・・・，ｒＭの初期値を設定する。
次に、ステップ９０４で、上記の初期値を持つＺコイルが撮影領域の球面上に生じる磁界（磁場強度）
【０１３０】
【数２７】
【０１３１】
ただし、
【０１３２】
【数２８】
【０１３３】
を計算する。磁界の計算はビオ・サバール（Ｂｉｏｔ−Ｓａｖａｒｔ）の法則を用いて行う。
次に、ステップ９０６で、上記の磁界を球面調和関数にフィッティングして
【０１３４】
【数２９】
【０１３５】
項を求める。
次に、ステップ９０８で、前述の残留磁場解析により特に
【０１３６】
【数３０】
【０１３７】
項が大きいことが判明しているとき、
【０１３８】
【数３１】
【０１３９】
となるように
【０１４０】
【数３２】
【０１４１】
を調節する。ただし、δは許容値である。
次に、ステップ９１０で、上記の条件を満足し、かつ、そのように調節したＺコイルが発生する磁界が所定の直線性を満足するか否かを判定する。
【０１４２】
満足しない場合は、ステップ９１２で、非線形計画法により新たに
【０１４３】
【数３３】
【０１４４】
を設定し、ステップ９０４以降の処理を繰り返す。
満足する場合は、ｒ１，ｒ２，・・・，ｒＭの値を確定する。そしてそのような半径の電気経路を持つＺコイルを製造する。そのようなＺコイルを備えたマグネットシステム１００は残留磁場特性が改善されたものとなる。したがって、そのようなマグネットシステムを用いる磁気共鳴撮影装置は、残留磁場の影響が少ない品質の良い画像を撮影することができる。
【０１４５】
【発明の効果】
以上詳細に説明したように、本発明によれば、勾配磁場の印加停止後に残留する磁場を測定する方法、測定によって判明した残留磁場の特性を反映させて勾配コイルを製造する方法、その方法によって製造された勾配コイル、および、そのような勾配コイルを有する磁気共鳴撮影装置を実現することができる。
【図面の簡単な説明】
【図１】本発明の実施の形態の一例の装置のブロック図である。
【図２】図１に示した装置が実行するパルスシーケンスの一例を示す図である。
【図３】図１に示した装置が実行するパルスシーケンスの一例を示す図である。
【図４】図１に示した装置におけるマグネットシステムの勾配コイル部付近の構成を示す模式図である。
【図５】勾配コイルの電流経路のパターンを示す略図である。
【図６】測定点Ｐの３次元座標を示す図である。
【図７】球面上の測定点Ｐを配置する経度および緯度を示す図である。
【図８】磁場測定のためのパルスシーケンスを示す図である。
【図９】勾配磁場と残留磁場の関係を示す図である。
【図１０】磁場測定のための勾配磁場印加のシーケンスを示す図である。
【図１１】測定値のヒステリシスを示す図である。
【図１２】磁場測定および関数フィッティングの手順を示すフロー図である。
【図１３】勾配コイル設計手順を示すフロー図である。
【符号の説明】
１００ マグネットシステム
１０２ 主磁場マグネット部
１０６ 勾配コイル部
１０８ ＲＦコイル部
１３０ 勾配駆動部
１４０ ＲＦ駆動部
１５０ データ収集部
１６０ 制御部
１７０ データ処理部
１８０ 表示部
１９０ 操作部
２０２ 磁極片
２０４ Ｘコイル
２０６ Ｙコイル
２０８ Ｚコイル
３００ 対象
５００ クレードル
【発明の属する技術分野】
本発明は、磁場測定方法、勾配コイル（ｃｏｉｌ）製造方法、勾配コイルおよび磁気共鳴撮影装置に関し、特に、勾配磁場の印加停止後に残留する磁場を測定する方法、残留磁場の高次項が小さい勾配コイルとその製造方法、並びに、そのような勾配コイルを有する磁気共鳴撮影装置に関する。
【０００２】
【従来の技術】
磁気共鳴撮影（ＭＲＩ：Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）装置では、マグネットシステム（ｍａｇｎｅｔ ｓｙｓｔｅｍ）の内部空間すなわち静磁場を形成した空間に撮影する対象を搬入し、勾配磁場および高周波磁場を印加して対象内に磁気共鳴信号を発生させ、その受信信号に基づいて断層像を生成（再構成）する。
【０００３】
静磁場を発生するのに永久磁石を用いるマグネットシステムでは、互いに対向する１対の永久磁石の先端に静磁場空間の磁束分布を均一化するための磁極片をそれぞれ設け、勾配磁場発生用の勾配コイルはこの磁極片の磁極面に沿って設けられる。
【０００４】
【発明が解決しようとする課題】
上記のようなマグネットシステムでは、勾配コイルが磁極片に近接しているので勾配磁場により磁極片が磁化され、その残留磁化が形成する磁場により、スピン（ｓｐｉｎ）の位相が、あたかも時定数が極めて長い渦電流が存在するかのような影響を受ける。このため、精密な位相制御が要求される例えばファースト・スピンエコー（ＦＳＥ：Ｆａｓｔ Ｓｐｉｎ Ｅｃｈｏ）法等による撮影に支障をきたす。
【０００５】
これまでは、残留磁化による磁場を測定する適切な方法がなかったので、その影響を除去する方策を立てることができなかった。
そこで、本発明の課題は、勾配磁場の印加停止後に残留する磁場を測定する方法、測定によって判明した残留磁場の特性を反映させて勾配コイルを製造する方法、その方法によって製造された勾配コイル、および、そのような勾配コイルを有する磁気共鳴撮影装置を実現することである。
【０００６】
【課題を解決するための手段】
（１）上記の課題を解決するための１つの観点での発明は、測定空間に想定した球の表面の測定点に磁気共鳴可能な試料を配置し、前記測定空間に勾配磁場を印加し、前記勾配磁場の印加停止後にＲＦ励起を行って前記試料が生じるＦＩＤ信号を測定し、前記ＦＩＤ信号の位相の微分値に基づいて前記測定点における磁場強度を求め、前記磁場強度を前記測定空間における磁場強度分布を表す球面調和関数にフィッティングすることを特徴とする磁場測定方法である。
【０００７】
この観点での発明では、残留磁化が生じる磁場を、球面上の測定点に配置した試料のＦＩＤ信号を利用して測定し、測定値を球面調和関数にフィッティングして、残留磁場の特性を表す球面調和関数を特定する。
【０００８】
（２）上記の課題を解決するための他の観点での発明は、測定空間に想定した球の表面の測定点に磁気共鳴可能な試料を配置し、前記球の中心に磁気共鳴可能な試料を配置し、前記測定空間に勾配磁場を印加し、前記勾配磁場の印加停止後にＲＦ励起を行って前記試料が生じるＦＩＤ信号を測定し、前記測定点および前記球の中心の試料がそれぞれ生じるＦＩＤ信号の位相の微分値の差に基づいて前記測定点における磁場強度を求め、前記磁場強度を前記測定空間における磁場強度分布を表す球面調和関数にフィッティングすることを特徴とする磁場測定方法である。
【０００９】
この観点での発明では、球の中心に配置した試料のＦＩＤ信号を測定してこれをリファレンスとする。
（３）上記の課題を解決するための他の観点での発明は、測定空間に想定した球の表面の測定点に磁気共鳴可能な試料を配置し、１つの極性における最大勾配から反対極性における最大勾配まで順次に勾配を変化させ次いで前記反対極性における最大勾配から前記１つの極性における最大勾配まで順次に勾配を変化させながら間欠的に勾配磁場を印加し、前記勾配磁場の間欠期間中にＲＦ励起を行って前記試料が生じるＦＩＤ信号を測定し、前記ＦＩＤ信号の位相の微分値を求め、前記勾配の一巡に伴う前記微分値のヒステリシスに基づいて前記測定点における磁場強度を求め、前記磁場強度を前記測定空間における磁場強度分布を表す球面調和関数にフィッティングすることを特徴とする磁場測定方法である。
【００１０】
この観点での発明では、勾配磁場を１つの極性における最大勾配と反対極性における最大勾配の間で往復変化させる。これにより、残留磁化の履歴を一定にして磁場測定を行うことができる。
【００１１】
（４）上記の課題を解決するための他の観点での発明は、測定空間に想定した球の表面の測定点に磁気共鳴可能な試料を配置し、前記球の中心に磁気共鳴可能な試料を配置し、１つの極性における最大勾配から反対極性における最大勾配まで順次に勾配を変化させ次いで前記反対極性における最大勾配から前記１つの極性における最大勾配まで順次に勾配を変化させながら間欠的に勾配磁場を印加し、前記勾配磁場の間欠期間中にＲＦ励起を行って前記試料が生じるＦＩＤ信号を測定し、前記測定点および前記球の中心の試料がそれぞれ生じるＦＩＤ信号の位相の微分値の差を求め、前記勾配の一巡に伴う前記微分値の差のヒステリシスに基づいて前記測定点における磁場強度を求め、前記磁場強度を前記測定空間における磁場強度分布を表す球面調和関数にフィッティングすることを特徴とする磁場測定方法である。
【００１２】
この観点での発明では、球の中心に配置した試料のＦＩＤ信号を測定してこれをリファレンスとする。
（５）上記の課題を解決するための他の観点での発明は、前記測定を複数の測定点について逐次に行うことを特徴とする（１）ないし（４）のうちのいずれか１つに記載の磁場測定方法である。
【００１３】
この観点での発明では、複数の測定点について逐次に測定を行うので、同一試料を使い回しすることができる。
（６）上記の課題を解決するための他の観点での発明は、勾配コイルを製造するにあたり、空間に想定した球面上に前記勾配コイルが発生する磁界を計算し、前記磁界を球面調和関数にフィッティングし、前記フィッティングした球面調和関数における高次項のうち、（１）ないし（５）のうちのいずれか１つに記載の方法で測定した磁場を表す球面調和関数における最大の高次項に相当する高次項が小さくなるように、コイルの電流経路を決定することを特徴とする勾配コイル製造方法である。
【００１４】
この観点での発明では、測定によって判明した残留磁場の最大の高次項に相当する高次項が小さくなるように、勾配コイルの電流経路を決めるので、影響が小さい残留磁化を生じる勾配コイルを製造することができる。
【００１５】
（７）上記の課題を解決するための他の観点での発明は、電流経路を流れる電流によって勾配磁場を発生する勾配コイルであって、前記電流経路が下記の手順で決定されたことを特徴とする勾配コイルである。
【００１６】
記
（イ）空間に想定した球面上に前記勾配コイルが発生する磁界を計算する。
（ロ）前記磁界を球面調和関数にフィッティングする。
【００１７】
（ハ）前記フィッティングした球面調和関数における高次項のうち、（１）ないし（５）のうちのいずれか１つに記載の方法で測定した磁場を表す球面調和関数における最大の高次項に相当する高次項が小さくなるように、コイルの電流経路を決定する。
【００１８】
この観点での発明では、測定によって判明した残留磁場の最大の高次項に相当する高次項が小さくなるように、勾配コイルの電流経路が決められているので、影響が小さい残留磁化を生じる勾配コイルとなる。
【００１９】
（８）上記の課題を解決するための他の観点での発明は、静磁場、勾配磁場および高周波磁場を用いて獲得した磁気共鳴信号に基づいて画像を構成する磁気共鳴撮影装置であって、前記勾配磁場を発生する勾配コイルとして、下記の手順で決定された電流経路を有する勾配コイルを具備することを特徴とする磁気共鳴撮影装置である。
【００２０】
記
（イ）空間に想定した球面上に前記勾配コイルが発生する磁界を計算する。
（ロ）前記磁界を球面調和関数にフィッティングする。
【００２１】
（ハ）前記フィッティングした球面調和関数における高次項のうち、（１）ないし（５）のうちのいずれか１つに記載の方法で測定した磁場を表す球面調和関数における最大の高次項に相当する高次項が小さくなるように、コイルの電流経路を決定する。
【００２２】
この観点での発明では、測定によって判明した残留磁場の最大の高次項に相当する高次項が小さくなるように勾配コイルの電流経路が決められているので、影響が小さい残留磁化を生じる勾配コイルを用いて適正な磁気共鳴撮影をすることができる。
【００２３】
【発明の実施の形態】
以下、図面を参照して本発明の実施の形態を詳細に説明する。なお、本発明は実施の形態に限定されるものではない。図１に磁気共鳴撮影装置のブロック（ｂｌｏｃｋ）図を示す。本装置は本発明の実施の形態の一例である。本装置の構成によって、本発明の装置に関する実施の形態の一例が示される。本装置の動作によって、本発明の方法に関する実施の形態の一例が示される。
【００２４】
図１に示すように、本装置はマグネットシステム１００を有する。マグネットシステム１００は、主磁場マグネット部１０２、勾配コイル部１０６およびＲＦ（ｒａｄｉｏ ｆｒｅｑｕｅｎｃｙ）コイル部１０８を有する。これら主磁場マグネット部１０２および各コイル部は、いずれも空間を挟んで互いに対向する１対のものからなる。また、いずれも概ね円板状の形状を有し中心軸を共有して配置されている。マグネットシステム１００の内部空間（ボア：ｂｏｒｅ）に、対象３００がクレードル（ｃｒａｄｌｅ）５００に搭載されて図示しない搬送手段により搬入および搬出される。
【００２５】
主磁場マグネット部１０２はマグネットシステム１００の内部空間に静磁場を形成する。静磁場の方向は概ね対象３００の体軸方向と直交する。すなわちいわゆる垂直磁場を形成する。主磁場マグネット部１０２は例えば永久磁石等を用いて構成される。なお、永久磁石に限らず超伝導電磁石あるいは常伝導電磁石等を用いて構成しても良いのはもちろんである。
【００２６】
勾配コイル部１０６は静磁場強度に勾配を持たせるための勾配磁場を生じる。発生する勾配磁場は、スライス（ｓｌｉｃｅ）勾配磁場、リードアウト（ｒｅａｄ ｏｕｔ）勾配磁場およびフェーズエンコード（ｐｈａｓｅ ｅｎｃｏｄｅ）勾配磁場の３種であり、これら３種類の勾配磁場に対応して勾配コイル部１０６は図示しない３系統の勾配コイルを有する。
【００２７】
３系統の勾配コイルは、互いに直交する３方向において静磁場にそれぞれ勾配を付与するための３つの勾配磁場をそれぞれ発生する。３方向のうちの１つは静磁場の方向（垂直方向）であり、通常これをｚ方向とする。他の１つは水平方向であり、通常これをｙ方向とする。残りの１つはｚ，ｙ方向に垂直な方向であり、通常これをｘ方向とする。x方向は垂直面内でｚ方向に垂直であり、水平面内でｙ方向に垂直である。以下、ｘ，ｙ，ｚを勾配軸ともいう。
【００２８】
ｘ，ｙ，ｚはいずれもスライス勾配の軸とすることができる。いずれか１つをスライス勾配の軸としたとき、残り２つのうちの１つをフェーズエンコード勾配の軸とし、他をリードアウト勾配の軸とする。３系統の勾配コイルについては後にあらためて説明する。
【００２９】
ＲＦコイル部１０８は静磁場空間に対象３００の体内のスピンを励起するためのＲＦ励起信号を送信する。ＲＦコイル部１０８は、また、励起されたスピンが生じる磁気共鳴信号を受信する。ＲＦコイル部１０８は図示しない送信用のコイルおよび受信用のコイルを有する。送信用のコイルおよび受信用のコイルは、同じコイルを兼用するかあるいはそれぞれ専用のコイルを用いる。
【００３０】
勾配コイル部１０６には勾配駆動部１３０が接続されている。勾配駆動部１３０は勾配コイル部１０６に駆動信号を与えて勾配磁場を発生させる。勾配駆動部１３０は、勾配コイル部１０６における３系統の勾配コイルに対応して、図示しない３系統の駆動回路を有する。
【００３１】
ＲＦコイル部１０８にはＲＦ駆動部１４０が接続されている。ＲＦ駆動部１４０はＲＦコイル部１０８に駆動信号を与えてＲＦ励起信号を送信し、対象３００の体内のスピンを励起する。
【００３２】
ＲＦコイル部１０８には、また、データ（ｄａｔａ）収集部１５０が接続されている。データ収集部１５０はＲＦコイル部１０８が受信した受信信号を取り込み、それをビューデータ（ｖｉｅｗ ｄａｔａ）として収集する。
【００３３】
勾配駆動部１３０、ＲＦ駆動部１４０およびデータ収集部１５０には制御部１６０が接続されている。制御部１６０は、勾配駆動部１３０ないしデータ収集部１５０をそれぞれ制御して撮影を遂行する。
【００３４】
データ収集部１５０の出力側はデータ処理部１７０に接続されている。データ処理部１７０は、例えばコンピュータ（ｃｏｍｐｕｔｅｒ）等を用いて構成される。データ処理部１７０は図示しないメモリ（ｍｅｍｏｒｙ）を有する。メモリはデータ処理部１７０用のプログラムおよび各種のデータを記憶している。本装置の機能は、データ処理部１７０がメモリに記憶されたプログラムを実行することによりを実現される。
【００３５】
データ処理部１７０は、データ収集部１５０から取り込んだデータをメモリに記憶する。メモリ内にはデータ空間が形成される。データ空間は２次元フ−リエ（Ｆｏｕｒｉｅｒ）空間を構成する。データ処理部１７０は、これら２次元フ−リエ空間のデータを２次元逆フ−リエ変換して対象３００の画像を生成（再構成）する。２次元フ−リエ空間をｋスペース（ｋ−ｓｐａｃｅ）ともいう。
【００３６】
データ処理部１７０は制御部１６０に接続されている。データ処理部１７０は制御部１６０の上位にあってそれを統括する。データ処理部１７０には、また、表示部１８０および操作部１９０が接続されている。表示部１８０は、グラフィックディスプレー（ｇｒａｐｈｉｃ ｄｉｓｐｌａｙ）等で構成される。操作部１９０はポインティングデバイス（ｐｏｉｎｔｉｎｇ ｄｅｖｉｃｅ）を備えたキーボード（ｋｅｙｂｏａｒｄ）等で構成される。
【００３７】
表示部１８０は、データ処理部１７０から出力される再構成画像および各種の情報を表示する。操作部１９０は、操作者によって操作され、各種の指令や情報等をデータ処理部１７０に入力する。操作者は表示部１８０および操作部１９０を通じてインタラクティブ（ｉｎｔｅｒａｃｔｉｖｅ）に本装置を操作する。
【００３８】
図２に、本装置で撮影を行うときのパルスシーケンス（ｐｕｌｓｅ ｓｅｑｕｅｎｃｅ）の一例を示す。このパルスシーケンスは、グラディエントエコー（ＧＲＥ：Ｇｒａｄｉｅｎｔ Ｅｃｈｏ）法のパルスシーケンスである。
【００３９】
すなわち、（１）はＧＲＥ法におけるＲＦ励起用のα°パルスのシーケンスであり、（２）、（３）、（４）および（５）は、同じくそれぞれ、スライス勾配Ｇｓ、リードアウト勾配Ｇｒ、フェーズエンコード勾配ＧｐおよびグラディエントエコーＭＲのシーケンスである。なお、α°パルスは中心信号で代表する。パルスシーケンスは時間軸ｔに沿って左から右に進行する。
【００４０】
同図に示すように、α°パルスによりスピンのα°励起が行われる。フリップアングル（ｆｌｉｐ ａｎｇｌｅ）α°は９０°以下である。このときスライス勾配Ｇｓが印加され所定のスライスについての選択励起が行われる。
【００４１】
α°励起後、フェーズエンコード勾配Ｇｐによりスピンのフェーズエンコードが行われる。次に、リードアウト勾配Ｇｒにより先ずスピンをディフェーズ（ｄｅｐｈａｓｅ）し、次いでスピンをリフェーズ（ｒｅｐｈａｓｅ）して、グラディエントエコーＭＲを発生させる。グラディエントエコーＭＲの信号強度は、α°励起からエコータイム（ｅｃｈｏ ｔｉｍｅ）ＴＥ後の時点で最大となる。グラディエントエコーＭＲはデータ収集部１５０によりビューデータとして収集される。
【００４２】
このようなパルスシーケンスが周期ＴＲ（ｒｅｐｅｔｉｔｉｏｎ ｔｉｍｅ）で６４〜５１２回繰り返される。繰り返しのたびにフェーズエンコード勾配Ｇｐを変更し、毎回異なるフェーズエンコードを行う。これによって、ｋスペースを埋める６４〜５１２ビューのビューデータが得られる。
【００４３】
磁気共鳴撮影用パルスシーケンスの他の例を図３に示す。このパルスシーケンスは、スピンエコー（ＳＥ：Ｓｐｉｎ Ｅｃｈｏ）法のパルスシーケンスである。
【００４４】
すなわち、（１）はＳＥ法におけるＲＦ励起用の９０°パルスおよび１８０°パルスのシーケンスであり、（２）、（３）、（４）および（５）は、同じくそれぞれ、スライス勾配Ｇｓ、リードアウト勾配Ｇｒ、フェーズエンコード勾配ＧｐおよびスピンエコーＭＲのシーケンスである。なお、９０°パルスおよび１８０°パルスはそれぞれ中心信号で代表する。パルスシーケンスは時間軸ｔに沿って左から右に進行する。
【００４５】
同図に示すように、９０°パルスによりスピンの９０°励起が行われる。このときスライス勾配Ｇｓが印加され所定のスライスについての選択励起が行われる。９０°励起から所定の時間後に、１８０°パルスによる１８０°励起すなわちスピン反転が行われる。このときもスライス勾配Ｇｓが印加され、同じスライスについての選択的反転が行われる。
【００４６】
９０°励起とスピン反転の間の期間に、リードアウト勾配Ｇｒおよびフェーズエンコード勾配Ｇｐが印加される。リードアウト勾配Ｇｒによりスピンのディフェーズが行われる。フェーズエンコード勾配Ｇｐによりスピンのフェーズエンコードが行われる。
【００４７】
スピン反転後、リードアウト勾配ＧｒでスピンをリフェーズしてスピンエコーＭＲを発生させる。スピンエコーＭＲの信号強度は、９０°励起からＴＥ後の時点で最大となる。スピンエコーＭＲはデータ収集部１５０によりビューデータとして収集される。このようなパルスシーケンスが周期ＴＲで６４〜５１２回繰り返される。繰り返しのたびにフェーズエンコード勾配Ｇｐを変更し、毎回異なるフェーズエンコードを行う。これによって、ｋスペースを埋める６４〜５１２ビューのビューデータが得られる。
【００４８】
なお、撮影に用いるパルスシーケンスはＧＲＥ法またはＳＥ法に限るものではなく、例えば、ＦＳＥ（Ｆａｓｔ Ｓｐｉｎ Ｅｃｈｏ）法、ファーストリカバリＦＳＥ（Ｆａｓｔ Ｒｅｃｏｖｅｒｙ Ｆａｓｔ Ｓｐｉｎ Ｅｃｈｏ）法、エコープラナー・イメージング（ＥＰＩ：Ｅｃｈｏ Ｐｌａｎａｒ Ｉｍａｇｉｎｇ）等、他の適宜の技法のものであって良い。
【００４９】
データ処理部１７０は、ｋスペースのビューデータを２次元逆フ−リエ変換して対象３００の断層像を再構成する。再構成した画像はメモリに記憶し、また、表示部１８０で表示する。
【００５０】
図４に、勾配コイル部１０６付近のマグネットシステム１００の構造を断面図により模式的に示す。同図において、Ｏは静磁場の中心すなわちマグネットセンター（ｍａｇｎｅｔ ｃｅｎｔｅｒ）であり、ｘ，ｙ，ｚは前述した３方向である。
【００５１】
マグネットセンターＯを中心とする半径Ｒの球形領域ＳＶ（ｓｐｈｅｒｉｃ ｖｏｌｕｍｅ）が撮影領域であり、マグネットシステム１００はこのＳＶにおいて静磁場および勾配磁場が所定の精度を持つように構成される。
【００５２】
１対の主磁場マグネット部１０２は互いに対向する１対の磁極片（ポールピース：ｐｏｌｅ ｐｉｅｃｅ）２０２を有する。磁極片２０２は例えば軟鉄等の高透磁率の磁性材料で構成され、静磁場空間における磁束分布を均一化する働きをする。
【００５３】
磁極片２０２は概ね円板形状を成すが、周縁部が板面に垂直な方向（ｚ方向）、すなわち、磁極片２０２同士が互いに対向する方向に突出している。突出部は磁極片２０２の周縁部における磁束密度の低下を補う働きをする。
【００５４】
突出した周縁部の内側に形成される磁極片２０２の凹部に、勾配コイル部１０６が設けられる。勾配コイル部１０６は、Ｘコイル２０４、Ｙコイル２０６およびＺ２０８コイルを有する。
【００５５】
Ｘコイル２０４、Ｙコイル２０６およびＺコイル２０８は、それぞれ、本発明の勾配コイルの実施の形態の一例である。各コイルはいずれも概ね円板形状を成し、磁極片２０２の磁極面に、図示しない適宜の取り付け手段により、順次に層を成すように取り付けられている。
【００５６】
図５に、Ｘコイル２０４、Ｙコイル２０６およびＺ２０８コイルの電流経路のパターン（ｐａｔｔｅｒｎ）を略図によって示す。同図に示すように、Ｘコイル２０４は、円の中心に近い部分ではｙ方向に平行な直線状の複数の主電流経路（メインパス：ｍａｉｎ ｐａｓｓ）を有する。これらメインパスは円の中心を通るｙ軸に関して対称である。メインパスの帰路（リターンパス：ｒｅｔｕｒｎ ｐａｓｓ）は円周に沿って形成される。最外側のリターンパスの半径すなわちＸコイル２０４の外径はｒ００である。
【００５７】
Ｙコイル２０６は、円の中心に近い部分ではｘ方向に平行な直線状の複数のメインパスを有する。これらメインパスは円の中心を通るｘ軸に関して対称である。メインパスのリターンパスは円周に沿って形成される。最外側のリターンパスの半径すなわちＹコイル２０６の外径はｒ００である。
【００５８】
Ｚコイル２０８は同心円を成す複数の電流経路を有する。これら電流経路は全てがメインパスである。各メインパスの半径は内側から順にｒ１，ｒ２，・・・，ｒＭである。ｒＭはＺコイル２０８の外径となる。
【００５９】
このようなマグネットシステム１００について、撮影領域すなわち球形領域ＳＶにおける残留磁場を測定する方法を説明する。先ず、球形領域ＳＶの表面に複数の測定点を設定する。測定点Ｐを、図６に示すように、３次元の極座標で表すとした場合、複数の測定点は、半径ｒを一定とし、角度θ，φの組み合わせを複数通りに設定することにより行われる。角度θはいわば緯度に相当し、角度φはいわば経度に相当する。
【００６０】
角度θは、図７（ａ）に示すように、例えば２２．５度刻みで設定する。これによって７つの緯度が設定される。角度φは、同図の（ｂ）に示すように、例えば３０度刻みで設定する。これによって１２の経度が設定される。
【００６１】
このような緯度と経度の交点を測定点Ｐとする。したがって、球面上の測定点数は８４となる。これらの測定点に加えて、マグネットセンターＯにも測定点Ｐ０を設定する。測定点Ｐ０はリファレンス（ｒｅｆｅｒｅｎｃｅ）用である。リファレンスが不要な場合は測定点Ｐ０は省略可能である。
【００６２】
測定点には、磁気共鳴信号を生じる試料に信号検出用のＲＦコイルを取り付けたプローブ（ｐｒｏｂｅ）が設置される。試料としては、例えば、硫酸銅（ＣｕＳＯ４）や塩化ニッケル（ＮｉＣｌ）等の水溶液を微小な容器に封入したものが用いられる。試料の体積は、撮影空間の１ボクセル（ｖｏｘｅｌ）に相当する大きさとする。
【００６３】
そのようなプローブを適宜の支持手段で支持して上記の測定点に順次に位置決めし、球面上の各測定点およびマグネットセンターＯにおける磁場を測定する。プローブは少なくとも２つ用い、球面上の測定点とマグネットセンターＯの磁場を同時に測定する。リファレンスを必要としないときは、球面上の測定点だけで測定する。
【００６４】
図８に、磁場測定用のパルスシーケンスを示す。同図において（１）は勾配磁場のシーケンス、（２）はＲＦ励起のシーケンス、（３）は磁気共鳴信号すなわちＦＩＤ（Ｆｒｅｅ Ｉｎｄｕｃｔｉｏｎ Ｄｅｃａｙ）信号のシーケンスである。
【００６５】
同図に示すように、１つの軸に勾配磁場Ｇを印加し、それが終了した後にＲＦ励起を行う。ＲＦ励起に基づきプローブ中の試料がＦＩＤ信号を生じる。このＦＩＤ信号を、プローブに備わるＲＦコイルで受信し、受信データをデータ収集部１５０で収集する。
【００６６】
勾配磁場Ｇが例えば図９に実線で示すようなものであるとしたとき、それによって磁極片２０２が磁化されることにより、勾配磁場の印加をやめた後にも、例えば一点鎖線で示すような残留磁場が存在する。このため、図８に斜線で示すように、あたかも勾配磁場が残留するかのような状態になる。
【００６７】
このような残留磁場の下で測定されるＦＩＤ信号は次式で表される。
【００６８】
【数１】
ここで、
ρ：スピン密度
γ：磁界回転比
Ｇｒｍ：磁場強度
（１）式に示すように、磁場強度ＧｒｍはＦＩＤ信号の位相に影響を与える。そこで、複素数で表現したＦＩＤ信号
【００７０】
【数２】
から、
【００７２】
【数３】
により位相を求め、その時間微分を用いて
【００７４】
【数４】
により磁場強度を求める。これによって、測定点Ｐにおける磁場の測定値が得られる。
実際には、測定空間に静磁場が存在するので、マグネットセンターＯで測定したＦＩＤ信号の位相φ０（ｔ）の微分値を求め、
【００７６】
【数５】
により磁場強度を求める。これにより、静磁場の強度変動による測定誤差を除去することができる。以上の計算は、データ処理部１７０によって行う。
勾配磁場の印加の経過を反映して、磁極片２０２では磁化の履歴が生じている。磁化の履歴はマイナー・ヒステリシスループ（ｍｉｎｏｒ ｈｙｓｔｅｒｅｓｉｓ ｌｏｏｐ）を描く。このような履歴を持つものについて常に同一の条件で測定を行うために、磁場測定は次のようにして行われる。
【００７８】
図１０に、測定時の１つの測定点あたりの勾配磁場印加のシーケンスを示す。同図に示すように、勾配磁場は、その強度を順次変化させながら間欠的に印加する。すなわち、例えば、シーケンスの前半は、本装置で発生可能な最大の負の勾配−Ｇｍａｘから最大の正の勾配＋Ｇｍａｘまで順次に変化させながら印加し、後半は、最大の正の勾配＋Ｇｍａｘから最大の負の勾配−Ｇｍａｘまで順次に変化させながら印加する。このように勾配の強度を一巡させることにより、シーケンスを開始する前と後における磁極片２０２の磁化状態を、同一にすることができる。
【００７９】
勾配印加の間欠期間中は斜線で概念的に示すように勾配磁場が残留する。この期間中に、前述のようなＲＦ励起とＦＩＤ測定を毎回行い、各回の測定信号に基づいて、上記のようにしてそれぞれ磁場強度Ｇｒｍを求める。これによって、１つの測定点あたり例えば１５個の磁場測定値が得られる。
【００８０】
このような磁場測定値について、印加した勾配との対応を求めると、例えば図１１に示すようなヒステリシスループが得られる。そこで、このヒステリシスループのＧｒｍ軸方向の幅、すなわち、ヒステリシスループが横切るＧｒｍ軸上の２箇所の値の差を測定点Ｐの磁場強度Ｂｒ（ｋ）とする。
【００８１】
このようにして求めた８４個の測定点Ｐの磁場測定値を球面調和関数にフィッティング（ｆｉｔｔｉｎｇ）して、残留磁場の空間成分を求める。フィッティングはデータ処理部１７０により行う。以下、磁場測定値を単に測定値ともいう。
【００８２】
球面調和関数展開式は次式で与えられる。
【００８３】
【数６】
ここで、
【００８５】
【数７】
はルジャンドル（Ｌｅｇｅｎｄｒｅ）の倍関数、
【００８７】
【数８】
は定数である。
フィッティングは、例えば次のようにして行われる。
１． ｍ＝０の場合
この場合、（６）式は次のようになる。
【００８９】
【数９】
ｚ座標が異なる測定点での測定値を
【００９１】
【数１０】
とすると次式が成り立つ。
【００９３】
【数１１】
この方程式を解くことにより、
【００９５】
【数１２】
成分が得られる。これらは通常
【００９７】
【数１３】
成分と呼ばれる。
２． ｍ≠０の場合
例えばｍ＝１のとき（６）式は次式のようになる。
【００９９】
【数１４】
球面上のあるｘｙ平面内の測定値を
【０１０１】
【数１５】
とすると、円周上の測定点における測定値は、
【０１０３】
【数１６】
このとき、
【０１０５】
【数１７】
【数１８】
各ｘｙ平面について
【０１０８】
【数１９】
を求めれば次式が成り立つ。
【０１１０】
【数２０】
この方程式を解けば、
【０１１２】
【数２１】
が求まる。これらは通常
【０１１４】
【数２２】
と呼ばれる。
３．ｍ＝２の場合
この場合は（１０）式において
【０１１６】
【数２３】
とすれば、同様に
【０１１８】
【数２４】
すなわち
【０１２０】
【数２５】
を求めることができる。
図１２に、以上のような複数の測定点における磁場測定とそれに基づく球面調和関数のフィッティングのフロー（ｆｌｏｗ）図を示す。
【０１２２】
同図に示すように、ステップ（ｓｔｅｐ）７０２で、球面上の測定点を設定する。測定点は前述の要領で設定される。測定点はＬ個である。
次に、ステップ７０４で、球の中心すなわちマグネットセンターＯと球面上の１つの測定点にそれぞれプローブをセット（ｓｅｔ）する。
【０１２３】
次に、ステップ７０６で、図１０に示したような測定シーケンスを実行して、２つのプローブのＦＩＤ信号を測定し、それらの位相の微分値の差に基づいて磁場測定値Ｇｒｍを求める。
【０１２４】
次に、ステップ７０８で、図１１に示したようなＧｒｍのヒステリシスに基づいて残留磁場Ｂｒ（ｋ）を求める。
次に、ステップ７１０で、全ての測定点について残留磁場Ｂｒ（ｋ）を求めたかどうかを判定し、否の場合はステップ７１２で測定点を変更し、ステップ７０６，７０８の動作を繰り返す。
【０１２５】
全ての測定点について残留磁場Ｂｒ（ｋ）を求めた終わったときは、ステップ７１４で、Ｂｒ（ｋ）を球面調和関数にフィッティングする。この結果、ステップ７１６で、残留磁場の空間成分
【０１２６】
【数２６】
が求まる。
これら残留磁場の空間成分は、球形領域ＳＶすなわち撮影空間における残留磁場の特性を表す項目であるから、以上の測定および関数フィッティングにより、マグネットシステム１００の残留磁場特性が解析（測定）できたことになる。
【０１２８】
そこで、このような解析結果を勾配コイルの設計に反映させれば、マグネットシステム１００の残留磁場特性を改善することができる。
図１３に、残留磁場解析結果を反映させた勾配コイルの設計手順の一例をフロー図によって示す。ここではＺコイルの設計手順を示す。
【０１２９】
同図に示すように、ステップ９０２で、Ｚコイルの電流経路の半径ｒ１，ｒ２，・・・，ｒＭの初期値を設定する。
次に、ステップ９０４で、上記の初期値を持つＺコイルが撮影領域の球面上に生じる磁界（磁場強度）
【０１３０】
【数２７】
ただし、
【０１３２】
【数２８】
を計算する。磁界の計算はビオ・サバール（Ｂｉｏｔ−Ｓａｖａｒｔ）の法則を用いて行う。
次に、ステップ９０６で、上記の磁界を球面調和関数にフィッティングして
【０１３４】
【数２９】
項を求める。
次に、ステップ９０８で、前述の残留磁場解析により特に
【０１３６】
【数３０】
項が大きいことが判明しているとき、
【０１３８】
【数３１】
となるように
【０１４０】
【数３２】
を調節する。ただし、δは許容値である。
次に、ステップ９１０で、上記の条件を満足し、かつ、そのように調節したＺコイルが発生する磁界が所定の直線性を満足するか否かを判定する。
【０１４２】
満足しない場合は、ステップ９１２で、非線形計画法により新たに
【０１４３】
【数３３】
を設定し、ステップ９０４以降の処理を繰り返す。
満足する場合は、ｒ１，ｒ２，・・・，ｒＭの値を確定する。そしてそのような半径の電気経路を持つＺコイルを製造する。そのようなＺコイルを備えたマグネットシステム１００は残留磁場特性が改善されたものとなる。したがって、そのようなマグネットシステムを用いる磁気共鳴撮影装置は、残留磁場の影響が少ない品質の良い画像を撮影することができる。
【０１４５】
【発明の効果】
以上詳細に説明したように、本発明によれば、勾配磁場の印加停止後に残留する磁場を測定する方法、測定によって判明した残留磁場の特性を反映させて勾配コイルを製造する方法、その方法によって製造された勾配コイル、および、そのような勾配コイルを有する磁気共鳴撮影装置を実現することができる。
【図面の簡単な説明】
【図１】本発明の実施の形態の一例の装置のブロック図である。
【図２】図１に示した装置が実行するパルスシーケンスの一例を示す図である。
【図３】図１に示した装置が実行するパルスシーケンスの一例を示す図である。
【図４】図１に示した装置におけるマグネットシステムの勾配コイル部付近の構成を示す模式図である。
【図５】勾配コイルの電流経路のパターンを示す略図である。
【図６】測定点Ｐの３次元座標を示す図である。
【図７】球面上の測定点Ｐを配置する経度および緯度を示す図である。
【図８】磁場測定のためのパルスシーケンスを示す図である。
【図９】勾配磁場と残留磁場の関係を示す図である。
【図１０】磁場測定のための勾配磁場印加のシーケンスを示す図である。
【図１１】測定値のヒステリシスを示す図である。
【図１２】磁場測定および関数フィッティングの手順を示すフロー図である。
【図１３】勾配コイル設計手順を示すフロー図である。
【符号の説明】
１００ マグネットシステム
１０２ 主磁場マグネット部
１０６ 勾配コイル部
１０８ ＲＦコイル部
１３０ 勾配駆動部
１４０ ＲＦ駆動部
１５０ データ収集部
１６０ 制御部
１７０ データ処理部
１８０ 表示部
１９０ 操作部
２０２ 磁極片
２０４ Ｘコイル
２０６ Ｙコイル
２０８ Ｚコイル
３００ 対象
５００ クレードル
Claims (7)
- 静磁場を発生する１対の永久磁石が測定空間を介して対向して設置され、前記１対の永久磁石における前記測定空間側に１対の磁極片が設置され、前記１対の磁極片における前記測定空間側に１対の勾配コイルが設置されたマグネットシステムにおいて、
前記測定空間に想定した球の表面の測定点に磁気共鳴可能な試料を配置し、
１つの極性における最大勾配から反対極性における最大勾配まで順次に勾配を変化させ次いで前記反対極性における最大勾配から前記１つの極性における最大勾配まで順次に勾配を変化させながら間欠的に勾配を印加し、
前記勾配磁場の間欠期間中にＲＦ励起を行って前記試料が生じるＦＩＤ信号を測定し、
前記ＦＩＤ信号の位相の微分値を求め、
前記勾配の一巡に伴う前記微分値のヒステリシスに基づいて前記測定点における磁場強度を求め、
前記磁場強度を前記測定空間における磁場強度分布を表す球面調和関数にフィッティングすることを特徴とする磁場測定方法。 - 静磁場を発生する１対の永久磁石が測定空間を介して対向して設置され、前記１対の永久磁石における前記測定空間側に１対の磁極片が設置され、前記１対の磁極片における前記測定空間側に１対の勾配コイルが設置されたマグネットシステムにおいて、
前記測定空間に想定した球の表面の測定点に磁気共鳴可能な試料を配置し、
１つの極性における最大勾配から反対極性における最大勾配まで順次に勾配を変化させ次いで前記反対極性における最大勾配から前記１つの極性における最大勾配まで順次に勾配を変化させながら間欠的に勾配を印加し、
前記勾配磁場の間欠期間中にＲＦ励起を行って前記試料が生じるＦＩＤ信号を測定し、
前記測定点および前記球の中心の試料がそれぞれ生じるＦＩＤ信号の位相の微分値の差を求め、
前記勾配の一巡に伴う前記微分値のヒステリシスに基づいて前記測定点における磁場強度を求め、
前記磁場強度を前記測定空間における磁場強度分布を表す球面調和関数にフィッティングすることを特徴とする磁場測定方法。 - 前記測定を複数の測定点について逐次行うこと特徴とする請求項１又は請求項２に記載の磁場測定方法。
- 前記ヒステリシスに基づく前記測定点における磁場強度として、ヒステリシスループがそのヒステリシスカーブがプロットされたグラフの縦軸上を横切った２箇所の値の差を用いたこと特徴とする請求項１ないし請求項３のうちのいずれか１つに記載の磁場測定方法。
- 勾配コイルを製造するにあたり、
空間に想定した球面上に前記勾配コイルが発生する磁界を計算し、
前記磁界を球面調和関数にフィッティングし、
前記フィッティングした球面調和関数における高次項のうち、請求項１ないし請求項４のうちのいずれか１つに記載の方法で測定した磁場を表す球面調和関数における最大の高次項に相当する高次項が小さくなるように、コイルの電流経路を決定すること特徴とする勾配コイル製造方法。 - 電流経路を流れる電流によって勾配磁場を発生する勾配コイルであって、
前記電流経路が下記の手順で決定されたことを特徴とする勾配コイル。
記
（イ）空間に想定した球面上に前記勾配コイルが発生する磁界を計算する。
（ロ）前記磁界を球面調和関数にフィッティングする。
（ハ）前記フィッティングした球面調和関数における高次項のうち、請求項１ないし請求 項４のうちのいずれか１つに記載の方法で測定した磁場を表す球面調和関数における最大の高次項に相当する高次項が小さくなるようにコイルの電流経路を決定する。 - 静磁場、勾配磁場および高周波磁場を用いて獲得した磁気共鳴信号に基づいて画像を再構成する磁気共鳴撮影装置であって、
前記勾配磁場を発生する勾配コイルとして、下記の手順で決定された電流経路を有する勾配コイルを具備することを特徴とする磁気共鳴撮影装置。
記
（イ）空間に想定した球面上に前記勾配コイルが発生する磁界を計算する。
（ロ）前記磁界を球面調和関数にフィッティングする。
（ハ）前記フィッティングした球面調和関数における高次項のうち、請求項１ないし請求項４のうちのいずれか１つに記載の方法で測定した磁場を表す球面調和関数における最大の高次項に相当する高次項が小さくなるようにコイルの電流経路を決定する。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000117880A JP3701540B2 (ja) | 2000-04-19 | 2000-04-19 | 磁場測定方法、勾配コイル製造方法、勾配コイルおよび磁気共鳴撮影装置 |
| US09/827,037 US6813513B2 (en) | 2000-04-19 | 2001-04-05 | Method for measurement of magnetic fields, method for production of gradient coils, gradient coil, and apparatus for magnetic resonance imaging |
| EP01303473A EP1148342A3 (en) | 2000-04-19 | 2001-04-12 | Method for measurement of magnetic fields, method for production of gradient coils, gradient coil, and apparatus for magnetic resonance imaging |
| CN01119654A CN1336558A (zh) | 2000-04-19 | 2001-04-18 | 测磁场的方法、产生梯度线圈方法、梯度线圈和磁共振成像设备 |
| KR10-2001-0020699A KR100458771B1 (ko) | 2000-04-19 | 2001-04-18 | 자계 측정 방법, 그레디언트 코일 제조 방법, 및 자기 공명 촬상 장치 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000117880A JP3701540B2 (ja) | 2000-04-19 | 2000-04-19 | 磁場測定方法、勾配コイル製造方法、勾配コイルおよび磁気共鳴撮影装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001299721A JP2001299721A (ja) | 2001-10-30 |
| JP3701540B2 true JP3701540B2 (ja) | 2005-09-28 |
Family
ID=18629122
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000117880A Expired - Fee Related JP3701540B2 (ja) | 2000-04-19 | 2000-04-19 | 磁場測定方法、勾配コイル製造方法、勾配コイルおよび磁気共鳴撮影装置 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US6813513B2 (ja) |
| EP (1) | EP1148342A3 (ja) |
| JP (1) | JP3701540B2 (ja) |
| KR (1) | KR100458771B1 (ja) |
| CN (1) | CN1336558A (ja) |
Families Citing this family (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7603156B2 (en) * | 2003-07-02 | 2009-10-13 | Ge Medical Systems Global Technology Co., Llc | Systems and methods for phase encode placement |
| US8016541B2 (en) * | 2003-09-10 | 2011-09-13 | Brooks Automation, Inc. | Substrate handling system for aligning and orienting substrates during a transfer operation |
| US9691651B2 (en) | 2005-01-28 | 2017-06-27 | Brooks Automation, Inc. | Substrate handling system for aligning and orienting substrates during a transfer operation |
| GB2408345B (en) * | 2003-11-18 | 2006-09-13 | Univ Queensland | Bi-planar coil assemblies for producing specified magnetic fields |
| CN101109720B (zh) * | 2006-07-19 | 2011-02-16 | 西门子（中国）有限公司 | 测量磁性材料磁感应强度对温度变化特性的方法及装置 |
| CN104181480B (zh) | 2013-05-21 | 2017-02-08 | 上海联影医疗科技有限公司 | 磁共振装置中成像磁场测量和校正的方法及系统 |
| US10156625B2 (en) * | 2013-08-12 | 2018-12-18 | Koninklijke Philips N.V. | MR imaging with B1 mapping |
| DE102014213413B4 (de) * | 2014-07-10 | 2018-12-20 | Siemens Healthcare Gmbh | Dynamische Felderfassung in einem MRT |
| CN104730476A (zh) * | 2015-02-28 | 2015-06-24 | 中国科学院电工研究所 | 一种磁共振显微成像用平面梯度线圈 |
| CN104833930B (zh) * | 2015-04-21 | 2017-11-17 | 中国科学院电工研究所 | 开放式磁共振系统梯度线圈磁场强度的计算方法 |
| CN107045112A (zh) * | 2016-02-09 | 2017-08-15 | 温伯格医学物理有限公司 | 操作用于磁共振成像和图像引导治疗的电永磁铁的方法和设备 |
| CN105572611B (zh) * | 2016-03-04 | 2018-09-28 | 中国海洋石油集团有限公司 | 一种静磁场核磁效应分析系统 |
| CN105806928B (zh) * | 2016-03-04 | 2019-02-26 | 中国海洋石油集团有限公司 | 一种静磁场核磁效应分析方法 |
| DE102016204863B3 (de) * | 2016-03-23 | 2017-06-14 | Siemens Healthcare Gmbh | Verfahren zu einer Bestimmung einer Abweichung einer Homogenität eines Magnetfeldes eines Magnetresonanzgerätes |
| CN109598004B (zh) * | 2017-09-30 | 2022-09-20 | 中国科学院长春光学精密机械与物理研究所 | 用于微尺度磁共振成像系统的横向梯度线圈及其设计方法 |
| CN107942272B (zh) * | 2018-01-11 | 2020-08-21 | 郑州航空工业管理学院 | 主磁场不均匀下的分数域磁共振成像方法 |
| AU2019417058A1 (en) * | 2018-12-28 | 2021-07-22 | Hyperfine Operations, Inc. | Correcting for hysteresis in magnetic resonance imaging |
| US11698430B2 (en) | 2019-08-15 | 2023-07-11 | Hyperfine Operations, Inc. | Eddy current mitigation systems and methods |
| CN112741612B (zh) * | 2019-10-31 | 2023-04-25 | 上海联影医疗科技股份有限公司 | 磁共振成像方法、装置、存储介质及计算机设备 |
| CN113933769B (zh) * | 2021-10-18 | 2024-08-09 | 上海复旦数字医疗科技股份有限公司 | 一种磁共振成像系统中剩磁的测量方法 |
| CN116359270B (zh) * | 2021-12-27 | 2026-04-21 | 上海复旦数字医疗科技股份有限公司 | 一种利用双回波测量磁共振成像系统中剩磁的方法 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2588997B1 (fr) * | 1985-10-18 | 1987-11-20 | Thomson Cgr | Procede de realisation d'une bobine de gradient et bobine obtenue par ce procede |
| DE3616078A1 (de) * | 1986-05-13 | 1987-11-19 | Bruker Analytische Messtechnik | Elektromagnetsystem fuer die kernspintomographie |
| JPH0690975B2 (ja) * | 1988-06-03 | 1994-11-14 | 三菱電機株式会社 | 磁場補正用磁性体シム |
| US5227728A (en) * | 1991-11-01 | 1993-07-13 | The Regents Of The University Of California | Gradient driver control in magnetic resonance imaging |
| DE4217496C2 (de) * | 1992-05-27 | 1994-06-16 | Bruker Analytische Messtechnik | Shim-Verfahren |
| US5313164A (en) * | 1992-11-27 | 1994-05-17 | Resonance Research Inc. | Apparatus for mapping static magnetic fields |
| EP1647831B1 (en) * | 1996-10-30 | 2008-12-03 | Hitachi Medical Corporation | Open superconducting magnet apparatus |
| US5770943A (en) * | 1996-12-30 | 1998-06-23 | General Electric Company | Method for measuring and compensating for spatially and temporally varying magnetic fields induced by eddy currents |
-
2000
- 2000-04-19 JP JP2000117880A patent/JP3701540B2/ja not_active Expired - Fee Related
-
2001
- 2001-04-05 US US09/827,037 patent/US6813513B2/en not_active Expired - Fee Related
- 2001-04-12 EP EP01303473A patent/EP1148342A3/en not_active Withdrawn
- 2001-04-18 CN CN01119654A patent/CN1336558A/zh active Pending
- 2001-04-18 KR KR10-2001-0020699A patent/KR100458771B1/ko not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| KR100458771B1 (ko) | 2004-12-03 |
| US6813513B2 (en) | 2004-11-02 |
| EP1148342A2 (en) | 2001-10-24 |
| KR20010098694A (ko) | 2001-11-08 |
| US20010041819A1 (en) | 2001-11-15 |
| EP1148342A3 (en) | 2003-09-03 |
| CN1336558A (zh) | 2002-02-20 |
| JP2001299721A (ja) | 2001-10-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3701540B2 (ja) | 磁場測定方法、勾配コイル製造方法、勾配コイルおよび磁気共鳴撮影装置 | |
| JPS5946546A (ja) | 核磁気共鳴による検査方法及び検査装置 | |
| JPH0224135B2 (ja) | ||
| JP2003111746A (ja) | Ｔ２／ｔ２＊の重み付けを低減した中位の磁場強度における水分および脂肪のｍｒ画像の区別および識別 | |
| USRE32701E (en) | NMR scanner with motion zeugmatography | |
| US20090253983A1 (en) | Image based measurement of contrast agents | |
| CN103969610A (zh) | 采集磁共振数据和确定b1磁场的方法及磁共振设备 | |
| EP0103372A2 (en) | Nuclear magnetic resonance method and apparatus | |
| JP4493763B2 (ja) | 磁気共鳴イメージング装置及び画像サイズ可変装置 | |
| CN103364747A (zh) | 用于控制磁共振系统的方法和控制设备 | |
| EP0161366B1 (en) | Nuclear magnetic resonance methods and apparatus | |
| CN103454606A (zh) | 用于确定磁共振技术中特定于对象的b1分布的方法 | |
| JP2023156431A (ja) | 組織を撮像するためのシステムおよび方法 | |
| JP3705995B2 (ja) | 勾配コイル製造方法および勾配コイル並びに磁気共鳴撮影装置 | |
| JPH0357774B2 (ja) | ||
| JP3705996B2 (ja) | 磁気共鳴撮影装置 | |
| JPS61100617A (ja) | 核磁気共鳴を用いた流体運動測定方法及び装置 | |
| JPS63216557A (ja) | 磁気共鳴映像装置 | |
| JPH03224538A (ja) | 一次の静磁場不均一を補正して計測する過程を備えたｍｒｉ装置 | |
| JP4247511B2 (ja) | 勾配磁場測定方法および装置並びに磁気共鳴撮影装置 | |
| JPH0685766B2 (ja) | 流体流速測定装置 | |
| JPS62129043A (ja) | 核磁気共鳴画像化方法ならびに装置 | |
| JPH04343833A (ja) | 静磁場不均一の一次項補正の可能な磁気共鳴イメージン グ装置 | |
| JP2004254884A (ja) | 磁気共鳴撮影装置 | |
| JP2005058252A (ja) | 磁気的状態表示方法および磁気共鳴撮影装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040914 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20041202 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050614 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050713 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090722 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090722 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100722 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100722 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100722 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110722 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110722 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120722 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120722 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120722 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120722 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130722 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130722 Year of fee payment: 8 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |