JPH1040365A

JPH1040365A - 磁気共鳴映像法および装置

Info

Publication number: JPH1040365A
Application number: JP9089747A
Authority: JP
Inventors: Haiying Liu; リウハイイン; James B Mcnally; ビーマックナリージェームズ; Gordon D Demeester; ディーディメースターゴードン
Original assignee: Picker International Inc
Current assignee: Philips Nuclear Medicine Inc
Priority date: 1996-04-26
Filing date: 1997-04-08
Publication date: 1998-02-13
Also published as: EP0803740A1; US5742163A

Abstract

(57)【要約】【課題】磁気共鳴映像における、二次元位相補正【解決手段】較正データ列は各データ収集期間中に、ゼ
ロ補充されて（８６）、完全なデータの組を生成し、そ
れぞれデータ収集期間の一つに対応する一連の低分解能
複素画像（９０₁、９０₂、・・・９０_n) にフーリエ
変換される（８８）。低分解能画像は正規化され（９
２）、それらの複素共役が取られる（９４）。映像デー
タ列はデータ収集期間により振分けられ、完全なデータ
の組を生成するためにゼロ補充される（１０４）。各デ
ータ・サンプリング期間に対応する完全なデータの組
が、部分画像表現（１０６₁、１０６₂、１０６_n）へ
とフーリエ変換される。各部分画像が正規化された位相
補正マップ（９６）の複素共役と乗算され（１０８）、
補正された部分画像を生成し、それを総計して（１１
２）、合成画像生成する（１１４）。合成画像は密度補
正される（１２０）。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は磁気共鳴映像方法お
よび装置に関する。本発明は特にエコー平面映像法（Ｅ
ＰＩ）に関連し、それにつき更に詳しく説明する。しか
しながら、本発明の技術は、例えば、エコー容積映像
（ＥＶＩ）法、急速スピンエコー（ＦＳＥ）法、ならび
に勾配スピンエコー（ＧＳＥ）法等の反復勾配エコー、
スピンエコーあるいは勾配およびスピンエコーの組合せ
からなる、他の急速連続映像法にも利用出来ることが理
解されるべきである。

【０００２】

【従来の技術】従来、磁気共鳴映像の被写体は、選択さ
れた双極子が磁場に好ましい方向に整列されるよう、一
時的定常磁場内に置かれていた。無線周波数パルスが印
加されて、好ましい方向に整列された双極子が共鳴し、
ある特性を有する共鳴無線周波数の磁気共鳴信号を放射
するようにされた。共鳴中の双極子からの無線周波数共
鳴信号が読み出されて、画像として再構成され、表現さ
れる。二次元フーリエ変換映像法では、読取り勾配が、
エコーの読み取りの間に、読み取り軸に沿って周波数エ
ンコーディングのために印加され、かつ位相エンコード
勾配が、エコーの間に位相エンコード軸に沿って、段階
的位相エンコーディングのためにパルスとして印加され
る。このようにして、各エコーはｋ空間中にデータ列を
生成させる。データ列の相対的位相エンコーディング
が、ｋ空間中のそれらの相対的位置を制御している。従
来例では、ゼロ位相エンコーディングのデータ列は、ｋ
空間中心を横切って延びている。段階に増加する位相エ
ンコーディング勾配を有するデータ列は、一般的にｋ空
間の中心線の上に描かれる。そして段階的に減少する位
相エンコーディング勾配でのデータ列は、ｋ空間の中心
線の下に描かれる。このようにして、２５６×２５６、
または５１２×５１２でのデータ値マトリックスのよう
なマトリックスがｋ空間中に生成される。これらの値の
フーリエ変換で、従来の磁気共鳴映像が形成される。

【０００３】受信した磁気共鳴信号を強化するために、
普通、最初の信号は再びあるエコーに集束される。これ
は磁場勾配の極性を反転させて、磁場あるいは勾配・エ
コーを誘導させることにより行なわれる。同様に、無線
周波数励振パルスの後に、１８０°パルスを加えて、信
号をスピン・エコーとして再び集束させてもよい。更
に、磁場勾配の極性の反転を反復して、一連の勾配エコ
ーを、各無線周波数励振パルスの後に生成させることも
出来る。また他の選択肢として、単一の無線周波数励振
パルスの後にスピンおよび勾配エコーを混合して加える
ことも出来る。例えば、ホーランド等の米国特許第4,８
３3,４０８号を参照されたい。２０ｍＴ／ｍ／ｍｓのオ
ーダのスルー・レートと、１５ｍＴ／ｍの勾配強度を有
する勾配システムを使用するシングルショットエコー平
面映像（ＥＰＩ）シークエンスでは、任意の磁化の傾斜
角を生じる単一無線周波数励振パルスに、十分な回数の
勾配反転を続けて、一組のデータ列全部を発生させるこ
とが出来る。被写体からの磁気共鳴データは、磁場勾配
の方向に画像被写体をエンコーディングする振動読取り
勾配を用いて一連のエコーの間に収集される。例えば、
Ｐ．マンスフィールド（Mansfield)著、J.Phys. Chemis
try Vol. 10. Ｌ５５−Ｌ５８頁、（１９７７年）参
照。更に、ｋ空間中のデータ列を段階的に得るために、
各エコーの前に、読取り勾配方向に垂直な一連の位相エ
ンコーディング勾配パルスが印加される。被写体の画像
がエコー・データの二つの一次元フーリエ変換により得
られることが好ましい。このシングル・ショットＥＰＩ
法は、秒以下の時間単位での真の動画映像法用の超高速
画像技術を提供する。

【０００４】

【発明が解決しようとする課題】しかしながら、シング
ル・ショットＥＰＩ法は、（特にシークエンスの最後に
近いエコーについて）比較的低いＳＮ比、低空間的分解
能、化学シフト・アーティファクト、位相エンコーディ
ング方向に沿った形状歪み等を含む欠点がある。内在的
に位相エンコーディング方向に沿ったデータ・サンプリ
ング帯域幅が狭いので、シングル・ショットＥＰＩ法
は、化学シフト・アーティファクトと同様、磁場の不均
一性により引起される形状歪みにさらされる傾向がある
が、これは重大な場合もありうる。マルチ・ショットＥ
ＰＩ法は、上記のシングル・ショットＥＰＩ法の欠点を
減少させる。従来のマルチ・ショットＥＰＩ法では、位
相−周波数空間、すなわちｋ空間は、例えば４〜１６区
分の複数の区分に分られている。共鳴励振の後、読取り
勾配が発振され、ｋ空間中の各区分で一連のエコーを、
それゆえデータ列を生成する。他の励振により各区分中
の複数のデータ列の内の異る列が生成される。ｋ空間が
この交互に繰返される方法で満たされるまで、この過程
が繰返される。

