JP2013009900A - 磁気共鳴装置およびプログラム - Google Patents

Abstract

【課題】代謝物スペクトルの周波数を十分な精度で補正する。
【解決手段】水の中心周波数の基準値ｆｗ０を求めるためのプリスキャンを行った後、本スキャンを行う。本スキャンでは、代謝物シーケンスＭＳ１〜ＭＳ４と水シーケンスＷＳとを有するシーケンス群Ｇ_ｉ（ｉ＝１〜ｎ）を実行する。水シーケンスＷＳにより得られた水の中心周波数と、水の中心周波数との基準値ｆｗ０との差に基づいて、代謝物シーケンスＭＳ１〜ＭＳ４により得られた代謝物スペクトルＭ_ｉ１〜Ｍ_ｉ４の周波数を補正する。その後、周波数補正後の代謝物スペクトルの位相を補正し、加算する。
【選択図】図１０

Description

本発明は、代謝物スペクトルを取得するための代謝物シーケンスと、水スペクトルを取得するための水シーケンスとを実行する磁気共鳴装置、およびプログラムに関する。

近年、磁気共鳴装置を用いて、生体内の代謝物のスペクトルを取得し、医学的な診断に使用することが行われている（特許文献１参照）。

特開2001-346779号公報

一般的に、代謝物の信号値は小さいので、代謝物のスペクトルを取得するためのシーケンスを繰り返し実行し、各シーケンスで得られた代謝物のスペクトルを加算している。スペクトルを加算することにより、ＳＮ比を大きくすることができる。しかし、代謝物の共鳴周波数は、勾配磁場を印加するときにコイルに発生する熱などが原因で、ドリフトする。したがって、代謝物のスペクトルに生じた周波数のドリフトを補正する必要があるが、周波数のドリフト量は時間とともに変化するので、スペクトルの周波数を十分な精度で補正することは難しいという問題がある。

したがって、代謝物のスペクトルの周波数を十分な精度で補正できることが望まれている。

本発明の第１の態様は、被検体の所定の領域に含まれている代謝物の周波数情報を有する代謝物スペクトルを取得するための代謝物シーケンスと、前記所定の領域に含まれている水の周波数情報を有する水スペクトルを取得するための水シーケンスとを実行する磁気共鳴装置であって、
前記代謝物シーケンスと前記水シーケンスとを有するシーケンス群をｎ回実行するスキャン手段と、
ｎ回のシーケンス群のうちのｉ（ｉ＝１〜ｎ）番目のシーケンス群により得られた水スペクトルに基づいて、前記ｉ番目のシーケンス群により得られた代謝物スペクトルの周波数を補正する補正手段と、
を有する、磁気共鳴装置である。

本発明の第２の態様は、被検体の所定の領域に含まれている代謝物の周波数情報を有する代謝物スペクトルを取得するための代謝物シーケンスと、前記所定の領域に含まれている水の周波数情報を有する水スペクトルを取得するための水シーケンスとを有するシーケンス群をｎ回実行する磁気共鳴装置のプログラムであって、
ｎ回のシーケンス群のうちのｉ（ｉ＝１〜ｎ）番目のシーケンス群により得られた水スペクトルに基づいて、前記ｉ番目のシーケンス群により得られた代謝物スペクトルの周波数を補正する補正処理、
を計算機に実行させるためのプログラムである。

ｎ回のシーケンス群のうちのｉ（ｉ＝１〜ｎ）番目のシーケンス群により得られた水スペクトルに基づいて、前記ｉ番目のシーケンス群により得られた代謝物スペクトルの周波数を補正している。したがって、各シーケンス群で生じた周波数のドリフト量に合わせて代謝物スペクトルの周波数が補正されるので、補正精度を向上させることができる。

