JP3042964B2 - Ultrasound Doppler diagnostic device - Google Patents
Ultrasound Doppler diagnostic deviceInfo
- 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
Description
【０００１】[0001]
【産業上の利用分野】本発明は、超音波ドプラ診断装
置、特に受信信号に含まれるクラッタ信号の除去に関す
る。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an ultrasonic Doppler diagnostic apparatus and, more particularly, to the removal of clutter signals contained in a received signal.
【０００２】[0002]
【従来の技術】超音波ドプラ診断装置は、ドプラ偏移周
波数に基づき血流等の生体内運動体の運動情報（速度
等）を画像表示する装置である。例えば、心臓のドプラ
画像を形成する場合、本来観察したい血流のほかに、動
きのある心臓壁等の目障りなクラッタが表示されてしま
う。心臓壁等の低速運動体からの不要反射波は、血流に
比べ非常に強く、このため受信信号中でクラッタ信号の
レベルは大きい。2. Description of the Related Art An ultrasonic Doppler diagnostic apparatus is an apparatus for displaying an image of motion information (speed, etc.) of a moving body in a living body such as a blood flow based on a Doppler shift frequency. For example, when a Doppler image of the heart is formed, an unsightly clutter, such as a moving heart wall, is displayed in addition to the blood flow that is originally desired to be observed. Unwanted reflected waves from a low-speed moving body such as a heart wall are much stronger than blood flows, and therefore the level of a clutter signal is large in a received signal.
【０００３】そこで、従来の超音波ドプラ診断装置で
は、いわゆるウォールフィルタと呼ばれるクラッタ除去
フィルタが設けられている。このウォールフィルタは、
ローカットフィルタであり、低速運動体のドプラ成分を
低減するものである。ここで図４には、ウォールフィル
タの遮断特性が示されている。Therefore, a conventional ultrasonic Doppler diagnostic apparatus is provided with a clutter removing filter called a so-called wall filter. This wall filter is
This is a low-cut filter that reduces the Doppler component of the low-speed moving body. Here, FIG. 4 shows the cutoff characteristics of the wall filter.
【０００４】[0004]
【発明が解決しようとする課題】従って、いま図４の低
い周波数のクラッタａが入力した場合、ウォールフィル
タによってクラッタのレベルはａ´のように低減され
る。Therefore, when the low frequency clutter a shown in FIG. 4 is input, the level of the clutter is reduced to a 'by the wall filter.
【０００５】しかし、クラッタｂのような高い周波数の
クラッタが入力した場合、ウォールフィルタでは低減で
きず、そのまま通過してしまう。そして、これが画像の
劣化を引き起こす。However, when a high-frequency clutter such as the clutter b is input, it cannot be reduced by the wall filter and passes through as it is. This causes image degradation.
【０００６】なお、ウォールフィルタの遮断特性は、ク
ラッタ信号及び血流信号のドプラ周波数に応じて選択で
きるが、カットオフ周波数を上げすぎると、血流信号が
損なわれる。また、血流信号への影響をできる限り少な
くするにはフィルタ特性を急峻にする必要があるが、超
音波ドプラ診断装置のような場合には、処理時間の制約
からそれには限界がある。The cutoff characteristic of the wall filter can be selected according to the Doppler frequency of the clutter signal and the blood flow signal. However, if the cutoff frequency is too high, the blood flow signal will be damaged. In order to minimize the influence on the blood flow signal, it is necessary to make the filter characteristics sharp. However, in the case of an ultrasonic Doppler diagnostic apparatus, there is a limit due to processing time restrictions.
【０００７】そこで、本発明は、効果的にクラッタの除
去を行うことを目的とする。Therefore, an object of the present invention is to effectively remove clutter.
【０００８】また、本発明は、クラッタの性質に応じて
適切なクラッタ除去を行うことを目的とする。Another object of the present invention is to perform appropriate clutter removal according to the nature of clutter.
【０００９】[0009]
【課題を解決するための手段】上記目的を達成するため
に、請求項１記載の発明は、生体内運動体からの受信信
号からドプラ偏移周波数を演算するドプラ演算部と、前
記ドプラ演算部への入力前に分岐取り出しされた受信信
号に含まれるクラッタ信号のドプラ偏移周波数であるク
ラッタ周波数を演算するクラッタ周波数演算部と、ドプ
ラ偏移周波数軸上に前記クラッタ周波数に応じてクラッ
タ判別範囲を適応的に設定し、前記受信信号のドプラ偏
移周波数が前記クラッタ判別範囲に入る場合に、クラッ
タ信号の存在を判別するクラッタ判別部と、前記クラッ
タ信号の存在が判別された場合に、前記ドプラ演算部出
力に対してクラッタ信号の除去を行うクラッタ除去部
と、を含み、前記ドプラ演算部に対し、前記クラッタ周
波数演算部及び前記クラッタ判別部が並列に設けられた
ことを特徴とする。According to a first aspect of the present invention, there is provided a Doppler operation unit for calculating a Doppler shift frequency from a signal received from a moving body in a living body, and the Doppler operation unit. A clutter frequency calculating unit for calculating a clutter frequency which is a Doppler shift frequency of a clutter signal included in a received signal that is branched and extracted before input to a clutter determination range according to the clutter frequency on a Doppler shift frequency axis When the Doppler shift frequency of the received signal falls within the clutter determination range, a clutter determination unit that determines the presence of a clutter signal, and when the presence of the clutter signal is determined, the It is seen containing a clutter removing unit for removing the clutter signals to the Doppler calculation section outputs, a, relative to the Doppler calculation section, the clutter circumferential
A wave number calculation unit and the clutter determination unit are provided in parallel .
【００１０】請求項２記載の発明は、生体内運動体から
の受信信号からドプラ偏移周波数を演算するドプラ演算
部と、前記ドプラ演算部への入力前に分岐取り出しされ
た受信信号に含まれるクラッタ信号のドプラ偏移周波数
であるクラッタ周波数を演算するクラッタ周波数演算部
と、前記受信信号のドプラ偏移成分のパワーを演算する
パワー演算部と、ドプラ偏移周波数軸とパワー軸とで構
成される二次元空間内に、前記クラッタ周波数に応じて
クラッタ判別範囲を適応的に設定し、前記受信信号のド
プラ偏移周波数及び前記パワーによる空間座標が前記ク
ラッタ判別範囲内に入る場合に、クラッタ信号の存在を
判別するクラッタ判別部と、前記クラッタ信号の存在が
判別された場合に、前記ドプラ演算部出力に対してクラ
ッタ信号の除去を行うクラッタ除去部と、を含み、前記
ドプラ演算部に対し、前記クラッタ周波数演算部、前記
パワー演算部及び前記クラッタ判別部が並列に設けられ
たことを特徴とする。According to a second aspect of the present invention, there is provided a Doppler operation unit for calculating a Doppler shift frequency from a received signal from a moving body in a living body, and a Doppler operation unit which branches and extracts the Doppler operation frequency before input to the Doppler operation unit.
A clutter frequency calculator that calculates a clutter frequency that is a Doppler shift frequency of a clutter signal included in the received signal, a power calculator that calculates the power of the Doppler shift component of the received signal, and a Doppler shift frequency axis. In a two-dimensional space composed of a power axis and a clutter discrimination range, the clutter discrimination range is adaptively set according to the clutter frequency, and the Doppler shift frequency of the received signal and the spatial coordinates of the power are within the clutter discrimination range. when entering the clutter discrimination unit for discriminating the presence of the clutter signal, when the presence of the clutter signal is determined, see containing and a clutter removing unit for removing the clutter signals to the Doppler calculation section outputs And said
For the Doppler calculation unit, the clutter frequency calculation unit,
A power calculation unit and the clutter determination unit are provided in parallel.
Characterized in that was.
【００１１】請求項３記載の発明は、前記受信信号中の
クラッタ信号の分散を演算するクラッタ分散演算器を含
み、前記クラッタ判別範囲は、更に前記分散に対して適
応的に設定され、前記ドプラ演算部に対し、前記クラッ
タ分散演算器も並列に設けられたことを特徴とする。ま
た、請求項４及び請求項５の発明は、低域通過フィルタ
を通過する前の受信信号からクラッタ周波数を演算する
ことを特徴とする。[0011] According to a third aspect, includes a clutter dispersion calculator for calculating the variance of the clutter signal in the received signal, the clutter determination range is adaptively set for further the dispersion, the Doppler The calculation unit
The distributed computing units are also provided in parallel . The invention according to claims 4 and 5 is characterized in that a clutter frequency is calculated from a received signal before passing through a low-pass filter.
【００１２】[0012]
【作用】上記構成によれば、受信信号に含まれるクラッ
タ周波数（例えば、クラッタスペクトルの平均周波数）
が演算され、そのクラッタ周波数に応じてクラッタ判別
範囲が設定される。そして、そのクラッタ判別範囲内
に、受信信号のドプラ偏移周波数が入った場合、クラッ
タの存在が判別され、ドプラ演算部出力に対してクラッ
タ除去処理が実行される。According to the above arrangement, the clutter frequency included in the received signal (for example, the average frequency of the clutter spectrum)
Is calculated, and a clutter determination range is set according to the clutter frequency. Then, when the Doppler shift frequency of the received signal falls within the clutter determination range, the presence of clutter is determined, and clutter removal processing is performed on the output of the Doppler calculation unit.
【００１３】ここで、クラッタ判別範囲は、ドプラ偏移
周波数上に設定できるほか、ドプラ偏移周波数軸とパワ
ー軸とで構成される二次元空間内に設定できる。二次元
空間内に設定すれば、ドプラ偏移成分のパワーを考慮し
て、より精度良く判別を行うことができる。更に、分散
に応じてクラッタ判別範囲を変化させて最適なクラッタ
除去処理を実現できる。Here, the clutter discrimination range can be set on the Doppler shift frequency and can be set in a two-dimensional space formed by the Doppler shift frequency axis and the power axis. If it is set in a two-dimensional space, the discrimination can be performed with higher accuracy in consideration of the power of the Doppler shift component. Further, an optimal clutter removal process can be realized by changing the clutter determination range according to the dispersion.
【００１４】[0014]
【実施例】以下、本発明の好適な実施例を図面に基づい
て説明する。DESCRIPTION OF THE PREFERRED EMBODIMENTS Preferred embodiments of the present invention will be described below with reference to the drawings.
【００１５】図１には、本発明に係る超音波ドプラ診断
装置の好適な実施例が示されており、図１はその全体構
成図である。図１において、送受信器１２の駆動の下、
超音波探触子１０によって生体１４に対して超音波パル
スが放射され、生体内の各運動反射体にてドプラシフト
を受けた反射波が超音波探触子１０にて受波される。こ
れにより得られる受信信号は、送受信器１２を介して直
交検波部１６に送られる。直交検波部１６は、参照信号
を用いて受信信号に対して直交検波を行い、直交検波後
の信号１００を出力する。FIG. 1 shows a preferred embodiment of an ultrasonic Doppler diagnostic apparatus according to the present invention, and FIG. 1 is an overall configuration diagram thereof. In FIG. 1, under the driving of the transceiver 12,
Ultrasonic pulses are emitted to the living body 14 by the ultrasonic probe 10, and reflected waves that have undergone Doppler shift at each motion reflector in the living body are received by the ultrasonic probe 10. The received signal thus obtained is sent to the quadrature detector 16 via the transceiver 12. The quadrature detection unit 16 performs quadrature detection on the received signal using the reference signal, and outputs a signal 100 after quadrature detection.
【００１６】この信号１００は、ドプラ演算部１８とク
ラッタ周波数演算部２０に送られている。つまり、ドプ
ラ演算部１８の手前から分岐取り出しされた信号がクラ
ッタ周波数演算部２０に入力されている。 The signal 100 is sent to a Doppler operation unit 18 and a clutter frequency operation unit 20. In other words, dope
The signal branched and taken out from just before the
The data is input to the frequency calculator 20.
【００１７】ドプラ演算部１８において、信号１００
は、まずクラッタ信号を除去するウォールフィルタ２２
に入力され、ここで図４に示した遮断特性に従いクラッ
タ信号が低減される。そして、その出力信号は周知の自
己相関器２４に入力され、自己相関演算が実行され、自
己相関信号１０２が得られる。In the Doppler operation section 18, the signal 100
First, a wall filter 22 for removing clutter signals
, Where the clutter signal is reduced according to the cutoff characteristics shown in FIG. Then, the output signal is input to a well-known autocorrelator 24, where an autocorrelation operation is executed, and an autocorrelation signal 102 is obtained.
【００１８】速度演算器２６は、前記自己相関信号１０
２の実数部と虚数部の逆正接演算によって速度を求める
ものであり、ドプラ偏移周波数ｆｄを表すドプラ演算部
出力信号１０４が出力されている。そして、この信号１
０４はクラッタ除去回路２８を介して表示器３０に送ら
れており、この表示器３０には二次元ドプラ画像が表示
される。The speed calculator 26 calculates the speed of the autocorrelation signal 10
The arctangent calculation of two real and imaginary parts are those seeking speed, Doppler operation unit output signal 104 representative of the Doppler shift frequency f d is outputted. And this signal 1
04 is sent to the display 30 via the clutter removing circuit 28, and the display 30 displays a two-dimensional Doppler image.
【００１９】ドプラ演算部１８に設けられたパワー演算
器３２は、自己相関器２４の出力信号１０２に基づいて
受信信号中のドプラ偏移成分のパワーを演算するもので
あり、演算されたパワー｜Ｒ（Ｔ）｜が出力されてい
る。The power calculator 32 provided in the Doppler calculator 18 calculates the power of the Doppler shift component in the received signal based on the output signal 102 of the autocorrelator 24. The calculated power | R (T) | is output.
【００２０】前記クラッタ周波数演算部２０において、
直交検波部１６の出力信号１００は、自己相関器３４に
入力され、自己相関演算が実行されている。この自己相
関器３４は、上述した自己相関器２４と同一の構成を有
するものであるが、その前段にウォールフィルタ２２が
設けられていないので、自己相関器３４の出力はほぼク
ラッタ信号に対する自己相関演算結果を示すことにな
る。すなわち、血流信号に対してクラッタ信号は数十倍
〜数百倍程度大きく、ウォールフィルタ２２を介さずに
自己相関演算を行えば、血流信号の寄与分を事実上無視
できる。In the clutter frequency calculating section 20,
The output signal 100 of the quadrature detector 16 is input to the autocorrelator 34, where an autocorrelation operation is performed. The auto-correlator 34 has the same configuration as the auto-correlator 24 described above, but since the wall filter 22 is not provided at the preceding stage, the output of the auto-correlator 34 is substantially equal to the auto-correlator for the clutter signal. The result of the operation will be shown. That is, the clutter signal is several tens to several hundred times larger than the blood flow signal, and if the autocorrelation calculation is performed without passing through the wall filter 22, the contribution of the blood flow signal can be ignored.
【００２１】自己相関器３４の出力信号１０６は、速度
演算器３６に入力され、ここでクラッタ周波数ｆｃが演
算されている。なお、この速度演算器３６は、ドプラ演
算部１８に設けられた速度演算器２６と同一の構成を有
する。The output signal 106 of the autocorrelator 34 is inputted to the speed calculator 36, where the clutter frequency f c is computed. The speed calculator 36 has the same configuration as the speed calculator 26 provided in the Doppler calculator 18.
【００２２】クラッタ周波数演算部２０に設けられた分
散演算器３８は、直交検波部１６の出力信号１００及び
自己相関器３４の出力１０６に基づいて、クラッタ信号
の分散σc ２を演算するものである。A variance calculator 38 provided in the clutter frequency calculator 20 calculates the variance σ c 2 of the clutter signal based on the output signal 100 of the quadrature detector 16 and the output 106 of the autocorrelator 34. is there.
【００２３】従って、ドプラ演算部１８からドプラ偏移
周波数ｆｄを表すドプラ演算部出力信号１０４とパワー
｜Ｒ（Ｔ）｜を示す信号とが出力され、クラッタ周波数
演算部２０からクラッタ周波数ｆｃを表す信号と分散σ
c ２を示す信号とが出力される。[0023] Thus, the Doppler calculation section output signal 104 and the power from the Doppler calculation section 18 representing the Doppler shift frequency f d | R (T) | a signal indicating a is outputted, the clutter frequency f c from the clutter frequency calculating unit 20 And the variance σ
signal are output to indicate c 2.
【００２４】これらの４つの信号は、クラッタ判別器４
０に入力されている。These four signals are output to a clutter discriminator 4
0 has been entered.
【００２５】図２には、クラッタ判別範囲２００が示さ
れている。すなわち、横軸がドプラ偏移周波数ｆｄであ
り縦軸がパワー｜Ｒ（Ｔ）｜である。FIG. 2 shows a clutter determination range 200. In other words, the horizontal axis represents the Doppler shift frequency f d and the vertical axis power | R (T) | is.
【００２６】クラッタ判別範囲２００は、この図２に示
される例において、クラッタ周波数ｆｃを中心としてド
プラ偏移周波数軸方向に一定の幅をもって設定されてお
り、本実施例においてその幅は分散σc ２により設定さ
れる。The clutter determining range 200, in the example shown in FIG. 2, are set to have a predetermined width in the Doppler shift frequency axis direction about the clutter frequency fc, the width in this embodiment variance sigma c 2 is set.
【００２７】すなわち、クラッタ判別器４０は、クラッ
タ周波数ｆｃを中心として分散σ２に基づいてクラッタ
判別領域２００を設定する。ここで、勿論クラッタ周波
数ｆｃが変動すればそれに追従してこのクラッタ判別領
域２００も左右にシフトする。これは、その幅について
も同様であり分散σ２によりその幅は適応的に可変す
る。[0027] That is, clutter discriminator 40 sets the clutter discrimination region 200 based on the variance sigma 2 around the clutter frequency f c. Here, even the clutter determination area 200 to follow it if change course clutter frequency f c is shifted to the left and right. It has a width by and variance sigma 2 also applies to the width variable adaptively.
【００２８】そして、ドプラ偏移周波数ｆｄとパワー｜
Ｒ（Ｔ）｜とで特定される座標上の点が、このクラッタ
判別領域２００内に入った場合に、クラッタ判別器４０
はクラッタ信号の存在を判別する。そして、除去制御信
号１０８を出力する。Then, Doppler shift frequency fd and power |
When a point on the coordinates specified by R (T) | enters this clutter determination area 200, the clutter discriminator 40
Determines the presence of a clutter signal. Then, a removal control signal 108 is output.
【００２９】なお、本実施例では、図２に示されるよう
に、クラッタ判別範囲２００はパワー｜Ｒ（Ｔ）｜が一
定値ｋ以上に設定されているが、これは図４に示したよ
うにクラッタ信号のレベルが大きいことに基づく。すな
わち、血流信号のみを通過させて、クラッタ信号に対し
て除去を行うために、このようなオフセットが設定され
ている。勿論、後に示すように他の形状を採用すること
も可能である。In this embodiment, as shown in FIG. 2, in the clutter determination range 200, the power | R (T) | is set to a constant value k or more, as shown in FIG. This is based on the fact that the level of the clutter signal is large. That is, such an offset is set in order to pass only the blood flow signal and remove the clutter signal. Of course, other shapes can be employed as shown below.
【００３０】図１において、クラッタ除去回路２８で
は、除去制御信号１０８が入力された場合、ドプラ出力
をゼロにし、除去制御信号１０８が入力されない場合の
みドプラ演算部出力信号１０４を表示器３０へ通過させ
る。従って、二次元ドプラ画像内において、クラッタの
存在している画素に対しては表示が行われないことにな
り、クラッタによる画像の劣化が防止される。勿論、こ
のような抜けの生じた画素に対しては補間処理などを行
うことも可能である。いずれにしてもこのクラッタ除去
回路２８によって、判別されたクラッタ信号に対する除
去処理が行われる。In FIG. 1, the clutter removal circuit 28 sets the Doppler output to zero when the removal control signal 108 is input, and passes the Doppler operation unit output signal 104 to the display 30 only when the removal control signal 108 is not input. Let it. Therefore, in the two-dimensional Doppler image, display is not performed on the pixel where the clutter exists, and the deterioration of the image due to the clutter is prevented. Of course, it is also possible to perform an interpolation process or the like on such a missing pixel. In any case, the clutter removing circuit 28 removes the determined clutter signal.
【００３１】図３には、クラッタ判別範囲の他の例が示
されている。（ａ）に示す例では、クラッタ判別範囲２
０２が逆三角形型に設定されている。そして、その幅が
分散σc ２により設定されている。（ｂ）に示すクラッ
タ判別範囲は逆ガウシアン型であり、その幅が分散σc
２により設定されている。FIG. 3 shows another example of the clutter determination range. In the example shown in FIG.
02 is set to an inverted triangle type. The width is set by the variance σ c 2 . The clutter determination range shown in (b) is of the inverse Gaussian type, and its width has a variance σ c
2 is set.
【００３２】クラッタ信号は、血流信号に比してパワー
が大きいため、このような形状に設定すれば、血流信号
を除去することなくクラッタ除去処理が実現できる。勿
論、このクラッタ除去範囲の設定に当たってはウォール
フィルタ２２の遮断特性を考慮して定めるのが望まし
い。Since the power of the clutter signal is greater than that of the blood flow signal, by setting such a shape, clutter removal processing can be realized without removing the blood flow signal. Needless to say, it is desirable to set the clutter removal range in consideration of the cutoff characteristics of the wall filter 22.
【００３３】以上の実施例によれば、ウォールフィルタ
で十分に遮断できないようなクラッタ信号が生じた場合
でも、クラッタ信号の平均周波数を中心として極めて限
定された周波数範囲内に対して減衰処理を実行して漏れ
出たクラッタ信号を完全に除去することが可能となる。
この場合、クラッタ信号と血流信号とが重なる範囲にお
いては血流信号についても除去されてしまうが、二次元
ドプラ画像全体としては目障りなクラッタを排除できる
という効果があり、その結果、全体としてきわめて良好
なドプラ画像を形成できる。According to the above-described embodiment, even when a clutter signal that cannot be sufficiently cut off by the wall filter is generated, the attenuation process is performed within a very limited frequency range centered on the average frequency of the clutter signal. Thus, it is possible to completely remove the leaked clutter signal.
In this case, the blood flow signal is also removed in a range where the clutter signal and the blood flow signal overlap, but there is an effect that unpleasant clutter can be eliminated as a whole of the two-dimensional Doppler image, and as a result, the overall effect is extremely low. A good Doppler image can be formed.
【００３４】本実施例においては、ドプラ周波数軸とパ
ワー軸とで構成される二次元座標空間内にクラッタ判別
範囲を設定したが、これには限られずドプラ周波数軸上
のみにすなわち一次元的にクラッタ判別範囲を設定する
ことも可能である。この場合、きわめて急峻な特性を有
する特定周波数内のみをカットするフィルタにより得ら
れる効果と同等の効果が得られる。In the present embodiment, the clutter determination range is set in the two-dimensional coordinate space composed of the Doppler frequency axis and the power axis. However, the present invention is not limited to this. It is also possible to set a clutter determination range. In this case, an effect equivalent to that obtained by a filter that cuts only within a specific frequency having extremely steep characteristics can be obtained.
【００３５】本実施例は、ドプラ演算部１８及びクラッ
タ周波数演算部２０においてそれぞれ別個に自己相関器
と速度演算器とを設けたが、時分割方式を採用すること
によって、一方を省略することも可能である。また、必
要に応じて他のパラメータに基づいて適応的にクラッタ
判別範囲を可変させることもできる。In this embodiment, the autocorrelator and the speed calculator are separately provided in the Doppler calculation unit 18 and the clutter frequency calculation unit 20, but one of them may be omitted by adopting the time division method. It is possible. Also, the clutter determination range can be adaptively varied based on other parameters as needed.
【００３６】[0036]
【発明の効果】以上説明したように、本発明によれば、
ウォールフィルタを通過してしまったクラッタ信号を効
果的に除去することができ、またそのクラッタの性質に
応じた適切なクラッタ除去を行うことができるという効
果がある。As described above, according to the present invention,
There is an effect that clutter signals that have passed through the wall filter can be effectively removed, and that clutter can be appropriately removed according to the nature of the clutter.
【図１】本発明に係る超音波ドプラ診断装置の全体構成
を示すブロック図である。FIG. 1 is a block diagram showing an overall configuration of an ultrasonic Doppler diagnostic apparatus according to the present invention.
【図２】クラッタ判別範囲２００を示す説明図である。FIG. 2 is an explanatory diagram showing a clutter determination range 200.
【図３】クラッタ判別範囲の他の例を示す説明図であ
る。FIG. 3 is an explanatory diagram illustrating another example of a clutter determination range.
【図４】ウォールフィルタの遮断特性とクラッタ信号と
の関係を示す図である。FIG. 4 is a diagram illustrating a relationship between a cutoff characteristic of a wall filter and a clutter signal.
１８ ドプラ演算部 ２０ クラッタ周波数演算部 ２２ ウォールフィルタ ２８ クラッタ除去回路 ３２ パワー演算器 ３８ 分散演算器 ４０ クラッタ判別器 18 Doppler operation unit 20 Clutter frequency operation unit 22 Wall filter 28 Clutter removal circuit 32 Power operation unit 38 Distributed operation unit 40 Clutter discrimination unit
フロントページの続き (56)参考文献 特開 平４−218143（ＪＰ，Ａ) 特開 平５−31111（ＪＰ，Ａ) 特開 平５−7588（ＪＰ，Ａ) 特開 平２−167151（ＪＰ，Ａ) (58)調査した分野(Int.Cl.7，ＤＢ名) A61B 8/00 - 8/15 Continuation of the front page (56) References JP-A-4-218143 (JP, A) JP-A-5-31111 (JP, A) JP-A-5-7588 (JP, A) JP-A-2-167151 (JP) , A) (58) Field surveyed (Int. Cl. 7 , DB name) A61B 8/00-8/15
Claims (5)
偏移周波数を演算するドプラ演算部と、 前記ドプラ演算部への入力前に分岐取り出しされた受信
信号に含まれるクラッタ信号のドプラ偏移周波数である
クラッタ周波数を演算するクラッタ周波数演算部と、 ドプラ偏移周波数軸上に前記クラッタ周波数に応じてク
ラッタ判別範囲を適応的に設定し、前記受信信号のドプ
ラ偏移周波数が前記クラッタ判別範囲に入る場合に、ク
ラッタ信号の存在を判別するクラッタ判別部と、 前記クラッタ信号の存在が判別された場合に、前記ドプ
ラ演算部出力に対してクラッタ信号の除去を行うクラッ
タ除去部と、 を含み、 前記ドプラ演算部に対し、前記クラッタ周波数演算部及
び前記クラッタ判別部が並列に設けられた ことを特徴と
する超音波ドプラ診断装置。1. A Doppler calculator for calculating a Doppler shift frequency from a received signal from a moving body in a living body, and a Doppler shift of a clutter signal included in a received signal branched and extracted before input to the Doppler calculator. A clutter frequency calculating unit that calculates a clutter frequency that is a frequency, and a clutter determination range is adaptively set on a Doppler shift frequency axis according to the clutter frequency, and the Doppler shift frequency of the received signal is the clutter determination range. A clutter discriminating unit that discriminates the presence of a clutter signal when entering, and a clutter removing unit that removes a clutter signal from the Doppler operation unit output when the presence of the clutter signal is discriminated. The clutter frequency calculation unit and the Doppler calculation unit
And an ultrasonic Doppler diagnostic apparatus, wherein the clutter determination unit is provided in parallel .
偏移周波数を演算するドプラ演算部と、 前記ドプラ演算部への入力前に分岐取り出しされた受信
信号に含まれるクラッタ信号のドプラ偏移周波数である
クラッタ周波数を演算するクラッタ周波数演算部と、 前記受信信号のドプラ偏移成分のパワーを演算するパワ
ー演算部と、 ドプラ偏移周波数軸とパワー軸とで構成される二次元空
間内に、前記クラッタ周波数に応じてクラッタ判別範囲
を適応的に設定し、前記受信信号のドプラ偏移周波数及
び前記パワーによる空間座標が前記クラッタ判別範囲内
に入る場合に、 クラッタ信号の存在を判別するクラッタ判別部と、 前記クラッタ信号の存在が判別された場合に、前記ドプ
ラ演算部出力に対してクラッタ信号の除去を行うクラッ
タ除去部と、 を含み、 前記ドプラ演算部に対し、前記クラッタ周波数演算部、
前記パワー演算部及び前記クラッタ判別部が並列に設け
られた ことを特徴とする超音波ドプラ診断装置。2. A Doppler operation unit for calculating a Doppler shift frequency from a received signal from a moving body in a living body, and a Doppler shift of a clutter signal included in a received signal branched and extracted before input to the Doppler operation unit. A clutter frequency calculator that calculates a clutter frequency that is a frequency, a power calculator that calculates the power of the Doppler shift component of the received signal, and a two-dimensional space that includes a Doppler shift frequency axis and a power axis. A clutter discriminating range is adaptively set according to the clutter frequency, and a clutter for discriminating the presence of a clutter signal when the Doppler shift frequency of the received signal and the spatial coordinates based on the power fall within the clutter discriminating range. A determining unit, and a clutter removing unit that removes a clutter signal from an output of the Doppler operation unit when the presence of the clutter signal is determined. , Only including, with respect to the Doppler calculation section, the clutter frequency calculating unit,
The power calculation unit and the clutter determination unit are provided in parallel.
An ultrasonic Doppler diagnostic apparatus, characterized in that:
タ分散演算器を含み、 前記クラッタ判別範囲は、更に前記分散に対して適応的
に設定され、 前記ドプラ演算部に対し、前記クラッタ分散演算器も並
列に設けられた ことを特徴とする超音波ドプラ診断装
置。3. The apparatus according to claim 1, further comprising a clutter dispersion calculator for computing a variance of a clutter signal in the received signal, wherein the clutter determination range is further set adaptively with respect to the variance. It is, with respect to the Doppler calculation section, also parallel the clutter dispersion calculator
An ultrasonic Doppler diagnostic apparatus provided in a row .
域成分を低減するフィルタと、 前記フィルタを通過したフィルタ通過後受信信号からド
プラ偏移周波数を演算するドプラ演算部と、 前記フィルタを通過する前のフィルタ通過前受信信号に
含まれるクラッタ信号のドプラ偏移周波数であるクラッ
タ周波数を演算するクラッタ周波数演算部と、 ドプラ偏移周波数軸上に前記クラッタ周波数に応じてク
ラッタ判別範囲を適応的に設定し、前記フィルタ通過後
受信信号のドプラ偏移周波数が前記クラッタ判別範囲に
入る場合に、クラッタ信号の存在を判別するクラッタ判
別部と、 前記クラッタ信号の存在が判別された場合に、前記ドプ
ラ演算部出力に対してクラッタ信号の除去を行うクラッ
タ除去部と、 を含み、 前記ドプラ演算部に対し、前記クラッタ周波数演算部及
び前記クラッタ判別部が並列に設けられた ことを特徴と
する超音波ドプラ診断装置。4. A filter for reducing a low-frequency component of a received signal from a moving body in a living body, a Doppler operation unit for calculating a Doppler shift frequency from a received signal after passing through the filter, and the filter. A clutter frequency calculation unit that calculates a clutter frequency that is a Doppler shift frequency of a clutter signal included in the received signal before passing through the filter before the filter, and a clutter determination range on the Doppler shift frequency axis according to the clutter frequency. Adaptively set, when the Doppler shift frequency of the received signal after passing through the filter falls within the clutter determination range, a clutter determination unit that determines the presence of a clutter signal, and when the presence of the clutter signal is determined. the a clutter removing unit for removing the clutter signals to the Doppler calculation section outputs, only including, with respect to the Doppler calculation section, before Clutter frequency calculating unit及
And an ultrasonic Doppler diagnostic apparatus, wherein the clutter determination unit is provided in parallel .
域成分を低減するフィルタと、 前記フィルタを通過したフィルタ通過後受信信号からド
プラ偏移周波数を演算するドプラ演算部と、 前記フィルタを通過する前のフィルタ通過前受信信号に
含まれるクラッタ信号のドプラ偏移周波数であるクラッ
タ周波数を演算するクラッタ周波数演算部と、 前記フィルタ通過後受信信号のドプラ偏移成分のパワー
を演算するパワー演算部と、 ドプラ偏移周波数軸とパワー軸とで構成される二次元空
間内に、前記クラッタ周波数に応じてクラッタ判別範囲
を適応的に設定し、前記フィルタ通過後受信信号のドプ
ラ偏移周波数及び前記パワーによる空間座標が前記クラ
ッタ判別範囲内に入る場合に、クラッタ信号の存在を判
別するクラッタ判別部と、 前記クラッタ信号の存在が判別された場合に、前記ドプ
ラ演算部出力に対してクラッタ信号の除去を行うクラッ
タ除去部と、 を含み、 前記ドプラ演算部に対し、前記クラッタ周波数演算部、
前記パワー演算部及び前記クラッタ判別部が並列に設け
られた ことを特徴とする超音波ドプラ診断装置。5. A filter for reducing low-frequency components of a received signal from a moving body in a living body, a Doppler operation unit for calculating a Doppler shift frequency from a received signal after passing through the filter, and the filter. A clutter frequency calculating unit that calculates a clutter frequency that is a Doppler shift frequency of a clutter signal included in the received signal before passing through the filter, and a power that calculates the power of the Doppler shifted component of the received signal after passing through the filter. An arithmetic unit, in a two-dimensional space composed of a Doppler shift frequency axis and a power axis, a clutter determination range is adaptively set according to the clutter frequency, and the Doppler shift frequency of the received signal after passing through the filter. And when a spatial coordinate based on the power falls within the clutter determination range, a clutter determination unit that determines the presence of a clutter signal. If the presence of the clutter signal is discriminated, the a clutter removing unit for removing the clutter signals to the Doppler calculation section outputs, only including, with respect to the Doppler calculation section, the clutter frequency calculating unit,
The power calculation unit and the clutter determination unit are provided in parallel.
An ultrasonic Doppler diagnostic apparatus, characterized in that:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6148885A JP3042964B2 (en) | 1994-06-30 | 1994-06-30 | Ultrasound Doppler diagnostic device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6148885A JP3042964B2 (en) | 1994-06-30 | 1994-06-30 | Ultrasound Doppler diagnostic device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0810258A JPH0810258A (en) | 1996-01-16 |
| JP3042964B2 true JP3042964B2 (en) | 2000-05-22 |
Family
ID=15462887
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP6148885A Expired - Fee Related JP3042964B2 (en) | 1994-06-30 | 1994-06-30 | Ultrasound Doppler diagnostic device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3042964B2 (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3887040B2 (en) * | 1996-09-05 | 2007-02-28 | 株式会社東芝 | Ultrasonic diagnostic equipment |
| KR101055581B1 (en) * | 2008-03-20 | 2011-08-23 | 삼성메디슨 주식회사 | How to set filters for ultrasonic system and clutter rejection |
-
1994
- 1994-06-30 JP JP6148885A patent/JP3042964B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0810258A (en) | 1996-01-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20050054931A1 (en) | Tracking clutter filter for spectral & audio doppler | |
| JP3397748B2 (en) | Color image display method and apparatus in color Doppler image system | |
| JP7405950B2 (en) | A method for performing high spatiotemporal resolution ultrasound imaging of microvasculature | |
| JP2002534186A (en) | Ultrasonic color flow display optimization method | |
| JP2001178721A (en) | Method and device for visualizing movement in ultrasonic flow imaging using packet data acquisition | |
| US20060030777A1 (en) | T-statistic method for suppressing artifacts in blood vessel ultrasonic imaging | |
| KR100232257B1 (en) | Ultrasound color doppler imaging system for minimizing transient response to a clutter signal | |
| JPH08201500A (en) | Underwater acoustic signal detector | |
| Lovstakken et al. | Blood flow imaging-a new real-time, flow imaging technique | |
| JP3042964B2 (en) | Ultrasound Doppler diagnostic device | |
| US8189427B2 (en) | Clutter signal filtering for doppler signal | |
| JP3773779B2 (en) | Radar signal processing device | |
| CN109363722B (en) | Method and device for suppressing motion artifact in color flow imaging | |
| JP2003529422A (en) | Ultrasound system using B color priority threshold calculation | |
| JP3188394B2 (en) | Filtering method and filter circuit for ultrasonic color Doppler imaging system | |
| JP2002532172A (en) | Ultrasound method and apparatus for determining absolute radii of flank walls and arteries in tissue | |
| JP4067914B2 (en) | Ultrasonic diagnostic equipment | |
| JP2723458B2 (en) | Ultrasound Doppler diagnostic device | |
| JP2720711B2 (en) | Ultrasound two-dimensional Doppler blood flow meter | |
| JP3408647B2 (en) | Ultrasound imaging device | |
| JP2003250802A (en) | Ultrasonic diagnostic apparatus | |
| JP2005318921A (en) | Ultrasonic diagnostic equipment | |
| JP6828017B2 (en) | Spectral Doppler processing with adaptive sample window size | |
| JPH08150142A (en) | Ultrasonic diagnostic device | |
| CN117462173A (en) | Blood flow image forming apparatus and computer-readable non-transitory storage medium |
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 |