JPH09122099A - Ｍｒｉパルスシーケンスの自動補正方法、ｍｒイメージング方法、及びｍｒｉ装置 - Google Patents

Abstract

(57)【要約】 【課題】高周波反転パルスのフリップ角の変動に高い耐
性を呈し、安定した動作のプリスキャンを実施でき、本
スキャン時のパルスシーケンスをより高精度に補正す
る。 【解決手段】高周波励起パルス及び複数の高周波反転パ
ルスを有しかつその複数の高周波反転パルスの印加位相
を基準位相値に設定したパルスシーケンスの第１のプリ
スキャン、及び高周波励起パルス及び複数の高周波反転
パルスを有しかつその複数の高周波反転パルスの内の偶
数番目の高周波反転パルスの印加位相を基準位相値に対
して１８０°の位相差を有する値に設定したパルスシー
ケンスの第２のプリスキャンを各別に実施する。第１の
プリスキャンにより得た複数のエコー信号から成る第１
のエコーデータ群と第２のプリスキャンにより得た複数
のエコー信号から成る第２のエコーデータ群とに基づき
本スキャンのパルスシーケンスを事前に補正する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】この発明は、ＭＲＩ（磁気共
鳴イメージング）に係り、とくに、ＣＰＭＧパルス系列
や位相反転ＣＰパルス系列など、複数の高周波反転パル
ス（ＲＦrefocusing pulse）を組み込んだパルスシーケ
ンス（代表的なシーケンス名としては例えば、ＲＡＲＥ
法（高速ＳＥ法）、ＧＲＡＳＥ法（Hybrid ＥＰＩ
法））を使う磁気共鳴イメージングに関する。

【０００２】

【従来の技術】原子核スピンの磁気共鳴現象を利用した
磁気共鳴イメージング（ＭＲＩ）は、非侵襲で被検体内
の画像を得ることができることから、医療技術の分野で
も益々盛んに実施されており、画像処理などの技術の進
歩、高度化に伴って、ＭＲ画像に対する質的要求及び高
速撮影の度合いも非常に高まっている。

【０００３】従来、これらの条件を満足させようとして
磁気共鳴イメージングのための各種のパルスシーケンス
が実施あるいは提案されている。その一つとして、ＲＡ
ＲＥ法と呼ばれるパルスシーケンス（例えば、文献
１）；"Magnetic Resonance inMedicine Vol.3, 823-83
3, 1986" 、文献２）；"Journal of MagneticResonanc
e Imaging Vol.78, 397-407, 1988" 参照）や、ＧＲＡ
ＳＥ法と呼ばれるパルスシーケンス（例えば“米国特許
第５２７０６５４号”参照）がある。

