JP4542258B2

JP4542258B2 - 超音波診断装置

Info

Publication number: JP4542258B2
Application number: JP2000388435A
Authority: JP
Inventors: 好一宮坂; 烈光原田
Original assignee: Aloka Co Ltd
Current assignee: Hitachi Ltd
Priority date: 2000-12-21
Filing date: 2000-12-21
Publication date: 2010-09-08
Anticipated expiration: 2020-12-21
Also published as: JP2002186615A

Description

【０００１】
【発明の属する技術分野】
本発明は超音波診断装置に関し、特に超音波画像に含まれるノイズの低減に関する。
【０００２】
【従来の技術及びその課題】
超音波診断装置は生体への超音波の送受波により超音波画像を表示する装置である。超音波画像としては、二次元白黒断層画像であるＢモード画像が周知である。かかるＢモード画像内には、生体内の構造物とは関係のない超音波波面の干渉によって生じる粒状模様のスペックルというアーチファクトが含まれる。そのようなノイズはＢモード画像の観察において障害となる無用なデータである。なお、そのようなスペックルあるいはアーチファクトはＭモード画像などの他の超音波画像においても認められる。
【０００３】
本発明は、上記従来の課題に鑑みなされたものであり、その目的は、超音波画像の画質を向上することにある。本発明の他の目的は、受信信号自体を利用してそれに含まれるノイズ成分を特定することにある。
【０００４】
【課題を解決するための手段】
（１）上記目的を達成するために、本発明は、超音波を送受波し、受信信号を出力する送受波手段と、前記受信信号から、第１周波数帯域の第１信号成分及び第２周波数帯域の第２信号成分を抽出する成分抽出手段と、前記第１信号成分から第１包絡線データを求める第１検波手段と、前記第２信号成分から第２包絡線データを求める第２検波手段と、前記第１包絡線データ及び前記第２包絡線データの相互相関演算を行って相互相関値を演算する相互相関演算手段と、前記相互相関値に応じて受信信号に対するフィルタリングを行うフィルタ手段と、を含むことを特徴とする。
【０００５】
上記構成によれば、検波後のエンベロープ信号としての第１包絡線データと第２包絡線データについて相互相関演算が実行され、その相互相関値に基づいて受信信号（エコーデータあるいは画像データ）に対するフィルタリングがなされる。これは生体組織からの真のエコーの場合には相互相関値が比較的大きくなり、ノイズの場合には相互相関値が比較的小さくなるという傾向を利用したものである。よって、本発明によれば、ノイズを選択的に除去、抑圧できるので、超音波画像の画質を高めることが可能となる。フィルタリングは一次元的に行うようにしてもよいが、二次元的に行うようにしてもよい。
【０００６】
望ましくは、前記第１包絡線データ及び前記第２包絡線データは時間軸上のデータであり、前記相互相関演算手段は時間軸上において第１包絡線データ及び第２包絡線データの相互相関演算を実行する。
【０００７】
望ましくは、前記第１包絡線データ及び前記第２包絡線データは周波数軸上のデータであり、前記相互相関演算手段は周波数軸上において第１包絡線データ及び第２包絡線データの相互相関演算を実行する。
【０００８】
望ましくは、前記フィルタ手段はローパスフィルタであり、前記相互相関値が大きい場合には前記ローパスフィルタのカットオフ周波数を高く設定し、前記相互相関値が小さい場合には前記ローパスフィルタのカットオフ周波数を低く設定するフィルタ特性設定手段が設けられる。この構成によれば、真のエコーデータに対するフィルタリングをできるだけ回避してノイズに対して選択的かつ能率的にフィルタリングを行える。もちろん、カットオフ周波数の可変の他、通過帯域の可変やゲイン調整を行うのもフィルタリングの一態様である。
【０００９】
（２）また、上記目的を達成するために、本発明は、超音波を送受波し、受信信号を出力する送受波手段と、前記受信信号から第１周波数帯域の第１信号成分及び第２周波数帯域の第２信号成分を抽出する成分抽出手段と、前記第１信号成分から第１包絡線データを求める第１検波手段と、前記第２信号成分から第２包絡線データを求める第２検波手段と、信号時系列に沿って一次元のウインドをスキャンさせながら、各ウインド位置において切り出される第１包絡線データ及び第２包絡線データの相互相関演算を行って相互相関値を演算する相互相関演算手段と、前記相互相関値に応じて受信信号に対するフィルタリングを行うフィルタ手段と、を含むことを特徴とする。
【００１０】
（３）また、上記目的を達成するために、本発明は、超音波ビームを走査し、各ビーム位置ごとに受信信号を出力する送受波手段と、前記受信信号から第１周波数帯域の第１信号成分及び第２周波数帯域の第２信号成分を抽出する成分抽出手段と、前記第１信号成分から第１包絡線データを求める第１検波手段と、前記第２信号成分から第２包絡線データを求める第２検波手段と、超音波ビーム上の深さ方向と超音波ビームの走査方向とによって定義される二次元座標系上において、二次元のウインドをスキャンさせながら、各ウインド位置において切り出される第１包絡線データ及び第２包絡線データの相互相関演算を行って相互相関値を演算する相互相関演算手段と、前記相互相関値に応じて受信信号に対するフィルタリングを行うフィルタ手段と、を含むことを特徴とする。この構成によれば、二次元的に相関演算を行えるので、よりノイズ判定精度を高めることができる。
【００１１】
（４）以下に、本発明の原理について説明する。超音波の送受波により得られた受信信号に含まれるノイズ（特にスペックル）を特定するため、その受信信号から互いに異なる複数の（望ましくは２つの）周波数帯域の信号成分が抽出される。
【００１２】
図１には、受信信号のスペクトル３００が示されており、符号３０２はその全体域を抽出する場合のバンドバスフィルタの通過帯域特性を示し、符号３０４及び符号３０６は、それぞれ上記のように２つの周波数帯域の信号成分を抽出するためのバンドパスフィルタの通過帯域を示している。符号３０４で示す通過帯域特性は受信信号のスペクトル３００における低域側に設定され、符号３０６で示す通過帯域特性は受信信号のスペクトル３００における高域側に設定されている。それらの通過帯域特性は実質的に分離され、一部の裾が相互にオーバーラップしている。もちろん、全体として２つの通過帯域特性に違いがあれば、各種の通過帯域特性を設定することが可能である。
【００１３】
さて、それぞれの周波数帯域の信号成分の包絡線データ相互間において相関演算を実行すると、生体内の構造体からの真のエコー成分の包絡線の形状は互いに似ているため相互相関値が大きくなり、一方、ノイズの包絡線の形状はランダムであるために相互相関値が小さくなる。よって、そのような相互相関値の大きさから、個々の受信信号が真のエコー成分（あるいはノイズ）である確率を求めることができ、あるいは、真のエコー成分かノイズかの判別を行うことができる。
【００１４】
より詳しくは、、例えば、まず受信信号から２つの周波数帯域の信号成分（データ列）が抽出され、各信号成分が検波されて包絡線（エンベロープ）データ（データ列）とされ、それぞれのデータ列が１ライン分（あるいは１フレーム分）を単位としてメモリに個別的に格納される。そして、それぞれのデータ列から、所定のウインド（図２参照）内のＮ個のデータからなる２組のデータ列（ｘ_i，ｙ_i）（但し、ｉ＝１，２，・・・，Ｎ）が切り出され、そのような２組のデータ列（ｘ_i，ｙ_i）について相互相関演算がなされる。ちなみに、図２には、超音波ビーム３０８を電子走査することにより形成される走査面３０６が示され、上記のデータ列の切り出しは、各超音波ビーム３０８上において一次元のウインド３１０を設定することに相当する。なお、信号成分の抽出をデータ列の切り出し後に行うようにしてもよい。
【００１５】
相関演算は、相互相関値をＲｘｙとして以下の数式によって表される。
【００１６】
【数１】

