JPH09299345A - 拡散強調イメージング方法およびｍｒｉ装置 - Google Patents

Abstract

(57)【要約】 【課題】 体動によるアーチファクトを十分に低減す
る。 【解決手段】 イメージング用データｅ１（ｍ）とナビ
ゲーション用データｅ２（ｍ）とを収集し（ステップＢ
３）、ナビゲーション用データｅ２（ｍ）の位相φ２
（ｍ）を求め（ステップＢ４）、その位相φ２（ｍ）が
予め求めた複数ビューのナビゲーション用データの平均
位相φ２avから所定の閾値δ以上離れているときは、当
該ビューについてのイメージング用データｅ１（ｍ）と
ナビゲーション用データｅ２（ｍ）を取り直す（ステッ
プＢ５）。そして、ナビゲーション用データｅ２（ｍ）
の位相φ２（ｍ）を基にイメージング用データｅ１
（ｍ）を補正し（ステップＢ６）、補正したイメージン
グ用データＥ１’（ｍ）に基づいて拡散強調画像を作成
する（ステップＢ９，Ｂ１０）。 【効果】 体動分が大き過ぎるか小さ過ぎる確率が高い
データを破棄し、画像作成に使わない。従って、適正に
補正されたイメージング用データのみから画像を作成で
き、良好な拡散強調画像が得られる。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、拡散強調イメージ
ング方法およびＭＲＩ（Magnetic Resonance Imaging）
装置に関する。さらに詳しくは、体動によるアーチファ
クトを低減できる拡散強調イメージング方法およびＭＲ
Ｉ装置に関する。

【０００２】

【従来の技術】イメージング用データの位相の体動分を
補正により除去し、体動によるアーチファクトを低減す
る従来の拡散強調イメージング方法として、特開平４−
３１４４２６号公報および特開平８−１０２３９号公報
に記載の拡散強調イメージング方法が知られている。

【０００３】特開平４−３１４４２６号公報に記載の拡
散強調イメージング方法では、ｋ−空間をラスタ・スキ
ャンのように埋める多数のリード軸方向の直線状軌跡の
それぞれ（これらをビューという）に沿うように位相エ
ンコード勾配およびリード勾配を印加して拡散情報を含
むイメージング用データを収集すると共に、位相巻き戻
し勾配を印加して補正用データを収集する。次いで、ビ
ュー毎に、イメージング用データおよび補正用データを
リード軸方向にフーリエ変換し、それぞれの中間データ
を得る。次に、ビュー毎に、イメージング用データの中
間データの位相を、補正用データの中間データの位相で
補正する。次に、全ビューの補正後のイメージング用デ
ータの中間データを位相エンコード軸方向にフーリエ変
換し、画像用データを生成する。そして、この画像用デ
ータから拡散強調画像を作成する。上記拡散強調イメー
ジング方法において、イメージング用データの位相は、
位相エンコード分と体動分（体動によって生じる位相成
分）の両方を含んでいる。一方、補正用データは、位相
巻き戻し勾配を印加してから収集するため、位相エンコ
ード量が“０”であり、位相エンコード分を含んでおら
ず、体動分のみを含んでいる。そこで、イメージング用
データの位相を補正用データの位相で補正することで、
体動分を除去でき、体動によるアーチファクトを低減で
きる。

【０００４】また、特開平８−１０２３９号公報に記載
の拡散強調イメージング方法では、ｋ−空間の中心部か
らｋ−空間の端部へと螺旋状に広がる多数の螺旋状軌跡
のそれぞれ（これらをビューという）に沿うように位相
エンコード勾配およびリード勾配を印加して拡散情報を
含むイメージング用データを収集する。次に、ビュー毎
に、ｋ−空間の中心部のイメージング用データの位相を
基にｋ−空間の中心部以外のイメージング用データの位
相を補正する。次に、その補正したイメージング用デー
タを２次元フーリエ変換し、画像用データを生成する。
そして、この画像用データから拡散強調画像を作成す
る。上記拡散強調イメージング方法において、ｋ−空間
の中心部以外のイメージング用データは、位相エンコー
ド分と体動分の両方を含んでいる。一方、ｋ−空間の中
心部のイメージング用データは、位相エンコード量が
“０”であるため、位相エンコード分を含んでおらず、
体動分のみを含んでいる。そこで、ｋ−空間の中心部以
外のイメージング用データの位相をｋ−空間の中心部の
イメージング用データで補正することで、体動分を除去
でき、体動によるアーチファクトを低減できる。

【０００５】

【発明が解決しようとする課題】上記特開平４−３１４
４２６号公報に記載の拡散強調イメージング方法におけ
る補正は、補正用データが含む体動分とイメージング用
データが含む体動分が同程度の時は適正になるが、補正
用データが含む体動分よりイメージング用データが含む
体動分が大き過ぎるか小さ過ぎる時には不適正になる。
しかし、従来は、かかる考慮をしていなかったため、不
適正な補正がなされてしまい、体動によるアーチファク
トを十分に低減できないことがある問題点があった。

