JP2024016804A - ｐＨ算出装置、磁気共鳴イメージング装置、およびｐＨ算出方法 - Google Patents

Abstract

【課題】算出されるｐＨ値の正確性（精度）を向上すること。【解決手段】本実施形態に係るｐＨ算出装置は、被検体に対しケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ撮像により生成され前記複数種の物質にそれぞれ対応する複数の信号値の低下ポイントを含む少なくとも一つのＺスペクトルと、前記被検体に関する情報とを取得する取得部と、前記複数の信号値の低下ポイントに対応する複数のＭＲ信号値を用いて、ｐＨ依存値を算出する算出部と、前記ｐＨ依存値に対するｐＨ値を示す検量線を、前記被検体に関する情報に応じて予め用意された複数の検量線から、前記被検体に関する情報に応じて選択する選択部と、前記選択された検量線と前記ｐＨ依存値とに基づいて、ｐＨ値を決定するｐＨ値決定部と、を備える。【選択図】図４

Description

本明細書及び図面に開示の実施形態は、ｐＨ算出装置、磁気共鳴イメージング装置、およびｐＨ算出方法に関する。

従来の技術として、化学交換飽和移動（以下、ＣＥＳＴ（Ｃｈｅｍｉｃａｌ Ｅｘｃｈａｎｇｅ Ｓａｔｕｒａｔｉｏｎ Ｔｒａｎｓｆｅｒ）と呼ぶ）を用いたＭＲＩによるイメージング手法が提唱されている。水に溶けている溶質のプロトンは、水のプロトンと化学的に交換される。このような化学交換の速度は、温度やｐＨに依存している。一方で、プロトンの共鳴周波数は、プロトンの状態によって変化する。このようなプロトンの共鳴周波数の変化は、ケミカルシフトと呼ばれている。これらのことから、溶質のプロトンに特異的な周波数に設定して飽和ＲＦパルスを送信することで、本来は飽和していないはずの水のプロトンが飽和状態になってしまう現象が発生する。このような現象を化学交換飽和移動（ＣＥＳＴ）と呼ぶ。ＣＥＳＴ現象を利用したイメージング手法は、ＣＥＳＴイメージングと呼ばれる。このとき、ＣＥＳＴ現象を実現する物質（ＣＥＳＴ物質）が、ＣＥＳＴイメージングにおいて用いられる。

ＣＥＳＴ現象は、上述のように、温度やｐＨなどの物質の特性を反映している。このため、例えば、２つのプロトンプールを持つＣＥＳＴ物質とレシオメトリック（Ｒａｔｉｏｍｅｔｒｉｃ）法という手法を組み合わせることで、ＣＥＳＴ物質の濃度に依らずにｐＨを撮像するｐＨイメージングが可能とされている。しかしながら、レシオメトリック法によるｐＨの算出は、ＣＥＳＴイメージングにおいて実行された飽和パルスの強度（実行飽和パルス強度と呼ぶ）および温度が一定等の制約がある。生体内では温度は概ね３７℃で一定であると言えるが、実行飽和パルス強度は、ＣＥＳＴイメージング中において一定であるとは言えない場合がある。このため、算出されたｐＨ値は正確ではない恐れがある。

特表２０１９－５０５２４８号公報

Ｋ．Ｍ． Ｗａｒｄ ａｎｄ Ｒ．Ｓ． Ｂａｌａｂａｎ Ｄｅｔｅｒｍｉｎａｔｉｏｎ ｏｆ ｐＨ Ｕｓｉｎｇ Ｗａｔｅｒ Ｐｒｏｔｏｎｓ ａｎｄ Ｃｈｅｍｉｃａｌ Ｅｘｃｈａｎｇｅ Ｄｅｐｅｎｄｅｎｔ Ｓａｔｕｒａｔｉｏｎ Ｔｒａｎｓｆｅｒ （ＣＥＳＴ）、 Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ ｉｎ Ｍｅｄｉｃｉｎｅ ４４：７９９－８０２ （２０００）

本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、算出されるｐＨ値の正確性（精度）を向上することである。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。

本実施形態に係るｐＨ算出装置は、取得部と、算出部と、選択部と、ｐＨ値決定部とを備える。取得部は、被検体に対しケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ（Ｃｈｅｍｉｃａｌ Ｅｘｃｈａｎｇｅ Ｓａｔｕｒａｔｉｏｎ Ｔｒａｎｓｆｅｒ）撮像により生成され前記複数種の物質にそれぞれ対応する複数の信号値の低下ポイントを含む少なくとも一つのＺスペクトルと、前記被検体に関する情報とを取得する。算出部は、前記複数の信号値の低下ポイントに対応する複数のＭＲ信号値を用いて、ｐＨ依存値を算出する。選択部は、前記ｐＨ依存値に対するｐＨ値を示す検量線を、前記被検体に関する情報に応じて予め用意された複数の検量線から、前記被検体に関する情報に応じて選択する。ｐＨ値決定部は、前記選択された検量線と前記ｐＨ依存値とに基づいて、ｐＨ値を決定する。

図１は、第１実施形態に係るＭＲＩ装置の構成を示すブロック図。 図２は、第１実施形態に係り、造影剤を用いたＣＥＳＴイメージングに関する概要の一例を示す概要図。 図３は、第１実施形態に係り、Ｂ_１値に応じた複数の検量線の一例を示す図。 図４は、第１実施形態に係り、ｐＨ値決定処理の手順の一例を示すフローチャート。 図５は、第２実施形態に係り、ｐＨ算出装置の構成の一例を示す図。 図６は、第３実施形態に係り、ファントムの温度が２５℃である場合の３つのｐＨ値に応じた３つのＺスペクトルの一例を示す図。 図７は、第３実施形態に係り、ファントムの温度が３６℃である場合の３つのｐＨ値に応じた３つのＺスペクトルの一例を示す図。 図８は、第３実施形態に係るｐＨ値決定処理の手順の一例を示すフローチャート。 図９は、第４実施形態に係るｐＨ値決定処理の手順の一例を示すフローチャート。

以下、図面を参照しながら、ｐＨ算出装置、当該ｐＨ算出装置を搭載した磁気共鳴イメージング装置（以下、ＭＲＩ（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）装置と呼ぶ）、およびｐＨ算出方法の実施形態について説明する。なお、実施形態は、以下の実施形態に限られるものではない。また、各実施形態において説明する内容は、原則として、他の実施形態においても同様に適用することができる。以下の実施形態では、同一の参照符号を付した部分は同様の動作をおこなうものとして、重複する説明は適宜省略する。

（第１実施形態）
図１は、第１実施形態に係るＭＲＩ装置１００の構成を示すブロック図である。図１に示すように、ＭＲＩ装置１００は、静磁場磁石１０１と、静磁場電源１０２と、傾斜磁場コイル１０３と、傾斜磁場電源１０４と、寝台１０５と、寝台制御回路１０６と、送信コイル１０７と、送信回路１０８と、受信コイル１０９と、受信回路１１０と、シーケンス制御回路１２０と、コンピューター１３０（画像処理装置とも称される）とを備える。なお、ＭＲＩ装置１００に、被検体Ｐ（例えば、人体）は含まれない。また、図１に示す構成は一例に過ぎない。例えば、シーケンス制御回路１２０及びコンピューター１３０内の各部は、適宜統合若しくは分離して構成されてもよい。

静磁場磁石１０１は、中空の略円筒形状に形成された磁石であり、内部の空間に静磁場を発生する。静磁場磁石１０１は、例えば、超伝導磁石等であり、静磁場電源１０２から電流の供給を受けて励磁する。静磁場電源１０２は、静磁場磁石１０１に電流を供給する。なお、静磁場磁石１０１は、永久磁石でもよく、この場合、ＭＲＩ装置１００は、静磁場電源１０２を備えなくてもよい。また、静磁場電源１０２は、ＭＲＩ装置１００とは別に備えられてもよい。

傾斜磁場コイル１０３は、中空の略円筒形状に形成されたコイルであり、静磁場磁石１０１の内側に配置される。傾斜磁場コイル１０３は、互いに直交するＸ、Ｙ、及びＺの各軸に対応する３つのコイルが組み合わされて形成されており、これら３つのコイルは、傾斜磁場電源１０４から個別に電流の供給を受けて、Ｘ、Ｙ、及びＺの各軸に沿って磁場強度が変化する傾斜磁場を発生する。傾斜磁場コイル１０３によって発生するＸ、Ｙ、及びＺの各軸の傾斜磁場は、例えば、スライス用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ、及びリードアウト用傾斜磁場Ｇｒである。傾斜磁場電源１０４は、傾斜磁場コイル１０３に電流を供給する。

寝台１０５は、被検体Ｐが載置される天板１０５ａを備え、寝台制御回路１０６による制御の下、天板１０５ａを、被検体Ｐが載置された状態で、傾斜磁場コイル１０３の空洞（撮像口）内へ挿入する。通常、寝台１０５は、長手方向が静磁場磁石１０１の中心軸と平行になるように設置される。寝台制御回路１０６は、コンピューター１３０による制御の下、寝台１０５を駆動して天板１０５ａを長手方向及び上下方向へ移動する。

送信コイル１０７は、傾斜磁場コイル１０３の内側に配置され、送信回路１０８からＲＦパルスの供給を受けて、高周波磁場を発生する。送信回路１０８は、対象とする原子の種類及び磁場強度で定まるラーモア（Ｌａｒｍｏｒ）周波数に対応するＲＦパルスを送信コイル１０７に供給する。

受信コイル１０９は、傾斜磁場コイル１０３の内側に配置され、高周波磁場の影響によって被検体Ｐから発せられる磁気共鳴信号（以下、ＭＲ（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ）信号と呼ぶ）を受信する。受信コイル１０９は、ＭＲ信号を受信すると、受信したＭＲ信号を受信回路１１０へ出力する。