【０００５】データ収集の過程で、スライディング・サ
ンプリング・ウインドウ技術を利用して、データ中の位
相と強さの不連続性を平滑化し、あるいは整合を助ける
ことも可能である。ブッツ（Butts)他、著「標準的ＭＲ
Ｉシステムにおける交互的エコー平面映像法」、Mag. R
es. Med. Vol. 31、６７−７２頁（１９９４）、または
マキノン著、「標準的スキャナーにおける超高速交互的
勾配エコー平面映像法」Mag. Res. Med. Vol. 30、６０
９−６１９頁、（１９９３年）参照。スライディング・
ウインドウ技術は、ｋ空間の連続した区分の各対の間の
位相と振幅の転移の平滑化を助ける。これは最終的な復
元された画像のリンギング・アーティファクトを減少さ
せる。ｋ−空間の区分の位相不連続性から生ずるアーテ
ィファクトを減少させる他の先行技術は空間的な一次元
位相補正に基いている。更に詳しくは、データ列につい
ての一次位相補正（first order phase correction) は
異なるエコーからなっている。マイヤー他の米国特許第
5,１５1,６５６号参照。

【０００６】

【課題を解決するための手段】本発明によれば、磁気共
鳴映像のためのシステムと方法が提供されている。本発
明によれば、このシステムは検査領域を通る時間的に一
定な磁場を生成する磁石と、検査領域にある双極子に共
鳴を誘起させ、無線周波数共鳴信号を発生させるための
無線周波数パルス制御器と送信器とを有している。また
このシステムは勾配磁場コイルと勾配磁場制御器とを有
し、少なくとも検査領域を直角方向に横断する位相およ
び読取り磁場勾配パルスを生成させ、かつ読取り勾配を
反復して反転させる。更に、このシステムは受信器を有
し、読み取り勾配が反転する毎にその後、無線周波数信
号を受信して復調し、一連のデータ列を作りだし、かつ
復元された画像表現を記憶する画像メモリーを有してい
る。更にこのシステムの特徴は複数の二次元位相補正マ
トリックスを生成させる二次元位相補正プロセッサーを
有することである。各二次元位相補正マトリックスは、
読取り勾配の反転の一つに対応する。二次元位相補正プ
ロセッサーは、読取り勾配の各反転後のデータ収集期間
に収集されたデータ列を受信して、各データ収集期間に
ついて、対応する二次元位相補正マトリックスを独立に
生成させる。一個の部分画像プロセッサーが各データ収
集期間に生成された画像データを受信する。この部分画
像プロセッサーは、各データ収集期間に収集された一組
の不完全なデータ列を、独立に対応する部分画像に再構
成する。一個のマトリックス乗算器が、各部分画像を同
じデータ収集期間に対応する二次元位相補正マトリック
スと乗算し、複数の二次元的に位相補正された部分画像
を生成する。画像加算器が二次元的に位相補正された部
分画像を加算して、合成画像を生成し、画像メモリーに
記憶する。

【０００７】本発明の他の側面によれば、エコーを形成
するよう誘起される双極子に、磁気共鳴が励起される磁
気共鳴映像法が提供されている。エコーは読取り勾配が
存在する時に、読取り軸に沿って読み取られ、一連のデ
ータ列を形成し、出力画像として再構成される。この方
法の特徴はエコーを位相エンコーディングして、ｋ空間
の複数の区分のそれぞれに沿っているデータ列を生成す
ることである。データ列はｋ空間の複数の区分の中に振
分けられる。位相補正マトリックスは、データ列から生
成される。各マトリックスはｋ空間の区分の一つに対応
する。複数の部分画像がデータ列から生成される。各部
分画像はｋ空間の区分の一つに対応する。複数の位相補
正された部分画像が、複数の部分画像と位相補正マトリ
ックスから生成される。位相補正された部分画像のそれ
ぞれは、ｋ空間の区分の一つに対応する。位相補正され
た部分画像が、総計されて合成画像を生成する。

【０００８】本発明の利点の一は、二次元位相補正の提
供にある。本発明の他の利点は、それにより多数ショッ
トＥＰＩ画像データが容易にえられることである。本発
明の他の利点は、機器の性能向上が不要なことである。
本発明の他の利点は画質の改善で、特に空間的分解能の
改善とギブス（Gibbs)リンギングと歪みとの低減であ
る。本発明の他の利点は、この技術による走査時間の大
きな増加がないことである。本発明の他の利点は、この
技術により更に安定した信頼性のある位相補正が得られ
ることである。

【０００９】本発明の更なる他の利点は、較正のために
一時的に保存されるデータの量を減少させることが出来
ることである。

【００１０】

【発明の実施の形態】本発明を添付図面を参照して、実
施例にもとづき更に説明する。図１についてみると、主
磁場制御１０が、超伝導または抵抗型磁石１２を制御し
て、ほぼ均一な時間的に一定な磁場を検査領域１４を通
るｚ軸に沿って形成する。磁気共鳴エコー手段が、磁気
共鳴映像生成と分光学的シーケンスを形成するために、
一連の高周波（ＲＦ）と磁場勾配パルスを印加し、磁気
スピンを反転、または励起し、磁気共鳴を誘導し、磁気
共鳴を再集束し、磁気共鳴を走査し、空間的およびその
他で磁気共鳴をエンコードし、スピンを飽和させる等を
行なう。更に詳しくは、勾配パルス増幅器２０が、全身
勾配コイル２２の選ばれたもの、あるいは対に電流パル
スを印加し、検査領域１４のｘ、ｙ、ｚ軸に沿って磁場
勾配を生みだす。デジタル無線周波数送信器２４が無線
周波数パルス、あるいはパルス・パケットを全身ＲＦコ
イル２６へと送信し、ＲＦパルスを検査領域へと送信す
る。典型的な無線周波数パルスは直接連続する短かい継
続期間のパルス区分のパケットからなり、それらは互に
かつ印加された勾配と相俟って、選択された磁気共鳴操
作を達成する。ＲＦパルスはスピンを飽和させ、共鳴を
励起し、磁化を反転させ、共鳴を再集束し、検査領域の
選択された部分での共鳴を操作する。全身撮影のために
は、共鳴信号は通常全身ＲＦコイル２６によりピックア
ップされる。

