JPH0241137A - ケミカルシフトイメージング装置 - Google Patents

Abstract

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

Description

【発明の詳細な説明】 （産業上の利用分野） 本発明は、ケミカルシフトを利用して体内の水と脂肪の
分離像を得る際の、静磁場平均−などに起因した位相エ
ラーを補正する方法に関する。更に詳しくは、この補正
における位相飛びエラーを除去したケミカルシフトイメ
ージングにおける位相補正方法に関する。

（従来の技術） Ｍ　ＲＩ　（Ｍａｇｎｅｔｉｃ　Ｒｅ５ｏｎａｎｃｅ　
ｌｌ１ｌａＱｉｎ（１）において、２成分の分子構造の
違いにより同一核種の共鳴周波数がずれるケミカルシフ
トを利用して、体内の同一断層面を、水のみのプロトン
画像と脂肪のみのプロトン画像に分離して表示する分離
イメージングがある。

始めに、水と脂肪の分離イメージング方法の１つである
従来例のデイクソン法を説明する。第５図は分離イメー
ジングのパルスシーケンスを表す図である。第５図にお
いて、ｔは時間軸であり、ＲＦは静磁場方向（２方向）
に向いたプロトンの磁化ベクトルを励起するＲ　Ｆ　（
Ｒａｄｉｏ−Ｆｒｅｑｕｅｎｃｙ　）波であり、回転角
に応じて９０°パルス、１８０゛パルスと呼ばれる。　
Ｇ５１ｉｃｅは９０’パルスと同時に印加され、断層面
の選択励起を行うスライス勾配、ＳＥは９０°パルス印
加後、静磁場方向に垂直な面（ＸＹ平面）内に、静磁場
平均−及び勾配磁場によってばらけな磁化ベクトルの位
相を、１８０°パルスで反転して再び収束した時に観測
されるスピンエコー信号（以後ＳＥ信号と呼ぶ）である
、　Ｇｐｈａｓｅは位相エンコード勾配磁場、Ｇｒｅａ
ｄは周波数エンコード勾配磁場であり、例えば各々Ｙ方
向を位相、Ｘ方向を周波数情報に対応させてＳＥ信号に
２次元の位置情報を与える。デイクソン法では、１８０
°パルスの印加するタイミングの異なる２種類のパルス
シーケンスを用いる。即ちＳ１スキャンでは１８０°パ
ルスの印加タイミングτを、９０゜パルスからスピンエ
コーが得られるまでの時間ＴＥの中間のτ＝ＴＥ／２と
し、Ｓ２スキャンではＳ１スキャン時よりε早いτ＝Ｔ
Ｅ／２−εとする。

Ｓ２スキャンもスピンエコーを得るタイミングは、Ｓ１
スキャンと同じＴ６とする。

ここで、前記εは次式を満たす。

ε＝Ｔ７７τ了　　　　　　　・・ｆ）σ：水と脂肪の
ケミカルシフト量 ｆ：プロトンの共鳴周波数 Ｓ１スキャンにおいては、時刻ＴＥ／２の１８０　”　
パルスによって位相が反転されるため、時刻Ｔ［で水と
脂肪の磁化ベクトルの位相が一致する。３２スキャンに
おいては、ε早＜１８０°パルスによって位相が反転さ
れるため、脂肪の磁化ベクトルが２ε時間分だけ余計に
回転し、ｔ＝ＴＥで水と脂肪の磁化ベクトルが１８０°
の位相ずれを生じる。

従って、各スキャンによって得られるローデータを画像
再構成した画像データは、 Ｓｌ　＝Ｗ＋Ｆ　　　　　　　　　　　　・・・（２）
Ｓｌ　＝Ｗ−Ｆ　　　　　　　　　　　　・・・（３）
になる。ここで、Ｗ（≧０）は水のプロトン密度、Ｆ（
≧０）は脂肪のプロトン密度である。そこで、Ｗ＝　（
Ｓｌ　＋Ｓ２　）／２　　　　　　　・・・（４）Ｆ＝
　（３１−Ｓ２　＞／２　　　　　　　・・・（５）を
各ピクセル毎に計算することによって、水と脂肪を分離
した画像を得ることができる。

