JP2791091B2 - 核磁気共鳴を用いた検査装置 - Google Patents

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は、核磁気共鳴（以下、「NMR」という）を用
いた検査対象内の水素あるいはリンなどの分布を映像化
する方法及び装置に係り、特に短時間で定量的な化学シ
フトイメージングを実行するのに必要なすべての信号を
測定するのに好適な、NMRを用いた検査方法及び装置に
関する。

〔従来の技術〕

静磁場、傾斜磁場、及び高周波磁場の各磁場を発生す
る手段と、検査対象からの核磁気共鳴信号を検出する信
号検出手段と、この信号検出手段による検出信号（核磁
気共鳴信号）の演算処理を行なう計算機とを有し、検査
対象の核スピンの分布情報を得る核磁気共鳴を用いた検
査装置は、従来技術として周知であり、例えば、特開昭
61−226648号公報、特開昭61−62851号公報、特開昭58
−41340号公報、特開昭61−205853号公報、特開昭61−1
3143号公報等に記載がある。検査物体からの信号を、物
体各部に対応させて分離識別する必要があるが、その一
つに、検査物体に傾斜磁場を印加し、物体各部の置かれ
た静磁場を異ならせ、これにより各部の共鳴周波数ある
いはフェーズエンコード量を異ならせることで位置の情
報を得る方法が広く知られている（例えば、特開昭61−
226648号公報、特開昭61−62851号公報、特開昭61−205
853号公報、特開昭61−13143号公報等）。さらに、これ
ら従来装置における、90゜、及び180゜高周波磁場、３
方向（ｘ、ｙ、ｚの各方向）の傾斜磁場の各磁場の印加
の順序は、それぞれ特徴あるパルスシーケンスにより、
装置の動作が規定されている（例えば、特開昭61−2266
48号公報、特開昭61−62851号公報、特開昭61−205853
号公報、特開昭61−13143号公報等）。これらのパルス
シーケンスは、周知のように、核磁気共鳴信号の演算処
理を行なう計算機により制御されるか（例えば、特開昭
61−226648号公報、特開昭58−41340号公報）、核磁気
共鳴信号の演算処理を行なう信号処理装置とは別の制御
装置により制御されるか（例えば、特開昭61−62851号
公報等）、あるいはCPUからの制御により動作するシー
ケンサ（例えば、特開昭61−205853号公報、特開昭61−
13143号公報等）により、制御されている。特開昭58−4
1340号公報の第２図には、180゜パルスを連続して印加
して、複数回にわたり、核磁気共鳴信号を検出するパル
ス（系列）シーケンスが記載されている。

従来、マルチエコーシーケンスを用いて化学シフトイ
メージングを行なう方法については、ラヂオロジー，第
159巻（1986年）第783頁から第786頁（Radiology,vol.1
59,（1986）PP783−786）において論じられている。

上記方法は、デイクソン（Dixon）法（ラジオロジイ
（Radiology）、第13巻（1984）第189頁〜第194頁）を
基にしており、実行にあたつては、検査対象に対してマ
ルチエコーシーケンスにより２回の測定を行ない、それ
ぞれの測定で第１エコー信号、第２エコー信号を取得
し、計４枚の画像より化学シフトイメージングによる分
離画像を得ている。

〔発明が解決しようとする課題〕

上記従来技術は、２回の測定が必要であるという点
で、検査対象からの信号測定に関する時間の短縮につい
て配慮がされておらず、生体を検査対象とする場合に
は、被検者に苦痛を与える。測定時間内における検査対
象の動きにより、最終的な結果である分離画像の定量性
が劣化する。分離画像を得るための画像処理が複雑であ
る、などの問題があつた。

本発明の目的は、化学シフトイメージングによる分離
画像を得るために必要な検査対象に対する信号測定に要
する時間を短縮し、上記問題をより軽減することにあ
る。また、本発明の他の目的は、分離画像を得るための
画像処理を容易にすることにある。

〔課題を解決するための手段〕