なお、上述した送信コイル１０７及び受信コイル１０９は一例に過ぎない。送信コイル１０７及び受信コイル１０９は、送信機能のみを備えたコイル、受信機能のみを備えたコイル、若しくは送受信機能を備えたコイルのうち、１つ若しくは複数を組み合わせることによって構成されてもよい。

受信回路１１０は、受信コイル１０９から出力されるＭＲ信号を検出し、検出したＭＲ信号に基づいてＭＲデータを生成する。具体的には、受信回路１１０は、受信コイル１０９から出力されるＭＲ信号をデジタル変換することによってＭＲデータを生成する。また、受信回路１１０は、生成したＭＲデータをシーケンス制御回路１２０へ送信する。なお、受信回路１１０は、静磁場磁石１０１や傾斜磁場コイル１０３等を備える架台装置側に備えられてもよい。

シーケンス制御回路１２０は、コンピューター１３０から送信されるシーケンス情報に基づいて、傾斜磁場電源１０４、送信回路１０８及び受信回路１１０を駆動することによって、被検体Ｐの撮像を行う。ここで、シーケンス情報は、撮像を行うための手順を定義した情報であって、シーケンス条件とも称される。シーケンス情報には、傾斜磁場電源１０４が傾斜磁場コイル１０３に供給する電流の強さや電流を供給するタイミング、送信回路１０８が送信コイル１０７に供給するＲＦパルスの強度やＲＦパルスを印加するタイミング、受信回路１１０がＭＲ信号を検出するタイミング等が定義される。

例えば、シーケンス制御回路１２０は、特定用途向け集積回路（ＡＳＩＣ：Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）、フィールドプログラマブルゲートアレイ（ＦＰＧＡ：Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）等の集積回路、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＭＰＵ（Ｍｉｃｒｏ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）等の電子回路である。シーケンス制御回路１２０は、シーケンス制御部に相当する。例えば、シーケンス制御回路１２０は、後述のＣＥＳＴ撮像とは異なる撮像（例えば、後述のマップイメージング）とＣＥＳＴ撮像とを実行する。

なお、シーケンス制御回路１２０は、傾斜磁場電源１０４、送信回路１０８及び受信回路１１０を駆動して被検体Ｐを撮像した結果、受信回路１１０からＭＲデータを受信すると、受信したＭＲデータをコンピューター１３０へ転送する。

コンピューター１３０は、ＭＲＩ装置１００の全体制御や、画像の生成等を行う。コンピューター１３０は、記憶回路１３２、入力装置１４１、ディスプレイ１４３、処理回路１５０を備える。処理回路１５０は、インタフェース機能１３１、制御機能１３３、画像生成機能１３４、取得機能１３６、特定機能１３８、算出機能１４０、選択機能１４２、及びｐＨ値決定機能１４４を備える。

インタフェース機能１３１、制御機能１３３、画像生成機能１３４、取得機能１３６、特定機能１３８、算出機能１４０、選択機能１４２、及びｐＨ値決定機能１４４にて行われる各処理機能は、コンピューター１３０によって実行可能なプログラムの形態で記憶回路１３２へ記憶されている。処理回路１５０はプログラムを記憶回路１３２から読み出し、実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路１５０は、図１の処理回路１５０内に示された各機能を有することになる。

なお、図１においては単一の処理回路１５０にて、インタフェース機能１３１、制御機能１３３、画像生成機能１３４、取得機能１３６、特定機能１３８、算出機能１４０、選択機能１４２、及びｐＨ値決定機能１４４にて行われる処理機能が実現されるものとして説明したが、複数の独立したプロセッサを組み合わせて処理回路１５０を構成し、各プロセッサがプログラムを実行することにより機能を実現するものとしても構わない。換言すると、上述のそれぞれの機能がプログラムとして構成され、１つの処理回路１５０が各プログラムを実行する場合であってもよいし、特定の機能が専用の独立したプログラム実行回路に実装される場合であってもよい。

上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ、ＧＰＵ（Ｇｒａｐｈｉｃａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）或いは、特定用途向け集積回路、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Ｓｉｍｐｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＳＰＬＤ）、複合プログラマブル論理デバイス（Ｃｏｍｐｌｅｘ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイＦＰＧＡ））等の回路を意味する。プロセッサは記憶回路１３２に保存されたプログラムを読み出し実行することで機能を実現する。

なお、記憶回路１３２にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、寝台制御回路１０６、送信回路１０８、受信回路１１０等も同様に、上記のプロセッサ等の電子回路により構成される。

処理回路１５０は、インタフェース機能１３１により、シーケンス情報をシーケンス制御回路１２０へ送信し、シーケンス制御回路１２０からＭＲデータを受信する。また、インタフェース機能１３１を有する処理回路１５０は、ＭＲデータを受信すると、受信したＭＲデータを記憶回路１３２に格納する。インタフェース機能１３１を実現する処理回路１５０は、インタフェース部に相当する。処理回路１５０における他の機能については、後ほど説明する。

記憶回路１３２は、インタフェース機能１３１を有する処理回路１５０によって受信されたＭＲデータや、取得機能１３６により取得された各種データ、画像生成機能１３４により生成された各種画像データ、算出機能１４０において用いられる演算処理、演算処理により算出されたｐＨ依存値、選択機能１４２において用いられる複数の検量線、ｐＨ値決定機能１４４により決定されたｐＨ値などを記憶する。複数の検量線各々は、ｐＨ依存値に対するｐＨ値を示すものであって、複数のＢ_１値に応じて異なる。記憶回路１３２は、例えば、複数の対応表（ＬＵＴ：Ｌｏｏｋ Ｕｐ Ｔａｂｌｅ）の形式で、複数の検量線を、複数のＢ_１値と対応付けて記憶する。

複数の検量線は、例えば、ファントムなどにより予め計測された実測値により予め用意されて、記憶回路１３２に記憶される。実測値は、例えば、Ｂ_１値、ｐＨ依存値、ｐＨ値などである。複数の検量線は、実測に基づくデータ、シミュレーションで計算されたデータ、または当該実測値から補間されたデータで、予め用意されたデータであってもよい。本実施形態における補間は、隣接する２つのＢ_１値の間を補間して検量線を算出することにある。

また、記憶回路１３２は、制御機能１３３によってｋ空間に配置されたＭＲデータ（ｋ空間データとも称される）を記憶する。これらの記憶される各種データについては、後ほど説明する。例えば、記憶回路１３２は、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等で実現される。

入力装置１４１は、ユーザからの各種指示や情報入力を受け付ける。入力装置１４１は、例えば、トラックボール、スイッチボタン、マウス、キーボード、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、及び音声入力回路等によって実現される。入力装置１４１は、処理回路１５０に電気的に接続されており、操作者から受け取った入力操作を電気信号へ変換し処理回路１５０へと出力する。

なお、本明細書において入力装置１４１は、マウス、キーボードなどの物理的な操作部品（入力インタフェース）を備えるものだけに限られない。例えば、ＭＲＩ装置１００とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を制御回路へ出力する電気信号の処理回路も入力装置１４１の例に含まれる。

ディスプレイ１４３は、制御機能１３３を有する処理回路１５０による制御の下、撮像条件の入力を受け付けるためのＧＵＩ（Ｇｒａｐｈｉｃａｌ Ｕｓｅｒ Ｉｎｔｅｒｆａｃｅ）や、画像生成機能１３４を有する処理回路１５０によって生成された画像等を表示する。ディスプレイ１４３は、例えば、ＣＲＴディスプレイや液晶ディスプレイ、有機ＥＬディスプレイ、ＬＥＤディスプレイ、プラズマディスプレイ、又は当技術分野で知られている他の任意のディスプレイ、モニタ等の表示デバイスにより実現される。

本実施形態に関する化学交換飽和移動（ＣＥＳＴ：Ｃｈｅｍｉｃａｌ Ｅｘｃｈａｎｇｅ Ｓａｔｕｒａｔｉｏｎ Ｔｒａｓｆｅｒ）効果について簡単に説明する。自由水（バルク水）におけるプロトンと化合物のプロトンとは、交換される。例えば、アミド基（－ＮＨ）、ヒドロキシル基（－ＯＨ）、及びアミノ基（－ＮＨ２）等の、化合物のプロトンは、交換性プロトンである。ＣＥＳＴ効果は、「Ｚスペクトル」や、「ＭＴＲａｓｙｍスペクトル」と呼ばれる量を用いて記述される。ＣＥＳＴ効果に関する撮像は、アミド基（－ＮＨ）、ヒドロキシル基（－ＯＨ）、及びアミノ基（－ＮＨ２）等のプロトンと自由水におけるプロトンとが交換されることを利用したＭＲイメージング手法（以下、ＣＥＳＴイメージングと呼ぶ）である。

ＣＥＳＴ効果に関する磁気共鳴イメージング手法（以下、ＣＥＳＴ撮像と呼ぶ）において、シーケンス制御回路１２０は、自由水の共鳴周波数から離れた周波数（ｏｆｆ－ｒｅｓｏｎａｎｃｅな周波数）であって、交換性プロトン（例えば、化合物のプロトン）の共鳴周波数で、周波数選択ＲＦ（Ｒａｄｉｏ Ｆｒｅｑｕｅｎｃｙ）パルスである飽和パルスを、ＭＲ信号の収集前において被検体Ｐに対して印加する。飽和パルスは、プリサチュレーションパルス（Ｐｒｅｓａｔｕｒａｔｉｏｎ ｐｕｌｓｅ）とも称される。具体的には、ＣＥＳＴ撮像は、複数の飽和パルスの印加と複数の飽和パルスに応じた複数の磁気共鳴信号（ＭＲ信号）の収集とを実行する磁気共鳴イメージング手法である。

本実施形態に関するＣＥＳＴイメージングについて、図２を用いて説明する。図２は、造影剤を用いたＣＥＳＴイメージングに関する概要の一例を示す概要図である。造影剤は、ケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む。以下、説明を具体的にするために、造影剤は、イオパミドールであるものとする。このとき、ケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質は、４．２ｐｐｍのケミカルシフトに対応する２つのアミド基と、５．６ｐｐｍのケミカルシフトに対応する１つのアミド基との２種の物質に相当する。なお、造影剤は、イオパミドールに限定されず、例えば、イオプロミドであってもよい。

