JP2572913B2

JP2572913B2 - 超音波診断装置

Info

Publication number: JP2572913B2
Application number: JP3277805A
Authority: JP
Inventors: 雅彦権藤
Original assignee: Olympus Optical Co Ltd
Current assignee: Olympus Corp
Priority date: 1991-10-24
Filing date: 1991-10-24
Publication date: 1997-01-16
Anticipated expiration: 2012-01-16
Also published as: JPH05115477A

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】この発明は、生体内の血流等の移
動物体の速度情報を映像化する超音波診断装置に関する
ものである。

【０００２】

【従来の技術】ドップラ効果を利用して血流速度の分布
を測定し、超音波反射の映像と重畳して表示する超音波
診断装置は、カラードップラ装置と呼ばれ、例えば特公
平３−23050 号公報、特開平１−99541 号公報、特開平
２−268747号公報に開示されていると共に、広く用いら
れている。この血流速度分布の測定原理は、振動子アレ
イを用いた超音波ビームによりパルス波を一定周期で送
信し、反射体より反射波がかえってくるまでの時間を測
定するとともに、受信信号の周波数変化を検出すること
で、反射体の位置と動きを測定するものである。このカ
ラードップラ装置は、はじめのうちは心臓等の循環器系
を対象としていたが、最近ではドップラのもつ情報量の
多さから腹部臓器の診断にも多用されるようになってき
た。

【０００３】図７は、従来のカラードップラ装置のドッ
プラ測定関連部の構成を示すものである。ここで、１は
振動子本体、２は振動素子、102 は遅延回路、103 は切
り換え回路、５は信号発生器、105 はパルス発生回路、
106 は加算回路、107 は受信増幅回路、10はバンドパス
フィルタ(BPF) 、11は90度移相器、12a,12bは乗算器、1
3a, 13bはローパスフィルタ(LPF) 、14a, 14bはＡ／Ｄ
コンバータ，113 は血流成分のみを取り出す MTIフィル
タ、114 は周波数成分を検出する自己相関回路をそれぞ
れ示す。

【０００４】このカラードップラ装置においては、信号
発生器５の出力をパルス発生器105に供給して中心周波
数ｆo のパルスを発生させ、このパルスを切り換え回路
103および遅延回路102 を経て振動子本体１に供給し
て、超音波ビーム100 が所定の方向に向くように個々の
振動素子２を所定のタイミングで駆動している。また、
超音波ビーム100 を横切る血管101 内の血球からの反射
波は、振動子本体１の各振動素子２で電気信号に変換し
て遅延回路102 および切り換え回路103 を経て加算回路
106 で加算し、その出力を受信増幅回路107 で適正な大
きさに増幅した後、BPF10 で帯域外の雑音をカットし
て、乗算器12a, 12bおよびLPF13a, 13b で構成される直
交検波回路に供給し、これにより信号発生器５と90度移
相器11の出力である中心周波数ｆo の基準信号からの差
分周波数成分Δｆを検出している。

【０００５】ここで、差分周波数成分Δｆは、血流のド
ップラ周波数で、このドップラ周波数は超音波ビーム10
0 と血管101 のなす角度に依存し、実際の血流速度Ｖの
ビーム方向速度成分Ｖｄに対応する。このため、血流速
度は振動子本体１と血管101の間の角度に大きく依存す
ることになる。

【０００６】直交検波回路を構成する LPF13a, 13bの出
力は、Ａ／Ｄコンバータ14a, 14bでそれぞれデジタル信
号に変換して MTIフィルタ113 に供給し、ここで直交検
波出力に含まれる静止組織からの信号に対応する直流成
分の信号 (クラッタ信号) を取り除いて血流反射信号の
みを抽出し、その出力を自己相関回路114 に供給して差
分周波数成分Δｆを測定し、これを図示しないデジタル
スキャンコンバータや表示器を用いて、組織の反射信号
の強さを表わすＢモード像に重畳させて、赤や青で色付
けして表示している。なお、血球からの反射波のドップ
ラ成分は、ばらつきが多いため、ある方向に対して10回
程度送信を繰り返し、その平均を求めることによりドッ
プラの測定精度を向上させるようにしている。

【０００７】

