JP3042964B2 - 超音波ドプラ診断装置 - Google Patents
超音波ドプラ診断装置Info
- Publication number
- JP3042964B2 JP3042964B2 JP6148885A JP14888594A JP3042964B2 JP 3042964 B2 JP3042964 B2 JP 3042964B2 JP 6148885 A JP6148885 A JP 6148885A JP 14888594 A JP14888594 A JP 14888594A JP 3042964 B2 JP3042964 B2 JP 3042964B2
- Authority
- JP
- Japan
- Prior art keywords
- clutter
- doppler
- frequency
- signal
- unit
- 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 - Fee Related
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、超音波ドプラ診断装
置、特に受信信号に含まれるクラッタ信号の除去に関す
る。
置、特に受信信号に含まれるクラッタ信号の除去に関す
る。
【０００２】
【従来の技術】超音波ドプラ診断装置は、ドプラ偏移周
波数に基づき血流等の生体内運動体の運動情報（速度
等）を画像表示する装置である。例えば、心臓のドプラ
画像を形成する場合、本来観察したい血流のほかに、動
きのある心臓壁等の目障りなクラッタが表示されてしま
う。心臓壁等の低速運動体からの不要反射波は、血流に
比べ非常に強く、このため受信信号中でクラッタ信号の
レベルは大きい。
波数に基づき血流等の生体内運動体の運動情報（速度
等）を画像表示する装置である。例えば、心臓のドプラ
画像を形成する場合、本来観察したい血流のほかに、動
きのある心臓壁等の目障りなクラッタが表示されてしま
う。心臓壁等の低速運動体からの不要反射波は、血流に
比べ非常に強く、このため受信信号中でクラッタ信号の
レベルは大きい。
【０００３】そこで、従来の超音波ドプラ診断装置で
は、いわゆるウォールフィルタと呼ばれるクラッタ除去
フィルタが設けられている。このウォールフィルタは、
ローカットフィルタであり、低速運動体のドプラ成分を
低減するものである。ここで図４には、ウォールフィル
タの遮断特性が示されている。
は、いわゆるウォールフィルタと呼ばれるクラッタ除去
フィルタが設けられている。このウォールフィルタは、
ローカットフィルタであり、低速運動体のドプラ成分を
低減するものである。ここで図４には、ウォールフィル
タの遮断特性が示されている。
【０００４】
【発明が解決しようとする課題】従って、いま図４の低
い周波数のクラッタａが入力した場合、ウォールフィル
タによってクラッタのレベルはａ´のように低減され
る。
い周波数のクラッタａが入力した場合、ウォールフィル
タによってクラッタのレベルはａ´のように低減され
る。
【０００５】しかし、クラッタｂのような高い周波数の
クラッタが入力した場合、ウォールフィルタでは低減で
きず、そのまま通過してしまう。そして、これが画像の
劣化を引き起こす。
クラッタが入力した場合、ウォールフィルタでは低減で
きず、そのまま通過してしまう。そして、これが画像の
劣化を引き起こす。
【０００６】なお、ウォールフィルタの遮断特性は、ク
ラッタ信号及び血流信号のドプラ周波数に応じて選択で
きるが、カットオフ周波数を上げすぎると、血流信号が
損なわれる。また、血流信号への影響をできる限り少な
くするにはフィルタ特性を急峻にする必要があるが、超
音波ドプラ診断装置のような場合には、処理時間の制約
からそれには限界がある。
ラッタ信号及び血流信号のドプラ周波数に応じて選択で
きるが、カットオフ周波数を上げすぎると、血流信号が
損なわれる。また、血流信号への影響をできる限り少な
くするにはフィルタ特性を急峻にする必要があるが、超
音波ドプラ診断装置のような場合には、処理時間の制約
からそれには限界がある。
【０００７】そこで、本発明は、効果的にクラッタの除
去を行うことを目的とする。
去を行うことを目的とする。
【０００８】また、本発明は、クラッタの性質に応じて
適切なクラッタ除去を行うことを目的とする。
適切なクラッタ除去を行うことを目的とする。
【０００９】
【課題を解決するための手段】上記目的を達成するため
に、請求項１記載の発明は、生体内運動体からの受信信
号からドプラ偏移周波数を演算するドプラ演算部と、前
記ドプラ演算部への入力前に分岐取り出しされた受信信
号に含まれるクラッタ信号のドプラ偏移周波数であるク
ラッタ周波数を演算するクラッタ周波数演算部と、ドプ
ラ偏移周波数軸上に前記クラッタ周波数に応じてクラッ
タ判別範囲を適応的に設定し、前記受信信号のドプラ偏
移周波数が前記クラッタ判別範囲に入る場合に、クラッ
タ信号の存在を判別するクラッタ判別部と、前記クラッ
タ信号の存在が判別された場合に、前記ドプラ演算部出
力に対してクラッタ信号の除去を行うクラッタ除去部
と、を含み、前記ドプラ演算部に対し、前記クラッタ周
波数演算部及び前記クラッタ判別部が並列に設けられた
ことを特徴とする。
に、請求項１記載の発明は、生体内運動体からの受信信
号からドプラ偏移周波数を演算するドプラ演算部と、前
記ドプラ演算部への入力前に分岐取り出しされた受信信
号に含まれるクラッタ信号のドプラ偏移周波数であるク
ラッタ周波数を演算するクラッタ周波数演算部と、ドプ
ラ偏移周波数軸上に前記クラッタ周波数に応じてクラッ
タ判別範囲を適応的に設定し、前記受信信号のドプラ偏
移周波数が前記クラッタ判別範囲に入る場合に、クラッ
タ信号の存在を判別するクラッタ判別部と、前記クラッ
タ信号の存在が判別された場合に、前記ドプラ演算部出
力に対してクラッタ信号の除去を行うクラッタ除去部
と、を含み、前記ドプラ演算部に対し、前記クラッタ周
波数演算部及び前記クラッタ判別部が並列に設けられた
ことを特徴とする。
【００１０】請求項２記載の発明は、生体内運動体から
の受信信号からドプラ偏移周波数を演算するドプラ演算
部と、前記ドプラ演算部への入力前に分岐取り出しされ
た受信信号に含まれるクラッタ信号のドプラ偏移周波数
であるクラッタ周波数を演算するクラッタ周波数演算部
と、前記受信信号のドプラ偏移成分のパワーを演算する
パワー演算部と、ドプラ偏移周波数軸とパワー軸とで構
成される二次元空間内に、前記クラッタ周波数に応じて
クラッタ判別範囲を適応的に設定し、前記受信信号のド
プラ偏移周波数及び前記パワーによる空間座標が前記ク
ラッタ判別範囲内に入る場合に、クラッタ信号の存在を
判別するクラッタ判別部と、前記クラッタ信号の存在が
判別された場合に、前記ドプラ演算部出力に対してクラ
ッタ信号の除去を行うクラッタ除去部と、を含み、前記
ドプラ演算部に対し、前記クラッタ周波数演算部、前記
パワー演算部及び前記クラッタ判別部が並列に設けられ
たことを特徴とする。
の受信信号からドプラ偏移周波数を演算するドプラ演算
部と、前記ドプラ演算部への入力前に分岐取り出しされ
た受信信号に含まれるクラッタ信号のドプラ偏移周波数
であるクラッタ周波数を演算するクラッタ周波数演算部
と、前記受信信号のドプラ偏移成分のパワーを演算する
パワー演算部と、ドプラ偏移周波数軸とパワー軸とで構
成される二次元空間内に、前記クラッタ周波数に応じて
クラッタ判別範囲を適応的に設定し、前記受信信号のド
プラ偏移周波数及び前記パワーによる空間座標が前記ク
ラッタ判別範囲内に入る場合に、クラッタ信号の存在を
判別するクラッタ判別部と、前記クラッタ信号の存在が
判別された場合に、前記ドプラ演算部出力に対してクラ
ッタ信号の除去を行うクラッタ除去部と、を含み、前記
ドプラ演算部に対し、前記クラッタ周波数演算部、前記
パワー演算部及び前記クラッタ判別部が並列に設けられ
たことを特徴とする。
【００１１】請求項３記載の発明は、前記受信信号中の
クラッタ信号の分散を演算するクラッタ分散演算器を含
み、前記クラッタ判別範囲は、更に前記分散に対して適
応的に設定され、前記ドプラ演算部に対し、前記クラッ
タ分散演算器も並列に設けられたことを特徴とする。ま
た、請求項４及び請求項５の発明は、低域通過フィルタ
を通過する前の受信信号からクラッタ周波数を演算する
ことを特徴とする。
クラッタ信号の分散を演算するクラッタ分散演算器を含
み、前記クラッタ判別範囲は、更に前記分散に対して適
応的に設定され、前記ドプラ演算部に対し、前記クラッ
タ分散演算器も並列に設けられたことを特徴とする。ま
た、請求項４及び請求項５の発明は、低域通過フィルタ
を通過する前の受信信号からクラッタ周波数を演算する
ことを特徴とする。
【００１２】
【作用】上記構成によれば、受信信号に含まれるクラッ
タ周波数（例えば、クラッタスペクトルの平均周波数）
が演算され、そのクラッタ周波数に応じてクラッタ判別
範囲が設定される。そして、そのクラッタ判別範囲内
に、受信信号のドプラ偏移周波数が入った場合、クラッ
タの存在が判別され、ドプラ演算部出力に対してクラッ
タ除去処理が実行される。
タ周波数（例えば、クラッタスペクトルの平均周波数）
が演算され、そのクラッタ周波数に応じてクラッタ判別
範囲が設定される。そして、そのクラッタ判別範囲内
に、受信信号のドプラ偏移周波数が入った場合、クラッ
タの存在が判別され、ドプラ演算部出力に対してクラッ
タ除去処理が実行される。
【００１３】ここで、クラッタ判別範囲は、ドプラ偏移
周波数上に設定できるほか、ドプラ偏移周波数軸とパワ
ー軸とで構成される二次元空間内に設定できる。二次元
空間内に設定すれば、ドプラ偏移成分のパワーを考慮し
て、より精度良く判別を行うことができる。更に、分散
に応じてクラッタ判別範囲を変化させて最適なクラッタ
除去処理を実現できる。
周波数上に設定できるほか、ドプラ偏移周波数軸とパワ
ー軸とで構成される二次元空間内に設定できる。二次元
空間内に設定すれば、ドプラ偏移成分のパワーを考慮し
て、より精度良く判別を行うことができる。更に、分散
に応じてクラッタ判別範囲を変化させて最適なクラッタ
除去処理を実現できる。
【００１４】
【実施例】以下、本発明の好適な実施例を図面に基づい
て説明する。
て説明する。
【００１５】図１には、本発明に係る超音波ドプラ診断
装置の好適な実施例が示されており、図１はその全体構
成図である。図１において、送受信器１２の駆動の下、
超音波探触子１０によって生体１４に対して超音波パル
スが放射され、生体内の各運動反射体にてドプラシフト
を受けた反射波が超音波探触子１０にて受波される。こ
れにより得られる受信信号は、送受信器１２を介して直
交検波部１６に送られる。直交検波部１６は、参照信号
を用いて受信信号に対して直交検波を行い、直交検波後
の信号１００を出力する。
装置の好適な実施例が示されており、図１はその全体構
成図である。図１において、送受信器１２の駆動の下、
超音波探触子１０によって生体１４に対して超音波パル
スが放射され、生体内の各運動反射体にてドプラシフト
を受けた反射波が超音波探触子１０にて受波される。こ
れにより得られる受信信号は、送受信器１２を介して直
交検波部１６に送られる。直交検波部１６は、参照信号
を用いて受信信号に対して直交検波を行い、直交検波後
の信号１００を出力する。
【００１６】この信号１００は、ドプラ演算部１８とク
ラッタ周波数演算部２０に送られている。つまり、ドプ
ラ演算部１８の手前から分岐取り出しされた信号がクラ
ッタ周波数演算部２０に入力されている。
ラッタ周波数演算部２０に送られている。つまり、ドプ
ラ演算部１８の手前から分岐取り出しされた信号がクラ
ッタ周波数演算部２０に入力されている。
【００１７】ドプラ演算部１８において、信号１００
は、まずクラッタ信号を除去するウォールフィルタ２２
に入力され、ここで図４に示した遮断特性に従いクラッ
タ信号が低減される。そして、その出力信号は周知の自
己相関器２４に入力され、自己相関演算が実行され、自
己相関信号１０２が得られる。
は、まずクラッタ信号を除去するウォールフィルタ２２
に入力され、ここで図４に示した遮断特性に従いクラッ
タ信号が低減される。そして、その出力信号は周知の自
己相関器２４に入力され、自己相関演算が実行され、自
己相関信号１０２が得られる。
【００１８】速度演算器２６は、前記自己相関信号１０
２の実数部と虚数部の逆正接演算によって速度を求める
ものであり、ドプラ偏移周波数ｆｄを表すドプラ演算部
出力信号１０４が出力されている。そして、この信号１
０４はクラッタ除去回路２８を介して表示器３０に送ら
れており、この表示器３０には二次元ドプラ画像が表示
される。
２の実数部と虚数部の逆正接演算によって速度を求める
ものであり、ドプラ偏移周波数ｆｄを表すドプラ演算部
出力信号１０４が出力されている。そして、この信号１
０４はクラッタ除去回路２８を介して表示器３０に送ら
れており、この表示器３０には二次元ドプラ画像が表示
される。
【００１９】ドプラ演算部１８に設けられたパワー演算
器３２は、自己相関器２４の出力信号１０２に基づいて
受信信号中のドプラ偏移成分のパワーを演算するもので
あり、演算されたパワー｜Ｒ（Ｔ）｜が出力されてい
る。
器３２は、自己相関器２４の出力信号１０２に基づいて
受信信号中のドプラ偏移成分のパワーを演算するもので
あり、演算されたパワー｜Ｒ（Ｔ）｜が出力されてい
る。
【００２０】前記クラッタ周波数演算部２０において、
直交検波部１６の出力信号１００は、自己相関器３４に
入力され、自己相関演算が実行されている。この自己相
関器３４は、上述した自己相関器２４と同一の構成を有
するものであるが、その前段にウォールフィルタ２２が
設けられていないので、自己相関器３４の出力はほぼク
ラッタ信号に対する自己相関演算結果を示すことにな
る。すなわち、血流信号に対してクラッタ信号は数十倍
〜数百倍程度大きく、ウォールフィルタ２２を介さずに
自己相関演算を行えば、血流信号の寄与分を事実上無視
できる。
直交検波部１６の出力信号１００は、自己相関器３４に
入力され、自己相関演算が実行されている。この自己相
関器３４は、上述した自己相関器２４と同一の構成を有
するものであるが、その前段にウォールフィルタ２２が
設けられていないので、自己相関器３４の出力はほぼク
ラッタ信号に対する自己相関演算結果を示すことにな
る。すなわち、血流信号に対してクラッタ信号は数十倍
〜数百倍程度大きく、ウォールフィルタ２２を介さずに
自己相関演算を行えば、血流信号の寄与分を事実上無視
できる。
【００２１】自己相関器３４の出力信号１０６は、速度
演算器３６に入力され、ここでクラッタ周波数ｆｃが演
算されている。なお、この速度演算器３６は、ドプラ演
算部１８に設けられた速度演算器２６と同一の構成を有
する。
演算器３６に入力され、ここでクラッタ周波数ｆｃが演
算されている。なお、この速度演算器３６は、ドプラ演
算部１８に設けられた速度演算器２６と同一の構成を有
する。
【００２２】クラッタ周波数演算部２０に設けられた分
散演算器３８は、直交検波部１６の出力信号１００及び
自己相関器３４の出力１０６に基づいて、クラッタ信号
の分散σc ２を演算するものである。
散演算器３８は、直交検波部１６の出力信号１００及び
自己相関器３４の出力１０６に基づいて、クラッタ信号
の分散σc ２を演算するものである。
【００２３】従って、ドプラ演算部１８からドプラ偏移
周波数ｆｄを表すドプラ演算部出力信号１０４とパワー
｜Ｒ（Ｔ）｜を示す信号とが出力され、クラッタ周波数
演算部２０からクラッタ周波数ｆｃを表す信号と分散σ
c ２を示す信号とが出力される。
周波数ｆｄを表すドプラ演算部出力信号１０４とパワー
｜Ｒ（Ｔ）｜を示す信号とが出力され、クラッタ周波数
演算部２０からクラッタ周波数ｆｃを表す信号と分散σ
c ２を示す信号とが出力される。
【００２４】これらの４つの信号は、クラッタ判別器４
０に入力されている。
０に入力されている。
【００２５】図２には、クラッタ判別範囲２００が示さ
れている。すなわち、横軸がドプラ偏移周波数ｆｄであ
り縦軸がパワー｜Ｒ（Ｔ）｜である。
れている。すなわち、横軸がドプラ偏移周波数ｆｄであ
り縦軸がパワー｜Ｒ（Ｔ）｜である。
【００２６】クラッタ判別範囲２００は、この図２に示
される例において、クラッタ周波数ｆｃを中心としてド
プラ偏移周波数軸方向に一定の幅をもって設定されてお
り、本実施例においてその幅は分散σc ２により設定さ
れる。
される例において、クラッタ周波数ｆｃを中心としてド
プラ偏移周波数軸方向に一定の幅をもって設定されてお
り、本実施例においてその幅は分散σc ２により設定さ
れる。
【００２７】すなわち、クラッタ判別器４０は、クラッ
タ周波数ｆｃを中心として分散σ２に基づいてクラッタ
判別領域２００を設定する。ここで、勿論クラッタ周波
数ｆｃが変動すればそれに追従してこのクラッタ判別領
域２００も左右にシフトする。これは、その幅について
も同様であり分散σ２によりその幅は適応的に可変す
る。
タ周波数ｆｃを中心として分散σ２に基づいてクラッタ
判別領域２００を設定する。ここで、勿論クラッタ周波
数ｆｃが変動すればそれに追従してこのクラッタ判別領
域２００も左右にシフトする。これは、その幅について
も同様であり分散σ２によりその幅は適応的に可変す
る。
【００２８】そして、ドプラ偏移周波数ｆｄとパワー｜
Ｒ（Ｔ）｜とで特定される座標上の点が、このクラッタ
判別領域２００内に入った場合に、クラッタ判別器４０
はクラッタ信号の存在を判別する。そして、除去制御信
号１０８を出力する。
Ｒ（Ｔ）｜とで特定される座標上の点が、このクラッタ
判別領域２００内に入った場合に、クラッタ判別器４０
はクラッタ信号の存在を判別する。そして、除去制御信
号１０８を出力する。
【００２９】なお、本実施例では、図２に示されるよう
に、クラッタ判別範囲２００はパワー｜Ｒ（Ｔ）｜が一
定値ｋ以上に設定されているが、これは図４に示したよ
うにクラッタ信号のレベルが大きいことに基づく。すな
わち、血流信号のみを通過させて、クラッタ信号に対し
て除去を行うために、このようなオフセットが設定され
ている。勿論、後に示すように他の形状を採用すること
も可能である。
に、クラッタ判別範囲２００はパワー｜Ｒ（Ｔ）｜が一
定値ｋ以上に設定されているが、これは図４に示したよ
うにクラッタ信号のレベルが大きいことに基づく。すな
わち、血流信号のみを通過させて、クラッタ信号に対し
て除去を行うために、このようなオフセットが設定され
ている。勿論、後に示すように他の形状を採用すること
も可能である。
【００３０】図１において、クラッタ除去回路２８で
は、除去制御信号１０８が入力された場合、ドプラ出力
をゼロにし、除去制御信号１０８が入力されない場合の
みドプラ演算部出力信号１０４を表示器３０へ通過させ
る。従って、二次元ドプラ画像内において、クラッタの
存在している画素に対しては表示が行われないことにな
り、クラッタによる画像の劣化が防止される。勿論、こ
のような抜けの生じた画素に対しては補間処理などを行
うことも可能である。いずれにしてもこのクラッタ除去
回路２８によって、判別されたクラッタ信号に対する除
去処理が行われる。
は、除去制御信号１０８が入力された場合、ドプラ出力
をゼロにし、除去制御信号１０８が入力されない場合の
みドプラ演算部出力信号１０４を表示器３０へ通過させ
る。従って、二次元ドプラ画像内において、クラッタの
存在している画素に対しては表示が行われないことにな
り、クラッタによる画像の劣化が防止される。勿論、こ
のような抜けの生じた画素に対しては補間処理などを行
うことも可能である。いずれにしてもこのクラッタ除去
回路２８によって、判別されたクラッタ信号に対する除
去処理が行われる。
【００３１】図３には、クラッタ判別範囲の他の例が示
されている。（ａ）に示す例では、クラッタ判別範囲２
０２が逆三角形型に設定されている。そして、その幅が
分散σc ２により設定されている。（ｂ）に示すクラッ
タ判別範囲は逆ガウシアン型であり、その幅が分散σc
２により設定されている。
されている。（ａ）に示す例では、クラッタ判別範囲２
０２が逆三角形型に設定されている。そして、その幅が
分散σc ２により設定されている。（ｂ）に示すクラッ
タ判別範囲は逆ガウシアン型であり、その幅が分散σc
２により設定されている。
【００３２】クラッタ信号は、血流信号に比してパワー
が大きいため、このような形状に設定すれば、血流信号
を除去することなくクラッタ除去処理が実現できる。勿
論、このクラッタ除去範囲の設定に当たってはウォール
フィルタ２２の遮断特性を考慮して定めるのが望まし
い。
が大きいため、このような形状に設定すれば、血流信号
を除去することなくクラッタ除去処理が実現できる。勿
論、このクラッタ除去範囲の設定に当たってはウォール
フィルタ２２の遮断特性を考慮して定めるのが望まし
い。
【００３３】以上の実施例によれば、ウォールフィルタ
で十分に遮断できないようなクラッタ信号が生じた場合
でも、クラッタ信号の平均周波数を中心として極めて限
定された周波数範囲内に対して減衰処理を実行して漏れ
出たクラッタ信号を完全に除去することが可能となる。
この場合、クラッタ信号と血流信号とが重なる範囲にお
いては血流信号についても除去されてしまうが、二次元
ドプラ画像全体としては目障りなクラッタを排除できる
という効果があり、その結果、全体としてきわめて良好
なドプラ画像を形成できる。
で十分に遮断できないようなクラッタ信号が生じた場合
でも、クラッタ信号の平均周波数を中心として極めて限
定された周波数範囲内に対して減衰処理を実行して漏れ
出たクラッタ信号を完全に除去することが可能となる。
この場合、クラッタ信号と血流信号とが重なる範囲にお
いては血流信号についても除去されてしまうが、二次元
ドプラ画像全体としては目障りなクラッタを排除できる
という効果があり、その結果、全体としてきわめて良好
なドプラ画像を形成できる。
【００３４】本実施例においては、ドプラ周波数軸とパ
ワー軸とで構成される二次元座標空間内にクラッタ判別
範囲を設定したが、これには限られずドプラ周波数軸上
のみにすなわち一次元的にクラッタ判別範囲を設定する
ことも可能である。この場合、きわめて急峻な特性を有
する特定周波数内のみをカットするフィルタにより得ら
れる効果と同等の効果が得られる。
ワー軸とで構成される二次元座標空間内にクラッタ判別
範囲を設定したが、これには限られずドプラ周波数軸上
のみにすなわち一次元的にクラッタ判別範囲を設定する
ことも可能である。この場合、きわめて急峻な特性を有
する特定周波数内のみをカットするフィルタにより得ら
れる効果と同等の効果が得られる。
【００３５】本実施例は、ドプラ演算部１８及びクラッ
タ周波数演算部２０においてそれぞれ別個に自己相関器
と速度演算器とを設けたが、時分割方式を採用すること
によって、一方を省略することも可能である。また、必
要に応じて他のパラメータに基づいて適応的にクラッタ
判別範囲を可変させることもできる。
タ周波数演算部２０においてそれぞれ別個に自己相関器
と速度演算器とを設けたが、時分割方式を採用すること
によって、一方を省略することも可能である。また、必
要に応じて他のパラメータに基づいて適応的にクラッタ
判別範囲を可変させることもできる。
【００３６】
【発明の効果】以上説明したように、本発明によれば、
ウォールフィルタを通過してしまったクラッタ信号を効
果的に除去することができ、またそのクラッタの性質に
応じた適切なクラッタ除去を行うことができるという効
果がある。
ウォールフィルタを通過してしまったクラッタ信号を効
果的に除去することができ、またそのクラッタの性質に
応じた適切なクラッタ除去を行うことができるという効
果がある。
【図１】本発明に係る超音波ドプラ診断装置の全体構成
を示すブロック図である。
を示すブロック図である。
【図２】クラッタ判別範囲２００を示す説明図である。
【図３】クラッタ判別範囲の他の例を示す説明図であ
る。
る。
【図４】ウォールフィルタの遮断特性とクラッタ信号と
の関係を示す図である。
の関係を示す図である。
１８ ドプラ演算部 ２０ クラッタ周波数演算部 ２２ ウォールフィルタ ２８ クラッタ除去回路 ３２ パワー演算器 ３８ 分散演算器 ４０ クラッタ判別器
フロントページの続き (56)参考文献 特開 平４−218143（ＪＰ，Ａ) 特開 平５−31111（ＪＰ，Ａ) 特開 平５−7588（ＪＰ，Ａ) 特開 平２−167151（ＪＰ，Ａ) (58)調査した分野(Int.Cl.7，ＤＢ名) A61B 8/00 - 8/15
Claims (5)
- 【請求項１】 生体内運動体からの受信信号からドプラ
偏移周波数を演算するドプラ演算部と、 前記ドプラ演算部への入力前に分岐取り出しされた受信
信号に含まれるクラッタ信号のドプラ偏移周波数である
クラッタ周波数を演算するクラッタ周波数演算部と、 ドプラ偏移周波数軸上に前記クラッタ周波数に応じてク
ラッタ判別範囲を適応的に設定し、前記受信信号のドプ
ラ偏移周波数が前記クラッタ判別範囲に入る場合に、ク
ラッタ信号の存在を判別するクラッタ判別部と、 前記クラッタ信号の存在が判別された場合に、前記ドプ
ラ演算部出力に対してクラッタ信号の除去を行うクラッ
タ除去部と、 を含み、 前記ドプラ演算部に対し、前記クラッタ周波数演算部及
び前記クラッタ判別部が並列に設けられた ことを特徴と
する超音波ドプラ診断装置。 - 【請求項２】 生体内運動体からの受信信号からドプラ
偏移周波数を演算するドプラ演算部と、 前記ドプラ演算部への入力前に分岐取り出しされた受信
信号に含まれるクラッタ信号のドプラ偏移周波数である
クラッタ周波数を演算するクラッタ周波数演算部と、 前記受信信号のドプラ偏移成分のパワーを演算するパワ
ー演算部と、 ドプラ偏移周波数軸とパワー軸とで構成される二次元空
間内に、前記クラッタ周波数に応じてクラッタ判別範囲
を適応的に設定し、前記受信信号のドプラ偏移周波数及
び前記パワーによる空間座標が前記クラッタ判別範囲内
に入る場合に、 クラッタ信号の存在を判別するクラッタ判別部と、 前記クラッタ信号の存在が判別された場合に、前記ドプ
ラ演算部出力に対してクラッタ信号の除去を行うクラッ
タ除去部と、 を含み、 前記ドプラ演算部に対し、前記クラッタ周波数演算部、
前記パワー演算部及び前記クラッタ判別部が並列に設け
られた ことを特徴とする超音波ドプラ診断装置。 - 【請求項３】 請求項１又は２記載の装置において、 前記受信信号中のクラッタ信号の分散を演算するクラッ
タ分散演算器を含み、 前記クラッタ判別範囲は、更に前記分散に対して適応的
に設定され、 前記ドプラ演算部に対し、前記クラッタ分散演算器も並
列に設けられた ことを特徴とする超音波ドプラ診断装
置。 - 【請求項４】 生体内運動体からの受信信号に対して低
域成分を低減するフィルタと、 前記フィルタを通過したフィルタ通過後受信信号からド
プラ偏移周波数を演算するドプラ演算部と、 前記フィルタを通過する前のフィルタ通過前受信信号に
含まれるクラッタ信号のドプラ偏移周波数であるクラッ
タ周波数を演算するクラッタ周波数演算部と、 ドプラ偏移周波数軸上に前記クラッタ周波数に応じてク
ラッタ判別範囲を適応的に設定し、前記フィルタ通過後
受信信号のドプラ偏移周波数が前記クラッタ判別範囲に
入る場合に、クラッタ信号の存在を判別するクラッタ判
別部と、 前記クラッタ信号の存在が判別された場合に、前記ドプ
ラ演算部出力に対してクラッタ信号の除去を行うクラッ
タ除去部と、 を含み、 前記ドプラ演算部に対し、前記クラッタ周波数演算部及
び前記クラッタ判別部が並列に設けられた ことを特徴と
する超音波ドプラ診断装置。 - 【請求項５】 生体内運動体からの受信信号に対して低
域成分を低減するフィルタと、 前記フィルタを通過したフィルタ通過後受信信号からド
プラ偏移周波数を演算するドプラ演算部と、 前記フィルタを通過する前のフィルタ通過前受信信号に
含まれるクラッタ信号のドプラ偏移周波数であるクラッ
タ周波数を演算するクラッタ周波数演算部と、 前記フィルタ通過後受信信号のドプラ偏移成分のパワー
を演算するパワー演算部と、 ドプラ偏移周波数軸とパワー軸とで構成される二次元空
間内に、前記クラッタ周波数に応じてクラッタ判別範囲
を適応的に設定し、前記フィルタ通過後受信信号のドプ
ラ偏移周波数及び前記パワーによる空間座標が前記クラ
ッタ判別範囲内に入る場合に、クラッタ信号の存在を判
別するクラッタ判別部と、 前記クラッタ信号の存在が判別された場合に、前記ドプ
ラ演算部出力に対してクラッタ信号の除去を行うクラッ
タ除去部と、 を含み、 前記ドプラ演算部に対し、前記クラッタ周波数演算部、
前記パワー演算部及び前記クラッタ判別部が並列に設け
られた ことを特徴とする超音波ドプラ診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6148885A JP3042964B2 (ja) | 1994-06-30 | 1994-06-30 | 超音波ドプラ診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6148885A JP3042964B2 (ja) | 1994-06-30 | 1994-06-30 | 超音波ドプラ診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0810258A JPH0810258A (ja) | 1996-01-16 |
| JP3042964B2 true JP3042964B2 (ja) | 2000-05-22 |
Family
ID=15462887
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP6148885A Expired - Fee Related JP3042964B2 (ja) | 1994-06-30 | 1994-06-30 | 超音波ドプラ診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3042964B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3887040B2 (ja) * | 1996-09-05 | 2007-02-28 | 株式会社東芝 | 超音波診断装置 |
| KR101055581B1 (ko) * | 2008-03-20 | 2011-08-23 | 삼성메디슨 주식회사 | 초음파 시스템 및 클러터 신호 제거용 필터 설정 방법 |
-
1994
- 1994-06-30 JP JP6148885A patent/JP3042964B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0810258A (ja) | 1996-01-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20050054931A1 (en) | Tracking clutter filter for spectral & audio doppler | |
| JP3397748B2 (ja) | カラードップラー映像システムにおけるカラー映像表示方法及び装置 | |
| JP7405950B2 (ja) | 微小脈管の高時空分解能超音波イメージングを行うための方法 | |
| JP2002534186A (ja) | 超音波カラーフロー表示最適化法 | |
| JP2001178721A (ja) | パケット・データ獲得を用いた超音波フロー撮像における運動の可視化のための方法および装置 | |
| US20060030777A1 (en) | T-statistic method for suppressing artifacts in blood vessel ultrasonic imaging | |
| KR100232257B1 (ko) | 클러터신호의 과도응답을 최소로 하는 초음파칼라도플러영상시스템 | |
| JPH08201500A (ja) | 水中音響信号検出装置 | |
| Lovstakken et al. | Blood flow imaging-a new real-time, flow imaging technique | |
| JP3042964B2 (ja) | 超音波ドプラ診断装置 | |
| US8189427B2 (en) | Clutter signal filtering for doppler signal | |
| JP3773779B2 (ja) | レーダ信号処理装置 | |
| CN109363722B (zh) | 彩色血流成像中运动伪像的抑制方法及设备 | |
| JP2003529422A (ja) | Ｂカラー優先閾値計算を用いた超音波システム | |
| JP3188394B2 (ja) | 超音波カラードップラー映像システムのためのフィルタリング方法及びそのフィルタ回路 | |
| JP2002532172A (ja) | 組織内の壁側面及び動脈の絶対半径の決定のための超音波方法及び装置 | |
| JP4067914B2 (ja) | 超音波診断装置 | |
| JP2723458B2 (ja) | 超音波ドプラ診断装置 | |
| JP2720711B2 (ja) | 超音波二次元ドプラ血流計 | |
| JP3408647B2 (ja) | 超音波イメージング装置 | |
| JP2003250802A (ja) | 超音波診断装置 | |
| JP2005318921A (ja) | 超音波診断装置 | |
| JP6828017B2 (ja) | 適応型サンプル窓サイズを用いるスペクトルドプラ処理 | |
| JPH08150142A (ja) | 超音波診断装置 | |
| JPH04250148A (ja) | 超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090310 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110310 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110310 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130310 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140310 Year of fee payment: 14 |
|
| LAPS | Cancellation because of no payment of annual fees |