JPS6288949A - 核磁気共鳴を用いた磁化分布の計測方法 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔発明の利用分野〕 本発明は核磁気共鳴を用い、対象物体中の該スピン分布
、緩和時間分布など（以下磁化分布と呼ぶ）を非破壊的
に求める検査装置に係り、特に対象の分布を高速で画像
とする装置に関する。

〔発明の背景〕

本発明に最も近い公知技術としては、例えはジャーナル
・オブ・ブイジックス・シー（Ｊ、　Ｐｈｙｓ。

Ｃ）、ソリッド・ステイＩ−・ンイジックス（Ｓｏｌｉ
ｄＳｔａｔｅ　Ｐｈｙｓ、　）第１０巻、第５５〜５８
頁に紹介されているＰ、マンスフイールド氏提案のエコ
ープラナ−法がある。この方法を実現するための典型的
なＲＴ−’および傾斜磁界の印加タイミングを第１図に
示す。図はＺ軸に垂直な断面を選択し計泪すする場合の
シーケンスであり、図中Ｇｙのように、１方向に定常的
に印加する傾斜磁界、及び図中Ｇｘのようにこれと直角
方向に振動する傾斜磁界を印加することによりスピン分
布を測定する。この方法では大振幅かつ高速の傾斜磁場
の切り変えが必要であるが、原理的には１回の９０°パ
ルス印加により８１１１定されるＮＭＲ信号を用い像を
構成できる。

しかし、特別な後処理を行わない限り、振動する傾斜磁
界の撮動波形として、矩形波を用いなければならず、こ
れはこの傾斜磁界発生用電源に対して大きな負担となり
、また、大きな周波数が必要とされる場合、傾斜磁場発
生用コイルも特別な工夫を必要とする。

〔発明の目的〕 本発明の目的は、核磁気共鳴を用い、対象物体中の磁化
分布を非破壊的に求める検査装置において、矩形的な傾
斜磁場を切り変えを用いずに、正弦波的な周期波形を用
いることにより、１）′６速で画像再構成に必要なデー
タを敗り込むことのできる装置６を提供することにある
。

〔発明の概要〕

２次元イメージングを仮定し、画像面を（χ＋　ｙ　）
面とする。傾斜磁場を Ｇｙ（ｔ）＝Ｇｏｓｉｎωｔ　＋ＧＯ（１１ｔｃｏｓω
ｔと駆動すると、得られる詠磁気共鳴信号は、５（ｔ）
＝、ｆｐ（ｘ＋ｙ）ｅｘｐ　　ｌ：１ｙ（Ｇｔ　　（ｘ
ｃｏｓωも＋ｙｓｉｎωｔ）　　３コｄｘｄｙとなる。

上式から得られるデータ列はρ（ｘ、ｙ）の空間周波数
領域を第２図に示すごとく螺線状にサンプリングしたも
のになる。このフーリエ領域のスパイラルサンプリング
については、ジャーナル・オブ・マグネティック・レゾ
ナンス（Ｊｏｕｒｎａｌｏｆ　ｍａｇｎｅｔｉｃ　ｒｅ
ｓｏｎａｎｃｅ）誌、第５４巻、３３８−３４３頁（１
９８３年）の中で示唆されているが、上記文献には、像
再生手法などの方法に関する示唆以上の議論は行われて
いない。

本発明は、螺線状に測定されたデータに刻する像再生手
法、この像再生を適用するためのデータ計測法および像
再生手法から導かれる傾斜磁場パｌ勅に対して負担の少
ないデータ１ｌｌ１１定法を特徴とする。

ます像再成について説明する。信号計＆ｌＩＩのサンプ
リング間隔、すなわち、」も体的には計？１ｌｌｌされ
た信号をディジタル化して計算機に与えるＡ　ｎ　ｇ換
器のサンプリング間隔Δｔを。

となるようにえらぶと、第２図に示されるごとく、ｋ８
軸と角度θをなす直線上に等間隔で測定点が゛ならぶよ
うになる。（同図中ではサンプリング点を丸印で示す。

） これらの直線を螺線が交差する時刻は、ω　　　　　　
　ω で表わされる。（Ｌ）式に代入すると ｓ　（ｔ　ｎ）＝、ｆｐ　（ｘ　ｅ　ｙ）　ｅｘｐ［ｉ
　ｙ　（Ｇｏｔｎ　’（ｘｃｏｓθ＋ｙｓｊｎθ）　）
］ｄｘｄｙとなる。ｔｎ′を連続変数とみなし、ｔ６と
表記すると。

