JP2000279390A - 磁気共鳴イメージング装置 - Google Patents
磁気共鳴イメージング装置Info
- Publication number
- JP2000279390A JP2000279390A JP11090128A JP9012899A JP2000279390A JP 2000279390 A JP2000279390 A JP 2000279390A JP 11090128 A JP11090128 A JP 11090128A JP 9012899 A JP9012899 A JP 9012899A JP 2000279390 A JP2000279390 A JP 2000279390A
- Authority
- JP
- Japan
- Prior art keywords
- diffusion
- image
- picture
- signal
- magnetic field
- 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)
Abstract
(57)【要約】
【課題】拡散強調画像を用いた計算処理の画質の劣化を
低減する。 【解決手段】拡散傾斜磁場の強度、印加時間等で決まる
b値を変えて複数の拡散画像データを得、これら複数の
拡散画像データから見かけの拡散係数等を求める拡散画
像処理において、複数の拡散画像データをそれぞれ基準
となる画像のパラメータと比較し画質の劣化を判定す
る。基準となる画像は、拡散傾斜磁場をかけていない画
像とする。判定の結果、画質の劣化が許容値以内である
画像データのみを用いて拡散画像処理計算を行う。高画
質な画像から拡散計算処理を行うことが可能となり、計
算処理後の画像の画質が向上する。
低減する。 【解決手段】拡散傾斜磁場の強度、印加時間等で決まる
b値を変えて複数の拡散画像データを得、これら複数の
拡散画像データから見かけの拡散係数等を求める拡散画
像処理において、複数の拡散画像データをそれぞれ基準
となる画像のパラメータと比較し画質の劣化を判定す
る。基準となる画像は、拡散傾斜磁場をかけていない画
像とする。判定の結果、画質の劣化が許容値以内である
画像データのみを用いて拡散画像処理計算を行う。高画
質な画像から拡散計算処理を行うことが可能となり、計
算処理後の画像の画質が向上する。
Description
【０００１】
【発明が属する技術分野】本発明は、被検体中の水素や
燐等からの核磁気共鳴(以下、「NMR」という)信号を測
定し、核の密度分布や緩和時間分布等を映像化する核磁
気共鳴撮影(MRI)装置に関し、特に拡散画像を作成する
機能を備えたMRI装置に関する。
燐等からの核磁気共鳴(以下、「NMR」という)信号を測
定し、核の密度分布や緩和時間分布等を映像化する核磁
気共鳴撮影(MRI)装置に関し、特に拡散画像を作成する
機能を備えたMRI装置に関する。
【０００２】
【従来の技術】MRI装置を用いた撮像手法の一つに拡散
イメージングがある。拡散イメージングは、磁場の不均
一を起こす強力な傾斜磁場（以下、拡散傾斜磁場とい
う）を用いることにより、分子レベルでの原子核スピン
の動き（拡散）による信号低下を強調して画像化する。
拡散傾斜磁場パルスの強度Ｇ、印加時間σおよび一対の
拡散傾斜磁場の間隔Δから計算される値（＝Ｇ2σ2γ2(
Δ-σ3）はb値（b value、b factor）と呼ばれ、この
b値が大きいほど、移動スピンの位相の変化が大きくな
り、拡散が強調される。この結果、得られる画像を拡散
強調画像という。また、異なるb値の拡散強調画像を複
数回撮影し、b値の変化と信号強度の変化から組織の見
かけの拡散係数ADC（Apparent Diffusion Coefficien
t）を求めることができ、またADCにより拡散のしやすさ
を示す拡散係数画像を算出できる。
イメージングがある。拡散イメージングは、磁場の不均
一を起こす強力な傾斜磁場（以下、拡散傾斜磁場とい
う）を用いることにより、分子レベルでの原子核スピン
の動き（拡散）による信号低下を強調して画像化する。
拡散傾斜磁場パルスの強度Ｇ、印加時間σおよび一対の
拡散傾斜磁場の間隔Δから計算される値（＝Ｇ2σ2γ2(
Δ-σ3）はb値（b value、b factor）と呼ばれ、この
b値が大きいほど、移動スピンの位相の変化が大きくな
り、拡散が強調される。この結果、得られる画像を拡散
強調画像という。また、異なるb値の拡散強調画像を複
数回撮影し、b値の変化と信号強度の変化から組織の見
かけの拡散係数ADC（Apparent Diffusion Coefficien
t）を求めることができ、またADCにより拡散のしやすさ
を示す拡散係数画像を算出できる。
【０００３】上記拡散強調イメージングでは、被検体の
わずかな動きにより画像が乱れることが知られている。
これを回避するために、被検体の動きが無視できるシン
グルショットEPIで撮影することが試みられている。し
かし、シングルショットEPIの空間分解能は必ずしも高
くはない。空間分解能を向上するにはマルチショットシ
ーケンスでの拡散強調イメージングが必須である。マル
チショットシーケンスでは、体動による取得信号の位相
乱れを補正するナビゲーションエコー法が公知である
(「正常及び虚血ヒト脳のナビゲートされた拡散イメー
ジング」Alex J. de Crespignyら、マグネチックレ
ゾナンスインメディスン33、p720-728（1995）、米国特
許4937526号、米国特許5254948号など）。
わずかな動きにより画像が乱れることが知られている。
これを回避するために、被検体の動きが無視できるシン
グルショットEPIで撮影することが試みられている。し
かし、シングルショットEPIの空間分解能は必ずしも高
くはない。空間分解能を向上するにはマルチショットシ
ーケンスでの拡散強調イメージングが必須である。マル
チショットシーケンスでは、体動による取得信号の位相
乱れを補正するナビゲーションエコー法が公知である
(「正常及び虚血ヒト脳のナビゲートされた拡散イメー
ジング」Alex J. de Crespignyら、マグネチックレ
ゾナンスインメディスン33、p720-728（1995）、米国特
許4937526号、米国特許5254948号など）。
【０００４】
【発明が解決しようとする課題】しかし拡散イメージン
グでは、ナビゲーションエコー法を用いても、体動によ
る取得信号の乱れを完全には除去できない場合がある。
これは拡散イメージングでは移動するスピンの位相をデ
ィフェイズする強い拡散傾斜磁場を印加するため、この
拡散傾斜磁場印加の間(100ms以下)に体動によるスピン
の移動があると本来的には停止しているスピンの位相が
回転し、擬似的にオフセット位相エンコードが付与され
るためと考えられる。また、強い拡散傾斜磁場を印加す
るために、傾斜磁場パルスが不安定な時にも上記と同様
の現象が生じる。さらに、撮影部位によって血流による
アーチファクトも生じる。このため、良好な画像を得ら
れるときと得られないときがあり、画質にばらつきが生
じる問題があり、これら複数の画像を利用して算出する
拡散係数などにも影響が生じる。そこで本発明は、画質
のばらつきがなく良好な拡散画像を得ることができるMR
I装置を提供することを目的とする。また精度よく拡散
係数を算出することができるMRI装置を目的とする。
グでは、ナビゲーションエコー法を用いても、体動によ
る取得信号の乱れを完全には除去できない場合がある。
これは拡散イメージングでは移動するスピンの位相をデ
ィフェイズする強い拡散傾斜磁場を印加するため、この
拡散傾斜磁場印加の間(100ms以下)に体動によるスピン
の移動があると本来的には停止しているスピンの位相が
回転し、擬似的にオフセット位相エンコードが付与され
るためと考えられる。また、強い拡散傾斜磁場を印加す
るために、傾斜磁場パルスが不安定な時にも上記と同様
の現象が生じる。さらに、撮影部位によって血流による
アーチファクトも生じる。このため、良好な画像を得ら
れるときと得られないときがあり、画質にばらつきが生
じる問題があり、これら複数の画像を利用して算出する
拡散係数などにも影響が生じる。そこで本発明は、画質
のばらつきがなく良好な拡散画像を得ることができるMR
I装置を提供することを目的とする。また精度よく拡散
係数を算出することができるMRI装置を目的とする。
【０００５】
【課題を解決するための手段】上記課題を解決するため
に、本発明のMRI装置は、被検体の組織を構成する原子
核スピンを励起する高周波パルスを印加する手段と、前
記原子核スピンのうち移動するスピンの位相をディフェ
イズする拡散傾斜磁場を印加する手段と、前記被検体か
ら発生する複数のエコー信号を計測する手段と、前記エ
コー信号に基づき前記移動スピンからの信号を抑制した
拡散画像を形成する手段と、前記拡散画像を用いて拡散
係数計算等の拡散画像処理を行なう計算手段とを備え、
前記計算手段は、複数の拡散画像の画質劣化を判定する
手段を備え、前記判定手段において画質劣化が許容値以
内と判定された画像を用いて拡散画像処理を行なうもの
である。ここで拡散画像処理とは拡散係数や非等方性拡
散の算出、拡散係数画像の作成を含む。
に、本発明のMRI装置は、被検体の組織を構成する原子
核スピンを励起する高周波パルスを印加する手段と、前
記原子核スピンのうち移動するスピンの位相をディフェ
イズする拡散傾斜磁場を印加する手段と、前記被検体か
ら発生する複数のエコー信号を計測する手段と、前記エ
コー信号に基づき前記移動スピンからの信号を抑制した
拡散画像を形成する手段と、前記拡散画像を用いて拡散
係数計算等の拡散画像処理を行なう計算手段とを備え、
前記計算手段は、複数の拡散画像の画質劣化を判定する
手段を備え、前記判定手段において画質劣化が許容値以
内と判定された画像を用いて拡散画像処理を行なうもの
である。ここで拡散画像処理とは拡散係数や非等方性拡
散の算出、拡散係数画像の作成を含む。
【０００６】本発明のMRI装置は、判定手段において拡
散係数計算等の拡散画像処理の基礎となる拡散画像の劣
化を判定し、劣化していない画像のみを用いて計算を行
なうので精度よく拡散係数や非等方性拡散等を求めるこ
とができ、また画質のよい拡散係数画像を得ることがで
きる。画質劣化の判定は、好適には拡散傾斜磁場をかけ
ずに取得した画像の背景の画素値を基準として、拡散画
像の背景の画素値を比較することにより行う。
散係数計算等の拡散画像処理の基礎となる拡散画像の劣
化を判定し、劣化していない画像のみを用いて計算を行
なうので精度よく拡散係数や非等方性拡散等を求めるこ
とができ、また画質のよい拡散係数画像を得ることがで
きる。画質劣化の判定は、好適には拡散傾斜磁場をかけ
ずに取得した画像の背景の画素値を基準として、拡散画
像の背景の画素値を比較することにより行う。
【０００７】
【発明の実施の形態】以下、本発明を図面に示す実施例
を参照して詳述する。図３は、本発明が適用される典型
的なMRI装置の構成を示す図で、このMRI装置は、被検体
401が置かれる空間に静磁場を発生する磁石402と、この
空間に傾斜磁場を発生する傾斜磁場コイル403と、この
領域に高周波磁場を発生するRFコイル404と、被検体401
が発生するNMR信号を検出するRFプローブ405とを備えて
いる。また被検体401を測定空間に搬送するためのベッ
ド412が備えられている。
を参照して詳述する。図３は、本発明が適用される典型
的なMRI装置の構成を示す図で、このMRI装置は、被検体
401が置かれる空間に静磁場を発生する磁石402と、この
空間に傾斜磁場を発生する傾斜磁場コイル403と、この
領域に高周波磁場を発生するRFコイル404と、被検体401
が発生するNMR信号を検出するRFプローブ405とを備えて
いる。また被検体401を測定空間に搬送するためのベッ
ド412が備えられている。
【０００８】傾斜磁場コイル403は、X,Y,Zの3方向の傾
斜磁場コイルで構成され、傾斜磁場電源409からの信号
に応じてそれぞれ傾斜磁場を発生する。この傾斜磁場の
かけ方によって被検体の断面を選択し、またNMR信号を
位相エンコードすることができる。さらに傾斜磁場コイ
ル403は、拡散イメージングにおいて移動スピンの位相
をディフェイズするための拡散傾斜磁場を発生する。RF
コイル404はRF送信部410の信号に応じて高周波磁場を発
生する。RFプローブ405が受信した信号は、信号検出部4
06で検出され、信号処理部407で信号処理され、また計
算により画像信号に変換される。信号処理部407は、拡
散傾斜磁場を用いた撮影の場合には拡散画像の再構成の
他、複数の拡散画像を用いて拡散係数および拡散係数画
像の算出を行う。画像は表示部408で表示される。傾斜
磁場電源409、RF送信部410、信号検出部406は制御部411
で制御される。制御のタイムチャートは一般にパルスシ
ーケンスと呼ばれ、拡散画像の撮影においては図２に示
すようなエコープレナー（EPI）法を基本とするパルス
シーケンスを実行する。
斜磁場コイルで構成され、傾斜磁場電源409からの信号
に応じてそれぞれ傾斜磁場を発生する。この傾斜磁場の
かけ方によって被検体の断面を選択し、またNMR信号を
位相エンコードすることができる。さらに傾斜磁場コイ
ル403は、拡散イメージングにおいて移動スピンの位相
をディフェイズするための拡散傾斜磁場を発生する。RF
コイル404はRF送信部410の信号に応じて高周波磁場を発
生する。RFプローブ405が受信した信号は、信号検出部4
06で検出され、信号処理部407で信号処理され、また計
算により画像信号に変換される。信号処理部407は、拡
散傾斜磁場を用いた撮影の場合には拡散画像の再構成の
他、複数の拡散画像を用いて拡散係数および拡散係数画
像の算出を行う。画像は表示部408で表示される。傾斜
磁場電源409、RF送信部410、信号検出部406は制御部411
で制御される。制御のタイムチャートは一般にパルスシ
ーケンスと呼ばれ、拡散画像の撮影においては図２に示
すようなエコープレナー（EPI）法を基本とするパルス
シーケンスを実行する。
【０００９】次に上記構成のMRI装置を用いた拡散イメ
ージングの手順、特に信号処理部407における処理手順
を図１を参照して説明する。まず拡散イメージングのた
めのパルスシーケンスを実行し、エコー信号(MR信号)を
計測する(101)。パルスシーケンスは例えば図２に示す
ようなマルチショット拡散EPIが採用される。図示する
パルスシーケンスは4ショットSE EPIシーケンスを基本
としており、90°パルス201と180°パルス202の照射
後、読み出し傾斜磁場Gr203を極性を反転させながら印
加してエコー信号205を発生させるとともに位相エンコ
ード傾斜磁場Gp204を印加し、各エコー信号205に異なる
位相エンコードを付与する。位相エンコード数を256と
した場合、１ショットで計測されるエコー信号205は64
である。また180°パルス202の前後に拡散強調傾斜磁場
パルス206が付加される。図では読み出し傾斜磁場Gr、
位相エンコード傾斜磁場Gp、スライス傾斜磁場Gsの全方
向に拡散強調傾斜磁場パルスが付加してあるが、これら
の内の１方向でも良く、通常は１方向が選ばれる。
ージングの手順、特に信号処理部407における処理手順
を図１を参照して説明する。まず拡散イメージングのた
めのパルスシーケンスを実行し、エコー信号(MR信号)を
計測する(101)。パルスシーケンスは例えば図２に示す
ようなマルチショット拡散EPIが採用される。図示する
パルスシーケンスは4ショットSE EPIシーケンスを基本
としており、90°パルス201と180°パルス202の照射
後、読み出し傾斜磁場Gr203を極性を反転させながら印
加してエコー信号205を発生させるとともに位相エンコ
ード傾斜磁場Gp204を印加し、各エコー信号205に異なる
位相エンコードを付与する。位相エンコード数を256と
した場合、１ショットで計測されるエコー信号205は64
である。また180°パルス202の前後に拡散強調傾斜磁場
パルス206が付加される。図では読み出し傾斜磁場Gr、
位相エンコード傾斜磁場Gp、スライス傾斜磁場Gsの全方
向に拡散強調傾斜磁場パルスが付加してあるが、これら
の内の１方向でも良く、通常は１方向が選ばれる。
【００１０】このようなシーケンスを4回繰り返すこと
により、１枚の画像再構成に必要なエコー信号（本計測
データ）が得られる。このシーケンスでは一対の強い拡
散傾斜磁場206を印加することによって、静止している
プロトンでは位相がリフェイズするのに対し、拡散して
いるプロトンでは位相がディフェイズし、それらからの
信号が低下する。このようなエコー信号を用いて画像再
構成を行うことにより移動プロトンからの信号を抑制し
た拡散強調画像が得られる。このようなステップ101をb
値を変えて複数回（例えば5回）実行する。既に述べた
ようにb値は、拡散傾斜磁場パルスの強度、印加時間お
よび２つの拡散傾斜磁場パルスの間隔によって決まるの
で、これらのファクタのいずれか１つ以上を変化させて
撮影することによりb値の異なる複数の拡散強調画像を
得ることができる。
により、１枚の画像再構成に必要なエコー信号（本計測
データ）が得られる。このシーケンスでは一対の強い拡
散傾斜磁場206を印加することによって、静止している
プロトンでは位相がリフェイズするのに対し、拡散して
いるプロトンでは位相がディフェイズし、それらからの
信号が低下する。このようなエコー信号を用いて画像再
構成を行うことにより移動プロトンからの信号を抑制し
た拡散強調画像が得られる。このようなステップ101をb
値を変えて複数回（例えば5回）実行する。既に述べた
ようにb値は、拡散傾斜磁場パルスの強度、印加時間お
よび２つの拡散傾斜磁場パルスの間隔によって決まるの
で、これらのファクタのいずれか１つ以上を変化させて
撮影することによりb値の異なる複数の拡散強調画像を
得ることができる。
【００１１】次に、複数枚のb値の異なる拡散強調画像
それぞれに対して拡散係数算出などに使用できるデータ
か否かの判定を行う（102）。判定ではまず任意の基準
画像を決めておき、基準画像と拡散強調画像とを比較す
る。基準画像としては複数の拡散強調画像の一つを用い
ることも可能であるが、好適には拡散傾斜磁場206をか
けずに撮影した画像、即ちb値が0の画像を用いる。これ
は強い拡散傾斜磁場を印加していないためにアーチファ
クトが拡散強調画像よりも少ないためである。
それぞれに対して拡散係数算出などに使用できるデータ
か否かの判定を行う（102）。判定ではまず任意の基準
画像を決めておき、基準画像と拡散強調画像とを比較す
る。基準画像としては複数の拡散強調画像の一つを用い
ることも可能であるが、好適には拡散傾斜磁場206をか
けずに撮影した画像、即ちb値が0の画像を用いる。これ
は強い拡散傾斜磁場を印加していないためにアーチファ
クトが拡散強調画像よりも少ないためである。
【００１２】画像と画像の比較は背景部分を用いて、そ
の部分の画素値を比較することにより行なう。比較する
背景部分105は流れアーチファクトが生じやすい撮影部
位の位相エンコード方向(Gp)の外側に設定するのが望ま
しい。例えば図４に示すように、被検体の頭部や腹部な
どの撮影部位の外側に存在する背景部分105で行なう。
画素値は、ある面積（例えば20×20ピクセル）を持つ任
意の背景部分の画素値の合計値、平均値或いは標準偏差
をパラメータとして基準画像のパラメータと比較する。
ある面積を待った背景部分を用いて判定することによ
り、１点の画素値のみで判定する場合に比べ、画素値の
ばらつきの影響を受けることなく正確な判定ができる。
の部分の画素値を比較することにより行なう。比較する
背景部分105は流れアーチファクトが生じやすい撮影部
位の位相エンコード方向(Gp)の外側に設定するのが望ま
しい。例えば図４に示すように、被検体の頭部や腹部な
どの撮影部位の外側に存在する背景部分105で行なう。
画素値は、ある面積（例えば20×20ピクセル）を持つ任
意の背景部分の画素値の合計値、平均値或いは標準偏差
をパラメータとして基準画像のパラメータと比較する。
ある面積を待った背景部分を用いて判定することによ
り、１点の画素値のみで判定する場合に比べ、画素値の
ばらつきの影響を受けることなく正確な判定ができる。
【００１３】具体的には、拡散強調画像の背景部分の所
定面積の画素値のパラメータＰ（合計値、平均値、標準
偏差など）が、基準画像のパラメータのＮ倍（Ｐ0×
Ｎ）よりも小さければ、その拡散強調画像の画質劣化は
許容値以内と判定し、基準画像のパラメータのＮ倍（Ｐ
0×Ｎ）よりも大きければ、許容値を超えていると判定
する。ここで、Ｎは許容値以内かどうか判定するための
係数である。通常、背景の信号値は０であるのに対し、
拡散強調画像のアーチファクトが存在すると信号値は０
よりも大きくなるので、Ｎの値としては１以上の値を設
定する。例えば2.0とする。このような判定の後、信号
処理部407は許容値以内であると判定された画像のみを
用いて拡散強調計算処理を行う（103）。画質劣化が許
容値を超える画像データは計算には使用しない。
定面積の画素値のパラメータＰ（合計値、平均値、標準
偏差など）が、基準画像のパラメータのＮ倍（Ｐ0×
Ｎ）よりも小さければ、その拡散強調画像の画質劣化は
許容値以内と判定し、基準画像のパラメータのＮ倍（Ｐ
0×Ｎ）よりも大きければ、許容値を超えていると判定
する。ここで、Ｎは許容値以内かどうか判定するための
係数である。通常、背景の信号値は０であるのに対し、
拡散強調画像のアーチファクトが存在すると信号値は０
よりも大きくなるので、Ｎの値としては１以上の値を設
定する。例えば2.0とする。このような判定の後、信号
処理部407は許容値以内であると判定された画像のみを
用いて拡散強調計算処理を行う（103）。画質劣化が許
容値を超える画像データは計算には使用しない。
【００１４】拡散強調計算処理では、b値の異なる複数
の画像の画素の信号強度から見かけの拡散係数（ADC）
を算出し、ADC計算画像を作成する。一般に画質劣化が
ない場合、拡散による信号低下は下式に示すように指数
関数的に減衰する。 Ｓ＝Ｓ0・exp（−b・D） 式中、Ｓは信号値、Ｓ0はb値が０の場合の信号値、Ｄは
拡散係数である。従ってＳとb値が与えられればＤを計
算によって求めることができる。具体的には異なるb値
の拡散傾斜磁場を印加して得られた複数の拡散画像から
b値と信号値との関係を求め、これを指数フィッティン
グすることにより見かけ拡散係数Ｄを求める。
の画像の画素の信号強度から見かけの拡散係数（ADC）
を算出し、ADC計算画像を作成する。一般に画質劣化が
ない場合、拡散による信号低下は下式に示すように指数
関数的に減衰する。 Ｓ＝Ｓ0・exp（−b・D） 式中、Ｓは信号値、Ｓ0はb値が０の場合の信号値、Ｄは
拡散係数である。従ってＳとb値が与えられればＤを計
算によって求めることができる。具体的には異なるb値
の拡散傾斜磁場を印加して得られた複数の拡散画像から
b値と信号値との関係を求め、これを指数フィッティン
グすることにより見かけ拡散係数Ｄを求める。
【００１５】図５は、５つの異なるb値の拡散傾斜磁場
を印加して拡散画像を撮影した場合におけるb値と信号
強度の関係をグラフにしたものである。ここでは、b値
がb2である画像（信号強度112に対応）は、画質の判定
において許容値以上と判定されたものとする。これをそ
のまま指数フィッティングすると、曲線121となる。し
かしここでは、画質の判定において許容値以内と判定さ
れたもののみ、即ち信号強度112を除いた4つの信号強度
（111,113,114,115）を用いて指数フィッティングす
る。その結果、曲線120となる。この曲線120は曲線121
に比べ、4つの信号強度にうまくフィッティングしてい
ることが分かる。これによりADC計算画像もより正確な
ものとなる。
を印加して拡散画像を撮影した場合におけるb値と信号
強度の関係をグラフにしたものである。ここでは、b値
がb2である画像（信号強度112に対応）は、画質の判定
において許容値以上と判定されたものとする。これをそ
のまま指数フィッティングすると、曲線121となる。し
かしここでは、画質の判定において許容値以内と判定さ
れたもののみ、即ち信号強度112を除いた4つの信号強度
（111,113,114,115）を用いて指数フィッティングす
る。その結果、曲線120となる。この曲線120は曲線121
に比べ、4つの信号強度にうまくフィッティングしてい
ることが分かる。これによりADC計算画像もより正確な
ものとなる。
【００１６】尚、以上の実施例ではステップ102で許容
値以上と判定された画像データを使用せずに、許容値以
内の画像データのみを用いて拡散強調計算処理（103）
を行う場合を説明したが、判定結果に基づき許容値以上
と判定された画像についてステップ101を再実行し、新
たな画像データを得ることも可能である。その場合、例
えばステップ102とステップ103との間で許容値以内と判
定された画像データ数をカウントするステップを挿入
し、その数が一定値以下の場合にのみ許容値以上のデー
タについてステップ101に戻り、新たにデータ取得（10
1）、画質の判定（102）を行うようにしてもよい。これ
によりデータ数の減少による計算精度の低下を防止する
ことができる。また図１の手順では、拡散強調画像エコ
ー計測と画像再構成から、拡散強調計算処理まで一連の
処理として行う場合を説明したが、拡散強調画像エコー
計測と画像再構成はMRI装置で行ない、データ採否の決
定と拡散強調計算処理をMRI装置にイーサネット（登録
商標）などでつながっているコンピュータに画像を転送
して行なうことも可能である。
値以上と判定された画像データを使用せずに、許容値以
内の画像データのみを用いて拡散強調計算処理（103）
を行う場合を説明したが、判定結果に基づき許容値以上
と判定された画像についてステップ101を再実行し、新
たな画像データを得ることも可能である。その場合、例
えばステップ102とステップ103との間で許容値以内と判
定された画像データ数をカウントするステップを挿入
し、その数が一定値以下の場合にのみ許容値以上のデー
タについてステップ101に戻り、新たにデータ取得（10
1）、画質の判定（102）を行うようにしてもよい。これ
によりデータ数の減少による計算精度の低下を防止する
ことができる。また図１の手順では、拡散強調画像エコ
ー計測と画像再構成から、拡散強調計算処理まで一連の
処理として行う場合を説明したが、拡散強調画像エコー
計測と画像再構成はMRI装置で行ない、データ採否の決
定と拡散強調計算処理をMRI装置にイーサネット（登録
商標）などでつながっているコンピュータに画像を転送
して行なうことも可能である。
【００１７】
【発明の効果】本発明によれば複数の拡散画像のうち画
質の劣化が許容値以内であると判定されたもののみを用
いて拡散計算処理を行うようにしたもので、精度の高い
拡散係数が可能であり、計算処理後の画像の画質が向上
する。
質の劣化が許容値以内であると判定されたもののみを用
いて拡散計算処理を行うようにしたもので、精度の高い
拡散係数が可能であり、計算処理後の画像の画質が向上
する。
【図１】本発明の拡散イメージングの一実施例を示すブ
ロック図。
ロック図。
【図２】本発明の拡散イメージングで採用するパルスシ
ーケンスの一実施例を示すタイミング図。
ーケンスの一実施例を示すタイミング図。
【図３】本発明が適用されるMRI装置の構成図。
【図４】画質の判定手法を説明する図。
【図５】拡散係数の算出方法を説明する図。
403・・・・・・傾斜磁場印加手段 404・・・・・・高周波パルス印加手段 405・・・・・・エコー信号計測手段 407・・・・・・計算手段 411・・・・・・画像形成手段
Claims (1)
- 【請求項１】 検体の組織を構成する原子核スピンを励
起する高周波パルスを印加する手段と、前記原子核スピ
ンのうち移動するスピンの位相をディフェイズする拡散
傾斜磁場を印加する手段と、前記被検体から発生する複
数のエコー信号を計測する手段と、前記エコー信号に基
づき前記移動スピンからの信号を抑制した拡散画像を形
成する手段と、前記拡散画像を用いて拡散係数計算等の
拡散画像処理を行なう計算手段とを備えた磁気共鳴イメ
ージング装置において、 前記計算手段は、複数の拡散画像の画質劣化を判定する
手段を備え、前記判定手段において画質劣化が許容値以
内と判定された画像のみを用いて拡散画像処理を行なう
ことを特徴とする磁気共鳴イメージング装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11090128A JP2000279390A (ja) | 1999-03-30 | 1999-03-30 | 磁気共鳴イメージング装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11090128A JP2000279390A (ja) | 1999-03-30 | 1999-03-30 | 磁気共鳴イメージング装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2000279390A true JP2000279390A (ja) | 2000-10-10 |
Family
ID=13989882
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP11090128A Pending JP2000279390A (ja) | 1999-03-30 | 1999-03-30 | 磁気共鳴イメージング装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2000279390A (ja) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005525208A (ja) * | 2002-05-15 | 2005-08-25 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Ｄｔｉｓｎｒを改善するための画像の遡及的な選択及び様々なタイプの画像アラインメント |
| JP2008183170A (ja) * | 2007-01-30 | 2008-08-14 | Ge Medical Systems Global Technology Co Llc | Ｍｒｉ装置 |
| JP2008284225A (ja) * | 2007-05-18 | 2008-11-27 | Toshiba Corp | 磁気共鳴画像診断装置 |
| JP2010042257A (ja) * | 2008-08-11 | 2010-02-25 | Siemens Medical Solutions Usa Inc | 灌流画像を生成する磁気共鳴方法及び装置 |
| JP2010527740A (ja) * | 2007-05-31 | 2010-08-19 | シーアール ディベロップメント アーベー | 拡散磁気共鳴画像化のための方法とシステム |
| JP2012529328A (ja) * | 2009-06-09 | 2012-11-22 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 記憶された画像を順序付けする装置及び方法 |
| JP2014210099A (ja) * | 2013-04-19 | 2014-11-13 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴装置およびプログラム |
| US9201129B2 (en) | 2006-09-13 | 2015-12-01 | Kabushiki Kaisha Toshiba | Magnetic-resonance image diagnostic apparatus and method of controlling the same |
| JP2017225688A (ja) * | 2016-06-23 | 2017-12-28 | 株式会社日立製作所 | 磁気共鳴イメージング装置、画像処理装置、及び拡散強調画像計算方法 |
| CN108490377A (zh) * | 2018-04-03 | 2018-09-04 | 上海东软医疗科技有限公司 | 一种磁共振弥散定量的获取方法和装置 |
-
1999
- 1999-03-30 JP JP11090128A patent/JP2000279390A/ja active Pending
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005525208A (ja) * | 2002-05-15 | 2005-08-25 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Ｄｔｉｓｎｒを改善するための画像の遡及的な選択及び様々なタイプの画像アラインメント |
| US9201129B2 (en) | 2006-09-13 | 2015-12-01 | Kabushiki Kaisha Toshiba | Magnetic-resonance image diagnostic apparatus and method of controlling the same |
| JP2008183170A (ja) * | 2007-01-30 | 2008-08-14 | Ge Medical Systems Global Technology Co Llc | Ｍｒｉ装置 |
| JP2008284225A (ja) * | 2007-05-18 | 2008-11-27 | Toshiba Corp | 磁気共鳴画像診断装置 |
| JP2010527740A (ja) * | 2007-05-31 | 2010-08-19 | シーアール ディベロップメント アーベー | 拡散磁気共鳴画像化のための方法とシステム |
| US8565854B2 (en) | 2007-05-31 | 2013-10-22 | Cr Development Ab | Method and system for diffusion magnetic resonance imaging |
| JP2010042257A (ja) * | 2008-08-11 | 2010-02-25 | Siemens Medical Solutions Usa Inc | 灌流画像を生成する磁気共鳴方法及び装置 |
| JP2012529328A (ja) * | 2009-06-09 | 2012-11-22 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 記憶された画像を順序付けする装置及び方法 |
| JP2014210099A (ja) * | 2013-04-19 | 2014-11-13 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴装置およびプログラム |
| JP2017225688A (ja) * | 2016-06-23 | 2017-12-28 | 株式会社日立製作所 | 磁気共鳴イメージング装置、画像処理装置、及び拡散強調画像計算方法 |
| WO2017221654A1 (ja) * | 2016-06-23 | 2017-12-28 | 株式会社日立製作所 | 磁気共鳴イメージング装置、画像処理装置、及び拡散強調画像計算方法 |
| CN108490377A (zh) * | 2018-04-03 | 2018-09-04 | 上海东软医疗科技有限公司 | 一种磁共振弥散定量的获取方法和装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| 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 | |
| Greenman et al. | Double inversion black‐blood fast spin‐echo imaging of the human heart: a comparison between 1.5 T and 3.0 T | |
| US9513358B2 (en) | Method and apparatus for magnetic resonance imaging | |
| EP0634664B1 (en) | Fast spin echo prescan for MRI system | |
| US6369568B1 (en) | Fast spin echo phase correction for MRI system | |
| WO2001017428A1 (fr) | Dispositif d'imagerie par resonance magnetique et procede associe | |
| US20070052417A1 (en) | Method and apparatus for diffusion magnetic resonance imaging with the effects of eddy currents compensated | |
| US20170102439A1 (en) | System and Method For Gradient Measurement Using Single-Point Imaging | |
| CN110720047B (zh) | 双回波Dixon型水/脂肪分离MR成像 | |
| US6127826A (en) | EPI image based long term eddy current pre-emphasis calibration | |
| JP3964110B2 (ja) | 磁気共鳴イメージング装置 | |
| JP2000279390A (ja) | 磁気共鳴イメージング装置 | |
| US10746832B2 (en) | DTI with correction of motion-induced diffusion gradient inconsistency | |
| JPH07327960A (ja) | 磁気共鳴イメージング装置 | |
| JP2003260037A (ja) | 磁気共鳴信号獲得装置および磁気共鳴撮影装置 | |
| JP2005111059A (ja) | 磁気共鳴イメージング装置 | |
| JPH09173315A (ja) | 磁気共鳴イメージング方法 | |
| US20020050816A1 (en) | MR imaging method, phase error measuring method, and MRI apparatus | |
| US7176686B2 (en) | Parallel MR imaging method | |
| US20220057467A1 (en) | Epi mr imaging with distortion correction | |
| JP2002085376A (ja) | 核磁気共鳴イメージング装置および方法 | |
| US7148687B2 (en) | Method for the acquisition of moving objects through nuclear magnetic resonance tomography | |
| JPH04208132A (ja) | 磁気共鳴映像装置 | |
| JP3450508B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3868754B2 (ja) | 磁気共鳴イメージング装置 |