上記目的を達成するために、本発明においては、通常
のマルチエコーシーケンスに対して化学シフトイメージ
ングを行なうための特定の時間差をもたせたタイミング
でNMR信号を検出するシーケンスにより、検査対象の
み、あるいは、検査対象ならびに一様な化学シフトをも
つ物質で満たされた一様フアントムを用いて第１エコー
信号，第２エコー信号までの測定を実行する。

即ち、本発明の特徴は、（Ａ） 静磁場、傾斜磁場、
及び高周波磁場の各磁場を発生し、静磁場が印加された
空間に置かれた検査対象からの核磁気共鳴信号を検出
し、検出された核磁気共鳴信号に対して演算処理を行な
い、検査対象の核スピンの分布情報を得る核磁気共鳴を
用いた検査装置において、（１） （ｉ）検査対象の核
スピンを励起するステップと、（ii）エンコード傾斜磁
場を印加するステップと、（iii）２つのケミカルシフ
トをもつ核スピンの間で所定の位相差θ_１がθ_１≠ｎπ
（ｎ＝０、±１、±２、…）となる時間差をもたせたタ
イミングで、検査対象に第１の180゜高周波磁場パルス
を印加するステップと、（iv）リードアウト傾斜磁場を
印加し、第１の核磁気共鳴信号を計測するステップと、
（ｖ）２つのケミカルシフトをもつ核スピンの間で位相
差θ_２がθ_２＝2nπ（ｎ＝±１、±２、…）となる時間
差をもたせたタイミングで、検査対象に第２の180゜高
周波磁場パルスを印加するステップと、（vi）リードア
ウト傾斜磁場を印加し、第２の核磁気共鳴信号を計測す
るステップと、及び（vii）（ｉ）から（vi）のステッ
プを、位相エンコード量を変化させてエンコード傾斜磁
場を印加し、繰返すステップと、を有するパルスシーケ
ンス制御を実行して、第１、第２の核磁気共鳴信号の計
測を行ない、第１、第２のデータセットを得る第１のス
テップと、（２）２つのケミカルシフトのうちの１つを
もつファントムを、静磁場が印加された空間に置いて、
上記のステップ（１）と同じパルスシーケンス制御を実
行して、第３、第４のデータセットを得る第２のステッ
プと、（３）第１、第２のデータセットをそれぞれフー
リエ変換して、第１、第２の画像データ（S₁、S₂）を求
める第３のステップと、（４）第３、第４のデータセッ
トをそれぞれフーリエ変換して、第３、第４の画像デー
タ（U₁、U₂）を求める第４のステップと、（５）第３、
第４の画像データの位相分布をそれぞれ補正データとし
て使用し、第１、第２の画像データからそれぞれ静磁場
の不均一に起因する位相誤差分布を除去した第５、第６
の画像データ（S₁′、S₂′）を求める第５のステップ
と、及び（６）第６の画像データから計算される位相分
布（Θ_２′）を補正データとして使用し、第５の画像デ
ータから、検査対象の透磁率分布（E_per）に起因する位
相誤差分布を除去することにより、２つのケミカルシフ
ト（σ_１、σ_２）の核スピンの分布情報を含む第７の画
像データを求める第６のステップと、を有する核磁気共
鳴を用いた検査方法に特徴がある。

