WO2002082994A1

WO2002082994A1 - Magnetic resonance imaging apparatus and image processing method

Info

Publication number: WO2002082994A1
Application number: PCT/JP2002/003614
Authority: WO
Inventors: Yumiko Yatsui; Tetsuhiko Takahashi
Original assignee: Hitachi Medical Corp
Current assignee: Hitachi Healthcare Manufacturing Ltd
Priority date: 2001-04-11
Filing date: 2002-04-11
Publication date: 2002-10-24
Anticipated expiration: 2003-10-11
Also published as: JP2002306445A; US8027713B2; US20040135577A1; JP3976515B2; EP1380257B1; EP1380257A4; EP1380257A1

Description

明細書磁気共嗚ィメ一ジング装置及び画像処理方法技術分野

本発明は、磁気共鳴イメージング装置（以下、 MRI装置という）に係り、特に 1つの画像中に 2つ以上の島状の撮影領域を含んでいる場合、各領域に対して正確に水 ·脂肪分離画像を取得する技術に関する。背景技術

磁気共鳴イメージング装置で画像を得る場合、エコー時間（TE) や繰り返し時間（TR)等のパラメータを変化させたり、画像演算を行うことによって、様々な組織コントラストを持つ画像を得ることができる。そこで、コントラストを持つ画像を得る方法の 1つに Dixon法がある。 Dixon法とは、ある 1つ領域において、水信号の MR信号と脂肪の MR信号の位相の変ィヒを利用して水 ·脂肪分離画像を得る方法である。

近年、臨床において、大腿部を左右同時に撮影したり、頭部脳底部を横断面で撮影したり、腹部を冠状断面で撮影し、水'脂肪分離画像で診断する機会が増えている。このように空間的に離れた複数の部分からなる被検体や、被検体中に空気等信号のない部分が存在する被検体など領域が島状にいくつかに分断されたときに各領域に対して正確に水 ·脂肪分離画像を取得する従来の技術として、「 Λ New Two - Dimensional Phase Unwrapping A丄 goritnm ror MRI Images ； M. Hedley and D. Rosenfeld ； Magnetic Resonance in medicine,

Vol. 24, 177- 181 (1992)」がある。この技術は 2つの分断された島状領域がある場合、第 1の領域と第 2の領域間の最近接点の位相のみを単純に比較することが開示されている。し力し、この技術では局所的に急激な位相変化が起こっている場所や、分断された島状領域間が離れている場合には、所望の画像を得ることが困難なため、水画像、脂肪画像を分離することができず、 1枚の画像の中に、水信号が描出された領域と脂肪信号が描出された領域が混在してしまう。発明の開示

本発明の目的は 1つの画像中に 2つ以上の領域が島状にいくつかに分断されたときでも、領域間でのアーチファタトの発生を抑制が可能な磁気共鳴ィメージング装置を提供することにある。本発明ではエコー時間の異なる複数枚の画像データを取得し、演算によって画像を得る制御手段を有する磁気共鳴ィメ一ジング装置において、前記制御手段は、複数の各領域の位相値に基づいた評価関数を用いて、前記各領域間の位相値の整合を取る磁気共鳴ィメ一ジング装置を開示する。図面の簡単な説明

図 1は MRI装置の構成を示すプロック図である。図 2は従来の F0V内の静磁場不均一分布を説明する図である。図 3は本発明の領域間アンラップアルゴリズムを示す図である。図 4は空間的に離れた複数の部分からなる被検体を撮影した画像を示す図である。図 5は空間的に離れた複数の部分からなる被検体を撮影した画像の位相のプロファイルを示す図である。図 6は本発明の空間的に離れた部分からなる被検体の F0V内の静磁場不均一分布を説明する図である。図 7は本発明の評価値の比較を示す図である。 101は被検体、 102は静磁場磁石、 103は傾斜磁場コイル、 104は RFコイル、 105は RFプローブ、 106は信号検出部、 107は信号処理部、 108は表示部、 109は傾斜磁場電源、 110は RF送信部、 111は制御部、 112はべッドを示す。発明を実施するための最良の形態

