"特定の物質からの信号の抑制又は強調効果を変えて被検体から複数の磁気共鳴信号の周波数スペクトルを収集するスペクトル収集手段と、前記複数の周波数スペクトル間における前記特定の物質又は他の物質からの信号の強度の相違を表す指標に基づいて、前記特定の物質又は前記他の物質の共鳴周波数を求める共鳴周波数取得手段と、中心周波数が前記特定の物質又は前記他の物質の共鳴周波数に設定された高周波パルスを用いてイメージングを行うイメージング手段と、を備える磁気共鳴イメージング装置。前記スペクトル収集手段は、脂肪抑制パルスのフリップアングル又は脂肪抑制パルスを印加するか否かを変えることによって脂肪抑制効果が異なる2つの周波数スペクトルを収集するように構成され、前記共鳴周波数取得手段は、脂肪からの信号の強度の相違を表す指標が閾値処理により大きいと判定される場合には、脂肪抑制効果が小さい側の周波数スペクトルを用いて脂肪の共鳴周波数を求め、更に求めた前記脂肪の共鳴周波数に基づいて水の共鳴周波数を求める一方、前記指標が前記閾値処理により小さいと判定される場合には、脂肪抑制効果が大きい側の周波数スペクトルを用いて水の共鳴周波数を求めるように構成される請求項1記載の磁気共鳴イメージング装置。前記共鳴周波数取得手段は、前記脂肪からの信号の強度を表す曲線の積分値の差又は比を前記指標とするように構成される請求項2記載の磁気共鳴イメージング装置。前記共鳴周波数取得手段は、前記脂肪からの信号強度のピークの差又は比を前記指標とするように構成される請求項2記載の磁気共鳴イメージング装置。前記脂肪抑制効果が大きい側の周波数スペクトルを用いて前記水の共鳴周波数が所定の精度で求められたか否かを判定し、前記水の共鳴周波数が前記所定の精度で求められたと判定された場合には、前記脂肪抑制効果が小さい側の周波数スペクトルが収集されないように前記スペクトル収集手段を制御する制御手段を更に備える請求項2乃至4のいずれか1項に記載の磁気共鳴イメージング装置。前記スペクトル収集手段は、前記特定の物質と他の物質との縦緩和時間の相違に基づいて、前記特定の物質からの信号の抑制又は強調効果の度合いの異なる複数の周波数スペクトルを収集するように構成される請求項1乃至5のいずれか1項に記載の磁気共鳴イメージング装置。前記共鳴周波数取得手段は、脂肪抑制効果が異なる複数の周波数スペクトルの各最大値に対応する複数の周波数帯域の前記複数の周波数スペクトル間における複数の信号変化率及び前記複数の周波数帯域と前記複数の周波数帯域からそれぞれ所定の周波数だけ離れた複数の周波数帯域との間における前記複数の信号変化率の比に基づいて水の共鳴周波数を求めるように構成される請求項1乃至6のいずれか1項に記載の磁気共鳴イメージング装置。特定の物質からの信号を抑制又は強調して被検体から磁気共鳴信号の周波数スペクトルを収集するスペクトル収集手段と、前記周波数スペクトルと参照用の周波数スペクトルとの間における一致度又は不一致度を表す指標に基づいて、前記特定の物質又は他の物質の共鳴周波数を求める共鳴周波数取得手段と、中心周波数が前記特定の物質又は前記他の物質の共鳴周波数に設定された高周波パルスを用いてイメージングを行うイメージング手段と、を備える磁気共鳴イメージング装置。しかし、特に胸部におけるイメージングを行う場合のように撮像領域に脂肪成分が多いと、脂肪抑制効果が不十分となる。この結果、周波数スペクトルの形状が乱れ、精度よくピークに対応する周波数を検出することが困難となる。場合によっては、脂肪の共鳴周波数において水よりも大きなピークが出現する場合がある。このため、誤って脂肪の共鳴周波数を水の共鳴周波数として認識してしまう恐れがある。一方、水の共鳴周波数の検出及びRFパルスの中心周波数の設定は自動化できることが望ましい。従って、ユーザが周波数スペクトルを目視して水の共鳴周波数を検出し、RFパルスの中心周波数を水の共鳴周波数に合わせる場合はもちろん、自動的に水の共鳴周波数の検出及びRFパルスの中心周波数の設定を行う場合においても、確実に水の共鳴周波数を十分な精度で検出することが重要である。これは、水以外の物質の共鳴周波数にRFパルスの中心周波数を設定し、特定の物質からのMR信号を発生させて映像化するmolecularimagingを行う場合においても同様である。本発明は、より高精度に水等の信号収集対象となる物質の共鳴周波数を検出することが可能な磁気共鳴イメージング装置を提供することを目的とする。本発明の実施形態に係る磁気共鳴イメージング装置は、スペクトル収集手段、共鳴周波数取得手段及びイメージング手段を備える。スペクトル収集手段は、特定の物質からの信号の抑制又は強調効果を変えて被検体から複数の磁気共鳴信号の周波数スペクトルを収集する。共鳴周波数取得手段は、前記複数の周波数スペクトル間における前記特定の物質又は他の物質からの信号の強度の相違を表す指標に基づいて、前記特定の物質又は前記他の物質の共鳴周波数を求める。イメージング手段は、中心周波数が前記特定の物質又は前記他の物質の共鳴周波数に設定された高周波パルスを用いてイメージングを行う。本発明の実施形態に係る磁気共鳴イメージング装置について添付図面を参照して説明する。第1の実施形態図1は本発明の第1の実施形態に係る磁気共鳴イメージング装置の構成図である。磁気共鳴イメージング装置20は、静磁場を形成する筒状の静磁場用磁石21、この静磁場用磁石21の内部に設けられたシムコイル22、傾斜磁場コイル23及びRFコイル24を備えている。また、磁気共鳴イメージング装置20には、制御系25が備えられる。制御系25は、静磁場電源26、傾斜磁場電源27、シムコイル電源28、送信器29、受信器30、シーケンスコントローラ31及びコンピュータ32を具備している。制御系25の傾斜磁場電源27は、X軸傾斜磁場電源27x、Y軸傾斜磁場電源27y及びZ軸傾斜磁場電源27zで構成される。また、コンピュータ32には、入力装置33、表示装置34、演算装置35及び記憶装置36が備えられる。静磁場用磁石21は静磁場電源26と接続され、静磁場電源26から供給された電流により撮像領域に静磁場を形成させる機能を有する。尚、静磁場用磁石21は超伝導コイルで構成される場合が多く、励磁の際に静磁場電源26と接続されて電流が供給されるが、一旦励磁された後は非接続状態とされるのが一般的である。また、静磁場用磁石21を永久磁石で構成し、静磁場電源26が設けられない場合もある。また、静磁場用磁石21の内側には、同軸上に筒状のシムコイル22が設けられる。シムコイル22はシムコイル電源28と接続され、シムコイル電源28からシムコイル22に電流が供給されて静磁場が均一化されるように構成される。傾斜磁場コイル23は、X軸傾斜磁場コイル23x、Y軸傾斜磁場コイル23y及びZ軸傾斜磁場コイル23zで構成され、静磁場用磁石21の内部において筒状に形成される。傾斜磁場コイル23の内側には寝台37が設けられて撮像領域とされ、寝台37には被検体Pがセットされる。RFコイル24にはガントリに内蔵されたRF信号の送受信用の全身用コイルや寝台37や被検体P近傍に設けられるRF信号の受信用の局所コイルなどがある。また、傾斜磁場コイル23は、傾斜磁場電源27と接続される。傾斜磁場コイル23のX軸傾斜磁場コイル23x、Y軸傾斜磁場コイル23y及びZ軸傾斜磁場コイル23zはそれぞれ、傾斜磁場電源27のX軸傾斜磁場電源27x、Y軸傾斜磁場電源27y及びZ軸傾斜磁場電源27zと接続される。そして、X軸傾斜磁場電源27x、Y軸傾斜磁場電源27y及びZ軸傾斜磁場電源27zからそれぞれX軸傾斜磁場コイル23x、Y軸傾斜磁場コイル23y及びZ軸傾斜磁場コイル23zに供給された電流により、撮像領域にそれぞれX軸方向の傾斜磁場Gx、Y軸方向の傾斜磁場Gy、Z軸方向の傾斜磁場Gzを形成することができるように構成される。RFコイル24は、送信器29及び受信器30の少なくとも一方と接続される。送信用のRFコイル24は、送信器29からRF信号を受けて被検体Pに送信する機能を有し、受信用のRFコイル24は、被検体P内部の原子核スピンのRF信号による励起に伴って発生したMR信号を受信して受信器30に与える機能を有する。一方、制御系25のシーケンスコントローラ31は、傾斜磁場電源27、送信器29及び受信器30と接続される。シーケンスコントローラ31は傾斜磁場電源27、送信器29及び受信器30を駆動させるために必要な制御情報、例えば傾斜磁場電源27に印加すべきパルス電流の強度や印加時間、印加タイミング等の動作制御情報を記述したシーケンス情報を記憶する機能と、記憶した所定のシーケンスに従って傾斜磁場電源27、送信器29及び受信器30を駆動させることによりX軸傾斜磁場Gx、Y軸傾斜磁場Gy,Z軸傾斜磁場Gz及びRF信号を発生させる機能を有する。また、シーケンスコントローラ31は、受信器30におけるMR信号の検波及びAD変換により得られた複素データである生データを受けてコンピュータ32に与えるように構成される。このため、送信器29には、シーケンスコントローラ31から受けた制御情報に基づいてRF信号をRFコイル24に与える機能が備えられる一方、受信器30には、RFコイル24から受けたMR信号を検波して所要の信号処理を実行するとともにA/D変換することにより、デジタル化された複素データである生データを生成する機能と生成した生データをシーケンスコントローラ31に与える機能とが備えられる。図2は、図1に示すコンピュータ32の機能ブロック図である。コンピュータ32の演算装置35は、記憶装置36に保存されたプログラムを実行することにより撮像条件設定部40及びデータ処理部41として、記憶装置36はk空間データ記憶部42及び画像データ記憶部43として、それぞれ機能する。また、撮像条件設定部40は、パルスシーケンス設定部40A、周波数スペクトル生成部40B、ピーク検出部40C、検出精度判定部40D、差分値計算部40E、共鳴周波数検出部40F及び中心周波数設定部40Gを有する。撮像条件設定部40は、パルスシーケンスを含む撮像条件を設定し、設定した撮像条件をシーケンスコントローラ31に出力して制御する機能を有する。特に、撮像条件設定部40は、イメージング用の撮像条件に加え、MR信号の周波数スペクトルを収集するためのプレスキャン用の撮像条件を設定する機能を備えている。周波数スペクトルは、通常シムコイル電源28からシムコイル22に供給すべき電流を決定するシミング用のプレスキャンにおいて収集される。そして、収集された周波数スペクトルは、イメージング用の撮像条件の1つであるRFパルスの中心周波数の設定に用いられる。また、撮像条件設定部40は、水や脂肪等の特定の物質からの信号の抑制効果又は強調効果を変えて収集された互いに異なる複数の周波数スペクトル間における特定の物質又は他の物質からのMR信号の強度の相違を表す指標に基づいて、特定の物質又は他の物質の共鳴周波数を求める機能を備えている。パルスシーケンス設定部40Aは、イメージングスキャン用のパルスシーケンスの他、特定の物質からの信号の抑制効果又は強調効果を変えて被検体Pから互いに異なる複数のMR信号の周波数スペクトルを収集するプレスキャン用のパルスシーケンスを設定する機能を有する。以下、2つの周波数スペクトルを収集する場合について説明する。図3は、図2に示すパルスシーケンス設定部40Aにおいて設定されるシーケンスの第1の例を示す図である。図3において横軸は時間を示す。図3に示すように、脂肪抑制パルスの印加後にデータ収集を行う第1のシーケンスと、図3に示すように脂肪抑制パルスを印加せずにデータ収集を行う第2のシーケンスを用いて、2回プレスキャンを行うと、脂肪抑制効果の異なる2セットのMRデータを収集することができる。従って第1のシーケンスにより収集された第1のMRデータから第1の周波数スペクトルを生成し、第2のシーケンスにより収集された第2のMRデータから第2の周波数スペクトルを生成すると、第1の周波数スペクトルにおける脂肪信号の抑制効果は第2の周波数スペクトルにおける脂肪信号の抑制効果と異なる効果となる。図4は、図2に示すパルスシーケンス設定部40Aにおいて設定されるシーケンスの第2の例を示す図である。図4において横軸は時間を示す。図4に示すように、フリップ角がα°の脂肪抑制パルスの印加後にデータ収集を行う第1のシーケンスと、図4に示すようにFAがα°と異なるβ°の脂肪抑制パルスの印加後にデータ収集を行う第2のシーケンスを用いて、2回プレスキャンを行うことによっても、脂肪抑制効果の異なる2セットのMRデータを収集することができる。従って、脂肪抑制効果が異なる2つの周波数スペクトルを取得することができる。尚、図4に示すβ°をゼロとした場合が図3のケースに相当する。つまり、脂肪抑制パルスのFA又は脂肪抑制パルスを印加するか否かを変えることによって脂肪抑制効果が異なる2つの周波数スペクトルを収集することができる。図5は、図2に示すパルスシーケンス設定部40Aにおいて設定されるシーケンスの第3の例を示す図である。図5において横軸は時間を示す。図5に示すように脂肪抑制パルスの印加後に第1のMRデータ及び第2のMRデータを連続収集するマルチエコーデータ収集シーケンスによっても脂肪抑制効果が異なる2つの周波数スペクトルを取得することができる。すなわち、図5に示すシーケンスでプレスキャンを行うと、第1のMRデータ及び第2のMRデータは、脂肪抑制パルスの印加時刻から互いに異なる経過時間後において収集される。従って、第1のMRデータ及び第2のMRデータは、脂肪抑制効果が互いに異なるデータとなる。このため、第1及び第2のMRデータからそれぞれ第1及び第2の周波数スペクトルを生成すると、第1及び第2の周波数スペクトルは互いに異なる脂肪抑制効果を呈する。すなわち、脂肪等の特定の物質からの信号抑制又は信号強調に影響を与える撮像条件のパラメータを変えて2回MRデータを収集すれば、異なる信号抑制効果又は信号強調効果を示す2つの周波数スペクトルを収集することができる。従って、上述の例では、プリパルスとして印加される脂肪抑制パルスに関する条件を変えたが、他の物質の信号を抑制する抑制パルス又は強調する励起パルスの条件を変えたり、抑制パルス又は励起パルスの印加を伴わずに特定の物質を選択的に抑制又は強調する撮像法の条件を変えてもよい。脂肪抑制パルスとしては、短時間反転回復パルス、化学シフト選択パルス、SPIRパルス、SPAIRパルスなどが知られている。また、複数の同一又は異なる脂肪抑制パルスを印加する脂肪抑制法も知られている。更に、プリパルスを印加せずに水からの信号を強調するシーケンスとしてPASTAが知られている。この他、シリコーンや水からの信号を抑制する撮像法も知られている。多くのケースでは、イメージングスキャンにおけるRFパルスの中心周波数が水の共鳴周波数に合わせられる。この場合、周波数スペクトルは、水の共鳴周波数を検出するために収集される。従って、以降では、RFパルスの中心周波数を水の共鳴周波数に合わせるために、周波数スペクトルから水の共鳴周波数を検出する場合について説明する。水の共鳴周波数と誤認識する可能性が高いのは、水の共鳴周波数に対して相対的に負方向に35ppmだけ離れている脂肪の共鳴周波数である。従って、従来、周波数スペクトルを収集するプレスキャンでは、周波数スペクトルにおいて脂肪信号のピークが抑制されるように脂肪抑制パルスが印加される。尚、データ収集領域にシリコーンが存在すると、水の共鳴周波数に対して約5ppmだけ負側にシフトした周波数においてシリコーンからの信号のピークが周波数スペクトルに現れる。この場合には、シリコーン抑制パルスを印加してもよい。これに対して、パルスシーケンス設定部40Aでは、収集される2つの周波数スペクトルにおいて、水信号に相当するピークが明瞭かつ同等となる一方、脂肪信号に相当する分布形状がなるべく異なる形状となるようにプレスキャン用の撮像条件が設定される。そのためには、静磁場の不均一性の影響が少ないSTIRパルスを脂肪抑制パルスとして印加する第1のシーケンスとSTIRパルスを印加しない第2のシーケンスをプレスキャン用の撮像条件とすることが好適である。STIRパルスの印加を伴うSTIRシーケンスは、180°RFパルスとしてSTIRパルスを印加することによってデータ収集部位における縦磁化を一旦反転させ、縦緩和によって脂肪の縦磁化が回復してゼロになるタイミングで90°RF励起パルスを印加してMRデータを収集するスピンエコー法によるシーケンスである。すなわち、水のT1緩和時間は、脂肪のT1緩和時間よりも長いため、脂肪の縦磁化がゼロとなるタイミングで励起パルスを印加してデータを収集すれば、水信号を選択的に収集することができる。同様に、物質ごとにT1緩和時間が異なるため、特定の物質と他の物質とのT1緩和時間の相違に基づいて、特定の物質からのMR信号を抑制又は強調することができる。従って、STIRパルスの印加の有無やFA等の条件を変えることによって、特定の物質と他の物質とのT1緩和時間の相違に基づいて、特定の物質からの信号の抑制又は強調効果の度合いの異なる複数の周波数スペクトルを収集することができる。周波数スペクトル生成部40Bは、プレスキャンによって収集されたMRデータに基づいて周波数スペクトルを生成する機能を有する。従って、2つの周波数スペクトルを生成するための第1及び第2のMRデータが収集された場合には、第1及び第2の周波数スペクトルが生成される。プレスキャンによって収集されたMRデータは、データ処理部41から与えられる。ピーク検出部40Cは、周波数スペクトル生成部40Bにおいて生成された周波数スペクトルのピーク及びピークに対応する周波数を検出する機能を有する。また、ピーク検出部40Cは、複数のピークが検出された場合には、ピーク間のシフト量及び最大値を求めるように構成される。尚、ピーク検出部40Cが、最初に検出されたピークに対応する周波数を35ppmだけ正方向又は負方向にシフトさせた周波数において別のピークが検出されるか否かを判定するようにしてもよい。検出精度判定部40Dは、脂肪抑制を伴って収集された周波数スペクトルから検出されたピークの検出精度を判定し、検出精度が所定の精度であり十分であると判定された場合には、脂肪抑制を伴わずに周波数スペクトルを収集するプレスキャンを実行しないように撮像条件設定部40を制御する機能を有する。例えば、検出精度判定部40Dがパルスシーケンス設定部40Aに脂肪抑制を伴わずに周波数スペクトルを収集するプレスキャンのシーケンスを生成しないように制御するようにしてもよい。すなわち、脂肪抑制を伴って収集された周波数スペクトルが水の共鳴周波数において明瞭なピークを呈し、十分な精度で水の共鳴周波数を検出できる場合には、更に周波数スペクトルを収集しなくてもRFパルスの中心周波数を設定することができる。そこで、無用なプレスキャンが実行されないように検出精度判定部40Dによって撮像条件設定部40が制御される。特に脂肪の少ない部位からデータが収集された場合には、十分な精度で水の共鳴周波数を検出できる場合が多い。このため、検出精度判定部40Dにおいて、脂肪抑制効果が大きい側の周波数スペクトルを用いて水の共鳴周波数が所定の精度で求められたか否かを判定し、水の共鳴周波数が所定の精度で求められたと判定された場合には、脂肪抑制効果が小さい側の周波数スペクトルが収集されないように周波数スペクトルの収集を制御することができる。差分値計算部40Eは、周波数スペクトル生成部40Bにおいて生成された2つの周波数スペクトルの脂肪信号及び水信号の一方又は双方に相当する分布形状の相違を表す指標を計算する機能を有する。この指標としては、例えば2つの周波数スペクトル間における脂肪信号又は水信号に相当するピーク値の差又は比、ある周波数の範囲における信号の積分値の差又は比を用いることができる。図6は、図2に示す周波数スペクトル生成部40Bにおいて生成された2つの周波数スペクトルの第1の例を示す図である。図6において横軸は周波数fを示し、縦軸は周波数fにおける各信号強度S1S2を示す。STIRパルスの印加を伴う第1のシーケンスにより第1のMRデータを収集して第1の周波数スペクトルを生成すると、脂肪信号が抑制されているため、図6に示すように1つのピークが第1の周波数スペクトルに現れる。従って、ピーク検出部40Cにより検出された第1の周波数スペクトルの1つのピークが水の共鳴周波数であると容易に判定することができる。そして、ピーク検出部40Cにより1つのピークに対応する周波数が十分な精度で検出された場合には、上述したように検出精度判定部40Dにおいてピークの検出精度が十分であると判定される。ピークの検出精度が十分であるか否かの判定は、例えば、信号強度が閾値以上となる周波数の範囲が閾値以下であるか否かの判定や標準偏差等の統計的指標の閾値判定によって行うことができる。ただし、図6のように、1つのピークがピーク検出部40Cにより検出されたとしても、精度が十分でないと検出精度判定部40Dにおいて判定された場合には、例えばSTIRパルスの印加を伴わない第2のシーケンスにより第2のMRデータが収集される。そうすると、第2のMRデータには脂肪抑制効果がないため、第2のMRデータから第2の周波数スペクトルを生成すると、図6に示すように脂肪の共鳴周波数及び水の共鳴周波数にそれぞれ対応する2つのピークが第2の周波数スペクトルに現れる。脂肪の共鳴周波数は水の共鳴周波数に対して負側に35ppmだけシフトしている。このため、2つのピークに対応する周波数のうち高い周波数が水の共鳴周波数に相当し、低い周波数が脂肪の共鳴周波数に相当する。図6に示す例では、高い周波数側のピークの方が、信号強度が大きい。従って、少なくともピーク検出部40Cにより第1及び第2の周波数スペクトルの各最大値を検出し、差分値計算部40Eにより各最大値に対応する周波数が同位置となるように第1及び第2の周波数スペクトル間における相対的な位置合わせを行えば、第2の周波数スペクトルにおける2番目のピークを検出しなくても、図6に示すように第1及び第2の周波数スペクトルにおける水の共鳴周波数の相対的な位置を合わせることができる。このため、式に示すように、脂肪信号に相当するピーク値の差Dpを、第1及び第2の周波数スペクトルの最大値に対応する周波数から負側に35ppmだけ低い周波数における信号強度P1P2の差として算出することができる。DpP2P1また、第1及び第2の周波数スペクトルの最大値に対応する周波数から負側に35ppmだけ低い周波数よりも更に低い周波数f1及び高い周波数f2を所定の周波数差Δfとなるように設定すれば、式に示すように、脂肪信号S1S2の積分値の差Diを算出することができる。尚、図6に示すような周波数スペクトルが複数の位置において収集される場合には、信号の積分値は位置方向軸も有する信号分布を表す局面で囲まれた範囲の体積ボリュームとなり、1つの位置において周波数スペクトルが収集される場合には、信号の積分値は信号分布を表す曲線で囲まれた範囲の面積となる。従って、より高精度に積分値を求めるためには、前処理として第1及び第2の周波数スペクトルのサーフェスフィッティング又はカーブフィッティングを行うことが好適である。また、水信号に相当するピーク値の差及び積分値の差についても第1及び第2の周波数スペクトルの最大値に対応する周波数を基準として同様に計算することができる。但し、STIRパルスを印加するか否か等の脂肪抑制条件によって水信号の強度は大きく変化しない。従って、積分範囲Δfを明らかに脂肪及び水の共鳴周波数を含む範囲となるように十分に大きく設定すれば、差分によって水信号がキャンセルされる。このため、積分値の差を脂肪信号の積分値の差とみなすことができる。この場合、第1及び第2の周波数スペクトルにおける周波数の相対的な位置を合わせ処理は不要となる。図7は、図2に示す周波数スペクトル生成部40Bにおいて生成された2つの周波数スペクトルの第2の例を示す図である。図7において横軸は周波数fを示し、縦軸は周波数fにおける各信号強度S1S2を示す。STIRパルスの印加を伴う第1のシーケンスにより胸部等の脂肪成分の多い部位から第1のMRデータを収集して第1の周波数スペクトルを生成すると、脂肪信号が十分に抑制されずに、図7に示すように脂肪及び水の共鳴周波数にそれぞれ対応する不明瞭な2つのピークが相対的に35ppmだけシフトした周波数において現れる場合がある。また、STIRパルスの印加を伴わない第2のシーケンスにより脂肪成分の多い部位から第2のMRデータを収集して第2の周波数スペクトルを生成すると、図7に示すように脂肪及び水の共鳴周波数にそれぞれ対応する2つのピークが現れ、かつ脂肪の共鳴周波数において信号強度が最大となる場合がある。更に、図7に示すように第1の周波数スペクトルが水の共鳴周波数において最大の信号強度を呈する場合に第1及び第2の周波数スペクトルの各最大値に対応する周波数の相対的な位置合わせを行うと、図7に示すように第1の周波数スペクトルにおける水の共鳴周波数と第1の周波数スペクトルにおける脂肪の共鳴周波数とが対応してしまう。従って、このような場合には、ピーク検出部40Cにより第1及び第2の周波数スペクトルの2つのピーク及びピーク間の相対的位置を検出し、かつ差分値計算部40Eにより式により脂肪信号に相当するピーク値の差Dpを算出する前に前処理として35ppmに相当する周波数の相対的なシフト補正処理を第1及び第2の周波数スペクトルに対して行うことが必要となる。但し、図7に示すように、水と脂肪間における共鳴周波数のシフト量の2倍を超えるような広い積分範囲Δfを設定すれば、脂肪抑制の条件によって水信号の強度が大きく変化しないため、周波数の相対的なシフト補正処理を行わなくても脂肪信号の積分値の差を算出することができる。従って、少なくとも第1及び第2の周波数スペクトルの各最大値を検出すればよい。一方、周波数の相対的なシフト補正処理を行えば、積分範囲を脂肪信号のピークに対応する周波数範囲程度に設定することによって、脂肪信号の積分値の差をより高精度に算出することができる。加えて、水信号の積分値の差を算出することもできる。共鳴周波数検出部40Fは、差分値計算部40Eにおける指標の計算結果に基づいて、2つの周波数スペクトルのいずれに基づいて水の共鳴周波数を検出するのかを決定する機能と、決定した一方の周波数スペクトルに基づいて、必要に応じて脂肪の共鳴周波数を検出することによって、水の共鳴周波数を求める機能を有する。また、検出精度判定部40Dにおいて、水の共鳴周波数が必要な精度で検出されたか否かの判定を行う場合には、共鳴周波数検出部40Fが最初に収集された周波数スペクトルから検出されたピークに基づいて水の共鳴周波数を検出するように構成される。図6に示すように、脂肪抑制を伴って脂肪成分が少ない部位から収集された周波数スペクトルの場合には、水信号に対応するピークの方が脂肪信号に対応するピークよりも明瞭に出現する。このため、水信号のピークから水の共鳴周波数を求めることが精度上望ましい。しかし、図7に示すように、脂肪抑制を伴わずに脂肪成分が多い部位から収集された周波数スペクトルの場合には、脂肪信号に対応するピークの方が水信号に対応するピークよりも明瞭となる場合がある。この場合、脂肪信号に対応するピークから脂肪の共鳴周波数を求め、求めた脂肪の共鳴周波数よりも35ppmだけ高い周波数を水の共鳴周波数とみなす方が精度上望ましいと考えられる。更に、脂肪成分の量が中間的な場合には、水信号及び脂肪信号のいずれのピークを用いて水の共鳴周波数を求めれば良いのかという判断をユーザが目視で行うことが困難となり得る。そこで、共鳴周波数検出部40Fは、差分値計算部40Eにおける指標の計算結果に基づいて、画一的に水信号及び脂肪信号のいずれのピークを用いて水の共鳴周波数を求めるのかを決定することができるように構成される。この決定は、指標に対する閾値処理によって行うことができる。例えば、脂肪成分が多い部位では、STIRパルスの印加によって脂肪抑制効果が不十分とはなるものの、脂肪成分が多いため、STIRパルスを印加するか否かによって脂肪信号の強度が大きく変化する。一方、脂肪成分が少ない部位では、STIRパルスの印加によって脂肪信号が抑制されるものの、脂肪成分が多い部位に比べると、STIRパルスを印加するか否かによる脂肪信号の強度の変化量は小さい。また、水信号の強度に着目すると、脂肪成分が少ない部位では、STIRパルスを印加するか否かによる水信号の絶対的な強度の変化量も小さい。そこで、STIRパルスを印加して収集された第1の周波数スペクトルとSTIRパルスを印加せずに収集された第2の周波数スペクトルとの間における脂肪信号又は水信号の強度の差又は比を表す指標が閾値よりも小さい場合には、脂肪成分が少ない部位と判定し、STIRパルスを印加して収集された第1の周波数スペクトルのピークに基づいて水の共鳴周波数を算出するように決定することができる。逆に、第1及び第2の周波数スペクトル間における脂肪信号の強度の差又は比を表す指標が閾値よりも大きい場合には、脂肪成分が多い部位と判定し、STIRパルスを印加せずに収集された第2の周波数スペクトルにおける脂肪信号のピークに基づいて脂肪の共鳴周波数を算出するように決定することができる。この場合、算出した脂肪の共鳴周波数から水の共鳴周波数が算出される。つまり、共鳴周波数検出部40Fは、脂肪からの信号の強度の相違を表す指標が閾値処理により大きいと判定される場合には、脂肪抑制効果が小さい側の周波数スペクトルを用いて脂肪の共鳴周波数を求め、更に求めた脂肪の共鳴周波数に基づいて水の共鳴周波数を求める一方、指標が閾値処理により小さいと判定される場合には、脂肪抑制効果が大きい側の周波数スペクトルを用いて水の共鳴周波数を求めるように構成される。また、閾値処理の対象となる指標は、脂肪からの信号の強度を表す曲線の積分値の差又は比とすることができる。或いは、脂肪からの信号強度のピークの差又は比を指標とすることもできる。尚、周波数スペクトルのピークから共鳴周波数を求める方法としては、信号の極大値に対応する周波数を共鳴周波数とする方法の他、周波数スペクトルにおける信号強度分布の重心に対応する周波数を共鳴周波数とする方法や信号強度分布の半値幅の中心に対応する周波数を共鳴周波数とする方法などがある。中心周波数設定部40Gは、イメージング用のRFパルスの中心周波数を共鳴周波数検出部40Fにより算出された水の共鳴周波数に設定する機能を有する。データ処理部41は、シーケンスコントローラ31から生データを受けてk空間データ記憶部42に形成されたk空間に配置する機能、シーケンスコントローラ31から周波数スペクトルの生成用のMRデータを受けて周波数スペクトル生成部40Bに与える機能、k空間データ記憶部42からイメージング用のk空間データを取得して、画像再構成処理を含む必要なデータ処理を行うことによって画像データを生成する機能、生成した画像データを画像データ記憶部43に書き込む機能、所望の画像データを画像データ記憶部43から読み込んで必要な画像処理を施して表示装置34に表示させる機能を有する。以上では、脂肪抑制を行う場合を例に説明したが、水励起を行って相対的に脂肪信号を抑制する場合においても同様に水の共鳴周波数を求めることができる。但し、水励起の場合には、水励起の有無によって水信号の強度が大きく変化する。このため、周波数スペクトル間における水信号の強度分布の相違を表す指標を定義し、水信号についての指標に対する閾値処理によって水の共鳴周波数を求めるための周波数スペクトルを決定するようにしてもよい。また、上述の例では、周波数スペクトルを2つ収集する場合について説明したが、3つ以上の周波数スペクトルを収集してもよい。この場合、所望の方法で信号強度分布の相違を表す指標を定義し、指標に基づいて共鳴周波数の算出に用いる単一又は複数の周波数スペクトルを決定することができる。次に磁気共鳴イメージング装置20の動作及び作用について説明する。図8は、図1に示す磁気共鳴イメージング装置20によりRFパルスの中心周波数を水の共鳴周波数に合わせてイメージングを行う際の流れの一例を示すフローチャートである。ここでは、水の共鳴周波数を求めるために、STIRパルスの印加を伴う第1のシーケンス及びSTIRパルスの印加を伴わない第2のシーケンスにより第1及び第2の周波数スペクトルを収集する場合を例に説明する。まず、予め寝台37に被検体Pがセットされ、静磁場電源26により励磁された静磁場用磁石21の撮像領域に静磁場が形成される。また、シムコイル電源28からシムコイル22に電流が供給されて撮像領域に形成された静磁場が均一化される。そして、ステップS1において、STIRパルスの印加を伴う第1のシーケンスにより第1の周波数スペクトルが収集される。すなわち、撮像条件設定部40のパルスシーケンス設定部40Aは、図3に示すようなSTIRパルスの印加後にデータ収集を行う第1のシーケンスをプレスキャン用の撮像条件として設定する。そして、撮像条件設定部40から第1のシーケンスを含むプレスキャン用の撮像条件がシーケンスコントローラ31に出力される。そうすると、シーケンスコントローラ31は、プレスキャン用の撮像条件に従って傾斜磁場電源27、送信器29及び受信器30を駆動させることにより被検体Pがセットされた撮像領域に傾斜磁場を形成させるとともに、RFコイル24からRF信号を発生させる。このため、被検体Pの内部における核磁気共鳴により生じたMR信号が、RFコイル24により受信されて受信器30に与えられる。受信器30は、RFコイル24からMR信号を受けてデジタルデータのMR信号を生成し、シーケンスコントローラ31に与える。シーケンスコントローラ31はMR信号をデータ処理部41に与え、データ処理部41はMR信号を周波数スペクトル生成部40Bに与える。そして、周波数スペクトル生成部40Bは、第1のシーケンスにより収集された第1のMRデータに基づいて第1の周波数スペクトルを生成する。次にステップS2において、ピーク検出部40Cは、第1の周波数スペクトルのピーク及びピークに対応する周波数を検出する。検出されたピークが1つである場合には、共鳴周波数検出部40Fがピークに対応する周波数を水の共鳴周波数と認識する。一方、検出されたピークが2つであり、かつピーク間における周波数シフト量が35ppmである場合には、共鳴周波数検出部40Fが高い周波数側のピークに対応する周波数を水の共鳴周波数と認識する。次にステップS3において、検出精度判定部40Dは、水の共鳴周波数が必要な精度で検出されたか否かの判定を行う。そして水の共鳴周波数が必要な精度で検出されていないと判定された場合には、ステップS4において、STIRパルスの印加を伴わない第2のシーケンスにより第2の周波数スペクトルが収集される。すなわち、検出精度判定部40Dがパルスシーケンス設定部40Aに第2のシーケンスの生成指示を与える。そうすると、パルスシーケンス設定部40Aは、図3に示すようなSTIRパルスの印加を伴わない第2のシーケンスをプレスキャン用の撮像条件として設定する。そして、第1の周波数スペクトルの収集と同様な流れで第2の周波数スペクトルが収集される。また、第2の周波数スペクトルのピーク及びピークに対応する周波数がピーク検出部40Cにより検出される。次にステップS5において、差分値計算部40Eは、第1及び第2の周波数スペクトル間における脂肪信号強度の相違を表す指標を計算する。指標としては、式に示す脂肪信号に相当するピーク値P1P2の差Dp又は比、式に示す脂肪信号S1S2の積分値の差Di又は比を用いることができる。また、図7に示すように、第2の周波数スペクトルにおいて脂肪信号が最大値を呈する場合には、必要に応じて第1及び第2の周波数スペクトル間における周波数シフトが補正される。次にステップS6において、共鳴周波数検出部40Fは、脂肪信号強度の相違を表す指標の閾値判定を行う。そして、指標が閾値以上であると判定された場合には、撮像部位に脂肪成分が多く、脂肪信号に対応するピークの方が水信号に対応するピークよりも明瞭となると考えられる。このため、ステップS7において、共鳴周波数検出部40Fは、脂肪信号に対応するピークから脂肪の共鳴周波数を求める。次にステップS8において、共鳴周波数検出部40Fは、脂肪の共鳴周波数よりも35ppmだけ高い周波数を水の共鳴周波数として算出する。これにより、水の共鳴周波数が得られる。また、ステップS3の判定において、検出精度判定部40Dにより、水の共鳴周波数が必要な精度で検出されていると判定された場合には、検出精度判定部40Dがパルスシーケンス設定部40Aに第2のシーケンスを生成しないよう指示を与える。このため、第2の周波数スペクトルの収集は行われない。この場合、ステップS2においてピーク検出部40Cにより検出されたピークに対応する周波数が水の共鳴周波数と認識される。更に、ステップS6の判定において、共鳴周波数検出部40Fにより、脂肪信号強度の相違を表す指標が閾値以上でないと判定された場合には、撮像部位に脂肪成分が少なく、水信号に対応するピークの方が脂肪信号に対応するピークよりも明瞭となると考えられる。この場合にも、ステップS2においてピーク検出部40Cにより検出されたピークに対応する周波数が水の共鳴周波数と認識される。次にステップS9において、中心周波数設定部40Gは、イメージング用のRFパルスの中心周波数を水の共鳴周波数に設定する。また、パルスシーケンスを含む他のイメージング用の撮像条件が撮像条件設定部40により設定される。次にステップS10において、中心周波数が特定の物質又は他の物質の一例に相当する水の共鳴周波数に設定されたRFパルスを用いてイメージングが実行される。すなわち、撮像条件設定部40からイメージング用の撮像条件がシーケンスコントローラ31に出力される。これにより、プレスキャンと同様な流れでイメージング用のMRデータが収集され、収集されたMRデータがデータ処理部41によりk空間データ記憶部42に形成されたk空間に配置される。そして、データ処理部41におけるMRデータに対する画像再構成処理及び画像処理によって表示用の画像データが生成される。さらに生成された画像データが表示装置34に表示される。また、必要な画像データは画像データ記憶部43に保存される。つまり以上のような磁気共鳴イメージング装置20は、脂肪等の特定の物質からの抑制又は強調効果が互いに異なる条件でMR信号に関する複数の周波数スペクトルを収集し、特定の物質からの信号強度の周波数スペクトル間における相違を表す指標に応じて選択した周波数スペクトルを用いて特定の物質又は他の物質の共鳴周波数をRFパルスの中心周波数の設定用に検出するようにしたものである。例えば、水信号の共鳴周波数を求める場合であれば、脂肪信号の抑制効果が大きい周波数スペクトルと脂肪信号の抑制効果が小さい周波数スペクトルが収集され、脂肪信号の強度の相違が大きければ脂肪信号の抑制効果が小さい周波数スペクトルから脂肪信号の共鳴周波数が検出される。そして、脂肪信号と水信号間におけるケミカルシフト量に基づいて脂肪信号の共鳴周波数から水の共鳴周波数が算出される。一方、脂肪信号の強度の相違が小さければ脂肪信号の抑制効果が大きい周波数スペクトルから水信号の共鳴周波数が検出される。このため、磁気共鳴イメージング装置20によれば、水の共鳴周波数を求める場合に撮像部位において脂肪成分が多く、脂肪抑制効果が不十分であっても、より高精度に水の共鳴周波数を求めることができる。従って、脂肪成分が多い胸部におけるイメージングにおいて特に有効である。また、他の物質の共鳴周波数を求める場合においても高精度な検出を行うことができる。この結果、RFパルスの中心周波数をより適切な値に設定することが可能となる。特に、ケミカルシフトを利用して脂肪成分を抑制した画像を収集する場合において、良好に脂肪信号が抑制された画像を得ることができる。また、磁気共鳴イメージング装置20によれば、周波数スペクトルにおいて検出されたピークが水の共鳴周波数に対応するのか脂肪の共鳴周波数に対応するのかの判定精度を向上させ、共鳴周波数の自動検出及びRFパルスの中心周波数の自動設定を容易に行うことができる。第2の実施形態第2の実施形態における磁気共鳴イメージング装置は、共鳴周波数検出部40Fの機能が第1の実施形態における磁気共鳴イメージング装置と相違する。他の構成及び機能は第1の実施形態における磁気共鳴イメージング装置と同等であるため同符号を付して説明を省略し、主として共鳴周波数検出部40Fの機能について説明する。具体的には、第2の実施形態における磁気共鳴イメージング装置の共鳴周波数検出部40Fは、脂肪抑制効果が異なる複数の周波数スペクトルの各最大値に対応する複数の周波数帯域の、複数の周波数スペクトル間における複数の信号変化率及び複数の周波数帯域と複数の周波数帯域からそれぞれ所定の周波数だけ離れた複数の周波数帯域との間における複数の信号変化率の比に基づいて水の共鳴周波数を求める機能を有しいている。すなわち、共鳴周波数検出部40Fは、複数の周波数スペクトル間における水又は脂肪からのMR信号の強度の相違を表す指標として信号量の変化率を用いるように構成されている。そして、共鳴周波数検出部40Fは、複数の周波数スペクトル間及び所定の周波数だけ離れた複数の周波数帯域間における信号変化率に基づいて水の共鳴周波数を求めるように構成されている。また、必要に応じて共鳴周波数検出部40Fには、ノイズ判定処理を実行する機能が備えられる。ノイズ判定処理は、脂肪抑制効果が異なる複数の周波数スペクトルの各最大値に対応する周波数帯域から35ppmだけ低い周波数帯域における信号量の周波数スペクトル間における変化量に基づいて、信号量の変化率の計算にノイズ成分の信号量が用いられたか否かを判定し、ノイズ成分の信号量が用いられたと判定された場合には、信号量の変化率を定数に補正する処理である。以下、脂肪抑制パルスの一例としてSTIRパルスを印加して収集された第1の周波数スペクトルと、STIRパルスを印加せずに収集された第2の周波数スペクトルに基づいて水の共鳴周波数を求める場合を例に説明する。図4に示すように脂肪抑制パルスのFAを変えたり、図5に示すように脂肪抑制パルスの印加時刻からの経過時間を変えた場合も同様である。図9は、本発明の第2の実施形態に係る磁気共鳴イメージング装置において実行される水の共鳴周波数の検出処理の流れを示すフローチャートである。まずステップS20において、STIRパルスを印加して第1の周波数スペクトルが収集される。すなわち、図8のステップS1と同様に、STIRパルスの印加を伴うプレスキャンによって第1の周波数スペクトルの生成用のMR信号が収集される。そして、周波数スペクトル生成部40Bにおいて第1の周波数スペクトルが生成される。次に、ステップS21において、STIRパルスを印加せずに第2の周波数スペクトルが収集される。すなわち、図8のステップS4と同様にSTIRパルスの印加を伴わないプレスキャンによって第2の周波数スペクトルの生成用のMR信号が収集される。そして、周波数スペクトル生成部40Bにおいて第2の周波数スペクトルが生成される。尚、図8と同様に、水の共鳴周波数の検出精度が十分であるか否かの判定を行い、水の共鳴周波数の検出精度が十分である判定された場合には、第2の周波数スペクトルを収集しないようにしてもよい。次に、ステップS22において、ピーク検出部40Cは、第1及び第2の周波数スペクトルの最大値に対応する周波数f0STIRf0を検出する。或いは、第1及び第2の周波数スペクトルの最大値よりも35ppmだけ負側にそれぞれ2番目のピークに相当する極大値が検出された場合には、極大値に対応する周波数f0STIRf0を検出する。最大値に対応する周波数f0STIRf0を以後の処理に採用するか極大値に対応する周波数f0STIRf0を以後の処理に採用するかについては、任意のアルゴリズムに従って決定することができる。例えば、最大値の1Y以下の極大値が最大値に対応する周波数f0STIRf0よりも35ppmだけ負側の周波数f0STIRf0において検出された場合には、最大値に対応する周波数f0STIRf0を水の共鳴周波数と仮定して以後の処理に採用する一方、最大値の1Yよりも大きい極大値が最大値に対応する周波数f0STIRf0よりも35ppmだけ負側の周波数f0STIRf0において検出された場合には、極大値に対応する周波数f0STIRf0を脂肪の共鳴周波数と仮定して以後の処理に採用するというアルゴリズムが考えられる。次に、ステップS23において、共鳴周波数検出部40Fは、第1及び第2の周波数スペクトルの周波数軸上における位置を合わせる。図10は、図9に示すステップS20及びステップS21において収集され得る第1の周波数スペクトルと第2の周波数スペクトルの組合せを示す図である。図10及びの各上段は、STIRパルスの印加を伴って収集された第1の周波数スペクトルの例であり、各下段はSTIRパルスを印加せずに収集された第2の周波数スペクトルの例である。従って、図10における各横軸は周波数fを示している。また、図10及びは、それぞれ第1及び第2の周波数スペクトルの組合せとなり得る各ケースを示している。頭部のように脂肪が比較的少ない撮像部位についてSTIRパルスを印加せずに周波数スペクトルを収集すると、水の共鳴周波数において最大値となり、脂肪の共鳴周波数において最大値よりも小さい極大値を呈する第2の周波数スペクトルが得られる。この場合、STIRパルスを印加して周波数スペクトルを収集すると、第2の周波数スペクトルよりも脂肪の共鳴周波数において強度が小さくなった第1の周波数スペクトルが得られる。従って、第1の周波数スペクトル及び第2の周波数スペクトルは図10に示すようなCASE1の関係となる。図10に示すような第1の周波数スペクトル及び第2の周波数スペクトルの各最大値に対応する周波数f0STIRf0をそれぞれ検出すると、最大値に対応する周波数f0STIRf0は、いずれも水の共鳴周波数となる。一方、脂肪が存在する撮像部位についてSTIRパルスを印加せずに周波数スペクトルを収集すると、脂肪の共鳴周波数において最大値となり、水の共鳴周波数において最大値よりも小さい極大値を呈する第2の周波数スペクトルが得られる場合が多い。この場合、STIRパルスの印加によって脂肪抑制効果が十分に得られた状態で周波数スペクトルを収集すると、水の共鳴周波数において最大値となり、脂肪の共鳴周波数において最大値よりも小さい極大値を呈する第1の周波数スペクトルが得られる。従って、第1の周波数スペクトル及び第2の周波数スペクトルは図10に示すようなCASE2の関係となる。図10に示すような第1の周波数スペクトル及び第2の周波数スペクトルの各最大値に対応する周波数f0STIRf0をそれぞれ検出すると、第1の周波数スペクトルの最大値に対応する周波数f0STIRは水の共鳴周波数となるが、第2の周波数スペクトルの最大値に対応する周波数f0は脂肪の共鳴周波数となる。図10に示すCASE2は典型的なケースである。また、比較的脂肪が多い撮像部位についてSTIRパルスを印加せずに周波数スペクトルを収集すると、脂肪の共鳴周波数において最大値となり、水の共鳴周波数において最大値よりも小さい極大値を呈する第2の周波数スペクトルが得られる。この場合、STIRパルスを印加しても十分な脂肪抑制効果が得られず、脂肪の共鳴周波数において最大値となり、水の共鳴周波数において最大値よりも小さい極大値を呈する第1の周波数スペクトルが得られる。従って、第1の周波数スペクトル及び第2の周波数スペクトルは図10に示すようなCASE3の関係となる。図10に示すような第1の周波数スペクトル及び第2の周波数スペクトルの各最大値に対応する周波数f0STIRf0をそれぞれ検出すると、水の共鳴周波数と仮定された最大値に対応する周波数f0STIRf0は、実際にはいずれも脂肪の共鳴周波数となる。この他、殆ど生じないCASE4として図10に示すケースが考えられる。CASE4は、STIRパルスを印加して収集された第1の周波数スペクトルの最大値に対応する周波数f0STIRが脂肪の共鳴周波数に相当し、STIRパルスを印加せずに収集された第2の周波数スペクトルの最大値に対応する周波数f0が水の共鳴周波数に相当するケースである。次に、ステップS24において、共鳴周波数検出部40Fは、STIRパルスの印加を伴って収集された第1の周波数スペクトルの最大値に対応する周波数f0STIRを含む帯域幅W1Hzにおける信号量S1STIRを算出する。一方、STIRパルスの印加を伴わずに収集された第2の周波数スペクトルの帯域幅W1における信号量S1を算出する。更に、STIRパルスの印加を伴わずに収集された第2の周波数スペクトルの最大値に対応する周波数f0を含む帯域幅W1における第2の周波数スペクトルの信号量S1が共鳴周波数検出部40Fにより算出される。一方、STIRパルスの印加を伴って収集された第1の周波数スペクトルの帯域幅W1における信号量S1STIRが共鳴周波数検出部40Fにより算出される。各信号量S1STIRS1S1STIRS1は、図10に示すようにそれぞれ第1及び第2の周波数スペクトルを表す曲線と周波数軸とで囲まれる部分の面積として積分処理によって算出することができる。この場合、帯域幅W1W1は、第1及び第2の周波数スペクトルの半値幅等に設定することができる。或いは、第1及び第2の周波数スペクトルの最大値又は半値幅の中心における信号強度を各信号量S1STIRS1S1STIRS1としてもよい。この場合には、帯域幅W1W1は、単位幅となる。但し、各信号量S1STIRS1S1STIRS1を面積とすることが精度上好適である。従って、以下、各信号量S1STIRS1S1STIRS1を面積として求める場合を例に説明する。次に、ステップS25において、共鳴周波数検出部40Fは、第1の周波数スペクトルの最大値に対応する周波数f0STIRから35ppmだけ負側の周波数を含む帯域幅W2における信号量S2STIRを算出する。一方、STIRパルスの印加を伴わずに収集された第2の周波数スペクトルの帯域幅W2における信号量S2を算出する。更に、STIRパルスの印加を伴わずに収集された第2の周波数スペクトルの最大値に対応する周波数f0から35ppmだけ負側の周波数を含む帯域幅W2における第2の周波数スペクトルの信号量S2が共鳴周波数検出部40Fにより算出される。一方、STIRパルスの印加を伴って収集された第1の周波数スペクトルの帯域幅W2における信号量S2STIRが共鳴周波数検出部40Fにより算出される。尚、各信号量S2STIRS2S2STIRS2については、35ppmだけ正側における各信号量S1STIRS1S1STIRS1と同様に算出することができる。帯域幅W2W2についても、35ppmだけ正側における帯域幅W1W1と同様に決定することができる。また、第1及び第2の周波数スペクトルの最大値に対応する周波数f0STIRf0の代わりに、第1及び第2の周波数スペクトルの2番目の極大値に対応する周波数f0STIRf0が検出されている場合には、2番目の極大値に対応する周波数f0STIRf0を基準として周波数が低い側における各信号量S2STIRS2S2STIRS2が算出される一方、極大値に対応する周波数f0STIRf0から35ppmだけ正側の周波数を基準として周波数が高い側における各信号量S1STIRS1S1STIRS1が算出される。次に、ステップS26において、共鳴周波数検出部40Fは、各帯域幅W1W1W2W2における第1及び第2の周波数スペクトル間の信号量の変化率R1R2R1R2を算出する。各変化率R1R2R1R2は、それぞれ式、式、式及び式により計算することができる。R1S1S1STIRR2S2S2STIRR1S1S1STIRR2S2S2STIR次に、ステップS27において、共鳴周波数検出部40Fは、周波数が低い側における各信号量S2STIRS2S2STIRS2の、第1及び第2の周波数スペクトル間における変化量ABを算出する。各信号量S2STIRS2S2STIRS2の変化量ABは、それぞれ式及び式により面積間の差分値として計算することができる。AS2STIRS2BS2STIRS2次に、ステップS28において、共鳴周波数検出部40Fは、ノイズ判定処理を実行する。ノイズ判定処理は、周波数が低い側における信号量の変化率R2R2の計算に用いられた各信号量S2STIRS2S2STIRS2がノイズ成分の信号量であるか否かを判定し、ノイズ成分の信号量が用いられたと判定された場合には、信号量の変化率R2R2を定数に設定する処理である。このノイズ判定処理は、周波数が低い側における各信号量S2STIRS2S2STIRS2の変化量ABに基づくノイズ判定アルゴリズムにより行うことができる。ノイズ判定アルゴリズムは、ノイズ成分の信号量を用いて算出されたと判定された信号量の変化率R2R2を強制的に定数C1とする処理である。このノイズ判定アルゴリズムは、例えば式のように定めることができる。R2C1ABεかつTεである場合には、第1及び第2の周波数スペクトルが図10に示すCASE1に該当すると判定することができる。従って、STIRパルスの印加を伴って収集された第1の周波数スペクトルの最大値に対応する周波数f0STIRを水の共鳴周波数として特定することができる。次に、閾値判定の結果が、TεかつTεであれば、Tの閾値判定を行わなくても、STIRパルスの印加を伴って収集された第1の周波数スペクトルの最大値に対応する周波数f0STIRを水の共鳴周波数として特定することが可能である。一方、閾値判定の結果が、Tεである場合には、第1及び第2の周波数スペクトルが図10に示すCASE4に該当すると判定することができる。尚、CASE4は、STIRパルスの印加を伴って収集された第1の周波数スペクトルでは脂肪の共鳴周波数において最大値となる一方、STIRパルスを印加せずに収集された第2の周波数スペクトルでは水の共鳴周波数において最大値となる稀なケースである。この場合には、STIRパルスの印加を伴って収集された第1の周波数スペクトルの最大値に対応する周波数f0STIRが脂肪の共鳴周波数に相当するため、35ppmだけ高い周波数を水の共鳴周波数として特定することができる。このような第2の実施形態における磁気共鳴イメージング装置によれば、撮像部位に脂肪が殆ど存在しない場合や水からのMR信号が殆ど収集されない場合など、様々なケースにおいて周波数スペクトルのピークが水の共鳴周波数に対応するのか或いは脂肪の共鳴周波数に対応するのかを高精度に判定することが可能である。すなわち、周波数スペクトルに基づいて、水及び脂肪の共鳴周波数を高精度で検出することができる。この結果、イメージング用のRFパルスの中心周波数を適切に設定するとともに、良好に脂肪抑制効果が得られたMR画像を収集することが可能となる。特に、水の共鳴周波数を特定するために用いられる判定パラメータTTの値が限りなく0に近い値になることをノイズ判定処理の実行によって回避することができる。このため、信号変化率を用いた判定パラメータTTを安定的な値とし、判定パラメータTTに基づく水の共鳴周波数の特定処理についてもエラーを回避して安定的に行うことができる。第3の実施形態第3の実施形態における磁気共鳴イメージング装置は、共鳴周波数検出部の機能が第1及び第2の実施形態における磁気共鳴イメージング装置と相違する。他の構成及び機能は第1及び第2の実施形態における磁気共鳴イメージング装置と同等であるため説明を省略し、共鳴周波数検出部の機能についてのみ説明する。上述した第1及び第2の実施形態では、実際に収集された周波数スペクトルに基づいて水や脂肪等の特定の物質の共鳴周波数を求める方法について説明したが、参照用の周波数スペクトルの形状を特定の物質からの信号の抑制効果又は強調効果の度合いごとに予めモデル化し、モデル化した参照用の周波数スペクトルと実際に収集された周波数スペクトルとに基づいて特定の物質の共鳴周波数を求めるようにすることもできる。具体的には、特定の物質からの信号を抑制又は強調して被検体からMR信号の周波数スペクトルを収集し、実際に収集された周波数スペクトルと参照用の周波数スペクトルとの一致度又は不一致度を表す指標に基づいて、特定の物質又は他の物質の共鳴周波数を求めることができる。換言すれば、カーブフィッティングによって確率論的に水や脂肪等の特定の物質の共鳴周波数を求めることができる。フィッティングの手法としては、相関係数を一致度の指標として相関係数がより1に近づくように周波数スペクトルをシフトさせる方法や2乗誤差を不一致度の指標として最小とする最小2乗法によって周波数スペクトルをシフトさせる方法などが挙げられる。すなわち、実際に収集された周波数スペクトルと参照用の周波数スペクトルとの一致度が最大となるように、実際に収集された周波数スペクトルを周波数方向にシフトさせることによって、実際に収集された周波数スペクトルから検出された各ピークがどの物質の共鳴周波数に相当するのかを自動判定することができる。そこで、図2に示す構成であれば、共鳴周波数検出部40Dに、上述したような参照用の周波数スペクトルと実際に収集された周波数スペクトルとの一致度又は不一致度に基づく共鳴周波数の自動判定機能を設ければよい。他の実施形態以上、特定の実施形態について記載したが、記載された実施形態は一例に過ぎず、発明の範囲を限定するものではない。ここに記載された新規な方法及び装置は、様々な他の様式で具現化することができる。また、ここに記載された方法及び装置の様式において、発明の要旨から逸脱しない範囲で、種々の省略、置換及び変更を行うことができる。添付された請求の範囲及びその均等物は、発明の範囲及び要旨に包含されているものとして、そのような種々の様式及び変形例を含んでいる。"