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

Abstract

(57)【要約】 【課題】 ノイズの影響を少なくして最高周波数を検出
することでその検出精度を向上させ、操作者の負担を軽
減でき、さらに、診断時間を著しく短縮できる超音波診
断装置を提供する。 【解決手段】 被検体内に超音波を送受波し、これによ
り得られた受信信号からドプラ信号を抽出し、このドプ
ラ信号を処理して周波数スペクトル信号を得、この周波
数スペクトル信号を表示するとともに、前記周波数スペ
クトル信号の各周波数についての最高周波数を自動的に
検出する超音波診断装置において、前記周波数スペクト
ルの各周波数の強度を示すデータ値を検索して閾値と比
較し、閾値を所定の周波数幅中で所定数越えた時、その
最初に越えた位置を最高周波数とすることを特徴とす
る。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、超音波のドプラ効
果を利用して血液等、被検体内の運動状態を診断する超
音波診断装置に係り、特に高精度でドプラ偏移の最高周
波数を検出する超音波診断装置に関する。

【０００２】

【従来の技術】従来より、超音波パルスドプラ法とパル
ス反射法とを併用することによって一つの超音波プロー
ブで血流情報と断層像（白黒Ｂモード像）情報を得、血
流情報をリアルタイムで表示するようにした超音波診断
装置が知られている。

【０００３】図１１は、従来の超音波診断装置の機能ブ
ロック図を示している。図１１に示すように、従来の超
音波診断装置１００は、超音波振動子１０１と、パルサ
１０３と、前置増幅器１０５と、ミキサ１０７と、ロー
パスフィルタ１０９と、サンプルホールド回路１１１
と、バンドパスフィルタ１１３と、高速フーリエ変換器
（以下ＦＦＴと記す）１１５と、表示器１１７と、レン
ジゲート回路１１９と、パルス発生回路１２１とを備え
ている。

【０００４】次に、従来の超音波診断装置１００を用い
て血流速度を測定する場合の動作原理を図１１と図１２
を用いて説明する。

【０００５】まず、パルス発生回路１２１は、クロック
パルスａを発振し、このクロックパルスａをレンジゲー
ト回路１１９に出力する。また、パルス発生回路１２１
は、分周回路およびゲート回路等を備え、このクロック
パルスａに基づき超音波繰り返し周波数に相当するレー
トパルスｂを発生し、パルサ１０３およびレンジゲート
回路１１９に出力する。

【０００６】パルサ１０３は、供給されたレートパルス
ｂから高電圧の駆動パルスを作り、超音波振動子１０１
を駆動する。超音波振動子１０１は、電気信号を機械振
動に変換し、生体表面Ｐを介して生体内へ超音波を送波
する。この超音波は生体内の血管壁および血管内の血流
（主に赤血球）により一部反射され、そのエコー信号は
同一の超音波振動子１０１で受信され、電気信号ｄに変
換される。

【０００７】すなわち、被検体である生体内を流れてい
る血流に対して超音波パルスを送波すると、この超音波
ビームの中心周波数ｆｃは、流動する血流によって散乱
され、ドプラ偏移を受けて周波数ｆｄだけ変化する。こ
のため、受信周波数ｆはｆ＝ｆｃ＋ｆｄとなる。このと
き、周波数ｆｃ，ｆｄは、次の式（１）により現され
る。 ｆｄ＝（２・Ｖ・ｃｏｓθ／Ｃ）・ｆｃ …（１） ここで、Ｖは血流速度、θは超音波ビームと血流との成
す角度、Ｃは生体内の音速（約１５３０ｍ／Ｓ）を表し
ている。

【０００８】したがって、ドプラ偏移ｆｄを検出するこ
とによって血流速度Ｖを得ることができる。

【０００９】前置増幅器１０５は、超音波振動子１０１
から入力する前記電気信号ｄを増幅した後、ミキサ１０
７に入力する。ミキサ１０７は前置増幅器１０５から入
力する電気信号ｄとクロックパルス発生部１２１から入
力するクロックパルスａとを混合する。そしてローパス
フィルタ１０９は、ミキサ１０７から入力する混合信号
の内、超音波周波数等の高周波成分を除去する。

【００１０】そして生体内の血流が流れている深さの位
置Ｏだけのドプラ信号を抽出するために、ローパスフィ
ルタ１０９からの信号をサンプルホールド回路１１１に
出力する。

【００１１】レンジゲート回路１１９は、遅延時間が任
意に設定できるが、ここでは、超音波が超音波振動子１
０１からレンジゲート位置（サンプリングポイント位
置）Ｏまでの往復する時間だけ、信号ｄより遅延させ、
設定された長さに対応する幅を有するサンプリングパル
スｃをサンプルホールド回路１１１に加える。

【００１２】前記レンジゲート位置Ｏは操作者が超音波
断層像中の血流信号を得たい血管に、トラックボールや
ジョイスティック等のポインティングデバイスを用いて
自由に位置を合わせることが可能である。

【００１３】サンプルホールド回路１１１は、サンプリ
ングパルスｃによりローパスフィルタ１０９の出力信号
をサンプルホールドする。バンドパスフィルタ１１３
は、サンプルホールド回路１１１でのサンプリングによ
り生じた高周波成分並びに血管等の固定反射信号または
比較的ゆっくりした動きによるドプラ偏移周波数を除去
し、血流によるドプラ信号のみを抽出する。次に周波数
分析器としてのＦＦＴ１１５を介して周波数スペクトル
パターンを表示器１１７に表示する。

【００１４】このようにして、超音波診断装置１００は
血流速度に対応するドプラ偏移周波数を検出、表示する
とともに、最高周波数を検出して表示している。

【００１５】次に、従来の最高周波数の検出アルゴリズ
ムを図１３に示す周波数スペクトルが得られた場合を例
に説明する。従来の超音波診断装置１００は、図１３
（ａ），（ｂ）に示すようにパワー値のピーク値を基準
にし、このピーク値より所定の比率低い値を閾値とし、
この閾値を最初に越える所のドプラ偏移周波数を最高周
波数としている。

【００１６】なお、最高周波数とは、流れの方向（順流
方向（超音波プローブに近付く方向）の場合は図１３
中、（＋）ｆ側、逆流方向（超音波プローブから遠ざか
る方向）の場合は図１３中、（−）ｆ側）から０Hz側へ
向かってパワー値をサーチしていき、パワー値が閾値を
最初に越えた所のドプラ偏移周波数としている。また、
前記流れの方向は、周波数スペクトラムにて最大周波数
ｆpeakまたは平均周波数ｆmeanがプラス側（（＋）ｆ
側）にあれば順流方向、マイナス側（（−）ｆ側）にあ
れば逆流方向となる。

【００１７】

