JP6581584B2 - 水／脂肪分離を用いた位相感応型反転回復ｍｒｉ - Google Patents

Description

本発明は、磁気共鳴（ＭＲ：magnetic resonance）撮像の分野に関する。本発明は、ＭＲデバイスの検査ボリューム内に配置された対象物をＭＲ撮像する方法に関する。また、本発明は、ＭＲデバイス、及びＭＲデバイス上で実行されるコンピュータプログラムに関する。

今日、特に医療診断の分野で、２次元又は３次元画像を形成するために磁場と核スピンとの相互作用を利用する画像形成ＭＲ法が広く使用されている。これは、そのような画像形成ＭＲ法が、軟組織の撮像に関して多くの点で他の撮像方法よりも優れており、電離放射線を必要とせず、且つ通常は侵襲性でないからである。

一般にＭＲ法によれば、検査対象の患者の身体が、強く且つ均一な磁場Ｂ_０内に配置され、磁場Ｂ_０の方向は、同時に、測定が基づく座標系の軸（通常はｚ軸）を定義する。磁場Ｂ_０は、磁場強度に依存して個々の核スピンに関して異なるエネルギーレベルを生じ、これらの核スピンは、規定の周波数（いわゆるラーモア（Larmor）周波数、又はＭＲ周波数）の交番電磁場（ＲＦ場）の印加によって励起（スピン共鳴）され得る。巨視的視点では、個々の核スピンの分布が全体的な磁化をもたらし、この磁化は、適切な周波数の電磁パルス（ＲＦパルス）の印加によって平衡状態からずれるように偏向され得、これにより磁化がｚ軸周りで歳差運動を行う。歳差運動は円錐の表面を表し、この円錐の開き角はフリップ角と呼ばれる。フリップ角の大きさは、印加される電磁パルスの強度及び持続時間に依存する。いわゆる９０°パルスの場合、スピンは、ｚ軸から横平面に偏向される（フリップ角９０°）。

ＲＦパルスの終了後、磁化は緩和して元の平衡状態に戻り、この平衡状態で、ｚ方向での磁化は、第１の時定数Ｔ_１（スピン格子又は縦緩和時間）で再び構築され、ｚ方向に垂直な方向での磁化は、第２の時定数Ｔ_２（スピン−スピン、又は横緩和時間）で緩和する。磁化の変化は、受信ＲＦコイルによって検出され得、受信ＲＦコイルは、磁化の変化がｚ軸に垂直な方向で測定されるようにＭＲデバイスの検査ボリューム内部に配置されて向きを定められる。横磁化の減衰に伴って、例えば９０°パルスの印加後、同位相での秩序状態から、全ての位相角が一様に分散される状態（デフェージング（dephasing））への、（局所磁磁場不均一性によって誘発される）核スピンの遷移が生じる。デフェージングは、再合焦パルス（例えば１８０°パルス）によって補償され得る。これは、受信コイルにおいてエコー信号（スピンエコー）をもたらす。

体内での空間分解能を実現するために、３つの主軸に沿って延びる一定の磁場勾配が、均一な磁場Ｂ_０に重畳されて、スピン共鳴周波数の線形空間依存性をもたらす。このとき、受信コイルでピックアップされる信号は、体内の異なる位置に関連付けられ得る異なる周波数の成分を含む。受信コイルによって取得される信号データは、空間周波数領域に対応し、ｋ空間データと呼ばれる。ｋ空間データは通常、異なる位相エンコードで取得された複数のラインを含む。各ラインは、複数の試料を収集することによってデジタル化される。ｋ空間データの組が、フーリエ変換によってＭＲ画像に変換される。

一般に、ＭＲ撮像における磁化準備技法として反転回復（ＩＲ：inversion recovery）が使用される。ＩＲ撮像において、まず、主磁場Ｂ_０に沿った縦磁化が、１８０°ＲＦパルスを使用して−ｚ方向に回転される。次いで、反転された磁化は、反転ＲＦパルスと励起ＲＦパルスとの間の反転回復時間（ＴＩ）中にＴ_１緩和によって回復される。