【０００４】図３５にはＲＡＲＥ法のパルスシーケンス
を、図３６にはＧＲＡＳＥ法のパルスシーケンスをそれ
ぞれ示す。これらの図において、パルスＲＦexは高周波
励起パルス(RF excitation pulse) を、パルスＲＦre1
〜ＲＦre5 は高周波反転パルス(RF refocusing pulse)
を示す。高周波励起パルスＲＦexに続いて高周波反転パ
ルスＲＦre1 〜ＲＦre5 が順次印加される。また、パル
スＧssはスライス方向傾斜磁場を、パルスＧpeは位相エ
ンコード方向傾斜磁場を、パルスＧroはリード方向傾斜
磁場を、さらに波形Ｅ(1) 〜Ｅ(5) （又はＥ(1,1) 〜Ｅ
(5,3）はＣＰＭＧパルス系列によって発生するＮＭＲエ
コー信号（位相検波後の波形）をそれぞれ示している。
高周波パルスＲＦex、ＲＦre1 〜ＲＦre5 を印加すると
きの位相φは、高周波励起パルスＲＦexに対してφ＝０
°、高周波反転パルスＲＦre1 〜ＲＦre5 に対してφ＝
９０°に夫々設定されている。また、高周波パルスによ
って磁化スピンを倒すフリップ角θは、高周波励起パル
スＲＦexに対してθ＝α（一般に９０°）、高周波反転
パルスＲＦre1 〜ＲＦre5 に対してθ＝β１〜β５（一
般に１８０°）である。図３５、３６のいずれのパルス
シーケンスもＣＰＭＧパルス系列を利用して複数のエコ
ー信号を発生させ、それぞのエコー信号に大きさ（パル
ス面積）の異なる位相エンコード（位相エンコード方向
傾斜磁場パルスＧpe）を施すことで、１フレームの画像
の再構成に必要なエコーデータを収集している。これに
より、従来のスピンエコー法に比べて、撮像時間が数分
の１から数百分の１で済む。

【０００５】このＣＰＭＧパルス系列、位相反転ＣＰパ
ルス系列の如く、複数の高周波反転パルスを用いたパル
ス系列の場合、その各高周波反転パルス間で収集される
それぞれのエコー信号は位相ダイヤグラム上の複数の経
路を経たエコー成分の総和であることが、明らかになっ
ている（例えば、文献２）；"Journal of MagneticReso
nance Imaging Vol.78, 397 -407, 1988" 、文献
３）；"MagneticResonance in Medicine Vol.30, 183-1
91, 1993"参照）。

【０００６】また、これらのパルスシーケンスを使う場
合、傾斜磁場パルスのスイッチングに起因してＭＲＩ装
置の導体部分に生じる渦電流、傾斜磁場コイルの製造誤
差、そのほかのＭＲＩ装置の校正の不完全性を主たる原
因として、各エコー成分に、互いに異なる空間的分布を
持った位相ずれが生じる。このことは、具体的には例え
ば、特許公報１）；特開平６−１２１７７７号、特許公
報２）；特開平６−５４８２７号、さらには、特許公報
３）；米国特許第５３７８９８５号に記載されている。
このため、再構成画像の位相エンコード方向にゴースト
アーチファクトが発生したり、画像の信号値が局所的に
低下したり、画像全体のＳ／Ｎ比が低下するなどの現象
が生じ、画質を著しく劣化させてしまうことがある。こ
の不具合は、傾斜磁場パルスを短いスイッチング時間の
間に強度を大きく変化させなけらばならない高空間分解
能の撮像や、エコー数を非常に多くとる撮像の場合に特
に顕著である。

【０００７】この画質劣化の不具合を解消する一つの方
法（第１の従来技術）が、前述した特許公報３）；米国
特許第５３７８９８５号に示されている。被検体に対す
るＭＲ画像取得のための実際の撮像（以下、「本スキャ
ン」という）に先立って、プリスキャンと呼ばれるスキ
ャンを実施する。このプリスキャンでは、位相エンコー
ド方向の傾斜磁場パルスの強度を零にした状態でスキャ
ンを実施し、これにより得た各エコー信号に１次元のフ
ーリエ変換を施す。これにより得た位相分布の０次、１
次の成分を補正データとして算出し、本スキャンの傾斜
磁場パルスの波形および高周波パルスを印加する位相
（以下、「印加位相」という）を補正している。この印
加位相は、図３８に示す如く、回転座標系Ｘ′，Ｙ′，
Ｚ′のＸ′，Ｙ′面における基準とする座標軸、例えば
Ｘ′に対する高周波パルスの印加時の磁化Ｍの位相φと
して表される。

【０００８】ところで、このＣＰＭＧパルス系列、位相
反転ＣＰパルス系列などの複数の高周波反転パルスを用
いた撮像においては、また別の物理的性質も分かってい
る。つまり、高周波励起パルスと第１番目の高周波反転
パルスとの間で横磁化が位相誤差を受けたり、またそれ
以降の隣接する高周波反転パルス間で横磁化が一律に同
じ大きさの位相誤差を受けると、各エコー信号（複数の
経路を経たエコー成分の総和）は２つのグループのエコ
ー群に分かれる、というものである。この性質は例え
ば、文献４）；"Magnetic Resonance in Medicine Vol.
30, 251-255, 1993"に示されている。

【０００９】ＲＡＲＥ法やＧＲＡＳＥ法などのパルスシ
ーケンスを使った場合、前述したように、ほぼ同じ形の
傾斜磁場パルスが繰り返して印加されるので、傾斜磁場
パルス自体は勿論のこと、傾斜磁場パルスの大きさに比
例する性質を有する渦磁場もそのパルス繰返しパターン
に基づく位相誤差を生じて、２つのエコー成分のグルー
プに分かれると考えられる。したがって、各エコー信号
の位相ずれのパターンと再構成画像の画質劣化は、近似
的には、２つのグループのエコー成分間の相互干渉に起
因していると、考えることができる。

【００１０】高周波励起パルスと第１番目の高周波反転
パルスとの間で原子核スピンの横磁化が受ける位相誤差
をΔφ１、それ以降の隣接する高周波反転パルス間でそ
の横磁化が一律に受ける位相誤差をΔφ２とし、第ｎ番
目と第ｎ＋１番目の高周波反転パルス間で観測される２
つのグループのエコー成分の内、全ての高周波反転パル
スで位相反転(refocus) される経路のエコー成分を含む
ものを主エコー成分Ｅmain(n) 、それ以外のエコー成分
を副エコー成分Ｅsub(n)とすると、両エコー成分Ｅmain
(n) 、Ｅsub(n)の位相は、

【数１】

【数２】 で表される。

【００１１】撮像領域の各々の点において生じるエコー
信号は、このような対応関係があるため、それらの総和
である（観測される）エコー信号の位相や位相のちらば
り程度を示す位相分散の大きさについても同様の対応関
係がある。図３８は各高周波反転パルス間でのリード方
向傾斜磁場に因る各エコー信号の位相分散の大きさの変
化を模式的に説明する位相ダイヤグラムである。同図の
縦軸はリード方向傾斜磁場に因るエコー信号の位相分散
の大きさ（相対値）であり、位相分散の等しいエコー成
分の同図上での変化を示す経路が分かり易くなるよう
に、故意に、高周波励起パルスと第１高周波反転パルス
との間のリード方向傾斜磁場の時間積分値を正規の値Ａ
（正規には同図中のパルスＢの１／２）よりもΔＡだけ
増やしている。前述した文献２）に示す如く、ある位相
分散を有する横磁化成分（エコー成分）は高周波反転パ
ルスによって、分散関係が反転する横磁化成分と、反転
しない横磁化成分と、縦磁化として保存される成分の３
つに分化することが知られている。その後の高周波反転
パルスの印加間隔、印加する傾斜磁場波形は規則的であ
ることから、位相ダイヤグラム上のいくつかの経路のエ
コー成分は等しい位相分散の大きさを持って重なるた
め、同図に示すように規則的な分散パターンになる。同
図中、太い実線は主エコー成分Ｅmainと同位相のエコー
成分のグループで、細い実線は副エコー成分Ｅsub と同
位相のエコー成分のグループである。またエコー収集の
期間において、位相分散の大きさが正又は負の方向に大
きく離れているエコー成分の殆どは、そのエコー成分が
属するグループ内のスピン同士の位相干渉に因ってエコ
ー信号としては観測できない、いわゆるスポイルされた
状態にある。したがって、実際にエコー信号として収集
される部分は主に、図中、矢印で示す部分（データ収集
中に位相分散が零になる瞬間がある成分）ということに
なる。

【００１２】一方、エコー信号の位相を見てみると、エ
コー信号の位相と高周波パルスの印加位相の対応関係か
ら、ＣＰＭＧパルス系列又は位相反転ＣＰパルス系列の
印加位相の条件を満たし、そのほかハードウエアの不完
全性に因るスピン位相の誤差が無いものとすると、主エ
コー成分Ｅmain(n) と副エコー成分Ｅsub(n)の位相は一
致する。実際に収集されるエコー信号は、両エコー成分
の和であり、

【数３】Ｅ(n) ＝Ｅmain(n) ＋Ｅsub(n) となる。未調整でも図３８中のΔＡは十分に小さいた
め、２つのエコー成分のピークは重なって見えることが
多い。

【００１３】また、両方のエコー成分Ｅmain(n) 、Ｅsu
b(n)の振幅は、高周波反転パルス数ｎ、高周波パルスの
スライス特性、フリップ角などに依存することが確認さ
れている（例えば、文献５）；"Society of Magnetic R
esonance in Medicine, Pro.of Annual Meeting in 199
2 No.4508" ，文献６）；" 同 in 1991 No.1025 "参
照）。

【００１４】このため、両方のエコー成分Ｅmain(n) 、
Ｅsub(n)の振幅比Ｒ(n) が分からない限り、収集される
エコー信号から位相誤差を抽出することはできない。

【００１５】このような状況下において、前述した特許
公報２）；特開平６−５４８２７号記載の磁気共鳴画像
化方法および装置には、上述した不具合を解消すべく別
の従来技術（第２の従来技術）が開示されている。この
開示による補正処理を上述と同一の記号で示すと、

【数４】｜Ｅsub(1)｜＝０ であり、第１番目の高周波反転パルスのフリップ角がほ
ぼ１８０°である場合、

【数５】｜Ｅsub(2)｜＝ほぼ零 である。このため、

【数６】 となる。これを利用して、前述と同様のプリスキャンか
ら得られた第１エコー信号および第２エコー信号の情報
のみから、本スキャンの傾斜磁場波形および高周波パル
スの印加位相を修正し、比較的簡便な補正を行おうとす
るものである。

【００１６】この第２の従来技術を図３９〜４１に基づ
いて詳述する。図３９にはプリスキャンのパルスシーケ
ンスを、図４０には本スキャンのパルスシーケンスを、
図４１には補正処理のフローチャートをそれぞれ示す。

【００１７】図４１に示すように、最初に図３９のパル
スシーケンスにしたがってプリスキャンが実施される
（ステップ１０１）。このプリスキャンは位相エンコー
ド方向の傾斜磁場が常に零である（Ｇpe＝０）。いま何
らかの原因に拠って位相ずれが在る場合、同図に示すよ
うに第１エコー信号Ｅ(1) を除く各エコー信号はそれぞ
れ２つのエコー成分Ｅmain(n) 、Ｅsub(n)に分かれる。
実際には、２つのエコー成分Ｅmain(n) 、Ｅsub(n)のず
れは小さく、エコーピークは２つに分かれず、エコーピ
ークの位置ずれｐ(n) と位相のずれφ(n) が各エコー毎
に交互に振動するように変化する。第２の従来技術によ
れば、第１、第２エコー信号Ｅ(1) 、Ｅ(2) のエコーピ
ークの位置ずれｐ(1) ，ｐ(2) および位相ずれφ(1) ，
φ(2) が測定される（ステップ１０２）。磁場中心から
選択されたスライス面までの距離をｘｓとすると、高周
波パルスの印加位相の補正量Δφ（あるいはスライス方
向傾斜磁場パルスＧssの補正量ΔＧssに対する「ΔＧss
・ｘｓ」）は位相ずれφ(1)，φ(2) の位相差に比例
し、ΔＧroは位置ずれｐ(1) ，ｐ(2) の差に比例する。
この関係に基づいて、求めたｐ(1) ，ｐ(2) ，φ(1) ，
φ(2) から本スキャン時の高周波パルスの印加位相の補
正量Δφ及びスライス方向傾斜磁場パルスＧssの補正量
ΔＧssのいずれか、及びリード方向傾斜磁場パルスＧro
の補正量ΔＧroが演算される（ステップ１０３）。つい
で、これらの補正量Δφ（又はΔＧss）およびΔＧroで
パルスシーケンスが修正され（図４０参照）、このパル
スシーケンスにより本スキャンが実行される（ステップ
１０４）。

【００１８】

【発明が解決しようとする課題】上述した第１、第２の
従来技術は、それぞれのエコー成分が相異なる空間的分
布を持った位相ずれを起こしたことに因る位相エンコー
ド方向のアーチファクトなどの不具合を改善しようとす
るものではあったが、以下の示す如く、未解決の課題が
あった。

【００１９】（１） 第１に、高周波反転パルスのフリ
ップ角が１８０°からずれている場合、本スキャンに対
する高周波パルスの印加位相及び傾斜磁場波形の補正精
度が低下するという問題である。

【００２０】例えば｜Ｅmain(n) ｜＝｜Ｅsub(n)｜の場
合、２つのエコー成分Ｅmain(n) 、Ｅsub(n)の位相差が
いくら大きくても、その２つのエコー成分Ｅmain(n) 、
Ｅsub(n)のベクトル和である第ｎ番目に収集されるエコ
ー信号の位相は変化しない。この理由に拠り、第１、第
２の従来技術のいずれの手法を実施したとしても補正が
的確に実施されず、ＭＲ画像の画質は一向に改善されな
いか、改善されたとしても不十分な程度に終わることが
多い。

【００２１】フリップ角が１８０°からずれる又はずら
しているという状況は決して特殊なものではない。例え
ば、前述した文献５）；"Society of Magnetic Resonan
cein Medicine, Pro. of Annual Meeting in 1992 No.4
508" でも述べられているように、現状のＦＳＥ（Fast
ＳＥ）法の場合でも、第４〜第８番目のエコーで主エコ
ー成分Ｅmain(n) と副エコー成分Ｅsub(n)の振幅がほぼ
等しくなっていることが報告されている。また、文献
３）；"Magnetic Resonance in Medicine Vol.30, 183-
191, 1993"では、各エコー信号の振幅を安定化させるた
めに高周波反転パルスのそれぞれのスライス特性、フリ
ップ角を個別に変化させている。さらに、この技術を利
用してＴ２緩和による各エコー信号の減衰を減らす試み
も報告されている（"Society of Magnetic Resonance,
Proc. of 2nd Annual Meeting in1994 No.27"参照）。
さらに、マルチスライス撮像の場合、高周波反転パルス
のフリップ角を１８０°以下の値に設定し、ＳＡＲ（Ｒ
Ｆ被爆）の低減、ＭＴＣ効果による組織コントラストの
変化を低減させることができるという報告もある（"Soc
iety of Magnetic Resonance in Medicine, Proc. of A
nnual Meeting in1993 No.1244"参照）。

【００２２】以上のように、高周波反転パルスのフリッ
プ角が１８０°からずれている場合も多く、かかる状況
下で前述した第１、第２の従来技術による補正はその精
度及び安定性の面で不足である。

【００２３】この問題は前述した特許公報２）；特開平
６−５４８２７号でも指摘されているところである（同
公報第１１ページ、１８−２７行目参照）。この問題に
対して、同特許公報の技術では、スライス方向の傾斜磁
場を駆使して一方のエコー成分をディフェーズさせるな
どの回避方法が示されているが、傾斜磁場パルスの波形
を変形させるので、システムの校正の不完全性の影響を
受けるなど、シーケンス調整の自動化には適さないこと
が多いという問題がある。

【００２４】（２） 第２の問題は、複数ある位相ずれ
の要因のうち、従来の方法で補正できるのは一部の要因
によるものであり、全てではないということである。前
述した第１、第２の従来技術のいずれにあっても、プリ
スキャンの位相エンコード方向の傾斜磁場パルスＧpeは
常に零である。このため、プリスキャンにより得られる
エコー信号は本スキャンにおける位相エンコード方向の
傾斜磁場パルスに起因する位相ずれ成分を含んでいな
い。つまり、ａ）；位相エンコード方向の傾斜磁場パル
ス自体のずれ、ｂ）；位相エンコード方向の傾斜磁場パ
ルスに起因する渦磁場に因る０次あるいは１次以上の空
間分布を持つ位相ずれ、などは原理的に補正できない。
スキャンの空間分解能の高度化、高速化が進んでいる昨
今において、パルスシーケンスの種類によっては、位相
エンコード方向の傾斜磁場に起因する位相ずれ成分を到
底無視できない状況にあり、かかる補正に対する制限は
何としても打破したいところである。

【００２５】（３） 第３の問題は、補正値の精度が撮
像対象の状態により変化してしまうことである。前述し
た第１、第２の従来技術では、撮像対象からのエコー信
号をもとに補正情報を求めているため、撮像対象の形
状、状態によっては正確な位相誤差を測定できないこと
もあり、補正の安定性の面で問題がある。

【００２６】例えば医用診断装置の場合、けい椎のアキ
シャル像のＴ２強調画像では主な信号源が断面に対して
ほぼ垂直な方向に流れる脳脊髄液(Cerebral Spinal Flu
id:CSF) であるため、パルスシーケンスで使用している
傾斜磁場パルスに因り流れ速度に関係した原子核スピン
の位相変化が生じる。この位相変化を従来の第１、第２
の従来技術に示された手法で補正する場合、かかる位相
変化はほぼそのままプリスキャンでも表れるから、プリ
スキャンから求められる補正量の正確性に欠け、信頼性
の低いものとなって、目的とする位相変化を誤って補正
してしまうこともあるなど、補正精度が低く、その安定
性に劣る。このため、原子核スピンが静止している部位
でも、ゴーストアーチファクトが現れたり、信号値が低
下してＳ／Ｎ比が低下するなどの不具合に帰着するし、
また、高いエコー信号が出る部位の動きが激しい場合、
画質が劣化する場合がある。これは、本スキャン時のパ
ルスシーケンスの補正情報（補正量）を形状、状態が未
知の撮像対象から予め得ることを特徴とする補正手法の
場合、ある程度避けられない問題ではあるものの、臨床
診断上問題となっている。

【００２７】この発明は上述した不具合、問題に鑑みて
なされたものである。（１）；プリスキャンを実施して
本スキャン時のパルスシーケンスに必要な補正量を求め
る処理を前提とする磁気共鳴イメージングにおいて、従
来よりも高周波反転パルスのフリップ角の変動に高い耐
性を呈し、安定した動作のプリスキャンを実施でき、本
スキャン時のパルスシーケンスをより高精度に補正で
き、再構成されるＭＲ画像の画質劣化を抑えることがで
きるＭＲＩパルスシーケンスの自動補正方法及びＭＲＩ
装置を提供することを、本発明の第１の目的とする。ま
た、（２）；プリスキャンを実施して本スキャン時のパ
ルスシーケンスに必要な補正量を求める処理を前提とす
る磁気共鳴イメージングにおいて、位相エンコード方向
の傾斜磁場パルスに起因する位相ずれ成分も精度良く補
正できるようにし、再構成されるＭＲ画像の画質向上を
図ることを、本発明の第２の目的とする。また、
（３）；プリスキャンを実施して本スキャン時のパルス
シーケンスに必要な補正量を求める処理を前提とする磁
気共鳴イメージングにおいて、主要な信号源が高速で動
いている場合でも、本スキャン時のパルスシーケンスの
補正精度の低下を防止し、再構成されるＭＲ画像の画質
劣化を抑えることを、本発明の第３の目的とする。さら
に、（４）；プリスキャンからの補正データに拠る本ス
キャンのパルスシーケンスの補正を行わずとも、結果と
して、上記第１〜第３の目的を同時に達成することでき
るＭＲイメージング法及びＭＲＩ装置を提供すること
を、本発明の第４の目的とする。

【００２８】以上を集約すると、上記第１〜第４の目的
の少なくとも１つを達成させて、ＣＰＭＧパルス系列、
位相反転ＣＰパルス系列などの複数の高周波反転パルス
を使うパルスシーケンスを用いた撮像において、必要な
補正情報の収集と補正を従来よりも系統的に実施して、
かかる補正の精度及び安定性を向上させ、又は、位相誤
差分布の影響をそれらの補正以外のアプローチによって
回避し、空間分解能、ゴーストアーチファクト、画像コ
ントラストなどの面で高画質のＭＲ画像を提供すること
である。

【００２９】

【課題を解決するための手段】上記目的を達成させるた
め、請求項１〜９，１７記載の発明によれば、被検体の
ＭＲ画像を得るための本スキャンの前に、高周波励起パ
ルス及び複数の高周波反転パルスを有しかつその複数の
高周波反転パルスの印加位相を基準位相値に設定したパ
ルスシーケンスの第１のプリスキャン、及び高周波励起
パルス及び複数の高周波反転パルスを有しかつその複数
の高周波反転パルスの内の偶数番目の高周波反転パルス
の印加位相を前記基準位相値に対して１８０°の位相差
を有する値に設定したパルスシーケンスの第２のプリス
キャンを各別に実施し、前記第１のプリスキャンにより
得た複数のエコー信号から成る第１のエコーデータ群と
前記第２のプリスキャンにより得た複数のエコー信号か
ら成る第２のエコーデータ群とに基づいて前記本スキャ
ンのパルスシーケンスを事前に補正する、ことを特徴と
したＭＲＩパルスシーケンスの自動補正方法である。

【００３０】例えば、前記第１、第２のプリスキャンの
パルスシーケンス及び前記本スキャンのパルスシーケン
スはＣＰＭＧパルス系列又は位相反転ＣＰパルス系列で
ある。また例えば、前記第１、第２のプリスキャンのパ
ルスシーケンス及び前記本スキャンのパルスシーケンス
はＧＲＡＳＥ法又は高速ＳＥ法に従うパルス系列であ
る。また例えば、前記第１、第２のプリスキャンのパル
スシーケンスの位相エンコード用傾斜磁場は常に零であ
る。また前記自動補正の処理は、典型的には、前記第１
のエコーデータ群と前記第２のエコーデータ群との対応
する順番のエコーデータ同士で加算又は減算して主エコ
ー成分から成る主エコー集合体とその主エコー成分以外
のエコー成分から成る副エコー集合体とを分離・抽出す
る処理と、前記主エコー集合体及び副エコー集合体の各
集合体の内の奇数番目及び偶数番目のエコーデータのエ
コーピークの位相ずれ及び位置ずれの少なくとも一方を
演算する処理と、この演算値から求めた補正データによ
り前記本スキャンのパルスシーケンスを自動的に補正す
る処理と、を含む。

【００３１】また請求項１０、１１記載の発明によれ
ば、前記第１、第２のプリスキャンのパルスシーケンス
の位相エンコード用傾斜磁場は、前記本スキャンにおけ
るｋ空間の中心位置を含むデータラインに配置するエコ
ーデータ収集時のショットと同一の位相エンコード用傾
斜磁場の波形であることを特徴とする。

【００３２】また請求項１２、１３記載の発明によれ
ば、前記第１、第２のプリスキャンのパルスシーケンス
の位相エンコード用傾斜磁場は、前記本スキャンにおけ
るｋ空間上の一部又は全面のエコーデータ収集時のショ
ットと同一の位相エンコード用傾斜磁場の波形であるこ
とを特徴とする。

【００３３】さらに請求項１４記載の発明によれば、前
記第１、第２のプリスキャンを実施する対象は、前記被
検体とは異なるファントムである。

【００３４】さらに請求項１５、１６記載の発明によれ
ば、その主な特徴は、前記第１、第２のエコーデータ群
とに基づく補正処理は、この第１、第２のエコーデータ
群に基づいてイメージングチャンネルと傾斜磁場コイル
の物理チャンネルとの組み合わせに対して固有の位相誤
差に関連する情報を求める処理と、本スキャン時にはそ
の位相誤差関連の情報と傾斜磁場に関連した撮像条件と
から本スキャンのパルスシーケンスの補正データを求め
る処理とを含むことである。

【００３５】本発明の別の側面として、請求項１８記載
の発明によれば、高周波励起パルス及び複数の高周波反
転パルスを有しかつその複数の高周波反転パルスの印加
位相を基準位相値に設定したパルスシーケンスの第１の
スキャン、及び高周波励起パルス及び複数の高周波反転
パルスを有しかつその複数の高周波反転パルスの内の偶
数番目の高周波反転パルスの印加位相を前記基準位相値
に対して１８０°の位相差を有する値に設定したパルス
シーケンスの第２のスキャンを対象物に各別に実施する
ステップと、前記第１のスキャンにより得た複数のエコ
ー信号から成る第１のエコーデータ群と前記第２のスキ
ャンにより得た複数のエコー信号から成る第２のエコー
データ群とを加算又は減算して主エコー成分の主エコー
集合体とその主エコー成分以外のエコー成分から成る副
エコー集合体とを分離・抽出するステップと、を含むＭ
Ｒイメージング方法が提供される。

【００３６】本発明のさらに別の側面として、請求項１
９〜２７記載の発明によれば、高周波励起パルス及び複
数の高周波反転パルスを有し、その複数の高周波反転パ
ルスの印加位相を基準位相値に設定し、かつｋ空間デー
タ収集用の位相エンコードを施した撮像用パルスシーケ
ンスの第１のスキャンと、高周波励起パルス及び複数の
高周波反転パルスを有し、その複数の高周波反転パルス
の内の偶数番目の高周波反転パルスの印加位相を前記基
準位相値に対して１８０°の位相差を有する値に設定
し、かつｋ空間データ収集用の位相エンコードを施した
撮像用パルスシーケンスの第２のスキャンとを対象物に
各別に実施するステップと、前記第１のスキャンにより
得た複数のエコー信号から成る第１のエコーデータ群と
前記第２のスキャンにより得た複数のエコー信号から成
る第２のエコーデータ群とを加減算して主エコー成分の
集合体からなるｋ空間データとその主エコー成分以外の
エコー成分の集合体からなるｋ空間データとを分離・抽
出するステップと、前記２組のｋ空間データ同士で、ｋ
空間上の複数のエコー配置区画の内の偶数番目又は奇数
番目のエコー配置区画のエコーデータを互いに交換する
ステップと、前記エコーデータ交換後の２組のｋ空間デ
ータを各別に再構成して２つの実空間画像を作るステッ
プと、前記２つの実空間画像から１つの実空間画像を形
成するステップと、を含むＭＲイメージング方法が提供
される。

【００３７】とくに、請求項２２記載の発明に係るＭＲ
イメージング方法では、前記対象物へのプリスキャンと
して、高周波励起パルス及び複数の高周波反転パルスを
有し、その複数の高周波反転パルスの印加位相を基準位
相値に設定したパルスシーケンスの第１のプリスキャン
と、高周波励起パルス及び複数の高周波反転パルスを有
し、その複数の前記高周波反転パルスの内の偶数番目の
高周波反転パルスの印加位相を前記基準位相値に対して
１８０°の位相差を有する値に設定たパルスシーケンス
の第２のプリスキャンとを対象物に各別に実施するステ
ップと、前記第１のスキャンにより得た複数のエコー信
号から成る第１のエコーデータ群と前記第２のスキャン
により得た複数のエコー信号から成る第２のエコーデー
タ群とを加減算して主エコー成分の集合体の振幅値情報
とその主エコー成分以外のエコー成分の集合体の振幅値
情報とを分離・抽出するステップと、前記エコーデータ
の交換ステップと前記２つの実空間画像の再構成ステッ
プとの間で実施するステップであって、前記２組の集合
体の振幅値情報を使って前記エコーデータ交換後の２組
のｋ空間データの振幅値を補正するステップと、をさら
に含むことを特徴とする。

【００３８】さらに本発明の別の側面として、請求項２
８記載の発明によって、被検体のＭＲ画像を得るための
本スキャンの前に、高周波励起パルス及び複数の高周波
反転パルスを有しかつその複数の高周波反転パルスの印
加位相を基準位相値に設定したパルスシーケンスの第１
のプリスキャン、及び高周波励起パルス及び複数の高周
波反転パルスを有しかつその複数の高周波反転パルスの
内の偶数番目の高周波反転パルスの印加位相を前記基準
位相値に対して１８０°の位相差を有する値に設定した
パルスシーケンスの第２のプリスキャンを各別に実施す
る手段と、前記第１のプリスキャンにより得た複数のエ
コー信号から成る第１のエコーデータ群と前記第２のプ
リスキャンにより得た複数のエコー信号から成る第２の
エコーデータ群とに基づいて前記本スキャンのパルスシ
ーケンスを事前に補正する手段とを備えたＭＲＩ装置が
提供される。

【００３９】さらに本発明の別の側面として、請求項２
９記載の発明により、高周波励起パルス及び複数の高周
波反転パルスを有し、その複数の高周波反転パルスの印
加位相を基準位相値に設定し、かつｋ空間データ収集用
の位相エンコードを施した撮像用パルスシーケンスの第
１のスキャンと、高周波励起パルス及び複数の高周波反
転パルスを有し、その複数の高周波反転パルスの内の偶
数番目の高周波反転パルスの印加位相を前記基準位相値
に対して１８０°の位相差を有する値に設定し、かつｋ
空間データ収集用の位相エンコードを施した撮像用パル
スシーケンスの第２のスキャンとを対象物に各別に実施
する手段と、前記第１のスキャンにより得た複数のエコ
ー信号から成る第１のエコーデータ群と前記第２のスキ
ャンにより得た複数のエコー信号から成る第２のエコー
データ群とを加減算して主エコー成分の集合体からなる
ｋ空間データとその主エコー成分以外のエコー成分の集
合体からなるｋ空間データとを分離・抽出する手段と、
前記２組のｋ空間データ同士で、ｋ空間上の複数のエコ
ー配置区画の内の偶数番目又は奇数番目のエコー配置区
画のエコーデータを互いに交換する手段と、前記エコー
データ交換後の２組のｋ空間データを各別に再構成して
２つの実空間画像を作る手段と、前記２つの実空間画像
から１つの実空間画像を形成する手段と、を備えたＭＲ
Ｉ装置が提供される。

【００４０】

【発明の実施の形態】以下、この発明の各種の実施の形
態を図面を参照して説明する。最初に、以下に詳述する
実施形態と本発明に係る上記第１〜第４の目的との対応
関係を示しておく。

【００４１】第１の実施形態は第１の目的を、第２の実
施形態及びその変形形態は第２の目的を、第３の実施形
態は同じく第２の目的を、第４の実施形態及びその変形
形態は第３の目的を、第５の実施形態は第４の目的を、
それぞれ達成しようとするものである。

【００４２】（第１の実施形態）この発明に係る第１の
実施形態を図１〜図８に基づき説明する。

【００４３】この実施形態に係る磁気共鳴イメージング
（ＭＲＩ）装置の概略構成を図１に示す。このＭＲＩ装
置は、静磁場発生用の磁石部と、静磁場に位置情報を付
加するための傾斜磁場部と、選択励起用及びＭＲ信号受
信用の送受信部と、システムコントロール及び画像再構
成を担う制御・演算部とを備えている。

【００４４】磁石部は、例えば超電導方式の磁石１と、
この磁石１に電流を供給する静磁場電源２とを備え、撮
影対象としての被検体Ｐが挿入される円筒状の診断空間
のＺ軸方向（長手方向）に静磁場Ｈ₀ を発生させる。

【００４５】傾斜磁場部は、磁石１に組み込まれたＸ，
Ｙ，Ｚ軸方向の３組の傾斜磁場コイル３ｘ〜３ｚと、こ
の傾斜磁場コイル３ｘ〜３ｚに電流を供給する傾斜磁場
電源４と、この電源４を制御する傾斜磁場シーケンサ５
ａとを備える。この傾斜磁場シーケンサ５ａはコンピュ
ータを備え、装置全体のコントローラ６（コンピュータ
を搭載）からプリスキャン及び本スキャンを実施するた
めのＦＳＥ法，ＧＲＡＳＥ法などの収集シーケンスを指
令する信号を受けると、この指令に応答して図２に示す
処理を実行する。これにより、傾斜磁場シーケンサ５ａ
は、指令されたシーケンスにしたがってＸ，Ｙ，Ｚ軸方
向の各傾斜磁場の印加及びその強度を制御し、それらの
傾斜磁場が静磁場Ｈ₀ に重畳可能になっている。この実
施例では、互いに直交する３軸の内のＺ軸方向の傾斜磁
場をスライス用傾斜磁場Ｇ_S とし、Ｘ軸方向のそれを読
出し用傾斜磁場Ｇ_R とし、さらにＹ軸方向のそれを位相
エンコード用傾斜磁場Ｇ_E とする。

【００４６】送受信部は、磁石１内の撮影空間にて被検
体Ｐの近傍に配設される高周波コイル７と、この高周波
コイル７に接続された送信機８Ｔ及び受信機８Ｒと、こ
の送信機８Ｔ及び受信機８Ｒの動作タイミングを制御す
るＲＦシーケンサ５ｂ（コンピュータを搭載）とを備え
る。ＲＦシーケンサ５ｂは前記傾斜磁場シーケンサ５ａ
と共にシーケンサ５を成す。ＲＦシーケンサ５ｂは傾斜
磁場シーケンサ５ａと同期した状態でＲＦパルスの印加
を指令できる。送信機８Ｔ及び受信機８Ｒは、ＲＦシー
ケンサ５ｂの制御のもと、核磁気共鳴（ＮＭＲ）を励起
させるためのラーモア周波数のＲＦ電流パルスを高周波
（ＲＦ）コイル７に供給する一方、高周波コイル７が受
信したＭＲ信号（高周波信号）に各種の信号処理を施し
てエコー信号を形成するようになっている。

【００４７】さらに、制御・演算部は、上述したコント
ローラ６のほか、受信機８Ｒで形成されたデジタル量の
エコー信号を受けるマルチプレクサ１１を備える。マル
チプレクサ１１はコントローラ６からの制御信号に応答
して、その出力経路をコントローラ側及び再構成ユニッ
ト側との間で択一的に切り換える。さらに、マルチプレ
クサ１１の一方の出力側には、フーリエ変換法により画
像再構成を行う再構成ユニット１２と、再構成した画像
データを保管する記憶ユニット１３と、画像を表示する
表示器１４と、入力器１５とを備えている。コントロー
ラ６は前述したようにコンピュータを内蔵し、システム
全体の動作内容及び動作タイミングを制御する。

【００４８】次に、この実施形態の動作を説明する。な
お、ここでは、パルスシーケンスとしてＦＳＥ（高速Ｓ
Ｅ）法が採用されているとする。

【００４９】コントローラ６は図２の処理を実施する。
この処理には２種類のプリスキャンＡ，Ｂが含まれてい
る。最初に、コントローラ６はマルチプレクサ１１のス
イッチ経路をコントローラ６側に切り換えた状態で、図
３に示すパルスシーケンスのプリスキャンＡをシーケン
サ５に実行させ、その結果得られたエコー信号を入力す
る（図２ステップＳ１）。

【００５０】このプリスキャンＡでは、最初に、スライ
ス用傾斜磁場パルスＧssとともに高周波励起パルスＲＦ
ex（印加位相φ＝０°、フリップ角θ＝α（ここでは９
０°））が印加される。このτ／２の時間経過後に、ス
ライス用傾斜磁場パルスＧssとともに第１番目の高周波
反転パルスＲＦre1 （印加位相φ＝９０°、フリップ角
θ＝β１（ここでは１８０°））が印加される。そし
て、最初の高周波励起パルスＲＦexの印加からτ時間後
に、リード用傾斜磁場パルスＧroを印加しながらエコー
信号Ｅ(1) 読み出される。以後、第１番目の高周波反転
パルスＲＦre1 の印加からτ時間毎に第２番目以降の高
周波反転パルスＲＦre2 ，ＲＦre3 …がスライス用傾斜
磁場パルスＧssとともに印加され、エコー信号Ｅ(2) ，
Ｅ(3) ，…が同様に読み出される。このプリンスキャン
Ａでは同図に示すように、位相エンコード用傾斜磁場パ
ルスＧpeは常に零となっている。

【００５１】このプリンスピンＡにおいては、何等かの
原因に因り原子核スピンに位相ずれがある場合、前述と
同様に、最初のエコー信号Ｅ(1) 以外のエコー信号Ｅ
(2) ，Ｅ(3) ，…のそれぞれは、２つのエコー成分、す
なわち主エコー成分Ｅmain(2)（Ｅmain(3) ，…）と副
エコー成分Ｅsub(2)（Ｅsub(3)，…）とに分かれる（図
３参照）。そのような位相ずれの原因が無いとき、主エ
コー成分Ｅmain(n) の位相arg{Ｅmain(n)}と副エコー成
分Ｅsub(n)の位相arg{Ｅsub(n)} は等しい（ｎ＝２、
３、…）。

【００５２】次いでコントローラ６は図４に示す２番目
のプリスピンＢをシーケンサ５に実行させ、その結果得
られたエコー信号を入力する（図２ステップＳ２）。

【００５３】このプリスキャンＢでも最初に、スライス
用傾斜磁場パルスＧssとともに高周波励起パルスＲＦex
（印加位相φ＝０°、フリップ角θ＝α（ここでは９０
°））が印加される。このτ／２の時間経過後に、スラ
イス用傾斜磁場パルスＧssとともに高周波反転パルスＲ
Ｆre1 （印加位相φ＝９０°、フリップ角θ＝α（ここ
では１８０°））が印加される。そして、最初の高周波
励起パルスＲＦexの印加からτ時間後に、リード用傾斜
磁場パルスＧroを印加しながらエコー信号Ｅ(1) が読み
出される。以後、第１番目の高周波反転パルスＲＦre1
の印加からτ時間毎に第２番目以降の高周波反転パルス
ＲＦre2 ，ＲＦre3 ，…がスライス用傾斜磁場パルスＧ
ssとともに印加され、エコー信号Ｅ(2) ，Ｅ(3) ，…が
同様に読み出される。

【００５４】このプリスキャンＢでは、偶数番目の高周
波反転パルスＲＦre2 ，ＲＦre4 ，…の印加位相φは図
４に記載の如く、奇数番目の高周波反転パルスＲＦre3
，ＲＥre5 ，…のそれをさらに１８０°だけ回転させ
た値（φ＝２７０°）に設定してある。また、位相エン
コード用傾斜磁場パルスＧpeはプリスキャンＡと同様に
常に零となっている。

【００５５】何等かの原因に因り原子核スピンに位相ず
れがある場合、前述と同様に、最初のエコー信号Ｅ(1)
以外のエコー信号Ｅ(2) ，Ｅ(3) ，…のそれぞれは、２
つのエコー成分、すなわち主エコー成分Ｅmain(2) （Ｅ
main(3) ，…）と副エコー成分Ｅsub(2)（Ｅsub(3)，
…）とに分かれ（図３参照）、しかも偶数番目の高周波
反転パルスＲＦre2 ，ＲＦre4 ，…の印加位相φを１８
０°だけ余計に回転させているため、副エコー成分Ｅsu
b(2)，Ｅsub(3)，…の位相arg{Ｅsub(n)} だけがプリン
スキャンＡ時の対応する位相に対して１８０°回転す
る。

【００５６】この理由を図８を用いて説明する。同図に
示すように、いま、第ｎ番目のエコー信号の位相を
θ_n ，第ｎ＋１番目、第ｎ＋２番目の高周波反転パルス
の位相をそれぞれφ_n+1 、φ_n+2 とする。第ｎ＋２番目
の高周波反転パルスで生じるＳＥ（ｓｐｉｎ ｅｃｈ
ｏ）成分の位相θ_n+2,SE、ＳＴＥ（ｓｔｉｍｕｌａｔｅ
ｄｅｃｈｏ）成分の位相θ_n+2,STE は、

【数７】 で表される。

【００５７】図３の場合、１以上の全ての整数ｎに対し
て、φ_n ＝９０°、θ_n ＝９０°、θ_n+2,SE＝９０°、
θ_n+2,STE ＝９０°、である。したがって、位相ダイヤ
グラム上の複数のエコー経路の合成であるＥmain(n) 、
Ｅsub(n)の位相も等しく、それぞれ９０°となる。

【００５８】図４の場合は２通りに分けて考える。

【００５９】ｎ＝２ｍ （ｍは１以上の整数）のとき。

【００６０】

【数８】

【００６１】第ｎ番目の主エコー成分Ｅmain(n) からφ
_n+1 ，φ_n+2 に因り生じるＳＴＥ成分はＥsub(n+2)であ
るから、上記により、ｎ＝２以上の全てのｎについて図
４における副エコー成分Ｅsub(n)の位相は２７０°、つ
まり図３の状態から１８０°反転していることが分か
る。

【００６２】図８には隣接する高周波反転パルス間で生
じる、ある限られた数の経路についての位相の対応関係
しか示していないが、前述した図３８に示す規則性がエ
コー信号の位相の対応関係にもあるため、プリスキャン
Ｂの副エコー成分の位相はプリスキャンＡのそれに対し
て、１８０°回転する。

【００６３】本実施例では、ＣＰＭＧパルス系列の印加
位相の組み合わせを基本に変更しているが、これに限定
されるものではなく、プリスキャンＡ，Ｂの高周波パル
スの印加位相の組み合わせは、上記（＊）式を満たし、
主エコー成分Ｅmain(n) 又は副エコー成分Ｅsub(n)が２
つのプリスキャンで１８０°の位相差を持つようにする
ならば、どのような位相の組み合わせでもよい。

【００６４】再び、図２に戻って説明する。上述のプリ
スキャンの結果、プリスキャンＡで得られたエコー信号
をＥａ(n) 、プリスキャンＢで得られたエコー信号をＥ
ｂ(n) とすると（ｎ＝２，３，…）、

【数９】 の関係が成立する。そこで、コントローラ６はプリスキ
ャンＡ，Ｂで得られたエコー信号同士の平均化の処理を
行う（図２ステップＳ３）。つまり、

【数１０】｛Ｅａ(n) ＋Ｅｂ(n) ｝／２ 、（ここでｎ＝１，２，…）の演算を行う。この平均化
演算の概念を図５（ａ），（ｂ）および図６の模式図に
示す。これらの図に示すように、２番目以降のエコー信
号Ｅ(2) ，Ｅ(3) ，…では副エコー成分Ｅsub(2)，Ｅsu
b(3)，…の位相が２回のプリスキャンＡ，Ｂにおいて互
いに１８０°異なるので適宜相殺され、主エコー成分Ｅ
main(n) （ｎ＝１，２，…）のみが抽出される。

【００６５】なお、図６（ｂ）には、

【数１１】｛Ｅａ(n) −Ｅｂ(n) ｝／２ の演算を行った場合の、副エコー成分Ｅsub(n)（ｎ＝
２，３，…）のみが抽出される様子を模式的に示す。

【００６６】さらにコントローラ６は第１、第２番目の
エコー信号に基づいて主エコー成分のピークの純粋な位
相ずれφmain(1) 、φmain(2) および位置ずれｐmain
(1) 、ｐmain(2) を演算する（図２ステップＳ４）。

【００６７】次いで、従来の場合と同様に、

【数１２】φmain(1) −φmain(2) および

【数１３】ｐmain(1) −ｐmain(2) の値に基づいて、本スキャン時における、高周波パルス
ＲＦの印加位相φの補正量Δφ又はスライス用傾斜磁場
パルスＧssの補正量ΔＧss、およびリード方向傾斜磁場
パルスＧroの補正量ΔＧroを演算する（同図ステップＳ
５）。

【００６８】最後に、コントローラ６はマルチプレクサ
１１のスイッチ経路を再構成ユニット１２側に切換え、
シーケンサ５に対して図７に示すパルスシーケンスに基
づく本スキャンを実行させる（図２ステップＳ６）。

【００６９】本スキャンでは同図に記載しているよう
に、上述の補正演算結果を反映させたパルスシーケンス
が用いられる。すなわち、最初にスライス用傾斜磁場パ
ルスＧssと共に印加される高周波励起パルスＲＦexの印
加位相φがφ＝Δφの値であるか、または、このスライ
ス用傾斜磁場パルスＧssの後に印加する、極性反転した
ディフェイズ防止用のスライス方向傾斜磁場パルスＧss
の強度がΔＧssだけ補正された値となっている。また同
時に、そのディフェイズ防止用のスライス方向傾斜磁場
パルスＧssと並行して印加するリード方向傾斜磁場パル
スＧroの強度はΔＧroだけ補正された値となっている。

【００７０】そして最初の励起からτ／２の時間経過後
には、スライス用傾斜磁場パルスＧssと共に最初の高周
波反転パルスＲＦre1 （印加位相φ＝９０°、フリップ
角θ＝β１（ここでは１８０°）が印加され、スピンの
位相反転がなされる。次いで、このシーケンス実行毎に
強度（パルス面積）が調整される位相エンコード用傾斜
磁場パルスＧpeが印加される。そして、最初の選択励起
からのτ時間の経過に合わせてリード用傾斜磁場パルス
Ｇroが印加され、これと並行して第１番目のエコー信号
Ｅ(1) が収集される。

【００７１】以下同様に、スライス用傾斜磁場パルスＧ
ssと共に高周波反転パルスＲＦre2，ＲＦre3 ，…（印
加位相φ＝９０°、フリップ角θ＝β２，β３，…（こ
こでは１８０°）が順次印加されるととも、その間に、
スピンの位相反転に伴うエコー信号Ｅ(2) ，Ｅ(3) ，…
が順次収集される。なお、図７において、各エコー収集
後には、極性反転させた位相エンコード方向傾斜磁場パ
ルスＧpeがリワインディングパルスとして印加され、こ
れにより本スキャン時の疑似エコーの発生防止が図られ
る。

【００７２】本スキャンにより収集されたエコー信号Ｅ
（n)は順次、受信機８Ｒにて受信処理され、エコーデー
タに変換される。この処理には、エコー信号Ｅ（n)の直
交検波やＡ／Ｄ変換も含まれる。エコーデータは再構成
ユニット１２にてｋ（フーリエ）空間を形成するメモリ
上にエンコード量に応じて並べられ、フーリエ変換法に
よりＭＲ画像に再構成される。