MRI装置の構成を図 1により説明する。被検体 101の周囲に静磁場を発生する磁石 102と、該空間に傾斜磁場を発生する傾斜磁場コイル 103と、この領域に高周波磁場を発生する RFコイル 104と被検体 101が発生する MR信号を検出する RFプロープ 105がある。傾斜磁場コイル 103は、 X、 Y、 Ζの 3方向の傾斜磁場コイルで構成され、傾斜磁場電源 109からの信号応じてそれぞれ傾斜磁場を発生する。 RFコイル 104は RF送信部 110の信号に応じて高周波磁場を発生する。 RFプロープ 105の信号は、信号検出部 106で検出され、信号処理部 107で信号処理され、また計算により画像信号に変換される。画像は表示部 108で表示される。傾斜磁場電源 109、 RF送信部 110、信号検出部 106は制御部 111で制御され、制御のタイムチャートは一般にノルスシーケンスと呼ばれている。べッド 112は被検体が横たわるためのものである。

図 2 (b)は被検体 208を撮影した場合の静磁場不均一マップであり、図 2 (a) は、図 2 (b)の M'の線上の位相値をプロットした図である。図 2(a)で、横軸は位置を示し、縦軸は静磁場不均一による位相の回転量を表している。そして、被検体内の静磁場不均一による位相回転量の分布 201を示す。被検体 208の外部の信号成分はほぼ 0なので、位相値は計測ノイズ（ランダム位相）となってしまうので、通常、信号の絶対値でしきい値を設け、信号値が一定値以下の領域の位相値は使用しない。

図 2 (a)では、位相値が + π以上となっている値 202は- πから + πの間に折り返された値 203を取る。また、 - π以下である値 204、 205も- πから + πの間に折り返された値 206、 207を取る。このように、折り返された部分では位相値が - πから + πと不連続な飛び（図中点線）が生じる。し力し、実際の磁場ではこのような不連続な飛びが生じることはありえない。そこで、不連続な飛ぴを除去し、滑らかな静磁場不均一マップを得るために、 USP5, 701， 074等に記載されているアンラップ（巻き戻し）処理を行う。

アンラップ対象が 2つの分断された島状領域となっている場合を図 6に示す。図 6 (b)は領域面積の大きいほうを第 1の領域 610、小さいほうを第 2の領域 611 としたときの静磁場不均一マップを示し、図 6 (a)は図 6 (b)の M'の線上の位相値をプロットした図である。図 6 (a)で、横軸は位置を示し、縦軸は静磁場不均一による位相の回転量を表しており、第 1の領域 610と第 2の領域 611の間に空隙 613が存在している。

このような処理対象に対して、図 3に示す領域間アンラップアルゴリズムを適用する。第 1の領域の領域の位相値の近似関数を求めるフィッティングを行い、第 2の領域の位相値を求め、その後アンラップをする方法を行う。

また、図 3に示す領域間アンラップアルゴリズムを適用する。まず、第 1の領域 610に対して領域内アンラップを行う（301)。次に、第 1の領域 610のァンラップ（301) 後の位相値のみに対して 2次元フィッティングを行う（302)。本実施形態ではフィッティング次数を例えば 4次とする。そして、第 2の領域 611に対してこのフイツティング結果に基づく評価関数によって評価値 f0を求める（303)。評価関数としては、例えば上記フィッティング結果と第 2の領域におけるフイツティング前位相値の差の 2乗和（SD： Standard Deviation標準偏差）を用いる。

続いて、第 2の領域 611に対して領域内アンラップを行う（304)。次に、第 1 の領域 610のアンラップ（301) 後の位相値と第 2の領域のアンラップ（304) 後の位相値に対して、フィッティングを行う（305)。そして、同様に評価関数によって評価値 flを求める（306)。ステップ 306はステップ 303と同じくフィッティング前後の位相値の差の 2乗和を用いて求める。