本発明の一形態の磁気共鳴装置の概略図である。 被検体１２の頭部の中の所定の領域Ｒを示す図である。 所定の領域Ｒから磁気共鳴信号を収集するためのスキャンの説明図である。 ＭＲ装置１００の動作フローを示す図である。 図３に示すシーケンス群Ｇ_１〜Ｇ_ｎを実行することにより得られたスペクトルを概略的に示す図である。 水スペクトルＷ_１の周波数を補正したときの様子を示す図である。 代謝物スペクトルＭ_１１〜Ｍ_１４の周波数を補正したときの様子を示す図である。 水スペクトルＷ_２の周波数を補正したときの様子を示す図である。 代謝物スペクトルＭ_２１〜Ｍ_２４の周波数を補正したときの様子を示す図である。 周波数補正後の水スペクトルＷ_ｎ′と、周波数補正後の代謝物スペクトルＭ_ｎ１′〜Ｍ_ｎ４′とを示す図である。 位相補正された代謝物スペクトルを示す図である。 位相補正後の代謝物スペクトルＭ_１１″〜Ｍ_ｎ４″を加算するときの説明図である。 ＭＲ装置の別の動作フローを示す図である。 加算により得られた水スペクトルＷ′と、加算により得られた代謝物スペクトルＭ′を示す図である。 シーケンス群Ｇ_ｉの変形例を説明する図である。 ｐ＝１、ｑ＝０の場合のシーケンス群Ｇ_ｉを示す図である。 代謝物シーケンスの数が異なる場合の一例を示す図である。

以下、発明を実施するための形態を説明するが、本発明は、以下の形態に限定されることはない。

図１は、本発明の一形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。

マグネット２は、被検体１２が収容されるボア２１と、超伝導コイル２２と、勾配コイル２３と、送信コイル２４とを有している。超伝導コイル２２は静磁場を印加し、勾配コイル２３は勾配磁場を印加し、送信コイル２４はＲＦパルスを送信する。尚、超伝導コイル２２の代わりに、永久磁石を用いてもよい。

テーブル３は、クレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１２はボア２１に搬送される。

受信コイル４は、被検体１２の頭部に取り付けられている。受信コイル４は、被検体１２からの磁気共鳴信号を受信する。

ＭＲ装置１００は、更に、シーケンサ５、送信器６、勾配磁場電源７、受信器８、中央処理装置９、操作部１０、および表示部１１を有している。

シーケンサ５は、中央処理装置９の制御を受けて、被検体１２を撮影するための情報を送信器６および勾配磁場電源７に送る。

送信器６は、シーケンサ５から送られた情報に基づいて、ＲＦコイル２４を駆動する駆動信号を出力する。

勾配磁場電源７は、シーケンサ５から送られた情報に基づいて、勾配コイル２３を駆動する駆動信号を出力する。

受信器８は、受信コイル４で受信された磁気共鳴信号を信号処理し、信号処理により得たれたデータを中央処理装置９に出力する。

中央処理装置９は、シーケンサ５および表示部１１に必要な情報を伝送したり、受信器８から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。中央処理装置９は、例えばコンピュータ（computer）によって構成される。中央処理装置９は、周波数検出手段９１、補正手段９２、および加算手段９３などを有している。

周波数検出手段９１は、水スペクトルに基づいて、水の共鳴周波数を検出する。
補正手段９２は、スペクトルの周波数の補正および位相の補正を行う。
加算手段９３は、スペクトルを加算する。

中央処理装置９は、周波数検出手段９１、補正手段９２、および加算手段９３の一例であり、所定のプログラムを実行することにより、これらの手段として機能する。

操作部１０は、オペレータ１３により操作され、種々の情報を中央処理装置９に入力する。表示部１１は種々の情報を表示する。

ＭＲ装置１００は、上記のように構成されている。尚、マグネット２、受信コイル４、送信器６、勾配磁場電源７、受信器８を合わせたものが、スキャン手段に相当する。

図２は、被検体１２の頭部の中の所定の領域Ｒを示す図、図３は、所定の領域Ｒから磁気共鳴信号を収集するためのスキャンの説明図である。
本形態では、プリスキャンと、本スキャンが実行される。

プリスキャンは、水の周波数情報（共鳴周波数）を有する水スペクトルを取得するためのスキャンである。

本スキャンは、水スペクトルと、代謝物の周波数情報を有する代謝物スペクトルを取得するためのスキャンである。本スキャンでは、シーケンス群Ｇ_ｉ（ｉ＝１〜ｎ）が実行される。以下に、シーケンス群Ｇ_ｉについて、説明する。尚、シーケンス群Ｇ_ｉは、ｉの値に関わらず同じであるので、以下では、ｉ＝１、すなわち、シーケンス群Ｇ_１を取り上げて説明する。