Ｓ　（ｔ’　）＝＜　１　＋ｓｇｎ　　（ｔ　−”　）
）Ｊ”／ｌ　（ＸＩ　　ｙ）　　ｅｘｐｃｙ　（Ｇｎｔ
、６（ｘｃｏＳＯ＋　ｙ　Ｓｉｎθ））］ｄｘｄｙ　　
　　　　　　　　　　　（３）ここで ｓｇｎ（ｔ）＝　　　　１　　　ｔ＞０−１　　ｔｏ。

である。ここでＳ　（ｔ’）は、第２図において、角度
りであるような動径上にならんだデータ群を指す。以後
ｔ＃に関する演算を行うことは、ラセン上のデータをθ
の異なる動径上のデータ群にグループ分けしたことにな
る。このことを明確にするために、角度Ｏなる動径上の
データを、　５ａ（ｔ　）と表記する。（３）式の両辺
のフーリエ変換をとると、 π↑ ここで、串はコンボリューションを意味する。またＰｅ
　（ｆ）＝Ｊ−ｐ　（ｘ＋　ｙ　）δ（ｆ　−−Ｇｏ（
ｘｃｏｓ　０２π ＋ｙｓｉｎθ））ｄｘｄｙ　　　　　　　　　　　（５
）であり、このＰθ（ｆ）は周波数軸ヒに与えられたρ
（ｘ＊ｙ）のＯ方向の射影である。（：３）および（４
）式により、 ２＋１１１！ Ｐ、　（ｆ　）＝　ＲＱ　（ｃ’　ζ（ｓ、（ｔ’）ｌ
　）　　（６）式を得る。したがって、（６）式より角
度０の射影を求め、０〜３６０″のＰ　（ｆ）からρ（
ｘ、ｙ）を再構成する。

再構成は以下の様に行う。

Ｐ＃（ｆ）をｆに関してフーリエ変換する。これを１ｒ
ａ（Ｒ）と表記すると。

Ｆ＊（Ｒ）＝Ｊ”’Ｐ＃（ｆ）ｅ”ｌＨ’ｄｆ　　　（
６Ｌ）となり、このＦｌＩ（Ｒ）に対して、以ドの様な
計算を行えばρ（ｘ、ｙ）を再生できる。すなわち、＋
ｙｓｉｎＯ）］　ＲｄＲｄ　ｆｌ　　　　　　　（６−
２）である。ところで１通常（６−２）式の計算は２つ
に分けて行われる。すなわち、 をます計算し、 を求めろ方法である。ここで特に（６−４）式の演算は
逆投影と呼ばれている。ここで、（６）。

としてもＷ、（ｕ）を求めることができる。この（６−
５）式を用いる方法は、スパイラル上のデータを各角度
の動径上のデータにグループ分けした後、動径方向へ重
みづけの後フーリエ変換し。

Ｗｅ（ｕ）を直接求める方法であり、（６）、（ｅ　−
１）ｓ（６−３）式を用いる方法に比ベフーリエ変換の
計算を１角度分のデータにつき２回節約できる。

〔発明の実施例〕

以下１本発明の一実施例を画面により説明する。

第３図は１本発明の一実施例による核磁共鳴を用いた検
査装！ｉ’ｔ　（以下、甲−に、「検査装置」と呼ぶ）
の概略構成図である。

第３図において、１は静磁場Ｈｏを発生させる電磁石、
２は対象物体、３は高周波磁場を発生させると同時に、
対象物体２から生ずる信号を検出するためのコイル、４
　ｘ　ｇ　４　ｙ＋および５は、それぞれＸ方向、Ｙ方
向、Ｚ方向の傾斜磁！　’；光発生せるための傾斜磁場
発生コイルである。傾斜磁場発生コイル５としては、互
いに逆向きに電流が流れるように配線された円線輪を用
いる。６，７゜８はそれぞれ上記各傾斜磁場発生コイル
４ｘ。

４ｙ、５に電流を供給するためのＶ＠装置である。

９は計算機、】０は静磁場発生用の１℃磁石］のための
電源、１１は対象物容積計量装置である。傾斜磁場発生
コイル４ｘ、４ｙ、５により発生する傾斜磁場の強度は
、上記対象物容量計量′！Ａ置１１からの指令により変
化させることができる。