【０００６】また、上記特開平８−１０２３９号公報に
記載の拡散強調イメージング方法における補正は、ｋ−
空間の中心部のイメージング用データが含む体動分とｋ
−空間の中心部以外のイメージング用データが含む体動
分が同程度の時は適正になるが、ｋ−空間の中心部のイ
メージング用データが含む体動分よりｋ−空間の中心部
以外のイメージング用データが含む体動分が大き過ぎる
か小さ過ぎる時は不適正になる。しかし、従来は、かか
る考慮をしていなかったため、不適正な補正がなされて
しまい、体動によるアーチファクトを十分に低減できな
いことがある問題点があった。

【０００７】そこで、本発明の目的は、不適正な補正が
なされてしまうことを回避し、体動によるアーチファク
トを十分に低減できるようにした拡散強調イメージング
方法およびＭＲＩ装置を提供することにある。

【０００８】

【課題を解決するための手段】第１の観点では、本発明
は、ｋ−空間をラスタ・スキャンのように埋める多数の
リード軸方向の直線状軌跡のそれぞれに沿うように位相
エンコード勾配およびリード勾配を印加して拡散情報を
含むイメージング用データを収集すると共に位相巻き戻
し勾配を印加して補正用データを収集し、各直線状軌跡
に対応するデータ毎に、その補正用データの位相を基に
イメージング用データの位相を補正し、その補正したイ
メージング用データに基づいて拡散情報を反映した画像
を作成する拡散強調イメージング方法において、複数の
直線状軌跡に対応する補正用データの平均位相を取得
し、各直線状軌跡に対応する補正用データの位相と前記
平均位相の差が所定の閾値以上なら、当該直線状軌跡に
ついてはイメージング用データと補正用データを収集し
直すことを特徴とする拡散強調イメージング方法を提供
する。上記第１の観点による拡散強調イメージング方法
では、ある補正用データの位相が、複数の補正用データ
の位相の平均から大きく外れた時は、データの取り直し
を行う。その理由は、かかる時は平均的でない体動があ
ったと考えられるから、補正用データが含む体動分より
イメージング用データが含む体動分が大き過ぎるか小さ
過ぎる確率が高く、そのデータを使わない方がよいから
である。これにより、不適正な補正を行ったイメージン
グ用データを画像の作成に用いることを回避できるよう
になり、体動によるアーチファクトを十分に低減できる
ようになる。なお、閾値は、ＭＲＩ装置や被検体や診断
部位などにより適正な値が異なるため、経験的に値を決
めるのが良い。

【０００９】第２の観点では、本発明は、ｋ−空間の中
心部からｋ−空間の端部へと螺旋状に広がる多数の螺旋
状軌跡のそれぞれに沿うように位相エンコード勾配およ
びリード勾配を印加して拡散情報を含むイメージング用
データを収集し、各螺旋状軌跡に対応するデータ毎に、
そのｋ−空間の中心部のイメージング用データの位相を
基にｋ−空間の中心部以外のイメージング用データの位
相を補正し、その補正したイメージング用データに基づ
いて拡散情報を反映した画像を作成する拡散強調イメー
ジング方法において、複数の螺旋状軌跡に対応するｋ−
空間の中心部のイメージング用データの平均位相を取得
し、各螺旋状軌跡に対応するｋ−空間の中心部のイメー
ジング用データの位相と前記平均位相の差が所定の閾値
以上なら、当該螺旋状軌跡についてはイメージング用デ
ータを収集し直すことを特徴とする拡散強調イメージン
グ方法を提供する。上記第２の観点による拡散強調イメ
ージング方法では、あるｋ−空間の中心部のイメージン
グ用データの位相が、複数のｋ−空間の中心部のイメー
ジング用データの位相の平均から大きく外れた時は、デ
ータの取り直しを行う。その理由は、かかる時は平均的
でない体動があったと考えられるから、ｋ−空間の中心
部のイメージング用データが含む体動分よりｋ−空間の
中心部以外のイメージング用データが含む体動分が大き
過ぎるか小さ過ぎる確率が高く、そのデータを使わない
方がよいからである。これにより、不適正な補正を行っ
たイメージング用データを画像の作成に用いることを回
避できるようになり、体動によるアーチファクトを十分
に低減できるようになる。なお、閾値は、ＭＲＩ装置や
被検体や診断部位などにより適正な値が異なるため、経
験的に値を決めるのが良い。