【発明が解決しようとする課題】しかしながら、従来の
超音波診断装置の最高周波数の検出アルゴリズムでは、
閾値を比率Ａに設定した場合、図１３（ａ）に示すよう
に低Ｓ／Ｎ比のときに検出される最高周波数（ｆmax-A
）は良いが、図１３（ｂ）に示すように高Ｓ／Ｎ比の
ときに検出される最高周波数（ｆmax-A ）は低めの値と
なり、閾値を比率Ｂに設定した場合、図１３（ｂ）に示
すように高Ｓ／Ｎ比のときに検出される最高周波数（ｆ
max-B ）は良いが、図１３（ａ）に示すように低Ｓ／Ｎ
比のときに検出される最高周波数（ｆmax-B ）は高めの
値になるという問題がある。

【００１８】また、従来の超音波診断装置の最高周波数
の検出アルゴリズムでは、閾値を最初に越える所のドプ
ラ偏移周波数を最高周波数としているので、ノイズ成分
でも最高周波数と検出する場合があった。

【００１９】さらに、従来の超音波診断装置の最高周波
数の検出アルゴリズムでは、平均処理をしていない周波
数スペクトルを用いて最高周波数の検出を行っているた
め、ノイズ成分でも最高周波数と検出する場合があっ
た。

【００２０】このため、正確な最高周波数の検出に時間
が掛かり、操作者の大きな負担となっていたと同時に、
診断に非常に多くの時間が掛かっていた。

【００２１】本発明は上記事情に鑑みて成されたもので
あり、その目的は、ノイズの影響を少なくして最高周波
数を検出することでその検出精度を向上させ、操作者の
負担を軽減でき、さらに、診断時間を著しく短縮できる
超音波診断装置を提供することにある。

【００２２】また、本発明の別の目的は、ノイズ等の影
響により少しの期間だけ最高周波数の時間経過に伴う推
移（最高周波数のトレース）を誤った場合に、操作者の
介入によりこの誤りを修正することが可能な超音波診断
装置を提供することにある。

【００２３】

【課題を解決するための手段】上記の目的を達成するた
めに本発明は、被検体内に超音波を送受波し、これによ
り得られた受信信号からドプラ信号を抽出し、このドプ
ラ信号を処理して周波数スペクトラムを得、この周波数
スペクトラムを表示するとともに、前記周波数スペクト
ラムの各周波数についての最高周波数を検出する超音波
診断装置において、前記周波数スペクトラムの各周波数
の信号強度を示すデータ値について、流れ方向側端部か
ら周波数０Ｈｚの方向に向かって順に検索し、該データ
値が予め設定された閾値を所定の周波数幅中で所定回数
以上越えた時、その最初に越えた位置を最高周波数とす
るか、或いは最も信号強度の大きい周波数ないし平均周
波数から流れ方向側端部の方向に向かって順に検索し、
該データ値が予め設定された閾値を所定の周波数幅中で
所定回数以上越えなかった時、その最初に越えなかった
位置を最高周波数とすることを要旨とする。

【００２４】すなわち、周波数スペクトラムを流れ方向
側端部から周波数０Ｈｚに向かって順に検索する場合、
検索開始当初は各データ値は予め設定された閾値より信
号強度が低く、ノイズによる不規則な変動を伴って徐々
に信号強度が高くなると考えられる。そして各周波数の
信号強度を示すデータ値が所定の周波数幅中で所定回数
以上閾値を越える時、その最初に閾値を越える位置を最
高周波数とする。

【００２５】これとは逆に、周波数スペクトラムを最も
信号強度の大きい周波数ないし平均周波数から流れ方向
側端部に向かって順に検索する場合、検索開始当初は各
データ値は予め設定された閾値より信号強度が高く、ノ
イズによる不規則な変動を伴って徐々に信号強度が低く
なると考えられる。そして各周波数の信号強度を示すデ
ータ値が所定の周波数幅中で所定回数以上閾値を越えな
い時、その最初に閾値を越えない位置を最高周波数とす
る。

【００２６】また本発明は、前記周波数スペクトラムを
時間方向、周波数方向の少なくとも一方向について平均
処理する平均処理手段と、前記平均処理手段により平均
処理された周波数スペクトラムの各周波数の信号強度を
示すデータ値について、流れ方向側端部から周波数０Ｈ
ｚの方向に向かって順に検索し、該データ値が予め設定
された閾値を所定の周波数幅中で所定回数以上越えた
時、その最初に越えた位置を最高周波数とするか、或い
は最も信号強度の大きい周波数ないし平均周波数から流
れ方向側端部の方向に向かって順に検索し、該データ値
が予め設定された閾値を所定の周波数幅中で所定回数以
上越えなかった時、その最初に越えなかった位置を最高
周波数とする最高周波数抽出手段と、を備えることを要
旨とする。

【００２７】また本発明は、前記周波数スペクトラムを
時間方向、周波数方向の少なくとも一方向について平均
処理する平均処理手段と、前記平均処理手段により平均
処理された周波数スペクトラムの各周波数について、そ
の信号値とノイズ値とを基に、最高周波数の閾値を演算
する閾値演算手段とをさらに備え、前記平均処理手段に
より平均処理された周波数スペクトラムに対して、前記
演算された閾値を用いて最高周波数を検索することがで
きる。

【００２８】さらに、前記最高周波数をリアルタイム表
示する場合は、所定の演算処理を省いて最高周波数を求
め、前記最高周波数を画像フリーズ後に表示させる命令
が出力されたときは全ての演算処理を行って最高周波数
を求める。さらに、前記周波数スペクトラムをフリーズ
して表示させる命令が出力されたときのみ前記最高周波
数を求める。これらにより、前記平均処理手段、前記最
高周波数抽出手段、前記閾値演算手段に掛かる負担を軽
減させる。

【００２９】さらに、前記閾値演算手段は、前記平均処
理手段により平均処理された周波数スペクトラムの各周
波数について、そのデータ中の最大値もしくは平均値を
その周波数の信号値とするとともに、前記各周波数につ
いて、そのデータ中の所定の領域をノイズとし、前記ノ
イズとする領域のデータの平均値をノイズ値として前記
最高周波数の閾値を演算することを特徴としている。

【００３０】さらに、流れ方向を検出する手段は、前記
周波数スペクトラムの周波数軸において、データ値が最
大となる最大周波数または周波数スペクトラムの平均周
波数が周波数０に対して位置する側を流れ方向として検
出する。

【００３１】また本発明は、被検体内に超音波を送受波
し、これにより得られた受信信号からドプラ信号を抽出
し、このドプラ信号を処理して周波数スペクトラムを
得、該周波数スペクトラムの各周波数についての最高周
波数を検出し、前記周波数スペクトラムと前記最高周波
数との時間経過に伴う推移を表示する超音波診断装置に
おいて、一定期間前記最高周波数が検出できなかったと
き、該期間の直前の期間における最高周波数値及び又は
該期間の直後の期間における最高周波数値を用いて該期
間の最高周波数値を補間することを特徴とする。