【発明が解決しようとする課題】しかしながら、上述し
た従来のカラードップラ装置にあっては、低速血流の観
察時において静止組織からのクラッタ信号と血流信号と
の分離が困難であるために、測定不能になるという問題
がある。すなわち、このクラッタ信号と血流信号との分
離には、MTI フィルタ113 が用いられているが、このMT
I フィルタ113 はデジタル回路によって構成され、その
通過特性は、図８に示すように、応答特性とのかねあい
もあって、あまり急峻にできない。このため、クラッタ
信号を抑圧しようとすると血流信号まで減衰させる結果
となって、血流速度の小さい信号の観測ができなくな
る。また、ドップラ効果は超音波ビーム方向について生
じるため、距離方向の血流速度は検出できるが、ビーム
と直交する方向の血流については、その血流速度を測定
することができないという本質的な問題がある。

【０００８】一方、超音波ビームは、各振動素子２から
の位相合成で作られるが、各振動素子２と反射体との関
係をミクロ的にみると、遠距離部の反射体を見込む角度
は、各振動素子２間で大幅な変化はないが、振動素子２
に近接したところでの反射体を見込む角度は、各振動素
子２間で大きく異なる。また、ドップラ効果は、振動素
子２から見た反射体までの見込み角度に大きく左右され
る。このため、近接した反射体の場合には、ドップラ周
波数に大きな偏差が生じ、測定精度が悪くなるという問
題がある。この問題は、ダイナミックフォーカス等を用
いて、近距離部の空間分解能を改善しても、解決するこ
とができない。

【０００９】さらに、上述した従来のカラードップラ装
置においては、Ｂモード像を得るための通常の送信シー
ケンス以外に、微弱な血流情報を検出するためのドップ
ラ専用の時間幅の大きいパルスを送信するシーケンスが
ある。また、良好なドップラ信号を得るには、10回程度
のシーケンス間の平均操作が必要となる。このため、従
来のカラードップラ装置では、Ｂモード像のみの場合と
比べて、大幅にフレームレイトが低下したり、空間分解
能が低下するという問題がある。

【００１０】この発明は、上述した従来の問題点に鑑み
てなされたもので、開口合成技術に空間周波数の概念を
取り入れることにより、Ｂモード像と同様に高い空間分
解能で、フレームレイトを低下させることなく、低速血
流も容易に測定できると共に、振動子面と平行な方向の
血流も測定できるよう適切に構成した超音波診断装置を
提供することを目的とする。

【００１１】

【課題を解決するための手段】上記目的を達成するた
め、この発明では、アレイ状に並べた振動素子群と、そ
れらを切り替えるためのマルチプレクサと、このマルチ
プレクサの出力を増幅する受信増幅回路と、その出力を
複素数信号としてデジタル信号に変換するＡ／Ｄコンバ
ータと、該デジタル信号を各振動素子に対応させて、複
素数の波面データとして格納する波面メモリと、これら
波面メモリに格納した複素数の波面データに基づいて空
間周波数を求めるスペクトル分布検出器とを具え、この
スペクトル分布検出器で求めた空間周波数の直流成分か
ら超音波反射強度を求め、空間周波数の直流以外の周波
数成分から血流速度を求めるよう構成する。

【００１２】

【作用】すなわち、この発明では、マルチプレクサによ
り振動素子を順次切り替えながら、送信・受信を繰り返
し、各振動素子によって得られた受信信号を波面データ
として波面メモリに格納し、これらの波面データから任
意の焦点を結ぶように時間調整してサンプリングし、こ
の一連のサンプリングデータを空間変数としてとらえて
スペクトル分布検出器により空間周波数を検出し、この
スペクトル分布により血流速度を求める。

【００１３】ここで、スペクトル分布検出器で検出され
る空間周波数は、各振動素子と時間と共に変化する血球
の相対位置関係によって関連づけられ、例えば血球が近
づく場合はプラスの周波数、遠ざかる場合はマイナスの
周波数、左右に動く場合はプラスとマイナスの周波数成
分を含むことになる。したがって、このスペクトル分布
の違いを検出すれば、超音波ビーム方向の血流速度のみ
ならず、ビームと直交する方向の血流速度も測定するこ
とができる。また、この空間周波数の直流成分は、単純
に反射強度を表すので、これをＢモード像のデータと
し、その他の周波数成分は血流成分を表すドップラデー
タとして表示することができる。

【００１４】

