JPH08215194A - Ultrasonic wave diagnosis device - Google Patents

Ultrasonic wave diagnosis device

Info

Publication number
JPH08215194A
JPH08215194A JP3084695A JP3084695A JPH08215194A JP H08215194 A JPH08215194 A JP H08215194A JP 3084695 A JP3084695 A JP 3084695A JP 3084695 A JP3084695 A JP 3084695A JP H08215194 A JPH08215194 A JP H08215194A
Authority
JP
Japan
Prior art keywords
value
filter
dispersion value
correction
component
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
Application number
JP3084695A
Other languages
Japanese (ja)
Other versions
JP3432627B2 (en
Inventor
Takao Jibiki
隆夫 地挽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Healthcare Japan Corp
Original Assignee
GE Yokogawa Medical System Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by GE Yokogawa Medical System Ltd filed Critical GE Yokogawa Medical System Ltd
Priority to JP3084695A priority Critical patent/JP3432627B2/en
Publication of JPH08215194A publication Critical patent/JPH08215194A/en
Application granted granted Critical
Publication of JP3432627B2 publication Critical patent/JP3432627B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Abstract

PURPOSE: To prevent a flow value from being underestimated lower than its true value when a part of flow components is removed by a MTI filter. CONSTITUTION: Ultrasonic waves are transmitted to a testee' body by an ultrasonic wave probe 1, ultrasonic wave echo signals from the testee' body are received, the low frequency components of the ultrasonic wave echo signals are removed by MTI filters 9 and 10 based on filter characteristics selected in advance, and a flow dispersion value σ is thereby obtained by a CFM operation section 12 based on the ultrasonic wave echo signals from which the low frequency components are removed. Besides, the ratio α of a reference band width to the band width of selected filter characteristics is obtained, a connection factor γ is computed by a correction factor computing section 18 by letting the band width ratio α be multiplied by an adjusting factor β, and a correction dispersion value σ" is thereby computed by a correction multiplying section 20 by letting the aforesaid flow dispersion value σ' be multiplied by the correction factor γ. By this constitution, the flow dispersion value can be prevented from being underestimated to a value parted from its true value because of the characteristics of the MTI filter, so that the reliability of diagnosis can thereby be enhanced.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】この発明は、超音波診断装置に関
する。さらに詳しくは、MTIフィルタの特性に起因し
て流れの分散値が真の値から離れた値に評価されること
を防止できるようにした超音波診断装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an ultrasonic diagnostic apparatus. More specifically, the present invention relates to an ultrasonic diagnostic apparatus capable of preventing a dispersion value of a flow from being evaluated to a value apart from a true value due to the characteristics of an MTI filter.

【0002】[0002]

【従来の技術】図5は、従来の超音波診断装置の一例の
ブロック図である。この超音波診断装置500では、超
音波探触子1および送受信器2により、各音線毎に、M
発(例えば15発)の超音波パルスを時間間隔T(例え
ば1ms)で送信する。そして、各音線の深さ方向のN
点(例えば1000点)のサンプリング点に対応した超
音波エコー信号を採取する。CFM(Color Flow Mappi
ng)画像を表示する場合には、上記超音波エコー信号を
直交検波部4に入力する。直交検波部4のミキサー4
a,4bは、参照信号発生器4c,4dからの参照信号
と前記超音波エコー信号とを乗算し、LPF(Low Pa
ss Filter)4g,4hを通して、直交成分Qと同相成
分Iを出力する。A/D変換部5,6は、直交成分Qと
同相成分IをそれぞれA/D変換し、メモリ7,8に書
き込む。
2. Description of the Related Art FIG. 5 is a block diagram of an example of a conventional ultrasonic diagnostic apparatus. In this ultrasonic diagnostic apparatus 500, the ultrasonic probe 1 and the transmitter / receiver 2 are used to set M for each sound ray.
Emitting (for example, 15) ultrasonic pulses are transmitted at time intervals T (for example, 1 ms). And N in the depth direction of each sound ray
Ultrasonic echo signals corresponding to sampling points (for example, 1000 points) are collected. CFM (Color Flow Mappi)
ng) When displaying an image, the ultrasonic echo signal is input to the quadrature detection unit 4. Mixer 4 of quadrature detector 4
a and 4b multiply the reference signal from the reference signal generators 4c and 4d by the ultrasonic echo signal to obtain an LPF (Low Pa).
The quadrature component Q and the in-phase component I are output through the ss Filter) 4g and 4h. The A / D converters 5 and 6 A / D-convert the quadrature component Q and the in-phase component I, respectively, and write them in the memories 7 and 8.