（Ｂ） 静磁場、傾斜磁場、及び高周波磁場の各磁場を
発生する手段と、静磁場が印加された空間に置かれた検
査対象からの核磁気共鳴信号を検出する信号検出手段
と、信号検出手段の検出信号の演算処理及びパルスシー
ケンスの制御を行なう信号処理制御手段とを有し、検査
対象の核スピンの分布情報を得る核磁気共鳴を用いた検
査装置において、信号処理制御手段は、（１）（ｉ）検
査対象の核スピンを励起するステップと、（ii）エンコ
ード傾斜磁場を印加するステップと、（iii）２つのケ
ミカルシフトをもつ核スピンの間で所定の位相差θ_１が
θ_１≠ｎπ（ｎ＝０、±１、±２、…）となる時間差を
もたせたタイミングで、検査対象に第１の180゜高周波
磁場パルスを印加するステップと、（iv）リードアウト
傾斜磁場を印加し、第１の核磁気共鳴信号を計測するス
テップと、（ｖ）２つのケミカルシフトをもつ核スピン
の間で位相差θ_２がθ_２＝2nπ（ｎ＝±１、±２、…）
となる時間差をもたせたタイミングで、検査対象に第２
の180゜高周波磁場パルスを印加するステップと、（v
i）リードアウト傾斜磁場を印加し、第２の核磁気共鳴
信号を計測するステップと、及び（vii）（ｉ）から（v
i）のステップを、位相エンコード量を変化させてエン
コード傾斜磁場を印加し、繰返すステップと、からな
り、第１、第２の核磁気共鳴信号の計測を行ない、第
１、第２のデータセットを得るための核磁気共鳴信号を
計測するパルスシーケンス制御と、（２）（ａ）第１、
第２のデータセットをそれぞれフーリエ変換し、第１、
第２の画像データ（S₁、S₂）を求めること、（ｂ）２つ
のケミカルシフトのうちの１つをもつファントムを、静
磁場が印加された空間に置いて、上記の（１）と同じパ
ルスシーケンス制御を実行して、第３、第４のデータセ
ットを得て、第３、第４のデータセットをそれぞれフー
リエ変換して、第３、第４の画像データ（U₁、U₂）を求
めること、（ｃ）第３、第４の画像データの位相分布を
それぞれ補正データとして使用し、第１、第２の画像デ
ータからそれぞれ静磁場の不均一に起因する位相誤差分
布を除去した第５、第６の画像データ（S₁′、S₂′）を
求めること、及び（ｄ）第６の画像データから計算され
る位相分布（Θ_２′）を補正データとして使用し、第５
の画像データから、検査対象の透磁率分布（E_per）に起
因する位相誤差分布を除去することにより、２つのケミ
カルシフト（σ_１、σ_２）を核スピンの分布情報を含む
第７の画像データを求めること、を含む演算処理を行な
う核磁気共鳴を用いた検査装置に特徴がある。

また、上記（Ａ）、（Ｂ）の特徴に加えて、上記ファ
ントムは前記検査対象の領域よりも大きい断面を有し、
上記２つのケミカルシフトのうちの１つのみからなる均
一な物質を含むこと、上記第７の画像データ（Ｒ）のピ
クセル位置に依存しないオフセット位相（θ_off）を除
去すること、オフセット位相（θ_off）は、上記第７の
画像データ（Ｒ）のピクセルのうち、上記２つのケミカ
ルシフトのうちの１つのみが存在するピクセルのデータ
の位相から検出されることにも特徴がある。

〔作用〕

検査対象として、化学シフトσ_１およびσ_２（単位は
［ppm］）の２種類の物質よりなる物体を考える。この
物体をマルチエコーシーケンスで測定する際に、第１エ
コーでは、化学シフトσ_１を持つスピンと化学シフトσ
_２をもつスピン間に θ_１≠ｎπ（ｎ＝0,±1,±２…） …（１） なる角度（位相差）を導入するような時間差をもたせ
たタイミングでNMR信号を検出する。

この信号より得られる画像の位相情報には、静磁場不
均一、つまり外部から印加した静磁場の不均一E_St（x,
y）と検査対象の透磁率分布に起因する磁束密度変化E
_per（x,y）の和の情報、および化学シフトの情報が含ま
れている。また、第２エコーでは、上記スピン間に θ_２＝2nπ（ｎ＝±1,±２…） …（２） なる角度（位相差）を導入するような時間差をもたせ
たタイミングでNMR信号を検出する。この信号より得ら
れる画像の位相情報には、静磁場不均一のみの情報が含
まれている。この情報を用いて第１エコー信号より得ら
れる画像より静磁場不均一の影響を除すれば、化学シフ
トイメージングによる分離画像を得ることができる。以
上のように検査対象に対するマルチエコーシーケンスに
よる１回の測定で化学シフトイメージングによる分離画
像を得ることができる。