【００１０】第３の観点では、本発明は、ｋ−空間をラ
スタ・スキャンのように埋める多数のリード軸方向の直
線状軌跡のそれぞれに沿うように位相エンコード勾配お
よびリード勾配を印加して拡散情報を含むイメージング
用データを収集すると共に位相巻き戻し勾配を印加して
補正用データを収集するラスタ・スキャン手段と、複数
の直線状軌跡に対応する補正用データの平均位相を取得
する平均位相取得手段と、補正用データの位相と前記平
均位相の差が所定の閾値以上の直線状軌跡についてはイ
メージング用データと補正用データを収集し直すように
前記ラスタ・スキャン手段を制御するデータ取り直し手
段と、補正用データの位相を基にイメージング用データ
の位相を補正する位相補正手段と、補正したイメージン
グ用データに基づいて拡散情報を反映した画像を作成す
る画像作成手段とを具備したことを特徴とするＭＲＩ装
置を提供する。上記第３の観点によるＭＲＩ装置では、
上記第１の観点による拡散強調イメージング方法を好適
に実施できる。従って、先述のように、体動によるアー
チファクトを十分に低減できるようになる。

【００１１】第４の観点では、本発明は、ｋ−空間の中
心部からｋ−空間の端部へと螺旋状に広がる多数の螺旋
状軌跡のそれぞれに沿うように位相エンコード勾配およ
びリード勾配を印加して拡散情報を含むイメージング用
データを収集するスパイラル・スキャン手段と、複数の
螺旋状軌跡に対応するｋ−空間の中心部のイメージング
用データの平均位相を取得する平均位相取得手段と、ｋ
−空間の中心部のイメージング用データの位相と前記平
均位相の差が所定の閾値以上の螺旋状軌跡についてはイ
メージング用データを収集し直すように前記スパイラル
・スキャン手段を制御するデータ取り直し手段と、ｋ−
空間の中心部のイメージング用データの位相を基にｋ−
空間の中心部以外のイメージング用データの位相を補正
する位相補正手段と、補正したイメージング用データに
基づいて拡散情報を反映した画像を作成する画像作成手
段とを具備したＭＲＩ装置を提供する。上記第４の観点
によるＭＲＩ装置では、上記第２の観点による拡散強調
イメージング方法を好適に実施できる。従って、先述の
ように、体動によるアーチファクトを十分に低減できる
ようになる。

【００１２】

【発明の実施の形態】以下、図に示す実施の形態により
本発明をさらに詳細に説明する。なお、これにより本発
明が限定されるものではない。

【００１３】−第１の実施形態− 図１は、本発明の第１の実施形態にかかるＭＲＩ装置の
構成図である。このＭＲＩ装置１００において、マグネ
ットアセンブリ１は、内部に被検体を挿入するための空
間部分（孔）を有し、この空間部分を取りまくようにし
て、被検体に一定の主磁場を印加する主磁場コイルと、
勾配磁場を発生するための勾配磁場コイル（勾配磁場コ
イルはｘ軸，ｙ軸，ｚ軸の各コイルを備えており、これ
らの組み合わせによりスライス軸，位相エンコード軸，
リード軸が決まる）と、被検体内の原子核のスピンを励
起するためのＲＦパルスを送信する送信コイルと、被検
体からのＮＭＲ信号を受信する受信コイル等が配置され
ている。主磁場コイル，勾配磁場コイル，送信コイルお
よび受信コイルは、それぞれ主磁場電源２，勾配磁場駆
動回路３，ＲＦ電力増幅器４および前置増幅器５に接続
されている。

【００１４】計算機７は、パルスシーケンスを作成し、
シーケンス記憶回路８に渡す。シーケンス記憶回路８
は、パルスシーケンスを記憶し、そのパルスシーケンス
に基づいて勾配磁場駆動回路３を操作し、マグネットア
センブリ１の勾配磁場コイルから勾配磁場を発生させる
と共に、ゲート変調回路９を操作し、ＲＦ発振回路１０
の搬送波出力信号を所定タイミング・所定包絡線形状の
パルス状信号に変調し、それをＲＦパルスとしてＲＦ電
力増幅器４に加え、ＲＦ電力増幅器４でパワー増幅した
後、前記マグネットアセンブリ１の送信コイルに印加す
る。

【００１５】前置増幅器５は、マグネットアセンブリ１
の受信コイルで受信したＮＭＲ信号を増幅し、位相検波
器１２に入力する。位相検波器１２は、ＲＦ発振回路１
０の搬送波出力信号を参照信号とし、ＮＭＲ信号を位相
検波して、Ａ／Ｄ変換器１１に与える。Ａ／Ｄ変換器１
１は、アナログ信号のＮＭＲ信号をディジタル信号のＭ
Ｒデータに変換し、計算機７に入力する。

【００１６】計算機７は、Ａ／Ｄ変換器１１からＭＲデ
ータを読み込み、位相補正演算や画像再構成演算を行
い、画像を作成する。この画像は、表示装置６にて表示
される。また、計算機７は、操作卓１３から入力された
情報を受け取るなどの全体的な制御を受け持つ。