シーケンス群Ｇ_１は、代謝物シーケンスＭＳ１〜ＭＳ４と水シーケンスＷＳとを有している。以下に、代謝物シーケンスＭＳ１〜ＭＳ４および水シーケンスＷについて、順に説明する。

（１）代謝物シーケンスＭＳ１〜ＭＳ４について
代謝物シーケンスＭＳ１〜ＭＳ４は、代謝物スペクトルを取得するためのシーケンスである。尚、代謝物シーケンスＭＳ１〜ＭＳ４は、同じシーケンスであるので、以下では、代表して代謝物シーケンスＭＳ１を取り上げて説明する。代謝物シーケンスＭＳ１は、ＲＦパルスＰ_１、Ｐ_２、およびＰ_３を有している。これらのＲＦパルスＰ_１、Ｐ_２、およびＰ_３の送信周波数ｆ_ｍは、観測したい代謝物の化学シフトを考慮して設定される。例えば、代謝物として、ＮＡＡ（N-acetyl-L-aspartate）およびコリン（choline）を観測したい場合、ＲＦパルスＰ_１、Ｐ_２、およびＰ_３の送信周波数ｆ_ｍは、ＮＡＡの化学シフト（≒２．００ｐｐｍ）とコリンの化学シフト（≒３．２ｐｐｍ）との中間値２．６ｐｐｍに対応する周波数に設定される。

尚、頭部は、脳室などに大量の水を含んでいるので、代謝物の信号と一緒に水の信号も収集してしまうと、代謝物スペクトルのピークが水スペクトルのピークの裾野に隠れてしまい、代謝物を観測できなくなる。そこで、代謝物シーケンスは、水の信号を抑圧するための水抑圧用シーケンス部ＳＵを有している。水抑圧用シーケンス部ＳＵに含まれる水抑圧用ＲＦパルスＰ_０の送信周波数ｆ_０は、水の化学シフト（≒４．７ｐｐｍ）に対応する周波数に設定されている。

水抑圧用シーケンス部ＳＵによって、水の信号を抑制することができ、所定の領域Ｒの代謝物スペクトルを観測することができる。

（２）水シーケンスＷＳについて
水シーケンスＷＳは、水スペクトルを取得するためのシーケンスである。水スペクトルは、代謝物スペクトルの位相を補正するために用いられる。一般的に、代謝物スペクトルは、渦電流などが原因で位相が変動するので、代謝物スペクトルの位相を補正するために、水スペクトルが用いられる。水スペクトルを用いて代謝物スペクトルの位相をどのように補正するかについては、例えば、文献「Klose U，In vivo proton spectroscopy in presence of eddy currents. Magn Reson
Med. 1990 Apr;14(1):26-30」に記載されているので、詳細な説明は省略する。

水シーケンスＷＳは、代謝物シーケンスＭＳ１から水抑圧用シーケンス部ＳＵを除くことにより得られるシーケンスである。水シーケンスＷＳは、水抑圧用シーケンス部ＳＵが除かれているので、水シーケンスＷＳを実行することにより、水スペクトルを取得することができる。したがって、取得した水スペクトルを用いて、代謝物スペクトルの位相を補正することができる。

以下に、本形態におけるＭＲ装置１００の動作について、図４のフローを参照しながら説明する。

先ず、ステップＳＴ１において、プリスキャン（図３参照）を実行する。プリスキャンでは、本スキャンで使用される水シーケンスＷＳと同じシーケンスが実行される。したがって、プリスキャンを実行することによって、水の周波数情報（共鳴周波数）を有する水スペクトルを取得することができる。周波数検出手段９１（図１参照）は、プリスキャンにより取得された水スペクトルを解析し、水の共鳴周波数ｆｗ０を検出する。検出された水の共鳴周波数ｆｗ０は、後述するステップＳＴ３において、水スペクトルの周波数のドリフト量を求めるための共鳴周波数の基準値として用いられる。プリスキャンが終了したら、ステップＳＴ２に進む。

ステップＳＴ２では、先ず、ステップＳＴ１で得られた水の共鳴周波数ｆｗ０に基づいて、本スキャンで使用されるシーケンス（代謝物シーケンスＭＳ１〜ＭＳ４および水シーケンスＷＳ）のＲＦパルスの送信周波数を設定する。具体的には、水抑圧用ＲＦパルスＰ_０の送信周波数ｆ_０は、ｆ_０＝ｆ_Ｗ０に設定される。また、ＲＦパルスＰ_１、Ｐ_２、およびＰ_３の送信周波数ｆ_ｍは、ｆ_Ｗ０からΔｆだけずれた値に設定される。Δｆの値は、水の化学シフトと代謝物の化学シフトなどを考慮して、事前に設定された値である。このように、送信周波数を設定したら、本スキャンを実行する。