【0003】MTIフィルタ9,10には、直交成分Q
と同相成分Iからクラッタ成分(比較的に動きの遅い心
臓壁などの組織からの成分)を除去するための幾つかの
フィルタ特性が記憶されている。操作者は、入力装置1
1を用いて、MTIフィルタ9,10に記憶されている
上記フィルタ特性の中から1つのフィルタ特性を選択で
きる。MTIフィルタ9,10は、メモリ7,8から直
交成分Qと同相成分Iとを呼び出し、選択されたフィル
タ特性に基づいて前記直交成分Qと同相成分Iからクラ
ッタ成分を除去する。
The MTI filters 9 and 10 have a quadrature component Q.
Several filter characteristics for removing clutter components (components from tissues such as a relatively slow moving heart wall) from the in-phase component I are stored. The operator operates the input device 1
1 can be used to select one of the filter characteristics stored in the MTI filters 9 and 10. The MTI filters 9 and 10 call the quadrature component Q and the in-phase component I from the memories 7 and 8, and remove the clutter component from the quadrature component Q and the in-phase component I based on the selected filter characteristic.

【0004】CMF演算部12は、前記不要なクラッタ
成分を除去した直交成分Qと同相成分Iとに対して、一
つの音線上の各サンプリング点n(n=1,2,…,
N)における流れの平均速度値vおよび分散値σ’をM
発の超音波パルス間の自己相関演算によって算出し、D
SC(Digital Scan Convertor)13に渡す。DSC1
3は、前記流れの平均速度値vおよび分散値σ’を表示
値(カラーおよび輝度)に変換し且つサンプリング点の
位置に応じて2次元マッピングし、CFM画像データを
生成する。なお、一般には、流れの平均速度値vは、そ
の符号により赤または青のカラーとし、絶対値の大きさ
により輝度を定める。また、分散値σ’は、緑色のカラ
ーとし、その大きさにより輝度を定める。CRT14
は、前記CFM画像データに基づいて、速度カラー画像
または分散カラー画像を画面に表示する。
The CMF calculator 12 has sampling points n (n = 1, 2, ...,) on one sound ray for the quadrature component Q and the in-phase component I from which the unnecessary clutter components have been removed.
The average velocity value v and the variance value σ ′ of the flow in N)
Calculated by the autocorrelation calculation between the emitted ultrasonic pulses, D
Pass to SC (Digital Scan Convertor) 13. DSC1
3 converts the average velocity value v and the variance value σ ′ of the flow into display values (color and brightness) and two-dimensionally maps them according to the positions of sampling points to generate CFM image data. In general, the average velocity value v of the flow is a red or blue color depending on its sign, and the brightness is determined by the magnitude of the absolute value. The variance value σ ′ is a green color, and the brightness is determined by its size. CRT14
Displays a velocity color image or a dispersed color image on the screen based on the CFM image data.

【0005】一方、Bモード画像を得る場合には、超音
波エコー信号をBモード処理部3に入力する。Bモード
処理部3は、超音波エコー信号の強さからBモードデー
タを生成し、DSC13に渡す。DSC13は、そのB
モードデータを表示値(カラーおよび輝度)に変換し且
つサンプリング点の位置に応じて2次元マッピングし、
Bモード画像データを生成する。なお、一般には、Bモ
ードデータは、グレイのカラーとし、その大きさにより
輝度を定める。CRT14は、前記Bモード画像データ
に基づいて、Bモード画像を画面に表示する。
On the other hand, when obtaining a B-mode image, an ultrasonic echo signal is input to the B-mode processing unit 3. The B-mode processing unit 3 generates B-mode data from the strength of the ultrasonic echo signal and transfers it to the DSC 13. DSC13 is the B
The mode data is converted into display values (color and brightness) and two-dimensionally mapped according to the position of the sampling point,
B-mode image data is generated. Generally, the B-mode data is a gray color, and the brightness is determined by its size. The CRT 14 displays a B-mode image on the screen based on the B-mode image data.

【0006】図6は、クラッタ成分と血流成分とMTI
フィルタ9,10のフィルタ特性の関係を示す説明図で
ある。クラッタ成分Cおよび血流成分Bの大きさは、そ
れぞれの直交成分Qと同相成分Iの2乗の和の平方根に
より算出したものである。パルス繰り返し周波数PRF
(Pulse Repetation Frequency)は、超音波パルスを送
信する時間間隔Tの逆数である。フィルタ特性Mは、そ
のカットオフ周波数fc1がクラッタ成分Cより高く且
つ血流成分Bよりも低くなるように選択される。
FIG. 6 shows clutter components, blood flow components, and MTI.
It is explanatory drawing which shows the relationship of the filter characteristics of the filters 9 and 10. The sizes of the clutter component C and the blood flow component B are calculated by the square root of the sum of squares of the orthogonal component Q and the in-phase component I. Pulse repetition frequency PRF
(Pulse Repetation Frequency) is the reciprocal of the time interval T for transmitting ultrasonic pulses. The filter characteristic M is selected so that its cutoff frequency fc1 is higher than the clutter component C and lower than the blood flow component B.

【0007】図7は、図6のフィルタ特性Mに基づいて
クラッタ成分Cが除去された後の信号成分を表わす説明
図である。血流成分Bに対応する信号成分bのみが取り
出されるため、正しい分散値σ’が得られる。
FIG. 7 is an explanatory diagram showing the signal component after the clutter component C has been removed based on the filter characteristic M of FIG. Since only the signal component b corresponding to the blood flow component B is extracted, the correct variance value σ ′ can be obtained.