ＩＲ撮像の既知の用途は、いわゆる位相感応型反転回復（ＰＳＩＲ：phase sensitive inversion recovery）法であり、これは、心筋梗塞の検出及び評価に特に良く適している。ＰＳＩＲ画像再構成での主な課題は、複素画像中の（信号の極性を決定する）固有信号位相を、ＭＲ画像において一般的である位相誤差（以下ではバックグランドフェージングと呼ぶ）から分離するための位相補正プロセスである。特に、バックグランドフェージングは、オフ共鳴（主磁場の空間的変化）による効果を含む。即ち、位相感応型再構成は、ＰＳＩＲにおいて、望ましい信号の極性を保ちながらバックグランドフェージングを除去するために使用される。

ＩＲを有さない又は異なるＴＩでのＩＲを有する別の画像（参照画像）の取得による位相誤差の較正を含め、ＰＳＩＲ画像再構成のための幾つかの手法が提案されている。しかし、これらの手法は、データ取得効率を低下させる。更に、患者の動きによる実際のスキャンと較正スキャンとの空間的な位置ずれが問題となり得る。

ＰＳＩＲ画像再構成のための代替手法は、適切な位相補正アルゴリズムを使用して、ＩＲ画像自体から位相誤差を決定することである。１つのそのような手法は、いわゆる、決定的再構成を用いた位相感応型反転回復の参照なしの取得（ＲＡＰＩＤ：reference-less acquisition of phase sensitive inversion recovery with decisive reconstruction）法であり、この方法は、参照画像に依拠せずに磁化の極性を確実に復元することが可能である（Jinnan Wang et al., Proc. ISMRM 2013, 2077参照）。

心筋トリグリセリド沈着は、ＰＳＩＲベースの心筋瘢痕検出及び定量化における誤差及び誤解釈の一般的な原因である。位相感応型再構成画像において、心筋瘢痕と心筋トリグリセリド沈着とはいずれも高信号をもたらし、このため、これらはほぼ区別不可能になる。心筋トリグリセリド沈着から発する脂肪信号は、最適な遅延時間（ＴＩ＝２００〜２５０ｍｓ）での自然な高速のＴ_１回復により、（例えば既知のＳＰＩＲ／ＳＰＡＩＲ技法による）スペクトル選択性パルスを使用しても容易には抑制することができない。

以上のことから、改良されたＭＲ撮像技法の必要性が容易に理解される。従って、本発明の目的は、心筋瘢痕と心筋トリグリセリド沈着との区別を可能にするＰＳＩＲベースのＭＲ撮像法を提供することである。

本発明によれば、異なるＭＲスペクトルを有する少なくとも２種の化学種をＭＲ撮像する方法が開示される。本発明の方法は、
ａ）ＭＲデバイスの検査ボリューム内に位置付けされた対象物を、ＲＦパルスとスイッチ磁場勾配との撮像シーケンスにさらすことによって、２つ以上の異なるエコー時間にエコー信号を発生させるステップであって、撮像シーケンスが、反転ＲＦパルスと、それに続く、反転回復時間後の励起ＲＦパルスとを含む反転回復シーケンスである、ステップと、
ｂ）エコー信号を取得するステップと、
ｃ）取得されたエコー信号に対する少なくとも２種の化学種の信号寄与を分離するステップと、
ｄ）化学種の少なくとも１種の信号寄与から、位相感応型ＭＲ画像を再構成するステップと
を含む。

本発明は、従来のＰＳＩＲ型ＭＲ技法の上述した限界に対処する。本発明は、水／脂肪のディクソン（Dixon）型撮像と、位相感応型反転回復技法との一体化を提案する。

ＭＲ撮像の分野の当業者には知られているように、全体の信号に対する水及び脂肪等の異なる化学種の相対的な寄与に関する情報を取得して、それらの幾つかの寄与を抑制すること、又はそれら全ての寄与を個々に若しくはまとめて分析することがしばしば望まれる。これらの寄与は、異なるエコー時間に取得される２つ以上の対応するエコーからの情報が組み合わされる場合に計算され得る。これは、化学シフトエンコードとみなされ得、ここで、追加の次元、即ち化学シフトの次元は、わずかに異なるエコー時間で画像の対を取得することによって定義されてエンコードされる。特に水と脂肪との分離に関して、これらのタイプの実験は、しばしばディクソン型の測定と呼ばれる。ディクソン撮像又はディクソン水／脂肪撮像により、異なるエコー時間に取得された２つ以上の対応するエコーから水及び脂肪の寄与を計算することによって水／脂肪分離が実現され得る。一般に、脂肪中と水中とで水素の歳差周波数の既知の差があるため、そのような分離が可能である。その最も単純な形態では、水画像と脂肪画像とは、「同位相（in phase）」データセットと「逆位相（out of phase）」データセットとの加算又は減算によって生成される。多くの既知のディクソン水／脂肪分離アルゴリズムは、水画像及び脂肪画像だけでなく、Ｂ_０磁場不均一性のマップ、即ち、いわゆるＢ_０マップも提供する。