次に本検査装置の動作を概略図に説明する。

対象物体２の核スピンを励振する高周波磁場は、シンセ
サイザ１２により発生された高周波を変調装置１３で波
形整形・電力増幅し、コイル３に電流を供給することに
より発生させる。対象物体２からの信号はコイル２によ
り受信され、増幅装置１４を通った後、検波器１５で直
交検波され計算機９に入力される。計算機９は、信号処
理後、核スピンの密度分布、あるいは緩和時間分布に対
応する画像をＣＲＴディスプレイ】６に表示する。

１７は計算の途中のデータあるいは最終データを格納す
るメモリーである。

第４図は本発明の実施例の計測シーケンスを示し、ＲＦ
は高周波磁場波形、Ｇｚ、Ｇｘ、ＧｙはそれぞれＺｙ　
Ｘｒ　’／力方向傾斜磁場発形、ＡＤは信号サンプリン
グの期間を示）−０ まず９０°高周波パルスと傾斜磁場Ｇ　ｚ　Ｑ）　ｔ４
Ｊ加により、被検体の所望スライス内のスピンを耐起す
る。次に９０℃高周波パルスのピークから時間τが経過
した時点で、［８０°高周波パルスを印加し、これによ
り再び時間τ経過した点（ｔ＝０）でのスピンの横磁化
信号が得られるようにする。

１＝０の時点から（１）式に示す波形で傾斜磁場Ｇｘ、
Ｇｙを印加し、信号のサンプリングを開始する。

前に述べたごと＜Ａｒ）のサンプリング間隔Δしされた
データを とグループ分けする。すなわち、角Ｊ皮の異なる動径上
のデータごとにグループ分けする。角度Ｏ１なる動径上
のデータはｂａｇと表記する。ここで’　ａ　　　　Ｉ
ｎ　Δ　０　　　　２　　ｎ　　πし、　　＝　　−＋
　　− ω　　　　　　　　　　　　　ω であり、ＡＤ　＝　２　ｔｃ　／　Ｍである。

各ｍについてデータ列Ｓθｍ（己”）を、ｎについて離
散フーリエ変換し、（６）式にしたがって、位相補正の
後、実部を取ってｐｇ、（ｆ）を求める。この計３）゛
を各々のｍについてそれぞれ行い、各角度におけるＰａ
、（ｆ）から（６−４）−（６−４）式、もしくは（６
−５）式にして説明した逆投影法により像を再構成する
。

以りのシーケンスにより、従来のエコープラナ−法のよ
うな矩形状のＭ斜磁場の振動を行うことなく、正弦波的
な周期波形の傾斜６Ｊ＆場の印加により、被検体のスピ
ン密度分布を示す磁化分布像が得られる。なお被検体の
スピンの縦緩和時間ＴＩの分布の効果を含む磁化分布像
（Ｔ’１強調像）を得るには、良く知られているように
９０°高周波パルスの印加から所定間前に１８０°高周
波パルスを印加して予めスピンの反転を行い、１＝０の
時点で縦緩和の効果を含んだ横磁化信号が１！）られる
の間隔でサンプリングした場合には、サンプリングされ
たデータ間の補間により予め角度Ｏｓ　（ｍ＝０．１．
・・・、Ｍ−１）の動径上のデータを求めてからフーリ
エ変換すれば良いに の像再生手法にはいくつかの変形手法が考えられる。ま
ず、（６）式を用いるかわりに５ａ（ｔ’）に対して重
みづけを行いコンボリューション後の射影を直接求めろ
方法である。

コンボリューションウェイトの周波数特性をψ（＋、）
とすると、コンボリューション後の射影Ｗ　ｅ　（ｆ　
）は Ｗ　　（ｆ）＝Ｒｅ（ｅ”””　　ζ　（Ｙ（ｔ）Ｓ、
（ｔ’））　　）　　　　　　　（’７）により求まる
。この（７）式は、前記（６−５）式に対応している。

ここで、’Ｐ　（ｔ、　）はスピンの横緩和（Ｔ　２　
ｄｅｃａｙ）の影響を考慮すれば、試料の平均的なＴｚ
を〒２として ’ｆ’　（ｔ、　）ωｔｅＴ２ とすることにより、横緩和の影響をある程度軽減できる
。