【００１１】被写体の限定された領域の画像の生成のた
めには、局所コイルが通常選択された領域の近傍に置か
れる。例えば挿入可能な頭部コイル３０が、挿入孔と同
心点の選択された頭脳領域を囲んで挿入される。挿入可
能な頭部コイル３０は、好ましくは局所勾配コイル３２
を有し、局所勾配コイル３２は勾配増幅器２０から電流
パルスを受けて、頭部コイル３０内の検査領域のｘ、
ｙ、ｚ軸に沿って磁場勾配を生みだす。局所無線周波数
コイル３４が、磁気共鳴を励振し、患者の頭部から発す
る磁気共鳴信号を受信するために使用される。それに代
り、受信専用局所無線周波数コイルが、身体−コイルＲ
Ｆ送信により誘導される共鳴信号を受信するために使用
出来る。一枚のＲＦスクリーン３６により、ＲＦ頭部コ
イルからのＲＦ信号が、勾配コイルと周囲の構造物中に
うず電流を誘起させることを阻止している。その結果生
じる無線周波数信号は全身ＲＦコイル２６、局所無線周
波数コイル３４、あるいは他の専用ＲＦコイルによって
ピックアップされ、受信器３８により復調される。

【００１２】シーケンス制御回路４０が、勾配パルス増
幅器２０と、送信器２４とを制御して、エコー平面映
像、エコー容積映像、勾配およびスピンエコー映像、急
速スピンエコー映像等を含む、複数の多エコーシーケン
スのどれかを発生させる。選択されたシーケンスについ
て、受信器３８は各ＲＦ励振パルスに続いて、急速に連
続して複数のデータ列を受信する。受信器３８はデジタ
ル受信器であるか、あるいはここに示されるように、各
データ列をデジタル形式に変換するため、Ａ−Ｄコンバ
ータ４２を伴なっていることが好ましい。図１および２
において、デジタルデータ列は先ず、較正データ列と画
像データ列とに区別される４４。較正データ列とｋ空間
中でエコー番号により区分１から区分ｎにまで４６で区
別される。画像データ・ラインは同様に４８でｋ空間の
ｎ個の区分中に分けられる。

【００１３】図３においてエコー平面映像シークエンス
は、任意に条件設定ルーチン５０、例えばスライス選択
勾配パルスを伴なう１８０°反転パルスで開始出来る。
適当な遅延を加えることにより、例えば反転回復画像が
生成される。一つの無線周波数励振パルス５２がスライ
ス選択勾配５４と同時に印加され、選択されたスライス
内で、選択された磁化傾斜角、例えばエルンスト（Erns
t)角で磁化を傾斜させる。位相エンコード勾配５６が、
ｋ空間のｎ個の区分の最初の一つに与えられるる。一個
の負極性の読取り勾配部分５８が読み取り、あるいは周
波数エンコード方向に沿って位相を乱れ（dephase)させ
る。その後、逆極性の読取り勾配部分６０が、磁化の位
相を勾配エコー６２の位相に再び合わせる（rephase)
が、勾配エコー６２は、位相エンコード勾配５６によ
り、設定された位相で位相エンコーディングされる。そ
の後、読取り勾配は複数回反転され（図示された実施例
では４回）、対応する数の勾配エコーを発生させる。位
相エンコード勾配パルス６４が、各勾配エコーの間に印
加され、次に、ｋ空間のｎ個の領域の何れかの後続の勾
配エコーからのデータ列を得る。それにより、読取り勾
配はｎ回反転されて、ｎ個の勾配エコーを発生させる
が、これはｋ空間のｎ個の区分のそれぞれについて一個
となる。図３はｎ＝４の実施例を図示している。図３に
示されたＥＰＩシーケンスのスピン・エコーもまた好ま
しく、そこでは再集束ＲＦパルスは、励振パルス５２と
読取り勾配６０の間に挟まれる。

【００１４】最後の勾配エコーが読まれた後、位相乱し
ルーチン（dephasing routine)６６が印加される。位相
乱しルーチンには、位相アンラップ（unwrap) 勾配６８
が含まれてもよく、位相アンラップ勾配６８は、位相エ
ンコード勾配５６と全ての位相エンコード勾配パルス６
４の合計の負数に等しい。位相乱し法はまた、残存位相
エンコーディングを乱すための、スポイラ勾配７０を含
んでもよい。これにより次の繰返しがほぼ直に開始する
ことが可能となる。これにかわり、マルチ・スライス法
では、スライス内での位相エンコーディングが時間の経
過中と共に位相を乱して、他の画像スライスの一つから
ｎデータ列を集めるために、シーケンスを直に繰り返し
ても良い。較正期間中、位相エンコーディング標点（bl
ips)６４は、各エコーが同じ位相エンコーディングを有
するようにゼロに設定される。各繰り返しにおいて、最
初の位相エンコーディング勾配５６が、ゼロ位相エンコ
ード原点の近くの少数の位相エンコード値のそれぞれに
与えられる。較正シーケンスは、その都度異なる位相エ
ンコーディングで約２〜６４回繰り返すことが好まし
く、その位相エンコーディングは全てのｎエコーに共通
である。

【００１５】更に図３および図２において、各エコーに
対応する較正データ列は、４６で一連の対応する較正デ
ータ列メモリ８０₁、８０₂、・・・、８０_nに振分け
られる。再度、好ましい実施例では、ｎ個のエコー位置
のそれぞれについて、約２〜６４のデータ列が発生され
る。好ましくは、ブラックマン（blackman）、ハニング
（Hanning)、ハミング（Hamming)、あるいはガウス（Ga
ussian) フィルタであるフィルタ８２が、各データ列を
読み取り方向に濾波する。濾波されたデータ列は、対応
する較正データ列メモリ８４₁、８４₂、・・・、８４
_nにロードされる。ゼロ−補充ルーチン（zero-filling
routine) ８６が、濾波後の較正データ列メモリ８
４₁、８４₂、・・・８４_nの中の空データ列を満たし
て、メモリ８４₁、８４₂、・・・８４_nの中に、完全
な矩形マトリックスを形成する。二次元フーリエ変換ル
ーチン８８が、濾波後の較正データ列メモリ８４₁、８
４₂、・・・８４_nのそれぞれの中のゼロ補充されたデ
ータを二次元フーリエ変換し、対応する低分解能複素画
像を生成して、対応する低分解能画像メモリ９０₁、９
０₂、・・・９０_nの中に記憶する。同じ二次元フーリ
エ変換ルーチンを直列に使用して、各データ列を処理す
ることも出来るし、あるいは複数のルーチンを備えて、
データ列を並列処理して高速化することも出来る。低分
解能複素画像はそれぞれ９２で正規化され、その大きさ
を１に等しく設定することが好ましい。一または複数の
複素共役ルーチン９４が、正規化された低分解能複素画
像のそれぞれの複素共役を計算するが、その画像の位相
は位相マップを表わす。このようにして、二次元位相補
正マトリックス９６₁、９６₂、・・・、９６_nが、ｋ
空間のｎ個の区分のそれぞれに、すなわち、エコー位置
のそれぞれに一個ずつ生成される。