本発明によれば、反転回復を使用して２つ以上のエコー画像（即ち各エコー時間に起因する１つのＭＲ画像）を取得することによって脂肪及び／又は水の位相感応型ＭＲ再構成を行うために、ディクソン型水／脂肪撮像の一般的な概念が採用される。本発明の手法は、水と脂肪との複素画像の分離を可能にし、ここで、水及び／又は脂肪画像の位相感応型再構成が行われる。この再構成は、手動入力が必要とされないように、完全に自動化されて実施され得る。追加の、例えば脂肪抑制された参照水画像の取得が避けられ得る。

本発明の好ましい実施形態では、信号寄与の分離は、主磁場Ｂ_０の空間的変化に従う、それぞれのエコー時間にわたる少なくとも２種の化学種の信号寄与の位相推移の予測と、（ｋ空間内又は画像空間内で）取得されたエコー信号の対応する復調とを含む。多くの既知のディクソン型の用途で一般的であるように、主磁場分布が空間にわたって滑らかに変化するという仮定に基づいて、化学シフトによって引き起こされる位相推移が、主磁場分布によって引き起こされる位相推移から分離され得る。

本発明に関連して、用語「化学種」は、ＭＲ特性を有する任意の種類の化学物質又は任意の種類の核として広範に解釈されなければならないことに留意されたい。単純な例では、２種の化学種のＭＲ信号が取得され、ここで、化学種は「化学組成物」である水及び脂肪中のプロトンである。より洗練された例では、マルチピークスペクトルモデルが、既知の相対量で生じる異なる化学組成物の組における核を実際に表す。

本発明の好ましい実施形態では、化学種の少なくとも１種の信号寄与から位相感応型ＭＲ画像を再構成するステップが、（ｉ）極性とバックグランドフェージングとの区別、及び（ｉｉ）バックグランドフェージングに従ったそれぞれの化学種の信号寄与の位相補正を含む。これらのステップは、（画像空間内又はｋ空間内で）取得されたエコー信号からバックグランドフェージングを決定するために使用される。この目的で、バックグランドフェージングの推定のために領域拡張技法を使用するそれ自体知られているＲＡＰＩＤ方法（上記参照）が採用され得る。ＲＡＰＩＤアルゴリズムは、別途取得される参照画像に依拠することなく、磁化の極性を確実に復元することが可能である。ディクソンアルゴリズム（上記参照）を使用してエコー信号から導出されたＢ_０マップは、有利には、位相感応型ＭＲ画像を再構成するステップにおいて、それぞれの化学種の信号寄与のバックグランドフェージングを予測するために使用され得る。即ち、水／脂肪分離ステップにより得られるＢ_０マップは、有利には、本発明に従ってＭＲ画像の位相感応型再構成を支援するために使用される。

上述した本発明の方法は、検査ボリューム内部で均一な定常磁場Ｂ_０を発生させるための少なくとも１つの主磁石コイルと、検査ボリューム内部で異なる空間方向でスイッチ磁場勾配を発生させるための複数の勾配コイルと、検査ボリューム内部でＲＦパルスを発生させるため及び／又は検査ボリューム内に位置付けされた患者の身体からＭＲ信号を受信するための少なくとも１つの身体ＲＦコイルと、時間的に連続するＲＦパルス及びスイッチ磁場勾配を制御するための制御ユニットと、受信されたＭＲ信号からＭＲ画像を再構成するための再構成ユニットとを含む、ＭＲデバイスによって実施され得る。本発明の方法は、ＭＲデバイスの再構成ユニット及び／又は制御ユニットの対応するプログラミングによって実装され得る。

本発明の方法は、有利には、現在臨床で使用されているほとんどのＭＲデバイスで実施され得る。このために、本発明の上述した方法ステップを実施するようにＭＲデバイスを制御するコンピュータプログラムを利用することのみが必要である。コンピュータプログラムは、データキャリア上に存在していてよく、又はＭＲデバイスの制御ユニットにインストールするためにダウンロードされるようにデータネットワーク上に存在していてもよい。