ところで、視野内における静磁場の不均一性を考慮する
場合、逆投影を工夫することによりこの影響を補正でき
る。以下これを説明する。

視野内における静磁場の不均一性をＥ（ｘ、ｙ）とする
と、これを考慮した場合、（１）式はＳ’　（ｔ　）＝
ｆ／）　（Ｘ　ｑ　ｙ）ｅｘｐ［ｉｙ　（Ｇｅｔ（ｘｃ
ｏｓωｔ＋ｙｓｊｎωｔ＋Ｅ（ｘ＋　ｙ）ｔ））コｄｘ
ｄｙ　　　　　　　　　　（１）’となり、同様な導出
により、 Ｐａ’　　（ｆ）＝＝Ｒｅ　　（ｅ”／＋ζ［８’　　
（ｔｎ’）コ）であり、 を得る。（３）’式はＥ　（ｘ　ｔ　ｙ　）が射影デー
タに周波数シフトを引き起こすことを示しており、逆投
影時にｈ　（ｘ　＋　ｙ　）を考慮してチー′夕を選ぶ
ことにより補正できろ。この補正について以下に説明す
る。

逆投影時における像面と周波数軸の関係を第６１閑に示
す６同図でｆｃは、視野中心における共鳴周波数である
。通常の逆投影では、図に示される通り、角度０から点
（ｘ＋ｙ）への逆投影データは次のｉＢで示される周波
数位置から取られる。

ＪＲ＝　γＧ（ｘｓｉｒ＋１？　＋　ｙｃｏｓＯ）／　
２７Ｃ＋ｆ。

けれども、静磁場の不均一性の影ｑ＞で点（ｘ、ｙ）の
共鳴周波数は上式のｆ８ではなく、ｆＢ’＝ｆｎ＋ｆｅ
で与えられろｆＢＩ　となる。したがって、磁場不均一
性が既知であれば、このＪ１１′　を逆投影の座標計算
において１口の代わりに用いることにより静磁場不均一
性の影響を補正できる。

さて、以ヒの説明ではＯ〜３６０°におけるＰ、（ｆ）
から像再生を行う方式を説明したが、Ｓ＃（ｔｎ”）と
Ｓ　ｌ＋＋、　（ｔ　ｎ　”’　）をつなぎ合わせてフ
ー工変換することによりｐ、（ｆ）を求めれば、この場
合にはＯ〜１８０°までのＰａ（ｆ）により像を再生で
きる。ただし、第２図に示される通り、５６（ｔ、ｎ’
）ｖ　５ｔｒａｙ　（ｔｎ”’　）とも原点付近で測定
点が不等間隔となるため、補間により等間隔データ点を
推定しなければならない。

さて、ω＝　２　ＦＣ／　ＴなるＴを導入し、視野の半
径をｒくとすると、 ２π　　　　　２Ｒ の関係がある。画像マトリックスを、Ｉ　Ｘ　、Ｊとす
ると、分解能Δγ＝２Ｒ／Ｊであり同時に、 Δγ＝　−（８） γＧｏｔｓａｘ であるので、これらの式より Ｊ ｔ、□＝　−ｒ　　　　　　　　　　　　　　　（９）
を得る。このときＧ８およびＧ、の最大振幅Ｇ＋＋ａｘ
は Ｇｍａｘ＝　ＧＯ（１１ｔ　ｍａｘ＝　　πＪＧｏ　　
　　　　　　　　　　　　　（１０）で与えられる。例
えばＲ＝　１５ｃｍ、　Ｔ　＝０．２５ｍ５゜Ｊ　＝　
１．２８とすると、Ｇｏ＝０．Ｏ］６ガウス／Ｃｍ。

Ｇ、、、＝６．４ガウス／ａｍとなる。このＴと６１１
８の値の組み合わせを現実の装置で実現するのは非常に
困難である。

この困難をさけるには、イ頃斜磁場Ｇｙ、Ｇｙを。

次式で表せるものとする。

ここで、この傾斜磁場波形の要素となる正弦波。

余弦波の位相φを等間隔で変更して複数回第４図のシー
ケンスの計ｉ！＋’ｌをくり返し、等間隔の複数本の螺
線上のデータの組を得る。次に第４図の実施例で説明し
たと同様に、各角度の動径上のデータごとにグループ分
けを行いそれぞれフーリエ変換して像を再構成すればよ
い。第５図にφが−ずつ異なる４本の螺線ヒの信号をサ
シブリングする例を示す。