【００１７】図２は、上記ＭＲＩ装置１００における平
均位相取得処理のフローチャートである。ステップＶ１
では、ビュー番号用カウンタｍを“１”に初期化する。
ステップＶ２では、図３に示すパルスシーケンスＱＲに
より補正用データｅ２（ｍ）を収集する。すなわち、こ
のパルスシーケンスＱＲでは、９０°のＲＦパルスＲ１
を印加すると共にスライス軸にスライス勾配Ｓ１を印加
する。次に、位相エンコード軸に位相エンコード勾配Ｐ
Ｅ（ｍ）を印加する。次に、拡散強調用の強力なＭＰ
（Motion Probing）勾配ＭＰＧ１を任意の勾配軸に印加
する。次に、１８０°のＲＦパルスＲ２を印加すると共
にスライス軸にスライス勾配Ｓ２を印加する。次に、拡
散強調用のＭＰ勾配ＭＰＧ２を任意の勾配軸に印加す
る。次に、リード軸にリード勾配ＲＤ１を印加しながら
エコーecho１をサンプリングして、イメージング用デー
タｅ１（ｍ）を収集する。次に、位相エンコード軸に前
記位相エンコード勾配ＰＥ（ｍ）と同じ大きさの位相巻
き戻しエンコード勾配ＰＲ（ｍ）を印加する。次に、リ
ード軸にリード勾配ＲＤ２を印加しながらエコーecho２
をサンプリングして、補正用データｅ２（ｍ）を収集す
る。図４の（ａ）に、第１ビュー（ｍ＝１）のイメージ
ング用データｅ１（１）および補正用データｅ２（１）
の直線状軌跡（ラスタ・トラジェクトリ）を例示する。
また、図４の（ｂ）に、第２ビュー（ｍ＝２）のイメー
ジング用データｅ１（２）および補正用データｅ２
（２）の直線状軌跡を示す。なお、上記平均位相取得処
理では、イメージング用データｅ１（ｍ）の収集を省い
てもよい。

【００１８】図２に戻って、ステップＶ３では、前記補
正用データｅ２（ｍ）をリード軸方向にフーリエ変換
し、中間データＥ２（ｍ）を求め、さらに、その位相φ
２（ｍ）を求める。ステップＶ４，Ｖ５では、ｍ＝２〜
Ｎ（Ｎは、例えば８〜１２）について上記ステップＶ２
およびステップＶ３を反復実行する。ステップＶ６で
は、位相φ２（１）〜φ２（Ｎ）の平均位相φ２avを算
出する。

【００１９】図５は、上記ＭＲＩ装置１００における拡
散強調イメージング処理のフローチャートである。ステ
ップＢ１では、操作者が閾値δ（δは、例えば１０°）
を設定する。ステップＢ２では、ビュー番号用カウンタ
ｍを“１”に初期化する。ステップＢ３では、図３に示
すパルスシーケンスＱＲによりイメージング用データｅ
１（ｍ）と補正用データｅ２（ｍ）を収集する。すなわ
ち、このパルスシーケンスＱＲでは、９０°のＲＦパル
スＲ１を印加すると共にスライス軸にスライス勾配Ｓ１
を印加する。次に、位相エンコード軸に位相エンコード
勾配ＰＥ（ｍ）を印加する。次に、拡散強調用のＭＰ勾
配ＭＰＧ１を任意の勾配軸に印加する。次に、１８０°
のＲＦパルスＲ２を印加すると共にスライス軸にスライ
ス勾配Ｓ２を印加する。次に、拡散強調用のＭＰ勾配Ｍ
ＰＧ２を任意の勾配軸に印加する。次に、リード軸にリ
ード勾配ＲＤ１を印加しながらエコーecho１をサンプリ
ングして、イメージング用データｅ１（ｍ）を収集す
る。次に、位相エンコード軸に前記位相エンコード勾配
ＰＥ（ｍ）と同じ大きさの位相巻き戻しエンコード勾配
ＰＲ（ｍ）を印加する。次に、リード軸にリード勾配Ｒ
Ｄ２を印加しながらエコーecho２をサンプリングして、
補正用データｅ２（ｍ）を収集する。

【００２０】図５に戻り、ステップＢ４では、前記補正
用データｅ２（ｍ）をリード軸方向にフーリエ変換して
中間データＥ２（ｍ）を求め、その中間データＥ２
（ｍ）の位相φ２（ｍ）を求める。ステップＢ５では、
｜φ２（ｍ）−φ２av｜＜δかを判定する。｜φ２
（ｍ）−φ２av｜＜δならステップＢ６に進む。｜φ２
（ｍ）−φ２av｜＜δでないなら上記ステップＢ３に戻
り、第ｍビューのイメージング用データｅ１（ｍ）と補
正用データｅ２（ｍ）を取り直す。

