(57)【特許請求の範囲】
【請求項１】
少なくとも２つの化学種を、該２つの化学種に関する信号分離を用いた磁気共鳴撮像を用いて撮像する方法であって：
−  相異なるエコー時間にて第１及び第２のエコーデータを収集して、第１及び第２の収集複素データセットを得るステップと、
−  前記第１及び第２の収集データセットをモデル化するステップであり、該モデル化は前記化学種のうちの少なくとも１つの化学種のスペクトル信号モデルを有するステップと、
−  前記第１及び第２の収集データセット内で、前記モデル化が前記信号分離に関して単一の一義的な数学解を生み出すボクセルを特定するステップと、
−  前記モデル化が２つ以上の数学解を生み出すボクセルが残っている場合に、該ボクセルの単一の数学解を決定するステップと、
を有し、
当該方法は更に、前記相異なるエコー時間を予め決定するステップを有し、該予め決定するステップは：
−  前記化学種のうちの少なくとも１つの化学種の適切なスペクトル信号モデルを選択あるいは校正するステップと、
−  時変的な複素重みによって記述される、前記化学種のうちの少なくとも１つの化学種の信号のエコー時間にわたる振幅及び位相の変化を計算するステップと、
−  これら複素重みが所定の制約を満足するエコー時間であり、それらエコー時間によって前記信号分離において前記モデル化が２つ以上の数学解を生み出すボクセルの個数が最小化されるエコー時間、を選択するステップと、
−  前記第１及び第２のエコーデータを収集して前記第１及び第２の収集複素データセットを得るために、選択されたエコー時間を使用するステップと
を有する、
方法。
【請求項２】
前記特定が為されていないボクセルに対し、先ず、前記モデル化によって、前記信号分離に関する全ての数学解が計算され、その後、前記特定が為されていないボクセルの直近に位置する前記特定が為されたボクセルに基づいて、１つの数学解が選択される、請求項１に記載の方法。
【請求項３】
前記第１及び第２の収集データセットの前記モデル化は、各ボクセルに関する第１及び第２の位相誤差を有し、１つの数学解の前記選択は、前記特定が為されていないボクセルの前記直近における、前記第１及び第２の位相誤差の、あるいは前記第１の位相誤差と前記第２の位相誤差との間の差の、空間的に滑らかな変化を仮定して行われる、請求項２に記載の方法。
【請求項４】
第１の化学種からの信号の振幅が第１のエコー時間から第２のエコー時間までに減衰するよりも大きく、第２の化学種からの信号の振幅が前記第１のエコー時間から前記第２のエコー時間までに減衰し、前記所定の制約は、
−  主として前記第２の化学種が識別される場合に、前記２つの化学種からの信号が、前記第１のエコー時間のときよりも、前記第２のエコー時間のときの方がより同相であるように、前記第１及び第２のエコー時間を選択し、
−  主として前記第１の化学種が識別される場合に、前記２つの化学種からの信号が、前記第２のエコー時間のときよりも、前記第１のエコー時間のときの方がより同相であるように、前記第１及び第２のエコー時間を選択し、
−  前記第１のエコー時間から前記第２のエコー時間までの前記第２の化学種からの信号の振幅の変化が最大化されるように、前記相異なるエコー時間を適応する、
ことを有する、請求項１に記載の方法。
【請求項５】
前記第１及び第２の収集データセット内での、前記モデル化が前記信号分離に関して単一の一義的な数学解を生み出すボクセルの前記特定を、予期されるＳＮＲ及び緩和による信号減衰に基づく範囲に制限すること、を更に有する請求項１に記載の方法。
【請求項６】
前記２つの化学種は水と脂肪であり、前記スペクトル信号モデルは脂肪のマルチピークスペクトルモデルである、請求項１に記載の方法。
【請求項７】
双方の化学種が識別される場合に、前記２つの化学種からの信号が双方のエコー時間において同相に近くなるように第１及び第２のエコー時間を選択する、請求項１に記載の方法。
【請求項８】
請求項１乃至７の何れかに記載の方法を実行するためのコンピュータ実行可能命令を有するコンピュータプログラム。
【請求項９】
少なくとも２つの化学種を撮像する磁気共鳴撮像装置であって、当該装置は、磁気共鳴画像データを収集する磁気共鳴撮像スキャナを有し、該スキャナは：
−  相異なるエコー時間にて第１及び第２のエコーデータを収集して、第１及び第２の収集複素データセットを取得し、
−  前記第１及び第２の収集データセットをモデル化し、該モデル化は前記化学種のうちの少なくとも１つの化学種のスペクトル信号モデルを有し、
−  前記第１及び第２の収集データセット内で、前記モデル化が信号分離に関して単一の一義的な数学解を生み出すボクセルを特定し、且つ
−  前記モデル化が２つ以上の数学解を生み出すボクセルが残っている場合に、該ボクセルの単一の数学解を決定する、
ように動作可能であり、
前記スキャナは更に、前記相異なるエコー時間を事前決定するように動作可能であり、該事前決定は：
−  前記化学種のうちの少なくとも１つの化学種の適切なスペクトル信号モデルを選択あるいは校正し、
−  時変的な複素重みによって記述される、前記化学種のうちの少なくとも１つの化学種の信号のエコー時間にわたる振幅及び位相の変化を計算し、
−  これら複素重みが所定の制約を満足するエコー時間であり、それらエコー時間によって前記信号分離において前記モデル化が２つ以上の数学解を生み出すボクセルの個数が最小化されるエコー時間、を選択し、且つ
−  前記第１及び第２のエコーデータを収集して前記第１及び第２の収集複素データセットを得るために、選択されたエコー時間を使用する
ことを有する、
磁気共鳴撮像装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、２つの化学種に対する信号分離を用いた磁気共鳴撮像を用いて少なくとも２つの化学種を撮像する方法、少なくとも２つの化学種を撮像するためのコンピュータプログラム及び磁気共鳴撮像装置に関する。
【背景技術】
【０００２】
磁場と核スピンとの間の相互作用を利用して２次元又は３次元の画像を形成する画像形成ＭＲ法が、特に医療診断の分野において、今日広く使用されている。何故なら、画像形成ＭＲ法は、軟組織の撮像において数多くの点で他の撮像法より優れており、電離放射線を必要とせず、且つ通常、非侵襲的であるからである。
【０００３】
ＭＲ法によれば一般的に、患者の体又は一般的には検査対象が、強く均一な磁場Ｂ０内に配置される。同時に、この磁場の方向が、測定が基礎とする座標系の軸（通常はｚ軸）を定める。
【０００４】
この磁場は、印加される磁場強度に応じて個々の核スピンに異なるエネルギーレベルを生成する。これらのスピンは、所謂ラーモア周波数又はＭＲ周波数である規定の周波数の交番磁場（ＲＦ場）の印加によって、励起（スピン共鳴）されることができる。巨視的な視点から、個々の核スピンの分布が全体的な磁化を生成し、該磁化は、磁場がｚ軸に垂直に広がりながら適切な周波数の電磁パルス（ＲＦパルス）が印加されることによって平衡状態から逸らされることが可能であり、結果として、該磁化はｚ軸の周りで歳差運動を行う。
【０００５】
磁化の如何なる変化も、受信ＲＦアンテナによって検出されることができる。受信ＲＦアンテナは、ｚ軸に垂直な方向での磁化の変化が測定されるように、ＭＲ装置の検査ボリューム（容積）内で方向付けて配置される。
【０００６】
人体内で空間分解能を実現するため、スピン共鳴周波数の線形な空間依存性をもたらすよう、３つの主軸に沿って延在する一定の磁場勾配が上記均一な磁場に重畳される。その結果、受信アンテナにて拾われる信号は、人体内の異なる位置に関連付けられることが可能な異なる周波数の成分を含むものとなる。
【０００７】
受信アンテナによって得られる信号データは、空間周波数ドメインに対応するものであり、ｋ空間データと呼ばれる。ｋ空間データは通常、異なる位相エンコーディングで収集された複数の線（ライン）を含む。多数のサンプルを収集することによって各ラインがデジタル化される。ｋ空間データの一組のサンプル（サンプルセット）が、例えばフーリエ変換によって、１枚のＭＲ画像に変換される。
【０００８】
ＭＲＩにおいては、しばしば、例えば水と脂肪などの２つの支配的な化学種について、それらの一方の寄与を抑制するために、あるいはそれら双方の寄与を別々若しくは一緒に分析するために、全体信号に対するそれら２つの化学種の相対的な寄与について情報を取得することが望ましい。それらの寄与分は、異なるエコー時間にて収集される２つ以上の対応するエコーからの情報が結合される場合に計算されることができる。
【０００９】
ＭＲ信号に対する水の寄与及び脂肪の寄与についての情報を同時に取得する一手法は化学（ケミカル）シフトエンコーディングである。化学シフトエンコーディングにおいては、化学シフト次元という追加の次元が定義され、僅かに異なるエコー時間にて２つの画像を収集することによってエンコードされる。
【００１０】
特に水−脂肪分離の場合、このような種類の実験はしばしば、ディクソン（Ｄｉｘｏｎ）型の測定と呼ばれている。Ｄｉｘｏｎ撮像又はＤｉｘｏｎ水／脂肪分離撮像により、異なるエコー時間にて収集される２つ以上の対応するエコーから水及び脂肪の寄与分を計算することによって、水−脂肪分離を得ることができる。Ｄｉｘｏｎ撮像は通常、水信号と脂肪信号とを分離するために少なくとも２つのエコーを収集することに頼っている。一般に、このような種類の分離は、脂肪内と水内とでの水素の既知の歳差運動周波数差の存在のおかげで可能である。その最も単純な形態においては、‘同相（インフェーズ）’データセット及び‘異相（アウトオブフェーズ）’データセットの加算又は減算の何れかによって水画像及び脂肪画像が生成されるが、この手法は主磁場の均一性の影響をかなり受けやすい。
【００１１】
水画像内に残留脂肪信号が存在しない高品質の水−脂肪分離は、脂肪スペクトルの複雑なモデルが水−脂肪分離プロセスに組み込まれる場合に得られ得る。これは、例えば非特許文献１にて三点Ｄｉｘｏｎ法に関して例証されている。
【００１２】
特に、例えば一度の息止めでの腹部撮像などの時間が重視される用途では、スキャン時間を可能な限り短縮するために二点法が好ましく使用される。しかしながら、二点法は、単一の支配的なピークによって脂肪スペクトルを近似するものであり、故に一般的に、より実効的な脂肪抑制を提供しない。また、脂肪抑制の品質は、画像データ収集におけるエコー時間の選択に大きく依存するものとなる。
【００１３】
水信号と脂肪信号との識別は、結果画像の正確なラベル付け又は表示にとって望ましいだけでなく、より良好な初期設定の恩恵を実質的に受けるものである水信号と脂肪信号とのロバストな分離にとっても望ましい。しかしながら、化学シフトと主磁場の不均一性との間の生来的な不明確性のため、これは難題である。二点Ｄｉｘｏｎ撮像に関し、２つの信号間の先行あるいは遅行する位相関係を活用する部分的逆位相収集に基づく水信号と脂肪信号との識別が提案されたことがある（非特許文献２）。この識別は、しかし、水及び脂肪の双方を含んだボクセル群のうちの通常はごく少数に制限されたものであった。
【００１４】
多点Ｄｉｘｏｎ撮像に関しては、分離に脂肪のシングルピークスペクトルモデル及びマルチピークスペクトルモデルを用いて残る残余の比較に基づく水信号と脂肪信号との識別が提案されたことがある（非特許文献３）。しかしながら、この手法は、より多点の場合にのみロバストである。
【先行技術文献】
【非特許文献】
【００１５】
【非特許文献１】Yu  H，Shimakawa  A，McKenzie  CA，Brodsky  E，Brittain  JH，Reeder  SB，“Multi-echo  water-fat  separation  and  simultaneous  R2＊  estimation  with  multi-frequency  fat  spectrum  modeling”，Magn  Reson  Med，2008；60：1122-1134
【非特許文献２】Xiang  QS，“Two-point  water-fat  imaging  with  partially-opposed-phase  (POP)  acquisition：an  asymmetric  Dixon  method”，Magn  Reson  Med，2006；56：572-584
【非特許文献３】Yu  H，Shimakawa  A，Brittain  JH，McKenzie  CA，Reeder  SB，“Exploiting  the  spectral  complexity  of  fat  for  robust  multi-point  water-fat  separation”，Proc  ISMRM  2010；771
【発明の概要】
【発明が解決しようとする課題】
【００１６】
以上から直ちに認識されるように、改善されたＭＲ撮像法が望まれる。従って、本発明の１つの目的は、高速且つ信頼性高く２つの支配的な化学種を高品質分離して、収集された全体信号に対するそれら２つの支配的化学種の相対寄与を決定するＭＲ撮像を実現することである。
【課題を解決するための手段】
【００１７】
本発明によれば、少なくとも２つの化学種を、該２つの化学種に関する信号分離を用いた磁気共鳴撮像を用いて撮像する方法が提供される。当該方法は、
−  相異なるエコー時間にて第１及び第２のエコーデータを収集して、第１及び第２の収集複素データセットを得るステップと、
−  前記第１及び第２の収集データセットをモデル化するステップであり、該モデル化は前記化学種のうちの少なくとも１つの化学種のスペクトル信号モデルを有するステップと、
−  前記第１及び第２の収集データセット内で、前記モデル化が前記信号分離に関して単一の一義的な数学解を生み出すボクセルを特定するステップと、
−  前記モデル化が２つ以上の数学解を生み出すボクセルが残っている場合に、該ボクセルの不明確性を解消するステップと、を有する。
【００１８】
換言すれば、本発明は、二点法による分離に、少なくとも２つの化学種のうちの少なくとも１つの化学種のスペクトルの、より複雑なモデルを組み入れるという洞察を用いる。そうすることは、２つの化学種の信号分離の精度と、化学種の抑制（サプレッション）の効率とを高めることを可能にする。このアプローチから始まって、本発明の実施形態は、複数の化学種の信号の一層ロバストな分離を、これら信号へのこれら化学種の正確な割付けにより可能にする。
【００１９】
なお、例えばマルチピークスペクトルモデルによって、化学種のうちのただ１つがモデル化される場合、その他の化学種は単に単一のラインと見なされ得る。従って、実際には、双方の化学種がモデル化され、それらモデルのうちの少なくとも１つがマルチピークスペクトルモデルを有する。
【００２０】
また、用語‘化学種’は、任意の種類の所定の化学物質、又は所定の磁気共鳴特性の任意の種類の核として、広く理解されなければならない。単純な一例において、２つの化学種は‘化学成分’水及び脂肪内のプロトンである。より洗練された一例において、マルチピークスペクトルモデルは実際に、既知の相対量にて生じる一組の‘化学成分’内の或る核を記述する。この場合、マルチピークスペクトルモデルは、例えば代謝過程といった化学反応の生成物などの、２つの信号成分を分離するために導入される。故に、特定の代謝過程の出力が、相異なる生成物の出力比に関してモデル化され得る。そして、相異なる化学種ではあるが、これらは単一の信号成分へと結合され得る。
【００２１】
また、ここでは明示的にはデュアルエコーデータ収集のみに言及するが、本提案に係る方法は、マルチエコーデータ収集の好適なサブセットにも等しく適用可能である。
【００２２】
本発明の一実施形態によれば、当該方法は更に、相異なるエコー時間を予め決定するステップを有し、該予め決定するステップは：
−  前記化学種のうちの少なくとも１つの化学種の適切なスペクトル信号モデルを選択あるいは校正するステップと、
−  時変的な複素重みによって記述される、前記化学種のうちの少なくとも１つの化学種の信号のエコー時間にわたる振幅及び位相の変化を計算するステップと、
−  これら複素重みが所定の制約を満足するエコー時間であり、それらエコー時間によって信号分離において前記モデル化が２つ以上の数学解を生み出すボクセルの個数が最小化されるエコー時間、を選択するステップと、
−  第１及び第２のエコーデータを収集して第１及び第２の収集複素データセットを得るために、選択されたエコー時間を使用するステップと、を有する。
【００２３】
これは、化学種の信号が一義的に識別されるボクセルの個数が最大化されるという利点を有する。故に、第１及び第２のエコーデータを用いて収集された信号の化学種への割付けが単純化され、ひいては、分離の堅牢性が高められる。
【００２４】
本発明の一実施形態によれば、前記特定が為されていないボクセルに対し、先ず、前記モデル化によって、信号分離に関する全ての数学解が計算され、その後、前記特定が為されていないボクセルの直近に位置する前記特定が為されたボクセルに基づいて、１つの数学解が選択される。例えば、第１及び第２の収集データセットの前記モデル化は、各ボクセルに関する第１及び第２の位相誤差を有し、１つの数学解の前記選択は、前記特定が為されていないボクセルの前記直近における、第１及び第２の位相誤差の、あるいは第１の位相誤差と第２の位相誤差との間の差の、空間的に滑らかな変化を仮定して行われる。
【００２５】
本発明の一実施形態によれば、第１の化学種からの信号の振幅が第１のエコー時間から第２のエコー時間までに減衰するより大きく、第２の化学種からの信号の振幅が第１のエコー時間から第２のエコー時間までに減衰し、前記所定の制約は、
−  主として第２の化学種が識別される場合に、２つの化学種からの信号が、第１のエコー時間において、第２のエコー時間においてよりも同相でないように、第１及び第２のエコー時間を選択し、
−  主として第１の化学種が識別される場合に、２つの化学種からの信号が、第１のエコー時間において、第２のエコー時間においてよりも同相であるように、第１及び第２のエコー時間を選択し、
−  第１のエコー時間から第２のエコー時間までの信号の振幅の相対的な低下が最大化されるように、前記相異なるエコー時間を適応する、ことを有する。
【００２６】
これは、化学種の信号が一義的に識別されるボクセルの個数が最大化されること、前もって確保することを可能にする。
【００２７】
本発明の一実施形態によれば、２つの化学種は水と脂肪であり、スペクトル信号モデルは脂肪のマルチピークスペクトルモデルである。従って、本発明は、スペクトル内に唯一の支配的な脂肪のスペクトルピークのみが存在すると単純に仮定するのではなく、脂肪のマルチピークスペクトルモデルを採用する。例えば脂肪である化学種のうちの１つに関し、相対的な共鳴周波数及び相対的な共鳴強さが前もって既知であり、例えば理論モデル又は実験モデルから、あるいは別個若しくは統合された校正（例えば、脂肪のような１つの化学種のみを含む画素の識別に基づく）から得られていると仮定する。
【００２８】
例えば脂肪である化学種のうちの１つのみがモデル化される場合、水はシングルピークスペクトルと見なされ得る。
【００２９】
２つの化学種が水と脂肪である場合、上述の所定の制約は、
−  水信号及び脂肪信号が第２のエコー時間においてよりも同相でないような第１のエコー時間を選択し、これは一般に、純粋な脂肪信号の一義的な識別を可能にし、
−  水信号及び脂肪信号が第２のエコー時間においてよりも同相であるような第１のエコー時間を選択し、これは一般に、純粋な水信号、又は更には水が支配的な信号、の一義的な識別を可能にし、
−  （使用するスペクトルモデルに従って）これら２つのエコー時間の間に純粋な脂肪信号の振幅が可能な限り大きく変化するようにこれらエコー時間を適応する、ことを有する。
【００３０】
本発明の更なる一実施形態によれば、当該方法は更に、第１及び第２の収集データセット内での、前記モデル化が信号分離に関して単一の一義的な数学解を生み出すボクセルの前記特定を、予期されるＳＮＲ及び緩和による信号減衰に基づく範囲に制限することを有する。故に、乏しいＳＮＲの場合、単一の数学解を有する信号領域内に位置しながら、２つ以上の数学解を有する信号領域の近くにあるボクセルが、２つの数学解を検討しなければならないボクセルとして扱われる。これにより、そのボクセルに関して２つの化学種の信号を識別するときに実際の数学解がうっかり排除されてしまう、ことがないことが確保される。
【００３１】
本発明の一実施形態によれば、記載のデータセットは、画像空間又はｋ空間のデータセットとし得る。画像データセットが好ましい場合、第１及び第２のエコーデータは、第１及び第２の収集データセットの再構成のために処理され得る。この場合、第１及び第２の収集データセットは画像データセットとなる。
【００３２】
本発明に係る方法は有利なことに、現在臨床的に使用されている大抵のＭＲ装置で実行されることが可能である。この目的のために、ＭＲ装置が本発明に係る上述の方法ステップを実行するようにＭＲ装置を制御するコンピュータプログラムを使用することを必要とするのみである。コンピュータプログラムは、データキャリア上に存在してもよいし、ＭＲ装置の制御ユニットへのインストールのためにダウンロードされるようにデータネットワーク内に存在してもよい。故に、本発明はまた、上述の方法を実行するためのコンピュータ実行可能命令を有するコンピュータプログラム製品にも関する。
【００３３】
さらに、本発明は、少なくとも２つの化学種を撮像する磁気共鳴撮像装置に関する。当該装置は、磁気共鳴画像データを収集する磁気共鳴撮像スキャナを有し、該スキャナは：
−  相異なるエコー時間にて第１及び第２のエコーデータを収集して、第１及び第２の収集複素データセットを取得し、
−  第１及び第２の収集データセットをモデル化し、該モデル化は前記化学種のうちの少なくとも１つの化学種のスペクトル信号モデルを有し、
−  第１及び第２の収集データセット内で、前記モデル化が信号分離に関して単一の一義的な数学解を生み出すボクセルを特定し、且つ
−  前記モデル化が２つ以上の数学解を生み出すボクセルが残っている場合に、該ボクセルの不明確性を解消する、ように動作可能である。
【図面の簡単な説明】
【００３４】
図面は、本発明の好適実施形態を開示するものである。しかしながら、理解されるように、図面は、単に例示目的で作成されたものであり、本発明の限定を定めるものとして作成されたものではない。
【図１】ＭＲ撮像システムを模式的に示す図である。
【図２】複素因子ｃの時間変動を例示する図である。
【図３】水信号及び脂肪信号の解を示す図である。
【図４】２つのエコー時間での信号振幅の比の分類を例示する図である。
【発明を実施するための形態】
【００３５】
図１を参照するに、ＭＲ撮像システム１が示されている。このシステムは、超電導又は抵抗性の主磁石コイル２を有しており、それにより、検査ボリュームを通るｚ軸に沿って、実質的に均一で時間的に一定の主磁場Ｂ０が作り出される。
【００３６】
ＭＲ撮像を行うため、磁気共鳴生成操作システムが、一連のＲＦパルス及び切換式の磁場勾配を印加して、核磁気スピンを反転あるいは励起し、磁気共鳴を誘起し、磁気共鳴をリフォーカスし、磁気共鳴を操作し、磁気共鳴を空間的あるいはその他手法でエンコードし、スピンを飽和させ、などする。
【００３７】
より具体的には、傾斜パルス増幅器３が、検査ボリュームのｘ軸、ｙ軸及びｚ軸に沿う全身傾斜コイル４、５及び６のうちの選択されたものに電流を与える。ＲＦ送信器７が送／受スイッチ９を介してＲＦアンテナ９にＲＦパルス又はパルスパケットを送ることで、ＲＦパルスが検査ボリューム内に送信される。典型的なＭＲ撮像シーケンスは、互いに共に且つ印加磁場勾配と共に用いられて選択された核磁気共鳴の操作を達成する短期間のＲＦパケットシーケンス群のパケットからなる。ＲＦパルスは、共鳴を飽和させ、共鳴を励起し、磁化を反転し、共鳴をリフォーカスし、或いは共鳴を操作し、そして検査ボリューム内に位置する人体１０の一部を選択するために使用される。ＭＲ信号もＲＦアンテナ９によって捕捉され得る。
【００３８】
例えばパラレルイメージング法による、人体又は一般的には対象１０の限られた領域のＭＲ画像の生成のため、撮像のために選択された領域に近接して一組の局部アレイＲＦコイル１１、１２及び１３が配置される。アレイコイル１１、１２及び１３は、ＲＦアンテナによって発せられるＲＦ送信によって誘起されたＭＲ信号を受信するために使用されることができる。しかしながら、検査ボリュームにＲＦ信号を送信することにアレイコイル１１、１２及び１３を使用することも可能である。
【００３９】
結果として得られるＭＲ信号は、ＲＦアンテナ９によって、且つ／或いはＲＦコイル１１、１２及び１３のアレイによって捕捉され、好ましくは前置増幅器（図示せず）を含む受信器１４によって復調される。受信器１４は、送／受スイッチ８を介してＲＦコイル９、１１、１２及び１３に接続される。
【００４０】
例えばエコープラナーイメージング（ＥＰＩ）、エコーボリュームイメージング、グラディエントスピンエコーイメージング、ファーストスピンエコーイメージングなどの複数の撮像シーケンスのうちの何れかを生成するように、ホストコンピュータ１５が傾斜パルス増幅器３及び送信器７を制御する。選択されたシーケンスで、受信器１４は、各ＲＦ励起パルスに素早く続いて単一又は複数のＭＲデータラインを受信する。データ収集システム１６が、受信信号のアナログ−デジタル変換を実行し、各ＭＲデータラインを更なる処理に適したデジタルフォーマットに変換する。最近のＭＲ装置においては、データ収集システム１６は、ロー（生）画像データの収集に特化された別個のコンピュータである。
【００４１】
最終的に、デジタルロー画像データは、フーリエ変換又は例えばＤｉｘｏｎ再構成のような適切な再構成アルゴリズムを適用する再構成プロセッサ１７によって、画像表現へと再構成される。ＭＲ画像は、患者を通る１つの平面スライス、平行な複数の平面スライスの配列、３次元ボリューム、又はこれらに類するものを表現し得る。そして、画像は画像メモリに格納され、そこで画像は、画像表現のスライス又はその他の部分を視覚化（例えば、結果として得られたＭＲ画像のヒト読み取り可能表示を提供するビデオモニタ１８による）に適したフォーマットに変換するためにアクセスされ得る。
【００４２】
再構成プロセッサ１７及びホストコンピュータ１５上で本発明に係る方法を用いることにより、水信号、脂肪信号、又は水信号と脂肪信号との双方の識別のために、個々のスペクトルピークの部分的なディフェージング及びリフェージングによる脂肪信号における振幅変調を利用することが可能である。脂肪のマルチピークスペクトルモデルを考えるに、脂肪信号の振幅の差を最大することを狙って、二点Ｄｉｘｏｎ撮像において２つのエコー時間を選定することが可能である。故に、２つのエコー時間にて信号方程式を解くことによって得られる一般的な２つの選択的な解は、広範囲の信号振幅比において、単一の一義的な解へと減少される。
【００４３】
以下、上述の方法を用いた例示的な画像再構成プロセスを更に詳細に説明する。
【００４４】
画像空間内の複合的な複素信号ｓが、
【００４５】
【数１】
0005914488.tif 000002
によってモデル化される。ただし、Ｗ及びＦは画像空間内の水信号及び脂肪信号を表し、ψ及びｅｉψは位相誤差及び対応するフェーザを表す。複素因子ｃが、
【００４６】
【数２】
0005914488.tif 000003
によって与えられる。ただし、ｗは合計で１になる重みを表し、θｍは、脂肪スペクトルのｍ番目のピークの共鳴周波数の水に対するオフセットをΔｆｍとして、２πΔｆｍＴＥに等しく、ＴＥはエコー時間である。
【００４７】
二点Ｄｉｘｏｎ撮像においては、Ｓが２つのエコー時間にて測定される。そして、これらのサンプルＳ１及びＳ２に基づいて、水信号と脂肪信号との分離及び識別が試みられる。因子ｃの時間変動、故に、純粋な脂肪信号の時間変動を図２に例示する。図２では、１．５Ｔ（ＦＩＧ．２）及び３．０Ｔ（ＦＩＧ．２ｃｏｎｔｉｎｕｅｄ）で、シングルピーク及び７ピークの脂肪スペクトルモデルを仮定し、純粋な脂肪信号の振幅及び位相をエコー時間（ＴＥ）の関数として概略的に示している。
【００４８】
次式：
【００４９】
【数３】
0005914488.tif 000004
を考えるに、Ｓ１とＳ２との間の振幅の差が主として、先ほど二点Ｄｉｘｏｎ撮像において仮定したように水信号と脂肪信号とのディフェージングに由来するのでなく、脂肪信号自体に由来し得ることは明らかである。これは特に、因子ｃの大きさが急速に減少する非常に短いエコー時間の場合に当てはまる。
【００５０】
故に、エコー時間の好適な選定により、この振幅変調を利用して、水信号と脂肪信号との識別における生来的な不明確性を解消することが可能である。
【００５１】
２つの非線形方程式の系を解くことによって得られる２つの可能性ある解：
【００５２】
【数４】
0005914488.tif 000005
を図３にて特徴付ける。第１のエコー時間での複合信号の振幅を｜Ｓ１｜＝１に固定したときの、第２のエコー時間での複合信号の振幅｜Ｓ２｜の関数としての水信号Ｗ及び脂肪信号Ｆの解が、ＴＥ１／ＴＥ２＝１．５／２．８ｍｓ（ＦＩＧ．３）及びＴＥ１／ＴＥ２＝１．８／３．２ｍｓ（ＦＩＧ．３ｃｏｎｔｉｎｕｅｄ）について示されている。下のグラフは上のグラフを拡大したものである。
【００５３】
｜Ｓ２｜の値（又は、一般的には比｜Ｓ１｜／｜Ｓ２｜）が非常に高い、あるいは低い場合、取り得る解は複素数である。この比が中間的な値の場合、取り得る解は実数であるが、負の場合もある。これら全ての特性を考慮して、｜Ｓ１｜／｜Ｓ２｜の値の範囲を図４にて分類する。
【００５４】
図４は、水信号及び脂肪信号の妥当解すなわち負でない（非負の）実数解の個数に従って、２つのエコー時間における信号振幅の比を３つのカテゴリーに分類したものを示している。故に、“０”は、その範囲内の信号振幅の比では、取り得る何れの解も負又は複素数の何れかであるので、水信号及び脂肪信号の妥当な解が存在しないことを指し示している。“１”は、水信号及び脂肪信号の取り得る解が１つ存在することを指し示し、“２”は、水信号及び脂肪信号の取り得る解が２つ存在することを指し示している。故に、領域“２”では不明確性（両義性）が残っているが、“１”として分類された小範囲では水信号と脂肪信号との識別が直ちに達成される。
【００５５】
故に、図４により、収集された画像データセット内の信号強度から、モデル化が水／脂肪信号分離の単一の一義的な数学解を生じさせるボクセルを特定することが可能である。
【００５６】
上述のように、信号分離においてモデル化が２つ以上の数学解を生じさせるボクセルの数を最小化するために、使用するエコー時間を前もって決定することが好ましい。これを、水信号及び脂肪信号の場合について、以下にて単純に例証する。
【００５７】
上述の２つの非線形方程式の系には２つの自明な解が存在する。第１の解は純粋な水信号：
｜Ｓ１｜＝Ｗ
｜Ｓ２｜＝Ｗ
であり、
ｒ＝｜Ｓ１｜／｜Ｓ２｜
で定義される信号振幅比ｒは
ｒ＝１
である。
【００５８】
第２の解は純粋な脂肪信号：
｜Ｓ１｜＝｜ｃ１｜Ｆ
｜Ｓ２｜＝｜ｃ２｜Ｆ
であり、
ｒ＝｜ｃ１｜／｜ｃ２｜
である。
【００５９】
その他の解は、例えば、Ｗに関して第１の方程式を、
【００６０】
【数５】
0005914488.tif 000006
と解き、この式を第２の方程式に代入して、
ａ１Ｆ４＋ａ２Ｆ２＋ａ３＝０
ただし、
【００６１】
【数６】
0005914488.tif 000007
を生成することによって得られる。明らかなように、ａ１は信号に無依存であり、且つ
ａ１≧０
ａ３≧０
である。ａ１≠０の場合、Ｆの２つの解が、
【００６２】
【数７】
0005914488.tif 000008
によって与えられ、２つの更なる解が、
【００６３】
【数８】
0005914488.tif 000009
によって与えられる。
【００６４】
後者は、負又は複素数の何れかであるので検討しない。
【００６５】
前者が実数である場合、それらは明らかに非負である。従って、２つの解の範囲を実数であるように構築すれば十分である。２つの必要条件は、
【００６６】
【数９】
0005914488.tif 000010
であり、ここで、最初の式は自明であり、第２の式は、次式の左辺が実数のままであるとして、
【００６７】
【数１０】
0005914488.tif 000011
に従っている。これらの条件下で、Ｆの２つの解はｒの同じ範囲で実数である。これにより、
ａ２２≧４ａ１ａ３
ａ２≦０
がもたらされ、ひいては、
ａ２≦−２√（ａ１ａ３）
がもたらされる。
【００６８】
単純化のため、幾つかの定数を、
【００６９】
【数１１】
0005914488.tif 000012
と定義する。Ｓ２≠０と仮定すると、
【００７０】
【数１２】
0005914488.tif 000013
となり、ｒ＞１の場合、
ｄ４ｒ２≦ｄ３
となる。故に、ｄ４＞０の場合、
ｒ≦√（ｄ３／ｄ４）
なる上限が得られる。ｄ３≧ｄ４であるので、これは１以上である。ｄ４≦０の場合には、
√ａ１≧ｄ１
であり、故にｄ３≧０であるので、制限は存在しない。ｒ＝１の場合にも制限は存在しない。ｒ＜１の場合には、
ｄ６ｒ２≦ｄ５
であり、ただし、ｄ６≦０である。ｄ５≧ｄ６であるので、ｄ６＝０の場合、制限は存在せず、ｄ６＜０の場合、
ｒ≧√（ｄ５／ｄ６）
の下限が得られる。これは最大で１である。
【００７１】
確立されたＦの実数の非負の解の限界内で、元々の２つの方程式のうちの一方への代入により、対応するＷの解が得られる。代替的に、これら２つの元々の方程式の間の差が用いられてもよく、それにより、
【００７２】
【数１３】
0005914488.tif 000014
がもたらされる。
【００７３】
この手法は、Ｆの各解に対してＷの一方の解のみが得られるという利点を有する。しかしながら、ｃ１Ｒ≠ｃ２Ｒ且つＦ＞０の場合にのみ妥当である。Ｆ＝０の場合は先ほどのものでカバーされているが、ｃ１Ｒ＝ｃ２Ｒの場合はなおも分析される必要がある。その場合、Ｆの該当する解は、
【００７４】
【数１４】
0005914488.tif 000015
であり、元々の第１の方程式への代入により、
【００７５】
【数１５】
0005914488.tif 000016
がもたらされる。ただし、ｃＲ＝ｃ１Ｒ且つｃＲ＝ｃ２Ｒである。以下ではＷ＋及びＷ−と表記するＷの両方の解も元々の第２の方程式を満足する。これらの解が実数である条件は、Ｆに関し、｜ｃ１Ｉ｜＞｜ｃ２Ｉ｜でｒ≧１、｜ｃ１Ｉ｜＜｜ｃ２Ｉ｜でｒ≦１であり、且つ、
｜ｃ１Ｉ｜＞｜ｃ２Ｉ｜で
ｒ≦｜ｃ１Ｉ｜／｜ｃ２Ｉ｜
｜ｃ１Ｉ｜＜｜ｃ２Ｉ｜で
ｒ≧｜ｃ１Ｉ｜／｜ｃ２Ｉ｜
である。Ｗが非負である条件は、
【００７６】
【数１６】
0005914488.tif 000017
であり、これは、ｃＲ≧０、｜ｃ１Ｉ｜＞｜ｃ２Ｉ｜、Ｓ２≠０且つ＋符号に関して、
ｒ≦｜ｃ１｜／｜ｃ２｜
をもたらす。
【００７７】
同様に、その他の場合を生じさせる可能性ある限界を得ることができる。全体として、ｃＲ≧０の場合、Ｗ＋が、｜ｃ１｜＞｜ｃ２｜の場合の１≦ｒ≦｜ｃ１｜／｜ｃ２｜、及び、｜ｃ１｜＜｜ｃ２｜の場合の｜ｃ１｜／｜ｃ２｜≦ｒ≦１の妥当な解であり、Ｗ−は妥当でない。ｃＲ＜０の場合には、Ｗ＋が、｜ｃ１｜＞｜ｃ２｜の場合の１≦ｒ≦｜ｃ１Ｉ｜／｜ｃ２Ｉ｜、及び、｜ｃ１｜＜｜ｃ２｜の場合の｜ｃ１Ｉ｜／｜ｃ２Ｉ｜≦ｒ≦１の妥当な解であり、Ｗ−が、｜ｃ１｜＞｜ｃ２｜の場合の｜ｃ１｜／｜ｃ２｜≦ｒ≦｜ｃ１Ｉ｜／｜ｃ２Ｉ｜、及び、｜ｃ１｜＜｜ｃ２｜の場合の｜ｃ１Ｉ｜／｜ｃ２Ｉ｜≦ｒ≦｜ｃ１｜／｜ｃ２｜の妥当な解である。故に、後者の範囲においては、２つの解が存在する。
【００７８】
ここで、ｃ１Ｒ≠ｃ２Ｒ且つＦ＞０を仮定すると、Ｗに関して解くことは、
【００７９】
【数１７】
0005914488.tif 000018
をもたらす。Ｆが実数であり、Ｗも実数である場合、これは、
【００８０】
【数１８】
0005914488.tif 000019
に等価である。
【００８１】
明らかなように、このような解はｒ＜１に関して存在しない。同様に、ｃ１Ｒ＜ｃ２Ｒ且つ｜ｃ１｜＜｜ｃ２｜の場合、このような解はｒ＞１に関して存在しない。対照的に、ｃ１Ｒ＜ｃ２Ｒ且つ｜ｃ１｜＞｜ｃ２｜の場合のｒ＜１に関して、及びｃ１Ｒ＞ｃ２Ｒ且つ｜ｃ１｜＜｜ｃ２｜の場合のｒ＞１に関しては、このような解が２つ存在する。Ｆの解を代入することにより、
【００８２】
【数１９】
0005914488.tif 000020
が得られる。これは、
【００８３】
【数２０】
0005914488.tif 000021
を用いると、
【００８４】
【数２１】
0005914488.tif 000022
と書き換えられることができ、｜Ｓ２｜≠０の場合、
【００８５】
【数２２】
0005914488.tif 000023
と書き換えられ得る。暫定的に、この不等式を等式と見なす。取り得る解は、二乗をとって、
【００８６】
【数２３】
0005914488.tif 000024
を得ることによって見出される。自明な解はｒ＝±１である。ａ１（ｒ２−１）で割ることにより、
【００８７】
【数２４】
0005914488.tif 000025
及び
【００８８】
【数２５】
0005914488.tif 000026
が得られる。これは、
【００８９】
【数２６】
0005914488.tif 000027
ひいては、ｒ＝±｜ｃ１｜／｜ｃ２｜に単純化されることができる。
【００９０】
従って、２つの正の解のみが関心あるものであるＦ＋及びＦ−の双方に対して、ｒの取り得る４つの解が存在する。残るは、それらのうちの何れが現実の解であるか、及びどの範囲内で元の不等式が成り立つかを決定することである。
【００９１】
ａ２１＋ａ２２＜０であるので、ｒ＝±１がＦ−のみの現実解であることを示すのは容易である。ｒのその他の取り得る解を元の不等式に代入すると、
【００９２】
【数２７】
0005914488.tif 000028
がもたらされる。この不等式を用いることで、ついに、以下の４つの異なるケース（場合）の各ｒの妥当解の個数を決定することが可能である。
【００９３】
ケース１：｜ｃ１｜＞｜ｃ２｜且つｃ１Ｒ＞ｃ２Ｒ
ａ）  ｒ＜１  ：解なし
ｂ）  １≦ｒ≦｜ｃ１／ｃ２｜  ：Ｆ−，Ｗ−（すなわち、単一の解）
ｃ）  ｒ＞｜ｃ１／ｃ２｜  ：ｃ１Ｒ｜ｃ２｜２−ｃ２Ｒ｜ｃ１｜２＞０の場合は、Ｆ−，Ｗ−、及びＦ＋，Ｗ＋（すなわち、２つの解）、その他の場合は解なし。
【００９４】
従って、ケース１では、如何なる測定も水／脂肪信号分離に関して単一の一義的な解をもたらすよう、条件ｃ１Ｒ｜ｃ２｜２−ｃ２Ｒ｜ｃ１｜２＜０で動作することが望ましい。それ以外では、なおも通常、水が支配的な信号に対して単一の一義的な解が得られる。
【００９５】
ケース２：｜ｃ１｜＞｜ｃ２｜且つｃ１Ｒ＜ｃ２Ｒ
ａ）  ｒ＜１  ：Ｆ−，Ｗ−、及びＦ＋，Ｗ＋（すなわち、２つの解）
ｂ）  １≦ｒ≦｜ｃ１／ｃ２｜  ：Ｆ＋，Ｗ＋（すなわち、単一の解）
ｃ）  ｒ＞｜ｃ１／ｃ２｜  ：ｃ１Ｒ｜ｃ２｜２−ｃ２Ｒ｜ｃ１｜２＞０の場合は、Ｆ−，Ｗ−、及びＦ＋，Ｗ＋（すなわち、２つの解）、その他の場合は解なし。
【００９６】
従って、ケース２でもやはり、純粋あるいは略純粋な脂肪信号に対して通常単一の一義的な解が得られるよう、条件ｃ１Ｒ｜ｃ２｜２−ｃ２Ｒ｜ｃ１｜２＜０で動作することが望ましい。しかしながら、ケース２は一般に、部分的な識別が可能であるのみであり、ケース１より望ましくない。
【００９７】
ケース３：｜ｃ１｜＜｜ｃ２｜且つｃ１Ｒ＞ｃ２Ｒ
ａ）  ｒ＜｜ｃ１／ｃ２｜  ：ｃ１Ｒ｜ｃ２｜２−ｃ２Ｒ｜ｃ１｜２＜０の場合は、Ｆ−，Ｗ−、及びＦ＋，Ｗ＋（すなわち、２つの解）、その他の場合は解なし
ｂ）  ｜ｃ１／ｃ２｜≦ｒ≦１  ：Ｆ＋，Ｗ＋（すなわち、単一の解）
ｃ）  ｒ＞１  ：Ｆ−，Ｗ−、及びＦ＋，Ｗ＋（すなわち、２つの解）。
【００９８】
従って、ケース３では、純粋あるいは略純粋な脂肪信号に対して通常単一の一義的な解が得られるよう、条件ｃ１Ｒ｜ｃ２｜２−ｃ２Ｒ｜ｃ１｜２＞０で動作することが望ましい。しかしながら、ケース３は一般に、部分的な識別が可能であるのみであり、ケース４より望ましくない。
【００９９】
ケース４：｜ｃ１｜＜｜ｃ２｜且つｃ１Ｒ＜ｃ２Ｒ
ａ）  ｒ＜｜ｃ１／ｃ２｜  ：ｃ１Ｒ｜ｃ２｜２−ｃ２Ｒ｜ｃ１｜２＜０の場合は、Ｆ−，Ｗ−、及びＦ＋，Ｗ＋（すなわち、２つの解）、その他の場合は解なし
ｂ）  ｜ｃ１／ｃ２｜≦ｒ≦１  ：Ｆ＋，Ｗ＋（すなわち、単一の解）
ｃ）  ｒ＞１  ：解なし。
【０１００】
従って、ケース４では、如何なる測定も水／脂肪信号分離に関して単一の一義的な解をもたらすよう、条件ｃ１Ｒ｜ｃ２｜２−ｃ２Ｒ｜ｃ１｜２＞０で動作することが望ましい。それ以外では、なおも通常、水が支配的な信号に対して単一の一義的な解が得られる。
【０１０１】
結論として、ケース１及び４は原理上ｒの全範囲で水−脂肪の識別を可能にするので、ケース１及び４が好ましい。驚くべきことに、ヒトの脂肪信号をモデル化するために一般的に適用されるスペクトルモデルにおいて、この条件は典型的に、双方のエコー時間において水及び脂肪からの信号が同相に近くなる２つのエコー時間が選択されるときに満足される。
【図１】
0005914488.tif 000029
【図２】
0005914488.tif 000030
【図３−１】
0005914488.tif 000031
【図３−２】
0005914488.tif 000032
【図４】
0005914488.tif 000033