JP6411776B2 - 磁気共鳴装置 - Google Patents

Description

本発明は、スペクトルを作成する磁気共鳴装置に関する。

近年、人体の関心部位に含まれる代謝物を測定する方法として、代謝物の情報を含むスペクトルを計測するＭＲスペクトロスコピー（ＭＲSpectroscopy）法が知られている（特許文献１参照）。

特開２０１４−００４０１２号公報

ＭＲスペクトロスコピー法の一例として、複数のコイルエレメントを有するコイルを用いてＭＲ信号を収集するためのスキャンを実行し、収集したＭＲ信号に基づいて、代謝物の情報を含むスペクトルを作成する方法が知られている。以下に、この方法の一例について説明する。

先ず、スキャンが行われる部位（関心部位）を設定する。
図４９に、設定された関心部位を概略的に示す。図４９では、被検体の頭部に、４つの領域Ｖ１〜Ｖ４を有する関心部位が設定された例が示されている。また、図４９では、ＭＲ信号を収集するための複数のコイルエレメントも概略的に示されている。ここでは、簡単のため、２つのコイルエレメント４ａおよび４ｂのみが示されている。

関心部位を設定した後、水信号を抑制するためのＲＦパルス（例えば、CHESSパルス）を含むシーケンスを実行し、コイルエレメント４ａおよび４ｂを用いてＭＲ信号を収集する。そして、収集したＭＲ信号を用いて、コイルエレメントごとに、各領域におけるスペクトルを作成する。図５０に、作成されたスペクトルを概略的に示す。図５０（ａ）は、コイルエレメント４ａが受信したＭＲ信号に基づいて作成されたスペクトルＡ１〜Ａ４を示す。図５０（ｂ）は、コイルエレメント４ｂが受信したＭＲ信号に基づいて作成されたスペクトルＢ１〜Ｂ４を示す。水信号を抑制するためのＲＦパルスを含むシーケンスが実行されているので、図５０のスペクトルは、水のピークが低減されている。したがって、代謝物のピークを含むスペクトルを得ることができる。

コイルエレメントごとにスペクトルを得た後、これらのスペクトルを合成する。図５１にスペクトルの合成方法の一例を示す。

先ず、コイルエレメント４ａにより得られた領域Ｖ１のスペクトルＡ１から、所定の代謝物（例えば、ＮＡＡ、Ｃｈｏ、Ｃｒｅ）のピーク値ＱａおよびノイズＮａを求め、ノイズで正規化された代謝物のピーク値Ｑａ／Ｎａを、スペクトルＡ１に乗算する。これにより、スペクトルの強度が増幅されたスペクトルＭＡが得られる。

次に、コイルエレメント４ｂにより得られた領域Ｖ１のスペクトルＢ１から、所定の代謝物（例えば、ＮＡＡ、Ｃｈｏ、Ｃｒｅ）のピーク値ＱｂおよびノイズＮｂを求め、ノイズで正規化された代謝物のピーク値Ｑｂ／Ｎｂを、スペクトルＢ１に乗算する。これにより、スペクトルの強度が増幅されたスペクトルＭＢが得られる。

スペクトルＭＡおよびＭＢを得た後、スペクトルＭＡおよびＭＢを加算することにより、領域Ｖ１の合成スペクトルＭ１を得ることができる。

図５１では、領域Ｖ１における合成スペクトルの作成方法について示したが、他の領域Ｖ２、Ｖ３、およびＶ４についても、同様の方法で、合成スペクトルを作成することができる。

しかし、一般的に、代謝物は濃度が低いので、代謝物のピーク値はあまり大きくならない。したがって、上記の方法では、ＳＮ比の大きい合成スペクトルを得ることが難しいという問題がある。

したがって、ＳＮ比の大きい合成スペクトルを容易に得ることができる技術が望まれている。

本発明の第１の観点は、水および代謝物を含む関心部位からＭＲ信号を発生させるためのシーケンスを実行するスキャン手段であって、前記関心部位に含まれる水の信号を抑制するためのＲＦパルスを含む第１のシーケンスと、前記ＲＦパルスを含まない第２のシーケンスとを実行するスキャン手段と、
前記第１のシーケンスにより発生した第１のＭＲ信号と、前記第２のシーケンスにより発生した第２のＭＲ信号とを受信する複数のコイルエレメントを備えたコイルと、
各コイルエレメントで受信された第１のＭＲ信号に基づいて、前記コイルエレメントごとに、前記代謝物のピークを含む第１のスペクトルを生成するとともに、各コイルエレメントで受信された第２のＭＲ信号に基づいて、前記コイルエレメントごとに、前記水のピークを含む第２のスペクトルを生成するスペクトル生成手段と、
前記第２のスペクトルの水のピークの特性値に基づいて、前記コイルエレメントごとに得られた前記第１のスペクトルの強度を調整するための係数を算出する係数算出手段と、
前記係数を用いて前記第１のスペクトルの各々の強度を調整し、強度が調整された前記第１のスペクトルを合成する合成手段と、
を有する、磁気共鳴装置である。

本発明の第２の観点は、水および代謝物を含む関心部位からＭＲ信号を発生させるためのシーケンスを実行するスキャン手段であって、前記関心部位に含まれる水の信号を抑制するためのＲＦパルスを含む第１のシーケンスと、前記ＲＦパルスを含まない第２のシーケンスとを実行するスキャン手段と、
前記第１のシーケンスにより発生した第１のＭＲ信号と、前記第２のシーケンスにより発生した第２のＭＲ信号とを受信する複数のコイルエレメントを備えたコイルと、
各コイルエレメントで受信された第１のＭＲ信号に基づいて、前記コイルエレメントごとに、前記代謝物のピークを含む第１のスペクトルを生成するとともに、各コイルエレメントで受信された第２のＭＲ信号に基づいて、前記コイルエレメントごとに、前記水のピークを含む第２のスペクトルを生成するスペクトル生成手段と、
を有する磁気共鳴装置に適用されるプログラムであって、
前記第２のスペクトルの水のピークの特性値に基づいて、前記コイルエレメントごとに得られた前記第１のスペクトルの強度を調整するための係数を算出する係数算出処理と、
前記係数を用いて前記第１のスペクトルの各々の強度を調整し、強度が調整された前記第１のスペクトルを合成する合成処理と、
をコンピュータに実行させるためのプログラムである。

第２のスペクトルの水のピーク値の特性値に基づいて算出した係数を用いて第１のスペクトルを合成するので、ＳＮ比の大きい合成スペクトルを得ることができる。

本発明の第１の形態の磁気共鳴装置の概略図である。 受信コイル４の説明図である。 プロセッサ９が実行する処理を示す図である。 第１の形態で実行されるスキャンの説明図である。 ＭＲ装置の動作フローを示す図である。 ローカライザスキャンＡ０により得られた画像の一例を示す図である。 オペレータが設定した関心部位ＶＯＩの一例を示す図である。 ＣＳＩスキャンＡ１で実行されるシーケンスａ１を示す図である。 ＣＳＩスキャンＡ１で実行されるシーケンスａ２を示す図である。 ＣＳＩスキャンＡ１で実行されるシーケンスａ３を示す図である。 ＣＳＩスキャンＡ１で実行されるシーケンスａ４を示す図である。 シーケンスａ１〜ａ４を実行することにより得られたＭＲ信号を概略的に示す図である。 ＭＲ信号Ｓ１〜Ｓ４が受信コイル４で受信され、受信器で信号処理されるときの様子を概略的に示す図である。 ＭＲ信号Ｓ１、Ｓ２、Ｓ３、およびＳ４から得られた受信データを概略的に示す図である。 ＭＲＳスキャンＡ２で実行されるシーケンスの一例を示す図である。 ＭＲ信号Ｓ５が受信コイル４で受信され、受信器で信号処理されるときの様子を概略的に示す図である。 受信データをスペクトルに変換するときの説明図である。 位相補正後の代謝物スペクトルを示す図である。 スペクトルを概略的に示す図である。 ピーク値ＰａおよびＰｂを示す図である。 ノイズＮａおよびＮｂを示す図である。 水のピーク値をノイズで正規化するときの説明図である。 位相補正後の代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａに、係数Ｐａ／Ｎａを乗算するときの説明図である。 位相補正後の代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算するときの説明図である。 スペクトルを加算するときの説明図である。 第２の形態においてプロセッサ９が実行する処理を示す図である。 第２の形態におけるＭＲ装置の動作フローを示す図である。 位相を算出するときの説明図である。 位相データＳＰ５ａに基づいて、コイルエレメント４ａにより得られた受信データＳ１ａ〜Ｓ４ａの位相を補正するときの説明図である。 位相データＳＰ５ｂに基づいて、コイルエレメント４ｂにより得られた受信データＳ１ｂ〜Ｓ４ｂの位相を補正するときの説明図である。 スペクトルを概略的に示す図である。 位相補正後の代謝物スペクトルＭ１ａ〜Ｍ４ａに、係数Ｐａ／Ｎａを乗算するときの説明図である。 位相補正後の代謝物スペクトルＭ１ｂ〜Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算するときの説明図である。 スペクトルを加算するときの説明図である。 第３の形態で実行されるスキャンの説明図である。 第３の形態においてプロセッサ９が実行する処理を示す図である。 第３の形態におけるＭＲ装置の動作フローを示す図である。 ＭＲＳスキャンＡ３で実行されるパルスシーケンスの一例を示す図である。 ＭＲ信号Ｓ６が受信コイル４で受信され、受信器で信号処理されるときの様子を概略的に示す図である。 スペクトルを概略的に示す図である。 差分スペクトルＷａ＿difを示す図である。 差分スペクトルＷｂ＿difを示す図である。 ピーク値ＰａおよびＰｂを示す図である。 ノイズＮａおよびＮｂを示す図である。 水のピーク値をノイズで正規化するときの説明図である。 位相補正後の代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａに、係数Ｐａ／Ｎａを乗算するときの説明図である。 位相補正後の代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算するときの説明図である。 実測により得られたスペクトルを示す図である。 設定された関心部位を概略的に示す図である。 作成されたスペクトルを概略的に示す図である。 スペクトルの合成方法を示す図である。

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

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