しかし、実際の８１スキャンデータには、静磁場不均一
等による位相のずれ成分が加わる。このような位相ずれ
成分があると、（４）、　（５）式での水と脂肪成分の
分離ができず、シェーディング等の分離エラーが生じる
。この静磁場不均一等による位相のずれには、以下のよ
うな要因がある。まず、静磁場の均一性はマグネットに
起因する不均一に加え、実際の撮影では人体が静磁場内
に入り、このため空気と人体各部位との透は率が異なり
、又、反磁界の影響で磁場不均一の状態が異なり、均一
・度が変化し、悪化する。このため、撮影面内の各ピク
セル毎にラーモア周波数がずれることにより位相ずれが
生じる。又、周波数エンコード勾配磁場印加時にマグネ
ットボア材等に誘起される渦電流によって、矩形の勾配
磁場Ｇｒｅａｄの立ち下りが第５図点線部のように鈍り
、残り勾配が生じる。

このため、１８０°パルス印加タイミングの異なる２ス
キャンでのＳＥ倍信号エコーセンターが点線のようにδ
すれ、その分データサンプリング点がずれるため、各ス
キャンデータには、更に１次位相エラー成分が含まれる
。又、中心位相のずれ等も１次位相エラー成分の原因と
なる。　そこで、第５図のように、　１８０°　パルス
の印加タイミングτを、Ｓ１スキャン時よりε遅いτ＝
ＴＥ／２±εとしなＳ３スキャンを行い、Ｓ２スキャン
と８３スキャンデータを２次元フーリエ変換して得なデ
ータにより静磁場不均一分布に関するデータマトリック
スを得て、これから８２又はＳ３のデータを位相補正し
、この補正されたデータ及びＳ１スキャンからのデータ
を用いて水と脂肪の分離像を得る方法が提案されている
。

第６図（ａ）、　（ｂ）、　［Ｃ）は３１．３２．３３
スキャンにおける水と脂肪の磁化ベクトルの位相関係を
表す図である。第６図において、Ｒｅａｌ−ＩＩＩｌａ
ａｉｎａｒｙ軸による座標系は、水の磁化ベクトルの位
相を基準にした位相平面を表す。前記のように、Ｓ１ス
キャンにおいては、ｔ＝ＴＥで水と脂肪の磁化ベクトル
の位相が一致し１、Ｓ２及びＳ３スキャンにおいては、
水と脂肪の磁化ベクトルが１８０°の位相ずれを生じる
。更に、Ｓ２及びＳ３スキャンでは、静磁場不均一等に
よる位相ずれ成分が各４士θ、−θだけ加わる。そこで
、Ｓ２スキャンと８３スキャンデータを２次元フーリエ
変換して得たデータにより各ピクセル毎にθを求めて、
この静磁場不均一分布によりＳ２又はＳ３のデータを位
相補正することにより、位相ずれに起因した分離エラー
を除去することができる。

（発明が解決しようとする課題） しかし、上記のようなケミカルシフトイメージングにお
ける位相補正方法においては以下の問題点がある。前記
不均一分布の位相値の検出範囲は、−πからπまでであ
るため、この範囲を越えると位相が折り返る位相飛びが
生じる。第７図（ａ）、（ｂ）は静磁場不均一による位
相分布く以ｆ＆静磁場不均−分布と呼ぶ）を表す図であ
る。第７図（ａ）は２５６ｘ２５６マトリツクスの静磁
場不均一分布を表す図であり、この分布図の１位相エン
コードラインのデータを周波数エンコード方向にプロッ
トしたグラフが第７図（ｂ）である。図のように、ノイ
ズピクセルや１次位相エラーがあると、位相値が±πを
越える位相飛びが生じる。この位相飛びが生じると２成
分の分離エラーが起こる。

