JPH02289237A - 超音波診断装置 - Google Patents

超音波診断装置

Info

Publication number
JPH02289237A
JPH02289237A JP2007986A JP798690A JPH02289237A JP H02289237 A JPH02289237 A JP H02289237A JP 2007986 A JP2007986 A JP 2007986A JP 798690 A JP798690 A JP 798690A JP H02289237 A JPH02289237 A JP H02289237A
Authority
JP
Japan
Prior art keywords
frequency
code
average
circuit
color doppler
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
JP2007986A
Other languages
English (en)
Other versions
JPH0572817B2 (ja
Inventor
Shinichi Amamiya
慎一 雨宮
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.)
Fujitsu Ltd
Original Assignee
Fujitsu 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 Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2007986A priority Critical patent/JPH02289237A/ja
Publication of JPH02289237A publication Critical patent/JPH02289237A/ja
Publication of JPH0572817B2 publication Critical patent/JPH0572817B2/ja
Granted legal-status Critical Current

Links

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)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • Chemical & Material Sciences (AREA)
  • Hematology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Biochemistry (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (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)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔目　次〕 概要 産業上の利用分野 従来の技術 発明が解決しようとする課題 課題を解決するための手段 作用 実施例 発明の効果 〔（既　　要〕 移動体のカラー表示を行う超音波診断装置に関し、特に
血流の流速及び方向に対応した周波数出力及びその符号
出力をカラードプラ解析器によって得るようにした超音
波診断装置に関し、ノイズによる画像の乱れの防止を目
的とし、移動体から反射した超音波受信信号に基づいて
移動体の動きの方向に対応した符号ｓ、移動体の速度に
対応した周波数ｆ及び超音波受信信号のパワーｐを解析
するカラードプラ解析器と、カラードプラ解析器によっ
て解析された所定期間における符号Ｓの平均を算出する
平均符号回路と、カラードプラ解析器によって解析され
た所定期間における周波数ｆの平均を算出する周波数平
均回路と、周波数平均回路で算出した平均周波数と第１
の闇値との大小を比較する周波数比較器と、カラードプ
ラ解析器によって解析されたパワーｐと第２の闇値との
大小を比較するパワー比較器と、周波数比較器の比較結
果とパワー比較器の比較結果の少なくとも一方が入力さ
れ、各比較結果に応じてカラードプラ解析器によって解
析された符号Ｓあるいは平均符号回路で算出した平均符
号を選択して出力する符号選択回路とを備えるように構
成する。
［産業上の利用分野］ 本発明は、移動体のカラー表示を行う超音波診断装置に
関し、特に血流の流速及び方向に対応した周波数出力及
びその符号出力をカラードプラ解析器によって得るよう
にした超音波診断装置に関するものである。
〔従来の技術〕
カラードプラ解析器を利用した超音波診断装置は、人体
の表面から超音波を照射して、この照射波と反射波の周
波数のずれを測定することにより血流の速度を測定する
ものであり、例えば心臓内の血流を容易に測定すること
が可能になる。
このような超音波診断装置で利用されるカラードプラ解
析器では、反射波から求めた受信信号を照射波に対応す
る送信信号及びこの送信信号と位相が９０度ずれた信号
とによって直交検波して、位相が９０度ずれた２つの出
力信号Ｘ、Ｙを得る。
そして、この２つの出力信号Ｘ、Ｙ中の心臓壁からの信
号のような不用信号を除去した後このＸ。
Ｙに基づいて自己相関から周波数Ｆを求め、この求めた
周波数Ｆの絶対値を周波数ｆ、正負を符号ｓ、２つの出
力信号Ｘ、　Ｙの強さをパワーＰとして出力する。この
ようにしてカラードプラ解析器から出力された周波数ｒ
は移動物体すなわち血流の速度に対応し、符号Ｓはその
速度方向に対応する。また、パワーＰは反射波の強度に
対応する。
カラードプラ解析器を利用した超音波診断装置では、カ
ラードプラ解析器から出力される符号Ｓが正である場合
には周波数ｆに応じた濃淡の例えば青色を表示し、符号
Ｓが負である場合には周波数ｆに応じた濃淡の例えば赤
色を表示して、血流の分布に対応した画像表示を行う。
また、カラードプラ解析器から出力されたパワーｐが所
定の閾値以下のときは、反射波のＳＮ比が悪化するため
周波数ｆの値を零に置き換えて、ノイズの影響を除去す
るようにしている。
超音波診断装置の従来技術としては、特開昭６２−２７
０１．３９号公報が知られている。
［発明が解決しようとする課題］ ところで、上述した従来方式にあっては、カラードプラ
解析器から出力されたパワーｐが所定の闇値以下の場合
には周波数ｆを零に置き換える処理を行っているが、そ
の闇値を余り大きく設定すると測定精度の低下につなが
るため、この測定精度が許容できる範囲で闇値を設定す
ることになる。
従って、カラードプラ解析器から闇値付近のパワーに対
応する周波数ｆが出力された場合にはわずかではあるが
ノイズが含まれることになり、その影響を完全に除去す
ることができない。
このようにしてノイズが含まれた場合には、周波数ｆの
値及び符号Ｓはランダムな値をとるため局部的に符号Ｓ
の値が反転することが生じる。そのため、血流が一定方
向に流れて例えば青色の表示がなされているときに、所
々赤色が混じるといった画像の乱れが発生するという問
題点があった。
このような画像の乱れが発生ずると、病院等において表
示画像を見て診断を行う場合に、誤診断につながること
になる。
また、ノイズよって符号Ｓが反転することによる画像の
乱れを除去する方法としては、所定領域内の符号Ｓの平
均を求めて局所的な符号Ｓの変動を均一化する補正を行
う方法が考えられる。しかし、一般には、低速の血流に
は逆流は存在しないが、乱流と称される高速の血流には
局部的に逆流が発生するという事実が確認されている。
このため、上述したように符号Ｓの平均を求めて画質を
補正する場合には、乱流が発生しているときの局部的な
逆流も均一化してしまうため、このような補正を適用す
ることはできない。
本発明は、このような点にかんがみて創作されたもので
あり、ノイズによる画像の乱れを防止することができる
超音波診断装置を提供することを目的としている。
〔課題を解決するための手段〕
ｉ　壬ン１の　日 請求項１の超音波診断装置において、カラードプラ解析
器は、移動体から反射した超音波受信信号に基づいて移
動体の動きの方向に対応した符号Ｓ及び移動体の速度に
対応した周波数ｒを解析する。
平均符号回路は、カラードプラ解析器によって解析され
た所定期間における符号Ｓの平均を算出する。
周波数平均回路は、カラードプラ解析器によって解析さ
れた所定期間における周波数ｆの平均を算出する。
周波数比較器は、周波数平均回路で算出した平均周波数
と所定の闇値との大小を比較する。
符号選択回路は、周波数比較器の比較結果が入力され、
周波数平均回路で算出した平均周波数が所定の闇値より
大きい場合にカラードプラ解析器によって解析された符
号Ｓを出力し、周波数平均回路で算出した平均周波数が
所定の闇値より小さい場合に平均符号回路で算出した平
均符号を出力する。
全体として、この符号選択回路の出力と周波数平均回路
で算出した平均周波数とに基づいて移動体の速度分布を
表示するように構成されている。
ｊ　　　量　）　１２の、−日 請求項２の超音波診断装置において、カラードプラ解析
器は、移動体から反射した超音波受信信号に基づいて移
動体の動きの方向に対応した符号Ｓ及び移動体の速度に
対応した周波数ｆを解析する。
平均符号回路は、カラードプラ解析器によって解析され
た所定期間における符号Ｓの平均を算出する。
周波数平均回路は、カラードプラ解析器によって解析さ
れた所定期間における周波数ｆの平均を算出する。
周波数比較器は、周波数平均回路で算出した平均周波数
と所定の闇値との大小を比較する。
符号選択回路は、周波数比較器の比較結果が入力され、
周波数平均回路で算出した平均周波数が所定の闇値より
大きい場合にカラードプラ解析器によって解析された符
号Ｓを出力し、周波数平均回路で算出した平均周波数が
所定の闇値より小さい場合に平均符号回路で算出した平
均符号を出力する。
全体として、この符号選択回路の出力とカラードプラ解
析器によって解析された周波数ｆとに基づいて移動体の
速度分布を表示するように構成されている。
利亘上ＪＪ尼りしλ免肌 請求項３の超音波診断装置において、カラードプラ解析
器は、移動体から反射した超音波受信信号に基づいて移
動体の動きの方向に対応した符号ｓ、移動体の速度に対
応した周波数ｆ及び超音波受信信号のパワーｐを解析す
る。
平均符号回路は、カラードプラ解析器によって解析され
た所定期間における符号Ｓの平均を算出する。
周波数平均回路は、カラードプラ解析器によって解析さ
れた所定期間における周波数ｆの平均を算出する。
パワー比較器は、カラードプラ解析器によって解析され
たパワーｐと所定の闇値との大小を比較する。
符号選択回路は、パワー比較器の比較結果が入力され、
カラードプラ解析器で解析されたパワーｐが所定の闇値
より大きい場合にカラードプラ解析器によって解析され
た符号Ｓを出力し、カラードプラ解析器で解析されたパ
ワーｐが所定の闇値より小さい場合に平均符号回路で算
出した平均符号を出力する。
全体として、この符号選択回路の出力と周波数平均回路
で算出した平均周波数とに基づいて移動体の速度分布を
表示するように構成されている。
ｉｖ　　云ン１４の　日 請求項４の超音波診断装置において、カラードプラ解析
器は、移動体から反射した超音波受信信号に基づいて移
動体の動きの方向に対応した符号ｓ、移動体の速度に対
応した周波数ｆ及び超音波受信信号のパワーｐを解析す
る。
平均符号回路は、カラードプラ解析器によって解析され
た所定期間における符号Ｓの平均を算出する。
パワー比較器は、カラードプラ解析器によって解析され
たパワーｐと所定の闇値との大小を比較する。
符号選択回路は、パワー比較器の比較結果が入力され、
カラードプラ解析器で解析されたパワーｐが所定の闇値
より大きい場合にカラードプラ解析器によって解析され
た符号Ｓを出力し、カラードプラ解析器で解析されたパ
ワーｐが所定の闇値より小さい場合に平均符号回路で算
出した平均符号を出力する。
全体として、この符号選択回路の出力とカラードプラ解
析器で解析された周波数ｒとに基づいて移動体の速度分
布を表示するように構成されている。
Ｖ　１ン１５の一■ 請求項５の超音波診断装置において、カラードプラ解析
器は、移動体から反射した超音波受信信号に基づいて移
動体の動きの方向に対応した符号ｓ、移動体の速度に対
応した周波数ｆ及び超音波受信信号のパワーｐを解析す
る。
平均符号回路は、カラードプラ解析器によって解析され
た所定期間における符号Ｓの平均を算出する。
周波数平均回路は、カラードプラ解析器によって解析さ
れた所定期間における周波数１の平均を算出する。
周波数比較器は、周波数平均回路で算出した平均周波数
と第１の闇値との大小を比較する。
パワー比較器は、カラードプラ解析器によって解析され
たパワーｐと第２の闇値との大小を比較する。
符号選択回路は、周波数比較器の比較結果とパワー比較
器の比較結果とが入力され、周波数平均回路で算出され
た平均周波数が第１の闇値より大きく、かつカラードプ
ラ解析器で解析されたパワーｐが第２の闇値より大きい
場合にカラードプラ解析器によって解析された符号Ｓを
出力し、それ以外の場合に平均符号回路で算出した平均
符号を出力する。
全体として、この符号選択回路の出力と周波数平均回路
で算出した平均周波数とに基づいて移動体の速度分布を
表示するように構成されている。
ｖｉ　　云ンゝＩ６の　■− 請求項６の超音波診断装置において、カラードプラ解析
器は、移動体から反射した超音波受信信号に基づいて移
動体の動きの方向に対応した符号ｓ、移動体の速度に対
応した周波数「及び超音波受信信号のパワーｐを解析す
る。
平均符号回路は、カラードプラ解析器によって解析され
た所定期間における符号Ｓの平均を算出する。
周波数平均回路は、カラードプラ解析器によって解析さ
れた所定期間における周波数ｆの平均を算出する。
周波数比較器は、ｒＴＪ波数平均回路で算出した平均周
波数と第１の闇値との大小を比較する。
パワー比較器は、カラードプラ解析器によって解析され
たパワーｐと第２の闇値との大小を比較する。
符号選択回路は、周波数比較器の比較結果とパワー比較
器の比較結果とが入力され、周波数平均回路で算出され
た平均周波数が第１の闇値より大きく、かつカラードプ
ラ解析器で解析されたパワーｐが第２の闇値より大きい
場合にカラードプラ解析器によって解析された符号Ｓを
出力し、それ以外の場合に平均符号回路で算出した平均
符号を出力する。
全体として、この符号選択回路の出力とカラードプラ解
析器によって解析された周波数ｆとに基づいて移動体の
速度分布を表示するように構成されている。
ｖｉ　ｆンＪ７の　口 請求項７の超音波診断装置において、カラードプラ解析
器は、移動体から反射した超音波受信信号に基づいて移
動体の動きの方向に対応した符号ｓ、移動体の速度に対
応した周波数ｆ及び超音波受信信号のパワーｐを解析す
る。
平均符号回路は、カラードプラ解析器によって解析され
た所定期間における符号Ｓの平均を算出する。
周波数平均回路は、カラードプラ解析器によって解析さ
れた所定期間における周波数ｆの平均を算出する。
周波数比較器は、周波数平均回路で算出した平均周波数
と第１の闇値との大小を比較する。
パワー比較器は、カラードプラ解析器によって解析され
たパワーｐと第２の闇値との大小を比較する。
符号選択回路は、周波数比較器の比較結果とパワー比較
器の比較結果とが入力され、周波数平均回路で算出され
た平均周波数が第１の闇値より小さく、かつカラードプ
ラ解析器で解析されたパワーｐが第２の闇値より小さい
場合に平均符号回路で算出された平均符号を出力し、そ
れ以外の場合にカラードプラ解析器によって解析された
符号Ｓを出力する。
全体として、この符号選択回路の出力と周波数平均回路
で算出した平均周波数とに基づいて移動体の速度分布を
表示するように構成されている。
Ｗ！ｊｌＵン゛工８の　日 請求項８の超音波診断装置において、カラードプラ解析
器は、移動体から反射した超音波受信信号に基づいて移
動体の動きの方向に対応した符号ｓ、移動体の速度に対
応した周波数ｆ及び超音波受信信号のパワーｐを解析す
る。
平均符号回路は、カラードプラ解析器によって解析され
た所定期間における符号Ｓの平均を算出する。
周波数平均回路は、カラードプラ解析器によって解析さ
れた所定期間における周波数ｆの平均を算出する。
周波数比較器は、周波数平均回路で算出した平均周波数
と第１の闇値との大小を比較する。
パワー比較器は、カラードプラ解析器によって解析され
たパワーｐと第２の闇値との大小を比較する。
符号選択回路は、周波数比較器の比較結果とパワー比較
器の比較結果とが入力され、周波数平均回路で算出され
た平均周波数が第１の闇値より小さく、かつカラードプ
ラ解析器で解析されたパワーｐが第２の闇値より小さい
場合に平均符号回路で算出された平均符号を出力し、そ
れ以外の場合にカラードプラ解析器によって解析された
符号Ｓを出力する。
全体として、この符号選択回路の出力とカラードプラ解
析器によって解析された周波数ｒとに基づいて移動体の
速度分布を表示するように構成されている。
〔作　用〕
＝（ｉ　　　ｉ　ン１　　の　−■ 周波数平均回路によって算出された平均周波数が所定の
闇値より大きい場合には、カラードプラ解析器によって
解析された符号Ｓを選択し、この選択した符号Ｓと周波
数平均回路で算出した平均周波数とに基づいて移動体の
速度分布の表示を行う。
一方、平均周波数が所定の闇値より小さい場合には、符
号平均回路によって算出した平均符号を選択し、この選
択した平均符号と周波数平均回路で算出した平均周波数
とに基づいて移動体の速度分布の表示を行う。
請求項１の発明にあっては、移動体の速度が遅く符号Ｓ
に含まれるノイズの割合が多い場合に、符号Ｓに代えて
平均符号が選択される。
ｉｉ　ｉン工２の　■ 周波数平均回路によって算出された平均周波数が所定の
闇値より大きい場合には、カラードプラ解析器によって
解析された符号Ｓを選択し、この選択した符号Ｓとカラ
ードプラ解析器によって解析された周波数ｆとに基づい
て移動体の速度分布の表示を行う。
一方、平均周波数が所定の闇値より小さい場合には、符
号平均回路によって算出した平均符号を選択し、この選
択した平均符号とカラードプラ解析器によって解析され
た周波数ｒとに基づいて移動体の速度分布の表示を行う
。
請求項２の発明にあっては、移動体の速度が遅く符号Ｓ
に含まれるノイズの割合が多い場合に、符号Ｓに代えて
平均符号が選択される。
ｉＩｉ　　ラ一ゝＪ３の　１１ カラードプラ解析器によって解析されたパワーｐが所定
の闇値より大きい場合には、カラードプラ解析器によっ
て解析された符号Ｓを選択し、この選択した符号Ｓと周
波数平均回路で算出した平均周波数とに基づいて移動体
の速度分布の表示を行う。
一方、パワーｐが所定の闇値より小さい場合には、符号
平均回路によって算出した平均符号を選択し、この選択
した平均符号と周波数平均回路で算出した平均周波数と
に基づいて移動体の速度分布の表示を行う。
請求項３の発明にあっては、反射してくる超音波受信信
号の強度が弱く符号Ｓに含まれるノイズの割合が多い場
合に、符号Ｓに代えて平均符号が選択される。
ｉｖ”、ｚ１４の　ロ カラードプラ解析器によって解析されたパワーｐが所定
の闇値より大きい場合には、カラードプラ解析器によっ
て解析された符号Ｓを選択し、この選択した符号Ｓとカ
ラードプラ解析器によって解析された周波数ｆとに基づ
いて移動体の速度分布の表示を行う。
一方、パワーｐが所定の闇値より小さい場合には、符号
平均回路によって算出した平均符号を選択し、この選択
した平均符号とカラードプラ解析器によって解析された
周波数ｒとに基づいて移動体の速度分布の表示を行う。
請求項４の発明にあっては、反射してくる超音波受信信
号の強度が弱く符号Ｓに含まれるノイズの割合が多い場
合に、符号Ｓに代えで平均符号が選Ｉ尺される。
Ｖ　ジンエ５の　■ 周波数平均回路によって算出された平均周波数が第１の
闇値より大きく、かつカラードプラ解析器によって解析
されたパワーｐが第２の闇値より大きい場合には、カラ
ードプラ解析器によって解析された符号Ｓを選択し、こ
の選択した符号Ｓと周波数平均回路で算出した平均周波
数とに基づいて移動体の速度分布の表示を行う。
一方、それ以外の場合には、符号平均回路によって算出
した平均符号を選択し、この選択した平均符号と周波数
平均回路で算出した平均周波数とに基づいて移動体の速
度分布の表示を行う。
請求項５の発明にあっては、移動体の速度が遅く、かつ
反射してくる超音波受信信号の強度が弱く符号Ｓに含ま
れるノイズの割合が多い場合に、符号Ｓに代えて平均符
号が選択される。
」■１箭末」１久λ班 周波数平均回路によって算出された平均周波数が第１の
闇値より大きく、かつカラードプラ解析器によって解析
されたパワーｐが第２の闇値より大きい場合には、カラ
ードプラ解析器によって解析された符号Ｓを選択し、こ
の選択した符号Ｓとカラードプラ解析器によって解析さ
れた周波数ｆとに基づいて移動体の速度分布の表示を行
なう。
一方、それ以外の場合には、符号平均回路によって算出
した平均符号を選択し、この選択した平均符号とカラー
ドプラ解析器によって解析された周波数ｒとに基づいて
移動体の速度分布の表示を行う。
請求項６の発明にあっては、移動体の速度が遅く、かつ
反射してくる超音波受信信号の強度が弱く符号Ｓに含ま
れるノイズの割合が多い場合に、符号Ｓに代えて平均符
号が選択される。
■ｉ　ｉン１の　口 周波数平均回路によって算出された平均周波数が第１の
閾値より小さく、かつカラードプラ解析器によって解析
されたパワーｐが第２の闇値より小さい場合には、符号
平均回路によって算出した平均符号を選択し、この選択
した平均符号と周波数平均回路で算出した平均周波数と
に基づいて移動体の速度分布の表示を行う。
一方、それ以外の場合には、カラードプラ解析器によっ
て解析された符号Ｓを選択し、この選択した符号Ｓと周
波数平均回路で算出した平均周波数とに基づいて移動体
の速度分布の表示を行う。
請求項７の発明にあっては、移動体の速度が遅いか、あ
るいは反射してくる超音波受信信号の強度が弱く符号Ｓ
に含まれるノイズの割合が多い場合に、符号Ｓに代えて
平均符号が選択される。
」血り許求頭ｌ傅光凱 周波数平均回路によって算出された平均周波数が第１の
闇値より小さく、かつカラードプラ解析器Ｇこよって解
析されたパワーｐが第２の闇値より小さい場合には、符
号平均回路によって算出した平均符号を選択し、この選
択した平均符号とカラードプラ解析器によって解析され
た周波数ｆとに基づいて移動体の速度分布の表示を行う
。
一方、それ以外の場合には、カラードプラ解析器によっ
て解析された符号Ｓを選択し、この選択した符号Ｓとカ
ラードプラ解析器によって解析された周波数ｆとに基づ
いて移動体の速度分布の表示を行う。
請求項８の発明にあっては、移動体の速度が遅いか、あ
るいは反射してくる超音波受信信号の強度が弱く符号Ｓ
に含まれるノイズの割合が多い場合に、符号Ｓに代えて
平均符号が選択される。
〔実施例〕
以下、図面に基づいて本発明の実施例について詳細に説
明する。
第１図に、本発明の超音波診断装置の一実施例の構成を
示す。
第１図において、１１１は超音波プローブ、１１３は送
信回路、１１５は受信回路、１２１はカラードプラ解析
器、１２３は雑音除去回路、１２５及び１４３はディジ
タル・スキャン・コンバータ、１２７及び１４５はメモ
リ、１２９及び１４７はルック・アップ・テーブル、１
３１は表示装置、１４１は振幅検出回路、１４９は合成
回路である。
超音波プローブ１１１は、例えばＰｈａｓｅｄ　Ａｒｒ
ａｙ型と称される構成を有しており、数十個程度の超音
波振動子を備えている。この超音波プローブ１１１は、
送信回路１１３及び受信回路１１５に接続されている。
送信回路１１３は、超音波プローブ１１１をドライブす
るためのものであり、超音波プローブ１１１に備わった
超音波振動子を励振することにより、所定の照射角をな
す超音波を出力させる。また、受信回路１１５は、被測
定物（移動体）によって反射された超音波が超音波プロ
ーブ１１１に入射されたときに、この反射波を検出する
ためのものであり、上述した超音波振動子を順次走査す
ることにより反射波の検出を行っている。
カラードプラ解析器１２１は、受信回路１１５の受信出
力に基づいて被測定物体の速度に関する解析を行うもの
である。解析の結果、カラードプラ解析器１．２１は、
被測定物体の速度に対応した周波数出力ｆと、この速度
の方向に対応した符号Ｓと、反射波の強度に対応したパ
ワーｐとを出力する。カラードプラ解析器の出力は雑音
除去回路１２３に入力される。
雑音除去回路１２３は、予め設定された所定のモードに
従って、カラードプラ解析器１２１から出力された符号
Ｓに補正を加えるものである。この補正は、ノイズによ
って符号が反転する現象を除去するための処理であり、
ノイズによって符号が反転する可能性がある状態を検出
したときに、符号Ｓに代えて所定時間の平均値ｇを出力
する。
この符号平均については、表示する画素の前後例えば連
続する１５画素の符号の平均をとり、その１５画素に含
まれる１つの画素の符号を決定する。
雑音除去回路１２３から出力される符号Ｓあるいはその
平均値ｉと周波数ｆはディジタル・スキャン・コンバー
タ１２５に入力される。
ディジタル・スキャン・コンバータ１２５は、雑音除去
回路１２３の出力をスキャンして順次取り込むためのも
のであり、この取り込んだデータはメモリ１２７に入力
される。このディジタル・スキャン・コンバータ１２５
によるスキャン動作は、超音波プローブ１１１から反射
波を取り込むスキャン動作に同期している。メモリ１２
７では、ディジタル・スキャン・コンバータ１２５から
出力されるデータを、被測定物体の位置に応じたアドレ
スに順次格納する。このメモリ！２７の格納データは順
次読み出されてルック・アップ・テーブル１２９に入力
される。
ルック・アップ・テーブル１２９は、メモリ１２７から
読み出されたデータに基づいて、カラーデータ（例えば
ＲＧＢデータ）を作成するためのものである。メモリ１
２７には、雑音除去回路１２３から出力された周波数ｆ
及び符号Ｓあるいはその平均値ｉが格納されており、こ
れらの格納データに基づいて被測定物体の速度に関する
カラーデータに変換する。例えば、符号Ｓあるいはその
平均値ｇが正である場合が赤色に対応しており、この赤
色の濃淡の度合いを周波数ｆに対応させたカラーデータ
への変換を行う。反対に、符号Ｓあるいはその平均値ｇ
が負である場合が青色に対応しており、この青色の濃淡
の度合いを周波数ｒに対応させたカラーデータへの変換
が行われる。この変換後のカラーデータは表示装置１３
１に入力される。
表示装置１３１は、ルック・アップ・テーブル１２９か
ら出力されるカラーデータに基づいて、被測定物体の速
度に関するカラー表示を行う。例えば、被測定物体とし
て血液を考えた場合、流れが一定方向の血流は赤色ある
いは青色の何れか一方に対応して表示され、乱流と称さ
れる高速流の場合は所々流れの方向が反転するので、青
色の所々に赤色が混じって表示されたり、赤色の所々に
青色が混じって表示されたりする。
上述した表示装置１３１による表示は血液の流れを表示
するような場合に用いられるが、超音波の反射波に基づ
いて被測定物体の形状を測定することもできる。以下、
形状を測定する動作について説明する。
振幅検出回路１４１は、受信回路１１５の受信出力の振
幅検波を行うものである。振幅検波回路１４１は、受信
信号の振幅に応じた出力、すなわち超音波プローブ１１
１に到達する反射波の強度に応じた検波出力が得られる
。この検波出力は、ディジタル・スキャン・コンバータ
１４３でスキャンされてメモリ１４５に格納され、その
後ルック・アップ・テーブル１４７に入力される。ルッ
ク・アップ・テーブル１４７は、メモリ１４５から読み
出されるデータに基づいて白黒データ（ＲＧＢが同一値
）を作成するためのものである。例えば、心臓の断面に
超音波を照射した場合には、心臓の壁面部分と内部の血
液が流れている部分とでは反射波の強度が異なるため振
幅検波回路１４Ｉの検波出力にも差異が生じ、この差異
が被測定物体の形状として表示装置１３１に表示される
。
合成回路１４９はルック・アップ・テーブル１２９．１
．４７の出力を合成して表示装置１３１に供給している
。
次に、上述した本実施例の詳細構成及び動作を第２図乃
至第１０図を参照して説明する。
第２図に、カラードプラ解析器１２１の詳細構成を示す
。２１１は直交検波回路、２１３はアナログ−ディジタ
ル（Ａ、／Ｄ）変換器、２１５はモーション・ターゲッ
ト・インジケータ（ＭＴｒ）フィルタ、２１７は自己相
関器である。
受信回路１１５の出力が直交検波回路２１１に入力され
ると、直交検波回路２１１は、照射波の出力周波数ｆＣ
と同一周波数の信号と、出力周波数ｆｃと９０度位相が
ずれた周波数の信号とによって入力に対する直交検波を
行う。この直交検波の結果、位相が９０度ずれた２つの
出力が得られ、Ａ／Ｄ変換器２１３に入力される。
Ａ、　／　Ｄ変換器２１３は、入力された２つの信号の
それぞれをディジタルデータに変換するものであり、変
換されたディジタルデータはＭＴＩフィルタ２１５に入
力される。
ＭＴＩフィルタ２１５は、Ａ／Ｄ変換器２１３から入力
される２つのディジタルデータのそれぞれに対して、バ
イパスフィルタとして動作する。
例えば、超音波プローブ１１１によって被測定物体の同
一位置に対する反射波を繰り返し受信した場合に、この
同一位置の反射波に対応するＡ／Ｄ変換器２１３の出力
データに対して、１次のバイパスフィルタとして動作す
る。従って、被測定物体が移動した場合には、この移動
速度に応じてＭＴｌフィルタ２１５の出力値が大きくな
る。ＭＴＩフィルタ２１５の２つの出力は、自己相関器
２１７に入力される。
自己相関器２１７は、ＭＴＩフィルタ２１５の２つの出
力に基づいて自己相関から周波数を求める。自己相関か
ら得られる周波数をＦ、ＭＴＩフィルタ２１５から出力
される２つの出力をＸ　（ｉ）、Ｙ（ｉ）（ｉは被測定
物体の同一位置に対して繰り返される超音波の送受信回
数）とすると、自己相関から得られる周波数Ｆは、 ・・・（１） で示される。ここで、Ｔは被測定物体の同一位置に対し
て繰り返される超音波の送信あるいは受信の周期である
。
自己相関器２１．７は、（１）式に従って求めた自己相
関による周波数Ｆの絶対値ＩＦ＋を周波数１として、自
己相関による周波数Ｆの正負を符号Ｓとして出力する。
また、自己相関器２１７は、被測定物体の同一位置に対
応する２入力Ｘ（ｉ）　、Ｙ（ｉ）の２乗和を累積した
値Σ（Ｘ　（ｉ）　２＋ｙ　（ｉ）　２）をパワーｐと
して出力する。このパワーｐは、反射波の強度に対応し
ている。自己相関器２１７のこれらの出力は、上述した
雑音除去回路１２３に入力される。
第３図に、雑音除去回路１２３の詳細構成を示す。第３
図において、３１１はパワー比較器、３１３は周波数選
択回路、３１．５は周波数平均回路、３１７は周波数比
較器、３２１は平均符号回路、３２３は等価遅延回路、
３２５は符号選択回路、３３１はモード制御回路である
。
パワー比較器３１１は、カラードプラ解析器１２１内の
自己相関器２１７から出力されたパワーｐが入力され、
このパワーｐと所定の閾値ｐ０との大小を比較するもの
である。比較結果は、周波数選択回路３１３及びモード
制御回路３３１に入力される。
周波数選択回路３１３には、カラードプラ解析器１２１
内の自己相関器２１７から出力された周波数［が入力さ
れる。周波数選択回路３１３は、パワー比較器３１１の
比較結果に応じて、この周波数ｆをそのまま出力するか
、あるいは零に置き換えて出力するかの択一的動作を行
っている。具体的には、パワー比較器３１１に入力され
たパワーＰが所定の閾値ｐ。以下のときは、周波数ｆは
ノイズの影響を大きく受けているものであるため、この
周波数ｆを零に置き換えて出力する。それ以外のときは
周波数ｆをそのまま出力する。周波数選択回路３１３の
出力は、周波数平均回路３１５に入力される。
周波数平均回路３１５は、周波数選択回路３１３の出力
に対して、所定の時間内の移動平均を求めて平均周波数
丁を出力するためのものである。
周波数ｆから平均周波数Ｔを求めることにより、画像の
相関を利用した補正を行って、後段の表示装置１３１に
よる表示を見易いものにしている。
周波数平均回路３１５の出力は、周波数比較器３１７に
入力されると共に、雑音除去回路１２３の外部に出力さ
れる。
周波数比較器３１７は、入力された平均周波数Ｔと所定
の閾値「。との大小を比較するものである。比較結果は
、モード制御回路３３１に入力される。
モード制御回路３３１は、パワー比較器３１１及び周波
数比較器３１７の各出力の少なくとも一方に基づいて、
符号選択回路３２５による選択動作を制御するものであ
る。モード制御回路３３１に入力されるモード選択信号
に応じて制御モードの設定が行われ、設定された制御モ
ードに従った制御信号が符号選択回路３２５に入力され
る。
平均符号回路３２１及び等価遅延回路３２３には、カラ
ードプラ解析器１２１内の自己相関器２１７から出力さ
れた符号Ｓが入力される。平均符号回路３２１は、この
符号Ｓを所定の時間内で平均した平均符号ｇを算出して
出力する。この平均符号は、表示する画素に対応する符
号につき、例えば連続する７画素の平均をとり、平均符
号ｉとして出力するものである。平均符号回路３２１が
ら出力された平均符号ｇは符号選択回路３２５に入力さ
れる。また、等価遅延回路３２３は、平均符号回路３２
１による平均符号算出に要する時間と同期をとるための
ものであり、入力された符号Ｓを所定時間遅延した後に
出力する。等価遅延回路３２３で遅延された符号Ｓは符
号選択回路３２５に入力される。
符号選択回路３２５は、平均符号回路３２１がら入力さ
れる平均符号ｉまたは等価遅延回路３２３から出力され
る符号Ｓの何れか一方を選択して出力する。この選択動
作は、モード制御回路３３１から入力される制御信号に
よって制御される。
符号選択回路３２５から出力される符号Ｓあるいは平均
符号ｉは、雑音除去回路１２３の外部に出力される。
第１表に、モード制御回路３３１において制御される制
御モードの詳細内容を示す。以下、各制御モードについ
て説明する。
第１表 一般に、血液が低速で流れている場合には、逆流がほと
んど存在しないことが知られている。この性質を用いた
制御モードがモード１である。
モード１に対応するモード選択信号が入力されたモード
制御回路３３１は、平均周波数Ｔが所定の閾値ｆ、以下
のときに平均符号ｇを選択し、それ以外のときに符号Ｓ
を選択する制御動作を行っている。周波数比較器３１７
が平均周波数Ｔと所定の閾値「。との比較を行っている
ので、モード制御回路３３１は、この比較結果に基づい
て符号選択回路３２５における選択動作を制御する。
このように、モード１では、雑音除去回路１２３は、被
測定物体が低速で移動している場合には符号Ｓに代えて
平均符号百を出力し、それ以外のときは符号Ｓを出力す
る。低速時には逆流がほとんど存在しないので、このよ
うな平均符号ｉへの置き換えを行っても不都合はない。
モード２に対応するモード選択信号が入力されたモード
制御回路３３１は、パワーｐが所定の闇値ｐｏ以下のと
きに平均符号ｇを選択し、それ以外のときに符号Ｓを選
択する制御動作を行っている。パワー比較器３１１がパ
ワーｐと所定の閾値ｐ０との比較を行っているので、モ
ード制御回路３３１は、この比較結果に基づいて符号選
択回路３２５における選択動作を制御する。
このように、モード２では、雑音除去回路１２３は、被
測定物体からの反射波が少なくパワーｐが小さいときに
平均符号ｉを出力し、それ以外のときは符号Ｓを出力す
る。パワーｐが小さいときは有効な反射波が得られなか
ったときに対応しているので、このような平均符号百へ
の置き換えを行っても不都合はない。
モード３に対応するモード選択信号が入力されたモード
制御回路３３１は、平均周波数Ｔが所定の閾値ｒ０以下
であり、かつパワーｐが所定の闇値２０以下のときに平
均符号ｉを選択し、それ以外のときに符号Ｓを選択する
制御動作を行っている。モード制御回路３３１は、周波
数比較器３１７及びパワー比較器３１１０両方の比較結
果に基づいて符号選択回路３２５の選択動作を制御する
。
モード４に対応するモード選択信号が入力されると、モ
ード制御Ｂ回路３３１は、平均周波数Ｔが所定の閾値ｆ
０以下であるか、あるいはパワーｐが所定の闇値２０以
下であるかの少なくとも一方が満たされたときに平均符
号ｇを選択し、それ以外のとき（平均周波数Ｔが所定の
閾値ｆ０より大きく、かつパワーｐが所定の閾値ｐ０よ
り大きいとき）に符号Ｓを選択する制御動作を行ってい
る。
モード制御回路３３１は、周波数比較器３１７及びパワ
ー比較器３１１の両方の比較結果に基づいて符号選択回
路３２５の選択動作を制御する。
このように、モード１では被測定物体が低速である場合
に、モード２では反射波が少ない場合に、モード３では
その両方を満たす場合に、モード４ではその少なくとも
一方を満たす場合に、符号Ｓに代えて平均符号ｉを選択
する。雑音除去回路１２３はこの選択した平均符号ｇと
平均周波数Ｔを出力し、後段の表示装置１３１はこれら
の各出力に基づいた表示を行う。表示装置１３１による
表示は平均符号ｇに基づいて行われるため、ノイズの影
舌が目立たない奇麗な表示が得られる。
また、表示装置１３１の表示におけるノイズを目立たせ
なくすることで、ノイズによって注意力が散漫によるこ
とがなくなり、病院等で診断を行う場合の誤診断を防止
することができる。
第４図に、平均符号回路３２１の詳細構成を示す。図に
おいて、５１１．５１２．５１３．５１４．５１５．５
１Ｇ、５１７はフリップフロップ（ＦＦ）、５２１．５
２３．５２５．５２７．５２９．５３１は加算器である
。
例えば符号Ｓは１ビツトデータであり、“ビが正に０”
が負にそれぞれ対応している。また、平均符号回路３２
１は、７サイクル分の符号Ｓを平均するものである。７
つのフリップフロップ５１１〜５１７は直列接続されて
いる。
フリップフロップ５１１及び５１２の各出力端子は加算
器５２１の２つの入力端子のそれぞれに接続されている
。フリップフロップ５１３及び５Ｉ４の各出力端子は加
算器５２３の２つの入力端子のそれぞれに接続されてい
る。フリップフロップ５１５及び５１６の各出力端子は
加算器５２７の２つの入力端子のそれぞれに接続されて
いる。
フリップフロップ５１７の出力端子及び加算器５２７の
出力端子は加算器５２９の２つの入力端子のそれぞれに
接続されている。加算器５２１及び５２３の各出力端子
は加算器５２５の２つの入力端子のそれぞれに接続され
ている。加算器５２５・及び５２９の各出力端子は加算
器５３１の２つの入力端子のそれぞれに接続されている
。
従って、フリップフロップ５１１に符号Ｓが順次入力さ
れると、加算器５３１からは、直列接続された７つのフ
リップフロップの各出力を加算した結果が得られ、その
最上位ビット（最下位ピントから４ビツト目）を平均符
号ｇとして出力する。
具体的には、加算結果が°“１０００”（十進数で４）
以上のときに平均符号ｇとして“１°′を、それ以下の
ときに平均符号ｉとして′０゛を出力する。
第５図に、等価遅延回路３２３の詳細構成を示す。図に
おいて、６１１．６１２．６１３．６１４はフリップフ
ロンプロ（ＦＦ）である。４つのフリップフロップ６１
１〜６１４は直列接続されており、フリップフロップ６
１１に入力された符号Ｓは４サイクル後にフリップフロ
ンプロ出力される。
第６図に、周波数平均回路３１５の詳細構成を示す。図
において、７１１、７１２、７】３、７１４、７１５、
７１６、７１７、７３１はフリップフロップ（ＦＦ）、
７２１及び７２３は加算器、７４１は除算器である。
周波数平均回路３１５に入力された周波数ｒは、フリッ
プフロップ７１１及び加算器７２１の一方の入力端子に
入力される。フリップフロップ７１１〜７１７は直列接
続されており、フリップフロップ７１７の出力は加算器
７２１の他方の入力端子に入力される。この他方の入力
端子は反転入力端子であり、加算器７２１は、一方の入
力端子に入力された周波数ｒから他方の入力端子に入力
されたフリップフロップ７１７の出力を減算して出力す
る。
加算器７２１の出力端子は加算器７２３の一方の入力端
子に接続されており、加算器７２３の出力端子は除算器
７４１の入力端子に接続されていると共に、フリップフ
ロップ７３１を介して自分自身の他方の入力端子に接続
されている。加算器７２３の出力はフリップフロップ７
３１にーら：；１格納された後、再度加算器７２３自身
に入力されるため、加算器７２３では入力される加算器
７２１の出力を積算することになる。このような構成に
よって、常にフリップフロップ７１１〜７１７の個数に
応じたサイクル分の周波数「の積算値を求めることがで
きる。
除算器７４１は、加算３７２３の出力をフリップフロッ
プ７１１〜７１７の個数に応じた数（この場合は７）で
除算し、除算結果として平均周波数Ｔを出力する。
なお、フリップフロップ７１１〜７１７の各出力を加算
し、この加算結果を所定の値で除算して平均周波数Ｔを
求めるようにしてもよいが、第６図に示す構成を用いる
ことにより、使用する加算器の数を減らすことができる
。また、直列接続するフリップフロップの数を２の倍数
に設定すると、除算器７４１における除算処理をビット
操作で行うことができるので、更に構成を簡略化するこ
とが可能になる。
第７図に、モード制御回路３３１の詳細構成を示す。図
において、８１１は選択回路、８２１はＡＮＤ回路、８
２３はＯＲ回路である。
周波数比較器３１７からは周波数ｆが所定の闇値ｆ０以
下のときにその旨を示す“１”が出力され、それ以外の
ときは“°０パが出力されるものとする。同様に、パワ
ー比較器３１２からはパワーｐが所定の閾値ｐ０以下の
ときにその旨を示す“１゛が、それ以外のときは“′０
°゛が出力されるものとする。
周波数比較器３１７の出力は、選択回路８１１の第１入
力端子に入力されると共に、Ａ．　Ｎ　Ｄ回路８２１の
一方の入力端子及びＯＲ回路８２３の一方の入力端子に
入力される。また、パワー比較器３１１の出力は、選択
回路８１１の第２入力端子に入力されると共に、ＡＮＤ
回路８２１の他方の入力端子及びＯＲ回路８２３の他方
の入力端子に入力される。更に、ＡＮＤ回路８２１の出
力は選択回路８１１の第３入力端子に入力され、ＯＲ回
路８２３の出力は選択回路８１１の第４入力端子に入力
される。
選択回路８１１は、上述した４つの入力端子の他にモー
ド選択信号を入力するだめの制御端子ををしている。選
択回路８１１は、モード１を示すモード選択信号が制御
端子に入力されると、第１入力端子に入力された信号を
選択して出力する。
同様に、モード２を示すモード選択信号が制御端子に入
力されると、第２入力端子に入力された信号を選択して
出力する。モード３を示すモード選択信号が制御端子に
入力されると、第３入力端子に入力された信号を選択し
て出力する。モード４を示すモード選択信号が制御端子
に入力されると、第４入力端子に入力された信号を選択
して出力する。
符号選択回路３２５は、上述した選択回路８１１の出力
に応じた選択動作を行っており、選択回路８１１の出力
が”　１　”であるときは平均符号回路３２１から出力
される平均符号ｇを選択して、それ以外のときは等価遅
延回路３２３から出力される符号Ｓを選択して出力する
。
第８図に、符号選択回路３２５及びモード制御回路３３
１内の選択回路８１１を共にＴＩ社製標準ロジック５Ｎ
７Ｊ１．５１５３（Ｄａｔａ　５ｅｌｅｃｔｏｒ）を用
いて実現した場合の構成を示す。この標準ロジックは、
４つの入力端子Ｃ０−Ｃ５を、入力端子Ａ及びＢの入力
の論理に従って選択することができる。入力端子Ａ及び
Ｂには、モード１〜４に対応して“００”から“１１”
°までの各２ビツトが入力される。なお、符号選択回路
３２５では２つの入力に対する選択を行うので、上述し
た４つの入力端子の中のＣＯ及びＣＩを使用し、入力端
子への入力論理に従って何れかを選択する（入力端子Ｂ
には固定的にローレベルの信号を入力しておく）。
また、パワー比較器３１１及び周波数比較器３１７は、
Ｔ１社製標準ロジック５Ｎ７４ＬＳ６８２〜６８５（Ｍ
ａｇｎｉｔｕｄｅ　Ｃｏｍｐａｒａｔｏｒ）を用いて実
現することができる。この場合、Ｐを周波数またはパワ
ーの入力値、Ｑを比較設定値（閾値ｐ０またはｆ、）と
し、Ｐ＞Ｑ（ＩＰＩＮ）を出力として得ることができる
。
第９図に、雑音除去回路１２３の別の構成例を示す。第
３図に示した構成が周波数平均回路３１５及び符号選択
回路３２５の出力を雑音除去回路１２３の出力としてい
たのに対し、第９図に示した構成は周波数選択回路３１
３及び符号選択回路３２５の出力を雑音除去回路１２３
の出力としている点で異なる。
第９図に示した構成では、周波数選択回路３１３から直
接的に出力を得ているため、出力される周波数ｒ内にノ
イズ成分が混じっている場合があるが、パワーｐが所定
の閾値ｐ０より大きな場合に対応しているためＳＮ比が
大きく、実際上差し支えない。従って、モード１〜モー
ド４のそれぞれにおいて平均符号ｇを選択し、この平均
符号ｇと周波数ｆに基づいてカラー表示を行う場合に、
ノイズの影響を除去して奇麗な画像を表示することがで
きる。また、病院等における誤診断を防止することがで
きる。
第１０図に、雑音除去回路１２３の更に別の構成例を示
す。第１０図に示した構成は、第３図の周波数比較器３
１７及びモード制御回路３３１を符号１１Ｊｅ論理ＲＯ
Ｍ３４１に置き換えたものである。この符号選択論理Ｒ
ＯＭ３４１に第１表に示した選択情報を格納しておいて
、各モードの指定及び符号の選択に必要な各種データを
アドレスとして入力し、対応する選択信号を出力する。
例えば、周波数平均回路３１５が出力する平均周波数Ｔ
が３ビツトデータとすると、この３ビツトデータとパワ
ー比較器３１１の比較結果に対応した１ビツトデータと
モート選択を指示する２ビツトデータの合計６ビツトデ
ータをアドレスとして符号選択論理ＲＯＭ３４１に入力
し、平均符号ｉあるいは符号Ｓの選択を指示する１ビツ
トデータを出力する。上述した符号選択論理ＲＯＭ３４
１はＲＡＭ等の他の半導体メモリで構成してもよい。
なお、上述した実施例では、モード制御回路３３１ある
いは符号選択論理ＲＯＭ３４１によってモード１〜４の
何れかを選択するようにして全モードに対応するように
したが、何れかのモードのみを固定的に実施する場合に
はモード制御回路３３１あるいは符号選択論理ＲＯＭ３
４１を省略してもよい。この場合、例えば固定的にモー
ド１の制御を行うには、周波数比較２Ｓ３１７の出力を
符号選択回路３２５に直接入力すればよい。また、固定
的にモード２の制御を行うには、パワー比較器３１１の
出力を符号選択回路３２５に直接入力すればよい。固定
的にモード３の制御を行うには、パワー比較器３１１の
出力と周波数比較器３１７の出力との論理積演算を行っ
た結果を符号選択回路３２５に入力すればよい。更に、
固定的にモード４の制御を行うには、パワー比較器３１
１の出力と周波数比較器３１７の出力との論理和演算を
行った結果を符号選択回路３２５に入力すればよい。
また、第１０図に示した構成の雑音除去回路１２３では
周波数平均回路３１５から出力される平均周波数Ｔを外
部に出力しているが、第９図と同様にして周波数選択回
路３１３の出力である周波数ｆを外部に出力するように
してもよい。
〔発明の効果］ 請求項１及び請求項２の発明によれば、移動体の速度が
遅く符号Ｓに含まれるノイズの割合が多い場合に、符号
Ｓに代えて平均符号を用いて移動体の速度分布の表示を
行うことにより、ノイズによる画像の乱れを防止するこ
とができる。
請求項３及び請求項４の発明によれば、反射してくる超
音波受信信号の強度が弱く符号Ｓに含まれるノイズの割
合が多い場合に、符号Ｓに代えて平均符号を用いて移動
体の速度分布の表示を行うことにより、ノイズによる画
像の乱れを防止することができる。
請求項５及び請求項６の発明にあっては、移動体の速度
が遅く、かつ反射してくる超音波受信信号の強度が弱く
符号Ｓに含まれるノイズの割合が多い場合に、符号Ｓに
代えて平均符号を用いて移動体の速度分布の表示を行う
ことにより、ノイズによる画像の乱れを防止することが
できる。
請求項７及び請求項８の発明によれば、移動体の速度が
遅いか、あるいは反射してくる超音波受信信号の強度が
弱く符号Ｓに含まれるノイズの割合が多い場合に、符号
Ｓに代えて平均符号を用いて移動体の速度分布の表示を
行うことにより、ノイズによる画像の乱れを防止するこ
とができる。
【図面の簡単な説明】
第１図は本発明の超音波診断装置の一実施例の構成を示
す図、 第２図は本実施例におけるカラードプラ解析器の詳細構
成を示す図、 第３図は本実施例における雑音除去回路の詳細構成を示
す図、 第４図は本実施例における平均符号回路の詳細構成を示
す図、 第５図は本実施例における等価遅延回路の詳細構成を示
す図、 第６図は本実施例における周波数平均回路の詳細構成を
示す図、 第７図は本実施例におけるモード制御回路の詳細構成を
示す図、 第８図は実施例における符号選択回路及びモード制御回
路内の選択回路の具体例の構成を示す図、第９図は本発
明の超音波診断装置の別の実施例の構成を示す図、 第１０図は本発明の超音波診断装置の更に別の実施例の
構成を示す図である。 図において、 １２１はカラードプラ解析器、 ３１１はパワー比較器、 ３１３は周波数選択回路、 ３１５は周波数平均回路、 ３１７は周波数比較器、 ３２１は平均符号回路、 ３２３は等価遅延回路、 ３２５は符号選択回路、 ３３１はモード制御回路である。 刀う−ドプラ角早才汀器／）詩阜咄或乞灯図第２図 事冷号回ヱｈ詳細構へ乞利図 第４図 慣音顯ｑ冒可語のＬ髄田矛動べＥ示ず図第３図 等１珀ノを起Ｉり路の６子細構メ濾εが１υｄ第５図 モービ制祥７７区は各め詳細Ｉ青へ↓胴囲第７図 絢’ｊｄｈ欠回詮及−ｒ”七子１り淳Ｖ口路円のユｉ巳
丁に回ヱ昏の」工１ネｆクリの横＊ｈ、示■図第８図 ン釦銹色）りＪの４１４〆”＜Ｌｆｉミず　μ円第９図

Claims (8)

    【特許請求の範囲】
  1. （１）移動体から反射した超音波受信信号に基づいて前
    記移動体の動きの方向に対応した符号ｓ及び前記移動体
    の速度に対応した周波数ｆを解析するカラードプラ解析
    器と、 前記カラードプラ解析器によって解析された所定期間に
    おける符号ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によって解析された所定期間に
    おける周波数ｆの平均を算出する周波数平均回路と、 前記周波数平均回路で算出した平均周波数と所定の閾値
    との大小を比較する周波数比較器と、前記周波数比較器
    の比較結果が入力され、前記周波数平均回路で算出した
    平均周波数が前記所定の閾値より大きい場合に前記カラ
    ードプラ解析器によって解析された符号ｓを出力し、前
    記周波数平均回路で算出した平均周波数が前記所定の閾
    値より小さい場合に前記平均符号回路で算出した平均符
    号を出力する符号選択回路と、 を備え、この符号選択回路の出力と前記周波数平均回路
    で算出した平均周波数とに基づいて前記移動体の速度分
    布を表示することを特徴とする超音波診断装置。
  2. （２）移動体から反射した超音波受信信号に基づいて前
    記移動体の動きの方向に対応した符号ｓ及び前記移動体
    の速度に対応した周波数ｆを解析するカラードプラ解析
    器と、 前記カラードプラ解析器によって解析された所定期間に
    おける符号ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によって解析された所定期間に
    おける周波数ｆの平均を算出する周波数平均回路と、 前記周波数平均回路で算出した平均周波数と所定の閾値
    との大小を比較する周波数比較器と、前記周波数比較器
    の比較結果が入力され、前記周波数平均回路で算出した
    平均周波数が前記所定の閾値より大きい場合に前記カラ
    ードプラ解析器によって解析された符号ｓを出力し、前
    記周波数平均回路で算出した平均周波数が前記所定の閾
    値より小さい場合に前記平均符号回路で算出した平均符
    号を出力する符号選択回路と、 を備え、この符号選択回路の出力と前記カラードプラ解
    析器によって解析された周波数ｆとに基づいて前記移動
    体の速度分布を表示することを特徴とする超音波診断装
    置。
  3. （３）移動体から反射した超音波受信信号に基づいて前
    記移動体の動きの方向に対応した符号ｓ、前記移動体の
    速度に対応した周波数ｆ及び前記超音波受信信号のパワ
    ーｐを解析するカラードプラ解析器と、 前記カラードプラ解析器によって解析された所定期間に
    おける符号ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によって解析された所定期間に
    おける周波数ｆの平均を算出する周波数平均回路と、 前記カラードプラ解析器によって解析されたパワーｐと
    所定の閾値との大小を比較するパワー比較器と、 前記パワー比較器の比較結果が入力され、前記カラード
    プラ解析器で解析されたパワーｐが前記所定の閾値より
    大きい場合に前記カラードプラ解析器によって解析され
    た符号ｓを出力し、前記カラードプラ解析器で解析され
    たパワーｐが前記所定の閾値より小さい場合に前記平均
    符号回路で算出した平均符号を出力する符号選択回路と
    、を備え、この符号選択回路の出力と前記周波数平均回
    路で算出した平均周波数とに基づいて前記移動体の速度
    分布を表示することを特徴とする超音波診断装置。
  4. （４）移動体から反射した超音波受信信号に基づいて前
    記移動体の動きの方向に対応した符号ｓ、前記移動体の
    速度に対応した周波数ｆ及び前記超音波受信信号のパワ
    ーｐを解析するカラードプラ解析器と、 前記カラードプラ解析器によって解析された所定期間に
    おける符号ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によって解析されたパワーｐと
    所定の閾値との大小を比較するパワー比較器と、 前記パワー比較器の比較結果が入力され、前記カラード
    プラ解析器で解析されたパワーｐが前記所定の閾値より
    大きい場合に前記カラードプラ解析器によって解析され
    た符号ｓを出力し、前記カラードプラ解析器で解析され
    たパワーｐが前記所定の閾値より小さい場合に前記平均
    符号回路で算出した平均符号を出力する符号選択回路と
    、を備え、この符号選択回路の出力と前記カラードプラ
    解析器で解析された周波数ｆとに基づいて前記移動体の
    速度分布を表示することを特徴とする超音波診断装置。
  5. （５）移動体から反射した超音波受信信号に基づいて前
    記移動体の動きの方向に対応した符号ｓ、前記移動体の
    速度に対応した周波数ｆ及び前記超音波受信信号のパワ
    ーｐを解析するカラードプラ解析器と、 前記カラードプラ解析器によって解析された所定期間に
    おける符号ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によって解析された所定期間に
    おける周波数ｆの平均を算出する周波数平均回路と、 前記周波数平均回路で算出した平均周波数と第１の閾値
    との大小を比較する周波数比較器と、前記カラードプラ
    解析器によって解析されたパワーｐと第２の閾値との大
    小を比較するパワー比較器と、 前記周波数比較器の比較結果と前記パワー比較器の比較
    結果とが入力され、前記周波数平均回路で算出された平
    均周波数が前記第１の閾値より大きく、かつ前記カラー
    ドプラ解析器で解析されたパワーｐが前記第２の閾値よ
    り大きい場合に前記カラードプラ解析器によって解析さ
    れた符号ｓを出力し、それ以外の場合に前記平均符号回
    路で算出した平均符号を出力する符号選択回路と、を備
    え、この符号選択回路の出力と前記周波数平均回路で算
    出した平均周波数とに基づいて前記移動体の速度分布を
    表示することを特徴とする超音波診断装置。
  6. （６）移動体から反射した超音波受信信号に基づいて前
    記移動体の動きの方向に対応した符号ｓ、前記移動体の
    速度に対応した周波数ｆ及び前記超音波受信信号のパワ
    ーｐを解析するカラードプラ解析器と、 前記カラードプラ解析器によって解析された所定期間に
    おける符号ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によって解析された所定期間に
    おける周波数ｆの平均を算出する周波数平均回路と、 前記周波数平均回路で算出した平均周波数と第１の閾値
    との大小を比較する周波数比較器と、前記カラードプラ
    解析器によって解析されたパワーｐと第２の閾値との大
    小を比較するパワー比較器と、 前記周波数比較器の比較結果と前記パワー比較器の比較
    結果とが入力され、前記周波数平均回路で算出された平
    均周波数が前記第１の闇値より大きく、かつ前記カラー
    ドプラ解析器で解析されたパワーｐが前記第２の閾値よ
    り大きい場合に前記カラードプラ解析器によって解析さ
    れた符号ｓを出力し、それ以外の場合に前記平均符号回
    路で算出した平均符号を出力する符号選択回路と、を備
    え、この符号選択回路の出力と前記カラードプラ解析器
    によって解析された周波数ｆとに基づいて前記移動体の
    速度分布を表示することを特徴とする超音波診断装置。
  7. （７）移動体から反射した超音波受信信号に基づいて前
    記移動体の動きの方向に対応した符号ｓ、前記移動体の
    速度に対応した周波数ｆ及び前記超音波受信信号のパワ
    ーｐを解析するカラードプラ解析器と、 前記カラードプラ解析器によって解析された所定期間に
    おける符号ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によって解析された所定期間に
    おける周波数ｆの平均を算出する周波数平均回路と、 前記周波数平均回路で算出した平均周波数と第１の閾値
    との大小を比較する周波数比較器と、前記カラードプラ
    解析器によって解析されたパワーｐと第２の閾値との大
    小を比較するパワー比較器と、 前記周波数比較器の比較結果と前記パワー比較器の比較
    結果とが入力され、前記周波数平均回路で算出された平
    均周波数が前記第１の閾値より小さく、かつ前記カラー
    ドプラ解析器で解析されたパワーｐが前記第２の閾値よ
    り小さい場合に前記平均符号回路で算出された平均符号
    を出力し、それ以外の場合に前記カラードプラ解析器に
    よって解析された符号ｓを出力する符号選択回路と、を
    備え、この符号選択回路の出力と前記周波数平均回路で
    算出した平均周波数とに基づいて前記移動体の速度分布
    を表示することを特徴とする超音波診断装置。
  8. （８）移動体から反射した超音波受信信号に基づいて前
    記移動体の動きの方向に対応した符号ｓ、前記移動体の
    速度に対応した周波数ｆ及び前記超音波受信信号のパワ
    ーｐを解析するカラードプラ解析器と、 前記カラードプラ解析器によって解析された所定期間に
    おける符号ｓの平均を算出する平均符号回路と、 前記カラードプラ解析器によって解析された所定期間に
    おける周波数ｆの平均を算出する周波数平均回路と、 前記周波数平均回路で算出した平均周波数と第１の閾値
    との大小を比較する周波数比較器と、前記カラードプラ
    解析器によって解析されたパワーｐと第２の閾値との大
    小を比較するパワー比較器と、 前記周波数比較器の比較結果と前記パワー比較器の比較
    結果とが入力され、前記周波数平均回路で算出された平
    均周波数が前記第１の閾値より小さく、かつ前記カラー
    ドプラ解析器で解析されたパワーｐが前記第２の閾値よ
    り小さい場合に前記平均符号回路で算出された平均符号
    を出力し、それ以外の場合に前記カラードプラ解析器に
    よって解析された符号ｓを出力する符号選択回路と、を
    備え、この符号選択回路の出力と前記カラードプラ解析
    器によって解析された周波数ｆとに基づいて前記移動体
    の速度分布を表示することを特徴とする超音波診断装置
    。
JP2007986A 1989-01-17 1990-01-16 超音波診断装置 Granted JPH02289237A (ja)

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 true JPH02289237A (ja) 1990-11-29
JPH0572817B2 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)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1993016399A1 (en) * 1992-02-06 1993-08-19 Furuno Electric Company, Limited Underwater detector
JP2009261436A (ja) * 2008-04-22 2009-11-12 Aloka Co Ltd 超音波診断装置