【００１６】画像データの収集の間に、第１から第ｎエ
コーまでのそれぞれから集められた画像データ列は、４
８で画像データ列メモリ１００₁、１００₂、・・・、
１００_nの中にそれぞれ振分けられる。データ列は、そ
れぞれが表わすｋ空間の区分に対応するメモリの適当な
部分に位置される。ゼロ補充ルーチン１０２が、各メモ
リをゼロ補充して、画像データ・メモリ１００₁、１０
０₂、・・・、１００_nの中に正方マトリックスを生成
する。二次元フーリエ変換プロセッサ１０４が、二次元
的に各ゼロ補充された正方マトリックスをフーリエ逆変
換して複素部分画像を生成し、対応する複素画像メモリ
１０６₁、１０６₂、・・・、１０６_nの中に記憶す
る。画像乗算ルーチン１０８₁、１０８₂、・・・１０
８_nが、対応する一つのメモリ１０６_xからの複素部分
画像と、対応する一つのメモリ９６_xからの正規化され
た低分解能画像の複素共役とを、各点ごとに複素乗算す
る。このようにして、各エコーに対応する複素部分画像
が同じエコーに対応する位相マップに従って補正され
る。画像総計ルーチン１１２が、位相補正された部分画
像１１０₁、１１０₂、・・・、１１０_nを総計し、画
像を生成して画像メモリ１１４に記憶する。

【００１７】画像メモリ１１４からの合成画像は、一次
元で位相エンコード方向に沿って、プロセッサ１１６に
より、フーリエ逆変換され、データをｋ空間に位相エン
コードのためにのみ帰還させる。ｋ空間データは一時的
メモリ１１８に記憶される。密度補正プロセッサ１２０
は、合成ｋ空間画像データに対して、密度補正マトリッ
クス・メモリ１２２から検索した密度補正マトリックス
で各点毎に演算し、密度補正されたｋ空間データを生成
する。ｋ空間区分が重畳した時には、ｋ空間サンプリン
グ密度は、単純な正規化をすることにより補正出来る。
複素共役プロセッサ１２４は、密度補正されたデータで
演算して、欠落したデータをその複素共役データから得
て補充し、完全な画像を生成して画像メモリ１２６に記
憶する。プロセッサ１３０は位相エンコード方向に沿っ
て、一次元フーリエ変換を行ない、メモリ１２６の密度
補正された画像を変換して画像空間に戻し、密度補正さ
れた合成画像を生成し、出力合成画像メモリ１３２に記
憶する。ビデオ・プロセッサ１３４が、合成画像を適当
なフォーマットに変換して、ビデオ・モニタ１３６ある
いはその他の人間が読み取れる出力ディスプレイ装置に
表示する。

【００１８】画像データ列は、ｋ空間を通る異なる軌跡
により集められてもよい。例えば、図３および４におい
て、位相エンコード勾配５６が、ｋ空間の頂部で位相エ
ンコード・ステップに付けられるが、それは第一区分１
４０の頂部であり、かつ位相エンコード標点６４のそれ
ぞれが、位相エンコーディングをｋ空間の次の最低区分
１４２の頂部データ列へ、次の最低区分１４４へ、そし
て最低部の区分１４６へとステップさせて行く。その後
の反復で、位相エンコード勾配が、当初位相エンコード
・ステップを、最初の区分の次第に低くなるデータ列に
設定する。区分１４８中のデータは、区分１４６のデー
タに対するその複素共役関係により生成される。それに
代り、ｎ個の区分の一つから開始しないで、位相エンコ
ード勾配５６は、中心区分の一つの中の、一ステップを
選択してもよい。位相エンコード勾配標点６４は、位相
エンコーディングが予め選択された順序で、中心、正、
および負の位相エンコーディング・ステップ間で、ステ
ップについて外方にステップするように、様々な正、負
の極性と振幅のものでもよい。様々な他の映像シークエ
ンスが考えられるが、それらでは、各エコー位置は常に
ｋ空間の同じ区分に一致させられる。

【００１９】任意的に、位相エンコード勾配５６は、各
区分が隣接する次の区分に（図４には図示せず）数デー
タ列伸長するよう、位相エンコーディングを制御するた
めに使用してもよい。このようにして、数ラインの重畳
が生み出され、ｋ空間を完全に満たすよう確保され、か
つ区分の境界において起り得る位相の不連続の補正を助
ける。図５、６Ａ、６Ｂ、６Ｃにおいて、ｋ空間の中央
セクション１５２と、ｋ空間の同じ側にある二つのｋ空
間セクション１５４と１５６とのデータから形成される
複素画像１１４は、少なくとも列、すなわち垂直な方向
に、１５８でフーリエ逆変換され、一連の一次元フーリ
エ変換されたデータ列を生成し、それらが、ｋ空間のセ
クション１５２、１５４、１５６を満たすが、ここで１
５２はｋ空間の中央にある。ゼロ補充ルーチン１６０
は、ｋ空間の二つの空区分１６２と１６４をゼロ補充す
る。図６Ａに図示されるフィルタ１６６は、ｋ空間の中
央のデータ列以外は全て除去する。更に詳しくは、フィ
ルタ機能はセクション１５２の中央を通る単位値（unit
y value)と、セクション１５２の外部のゼロ値への円滑
な転移を生みだす。

【００２０】ｋ空間の中央セクションからのデータは、
１７０で列方向フーリエ変換をされ、低空間分解能画像
をｋ空間の中央セクションのデータから生成する。複素
低分解能画像は、その位相を維持しながら単位値の偏角
を有するように正規化される。このようにして、低分解
能画像位相マップが生成され、低分解能画像位相マップ
・メモリ１７２に記憶される。フィルタ１７４は、図６
Ｂのフィルタでデータを濾波する。図６Ｂのフィルタ
は、ｋ空間のセクション１５４と１５６のデータ全て
と、セクション１５２のデータの殆どを保持している。
しかしながら、セクション１５２のデータと、セクショ
ン１６２と１６４のゼロ補充されたデータとの間に平滑
な遷移が行なわれる。

