JP3405578B2 - 超音波診断装置 - Google Patents
超音波診断装置Info
- Publication number
- JP3405578B2 JP3405578B2 JP32348493A JP32348493A JP3405578B2 JP 3405578 B2 JP3405578 B2 JP 3405578B2 JP 32348493 A JP32348493 A JP 32348493A JP 32348493 A JP32348493 A JP 32348493A JP 3405578 B2 JP3405578 B2 JP 3405578B2
- Authority
- JP
- Japan
- Prior art keywords
- time
- interest
- change
- region
- luminance signal
- 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.)
- Expired - Lifetime
Links
Description
【０００１】
【産業上の利用分野】本発明は、被検体内を超音波ビー
ムで走査することにより得られた受信信号の包絡線を検
波し、この検波信号を用いて断層像を表示する超音波診
断装置に関する。
ムで走査することにより得られた受信信号の包絡線を検
波し、この検波信号を用いて断層像を表示する超音波診
断装置に関する。
【０００２】
【従来の技術】近年、超音波診断の分野では、Ｘ線診断
における造影剤と同様の機能としての超音波反射特性が
良好なコントラスト剤の開発が進み、実用域に入ってい
る。このコントラスト剤を例えば血管に注入すると、Ｂ
モード画像に血流の様子がその周辺組織と異なるコント
ラストで鮮明に表れる。
における造影剤と同様の機能としての超音波反射特性が
良好なコントラスト剤の開発が進み、実用域に入ってい
る。このコントラスト剤を例えば血管に注入すると、Ｂ
モード画像に血流の様子がその周辺組織と異なるコント
ラストで鮮明に表れる。
【０００３】ところで、血流が関心部分に流入してから
流出するまでの時間的変化には貴重な情報が含まれてい
る。しかし従来の超音波診断装置では、上記時間的変化
を計算する機能を備えていないため、検査終了後に複数
枚のＢモード画像データを時間的変化を計算する処理装
置に転送して、そこで計算するようにしていた。このた
め、診断結果が出されるまでに長時間必要であった。
流出するまでの時間的変化には貴重な情報が含まれてい
る。しかし従来の超音波診断装置では、上記時間的変化
を計算する機能を備えていないため、検査終了後に複数
枚のＢモード画像データを時間的変化を計算する処理装
置に転送して、そこで計算するようにしていた。このた
め、診断結果が出されるまでに長時間必要であった。
【０００４】
【発明が解決しようとする課題】そこで、本発明は、上
述した事情に対処するべく、Ｂモード画像の生成に並行
して輝度の時間的変化をリアルタイムで計算する超音波
診断装置を提供することを目的とする。
述した事情に対処するべく、Ｂモード画像の生成に並行
して輝度の時間的変化をリアルタイムで計算する超音波
診断装置を提供することを目的とする。
【０００５】
【課題を解決するための手段】本発明の第１局面は、コ
ントラスト剤を注入した状態で被検体内を超音波ビーム
で走査することにより得られた受信信号の包絡線を検波
し、この検波結果を輝度信号として断層像を表示する超
音波診断装置において、前記断層像の関心領域内の輝度
信号の最大値及び最小値を選択し、前記最大値の時間変
化を表すグラフと前記最小値の時間変化を表すグラフと
を生成する手段と、前記最大値の時間変化を表すグラフ
と前記最小値の時間変化を表すグラフとを前記断層像と
共に表示する手段とを具備する。本発明の第２局面は、
コントラスト剤を注入した状態で被検体内を超音波ビー
ムで走査することにより得られた受信信号の包絡線を検
波し、この検波結果を輝度信号として断層像を表示する
超音波診断装置において、前記断層像の関心領域内の輝
度信号の平均値とこの平均値の時間変化率とを計算し、
前記平均値の時間変化率に応じた線分長を有する線分を
前記平均値と時間に応じて配列させたグラフを生成する
手段と、前記グラフを前記断層像と共に表示する手段と
を具備する。
ントラスト剤を注入した状態で被検体内を超音波ビーム
で走査することにより得られた受信信号の包絡線を検波
し、この検波結果を輝度信号として断層像を表示する超
音波診断装置において、前記断層像の関心領域内の輝度
信号の最大値及び最小値を選択し、前記最大値の時間変
化を表すグラフと前記最小値の時間変化を表すグラフと
を生成する手段と、前記最大値の時間変化を表すグラフ
と前記最小値の時間変化を表すグラフとを前記断層像と
共に表示する手段とを具備する。本発明の第２局面は、
コントラスト剤を注入した状態で被検体内を超音波ビー
ムで走査することにより得られた受信信号の包絡線を検
波し、この検波結果を輝度信号として断層像を表示する
超音波診断装置において、前記断層像の関心領域内の輝
度信号の平均値とこの平均値の時間変化率とを計算し、
前記平均値の時間変化率に応じた線分長を有する線分を
前記平均値と時間に応じて配列させたグラフを生成する
手段と、前記グラフを前記断層像と共に表示する手段と
を具備する。
【０００６】
【０００７】
【作用】本発明によれば、関心領域の血流状態を良好に
診断することができる。
診断することができる。
【０００８】
【０００９】
【実施例】図面を参照して本発明の一実施例について説
明する。図１は第１実施例のブロック図である。セクタ
電子走査型超音波プローブ１は複数の圧電振動子が一次
元に配列される。このプローブ１には送受信制御と受信
信号の検波とを兼ねたアナログ部２が接続される。この
アナログ部２にはパルサ２２が含まれる。
明する。図１は第１実施例のブロック図である。セクタ
電子走査型超音波プローブ１は複数の圧電振動子が一次
元に配列される。このプローブ１には送受信制御と受信
信号の検波とを兼ねたアナログ部２が接続される。この
アナログ部２にはパルサ２２が含まれる。
【００１０】パルサ２２はディレーライン２４から供給
される各圧電振動子毎の遅延情報にしたがって異なるタ
イミングで各圧電振動子に駆動パルスを供給する。この
駆動パルスによりプローブ１から超音波ビームが上記遅
延情報に応じた方向に送信される。この遅延情報をラス
タ毎に変化させることにより超音波ビームは扇状に移動
する。超音波ビームは被検体内の音響インピーダンスの
境界で反射する。この反射波は、プローブ１の各振動子
で受信される。各受信信号はプリアンプ２１を介してデ
ィレーライン２４に個別に送られ、そこで所定の焦点深
度に応じた遅延時間を与えられた後、加算器２５で加算
される。この加算信号は、検波器２６に送られ、そこで
対数増幅および包絡線が検波される。この検波結果は輝
度信号として一般にはディジタルスキャンコンバータと
呼ばれる走査変換手段３と共に、信号レベル時間変化抽
出器６に送られる。
される各圧電振動子毎の遅延情報にしたがって異なるタ
イミングで各圧電振動子に駆動パルスを供給する。この
駆動パルスによりプローブ１から超音波ビームが上記遅
延情報に応じた方向に送信される。この遅延情報をラス
タ毎に変化させることにより超音波ビームは扇状に移動
する。超音波ビームは被検体内の音響インピーダンスの
境界で反射する。この反射波は、プローブ１の各振動子
で受信される。各受信信号はプリアンプ２１を介してデ
ィレーライン２４に個別に送られ、そこで所定の焦点深
度に応じた遅延時間を与えられた後、加算器２５で加算
される。この加算信号は、検波器２６に送られ、そこで
対数増幅および包絡線が検波される。この検波結果は輝
度信号として一般にはディジタルスキャンコンバータと
呼ばれる走査変換手段３と共に、信号レベル時間変化抽
出器６に送られる。
【００１１】信号レベル時間変化抽出器６には、超音波
ビームによる走査の繰り返しに応じて図示しないマウス
やトラックボール等の関心領域設定部により設定された
関心領域に相当する輝度信号が２次元走査面の走査毎に
検波器２６から連続的に供給されてくる。信号レベル時
間変化抽出器６には縦軸を輝度レベル、横軸を時間とし
たグラフ領域を内部メモリに備えている。関心領域の輝
度信号が供給される毎に、その輝度信号の輝度レベルお
よび走査時刻に応じたグラフ領域の位置に点が分布され
る。このグラフ情報は、次の輝度信号が供給されるまで
に作成され、走査変換手段３に出力される。
ビームによる走査の繰り返しに応じて図示しないマウス
やトラックボール等の関心領域設定部により設定された
関心領域に相当する輝度信号が２次元走査面の走査毎に
検波器２６から連続的に供給されてくる。信号レベル時
間変化抽出器６には縦軸を輝度レベル、横軸を時間とし
たグラフ領域を内部メモリに備えている。関心領域の輝
度信号が供給される毎に、その輝度信号の輝度レベルお
よび走査時刻に応じたグラフ領域の位置に点が分布され
る。このグラフ情報は、次の輝度信号が供給されるまで
に作成され、走査変換手段３に出力される。
【００１２】走査変換手段３には表示手段５の表示画面
に対応する２次元メモリが備えられている。この２次元
メモリには検波器２６からの輝度信号を記憶する領域
と、信号レベル時間変化抽出器６からのグラフ情報を記
憶する領域とが予め区別して割り当てられている。走査
変換手段３は輝度信号やグラフ情報を一旦ディジタル信
号として記憶する。走査変換手段３から表示手段５の走
査方式にしたがって読み出された輝度信号やグラフ情報
は、Ｄ／Ａ変換手段４を介してモニタとしての表示手段
５に送られ、そこで関心領域に関する輝度信号の強度
（レベル）の時間的変化がグラフとしてＢモード画像
（断層像）と別々に表示される。
に対応する２次元メモリが備えられている。この２次元
メモリには検波器２６からの輝度信号を記憶する領域
と、信号レベル時間変化抽出器６からのグラフ情報を記
憶する領域とが予め区別して割り当てられている。走査
変換手段３は輝度信号やグラフ情報を一旦ディジタル信
号として記憶する。走査変換手段３から表示手段５の走
査方式にしたがって読み出された輝度信号やグラフ情報
は、Ｄ／Ａ変換手段４を介してモニタとしての表示手段
５に送られ、そこで関心領域に関する輝度信号の強度
（レベル）の時間的変化がグラフとしてＢモード画像
（断層像）と別々に表示される。
【００１３】次に本実施例の動作について説明する。Ｘ
線診断における造影剤と同様の機能としての超音波反射
特性が良好なコントラスト剤を例えば血管に注入した
後、以下のような動作が開始される。
線診断における造影剤と同様の機能としての超音波反射
特性が良好なコントラスト剤を例えば血管に注入した
後、以下のような動作が開始される。
【００１４】アナログ部２のパルサ２２からプローブ１
１の各圧電振動子にディレーライン２４の遅延制御にし
たがって駆動電圧がされる。これにより超音波ビームが
被検体内に送信される。生体中で反射してきた反射波
は、プローブ１１で受信される。この受信信号は、アナ
ログ部２のプリアンプ２１、ディレーライン２４、加算
器２５を順に介して検波器２６に送られる。受信信号は
その包絡線を検波器２６で検波され、図示しないアナロ
グディジタル変換器を介して輝度信号として走査変換手
段３に送られる。走査変換手段３では輝度信号は当該ラ
スタの方位角にしたがってメモリ領域に記憶される。
１の各圧電振動子にディレーライン２４の遅延制御にし
たがって駆動電圧がされる。これにより超音波ビームが
被検体内に送信される。生体中で反射してきた反射波
は、プローブ１１で受信される。この受信信号は、アナ
ログ部２のプリアンプ２１、ディレーライン２４、加算
器２５を順に介して検波器２６に送られる。受信信号は
その包絡線を検波器２６で検波され、図示しないアナロ
グディジタル変換器を介して輝度信号として走査変換手
段３に送られる。走査変換手段３では輝度信号は当該ラ
スタの方位角にしたがってメモリ領域に記憶される。
【００１５】検波器２６からの輝度信号は信号レベル時
間変化抽出器６にも供給される。信号レベル時間変化抽
出器６では、関心領域設定部の操作により設定された少
なくとも１点（１画素）の関心領域の輝度信号の時間経
過に伴うレベル変化をグラフとして作成するが、複数点
（複数画素）を含む関心領域が設定された場合、本実施
例では、関心領域内全点の平均値の時間的変化を求める
こと、関心領域内の最大値及び／又は最小値の時間的変
化を求めること、関心領域内全点の平均値の所定の短時
間間隔における変化率の時間的変化を求めることのいず
れかにより時間経過に伴う輝度信号のレベル変化を求め
る。複数点を含む関心領域内全点の平均値を求める場
合、各点に対し窓関数による重み付けを行うこともあ
る。
間変化抽出器６にも供給される。信号レベル時間変化抽
出器６では、関心領域設定部の操作により設定された少
なくとも１点（１画素）の関心領域の輝度信号の時間経
過に伴うレベル変化をグラフとして作成するが、複数点
（複数画素）を含む関心領域が設定された場合、本実施
例では、関心領域内全点の平均値の時間的変化を求める
こと、関心領域内の最大値及び／又は最小値の時間的変
化を求めること、関心領域内全点の平均値の所定の短時
間間隔における変化率の時間的変化を求めることのいず
れかにより時間経過に伴う輝度信号のレベル変化を求め
る。複数点を含む関心領域内全点の平均値を求める場
合、各点に対し窓関数による重み付けを行うこともあ
る。
【００１６】また、信号レベル時間変化抽出器６は、１
心拍期間のグラフから輝度レベルが最大値を示す時刻か
らその最大値の半値になる時刻までの時間幅をレベル変
化を示す数値情報として求める。
心拍期間のグラフから輝度レベルが最大値を示す時刻か
らその最大値の半値になる時刻までの時間幅をレベル変
化を示す数値情報として求める。
【００１７】図２には関心領域内全点の平均値の時間的
変化を求める場合の信号レベル時間変化抽出器６による
２次元走査面の１回の走査動作に対応する動作手順が順
番に示されている。なお、この動作は２次元走査面を一
通り走査する毎に繰り返される。
変化を求める場合の信号レベル時間変化抽出器６による
２次元走査面の１回の走査動作に対応する動作手順が順
番に示されている。なお、この動作は２次元走査面を一
通り走査する毎に繰り返される。
【００１８】ステップ（Ｓ１）で、１回の２次元走査面
分の輝度信号が検波器２６から信号レベル時間変化抽出
器６に供給される。ステップ（Ｓ２）では、１回の２次
元走査面分の輝度信号の中の表示手段５に表示されたＢ
モード画像上に関心領域設定部の操作により設定された
少なくとも１点（１画素）の関心領域の輝度信号だけが
信号レベル時間変化抽出器６の内部メモリに格納され、
関心領域以外の輝度信号は破棄される。
分の輝度信号が検波器２６から信号レベル時間変化抽出
器６に供給される。ステップ（Ｓ２）では、１回の２次
元走査面分の輝度信号の中の表示手段５に表示されたＢ
モード画像上に関心領域設定部の操作により設定された
少なくとも１点（１画素）の関心領域の輝度信号だけが
信号レベル時間変化抽出器６の内部メモリに格納され、
関心領域以外の輝度信号は破棄される。
【００１９】ステップ（Ｓ３）では、内部メモリに格納
されている関心領域内の輝度信号のレベルの平均値が以
下の（１）式にしたがって計算される。なお、Ｉ´を輝
度レベルの平均値とし、ｘ，ｙそれぞれｘ方向，ｙ方向
の画素の座標（関心領域ｘ0〜ｘ1 ，ｙ0 〜ｙ1 の範囲
内の座標）とし、Ｉ(x,y) は座標（ｘ，ｙ）の輝度レベ
ルとする。
されている関心領域内の輝度信号のレベルの平均値が以
下の（１）式にしたがって計算される。なお、Ｉ´を輝
度レベルの平均値とし、ｘ，ｙそれぞれｘ方向，ｙ方向
の画素の座標（関心領域ｘ0〜ｘ1 ，ｙ0 〜ｙ1 の範囲
内の座標）とし、Ｉ(x,y) は座標（ｘ，ｙ）の輝度レベ
ルとする。
【００２０】
【数１】
【００２１】この平均値はステップ（Ｓ４）で内部メモ
リに記憶された後、ステップ（Ｓ５）で、所定時間前か
ら直前までにすでに作成されているグラフ情報に、当該
平均値のレベル及び走査時刻に応じた位置に点が追加さ
れて、最新のグラフ情報が作成される。この最新のグラ
フ情報は、ステップ（Ｓ６）で走査変換手段３に送ら
れ、走査変換手段３のメモリ領域における輝度信号を記
憶する領域以外の専用領域に記憶される。
リに記憶された後、ステップ（Ｓ５）で、所定時間前か
ら直前までにすでに作成されているグラフ情報に、当該
平均値のレベル及び走査時刻に応じた位置に点が追加さ
れて、最新のグラフ情報が作成される。この最新のグラ
フ情報は、ステップ（Ｓ６）で走査変換手段３に送ら
れ、走査変換手段３のメモリ領域における輝度信号を記
憶する領域以外の専用領域に記憶される。
【００２２】走査変換手段３では、同じ２次元走査によ
る受信信号から作成されたＢモード画像とグラフ情報を
順次走査によりＤ／Ａ変換手段４に読出し、図３（ａ）
に示すように表示手段５の１画面内にＢモード画像とグ
ラフ情報を別々に表示させる。なお図３（ｂ）には図３
（ａ）のグラフ情報の拡大図である。
る受信信号から作成されたＢモード画像とグラフ情報を
順次走査によりＤ／Ａ変換手段４に読出し、図３（ａ）
に示すように表示手段５の１画面内にＢモード画像とグ
ラフ情報を別々に表示させる。なお図３（ｂ）には図３
（ａ）のグラフ情報の拡大図である。
【００２３】このように関心領域の輝度信号の時間変化
曲線がＢモード画像と共にリアルタイムで表示される。
図４には関心領域内の最大値及び／又は最小値の時間的
変化を求める場合の信号レベル時間変化抽出器６による
２次元走査面の１回の走査動作に対応する動作手順が順
番に示されている。なお、この動作は２次元走査面を一
通り走査する毎に繰り返される。図４において図２と同
じステップは同符号を付して説明は省略する。ステップ
（Ｓ１３）では、内部メモリに格納されている関心領域
内の全ての点の輝度信号の中の最大値及び最小値が選択
される。これらの最大値及び最小値はステップ（Ｓ
４）、ステップ（Ｓ５）、ステップ（Ｓ６）を経過し
て、最大値及び最小値の時間変化を現すグラフ情報とし
て走査変換手段３に送られ、図３（ａ）と同様にＢモー
ド画像と共に同時表示される。図５には表示手段５に表
示されるグラフ情報を示している。
曲線がＢモード画像と共にリアルタイムで表示される。
図４には関心領域内の最大値及び／又は最小値の時間的
変化を求める場合の信号レベル時間変化抽出器６による
２次元走査面の１回の走査動作に対応する動作手順が順
番に示されている。なお、この動作は２次元走査面を一
通り走査する毎に繰り返される。図４において図２と同
じステップは同符号を付して説明は省略する。ステップ
（Ｓ１３）では、内部メモリに格納されている関心領域
内の全ての点の輝度信号の中の最大値及び最小値が選択
される。これらの最大値及び最小値はステップ（Ｓ
４）、ステップ（Ｓ５）、ステップ（Ｓ６）を経過し
て、最大値及び最小値の時間変化を現すグラフ情報とし
て走査変換手段３に送られ、図３（ａ）と同様にＢモー
ド画像と共に同時表示される。図５には表示手段５に表
示されるグラフ情報を示している。
【００２４】図６には関心領域内全点の平均値の所定の
短時間間隔における変化率の時間的変化を求める場合の
信号レベル時間変化抽出器６による２次元走査面の１回
の走査動作に対応する動作手順が順番に示されている。
なお、この動作は２次元走査面を一通り走査する毎に繰
り返される。図６において図２と同じステップは同符号
を付して説明は省略する。
短時間間隔における変化率の時間的変化を求める場合の
信号レベル時間変化抽出器６による２次元走査面の１回
の走査動作に対応する動作手順が順番に示されている。
なお、この動作は２次元走査面を一通り走査する毎に繰
り返される。図６において図２と同じステップは同符号
を付して説明は省略する。
【００２５】ステップ（Ｓ２５）では、以下の（２）に
したがって内部メモリに格納されている関心領域内の全
ての点の輝度信号の平均値と前回または所定時間前の走
査による輝度信号の平均値との差ｄＩが求められ、この
ｄＩが今回から前回または所定時間前の走査時刻までの
所定の短時間間隔ｄｔで割り算されることにより、その
変化率△Ｉが計算される。
したがって内部メモリに格納されている関心領域内の全
ての点の輝度信号の平均値と前回または所定時間前の走
査による輝度信号の平均値との差ｄＩが求められ、この
ｄＩが今回から前回または所定時間前の走査時刻までの
所定の短時間間隔ｄｔで割り算されることにより、その
変化率△Ｉが計算される。
【００２６】
△Ｉ＝ｄＩ／ｄｔ ・・・（２）
この変化率はステップ（Ｓ４）と同じステップ（Ｓ２
６）、ステップ（Ｓ５）、ステップ（Ｓ６）を経過し
て、所定の短時間間隔における変化率の大きさに応じた
長さの線分で、その線分の上点が今回の平均値に一致す
るグラフ情報として、図７に示すように表示手段５に表
示される。勿論、変化率の時間変化も図３（ａ）と同様
にＢモード画像と共に同時表示される。
６）、ステップ（Ｓ５）、ステップ（Ｓ６）を経過し
て、所定の短時間間隔における変化率の大きさに応じた
長さの線分で、その線分の上点が今回の平均値に一致す
るグラフ情報として、図７に示すように表示手段５に表
示される。勿論、変化率の時間変化も図３（ａ）と同様
にＢモード画像と共に同時表示される。
【００２７】また、信号レベル時間変化抽出器６は、１
心拍期間のグラフから輝度レベルが最大値を示す時刻か
らその最大値の半値になる時刻までの半値時間幅をレベ
ル変化を示す数値情報として求める。図８はこの半値時
間幅を説明する図であり、１心拍期間のグラフを示して
いる。信号レベル時間変化抽出器６は、１心拍期間が経
過する毎に、そのグラフの中の輝度レベルが最大値を示
す時刻ｔ1 と、その最大値の半値になる時刻ｔ2 を求
め、時刻ｔ1 からｔ2 に至る半値時間幅ｔw を数値情報
として計算する。このように求めた半値時間幅ｔw は、
走査変換手段３を介して、Ｂモード画像及びグラフ情報
と共に表示手段５に表示される。
心拍期間のグラフから輝度レベルが最大値を示す時刻か
らその最大値の半値になる時刻までの半値時間幅をレベ
ル変化を示す数値情報として求める。図８はこの半値時
間幅を説明する図であり、１心拍期間のグラフを示して
いる。信号レベル時間変化抽出器６は、１心拍期間が経
過する毎に、そのグラフの中の輝度レベルが最大値を示
す時刻ｔ1 と、その最大値の半値になる時刻ｔ2 を求
め、時刻ｔ1 からｔ2 に至る半値時間幅ｔw を数値情報
として計算する。このように求めた半値時間幅ｔw は、
走査変換手段３を介して、Ｂモード画像及びグラフ情報
と共に表示手段５に表示される。
【００２８】このように本実施例によれば、Ｂモード画
像と同じ輝度信号を用いて、関心領域の輝度信号の輝度
レベルの時間的変化がグラフ情報として、また半値時間
幅が数値情報としてＢモード画像と同じ受信信号から求
められ、Ｂモード画像と共にリアルタイムで表示され
る。したがって、超音波反射特性が良好なコントラスト
剤を血管に注入した場合、コントラスト剤が関心領域に
流入するに応じて輝度レベルが序々に高くなり、コント
ラスト剤が関心領域から流出するに応じて輝度レベルが
序々に低くなる様子を、リアルタイムで且つ定量的に診
断することができる。
像と同じ輝度信号を用いて、関心領域の輝度信号の輝度
レベルの時間的変化がグラフ情報として、また半値時間
幅が数値情報としてＢモード画像と同じ受信信号から求
められ、Ｂモード画像と共にリアルタイムで表示され
る。したがって、超音波反射特性が良好なコントラスト
剤を血管に注入した場合、コントラスト剤が関心領域に
流入するに応じて輝度レベルが序々に高くなり、コント
ラスト剤が関心領域から流出するに応じて輝度レベルが
序々に低くなる様子を、リアルタイムで且つ定量的に診
断することができる。
【００２９】次に第２実施例について説明する。本実施
例は、輝度信号の時間的変化をＢモード画像の各点毎に
計算し、この計算した結果（これを以下「インデック
ス」という）をそのレベルに応じてカラー情報に変換
し、カラーの２次元画像としてＢモード画像に重ねて表
示することを特徴とする。インデックスを得る手段は前
述の通り、複数画素の平均をその中心画素のインデック
スとすることも可能である。
例は、輝度信号の時間的変化をＢモード画像の各点毎に
計算し、この計算した結果（これを以下「インデック
ス」という）をそのレベルに応じてカラー情報に変換
し、カラーの２次元画像としてＢモード画像に重ねて表
示することを特徴とする。インデックスを得る手段は前
述の通り、複数画素の平均をその中心画素のインデック
スとすることも可能である。
【００３０】図９は本実施例のブロック図を示す。な
お、図９において図１と同じ部分は同符号を付して説明
は省略する。検波器２６からの輝度信号はディジタル信
号としてインデックス演算器１０に供給される。インデ
ックス演算器１０では、輝度レベルと前回または所定時
間前の走査による輝度信号の輝度レベルとの差が２次元
走査面の全ての点について各々求められ、この差が今回
から前回または所定時間前の走査時刻までの所定の短時
間間隔で割り算されることにより、その変化率が２次元
走査面の全ての点について各々計算される。
お、図９において図１と同じ部分は同符号を付して説明
は省略する。検波器２６からの輝度信号はディジタル信
号としてインデックス演算器１０に供給される。インデ
ックス演算器１０では、輝度レベルと前回または所定時
間前の走査による輝度信号の輝度レベルとの差が２次元
走査面の全ての点について各々求められ、この差が今回
から前回または所定時間前の走査時刻までの所定の短時
間間隔で割り算されることにより、その変化率が２次元
走査面の全ての点について各々計算される。
【００３１】このインデックス演算器１０で演算された
２次元走査面内の各点の変化率は、走査変換手段３と同
じ走査変換手段３２に送られ、一旦記憶される。そし
て、走査変換手段３の読出しと同期して走査変換手段３
２から読み出され、カラー情報変換手段１２に送られ、
そこで各点毎に又は微小領域毎にその変化率（微小領域
内の各点の変化率合計）の相違に応じて異なる色彩又は
彩度のカラー情報に変換された後、マルチプレクサ（Ｍ
ＰＸ）７に送られる。マルチプレクサ７には、走査変換
手段３の出力も供給されている。図示しない表示選択入
力手段の操作による変化率の２次元カラー画像とＢモー
ド画像を重ねて表示するモードと、２次元カラー画像ま
たはＢモード画像だけを表示するモードのいずれかのモ
ードにしたがって、変化率の２次元カラー画像とＢモー
ド画像とが１枚の画像に合成されて出力されるか、また
は変化率の２次元カラー画像またはＢモード画像だけが
出力され、Ｄ／Ａ変換手段４を介して表示手段５に表示
される。図１０は変化率の２次元カラー画像とＢモード
画像との合成画像の表示例を示している。
２次元走査面内の各点の変化率は、走査変換手段３と同
じ走査変換手段３２に送られ、一旦記憶される。そし
て、走査変換手段３の読出しと同期して走査変換手段３
２から読み出され、カラー情報変換手段１２に送られ、
そこで各点毎に又は微小領域毎にその変化率（微小領域
内の各点の変化率合計）の相違に応じて異なる色彩又は
彩度のカラー情報に変換された後、マルチプレクサ（Ｍ
ＰＸ）７に送られる。マルチプレクサ７には、走査変換
手段３の出力も供給されている。図示しない表示選択入
力手段の操作による変化率の２次元カラー画像とＢモー
ド画像を重ねて表示するモードと、２次元カラー画像ま
たはＢモード画像だけを表示するモードのいずれかのモ
ードにしたがって、変化率の２次元カラー画像とＢモー
ド画像とが１枚の画像に合成されて出力されるか、また
は変化率の２次元カラー画像またはＢモード画像だけが
出力され、Ｄ／Ａ変換手段４を介して表示手段５に表示
される。図１０は変化率の２次元カラー画像とＢモード
画像との合成画像の表示例を示している。
【００３２】このように本実施例によれば、Ｂモード画
像と同じ輝度信号を用いて、２次元走査面全域の輝度レ
ベルの時間的変化の空間的な様子が色彩又は彩度の相違
としてＢモード画像と共にリアルタイムで表示される。
したがって、超音波反射特性が良好なコントラスト剤を
血管に注入した場合、コントラスト剤がある部分に流入
してから流出する様子をリアルタイムで且つ空間的に診
断することができる。インデックスをそのレベルに応じ
た輝度信号として画像化しても同様の効果が得られるこ
とは勿論である。
像と同じ輝度信号を用いて、２次元走査面全域の輝度レ
ベルの時間的変化の空間的な様子が色彩又は彩度の相違
としてＢモード画像と共にリアルタイムで表示される。
したがって、超音波反射特性が良好なコントラスト剤を
血管に注入した場合、コントラスト剤がある部分に流入
してから流出する様子をリアルタイムで且つ空間的に診
断することができる。インデックスをそのレベルに応じ
た輝度信号として画像化しても同様の効果が得られるこ
とは勿論である。
【００３３】本発明による関心領域の輝度レベルを得る
とき、操作者が超音波プローブを誤って僅かに動かして
しまった場合や、被検体の呼吸や心拍動による臓器の移
動がある場合、Ｂモード画像又はそれに代わる２次元又
は３次元情報に基づいてパターンマッチング技法等によ
り、関心領域の動きの補正をすることが考えられる。こ
れにより精度向上を図ることができる。
とき、操作者が超音波プローブを誤って僅かに動かして
しまった場合や、被検体の呼吸や心拍動による臓器の移
動がある場合、Ｂモード画像又はそれに代わる２次元又
は３次元情報に基づいてパターンマッチング技法等によ
り、関心領域の動きの補正をすることが考えられる。こ
れにより精度向上を図ることができる。
【００３４】本発明は上述した実施例に限定されず、種
々変形して実施可能である。例えば上述の説明では、輝
度レベルの時間的変化をグラフや数値情報として表示す
ることと説明したが、一般的な超音波ドプラ装置で計算
される周波数シフト（偏移周波数）のパワーを輝度レベ
ルに置換して、このパワーの時間的変化を求めるように
してもよい。また上述の説明では、検波結果としての輝
度信号を用いて時間的変化を求めているが、検波器に入
力される加算器の出力信号を用いて、その振幅値の時間
的変化を求めるようにしてもよい。
々変形して実施可能である。例えば上述の説明では、輝
度レベルの時間的変化をグラフや数値情報として表示す
ることと説明したが、一般的な超音波ドプラ装置で計算
される周波数シフト（偏移周波数）のパワーを輝度レベ
ルに置換して、このパワーの時間的変化を求めるように
してもよい。また上述の説明では、検波結果としての輝
度信号を用いて時間的変化を求めているが、検波器に入
力される加算器の出力信号を用いて、その振幅値の時間
的変化を求めるようにしてもよい。
【００３５】
【発明の効果】本発明によれば、関心領域の血流状態を
良好に診断することができる。
良好に診断することができる。
【００３６】
【図１】第１実施例のブロック図。
【図２】関心領域内全点の平均値の時間的変化を求める
場合の図１の信号レベル時間変化抽出器による２次元走
査面の１回の走査動作に対応する動作手順を示す流れ
図。
場合の図１の信号レベル時間変化抽出器による２次元走
査面の１回の走査動作に対応する動作手順を示す流れ
図。
【図３】図２の場合の表示手段の表示画面を示す図。
【図４】関心領域内の最大値及び／又は最小値の時間的
変化を求める場合の図１の信号レベル時間変化抽出器に
よる２次元走査面の１回の走査動作に対応する動作手順
を示す流れ図。
変化を求める場合の図１の信号レベル時間変化抽出器に
よる２次元走査面の１回の走査動作に対応する動作手順
を示す流れ図。
【図５】図４の場合の表示手段の表示画面を示す図。
【図６】関心領域内全点の平均値の所定の短時間間隔に
おける変化率の時間的変化を求める場合の図１の信号レ
ベル時間変化抽出器による２次元走査面の１回の走査動
作に対応する動作手順を示す流れ図。
おける変化率の時間的変化を求める場合の図１の信号レ
ベル時間変化抽出器による２次元走査面の１回の走査動
作に対応する動作手順を示す流れ図。
【図７】図６の場合の表示手段の表示画面を示す図。
【図８】図１の信号レベル時間変化抽出器で計算される
半値時間幅を説明する図。
半値時間幅を説明する図。
【図９】第２実施例のブロック図。
【図１０】第２実施例の表示画面を示す図。
１…超音波プローブ、２…アナログ部、３…走査変換手
段、４…Ｄ／Ａ変換手段、５…表示手段、６…信号レベ
ル時間変化抽出器、２１…プリアンプ、２２…パルサ、
２４…ディレーライン、２５…加算器、２６…検波器。
段、４…Ｄ／Ａ変換手段、５…表示手段、６…信号レベ
ル時間変化抽出器、２１…プリアンプ、２２…パルサ、
２４…ディレーライン、２５…加算器、２６…検波器。
─────────────────────────────────────────────────────
フロントページの続き
(72)発明者 内海 勲
栃木県大田原市下石上1385番の１ 株式
会社東芝那須工場内
(72)発明者 上木 康至
栃木県大田原市下石上1385番の１ 株式
会社東芝那須工場内
(56)参考文献 特開 平１−254151（ＪＰ，Ａ)
特開 平４−117950（ＪＰ，Ａ)
特開 平４−208143（ＪＰ，Ａ)
特開 平７−59781（ＪＰ，Ａ)
特開 平７−75638（ＪＰ，Ａ)
特開 平７−79974（ＪＰ，Ａ)
米国特許5235984（ＵＳ，Ａ)
Ｆｌｏｒｄｅｌｉｚａ Ｓ．Ｖｉｌｌ
ａｎｕｅｖａ，ｅｔ ａｌ，Ｃｉｒｃｕ
ｌａｔｉｏｎ，1992年 ４月，Ｖｏｌ.
85，Ｎｏ．４，ｐｐ1557−1564
Ｃ．Ｂｏｙｌｅ，ｅｔ ａｌ，Ｃｏｍ
ｐｕｔｅｒｓ ｉｎ Ｃａｒｄｉｏｌｏ
ｇｙ，1988．Ｐｒｏｃｅｅｄｉｎｇｓ,
1989年，ｐｐ203−206
Ｂｙｒｏｎ Ｆ．Ｖａｎｄｅｎｂｅｒ
ｇ，ｅｔ ａｌ，Ｊｏｕｒｎａｌ ｏｆ
ｔｈｅ Ａｍｅｒｉｃａｎ Ｃｏｌｌ
ｅｇｅ ｏｆ Ｃａｒｄｉｏｌｏｇｙ,
1989年 １月，ｖｏｌ．13，ｎｏ．１,
ｐｐ200−206
(58)調査した分野(Int.Cl.7，ＤＢ名)
A61B 8/00
Claims (2)
- 【請求項１】 コントラスト剤を注入した状態で被検体
内を超音波ビームで走査することにより得られた受信信
号の包絡線を検波し、この検波結果を輝度信号として断
層像を表示する超音波診断装置において、 前記断層像の関心領域内の輝度信号の最大値及び最小値
を選択し、前記最大値の時間変化を表すグラフと前記最
小値の時間変化を表すグラフとを生成する手段と、 前記最大値の時間変化を表すグラフと前記最小値の時間
変化を表すグラフとを前記断層像と共に表示する手段と
を具備することを特徴とする超音波診断装置。 - 【請求項２】 コントラスト剤を注入した状態で被検体
内を超音波ビームで走査することにより得られた受信信
号の包絡線を検波し、この検波結果を輝度信号として断
層像を表示する超音波診断装置において、 前記断層像の関心領域内の輝度信号の平均値とこの平均
値の時間変化率とを計算し、前記平均値の時間変化率に
応じた線分長を有する線分を前記平均値と時間に応じて
配列させたグラフを生成する手段と、 前記グラフを前記断層像と共に表示する手段とを具備す
ることを特徴とする超音波診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP32348493A JP3405578B2 (ja) | 1993-12-22 | 1993-12-22 | 超音波診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP32348493A JP3405578B2 (ja) | 1993-12-22 | 1993-12-22 | 超音波診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH07178095A JPH07178095A (ja) | 1995-07-18 |
| JP3405578B2 true JP3405578B2 (ja) | 2003-05-12 |
Family
ID=18155210
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP32348493A Expired - Lifetime JP3405578B2 (ja) | 1993-12-22 | 1993-12-22 | 超音波診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3405578B2 (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4851078B2 (ja) * | 2004-10-14 | 2012-01-11 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波イメージング装置 |
| JP4679141B2 (ja) * | 2004-12-27 | 2011-04-27 | 株式会社東芝 | 超音波診断装置および超音波診断画像の表示方法 |
| JP2009028194A (ja) * | 2007-07-26 | 2009-02-12 | Ge Medical Systems Global Technology Co Llc | 超音波撮像装置 |
| JP6139186B2 (ja) | 2013-03-11 | 2017-05-31 | 東芝メディカルシステムズ株式会社 | 超音波診断装置、画像処理装置及び画像処理プログラム |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH01254151A (ja) * | 1988-04-04 | 1989-10-11 | Fujitsu Ltd | 超音波診断装置 |
| JP2884184B2 (ja) * | 1990-09-07 | 1999-04-19 | 謙 石原 | 超音波診断装置 |
| JP3108879B2 (ja) * | 1990-10-12 | 2000-11-13 | 株式会社日立メディコ | 超音波診断装置 |
| US5235984A (en) * | 1992-03-30 | 1993-08-17 | Hewlett-Packard Company | On-line acoustic densitometry tool for use with an ultrasonic imaging system |
| JP3403809B2 (ja) * | 1993-06-14 | 2003-05-06 | 株式会社東芝 | 超音波診断装置 |
| JP3410821B2 (ja) * | 1993-07-12 | 2003-05-26 | 株式会社東芝 | 超音波診断装置 |
| JPH0779974A (ja) * | 1993-09-13 | 1995-03-28 | Toshiba Corp | 超音波診断装置 |
-
1993
- 1993-12-22 JP JP32348493A patent/JP3405578B2/ja not_active Expired - Lifetime
Non-Patent Citations (3)
| Title |
|---|
| Ｂｙｒｏｎ Ｆ．Ｖａｎｄｅｎｂｅｒｇ，ｅｔ ａｌ，Ｊｏｕｒｎａｌ ｏｆ ｔｈｅ Ａｍｅｒｉｃａｎ Ｃｏｌｌｅｇｅ ｏｆ Ｃａｒｄｉｏｌｏｇｙ，１９８９年 １月，ｖｏｌ．１３，ｎｏ．１，ｐｐ２００−２０６ |
| Ｃ．Ｂｏｙｌｅ，ｅｔ ａｌ，Ｃｏｍｐｕｔｅｒｓ ｉｎ Ｃａｒｄｉｏｌｏｇｙ，１９８８．Ｐｒｏｃｅｅｄｉｎｇｓ，１９８９年，ｐｐ２０３−２０６ |
| Ｆｌｏｒｄｅｌｉｚａ Ｓ．Ｖｉｌｌａｎｕｅｖａ，ｅｔ ａｌ，Ｃｉｒｃｕｌａｔｉｏｎ，１９９２年 ４月，Ｖｏｌ．８５，Ｎｏ．４，ｐｐ１５５７−１５６４ |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH07178095A (ja) | 1995-07-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9241689B2 (en) | Ultrasonic diagnostic equipment and imaging processing apparatus | |
| EP1614387B1 (en) | Ultrasonic diagnostic apparatus, image processing apparatus and image processing method | |
| US5873829A (en) | Diagnostic ultrasound system using harmonic echo imaging | |
| US6638228B1 (en) | Contrast-agent enhanced color-flow imaging | |
| JP3410843B2 (ja) | 超音波診断装置 | |
| CN103648400B (zh) | 超声波诊断装置以及方法 | |
| JP5134787B2 (ja) | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム | |
| JP2007513726A (ja) | 浸透、解像度及びフレームレートの自動制御を有する超音波画像診断システム | |
| JP7313902B2 (ja) | 超音波診断装置 | |
| JP2002085409A (ja) | 超音波診断装置 | |
| JP2007222253A (ja) | 血流画像表示装置 | |
| JPH0924047A (ja) | 超音波診断装置 | |
| US20160220232A1 (en) | Ultrasonic diagnostic apparatus | |
| JPH11327A (ja) | 超音波診断装置 | |
| JP4764209B2 (ja) | 超音波信号解析装置、超音波信号解析方法、超音波解析プログラム、超音波診断装置、及び超音波診断装置の制御方法 | |
| US6322510B1 (en) | Ultrasonic imaging method and apparatus | |
| JPH10118061A (ja) | 超音波診断装置 | |
| JPH10262970A (ja) | 超音波カラードプラ断層装置 | |
| US20120203111A1 (en) | Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image acquisition method | |
| JP3405578B2 (ja) | 超音波診断装置 | |
| JP5606025B2 (ja) | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム | |
| CN100383554C (zh) | 心壁应变成像 | |
| JP4427139B2 (ja) | 超音波診断装置及び超音波送信方法 | |
| JP2008264314A (ja) | 超音波診断装置及び超音波画像データ生成プログラム | |
| US20050075566A1 (en) | Ultrasonice diagnosing apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080307 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090307 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100307 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100307 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110307 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120307 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130307 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140307 Year of fee payment: 11 |
|
| EXPY | Cancellation because of completion of term |