本発明は、前記問題点を解消し、静磁場不均一分布上の
位相飛びを最少限に抑え、分離エラーを低減したケミカ
ルシフトイメージングにおける位相補正方法を提供する
ことにある。

（課題を解決するための手段） ９０”のＲＦパルス印加がら１８０°のＲＦパルス印加
までの時間τをＴＥ／２としたＳ１スキャンと、前記時
間τをＴＥ／２−εとしたＳ２スキャンと、前記時間τ
をＴＥ／２＋εとしたＳ３スキャンを行い、Ｓ２スキャ
ン及びＳ３スキャンからのデータにより静磁場不均一分
布を求め、位相ずれ防止のための位相補正を行うゲミカ
ルシフトイメージングにおける位相補正方法において、
前記静磁場不均一分布に対して位相飛び防止のための位
相補正を施し、この位相補正後の静磁場不均一分布によ
りＳｌ又はＳ３スキャンからのデータに対して位相ずれ
防止のための位相補正を行い、この位相補正されたＳｌ
又はＳ３スキャンのデータに対して、前記静磁場不均一
分布に対する位相補正分を補うための位相補正を行うこ
とを特徴とする。

（作用） まず、静磁場不均一分布に対して位相値を±π内に抑え
るために位相補正を行い、この補正後の不向−分を８２
又はＳ３スキャンからのデータへ位相補正し、本来静磁
場不均一分布に含まれていた前記補正分の位相補正を８
２又はＳ３スキャンからのデータに対して行うため、位
相飛びに起因した分離エラーを除去できると共に、位相
補正が正確になる。

（実施例） 以下、図面を参照して本発明について詳細に説明する。

スキャンシーケンスは、デイクソン法を改良した第５図
の３１．　Ｓ　２．　Ｓ　３スキャンを用いる。

第１図は本願発明の一実施例の位相補正方法を表すフロ
ーチャートである。以下、第１図を用いて、本実施例の
位相補正方法を説明する。始めに、Ｓｌ、　Ｓ　２．　
Ｓ　３スキャンで得られた複素数のローデータを２次元
フーリエ変換する。各データＳｌ、８２゜Ｓ３は、次式
で表される。

Ｓｌ　＝Ｗ＋Ｆ　　　　　　　　　　　　・・・（６）
Ｓｌ　＝　（Ｗ＋Ｆ−ｅｘｐ　　（ｉ　π）　１１３Ｘ
Ｅｌ　　（ｉθ）・・（７） Ｓ３　＝　（Ｗ＋Ｆ’−ｅＸＤ　　（ｉ　π）　）ｅｘ
ｐ（−ｉθ）・・・（８） ここで、１２＝−１であり、ｅｘｐ（ｉθ）は、静磁場
不均一等による位相ずれ成分である。

次に、（７）（８）式から各ピクセル毎に、静磁場不均
一分布θｏ（ｒ、ｗ）　（１≦「、ｖ４≦２５６）を次
式に従って求める。

θｏ（ｒ、ｗ）＝ｔａｎ　−”　５２）−（９）Ｔ丁丁
丁 ここで、ｒは周波数エンコード軸、Ｗは位相エンコード
軸である。

次に、ローデータからノイズレベルＮｔ、ｃｖ−を設定
する。第２図はデータ平面上のローデータを表す図であ
る。第２図における座標は、横軸が時間、縦軸が位相エ
ンコード量に対応している０位相エンコード量の大きい
所のローデータは、はとんどがノイズであることからこ
の領域のデータを用い、周波数エンコード方向に１次元
逆フーリエ変換を施し、その結果を絶対値処理したデー
タの平均値Ｍｅａｎと分散σ０を求める。前記ノイズレ
ベルＮＬｉｖｇｔは Ｎ　Ｌｌ！ＶＥＬ　　＝　Ｈｅａｎ＋　　３．５σＯ−
（１０）とする、σ０の係数を３．５とすることで、ノ
イズの９９．８％がＮＬｌ、ＶＩＩＬに含まれる。