以下、説明を具体的にするために、ＣＥＳＴイメージングに用いられる複数の飽和パルスの周波数は、－８ｐｐｍから８ｐｐｍの範囲において、０．１ｐｐｍ刻みであるものとする。例えは、静磁場強度が３Ｔである場合、０ｐｐｍに対応する飽和パルスの周波数は、静磁場強度に基づく自由水の共鳴周波数（以下、中心周波数と呼ぶ）であって、１２８ＭＨｚである。このとき、＋８ｐｐｍに対応する飽和パルスの周波数は、（１２８ＭＨｚ＋１２８×８Ｈｚ）となる。一方、－８ｐｐｍに対応する飽和パルスの周波数は、（１２８ＭＨｚ－１２８×８Ｈｚ）となる。

図２に示すように、ＣＥＳＴイメージングによれば、中心周波数を０ｐｐｍとして－８ｐｐｍから８ｐｐｍまで０．１ｐｐｍ刻みの複数の飽和パルスの印加に伴ってＭＲ画像を取得する。１６１枚のＭＲ画像に基づいて、Ｂ_０（静磁場）不均一の補正および当該補正後の正確なケミカルシフトにおける化学交換飽和移動（ＣＥＳＴ）の効果の影響を示すＺスペクトルが生成される。

図２に示す領域ＮＣＲは、造影剤では造影されていない領域（以下、非造影ＲＯＩ（Ｒｅｇｉｏｎ ｏｆ Ｉｎｔｅｒｅｓｔ：関心領域）とよぶ）を示している。非造影ＲＯＩにおけるＺスペクトルＮＣＺでは、ＭＲ信号の低下は、飽和パルスの周波数（０ｐｐｍ）を中心とする所定の周波数の範囲で生じている。一方、図２に示す領域ＣＲは、造影剤により造影された領域（以下、造影ＲＯＩとよぶ）を示している。図２に示すＺスペクトルＣＺは、ケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴイメージング（以下、造影ＣＥＳＴイメージングと呼ぶ）により生成されたＺスペクトル（以下、造影Ｚスペクトルと呼ぶ）を示している。造影ＲＯＩにおけるＺスペクトルＣＺでは、ＭＲ信号の低下は、飽和パルスの周波数（０ｐｐｍ）を中心するとする所定の周波数の範囲に加えて、４．２ｐｐｍと５．６ｐｐｍとの位置で生じている。

図２において、造影ＲＯＩに関する造影ＺスペクトルＣＺにおける点線ＮＣＣは、造影ＲＯＩと同一のＲＯＩに対して造影剤を用いずにＣＥＳＴイメージング（以下、非造影ＣＥＳＴイメージングと呼ぶ）を実行することで生成されたＺスペクトル（以下、非造影Ｚスペクトルと呼ぶ）を示している。図２に示すように、造影ＲＯＩに関するＺスペクトルＣＺと点線ＮＣＣとの差は、両矢印で示されている。造影ＲＯＩに関するＺスペクトルＣＺに関して、４．２ｐｐｍと５．６ｐｐｍとの位置における両矢印は、造影剤の成分のイオパミドールにおける２種類のアミド基に対応する信号値の低下を示している。ケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤としてイオパミドールが用いられる場合、ｐＨ依存値は、図２に示すような４．２ｐｐｍと５．６ｐｐｍとの位置における信号値の低下の成分を用いて、算出機能１４０により算出される。

図３は、Ｂ_１値に応じた複数の検量線の一例を示す図である。検量線（標準曲線：Ｓｔａｎｄａｒｄ Ｃｕｒｖｅ）は、校正曲線：Ｃａｌｉｂｒａｔｉｏｎ Ｃｕｒｖｅと称されてもよい。Ｂ_１値は、例えば、シーケンス情報として設定されたＲＦパルスの強度を基準としたずれを、マイクロテスラ（μＴ）の単位で示した値である。すなわち、Ｂ_１値は、被検体ＰにＲＦパルスを照射して得られる。このため、Ｂ_１値は、被検体Ｐから取得された情報に対応する。図３に示すように、複数の検量線各々は、ｐＨ依存値に対するｐＨ値を示している。複数の検量線各々は、既知のｐＨを有するファントムを撮影して算出されたｐＨ依存値と当該ｐＨとを用いて決定され、記憶回路１３２に予め記憶される。Ｂ_１値を撮像領域に亘って配置したＢ_１マップの生成については後ほど説明する。

図３に示すように、縦軸をｐＨ依存値として横軸をｐＨとする座標系において、検量線の傾きは、Ｂ_１値が増加するにつれて減少する。換言すれば、検量線の傾きは、Ｂ_１値が減少するにつれて増加する。図３では、４つのＢ_１値にそれぞれ対応する４つの検量線が示されているが、記憶回路１３２に記憶される検量線の数は４つに限定されない。また、複数の検量線に関する複数のＢ_１値は、図３では０．５μＴの等間隔で示されているが、複数のＢ_１値の間隔は、不等間隔であってよい。

処理回路１５０は、制御機能１３３により、ＭＲＩ装置１００の全体制御を行い、撮像や画像の生成、画像の表示等を制御する。例えば、制御機能１３３を有する処理回路１５０は、撮像条件（撮像パラメータ等）の入力をＧＵＩ上で受け付け、受け付けた撮像条件により設定された飽和パルスの条件に従ってシーケンス情報を生成する。また、制御機能１３３を有する処理回路１５０は、生成したシーケンス情報をシーケンス制御回路１２０へ送信する。制御機能１３３を実現する処理回路１５０は、制御部に相当する。本実施形態では、制御機能１３３は、非造影ＣＥＳＴイメージングのシーケンス情報と、Ｂ_０マップおよびＢ_１マップの生成に関するＭＲデータを収集するＭＲ撮像（以下、マップイメージングと呼ぶ）のシーケンス情報と、造影ＣＥＳＴイメージングのシーケンス情報をそれぞれシーケンス制御回路１２０へ送信する。

処理回路１５０は、画像生成機能１３４により、ｋ空間データを記憶回路１３２から読み出し、読み出したｋ空間データにフーリエ変換等の再構成処理を施すことで、画像を生成する。例えば、画像生成機能１３４は、マップイメージングにより収集されたＭＲデータ（以下、マップＭＲデータと呼ぶ）に基づいて、Ｂ_０マップおよびＢ_１マップを生成する。Ｂ_０マップは、撮像領域における静磁場（Ｂ_０）の不均一性を示すマップである。Ｂ_１マップは、撮像領域亘ってＢ_１値を配置したマップである。Ｂ_１マップは、撮像領域における飽和パルスの強度（実行飽和パルス強度ともいう）の不均一性を示している。Ｂ_０マップおよびＢ_１マップの生成は、既知の手法を適宜利用可能であるため、説明は省略する。

処理回路１５０は、画像生成機能１３４により、非造影ＣＥＳＴイメージングにより収集されたＭＲデータ（以下、非造影データと呼ぶ）に基づいて、複数のＭＲ画像（以下、非造影ＭＲ画像と呼ぶ）を生成する。画像生成機能１３４は、複数の非造影ＭＲ画像とＢ_０マップとに基づいて、飽和パルスの位置の補正（以下、Ｂ_０補正と呼ぶ）を伴って、非造影Ｚスペクトルを生成する。なお、Ｂ_０補正は、非造影ＭＲ画像に対して実行されてもよい。このとき、画像生成機能１３４は、Ｂ_０補正が実施された複数の非造影ＭＲ画像に基づいて、非造影Ｚスペクトルを生成する。

処理回路１５０は、画像生成機能１３４により、造影ＣＥＳＴイメージングにより収集されたＭＲデータ（以下、造影データと呼ぶ）に基づいて、複数のＭＲ画像（以下、造影ＭＲ画像と呼ぶ）を生成する。画像生成機能１３４は、複数の造影ＭＲ画像とＢ_０マップとに基づいて、Ｂ_０補正を伴って、造影Ｚスペクトルを生成する。なお、Ｂ_０補正は、造影ＭＲ画像に対して実行されてもよい。このとき、画像生成機能１３４は、Ｂ_０補正が実施された複数の造影ＭＲ画像に基づいて、造影Ｚスペクトルを生成する。以上のことから、実施形態におけるＺスペクトルは、異なる撮像により生成されたＢ_０マップに基づいてＣＥＳＴ撮像における飽和パルスの位置が補正されたスペクトルデータである。画像生成機能１３４を実現する処理回路１５０は、画像生成部に相当する。

処理回路１５０は、取得機能１３６により、ケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ撮像により生成され当該複数種の物質にそれぞれ対応する複数の信号値の低下ポイント（または極小値）を含む少なくとも一つのＺスペクトルと、被検体Ｐに関する情報とを取得する。本実施形態において、被検体Ｐに関する情報は、例えば、Ｂ_１マップである。すなわち、Ｂ_１マップは、被検体Ｐから取得された情報である。このとき、取得機能１３６は、少なくとも一つのＺスペクトルとＢ_１マップとを取得する。具体的には、取得機能１３６は、画像生成機能１３４から、例えば、非造影Ｚスペクトルと、造影Ｚスペクトルと、Ｂ_１マップとを取得する。被検体Ｐに関する情報のうちＢ_１マップとは異なる情報については、第３実施形態において説明する。

処理回路１５０がＭＲＩ装置１００とは別個のｐＨ算出装置などに搭載される場合、処理回路１５０は、取得機能１３６により、他のＭＲＩ装置から、例えば、非造影Ｚスペクトルと、造影Ｚスペクトルと、Ｂ_１マップとを取得する。この場合においても、造影Ｚスペクトルと非造影Ｚスペクトルとは、Ｂ_０マップを用いたＢ_０補正が実施されたＺスペクトルである。取得機能１３６を実現する処理回路１５０は、取得部に相当する。

