関心領域の各ボクセルごとに、水の共鳴周波数を含む第1の水スペクトルと、代謝物の共鳴周波数を含む第1の代謝物スペクトルとを取得するための第1のスキャンと、前記関心領域の各ボクセルごとに、水の共鳴周波数を含む第2の水スペクトルを取得するとともに、前記関心領域の全領域に対して、代謝物の共鳴周波数を含む第2の代謝物スペクトルを取得するための第2のスキャンと、を実行する手段と、前記第1の代謝物スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数を算出する算出手段と、前記第1のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数と、前記第1のスキャンにおける前記関心領域の各ボクセルの代謝物の共鳴周波数とに基づいて、前記第1のスキャンにおける前記関心領域の温度分布画像データを作成する温度分布画像データ作成手段と、を有し、前記温度分布画像データ作成手段は、A前記第1のスキャンにおける前記関心領域の温度分布画像データと、B前記第1のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数と、C前記第1のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、D前記第2のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数と、E前記第2のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、に基づいて、前記第2のスキャンにおける前記関心領域の温度分布画像データを作成する、磁気共鳴イメージング装置。前記温度分布画像データ作成手段は、A前記第1の水スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の各ボクセルごとに、水の共鳴周波数を検出し、B前記第1の代謝物スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の各ボクセルごとに、代謝物の共鳴周波数を検出し、C前記第2の水スペクトルに基づいて、前記第2のスキャンにおける前記関心領域の各ボクセルごとに、水の共鳴周波数を検出し、D前記第2の代謝物スペクトルに基づいて、前記第2のスキャンにおける前記関心領域の全領域に対して、代謝物の共鳴周波数を検出する、請求項1に記載の磁気共鳴イメージング装置。前記温度分布画像データ作成手段は、前記第1の水スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数を表す第1の共鳴周波数テーブルを作成し、前記第2の水スペクトルに基づいて、前記第2のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数を表す第2の共鳴周波数テーブルを作成する、請求項2に記載の磁気共鳴イメージング装置。前記温度画像分布データ作成手段は、前記第1の共鳴周波数テーブルと、前記第2の共鳴周波数テーブルと、前記第1のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、前記第2のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、に基づいて、前記関心領域の各ボクセルにおける温度の変化量を表す温度変化マップを作成し、前記第1のスキャンにおける前記関心領域の温度分布画像データに、前記温度変化マップを加算することによって、前記第2のスキャンにおける前記関心領域の温度分布画像データを作成する、請求項3に記載の磁気共鳴イメージング装置。関心領域の各ボクセルごとに、水の共鳴周波数を含む第1の水スペクトルと、代謝物の共鳴周波数を含む第1の代謝物スペクトルとを取得するための第1のスキャンと、前記関心領域の各ボクセルごとに、水の共鳴周波数を含む第2の水スペクトルを取得するとともに、前記関心領域の全領域に対して、代謝物の共鳴周波数を含む第2の代謝物スペクトルを取得するための第2のスキャンと、を実行する手段と、前記第1の代謝物スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数を算出する算出手段と、前記第2のスキャンにおける前記関心領域の温度分布画像データを作成する温度分布画像データ作成手段と、を有し、前記温度分布画像データ作成手段は、A前記第1のスキャンにおける前記関心領域の各ボクセルの代謝物の共鳴周波数と、B前記第1のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、C前記第2のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数と、D前記第2のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、に基づいて、前記第2のスキャンにおける前記関心領域の温度分布画像データを作成する、磁気共鳴イメージング装置。前記温度分布画像データ作成手段は、A前記第1の水スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の各ボクセルごとに、水の共鳴周波数を検出し、B前記第1の代謝物スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の各ボクセルごとに、代謝物の共鳴周波数を検出し、C前記第2の水スペクトルに基づいて、前記第2のスキャンにおける前記関心領域の各ボクセルごとに、水の共鳴周波数を検出し、D前記第2の代謝物スペクトルに基づいて、前記第2のスキャンにおける前記関心領域の全領域に対して、代謝物の共鳴周波数を検出する、請求項5に記載の磁気共鳴イメージング装置。前記温度分布画像データ作成手段は、前記第1の水スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数を表す第1の共鳴周波数テーブルを作成し、前記第2の水スペクトルに基づいて、前記第2のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数を表す第2の共鳴周波数テーブルを作成し、前記第1の代謝物スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の各ボクセルの代謝物の共鳴周波数を表す第3の共鳴周波数テーブルを作成する、請求項6に記載の磁気共鳴イメージング装置。前記温度画像分布データ作成手段は、前記第3の共鳴周波数テーブルと、前記第1のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、前記第2のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、に基づいて、前記第2のスキャンにおける前記関心領域の各ボクセルの代謝物の共鳴周波数の計算値を表す第4の共鳴周波数テーブルを作成し、前記第2の共鳴周波数テーブルに、前記第4の共鳴周波数テーブルを加算することにより、前記第2のスキャンにおける前記関心領域の温度分布画像データを作成する、請求項7に記載の磁気共鳴イメージング装置。前記第1のスキャンは、前記第1の水スペクトルを取得するための第1の水シーケンスと、前記第1の代謝物スペクトルを取得するための第1の代謝物シーケンスとを有し、前記第2のスキャンは、前記第2の水スペクトルを取得するための第2の水シーケンスと、前記第2の代謝物スペクトルを取得するための第2の代謝物シーケンスとを有する、請求項1~8のうちのいずれか一項に記載の磁気共鳴イメージング装置。前記第1の代謝物シーケンスおよび前記第2の代謝物シーケンスは、水信号を抑制するための水抑制部を有する、請求項9に記載の磁気共鳴イメージング装置。前記第1のスキャンにおいて、前記第1の水シーケンスの繰り返し時間は、前記第1の代謝物シーケンスの繰り返し時間に等しい、請求項10に記載の磁気共鳴イメージング装置。前記代謝物は、NAAである、請求項1~11のうちのいずれか一項に記載の磁気共鳴イメージング装置。水の共鳴周波数は、生体の温度によって変化するという特性がある。したがって、水の共鳴周波数の変化量を求めることによって、生体の温度を計測することが可能となる。しかし、水の共鳴周波数は、生体の温度だけでなく、様々な要因で変化する。生体の温度以外に水の共鳴周波数を変化させる大きな要因の一つとして、勾配磁場を印加するときにコイルに発生する熱がある。したがって、温度の時間変化の計測精度を高めるためには、勾配磁場を印加するときにコイルに発生する熱が原因で生じる水の共鳴周波数の変化量をできるだけ正確に知る必要がある。そこで、温度の時間変化の計測精度を高める方法として、以下のような方法が考えられる。図22および図23は、温度の計測精度を高める方法の一例の説明図である。図22は、温度の計測が行われる関心領域の一例を示す図である。関心領域Rは、n×nのボクセルによって規定されているとする。図22では、n=4の場合が示されているので、ボクセルの総数は、16個である。ボクセルは、符号V11~V44で示されている。尚、図22では、説明の便宜上、ボクセルの総数は16個であるが、16個に限られることは無く、ボクセルの総数は、必要に応じて、任意の値例えば64個に設定することができる。図23は、関心領域Rの温度を計測するときに実行されるスキャンを示す図である。スキャンAは、水シーケンスWおよび代謝物シーケンスMaを有している。水シーケンスWは、関心領域Rの各ボクセルごとに、水の共鳴周波数を含む水スペクトルを取得するためのシーケンスである。代謝物シーケンスMaは、関心領域Rの各ボクセルごとに、代謝物の共鳴周波数を含む代謝物スペクトルを取得するためのシーケンスである。水シーケンスWは、例えばPRESSである。代謝物シーケンスMaは、水シーケンスWに、水信号を抑制するための水抑制部WSを追加することによって得られるものである。関心領域Rをn×n個のボクセルとしてスキャンを実行する場合、水シーケンスWは、位相エンコードのステップを変えながら、繰り返し時間TRwでn×n回実行される。また、代謝物シーケンスMaは、位相エンコードのステップを変えながら、繰り返し時間TRmで、n×n回実行される。水シーケンスWをn×n回実行することによって、関心領域Rの各ボクセルごとに水スペクトルが得られ、代謝物シーケンスMaをn×n回実行することによって、関心領域Rの各ボクセルごとに代謝物スペクトルが得られる。したがって、各ボクセルごとに、水の共鳴周波数および代謝物の共鳴周波数が得られる。水の共鳴周波数は、関心領域Rの温度が原因で変化するだけでなく、勾配磁場を印加するときにコイルに発生する熱も原因となって変化する。一方、代謝物の共鳴周波数は、関心領域Rの温度によってはほとんど変化しないので、代謝物の共鳴周波数の変化の大部分は、勾配磁場を印加するときにコイルに発生する熱が原因であると考えることができる。したがって、各ボクセルごとに、水の共鳴周波数と代謝物の共鳴周波数との差を求めることによって、計測誤差の小さい温度分布画像データを作成することができる。以下同様に、2回目以降のスキャンAでも、1回目のスキャンAと同様に、水シーケンスをn×n回実行し、代謝物シーケンスをn×n回実行し、温度分布画像データを作成する。したがって、スキャンAを実行するたびに、関心領域Rの温度分布画像データを作成することができるので、関心領域Rの温度の時間変化を知ることができる。しかし、図23に示す方法では、1枚の温度分布画像データを取得する場合、水シーケンスをn×n回実行し、更に代謝物シーケンスをn×n回実行しなければならないので、スキャン時間がかかるという問題がある。例えば、n=4の場合、1回のスキャンAで、水シーケンスを4×4=16回実行し、代謝物シーケンスを4×4=16回実行する必要がある。したがって、水シーケンスの繰り返し時間および代謝物シーケンスの繰り返し時間が1秒の場合、1回のスキャンAで、16秒16秒=32秒の時間が必要となり、1枚の温度分布画像データを取得するのに32秒かかる。また、解像度を高くする場合、例えば、n=8とした場合、1回のスキャンAで、水シーケンスを8×8=64回実行し、代謝物シーケンスを8×8=64回実行する必要がある。したがって、水シーケンスの繰り返し時間および代謝物シーケンスの繰り返し時間を1秒とすると、1回のスキャンで、64秒64秒=128秒の時間が必要となり、1枚の温度分布画像データを取得するのに128秒かかる。本発明は、上記の事情に鑑み、短いスキャン時間で温度分布画像データを取得することが可能な磁気共鳴イメージング装置を提供することを目的とする。本発明の第1の態様は、関心領域の各ボクセルごとに、水の共鳴周波数を含む第1の水スペクトルと、代謝物の共鳴周波数を含む第1の代謝物スペクトルとを取得するための第1のスキャンと、前記関心領域の各ボクセルごとに、水の共鳴周波数を含む第2の水スペクトルを取得するとともに、前記関心領域の全領域に対して、代謝物の共鳴周波数を含む第2の代謝物スペクトルを取得するための第2のスキャンと、を実行する手段と、前記第1の代謝物スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数を算出する算出手段と、前記第1のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数と、前記第1のスキャンにおける前記関心領域の各ボクセルの代謝物の共鳴周波数とに基づいて、前記第1のスキャンにおける前記関心領域の温度分布画像データを作成する温度分布画像データ作成手段と、を有し、前記温度分布画像データ作成手段は、A前記第1のスキャンにおける前記関心領域の温度分布画像データと、B前記第1のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数と、C前記第1のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、D前記第2のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数と、E前記第2のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、に基づいて、前記第2のスキャンにおける前記関心領域の温度分布画像データを作成する、磁気共鳴イメージング装置である。本発明の第2の態様は、関心領域の各ボクセルごとに、水の共鳴周波数を含む第1の水スペクトルと、代謝物の共鳴周波数を含む第1の代謝物スペクトルとを取得するための第1のスキャンと、前記関心領域の各ボクセルごとに、水の共鳴周波数を含む第2の水スペクトルを取得するとともに、前記関心領域の全領域に対して、代謝物の共鳴周波数を含む第2の代謝物スペクトルを取得するための第2のスキャンと、を実行する手段と、前記第1の代謝物スペクトルに基づいて、前記第1のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数を算出する算出手段と、前記第2のスキャンにおける前記関心領域の温度分布画像データを作成する温度分布画像データ作成手段と、を有し、前記温度分布画像データ作成手段は、A前記第1のスキャンにおける前記関心領域の各ボクセルの代謝物の共鳴周波数と、B前記第1のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、C前記第2のスキャンにおける前記関心領域の各ボクセルの水の共鳴周波数と、D前記第2のスキャンにおける前記関心領域の全領域に対する代謝物の共鳴周波数と、に基づいて、前記第2のスキャンにおける前記関心領域の温度分布画像データを作成する、磁気共鳴イメージング装置である。第2のスキャンでは、関心領域の全領域に対して第2の代謝物スペクトルを取得すればよい。つまり、第2のスキャンでは、関心領域の各ボクセルごとに代謝物スペクトルを取得しなくてもよいので、スキャン時間を短縮することができる。以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることは無い。1第1の形態図1は、本発明の第1の形態の磁気共鳴イメージング装置の概略図である。磁気共鳴イメージング装置以下、MRI装置と呼ぶ。MRI:MagneticResonanceImaging100は、磁場発生装置2、テーブル3、受信コイル4などを有している。磁場発生装置2は、被検体12が収容されるボア21と、超伝導コイル22と、勾配コイル23と、送信コイル24とを有している。超伝導コイル22は静磁場B0を印加し、勾配コイル23は勾配磁場を印加し、送信コイル24はRFパルスを送信する。尚、超伝導コイル22の代わりに、永久磁石を用いてもよい。テーブル3は、クレードル31を有している。クレードル31は、ボア21内に移動できるように構成されている。クレードル31によって、被検体12はボア21に搬送される。受信コイル4は、被検体12からの磁気共鳴信号を受信する。MRI装置100は、更に、シーケンサ5、送信器6、勾配磁場電源7、受信器8、中央処理装置9、操作部10、および表示部11を有している。シーケンサ5は、中央処理装置9の制御を受けて、被検体12を撮影するための情報を送信器6および勾配磁場電源7に送る。送信器6は、シーケンサ5から送られた情報に基づいて、RFコイル24を駆動する駆動信号を出力する。勾配磁場電源7は、シーケンサ5から送られた情報に基づいて、勾配コイル23を駆動する駆動信号を出力する。受信器8は、受信コイル4で受信された磁気共鳴信号を信号処理し、中央処理装置9に出力する。中央処理装置9は、シーケンサ5および表示部11に必要な情報を伝送したり、受信器8から受け取った信号に基づいて画像を再構成するなど、MRI装置100の各種の動作を実現するように、MRI装置100の各部の動作を制御する。中央処理装置9は、例えばコンピュータcomputerによって構成される。中央処理装置9は、温度分布画像データ作成手段91および算出手段92などを有している。温度分布画像データ作成手段91は、スキャンAおよびB1~Bz後述する図2~図4参照により取得されたスペクトルに基づいて、各スキャンにおける関心領域の温度分布を表す温度分布画像データを作成する。算出手段92は、スキャンAの代謝物シーケンスMa後述する図3参照により取得された代謝物スペクトルに基づいて、スキャンAにおける関心領域の全領域に対する代謝物の共鳴周波数を算出する。操作部10は、オペレータ13により操作され、種々の情報を中央処理装置9に入力する。表示部11は種々の情報を表示する。MRI装置100は、上記のように構成されている。次に、被検体12を撮影するときに実行されるスキャンについて説明する。図2~図4は、第1の形態で実行されるスキャンの説明図である。図2は、第1の形態で実行されるスキャンの概略を示す図である。第1の形態では、スキャンAが1回だけ実行され、スキャンBii=1~zが繰り返し実行される。図3は、スキャンAの説明図である。スキャンAは、図23に示すスキャンAと同じであり、水シーケンスWおよび代謝物シーケンスMaを有している。水シーケンスWは、関心領域R図22参照の各ボクセルごとに、水の共鳴周波数を含む水スペクトルを取得するためのシーケンスである。代謝物シーケンスMaは、関心領域Rの各ボクセルごとに、代謝物の共鳴周波数を含む代謝物スペクトルを取得するためのシーケンスである。水シーケンスWは、例えばPRESSである。代謝物シーケンスMaは、水シーケンスWに、水信号を抑制するための水抑制部WSを追加することによって得られるものである。図4は、スキャンBii=1~zの説明図である。スキャンBiは、水シーケンスWおよび代謝物シーケンスMbを有している。水シーケンスWは、関心領域Rの各ボクセルごとに、水の共鳴周波数を含む水スペクトルを取得するためのシーケンスである。代謝物シーケンスMbは、関心領域Rの全領域に対して、代謝物の共鳴周波数を含む代謝物スペクトルを取得するためのシーケンスである。水シーケンスWは、図3に示す水シーケンスWと同じであり、n×n回実行される。しかし、代謝物シーケンスMbは、図3に示す代謝物シーケンスMaとは異なり、位相エンコードのための勾配磁場は与えないシーケンスであり、1回のみ実行される。代謝物シーケンスMbは、位相エンコードのための勾配磁場は与えないシーケンスであるので、n×nの各ボクセルごとに代謝物スペクトルを取得するのではなく、関心領域Rの全領域に対して代謝物スペクトルを取得するシーケンスである。尚、図4では、スキャンB1について説明されているが、他のスキャンB2~Bzについても、スキャンB1と同様に、水シーケンスWがn×n回実行され、代謝物シーケンスMbは1回のみ実行される。第1の形態では、スキャンAおよびB1~Bzを実行し、各スキャンが行われたときの関心領域Rの温度分布画像データを取得する。以下では、先ず、スキャンAが行われたときの関心領域Rの温度分布画像データを取得する方法について説明し、次に、スキャンB1~Bzが行われたときの関心領域Rの温度分布画像データを取得する方法について説明する。図5は、スキャンAにより得られたスペクトルを示す図である。スキャンAの水シーケンスWをn×n回実行することによって、関心領域Rの各ボクセルごとに水スペクトルWFが得られる。また、スキャンAの代謝物シーケンスMaをn×n回実行することによって、関心領域Rの各ボクセルごとに代謝物スペクトルMFaが得られる。図5では、説明の便宜上、n=4、つまり、関心領域Rのボクセル数が4×4=16個としている。したがって、水スペクトルWFおよび代謝物スペクトルMFaは、それぞれ16個得られる。そして、温度分布画像データ作成手段91図1参照は、水スペクトルWFに基づいて、スキャンAにおける関心領域Rの各ボクセルごとに、水の共鳴周波数を検出する。温度分布画像データ作成手段91は、更に、代謝物スペクトルMFaに基づいて、スキャンAにおける関心領域Rの各ボクセルごとに、代謝物の共鳴周波数を検出する図6参照。図6は、温度分布画像データ作成手段91により検出された共鳴周波数を示す図である。共鳴周波数テーブルwatermap0は、関心領域Rの各ボクセルの水の共鳴周波数を表している。共鳴周波数テーブルwatermap0において、各ボクセルの水の共鳴周波数は、fpap=1~4の整数、q=1~4の整数で表されている。また、共鳴周波数テーブルmetabmap0は、関心領域Rの各ボクセルの代謝物の共鳴周波数を表している。共鳴周波数テーブルmetabmap0において、各ボクセルの代謝物の共鳴周波数は、frsr=1~4の整数、s=1~4の整数で表されている。尚、代謝物スペクトルMFaには、複数の代謝物の共鳴周波数のピークが現れるが、全ての代謝物の共鳴周波数を検出する必要は無く、いずれか一つの代謝物の共鳴周波数を検出すればよい。一般的には、検出精度を高めるという観点から、スペクトルの中に大きなピークとして現れやすい代謝物の共鳴周波数を検出する。このような代謝物の共鳴周波数としては、例えば、NAANacetylLaspartateがある。そこで、第1の形態では、NAAの共鳴周波数を検出している。ただし、コリンなど、NAA以外の代謝物の共鳴周波数を検出してもよい。共鳴周波数テーブルwatermap0およびmetabmap0を作成した後、温度分布画像データ作成手段91図1参照は、共鳴周波数テーブルwatermap0およびmetabmap0に基づいて、関心領域Rの絶対温度の分布を表す温度分布画像データtempmap0を作成する図7参照。図7は、温度分布画像データtempmap0の一例を示す図である。共鳴周波数テーブルwatermap0が表す水の共鳴周波数fpaと、共鳴周波数テーブルmetabmap0が表す代謝物の共鳴周波数frsとの差は、温度の関数になる。したがって、温度分布画像データ作成手段91は、水の共鳴周波数fpaと代謝物の共鳴周波数frsとの差に基づいて、関心領域Rの絶対温度の分布を表す温度分布画像データtempmap0を作成することができる。温度分布画像データtempmap0の各ボクセルは、温度が高いほど白色に近い色で表されており、一方、温度が低いほど黒に近い色で表されている。温度分布画像データtempmap0によって、スキャンAにおける関心領域Rの温度分布を知ることができる。また、第1の形態では、関心領域Rの各ボクセルごとに得られた代謝物スペクトルMFaを加算し、スキャンAにおける関心領域Rの全領域に対するNAAの共鳴周波数を算出する図8参照。図8は、スキャンAにおける関心領域Rの全領域に対するNAAの共鳴周波数fM0を概略的に示す図である。第1の形態では、算出手段92図1参照が、関心領域Rの各ボクセルの代謝物スペクトルMFaを加算することにより加算スペクトルMFaddを求め、加算スペクトルMFaddから、スキャンAにおける関心領域Rの全領域に対するNAAの共鳴周波数fM0を決定する。NAAの共鳴周波数fM0は、スキャンBii=1~zおける温度分布画像データを作成するときに使用されるものである。スキャンBii=1~zにおける温度分布画像データを作成するときに、NAAの共鳴周波数fM0がどのように使用されるかについては、後述する。次に、スキャンBii=1~zにおける温度分布画像データを作成する手順について説明する。先ず、i=1、すなわち、スキャンB1における温度分布画像データを作成する手順について、図9~図14を参照しながら説明する。図9は、スキャンB1により得られた水スペクトルおよび代謝物スペクトルを概略的に示す図である。スキャンB1の水シーケンスWをn×n回実行することによって、関心領域Rの各ボクセルごとに水スペクトルWFが得られる。温度分布画像データ作成手段91は、水スペクトルWFに基づいて、スキャンB1における関心領域Rの各ボクセルごとに、水の共鳴周波数を検出し、水の共鳴周波数テーブルwatermap1を作成する。図9では、各ボクセルの水の共鳴周波数は、ftut=1~4の整数、u=1~4の整数で表されている。また、スキャンB1では、代謝物シーケンスMbは1回のみ実行される。代謝物シーケンスMbは、位相エンコードのための勾配磁場は与えないシーケンスであるので、代謝物シーケンスMbを実行することによって、関心領域Rの全領域に対して代謝物スペクトルMFbを取得することができる。温度分布画像データ作成手段91は、代謝物スペクトルMFbに基づいて、スキャンB1における関心領域Rの全領域に対するNAAの共鳴周波数fM1を求めることができる。共鳴周波数テーブルwatermap1およびNAAの共鳴周波数fM1を求めた後、スキャンB1における温度分布画像データを作成する。ただし、スキャンB1では、NAAの共鳴周波数fM1は、関心領域Rの各ボクセルごとに求められているわけではないので、水の共鳴周波数ftuとNAAの共鳴周波数fM1との差を求めるだけでは、十分な精度の温度分布画像データを求めることはできない。そこで、第1の形態では、以下のようにして、スキャンB1における温度分布画像データを作成する。第1の形態では、スキャンB1における温度分布画像データを作成する場合、先ず、スキャンAからスキャンB1までの間に生じた温度変化を表す温度変化マップを作成する。以下に、温度変化マップの作成方法について説明する。図10~図13は、温度変化マップの作成方法の一例の説明図である。図10は、スキャンAにより得られた共鳴周波数テーブルwatermap0、NAAの共鳴周波数fM0、および温度分布画像データtempmap0と、スキャンB1により得られた共鳴周波数テーブルwatermap1およびNAAの共鳴周波数fM1とを示す図である。第1の形態では、温度変化マップを作成するに当たり、先ず、温度分布画像データ作成手段91が、スキャンB1により得られた共鳴周波数テーブルwatermap1から、スキャンAにより得られた共鳴周波数テーブルwatermap0を減算する図11参照。図11は、共鳴周波数テーブルwatermap1から共鳴周波数テーブルwatermap0を減算することにより得られた周波数変化テーブルwatermapsub1を示す図である。周波数変化テーブルwatermapsub1は、関心領域Rの各ボクセルにおける水の共鳴周波数の変化量Δfvwv=1~4、w=1~4表している。水の共鳴周波数は、温度によって変化するという特性がある。したがって、周波数変化テーブルwatermapsub1の水の共鳴周波数の変化量Δfvwは、スキャンAからスキャンB1までの間に関心領域Rの温度がどれだけ変化したのかを求めるための指標として使用することができる。ただし、水の共鳴周波数は、関心領域Rの温度が原因で変化するだけでなく、勾配磁場を印加するときにコイルに発生する熱も原因となって変化する。したがって、周波数変化テーブルwatermapsub1が表す水の共鳴周波数の変化量Δfvwには、温度による変化分だけでなく、勾配磁場を印加するときにコイルに発生する熱が原因で生じる変化分も含まれている。つまり、水の共鳴周波数の変化量Δfvwは、以下の式1で近似することができる。Δfvw=ΔftempΔfMRI1ここで、Δftemp:関心領域Rの温度による水の共鳴周波数の変化分ΔfMRI:勾配磁場を印加するときにコイルに発生する熱が原因で生じる水の共鳴周波数の変化分したがって、周波数変化テーブルwatermapsub1だけでは、スキャンAからスキャンB1までの間に生じた温度変化を精度よく求めることができない。そこで、温度分布画像データ作成手段91は、スキャンAとスキャンB1との間に生じたNAAの共鳴周波数の変化量ΔfM1も求める図12参照。図12は、NAAの共鳴周波数の変化量ΔfM1を示す図である。NAAの共鳴周波数の変化量ΔfM1は、以下の式によって求められる。ΔfM1=fM1-fM02ここで、fM0:スキャンAにおけるNAAの共鳴周波数fM1:スキャンB1におけるNAAの共鳴周波数NAAが関心領域Rの中のどのボクセルに存在していても、NAAの共鳴周波数の変化量ΔfM1にはそれほど違いが無いことが実験的にわかっている。例えば、関心領域Rの中のボクセルV11図22参照に存在するNAAの共鳴周波数がΔfM1だけ変化した場合、他のボクセルV12~V44に存在するNAAの共鳴周波数も、ほぼΔfM1だけ変化することが実験的にわかっている。したがって、NAAの共鳴周波数の変化量ΔfM1は、関心領域Rの各ボクセルごとに別個に求める必要は無く、関心領域Rの全領域に対して1つの値ΔfM1だけ求めればよい。NAAの共鳴周波数の変化量ΔfM1を求めたら、温度分布画像データ作成手段91は、周波数変化テーブルwatermapsub1と、NAAの共鳴周波数の変化量ΔfM1とに基づいて、スキャンAからスキャンB1までの間に生じた温度変化を表す温度変化マップを作成する図13参照。図13は、温度変化マップtempsub1を示す図である。温度分布画像データ作成手段91は、周波数変化テーブルwatermapsub1の各ボクセルの水の共鳴周波数の変化量Δfvwから、NAAの共鳴周波数の変化量ΔfM1を減算する。ΔfvwとΔfM1との差分は温度変化の関数であるので、ΔfvwからΔfM1を減算することによって、温度変化マップtempsub1を作成することができる。温度変化マップtempsub1は、スキャンAからスキャンB1までの間に関心領域Rの各ボクセルに生じた温度変化量を表している。温度変化マップtempsub1の各ボクセルの値Δfvw′i=1~4、j=1~4は、以下の式で表される。Δfvw′=Δfvw-ΔfM13ここで、Δfvw:周波数変化テーブルwatermapsub1の各ボクセルの水の共鳴周波数の変化量ΔfM1:NAAの共鳴周波数の変化量NAAの共鳴周波数は、関心領域Rの温度によってはほとんど変化しないので、式3におけるNAAの共鳴周波数の変化量ΔfM1の大部分は、勾配磁場を印加するときにコイルに発生する熱が原因で生じる変化分と考えることができる。したがって、式3に従って、ΔfvwからΔfM1を減算することによって、Δfvwに含まれるΔfMRIの成分式1参照をほぼキャンセルすることができるので、誤差の少ない温度変化マップtempsub1を作成することができる。温度変化マップtempsub1を作成した後、温度分布画像データ作成手段91は、図14に示すように、温度変化マップtempsub1を、スキャンAにおける温度分布画像データtempmap0に加算する。このようにして、スキャンB1における温度分布画像データtempmap1が求められる。上記の説明では、スキャンB1における温度分布画像データtempmap1の作成方法について説明されている。しかし、スキャンB2~Bzについても、同様の手順で、温度分布画像データが作成される。例えば、スキャンBkにおける温度分布画像データtempmapkを作成する場合は、図15に示すように、先ず、スキャンBkにおける共鳴周波数テーブルwatermapkを作成し、NAAの共鳴周波数fMkを求める。そして、周波数変化テーブルwatermapsubkを作成し、NAAの共鳴周波数の変化量ΔfMkを求め、温度変化マップtempsubkを作成する。温度変化マップtempsubkを作成したら、温度変化マップtempsubkを温度分布画像データtempmap0に加算する。これにより、スキャンBkにおける温度分布画像データtempmapkが作成される。第1の形態では、スキャンBiにおいて、スキャンAのように関心領域Rの各ボクセルごとに代謝物スペクトルを取得する必要は無く、関心領域Rの全領域に対して代謝物スペクトルMFb図9参照を取得すればよい。したがって、スキャンBiで実行される代謝物シーケンスMbには、位相エンコードのための勾配磁場が不要となるので図4参照、スキャンBiは、スキャンAよりも、スキャン時間を大幅に短縮することができる。図16は、スキャンAとスキャンBiのスキャン時間の違いを説明する図である。例えば、関心領域Rを4×4のマトリクスと考えると、スキャンAでは、水シーケンスを16回、代謝物シーケンスを16回行うので、各シーケンスの繰り返し時間を1秒とすると、スキャン時間は32秒となる。一方、スキャンBiでは、水シーケンスを16回行うが、代謝物シーケンスは1回のみでよいので、各シーケンスの繰り返し時間を1秒とすると、スキャン時間は17秒となる。つまり、スキャンAについては32秒かかるが、スキャンA以降に実行されるスキャンBii=1~zについては、17秒ですむ。したがって、第1の形態の方法は、スキャンAを繰り返し実行する方法図23参照と比較して、撮影を高速に行うことができる。尚、第1の形態では、スキャンAは、水シーケンスWをn×n回実行した後で、代謝物シーケンスMaをn×n回実行している図3参照。しかし、水シーケンスWと代謝物シーケンスMaとを交互に実行してもよい図17参照。図17は、スキャンAにおいて、水シーケンスWと代謝物シーケンスMaとを交互に実行する場合の例を示す図である。水シーケンスWと代謝物シーケンスMaとを交互に実行することにより、動きアーチファクトを低減することができる。尚、スキャンAにおいて、代謝物シーケンスMaは、水抑制部WSが必要となるが、水シーケンスWは、水抑制部WSは不要である図3参照。したがって、水シーケンスWの繰り返し時間TRwは、代謝物シーケンスMaの繰り返し時間TRmよりも短くすることが可能である。ただし、水シーケンスWを実行しているときに流れる渦電流の影響と、代謝物シーケンスMaを実行しているときに流れる渦電流の影響を、できるだけ均一になるようにするためには、水シーケンスWの繰り返し時間TRwは、代謝物シーケンスMの繰り返し時間TRmと同じにすることが望ましい。また、第1の形態では、スキャンBiは、水シーケンスWをn×n回実行した後で、代謝物シーケンスMbを1回実行している図4参照。しかし、代謝物シーケンスMbは、水シーケンスWよりも先に実行してもよいし、水シーケンスWと次の水シーケンスWとの間に実行してもよい。2第2の形態第1の形態では、スキャンBiにおける温度分布画像データtempmapiを作成する場合、スキャンAにより得られた以下のi~iiiの組み合わせSet1が使用されている。i水の共鳴周波数テーブルwatermap0iiNAAの共鳴周波数fM0iii温度分布画像データtempmap0しかし、上記のiおよびiiiの代わりに、スキャンAにおいて得られた代謝物の共鳴周波数テーブルmetabmap0図6参照を用いてもよい。つまり、上記のi~iiiの組み合わせSet1の代わりに、以下のiiおよびivの組み合わせSet2を用いて、スキャンBiにおける温度分布画像データtempmapiを作成してもよい。iiNAAの共鳴周波数fM0iv代謝物の共鳴周波数テーブルmetabmap0第2の形態では、温度分布画像データ作成手段91が、i~iiiの組み合わせSet1の代わりに、iiおよびivの組み合わせSet2を用いて、スキャンBiにおける温度分布画像データtempmapiを作成する例について説明する。図18~図21は、上記の組み合わせSet2を用いてスキャンBiにおける温度分布画像データtempmapiを作成する一例の説明図である。尚、以下の説明では、i=1、すなわち、スキャンB1における温度分布画像データtempmap1を作成する一例について説明するが、i≠1の場合でも、同じ方法で温度分布画像データを求めることができる。図18は、スキャンAにより得られたNAAの共鳴周波数fM0および代謝物の共鳴周波数テーブルmetabmap0と、スキャンB1により得られた水の共鳴周波数テーブルwatermap1およびNAAの共鳴周波数fM1とを示す図である。第2の形態では、温度分布画像データ作成手段91は、スキャンAとスキャンB1との間に生じたNAAの共鳴周波数の変化量ΔfM1を求める図19参照。図19は、NAAの共鳴周波数の変化量ΔfM1を示す図である。NAAの共鳴周波数の変化量ΔfM1の求め方は、第1の形態と同じである。NAAの共鳴周波数の変化量ΔfM1を求めた後、NAAの共鳴周波数の変化量ΔfM1を、スキャンAにおけるNAAの共鳴周波数テーブルmetabmap0に加算する図20参照。図20は、スキャンAにおけるNAAの共鳴周波数テーブルmetabmap0にNAAの共鳴周波数の変化量ΔfM1を加算した後の様子を示す図である。NAAの共鳴周波数の変化量ΔfM1を、スキャンAにおけるNAAの共鳴周波数テーブルmetabmap0に加算することにより、スキャンB1におけるNAAの共鳴周波数テーブルmetabmap1′が得られる。NAAの共鳴周波数テーブルmetabmap1′は、スキャンB1における関心領域Rの各ボクセルのNAAの共鳴周波数の計算値Δfxy′x=1~4、y=1~4を表している。NAAの共鳴周波数テーブルmetabmap1′を作成したら、温度分布画像データ作成手段91は、図21に示すように、NAAの共鳴周波数テーブルmetabmap1′を、スキャンB1における水の共鳴周波数テーブルwatermap1に加算する。これにより、スキャンB1における温度分布画像データtempmap1が作成される。第2の形態においても、第1の形態と同様に、スキャンBiにおいて、スキャンAのように関心領域Rの各ボクセルごとに代謝物スペクトルを取得する必要は無く、関心領域Rの全領域に対して代謝物スペクトルMFb図9参照を取得すればよい。したがって、スキャンBiで実行される代謝物シーケンスMbには、位相エンコードのための勾配磁場が不要となるので図4参照、スキャンBiは、スキャンAよりも、スキャン時間を大幅に短縮することができる。