【００７３】このように本実施形態では、従来の手法で
は除去しきれなかった、第２番目以降のエコー信号の副
エコー成分（第２番目のエコー信号に限って言及する
と、スティミュレーティッドエコーの成分）をほぼ完全
に除去でき、主エコー成分のみを的確に抽出できるの
で、本スキャンのパルスシーケンスに対する補正精度が
格段に向上し、安定した本スキャンを実施できるという
利点がある。

【００７４】なお、本実施形態における補正量Δφ、Δ
Ｇss、ΔＧroの算出法は前述したものに限定されること
なく、例えば、前記（１）式から明らかなように、第
１、２番目以外、すなわち第３番目以降の偶数番目、奇
数番目の主エコー成分の位相差およびエコーピーク位置
の差を、かかる補正量Δφ、ΔＧss、ΔＧroの算出に用
いてもよい。また、プリスキャンで得たエコーデータに
１次元のフーリエ変換をｋｒ方向について施し、読出し
方向の位相分布曲線を得て、その曲線の傾きと切片、つ
まり位相分布曲線の読出し方向の１次成分および０次成
分を求め、この１次成分および０次成分に基づいてエコ
ーピークの位置ずれ、位相ずれを求めるようにしてもよ
い。

【００７５】（第２の実施の形態）この発明に係る第２
の実施形態を図９〜図１３に基づき説明する。なお、こ
の実施形態を含め、以下の実施形態において上記第１の
実施形態と同一又は同等の構成要素には同一符号を用い
てその説明を省略又は簡略化する。