Families Citing this family (9)

* Cited by examiner, † Cited by third party
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 株式会社東芝 超音波診断装置
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 深圳迈瑞生物医疗电子股份有限公司 消除彩色血流图像中速度异常点的方法
CN102429684B (zh) 2010-09-28 2013-10-09 深圳迈瑞生物医疗电子股份有限公司 一种多普勒彩色血流成像方法和装置

Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 超音波血流画像形成装置
EP0163664A4 (en) * 1983-11-10 1987-01-10 Acoustec Partners ULTRASONIC DIAGNOSTIC DEVICE.
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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 超音波診断装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1993016399A1 (en) * 1992-02-06 1993-08-19 Furuno Electric Company, Limited Underwater detector
JP2009261436A (ja) * 2008-04-22 2009-11-12 Aloka Co Ltd 超音波診断装置

Also Published As

Publication number Publication date
EP0379139B1 (en) 1994-07-27
JPH0572817B2 (ja) 1993-10-13
US5042491A (en) 1991-08-27
DE69010923T2 (de) 1994-11-10
DE69010923D1 (de) 1994-09-01
EP0379139A1 (en) 1990-07-25

Similar Documents

Publication Publication Date Title
US7713204B2 (en) Image data processing method and apparatus for ultrasonic diagnostic apparatus, and image processing apparatus
US5931784A (en) Ultrasonic diagnostic apparatus
JP3406106B2 (ja) 超音波画像表示方法および超音波診断装置
EP1609423B1 (en) Ultrasonic diagnostic apparatus
US5105817A (en) Ultrasonic bloodstream imaging apparatus
EP2982305A1 (en) Diagnostic ultrasound apparatus and elasticity evaluation method
JPS60119929A (ja) 超音波診断装置
KR100742467B1 (ko) 촬상 시스템, 촬상 방법 및 혈액 움직임 촬상 시스템
US5188113A (en) Ultrasonic diagnosis apparatus
JP2001178720A (ja) 連続データ獲得を用いた超音波フロー・イメージングにおける運動の可視化のための方法および装置
US8708911B2 (en) Ultrasonic imaging apparatus and a method for displaying diagnostic images
JPH02289237A (ja) 超音波診断装置
US6402694B1 (en) Color imaging display method and apparatus for color doppler imaging system
JPH0614930B2 (ja) 超音波診断装置
JPH0548692B2 (ja)
US5072734A (en) Pulse doppler mti system
US20070161900A1 (en) Ultrasonic doppler blood flow measuring device
CN112672694B (zh) 确定超声波发射角度的方法以及超声设备
JP6677042B2 (ja) 画像処理装置、超音波診断装置及び画像処理プログラム
JP3260037B2 (ja) 超音波ドプラ血流計
JP2950600B2 (ja) カラードプラ超音波診断装置
JP2008173448A (ja) 超音波ドプラ診断装置
JP2602898B2 (ja) 超音波血流イメージング装置
JP2565917B2 (ja) 画像処理表示装置
JPH03162841A (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