JPH03295537A - Ｍｒｉ装置 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 （産業上の利用分野） 本発明はＭＲＩ（磁気共鳴画像撮像装置）の３次元イメ
ージング方法に関し、更に詳しくは３次元スキャンにお
いて発生するエリアジングイメージを分離除去し、有効
なイメージを効率的に得ることのできるＭＲＩの３次元
イメージング方法に関する。

（従来の技術） 原子核を静磁場中におくと、原子核は磁界の強さと原子
核の種類によって異なる定数に比例した角速度で歳差運
動をする。この静磁場に垂直な軸に前記の周波数の高周
波回転磁場を印加すると磁気共鳴が起こり、前記定数を
有する特定の原子核の集団は共鳴条件を満足する高周波
磁場によって準位間の遷移を生じ、エネルギー準位の高
い方の準位に遷移する。共鳴後高い準位に励起された原
子核は低い準位に戻ってエネルギーの放射を行う。

ＭＲＩはこの特定の原子核によるＮＭＲ（核磁気共鳴）
現象を観察して被検体の断層像を撮像する装置である。

ＭＲＩにおいてフーリエ変換法に用いる高周波磁場及び
勾配磁場印加のパルスシーケンスを第８図に示す。図に
おいて、（イ）図はそれぞれリード軸、ワーブ軸、スラ
イス軸であるｘｒ　　Ｙｒ　　Ｚ軸にＧｘ、Ｇｙ、Ｇｚ
の勾配磁場を与え、高周波磁場をＸ軸に印加する状態を
示す図で、（ロ）図はそれぞれの磁場を印加するタイミ
ングを示す図である。期間１において、９０°パルス１
とスライス勾配２によりｚ−０を中心とする２方向に垂
直なスライス面内のスピンが選択的に励起される。

期間２のリフェーズ勾配３はスライス勾配２により乱れ
たスピンの位相を元に戻すためのものである。同じ期間
２のデイフェーズ勾配４はデータ読み出し期間４の時間
的中心にＳＥ信号７の中心が一致するようにスピンに場
所に応じた位相差を与えるためのものである。期間２で
は更にｙ方向の位置に比例してスピンの位相をずらせて
やるためのワーブ勾配６を印加しており、ワーブ勾配６
は毎周期その強度を変えて印加されている。その後１８
０”パルス５を与えて磁気モーメントを揃え、その後に
現れるＳＥ信号７を観察する。期間４ではＸ軸にリード
勾配８を印加する。これにより、デイフェーズ勾配４で
与えられた位相差は、期間４のリード勾配８の時間的中
心で相殺されＳＥ信号７が現れる。このシーケンスをビ
ューといい、パルス繰り返し周期ＴＲ後に再び９０°パ
ルス１を加えて、次のビューを開始する。

（発明が解決しようとする課題） 上記のパルスシーケンスは２次元画像データを採取する
ためのスピンエコー法パルスシーケンスであったが、３
次元画像データを取るために従来第９図のようなパルス
シーケンスで行っていた。

図において、第８図と同等な部分には同一の符号を付し
である。図中、９はスライス軸に印加した第２ワーブ勾
配で、第１ワーブ勾配６と同様に毎ビューその強度を変
えて印加されている。そして、第１ワーブ勾配６は２次
元平面を掃引し、第２ワーブ勾配９はスライス方向即ち
厚み方向を掃引して、多数の薄いスライス像を作ってい
るが、読み出し方向へはＬＰＦ等を用いて帯域制限をす
ることにより必要な範囲外のデータを取り除くことがで
きるが、ワーブ方向ではデータ収集をしながらハードウ
ェアで帯域制限をすることができないため、必要範囲外
のデータが廻り込んでエリアジングとなってしまう。こ
れを除くために第１０図や第１１図に示すような方法を
取ることにより、余分な領域のデータを取り除いている
。第１０図はＲＦの感度分布を利用した例である。図に
おいて、１４はＲＦコイルの感度分布を示す曲線で、コ
イルの端末における感度低下を利用して余分な領域のデ
ータを取り除こうとした状態を示しているが、ＲＦコイ
ルの感度分布曲線１４はなだらかなため必要な第２ワー
プ方向視野（以下ＦＯＶという）１５の領域内に、領域
外のデータが廻り込んで、両端部がエリアジング１６で
潰されてしまう。第１１図はＲＦの選択励起による帯域
制限を利用した例である。図において、第１０図と同等
の部分には同一の符号を付しである。１７はＲＦの選択
励起によるスライス形状である。この場合はスライス形
状１７を矩形に切り出すことができればエリアジング１
６を取り除くことができるが、矩形に切り出すことは難
しいという問題がある。

