"被検体の所定の領域に含まれている代謝物の周波数情報を有する代謝物スペクトルを取得するための代謝物シーケンスと、前記所定の領域に含まれている水の周波数情報を有する水スペクトルを取得するための水シーケンスとを実行する磁気共鳴装置であって、前記代謝物シーケンスと前記水シーケンスとを有するシーケンス群をn回実行するスキャン手段と、n回のシーケンス群のうちのii=1~n番目のシーケンス群により得られた水スペクトルに基づいて、前記i番目のシーケンス群により得られた代謝物スペクトルの周波数を補正する補正手段と、を有する、磁気共鳴装置。前記シーケンス群は、前記代謝物シーケンスを複数有している、請求項1に記載の磁気共鳴装置。前記シーケンス群は、前記代謝物シーケンスを1個有している、請求項1に記載の磁気共鳴装置。前記スキャン手段は、水の共鳴周波数の基準値を求めるためのスキャンを実行する、請求項1~3のうちのいずれか一項に記載の磁気共鳴装置。前記水スペクトルに基づいて水の共鳴周波数を検出する周波数検出手段を有する、請求項1~4のうちのいずれか一項に記載の磁気共鳴装置。前記補正手段は、前記i番目のシーケンス群により得られた水スペクトルの水の共鳴周波数と、前記水の共鳴周波数の基準値とに基づいて、前記i番目のシーケンス群により得られた代謝物スペクトルの周波数を補正する、請求項4又は5に記載の磁気共鳴装置。前記補正手段は、周波数が補正された前記代謝物スペクトルの位相を補正する、請求項6に記載の磁気共鳴装置。被検体の所定の領域に含まれている代謝物の周波数情報を有する代謝物スペクトルを取得するための代謝物シーケンスと、前記所定の領域に含まれている水の周波数情報を有する水スペクトルを取得するための水シーケンスとを有するシーケンス群をn回実行する磁気共鳴装置のプログラムであって、n回のシーケンス群のうちのii=1~n番目のシーケンス群により得られた水スペクトルに基づいて、前記i番目のシーケンス群により得られた代謝物スペクトルの周波数を補正する補正処理、を計算機に実行させるためのプログラム。一般的に、代謝物の信号値は小さいので、代謝物のスペクトルを取得するためのシーケンスを繰り返し実行し、各シーケンスで得られた代謝物のスペクトルを加算している。スペクトルを加算することにより、SN比を大きくすることができる。しかし、代謝物の共鳴周波数は、勾配磁場を印加するときにコイルに発生する熱などが原因で、ドリフトする。したがって、代謝物のスペクトルに生じた周波数のドリフトを補正する必要があるが、周波数のドリフト量は時間とともに変化するので、スペクトルの周波数を十分な精度で補正することは難しいという問題がある。したがって、代謝物のスペクトルの周波数を十分な精度で補正できることが望まれている。本発明の第1の態様は、被検体の所定の領域に含まれている代謝物の周波数情報を有する代謝物スペクトルを取得するための代謝物シーケンスと、前記所定の領域に含まれている水の周波数情報を有する水スペクトルを取得するための水シーケンスとを実行する磁気共鳴装置であって、前記代謝物シーケンスと前記水シーケンスとを有するシーケンス群をn回実行するスキャン手段と、n回のシーケンス群のうちのii=1~n番目のシーケンス群により得られた水スペクトルに基づいて、前記i番目のシーケンス群により得られた代謝物スペクトルの周波数を補正する補正手段と、を有する、磁気共鳴装置である。本発明の第2の態様は、被検体の所定の領域に含まれている代謝物の周波数情報を有する代謝物スペクトルを取得するための代謝物シーケンスと、前記所定の領域に含まれている水の周波数情報を有する水スペクトルを取得するための水シーケンスとを有するシーケンス群をn回実行する磁気共鳴装置のプログラムであって、n回のシーケンス群のうちのii=1~n番目のシーケンス群により得られた水スペクトルに基づいて、前記i番目のシーケンス群により得られた代謝物スペクトルの周波数を補正する補正処理、を計算機に実行させるためのプログラムである。n回のシーケンス群のうちのii=1~n番目のシーケンス群により得られた水スペクトルに基づいて、前記i番目のシーケンス群により得られた代謝物スペクトルの周波数を補正している。したがって、各シーケンス群で生じた周波数のドリフト量に合わせて代謝物スペクトルの周波数が補正されるので、補正精度を向上させることができる。以下、発明を実施するための形態を説明するが、本発明は、以下の形態に限定されることはない。図1は、本発明の一形態の磁気共鳴装置の概略図である。磁気共鳴装置以下、MR装置と呼ぶ。MR:MagneticResonance100は、マグネット2、テーブル3、受信コイル4などを有している。マグネット2は、被検体12が収容されるボア21と、超伝導コイル22と、勾配コイル23と、送信コイル24とを有している。超伝導コイル22は静磁場を印加し、勾配コイル23は勾配磁場を印加し、送信コイル24はRFパルスを送信する。尚、超伝導コイル22の代わりに、永久磁石を用いてもよい。テーブル3は、クレードル3aを有している。クレードル3aは、ボア21内に移動できるように構成されている。クレードル3aによって、被検体12はボア21に搬送される。受信コイル4は、被検体12の頭部に取り付けられている。受信コイル4は、被検体12からの磁気共鳴信号を受信する。MR装置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から受け取ったデータに基づいて画像を再構成するなど、MR装置100の各種の動作を実現するように、MR装置100の各部の動作を制御する。中央処理装置9は、例えばコンピュータcomputerによって構成される。中央処理装置9は、周波数検出手段91、補正手段92、および加算手段93などを有している。周波数検出手段91は、水スペクトルに基づいて、水の共鳴周波数を検出する。補正手段92は、スペクトルの周波数の補正および位相の補正を行う。加算手段93は、スペクトルを加算する。中央処理装置9は、周波数検出手段91、補正手段92、および加算手段93の一例であり、所定のプログラムを実行することにより、これらの手段として機能する。操作部10は、オペレータ13により操作され、種々の情報を中央処理装置9に入力する。表示部11は種々の情報を表示する。MR装置100は、上記のように構成されている。尚、マグネット2、受信コイル4、送信器6、勾配磁場電源7、受信器8を合わせたものが、スキャン手段に相当する。図2は、被検体12の頭部の中の所定の領域Rを示す図、図3は、所定の領域Rから磁気共鳴信号を収集するためのスキャンの説明図である。本形態では、プリスキャンと、本スキャンが実行される。プリスキャンは、水の周波数情報共鳴周波数を有する水スペクトルを取得するためのスキャンである。本スキャンは、水スペクトルと、代謝物の周波数情報を有する代謝物スペクトルを取得するためのスキャンである。本スキャンでは、シーケンス群Gii=1~nが実行される。以下に、シーケンス群Giについて、説明する。尚、シーケンス群Giは、iの値に関わらず同じであるので、以下では、i=1、すなわち、シーケンス群G1を取り上げて説明する。シーケンス群G1は、代謝物シーケンスMS1~MS4と水シーケンスWSとを有している。以下に、代謝物シーケンスMS1~MS4および水シーケンスWについて、順に説明する。1代謝物シーケンスMS1~MS4について代謝物シーケンスMS1~MS4は、代謝物スペクトルを取得するためのシーケンスである。尚、代謝物シーケンスMS1~MS4は、同じシーケンスであるので、以下では、代表して代謝物シーケンスMS1を取り上げて説明する。代謝物シーケンスMS1は、RFパルスP1、P2、およびP3を有している。これらのRFパルスP1、P2、およびP3の送信周波数fmは、観測したい代謝物の化学シフトを考慮して設定される。例えば、代謝物として、NAANacetylLaspartateおよびコリンcholineを観測したい場合、RFパルスP1、P2、およびP3の送信周波数fmは、NAAの化学シフト≒2.00ppmとコリンの化学シフト≒3.2ppmとの中間値2.6ppmに対応する周波数に設定される。尚、頭部は、脳室などに大量の水を含んでいるので、代謝物の信号と一緒に水の信号も収集してしまうと、代謝物スペクトルのピークが水スペクトルのピークの裾野に隠れてしまい、代謝物を観測できなくなる。そこで、代謝物シーケンスは、水の信号を抑圧するための水抑圧用シーケンス部SUを有している。水抑圧用シーケンス部SUに含まれる水抑圧用RFパルスP0の送信周波数f0は、水の化学シフト≒4.7ppmに対応する周波数に設定されている。水抑圧用シーケンス部SUによって、水の信号を抑制することができ、所定の領域Rの代謝物スペクトルを観測することができる。2水シーケンスWSについて水シーケンスWSは、水スペクトルを取得するためのシーケンスである。水スペクトルは、代謝物スペクトルの位相を補正するために用いられる。一般的に、代謝物スペクトルは、渦電流などが原因で位相が変動するので、代謝物スペクトルの位相を補正するために、水スペクトルが用いられる。水スペクトルを用いて代謝物スペクトルの位相をどのように補正するかについては、例えば、文献KloseU,InvivoprotonspectroscopyinpresenceofeddycurrentsMagnResonMed1990Apr142630に記載されているので、詳細な説明は省略する。水シーケンスWSは、代謝物シーケンスMS1から水抑圧用シーケンス部SUを除くことにより得られるシーケンスである。水シーケンスWSは、水抑圧用シーケンス部SUが除かれているので、水シーケンスWSを実行することにより、水スペクトルを取得することができる。したがって、取得した水スペクトルを用いて、代謝物スペクトルの位相を補正することができる。以下に、本形態におけるMR装置100の動作について、図4のフローを参照しながら説明する。先ず、ステップST1において、プリスキャン図3参照を実行する。プリスキャンでは、本スキャンで使用される水シーケンスWSと同じシーケンスが実行される。したがって、プリスキャンを実行することによって、水の周波数情報共鳴周波数を有する水スペクトルを取得することができる。周波数検出手段91図1参照は、プリスキャンにより取得された水スペクトルを解析し、水の共鳴周波数fw0を検出する。検出された水の共鳴周波数fw0は、後述するステップST3において、水スペクトルの周波数のドリフト量を求めるための共鳴周波数の基準値として用いられる。プリスキャンが終了したら、ステップST2に進む。ステップST2では、先ず、ステップST1で得られた水の共鳴周波数fw0に基づいて、本スキャンで使用されるシーケンス代謝物シーケンスMS1~MS4および水シーケンスWSのRFパルスの送信周波数を設定する。具体的には、水抑圧用RFパルスP0の送信周波数f0は、f0=fW0に設定される。また、RFパルスP1、P2、およびP3の送信周波数fmは、fW0からΔfだけずれた値に設定される。Δfの値は、水の化学シフトと代謝物の化学シフトなどを考慮して、事前に設定された値である。このように、送信周波数を設定したら、本スキャンを実行する。図5は、図3に示すシーケンス群G1~Gnを実行することにより得られたスペクトルを概略的に示す図である。縦軸はスペクトルが取得された時間を表しており、横軸は周波数を表している。尚、図5には、ステップST1で得られた水の共鳴周波数の基準値fw0も示されている。図5において、シーケンス群Gii=1~nの代謝物シーケンスMS1~MS4により得られたスペクトルは、それぞれ、符号Mi1、Mi2、Mi3、Mi4で表されている。また、シーケンス群Giの水シーケンスWSにより得られたスペクトルは、符号Wiで表されている。例えば、シーケンス群G1の場合i=1、代謝物シーケンスMS1~MS4により得られたスペクトルは、符号M11、M12、M13、M14で表され、水シーケンスWSにより得られたスペクトルは、符号W1で表されている。図5を参照すると、時間とともに、スペクトルの周波数がドリフトすることがわかる。図5では、代謝物スペクトルM11~Mn4の周波数のドリフトの時間変化を破線Lで示してある。シーケンス群G1~Gnを実行した後、ステップST3に進む。ステップST3では、先ず、シーケンス群G1で取得されたスペクトルの周波数の補正処理を行う。以下に、この補正方法について説明する。補正手段92図1参照は、先ず、シーケンス群G1で取得された水スペクトルW1の周波数を補正する図6参照。図6は、水スペクトルW1の周波数を補正したときの様子を示す図である。周波数を補正する前の水スペクトルW1破線で示されているでは、水の共鳴周波数は、基準値fw0からΔfw1だけずれている。そこで、補正手段92は、水の共鳴周波数が基準値fw0に一致するように、水スペクトルW1の周波数を補正する。図5では、周波数補正後の水スペクトルを符号W1′で示してある。このようにして、水スペクトルW1の周波数が補正される。水スペクトルW1の周波数を補正したら、代謝物スペクトルM11~M14の周波数を補正する図7参照。図7は、代謝物スペクトルM11~M14の周波数を補正したときの様子を示す図である。補正手段92は、シーケンス群G1により得られた代謝物スペクトルM11~M14の周波数を、同じシーケンス群G1により得られた水スペクトルW1と同様に、Δfw1だけ補正する。図7では、周波数補正後の代謝物スペクトルを、符号M11′、M12′、M13′、M14′で示してある。このように、水スペクトルW1を用いて周波数のドリフト量Δfw1を求めることによって、代謝物スペクトルM11~M14の周波数のドリフトを補正することができる。次に、補正手段92は、シーケンス群G2により得られたスペクトルの周波数の補正処理を行う。補正手段92は、先ず、シーケンス群G2により得られた水スペクトルW2の周波数を補正する図8参照。図8は、水スペクトルW2の周波数を補正したときの様子を示す図である。周波数を補正する前の水スペクトルW2破線で示されているでは、水の共鳴周波数は、基準値fw0からΔfw2だけずれている。そこで、補正手段92は、水の共鳴周波数が基準値fw0に一致するように、水スペクトルW2の周波数を補正する。図8では、周波数補正後の水スペクトルを符号W2′で示してある。このようにして、水スペクトルW1の周波数が補正される。水スペクトルW2の周波数を補正したら、代謝物スペクトルM21~M24の周波数を補正する図9参照。図9は、代謝物スペクトルM21~M24の周波数を補正したときの様子を示す図である。補正手段92は、シーケンス群G2により得られた代謝物スペクトルM21~M24の周波数を、同じシーケンス群G2により得られた水スペクトルW2と同様に、Δfw2だけ補正する。図9では、周波数補正後の代謝物スペクトルを、符号M21′、M22′、M23′、M24′で示してある。このように、水スペクトルW2を用いて周波数のドリフト量Δfw2を求めることにより、代謝物スペクトルM21~M24の周波数のドリフトを補正することができる。以下同様に、補正手段92は、シーケンス群Giにより得られたスペクトルの周波数の補正処理を行う。例えば、i=n、つまり、シーケンス群Gnにより得られたスペクトルの補正処理を行うと、図10に示すように、周波数補正後の水スペクトルWn′と、周波数補正後の代謝物スペクトルMn1′~Mn4′とが得られる。このようにして、全ての代謝物スペクトルの周波数が補正される。図10には、2つの破線LおよびL′が示されている。破線Lは、周波数の補正をする前の代謝物シーケンスの周波数のドリフトの時間変化を表している。一方、破線L′は、周波数の補正をした後の代謝物シーケンスの周波数のドリフトの時間変化を表している。2つの破線LおよびL′を比較すると、代謝物スペクトルの周波数を補正することによって、代謝物スペクトルの周波数のドリフト量のばらつきが小さくなっていることがわかる。上記のようにして、全ての代謝物スペクトルの周波数を補正することができる。スペクトルの周波数を補正したら、ステップST4に進む。ステップST4では、補正手段92は、周波数補正後の代謝物スペクトルM11′~Mn4′の位相を補正する。補正手段92は、先ず、シーケンス群G1の周波数補正後の水スペクトルW1′に基づいて、同じシーケンス群G1の周波数補正後の代謝物スペクトルM11′~M14′の位相を補正する。次に、シーケンス群G2の周波数補正後の水スペクトルW2′に基づいて、同じシーケンス群G2の周波数補正後の代謝物スペクトルM21′~M24′の位相を補正する。以下同様に、シーケンス群Giの周波数補正後の水スペクトルWi′に基づいて、同じシーケンス群Giの周波数補正後の代謝物スペクトルMi1′~Mi4′の位相を補正する。例えば、i=nの場合、シーケンス群Gnの周波数補正後の水スペクトルWn′に基づいて、同じシーケンス群Gnの周波数補正後の代謝物スペクトルMn1′~Mn4′の位相を補正する。このようにして、全ての代謝物シーケンスの位相が補正される。図11に、位相補正された代謝物スペクトルを示す。図11では、位相補正後の代謝物スペクトルを符号M11′′~Mn4′′で示してある。位相補正をした後、ステップST5に進む。ステップST5では、加算手段93図1参照が、位相補正後の代謝物スペクトルM11′′~Mn4′′を加算する図12参照。図12は、位相補正後の代謝物スペクトルM11′′~Mn4′′を加算するときの説明図である。一般的に、代謝物の信号値は小さいので、代謝物スペクトルはSN比が小さくなる傾向がある。そこで、本形態では、加算手段93によって、位相補正後の代謝物スペクトルM11′′~Mn4′′を加算する。スペクトルを加算することによって、SN比の大きい代謝物スペクトルを得ることができる。図12では、位相補正後の代謝物スペクトルM11′′~Mn4′′を加算することにより得られた代謝物スペクトルを符号Maddで示してある。代謝物のペクトルMaddを得た後、代謝物スペクトルMaddを解析し、代謝物に関するデータを算出し、フローを終了する。本形態では、シーケンス群Giを実行することによって、代謝物スペクトルMi1~Mi4の他に、周波数のドリフト量を求めるための水スペクトルWiを取得している。そして、シーケンス群Giの水スペクトルWiにより求められた周波数のドリフト量を用いて、同じシーケンス群Giで得られた代謝物スペクトルMi1~Mi4の周波数を補正する。したがって、周波数のドリフト量が時間とともに変化しても、代謝物スペクトルの周波数を十分な精度で補正することができる。また、このように周波数補正をすることによって、代謝物スペクトルの周波数のドリフト量のばらつきを小さくすることができるので図10参照、ステップST5で得られた代謝物スペクトルMaddの半値幅の広がりを抑えることができる。尚、ステップST2を実行している間に、水スペクトルWiの共鳴周波数と、水の共鳴周波数の基準値fw0との周波数差を求め、求めた周波数差に基づいて、RFパルスの送信周波数を変更してもよい。このようにRFパルスの送信周波数を変更することによって、励起領域の位置ずれを補正することができ、良好な水抑制を行うこともできる。また、本形態では、図4のフローに従って代謝物スペクトルMaddを得ているが、図4とは別のフローに従って代謝物スペクトルMaddを得てもよい図13参照。図13は、別のフローを示す図である。ステップST1~ST3は、図4のフローと同じであるので説明は省略する。ステップST1~ST3を実行することによって、図10に示すように、周波数補正後の代謝物スペクトルM11′~Mn4′を得ることができる。ステップST3を実行したら、ステップST4に進む。ステップST4では、加算手段93図1参照が、周波数補正後の水スペクトルW1′~Wn′を加算し、更に、周波数補正後の代謝物スペクトルM11′~Mn4′を加算する。図14に、加算により得られた水スペクトルW′と、加算により得られた代謝物スペクトルM′を示す。スペクトルを加算した後、ステップST5に進む。ステップST5では、補正手段92が、水スペクトルW′を用いて、代謝物スペクトルM′の位相を補正する。これによって、位相補正後の代謝物スペクトルMaddを得ることができる。代謝物スペクトルMaddを得た後、代謝物スペクトルMaddを解析し、代謝物に関するデータを求め、フローを終了する。このように、スペクトルを加算した後で、位相を補正してもよい。尚、本形態では、シーケンス群Giは、図3に示すように、4つの代謝物シーケンスと1つの水シーケンスとを有している。しかし、シーケンス群Giは、図3に限定されることはなく、種々の変形が可能である図15参照。図15は、シーケンス群Giの変形例を説明する図である。図15では、シーケンス群Giは、水シーケンスWSの前にp個の代謝物シーケンスを有し、水シーケンスWSの後にq個の代謝物シーケンスを有している。図3では、p=q=2であったが、pおよびqは、2に限定されることはない。また、p=qである必要はなく、p≠qであってもよい。更に、p=0又はq=0にしてもよい。例えば、p=1、q=0の場合、シーケンス群Giは、以下のように表される。図16は、p=1、q=0の場合のシーケンス群Giを示す図である。図16では、1個の代謝物シーケンスMSごとに、周波数のドリフト量を求めるための水シーケンスWSが実行されるので、代謝物スペクトルの補正精度を更に向上させることができる。ただし、図16に示すシーケンス群を用いる場合、代謝物シーケンスの実行回数と同じ回数だけ水シーケンスを実行する必要があるので、スキャン時間が長くなる。したがって、スキャン時間の短縮化を優先したい場合は、1個のシーケンス群が、複数の代謝物シーケンスを有するようにすることが望ましい。尚、シーケンス群G1~Gnの各々が有する代謝物シーケンスの数は、必ずしも同じ数である必要はない図17参照。図17は、代謝物シーケンスの数が異なる場合の一例を示す図である。図17では、シーケンス群G1は、1個の代謝物シーケンスを有しているが、シーケンス群G2は、2個の代謝物シーケンスを有している。このように、代謝物シーケンスの数は異なっていてもよい。"