上記した、第２エコー信号より得られる画像から抽出
した静磁場不均一の情報は、そのままでは第１エコー信
号より得られる画像に適用することはできない。すなわ
ち、第２エコー信号からの静磁場不均一の情報を、第１
エコー信号の測定されたタイミングにおける静磁場不均
一の情報に変換する必要がある。

この操作は一般に複雑かつ困難であるという問題があ
る。この操作を容易かつ確実なものとするためには、検
査対象に適用したのと同一のシーケンスにより、一様フ
アントムを測定する。この測定より得られた画像には、
第1,第２エコーそれぞれのタイミングにおける外部から
印加した静磁場の不均一E_St（x,y）の情報が含まれてい
ると考えてよい。一般に静磁場不均一Ｅ（x,y）は、E_St
（x,y）と、検査対象の透磁率分布に起因する磁束密度E
_per（x,y）の和で表わされ、 E_St（x,y）≫E_per（x,y） …（３） なる関係が成立つ。従がつて、検査対象から得られた画
像から、一様フアントムで得られた画像を用いてE
_St（x,y）を除去することにより、静磁場不均一除去の
操作が容易となる。これにより、分離画像を得るための
画像処理を容易とすることができる。

〔実施例〕

以下、本発明の一実施例を第１図により説明する。第
１図は、マルチエコーシーケンスにより化学シフトイメ
ージングを行なうためのシーケンスの一例である。第１
図に点線で示すように通常のマルチエコーシーケンスに
対して180゜パルスとスライス選択傾斜磁場をτ₁,τ_２
だけ移動して印加しシーケンスである。

検査対象として、化学シフトσ_１およびσ_２（単位は
［ppm］）の２種類の物質よりなる物体を考える。第１
図のシーケンスで測定したときに第1,第２エコーからそ
れぞれ得られる再構成画像S₁（x,y）,S₂（x,y）は式
（４）で示すことができる。

ここでρ₁,ρ_１′は化学シフトσ_１を持つ物質の密度
（ρ_１＞ρ_１′），ρ₂,ρ_２′は化学シフトσ_２をもつ
物質の密度（ρ_２＞ρ_２′），γは磁気回転比，
θ_１′，θ_２′は座標に依存しない一定位相,E（x,y）
は静磁場不均一でＥ（x,y）＝E_St（x,y）＋E_per（x,y）
と表わすことができる。またB₀は印加静磁場の空間的な
平均値である。

τ₁,τ_２を調整して、θ₁,θ_２が式（１），（２）を
満足するようにする。具体的な値として、θ_１＝π/2、
θ_２＝２πとする。この場合、式（４）は、式（５）の
ようになる。式（５）S₂（x,y）の位相 項Θ_２（x,y）は、式（６）で表わすことができ Θ_２（x,y）＝−γ〔Ｅ（x,y）＋σ_１〕τ_２＋θ_２′ …（６） る。式（６）より、新たに式（７）の画像を作る。

S₁（x,y）と、Ｓ（x,y）より、静磁場不均一 Ｓ（x,y）＝exp｛ｊ（τ₁/τ_２）Θ_２（x,y）｝ ＝exp｛−ｊγ〔Ｅ（x,y）＋σ_１〕τ_１＋
（τ₁/τ_２）θ_２′｝ …（７） Ｒ（x,y）＝S₁（x,y）×Ｓ＊（x,y） ＝｛ρ_１（x,y）×ｊρ_２（x,y）｝ ×exp｛ｊ（θ_１′−（τ₁/τ_２）
θ_２′）｝ …（８） が除去された画像式（８）が得られる。なお、＊は複素
共役を表わす。例えば、Ｓ＊（x,y）は、Ｓ（x,y）の複
素共役を表わす。式（８）の下線部のオフセツト位相を
除去すれば目的とする分離画像Ｔ（x,y）が得られる。
そのためには、ρ_１＝0,ρ_２≠０、あるいは、ρ_１≠0,
ρ_２＝０ Ｔ（x,y）＝ρ_１（x,y）＋ｊρ_２（x,y） …（９） なる（x,y）における位相θ_offを調べればよい。すなわ
ち、式（８）よりρ_１≠0,ρ_２＝０なる（x,y）では、
式（10）となり、式（10）の位相項Θ_Ｒ′は Ｒ′（x,y）＝ρ_１（x,y）exp｛ｊ（θ_１′−（τ₁/
τ_２）θ_２′｝ …（10） 式（11）のようになる。式（８），（11）より、式
（12）の演算を行なうことによつて式（９）が得られ
る。