【００７６】この実施形態に係るＭＲＩ装置は上述した
ものと同一であり、コントローラ６は図９に示す一連の
処理を実行する。

【００７７】まず、コントローラ６はマルチプレクサ１
１のスイッチ経路をコントローラ６側に切り換えた状態
で、ＦＳＥ法に係る２種類のプリスキャンＣ，Ｄの実行
を順次、シーケンサ５に指令し、収集されたエコー信号
を入力する（図９ステップＳ１１，Ｓ１２）。最初のプ
リスキャンＣについては図１０に、後のプリスキャンＤ
については図１１に、それぞれ、パルスシーケンスを示
す。

【００７８】最初のプリスキャンＣでは、１つの高周波
励起パルスＲＦex（印加位相φ＝０°、フリップ角θ＝
α（ここでは９０°））に続いて複数の高周波反転パル
スＲＦre1 ，ＲＦre2 ，…（印加位相φ＝９０°、フリ
ップ角θ＝β１、β２、…（ここでは１８０°））が順
次、スライス用傾斜磁場パルスＧssと共に印加され、読
出し方向の傾斜磁場パルスＧroの印加と共に各エコー信
号Ｅｃ(n) が収集される。位相エンコード用傾斜磁場パ
ルスＧpeは、本スキャンでのｋ空間の中心を含むデータ
ラインを収集するときのショットと同じ強度変化に設定
されている。この結果、プリスキャンＣにおける高周波
パルスおよび傾斜磁場パルスは、本スキャンでのｋ空間
の中心を含むデータラインを収集するときのショットと
同じ波形に設定されている（この点は、従来および第１
の実施形態に係るプリスキャンのパルスシーケンスと異
なる）。

【００７９】いま、何等かの原因に因り原子核スピンに
位相ずれがあるとすると、前述と同様に、最初のエコー
信号Ｅ(1) 以外のエコー信号Ｅ(2) ，Ｅ(3) ，…のそれ
ぞれは、２つのエコー成分、すなわち主エコー成分Ｅma
in(2) （Ｅmain(3) ，…）と副エコー成分Ｅsub(2)（Ｅ
sub(3)，…）とに分かれる（図１０参照）。

【００８０】後のプリスキャンＤも、１つの高周波励起
パルスＲＦex（印加位相φ＝０°、フリップ角θ＝α
（ここでは９０°））に続いて複数の高周波反転パルス
ＲＦre1 ，ＲＦre2 ，…（印加位相φ＝９０°又は２７
０°、フリップ角θ＝β１、β２、…（ここでは１８０
°））が順次、スライス用傾斜磁場パルスＧssと共に印
加され、読出し方向の傾斜磁場パルスＧroの印加と共に
各エコー信号Ｅｄ(n) が収集される。位相エンコード用
傾斜磁場パルスＧpeは、本スキャンでのｋ空間の中心を
含むデータラインを収集するときのショットと同じ強度
変化に設定されている。この結果、プリスキャンＤにお
ける傾斜磁場パルスは、本スキャンでのｋ空間の中心を
含むデータラインを収集するときのショットと同じ波形
に設定されている（この点は、従来および第１の実施形
態に係るプリスキャンのパルスシーケンスと異なる）。
しかし、高周波パルスに関しては、前述した第１の実施
形態におけるプリスキャンＢと同様に、偶数番目の高周
波反転パルスＲＦre2 ，ＲＦre4 ，…の印加位相をさら
に１８０°だけ回転させ、φ＝２７０°に設定されてい
る。

【００８１】ここでも、何等かの原因に因り原子核スピ
ンに位相ずれがある場合、最初のエコー信号Ｅ(1) 以外
のエコー信号Ｅ(2) ，Ｅ(3) ，…のそれぞれは、２つの
エコー成分、すなわち主エコー成分Ｅmain(2) （Ｅmain
(3) ，…）と副エコー成分Ｅsub(2)（Ｅsub(3)，…）と
に分かれ（図１１参照）、しかも偶数番目の高周波反転
パルスＲＦre2 ，ＲＦre4 ，…の印加位相φを１８０°
だけ余計に回転させているため、偶数番目の副エコー成
分Ｅsub(2)，Ｅsub(3)，…の位相arg{Ｅsub(n)} だけが
プリンスキャンＣ時のそれに対して１８０°回転する。

【００８２】いま、プリスキャンＣで得られたエコー信
号をＥｃ(n) 、プリスキャンＤで得られたエコー信号を
Ｅｄ(n) とする。両方のプリスキャンにおいて実効エコ
ー時間に相当する位置、すなわち第ｎ_eff 番目のエコー
信号Ｅ（ｎ_eff ）（図１０、１１の例では第３番目のエ
コー信号Ｅ(3) ）を除くエコー信号には位相エンコード
が施されている。そこで、本スキャン時の補正量演算に
使用できるエコー信号は、実効エコー時間に相当する位
置、すなわち第３番目のエコー信号Ｅｃ（３）、Ｅｄ
（３）のみとなる。

【００８３】次いでコントローラ６は、第ｎ_eff 番目の
エコー信号Ｅ（ｎ_eff ）に相当する第３番目のエコー信
号Ｅ(3) の主エコー成分Ｅmain(3) 及び副エコー成分Ｅ
sub(3)を、

【数１４】 に基づく、

【数１５】 の演算を行って求める（図９ステップＳ１３）。この演
算により抽出された主エコー成分Ｅmain(3) 及び副エコ
ー成分Ｅsub(3)を図１２に示す。本実施形態では第ｎ
_eff 番目のエコー信号Ｅ（ｎ_eff ）に相当する第３番目
のエコー信号Ｅ(3)の内、副エコー信号Ｅsub(3)同士の
位相が１８０°異なるため、上述のように、主エコー信
号Ｅmain(3) は加算平均により抽出され、副エコー信号
Ｅsub(3)は引算平均により抽出される。

【００８４】さらにコントローラ６はこれらのエコーデ
ータに基づいて第３番目のエコー信号Ｅ(3) の主エコー
成分のピークの位相ずれφmain(3) 及び位置ずれｐmain
(3)ならびに副エコー成分のピークの位相ずれφsub(3)
及び位置ずれｐsub(3)を演算する（図９ステップＳ１
４）。

【００８５】次いで、従来手法と同様に、同一番目ｎの
主エコー成分と副エコー成分の位相差

【数１６】φmain（ｎ_eff ）−φsub （ｎ_eff ） 、および、ピーク位置の差

【数１７】ｐmain（ｎ_eff ）−ｐsub （ｎ_eff ） の値、すなわち、ここでは、

【数１８】 の値に基づいて、本スキャン時における、高周波パルス
ＲＦの印加位相φの補正量Δφ又はスライス用傾斜磁場
パルスＧssの補正量ΔＧss、およびリード方向傾斜磁場
パルスＧroの補正量ΔＧroを演算する（同図ステップＳ
１５）。なお、第ｎ_eff 番目が奇数か偶数かにより上記
補正のための演算結果の極性が反転するので、常に同一
極性となるように減算の順を統一する。

【００８６】最後に、コントローラ６はマルチプレクサ
１１のスイッチ経路を再構成ユニット１２側に切換え、
シーケンサ５に対して図１２に示すパルスシーケンスに
基づく本スキャンを実行させる（図９ステップＳ１
６）。

【００８７】本スキャンでは同図に記載しているよう
に、上述の補正演算結果を反映させたパルスシーケンス
が用いられる。すなわち、最初にスライス用傾斜磁場パ
ルスＧssと共に印加される高周波励起パルスＲＦex（フ
リップ角θ＝α（ここでは９０°）の印加位相φがφ＝
Δφの値であるか、または、このスライス用傾斜磁場パ
ルスＧssの後に印加する、極性反転したディフェイズ防
止用のスライス方向傾斜磁場パルスＧssの強度がΔＧss
だけ補正された値となっている。また同時に、そのディ
フェイズ防止用のスライス方向傾斜磁場パルスＧssと並
行して印加するリード方向傾斜磁場パルスＧroの強度は
ΔＧroだけ補正された値となっている。

【００８８】以下、前述した図７と同様に、複数の高周
波反転パルスＲＦre1 ，ＲＦre2 ，…を使って、原子核
スピンの位相反転に伴うエコー信号Ｅ(1) ，Ｅ(2) ，…
が順次収集される。この本スキャンにより収集されたエ
コー信号Ｅ（n)は順次、受信機８Ｒを介してエコーデー
タとして再構成ユニット１２に入力し、フーリエ変換法
によりＭＲ画像に再構成される。

【００８９】このように本実施形態によれば、前述した
第１の実施形態における利点に加えて、位相エンコード
方向の傾斜磁場パルスに起因する位相ずれ成分の内、空
間的に一様な位相ずれ成分、つまり０次成分を好適に補
正することができ、本スキャンのパルスシーケンスに対
する補正精度がより一層向上するという利点が得られ
る。

【００９０】なお、本実施形態における補正量Δφ、Δ
Ｇss、ΔＧroの算出法は前述したものに限定されること
なく、例えば、プリスキャンで得たエコーデータに１次
元のフーリエ変換をｋｒ方向について施し、読出し方向
の位相分布曲線を得て、その曲線の傾きと切片、つまり
位相分布曲線の読出し方向の１次成分および０次成分を
求め、この１次成分および０次成分に基づいてエコーピ
ークの位置ずれ、位相ずれを求めるようにしてもよい。

【００９１】（第２の実施の形態の変形）上記第２の実
施形態の変形形態を図１４〜図１８に基づき説明する。
この変形形態は、上述したＦＳＥ法に代えて、ＧＲＡＳ
Ｅ法を適用したものであり、シーケンサ５はコントロー
ラ６からの指令に応答してプリスキャン時及び本スキャ
ン時ともにＧＲＡＳＥ法のパルスシーケンスを指令でき
るようになっている。

【００９２】コントローラ６は図１４に示す一連の処理
を実行する。

【００９３】まず、コントローラ６はマルチプレクサ１
１のスイッチ経路をコントローラ６側に切り換えた状態
で、ＧＲＡＳＥ法に係る２種類のプリスキャンＥ，Ｆの
実行を順次、シーケンサ５に指令し、収集されたエコー
信号を入力する（図１４ステップＳ２１，Ｓ２２）。最
初のプリスキャンＥについては図１５に、後のプリスキ
ャンＦについては図１６に、パルスシーケンスをそれぞ
れ示す。

【００９４】この２種類のプリスキャンＥ，Ｆのパルス
シーケンスは、第２の実施形態に係る図１０、１１のそ
れと同一の１つの高周波励起パルスＲＦexおよび複数の
高周波反転パルスＲＦre1 ，ＲＦre2 ，…を夫々使用
し、かつ、高周波反転パルス間において読出し用傾斜磁
場パルスＧroを３回極性反転してエコー信号Ｅe(m,n)、
Ｅf(m,n)を夫々収集する。このときの位相エンコード用
傾斜磁場パルスＧpeは本スキャンでのｋ空間の中心を含
むデータラインを収集するときのショットと同じ強度変
化に夫々設定されている。

【００９５】この両方のプリスキャンＥ，Ｆにおいて、
何等かの原因に因り原子核スピンに位相ずれがある場
合、前述と同様に、最初のエコー信号Ｅ(1,1) 〜Ｅ(1,
3) 以外のエコー信号Ｅ(2,1) ，Ｅ(2,2) ，(2,3) ，Ｅ
(3,1) ，…のそれぞれは、２つのエコー成分、すなわち
主エコー成分Ｅmain(2,1) （Ｅmain(2,2) ，Ｅmain(2,
3)，Ｅmain(3,1) ，…）と副エコー成分Ｅsub(2,1)（Ｅ
sub(2,2)，Ｅsub(2,3)，Ｅsub(3,1)，…）とに分かれる
（その一例を、エコー信号Ｅ(3,2) について図１５〜１
７に示す）。しかも、後の方のプリスキャンＦでは、偶
数番目の高周波反転パルスＲＦre2 ，ＲＦre4 ，…の印
加位相φを１８０°だけ余計に回転させ、φ＝２７０°
としているため、偶数番目の副エコー成分の位相だけが
プリスキャンＥ時のそれに対して１８０°回転する。