マグネット２は、被検体１３が収容されるボア２１を有している。また、マグネット２には、静磁場を発生させるための超伝導コイル、勾配パルスを印加するための勾配コイル、およびＲＦパルスを送信するとともにＭＲ信号を受信するための送受信コイルなどが内蔵されている。

テーブル３は、被検体１３を支持するクレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１３はボア２１に搬送される。
受信コイル４は、被検体１３からの磁気共鳴信号を受信する。

図２は、受信コイル４の説明図である。図２（ａ）は、受信コイル４の斜視図、図２（ｂ）は、図２（ａ）のＡ−Ａ断面図である。
受信コイル４は、被検体からの磁気共鳴信号を受信するためのｎ個のコイルエレメントを有している。本形態では、説明の便宜上、ｎ＝２、すなわち、受信コイル４は２つのコイルエレメント４ａおよび４ｂを有する例が示されている。
図１に戻って説明を続ける。

ＭＲ装置１００は、更に、送信器５、勾配磁場電源６、受信器７、コンピュータ８、操作部１１、および表示部１２などを有している。

送信器５は、マグネット２に備えられた送受信コイルに電流を供給する。勾配磁場電源６はマグネット２に備えられた勾配コイルに電流を供給する。受信器７は、受信コイル４が受信した信号に対して、検波などの信号処理を行う。ＭＲ装置１００は、複数個の受信器７を有しているが、図１では、１個の受信器７のみが示されている。尚、マグネット２、送信器５、勾配磁場電源６を合わせたものがスキャン手段に相当する。また、受信器７は、データ生成手段に相当する。

コンピュータ８は、表示部１２に必要な情報を伝送したり、画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。コンピュータ８は、プロセッサ９およびメモリ１０などを有している。

メモリ１０には、プロセッサ９により実行されるプログラムなどが記憶されている。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出し、プログラムに記述されている処理を実行する。図３に、プロセッサ９が実行する処理を示す。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出すことにより、スペクトル作成手段９１〜画像作成手段９４などを構成する。

スペクトル作成手段９１は、受信器７から得られたデータに基づいて、スペクトルを作成する。スペクトル作成手段９１は、変換手段９１ａおよび補正手段９１ｂを有している。変換手段９１ａは、受信器７から得られたデータをスペクトルに変換する。補正手段９１ｂは、変換手段９１ａにより得られたスペクトルの位相を補正する。
係数算出手段９２は、スペクトルの強度を調整するための係数を算出する。
合成手段９３は、係数算出手段９２により算出された係数を用いてスペクトルの強度を調整し、強度が調整されたスペクトルを合成する。
画像作成手段９４は、合成手段９３により合成されたスペクトルに基づいて画像を作成する。

尚、受信器７およびスペクトル作成手段９１を合わせたものが、スペクトル生成手段に相当する。
図１に戻って説明を続ける。

操作部１１は、オペレータにより操作され、種々の情報をコンピュータ８に入力する。表示部１２は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
次に、被検体１３を撮影するときに実行されるスキャンについて説明する。

図４は、第１の形態で実行されるスキャンの説明図である。
第１の形態では、ローカライザスキャンＡ０、ＣＳＩ（Chemical Shift Imaging）スキャンＡ１、およびＭＲＳ（Magnetic Resonance Spectroscopy）スキャンＡ２が実行される。

ローカライザスキャンＡ０では、関心部位ＶＯＩ（後述する図７参照）を設定するときに使用される画像を取得するためのシーケンスが実行される。

ＣＳＩスキャンＡ１では、水の信号を抑制するためのＲＦパルスを含むシーケンスが実行される。

ＭＲＳスキャンＡ２では、水の信号を抑制するためのＲＦパルスを含まないシーケンスが実行される。

第１の形態では、図４に示すスキャンを実行し、ＣＳＩ画像を作成する。以下に、図４に示すスキャンを実行しＣＳＩ画像を作成するときのＭＲ装置の動作フローについて説明する。

図５はＭＲ装置の動作フローを示す図である。
ステップＳＴ１では、ローカライザスキャンＡ０を実行する。図６に、ローカライザスキャンＡ０により得られた画像の一例を示す。オペレータは、ローカライザスキャンＡ０により得られた画像に基づいて、ＣＳＩスキャンＡ１を実行するときの関心部位ＶＯＩを設定する（図７参照）。

図７は、オペレータが設定した関心部位ＶＯＩの一例を示す図である。図７では、説明の便宜上、４つの領域Ｖ１〜Ｖ４を有する関心部位ＶＯＩが設定されたとする。関心部位ＶＯＩを設定した後、ステップＳＴ２に進む。

ステップＳＴ２では、ＣＳＩスキャンＡ１およびＭＲＳスキャンＡ２が実行される。以下、ＣＳＩスキャンＡ１およびＭＲＳスキャンＡ２について順に説明する。

（１）ＣＳＩスキャンＡ１について
図８は、ＣＳＩスキャンＡ１で実行されるパルスシーケンスの一例を示す図である。
ＣＳＩスキャンＡ１は、関心部位ＶＯＩからＭＲ信号を発生させるためのシーケンスがｎ回実行される。第１の形態では、関心部位ＶＯＩは、４つの領域Ｖ１〜Ｖ４を有しているので、ｎ＝４、すなわち、シーケンスが４回実行される。図８では、実行される４回のシーケンスを、符号ａ１、ａ２、ａ３、およびａ４で示してある。また、図８には、シーケンスａ１の具体例が示されている。シーケンスａ１は、水信号を抑制するためのCHESS（Chemical Shift Selective）パルスを有する水抑制部ｐ１と、関心部位ＶＯＩからＭＲ信号Ｓ１を収集するためのデータ収集部ｐ２を有している。シーケンスａ１を実行することにより、ＭＲ信号Ｓ１が得られる。

データ収集部ｐ２は、ＲＦパルスとして、９０°パルスと２つの１８０°パルスとを有している。これらのＲＦパルスの中心周波数は、水の共鳴周波数に設定されている。また、データ収集部ｐ２は、ｘ軸に印加される位相エンコード勾配パルスＰｘと、ｙ軸に印加される位相エンコード勾配パルスＰｙとを有している。位相エンコード勾配パルスＰｘの磁場強度をｘとし、位相エンコード勾配パルスＰｙの磁場強度をｙとすると、シーケンスａ１では、（ｘ，ｙ）＝（ｘ１，ｙ１）に設定されている。

シーケンスａ１を実行した後、シーケンスａ２、ａ３、ａ４が順に実行される。シーケンスａ２、ａ３、およびａ４は、位相エンコード勾配パルスの磁場強度を除いて、シーケンスａ１と同じである。図９、図１０、図１１に、それぞれ、シーケンスａ２、ａ３、およびａ４を示す。

シーケンスａ２（図９参照）では、位相エンコード勾配パルスＰｘの磁場強度ｘは、ｘ＝ｘ１に設定されており、位相エンコード勾配パルスＰｙの磁場強度ｙは、ｙ＝ｙ２（＜ｙ１）に設定されている。したがって、（ｘ，ｙ）＝（ｘ１，ｙ２）である。シーケンスａ２は、シーケンスａ１と比較すると、位相エンコード勾配パルスＰｙの磁場強度が異なっているが、位相エンコード勾配パルスＰｘの磁場強度は同じである。シーケンスａ２を実行することにより、ＭＲ信号Ｓ２が得られる。

シーケンスａ３（図１０参照）では、位相エンコード勾配パルスＰｘの磁場強度ｘは、ｘ＝ｘ２（＜ｘ１）に設定されており、位相エンコード勾配パルスＰｙの磁場強度ｙは、ｙ＝ｙ１に設定されている。したがって、（ｘ，ｙ）＝（ｘ２，ｙ１）である。シーケンスａ３は、シーケンスａ１と比較すると、位相エンコード勾配パルスＰｘの磁場強度が異なっているが、位相エンコード勾配パルスＰｙの磁場強度は同じである。シーケンスａ３を実行することにより、ＭＲ信号Ｓ３が得られる。

シーケンスａ４（図１１参照）では、位相エンコード勾配パルスＰｘの磁場強度ｘは、ｘ＝ｘ２（＜ｘ１）に設定されており、位相エンコード勾配パルスＰｙの磁場強度ｙは、ｙ＝ｙ２（＜ｙ１）に設定されている。したがって、（ｘ，ｙ）＝（ｘ２，ｙ２）である。シーケンスａ４は、シーケンスａ１と比較すると、位相エンコード勾配パルスＰｘの磁場強度が異なっており、また、位相エンコード勾配パルスＰｙの磁場強度も異なっている。シーケンスａ４を実行することにより、ＭＲ信号Ｓ４が得られる。

したがって、ＣＳＩスキャンＡ１のシーケンスａ１〜ａ４を実行することにより、ＭＲ信号Ｓ１〜Ｓ４が得られる。図１２に、シーケンスａ１〜ａ４を実行することにより得られたＭＲ信号を概略的に示す。これらのＭＲ信号Ｓ１〜Ｓ４は受信コイル４で受信され、受信器７で信号処理される（図１３参照）。

図１３は、ＭＲ信号Ｓ１〜Ｓ４が受信コイル４で受信され、受信器７で信号処理されるときの様子を概略的に示す図である。
ＣＳＩスキャンＡ１のシーケンスａ１により得られたＭＲ信号Ｓ１は、受信コイル４で受信され、受信器７で信号処理（例えば、検波）される。受信コイル４は、２つのコイルエレメント４ａ及び４ｂを有しているので、コイルエレメントごとにＭＲ信号Ｓ１が受信される。各コイルエレメント４ａ，４ｂはＭＲ信号Ｓ１を電気信号Ｓａ，Ｓｂに変換し、受信器７に供給する。受信器７は、各コイルエレメントから受け取った電気信号Ｓａ，Ｓｂに対して、検波などの信号処理を行う。したがって、一方の受信器７は、コイルエレメント４ａによって得られた受信データＳ１ａを出力し、他方の受信器７は、コイルエレメント４ｂにより得られた受信データＳ１ｂを出力する。受信データＳ１ａは、コイルエレメント４ａで受信されたＭＲ信号Ｓ１の振幅と位相の情報を含んでおり、一方、受信データＳ１ｂは、コイルエレメント４ｂで受信されたＭＲ信号Ｓ１の振幅と位相の情報を含んでいる。

図１３では、ＭＲ信号Ｓ１から得られた受信データＳ１ａ，Ｓ１ｂを示したが、他のシーケンスａ２、ａ３、およびａ４により得られたＭＲ信号Ｓ２、Ｓ３、およびＳ４についても、ＭＲ信号Ｓ１と同じ処理が実行される。したがって、ＭＲ信号Ｓ２、Ｓ３、およびＳ４についても、２つの受信データが得られる。

図１４に、ＭＲ信号Ｓ１、Ｓ２、Ｓ３、およびＳ４から得られた受信データを概略的に示す。受信データＳ１ａ、Ｓ２ａ、Ｓ３ａ、およびＳ４ａは、コイルエレメント４ａにより得られた受信データであり、一方、受信データＳ１ｂ、Ｓ２ｂ、Ｓ３ｂ、およびＳ４ｂは、コイルエレメント４ｂにより得られた受信データである。これらの受信データは、プロセッサ９（図１参照）に伝送される。
次に、ＭＲＳスキャンＡ２について説明する。

（２）ＭＲＳスキャンＡ２について
図１５は、ＭＲＳスキャンＡ２で実行されるシーケンスの一例を示す図である。
ＭＲＳスキャンＡ２は、関心部位ＶＯＩを一つの領域（シングルボクセル）と見なして、関心部位ＶＯＩからＭＲ信号を発生させるためのスキャンである。ＭＲＳスキャンＡ２では、関心部位ＶＯＩからＭＲ信号を発生させるためのシーケンスｂが１回実行される。シーケンスｂは、CHESSパルスを備えていない点を除いて、図８に示すシーケンスａ１と同じである。

ＭＲＳスキャンＡ２を実行することにより、ＭＲ信号Ｓ５が得られる。ＭＲ信号Ｓ５は、受信コイル４で受信され、受信器で信号処理される（図１６参照）。

図１６は、ＭＲ信号Ｓ５が受信コイル４で受信され、受信器で信号処理されるときの様子を概略的に示す図である。
ＭＲＳスキャンＡ２のシーケンスｂにより得られたＭＲ信号Ｓ５は、受信コイル４で受信され、受信器７で信号処理（例えば、検波）される。受信コイル４は、２つのコイルエレメント４ａ及び４ｂを有しているので、コイルエレメントごとにＭＲ信号Ｓ５が受信される。各コイルエレメント４ａ，４ｂはＭＲ信号Ｓ５を電気信号Ｓａ，Ｓｂに変換し、受信器７に供給する。受信器７は、各コイルエレメントから受け取った電気信号Ｓａ，Ｓｂに対して、検波などの信号処理を行う。したがって、一方の受信器７は、コイルエレメント４ａによって得られた受信データＳ５ａを出力し、他方の受信器７は、コイルエレメント４ｂにより得られた受信データＳ５ｂを出力する。受信データＳ５ａおよびＳ５ｂは、プロセッサ９（図１参照）に伝送される。

したがって、ＣＳＩスキャンＡ１を実行することにより、受信データＳ１ａ〜Ｓ４ａおよびＳ１ｂ〜Ｓ４ｂ（図１４参照）を得ることができ、ＭＲＳスキャンＡ２を実行することにより、受信データＳ５ａおよびＳ５ｂ（図１６参照）を得ることができる。ＣＳＩスキャンＡ１およびＭＲＳスキャンＡ２を実行した後、ステップＳＴ３に進む。

ステップＳＴ３では、スペクトル作成手段９１（図３参照）が、ＣＳＩスキャンＡ１により得られた受信データ（図１４参照）に基づいて、関心部位ＶＯＩの各領域におけるスペクトルを作成する。ステップＳＴ３は、ステップＳＴ３１およびＳＴ３２を有している。以下、ステップＳＴ３１およびＳＴ３２について順に説明する。

ステップＳＴ３１では、変換手段９１ａ（図３参照）が、受信データをスペクトルに変換する（図１７参照）。

図１７は、受信データをスペクトルに変換するときの説明図である。
変換手段９１ａは、フーリエ変換を用いて、コイルエレメントごとに得られた受信データをスペクトルに変換する。

図１７（ａ）は、コイルエレメント４ａによる受信データＳ１ａ、Ｓ２ａ、Ｓ３ａ、およびＳ４ａから得られたスペクトルＭ１ａ、Ｍ２ａ、Ｍ３ａ、およびＭ４ａを示す。一方、図１７（ｂ）は、コイルエレメント４ｂによる受信データＳ１ｂ、Ｓ２ｂ、Ｓ３ｂ、およびＳ４ｂから得られたスペクトルＭ１ｂ、Ｍ２ｂ、Ｍ３ｂ、およびＭ４ｂを示す。ＣＳＩスキャンＡ１では、水抑制を行うためのCHESSパルスが使用されているので（図８参照）、スペクトルＭ１ａ〜Ｍ４ａおよびＭ１ｂ〜Ｍ４ｂは、水のピークが抑制されている。また、スペクトルＭ１ａ〜Ｍ４ａおよびＭ１ｂ〜Ｍ４ｂは、水のピークが抑制されているので、水よりも濃度の小さい代謝物のピークを得ることができる。したがって、ＣＳＩスキャンＡ１を実行することにより、代謝物のピークを含むスペクトル（以下、「代謝物スペクトル」と呼ぶ）Ｍ１ａ〜Ｍ４ａおよびＭ１ｂ〜Ｍ４ｂが得られる。
代謝物スペクトルを得た後、ステップＳＴ３２に進む。

ステップＳＴ３２では、補正手段９１ｂ（図３参照）が、図１７に示す代謝物スペクトルＭ１ａ〜Ｍ４ａおよびＭ１ｂ〜Ｍ４ｂの位相を補正する。図１８に、位相補正後の代謝物スペクトルを示す。

図１８（ａ１）は位相補正前の代謝物スペクトルＭ１ａ〜Ｍ４ａを示し、図１８（ａ２）は、代謝物スペクトルＭ１ａ〜Ｍ４ａを位相補正することにより得られた代謝物スペクトルを示す。図１８（ａ２）では、位相補正後の代謝物スペクトルには、符号「Ｐ−Ｍ１ａ」、「Ｐ−Ｍ２ａ」、「Ｐ−Ｍ３ａ」、および「Ｐ−Ｍ４ａ」が付されている。

一方、図１８（ｂ１）は位相補正前の代謝物スペクトルＭ１ｂ〜Ｍ４ｂを示し、図１８（ｂ２）は、代謝物スペクトルＭ１ｂ〜Ｍ４ｂを位相補正することにより得られた代謝物スペクトルを示す。図１８（ｂ２）では、位相補正後の代謝物スペクトルには、符号「Ｐ−Ｍ１ｂ」、「Ｐ−Ｍ２ｂ」、「Ｐ−Ｍ３ｂ」、および「Ｐ−Ｍ４ｂ」が付されている。

位相補正では、例えば、０次の位相補正および１次の位相補正を行うことができる。
位相補正を行った後、ステップＳＴ４に進む。

ステップＳＴ４では、変換手段９１ａが、ＭＲＳスキャンＡ２により得られた受信データＳ５ａおよびＳ５ｂ（図１６参照）をスペクトルに変換する。変換手段９１ａは、フーリエ変換を用いて、受信データＳ５ａおよびＳ５ｂをスペクトルに変換する。図１９に、スペクトルを概略的に示す。図１９（ａ）は、コイルエレメント４ａによる受信データＳ５ａから得られたスペクトルＷａを示す。一方、図１９（ｂ）は、コイルエレメント４ｂによる受信データＳ５ｂから得られたスペクトルＷｂを示す。ＭＲＳスキャンＡ２では、水抑制を行うためのCHESSパルスが使用されていないので（図１５参照）、スペクトルＷａおよびＷｂは、大きい水のピークを有している。したがって、ＭＲＳスキャンＡ２を実行することにより、水のピークを含むスペクトル（以下、「水スペクトル」と呼ぶ）ＷａおよびＷｂを作成することができる。
水スペクトルＷａおよびＷｂを作成した後、ステップＳＴ５に進む。

ステップＳＴ５では、係数算出手段９２（図３参照）が、水スペクトルＷａおよびＷｂに基づいて、後述するステップＳＴ６においてスペクトルを合成するときに使用する係数を求める。以下、係数の求め方について説明する。尚、ステップＳＴ５は、ステップＳＴ５１〜ＳＴ５３を有しているので、各ステップについて順に説明する。

ステップＳＴ５１では、係数算出手段９２が、水スペクトルＷａおよびＷｂに現れる水のピークの特性値を求める。水のピークの特性値とは、水のピークの大きさを表す値であり、例えば、水のピークのピーク値や、水のピークの面積である。本形態では、水のピークの特性値として、水のピークのピーク値を求める。図２０に、水スペクトルＷａに現れる水のピークのピーク値Ｐａと、水スペクトルＷｂに現れる水のピークのピーク値Ｐｂとを概略的に示す。ピーク値ＰａおよびＰｂを求めた後、ステップＳＴ５２に進む。

ステップＳＴ５２では、係数算出手段９２が、水スペクトルＷａおよびＷｂのノイズを求める。図２１に、求められたノイズＮａおよびＮｂを示す。ノイズＮａおよびＮｂは、例えば、水のピークから離れた周波数範囲の信号値を平均することにより求めることができる。ノイズＮａおよびＮｂを求めた後、ステップＳＴ５３に進む。

ステップＳＴ５３では、係数算出手段９２が、水のピーク値ＰａおよびＰｂをノイズＮａおよびＮｂで正規化する。図２２は、水のピーク値をノイズで正規化するときの説明図である。係数算出手段９２は、先ず、水のピーク値ＰａをノイズＮａで除算する。これにより、ノイズで正規化されたピーク値Ｐａ／Ｎａを求めることができる。また、係数算出手段９２は、水のピーク値ＰｂをノイズＮｂで除算する。これにより、水のピーク値ＰｂがノイズＮｂで正規化されたピーク値Ｐｂ／Ｎｂを求めることができる。正規化されたピーク値Ｐａ／ＮａおよびＰｂ／Ｎｂが、ステップＳＴ６においてスペクトルを合成するときの係数として用いられる。係数Ｐａ／ＮａおよびＰｂ／Ｎｂを求めた後、ステップＳＴ６に進む。

ステップＳＴ６では、合成手段９３（図３参照）が、係数Ｐａ／ＮａおよびＰｂ／Ｎｂに基づいて、位相補正後の代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａ（図１８（ａ２）参照）と、位相補正後の代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂ（図１８（ｂ２）参照）とを合成する。以下に、合成方法について説明する。

先ず、合成手段９３は、位相補正後の代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａに、係数Ｐａ／Ｎａを乗算する（図２３参照）。

図２３は、位相補正後の代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａに、係数Ｐａ／Ｎａを乗算するときの説明図である。
合成手段９３は、位相補正後の代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａに、係数Ｐａ／Ｎａを乗算する。これにより、Ｐａ／Ｎａの値に応じて、代謝物スペクトルの強度を調整することができる。図２３では、代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａに係数Ｐａ／Ｎａを乗算することにより得られた代謝物スペクトルを、符号ＭＡ１〜ＭＡ４で示してある。

次に、合成手段９３は、位相補正後の代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算する（図２４参照）。

図２４は、位相補正後の代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算するときの説明図である。
合成手段９３は、位相補正後の代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算する。これにより、Ｐｂ／Ｎｂの値に応じて、代謝物スペクトルの強度を調整することができる。図２４では、代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂに係数Ｐｂ／Ｎｂを乗算することにより得られた代謝物スペクトルを、符号ＭＢ１〜ＭＢ４で示してある。

次に、合成手段９３は、代謝物スペクトルＭＡ１〜ＭＡ４（図２３参照）と代謝物スペクトルＭＢ１〜ＭＢ４（図２４参照）とを加算する。図２５は、スペクトルを加算するときの説明図である。合成手段９３は、領域Ｖ１〜Ｖ４ごとにスペクトルを加算する。これにより、領域Ｖ１〜Ｖ４ごとに合成スペクトルＭ１〜Ｍ４を得ることができる。

第１の形態では、ステップＳＴ５において、水スペクトルＷａおよびＷｂから水のピーク値ＰａおよびＰｂを検出し、水のピーク値をノイズで正規化することにより、代謝物スペクトルの強度を調整するための係数Ｐａ／ＮａおよびＰｂ／Ｎｂを求めている。水は代謝物よりも濃度がかなり大きいので、水のピーク値は代謝物のピーク値よりもかなり大きくなる。したがって、係数Ｐａ／Ｎａを用いて代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａの強度を調整するとともに、係数Ｐｂ／Ｎｂを用いて代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂの強度を調整し、強度が調整された代謝物スペクトルを加算することにより、ＳＮ比の大きい合成スペクトルＭ１〜Ｍ４を得ることができる。
スペクトルを合成した後、ステップＳＴ７に進む。

ステップＳＴ７では、画像作成手段９４（図３参照）は、合成スペクトルに現れているピークを検出し、ピーク値を求める。そして、画像作成手段９４は、ピーク値に基づいて、代謝物ごとにＣＳＩ画像を作成する。ＣＳＩ画像を作成したら、フローを終了する。

第１の形態では、水のピーク値を用いて求めた係数Ｐａ／ＮａおよびＰｂ／Ｎｂで代謝物スペクトルの強度を調整し、強度が調整された後のスペクトルを加算するので、ＳＮ比の大きい合成スペクトルを得ることができる。

尚、生体内に高濃度で含まれている物質として、水の他に、脂質がある。したがって、水スペクトルに、脂質のピークが含まれている場合は、水のピーク値の代わりに、脂質のピーク値を検出し、ノイズで正規化された脂質のピーク値を用いて代謝物スペクトルを合成することも考えられる。しかし、代謝物は、一般的に、細胞膜近傍の束縛水からの信号であるので、代謝物の濃度分布は、脂質よりも水の濃度分布に近いと推定することができる。したがって、脂質のピーク値よりも水のピーク値を用いた方が、代謝物の濃度分布を反映した合成スペクトルを得ることができるので、水のピーク値を検出することが望ましい。

（２）第２の形態
第２の形態では、第１の形態と同様に、図４に示すスキャンを実行し、ＣＳＩ画像を作成する。尚、第２の形態のＭＲ装置の構成は、第１の形態のＭＲ装置と比較すると、プロセッサ９が実行する処理が異なるが、プロセッサ９以外の部分は、第１の形態のＭＲ装置と同じ構成である。したがって、第２のＭＲ装置の説明にあたっては、プロセッサ９について主に説明する。

図２６は、第２の形態においてプロセッサ９が実行する処理を示す図である。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出すことにより、スペクトル作成手段９１１〜画像作成手段９１５などを構成する。

スペクトル作成手段９１１は、受信器７から得られたデータに基づいて、スペクトルを作成する。スペクトル作成手段９１１は、位相算出手段９１１ａ、補正手段９１１ｂ、および変換手段９１１ｃを有している。

位相算出手段９１１ａは、ＭＲＳスキャンＡ２により得られた受信データＳ５ａおよびＳ５ｂに基づいて、ＭＲ信号の位相を算出する（後述する図２８参照）。補正手段９１１ｂは、位相算出手段９１１ａにより算出された位相に基づいて、ＣＳＩスキャンＡ１により得られた受信データの位相を補正する（後述する図２９および図３０参照）。変換手段９１１ｃは、補正手段９１１ｂにより得られた位相補正後のデータを、スペクトルに変換する。
係数算出手段９１２は、スペクトルの強度を調整するための係数を算出する。
合成手段９１３は、係数算出手段９１２により算出された係数を用いてスペクトルの強度を調整し、強度が調整されたスペクトルを合成する。
補正手段９１４は、合成手段９１３により得られた合成スペクトルの位相を補正する。
画像作成手段９１５は、補正手段９１４により得られた位相補正後のスペクトルに基づいて画像を作成する。
第２の形態のプロセッサ９は、上記のように構成されている。