【0008】[0008]

【発明が解決しようとする課題】一般に、心臓領域で
は、図8に示すように、クラッタ成分Cの周波数帯域が
広く、クラッタ成分Cと血流成分Bの周波数帯域があま
り離れていない。このため、MTIフィルタ9,10の
フィルタ特性として比較的低いカットオフ周波数fc1
のフィルタ特性Mを選択した場合は、クラッタ成分Cが
うまく除去されず、図9に示すように、クラッタ成分C
に対応する信号成分c’が残ってしまう。しかし、クラ
ッタ成分Cに対応する信号成分c’が残ると、モーショ
ンアーチファクトを生じてしまう。また、クラッタ成分
Cに対応する信号成分c’の影響で真の値より過大に分
散値σ1’が算出されてしまう。一方、図10に示すよ
うに、MTIフィルタ9,10のフィルタ特性として比
較的高いカットオフ周波数fc2のフィルタ特性Mを選
択した場合は、クラッタ成分Cが完全に除去されるた
め、モーションアーチファクトを生じない。ところが、
血流成分Bも部分的にカットされてしまうため、図11
に示すように、血流成分Bに対応する信号成分b2が小
さくなり、これの影響で真の値より過小に分散値σ2’
が算出されてしまう。そこで、この発明の目的は、MT
Iフィルタのフィルタ特性として比較的高いカットオフ
周波数のフィルタ特性を選択した場合でも、流れの分散
値が真の値に近く評価されるようにした超音波診断装置
を提供することにある。
Generally, in the heart region, as shown in FIG. 8, the frequency band of the clutter component C is wide and the frequency bands of the clutter component C and the blood flow component B are not so far apart from each other. Therefore, the cut-off frequency fc1 is relatively low as the filter characteristic of the MTI filters 9 and 10.
When the filter characteristic M of is selected, the clutter component C is not removed well, and as shown in FIG.
The signal component c'corresponding to is left. However, if the signal component c ′ corresponding to the clutter component C remains, motion artifact will occur. In addition, the variance value σ1 ′ is excessively calculated from the true value due to the influence of the signal component c ′ corresponding to the clutter component C. On the other hand, as shown in FIG. 10, when the filter characteristic M of the relatively high cutoff frequency fc2 is selected as the filter characteristic of the MTI filters 9 and 10, the clutter component C is completely removed, so that a motion artifact occurs. Absent. However,
Since the blood flow component B is also partially cut, FIG.
As shown in, the signal component b2 corresponding to the blood flow component B becomes smaller, and due to this, the variance value σ2 ′ is smaller than the true value.
Is calculated. Therefore, an object of the present invention is MT
An object of the present invention is to provide an ultrasonic diagnostic apparatus in which a flow dispersion value is evaluated close to a true value even when a filter characteristic having a relatively high cutoff frequency is selected as the filter characteristic of the I filter.

【0009】[0009]

【課題を解決するための手段】第1の観点では、この発
明は、超音波探触子から被検体内に超音波を送信し、被
検体内からの超音波エコー信号を受信し、予め選択した
MTIフィルタを通して前記超音波エコー信号の低域周
波数成分を除去し、低域周波数成分を除去した前記超音
波エコー信号に基づいて少なくとも流れの分散値を取得
して表示する超音波診断装置において、選択したMTI
フィルタのフィルタ特性を取得しそのフィルタ特性に基
づいて前記流れの分散値を補正する分散値補正手段を具
備したことを特徴とする超音波診断装置を提供する。
According to a first aspect of the present invention, an ultrasonic probe transmits an ultrasonic wave into a subject, receives an ultrasonic echo signal from the inside of the subject, and preselects the ultrasonic echo signal. In the ultrasonic diagnostic apparatus, the low frequency component of the ultrasonic echo signal is removed through the MTI filter, and at least the dispersion value of the flow is acquired and displayed based on the ultrasonic echo signal from which the low frequency component is removed. Selected MTI
There is provided an ultrasonic diagnostic apparatus comprising a dispersion value correcting means for acquiring a filter characteristic of a filter and correcting the dispersion value of the flow based on the filter characteristic.

【0010】第2の観点では、この発明は、上記構成の
超音波診断装置において、前記分散値補正手段が、選択
したMTIフィルタの帯域幅と基準帯域幅の帯域幅比を
求める帯域幅比算出手段と、前記帯域幅比に調整係数を
乗じて補正係数を算出する補正係数算出手段と、前記流
れの分散値に前記補正係数を乗じて補正分散値を算出す
る補正乗算手段とを備えたことを特徴とする超音波診断
装置を提供する。
According to a second aspect of the present invention, in the ultrasonic diagnostic apparatus having the above structure, the dispersion value correcting means calculates a bandwidth ratio for obtaining a bandwidth ratio between the bandwidth of the selected MTI filter and the reference bandwidth. Means, correction coefficient calculation means for calculating a correction coefficient by multiplying the bandwidth ratio by an adjustment coefficient, and correction multiplication means for calculating a correction dispersion value by multiplying the dispersion value of the flow by the correction coefficient. An ultrasonic diagnostic apparatus is provided.