【００２１】最初の反復では、一次元フーリエ変換プロ
セッサ１７６が、フィルタ１７４からのデータに、列方
向に一次元フーリエ変換を行ない画像を生成する。絶対
値ルーチン１７８が、画像から位相データを除去する。
乗算器１８０が、絶対値ルーチン１７８からの画像の各
画素の強度の値を、低分解能画像位相マップメモ１７２
からの対応する位相と結合し、複素画像を生成する。乗
算器１８０により形成された複素画像データは、１８２
で列方向に一次元フーリエ逆変換をされ、ｋ空間を補充
する行方向に変換されたデータ列の他の系列を生成す
る。この過程により、ｋ空間のセクション１６２と１６
４内のゼロデータ列は、第１近似データ列に置き換えら
れている。フィルタ１８４は、図６Ｃのフィルタでデー
タ列を濾波し、セクション１６２と１６４内のデータ以
外のデータの全ての部分を除去し、マトリックスの残り
の部分について、ゼロ値への平滑な遷移を達成させる。
フィルタ１７４と１８４との間には、データの完全な集
合が規定されていることに注目すべきである。これらの
二つのデータの部分集合が１８６で総計され、この過程
が反復される。すなわち、セクション１６２と１６４に
ついての第１近似データと、ｋ空間のセクション１５
２、１５４、１５６からの実データとを有する完全なデ
ータの集合が、１７６でフーリエ変換されて、複素画像
を規定する。位相情報は１７８で除去され、１８０で低
分解能画像位相マップメモリ１７２からの位相情報で置
き換えられる。その結果得られる複素画像が１８２で列
方向にフーリエ逆変換されると、ｋ空間のセクション１
６２と１６４についてのデータ列の第２近似値でデータ
の集合を再生成する。この過程はセクション１６２と１
６４についての第２近似データ列が第３近似値を生成さ
せる等々により反復される。結果として得られる画像が
最後に反復されたものとほぼ同じであれば、それが正確
なものであると決定される等により、セクション１６２
と１６４からのデータが、正確さで満足すべきものであ
れば、その結果得られる画像は最終画像メモリ１３２へ
と送られる。

【００２２】図７において、他の駆動平衡シーケンスも
構想される。図７では図３と同様な要素は同じ参照番号
を付されているが、ダッシュ（’）がついている。例え
ば、一個のα°の励振パルス５２’に、１８０°反転パ
ルス１９０が続いている。勾配エコー５８’、６０’の
タイミングは、１８０°再集束パルスと励振パルスに関
連して調整され、生成されたエコーの一つ６２_s’が、
スピン・エコーで残りは勾配・エコーであるようになさ
れる。一以上の追加の１８０°の再集束パルスが印加さ
れて、共鳴をそれに続くスピン・エコーに再集束させ
る。図８に図示ささるように、図３と同様な要素は二重
ダッシュ（”）がついた同じ参照番号を付されており、
１８０°再集束パルス１９２の数が、エコーの数に合う
ように増加されると、一連のスピン・エコー６２_s”が
生成される。任意的に、勾配パルスが、共鳴励振パルス
と第一の反転パルスとの間に生成されてもよい。較正デ
ータ列を生成させるために、勾配エコーを有利に使用
し、かつスピン・エコーに関する帯域幅を減少させるこ
とが出来る。

【００２３】図９において、図３と同様な要素は三重ダ
ッシュ（''')がついた同じ参照番号を付されており、位
相エンコード勾配２００が、「スライス選択」軸上に印
加されて、その方向に沿って位相をステップさせる。こ
の容積映像法では、ｋ空間は三次元矩形プリズムとして
想定される。標点６４''' は、ｋ空間中を位相エンコー
ド軸に沿ってステップする。スラブ選択勾配２０２が共
鳴励振パルスと同時に印加されて、励振をスラブ一個に
制限することが好ましい。発明の原理を更に明確に説明
するために、一般的に図３と４において、多重エコーを
有する二次元磁気共鳴画像について考慮する。スライス
選択励振５２の後に、その結果得られる水素原子核スピ
ン密度ρ（ｘ，ｙ）の三次元物体の磁気共鳴磁化誘導信
号は、下記の式によって与えられる。

【００２４】

【数１】

【００２５】ここで、ｑはエコーの数、ｋ_rd(t) とｋ_pe
(t) は読み取りと位相エンコーディング勾配のｋ空間軌
跡をそれぞれ表わす。

【００２６】

【数２】

【００２７】ここで、Ｇ_rdとＧ_peは信号読み取り時の時
間依存磁場読み取り勾配と位相エンコーディング勾配
で、△Φ（ｘ，ｙ）は磁石の主磁場の磁化率の両者によ
る影響を含む主磁場の不均一性による位相誤差を表わ
す。横方向スピン緩和はＴ₂を含む項により考慮に入れ
られている。勾配うず電流の磁界効果は、計算を簡略化
するためにこの式には含まれていない。スライス選択信
号取り込み（slice selection acquisition)速度を早め
るには、図３、５、６、７に見られるように、振動読み
取り勾配あるいは多重高周波集束パルスを用いて、いく
つかの勾配および／またはスピン・エコーが、各スピン
励振毎に得られる。位相補正に用いられる較正データ列
は、全てのエコーと全ての空間スライスについて得られ
る。較正データ列は画像データと較正データとの間の解
剖図の動揺の可能性を最小にするために、区分化された
画像取り込みの前、後、あるいはその間に挟まれて得る
ようにしてもよい。位相エンコーディング標点６４につ
いて、ゼロに近い振幅での映像と同じシークエンスを用
いて、ｋ空間１８０の原点の周囲の一組の共通位相エン
コーディング・ステップのもとに、各エコー毎に較正デ
ータが得られる。位相エンコーディング・ステップの数
は、比較的少なくてもよく、あるいは区分の一つ１６２
のサイズでもよい。もし位相エンコーディング・ステッ
プが、区分の一つのサイズである場合には、較正データ
は画像データとしても使用可能で、その区分から再びデ
ータを集める必要がなくなる。