今、１〕本のＴ！！Ａ線を用いて像を再構成したとする
と１ 、■ ｔ　ｗａｘ　”　　　Ｔ’ Ｄ Ｇｓａｘ＝−πＪＧ。

となり、Ｇｍａｘに対する条件がＤ倍緩和されているこ
とがわかる。

また、本発明け、ＮＭＲ顕′ｐＩ鏡などのように。

被検体を回転させることができる場合に、特に容易に実
現できろ。

すなわち、このような場合を第７図に示す。同図におい
て、被検体は視野内で同転しており、回転中心と視野中
心は一致している。実験室に固定した座標系を（ｘ、ｙ
）とする。第４図に示す回転磁界Ｇｘ、Ｇｙに変え、第
８図に示すように、で示されるに１斜磁場Ｇ　Ｋ　、　
Ｇ　ｙを印加する。

すなわち、江いに直交する２つの傾斜磁界を。

１つは定常的に、１つは時間に比例して、印加し。

その比例定数を定祁゛的に印加している傾斜磁界の大き
さでおった１１′Ｈに！１１当する角速度Ｓ彼倹体を回
転させろ。このようにすると１回転している被検体に固
定されたＰＰ、漂糸を（ｘ　ｒ　＋　ｙ　ｒ　）とする
と、となる。こ、：で、 を（１４）式に代入すると ５（ｔ）＝ｆｐ（ｘｒ、ｙｒ）ｃｘｐ［ｊｙ（Ｇｎｘｒ
、ｆ　’（ｃｏｓｗｔ、’−ｕｔ’　ｓｉｎωｔ’　）
ｄｔ’ ＋ＧｏｙＪ”　（ｓｉｎωｔ’ ＋　ｔａｓ　ｔｃｏｓ　（１１ｔ’　　）ｄｔ）コｄ　
ｘ　ｒｄ　ｙ　ｒ　　　　　　　　　　　　　　　（１
６）であり、積分を計算すると。

５（ｔ）＝Ｊ’ρ（ｘｒｐｙｒ）ｅｘｐ［１ｙ（ｘｒＧ
ｏｔｃｏｓωｔ＋ｙｒＧｏｔｓｉｎωｔ）ｌｄｘｒｄｙ
ｒ　　　（１７）となる。すなわち、被検体の回転と第
８図の傾斜磁場の印加により第２図に示すと同様なスパ
イラルサンプリングが実行できる。

したがって、信号のサンプリング間隔を被検体の回転周
期の′Ｍ数分の１とし、得られたサンプリング信号のグ
ルービング、及び像再構成処理を第４図にて説明したの
と同様に行うことにより被検体のスピン密度分布が得ら
れる。また信号のサンプリング間隔が被検体の回転周期
の整数分の１のでない場合は補間により回転周期の整数
分の１に対応する信号を算出してからグルービング、及
び像再構成処理を行えば良い。

さらに、第８図に示す傾斜磁場Ｇ　ｘ　、　Ｇ　ｙの印
加タイミングを被検体の回転角に対して変化させて複数
回の測定を行えば、第５図に示すように、空間周波数領
域において複数の螺線上のサンプリングデータが得られ
る。このデータからも、像面構成が行えることも先に示
した実施例と同様である６〔発明の効果〕 以上に示した本発明によれば、フーリエ空間を螺線状に
サンプリングするとの新たな手法を取り入れたことによ
り、矩形的な傾斜磁場の高速振動乞行うことなく、高速
で画像再構成の必要なデータを取り込むことができ、対
象物体中の該スピン分布もしくは緩和時間分布などを非
破壊的に求める装置の装置上の制約を緩和できる。

図面のｆｆ？１＠な説明 第１図は本発明に近い公知技術であるエコープレチー法
を実現するパルスシーフェンス、第２図は本発明による
測定点のならび方の一例、第３図は本発明の実施例の装
置４９成例、第４図は本発明の実施例のＲＦおよび傾斜
磁場の印加シーフェンス、第５図は本発明による測定点
のならび方の別の一例、第６図は逆投影時における像面
と周波数軸の関係、第７図は被検体を同転する本発明の
他の実施例の座部系、第８図は第７図の実施例のパルス
シーケンスである。