そして、領域内アンラップを行った第 2の領域の位相値に +2 πのオフセットを与え（307)、同様にフイツティングを行い（308)、評価値 f 2を求める（309)。さらに、領域内アンラップを行った第 2の領域の位相値に - 2 πのオフセットを与え（310)、フィッティングを行い（311)、この状態において評価値 f3を求める（312)。ここで、それぞれ求めた評価値 f0， fl, f2, ： f3を比較し、どの状態が最適であるかを判定する（313)。

また、他の方法として、最初に第 1 の領域に対してアンラップをせず、第 1 の領域の位相値に対してフィッティングを行い、第 2の領域に対して評価関数を求める。その後、第 1の領域及ぴ第 2の領域の位相値のアンラップを行い、それぞれの評価関数値について、比較を行う。そして、どの状態が最適である力判断する。比較の行い方は、図 3に示した方法と同様である。

この判定は、例えば最も小なる 2乗和を選択する。このことを図 4、図 5 'を用いて説明する。図 4は空間的に離れた 2つの部分 401、 402からなる被検体を撮影した画像を示し、図 5はこの画像の位相の B- B'で示したラインでのプロファィルを示す。図 4における範囲 403、 404が、図 5の範囲 501、 502に対応している。曲線 503は範囲 501で位相アンラップした結果を示し、曲線 504は先に示した領域内アンラップ処理を用いて、範囲 502を領域内アンラップした結果を示す。この曲線 504は 2ηπのオフセットを持ち、実際の位相は曲線 505であることを示している。曲線 503'と曲線 504、もしくは、曲線 503と曲線 505を使つて、それぞれフイツティングした結果を曲線 506、曲線 507で示す。

ここで、曲線 504と曲線 506、曲線 505と曲線 507を比較すると、曲線 505 と曲線 507の軌跡はほぼ一致しているのに対し、曲線 504と曲線 506の軌跡は範囲 502 内において、値の差が大きくなつている。このことから、領域内アンラップ後の位相とフイツティング結果の差の 2乗和が最小になる状態が最も適切と判定される。

ステップ 313で判定の結果、評価値 f0、すなわち第 2の領域において領域内アンラップ前の位相値が最適とみなされた場合は、第 2の領域の領域内アンラップを行う前の値に戻す（314)。評価値 fl、すなわち第 2の領域において領域内アンラップ後の位相値が最適とみなされた場合は、第 2の領域の領域内アンラップ後の値を採用する（315)。評価値 f2、すなわち第 2の領域において領域内アンラップ後の位相値に +2 %のオフセットを与えたものが最適とみなされた場合は、さらに +4 π，+6 π，…とオフセットの値を変化させて評価値： fnを求めていき、最適なオフセット値を求め、その値を採用する（316)。同様に、評価値 f3、すなわち位相値に- 2 πのオフセットを与えた場合が最適と判断された場合も、 - 4 π，-6 π，…とオフセットの値を変化させて評価関数の値を求めていき、最適なオフセット値を求め、その値を採用する（316)。

以上の判定の手順を図 7を使って模式的に示す。（a)はステップ 301、 302の結果である。区間 701の太線（実際の位相値）と細線 702 (フィッティング結果）のずれが大きいと、評価値 fOは大きくなる。（b)はステップ 304、 305の結果である。区間 703の太線と細線 704は非常に近い値をとり、評価値 flは小さい値をとる。（c)、（d)はそれぞれ、ステップ 307、 308、及ぴステップ 310、 311 の結果である。区間 705の太線と細線 706、区間 707の太線と細線 708のずれは大きく、評価値 f2、 f3も大きな値をとる。図 7の例では、評価値 flが最小値をとり、ステップ 313で比較、判定の結果、ステップ 315が最適となる。実際にこの処理を大腿部の左右同時撮影（横断面、冠状断面）や腹部横断面同冠状断面などのマルチスライス撮影、大 F0Vでの撮影、頭部横断面の撮影、足首の撮影などに適用し、 300枚以上の撮影を行ったところ、全画像において、領域間アンラップが正しく行われ、アーチファタトのない水'脂肪分離画像が得られた。また、本手法は、グラジェントエコータイプの 3点 Dixon、スピンェコータイプの 3点 Dixonの双方に適用できた。