【００２８】様々なデータ時点において得られた、異な
るエコーから取り込まれた複数の位相エンコードされた
データ列を結合して、一つのｋ空間データ領域とし、ア
ーティファクトのない画像を形成するためには、隣接す
るｋ空間区分の間の位相不連続性を最小化するために、
下記に説明する位相補正がデータについて行なわれる。
二次元フーリエ逆変換で画像を再構成するには、異なる
エコーからのデータは、通常読み取りおよび位相の両方
向に、タイム・シフトあるいは位相補正される必要があ
る。各エコーについての磁場誘導位相誤差マップＣ
（ｘ，ｙ，ｑ）は、各エコーについての一組の二次元較
正データから、ｋ空間の中央の周囲に対称的に少数の位
相エンコーディング・ステップを用いて決定される。

【００２９】

【数３】

【００３０】ここで、ｒ（ｔ，ｖ，ｑ）は、視野ｖとエ
コーｑのサンプリング時間ｔの関数としての生データ列
を表わす。

【００３１】

【数４】

【００３２】ここで、ｆ_l（ｔ，ｖ）は、位相と読み取
りの両方向に沿った、二次元平滑フィルタ関数を表わ
し、Φ（ｘ，ｙ，ｑ）は、エコー時ｑでの、磁場誤差に
ついての低空間分解能位相推定（low spatial resoluti
on phase estimate)を表わす。較正データのための視野
の数は、いくつでもよいが、８、１６、または３２が好
ましい。エコーｑの低空間分解能位相誤差を表わす単位
複素数的二次元画像マトリックスが、較正データの第ｑ
エコーから下記のように得られる。

【００３３】

【数５】

【００３４】収集した映像データ列は、サンプリング時
間内に先ずフィルタされ、それらが受取られた順序で一
次元的にフーリエ逆変換される。この段階で一次元的に
フーリエ変換されたデータ・ラインは、ｋ空間の連続す
る一組の区分を形成する同じエコー・インデックスに対
応するデータ列に選別される。インデックス整数ｑによ
り表わされるあるエコー区分が完全に満たされると直に
その区分は取出され、対応する視野位置の二次元複素画
像データ・バッファ内に置かれる。残余のエコー区分
は、それから１０２によりゼロ補充される。一次元的な
列方向のフーリエ変換が行なわれる。これに代り、二次
元フーリエ変換が、ゼロ補充の後に行なわれてもよい。
二次元フーリエ変換を使用する場合には、エコーｑの中
間的画像が下記のように得られる。

【００３５】

【数６】

【００３６】ここで、Ｖ_segは各ｋ空間の区分の幅で、
Ｖ_Oはその区分の中央視野インデックスを表わし、２Ｄ
−ＦＴはデータ・サンプリング時間（ｔ）と視野（ｖ）
方向に沿った、離散的二次元フーリエ変換を表わし、Π
（ｘ）はｋ空間をエコー位置または時間により複数区分
に分割する標準的矩形ウインドウ関数を表わす。このよ
うにして、代表的４エコー複数ショットＥＰＩ撮影の各
画像スライスは、対応する４個のｋ空間区分から４個の
中間画像を生みだす。これらの中間二次元画像は、複素
画像メモリ１０６に記憶される。二次元位相補正は、各
ｋ空間データ区分に対応する各中間画像について独立し
て行なわれる。各部分またはサブ（sub)画像Ｒ（ｘ，
ｙ，ｑ）は画像領域で下記により補正される。

【００３７】

【数７】

【００３８】ここで、^*は複素共役演算を表わす。最終
補正済み複素画像は、１１２で異なるエコーに対応する
全ての個別位相補正済み中間画像を総合することで生成
される。

【００３９】

【数８】

【００４０】ここで、１ｍ（ｘ，ｙ，Ｅ_n）は、疑似エ
コー時間Ｅ_nの最終画像を表わし、Ｅ _nは使用された中
央ｋ空間データのエコー時間を表わす。どの二つの隣接
するｋ空間区分も、二次元位相方向に、間隙があっても
重複してもよい。これは二次元位相補正後、ｋ空間の軌
跡を検証して決定される。この結果生じるｋ空間サンプ
リング密度の変化は、ｋ空間データ中で容易に補正され
る。周波数エンコーディングのための読み取り勾配パル
スの波形は、台形でも、正弦波でも、あるいはその他の
適当な波形でもよい。データは読取り勾配の発振期間全
体を通して、連続的にサンプルすることが出来る。もし
サンプル・データがｋ空間で等距離でなければ、データ
は等距離離散的表示に再等間隔化されることが望まし
い。

【００４１】ガラス、ハニング、ブラックマン等の適当
な低域フィルタが、生データの信号ピーク位置を中央に
位置させる。低域フィルタは収集されたデータの性格に
最適化された、ｎ項コサイン拡張形式のような一般的な
ものでよい。例えば、副ローブ出力は、分解能に大きな
損失なしに最小化出来る。較正データはこのフィルタで
平滑化され、決定される位相値でのノイズ・レベルを減
少させる。その結果生じる補正二次元位相アレイは、空
間的ＤＣ、線形、および高次（high order) 成分を含ん
でいる。離散的形式では、画素位置（ｉ，ｊ）での単位
複素数補正ベクトルは、下記のとおりである。

【００４２】

【数９】

【００４３】ここでΦ_ijは一組の整数（ｉ，ｊ）により
表示された画素位置での位相誤差を表し、一般的に空間
内で比較的平滑な関数である。位相誤差は通常相当に平
滑な関数であるので、複素数マトリックスの間引きされ
た一組のみを保存すればよい。マトリックスは乗算段階
の前に、補完によって容易に拡張出来る。全てのデータ
区分を総合して実効総ｋ空間データ・サンプリング密度
は、全ての区分についての寄与の総計によって得られ
る。

【００４４】

【数１０】

【００４５】ここで、異なるｋ空間の区分は、当初のデ
ータ取り込みにおいて互いに重なり合ってもよい。実際
でのｋ空間におけるこの不均一サンプリング密度に対す
る補償のために、ｋ空間で二次元位相補正済みサブ画像
が、結合されるや否や、結合された画像は画像空間間か
らｋ空間へと、フーリエ逆変換されて戻される。すなわ
ち：

【００４６】

【数１１】

【００４７】ここで、スーパースクリプト「^com」は結
合された画像を、ｑは区分、すなわちエコー番号インデ
ックスを表わす。本記述での約束として、Ｒは画像表現
を、ｒは生表現を表わすものとする。結合されたｋ空間
データの信号強度は、各点毎にサンプリング密度に従
い、尺度（スケール）が直される。すなわち：

【００４８】

【数１２】

【００４９】最終の補正された複素画像が、標準的二次
元フーリエ変換を使用して、再構成される。すなわち

【００５０】

【数１３】

