JPH1133024A - 超音波ドプラ診断装置 - Google Patents
超音波ドプラ診断装置Info
- Publication number
- JPH1133024A JPH1133024A JP9194041A JP19404197A JPH1133024A JP H1133024 A JPH1133024 A JP H1133024A JP 9194041 A JP9194041 A JP 9194041A JP 19404197 A JP19404197 A JP 19404197A JP H1133024 A JPH1133024 A JP H1133024A
- Authority
- JP
- Japan
- Prior art keywords
- doppler
- frequency
- ultrasonic
- signal
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
最大周波数や重心周波数のトレース精度を向上し得る超
音波ドプラ診断装置を提供すること。 【解決手段】本発明に係る被検体内の運動流体を含む診
断部位との間で超音波ビームを送受信し、得られた受信
信号から所望のレンジゲート位置の流体に起因したドプ
ラ信号を抽出し、このドプラ信号に基づいてドプラスペ
クトラムを演算し、このドプラスペクトラムの時間変化
を縦軸を周波数、横軸を時間、信号強度又はパワーを輝
度として表すドプラ画像を作成する超音波ドプラ診断装
置は、ドプラ画像の中のノイズが支配的な特定領域の輝
度平均に基づいて閾値を決定し、ドプラ画像に対して輝
度が閾値を越える周波数成分の中の最大周波数を時間ご
とに検出する画像処理演算部301を具備する。
Description
果を利用して、血流などの体内の流体の運動状態を診断
する超音波ドプラ診断装置に関する。
置があるが、その主流は超音波パルス反射法を用いて生
体の軟部組織の断層像を超音波診断装置である。この超
音波診断装置は無侵襲検査法で、組織の断層像を表示す
るものであり、X線診断装置、X線CT装置、MRIお
よび核医学診断装置などの他の診断装置に比べて、リア
ルタイム表示が可能、装置が小型で安価、X線などの被
曝がなく安全性が高く、さらに超音波ドプラ法により血
流イメージングが可能であるなどの独自の特徴を有して
いる。
び産婦人科などでその活用範囲は広い。特に、超音波プ
ローブを体表から割り当てるだけの簡単な操作で心臓の
拍動や胎児の動きの様子がリアルタイム表示で得られ、
かつ安全性が高いため繰り返して検査が行えるほか、ベ
ッドサイドへ移動していっての検査も容易に行えるなど
簡便である。
超音波パルス反射法とを併用し、一つの超音波プローブ
で断層像(白黒Bモード像)と血流情報とを得るととも
に、少なくともその血流情報をリアルタイムで表示でき
るようにした超音波ドプラ診断装置が普及している。
して血流速度を計測するものであり、送信制御及び受信
制御のためのパルス発生回路は、クロック回路、分周回
路などを備えており、所定周波数のクロックパルスに基
づいて超音波繰り返し周波数に相当するレートパルスを
生成し、そのレートパルスをパルサに供給している。パ
ルサは、供給されたレートパルスに同期して高電圧の駆
動電圧パルスで超音波振動子を励振する。この励振によ
り、超音波振動子から超音波パルスが生体内に送波さ
れ、この送波された超音波パルスの一部は、生体内の血
管壁及び血管内の血流(主に赤血球)で反射して超音波
エコーとなる。このエコーは、超音波のドプラ効果を受
けている。つまり、生体内を流れている血流に対して超
音波パルスを走波すると、流動する血流によって散乱さ
れ、ドプラ偏移を受ける。このため、超音波ビームの中
心周波数fcがfdだけ変化し、受信周波数fはf=f
c+fdとなる。このドプラ偏移周波数fdは、血流速
度v、超音波ビームと血管の成す角度θ、音速cとし
て、およそ以下のように表される。
ることにより血流速度vを知ることができる。このため
に、まず前置増幅器でエコー信号を増幅し、その増幅信
号をミキサでクロックパルスと混合し、その混合信号を
次段のローパスフィルタで高周波成分(2・fc)を除
去し、ドプラ偏移周波数fdを中心とする低周波分のみ
を抽出する。この信号をここではドプラ偏移信号と称す
る。
回路からのサンプリングパルスに従ってドプラ偏移信号
にゲートを掛けて、所望の観測位置、すなわちサンプリ
ング・ラスタ上に設定されたレンジゲート(サンプリン
グポイント、サンプリングボリュームともいう)の位置
に対応するドプラ偏移信号をサンプリングする。レンジ
ゲート回路は遅延時間を任意に設定できる回路で、超音
波パルスが振動子とレンジゲート位置との間を超音波が
往復伝搬するに等しい時間だけレートパルスよりも遅延
させ、且つ、設定されたパルス幅でサンプリングパルス
を形成する。なお、レンジゲート位置は、オペレータに
より、Bモード断層上の血流速度を得たい血管の位置
に、トラックボールやジョイスティックで任意に設定さ
れる。
ド回路のサンプリングで生じた高周波成分や血管などの
固定反射信号及び比較的遅い生体内の動きに拠る比較的
低周波のクラッタ成分を除去し、血流のドプラ偏移周波
数成分のみを抽出する。この抽出信号を周波数分析器で
高速フーリエ変換などの周波数分析にかけて、周波数ス
ペクトルパターン(ドプラスペクトラム)を求める。
得、これらを信号強度又はパワーで輝度を割り当てて、
縦軸を周波数として、横軸の時間軸に沿って配列するこ
とにより、図12に示すようなドプラ偏移周波数の時間
変化を作成し、これをリアルタイムで表示している。な
お、この時間変化を表す画像を、ここでは“ドプラ画
像”と称するものとする。
て、最高血流速度の時間変化を、図13中の太線で図示
した如く、最大周波数(最高周波数)のトレースにより
見ることがある。
信号強度がかならずしも十分とれない場合もある。特
に、最大周波数を示す血流信号の強度又はパワーの値は
小さく、ノイズとの分別が重要である。通常ノイズと信
号の分別は信号強度又はパワーに対する閾値処理によっ
て行われる。閾値を決定する方法としては、ある位置
(時刻)の血流成分の最大輝度値をもとに算出した値を
基準値として用いるトレース方法(特開平7−3036
41号公報)、また、ドプラ波形の存在確率の低い領域
の輝度分布から基準ノイズレベルを検出してトレースす
る方法(特開平7−241291号公報)等が提案され
ている。
の太線で図示した如く、重心周波数のトレースにより見
ることがある。さらに、これら最高速度や重心速度のト
レース波形を使って、指定された区間(時間)内におけ
る平均流速値、最高流速値、最低流速値などを計測した
り、これらからさらに様々な特徴量を算出して診断を支
援するようなことも行われている。
ような問題がある。まず、最高周波数のトレースに関し
て、信号(血流)とノイズとを弁別するための閾値を、
血流信号の最高輝度値から予め設定したdBだけ下がっ
た所に求めているが、ドプラ信号のS/Nによっては閾
値が適切とならない場合がある。例えば、S/Nの良い
場合にはドプラ波形の内側(低流速側)をトレースして
しまうし、S/Nの悪い場合にはノイズに引っかかって
局所的に外側(高流速側)をトレースしてしまう。この
様子を図15に示す。
合、算出する領域に血流成分の信号が含まれた場合に
は、正確なトレースができない。特開平7−24129
1号公報においては、血流信号の存在確率の低い所を用
いるとしているが、血流信号が領域に入ることを防ぐの
は難しく、特に拍動性に富む波形では血流信号が当該領
域に入り込むことが容易に予想される。臨床的にもこの
ような波形は極めて一般的である。
ドプラスペクトラム上での重心周波数を自動的にトレー
スする場合に、超音波画像に含まれるノイズの影響によ
りS/Nが低いドプラ画像では正確に重心の位置を求め
ることは難しく、図16に示すように、特に低周波数
(低速)側でエラーが生じやすい。
々な特徴量を算出して診断を行う場合、この算出のため
のパラメータとして必要なのは、平均流速値、最高流速
値、最低流速値、さらにこの他に収縮期最大流速値、拡
張末期流速値が必要となる場合があるが、収縮期最大流
速値、拡張末期流速値に関しては従来では自動検出でき
ない。
上述のようなパラメータを検出することができないた
め、特徴量も心周期に対応させることができない。本発
明は、このような従来技術の問題に鑑みてなされたもの
で、超音波画像に含まれるノイズの影響を抑制して最大
周波数や重心周波数のトレース精度を向上し得る超音波
ドプラ診断装置を提供することを目的としている。
め、この発明に係る超音波ドプラ診断装置は、被検体内
の運動流体を含む診断部位との間で超音波ビームを送受
信する送受信手段と、この送受信手段により得られた受
信信号から所望のレンジゲート位置の流体に起因したド
プラ信号を抽出する抽出手段と、この抽出手段により抽
出されたドプラ信号に基づいてドプラスペクトラムを演
算するスペクトラム演算手段と、このドプラスペクトラ
ムの時間変化を縦軸を周波数、横軸を時間、信号強度又
はパワーを輝度として表すドプラ画像を前記ドプラスペ
クトラムから作成する手段と、前記ドプラ画像の中のノ
イズが支配的な特定領域の輝度平均に基づいて閾値を決
定する手段と、前記ドプラ画像に対して輝度が前記閾値
を越える周波数成分の中の最大周波数を時間ごとに検出
する手段と、前記最大周波数の時間変化を表示する手段
とを具備したことを特徴とする。
は、被検体内の運動流体を含む診断部位との間で超音波
ビームを送受信する送受信手段と、この送受信手段によ
り得られた受信信号から所望のレンジゲート位置の流体
に起因したドプラ信号を抽出する抽出手段と、この抽出
手段により抽出されたドプラ信号に基づいてドプラスペ
クトラムを演算するスペクトラム演算手段と、このドプ
ラスペクトラムの時間変化を縦軸を周波数、横軸を時
間、信号強度又はパワーを輝度として表すドプラ画像を
前記ドプラスペクトラムから作成する手段と、前記ドプ
ラ画像に対して血流信号存在領域を時間ごとに限定する
手段と、前記血流信号存在領域の範囲内で重心周波数を
時間ごとに計算する手段と、前記重心周波数の時間変化
を表示する手段とを具備したことを特徴とする。 (作用)最大周波数を自動的にトレースする場合、血流
成分の信号が存在する割合が低く、ノイズ成分を多く含
む特定の領域の平均輝度から閾値を決めているので、高
精度にノイズを除去して、最大周波数のトレースを高精
度で行うことができる。
存在領域を限定してから、重心周波数をトレースするの
で、ノイズの影響を抑えて、高精度に重心周波数のトレ
ースを行うことができる。また、インデックスを簡便に
求める求めることで、計測操作全体を簡便にすることが
できる。
好ましい実施形態により説明する。図1にこの実施形態
に係る超音波ドプラ診断装置の構成を示している。この
診断装置は血流情報として血流速度を計測するものであ
る。
と、超音波に内部情報を付与する被検体側との間を媒介
するために、配列された複数の微小圧電素子を先端部分
に有している。この振動子201の形態としては、セク
タ対応、リニア対応、コンベックス対応等の中から任意
に選択される。
分周回路などを有しており、クロック回路で発生した図
2に示す所定周波数のクロックパルス(a)を分周し
て、超音波繰り返し周波数に相当するレートパルス
(b)を生成し、そのレートパルスbをパルサ202に
供給する。
ら供給されたレートパルス(b)に同期して高電圧の駆
動電圧パルスを生成し、その駆動電圧パルスにより超音
波振動子201を励振する。この励振により、超音波振
動子201は超音波パルスを生体(P)内に送波する。
送波された超音波パルスの一部は、生体(P)内の血管
壁(BW)及び血管内の血流(B)の主に赤血球で反射
する。この超音波エコーは超音波振動子201により受
信され、電圧信号としてのエコー信号(d)変換され
る。
ンスの差に関する情報の他に、ドプラ効果による周波数
偏移に関する情報を有している。つまり、生体(P)内
を流れている血流(血球)に超音波パルスが当たると、
その流動する血流によって散乱され、周波数偏移を受け
る。このため、超音波ビームの中心周波数fcがfdだ
け変化し、受信周波数fは、(fc+fd)になる。こ
のドプラ偏移周波数fdは、血流速度をv、超音波ビー
ムと血管の成す角度をθ、音速をcとすると、およそ以
下のように表される。
血流速度のビーム方向成分(v・ cosθ)を知ることが
できる。
コー信号(d)を増幅し、その増幅信号をミキサ204
で参照周波数(fc)と混合し、その混合信号(2・f
c+fd)の高周波成分(2・fc)を次段のローパス
フィルタ205で除去する。
とする低周波分のみからなる信号(以下、“ドプラ信
号”と称する)の一部分が、レンジゲート回路211か
らのサンプリングパルスに従ってサンプルホールド回路
206でサンプルされ、ホールドされる。
意に設定できる回路で、超音波パルスが振動子201
と、オペレータにより任意の深さに設定されたレンジゲ
ート(サンプリングポイント、サンプリングボリューム
ともいう)との間の距離を超音波が往復伝搬するに等し
い時間だけ、レートパルス(b)よりも遅延させたタイ
ミングで、所定のパルス幅のサンプリングパルス(c)
を発生する。このサンプリングパルス(c)により、血
流(B)の速度の観測位置、すなわちサンプリング・ラ
スタ上の血流(B)に対して設定されたレンジゲートの
位置で周波数偏移を受けたドプラ信号を抽出することが
できる。
ホールド回路206のサンプリングで生じた高周波成分
や、血管などの固定反射信号及び比較的遅い生体内の動
きによる比較的低周波のクラッタ成分を除去し、血流
(B)のドプラ偏移周波数のみを抽出する。
に送られ、高速フーリエ変換(FFT)などの周波数分
析によってドプラ偏移周波数の周波数スペクトラムが演
算される。この周波数スペクトラムは、表示器209に
送られ、ここで縦軸を周波数、横軸を時間、信号強度又
はパワーを輝度として図12に示したようなドプラ画像
(ドプラスペクトラムの時間変化)として表示される。
像処理演算部301と、画像処理のためにフリーズした
ドプラ画像を記憶するためのフリーズ用メモリ302
と、操作パネル303とを特徴的に備えている。
有し、操作パネル303を介して入力されたフリーズ指
令信号により、フリーズ用メモリ302からフリーズに
係るドプラ画像のデータを読み込み、そしてこのドプラ
画像に対して最大周波数トレースや重心周波数トレース
を実施し、これらトレース結果を利用して様々な特徴的
な速度を検出して、この検出結果を表示器209に供給
する。
操作可能なトラックボールやキーボードを備えており、
この操作パネル303の操作を介して前述したレンジゲ
ートが所望の位置(深さ)に設定され、またフリーズ指
令等が入力され得る。また、操作パネル303は表示器
209に表示されたドプラスペクトラム上の演算位置
(時間位置)信号を画像処理演算部301に出力する。
なお、レンジゲート位置信号及び演算位置信号に対応し
て所定のマーカーが表示器209に表示される。
ート回路211にも送られ、指定したレンジゲートの位
置(深さ)に対応した走査線上に受信超音波ビームを絞
ることができるようになっている。
て図3を参照して説明する。まず、ステップS11にお
いて、まず、画像処理演算部301は、フリーズが指令
されたときに表示されているドプラ画像のデータを、最
大周波数トレースや重心周波数トレース、またこれらト
レース結果を利用した様々な特徴的な速度の検出の対象
として、フリーズ用メモリ302から読み込む。
とノイズ成分とを弁別するための閾値を決定するために
最適な閾値演算領域を選択する。この選択にあたって
は、図4と図5に示すように、複数、個々では6つの閾
値演算領域の候補(以下、“サンプル領域”と称する)
が予め画像上に規定されている。このサンプル領域の中
から、上述した最適な1つの閾値演算領域を選択するの
は、輝度の平均値を領域ごとに演算し、この平均値の比
較により行うのであり、具体的には平均値が最小のサン
プル領域を閾値演算領域として選出する(図6参照)。
上述したように、ドプラ画像の輝度は各周波数での信号
強度又はパワーによって決められており、一般的には、
ランダムに発生するノイズの信号強度又はパワーは、血
流信号の信号強度又はパワーより小さくなる。従って、
平均値が最小のサンプル領域とは、その領域内で血流信
号成分の占める割合が最小で、しかもノイズ成分の占め
る割合が最大となる傾向にある。勿論、予めサンプル領
域を規定せずに、オペレータが血流信号成分の割合が低
く、理想的にはこの割合がゼロで、しかもノイズ成分の
割合が高いと判断した領域を任意の大きさで、任意の形
状に設定できるようにしてもよい。
S13において、図6に示すように、その領域内の輝度
値の標準偏差値を演算し、この標準偏差値とステップS
12で計算済みの輝度平均値とで、次のように閾値を演
算する。
て、その精度から判断して任意の値に変更可能とするこ
とが好適である。
4において、最大周波数を時間軸に沿ってトレースす
る。このトレースの処理の流れを図7に示している。ス
テップS21では、時間軸上で時刻ごとに各画素の輝度
値(信号強度又はパワー)をステップS13で出力され
た閾値に対して、周波数軸に沿って上方(高周波数側)
から順番に比較していき、図8に示すように、ステップ
S22で最初に輝度値が閾値以上と判定された画素の位
置(周波数)を、ステップS23において最大周波数値
として検出する。また、時間軸上で時刻ごとに各画素の
輝度値(信号強度又はパワー)をステップS13で出力
された閾値に対して、周波数軸に沿って下方(低周波数
側)から順番に比較していき、最初に輝度値が閾値以上
と判定された画素の位置(周波数)を、最小周波数値と
して検出する。
間軸方向に連続して繰り返すことにより、最大周波数と
最小周波数を時間軸に沿ってトレースしていき、最終的
に図13に示したような最大周波数(最高速度)の時間
曲線と最小周波数(最低速度)の時間曲線を得る。
ップS15)について説明する。この重心速度トレース
処理の流れを、図9に示している。この処理でも、最大
周波数トレース処理で求めた閾値を使って、画素の輝度
を上方(高周波数側)から順番に当該閾値と比較し、最
初に出現した閾値以上の輝度値を示す画素の位置を最大
周波数として検出し(ステップS31)、また画素の輝
度を、周波数軸に沿って下方(低周波数側)から順番に
当該閾値と比較し、最初に出現した閾値以上の輝度値を
示す画素の位置を最小周波数として検出する(ステップ
S32)。
0に示すように、血流成分が存在する周波数帯域、つま
り血流信号存在領域が限定され得る。ステップS33で
は、この血流信号存在領域に限定して、重心周波数を演
算する。このような血流信号存在領域の限定と、その範
囲内での重心周波数の演算との処理を、ステップS34
を介して時刻ごとに繰り返し行い、重心周波数を時間軸
に沿ってトレースしていき、最終的に図14に示したよ
うな重心周波数(重心速度)の時間曲線を得る。
分が存在する周波数帯域)を限定せずに、周波数軸全域
を対象として重心周波数を演算したとき、S/Nが悪い
ドプラ画像では、図16に示したように特に流速が低い
領域でのトレース精度が悪化する。
数のトレースの方法を用いて血流信号存在領域を周波数
軸上で限定したが、図3のステップS13で求めた閾値
を用いて、閾値以上の輝度を持つ画素群を抽出し、この
抽出した画素群を対象にして重心周波数を演算してもよ
い。また、周波数軸方向に重み付けを行って重心を演算
してもよい。
スペクトラム上の上記ステップS15で出力された最大
周波数トレースまたはステップS16で出力された重心
周波数トレース結果を用いて、図11に示すように指定
範囲内に含まれる数心拍中の個々の心拍周期内での収縮
期最高流速点(PSV)や拡張末期速度点(EDV)、
また個々の心拍周期内での最低流速点(Vmin )を同定
する。各点が満たす条件を以下に説明する。
する。図11に示すように、最大周波数でトレースした
時間曲線の全領域(ドプラ画像の時間幅)からその中の
最も高い周波数に、所定の係数(例えば0.8)を乗算
することにより、収縮期間を判定するためのPSV判定
閾値とする。PSVは、最大周波数が当該PSV判定閾
値を越える時間範囲以内で、最大周波数の時間曲線が凸
部分の最大値(極大値)の位置を収縮期最高流速点(P
SV)として同定する。
定するのに用いたPSV判定閾値を最大周波数が越える
時間範囲内を対象に、その範囲の最も右(時間軸で過
去)側の点から時間軸に沿って検索していき、最初に表
れる凹部分の頂点(極小値)の点である。
検索していくような方法も可能である。しかし、このよ
うな方法では、PSV点付近の上下変動が大きいような
場合、EDVをPSV付近で誤認識してしまうことが予
想されるが、上述したように拡張期間を限定してから、
時間軸上で過去の側から時間軸に沿って検索していくこ
とにより、EDVを正確に同定することができる。最低
流速点(Vmin )は2つのEDV点間のトレース曲線領
域中で最小速度値点である。
プS17で同定された点の時刻やその血流速度(周波
数)に基づいて、ドプラ波形の形状的特徴をインデック
ス値として数値化することができる。
5で得られた最大周波数トレース結果と、ステップS1
6で得られた重心周波数トレース結果と、ステップS1
7で同定された点と、ステップS18で得られたインデ
ックス値を図1の表示器209に出力し、ドプラ画像中
に表示する。
を用いて構成するものと想定して説明したが、それぞれ
の機能を実現するように設計された回路により当該画像
信号処理部をハードウェアで構成しても良い。その他、
本発明は、上述した実施形態に限定されることなく、種
々変形して実施可能である。
トレースする場合、血流成分の信号が存在する割合が低
く、ノイズ成分を多く含む特定の領域の平均輝度から閾
値を決めているので、高精度にノイズを除去して、最大
周波数のトレースを高精度で行うことができる。
存在領域を限定してから、重心周波数をトレースするの
で、ノイズの影響を抑えて、高精度に重心周波数のトレ
ースを行うことができる。また、インデックスを簡便に
求める求めることで、計測操作全体を簡便にすることが
できる。
を示すブロック図。
図。
ル領域の一例を示す図。
ル領域の一例を示す図。
領域の周波数スペクトラムから閾値を演算する処理の説
明図。
ーを示す図。
の周波数スペクトラムから閾値を使って最大周波数を同
定する処理の説明図。
ーを示す図。
間の周波数スペクトラムから閾値を使って血流信号存在
領域を限定する処理の説明図。
数や重心周波数のトレース結果から同定するパラメータ
点を示す図
図。
を示す図。
を示す図。
Claims (16)
- 【請求項1】 被検体内の運動流体を含む診断部位との
間で超音波ビームを送受信する送受信手段と、 この送受信手段により得られた受信信号から所望のレン
ジゲート位置の流体に起因したドプラ信号を抽出する抽
出手段と、 この抽出手段により抽出されたドプラ信号に基づいてド
プラスペクトラムを演算するスペクトラム演算手段と、 このドプラスペクトラムの時間変化を縦軸を周波数、横
軸を時間、信号強度又はパワーを輝度として表すドプラ
画像を前記ドプラスペクトラムから作成する手段と、 前記ドプラ画像の中のノイズが支配的な特定領域の輝度
平均に基づいて閾値を決定する手段と、 前記ドプラ画像に対して輝度が前記閾値を越える周波数
成分の中の最大周波数を時間ごとに検出する手段と、 前記最大周波数の時間変化を表示する手段とを具備した
ことを特徴とする超音波ドプラ診断装置。 - 【請求項2】 前記閾値決定手段は、前記ドプラ画像内
の複数の領域の中から輝度平均が最小になる領域を前記
特定領域として選択する手段を有することを特徴とする
請求項1記載の超音波ドプラ診断装置。 - 【請求項3】 前記閾値決定手段は、前記特定領域内の
輝度値の標準偏差と前記輝度平均とに基づいて前記閾値
を演算する手段を有することを特徴とする請求項1記載
の超音波ドプラ診断装置。 - 【請求項4】 前記最大周波数の時間変化に基づいて収
縮期間を判定する手段をさらに備えることを特徴とする
請求項1記載の超音波ドプラ診断装置。 - 【請求項5】 前記収縮期間を、最大周波数の時間変化
の中で最も高い周波数に所定の係数を乗算することによ
り、収縮期間を判定するための閾値を求め、最大周波数
が当該閾値を越える時間範囲として判定することを特徴
とする請求項4記載の超音波ドプラ診断装置。 - 【請求項6】 前記判定した収縮期間を検索範囲とし
て、収縮期最高流速点、拡張末期速度点、最低流速点の
少なくとも1つを同定する手段をさらに備えることを特
徴とする請求項4記載の超音波ドプラ診断装置。 - 【請求項7】 前記収縮期最高流速点を、前記検索範囲
で前記最大周波数の時間変化の極大点に同定することを
特徴とする請求項6記載の超音波ドプラ診断装置。 - 【請求項8】 前記拡張末期速度点を、前記検索範囲内
で時間軸に沿って検索していき、最初に出現する極小点
に同定することを特徴とする請求項6記載の超音波ドプ
ラ診断装置。 - 【請求項9】 被検体内の運動流体を含む診断部位との
間で超音波ビームを送受信する送受信手段と、 この送受信手段により得られた受信信号から所望のレン
ジゲート位置の流体に起因したドプラ信号を抽出する抽
出手段と、 この抽出手段により抽出されたドプラ信号に基づいてド
プラスペクトラムを演算するスペクトラム演算手段と、 このドプラスペクトラムの時間変化を縦軸を周波数、横
軸を時間、信号強度又はパワーを輝度として表すドプラ
画像を前記ドプラスペクトラムから作成する手段と、 前記ドプラ画像に対して血流信号存在領域を時間ごとに
限定する手段と、 前記血流信号存在領域の範囲内で重心周波数を時間ごと
に計算する手段と、 前記重心周波数の時間変化を表示する手段とを具備した
ことを特徴とする超音波ドプラ診断装置。 - 【請求項10】 前記血流信号存在領域限定手段は、前
記ドプラ画像の中のノイズが支配的な特定領域の輝度平
均に基づいて閾値を決定する手段と、前記ドプラ画像に
対して輝度が前記閾値を越える周波数成分の中の最小周
波数と最大周波数を時間ごとに検出する手段とを有する
ことを特徴とした請求項9記載の超音波ドプラ診断装
置。 - 【請求項11】 前記閾値決定手段は、前記ドプラ画像
内の複数の領域の中から輝度平均が最小になる領域を前
記特定領域として選択する手段を有することを特徴とす
る請求項10記載の超音波ドプラ診断装置。 - 【請求項12】 前記閾値決定手段は、前記特定領域内
の輝度値の標準偏差と前記輝度平均とに基づいて前記閾
値を演算する手段を有することを特徴とする請求項10
記載の超音波ドプラ診断装置。 - 【請求項13】 前記重心周波数の時間変化に基づいて
収縮期間を判定する手段をさらに備えることを特徴とす
る請求項9記載の超音波ドプラ診断装置。 - 【請求項14】 前記判定した収縮期間を検索範囲とし
て、収縮期最高流速点、拡張末期速度点、最低流速点の
少なくとも1つを同定する手段をさらに備えることを特
徴とする請求項13記載の超音波ドプラ診断装置。 - 【請求項15】 前記収縮期最高流速点を、前記検索範
囲で前記重心周波数の時間変化の極大点に同定すること
を特徴とする請求項14記載の超音波ドプラ診断装置。 - 【請求項16】 前記拡張末期速度点を、前記検索範囲
内で時間軸に沿って検索していき、最初に出現する極小
点に同定することを特徴とする請求項14記載の超音波
ドプラ診断装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP19404197A JP3892538B2 (ja) | 1997-07-18 | 1997-07-18 | 超音波ドプラ診断装置 |
US09/116,963 US6050948A (en) | 1997-07-18 | 1998-07-17 | Ultrasound Doppler diagnostic apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP19404197A JP3892538B2 (ja) | 1997-07-18 | 1997-07-18 | 超音波ドプラ診断装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH1133024A true JPH1133024A (ja) | 1999-02-09 |
JP3892538B2 JP3892538B2 (ja) | 2007-03-14 |
Family
ID=16317964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP19404197A Expired - Lifetime JP3892538B2 (ja) | 1997-07-18 | 1997-07-18 | 超音波ドプラ診断装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US6050948A (ja) |
JP (1) | JP3892538B2 (ja) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000229082A (ja) * | 1999-02-10 | 2000-08-22 | Toshiba Corp | 超音波ドプラ診断装置 |
JP2003245279A (ja) * | 2002-02-19 | 2003-09-02 | Ge Medical Systems Global Technology Co Llc | スペクトル・ドプラ・イメージングの自動制御のための方法及び装置 |
JP2003284718A (ja) * | 2002-03-28 | 2003-10-07 | Toshiba Corp | 超音波診断装置 |
JP2004127064A (ja) * | 2002-10-04 | 2004-04-22 | Konica Minolta Holdings Inc | 画像処理方法、画像処理装置、画像処理プログラム及び画像記録装置 |
WO2005115248A1 (ja) * | 2004-05-26 | 2005-12-08 | Hitachi Medical Corporation | 超音波撮像装置 |
JP2006326291A (ja) * | 2005-05-26 | 2006-12-07 | Medison Co Ltd | 超音波スペクトル映像を処理する装置及び方法 |
JP2007152111A (ja) * | 2005-12-01 | 2007-06-21 | General Electric Co <Ge> | スペクトルドプラ利得を自動調整するための方法及び装置 |
JP2008511367A (ja) * | 2004-08-30 | 2008-04-17 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | ドップラ速度スペクトルにおける流速の調節可能トレース方法及び装置 |
JP2008512167A (ja) * | 2004-09-13 | 2008-04-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 超音波を使用して体液の流動性を測定および/又は検出する方法および装置 |
JP2008514263A (ja) * | 2004-09-28 | 2008-05-08 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 超音波により外部的に測定された体液の流れ挙動に関する情報を提示する方法及び装置 |
WO2008123596A1 (en) * | 2007-03-29 | 2008-10-16 | Aloka Co., Ltd. | Methods and apparatus for ultrasound imaging |
JP2009532089A (ja) * | 2006-03-30 | 2009-09-10 | アロカ株式会社 | 超音波受信ビームフォーマ用の遅延コントローラ |
WO2012160873A1 (ja) * | 2011-05-25 | 2012-11-29 | 日立アロカメディカル株式会社 | 超音波画像処理装置およびプログラム |
JP2013526890A (ja) * | 2009-10-27 | 2013-06-27 | パルティ、ヨーラム | 経胸壁肺ドップラー超音波 |
JP2016025958A (ja) * | 2015-10-07 | 2016-02-12 | パルティ、ヨーラム | 経胸壁肺ドップラー超音波 |
US20160058410A1 (en) * | 2014-09-01 | 2016-03-03 | Samsung Medison Co., Ltd. | Ultrasound diagnosis apparatus and method of operating the same |
KR20170075363A (ko) | 2015-12-23 | 2017-07-03 | 지멘스 메디컬 솔루션즈 유에스에이, 인크. | 도플러 스펙트럼 영상을 표시하는 초음파 시스템 및 방법 |
JP2018157929A (ja) * | 2017-03-22 | 2018-10-11 | 株式会社日立製作所 | 超音波診断装置及び波形処理方法 |
JP2020156810A (ja) * | 2019-03-27 | 2020-10-01 | 株式会社日立製作所 | 超音波診断装置、トレース方法及びプログラム |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6221020B1 (en) * | 1999-04-22 | 2001-04-24 | G.E. Medical Systems Global Technology Company, Llc | System and method for providing variable ultrasound analyses in a post-storage mode |
KR100381875B1 (ko) * | 2000-07-08 | 2003-04-26 | 주식회사 메디슨 | 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법 |
JP4245428B2 (ja) * | 2003-07-29 | 2009-03-25 | 東芝医用システムエンジニアリング株式会社 | 超音波ドプラ診断装置及び画像データ生成方法 |
DE10338940B3 (de) * | 2003-08-22 | 2005-02-10 | Fresenius Medical Care Deutschland Gmbh | Verfahren zum Messen einer Signallaufzeit in Blut und Vorrichtung zur Anwendung des Verfahrens |
CN100443056C (zh) * | 2003-10-17 | 2008-12-17 | 松下电器产业株式会社 | 超声波多普勒血流测定装置 |
US7248208B2 (en) * | 2004-06-10 | 2007-07-24 | Honeywell International Inc. | Methods and systems for maintaining a position during hovering operations |
US8303507B2 (en) * | 2004-09-07 | 2012-11-06 | Kabushiki Kaisha Toshiba | Ultrasonic doppler diagnostic apparatus and measuring method of diagnostic parameter |
US7513872B2 (en) * | 2004-10-18 | 2009-04-07 | Kabushiki Kaisha Toshiba | Ultrasonic doppler measuring apparatus and control method therefor |
WO2006098354A1 (ja) * | 2005-03-15 | 2006-09-21 | Kabushiki Kaisha Toshiba | 超音波診断装置及びその制御方法 |
CN100525713C (zh) * | 2005-06-10 | 2009-08-12 | 深圳迈瑞生物医疗电子股份有限公司 | 声谱图包络曲线提取方法 |
KR20080091350A (ko) * | 2006-01-27 | 2008-10-10 | 코닌클리케 필립스 일렉트로닉스 엔.브이. | 자동 초음파 도플러 측정 |
JP5300188B2 (ja) * | 2006-09-11 | 2013-09-25 | 株式会社東芝 | 超音波診断装置及び超音波診断装置の制御プログラム |
CN101301212B (zh) * | 2007-05-11 | 2011-08-17 | 深圳迈瑞生物医疗电子股份有限公司 | 实时估计多普勒参数的方法及装置 |
JP4960162B2 (ja) * | 2007-07-17 | 2012-06-27 | 株式会社東芝 | 超音波診断装置 |
KR100951586B1 (ko) | 2008-10-16 | 2010-04-09 | 주식회사 메디슨 | 도플러 스펙트럼 영상을 처리하는 방법 |
CA3103546A1 (en) * | 2011-02-03 | 2012-08-09 | Yoram Palti | Transthoracic cardio-pulmonary monitor |
US9783929B2 (en) | 2011-04-28 | 2017-10-10 | Abb Schweiz Ag | Determination of CD and/or MD variations from scanning measurements of a sheet of material |
CN102499718B (zh) * | 2011-11-25 | 2013-04-03 | 深圳市理邦精密仪器股份有限公司 | 一种声谱图最大频率提取的方法和装置 |
EP3323350A3 (en) * | 2016-11-16 | 2018-06-20 | Huntleigh Technology Limited | Ultrasound doppler monitoring system |
CN106943158A (zh) * | 2017-03-27 | 2017-07-14 | 深圳安盛生物医疗技术有限公司 | 一种频谱参数实时计算方法 |
WO2019096031A1 (en) * | 2017-11-14 | 2019-05-23 | Edan Instruments, Inc. | Systems and methods for multi-resolution discriminant analysis for ultrasound imaging |
US11207055B2 (en) | 2018-10-08 | 2021-12-28 | General Electric Company | Ultrasound Cardiac Doppler study automation |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4373532A (en) * | 1980-07-07 | 1983-02-15 | Palo Alto Medical Research Foundation | Ultrasonic marker for physiologic diagnosis and method of using same |
US5146414A (en) * | 1990-04-18 | 1992-09-08 | Interflo Medical, Inc. | Method and apparatus for continuously measuring volumetric flow |
US5287753A (en) * | 1992-05-02 | 1994-02-22 | Advanced Technology Laboratories, Inc. | Continuous display of peak and mean blood flow velocities |
JPH06327672A (ja) * | 1993-05-18 | 1994-11-29 | Toshiba Medical Eng Co Ltd | 超音波ドプラ診断装置 |
JP2664633B2 (ja) * | 1994-03-07 | 1997-10-15 | アロカ株式会社 | 超音波ドプラ診断装置 |
US5634465A (en) * | 1995-06-09 | 1997-06-03 | Advanced Technology Laboratories, Inc. | Continuous display of cardiac blood flow information |
US5647366A (en) * | 1996-09-17 | 1997-07-15 | Siemens Medical Systems, Inc. | Method and system for automatic measurements of doppler waveforms |
US5935074A (en) * | 1997-10-06 | 1999-08-10 | General Electric Company | Method and apparatus for automatic tracing of Doppler time-velocity waveform envelope |
-
1997
- 1997-07-18 JP JP19404197A patent/JP3892538B2/ja not_active Expired - Lifetime
-
1998
- 1998-07-17 US US09/116,963 patent/US6050948A/en not_active Expired - Lifetime
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000229082A (ja) * | 1999-02-10 | 2000-08-22 | Toshiba Corp | 超音波ドプラ診断装置 |
JP2003245279A (ja) * | 2002-02-19 | 2003-09-02 | Ge Medical Systems Global Technology Co Llc | スペクトル・ドプラ・イメージングの自動制御のための方法及び装置 |
JP2003284718A (ja) * | 2002-03-28 | 2003-10-07 | Toshiba Corp | 超音波診断装置 |
JP2004127064A (ja) * | 2002-10-04 | 2004-04-22 | Konica Minolta Holdings Inc | 画像処理方法、画像処理装置、画像処理プログラム及び画像記録装置 |
JP4996247B2 (ja) * | 2004-05-26 | 2012-08-08 | 株式会社日立メディコ | 超音波診断装置 |
WO2005115248A1 (ja) * | 2004-05-26 | 2005-12-08 | Hitachi Medical Corporation | 超音波撮像装置 |
US8353836B2 (en) | 2004-05-26 | 2013-01-15 | Hitachi Medical Corporation | Ultrasonographic device for performing luminance correction |
JPWO2005115248A1 (ja) * | 2004-05-26 | 2008-03-27 | 株式会社日立メディコ | 超音波撮像装置 |
JP2008511367A (ja) * | 2004-08-30 | 2008-04-17 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | ドップラ速度スペクトルにおける流速の調節可能トレース方法及び装置 |
JP2008512167A (ja) * | 2004-09-13 | 2008-04-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 超音波を使用して体液の流動性を測定および/又は検出する方法および装置 |
JP2008514263A (ja) * | 2004-09-28 | 2008-05-08 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 超音波により外部的に測定された体液の流れ挙動に関する情報を提示する方法及び装置 |
JP2006326291A (ja) * | 2005-05-26 | 2006-12-07 | Medison Co Ltd | 超音波スペクトル映像を処理する装置及び方法 |
JP2007152111A (ja) * | 2005-12-01 | 2007-06-21 | General Electric Co <Ge> | スペクトルドプラ利得を自動調整するための方法及び装置 |
JP2009532089A (ja) * | 2006-03-30 | 2009-09-10 | アロカ株式会社 | 超音波受信ビームフォーマ用の遅延コントローラ |
WO2008123596A1 (en) * | 2007-03-29 | 2008-10-16 | Aloka Co., Ltd. | Methods and apparatus for ultrasound imaging |
JP2010522581A (ja) * | 2007-03-29 | 2010-07-08 | アロカ株式会社 | 超音波画像処理方法および超音波画像処理装置 |
JP2013526890A (ja) * | 2009-10-27 | 2013-06-27 | パルティ、ヨーラム | 経胸壁肺ドップラー超音波 |
US9817119B2 (en) | 2011-05-25 | 2017-11-14 | Hitachi, Ltd. | Ultrasound image processing apparatus and program |
JP2012245049A (ja) * | 2011-05-25 | 2012-12-13 | Hitachi Aloka Medical Ltd | 超音波画像処理装置およびプログラム |
CN103561661A (zh) * | 2011-05-25 | 2014-02-05 | 日立阿洛卡医疗株式会社 | 超声图像处理装置和程序 |
WO2012160873A1 (ja) * | 2011-05-25 | 2012-11-29 | 日立アロカメディカル株式会社 | 超音波画像処理装置およびプログラム |
US20160058410A1 (en) * | 2014-09-01 | 2016-03-03 | Samsung Medison Co., Ltd. | Ultrasound diagnosis apparatus and method of operating the same |
US10743835B2 (en) * | 2014-09-01 | 2020-08-18 | Samsung Medison Co., Ltd. | Ultrasound diagnosis apparatus and method of operating the same |
JP2016025958A (ja) * | 2015-10-07 | 2016-02-12 | パルティ、ヨーラム | 経胸壁肺ドップラー超音波 |
KR20170075363A (ko) | 2015-12-23 | 2017-07-03 | 지멘스 메디컬 솔루션즈 유에스에이, 인크. | 도플러 스펙트럼 영상을 표시하는 초음파 시스템 및 방법 |
JP2018157929A (ja) * | 2017-03-22 | 2018-10-11 | 株式会社日立製作所 | 超音波診断装置及び波形処理方法 |
JP2020156810A (ja) * | 2019-03-27 | 2020-10-01 | 株式会社日立製作所 | 超音波診断装置、トレース方法及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
JP3892538B2 (ja) | 2007-03-14 |
US6050948A (en) | 2000-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3892538B2 (ja) | 超音波ドプラ診断装置 | |
Schäberle | Ultrasonography in vascular diagnosis: a therapy-oriented textbook and atlas | |
US11635514B2 (en) | Imaging methods and apparatuses for performing shear wave elastography imaging | |
US7749165B2 (en) | Instantaneous ultrasonic echo measurement of bladder volume with a limited number of ultrasound beams | |
JP3746115B2 (ja) | 超音波診断装置 | |
US6884216B2 (en) | Ultrasound diagnosis apparatus and ultrasound image display method and apparatus | |
US8597191B2 (en) | Ultrasonic imaging apparatus and a method of generating ultrasonic images | |
EP2437666B1 (en) | Method and apparatus for measuring heart contractility | |
US7871379B2 (en) | Ultrasonic diagnostic apparatus and method of ultrasonic measurement | |
KR100748178B1 (ko) | 임의 m-모드 영상을 디스플레이하는 초음파 진단 시스템및 방법 | |
CN110415248B (zh) | 一种基于超声的血管监测方法、装置、设备及存储介质 | |
JPH08308842A (ja) | 超音波診断装置 | |
JP2003010183A (ja) | 超音波診断装置 | |
JP2008528167A (ja) | 非3d用途における心電図を使用せずに心拍数を導出する方法及びシステム | |
US20060079783A1 (en) | Method and system for deriving a fetal heart rate without the use of an electrocardiogram in non-3D imaging applications | |
Atkinson et al. | Pulse-Doppler ultrasound and its clinical application | |
JP3668687B2 (ja) | 脈波伝播速度計測装置及び超音波診断装置 | |
JP3238467B2 (ja) | 超音波ドップラ診断装置 | |
JP4679141B2 (ja) | 超音波診断装置および超音波診断画像の表示方法 | |
JP4373718B2 (ja) | 血管計測機能を有する超音波診断装置 | |
Bonnefous | Blood flow and tissue motion with ultrasound for vascular applications | |
JP3349695B2 (ja) | 超音波診断装置 | |
JPWO2013161228A1 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
JP2004173753A (ja) | 超音波診断装置 | |
Merritt | Spectral analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050906 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050913 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20051109 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20060822 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061017 |
|
A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20061030 |
|
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: 20061205 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20061207 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091215 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101215 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101215 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111215 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121215 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131215 Year of fee payment: 7 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313114 Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 Free format text: JAPANESE INTERMEDIATE CODE: R313114 Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
EXPY | Cancellation because of completion of term |