なお、評価関数と判定式の他の例としては、次のようなものも考えられる。

(a)静磁場不均一による位相回転量分布図を多項式でブイッティングしたときの、最高次数項の係数がすべて定めたしきい値以下である力否かを判定し、しきい値以下なら正しいとする。

(b) 上記フィッティングにおいて、最高次数項の係数の 2乗和を求め、その最小値を採用する。

ここで、（a) は、「磁石の作る静磁場は緩やかに変化し、位相マップは局所的な変化をもたず、フィッティング結果において、高次項の係数は大きな値になることはない」ということを根拠としている。実際の処理では 4次程度のフイツティングを行い、 4次項 (x4_s x3y, x2y2, xy3, y4) の係数がある値以下であればよいとした。実験の結果、この評価関数である程度の処理の安定化ができた。なお、適切なしきい値は F0Vや部位に依存して変化した。

(b) は（a) に改良を加えたもので、あらかじめしきい値を与えるのでなく、各状態での評価関数の値を相対的に比較し、どの状態が最適かを判断する。この場合、 (a) と同じ理由により、（b) の値が最小になる状態が最適であると考えられた。この手法は（a) よりも安定であった。なお、（a) と同様に撮影条件が異なると処理結果に違いが生じた。

上記実施形態において、島状領域のうち、第 1の領域と第 2の領域の決定は、面積の大小以外の方法で決定しても良い。例えば、画像のより中心部に存在する島状領域を第 1の領域とし、周辺部に存在する島状領域を第 2の領域とする。すなわち、処理の手順として面積が多く中心部に位置する領域からアンラップをスタートすることが、早く安定にアンラップを行うために重要である。

評価関数の求め方について、第 2領域についてフイツティング前後の位相差の 2乗和を求めたが、他の方法として第 1と第 2の合計の領域に対してフィッティング前後の位相差の 2乗和を求めても良い。この場合、評価は局所的なフィッティングの状況ではなく被検体全体のフィッティング効果を評価することになる。ただし、対象面積が増え評価値の計算に時間を要するため、より高速な処理が必要となる。