【００５１】強度と位相画像は、単に複素画像の絶対値
と位相値とを採ることで得られる。すなわち

【００５２】

【数１４】

【００５３】これに代り、ｋ空間中の不等距離サンプリ
ング・データに対する標準的等間隔化法が、区分された
ｋ空間データの結合と補正に使用出来る。較正エコー・
データは、画像再構成について所望のＴ₂重みずけを実
現するために使用出来る。位相アレイ・コイルと共に使
用される場合には、較正データ・ラインは、全てのＲＦ
受信チャネルから同時に得られる。個別の受信コイルに
対応する画像は別々に補正され、再構成され、最後で結
合されて、最終合成画像を生成する。較正データと映像
データとは共に、鋭い上昇および下降ランプの台形読み
取り勾配の間に集めることが出来る。これに代り、正弦
波形の読み取り勾配が使用出来るが、但し、周波数方向
での非直線性にたいする対応する調整をしなければなら
ない。

【００５４】部分フーリエ再構成についての反復法で、
一次元において画像再構成のために、部分データのみし
か得られないものとすれば：

【００５５】

【数１５】

【００５６】ここで、ｓ（ｋ）は部分位相エンコーディ
ング（ＰＥ）視野について規定されているのみで、数学
的には下記のように表される。

【００５７】

【数１６】

【００５８】反復法は三つの段階を含んでいる。先ず中
央対称部分データで、完全サイズへのゼロ補充とフィル
タされたフーリエ変換により、低周波数複素画像が得ら
れ、それから低周波数位相マップが構成される。画像再
構成において起り得る周辺部でのぼやけや、リンギング
を減少させるために、位相エンコーディング方向に沿っ
てフィルタ（図６Ｂ）が、下記のように定義される。

【００５９】

【数１７】

【００６０】ここで、ｎ_wはフィルタの遷移領域の幅を
表わしている。低空間分解能画像は先ず、生ｋ空間デー
タの中央部分から下記のように決定される。

【００６１】

【数１８】

【００６２】同時に、信頼出来る低周波数位相マップが
規定され、それが位相制約として使用される。この位相
は主磁場不均一性の位相誤差の尺度である。第二に反復
再構成のための最初の画像が、下記のように生成され
る。

【００６３】

【数１９】

【００６４】ここで、フィルタ関数ｗ_o（ｋ）（図６
Ａ）は、下記となる。

【００６５】

【数２０】

【００６６】第三に、反復は下記のように行なわれる。
第ｑ反復段階で、新しい複素画像が、その前のものから
下記のように構成される。

【００６７】

【数２１】

【００６８】逆ＦＴを取れば、新しい生データが下記の
ように取られる。

【００６９】

【数２２】

【００７０】次の反復の新しい生データが下記のように
構成される。

【００７１】

【数２３】

【００７２】ここで、ｗ（ｋ）は、重みづけ関数（図６
Ｃ）を表わし、それは下記のように規定される。

【００７３】

【数２４】

【００７４】ここで、ｎ_tは混合関数（mixing functio
n)についての、遷移領域の幅を規定する。

【図面の簡単な説明】

【図１】図２とともに本発明による磁気共鳴映像システ
ムを図示する略線図である。

【図２】図１とともに本発明による磁気共鳴映像システ
ムを図示する略線図である。

【図３】好ましいエコー平面映像シーケンスの一回分を
示す図である。

【図４】ｋ空間区分の代表的補充例を示す図である。

【図５】部分的フーリエ再構成のための反復法を示す略
線図である。

【図６】図６Ａ、図６Ｂ、図６Ｃは、図４の方法で使用
されるフィルタの関数を示す図である。

【図７】好ましい勾配スピンエコー（ＧＳＥ）映像シー
ケンスを示す図である。

【図８】好ましい急速スピンエコー（ＦＳＥ）映像シー
ケンスを示す図である。

【図９】好ましいエコー容積映像シーケンスを示す図で
ある。

【符号の説明】

１０主磁場制御器１２超伝導または抵抗性磁石１４検査領域２０勾配パルス増幅器２２全身勾配コイル２４デジタル無線周波数送信器２６全身ＲＦＲコイル３０挿入可能な頭部コイル３２局所勾配コイル３４局所無線周波数コイル３６ＲＦスクリーン３８受信器４０シーケンス制御回路４２Ａ−Ｄコンバータ４４選別器

───────────────────────────────────────────────────── フロントページの続き (72)発明者ジェームズビーマックナリーアメリカ合衆国オハイオ州 44022 シャグリンフォールズラークスパーレーン 7919 (72)発明者ゴードンディーディメースターアメリカ合衆国オハイオ州 44092 ウィクリフデニスドライヴ 1613

Claims

【特許請求の範囲】