【００９６】いま、プリスキャンＥで得られたエコー信
号をＥe(n,m)、プリスキャンＦで得られたエコー信号を
Ｅf(n,m)とする。両方のプリスキャンにおいて実効エコ
ー時間に相当する位置、すなわち第ｎ_eff 番目のエコー
信号Ｅ（ｎ_eff ,m' ）（図１５、１６の例では第３番目
のエコー信号中の信号Ｅ(3,2) ）を除くエコー信号には
位相エンコードが施されている。そこで、本スキャン時
の補正量演算に使用できるエコー信号は、実効エコー時
間に相当する位置、すなわち第３番目のエコー信号中の
信号Ｅ(3,2) のみとなる。

【００９７】次いでコントローラ６は、ここではエコー
信号Ｅ(3,2) の主エコー成分Ｅmain(3.2) 及び副エコー
成分Ｅsub(3,2)を、

【数１９】 の演算を行って求める（図１４ステップＳ２３）。この
演算により抽出された主エコー成分Ｅmain(3,2) 及び副
エコー成分Ｅsub(3,2)を図１７に示す。

【００９８】さらにコントローラ６はこれらのエコーデ
ータに基づいてエコー信号Ｅ(3,2)の主エコー成分のピ
ークの位相ずれφmain(3,2) 及び位置ずれｐmain(3,2)
ならびに副エコー成分のピークの位相ずれφsub(3,2)及
び位置ずれｐsub(3,2)を演算する（図１４ステップＳ２
４）。

【００９９】次いで、従来手法と同様に、同一番目ｎの
主エコー成分と副エコー成分の位相差の値、すなわち、
ここでは、

【数２０】 の値に基づいて、本スキャン時における、高周波パルス
ＲＦの印加位相φの補正量Δφ又はスライス用傾斜磁場
パルスＧssの補正量ΔＧss、およびリード方向傾斜磁場
パルスＧroの補正量ΔＧroを演算する（同図ステップＳ
２５）。

【０１００】最後に、コントローラ６はマルチプレクサ
１１のスイッチ経路を再構成ユニット１２側に切換え、
シーケンサ５に対して図１８に示すパルスシーケンスに
基づく本スキャンを実行させる（図１４ステップＳ２
６）。本スキャンでは同図に記載しているように、最初
にスライス用傾斜磁場パルスＧssと共に印加される高周
波励起パルスＲＦexの印加位相φがφ＝Δφの値である
か、または、このスライス用傾斜磁場パルスＧssの後に
印加する、極性反転したディフェイズ防止用のスライス
方向傾斜磁場パルスＧssの強度がΔＧssだけ補正された
値となっている。また同時に、そのディフェイズ防止用
のスライス方向傾斜磁場パルスＧssと並行して印加する
リード方向傾斜磁場パルスＧroの強度はΔＧroだけ補正
された値となっている。

【０１０１】このように本実施形態によれば、前述した
第１及び第２の実施形態における利点が得られる。さら
に、ＧＲＡＳＥ法はＦＳＥ法とは異なり、位相エンコー
ド用傾斜磁場パルスの波形の形状は高周波反転パルス間
であまり変化させないことから、位相エンコードに拠る
位相誤差の発生パターンが規則的となり、したがって図
１８に斜線で示すΔＧss、ΔＧroの補正によってアーチ
ファクトを効果的に低減させることができる。

【０１０２】（第３の実施の形態）さらに、第３の実施
形態に係るＭＲＩ装置を図１９〜図２３を参照して説明
する。この実施形態においても２種類のプリスキャン
Ｇ，Ｈを使用するが、そのスキャン時に位相エンコード
量を本スキャンと同様に各エコーに施すようにしたもの
である。ここでは、パルスシーケンスがＧＲＡＳＥ法で
ある場合について説明する。

【０１０３】このＭＲＩ装置も前述した各実施形態のも
のと同一の構成であり、コントローラ６は図１９に示す
処理を実行する。

【０１０４】コントローラ６は、マルチプレクサ１１の
スイッチ経路をコントローラ６側に切り換え、ＧＲＡＳ
Ｅ法に係る２種類のプリスキャンＧ，Ｈの実行を順次、
シーケンサ５に指令し、収集されたエコー信号を入力す
る（図１９ステップＳ３１，Ｓ３２）。

【０１０５】図２０に示す最初のプリスキャンＧ及び図
２１に示す次のプリスキャンＨは共にこれまでのものと
は異なり、位相エンコード用傾斜磁場パルスＧpeを本ス
キャン時と同じように各エコーＥg(n,m)又はＥh(n,m)に
施す位相エンコード量を各ショット毎に変化させる。こ
の位相エンコードの付与は、ｋ空間上の全面又は一部に
当たる未調整状態の２次元のエコーデータを収集するた
めである。高周波パルスＲＦex、ＲＦre1,ＲＦre2,…の
印加位相は、前述した第１の実施形態の２種類のプリス
キャンＡ，Ｂの場合と同じに変化させる。このため、２
次元のエコーデータは高周波パルスの印加位相の組み合
わせを変えて２組、すなわちＥg(n,m)又はＥh(n,m)、収
集される。

【０１０６】上記プリスキャンＧ，Ｈでは、一例とし
て、ＧＲＡＳＥ法の高周波反転パルス数＝５、隣接する
高周波反転パルス間のグラジェントエコー数＝３とし、
１ショット当たりの合計エコー数＝１５のパルスシーケ
ンスの場合を示している。図面上には、表示の関係で、
２ショット、画像の位相エンコード方向のマトリクス＝
３０のときの傾斜磁場波形のみを表してあるが、ショッ
ト数は、実効エコー時間に相当する第ｎ_eff 番目のエコ
ー信号のみで、画面の位相エンコード方向の１次の位相
変化を安定して測定できるだけのショット数、例えば８
ショット以上（これをｋ空間上のデータで表現すると、
ｋ空間の中心位置を含むデータラインを中心として、隣
接するエコーデータライン数が８ライン以上）であるこ
とが望ましい。

【０１０７】いま、何等かの原因に因り原子核スピンに
位相ずれがある場合、最初のエコー信号Ｅ(1,m) 以外の
エコー信号Ｅ(2,m) ，Ｅ(3,m) ，…のそれぞれは、２つ
のエコー成分、すなわち主エコー成分Ｅmain(2,m) （Ｅ
main(3,m) ，…）と副エコー成分Ｅsub(2,m)（Ｅsub(3,
m)，…）とに分かれる（例えば、図２０、２１中のエコ
ー信号Ｅg(3,2)、Ｅh(3,2)参照）。しかも後の方のプリ
スキャンＨでは、偶数番目の高周波反転パルスＲＦre2
，ＲＦre4 ，…の印加位相φを１８０°だけ余計に回
転させているため、偶数番目の副エコー成分Ｅsub(2,
m)，Ｅsub(4,m)，…の位相だけがプリンスキャンＧ時の
それに対して１８０°回転する。

【０１０８】前述した各実施形態のときの同様に、図２
２に示す如く、最初のプリスキャンＧで得られた２次元
のエコー信号をＥg(n,m)、次のプリスキャンＨで得られ
たエコー信号をＥh(n,m)とすると、ｋ平面上の２次元の
主エコー成分Ｅmain(n,m) 、２次元の副エコー成分Ｅsu
b(n,m)は、

【数２１】 で抽出される。

【０１０９】このため、両方のプリスキャンにおいて実
効エコー時間に相当する位置、すなわち第ｎ_eff 番目の
エコー信号Ｅ(n,m) を例えばＥ(3,2) とすると、

【数２２】 により主エコー成分Ｅmain(3,2) 及び副エコー成分Ｅsu
b (3,2) が抽出される。（図１９ステップＳ３３）。抽
出されたエコー成分Ｅmain(3,2) 、Ｅsub (3,2)を図２
２（ｃ），（ｄ）にそれぞれ例示する。第ｎ_eff 番目に
相当する２組の２次元エコー信号Ｅ(3,2) の内、副エコ
ー信号Ｅsub(3,2)同士の位相が１８０°異なるため（図
２１、２２では一方を点線で示す）、上述のように、主
エコー信号Ｅmain(3,2) は加算平均により抽出され、副
エコー信号Ｅsub(3,2)は引算平均により抽出される。

【０１１０】さらにコントローラ６はこれらのｋ平面上
の２次元のエコーデータＥmain(3,2) 、Ｅsub(3,2)に基
づいて第(3,2) 番目のエコー信号Ｅ(3,2) の主エコー成
分のピークの位相ずれφmain(3.2) 及び位置ずれｐmain
(3,2) ならびに副エコー成分のピークの位相ずれφsub
(3,2)及び位置ずれｐsub(3,2)を演算する（図１９ステ
ップＳ３４）。ここで、位置ずれｐmain(3,2) 、ｐsub
(3,2)は共に２次元のベクトル量である。

【０１１１】次いで、前記（１），（２）式に基づいて
同一(3,2) 番目の主エコー成分と副エコー成分の位相差

【数２３】φmain(3,2) −φsub(3,2) 、および、ピーク位置ベクトルの差

【数２４】ｐmain(3,2) −ｐsub(3,2) を演算し、この値に基づいて本スキャン時における、高
周波パルスＲＦの印加位相φの補正量Δφ、スライス用
傾斜磁場パルスＧssの補正量ΔＧss、リード用傾斜磁場
パルスＧroの補正量ΔＧro、位相エンコード用傾斜磁場
パルスＧpeの補正量ΔＧpeをそれぞれ演算する（同図ス
テップＳ３５）。なお、第ｎ_eff 番目が奇数か偶数かに
より上記補正のための演算結果の極性が反転するので、
常に同一極性となるように減算の順を統一する。

【０１１２】最後に、コントローラ６はマルチプレクサ
１１のスイッチ経路を再構成ユニット１２側に切換え、
シーケンサ５に対して図２３に示すパルスシーケンスに
基づく本スキャンを実行させる（図１９ステップＳ３
６）。

【０１１３】本スキャンでは同図に記載しているよう
に、上述の補正演算結果を反映させたパルスシーケンス
が用いられる。すなわち、最初にスライス用傾斜磁場パ
ルスＧssと共に印加される高周波励起パルスＲＦex（フ
リップ角θ＝α（ここでは90°）の印加位相φがφ＝Δ
φの値であるか、または、このスライス用傾斜磁場パル
スＧssの後に印加する、極性反転したディフェイズ防止
用のスライス用傾斜磁場パルスＧssの強度がΔＧssだけ
補正された値となっている。また同時に、そのディフェ
イズ防止用のスライス方向傾斜磁場パルスＧssと並行し
て印加するリード用傾斜磁場パルスＧroの強度はΔＧro
だけ補正された値で、かつ位相エンコード用傾斜磁場パ
ルスＧpeの強度はΔＧpeだけ補正された値となってい
る。

【０１１４】以下、前述した実施形態と同様に、複数の
高周波反転パルスＲＦre1 ，ＲＦre2 ，…を使って、原
子核スピンの位相反転に伴うエコー信号Ｅ(n,m) が順次
収集され、ＭＲ画像に再構成される。

【０１１５】このように本実施形態によれば、第１、第
２の実施形態における利点に加えて、位相エンコード方
向の傾斜磁場パルスに起因する位相ずれ成分の内、位相
エンコード方向に１次の空間分布を持つ成分までをも好
適に補正することができ、本スキャンのパルスシーケン
スに対する補正精度を著しく向上させ、高画質のＭＲ画
像を提供できるという利点が得られる。

【０１１６】なお、本実施形態における補正量Δφ、Δ
Ｇ（ΔＧss、ΔＧro、ΔＧpe）の算出法は前述したもの
に限定されない。例えば、プリスキャンで得た２次元の
エコーデータから一度、２次元画像を再構成し、読出し
方向及び位相エンコード方向の位相分布曲線の傾きと切
片、つまり位相の読出し方向及び位相エンコード方向の
１次成分及び０次成分を求める。そして、この１次成分
および０次成分に基づいてエコーピークの位置ずれ、位
相ずれを求め、本スキャンのシーケンスの補正量Δφ、
ΔＧを算出するようにしてもよい。

【０１１７】また、上記各実施形態は２次元スキャンで
説明してきたが、３次元スキャンやそれ以上の次元のス
キャンについても、ベクトルの次元を増加させることで
同様に補正できる。

【０１１８】（第４の実施の形態）さらに、第４の実施
形態に係るＭＲＩ装置を図２４、２５に基づき説明す
る。この実施形態は、２種類のプリスキャンを行うこと
は従前と同じであるが、プリスキャンの特質を考慮し、
プリスキャン時にはその特質に合致した撮像対象を使用
するようにしたものである。

【０１１９】使用するＭＲＩ装置は第１の実施形態のも
のと同一の構成であり、コントローラ６は図２４及び図
２５に示す一連の処理を実行する。図２４はプリスキャ
ン時の処理を、図２５は本スキャン時の処理をそれぞれ
示している。

【０１２０】プリスキャンは、傾斜磁場コイルのような
ハードウエアの不完全性に因る主エコー成分と副エコー
成分の位相誤差を測定することが目的である。このた
め、プリスキャン時の被検体（撮像対象）は実際のＭＲ
画像を得る本スキャン時の被検体と同一である必要な無
い。そこで、この実施形態ではプリスキャン時に実際の
被検体の代わりにファントムを使用する。ファントムと
しては、測定精度を確保する必要上、なるべく信号値分
布が一様なものが望ましく、例えばアクリル樹脂容器に
硫酸銅水溶液を満たしたファントムが良い。

【０１２１】このようなファントムを磁石１の診断用空
間に設置した後、コントローラ６に指令を与えること
で、コントローラ６は図２４のプリスキャン時の処理を
開始する。

【０１２２】コントローラ６はまず、図２４のステップ
Ｓ４１で、２種類のプリスキャンＧ，Ｈの撮像条件であ
る分解能及びスライス厚を決定する。プリスキャンＧ，
Ｈは前述した第３の実施形態で説明したものと同じで、
ＧＲＡＳＥ法に係るパルスシーケンスを使用する。次い
でステップＳ４２にて、決定した撮像条件にしたがった
２種類のプリスキャンＧ，Ｈを順次実行させる。次いで
ステップＳ４３にて、第３の実施形態のときと同様に演
算処理を行い、第ｎ_eff 番目（＝(n,m) 番目）に相当す
るエコー信号Ｅ(n,m) の主エコー成分Ｅmain(n,m) 、副
エコー成分Ｅsub (n,m) を、

【数２５】 の式に基づき抽出する。

【０１２３】さらにステップＳ４４に移行し、主エコー
成分Ｅmain(n,m) 、副エコー成分Ｅsub (n,m) それぞれ
の位相誤差の空間分布の０次、１次の成分φmain(n,m)
、ｐmain(n,m) 、φsub(n,m)、ｐsub(n,m)を求める。
さらにステップ４５で、前記（１），（２）式に基づく
主エコー成分と副エコー成分の位相差、および、ピーク
位置ベクトルの差を演算する。ステップ４６では、これ
らの差情報を、プリスキャンの空間分解能、スライス厚
の情報とともに補正用データとして例えば、記憶ユニッ
ト１３に格納する。

【０１２４】本スキャンが指令されると、コントローラ
６は図２５の処理を開始する。最初に図２５のステップ
Ｓ５１にて、オペレータから与えられた情報を参照する
などの手法で、本スキャンの撮像条件である分解能、ス
ライス厚を決める。次いでステップＳ５２に移行し、予
めプリスキャン時に記憶した、プリスキャンＧ，Ｈの撮
像条件（分解能、スライス厚）並びに主エコー成分と副
エコー成分の位相差およびピーク位置ベクトルの差を補
正用データとして読み出す。

【０１２５】そして、ステップＳ５３で、読み出して補
正用データに基づいて前述と同様に補正量Δφ、ΔＧを
算出する。このとき、本スキャンの撮影条件、画像の各
方向の空間分解能をプリスキャン時のそれらと比較し、
それらの比率に応じて補正量Δφ、ΔＧを求めるように
する。さらに、ステップＳ５４にて、補正量Δφ、ΔＧ
を反映させたパルスシーケンスにより本スキャンの実行
をシーケンサ５に指令する。