又、第２７−ブ方向のマトリックス（スライス画像枚数
）とＦＯＶを大きくしてデータを得ることによりエリア
ジングを除くこともできるが、スキャン時間が増大し、
更に、不必要なデータを取ることになって効率が低下し
てしまう。

本発明は上記の点に鑑みてなされたもので、その目的は
、スキャン時間の増大、分解能の劣化、ＳＮ比の劣化等
を生ずることなくエリアジングイメージを分離除去し、
効率良く有効なイメージを得るＭＲＩの３次元イメージ
ング方法を実現することにある。

（課題を解決するための手段） 前記の課題を解決する本発明は、第１ワープ勾配をワー
プ軸に印加し、第２ワープ勾配をスライス軸に印加して
スキャンを行うＭＲＩの３次元イメージング方法におい
て、使用者指定の厚み方向の視野ＦＯＶｚ、厚み方向マ
トリスクＮ２及びスキャン平均回数Ｎａを入力する条件
設定の段階と、実スキャンにおいて厚み方向視野を増加
してエリアジングイメージングを有効視野がら追い出す
ために次式により実スキャン厚み方向視野ＦｏＶ２、実
スキャン厚み方向マトリクスＮ２　、スキャン平均回数
Ｎａ／及びＮＡ　、厚み方向部分平均ビュー数Ｎｆｒを
演算により求める段階と、ＮＺ　　−ＩＮＴ（ＮＺ　”
ｃ／２．０＋０．９９９）ｘ２ＣはＣ≧１．２ ＦＯＶＺ　　−ＦＯＶｚ　Ｘ　ＮＺ　’　／　　ＮＺＮ
ａ　　＝　　Ｎ−ＮＺ／　　ＮＺ ＮＡ　−ＩＮＴ（Ｎａ　’　） Ｎｔ、＝　　Ｎａ　　　ＮＺ−ＮＡ　　　　ＮＺ実スス
キャン第２ワーブステツプＧＷＳ２次式により求めて決
定し、そのＧＷＳ２の間隔を持つワープ勾配を印加する
段階と、 ＧＶＳ２−１バγ・ＴＳｗ・ＦＯｖｚ′）ＴＳＷ−・ワ
ープ時間 厚み方向マトリクスＮＺ　　＋　スキャン平均回数Ｎ９
′及び厚み方向部分平均ビュー数Ｎａ、に基づいて厚み
方向マトリクスＮ２′−で部分的に異なる実スキャン平
均回数ＮＡ′のスキャンを行う段階と、厚み方向マトリ
クスＮＺ′におけるデータの離散フーリエ変換演算を行
ってデータ処理を行う段階と、前記段階においてフーリ
エ変換を行ったデータから不必要なエリアジングデータ
を捨て有効データのみを保存する段階と、前記段階で得
た有効データにより２次元画像を再構成する段階とから
成ることを特徴とするものである。

（作用） 指定された厚み方向の視野、マトリクス数及びスキャン
平均回数を入力し、その入力データに基づき一定の法則
に従って視野を拡大し、拡大した視野に応じたマトリク
ス数を決定し、合計のマトリクス数を指定された厚み方
向マトリクス数とスキャン平均回数によるマトリクス数
に合わせるように残りのマトリクス数を決める。

（実施例） 以下、図面を参照して本発明の方法を詳細に説明する。

第２図は本発明の方法を実施するＭＲＩのパルスシーケ
ンスの図である。このパルスシーケンスを説明する前に
スキャンの時間短縮のための手法を説明する。第３図は
第２図のパルスシーケンスによって行おうとする手法の
説明図である。図において、（イ）図はデータを取ろう
とする被検体１０のスライス方向の厚みを示す図で、こ
の厚みを２５６ビユーの走査でデータを採取するもので
ある。（ロ）図は（イ）図の被検体１０の厚みを５等分
して、１０ａ、１０ｂ、１０ｃ、１０ｄ。

１０ｅに分け、１０ｂ、１０ｄは関心の無い領域として
データを取ることなく、１０ａ、ＩＯＣ。