次に、第２の形態において、図４に示すスキャンを実行しＣＳＩ画像を作成するときのＭＲ装置の動作フローについて説明する。

図２７は、第２の形態におけるＭＲ装置の動作フローを示す図である。
ステップＳＴ１およびＳＴ２は、第１の形態と同じであるので、説明は省略する。ステップＳＴ２においてＣＳＩスキャンＡ１およびＭＲＳスキャンＡ２を実行した後、ステップＳＴ２１に進む。

ステップＳＴ２１では、位相算出手段９１１ａ（図２６参照）が、ＭＲＳスキャンＡ２により得られた受信データＳ５ａおよびＳ５ｂ（図１６参照）に基づいて、ＭＲ信号Ｓ５の位相を算出する。図２８は、位相を算出するときの説明図である。位相算出手段９１１ａは、受信データＳ５ａからＭＲ信号Ｓ５の位相情報を検出し、各時刻における位相を算出する。したがって、受信データＳ５ａから、ＭＲ信号の位相の時間変化を表す位相データＳＰ５ａを得ることができる。また、位相算出手段９１１ａは、受信データＳ５ｂからＭＲ信号Ｓ５の位相情報を検出し、各時刻における位相を算出する。したがって、受信データＳ５ｂから、ＭＲ信号の位相の時間変化を表す位相データＳＰ５ｂを得ることができる。位相データＳＰ５ａおよびＳＰ５ｂを求めた後、ステップＳＴ２２に進む。

ステップＳＴ２２では、補正手段９１１ｂ（図２６参照）が、ステップＳＴ２１で求めた位相データＳＰ５ａおよびＳＰ５ｂに基づいて、ＣＳＩスキャンにより得られた受信データＳ１ａ〜Ｓ４ａおよびＳ１ｂ〜Ｓ４ｂ（図１４参照）の位相を補正する。位相データＳＰ５ａは、コイルエレメント４ａにより得られた受信データＳ１ａ〜Ｓ４ａの位相を補正するために使用される。一方、位相データＳＰ５ｂは、コイルエレメント４ｂにより得られた受信データＳ１ｂ〜Ｓ４ｂの位相を補正するために使用される。以下、補正方法について図２９および図３０を参照しながら説明する。

図２９は、位相データＳＰ５ａに基づいて、コイルエレメント４ａにより得られた受信データＳ１ａ〜Ｓ４ａの位相を補正するときの説明図である。
補正手段９１１ｂは、先ず、受信データＳ１ａの各時刻における位相を求める。次に、受信データＳ１ａの各時刻における位相を、位相データＳＰ５ａの各時刻における位相で補正する。具体的には、受信データＳ１ａの各時刻における位相から、位相データＳ５ａの各時刻における位相を減算する。したがって、位相データＳ５ａに基づいて受信データＳ１ａの位相を補正することができる。図２９では、位相補正された受信データＳ１ａを、符号「Ｓ１ａ＿ｐ」で示してある。

上記のように、受信データＳ１ａの各時刻における位相から、位相データＳ５ａの各時刻における位相を減算することにより、受信データＳ１ａの位相から、渦電流により生じる位相成分と、ＲＦパルスの中心周波数とは異なる周波数で共鳴すること（オフレゾナンス）により生じる位相成分とを除去することができる。以下に、受信データＳ１ａの位相から、渦電流により生じる位相成分と、オフレゾナンスにより生じる位相成分とを除去することができる理由について説明する。

先ず、受信データＳ１ａの位相について考える。受信データＳ１ａは、ＣＳＩスキャンＡ１のシーケンスａ１（図８参照）により得られたＭＲ信号Ｓ１をコイルエレメント４ａで受信したときの位相を表している。シーケンスａ１は水抑制を行うので、シーケンスａ１を実行すると、水のスピンによるＭＲ信号は抑制される。したがって、シーケンスａ１を実行することにより、濃度の低い代謝物のスピンによるＭＲ信号の位相情報を取得することができる。このため、受信データＳ１ａの位相をφｓ１とすると、位相φｓ１は、代謝物の化学シフトによる位相成分ωｍを用いて、以下の式で表すことができる。
φｓ１＝ωｏ＋ωｍ＋ωｐ＋ωｅ ・・・（１）
ここで、ωｏ：オフレゾナンスにより生じる位相成分
ωｍ：代謝物の化学シフトによる位相成分
ωｐ：位相エンコード勾配パルスにより生じる位相成分
ωｅ：渦電流により生じる位相成分

次に、位相データＳ５ａの位相について考える。位相データＳ５ａは、ＭＲＳスキャンＡ２のシーケンスｂ（図１５参照）により得られたＭＲ信号Ｓ５をコイルエレメント４ａで受信したときの位相を表している。ＭＲＳスキャンＡ２のシーケンスｂは水抑制が行われていないので、シーケンスｂを実行すると、代謝物のスピンによるＭＲ信号よりも、水のスピンによるＭＲ信号が大きくなる。したがって、ＭＲ信号Ｓ５の位相を考える場合、代謝物の化学シフトによる位相成分ωｍは無視することができ、水の化学シフトによる位相成分ωｗを考えればよい。ここで、位相データＳ５ａの位相をφｓ５とすると、位相φｓ５は、水の化学シフトによる位相成分ωｗを用いて、以下の式で表すことができる。
φｓ５＝ωｏ＋ωｗ＋ωｐ＋ωｅ ・・・（２）
ここで、ωｏ：オフレゾナンスにより生じる位相成分
ωｗ：水の化学シフトによる位相成分
ωｐ：位相エンコード勾配パルスにより生じる位相成分
ωｅ：渦電流により生じる位相成分

ただし、位相データＳ５ａは、ＭＲＳスキャンＡ２のシーケンスｂ（図１５参照）により得られたＭＲ信号Ｓ５ａの位相を表している。ここで、ＭＲＳスキャンＡ２のシーケンスｂで使用されているＲＦパルスについて考えると、シーケンスｂのＲＦパルス（９０°パルス、１８０°パルス）の中心周波数は、第１の形態で述べたように、水の共鳴周波数に設定されている。したがって、式（２）の右辺第２項ωｗはゼロと考えることができる。また、ＭＲＳスキャンＡ２では、シーケンスｂが１回のみ実行されているので、位相エンコード勾配パルスにより生じる位相変化は発生しないと考えることができる。したがって、式（２）の右辺第３項ωｐもゼロと考えることができる。このため、式（２）は、以下の式で表すことができる。
φｓ５＝ωｏ＋ωｅ ・・・（３）

次に、φｓ１からφｓ５を減算することにより、φｓ１の位相を補正する。ここで、位相補正後のφｓ１を「φｓ１＿ｐ」で表すと、φｓ１＿ｐは、式（１）および式（３）を用いて、以下の式で表すことができる。

φｓ１＿ｐ＝φｓ１−φｓ５
＝（ωｏ＋ωｃ＋ωｐ＋ωｅ）−（ωｏ＋ωｅ）
＝ωｃ＋ωｐ ・・・（４）

式（４）を参照すると、位相補正後のφｓ１＿ｐには、オフレゾナンスにより生じる位相成分ωｏと、渦電流により生じる位相成分ωｅが含まれていないことがわかる。したがって、受信データＳ１ａの位相φｓ１から位相データＳ５ａの位相φｓ５を減算することにより、受信データＳ１ａの位相φｓ１から、オフレゾナンスにより生じる位相成分ωｏと、渦電流により生じる位相成分ωｅを除去できることがわかる。

尚、上記の説明では、受信データＳ１ａの位相を補正する方法が示されているが、他の受信データＳ２ａ、Ｓ３ａ、Ｓ４ａの位相も、位相データＳＰ５ａを用いて補正される。図２９では、位相補正された受信データＳ２ａ、Ｓ３ａ、およびＳ４ａを、それぞれ、符号「Ｓ２ａ＿ｐ」、「Ｓ３ａ＿ｐ」、および「Ｓ４ａ＿ｐ」で示してある。

したがって、受信データＳ１ａ〜Ｓ４ａの位相から、オフレゾナンスにより生じる位相成分ωｏと、渦電流により生じる位相成分ωｅを除去することができる。

次に、補正手段９１１ｂは、位相データＳＰ５ｂに基づいて、コイルエレメント４ｂにより得られた受信データＳ１ｂ〜Ｓ４ｂに含まれる位相を補正する（図３０参照）。

図３０は、位相データＳＰ５ｂに基づいて、コイルエレメント４ｂにより得られた受信データＳ１ｂ〜Ｓ４ｂの位相を補正するときの説明図である。
補正手段９１１ｂは、先ず、受信データＳ１ｂの各時刻における位相を求める。次に、受信データＳ１ｂの各時刻における位相を、位相データＳＰ５ｂの各時刻における位相で補正する。具体的には、受信データＳ１ｂの各時刻における位相から、位相データＳ５ｂの各時刻における位相を減算する。したがって、位相データＳ５ｂに基づいて受信データＳ１ｂの位相を補正することができる。図３０では、位相補正された受信データＳ１ｂを、符号「Ｓ１ｂ＿ｐ」で示してある。