【０１２６】このように理想的な信号値分布のファント
ムを使ってプリスキャンを実施することで、正確度を増
した補正データを得ることができる。プリスキャン時に
撮像対象が動いて、位相誤差の補正データの正確度が下
がることも無く、撮像対象からの信号の変化に影響を受
けることもなく、安定かつ正確な補正データを収集し
て、本スキャンを実施できる。患者にとっては、位相誤
差の補正データを収集するためのプリスキャンが不要に
なるから、本スキャンのみで済み、患者の拘束時間が減
り、診断を受け易くもなる。

【０１２７】なお、本実施形態において、位相エンコー
ド方向の位相誤差分布が無視できるほど小さい場合、第
１の実施形態におけるプリスキャンＡ，Ｂのように位相
エンコードを零としたプリスキャンを行い、位相エンコ
ード方向以外の方向の位相誤差分布だけを測定し、この
位相誤差情報に基づいて補正を実施するようにしてもよ
い。また必要に応じて、第２の実施形態のプリスキャン
Ｃ，Ｄでみられたように、位相エンコード方向の傾斜磁
場に因る０次の位相誤差成分も含めた位相誤差分布の得
るプリスキャンを実施してもよい。

【０１２８】（第４の実施の形態の変形）上記第４の実
施形態をさらに変形した例を図２６、図２７に基づき説
明する。この変形形態に係るＭＲＩ装置は、傾斜磁場の
各物理チャンネルの位相誤差特性を予め２種類のプリス
キャンＧ，Ｈ（第４の実施形態と同一のプリスキャン）
により収集しておき、本スキャンのときには、断面方
向、分解能、スライス厚などの傾斜磁場に関連した撮影
条件からパルスシーケンスの補正量Δφ、ΔＧを算出
し、本スキャンを実施するものである。これにより、後
述するように、本スキャンと同じ断面方向のプリスキャ
ンを実行しなくても任意方向の断面の撮像時におけるシ
ーケンスの補正量Δφ、ΔＧを求めることができる。

【０１２９】コントローラ６は、前述の図２４の処理に
代えて、図２６、２７の処理を実施するとともに、前述
した図２５の処理を実施する。

【０１３０】図２６の処理は、プリスキャンを含むもの
で、補正用傾斜磁場係数マトリクスＫp,ro、Ｋp,pe，Ｋ
p,ssの各要素（配列αで表示）の算出手順を示す。図２
７の処理は、補正用高周波位相係数ベクトルＫφ,ro 、
Ｋφ,pe 、Ｋφ,ss の各要素（配列βで表示）の算出手
順を示す。

【０１３１】また、補正用傾斜磁場係数マトリクスＫp,
ro、Ｋp,pe，Ｋp,ssの各要素と配列αの各要素との対応
関係は、以下の式（９）により表される。

【０１３２】

【数２６】 補正用高周波位相係数ベクトルＫφ,ro 、Ｋφ,pe 、Ｋ
φ,ss の各要素と配列βの各要素との対応関係は、以下
の式（１０）により表される。

【０１３３】

【数２７】 補正用傾斜磁場係数マトリクスＫp,ro、Ｋp,pe，Ｋp,ss
と補正用傾斜磁場ΔＧとの対応関係は、以下の式（１
１）により表される。

【０１３４】

【数２８】 補正用高周波位相係数ベクトルＫφ,ro 、Ｋφ,pe 、Ｋ
φ,ss と高周波パルスの印加位相の補正量Δφとの対応
関係は、以下の式（１２）により表される。

【０１３５】

【数２９】 さらに、撮像断面方向の回転を表すマトリクスＲの定義
を式（１３）に、補正用傾斜磁場ΔＧの各要素を式（１
４）に夫々示す。

【０１３６】

【数３０】

【数３１】

【０１３７】また、Ｍ,imgＣh をイメージングチャンネ
ル imgＣh 方向（パルスシーケンスに関係するRO,PE,SS
の方向）の空間周波数の大きさとし、Ｍ, phＣh を物理
的なチャンネルphＣh 方向（Ｘ，Ｙ，Ｚの傾斜磁場コイ
ルに関係した物理的なチャンネルの方向）の空間周波数
の大きさとすると、それらと撮像断面の方向を示す回転
マトリクスＲとは前述した式（１３）で結合されている
とする。

【０１３８】式（９），（１１）に示す如く、

【数３２】α［ imgＣh ］［phＣh1］［phＣh2］ は、 imgＣh をphＣh1に対応させると、phＣh2方向の１
次の位相分布をどの程度補正するかを示す係数である。
式（１１）中のＴは補正用傾斜磁場の時間幅であり、同
式中のＦ,imgＣh は添字のイメージングチャンネル img
Ｃh について、基準とする撮像条件を得るために必要な
強度の傾斜磁場に対する実際のスキャンの傾斜磁場強度
の比であり、同式中のＭo, imgＣh はその imgＣh 方向
の基準とする空間周波数の大きさである。

【０１３９】また式（１０），（１２）に示す如く、

【数３３】β［ imgＣh ］［phＣh ］ は、パルスシーケンスに関係するRO,PE,SSのイメージン
グチャンネル imgＣh を物理チャンネルphＣh に対応さ
せると、０次の位相分布をどの程度補正させるかを示す
係数である。式（１２）中のＦ,imgＣh は式（１１）と
同じ値である。

【０１４０】続いて、コントローラ６によって実施され
る、図２６に示す補正用傾斜磁場係数マトリクスの各要
素の算出手順を説明する。まず、マトリクスの内のどの
要素を求めるかを選択する（同図ステップＳ６１）。次
いで、その要素に関係する断面方向を選択する（ステッ
プＳ６２）。マトリクスのある要素を求めるには、撮像
条件をいくつか変え、それぞれの撮像条件におけるスピ
ンの位相ずれ分布を求め、比較する必要がある。

【０１４１】そこで、収集回数カウント用の変数ｉ＝０
に初期設定した後（ステップＳ６３）、目的とするマト
リクス要素が抽出できるように空間分解能、スライス厚
などの撮像条件を決める（ステップＳ６４）。傾斜磁場
コイルあるいはシーケンスの特性上、いくつかのマトリ
クス要素が明らかに零と見做せる場合、収集回数を減ら
すことも可能である。

【０１４２】次いで、前述した第４の実施形態に係る図
２４のプリスキャンＧ，Ｈ（図２４のステップＳ４２〜
Ｓ４４を含む一連の処理）を実施する（図２６ステップ
Ｓ６５）。次いで、プリスキャンにより求めた主エコー
成分及び副エコーせいぶんの位置ずれを示す２次元ベク
トルｐmain及びｐsub から、それらのピーク位置ベクト
ルの差を、変数ｉ毎に、

【数３４】ΔＰi ＝ｐmain−ｐsub より演算する（ステップＳ６６）。

【０１４３】次いでコントローラ６はステップＳ６７に
処理を移行させて、撮像条件が残っているかどうかを判
断し、未だ撮像条件があるならば変数ｉをインクリメン
トし（ｉ＝ｉ＋１）、ステップＳ６４に戻る。これによ
り、ステップＳ６４〜Ｓ６６の処理が１回又は複数回実
行され、異なる（又は一つの）撮像条件で上記ピーク位
置ベクトルの差ΔＰｉが求められる。

【０１４４】上記ステップＳ６７の判断でＮＯとなり、
予め定めた撮像条件におけるプリスキャンが全部終了し
たと判断すると、コントローラ６は次いで、撮像条件の
情報とピーク位置ベクトルの差ΔＰｉとに基づいて対応
するマトリクス要素の補正量

【数３５】Δα［ imgＣh ］［phＣh1］［phＣh2］ を決める（ステップＳ６９）。この補正量Δα［ imgＣ
h ］［phＣh1］［phＣh2］を用いて、対応するマトリク
ス要素

【数３６】 を求める（ステップＳ７０）。この一連の処理を、ステ
ップＳ６１で予め選択した要素α全てについて実施する
（ステップＳ７１）。

【０１４５】一方、図２７に示す補正用高周波位相係数
ベクトルの各要素も上述の図２６の場合と同様の流れで
算出される（図２７ステップＳ８１〜９１参照）。ステ
ップＳ８１では補正用高周波位相係数ベクトルの求める
要素β［ imgＣh ］［phＣh］が予め選択される。ステ
ップＳ８６では、前述したプリスキャンＧ，Ｈ及びその
後処理演算に基づいて主エコー成分と副エコー成分の位
相差

【数３７】Δφｉ＝φmain−φsub が撮像条件語とに演算される。さらにステップＳ８９
で、撮像条件の情報と位相差Δφｉとに基づいて対応す
るベクトル要素の補正量

【数３８】Δβ［ imgＣh ］［phＣh ］ が決められ、ステップＳ９０で、対応するベクトル要素

【数３９】 が求められる。

【０１４６】このように求められた補正用傾斜磁場係数
マトリクス及び補正用高周波位相係数ベクトルは、記憶
ユニット１３又は内部メモリに格納される。そして、本
スキャンの際、図２５で説明した処理が実行される。つ
まり、本スキャンの撮像条件を決定し（同図ステップＳ
５１）、格納していた補正用傾斜磁場係数マトリクス及
び補正用高周波位相係数ベクトルを読出し（同図ステッ
プＳ５２）、式（１１）及び（１２）に基づいて印加位
相の補正量Δφ及び傾斜磁場の補正量ΔＧが算出される
（ステップＳ５３）。これらの補正量Δφ、ΔＧを加味
したパルスシーケンスの本スキャンが実行される（ステ
ップＳ５４）。

【０１４７】このように処理することで、本スキャンと
同じ断面方向のプリスキャンが不要になり、そのような
プリスキャンを行うことなく、任意方向の断面を撮像
（本スキャン）できる。

【０１４８】この変形形態に係る撮像手法の特質を従来
技術との対比で詳述する。従来技術の項で述べた米国特
許第５３７８９８５号、及び特開平６−５４８２７号公
報の記載手法において補正を行っている要素は、前述の
補正用傾斜磁場係数マトリクスの内、Ｋp,ro、Ｋp,ssの
要素、及び補正用高周波位相係数ベクトルの内、Ｋφ,r
o 、Ｋφ,ss の要素に相当すると考えてよい。各要素の
分離も明確ではないため、プリスキャンは本スキャンの
直前に実施し、しかも本スキャンと同じ断面方向、空間
分解能、スライス厚に設定する必要があった。

【０１４９】これに対し、この変形形態では、従来の要
素に加えて、従来では補正不可能であった補正用傾斜磁
場係数マトリクスのＫp,peの対角要素、及び補正用高周
波位相係数ベクトルＫφ,pe も補正される。この補正さ
れる要素の増加に加えて、本変形形態では、全ての補正
用傾斜磁場係数マトリクスの非対角要素がそれぞれ独立
に測定及び補正できるため、本スキャンとプリスキャン
との撮像断面方向が異なる場合でも補正量が算出され
る。しかも、高周波パルスの位相と傾斜磁場を補正する
ことで、かかるマトリクス及びベクトルのほとんどの要
素が補正される。ゆえに、同一の傾斜磁場コイルのＭＲ
Ｉ装置であっても、本変形形態にかかるＭＲＩ装置の方
が従来よりも格段に安定した高画質のＭＲ画像を再構成
できる。

【０１５０】なお、上述した式（１１），（１２）など
の演算モデルは、傾斜磁場システムの特性に応じて適宜
変更可能なものである。

【０１５１】（第５の実施の形態）さらに第５の実施形
態の係るＭＲＩ装置を図２８〜図３３に基づいて説明す
る。この実施形態は、従前のものとは異なり、本スキャ
ンにおけるパルスシーケンスの高周波パルスの印加位相
や傾斜磁場の波形を補正することはせずに、２次以上の
高次の位相分布誤差を含む位相誤差を的確に補正したＭ
Ｒ画像を得るようにしたものである。

【０１５２】ＭＲＩ装置のハードウエアはこれまでのも
のと同一の構成であり、コントローラ６は図２８に示
す、プリスキャン及び本スキャンを組み合わせた一連の
処理を実行する。ここではＦＳＥ法を使っている場合を
例示する。

【０１５３】まず、コントローラ６はマルチプレクサ１
１をコントローラ側に切り換えた状態で、２種類のプリ
スキャンＡ，Ｂをシーケンサ５に順次実行させる（同図
ステップＳ９１）。このプリスキャンＡ，Ｂは第１の実
施形態で説明したものと同一に設定してある。

【０１５４】そして、プリスキャンＡ，Ｂで得られたエ
コー信号Ｅa(n)、Ｅb(n)から前述と同様に、主エコー成
分Ｅmain(n) 、副エコー成分Ｅsub(n)を

【数４０】 を抽出するとともに、その絶対値｜｛Ｅa(n)＋Ｅb(n)｝
／２｜、｜｛Ｅa(n)−Ｅb(n)｝／２｜を演算し、その絶
対値を主エコー成分Ｅmain(n) 、副エコー成分Ｅsub(n)
の振幅情報ＥＰmain、ＥＰsub として設定する（ステッ
プＳ９２；図３２参照）。

【０１５５】次いでコントローラ６は、マルチプレクサ
１１を再構成ユニット側に切り換え、２種類の本スキャ
ンＩ，Ｊをシーケンサ５に順次実行させる（ステップＳ
９３）。この本スキャンＩ，Ｊのパルスシーケンスは図
２９、３０に示す。これらの図に示すように、高周波パ
ルスの印加位相のパターン以外は通常のＦＳＥ法と同一
であり、前述した実施形態のように、高周波パルスの印
加位相や傾斜磁場波形を補正するようにはしていない。
最初の本スキャンＩにおける高周波パルスの印加位相パ
ターンは、通常のＦＳＥ法のときと同一であるが、続く
本スキャンＪにおけるそれは、偶数番目の高周波反転パ
ルスのときに印加位相φを１８０°だけ余計に回転させ
た値に設定してある。つまり、φ＝２７０°であって、
この印加位相のパターンは例えば前述したプリスキャン
Ｂのときと同一である。

【０１５６】次いで、コントローラ６は再構成ユニット
１１にステップＳ９４以降の処理を自動的に行うよう指
令する。

【０１５７】この指令に応答して、再構成ユニット１１
は、本スキャンＩ，Ｊにより得られるエコーデータＥi
(n)、Ｅj(n)を用いてｋ空間上の２次元分布の主エコー
成分Ｅmain、副エコー成分Ｅsub を抽出分離する（ステ
ップＳ９４）。この処理は前述してきた手法と同一で、

【数４１】 の式にしたがって行われる。分離抽出された主エコー成
分Ｅmain、副エコー成分Ｅsub の模式図を図３１、３２
に示す。主エコー成分Ｅmainは実線で、副エコー成分Ｅ
sub は点線で示してある。

【０１５８】前記式（１），（２）から明らかなよう
に、主エコー成分Ｅmain、副エコー成分Ｅsub のいずれ
においても、位相の値は、奇数番目においてある値φ
１、偶数番目において別のある値φ２というように２つ
の値φ１、φ２を交互に繰り返す。しかも、主エコー成
分Ｅmainにおける奇数番目のエコー信号の位相値と副エ
コー成分Ｅsub 成分における偶数番目のエコー信号の位
相値は、傾斜磁場の規則性を考慮すると、ほとんど場合
ほぼ同一になる。

【０１５９】この原理にしたがって、再構成ユニット１
１は図３２に示す如く、ｋ空間の中心位置に配置するエ
コー信号が偶数番目のときは奇数番目のエコーデータ
を、反対にそのエコー信号が奇数番目のときは偶数番目
のエコーデータを、主エコー成分Ｅmain、副エコー成分
Ｅsub の各ｋ空間上で、エコーブロック毎にそっくり入
れ替えて新たな主エコー成分Ｅmain′、副エコー成分Ｅ
sub ′を作る（ステップＳ９５）。この結果、新たなｋ
空間上の主エコー成分Ｅmain′、副エコー成分Ｅsub ′
ではエコーブロック毎の位相ずれが殆ど無くなって、図
３２に示すようにｋ空間上の配置データの繋がりがほぼ
完全になる。