【実施例】図１は、この発明の一実施例を示すものであ
る。図１において、１は振動子本体、2-1 〜2-7 は振動
素子、３は超音波パルス波面、４は振動素子を切り換え
るための先端部マルチプレクサ、５は信号発生器、６は
パルス発生器、７は送信アンプ、８は受信増幅回路、９
は受信増幅回路８の増幅度を制御するための STCコント
ロール回路、10はＢＰＦ、11は90度移相器、12a, 12bは
乗算器、13a, 13bはＬＰＦ、14a, 14bはＡ／Ｄコンバー
タ、15a, 15bはメモリを切り換えるためのメモリ用マル
チプレクサ、16は先端部マルチプレクサ４およびメモリ
用マルチプレクサ15a, 15bを送信毎に切り換えるための
切り換えコントロール回路、17は波面メモリ、18はアド
レス制御回路、19は高速フーリエ変換回路(FFT) 、20は
速度演算回路をそれぞれ示す。なお、波面メモリ17は、
振動素子2-1〜2-7 に対応して、実数部波面メモリA1〜A
7および虚数部波面メモリB1〜B7を具える。

【００１５】この実施例では、まず、切り換えコントロ
ール回路16の指示により先端部マルチプレクサ４とメモ
リ用マルチプレクサ15a, 15bを、振動素子2-1 および波
面メモリA1, B1に切り換え、この状態で信号発生器５の
出力によりパルス発生回路６、送信アンプ７および先端
部マルチプレクサ４を経て振動素子2-1 を駆動し、超音
波パルス３を発生させる。この超音波パルス３による成
体組織および血球からの反射波は、振動素子2-1 で電気
信号に変換し、その受信信号を先端部マルチプレクサ４
を経て受信増幅回路８に供給して適正な大きさに増幅す
る。この時、遠距離部になるほど信号が小さくなるの
で、STC コントロール回路９の制御により時間と共に増
幅度を大きくする。

【００１６】受信増幅回路８の出力は、BPF10 により余
分な雑音を取り除いた後、次段の乗算器12a, 12b、LPF1
3a, 13b および移相器11から成る直交検波回路に供給し
て、入力信号を(1) 式で示すベースバンド帯の複素数信
号ｇ(t, x)に変換する。これら実数部と虚数部とに分離
された複素数信号は、それぞれＡ／Ｄコンバータ14aお
よび14b によってデジタル信号に変換した後、マルチプ
レクサ15a, 15bを経て実数部波面メモリA1および虚数部
波面メモリB1に時系列データとしてそれぞれ格納する。

【数１】ｇ(t,x) ＝ａ(t,x) ＋ｊｂ(t,x) ------- (1) ここに、ａ(t,x) は実数部データ、ｂ(t,x) は虚数部デ
ータを示し、ｘは振動素子の空間位置を示す。

【００１７】次に切り換えコントロール回路16の制御下
に、先端部マルチプレクサ４およびメモリ用マルチプレ
クサ15a, 15bにより、振動素子2-1 から振動素子2-2
へ、波面メモリA1, B1から波面メモリA2, B2へ切り換え
て、同様にして送信・受信を行い、その組織および血流
からの反射信号を時系列データとして波面メモリA2, B2
へ格納する。

【００１８】以上の処理を繰り返して、全ての振動素子
2-1 〜2-7からの波面データを波面メモリ17に格納した
ら、その格納した波面データを開口合成法により処理し
て、各空間における超音波反射像を再現する。図２は、
開口合成法を説明するための図で、21a, 21bは超音波反
射体、22は波面メモリ、23は波面合成した結果を示すも
のである。ここで、振動素子2-1 から超音波パルスが送
波されると、反射体21a, 21bによって超音波パルスが反
射され、波面信号24-1として波面メモリ22に格納され
る。同様に、他の振動素子 (例えば2-7)から超音波パル
スが送波されると、その波面信号 (例えば24-7) が波面
メモリ22に格納される。これらの波面信号24-1や24-7
は、各振動素子2-1, 2-7と反射体21a, 21bの空間相対位
置関係が異なるために、それぞれ異なる波形となる。

