(57)【特許請求の範囲】
【請求項１】
磁気共鳴画像法（ＭＲＩ）システムであって、
少なくとも一つのエコー時間を使用して解剖学的領域を撮像することにより第一のフリップ角に対応する第一の磁気共鳴（ＭＲ）データセットおよび第二のフリップ角に対応する第二のＭＲデータセットを取得するように構成されたＭＲＩスキャナと、
少なくとも一つのプロセッサと、
メモリであって、そこに格納されたコンピュータコード命令を備え、前記少なくとも一つのプロセッサにより実行されると、前記コンピュータコード命令により前記少なくとも一つのプロセッサは、
前記第一のＭＲデータセット、前記第二のＭＲデータセット、および前記解剖学的領域内の送信無線周波数（ＲＦ）磁場に対する定数値を使用する前記解剖学的領域内の見かけ上の縦緩和時間（Ｔ1app）の空間分布を表すＴ1appマップを生成し、
前記解剖学的領域内の第一の組織タイプに関連付けられた縦緩和時間（Ｔ1）の第一の定数値に基づいて前記Ｔ1appマップをスケーリングすることにより第一の送信ＲＦ磁場マップを推定し、
前記解剖学的領域内の第二の組織タイプに関連付けられたＴ1の第二の定数値に基づいて前記Ｔ1appマップをスケーリングすることにより第二の送信ＲＦ磁場マップを推定し、
推定された第一の送信ＲＦ磁場マップおよび推定された第二の送信ＲＦ磁場マップを使用して前記解剖学的領域内の送信ＲＦ磁場の空間分布を表す第三の送信ＲＦ磁場マップを生成する、
メモリと、
を備える、ＭＲＩシステム。
【請求項２】
前記解剖学的領域が人間の脳であり、前記第一の組織タイプが白質であり、前記第二の組織タイプが灰白質である、請求項１に記載のＭＲＩシステム。
【請求項３】
前記少なくとも一つのプロセッサにより実行されると、前記コンピュータコード命令により前記少なくとも一つのプロセッサがさらに、
前記解剖学的領域内の第三の組織タイプに関連付けられたＴ1の第三の定数値に基づいて前記Ｔ1appマップをスケーリングすることにより第四の送信ＲＦ磁場マップを推定し、
前記第三の送信ＲＦ磁場マップを生成することは、推定された第一の送信ＲＦ磁場マップ、前記された第二の送信ＲＦ磁場マップ、および推定された第四の送信ＲＦ磁場マップを使用することを含む、
請求項１に記載のＭＲＩシステム。
【請求項４】
前記解剖学的領域が人間の脳であり、前記第一の組織タイプが白質であり、前記第二の組織タイプが灰白質であり、前記第三の組織タイプが脳脊髄液である、請求項３に記載のＭＲＩシステム。
【請求項５】
前記少なくとも一つのプロセッサにより実行されると、前記コンピュータコード命令により前記少なくとも一つのプロセッサがさらに、
前記Ｔ1appマップを前記第三の送信ＲＦ磁場マップの二乗で除算することにより前記解剖学的領域内の前記Ｔ1の空間分布を表す第一の縦緩和時間（Ｔ1）マップを生成する、
請求項１に記載のＭＲＩシステム。
【請求項６】
前記少なくとも一つのプロセッサにより実行されると、前記コンピュータコード命令により前記少なくとも一つのプロセッサがさらに、
前記第一のＭＲデータセット、前記第二のＭＲデータセット、および前記送信ＲＦ磁場に対する前記定数値を使用して前記解剖学的領域内の見かけ上のスピン密度（ＳＤapp）の空間分布を表すＳＤappマップを生成し、
前記第三の送信ＲＦ磁場マップに基づいて前記ＳＤappマップをスケーリングすることにより第一のスピン密度（ＳＤ）マップを生成し、
少なくとも前記第一の組織タイプに対応する前記解剖学的領域の第一の副領域と前記第二の組織タイプに対応する前記解剖学的領域の第二の副領域とが等強度であるように前記第三の送信ＲＦ磁場マップ、第一のＴ1マップ、および前記第一のＳＤマップを使用して第三のフリップ角に対応する第三のＭＲデータセットを合成し、
合成された第三のＭＲデータセットを使用して前記解剖学的領域内の受信ＲＦ磁場の空間分布を表す受信ＲＦ磁場マップを推定し、
推定された受信ＲＦ磁場マップに基づいて前記第一のＳＤマップをスケーリングすることにより第二のＳＤマップを生成する、
請求項１に記載のＭＲＩシステム。
【請求項７】
前記少なくとも一つのプロセッサにより実行されると、前記コンピュータコード命令により前記少なくとも一つのプロセッサが、
複数の被験体に対する複数の送信ＲＦ磁場マップを推定し、
前記複数の被験体に対する推定された複数の送信ＲＦ磁場マップの平均化を使用して送信ＲＦ磁場テンプレートを生成する、
請求項１に記載のＭＲＩシステム。
【請求項８】
前記少なくとも一つのプロセッサにより実行されると、前記コンピュータコード命令により前記少なくとも一つのプロセッサが、前記Ｔ1appマップおよび前記送信ＲＦ磁場テンプレートを使用して第二のＴ1マップを生成する、請求項７に記載のＭＲＩシステム。
【請求項９】
前記少なくとも一つのプロセッサにより実行されると、前記コンピュータコード命令により前記少なくとも一つのプロセッサが、推定されたＳＤappマップ、前記送信ＲＦ磁場テンプレート、および推定された受信ＲＦ磁場マップを使用してスピン密度（ＳＤ）マップを生成する、請求項７に記載のＭＲＩシステム。
【請求項１０】
前記少なくとも一つのプロセッサにより実行されると、前記コンピュータコード命令により前記少なくとも一つのプロセッサが、
前記受信ＲＦ磁場マップに基づいて前記第一のＭＲデータセットをスケーリングし、
前記受信ＲＦ磁場マップに基づいて前記第二のＭＲデータセットをスケーリングし、
スケーリングされた第二のＭＲデータセットからのスケーリングされた第一のＭＲデータセットの加重減算を表す画像を生成する、
請求項６に記載のＭＲＩシステム。
【請求項１１】
前記少なくとも一つのプロセッサにより実行されると、前記コンピュータコード命令により前記少なくとも一つのプロセッサが、前記第二のＭＲデータセットからの前記第一のＭＲデータセットの加重減算を表す画像を生成する、請求項１に記載のＭＲＩシステム。
【請求項１２】
磁気共鳴画像法（ＭＲＩ）のための方法であって、
少なくとも一つのプロセッサにより、第一のフリップ角に対応する第一の磁気共鳴（ＭＲ）データセットおよび第二のフリップ角に対応する第二のＭＲデータセットを受信することであって、前記第一のＭＲデータセットおよび前記第二のＭＲデータセットは少なくとも一つのエコー時間を使用して解剖学的領域を撮像することにより取得される、受信することと、
少なくとも一つのプロセッサにより、見かけ上の縦緩和時間（Ｔ1app）マップを
生成することであって、Ｔ1appマップは前記第一のＭＲデータセット、前記第二のＭＲデータセット、および前記解剖学的領域内の送信無線周波数（ＲＦ）磁場に対する定数値を使用して前記解剖学的領域内のＴ1appの空間分布を表す、生成することと、
前記少なくとも一つのプロセッサにより、縦緩和時間（Ｔ1）の第一の定数値に基づいて前記Ｔ1appマップをスケーリングすることにより第一の送信ＲＦ磁場マップを推定することであって、前記Ｔ1の第一の定数値は前記解剖学的領域内の第一の組織タイプに関連付けられる、推定することと、
前記少なくとも一つのプロセッサにより、Ｔ1の第二の定数値に基づいて前記Ｔ1appマップをスケーリングすることにより第二の送信ＲＦ磁場マップを推定することであって、前記Ｔ1の第二の定数値は前記解剖学的領域内の第二の組織タイプに関連付けられる、推定することと、
前記少なくとも一つのプロセッサにより、推定された第一の送信ＲＦ磁場マップおよび推定された第二の送信ＲＦ磁場マップを使用して第三の送信ＲＦ磁場マップを
生成することであって、前記第三の送信ＲＦ磁場マップは前記解剖学的領域内の前記送信ＲＦ磁場の空間分布を表す、生成することと、
を含む、方法。
【請求項１３】
前記少なくとも一つのプロセッサにより、Ｔ1の第三の定数値に基づいて前記Ｔ1appマップをスケーリングすることにより第四の送信ＲＦ磁場マップを推定することであって、前記Ｔ1の第三の定数値は前記解剖学的領域内の第三の組織タイプに関連付けられる、推定すること、
をさらに含み、
前記第三の送信ＲＦ磁場マップを生成することは、推定された第一の送信ＲＦ磁場マップ、推定された第二の送信ＲＦ磁場マップ、および推定された第四の送信ＲＦ磁場マップを使用することを含む、
請求項１２に記載の方法。
【請求項１４】
前記少なくとも一つのプロセッサにより、前記Ｔ1appマップを前記第三の送信ＲＦ磁場マップの二乗で除算することにより第一の縦緩和時間（Ｔ1）を生成することであって、第一のＴ1マップは前記解剖学的領域内のＴ1の空間分布を表す、生成すること、
をさらに含む、請求項１２に記載の方法。
【請求項１５】
前記少なくとも一つのプロセッサにより、前記第一のＭＲデータセット、前記第二のＭＲデータセット、および前記送信ＲＦ磁場に対する前記定数値を使用して見かけ上のスピン密度（ＳＤapp）マップを生成することであって、前記ＳＤappマップは前記解剖学的領域内のＳＤappの空間分布を表す、生成することと、
前記少なくとも一つのプロセッサにより、前記第三の送信ＲＦ磁場マップに基づいて前記ＳＤappマップをスケーリングすることにより第一のスピン密度（ＳＤ）マップを生成することと、
前記少なくとも一つのプロセッサにより、少なくとも前記第一の組織タイプに対応する前記解剖学的領域の第一の副領域と前記第二の組織タイプに対応する前記解剖学的領域の第二の副領域とが等強度であるように、第三のＲＦ磁場マップ、前記第一のＴ1マップ、および前記第一のＳＤマップを使用して第三のフリップ角に対応する第三のＭＲデータセットを合成することと、
前記少なくとも一つのプロセッサにより、合成された第三のＭＲデータセットを使用して受信ＲＦ磁場マップを推定することであって、前記受信ＲＦ磁場マップは前記解剖学的領域内の受信ＲＦ磁場の空間分布を表す、推定することと、
前記少なくとも一つのプロセッサにより、推定された受信ＲＦ磁場マップに基づいて前記第一のＳＤマップをスケーリングすることにより第二のＳＤマップを生成することと、
をさらに含む、請求項１４に記載の方法。
【請求項１６】
前記少なくとも一つのプロセッサにより、複数の被験体に対する複数の送信ＲＦ磁場マップを推定することと、
前記少なくとも一つのプロセッサにより、前記複数の被験体に対する推定された複数の送信ＲＦ磁場マップの平均化を使用して送信ＲＦ磁場テンプレートを生成することと、
をさらに含む、請求項１２に記載の方法。
【請求項１７】
前記少なくとも一つのプロセッサにより、推定されたＳＤappマップ、前記送信ＲＦ磁場テンプレート、および推定された受信ＲＦ磁場マップを使用してスピン密度（ＳＤ）マップを生成すること、
をさらに含む、請求項１６に記載の方法。
【請求項１８】
前記少なくとも一つのプロセッサにより、前記Ｔ1appマップおよび前記送信ＲＦ磁場テンプレートを使用して第二のＴ1マップを生成すること、
をさらに含む、請求項１６に記載の方法。
【請求項１９】
前記少なくとも一つのプロセッサにより、前記受信ＲＦ磁場マップに基づいて前記第一のＭＲデータセットをスケーリングすることと、
前記少なくとも一つのプロセッサにより、前記受信ＲＦ磁場マップに基づいて前記第二のＭＲデータセットをスケーリングすることと、
前記少なくとも一つのプロセッサにより、スケーリングされた第二のＭＲデータセットからのスケーリングされた第一のＭＲデータセットの加重減算を表す画像を生成することと、
をさらに含む、請求項１５に記載の方法。
【請求項２０】
非一時的なコンピュータ可読媒体であって、そこに格納されたコンピュータコード命令を備え、少なくとも一つのプロセッサにより実行されると、前記コンピュータコード命令により前記少なくとも一つのプロセッサは、
少なくとも一つのエコー時間を使用して解剖学的領域を撮像することにより取得される、第一のフリップ角に対応する第一の磁気共鳴（ＭＲ）データセットおよび第二のフリップ角に対応する第二のＭＲデータセットを受信し、
前記第一のＭＲデータセット、前記第二のＭＲデータセット、および前記解剖学的領域内の送信無線周波数（ＲＦ）磁場に対する定数値を使用する前記解剖学的領域内の見かけ上の縦緩和時間（Ｔ1app）の空間分布を表すＴ1appマップを生成し、
縦緩和時間（Ｔ1）の第一の定数値に基づいて前記Ｔ1appマップをスケーリングすることにより前記解剖学的領域内の第一の組織タイプに関連付けられた第一の送信ＲＦ磁場マップを推定し、
Ｔ1の第二の定数値に基づいて前記Ｔ1appマップをスケーリングすることにより前記解剖学的領域内の第二の組織タイプに関連付けられた第二の送信ＲＦ磁場マップを推定し、
推定された第一の送信ＲＦ磁場マップおよび推定された第二の送信ＲＦ磁場マップを使用して前記解剖学的領域内の前記送信ＲＦ磁場の空間分布を表す第三の送信ＲＦ磁場マップを生成する、
非一時的なコンピュータ可読媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は概して磁気共鳴画像法（ＭＲＩ）の分野に関する。より具体的には、本開示は可変フリップ角（ＶＦＡ）技術を使用した定量磁気共鳴画像法（ｑＭＲＩ）の方法およびシステムに関する。
【背景技術】
【０００２】
磁気共鳴画像法（ＭＲＩ）はスキャンされた関心対象の構造を再構成するために磁場を使用するイメージングモダリティである。ＭＲＩスキャナは０．１０テスラ（Ｔ）から７Ｔの範囲の磁場のような強力な静磁場を生成するための磁石ならびにＲＦ信号を送信および／または受信するための無線周波数（ＲＦ）トランシーバを含む。生成された静磁場に体が配置されると、体内の水素の陽子が磁場に沿って整列する。横磁場に残存するバルク磁化が存在するようにスピンを傾けるための振動Ｂ1磁場の存在下でＲＦパルスが印加される。ＲＦパルスがオフにされると、縦成分が増加し、横成分が減少して水素の陽子は静磁場に沿った列に戻る。サンプリング時間、またはエコー時間、もしくはグラジェントエコー時間と呼ばれる選択された時点において、データが収集され、受信された信号がスキャンされた体またはその一部の画像を再構成するために使用される。本開示では、収集されたＭＲデータセットに基づいてＭＲ画像を再構成するための様々なＶＦＡ技術が説明される。
【発明の概要】
【０００３】
少なくとも一つの態様によると、磁気共鳴画像法（ＭＲＩ）システムはＭＲＩスキャナ、少なくとも一つのプロセッサ、およびそこに格納されたコンピュータコード命令を備えるメモリを含むことができる。ＭＲＩスキャナは少なくとも一つのエコー時間を使用して解剖学的領域を撮像することにより、第一のフリップ角に対応する第一の磁気共鳴（ＭＲ）データセットおよび第二のフリップ角に対応する第二のＭＲデータセットを取得するように構成することができる。コンピュータコード命令は少なくとも一つのプロセッサにより実行されると少なくとも一つのプロセッサに見かけ上の縦緩和時間（Ｔ1app）を生成させ、このＴ1appマップは第一のＭＲデータセットおよび第二のＭＲデータセットを使用する解剖学的領域内のＴ1appの空間分布を表している。少なくとも一つのプロセッサは縦緩和時間（Ｔ1）の第一の定数値に基づいてＴ1appマップをスケーリングする（ｓｃａｌｅ）ことにより第一の送信ＲＦ磁場マップを推定することができる。Ｔ1の第一の定数値は解剖学的領域内の第一の組織タイプに関連付けることができる。少なくとも一つのプロセッサはＴ1の第二の定数値に基づいてＴ1appマップをスケーリングすることにより第二の送信ＲＦ磁場マップを推定することができる。Ｔ1の第二の定数値は解剖学的領域内の第二の組織タイプに関連付けることができる。プロセッサは推定された第一の送信ＲＦ磁場マップおよび推定された第二の送信ＲＦ磁場マップを使用して第三の送信ＲＦ磁場マップを生成することができる。第三の送信ＲＦ磁場マップは解剖学的領域内の送信ＲＦ磁場の空間分布を表すことができる。
【０００４】
解剖学的領域は人間の脳とすることができ、第一の組織タイプは白質、第二の組織タイプは灰白質とすることができる。少なくとも一つのプロセッサはＴ1の第三の定数値に基づいてＴ1appマップをスケーリングすることにより第四の送信ＲＦ磁場マップを推定することができる。Ｔ1の第三の定数値は解剖学的領域内の第三の組織タイプに関連付けることができる。第三の送信ＲＦ磁場マップを生成することは、推定された第一の送信ＲＦ磁場マップ、推定された第二の送信ＲＦ磁場マップ、および推定された第四の送信ＲＦ磁場マップを使用することを含むことができる。解剖学的領域は人間の脳とすることができ、第一の組織タイプは白質、第二の組織タイプは灰白質、第三の組織タイプは脳脊髄液とすることができる。少なくとも一つのプロセッサはＴ1appマップを第三の送信ＲＦ磁場マップの二乗で除算することにより第一の縦緩和時間（Ｔ1）マップを生成することができる。第一のＴ1マップは解剖学的領域内のＴ1の空間分布を表すことができる。
【０００５】
少なくとも一つのプロセッサは第一のＭＲデータセットおよび第二のＭＲデータセットを使用して見かけ上のスピン密度（ＳＤapp）マップを生成することができる。ＳＤappマップは解剖学的領域内のＳＤappの空間分布を表すことができる。少なくとも一つのプロセッサは第三の送信ＲＦ磁場マップに基づいてＳＤappマップをスケーリングすることにより第一のスピン密度（ＳＤ）マップを生成することができる。少なくとも一つのプロセッサは、少なくとも第一の組織タイプに対応する解剖学的領域の第一の副領域と第二の組織タイプに対応する解剖学的領域の第二の副領域とが等強度であるように、第三の送信ＲＦ磁場マップ、第一のＴ１マップ、および第一のＳＤマップを使用して第三のフリップ角に対応する第三のＭＲデータセットを合成することができる。少なくとも一つのプロセッサは合成された第三のＭＲデータセットを使用して受信ＲＦ磁場マップを推定することができる。受信ＲＦ磁場マップは解剖学的領域内の受信ＲＦ磁場の空間分布を表すことができる。少なくとも一つのプロセッサは推定された受信ＲＦ磁場マップに基づいて第一のＳＤマップをスケーリングすることにより第二のＳＤマップを生成することができる。
【０００６】
少なくとも一つのプロセッサは受信ＲＦ磁場マップに基づいて第一のＭＲデータセットをスケーリングし、受信ＭＲデータセットに基づいて第二のＭＲデータセットをスケーリングし、スケーリングされた第二のＭＲデータセットからのスケーリングされた第一のＭＲデータセットの加重減算を表す画像を生成することができる。少なくとも一つのプロセッサは第二のＭＲデータセットからの第一のＭＲデータセットの加重減算を表す画像を生成することができる。
【０００７】
少なくとも一つのプロセッサは複数の被験体に対する複数の送信ＲＦ磁場マップを推定し、複数の被験体に対する推定された複数の送信ＲＦ磁場マップの平均化を使用して送信ＲＦ磁場テンプレートを生成することができる。少なくとも一つのプロセッサはＴ1appマップおよび送信ＲＦ磁場テンプレートを使用して第二のＴ1マップを生成することができる。少なくとも一つのプロセッサは推定されたＳＤappマップ、送信ＲＦ磁場テンプレート、および推定された受信ＲＦ磁場マップを使用してスピン密度（ＳＤ）マップを生成することができる。
【０００８】
少なくとも一つの態様によると、磁気共鳴画像法（ＭＲＩ）のための方法は第一のフリップ角に対応する第一の磁気共鳴（ＭＲ）データセットおよび第二のフリップ角に対応する第二のＭＲデータセットを受信する少なくとも一つのプロセッサを含むことができる。第一のＭＲデータセットおよび第二のＭＲデータセットは少なくとも一つのエコー時間を使用して解剖学的領域を撮像することにより取得することができる。本方法は見かけ上の縦緩和時間（Ｔ1app）マップを生成する少なくとも一つのプロセッサを含むことができ、Ｔ1appマップは解剖学的領域内の第一のＭＲデータセットおよび第二のＭＲデータセットを使用する解剖学的領域内のＴ1appの空間分布を表している。本方法は縦緩和時間（Ｔ1）の第一の定数値に基づいてＴ1appマップをスケーリングすることにより第一の送信ＲＦ磁場マップを推定する少なくとも一つのプロセッサを含むことができる。Ｔ1の第一の定数値は解剖学的領域内の第一の組織タイプに関連付けることができる。本方法はＴ1の第二の定数値に基づいてＴ1appマップをスケーリングすることにより第二の送信ＲＦ磁場マップを推定する少なくとも一つのプロセッサを含むことができる。Ｔ1の第二の定数値は解剖学的領域内の第二の組織タイプに関連付けることができる。本方法は推定された第一の送信ＲＦ磁場マップおよび推定された第二の送信ＲＦ磁場マップを使用して第三の送信ＲＦ磁場マップを生成する少なくとも一つのプロセッサを含むことができる。第三の送信ＲＦ磁場マップは解剖学的領域内の送信ＲＦ磁場の空間分布を表すことができる。
【０００９】
本方法はＴ1の第三の定数値に基づいてＴ1appマップをスケーリングすることにより第四の送信ＲＦ磁場マップを推定する少なくとも一つのプロセッサをさらに含むことができる。Ｔ1の第三の定数値は解剖学的領域内の第三の組織タイプに関連付けることができる。第三の送信ＲＦ磁場マップを生成することは、推定された第一の送信ＲＦ磁場マップ、推定された第二の送信ＲＦ磁場マップ、および推定された第四の送信ＲＦ磁場マップを使用することを含むことができる。解剖学的領域は人間の脳とすることができ、第一の組織タイプは白質、第二の組織タイプは灰白質、第三の組織タイプは脳脊髄液とすることができる。本方法はＴ1appマップを第三の送信ＲＦ磁場マップの二乗で除算することにより第一の縦緩和時間（Ｔ1）マップを生成する少なくとも一つのプロセッサをさらに含むことができる。第一のＴ1マップは解剖学的領域内のＴ1の空間分布を表すことができる。
【００１０】
本方法は第一のＭＲデータセットおよび第二のＭＲデータセットを使用して見かけ上のスピン密度（ＳＤapp）マップを生成する少なくとも一つのプロセッサをさらに含むことができる。ＳＤappマップは解剖学的領域内のＳＤappの空間分布を表すことができる。少なくとも一つのプロセッサは第三の送信ＲＦ磁場マップに基づいてＳＤappマップをスケーリングすることにより第一のスピン密度（ＳＤ）マップを生成することができる。少なくとも一つのプロセッサは、少なくとも第一の組織タイプに対応する解剖学的領域の第一の副領域と第二の組織タイプに対応する解剖学的領域の第二の副領域とが等強度であるように、第三の送信ＲＦ磁場マップ、第一のＴ1マップ、および第一のＳＤマップを使用して第三のフリップ角に対応する第三のＭＲデータセットを合成することができる。少なくとも一つのプロセッサは合成された第三のＭＲデータセットを使用して受信ＲＦ磁場マップを推定することができる。受信ＲＦ磁場マップは解剖学的領域内の受信ＲＦ磁場の空間分布を表すことができる。少なくとも一つのプロセッサは推定された受信ＲＦ磁場マップに基づいて第一のＳＤマップをスケーリングすることにより第二のＳＤマップを生成することができる。
【００１１】
本方法は受信ＲＦ磁場マップに基づいて第一のＭＲデータセットをスケーリングし、受信ＭＲデータセットに基づいて第二のＭＲデータセットをスケーリングし、スケーリングされた第二のＭＲデータセットからのスケーリングされた第一のＭＲデータセットの加重減算を表す画像を生成する少なくとも一つのプロセッサをさらに含むことができる。本方法は第二のＭＲデータセットからの第一のＭＲデータセットの加重減算を表す画像を生成する少なくとも一つのプロセッサをさらに含むことができる。
【００１２】
本方法は複数の被験体に対する複数の送信ＲＦ磁場マップを推定し、複数の被験体に対する推定された複数の送信ＲＦ磁場マップの平均化を使用して送信ＲＦ磁場テンプレートを生成する少なくとも一つのプロセッサをさらに含むことができる。本方法はＴ1appマップおよび送信ＲＦ磁場テンプレートを使用して第二のＴ1マップを生成する少なくとも一つのプロセッサをさらに含むことができる。本方法は推定されたＳＤappマップ、送信ＲＦ磁場テンプレート、および推定された受信ＲＦ磁場マップを使用してスピン密度（ＳＤ）マップを生成する少なくとも一つのプロセッサをさらに含むことができる。
【００１３】
少なくとも一つの態様によると、コンピュータ可読媒体はそれに格納されたコンピュータコード命令を含む。コンピュータコード命令は少なくとも一つのプロセッサにより実行されると磁気共鳴画像法（ＭＲＩ）のための方法を実行させる。本方法は第一のフリップ角に対応する第一の磁気共鳴（ＭＲ）データセットおよび第二のフリップ角に対応する第二のＭＲデータセットを受信する少なくとも一つのプロセッサを含むことができる。第一のＭＲデータセットおよび第二のＭＲデータセットは少なくとも一つのエコー時間を使用して解剖学的領域を撮像することにより取得することができる。本方法は見かけ上の縦緩和時間（Ｔ1app）マップを生成する少なくとも一つのプロセッサを含むことができ、Ｔ1appマップは第一のＭＲデータセット、第二のＭＲデータセット、および解剖学的領域内の送信無線周波数（ＲＦ）磁場に対する定数値を使用する解剖学的領域内のＴ1appの空間分布を表す。本方法は縦緩和時間（Ｔ1）の第一の定数値に基づいてＴ1appマップをスケーリングすることにより第一の送信ＲＦ磁場マップを推定する少なくとも一つのプロセッサを含むことができる。Ｔ1の第一の定数値は解剖学的領域内の第一の組織タイプに関連付けることができる。本方法はＴ1の第二の定数値に基づいてＴ1appマップをスケーリングすることにより第二の送信ＲＦ磁場マップを推定する少なくとも一つのプロセッサを含むことができる。Ｔ1の第二の定数値は解剖学的領域内の第二の組織タイプに関連付けることができる。本方法は推定された第一の送信ＲＦ磁場マップおよび推定された第二の送信ＲＦ磁場マップを使用して第三の送信ＲＦ磁場マップを生成する少なくとも一つのプロセッサを含むことができる。第三の送信ＲＦ磁場マップは解剖学的領域内の送信ＲＦ磁場の空間分布を表すことができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の発明概念による、磁気共鳴画像法（ＭＲＩ）システムを例示するブロック図である。
【図２】本開示の発明概念による、磁気共鳴画像法（ＭＲＩ）のための方法２００を例示するフローチャートである。
【図３Ａ−３Ｂ】二つの異なるフリップ角に対応する脳の二つのＭＲ画像を示す。
【図４Ａ−４Ｂ】図４Ａは、見かけ上の縦緩和時間（Ｔ1app）マップの例を描写する画像を示し、図４Ｂは、図４Ａに示されるＴ1appマップを使用して生成される送信ＲＦ磁場（Ｂ1t）マップの第一の推定の画像を示す。
【図５Ａ−５Ｂ】白質および灰白質に対するマスクを表す画像を示す。
【図６Ａ】解剖学的領域内の線におけるＢ1tの第一の推定およびＢ1tの第二の推定のプロットを示す。
【図６Ｂ】ローカル二次フィッティングを用いた場合と用いない場合の統合されたＢ1tマップに関連付けられたプロットを示す。
【図７Ａ−７Ｂ】統合されたＢ1tマップの例および統合されたＢ1tマップにローカル二次フィッティングを適用した後に得られるＢ1tマップの例を例示する画像を示す。
【図８Ａ−８Ｃ】それぞれＴ1マップの推定、見かけ上のスピン密度（ＳＤapp）マップの推定、およびスピン密度（ＳＤ）マップの推定の例を例示する画像を示す。
【図９】ＦＡの範囲にわたる異なる組織タイプに対応する信号の例を示す。
【図１０Ａ−１０Ｃ】それぞれ合成されたＭＲデータセットの例、受信ＲＦ磁場の例、およびスケーリングされたスピン密度マップの例を例示する画像を示す。
【図１１】Ｂ1tおよびＢ1rマップに対する脳テンプレートの例を示す。
【図１２Ａ−１２Ｂ】スケーリングされたＭＲデータセットおよび減算ＭＲ画像を例示する画像を示す。
【図１２Ｃ−１２Ｄ】スケーリングされたＭＲデータセットおよび減算ＭＲ画像を例示する画像を示す。
【発明を実施するための形態】
【００１５】
定量磁気共鳴画像法（ｑＭＲＩ）は異なるタイプの組織間を区別するため、または異常な組織を識別するために信号強度または相対信号強度を測定および／または使用することを含む。組織特性を定量化することにより、正常および異常な組織タイプの区分および分類の両方をするための方法を開発することがより容易になる。定量スピン密度（ＳＤ）、縦緩和時間（Ｔ1）、横磁化減衰レート（Ｔ2*）、および定量磁化率マッピング（ＱＳＭ）の使用は、今日、重要な臨床での用途を有する。例えば、定量Ｔ1イメージングは心臓血管イメージングにおけるアテローム性動脈硬化、腫瘍、脳卒中、および多発性硬化症の研究において潜在的な実用性を有する。また、定量ＳＤマッピングは浮腫および治療後の組織水分量の変化を研究するために重要な役割を果たし、また使用されている。一般に、このタイプの組織定量化によってより厳密な方法での治療に対する組織の応答を追跡することが可能になる。
【００１６】
定量Ｔ1イメージングおよび定量ＳＤマッピングは可変フリップ角（ＶＦＡ）技術を使用して達成することができる。ＶＦＡ法を使用することは複数の異なるフリップ角に対するＭＲデータの取得を含む。具体的には、ＭＲＩスキャナは異なる名目上のフリップ角を有する複数の損なわれたグラジェントエコーデータセットを取得し、Ｔ1および／またはＳＤを適切に定量化するために取得されたデータセットを使用することができる。ＶＦＡ技術がＭＲデータ収集の容易さにより特徴付けられる一方で、ＶＦＡベースの技術は位置の関数として無線周波数（ＲＦ）励起（または送信）磁場（Ｂ1t）またはフリップ角（ＦＡ）の正確な情報を要求する。実際には、ＶＦＡデータセットに基づいてＴ1マップを正確に計算するための能力は送信ＲＦ磁場（Ｂ1t）マップの空間分布の正確な情報に依存する。また、ＳＤマップの再構成には受信ＲＦ磁場（Ｂ1r）マップの使用を含む。いくつかの実施形態では、ＭＲＩスキャナまたは各プロセッサは例えば他の二つのフリップ角（ＦＡ）と比較される相対的に大きなＦＡに関連付けられたデータセットを含む少なくとも三つの異なるフリップ角に対して取得されたデータセットを使用することにより送信ＲＦ磁場Ｂ1tを算出することができる。三未満のＦＡまたは相対的に小さなＦＡに対応するデータセットを使用することより不適切に送信ＲＦ磁場（Ｂ1t）を再構成する問題が生じる可能性があり、またその問題の不適切さを改善するために追加の制約が要求される可能性がある。不適切さに対処するためのこのような制約の例は胸部イメージングに対するＢ1t磁場分布を決定するために脂肪等の単一の組織を使用すること、または白質もしくは灰白質に対するＴ1とＳＤとの間の関係に制約を使用することを含むことができる。
【００１７】
本開示では、改善された画像品質および定量データを用いる戦略的に取得されたグラジェントエコー（ＳＴＡＧＥ）イメージングのための方法およびシステムが開示される。ＳＴＡＧＥイメージングにより二つのＦＡに対応する取得されたデータセットを使用した包括的で迅速なイメージングが可能になる。また、本明細書に記載のＳＴＡＧＥイメージングのためのシステムおよび方法によりＴ1およびＳＤマップの信頼性が高く正確な再構成が可能になる。
【００１８】
図１は本開示の発明概念による、磁気共鳴画像法（ＭＲＩ）システム１００を例示するブロック図である。概要では、ＭＲＩシステム１００はＭＲＩスキャナ１０２、プロセッサ１０４、メモリ１０６、およびディスプレイデバイス１０８を含むことができる。プロセッサ１０４はＭＲＩスキャナ１０２、メモリ１０６およびディスプレイデバイス１０８と通信可能に結合することができる。
【００１９】
ＭＲＩスキャナ１０２は０．１テスラ（Ｔ）から７Ｔの範囲の磁場のような強力な静磁場を生成するための磁石（図１に示されていない）、ならびにＲＦ信号を送信および／または受信するための複数の無線周波数（ＲＦ）コイル（図１に示されていない）を含むことができる。ＲＦコイルは送信ＲＦコイルおよび受信ＲＦコイルを含むことができる。ＲＦ送信コイルはＭＲＩパルスシーケンスに従って患者の解剖学的領域等の被験体を励起するためのＲＦパルスを放出することができる。受信ＲＦコイルはＲＦパルスの放出を追跡する被験体により生成されるＭＲＩ信号を記録することができる。ＲＦコイルはＲＦ信号を交互に送信および受信することのできるＲＦトランシーバを含み得る。ＲＦコイルは損なわれたＲＦグラジェントデータ取得に従ってＭＲＩデータを取得することができる。記録されたＭＲＩ信号は二つの異なるＦＡに関連付けることができる。
【００２０】
イメージングシステム１００は一または二以上のプロセッサ１０４を含むことができる。一または二以上のプロセッサ１０４はＭＲＩスキャナ１０２内部に集積されたプロセッサ、ＭＲＩスキャナ１０２に通信可能に結合されたコンピューティングデバイスのプロセッサ、またはその組合せを含むことができる。メモリ１０６はＭＲＩスキャナ１０２のメモリコンポーネント、ＭＲＩスキャナ１０２に通信可能に結合されたコンピューティングデバイスのメモリコンポーネント、またはその組合せを含むことができる。メモリ１０６はコンピュータ実行可能命令を含むことができ、この命令は一または二以上のプロセッサ１０４により実行されると一または二以上のプロセッサ１０４に本明細書に記載のＳＴＡＧＥイメージングのための方法を実行させることができる。メモリ１０６はＭＲＩスキャナ１０２により取得されるＭＲＩデータを格納することができ、プロセッサ１０４はメモリ１０６からのこのようなデータにアクセスすることができる。メモリ１０６はスキャナ１０２により取得されるＭＲＩデータに基づいてプロセッサ１０４により生成される画像を受信および格納することができる。
【００２１】
ディスプレイデバイス１０８はカソードレイチューブ（ＣＲＴ）ディスプレイ、発光ダイオード（ＬＥＤ）ディスプレイ、有機発光ダイオード（ＯＬＥＤ）ディスプレイ、プラズマディスプレイパネル（ＰＤＰ）、液晶ディスプレイ（ＬＣＤ）、または当該技術分野における当業者に知られている他のディスプレイを含むことができる。ディスプレイデバイス１０８はスタンドアロンデバイスまたはＭＲＩスキャナ１０２に通信可能に結合されたコンピューティングデバイス（例えば、デスクトップ、ラップトップ、またはタブレット）のディスプレイであり得る。ディスプレイデバイス１０８はタッチスクリーンを含むことができる。ディスプレイデバイス１０８はプロセッサ１０４またはメモリ１０６から画像データを受信し、その受信した画像データを表示することができる。例えば、ＭＲＩスキャナ１０２により取得されるデータに基づいてＭＲＩ画像を再構成する際、プロセッサ１０４はディスプレイデバイス１０８で表示するための再構成された画像を提供することができる。フリップ角θの関数としての損なわれたＲＦグラジェントエコーデータ取得に対する信号強度は以下のように表すことができる。
【数１】
0006872079.tif 000002
Ｔ1は縦緩和時間、ＳＤはスピン密度、ＴＲは繰返し時間、ＴＥはエコー時間、Ｂ1tは送信ＲＦ磁場、Ｂ1rは受信ＲＦ磁場（ＲＦコイル感受性またはバイアス磁場とも呼ばれる）であり、Ｅ１＝ｅｘｐ（−ＴＲ／Ｔ1）である。ここで、Ｂ1tの項は補正フリップ角が得られる時に１に対して正規化されると仮定されている。（１−Ｅ１ｃｏｓ（Ｂ1t・θ））／ｓｉｎ（Ｂ1t・θ）を式（１）の両側にかけると、式（１）は以下のように書き直すことができる。
【数２】
0006872079.tif 000003
ここで、ＳＤeff＝ＳＤ・Ｂ1r・e-TE/T2*は有効スピン密度を表す。
【００２２】
式（２）に従うと、Ｓ（θ）/ｓｉｎ（Ｂ1tθ）はＳ（θ）/ｔａｎ（Ｂ1tθ）の線形関数とみなすことができる。具体的には、異なる名目上の励起フリップ角に対して収集されるデータに対し、変換されたデータをＳＤeff・（１−Ｅ1）において傾きＥ1および縦軸切片（またはｘ軸交差）を有する線にフィッティングすることができる。このように、（Ｓ（θ）/ｔａｎ（Ｂ1tθ），Ｓ（θ）/ｓｉｎ（Ｂ1tθ））座標系における少なくとも二つのＦＡに対応する少なくとも二つのデータ点により定義される線の傾きを算出することによりＥ1の値を決定することができる。また、決定されたＥ1値およびｘ軸と少なくとも二つのデータ点により形成される線との間の交点を使用して、ＳＤeffの値を決定することができる。しかしながら、Ｂ1t磁場の不均一性に起因して、測定されたフリップ角はスキャンを動作させるために選択された実際のフリップ角の値とは異なる可能性があり、これにより特に高磁場でのＴ1の推定において著しい誤差が導かれる可能性がある。したがって、Ｔ1の正確な推定では送信ＲＦ磁場Ｂ1tの正確な情報が要求される。
【００２３】
ＦＡが小さく、またＴＲ＜＜Ｔ1である場合、式（１）は以下のように近似することができる。
【数３】
0006872079.tif 000004
ここで角度θEはｃｏｓθE＝ｅｘｐ（−ＴＲ/Ｔ1）と定義される。式（３）を使用すると、見かけ上のスピン密度（ＳＤapp）および見かけ上の縦緩和時間（Ｔ1app）を以下のように導出することができる。
【数４】
0006872079.tif 000005
【数５】
0006872079.tif 000006
式（３）−（５）を考慮すると、見かけ上のスピン密度ＳＤappはＢ1tに線形に比例し、見かけ上の縦緩和時間Ｔ1appはＢ1t2に線形に比例する。したがってＢ1tにおける１０％の誤差（または実際のＦＡがθ’＝Ｂ1tθである場合のＦＡにおける１０％）がＴ1推定における２０％の誤差およびＳＤ推定における１０％の誤差をもたらす可能性がある。そのため、Ｔ1および／またはＳＤマップの正確な再構成ではＢ1tマップの正確な推定が要求される。
【００２４】
図２は磁気共鳴画像法（ＭＲＩ）のための方法２００を例示するフローチャートである。概要を述べると、方法２００は第一のフリップ角に対応する第一の磁気共鳴（ＭＲ）データセットおよび第二のフリップ角に対応する第二のＭＲデータセットを得る（または受信する）こと（ステップ２０２）、ならびに第一のＭＲデータ、第二のＭＲデータ、および解剖学的領域内の送信無線周波数（ＲＦ）磁場に対する定数値を使用して見かけ上の縦緩和時間（Ｔ1app）マップを生成すること（ステップ２０４）を含むことができる。方法２００は縦緩和時間（Ｔ1）の第一の定数値に基づいてＴ1appマップをスケーリングすることにより第一の送信ＲＦ磁場マップを推定すること（ステップ２０６）、Ｔ1の第二の定数値に基づいてＴ1appマップをスケーリングすることにより第二の送信ＲＦ磁場マップを推定すること（ステップ２０８）、ならびに推定された第一の送信ＲＦ磁場マップおよび推定された第二の送信ＲＦ磁場マップを使用して第三の送信ＲＦ磁場マップを生成すること（ステップ２１０）も含むことができる。
【００２５】
方法２００は第一のフリップ角に対応する第一の磁気共鳴（ＭＲ）データセットおよび第二のフリップ角に対応する第二のＭＲデータセットを得ることを含むことができる（ステップ２０２）。スキャナ１０２は少なくとも一つのエコー時間ＴＥを使用して解剖学的領域を撮像することにより第一のＭＲデータセットおよび第二のＭＲデータセットを取得することができる。プロセッサ１０４はユーザがＭＲデータ取得に対する設定を選択することを可能にするために、例えば、ディスプレイデバイス１０８にユーザインタフェース（ＵＩ）を設けることができる。プロセッサ１０２はＭＲスキャナ１０２に選択された設定に従って解剖学的領域を撮像させることができる。選択された設定は二つの別々のＦＡを使用するＭＲデータの取得を表示することができる。選択された設定は（例えば一または二以上の損なわれたグラジェントエコーＭＲシーケンスを使用する）損なわれたグラジェントエコーＭＲデータ取得の表示であり得る。撮像される解剖学的領域は人間の脳もしくは他の臓器または患者の体の一部を含むことができる。
【００２６】
プロセッサ１０４はＭＲスキャナ１０２（または各ＲＦコイル）にＲＦを用いて解剖学的領域を励起すること、および選択された設定に従って解剖学的領域により生成されるＭＲ信号を記録させることができる。具体的には、ＭＲスキャナ１０２の受信ＲＦコイルは第一のＦＡに関連付けられたＭＲ信号の第一の組および第二のＦＡに関連付けられたＭＲ信号の第二の組を記録することができる。記録されたＭＲ信号と対応するＦＡとの間の関係は式（１）および／または式（２）を満たす。第一および第二のＭＲデータセットを得ることは、ＭＲスキャナ１０２から記録された信号を受信し、各ＭＲデータセット（またはＭＲ信号の組）に対する各ＭＲ画像を生成するプロセッサ１０４を含むことができる。
【００２７】
図３Ａおよび図３Ｂを参照すると、二つの異なるＦＡに対応する脳の二つのＭＲ画像が示されている。例えば、図３Ａは６°に等しいＦＡを使用して取得されたＭＲデータの画像３０２を示しており、図３Ｂは２４°に等しいＦＡを使用して取得されたＭＲデータの画像３０４を示している。図３Ａおよび図３Ｂに例示されるＭＲデータの取得において使用される繰返し時間ＴＲは２５ミリ秒（ｍｓ）に等しい。図３Ａおよび図３Ｂの画像に関連付けられたＦＡの値は例示の目的のために選択されたものであり、限定するものと解釈されるべきではない。例えば、他の角度の値（６°および／または２４°以外）は（二つのＦＡの中で）小さい方のＦＡがエルンスト角未満であり大きい方のＦＡがエルンスト角よりも大きい場合にＭＲデータ取得のために使用することができる。プロセッサ１０４は対応するＦＡに関連付けられた記録された信号の逆フーリエ変換を行うことにより図３Ａおよび図３Ｂに示される画像をそれぞれ再構成することができる。
【００２８】
再び図２を参照すると、方法２００は第一のＭＲデータセット、第二のＭＲデータセット、および解剖学的領域内の送信ＲＦ磁場Ｂ1tに対する定数値を使用して見かけ上の縦緩和時間（Ｔ1app）マップを生成することを含むことができる（ステップ２０４）。Ｔ1appマップは解剖学的領域内のＴ1appの空間分布を表す。一般に、また式（５）を考慮すると、Ｔ1appはＢ1tの空間不均一性に起因して単一の組織に関連付けられた領域の内部であっても空間的に変化する可能性がある。例えば、Ｔ1appの値は一つの組織を表す領域の内部で比較的大きな分散を有する可能性がある。全ての人に対して利用可能な標準Ｂ1tがあれば、次にそれを人によって異なるローカルＴ1値を発見するために使用することができる。しかしながら、Ｂ1tのこのような標準値（または標準マップ）が無いため、Ｂ1tマップを正確に推定するための方法または技術が要求される。Ｔ1appマップを生成すること（ステップ２０４）は、Ｂ1tマップを再構成するために使用されるＴ1appマップの第一の推定を決定することとみなすことができる。
【００２９】
プロセッサ１０４は式（２）を使用して、また解剖学的領域内のＢ1tに対する正規化された定数値を例えば１に等しいと仮定して、画素単位でＴ1appを計算することができる。具体的には、プロセッサ１０４は異なるＦＡに関連付けられたデータ点に基づいて定義される線の式（２）における傾きとして各画素に対するＥ１を最初に決定することができる。具体的には、ＦＡθ1およびθ2にそれぞれ対応するデータ点Ｓ（θ1）およびＳ（θ2）を考慮すると、Ｅ１は以下のように算出することができる。
【数６】
0006872079.tif 000007
Ｅ１＝ｅｘｐ（−ＴＲ／Ｔ1）であるため、プロセッサ１０４は算出されたＥ１および繰返し時間ＴＲを使用してＴ1の推定を決定することができる。プロセッサ１０４は次にＴ1の推定された値およびＢ1tの定数値を使用する式（５）に従って各画素に対するＴ1appを算出することができ、したがってＴ1appマップが生成される。図４Ａを参照すると、Ｔ1appマップの画像４０２が描写されている。図４Ａに示されているＴ1appマップは図３Ａおよび図３Ｂに示されるＭＲデータおよび１に等しいＢ1tの正規化された定数値を使用して、上記で説明したように画素毎に算出される。ＭＲデータは２５ｍｓに等しいＴＲの値を使用して取得される。
【００３０】
方法２００はＴ1の第一の定数値に基づいてＴ1appマップをスケーリングすることにより第一の送信ＲＦ磁場マップを推定し（ステップ２０６）、Ｔ1の第二の定数値に基づいてＴ1appマップをスケーリングすることにより第二の送信ＲＦ磁場マップを推定する（ステップ２０８）プロセッサ１０４を含むことができる。Ｔ1の第一の定数値は解剖学的領域内の第一の組織タイプに対応することができ、Ｔ1の第二の定数値は解剖学的領域内の第二の組織タイプに対応することができる。例えば、第一の組織タイプは白質、第二の組織タイプは灰白質とすることができる。一般に、人間の脳のＴ1測定値を確認するための三次元（３Ｄ）の信頼性が高く正確な標準は存在しない。また、Ｔ1は温度、化学交換およびかん流のような多くの要因から影響を受ける可能性がある。様々な研究者によって記録された白質および灰白質の両方に対するＴ1値の値は実質的なばらつきを示す。しかしながら、このようなばらつきにもかかわらず、白質に対するＴ1値で除算される灰白質に対するＴ1値の比率は異なる研究者による様々な測定値において極めて好適であるように見受けられる。具体的には、比率はおよそ１．７である。したがって、Ｔ1の第一の定数値は１０００ｍｓに設定することができ、一方Ｔ1の第二の定数値は１７００ｍｓに設定することができる。いくつかの実施形態では、白質および灰白質に対してＴ1の他の値を選択することができる。例えば、白質に対するＴ1値は８００から１２００ｍｓの範囲の任意の値に設定することができ、灰白質に対するＴ1値は白質に対するＴ1値の１．７倍に等しい値に設定することができる。また、他のタイプの組織（例えば、灰白質および／または白質以外）に対して、異なるＴ1値を使用することができる。
【００３１】
プロセッサ１０４はＴ1の第一の定数値に基づいて各画素においてＴ1app値をスケーリングする（または除算する）ことによりＢ1tマップの第一の推定を生成することができる。図４Ｂを参照すると、図４Ａに示されるＴ1appマップを使用して生成されるＢ1tマップの第一の推定の画像４０４。図４Ｂに示されるＢ1tマップの第一の推定は１０００ｍｓに等しく、また白質に対応するＴ1の定数値を使用して算出される。プロセッサ１０４はＴ1の第二の定数値に基づいて各画素においてＴ1app値をスケーリングする（または除算する）ことによりＢ1tの第二の推定を生成することができる。プロセッサ１０４はＢ1tマップの第一の推定および第二の推定のうち少なくとも一つを使用して解剖学的領域内の第一の組織タイプおよび第一の組織タイプに対するマスクも生成することができる。例えば、プロセッサ１０４は白質のマスクおよび灰白質のマスクを生成するために白質に対応するＴ1の定数値を使用して決定されるＢ1tマップの第一の推定を使用することができる。Ｂ1tマップの推定は低空間周波数成分により特徴付けられ、解剖学的領域を異なる組織タイプに対応する様々な領域に区分するために適切である可能性がある。プロセッサ１０４は白質領域と灰白質領域との間を区別するためにＢ1tマップの第一の推定に対してハイパスフィルタを適用することができる。プロセッサ１０４は白質領域を決定するためにＢ1tマップの第一の推定に対してローパスフィルタを適用することができ、灰白質領域を決定するためにＢ1tマップの第一の推定に対してハイパスフィルタを適用することができる。プロセッサ１０４は異なる組織タイプに関連付けられた領域間の干渉または重複を回避するために各マスクに対してエロ―ジョンアルゴリズム（ｅｒｏｓｉｏｎ ａｌｇｏｒｉｔｈｍ）を適用することもできる。プロセッサ１０４により生成されるマスクの各々に対して、対応する組織タイプに関連付けられた画素は１に等しい値（またはゼロでない値）を有する可能性があり、一方他の画素はゼロに等しい値を有する可能性がある。
【００３２】
図５Ａおよび図５Ｂを参照すると、白質および灰白質に対するマスクを表す画像が示されている。画像５０２は白質のマスクを表しており、画像５０４は灰白質のマスクを表している。両方のマスクは図４Ｂに示されるＢ1tマップの第一の推定を使用して生成される。
【００３３】
ステップ２１０において、プロセッサ１０４は推定された第一のＢ1tマップおよび推定された第二のＢ1tマップを使用して第三のＢ1tマップを生成することができる。図６Ａを参照すると、解剖学的領域内の線におけるＢ1tの第一の推定およびＢ1tの第二の推定のプロットが示されている。プロット６０２は解剖学的領域内の線におけるＢ1tの第一の推定（または白質に対応するＴ1の定数値を使用するＢ1tの推定）を表している。プロット６０４は解剖学的領域内の線におけるＢ1tの第二の推定（または灰白質に対応するＴ1の定数値を使用するＢ1tの推定）を表している。二つのプロット６０２と６０４とを比較すると、Ｂ1tの第一の推定およびＢ1tの第二の推定が実質的に互いのシフトしたバージョンであることがわかる。第三のＢ1tマップを生成することはＢ1tの第一および第二の推定されたマップの両方が実質的に重複するまでＢ1tの第一および第二の推定されたマップのうち少なくとも一つをシフトさせるプロセッサ１０４を含むことができる。例えば、プロセッサ１０４はＢ1tのシフトされた第一の推定されたマップとＢ1tの第二の推定されたマップとの間の平均二乗誤差が最小化されるようにシフト値に基づいてＢ1tの第一の推定されたマップをシフトさせることができる。
【００３４】
プロセッサ１０４はＢ1tのシフトされた第一の推定されたマップ、Ｂ1tの第二の推定されたマップ、ならびに第一および第二の組織タイプに対応するマスクを使用して第三のＢ1tマップを生成することができる。具体的には、プロセッサは（ｉ）第一の組織タイプ（例えば、白質）に対応するマスクにＢ1tのシフトされた第一の推定されたマップを乗算（すなわち、画素単位で乗算）すること、（ｉｉ）第二の組織タイプ（例えば、灰白質）に対応するマスクにＢ1tの第二の推定されたマップを乗算（すなわち、画素単位で乗算）すること、および（ｉｉｉ）Ｂ1tの統合されたマップ（または第三のＢ1tマップ）を形成するためにこれらの乗算の結果を統合することができる。プロセッサ１０４は生成されたマスクのいずれにも属さない任意の画素を埋める（またはそれに値を割り当てる）ためにＢ1tの統合されたマップ（または第三のＢ1tマップ）にローカル二次フィッティング（ｌｏｃａｌ ｑｕａｄｒａｔｉｃ ｆｉｔｔｉｎｇ）を適用し得る。
【００３５】
図６Ｂを参照すると、ローカル二次フィッティングを用いた場合と用いない場合の統合されたＢ1tマップに関連付けられたプロットが描写されている。プロット６０６は解剖学的領域内の線における統合されたＢ1tの値を表している。プロット６０８はローカル二次フィッティングを適用した後の解剖学的領域内の線における統合されたＢ1tの同じ値を表している。ローカル二次フィッティングはサイズがｍ×ｍのローカル二次フィッティング行列を使用して達成することができ、ここでｍは整数である。ｎ×ｎ画像の全体に滑らかでよりノイズの少ない結果を創造するために使用される。図６Ｂに例示されるように、ローカル二次フィッティング操作はＢ1tの滑らかでよりノイズの少ないマップを導く。
【００３６】
図７Ａを参照すると、統合されたＢ1tマップの例を例示する画像７０２が示されている。具体的には、画像７０２は白質に対応するマスク５０２とＢ1tのシフトされた第一の推定されたマップとの積と、灰白質に対応するマスク５０４とＢ1tの第二の推定されたマップとの積とを統合して得られる画像を表している。画像７０２における黒色領域はマスク５０２またはマスク５０４のいずれにも属さない画素または領域を表している。
【００３７】
図７Ｂを参照すると、画像７０２にローカル二次フィッティングを適用した後に得られるＢ1tマップの画像７０４が示されている。図７０４に例示されるように、ローカル二次フィッティングにより図７０２に示される統合されたＢ1tにおける間隙を埋めることが可能になる。画像７０４に示されるＢ1tマップにより解剖学的領域内のＢ1tの信頼性の高い推定が提供され、またＴ1マップおよび／またはＳＤマップの正確な再構成が可能になる。
【００３８】
上記において方法２００は二つの組織タイプに関して説明されたが、解剖学的領域は二よりも多くの異なる組織タイプを含むことができる。例えば、解剖学的領域は三つの区別できる組織タイプを含むことができる。このような例では、プロセッサ１０４は解剖学的領域内の第三の組織タイプに関連付けられたＴ1の第三の定数値に基づいて（ステップ２０４において生成された）Ｔ1appマップをスケーリングすることによりＢ1tマップの別の推定をさらに決定することができる。プロセッサ１０４は解剖学的領域内の第三の組織タイプにより占有される領域を表す第三のマスクを生成することもできる。例えば、プロセッサ１０４は三つのマスクの各々を生成するために異なるフィルタ（例えば、ローパスフィルタ、バンドパスフィルタ、および／またはハイパスフィルタ）を使用することができる。プロセッサ１０４は（例えば、第三のマスクに従って）Ｔ1の第三の定数値に基づいて推定されるＢ1tのマップを使用して第三のＢ1tマップを更新することができる。あるいは、（例えば、三つの組織タイプに対応する三つのマスクに基づく）Ｂ1tマップの三つの推定を使用して第三のＢ1tマップを最初に生成することができる。例えば、第三の組織タイプは脳脊髄液とすることができ、プロセッサ１０４は（ｉ）白質に対する一定のＴ1値に基づいて推定される第一のＢ1tマップに第一のマスクを乗算すること、（ｉｉ）灰白質に対する一定のＴ1値を使用して推定される第二のＢ1tマップに第二のマスクを乗算すること、（ｉｉｉ）第三のマスクにより脳脊髄液に対する一定のＴ1値を使用して推定されるＢ1tマップを推定すること、また（ｉｖ）他の組織に対しても同様に推定することができる。脳脊髄液に対する一定のＴ1値は４５００ｍｓに等しい値とすることができる。プロセッサ１０４は次に単一のより正確なＢ1tマップ画像を創造するために（マスクを用いた乗算から）得られた画像を統合することができる。
【００３９】
いくつかの実施形態では、プロセッサ１０４は窓における平均および標準偏差を計算し、平均を上回るかまたは下回る所与の閾値（例えば、ノイズの三つの標準偏差）を超えているこれらの点（または画素）を除去するために、例えばスライドする窓を使用することによりローカルノイズスパイクを除去し得る。プロセッサ１０４は滑らかでよりノイズの少ないＢ1tマップ画像を提供するために例えばサイズがｍ×ｍのローカル二次フィッティング行列を使用して統合されたＢ1tマップに図６Ｂに関して上記で説明されたようなローカル二次フィッティングを適用することができ、ここでｍは整数である。いくつかの実施形態では、整数ｍは２０に等しい値とすることができる。しかしながら、ｍに対して他の値を使用してもよい。人間の脳のケースでは組織タイプは白質、灰白質、および脳脊髄液を含むことができるが、本方法およびシステムは人間の脳またはこれらの組織タイプに限定されるものと解釈される。一般に、本明細書に記載の方法は固有の解剖学的領域または解剖学的領域内の固有の数の組織タイプに制限されるべきではない。
【００４０】
方法２００はＴ1マップの第一の推定を生成するためにステップ２１０において決定されたＢ1tマップを使用するプロセッサ１０４をさらに含むことができる。プロセッサ１０４は式（５）を使用すること、およびステップ２０４において決定されたＴ1appマップをステップ２０１において決定されたＢ1tマップの二乗で除算することができる。具体的には、解剖学的領域の各画素またはボクセルｘに対して、プロセッサ１０４は対応するＴ1値をＴ1（ｘ）＝Ｔ1app（ｘ）／Ｂ1t2（ｘ）として算出することができる。そのため、このＴ1マップの第一の推定は解剖学的領域内のＢ1tの分布におけるばらつきを説明し、また補正されたＴ1マップと呼ばれ得る。
【００４１】
プロセッサ１０４は第一のＭＲデータセットおよび第二のＭＲデータセットを使用して見かけ上のスピン密度（ＳＤapp）マップをさらに生成することができる。ＳＤappマップは解剖学的領域内のＳＤappの空間分布を表す。プロセッサ１０４は式（３）およびステップ２０４において仮定されたＢ1tの定数値を使用してＳＤappマップを決定することができる。例えば、プロセッサ１０４はＦＡθ1またはθ2に対応するデータセットからのデータ点およびステップ２０４において仮定されたＢ1tの定数値を使用して各画素においてＳＤappに対して式（３）を解くことができる。あるいは、式（２）に基づいて、また再びステップ２０４を参照すると、プロセッサ１０４はＢ1tの正規化された定数値（例えば、Ｂ1t＝１．０）を仮定する時に（Ｓ（θ）／ｔａｎ(Ｂ1tθ)，Ｓ（θ）／ｓｉｎ（Ｂ1tθ））座標系におけるｘ軸を用いてデータ点（Ｓ（θ1）／ｔａｎ（Ｂ1tθ1），Ｓ（θ1）／ｓｉｎ（Ｂ1tθ1））および（Ｓ（θ2）／ｔａｎ（Ｂ1tθ2），Ｓ（θ2）／ｓｉｎ（Ｂ1tθ2））により定義される線の交点に基づいて各画素においてＳＤeff・（１−Ｅ１）を決定することができる。ステップ２０４においてＥ1が既に決定されているため、プロセッサ１０４は各画素においてｘ軸との交点に対応する値を１−Ｅ1で除算することによりＳＤeffのマップを生成すること、およびＳＤappマップを生成するためにＳＤeffマップにステップ２０４において仮定されたＢ1tの定数値を乗算することができる。
【００４２】
プロセッサ１０４は決定されたＳＤappマップをステップ２１０において決定されたＢ1tマップで除算することによりＳＤマップの第一の推定を生成することができる。式（４）に従って、解剖学的領域内の各画素またはボクセルに対し、プロセッサ１０４は（ステップ２１０において生成されたＢ1tマップからの）同じ画素におけるＢ1t値に基づいてその画素におけるＳＤapp値をスケーリングする（または除算する）ことができる。ＳＤeff＝ＳＤ・Ｂ1r・ｅ-TE/T2*であるため、ＳＤマップの生成された第一の推定はＳＤeffマップの推定である。このＳＤマップの第一の推定によりＢ1tの分布におけるばらつき（ただしＢ1rの分布におけるばらつきではない）が説明される。Ｂ1rデータが得られると、補正されたＳＤマップはＳＤcor＝ＳＤ・ｅ-TE/T2*を介して得ることができ、また補正されたＳＤマップと呼ばれ得る。最後に、Ｔ２*を計算するために多重エコーを使用して、絶対ＳＤを式ＳＤ＝ＳＤcor・ｅTE/T2*から得ることができる。
【００４３】
図８Ａ−図８Ｃを参照すると、画像８０２、８０４および８０６はそれぞれＴ1マップの推定、ＳＤappマップの推定、およびＳＤマップの推定の例を示している。画像８０２に示されるＴ1マップは画像７０４（図７Ｂ）に示されるＢ1tの二乗に基づいて画像４０２（図４Ａに示されている）におけるＴ1appマップを（画素毎に）スケーリングすることにより生成される。画像８０４におけるＳＤappマップの推定は式（３）、画像３０２または３０４（図３Ａおよび図３Ｂ）に例示されるデータセットのうちの一つを使用して生成される。画像８０６のＳＤマップ（またはＳＤeffマップ）の推定は画像７０４のＳＤappマップを画像７０４（図７Ｂ）のＢ1tマップで除算することにより生成される。
【００４４】
プロセッサ１０４はＢ1rのマップを推定する（または生成する）ために（ステップ２１０において決定された）推定されたＢ1tマップ、補正されたＴ1マップおよび補正されたＳＤマップを使用することができる。具体的には、プロセッサ１０４は第一の組織タイプ（例えば、白質）に関連付けられた第一の領域および第二の組織タイプ（例えば、灰白質）に関連付けられた第二の領域を等強度にすることのできる固有のフリップ角に対応するＭＲデータセットを合成するために推定されたＢ1tマップだけでなく補正されたＴ1およびＳＤマップを用いることができる。プロセッサ１０４はＭＲＩスキャナ１０２を介した追加のＭＲデータ取得を用いるのではなく、式（１）を使用することによりこのようなＭＲデータセットを合成することができる。プロセッサ１０４は三つの異なる組織タイプ（例えば、白質、灰白質、および脳脊髄液）に対応する三つの領域が等強度であるＭＲデータセットを創造するためにエコー時間ＴＥをさらに変え得る。プロセッサ１０４は解剖学的領域内のＢ1rの空間ばらつきの影響を緩和するために推定されたＢ1rマップによりＳＤマップ推定（例えば、図８Ｃの画像８０６に示されるＳＤマップ）をスケーリングすること、およびＳＤマップの改善された推定を達成することができる。
【００４５】
図９を参照すると、ＦＡの範囲にわたる異なる組織タイプに対応する信号の例が示されている。信号９０２は白質（またはＷＭ）に対応し、信号９０４は灰白質（またはＧＭ）に対応し、信号９０６は脳脊髄液（またはＣＳＦ）に対応する。信号９０８は信号９０２と９０４との間の差を表し、信号９１０は信号９０４と９０６との間の差を表し、信号９１２は信号９０２と９０６との間の差を表している。信号９０２、９０４、および９０６は白質内の１０００ｍｓ、灰白質内の１７００ｍｓ、そして脳脊髄液内の４５００ｍｓに等しい一定のＴ1値を使用して生成された。繰返し時間ＴＲは３０ｍｓに等しく、白質ＳＤは０．６８に等しく、灰白質のＳＤは０．８２に等しい。また、４４ｍｓおよび５１ｍｓのＴ2*値がそれぞれ白質および灰白質に対して使用された。信号９０２および９０４はおよそ１０°（またはそれをわずかに上回る）ＦＡにおいて交差する。異なる組織を等強度にすることにより、プロセッサ１０４はローカルスピン密度を決定する必要性を回避することができる。全ての信号が等強度であるため、合成されたデータセットにおける増幅のばらつきはＢ1r磁場におけるばらつきを表している。したがって、プロセッサ１０４は合成されたデータセットに基づいてＢ1rマップを決定する（または生成する）ことができる。プロセッサ１０４はＢ1rマップの推定においてローカル二次フィッティングを用いることもできる。当初の画像を補正するために使用されるＢ1rマップは画像の中心におけるＢ1rの値に対して正規化される等強度画像からの信号と定義することができる。
【００４６】
図１０Ａ−図１０Ｃを参照すると、合成されたＭＲデータセットの例、受信ＲＦ磁場の例、およびスケーリングされたスピン密度マップの例を例示する画像が示されている。画像１００２は白質と灰白質とが等強度であるようにＦＡが選択される二重の等強度ＭＲデータセットを表している。画像１００４は画像１００２の等強度ＭＲデータセットに基づいて生成（または推定）されるＢ1rマップを表している。画像１００２における等強度ＭＲデータセットは間隙（黒色領域または画素）を含む。このような間隙はローカル二次フィッティングを使用することにより画像１００４のＢ1rマップにおいて除去される。いくつかの実施形態では、プロセッサ１０４は推定されたＳＤマップ上の解剖学的領域内のＢ1rの空間ばらつきの影響を緩和するために合成されたデータセットに基づいて推定されたＢ1rマップによりＳＤマップ推定（例えば、図８Ｃの画像８０６に示されるＳＤマップ）をスケーリングし得る。画像１００６は画像１００４のＢ1rマップによる画像８０６におけるＳＤマップのスケーリング（画素毎のスケーリング）を表している。画像１００６と８０６とを比較すると、画像１００６におけるＳＤマップが画像８０６と比較される解剖学的領域内の様々な副領域の改善された表現を提供することが見てとれる。
【００４７】
いくつかの実施形態では、プロセッサ１０４（または他のコンピューティングデバイス）は複数の被験体（例えば、患者またはボランティア）に対する複数のＢ1tマップを推定することができる。プロセッサ１０４は方法２００に関して上記で説明されたような各被験体に対するＢ1tマップを推定することができる。各被験体に対して、プロセッサ１０４は各被験体に対する対応するＢ1tマップを解剖学的領域（例えば、脳）の各テンプレートに変換することができる。プロセッサ１０４は（例えば、Ｉｎｓｉｇｈｔ Ｓｅｇｍｅｎｔａｔｉｏｎ， ａｎｄ Ｒｅｇｉｓｔｒａｔｉｏｎ Ｔｏｏｌｋｉｔ（ＩＴＫ）を使用する）テンプレート脳への第一のデータセットの強度画像（ｍａｇｎｉｔｕｄｅ ｉｍａｇｅ）の第一の変換によりこれを達成することができる。プロセッサ１０４はＢ1t磁場をテンプレート空間にマッピングするために得られた変換を使用することができる。プロセッサ１０４は次に複数の被験体に関連付けられたテンプレートにおけるＢ1t値を平均化して共通のＢ1tテンプレートにすることができる。いくつかの実施形態では、この平均化は加重平均とすることができる。例えば、いくつかの被験体に関連付けられたいくつかのテンプレートは他のテンプレートよりも多く、または少なく重み付けされ得る。プロセッサ１０４は患者空間へ戻す逆変換を実行することにより様々な患者に対するＢ1tマップの推定として共通（または平均）Ｂ1tテンプレートを使用することができる。例えば、共通（または平均）Ｂ1tテンプレートはオフラインで算出され、様々な患者に対する（ステップ２０１において生成されたＢ1tマップに関して上記で説明されたような）Ｔ1マップまたはＳＤマップの再構成において使用することができる。共通（または平均）Ｂ1tテンプレートを使用することにより、プロセッサ１０４は各新しくスキャンされた被験体に対するＢ1tマップを生成することを回避し、代わりに予め算出された共通（または平均）Ｂ1tテンプレートを使用することができる。いくつかの実施形態では、プロセッサは解剖学的領域の別々の画像（例えば、Ｔ1マップの別々の画像）を生成するために予め算出された共通（または平均）Ｂltテンプレートおよび患者に固有のＢ1tマップの推定の両方を使用することができる。
【００４８】
プロセッサ１０４（または他のコンピューティングデバイス）は共通（または平均）Ｂ1rテンプレートを生成することもできる。プロセッサ１０４は例えば上記で説明されたような別々の被験体に対して等強度領域を有する合成されたＭＲデータセットに基づいて、複数の被験体（例えば、患者またはボランティア）に対する複数のＢ1rマップを推定することができる。各被験体に対して、プロセッサ１０４は対応するＢ1rマップを解剖学的領域（例えば、脳）の各テンプレートに変換すること、および共通（または平均）Ｂ1rテンプレートを生成するために様々な被験体に対するテンプレートを平均化（例えば、重み付けされたかまたはされない平均化）することができる。共通（または平均）Ｂ1rテンプレートはオフラインで算出することができ、プロセッサ１０４は様々な患者に対するＢ1rマップの推定として共通（または平均）Ｂ1rテンプレートを使用することができる。予め算出された共通（または平均）Ｂ1rテンプレートが利用できるため、各新しくスキャンされた患者に対する各Ｂ1rマップの構成（または計算）を回避することが可能になる。いくつかの実施形態では、プロセッサは予め算出された共通（または平均）Ｂ1rテンプレートおよび解剖学的領域の別々の画像を生成するための患者に固有のＢ1rマップの推定の両方を使用することができる。
【００４９】
図１１を参照すると、Ｂ1tおよびＢ1rに対する脳テンプレートの例が示されている。画像の上列はＢ1tマップに対する共通（または平均）脳テンプレートの例を示している。この共通（または平均）テンプレートは様々な被験体に対して使用することのできるＢ1tマップの推定を表している。画像の下列はＢ1rマップに対する例の共通（または平均）脳テンプレートを例示している。この共通（または平均）テンプレートは様々な被験体に対して使用することのできるＢ1rマップの推定を表している。
【００５０】
各スキャンされた被験体に対して、プロセッサ１０４はその固有の被験体に対する（例えば、ＩＴＫを使用する）平均Ｂ1rテンプレートおよび／または平均Ｂ1rテンプレートを逆変換することができる。例えば、プロセッサは固有の患者に対する解剖学的領域（例えば、脳）のサイズに一致させるために共通（もしくは平均）Ｂ1tテンプレートおよび／または共通（もしくは平均）Ｂ1rテンプレートを調整することができる。プロセッサ１０４はその固有の患者に対するＴ1マップを推定するために逆変換された共通（または平均）Ｂ1tテンプレートを使用し得る。具体的には、プロセッサ１０４はＴ1マップの推定を生成するために逆変換された共通（または平均）Ｂ1tテンプレートの二乗に基づいて患者に対するＴ1appマップをスケーリング（または画素単位もしくはボクセル単位で除算）することができる。このような推定は患者に固有の推定されたＢ1tマップに基づいて実質的に推定されるか、またはＴ1マップの追加の（もしくは第二の）推定とすることができる。プロセッサ１０４は患者に対するＳＤ（またはＳＤeff）マップの推定の時に共通（または平均）Ｂ1tを使用することもできる。具体的には、プロセッサ１０４はＳＤマップの推定を生成するために逆変換された共通（または平均）Ｂ1tテンプレートに基づいて患者に対するＳＤappマップをスケーリング（または画素単位もしくはボクセル単位で除算）し得る。
【００５１】
プロセッサ１０４は合成されたデータセットに基づいて推定されたＢ1rマップに基づいてステップ２０２において取得された（ＦＡθ1およびθ2に対応する）当初のデータセットをスケーリングすること、およびスケーリングされたデータセットを減算することにより解剖学的領域のＭＲ画像を生成することもできる。生成されたＭＲ画像は各Ｔ1値に基づいて解剖学的領域内の様々な組織タイプを示すことができる。プロセッサ１０４は第一のＦＡに対応する第一のスケーリングされたＭＲデータセットを第二のＦＡに対応する他のスケーリングされたデータセットから減算する場合に加重減算を用いることができる。いくつかの実施形態では、プロセッサ１０４はオフラインで算出または再構成された共通（または平均）Ｂ1rテンプレート（またはその逆変換）に基づいてステップ２０２において取得された（ＦＡθ1およびθ2に対応する）当初のデータセットをスケーリングすることができる。いくつかの実施形態では、プロセッサ１０４は推定されたＳＤマップ上の解剖学的領域内のＢ1rの空間ばらつきの影響を緩和するために共通（または平均）Ｂ1rテンプレートに基づいてＳＤマップ推定（例えば、図８Ｃの画像８０６に示されるＳＤマップ）をスケーリングし得る。
【００５２】
いくつかの実施形態では、プロセッサ１０４は解剖学的領域の別の（または代わりの）ＭＲ画像を生成するためにステップ２０２において取得された（ＦＡθ1およびθ2に対応する）当初のデータセットに対応する画像を（スケーリングすることなく）減算することができる。このようなＭＥ画像は各Ｔ1値に基づいて解剖学的領域内の様々な組織タイプを例示することができる。
【００５３】
図１２Ａ−図１２Ｄを参照すると、スケーリングされたＭＲデータセットの画像および減算ＭＲ画像が例示されている。画像１２０２および１２０４はそれぞれ画像１２０２内のＢ1rにおけるマップによる画像３０２および３０４のスケーリングを表している。Ｂ1rマップによるスケーリングが解剖学的領域内のＢ1rの空間ばらつきの影響を緩和するため、画像１２０２および１２０４は対応する当初の画像３０２および３０４に対する改善を表している。しかしながら、画像１２０２および１２０４は依然として解剖学的領域内のＢ1tの空間ばらつきの影響を受けている。画像１２０６は二つのＦＡに対応する当初の取得されたデータセットに対応する画像３０２と３０４との間の減算（例えば、大きい方のＦＡに対応する画像から小さい方のＦＡに対応する画像の減算）を表している。画像１２０８はスケーリングされた画像１２０２と１２０４との間の減算（例えば、大きい方のＦＡに対応する画像から小さい方のＦＡに対応する画像の減算）を表している。画像１２０６および１２０８の両方が当初の取得されたデータに対応する画像３０２および３０４と比較される解剖学的領域内の様々な組織の改善された表現を提供する一方で、画像１２０８は解剖学的領域内のＢ1rの空間ばらつきの影響が元の画像１２０８に内で除去（または緩和）されるため、画像１２０６に対する改善を表現している。
【００５４】
本明細書に記載の方法およびシステムは二つのフリップ角を使用してスキャンされる解剖学的領域の改善された画像を生成するための様々な技術を提供する。これらの方法およびシステムは人間の脳に限定されると解釈されるべきではなく、他の解剖学的領域に対して使用することができる。また、本開示は（例えば、推定されたＴ1マップ、推定されたＳＤマップ、スケーリングされたＳＤマップ、減算ＭＲ画像、スケーリングされた減算ＭＲ画像、推定されたＢ1tマップ、推定されたＢ1rマップ、共通Ｂ1tテンプレート、または共通Ｂ1rテンプレートの）様々なＭＲ画像を生成するために説明される様々な技術を説明する一方で、本開示はこのような技術または対応するＭＲ画像の様々な組合せを包含するものと解釈されるべきである。さらに、プロセス１０４は複数のエコー時間に対する本明細書に記載の様々なＭＲ画像を生成するための技術を適用することができる。例えば、プロセッサ１０４は解剖学的領域の改善された画像を提供するために様々なエコー時間に関連付けられたデータセットに基づいて生成されるＭＲ画像（例えば、Ｂ1tマップ、Ｂ1rマップ、Ｔ1appマップ、Ｔ1マップ、ＳＤappマップ、ＳＤマップ等）に平均化または加重平均を適用することができる。
【００５５】
当業者においては本開示において説明されるプロセスはプロセッサにより実行可能なコンピュータコード命令を使用して実装されることができることを理解されるべきである。コンピュータコード命令はランダムアクセスメモリ（ＲＡＭ）、読取り専用メモリ（ＲＯＭ）、キャッシュメモリ、ディスクメモリ、任意の他のメモリ、または任意の他のコンピュータ可読媒体のような非一時的または有形のコンピュータ可読媒体に格納することができる。少なくとも一つのプロセッサにより実行される時のコード命令は本開示において説明されるプロセスまたは操作のうち少なくとも一つを実行することを引き起こすことができる。装置は、例えばＭＲＩスキャナ、コンピュータデバイス、またはＭＲＩスキャナに関連付けられた他の電子デバイスとすることができる。
【図１】
0006872079.tif 000008
【図２】
0006872079.tif 000009
【図３Ａ】
0006872079.tif 000010
【図３Ｂ】
0006872079.tif 000011
【図４Ａ】
0006872079.tif 000012
【図４Ｂ】
0006872079.tif 000013
【図５Ａ】
0006872079.tif 000014
【図５Ｂ】
0006872079.tif 000015
【図６Ａ】
0006872079.tif 000016
【図６Ｂ】
0006872079.tif 000017
【図７Ａ】
0006872079.tif 000018
【図７Ｂ】
0006872079.tif 000019
【図８Ａ】
0006872079.tif 000020
【図８Ｂ】
0006872079.tif 000021
【図８Ｃ】
0006872079.tif 000022
【図９】
0006872079.tif 000023
【図１０Ａ】
0006872079.tif 000024
【図１０Ｂ】
0006872079.tif 000025
【図１０Ｃ】
0006872079.tif 000026
【図１１】
0006872079.tif 000027
【図１２Ａ】
0006872079.tif 000028
【図１２Ｂ】
0006872079.tif 000029
【図１２Ｃ】
0006872079.tif 000030
【図１２Ｄ】
0006872079.tif 000031