JPH1189839A - 灌流パラメータ測定方法、灌流画像生成方法および超音波診断装置 - Google Patents
灌流パラメータ測定方法、灌流画像生成方法および超音波診断装置Info
- Publication number
- JPH1189839A JPH1189839A JP25959897A JP25959897A JPH1189839A JP H1189839 A JPH1189839 A JP H1189839A JP 25959897 A JP25959897 A JP 25959897A JP 25959897 A JP25959897 A JP 25959897A JP H1189839 A JPH1189839 A JP H1189839A
- Authority
- JP
- Japan
- Prior art keywords
- perfusion
- image
- signal intensity
- pixel
- time
- 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
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
(57)【要約】 （修正有）
【課題】 灌流の活発さの度合いを表す灌流パラメータ
を定量的に測定する。 【解決手段】 造影剤としての気泡が存在する状態で信
号強度ｇk -を測定することによって気泡を破壊する。次
に、時間Ｔk を経過した後、再び、信号強度ｇk +を測定
する。もし、灌流が活発であるなら、気泡を含んだ血液
の流入によって速やかに気泡が補填されるから、比較的
短い時間Ｔk で信号強度ｇk +は信号強度ｇk-まで回復す
る。一方、灌流が不活発なら、比較的長い時間Ｔk をか
けないと信号強度ｇk +は信号強度ｇk -まで回復しない。
換言すれば、信号強度ｇk +が信号強度ｇk -の例えば９０
％まで回復してくる最短時間Ｔkminに基づいて灌流の活
発さの度合いを表す灌流パラメータの値を決める。これ
を画像の全画素について行い、灌流パラメータの値に応
じて色分けし、表示する。
を定量的に測定する。 【解決手段】 造影剤としての気泡が存在する状態で信
号強度ｇk -を測定することによって気泡を破壊する。次
に、時間Ｔk を経過した後、再び、信号強度ｇk +を測定
する。もし、灌流が活発であるなら、気泡を含んだ血液
の流入によって速やかに気泡が補填されるから、比較的
短い時間Ｔk で信号強度ｇk +は信号強度ｇk-まで回復す
る。一方、灌流が不活発なら、比較的長い時間Ｔk をか
けないと信号強度ｇk +は信号強度ｇk -まで回復しない。
換言すれば、信号強度ｇk +が信号強度ｇk -の例えば９０
％まで回復してくる最短時間Ｔkminに基づいて灌流の活
発さの度合いを表す灌流パラメータの値を決める。これ
を画像の全画素について行い、灌流パラメータの値に応
じて色分けし、表示する。
Description
【０００１】
【発明の属する技術分野】本発明は、灌流パラメータ測
定方法、灌流画像生成方法および超音波診断装置に関
し、さらに詳しくは、灌流（臓器の血液の入れ換わり）
の活発さの度合いを表す灌流パラメータを測定する灌流
パラメータ測定方法、灌流の活発さの度合いの２次元分
布を表す灌流画像を生成する灌流画像生成方法、およ
び、前記灌流画像生成方法を実施しうる超音波診断装置
に関する。
定方法、灌流画像生成方法および超音波診断装置に関
し、さらに詳しくは、灌流（臓器の血液の入れ換わり）
の活発さの度合いを表す灌流パラメータを測定する灌流
パラメータ測定方法、灌流の活発さの度合いの２次元分
布を表す灌流画像を生成する灌流画像生成方法、およ
び、前記灌流画像生成方法を実施しうる超音波診断装置
に関する。
【０００２】
【従来の技術】造影剤として気泡を被検体内に注入し、
超音波エコーの信号強度画像を撮像する技術が知られて
いる。
超音波エコーの信号強度画像を撮像する技術が知られて
いる。
【０００３】
【発明が解決しようとする課題】しかし、従来、造影剤
としての気泡を利用して、灌流の活発さの度合いを測定
する技術は知られていない。また、従来、造影剤として
の気泡を利用して、灌流の活発さの度合いの２次元分布
を表す灌流画像を生成する技術は知られていない。そこ
で、本発明の目的は、灌流の活発さの度合いを表す灌流
パラメータを測定する灌流パラメータ測定方法、灌流の
活発さの度合いの２次元分布を表す灌流画像を生成する
灌流画像生成方法、および、前記灌流画像生成方法を実
施しうる超音波診断装置を提供することにある。
としての気泡を利用して、灌流の活発さの度合いを測定
する技術は知られていない。また、従来、造影剤として
の気泡を利用して、灌流の活発さの度合いの２次元分布
を表す灌流画像を生成する技術は知られていない。そこ
で、本発明の目的は、灌流の活発さの度合いを表す灌流
パラメータを測定する灌流パラメータ測定方法、灌流の
活発さの度合いの２次元分布を表す灌流画像を生成する
灌流画像生成方法、および、前記灌流画像生成方法を実
施しうる超音波診断装置を提供することにある。
【０００４】
【課題を解決するための手段】第１の観点では、本発明
は、被検体内の測定点に造影剤としての気泡が十分存在
する状態で当該測定点からの超音波エコーの信号強度ｇ
k -を測定し、その測定のために送信した超音波ビームに
より前記気泡を破壊し、時間Ｔk 後に前記測定点からの
超音波エコーの信号強度ｇk +を再び測定し、比較係数ｗ
を０＜ｗ＜１とするとき、 ｇk +≧ｗ・ｇk - となるか否かを判定することを、前記時間Ｔk を変えな
がら繰り返して、前記不等式を満たす最短時間Ｔkminを
求め、その最短時間Ｔkminに基づいて灌流の活発さの度
合いを表すパラメータを決定することを特徴とする灌流
パラメータ測定方法を提供する。上記第１の観点による
灌流パラメータ測定方法では、気泡が存在する状態で信
号強度ｇk -を測定することによって気泡を破壊する。次
に、時間Ｔk を経過した後、再び、信号強度ｇk +を測定
する。もし、灌流が活発であるなら、気泡を含んだ血液
の流入によって速やかに気泡が補填されるから、比較的
短い時間Ｔk で信号強度ｇk +は信号強度ｇk -まで回復す
る。一方、灌流が不活発なら、気泡の補填が遅いから、
比較的長い時間Ｔk をかけないと信号強度ｇk +は信号強
度ｇk -まで回復しない。換言すれば、信号強度ｇk +が信
号強度ｇk -の（ｗ×１００）％まで回復してくる最短時
間Ｔkminは、臓器の血液の入れ換わりが活発なら短くな
り、臓器の血液の入れ換わりが不活発なら長くなる。よ
って、この最短時間Ｔkminに基づいて灌流の活発さの度
合いを表すパラメータを決めることが出来る。
は、被検体内の測定点に造影剤としての気泡が十分存在
する状態で当該測定点からの超音波エコーの信号強度ｇ
k -を測定し、その測定のために送信した超音波ビームに
より前記気泡を破壊し、時間Ｔk 後に前記測定点からの
超音波エコーの信号強度ｇk +を再び測定し、比較係数ｗ
を０＜ｗ＜１とするとき、 ｇk +≧ｗ・ｇk - となるか否かを判定することを、前記時間Ｔk を変えな
がら繰り返して、前記不等式を満たす最短時間Ｔkminを
求め、その最短時間Ｔkminに基づいて灌流の活発さの度
合いを表すパラメータを決定することを特徴とする灌流
パラメータ測定方法を提供する。上記第１の観点による
灌流パラメータ測定方法では、気泡が存在する状態で信
号強度ｇk -を測定することによって気泡を破壊する。次
に、時間Ｔk を経過した後、再び、信号強度ｇk +を測定
する。もし、灌流が活発であるなら、気泡を含んだ血液
の流入によって速やかに気泡が補填されるから、比較的
短い時間Ｔk で信号強度ｇk +は信号強度ｇk -まで回復す
る。一方、灌流が不活発なら、気泡の補填が遅いから、
比較的長い時間Ｔk をかけないと信号強度ｇk +は信号強
度ｇk -まで回復しない。換言すれば、信号強度ｇk +が信
号強度ｇk -の（ｗ×１００）％まで回復してくる最短時
間Ｔkminは、臓器の血液の入れ換わりが活発なら短くな
り、臓器の血液の入れ換わりが不活発なら長くなる。よ
って、この最短時間Ｔkminに基づいて灌流の活発さの度
合いを表すパラメータを決めることが出来る。
【０００５】第２の観点では、本発明は、被検体内に造
影剤としての気泡が十分存在する状態で超音波エコーの
信号強度画像ｉk -を撮像し、その撮像のために送信した
超音波ビームにより前記気泡を破壊し、時間Ｔk 後に超
音波エコーの信号強度画像ｉk+を再び撮像し、前記信号
強度画像ｉk -の各画素の信号強度をｇk -とし、前記信号
強度画像ｉk +の各画素の信号強度をｇk +とし、比較係数
ｗを０＜ｗ＜１とするとき、 ｇk +≧ｗ・ｇk - となるか否かを判定することを、前記時間Ｔk を変えな
がら繰り返して、各画素の信号強度が前記不等式を満た
す最短時間Ｔkminを求め、その最短時間Ｔkminに基づい
て灌流画像の各画素の画素値を決定することを特徴とす
る灌流画像生成方法を提供する。上記第２の観点による
灌流画像生成方法では、信号強度画像の各画素毎に上記
第１の観点による灌流パラメータ測定方法を実施するこ
とになるから、灌流の活発さの度合いの２次元分布を表
す灌流画像を生成できる。
影剤としての気泡が十分存在する状態で超音波エコーの
信号強度画像ｉk -を撮像し、その撮像のために送信した
超音波ビームにより前記気泡を破壊し、時間Ｔk 後に超
音波エコーの信号強度画像ｉk+を再び撮像し、前記信号
強度画像ｉk -の各画素の信号強度をｇk -とし、前記信号
強度画像ｉk +の各画素の信号強度をｇk +とし、比較係数
ｗを０＜ｗ＜１とするとき、 ｇk +≧ｗ・ｇk - となるか否かを判定することを、前記時間Ｔk を変えな
がら繰り返して、各画素の信号強度が前記不等式を満た
す最短時間Ｔkminを求め、その最短時間Ｔkminに基づい
て灌流画像の各画素の画素値を決定することを特徴とす
る灌流画像生成方法を提供する。上記第２の観点による
灌流画像生成方法では、信号強度画像の各画素毎に上記
第１の観点による灌流パラメータ測定方法を実施するこ
とになるから、灌流の活発さの度合いの２次元分布を表
す灌流画像を生成できる。
【０００６】第３の観点では、本発明は、上記構成の灌
流画像生成方法において、実質的に灌流がされない程度
に前記時間Ｔk を短くして信号強度画像ｉk +’を撮像
し、その信号強度画像ｉk +’の各画素の画素値と前記灌
流画像の画素値とに基づいて灌流画像の各画素の表示値
を決定することを特徴とする灌流画像生成方法を提供す
る。実質的に灌流がされない程度に時間Ｔk を短くする
と、造影剤としての気泡が破壊されたままであるから、
撮像された信号強度画像ｉk +’の画素値は、臓器の構造
を表した値となる。一方、灌流画像の画素値は、その部
分の灌流の活発さを表している。そこで、両者を基に当
該画素の表示値を決定すれば、臓器の構造と灌流の活発
さの度合いとを同時に表した画像表示が得られ、患部を
的確に視認できるようになる。
流画像生成方法において、実質的に灌流がされない程度
に前記時間Ｔk を短くして信号強度画像ｉk +’を撮像
し、その信号強度画像ｉk +’の各画素の画素値と前記灌
流画像の画素値とに基づいて灌流画像の各画素の表示値
を決定することを特徴とする灌流画像生成方法を提供す
る。実質的に灌流がされない程度に時間Ｔk を短くする
と、造影剤としての気泡が破壊されたままであるから、
撮像された信号強度画像ｉk +’の画素値は、臓器の構造
を表した値となる。一方、灌流画像の画素値は、その部
分の灌流の活発さを表している。そこで、両者を基に当
該画素の表示値を決定すれば、臓器の構造と灌流の活発
さの度合いとを同時に表した画像表示が得られ、患部を
的確に視認できるようになる。
【０００７】第４の観点では、本発明は、上記構成の灌
流画像生成方法において、位置の異なる複数の灌流画像
を基に灌流３次元画像を生成することを特徴とする灌流
画像生成方法を提供する。灌流３次元画像を利用すれ
ば、灌流の活発さの度合いの３次元分布を認識できるよ
うになる。
流画像生成方法において、位置の異なる複数の灌流画像
を基に灌流３次元画像を生成することを特徴とする灌流
画像生成方法を提供する。灌流３次元画像を利用すれ
ば、灌流の活発さの度合いの３次元分布を認識できるよ
うになる。
【０００８】第５の観点では、本発明は、被検体内に造
影剤として存在する気泡を破壊しうる超音波ビームによ
り超音波エコーの第１の信号強度画像ｉk -を撮像し続い
て時間Ｔk 後に第２の信号強度画像ｉk +を撮像する撮像
手段と、前記第１の信号強度画像ｉk -の各画素の信号強
度をｇk -とし，前記第２の信号強度画像ｉk +の各画素の
信号強度をｇk +とし，比較係数ｗを０＜ｗ＜１とすると
き，ｇk +≧ｗ・ｇk -となるか否かを判定する判定手段
と、前記撮像手段と前記判定手段とにより前記時間Ｔk
を変えながら前記撮像と前記判定とを繰り返して，各画
素の信号強度が前記不等式を満たす最短時間Ｔkminを求
め，その最短時間Ｔkminに基づいて灌流画像の各画素の
画素値を決定する画素値決定手段とを具備したことを特
徴とする超音波診断装置を提供する。上記第５の観点に
よる超音波診断装置では、上記第２の観点による灌流画
像生成方法を好適に実施することが出来る。
影剤として存在する気泡を破壊しうる超音波ビームによ
り超音波エコーの第１の信号強度画像ｉk -を撮像し続い
て時間Ｔk 後に第２の信号強度画像ｉk +を撮像する撮像
手段と、前記第１の信号強度画像ｉk -の各画素の信号強
度をｇk -とし，前記第２の信号強度画像ｉk +の各画素の
信号強度をｇk +とし，比較係数ｗを０＜ｗ＜１とすると
き，ｇk +≧ｗ・ｇk -となるか否かを判定する判定手段
と、前記撮像手段と前記判定手段とにより前記時間Ｔk
を変えながら前記撮像と前記判定とを繰り返して，各画
素の信号強度が前記不等式を満たす最短時間Ｔkminを求
め，その最短時間Ｔkminに基づいて灌流画像の各画素の
画素値を決定する画素値決定手段とを具備したことを特
徴とする超音波診断装置を提供する。上記第５の観点に
よる超音波診断装置では、上記第２の観点による灌流画
像生成方法を好適に実施することが出来る。
【０００９】
【発明の実施の形態】以下、図に示す発明の実施の形態
により本発明をさらに詳細に説明する。なお、これによ
り本発明が限定されるものではない。
により本発明をさらに詳細に説明する。なお、これによ
り本発明が限定されるものではない。
【００１０】図１は、本発明の一実施形態にかかる超音
波診断装置の構成図である。この超音波診断装置１００
は、超音波パルスを送信し超音波エコーを受信する超音
波探触子１と、走査平面を電子走査して音線信号を取得
する送受信制御部２と、前記超音波エコーの信号強度に
基づく画像データを生成するＢモード処理部３と、前記
超音波エコーのドプラ成分の位相に基づく画像データを
生成するＣＦ（Ｃolor Ｆlow）モード処理部４と、前記
超音波エコーのドプラ成分のパワーに基づく画像データ
を生成するＰＤ（Ｐower Doppler）モード処理部５と、
全体的な動作の制御を行うと共に前記画像データにより
表示画像を生成する中央処理部６と、前記表示画像を表
示するＣＲＴ７と、操作者が指示を入力する入力装置８
とを有している。また、前記中央処理部６は、灌流の活
発さの度合いの２次元分布を表す灌流画像を生成する灌
流画像生成部９と、灌流画像の画素値を表示値に変換し
て表示する灌流画像表示部１０と、灌流３次元画像を生
成し表示する灌流３Ｄ画像表示部１１とを含んでいる。
波診断装置の構成図である。この超音波診断装置１００
は、超音波パルスを送信し超音波エコーを受信する超音
波探触子１と、走査平面を電子走査して音線信号を取得
する送受信制御部２と、前記超音波エコーの信号強度に
基づく画像データを生成するＢモード処理部３と、前記
超音波エコーのドプラ成分の位相に基づく画像データを
生成するＣＦ（Ｃolor Ｆlow）モード処理部４と、前記
超音波エコーのドプラ成分のパワーに基づく画像データ
を生成するＰＤ（Ｐower Doppler）モード処理部５と、
全体的な動作の制御を行うと共に前記画像データにより
表示画像を生成する中央処理部６と、前記表示画像を表
示するＣＲＴ７と、操作者が指示を入力する入力装置８
とを有している。また、前記中央処理部６は、灌流の活
発さの度合いの２次元分布を表す灌流画像を生成する灌
流画像生成部９と、灌流画像の画素値を表示値に変換し
て表示する灌流画像表示部１０と、灌流３次元画像を生
成し表示する灌流３Ｄ画像表示部１１とを含んでいる。
【００１１】図２は、前記灌流画像生成部９により実行
される灌流画像生成処理を示すフロー図である。なお、
灌流画像生成処理を実行させる前提として、造影剤とし
ての気泡が被検体内に注入されているものとする。ステ
ップＳ１では、灌流画像の全画素値を“０”に初期化す
る。ステップＳ２では、撮像対象の臓器に適したＴｋテ
ーブルを選択する。図３の（ａ）〜（ｃ）に３種類のＴ
ｋテーブルを例示する。これらＴｋテーブルは、ｋ＝
０，１，２，…，ｅ（ｋの定義の終値）に対する時間Ｔ
k を定義した表である。時間Ｔ0 は、当該Ｔｋテーブル
に対応する臓器が正常なときに灌流が十分に行われる時
間である。また、ｋが大きくなるほど、時間Ｔk が長く
なっている。例えば、図３の（ａ）は、肝臓に適したＴ
ｋテーブルであり、Ｔ0 ＝１０秒，Ｔ1＝０．５秒，Ｔ2
＝１秒，Ｔ3 ＝２秒，Ｔ4 ＝４秒，…というように時
間間隔が等比的に長くなるように定義されている。ま
た、図３の（ｂ）のＴｋテーブルでは、時間間隔が等差
的に長くなるように定義されている。また、図３の
（ｃ）のＴｋテーブルでは、時間間隔が等差的に長くな
るが、その差が最初は小さく後は大きくなるように２段
階に定義されている。一般的には、灌流が速い臓器に対
しては時間間隔が短いＴｋテーブルを選択し、灌流が遅
い臓器に対しては時間間隔が長いＴｋテーブルを選択す
る。ステップＳ３では、シーケンスカウンタｋの値を
“１”に初期化する。
される灌流画像生成処理を示すフロー図である。なお、
灌流画像生成処理を実行させる前提として、造影剤とし
ての気泡が被検体内に注入されているものとする。ステ
ップＳ１では、灌流画像の全画素値を“０”に初期化す
る。ステップＳ２では、撮像対象の臓器に適したＴｋテ
ーブルを選択する。図３の（ａ）〜（ｃ）に３種類のＴ
ｋテーブルを例示する。これらＴｋテーブルは、ｋ＝
０，１，２，…，ｅ（ｋの定義の終値）に対する時間Ｔ
k を定義した表である。時間Ｔ0 は、当該Ｔｋテーブル
に対応する臓器が正常なときに灌流が十分に行われる時
間である。また、ｋが大きくなるほど、時間Ｔk が長く
なっている。例えば、図３の（ａ）は、肝臓に適したＴ
ｋテーブルであり、Ｔ0 ＝１０秒，Ｔ1＝０．５秒，Ｔ2
＝１秒，Ｔ3 ＝２秒，Ｔ4 ＝４秒，…というように時
間間隔が等比的に長くなるように定義されている。ま
た、図３の（ｂ）のＴｋテーブルでは、時間間隔が等差
的に長くなるように定義されている。また、図３の
（ｃ）のＴｋテーブルでは、時間間隔が等差的に長くな
るが、その差が最初は小さく後は大きくなるように２段
階に定義されている。一般的には、灌流が速い臓器に対
しては時間間隔が短いＴｋテーブルを選択し、灌流が遅
い臓器に対しては時間間隔が長いＴｋテーブルを選択す
る。ステップＳ３では、シーケンスカウンタｋの値を
“１”に初期化する。
【００１２】ステップＳ４では、時間Ｔ0 が経過するの
を待ってから、ステップＳ５へ進む。これは、灌流を十
分に行わせて、造影剤としての気泡を臓器に十分補填す
るためである。ステップＳ５では、第ｋ前画像ｉk -（Ｂ
モード画像）を撮像する。ステップＳ６では、時間Ｔk
が経過するのを待ってから、ステップＳ７へ進む。ステ
ップＳ７では、第ｋ後画像ｉk +（Ｂモード画像）を撮像
する。
を待ってから、ステップＳ５へ進む。これは、灌流を十
分に行わせて、造影剤としての気泡を臓器に十分補填す
るためである。ステップＳ５では、第ｋ前画像ｉk -（Ｂ
モード画像）を撮像する。ステップＳ６では、時間Ｔk
が経過するのを待ってから、ステップＳ７へ進む。ステ
ップＳ７では、第ｋ後画像ｉk +（Ｂモード画像）を撮像
する。
【００１３】ステップＳ８では、第ｋ前画像ｉk -の各画
素の信号強度をｇk -とし、第ｋ後画像ｉk +の各画素の信
号強度をｇk +とし、比較係数ｗ＝０．９とするとき、ｇ
k +≧ｗ・ｇk -となるか否かを全画素について判定し、前
記不等式を満たす画素が見つかれば、その画素に対応す
る灌流画像の画素の画素値を“ｋ”にする。ステップＳ
９では、ｋ＝ｅなら処理を終了し、そうでないならステ
ップＳ１０に進む。ステップＳ１０では、ｋの値を
“１”増加させ、前記ステップＳ４に戻る。
素の信号強度をｇk -とし、第ｋ後画像ｉk +の各画素の信
号強度をｇk +とし、比較係数ｗ＝０．９とするとき、ｇ
k +≧ｗ・ｇk -となるか否かを全画素について判定し、前
記不等式を満たす画素が見つかれば、その画素に対応す
る灌流画像の画素の画素値を“ｋ”にする。ステップＳ
９では、ｋ＝ｅなら処理を終了し、そうでないならステ
ップＳ１０に進む。ステップＳ１０では、ｋの値を
“１”増加させ、前記ステップＳ４に戻る。
【００１４】以上により、図４に示すように第１前画像
ｉ1 -〜第ｅ後画像ｉe +が順に撮像されながら、灌流画像
の画素値が決定されていく。灌流画像の画素値は、当該
画素における信号強度が前記不等式を満たす最短時間Ｔ
kminに応じた値であり、灌流の活発さの度合いを示す灌
流パラメータになっている。画素値“１”の画素は最も
灌流が活発な部位に対応し、画素値が“１”より大きく
なるほど灌流が不活発な部位になる。そして、画素値
“０”の画素は最も灌流が不活発な部位に当る。
ｉ1 -〜第ｅ後画像ｉe +が順に撮像されながら、灌流画像
の画素値が決定されていく。灌流画像の画素値は、当該
画素における信号強度が前記不等式を満たす最短時間Ｔ
kminに応じた値であり、灌流の活発さの度合いを示す灌
流パラメータになっている。画素値“１”の画素は最も
灌流が活発な部位に対応し、画素値が“１”より大きく
なるほど灌流が不活発な部位になる。そして、画素値
“０”の画素は最も灌流が不活発な部位に当る。
【００１５】なお、前記ステップＳ５，Ｓ７が撮像手段
に相当し、前記ステップＳ８が判定手段および画素値決
定手段に相当する。
に相当し、前記ステップＳ８が判定手段および画素値決
定手段に相当する。
【００１６】図５は、前記灌流画像表示部１０により実
行される灌流画像表示処理を示すフロー図である。ステ
ップＤ１では、灌流画像の画素ｋを、カラー表示値Ｍ
R ，ＭG ，ＭB に変換する。例えば、ｋ＝０ならカラー
表示値“緑”，ｋ＝１ならカラー表示値“青”，ｋ＝２
ならカラー表示値“黄”，ｋ＝３ならカラー表示値
“赤”，ｋ＝４ならカラー表示値“橙”，…というよう
に変換する。ステップＤ２では、次式により灌流画像の
表示値ＭR ，ＭG ，ＭB を、合成表示値ＮR ，ＮG ，Ｎ
B に変換する。 ＮR＝（１−α）ＩR＋α・ＭR ＮG＝（１−α）ＩG＋α・ＭG ＮB＝（１−α）ＩB＋α・ＭB ここで、αは重み係数であり、０≦α≦１である。ま
た、ＩR ，ＩG ，ＩB は、造影剤が存在しない状態での
画像であり、簡易には第１後画像ｉ1 +の各画素値ｇ1+を
グレースケールにより変換した白黒表示値である。ステ
ップＤ３では、前記合成表示値ＮR ，ＮG ，ＮB により
灌流画像を表示する。
行される灌流画像表示処理を示すフロー図である。ステ
ップＤ１では、灌流画像の画素ｋを、カラー表示値Ｍ
R ，ＭG ，ＭB に変換する。例えば、ｋ＝０ならカラー
表示値“緑”，ｋ＝１ならカラー表示値“青”，ｋ＝２
ならカラー表示値“黄”，ｋ＝３ならカラー表示値
“赤”，ｋ＝４ならカラー表示値“橙”，…というよう
に変換する。ステップＤ２では、次式により灌流画像の
表示値ＭR ，ＭG ，ＭB を、合成表示値ＮR ，ＮG ，Ｎ
B に変換する。 ＮR＝（１−α）ＩR＋α・ＭR ＮG＝（１−α）ＩG＋α・ＭG ＮB＝（１−α）ＩB＋α・ＭB ここで、αは重み係数であり、０≦α≦１である。ま
た、ＩR ，ＩG ，ＩB は、造影剤が存在しない状態での
画像であり、簡易には第１後画像ｉ1 +の各画素値ｇ1+を
グレースケールにより変換した白黒表示値である。ステ
ップＤ３では、前記合成表示値ＮR ，ＮG ，ＮB により
灌流画像を表示する。
【００１７】例えばα＝１とすれば、図６の（ａ）
（ｂ）に示すように、灌流の活発さの度合いの２次元分
布を表すカラー灌流画像が表示される。なお、（ａ）は
臓器が正常な灌流能力を有する場合、（ｂ）は臓器が正
常な灌流能力を有さない場合である。また、α＝０．５
とすれば、図７の（ａ）（ｂ）に示すように、灌流の活
発さの度合いの２次元分布および臓器の構造が重畳され
たカラー灌流画像が表示される。なお、（ａ）は臓器が
正常な灌流能力を有する場合、（ｂ）は臓器が正常な灌
流能力を有さない場合である。また、ｋ＝０ならα＝０
かつｋ＞０ならα＝１とすれば、灌流がある程度行われ
た部分では灌流の活発さの度合いの２次元分布を表すカ
ラー灌流画像が表示され、灌流がほとんど行われていな
い部分では臓器の構造を表す白黒画像が表示される。
（ｂ）に示すように、灌流の活発さの度合いの２次元分
布を表すカラー灌流画像が表示される。なお、（ａ）は
臓器が正常な灌流能力を有する場合、（ｂ）は臓器が正
常な灌流能力を有さない場合である。また、α＝０．５
とすれば、図７の（ａ）（ｂ）に示すように、灌流の活
発さの度合いの２次元分布および臓器の構造が重畳され
たカラー灌流画像が表示される。なお、（ａ）は臓器が
正常な灌流能力を有する場合、（ｂ）は臓器が正常な灌
流能力を有さない場合である。また、ｋ＝０ならα＝０
かつｋ＞０ならα＝１とすれば、灌流がある程度行われ
た部分では灌流の活発さの度合いの２次元分布を表すカ
ラー灌流画像が表示され、灌流がほとんど行われていな
い部分では臓器の構造を表す白黒画像が表示される。
【００１８】そして、灌流パラメータを画素値とする灌
流画像および表示値を画素値とする灌流画像の両方また
は一方は、ディスクやテープに記録され、所望によりプ
リントアウトされる。
流画像および表示値を画素値とする灌流画像の両方また
は一方は、ディスクやテープに記録され、所望によりプ
リントアウトされる。
【００１９】図８は、前記灌流３Ｄ画像表示部１１によ
り実行される灌流３Ｄ画像表示処理を示すフロー図であ
る。ステップＱ１では、連続した複数の異なる断面位置
で、それぞれ灌流画像を得る。ステップＱ２では、前記
複数の灌流画像から灌流３Ｄモデルを作成する。なお、
ここで用いる灌流画像は、灌流パラメータを画素値とす
る灌流画像または表示値を画素値とする灌流画像のいず
れでもよい。ステップＱ３では、灌流３Ｄモデルに対し
て視線方向に応じたＭＩＰ処理またはレンダリング処理
を施し、灌流３Ｄ画像を作成する。ステップＱ４では、
前記灌流３Ｄ画像を表示する。なお、前記ステップＱ２
で灌流パラメータを画素値とする灌流画像を用いた場合
には、ここで表示値を設定する必要がある。また、前記
ステップＱ２で表示値を画素値とする灌流画像を用いた
場合でも、ここで表示値を改めて設定し直してもよい。
例えば、前記ステップＱ３でＭＩＰ処理を施した場合に
は、表示値を白黒画像の表示値に設定した方が見やすく
なることがある。
り実行される灌流３Ｄ画像表示処理を示すフロー図であ
る。ステップＱ１では、連続した複数の異なる断面位置
で、それぞれ灌流画像を得る。ステップＱ２では、前記
複数の灌流画像から灌流３Ｄモデルを作成する。なお、
ここで用いる灌流画像は、灌流パラメータを画素値とす
る灌流画像または表示値を画素値とする灌流画像のいず
れでもよい。ステップＱ３では、灌流３Ｄモデルに対し
て視線方向に応じたＭＩＰ処理またはレンダリング処理
を施し、灌流３Ｄ画像を作成する。ステップＱ４では、
前記灌流３Ｄ画像を表示する。なお、前記ステップＱ２
で灌流パラメータを画素値とする灌流画像を用いた場合
には、ここで表示値を設定する必要がある。また、前記
ステップＱ２で表示値を画素値とする灌流画像を用いた
場合でも、ここで表示値を改めて設定し直してもよい。
例えば、前記ステップＱ３でＭＩＰ処理を施した場合に
は、表示値を白黒画像の表示値に設定した方が見やすく
なることがある。
【００２０】以上の超音波診断装置１００によれば、臓
器の灌流の活発さの度合い、すなわち、灌流パラメータ
を定量的に測定でき、その２次元分布を表す灌流画像お
よび３次元分布を表す灌流３Ｄ画像を生成し、表示でき
る。従って、臓器の疾病を従来と異なった観点で診断で
きるようになる。
器の灌流の活発さの度合い、すなわち、灌流パラメータ
を定量的に測定でき、その２次元分布を表す灌流画像お
よび３次元分布を表す灌流３Ｄ画像を生成し、表示でき
る。従って、臓器の疾病を従来と異なった観点で診断で
きるようになる。
【００２１】
【発明の効果】本発明の灌流パラメータ測定方法によれ
ば、灌流の活発さの度合いを表す灌流パラメータを定量
的に測定することが出来る。また、本発明の灌流画像生
成方法および超音波診断装置によれば、灌流の活発さの
度合いの２次元分布を表す灌流画像を生成することが出
来る。
ば、灌流の活発さの度合いを表す灌流パラメータを定量
的に測定することが出来る。また、本発明の灌流画像生
成方法および超音波診断装置によれば、灌流の活発さの
度合いの２次元分布を表す灌流画像を生成することが出
来る。
【図１】本発明の一実施形態にかかる超音波診断装置を
示す構成図である。
示す構成図である。
【図２】図１の超音波診断装置の灌流画像生成部により
実行される灌流画像生成処理を示すフロー図である。
実行される灌流画像生成処理を示すフロー図である。
【図３】種々のＴｋテーブルの例示図である。
【図４】信号強度画像の撮像タイミングを示すタイムチ
ャートである。
ャートである。
【図５】図１の超音波診断装置の灌流画像表示部により
実行される灌流画像表示処理を示すフロー図である。
実行される灌流画像表示処理を示すフロー図である。
【図６】表示された灌流画像の例示図である。
【図７】表示された灌流画像の別の表示態様の例示図で
ある。
ある。
【図８】図１の超音波診断装置の灌流３Ｄ画像表示部に
より実行される灌流３Ｄ画像表示処理を示すフロー図で
ある。
より実行される灌流３Ｄ画像表示処理を示すフロー図で
ある。
１ 超音波探触子 ２ 送受信制御部 ３ Ｂモード処理部 ６ 中央処理部 ７ ＣＲＴ ９ 灌流画像生成部 １０ 灌流画像表示部 １１ 灌流３Ｄ画像表示部 １００ 超音波診断装置
Claims (5)
- 【請求項１】 被検体内の測定点に造影剤としての気泡
が十分存在する状態で該測定点からの超音波エコーの信
号強度ｇk -を測定し、その測定のために送信した超音波
ビームにより前記気泡を破壊し、時間Ｔk 後に前記測定
点からの超音波エコーの信号強度ｇk +を再び測定し、比
較係数ｗを０＜ｗ＜１とするとき、 ｇk +≧ｗ・ｇk - となるか否かを判定することを、前記時間Ｔk を変えな
がら繰り返して、前記不等式を満たす最短時間Ｔkminを
求め、その最短時間Ｔkminに基づいて灌流の活発さの度
合いを表すパラメータを決定することを特徴とする灌流
パラメータ測定方法。 - 【請求項２】 被検体内に造影剤としての気泡が十分存
在する状態で超音波エコーの信号強度画像ｉk -を撮像
し、その撮像のために送信した超音波ビームにより前記
気泡を破壊し、時間Ｔk 後に超音波エコーの信号強度画
像ｉk +を再び撮像し、前記信号強度画像ｉk -の各画素の
信号強度をｇk -とし、前記信号強度画像ｉk +の各画素の
信号強度をｇk +とし、比較係数ｗを０＜ｗ＜１とすると
き、 ｇk +≧ｗ・ｇk - となるか否かを判定することを、前記時間Ｔk を変えな
がら繰り返して、各画素の信号強度が前記不等式を満た
す最短時間Ｔkminを求め、その最短時間Ｔkminに基づい
て灌流画像の各画素の画素値を決定することを特徴とす
る灌流画像生成方法。 - 【請求項３】 請求項２に記載の灌流画像生成方法にお
いて、実質的に灌流がされない程度に前記時間Ｔk を短
くして信号強度画像ｉk +’を撮像し、その信号強度画像
ｉk +’の各画素の画素値と前記灌流画像の画素値とに基
づいて灌流画像の各画素の表示値を決定することを特徴
とする灌流画像生成方法。 - 【請求項４】 請求項２または請求項３に記載の灌流画
像生成方法において、位置の異なる複数の灌流画像を基
に灌流３次元画像を生成することを特徴とする灌流画像
生成方法。 - 【請求項５】 被検体内に造影剤として存在する気泡を
破壊しうる超音波ビームにより超音波エコーの第１の信
号強度画像ｉk -を撮像し続いて時間Ｔk 後に第２の信号
強度画像ｉk +を撮像する撮像手段と、前記第１の信号強
度画像ｉk -の各画素の信号強度をｇk -とし，前記第２の
信号強度画像ｉk +の各画素の信号強度をｇk +とし，比較
係数ｗを０＜ｗ＜１とするとき，ｇk +≧ｗ・ｇk -となる
か否かを判定する判定手段と、前記撮像手段と前記判定
手段とにより前記時間Ｔk を変えながら前記撮像と前記
判定とを繰り返して，各画素の信号強度が前記不等式を
満たす最短時間Ｔkminを求め，その最短時間Ｔkminに基
づいて灌流画像の各画素の画素値を決定する画素値決定
手段とを具備したことを特徴とする超音波診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25959897A JP3712510B2 (ja) | 1997-09-25 | 1997-09-25 | 超音波診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25959897A JP3712510B2 (ja) | 1997-09-25 | 1997-09-25 | 超音波診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH1189839A true JPH1189839A (ja) | 1999-04-06 |
| JP3712510B2 JP3712510B2 (ja) | 2005-11-02 |
Family
ID=17336332
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP25959897A Expired - Fee Related JP3712510B2 (ja) | 1997-09-25 | 1997-09-25 | 超音波診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3712510B2 (ja) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7302850B2 (en) | 2003-04-28 | 2007-12-04 | Kabushiki Kaisha Toshiba | Ultrasonic diagnostic apparatus and image processing apparatus |
| JP2008295691A (ja) * | 2007-05-31 | 2008-12-11 | Ge Medical Systems Global Technology Co Llc | 超音波撮像装置 |
| JP2009233408A (ja) * | 2009-07-21 | 2009-10-15 | Toshiba Corp | 超音波診断装置 |
-
1997
- 1997-09-25 JP JP25959897A patent/JP3712510B2/ja not_active Expired - Fee Related
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7302850B2 (en) | 2003-04-28 | 2007-12-04 | Kabushiki Kaisha Toshiba | Ultrasonic diagnostic apparatus and image processing apparatus |
| US7886603B2 (en) | 2003-04-28 | 2011-02-15 | Kabushiki Kaisha Toshiba | Ultrasonic diagnostic apparatus and image processing apparatus |
| US8499636B2 (en) | 2003-04-28 | 2013-08-06 | Kabushiki Kaisha Toshiba | Ultrasonic diagnostic apparatus and image processing apparatus |
| JP2008295691A (ja) * | 2007-05-31 | 2008-12-11 | Ge Medical Systems Global Technology Co Llc | 超音波撮像装置 |
| JP2009233408A (ja) * | 2009-07-21 | 2009-10-15 | Toshiba Corp | 超音波診断装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3712510B2 (ja) | 2005-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7302850B2 (en) | Ultrasonic diagnostic apparatus and image processing apparatus | |
| US5873829A (en) | Diagnostic ultrasound system using harmonic echo imaging | |
| US7713208B2 (en) | Ultrasound diagnostic apparatus | |
| CN101264023B (zh) | 超声波诊断装置、图像处理装置和图像处理方法 | |
| US4790321A (en) | Method of displaying stream lines of an inhomogeneous flowing medium and a device therefor | |
| US20040126007A1 (en) | System and method for improved multiple-dimension image displays | |
| US6413218B1 (en) | Medical diagnostic ultrasound imaging system and method for determining an acoustic output parameter of a transmitted ultrasonic beam | |
| CN102274046A (zh) | 超声波诊断装置、超声波图像处理装置以及医用图像诊断装置 | |
| US11701091B2 (en) | Ultrasound analysis apparatus and method for tissue elasticity and viscosity based on the hormonic signals | |
| JP2007044499A (ja) | 超音波診断装置及び超音波画像処理プログラム | |
| JP2014028029A (ja) | 超音波診断装置及び超音波診断装置制御プログラム | |
| JP2004202229A (ja) | 造影剤時間強度曲線分析のための方法及び装置 | |
| JP2012176232A (ja) | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム | |
| JPH0779981A (ja) | 超音波診断装置 | |
| JP3961209B2 (ja) | 超音波撮影装置 | |
| JP2021529035A (ja) | 経弁圧の定量化を行うための方法及びシステム | |
| JP4521204B2 (ja) | 超音波診断装置、画像処理装置及び超音波画像撮影方法 | |
| WO2020082270A1 (zh) | 一种成像方法以及成像*** | |
| JP3712510B2 (ja) | 超音波診断装置 | |
| JP2002301070A (ja) | 超音波撮像方法及び装置 | |
| JP3405578B2 (ja) | 超音波診断装置 | |
| JP2013244138A (ja) | 超音波診断装置および音速表示方法 | |
| JP7297485B2 (ja) | 超音波診断装置、医用画像処理装置および医用画像処理プログラム | |
| JP2023048353A (ja) | 超音波診断装置及びプログラム | |
| JP5454901B2 (ja) | 超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050426 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050502 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050624 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050809 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050817 |
|
| R150 | Certificate of patent (=grant) or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |