JPH0572817B2 - - Google Patents
Info
- Publication number
- JPH0572817B2 JPH0572817B2 JP2007986A JP798690A JPH0572817B2 JP H0572817 B2 JPH0572817 B2 JP H0572817B2 JP 2007986 A JP2007986 A JP 2007986A JP 798690 A JP798690 A JP 798690A JP H0572817 B2 JPH0572817 B2 JP H0572817B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- code
- average
- color doppler
- analyzed
- 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
- 238000012935 Averaging Methods 0.000 claims description 73
- 238000004458 analytical method Methods 0.000 claims description 30
- 230000033001 locomotion Effects 0.000 claims description 18
- 238000001514 detection method Methods 0.000 description 14
- 230000017531 blood circulation Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 11
- 239000000523 sample Substances 0.000 description 9
- 230000015654 memory Effects 0.000 description 8
- 238000002604 ultrasonography Methods 0.000 description 7
- 238000012937 correction Methods 0.000 description 4
- 238000000034 method Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000010992 reflux Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0609—Display arrangements, e.g. colour displays
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
- G01S15/52—Discriminating between fixed and moving objects or between objects moving at different speeds
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
- G01S15/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S15/582—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse-modulated waves and based upon the Doppler effect resulting from movement of targets
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Hematology (AREA)
- Immunology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Description
〔目次〕
概 要
産業上の利用分野
従来の技術
発明が解決しようとする課題
課題を解決するための手段
作 用
実施例
発明の効果
〔概 要〕
移動体のカラー表示を行う超音波診断装置に関
し、特に血流の流速及び方向に対応した周波数出
力及びその符号出力をカラードプラ解析器によつ
て得るようにした超音波診断装置に関し、 移動体から反射した超音波受信信号に基づいて
前記移動体の速度に対応した周波ｆ、及び該周波
数解析時刻における前記移動体の移動に伴う超音
波受信信号の周波数偏移の正負を表す符号ｓ及び
超音波受信信号のパワーｐを解析するカラードプ
ラ解析器と、カラードプラ解析器によつて解析さ
れた所定期間内の各周波数解析時刻に出力される
各符号ｓの平均を算出する平均符号回路と、カラ
ードプラ解析器によつて解析された所定期間にお
ける周波数ｆの平均を算出する周波数平均回路
と、周波数平均回路で算出した平均周波数と第１
の閾値との大小を比較する周波数比較器と、カラ
ードプラ解析器によつて解析されたパワーｐと第
２の閾値との大小を比較するパワー比較器と、周
波数比較器の比較結果とパワー比較器の比較結果
の少なくとも一方が入力され、各比較結果に応じ
てカラードプラ解析器によつて解析された符号ｓ
あるいは平均符号回路で算出した平均符号を選択
して出力する符号選択回路とを備えるように構成
した。 〔産業上の利用分野〕 本発明は、移動体のカラー表示を行う超音波診
断装置に関し、特に血流の流速及び方向に対応し
た周波数出力及びその符号出力をカラードプラ解
析器によつて得るようにした超音波診断装置に関
するものである。 〔従来の技術〕 カラードプラ解析器を利用した超音波診断装置
は、人体の表示から超音波を照射して、この照射
波と反射波の周波数のずれを測定することにより
血流の速度を測定するものであり、例えば心臓内
の血流を容易に測定することが可能になる。 このような超音波診断装置で利用されるカラー
ドプラ解析器では、反射波から求めた受信信号を
照射波に対応する送信信号及びこの送信信号と位
相が90度ずれた信号とによつて直交検波して、位
相が90度ずれた２つの出力信号Ｘ，Ｙを得る。そ
して、この２つの出力信号Ｘ，Ｙ中の心臓壁から
の信号のような不用信号を除去した後このＸ，Ｙ
に基づいて自己相関から周波数Ｆを求め、この求
めた周波数Ｆの絶対値を周波数ｆ、正負を符号
ｓ、２つの出力信号Ｘ，Ｙの強さをパワーｐとし
て出力する。このようにしてカラードプラ解析器
から出力された周波数ｆは移動物体すなわち血流
の速度に対応し、符号ｓはその速度方向に対応す
る。また、パワーｐは反射波の強度に対応する。 カラードプラ解析器を利用した超音波診断装置
では、カラードプラ解析器から出力される符号ｓ
が正である場合には周波数ｆに応じた濃淡の例え
ば青色を表示し、符号ｓが負である場合には周波
数ｆに応じた濃淡の例えば赤色を表示して、血流
の分布に対応した画像表示を行う。また、カラー
ドプラ解析器から出力されたパワーｐが所定の閾
値以下のときは、反射波のSN比が悪化するため
周波数ｆの値を零に置き換えて、ノイズの影響を
除去するようにしている。 超音波診断装置の従来技術としては、特開昭62
−270139号公報が知られている。 〔発明が解決しようとする課題〕 ところで、上述した従来方式にあつては、カラ
ードプラ解析器から出力されたパワーｐが所定の
閾値以下の場合には周波数ｆを零に置き換える処
理を行つているが、その閾値を余り大きく設定す
ると測定精度の低下につながるため、この測定精
度が許容できる範囲で閾値を設定することにな
る。従つて、カラードプラ解析器から閾値付近の
パワーに対応する周波数ｆが出力された場合には
わずかではあるがノイズが含まれることになり、
その影響を完全に除去することができない。 このようにしてノイズが含まれた場合には、周
波数ｆの値及び符号ｓはランダムな値をとるため
局部的に符号ｓの値が反転することが生じる。そ
のため、血流が一定方向に流れて例えば青色の表
示がなされているときに、所々赤色が混じるとい
つた画像の乱れが発生するという問題点があつ
た。 このような画像の乱れが発生すると、病院等に
おいて表示画像を見て診断を行う場合に、誤診断
につながることになる。 また、ノイズによつて符号ｓが反転することに
よる画像の乱れを除去する方法としては、所定領
域内の符号ｓの平均を求めて局所的な符号ｓの変
動を均一化する補正を行う方法が考えられる。し
かし、一般には、低速の血流には逆流は存在しな
いが、乱流と称される高速の血流には局部的に逆
流が発生するという事実が確認されている。この
ため、上述したように符号ｓの平均を求めて画質
を補正する場合には、乱流が発生しているときの
局部的な逆流も均一化してしまうため、このよう
な補正を適用することはできない。 本発明は、このような点にかんがみて創作され
たものであり、ノイズによる画像の乱れを防止す
ることができる超音波診断装置を提供することを
目的としている。 〔課題を解決するための手段〕 (i) 請求項１の発明 請求項１の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、及び該周波数解析時刻における前記移動体の
移動に伴う超音波受信信号の周波数偏移の正負を
表す符号ｓを解析して出力する。 平均符号回路は、平均カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と所定の閾値との大小を比較する。 符号選択回路は、前記周波数比較器の比較結果
が入力され、前記周波数平均回路で算出した平均
周波数が前記所定の閾値より大きい場合に前記カ
ラードプラ解析器によつて解析された符号ｓを出
力し、前記周波数平均回路で算出した平均周波数
が前記所定の閾値より小さい場合に前記平均符号
回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記周
波数平均回路で算出した平均周波数とに基づいて
前記移動体の速度分布を表示するように構成され
ている。 (ii) 請求項２の発明 請求項２の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、及び該周波数解析時刻における前記移動体の
移動に伴う超音波受信信号の周波数偏移の正負を
表す符号ｓを解析して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と所定の閾値との大小を比較する。 符号選択回路は、前記周波数比較器の比較結果
が入力され、前記周波数平均回路で算出した平均
周波数が前記所定の閾値より大きい場合に前記カ
ラードプラ解析器によつて解析された符号ｓを出
力し、前記周波数平均回路で算出した平均周波数
が前記所定の閾値より小さい場合に前記平均符号
回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記カ
ラードプラ解析器によつて解析された周波数ｆと
に基づいて前記移動体の速度分布を表示するよう
に構成されている。 (iii) 請求項３の発明 請求項３の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと所定の閾値との大小を
比較する。 符号選択回路は、前記パワー比較器の比較結果
が入力され、前記カラードプラ解析器で解析され
たパワーｐが前記所定の闘値より大きい場合に前
記カラードプラ解析器によつて解析された符号ｓ
を出力し、前記カラードプラ解析器で解析された
パワーｐが前記所定の閾値より小さい場合に前記
平均符号回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記周
波数平均回路で算出した平均周波数とに基づいて
前記移動体の速度分布を表示するように構成され
ている。 (iv) 請求項４の発明 請求項４の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと所定の閾値との大小を
比較する。 符号選択回路は、前記パワー比較器の比較結果
が入力され、前記カラードプラ解析器で解析され
たパワーｐが前記所定の閾値より大きい場合に前
記カラードプラ解析器によつて解析された符号ｓ
を出力し、前記カラードプラ解析器で解析された
パワーｐが前記所定の閾値より小さい場合に前記
平均符号回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記カ
ラードプラ解析器によつて解析された周波数ｆと
に基づいて前記移動体の速度分布を表示するよう
に構成されている。 (v) 請求項５の発明 請求項５の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と第１の閾値との大小を比較する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと第２の閾値との大小を
比較する。 符号選択回路は、前記周波数比較器の比較結果
と前記パワー比較器の比較結果とが入力され、前
記周波数平均回路で算出された平均周波数が前記
第１の閾値より大きく、かつ前記カラードプラ解
析器で解析されたパワーｐが前記第２の閾値より
大きい場合に前記カラードプラ解析器によつて解
析された符号ｓを出力し、それ以外の場合に前記
平均符号回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記周
波数平均回路で算出した平均周波数とに基づいて
前記移動体の速度分布を表示するように構成され
ている。 (vi) 請求項６の発明 請求項６の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と第１の閾値との大小を比較する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと第２の閾値との大小を
比較する。 符号選択回路は、前記周波数比較器の比較結果
と前記パワー比較器の比較結果とが入力され、前
記周波数平均回路で算出された平均周波数が前記
第１の閾値より大きく、かつ前記カラードプラ解
析器で解析されたパワーｐが前記第２の閾値より
大きい場合に前記カラードプラ解析器によつて解
析された符号ｓを出力し、それ以外の場合に前記
平均符号回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記カ
ラードプラ解析器によつて解析された周波数ｆと
に基づいて前記移動体の速度分布を表示するよう
に構成されている。 (vii) 請求項７の発明 請求項７の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と第１の閾値との大小を比較する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと第２の閾値との大小を
比較する。 符号選択回路は、前記周波数比較器の比較結果
と前記パワー比較器の比較結果とが入力され、前
記周波数平均回路で算出された平均周波数が前記
第１の閾値より小さく、かつ前記カラードプラ解
析器で解析されたパワーｐが前記第２の閾値より
小さい場合に前記平均符号回路で算出された平均
符号を出力し、それ以外の場合に前記カラードプ
ラ解析器で解析された符号ｓを出力する。 全体として、前記符号選択回路の出力と前記周
波数平均回路で算出した平均周波数とに基づいて
前記移動体の速度分布を表示するように構成され
ている。 (viii) 請求項８の発明 請求項８の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と第１の閾値との大小を比較する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと第２の閾値との大小を
比較する。 符号選択回路は、前記周波数比較器の比較結果
と前記パワー比較器の比較結果とが入力され、前
記周波数平均回路で算出された平均周波数が前記
第１の閾値より小さく、かつ前記カラードプラ解
析器で解析されたパワーｐが前記第２の閾値より
小さい場合に前記平均符号回路で算出された平均
符号を出力し、それ以外の場合に前記カラードプ
ラ解析器で解析された符号ｓを出力する。 全体として、前記符号選択回路の出力と前記カ
ラードプラ解析器によつて解析された周波数ｆと
に基づいて前記移動体の速度分布を表示するよう
に構成されている。 〔作 用〕 (i) 請求項１の発明 周波数平均回路によつて算出された平均周波数
が所定の閾値より大きい場合には、カラードプラ
解析器によつて解析された符号ｓを選択し、この
選択した符号ｓと周波数平均回路で算出した平均
周波数とに基づいて移動体の速度分布の表示を行
う。 一方、平均周波数が所定の閾値より小さい場合
には、符号平均回路によつて算出した平均符号を
選択し、この選択した平均符号と周波数平均回路
で算出した平均周波数とに基づいて移動体の速度
分布の表示を行う。 請求項１の発明にあつては、移動体の速度が遅
く符号ｓに含まれるノイズの割合が多い場合に、
符号ｓに代えて平均符号が選択される。 (ii) 請求項２の発明 周波数平均回路によつて算出された平均周波数
が所定の閾値より大きい場合には、カラードプラ
解析器によつて解析された符号ｓを選択し、この
選択した符号ｓとカラードプラ解析器によつて解
析された周波数ｆとに基づいて移動体の速度分布
の表示を行う。 一方、平均周波数が所定の閾値より小さい場合
には、符号平均回路によつて算出した平均符号を
選択し、この選択した平均符号とカラードプラ解
析器によつて解析された周波数ｆとに基づいて移
動体の速度分布の表示を行う。 請求項２の発明にあつては、移動体の速度が遅
く符号ｓに含まれるノイズの割合が多い場合に、
符号ｓに代えて平均符号が選択される。 (iii) 請求項３の発明 カラードプラ解析器によつて解析されたパワー
ｐが所定の閾値より大きい場合には、カラードプ
ラ解析器によつて解析された符号ｓを選択し、こ
の選択した符号ｓと周波数平均回路で算出した平
均周波数とに基づいて移動体の速度分布の表示を
行う。 一方、パワーｐが所定の閾値より小さい場合に
は、符号平均回路によつて算出した平均符号を選
択し、この選択した平均符号と周波数平均回路で
算出した平均周波数とに基づいて移動体の速度分
布の表示を行う。 請求項３の発明にあつては、反射してくる超音
波受信信号の強度が弱く符号ｓに含まれるノイズ
の割合が多い場合に、符号ｓに代えて平均符号が
選択される。 (iv) 請求項４の発明 カラードプラ解析器によつて解析されたパワー
ｐが所定の閾値より大きい場合には、カラードプ
ラ解析器によつて解析された符号ｓを選択し、こ
の選択した符号ｓとカラードプラ解析器によつて
解析された周波数ｆとに基づいて移動体の速度分
布の表示を行う。 一方、パワーｐが所定の閾値より小さい場合に
は、符号平均回路によつて算出した平均符号を選
択し、この選択した平均符号とカラードプラ解析
器によつて解析された周波数ｆとに基づいて移動
体の速度分布の表示を行う。 請求項４の発明にあつては、反射してくる超音
波受信信号の強度が弱く符号ｓに含まれるノイズ
の割合が多い場合に、符号ｓに代えて平均符号が
選択される。 (v) 請求項５の発明 周波数平均回路によつて算出された平均周波数
が第１の閾値より大きく、かつカラードプラ解析
器によつて解析されたパワーｐが第２の閾値より
大きい場合には、カラードプラ解析器によつて解
析された符号ｓを選択し、この選択した符号ｓと
周波数平均回路で算出した平均周波数とに基づい
て移動体の速度分布の表示を行う。 一方、それ以外の場合には、符号平均回路によ
つて算出した平均符号を選択し、この選択した平
均符号と周波数平均回路で算出した平均周波数と
に基づいて移動体の速度分布の表示を行う。 請求項５の発明にあつては、移動体の速度が遅
く、かつ反射してくる超音波受信信号の強度が弱
く符号ｓに含まれるノイズの割合が多い場合に、
符号ｓに代えて平均符号が選択される。 (vi) 請求項６の発明 周波数平均回路によつて算出された平均周波数
が第１の閾値より大きく、かつカラードプラ解析
器によつて解析されたパワーｐが第２の閾値より
大きい場合には、カラードプラ解析器によつて解
析された符号ｓを選択し、この選択した符号ｓと
カラードプラ解析器によつて解析された周波数ｆ
とに基づいて移動体の速度分布の表示を行なう。 一方、それ以外の場合には、符号平均回路によ
つて算出した平均符号を選択し、この選択した平
均符号とカラードプラ解析器によつて解析された
周波数ｆとに基づいて移動体の速度分布の表示を
行う。 請求項６の発明にあつては、移動体の速度が遅
く、かつ反射してくる超音波受信信号の強度が弱
く符号ｓに含まれるノイズの割合が多い場合に、
符号ｓに代えて平均符号が選択される。 (vii) 請求項７の発明 周波数平均回路によつて算出された平均周波数
が第１の閾値より小さく、かつカラードプラ解析
器によつて解析されたパワーｐが第２の閾値より
小さい場合には、符号平均回路によつて算出した
平均符号を選択し、この選択した平均符号と周波
数平均回路で算出した平均周波数とに基づいて移
動体の速度分布の表示を行う。 一方、それ以外の場合には、カラードプラ解析
器によつて解析された符号ｓを選択し、この選択
した符号ｓと周波数平均回路で算出した平均周波
数とに基づいて移動体の速度分布の表示を行う。 請求項７の発明にあつては、移動体の速度が遅
いか、あるいは反射してくる超音波受信信号の強
度が弱く符号ｓに含まれるノイズの割合が多い場
合に、符号ｓに代えて平均符号が選択される。 (viii) 請求項８の発明 周波数平均回路によつて算出された平均周波数
が第１の閾値より小さく、かつカラードプラ解析
器によつて解析されたパワーｐが第２の閾値より
小さい場合には、符号平均回路によつて算出した
平均符号を選択し、この選択した平均符号とカラ
ードプラ解析器によつて解析された周波数ｆとに
基づいて移動体の速度分布の表示を行う。 一方、それ以外の場合には、カラードプラ解析
器によつて解析された符号ｓを選択し、この選択
した符号ｓとカラードプラ解析器によつて解析さ
れた周波数ｆとに基づいて移動体の速度分布の表
示を行う。 請求項８の発明にあつては、移動体の速度が遅
いか、あるいは反射してくる超音波受信信号の強
度が弱く符号ｓに含まれるノイズの割合が多い場
合に、符号ｓに代えて平均符号が選択される。 〔実施例〕 以下、図面に基づいて本発明の実施例について
詳細に説明する。 第１図に、本発明の超音波診断装置の一実施例
の構成を示す。 第１図において、１１１は超音波プローブ、１
３は送信回路、１１５は受信回路、１２１はカラ
ードプラ解析器、１２３は雑音除去回路、１２５
及び１４３はデイジタル・スキヤン・コンバー
タ、１２７及び１４５はメモリ、１２９及び１４
７はルツク・アツプ・テーブル、１３１は表示装
置、１４１は振幅検出回路、１４９は合成回路で
ある。 超音波プローブ１１１は、例えばPhased
Array型と称される構成を有しており、数十個程
度の超音波振動子を備えている。この超音波プロ
ーブ１１１は、送信回路１１３及び受信回路１１
５に接続されている。 送信回路１１３は、超音波プローブ１１１をド
ライブするためのものであり、超音波プローブ１
１１に備わつた超音波振動子を励振することによ
り、所定の照射角をなす超音波を出力させる。ま
た、受信回路１１５は、被測定物（移動体）によ
つて反射された超音波が超音波プローブ１１１に
入射されたときに、この反射波を検出するための
ものであり、上述した超音波振動子を順次走査す
ることにより反射波の検出を行つている。 カラードプラ解析器１２１は、受信回路１１５
の受信出力に基づいて被測定物体の速度に関する
解析を行うものである。解析の結果、カラードプ
ラ解析器１２１は、被測定物体の速度に対応した
周波数出力ｆと、この速度の方向に対応した符号
ｓと、反射波の強度に対応したパワーｐとを出力
する。カラードプラ解析器の出力は雑音除去回路
１２３に入力される。 雑音除去回路１２３は、予め設定された所定の
モードに従つて、カラードプラ解析器１２１から
出力された符号ｓに補正を加えるものである。こ
の補正は、ノイズによつて符号が反転する現象を
除去するための処理であり、ノイズによつて符号
が反転する可能性がある状態を検出したときに、
符号ｓに代えて所定時間の平均値を出力する。
この符号平均については、表示する画素の前後例
えば連続する15画素の符号の平均をとり、その15
画素に含まれる１つの画素の符号を決定する。雑
音除去回路１２３から出力される符号ｓあるいは
その平均値と周波数ｆはデイジタル・スキヤ
ン・コンバータ１２５に入力される。 デイジタル・スキヤン・コンバータ１２５は、
雑音除去回路１２３の出力をスキヤンして順次取
り込むためのものであり、この取り込んだデータ
はメモリ１２７に入力される。このデイジタル・
スキヤン・コンバータ１２５によるスキヤン動作
は、超音波プローブ１１１から反射波を取り込む
スキヤン動作に同期している。メモリ１２７で
は、デイジタル・スキヤン・コンバータ１２５か
ら出力されるデータを、被測定物体の位置に応じ
たアドレスに順次格納する。このメモリ１２７の
格納データは順次読み出されてルツク・アツプ・
テーブル１２９に入力される。 ルツク・アツプ・テーブル１２９は、メモリ１
２７から読み出されたデータに基づいて、カラー
データ（例えばRGBデータ）を作成するための
ものである。メモリ１２７には、雑音除去回路１
２３から出力された周波数ｆ及び符号ｓあるいは
その平均値が格納されており、これらの格納デ
ータに基づいて被測定物体の速度に関するカラー
データに変換する。例えば、符号ｓあるいはその
平均値が正である場合が赤色に対応しており、
この赤色の濃淡の度合いを周波数ｆに対応させた
カラーデータへの変換を行う。反対に、符号ｓあ
るいはその平均値が負である場合が青色に対応
しており、この青色の濃淡の度合いを周波数ｆに
対応させたカラーデータへの変換が行われる。こ
の変換後のカラーデータは表示装置１３１に入力
される。 表示装置１３１は、ルツク・アツプ・テーブル
１２９から出力されるカラーデータに基づいて、
被測定物体の速度に関するカラー表示を行う。例
えば、被測定物体として血液を考えた場合、流れ
が一定方向の血流は赤色あるいは青色何れか一方
に対応して表示され、乱流と称される高速流の場
合は所々流れの方向が反転するので、青色の所々
に赤色が混じて表示されたり、赤色の所々に青色
が混じつて表示されたりする。 上述した表示装置１３１による表示は血液の流
れを表示するような場合に用いられるが、超音波
の反射波に基づいて被測定物体の形状を測定する
こともできる。以下、形状を測定する動作につい
て説明する。 振幅検出回路１４１は、受信回路１１５の受信
出力の振幅検波を行うものである。振幅検波回路
１４１は、受信信号の振幅に応じた出力、すなわ
ち超音波プローブ１１１に到達する反射波の強度
に応じた検波出力が得られる。この検波出力は、
デイジタル・スキヤン・コンバータ１４３でスキ
ヤンされてメモリ１４５に格納され、その後ルツ
ク・アツプ・テーブル１４７に入力される。ルツ
ク・アツプ・テーブル１４７は、メモリ１４５か
ら読み出されるデータに基づいて白黒データ
（RGBが同一値）を作成するためのものである。
例えば、心臓の断面に超音波を照射した場合に
は、心臓の壁面部分と内部の血液が流れている部
分とでは反射波の強度が異なるため振幅検波回路
１４１の検波出力にも差異が生じ、この差異が被
測定物体の形状として表示装置１３１に表示され
る。合成回路１４９はルツク・アツプ・テーブル
１２９，１４７の出力を合成して表示装置１３１
に供給している。 次に、上述した本実施例の詳細構成及び動作を
第２図乃至第１０図を参照して説明する。 第２図に、カラードプラ解析器１２１の詳細構
成を示す。２１１は直交検波回路、２１３はアナ
ログ−デイジタル（Ａ／Ｄ）変換器、２１５はモ
ーシヨン・ターゲツト・インジケータ（MTI）
フイルタ、２１７は自己相関器である。 受信回路１１５の出力が直交検波回路２１１に
入力されると、直交検波回路２１１は、照射波の
出力周波数fcと同一周波数の信号と、出力周波数
fcと90度位相がずれた周波数の信号とによつて入
力に対する直交検波を行う。この直交検波の結
果、位相が90度ずれた２つの出力が得られ、Ａ／
Ｄ変換器２１３に入力される。 Ａ／Ｄ変換器２１３は、入力された２つの信号
のそれぞれをデイジタルデータに変換するもので
あり、変換されたデイジタルデータはMTIフイ
ルタ２１５に入力される。 MTIフイルタ２１５は、Ａ／Ｄ変換器２１３
から入力される２つのデイジタルデータのそれぞ
れに対して、ハイパスフイルタとして動作する。
例えば、超音波プローブ１１１によつて被測定物
体の同一位置に対する反射波を繰り返し受信した
場合に、この同一位置の反射波に対応するＡ／Ｄ
変換器２１３の出力データに対して、１次のハイ
パスフイルタとして動作する。従つて、被測定物
体が移動した場合には、この移動速度に応じて
MTIフイルタ２１５の出力値が大きくなる。
MTIフイルタ２１５の２つの出力は、自己相関
器２１７に入力される。 自己相関器２１７は、MTIフイルタ２１５の
２つの出力に基づいて自己相関から周波数を求め
る。自己相関から得られる周波数をＦ、MTIフ
イルタ２１５から出力される２つの出力をＸ(i)、
Ｙ(i)（ｉは被測定物体の同一位置に対して繰り返
される超音波の送受信回数）とすると、自己相関
から得られる周波数Ｆは、 Ｆ＝１／Ｔtan-1｛Σ（Ｘ(i)Ｙ（ｉ−１
）−Ｙ(i)Ｘ（ｉ−１））／Σ（Ｘ(i)Ｘ（ｉ−１）＋Ｙ
(i)Ｙ（ｉ−１））｝…(1) で示される。ここで、Ｔは被測定物体の同一位置
に対して繰り返される超音波の送信あるいは受信
の周期である。 自己相関器２１７は、(1)式に従つて求めた自己
相関による周波数Ｆの絶対値｜Ｆ｜を周波数ｆと
して、自己相関による周波数Ｆの正負を符号ｓと
して出力する。また、自己相関器２１７は、被測
定物体の同一位置に対応する２入力Ｘ(i)、Ｙ(i)の
２乗和を累積した値Σ（Ｘ(i)2＋Ｙ(i)2）をパワー
ｐとして出力する。このパワーｐは、反射波の強
度に対応している。自己相関器２１７のこれらの
出力は、上述した雑音除去回路１２３に入力され
る。 第３図に、雑音除去回路１２３の詳細構成を示
す。第３図において、３１１はパワー比較器、３
１３は周波数選択回路、３１５は周波数平均回
路、３１７は周波数比較、３２１は平均符号回
路、３２３は等価遅延回路、３２５は符号選択回
路、３３１はモード制御回路である。 パワー比較器３１１は、カラードプラ解析器１
２１内の自己相関器２１７から出力されたパワー
ｐが入力され、このパワーｐと所定の閾値p0との
大小を比較するものである。比較結果は、周波数
選択回路３１３及びモード制御回路３３１に入力
される。 周波数選択回路３１３は、カラードプラ解析器
１２１内の自己相関器２１７から出力された周波
数ｆが入力される。周波数選択回路３１３は、パ
ワー比較器３１１の比較結果に応じて、この周波
数ｆをそのまま出力するか、あるいは零に置き換
えて出力するかの択一的動作を行つている。具体
的には、パワー比較器３１１に入力されたパワー
ｐが所定の閾値p0以下のときは、周波数ｆはノイ
ズの影響を大きく受けているものであるため、こ
の周波数ｆを零に置き換えて出力する。それ以外
のときは周波数ｆをそのまま出力する。周波数選
択回路３１３の出力は、周波数平均回路３１５に
入力される。 周波数平均回路３１５は、周波数選択回路３１
３の出力に対して、所定の時間内の移動平均を求
めて平均周波数を出力するためのものである。
周波数ｆから平均周波数を求めることにより、
画像の相関を利用した補正を行つて、後段の表示
装置１３１による表示を見易いものにしている。
周波数平均回路３１５の出力は、周波数比較器３
１７に入力されると共に、雑音除去回路１２３の
外部に出力される。 周波数比較器３１７は、入力された平均周波数
ｆと所定の閾値f0との大小を比較するものであ
る。比較結果は、モード制御回路３３１に入力さ
れる。 モード制御回路３３１は、パワー比較器３１１
及び周波数比較器３１７の各出力の少なくとも一
方に基づいて、符号選択回路３２５による選択動
作を制御するものである。モード制御回路３３１
に入力されるモード選択信号に応じて制御モード
の設定が行われ、設定された制御モードに従つた
制御信号が符号選択回路３２５に入力される。 平均符号回路３２１及び等価遅延回路３２３に
は、カラードプラ解析器１２１内の自己相関器２
１７から出力された符号ｓが入力される。平均符
号回路３２１は、この符号ｓを所定の時間内で平
均した平均符号を算出して出力する。この平均
符号は、表示する画素に対応する符号につき、例
えば連続する７画素の平均をとり、平均符号と
して出力するものである。平均符号回路３２１か
ら出力された平均符号は符号選択回路３２５に
入力される。また、等価遅延回路３２３は、平均
符号回路３２１による平均符号算出に要する時間
と同期をとるためのものであり、入力された符号
ｓを所定時間遅延した後に出力する。等価遅延回
路３２３で遅延された符号ｓは符号選択回路３２
５に入力される。 符号選択回路３２５は、平均符号回路３２１か
ら入力される平均符号または等価遅延回路３２
３から出力される符号ｓの何れか一方を選択して
出力する。この選択動作は、モード制御回路３３
１から入力される制御信号によつて制御される。
符号選択回路３２５から出力される符号ｓあるい
は平均符号は、雑音除去回路１２３の外部に出
力される。 第１表に、モード制御回路３３１において制御
される制御モードの詳細内容を示す。以下、各制
御モードについて説明する。
し、特に血流の流速及び方向に対応した周波数出
力及びその符号出力をカラードプラ解析器によつ
て得るようにした超音波診断装置に関し、 移動体から反射した超音波受信信号に基づいて
前記移動体の速度に対応した周波ｆ、及び該周波
数解析時刻における前記移動体の移動に伴う超音
波受信信号の周波数偏移の正負を表す符号ｓ及び
超音波受信信号のパワーｐを解析するカラードプ
ラ解析器と、カラードプラ解析器によつて解析さ
れた所定期間内の各周波数解析時刻に出力される
各符号ｓの平均を算出する平均符号回路と、カラ
ードプラ解析器によつて解析された所定期間にお
ける周波数ｆの平均を算出する周波数平均回路
と、周波数平均回路で算出した平均周波数と第１
の閾値との大小を比較する周波数比較器と、カラ
ードプラ解析器によつて解析されたパワーｐと第
２の閾値との大小を比較するパワー比較器と、周
波数比較器の比較結果とパワー比較器の比較結果
の少なくとも一方が入力され、各比較結果に応じ
てカラードプラ解析器によつて解析された符号ｓ
あるいは平均符号回路で算出した平均符号を選択
して出力する符号選択回路とを備えるように構成
した。 〔産業上の利用分野〕 本発明は、移動体のカラー表示を行う超音波診
断装置に関し、特に血流の流速及び方向に対応し
た周波数出力及びその符号出力をカラードプラ解
析器によつて得るようにした超音波診断装置に関
するものである。 〔従来の技術〕 カラードプラ解析器を利用した超音波診断装置
は、人体の表示から超音波を照射して、この照射
波と反射波の周波数のずれを測定することにより
血流の速度を測定するものであり、例えば心臓内
の血流を容易に測定することが可能になる。 このような超音波診断装置で利用されるカラー
ドプラ解析器では、反射波から求めた受信信号を
照射波に対応する送信信号及びこの送信信号と位
相が90度ずれた信号とによつて直交検波して、位
相が90度ずれた２つの出力信号Ｘ，Ｙを得る。そ
して、この２つの出力信号Ｘ，Ｙ中の心臓壁から
の信号のような不用信号を除去した後このＸ，Ｙ
に基づいて自己相関から周波数Ｆを求め、この求
めた周波数Ｆの絶対値を周波数ｆ、正負を符号
ｓ、２つの出力信号Ｘ，Ｙの強さをパワーｐとし
て出力する。このようにしてカラードプラ解析器
から出力された周波数ｆは移動物体すなわち血流
の速度に対応し、符号ｓはその速度方向に対応す
る。また、パワーｐは反射波の強度に対応する。 カラードプラ解析器を利用した超音波診断装置
では、カラードプラ解析器から出力される符号ｓ
が正である場合には周波数ｆに応じた濃淡の例え
ば青色を表示し、符号ｓが負である場合には周波
数ｆに応じた濃淡の例えば赤色を表示して、血流
の分布に対応した画像表示を行う。また、カラー
ドプラ解析器から出力されたパワーｐが所定の閾
値以下のときは、反射波のSN比が悪化するため
周波数ｆの値を零に置き換えて、ノイズの影響を
除去するようにしている。 超音波診断装置の従来技術としては、特開昭62
−270139号公報が知られている。 〔発明が解決しようとする課題〕 ところで、上述した従来方式にあつては、カラ
ードプラ解析器から出力されたパワーｐが所定の
閾値以下の場合には周波数ｆを零に置き換える処
理を行つているが、その閾値を余り大きく設定す
ると測定精度の低下につながるため、この測定精
度が許容できる範囲で閾値を設定することにな
る。従つて、カラードプラ解析器から閾値付近の
パワーに対応する周波数ｆが出力された場合には
わずかではあるがノイズが含まれることになり、
その影響を完全に除去することができない。 このようにしてノイズが含まれた場合には、周
波数ｆの値及び符号ｓはランダムな値をとるため
局部的に符号ｓの値が反転することが生じる。そ
のため、血流が一定方向に流れて例えば青色の表
示がなされているときに、所々赤色が混じるとい
つた画像の乱れが発生するという問題点があつ
た。 このような画像の乱れが発生すると、病院等に
おいて表示画像を見て診断を行う場合に、誤診断
につながることになる。 また、ノイズによつて符号ｓが反転することに
よる画像の乱れを除去する方法としては、所定領
域内の符号ｓの平均を求めて局所的な符号ｓの変
動を均一化する補正を行う方法が考えられる。し
かし、一般には、低速の血流には逆流は存在しな
いが、乱流と称される高速の血流には局部的に逆
流が発生するという事実が確認されている。この
ため、上述したように符号ｓの平均を求めて画質
を補正する場合には、乱流が発生しているときの
局部的な逆流も均一化してしまうため、このよう
な補正を適用することはできない。 本発明は、このような点にかんがみて創作され
たものであり、ノイズによる画像の乱れを防止す
ることができる超音波診断装置を提供することを
目的としている。 〔課題を解決するための手段〕 (i) 請求項１の発明 請求項１の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、及び該周波数解析時刻における前記移動体の
移動に伴う超音波受信信号の周波数偏移の正負を
表す符号ｓを解析して出力する。 平均符号回路は、平均カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と所定の閾値との大小を比較する。 符号選択回路は、前記周波数比較器の比較結果
が入力され、前記周波数平均回路で算出した平均
周波数が前記所定の閾値より大きい場合に前記カ
ラードプラ解析器によつて解析された符号ｓを出
力し、前記周波数平均回路で算出した平均周波数
が前記所定の閾値より小さい場合に前記平均符号
回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記周
波数平均回路で算出した平均周波数とに基づいて
前記移動体の速度分布を表示するように構成され
ている。 (ii) 請求項２の発明 請求項２の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、及び該周波数解析時刻における前記移動体の
移動に伴う超音波受信信号の周波数偏移の正負を
表す符号ｓを解析して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と所定の閾値との大小を比較する。 符号選択回路は、前記周波数比較器の比較結果
が入力され、前記周波数平均回路で算出した平均
周波数が前記所定の閾値より大きい場合に前記カ
ラードプラ解析器によつて解析された符号ｓを出
力し、前記周波数平均回路で算出した平均周波数
が前記所定の閾値より小さい場合に前記平均符号
回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記カ
ラードプラ解析器によつて解析された周波数ｆと
に基づいて前記移動体の速度分布を表示するよう
に構成されている。 (iii) 請求項３の発明 請求項３の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと所定の閾値との大小を
比較する。 符号選択回路は、前記パワー比較器の比較結果
が入力され、前記カラードプラ解析器で解析され
たパワーｐが前記所定の闘値より大きい場合に前
記カラードプラ解析器によつて解析された符号ｓ
を出力し、前記カラードプラ解析器で解析された
パワーｐが前記所定の閾値より小さい場合に前記
平均符号回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記周
波数平均回路で算出した平均周波数とに基づいて
前記移動体の速度分布を表示するように構成され
ている。 (iv) 請求項４の発明 請求項４の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと所定の閾値との大小を
比較する。 符号選択回路は、前記パワー比較器の比較結果
が入力され、前記カラードプラ解析器で解析され
たパワーｐが前記所定の閾値より大きい場合に前
記カラードプラ解析器によつて解析された符号ｓ
を出力し、前記カラードプラ解析器で解析された
パワーｐが前記所定の閾値より小さい場合に前記
平均符号回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記カ
ラードプラ解析器によつて解析された周波数ｆと
に基づいて前記移動体の速度分布を表示するよう
に構成されている。 (v) 請求項５の発明 請求項５の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と第１の閾値との大小を比較する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと第２の閾値との大小を
比較する。 符号選択回路は、前記周波数比較器の比較結果
と前記パワー比較器の比較結果とが入力され、前
記周波数平均回路で算出された平均周波数が前記
第１の閾値より大きく、かつ前記カラードプラ解
析器で解析されたパワーｐが前記第２の閾値より
大きい場合に前記カラードプラ解析器によつて解
析された符号ｓを出力し、それ以外の場合に前記
平均符号回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記周
波数平均回路で算出した平均周波数とに基づいて
前記移動体の速度分布を表示するように構成され
ている。 (vi) 請求項６の発明 請求項６の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と第１の閾値との大小を比較する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと第２の閾値との大小を
比較する。 符号選択回路は、前記周波数比較器の比較結果
と前記パワー比較器の比較結果とが入力され、前
記周波数平均回路で算出された平均周波数が前記
第１の閾値より大きく、かつ前記カラードプラ解
析器で解析されたパワーｐが前記第２の閾値より
大きい場合に前記カラードプラ解析器によつて解
析された符号ｓを出力し、それ以外の場合に前記
平均符号回路で算出した平均符号を出力する。 全体として、前記符号選択回路の出力と前記カ
ラードプラ解析器によつて解析された周波数ｆと
に基づいて前記移動体の速度分布を表示するよう
に構成されている。 (vii) 請求項７の発明 請求項７の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と第１の閾値との大小を比較する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと第２の閾値との大小を
比較する。 符号選択回路は、前記周波数比較器の比較結果
と前記パワー比較器の比較結果とが入力され、前
記周波数平均回路で算出された平均周波数が前記
第１の閾値より小さく、かつ前記カラードプラ解
析器で解析されたパワーｐが前記第２の閾値より
小さい場合に前記平均符号回路で算出された平均
符号を出力し、それ以外の場合に前記カラードプ
ラ解析器で解析された符号ｓを出力する。 全体として、前記符号選択回路の出力と前記周
波数平均回路で算出した平均周波数とに基づいて
前記移動体の速度分布を表示するように構成され
ている。 (viii) 請求項８の発明 請求項８の超音波診断装置において、カラード
プラ解析器は、移動体から反射した超音波受信信
号に基づいて前記移動体の速度に対応した周波数
ｆ、該周波数解析時刻における前記移動体の移動
に伴う超音波受信信号の周波数偏移の正負を表す
符号ｓ及び前記超音波受信信号のパワーｐを解析
して出力する。 平均符号回路は、前記カラードプラ解析器によ
つて解析された所定期間内の各周波数解析時刻に
出力される各符号ｓの平均を算出する。 周波数平均回路は、前記カラードプラ解析器に
よつて解析された所定期間における周波数ｆの平
均を算出する。 周波数比較器は、前記周波数平均回路で算出し
た平均周波数と第１の閾値との大小を比較する。 パワー比較器は、前記カラードプラ解析器によ
つて解析されたパワーｐと第２の閾値との大小を
比較する。 符号選択回路は、前記周波数比較器の比較結果
と前記パワー比較器の比較結果とが入力され、前
記周波数平均回路で算出された平均周波数が前記
第１の閾値より小さく、かつ前記カラードプラ解
析器で解析されたパワーｐが前記第２の閾値より
小さい場合に前記平均符号回路で算出された平均
符号を出力し、それ以外の場合に前記カラードプ
ラ解析器で解析された符号ｓを出力する。 全体として、前記符号選択回路の出力と前記カ
ラードプラ解析器によつて解析された周波数ｆと
に基づいて前記移動体の速度分布を表示するよう
に構成されている。 〔作 用〕 (i) 請求項１の発明 周波数平均回路によつて算出された平均周波数
が所定の閾値より大きい場合には、カラードプラ
解析器によつて解析された符号ｓを選択し、この
選択した符号ｓと周波数平均回路で算出した平均
周波数とに基づいて移動体の速度分布の表示を行
う。 一方、平均周波数が所定の閾値より小さい場合
には、符号平均回路によつて算出した平均符号を
選択し、この選択した平均符号と周波数平均回路
で算出した平均周波数とに基づいて移動体の速度
分布の表示を行う。 請求項１の発明にあつては、移動体の速度が遅
く符号ｓに含まれるノイズの割合が多い場合に、
符号ｓに代えて平均符号が選択される。 (ii) 請求項２の発明 周波数平均回路によつて算出された平均周波数
が所定の閾値より大きい場合には、カラードプラ
解析器によつて解析された符号ｓを選択し、この
選択した符号ｓとカラードプラ解析器によつて解
析された周波数ｆとに基づいて移動体の速度分布
の表示を行う。 一方、平均周波数が所定の閾値より小さい場合
には、符号平均回路によつて算出した平均符号を
選択し、この選択した平均符号とカラードプラ解
析器によつて解析された周波数ｆとに基づいて移
動体の速度分布の表示を行う。 請求項２の発明にあつては、移動体の速度が遅
く符号ｓに含まれるノイズの割合が多い場合に、
符号ｓに代えて平均符号が選択される。 (iii) 請求項３の発明 カラードプラ解析器によつて解析されたパワー
ｐが所定の閾値より大きい場合には、カラードプ
ラ解析器によつて解析された符号ｓを選択し、こ
の選択した符号ｓと周波数平均回路で算出した平
均周波数とに基づいて移動体の速度分布の表示を
行う。 一方、パワーｐが所定の閾値より小さい場合に
は、符号平均回路によつて算出した平均符号を選
択し、この選択した平均符号と周波数平均回路で
算出した平均周波数とに基づいて移動体の速度分
布の表示を行う。 請求項３の発明にあつては、反射してくる超音
波受信信号の強度が弱く符号ｓに含まれるノイズ
の割合が多い場合に、符号ｓに代えて平均符号が
選択される。 (iv) 請求項４の発明 カラードプラ解析器によつて解析されたパワー
ｐが所定の閾値より大きい場合には、カラードプ
ラ解析器によつて解析された符号ｓを選択し、こ
の選択した符号ｓとカラードプラ解析器によつて
解析された周波数ｆとに基づいて移動体の速度分
布の表示を行う。 一方、パワーｐが所定の閾値より小さい場合に
は、符号平均回路によつて算出した平均符号を選
択し、この選択した平均符号とカラードプラ解析
器によつて解析された周波数ｆとに基づいて移動
体の速度分布の表示を行う。 請求項４の発明にあつては、反射してくる超音
波受信信号の強度が弱く符号ｓに含まれるノイズ
の割合が多い場合に、符号ｓに代えて平均符号が
選択される。 (v) 請求項５の発明 周波数平均回路によつて算出された平均周波数
が第１の閾値より大きく、かつカラードプラ解析
器によつて解析されたパワーｐが第２の閾値より
大きい場合には、カラードプラ解析器によつて解
析された符号ｓを選択し、この選択した符号ｓと
周波数平均回路で算出した平均周波数とに基づい
て移動体の速度分布の表示を行う。 一方、それ以外の場合には、符号平均回路によ
つて算出した平均符号を選択し、この選択した平
均符号と周波数平均回路で算出した平均周波数と
に基づいて移動体の速度分布の表示を行う。 請求項５の発明にあつては、移動体の速度が遅
く、かつ反射してくる超音波受信信号の強度が弱
く符号ｓに含まれるノイズの割合が多い場合に、
符号ｓに代えて平均符号が選択される。 (vi) 請求項６の発明 周波数平均回路によつて算出された平均周波数
が第１の閾値より大きく、かつカラードプラ解析
器によつて解析されたパワーｐが第２の閾値より
大きい場合には、カラードプラ解析器によつて解
析された符号ｓを選択し、この選択した符号ｓと
カラードプラ解析器によつて解析された周波数ｆ
とに基づいて移動体の速度分布の表示を行なう。 一方、それ以外の場合には、符号平均回路によ
つて算出した平均符号を選択し、この選択した平
均符号とカラードプラ解析器によつて解析された
周波数ｆとに基づいて移動体の速度分布の表示を
行う。 請求項６の発明にあつては、移動体の速度が遅
く、かつ反射してくる超音波受信信号の強度が弱
く符号ｓに含まれるノイズの割合が多い場合に、
符号ｓに代えて平均符号が選択される。 (vii) 請求項７の発明 周波数平均回路によつて算出された平均周波数
が第１の閾値より小さく、かつカラードプラ解析
器によつて解析されたパワーｐが第２の閾値より
小さい場合には、符号平均回路によつて算出した
平均符号を選択し、この選択した平均符号と周波
数平均回路で算出した平均周波数とに基づいて移
動体の速度分布の表示を行う。 一方、それ以外の場合には、カラードプラ解析
器によつて解析された符号ｓを選択し、この選択
した符号ｓと周波数平均回路で算出した平均周波
数とに基づいて移動体の速度分布の表示を行う。 請求項７の発明にあつては、移動体の速度が遅
いか、あるいは反射してくる超音波受信信号の強
度が弱く符号ｓに含まれるノイズの割合が多い場
合に、符号ｓに代えて平均符号が選択される。 (viii) 請求項８の発明 周波数平均回路によつて算出された平均周波数
が第１の閾値より小さく、かつカラードプラ解析
器によつて解析されたパワーｐが第２の閾値より
小さい場合には、符号平均回路によつて算出した
平均符号を選択し、この選択した平均符号とカラ
ードプラ解析器によつて解析された周波数ｆとに
基づいて移動体の速度分布の表示を行う。 一方、それ以外の場合には、カラードプラ解析
器によつて解析された符号ｓを選択し、この選択
した符号ｓとカラードプラ解析器によつて解析さ
れた周波数ｆとに基づいて移動体の速度分布の表
示を行う。 請求項８の発明にあつては、移動体の速度が遅
いか、あるいは反射してくる超音波受信信号の強
度が弱く符号ｓに含まれるノイズの割合が多い場
合に、符号ｓに代えて平均符号が選択される。 〔実施例〕 以下、図面に基づいて本発明の実施例について
詳細に説明する。 第１図に、本発明の超音波診断装置の一実施例
の構成を示す。 第１図において、１１１は超音波プローブ、１
３は送信回路、１１５は受信回路、１２１はカラ
ードプラ解析器、１２３は雑音除去回路、１２５
及び１４３はデイジタル・スキヤン・コンバー
タ、１２７及び１４５はメモリ、１２９及び１４
７はルツク・アツプ・テーブル、１３１は表示装
置、１４１は振幅検出回路、１４９は合成回路で
ある。 超音波プローブ１１１は、例えばPhased
Array型と称される構成を有しており、数十個程
度の超音波振動子を備えている。この超音波プロ
ーブ１１１は、送信回路１１３及び受信回路１１
５に接続されている。 送信回路１１３は、超音波プローブ１１１をド
ライブするためのものであり、超音波プローブ１
１１に備わつた超音波振動子を励振することによ
り、所定の照射角をなす超音波を出力させる。ま
た、受信回路１１５は、被測定物（移動体）によ
つて反射された超音波が超音波プローブ１１１に
入射されたときに、この反射波を検出するための
ものであり、上述した超音波振動子を順次走査す
ることにより反射波の検出を行つている。 カラードプラ解析器１２１は、受信回路１１５
の受信出力に基づいて被測定物体の速度に関する
解析を行うものである。解析の結果、カラードプ
ラ解析器１２１は、被測定物体の速度に対応した
周波数出力ｆと、この速度の方向に対応した符号
ｓと、反射波の強度に対応したパワーｐとを出力
する。カラードプラ解析器の出力は雑音除去回路
１２３に入力される。 雑音除去回路１２３は、予め設定された所定の
モードに従つて、カラードプラ解析器１２１から
出力された符号ｓに補正を加えるものである。こ
の補正は、ノイズによつて符号が反転する現象を
除去するための処理であり、ノイズによつて符号
が反転する可能性がある状態を検出したときに、
符号ｓに代えて所定時間の平均値を出力する。
この符号平均については、表示する画素の前後例
えば連続する15画素の符号の平均をとり、その15
画素に含まれる１つの画素の符号を決定する。雑
音除去回路１２３から出力される符号ｓあるいは
その平均値と周波数ｆはデイジタル・スキヤ
ン・コンバータ１２５に入力される。 デイジタル・スキヤン・コンバータ１２５は、
雑音除去回路１２３の出力をスキヤンして順次取
り込むためのものであり、この取り込んだデータ
はメモリ１２７に入力される。このデイジタル・
スキヤン・コンバータ１２５によるスキヤン動作
は、超音波プローブ１１１から反射波を取り込む
スキヤン動作に同期している。メモリ１２７で
は、デイジタル・スキヤン・コンバータ１２５か
ら出力されるデータを、被測定物体の位置に応じ
たアドレスに順次格納する。このメモリ１２７の
格納データは順次読み出されてルツク・アツプ・
テーブル１２９に入力される。 ルツク・アツプ・テーブル１２９は、メモリ１
２７から読み出されたデータに基づいて、カラー
データ（例えばRGBデータ）を作成するための
ものである。メモリ１２７には、雑音除去回路１
２３から出力された周波数ｆ及び符号ｓあるいは
その平均値が格納されており、これらの格納デ
ータに基づいて被測定物体の速度に関するカラー
データに変換する。例えば、符号ｓあるいはその
平均値が正である場合が赤色に対応しており、
この赤色の濃淡の度合いを周波数ｆに対応させた
カラーデータへの変換を行う。反対に、符号ｓあ
るいはその平均値が負である場合が青色に対応
しており、この青色の濃淡の度合いを周波数ｆに
対応させたカラーデータへの変換が行われる。こ
の変換後のカラーデータは表示装置１３１に入力
される。 表示装置１３１は、ルツク・アツプ・テーブル
１２９から出力されるカラーデータに基づいて、
被測定物体の速度に関するカラー表示を行う。例
えば、被測定物体として血液を考えた場合、流れ
が一定方向の血流は赤色あるいは青色何れか一方
に対応して表示され、乱流と称される高速流の場
合は所々流れの方向が反転するので、青色の所々
に赤色が混じて表示されたり、赤色の所々に青色
が混じつて表示されたりする。 上述した表示装置１３１による表示は血液の流
れを表示するような場合に用いられるが、超音波
の反射波に基づいて被測定物体の形状を測定する
こともできる。以下、形状を測定する動作につい
て説明する。 振幅検出回路１４１は、受信回路１１５の受信
出力の振幅検波を行うものである。振幅検波回路
１４１は、受信信号の振幅に応じた出力、すなわ
ち超音波プローブ１１１に到達する反射波の強度
に応じた検波出力が得られる。この検波出力は、
デイジタル・スキヤン・コンバータ１４３でスキ
ヤンされてメモリ１４５に格納され、その後ルツ
ク・アツプ・テーブル１４７に入力される。ルツ
ク・アツプ・テーブル１４７は、メモリ１４５か
ら読み出されるデータに基づいて白黒データ
（RGBが同一値）を作成するためのものである。
例えば、心臓の断面に超音波を照射した場合に
は、心臓の壁面部分と内部の血液が流れている部
分とでは反射波の強度が異なるため振幅検波回路
１４１の検波出力にも差異が生じ、この差異が被
測定物体の形状として表示装置１３１に表示され
る。合成回路１４９はルツク・アツプ・テーブル
１２９，１４７の出力を合成して表示装置１３１
に供給している。 次に、上述した本実施例の詳細構成及び動作を
第２図乃至第１０図を参照して説明する。 第２図に、カラードプラ解析器１２１の詳細構
成を示す。２１１は直交検波回路、２１３はアナ
ログ−デイジタル（Ａ／Ｄ）変換器、２１５はモ
ーシヨン・ターゲツト・インジケータ（MTI）
フイルタ、２１７は自己相関器である。 受信回路１１５の出力が直交検波回路２１１に
入力されると、直交検波回路２１１は、照射波の
出力周波数fcと同一周波数の信号と、出力周波数
fcと90度位相がずれた周波数の信号とによつて入
力に対する直交検波を行う。この直交検波の結
果、位相が90度ずれた２つの出力が得られ、Ａ／
Ｄ変換器２１３に入力される。 Ａ／Ｄ変換器２１３は、入力された２つの信号
のそれぞれをデイジタルデータに変換するもので
あり、変換されたデイジタルデータはMTIフイ
ルタ２１５に入力される。 MTIフイルタ２１５は、Ａ／Ｄ変換器２１３
から入力される２つのデイジタルデータのそれぞ
れに対して、ハイパスフイルタとして動作する。
例えば、超音波プローブ１１１によつて被測定物
体の同一位置に対する反射波を繰り返し受信した
場合に、この同一位置の反射波に対応するＡ／Ｄ
変換器２１３の出力データに対して、１次のハイ
パスフイルタとして動作する。従つて、被測定物
体が移動した場合には、この移動速度に応じて
MTIフイルタ２１５の出力値が大きくなる。
MTIフイルタ２１５の２つの出力は、自己相関
器２１７に入力される。 自己相関器２１７は、MTIフイルタ２１５の
２つの出力に基づいて自己相関から周波数を求め
る。自己相関から得られる周波数をＦ、MTIフ
イルタ２１５から出力される２つの出力をＸ(i)、
Ｙ(i)（ｉは被測定物体の同一位置に対して繰り返
される超音波の送受信回数）とすると、自己相関
から得られる周波数Ｆは、 Ｆ＝１／Ｔtan-1｛Σ（Ｘ(i)Ｙ（ｉ−１
）−Ｙ(i)Ｘ（ｉ−１））／Σ（Ｘ(i)Ｘ（ｉ−１）＋Ｙ
(i)Ｙ（ｉ−１））｝…(1) で示される。ここで、Ｔは被測定物体の同一位置
に対して繰り返される超音波の送信あるいは受信
の周期である。 自己相関器２１７は、(1)式に従つて求めた自己
相関による周波数Ｆの絶対値｜Ｆ｜を周波数ｆと
して、自己相関による周波数Ｆの正負を符号ｓと
して出力する。また、自己相関器２１７は、被測
定物体の同一位置に対応する２入力Ｘ(i)、Ｙ(i)の
２乗和を累積した値Σ（Ｘ(i)2＋Ｙ(i)2）をパワー
ｐとして出力する。このパワーｐは、反射波の強
度に対応している。自己相関器２１７のこれらの
出力は、上述した雑音除去回路１２３に入力され
る。 第３図に、雑音除去回路１２３の詳細構成を示
す。第３図において、３１１はパワー比較器、３
１３は周波数選択回路、３１５は周波数平均回
路、３１７は周波数比較、３２１は平均符号回
路、３２３は等価遅延回路、３２５は符号選択回
路、３３１はモード制御回路である。 パワー比較器３１１は、カラードプラ解析器１
２１内の自己相関器２１７から出力されたパワー
ｐが入力され、このパワーｐと所定の閾値p0との
大小を比較するものである。比較結果は、周波数
選択回路３１３及びモード制御回路３３１に入力
される。 周波数選択回路３１３は、カラードプラ解析器
１２１内の自己相関器２１７から出力された周波
数ｆが入力される。周波数選択回路３１３は、パ
ワー比較器３１１の比較結果に応じて、この周波
数ｆをそのまま出力するか、あるいは零に置き換
えて出力するかの択一的動作を行つている。具体
的には、パワー比較器３１１に入力されたパワー
ｐが所定の閾値p0以下のときは、周波数ｆはノイ
ズの影響を大きく受けているものであるため、こ
の周波数ｆを零に置き換えて出力する。それ以外
のときは周波数ｆをそのまま出力する。周波数選
択回路３１３の出力は、周波数平均回路３１５に
入力される。 周波数平均回路３１５は、周波数選択回路３１
３の出力に対して、所定の時間内の移動平均を求
めて平均周波数を出力するためのものである。
周波数ｆから平均周波数を求めることにより、
画像の相関を利用した補正を行つて、後段の表示
装置１３１による表示を見易いものにしている。
周波数平均回路３１５の出力は、周波数比較器３
１７に入力されると共に、雑音除去回路１２３の
外部に出力される。 周波数比較器３１７は、入力された平均周波数
ｆと所定の閾値f0との大小を比較するものであ
る。比較結果は、モード制御回路３３１に入力さ
れる。 モード制御回路３３１は、パワー比較器３１１
及び周波数比較器３１７の各出力の少なくとも一
方に基づいて、符号選択回路３２５による選択動
作を制御するものである。モード制御回路３３１
に入力されるモード選択信号に応じて制御モード
の設定が行われ、設定された制御モードに従つた
制御信号が符号選択回路３２５に入力される。 平均符号回路３２１及び等価遅延回路３２３に
は、カラードプラ解析器１２１内の自己相関器２
１７から出力された符号ｓが入力される。平均符
号回路３２１は、この符号ｓを所定の時間内で平
均した平均符号を算出して出力する。この平均
符号は、表示する画素に対応する符号につき、例
えば連続する７画素の平均をとり、平均符号と
して出力するものである。平均符号回路３２１か
ら出力された平均符号は符号選択回路３２５に
入力される。また、等価遅延回路３２３は、平均
符号回路３２１による平均符号算出に要する時間
と同期をとるためのものであり、入力された符号
ｓを所定時間遅延した後に出力する。等価遅延回
路３２３で遅延された符号ｓは符号選択回路３２
５に入力される。 符号選択回路３２５は、平均符号回路３２１か
ら入力される平均符号または等価遅延回路３２
３から出力される符号ｓの何れか一方を選択して
出力する。この選択動作は、モード制御回路３３
１から入力される制御信号によつて制御される。
符号選択回路３２５から出力される符号ｓあるい
は平均符号は、雑音除去回路１２３の外部に出
力される。 第１表に、モード制御回路３３１において制御
される制御モードの詳細内容を示す。以下、各制
御モードについて説明する。
【表】
請求項１及び請求項２の発明によれば、移動体
の速度が遅く符号ｓに含まれるノイズの割合が多
い場合に、符号ｓに代えて平均符号を用いて移動
体の速度分布の表示を行うことにより、ノイズに
よる画像の乱れを防止することができる。 請求項３及び請求項４の発明によれば、反射し
てくる超音波受信信号の強度が弱く符号ｓに含ま
れるノイズの割合が多い場合に、符号ｓに代えて
平均符号を用いて移動体の速度分布の表示を行う
ことにより、ノイズによる画像の乱れを防止する
ことができる。 請求項５及び請求項６の発明にあつては、移動
体の速度が遅く、かつ反射してくる超音波受信信
号の強度が弱く符号ｓに含まれるノイズの割合が
多い場合に、符号ｓに代えて平均符号を用いて移
動体の速度分布の表示を行うことにより、ノイズ
による画像の乱れを防止することができる。 請求項７及び請求項８の発明によれば、移動体
の速度が遅いか、あるいは反射してくる超音波受
信信号の強度が弱く符号ｓに含まれるノイズの割
合が多い場合に、符号ｓに代えて平均符号を用い
て移動体の速度分布の表示を行うことにより、ノ
イズによる画像の乱れを防止することができる。
の速度が遅く符号ｓに含まれるノイズの割合が多
い場合に、符号ｓに代えて平均符号を用いて移動
体の速度分布の表示を行うことにより、ノイズに
よる画像の乱れを防止することができる。 請求項３及び請求項４の発明によれば、反射し
てくる超音波受信信号の強度が弱く符号ｓに含ま
れるノイズの割合が多い場合に、符号ｓに代えて
平均符号を用いて移動体の速度分布の表示を行う
ことにより、ノイズによる画像の乱れを防止する
ことができる。 請求項５及び請求項６の発明にあつては、移動
体の速度が遅く、かつ反射してくる超音波受信信
号の強度が弱く符号ｓに含まれるノイズの割合が
多い場合に、符号ｓに代えて平均符号を用いて移
動体の速度分布の表示を行うことにより、ノイズ
による画像の乱れを防止することができる。 請求項７及び請求項８の発明によれば、移動体
の速度が遅いか、あるいは反射してくる超音波受
信信号の強度が弱く符号ｓに含まれるノイズの割
合が多い場合に、符号ｓに代えて平均符号を用い
て移動体の速度分布の表示を行うことにより、ノ
イズによる画像の乱れを防止することができる。
第１図は本発明の超音波診断装置の一実施例の
構成を示す図、第２図は本実施例におけるカラー
ドプラ解析器の詳細構成を示す図、第３図は本実
施例における雑音除去回路の詳細構成を示す図、
第４図は本実施例における平均符号回路の詳細構
成を示す図、第５図は本実施例における等価遅延
回路の詳細構成を示す図、第６図は本実施例にお
ける周波数平均回路の詳細構成を示す図、第７図
は本実施例におけるモード制御回路の詳細構成を
示す図、第８図は実施例における符号選択回路及
びモード制御回路内の選択回路の具体例の構成を
示す図、第９図は本発明の超音波診断装置の別の
実施例の構成を示す図、第１０図は本発明の超音
波診断装置の更に別の実施例の構成を示す図であ
る。 図において、１２１はカラードプラ解析器、３
１１はパワー比較器、３１３は周波数選択回路、
３１５は周波数平均回路、３１７は周波数比較
器、３２１は平均符号回路、３２３は等価遅延回
路、３２５は符号選択回路、３３１はモード制御
回路である。
構成を示す図、第２図は本実施例におけるカラー
ドプラ解析器の詳細構成を示す図、第３図は本実
施例における雑音除去回路の詳細構成を示す図、
第４図は本実施例における平均符号回路の詳細構
成を示す図、第５図は本実施例における等価遅延
回路の詳細構成を示す図、第６図は本実施例にお
ける周波数平均回路の詳細構成を示す図、第７図
は本実施例におけるモード制御回路の詳細構成を
示す図、第８図は実施例における符号選択回路及
びモード制御回路内の選択回路の具体例の構成を
示す図、第９図は本発明の超音波診断装置の別の
実施例の構成を示す図、第１０図は本発明の超音
波診断装置の更に別の実施例の構成を示す図であ
る。 図において、１２１はカラードプラ解析器、３
１１はパワー比較器、３１３は周波数選択回路、
３１５は周波数平均回路、３１７は周波数比較
器、３２１は平均符号回路、３２３は等価遅延回
路、３２５は符号選択回路、３３１はモード制御
回路である。
Claims (1)
- 【特許請求の範囲】 １ 移動体から反射した超音波受信信号に基づい
て前記移動体の速度に対応した周波数ｆ、及び該
周波数解析時刻における前記移動体の移動に伴う
超音波受信信号の周波数偏移の正負を表す符号ｓ
を解析して出力するカラードプラ解析器と、 前記カラードプラ解析器によつて解析された所
定期間内の各周波数解析時刻に出力される各符号
ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によつて解析された所
定期間における周波数ｆの平均を算出する周波数
平均回路と、 前記周波数平均回路で算出した平均周波数と所
定の閾値との大小を比較する周波数比較器と、 前記周波数比較器の比較結果が入力され、前記
周波数平均回路で算出した平均周波数が前記所定
の閾値より大きい場合に前記カラードプラ解析器
によつて解析された符号ｓを出力し、前記周波数
平均回路で算出した平均周波数が前記所定の閾値
より小さい場合に前記平均符号回路で算出した平
均符号を出力する符号選択回路と、 を備え、この符号選択回路の出力と前記周波数平
均回路で算出した平均周波数とに基づいて前記移
動体の速度分布を表示することを特徴とする超音
波診断装置。 ２ 移動体から反射した超音波受信信号に基づい
て前記移動体の速度に対応した周波数ｆ、及び該
周波数解析時刻における前記移動体の移動に伴う
超音波受信信号の周波数偏移の正負を表す符号ｓ
を解析して出力するカラードプラ解析器と、 前記カラードプラ解析器によつて解析された所
定期間内の各周波数解析時刻に出力される各符号
ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によつて解析された所
定期間における周波数ｆの平均を算出する周波数
平均回路と、 前記周波数平均回路で算出した平均周波数と所
定の閾値との大小を比較する周波数比較器と、 前記周波数比較器の比較結果が入力され、前記
周波数平均回路で算出した平均周波数が前記所定
の閾値より大きい場合に前記カラードプラ解析器
によつて解析された符号ｓを出力し、前記周波数
平均回路で算出した平均周波数が前記所定の閾値
より小さい場合に前記平均符号回路で算出した平
均符号を出力する符号選択回路と、 を備え、この符号選択回路の出力と前記カラード
プラ解析器によつて解析された周波数ｆとに基づ
いて前記移動体の速度分布を表示することを特徴
とする超音波診断装置。 ３ 移動体から反射した超音波受信信号に基づい
て前記移動体の速度に対応した周波数ｆ、該周波
数解析時刻における前記移動体の移動に伴う超音
波受信信号の周波数偏移の正負を表す符号ｓ及び
前記超音波受信信号のパワーｐを解析して出力す
るカラードプラ解析器と、 前記カラードプラ解析器によつて解析された所
定期間内の各周波数解析時刻に出力される各符号
ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によつて解析された所
定期間における周波数ｆの平均を算出する周波数
平均回路と、 前記カラードプラ解析器によつて解析されたパ
ワーｐと所定の閾値との大小を比較するパワー比
較器と、 前記パワー比較器の比較結果が入力され、前記
カラードプラ解析器で解析されたパワーｐが前記
所定の閾値より大きい場合に前記カラードプラ解
析器によつて解析された符号ｓを出力し、前記カ
ラードプラ解析器で解析されたパワーｐが前記所
定の閾値より小さい場合に前記平均符号回路で算
出した平均符号を出力する符号選択回路と、 を備え、この符号選択回路の出力と前記周波数平
均回路で算出した平均周波数とに基いて前記移動
体の速度分布を表示することを特徴とする超音波
診断装置。 ４ 移動体から反射した超音波受信信号に基づい
て前記移動体の速度に対応した周波数ｆ、該周波
数解析時刻における前記移動体の移動に伴う超音
波受信信号の周波数偏移の正負を表す符号ｓ及び
前記超音波受信信号のパワーｐを解析して出力す
るカラードプラ解析器と、 前記カラードプラ解析器によつて解析された所
定期間内の各周波数解析時刻に出力される各符号
ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によつて解析されたパ
ワーｐと所定の閾値との大小を比較するパワー比
較器と、 前記パワー比較器の比較結果が入力され、前記
カラードプラ解析器で解析されたパワーｐが前記
所定の閾値より大きい場合に前記カラードプラ解
析器によつて解析された符号ｓを出力し、前記カ
ラードプラ解析器で解析されたパワーｐが前記所
定の閾値より小さい場合に前記平均符号回路で算
出した平均符号を出力する符号選択回路と、 を備え、この符号選択回路の出力と前記カラード
プラ解析器によつて解析された周波数ｆとに基づ
いて前記移動体の速度分布を表示することを特徴
とする超音波診断装置。 ５ 移動体から反射した超音波受信信号に基づい
て前記移動体の速度に対応した周波数ｆ、該周波
数解析時刻における前記移動体の移動に伴う超音
波受信信号の周波数偏移の正負を表す符号ｓ及び
前記超音波受信信号のパワーｐを解析して出力す
るカラードプラ解析器と、 前記カラードプラ解析器によつて解析された所
定期間内の各周波数解析時刻に出力される各符号
ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によつて解析された所
定移動体における周波数ｆの平均を算出する周波
数平均回路と、 前記周波数平均回路で算出した平均周波数と第
１の閾値との大小を比較する周波数比較器と、 前記カラードプラ解析器によつて解析されたパ
ワーｐと第２の閾値との大小を比較するパワー比
較器と、 前記周波数比較器の比較結果と前記パワー比較
器の比較結果とが入力道され、前記周波数平均回
路で算出された平均周波数が前記第１の閾値より
大きく、かつ前記カラードプラ解析器で解析され
たパワーｐが前記第２の閾値より大きい場合に前
記カラードプラ解析器によつて解析された符号ｓ
を出力し、それ以外の場合に前記平均符号回路で
算出した平均符号を出力する符号選択回路と、 を備え、この符号選択回路の出力と前記周波数平
均回路で算出した平均周波数とに基づいて前記移
動体の速度分布を表示することを特徴とする超音
波診断装置。 ６ 移動体から反射した超音波受信信号に基づい
て前記移動体の速度に対応した周波数ｆ、該周波
数解析時刻における前記移動体の移動に伴う超音
波受信信号の周波数偏移の正負を表す符号ｓ及び
前記超音波受信信号のパワーｐを解析して出力す
るカラードプラ解析器と、 前記カラードプラ解析器によつて解析された所
定期間内の各周波数解析時刻に出力される各符号
ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によつて解析された所
定期間における周波数ｆの平均を算出する周波数
平均回路と、 前記周波数平均回路で算出した平均周波数と第
１の閾値との大小を比較する周波数比較器と、 前記カラードプラ解析器によつて解析されたパ
ワーｐと第２の閾値との大小を比較するパワー比
較器と、 前記周波数比較器の比較結果と前記パワー比較
器の比較結果とが入力され、前記周波数平均回路
で算出された平均周波数が前記第１の閾値より大
きく、かつ前記カラードプラ解析器で解析された
パワーｐが前記第２の閾値より大きい場合に前記
カラードプラ解析器によつて解析された符号ｓを
出力し、それ以外の場合に前記平均符号回路で算
出した平均符号を出力する符号選択回路と、 を備え、この符号選択回路の出力と前記カラード
プラ解析器によつて解析された周波数ｆとに基づ
いて前記移動体の速度分布を表示することを特徴
とする超音波診断装置。 ７ 移動体から反射した超音波受信信号に基づい
て前記移動体の速度に対応した周波数ｆ、該周波
数解析時刻における前記移動体の移動に伴う超音
波受信信号の周波数偏移の正負を表す符号ｓ及び
前記超音波受信信号のパワーｐを解析して出力す
るカラードプラ解析器と、 前記カラードプラ解析器によつて解析された所
定期間内の各周波数解析時刻に出力される各符号
ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によつて解析された所
定期間における周波数ｆの平均を算出する周波数
平均回路と、 前記周波数平均回路で算出した平均周波数と第
１の閾値との大小を比較する周波数比較器と、 前記カラードプラ解析器によつて解析されたパ
ワーｐと第２の閾値との大小を比較するパワー比
較器と、 前記周波数比較器の比較結果と前記パワー比較
器の比較結果とが入力され、前記周波数平均回路
で算出された平均周波数が前記第１の閾値より小
さく、かつ前記カラードプラ解析器で解析された
パワーｐが前記第２の閾値より小さい場合に前記
平均符号回路で算出された平均符号を出力し、そ
れ以外の場合に前記カラードプラ解析器で解析さ
れた符号ｓを出力する符号選択回路と、 を備え、この符号選択回路の出力と前記周波数平
均回路で算出した平均周波数とに基づいて前記移
動体の速度分布を表示することを特徴とする超音
波診断装置。 ８ 移動体から反射した超音波受信信号に基づい
て前記移動体の速度に対応した周波数ｆ、該周波
数解析時刻における前記移動体の移動に伴う超音
波受信信号の周波数偏移の正負を表す符号ｓ及び
前記超音波受信信号のパワーｐを解析して出力す
るカラードプラ解析器と、 前記カラードプラ解析器によつて解析された所
定期間内の各周波数解析時刻に出力される各符号
ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によつて解析された所
定期間における周波数ｆの平均を算出する周波数
平均回路と、 前記周波数平均回路で算出した平均周波数と第
１の閾値との大小を比較する周波数比較器と、 前記カラードプラ解析器によつて解析されたパ
ワーｐと第２の閾値との大小を比較するパワー比
較器と、 前記周波数比較器の比較結果と前記パワー比較
器の比較結果とが入力され、前記周波数平均回路
で算出された平均周波数が前記第１の閾値より小
さく、かつ前記カラードプラ解析器で解析された
パワーｐが前記第２の閾値より小さい場合に前記
平均符号回路で算出された平均符号を出力して、
それ以外の場合に前記カラードプラ解析器で解析
された符号ｓを出力する符号選択回路と、 を備え、この符号選択回路の出力と前記カラード
プラ解析器によつて解析された周波数ｆとに基づ
いて前記移動体の速度分布を表示することを特徴
とする超音波診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007986A JPH02289237A (ja) | 1989-01-17 | 1990-01-16 | 超音波診断装置 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP806589 | 1989-01-17 | ||
| JP1-8065 | 1989-01-17 | ||
| JP2007986A JPH02289237A (ja) | 1989-01-17 | 1990-01-16 | 超音波診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH02289237A JPH02289237A (ja) | 1990-11-29 |
| JPH0572817B2 true JPH0572817B2 (ja) | 1993-10-13 |
Family
ID=11682945
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007986A Granted JPH02289237A (ja) | 1989-01-17 | 1990-01-16 | 超音波診断装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5042491A (ja) |
| EP (1) | EP0379139B1 (ja) |
| JP (1) | JPH02289237A (ja) |
| DE (1) | DE69010923T2 (ja) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5216602A (en) * | 1989-11-06 | 1993-06-01 | The Board Of Trustees Of The University Of Illinois | Color imaging system |
| DE4134724C2 (de) * | 1990-10-24 | 1995-11-16 | Hitachi Medical Corp | Einrichtung zur farbigen Strömungsaufzeichnung mit Ultraschall |
| JPH0793927B2 (ja) * | 1990-11-02 | 1995-10-11 | 富士通株式会社 | 超音波カラードプラ診断装置 |
| JP2931707B2 (ja) * | 1991-10-04 | 1999-08-09 | 株式会社東芝 | 超音波診断装置 |
| JP3088174B2 (ja) * | 1992-02-06 | 2000-09-18 | 古野電気株式会社 | 水中探知装置 |
| US5373847A (en) * | 1992-11-12 | 1994-12-20 | University Of Pittsburgh Of The Commonwealth System Of Higher Education | Method of color doppler mapping of blood flow in a patient |
| US5857973A (en) * | 1997-09-30 | 1999-01-12 | Siemens Medical Systems, Inc. | Fuzzy logic tissue flow determination system |
| US6066095A (en) * | 1998-05-13 | 2000-05-23 | Duke University | Ultrasound methods, systems, and computer program products for determining movement of biological tissues |
| CN100469320C (zh) * | 2005-09-29 | 2009-03-18 | 深圳迈瑞生物医疗电子股份有限公司 | 消除彩色血流图像中速度异常点的方法 |
| JP2009261436A (ja) * | 2008-04-22 | 2009-11-12 | Aloka Co Ltd | 超音波診断装置 |
| CN102429684B (zh) | 2010-09-28 | 2013-10-09 | 深圳迈瑞生物医疗电子股份有限公司 | 一种多普勒彩色血流成像方法和装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63315040A (ja) * | 1987-06-19 | 1988-12-22 | Toshiba Corp | 超音波診断装置 |
| JPS63317142A (ja) * | 1987-06-22 | 1988-12-26 | Toshiba Corp | 超音波診断装置 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA1135827A (en) * | 1978-12-04 | 1982-11-16 | Rainer Fehr | Determination of flow velocities by measuring phase difference between the doppler signals |
| US4324258A (en) * | 1980-06-24 | 1982-04-13 | Werner Huebscher | Ultrasonic doppler flowmeters |
| JPS5920820A (ja) * | 1982-07-28 | 1984-02-02 | Aloka Co Ltd | 超音波血流画像形成装置 |
| WO1985002105A1 (en) * | 1983-11-10 | 1985-05-23 | Acoustec Partners | Ultrasound diagnostic apparatus |
| JPS61141347A (ja) * | 1984-12-12 | 1986-06-28 | 株式会社 日立メデイコ | 超音波診断装置 |
| JP2639391B2 (ja) * | 1986-03-06 | 1997-08-13 | ヒューレット・パッカード・カンパニー | 流体マップの画像形成装置 |
| JPS63317137A (ja) * | 1987-06-20 | 1988-12-26 | Toshiba Corp | 超音波血流イメ−ジング装置 |
| US4850364A (en) * | 1987-11-12 | 1989-07-25 | Hewlett-Packard Company | Medical ultrasound imaging system with velocity-dependent rejection filtering |
| US4932415A (en) * | 1988-11-14 | 1990-06-12 | Vingmed Sound A/S | Method of color coding two dimensional ulltrasonic doppler velocity images of blood flow on a display |
-
1990
- 1990-01-16 JP JP2007986A patent/JPH02289237A/ja active Granted
- 1990-01-16 EP EP90100814A patent/EP0379139B1/en not_active Expired - Lifetime
- 1990-01-16 DE DE69010923T patent/DE69010923T2/de not_active Expired - Fee Related
- 1990-01-17 US US07/466,662 patent/US5042491A/en not_active Expired - Lifetime
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63315040A (ja) * | 1987-06-19 | 1988-12-22 | Toshiba Corp | 超音波診断装置 |
| JPS63317142A (ja) * | 1987-06-22 | 1988-12-26 | Toshiba Corp | 超音波診断装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH02289237A (ja) | 1990-11-29 |
| EP0379139A1 (en) | 1990-07-25 |
| DE69010923D1 (de) | 1994-09-01 |
| EP0379139B1 (en) | 1994-07-27 |
| US5042491A (en) | 1991-08-27 |
| DE69010923T2 (de) | 1994-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6277075B1 (en) | Method and apparatus for visualization of motion in ultrasound flow imaging using continuous data acquisition | |
| US7713204B2 (en) | Image data processing method and apparatus for ultrasonic diagnostic apparatus, and image processing apparatus | |
| US8036856B2 (en) | Method and apparatus for automatically adjusting spectral doppler gain | |
| US6478742B1 (en) | PRF adjustment method and apparatus, and ultrasonic wave imaging apparatus | |
| JPH0613031B2 (ja) | 超音波血流イメ−ジング装置 | |
| KR100742467B1 (ko) | 촬상 시스템, 촬상 방법 및 혈액 움직임 촬상 시스템 | |
| US5301670A (en) | Ultrasonic diagnosis apparatus | |
| US5188113A (en) | Ultrasonic diagnosis apparatus | |
| JPS6244494B2 (ja) | ||
| KR100352639B1 (ko) | 칼라 도플러 영상화 시스템을 위한 칼라 영상 표시방법 및 장치 | |
| US5148808A (en) | Ultrasonic diagnostic apparatus detecting Doppler shift | |
| JPH062134B2 (ja) | 超音波診断装置 | |
| JPH03188841A (ja) | 超音波診断装置 | |
| JPH0548692B2 (ja) | ||
| JPH0572817B2 (ja) | ||
| JP2000149015A (ja) | 画像のエッジ強調方法及びイメ―ジング装置 | |
| JP5455567B2 (ja) | 超音波診断装置 | |
| JP2823252B2 (ja) | 超音波診断装置 | |
| JP2987109B2 (ja) | 超音波ドプラ診断装置 | |
| JPH0549639A (ja) | 超音波カラードツプラー診断装置 | |
| JP2602898B2 (ja) | 超音波血流イメージング装置 | |
| JPH0556973A (ja) | 超音波血流イメージング装置 | |
| JPH0216138B2 (ja) | ||
| JPH04250148A (ja) | 超音波診断装置 | |
| JPH064073B2 (ja) | 超音波血流イメージング装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 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: 20091013 Year of fee payment: 16 |
|
| LAPS | Cancellation because of no payment of annual fees |