JPH06296600A - 二次元ｐｃ法を用いたｍｒアンギオグラフィー方法及び装置 - Google Patents
二次元ｐｃ法を用いたｍｒアンギオグラフィー方法及び装置Info
- Publication number
- JPH06296600A JPH06296600A JP5089912A JP8991293A JPH06296600A JP H06296600 A JPH06296600 A JP H06296600A JP 5089912 A JP5089912 A JP 5089912A JP 8991293 A JP8991293 A JP 8991293A JP H06296600 A JPH06296600 A JP H06296600A
- Authority
- JP
- Japan
- Prior art keywords
- pulse
- angiography
- dimensional
- gradient
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Abstract
(57)【要約】
【目的】 血流の流速と方向をイメージングすることが
可能な二次元ＰＣ法を用いたＭＲアンギオグラフィー方
法及び装置を提供する。 【構成】 ＭＴ法によるポリノミアルパルスＰＰおよび
スポイラ勾配ＳＰを印加して静止部の磁化強度を減少さ
せてから、ＲＦパルスαおよびスライス勾配Ｓを印加し
て所定のスライスを選択励起し、異なる２種類のバイポ
ーラ勾配パルスＧ１，Ｇ２を印加して２種類のデータを
収集し、それら２種類のデータの持つ位相成分の違いか
ら血流のイメージを取得する。プロジェクションディフ
ェーズ勾配は印加しない。 【効果】 ボクセルがスライス方向に長い場合でも流れ
の方向が判り且つ消え残り誤差を小さくすることが出来
る。
可能な二次元ＰＣ法を用いたＭＲアンギオグラフィー方
法及び装置を提供する。 【構成】 ＭＴ法によるポリノミアルパルスＰＰおよび
スポイラ勾配ＳＰを印加して静止部の磁化強度を減少さ
せてから、ＲＦパルスαおよびスライス勾配Ｓを印加し
て所定のスライスを選択励起し、異なる２種類のバイポ
ーラ勾配パルスＧ１，Ｇ２を印加して２種類のデータを
収集し、それら２種類のデータの持つ位相成分の違いか
ら血流のイメージを取得する。プロジェクションディフ
ェーズ勾配は印加しない。 【効果】 ボクセルがスライス方向に長い場合でも流れ
の方向が判り且つ消え残り誤差を小さくすることが出来
る。
Description
【０００１】
【産業上の利用分野】この発明は、二次元ＰＣ（Ｐhase
Ｃontrast）法を用いたＭＲアンギオグラフィー方法及
び装置に関し、さらに詳しくは、ボクセルがスライス方
向に長い場合でも流れの方向が判り且つ消え残り誤差を
小さくすることが出来る二次元ＰＣ法を用いたＭＲアン
ギオグラフィー方法及び装置に関する。
Ｃontrast）法を用いたＭＲアンギオグラフィー方法及
び装置に関し、さらに詳しくは、ボクセルがスライス方
向に長い場合でも流れの方向が判り且つ消え残り誤差を
小さくすることが出来る二次元ＰＣ法を用いたＭＲアン
ギオグラフィー方法及び装置に関する。
【０００２】
【従来の技術】図６は、二次元ＰＣ法を用いたＭＲアン
ギオグラフィー方法に用いられる従来のパルスシーケン
スである。このパルスシーケンスＢでは、ＲＦパルスα
およびスライス勾配Ｓを印加してスライスを選択励起
し、実線で示すバイポーラ勾配パルスＧ１と、リード勾
配Ｒと、プロジェクションディフェーズ勾配ＰＤと、位
相エンコード勾配ＰＥを印加して、第１のデータを収集
する。次に、実線で示すバイポーラ勾配パルスＧ１を点
線で示すバイポーラ勾配パルスＧ２に替えて、第２のデ
ータを収集する。そして、前記第１のデータと前記第２
のデータの差をとって画像化するか、または、前記第１
のデータによる画像と前記第２のデータによる画像の差
をとって、流動部のイメージを取得する。この流動部の
イメージでは、流速が輝度で表現されている。
ギオグラフィー方法に用いられる従来のパルスシーケン
スである。このパルスシーケンスＢでは、ＲＦパルスα
およびスライス勾配Ｓを印加してスライスを選択励起
し、実線で示すバイポーラ勾配パルスＧ１と、リード勾
配Ｒと、プロジェクションディフェーズ勾配ＰＤと、位
相エンコード勾配ＰＥを印加して、第１のデータを収集
する。次に、実線で示すバイポーラ勾配パルスＧ１を点
線で示すバイポーラ勾配パルスＧ２に替えて、第２のデ
ータを収集する。そして、前記第１のデータと前記第２
のデータの差をとって画像化するか、または、前記第１
のデータによる画像と前記第２のデータによる画像の差
をとって、流動部のイメージを取得する。この流動部の
イメージでは、流速が輝度で表現されている。
【０００３】図７に、二次元ＰＣ法によるボクセルＶを
示す。ｙは位相エンコード方向、ｘは周波数エンコード
方向、ｚはスライス方向である。二次元ＰＣ法では、ス
ライス方向には空間分解能がなく（それ故、スライス方
向に投影され、二次元となる）、ボクセルＶはスライス
方向に長い直方体になる。しかし、スライス方向に長く
なると、各ボクセルＶに占める血流Ｋの割合が小さくな
り、大部分は静止部になる。そこで、スライス方向の深
さに応じて位相を変化させるプロジェクションディフェ
ーズ勾配ＰＤを印加して、静止部からの信号を抑制して
いる。
示す。ｙは位相エンコード方向、ｘは周波数エンコード
方向、ｚはスライス方向である。二次元ＰＣ法では、ス
ライス方向には空間分解能がなく（それ故、スライス方
向に投影され、二次元となる）、ボクセルＶはスライス
方向に長い直方体になる。しかし、スライス方向に長く
なると、各ボクセルＶに占める血流Ｋの割合が小さくな
り、大部分は静止部になる。そこで、スライス方向の深
さに応じて位相を変化させるプロジェクションディフェ
ーズ勾配ＰＤを印加して、静止部からの信号を抑制して
いる。
【０００４】
【発明が解決しようとする課題】上記従来の二次元ＰＣ
法を用いたＭＲアンギオグラフィー方法では、プロジェ
クションディフェーズ勾配ＰＤを印加しているため、投
影方向に位相が乱されてしまい、流れの方向を判定する
ことが出来なくなる問題点がある。しかし、プロジェク
ションディフェーズ勾配ＰＤを印加しなければ、前述の
ようにボクセルＶがスライス方向に長いときに静止部か
らの信号が支配的になって、いわゆる「桁落ち」を引き
起こし、消え残り誤差を生じてしまう問題点がある。そ
こで、この発明の目的は、ボクセルＶがスライス方向に
長い場合でも流れの方向が判り且つ消え残り誤差を小さ
くすることが出来る二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法及び装置を提供することにある。
法を用いたＭＲアンギオグラフィー方法では、プロジェ
クションディフェーズ勾配ＰＤを印加しているため、投
影方向に位相が乱されてしまい、流れの方向を判定する
ことが出来なくなる問題点がある。しかし、プロジェク
ションディフェーズ勾配ＰＤを印加しなければ、前述の
ようにボクセルＶがスライス方向に長いときに静止部か
らの信号が支配的になって、いわゆる「桁落ち」を引き
起こし、消え残り誤差を生じてしまう問題点がある。そ
こで、この発明の目的は、ボクセルＶがスライス方向に
長い場合でも流れの方向が判り且つ消え残り誤差を小さ
くすることが出来る二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法及び装置を提供することにある。
【０００５】
【課題を解決するための手段】この発明の二次元ＰＣ法
を用いたＭＲアンギオグラフィー方法は、ＲＦパルスお
よびスライス勾配を印加して所定のスライスを選択励起
し、異なる２種類のバイポーラ勾配パルスを印加して２
種類のデータを収集し、それら２種類のデータの持つ位
相成分の違いから流動部のイメージを取得する二次元Ｐ
Ｃ法を用いたＭＲアンギオグラフィー方法において、Ｒ
Ｆパルスを印加する直前に、マグネタイゼーショントラ
ンスファ法により静止部の磁化を飽和させ、且つ、プロ
ジェクションディフェーズ勾配を印加しないことを構成
上の特徴とするものである。
を用いたＭＲアンギオグラフィー方法は、ＲＦパルスお
よびスライス勾配を印加して所定のスライスを選択励起
し、異なる２種類のバイポーラ勾配パルスを印加して２
種類のデータを収集し、それら２種類のデータの持つ位
相成分の違いから流動部のイメージを取得する二次元Ｐ
Ｃ法を用いたＭＲアンギオグラフィー方法において、Ｒ
Ｆパルスを印加する直前に、マグネタイゼーショントラ
ンスファ法により静止部の磁化を飽和させ、且つ、プロ
ジェクションディフェーズ勾配を印加しないことを構成
上の特徴とするものである。
【０００６】また、この発明のＭＲアンギオグラフィー
装置は、ＲＦパルスおよびスライス勾配を印加して所定
のスライスを選択励起し、異なる２種類のバイポーラ勾
配パルスを印加して２種類のデータを収集し、それら２
種類のデータの持つ位相成分の違いから流動部のイメー
ジを取得する二次元ＰＣ法を用いたＭＲアンギオグラフ
ィー方法を実施するＭＲアンギオグラフィー装置におい
て、ＲＦパルスを印加する直前に、マグネタイゼーショ
ントランスファ法により静止部の磁化強度を減少させ、
且つ、プロジェクションディフェーズ勾配を印加しない
パルスシーケンスを実施する制御手段を具備したことを
構成上の特徴とするものである。
装置は、ＲＦパルスおよびスライス勾配を印加して所定
のスライスを選択励起し、異なる２種類のバイポーラ勾
配パルスを印加して２種類のデータを収集し、それら２
種類のデータの持つ位相成分の違いから流動部のイメー
ジを取得する二次元ＰＣ法を用いたＭＲアンギオグラフ
ィー方法を実施するＭＲアンギオグラフィー装置におい
て、ＲＦパルスを印加する直前に、マグネタイゼーショ
ントランスファ法により静止部の磁化強度を減少させ、
且つ、プロジェクションディフェーズ勾配を印加しない
パルスシーケンスを実施する制御手段を具備したことを
構成上の特徴とするものである。
【０００７】
【作用】この発明の二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法及び装置は、次の特徴を持っている。 マグネタイゼーショントランスファ法（Ｍagnetizati
on Ｔransfer法、以下、ＭＴ法という）により静止部の
磁化強度を減少させてから、ＲＦパルス以下を印加し、
データを収集する。
グラフィー方法及び装置は、次の特徴を持っている。 マグネタイゼーショントランスファ法（Ｍagnetizati
on Ｔransfer法、以下、ＭＴ法という）により静止部の
磁化強度を減少させてから、ＲＦパルス以下を印加し、
データを収集する。
【０００８】プロジェクションディフェーズ勾配を印
加しない。
加しない。
【０００９】上記により、プロジェクションディフェ
ーズ勾配を印加しなくても、静止部からの信号を抑制可
能となる。上記により、投影方向に位相が乱されるこ
とがなくなるので、流れの方向の情報を失わない。以上
のため、ボクセルＶがスライス方向に長い場合でも、流
れの方向が判り、且つ、消え残り誤差を小さくすること
が出来るようになる。
ーズ勾配を印加しなくても、静止部からの信号を抑制可
能となる。上記により、投影方向に位相が乱されるこ
とがなくなるので、流れの方向の情報を失わない。以上
のため、ボクセルＶがスライス方向に長い場合でも、流
れの方向が判り、且つ、消え残り誤差を小さくすること
が出来るようになる。
【００１０】
【実施例】以下、図に示す実施例に基づいてこの発明を
さらに詳細に説明する。なお、これによりこの発明が限
定されるものではない。図１は、この発明の二次元ＰＣ
法を用いたＭＲアンギオグラフィー方法を実施するＭＲ
Ｉ装置（すなわち、ＭＲアンギオグラフィー装置として
機能する）のブロック図である。計算機２は、操作卓１
３からの指示に基づき、全体の作動を制御する。シーケ
ンスコントローラ３は、記憶しているシーケンスに基づ
いて、勾配磁場駆動回路４を作動させ、マグネットアセ
ンブリ５の勾配磁場コイルで勾配磁場を発生させる。ま
た、ゲート変調回路７を制御し、ＲＦ発振回路６で発生
したＲＦパルスを所定の波形に変調して、ＲＦ電力増幅
器８からマグネットアセンブリ５の送信コイルに加え
る。
さらに詳細に説明する。なお、これによりこの発明が限
定されるものではない。図１は、この発明の二次元ＰＣ
法を用いたＭＲアンギオグラフィー方法を実施するＭＲ
Ｉ装置（すなわち、ＭＲアンギオグラフィー装置として
機能する）のブロック図である。計算機２は、操作卓１
３からの指示に基づき、全体の作動を制御する。シーケ
ンスコントローラ３は、記憶しているシーケンスに基づ
いて、勾配磁場駆動回路４を作動させ、マグネットアセ
ンブリ５の勾配磁場コイルで勾配磁場を発生させる。ま
た、ゲート変調回路７を制御し、ＲＦ発振回路６で発生
したＲＦパルスを所定の波形に変調して、ＲＦ電力増幅
器８からマグネットアセンブリ５の送信コイルに加え
る。
【００１１】マグネットアセンブリ５の受信コイルで得
られたＮＭＲ信号は、前置増幅器９を介して位相検波器
１０に入力され、さらにＡＤ変換器１１を介して計算機
２に入力される。計算機２は、ＡＤ変換器１１から得た
ＮＭＲ信号のデータに基づき、イメージを再構成し、表
示装置１２で表示する。この発明の二次元ＰＣ法を用い
たＭＲアンギオグラフィー方法は、計算機２およびシー
ケンスコントローラ３に記憶された手順により実施され
る。
られたＮＭＲ信号は、前置増幅器９を介して位相検波器
１０に入力され、さらにＡＤ変換器１１を介して計算機
２に入力される。計算機２は、ＡＤ変換器１１から得た
ＮＭＲ信号のデータに基づき、イメージを再構成し、表
示装置１２で表示する。この発明の二次元ＰＣ法を用い
たＭＲアンギオグラフィー方法は、計算機２およびシー
ケンスコントローラ３に記憶された手順により実施され
る。
【００１２】図２は、この発明の二次元ＰＣ法を用いた
ＭＲアンギオグラフィー方法の一実施例のフロー図であ
る。ユーザが、マグネットアセンブリ５に被検体をセッ
ティングしたのち、操作卓１３を用いて二次元ＰＣ法を
用いたＭＲアンギオグラフィー方法の実施の指示を与え
ると、計算機２が、図２の処理を実行する。
ＭＲアンギオグラフィー方法の一実施例のフロー図であ
る。ユーザが、マグネットアセンブリ５に被検体をセッ
ティングしたのち、操作卓１３を用いて二次元ＰＣ法を
用いたＭＲアンギオグラフィー方法の実施の指示を与え
ると、計算機２が、図２の処理を実行する。
【００１３】ステップＳ１では、図３に示すようなパル
スシーケンスに基づいてデータ収集を行う。このパルス
シーケンスＡでは、まず、プロトンのラーモア周波数に
合わせたＲＦ周波数のポリノミアルパルス（Ｐolynomia
l Ｐulse）ＰＰを印加する。これにより、マグネタイゼ
ーショントランスファ効果を生じ、血流には影響を及ぼ
さずに、静止部（軟部組織）を選択的に励起する。な
お、図２のポリノミアルパルスＰＰは、いわゆる“１，
１パルス”であるが、“１，２，１パルス”や“１，
３，３，１パルス”などを用いてもよい。次に、スポイ
ラ勾配ＳＰを印加して、前記ポリノミアルパルスＰＰに
より励起された静止部の磁化を飽和状態にする。
スシーケンスに基づいてデータ収集を行う。このパルス
シーケンスＡでは、まず、プロトンのラーモア周波数に
合わせたＲＦ周波数のポリノミアルパルス（Ｐolynomia
l Ｐulse）ＰＰを印加する。これにより、マグネタイゼ
ーショントランスファ効果を生じ、血流には影響を及ぼ
さずに、静止部（軟部組織）を選択的に励起する。な
お、図２のポリノミアルパルスＰＰは、いわゆる“１，
１パルス”であるが、“１，２，１パルス”や“１，
３，３，１パルス”などを用いてもよい。次に、スポイ
ラ勾配ＳＰを印加して、前記ポリノミアルパルスＰＰに
より励起された静止部の磁化を飽和状態にする。
【００１４】次に、プロジェクションディフェーズ勾配
ＰＤを印加しないという制限以外は従来のパルスシーケ
ンスＢ（図７）と同様にして、データを収集する。すな
わち、ＲＦパルスαおよびスライス勾配Ｓを印加してス
ライスを選択励起し、実線で示すバイポーラ勾配パルス
Ｇ１と、リード勾配Ｒと、位相エンコード勾配ＰＥを印
加して、第１のＲａｗデータＲ１を収集する。次に、実
線で示すバイポーラ勾配パルスＧ１を点線で示すバイポ
ーラ勾配パルスＧ２に替えて、第２のＲａｗデータＲ２
を収集する。
ＰＤを印加しないという制限以外は従来のパルスシーケ
ンスＢ（図７）と同様にして、データを収集する。すな
わち、ＲＦパルスαおよびスライス勾配Ｓを印加してス
ライスを選択励起し、実線で示すバイポーラ勾配パルス
Ｇ１と、リード勾配Ｒと、位相エンコード勾配ＰＥを印
加して、第１のＲａｗデータＲ１を収集する。次に、実
線で示すバイポーラ勾配パルスＧ１を点線で示すバイポ
ーラ勾配パルスＧ２に替えて、第２のＲａｗデータＲ２
を収集する。
【００１５】ステップＳ２では、ＲａｗデータＲ１から
複素画像ｆ１（ｉ，ｊ）を得、ＲａｗデータＲ２から複
素画像ｆ２（ｉ，ｊ）を得る。ここで、ｉは位相エンコ
ード方向のピクセル番号であり、ｊは周波数エンコード
方向のピクセル番号である。ステップＳ３では、複素画
像ｆ１（ｉ，ｊ），ｆ２（ｉ，ｊ）について、静止部の
位相をゼロに一致させる０次の位相補正と、渦電流によ
る位相ズレを補正する１次の位相補正を施す。０次の位
相補正と１次の位相補正を施したのちの複素画像をＦ１
（ｉ，ｊ），Ｆ２（ｉ，ｊ）で示す。
複素画像ｆ１（ｉ，ｊ）を得、ＲａｗデータＲ２から複
素画像ｆ２（ｉ，ｊ）を得る。ここで、ｉは位相エンコ
ード方向のピクセル番号であり、ｊは周波数エンコード
方向のピクセル番号である。ステップＳ３では、複素画
像ｆ１（ｉ，ｊ），ｆ２（ｉ，ｊ）について、静止部の
位相をゼロに一致させる０次の位相補正と、渦電流によ
る位相ズレを補正する１次の位相補正を施す。０次の位
相補正と１次の位相補正を施したのちの複素画像をＦ１
（ｉ，ｊ），Ｆ２（ｉ，ｊ）で示す。
【００１６】ステップＳ４にて、複素画像Ｆ１（ｉ，
ｊ）と複素画像Ｆ２（ｉ，ｊ）との差Ｄ（ｉ，ｊ）をと
る。この差Ｄ（ｉ，ｊ）の絶対値｜Ｄ（ｉ，ｊ）｜を画
像化したものが、従来の二次元ＰＣ法を用いたＭＲアン
ギオグラフィー方法による血流イメージに相当する。
ｊ）と複素画像Ｆ２（ｉ，ｊ）との差Ｄ（ｉ，ｊ）をと
る。この差Ｄ（ｉ，ｊ）の絶対値｜Ｄ（ｉ，ｊ）｜を画
像化したものが、従来の二次元ＰＣ法を用いたＭＲアン
ギオグラフィー方法による血流イメージに相当する。
【００１７】ステップＳ５では、差Ｄ（ｉ，ｊ）の偏角
arg［Ｄ（ｉ，ｊ）］（＝arctan［虚数部／実数部］）
から血流の方向を判定する。例えば、 −９０゜≦arg［Ｄ（ｉ，ｊ）］＜９０゜のとき Ｓ
（ｉ，ｊ）＝＋１ arg［Ｄ（ｉ，ｊ）］＜−９０゜ のとき Ｓ
（ｉ，ｊ）＝−１ ９０゜≦arg［Ｄ（ｉ，ｊ）］ のとき Ｓ
（ｉ，ｊ）＝−１ とする。Ｓ（ｉ，ｊ）は、血流の方向を表わす。
arg［Ｄ（ｉ，ｊ）］（＝arctan［虚数部／実数部］）
から血流の方向を判定する。例えば、 −９０゜≦arg［Ｄ（ｉ，ｊ）］＜９０゜のとき Ｓ
（ｉ，ｊ）＝＋１ arg［Ｄ（ｉ，ｊ）］＜−９０゜ のとき Ｓ
（ｉ，ｊ）＝−１ ９０゜≦arg［Ｄ（ｉ，ｊ）］ のとき Ｓ
（ｉ，ｊ）＝−１ とする。Ｓ（ｉ，ｊ）は、血流の方向を表わす。
【００１８】ステップＳ６では、表示画面上のピクセル
の輝度Ｉ（ｉ，ｊ）を、 Ｉ（ｉ，ｊ）＝Ｉｓ＋Ｓ（ｉ，ｊ）・｜Ｄ（ｉ，ｊ）｜ のように算出して、血流イメージを表示する。その血流
イメージでは、差Ｄ（ｉ，ｊ）の絶対値｜Ｄ（ｉ，ｊ）
｜がゼロであるピクセル、即ち、静止部に対応するピク
セルの輝度がＩｓとなり、一方向に流れる血流がＩｓよ
り明るい輝度となり、反対方向に流れる血流がＩｓより
暗い輝度となる。また、その明るさ，暗さの度合が流速
を表わすようになる。なお、超音波診断装置のＣＦＭの
ように、血流の向きに応じて色を変えてカラー表示とす
ることも可能である。
の輝度Ｉ（ｉ，ｊ）を、 Ｉ（ｉ，ｊ）＝Ｉｓ＋Ｓ（ｉ，ｊ）・｜Ｄ（ｉ，ｊ）｜ のように算出して、血流イメージを表示する。その血流
イメージでは、差Ｄ（ｉ，ｊ）の絶対値｜Ｄ（ｉ，ｊ）
｜がゼロであるピクセル、即ち、静止部に対応するピク
セルの輝度がＩｓとなり、一方向に流れる血流がＩｓよ
り明るい輝度となり、反対方向に流れる血流がＩｓより
暗い輝度となる。また、その明るさ，暗さの度合が流速
を表わすようになる。なお、超音波診断装置のＣＦＭの
ように、血流の向きに応じて色を変えてカラー表示とす
ることも可能である。
【００１９】以上の二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法では、ＭＴ法により、血流には影響を及
ぼすことなく，静止部からの信号を抑制している。この
ため、血流の流速だけではなく、方向も正確に算出でき
るようになり、血流の流速および方向を好適に画像化す
ることが出来る。
グラフィー方法では、ＭＴ法により、血流には影響を及
ぼすことなく，静止部からの信号を抑制している。この
ため、血流の流速だけではなく、方向も正確に算出でき
るようになり、血流の流速および方向を好適に画像化す
ることが出来る。
【００２０】図４は、この発明の二次元ＰＣ法を用いた
ＭＲアンギオグラフィー方法の他の実施例のフロー図で
ある。ステップＳ１〜Ｓ３は、図２と同様である。ステ
ップＳ１４では、複素画像Ｆ１（ｉ，ｊ）の偏角と複素
画像Ｆ２（ｉ，ｊ）の偏角との差Ｐ（ｉ，ｊ）をとる。
ステップＳ１５では、差Ｐ（ｉ，ｊ）から血流の方向を
判定する。例えば、 −９０゜≦Ｐ（ｉ，ｊ）＜９０゜のとき Ｓ（ｉ，ｊ）
＝＋１ Ｐ（ｉ，ｊ）＜−９０゜ のとき Ｓ（ｉ，ｊ）
＝−１ ９０゜≦Ｐ（ｉ，ｊ） のとき Ｓ（ｉ，ｊ）
＝−１ とする。ステップＳ１６では、表示画面上の信号強度
Ｉ’（ｉ，ｊ）を、 Ｉ'（ｉ，ｊ）＝Ｉｓ＋Ｓ（ｉ，ｊ）・｜Ｆ１（ｉ，ｊ）
｜・｜Ｐ（ｉ，ｊ）｜ のように算出して、絶対値｜Ｆ１（ｉ，ｊ）｜でマスク
した血流イメージを表示する。なお、上式の｜Ｆ１
（ｉ，ｊ）｜に代えて、｜Ｆ２（ｉ，ｊ）｜を用いても
よい。
ＭＲアンギオグラフィー方法の他の実施例のフロー図で
ある。ステップＳ１〜Ｓ３は、図２と同様である。ステ
ップＳ１４では、複素画像Ｆ１（ｉ，ｊ）の偏角と複素
画像Ｆ２（ｉ，ｊ）の偏角との差Ｐ（ｉ，ｊ）をとる。
ステップＳ１５では、差Ｐ（ｉ，ｊ）から血流の方向を
判定する。例えば、 −９０゜≦Ｐ（ｉ，ｊ）＜９０゜のとき Ｓ（ｉ，ｊ）
＝＋１ Ｐ（ｉ，ｊ）＜−９０゜ のとき Ｓ（ｉ，ｊ）
＝−１ ９０゜≦Ｐ（ｉ，ｊ） のとき Ｓ（ｉ，ｊ）
＝−１ とする。ステップＳ１６では、表示画面上の信号強度
Ｉ’（ｉ，ｊ）を、 Ｉ'（ｉ，ｊ）＝Ｉｓ＋Ｓ（ｉ，ｊ）・｜Ｆ１（ｉ，ｊ）
｜・｜Ｐ（ｉ，ｊ）｜ のように算出して、絶対値｜Ｆ１（ｉ，ｊ）｜でマスク
した血流イメージを表示する。なお、上式の｜Ｆ１
（ｉ，ｊ）｜に代えて、｜Ｆ２（ｉ，ｊ）｜を用いても
よい。
【００２１】上記実施例では、ポリノミアルパルスＰＰ
を用いたＭＴ法を採用したが、キャリア周波数をラーモ
ア周波数から数ｋHzずらしたＲＦ周波数を一定期間照射
したり、これをｓｉｎｃ関数で振幅変調したパルスを用
いるなどの他のＭＴ法を採用してもよい。
を用いたＭＴ法を採用したが、キャリア周波数をラーモ
ア周波数から数ｋHzずらしたＲＦ周波数を一定期間照射
したり、これをｓｉｎｃ関数で振幅変調したパルスを用
いるなどの他のＭＴ法を採用してもよい。
【００２２】さて、ＭＴ法は、筋肉や脳実質の信号強度
を落とすには有効であるが、脂肪に対しては有効でな
い。このため、静止部に脂肪が多く含まれている腹部な
どを診断部位とする場合には、ＭＴ法に加えて、水（血
流）からの信号と脂肪からの信号との位相が反転するよ
うにエコー時間を設定したり、ＣＨＥＳＳ法などの脂肪
抑制パルスを用いるのが望ましい。図５に、ＭＴ法のポ
リノミアルパルスＰＰおよびスポイラ勾配ＳＰ（実線）
と、ＣＨＥＳＳ法などの脂肪抑制パルスＣおよびスポイ
ラ勾配Ｄと（点線）を、ビュースキャンごとに交互に印
加するパルスシーケンスＡ’を示す。これによって、脂
肪が多く含まれる場合でも、静止部からの信号を好適に
減衰させることが出来る。
を落とすには有効であるが、脂肪に対しては有効でな
い。このため、静止部に脂肪が多く含まれている腹部な
どを診断部位とする場合には、ＭＴ法に加えて、水（血
流）からの信号と脂肪からの信号との位相が反転するよ
うにエコー時間を設定したり、ＣＨＥＳＳ法などの脂肪
抑制パルスを用いるのが望ましい。図５に、ＭＴ法のポ
リノミアルパルスＰＰおよびスポイラ勾配ＳＰ（実線）
と、ＣＨＥＳＳ法などの脂肪抑制パルスＣおよびスポイ
ラ勾配Ｄと（点線）を、ビュースキャンごとに交互に印
加するパルスシーケンスＡ’を示す。これによって、脂
肪が多く含まれる場合でも、静止部からの信号を好適に
減衰させることが出来る。
【００２３】なお、上記実施例のステップＳ３におい
て、Ｒａｗデータから複素画像を得る際にｋ空間上の全
てのＲａｗデータを複素加算した和の位相角で全てのＲ
ａｗデータについて位相回転させて０次の位相補正を施
してもよい。また、うず電流を補正できるＳＧＣ（Ｓhi
elded Ｇradient Ｃoil)を備えたＭＲＩ装置では、１次
の位相補正を施さなくてもよい。
て、Ｒａｗデータから複素画像を得る際にｋ空間上の全
てのＲａｗデータを複素加算した和の位相角で全てのＲ
ａｗデータについて位相回転させて０次の位相補正を施
してもよい。また、うず電流を補正できるＳＧＣ（Ｓhi
elded Ｇradient Ｃoil)を備えたＭＲＩ装置では、１次
の位相補正を施さなくてもよい。
【００２４】また、上記実施例のステップＳ５では、差
Ｄ（ｉ，ｊ）の偏角が±９０゜以内か否かにより血流の
方向を判定しているが、例えば、 −４５゜≦arg［Ｄ（ｉ，ｊ）］＜４５゜ のとき
Ｓ（ｉ，ｊ）＝＋１ ４５゜≦arg［Ｄ（ｉ，ｊ）］＜１３５゜ のとき
Ｓ（ｉ，ｊ）＝０ １３５゜≦arg［Ｄ（ｉ，ｊ）］ のとき
Ｓ（ｉ，ｊ）＝−１ −１３５゜≦arg［Ｄ（ｉ，ｊ）］＜−４５゜のとき
Ｓ（ｉ，ｊ）＝０ arg［Ｄ（ｉ，ｊ）］＜−１３５゜ のとき
Ｓ（ｉ，ｊ）＝−１ のように判定し、前記偏角の大きいところ（信頼性が低
いところ）では血流を表示させないようにしてもよい。
Ｄ（ｉ，ｊ）の偏角が±９０゜以内か否かにより血流の
方向を判定しているが、例えば、 −４５゜≦arg［Ｄ（ｉ，ｊ）］＜４５゜ のとき
Ｓ（ｉ，ｊ）＝＋１ ４５゜≦arg［Ｄ（ｉ，ｊ）］＜１３５゜ のとき
Ｓ（ｉ，ｊ）＝０ １３５゜≦arg［Ｄ（ｉ，ｊ）］ のとき
Ｓ（ｉ，ｊ）＝−１ −１３５゜≦arg［Ｄ（ｉ，ｊ）］＜−４５゜のとき
Ｓ（ｉ，ｊ）＝０ arg［Ｄ（ｉ，ｊ）］＜−１３５゜ のとき
Ｓ（ｉ，ｊ）＝−１ のように判定し、前記偏角の大きいところ（信頼性が低
いところ）では血流を表示させないようにしてもよい。
【００２５】
【発明の効果】この発明の二次元ＰＣ法を用いたＭＲア
ンギオグラフィー方法及び装置によれば、ＭＴ法により
流動部には影響を及ぼすことなく静止部からの信号を抑
制するため、流れの大きさだけでなく、方向をも正確に
算出し、画像化することが出来る。
ンギオグラフィー方法及び装置によれば、ＭＴ法により
流動部には影響を及ぼすことなく静止部からの信号を抑
制するため、流れの大きさだけでなく、方向をも正確に
算出し、画像化することが出来る。
【図１】この発明の二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法を実施するためのＭＲＩ装置のブロック
図である。
グラフィー方法を実施するためのＭＲＩ装置のブロック
図である。
【図２】この発明の二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法の一実施例にかかるフロー図である。
グラフィー方法の一実施例にかかるフロー図である。
【図３】この発明の二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法の一実施例にかかるパルスシーケンスの
例示図である。
グラフィー方法の一実施例にかかるパルスシーケンスの
例示図である。
【図４】この発明の二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法の他の実施例にかかるフロー図である。
グラフィー方法の他の実施例にかかるフロー図である。
【図５】この発明の二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法のさらに他の実施例にかかるパルスシー
ケンスの例示図である。
グラフィー方法のさらに他の実施例にかかるパルスシー
ケンスの例示図である。
【図６】従来の二次元ＰＣ法を用いたＭＲアンギオグラ
フィー方法にかかるパルスシーケンスの例示図である。
フィー方法にかかるパルスシーケンスの例示図である。
【図７】二次元ＰＣ法により選択励起されるスライスの
例示図である。
例示図である。
１ ＭＲＩ装置 ２ 計算機 Ａ，Ｂ パルスシーケンス Ｇ１ バイポーラ勾配パルス Ｋ 血流 ＰＤ プロジェクションディフェーズ勾配 ＰＰ ポリノミアルパルス ＳＰ スポイラ勾配
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.5 識別記号 庁内整理番号 ＦＩ 技術表示箇所 9219−2Ｊ Ｇ０１Ｎ 24/08 Ｙ
Claims (4)
- 【請求項１】 ＲＦパルスおよびスライス勾配を印加し
て所定のスライスを選択励起し、異なる２種類のバイポ
ーラ勾配パルスを印加して２種類のデータを収集し、そ
れら２種類のデータの持つ位相成分の違いから流動部の
イメージを取得する二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法において、 ＲＦパルスを印加する直前に、マグネタイゼーショント
ランスファ法により静止部の磁化強度を減少させ、且
つ、プロジェクションディフェーズ勾配を印加しないこ
とを特徴とする二次元ＰＣ法を用いたＭＲアンギオグラ
フィー方法。 - 【請求項２】 請求項１に記載の二次元ＰＣ法を用いた
ＭＲアンギオグラフィー方法において、 ＣＨＥＳＳ法の如き脂肪からの信号強度を抑制する方法
を併用することを特徴とする二次元ＰＣ法を用いたＭＲ
アンギオグラフィー方法。 - 【請求項３】 ＲＦパルスおよびスライス勾配を印加し
て所定のスライスを選択励起し、異なる２種類のバイポ
ーラ勾配パルスを印加して２種類のデータを収集し、そ
れら２種類のデータの持つ位相成分の違いから流動部の
イメージを取得する二次元ＰＣ法を用いたＭＲアンギオ
グラフィー方法を実施するＭＲアンギオグラフィー装置
において、 ＲＦパルスを印加する直前に、マグネタイゼーショント
ランスファ法により静止部の磁化強度を減少させ、且
つ、プロジェクションディフェーズ勾配を印加しないパ
ルスシーケンスを実施する制御手段を具備したことを特
徴とするＭＲアンギオグラフィー装置。 - 【請求項４】 請求項３に記載のＭＲアンギオグラフィ
ー装置において、 前記制御手段は、ＣＨＥＳＳ法の如き脂肪からの信号強
度を抑制する方法を前記パルスシーケンスに加えたパル
スシーケンスを実施することを特徴とするＭＲアンギオ
グラフィー装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP08991293A JP3343392B2 (ja) | 1993-04-16 | 1993-04-16 | Ｍｒアンギオグラフィー装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP08991293A JP3343392B2 (ja) | 1993-04-16 | 1993-04-16 | Ｍｒアンギオグラフィー装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06296600A true JPH06296600A (ja) | 1994-10-25 |
| JP3343392B2 JP3343392B2 (ja) | 2002-11-11 |
Family
ID=13983930
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP08991293A Expired - Lifetime JP3343392B2 (ja) | 1993-04-16 | 1993-04-16 | Ｍｒアンギオグラフィー装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3343392B2 (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11299753A (ja) * | 1998-04-17 | 1999-11-02 | Toshiba Corp | Ｍｒｉ装置およびｍｒ撮像方法 |
| KR100335783B1 (ko) * | 1999-09-08 | 2002-05-09 | 이민화 | 듀얼밴드의 사전포화를 이용하여 영상영역의 부분포화를 방지하는 방법 |
| KR100459097B1 (ko) * | 2000-10-31 | 2004-12-03 | 지이 메디컬 시스템즈 글로발 테크놀러지 캄파니 엘엘씨 | Ｍｒ 촬상 방법, 잔류 자화량 측정 방법 및 ｍｒｉ 장치 |
| US8131338B2 (en) | 1998-04-17 | 2012-03-06 | Kabushiki Kaisha Toshiba | MR imaging providing tissue/blood contrast image |
| JP2014514112A (ja) * | 2011-04-21 | 2014-06-19 | コーニンクレッカ フィリップス エヌ ヴェ | 交互的ブラック及びブライトブラッド動的コントラスト促進（ｄｃｅ）ｍｒｉ |
-
1993
- 1993-04-16 JP JP08991293A patent/JP3343392B2/ja not_active Expired - Lifetime
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11299753A (ja) * | 1998-04-17 | 1999-11-02 | Toshiba Corp | Ｍｒｉ装置およびｍｒ撮像方法 |
| US8131338B2 (en) | 1998-04-17 | 2012-03-06 | Kabushiki Kaisha Toshiba | MR imaging providing tissue/blood contrast image |
| US8560049B2 (en) | 1998-04-17 | 2013-10-15 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging system |
| KR100335783B1 (ko) * | 1999-09-08 | 2002-05-09 | 이민화 | 듀얼밴드의 사전포화를 이용하여 영상영역의 부분포화를 방지하는 방법 |
| KR100459097B1 (ko) * | 2000-10-31 | 2004-12-03 | 지이 메디컬 시스템즈 글로발 테크놀러지 캄파니 엘엘씨 | Ｍｒ 촬상 방법, 잔류 자화량 측정 방법 및 ｍｒｉ 장치 |
| JP2014514112A (ja) * | 2011-04-21 | 2014-06-19 | コーニンクレッカ フィリップス エヌ ヴェ | 交互的ブラック及びブライトブラッド動的コントラスト促進（ｄｃｅ）ｍｒｉ |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3343392B2 (ja) | 2002-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1060706B1 (en) | Magnetic resonance imaging device | |
| US8588890B2 (en) | Magnetic resonance imaging apparatus | |
| US5565777A (en) | Method/apparatus for NMR imaging using an imaging scheme sensitive to inhomogeneity and a scheme insensitive to inhomogeneity in a single imaging step | |
| US5115812A (en) | Magnetic resonance imaging method for moving object | |
| US6259250B1 (en) | Method and apparatus for reducing artifacts in echo planar imaging | |
| US6285187B1 (en) | Method and apparatus for reducing artifacts in echo planar imaging | |
| US20010004211A1 (en) | Magnetic resonance imaging apparatus | |
| US7145334B2 (en) | Method and magnetic resonance tomography apparatus for producing phase-coded flow images | |
| US20030160616A1 (en) | MRI apparatus | |
| JP3877926B2 (ja) | 核磁気共鳴撮影装置 | |
| JP2002306445A (ja) | 磁気共鳴イメージング装置及び画像処理方法 | |
| JPH04364829A (ja) | 磁気共鳴画像処理方法及び装置 | |
| JP3431249B2 (ja) | 磁気共鳴映像装置 | |
| JP3472620B2 (ja) | Ｍｒｉ装置 | |
| JP3343392B2 (ja) | Ｍｒアンギオグラフィー装置 | |
| JP3884227B2 (ja) | 磁気共鳴撮影装置 | |
| JP2001511058A (ja) | 磁気共鳴画像用の方法及び装置 | |
| JP4006080B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3588690B2 (ja) | 磁気共鳴装置 | |
| JP3690874B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3907944B2 (ja) | 磁気共鳴イメージング方法及び装置 | |
| JPH07108288B2 (ja) | Ｎｍｒイメ−ジング方法 | |
| JPH0956695A (ja) | 拡散強調イメージング方法、動的イメージング方法およびｍｒｉ装置 | |
| JPH06169A (ja) | Ｍｒによるディフュージョン測定方法および装置 | |
| JP2001340316A (ja) | 磁気共鳴イメージング装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080823 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090823 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090823 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090823 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100823 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100823 Year of fee payment: 8 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100823 Year of fee payment: 8 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100823 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110823 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110823 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120823 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120823 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130823 Year of fee payment: 11 |
|
| EXPY | Cancellation because of completion of term |