"静磁場、高周波磁場および傾斜磁場をそれぞれ発生する磁場発生手段と、前記静磁場中に置かれた被検体から発生する磁気共鳴信号を検出する検出手段と、前記磁場発生手段および検出手段を制御する計測制御手段と、前記核磁気共鳴信号を用いて磁気共鳴スペクトル画像を作成し、表示する演算手段とを備え、前記検出手段が複数の小型コイルからなる受信コイルを有する磁気共鳴撮影装置であって、前記計測制御手段は、水からの核磁気共鳴信号を抑圧することなく核磁気共鳴スペクトルを計測する第一の計測シーケンス手段と、水からの核磁気共鳴信号を抑圧して核磁気共鳴スペクトルを計測する第二の計測シーケンス手段とを備え、前記演算手段は、前記第一の計測シーケンス手段により前記複数の小型コイル各々で計測された各画素の核磁気共鳴信号を用いて、前記複数の小型コイル各々の画素毎の信号位相値を算出し、前記第二の計測シーケンス手段により前記複数の小型コイル各々で計測した各画素の核磁気共鳴信号を前記信号位相値を用いて各々位相補正した後、前記小型コイル毎に作成される画像を合成することを特徴とする磁気共鳴撮影装置。請求項1記載の磁気共鳴撮影装置であって、前記演算手段は、前記第一の計測シーケンス手段により前記複数の小型コイル各々で計測した各画素の磁気共鳴信号を用いて、前記複数の小型コイル各々の画素毎に前記複数の小型コイル各々の信号雑音比に応じた加算係数を算出し、前記小型コイル毎に作成される画像を合成する際、前記加算係数を前記複数の小型コイル各々の画素毎に乗算することを特徴とする磁気共鳴撮影装置。請求項1記載の磁気共鳴撮影装置であって、前記演算手段は、前記位相補正後の合成により得られた各画素の核磁気共鳴信号から算出される合成後の各画素の核磁気共鳴スペクトルに対し、スペクトル軸方向の位相補正を行うことを特徴とする磁気共鳴撮影装置。請求項1記載の磁気共鳴撮影装置であって、前記計測制御手段は、前記第一の計測シーケンスにおける前記磁場発生手段が発生する傾斜磁場のうち少なくとも1軸の傾斜磁場であって位相エンコードに用いられる傾斜磁場の印加ステップ数が、前記第二の計測シーケンスにおける前記磁場発生手段が発生する前記少なくとも1軸の傾斜磁場の印加ステップ数よりも少ないことを特徴とする磁気共鳴撮影装置。請求項1記載の磁気共鳴撮影装置であって、前記計測制御手段は、前記第一の計測シーケンスにおいて、核磁気共鳴信号計測を繰り返す際の繰り返し時間を、前記第二の計測シーケンスで用いる繰り返し時間よりも短くすることを特徴とする磁気共鳴撮影装置。請求項5記載の磁気共鳴撮影装置であって、前記計測制御手段は、前記第一の計測シーケンスにおいて、前記磁場発生手段が発生する励起用の高周波磁場のフリップ角を、前記第二の計測シーケンスで用いるフリップ角よりも小さくすることを特徴とする磁気共鳴撮影装置。請求項1記載の磁気共鳴撮影装置であって、前記計測制御手段は、前記第一の計測シーケンスにおいて、前記被検体から発生する磁気共鳴信号を前記検出手段が検出し始める検出開始時刻を、前記第二の計測シーケンスにおける検出開始時刻よりも早くすることを特徴とする磁気共鳴撮影装置。請求項1記載の磁気共鳴撮影装置であって、前記計測制御手段は、前記第一の計測シーケンスにおいて、前記被検体から発生する磁気共鳴信号を前記検出手段が検出する検出時間を、前記第二の計測シーケンスで用いる検出時間よりも短くすることを特徴とする磁気共鳴撮影装置。静磁場中に置かれた検査対象に高周波磁場を少なくとも1回照射し、前記高周波磁場の照射の後に前記検査対象に傾斜磁場を少なくとも1回印加し、前記傾斜磁場の印加の後に前記検査対象から発生する磁気共鳴信号を受信コイルで検出し、検出した前記磁気共鳴信号から磁気共鳴スペクトル画像を算出する磁気共鳴スペクトル画像算出方法であって、前記受信コイルは複数の小型コイルLからなるマルチアレイコイルであり、2次元k空間もしくは3次元k空間,ky,kzの非水抑圧時のスペクトル情報を計測する第一の計測ステップと、前記第一の計測ステップにおいて、前記各小型コイルで得られたk空間の磁気共鳴信号Sw1,ky,kz方向のデータ番号を表すJ以下の整数に対して、各k空間軸方向にk空間成分と実空間成分とを変換し、2次元実空間もしくは3次元実空間,y,zの磁気共鳴信号を算出する第一のk空間-実空間変換ステップと、前記第一のk空間-実空間変換ステップで算出された実空間の磁気共鳴信号Sw2,y,zに対して、各小型コイルの実空間点毎に、先頭からN番目の信号Sw2,y,zの位相値φw,y,zを算出し、算出結果から補正値を求める位相値算出ステップと、2次元k空間もしくは3次元k空間,ky,kzの水抑圧時のスペクトル情報を計測する第二の計測ステップと、前記第二の計測ステップにおいて、前記各小型コイルで得られたk空間の磁気共鳴信号Sm1,ky,kzに対して、各k空間軸方向にk空間成分と実空間成分とを変換し、2次元実空間もしくは3次元実空間,y,zの磁気共鳴信号を算出する第二のk空間-実空間変換ステップと、前記第二のk空間-実空間変換ステップで算出された実空間の磁気共鳴信号Sm2,y,zに対して、各小型コイルの実空間点毎に、全点の位相値φm,y,zを前記位相値算出ステップで算出した補正値を用いて位相補正し、位相補正後の磁気共鳴信号Sm3,y,zを算出する位相値補正ステップと、前記位相値補正ステップで算出された各小型コイルの実空間点毎の磁気共鳴信号Sm3,y,zを、実空間点毎に複素数加算し、加算後の信号Sm4,y,zを、MACで得られた複数の画像を足し合わせるMAC合成を行った磁気共鳴信号とするMAC合成ステップと、前記MAC合成ステップで算出されたMAC合成後磁気共鳴信号の実空間点毎の磁気共鳴信号Sm4,y,zに対して、時間軸方向に時間成分と周波数成分とを変換し、周波数軸方向の情報を有する2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm5,y,zを算出する時間-周波数変換ステップと、を含むことを特徴とする磁気共鳴スペクトル画像算出方法。請求項9記載の磁気共鳴スペクトル画像算出方法であって、前記時間-周波数変換ステップで算出した周波数軸方向の情報を有する2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm5,y,zに対してスペクトル軸方向の位相補正を施し、2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm6,y,zを算出する第二の位相値補正ステップをさらに備えることを特徴とする磁気共鳴スペクトル画像算出方法。静磁場中に置かれた検査対象に高周波磁場を少なくとも1回照射し、前記高周波磁場の照射の後に、前記検査対象に傾斜磁場を少なくとも1回印加し、前記傾斜磁場の印加の後に、前記検査対象から発生する磁気共鳴信号を受信コイルで検出し、検出された前記磁気共鳴信号から磁気共鳴スペクトル画像を算出する磁気共鳴スペクトル画像算出方法であって、前記受信コイルは複数の小型コイルLからなるマルチアレイコイルであり、2次元k空間もしくは3次元k空間,ky,kzの非水抑圧時のスペクトル情報を計測する第一の計測ステップと、前記第一の計測ステップにおいて、前記各小型コイルで得られたk空間の磁気共鳴信号Sw1,ky,kz方向のデータ番号を表すJ以下の整数に対して、各k空間軸方向にk空間成分と実空間成分を変換する演算を施し、2次元実空間もしくは3次元実空間,y,zの磁気共鳴信号を算出する第一のk空間-実空間変換ステップと、前記第一のk空間-実空間変換ステップで算出された実空間の磁気共鳴信号Sw2,y,zに対して、各小型コイルの実空間点毎に、先頭からN番目の信号Sw2,y,zの絶対値Aw,y,zと位相値φw,y,zとを算出し、両算出結果をそれぞれ補正用絶対値および補正用位相値とする絶対値位相値算出ステップと、前記第一のk空間-実空間変換ステップで算出された実空間の磁気共鳴信号Sw2,y,zに対して、所定数の点の絶対値信号から、各小型コイルの実空間点毎に、ノイズ平均値Na,y,zとノイズ標準偏差Nd,y,zとを算出するノイズ算出ステップと、前記絶対値位相値算出ステップで算出した補正用絶対値と、前記ノイズ算出ステップで算出したノイズ平均値Na,y,zおよびノイズ標準偏差Nd,y,zとから、各小型コイルの実空間点毎に、信号雑音比SNR,y,zを算出する信号雑音比算出ステップと、2次元k空間もしくは3次元k空間,ky,kzの水抑圧時のスペクトル情報を計測する第二の計測ステップと、前記第二の計測ステップにおいて、前記各小型コイルで得られたk空間の磁気共鳴信号Sm1,ky,kzに対して、各k空間軸方向にk空間成分と実空間成分とを変換し、2次元実空間もしくは3次元実空間の磁気共鳴信号を算出する第2,y,zのk空間-実空間変換ステップと、前記第二のk空間-実空間変換ステップで算出された実空間の磁気共鳴信号Sm2,y,zに対して、各小型コイルの実空間点毎に、全点の位相値φm,y,zを前記絶対値位相値算出ステップで算出した補正用位相値を用いて位相補正し、位相補正後の磁気共鳴信号Sm3,y,zを算出する位相値補正ステップと、前記位相値補正ステップで算出された各小型コイルの実空間点毎の磁気共鳴信号Sm3,y,zに対して、前記信号雑音比算出ステップで算出した前記各小型コイルの実空間点毎の信号雑音比SNR,y,zを用いて算出した加算係数を乗算する信号雑音比乗算ステップと、前記信号雑音比乗算ステップで算出された各小型コイルの実空間点毎の磁気共鳴信号Sm4,y,zを、実空間点毎に複素数加算し、加算後の信号Sm5,y,zをMACで得られた複数の画像を足し合わせるMAC合成を行った磁気共鳴信号とするMAC合成ステップと、前記MAC合成ステップで算出されたMAC合成後磁気共鳴信号の実空間点毎の磁気共鳴信号Sm5,y,zに対して、時間軸方向に時間成分と周波数成分を変換する演算を施し、周波数軸方向の情報を有する2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm6,y,zを算出する時間-周波数変換ステップと、を含むことを特徴とする磁気共鳴スペクトル画像算出方法。請求項11記載の磁気共鳴スペクトル画像算出方法であって、前記ノイズ算出ステップでは、前記第一のk空間-実空間変換ステップで算出された実空間の磁気共鳴信号Sw2,y,zに対して、後部のM点分の絶対値信号を用いて前記ノイズ平均値Na,y,zと前記ノイズ標準偏差Nd,y,zとを算出することを特徴とする磁気共鳴スペクトル画像算出方法。請求項11記載の磁気共鳴スペクトル画像算出方法であって、前記時間-周波数変換ステップで算出した周波数軸方向の情報を有する2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm6,y,zに対してスペクトル軸方向の位相補正を施し、2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm7,y,zを算出する第二の位相値補正ステップをさらに備えることを特徴とする磁気共鳴スペクトル画像算出方法。請求項9記載の磁気共鳴スペクトル画像算出方法であって、前記第一の計測ステップにおいて各小型コイルで得られた磁気共鳴信号を用いて各小型コイルの感度分布を算出する感度算出ステップをさらに含み、前記MAC合成ステップに代えて、前記位相値補正ステップで算出された各小型コイルの実空間点毎の磁気共鳴信号と、前記感度算出ステップで算出された各小型コイルの感度分布とを用いて、画像の折り返しを除去する演算を行ない、各小型コイルで得られた複数の画像を合成する合成ステップを含むことを特徴とする磁気共鳴スペクトル画像算出方法。請求項11記載の磁気共鳴スペクトル画像算出方法であって、前記第一の計測ステップにおいて各小型コイルで得られた磁気共鳴信号を用いて各小型コイルの感度分布を算出する感度算出ステップをさらに含み、前記MAC合成ステップに代えて、前記位相値補正ステップで算出された各小型コイルの実空間点毎の磁気共鳴信号と、前記感度算出ステップで算出された各小型コイルの感度分布とを用いて、画像の折り返しを除去する演算を行ない、各小型コイルで得られた複数の画像を合成する合成ステップを含むことを特徴とする磁気共鳴スペクトル画像算出方法。MRIにおいて計測される磁気共鳴信号は、実部と虚部とから構成される。MRIにおけるMAC合成の方法には、各コイルで得られた画像信号をそのままで足し合わせる方法、各コイルで得られる画像信号の二乗和の平方根を足し合わせる方法、各コイルで得られた画像信号の位相補正後に実部毎、虚部毎で足し合わせる方法等がある。上記~の方式について、の方法では、合成後スペクトルの実部のみ/虚部のみを取り出すことが可能だが、各受信コイルの空間的な位相分布に不均一がある場合には合成効果が低下する。の方法では、空間的な位相不均一がある場合にも合成効果は低下しないが、合成後スペクトルの実部のみ/虚部のみを取り出すことは出来ない。の方法では、合成後スペクトルの実部のみ/虚部のみを取り出すことが可能であり、かつ、得られた信号のSNRが高い場合、空間的な位相不均一があっても合成効果は低下しない。スペクトル解析を行う際には、MAC合成後スペクトルの実部のみを表示する必要があるため、MAC合成後のスペクトルの実部のみを取り出すことが出来る必要がある。従って、これが可能であり、かつ、各受信コイルの空間的位相不均一があっても合成効果が低下しないの方法は、MRSIのMAC合成には優れている。しかし、の方法を用いたとしても、MRSI計測では、代謝物質の生体内濃度が非常に低いため、水抑圧計測で得られる信号のSNRは非常が低くなり、位相補正の補正精度が低下する場合がある。また、の方法では、位相補正の計算量がコイル数に比例して増加するため処理時間が長くなる。本発明は、上記事情に鑑みてなされたもので、MACを備えた磁気共鳴撮影装置を用いたMRSI計測において、MACの位相特性分布に空間的不均一がある場合にも、短い処理時間で高精度にMAC合成された磁気共鳴スペクトル画像を得ることを目的とする。本発明では、上記の方式を改良したMAC合成アルゴリズムを提供する。本発明によれば、水抑圧を行わない非水抑圧計測で取得されるSNRの高い非水抑圧画像信号を用いて、コイル毎および画素毎に、MAC合成のための位相歪みを揃えるための補正値を算出し、水抑圧を行って計測した本計測画像信号に対して、前記補正値を適用して、コイル毎および画素毎の位相を揃えた後で、信号の足し合わせ処理を行う。そして、合成後のスペクトル信号に対して、スペクトル軸方向の位相補正を行う。また、上記リファレンス計測を行う際、計測の繰り返し時間や繰り返し回数を減じることにより、計測時間を短縮する。また、各コイルで得られる信号を足し合わせる際、従来は足し合わせ係数1の等価合成が行われていた。これに対し、本発明では、上記リファレンス計測で取得される画像信号を用いてコイル毎および画素毎のSNRを算出し、本計測信号を足し合わせる際、前記コイル毎および画素毎のSNRから算出した係数を本計測信号に乗じた後に合成する。具体的には、静磁場、高周波磁場および傾斜磁場をそれぞれ発生する磁場発生手段と、前記静磁場中に置かれた被検体から発生する磁気共鳴信号を検出する検出手段と、前記磁場発生手段および検出手段を制御する計測制御手段と、前記核磁気共鳴信号を用いて磁気共鳴スペクトル画像を作成し、表示する演算手段とを備え、前記検出手段が複数の小型コイルからなる受信コイルを有する磁気共鳴撮影装置であって、前記計測制御手段は、水からの核磁気共鳴信号を抑圧することなく核磁気共鳴スペクトルを計測する第一の計測シーケンス手段と、水からの核磁気共鳴信号を抑圧して核磁気共鳴スペクトルを計測する第二の計測シーケンス手段とを備え、前記演算手段は、前記第一の計測シーケンス手段により前記複数の小型コイル各々で計測された各画素の核磁気共鳴信号を用いて、前記複数の小型コイル各々の画素毎の信号位相値を算出し、前記第二の計測シーケンス手段により前記複数の小型コイル各々で計測した各画素の核磁気共鳴信号を前記信号位相値を用いて各々位相補正した後、前記小型コイル毎に作成される画像を合成することを特徴とする磁気共鳴撮影装置を提供する。また、静磁場中に置かれた検査対象に高周波磁場を少なくとも1回照射し、前記高周波磁場の照射の後に前記検査対象に傾斜磁場を少なくとも1回印加し、前記傾斜磁場の印加の後に前記検査対象から発生する磁気共鳴信号を受信コイルで検出し、検出した前記磁気共鳴信号から磁気共鳴スペクトル画像を算出する磁気共鳴スペクトル画像算出方法であって、前記受信コイルは複数の小型コイルLからなるマルチアレイコイルであり、2次元k空間もしくは3次元k空間,ky,kzの非水抑圧時のスペクトル情報を計測する第一の計測ステップと、前記第一の計測ステップにおいて、前記各小型コイルで得られたk空間の磁気共鳴信号Sw1,ky,kz方向のデータ番号を表すJ以下の整数に対して、各k空間軸方向にk空間成分と実空間成分とを変換し、2次元実空間もしくは3次元実空間,y,zの磁気共鳴信号を算出する第一のk空間-実空間変換ステップと、前記第一のk空間-実空間変換ステップで算出された実空間の磁気共鳴信号Sw2,y,zに対して、各小型コイルの実空間点毎に、先頭からN番目の信号Sw2,y,zの位相値φw,y,zを算出し、算出結果から補正値を求める位相値算出ステップと、2次元k空間もしくは3次元k空間,ky,kzの水抑圧時のスペクトル情報を計測する第二の計測ステップと、前記第二の計測ステップにおいて、前記各小型コイルで得られたk空間の磁気共鳴信号Sm1,ky,kzに対して、各k空間軸方向にk空間成分と実空間成分とを変換し、2次元実空間もしくは3次元実空間,y,zの磁気共鳴信号を算出する第二のk空間-実空間変換ステップと、前記第二のk空間-実空間変換ステップで算出された実空間の磁気共鳴信号Sm2,y,zに対して、各小型コイルの実空間点毎に、全点の位相値φm,y,zを前記位相値算出ステップで算出した補正値を用いて位相補正し、位相補正後の磁気共鳴信号Sm3,y,zを算出する位相値補正ステップと、前記位相値補正ステップで算出された各小型コイルの実空間点毎の磁気共鳴信号Sm3,y,zを、実空間点毎に複素数加算し、加算後の信号Sm4,y,zを、MACで得られた複数の画像を足し合わせるMAC合成を行った磁気共鳴信号とするMAC合成ステップと、前記MAC合成ステップで算出されたMAC合成後磁気共鳴信号の実空間点毎の磁気共鳴信号Sm4,y,zに対して、時間軸方向に時間成分と周波数成分とを変換し、周波数軸方向の情報を有する2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm5,y,zを算出する時間-周波数変換ステップと、を含むことを特徴とする磁気共鳴スペクトル画像算出方法を提供する。本発明によれば、MACを備えた磁気共鳴撮影装置を用いたMRSI計測において、MACの位相特性分布に空間的不均一がある場合にも、短い処理時間で高精度にMAC合成された磁気共鳴スペクトル画像を得ることができる。<<第一の実施形態>>以下、本発明を適用した第一の実施の形態を図面を参照して説明する。図1~は、それぞれ本実施形態の磁気共鳴撮影装置の外観図である。図1はソレノイドコイルで静磁場を発生するトンネル型磁石を用いた水平磁場方式の磁気共鳴撮影装置であり、図1は開放感を高めるために磁石を上下に分離したハンバーガー型の垂直磁場方式の磁気共鳴撮影装置である。また、図1は、図1と同じトンネル型の磁気共鳴撮影装置であるが、磁石の奥行を短くし且つ斜めに傾けることによって、開放感を高めている。なお、本発明は、これら磁気共鳴撮影装置を含む公知の構造の磁気共鳴撮影装置に適用することができる。図2は、本実施形態の磁気共鳴撮影装置の構成を示すブロック図である。この磁気共鳴撮影装置は、被検体1が置かれる空間に、静磁場を発生する静磁場コイル2と、互いに直交する3方向の傾斜磁場を与えるための傾斜磁場コイル3と、被検体1に対し高周波磁場を照射する送信用高周波コイル5と、被検体1から発生する磁気共鳴信号を受信する受信用高周波コイル6とを備えている。また静磁場均一度を調整できるシムコイル4を備えている場合もある。静磁場コイル2は、図1に示した装置の構造に応じて、種々の形態のものが採用される。傾斜磁場コイル3及びシムコイル4は、それぞれ傾斜磁場用電源部12及びシム用電源部13により駆動される。送信コイル5が照射する高周波磁場は、送信機7により生成され、静磁場中に置かれた被検体に印加される。受信コイル6は、複数のコイルを平面的あるいは立体的に配置したマルチアレイコイルからなる。マルチアレイコイルとしては、MRIにおいて用いられている公知のマルチアレイコイルを用いることができる。図3にその一例を示す。図3に示すマルチアレイコイルは、4個の平面型コイルを円筒状に配置したもので、測定対象とする原子核の共鳴周波数に合わせた共振特性を備えている。図では省略されているが、この共振回路は電気容量CのコンデンサとインダクタンスLのコイルの組み合わせにより構成される。また必要に応じて、各小型コイル間の磁気的な結合状態を回避するためのデカップリング回路を設けてもよい。なお、図3に示した受信コイルの例では、小型コイルの数が4個であるが、小型コイルの数は複数であればよく4個に限定されない。図2では、送信コイル5と受信コイル6が別個に設けられた構成を示しているが、送信用と受信用を兼用する一つの高周波コイルのみを用いる構成であってもよい。すなわち、図3に示すようなマルチアレイコイルは送信用を兼ねることもできる。この受信コイル6が検出した磁気共鳴信号は、受信機8を通して計算機9に送られる。計算機9は、予め登録されているプログラム、または、ユーザからの指示に従って、磁気共鳴信号に対して様々な演算処理を行ないスペクトル情報や画像情報を生成する。本実施形態の磁気共鳴撮影装置では、受信コイル6を構成する各小型コイルが検出した磁気共鳴信号の補正計算やMAC合成に必要な補正値や合成係数の計算、および補正計算やMAC合成の演算を計算機9が行う。計算機9には、ディスプレイ10、記憶装置11、シーケンス制御装置14、入力装置15などが接続されており、上述した生成したスペクトル情報や画像情報をディスプレイ10に表示したり記憶装置11に記録したりする。入力装置15は、測定条件や演算処理に必要な条件などを入力するためのもので、これら測定条件等も必要に応じて記憶装置11に記録される。シーケンス制御装置14は、予め登録されているプログラム、または、ユーザからの指示に従って、傾斜磁場用電源部12、シム用電源部13、送信機7及び受信機9を制御する。制御のタイムチャートは撮影方法によって予め設定されており、記憶装置11に格納されている。本実施形態の磁気共鳴撮影装置では、シーケンス制御装置14は、水からの核磁気共鳴信号を抑圧しないでMRSI計測を行うパルスシーケンスと、水からの核磁気共鳴信号を抑圧してMRSI計測を行うパルスシーケンスとを備え、これら2種のパルスシーケンスを組み合わせて実行する。これらパルスシーケンスとしては公知のパルスシーケンスを採用できる。その一例を図4および図5に示す。図4は、MRSI計測用のパルスシーケンスの一例を示す図である。本実施形態では、このパルスシーケンスは非水抑圧計測および水抑圧計測の両方で実行される。また、図5は、水抑圧計測の前に実行されるプリパルスシーケンスの一例を示す図である。このMRSIパルスシ-ケンスに従って、シーケンス制御装置14は、以下の制御を行う。まず初めに、第1スラブ選択用の第1の傾斜磁場Gs1と90°パルスと呼ばれる第1の高周波磁場RF1とを同時に印加し、第1スラブ内の核磁化を励起状態にする。ここで、TEはエコー時間、TRは繰返し時間である。次に、RF1の照射からTE/4後に、第2スラブ選択用の第2の傾斜磁場Gs2と180°パルスと呼ばれる第2の高周波磁場RF2とを同時に印加し、RF1によって励起されていた第1スラブ内の核磁化のうち、第2スラブにも含まれる核磁化を180°反転する。さらに、RF2の照射からTE/2後に、第3スラブ選択用の第3の傾斜磁場Gs3と180°パルスと呼ばれる第3の高周波磁場RF3とを同時に印加し、RF2によって反転された第1スラブと第2スラブとの交差領域内にある核磁化のうち、第3スラブにも含まれる核磁化を再度180°反転する。上記の3組の高周波磁場及び傾斜磁場の印加により、RF3の照射からTE/4後の時点をエコ-タイムとする磁気共鳴信号Sig1が発生する。発生した磁気共鳴信号Sig1は、所定のサンプリング間隔でj点分計測される。そして、3次元の空間情報を付与する位相エンコード用の傾斜磁場Gp1、Gp2およびGp3の印加強度を段階的に変化させて、上記Sig1の計測を繰り返し行う。この時間軸方向の信号変化を有する磁気共鳴信号Sig1には、上述したケミカルシフトの情報が含まれているため、後述する通り、時間軸方向のフーリエ変換を施すことにより、磁気共鳴スペクトル信号を得ることができる。なお、Gs1の印加の直後に印加されるGs1’は、Gs1に対するリフェイズ用の傾斜磁場である。また、RF2の印加の前後で印加されるGd1とGd1’、Gd2とGd2’及びGd3とGd3’は、RF1の照射により励起された核磁化の位相は乱さず、RF2の照射により励起された核磁化をディフェイズするための傾斜磁場である。さらに、RF3の印加の前後で印加されるGd4とGd4’、Gd5とGd5’及びGd6とGd6’は、RF1の照射により励起された核磁化の位相は乱さず、RF3の照射によって励起された核磁化をディフェイズするための傾斜磁場である。上述した図4のパルスシーケンスを実行することにより、上記の3つのスラブが交差する領域に含まれる核磁化のみを選択的に励起することが出来る。そして、上記位相エンコード用傾斜磁場Gp1、Gp2およびGp3の印加強度を、Gp1はa段階、Gp2はb段階およびGp3はc段階、それぞれ変化させて励起と計測を繰り返してa×b×c個の磁気共鳴信号を検出し、周波数空間と呼ばれる計測空間を埋める一連のデータを計測する。なお、ここで、kxはk空間のx軸成分を意味し、aはkx軸方向の計測番号を表すA以下の整数を意味する。同様に、kyおよびkzはそれぞれk空間のy軸およびz軸成分を表し、bおよびcはそれぞれky軸およびkz方向の計測番号を表すBおよびC以下の整数を意味する。なお、2次元k空間の場合はA,B,Cのいずれか1つが1で他の2つは2以上の整数、3次元k空間の場合はA,B,Cのいずれもが2以上の整数となる。次に、得られたa×b×c個の磁気共鳴信号に対して、各k空間軸方向にk空間成分と実空間成分とを変換する逆フーリエ変換を施し、更に、時間軸方向に時間成分と周波数成分とを変換するフーリエ変換を施し、a×b×c個の各空間点から発生した磁気共鳴スペクトルを得る。一方、水抑圧計測を行う場合、シーケンス制御装置14は、前記図4のMRSIシーケンスに従った励起検出を行う直前に、図5に示す水信号を抑圧するためのプリパルスシーケンスに従った制御を行う。すなわち、水抑圧計測は、図5に示すパルスシーケンスと図4に示すMRSIシーケンスとを一組とするパルスシーケンスが実行される。プリパルスシーケンスでは、まず初めに、水分子にのみ含まれている核磁化を励起させるために、送信周波数Ftを水の共鳴周波数Fwに設定し、且つ励起周波数帯域ΔFtを水ピーク幅ΔFw程度に設定した高周波磁場RFw1の照射を行う。次に、励起状態にある水の核磁化の位相をバラバラにして、水の核磁化のベクトル和をゼロとするために、ディフェイズ用傾斜磁場Gdw1の印加を行う。更に水信号の抑圧効果を増すために、水励起用高周波磁場RFw1及びディフェイズ用傾斜磁場Gdw1と同様の高周波磁場及びディフェイズ用傾斜磁場の印加を複数回繰り返して行う。高周波磁場RFwには、狭帯域の励起周波数特性を有するガウス波形が用いられる場合が多い。図5に示す例は、ディフェイズ用傾斜磁場としてGx、Gy、Gzのうちいずれか1軸の傾斜磁場を印加する例である。しかし、Gx、Gy、Gzの3軸全ての傾斜磁場を同時に印加しても良いし、いずれか2軸を同時に印加しても構わない。そして、この水磁化の疑似飽和状態が続いている間に、図4のMRSIシーケンスを行い、微弱な代謝物質の信号を測定する。なお通常、水励起用高周波磁場RFwのフリップ角は90°前後に設定する場合が多いが、ディフェイズ用傾斜磁場Gdwについては、印加軸数や印加強度として様々な組合せや数値を用いることができる。なお、非水抑圧計測時にも全ての水励起用高周波磁場のRFwのフリップ角を0°に設定して図5に示すプリパルスシーケンスを図4のMRSIシーケンスの直前に行うことができる。このようにすることで、非水抑圧計測時のディフェイズ用傾斜磁場Gdwにより引き起こされる渦電流の影響を水抑圧計測時と同様にすることができる。すなわち、非水抑圧計測も水抑圧計測もともに図5に示すプリパルスシーケンスを図4のMRSIシーケンスの前に行うよう構成し、図5に示す全ての水励起用高周波磁場RFwのフリップ角を0°と90°との間で切り替えることによって、非水抑圧計測と水抑圧計測とを切り替える。また、生体内から検出できる代謝物質の信号は、非常に微弱である場合が多いため、得られるスペクトルのSNRを向上させることを目的に、同一条件での計測を複数回繰り返し行ない、得られた信号を足し合わせる処理を行ってもよい。次に、MRSIで行われる位相補正の概要について説明する。MRSIにおいて、位相検波により計測される磁気共鳴信号は、実部と虚部とから構成され、時刻tにおける磁気共鳴信号Fの信号位相φは以下ので表される。φ=tan-1/Reここで、tan-1はアークタンジェント関数、Imは複素数Xの虚部、Reは複素数Xの実部を表す。上記磁気共鳴信号Fに対して、時間軸方向に時間成分と周波数成分とを変換する演算を施すと、ケミカルシフト情報を有するスペクトル信号Sが得られる。位相補正の対象はこのスペクトル信号Sである。このスペクトル信号Sにおける第i番目のデータSの信号位相φは、以下ので表される。φ=tan-1/Reスペクトル信号Sの信号位相の補正すべき成分には、1次成分C1と0次成分C0とが含まれる。1次成分C1は、データの番号iに比例して変化する信号位相φの変動であり、ケミカルシフトの大きさに応じて生じる位相歪みである。また、0次成分C0は、データの番号iに依存しない信号位相φのオフセット変動成分である。両成分から成る位相歪みPは、以下ので表される。P=C1×iC0位相補正は、で求められる信号位相φから、で求められる位相歪みP分を減算することにより行う。本実施形態では、MAC合成前の位相補正については、影響の大きい位相歪みのオフセット変動成分に対してのみ行い、MAC合成後に1次成分まで補正する。また、位相歪みのオフセット変動成分は、SNRの高い非水抑圧計測でえられた磁気共鳴画像信号から得る。次に、本実施形態のMRSI計測を説明する。図6は本実施形態の計測手順の一例を示す図である。本実施形態では、まず、非水抑圧計測を行い、小型コイル毎の各画素の、位相歪みの0次成分をMAC合成用の位相補正値として得る。そして、得られたMAC合成用補正値を用いて、水抑圧計測で得られる小型コイル毎の各画素の磁気共鳴信号の位相を補正し、MAC合成を行う。MAC合成後、各画素毎の磁気共鳴信号に対し、スペクトル軸方向の位相補正を行う。以下の各ステップにおいて、シーケンスの制御はシーケンス制御装置14が、得られたデータの処理は計算機9が行う。まず、複数の小型コイルLから構成される受信コイルを用いて上述した非水抑圧計測を実行することにより、小型コイルL毎に、2次元k空間もしくは3次元k空間,ky,kzの非水抑圧時のスペクトル情報を取得する。ここで、取得した各小型コイルLのk空間のスペクトル情報を、Sw1,ky,kz方向のデータ番号を表すJ以下の整数とする。次に、ステップ701において得られた各小型コイルLのk空間の磁気共鳴信号Sw1,ky,kzに対して、各k空間軸方向にk空間成分と実空間成分とを変換する演算を施すことにより、非水抑圧時の2次元実空間もしくは3次元実空間,y,zの磁気共鳴信号Sw2,y,zを算出する。次に、ステップ702で算出された実空間の磁気共鳴信号Sw2,y,zに対して、各小型コイルLの実空間点毎に、信号位相値を算出する。信号位相値は、t軸方向のデータ数Jの信号の全て或いは複数を用いて算出しても良いが、信号値の絶対値が比較的大きい一つのデータを用いて算出してもよい。本実施形態では、小型コイルL毎に先頭からN点目の信号Sw2,y,zを用いて、それぞれの信号位相値φw,y,zを算出する。Nとして、信号の絶対値が大きくなる値を選ぶことが好ましいが、所定の範囲のいずれか1点を用いてもよい。ここで、Sw2,y,zの位相値φw,y,zは、以下のに従って導出する。φw,y,z=tan-1{Im,y,z/Re,y,z}なお、求めた小型コイルL毎の信号位相値φw,y,zは、本計測で計測する信号を補正する補正値、すなわち、MAC合成用補正値として記憶装置11に記憶する。この信号位相値φw,y,zは、上述のデータの番号iに依存しない信号位相φのオフセット変動成分C0に相当する。次に、上述した水抑圧計測を実行することにより、小型コイルL毎に、2次元k空間もしくは3次元k空間の水抑圧時のスペクトル情報を取得する。取得した各小型コイルLのk空間のスペクトル情報をSm1,ky,kzとする。次に、ステップ704において得られた各小型コイルLのk空間の磁気共鳴信号Sm1,ky,kzに対して、各k空間軸方向にk空間成分と実空間成分を変換する演算を施すことにより、水抑圧時の2次元実空間もしくは3次元実空間,y,zの磁気共鳴信号Sm2,y,zを算出する。次に、ステップ703で算出した非水抑圧信号の信号位相値φw,y,zを用いて、ステップ705で算出した実空間の磁気共鳴信号Sm2,y,zに対して、各小型コイルLの実空間点毎に、コイル間の位相を揃える位相補正を行い、各小型コイルの実空間点毎の磁気共鳴信号Sm3,y,zを得る。なお、この位相補正では、t軸方向のデータ数Jの信号の全てを対象とし、各データtの実部および虚部を、以下のおよびに従って導出する。Re,y,z=cos×π×φw,y,z/180×Re,y,z-sin×π×φw,y,z/180×Im,y,zIm,y,z=sin×π×φw,y,z/180×Re,y,zcos×π×φw,y,z/180×Im,y,zそして、ステップ706で算出した各小型コイルの実空間点毎の磁気共鳴信号Sm3,y,zを、以下のおよびに従って複素数加算し、加算後の実空間点毎の磁気共鳴信号Sm4,y,zを得る。さらに、ステップ707で算出した加算後の実空間点毎の磁気共鳴信号Sm4,y,zに対して、時間軸方向に時間成分と周波数成分とを変換する演算を施すことにより、周波数軸方向の情報を有する2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm5,y,zを算出する。そして、ステップ708で算出された2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm5,y,zの各画素の磁気共鳴スペクトルに対して、スペクトル軸方向の位相補正を施す。ここで、スペクトル軸方向の位相補正の手順について説明する。図7は、スペクトル信号の絶対値と位相値との一例を示す図である。図7において、横軸は周波数を、スペクトルデータの縦軸は信号強度の絶対値を、位相データの縦軸は位相値を示す。図7に示すように、位相補正の対象となるスペクトル信号Sには、いくつかの信号ピークが離散的に存在する。これらの信号ピーク以外はノイズ領域であり、ノイズ領域の位相成分は意味の無い値である。このため、位相歪みを算出する際は、まず初めに信号ピークの検出を行う。信号ピークは、多項式適合法等を用いて2次微分値および3次微分値を算出し、適当な閾値を設定することによりピーク判定を行い、検出する。次に、検出した信号ピークの位相値をに従って算出する。更に、この算出した離散的な位相値についてへのフィッティング処理を行い、近似直線の算出を行う。近時直線の算出には、例えば、最小自乗近似法等が用いられる。フィッティング処理で得られた近似直線の傾きが上記C1、切片が上記C0に相当する。これらの値を用いて、磁気共鳴スペクトル画像Sm5,y,zのj番目のデータの位相歪をP=C1×jC0を算出し、スペクトル軸方向の位相補正を行い、2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm6,y,zを算出する。そして、磁気共鳴スペクトルの2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm6,y,zの実部成分Re,y,zを抽出する。以上説明したように、本実施形態によれば、MACを用いてMRSI計測を行う場合、水抑圧を行わない第1の計測で取得されるSNRの高いリファレンス画像信号を用いて、小型コイル毎および画素毎の位相を揃えるための補正値を算出する。そして、水抑圧を伴う本計測画像信号を、小型コイル毎および画素毎にMAC合成用補正値で補正する。その後、MAC合成し、合成後の信号について、画素毎にスペクトル軸方向の位相補正を行う。すなわち、MAC合成を行う前に、小型コイル毎および画素毎に位相歪みの0次成分のみを揃えるための補正を行い、位相歪みの1次成分を含めたスペクトル軸方向の位相補正は、MAC合成後のスペクトル信号に対して画素毎に行う。位相歪みの0次成分は、信号位相のオフセット成分であり、データの番号に依存しないものである。本実施形態によれば、各受信コイルの空間的な位相分布に不均一がある場合に合成効果を低下させる影響の大きいオフセット成分のみ補正して合成し、スペクトル軸方向の位相補正は合成後に行う。本実施形態によれば、MAC合成の前に位相歪みの0次成分を補正するため、位相の空間的不均一がある場合でもMAC合成効果の高い信号が得られる。さらに、位相歪みの1次成分も含むスペクトル軸方向の位相補正は、MAC合成を行った後の1組の磁気共鳴スペクトル画像に対してのみ行えば良いため、この補正に要する計算量がコイル数に比例しない。従って、小型コイル毎にスペクトル軸方向の全ての位相補正を行う従来の手法に較べて、処理時間が短くてすむ。また、小型コイル毎に行う位相補正に使用する補正値として、水抑圧を行わないSNRの高い計測データを用いるため、小型コイル毎の位相補正の精度も高い。従って、短時間に高精度にMAC合成された磁気共鳴スペクトル画像を得ることが可能となる。<<第二の実施形態>>次に本発明の第二の実施の形態を説明する。本実施形態においても、第一の実施形態と同様に、非水抑圧計測で得られた小型コイル毎の位相補正値を用いて、その後に行われる水抑圧計測で得られる信号を小型コイル毎に位相補正する。ただし、本実施形態では、水抑圧計測による信号を足し合わせる際、各小型コイルのSNRに応じた異なる加算係数を乗算する。なお、第一の実施形態は、本実施形態の全ての小型コイルの加算係数を常に1にしたものといえる。以下、第一の実施形態と異なる点についてのみ説明する。以下、本実施形態のMRSI計測を説明する。図8は本実施形態の計測手順の一例を示す図である。以下の各ステップは、計算機9が行う処理である。本実施形態では、まず、複数の小型コイルLから構成される受信コイルを用いて、上述した非水抑圧計測を実行することにより、小型コイルL毎に、2次元k空間もしくは3次元k空間,ky,kzの非水抑圧時のスペクトル情報Sw1,ky,kz方向のデータ番号を表すJ以下の整数を取得する。次に、ステップ801において取得した各小型コイルLのk空間の磁気共鳴信号Sw1,ky,kzに対して、各k空間軸方向にk空間成分と実空間成分とを変換する演算を施すことにより、非水抑圧時の2次元実空間もしくは3次元実空間,y,zの磁気共鳴信号Sw2,y,zを算出する。次に、ステップ802で算出された実空間の磁気共鳴信号Sw2,y,zに対して、各小型コイルLの実空間点毎に、信号の絶対値と位相値とを算出する。絶対値と位相値は、t軸方向のデータ数Jの信号の全て或いは複数を用いて算出しても良いが、信号の絶対値が比較的大きい一つのデータを用いて算出してもよい。本実施形態では、小型コイルL毎に先頭からN点目の信号Sw2,y,zを用いて、それぞれの信号の絶対値Aw,y,zと位相値φw,y,zとを算出する。Nとして、信号の絶対値が大きくなる値を選ぶことが好ましいが、所定の範囲のいずれか1点を用いてもよい。ここで、Sw2,y,zの絶対値Aw,y,zと位相値φw,y,zとは、以下のおよびに従って導出する。なお、求めた小型コイルL毎のN点目の信号の絶対値Aw,y,zと位相値φw,y,zとは、本計測で計測する信号を補正する補正値としてそれぞれ記憶装置11に記憶される。また、ステップ802で算出された実空間の磁気共鳴信号Sw2,y,zに対して、各小型コイルLの実空間点毎に、ノイズ平均値とノイズ標準偏差とを求める。ノイズ平均値とノイズ標準偏差とは一般に絶対値Aw,y,zの小さい信号から求めることが好ましい。本実施形態では、信号Sw2,y,zを構成するデータのうち、最後尾のM点分の絶対値Aw,y,z~Aw,y,zから、下記のおよびに従い、ノイズ平均値Na,y,zとノイズ標準偏差Nd,y,zとを算出する。なお、上述したようにMとして、Sw2,y,zの絶対値Aw,y,zが小さくなる値を選ぶことが好ましいが、例えば、最後尾から16~32点の範囲で固定された値を用いてもよい。さらに、ステップ803で算出した信号の絶対値Aw,y,zと、ステップ804で算出したノイズ平均値Na,y,zおよびノイズ標準偏差Nd,y,zと、を用いて各小型コイルLの信号雑音比SNR,y,zを以下のにより算出する。SNR,y,z=,y,z-Na,y,z/Nd,y,zそして、各空間点について、各小型コイルLで検出したその空間点の信号雑音比のコイル間の比R,y,zの合計値を一定する係数,y,zを、小型コイルL毎に算出する。一例として、4個の小型コイルL、L、LおよびLで検出した空間点,y,zにおける信号雑音比が、下記の、、およびで示される場合についての計算式を以下に示す。SNR,y,z=eSNR,y,z=fSNR,y,z=gSNR,y,z=hこの場合、信号雑音比のコイル間の比R,y,zは、下記の、、およびで示される。R,y,z=e/max=ER,y,z=f/max=FR,y,z=g/max=GR,y,z=h/max=Hなお、maxは、e,f,g,hのうちの最大の値を取る関数である。空間点,y,zについて、信号雑音比のコイル間の比R,y,zの全小型コイルLの合計値を常にPとする係数W,y,zは、4個の小型コイルL、L、LおよびLのそれぞれについて下記の、、およびで示される。W,y,z=E×P/W,y,z=F×P/W,y,z=G×P/W,y,z=H×P/なお、算出した係数W,y,zは、以下で行う加算処理における加算係数として用いるため、記憶装置11に保存する。なお、Pの値は、1組の磁気共鳴スペクトル画像信号の処理において同一の値であれば、磁気共鳴スペクトル画像内の画素間のピクセル濃度値の比は一定となるため、どのような値を設定してもよい。例えば、常に、小型受信コイルの数を設定しても良い。次に、上述した水抑圧計測を実行することにより、小型コイル毎に、2次元k空間もしくは3次元k空間の水抑圧時のスペクトル情報Sm1,ky,kzを取得する。次に、ステップ807において、前記各小型コイルで得られたk空間の磁気共鳴信号Sm1,ky,kzに対して、各k空間軸方向にk空間成分と実空間成分を変換する演算を施すことにより、水抑圧時の2次元実空間もしくは3次元実空間,y,zの磁気共鳴信号Sm2,y,zを算出する。次に、ステップ803で算出した非水抑圧信号の位相値φw,y,zを用いて、ステップ808で算出された実空間の磁気共鳴信号Sm2,y,zに対して、各小型コイルの実空間点毎に、コイル間の位相を揃える位相補正を行い、各小型コイルの実空間点毎の磁気共鳴信号Sm3,y,zを得る。なお、この位相補正では、t軸方向のデータ数Jの信号の全てを対象とし、各データtの実部および虚部を、以下のおよびに従って導出する。Re,y,z=cos×π×φw,y,z/180×Re,y,z-sin×π×φw,y,z/180×Im,y,zIm,y,z=sin×π×φw,y,z/180×Re,y,zcos×π×φw,y,z/180×Im,y,zそして、ステップ809で算出した各小型コイルの実空間点毎の磁気共鳴信号Sm3,y,zについて、ステップ806で算出した加算係数W,y,zを乗算し、磁気共鳴信号Sm4,y,zを算出する。この処理では、各データtの実部および虚部を、以下のおよびに従って導出する。Re,y,z=Re,y,z×W,y,zIm,y,z=Re,y,z×W,y,zそして、ステップ810で算出された各小型コイルの実空間点毎の磁気共鳴信号Sm4,y,zを、以下のおよびに従って複素数加算し、実空間点毎の磁気共鳴信号Sm5,y,zを算出する。さらに、前記ステップ811で算出された加算後の実空間点毎の磁気共鳴信号Sm5,y,zに対して、時間軸方向に時間成分と周波数成分とを変換する演算を施し、周波数軸方向の情報を有する2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm6,y,zを算出する。そして、ステップ812で算出された2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm6,y,zの各画素の磁気共鳴スペクトルに対して、第一の実施形態同様、スペクトル軸方向の位相補正を施し、2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像の実部成分Re,y,zを算出する。以上説明したように、各コイルで得られる信号を足し合わせる際、従来は足し合わせ係数1の等価合成が行われていたが、本実施形態では、非水抑圧により行われるリファレンス計測で取得される画像信号からコイル毎および画素毎のSNRを算出し、水抑圧により行われる本計測で取得される信号の合成を、前記コイル毎および画素毎のSNRから算出した係数を乗じた後に行う。すなわち、本実施形態によれば、MACを用いてMRSI計測を行う場合において、水抑圧を行わない第1の計測で取得されるSNRの高いリファレンス画像信号を用いて、小型コイル毎および画素毎の位相を揃えるための位相補正値およびMAC合成時の加算係数を算出する。そして、水抑圧を伴う本計測画像信号について、MAC合成を行う前に算出した位相補正値で位相補正を行い、MAC合成を行う際は、算出した加算係数を適用する。すなわち、MAC合成を行う前には、位相歪の0次成分のみを揃えるための補正を行い、位相歪の1次成分の補正については、MAC合成後のスペクトル信号に対して行う。このため、本実施形態によれば、第一の実施の形態と同様に、位相の空間的不均一がある場合でもMAC合成効果の高い信号が得られるとともに、MAC合成を行った後の1組の磁気共鳴スペクトル画像に対してのみスペクトル軸方向の位相補正を行えば良いため、この補正に要する計算量がコイル数に比例せず、小型コイル毎にスペクトル軸方向の全ての位相補正を行う従来の手法に較べて、処理時間が短くてすむ。また、小型コイル毎に行う位相補正に使用する補正値として、水抑圧を行わない計測データを用いるため、正確な位相補正を行うことができる。さらに、加算係数を採用しているため、受信コイルを構成する小型コイルのSNRにばらつきがある場合にも、ばらつきによるMAC合成信号の劣化を防止することができる。従って、本実施形態においても、短時間に高精度にMAC合成された磁気共鳴スペクトル画像を得ることが可能となる。以上、本発明を適用した磁気共鳴撮影装置の各実施形態を説明した。本発明は上記の各実施形態に限定されず種々の変更や応用が可能である。例えば、受信コイルを構成する小型コイルの配列は、平面的な配置であってもよい。また各小型コイルの形状は平面に限定されず、検査部位に応じた曲面状に配置されていてもよい。また、上記の各実施形態では、第1の計測として非水抑圧画像計測を行い、第2の計測として水抑圧画像計測を行っている。しかし、第1の計測と第2の計測との計測順番を入れ替え、第1の計測として水抑圧画像計測を行い、第2の計測として非水抑圧画像計測を行っても良い。また、上記の各実施形態では、図4に示す領域選択スピンエコー型のパルスシーケンスを用いてMRSI計測を行う場合を例にあげて記述した。しかし、パルスシーケンスはこれに限られない。MRSI計測が可能なパルスシーケンスならば、全て同様の効果を得ることができる。例えば、3D-CSIや4D-CSIと呼ばれる領域選択を行わない計測シーケンスや、EPSIと呼ばれる振動傾斜磁場を用いた高速計測シーケンス等を用いた場合においても、同等の効果が得られる。また、上記各実施形態では、初めに、補正値を求めるためのリファレンス計測として、2次元k空間もしくは3次元k空間の非水抑圧時のスペクトル情報を取得している。しかし、図4に示すパルスシーケンスを用いてこの非水抑圧計測を実行すると、計測条件によっては計測時間が長くなるため、下記の方法により、計測時間の増大を軽減してもよい。例えば、上記位相エンコード用傾斜磁場Gp1、Gp2およびGp3の印加強度を、Gp1は16段階、Gp2は16段階、Gp3は1段階と、それぞれ変化させ、256回の計測を繰り返した場合、繰り返し時間を2秒とすると、リファレンス計測に要する時間は512秒となる。従って、非水抑圧計測時間と水抑圧計測時間を併せた総計測時間は2倍の1024秒となり、臨床に適用するにはやや長すぎる。そこで、リファレンス計測に要する時間を短縮するために、計測の繰り返し回数を減らす。すなわち、位相エンコード用傾斜磁場の印加強度のステップ幅は変えずに、ステップ数だけを減らす。例えば、上記例の位相エンコード用傾斜磁場Gp1、Gp2およびGp3の印加強度のステップ幅は変えずに、Gp1、Gp2のステップ数だけをそれぞれ半分に減らし、Gp1は8段階、Gp2は8段階、Gp3は1段階とすると、計測時間を1/4の128秒に短縮できる。ここで、非水抑圧計測の位相エンコードのステップ数を減らした計測データをそのままリファレンス画像として用る場合、非水抑圧計測画像の空間分解能と水抑圧計測画像の空間分解能が異なり、そのままでは上述の位相補正等を行うことができない。このため、通常、下記のゼロフィリングと呼ばれる手法を用いて、空間分解能を揃える。以下、ゼロフィリングの具体的な処理手順について、図9を用いて簡単に説明する。例えば、上述の計測条件と同様に水抑圧計測時の位相エンコードのステップ数が16×16×1である場合、非水抑圧計測時に位相エンコードのステップ幅を変えずに位相エンコードのステップ数を8×8×1に減じ、リファレンス計測に要する時間を1/4に短縮すると、再構成処理を経て得られる画像の画素サイズは4倍となる。図9に示すように、リファレンス計測で得た8×8×1のデータを、16×16×1のマトリクスの中央に配置し、周囲のデータ欠損部分にゼロデータを詰める。そして、各位相エンコード軸方向に逆フーリエ変換を施し、画素サイズを、位相エンコード数が16×16×1の場合と同一にする。なお、ゼロフィリング処理に関して、計測データを再構成マトリクスに配置する際に、計測データの真のゼロエンコード点を検出して中心を合わせる、計測データに空間分布上のハミング関数等を乗算することにより計測中心のデータが最大信号強度を有するように補正を行う、等の処理を加えてもよい。また、リファレンス計測に要する時間を短縮する別の方法として、計測の繰り返し時間を短くする方法がある。例えば、上記例の繰り返し時間2秒を半分に減じた場合、計測時間を1/2の256秒に短縮できる。しかし、単純に繰り返し時間を短くすると、縦緩和時間の影響で核磁化の初期状態に変化が生じ、正しい位相情報が得られなくなる場合がある。このような場合は、非水抑圧計測で用いる励起用高周波磁場のフリップ角を水抑圧計測のフリップ角より小さく設定する。こうして、繰り返し時間を短くしても、縦緩和時間の影響を減じることができる。また、エコー時間もしくは信号検出時間が長いため、繰り返し時間を単純に短くできない場合がある。このような場合には、非水抑圧計測で用いるエコー時間を水抑圧計測のエコー時間より短く設定する。すなわち、非水抑圧計測時の被検体から発生する磁気共鳴信号の検出開始時刻を水抑圧計測時より早める。こうして、繰り返し時間を短くすることが可能となる。また、非水抑圧計測で用いる信号検出時間を水抑圧計測の信号検出時間より短く設定し、繰り返し時間を短くすることができる。なお、非水抑圧計測で用いる信号検出時間を極端に短くした場合、非水抑圧計測信号の信号検出時間内での時間変化をリファレンスとして用いる渦電流補正を行うことができなくなる。このような場合は、水抑圧計測画像において、信号強度が最も強い代謝物質のピーク位置を揃えるといった後処理を追加し、渦電流補正と同等な補正を行う。また、マルチアレイコイルを備えた磁気共鳴撮影装置を用いたMRSI計測において、各受信コイルの感度分布画像を用いて、位相エンコード用傾斜磁場強度の印加ステップ幅の低下により生じる画像の折り返しを除去し、撮影時間を短縮するパラレルイメージングと呼ばれる計測に本発明を適用することができる。感度分布画像に水抑圧を行わないリファレンス計測で取得されるSNRの高い磁気共鳴スペクトル画像を用いればよい。例えば、2つの受信コイルを用いて上述した図4のパルスシーケンスによるMRSI計測を行う際に、X軸方向の位相エンコード用傾斜磁場Gp1の印加強度のステップ幅を倍にしてステップ数を半分に減らすと、計測時間を1/2に短縮できるもののX軸方向の撮影視野も1/2になる。撮像視野が1/2の画像は、被写体の折り返しの入った画像となる。このため、パラレルMRSI計測では、2つのコイルの感度分布の違いを利用して連立方程式を解き、折り返しを分離し、撮像視野が1/1の画像を再構成する。参照。。ここでは、感度分布画像に水抑圧を行わないリファレンス計測で取得されるSNRの高い磁気共鳴スペクトル画像を用いるため、折り返しの分離の精度を高めることができる。以下、第一の実施形態および第二の実施形態の実施例を示す。図3に示すマルチアレイコイルを実装した図1に示す型の静磁場強度1.5テスラの磁気共鳴撮影装置を用い、図5および図4のパルスシーケンスを実行し、MRSI計測を行った。このとき、対象核種はプロトン、測定対象物はN-アセチルアラニン水溶液ファントムとした。その結果のMAC合成後のスペクトル画像を図10に示す。図10は、第一の実施形態の図6の実施手順に従ってMAC合成した結果で、図10は、第二の実施形態の図8の実施手順に従ってMAC合成した結果である。両画像とも、概ね、選択励起領域内が均等な信号領域として描出されており、本来のファントム内水溶液の均質な濃度分布を反映した均一な輝度分布を有する合成画像が得られていることが分かる。ただ、図10の画像では、選択励起領域内の左下部分がやや高信号の領域として描出されているが、これは受信コイル間の受信感度のばらつきに起因するものと推測され、図10の画像では、この不均一も解消されていることが分かる。"