１０ｅの部分に対して各５１ビユーずつでデータを採取
する場合を示している。この１０ｇ、１０ｃ、１０ｅの
部分は同−ＴＲ時間内に３回の走査をそれぞれ行って走
査時間を１１５にしようとするものである。

次に、第２図のパルスシーケンスについて説明する。図
において、第８図、第９図と同等の部分には同一の符号
を付しである。図中、１１は９０°パルス１を印加した
時スライス勾配２によって特定の部位（以下スラブとい
う）、例えば、第３図の１０ａを励起し、次に１８０°
パルス５を印加して反転させる時、他のスラブＩＯＣ，
１０ｅにも及ぶ範囲を励起すれば、その回復を待つ必要
があるので、他のスラブに影響を及ぼさないようにして
第１のスラブ１０ａのみに限定するためのスライス勾配
である。ここでスラブというのはデータを取る厚み範囲
の中、関心のない領域１０ｂ。

１０ｄはデータを取らないで、関心領域１０ａ。

１０ｃ、１０ｅのデータを取るようにした各範囲１０ｇ
、１０ｃ、１０ｅのことである。１２はスライス勾配置
１により乱れたスピンの位相を元に戻すためのリフェー
ズ勾配である。

２１は同−ＴＲ内に第２スラブ１０ｃを励起するための
９０１′パルスで、第１スラブ１０ａを励起する９０°
パルス１とは周波数が異なっている。

２２は９０°パルス２１と共に第２スラブ１０ｃを共鳴
させるためのスライス勾配である。２４〜３２はそれぞ
れ４〜１２の第１スラブ１０ａに対する各ＲＦパルス勾
配、信号等に相当する第２スラブ１０ｃに対するもので
、２〜１２に対して２０を加えた番号のものか対応して
いる。第２スラブ１０ｃの後に、第３スラブ１０ｅに対
してＲＦパルス、勾配等を印加するが、前記の場合と同
様なので省略する。この場合も９０″パルス、１８０°
パルスの周波数は第３スラブ１０ｅの位置を指定する周
波数となっている。

３Ｄスキヤンはこの各スラブについて行うもので、第１
スラブ１０ａについて見ると、第４図に示すようにこの
スラブ１０ａの範囲内を複数のスライス面３０で構成さ
れるように３次元データを得るものである。図では６枚
のスライス面３０によるデータを得る場合を示している
。

次に上記のように構成された実施例のパルスシーケンス
によるスキャンを実施するためのＭＲＩの要部構成図を
第５図に示す。

図において、４１は内部に被検体を挿入するための空間
部分（孔）を有し、この空間部分を取巻くようにして、
被検体に一定の静磁場を印加する静磁場コイルと勾配磁
場を発生する勾配磁場コイル（勾配磁場コイルはＸ＊　
　Ｖ＋　　Ｚの３軸のコイルを備えている。）と被検体
内の原子核のスピンを励起するためのＲＦパルスを与え
るＲＦ送信コイルと被検体からのＮＭＲ信号を検出する
受信コイル等が配置されているマグネットアセンブリで
ある。静磁場コイル、勾配磁場コイル、ＲＦ送信コイル
及び受信コイルは、それぞれ静磁場電源４２、勾配磁場
駆動回路４３、ＲＦ電力増幅器４４及び前置増幅器４５
に接続されている。シーケンス記憶回路４６は計算機４
７からの指令に従って任意のビューで、ゲート変調回路
４８を操作（所定のタイミングでＲＦ発振回路４９のＲ
Ｆ出力信号を変調）シ、第１図のパルスシーケンスに基
づくＲＦパルス信号をＲＦ電力増幅器４４からＲＦ送信
コイルに印加する。又、シーケンス記憶回路４６は、同
じく第１図のパルスシーケンスに基づくシーケンス信号
によって勾配磁場駆動回路４３を操作して、第３図に示
すようにｘ、　　ｙ、　　ｚの３軸にそれぞれ勾配磁場
を供給する。５０はＲＦ発振回路４９の出力を参照信号
として、前置増幅器４５の受信信号出力を位相検波する
位相検波器である。

この出力信号はＡＤ変換器５１においてディジタル信号
に変換され、計算機４７に入力される。５２は計算機４
７に種々のパルスシーケンスの実現のための指示及び種
々の設定値等の入力をするための操作コンソール、５３
は計算機４７で再構成された画像を表示する表示装置で
ある。