【００２１】ステップＢ６では、前記イメージング用デ
ータｅ１（ｍ）をリード軸方向にフーリエ変換して中間
データＥ１（ｍ）を求め、それに対して、 exp｛−ｉφ２（ｍ）｝ を乗算し、第ｍビューの位相補正後の中間データＥ１’
（ｍ）を得る。ステップＢ７，Ｂ８では、ｍ＝２〜Ｍ
（Ｍは、例えば１２６）について上記ステップＢ３〜ス
テップＢ６を反復実行する。ステップＢ９では、得られ
た位相補正後の中間データＥ１’（１）〜Ｅ１’（Ｍ）
を位相エンコード軸方向にフーリエ変換し、画像データ
を生成する。ステップＢ１０では、前記画像データから
拡散強調画像を作成する。

【００２２】図６は、上記拡散強調イメージング処理
（図５）におけるｍの進行の説明図である。ｍ＝６，８
では、φ２（ｍ）と平均位相φ２avの差が閾値δ内にな
るまでにデータを２回取り直している。

【００２３】上記第１の実施形態のＭＲＩ装置１００に
よれば、補正用データｅ２が含む体動分よりイメージン
グ用データｅ１が含む体動分が大き過ぎるか小さ過ぎる
確率が高いデータを使わないから、常に適正に補正され
たイメージング用データのみを画像の作成に用いること
が出来るようになり、体動によるアーチファクトを十分
に低減できるようになる。

【００２４】なお、図５の拡散強調イメージング処理で
は、｜φ２（ｍ）−φ２av｜＜δとなるまでｍをインク
リメントせずにデータを取り直しているが、データを取
り直さずにｍをインクリメントしながらデータを収集
し、その後、｜φ２（ｍ）−φ２av｜＜δとならなかっ
たｍについてデータを取り直すようにしてもよい。さら
に、この方式の拡張として、図２の平均位相取得処理
で、Ｎ＝Ｍとしてイメージング用データｅ１（ｍ）およ
び補正用データｅ２（ｍ）を収集しておき、補正用デー
タｅ２（ｍ）の全部または一部から平均位相φ２avを求
め、｜φ２（ｍ）−φ２av｜＜δであるｍについてのデ
ータは画像作成に使い、｜φ２（ｍ）−φ２av｜＜δで
ないｍについてはデータを取り直すようにしてもよい。

【００２５】また、図５の拡散強調イメージング処理で
は、データの取り直しの回数を制限していないが、同一
ビューでの取り直し回数を制限し（例えば１０回）、そ
の回数だけデータの取り直しを行っても｜φ２（ｍ）−
φ２av｜＜δにならなかった場合は、最後に得られたデ
ータを画像作成に使うようにしてもよい。

【００２６】−第２の実施形態− 本発明の第２の実施形態のＭＲＩ装置の構成は、図１と
同様である。図７，図８は、本発明の第２の実施形態の
ＭＲＩ装置における拡散強調イメージング処理のフロー
チャートである。ステップＣ１では、操作者が閾値δ
（δは、例えば１０°）を設定する。ステップＣ２で
は、ビュー番号用カウンタｍを“１”に初期化する。ス
テップＣ３では、図３に示すパルスシーケンスＱＲによ
りイメージング用データｅ１（ｍ）と補正用データｅ２
（ｍ）を収集する。すなわち、このパルスシーケンスＱ
Ｒでは、９０°のＲＦパルスＲ１を印加すると共にスラ
イス軸にスライス勾配Ｓ１を印加する。次に、位相エン
コード軸に位相エンコード勾配ＰＥ（ｍ）を印加する。
次に、拡散強調用のＭＰ勾配ＭＰＧ１を任意の勾配軸に
印加する。次に、１８０°のＲＦパルスＲ２を印加する
と共にスライス軸にスライス勾配Ｓ２を印加する。次
に、拡散強調用のＭＰ勾配ＭＰＧ２を任意の勾配軸に印
加する。次に、リード軸にリード勾配ＲＤ１を印加しな
がらエコーecho１をサンプリングして、イメージング用
データｅ１（ｍ）を収集する。次に、位相エンコード軸
に前記位相エンコード勾配ＰＥ（ｍ）と同じ大きさの位
相巻き戻しエンコード勾配ＰＲ（ｍ）を印加する。次
に、リード軸にリード勾配ＲＤ２を印加しながらエコー
echo２をサンプリングして、補正用データｅ２（ｍ）を
収集する。

【００２７】図７に戻り、ステップＣ４では、前記イメ
ージング用データｅ１（ｍ）をリード軸方向にフーリエ
変換して中間データＥ１（ｍ）を求め、それに対して、 exp｛−ｉφ２（ｍ）｝ を乗算し、第ｍビューの位相補正後の中間データＥ１’
（ｍ）を得る。ステップＣ５，Ｃ６では、ｍ＝２〜Ｎ
（Ｎは、例えば８〜１２）について上記ステップＣ３を
反復実行する。