以下同様に、他の受信データＳ２ｂ、Ｓ３ｂ、Ｓ４ｂの位相も、位相データＳＰ５ｂを用いて補正される。図３０では、位相補正された受信データＳ２ｂ、Ｓ３ｂ、およびＳ４ｂを、それぞれ、符号「Ｓ２ｂ＿ｐ」、「Ｓ３ｂ＿ｐ」、および「Ｓ４ｂ＿ｐ」で示してある。

したがって、受信データＳ１ｂ〜Ｓ４ｂの位相から、オフレゾナンスにより生じる位相成分ωｏと、渦電流により生じる位相成分ωｅを除去することができる。
受信データの位相を補正した後、ステップＳＴ２３に進む。

ステップＳＴ２３では、変換手段９１１ｃ（図２６参照）が、ステップＳＴ２２により得られた位相補正後の受信データ（図２９および図３０参照）を、関心部位ＶＯＩの各領域におけるスペクトルに変換する。図３１に、得られたスペクトルを概略的に示す。図３１（ａ）は、位相補正後の受信データＳ１ａ＿ｐ、Ｓ２ａ＿ｐ、Ｓ３ａ＿ｐ、およびＳ４ａ＿ｐにより得られた代謝物スペクトルＭ１ａ、Ｍ２ａ、Ｍ３ａ、およびＭ４ａを示す。一方、図３１（ｂ）は、位相補正後の受信データＳ１ｂ＿ｐ、Ｓ２ｂ＿ｐ、Ｓ３ｂ＿ｐ、およびＳ４ｂ＿ｐにより得られた−０代謝物スペクトルＭ１ｂ、Ｍ２ｂ、Ｍ３ｂ、およびＭ４ｂを示す。
スペクトルを求めた後、ステップＳＴ４進む。

ステップＳＴ４およびＳＴ５は、第１の形態と同じであるので説明は省略する。ステップＳＴ４において水スペクトルＷａおよびＷｂを作成し（図１９参照）、ステップＳＴ５において係数算出手段９１２（図２６参照）が係数Ｐａ／ＮａおよびＰｂ／Ｎｂ（図２２参照）を求めた後、ステップＳＴ６に進む。

ステップＳＴ６では、合成手段９１３（図２６参照）が、係数Ｐａ／ＮａおよびＰｂ／Ｎｂに基づいて、代謝物スペクトルＭ１ａ〜Ｍ４ａ（図３１（ａ）参照）と、代謝物スペクトルＭ１ｂ〜Ｍ４ｂ（図３１（ｂ）参照）とを合成する。以下に、合成方法について説明する。

先ず、合成手段９１３は、位相補正後の代謝物スペクトルＭ１ａ〜Ｍ４ａに、係数Ｐａ／Ｎａを乗算する（図３２参照）。

図３２は、位相補正後の代謝物スペクトルＭ１ａ〜Ｍ４ａに、係数Ｐａ／Ｎａを乗算するときの説明図である。
合成手段９１３は、位相補正後の代謝物スペクトルＭ１ａ〜Ｍ４ａに、係数Ｐａ／Ｎａを乗算する。これにより、Ｐａ／Ｎａの値に応じて、代謝物スペクトルの強度を調整することができる。図３２では、代謝物スペクトルＭ１ａ〜Ｍ４ａに係数Ｐａ／Ｎａを乗算することにより得られた代謝物スペクトルを、符号ＭＡ１〜ＭＡ４で示してある。

次に、合成手段９１３は、位相補正後の代謝物スペクトルＭ１ｂ〜Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算する（図３３参照）。

図３３は、位相補正後の代謝物スペクトルＭ１ｂ〜Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算するときの説明図である。
合成手段９１３は、位相補正後の代謝物スペクトルＭ１ｂ〜Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算する。これにより、Ｐｂ／Ｎｂの値に応じて、代謝物スペクトルの強度を調整することができる。図３３では、代謝物スペクトルＭ１ｂ〜Ｍ４ｂに係数Ｐｂ／Ｎｂを乗算することにより得られた代謝物スペクトルを、符号ＭＢ１〜ＭＢ４で示してある。

次に、合成手段９１３は、代謝物スペクトルＭＡ１〜ＭＡ４（図３２参照）と代謝物スペクトルＭＢ１〜ＭＢ４（図３３参照）とを加算する。図３４は、スペクトルを加算するときの説明図である。合成手段９１３は、領域Ｖ１〜Ｖ４ごとにスペクトルを加算する。これにより、領域Ｖ１〜Ｖ４ごとに合成スペクトルＭ１〜Ｍ４を得ることができる。

尚、合成スペクトルＭ１〜Ｍ４を比較すると、合成スペクトルＭ１はピークが上向きであるが、合成スペクトルＭ３およびＭ４はピークが下向きであり、更に、合成スペクトルＭ２のピークはスプリットしている。したがって、合成スペクトルのピークを上向きに揃える必要がある。そこで、第２の形態では、合成スペクトルＭ１〜Ｍ４のピークを上向きの揃えるための位相補正を行う。この位相補正を行うために、ステップＳＴ６１に進む。

ステップＳＴ６１では、補正手段９１４（図２６参照）が、合成スペクトルＭ１〜Ｍ４のピークが上向きに揃うように、合成スペクトルＭ１〜Ｍ４の位相を補正する。合成スペクトルＭ１〜Ｍ４のピークを上向きに揃えた後、ステップＳＴ７に進み、画像作成手段９１５がＣＳＩ画像を作成し、フローを終了する。

第２の形態でも、第１の形態と同様にＳＮ比の大きい合成スペクトルを得ることができる。また、第２の形態では、渦電流による位相成分ωｅおよびオフレゾナンスによる位相成分ωｏを除去した後で、各領域のスペクトルを作成するので、位相補正の精度を高めることができる。

（３）第３の形態
図３５は、第３の形態で実行されるスキャンの説明図である。
第３の形態では、ローカライザスキャンＡ０、ＣＳＩスキャンＡ１、ＭＲＳスキャンＡ２、およびＭＲＳスキャンＡ３が実行される。

スキャンＡ０、Ａ１、およびＡ２は第１の形態と同じであるので、説明は省略する。
ＭＲＳスキャンＡ３では、水の信号を抑制するためのＲＦパルスを含むシーケンスが実行される。

第３の形態では、図３５に示すスキャンを実行し、ＣＳＩ画像を作成する。以下に、図３５に示すスキャンを実行しＣＳＩ画像を作成するときのＭＲ装置の動作フローについて説明する。尚、第３の形態のＭＲ装置の構成は、第１の形態のＭＲ装置と比較すると、プロセッサ９が実行する処理が異なるが、プロセッサ９以外の部分は、第１の形態のＭＲ装置と同じ構成である。したがって、第３のＭＲ装置の説明にあたっては、プロセッサ９について主に説明する。

図３６は、第３の形態においてプロセッサ９が実行する処理を示す図である。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出すことにより、スペクトル作成手段９１〜画像作成手段９４などを構成する。

スペクトル作成手段９１は、受信器７から得られたデータに基づいて、スペクトルを作成する。スペクトル作成手段９１は、変換手段９１ａおよび補正手段９１ｂを有している。変換手段９１ａは、受信器７から得られたデータをスペクトルに変換する。補正手段９１ｂは、変換手段９１ａにより得られたスペクトルの位相を補正する。差分手段９１ｃは、スペクトルを差分することにより、差分スペクトルを求める。
係数算出手段９２は、差分手段９１ｃにより求められた差分スペクトルに基づいて、スペクトルの強度を調整するための係数を算出する。
合成手段９３は、係数算出手段９２により算出された係数を用いてスペクトルの強度を調整し、強度が調整されたスペクトルを合成する。
画像作成手段９４は、合成手段９３により合成されたスペクトルに基づいて画像を作成する。
第３の形態のプロセッサ９は、上記のように構成されている。

次に、第３の形態において、スキャンを実行しＣＳＩ画像を作成するときのＭＲ装置の動作フローについて説明する。

図３７は、第３の形態におけるＭＲ装置の動作フローを示す図である。
ステップＳＴ１は、第１の形態と同じであるので、説明は省略する。ステップＳＴ１においてローカライザスキャンＡ０実行した後、ステップＳＴ２に進む。

ステップＳＴ２では、ＣＳＩスキャンＡ１、ＭＲＳスキャンＡ２、およびＭＲＳスキャンＡ３が実行される。以下、ＣＳＩスキャンＡ１およびＭＲＳスキャンＡ２は、第１の形態と同じであるので説明は省略し、以下では、ＭＲＳスキャンＡ３について詳細に説明する。

ＣＳＩスキャンＡ１およびＭＲＳスキャンＡ２を実行した後、ＭＲＳスキャンＡ３が実行される（図３８参照）。

図３８は、ＭＲＳスキャンＡ３で実行されるパルスシーケンスの一例を示す図である。
ＭＲＳスキャンＡ３は、関心部位ＶＯＩを一つの領域（シングルボクセル）と見なして、関心部位ＶＯＩからＭＲ信号を発生させるためのスキャンである。ＭＲＳスキャンＡ３では、関心部位ＶＯＩからＭＲ信号を発生させるためのシーケンスｃが１回実行される。シーケンスｃは、水を抑制するためのCHESSパルスを備えている点を除いて、図１５に示すシーケンスｂと同じである。