このノイズレベルＮＬＥＶＥＬと各ピクセルのイメージ
データの絶対値とを比較して、ノイズレベルＮＬＥＶＥ
Ｌ以上を１、ノイズレベルＮＬＥＶＩ！Ｌ未満を０とし
たマスクファイルＶ（ｒ、ｗ）　（１≦ｒ、ｗ≦２５６
）を作成する。このマスクファイルは以後の位相係数算
出において、ノイズピクセルを排除するために必要にな
る０次に、このマスクファイルを用いて、ノイズとイメ
ージの境界を繋いでイメージエリヤを設定する。第３図
は、このようにして設定したイメージエリヤを表す図で
ある。第３図において、ｒは周波数エンコード方向、Ｗ
は位相エンコード方向を表□し、ｒｏ、ｗｏは各々周波
数エンコード方向及び位相エンコード方向のイメージの
中心座標である。次に、イメージエリヤ内のイメージ中
心付近の平均位相α０を次式に従って算出する。

ここで、Ａは平均位相αＯ算出に用いた有効ピクセル数
である。この平均位相α０を用いて、次式に従って静磁
場不均一分布θｏ（ｒ、ｖ）に対して位相補正する。

θ１（ｒ、ｗ）＝θｏ（ｒ、ｗ）ｘ　ｅｘｐ（−ｉαｏ
）　　　　　−（１２）通常のイメージングでは、最も
感心のある部位をイメージ中心となるようにスキャンを
設定するため、この処理によりイメージ中心での位相を
０とする、ここで、平均位相α０算出に用いた領域内に
位相飛びによる位相段差かある恐れがあるために、以下
の処理を行う、始めに、イメージ中心付近の分散σを次
式に従って算出する。

σ　＝ ・・・　（１３） ここで、平均位相α０が０以上のときは、α１＝α０＋
σ、０未満のときはα１−αＯ−σとして、（１２）式
に従ってθｏ（ｒ、ｗ）に対して位相補正したθ１ａ（
ｒ、ｖ）を用い、再び（１１）、　（１３）式から平均
位相α２及び分散σ２を求める。そして、分散σ２が０
．８５σより大きいときは位相段差がないと見なし、（
１２）式の補正を有効とし、分散σ２か０．８５σ以下
のときは、（１２）式の補正を無効とし、α２を平均位
相として、次式に従って○次位相補正し直す。

θＢｒ、Ｗ）＝θＯ（ｒ、Ｗ）Ｘ　ｅｘｐ（−ｉα２）
　　　−（１４）次に、周波数エンコード方向での一次
位相係数を求める。第４図は、本実施例の静磁場不均一
分布に対する位相補正方法を表す図である。第４図（ａ
）、　（ｂ）、　（ｃ）において、各グラフは静磁場不
均一による１位相エンコードラインの位相を周波数エン
コード方向にプロットしたものである。第４図（ａ）の
ように、渦電流に起因したエコーセンターのずれ等によ
り、不均一分布上の周波数エンコード量にほぼ比例して
大きくなる１次の位相ずれが生じている。この位相ずれ
による位相飛びを減らすために、１次位相係数を求めて
補正する。まず、周波数エンコード方向の隣接ピクセル
間の位相差を次式に従って求める。

Δ　θ　ｆｒ、Ｗ）　　＝　　（θ　１（ｒ＋１．ｗ）
−θ　１（ｒ、ｗ））ｘ　Ｖ　　ｆｒ＋１．ｗ）　　ｘ
　　Ｖ　（ｒ、ｗ）（ｒ＝１．３，５．・・・、２５５
．　Ｗ＝１．２，３．・・・、　２５６　）・・ｆ５）
ここで、前記隣接側ピクセルの内どちらか一方がノイズ
ピクセルの位相差は、マスクファイルにより０として除
外される。ここで、同時に位相飛びレベルθＬＥＶＥＬ
を設定し、位相飛びピクセルが関与した位相差を除外す
る。例えば、位相飛びレベルθＬｌ！Ｖｌ！Ｌを１．５
πとすると、隣接ピクセル間の位相差が１．５πを越え
るものは、位相飛び部分として除外する。このようにし
て算出された位相差に基つき、全位相差の加算平均値を
求め、これを１次位相係数α１とし、次式に従って不均
一分布に対して補正する。