【００２８】ステップＣ７では、前記補正用データｅ２
（１）〜ｅ２（Ｎ）をリード軸方向にフーリエ変換して
中間データＥ２（１）〜Ｅ２（Ｎ）を求め、それらの中
間データＥ２（１）〜Ｅ２（Ｎ）の位相φ２（１）〜φ
２（Ｎ）を求め、さらに、それら位相φ２（１）〜φ２
（Ｎ）の平均位相φ２avを算出する。

【００２９】図８へ進み、ステップＣ８では、図３に示
すパルスシーケンスＱＲによりイメージング用データｅ
１（ｍ）と補正用データｅ２（ｍ）を収集する。ここで
は、Ｎ＜ｍである。ステップＣ９では、前記補正用デー
タｅ２（ｍ）をリード軸方向にフーリエ変換して中間デ
ータＥ２（ｍ）を求め、その中間データＥ２（ｍ）の位
相φ２（ｍ）を求める。ステップＣ１０では、｜φ２
（ｍ）−φ２av｜＜δかを判定する。｜φ２（ｍ）−φ
２av｜＜δならステップＣ１１に進む。｜φ２（ｍ）−
φ２av｜＜δでないなら上記ステップＣ８に戻り、第ｍ
ビューのイメージング用データｅ１（ｍ）と補正用デー
タｅ２（ｍ）を取り直す。

【００３０】ステップＣ１１では、前記イメージング用
データｅ１（ｍ）をリード軸方向にフーリエ変換して中
間データＥ１（ｍ）を求め、それに対して、 exp｛−ｉφ２（ｍ）｝ を乗算し、第ｍビューの位相補正後の中間データＥ１’
（ｍ）を得る。ステップＣ１２，Ｃ１３では、ｍ＝Ｎ＋
２〜Ｍ（Ｍは、例えば１２６）について上記ステップＣ
８〜ステップＣ１１を反復実行する。ステップＣ１４で
は、得られた位相補正後の中間Ｅ１’（１）〜Ｅ１’
（Ｍ）を位相エンコード軸方向にフーリエ変換し、画像
データを生成する。ステップＣ１５では、前記画像デー
タから拡散強調画像を作成する。

【００３１】上記第２の実施形態のＭＲＩ装置でも、体
動によるアーチファクトを十分に低減できるようにな
る。また、図２に示す平均位相取得処理のプリスキャン
を省くことが出来る。

【００３２】−第３の実施形態− 本発明の第３の実施形態のＭＲＩ装置の構成は図１と同
様である。図９は、本発明の第３の実施形態のＭＲＩ装
置における平均位相取得処理のフローチャートである。
ステップＦ１では、ビュー番号用カウンタｍを“１”に
初期化する。ステップＦ２では、図１０に示すパルスシ
ーケンスＱＳによりイメージング用データＤ（ｍ，ｋ
ｘ，ｋｙ）を収集する。このパルスシーケンスＱＳで
は、９０°のＲＦパルスＲ１を印加すると共にスライス
軸にスライス勾配Ｓ１を印加する。次に、拡散強調用の
ＭＰ勾配ＭＰＧ１を任意の勾配軸に印加する。次に、１
８０°のＲＦパルスＲ２を印加すると共にスライス軸に
スライス勾配Ｓ２を印加する。次に、拡散強調用のＭＰ
勾配ＭＰＧ２を任意の勾配軸に印加する。次に、図１１
に示すようにｋ−空間Ｓの中心部から端部へと螺旋状に
広がる螺旋状軌跡（スパイラル・トラジェクトリ）を形
成するように位相エンコード勾配ＲＤ３とリード勾配Ｒ
Ｄ４を印加しながら、エコーechoをサンプリングして、
イメージング用データＤ（ｍ，ｋｘ，ｋｙ）を収集す
る。図１１の（ａ）に、第１ビュー（ｍ＝１）のイメー
ジング用データＤ（１，ｋｘ，ｋｙ）のスパイラル・ト
ラジェクトリを示す。Ｄ（１，０，０）は、ｋ−空間Ｓ
の中心部のイメージング用データである。図１１の
（ｂ）に、第２ビュー（ｍ＝２）のイメージング用デー
タＤ（２，ｋｘ，ｋｙ）のスパイラル・トラジェクトリ
を示す。なお、Ｄ（２，０，０）は、ｋ−空間Ｓの中心
部のイメージング用データである。

【００３３】図９に戻って、ステップＦ３では、ｋ−空
間の中心部のイメージング用データＤ（ｍ，０，０）の
位相φ（ｍ）を次式により求める。 φ（ｍ）＝arg｛Ｄ（ｍ，０，０）｝ ステップＦ４，Ｆ５では、ｍ＝２〜Ｎ（Ｎは、例えば４
０）について上記ステップＦ２およびステップＦ３を反
復実行する。ステップＦ６では、位相φ（１）〜φ
（Ｎ）の平均位相φ２avを求める。