Claims (1)

1. 【特許請求の範囲】 １、ほぼ均一の静磁場中におかれた被検体の所望領域の
   横磁化信号を準備する第１の過程 前記横磁化信号と領斜磁場の印加のもとに計測する第２
   の過程 計測された信号から磁化分布像を再構成する第３の過程
   を含み、 前記傾斜磁場は正弦波あるいは余弦波関数の組合せによ
   り構成される振動波形であることを特徴とする核磁気共
   鳴を用いた磁化分布の計測方法。 ２、前記振動波形の積分値が空間周波数領域において螺
   線を描くことを特徴とする特許請求の範囲第１項に記載
   の核磁気共鳴を用いた磁化分布の計測方法。 ３、前記振動波形の構成要素である正弦波もしくは余弦
   波関数の周期を、前記第２の過程における信号のサンプ
   リング間隔で割つた値が整数であることを特徴とする特
   許請求の範囲第１項に記載の核磁気共鳴を用いた磁化分
   布の計測方法。 ４、前記第３の過程は、前記第２の過程においてサンプ
   リングされた信号からの補間により、前記正弦波あるい
   は余弦波関数の周期との比が整数である間隔の測定点の
   信号を算出する過程を含むことを特徴とする特許請求の
   範囲第１項に記載の核磁気共鳴を用いた磁化分布の計測
   方法。 ５、前記第３の過程は、前記螺線の上に分布する計測点
   のデータを動径方向ごとの１次元データにグループ分け
   する過程、該１次元データに重み付けをした後にフーリ
   エ変換する過程、フーリエ変換の結果を逆投影すること
   により画像を再構成する過程を含む特許請求の範囲第２
   項に記載の核磁気共鳴を用いた磁化分布の計測方法。 ６、前記逆投影に際して、逆投影データの座標計算に視
   野内の静磁場分布に関するデータを用いることを特徴と
   する特許請求の範囲第１項に記載の検査装置。 ７、前記振動波形の構成要素である正弦波あるいは余弦
   波関数の位相を変更して前記第１、第２の過程を複数回
   くり返し行うことを特徴とする特許請求の範囲第１項も
   しくは第２項に記載の核磁気共鳴を用いた磁化分布の計
   測方法。 ８、ほぼ均一の静磁場中におかれた被検体の所望領域の
   横磁化信号を準備する第１の過程、 前記被検体を回転させ、所定の傾斜磁場の印加のもとに
   前記横磁化信号を計測する第２の過程、 計測された信号から磁化分布像を再構成する第３の過程
   を含み、 前記傾斜磁場は一定の第１傾斜磁場成分と、該第１傾斜
   磁場成分と直交する方向に印加し、一定の増分で大きさ
   の変化する第２傾斜磁場成分からなることを特徴とする
   核磁気共鳴を用いた検査装置。 ９、前記傾斜磁場は前記回転する被検体に関する空間周
   波数領域において螺線を抽くことを特徴とする特許請求
   の範囲第８項に記載の核磁気共鳴を用いた磁化分布の計
   測方法。 １０、上記被検体の回転周期を上記第２の過程における
   信号のサンプリング間隔で割つた値が整数であることを
   特徴とする特許請求の範囲第８項に記載の核磁気共鳴を
   用いた磁化分布の計測方法。 １１、前記第３の過程は、前記第２の過程においてサン
   プリングされた信号からの補間により、前記被検体の回
   転周期との比が整数である間隔の測定点の信号を算出す
   る過程を含むことを特徴とする特許請求の範囲第８項記
   載の核磁気共鳴を用いた磁化分布の計測方法。 １２、前記第３の過程は、前記螺線の上に分布する計測
   点のデータを動径方向ごとの１次元データにグループ分
   けする過程、該１次元データに重み付けをした後にフー
   リエ変換する過程、フーリエ変換の結果を逆投影するこ
   とにより画像を再構成する過程を含む特許請求の範囲第
   ９項に記載の核磁気共鳴を用いた磁化分布の計測方法。 １３、前記傾斜磁場の印加タイミングを前記被検体の回
   転に対して変化させて前記第１、第２の過程を複数回く
   り返し行うことを特徴とする特許請求の範囲第８項もし
   くは第９項に記載の核磁気共鳴を用いた磁化分布の計測
   方法。