θ２（ｒ、ｖ）＝θ１（ｒ、ｗ）　ｘ　ｅ×ｐ’（ｆ−
ｉ　α１）・（ｒ−ｒｏ））・・ｆ６） このようにして、１次の位相補正が行われた結果が第４
図（ｂ）である６以上の補正により不均一分布上で低次
（０，１次）の位相ずれ成分による位相飛び領域が抑え
られる。

ここで、更に不均一分布上での高次の位相ずれ成分によ
る位相飛び領域を抑えるために、高次関数による位相補
正を、次式に従って各位相エンコードライン毎に周波数
エンコード方向に対して行う。

θ３（ｒ　ｗ）＝０２　（ｒ、　Ｗ）　Ｘ　ｅｘｐ（（
−ｉΣγにｗ（ｒ−ｒｏ）Ｋ）・・・　（１７） ここで、γに１はＷ番目の位相エンコードラインでのに
次の係数であり、前記１次係数算出時と同様に周波数エ
ンコード方向のに階の位相差の１７に１から求めること
ができる。

更に、次式に従って各周波数エンコードライン毎に位相
エンコード方向に対して高次関数による位相補正を行う
。

θ４（ｒ、ｗ）＝０３（ｒ、ｖ）ｘ　ｅｘｐ（［−ｉ　
　Σβ、ｒ（ＩＩＩ−１ＩＩＯ）Ｋ）・・・　（１８） ここで、βｋｒはＷ番目の周波数エンコードラインでの
に次の係数であり、前記と同様に位相エンコード方向の
に階の位相差の１７に！から求めることができる。この
ようにして、例えば２次の位相補正が行われた結果が第
４図（Ｃ）である。この補正により、更に不均一分布上
で高次の位相すれ成分による位相飛び領域が抑えられる
。即ち、（１７）、　（１８）式による補正は不均一成
分を高次の多項式でフィッティングすることに相当して
いる。通常のＭ　Ｒ■装置では、シムコイルなどにより
低次の静磁場不均一成分は補正しであるが、人体がマグ
ネット内に入ることにより静磁場不均一分布が変化する
。

従って、低次及び高次のフィッティングを、不均−分布
に対して行うことにより、人体がある場合の不均一状態
を概略近似することができ、位相飛びを防ぐことができ
る。

次に、このようにして位相飛びを抑えるために補正され
た不均一分布θ４（ｒ、ｗ）を用いて、Ｓ２スキャンか
らのデータに対して、次式に従って位相補正する。

Ｓ　２ａ＝　Ｓ　２　ｘ　ｅｘｔ）（−ｉθ４（ｒ、ｗ
）／２）　　　−（１９）次に、位相飛びを抑えるため
に不均一分布θ（「、Ｗ）を補正した位相量は、本来に
８２スキャンからのデータに対して行われるべきもので
ある。従って、前記補正された位相量分をＳ２ａデータ
に対して、次式に従って順次補正する。

５２ｂ＝　５２ａｘｅｘｐ（−ｉαｏ／２）　　　　　
　　−（２０）（但し、（１４）式を用いたときはαＯ
→α２とする）Ｓ　２ｃ＝　Ｓ　２ｂｘ　ｅｘｐ（（−
ｉ　ａ　１）　−（ｒ−ｒｏ）／２）・・・　（２１） ・・・　（２２） Ｓ　２ｅ＝　Ｓ　２ｄｘ　ｅｘｐ（（−ｉ　　Σβｋｒ
（Ｗ−ＷＯ）”　　／２１・・・　（２３） 以上で８２データに対する位相補正を終了し、最後に、
このようにして位相補正されたＳ２データ及びＳ１デー
タを用いて、次式に従って、複素演算及び絶対値処理を
行い水と脂肪の分離像を得る。