【００３４】図１２は、第３の実施形態のＭＲＩ装置に
おける拡散強調イメージング処理のフローチャートであ
る。ステップＨ１では、操作者が閾値δ（δは、例えば
１０°）を設定する。ステップＨ２では、ビュー番号用
カウンタｍを“１”に初期化する。ステップＨ３では、
図１０に示すパルスシーケンスＱＳによりイメージング
用データＤ（ｍ，ｋｘ，ｋｙ）を収集する。ステップＨ
４では、ｋ−空間の中心部のイメージング用データＤ
（ｍ，０，０）の位相φ（ｍ）を次式により求める。 φ（ｍ）＝arg｛Ｄ（ｍ，０，０）｝ ステップＨ５では、｜φ（ｍ）−φ２av｜＜δか否かを
判定する。｜φ（ｍ）−φ２av｜＜δならステップＨ６
に進む。｜φ（ｍ）−φ２av｜＜δでないなら上記ステ
ップＨ３に戻り、第ｍビューのイメージング用データＤ
（ｍ，ｋｘ，ｋｙ）を取り直す。

【００３５】ステップＨ６では、イメージング用データ
Ｄ（ｍ，ｋｘ，ｋｙ）に対して、 exp｛−ｉφ（ｍ）｝ を乗算し、これを第ｍビューの位相補正後のイメージン
グ用データＤ’（ｍ，ｋｘ，ｋｙ）とする。ステップＨ
７，Ｈ８では、ｍ＝２〜Ｍ（Ｍは、例えば１６）につい
て上記ステップＨ３〜ステップＨ６を反復実行する。ス
テップＨ９では、得られた位相補正後のイメージング用
データＤ’（１，ｋｘ，ｋｙ）〜Ｄ’（Ｍ，ｋｘ，ｋ
ｙ）を２次元フーリエ変換し、画像用データを生成す
る。ステップＨ１０では、前記画像データから拡散強調
画像を作成する。

【００３６】上記第３の実施形態のＭＲＩ装置によれ
ば、ｋ−空間の中心部のイメージング用データが含む体
動分よりｋ−空間の中心部以外のイメージング用データ
が含む体動分が大き過ぎるか小さ過ぎる確率が高いデー
タを使わないから、常に適正に補正されたイメージング
用データのみを画像の作成に用いることが出来るように
なり、体動によるアーチファクトを十分に低減できるよ
うになる。

【００３７】なお、図１２の拡散強調イメージング処理
では、｜φ（ｍ）−φav｜＜δとなるまでｍをインクリ
メントせずにデータを取り直しているが、データを取り
直さずにｍをインクリメントしながらデータを収集し、
その後、｜φ（ｍ）−φav｜＜δとならなかったｍにつ
いてデータを取り直すようにしてもよい。さらに、この
方式の拡張として、図９の平均位相取得処理で、Ｎ＝Ｍ
としてイメージング用データＤ（ｍ）を収集しておき、
イメージング用データＤ（ｍ）の全部または一部から平
均位相φavを求め、｜φ（ｍ）−φav｜＜δであるｍに
ついてのデータは画像作成に使い、｜φ（ｍ）−φav｜
＜δでないｍについてはデータを取り直すようにしても
よい。

【００３８】また、図１２の拡散強調イメージング処理
では、データの取り直しの回数を制限していないが、同
一ビューでの取り直し回数を制限し（例えば１０回）、
その回数だけデータの取りを行っても｜φ（ｍ）−φav
｜＜δにならなかった場合は、最後に得られたデータを
画像作成に使うようにしてもよい。

【００３９】さらに、前記第２の実施形態と同様に、ｍ
＝１〜Ｎのイメージング用データＤ（ｍ）を収集し、そ
れらから平均位相φavを求め、ｍ＝Ｎ＋１〜Ｍのイメー
ジング用データＤ（ｍ）については｜φ（ｍ）−φav｜
＜δにならなかった場合にデータを取り直すようにして
もよい。

【００４０】

【発明の効果】本発明の拡散強調イメージング方法およ
びＭＲＩ装置によれば、イメージング用データが含む体
動分が大き過ぎるか小さ過ぎる確率が高いデータを破棄
し、画像作成に使わないから、常に適正に補正されたイ
メージング用データのみから画像を作成できるようにな
り、体動によるアーチファクトを十分に低減した良好な
拡散強調画像が得られるようになる。

【図面の簡単な説明】

【図１】本発明の第１の実施形態にかかるＭＲＩ装置の
構成図である。

【図２】本発明の第１の実施形態に係る平均位相取得処
理のフローチャートである。

【図３】本発明の第１の実施形態にかかるパルスシーケ
ンスの例示図である。

【図４】ラスタ・トラジェクトリの説明図である。

【図５】本発明の第１の実施形態にかかる拡散強調イメ
ージング処理のフローチャートである。

【図６】データの取り直しの説明図である。