【００１９】したがって、ある任意の空間において、各
振動素子2-1 〜2-7 で焦点を作るように、波面メモリ22
内で伝搬時間の遅れ進みを調整して加算すれば、従来の
遅延回路によるビーム合成と同様に、方位指向性をもた
せることができる。これは一般に“Dilay and Sum ”と
呼ばれている。この合成する波面を、全ての空間を対象
にして合成すれば、反射体の分布の様子が得られる。

【００２０】ここで、特筆すべきことは各振動素子2-1
〜2-7 単体での指向特性が広く、十分な方位分解能が得
られなくても、波面合成することにより振動子本体１の
もつ指向特性が鋭くなり、方位分解能が改善されること
である。これは、見方をかえれば、開口の小さな振動子
を大きい空間を移動しながら波面信号を検出し、それら
を全て用いて合成すれば、あたかも開口の大きい振動子
を用いたときと同じように鋭い指向特性が得られるもの
で、このことから開口合成といわれるゆえんである。

【００２１】次に、反射体の動きを検出する原理を、図
３を用いて説明する。図３(a) は反射体21が静止してい
る場合を示すもので、前述した送受信のシーケンスによ
り波面メモリ22に一連の時系列データが格納される。そ
の時系列データを、反射体21の空間位置に応じた遅延時
間で・印25に示すようなサンプリングポイントでサンプ
リングすると、実線26で示すようなある瞬間における空
間データが得られる。次に、受信波の中間周波数におけ
る位相が90度遅れたところでサンプリングすると、破線
27で示すような空間データが得られる。この場合、これ
ら実線26および破線27のデータは、一組のペアとして複
素数信号として表すことができる。この空間データの周
期的な変化は、空間周波数と呼ばれるもので、この場合
はほとんど直流成分のみとなる。この空間データの空間
周波数成分を検出するために、図１に示す FFT19でフー
リエ変換すれば、Ｐ(f) に示すような空間スペクトルが
得られる。このスペクトルは、ある空間の１ポイントの
データを示すものである。

【００２２】図３(b) は反射体21が近づいている場合を
示したもので、送受信のシーケンスが振動素子2-1 から
振動素子2-7 へ変わるのに従い、反射体21と振動素子2-
1 〜2-7 の距離が図３(a) の場合に比べて、だんだんと
短くなるため、波面メモリ22のデータも時間的にわずか
ずつずれた信号となる。この波面データを図３(a) の場
合と同じサンプリング時刻（・印25で示す) でサンプリ
ングすると、周期成分を持った信号が検出される。これ
を図１の FFT19でフーリエ変換すると、正の周波数の所
にピークを持った空間スペクトルが得られる。

【００２３】図４は、この波面データが少しずつずれる
現象を分かりやすく説明したもので、上から順に反射体
21が振動子２に少しずつ近づいている状態を示したもの
である。反射体21と振動子２の距離が短くなるに従い、
受信波が早く到達するようになる。したがって、ある時
刻でサンプリングすれば、反射体21の移動の様子が位相
の変化となって現れてくるので、このサンプリングを90
度位相が異なる２箇所で行えば、複素数信号として表す
ことができる。

【００２４】図３(c) は反射体21が振動素子列 (2-1 〜
2-7)と、すなわち超音波ビーム方向と直交する方向に平
行に動いている場合を示したものである。この場合に
は、送受信のシーケンスが、振動素子2-1 から振動素子
2-3 までは、反射体21の動きと共に空間距離が短くな
る。これに対し、振動素子2-4 から振動素子2-7 までの
送受信のシーケンスでは、反射体の動きと共に空間距離
が長くなる。このため、図３(a) の場合と同じ時刻でサ
ンプリングすると、正および負の周波数成分を持った信
号成分が検出される。これら正および負の周波数成分
は、複素数化したサンプリングデータによって見分けが
つくので、このデータをもとに図１の FFT19でフーリエ
変換を行えば、正の成分と負の成分の両方を持ったスペ
クトルが検出される。

【００２５】なお、反射体21が振動素子列から遠ざかる
場合は、図３(b) の場合と反対となり、フーリエ変換の
結果、負のスペクトル成分が検出される。また、反射体
21が図３(c) の場合と反対方向に振動素子列と平行に動
く場合は、フーリエ変換の出力である移相スペクトルが
図３(c) の場合と異なるため、その方向も検出できる。