また、上記実施形態では、アンラップすべき島状領域が 2つの場合で説明したが、島状領域は 3つ以上であっても良い。この場合、第 1の領域と第 2の領域の領域間アンラップ処理が完成してから、第 3の領域の領域間アンラップを行う方法で、短時間で安定な結果を得ることができる。なお、第 3の領域に対する評価値は、第 1の領域と第 2の領域の領域間アンラップ、つまり、第 1の領域と第 2の領域との間で整合された位相値に基づく評価関数にて求められる。さらに、第 4の領域がある場合には第 1〜第 3の領域間で整合された位相値に基づく評価関数にて評価値が求められる。また、第 1 と第 2の合計の領域のフィティング効果を評価する手段と同様に、第 1と第 a (aは 3以上の整数）の合計の領域に対してフィッティング前後の位相差の 2乗和を求め、被検体全体のフィッティング効果を評価しても良い。

アンラップ対象となる島状領域の大きさをどの程度までにするかは、本処理を安定に動作させるために十分な検討が必要である。小さくて画像の端にある島状領域を静磁場マップの計算に組み込むことがどの程度のメリットがあるかを考えれば、この判断は自然につく。典型的には、 256 X 256 の画像において、画素数 100程度以下の島状領域については、無視して静磁場マップの確定に使用しないことが望ましい。その理由は、あまり小さな領域は例えば血流アーチファタトゃ体動アーチファタトである可能性があり、その位相値が静磁場分布を反映している可能性が低いからである。また、小さな島状領域で評価値（2 乗和）を求めても、その精度が必ずしも高い保証が無い。さらに、小さな島状領域が、フィッティング関数の全体に与える影響は小さいので、静磁場分布の結果に大きな影響を与えない場合が多い。また、島状領域の数が増えると、演算時間が延長することの考慮も必要である。また画像の端部に臨床上重要な部位が少ないことも判断の基準になる。アンラップ処理の対象とする島状領域の数を制限する方法も有効である。例えば島状領域を画素数の大きい順に 5個のみ演算対象とする。このような制限を用いれば、演算時間をある範囲以内に制限できるメリットがある。

さらに、上記実施形態では、フィッティング関数を 4次としたが、その次数は変更可能である。この場合、次数の決定には下記の配慮が必要である。次数が低いと得られるフィッティング結果の静磁場不均一マップはなだらかになる。従って、処理結果は安定になるが局所的な磁場の変動に追随できない欠点がある。例えば、オープン型磁石を使った MRI装置では一般に磁場中心から離れると急激に磁場均一性が低下することが知られている。従って、撮影視野（F0V) を大きくしたり、磁場中心以外の位置で撮影を行ったりするときに、静磁場の変化を正しく求められなくなる。その結果として、静磁場補正を施した水脂肪分離画像に水と脂肪が混在することになり、臨床診断の制度は著しく低下する。また、生体起因の静磁場乱れとして、例えば鼻腔、副鼻腔など空隙が複雑に入り込んでいる部位は局所的に静磁場が乱れる事が知られている。このような部位についても、次数を下げると正しくフィッティングが出来なくなる。他方、フィッティング次数を増やすとこれらの問題はなくなるが、被検体の周辺部でフィッティング結果が発散し不安定になる場合がある。また、フィッティングの演算時間が延びる傾向もある。発明者らはこのような経験から、フイツティング次数を 4次として良好な結果を得た。具体的には、オープン MRI装置において、 F0V320雇程度まで安定な処理が行えた。

島状のアンラップ対象領域が生まれる理由は、前述したように、信号の絶対値が小さい画素をノイズとみなして、演算対象からはずすためである。この処理は、同業者にとっては当然の処理であるが、その閾値をどこにするかは、処理を安定に行うために十分検討する必要がある。典型的には、画像の最大値の 5%から 40%程度が望ましい。あまり閾値が大きいと各島状領域の面積が小さくなり、少ないデータで静磁場不均一マップを作るため、求める静磁場不均一マップの精度が下がる。また島状領域の数が増える傾向がある。閾値が小さいと島状領域が分断される数は減るがノィズ成分が入り込むので、やはり静磁場不均一マップの精度が下がる。発明者らの検討によれば、閾値は 10%から 30% の範囲が安定に処理を行うために適当であった。