ＭＲＳスキャンＡ３を実行することにより、ＭＲ信号Ｓ６が得られる。ＭＲ信号Ｓ６は、受信コイル４で受信され、受信器で信号処理される（図３９参照）。

図３９は、ＭＲ信号Ｓ６が受信コイル４で受信され、受信器で信号処理されるときの様子を概略的に示す図である。
受信コイル４は、２つのコイルエレメント４ａ及び４ｂを有しているので、コイルエレメントごとにＭＲ信号Ｓ６が受信される。各コイルエレメント４ａ，４ｂはＭＲ信号Ｓ６を電気信号Ｓａ，Ｓｂに変換し、受信器７に供給する。受信器７は、各コイルエレメントから受け取った電気信号Ｓａ，Ｓｂに対して、検波などの信号処理を行う。したがって、一方の受信器７は、コイルエレメント４ａによって得られた受信データＳ６ａを出力し、他方の受信器７は、コイルエレメント４ｂにより得られた受信データＳ６ｂを出力する。受信データＳ６ａおよびＳ６ｂは、プロセッサ９（図１参照）に伝送される。
ＣＳＩスキャンＡ１〜ＭＲＳスキャンＡ３を実行した後、ステップＳＴ３に進む。

ステップＳＴ３およびＳＴ４は、第１の形態と同じであるので説明は省略する。ステップＳＴ３において位相補正後の代謝物スペクトル（図１８（ａ２）および（ｂ２））を作成し、ステップＳＴ４において水スペクトルＷａおよびＷｂを作成した後（図１９参照）、ステップＳＴ４１に進む。

ステップＳＴ４１では、変換手段９１ａ（図３６参照）が、ＭＲＳスキャンＡ３により得られた受信データＳ６ａおよびＳ６ｂ（図３９参照）を、スペクトルに変換する。図４０に、スペクトルを概略的に示す。図４０（ａ）は、コイルエレメント４ａの受信データＳ６ａから得られたスペクトルＭａを示しており、図４０（ｂ）は、コイルエレメント４ｂの受信データＳ６ｂから得られたスペクトルＭｂを示している。ＭＲＳスキャンＡ３では、水抑制を行うためのCHESSパルスが使用されているので（図３８参照）、スペクトルＭａおよびＭｂは、代謝物のピークを含む代謝物スペクトルである。
代謝物スペクトルＭａおよびＭｂを作成した後、ステップＳＴ４２に進む。

ステップＳＴ４２では、差分手段９１ｃ（図３６参照）が、ステップＳＴ４で得られた水スペクトルＷａおよびＷｂ（図１９参照）と、ステップＳＴ４１で得られた代謝物スペクトルＭａおよびＭｂ（図４０参照）との差分スペクトルを求める。以下に、差分スペクトルの求め方について、図４１および図４２を参照しながら説明する。

差分手段９１ｃは、先ず、図４１に示すように、コイルエレメント４ａにより得られた水スペクトルＷａから、コイルエレメント４ａにより得られた代謝物スペクトルＭａを減算する。水スペクトルＷａから代謝物スペクトルＭａを減算することにより、差分スペクトルＷａ＿difを得ることができる。この減算により、水スペクトルＷａに現れる水のピークを残したまま、水スペクトルＷａに現れる代謝物のピークやその他の不要なピーク（例えば、脂質のピーク）を小さくすることができる。

次に、差分手段９１ｃは、図４２に示すように、コイルエレメント４ｂにより得られた水スペクトルＷｂから、コイルエレメント４ｂにより得られた代謝物スペクトルＭｂを減算する。水スペクトルＷｂから代謝物スペクトルＭｂを減算することにより、差分スペクトルＷｂ＿difを得ることができる。この減算により、水スペクトルＷｂに現れる水のピークを残したまま、水スペクトルＷｂに現れる代謝物のピークやその他の不要なピーク（例えば、脂質のピーク）を小さくすることができる。
差分スペクトルＷａ＿difおよびＷｂ＿difを求めた後、ステップＳＴ５に進む。

ステップＳＴ５では、係数算出手段９２（図３６参照）が、差分スペクトルＷａ＿difおよびＷｂ＿difに基づいて、後述するステップＳＴ６においてスペクトルを合成するときに使用する係数を求める。以下、係数の求め方について説明する。尚、ステップＳＴ５は、ステップＳＴ５１〜ＳＴ５３を有しているので、各ステップについて順に説明する。

ステップＳＴ５１では、係数算出手段９２が、差分スペクトルＷａ＿difおよびＷｂ＿difに現れる水のピーク値を求める。図４３に、求められたピーク値ＰａおよびＰｂを示す。ピーク値ＰａおよびＰｂを求めた後、ステップＳＴ５２に進む。

ステップＳＴ５２では、係数算出手段９２が、水スペクトルＷａおよびＷｂ（図１９参照）を求める。図４４に、求められたノイズＮａおよびＮｂを示す。ノイズＮａおよびＮｂは、例えば、水のピークから離れた周波数範囲の信号値を平均することにより求めることができる。ノイズＮａおよびＮｂを求めた後、ステップＳＴ５３に進む。

ステップＳＴ５３では、係数算出手段９２が、差分スペクトルから求めた水のピーク値を、水スペクトルから求めたノイズで正規化する。図４５は、水のピーク値をノイズで正規化するときの説明図である。係数算出手段９２は、先ず、水のピーク値ＰａをノイズＮａで除算する。これにより、ノイズで正規化されたピーク値Ｐａ／Ｎａを求めることができる。また、係数算出手段９２は、水のピーク値ＰｂをノイズＮｂで除算する。これにより、水のピーク値ＰｂがノイズＮｂで正規化されたピーク値Ｐｂ／Ｎｂを求めることができる。正規化されたピーク値Ｐａ／ＮａおよびＰｂ／Ｎｂが、ステップＳＴ６においてスペクトルを合成するときの係数として用いられる。係数Ｐａ／ＮａおよびＰｂ／Ｎｂを求めた後、ステップＳＴ６に進む。

ステップＳＴ６では、合成手段９３（図３６参照）が、係数Ｐａ／ＮａおよびＰｂ／Ｎｂに基づいて、ステップＳＴ３で求めた位相補正後の代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａおよびＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂ（図１８（ａ２）および（ｂ２）参照）を合成する。以下に、合成方法について説明する。

先ず、合成手段９３は、位相補正後の代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａに、係数Ｐａ／Ｎａを乗算する（図４６参照）。

図４６は、位相補正後の代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａに、係数Ｐａ／Ｎａを乗算するときの説明図である。
合成手段９３は、位相補正後の代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａに、係数Ｐａ／Ｎａを乗算する。これにより、Ｐａ／Ｎａの値を考慮して、代謝物スペクトルの強度を調整することができる。図４６では、代謝物スペクトルＰ−Ｍ１ａ〜Ｐ−Ｍ４ａに係数Ｐａ／Ｎａを乗算することにより得られた代謝物スペクトルを符号ＭＡ１〜ＭＡ４で示してある。

次に、合成手段９３は、位相補正後の代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算する（図４７参照）。

図４７は、位相補正後の代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算するときの説明図である。
合成手段９３は、位相補正後の代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂに、係数Ｐｂ／Ｎｂを乗算する。これにより、Ｐｂ／Ｎｂの値を考慮して、代謝物スペクトルの強度を調整することができる。図４７では、代謝物スペクトルＰ−Ｍ１ｂ〜Ｐ−Ｍ４ｂに係数Ｐｂ／Ｎｂを乗算することにより得られた代謝物スペクトルを符号ＭＢ１〜ＭＢ４で示してある。

次に、合成手段９３は、代謝物スペクトルＭＡ１〜ＭＡ４（図４６参照）と代謝物スペクトルＭＢ１〜ＭＢ４（図４７参照）とを加算する。これにより、領域ごとに合成スペクトルＭ１〜Ｍ４（図２５参照）を得ることができる。スペクトルを合成した後、ステップＳＴ７に進む。

ステップＳＴ７では、画像作成手段９４（図３６参照）が、合成スペクトルＭ１〜Ｍ４に現れているピークを検出し、ピーク値を求める。そして、画像作成手段９４は、ピーク値に基づいて、代謝物ごとにＣＳＩ画像を作成する。ＣＳＩ画像を作成したら、フローを終了する。

第３の形態でも、第１の形態と同様に、ＳＮ比の大きい合成スペクトルを得ることができる。また、第３の形態では、ステップＳＴ４２において水スペクトルから代謝物スペクトルを減算しているので、水スペクトルに現れる水のピークを残したまま、水スペクトルに現れる代謝物のピークや脂質のピークを小さくすることができる。したがって、差分スペクトルは、水スペクトルよりも、代謝物のピークが小さいので、差分スペクトルから水のピークを検出することにより、代謝物のピークや脂質のピークを水のピークとして誤検出する可能性を更に低減することができる。このため、ステップＳＴ６において代謝物スペクトルの強度を調整するために使用される係数の信頼性を高めることができる。