この相互相関値Ｒｘｙは、互いのデータ列がまったく同一であれば１となり、まったく異なれば０となる。つまり、違いが大きい程、相互相関値は小さくなる。ここで、スペックル（ノイズ）の場合には、相互相関値は非常に小さくなる。
一方、生体内の構造体からのエコー（真のエコー）の場合には、相互相関値は大きくなる。よって、相互相関値を基準としてスペックルか真のエコーかを弁別することが可能であり、更に、その相互相関値にフィルタ特性を連動させれば、スペックルに対して選択的に抑圧処理を行うことができる。例えば、真のエコーである場合には、当該データが通過するローパスフィルタ（ＬＰＦ）のカットオフ周波数を高くして当該データをそのまま通過させ、一方、スペックルである場合には、ＬＰＦのカットオフ周波数を低くして当該データが抑圧、除去されるようにする。
【００１７】
あるウインド位置において上記相互相関演算が実行されると、ウインド位置がシフトされて上記同様の相互相関演算が実行される。そして、これが繰り返されると、全データに対してノイズ低減処理を行える。
【００１８】
上記のウインドは、深さ方向（超音波ビーム方向）に沿った１次元のウインドであったが、図３に示すように、二次元のウインドであってもよい。図３において、Ｊは深さ方向（超音波ビーム方向）、Ｉは超音波ビームの走査方向あるいはＢモード画像の水平方向に対応している。二次元マトリクス３１２上において、二次元のウインド３１４がラスタースキャンされ、各スキャン位置においてデータ列が切り出される。ウインド３１４のサイズはＷＩ×ＷＪであり、その場合、以下の演算式によって相互相関値Ｒｘｙが求められる。
【００１９】
【数２】