【００３２】ここでいう最高周波数が“検出できなかっ
た”とは、最高周波数が実際に検出されず、検出出力が
０Ｈｚになる場合と、直前の期間で検出された最高周波
数値から当該期間に検出された最高数周波数値が所定以
上離れている場合を含むものとする。

【００３３】この補間方法には、直線補間や曲線補間が
ある。なお、実時間で周波数スペクトラム及び最高周波
数を検出して表示するリアルタイムオートトレースを行
う場合には、最後に検出された最高周波数値を検出でき
なかった期間の最高周波数値とするホールドを行う。

【００３４】また本発明は、被検体内に超音波を送受波
し、これにより得られた受信信号からドプラ信号を抽出
し、このドプラ信号を処理して周波数スペクトラムを
得、該周波数スペクトラムの各周波数についての最高周
波数を検出し、前記周波数スペクトラムと前記最高周波
数との時間経過に伴う推移を表示する超音波診断装置に
おいて、第１の期間において第１の最高周波数値が検出
されたとき、第１の期間に引き続く第２の期間における
最高周波数検出窓を前記第１の最高周波数値から所定の
範囲内とし、該最高周波数検出窓の範囲以外の最高周波
数値が検出されたときは、最高周波数が検出できなかっ
たとすることを特徴とする。

【００３５】この最高周波数検出窓ｆｗは、前記第１の
最高周波数値をｆ１とすると、定数ａ，ｂ（０＜ｂ＜
１）を使用して、 （ｆ１−ａ）≦ｆｗ≦（ｆ１＋ａ） …（２） または、

【数１】 ｆ１×（１−ｂ）≦ｆｗ≦ｆ１×（１＋ｂ） …（３） 式（２）、式（３）としてもよく、また等号を含まなく
ても良い。

【００３６】また本発明は、被検体内に超音波を送受波
し、これにより得られた受信信号からドプラ信号を抽出
し、このドプラ信号を処理して周波数スペクトラムを
得、該周波数スペクトラムの各周波数についての最高周
波数を検出し、前記周波数スペクトラムと前記最高周波
数との時間経過に伴う推移を表示する超音波診断装置に
おいて、前記周波数スペクトラム及び前記最高周波数が
表示された画面上に第１の時刻及び第２の時刻を指定す
ることが可能な入力手段と、第１の時刻及び第２の時刻
で挟まれた期間に検出された最高周波数値を除外し、該
期間以外に検出された最高周波数値を用いて該期間内の
最高周波数値を補間する補間手段と、を備えたことを特
徴とする。この補間方法には、直線補間や曲線補間があ
る。

【００３７】また本発明は、被検体内に超音波を送受波
し、これにより得られた受信信号からドプラ信号を抽出
し、このドプラ信号を処理して周波数スペクトラムを
得、該周波数スペクトラムの各周波数についての最高周
波数を検出し、前記周波数スペクトラムと前記最高周波
数との時間経過に伴う推移を表示する超音波診断装置に
おいて、前記周波数スペクトラム及び前記最高周波数が
表示された画面上に第１の時刻及び第２の時刻を指定す
るととももに、任意の最高周波数値を入力することが可
能な入力手段と、第１の時刻及び第２の時刻で挟まれた
期間に検出された最高周波数値を除外し、前記入力手段
より入力された最高周波数値に基づいて、該期間内の最
高周波数値を修正する修正手段と、を備えたことを特徴
とする。

【００３８】この修正手段による修正には、入力手段に
よる第１の時刻及び第２の時刻で挟まれた修正期間を指
定する段階と、入力手段より修正期間中の任意の時刻の
最高周波数値を入力する段階と、入力された最高周波数
値に基づいて補間または近似により修正する段階を含ん
でもよい。

【００３９】例えば、第１の時刻及び第２の時刻をそれ
ぞれ指定し、第１の時刻及び第２の時刻で挟まれた修正
期間中の任意の時刻と最高周波数値との組からなる点を
１点、または数点入力し、期間両端の２つの最高周波数
値を持つ点とこれら入力された点の値を折れ線や２次曲
線、３次曲線またはスプライン曲線で結んでも良い。ス
プライン曲線には、入力された点列を通過する内挿スプ
ライン曲線と、必ずしも点列を通らない近似スプライン
曲線とがあり、いずれを使用してもよい。さらに、グラ
フイック入力手段により、修正期間中の最高周波数値を
示す曲線そのものを入力してもよいし、この入力された
曲線に対して適当な平均化または平滑化を施してもよ
い。

【００４０】

【発明の実施の形態】次に図面を参照して本発明の実施
の形態を詳細に説明する。図１は、本発明に係る超音波
診断装置の実施の形態を示す機能ブロック図である。

【００４１】図１に示すように、本実施例の超音波診断
装置１０は、電子走査型超音波プローブ（以下超音波プ
ローブと記す）１１と、電子走査装置アナログ部１３
と、走査変換手段（ディジタル・スキャン・コンバー
タ，以下ＤＳＣと記す）３０と、Ｄ／Ａ変換器３１と、
表示器３３と、入力手段としてのトラックボール３５
と、ミキサ４０と、９０゜位相器４１と、ローパスフィ
ルタ４３と、スペクトラムドプラ演算部５０と、最高周
波数検出部８０とを備えている。

【００４２】また、電子走査装置アナログ部１３は、基
準信号発生器２０と、ディレーライン２１と、パルサ２
３と、プリアンプ２５と、加算器２７と、検波器２９と
から成る。

【００４３】まず、電子走査装置アナログ部１３では、
基準信号発生器２０から出力される基準信号を基に、デ
ィレーライン２１により超音波プローブ１１の振動子の
駆動タイミングが制御された後、パルサ２３により超音
波プローブ１１の振動子が駆動される。これにより超音
波プローブ１１の振動子は生体内に超音波を送波する。
また、その反射波信号がプリアンプ２３により増幅さ
れ、ディレーライン２１により受信タイミングが制御さ
れて加算器２５により加算される。

【００４４】次いで加算器２５から出力された信号の一
方は検出器２９に供給され、そこで対数増幅処理、包絡
線検波処理を経てＤＳＣ３０へ送られＢモード像として
形成され、Ｄ／Ａ変換器３１を介して表示器３３に表示
される。

【００４５】他方、加算器２５からの出力はミキサ４０
にも加えられる。また、９０゜位相器４１よって基準信
号発生器２０からの基準信号ｆ0 が９０゜の位相差でミ
キサ４０に加えられる。

【００４６】この結果、ローパスフィルタ４３にはドプ
ラ偏移信号ｆｄと（２ｆ0 ＋ｆｄ）信号が入力され、ロ
ーパスフィルタ４３によって高周波成分が除去されてド
プラ偏移信号ｆｄのみが得られる。これは血流情報演算
のための位相検波出力信号となる。