尚、第３の形態では、水スペクトルＷａおよびＷｂのノイズを求めているが、水スペクトルの代わりに、代謝物スペクトルＭａおよびＭｂのノイズを求め、代謝物スペクトルＭａおよびＭｂから求めたノイズを用いて差分スペクトルの水のピーク値を正規化してもよい。

最後に、本発明により、ＳＮ比が高い合成スペクトルが得られることを検証するために、スペクトルを実測した。以下に、実測により得られたスペクトルを示す。

図４８は実測により得られたスペクトルを示す図である。
図４８（ａ）は、第１の形態の方法（水スペクトルから水のピーク値とノイズを求め、ピーク値をノイズで正規化することにより得られた係数を用いてスペクトルを合成する方法）により得られた合成スペクトルを示す図である。
図４８（ｂ）は、第３の形態の方法（差分スペクトルから水のピーク値を求めるとともに、水スペクトルからノイズを求め、ピーク値をノイズで正規化することにより得られた係数を用いてスペクトルを合成する方法）により得られた合成スペクトルを示す図である。
図４８（ｃ）は、本発明とは別の方法（代謝物スペクトルのピーク値をノイズで正規化することにより得られた係数を用いてスペクトルを合成する方法）により得られた合成スペクトルを示す図である。

図４８を参照すると、図４８（ａ）および図４８（ｂ）の合成スペクトルは、図４８（ｃ）の合成スペクトルよりも、ＳＮ比が大きくなっている。したがって、本発明の方法により、ＳＮ比の大きい合成スペクトルが得られることがわかる。

２ マグネット
３ テーブル
３ａ クレードル
４ 受信コイル
５ 送信器
６ 勾配磁場電源
７ 受信器
８ コンピュータ
９ プロセッサ
１０ メモリ
１１ 操作部
１２ 表示部
１３ 被検体
２１ ボア
１００ ＭＲ装置

Claims (14)

1. 水および代謝物を含む関心部位からＭＲ信号を発生させるためのシーケンスを実行するスキャン手段であって、前記関心部位に含まれる水の信号を抑制するためのＲＦパルスを含む第１のシーケンスと、前記ＲＦパルスを含まない第２のシーケンスとを実行するスキャン手段と、
   前記第１のシーケンスにより発生した第１のＭＲ信号と、前記第２のシーケンスにより発生した第２のＭＲ信号とを受信する複数のコイルエレメントを備えたコイルと、
   各コイルエレメントで受信された第１のＭＲ信号に基づいて、前記コイルエレメントごとに、前記代謝物のピークを含む第１のスペクトルを生成するとともに、各コイルエレメントで受信された第２のＭＲ信号に基づいて、前記コイルエレメントごとに、前記水のピークを含む第２のスペクトルを生成するスペクトル生成手段と、
   前記第２のスペクトルの水のピークの特性値に基づいて、前記コイルエレメントごとに得られた前記第１のスペクトルの強度を調整するための係数を算出する係数算出手段と、
   前記係数を用いて前記第１のスペクトルの各々の強度を調整し、強度が調整された前記第１のスペクトルを合成する合成手段と、
   を有し、
   前記係数算出手段は、
   前記第２のスペクトルのノイズを求め、前記ノイズで前記特性値を正規化し、正規化された前記特性値に基づいて、前記係数を算出する、磁気共鳴装置。
2. 前記スペクトル生成手段は、
   前記各コイルエレメントで受信された第１のＭＲ信号に基づいて、前記コイルエレメントごとに、前記第１のＭＲ信号の位相と振幅の情報を含む第１のデータを生成するデータ生成手段と、
   前記コイルエレメントごとに生成された第１のデータに基づいて、前記第１のスペクトルを作成する作成手段と、
   を有する、請求項１に記載の磁気共鳴装置。
3. 前記作成手段は、
   前記第１のデータを、代謝物のピークを含む第３のスペクトルに変換する変換手段と、
   前記第３のスペクトルの位相を補正することにより、前記第１のスペクトルを得る補正手段と、
   を有する請求項２に記載の磁気共鳴装置。
4. 前記第１のシーケンスは、第１の軸に印加される第１の位相エンコード勾配パルスと、第２の軸に印加される第２の位相エンコード勾配パルスとを有する、請求項３に記載の磁気共鳴装置。
5. 前記関心部位は複数の領域を有しており、
   前記スキャン手段は、前記第１の位相エンコード勾配パルスの磁場強度と前記第２の位相エンコード勾配パルスの磁場強度との組み合わせが異なる複数の第１のシーケンスを実行し、
   前記コイルは、前記複数の第１のシーケンスにより発生した複数の第１のＭＲ信号を受信し、
   前記データ生成手段は、前記各コイルエレメントで受信された複数の第１のＭＲ信号に基づいて、前記コイルエレメントごとに複数の第１のデータを生成し、
   前記変換手段は、前記コイルエレメントごとに生成された複数の第１のデータを、前記関心部位の各領域における前記第３のスペクトルに変換する、請求項４に記載の磁気共鳴装置。
6. 前記データ生成手段は、前記各コイルエレメントで受信された第２のＭＲ信号に基づいて、前記コイルエレメントごとに、前記第２のＭＲ信号の位相と振幅の情報を含む第２のデータを生成し、
   前記作成手段は、前記コイルエレメントごとに生成された前記第２のデータに基づいて、前記第２のスペクトルを作成する、請求項３〜５のうちのいずれか一項に記載の磁気共鳴装置。
7. 前記スキャン手段は、前記ＲＦパルスを含む第３のシーケンスを実行し、
   前記コイルは、前記第３のシーケンスにより発生した第３のＭＲ信号を受信し、
   前記データ生成手段は、前記各コイルエレメントで受信された第３のＭＲ信号に基づいて、前記コイルエレメントごとに、前記第３のＭＲ信号の位相と振幅の情報を含む第３のデータを生成し、
   前記作成手段は、前記第２のデータと前記第３のデータとに基づいて、前記第２のスペクトルを作成する、請求項６に記載の磁気共鳴装置。
8. 前記作成手段は、
   前記第２のデータを、水のピークを含む第４のスペクトルに変換するとともに、前記第３のデータを、代謝物のピークを含む第５のスペクトルに変換し、前記第４のスペクトルと前記第５のスペクトルとを差分することにより、前記第２のスペクトルを作成する、請求項７に記載の磁気共鳴装置。
9. 前記作成手段は、
   前記第２のデータに含まれる第２のＭＲ信号の位相を算出する位相算出手段と、
   前記位相算出手段により算出された位相に基づいて、前記第１のデータに含まれる第１のＭＲ信号の位相を補正する補正手段と、
   位相が補正された前記第１のデータを、前記第１のスペクトルに変換する変換手段と、
   を有する、請求項６〜８のうちのいずれか一項に記載の磁気共鳴装置。
10. 前記スキャン手段は、前記関心部位を単一の領域と見なした１回の第２のシーケンスを実行する、請求項１〜９のうちのいずれか一項に記載の磁気共鳴装置。
11. 前記係数算出手段は、
    前記第４のスペクトル又は前記第５のスペクトルのノイズを求め、前記ノイズで前記特性値を正規化し、正規化された特性値に基づいて前記係数を算出する、請求項８に記載の磁気共鳴装置。
12. 前記合成手段は、
    前記コイルエレメントごとに得られた前記第１のスペクトルに前記係数を乗算し、係数が乗算された前記第１のスペクトルを加算する、請求項１〜１１のうちのいずれか一項に記載の磁気共鳴装置。
13. 前記特性値は、水のピークのピーク値、又は水のピークの面積である、請求項１〜１２のうちのいずれか一項に記載の磁気共鳴装置。
14. 水および代謝物を含む関心部位からＭＲ信号を発生させるためのシーケンスを実行するスキャン手段であって、前記関心部位に含まれる水の信号を抑制するためのＲＦパルスを含む第１のシーケンスと、前記ＲＦパルスを含まない第２のシーケンスとを実行するスキャン手段と、
    前記第１のシーケンスにより発生した第１のＭＲ信号と、前記第２のシーケンスにより発生した第２のＭＲ信号とを受信する複数のコイルエレメントを備えたコイルと、
    各コイルエレメントで受信された第１のＭＲ信号に基づいて、前記コイルエレメントごとに、前記代謝物のピークを含む第１のスペクトルを生成するとともに、各コイルエレメントで受信された第２のＭＲ信号に基づいて、前記コイルエレメントごとに、前記水のピークを含む第２のスペクトルを生成するスペクトル生成手段と、
    を有する磁気共鳴装置に適用されるプログラムであって、
    前記第２のスペクトルの水のピークの特性値に基づいて、前記コイルエレメントごとに得られた前記第１のスペクトルの強度を調整するための係数を算出する係数算出処理と、
    前記係数を用いて前記第１のスペクトルの各々の強度を調整し、強度が調整された前記第１のスペクトルを合成する合成処理と、
    をコンピュータに実行させるためのプログラムであり、
    前記係数算出処理は、
    前記第２のスペクトルのノイズを求め、前記ノイズで前記特性値を正規化し、正規化された前記特性値に基づいて、前記係数を算出する、プログラム。