【００２６】以上の説明から明らかなように、ある焦点
の反射信号を構成する複素信号は、図１に示すように入
力信号を乗算器12a, 12b、LPF13a, 13b 、移相器11を用
いた直交検波回路で検出できる他、相対的に位相が90度
ずれたところでサンプリングして検出することもでき
る。

【００２７】図５は、開口合成法の波面合成により、超
音波ビームの指向特性が変化する様子を概念的に説明す
るもので、(a) は真正面方向、(b)および(c) の左右の
方向に指向特性を作る場合をそれぞれ示したものであ
る。28a, 28bおよび28c は、データをサンプリングする
ための時刻を示す整相波面である。また、29a, 29bおよ
び29c は、サンプリングデータを距離に応じて格納する
ためのメモリで、この図では整相波面の数に応じてメモ
リが用意されているが、実際にはもっとたくさんにな
る。メモリ22内の整相波面をこのようにいろいろと変え
ることで、任意の方向に指向特性を得ることができる。

【００２８】この実施例では、図１に示すアドレス制御
回路18により波面メモリ17のアドレスを制御して上記の
整相波面を得、その波面データをもとに反射映像再生処
理を行う。すなわち、波面メモリ17に振動素子2-1 〜2-
7 に対応する全ての波面データを格納したら、アドレス
制御回路18で、ある方向の整相波面を構成するように時
間軸を調整して波面メモリ17のアドレスを制御し、これ
によりその整相波面の波面データｇ(x) を読み出す。こ
の波面データｇ(x) は、FFT19 に供給し、ここで、次式
により空間周波数を示すスペクトル分布Ｇ(w) を求め
る。

【数２】ここで、l は、振動素子で構成するアレイ長を、Ｐ(w)
はパワースペトクルを、θ(w) は位相スペクトルを示
す。

【００２９】波面データｇ(x) は、上述した説明から明
らかなように、実数部データA1〜A7（ａ(x) ）と、虚数
部データB1〜B7（ｂ(x) ）とに分かれているので、FFT1
9 によりスペクトル分布Ｇ(w) を求めることにより、こ
のスペクトル分布Ｇ(w) の実数部データＡ(w) と、虚数
部データＢ(w) とから、上記の(4) 式および(5) 式をも
とにパワースペクトルＰ(w) および位相スペクトルθ
(w) を求めることができる。

【００３０】FFT19 の出力であるスペクトル分布は、図
６に示すように、静止組織からの反射信号（クラッタ信
号）と血流からの反射信号とが重なったものとなる。こ
こで、静止組織からの反射信号は常に直流成分となるの
で、これをＢモード像の信号源として出力する。また、
このスペクトル分布は、速度演算回路20に供給して周波
数弁別により血流信号を取り出したり、周波数スペクト
ル分布の形状により血流の大きさおよび方向を求めて、
ドップラ信号として出力する。ここで、血流速度の大き
さは、スペクトルの平均周波数に対応し、血流の方向は
図３を用いて説明したように、そのパワースペクトル分
布Ｐ(W) や位相スペクトルθ(W) の形状によって決定さ
れる。なお、速度演算回路20でのスペクトル分布の周波
数弁別は、容易に行うことができるので、血流の信号成
分を減算させることなく、クラッタ信号を含まない血流
信号のみを取り出すことができる。

【００３１】速度演算回路20の出力であるドップラ信号
は、その情報に基づいて彩色して図示しないデジタルス
キャンコンバータ等により画像処理した後、Ｂモード像
に重畳させて表示する。

【００３２】なお、この発明は上述した実施例にのみ限
定されるものではなく、幾多の変形または変更が可能で
ある。例えば、上述した実施例では、振動子本体の振動
素子の数を７個として説明したが、実際にはもっとたく
さんの振動素子数（例えば、64個程度) で振動子本体を
構成する。このように、振動素子数を増加すれば、ビー
ム指向特性が鋭くなるので、空間分解能を向上すること
ができる。また、あるポイントの反射映像を作るための
全体の送受信シーケンスに要する時間が長くなることか
ら、血球のゆっくりした動きを検出することができる共
に、送受信シーケンスが増えることで、データ数も増え
るので、ドップラの測定精度も向上することができる。