【0011】なお、上記構成において、前記分散値補正
手段の有効化(補正する)/無効化(補正しない)を操
作者が選択するための切換手段をさらに具備することが
好ましい。
In the above structure, it is preferable to further include a switching unit for the operator to select whether the dispersion value correcting unit is enabled (corrected) or disabled (not corrected).

【0012】[0012]

【作用】上記第1の観点による超音波診断装置では、M
TIフィルタを通して超音波エコー信号の低域周波数成
分を除去した後で算出した分散値を、MTIフィルタの
フィルタ特性に基づいて補正する。例えば、MTIフィ
ルタのフィルタ特性として比較的低いカットオフ周波数
のフィルタ特性を選択すると、上述のように、クラッタ
成分を完全に除去できずに過大に分散値が算出されてし
まうことがある。そこで、比較的低いカットオフ周波数
のフィルタ特性の場合には、“1”より小さな補正係数
を乗算すれば、真の値に近い分散値を得られる。一方、
MTIフィルタのフィルタ特性として比較的高いカット
オフ周波数のフィルタ特性を選択すると、上述のよう
に、クラッタ成分を完全に除去できるが、血流成分Bも
部分的にカットされてしまうため過小に分散値が算出さ
れてしまうことがある。そこで、比較的高いカットオフ
周波数のフィルタ特性の場合には、“1”より大きな補
正係数を乗算すれば、真の値に近い分散値を得られる。
In the ultrasonic diagnostic apparatus according to the first aspect, M
The variance value calculated after removing the low frequency components of the ultrasonic echo signal through the TI filter is corrected based on the filter characteristic of the MTI filter. For example, if a filter characteristic with a relatively low cutoff frequency is selected as the filter characteristic of the MTI filter, the clutter component may not be completely removed, and the variance value may be excessively calculated, as described above. Therefore, in the case of a filter characteristic having a relatively low cutoff frequency, a dispersion value close to a true value can be obtained by multiplying by a correction coefficient smaller than "1". on the other hand,
If a filter characteristic having a relatively high cutoff frequency is selected as the filter characteristic of the MTI filter, the clutter component can be completely removed as described above, but the blood flow component B is also partially cut, so that the dispersion value is too small. May be calculated. Therefore, in the case of a filter characteristic having a relatively high cutoff frequency, a dispersion value close to a true value can be obtained by multiplying by a correction coefficient larger than "1".

【0013】上記第2の観点による超音波診断装置で
は、選択したMTIフィルタの帯域幅と基準帯域幅の帯
域幅比を求め、その帯域幅比に調整係数を乗じて補正係
数を算出し、その補正係数を流れの分散値に乗じて補正
分散値を算出する。例えば、MTIフィルタのフィルタ
特性として比較的低いカットオフ周波数fc1と比較的
高いカットオフ周波数fc2とがあるとき、比較的低い
カットオフ周波数fc1に対応する帯域幅を基準帯域幅
とする。すると、比較的低いカットオフ周波数fc1の
フィルタ特性を選択した場合、帯域幅比は“1”にな
る。ここで、調整係数を“1”より小さな値としておく
と、“1”より小さな補正係数が得られるので、これを
乗算すれば、真の値に近い分散値を得られる。一方、比
較的高いカットオフ周波数fc2のフィルタ特性を選択
した場合、帯域幅比は“1”より大きな値になる。ここ
で、この帯域幅比と乗算した積が“1”より大きな値に
なり且つ“1”より小さな値の調整係数を選んでおけ
ば、補正係数は“1”より大きくなるので、これを乗算
すれば、真の値に近い分散値を得られる。
In the ultrasonic diagnostic apparatus according to the second aspect, the bandwidth ratio between the bandwidth of the selected MTI filter and the reference bandwidth is obtained, the bandwidth ratio is multiplied by the adjustment coefficient, and the correction coefficient is calculated. The correction variance value is calculated by multiplying the flow variance value by the correction coefficient. For example, when the MTI filter has a relatively low cutoff frequency fc1 and a relatively high cutoff frequency fc2 as filter characteristics, the bandwidth corresponding to the relatively low cutoff frequency fc1 is used as the reference bandwidth. Then, when the filter characteristic of the comparatively low cutoff frequency fc1 is selected, the bandwidth ratio becomes "1". Here, if the adjustment coefficient is set to a value smaller than "1", a correction coefficient smaller than "1" is obtained. Therefore, when this is multiplied, a variance value close to a true value can be obtained. On the other hand, when the filter characteristic of the relatively high cutoff frequency fc2 is selected, the bandwidth ratio becomes a value larger than "1". Here, if the product multiplied by this bandwidth ratio becomes a value larger than "1" and an adjustment coefficient having a value smaller than "1" is selected, the correction coefficient becomes larger than "1". Then, a variance value close to the true value can be obtained.

【0014】なお、前記分散値補正手段の有効化(補正
する)/無効化(補正しない)を操作者が選択するため
の切換手段をさらに具備すれば、クラッタ成分の帯域幅
が比較的広い心臓部診断では前記分散値補正手段を有効
化し,クラッタ成分の帯域幅が比較的狭い腹部診断では
前記分散値補正手段を無効化する、といった切り換えが
可能となる。
If the switching means for the operator to select whether the dispersion value correcting means is to be enabled (corrected) or disabled (not corrected) is provided, the heart having a relatively wide clutter component bandwidth. It is possible to perform switching such that the dispersion value correction means is enabled in the partial diagnosis, and the dispersion value correction means is disabled in the abdominal diagnosis in which the clutter component bandwidth is relatively narrow.

【0015】[0015]

【実施例】以下、図に示す実施例によりこの発明をさら
に詳しく説明する。なお、これによりこの発明が限定さ
れるものではない。図1は、この発明の超音波診断装置
の一実施例のブロック図である。この超音波診断装置1
00では、超音波探触子1および送受信器2により、各
音線毎に、M発(例えば15発)の超音波パルスを時間
間隔T(例えば1ms)で送信する。そして、各音線の
深さ方向のN点(例えば1000点)のサンプリング点
に対応した超音波エコー信号を採取する。
The present invention will be described in more detail with reference to the embodiments shown in the drawings. The present invention is not limited to this. FIG. 1 is a block diagram of an embodiment of the ultrasonic diagnostic apparatus of the present invention. This ultrasonic diagnostic device 1
At 00, the ultrasonic probe 1 and the transceiver 2 transmit M (for example, 15) ultrasonic pulses at time intervals T (for example, 1 ms) for each sound ray. Then, ultrasonic echo signals corresponding to N (for example, 1000) sampling points in the depth direction of each sound ray are collected.

【0016】CFM画像を表示する場合には、上記超音
波エコー信号を直交検波部4に入力する。直交検波部4
のミキサー4a,4bは、参照信号発生器4c,4dか
らの参照信号と前記超音波エコー信号とを乗算し、LP
F4g,4hを通して、直交成分Qと同相成分Iを出力
する。A/D変換部5,6は、直交成分Qと同相成分I
をそれぞれA/D変換し、メモリ7,8に書き込む。
When displaying a CFM image, the ultrasonic echo signal is input to the quadrature detection section 4. Quadrature detector 4
Mixers 4a, 4b of the reference signal generators 4c, 4d multiply the ultrasonic echo signal by the reference signal,
The quadrature component Q and the in-phase component I are output through F4g and 4h. The A / D converters 5 and 6 have a quadrature component Q and an in-phase component I.
Are A / D converted and written in the memories 7 and 8.

【0017】MTIフィルタ9,10には、直交成分Q
と同相成分Iからクラッタ成分を除去するための4段階
のフィルタ特性が記憶されている。図2に、前記フィル
タ特性を示す。フィルタ特性M0は、最も低いカットオ
フ周波数fc0を有する。すなわち、最も広い帯域幅B
0を有する。フィルタ特性M1は、2番目に低いカット
オフ周波数fc1を有する。すなわち、2番目に広い帯
域幅B1を有する。フィルタ特性M2は、2番目に高い
カットオフ周波数fc2を有する。すなわち、2番目に
狭い帯域幅B2を有する。フィルタ特性M3は、最も高
いカットオフ周波数fc3を有する。すなわち、最も狭
い帯域幅B3を有する。図1に戻り、操作者は、入力装
置11を用いて、MTIフィルタ9,10に記憶されて
いる上記フィルタ特性M0〜M3の中から1つのフィル
タ特性Miを選択することが出来る。
The MTI filters 9 and 10 have a quadrature component Q.
And four-stage filter characteristics for removing the clutter component from the in-phase component I are stored. FIG. 2 shows the filter characteristics. The filter characteristic M0 has the lowest cutoff frequency fc0. That is, the widest bandwidth B
Has zero. The filter characteristic M1 has the second lowest cutoff frequency fc1. That is, it has the second largest bandwidth B1. The filter characteristic M2 has the second highest cutoff frequency fc2. That is, it has the second narrowest bandwidth B2. The filter characteristic M3 has the highest cutoff frequency fc3. That is, it has the narrowest bandwidth B3. Returning to FIG. 1, the operator can use the input device 11 to select one filter characteristic Mi from the filter characteristics M0 to M3 stored in the MTI filters 9 and 10.