添付図面は、本発明の好ましい実施形態を開示する。しかし、図面は、例示の目的を意図されているに過ぎず、本発明の範囲を定義するものとは意図されていないことを理解すべきである。

本発明の方法を実施するためのＭＲデバイスを示す。 本発明の方法を例示するフローチャートを示す。 本発明の方法に従って取得されたＭＲ画像を示す。

図１を参照すると、ＭＲデバイス１が示されている。このデバイスは、超伝導性又は抵抗性の主磁石コイル２を備え、これにより、実質的に均一であり時間的に一定の主磁場Ｂ_０が、検査ボリュームを通ってｚ軸に沿って生成される。このデバイスは、（１次、２次、及び適用可能であれば３次）シムコイル２’の組を更に備え、その組の個々のシムコイル２’を通る電流は、検査ボリューム内部でのＢ_０偏差を最小限にする目的で制御可能である。

磁気共鳴発生及び操作システムは、一連のＲＦパルスとスイッチ磁場勾配とを印加して、核磁気スピンの反転又は励起、磁気共鳴の誘発、磁気共鳴の再合焦、磁気共鳴の操作、磁気共鳴の空間的及びその他のエンコード、スピンの飽和等を行って、ＭＲ撮像を実施する。

より具体的には、勾配パルス増幅器３が、検査ボリュームのｘ、ｙ、及びｚ軸に沿って、全身勾配コイル４、５、及び６のうちの選択されたコイルに電流パルスを印加する。デジタルＲＦ周波数送信機７が、ＲＦパルス又はパルスパケットを、送受信スイッチ８を介して身体ＲＦコイル９に送信して、ＲＦパルスを検査ボリューム内に送信する。典型的なＭＲ撮像シーケンスは、短い持続時間のＲＦパルスセグメントのパケットから構成され、これは、任意の印加される磁場勾配と共に、核磁気共鳴の選択された操作を実現する。ＲＦパルスは、飽和、共鳴の励起、磁化の反転、共鳴の再合焦、又は共鳴の操作、及び検査ボリューム内に位置付けされた身体１０の一部分の選択を行うために使用される。また、ＭＲ信号が、身体ＲＦコイル９によってピックアップされる。

身体１０の限定された領域のＭＲ画像を生成するために、局所アレイＲＦコイル１１、１２、１３の組が、撮像のために選択された領域に隣接して配置される。アレイコイル１１、１２、１３は、身体コイルＲＦ送信によって誘発されたＭＲ信号を受信するために使用され得る。

生じたＭＲ信号は、身体ＲＦコイル９によって及び／又はアレイＲＦコイル１１、１２、１３によってピックアップされ、好ましくは前置増幅器（図示せず）を含む受信機１４によって復調される。受信機１４は、送受信スイッチ８を介してＲＦコイル９、１１、１２、及び１３に接続される。

ホストコンピュータ１５が、シムコイル２’、並びに勾配パルス増幅器３及び送信機７を制御して、本発明によるＭＲ撮像シーケンスをもたらす。受信機１４は、各ＲＦ励起パルスに続いて直ちに、単一の又は複数のＭＲデータラインを受信する。データ取得システム１６は、受信された信号のアナログ／デジタル変換を行い、各ＭＲデータラインを、更なる処理に適したデジタルフォーマットに変換する。現代のＭＲデバイスでは、データ取得システム１６は、生の画像データの取得に特化された独立したコンピュータである。

最終的に、生のデジタル画像データは、再構成処理装置１７によって再構成されて画像表現となり、再構成処理装置１７は、フーリエ変換又は他の適切な再構成アルゴリズムを適用する。ＭＲ画像は、患者を通る平面状のスライス、平行な平面状スライスのアレイ、３次元ボリューム等を表現することがある。次いで画像は画像メモリに記憶され、そこで画像はアクセスされ得、スライス、投影、又は画像表現の他の部分を、例えばビデオモニタ１８による視覚化のための適切なフォーマットに変換する。ビデオモニタ１８は、得られたＭＲ画像の人間可読の表示を提供する。

