JPH0654830A - 拡散係数分布画像の表示方法 - Google Patents
拡散係数分布画像の表示方法Info
- Publication number
- JPH0654830A JPH0654830A JP4210070A JP21007092A JPH0654830A JP H0654830 A JPH0654830 A JP H0654830A JP 4210070 A JP4210070 A JP 4210070A JP 21007092 A JP21007092 A JP 21007092A JP H0654830 A JPH0654830 A JP H0654830A
- Authority
- JP
- Japan
- Prior art keywords
- diffusion coefficient
- coefficient distribution
- distribution image
- image
- displaying
- 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)【要約】
【目的】 定量的な評価が可能で、拡散係数以外の量を
同時に表示する拡散係数分布画像の表示方法を提供す
る。 【構成】 拡散係数値を画像の表示色に対応させ、他の
量を画像の濃淡値に対応させて表示し、画像の任意の領
域における拡散係数分布のヒストグラムを１次元表示す
る。ヒストグラムの横軸上の任意の範囲を指定して、範
囲内の拡散係数を有するボクセルのみを画像として表示
する。さらに、ヒストグラムの横軸上の任意の範囲を２
か所指定して、任意の領域において、２つの範囲の拡散
係数値を有するボクセル数の比率を表示する。また、異
なる拡散係数分布画像上の任意の領域において、それぞ
れの画像に対応するヒストグラムと、両者の差を表示す
る。 【効果】 拡散係数分布を１次元表示することにより、
分布状態のずれが容易に区別できる。
同時に表示する拡散係数分布画像の表示方法を提供す
る。 【構成】 拡散係数値を画像の表示色に対応させ、他の
量を画像の濃淡値に対応させて表示し、画像の任意の領
域における拡散係数分布のヒストグラムを１次元表示す
る。ヒストグラムの横軸上の任意の範囲を指定して、範
囲内の拡散係数を有するボクセルのみを画像として表示
する。さらに、ヒストグラムの横軸上の任意の範囲を２
か所指定して、任意の領域において、２つの範囲の拡散
係数値を有するボクセル数の比率を表示する。また、異
なる拡散係数分布画像上の任意の領域において、それぞ
れの画像に対応するヒストグラムと、両者の差を表示す
る。 【効果】 拡散係数分布を１次元表示することにより、
分布状態のずれが容易に区別できる。
Description
【０００１】
【産業上の利用分野】本発明は、核磁気共鳴を用いたイ
メージング装置のデータ表示方法に係り、特に拡散係数
分布画像の表示方法に係る。
メージング装置のデータ表示方法に係り、特に拡散係数
分布画像の表示方法に係る。
【０００２】
【従来の技術】従来、核磁気共鳴等により得られた拡散
係数の分布画像を表示する際には、ラジオロジー、168
(1988年)第497頁から505頁(Radiology 168,(1988)pp497
-505)に見られるように、画像の濃淡を拡散係数値に対
応づけたモノクロ画像が用いられていた。従って、これ
を生体計測データに適用すれば、拡散係数値が高い、す
なわち可動性の大きい分子からなるボクセルは明るく、
逆に拡散係数値が低い、すなわち可動性の小さい分子か
らなるボクセルは暗く表示される。例えば、毛細血管等
が発達した部分は明るく表示されるが、腫瘍などで石灰
化した部分は暗く表示される。この方法によれば、異な
る部位における分子の可動性の違いが明瞭である。
係数の分布画像を表示する際には、ラジオロジー、168
(1988年)第497頁から505頁(Radiology 168,(1988)pp497
-505)に見られるように、画像の濃淡を拡散係数値に対
応づけたモノクロ画像が用いられていた。従って、これ
を生体計測データに適用すれば、拡散係数値が高い、す
なわち可動性の大きい分子からなるボクセルは明るく、
逆に拡散係数値が低い、すなわち可動性の小さい分子か
らなるボクセルは暗く表示される。例えば、毛細血管等
が発達した部分は明るく表示されるが、腫瘍などで石灰
化した部分は暗く表示される。この方法によれば、異な
る部位における分子の可動性の違いが明瞭である。
【０００３】
【発明が解決しようとする課題】しかしながら、上記従
来技術には、異なる部位間での拡散係数の相対的な比較
には適するが、定量的な評価には不適であるという問題
があった。また、上記従来技術にはスピン密度等、拡散
係数以外の量を同時に表示できない、という問題があっ
た。本発明は、上記欠点に鑑みてなされたもので、その
目的とするところは、従来の技術における上述の如き問
題を解消し、定量的な評価が可能であり、拡散係数以外
の量を同時に表示可能な拡散係数分布画像の表示方法を
提供することにある。
来技術には、異なる部位間での拡散係数の相対的な比較
には適するが、定量的な評価には不適であるという問題
があった。また、上記従来技術にはスピン密度等、拡散
係数以外の量を同時に表示できない、という問題があっ
た。本発明は、上記欠点に鑑みてなされたもので、その
目的とするところは、従来の技術における上述の如き問
題を解消し、定量的な評価が可能であり、拡散係数以外
の量を同時に表示可能な拡散係数分布画像の表示方法を
提供することにある。
【０００４】
【課題を解決するための手段】上記課題を解決するため
に、本発明では、拡散係数に画像の表示色を対応させ
る。また、拡散係数の分布画像上の任意の領域内の拡散
係数分布を１次元表示させる。この１次元表示した拡散
係数分布のヒストグラムの横軸上の任意の範囲を指定す
ることにより、拡散係数の分布画像上において、上記範
囲の拡散係数を有するボクセルのみ画像として表示させ
たる。また、２つの拡散係数分布画像に対して上記の１
次元表示を行い、両者の差を１次元表示させる。２つの
拡散分布画像を用いた上記のような表示を行う場合に
は、拡散係数の計測誤差を補正した後に、上記の表示を
行う。さらに、拡散係数以外の量を同時に表示するため
に、これらの量を画像の濃淡に対応させる。
に、本発明では、拡散係数に画像の表示色を対応させ
る。また、拡散係数の分布画像上の任意の領域内の拡散
係数分布を１次元表示させる。この１次元表示した拡散
係数分布のヒストグラムの横軸上の任意の範囲を指定す
ることにより、拡散係数の分布画像上において、上記範
囲の拡散係数を有するボクセルのみ画像として表示させ
たる。また、２つの拡散係数分布画像に対して上記の１
次元表示を行い、両者の差を１次元表示させる。２つの
拡散分布画像を用いた上記のような表示を行う場合に
は、拡散係数の計測誤差を補正した後に、上記の表示を
行う。さらに、拡散係数以外の量を同時に表示するため
に、これらの量を画像の濃淡に対応させる。
【０００５】すなわち、静磁場と傾斜磁場および高周波
磁場の発生手段と、検査対象からの核磁気共鳴信号を検
出する信号検出手段と、信号検出手段の検出信号の処理
を行う計算機と、検出信号から得られる画像を表示する
表示手段とを備えた核磁気共鳴検査装置における拡散係
数分布画像の表示方法において、検出信号から得られる
画像の表示色を検査対象から得られた拡散係数に対応さ
せることに特徴を有する。また、画像の表示色を検査対
象から得られた拡散係数に対応させ、他の量を画像の濃
淡に対応させることを特徴とする。
磁場の発生手段と、検査対象からの核磁気共鳴信号を検
出する信号検出手段と、信号検出手段の検出信号の処理
を行う計算機と、検出信号から得られる画像を表示する
表示手段とを備えた核磁気共鳴検査装置における拡散係
数分布画像の表示方法において、検出信号から得られる
画像の表示色を検査対象から得られた拡散係数に対応さ
せることに特徴を有する。また、画像の表示色を検査対
象から得られた拡散係数に対応させ、他の量を画像の濃
淡に対応させることを特徴とする。
【０００６】
【作用】拡散係数分布画像を表示する際に、拡散係数値
に画像の表示色を対応させることにより、拡散係数の分
布状態が視覚的に捕らえやすくなる。また、画像表示に
用いる色の波長範囲と表示する拡散係数の範囲を可変と
することにより、様々な計測部位に対応可能となる。ま
た上記の画像の任意の領域内の拡散係数分布を１次元表
示することにより、異なる部位間において拡散係数の定
量的な評価が可能である。また上記の１次元表示した拡
散係数分布のヒストグラムの横軸上の任意の範囲を指定
することにより、拡散係数の分布画像上において、この
範囲の拡散係数を有するボクセルのみ画像として表示す
るので、特定の範囲の分子可動性を有する部位を抽出す
ることができる。また、上記の１次元表示した拡散係数
分布のヒストグラムの横軸上の任意の範囲を２か所指定
することにより、拡散係数分布画像上の任意の領域にお
いて上記の２つの範囲の拡散係数を有するボクセル数の
比率を表示するので、異なる部位間、あるいは異なる拡
散係数分布画像間における拡散係数の定量的な評価が可
能である。また、２つの拡散係数分布画像に対して上記
の１次元表示を行い、両者の差を１次元表示するので、
異なる拡散係数分布画像間において定量的な評価が可能
である。また、２つの拡散分布画像を用いた上記のよう
な表示を行う場合には、拡散係数の計測誤差を補正した
後に上記表示を行うので、計測条件等の違いによる計測
誤差のために定量性を失うこと無く表示を行うことがで
きる。また、拡散係数以外の量を画像の濃淡に対応させ
ることができるので、拡散係数分布と同時に表示するこ
とが可能である。
に画像の表示色を対応させることにより、拡散係数の分
布状態が視覚的に捕らえやすくなる。また、画像表示に
用いる色の波長範囲と表示する拡散係数の範囲を可変と
することにより、様々な計測部位に対応可能となる。ま
た上記の画像の任意の領域内の拡散係数分布を１次元表
示することにより、異なる部位間において拡散係数の定
量的な評価が可能である。また上記の１次元表示した拡
散係数分布のヒストグラムの横軸上の任意の範囲を指定
することにより、拡散係数の分布画像上において、この
範囲の拡散係数を有するボクセルのみ画像として表示す
るので、特定の範囲の分子可動性を有する部位を抽出す
ることができる。また、上記の１次元表示した拡散係数
分布のヒストグラムの横軸上の任意の範囲を２か所指定
することにより、拡散係数分布画像上の任意の領域にお
いて上記の２つの範囲の拡散係数を有するボクセル数の
比率を表示するので、異なる部位間、あるいは異なる拡
散係数分布画像間における拡散係数の定量的な評価が可
能である。また、２つの拡散係数分布画像に対して上記
の１次元表示を行い、両者の差を１次元表示するので、
異なる拡散係数分布画像間において定量的な評価が可能
である。また、２つの拡散分布画像を用いた上記のよう
な表示を行う場合には、拡散係数の計測誤差を補正した
後に上記表示を行うので、計測条件等の違いによる計測
誤差のために定量性を失うこと無く表示を行うことがで
きる。また、拡散係数以外の量を画像の濃淡に対応させ
ることができるので、拡散係数分布と同時に表示するこ
とが可能である。
【０００７】
【実施例】以下、本発明の実施例を図面に基づいて説明
する。図２は核磁気共鳴を用いた検査装置の一例を示す
構成図である。図において、１は静磁場を発生するコイ
ル、２は傾斜磁場を発生するコイル、３は検査対象であ
り、この検査対象はコイル１、２内に配置される。シー
ケンサ４は傾斜磁場電源５、高周波発振器６に命令を送
り、傾斜磁場および高周波磁場を印加する。高周波磁場
は、高周波変調器７、高周波増幅器８を経て高周波送信
器９により、検査対象３に印加される。検査対象から発
生した信号は受信器１０によって受波され、増幅器１
１、位相検波器１２、ＡＤ変換器１３を通ってＣＰＵ１
４に送られ、ここで信号処理が行われた後、表示装置１
５に表示される。必要に応じて、記憶媒体１６に信号や
測定条件を記憶させることもできる。ここで、核磁気共
鳴を用いた拡散係数分布の計測法について簡単に述べて
おく。拡散とは分子のブラウン運動であり、拡散係数は
分子の動きやすさを表している。例えば４０℃の純水で
は０．００２５ｍｍ2／ｓ程度であるが、生体中の水で
は０．００２ｍｍ2／ｓ程度である。通常のイメージン
グでは、この拡散現象による核磁気共鳴信号の減少量は
非常にわずかであまり問題にならないが、リードアウト
傾斜磁場方向に強い傾斜磁場を印加すると信号減少量が
増加する。図３に示すようなスピンエコーシーケンスに
おいて、１８０°パルスの両側でこのような強い傾斜磁
場を印加すると、静止している部位では第１の傾斜磁場
と第２の傾斜磁場による影響を打消し合うが、動きの大
きい部位ほど両者の影響を受けるため、信号の減少量が
大きくなる。また、信号の減少量は傾斜磁場の印加量、
すなわち印加時間あるいは傾斜磁場強度に応じて変化す
る。ここで、ｂを傾斜磁場の印加量を表す傾斜磁場因
子、Ｄを拡散係数とすると、自己拡散による信号減少
は、ｅｘｐ（−ｂ・Ｄ）で表される。従って傾斜磁場因
子ｂの異なる複数の画像から、自己拡散係数Ｄを求める
ことができる。詳しくはマグネティック レゾナンス
イン メディスン,19,(1991年)第221頁から227頁(Magne
tic Resonance inMedicine 19,(1991)pp221-227)の記載
が参考になる。
する。図２は核磁気共鳴を用いた検査装置の一例を示す
構成図である。図において、１は静磁場を発生するコイ
ル、２は傾斜磁場を発生するコイル、３は検査対象であ
り、この検査対象はコイル１、２内に配置される。シー
ケンサ４は傾斜磁場電源５、高周波発振器６に命令を送
り、傾斜磁場および高周波磁場を印加する。高周波磁場
は、高周波変調器７、高周波増幅器８を経て高周波送信
器９により、検査対象３に印加される。検査対象から発
生した信号は受信器１０によって受波され、増幅器１
１、位相検波器１２、ＡＤ変換器１３を通ってＣＰＵ１
４に送られ、ここで信号処理が行われた後、表示装置１
５に表示される。必要に応じて、記憶媒体１６に信号や
測定条件を記憶させることもできる。ここで、核磁気共
鳴を用いた拡散係数分布の計測法について簡単に述べて
おく。拡散とは分子のブラウン運動であり、拡散係数は
分子の動きやすさを表している。例えば４０℃の純水で
は０．００２５ｍｍ2／ｓ程度であるが、生体中の水で
は０．００２ｍｍ2／ｓ程度である。通常のイメージン
グでは、この拡散現象による核磁気共鳴信号の減少量は
非常にわずかであまり問題にならないが、リードアウト
傾斜磁場方向に強い傾斜磁場を印加すると信号減少量が
増加する。図３に示すようなスピンエコーシーケンスに
おいて、１８０°パルスの両側でこのような強い傾斜磁
場を印加すると、静止している部位では第１の傾斜磁場
と第２の傾斜磁場による影響を打消し合うが、動きの大
きい部位ほど両者の影響を受けるため、信号の減少量が
大きくなる。また、信号の減少量は傾斜磁場の印加量、
すなわち印加時間あるいは傾斜磁場強度に応じて変化す
る。ここで、ｂを傾斜磁場の印加量を表す傾斜磁場因
子、Ｄを拡散係数とすると、自己拡散による信号減少
は、ｅｘｐ（−ｂ・Ｄ）で表される。従って傾斜磁場因
子ｂの異なる複数の画像から、自己拡散係数Ｄを求める
ことができる。詳しくはマグネティック レゾナンス
イン メディスン,19,(1991年)第221頁から227頁(Magne
tic Resonance inMedicine 19,(1991)pp221-227)の記載
が参考になる。
【０００８】網目状の微細構造を持つ毛細血管中の血流
等も、巨視的には疑似拡散と考えられる。上記のような
方法で生体の拡散係数を計測する場合、これらの影響と
真の拡散とを区別することはできない。両者の影響を含
めた疑似的な拡散係数を見かけの拡散係数（ＡＤＣ）と
いう。上記のような毛細血管を有する部位では、ＡＤＣ
は真の拡散係数Ｄの１０倍程度の値を示す。さて、本発
明の実施例として核磁気共鳴を用いたイメージング装置
により計測した生体の見かけの拡散係数分布を表示する
場合について述べる。例えば人体頭部のＡＤＣを上記の
ような方法で計測すると、骨や石灰質化した部位、灰白
質や白質等、純水の拡散係数に近いＡＤＣ値を示す部
位、毛細血管床や脳脊髄流等の部位ではＡＤＣ値に非常
に大きな幅があるので、これらを同一の画像上に表示す
る場合には、対数で表示すると便利である。計測したＡ
ＤＣ分布を表示するために、まずＡＤＣの対数を色に対
応させる。表示するＡＤＣ値の範囲を、Ｄ１〜Ｄ２と
し、表示色の波長の範囲を、Ｗ１〜Ｗ２とすると、ＡＤ
Ｃ値がＤｘであるボクセルの表示色の波長は、Ｗｘ＝Ｄ
ｘ（Ｗ２−Ｗ１）／（Ｄ２−Ｄ１）＋Ｗ１で与えられ
る。Ｄ１、Ｄ２、Ｗ１、Ｗ２の値としては、それぞれ例
えば０ｍｍ2／ｓ、０．０３ｍｍ2／ｓ、４０００Å、７
０００Å等とし、これらは計測データに応じて可変とす
る。上記の表示例では、拡散係数の高い部位ほど赤に近
い色で表示される。図１に示すように、表示色と拡散係
数値との対応関係を明確にするために、画像の端にカラ
ーバーを設け、対応する拡散係数値を表示する。拡散係
数値の他にスピン密度等、別の計測値を同時に表示した
い場合には、画像の濃淡値に上記の計測値を対応させ
る。
等も、巨視的には疑似拡散と考えられる。上記のような
方法で生体の拡散係数を計測する場合、これらの影響と
真の拡散とを区別することはできない。両者の影響を含
めた疑似的な拡散係数を見かけの拡散係数（ＡＤＣ）と
いう。上記のような毛細血管を有する部位では、ＡＤＣ
は真の拡散係数Ｄの１０倍程度の値を示す。さて、本発
明の実施例として核磁気共鳴を用いたイメージング装置
により計測した生体の見かけの拡散係数分布を表示する
場合について述べる。例えば人体頭部のＡＤＣを上記の
ような方法で計測すると、骨や石灰質化した部位、灰白
質や白質等、純水の拡散係数に近いＡＤＣ値を示す部
位、毛細血管床や脳脊髄流等の部位ではＡＤＣ値に非常
に大きな幅があるので、これらを同一の画像上に表示す
る場合には、対数で表示すると便利である。計測したＡ
ＤＣ分布を表示するために、まずＡＤＣの対数を色に対
応させる。表示するＡＤＣ値の範囲を、Ｄ１〜Ｄ２と
し、表示色の波長の範囲を、Ｗ１〜Ｗ２とすると、ＡＤ
Ｃ値がＤｘであるボクセルの表示色の波長は、Ｗｘ＝Ｄ
ｘ（Ｗ２−Ｗ１）／（Ｄ２−Ｄ１）＋Ｗ１で与えられ
る。Ｄ１、Ｄ２、Ｗ１、Ｗ２の値としては、それぞれ例
えば０ｍｍ2／ｓ、０．０３ｍｍ2／ｓ、４０００Å、７
０００Å等とし、これらは計測データに応じて可変とす
る。上記の表示例では、拡散係数の高い部位ほど赤に近
い色で表示される。図１に示すように、表示色と拡散係
数値との対応関係を明確にするために、画像の端にカラ
ーバーを設け、対応する拡散係数値を表示する。拡散係
数値の他にスピン密度等、別の計測値を同時に表示した
い場合には、画像の濃淡値に上記の計測値を対応させ
る。
【０００９】さて、同一画像上の異なる部位における拡
散係数や、異なる画像上の同一部位における拡散係数を
比較する場合、色の違いだけでは評価が難しい。そこ
で、図４に示すように、矩形で評価したい領域を画像上
で指定し、この領域に含まれるボクセルの拡散係数分布
のヒストグラムを１次元表示する。この場合、横軸は拡
散係数値、あるいは拡散係数値の対数とし、縦軸はボク
セル数とする。画像上で指定する領域の形状は、たとえ
ば矩形、円形、楕円形等とし、マウス等で指定する。毛
細血管が発達している部位等では高いＡＤＣ値を示すボ
クセルが多いので、拡散係数分布のヒストグラムを１次
元表示すると、他の部位に比べて分布状態がずれるので
容易に区別できる。また、図５に示すように上記の１次
元表示した拡散係数分布のヒストグラムにおいて、横軸
上の任意の範囲を指定することにより、拡散係数画像上
で特定範囲の拡散係数値を示す部位を選択的に表示す
る。たとえば、０．０１〜０．０３ｍｍ2／ｓの範囲に
ある拡散係数を選択すると、比較的可動性の高い部位だ
けが表示される。範囲指定方法は、２本のラインカーソ
ルを用いて１次元拡散係数分布上の任意の範囲を選択し
ても良いし、直接数値を入力しても良い。上記特定範囲
の拡散係数値を示す部位が、計測部位全体においてどこ
に位置するのかが不明な場合には、上記部位と同時に計
測部位全体の輪郭を表示するか、あるいは残りの部位を
薄く表示することにより、上記部位の位置関係を明確に
することができる。
散係数や、異なる画像上の同一部位における拡散係数を
比較する場合、色の違いだけでは評価が難しい。そこ
で、図４に示すように、矩形で評価したい領域を画像上
で指定し、この領域に含まれるボクセルの拡散係数分布
のヒストグラムを１次元表示する。この場合、横軸は拡
散係数値、あるいは拡散係数値の対数とし、縦軸はボク
セル数とする。画像上で指定する領域の形状は、たとえ
ば矩形、円形、楕円形等とし、マウス等で指定する。毛
細血管が発達している部位等では高いＡＤＣ値を示すボ
クセルが多いので、拡散係数分布のヒストグラムを１次
元表示すると、他の部位に比べて分布状態がずれるので
容易に区別できる。また、図５に示すように上記の１次
元表示した拡散係数分布のヒストグラムにおいて、横軸
上の任意の範囲を指定することにより、拡散係数画像上
で特定範囲の拡散係数値を示す部位を選択的に表示す
る。たとえば、０．０１〜０．０３ｍｍ2／ｓの範囲に
ある拡散係数を選択すると、比較的可動性の高い部位だ
けが表示される。範囲指定方法は、２本のラインカーソ
ルを用いて１次元拡散係数分布上の任意の範囲を選択し
ても良いし、直接数値を入力しても良い。上記特定範囲
の拡散係数値を示す部位が、計測部位全体においてどこ
に位置するのかが不明な場合には、上記部位と同時に計
測部位全体の輪郭を表示するか、あるいは残りの部位を
薄く表示することにより、上記部位の位置関係を明確に
することができる。
【００１０】また、拡散係数分布画像上で複数あるいは
１つの領域を選択し、これらの領域における高可動性部
位と低可動性部位との存在比率を知りたい場合には、図
６に示すように上記１次元した拡散係数分布のヒストグ
ラムにおいて、横軸上の任意の範囲を２か所指定するこ
とにより、これらの範囲の拡散係数値を示すボクセル数
の比率を数値あるいは画像の濃淡値として表示する。図
６は画像の濃淡値による表示例を示しているが、数値で
表示する場合には、この濃淡値に対応する値を各領域上
等に表示すれば良い。また、数値と画像の濃淡値を同時
に表示しても良い。拡散係数分布画像上で選択した領域
が１つである場合、または正確な値が知りたい場合には
数値表示が適しているが、複数の領域について上記の存
在比率を相対的に比較したい場合には、存在比率を濃淡
値に対応させて画像として表示する方がわかりやすい。
領域の指定方法は、マウス等で矩形、円形、楕円形等、
任意の形状の領域を指定しても良いし、拡散係数分布画
像全体を適当な大きさに等分しても良い。
１つの領域を選択し、これらの領域における高可動性部
位と低可動性部位との存在比率を知りたい場合には、図
６に示すように上記１次元した拡散係数分布のヒストグ
ラムにおいて、横軸上の任意の範囲を２か所指定するこ
とにより、これらの範囲の拡散係数値を示すボクセル数
の比率を数値あるいは画像の濃淡値として表示する。図
６は画像の濃淡値による表示例を示しているが、数値で
表示する場合には、この濃淡値に対応する値を各領域上
等に表示すれば良い。また、数値と画像の濃淡値を同時
に表示しても良い。拡散係数分布画像上で選択した領域
が１つである場合、または正確な値が知りたい場合には
数値表示が適しているが、複数の領域について上記の存
在比率を相対的に比較したい場合には、存在比率を濃淡
値に対応させて画像として表示する方がわかりやすい。
領域の指定方法は、マウス等で矩形、円形、楕円形等、
任意の形状の領域を指定しても良いし、拡散係数分布画
像全体を適当な大きさに等分しても良い。
【００１１】また、上記の操作は同一の拡散係数分布画
像に限らず、複数の画像に対して行っても良い。たとえ
ば、同一人物における薬物治療効果等を拡散係数分布を
用いて評価したい場合には、治療前後に計測した拡散係
数分布画像上で比較したい領域を選択し、上記のような
操作を行えば良い。同じ目的で、複数の拡散係数分布画
像上の同一領域における上記１次元表示した拡散係数分
布のヒストグラムを比較する場合、図７に示すように、
それぞれの画像に対応する上記１次元した拡散係数分布
のヒストグラムの他に、両者の差を表示すると、拡散係
数の経時変化がわかりやすい。領域の指定方法は、マウ
ス等で矩形、円形、楕円形等、任意の形状の領域を指定
しても良いし、拡散係数分布画像全体を適当な大きさに
等分しても良い。上記図６、図７に示す表示法等のよう
に複数の拡散係数分布画像上の拡散係数を比較する場合
には、計測条件が異なることにより拡散係数の計測値に
誤差が生じる場合がある。このような場合には、上記の
誤差を補正する必要がある。このような誤差を補正する
方法の例を述べる。特定の部位における拡散係数の計測
値を標準値として、上記複数の拡散係数分布画像間の規
格化を行う。この標準値としては、拡散係数分布画像の
計測時に同時に計測した標準試料の拡散係数を用いる
か、あるいは個体差による値の変動が小さい部位の拡散
係数値を用いる。上記標準試料しては、例えば水を用い
る。また、上記個体差による値の変動が小さい部位とし
ては、例えば骨を用いる。例えば２枚の拡散係数分布画
像Ｓ1（ｘ，ｙ，ｚ）、Ｓ2（ｘ，ｙ，ｚ）の規格化を行
う場合、上記標準値の比がα＝Ｓ1（ａ，ｂ，ｃ）／Ｓ2
（ａ，ｂ，ｃ）であれば、Ｓ2（ｘ，ｙ，ｚ）の代わり
に、Ｓ2’（ｘ，ｙ，ｚ）＝α・Ｓ2（ｘ，ｙ，ｚ）を用
いれば良い。ただし、（ａ，ｂ，ｃ）は拡散係数の標準
値として用いた部位の位置座標である。
像に限らず、複数の画像に対して行っても良い。たとえ
ば、同一人物における薬物治療効果等を拡散係数分布を
用いて評価したい場合には、治療前後に計測した拡散係
数分布画像上で比較したい領域を選択し、上記のような
操作を行えば良い。同じ目的で、複数の拡散係数分布画
像上の同一領域における上記１次元表示した拡散係数分
布のヒストグラムを比較する場合、図７に示すように、
それぞれの画像に対応する上記１次元した拡散係数分布
のヒストグラムの他に、両者の差を表示すると、拡散係
数の経時変化がわかりやすい。領域の指定方法は、マウ
ス等で矩形、円形、楕円形等、任意の形状の領域を指定
しても良いし、拡散係数分布画像全体を適当な大きさに
等分しても良い。上記図６、図７に示す表示法等のよう
に複数の拡散係数分布画像上の拡散係数を比較する場合
には、計測条件が異なることにより拡散係数の計測値に
誤差が生じる場合がある。このような場合には、上記の
誤差を補正する必要がある。このような誤差を補正する
方法の例を述べる。特定の部位における拡散係数の計測
値を標準値として、上記複数の拡散係数分布画像間の規
格化を行う。この標準値としては、拡散係数分布画像の
計測時に同時に計測した標準試料の拡散係数を用いる
か、あるいは個体差による値の変動が小さい部位の拡散
係数値を用いる。上記標準試料しては、例えば水を用い
る。また、上記個体差による値の変動が小さい部位とし
ては、例えば骨を用いる。例えば２枚の拡散係数分布画
像Ｓ1（ｘ，ｙ，ｚ）、Ｓ2（ｘ，ｙ，ｚ）の規格化を行
う場合、上記標準値の比がα＝Ｓ1（ａ，ｂ，ｃ）／Ｓ2
（ａ，ｂ，ｃ）であれば、Ｓ2（ｘ，ｙ，ｚ）の代わり
に、Ｓ2’（ｘ，ｙ，ｚ）＝α・Ｓ2（ｘ，ｙ，ｚ）を用
いれば良い。ただし、（ａ，ｂ，ｃ）は拡散係数の標準
値として用いた部位の位置座標である。
【００１２】
【発明の効果】以上、詳細に説明したとおり、本発明に
よれば、拡散係数分布画像を表示する際に、拡散係数に
画像の表示色を対応させ、また上記画像の任意の領域内
の拡散係数分布のヒストグラムを１次元表示し、また上
記１次元表示した拡散係数分布のヒストグラムの横軸上
の任意の範囲を指定することにより、拡散係数の分布画
像上において、上記範囲の拡散係数を有するボクセルの
み画像として表示するので、異なる部位間において拡散
係数の定量的な評価が可能である。また、上記１次元表
示した拡散係数分布のヒストグラムの横軸上の任意の範
囲を２か所指定することにより、拡散係数分布画像上の
任意の領域において上記２範囲の拡散係数を有するボク
セル数の比率を表示するので、異なる部位間、あるいは
異なる拡散係数分布画像間における拡散係数の定量的な
評価が可能である。また、２つの拡散係数分布画像に対
して上記１次元表示を行い、両者の差を１次元表示する
ので、異なる拡散係数分布画像間において、定量的な評
価が可能である。また、２つの拡散分布画像を用いた上
記のような表示を行う場合には、拡散係数の計測誤差を
補正した後に上記表示を行うので、計測条件等の違いに
よる計測誤差のために定量性を失うこと無く上記表示を
行うことができる。また、拡散係数以外の量を画像の濃
淡に対応させることができるので、拡散係数分布と同時
に表示することが可能である。
よれば、拡散係数分布画像を表示する際に、拡散係数に
画像の表示色を対応させ、また上記画像の任意の領域内
の拡散係数分布のヒストグラムを１次元表示し、また上
記１次元表示した拡散係数分布のヒストグラムの横軸上
の任意の範囲を指定することにより、拡散係数の分布画
像上において、上記範囲の拡散係数を有するボクセルの
み画像として表示するので、異なる部位間において拡散
係数の定量的な評価が可能である。また、上記１次元表
示した拡散係数分布のヒストグラムの横軸上の任意の範
囲を２か所指定することにより、拡散係数分布画像上の
任意の領域において上記２範囲の拡散係数を有するボク
セル数の比率を表示するので、異なる部位間、あるいは
異なる拡散係数分布画像間における拡散係数の定量的な
評価が可能である。また、２つの拡散係数分布画像に対
して上記１次元表示を行い、両者の差を１次元表示する
ので、異なる拡散係数分布画像間において、定量的な評
価が可能である。また、２つの拡散分布画像を用いた上
記のような表示を行う場合には、拡散係数の計測誤差を
補正した後に上記表示を行うので、計測条件等の違いに
よる計測誤差のために定量性を失うこと無く上記表示を
行うことができる。また、拡散係数以外の量を画像の濃
淡に対応させることができるので、拡散係数分布と同時
に表示することが可能である。
【図１】本発明の実施例を示す拡散係数分布画像の表示
方法。
方法。
【図２】本発明の実施例を示す核磁気共鳴を用いた検査
装置の構成図。
装置の構成図。
【図３】核磁気共鳴を用いた拡散係数分布の計測方法。
【図４】本発明の実施例を示す１次元拡散係数分布の表
示方法。
示方法。
【図５】本発明の実施例を示す拡散係数分布画像の他の
表示方法。
表示方法。
【図６】本発明の実施例を示す拡散係数分布画像の他の
表示方法。
表示方法。
【図７】本発明の実施例を示す１次元拡散係数分布の他
の表示方法。
の表示方法。
１…静磁場発生コイル、２…傾斜磁場発生コイル、３…
検査対象、４…シーケンサ、５…傾斜磁場電源、６…高
周波発振器、７…高周波変調器、８…高周波増幅器、９
…高周波送信器、１０…受信器、１１…増幅器、１２…
位相検波器、１３…ＡＤ変換器、１４…ＣＰＵ、１５…
表示装置、１６…記憶媒体。
検査対象、４…シーケンサ、５…傾斜磁場電源、６…高
周波発振器、７…高周波変調器、８…高周波増幅器、９
…高周波送信器、１０…受信器、１１…増幅器、１２…
位相検波器、１３…ＡＤ変換器、１４…ＣＰＵ、１５…
表示装置、１６…記憶媒体。
Claims (14)
- 【請求項１】静磁場と傾斜磁場および高周波磁場の発生
手段と、検査対象からの核磁気共鳴信号を検出する信号
検出手段と、この信号検出手段の検出信号の処理を行う
計算機と、前記検出信号から得られる画像を表示する表
示手段とを備えた核磁気共鳴検査装置における拡散係数
分布画像の表示方法において、前記画像の表示色を前記
検査対象から得られた拡散係数に対応させたことを特徴
とする、拡散係数分布画像の表示方法。 - 【請求項２】静磁場と傾斜磁場および高周波磁場の発生
手段と、検査対象からの核磁気共鳴信号を検出する信号
検出手段と、この信号検出手段の検出信号の処理を行う
計算機と、前記検出信号から得られる画像を表示する表
示手段とを備えた核磁気共鳴検査装置における拡散係数
分布画像の表示方法において、前記画像の表示色を前記
検査対象から得られた拡散係数に対応させ、他の量を画
像の濃淡に対応させることを特徴とする、拡散係数分布
画像の表示方法。 - 【請求項３】請求項１に記載の拡散係数分布画像の表示
方法は、前記拡散係数の対数と前記表示色の波長が１対
１に対応することを特徴とする、拡散係数分布画像の表
示方法。 - 【請求項４】請求項１に記載の拡散係数分布画像の表示
方法は、カラー表示の前記表示色と前記拡散係数の数値
とを対応づける手段を表示画像上に有することを特徴と
する、拡散係数分布画像の表示方法。 - 【請求項５】請求項３に記載の拡散係数分布画像の表示
方法は、前記拡散係数の対数と前記表示色の波長との対
応関係が１次式で記述できることを特徴とする、拡散係
数分布画像の表示方法。 - 【請求項６】請求項５に記載の拡散係数分布画像の表示
方法は、上記１次式の傾きおよび切片を任意の値に設定
することができることを特徴とする、拡散係数分布画像
の表示方法。 - 【請求項７】請求項１に記載の拡散係数分布画像の表示
方法は、前記拡散係数の分布画像上の任意の領域を指定
することにより、前記領域内の拡散係数分布のヒストグ
ラムを１次元表示することを特徴とする拡散係数分布画
像の表示方法。 - 【請求項８】請求項７に記載の拡散係数分布画像の表示
方法は、前記１次元表示の横軸を、前記拡散係数値、あ
るいは拡散係数の対数とし、縦軸をボクセル数とするこ
とを特徴とする、拡散係数分布画像の表示方法。 - 【請求項９】請求項７に記載の拡散係数分布画像の表示
方法は、前記１次元表示した拡散係数分布のヒストグラ
ム上の横軸上の任意の範囲を指定することにより、前記
拡散係数の分布画像上において前記範囲の拡散係数を有
するボクセルのみ表示するか、あるいはそれと共に残り
の部分については輪郭のみあるいは薄く表示することを
特徴とする、拡散係数分布画像の表示方法。 - 【請求項１０】請求項７に記載の拡散係数分布画像の表
示方法は、前記１次元表示した拡散係数分布のヒストグ
ラムの横軸上の任意の２つの範囲を指定することによ
り、１つあるいは複数の前記拡散係数の分布画像上にお
いて指定した１つあるいは複数の領域における、前記２
つの範囲の拡散係数を有するボクセル数の比率を数値あ
るいは画像の濃淡で表示することを特徴とする、拡散係
数分布画像の表示方法。 - 【請求項１１】請求項７に記載の拡散係数分布画像の表
示方法は、前記拡散係数分布の画像上において任意の領
域を指定することにより、前記領域内の拡散係数分布の
ヒストグラムを１次元表示し、他の拡散係数の分布画像
上において、前記領域と同一の領域内の拡散係数分布の
ヒストグラムを１次元表示し、前記２つの１次元表示し
た拡散係数分布のヒストグラムの差を１次元表示するこ
とを特徴とする拡散係数分布画像の表示方法。 - 【請求項１２】請求項１０あるいは請求項１１に記載の
拡散係数分布画像の表示方法は、前記拡散係数の分布画
像に対して、計測条件が異なることに由来する拡散係数
の誤差補正を行った後、前記表示を行うことを特徴とす
る、拡散係数分布画像の表示方法。 - 【請求項１３】請求項１２に記載の拡散係数分布画像の
表示方法は、標準試料を用いて拡散係数値を規格化する
ことにより、拡散係数の誤差補正を行うことを特徴とす
る、拡散係数分布画像の表示方法。 - 【請求項１４】請求項１２に記載の拡散係数分布画像の
表示方法は、個体により値が変化しないと考えられる部
位における拡散係数値を標準値として、該標準値を用い
て拡散係数値を規格化することにより、拡散係数の誤差
補正を行うことを特徴とする、拡散係数分布画像の表示
方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4210070A JPH0654830A (ja) | 1992-08-06 | 1992-08-06 | 拡散係数分布画像の表示方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4210070A JPH0654830A (ja) | 1992-08-06 | 1992-08-06 | 拡散係数分布画像の表示方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0654830A true JPH0654830A (ja) | 1994-03-01 |
Family
ID=16583326
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP4210070A Pending JPH0654830A (ja) | 1992-08-06 | 1992-08-06 | 拡散係数分布画像の表示方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0654830A (ja) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007203106A (ja) * | 2007-05-15 | 2007-08-16 | Toshiba Corp | Ｍｒｉ装置 |
| US7271838B2 (en) | 2002-05-08 | 2007-09-18 | Olympus Corporation | Image pickup apparatus with brightness distribution chart display capability |
| JP2008012172A (ja) * | 2006-07-07 | 2008-01-24 | Toshiba Corp | 磁気共鳴イメージング装置および画像処理装置 |
| JP2008059274A (ja) * | 2006-08-31 | 2008-03-13 | Toyota Motor Corp | 振動解析システム |
| KR100938863B1 (ko) * | 2008-05-28 | 2010-01-27 | 재단법인 아산사회복지재단 | 자동화된 목적 장기의 체절 및 암 조직의 시각화를 위한확산 강조 영상의 처리장치 및 그 방법 |
| JP2011062508A (ja) * | 2009-08-20 | 2011-03-31 | Toshiba Corp | 磁気共鳴イメージング装置および高周波コイルユニット |
| JP2012130759A (ja) * | 2006-01-24 | 2012-07-12 | Schlumberger Holdings Ltd | 骨構造をキャラクタライズするための拡散型磁気共鳴法の手順 |
| WO2015133363A1 (ja) * | 2014-03-05 | 2015-09-11 | 国立大学法人京都大学 | 画像処理装置、磁気共鳴イメージング装置、画像処理方法、および記録媒体 |
| JP2016000102A (ja) * | 2014-06-11 | 2016-01-07 | ザイオソフト株式会社 | 画像処理装置、画像処理方法、及び画像処理プログラム |
| US9505159B2 (en) | 2014-04-25 | 2016-11-29 | Medtronic Ablation Frontiers Llc | Methods of dimensionally stabilizing a lumen of a multi-lumen device during manufacture |
-
1992
- 1992-08-06 JP JP4210070A patent/JPH0654830A/ja active Pending
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7271838B2 (en) | 2002-05-08 | 2007-09-18 | Olympus Corporation | Image pickup apparatus with brightness distribution chart display capability |
| US7649564B2 (en) | 2002-05-08 | 2010-01-19 | Olympus Corporation | Image pickup apparatus with brightness distribution chart display capability |
| JP2012130759A (ja) * | 2006-01-24 | 2012-07-12 | Schlumberger Holdings Ltd | 骨構造をキャラクタライズするための拡散型磁気共鳴法の手順 |
| US7701210B2 (en) | 2006-07-07 | 2010-04-20 | Kabushiki Kaisha Toshiba | MRI apparatus, method and process generating diffusion weighted images using apparent diffusion coefficients and threshold values |
| JP2008012172A (ja) * | 2006-07-07 | 2008-01-24 | Toshiba Corp | 磁気共鳴イメージング装置および画像処理装置 |
| JP2008059274A (ja) * | 2006-08-31 | 2008-03-13 | Toyota Motor Corp | 振動解析システム |
| JP2007203106A (ja) * | 2007-05-15 | 2007-08-16 | Toshiba Corp | Ｍｒｉ装置 |
| KR100938863B1 (ko) * | 2008-05-28 | 2010-01-27 | 재단법인 아산사회복지재단 | 자동화된 목적 장기의 체절 및 암 조직의 시각화를 위한확산 강조 영상의 처리장치 및 그 방법 |
| JP2011062508A (ja) * | 2009-08-20 | 2011-03-31 | Toshiba Corp | 磁気共鳴イメージング装置および高周波コイルユニット |
| WO2015133363A1 (ja) * | 2014-03-05 | 2015-09-11 | 国立大学法人京都大学 | 画像処理装置、磁気共鳴イメージング装置、画像処理方法、および記録媒体 |
| JP2015167610A (ja) * | 2014-03-05 | 2015-09-28 | 国立大学法人京都大学 | 画像処理装置、磁気共鳴イメージング装置、画像処理方法、およびプログラム |
| US10401457B2 (en) | 2014-03-05 | 2019-09-03 | Kyoto University | Image processing apparatus, magnetic resonance imaging apparatus, image processing method, and storage medium |
| US9505159B2 (en) | 2014-04-25 | 2016-11-29 | Medtronic Ablation Frontiers Llc | Methods of dimensionally stabilizing a lumen of a multi-lumen device during manufacture |
| JP2016000102A (ja) * | 2014-06-11 | 2016-01-07 | ザイオソフト株式会社 | 画像処理装置、画像処理方法、及び画像処理プログラム |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5311131A (en) | Magnetic resonance imaging using pattern recognition | |
| US6310477B1 (en) | MR imaging of lesions and detection of malignant tumors | |
| KR100828220B1 (ko) | 통계적 영상들을 이용한 단층촬영 측정의 슬라이스 위치 계획을 위한 방법 | |
| RU2216751C2 (ru) | Способ и устройство для формирования изображений магнитного резонанса | |
| US5204625A (en) | Segmentation of stationary and vascular surfaces in magnetic resonance imaging | |
| US8874189B2 (en) | Method of visualizing MR images | |
| US20030006770A1 (en) | Quantitation and standardization of magnetic resonance measurements | |
| JP2004527301A (ja) | 連続的テーブル移動を使用して大視野からｍｒｉデータを取得する方法 | |
| JP2006527636A (ja) | 時間分解された磁気共鳴血管造影法のためのバックグラウンド抑制法 | |
| JPH0871058A (ja) | 磁気共鳴映像装置 | |
| US7145334B2 (en) | Method and magnetic resonance tomography apparatus for producing phase-coded flow images | |
| US10451700B2 (en) | System and method for reducing partial voluming artifacts in quantitative myocardial tissue characterization | |
| JPH0654830A (ja) | 拡散係数分布画像の表示方法 | |
| US8995738B2 (en) | System and method for magnetic resonance imaging parametric mapping using confidence maps | |
| JPH1066685A (ja) | 磁気共鳴画像装置とその方法 | |
| JP2002511330A (ja) | 磁気共鳴エコー平面イメージングにおける磁場の不均質性から生じる信号欠如の回復 | |
| US7190163B2 (en) | Method for producing multiple MR images with different contrast from a single image acquisition | |
| US7317940B2 (en) | MRT apparatus, method and computer program product for determining the speed of a flowing medium | |
| US20040162482A1 (en) | MRT apparatus, method and computer program product for speed-resolved flow measurement | |
| JPH06327649A (ja) | 核磁気共鳴イメージング装置における画像再構成方法 | |
| US10908247B2 (en) | System and method for texture analysis in magnetic resonance fingerprinting (MRF) | |
| JPH07108288B2 (ja) | Ｎｍｒイメ−ジング方法 | |
| US6342785B1 (en) | Method for quantification and representation of in vivo fluid flow using magnetic resonance | |
| JPS61106140A (ja) | 高分解能核磁気共鳴イメ−ジング方式 | |
| JP3104985B2 (ja) | 磁気共鳴診断装置 |