【0018】MTIフィルタ9,10は、メモリ7,8
から直交成分Qと同相成分Iとを呼び出し、選択された
フィルタ特性Miに基づいて前記直交成分Qと同相成分
Iからクラッタ成分を除去する。CMF演算部12は、
前記不要なクラッタ成分を除去した直交成分Qと同相成
分Iとに対して、一つの音線上の各サンプリング点n
(n=1,2,…,N)における流れの平均速度値vお
よび分散値σ’をM発の超音波パルス間の自己相関演算
によって算出し、分散値補正部15の補正乗算部20に
渡す。
The MTI filters 9 and 10 include memories 7 and 8, respectively.
To call the quadrature component Q and the in-phase component I, and remove the clutter component from the quadrature component Q and the in-phase component I based on the selected filter characteristic Mi. The CMF calculator 12
For each of the quadrature component Q and the in-phase component I from which the unnecessary clutter component has been removed, each sampling point n on one sound ray
The average velocity value v and the dispersion value σ ′ of the flow at (n = 1, 2, ... hand over.

【0019】分散値補正部15は、帯域幅比算出部16
と、調整係数設定部17と、補正係数算出部18と、有
効/無効切換部19と、前記補正乗算部20とを備えて
いる。前記帯域幅比算出部16は、MTIフィルタ9,
10のフィルタ特性M0〜M1に対応した帯域幅B0〜
M3が記憶されている。操作者が入力装置11を用いて
MTIフィルタ9,10のフィルタ特性Miを選択する
と、帯域幅比算出部16は、最も広い帯域幅B0を基準
帯域幅とし、選択したMTIフィルタの帯域幅Biと基
準帯域幅の帯域幅比α=B0/Biを算出する。前記調
整係数設定部17は、操作者が入力装置11を用いて設
定した調整係数βを保持する。この調整係数βは、補正
が過剰にならないようにするためのもので、操作者が経
験的に設定する。前記補正係数算出部18は、前記帯域
幅比αに調整係数βを乗じて補正係数γ=α・βを算出
する前記有効/無効切換部19は、操作者が入力装置1
1を用いて“有効”を選択すると前記補正係数γを前記
補正乗算部20に入力し、“無効”を選択すると“1”
を前記補正乗算部20に入力する。前記補正乗算部20
は、分散値σ’と前記補正係数γまたは“1”とを乗算
し、補正した分散値σ”をDSC13に渡す。DSC1
3は、前記流れの平均速度値vおよび分散値σ”を表示
値(カラーおよび輝度)に変換し且つサンプリング点の
位置に応じて2次元マッピングし、CFM画像データを
生成する。CRT14は、前記CFM画像データに基づ
いて、速度カラー画像または分散カラー画像を画面に表
示する。
The dispersion value correction unit 15 includes a bandwidth ratio calculation unit 16
An adjustment coefficient setting unit 17, a correction coefficient calculation unit 18, a valid / invalid switching unit 19, and the correction multiplication unit 20. The bandwidth ratio calculator 16 includes the MTI filter 9,
Bandwidth B0 corresponding to 10 filter characteristics M0 to M1
M3 is stored. When the operator uses the input device 11 to select the filter characteristics Mi of the MTI filters 9 and 10, the bandwidth ratio calculation unit 16 sets the widest bandwidth B0 as the reference bandwidth and sets the bandwidth Bi of the selected MTI filter. The bandwidth ratio α = B0 / Bi of the reference bandwidth is calculated. The adjustment coefficient setting unit 17 holds the adjustment coefficient β set by the operator using the input device 11. This adjustment coefficient β is for preventing excessive correction, and is set empirically by the operator. The correction coefficient calculation unit 18 calculates the correction coefficient γ = α · β by multiplying the bandwidth ratio α by the adjustment coefficient β.
When 1 is used to select “valid”, the correction coefficient γ is input to the correction multiplication unit 20, and when “invalid” is selected, “1” is input.
Is input to the correction multiplication unit 20. The correction multiplication unit 20
Multiplies the variance value σ ′ by the correction coefficient γ or “1” and passes the corrected variance value σ ”to the DSC 13.
3 converts the average velocity value v and the variance value σ ″ of the flow into display values (color and luminance) and performs two-dimensional mapping according to the position of the sampling point to generate CFM image data. A velocity color image or a dispersed color image is displayed on the screen based on the CFM image data.

【0020】さて、図3に示すように、MTIフィルタ
9,10のフィルタ特性として比較的高いカットオフ周
波数fc2のフィルタ特性Mを選択した場合、クラッタ
成分Cが完全に除去されるため、モーションアーチファ
クトを生じない。ところが、血流成分Bも部分的にカッ
トされてしまうため、図4に示すように、血流成分Bに
対応する信号成分b2が小さくなり、これの影響で真の
値より過小に分散値σ2’が算出されてしまう。しか
し、例えば、帯域幅比α=2.0,調整係数β=0.7
とするとき、補正係数γ=α・β=1.4が分散値σ
2’に乗算されるため、真の値に近い分散値σ2”が算
出される結果となる。Bモード画像を得る場合の動作は
従来と同じであり、説明を省略する。
As shown in FIG. 3, when the filter characteristic M of the relatively high cutoff frequency fc2 is selected as the filter characteristic of the MTI filters 9 and 10, the clutter component C is completely removed, so that the motion artifact is generated. Does not occur. However, since the blood flow component B is also partially cut, the signal component b2 corresponding to the blood flow component B becomes small as shown in FIG. 4, and due to this, the variance value σ2 is too small from the true value. 'Is calculated. However, for example, the bandwidth ratio α = 2.0 and the adjustment coefficient β = 0.7
And the correction coefficient γ = α · β = 1.4 is the variance value σ
Since it is multiplied by 2 ′, the result is that the variance value σ2 ″ close to the true value is calculated. The operation for obtaining the B-mode image is the same as the conventional one, and the description thereof is omitted.

【0021】以上の超音波診断装置100によれば、心
臓部診断では前記分散値補正部15を有効化し,腹部診
断では前記分散値補正部15を無効化することで、流れ
の分散値が真の値に近く評価されるようになる。すなわ
ち、MTIフィルタの特性に起因して流れの分散値が真
の値から離れた値に評価されることを防止できる。
According to the above ultrasonic diagnostic apparatus 100, the variance value correction unit 15 is enabled in the heart part diagnosis, and the variance value correction unit 15 is disabled in the abdominal part diagnosis, so that the variance value of the flow is true. It comes to be evaluated close to the value of. That is, it is possible to prevent the variance value of the flow from being evaluated to a value apart from the true value due to the characteristics of the MTI filter.

【0022】[0022]

【発明の効果】この発明の超音波診断装置によれば、M
TIフィルタの特性に起因して流れの分散値が真の値か
ら離れた値に評価されることを防止することが出来る。
このため、診断の信頼性を向上できる。
According to the ultrasonic diagnostic apparatus of the present invention, M
It is possible to prevent the variance value of the flow from being evaluated to a value apart from the true value due to the characteristics of the TI filter.
Therefore, the reliability of diagnosis can be improved.

【図面の簡単な説明】[Brief description of drawings]

【図1】この発明の超音波診断装置の一実施例のブロッ
ク図である。
FIG. 1 is a block diagram of an embodiment of an ultrasonic diagnostic apparatus of the present invention.

【図2】フィルタ特性の説明図である。FIG. 2 is an explanatory diagram of filter characteristics.

【図3】クラッタ成分と血流成分とMTIフィルタのフ
ィルタ特性の関係を示す説明図である。
FIG. 3 is an explanatory diagram showing a relationship among clutter components, blood flow components, and filter characteristics of an MTI filter.

【図4】分散値の補正の説明図である。FIG. 4 is an explanatory diagram of correction of a variance value.

【図5】従来の超音波診断装置の一例のブロック図であ
る。
FIG. 5 is a block diagram of an example of a conventional ultrasonic diagnostic apparatus.

【図6】MTIフィルタによるクラッタ成分の除去の説
明図である。
FIG. 6 is an explanatory diagram of clutter component removal by an MTI filter.

【図7】クラッタ成分が除去された後の信号成分を表わ
す説明図である。
FIG. 7 is an explanatory diagram showing a signal component after the clutter component is removed.

【図8】比較的低いカットオフ周波数のフィルタ特性に
よるクラッタ成分の除去の説明図である。
FIG. 8 is an explanatory diagram of clutter component removal by a filter characteristic of a relatively low cutoff frequency.

【図9】クラッタ成分が不完全に除去された後の信号成
分を表わす説明図である。
FIG. 9 is an explanatory diagram illustrating a signal component after the clutter component is incompletely removed.

【図10】比較的高いカットオフ周波数のフィルタ特性
によるクラッタ成分の除去の説明図である。
FIG. 10 is an explanatory diagram of clutter component removal by a filter characteristic of a relatively high cutoff frequency.

【図11】血流成分も部分的に除去された信号成分を表
わす説明図である。
FIG. 11 is an explanatory diagram showing a signal component in which a blood flow component is also partially removed.

【符号の説明】[Explanation of symbols]

100 超音波診断装置 1 超音波探触子 9,10 MTIフィルタ 11 入力装置 15 分散値補正部 16 帯域幅比算出部 17 調整係数設定部 18 補正係数算出部 19 有効/無効切換部 20 補正乗算部 σ,σ’,σ” 分散値 α 帯域幅比 β 調整係数 γ 補正係数 100 ultrasonic diagnostic apparatus 1 ultrasonic probe 9,10 MTI filter 11 input device 15 dispersion value correction unit 16 bandwidth ratio calculation unit 17 adjustment coefficient setting unit 18 correction coefficient calculation unit 19 valid / invalid switching unit 20 correction multiplication unit σ, σ ', σ ”variance value α bandwidth ratio β adjustment coefficient γ correction coefficient

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】 超音波探触子から被検体内に超音波を送
信し、被検体内からの超音波エコー信号を受信し、予め
選択したMTI(Moving Target Indication)フィルタ
を通して前記超音波エコー信号の低域周波数成分を除去
し、低域周波数成分を除去した前記超音波エコー信号に
基づいて少なくとも流れの分散値を取得して表示する超
音波診断装置において、 選択したMTIフィルタのフィルタ特性を取得しそのフ
ィルタ特性に基づいて前記流れの分散値を補正する分散
値補正手段を具備したことを特徴とする超音波診断装
置。
1. An ultrasonic echo signal is transmitted from an ultrasonic probe into a subject, an ultrasonic echo signal is received from the subject, and the ultrasonic echo signal is passed through a preselected MTI (Moving Target Indication) filter. In the ultrasonic diagnostic apparatus that removes the low frequency component of, and acquires and displays at least the dispersion value of the flow based on the ultrasonic echo signal from which the low frequency component has been removed, the filter characteristic of the selected MTI filter is acquired. An ultrasonic diagnostic apparatus is provided with a dispersion value correcting means for correcting the dispersion value of the flow based on the filter characteristic.
【請求項2】 請求項1に記載の超音波診断装置におい
て、前記分散値補正手段が、選択したMTIフィルタの
帯域幅と基準帯域幅の帯域幅比を求める帯域幅比算出手
段と、前記帯域幅比に調整係数を乗じて補正係数を算出
する補正係数算出手段と、前記流れの分散値に前記補正
係数を乗じて補正分散値を算出する補正乗算手段とを備
えてなることを特徴とする超音波診断装置。
2. The ultrasonic diagnostic apparatus according to claim 1, wherein the dispersion value correcting means determines a bandwidth ratio between the bandwidth of the selected MTI filter and the reference bandwidth, and the bandwidth ratio calculating means. It is characterized by comprising a correction coefficient calculation means for calculating a correction coefficient by multiplying the width ratio by an adjustment coefficient, and a correction multiplication means for calculating a correction dispersion value by multiplying the dispersion value of the flow by the correction coefficient. Ultrasonic diagnostic equipment.
JP3084695A 1995-02-20 1995-02-20 Ultrasound diagnostic equipment Expired - Fee Related JP3432627B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP3084695A JP3432627B2 (en) 1995-02-20 1995-02-20 Ultrasound diagnostic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP3084695A JP3432627B2 (en) 1995-02-20 1995-02-20 Ultrasound diagnostic equipment

Publications (2)

Publication Number Publication Date
JPH08215194A true JPH08215194A (en) 1996-08-27
JP3432627B2 JP3432627B2 (en) 2003-08-04

Family

ID=12315077

Family Applications (1)

Application Number Title Priority Date Filing Date
JP3084695A Expired - Fee Related JP3432627B2 (en) 1995-02-20 1995-02-20 Ultrasound diagnostic equipment

Country Status (1)

Country Link
JP (1) JP3432627B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006255322A (en) * 2005-03-18 2006-09-28 Matsushita Electric Ind Co Ltd Ultrasonic doppler blood flowmeter
CN100342719C (en) * 2004-01-23 2007-10-10 三洋电机株式会社 Image signal processor
JP2011125690A (en) * 2009-11-18 2011-06-30 Toshiba Corp Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100342719C (en) * 2004-01-23 2007-10-10 三洋电机株式会社 Image signal processor
JP2006255322A (en) * 2005-03-18 2006-09-28 Matsushita Electric Ind Co Ltd Ultrasonic doppler blood flowmeter
JP2011125690A (en) * 2009-11-18 2011-06-30 Toshiba Corp Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program

Also Published As

Publication number Publication date
JP3432627B2 (en) 2003-08-04

Similar Documents

Publication Publication Date Title
US6520915B1 (en) Ultrasound imaging system with intrinsic doppler capability
JP3696763B2 (en) Ultrasound imaging device
US6814703B2 (en) Apparatus and method for ultrasonic diagnostic imaging using a contrast medium
US9107602B2 (en) Adaptive clutter filtering in an ultrasound system
JP2003527905A (en) Ultrasound B-mode and Doppler flow imaging
WO2011036891A1 (en) Ultrasonic diagnostic device
JPH0975348A (en) Ultrasonic diagnostic equipment
US6146331A (en) Method for improved clutter suppression for ultrasonic color doppler imaging
JP3847124B2 (en) Doppler signal processing method and apparatus and ultrasonic diagnostic apparatus
JP4481386B2 (en) Ultrasonic diagnostic equipment
JPH0810257A (en) Ultrasonic diagnostic device
JPH1075955A (en) Ultrasonic diagnostic device
JP3432627B2 (en) Ultrasound diagnostic equipment
JP6245853B2 (en) Ultrasonic diagnostic equipment
JP3451025B2 (en) Ultrasound Doppler diagnostic equipment
JP3691825B2 (en) Ultrasonic diagnostic equipment
JP2012010845A (en) Ultrasonograph
JP2807131B2 (en) Ultrasound diagnostic equipment
JP3447150B2 (en) Ultrasound imaging device
JP3353512B2 (en) Ultrasound diagnostic equipment
JPS6125536A (en) Ultrasonic diagnostic apparatus
JP2000237187A (en) Ultrasonic diagnostic device
JP3329589B2 (en) Ultrasonic imaging display method and ultrasonic diagnostic apparatus
JP3332090B2 (en) Ultrasound diagnostic equipment
JP3187169B2 (en) Ultrasound diagnostic equipment

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

LAPS Cancellation because of no payment of annual fees