【００４７】ローパスフィルタ４３より出力された位相
検波出力信号はスペクトラムドプラ演算部５０に加えら
れる。なお、血流が流れている深さの位置だけのドプラ
信号を抽出するために、前記ローパスフィルタ４３から
の位相検波出力信号はサンプルホールド回路６１に出力
する。

【００４８】レンジゲート回路６０は、遅延時間が任意
に設定できるが、ここでは、超音波が超音波プローブ１
１の振動子からレンジゲート位置（サンプリングポイン
ト位置）までの往復する時間だけレートパルスより遅延
させ、設定された長さに対応する幅を有するサンプリン
グパルスをサンプリングホールド回路６１に加える。サ
ンプルホールド回路６１は、サンプリングパルスにより
ローパスフィルタ４３の出力信号をサンプルホールドす
る。バンドパスフィルタ６３は、サンプルホールド回路
６１でのサンプリングにより生じた高周波成分並びに血
管等の固定反射信号または比較的ゆっくりした動きによ
るドプラ偏移周波数を除去し、血流によるドプラ信号の
みを抽出する。バンドパスフィルタ６３からの出力は、
Ａ／Ｄ変換器６５でディジタル信号に変換された後、周
波数分析器としてのＦＦＴ６７を介して周波数スペクト
ルに変換され、ＤＳＣ３０、Ｄ／Ａ変換器３１を介して
表示器３３に表示される。

【００４９】また、前記周波数スペクトルに変換された
信号は、最高周波数検出部８０に送られる。これを基に
最高周波数検出部８０では最高周波数が検出される。こ
の検出された最高周波数も表示器３３に表示される。

【００５０】表示器３３は、ドプラ周波数の周波数スペ
クトルとともにドプラ周波数の最高周波数の時間経過に
伴う推移を表示し、その表示は、横軸（Ｘ軸）に時間、
縦軸（Ｙ軸）にドプラ偏移周波数、輝度軸（Ｚ軸）に信
号強度または信号パワーを表示しており、この表示は通
常ソナグラム表示と呼ばれている。

【００５１】トラックボール３５は、周波数スペクトラ
ム及び最高周波数の時間経過に伴う推移が表示された表
示器３３の画面上で、ノイズなどの影響により異常な最
高周波数が表示されている期間をオペレータが指定して
補間させたり、この期間の推測される最高周波数値をオ
ペレータが入力する際に使用される入力手段である。こ
の入力手段は、トラックボールに限らず、マウス、ライ
トペン、キーボード等の他の入力手段であってもかまわ
ない。

【００５２】図２は、本実施例の超音波診断装置１０の
特徴である最高周波数検出部８０を示す機能ブロック図
である。最高周波数検出部８０は、並べ変え部８１と、
リニア軸変換部８３と、スペクトラム平均部８５と、ノ
イズ検出部８７と、信号値検出部８９と、閾値演算部９
１と、流れ方向検出部９３と、最高周波数抽出部９５
と、最高周波数補間部９６と、平均処理部９７と、最高
周波数補正部９８とから成る。

【００５３】まず、ＦＦＴ６７により周波数スペクトル
に変換された信号は最高周波数検出手段８０の並べ変え
部８１に供給される。並べ変え部８１は、ＦＦＴ６７か
ら供給された周波数スペクトルの各周波数軸について平
均周波数もしくは最大周波数が中心に位置にするように
データ列の並べ換えを行う。

【００５４】例えば、図３に示すような周波数スペクト
ルの時系列がＦＦＴ６７から供給されたとする。このと
き、図３に示すある時刻における一周波数軸について、
周波数に対するパワー値（輝度値）を示すと図４（ａ）
のようになり、ドプラ偏移周波数０が中心となってい
る。なお、図３、図４中、ｔは時間、ｆは周波数、Ｐは
パワー値を表し、以下図５〜図７についても同様とす
る。

【００５５】この状態で平均処理等を行うとデータの折
り返りにより正確な処理が行えないため、図４（ｂ）に
示すように平均周波数ｆmean（または最大周波数）が中
心にくるようにデータ列の並べ変えを行う。なお、デー
タの折り返りが生じない範囲の周波数帯である場合、こ
の処理は省略しても良い。

【００５６】リニア軸変換部８３は、ＦＦＴ６７により
ログ変換されたデータをリニア軸の数に変換する。な
お、前記ログ変換されたデータを用いて演算処理を行っ
ても誤差が少ない場合、この処理は省略しても良いし、
上端、下端どちらか一方でも良い。

【００５７】スペクトラム平均部８５は、前記並べ変え
部８１によりデータ列の並べ変え、前記リニア軸変換部
８３によりリニア軸の数に変換された周波数スペクトラ
ムについて、周波数軸方向、時間軸方向に平均処理を行
う。なお、前記周波数軸方向、時間軸方向の内、いずれ
か一方向について平均処理するようにしても良い。

【００５８】ノイズ検出部８７は、スペクトラム平均部
８５により平均処理された周波数スペクトラムの各周波
数軸について、上端側、下端側それぞれ所定の領域をノ
イズ領域Ｎｒとし、このノイズ領域Ｎｒ内のパワー値の
平均をノイズ値Ｎとする。例えば、図５（ａ）に示す一
周波数軸についてのノイズ領域Ｎｒは、図５（ｂ）に示
すような領域となる。なお、ここでは各周波数軸につい
て、上端側、下端側それぞれ所定の領域をノイズ領域Ｎ
ｒとしてノイズ値を求めているが、これに限らず、複数
の周波数軸について、上端側、下端側それぞれ所定の領
域をノイズ領域Ｎｒとしてノイズ値を求めるようにして
も良い。

【００５９】また、最高周波数の検出（表示）を、ある
一定時間幅を指示して行う場合には、その時間幅内の所
定領域をＮｒとしても良い。

【００６０】信号値検出部８９は、スペクトラム平均部
８５により平均処理された周波数スペクトラムの各周波
数軸について、パワー値の最大値を信号値Ｓとする。例
えば図６に示すような周波数スペクトルがスペクトラム
平均部８５により供給されたとき、図６に示す一周波数
軸については図６（ｂ）に示すように、パワー値の最大
値（Ｐｅａｋ）が信号値Ｓとなる。

【００６１】なお、前記各周波数軸について、最大値の
前後所定数の周波数に対応するパワー値の平均もしくは
平均周波数に対応するパワー値を信号値Ｓとするように
しても良い。

【００６２】閾値演算部９１は、ノイズ検出部８７によ
りノイズ値、信号値検出部８９により信号値が得られる
と、前記得られたノイズ値と信号値を基に以下に示す式
（４）ないし式（６）のいずれかを用いて閾値Ｔｈを演
算する。 Ｔｈ＝（Ｓ＋Ｎ）／２ …（４） またはＴｈ＝（Ｓ／Ｎ）^1/2 ＊Ｎ …（５） またはＴｈ＝Ｎ＊Ｃ（ｃｏｎｓｔ．） …（６）