Ｗａｔｅｒ−ｌ　１／２（Ｓ　１　＋　５２ｅ）Ｆａｔ
　　−ｌ　１／２（Ｓｌ　−３２ｅ）　　ｌ　　　　　
−（２４）上記のように、本実施例のゲミカルシフトイ
メージングにおける位相補正方法においては、静磁場不
均一分布θの位相飛びを抑えるために、０．１、高次の
関数による位相補正を順次行い、この位相補正された静
磁場不均一分布により１８０°パルスをずらしたＳ２ス
キャンデータに対して、前記関数から外れた分の不均一
分布による位相補正を施し、その後Ｓ２スキャンデータ
に対して、前記関数分の位相補正を行うので、位相飛び
によるエラーを極力抑えることができ、分離が正確にな
る。

尚、本発明は上記実施例に限定するものではなく、特許
請求の範囲内で種々の変形が可能である。

本実施例では高次関数のフィッティングを位相エンコー
ド又は周波数エンコードの各ライン毎に行ったが、不均
一分布全体に対して、次式による２変数の多項式からな
る関数によるフィッティングを行って補正しても良い。

（ａｍ　＋　ｂ＊　＋　Ｃｋ、ｄｔ、”・定数）　　　
−＜２４＞更に、Ｓ２スキャンデータに対して位相補正
を行い、（２３）式の分離に用いたが、Ｓ３スキャンデ
ータを用いても良い、更に、（１５）式による１次位相
係数算出は、各位相エンコードライン毎に行わず、位相
エンコードライン方向の平均を取り、この平均ラインに
対して周波数エンコード方向に１次係数を求めて用いて
も良い。

（発明の効果） 以上の説明の通り、本発明のゲミカルシフトイメージン
グにおける位相補正方法によれば、以下の効果が得られ
る。

（１）Ｓ２又はＳ３スキャンデータがら静磁場不均一分
布を求め、位相飛び防止のための位相補正された静磁場
不均一分布により、Ｓ２又はＳ３スキャンデータに対し
て位相ずれ防止のための位相補正を施し、その後このス
キャンデータに対して前記静磁場不均一分布への位相補
正分を補うための位相補正を行うので、前記分布上の位
相飛びが最少限に抑えられ、この分布による補正時のエ
ラーを防止し、なおかつ、Ｓ２又はＳ３スキャンデータ
に含まれる位相ずれに起因した分離エラーを除去できる
。又、Ｓ２スキャン及び３３４ヤンデータによって得ら
れた、対象とする人体内での静磁場不均一情報を余すと
ころなく利用して補正できるため、位相補正が正確にな
る。

（２）静磁場不均一分布に対して０．１次の関数による
位相補正を順次行う、従って、静磁場不均一成分を含む
中心周波数のズレ及び渦電流に起因しな１次位相エラー
を不均一分布から取り除くことができるため、不均一分
布での位相飛び領域を低減することができる。更に、静
磁場不均一分布に対して高次の関数による位相補正を順
次行う。従って、人体が静磁場内に入ることによる生じ
る静磁場不均一に基因した位相エラーを不均一分布から
収り除くことができるため、不均一分布での位相飛び領
域を更に低減することができる。これにより、位相飛び
に起因した分離エラーを除去できる。

【図面の簡単な説明】