上記の構成を採用する場合、受信信号を処理するフィルタ手段として、二次元のＬＰＦを利用するようにしてもよい。相互相関演算は時間軸上の他、周波数軸上で行うようにしてもよい。
【００２０】
なお、受信信号スペクトルは、周波数依存減衰などの影響により、反射点の深さに依存して変化するため、それに応じて図１の通過帯域特性３０４及び３０６を動的に変化させるようにするのが望ましい。
【００２１】
【発明の実施の形態】
以下、本発明の好適な実施形態を図面に基づいて説明する。
【００２２】
図４には、本発明に係る超音波診断装置の好適な実施形態が示されており、図４はその全体構成を示すブロック図である。
【００２３】
図４において、アレイ振動子１０は、図示されていない超音波探触子内に設けられるものであり、そのアレイ振動子は複数の振動素子によって構成される。そのアレイ振動子上において超音波ビームが電子走査され、これによって図２に示したような走査面３０６が形成される。
【００２４】
送信器１２はいわゆる送信ビームフォーマーとして機能するものであり、アレイ振動子を構成する複数の振動素子に対して送信信号を供給する回路である。受信器１３は、いわゆる受信ビームフォーマーとして機能するものであり、アレイ振動子を構成する複数の振動素子からの複数の受信信号に対して整相加算処理を実行し、整相加算後の受信信号を出力する回路である。この受信信号は本実施形態において３つの信号処理部１４，１６，１８に並列的に入力されている。
【００２５】
まず、信号処理部１４について説明すると、この信号処理部１４は従来の超音波診断装置においても一般的に設けられているものであり、受信信号に対して直交検波などの検波処理を実行する検波器２０と、その検波後の受信信号（複素信号）に対して振幅演算を実行する振幅演算器２２と、演算された振幅値すなわちエコーデータを超音波ビーム１本分ごとに格納するラインメモリ部２４と、で構成されている。検波器２０の具体的な構成例については後に図５及び図６を用いて説明する。また、振幅演算器２２については後に図７を用いて説明する。ラインメモリ部２４は、本実施形態において２つのラインメモリによって構成され、いわゆるデータの書き込み及び読み出しに関してそれらの２つのメモリを交互に使用するピンポン動作を行わせている。
【００２６】
フィルタ４２は本実施形態においてローパスフィルタを構成しており、カットオフ周波数が後に詳述する相互相関値によって動的に設定される。フィルタ４２から出力されたデータはスキャン変換器４４に入力され、ここで計測座標系から表示座標系への座標変換などの処理を経た後、そのデータが表示器４６に表示される。具体的には、表示器４６にはＢモード画像などの超音波画像が表示される。ちなみに、スキャン変換器４４は例えばデジタルスキャンコンバータなどによって構成される。
【００２７】
次に、信号処理部１６，１８について説明する。これらの信号処理部１６，１８は上述した信号処理部１４と同様の構成を有しており、すなわち検波器２６，２８と、振幅演算器３０，３２と、ラインメモリ部３４，３６と、によって構成されている。ただし、検波器２０，２６，２８はそれぞれ互いに異なる通過帯域特性をもったフィルタを有している。具体的に説明すると、検波器２０は、図１において符号３０２で示した全周波数帯域をカバーする通過帯域特性をもったフィルタを有している。検波器２６は図１において符号３０４で示した低域側の通過帯域特性をもったフィルタを備えている。また、検波器２８は図１において符号３０６で示した高域側の通過帯域特性をもったフィルタを備えている。各検波器はＲＦ信号を包絡線信号に変換する機能を有する。
【００２８】
したがって、振幅演算器２２では、受信信号の全周波数帯域における振幅が演算され、一方、振幅演算器３０においては低域側の信号成分の振幅が演算され、振幅演算器３２においては高域側の信号成分の振幅が演算されることになる。ラインメモリ部３４，３６には以上のようにして演算された振幅値をもったデータ列が超音波ビーム１本分ごとに格納される。
【００２９】
そして、本実施形態においては、ラインメモリ部３４，３６から、上述したようにＮ個のデータからなるデータ列が抽出され、それらの２組のデータ列に対して相互相関器４０において相互相関演算が実行される。ここで、ラインメモリ部３４，３６からのＮ個のデータの切り出しは上述したように所定のウインドの設定によってなされており、そのウインドをデータ列方向に沿ってスキャンすることにより、各ウインド位置において相互相関演算が実行される。
【００３０】
上述の説明から理解されるように、図４に示す相互相関器４０は、時間軸上において２組のデータ列間で相互相関演算を実行するものであり、その相互相関の演算式としては上述した式（１）などの数式が用いられる。その相互相関演算結果である相互相関値はフィルタ４２に出力され、その相互相関値にフィルタ４２におけるフィルタ特性が適応的に設定されることになる。これにより、ある深さにおけるエコーデータに対して、相互相関値に応じてフィルタリングを行うことができ、そのエコーデータがノイズであればそのフィルタリングによって当該ノイズを効果的に除去・低減することが可能となる。
【００３１】
図５には、図４に示した検波器２０，２６，２８の具体的な構成例が示されている。この図５に示す構成例はいわゆる直交検波器を示すものである。
【００３２】
２つのミキサ７０，７２の一方の入力端子には受信信号が入力され、他方の入力端子には参照信号が入力されている。ミキサ７２はπ／２位相シフタ６８を介して位相角が９０度ずれた参照信号が入力されており、一方、ミキサ７２は入力される参照信号がそのまま入力されている。このような２つの参照信号の位相関係によって直交検波を行うことができ、その直交検波後の信号のうちベースバンド領域の信号成分が低域通過フィルタ７４，７６によって抽出されている。この場合において、図４に示した各検波器２０，２６，２８ごとに低域通過フィルタ７４，７６の通過帯域特性が設定されており、すなわちそれぞれの通過帯域特性は図１において符号３０２，３０４，３０６で示したものが設定されている。
【００３３】
図６には、検波器２０，２６，２８の他の構成例が示されており、この図６に示す構成例はいわゆる直交サンプリング回路である。帯域通過フィルタ７８には受信信号が入力され、その帯域通過フィルタ７８を通過した受信信号がサンプリング回路８２，８４に並列的に入力される。ここで、サンプリング回路８２は、サンプリングクロックに従ってサンプリングを行うものであり、サンプリング８４はπ／２シフタ８０を介して９０度位相がずれたサンプリングクロックに従ってサンプリングを実行している。このような直交サンプリング方式により周知のように直交検波と同様の検波信号を得ることが可能となる。もちろん、検波器２０，２６，２８の相互間においては帯域通過フィルタ７８の帯域通過特性としてそれぞれ異なる特性が設定されており、具体的にはそれぞれの通過帯域特性として図１において符号３０２，３０４，３０６で示したものが設定されている。
【００３４】
図７には、振幅演算器２２，３０，３２の構成例が示されている。絶対値演算器８６は、検波器２０から出力される複素信号における実数部の二乗及び虚数部の二乗を加算し、その加算結果の平方根をとることによって絶対値を求めている。これによってエコーデータの振幅を演算することが可能となる。
【００３５】
その振幅を表すデータは対数変換器８８に入力され、その対数変換器８８において対数圧縮処理がなされる。そして、その対数圧縮後のデータに対して加算器９０に係数βが加算され、また乗算器９２において係数αが乗算される。ここで加算器９０はゲイン調整回路として機能し、乗算器９２はコントラスト回路として機能するものである。乗算器９２から出力されるデータはローパスフィルタ９４に入力され、ローパスフィルタ９４直後のデシメータ９６によるリサンプリング処理のために帯域が制限され、ローパスフィルタ９４から出力されたデータがデシメータ９６に入力され、そのデシメータ９６において表示ピクセルレートに従ったリサンプリング処理が実行される。
【００３６】
図８には、ラインメモリ部３４，３６に格納されるデータ列が概念的に示されている。ここで符号３０４Ａはラインメモリ部３４に格納されたデータ列であり、このデータ列は図１において通過帯域特性３０４によって抽出された信号成分に相当する。また図８において符号３０６Ａは、ラインメモリ部３６に格納されたデータ列を示しており、これは図１において通過帯域特性３０６によって抽出された信号成分に相当する。
【００３７】
上述したように、このような２つのデータ列からウインドによって２組のデータ列が部分的に切り出され、そのウインド内のデータ列に対して相互相関演算が実行される。
【００３８】
図９には、図４に示したフィルタ４２の構成例が示されている。
【００３９】
フィルタ特性選択テーブル９８には、相互相関値が入力されており、その相互相関値に応じたフィルタ特性を表す選択信号がフィルタ係数テーブル１００に出力される。フィルタ係数テーブル１００は例えば図１０に示すようなテーブル内容を有しており、具体的には、ローパスフィルタとしての当該フィルタ４２のカットオフ周波数を相互相関値に応じて変化させるために、複数のフィルタ特性に対応付けられたフィルタ係数列が格納されている。上記のフィルタ特性の選択信号が入力されると、フィルタ係数テーブル１００から当該選択信号によって選択選択されるフィルタ特性を表すフィルタ係数列ａ₁〜ａ_nが出力され、それぞれのフィルタ係数がそれに対応する乗算器１０４−１〜１０４−ｎの一方の入力端子に入力される。受信信号は直列接続された複数のディレイライン１０２−１〜１０２−ｎに順番に入力されており、それぞれのディレイラインの前後及び中間のポイントから出力されたデータが各乗算器１０４−１〜１０４−ｎの他方の入力端子に入力されている。そして、各乗算器１０４−１〜１０４−ｎはそれぞれのデータに対してフィルタ係数を乗算し、その乗算結果が加算器１０６において加算されている。
【００４０】
図１１には、相互相関値とローパスフィルタのカットオフ周波数との関係が示されており、本実施形態においては相互相関値が小さくなるほどカットオフ周波数が低められており、一方、相互相関値が高くなるほどカットオフ周波数が高められている。図においてはＳ字関数のような周波数変化特性が示されている。
【００４１】
図１２には、他の実施形態に係る超音波診断装置の全体構成がブロック図として示されている。なお、図４に示した構成と同様の構成には同一符号を付しその説明を省略する。
【００４２】
この図１２に示す構成例においては各信号処理部５０，５２，５４において振幅演算器２２，３０，３２の後段にスキャン変換器５６，６０，６２及びフレームメモリ部５８，６４，６６が設けられている。ここで、各スキャン変換器５６，６０，６２はそれぞれ図４に示したスキャン変換器４４に相当するものであり、すなわちスキャン変換器５６，６０，６２において送受波座標系から表示座標系への座標変換がなされる。そして、そのデータはそれぞれフレームメモリ部５８，６４，６６に１フレームごとに格納されることになる。
【００４３】
このようにフレームメモリ部６４，６６へのデータの格納がなされると図３に示した二次元のウインド３１４が利用され、そのウインド３１４をラスタースキャンさせながら各ウインド位置においてデータ列が抽出され、相互相関器４０においては２組の二次元のデータ列の相互相関値が演算される。そしてその相関値に従って、フィルタ４２がフレームメモリ部５８から出力されるデータに対してフィルタリングを実行している。
【００４４】
図１３には、図４の相互相関器４０の周波数軸上での処理構成例を示す。この構成例では、２つの一次元データ列に対してそれぞれＦＦＴ回路１１０，１１２においてＦＦＴ演算がなされ、これによりそれらのデータ列が周波数軸上のデータ列に変換される。それらのデータ列は乗算器１１３で乗算され、その乗算結果（相互相関値に相当）がＩＦＦＴ回路１１４に入力され、そこで逆ＦＦＴ演算が実行される。
【００４５】
図１４には図１２の相互相関器４０の周波数軸上での処理構成例を示す。この構成例では、２つの二次元データ列に対してそれぞれＦＦＴ回路１２８，１３０において二次元のＦＦＴ演算がなされ、これによりそれらのデータ列が周波数軸上のデータ列に変換される。それらのデータ列は乗算器１３１で乗算され、その乗算結果（相互相関値に相当）が二次元ＩＦＦＴ回路１３２に入力され、そこで二次元の逆ＦＦＴ演算が実行される。
【００４６】
以上のように、上記の各実施形態によれば、受信信号における２つの周波数成分から相互相関値を求め、その相互相関値に応じて受信信号に対するフィルタリングを実行することによりノイズ（特にスペックル）に対して効果的な除去処理を行うことが可能となる。上記の実施形態においては、２つの通過帯域特性を設定したが、それぞれの通過帯域特性が部分的にオーバーラップしていてもよい。
【００４７】
【発明の効果】
以上説明したように、本発明によれば、超音波画像内におけるノイズを効果的に低減して超音波画像の画質を向上できる。
【図面の簡単な説明】
【図１】受信信号のスペクトルと通過帯域特性の関係を示す図である。
【図２】走査面及び一次元のウインドを示す図である。
【図３】二次元のウインドを示す図である。
【図４】実施形態に係る超音波診断装置の全体構成を示すブロック図である。
【図５】検波器の構成例を示すブロック図である。
【図６】検波器の他の構成例を示すブロック図である。
【図７】振幅演算器の構成例を示すブロック図である。
【図８】ラインメモリ部に格納されるデータ列を示す図である。
【図９】フィルタの構成例を示すブロック図である。
【図１０】フィルタ係数テーブルの構成例を示す概念図である。
【図１１】相互相関値とカットオフ周波数との関係を示す図である。
【図１２】他の実施形態に係る超音波診断装置の全体構成を示すブロック図である。
【図１３】周波数軸上での一次元相互相関処理の構成例を示す図である。
【図１４】周波数軸上での二次元相互相関処理の構成例を示す図である。
【符号の説明】
１０アレイ振動子、１２送信器、１３受信器、１４，１６，１８信号処理部、２０，２６，２８検波器、２２，３０，３２振幅演算器、２４，３４，３６ラインメモリ部、４０相互相関器、４２フィルタ、４４スキャン変換器、４６表示器。