処理回路１５０は、特定機能１３８により、Ｚスペクトルに対応する、Ｂ_１マップ上の位置でのＢ_１値を特定する。具体的には、特定機能１３８は、造影Ｚスペクトルと非造影Ｚスペクトルとが生成されたＢ_１マップ上の位置を特定し、当該特定された位置でのＢ_１値を特定する。例えば、造影Ｚスペクトルと非造影Ｚスペクトルとが、図２に示す造影ＲＯＩについて１つずつ生成された場合、特定機能１３８は、造影ＲＯＩに含まれる複数の画素の位置を特定する。次いで、特定機能１３８は、特定された複数の画素に対応する複数のＢ_１値に対して平均を計算することにより、造影Ｚスペクトルと非造影Ｚスペクトルとに対応するＢ_１値を特定する。

また、造影Ｚスペクトルと非造影Ｚスペクトルとが、図２に示す造影ＲＯＩに含まれる複数の画素各々に対して生成された場合についても、特定機能１３８は、造影ＲＯＩに含まれる複数の画素の位置を特定する。次いで、特定機能１３８は、画素の位置に従って造影Ｚスペクトルと非造影Ｚスペクトルとに対応付けて、特定された複数の画素に対応する複数のＢ_１値各々を特定する。特定機能１３８により、造影Ｚスペクトルと非造影Ｚスペクトルとは、Ｂ_１値と対応付けられる。特定機能１３８を実現する処理回路１５０は、特定部に相当する。

処理回路１５０は、算出機能１４０により、複数の信号値の低下ポイント（または極小値）に対応する複数のＭＲ信号値を用いて、ｐＨ依存値を算出する。例えば、算出機能１４０は、複数のＭＲ信号値の低下ポイントのうち少なくとも２つのＭＲ信号値を用いて、ｐＨ依存値を算出する。このとき、ｐＨ依存値は、上記２つのＭＲ信号値の比を含む。ｐＨ依存値は、ｐＨに依存する値である。より詳細には、算出機能１４０は、レシオメトリック（Ｒａｔｉｏｍｅｔｒｉｃ）法を用いて、ｐＨ依存値を算出する。なお、ｐＨ依存値の算出方法は、レシオメトリック法に限定されず、既知の手法が適宜利用可能である。

具体的には、算出機能１４０は、造影Ｚスペクトルと非造影Ｚスペクトルとにより、複数の信号値の低下ポイントにおいて、複数種の物質に対応する複数の信号値の低下を算出する。例えば、造影剤がイオパミドールである場合、複数種の物質として２種のアミド基に関する周波数（４．２ｐｐｍと５．６ｐｐｍ）が、複数の信号値の低下に関する複数の周波数に該当する。すなわち、当該複数の周波数（４．２ｐｐｍと５．６ｐｐｍ）は、図２に示すように、イオパミドールに関する低下ポイントに対応する。このとき、算出機能１４０は、図２に示すように４．２ｐｐｍと５．６ｐｐｍとの位置における信号値の低下を、造影Ｚスペクトルと非造影ＺスペクトルとにおけるＭＲ信号値の差分（以下、低下信号値と呼ぶ）によりそれぞれ算出する。次いで、算出機能１４０は、４．２ｐｐｍに対応する低下信号値と５．６ｐｐｍに対応する低下信号値との比を用いて、例えば、レシオメトリック法により、ｐＨ依存値を算出する。算出機能１４０は、ｐＨ依存値の算出で用いられたＺスペクトルに関して特定されたＢ_１値を、当該算出されたｐＨ依存値と関連付ける。算出機能１４０を実現する処理回路１５０は、算出部に相当する。

処理回路１５０は、選択機能１４２により、ｐＨ依存値に対するｐＨ値を示す検量線を、被検体Ｐに関する情報に応じて予め用意された複数の検量線から、被検体Ｐに関する情報に応じて選択する。本実施形態においては、被検体Ｐに関する情報は、Ｂ_１マップである。このとき、選択機能１４２は、検量線を、Ｂ_１値に応じて予め用意された複数の検量線から、Ｂ_１マップに基づくＢ_１値（例えば、特定機能１３８により特定されたＢ_１値）に応じて選択する。具体的には、選択機能１４２は、ｐＨ依存値に関連するＢ_１値と同一のＢ_１値に関する検量線を、記憶回路１３２に記憶された複数の検量線から選択する。選択機能１４２による処理は、算出されたｐＨ依存値の数に対応するＢ_１値に応じて繰り返し選択される。選択機能１４２を実現する処理回路１５０は、選択部に相当する。

処理回路１５０は、ｐＨ値決定機能１４４により、選択された検量線とｐＨ依存値とに基づいて、ｐＨ値を決定する。具体的には、ｐＨ値決定機能１４４は、選択された検量線におけるｐＨ依存値と、算出されたｐＨ依存値とを照合することで、照合されたｐＨ依存値に対応するｐＨ値を決定する。ｐＨ値決定機能１４４による処理は、ｐＨ依存値の数に応じて繰り返し決定される。例えば、ｐＨ依存値の総数が、造影ＲＯＩまたは撮像領域全域に含まれる複数の画素の総数に対応する場合、ｐＨ値決定機能１４４は、Ｂ_１値に応じて選択された複数の検量線と複数のｐＨ依存値との照合により、複数のｐＨ依存値にそれぞれ対応する複数のｐＨ値を決定する。ｐＨ値決定機能１４４は、決定されたｐＨ値を、ｐＨ依存値に対応する画素の位置と関連づけて、記憶回路１３２に記憶させる。ｐＨ値決定機能１４４を実現する処理回路１５０は、ｐＨ値決定部に相当する。

以上、第１実施形態に係るＭＲＩ装置１００の全体構成について説明した。かかる構成のもと、第１実施形態に係るＭＲＩ装置１００は、非造影ＣＥＳＴイメージング、マップイメージング、造影ＣＥＳＴイメージングを実行し、各々のイメージングで収集されたＭＲデータと、Ｂ_１値を用いて決定された検量線とを用いてｐＨ値を決定する処理（以下、ｐＨ値決定処理と呼ぶ）実行する。以下、ｐＨ値決定処理に関する手順について、図４を参照して説明する。図４は、ｐＨ値決定処理の手順の一例を示すフローチャートである。

以下、説明を具体的にするために、静磁場強度は３Ｔであって、造影および非造影ＣＥＳＴイメージングにおいて複数の飽和パルスが印加される周波数の範囲は、－８ｐｐｍから８ｐｐｍの範囲であるものとする。このとき、０ｐｐｍに対応する飽和パルスの周波数、すなわち中心周波数は、１２８ＭＨｚである。また、非造影ＣＥＳＴイメージング、マップイメージング、および造影ＣＥＳＴイメージングに関するシーケンス情報は、予め設定されているものとする。上記シーケンス情報の設定としては、既知の手法が適用可能であるため、説明は省略する。また、ｐＨ値決定処理において、非造影ＣＥＳＴイメージング、マップイメージング、および造影ＣＥＳＴイメージングは、同一の対象を撮像対象としている。すなわち、非造影ＣＥＳＴイメージング、マップイメージング、および造影ＣＥＳＴイメージングにおける撮像対象部位は、同一である。

（ｐＨ値決定処理）
（ステップＳ４０１）
処理回路１５０は、インタフェース機能１３１により、非造影ＣＥＳＴイメージング、マップイメージング、および造影ＣＥＳＴイメージングに関するシーケンス情報を、シーケンス制御回路１２０に送信する。シーケンス制御回路１２０は、非造影ＣＥＳＴイメージングに関するシーケンス情報に従って、被検体Ｐに対して非造影ＣＥＳＴイメージングを実行する。処理回路１５０は、インタフェース機能１３１により、シーケンス制御回路１２０から非造影データを受信する。処理回路１５０は、画像生成機能１３４により、非造影データに基づいて、複数の非造影ＭＲ画像を生成する。画像生成機能１３４は、生成された複数の非造影ＭＲ画像を、記憶回路１３２に記憶させる。

（ステップＳ４０２）
シーケンス制御回路１２０は、マップイメージングに関するシーケンス情報に従って、被検体Ｐに対してマップイメージングを実行する。処理回路１５０は、インタフェース機能１３１により、シーケンス制御回路１２０からマップＭＲデータを受信する。処理回路１５０は、画像生成機能１３４により、マップＭＲデータに基づいて、Ｂ_０マップとＢ_１マップとを生成する。画像生成機能１３４は、生成されたＢ_０マップとＢ_１マップとを、記憶回路１３２に記憶させる。なお、本ステップの処理は、以下のステップＳ４０３の後に実行されてもよい。

（ステップＳ４０３）
本ステップに先立って、被検体Ｐに造影剤が注入される。このとき、注入される造影剤は、ケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤である。以下、説明を具体的にするために、造影剤は、イオパミドールであるものとする。なお、当該造影剤は、イオプロミドであってもよい。シーケンス制御回路１２０は、造影ＣＥＳＴイメージングに関するシーケンス情報に従って、被検体Ｐに対して造影ＣＥＳＴイメージングを実行する。処理回路１５０は、インタフェース機能１３１により、シーケンス制御回路１２０から造影データを受信する。処理回路１５０は、画像生成機能１３４により、造影データに基づいて複数の造影ＭＲ画像を生成する。画像生成機能１３４は、生成された複数の造影ＭＲ画像を、記憶回路１３２に記憶させる。

（ステップＳ４０４）
ディスプレイ１４３は、制御機能１３３による制御の元で、複数の非造影画像と複数の造影ＭＲ画像とのうち少なくとも一つを、表示する。入力装置１４１は、ユーザの指示により、ディスプレイ１４３に表示された画像に対して複数の非造影画像のうち少なくとも一つの非造影ＭＲ画像に対してＲＯＩを入力する。これにより、複数の非造影ＭＲ画像と複数の造影ＭＲ画像とのうち少なくとも一つに対してＲＯＩが設定される。なお、ｐＨ値の決定の対象となる画素が、撮像領域全域にわたる場合、本ステップにおける処理は不要となる。