【０１６０】しかし、この段階の主エコー成分Ｅmai
n′、副エコー成分Ｅsub ′は共に、ｋ空間上でエコー
ブロック毎に交互に並んでいるだけであるから、図３２
中の振幅値のｋｅ方向の変化グラフＡmain、Ａsub に示
すように、エコーブロック毎に振幅値のばらつきが存在
している。

【０１６１】そこで、再構成ユニット１１はコントロー
ラ６から前述したプリスキャンＡ，Ｂに基づいて作成し
た振幅情報ＥＰmain、ＥＰsub を入力し、これ振幅補正
データとして、主エコー成分Ｅmain′、副エコー成分Ｅ
sub ′の振幅値がｋ空間全体として一定となるようにエ
コーブロック毎に振幅補正を施す（ステップＳ９６）。
この振幅補正により、さらに新たな２組の主エコー成分
Ｅmain”、副エコー成分Ｅsub ”がｋ空間上で作成され
る（図３２参照）。

【０１６２】次いで再構成ユニット１１は、２組の主エ
コー成分Ｅmain”、副エコー成分Ｅsub ”を個別に２次
元フーリエ変換により再構成し、実時間の２枚のＭＲ画
像ＭＲmain、ＭＲsub を夫々作る（ステップＳ９７）。
いずれにｋ空間のエコーデータＥmain”、Ｅsub ”も位
相ずれが殆ど無いので、調整がずれていても再構成され
た画像ＭＲmain、ＭＲsub にゴーストアーチファクトが
発生することは殆ど回避できる。

【０１６３】再構成ユニット１１はさらに、２枚の画像
の同一位置のピクセル信号の位相を互いに揃える処理を
行う（ステップＳ９８）。この処理としては例えば、２
画像の同一位置のピクセル値（複素データ）をＩ_A (r)
，Ｉ_B (r) とすると（図３３参照）、その位相差Δθ
(r) を、

【数４２】Δθ(r) ＝arg{Ｉ_B (r)}−arg{Ｉ_A (r)} を求め、

【数４３】Ｉ_B '(r)＝Ｉ_B (r) ・exp{-jΔθ(r)} の関係を使って、最終画像Ｉ(r) をピクセル毎に、

【数４４】Ｉ(r) ＝｛Ｉ_A (r) ＋Ｉ_B '(r)｝／２ として求める。

【０１６４】あるいは、近似的に

【数４５】 ここで（Ｉ_B (r) ／Ｉ_A (r) ）^＊は（Ｉ_B (r) ／Ｉ
_A (r) ）の共役複素数、として最終画像Ｉ(r) をピクセ
ル毎に求めてもよい。

【０１６５】さらに、

【数４６】 Ｉ(r) ＝｛｜Ｉ_A (r) ｜＋｜Ｉ_B (r) ｜｝／２ として最終画像Ｉ(r) をピクセル毎に求めてもよい。

【０１６６】このように互いに同一ピクセル位置の信号
の位相が揃うと、最後に、両方のＭＲ画像ＭＲmain、Ｍ
Ｒsub をピクセル毎に加算し、１枚の最終画像ＭＲ
_image を作る（ステップＳ９９）。

【０１６７】このように本実施形態によれば、撮像には
最低２回の本スキャン分の時間が必要になるものの、高
周波パルスの印加位相や傾斜磁場波形の補正では原理的
に補正不可能な２次以上の高次の位相誤差分布が存在す
る場合でも、そのような高次成分を含む位相誤差分布を
確実に解消することができ、画質劣化の無い、安定した
高画質のＭＲ画像を提供できる。また、プリスキャンの
位相情報を使用しない方法を採用しているから、撮像対
象の動きなどに耐性を持ったＭＲイメージング法を提供
できる。

【０１６８】なお、上述した第５の実施形態に係る手法
において、振幅補正は状況により不要の場合があるか
ら、プリスキャンＡ，Ｂは必ずしも必須要件ではない。

【０１６９】また高周波反転パルスのフリップ角を下げ
ることで、主エコー成分と副エコー成分の振幅比を変化
させ、ｋ空間上の振幅のばたつきを抑えることもでき
る。かかるフリップ角低下に拠る変形手法の場合、本ス
キャンのパルスシーケンス自体をプリスキャン情報によ
り補正、変更しないので（第５実施形態でもそうであ
る）、プリスキャンに相当するスキャンよりも本スキャ
ンを時間的に先に収集してもよい。

【０１７０】さらに、第５の実施形態及びその変形形態
に係る手法は、ＧＲＡＳＥ法など、ＣＰＭＧパルス系
列、位相反転ＣＰパルス系列にみられる複数の高周波反
転パルスを用いた撮像法に実施できる。

【０１７１】ただし、ＧＲＡＳＥ法の場合、ｋ空間上の
位相エンコード方向に並ぶ複数のエコーブロック（エコ
ー配置区画）の偶数番目又は奇数番目のエコーデータを
単純に変換すればよいという順序性は維持されないこと
がある。その一例を図３４に示す。同図に示す２つのｋ
空間データＥmain，Ｅsub はＧＲＡＳＥ法であって高周
波反転パルス数が奇数（＝Ｎ_RF＝３）、高周波反転パル
ス間のグラジェントエコー数Ｎ_GE＝３でエコーデータＥ
（Ｎ_RF，Ｎ_GE）を各々収集したものである。このように
高周波反転パルス数が奇数のときは、各高周波反転パル
ス間における奇数番目（又は偶数番目）のエコーブロッ
クのエコーデータＥ（Ｎ_RF，Ｎ_GE）を、前述した図２８
のステップＳ２５に相当する処理において変換して位相
が揃った２つの新たなｋ空間データＥmain′，Ｅsub ′
を形成すればよい。このように、位相ずれの規則性を利
用しているので、ＧＲＡＳＥ法のような位相エンコード
波形が各高周波反転パルス間であまり形状的に変化しな
い撮像法ではとくに効果的にアーチファクトを低減させ
ることができる。

【０１７２】

【発明の効果】以上説明したように、本発明のＭＲＩパ
ルスシーケンスの自動補正方法及びＭＲＩ装置によれ
ば、複数の高周波反転パルスの印加位相を基準位相値に
設定したパルスシーケンスの第１のプリスキャン、及び
複数の前記高周波反転パルスの内の偶数番目の高周波反
転パルスの印加位相を基準位相値に対して１８０°の位
相差を有する値に設定したパルスシーケンスの第２のプ
リスキャンを各別に実施し、第１のプリスキャンにより
得た複数のエコー信号から成る第１のエコーデータ群と
第２のプリスキャンにより得た複数のエコー信号から成
る第２のエコーデータ群とに基づいて本スキャンのパル
スシーケンスを事前に補正するため、本スキャンのパル
スシーケンスの補正が、従来とは異なり、十分に分離さ
れた主エコー成分及び副エコー成分に基づいて実施され
ることから、従来よりも高周波反転パルスのフリップ角
の変化に高い耐性を発揮でき、安定したプリスキャンと
なり、強いては補正精度が高く、高品質のＭＲ画像を提
供することができる。

【０１７３】とくに、上記第１、第２のプリスキャンの
パルスシーケンスの位相エンコード用傾斜磁場は、本ス
キャンにおけるｋ空間の中心位置を含むデータラインに
配置するエコーデータ収集時のショットと同一の位相エ
ンコード用傾斜磁場の波形とすることで、位相エンコー
ド方向の傾斜磁場に起因する位相ずれ成分の内、空間的
に一様な位相ずれ成分（０次成分）を補正できる。また
上記第１、第２のプリスキャンのパルスシーケンスの位
相エンコード用傾斜磁場は、本スキャンにおけるｋ空間
上の一部又は全面のエコーデータ収集時のショットと同
一の位相エンコード用傾斜磁場の波形とすることで、位
相エンコード方向の傾斜磁場に起因する位相ずれ成分の
内、位相エンコード方向に１次の空間分布を持つ成分も
補正できる。このように、従来では測定不可能であった
位相エンコード方向の傾斜磁場に起因する位相ずれ成分
を好適に補正でき、補正精度がさらに向上し、高画質の
ＭＲ画像を提供できる。

【０１７４】またとくに、前記第１、第２のプリスキャ
ンを実施する対象を、被検体とは異なるファントムとす
ることで、動きの無い一様なＭＲ信号強度分布を持つ理
想的なファントムを予め用意できるから、被検体内の主
なＭＲ信号源が激しく動くような場合でも、正確な補正
データを得て補正精度を向上させることができる一方、
スキャンの際、被検体に対してプリスキャンを行わなく
ても済むため、被検体の拘束時間が少なくなり、ＭＲ診
断の容易化を図ることもできる。

【０１７５】さらに、本発明に係るＭＲイメージング方
法及びＭＲＩ装置によれば、複数の高周波反転パルスの
印加位相を基準位相値に設定し、かつｋ空間データ収集
用の位相エンコードを施した撮像用パルスシーケンスの
第１のスキャンと、複数の高周波反転パルスの内の偶数
番目の高周波反転パルスの印加位相を基準位相値に対し
て１８０°の位相差を有する値に設定し、かつｋ空間デ
ータ収集用の位相エンコードを施した撮像用パルスシー
ケンスの第２のスキャンとを対象物に各別に実施するス
テップと、第１のスキャンにより得た複数のエコー信号
から成る第１のエコーデータ群と第２のスキャンにより
得た複数のエコー信号から成る第２のエコーデータ群と
を加減算して主エコー成分の集合体からなるｋ空間デー
タとその主エコー成分以外のエコー成分の集合体からな
るｋ空間データとを分離・抽出するステップと、２組の
ｋ空間データ同士で、ｋ空間上の複数のエコー配置区画
の内の偶数番目又は奇数番目のエコー配置区画のエコー
データを互いに交換するステップと、エコーデータ交換
後の２組のｋ空間データを各別に再構成して２つの実空
間画像を作るステップと、２つの実空間画像から１つの
実空間画像を形成するステップと、を含むようにしたの
で、例えば、同じ位相誤差分布を持つ奇数番目の主エコ
ー成分と偶数番目の副エコー成分を用いて画像を再構成
し、一方、同じ位相誤差分布を持つ偶数番目の主エコー
成分と奇数番目の副エコー成分を用いて画像を再構成
し、両方の再構成画像を、例えば位相補正を行って位相
干渉が無いようにした状態で加算することで、最終の実
空間画像を得ることができる。この結果、高周波パルス
の印加位相の補正及び／又は傾斜磁場波形の補正だけで
は原理的に補正不可能であった２次以上の高次の位相誤
差分布が在る場合でも、それらまでも含めたトータルの
位相誤差分布の影響を回避でき、画質劣化が極めて少な
い、安定した高画質のＭＲ画像を提供できる。

【０１７６】このように、本発明によれば、ＣＰＭＧパ
ルス系列、位相反転ＣＰパルス系列などの複数の高周波
反転パルスを使うパルスシーケンスを用いた撮像におい
て、必要な補正情報の収集と補正を従来よりも系統的に
実施して、かかる補正の精度及び安定性を向上させ、又
は、それらの補正をせずとも位相誤差分布の影響を回避
し、空間分解能、ゴーストアーチファクト、画像コント
ラストなどの面で高画質のＭＲ画像を提供することがで
きる。

【図面の簡単な説明】

【図１】本発明の実施形態に係るＭＲＩ装置の概略を示
すブロック図。

【図２】第１の実施形態におけるコントローラの処理を
示す概略フローチャート。

【図３】第１の実施形態の一方のプリスキャンＡを示す
パルスシーケンス。

【図４】第１の実施形態のもう一方のプリスキャンＢを
示すパルスシーケンス。

【図５】主エコー成分の分離抽出を示す説明図。

【図６】主エコー成分、副エコー成分の分離抽出を示す
説明図。

【図７】第１の実施形態における本スキャンのパルスシ
ーケンス。

【図８】副エコー信号の１８０°の位相回転を説明する
ための図。

【図９】第２の実施形態におけるコントローラの処理を
示す概略フローチャート。

【図１０】第２の実施形態の一方のプリスキャンＣを示
すパルスシーケンス。

【図１１】第２の実施形態のもう一方のプリスキャンＤ
を示すパルスシーケンス。

【図１２】主エコー成分、副エコー成分の分離抽出を示
す説明図。

【図１３】第２の実施形態における本スキャンのパルス
シーケンス。

【図１４】第２の実施形態の変形形態におけるコントロ
ーラの処理を示す概略フローチャート。

【図１５】第２の実施形態の変形形態における一方のプ
リスキャンＥを示すパルスシーケンス。

【図１６】第２の実施形態の変形形態におけるもう一方
のプリスキャンＦを示すパルスシーケンス。

【図１７】主エコー成分、副エコー成分の分離抽出を示
す説明図。

【図１８】第２の実施形態の変形形態における本スキャ
ンのパルスシーケンス。

【図１９】第３の実施形態におけるコントローラの処理
を示す概略フローチャート。

【図２０】第３の実施形態の一方のプリスキャンＧを示
すパルスシーケンス。

【図２１】第３の実施形態のもう一方のプリスキャンＨ
を示すパルスシーケンス。

【図２２】主エコー成分、副エコー成分の分離抽出を示
す説明図。

【図２３】第３の実施形態における本スキャンのパルス
シーケンス。

【図２４】第４の実施形態におけるコントローラのプリ
スキャン時の処理を示す概略フローチャート。

【図２５】第４の実施形態におけるコントローラの本ス
キャン時の処理を示す概略フローチャート。

【図２６】第４の実施形態の変形形態におけるプリスキ
ャン時の補正用傾斜磁場係数マトリクスの算出処理を示
す概略フローチャート。

【図２７】第４の実施形態の変形形態におけるプリスキ
ャン時の補正用高周波位相係数ベクトルの算出処理を示
す概略フローチャート。

【図２８】第５の実施形態における画像取得の処理を示
す概略フローチャート。

【図２９】第５の実施形態の一方の本スキャンＩを示す
パルスシーケンス。

【図３０】第５の実施形態のもう一方の本スキャンＪを
示すパルスシーケンス。

【図３１】２次元ｋ空間上の主エコー成分、副エコー成
分の分離抽出を示す説明図。

【図３２】２つの２次元ｋ空間上の主エコー成分、副エ
コー成分から１つの実空間画像までの生成過程を説明す
る図。

【図３３】ピクセル同士の位相合わせを説明するための
図。

【図３４】エコーデータの変換を示すその他の例の模式
図。

【図３５】従来のＦＳＥ法の一例を示すパルスシーケン
ス。

【図３６】従来のＧＲＡＳＥ法の一例を示すパルスシー
ケンス。

【図３７】高周波パルスの印加位相φを説明する模式
図。

【図３８】位相分散ダイヤグラムの一例の図。

【図３９】従来法の一つに係るプリスキャンのパルスシ
ーケンス。

【図４０】従来法の一つに係る本スキャンのパルスシー
ケンス。

【図４１】従来法のプリスキャンから本スキャンまでの
処理の概略フローチャート。

【符号の説明】

１ 磁石 ２ 静磁場電源 ３ｘ〜３ｙ 傾斜磁場コイル ５ シーケンサ ６ コントローラ ７ 高周波コイル ８Ｔ，８Ｒ 送信機、受信機 １１ マルチプレクサ １２ 再構成ユニット １３ 記憶ユニット

Claims (29)