Claims

超音波を送受波し、受信信号を出力する送受波手段と、
前記受信信号から、第１周波数帯域の第１信号成分及び第２周波数帯域の第２信号成分を抽出する成分抽出手段と、
前記第１信号成分から第１包絡線データを求める第１検波手段と、
前記第２信号成分から第２包絡線データを求める第２検波手段と、
前記第１包絡線データ及び前記第２包絡線データの相互相関演算を行って相互相関値を演算する相互相関演算手段と、
前記相互相関値に応じて受信信号に対するフィルタリングを行うフィルタ手段と、
を含むことを特徴とする超音波診断装置。
請求項１記載の装置において、
前記第１包絡線データ及び前記第２包絡線データは時間軸上のデータであり、前記相互相関演算手段は時間軸上において第１包絡線データ及び第２包絡線データの相互相関演算を実行することを特徴とする超音波診断装置。
請求項１記載の装置において、
前記第１包絡線データ及び前記第２包絡線データは周波数軸上のデータであり、
前記相互相関演算手段は周波数軸上において第１包絡線データ及び第２包絡線データの相互相関演算を実行することを特徴とする超音波診断装置。
請求項１記載の装置において、
前記フィルタ手段はローパスフィルタであり、
前記相互相関値が大きい場合には前記ローパスフィルタのカットオフ周波数を高く設定し、前記相互相関値が小さい場合には前記ローパスフィルタのカットオフ周波数を低く設定するフィルタ特性設定手段が設けられたことを特徴とする超音波診断装置。
超音波を送受波し、受信信号を出力する送受波手段と、
前記受信信号から第１周波数帯域の第１信号成分及び第２周波数帯域の第２信号成分を抽出する成分抽出手段と、
前記第１信号成分から第１包絡線データを求める第１検波手段と、
前記第２信号成分から第２包絡線データを求める第２検波手段と、
信号時系列に沿って一次元のウインドをスキャンさせながら、各ウインド位置において切り出される第１包絡線データ及び第２包絡線データの相互相関演算を行って相互相関値を演算する相互相関演算手段と、
前記相互相関値に応じて受信信号に対するフィルタリングを行うフィルタ手段と、
を含むことを特徴とする超音波診断装置。
超音波ビームを走査し、各ビーム位置ごとに受信信号を出力する送受波手段と、
前記受信信号から第１周波数帯域の第１信号成分及び第２周波数帯域の第２信号成分を抽出する成分抽出手段と、
前記第１信号成分から第１包絡線データを求める第１検波手段と、
前記第２信号成分から第２包絡線データを求める第２検波手段と、
超音波ビーム上の深さ方向と超音波ビームの走査方向とによって定義される二次元座標系上において、二次元のウインドをスキャンさせながら、各ウインド位置において切り出される第１包絡線データ及び第２包絡線データの相互相関演算を行って相互相関値を演算する相互相関演算手段と、
前記相互相関値に応じて受信信号に対するフィルタリングを行うフィルタ手段と、
を含むことを特徴とする超音波診断装置。