Θ_Ｒ′＝θ_off＝θ_１′−（τ₁/τ_２）θ_２′ …（11） Ｔ（x,y）＝Ｒ（x,y）×exp（−ｊΘ_Ｒ′） …（12） 以上が分離画像を得るための画像処理の流れである。
以上の処理で問題になるのは、式（７）で示したＳ（x,
y）を作る操作である。すなわち、Ｅ（x,y）が小さく、
式（６）が実際に（13）を満足する場合には単純に式
（７）の計算でよいが、Ｅ（x,y） −π≦Θ_２（x,y）≦＋π …（13） が大きくなると、Θ_２（x,y）は、見かけ上、式（13）
を満足していても、実際には式（13）を満足しなくなつ
てくる。なぜなら、Θ_２（x,y）は式（14）で計算さ
れ、−π≦Θ_２（x,y）≦＋π、あるいは０≦Θ_２（x,
y）≦２π以外の値はとり得ないから Θ_２（x,y）＝tan^-1（Im（S₂（x,y）） /Re（S₂（x,y））） …（14） である。式（14）でImは（ ）内の虚数部、Reは（ ）
内の実数部をそれぞれ表わしている。静磁場不均一が大
きく、式（13）が満足されない場合には、S₂（x,y）の
位相像すなわちΘ_２（x,y）を画像表示すると隣り合つ
た座標点において−π→π，π→−πへの急激な位相変
化が表われているのがわかる。

このような場合のΘ_２（x,y）のある一次元の値を表
示すると第２図のようになる。図中で＊で示したところ
で＋π→−πへの変化が生じている式（13）が満足され
ない場合、すなわちΘ_２（x,y）が第２図のようになる
場合には、Ｅ（x,y）の空間的な変化がなだらかである
こと、つまり連続的に変化することを考慮し第３図のよ
うにΘ_２をΘ_２′に変換して位相連続性をとつてから、
式（７）の計算を行う。すなわち、式（７）のかわりに
式（15）のＳ（x,y）を用いる。

Ｓ（x,y）＝exp｛ｊ（τ₁/τ_２）Θ_２′（x,y）｝ …（15） 第２図，第３図は一次元での位相変化を示したもので
あるが、実際にはΘ_２（x,y）は二次元であるから、当
然のことながら二次元での位相連続性をとる計算処理が
必要となる。

この二次元での位相連続性をとる計算処理は、例えば
検査対象が人体のような複雑な構造の場合には計算処理
が複雑になるばかりでなく、良好に実行することができ
ない場合が多い。この計算処理を避けるためには、検査
対象を測定したのと同一のシーケンスにより、検査対象
領域よりも大きな断面をもつ一様フアントムを測定す
る。この一様フアントムが化学シフトσ_ｒをもつ物質の
みで満たされているとすると、第１エコー，第２エコー
より得られる画像U₁（x,y）,U₂（x,y）は式（５）に対
応して式（16）のようになる。

ここでρ_r,ρ_ｒ′（ρ_ｒ＞ρ_ｒ′）は化学シフトσ_ｒ
をもつ物質の密度、θ_r,θ_ｒ′は座標に依存しない一定
位相である。

U₁,U₂を用いて、式（５）のS₁,S₂に次のような演算を
行なう。S₂′（x,y）の位相項Θ_２′（x,y） は、式（18）のようになる。式（６）のΘ_２（x,y） Θ_２′（x,y）＝−γ〔E_per（x,y）＋（σ_１−
σ_ｒ）〕τ_２ ＋（θ_２′−θ_ｒ′）｝ …（18） と式（18）とを比較して、Ｅ（x,y）が、E_per（x,y）に
置き換つていることに注意されたい。