（ステップＳ４０５）
処理回路１５０は、画像生成機能１３４により、複数の非造影ＭＲ画像とＢ_０マップとに基づいて、Ｂ_０補正を適用して、非造影Ｚスペクトルを生成する。具体的には、画像生成機能１３４は、設定されたＲＯＩに関して、非造影Ｚスペクトルを生成する。このとき、ＲＯＩに含まれる複数の画素各々に対してｐＨ値を決定すること（例えば、ｐＨマップの生成など）が入力装置１４１を介したユーザの指示により入力されていれば、画像生成機能１３４は、設定されたＲＯＩに含まれる複数の画素にそれぞれ対応する複数の非造影Ｚスペクトルを生成する。以下、説明を簡単にするために、設定されたＲＯＩに対して、一つの非造影Ｚスペクトルが生成されるものとする。画像生成機能１３４は、生成された非造影Ｚスペクトルを、記憶回路１３２に記憶させる。

（ステップＳ４０６）
処理回路１５０は、画像生成機能１３４により、複数の造影ＭＲ画像とＢ_０マップとに基づいて、Ｂ_０補正を適用して、複数種の物質に対応する複数の信号値の低下をそれぞれ含む造影Ｚスペクトルを生成する。具体的には、画像生成機能１３４は、設定されたＲＯＩに関して、造影Ｚスペクトルを生成する。このとき、ＲＯＩに含まれる複数の画素各々に対してｐＨ値を決定すること（例えば、ｐＨマップの生成など）が入力装置１４１を介したユーザの指示により入力されていれば、画像生成機能１３４は、設定されたＲＯＩに含まれる複数の画素にそれぞれ対応する複数の造影Ｚスペクトルを生成する。設定されたＲＯＩに対して一つの非造影Ｚスペクトルが生成され場合、画像生成機能１３４は、設定されたＲＯＩに対して一つの造影Ｚスペクトルを生成する。画像生成機能１３４は、生成された造影Ｚスペクトルを、記憶回路１３２に記憶させる。なお、ステップＳ４０５の処理とステップＳ４０６の処理との順番は逆であってもよい。

（ステップＳ４０７）
処理回路１５０は、取得機能１３６により、記憶回路１３２から、非造影Ｚスペクトルと、造影Ｚスペクトルと、被検体Ｐに関する情報に対応するＢ_１マップとを取得する。なお、処理回路１５０がＭＲＩ装置１００とは別個のｐＨ算出装置などに搭載される場合、取得機能１３６は、例えば他のＭＲＩ装置から、非造影Ｚスペクトルと、造影Ｚスペクトルと、Ｂ_１マップとを取得する。

（ステップＳ４０８）
処理回路１５０は、特定機能１３８により、非造影Ｚスペクトルおよび造影Ｚスペクトルに関して、Ｂ_１マップ上での位置でのＢ_１値を特定する。ＲＯＩに対して１つの非造影Ｚスペクトルおよび１つの造影Ｚスペクトルが生成された場合、特定機能１３８は、Ｂ_１マップ上でのＲＯＩに対応する領域に含まれる複数のＢ_１値の和を、ＲＯＩに含まれる画素数で除算することにより、非造影Ｚスペクトルおよび造影Ｚスペクトルに対応するＢ_１値を特定する。すなわち、１つの非造影Ｚスペクトルの生成および１つの造影Ｚスペクトルの生成において複数の画素が関連している場合、特定機能１３８は、当該複数の画素の位置に対応する複数のＢ_１値の平均を演算することで、１つの非造影Ｚスペクトルおよび１つの造影Ｚスペクトルに対応するＢ_１値を特定する。なお、本ステップにおける処理は、後述のステップＳ４０９の後に実行されてもよい。

なお、複数の画素に応じて複数の非造影Ｚスペクトルおよび複数の造影Ｚスペクトルが生成された場合、特定機能１３８は、当該複数の画素各々に対応するＢ_１マップ上での位置でのＢ_１値を、特定する。例えば、撮像領域の全域に亘る複数の画素ごとに非造影Ｚスペクトルおよび造影Ｚスペクトルが生成された場合、特定機能１３８は、非造影Ｚスペクトルの生成および造影Ｚスペクトルの生成に関する画素の位置を用いて、非造影Ｚスペクトルおよび造影Ｚスペクトルに対応するＢ_１値を、Ｂ_１マップから特定する。

（ステップＳ４０９）
処理回路１５０は、算出機能１４０により、造影Ｚスペクトルと非造影Ｚスペクトルとに基づいて、ｐＨ依存値を算出する。具体的には、算出機能１４０は、造影スペクトルにおける複数の信号値の低下ポイントに対応する複数のＭＲ信号値と、非造影Ｚスペクトルにおける当該複数の周波数に対応する複数のＭＲ信号値とを用いて、４．２ｐｐｍと５．６ｐｐｍとにおける２つの低下信号値を算出する。次いで、算出機能１４０は、２つの低下信号値の比の演算を含む所定の演算処理により、ｐＨ依存値を算出する。所定の演算処理は、例えば、レシオメトリック（Ｒａｔｉｏｍｅｔｒｉｃ）法における演算の処理である。

（ステップＳ４１０）
処理回路１５０は、選択機能１４２により、記憶回路１３２に記憶された複数の検量線から、算出されたｐＨ依存値に対するｐＨ値を示す検量線を、特定されたＢ_１値に応じて選択する。検量線の選択は、例えば、特定されたＢ_１値の数に応じて複数回実行される。なお、特定されたＢ_１値に対応する検量線が記憶回路１３２に記憶されていない場合、選択機能１４２は、特定されたＢ_１値に最も近いＢ_１値に対応する検量線を、複数の検量線から選択してもよい。

また、特定されたＢ_１値に対応する検量線が記憶回路１３２に記憶されていない場合、選択機能１４２は、特定されたＢ_１値に近い２つのＢ_１値に対応する２つの検量線を、複数の検量線から選択してもよい。このとき、算出機能１４０は、選択された２つの検量線に関する２つのＢ_１値と特定されたＢ_１値との差に基づいて、選択された２つの検量線から、特定されたＢ_１値に対応する検量線を、例えば重みづけ加算、平均などの補間処理により計算する。これらにより、選択機能１４２は、計算された検量線を、特定されたＢ_１値に対応する検量線として、複数の検量線から選択する。

（ステップＳ４１１）
処理回路１５０は、ｐＨ値決定機能１４４により、選択された検量線とｐＨ依存値とに基づいて、ｐＨ値を決定する。具体的には、ｐＨ値決定機能１４４は、選択された検量線において、算出されたｐＨ依存値に対応するｐＨ値を、ｐＨ依存値の算出に用いられた造影Ｚスペクトルに関する画素の位置におけるｐＨ値として決定する。算出されたｐＨ依存値が複数である場合、ｐＨ値決定機能１４４は、ｐＨ依存値の決定を、造影Ｚスペクトルの数に対応する回数に亘って繰り返し実行する。

（ステップＳ４１２）
処理回路１５０は、制御機能１３３による制御の元で、決定されたｐＨ値をディスプレイ１４３に表示する。例えば、処理回路１５０は、画像生成機能１３４により、１つの造影ＭＲ画像または１つの非造影画像などにおけるＭＲ画像において、ｐＨ依存値が算出された画素の位置に、決定されたｐＨ値を重畳した重畳画像を生成する。重畳画像におけるｐＨ値は、例えば、ｐＨの値に応じたカラースケールまたはグレースケールで表される。処理回路１５０は、制御機能１３３により、重畳画像をディスプレイ１４３に表示する。

また、複数の画素にそれぞれ対応する複数のｐＨ値が決定された場合、処理回路１５０は、画像生成機能１３４により、複数のｐＨ値に対応する複数の画素の位置に当該複数のｐＨ値を配列することで、ｐＨマップを生成してもよい。ｐＨマップにおけるｐＨ値は、例えば、ｐＨの値に応じたカラースケールまたはグレースケールで表される。また、画像生成機能１３４は、上記ＭＲ画像にｐＨマップを重畳してもよい。処理回路１５０は、制御機能１３３により、ｐＨマップが重畳された重畳画像および／またはｐＨマップを、ディスプレイ１４３に表示する。

以上に述べた第１実施形態に係るＭＲＩ装置１００は、被検体Ｐに対しケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ撮像と、ＣＥＳＴ撮像とは異なる撮像とを実行し、異なる撮像により収集されたＭＲデータに基づいてＢ_１マップを生成し、ＣＥＳＴ撮像により収集されたＭＲデータに基づいてＺスペクトルを生成し、複数の信号値の低下ポイントに対応する複数のＭＲ信号値を用いてｐＨ依存値を算出し、ｐＨ依存値に対するｐＨ値を示す検量線を、Ｂ_１値に応じて予め用意された複数の検量線から、Ｂ_１マップに基づくＢ_１値に応じて選択し、選択された検量線とｐＨ依存値とに基づいて、ｐＨ値を決定する。

また、第１実施形態に係るＭＲＩ装置１００は、複数のＭＲ信号値のうち少なくとも２つのＭＲ信号値を用いてｐＨ依存値を算出し、ｐＨ依存値は、２つのＭＲ信号値の比を含む。また、第１実施形態に係るＭＲＩ装置１００は、レシオメトリック（Ｒａｔｉｏｍｅｔｒｉｃ）法を用いてｐＨ依存値を算出する。また、第１実施形態に係るＭＲＩ装置１００におけるＺスペクトルは、異なる撮像により生成されたＢ_０マップに基づいてＣＥＳＴ撮像における飽和パルスの位置が補正されたスペクトルデータである。