【００６３】流れ方向検出部９３は、スペクトラム平均
部８５により平均処理された周波数スペクトラムの各周
波数軸について、最大周波数ｆpeak（または平均周波数
ｆmean）がプラス側にある場合は、流れの方向を順流方
向（超音波プローブに近付く方向）とし、最大周波数ｆ
peak（または平均周波数ｆmean）がマイナス側にある場
合は、流れの方向を逆流方向（超音波プローブから遠ざ
かる方向）とする。また、スペクトラム平均部８５によ
り平均処理された周波数スペクトラムを最高周波数抽出
部９５に供給する。

【００６４】例えば図６に示すような周波数スペクトル
がスペクトラム平均部８５により供給され、このとき、
図６に示す一周波数軸について、周波数に対するパワー
値を示すと図６（ｂ）に示すようになる。この場合、最
大周波数はプラス側にあるので流れの方向は順流方向と
検出される。なお、前記流れの方向および閾値は、操作
者が入力するようにしても良い。

【００６５】最高周波数抽出部９５は、閾値演算部９１
により閾値が演算されるとともに、流れ方向検出部９３
により流れ方向が検出され、かつ、流れ方向検出部９３
から周波数スペクトラムが供給されると、前記周波数ス
ペクトラムの各周波数軸について、前記検出された流れ
の方向の端部から０Hz側に向かってパワー値をサーチし
て行き、前記パワー値が閾値を連続してＮ個越えた最初
の位置のドプラ偏移周波数を最高周波数ｆmax とする。

【００６６】例えば、閾値を連続して４個越えた最初の
位置のドプラ偏移周波数を最高周波数ｆmax とするよう
に最高周波数抽出部９５に設定し、そして、図７（ａ）
に示すような周波数スペクトラムがスペクトラム平均部
８５から流れ方向検出部９３を介して供給され、図７
（ａ）に示す一周波数について、周波数に対するパワー
値が図７（ｂ）に示すように得られたとする。このと
き、流れの方向は順流方向であるので、最高周波数抽出
部９５はプラス側の端部から０Hz側に向かってサーチし
ていき、図７（ｂ）に示すように閾値を連続して４個越
えた最初の位置のドプラ偏移周波数を最高周波数ｆmax
とする。

【００６７】なお、ここでは前記パワー値が閾値を連続
してＮ個越えた最初の位置のドプラ偏移周波数を最高周
波数ｆmax としているが、連続して越えていなくても、
Ｍ個中Ｒ個越えた場合は、その最初に越えた位置のドプ
ラ偏移周波数を最高周波数ｆmax とするようにしても良
い。さらに、スペクトラム平均部８５により時間軸方向
に平均処理を行った場合は、その平均幅の分、元に戻す
ように最高周波数を補正しても良い。

【００６８】最高周波数抽出部９５による周波数スペク
トラムのサーチにおいて、上記のように流れ方向端部か
らサーチする他に、パワー値が最大となる最大周波数ま
たはパワー値が平均となる平均周波数から流れ方向端部
に向かってサーチしても良い。この場合には、閾値より
十分大きいパワー値を持つ周波数から検索を開始して、
各周波数のパワー値が閾値を下回る周波数を求めること
となるので、例えば、閾値を連続して４個下回った最初
の位置のドプラ偏移周波数を最高周波数ｆmaxとする。

【００６９】また最高周波数抽出部９５による最高周波
数抽出において、最高周波数検出窓という概念を導入
し、最高周波数抽出の耐ノイズ性を高めることも考えら
れる。即ち、ある時刻ｔ１の周波数スペクトルから最高
周波数ｆ１が抽出されたとしたとき、次の時刻ｔ２の周
波数スペクトルから抽出される最高周波数はｆ１から大
きく離れているとは考えられないので、ｆ１に基づいて
妥当な範囲であるｆｗを定数ａ，ｂを使用して、 （ｆ１−ａ）≦ｆｗ≦（ｆ１＋ａ） …（２） または、

【数２】 ｆ１×（１−ｂ）≦ｆｗ≦ｆ１×（１＋ｂ） …（３） 式（２）または式（３）によって定める。

【００７０】そして、時刻ｔ２の周波数スペクトルから
抽出される最高周波数は、このｆｗの範囲に入る最高周
波数値のみを有効な値とし、この範囲に入らない値はノ
イズ等による誤りとして最高周波数値が検出されなかっ
たとする。

【００７１】また、従来は図８（ａ）に示すように、最
高周波数値が実際に検出されず０Ｈｚとなった場合、そ
の期間の最高周波数値の軌跡は、０Ｈｚとなり不自然な
表示となっていたが、本発明を適用すれば、最高周波数
値が実際に検出されず０Ｈｚとなった場合、または上記
最高周波数検出窓の範囲に入らなかった場合に、その期
間が一定範囲内であれば、この期間の最高周波数値を補
間することができる。この補間機能は、最高周波数補間
部９６のソフトウェアにより実現され、最高周波数値が
検出されなかった期間の前後の期間に抽出された最高周
波数値を用いて、最高周波数値を補間するものであり、
その表示は図８（ｂ）に示すようになる。補間方法に
は、直線補間、２次曲線補間、３次曲線補間等があり、
いずれの補間方法を選んでも良い。

【００７２】なお、実時間で周波数スペクトラム及び最
高周波数を検出して表示するリアルタイムオートトレー
スを行う場合には、最後に検出された最高周波数値をホ
ールドし、この値を最高周波数が検出されなかった期間
の最高周波数値として表示し、これ以後新たに最高周波
数値が検出されるとその値を表示するようにしてもよ
い。

【００７３】平均処理部９７は、最高周波数抽出部９５
により最高周波数が求められると、この最高周波数に対
し、時間軸方向に平均処理を施す。なお、平均処理部９
７による平均処理は省略しても良い。こうして得られた
最高周波数はＤＳＣ３０に供給され、Ｄ／Ａ変換器３１
を介して表示器３３に表示される。

【００７４】このように本実施例の超音波診断装置１０
では、周波数分析器であるＦＦＴ６７から供給される周
波数スペクトラムを各周波数軸について、データ列の並
べ変え、リニア軸の数への変換および平均処理を行った
後、ノイズ値と信号値を求め、それを基に閾値を演算す
るとともに流れ方向を検出し、前記検出された流れの方
向の端部から０Hz側に向かって、または最大のパワー値
を持つ最大周波数ないし平均のパワー値を持つ平均周波
数から流れ方向端部へ向かって、パワー値をサーチして
行き、前記パワー値が閾値を連続してＮ個越えた最初の
位置のドプラ偏移周波数を最高周波数ｆmax としている
ので、ノイズの影響を少なくすることができ、そのた
め、最高周波数の検出精度を向上させることが可能とな
る。また、特別な操作を行うこと無く前記検出精度を向
上させることができるので、操作者の負担が軽減し、さ
らに、診断時間を著しく短縮することができる。