第１図は本願発明の一実施例の位相補正方法を表すフロ
ーチャート、第２図は本願発明の一実施例のデータ平面
上のローデータを表す図、第３図は木琴発明の一実施例
のイメージエリヤの設定を表す図、第４図（ａ）、　（
ｂ）、　（ｃ）は本願発明の一実施例の静磁場不均一分
布に対する位相補正方法を表す図、第５図は従来例の分
離イメージングのパルスシーゲンスを表す図、第６図（
ａ）、　（ｂ）、　（ｃ）は従来例の３１．Ｓ２　、Ｓ
３スキャンにおける水と脂肪の磁化ベクトルの位相関係
を表す図、第７図（ａ）（ｂ）は従来例の静磁場不均一
による位相分布を表す図である。 第２図 第３図 （Ｇ） 第６図 （ｂ） （Ｃ） 第７図 （Ｇ） （ｂ）

Claims (7)

【特許請求の範囲】
1. （１）Ｚ軸方向に略一様な静磁場を発生させ、この静磁
   場中の被検体にＺ軸に垂直な方向から９０°及び１８０
   °の周波数ｆのＲＦパルスを印加し、前記静磁場に周波
   数エンコード勾配磁場と位相エンコード勾配磁場を重畳
   し、９０°のＲＦパルスからＴＥ秒後に被検体からのス
   ピンエコー信号を受信し、このスピンエコー信号からの
   スキャンデータを２次元フーリエ変換したデータを用い
   て水と脂肪の分離像を得る方法であって、前記９０°の
   ＲＦパルス印加から１８０°のＲＦパルス印加までの時
   間τをＴＦ／２としたＳ１スキャンと、前記時間τをＴ
   Ｅ／２−εとしたＳ２スキャンと、前記時間τをＴＥ／
   ２＋εとしたＳ３スキャンを行い、Ｓ２スキャン及びＳ
   ３スキャンからのデータにより静磁場不均一分布を求め
   、この静磁場不均一による位相ずれの補正を行うケミカ
   ルシフトイメージングにおける位相補正方法において、
   前記静磁場不均一分布に対して位相飛び防止のための位
   相補正を施し、この位相補正後の静磁場不均一分布によ
   りＳ２又はＳ３スキャンからのデータに対して位相ずれ
   防止のための位相補正を行い、この位相補正されたＳ２
   又はＳ３スキャンのデータに対して、前記静磁場不均一
   分布に対する位相補正分を補うための位相補正を行うこ
   とを特徴とするケミカルシフトイメージングにおける位
   相補正方法。 但し、前記εは次式を満たす。 ε＝１／（４・σ・ｆ）…（１） σ：水と脂肪のケミカルシフト量
2. （２）前記静磁場不均一分布に対する位相補正は、ノイ
   ズレベルを設定し、ノイズレベル以上のノイズピクセル
   を除去する過程を含むことを特徴とする請求項（１）記
   載のケミカルシフトイメージングにおける位相補正方法
   。
3. （３）前記静磁場不均一分布に対する位相補正は、前記
   静磁場不均一分布のイメージ中心付近の平均位相により
   ０次位相補正を施す過程を含むことを特徴とする請求項
   （１）記載のケミカルシフトイメージングにおける位相
   補正方法。
4. （４）前記静磁場不均一分布に対する位相補正は、周波
   数エンコード方向の位相データから１次係数を求め、１
   次位相補正を施す過程を含むことを特徴とする請求項（
   １）記載のケミカルシフトイメージングにおける位相補
   正方法。
5. （５）前記１次係数は、周波数エンコード方向の隣接ピ
   クセル間の位相差平均から求めることを特徴とする請求
   項（４）記載のケミカルシフトイメージングにおける位
   相補正方法。
6. （６）前記静磁場不均一分布に対する位相補正は、周波
   数エンコード方向又は位相エンコード方向に高次関数に
   よる位相補正を施す過程を含むことを特徴とする請求項
   （１）記載のケミカルシフトイメージングにおける位相
   補正方法。
7. （７）前記静磁場不均一分布に対する位相補正は、２変
   数の多項式からなる関数を用いて、前記静磁場不均一分
   布全体に対してフィッティングを行い、この関数を用い
   てＳ２又はＳ３スキャンのデータに位相補正を行う過程
   を含むことを特徴とする請求項（１）記載のケミカルシ
   フトイメージングにおける位相補正方法。