一般に静磁場不均一Ｅ（x,y）は、E_St（x,y）とE_per
（x,y）の和で表わされ、式（３）のような関係があ
る。従つてE_St（x,y）を除去したΘ_２′（x,y）はすべ
ての（x,y）においてΘ_２′にある一定値を加算あるい
は減算することにより、式（13）の条件を実際に満足す
ることが可能となる。すなわち、先に述べた位相連続性
をとる処理が不要となる。

S₁′（x,y）,S₂′（x,y）より分離画像を得る処理は
式（７）〜式（12）と同様であるので省略する。

先に述べたように、オフセツト位相θ_offを除去する
ためにはρ_１＝0,ρ_２≠0,ρ_２＝０あるいはρ_１≠0,ρ
_２＝０なる（x,y）における位相を調べる必要がある。
この操作を省略し、θ_offを除去するためには、θ_１＝
0,θ_２＝０とするタイミング、すなわち、通常のマルチ
エコーシーケンスにより得られた再構成画像を用いるこ
とができる。τ_１＝τ_２＝０つまりθ_１＝θ_２＝０とし
た場合、式（５）は式（19）に示すようになる。

このS₁₀,S₂₀を用いて、式（５）で示した、化学シフ
トイメージング用の画像に処理を行ない式（20）の画像
P₁,P₂を得る。

このP₁（x,y）,P₂（x,y）を用いて（６）式以降の処
理を行なえば、式（８）の下線部が１、となり、式（1
0）位降の処理なくして目的とする分離画像Ｔ（x,y）が
得られる。

以上の説明からわかるように、この方法でθ_offを除
去するためには化学シフトイメージング用データの他
に、通常のマルチエコーシーケンスによる測定が必要と
なり、測定時間が増大することになる。しかし、この通
常のマルチエコーシーケンスによる測定は、化学シフト
イメージング用データの測定に比して画像数を少なく、
具体的にはエンコード数を少なくして測定を行なうこと
が可能であり、通常のマルチエコーシーケンスによる、
検査対象に対する測定時間全体の増加を不必要に大きく
することはない。

この通常のマルチエコーシーケンスによるオフセツト
位相の除去法は前記の一様フアントムを用いる方法につ
いても同様に適用できる。この場合には検査対象だけで
なく一様フアントムについても通常のマルチエコーシー
ケンスによる測定を行ない、式（５）のθ_１′，
θ_２′、式（16）のθ_r,θ_ｒ′を除去した上で式（17）
以降を実行すればよい。

以上では、マルチエコーシーケンスにより第1,第２エ
コーを取得し、分離画像を得る方法について示したが、
マルチエコーシーケンスによらず、通常のスピンエコー
シーケンス、すなわち、第１図の第１エコーまでのシー
ケンスで化学シフトイメージングを実現するためには、
測定を２回実施する必要がある。つまり、１回目の測定
で第１図のτ_１に相当するタイミングで信号計測を行な
い、２回目の測定でτ_２に相当するタイミングで信号計
測を行なう。この際、２回目の測定では、フアントムを
使う場合、あるいは、フアントムを使用しない場合でも
静磁場不均一の空間的な変化の割合が小さければ１回目
の測定と比較して画素数を少なく、具体的にはエンコー
ド数を少なくして測定全体に要する時間を短縮すること
が可能である。

マルチエコーシーケンスでも、第２エコーを測定する
際にエンコード数を少なくすることが可能であるが、測
定時間の短縮にはあまり効果がない。

マルチエコーシーケンスによる測定では、一般に第１
エコーより得られる画像に比べて第２エコーより得られ
る画像のS/Nが低下するという欠点があるが、フアント
ムを使う場合、あるいは、フアントムを使用しない場合
でも第２エコー信号に適当なフイルター処理を施すこと
によつてS/N向上を計ることが可能である。

〔発明の効果〕

本発明によれば、検査対象に対するマルチエコーシー
ケンスによるほぼ１回の測定時間で化学シフトイメージ
ングによる分離画像を得ることができるので、検査対象
からの信号測定に要する時間を短縮でき、良好な分離画
像を得ることができるという効果がある。