図５は、図３に示すシーケンス群Ｇ_１〜Ｇ_ｎを実行することにより得られたスペクトルを概略的に示す図である。縦軸はスペクトルが取得された時間を表しており、横軸は周波数を表している。尚、図５には、ステップＳＴ１で得られた水の共鳴周波数の基準値ｆｗ０も示されている。

図５において、シーケンス群Ｇ_ｉ（ｉ＝１〜ｎ）の代謝物シーケンスＭＳ１〜ＭＳ４により得られたスペクトルは、それぞれ、符号「Ｍ_ｉ１」、「Ｍ_ｉ２」、「Ｍ_ｉ３」、「Ｍ_ｉ４」で表されている。また、シーケンス群Ｇ_ｉの水シーケンスＷＳにより得られたスペクトルは、符号「Ｗ_ｉ」で表されている。例えば、シーケンス群Ｇ_１の場合（ｉ＝１）、代謝物シーケンスＭＳ１〜ＭＳ４により得られたスペクトルは、符号「Ｍ_１１」、「Ｍ_１２」、「Ｍ_１３」、「Ｍ_１４」で表され、水シーケンスＷＳにより得られたスペクトルは、符号「Ｗ_１」で表されている。

図５を参照すると、時間とともに、スペクトルの周波数がドリフトすることがわかる。図５では、代謝物スペクトルＭ_１１〜Ｍ_ｎ４の周波数のドリフトの時間変化を破線Ｌで示してある。シーケンス群Ｇ_１〜Ｇ_ｎを実行した後、ステップＳＴ３に進む。

ステップＳＴ３では、先ず、シーケンス群Ｇ_１で取得されたスペクトルの周波数の補正処理を行う。以下に、この補正方法について説明する。

補正手段９２（図１参照）は、先ず、シーケンス群Ｇ_１で取得された水スペクトルＷ_１の周波数を補正する（図６参照）。

図６は、水スペクトルＷ_１の周波数を補正したときの様子を示す図である。
周波数を補正する前の水スペクトルＷ_１（破線で示されている）では、水の共鳴周波数は、基準値ｆｗ０からΔｆｗ１だけずれている。そこで、補正手段９２は、水の共鳴周波数が基準値ｆｗ０に一致するように、水スペクトルＷ_１の周波数を補正する。図５では、周波数補正後の水スペクトルを符号「Ｗ_１′」で示してある。このようにして、水スペクトルＷ_１の周波数が補正される。

水スペクトルＷ_１の周波数を補正したら、代謝物スペクトルＭ_１１〜Ｍ_１４の周波数を補正する（図７参照）。

図７は、代謝物スペクトルＭ_１１〜Ｍ_１４の周波数を補正したときの様子を示す図である。
補正手段９２は、シーケンス群Ｇ_１により得られた代謝物スペクトルＭ_１１〜Ｍ_１４の周波数を、同じシーケンス群Ｇ_１により得られた水スペクトルＷ_１と同様に、Δｆｗ１だけ補正する。図７では、周波数補正後の代謝物スペクトルを、符号「Ｍ_１１′」、「Ｍ_１２′」、「Ｍ_１３′」、「Ｍ_１４′」で示してある。

このように、水スペクトルＷ_１を用いて周波数のドリフト量Δｆｗ１を求めることによって、代謝物スペクトルＭ_１１〜Ｍ_１４の周波数のドリフトを補正することができる。

次に、補正手段９２は、シーケンス群Ｇ_２により得られたスペクトルの周波数の補正処理を行う。補正手段９２は、先ず、シーケンス群Ｇ_２により得られた水スペクトルＷ_２の周波数を補正する（図８参照）。