【請求項１】検査領域（１４）を通る時間的に一定な
磁場を発生させる磁石（１２）と、前記検査領域にある
双極子に共鳴を誘起して、無線周波数共鳴信号を発生さ
せるための無線周波数パルス制御器（４０）と送信器
（２４）と、少なくとも前記検査領域を直角方向に横断
する位相および読取り磁場勾配パルスを生成させ、かつ
前記読取り勾配を反復して反転させる勾配磁場コイル
（２２）と勾配磁場制御器と、前記読取り勾配が反転す
る毎にその後、前記無線周波数磁気共鳴信号を受信して
復調し、一連のデータ列を作りだす受信器（３８）と、
再構成された画像表現を記憶する画像メモリーを有して
いる磁気共鳴映像システムにおいて、前記読取り勾配の反転の一つに各々対応する複数の二次
元位相補正マトリックスを生成し、この生成が、前記読
取り勾配の各反転後のデータ収集期間に収集されたデー
タ列を受信して各データ収集期間に対応して前記二次元
位相補正マトリックスを独立に生成することにより行わ
れる二次元位相補正プロセッサ（８８、９０、９２、９
４）と、各前記データ収集期間中に生成された画像データを受信
し、各前記データ収集期間に収集された一組の不完全な
データ列を対応する部分画像に独立に再構成する部分画
像プロセッサ（１０４）と、各部分画像を同じ前記データ収集期間に対応する前記二
次元位相補正マトリックスに乗算し複数の二次元的に位
相補正された部分画像を生成するマトリックス乗算器
（１０８）と、前記二次元的に位相補正された部分画像を加算して合成
画像を生成し前記画像メモリ（１１４）に記憶する画像
加算器（１１２）を有することを特徴とする磁気共鳴映
像システム。
【請求項２】各データ収集期間に対応する較正データ
列を処理して位相が位相誤差マップを表す対応する低空
間分解能複素画像を生成させるフーリエ変換プロセッサ
と、各低空間分解能複素画像を正規化するためのプロセ
ッサと、各正規化された低分解能複素画像の複素共役を
計算して各データ収集期間に対応する前記二次元位相補
正マトリックスを生成させるための複素共役プロセッサ
を有する前記二次元位相補正プロセッサを更に備えるこ
とを特徴とする請求項１に記載の磁気共鳴映像システ
ム。
【請求項３】ｋ空間のある区分に対応する各データ収
集期間に不完全なデータの組が収集されるようにデータ
列を受信する前記部分画像プロセッサと、更に収集され
たデータがないｋ空間の部分をゼロで満たして前記部分
画像にと再構成するための完全なデータの組を生成させ
るゼロ補充手段とを更に有することを特徴とする請求項
１に記載の磁気共鳴映像システム。
【請求項４】各点毎に前記合成画像を密度補正マトリ
ックスで補正し密度補正済み合成画像を生成させるため
のｋ空間サンプリング密度補正プロセッサを更に有する
ことを特徴とする請求項１に記載の磁気共鳴映像システ
ム。
【請求項５】ｋ空間の中央区分と一方の側にある第一
の区分に隣接するｋ空間の第二の区分との中の一連のデ
ータ列を規定するために前記合成画像をフーリエ逆変換
するフーリエ変換プロセッサと、ｋ空間の前記中央セク
ションに他方の側で隣接する第三のｋ空間区分をゼロ補
充するためのゼロ補充ルーチンと、ｋ空間の前記中央区
分から前記データ列を抽出するためのフィルタと、前記
中央ｋ空間データ列を低分解能画像位相マップにとフー
リエ変換するフーリエ変換プロセッサと、ｋ空間の中央
と第二の区分とからの前記データ列を通過させるフィル
タと、前記第二のフィルタからのｋ空間の前記中央と第
二の区分とからの前記データ列を変換するための第二の
フーリエ変換プロセッサと、前記第二のフーリエ変換プ
ロセッサからの画像表現の大きさを決定する絶対値ルー
チンと、前記絶対値ルーチンからの対応する大きさを前
記低分解能画像位相マップからの位相と結合し複素画像
を生成する乗算器と、前記複素画像をフーリエ変換しｋ
空間の前記中央と第二と第三の区分中のデータ列を含む
一組のデータ列を生成するフーリエ逆変換手段と、ｋ空
間の前記第三区分からの前記データ列を通過させる第三
のフィルタと、前記第二と第三のフィルタからの前記デ
ータ列を結合し、前記第二のフーリエ変換プロセッサに
前記データ列を供給する加算器とを有し、それにより前
記複素画像が前記乗算器からの前記複素画像により反復
して改良され選択された回数の反復の後、最終画像とさ
れることを更に特徴とする請求項１に記載の磁気共鳴映
像システム。
【請求項６】双極子中に磁気共鳴が励起され、かつ一
連の磁気共鳴エコーが誘起され、エコーを読取り勾配の
存在中に読み取り軸に沿って読み取って一連のデータ列
を形成し、出力画像表現として再構成する磁気共鳴映像
法において、前記エコーを位相エンコードしてｋ空間の
複数の区分のそれぞれに沿っているデータ列を生成し、
ｋ空間の前記複数の区分の間に前記データ列を振分け、
前記データ列から位相補正マトリックスを生成し、各マ
トリックスはｋ空間の前記区分の一つに対応し、前記デ
ータ列から複数の部分画像を生成し、各部分画像はｋ空
間の前記区分の一つに対応し、複数の部分画像と位相補
正マトリックスから複数の位相補正された部分画像を生
成させ、位相補正された部分画像のそれぞれはｋ空間の
前記区分の一に対応し、かつ前記位相補正された部分画
像を総計して合成画像を生成することを特徴とする磁気
共鳴映像方法。
【請求項７】前記位相補正マトリックスを生成する段
階において、ｋ空間の各前記区分に対応する位相補充さ
れた較正データ列をフーリエ変換して一組の低分解能複
素画像を生成させ、かつ各低分解能画像を正規化して前
記位相補正マトリックスを生成させることを更に特徴と
する請求項６に記載の磁気共鳴映像方法。
【請求項８】前記位相補正マトリックスを生成する段
階において、ｋ空間の各区分について、ｋ空間の空部分
をゼロ補充してゼロ補充された位相エンコードされた較
正データの組を生成させ、空部分はｋ空間の他の区分に
対応し、各ゼロ補充されたデータの組をフーリエ変換し
て前記位相補正マトリックスの一つを生成させることを
特徴とする請求項７に記載の磁気共鳴映像方法。
【請求項９】前記合成画像を密度補正マトリックスで
演算して密度補正された合成画像を生成させることを更
に特徴とする請求項６に記載の磁気共鳴映像方法。
【請求項１０】（ａ）少なくとも列方向に前記合成画
像をフーリエ逆変換し、一組のデータ列を生成させ、そ
のデータ列の組がｋ空間の中央部分と前記中央区分の一
方の側にあるｋ空間の第二の区分とからのデータ列を含
み、（ｂ）前記中央区分の他方の側にあるｋ空間の第三
の区分をゼロ補充し、（ｃ）前記データ列をフィルタし
てｋ空間の前記中央区分からの前記データ列を抽出し、
前記抽出されたデータ列をフーリエ変換し、低分解能画
像を生成し、前記低分解能画像を正規化して低分解能画
像位相マップを生成させ、（ｄ）前記データ列をフィル
タしてｋ空間の前記中央と第二の区分とからの前記デー
タ列を抽出し、（ｅ）前記抽出されたデータ列をフーリ
エ変換して画像を生成し、大きさ画像を生成するために
前記画像の絶対値を取り、前記大きさ画像と前記低分解
能位相マップとを結合複素画像を生成し、（ｆ）少なく
とも前記列方向に前記複素画像をフーリエ逆変換し、ｋ
空間の前記中央と第二と第三の区分にデータ列を生成さ
せ、かつ前記データ列をフィルタして、ｋ空間の前記第
三区分からの前記データ列を抽出し、（ｇ）ｋ空間の前
記中央と第二の区分とから抽出された前記データ列をｋ
空間の前記第三区分から抽出された前記データと結合
し、（ｅ）段階と（ｆ）段階を反復して前記生成された
複素画像が反復して改善されることを更に特徴とする請
求項６に記載の磁気共鳴映像方法。