次に、上記のように構成された装置を用い、第２図のパ
ルスシーケンスで行われる本実施例の３次元イメージン
グ方法の手順を第１図のフローチャートを用いて説明す
る。

ステップ１ 使用者指定の厚み方向の視野（スラブの幅に相当）ＦＯ
Ｖｚ、厚み方向のスライス枚数であるマトリクスＮＺ（
≧４の偶数）、第１スラブ１０ａのデータ採取のために
行うスキャンの繰り返しの平均回数Ｎａを操作コンソー
ル５２により計算機４７に入力する。

ステップ２ エリアジングイメージが有効スライス以外の厚み方向の
スライス位置に分離できるように厚み方向の視野、厚み
方向のマトリクスを大きく取る。

この場合、厚み方向の１枚当たりのスライス厚を変化さ
せないように厚み方向の視野とマトリクスを決定する。

このように大きく取った視野をＦＯＶＺ′、′　とし、
マトリクスをＮ２′　とする。このステップではステッ
プ１で入力された各データに基づき計算機４７は次式の
計算を行い、マトリクスＮ２′、視野ＦＯＶｚ’、スキ
ャン平均回数Ｎ。

′を求める。平均回数Ｎａｒ　は全範囲のスキャン回数
即ち整数スキャン回数ＮＡ′と厚み方向部分平均ビュー
数Ｎｆｒとに分けられる。

Ｎ２−ＩＮＴ（ＮＺ　”ｃ／２’、０＋０．９９９）ｘ
２−（１）式中、ＩＮＴは整数化を意味し、０．９９９
を加えて整数化するのは切り上げを意味している。又、
ＩＮＴ中で１／２にし整数化後２倍するのは０倍した範
囲内の偶数マトリクスの数を得るためである。尚、平均
回数ＮａはＮａ≧１．５で指定される。又、Ｃは厚み方
向マトリクスＮＺに対する倍数で、Ｃ≧１．２である。

ＦＯＶｚ　’　−ＦＯＶｚ　Ｘ　ＮＺ　’　／　　ＮＺ
　　　−（２）この式は（１）式によって増やされたマ
トリクスに対応して増加したＦＯＶである。

スキャン時間を使用者指定のものと等しくするために実
スキャンパラメータとなるＮ２′で行う平均回数ＮａＬ
を決定する。

Ｎ−”　　Ｎ−ＮＺ／　　ＮＺ’　　　　・・・（３）
（３）式で決定されるＮ　％から実スキャンにおいて全
厚み方向のビューに亘るスキャンの平均回数ＮＡ′及び
部分的に平均を行う厚み方向のビュー数Ｎ　ＩＩを求め
ると、 ＮＡ　　−ＩＮＴ（Ｎ−’　）　　　　　　　　・・・
（４）（４）式のＩＮＴは小数部分を切り捨てた整数部
分を意味する。

Ｎａ、−Ｎａ　　　ＮＺ−Ｎ＾′ＮＺ　　・・・（５）
Ｎｌ＋はＦＯＶＺ′、′の整数部分を除いた部分的なス
キャンの平均マトリクス数である。以上の演算はすべて
計算機４７て行われる。

ステップ３ （１）式、　（２）式で得られたマトリクスＮＺ′及び
ＦＯＶＺ’　を実スキャンのパラメータとして厚み方向
のワープステップＧＷＳ２を次式により求める。

ＧＷＳ２−１バフ　・ＴＳＷ−ＦＯＶｚ　’　）　　−
（６）但し；γ・・・核磁気回転比 ＴＳＷ・・・ワープ時間 このＧＷＳ２はシーケンス記憶回路４６に与えられ、勾
配磁場駆動回路４３によりＧＷＳ２の間隔で変化するワ
ープ勾配を印加する。