【特許請求の範囲】
1. 【請求項１】 被検体のＭＲ画像を得るための本スキャ
   ンの前に、高周波励起パルス及び複数の高周波反転パル
   スを有しかつその複数の高周波反転パルスの印加位相を
   基準位相値に設定したパルスシーケンスの第１のプリス
   キャン、及び高周波励起パルス及び複数の高周波反転パ
   ルスを有しかつその複数の高周波反転パルスの内の偶数
   番目の高周波反転パルスの印加位相を前記基準位相値に
   対して１８０°の位相差を有する値に設定したパルスシ
   ーケンスの第２のプリスキャンを各別に実施し、 前記第１のプリスキャンにより得た複数のエコー信号か
   ら成る第１のエコーデータ群と前記第２のプリスキャン
   により得た複数のエコー信号から成る第２のエコーデー
   タ群とに基づいて前記本スキャンのパルスシーケンスを
   事前に補正する、ことを特徴としたＭＲＩパルスシーケ
   ンスの自動補正方法。
2. 【請求項２】 前記第１、第２のプリスキャンのパルス
   シーケンス及び前記本スキャンのパルスシーケンスはＣ
   ＰＭＧパルス系列又は位相反転ＣＰパルス系列である請
   求項１記載の自動補正方法。
3. 【請求項３】 前記第１、第２のプリスキャンのパルス
   シーケンス及び前記本スキャンのパルスシーケンスはＧ
   ＲＡＳＥ法又は高速ＳＥ法に従うパルス系列である請求
   項１又は２記載の自動補正方法。
4. 【請求項４】 前記第１、第２のプリスキャンのパルス
   シーケンスの位相エンコード用傾斜磁場は常に零である
   請求項１、２又は３記載の自動補正方法。
5. 【請求項５】 前記自動補正の処理は、 前記第１のエコーデータ群と前記第２のエコーデータ群
   との対応する順番のエコーデータ同士で加算又は減算し
   て主エコー成分から成る主エコー集合体とその主エコー
   成分以外のエコー成分から成る副エコー集合体とを分離
   ・抽出する処理と、 前記主、副エコー集合体に基づいて求めた補正データに
   より前記本スキャンのパルスシーケンスを自動的に補正
   する処理と、を含む請求項４記載の自動補正方法。
6. 【請求項６】 前記自動補正の処理は、 前記第１のエコーデータ群と前記第２のエコーデータ群
   との対応する順番のエコーデータ同士で加算又は減算し
   て主エコー成分から成る主エコー集合体とその主エコー
   成分以外のエコー成分から成る副エコー集合体とを分離
   ・抽出する処理と、 前記主エコー集合体及び副エコー集合体の各集合体の内
   の奇数番目及び偶数番目のエコーデータのエコーピーク
   の位相ずれ及び位置ずれの少なくとも一方を演算する処
   理と、 この演算値から求めた補正データにより前記本スキャン
   のパルスシーケンスを自動的に補正する処理と、を含む
   請求請４記載の自動補正方法。
7. 【請求項７】 前記自動補正の処理は、 前記第１のエコーデータ群と前記第２のエコーデータ群
   との対応する順番のエコーデータ同士で加算又は減算し
   て主エコー成分から成る主エコー集合体とその主エコー
   成分以外のエコー成分から成る副エコー集合体とを分離
   ・抽出する処理と、 前記主エコー集合体及び副エコー集合体の各集合体の内
   の少なくとも一部の０次又は１次の位相分布を比較する
   処理と、 この比較結果に応じて求めた補正データにより前記本ス
   キャンのパルスシーケンスを自動的に補正する処理と、
   を含む請求請４記載の自動補正方法。
8. 【請求項８】 前記補正データは、前記本スキャンのパ
   ルスシーケンスの高周波励起パルスの印加位相の補正
   値、スライス用傾斜磁場の補正値、リード用傾斜磁場の
   補正値の少なくとも一つである請求項５、６又は７記載
   の自動補正方法。
9. 【請求項９】 前記補正データは、前記高周波励起パル
   スの印加位相の補正値又は前記スライス用傾斜磁場の補
   正値、及び前記リード様傾斜磁場の補正値である請求項
   ８記載の自動補正方法。
10. 【請求項１０】 前記第１、第２のプリスキャンのパル
    スシーケンスの位相エンコード用傾斜磁場は、前記本ス
    キャンにおけるｋ空間の中心位置を含むデータラインに
    配置するエコーデータ収集時のショットと同一の位相エ
    ンコード用傾斜磁場の波形である請求項１、２又は３記
    載の自動補正方法。
11. 【請求項１１】 前記自動補正の処理は、 前記第１のエコーデータ群及び前記第２のエコーデータ
    群の内の時間軸上の実効エコー時間に相当するエコーデ
    ータ同士で加算又は減算して実効エコー時間相当のエコ
    ーデータに関する主エコー成分から成る主エコー集合体
    とその主エコー成分以外のエコー成分から成る副エコー
    集合体とを分離・抽出する処理と、 前記主エコー集合体と副エコー集合体との実効エコー時
    間相当のエコーデータのエコーピークの位相ずれ及び位
    置ずれの少なくとも一方を演算する処理と、 この演算値から求めた補正用データで前記本スキャンの
    パルスシーケンスを自動的に補正する処理とから成る請
    求請１０記載の自動補正方法。
12. 【請求項１２】 前記第１、第２のプリスキャンのパル
    スシーケンスの位相エンコード用傾斜磁場は、前記本ス
    キャンにおけるｋ空間上の一部又は全面のエコーデータ
    収集時のショットと同一の位相エンコード用傾斜磁場の
    波形である請求項１、２又は３記載の自動補正方法。
13. 【請求項１３】 前記自動補正の処理は、 前記第１のエコーデータ群及び前記第２のエコーデータ
    群のエコーデータ同士で加算又は減算してｋ空間上の２
    次元の主エコー成分から成る主エコー集合体とその主エ
    コー成分以外のエコー成分から成るｋ空間上の２次元の
    副エコー集合体とを分離・抽出する処理と、 前記主エコー集合体及び副エコー集合体から所望の前記
    主エコー成分のエコーピークの位相ずれと２次元位置ず
    れベクトル及び所望の前記副エコー成分のエコーピーク
    の位相ずれと２次元位置ずれベクトルの少なくとも一方
    を演算する処理と、 この演算値から求めた補正用データで前記本スキャンの
    パルスシーケンスを自動的に補正する処理とから成る請
    求請１２記載の自動補正方法。
14. 【請求項１４】 前記第１、第２のプリスキャンを実施
    する対象は、前記被検体とは異なるファントムである請
    求項１乃至１３のいずれか一項記載の自動補正方法。
15. 【請求項１５】 前記第１、第２のエコーデータ群とに
    基づく補正処理は、この第１、第２のエコーデータ群に
    基づいてイメージングチャンネルと傾斜磁場コイルの物
    理チャンネルとの組み合わせに対して固有の位相誤差に
    関連する情報を求める処理と、本スキャン時にはその位
    相誤差関連の情報と傾斜磁場に関連した撮像条件とから
    本スキャンのパルスシーケンスの補正データを求める処
    理とを含む請求項１、２又は３記載の自動補正方法。
16. 【請求項１６】 前記撮像条件は、本スキャン時に指定
    される撮像断面の方向、分解能、及びスライス厚を含む
    請求項１５記載の自動補正方法。
17. 【請求項１７】 前記高周波反転パルスの印加位相の基
    準位相値は９０°であり、前記高周波反転パルスの１８
    ０°の位相差を有する印加位相値は２７０°である請求
    項１ないし１６の何れか一項記載の自動補正方法。
18. 【請求項１８】 高周波励起パルス及び複数の高周波反
    転パルスを有しかつその複数の高周波反転パルスの印加
    位相を基準位相値に設定したパルスシーケンスの第１の
    スキャン、及び高周波励起パルス及び複数の高周波反転
    パルスを有しかつその複数の前記高周波反転パルスの内
    の偶数番目の高周波反転パルスの印加位相を前記基準位
    相値に対して１８０°の位相差を有する値に設定したパ
    ルスシーケンスの第２のスキャンを対象物に各別に実施
    するステップと、 前記第１のスキャンにより得た複数のエコー信号から成
    る第１のエコーデータ群と前記第２のスキャンにより得
    た複数のエコー信号から成る第２のエコーデータ群とを
    加算又は減算して主エコー成分の主エコー集合体とその
    主エコー成分以外のエコー成分から成る副エコー集合体
    とを分離・抽出するステップと、を含むＭＲイメージン
    グ方法。
19. 【請求項１９】 高周波励起パルス及び複数の高周波反
    転パルスを有し、その複数の高周波反転パルスの印加位
    相を基準位相値に設定し、かつｋ空間データ収集用の位
    相エンコードを施した撮像用パルスシーケンスの第１の
    スキャンと、高周波励起パルス及び複数の高周波反転パ
    ルスを有し、その複数の高周波反転パルスの内の偶数番
    目の高周波反転パルスの印加位相を前記基準位相値に対
    して１８０°の位相差を有する値に設定し、かつｋ空間
    データ収集用の位相エンコードを施した撮像用パルスシ
    ーケンスの第２のスキャンとを対象物に各別に実施する
    ステップと、 前記第１のスキャンにより得た複数のエコー信号から成
    る第１のエコーデータ群と前記第２のスキャンにより得
    た複数のエコー信号から成る第２のエコーデータ群とを
    加減算して主エコー成分の集合体からなるｋ空間データ
    とその主エコー成分以外のエコー成分の集合体からなる
    ｋ空間データとを分離・抽出するステップと、 前記２組のｋ空間データ同士で、ｋ空間上の複数のエコ
    ー配置区画の内の偶数番目又は奇数番目のエコー配置区
    画のエコーデータを互いに交換するステップと、 前記エコーデータ交換後の２組のｋ空間データを各別に
    再構成して２つの実空間画像を作るステップと、 前記２つの実空間画像から１つの実空間画像を形成する
    ステップと、を含むＭＲイメージング方法。
20. 【請求項２０】 前記エコーデータの交換ステップは、
    前記２組のｋ空間データ同士で、ｋ空間上の複数のエコ
    ー配置区画の内の、一つの高周波反転パルスとその次の
    高周波反転パルスの間で生じるエコー信号を配置する偶
    数番目又は奇数番目のエコー配置区画のエコーデータを
    互いに交換するステップである請求項１９記載のＭＲイ
    メージング方法。
21. 【請求項２１】 前記エコーデータの交換ステップは、
    前記２組のｋ空間データ同士で、ｋ空間上の位相エンコ
    ード方向に並ぶ複数のエコー配置区画の内の偶数番目又
    は奇数番目のエコー配置区画のエコーデータを互いに交
    換するステップである請求項１９記載のＭＲイメージン
    グ方法。
22. 【請求項２２】 前記対象物へのプリスキャンとして、
    高周波励起パルス及び複数の高周波反転パルスを有し、
    その複数の高周波反転パルスの印加位相を基準位相値に
    設定したパルスシーケンスの第１のプリスキャンと、高
    周波励起パルス及び複数の高周波反転パルスを有し、そ
    の複数の高周波反転パルスの内の偶数番目の高周波反転
    パルスの印加位相を前記基準位相値に対して１８０°の
    位相差を有する値に設定たパルスシーケンスの第２のプ
    リスキャンとを対象物に各別に実施するステップと、 前記第１のスキャンにより得た複数のエコー信号から成
    る第１のエコーデータ群と前記第２のスキャンにより得
    た複数のエコー信号から成る第２のエコーデータ群とを
    加減算して主エコー成分の集合体の振幅値情報とその主
    エコー成分以外のエコー成分の集合体の振幅値情報とを
    分離・抽出するステップと、 前記エコーデータの交換ステップと前記２つの実空間画
    像の再構成ステップとの間で実施するステップであっ
    て、前記２組の集合体の振幅値情報を使って前記エコー
    データ交換後の２組のｋ空間データの振幅値を補正する
    ステップと、をさらに含む請求項１９、２０又は２１記
    載のＭＲイメージング方法。
23. 【請求項２３】 前記２つの実空間画像から１つの実空
    間画像を形成するステップは、その２つの実空間画像上
    の少なくとも一部領域の同一ピクセル位置同士のピクセ
    ル信号の位相をその絶対値を変えずに揃える処理を含む
    請求項１９乃至２２のいずれか一項に記載のＭＲイメー
    ジング方法。
24. 【請求項２４】 前記位相を揃える処理は、低空間周波
    数のデータのみの位相分布を利用して位相を補正するこ
    とで行う処理である請求項２３記載のＭＲイメージング
    方法。
25. 【請求項２５】 前記位相を揃える処理は、前記２つの
    実空間画像上の全ての領域の同一ピクセル位置同士のピ
    クセル信号の位相をその絶対値を変えずに揃える処理で
    あり、前記２つの実空間画像から１つの実空間画像を形
    成するステップはさらに、位相を揃えた前記２つの実空
    間画像の同一ピクセル位置のピクセル信号同士を加算す
    る処理を含む請求項２３記載のＭＲイメージング方法。
26. 【請求項２６】 前記２つの実空間画像から１つの実空
    間画像を形成するステップは、この２つの実空間画像の
    各ピクセル位置のピクセル信号の絶対値を各画像別に求
    めて２つの絶対値画像を作る処理と、この２つの絶対値
    画像をピクセル毎に加算して前記１つの実空間画像を形
    成する処理とを含む請求項１９、２０又は２１記載のＭ
    Ｒイメージング方法。
27. 【請求項２７】 前記高周波反転パルスの印加位相の基
    準位相値は９０°であり、前記高周波反転パルスの１８
    ０°の位相差を有する印加位相値は２７０°である請求
    項１９ないし２６の何れか一項記載のＭＲイメージング
    方法。
28. 【請求項２８】 被検体のＭＲ画像を得るための本スキ
    ャンの前に、高周波励起パルス及び複数の高周波反転パ
    ルスを有しかつその複数の高周波反転パルスの印加位相
    を基準位相値に設定したパルスシーケンスの第１のプリ
    スキャン、及び高周波励起パルス及び複数の高周波反転
    パルスを有しかつその複数の高周波反転パルスの内の偶
    数番目の高周波反転パルスの印加位相を前記基準位相値
    に対して１８０°の位相差を有する値に設定したパルス
    シーケンスの第２のプリスキャンを各別に実施する手段
    と、 前記第１のプリスキャンにより得た複数のエコー信号か
    ら成る第１のエコーデータ群と前記第２のプリスキャン
    により得た複数のエコー信号から成る第２のエコーデー
    タ群とに基づいて前記本スキャンのパルスシーケンスを
    事前に補正する手段とを備えたことを特徴としたＭＲＩ
    装置。
29. 【請求項２９】 高周波励起パルス及び複数の高周波反
    転パルスを有し、その複数の高周波反転パルスの印加位
    相を基準位相値に設定し、かつｋ空間データ収集用の位
    相エンコードを施した撮像用パルスシーケンスの第１の
    スキャンと、高周波励起パルス及び複数の高周波反転パ
    ルスを有し、その複数の高周波反転パルスの内の偶数番
    目の高周波反転パルスの印加位相を前記基準位相値に対
    して１８０°の位相差を有する値に設定し、かつｋ空間
    データ収集用の位相エンコードを施した撮像用パルスシ
    ーケンスの第２のスキャンとを対象物に各別に実施する
    手段と、 前記第１のスキャンにより得た複数のエコー信号から成
    る第１のエコーデータ群と前記第２のスキャンにより得
    た複数のエコー信号から成る第２のエコーデータ群とを
    加減算して主エコー成分の集合体からなるｋ空間データ
    とその主エコー成分以外のエコー成分の集合体からなる
    ｋ空間データとを分離・抽出する手段と、 前記２組のｋ空間データ同士で、ｋ空間上の複数のエコ
    ー配置区画の内の偶数番目又は奇数番目のエコー配置区
    画のエコーデータを互いに交換する手段と、 前記エコーデータ交換後の２組のｋ空間データを各別に
    再構成して２つの実空間画像を作る手段と、 前記２つの実空間画像から１つの実空間画像を形成する
    手段と、を備えることを特徴としたＭＲＩ装置。