【００３３】さらに、上述した実施例では、FFT を用い
て空間周波数のスペクトルを求めたが、他のスペクトル
分布検出器を用いて空間周波数のスペクトルを求めるこ
ともできる。また、振動素子の配列は、直線上に限ら
ず、円周配列や、直線配列と円周配列との中間である通
称“コンベックス”と呼ばれているアレイ型配列とする
こともできる。このように振動素子の配列が変わった場
合には、その振動素子配列が所望の焦点を構成するよう
に、図１のアドレス制御回路18によって時間調整による
波面整合を行えばよい。

【００３４】

【発明の効果】以上のように、この発明によれば、各振
動素子について得られた波面データをメモリに蓄え、こ
れらの波面データに基づいて空間周波数のスペクトル分
布を求めるようにしたので、その空間周波数の直流成分
から組織の反射強度を、それ以外の周波数成分から血流
速度を求めることができる。また、組織からの反射信号
と血流速度信号は、単にフーリエ変換の結果を用いて周
波数弁別することにより容易に分離することができるの
で、低周波成分の減衰もなく、かなり低速血流まで検出
することができる。さらに、空間周波数のスペクトル分
布の形状を求めることにより、超音波ビーム方向のみの
血流成分だけでなく、ビームと直交する方向についても
血流速度を求めることができる。

【００３５】また、従来の超音波ビーム合成法では、各
振動素子の出力を用いてビーム合成し一括してそのビー
ム内の血流成分を求めていたので、特に至近距離におい
て測定精度が劣化するという問題があったが、この発明
では１個１個の振動素子と血球の動きの相対関係によっ
て空間周波数が決定されるので、至近距離においても血
球の動きが厳密に空間周波数に反映され、精度よく血流
速度を求めることができる。

【００３６】さらに、この発明では、スペクトル分布検
出器の出力である直流成分からＢモード像を求め、それ
以外の周波数成分から血流速度成分を求めるようにして
いるので、特別に血流速度を求めるための送信シーケン
スを設ける必要がない。また、通常の構成では、あるポ
イントの血球に対して数10回の反射エコーをもとに波面
合成して血流速度成分を検出するようにしているが、こ
の発明ではこのような平均操作が不要になる。したがっ
て、従来のＢモード像と比べてフレームレイトが低下す
ることがないと共に、カラードップラのモードにおいて
も応答性のよい血流速度分布を映像化することができ
る。さらに、また、多数の送信シーケンスによって血流
速度を求めるようにしているので、最低測定速度が小さ
くなり、流速の小さい血管の観察も容易にできる。この
ように、流速の小さい血流も測定でき、また振動素子列
と平行（プローブの振動面と平行）な血流も測定できる
ことから、腹部臓器を対象とした診断にきわめて有用に
なる。

【図面の簡単な説明】

【図１】この発明の一実施例を示す図である。

【図２】開口合成の概念を説明するための図である。

【図３】反射体の動きを検出する原理を説明するための
図である。

【図４】反射体の動きによって位相波面がずれることを
説明するための図である。

【図５】指向特性を変化させる概念を説明するための図
である。

【図６】周波数弁別によるクラッタ信号と血流信号の分
離を説明するための図である。

【図７】従来の技術を説明するための図である。

【図８】従来例におけるクラッタ信号と血流信号の分離
を説明するための図である。

【符号の説明】

１振動子本体 2-1 〜2-7 振動素子４先端部マルチプレクサ８受信増幅回路 11 90度移相器 12a, 12b 乗算器 13a, 13b ローパスフィルタ(LPF) 14a 14b Ａ／Ｄコンバータ 15a 15b メモリ用マルチプレクサ 17 波面メモリ 18 アドレス制御回路 19 高速フーリエ変換器(FFT) 20 速度演算回路

Claims

(57)【特許請求の範囲】

【請求項１】アレイ状に並べた振動素子群と、それら
を切り替えるためのマルチプレクサと、このマルチプレ
クサの出力を増幅する受信増幅回路と、その出力を複素
数信号としてデジタル信号に変換するＡ／Ｄコンバータ
と、該デジタル信号を各振動素子に対応させて、複素数
の波面データとして格納する波面メモリと、これら波面
メモリに格納した複素数の波面データに基づいて空間周
波数を求めるスペクトル分布検出器とを具え、このスペ
クトル分布検出器で求めた空間周波数の直流成分から超
音波反射強度を求め、空間周波数の直流以外の周波数成
分から血流速度を求めるよう構成したことを特徴とする
超音波診断装置。