分断される島状領域の数は、測定するエコーの特質によっても変わる。水と脂肪が同相になるエコー（インフェイズエコー）同士の演算で位相マップを作成すると、島状領域の数は少なくなる。しかし、水と脂肪が逆相になるエコー (アウトフェイズエコー）を位相マップを作成する演算に使用すると、島状領域の分断数は多くなる。その理由は、アウトフェイズ画像では、水と脂肪が混在する画素においてそれぞれの信号成分が打ち消し合い、信号値が小さくなるためである。従って、静磁場不均一マップを安定に作成するには、インフェイズエコーと本発明の島状領域間の位相アンラップを併用することが最も望ましい。なお、本発明の効果は、アウトフェイズエコーとの併用であっても、従来の処理に比べ大幅に処理の安定性と精度が向上することは、前述の説明で明らかである。また、インフェイズエコーとアウトフェイズエコーの位相差で静磁場不均一マップを求める場合にも本発明は適用できる。実際にこのようなデータに本処理を適用した結果、脳底部の画像や、腹部 TRS断面、腹部 C0R断面のアンラップ処理が著しく向上した。

分断された島状領域をどの順番でアンラップしていくのかも重要である。経験的には、画像中心部に近い領域から行うのが望ましい。その理由は、画像中心は磁場中心に一致することが多いので、静磁場均一度が高く空隙内でのアンラップ数も少なく、得られる結果が正しい確率が高いからである。または、大きな島状領域から行うのがよい。その理由は、大きな島状領域では前述のように評価値の精度が高いので正確な判定が出来るため、得られる結果が正しい確率が高いからである。さらには、アンラップ済みの領域との空隙が狭い島状領域から行うのが良い。これも空隙内でのアンラップ数も少なく、得られる結果が正しい確率が高いからである。すなわち正しい結果が得られやすい領域から順にアンラップを適用していく。また、これらの処理順番は、すべて簡単な数式のみで記述できるため、処理の自動化を妨げることは無い。

上記実施形態では、評価値を求める準備段階として、フィッティングを用いている。フィッティング結果と元データを比較する方法は適切な方法の一つではあるが、上記実施形態に限定されるものではない。すなわち、分断された島状領域間のアンラップを行う際に、複数のアンラップ状態を想定し、それぞれの結果に対して評価関数を使って判定を行い、最も適切と思われる評価値、つまりアンラップ条件を採用することができればよい。

上記実施例では、 2次元の画像に対して説明を行ったが、 3次元画像への拡張は当然可能であることは、言うまでも無い。

本発明は、以上のように構成されたので、アンラップすべき領域が島状にいくつかに分断されたときでも、適切に領域間アンラップ行うことができる。よつて、位相回転量を正確に求めることができるため、 1枚の画像の中に水信号が描出された領域と脂肪信号が描出された領域が混在するようなアーチファタトの発生を抑制することが可能となる。

Claims

求の

1. エコー時間の異なる複数枚の画像データを取得し、演算によって画像を得る制御手段を有する磁気共鳴ィメ一ジング装置において、前記制御手段は、複数の各領域の位相値に基づいた評価関数を用いて、前記各領域間の位相値の整合を取ることを特徴とする磁気共鳴ィメ一ジング装置。

2. 前記制御手段は、位相回転量の主値周りを除去するアンラップ処理を複数領域に対し実行すると共に前記各領域間の位相値の整合を取ることを特徴とする請求項 1記載の磁気共鳴ィメ一ジング装置。

3. 前記制御手段は、複数の各領域間の位相値に基づいた評価値を複数求め、前記複数の評価値から所望の評価値を選択し、前記所望の評価値に基づき前記各領域間の位相値の整合を取ることを特徴とする請求項 2記載の磁気共鳴イメージング装置。

4. 前記制御手段は、前記評価値を前記位相回転量の分布図をアンラップ後にフィッティングした位相値とフイツティング前の位相値とに基づいて求めることを特徴とする請求項 3記載の磁気共鳴ィメ一ジング装置。

5. 前記制御手段は、前記複数の評価値のうち最小となる評価値を選択することを特徴とする請求項 3又は請求項 4記載の磁気共鳴ィメ一ジング装置。

6. 被検体の置かれた空間に静磁場を発生するための静磁場発生手段と、前記被検体に対し核磁気共鳴現象を起こすための高周波パルスを繰り返し印加する高周波パルス印加手段と、前記被検体にスライス、位相エンコード、読み出し方向の各傾斜磁場を印加する傾斜磁場発生手段と、前記被検体からの複数のエコー信号を検出する受信手段と、前記高周波パルスの 1回の繰り返し時間内にエコー時間を異ならせた複数のエコー信号を発生させると共に前記各エコー信号の発生と共に読み出し傾斜磁場を印加し、繰り返し毎に位相ェンコ一ドを変化させるよう各手段を制御する制御手段を有する磁気共鳴イメージング装置において、前記制御手段は、分断された領域を持つ被検体に対し各領域毎に位相値の整合を行う第 1のアンラップ処理を実行し、各領域間の位相値の整合を行う第 2のアンラップ処理を実行すると共に、前記第 2のアンラップ処理は複数の位相値整合のパターンを求め、この複数のパターンから所望のパターンを選択して位相値整合を行うことを特徴とする磁気共鳴イメージング装置。

7. ェコ一時間の異なる複数枚の画像に対し水と脂肪に分離した画像を得るための画像処理をするための画像処理方法において、

(a) 分断された領域が複数存在する画像に対し各領域毎に静磁場不均一による位相回転量の分布図を求めるためのアンラップ処理を行うステップ、

(b) 各領域間の位相値の整合を取るため複数のパラメータを求め、複数のパラメータから所望のパラメータを選択し、選択されたパラメータに基づいて各領域間のアンラップ処理を行うステップ、を含むことを特徴とする画像処理方法。

8. ェコ一時間の異なる複数枚の画像に対し水と脂肪に分離した画像を得るための画像処理をするための画像処理方法において、

(a) 分断された領域が複数存在する画像に対し基準となる領域に対し静磁場不均一による位相回転量の分布図を求めるためのアンラップ処理を行うステップ、

(b) 前記基準領域の静磁場不均一による位相回転量の分布図に対しフィッティング処理を行うステップ、

(c) 前記基準領域と隣り合う第 1の領域の位相値と、ステップ（b) により得られた基準領域の位相値とにより評価値を求めるステップ、

