JPH07265277A - Ｍｒ血流情報収集方法およびｍｒｉ装置 - Google Patents
Ｍｒ血流情報収集方法およびｍｒｉ装置Info
- Publication number
- JPH07265277A JPH07265277A JP6059266A JP5926694A JPH07265277A JP H07265277 A JPH07265277 A JP H07265277A JP 6059266 A JP6059266 A JP 6059266A JP 5926694 A JP5926694 A JP 5926694A JP H07265277 A JPH07265277 A JP H07265277A
- Authority
- JP
- Japan
- Prior art keywords
- blood flow
- velocity
- profile
- pulse
- distribution
- 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.)
- Pending
Links
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】
【目的】 フリップ角を変化させるプロファイルを実際
の血流速度に合せて最適に調整する。 【構成】 励起スラブ領域を複数のスライスに分割し位
相コントラスト法により各スライスについてスライス厚
方向速度の２次元分布を取得し（ステップＶ２）、全て
のスライスの前記２次元分布を合算して励起スラブ領域
内の血流の速度分布を取得し（ステップＶ３）、その速
度分布を表示し（ステップＶ４）、操作者に対象血流速
度を設定させ（ステップＶ５）、その対象血流速度から
励起スラブ領域の厚さ方向にフリップ角を変化させる変
化関数を算出しシンクパルスをフーリエ変換した基本関
数に前記変化関数を乗算してＲＦパルスのプロファイル
を決定し（ステップＶ６）、そのプロファイルを逆フー
リエ変換したＲＦパルスを用いて３次元ＴＯＦアンギオ
グラフィを行う（ステップＶ７）。 【効果】 被検体や撮像部位により血流速度Ｖｏが変わ
っても、常に最適のプロファイルのＲＦパルスを用いる
ことが出来る。
の血流速度に合せて最適に調整する。 【構成】 励起スラブ領域を複数のスライスに分割し位
相コントラスト法により各スライスについてスライス厚
方向速度の２次元分布を取得し（ステップＶ２）、全て
のスライスの前記２次元分布を合算して励起スラブ領域
内の血流の速度分布を取得し（ステップＶ３）、その速
度分布を表示し（ステップＶ４）、操作者に対象血流速
度を設定させ（ステップＶ５）、その対象血流速度から
励起スラブ領域の厚さ方向にフリップ角を変化させる変
化関数を算出しシンクパルスをフーリエ変換した基本関
数に前記変化関数を乗算してＲＦパルスのプロファイル
を決定し（ステップＶ６）、そのプロファイルを逆フー
リエ変換したＲＦパルスを用いて３次元ＴＯＦアンギオ
グラフィを行う（ステップＶ７）。 【効果】 被検体や撮像部位により血流速度Ｖｏが変わ
っても、常に最適のプロファイルのＲＦパルスを用いる
ことが出来る。
Description
【０００１】
【産業上の利用分野】この発明は、ＭＲ（Magnetic Res
onance）血流情報収集方法（アンギオグラフィ）および
ＭＲＩ（Magnetic Resonance Imaging）装置に関する。
さらに詳しくは、励起スラブ（Slab）領域の厚さ方向に
フリップ角が徐々に変化するようなプロファイルを持つ
ＲＦ（高周波）パルスを照射してその励起スラブ領域内
の血流情報を収集するＭＲ血流情報収集方法およびＭＲ
Ｉ装置に関する。
onance）血流情報収集方法（アンギオグラフィ）および
ＭＲＩ（Magnetic Resonance Imaging）装置に関する。
さらに詳しくは、励起スラブ（Slab）領域の厚さ方向に
フリップ角が徐々に変化するようなプロファイルを持つ
ＲＦ（高周波）パルスを照射してその励起スラブ領域内
の血流情報を収集するＭＲ血流情報収集方法およびＭＲ
Ｉ装置に関する。
【０００２】
【従来の技術】従来、３次元ＴＯＦ（Time Of Flight）
アンギオグラフィにおいて、図９に示すように、被検体
Ｈの励起スラブ領域Ｓの厚さ方向にフリップ角θがラン
プ（Ramp）状に変化する（血流方向に大きくなる）よう
なプロファイルを持つＲＦパルスＲを照射することが行
なわれている。このようなプロファイルを持つＲＦパル
スＲを用いるのは、次の理由による。
アンギオグラフィにおいて、図９に示すように、被検体
Ｈの励起スラブ領域Ｓの厚さ方向にフリップ角θがラン
プ（Ramp）状に変化する（血流方向に大きくなる）よう
なプロファイルを持つＲＦパルスＲを照射することが行
なわれている。このようなプロファイルを持つＲＦパル
スＲを用いるのは、次の理由による。
【０００３】図１０に示すように、血流速度Ｖｏの血流
が断面積ａの血管Ｂを図の左から右へ流れているとす
る。また、励起の繰り返し時間をＴＲとする。また、励
起スラブ領域Ｓの厚さをＬ＝３・Ｖｏ・ＴＲとする。ま
た、励起スラブ領域Ｓをその右端から順に厚さＶｏ・Ｔ
ＲずつのスライスＳ１，Ｓ２，Ｓ３に分割する。さら
に、スライスＳ１，Ｓ２，Ｓ３のそれぞれに対するフリ
ップ角をθ１，θ２，θ３とする。すると、単位体積当
りの磁化ベクトルの大きさをＭｏとし、エコー時間をＴ
Ｅとし、緩和時間をＴ１，Ｔ２とするとき、スライスＳ
１がフリップ角θ１で１回励起されて生じる信号強度Ｉ
１は、横磁化を１ＴＲ毎に消失させる高速シーケンスに
おいて、 Ｉ１＝Mo・sin{θ１}・exp{−TE／T2}・ａ・Vo・TR である。このとき、スライスＳ２は、スライスＳ１の位
置でフリップ角θ１で励起され、且つ、現在の位置でフ
リップ角θ２で励起されているから、その信号強度Ｉ２
は、 Ｉ２＝Mo・[(１−exp{−TR／T1})＋cos{θ１}・exp{−TR
／T1｝]・sin{θ２}・exp{−TE／T2}・ａ・Vo・TR となる。同様に、スライスＳ３の信号強度Ｉ３は、 Ｉ３＝Mo・[(１−exp{−TR／T1})＋[(１−exp{−TR／T
1})cos{θ１}＋cos{θ１}・cos{θ２}・exp{−TR／T1｝]・
sin{θ３}・exp{−TE／T2}・ａ・Vo・TR となる。もし、θ１＝θ２＝θ３なら、Ｉ１＞Ｉ２＞Ｉ
３となり、信号強度が減衰してしまうことになる。そこ
で、所望の信号強度Ｉｏを定めると共にフリップ角θ
１，θ２，θ３以外の変数を定めて、 Ｉ１＝Ｉｏ Ｉ２＝Ｉｏ Ｉ３＝Ｉｏ 且つ、 ０≦θ１，θ２，θ３≦π／２ なる条件式を解いてフリップ角θ１，θ２，θ３を求め
ると、適当な傾きでフリップ角θ１，θ２，θ３をラン
プ状に変化させれば、励起スラブ領域Ｓ内で一定の信号
強度Ｉｏとなるように出来ることが判る。このような理
由により、フリップ角θがランプ状に変化する（血流方
向に大きくなる）ようなプロファイルを持つＲＦパルス
Ｒを照射するのである。
が断面積ａの血管Ｂを図の左から右へ流れているとす
る。また、励起の繰り返し時間をＴＲとする。また、励
起スラブ領域Ｓの厚さをＬ＝３・Ｖｏ・ＴＲとする。ま
た、励起スラブ領域Ｓをその右端から順に厚さＶｏ・Ｔ
ＲずつのスライスＳ１，Ｓ２，Ｓ３に分割する。さら
に、スライスＳ１，Ｓ２，Ｓ３のそれぞれに対するフリ
ップ角をθ１，θ２，θ３とする。すると、単位体積当
りの磁化ベクトルの大きさをＭｏとし、エコー時間をＴ
Ｅとし、緩和時間をＴ１，Ｔ２とするとき、スライスＳ
１がフリップ角θ１で１回励起されて生じる信号強度Ｉ
１は、横磁化を１ＴＲ毎に消失させる高速シーケンスに
おいて、 Ｉ１＝Mo・sin{θ１}・exp{−TE／T2}・ａ・Vo・TR である。このとき、スライスＳ２は、スライスＳ１の位
置でフリップ角θ１で励起され、且つ、現在の位置でフ
リップ角θ２で励起されているから、その信号強度Ｉ２
は、 Ｉ２＝Mo・[(１−exp{−TR／T1})＋cos{θ１}・exp{−TR
／T1｝]・sin{θ２}・exp{−TE／T2}・ａ・Vo・TR となる。同様に、スライスＳ３の信号強度Ｉ３は、 Ｉ３＝Mo・[(１−exp{−TR／T1})＋[(１−exp{−TR／T
1})cos{θ１}＋cos{θ１}・cos{θ２}・exp{−TR／T1｝]・
sin{θ３}・exp{−TE／T2}・ａ・Vo・TR となる。もし、θ１＝θ２＝θ３なら、Ｉ１＞Ｉ２＞Ｉ
３となり、信号強度が減衰してしまうことになる。そこ
で、所望の信号強度Ｉｏを定めると共にフリップ角θ
１，θ２，θ３以外の変数を定めて、 Ｉ１＝Ｉｏ Ｉ２＝Ｉｏ Ｉ３＝Ｉｏ 且つ、 ０≦θ１，θ２，θ３≦π／２ なる条件式を解いてフリップ角θ１，θ２，θ３を求め
ると、適当な傾きでフリップ角θ１，θ２，θ３をラン
プ状に変化させれば、励起スラブ領域Ｓ内で一定の信号
強度Ｉｏとなるように出来ることが判る。このような理
由により、フリップ角θがランプ状に変化する（血流方
向に大きくなる）ようなプロファイルを持つＲＦパルス
Ｒを照射するのである。
【０００４】
【発明が解決しようとする課題】上記のように、３次元
ＴＯＦアンギオグラフィにおいて、フリップ角θがラン
プ状に変化するようなプロファイルを持つＲＦパルスＲ
を照射することが行なわれているが、被検体や撮像部位
が異なると血流速度Ｖｏが異なるため、それに合せてフ
リップ角θを変化させるプロファイルを調整する必要が
ある。しかし、従来は、操作者（例えば、医師）が経験
によってプロファイルを調整していたため、必ずしも最
適のプロファイルの設定がなされるとは限らない問題点
があった。そこで、この発明の目的は、フリップ角θを
変化させるプロファイルを、実際の血流速度に合せて最
適に調整することが出来るようにしたＭＲ血流情報収集
方法およびＭＲＩ装置を提供することにある。
ＴＯＦアンギオグラフィにおいて、フリップ角θがラン
プ状に変化するようなプロファイルを持つＲＦパルスＲ
を照射することが行なわれているが、被検体や撮像部位
が異なると血流速度Ｖｏが異なるため、それに合せてフ
リップ角θを変化させるプロファイルを調整する必要が
ある。しかし、従来は、操作者（例えば、医師）が経験
によってプロファイルを調整していたため、必ずしも最
適のプロファイルの設定がなされるとは限らない問題点
があった。そこで、この発明の目的は、フリップ角θを
変化させるプロファイルを、実際の血流速度に合せて最
適に調整することが出来るようにしたＭＲ血流情報収集
方法およびＭＲＩ装置を提供することにある。
【０００５】
【課題を解決するための手段】第１の観点では、この発
明は、励起スラブ領域の厚さ方向にフリップ角が徐々に
変化するようなプロファイルを持つＲＦパルスを照射し
てその励起スラブ領域内の血流情報を収集するＭＲ血流
情報収集方法において、前記励起スラブ領域内の血流の
厚さ方向の速度分布を求め、その速度分布中に含まれる
速度を対象血流速度として選択し、その対象血流速度に
対して前記ＲＦパルスのプロファイルを最適化すること
を特徴とするＭＲ血流情報収集方法を提供する。
明は、励起スラブ領域の厚さ方向にフリップ角が徐々に
変化するようなプロファイルを持つＲＦパルスを照射し
てその励起スラブ領域内の血流情報を収集するＭＲ血流
情報収集方法において、前記励起スラブ領域内の血流の
厚さ方向の速度分布を求め、その速度分布中に含まれる
速度を対象血流速度として選択し、その対象血流速度に
対して前記ＲＦパルスのプロファイルを最適化すること
を特徴とするＭＲ血流情報収集方法を提供する。
【０００６】第２の観点では、この発明は、励起スラブ
領域の厚さ方向にフリップ角が徐々に変化するようなプ
ロファイルを持つＲＦパルスを照射してその励起スラブ
領域内の血流情報を収集するＭＲＩ装置において、前記
励起スラブ領域内の血流の厚さ方向の速度分布を取得す
る速度分布取得手段と、その速度分布中に含まれる速度
を対象血流速度として選択する対象血流速度選択手段
と、その対象血流速度に対して前記ＲＦパルスのプロフ
ァイルを最適化するプロファイル最適化手段とを具備し
たことを特徴とするＭＲＩ装置を提供する。
領域の厚さ方向にフリップ角が徐々に変化するようなプ
ロファイルを持つＲＦパルスを照射してその励起スラブ
領域内の血流情報を収集するＭＲＩ装置において、前記
励起スラブ領域内の血流の厚さ方向の速度分布を取得す
る速度分布取得手段と、その速度分布中に含まれる速度
を対象血流速度として選択する対象血流速度選択手段
と、その対象血流速度に対して前記ＲＦパルスのプロフ
ァイルを最適化するプロファイル最適化手段とを具備し
たことを特徴とするＭＲＩ装置を提供する。
【０００７】上記構成のＭＲＩ装置において、前記速度
分布取得手段は、例えば、前記励起スラブ領域を複数の
スライスに分割し、位相コントラスト法により各スライ
スについてスライス厚方向速度のスライス面内での２次
元分布を取得し、全てのスライスの前記２次元分布を合
算して、前記速度分布を取得するものである。また、上
記構成のＭＲＩ装置において、前記対象血流速度選択手
段は、例えば、前記速度分布を表示し、操作者に対象血
流速度を設定させるものである。また、上記構成のＭＲ
Ｉ装置において、前記プロファイル最適化手段は、例え
ば、励起スラブ領域の厚さ方向にフリップ角を変化させ
る変化関数を前記対象血流速度から算出し、シンクパル
スをフーリエ変換した基本関数に前記変化関数を乗算し
てＲＦパルスのプロファイルを決定するものである。
分布取得手段は、例えば、前記励起スラブ領域を複数の
スライスに分割し、位相コントラスト法により各スライ
スについてスライス厚方向速度のスライス面内での２次
元分布を取得し、全てのスライスの前記２次元分布を合
算して、前記速度分布を取得するものである。また、上
記構成のＭＲＩ装置において、前記対象血流速度選択手
段は、例えば、前記速度分布を表示し、操作者に対象血
流速度を設定させるものである。また、上記構成のＭＲ
Ｉ装置において、前記プロファイル最適化手段は、例え
ば、励起スラブ領域の厚さ方向にフリップ角を変化させ
る変化関数を前記対象血流速度から算出し、シンクパル
スをフーリエ変換した基本関数に前記変化関数を乗算し
てＲＦパルスのプロファイルを決定するものである。
【０００８】
【作用】この発明のＭＲ血流情報収集方法およびＭＲＩ
装置では、励起スラブ領域内の血流の速度分布を実測
し、その速度分布中から対象血流速度を選択し、その対
象血流速度に対してＲＦパルスのプロファイルを最適化
する。従って、被検体や撮像部位により血流速度Ｖｏが
変わっても、常に最適のプロファイルのＲＦパルスを用
いることが出来るようになる。
装置では、励起スラブ領域内の血流の速度分布を実測
し、その速度分布中から対象血流速度を選択し、その対
象血流速度に対してＲＦパルスのプロファイルを最適化
する。従って、被検体や撮像部位により血流速度Ｖｏが
変わっても、常に最適のプロファイルのＲＦパルスを用
いることが出来るようになる。
【０００９】なお、励起スラブ領域を複数のスライスに
分割し、位相コントラスト法により各スライスについて
スライス厚方向速度のスライス面内での２次元分布を取
得し、全てのスライスの前記２次元分布を合算して、前
記速度分布を取得し、その速度分布を表示し、操作者に
対象血流速度を設定させ、その対象血流速度を基にして
励起スラブ領域の厚さ方向にフリップ角を変化させる変
化関数を算出し、シンクパルスをフーリエ変換した基本
関数に前記変化関数を乗算してＲＦパルスのプロファイ
ルを決定すれば、最適のプロファイルのＲＦパルスを好
適に得ることが出来る。
分割し、位相コントラスト法により各スライスについて
スライス厚方向速度のスライス面内での２次元分布を取
得し、全てのスライスの前記２次元分布を合算して、前
記速度分布を取得し、その速度分布を表示し、操作者に
対象血流速度を設定させ、その対象血流速度を基にして
励起スラブ領域の厚さ方向にフリップ角を変化させる変
化関数を算出し、シンクパルスをフーリエ変換した基本
関数に前記変化関数を乗算してＲＦパルスのプロファイ
ルを決定すれば、最適のプロファイルのＲＦパルスを好
適に得ることが出来る。
【００１０】
【実施例】以下、図に示す実施例によりこの発明をさら
に詳しく説明する。なお、これによりこの発明が限定さ
れるものではない。図１は、この発明のＭＲＩ装置の一
実施例のブロック図である。このＭＲＩ装置１００にお
いて、マグネットアセンブリ１は、内部に被検体を挿入
するための空間部分（孔）を有し、この空間部分を取り
まくようにして、被検体に一定の静磁場を印加する静磁
場コイルと、勾配磁場を発生するための勾配磁場コイル
（勾配コイルは、読み出し軸，位相軸，スライス軸のコ
イルを備えている）と、被検体内の原子核のスピンを励
起するためのＲＦパルスを与える送信コイルと、被検体
からのＮＭＲ信号を検出する受信コイル等が配置されて
いる。静磁場コイル，勾配磁場コイル，送信コイルおよ
び受信コイルは、それぞれ主磁場電源２，勾配磁場駆動
回路３，ＲＦ電力増幅器４および前置増幅器５に接続さ
れている。
に詳しく説明する。なお、これによりこの発明が限定さ
れるものではない。図１は、この発明のＭＲＩ装置の一
実施例のブロック図である。このＭＲＩ装置１００にお
いて、マグネットアセンブリ１は、内部に被検体を挿入
するための空間部分（孔）を有し、この空間部分を取り
まくようにして、被検体に一定の静磁場を印加する静磁
場コイルと、勾配磁場を発生するための勾配磁場コイル
（勾配コイルは、読み出し軸，位相軸，スライス軸のコ
イルを備えている）と、被検体内の原子核のスピンを励
起するためのＲＦパルスを与える送信コイルと、被検体
からのＮＭＲ信号を検出する受信コイル等が配置されて
いる。静磁場コイル，勾配磁場コイル，送信コイルおよ
び受信コイルは、それぞれ主磁場電源２，勾配磁場駆動
回路３，ＲＦ電力増幅器４および前置増幅器５に接続さ
れている。
【００１１】シーケンス記憶回路８は、計算機７からの
指令に従い、記憶されているデータ収集パルスシーケン
スに基づいて勾配磁場駆動回路３を操作し、前記マグネ
ットアセンブリ１の勾配磁場コイルから勾配磁場を発生
させると共に、ゲート変調回路９を操作し、ＲＦ発振回
路１０の高周波出力信号を所定タイミング・所定包絡線
のパルス状信号に変調し、それをＲＦパルスとしてＲＦ
電力増幅器４に加え、ＲＦ電力増幅器４でパワー増幅し
た後、前記マグネットアセンブリ１の送信コイルに印加
し、目的の領域を選択励起する。
指令に従い、記憶されているデータ収集パルスシーケン
スに基づいて勾配磁場駆動回路３を操作し、前記マグネ
ットアセンブリ１の勾配磁場コイルから勾配磁場を発生
させると共に、ゲート変調回路９を操作し、ＲＦ発振回
路１０の高周波出力信号を所定タイミング・所定包絡線
のパルス状信号に変調し、それをＲＦパルスとしてＲＦ
電力増幅器４に加え、ＲＦ電力増幅器４でパワー増幅し
た後、前記マグネットアセンブリ１の送信コイルに印加
し、目的の領域を選択励起する。
【００１２】前置増幅器５は、マグネットアセンブリ１
の受信コイルで検出された被検体からのＮＭＲ信号を増
幅し、位相検波器１２に入力する。位相検波器１２は、
ＲＦ発振回路１０の出力を参照信号とし、前置増幅器５
からのＮＭＲ信号を位相検波して、Ａ／Ｄ変換器１１に
与える。Ａ／Ｄ変換器１１は、位相検波後のアナログ信
号をディジタル信号に変換して、計算機７に入力する。
計算機７は、Ａ／Ｄ変換器１１からのデジタル信号に対
する画像再構成演算を行い、目的の領域のイメージ（プ
ロトン密度像）を生成する。このイメージは、表示装置
６にて表示される。また、計算機７は、操作卓１３から
入力された情報を受け取るなどの全体的な制御を受け持
つ。
の受信コイルで検出された被検体からのＮＭＲ信号を増
幅し、位相検波器１２に入力する。位相検波器１２は、
ＲＦ発振回路１０の出力を参照信号とし、前置増幅器５
からのＮＭＲ信号を位相検波して、Ａ／Ｄ変換器１１に
与える。Ａ／Ｄ変換器１１は、位相検波後のアナログ信
号をディジタル信号に変換して、計算機７に入力する。
計算機７は、Ａ／Ｄ変換器１１からのデジタル信号に対
する画像再構成演算を行い、目的の領域のイメージ（プ
ロトン密度像）を生成する。このイメージは、表示装置
６にて表示される。また、計算機７は、操作卓１３から
入力された情報を受け取るなどの全体的な制御を受け持
つ。
【００１３】図２は、ＭＲＩ装置１００における血流情
報収集処理を示す要部フローチャートである。ステップ
Ｖ１では、操作者が、例えば図３に示す励起スラブ領域
Ｓの位置や，各種のパラメータを設定する。ステップＶ
２では、図３に示すように、励起スラブ領域Ｓを複数の
スライスＳｉに分割する。そして、各スライスＳｉにつ
いて、位相コントラスト（ＰＣ）法により、図４に示す
ように、スライス厚方向速度のスライス面内での２次元
分布Ｖｉ（ｘ，ｙ）を取得する。ステップＶ３では、全
てのスライスＳｉの２次元分布を合算して、励起スラブ
領域Ｓ内の血流速度分布を求める。すなわち、ある血流
速度の画素数を全てのスライスＳｉについてカウントし
て、図５に示すヒストグラムを求める。ステップＶ４で
は、血流速度のヒストグラムをディスプレイに表示す
る。ステップＶ５では、操作者に対象血流速度を設定さ
せる。操作者は、図６に示すように、所望の血流速度を
対象血流速度として設定する。
報収集処理を示す要部フローチャートである。ステップ
Ｖ１では、操作者が、例えば図３に示す励起スラブ領域
Ｓの位置や，各種のパラメータを設定する。ステップＶ
２では、図３に示すように、励起スラブ領域Ｓを複数の
スライスＳｉに分割する。そして、各スライスＳｉにつ
いて、位相コントラスト（ＰＣ）法により、図４に示す
ように、スライス厚方向速度のスライス面内での２次元
分布Ｖｉ（ｘ，ｙ）を取得する。ステップＶ３では、全
てのスライスＳｉの２次元分布を合算して、励起スラブ
領域Ｓ内の血流速度分布を求める。すなわち、ある血流
速度の画素数を全てのスライスＳｉについてカウントし
て、図５に示すヒストグラムを求める。ステップＶ４で
は、血流速度のヒストグラムをディスプレイに表示す
る。ステップＶ５では、操作者に対象血流速度を設定さ
せる。操作者は、図６に示すように、所望の血流速度を
対象血流速度として設定する。
【００１４】ステップＶ６では、設定された対象血流速
度を前記血流速度Ｖｏとして前記条件式を解いて、図７
に示すように、フリップ角θｉを励起スラブ領域Ｓの厚
さ方向に変化させる傾き関数Ｒ（ω）を算出する。な
お、励起スラブ領域Ｓの厚さ方向の位置ｚは角周波数ω
に換算する。次に、シンク（sinc）パルスをフーリエ変
換した図８に示す如き基本関数Ｆ（ω）に前記変化関数
Ｒ（ω）を乗算し、その乗算結果Ｆ（ω）・Ｒ（ω）を
ＲＦパルスのプロファイルとして決定する。ステップＶ
７では、前記決定したプロファイルをフーリエ逆変換し
てＲＦパルスを算出し、そのＲＦパルスを励起スラブ領
域Ｓに照射し、３次元ＴＯＦアンギオグラフィのパルス
シーケンスにより、励起スラブ領域Ｓ内の血流情報を収
集する。
度を前記血流速度Ｖｏとして前記条件式を解いて、図７
に示すように、フリップ角θｉを励起スラブ領域Ｓの厚
さ方向に変化させる傾き関数Ｒ（ω）を算出する。な
お、励起スラブ領域Ｓの厚さ方向の位置ｚは角周波数ω
に換算する。次に、シンク（sinc）パルスをフーリエ変
換した図８に示す如き基本関数Ｆ（ω）に前記変化関数
Ｒ（ω）を乗算し、その乗算結果Ｆ（ω）・Ｒ（ω）を
ＲＦパルスのプロファイルとして決定する。ステップＶ
７では、前記決定したプロファイルをフーリエ逆変換し
てＲＦパルスを算出し、そのＲＦパルスを励起スラブ領
域Ｓに照射し、３次元ＴＯＦアンギオグラフィのパルス
シーケンスにより、励起スラブ領域Ｓ内の血流情報を収
集する。
【００１５】以上のＭＲＩ装置１００によれば、励起ス
ラブ領域Ｓ内の血流の速度分布を実測し、その速度分布
中から対象血流速度を選択し、その対象血流速度に対し
てＲＦパルスのプロファイルを最適化するから、被検体
Ｈや撮像部位により支配的な血流速度が変わっても、常
に最適のプロファイルのＲＦパルスを用いることが出来
るようになる。
ラブ領域Ｓ内の血流の速度分布を実測し、その速度分布
中から対象血流速度を選択し、その対象血流速度に対し
てＲＦパルスのプロファイルを最適化するから、被検体
Ｈや撮像部位により支配的な血流速度が変わっても、常
に最適のプロファイルのＲＦパルスを用いることが出来
るようになる。
【００１６】なお、上記実施例では操作者が所望の血流
速度を対象血流速度として設定したが、ヒストグラムの
最も頻度の高い血流速度を自動的に対象血流速度として
設定するようにしてもよい。また、上記実施例では励起
スラブ領域Ｓの全体の血流の速度分布を実測したが、励
起スラブ領域Ｓの一部（例えば操作者が指定した部分）
の血流の速度分布のみを実測してもよい。この場合は、
実測時間を短縮できる。さらに、上記実施例では、演算
によりプロファイルを最適化しＲＦパルスを算出した
が、複数種類の対象血流速度（および他のパラメータ）
とＲＦパルスの組合せを予め登録しておき、設定された
対象血流速度（および他のパラメータ）に対応するＲＦ
パルスを検索して取得し、それを用いて３次元ＴＯＦア
ンギオグラフィを行うようにしてもよい。
速度を対象血流速度として設定したが、ヒストグラムの
最も頻度の高い血流速度を自動的に対象血流速度として
設定するようにしてもよい。また、上記実施例では励起
スラブ領域Ｓの全体の血流の速度分布を実測したが、励
起スラブ領域Ｓの一部（例えば操作者が指定した部分）
の血流の速度分布のみを実測してもよい。この場合は、
実測時間を短縮できる。さらに、上記実施例では、演算
によりプロファイルを最適化しＲＦパルスを算出した
が、複数種類の対象血流速度（および他のパラメータ）
とＲＦパルスの組合せを予め登録しておき、設定された
対象血流速度（および他のパラメータ）に対応するＲＦ
パルスを検索して取得し、それを用いて３次元ＴＯＦア
ンギオグラフィを行うようにしてもよい。
【００１７】
【発明の効果】この発明のＭＲ血流情報収集方法および
ＭＲＩ装置によれば、励起スラブ領域内の血流速度分布
を実測し、その血流速度分布中から選択した対象血流速
度に対してＲＦパルスのプロファイルを最適化するか
ら、被検体や撮像部位により血流速度が変わっても、常
に最適のプロファイルのＲＦパルスを用いて血流情報を
収集できるようになる。
ＭＲＩ装置によれば、励起スラブ領域内の血流速度分布
を実測し、その血流速度分布中から選択した対象血流速
度に対してＲＦパルスのプロファイルを最適化するか
ら、被検体や撮像部位により血流速度が変わっても、常
に最適のプロファイルのＲＦパルスを用いて血流情報を
収集できるようになる。
【図１】この発明の一実施例のＭＲＩ装置を示すブロッ
ク図である。
ク図である。
【図２】図１のＭＲＩ装置の動作を説明するフローチャ
ートである。
ートである。
【図３】励起スラブ領域とスライスの概念図である。
【図４】スライス面内の血流の２次元分布の概念図であ
る。
る。
【図５】励起スラブ領域内での血流速度のヒストグラム
である。
である。
【図６】対象血流速度の設定を説明する概念図である。
【図７】傾き関数の概念図である。
【図８】基本関数の概念図である。
【図９】フリップ角がランプ状に変化するプロファイル
を持つＲＦパルスを照射する３次元ＴＯＦアンギオグラ
フィの説明図である。
を持つＲＦパルスを照射する３次元ＴＯＦアンギオグラ
フィの説明図である。
【図１０】フリップ角をランプ状に変化させる理由の説
明図である。
明図である。
１００ ＭＲ装置 １ マグネットアセンブリ ３ 勾配磁場駆動回路 ７ 計算機 ８ シーケンス記憶回路 ９ ゲート変調回路 １０ ＲＦ発振回路 １３ 操作卓 Ｓ 励起スラブ領域
フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 ＦＩ 技術表示箇所 Ｇ０６Ｆ 15/62 ３９０ Ｃ
Claims (5)
- 【請求項１】 励起スラブ領域の厚さ方向にフリップ角
が徐々に変化するようなプロファイルを持つＲＦパルス
を照射してその励起スラブ領域内の血流情報を収集する
ＭＲ血流情報収集方法において、 前記励起スラブ領域内の血流の厚さ方向の速度分布を求
め、その速度分布中に含まれる速度を対象血流速度とし
て選択し、その対象血流速度に対して前記ＲＦパルスの
プロファイルを最適化することを特徴とするＭＲ血流情
報収集方法。 - 【請求項２】 励起スラブ領域の厚さ方向にフリップ角
が徐々に変化するようなプロファイルを持つＲＦパルス
を照射してその励起スラブ領域内の血流情報を収集する
ＭＲＩ装置において、 前記励起スラブ領域内の血流の厚さ方向の速度分布を取
得する速度分布取得手段と、その速度分布中に含まれる
速度を対象血流速度として選択する対象血流速度選択手
段と、その対象血流速度に対して前記ＲＦパルスのプロ
ファイルを最適化するプロファイル最適化手段とを具備
したことを特徴とするＭＲＩ装置。 - 【請求項３】 請求項２に記載のＭＲＩ装置において、
前記速度分布取得手段は、前記励起スラブ領域を複数の
スライスに分割し、位相コントラスト法により各スライ
スについてスライス厚方向速度のスライス面内での２次
元分布を取得し、全てのスライスの前記２次元分布を合
算して、前記速度分布を取得することを特徴とするＭＲ
Ｉ装置。 - 【請求項４】 請求項２または請求項３に記載のＭＲＩ
装置において、前記対象血流速度選択手段は、前記速度
分布を表示し、操作者に対象血流速度を設定させること
を特徴とするＭＲＩ装置。 - 【請求項５】 請求項２から請求項４のいずれかに記載
のＭＲＩ装置において、前記プロファイル最適化手段
は、励起スラブ領域の厚さ方向にフリップ角を変化させ
る変化関数を前記対象血流速度から算出し、シンクパル
スをフーリエ変換した基本関数に前記変化関数を乗算し
てＲＦパルスのプロファイルを決定することを特徴とす
るＭＲＩ装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6059266A JPH07265277A (ja) | 1994-03-29 | 1994-03-29 | Ｍｒ血流情報収集方法およびｍｒｉ装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6059266A JPH07265277A (ja) | 1994-03-29 | 1994-03-29 | Ｍｒ血流情報収集方法およびｍｒｉ装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH07265277A true JPH07265277A (ja) | 1995-10-17 |
Family
ID=13108407
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP6059266A Pending JPH07265277A (ja) | 1994-03-29 | 1994-03-29 | Ｍｒ血流情報収集方法およびｍｒｉ装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH07265277A (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002079791A1 (en) * | 2001-03-30 | 2002-10-10 | Koninkliijke Philips Electronics Nv | Flow adaptive guided determination of imaging parameters for 3d magnetic resonance angiography |
| JP2007029763A (ja) * | 2006-11-13 | 2007-02-08 | Toshiba Corp | 磁気共鳴映像化装置 |
| JP2008278961A (ja) * | 2007-05-08 | 2008-11-20 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置およびスラブ領域設定方法 |
| WO2016003070A1 (en) * | 2014-06-30 | 2016-01-07 | Samsung Electronics Co., Ltd. | Method of measuring blood flow velocity performed by medical imaging apparatus, and the medical imaging apparatus |
| KR20160002294A (ko) * | 2014-06-30 | 2016-01-07 | 삼성전자주식회사 | 의료 영상 장치에서의 혈류 속도 산출 방법 및 그 의료 영상 장치 |
| CN106539585A (zh) * | 2015-09-18 | 2017-03-29 | 三星电子株式会社 | 扫描血管的磁共振成像设备和方法 |
-
1994
- 1994-03-29 JP JP6059266A patent/JPH07265277A/ja active Pending
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002079791A1 (en) * | 2001-03-30 | 2002-10-10 | Koninkliijke Philips Electronics Nv | Flow adaptive guided determination of imaging parameters for 3d magnetic resonance angiography |
| JP2007029763A (ja) * | 2006-11-13 | 2007-02-08 | Toshiba Corp | 磁気共鳴映像化装置 |
| JP2008278961A (ja) * | 2007-05-08 | 2008-11-20 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置およびスラブ領域設定方法 |
| WO2016003070A1 (en) * | 2014-06-30 | 2016-01-07 | Samsung Electronics Co., Ltd. | Method of measuring blood flow velocity performed by medical imaging apparatus, and the medical imaging apparatus |
| KR20160002294A (ko) * | 2014-06-30 | 2016-01-07 | 삼성전자주식회사 | 의료 영상 장치에서의 혈류 속도 산출 방법 및 그 의료 영상 장치 |
| CN106539585A (zh) * | 2015-09-18 | 2017-03-29 | 三星电子株式会社 | 扫描血管的磁共振成像设备和方法 |
| EP3144691A3 (en) * | 2015-09-18 | 2017-06-14 | Samsung Electronics Co., Ltd. | Magnetic resonance imaging apparatus and method of scanning blood vessel |
| US10151819B2 (en) | 2015-09-18 | 2018-12-11 | Samsung Electronics Co., Ltd. | Magnetic resonance imaging apparatus and method of scanning blood vessel |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4471305A (en) | Method and apparatus for rapid NMR imaging of nuclear parameters with an object | |
| JP4141147B2 (ja) | 磁気共鳴イメージング装置 | |
| JP4416221B2 (ja) | 磁気共鳴画像診断装置 | |
| USRE33259E (en) | Method and apparatus for rapid NMR imaging of nuclear parameters with an object | |
| US6618609B2 (en) | Flow adaptive guided determination of imaging parameters for 3D magnetic resonance angiography | |
| JPH01101970A (ja) | 磁気共鳴イメージング装置 | |
| JPH05253208A (ja) | 磁気共鳴イメージング装置 | |
| JP3538249B2 (ja) | Ｍｒｉｓ装置 | |
| JP3987223B2 (ja) | 核磁気共鳴イメージング装置 | |
| JP3869056B2 (ja) | 医用画像診断装置における画像表示方法 | |
| JP3137366B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3519794B2 (ja) | Ｓａｔ設定方法及びｓａｔ設定装置 | |
| JP4229254B2 (ja) | 磁気共鳴イメージング装置 | |
| JP4297541B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3884282B2 (ja) | Ｍｒｉ装置 | |
| JP3526347B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3439845B2 (ja) | Ｍｒｉ用ランプドｒｆパルスの生成方法およびｍｒｉ装置 | |
| CN112824922B (zh) | 磁共振拍摄装置及其控制方法 | |
| JP3447099B2 (ja) | Ｍｒｉ装置 | |
| JPH08591A (ja) | 磁気共鳴イメージング装置 | |
| JP3454865B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3124789B2 (ja) | Ｍｒイメージング方式 | |
| JPH05103768A (ja) | 磁気共鳴イメージング方法 | |
| JP4067207B2 (ja) | 核磁気共鳴撮像装置 | |
| JP2000316831A (ja) | 磁気共鳴イメージング装置 |