本発明の例示的実施形態では、第１及び第２のエコー信号の発生のために反転回復デュアル勾配エコーシーケンスが使用され、第１のエコー（ＴＥ１）のエコー時間は１．８ｍｓであり、第２のエコーのエコー時間（ＴＥ２）は３．１ｍｓである。１．５テスラの主磁場において、第１のエコー信号に対する水スピンと脂肪スピンとの寄与は、より位相がずれたものとなり、一方、第２のエコー信号の時点ではより位相が合ったものとなる。望ましい視野の完全なＭＲ画像を再構成することが可能であるように、一般的な方法で適切な位相エンコードを用いて、複数の第１及び第２のエコー信号が発生及び取得される。撮像シーケンスは、最初の反転ＲＦパルスと、それに続く励起ＲＦパルスとを含む。最初の反転ＲＦパルス（１８０°）は、主磁場に沿った縦磁場を−ｚ方向に回転させる。次いで、反転された磁化は、反転ＲＦパルスと励起ＲＦパルスとの間の反転回復時間（ＴＩ）中にＴ_１緩和によって回復される。励起ＲＦパルスの後に、適切な勾配切り替えによってエコー信号が発生される。心筋梗塞の診断及び定量化のための典型的な用途において、ＴＩの持続時間は、２３０〜２７０ｍｓの範囲内でよい。

図２は、それぞれＴＥ１及びＴＥ２で第１及び第２のエコー信号から再構成される２つの単一エコー画像２１、２２を概略的に示す。取得されたエコー信号データ２１、２２に対する水と脂肪との信号寄与を分離するために、一般的な２点ディクソンアルゴリズム２３（例えば、EggersらによってMagnetic Resonance in Medicine, 65, 96-107, 2011で述べられているもの等）が採用される。その結果、水の（複素）中間画像２４と脂肪の（複素）中間画像２５とが得られる。更に、ディクソンアルゴリズムは、Ｂ_０磁場不均一性のマップ、即ちＢ_０マップ２６を提供する。次のステップで、それぞれ中間画像２４及び２５からの水画像２８及び脂肪画像２９の位相感応型再構成のために、上述したＲＡＰＩＤアルゴリズム２７が採用される。エコー信号２１、２２から導出されるＢ_０マップ２６は、有利には、「真の位相（true phase）」ＭＲ画像２８、２９の再構成における水と脂肪との信号寄与のバックグランドフェージングを予測するために使用される。このようにして、Ｂ_０マップ２６は、位相感応型再構成２７を支援するために使用される。

上述した技法は、例えば、診断の信頼性を改良するために脂肪（トリグリセリド沈着）を心筋瘢痕から分離するための位相感応型ガドリニウム遅延造影（ＬＧＥ：late gadolinium enhancement）で利用することができる。また、位相感応型であり且つロバストな水／脂肪分離が必要とされる他の用途にも使用され得る。

図３は試料画像を示し、ここでは、本発明の方法が、身体撮像における２点ディクソンＰＳＩＲ取得に適用されている（ＴＩ＝５００ｍｓ）。図３ａは水画像であり、図３ｂは脂肪画像である。図３ａ及び図３ｂの画像は、本発明に従って分離されている。図３ｃに示されるような位相感応型の水画像が、追加の位相参照スキャンなしで、図３ａの画像から再構成される。図３ｃの位相感応型の水画像で見られ得るように、異なる臓器が異なる信号極性を表す。即ち、肝臓組織（血管を除く）３１が正の信号極性を有し、腎臓組織３２が負の磁化を示す。このようにして、かなり強調された画像コントラストが実現される。それぞれの臓器の異なるＴ_１値によって、非常に異なる信号極性が引き起こされる。これらは、本発明のアルゴリズムによって適正に識別される。

Claims (8)

1. 異なるＭＲスペクトルを有する少なくとも２種の化学種をＭＲ撮像する方法であって、
   ａ）ＭＲデバイスの検査ボリューム内に位置付けされた対象物を、ＲＦパルスとスイッチ磁場勾配との撮像シーケンスにさらすことによって、２つ以上の異なるエコー時間にエコー信号を発生させるステップであって、前記撮像シーケンスが、反転ＲＦパルスと、それに続く、反転回復時間後の励起ＲＦパルスとを含む反転回復シーケンスである、ステップと、
   ｂ）前記エコー信号を取得するステップと、
   ｃ）取得された前記エコー信号に対する前記少なくとも２種の化学種の信号寄与を分離するステップと、
   ｄ）前記化学種の少なくとも１種の前記信号寄与から、位相感応型ＭＲ画像を再構成するステップと
   を含み、
   ｅ）前記化学種の少なくとも１種の前記信号寄与から前記位相感応型ＭＲ画像を再構成する前記ステップが、（ｉ）極性とバックグランドフェージングとを区別すること、及び（ｉｉ）前記バックグランドフェージングに従った前記それぞれの化学種の前記信号寄与の位相補正を行うことを含む、
   方法。