これらのことから、第１実施形態に係るＭＲＩ装置１００によれば、複数のＢ_１値に対応する複数の検量線を予め記憶回路１３２に記憶することで、Ｂ_１マップにおいて特定されたＢ_１値に対応する検量線を用いて、ｐＨ依存値からｐＨ値を決定することができる。このため、第１実施形態に係るＭＲＩ装置１００によれば、Ｂ_１不均一性によるｐＨ値のズレを小さく補正することができる。

以上のことから、第１実施形態に係るＭＲＩ装置１００によれば、Ｂ_１不均一性によらずに、より正確なｐＨ値を決定することができるため、被検体Ｐに対する診断の正確性を向上させることができる。

（第２実施形態）
本実施形態は、ｐＨ値決定処理におけるステップＳ４０７以降の処理を、ｐＨ算出装置で実現することにある。図５は、ｐＨ算出装置１の構成の一例を示す図である。第１実施形態との相違は、制御機能１３３における処理と取得機能１３６における処理とにある。制御機能１３３と取得機能１３６とを除く他の構成は、ＭＲＩ装置３特有の処理を除き第１実施形態と同様なため、説明は省略する。制御機能１３３は、ｐＨ算出装置１の全体を制御する。また、取得機能１３６は、例えば他のＭＲＩ装置３から、非造影Ｚスペクトルと、造影Ｚスペクトルと、被検体Ｐに関する情報とを取得する。被検体Ｐに関する情報は、例えば、Ｂ_１マップである。また、本実施形態により実現されるｐＨ値決定処理は、ステップＳ４０７以降の処理が実行されることにある。ｐＨ算出装置１によるステップＳ４０７以降の処理は、第１実施形態と同様なため、説明は算出する。

第２実施形態に係るｐＨ算出装置１は、被検体Ｐに対しケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ撮像により生成され複数種の物質に対応する複数の信号値の低下ポイントを含む少なくとも一つのＺスペクトルと、ＣＥＳＴ撮像とは異なる撮像により生成されたＢ_１マップとを取得し、複数の信号値の低下ポイントに対応する複数のＭＲ信号値を用いてｐＨ依存値を算出し、ｐＨ依存値に対するｐＨ値を示す検量線を、Ｂ_１値に応じて予め用意された複数の検量線から、前記Ｂ_１マップに基づくＢ_１値に応じて選択し、選択された検量線とｐＨ依存値とに基づいて、ｐＨ値を決定する。ｐＨ算出装置１の各種構成要素のおける処理およびｐＨ算出装置１により実現されるｐＨ値決定処理による効果等は、第１実施形態と同様なため、説明は省略する。

（第３実施形態）
本実施形態は、被検体Ｐに関する情報として、Ｂ_１マップと被検体Ｐの温度とを用いることにある。第１実施形態および第２実施形態では、被検体Ｐに関する情報として、Ｂ_１マップを用いて説明した。第１実施形態および第２実施形態では、被検体Ｐの温度は、例えば、概ね３７℃と仮定している。一方、被検体Ｐの温度（体温）は、個体差がある。加えて、同一の被検体Ｐであっても被検体Ｐの体調により、体温が異なることがある。また、ＲＯＩ内に癌が存在する場合、癌の温度は、被検体Ｐの体温と異なること（例えば、体温より低温または体温より高温）がある。ＣＥＳＴ効果におけるプロトンの交換速度は、温度とｐＨとに依存する。

図６は、ファントムの温度が２５℃である場合の３つのｐＨ値（ｐＨ７、ｐＨ７．２、ｐＨ７．４）に応じた３つのＺスペクトルの一例を示す図である。図７は、ファントムの温度が３６℃である場合の３つのｐＨ値（ｐＨ６．４、ｐＨ６．６、ｐＨ６．８）に応じた３つのＺスペクトルの一例を示す図である。図６および図７における複数のＺスペクトルは、温度およびｐＨを除いて、略同一の撮像条件によりＭＲ信号に基づいて、プロットされている。

図６および図７に示すように、図６におけるｐＨ７、ｐＨ７．２、ｐＨ７．４に対応するＺスペクトルの形状は、図７におけるｐＨ６．４、ｐＨ６．６、ｐＨ６．８に対応するＺスペクトルの形状と酷似している。Ｚスペクトルの形状の酷似は、水素原子の化学交換速度の酷似に対応する。例えば、図６および図７に示すように、２５℃のときのｐＨ７の時の水素原子の化学交換速度と、３６℃のときのｐＨ６．４の時の水素原子の化学交換速度とは酷似している。このため、例えば、図６における２５℃、ｐＨ７に対応するＺスペクトルがＣＥＳＴ撮像により取得された場合、被検体Ｐの温度が３６℃であればｐＨ値は６．４であるにもかかわらず、２５℃の検量線を用いてｐＨ値を決定するとｐＨ値は７となる。

これらのことから、被検体Ｐの温度（体温）を一定と仮定した第１実施形態および第２実施形態では、ｐＨ値の精度が低下する可能性がある。このため、本実施形態では、被検体Ｐに関する情報として、Ｂ_１マップに加えて被検体Ｐの温度を用いる。すなわち、本実施形態では、Ｂ_１マップと被検体Ｐの温度とを取得し、被検体Ｐの温度とＢ_１値とに応じた検量線を用いてｐＨ値を決定する。

記憶回路１３２は、複数のＢ_１値および複数の温度に応じた複数の検量線を記憶する。記憶回路１３２は、例えば、複数の対応表（ＬＵＴ：Ｌｏｏｋ Ｕｐ Ｔａｂｌｅ）の形式で、複数の検量線を、複数のＢ_１値と複数の温度と対応付けて記憶する。複数の検量線は、例えば、実施形態と同様にファントムなどにより予め計測された実測値により予め用意されて、記憶回路１３２に記憶される。本実施形態における実測値は、例えば、Ｂ_１値、温度、ｐＨ依存値、ｐＨ値などである。複数の検量線は、実測に基づくデータ、シミュレーションで計算されたデータ、または当該実測値から補間されたデータで、予め用意されたデータであってもよい。本実施形態における補間は、隣接する２つのＢ_１値の間および隣接する２つの温度の間を補間して検量線を算出することにある。複数の検量線は、複数のＢ_１値および複数の温度により分類された検量線のデータベースに相当する。

処理回路１５０は、取得機能１３６により、少なくとも一つのＺスペクトルと、被検体Ｐに関する情報として被検体Ｐの温度とＢ_１マップとを取得する。被検体Ｐの温度は、ｐＨ値を決定する処理の実行中にモニタリングされてもよいし、当該処理の実行前に取得されてもよい。

処理回路１５０は、選択機能１４２により、ｐＨ依存値に対するｐＨ値を示す検量線を、被検体Ｐに関する情報に応じて予め用意された複数の検量線から、被検体Ｐに関する情報に応じて選択する。本実施形態においては、被検体Ｐに関する情報は、Ｂ_１マップと被検体Ｐの温度である。このとき、選択機能１４２は、検量線を、Ｂ_１値と被検体Ｐの温度とに応じて予め用意された複数の検量線から、Ｂ_１マップに基づくＢ_１値（例えば、特定機能１３８により特定されたＢ_１値）と被検体Ｐの温度とに応じて選択する。具体的には、選択機能１４２は、ｐＨ依存値に関連するＢ_１値と同一のＢ_１値に関する検量線であって、被検体Ｐの温度と同一の温度に対応する検量線を、記憶回路１３２に記憶された複数の検量線から選択する。

本実施形態に係るＭＲＩ装置１００は、非造影ＣＥＳＴイメージング、マップイメージング、造影ＣＥＳＴイメージングを実行し、各々のイメージングで収集されたＭＲデータとＢ_１値と被検体Ｐの温度とを用いて決定された検量線とを用いてｐＨ値を決定するｐＨ値決定処理を実行する。以下、本実施形態に係るｐＨ値決定処理に関する手順について、図８を用いて説明する。

図８は、本実施形態に係るｐＨ値決定処理の手順の一例を示すフローチャートである。図８に示すフローチャートは、第１実施形態および第２実施形態におけるｐＨ値決定処理のステップＳ４０６に続く処理の手順に相当する。すなわち、図８におけるステップＳ８０１の処理は、ステップＳ４０６の後に実行される。

（ｐＨ値決定処理）
（ステップＳ８０１）
被検体Ｐの温度が計測される。被検体Ｐの温度は、例えば体温である。被検体Ｐの温度は、例えば、既知の温度計、ＣＥＳＴイメージングまたはマップイメージングにより収集された位相画像における位相差に基づく既知の温度計測手法などにより、計測される。被検体Ｐの温度の計測は、ｐＨ値決定処理におけるＣＥＳＴイメージングの実施前に実行されてもよいし、ＣＥＳＴイメージング中に実施されてもよい。

（ステップＳ８０２）
処理回路１５０は、取得機能１３６により、記憶回路１３２から、非造影Ｚスペクトルと、造影Ｚスペクトルと、被検体Ｐに関する情報に含まれるＢ_１マップおよび被検体Ｐの温度とを取得する。なお、処理回路１５０がＭＲＩ装置１００とは別個のｐＨ算出装置１などに搭載される場合、取得機能１３６は、例えば他のＭＲＩ装置から、非造影Ｚスペクトルと、造影Ｚスペクトルと、Ｂ_１マップと、被検体Ｐの温度とを取得する。なお、被検体Ｐの温度は、他のＭＲＩ装置の代わりに既知の温度計から取得されてもよい。

（ステップＳ８０３）
処理回路１５０は、特定機能１３８により、非造影Ｚスペクトルおよび造影Ｚスペクトルに関して、Ｂ_１マップ上での位置でのＢ_１値を特定する。本ステップの処理は、ステップＳ４０８と同様なため、説明は省略する。

（ステップＳ８０４）
処理回路１５０は、算出機能１４０により、造影スペクトルにと非造影Ｚスペクトルとに基づいて、ｐＨ依存値を算出する。本ステップの処理は、ステップＳ４０９と同様なため、説明は省略する。

（ステップＳ８０５）
処理回路１５０は、選択機能１４２により、記憶回路１３２に記憶された複数の検量線から、算出されたｐＨ依存値に対するｐＨ値を示す検量線を、特定されたＢ_１値および取得された被検体Ｐの体温に応じて選択する。検量線の選択は、例えば、特定されたＢ_１値の数および取得された被検体Ｐの体温の数に応じて複数回実行される。なお、特定されたＢ_１値および取得された被検体Ｐの体温に対応する検量線が記憶回路１３２に記憶されていない場合、選択機能１４２は、特定されたＢ_１値に最も近いＢ_１値および取得された被検体Ｐの体温に最も近い温度に対応する検量線を、複数の検量線から選択してもよい。

また、特定されたＢ_１値および取得された被検体Ｐの温度に対応する検量線が記憶回路１３２に記憶されていない場合、選択機能１４２は、特定されたＢ_１値に近い２つのＢ_１値に対応する２つの検量線と、取得された被検体Ｐの温度に近い２つの温度に対応する２つの検量線とを、複数の検量線から選択してもよい。このとき、算出機能１４０は、例えば、選択された４つの検量線に関する４つのＢ_１値と特定されたＢ_１値との差および４つの検量線に関する４つの温度と取得された温度との差に基づいて、選択された４つの検量線から、特定されたＢ_１値および取得された温度に対応する検量線を、例えば重みづけ加算、平均などの補間処理により計算する。これらにより、選択機能１４２は、計算された検量線を、特定されたＢ_１値および取得された温度に対応する検量線に対応する検量線として、複数の検量線から選択してもよい。

（ステップＳ８０６）
処理回路１５０は、ｐＨ値決定機能１４４により、選択された検量線とｐＨ依存値とに基づいて、ｐＨ値を決定する。本ステップの処理は、ステップＳ４１１と同様なため、説明は省略する。

（ステップＳ８０７）
処理回路１５０は、制御機能１３３による制御の元で、決定されたｐＨ値をディスプレイ１４３に表示する。本ステップの処理は、ステップＳ４１２と同様なため、説明は省略する。

以上に述べた第３実施形態に係るＭＲＩ装置１００は、被検体Ｐに対しケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ撮像と、ＣＥＳＴ撮像とは異なる撮像とを実行し、異なる撮像により収集されたＭＲデータに基づいてＢ_１マップを生成し、被検体Ｐの温度を取得し、ＣＥＳＴ撮像により収集されたＭＲデータに基づいてＺスペクトルを生成し、複数の信号値の低下ポイントに対応する複数のＭＲ信号値を用いてｐＨ依存値を算出し、ｐＨ依存値に対するｐＨ値を示す検量線を、Ｂ_１値および温度に応じて予め用意された複数の検量線から、Ｂ_１マップに基づくＢ_１値および被検体Ｐの温度に応じて選択し、選択された検量線とｐＨ依存値とに基づいて、ｐＨ値を決定する。

第３実施形態に係るＭＲＩ装置１００によれば、複数のＢ_１値および複数の温度に対応する複数の検量線を予め記憶回路１３２に記憶することで、Ｂ_１マップにおいて特定されたＢ_１値および被検体Ｐから取得された被検体Ｐの体温に対応する検量線を用いて、ｐＨ依存値からｐＨ値を決定することができる。このため、第３実施形態に係るＭＲＩ装置１００によれば、Ｂ_１不均一性によるｐＨ値のズレを小さく補正し、かつ温度により生じるｐＨ値のズレを補正することができる。

以上のことから、第３実施形態に係るＭＲＩ装置１００によれば、Ｂ_１不均一性および被検体Ｐの温度によらずに、より正確なｐＨ値を決定することができる。このため、第３実施形態に係るＭＲＩ装置１００によれば、被検体Ｐの体調および癌の有無に依らず、ｐＨの測定精度をより向上させ、被検体Ｐに対する診断の正確性を向上させることができる。

なお、本実施形態の変形例として、本実施形態における技術的特徴は、ｐＨ算出装置１により実現されてもよい。ｐＨ算出装置１の各種構成要素のおける処理およびｐＨ算出装置１により実現されるｐＨ値決定処理による効果等は、第３実施形態と同様なため、説明は省略する。

（第４実施形態）
本実施形態は、被検体Ｐに関する情報として、被検体Ｐの温度を用いることにある。記憶回路１３２は、複数の温度に応じた複数の検量線を記憶する。記憶回路１３２は、例えば、複数の対応表（ＬＵＴ：Ｌｏｏｋ Ｕｐ Ｔａｂｌｅ）の形式で、複数の検量線を、複数の温度と対応付けて記憶する。複数の検量線は、例えば、実施形態と同様にファントムなどにより予め計測された実測値により予め用意されて、記憶回路１３２に記憶される。本実施形態における実測値は、例えば、温度、ｐＨ依存値、ｐＨ値などである。複数の検量線は、実測に基づくデータ、シミュレーションで計算されたデータ、または当該実測値から補間されたデータで、予め用意されたデータであってもよい。本実施形態における補間は、隣接する２つの温度の間を補間して検量線を算出することにある。複数の検量線は、複数の温度により分類された検量線のデータベースに相当する。

処理回路１５０は、取得機能１３６により、少なくとも一つのＺスペクトルと、被検体Ｐに関する情報として被検体Ｐの温度を取得する。被検体Ｐの温度は、ｐＨ値を決定する処理の実行中にモニタリングされてもよいし、当該処理の実行前に取得されてもよい。

処理回路１５０は、選択機能１４２により、ｐＨ依存値に対するｐＨ値を示す検量線を、被検体Ｐに関する情報に応じて予め用意された複数の検量線から、被検体Ｐの温度に応じて選択する。このとき、選択機能１４２は、検量線を、被検体Ｐの温度に応じて予め用意された複数の検量線から、被検体Ｐの温度に応じて選択する。具体的には、選択機能１４２は、被検体Ｐの温度と同一の温度に対応する検量線を、記憶回路１３２に記憶された複数の検量線から選択する。

本実施形態に係るＭＲＩ装置１００は、非造影ＣＥＳＴイメージング、マップイメージング、造影ＣＥＳＴイメージングを実行し、各々のイメージングで収集されたＭＲデータと被検体Ｐの温度を用いて決定された検量線とを用いてｐＨ値を決定するｐＨ値決定処理を実行する。以下、本実施形態に係るｐＨ値決定処理に関する手順について、図９を用いて説明する。

図９は、本実施形態に係るｐＨ値決定処理の手順の一例を示すフローチャートである。図９に示すフローチャートは、第１実施形態および第２実施形態におけるｐＨ値決定処理のステップＳ４０６に続く処理の手順に相当する。すなわち、図９におけるステップＳ９０１の処理は、ステップＳ４０６の後に実行される。なお、本実施形態において、ステップＳ４０２におけるＢ_１マップの生成は不要である。

（ｐＨ値決定処理）
（ステップＳ９０１）
被検体Ｐの温度が計測される。本ステップにおける処理は、ステップＳ８０１と同様なため、説明は省略する。

（ステップＳ９０２）
処理回路１５０は、取得機能１３６により、記憶回路１３２から、非造影Ｚスペクトルと、造影Ｚスペクトルと、被検体Ｐに関する情報としての被検体Ｐの温度を取得する。なお、処理回路１５０がＭＲＩ装置１００とは別個のｐＨ算出装置１などに搭載される場合、取得機能１３６は、例えば他のＭＲＩ装置から、非造影Ｚスペクトルと、造影Ｚスペクトルと、被検体Ｐの温度とを取得する。なお、被検体Ｐの温度は、他のＭＲＩ装置の代わりに既知の温度計から取得されてもよい。

（ステップＳ９０３）
処理回路１５０は、算出機能１４０により、造影スペクトルにと非造影Ｚスペクトルとに基づいて、ｐＨ依存値を算出する。本ステップの処理は、ステップＳ４０９と同様なため、説明は省略する。

（ステップＳ９０４）
処理回路１５０は、選択機能１４２により、記憶回路１３２に記憶された複数の検量線から、算出されたｐＨ依存値に対するｐＨ値を示す検量線を、取得された被検体Ｐの体温に応じて選択する。検量線の選択は、例えば、取得された被検体Ｐの体温の数に応じて複数回実行される。なお、取得された被検体Ｐの体温に対応する検量線が記憶回路１３２に記憶されていない場合、選択機能１４２は、取得された被検体Ｐの体温に最も近い温度に対応する検量線を、複数の検量線から選択してもよい。

なお、取得された被検体Ｐの温度に対応する検量線が記憶回路１３２に記憶されていない場合、選択機能１４２は、取得された被検体Ｐの温度に近い２つの温度に対応する２つの検量線を、複数の検量線から選択してもよい。このとき、算出機能１４０は、例えば、選択された２つの検量線に関する２つの温度と取得された温度との差に基づいて、選択された２つの検量線から、取得された温度に対応する検量線を、例えば重みづけ加算、平均などの補間処理により計算する。これらにより、選択機能１４２は、計算された検量線を、取得された温度に対応する検量線に対応する検量線として、複数の検量線から選択してもよい。

（ステップＳ９０５）
処理回路１５０は、ｐＨ値決定機能１４４により、選択された検量線とｐＨ依存値とに基づいて、ｐＨ値を決定する。本ステップの処理は、ステップＳ４１１と同様なため、説明は省略する。

（ステップＳ９０６）
処理回路１５０は、制御機能１３３による制御の元で、決定されたｐＨ値をディスプレイ１４３に表示する。本ステップの処理は、ステップＳ４１２と同様なため、説明は省略する。

以上に述べた第４実施形態に係るＭＲＩ装置１００は、被検体Ｐに対しケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ撮像を実行し、被検体Ｐの温度を取得し、ＣＥＳＴ撮像により収集されたＭＲデータに基づいてＺスペクトルを生成し、複数の信号値の低下ポイントに対応する複数のＭＲ信号値を用いてｐＨ依存値を算出し、ｐＨ依存値に対するｐＨ値を示す検量線を、温度に応じて予め用意された複数の検量線から、被検体Ｐの温度に応じて選択し、選択された検量線とｐＨ依存値とに基づいて、ｐＨ値を決定する。