また、一様フアントムを検査対象と同一のシーケンス
で測定し、得られた画像を利用することにより、二次元
での位相連続性をとる処理が不要となるので、分離画像
を得るための画像処理が容易になるという効果がある。

【図面の簡単な説明】

第１図は本発明の一実施例である。変形マルチエコーシ
ーケンスを示す図、第２図は位相連続性をとる処理が必
要な場合の一次元での位相変化を示す図、第３図は第２
図より位相連続性をとる処理を行なつた場合の一次元で
の位相変化を示した図である。

フロントページの続き (58)調査した分野(Int.Cl.⁶，ＤＢ名) A61B 5/055

Claims (7)

(57)【特許請求の範囲】
1. 【請求項１】静磁場、傾斜磁場、及び高周波磁場の各磁
   場を発生する手段と、前記静磁場が印加された空間に置
   かれた検査対象からの核磁気共鳴信号を検出する信号検
   出手段と、該信号検出手段の検出信号の演算処理及びパ
   ルスシーケンスの制御を行なう信号処理制御手段とを有
   し、前記検査対象の核スピンの分布情報を得る核磁気共
   鳴を用いた検査装置において、前記信号処理制御手段
   は、 （１）（ｉ）前記検査対象の核スピンを励起すること； （ii）エンコード傾斜磁場を印加すること； （iii）２つのケミカルシフトをもつ核スピンの間で、
   所定の位相差θ_１がθ_１≠ｎπ（ｎ＝０、±１、±２、
   …）となる時間差をもたせたタイミングで、第１の180
   ゜高周波磁場パルスを、前記検査対象に印加すること； （iv）リードアウト傾斜磁場を印加し、第１の核磁気共
   鳴信号を計測すること； （ｖ）前記２つのケミカルシフトをもつ核スピンの間
   で、位相差θ_２がθ_２＝2nπ（ｎ＝±１、±２、…）と
   なる時間差をもたせたタイミングで、第２の180゜高周
   波磁場パルスを、前記検査対象に印加すること； （vi）前記リードアウト傾斜磁場を印加し、第２の核磁
   気共鳴信号を計測すること；及び （vii）位相エンコード量を変化させて前記エンコード
   傾斜磁場を印加し、（ｉ）から（vi）を繰返すこと；か
   らなり、前記第１、第２の核磁気共鳴信号の計測を行な
   い、第１、第２のデータセットを得るための前記核磁気
   共鳴信号を計測するパルスシーケンス制御と、 （２）（ａ）前記第１、第２のデータセットをそれぞれ
   フーリエ変換して、第１、第２の画像データを求めるこ
   と； （ｂ）前記２つのケミカルシフトのうちの１つをもつフ
   ァントムを、前記静磁場が印加された前記空間に置い
   て、前記の（１）と同じ前記パルスシーケンス制御を実
   行して、第３、第４のデータセットを得て、前記第３、
   第４のデータセットをそれぞれフーリエ変換して、第
   ３、第４の画像データを求めること； （ｃ）前記第３、第４の画像データの位相分布をそれぞ
   れ補正データとして使用し、前記第１、第２の画像デー
   タからそれぞれ前記静磁場の不均一に起因する位相誤差
   分布を除去した第５、第６の画像データを求めること；
   及び （ｄ）前記第６の画像データから計算される位相分布を
   補正データとして使用し、前記第５の画像データから、
   前記検査対象の透磁率分布に起因する位相誤差分布を除
   去することにより、前記２つのケミカルシフトの核スピ
   ンの分布情報を含む第７の画像データを求めること；を
   含む演算処理を行なうことを特徴とする核磁気共鳴を用
   いた検査装置。
2. 【請求項２】請求項１に記載の核磁気共鳴を用いた検査
   装置において、 前記ファントムは前記検査対象の領域よりも大きい断面
   を有し、前記２つのケミカルシフトのうちの１つのみか
   らなる均一な物質を含むことを特徴とする核磁気共鳴を
   用いた検査装置。