【図７】本発明の第２の実施形態にかかる拡散強調イメ
ージング処理の前半のフローチャートである。

【図８】本発明の第２の実施形態にかかる拡散強調イメ
ージング処理の後半のフローチャートである。

【図９】本発明の第３の実施形態に係る平均位相取得処
理のフローチャートである。

【図１０】本発明の第３の実施形態にかかるパルスシー
ケンスの例示図である。

【図１１】スパイラル・トラジェクトリの説明図であ
る。

【図１２】本発明の第３の実施形態にかかる拡散強調イ
メージング処理のフローチャートである。

【符号の説明】

１００ ＭＲＩ装置 １ マグネットアセンブリ ３ 勾配磁場駆動回路 ７ 計算機 ８ シーケンス記憶回路 ＱＲ，ＱＳ パルスシーケンス ｅ１，Ｄ イメージング用データ ｅ２ 補正用データ ＭＰＧ１，ＭＰＧ２ ＭＰ勾配

Claims (4)

【特許請求の範囲】
1. 【請求項１】 ｋ−空間をラスタ・スキャンのように埋
   める多数の周波数軸方向の直線状軌跡のそれぞれに沿う
   ように位相エンコード勾配およびリード勾配を印加して
   拡散情報を含むイメージング用データを収集すると共に
   位相巻き戻し勾配を印加して補正用データを収集し、各
   直線状軌跡に対応するデータ毎に、その補正用データの
   位相を基にイメージング用データの位相を補正し、その
   補正したイメージング用データに基づいて拡散情報を反
   映した画像を作成する拡散強調イメージング方法におい
   て、 複数の直線状軌跡に対応する補正用データの平均位相を
   取得し、各直線状軌跡に対応する補正用データの位相と
   前記平均位相の差が所定の閾値以上なら、当該直線状軌
   跡についてはイメージング用データと補正用データを収
   集し直すことを特徴とする拡散強調イメージング方法。
2. 【請求項２】 ｋ−空間の中心部からｋ−空間の端部へ
   と螺旋状に広がる多数の螺旋状軌跡のそれぞれに沿うよ
   うに位相エンコード勾配およびリード勾配を印加して拡
   散情報を含むイメージング用データを収集し、各螺旋状
   軌跡に対応するデータ毎に、そのｋ−空間の中心部のイ
   メージング用データの位相を基にｋ−空間の中心部以外
   のイメージング用データの位相を補正し、その補正した
   イメージング用データに基づいて拡散情報を反映した画
   像を作成する拡散強調イメージング方法において、 複数の螺旋状軌跡に対応するｋ−空間の中心部のイメー
   ジング用データの平均位相を取得し、各螺旋状軌跡に対
   応するｋ−空間の中心部のイメージング用データの位相
   と前記平均位相の差が所定の閾値以上なら、当該螺旋状
   軌跡についてはイメージング用データを収集し直すこと
   を特徴とする拡散強調イメージング方法。
3. 【請求項３】 ｋ−空間をラスタ・スキャンのように埋
   める多数のリード軸方向の直線状軌跡のそれぞれに沿う
   ように位相エンコード勾配およびリード勾配を印加して
   拡散情報を含むイメージング用データを収集すると共に
   位相巻き戻し勾配を印加して補正用データを収集するラ
   スタ・スキャン手段と、複数の直線状軌跡に対応する補
   正用データの平均位相を取得する平均位相取得手段と、
   補正用データの位相と前記平均位相の差が所定の閾値以
   上の直線状軌跡についてはイメージング用データと補正
   用データを収集し直すように前記ラスタ・スキャン手段
   を制御するデータ取り直し手段と、補正用データの位相
   を基にイメージング用データの位相を補正する位相補正
   手段と、補正したイメージング用データに基づいて拡散
   情報を反映した画像を作成する画像作成手段とを具備し
   たことを特徴とするＭＲＩ装置。
4. 【請求項４】 ｋ−空間の中心部からｋ−空間の端部へ
   と螺旋状に広がる多数の螺旋状軌跡のそれぞれに沿うよ
   うに位相エンコード勾配およびリード勾配を印加して拡
   散情報を含むイメージング用データを収集するスパイラ
   ル・スキャン手段と、複数の螺旋状軌跡に対応するｋ−
   空間の中心部のイメージング用データの平均位相を取得
   する平均位相取得手段と、ｋ−空間の中心部のイメージ
   ング用データの位相と前記平均位相の差が所定の閾値以
   上の螺旋状軌跡についてはイメージング用データを収集
   し直すように前記スパイラル・スキャン手段を制御する
   データ取り直し手段と、ｋ−空間の中心部のイメージン
   グ用データの位相を基にｋ−空間の中心部以外のイメー
   ジング用データの位相を補正する位相補正手段と、補正
   したイメージング用データに基づいて拡散情報を反映し
   た画像を作成する画像作成手段とを具備したことを特徴
   とするＭＲＩ装置。