第４実施形態に係るＭＲＩ装置１００によれば、複数の温度に対応する複数の検量線を予め記憶回路１３２に記憶することで、被検体Ｐから取得された被検体Ｐの体温に対応する検量線を用いて、ｐＨ依存値からｐＨ値を決定することができる。このため、第４実施形態に係るＭＲＩ装置１００によれば、温度により生じるｐＨ値のズレを補正することができる。

以上のことから、第４実施形態に係るＭＲＩ装置１００によれば、Ｂ_１不均一性が極めて小さい場合であっても、被検体Ｐの温度によらずに、正確なｐＨ値を決定することができる。このため、第４実施形態に係るＭＲＩ装置１００によれば、被検体Ｐの体調および癌の有無に依らず、ｐＨの測定精度を向上させ、被検体Ｐに対する診断の正確性を向上させることができる。

なお、本実施形態の変形例として、本実施形態における技術的特徴は、ｐＨ算出装置１により実現されてもよい。ｐＨ算出装置１の各種構成要素のおける処理およびｐＨ算出装置１により実現されるｐＨ値決定処理による効果等は、第４実施形態と同様なため、説明は省略する。

実施形態における技術的思想をｐＨ算出方法で実現する場合、当該ｐＨ算出方法は、被検体Ｐに対しケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ（Ｃｈｅｍｉｃａｌ Ｅｘｃｈａｎｇｅ Ｓａｔｕｒａｔｉｏｎ Ｔｒａｎｓｆｅｒ）撮像により生成され複数種の物質にそれぞれ対応する複数の信号値の低下ポイントを含む少なくとも一つのＺスペクトルと、被検体Ｐに関する情報とを取得し、複数の信号値の低下ポイントに対応する複数のＭＲ信号値を用いてｐＨ依存値を算出し、ｐＨ依存値に対するｐＨ値を示す検量線を、被検体Ｐに関する情報に応じて予め用意された複数の検量線から、被検体Ｐに関する情報に応じて選択し、選択された検量線とｐＨ依存値とに基づいて、ｐＨ値を決定する。ｐＨ算出方法により実行されるｐＨ値決定処理の手順および効果は、第１乃至第４実施形態と同様なため、説明は省略する。

以上説明した少なくとも一つの実施形態等によれば、算出されるｐＨ値の正確性（精度）を向上することができる。

いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

以上の実施形態等に関し、発明の一側面および選択的な特徴として以下の付記を開示する。
（付記１）
被検体に対しｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ（Ｃｈｅｍｉｃａｌ Ｅｘｃｈａｎｇｅ Ｓａｔｕｒａｔｉｏｎ Ｔｒａｎｓｆｅｒ）撮像により生成され前記複数種の物質にそれぞれ対応する複数の信号値の低下ポイントまたは極小値を含む少なくとも一つのＺスペクトルと、前記被検体に関する情報とを取得する取得部と、
前記複数の信号値の低下ポイントまたは極小値に対応する複数のＭＲ信号値を用いて、ｐＨ依存値を算出する算出部と、
前記ｐＨ依存値に対するｐＨ値を示す検量線を、前記被検体に関する情報に応じて予め用意された複数の検量線から、前記被検体に関する情報に応じて選択する選択部と、
前記選択された検量線と前記ｐＨ依存値とに基づいて、ｐＨ値を決定するｐＨ値決定部と、
を備えるｐＨ算出装置。
（付記２）
前記被検体に関する情報は、前記ＣＥＳＴ撮像とは異なる撮像により生成されたＢ_１マップと前記被検体の温度とのうち少なくとも一つであってもよい。
（付記３）
前記算出部は、前記複数のＭＲ信号値のうち少なくとも２つのＭＲ信号値を用いて、前記ｐＨ依存値を算出してもよく、
前記ｐＨ依存値は、前記２つのＭＲ信号値の比を含んでもよい。
（付記４）
前記算出部は、レシオメトリック（Ｒａｔｉｏｍｅｔｒｉｃ）法を用いて前記ｐＨ依存値を算出してもよい。
（付記５）
前記Ｚスペクトルは、前記ＣＥＳＴ撮像とは異なる撮像により生成されたＢ_０マップに基づいて前記ＣＥＳＴ撮像における飽和パルスの位置が補正されたスペクトルデータであってもよい。
（付記６）
前記造影剤は、イオパミドールまたはイオプロミドであってもよい。
（付記７）
前記複数の検量線は、シミュレーションで計算されたデータ、または実測値から補間されたデータであってもよい。
（付記８）
付記１に記載のｐＨ算出装置と、
前記ＣＥＳＴ撮像を実行する制御部と、
を備える磁気共鳴イメージング装置。
（付記９）
被検体に対しｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ（Ｃｈｅｍｉｃａｌ Ｅｘｃｈａｎｇｅ Ｓａｔｕｒａｔｉｏｎ Ｔｒａｎｓｆｅｒ）撮像により生成され前記複数種の物質にそれぞれ対応する複数の信号値の低下ポイントまたは極小値を含む少なくとも一つのＺスペクトルと、前記被検体に関する情報とを取得し、
前記複数の信号値の低下ポイントまたは極小値に対応する複数のＭＲ信号値を用いて、ｐＨ依存値を算出し、
前記ｐＨ依存値に対するｐＨ値を示す検量線を、前記被検体に関する情報に応じて予め用意された複数の検量線から、前記被検体に関する情報に応じて選択し、
前記選択された検量線と前記ｐＨ依存値とに基づいて、ｐＨ値を決定すること、
を備えるｐＨ算出方法。

１ ｐＨ算出装置
１００ 磁気共鳴イメージング装置
１０１ 静磁場磁石
１０２ 静磁場電源
１０３ 傾斜磁場コイル
１０４ 傾斜磁場電源
１０５ 寝台
１０５ａ 天板
１０６ 寝台制御回路
１０７ 送信コイル
１０８ 送信回路
１０９ 受信コイル
１２０ シーケンス制御回路
１３０ コンピューター
１３１ インタフェース機能
１３２ 記憶回路
１３３ 制御機能
１３４ 画像生成機能
１３６ 取得機能
１３８ 特定機能
１４０ 算出機能
１４１ 入力装置
１４２ 選択機能
１４３ ディスプレイ
１４４ ｐＨ値決定機能
１５０ 処理回路

Claims (9)

1. 被検体に対しケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ（Ｃｈｅｍｉｃａｌ Ｅｘｃｈａｎｇｅ Ｓａｔｕｒａｔｉｏｎ Ｔｒａｎｓｆｅｒ）撮像により生成され前記複数種の物質にそれぞれ対応する複数の信号値の低下ポイントを含む少なくとも一つのＺスペクトルと、前記被検体に関する情報とを取得する取得部と、
   前記複数の信号値の低下ポイントに対応する複数のＭＲ信号値を用いて、ｐＨ依存値を算出する算出部と、
   前記ｐＨ依存値に対するｐＨ値を示す検量線を、前記被検体に関する情報に応じて予め用意された複数の検量線から、前記被検体に関する情報に応じて選択する選択部と、
   前記選択された検量線と前記ｐＨ依存値とに基づいて、ｐＨ値を決定するｐＨ値決定部と、
   を備えるｐＨ算出装置。
2. 前記被検体に関する情報は、前記ＣＥＳＴ撮像とは異なる撮像により生成されたＢ_１マップと前記被検体の温度とのうち少なくとも一つである、
   請求項１に記載のｐＨ算出装置。
3. 前記算出部は、前記複数のＭＲ信号値のうち少なくとも２つのＭＲ信号値を用いて、前記ｐＨ依存値を算出し、
   前記ｐＨ依存値は、前記２つのＭＲ信号値の比を含む、
   請求項１に記載のｐＨ算出装置。
4. 前記算出部は、レシオメトリック（Ｒａｔｉｏｍｅｔｒｉｃ）法を用いて前記ｐＨ依存値を算出する、
   請求項３に記載のｐＨ算出装置。
5. 前記Ｚスペクトルは、前記ＣＥＳＴ撮像とは異なる撮像により生成されたＢ_０マップに基づいて前記ＣＥＳＴ撮像における飽和パルスの位置が補正されたスペクトルデータである、
   請求項１乃至４のいずれか一項に記載のｐＨ算出装置。
6. 前記造影剤は、イオパミドールまたはイオプロミドである、
   請求項１に記載のｐＨ算出装置。
7. 前記複数の検量線は、シミュレーションで計算されたデータ、または実測値から補間されたデータである、
   請求項１に記載のｐＨ算出装置。
8. 請求項１に記載のｐＨ算出装置と、
   前記ＣＥＳＴ撮像を実行するシーケンス制御部と、
   を備える磁気共鳴イメージング装置。
9. 被検体に対しケミカルシフトイメージングにおけるｐＨに対する応答性が互いに異なる複数種の物質を含む造影剤を用いたＣＥＳＴ（Ｃｈｅｍｉｃａｌ Ｅｘｃｈａｎｇｅ Ｓａｔｕｒａｔｉｏｎ Ｔｒａｎｓｆｅｒ）撮像により生成され前記複数種の物質にそれぞれ対応する複数の信号値の低下ポイントを含む少なくとも一つのＺスペクトルと、前記被検体に関する情報とを取得し、
   前記複数の信号値の低下ポイントに対応する複数のＭＲ信号値を用いて、ｐＨ依存値を算出し、
   前記ｐＨ依存値に対するｐＨ値を示す検量線を、前記被検体に関する情報に応じて予め用意された複数の検量線から、前記被検体に関する情報に応じて選択し、
   前記選択された検量線と前記ｐＨ依存値とに基づいて、ｐＨ値を決定すること、
   を備えるｐＨ算出方法。

Images