(d) 第 1の領域に対しアンラップ処理を行うステップ、

(e) 基準領域と第 1の領域に対しフイツティング処理を行うステップ、

(f) 基準領域と第 1の領域の位相値と、ステップ（e) により得られた基準領域と第 1の領域の位相値とにより評価値を求めるステップ、

(g) 第 ηの領域に対し所望のオフセットを与えるステップ、

(h) 基準領域と第 1の領域に対しフイツティング処理を行うステップ、

(i) 基準領域と第 1の領域の位相値と、ステップ（h) により得られた基準領域と第 1の領域の位相値とにより評価値を求めるステップ、

(j) オフセット値を変えて（g) 〜（i) を所望の回数繰り返すステップ、 (k)ステップ（c) , (f)， (i)及びステップ（j)の繰り返しにて得られた複数の評価関数を比較するステップ、

(1) ステップ（k) の比較結果に基づいて所望の評価値を選択するステツプ、

On) ステップ（1) にて選択された評価値に基づいて基準領域と第 1の領域との位相値の整合を行うステップ、を含むことを特徴とする画像処理方法。

9. 前記画像処理方法において、分断された領域が 3個以上存在する場合において、さらに

(n) ステップ（m) で位相値整合された基準領域と第 1 の領域とを新たな基準領域とし、ステップ（a)〜（n) を繰り返すステップ、を含むことを特徴とする請求項 8記載の画像処理方法。

10. 前記制御手段は、少なくとも 2つ以上に分断された領域を含む画像を撮り、第 1の領域の位相値に対して近似関数を求め、前記第 1の領域と対向する第 2の領域の位相値は前記近似関数に適応させ、水 ·脂肪分離画像を得ることを特徴とする請求項 1曾己載の磁気共鳴ィメ一ジング装置。

11. 前記制御手段は、前記第 1の領域に対して静磁場不均一による位相回転を補正し、前記近似関数を求め、前記第 2の領域の位相値は前記近似関数に適応させることを特徴とする請求項 10記載の磁気共鳴ィメ一ジング装置。

12. 前記第 1 の領域は前記画像の中心付近に位置することを特徴とする請求項 10記載の磁気共鳴ィメ一ジング装置。