図８は、水スペクトルＷ_２の周波数を補正したときの様子を示す図である。
周波数を補正する前の水スペクトルＷ_２（破線で示されている）では、水の共鳴周波数は、基準値ｆｗ０からΔｆｗ２だけずれている。そこで、補正手段９２は、水の共鳴周波数が基準値ｆｗ０に一致するように、水スペクトルＷ_２の周波数を補正する。図８では、周波数補正後の水スペクトルを符号「Ｗ_２′」で示してある。このようにして、水スペクトルＷ_１の周波数が補正される。

水スペクトルＷ_２の周波数を補正したら、代謝物スペクトルＭ_２１〜Ｍ_２４の周波数を補正する（図９参照）。

図９は、代謝物スペクトルＭ_２１〜Ｍ_２４の周波数を補正したときの様子を示す図である。
補正手段９２は、シーケンス群Ｇ_２により得られた代謝物スペクトルＭ_２１〜Ｍ_２４の周波数を、同じシーケンス群Ｇ_２により得られた水スペクトルＷ_２と同様に、Δｆｗ２だけ補正する。図９では、周波数補正後の代謝物スペクトルを、符号「Ｍ_２１′」、「Ｍ_２２′」、「Ｍ_２３′」、「Ｍ_２４′」で示してある。

このように、水スペクトルＷ_２を用いて周波数のドリフト量Δｆｗ２を求めることにより、代謝物スペクトルＭ_２１〜Ｍ_２４の周波数のドリフトを補正することができる。

以下同様に、補正手段９２は、シーケンス群Ｇ_ｉにより得られたスペクトルの周波数の補正処理を行う。例えば、ｉ＝ｎ、つまり、シーケンス群Ｇ_ｎにより得られたスペクトルの補正処理を行うと、図１０に示すように、周波数補正後の水スペクトルＷ_ｎ′と、周波数補正後の代謝物スペクトルＭ_ｎ１′〜Ｍ_ｎ４′とが得られる。このようにして、全ての代謝物スペクトルの周波数が補正される。

図１０には、２つの破線ＬおよびＬ′が示されている。破線Ｌは、周波数の補正をする前の代謝物シーケンスの周波数のドリフトの時間変化を表している。一方、破線Ｌ′は、周波数の補正をした後の代謝物シーケンスの周波数のドリフトの時間変化を表している。２つの破線ＬおよびＬ′を比較すると、代謝物スペクトルの周波数を補正することによって、代謝物スペクトルの周波数のドリフト量のばらつきが小さくなっていることがわかる。

上記のようにして、全ての代謝物スペクトルの周波数を補正することができる。スペクトルの周波数を補正したら、ステップＳＴ４に進む。

ステップＳＴ４では、補正手段９２は、周波数補正後の代謝物スペクトルＭ_１１′〜Ｍ_ｎ４′の位相を補正する。

補正手段９２は、先ず、シーケンス群Ｇ_１の周波数補正後の水スペクトルＷ_１′に基づいて、同じシーケンス群Ｇ_１の周波数補正後の代謝物スペクトルＭ_１１′〜Ｍ_１４′の位相を補正する。次に、シーケンス群Ｇ_２の周波数補正後の水スペクトルＷ_２′に基づいて、同じシーケンス群Ｇ_２の周波数補正後の代謝物スペクトルＭ_２１′〜Ｍ_２４′の位相を補正する。以下同様に、シーケンス群Ｇ_ｉの周波数補正後の水スペクトルＷ_ｉ′に基づいて、同じシーケンス群Ｇ_ｉの周波数補正後の代謝物スペクトルＭ_ｉ１′〜Ｍ_ｉ４′の位相を補正する。例えば、ｉ＝ｎの場合、シーケンス群Ｇ_ｎの周波数補正後の水スペクトルＷ_ｎ′に基づいて、同じシーケンス群Ｇ_ｎの周波数補正後の代謝物スペクトルＭ_ｎ１′〜Ｍ_ｎ４′の位相を補正する。このようにして、全ての代謝物シーケンスの位相が補正される。図１１に、位相補正された代謝物スペクトルを示す。図１１では、位相補正後の代謝物スペクトルを符号「Ｍ_１１″」〜「Ｍ_ｎ４″」で示してある。位相補正をした後、ステップＳＴ５に進む。

ステップＳＴ５では、加算手段９３（図１参照）が、位相補正後の代謝物スペクトルＭ_１１″〜Ｍ_ｎ４″を加算する（図１２参照）。

図１２は、位相補正後の代謝物スペクトルＭ_１１″〜Ｍ_ｎ４″を加算するときの説明図である。