ステップ４ 計算機４７はステップ２とステップ３で得られた演算結
果に基づき、シーケンス記憶回路４６に視野ＦＯＶ２　
、厚み方向のマトリクスＮＺワープステップＧＷＳ　２
．スキャン回数ＮＡ′と厚み方向部分スキャンマトリク
ス数Ｎｅｔとを出力する。シーケンス記憶回路４６は与
えられた条件により第５図に示す勾配磁場駆動回路４３
．ゲート変調回路４８．ＲＦ発振回路４９．ＲＦ電力増
幅器４４を動作させてスキャンを行う。この場合、スキ
ャンの回数は与えられたＮＡ　　、　　Ｎｔ−、ＮＺ′
から厚み方向マトリクスＮ２′で部分的に異なる平均回
数をもつスキャンを行う（Ｎ　ｒ−はワープ勾配の低周
波部を中心に行う）。このようにして行ったスキャンの
方法と従来の方法との比較の一例を第６図に示す。図に
おいて、（イ）は従来の方法によるスキャンの図で、厚
み方向のマトリクスをＮｏ、１〜Ｎｏ、４のＮＺ−４の
スライスによって行い、スキャンの平均回数を２回とし
た図、（ロ）は本実施例による３次元イメージング方法
により平均回数Ｎ＾　−１，Ｎｌ、−２″ｔ−行うスキ
ャンで、いずれの場合もそのスキャンの数は（イ）では
４Ｘ２−８．　　（ロ）では６＋２−８で同一回数のス
キャンを行っている。

ステップ５ 厚み方向へのマトリクスＮＺ′点のデータを計算機４７
により離散フーリエ変換（Ｄ　Ｆ　Ｔ）処理する。

ステップ６ 第６図（ロ）のエリアジングデータＮＯ６１とＮｏ、６
を捨て、有効データＮｏ、２〜Ｎｏ、５を計算機４７の
メモリに保存する。この状態を第７図に示す。図におい
て、（イ）図は第６図（ロ）に示したスキャンにより得
た第１スラブ１０ａのデータで＃３．　＃４は部分的平
均を行ったデータである。（ロ）は（イ）のマトリクス
ＮＺ′のデータをフーリエ変換したデータで、６１はＮ
Ｚのうち両端のスライスにおける＃１．　　＃６のエリ
アジングデータとして捨てるデータ、６２は残しておく
スライスＮＺに相当する有効データである。

ステップ７ 第７図（ロ）のＮ２に相当する有効データについて計算
機４７は２次元画像を再構成する。

尚、ここでｗ４６図の例について、ステップ２に示した
式に数値を代入して検討してみる。この例では、Ｎ−＝
２．ＮＺ−４が与えられている。Ｃ−１，２として（１
）式から厚み方向マトリクスＮ２／は ＮＺ　　　＝ｌＮＴ（４Ｘ　１．２７２．０＋０．９９
９）　Ｘ　２−１３（３）式からスキャン平均回数Ｎ　
％　はスキャンを行うビュー数は ■ し　　−ＩＮＴ（１−）−１ Ｎ＋、−２Ｘ４−　１ｘ６−２ 以上の結果から第６図（ロ）のスキャンが行われること
が分る。

以上詳細に説明したように本実施例によれば、使用者側
に指定されたパラメータであるスキャン時間を延長する
ことなく、又、ＳＮ比や分解能の劣化なしにエリアジン
グイメージを分離除去し、効率良く有効イメージを得る
ことができる。

従来の選択励起によるマルチスラブ法において連続する
マルチスラブの３次元イメージングでは、厚み方向視野
の両端のスライスがエリアジングにより潰れるために一
度に関心部位の有効イメージが得られず、２回に亘るス
キャンを行ってエリアジングによる欠落イメージを埋め
ていたが、本実施例により厚み方向ＦＯＶを互いに重な
るよう大きく取ることにより、エリアジングイメージを
除去すれば、欠落イメージなしに連続した関心部位の有
効イメージが１回のスキャンで得られ、効率の良いイメ
ージングが行える。

尚、本発明は本実施例に限定されるものではない。この
自動調整法は他の３Ｄスピンワーブ法のパルスシーケン
スを用いてもよい。

又、深さ方向のフーリエ変換は離散フーリエ変換に限る
必要はなく、これと等価なアルゴリズムを用いてもよい
。

（発明の効果） 以上詳細に説明したように本発明によれば、スキャン時
間の増大、分解能の劣化、ＳＮ比の劣化等の不具合を生
ずることなくエリアジングイメージを分離除去し、効率
良く有効なイメージを得ることができて、実用上の効果
は大きい。

【図面の簡単な説明】