13. 前記第 1 の領域は画像の最も広い領域であることを特徴とする請求項 10 記載の磁気共鳴ィメ一ジング装置。

14. 前記第 2の領域は前記第 1の領域と空隙が狭い領域から選択されることを特徴とする請求項 10記載の磁気共鳴ィメ一ジング装置。

15. 第 n ( nは 3以上の整数）の領域の位相値を求める際、前記第 1の領域から第 n- 1の領域の位相値の近似関数を求め、前記第 nの領域の位相値は前記第 1の領域から前記第 n-lの領域の位相値の前記近似関数に適応させることを特徴とする請求項 10記載の磁気共鳴ィメ一ジング装置。

16. 被検体から画像データを取得し、演算によって画像を得る制御手段を有する磁気共鳴ィメ一ジング装置において、前記制御手段は少なくとも 2つ以上に分断された領域を含む画像を取り、第 1の領域及び前記第 1の領域と対向する第 2の領域の合計領域に対してそれぞれの位相値の評価関数を求め、前記評価関数から前記第 1の領域の位相値に前記第 2の領域の位相値を適応させ、第 1の領域及び第 2の領域の水 ·脂肪分離画像を得ること特徴とする磁気共鳴イメージング装置。

17. 前記制御手段は前記第 1の領域及び前記第 2の領域に対して静磁場不均一による位相回転を補正し、それぞれの位相値の評価関数を求めること特徴とする請求項 16記載の磁気共鳴ィメ一ジング装置。

18. 前記第 1 の領域は前記画像の中心付近に位置することを特徴とする請求項 16記載の磁気共鳴ィメ一ジング装置。

19. 前記第 1 の領域は画像の中で最も広い領域であることを特徴とする請求項 16記載の磁気共鳴ィメ一ジング装置。

20. 前記第 2の領域は前記第 1の領域と空隙が狭い領域から選択されることを特徴とする請求項 16記載の磁気共鳴ィメ一ジング装置。

21. 第 n ( nは 3以上の整数）の領域の位相値を求める際、前記第 nの領域の位相値は第 1の領域から第 nの領域の合計領域の前記近似関数に適応させることを特徴とする請求項 16記載の磁気共鳴ィメ一ジング装置。

22. 複数枚の画像に対し水と脂肪に分離した画像を得るための画像処理をするための画像処理方法において、

(a)分断された領域が複数存在する画像に対し基準となる第 1の領域に対し静磁場不均一による位相回転量の分布図を求めるためのアンラップ処理を行うステップと、

(b)前記第 1の領域の静磁場不均一による位相回転量の分布図に対しフィッティング処理を行うステップと、

(c) 前記第 1の領域と隣り合う第 2の領域の位相値及びステップ（b ) により得られる前記第 1の領域の位相値により評価値を求めるステップとを含むことを特徴とする画像処理方法。

23. (d) 前記第 2の領域に対しアンラップ処理を行うステップと、 (e)前記第 1の領域と前記第 2の領域に対しフイツティング処理を行うステツプと、

(f) 前記第 1の領域と前記第 2の領域の位相値及ぴステップ（e) により得られる前記第 1の領域と前記第 2の領域の位相値とにより評価値を求めるステップとを含むことを特徴とする請求項 22記載の画像処理方法。

24. (g) 第 2の領域の位相値に ± 2m π (mは自然数）のオフセットを与え評価値を求めるステップを含むことを特徴とする請求項 23記載の画像処理方法。

25. (j)ステップ（c) , (f)，（g)にて得られた複数の評価値に対する評価関数を比較し、所望の評価値を選択するステップを含むことを特徴とする請求項 24記載の画像処理方法。

26. ( ) 前記第 1の領域と前記第 2の領域の位相値と第 X ( Xは 3以上の整数）の領域の位相値に対し、フッティング処理を行うステップを含むことを特徴とする請求項 25記載の画像処理方法。