一般的に、代謝物の信号値は小さいので、代謝物スペクトルはＳＮ比が小さくなる傾向がある。そこで、本形態では、加算手段９３によって、位相補正後の代謝物スペクトルＭ_１１″〜Ｍ_ｎ４″を加算する。スペクトルを加算することによって、ＳＮ比の大きい代謝物スペクトルを得ることができる。図１２では、位相補正後の代謝物スペクトルＭ_１１″〜Ｍ_ｎ４″を加算することにより得られた代謝物スペクトルを符号「Ｍ_ａｄｄ」で示してある。代謝物のペクトルＭ_ａｄｄを得た後、代謝物スペクトルＭ_ａｄｄを解析し、代謝物に関するデータを算出し、フローを終了する。

本形態では、シーケンス群Ｇ_ｉを実行することによって、代謝物スペクトルＭ_ｉ１〜Ｍ_ｉ４の他に、周波数のドリフト量を求めるための水スペクトルＷ_ｉを取得している。そして、シーケンス群Ｇ_ｉの水スペクトルＷ_ｉにより求められた周波数のドリフト量を用いて、同じシーケンス群Ｇ_ｉで得られた代謝物スペクトルＭ_ｉ１〜Ｍ_ｉ４の周波数を補正する。したがって、周波数のドリフト量が時間とともに変化しても、代謝物スペクトルの周波数を十分な精度で補正することができる。また、このように周波数補正をすることによって、代謝物スペクトルの周波数のドリフト量のばらつきを小さくすることができるので（図１０参照）、ステップＳＴ５で得られた代謝物スペクトルＭ_ａｄｄの半値幅の広がりを抑えることができる。

尚、ステップＳＴ２を実行している間に、水スペクトルＷ_ｉの共鳴周波数と、水の共鳴周波数の基準値ｆｗ０との周波数差を求め、求めた周波数差に基づいて、ＲＦパルスの送信周波数を変更してもよい。このようにＲＦパルスの送信周波数を変更することによって、励起領域の位置ずれを補正することができ、良好な水抑制を行うこともできる。

また、本形態では、図４のフローに従って代謝物スペクトルＭ_ａｄｄを得ているが、図４とは別のフローに従って代謝物スペクトルＭ_ａｄｄを得てもよい（図１３参照）。

図１３は、別のフローを示す図である。
ステップＳＴ１〜ＳＴ３は、図４のフローと同じであるので説明は省略する。ステップＳＴ１〜ＳＴ３を実行することによって、図１０に示すように、周波数補正後の代謝物スペクトルＭ_１１′〜Ｍ_ｎ４′を得ることができる。ステップＳＴ３を実行したら、ステップＳＴ４に進む。

ステップＳＴ４では、加算手段９３（図１参照）が、周波数補正後の水スペクトルＷ_１′〜Ｗ_ｎ′を加算し、更に、周波数補正後の代謝物スペクトルＭ_１１′〜Ｍ_ｎ４′を加算する。図１４に、加算により得られた水スペクトルＷ′と、加算により得られた代謝物スペクトルＭ′を示す。スペクトルを加算した後、ステップＳＴ５に進む。

ステップＳＴ５では、補正手段９２が、水スペクトルＷ′を用いて、代謝物スペクトルＭ′の位相を補正する。これによって、位相補正後の代謝物スペクトルＭ_ａｄｄを得ることができる。代謝物スペクトルＭ_ａｄｄを得た後、代謝物スペクトルＭ_ａｄｄを解析し、代謝物に関するデータを求め、フローを終了する。
このように、スペクトルを加算した後で、位相を補正してもよい。

尚、本形態では、シーケンス群Ｇ_ｉは、図３に示すように、４つの代謝物シーケンスと１つの水シーケンスとを有している。しかし、シーケンス群Ｇ_ｉは、図３に限定されることはなく、種々の変形が可能である（図１５参照）。

図１５は、シーケンス群Ｇ_ｉの変形例を説明する図である。
図１５では、シーケンス群Ｇ_ｉは、水シーケンスＷＳの前にｐ個の代謝物シーケンスを有し、水シーケンスＷＳの後にｑ個の代謝物シーケンスを有している。図３では、ｐ＝ｑ＝２であったが、ｐおよびｑは、２に限定されることはない。また、ｐ＝ｑである必要はなく、ｐ≠ｑであってもよい。更に、ｐ＝０又はｑ＝０にしてもよい。例えば、ｐ＝１、ｑ＝０の場合、シーケンス群Ｇ_ｉは、以下のように表される。

