"被検体が置かれる空間に静磁場を発生する静磁場発生手段と、前記被検体に高周波磁場パルスを照射する照射手段と、高周波磁場パルスの照射により前記被検体から生じる核磁気共鳴信号を受信する受信手段と、前記核磁気共鳴信号に位置情報を付加するための傾斜磁場を印加する傾斜磁場印加手段と、前記照射手段、前記傾斜磁場発生手段及び前記受信手段を所定のパルスシーケンスに従って動作させて異なる2つのエコー時間の撮影を実行する制御手段と、前記異なるエコー時間の核磁気共鳴信号から原画像を再構成する画像再構成手段と、を備える磁気共鳴撮影装置であって、前記原画像から、静磁場不均一による位相差マップを算出する位相差マップ算出手段と、前記位相差マップから前記原画像のうちの1画像を位相補正して位相補正画像を算出する位相補正画像算出手段と、前記位相補正画像と位相補正を実施しない前記原画像の2画像から水と脂肪を分離した画像を算出する水脂肪分離画像算出手段とを有することを特徴とする磁気共鳴撮影装置。請求項1記載の磁気共鳴撮影装置であって、前記制御手段は、水と脂肪との位相差が正負逆とならず、かつ、水と脂肪との位相差がπの整数倍とならないよう、前記異なる2つのエコー時間を設定するエコー時間設定手段を有することを特徴とする磁気共鳴撮影装置。請求項2記載の磁気共鳴撮影装置であって、前記位相差マップ算出手段は、前記原画像から、水と脂肪の比を算出する水対脂肪比算出手段と、前記水対脂肪比から、位相マップを算出する位相マップ算出手段と、前記位相マップの基準点から隣り合う点の位相差が最も小さくなる点を、領域を拡大しながら探索し、最小位相差マップを算出する最小位相差マップ算出手段と、前記最小位相差マップから、静磁場不均一による位相差マップを判定する位相差マップ判定手段と、を有することを特徴とする磁気共鳴撮影装置。前記水対脂肪比算出手段は、水と脂肪との割合の和が一定となる条件を付加して代数的に算出する算出手段を有することを特徴とする請求項3記載の磁気共鳴撮影装置。前記水対脂肪比算出手段は、水と脂肪の割合の和が一定となる条件を付加して数値解析的に算出する算出手段を有することを特徴とする請求項3記載の磁気共鳴撮影装置。前記位相差マップ判定手段は、前記最小位相差マップの微分値から、標準偏差が最小となる最小位相差マップを選択する選択手段を有することを特徴とする請求項3記載の磁気共鳴撮影装置。上述した非特許文献2は、3回の計測を実施する必要があるため、計測時間が延長してしまうという課題がある。また、非特許文献3の方法は、傾斜磁場パルスを反転することで、水と脂肪が同位相の画像と逆位相の画像を1回の計測で取得している。しかし、tinとtoutの時間間隔dtは、磁場強度に反比例する。nは整数としたとき、たとえば、磁場強度が3.0テスラのとき、dtは1.122.24×nms、磁場強度が1.5テスラのとき、dtは2.244.48×nms、磁場強度が0.5テスラのとき、dtは6.7013.39×nmsとなる。磁場強度が3.0テスラや1.5テスラのとき、dtが狭いため、計測帯域を広げたり、計測点数を減少させたりするなど計測上の制限が生じる。計測帯域を広げると、信号対雑音比以下、SNRというが低下し画像が劣化する。また、計測点数を減らすと、画像の空間分解能が低下し、微小な疾患の診断が困難になる。整数nの数を調整することで時間間隔dtを広げられるため、計測制限はある程度緩和される。しかし、取得する2枚の画像の両方で任意のコントラストの画像を取得することは困難である。磁場強度が0.5テスラのとき、dtが広いため、1スライスあたりの計測時間が延長してしまう。すなわち、2次元マルチスライス計測においては、所望の繰り返し時間TRあたりのスライス数を減らさなければならないため、取得したいスライス数が減少してしまう。また、3次元計測においては、単に計測時間が延長してしまう。計測時間が延長すると、被検体の体動アーチファクトが増大する。アーチファクトが増大すると、演算処理に対する誤差が大きくなり、水と脂肪の分離が困難となる。また、水と脂肪が混在する人体組織では、水信号と脂肪信号が打ち消しあうため、逆位相の画像のSNRは低下する。そのため、静磁場不均一による位相まわりがノイズによって正確に計算できなくなる。その結果、演算処理に対する誤差が大きくなり、水と脂肪の分離が困難となる。更に、得られた2つの分離画像が水画像、脂肪画像のどちらかであるかの判別ができないという課題がある。本発明では、水脂肪の分離撮影において、水と脂肪との間の位相が同位相および逆位相となるようなエコー時間で計測しなくても良好な水脂肪分離を実現することを課題とする。さらに、得られた2つの分離画像ついて、水画像か脂肪画像かの判別が可能となる装置および方法を実現することを課題とする。本発明は、水と脂肪の位相差が、正負逆とならず、かつ、πの整数倍とならない2つのエコー時間で2つの原画像を計測する。例えば、あるエコー時間t1における水と脂肪の位相差をP1、エコー時間t1からdt後のエコー時間t1dt=t2における水と脂肪の位相差をP2としたとき、P1≠n1πn1は整数5かつP2≠m1πm1は整数6かつ2n2πP1≠-2m2πP2n2およびm2は整数7を満たすようなエコー時間t1およびt2で計測し2つの原画像を得る。前記条件で取得した2つの原画像から、代数的あるいは数値解析的に2つの水対脂肪比マップを算出する。算出された2つの水対脂肪比マップから、2つの位相マップを算出する。算出された2つの位相マップを組み合わせて、空間的な位相差分がそれぞれ最小になるような2つの最小位相差マップを算出する。最小位相差マップを空間的に微分した微分マップのばらつきから、正しい最小位相差マップかどうかを判別する。判別された位相差マップを用いて、原画像を位相補正する。補正された原画像から水脂肪分離演算処理を行う。具体的には、被検体が置かれる空間に静磁場を発生する静磁場発生手段と、前記被検体に高周波磁場パルスを照射する照射手段と、高周波磁場パルスの照射により前記被検体から生じる核磁気共鳴信号を受信する受信手段と、前記核磁気共鳴信号に位置情報を付加するための傾斜磁場を印加する傾斜磁場印加手段と、前記照射手段、前記傾斜磁場発生手段及び前記受信手段を所定のパルスシーケンスに従って動作させて異なる2つのエコー時間の撮影を実行する制御手段と、前記異なるエコー時間の核磁気共鳴信号から原画像を再構成する画像再構成手段と、を備える磁気共鳴撮影装置であって、前記原画像から、静磁場不均一による位相差マップを算出する位相差マップ算出手段と、前記位相差マップから前記原画像のうちの1画像を位相補正して位相補正画像を算出する位相補正画像算出手段と、前記位相補正画像と位相補正を実施しない前記原画像の2画像から水と脂肪を分離した画像を算出する水脂肪分離画像算出手段とを有することを特徴とする磁気共鳴撮影装置を提供する。上記原画像計測手段は、例えば、前記式5、式6、式7を満たす2つのエコー時間で原画像を計測する。また、上記位相差マップ算出手段は、例えば、前記2つの原画像から、代数的あるいは数値解析的に2つの水対脂肪比マップを算出する水対脂肪比マップ算出手段と、前記2つの水対脂肪比マップから、2つの位相マップを算出する位相マップ算出手段と、前記2つの位相マップの各点の組み合わせから、空間的な位相差分がそれぞれ最小となるように2つの最小位相差マップを算出する最小位相差マップ算出手段と、前記最小位相差マップを空間的に微分した微分マップのばらつきからから、正しい最小位相差マップかどうかを判定する位相差マップ判定手段と、を有する。本発明によれば、水脂肪分離のために異なるエコー時間に画像を取得する手法において計測上の制限が緩和され、かつ良好な水脂肪分離が実現できる。さらに、得られた2つの分離画像ついて、水画像か脂肪画像かの判別が可能となる。以下、本発明のMRI装置の実施の形態を、図面を参照して説明する。図1a~cに、本発明が適用される磁気共鳴撮影装置の全体構成と外観図を例示する。図1aはソレノイドコイルで静磁場を発生するトンネル型磁石を用いた水平磁場方式のMRI装置である。図1bは開放感を高めるために磁石を上下に分離したハンバーガー型オープン型の垂直磁場方式のMRI装置である。また、図1cは、図1aと同じトンネル型のMRI装置である。しかし、磁石の奥行を短くし且つ斜めに傾けることによって、開放感を高めている。本発明は、これらの形態のMRI装置を含む公知の構造のMRI装置に適用することができるが、これに限定されるものではなく、MRI装置の形態やタイプを問わず適用可能である。本実施形態のMRI装置の構成の一例を図2に示す。このMRI装置は、被検体1が置かれる空間に、静磁場を発生する静磁場コイル2と、静磁場に直交する三方向の傾斜磁場を与えるための傾斜磁場コイル3と、静磁場均一度を調整できるシムコイル4と、被検体1の計測領域に対し高周波磁場を照射する計測用高周波コイル5以下、単に送信コイルというと、被検体1から発生する核磁気共鳴信号を受信する受信用高周波コイル6以下、単に受信コイルというとを備えている。静磁場コイル2は、図1に示した装置の構造に応じて、種々の形態のものが採用される。傾斜磁場コイル3およびシムコイル4は、それぞれ傾斜磁場用電源部12およびシム用電源部13により駆動される。図2では、送信コイル5と受信コイル6とを別個に示しているが、送信用と受信用を兼用する一つの高周波コイルのみを用いる構成もある。送信コイル5が照射する高周波磁場は、送信機7により生成される。受信コイル6が検出した核磁気共鳴信号は、受信機8を通して計算機9に送られる。計算機9は、核磁気共鳴信号に対して様々な演算処理を行い、画像情報などを生成する。上記演算処理には、水と脂肪の分離画像を算出するための各種演算処理が含まれる。計算機9には、ディスプレイ10、記憶装置11、入力装置15などが接続されている。上述した画像情報や、撮影領域の情報は、ディスプレイ10に表示したり記憶装置11に記録されたりする。入力装置15は、測定条件や、撮影領域情報の設定、演算処理に必要な条件などを入力するためのもので、これら測定条件等も必要に応じて記憶装置11に記録される。シーケンス制御装置14は、傾斜磁場コイル3を駆動する傾斜磁場用電源部12、シムコイル4を駆動するシム用電源部13、送信機7および受信機8を制御する。制御のタイムチャートパルスシーケンスは撮影方法によって予め設定されており、記憶装置11に格納されている。次に本実施形態のMRI装置で実行されるパルスシーケンスについて説明する。図3に、本実施形態において採用されるパルスシーケンスのタイムチャートを示す。このパルスシーケンスでは、まずRFパルス301を照射し、被検体のスピンを励起する。この際、被検体の特定のスライスを選択するためにスライス選択傾斜磁場GsをRFパルス301と同時に印加する。続いて核磁気共鳴信号に位相エンコードするための位相エンコード傾斜磁場Gpを印加し、その後、読み出し傾斜磁場Gr302を印加して、最初のRFパルス301照射から時間TE1後にエコー信号第1エコー303を計測する。更に極性の反転した読み出し傾斜磁場Gr304を印加して、第1エコー303計測から時間dt後にエコー信号第2エコー305を計測する。このときのエコー時間TE2は、TE2=TE1dtとなる。第1エコー303を計測する際のエコー時間TE1および、第2エコー305を計測する際のエコー時間TE2は、以下の式5、式6、および式7の条件を満たすよう選択される。すなわち、エコー時間TE1における水と脂肪の位相差をP1、エコー時間TE1からdt後のエコー時間TE1dt=TE2における水と脂肪の位相差をP2としたとき、P1≠n1πn1は整数5かつP2≠m1πm1は整数6かつ2n2πP1≠-2m2πP2n2およびm2は整数7を満たすエコー時間TE1およびエコー時間TE2が選択される。このように、本実施形態では、エコー時間TE1とエコー時間TE2とを、式7に示すように、水と脂肪の位相差が正負逆とならないように設定する。また、式5および式6に示すように、水と脂肪との位相差P1およびP2が、それぞれ、πの整数倍とならないよう各エコー時間TE1とエコー時間TE2を設定する。このシーケンスを位相エンコード傾斜磁場Gpの強度を変化させながら所定回数繰り返す。例えば128回、256回等、画像再構成に必要な数のエコー信号を繰り返し取得する。そして、繰り返し回数分の第1エコー303から1つの画像第1エコー画像が再構成され、繰り返し回数分の第2エコー305からもう1つの画像第2エコー画像が再構成される。これらは後述する水画像および脂肪画像を算出するための演算用の画像として使用される。なお、図3では、グラディエントエコー型のパルスシーケンスを示したが、図4に示すようなスピンエコー型パルスシーケンスでもよい。図4に、スピンエコー型パルスシーケンスを示す。このパルスシーケンスでは、まずRFパルス401を照射し、被検体のスピンを励起する。この際、被検体の特定のスライスを選択するためにスライス選択傾斜磁場GsをRFパルス401と同時に印加する。続いて核磁気共鳴信号に位相エンコードするための位相エンコード傾斜磁場Gpを印加し、更にスピンを反転させるためのRFパルス402をスライス選択傾斜磁場Gsとともに照射する。その後、読み出し傾斜磁場Gr403を印加して、最初のRFパルス401照射から時間TE1後にエコー信号第1エコー404を計測する。更に極性の反転した読み出し傾斜磁場Gr405を印加して、第1エコー404計測から時間dt後にエコー信号第2エコー406を計測する。第1エコー404を計測する際のエコー時間TE1および、第2エコー406を計測する際のエコー時間TE2は、グラディエント型パルスシーケンスと同様、式5および式6、式7の条件を満たすように選択される。グラディエントエコー型の場合にもスピンエコー型の場合にも、図3および図4に示すように1回の繰り返し時間内でエコー時間の異なる2つの信号を計測してもよいし、2回の計測でエコー時間の異なる信号を計測するようにしてもよい。本実施形態では、図5aに示すように、水と脂肪とを分離した画像を算出する。すなわち、シーケンス制御装置14により各部を制御し、上記図3または図4に示すパルスシーケンスを実行し、異なる時間TE1およびTE2でそれぞれエコー信号を得る撮影;ステップ511。そして、計算機9では、異なる時間TE1およびTE2で得たエコー信号から、それぞれ原画像を再構成する画像再構成;ステップ512。次に、計算機9では、ステップ512で得た2つの原画像から、静磁場不均一による位相差マップを算出し位相差マップ算出;ステップ513、得られた位相差マップを用いて、ステップ512で得た2つの原画像のうちの一方を補正し位相補正画像算出;ステップ514、補正後の画像と、ステップ514で補正しなかった原画像とから、水と脂肪とを分離した画像を算出する水脂肪分離画像算出;ステップ515。次に、上記ステップ513の、本実施形態の磁気共鳴撮影装置における静磁場不均一による位相差マップを算出する手順について、図5bを用いて説明する。ステップ501において、2つの水対脂肪比マップを算出する。まず、水対脂肪比マップを算出するための算出式を求める。第1エコーを計測する際のエコー時間t1における水の位相をσw1、脂肪の位相をσF1、静磁場不均一による位相まわりをδ1、第2エコーを計測する際のエコー時間t2における水の位相をσw2、脂肪の位相をσF2、静磁場不均一による位相まわりをδ2、初期位相をφ、第1エコーとして得られた信号の絶対値をr1、偏角をθ1、第2エコーとして得られた信号の絶対値をr2、偏角をθ2、ボクセルに含まれる水信号の割合をRW、脂肪信号の割合をRF、比例定数をkとすると、得られた第1エコーと第2エコーはそれぞれ式8、式9で表現できる。次に、式9を式8で除算し、式10を得る。ただし、δ=δ2-δ1、r=r2/r1、θ=θ2-θ1である。次に、式10について両辺の絶対値の2乗をとる。ここで、ボクセルに含まれる水信号の割合RW、と脂肪信号の割合RFとに、両者の和が一定となる条件を付加する。例えば、RWとRFとには、式12の関係があると仮定する。RWRF=112式11および式12から、以下の式13および式14を得る。を得る。ここで、AおよびBは、以下の式15および式16で表される。σw1、σF1、σw2、σF2、は全て既知の値であり、rは計測された2つの原画像から算出できる。このため、式13~式16より2つの水対脂肪比マップRW:RF、RW-:RF-を、以下の式17および式18のように算出できる。ここで、RWの符号が正のときの水信号の割合をRWと定義したときに式12から求まるRFをRFと定義する。また、RWの符号が負のときの水信号の割合をRW-と定義したときに式12から求まるRFをRF-と定義する。ここで、式11および式12から数値解析的にRW、RF、RW-、RF-を求めてもよい。すなわち、式12の関係を保ちつつ、式11の左辺のRWを0から1まで変化させたときRFを1から0まで変化させたとき、複素除算後の絶対値rとの差が最も小さくなる組み合わせを探索することによってRW、RF、RW-、RF-を求めてもよい。次にステップ502において、ステップ501で算出した2つの水対脂肪比マップから2つの位相マップδ、δ-を算出する。ここで、RWとRFから算出した位相マップをδと定義する。また、RW-とRF-から算出した位相マップをδ-と定義する。以下のように、σF1が0~πの範囲にあるとき、δは水信号が主成分の領域で正しい静磁場マップを表し、δ-は、脂肪信号が主成分の領域で正しい静磁場マップを表している。一方、σF1が-π~0の範囲にあるとき、δは脂肪信号が主成分の領域で正しい静磁場マップを表しており、δ-は、水信号が主成分の領域で正しい静磁場マップを表す。式10にσw1、σF1、σw2、σF2、RW、RFをそれぞれ代入し、式19~式21を得る。また、式10にσw1、σF1、σw2、σF2、RW-、RF-をそれぞれ代入し、式22~式24を得る。ここで、s=s-=rであるため、求める位相マップδ、δ-は、以下の式25および式26のようにそれぞれ算出できる。δ=θ-ω25δ-=θ-ω-26続いて、ステップ503において、2つの位相マップの各点において、空間的な位相変化が最小となる2つの最小位相差マップを算出する。その手順を、モデルを用いて説明する。まず、水の化学シフトfWを0Hz、脂肪の化学シフトfFを224Hzとし、簡便のため静磁場不均一は無いものとする。そして、水と脂肪の割合が0:100、25:75、50:50、75:25、100:0となる領域が連続的に並んだ1次元モデルを考える。ここで、第1エコーのエコー時間t1と、第2エコーのエコー時間t2の条件を以下のA~Dの4つのパターンで比較する。At1を4.5ms同位相、t2を6.0msBt1を6.7ms逆位相、t2を8.2msCt1を4.5ms同位相、t2を6.7ms逆位相Dt1を6.7ms逆位相、t2を8.9ms同位相上記各条件において、同位相とは、そのエコー時間において、水脂肪位相差が0であることを意味し、逆位相とは、そのエコー時間において、水脂肪位相差がπであることを示す。条件AおよびBは本実施形態の条件であり、条件CおよびDは従来用いられている条件である。なお、条件Aのエコー時間がt1の水脂肪位相差は0であるが、式5~式7の条件を満たせば位相差は0でなくてもよい。また、条件Bのエコー時間がt1の水脂肪位相差はπであるが、式5~式7の条件を満たせば位相差はπでなくてもよい。ここで、条件Aにおける位相マップδ、δ-を図6aに、条件Bにおける位相マップδ、δ-を図6bに、それぞれ示す。また、条件Cにおける位相マップδ、δ-を図7aに、条件Dにおける位相マップδ、δ-を図7bに、それぞれ示す。図6a601~605および図7a701~705は、水と脂肪の割合が0:100、25:75、50:50、75:25、100:0の領域をそれぞれ示している。また、図6606:黒丸および図7706:黒丸は、各条件におけるδのグラフ、図6607:白枠および図7707:白枠は、各条件におけるδ-のグラフをそれぞれ示している。図6aに示すように、本実施形態における条件では、水と脂肪の割合によってδとδ-の位相差608が変化していることがわかる。水と脂肪の割合によりδとδ-の位相差608が急激に変化することを位相とびと呼ぶ。静磁場不均一による位相まわりは、金属等の磁性体が近傍にない限り、空間的になだらかに変化する。このため、プロトンが連続的に存在している領域では、図6aに示す位相差608のような位相とびは生じない。したがって、脂肪が主要な成分を占めている601および602の領域ではδ-が静磁場不均一による位相まわりを表し、水が主要な成分を占めている604および605の領域では、δが静磁場不均一による位相まわりを表していることになる。なお603の領域では、水と脂肪の割合が同じであるため、δ、δ-ともに静磁場不均一による位相まわりを表している。なお、図6bに示すように、位相差609は、第1エコーのエコー時間t1が逆位相のときは、図6aに示す位相差608と逆の関係になる。すなわち、脂肪が主要な成分を占めている601および602の領域ではδが静磁場不均一による位相まわりを表し、水が主要な成分を占めている604および605の領域では、δ-が静磁場不均一による位相まわりを表している。一方、図7aに示すように、従来の計測条件では、水と脂肪の割合によらずδとδ-の位相差708は、πあるいは-πとなる。そのため、δとδ-のどの領域が正しい位相まわりを反映しているかが判断できない。なお、図7bに示すように、位相差709は、第1エコーのエコー時間t1が逆位相のときは、図6bに示す位相差609と同様、図7aに示す位相差708と関係が逆になる。本実施形態では、水画像と脂肪画像を区別するために、図6aに示す位相差608や図6bに示す位相差609にあるような水と脂肪の割合が変化することによって生じる位相とびを利用する。上述したように、静磁場不均一による位相まわりは、金属等の磁性体が近傍にない限り、空間的になだらかに変化する。したがって、算出した位相マップδ、δ-の任意の点を基準点とし、基準点と隣り合う点との位相差がもっとも小さい点を、領域を拡大しながら探索し、算出していき、最終的に最小位相差マップδ1およびδ-1を算出する。この最小位相差マップδ1およびδ-1のうち,どちらかが静磁場不均一による位相まわりを正しく表現していることになる。ここで、2つの最小位相差マップδ1およびδ-1を算出する手順を、図6aを使って説明する。[手順1]δおよびδ-での基準点x0を定める。ここでは、位相差|δ-δ-|が最も大きい点を基準点x0として設定する。たとえば、図6aの点610を基準点x0として設定する。[手順2]最小位相差マップδ1およびδ-1を算出するために、基準点x0におけるδ1およびδ-1の値を定める。ここでは、δ1x0=δx0、δ-1x0=δ-x0とする。[手順3]δ1およびδ-1を算出するために、δ1x0に対して、隣り合う点δx01およびδ-x01との位相差分のうち、最も小さくなる点をδ1x01に設定する。例えば、図6aにおいて、601、602の領域では、δ1x01にδx01が設定され、604、605の領域では、δ1x01にδ-x01が設定される。また、δ-1x0に対して、隣り合う点δx01およびδ-x01との位相差分のうち、最も小さくなる点をδ-1x01に設定する。[手順4]手順3をデータ全体で処理が終了するまでを繰り返し、δ1、δ-1を算出する。図8に最小位相差マップδ1およびδ-1をそれぞれ示す。801:白丸はδ1、802:米印はδ-1をそれぞれ示している。最小位相差マップδ1およびδ-1のそれぞれが隣り合う点との差分が最小になっていることがわかる。なお、最小位相差マップδ1上には位相とびが残っているが、それは、基準点x0の位置によって変わる。本実施形態では、基準点x0を図6aの点610に設定したため、位相とびが残る領域が存在する。ここで求まる2つの最小位相差マップδ1およびδ-1のうちいずれかが正しい静磁場不均一による位相差マップを表していることになる。次に、この最小位相差マップから、静磁場不均一による位相差マップを判定する。まず、図5ステップ504において、ステップ503で算出した2つの最小位相差マップの微分マップDおよびD-を算出する。微分をとる代わりに、式27、式28で示すように、隣り合う点との複素除算から求めてもよい。ただし、nは整数である。図5ステップ505において、算出した微分マップDおよびD-のばらつきを算出し、そのばらつきが最小の最小位相差マップを選択する。このばらつきは、たとえば標準偏差の大小から求めてもよいし、微分値のヒストグラムから検査者が判断してもよい。図8801に示すδ1は、水と脂肪の割合によって生じる位相とびのぶんだけ微分値のばらつきが増す。したがって、正しい位相差マップδは、802に示すδ-1となる。以上のステップによって、正しい静磁場不均一による位相差マップδが算出される。次に、上記ステップ514の位相補正画像算出方法、および、ステップ515の水画像および脂肪画像の算出方法について説明する。上記ステップ514では、第1エコーで計測した画像をI1、第2エコーで計測した画像をI2としたとき、上述した位相差マップδを使って画像I2を式29によって位相補正しI2’を算出する。なお、画像I1を式30によって位相補正しI1’を算出してもよい。以下、ここでは、ステップ514においてI2’を算出し、用いる場合を例にあげて説明する。次に、ステップ515において、画像I1とI2’とを使って水画像と脂肪画像を算出する。上述のように、第1エコー時間における水の位相をσw1、脂肪の位相をσF1、第2エコー時間における水の位相をσw2、脂肪の位相をσF2、水画像をW、脂肪画像をFとしたとき、I1とI2’は以下の連立方程式で表現できる。式33を、転置行列Atおよび逆行列を用いて変形し、得られた式35を用い、水画像および脂肪画像を算出する。ここで、本発明の効果を計算機シミュレーションおよびファントム実験にて確認する。計算機シミュレーション結果を図9に示す。図9aは水画像、図9bは脂肪画像である。なお、ここで使用したモデルは、上述した水と脂肪の割合が0:100、25:75、50:50、75:25、100:0となる領域が連続的に並んだモデルである。また、第1エコー時間を4.5ms、第2エコー時間を6.0msとした。図9a、bに示すように、設定した水と脂肪の割合に合わせて、水画像と脂肪画像とが分離できていることがわかる。次にファントム実験結果を図10に示す。図101001は塩化ニッケル水溶液ファントム、図101002は脂肪マットである。また、計測条件は、繰り返し時間TRが80ms、フリップ角は15度、撮影視野は360mm、計測マトリクスは128×128、積算回数は2回とした。図10aはエコー時間4.5msで計測した第1エコー画像、図10bはエコー時間6.0msで計測した第2エコー画像である。また、図10cおよびdは本発明によって分離された水画像および脂肪画像をそれぞれ示している。図10cおよびdに示すように、それぞれ水と脂肪とが同位相および逆位相になるエコー時間で取得した画像でなくても水画像と脂肪画像とが分離できていることがわかる。また、水画像と脂肪画像とが、間違えることなく判別されていることがわかる。本発明によって、以上のように、水脂肪の分離撮影において、水と脂肪との間の位相が、同位相および逆位相となるようなエコー時間で計測しなくても良好な水脂肪分離を実現する。さらに、得られた2つの分離画像ついて、水画像か脂肪画像かの判別が可能となる。"