3. 【請求項３】請求項１に記載の核磁気共鳴を用いた検査
   装置において、 前記演算処理は、前記第７の画像データのピクセル位置
   に依存しないオフセット位相を除去することを、さらに
   含むことを特徴とする核磁気共鳴を用いた検査装置。
4. 【請求項４】請求項３に記載の核磁気共鳴を用いた検査
   装置において、 前記オフセット位相は、前記第７の画像データのピクセ
   ルのうち、前記２つのケミカルシフトのうちの１つのみ
   が存在するピクセルのデータの位相から検出されること
   を特徴とする核磁気共鳴を用いた検査装置。
5. 【請求項５】請求項２に記載の核磁気共鳴を用いた検査
   装置において、 前記演算処理は、前記第７の画像データのピクセル位置
   に依存しないオフセット位相を除去することを、さらに
   含むことを特徴とする核磁気共鳴を用いた検査装置。
6. 【請求項６】請求項５に記載の核磁気共鳴を用いた検査
   装置において、 前記オフセット位相は、前記第７の画像データのピクセ
   ルのうち、前記２つのケミカルシフトのうちの１つのみ
   が存在するピクセルのデータの位相から検出されること
   を特徴とする核磁気共鳴を用いた検査装置。
7. 【請求項７】静磁場、傾斜磁場、及び高周波磁場の各磁
   場を発生する手段と、前記静磁場が印加された空間に置
   かれた検査対象からの核磁気共鳴信号を検出する信号検
   出手段と、該信号検出手段の検出信号の演算処理及びパ
   ルスシーケンスの制御を行なう信号処理制御手段とを有
   し、前記検査対象の核スピンの分布情報を得る核磁気共
   鳴を用いた検査装置において、前記信号処理制御手段
   は、 （１）（ｉ）前記検査対象の核スピンを励起すること； （ii）エンコード傾斜磁場を印加すること； （iii）第１、第２の180゜高周波磁場パルスを、前記検
   査対象に順次印加すること； （iv）２つのケミカルシフトをもつ核スピンの間で、所
   定の位相差θ_１がθ_１≠ｎπ（ｎ＝０、±１、±２、
   …）となる時間差をもたせたタイミングで、リードアウ
   ト傾斜磁場を印加し、第１の核磁気共鳴信号を計測する
   こと； （ｖ）前記２つのケミカルシフトをもつ核スピンの間
   で、位相差θ_２がθ_２＝2nπ（ｎ＝±１、±２、…）と
   なる時間差をもたせたタイミングで、前記リードアウト
   傾斜磁場を印加し、第２の核磁気共鳴信号を計測するこ
   と；及び （vi）位相エンコード量を変化させて前記エンコード傾
   斜磁場を印加し、（ｉ）から（ｖ）を繰返すこと；から
   なり、前記第１、第２の核磁気共鳴信号の計測を行な
   い、第１、第２のデータセットを得るための前記核磁気
   共鳴信号を計測するパルスシーケンス制御と、 （２）（ａ）前記第１、第２のデータセットをそれぞれ
   フーリエ変換して、第１、第２の画像データを求めるこ
   と； （ｂ）前記２つのケミカルシフトのうちの１つをもつフ
   ァントムを、前記静磁場が印加された前記空間に置い
   て、前記の（１）と同じ前記パルスシーケンス制御を実
   行して、第３、第４のデータセットを得て、前記第３、
   第４のデータセットをそれぞれフーリエ変換して、第
   ３、第４の画像データを求めること； （ｃ）前記第３、第４の画像データの位相分布をそれぞ
   れ補正データとして使用し、前記第１、第２の画像デー
   タからそれぞれ前記静磁場の不均一に起因する位相誤差
   分布を除去した第５、第６の画像データを求めること；
   及び （ｄ）前記第６の画像データから計算される位相分布を
   補正データとして使用し、前記第５の画像データから、
   前記検査対象の透磁率分布に起因する位相誤差分布を除
   去することにより、前記２つのケミカルシフトの核スピ
   ンの分布情報を含む第７の画像データを求めること；を
   含む演算処理を行なうことを特徴とする核磁気共鳴を用
   いた検査装置。