【００７５】なお、最高周波数をリアルタイム表示する
場合は、演算処理時間を短縮させるために所定の演算処
理（並べ変え部８１によるデータ列の並べ変え、リニア
軸変換部８３によるリニア軸の数への変換、スペクトラ
ム平均部８５による周波数軸方向の平均処理または時間
軸方向の平均処理、平均処理部９７による平均処理の
内、少なくとも一つ）を省いて最高周波数を求め、ま
た、最高周波数をフリーズ後に表示させる命令が出力さ
れたときは全ての演算処理を行って最高周波数を求める
ようにしても良い。さらに、本実施例の超音波診断装置
１０では、最高周波数をリアルタイムで求める場合を例
に説明したが、これに限らず、例えば、周波数スペクト
ルをフリーズさせる命令が図示しない入力手段を介して
オペレータにより出されたときのみ前記最高周波数を求
めるようにしても良い。これらにより、スペクトラム平
均部８５、閾値演算部９１、最高周波数抽出部９５に掛
かる負担を軽減させることができる。

【００７６】また、並べ変え部８１と、リニア軸変換部
８３とスペクトラム平均部８５とノイズ検出部８７と信
号値検出部８９と流れ方向検出部９３と閾値演算部９１
と最高周波数抽出部９５と平均処理部９７の各処理につ
いては、その一部または全てをソフトウェアで行うよう
にしても良い。

【００７７】また、本発明によれば、オペレータの判断
により最高周波数の時間経過をマニュアルで修正するこ
とができる。このマニュアル修正機能は、最高周波数補
間部９６及び最高周波数修正部９８のソフトウェアによ
って実現され、ノイズなどの影響を受けた周波数スペク
トラムにより上記の自動最高周波数検出を行った結果、
ほんの僅かの最高周波数のミストレースが生じた場合、
この僅かのミストレースのために再度データを収集、分
析するのではなく、ノイズによりミストレースした部分
をマニュアル修正するものである。

【００７８】図９（ａ）ないし（ｄ）は、マニュアル修
正機能による最高周波数トレースの修正の様子を示す周
波数スペクトラム及び最高周波数を表示する画面の変遷
である。また図１０は、マニュアル修正の操作を示すフ
ローチャートである。

【００７９】図１０のフローチャートにおいて、まず、
ドプラ周波数偏移の周波数スペクトラム及び自動検出さ
れた最高周波数の時間経過であるトレースが表示器３３
に表示される（ステップＳ１０）。この表示が例えば、
図９（ａ）に示すようにノイズ等の影響により、ドプラ
周波数偏移トレースの一部に最高周波数が高い値へ突出
していたとする。

【００８０】オペレータはノイズの影響の程度を判断し
て、再度データを取るか、マニュアル修正機能を使用し
てデータを修正するかを決め、マニュアル修正をする場
合、修正作業開始を入力する（ステップＳ１２）。修正
作業開始の入力には、修正作業開始スイッチを別途設け
て認識してもよいし、トラックボールにより修正作業開
始のアイコンを指定するか又はトラックボール３５の作
動を検知して入力作業開始を認識してもよい。

【００８１】次いで、最高周波数トレースを修正すべき
修正期間の始点ｔ１にトラックボールによりカーソルを
合わせて始点ｔ１を指定する（ステップＳ１４）。次い
で、修正期間の終点ｔ２にトラックボールによりカーソ
ルを合わせて終点ｔ２を指定する（ステップＳ１６）。
始点と終点とが指定された状態を図９（ｂ）に示す。な
お、始点と終点の指定は順序が入れ替わってもよい。こ
こで、修正期間が指定されると同時に修正期間の最高周
波数の自動トレースを消してもよいが、補間によるトレ
ースまたはマニュアルトレースにより最高周波数トレー
スが確定するまで異なる色で２つのトレースを表示して
もよい。

【００８２】次いで、修正期間の最高周波数値をマニュ
アルで入力するか、入力しないで修正期間の前後の期間
の値から補間するかを指定する（ステップＳ１８）。こ
れは極めて僅かのミストレースであれば、人手により修
正値を入力しなくても補間により十分修正できるからで
ある。

【００８３】マニュアルトレースしない場合には（ステ
ップＳ１８でＮｏの場合）、修正期間の前後の期間の最
高周波数値から補間を行い（ステップＳ３０）、修正期
間の最高周波数値を補間曲線（または補間直線）により
代替して（ステップＳ３２）、終了する。

【００８４】一方マニュアルトレースする場合には（ス
テップＳ１８でＹｅｓの場合）、トラックボールにより
カーソルを移動させて、時間軸座標と周波数軸座標とか
らなる修正点を入力する（ステップＳ２０）。修正点を
入力した状態を図９（ｃ）に示す。修正点は複数入力し
てもよい。修正点の入力が終了すれば（ステップＳ２２
でＹｅｓ）、修正期間両端の２つの最高周波数値と、入
力された修正点（または修正点列）よりスプライン曲線
を計算して表示器の画面に描画する（ステップＳ２
４）。次いで、修正期間の最高周波数値をスプライン曲
線により代替して（ステップＳ２６）、マニュアル修正
作業を終了する。

【００８５】なお、図１０のフローチャートでは、１つ
または複数の修正点を入力してスプライン曲線により最
高周波数トレースを修正するとしたが、トレース（修正
曲線）そのものを入力してもよい。この場合には、入力
されたトレースを自動的に滑らかにする平均化ないし平
滑化をさらに行ってもよい。

【００８６】以上好ましい実施の形態をパルスドプラを
用いた場合を例にして説明したが、これに限らず、連続
波ドプラを用いた場合でも、本発明を適用できることは
明らかである。

【００８７】

【発明の効果】以上説明したように本発明によれば、周
波数スペクトルの各周波数について、そのデータ値の
内、流れ方向側端部から順にサーチしていき、該値が予
め設定されている閾値を所定のデータ幅中で所定数越え
た時、その最初に越えた位置を最高周波数としているの
で、ノイズの影響を少なくすることができ、そのため、
最高周波数の検出精度を向上させることが可能となる。
また、特別な操作を行うこと無く前記検出精度を向上さ
せることができるので、操作者の負担が軽減し、さら
に、診断時間を著しく短縮することができる。

【００８８】また、周波数スペクトラムを時間方向、周
波数方向の少なくとも一方向について平均処理した後、
その周波数スペクトルの各周波数について、そのデータ
値の内、流れ方向側端部から順にサーチしていき、該値
が予め設定されている閾値を所定のデータ幅中で所定数
越えた時、その最初に越えた位置を最高周波数としてい
るので、ノイズの影響を少なくすることができ、そのた
め、最高周波数の検出精度を向上させることが可能とな
る。また、特別な操作を行うこと無く前記検出精度を向
上させることができるので、操作者の負担が軽減し、さ
らに、診断時間を著しく短縮することができる。