2. ステップｃ）での前記信号寄与の前記分離が、主磁場Ｂ０の少なくとも空間的変化に従った、それぞれの前記エコー時間にわたる前記少なくとも２種の化学種の前記信号寄与の位相推移の予測を含む、請求項１に記載の方法。
3. 前記それぞれのエコー時間にわたる前記少なくとも２種の化学種の前記信号寄与の前記位相推移の前記予測が、前記信号寄与の位相を復調するために使用される、請求項２に記載の方法。
4. 前記少なくとも２種の化学種の前記信号寄与を分離する前記ステップにおいてディクソン技法が使用される、請求項１乃至３の何れか一項に記載の方法。
5. 前記バックグランドフェージングの推定のために領域拡張技法が使用される、請求項１に記載の方法。
6. Ｂ０マップが、ディクソンアルゴリズムを使用して前記エコー信号から導出され、前記Ｂ０マップが、前記位相感応型ＭＲ画像を再構成する前記ステップにおいて、前記それぞれの化学種の前記信号寄与の前記バックグランドフェージングを予測するために使用される、請求項１又は５に記載の方法。
7. 異なるＭＲスペクトルを有する少なくとも２種の化学種をＭＲ撮像するためのＭＲデバイスであって、
   検査ボリューム内部で均一な定常磁場Ｂ０を発生させるための少なくとも１つの主磁石コイルと、
   前記検査ボリューム内部で異なる空間方向でスイッチ磁場勾配を発生させるための複数の勾配コイルと、
   前記検査ボリューム内部でＲＦパルスを発生させるため及び／又は前記検査ボリューム内に位置付けされた対象物からＭＲ信号を受信するための少なくとも１つのＲＦコイルと、
   時間的に連続するＲＦパルス及びスイッチ磁場勾配を制御するための制御ユニットと、
   前記受信されたＭＲ信号からＭＲ画像を再構成するための再構成ユニットとを含み、
   ａ）ＭＲデバイスの検査ボリューム内に位置付けされた対象物を、ＲＦパルスとスイッチ磁場勾配との撮像シーケンスにさらすことによって、第１のエコー時間に一連の第１のエコー信号を発生させ、且つ第２のエコー時間に一連の第２のエコー信号を発生させることであって、前記撮像シーケンスが、反転ＲＦパルスと、それに続く、反転回復時間後の励起ＲＦパルスとを含む反転回復シーケンスであることと、
   ｂ）前記エコー信号を取得することと、
   ｃ）取得された前記エコー信号に対する前記少なくとも２種の化学種の信号寄与を分離することと、
   ｄ）前記化学種の少なくとも１種の前記信号寄与から、位相感応型ＭＲ画像を再構成することと、
   ｅ）前記化学種の少なくとも１種の前記信号寄与から前記位相感応型ＭＲ画像を再構成することにおいて、（ｉ）極性とバックグランドフェージングとを区別すること、及び（ｉｉ）前記バックグランドフェージングに従った前記それぞれの化学種の前記信号寄与の位相補正を行うことと
   を実施する、ＭＲデバイス。
8. ＭＲデバイス上で実行されるコンピュータプログラムであって、
   ａ）ＲＦパルスとスイッチ磁場勾配との撮像シーケンスであって、反転ＲＦパルスと、それに続く、反転回復時間後の励起ＲＦパルスとを含む反転回復シーケンスである撮像シーケンスによって、第１のエコー時間に一連の第１のエコー信号を発生させ、且つ第２のエコー時間に一連の第２のエコー信号を発生させ、
   ｂ）前記エコー信号を取得し、
   ｃ）取得された前記エコー信号に対する前記少なくとも２種の化学種の信号寄与を分離し、
   ｄ）前記化学種の少なくとも１種の前記信号寄与から、位相感応型ＭＲ画像を再構成し、及び
   ｅ）前記化学種の少なくとも１種の前記信号寄与から前記位相感応型ＭＲ画像を再構成する際に、（ｉ）極性とバックグランドフェージングとを区別すること、及び（ｉｉ）前記バックグランドフェージングに従った前記それぞれの化学種の前記信号寄与の位相補正を行う
   ための命令を含む、コンピュータプログラム。

Images