第１図は本発明の一実施例の方法実施の手順を示すフロ
ーチャート、 第２図は本発明の方法に用いるパルスシーケンスの図、 第３図は３次元イメージング方法において、厚み方向に
スラブを設けて厚み方向のスライスを分割する方法を示
す図、 第４図は１つのスラブをスライス面で分割した状態を示
す図、 第５図は本発明の実施例の方法を実施するための装置の
図、 第６図は本発明の実施例の方法を用いて行ったスキャン
の説明図、 第７図は第６図の方法によりスキャンして得たデータの
図、 第８図は従来のＭＲＩのパルスシーケンスの図、第９図
は従来のＭＲＩによる３次元イメージング方法のパルス
シーケンスの図、 第１０図はＲＦコイルの感度分布を利用した視野制限の
概念図、 第１１図はＲＦの選択励起を利用した視野制限の概念図
である。 １．２１・・・９０″パルス ２．１１，２２．３１・・・スライス勾配５．２５・・
・１８０ｍパルス ６．２６・・・第１ワーブ勾配 ７．２７・・・ＳＥ倍信 号、２９・・・第２ワーブ勾配 置０・・・被検体　　　　１０ａ・・・第１スラブ１０
ｃ・・・第２スラブ　１０ｅ・・・第３スラブ４１・・
・マグネットアセンブリ ４４・・・ＲＦ電力増幅器 ４６・・・シーケンス記憶回路 ４７・・・計算機　　　　５２・・・操作コンソールＮ
ａ・・・指定されたスキャン繰り返し平均回数Ｎａ　・
・・実スキャン繰り返し平均回数ＮＺ・・・指定された
厚み方向マトリクスＮＺ　・・・実スキャン厚み方向マ
トリクスＦＯｖｚ・・・厚み方向視野 ＦＯｖｚ　・・・実スキャン厚み方向視野ＮＡ　・・・
実スキャン整数平均回数

Claims (1)

1. 【特許請求の範囲】 第１ワープ勾配（６）をワープ軸に印加し、第２ワープ
   勾配（９）をスライス軸に印加してスキャンを行うＭＲ
   Ｉの３次元イメージング方法において、 使用者指定の厚み方向の視野ＦＯＶ＿Ｚ、厚み方向マト
   リスクＮ＿Ｚ及びスキャン平均回数Ｎ＿ａを入力する条
   件設定の段階と、 実スキャンにおいて厚み方向視野を増加してエリアジン
   グイメージを有効視野から追い出すために次式により実
   スキャン厚み方向視野ＦＯＶ＿Ｚ′、実スキャン厚み方
   向マトリクスＮ＿Ｚ′、スキャン平均回数Ｎ＿ａ′及び
   Ｎ＿Ａ′、厚み方向部分平均ビュー数Ｎ＿ｆ＿ｒを演算
   により求める段階と、Ｎ＿Ｚ′＝ＩＮＴ（Ｎ＿Ｚ・ｃ／
   ２．０＋０．９９９）×２ｃはｃ≧１．２ ＦＯＶ＿Ｚ′＝ＦＯＶ＿Ｚ×Ｎ＿Ｚ′／Ｎ＿ＺＮ＿ａ′
   ＝Ｎ＿ａ・Ｎ＿Ｚ／Ｎ＿Ｚ′ Ｎ＿Ａ′＝ＩＮＴ（Ｎ＿ａ′） Ｎ＿ｆ＿ｒ＝Ｎａ・Ｎ＿Ｚ−Ｎ＿Ａ′・Ｎ＿Ｚ′実スキ
   ャン第２ワープステップＧＷＳ２を次式により求めて決
   定し、そのＧＷＳ２の間隔を持つワープ勾配を印加する
   段階と、 ＧＷＳ２＝１／（γ・ＴＳＷ・ＦＯＶ＿Ｚ′）ＴＳＷ・
   ・・・ワープ時間 厚み方向マトリクスＮ＿Ｚ′、スキャン平均回数Ｎ＿Ａ
   ′及び厚み方向部分平均ビュー数Ｎ＿ｆ＿ｒに基づいて
   厚み方向マトリクスＮ＿Ｚ′で部分的に異なる実スキャ
   ン平均回数Ｎ＿Ａ′のスキャンを行う段階と、 厚み方向マトリクスＮ＿Ｚ′におけるデータの離散フー
   リエ変換演算を行ってデータ処理を行う段階と、 前記段階においてフーリエ変換を行ったデータから不必
   要なエリアジングデータを捨て有効データのみを保存す
   る段階と、 前記段階で得た有効データにより２次元画像を再構成す
   る段階とから成ることを特徴とするＭＲＩの３次元イメ
   ージング方法。