【００８９】さらに、周波数スペクトラムを時間方向、
周波数方向の少なくとも一方向について平均処理した
後、その周波数スペクトルの各周波数について、その信
号値とノイズ値とを基に、最高周波数の閾値を演算する
とともに、前記平均処理された周波数スペクトルの各周
波数について、そのデータ値の内、流れ方向側端部から
順にサーチしていき、該値が前記閾値演算手段により演
算された閾値を所定のデータ幅中で所定数越えた時、そ
の最初に越えた位置を最高周波数としているので、ノイ
ズの影響を少なくすることができ、そのため、最高周波
数の検出精度を向上させることが可能となる。また、特
別な操作を行うこと無く前記検出精度を向上させること
ができるので、操作者の負担が軽減し、さらに、診断時
間を著しく短縮することができる。

【００９０】さらに、最高周波数をリアルタイム表示す
る場合は、所定の演算処理を省いて最高周波数を求め、
前記最高周波数をフリーズ後に表示させる命令が出力さ
れたときは全ての演算処理を行って最高周波数を求める
ようにしているので、前記演算処理に掛かる負担を軽減
させることができる。

【００９１】さらに、周波数スペクトラムをフリーズさ
せる命令が出力されたときのみ前記最高周波数を求める
ようにしているので、前記演算処理に掛かる負担を軽減
させることができる。

【００９２】さらに、ノイズ等の影響により、最高周波
数トレースが少しだけ誤った場合には、マニュアルで誤
った期間を修正期間と指定し、修正期間の前後の期間の
最高周波数値から補間したり、必要ならば推測される正
しい値である修正点ないし修正トレースそのものをマニ
ュアルで入力することにより、誤ったデータを簡単に修
正することができるので、再度データ収集する必要がな
くなり、診断時間を短縮し、操作者の負担を軽減するこ
とができる。

【図面の簡単な説明】

【図１】本発明に係る超音波診断装置の一実施例を示す
機能ブロック図である。

【図２】図１に示す超音波診断装置の最高周波数検出部
８０を示す機能ブロック図である。

【図３】周波数スペクトルの一例を示す説明図である。

【図４】図３に示す周波数スペクトルの一周波数軸につ
いて、周波数に対するパワー値を示す説明図である。

【図５】周波数スペクトル上のノイズ領域Ｎｒ（図５
（ａ））と、周波数スペクトルの一周波数軸についての
ノイズ領域Ｎｒ（図５（ｂａ））を示す説明図である。

【図６】周波数スペクトル（図６（ａ））の一周波数軸
について、パワー値の最大値（図６（ｂ））を示す説明
図である。

【図７】周波数スペクトル（図７（ａ））の一周波数軸
について、最高周波数ｆmax （図７（ｂ））を示す説明
図である。

【図８】従来の最高周波数値が検出できなかった期間を
含むトレース（図８（ａ））、及び本発明による最高周
波数値の補間を含むトレース（図８（ｂ））を示す説明
図である。

【図９】マニュアル修正機能による最高周波数トレース
の修正の様子を示す周波数スペクトラム及び最高周波数
値を表示する画面の変遷を示す説明図である。

【図１０】マニュアル修正の操作を示すフローチャート
である。

【図１１】従来の超音波診断装置を示す機能ブロック図
である。

【図１２】従来の超音波診断装置により発生されるクロ
ックパルスと、レートパルスと、サンプリングパルスお
よび超音波診断装置により受信されたエコー信号を示す
説明図である。

【図１３】従来の最高周波数の検出アルゴリズムの例を
示す説明図である。

【符号の説明】

１０ 超音波診断装置、 １１ 超音波プローブ、 １
３ 電子走査装置アナログ部、 ２０ 基準信号発生
器、 ２１ ディレーライン、 ２３ パルサ、２５
プリアンプ、 ２７ 加算器、 ２９ 検波器、 ３０
ＤＳＣ、 ３１ Ｄ／Ａ変換器、 ３３ 表示器、
３５ トラックボール、 ４０ ミキサ、 ４１ ９０
゜位相器、 ４３ ローパスフィルタ、 ５０ スペク
トラムドプラ演算部、 ６０ レンジゲート回路、 ６
１ サンプルホールド回路、 ６３ バンドパスフィル
タ、 ６５ Ａ／Ｄ変換器、 ６７ ＦＦＴ、 ８０
最高周波数検出部、 ８１ 並べ変え部、 ８３ リニ
ア軸変換部、 ８５ スペクトラム平均部、 ８７ ノ
イズ検出部、 ８９ 信号値検出部、 ９１ 閾値演算
部、 ９３ 流れ方向検出部、 ９５ 最高周波数抽出
部、 ９６ 最高周波数補間部、 ９７ 平均処理部、
９８ 最高周波数修正部。

Claims (12)

【特許請求の範囲】
1. 【請求項１】 被検体内に超音波を送受波し、これによ
   り得られた受信信号からドプラ信号を抽出し、このドプ
   ラ信号を処理して周波数スペクトラムを得、この周波数
   スペクトラムを表示するとともに、前記周波数スペクト
   ラムの各周波数についての最高周波数を検出する超音波
   診断装置において、 前記周波数スペクトラムの各周波数の信号強度を示すデ
   ータ値について、 流れ方向側端部から周波数０Ｈｚの方向に向かって順に
   検索し、該データ値が予め設定された閾値を所定の周波
   数幅中で所定回数以上越えた時、その最初に越えた位置
   を最高周波数とするか、 或いは最も信号強度の大きい周波数ないし平均周波数か
   ら流れ方向側端部の方向に向かって順に検索し、該デー
   タ値が予め設定された閾値を所定の周波数幅中で所定回
   数以上越えなかった時、その最初に越えなかった位置を
   最高周波数とすることを特徴とする超音波診断装置。
2. 【請求項２】 被検体内に超音波を送受波し、これによ
   り得られた受信信号からドプラ信号を抽出し、このドプ
   ラ信号を処理して周波数スペクトラムを得、この周波数
   スペクトラムを表示するとともに、前記周波数スペクト
   ラムの各周波数についての最高周波数を検出する超音波
   診断装置において、 前記周波数スペクトラムを時間方向、周波数方向の少な
   くとも一方向について平均処理する平均処理手段と、 前記平均処理手段により平均処理された周波数スペクト
   ラムの各周波数の信号強度を示すデータ値について、流
   れ方向側端部から周波数０Ｈｚの方向に向かって順に検
   索し、該データ値が予め設定された閾値を所定の周波数
   幅中で所定回数以上越えた時、その最初に越えた位置を
   最高周波数とするか、或いは最も信号強度の大きい周波
   数ないし平均周波数から流れ方向側端部の方向に向かっ
   て順に検索し、該データ値が予め設定された閾値を所定
   の周波数幅中で所定回数以上越えなかった時、その最初
   に越えなかった位置を最高周波数とする最高周波数抽出
   手段と、 を備えることを特徴とする超音波診断装置。
