(57)【特許請求の範囲】
【請求項１】
プロトンを含む第１の物質とプロトンを含む第２の物質との間においてCEST（化学交換飽和移動）により生じる磁化の移動を反映した情報を取得するための磁気共鳴装置であって、
ＲＦパルスを有する複数のシーケンスを実行するスキャン手段であって、前記ＲＦパルスの周波数が前記シーケンスごとに異なるように設定された複数のシーケンスを実行するスキャン手段と、
前記複数のシーケンスにより得られたデータに基づいて、CESTの影響を受けた信号成分を表すCEST成分と、CESTの影響を受けていない信号成分を表すベースライン成分とを含むＺスペクトルを作成するスペクトル作成手段と、
Ｚスペクトルを第１のスペクトルに変換するスペクトル変換手段であって、CESTの影響が現れる周波数において、前記第１のスペクトルのCEST成分の信号値とベースライン成分の信号値との比が、前記ＺスペクトルのCEST成分の信号値とベースライン成分の信号値との比よりも大きい前記第１のスペクトルが得られるように、前記Ｚスペクトルを前記第１のスペクトルに変換するスペクトル変換手段と、
前記第１のスペクトルに含まれるCEST成分の情報を求める手段と、
を有する磁気共鳴装置。
【請求項２】
プロトンを含む第１の物質とプロトンを含む第２の物質との間においてCEST（化学交換飽和移動）により生じる磁化の移動を反映した情報を取得するための磁気共鳴装置であって、
複数のＲＦパルスを含むパルスセットを有する複数のシーケンスを実行するスキャン手段であって、前記複数のＲＦパルスのうちのｐ番目のＲＦパルスの位相とｐ＋１番目のＲＦパルスの位相との位相差が前記シーケンスごとに異なるように、前記複数のＲＦパルスの位相をサイクルさせ、前記複数のシーケンスを実行するスキャン手段と、
前記複数のシーケンスにより得られたデータに基づいて、CESTの影響を受けた信号成分を表すCEST成分と、CESTの影響を受けていない信号成分を表すベースライン成分とを含むＺスペクトルを作成するスペクトル作成手段と、
Ｚスペクトルを第１のスペクトルに変換するスペクトル変換手段であって、CESTの影響が現れる位相差において、前記第１のスペクトルのCEST成分の信号値とベースライン成分の信号値との比が、前記ＺスペクトルのCEST成分の信号値とベースライン成分の信号値との比よりも大きい前記第１のスペクトルが得られるように、前記Ｚスペクトルを前記第１のスペクトルに変換するスペクトル変換手段と、
前記第１のスペクトルに含まれるCEST成分の情報を求める手段と、
を有する磁気共鳴装置。
【請求項３】
前記位相差を周波数に変換する周波数変換手段を有する、請求項２に記載の磁気共鳴装置。
【請求項４】
前記CEST成分の情報を求める手段は、
CEST成分を表すCEST項と、ベースライン成分を表すベースライン項とを含む第１の関数に基づいて、前記第１のスペクトルに含まれるCEST成分の情報を求める、請求項１又は３に記載の磁気共鳴装置。
【請求項５】
前記CEST項は複数の係数を含んでおり、
前記CEST成分の情報を求める手段は、
前記第１の関数を用いて前記第１のスペクトルをフィッティングすることにより、前記CEST項に含まれる複数の係数の各々の値を計算する第１のフィッティング手段を有する、請求項４に記載の磁気共鳴装置。
【請求項６】
前記第１のフィッティング手段は、
前記CEST項に含まれる複数の係数の各々の初期値を設定し、前記初期値を基準にして係数の値を変更することによりフィッティングを行う、請求項５に記載の磁気共鳴装置。
【請求項７】
前記第１のスペクトルから、CESTの影響を受けた信号成分が現れる周波数を検出する検出手段を有する、請求項６に記載の磁気共鳴装置。
【請求項８】
前記CEST項に含まれる複数の係数のうちの一つの係数は、CESTの影響を受けた信号成分が現れる周波数を表しており、
前記第１のフィッティング手段は、
前記検出手段により検出された周波数の値を、前記一つの係数の初期値として設定する、請求項７に記載の磁気共鳴装置。
【請求項９】
前記ＺスペクトルからCEST成分が除去された第２のスペクトルを求め、前記第２のスペクトルに基づいて、前記第１のスペクトルに含まれるベースライン成分の情報を求める手段を有する、請求項８に記載の磁気共鳴装置。
【請求項１０】
前記ベースライン項は複数の係数を含み、
前記ベースライン成分の情報を求める手段は、前記ベースライン項に含まれる複数の係数の各々の値を計算する、請求項９に記載の磁気共鳴装置。
【請求項１１】
前記ベースライン成分の情報を求める手段は、
前記第２のスペクトルを作成する他のスペクトル作成手段と、
複数の係数を含む第２の関数を用いて前記第２のスペクトルをフィッティングすることにより、前記第２の関数に含まれる複数の係数の各々の値を計算する第２のフィッティング手段と、
前記第２のフィッティング手段により計算された係数の値に基づいて、前記ベースライン項に含まれる複数の係数の各々の値を計算する係数値計算手段と、
を有する、請求項１０に記載の磁気共鳴装置。
【請求項１２】
前記第２の関数は、定数項に含まれる係数と、ローレンツ関数の項に含まれる係数とを有している、請求項１１に記載の磁気共鳴装置。
【請求項１３】
前記第２のフィッティング手段は、
前記第２の関数に含まれる複数の係数の各々の初期値を設定し、前記初期値を基準にして係数の値を変更することによりフィッティングを行う、請求項１１又は１２に記載の磁気共鳴装置。
【請求項１４】
前記CEST項と前記ベースライン項との和で表される第３のスペクトルを計算するスペクトル計算手段を有する、請求項１３に記載の磁気共鳴装置。
【請求項１５】
前記第３のスペクトルに基づいて、前記第１のスペクトルに他のCEST成分が含まれているか否かを判断する判断手段を有する、請求項１４に記載の磁気共鳴装置。
【請求項１６】
前記判断手段が、前記第１のスペクトルに他のCEST成分が含まれていると判断した場合、前記第１のフィッティング手段は、前記第１の関数に他のCEST成分に対応するCEST項を追加し、CEST項が追加された第１の関数を用いて前記第１のスペクトルをフィッティングすることにより、追加されたCEST項に含まれる複数の係数の各々の値を計算する、請求項１５に記載の磁気共鳴装置。
【請求項１７】
前記第１の関数に含まれるCEST項の数を表すｎの値を設定する設定手段を有し、
前記判断手段が、前記第１のスペクトルに他のCEST成分が含まれていると判断した場合、前記設定手段は、CEST項の数を表すｎの値をインクリメントする、請求項１６に記載の磁気共鳴装置。
【請求項１８】
前記第２のフィッティング手段は、
ｎ＝ｉ−１において、前記第２の関数に含まれる複数の係数の各々の値を計算するための第１のフィッティングを行い、
ｎ＝ｉにおいて、前記第２の関数に含まれる複数の係数の各々の値を計算するための第２のフィッティングを行う、請求項１７に記載の磁気共鳴装置。
【請求項１９】
前記第２のフィッティング手段は、
前記第２のフィッティングを行う場合、前記第１のフィッティングにより計算された各係数の値を、前記第２のフィッティングにおける各係数の初期値として設定する、請求項１８に記載の磁気共鳴装置。
【請求項２０】
前記第１のスペクトルに含まれるCEST成分の総数をカウントするカウント手段を有する、請求項１９に記載の磁気共鳴装置。
【請求項２１】
前記第１のフィッティング手段は、
前記ベースライン項に含まれる複数の係数の値を固定し、CEST項に含まれる係数の値を再計算するための第３のフィッティングを行う、請求項２０に記載の磁気共鳴装置。
【請求項２２】
前記第１のフィッティング手段は、
前記第３のフィッティングを行う場合、前記第３のフィッティングを行う前に計算されたCEST項の各係数の値を、CEST項の係数の初期値として設定する、請求項２１に記載の磁気共鳴装置。
【請求項２３】
前記第１のフィッティング手段が、CEST項に含まれる係数の値を再計算した後、前記カウント手段は、前記第１のスペクトルに含まれるCEST成分の総数をカウントする、請求項２２に記載の磁気共鳴装置。
【請求項２４】
CEST項に含まれる係数の値を再計算した後にカウントされたCEST成分の総数が、CEST項に含まれる係数の値を再計算する前にカウントされたCEST成分の総数と異なる場合、CEST項に含まれる係数の値の計算と、ベースライン項に含まれる係数の値の計算をやり直す、請求項２３に記載の磁気共鳴装置。
【請求項２５】
前記CEST項に含まれる複数の係数の値と、前記ベースライン項に含まれる複数の係数の値とに基づいて、理想的なＺスペクトルを推定するスペクトル推定手段と、
前記理想的なＺスペクトルと前記Ｚスペクトルをと比較し、前記理想的なＺスペクトルによって前記Ｚスペクトルが再現されているか否かを判断するスペクトル比較手段と、
を有する、請求項２４に記載の磁気共鳴装置。
【請求項２６】
前記スペクトル比較手段が、前記理想的なＺスペクトルによって前記Ｚスペクトルが再現されていると判断した場合、前記カウント手段は、前記第１のスペクトルに含まれるCEST成分の総数をカウントする、請求項２５に記載の磁気共鳴装置。
【請求項２７】
前記第１の関数のベースライン項に含まれる係数の個数は、前記第２の関数に含まれる係数の個数よりも少ない、請求項１１〜請求項２６のうちのいずれか一項に記載の磁気共鳴装置。
【請求項２８】
前記ＲＦパルスは、連続波のＲＦパルスである、請求項１に記載の磁気共鳴装置。
【請求項２９】
前記シーケンスは、複数のプリパレーションパルスを有しており、
前記複数のプリパレーションパルスの各々は、前記ＲＦパルスと、縦磁化を定常状態にするためのキラー勾配パルスとを含む、請求項１に記載の磁気共鳴装置。
【請求項３０】
前記スペクトル変換手段は、
Ｚスペクトルを第１のスペクトルに変換する変換式に基づいて、第１のスペクトルを得る、請求項１、３〜２９のうちのいずれか一項に記載の磁気共鳴装置。
【請求項３１】
前記第１の物質は水であり、
前記変換式は、水の共鳴周波数からのずれを表すオフセット周波数の値に応じて変化する偶関数を含む、請求項３０に記載の磁気共鳴装置。
【請求項３２】
前記オフセット周波数の値がゼロの場合、前記偶関数の値はゼロである、請求項３１に記載の磁気共鳴装置。
【請求項３３】
プロトンを含む第１の物質とプロトンを含む第２の物質との間においてCEST（化学交換飽和移動）により生じる磁化の移動を反映した情報を取得するための磁気共鳴装置であって、ＲＦパルスを有する複数のシーケンスを実行するスキャン手段を有し、前記ＲＦパルスの周波数が前記シーケンスごとに異なるように設定された複数のシーケンスを実行する磁気共鳴装置に適用されるプログラムであって、
前記複数のシーケンスにより得られたデータに基づいて、CESTの影響を受けた信号成分を表すCEST成分と、CESTの影響を受けていない信号成分を表すベースライン成分とを含むＺスペクトルを作成するスペクトル作成処理と、
Ｚスペクトルを第１のスペクトルに変換するスペクトル変換処理であって、CESTの影響が現れる周波数において、前記第１のスペクトルのCEST成分の信号値とベースライン成分の信号値との比が、前記ＺスペクトルのCEST成分の信号値とベースライン成分の信号値との比よりも大きい前記第１のスペクトルが得られるように、前記Ｚスペクトルを前記第１のスペクトルに変換するスペクトル変換処理と、
前記第１のスペクトルに含まれるCEST成分の情報を求める処理と、
をコンピュータに実行させるためのプログラム。
【請求項３４】
プロトンを含む第１の物質とプロトンを含む第２の物質との間においてCEST（化学交換飽和移動）により生じる磁化の移動を反映した情報を取得するための磁気共鳴装置であって、複数のＲＦパルスを含むパルスセットを有する複数のシーケンスを実行するスキャン手段を有し、前記複数のＲＦパルスのうちのｐ番目のＲＦパルスの位相とｐ＋１番目のＲＦパルスの位相との位相差が前記シーケンスごとに異なるように、前記複数のＲＦパルスの位相をサイクルさせ、前記複数のシーケンスを実行する磁気共鳴装置に適用されるプログラムであって、
前記複数のシーケンスにより得られたデータに基づいて、CESTの影響を受けた信号成分を表すCEST成分と、CESTの影響を受けていない信号成分を表すベースライン成分とを含むＺスペクトルを作成するスペクトル作成処理と、
Ｚスペクトルを第１のスペクトルに変換するスペクトル変換処理であって、CESTの影響が現れる位相差において、前記第１のスペクトルのCEST成分の信号値とベースライン成分の信号値との比が、前記ＺスペクトルのCEST成分の信号値とベースライン成分の信号値との比よりも大きい前記第１のスペクトルが得られるように、前記Ｚスペクトルを前記第１のスペクトルに変換するスペクトル変換処理と、
前記第１のスペクトルに含まれるCEST成分の情報を求める処理と、
をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｚスペクトルを得る磁気共鳴装置、およびその磁気共鳴装置に適用されるプログラムに関する。
【背景技術】
【０００２】
近年、低濃度の化合物を観測する方法として、化学交換により生じる信号減衰を利用したCEST（Chemical Exchange Saturation Transfer）イメージング法が注目されている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特表２０１２−５１３２３９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
CESTイメージングでは、ＲＦパルスの周波数を変更しながらシーケンスを実行し、シーケンスにより得られたデータに基づいて、Ｚスペクトルを作成する。Ｚスペクトルの信号値は、CEST効果が現れる周波数の付近において減衰する。したがって、Ｚスペクトルから、どの周波数において信号が減衰しているかを特定することにより、化合物の濃度や磁化交換速度に依存したCESTの情報を取得することが可能となる。
【０００５】
しかし、Ｚスペクトルは、水の共鳴周波数の付近において、自由水のローレンツ分布を表す下向きのピークが支配的である。このピークのピーク幅は送信磁場Ｂ１に比例しており、ピークの値は０近くにまで低下する。したがって、CESTにより生じる信号の低下が、自由水のローレンツ分布を表す下向きのピークに埋もれてしまい、Ｚスペクトルから、CESTの影響を反映した情報を取り出すことが難しいという問題がある。
【０００６】
したがって、CESTの情報を取得するのに適したスペクトルを得るための技術が望まれている。
【課題を解決するための手段】
【０００７】
本発明の第１の観点は、プロトンを含む第１の物質とプロトンを含む第２の物質との間においてCEST（化学交換飽和移動）により生じる磁化の移動を反映した情報を取得するための磁気共鳴装置であって、
ＲＦパルスを有する複数のシーケンスを実行するスキャン手段であって、前記ＲＦパルスの周波数が前記シーケンスごとに異なるように設定された複数のシーケンスを実行するスキャン手段と、
前記複数のシーケンスにより得られたデータに基づいて、CESTの影響を受けた信号成分を表すCEST成分と、CESTの影響を受けていない信号成分を表すベースライン成分とを含むＺスペクトルを作成するスペクトル作成手段と、
Ｚスペクトルを第１のスペクトルに変換するスペクトル変換手段であって、CESTの影響が現れる周波数において、前記第１のスペクトルのCEST成分の信号値とベースライン成分の信号値との比が、前記ＺスペクトルのCEST成分の信号値とベースライン成分の信号値との比よりも大きくなるように、前記Ｚスペクトルを前記第１のスペクトルに変換するスペクトル変換手段と、
前記第１のスペクトルに含まれるCEST成分の情報を求める手段と、
を有する磁気共鳴装置である。
【０００８】
本発明の第２の観点は、プロトンを含む第１の物質とプロトンを含む第２の物質との間においてCEST（化学交換飽和移動）により生じる磁化の移動を反映した情報を取得するための磁気共鳴装置であって、
前記複数のＲＦパルスを含むパルスセットを有する複数のシーケンスを実行するスキャン手段であって、前記複数のＲＦパルスのうちのｐ番目のＲＦパルスの位相とｐ＋１番目のＲＦパルスの位相との位相差が前記シーケンスごとに異なるように、前記複数のＲＦパルスの位相をサイクルさせ、前記複数のシーケンスを実行するスキャン手段と、
複数のシーケンスにより得られたデータに基づいて、CESTの影響を受けた信号成分を表すCEST成分と、CESTの影響を受けていない信号成分を表すベースライン成分とを含むＺスペクトルを作成するスペクトル作成手段と、
Ｚスペクトルを第１のスペクトルに変換するスペクトル変換手段であって、CESTの影響が現れる位相差において、前記第１のスペクトルのCEST成分の信号値とベースライン成分の信号値との比が、前記ＺスペクトルのCEST成分の信号値とベースライン成分の信号値との比よりも大きくなるように、前記Ｚスペクトルを前記第１のスペクトルに変換するスペクトル変換手段と、
前記第１のスペクトルに含まれるCEST成分の情報を求める手段と、
を有する磁気共鳴装置である。
【０００９】
本発明の第３の観点は、プロトンを含む第１の物質とプロトンを含む第２の物質との間においてCEST（化学交換飽和移動）により生じる磁化の移動を反映した情報を取得するための磁気共鳴装置であって、ＲＦパルスを有する複数のシーケンスを実行するスキャン手段を有し、前記ＲＦパルスの周波数が前記シーケンスごとに異なるように設定された複数のシーケンスを実行する磁気共鳴装置に適用されるプログラムであって、
前記複数のシーケンスにより得られたデータに基づいて、CESTの影響を受けた信号成分を表すCEST成分と、CESTの影響を受けていない信号成分を表すベースライン成分とを含むＺスペクトルを作成するスペクトル作成処理と、
Ｚスペクトルを第１のスペクトルに変換するスペクトル変換処理であって、CESTの影響が現れる周波数において、前記第１のスペクトルのCEST成分の信号値とベースライン成分の信号値との比が、前記ＺスペクトルのCEST成分の信号値とベースライン成分の信号値との比よりも大きくなるように、前記Ｚスペクトルを前記第１のスペクトルに変換するスペクトル変換処理と、
前記第１のスペクトルに含まれるCEST成分の情報を求める処理と、
をコンピュータに実行させるためのプログラムである。
【００１０】
本発明の第４の観点は、プロトンを含む第１の物質とプロトンを含む第２の物質との間においてCEST（化学交換飽和移動）により生じる磁化の移動を反映した情報を取得するための磁気共鳴装置であって、複数のＲＦパルスを含むパルスセットを有する複数のシーケンスを実行するスキャン手段を有し、前記複数のＲＦパルスのうちのｐ番目のＲＦパルスの位相とｐ＋１番目のＲＦパルスの位相との位相差が前記シーケンスごとに異なるように、前記複数のＲＦパルスの位相をサイクルさせ、前記複数のシーケンスを実行する磁気共鳴装置に適用されるプログラムであって、
前記複数のシーケンスにより得られたデータに基づいて、CESTの影響を受けた信号成分を表すCEST成分と、CESTの影響を受けていない信号成分を表すベースライン成分とを含むＺスペクトルを作成するスペクトル作成処理と、
Ｚスペクトルを第１のスペクトルに変換するスペクトル変換処理であって、CESTの影響が現れる位相差において、前記第１のスペクトルのCEST成分の信号値とベースライン成分の信号値との比が、前記ＺスペクトルのCEST成分の信号値とベースライン成分の信号値との比よりも大きくなるように、前記Ｚスペクトルを前記第１のスペクトルに変換するスペクトル変換処理と、
前記第１のスペクトルに含まれるCEST成分の情報を求める処理と、
をコンピュータに実行させるためのプログラムである。
【発明の効果】
【００１１】
Ｚスペクトルは、CEST成分とベースライン成分との比が大きい第１のスペクトルに変換される。したがって、第１のスペクトルを用いることによりCESTの情報を取得することが容易となる。
【図面の簡単な説明】
【００１２】
【図１】本発明の第１の形態の磁気共鳴装置の概略図である。
【図２】第１の形態においてプロセッサ９が実現する手段の説明図である。
【図３】第１の形態で実行されるスキャンの説明図である。
【図４】第１の形態におけるシーケンスＳＥｋを具体的に示す図である。
【図５】Ｚスペクトルの説明図である。
【図６】ＺスペクトルとCPEスペクトルとの違いを説明するための図である。
【図７】CEST画像を作成するためのフローを示す図である。
【図８】Ｚスペクトルを概略的に示す図である。
【図９】CPEスペクトルＦCPE（Δωａ）を概略的に示す図である。
【図１０】フィッティングにより計算された各係数の値を示す図である。
【図１１】フィッティングにより計算された係数（ｂ０，ｂ１，ｂ２，Δω０）の値を示す図である。
【図１２】計算された係数（ｃ０，ｃＭＴ）の値を示す図である。
【図１３】スペクトルＦCPE＿１（Δωａ）を概略的に示す図である。
【図１４】CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれているか否かを判断する方法の説明図である。
【図１５】フィッティングにより計算された各係数の値を示す図である。
【図１６】ｎ＝２において、フィッティングにより計算された係数（ｂ０，ｂ１，ｂ２，Δω０）の値を示す図である。
【図１７】計算された係数（ｃ０，ｃＭＴ）の値を示す図である。
【図１８】スペクトルＦCPE＿２（Δωａ）を概略的に示す図である。
【図１９】CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれているか否かを判断する方法の説明図である。
【図２０】差分スペクトルＤ（Δωａ）が閾値ＴＨ１を超えた例を示す図である。
【図２１】ｎ＝ｉにおいて計算された係数の値を示す図である。
【図２２】ｎ＝ｉにおいて、フィッティングにより計算された係数（ｂ０，ｂ１，ｂ２，Δω０）の値を示す図である。
【図２３】計算された係数（ｃ０，ｃＭＴ）の値を示す図である。
【図２４】スペクトルＦCPE＿ｉ（Δωａ）を概略的に示す図である。
【図２５】ステップＳＴ１５のフローを示す図である。
【図２６】第２の形態で使用されるシーケンスＳＥｋを具体的に示す図である。
【図２７】第３の形態においてプロセッサ９が実現する手段の説明図である。
【図２８】第３の形態で使用されるシーケンスＳＥｋを具体的に示す図である。
【発明を実施するための形態】
【００１３】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
（１）第１の形態
図１は、本発明の第１の形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１は、マグネット２、テーブル３、受信ＲＦコイル（以下、「受信コイル」と呼ぶ）４などを有している。
【００１４】
マグネット２は、被検体１３が収容される収容空間２１を有している。また、マグネット２は、静磁場を発生させるための超伝導コイル２２、勾配パルスを印加するための勾配コイル２３、およびＲＦパルスを印加するためのＲＦコイル２４などを有している。超伝導コイルの代わりに、永久磁石を用いてもよい。
【００１５】
テーブル３は、クレードル３ａを有している。クレードル３ａは、収容空間２１内に移動できるように構成されている。クレードル３ａによって、被検体１３は収容空間２１に搬送される。
【００１６】
受信コイル４は、被検体１３の腹部および胸部を含む部位に取り付けられている。受信コイル４は、被検体１３からの磁気共鳴信号を受信する。
【００１７】
ＭＲ装置１は、更に、送信器５、勾配磁場電源６、受信器７、コンピュータ８、操作部１１、および表示部１２などを有している。
【００１８】
送信器５はＲＦコイル２４に電流を供給し、勾配磁場電源６は勾配コイル２３に電流を供給する。受信器７は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。尚、マグネット２、受信コイル４、送信器５、勾配磁場電源６、および受信器７を合わせたものがスキャン手段に相当する。
【００１９】
コンピュータ８は、表示部１２に必要な情報を伝送したり、画像を再構成するなど、ＭＲ装置１の各種の動作を実現するように、ＭＲ装置１の各部の動作を制御する。コンピュータ８は、プロセッサ９および記憶部１０などを有している。
【００２０】
記憶部１０には、プロセッサ９により実行されるプログラムなどが記憶されている。尚、記憶部１０は、コンピュータで読取り可能な非一過性の記録媒体であってもよい。非一過性の記録媒体としては、例えば、ＣＤ−ＲＯＭを用いることができる。プロセッサ９は、記憶部１０に記憶されているプログラムを読み出し、プログラムに記述されている処理を実行するための手段を実現する。図２は、プロセッサ９が実現する手段の説明図である。プロセッサ９は、記憶部１０に記憶されているプログラムを読み出すことにより、画像作成手段９０〜カウント手段１０３などを実現する。
【００２１】
画像作成手段９０は、後述するシーケンスＳＥ１〜ＳＥ１６（図１１参照）により得られたデータに基づいて画像を作成する。
【００２２】
Ｚスペクトル作成手段９１は、画像作成手段９０により得られた画像に基づいて、Ｚスペクトルを作成する。
【００２３】
スペクトル変換手段９２は、Ｚスペクトルを、後述するCPEスペクトルＦCPE（Δωａ）（例えば、図９参照）に変換する。CPEスペクトルＦCPE（Δωａ）は第１のスペクトルに相当する。
【００２４】
検出手段９３は、CPEスペクトルに基づいて、CESTの影響を受けた信号成分が現れるオフセット周波数を検出する。
【００２５】
設定手段９４は、CEST項の数を表すｎの値を設定する。
第１のフィッティング手段９５は、後述する式（１６）を用いてフィッティングを行い、CEST項に含まれる係数の値を計算する。CEST項については後述する。
【００２６】
CRZスペクトル作成手段９６は、ＺスペクトルからCEST成分が除去されたスペクトル（後述するCRZスペクトル）を作成する。CRZスペクトルは第２のスペクトルに相当する。
【００２７】
第２のフィッティング手段９７は、後述する式（２０）を用いてフィッティングを行い、式（２０）に含まれる係数（ｂ０，ｂ１，ｂ２，Δω０）の値を計算する。
【００２８】
係数値計算手段９８は、第２のフィッティング手段９７により計算された係数の値に基づいて、ベースライン項に含まれる係数の値を計算する。ベースライン項については後述する。
【００２９】
スペクトル計算手段９９は、後述するスペクトルＦCPE＿ｎ（Δωａ）（式２１参照）を計算する。スペクトルＦCPE＿ｎ（Δωａ）は第３のスペクトルに相当する。
【００３０】
CEST判断手段１００は、スペクトル変換手段９２により得られたCPEスペクトルＦCPE（Δωａ）に、CESTの影響を受けた他の信号成分が含まれているか否かを判断する。
【００３１】
スペクトル推定手段１０１は、ベースライン項とCEST項との和で表されるＺスペクトル（理想的なＺスペクトル）を推定する。
【００３２】
スペクトル比較手段１０２は、スペクトル推定手段１０１により推定された理想的なＺスペクトルと、Ｚスペクトル作成手段９１により作成されたＺスペクトルとを比較し、理想的なＺスペクトルによってＺスペクトルが再現されているか否かを判断する。
【００３３】
カウント手段１０３は、検出手段９３により検出されたオフセット周波数の総数をカウントする。
【００３４】
プロセッサ９は、画像作成手段９０〜カウント手段１０３を構成する一例であり、記憶部１０に記憶されたプログラムを実行することによりこれらの手段として機能する。尚、第１のフィッティング手段９５がCEST成分の情報を求める手段に相当し、CRZスペクトル作成手段９６、第２のフィッティング手段９７、および係数値計算手段９８を合わせたものがベースライン成分の情報を求める手段に相当する。
【００３５】
操作部１１は、オペレータにより操作され、種々の情報をコンピュータ８に入力する。表示部１２は種々の情報を表示する。
ＭＲ装置１は、上記のように構成されている。
【００３６】
図３は第１の形態で実行されるスキャンの説明図である。
スキャンＳＣは、CESTイメージング法を用いてスライスＳＬの画像を取得するためのスキャンである。スキャンＳＣでは、スライスＳＬの画像Ｄｋを取得するためのシーケンスＳＥｋ（ｋ＝１〜ｒ）が実行される。第１の形態では、シーケンスＳＥｋはｒ回実行されるので、スキャンＳＣを実行することにより、ｒ個の画像Ｄ１〜Ｄｒを取得することができる。
【００３７】
図４は、第１の形態におけるシーケンスＳＥｋを具体的に示す図である。
ｋ回目のシーケンスＳＥｋは、連続波のＲＦパルスＣＷと、横磁化を消失させるためのキラー勾配パルスＧｃと、シングルショット法によりデータを収集するためのデータ収集セグメントＤＡＱとを含んでいる。ＲＦパルスＣＷの周波数ｆ[Hz]は、ｆ＝ｆｋに設定されている。連続波のＲＦパルスＣＷを印加した後、キラー勾配パルスＧｃが印加され、キラー勾配パルスＧｃを印加した後にデータ収集セグメントＤＡＱが実行される。
【００３８】
ｋ回目のシーケンスＳＥｋは、上記のように構成されている。シーケンスＳＥ１、ＳＥ２、・・・ＳＥｒのＲＦパルスＣＷの周波数を、それぞれｆ１、ｆ２、・・・ｆｒで表すと、これらの周波数ｆ１、ｆ２、・・・ｆｒは、互いに異なる値に設定されている。
【００３９】
第１の形態では、シーケンスＳＥ１〜ＳＥｒを実行することにより、画像Ｄ１〜Ｄｒを取得し、これらの画像Ｄ１〜Ｄｒに基づいて、Ｚスペクトルを作成する。
【００４０】
図５は、Ｚスペクトルの説明図である。
図５（ａ１）は、Ｚスペクトルの波形を概略的に示す図である。Ｚスペクトルの横軸は、水の共鳴周波数からのずれを表すオフセット周波数Δωａである。Δωａは、Δωａ＝２π（ｆ−ｆｗ）［rad/sec］で計算される。ここで、ｆｗは、水の共鳴周波数である。
【００４１】
図５（ａ１）に示すように、Ｚスペクトルでは、或るオフセット周波数Δωｃにおいて、CESTによる信号減衰が現れる。したがって、Ｚスペクトルから、CESTの影響を受けた信号成分を分離することによって、CESTを定量的に評価することが可能となる。
【００４２】
図５（ａ２）は、Ｚスペクトルを、CESTの影響を受けた信号成分（以下、「CEST成分」と呼ぶ）Ｐ１と、CESTの影響を受けていない信号成分（以下、「ベースライン成分」と呼ぶ）Ｐ２とに分けて示した図である。尚、CEST成分Ｐ１は、実際には周波数Δωｃにおいて下向きのピークを有しているが、図５（ａ２）では、説明の便宜上、CEST成分Ｐ１のピークを上向きに反転させて示してある。
【００４３】
Ｚスペクトルから、CEST成分Ｐ１を分離することにより、CESTを定量的に評価することが可能となる。しかし、Ｚスペクトルは、CESTの影響を受けた信号成分（CEST成分）Ｐ１の他に、CESTの影響を受けていない信号成分（ベースライン成分）Ｐ２を含んでいる。一般的に、CEST成分Ｐ１およびベースライン成分Ｐ２は、ローレンツ関数で近似できる。しかし、ベースライン成分Ｐ２は、CEST成分Ｐ１よりも大きなピークを持つので、オフセット周波数Δωｃの近傍において、CEST成分Ｐ１とベースライン成分Ｐ２との比Ｒ（＝Ｐ１／Ｐ２）は小さい値になる。したがって、CEST成分Ｐ１のピークがベースライン成分Ｐ２に埋もれてしまい、ＺスペクトルからCEST成分Ｐ１を分離することが難しい場合がある。そこで、第１の形態では、Ｚスペクトルを、CEST成分の抽出に適したスペクトルに変換する。以下に、Ｚスペクトルを、CEST成分の抽出に適したスペクトルに変換する方法について説明する。
【００４４】
Ｚスペクトルは以下の式（１）で表すことができる。
【数１】
0006348449.tif 000002
【００４５】
ここで、Δωａ：水の共鳴周波数からのずれを表すオフセット周波数
Ｍｚａ：シーケンスＳＥｋのデータ収集セグメントＤＡＱ（図４参照）の直前における縦磁化の大きさ
Ｍ０：ＲＦパルスＷＣおよびキラー勾配パルスＧｃを印加せずにデータ収集セグメントＤＡＱを実行する場合において、データ収集セグメントＤＡＱの直前における磁化の大きさ、である。
尚、文字の添え字「a」は、自由水（free water）に起因することを表している。
【００４６】
また、Ｚスペクトルは、Zaiss等により以下の式（２）で近似できることが示されている（Zaiss M, et al. NMR Biomed 2013;26:507-18.）。
【数２】
0006348449.tif 000003
【００４７】
式（２）のＲ１ρは、ＲＦパルス印加中のＴ１回復の時定数であり、Trott等により以下の式（３）で近似できることが示されている（Trott O, et al. J Magn Reson 2002;154:157-60）。
【数３】
0006348449.tif 000004
【００４８】
式（３）のcos２θ、sin２θ、およびＲｅｘは以下の式で表される。
【数４】
0006348449.tif 000005
【００４９】
ω１：ＲＦパルスＷＣの送信磁場強度（Ｂ１強度）から求められる周波数[radian/sec]
Δωｃ：CESTの影響を受けた信号成分が現れるオフセット周波数[radian/sec]
ｋａ：自由水プールからCESTプールへの磁化移動（Magnetization Transfer）の時定数[Hz]
ｋ：CESTプールから自由水プールへの磁化移動の時定数[Hz]
【００５０】
ここで、Δωａ２を表す以下の関数Ｆ（Δωａ）について考える。
【数５】
0006348449.tif 000006
【００５１】
式（４）のＦ（Δωａ）は偶関数であり、Δωａ＝０の場合、Ｆ（Δωａ）＝０となる。式（４）を、式（３ａ）、（３ｂ）、および（３ｃ）に代入すると、以下の式が得られる。
【数６】
0006348449.tif 000007
【００５２】
次に、式（１）で表されるスペクトルＺの逆数のスペクトル１／Ｚについて考える。式（１）、（２）、（３）、および（５ｃ）から、１／Ｚは、以下の式で表すことができる。
【数７】
0006348449.tif 000008
【００５３】
式（６）を整理すると、以下の式（７）が得られる。
【数８】
0006348449.tif 000009
【００５４】
ここで、Δωａがω１に比べて十分に小さい場合、即ち、以下の関係が成り立つ場合について考える。
【数９】
0006348449.tif 000010
【００５５】
式（８）が成り立つ場合、式（５ａ）は、以下の式に近似できる。
【数１０】
0006348449.tif 000011
【００５６】
したがって、式（９ａ）から、cosθは、以下の式に近似できる。
【数１１】
0006348449.tif 000012
【００５７】
また、式（８）が成り立つ場合、式（５ｂ）は、以下の式に近似できる。
【数１２】
0006348449.tif 000013
【００５８】
式（９ａ−１）および（９ｂ）を用いて式（７）を整理すると以下の近似式が得られる。
【数１３】
0006348449.tif 000014
【００５９】
Δωａ≫ω１が成り立つ場合（式８参照）、式（１０）の近似式を得ることができる。右辺第１項Ｒ２ａω１２／Ｒ１ａは、自由水プールにおける緩和時間を表す項である。また、右辺第２項は、CEST成分（CESTの影響を受けた信号成分）を表す項であり、以下では、この項を、「CEST項」と呼ぶことにする。CEST項は、強度がａ１、幅がａ２、位置がΔωｃのローレンツ関数で表されている。したがって、CEST成分を、ローレンツ関数で表されるピークとして抽出できることがわかる。そこで、第１の形態では、式（１０）で表されるスペクトルをCPE（CEST Peak Extraction）スペクトルと呼ぶことにする。
【００６０】
CPEスペクトルを、ＦCPE（Δωａ）とすると、ＦCPE（Δωａ）は以下の式で表すことができる。
【数１４】
0006348449.tif 000015
【００６１】
また、式（１０）および（１１）から、CPEスペクトルＦCPE（Δωａ）の近似式は、以下の式で表すことができる。
【数１５】
0006348449.tif 000016
【００６２】
尚、上記の説明で使用された式（３）は、Trott等のモデルを仮定している。Trott等のモデルは、２つのプール（例えば、自由水およびＮＨ）の間で生じるCESTを考慮したTwo Poolモデルである。しかし、実際の生体組織では、結合水と自由水との間で生じる磁化移動（ＭＴ：Magnetization Transfer）も考慮する必要がある。ここで、結合水と自由水との間で生じるＭＴ（磁化移動）の影響を表すスペクトルをＺＭＴとし、スペクトルＺＭＴは、定数からローレンツ関数を引いた式で表現できると仮定する。この場合、スペクトルＺＭＴは以下の式で表すことができる。
【数１６】
0006348449.tif 000017
【００６３】
式（１３）において、Δω０は、周波数の誤差を表している。式（１３）の右辺第１項（ｂ０）は定数項であり、右辺第２項はローレンツ関数の項である。ここで、式（１１）の右辺のＺをＺＭＴで置き換え、更に、式（１３）のΔω０が十分に小さくΔω０＝０で表せると仮定する。この場合、ＺをＺＭＴで置き換えた後の式（１１）と、Δω０＝０と仮定された後の式（１３）とを用いて、以下の式が得られる。
【数１７】
0006348449.tif 000018
【００６４】
式（１４）の右辺第３項は十分に小さいので無視できる。したがって、式（１４）は、以下の式で近似できる。
【数１８】
0006348449.tif 000019
【００６５】
式（１５）の右辺第２項は、自由水と結合水との間で生じるMTの影響を受けた信号成分を表す項（以下、「ＭＴ項」と呼ぶ）である。
【００６６】
また、Two poolモデルでは、CESTピークは１つしか考えていないが、CESTピークが複数現れる場合もある。そこで、ｎ個のCEST項を考え、各CEST項とＭＴ項との間に一次結合が成り立つと仮定すると、式（１２）および式（１５）に基づいて、CPEスペクトルＦCPE（Δωａ）は、以下の式に近似することができる。
【数１９】
0006348449.tif 000020
【００６７】
ただし、Δωａ≫ω１である。式（１６）の右辺第１項および右辺第２項の和は、ベースライン成分（CESTの影響を受けていない信号成分）を表す項である。以下では、この項を、ベースライン項と呼ぶことにする。また、式（１６）の右辺第３項のＦＬ，ｎ（Δωａ）は、第ｎのCEST項を表している。したがって、式（１６）から、CPEスペクトルＦCPE（Δωａ）は、ベースライン項とｎ個のCEST項とを含む関数で近似できることが分かる。ベースライン項とｎ個のCEST項とを含む関数は、第１の関数に相当する。
【００６８】
式（１６）のベースライン項は、定数項ｃ０と、ＭＴ項ｃＭＴＦ（Δωａ）の和で表されている。ＭＴ項のＦ（Δωａ）は、ローレンツ関数ではなく、式（４）で定義される偶関数である。したがって、CPEスペクトルＦCPE（Δωａ）に含まれるベースライン成分（CESTの影響を受けていない信号成分）は、偶関数で近似できることが分かる。
【００６９】
図６は、ＺスペクトルとCPEスペクトルとの違いを説明するための図である。
図６（ａ１）はＺスペクトルの波形の概略図、図６（ａ２）はＺスペクトルをCEST成分Ｐ１とベースライン成分Ｐ２とに分けて示した図である。
【００７０】
図５を参照しながら説明したように、Ｚスペクトルのベースライン成分Ｐ２は、大きなピークを持つローレンツ関数で近似される。したがって、周波数Δωｃの近傍において、CEST成分Ｐ１とベースライン成分Ｐ２との比Ｒ（＝Ｐ１／Ｐ２）は小さい値になるので、CEST成分Ｐ１のピークがベースライン成分Ｐ２に埋もれてしまい、ＺスペクトルからCESTの影響を受けた信号成分Ｐ１を分離することが難しい場合がある。
【００７１】
一方、図６（ｂ１）はCPEスペクトルの波形の概略図、図６（ｂ２）はCPEスペクトルをCEST成分Ｑ１とベースライン成分Ｑ２とに分けて示した図である。
式（１６）の説明で言及したように、CPEスペクトルのベースライン成分は、偶関数で近似することができる。したがって、CPEスペクトルのベースライン成分Ｑ２はローレンツ関数による大きなピークを持たないので、周波数Δωｃの近傍において、CEST成分Ｑ１とべースライン成分Ｑ２との比Ｒ（＝Ｑ１／Ｑ２）を大きくすることができる。このため、CPEスペクトルは、Zスペクトルよりも、大きな比Ｒが得られるので、CPEスペクトルから、CEST成分Ｑ１を分離し易くすることができる。
【００７２】
尚、図６では、説明を簡単にするため、Ｚスペクトルには、CEST成分が一つしか含まれていない例が示されているが、Ｚスペクトルに複数のCEST成分が含まれている場合もある。しかし、Ｚスペクトルに複数のCEST成分が含まれている場合であっても、ＺスペクトルをCPEスペクトルに変換することによりベースライン成分の影響を軽減することができる。したがって、Ｚスペクトルに複数のCEST成分が含まれている場合であっても、ＺスペクトルをCPEスペクトルに変換することにより、CPEスペクトルから複数のCEST成分の各々を分離し易くすることができる。
【００７３】
第１の形態では、ＺスペクトルをCPEスペクトルに変換し、CPEスペクトルを用いてCEST画像を作成している。以下に、CPEスペクトルを用いてCEST画像を作成する方法について具体的に説明する。
【００７４】
図７は、CEST画像を作成するためのフローを示す図である。
ステップＳＴ１では、スキャンＳＣ（図４参照）が実行される。スキャンＳＣでは、シーケンスＳＥ１〜ＳＥｒが実行される。画像作成手段９０（図２参照）は、シーケンスＳＥ１〜ＳＥｒにより得られたデータに基づいて、スライスＳＬの画像Ｄ１〜Ｄｒ（図３参照）を作成する。シーケンスＳＥ１〜ＳＥｒのＲＦパルスＣＷの周波数は互いに異なる値に設定されているので、シーケンスＳＥ１〜ＳＥｒを実行することにより、ＲＦパルスの周波数をｒ通りに変化させたときの画像Ｄ１〜Ｄｒを得ることができる。シーケンスＳＥ１〜ＳＥｒを実行した後、ステップＳＴ２に進む。
【００７５】
ステップＳＴ２では、Ｚスペクトル作成手段９１（図２参照）がＺスペクトルを作成する。図８に、Ｚスペクトルを概略的に示す。Ｚスペクトル作成手段９１は、画像Ｄ１〜Ｄｒの同じ位置のピクセルを抽出し、水の共鳴周波数からの周波数のずれを表すオフセット周波数Δωａと信号値との関係を表すＺスペクトルを作成する。図８では、画像Ｄ１〜Ｄｒの同じ位置のピクセルｇ１におけるＺスペクトルが示されているが、他のピクセルにおけるＺスペクトルも作成する。
Ｚスペクトルを作成した後、ステップＳＴ３に進む。
【００７６】
ステップＳＴ３では、スペクトル変換手段９２（図２参照）が、式（１１）を用いて、ＺスペクトルをCPEスペクトルＦCPE（Δωａ）に変換する。図９に、CPEスペクトルＦCPE（Δωａ）を概略的に示す。ここでは、CPEスペクトルＦCPE（Δωａ）は、オフセット周波数Δωｃ，１およびΔωｃ，２においてCESTの影響を受けた信号成分（CEST成分）を含んでいるとする。尚、オフセット周波数の単位は［rad/sec］であるが、オフセット周波数の単位は、［rad/sec］から［ppm］に変換することができる。ここでは、オフセット周波数の単位は［ppm］に変換されているとする。ただし、説明の便宜上、オフセット周波数の単位を［rad/sec］から［ppm］に変換した後も、オフセット周波数は符号「Δωａ」で表してある。ＺスペクトルをCPEスペクトルＦCPE（Δωａ）に変換した後、ステップＳＴ４に進む。
【００７７】
ステップＳＴ４では、検出手段９３（図２参照）が、CPEスペクトルＦCPE（Δωａ）の中から、CESTの影響を受けた信号成分（CEST成分）が現れるオフセット周波数を検出する。以下に、CEST成分が現れるオフセット周波数の検出方法について説明する。
【００７８】
式（１６）のベースライン項に含まれるＦ（Δωａ）は、Δωａの二次関数で表されるので（式（４）参照）、CPEスペクトルＦCPE（Δωａ）のベースライン成分は二次関数で近似できることが分かる。したがって、CPEスペクトルＦCPE（Δωａ）と二次関数とを比較し、CPEスペクトルＦCPE（Δωａ）と二次関数とのずれが大きくなるときのオフセット周波数を求めることにより、CEST成分が現れるオフセット周波数を検出することができる。ここでは、CPEスペクトルＦCPE（Δωａ）は、オフセット周波数Δωｃ，１において、二次関数から最もずれているとする。したがって、検出手段９３は、CEST成分が現れるオフセット周波数として、オフセット周波数Δωｃ，１を検出する。ここでは、Δωｃ，１の検出値はΔωａ１であるとする。したがって、Δωｃ，１の検出値は、以下の式で表される。
【数２０】
0006348449.tif 000021
【００７９】
オフセット周波数Δωｃ，１＝Δωａ１を検出した後、ステップＳＴ５に進む。
ステップＳＴ５では、設定手段９４（図２参照）が、CPEスペクトルＦCPE（Δωａ）の近似式（１６）に含まれるCEST項の数を表すｎを初期値１に設定する。ｎ＝１に設定された場合、式（１６）は、以下の式で表される。
【数２１】
0006348449.tif 000022
【００８０】
式（１７）のベースライン項の係数（ｃ０，ｃＭＴ）、および式（１７ｃ）のCEST項ＦＬ，１（Δωａ）の係数（ａ１,１，ａ２,１，Δωｃ，１）は、未知の係数である。ｎ＝１に設定した後、ステップＳＴ６に進む。
【００８１】
ステップＳＴ６では、第１のフィッティング手段９５（図２参照）が、式（１１）により求められたCPEスペクトルＦCPE（Δωａ）と式（１７）との誤差が最小になるようにフィッティングを行い、誤差が最小になるときの式（１７）のCEST項ＦＬ，１（Δωａ）の係数（ａ１,１，ａ２,１，Δωｃ，１）の値と、ベースライン項の係数（ｃ０，ｃＭＴ）の値を計算する。フィッティングを行う場合、第１のフィッティング手段９５は、先ず、係数（ａ１,１，ａ２,１，Δωｃ，１）の初期値と、係数（ｃ０，ｃＭＴ）の初期値を設定する。例えば、係数Δωｃ，１の初期値は、ステップＳＴ４で検出したオフセット周波数Δωｃ，１の値、即ち、Δωｃ，１＝Δωａ１（式１６ｄ参照）に設定される。また、他の係数ａ１,１、ａ２,１、ｃ０、ｃＭＴの初期値は、CPEスペクトルＦCPE（Δωａ）の特徴量（CPEスペクトルの最大値、最小値など）を用いて計算することができる。係数の初期値を設定した後、第１のフィッティング手段９５は、初期値を基準にして各係数の値を変更し、式（１１）により求められたCPEスペクトルＦCPE（Δωａ）と式（１７）との誤差が最小になるときの係数（ａ１,１，ａ２,１，Δωｃ，１）の値、および係数（ｃ０，ｃＭＴ）の値を計算する。図１０に、フィッティングにより計算された各係数の値を示す。図１０では、（ｃ０，ｃＭＴ）＝（ｃ０(1)，ｃＭＴ(1)）、（ａ１,１，ａ２,１，Δωｃ，１）＝（ａ１,１(1)，ａ２,１(1)，Δωｃ，１(1)）で示されている。
【００８２】
尚、フィッティングにより、CEST項ＦＬ，１（Δωａ）の係数の値の他に、ベースライン項の係数（ｃ０，ｃＭＴ）の値も計算することができる。ただし、CPEスペクトルＦCPE（Δωａ）は、Ｚスペクトルよりも、ベースライン成分が抑制されている（図６参照）。したがって、CPEスペクトルＦCPE（Δωａ）を近似式（１７）でフィッティングすることによりベースライン項の係数（ｃ０，ｃＭＴ）を求めた場合、係数（ｃ０，ｃＭＴ）の推定誤差が大きくなる可能性がある。そこで、第１の形態では、係数（ｃ０，ｃＭＴ）の推定誤差が小さくなるように、係数（ｃ０，ｃＭＴ）を計算し直す。係数（ｃ０，ｃＭＴ）を計算し直すために、ステップＳＴ７に進む。
【００８３】
ステップＳＴ７は、２つのステップＳＴ７１およびＳＴ７２を有している。以下、各ステップＳＴ７１およびＳＴ７２について説明する。
【００８４】
ステップＳＴ７１では、CRZスペクトル作成手段９６（図２参照）が、ＺスペクトルからCEST成分が除去されたスペクトルを作成する。以下では、ＺスペクトルからCEST成分が除去されたスペクトルをCRZスペクトル（CEST Removed Z-spectrum）と呼ぶことにする。CRZスペクトルを、「ＺＣＲＺ」で表すと、CRZスペクトルＺＣＲＺは、Ｚスペクトルを用いて、以下の式で表すことができる。
【数２２】
0006348449.tif 000023
【００８５】
尚、δ（Δωａ）は、Δωａ＝０において式（１８）の右辺第２項の分母がゼロにならないようにするために導入された関数である。ここでは、ｎ＝１であるので、式（１８）は、以下の式で表される。
【数２３】
0006348449.tif 000024
【００８６】
式（１９）のＺは、ステップＳＴ２で求められている。また、式（１９ａ）のＦＬ，１（Δωａ）の係数（ａ１,１，ａ２,１，Δωｃ，１）は、図１０に示すように、ステップＳＴ６で（ａ１,１，ａ２,１，Δωｃ，１）＝（ａ１,１(1)，ａ２,１(1)，Δωｃ，１(1)）と計算されている。したがって、式（１９）から、CESTの影響を受けた信号成分（CEST成分）が除去されたCRZスペクトルＺＣＲＺを作成することができる。CRZスペクトルＺＣＲＺを作成した後、ステップＳＴ７２に進む。
【００８７】
ステップＳＴ７２では、ステップＳＴ７１で作成したCRZスペクトルＺＣＲＺに基づいて、式（１７）のベースライン項の係数（ｃ０，ｃＭＴ）の値を計算する。以下に、係数（ｃ０，ｃＭＴ）の値の求め方について説明する。
【００８８】
CRZスペクトルＺＣＲＺは、ＺスペクトルからCEST成分が除去されたスペクトルを表している。したがって、CRZスペクトルＺＣＲＺは、CESTの影響を受けた信号成分ではなく、自由水と結合水との間で生じるMTの影響を受けた信号成分を主に表していると考えることができる。自由水と結合水との間で生じるMTの影響を受けた信号成分を表すスペクトルは、式（１３）のスペクトルＺＭＴで表されている。したがって、CRZスペクトルＺＣＲＺは、スペクトルＺＭＴを用いて、以下の式で近似することができる。
【数２４】
0006348449.tif 000025
【００８９】
式（２０）から、CRZスペクトルＺＣＲＺは、定数項（ｂ０）とローレンツ関数の項との和で表される関数で近似できることが分かる。定数項（ｂ０）とローレンツ関数の項との和で表される関数は、第２の関数に相当する。
【００９０】
第２のフィッティング手段９７（図２参照）は、CRZスペクトルＺＣＲＺと式（２０）との誤差が最小になるようにフィッティングを行い、誤差が最小になるときの式（２０）の係数（ｂ０，ｂ１，ｂ２，Δω０）の値を計算する。フィッティングを行う場合、第２のフィッティング手段９７は、先ず、係数（ｂ０，ｂ１，ｂ２，Δω０）の初期値を計算する。係数（ｂ０，ｂ１，ｂ２，Δω０）の初期値は、例えば、ステップＳＴ６で計算したベースライン項の値（ｃ０，ｃＭＴ）＝（ｃ０(1)，ｃＭＴ(1)）に基づいて計算することができる。係数（ｂ０，ｂ１，ｂ２，Δω０）の初期値を計算した後、第２のフィッティング手段９７は、初期値を基準にして係数（ｂ０，ｂ１，ｂ２，Δω０）の値を変更し、CRZスペクトルＺＣＲＺと式（２０）との誤差が最小になるときの係数（ｂ０，ｂ１，ｂ２，Δω０）の値を計算する。図１１に、フィッティングにより計算された係数（ｂ０，ｂ１，ｂ２，Δω０）の値を示す。図１１では、計算された係数（ｂ０，ｂ１，ｂ２，Δω０）の値は、（ｂ０，ｂ１，ｂ２，Δω０）＝（ｂ０(1)，ｂ１(1)，ｂ２(1)，Δω０(1)）で示されている。
【００９１】
これらの係数の値を求めた後、係数値計算手段９８（図２参照）は、（ｂ０，ｂ１）＝（ｂ０(1)，ｂ１(1)）を式（１７ａ）に代入し、ｃ０を計算する。また、係数値計算手段９８は、ｂ０＝ｂ０(1)を式（１７ｂ）に代入し、ｃＭＴを計算する。したがって、式（１７）のベースライン項の係数（ｃ０，ｃＭＴ）の値を計算することができる。図１２に、計算された係数（ｃ０，ｃＭＴ）の値を示す。図１２では、（ｃ０，ｃＭＴ）＝（ｃ０(1)’，ｃＭＴ(1)’）で示されている。これらの値を求めた後、ステップＳＴ８に進む。
【００９２】
ステップＳＴ８では、スペクトル計算手段９９（図２参照）が、ベースライン項ｃ０＋ｃＭＴＦ（Δωａ）とCEST項ΣｎＦＬ，ｎ（Δωａ）との和で表されるスペクトルＦCPE＿ｎ（Δωａ）を計算する。このスペクトルＦCPE＿ｎは、以下の式で定義することができる。
【数２５】
0006348449.tif 000026
【００９３】
尚、ステップＳＴ５において、ｎ＝１に設定されているので、式（２１）は、以下の式で表される。
【数２６】
0006348449.tif 000027
【００９４】
式（２２）のベースライン項の係数（ｃ０，ｃＭＴ）は、ステップＳＴ７２において（ｃ０，ｃＭＴ）＝（ｃ０(1)’，ｃＭＴ(1)’）と計算されている（図１２参照）。また、式（２２ｃ）のCEST項ＦＬ，１（Δωａ）の係数（ａ１,１，ａ２,１，Δωｃ，１）は、ステップＳＴ６において（ａ１,１，ａ２,１，Δωｃ，１）＝（ａ１,１(1)，ａ２,１(1)，Δωｃ，１(1)）と求められている（図１２参照）。したがって、式（２２）および（２２ｃ）に、これらの係数の値を代入することにより、スペクトルＦCPE＿１（Δωａ）を計算することができる。図１３に、スペクトルＦCPE＿１（Δωａ）を概略的に示す。尚、図１３では、比較のため、スペクトルＦCPE＿１（Δωａ）の他に、CPEスペクトルＦCPE（Δωａ）も示してある。スペクトルＦCPE＿１（Δωａ）は、オフセット周波数Δωｃ，１の付近において、CPEスペクトルＦCPE（Δωａ）に十分に近い波形を有していることが分かる。スペクトルＦCPE＿１（Δωａ）を計算した後、ステップＳＴ９に進む。
【００９５】
ステップＳＴ９では、CEST判断手段１００（図２参照）が、CPEスペクトルＦCPE（Δωａ）に、CEST項ＦＬ，１（Δωａ）で表されるCEST成分（式２２ｃ参照）とは異なる他のCEST成分が含まれているか否かを判断する。
【００９６】
図１４は、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれているか否かを判断する方法の説明図である。
図１４の上側には、CPEスペクトルＦCPE（Δωａ）およびスペクトルＦCPE＿１（Δωａ）が示されており、図１４の下側には、CPEスペクトルＦCPE（Δωａ）とスペクトルＦCPE＿１（Δωａ）との差分スペクトルＤ（Δωａ）が示されている。
【００９７】
先ず、CEST判断手段１００は、CPEスペクトルＦCPE（Δωａ）から、スペクトルＦCPE＿１（Δωａ）を差分し、差分スペクトルＤ（Δωａ）を求める。
【００９８】
次に、CEST判断手段１００は、差分スペクトルＤ（Δωａ）に基づいて、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれているか否かを判断する。以下に、この判断方法について説明する。
【００９９】
スペクトルＦCPE＿１（Δωａ）を求めるための式（２２）は、オフセット周波数Δωｃ,１に対応したCEST項ＦＬ，１（Δωａ）を含んでいる。したがって、オフセット周波数Δωｃ,１の付近では、スペクトルＦCPE＿１（Δωａ）は、CPEスペクトルＦCPE（Δωａ）に十分に近い値を有している。このため、差分スペクトルＤ（Δωａ）の信号値は、オフセット周波数Δωｃ,１の付近ではゼロに近い値になる。
【０１００】
しかし、ＦCPE＿１（Δωａ）の式（２２）は、オフセット周波数Δωｃ,２に対応したCEST項を含んでいない。したがって、オフセット周波数Δωｃ,２の付近では、スペクトルＦCPE＿１（Δωａ）と、CPEスペクトルＦCPE（Δωａ）との間に、ある程度の信号値の差が生じる。このため、差分スペクトルＤ（Δωａ）には、オフセット周波数Δωｃ,２の付近に、CESTの影響を受けた信号成分を表すピークＰ２が現れる。
【０１０１】
したがって、差分スペクトルＤ（Δωａ）にピークＰ２が現れているかを判断することにより、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれているか否かを判断することができる。第１の形態では、差分スペクトルＤ（Δωａ）にピークＰ２が現れているかを判断するために、２つの閾値ＴＨ１およびＴＨ２が用いられる。CEST判断手段１００は、２つの閾値ＴＨ１およびＴＨ２と、差分スペクトルＤ（Δωａ）とを比較し、差分スペクトルＤ（Δωａ）が、閾値ＴＨ１又は閾値ＴＨ２を横切っているか否かを判断する。差分スペクトルＤ（Δωａ）が、閾値ＴＨ１又は閾値ＴＨ２を横切っている場合、CEST判断手段１００は、差分スペクトルＤ（Δωａ）にピークＰ２が現れていると判断する。一方、差分スペクトルＤ（Δωａ）が、閾値ＴＨ１又は閾値ＴＨ２を横切っていない場合、CEST判断手段１００は、差分スペクトルにピークＰ２が現れていないと判断する。
【０１０２】
図１４を参照すると、差分スペクトルＤ（Δωａ）は、オフセット周波数Δωｃ,２の付近において、閾値ＴＨ１を超えるピークＰ２を含んでいる。したがって、CEST判断手段１００は、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれていると判断する。CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれていると判断された場合、ステップＳＴ１０に進む。
【０１０３】
ステップＳＴ１０では、検出手段９３が、CPEスペクトルＦCPE（Δωａ）の他のCEST成分が現れるオフセット周波数Δωｃ，２を検出する。ここでは、Δωｃ，２の検出値はΔωａ２であるとする。したがって、Δωｃ，２の検出値は、以下の式で表される。
【数２７】
0006348449.tif 000028
【０１０４】
オフセット周波数Δωｃ，２＝Δωａ２を検出した後、ステップＳＴ１１に進む。
ステップＳＴ１１では、設定手段９４が、CEST項の数を表すｎをインクリメントする。したがって、ｎは、ｎ＝１からｎ＝２に設定される。ｎ＝２に設定された場合、CPEスペクトルの近似式（１６）は、以下の式で表される。
【数２８】
0006348449.tif 000029
【０１０５】
ｎ＝２に設定した後、ステップＳＴ６に戻る。
ステップＳＴ６では、第１のフィッティング手段９５が、フィッティングにより、式（２３ｃ＿２）のCEST項ＦＬ，２（Δωａ）の係数（ａ１,２，ａ２,２，Δωｃ，２）の値を計算する。以下に、係数（ａ１,２，ａ２,２，Δωｃ，２）の求め方について説明する。
【０１０６】
第１の形態では、既に、CEST項ＦＬ，１（Δωａ）の係数（ａ１,１，ａ２,１，Δωｃ，１）は、（ａ１,１(1)，ａ２,１(1)，Δωｃ，１(1)）と計算されており、ベースライン項の係数（ｃ０，ｃＭＴ）は、（ｃ０(1)’，ｃＭＴ(1)’）と計算されている（図１２参照）。したがって、これらの値を式（２３）および式（２３ｃ＿１）に代入すると、式（２３ｃ＿２）で表されるCEST項ＦＬ，２（Δωａ）の３つの係数（ａ１,２，ａ２,２，Δωｃ，２）のみが、未知の係数となる。この場合、CPEスペクトルＦCPE（Δωａ）と、３つの未知の係数（ａ１,２，ａ２,２，Δωｃ，２）を含む近似式（２３）との誤差が最小になるように、CPEスペクトルＦCPE（Δωａ）をフィッティングすることにより、３つの係数（ａ１,２，ａ２,２，Δωｃ，２）の値を求めることができる。
【０１０７】
しかし、ベースライン項の係数（ｃ０，ｃＭＴ）＝（ｃ０(1)’，ｃＭＴ(1)’）は、CEST項が１項（ＦＬ，１（Δωａ））しか含まれていない近似式（１７）に基づいて求められた値である。一方、近似式（２３）は、CEST項ＦＬ，１（Δωａ）の他に、新たなCEST項ＦＬ，２（Δωａ）が追加されている。したがって、ベースライン項の係数（ｃ０，ｃＭＴ）を（ｃ０(1)’，ｃＭＴ(1)’）に固定してフィッティングを行うと、CEST項ＦＬ，２（Δωａ）の係数（ａ１,２，ａ２,２，Δωｃ，２）の推定誤差が大きくなる恐れがある。そこで、第１の形態では、CEST項ＦＬ，２（Δωａ）の係数（ａ１,２，ａ２,２，Δωｃ，２）の推定誤差を小さくするために、CEST項ＦＬ，２（Δωａ）の係数（ａ１,２，ａ２,２，Δωｃ，２）だけでなく、ベースライン項の係数（ｃ０、ｃＭＴ）も未知の係数として、フィッティングを行う。したがって、５つの係数が未知の係数となる。第１のフィッティング手段９５は、５つの未知の係数を含む近似式（２３）を用いて、CPEスペクトルＦCPE（Δωａ）のフィッティングを行う。フィッティングを行う場合、第１のフィッティング手段９５は、先ず、CEST項ＦＬ，２（Δωａ）の係数（ａ１,２，ａ２,２，Δωｃ，２）の初期値と、ベースライン項の係数（ｃ０，ｃＭＴ）の初期値とを設定する。係数Δωｃ，２の初期値は、ステップＳＴ１０で検出したオフセット周波数Δωｃ，２の値、即ち、Δωｃ，２＝Δωａ２（式２２ｄ参照）に設定される。また、係数ａ１,２、ａ２,２の初期値は、差分スペクトルＤ（Δωａ）のピークＰ２（図１４参照）の特徴量（ピークの高さ、半値幅など）に基づいて計算することができる。一方、ベースライン項の係数（ｃ０，ｃＭＴ）の初期値は、ｎ＝１のときに求めた値（ｃ０，ｃＭＴ）＝（ｃ０(1)’，ｃＭＴ(1)’）（図１２参照）に設定することができる。初期値を設定した後、第１のフィッティング手段９５は、初期値を基準にして係数の値を変更し、CPEスペクトルＦCPE（Δωａ）と式（２３）との誤差が最小になるときの係数（ａ１,２，ａ２,２，Δωｃ，２）および係数（ｃ０，ｃＭＴ）の値を計算する。図１５に、フィッティングにより計算された各係数の値を示す。図１５では、（ｃ０，ｃＭＴ）＝（ｃ０(2)，ｃＭＴ(2)）、（ａ１,２，ａ２,２，Δωｃ，２）＝（ａ１,２(2)，ａ２,２(2)，Δωｃ，２(2)）で示されている。
【０１０８】
尚、フィッティングにより、CEST項ＦＬ，２（Δωａ）の係数（ａ１,２，ａ２,２，Δωｃ，２）の他に、ベースライン項の係数（ｃ０，ｃＭＴ）も計算される。ただし、先に説明したように、ステップＳＴ６で計算されたベースライン項の係数（ｃ０，ｃＭＴ）の値は、推定誤差が大きい可能性がある。そこで、ベースライン項の係数（ｃ０，ｃＭＴ）を計算し直すために、ステップＳＴ７に進む。
【０１０９】
ステップＳＴ７では、２つのステップＳＴ７１およびＳＴ７２が順に実行される。
ステップＳＴ７１では、CRZスペクトル作成手段９６が、式（１８）を用いて、ＺスペクトルからCEST成分が除去されたCRZスペクトルＺＣＲＺを作成する。ただし、ｎ＝２に設定されているので、式（１８）は、以下の式で表される。
【数２９】
0006348449.tif 000030
【０１１０】
CRZスペクトル作成手段９６は、式（２４）を用いて、ＺスペクトルからCEST成分が除去されたCRZスペクトルＺＣＲＺを作成する。CRZスペクトルＺＣＲＺを作成した後、ステップＳＴ７２に進む。
【０１１１】
ステップＳＴ７２では、第２のフィッティング手段９７が、式（２４）を用いて作成されたCRZスペクトルＺＣＲＺと式（２０）との誤差が最小になるようにフィッティングを行い、誤差が最小になるときの式（２０）の係数（ｂ０，ｂ１，ｂ２，Δω０）の値を計算する。フィッティングを行う場合、第２のフィッティング手段９７は、先ず、係数（ｂ０，ｂ１，ｂ２，Δω０）の初期値を設定する。ここでは、第２のフィッティング手段９７は、ｎ＝１のときに計算された係数（ｂ０，ｂ１，ｂ２，Δω０）＝（ｂ０(1)，ｂ１(1)，ｂ２(1)，Δω０(1)）（図１２参照）を、ｎ＝２における係数（ｂ０，ｂ１，ｂ２，Δω０）の初期値として設定する。ｎ＝２における係数（ｂ０，ｂ１，ｂ２，Δω０）の初期値を設定した後、第２のフィッティング手段９７は、初期値を基準にして係数（ｂ０，ｂ１，ｂ２，Δω０）の値を変更し、式（２４）を用いて作成されたCRZスペクトルＺＣＲＺと式（２０）との誤差が最小になるときの係数（ｂ０，ｂ１，ｂ２，Δω０）の値を計算する。図１６に、ｎ＝２において、フィッティングにより計算された係数（ｂ０，ｂ１，ｂ２，Δω０）の値を示す。図１６では、計算された係数（ｂ０，ｂ１，ｂ２，Δω０）の値は、（ｂ０，ｂ１，ｂ２，Δω０）＝（ｂ０(2)，ｂ１(2)，ｂ２(2)，Δω０(2)）で示されている。
【０１１２】
これらの係数の値を求めた後、係数値計算手段９８は、（ｂ０，ｂ１）＝（ｂ０(2)，ｂ１(2)）を式（２３ａ）に代入し、ｃ０を計算する。また、係数値計算手段９８は、ｂ０＝ｂ０(2)を式（２３ｂ）に代入し、ｃＭＴを計算する。したがって、式（２３）のベースライン項の係数（ｃ０，ｃＭＴ）の値を計算することができる。図１７に、計算された係数（ｃ０，ｃＭＴ）の値を示す。図１７は、（ｃ０，ｃＭＴ）＝（ｃ０(2)’，ｃＭＴ(2)’）で示されている。これらの値を求めた後、ステップＳＴ８に進む。
【０１１３】
ステップＳＴ８では、スペクトル計算手段９９が、式（２１）を用いてスペクトルＦCPE＿ｎ（Δωａ）を計算する。ただし、ｎ＝２に設定されているので、式（２１）は、以下の式で表される。
【数３０】
0006348449.tif 000031
【０１１４】
図１７に示すように、係数（ｃ０，ｃＭＴ）は、（ｃ０，ｃＭＴ）＝（ｃ０(2)’，ｃＭＴ(2)’）と計算されており、係数（ａ１,２，ａ２,２，Δωｃ，２）は、（ａ１,２，ａ２,２，Δωｃ，２）＝（ａ１,２(2)，ａ２,２(2)，Δωｃ，２(2)）と計算されている。また、図１２に示すように、係数（ａ１,１，ａ２,１，Δωｃ，１）は、（ａ１,１，ａ２,１，Δωｃ，１）＝（ａ１,１(1)，ａ２,１(1)，Δωｃ，１(1)）と計算されている。したがって、式（２５）、（２５ｃ＿１）、および（２５ｃ＿２）に、これらの係数の値を代入することにより、スペクトルＦCPE＿２（Δωａ）を計算することができる。図１８に、スペクトルＦCPE＿２（Δωａ）を概略的に示す。尚、図１８では、比較のため、スペクトルＦCPE＿２（Δωａ）の他に、CPEスペクトルＦCPE（Δωａ）も示してある。スペクトルＦCPE＿２（Δωａ）は、オフセット周波数Δωｃ，１およびΔωｃ，２の付近において、CPEスペクトルＦCPE（Δωａ）に十分に近い波形を有していることが分かる。スペクトルＦCPE＿２（Δωａ）を計算した後、ステップＳＴ９に進む。
【０１１５】
ステップＳＴ９では、CEST判断手段１００が、CPEスペクトルＦCPE（Δωａ）に、CEST項ＦＬ，１（Δωａ）およびＦＬ，２（Δωａ）で表されるCEST成分（式２３ｃ＿１および２３ｃ＿２参照）とは異なる他のCEST成分が含まれているか否かを判断する。
【０１１６】
図１９は、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれているか否かを判断する方法の説明図である。
CEST判断手段１００は、CPEスペクトルＦCPE（Δωａ）から、スペクトルＦCPE＿２（Δωａ）を差分し、差分スペクトルＤ（Δωａ）を求め、差分スペクトルＤ（Δωａ）と閾値ＴＨ１およびＴＨ２とを比較する。
【０１１７】
差分スペクトルＤ（Δωａ）は、閾値ＴＨ１およびＴＨ２を超えていないので、スペクトルＦCPE＿２（Δωａ）は、CPEスペクトルＦCPE（Δωａ）に十分に近い波形を有していると考えることができる。この場合、CPEスペクトルＦCPE（Δωａ）に含まれるCEST成分は、スペクトルＦCPE＿２（Δωａ）の式（２５）に含まれている２つのCEST項ＦＬ，１（Δωａ）およびＦＬ，２（Δωａ）で十分に表すことができたと考えられる。したがって、差分スペクトルＤ（Δωａ）が閾値ＴＨ１およびＴＨ２を超えていない場合、CEST判断手段１００は、CPEスペクトルＦCPE（Δωａ）に他のCEST成分は含まれていないと判断する。他のCEST成分が存在していないと判断された場合、ステップＳＴ１２に進む。
【０１１８】
尚、図１９では、差分スペクトルＤ（Δωａ）が閾値ＴＨ１およびＴＨ２を超えていない例が示されている。しかし、差分スペクトルＤ（Δωａ）が閾値ＴＨ１又はＴＨ２を超える場合もある。以下に、差分スペクトルＤ（Δωａ）が閾値ＴＨ１又はＴＨ２を超える場合について説明する。
【０１１９】
図２０は、差分スペクトルＤ（Δωａ）が閾値ＴＨ１を超えた例を示す図である。
図２０では、差分スペクトルＤ（Δωａ）には、オフセット周波数Δωｃ，３の付近に、閾値ＴＨ１を超えたピークＰ３が現れている。したがって、CEST判断手段１００は、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれていると判断する。CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれていると判断された場合、ステップＳＴ１０に進み、他のCEST成分が現れるオフセット周波数Δωｃ，３を検出する。オフセット周波数Δωｃ，３を検出した後、ステップＳＴ１１に進み、ｎがインクリメントされ、CPEスペクトルＦCPE（Δωａ）の近似式（２３）に新たなCEST項ＦＬ，３（Δωａ）が追加される。そして、ステップＳＴ６〜ＳＴ９が実行される。したがって、ステップＳＴ９において、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれていると判断されるたびに、CPEスペクトルＦCPE（Δωａ）の近似式に新たなCEST項が追加され、ステップＳＴ６〜ＳＴ９が実行される。例えば、ステップＳＴ９において、CPEスペクトルＦCPE（Δωａ）にｉ個目のCEST成分が含まれていると判断された場合について考えてみる。この場合、ステップＳＴ１１においてｎ＝ｉに設定されるので、CPEスペクトルＦCPE（Δωａ）の近似式は、以下の式で表される。
【数３１】
0006348449.tif 000032
【０１２０】
上記の近似式（２６）では、ＦＬ，ｉ（Δωａ）が、新たに追加されたCEST項を表している。CPEスペクトルが近似式（２６）で表される場合、式（２６ｃ＿１）〜式（２６ｃ＿ｉ−１）のCEST項の係数は、既に計算されている。したがって、係数（ｃ０，ｃＭＴ）および係数（ａ１,ｉ，ａ２,ｉ，Δωｃ，ｉ）の５つの係数が未知の係数となる。ステップＳＴ１１においてｎ＝ｉに設定されたら、ステップＳＴ６に戻る。ステップＳＴ６では、第１のフィッティング手段９５が、近似式（２６）を用いて５つの係数（ｃ０，ｃＭＴ）および（ａ１,ｉ，ａ２,ｉ，Δωｃ，ｉ）の値を計算する。図２１に、ｎ＝ｉにおいて計算された係数の値を示す。
【０１２１】
CEST項の係数の値を計算した後、ステップＳＴ７１に進む。ステップＳＴ７１では、CRZスペクトル作成手段９６が、式（１８）を用いて、ＺスペクトルからCEST成分が除去されたCRZスペクトルＺＣＲＺを作成する。尚、ｎ＝ｉであるので、式（１８）は、以下の式で表される。
【数３２】
0006348449.tif 000033
【０１２２】
CRZスペクトルＺＣＲＺを求めた後、ステップＳＴ７２に進む。
ステップＳＴ７２では、第２のフィッティング手段９７が、式（２７）により作成されたCRZスペクトルＺＣＲＺと式（２０）との誤差が最小になるようにフィッティングを行い、誤差が最小になるときの式（２０）の係数（ｂ０，ｂ１，ｂ２，Δω０）の値を計算する。フィッティングを行う場合、第２のフィッティング手段９７は、先ず、係数（ｂ０，ｂ１，ｂ２，Δω０）の初期値を設定する。ここでは、第２のフィッティング手段９７は、ｎ＝ｉ−１のときに計算された係数（ｂ０，ｂ１，ｂ２，Δω０）の値（図示せず）を、ｎ＝ｉにおける係数（ｂ０，ｂ１，ｂ２，Δω０）の初期値として設定する。ｎ＝ｉにおける係数（ｂ０，ｂ１，ｂ２，Δω０）の初期値を設定した後、第２のフィッティング手段９７は、初期値を基準にして係数（ｂ０，ｂ１，ｂ２，Δω０）の値を変更し、式（２７）を用いて作成されたCRZスペクトルＺＣＲＺと式（２０）との誤差が最小になるときの係数（ｂ０，ｂ１，ｂ２，Δω０）の値を計算する。図２２に、ｎ＝ｉにおいて、フィッティングにより計算された係数（ｂ０，ｂ１，ｂ２，Δω０）の値を示す。図２２では、（ｂ０，ｂ１，ｂ２，Δω０）＝（ｂ０(i)，ｂ１(i)，ｂ２(i)，Δω０(i)）で示されている。
【０１２３】
これらの係数の値を求めた後、係数値計算手段９８は、（ｂ０，ｂ１）＝（ｂ０(i)，ｂ１(i)）を式（２６ａ）に代入し、ｃ０を計算する。また、係数値計算手段９８は、ｂ０＝ｂ０(i)を式（２６ｂ）に代入し、ｃＭＴを計算する。したがって、式（２６）のベースライン項の係数（ｃ０，ｃＭＴ）の値を計算することができる。図２３に、計算された係数（ｃ０，ｃＭＴ）の値を示す。図２３では、（ｃ０，ｃＭＴ）＝（ｃ０(i)’，ｃＭＴ(i)’）で示されている。これらの値を計算した後、ステップＳＴ８に進む。
【０１２４】
ステップＳＴ８では、スペクトル計算手段９９が、式（２１）を用いてスペクトルＦCPE＿ｎ（Δωａ）を計算する。ただし、ｎ＝ｉであるので、式（２１）は、以下の式で表される。
【数３３】
0006348449.tif 000034
【０１２５】
式（２８）の各係数は既に求められているので、これらの係数の値を代入することにより、スペクトルＦCPE＿ｉ（Δωａ）を計算することができる。図２４に、スペクトルＦCPE＿ｉ（Δωａ）を概略的に示す。尚、図２４では、比較のため、スペクトルＦCPE＿ｉ（Δωａ）の他に、CPEスペクトルＦCPE（Δωａ）も示してある。スペクトルＦCPE＿ｉ（Δωａ）を計算した後、ステップＳＴ９に進む。
【０１２６】
ステップＳＴ９では、CEST判断手段１００が、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれているか否かを判断する。したがって、ステップＳＴ９において、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれていないと判断されるまでは、ステップＳＴ１０、ステップＳＴ１１、ステップＳＴ６、ステップＳＴ７、ステップＳＴ８、およびステップＳＴ９のループが繰り返し実行される。そして、ステップＳＴ９において、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれていないと判断されると、ステップＳＴ１２に進む。
【０１２７】
ステップＳＴ１２では、スペクトル推定手段１０１（図２参照）が、ベースライン項とCEST項との和で表されるＺスペクトル（以下、「理想的なＺスペクトル」と呼ぶ）Ｚidealを推定する。理想的なＺスペクトルＺidealは、式（１１）および式（１６）を用いて、以下の式で表すことができる。
【数３４】
0006348449.tif 000035
【０１２８】
スペクトル推定手段１０１は、式（２９）を用いて、理想的なＺスペクトルＺidealを推定する。理想的なＺスペクトルＺidealを求めた後、ステップＳＴ１３に進む。
【０１２９】
ステップＳＴ１３では、スペクトル比較手段１０２（図２参照）が、理想的なＺスペクトルＺidealとステップＳＴ２で作成されたＺスペクトルとを比較し、理想的なＺスペクトルＺidealによってＺスペクトルが再現されているか否かを判断する。この判断は、以下のように行う。
【０１３０】
先ず、スペクトル比較手段１０２は、理想的なＺスペクトルＺidealと、ステップＳＴ２で作成されたＺスペクトルとを比較し、オフセット周波数ごとに、理想的なＺスペクトルＺidealの信号値と、Ｚスペクトルの信号値との差を求める。次に、スペクトル比較手段１０２は、差の平方和が十分に小さいか否かを判断する。差の平方和が大きいか否かの判断は、例えば、差の平方和が大きいか小さいかを判断するための閾値を予め決定しておき、差の平方和と閾値とを比較することにより、行うことができる。スペクトル比較手段１０２は、差の平方和が閾値以下の場合、差の平方和が小さく、閾値よりも大きい場合、差の平方和は大きいと判断することができる。
【０１３１】
差の平方和が小さい場合、スペクトル比較手段１０２は、理想的なＺスペクトルＺidealによってＺスペクトルが再現されていると判断する。この場合、ステップＳＴ１３からステップＳＴ１４に進む。
【０１３２】
一方、差の平方和が大きい場合、スペクトル比較手段１０２は、理想的なＺスペクトルＺidealによってＺスペクトルが再現されていないと判断する。この場合、ステップＳＴ７に戻り、ベースライン項が再計算される。したがって、ステップＳＴ１３において、理想的なＺスペクトルＺidealによりＺスペクトルが再現されていると判断されるまで、ベースライン項の係数が再計算される。ステップＳＴ１３において、理想的なＺスペクトルＺidealによりＺスペクトルが再現されていると判断されたら、ステップＳＴ１４に進む。
【０１３３】
ステップＳＴ１４では、カウント手段１０３（図２参照）が、CPEスペクトルＦCPE（Δωａ）に含まれているCEST成分の総数ＴＮをカウントする。第１の形態では、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれていると判断されるたびにｎがインクリメントされるので（ステップＳＴ１１参照）、ｎの値が、CPEスペクトルＦCPE（Δωａ）に含まれるCEST成分の総数ＴＮを表している。つまり、CEST成分の総数ＴＮは、ＴＮ＝ｎとなる。したがって、CPEスペクトルＦCPE（Δωａ）に含まれているCEST成分の総数ＴＮをカウントすることができる。ここでは、説明の便宜上、ｎ＝２であるとする。したがって、CPEスペクトルＦCPE（Δωａ）に含まれているCEST成分の総数ＴＮは、ＴＮ＝２とカウントされる。ＴＮ＝２とカウントした後、ステップＳＴ１５に進む。
【０１３４】
ステップＳＴ１５では、CPEスペクトルＦCPE（Δωａ）の近似式（１６）で用いられるベースライン項の係数（ｃ０，ｃＭＴ）を、最終的に得られたベースライン項の係数の値に固定し、CEST項の係数を再計算する。ここでは、ｎ＝２であるとしているので、最終的に得られたベースライン項の係数の値は、近似式（２３）を用いて求められた値（ｃ０(2)’，ｃＭＴ(2)’）である（図１７参照）。したがって、ベースライン項の係数の値は、（ｃ０，ｃＭＴ）＝（ｃ０(2)’，ｃＭＴ(2)’）に固定され、CEST項の係数が再計算される。以下に、ステップＳＴ１５について、説明する。
【０１３５】
図２５は、ステップＳＴ１５のフローを示す図である。
ステップＳＴ１５１では、設定手段９４が、CEST項の数を表すｎを、初期値（ｎ＝１）に設定する。ｎ＝１に設定した後、ステップＳＴ１５２に進む。
【０１３６】
ステップＳＴ１５２では、第１のフィッティング手段９５が、CPEスペクトルの近似式（１７）のCEST項ＦＬ，１（Δωａ）に含まれる３つの係数（ａ１,１，ａ２,１，Δωｃ，１）の値を計算する。係数の値を計算する場合、第１のフィッティング手段９５は、先ず、式（１７）のベースライン項の係数（ｃ０，ｃＭＴ）を、（ｃ０，ｃＭＴ）＝（ｃ０(2)’，ｃＭＴ(2)’）に固定する。したがって、近似式（１７）では、CEST項ＦＬ，１（Δωａ）の３つの係数（ａ１,１，ａ２,１，Δωｃ，１）のみが未知の係数となる。第１のフィッティング手段９５は、CPEスペクトルＦCPE（Δωａ）と３つの未知の係数を含む近似式（１７）との誤差が最小になるように、フィッティングする。フィッティングを行う場合、第１のフィッティング手段９５は、先ず、CEST項ＦＬ，１（Δωａ）の係数（ａ１,１，ａ２,１，Δωｃ，１）の初期値を設定する。CEST項ＦＬ，１（Δωａ）の係数（ａ１,１，ａ２,１，Δωｃ，１）の初期値は、ステップＳＴ６で計算されたCEST項の係数の値を使用することができる。初期値を設定した後、第１のフィッティング手段は、初期値を基準にして係数の値を変更し、CPEスペクトルＦCPE（Δωａ）と式（１７）との誤差が最小になるときのCEST項ＦＬ，１（Δωａ）の３つの係数（ａ１,１，ａ２,１，Δωｃ，１）の値を計算する。これらの係数の値を計算した後、ステップＳＴ１５３に進む。
【０１３７】
ステップＳＴ１５３では、スペクトル計算手段９９が、式（２２）および式（２２ｃ）を用いて、スペクトルＦCPE＿１（Δωａ）を計算する。
【０１３８】
式（２２）のベースライン項の係数（ｃ０，ｃＭＴ）の値は、（ｃ０，ｃＭＴ）＝（ｃ０(2)’，ｃＭＴ(2)’）に固定されている（図１７参照）。また、式（２２ｃ）のCEST項の係数（ａ１,１，ａ２,１，Δωｃ，１）は、ステップＳＴ１５２において計算されている。したがって、式（２２）および（２２ｃ）に、これらの係数の値を代入することにより、スペクトルＦCPE＿１（Δωａ）を計算することができる。スペクトルＦCPE＿１（Δωａ）を計算した後、ステップＳＴ１５４に進む。
【０１３９】
ステップＳＴ１５４では、CEST判断手段１００が、CPEスペクトルＦCPE（Δωａ）に、CEST項ＦＬ，１（Δωａ）で表されるCEST成分（式１７ｃ参照）とは異なる他のCEST成分が含まれているか否かを判断する。他のCEST成分が含まれていないと判断された場合は、ステップＳＴ１５７に進む。一方、他のCEST成分が含まれていると判断された場合はステップＳＴ１５５に進む。
【０１４０】
ステップＳＴ１５５では、検出手段９３が、CPEスペクトルＦCPE（Δωａ）の他のCEST成分が現れるオフセット周波数を検出する。オフセット周波数を検出した後、ステップＳＴ１５６に進む。
【０１４１】
ステップＳＴ１５６では、設定手段９４が、ｎを、ｎ＝１からｎ＝２にインクリメントする。ｎ＝２に設定された場合、CPEスペクトルＦCPE（Δωａ）の近似式は、式（２３）で表される。ｎをインクリメントした後、ステップＳＴ１５２に戻る。
【０１４２】
したがって、ステップＳＴ１５４において、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれていると判断されるたびに、CPEスペクトルＦCPE（Δωａ）の近似式にCEST項が追加され、追加されたCEST項に含まれる係数の値が計算される。したがって、ベースライン項の係数（ｃ０，ｃＭＴ）を、（ｃ０，ｃＭＴ）＝（ｃ０(2)’，ｃＭＴ(2)’）に固定した場合のCEST項の係数の値を計算することができる。そして、ステップＳＴ１５４において、CPEスペクトルＦCPE（Δωａ）に他のCEST成分が含まれていないと判断されたら、ステップＳＴ１５７に進む。
【０１４３】
ステップＳＴ１５７では、カウント手段１０３が、CPEスペクトルＦCPE（Δωａ）に含まれているCEST成分の総数ＴＮをカウントする。CEST成分の総数ＴＮをカウントした後、ステップＳＴ１５８に進む。
【０１４４】
ステップＳＴ１５８では、カウント手段１０３が、ステップＳＴ１５７でカウントされたCEST成分の総数ＴＮが、ステップＳＴ１４でカウントされたCEST成分の総数ＴＮ（＝２）に等しいか否かを判断する。CEST成分の総数ＴＮが異なっていると判断された場合、CEST項の係数又はベースラインの係数の推定誤差が大きいと考えられる。そこで、CEST成分の総数ＴＮが異なっていると判断された場合、ステップＳＴ５（図７参照）に戻る。そして、テップＳＴ１５８においてCEST成分の総数が等しいと判断されるまで、ステップＳＴ５〜ＳＴ１５が繰り返し実行される。ステップＳＴ１５８において、CEST成分の総数ＴＮは等しいと判断されたら、ステップＳＴ１５を抜け出し、ステップＳＴ１６に進む。
【０１４５】
ステップＳＴ１６では、画像作成手段９０が、ステップＳＴ１〜ＳＴ１５の処理により得られたCEST項の係数の値に基づいてCEST画像を作成する。このようにして、フローが終了する。
【０１４６】
第１の形態では、ＺスペクトルをCPEスペクトルＦCPE（Δωａ）に変換する。CPEスペクトルＦCPE（Δωａ）のベースライン成分（CESTの影響を受けていない信号成分）は、偶関数で近似することができるので、ローレンツ関数による大きなピークを持たない。したがって、ＺスペクトルをCPEスペクトルＦCPE（Δωａ）に変換することにより、CEST成分を分離し易くすることができるので、CEST項に含まれる係数の値の推定誤差を小さくすることができる。
【０１４７】
また、第１の形態では、ステップＳＴ１３において、理想的なＺスペクトルＺidealとＺスペクトルとの差を求め、差の平方和が閾値を超えている場合は、ステップＳＴ７に戻ってベースライン項の係数を再計算する。したがって、ベースライン項の係数の推定誤差を更に小さくすることができる。
【０１４８】
また、第１の形態では、ステップＳＴ１５において、ベースライン項の係数の値を固定して、CEST項の係数の再計算している。したがって、CEST項の係数の推定誤差を更に小さくすることができる。
【０１４９】
また、ステップＳＴ１５８では、CEST成分の総数ＴＮが異なっていると判断した場合、ステップＳＴ５に戻り、CEST項の係数およびベースライン項の係数の値を再計算する。したがって、ベースライン項およびCEST項の係数の推定誤差を更に小さくすることができる。
【０１５０】
また、第１の形態では、CRZスペクトルＺＣＲＺの近似式（２０）は、定数項（ｂ０）とローレンツ関数の項との和で表されている。ローレンツ関数の項は３つの係数（ｂ１，ｂ２，Δω０）を含んでいるので、CRZスペクトルＺＣＲＺの近似式（２０）は、合計４つの係数（ｂ０，ｂ１，ｂ２，Δω０）を含んでいる。一方、CPEスペクトルＦCPE（Δωａ）の近似式（１６）で使用されるベースライン項は、定数項ｃ０と、ＭＴ項ｃＭＴＦ（Δωａ）の和で表されている。ＭＴ項のＦ（Δωａ）は、ローレンツ関数ではなくΔωａの二次関数（式（４）参照）であるので、ＭＴ項に含まれる係数は１個（ｃＭＴ）で済む。したがって、CPEスペクトルＦCPE（Δωａ）の近似式（１６）では、ベースライン項に含まれる係数の合計は２個（ｃ０，ｃＭＴ）で済む。このため、ＺスペクトルをCPEスペクトルＦCPE（Δωａ）に変換した場合、２つの係数（ｃ０，ｃＭＴ）の値を求めるだけで、CESTの影響を受けていない信号成分（ベースライン成分）を特定することができるので、フィッティングの精度を高めることができる。
【０１５１】
（２）第２の形態
第１の形態では、連続波のＲＦパルスを有するシーケンスを用いた例について説明したが、第２の形態では、複数のプリパレーションパルスを有するシーケンスを用いた例について説明する。
尚、第２の形態のＭＲ装置は、第１の形態のＭＲ装置１と同じである。
【０１５２】
図２６は、第２の形態におけるシーケンスＳＥｋを具体的に示す図である。
ｋ回目のシーケンスＳＥｋは、ｍ個のプリパレーションパルスと、データ収集セグメントＤＡＱとを有している。各プリパレーションパルスは、ＲＦパルスＸと、縦磁化を定常状態にするためのキラー勾配パルスＫとを含んでいる。ＲＦパルスＸの周波数ｆは、ｆ＝ｆｋに設定されている。プリパレーションパルスを繰り返し実行し、第ｍのプリパレーションパルスが実行された後に、シングルショット法によりデータを収集するためのデータ収集セグメントＤＡＱが実行される。
【０１５３】
図２６に示すシーケンスを用いた場合でも、連続波のＲＦを用いたシーケンス（図４参照）と同様に、式１１を用いてCPEスペクトルＦCPE（Δωａ）を得ることができる。したがって、第２の形態でも、第１の形態と同様の効果を得ることができる。
【０１５４】
（３）第３の形態
第３の形態では、ＲＦパルスの位相をサイクルさせるフェーズサイクリング法によりデータを収集するシーケンスを用いた例について説明する。
【０１５５】
尚、第３の形態のＭＲ装置は、第１の形態のＭＲ装置１と比較すると、プロセッサ９が実現する手段に違いがあるが、ハードウェアの構成は、第１の形態のＭＲ装置１と同じである。したがって、第３の形態のＭＲ装置の説明に当たっては、ハードウェア構成の説明は省略し、プロセッサについて主に説明する。
【０１５６】
図２７は、第３の形態においてプロセッサ９が実現する手段の説明図である。
第３の形態におけるプロセッサは、第１の形態におけるプロセッサと比較すると、以下の点（１）および（２）が異なっている。
【０１５７】
（１）Ｚスペクトル作成手段９１はＺスペクトルを作成する。ただし、第１の形態では、Ｚスペクトルの横軸はオフセット周波数であるが、第３の形態では、Ｚスペクトルの横軸は後述する位相差である。
（２）第３の形態では、プロセッサは周波数変換手段１０４を有している。周波数変換手段１０４は、位相差を周波数に変換する。
【０１５８】
尚、その他の点については、第３の形態におけるプロセッサは、第１の形態におけるプロセッサと同じであるので、説明は省略する。プロセッサ９は、記憶部１０に記憶されているプログラムを読み出し、プログラムに記述されている処理を実行するための手段９０〜１０４を実現する。
次に、第３の形態で使用されるシーケンスについて説明する。
【０１５９】
図２８は、第３の形態で使用されるシーケンスＳＥｋを具体的に示す図である。
ｋ回目のシーケンスＳＥｋは、第１〜第ｍのパルスセットSet１〜Setｍ、キラー勾配パルス、およびデータ収集セグメントＤＡＱを有している。以下では、先ず、第１〜第ｍのパルスセットSet１〜Setｍについて説明する。尚、第１〜第ｍのパルスセットSet１〜Setｍは同じ構成であるので、第１〜第ｍのパルスセットSet１〜Setｍの説明に当たっては、代表して第１のパルスセットSet１を取り上げて説明する。
【０１６０】
図２８には、第１のパルスセットSet１が拡大して示されている。
第１のパルスセットSet１は、ｒ個のＲＦパルスＸ１〜Ｘｒを有している。ＲＦパルスＸ１〜Ｘｒは、正のＲＦパルスと負のＲＦパルスが交互に現れるように構成されている。ＲＦパルスＸ１〜Ｘｒは、繰り返し時間Ｔ_iterで印加される。符号「Ｘ１」〜「Ｘｒ」の下に記載されている「φ１」〜「φｒ」は、ＲＦパルスの位相を表している。
【０１６１】
次に、ｒ個のＲＦパルスＸ１〜Ｘｒの位相φ１〜φｒについて説明する。先ず、ｒ個のＲＦパルスＸ１〜Ｘｒの中で、ｐ番目のＲＦパルスＸpと、ｐ＋１番目のＲＦパルスＸp+1について考える（尚、ｐは、１≦ｐ≦ｒ−１である）。ｐ番目のＲＦパルスＸpの位相を「φp」で表し、ｐ＋１番目のＲＦパルスＸp+1の位相を「φp+1」で表すと、ｋ回目のシーケンスＳＥｋにおけるＲＦパルスの位相差Δφ（ｋ）＝φp+1−φpは、以下の式を満たすように設定されている。
【数３５】
0006348449.tif 000036
【０１６２】
式（３０）から、位相差Δφ（ｋ）は、ｋの値に応じて変化するように設定されていることが分かる。
【０１６３】
図２８では、第１のパルスセットSet１について示されているが、第２〜第ｍのパルスセットSet２〜Setｍも、第１のパルスセットSet１と同じ構成である。したがって、どのパルスセットも、ｒ個のＲＦパルスＸ１〜Ｘｒを有しており、ＲＦパルスの位相差Δφ（ｋ）は式（３０）を満たすように設定されている。
【０１６４】
第１〜第ｍのパルスセットSet１〜Setｍを印加した後、横磁化を消失させるためのキラー勾配パルスを印加する。そして、キラー勾配パルスを印加した後、データを収集するためのデータ収集セグメントＤＡＱが実行される。ここでは、データ収集セグメントＤＡＱは、シングルショット法でデータを収集するとする。
【０１６５】
ｋ回目のシーケンスＳＥｋは、上記のように構成されている。第３の形態では、シーケンスＳＥｋがｒ回実行される。尚、実行されるシーケンスの回数ｒが多いほど、周波数分解能の高いＺスペクトルが得られるので、ｒはある程度大きい値であることが望ましい。一般的には、ｒ＝１６〜３２に設定することが考えられる。
【０１６６】
ここで、フェーズサイクリング法のシーケンスを使用する場合のＦ（Δωａ）（式（４）参照）について考察する。
Δωａは、Miyoshi等により、以下の周期関数で置き換えられることが示されている（Miyoshi M, et al., Proceedings of ISMRM2014, #3299）。
【数３６】
0006348449.tif 000037
【０１６７】
そこで、式（４）において、Δωａ２を２（１−cosΔφａ）／Ｔｉｔｅｒ２に置き換えると、以下の式が得られる。
【数３７】
0006348449.tif 000038
【０１６８】
式（３２）で定義されたＦ（Δφａ）は、式（４）で定義されたＦ（Δωａ）と同様に偶関数である。したがって、式（４）の代わりに式（３２）を用いても、図６に示すように、周波数ΔωＣの近傍において、CPEスペクトルのCEST成分Ｑ１とベースライン成分Ｑ２との比Ｒ（＝Ｑ１／Ｑ２）を大きくすることができる。このため、CPEスペクトルから、CEST成分Ｑ１を分離し易くすることができる。
【０１６９】
尚、第１の形態では、周波数ΔωａおよびΔωｃ，ｎを変数として含む式を用いて、各種スペクトルが算出されている。しかし、第３の形態では、周波数の代わりに位相差を用いているので（式３２参照）、周波数を位相差に置き換えた式を用いて各スペクトル（CPEスペクトルなど）を求める必要がある。具体的には、周波数ΔωａおよびΔωｃ，ｎは、それぞれ、以下の位相差ΔφａおよびΔφｃ，ｎに置き換えればよい。
【数３８】
0006348449.tif 000039
【０１７０】
位相差ΔφａおよびΔφｃ，ｎは、以下の式で表される。
【数３９】
0006348449.tif 000040
【０１７１】
第３の形態のＭＲ装置では、プロセッサは、位相差を周波数に変換する周波数変換手段１０４（図２７参照）を備えている。周波数変換手段１０４は、式（３３）および（３４）に基づいて、位相差を周波数に変換することができる。したがって、フェーズサイクル法を用いる場合であっても、各スペクトルの位相差を周波数に変換できることが分かる。
【符号の説明】
【０１７２】
１  ＭＲ装置
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  送信器
６  勾配磁場電源
７  受信器
８  コンピュータ
９  プロセッサ
１０  記憶部
１１  操作部
１２  表示部
１３  被検体
２１  収容空間
９０  画像作成手段
９１  Ｚスペクトル作成手段
９２  スペクトル変換手段
９３  検出手段
９４  設定手段
９５  第１のフィッティング手段
９６  CRZスペクトル作成手段
９７  第２のフィッティング手段
９８  係数値計算手段
９９  スペクトル計算手段
１００  CEST判断手段
１０１  スペクトル推定手段
１０２  スペクトル比較手段
１０３  カウント手段
【図１】
0006348449.tif 000041
【図２】
0006348449.tif 000042
【図３】
0006348449.tif 000043
【図４】
0006348449.tif 000044
【図５】
0006348449.tif 000045
【図６】
0006348449.tif 000046
【図７】
0006348449.tif 000047
【図８】
0006348449.tif 000048
【図９】
0006348449.tif 000049
【図１０】
0006348449.tif 000050
【図１１】
0006348449.tif 000051
【図１２】
0006348449.tif 000052
【図１３】
0006348449.tif 000053
【図１４】
0006348449.tif 000054
【図１５】
0006348449.tif 000055
【図１６】
0006348449.tif 000056
【図１７】
0006348449.tif 000057
【図１８】
0006348449.tif 000058
【図１９】
0006348449.tif 000059
【図２０】
0006348449.tif 000060
【図２１】
0006348449.tif 000061
【図２２】
0006348449.tif 000062
【図２３】
0006348449.tif 000063
【図２４】
0006348449.tif 000064
【図２５】
0006348449.tif 000065
【図２６】
0006348449.tif 000066
【図２７】
0006348449.tif 000067
【図２８】
0006348449.tif 000068