図１６は、ｐ＝１、ｑ＝０の場合のシーケンス群Ｇ_ｉを示す図である。
図１６では、１個の代謝物シーケンスＭＳごとに、周波数のドリフト量を求めるための水シーケンスＷＳが実行されるので、代謝物スペクトルの補正精度を更に向上させることができる。

ただし、図１６に示すシーケンス群を用いる場合、代謝物シーケンスの実行回数と同じ回数だけ水シーケンスを実行する必要があるので、スキャン時間が長くなる。したがって、スキャン時間の短縮化を優先したい場合は、１個のシーケンス群が、複数の代謝物シーケンスを有するようにすることが望ましい。

尚、シーケンス群Ｇ_１〜Ｇ_ｎの各々が有する代謝物シーケンスの数は、必ずしも同じ数である必要はない（図１７参照）。

図１７は、代謝物シーケンスの数が異なる場合の一例を示す図である。
図１７では、シーケンス群Ｇ_１は、１個の代謝物シーケンスを有しているが、シーケンス群Ｇ_２は、２個の代謝物シーケンスを有している。このように、代謝物シーケンスの数は異なっていてもよい。

２ マグネット
３ テーブル
３ａ クレードル
４ 受信コイル
５ シーケンサ
６ 送信器
７ 勾配磁場電源
８ 受信器
９ 中央処理装置
１０ 操作部
１１ 表示部
１２ 被検体
２１ ボア
２２ 超伝導コイル
２３ 勾配コイル
２４ 送信コイル
９１ 周波数検出手段
９２ 補正手段
９３ 加算手段
１００ ＭＲ装置

Claims (8)

1. 被検体の所定の領域に含まれている代謝物の周波数情報を有する代謝物スペクトルを取得するための代謝物シーケンスと、前記所定の領域に含まれている水の周波数情報を有する水スペクトルを取得するための水シーケンスとを実行する磁気共鳴装置であって、
   前記代謝物シーケンスと前記水シーケンスとを有するシーケンス群をｎ回実行するスキャン手段と、
   ｎ回のシーケンス群のうちのｉ（ｉ＝１〜ｎ）番目のシーケンス群により得られた水スペクトルに基づいて、前記ｉ番目のシーケンス群により得られた代謝物スペクトルの周波数を補正する補正手段と、
   を有する、磁気共鳴装置。
2. 前記シーケンス群は、前記代謝物シーケンスを複数有している、請求項１に記載の磁気共鳴装置。
3. 前記シーケンス群は、前記代謝物シーケンスを１個有している、請求項１に記載の磁気共鳴装置。
4. 前記スキャン手段は、水の共鳴周波数の基準値を求めるためのスキャンを実行する、請求項１〜３のうちのいずれか一項に記載の磁気共鳴装置。
5. 前記水スペクトルに基づいて水の共鳴周波数を検出する周波数検出手段を有する、請求項１〜４のうちのいずれか一項に記載の磁気共鳴装置。
6. 前記補正手段は、
   前記ｉ番目のシーケンス群により得られた水スペクトルの水の共鳴周波数と、前記水の共鳴周波数の基準値とに基づいて、前記ｉ番目のシーケンス群により得られた代謝物スペクトルの周波数を補正する、請求項４又は５に記載の磁気共鳴装置。
7. 前記補正手段は、
   周波数が補正された前記代謝物スペクトルの位相を補正する、請求項６に記載の磁気共鳴装置。
8. 被検体の所定の領域に含まれている代謝物の周波数情報を有する代謝物スペクトルを取得するための代謝物シーケンスと、前記所定の領域に含まれている水の周波数情報を有する水スペクトルを取得するための水シーケンスとを有するシーケンス群をｎ回実行する磁気共鳴装置のプログラムであって、
   ｎ回のシーケンス群のうちのｉ（ｉ＝１〜ｎ）番目のシーケンス群により得られた水スペクトルに基づいて、前記ｉ番目のシーケンス群により得られた代謝物スペクトルの周波数を補正する補正処理、
   を計算機に実行させるためのプログラム。