3. 【請求項３】 被検体内に超音波を送受波し、これによ
   り得られた受信信号からドプラ信号を抽出し、このドプ
   ラ信号を処理して周波数スペクトラムを得、この周波数
   スペクトラムを表示するとともに、前記周波数スペクト
   ラムの各周波数についての最高周波数を検出する超音波
   診断装置において、 前記周波数スペクトラムを時間方向、周波数方向の少な
   くとも一方向について平均処理する平均処理手段と、 前記平均処理手段により平均処理された周波数スペクト
   ラムの各周波数について、その信号値とノイズ値とを基
   に、最高周波数の閾値を演算する閾値演算手段と、 前記平均処理手段により平均処理された周波数スペクト
   ラムの各周波数の信号強度を示すデータ値について、流
   れ方向側端部から周波数０Ｈｚの方向に向かって順に検
   索し、該データ値が予め設定された閾値を所定の周波数
   幅中で所定回数以上越えた時、その最初に越えた位置を
   最高周波数とするか、或いは最も信号強度の大きい周波
   数ないし平均周波数から流れ方向側端部の方向に向かっ
   て順に検索し、該データ値が予め設定された閾値を所定
   の周波数幅中で所定回数以上越えなかった時、その最初
   に越えなかった位置を最高周波数とする最高周波数抽出
   手段と、 を備えることを特徴とする超音波診断装置。
4. 【請求項４】 前記閾値演算手段は、前記平均処理手段
   により平均処理された周波数スペクトラムの各周波数軸
   でその信号値とノイズ値との比を基に、最高周波数の閾
   値を演算することを特徴とする請求項３記載の超音波診
   断装置。
5. 【請求項５】 前記閾値演算手段は、前記平均処理手段
   により平均処理された周波数スペクトラムの各周波数に
   ついて、そのデータ中の最大値もしくは平均値をその周
   波数の信号値とするとともに、前記各周波数について、
   そのデータ中の所定の領域をノイズとし、前記ノイズと
   する領域のデータの平均値をノイズ値として前記最高周
   波数の閾値を演算することを特徴とする請求項３または
   請求項４いずれか記載の超音波診断装置。
6. 【請求項６】 被検体内に超音波を送受波し、これによ
   り得られた受信信号からドプラ信号を抽出し、このドプ
   ラ信号を処理して周波数スペクトラムを得、この周波数
   スペクトラムを表示するとともに、前記周波数スペクト
   ラムの各周波数についての最高周波数を検出する超音波
   診断装置において、 前記最高周波数をリアルタイム表示する場合は、所定の
   演算処理を省いて最高周波数を求め、前記最高周波数を
   画像フリーズ後に表示させる命令が出力されたときは全
   ての演算処理を行って最高周波数を求めることを特徴と
   する請求項１乃至請求項５いずれか１項記載の超音波診
   断装置。
7. 【請求項７】 被検体内に超音波を送受波し、これによ
   り得られた受信信号からドプラ信号を抽出し、このドプ
   ラ信号を処理して周波数スペクトラムを得、この周波数
   スペクトラムを表示するとともに、前記周波数スペクト
   ラムの各周波数についての最高周波数を検出する超音波
   診断装置において、 前記周波数スペクトラムをフリーズさせる命令が出力さ
   れたときのみ前記最高周波数を求めることを特徴とする
   請求項１乃至請求項５いずれか１項記載の超音波診断装
   置。
8. 【請求項８】 前記周波数スペクトラムの周波数軸にお
   いて、データ値が最大となる最大周波数または周波数ス
   ペクトラムの平均周波数が周波数０に対して位置する側
   を流れ方向として検出する手段を備えたことを特徴とす
   る請求項１乃至請求項７いずれか１項記載の超音波診断
   装置。
9. 【請求項９】 被検体内に超音波を送受波し、これによ
   り得られた受信信号からドプラ信号を抽出し、このドプ
   ラ信号を処理して周波数スペクトラムを得、該周波数ス
   ペクトラムの各周波数についての最高周波数を検出し、
   前記周波数スペクトラムと前記最高周波数との時間経過
   に伴う推移を表示する超音波診断装置において、 一定期間前記最高周波数が検出できなかったとき、該期
   間の直前の期間における最高周波数値及び又は該期間の
   直後の期間における最高周波数値を用いて該期間の最高
   周波数値を補間することを特徴とする超音波診断装置。
10. 【請求項１０】 被検体内に超音波を送受波し、これに
    より得られた受信信号からドプラ信号を抽出し、このド
    プラ信号を処理して周波数スペクトラムを得、該周波数
    スペクトラムの各周波数についての最高周波数を検出
    し、前記周波数スペクトラムと前記最高周波数との時間
    経過に伴う推移を表示する超音波診断装置において、 第１の期間において第１の最高周波数値が検出されたと
    き、第１の期間に引き続く第２の期間における最高周波
    数検出窓を前記第１の最高周波数値から所定の範囲内と
    し、該最高周波数検出窓の範囲以外の最高周波数値が検
    出されたときは、最高周波数が検出できなかったとする
    ことを特徴とする超音波診断装置。
11. 【請求項１１】 被検体内に超音波を送受波し、これに
    より得られた受信信号からドプラ信号を抽出し、このド
    プラ信号を処理して周波数スペクトラムを得、該周波数
    スペクトラムの各周波数についての最高周波数を検出
    し、前記周波数スペクトラムと前記最高周波数との時間
    経過に伴う推移を表示する超音波診断装置において、 前記周波数スペクトラム及び前記最高周波数が表示され
    た画面上に第１の時刻及び第２の時刻を指定することが
    可能な入力手段と、 第１の時刻及び第２の時刻で挟まれた期間に検出された
    最高周波数値を除外し、該期間以外に検出された最高周
    波数値を用いて該期間内の最高周波数値を補間する補間
    手段と、 を備えたことを特徴とする超音波診断装置。
12. 【請求項１２】 被検体内に超音波を送受波し、これに
    より得られた受信信号からドプラ信号を抽出し、このド
    プラ信号を処理して周波数スペクトラムを得、該周波数
    スペクトラムの各周波数についての最高周波数を検出
    し、前記周波数スペクトラムと前記最高周波数との時間
    経過に伴う推移を表示する超音波診断装置において、 前記周波数スペクトラム及び前記最高周波数が表示され
    た画面上に第１の時刻及び第２の時刻を指定するととも
    もに、任意の最高周波数値を入力することが可能な入力
    手段と、 第１の時刻及び第２の時刻で挟まれた期間に検出された
    最高周波数値を除外し、前記入力手段より入力された最
    高周波数値に基づいて該期間内の最高周波数値を修正す
    る修正手段と、 を備えたことを特徴とする超音波診断装置。