JP4312202B2

JP4312202B2 - 超音波診断装置および超音波診断装置のデータ処理方法

Info

Publication number: JP4312202B2
Application number: JP2005506588A
Authority: JP
Inventors: 浩仲宮木
Original assignee: Olympus Corp
Current assignee: Olympus Corp
Priority date: 2003-05-29
Filing date: 2004-05-31
Publication date: 2009-08-12
Anticipated expiration: 2024-05-31
Also published as: EP1627602B1; EP1627602A4; US20060122510A1; US7803114B2; JPWO2004105614A1; EP1627602A1; WO2004105614A1

Description

本発明は、生体組織に対して超音波を送受波して生体内部の断層像を得る超音波診断装置のデータ処理方法及び超音波診断装置に関する。

従来、超音波振動子から生体組織に対して超音波を送受信して生体組織の断層像を得る超音波診断装置が種々実用化されている。
さらに、この種の超音波診断装置においては、超音波を用いて生体組織における血流の有無を感度良く表示する方法の１つにパワードプラ処理が知られている。このパワードプラ処理を行う超音波診断装置の一例としては、受信した反射エコー信号に直交検波を施し、この検波信号をＭＴＩ（ＭｏｖｉｎｇＴａｒｇｅｔＩｎｄｉｃａｔｏｒ）フィルタに通してからＩ信号とＱ信号の２乗和をとることで血流のパワーを求め、２次元画像として表示するものが知られている（例えば、特開昭６１−２５７６３１公報（第２−４頁、第１−５図））
上記ＭＴＩフィルタは、ＭＴＩ（ＭｏｖｉｎｇＴａｒｇｅｔＩｎｄｉｃａｔｏｒ）レーダーの原理を用いるもので、主として低い次数（ｏｒｄｅｒ）の無限衝撃応答（ＩＩＲ）が使われている。
本発明は、血流可視化を行うとともにノイズの抑制効果を向上することができる超音波診断装置のデータ処理方法及び超音波診断装置を提供することを目的とする。

簡略に本発明の超音波診断装置のデータ処理方法は、所定間隔で同一の音線方向に対する複数回の超音波送受信を行って得られた複数の超音波データの中から所定の同一の深さに対応するデータを得、得られた所定の同一の深さに対応するデータをフィルタリングして体内の運動部位からのエコーデータを抽出し、このエコーデータを第１所定時間遅延させた場合の第１自己相関の振幅値を得ると共に、同エコーデータを第２所定時間遅延させた場合の第２自己相関の振幅値を得て、さらに、上記第１、第２の自己相関の振幅値を掛け合わせて得られた乗算値を所定色相の輝度に対応させる。

図１：本発明の第１実施形態の超音波診断装置の全体構成を示す説明図。
図２：第１実施形態の超音波診断装置のデータ処理方法を示すフローチャート。
図３：第１実施形態の超音波診断装置における第１の自己相関演算部の動作を示すタイミングチャート。
図４：第１実施形態の超音波診断装置における第２の自己相関演算部の動作を示すタイミングチャート。
図５：第１実施形態の超音波診断装置における自己相関振幅値とホワイトノイズの関係を示す説明図。
図６：第１実施形態の超音波診断装置における超音波診断装置のデータ処理方法を更に具体的に示すフローチャート。
図７：第１実施形態の超音波診断装置における上記図６のステップＳ１５を詳細に示すフローチャート。
図８：第１実施形態の超音波診断装置における上記図６のステップＳ１６を詳細に示すフローチャート。
図９：本発明の第２実施形態の超音波診断装置の全体構成を示す説明図。
図１０：第２実施形態の超音波診断装置のデータ処理方法を示すフローチャート。
図１１：第２実施形態の超音波診断装置のデータ処理方法を更に具体的に示すフローチャート。
図１２：本発明の第３実施形態の超音波診断装置の全体構成を示す説明図。
図１３：第３実施形態の超音波診断装置のデータ処理方法を具体的に示すフローチャート。

以下、図面を参照して本発明の実施の形態を説明する。
（第１の実施形態）
（構成）
図１〜図８は本発明の第１実施形態の超音波診断装置を示したものである。
図１に示すように、超音波診断装置１は、超音波探触子１１と、超音波送受信部１２と、アナログ／デジタル変換器（以下、Ａ／Ｄ変換器と呼ぶ）１３と、直交変換器１４と、遅延器１５と、ＭＴＩフィルタ１６，１７と、第１及び第２の自己相関演算部１８，１９と、乗算器２０と、Ｂモード像演算部２１と、デジタルスキャンコンバータ（以下、ＤＳＣと呼ぶ）２２と、モニタ２３とを含んで構成されている。
超音波探触子１１は、超音波送受信部１２に接続され、超音波の送信及び受信を行う。
ここで、超音波探触子１１は、電気的な走査を行うアレイ振動子、または機械的な走査を行う単板あるいはアレイ振動子のいずれでも良い。
超音波送受信部１２は、超音波探触子１１にパルス信号を供給するようになっている。超音波探触子１１は超音波送受信部１２からパルス信号を供給されると超音波パルスの送信を行い、生体から戻ってくる反射エコーを受信する。
ここで、超音波探触子１１にアレイ振動子を用いた場合、超音波送受信部１２は、ビームフォーミング処理を行い、音線データを生成する。
超音波送受信部１２は、血流情報を取得する場合に、同一の音線に対して所定の時間間隔で複数回送受信を行うように動作する。本第１実施形態の場合、同一の音線に対する送受信の回数は６回程度である。各音線では、所定の同一深さのデータが送受信の回数分得られる。
超音波送受信部１２で生成された音線データは、Ａ／Ｄ変換器１３に入力され、所定のサンプリング周波数にてデジタルデータ化される。なお、超音波送受信部１２とＡ／Ｄ変換器１３の間には図示しないアンチエリアジングフィルタを備えており、Ａ／Ｄ変換時にエリアジングを防ぐ役割を果たす。
ここで、エリアジングとは、アナログ信号をＡ／Ｄ変換する際に、本来の信号の周波数とは違う偽物の信号が発生してしまう現象のことである。エリアジングが発生するのは、Ａ／Ｄ変換器のサンプリング周波数の１／２を越える周波数の信号がＡ／Ｄ変換器に入力された場合である。例えば、４ＭＨｚのＡ／Ｄ変換器であれば２ＭＨｚ以上の信号を入力するとエリアジングが発生する。このため、本第１実施形態では、サンプリング周波数の１／２を超える周波数成分をカットするアンチエリアジングフィルタ（ローパスフィルタ）をＡ／Ｄ変換器１３の前に設けている。
Ａ／Ｄ変換器１３には直交変換器１４が接続されている。該直交変換器１４は直交変換を行う部分であってヒルベルト変換フィルタで構成されており、入力された信号の位相を９０度遅らせる働きを持つ。ビルベルト変換フィルタは所定のフィルタ係数を持つデジタルＦＩＲフィルタで実現することができる。
遅延器１５は、直交変換器１４にて発生するデータ遅延量と等価な遅延を与える働きを持つ。直交変換器１４及び遅延器１５の出力データは互いに９０度の位相差を持つペア信号となる。
ＭＴＩフィルタ１６，１７は、それぞれ直交変換器１４及び遅延器１５からのデータに対して、血流以外の組織からのゆっくりとしたドプラ信号を除去するためのハイパスフィルタであり、デジタルＩＩＲフィルタで構成する。フィルタ係数は、図示しない外部入力部から任意に設定可能とする。
第１及び第２の自己相関演算部１８，１９は、本発明の特徴であるデータ処理を行う部分であり、ＭＴＩフィルタ１６，１７の出力に対してそれぞれ異なる遅延時間の自己相関の振幅値を得ことで、血流情報データを演算する演算処理を行い、乗算器２０に出力する。
乗算器２０は、第１及び第２の自己相関演算部１８，１９がそれぞれ出力する第１の血流情報データと第２の血流情報データを乗算して乗算血流情報データを出力する。
Ｂモード像演算部２１は、直交変換器１４及び遅延器１５の出力データを基にＢモード画像データを生成する部分である。直交変換器１４及び遅延器１５の出力データは互いに直交する関係にあり、両者の２乗和の平方根を計算することで、反射エコーの振幅情報を得ることができる。得られた振幅情報は、対数圧縮、ゲイン調整、コントラスト調整などの処理を行いＢモード画像データとなる。
ＤＳＣ２２は、乗算器２０で演算された乗算血流情報データ及びＢモード像演算部２１で演算されたＢモード画像データを基に超音波探触子１１のスキャン形状にあった形に座標変換及び補間処理し、乗算血流情報データ及びＢモード画像データが合成された超音波画像データを生成する。
ここで、例えば、超音波探触子１１として曲率６ｍｍで画角１８０度のコンベツクス探触子を用いた場合は、放射状の半円の２次元音線データが得られる。２次元音線データは、円の中心からの距離と音線の角度によって決まる極座標（ｒ，θ）形式のデータである。この場合、ＤＳＣ２２は、座標変換の処理として、極座標のデータをモニター表示に適した直交座標（ｘ，ｙ座標）に変換する。具体的には、以下の関係式（１）、（２）を利用して極座標（ｒ，０）と直交座標（ｘ，ｙ）の変換を行う。
ｘ＝ｒ×ｃｏｓθ … （１）
ｙ＝ｒ×ｓｉｎθ … （２）
その際、音腺データの密度が荒いとデータとデータの間に隙間ができてしまうため、ＤＳＣ２２は、周囲のデータから補間処理によってデータを作り出す。例えば、ＤＳＣ２２は、周囲の４点のデータを利用して補間する４点補間を行う。
モニタ２３は、ＤＳＣ２２で生成された超音波画像データを映像表示する。
（作用）
以下、第１実施形態の超音波診断装置の作用を図２〜図５を参照して説明する。なお、この作用は当該超音波診断装置の観察動作として説明する。
まず、図２に示すステップＳ１において、超音波探触子１１と超音波送受信部１２は、所定の時間間隔で同一の音線方向に対して複数回の超音波送受信を行って複数の超音波データを得る。これにより、各音線では、所定の同一深さのデータが送受信の回数分得られる。
次に、超音波送受信部１２は、図２に示すステップＳ２において、ステップＳ１で得られた複数の超音波データの中から所定の同一の深さに対応するデータを抽出する。
このようなステップＳ１，Ｓ２の処理において、例えば、同一の音線に対して６回送受信を繰り返した場合は、所定の同一深さのデータが６個得られる。
次に、図２に示すステップＳ３において、超音波診断装置１は、ＭＴＩフィルタ１６，１７に、直交変換器１４及び遅延器１５を介したデータを出力する。これにより、ＭＴＩフィルタ１６，１７は、エコーデータを出力する。
以下、図２に示すステップＳ４における第１の自己相関演算部１８及び図２に示すステップＳ５における第２の自己相関演算部１９内部での演算について説明する。
前述したように、各音線では、所定の同一深さのデータが送受信の回数分得られる。例えば、同一の音線に対して６回送受信を繰り返した場合は、所定の同一深さのデータが図３Ｂに示す元データＺ（ｉ）のように６個得られている。データａ〜データｆは、ＭＴＩフィルタ１６からのデータを虚数成分、ＭＴＩフイルタ１７を実数成分として組み合わせた複素データである。一般に同一の音線に対してＮ回送受信を繰り返した場合の元データＺ（ｉ）は以下の式（３）で表現される。
Ｚ（ｉ）＝Ｘ（ｉ）＋ｊＹ（ｉ） …（３）
但し、ｉ＝１〜Ｎの整数である。
第１の自己相関演算部１８では、前記元データＺ（ｉ）を図３Ａに示すクロック信号ＣＬＫのＭクロック分だけ遅延させ、この遅延させたデータを図３Ｃに示す遅延データＺ_Ｍ（ｉ）とする。Ｍは１からＮ−１までの任意の整数であるが、Ｍを大きくするとベクトル偏移データをとれるデータ数が減るのでＭ＝１〜２程度が望ましい。
このとき、図３Ｄのように、各時刻においてベクトル偏移データＲ_Ｍ（ｉ）を計算することができる。Ｒ_Ｍ（ｉ）は一般的に式（４）で表現される。
Ｒ_Ｍ（ｉ）＝Ｚ_Ｍ（ｉ）＊×Ｚ（ｉ） …（４）
但し、＊は共役複素数を表す。
図３Ｅにおいて、ゲートパルスはベクトル偏移データＲ_Ｍ（ｉ）から有効データが出力される期間と等価な時間幅を持つパルスである。
ベクトル偏移データＲ_Ｍ（ｉ）は、第１の自己相関演算部１８内の図示しない積分器においてゲートパルスがＨｉｇｈの期間中だけ積分され、第１の自己相関値が演算される。
図３Ｅに示すゲートパルスは、図３Ｂに示す元データと図３Ｃに示す遅延データが重なり合う時間の間だけＨｉｇｈになるので、その期間中だけ図示しない積分器においてベクトル偏移データＲ_Ｍ（ｉ）が積分される。次にゲートパノレスがＨｉｇｈになる前に積分値は０にリセットされる。
図３Ａ〜図３Ｅに示すデータの例で考えると、図３Ｄに示す第１の自己相関データに対応する第１の自己相関値Ｒ_Ｍは以下の式（５）で表現される。
Ｒ_Ｍ＝ｃａ＊＋ｄｂ＊＋ｅｃ＊＋ｆｄ＊ …（５）
最終的には、複素数である第１の自己相関値Ｒ_Ｍの振幅値をとったものを第１の血流情報データとする。
同様に、第２の自己相関演算部１９では、図４Ｂに示す元データＺ（ｉ）を図４Ａに示すクロック信号ＣＬＫのＫクロック分だけ遅延させ、この遅延させたデータを図４Ｃに示す遅延データＺ_Ｋ（ｉ）とする。Ｋは１からＮ−１までの任意の整数であり、Ｍと同じ値であっても良い。Ｋの値もＫ＝１〜２が望ましい。
このとき、図４Ｄのように、ベクトル偏移データＲ_Ｋ（ｉ）は一般的に式（６）で表現される。
Ｒ_Ｋ（ｉ）＝Ｚ_Ｋ（ｉ）＊×Ｚ（ｉ） …（６）
図４Ｅにおいて、ゲートパルスはベクトル偏移データＲ_Ｋ（ｉ）から有効データが出力される期間と等価な時間幅を持つパルスである。
第２の自己相関演算部１９でも第１の自己相関演算部１８と同様に自己相関値が計算され、図４Ａ〜図４Ｅに示すデータの例で考えると、第２の自己相関演算部１９における図４Ｄに示す第２の自己相関データに対応する第２の自己相関値Ｒ_Ｋは式（７）で表現される。
Ｒ_Ｋ＝ｂａ＊＋ｃｂ＊＋ｄｃ＊＋ｅｄ＊＋ｆｅ＊ …（７）
最終的には、複素数である第２の自己相関値Ｒ_Ｋの振幅値をとったものを第２の血流情報データとする。
第１の血流情報データと第２の血流情報データに共通して、ホワイトノイズは各ベクトル偏移データがランダムな位相を持つため、それらを加算した自己相関の振幅値は０に近づく。
一方、血流からのエコーは各ベクトル偏移データがほぼ同じ位相を持つため、それらを加算した自己相関の振幅値は特定の位相の向きに加算され、振幅はホワイトノイズの場合と比べて大きくなる。生体組織からのエコーは各ベクトル偏移データがほぼ同じ位相を持つうえに各ベクトル偏移データの振幅が大きいため、自己相関の振幅値も大きくなる。
従って、第１の血流情報データ、第２の血流情報データともに、図５に示すように値が小さい領域はホワイトノイズを表し、値が大きい領域は生体組織からのエコーを表し、その中間の値が血流からのエコーを表すことになる．
第１の血流情報データと第２の血流情報データは、図２に示すステップＳ６において、乗算器２０で乗算され、乗算血流情報データが生成される。乗算器２０で乗算を行うことで、乗算血流情報データのうち、ホワイトノイズがとる値はより小さくなり、生体組織からのエコーがとる値はより大きくなるため、結果として、ホワイトノイズ、血流からのエコー、生体組織からのエコーそれぞれの比がより大きくなる。従って、ホワイトノイズ、血流からのエコー、生体組織からのエコーそれぞれの弁別がより容易となる。
そのため、ＤＳＣ２２が特定の大きさの乗算血流情報データを抽出してモニタ２３に画面表示することで、血流からのエコーのみを表示することが可能となる。例えば、ホワイトノイズと血流エコーとの境界を定める閾値Ｂ１、及び生体組織からのエコーと血流からのエコーの境界を定める閾値Ｂ２を設け、ＤＳＣ２２が閾値Ｂ１以上でかつ閾値Ｂ２以下の値を持つデータのみをモニタ２３に画面表示することで、血流からのエコーのみを表示できる。
閾値Ｂ１，Ｂ２は、超音波送受信部１２の送信電圧と超音波送受信部１２内の受信回路の増幅率により一意的に決まる定数である。閾値Ｂ１，Ｂ２は、実際に使用する送信電圧と受信回路を決めてから、実験的に求めて、一番効果的な値に設定する。
ＤＳＣ２２は、乗算血流情報データに対して、振幅値を所定の色相の輝度に対応させて表示、もしくは振幅値を所定の色相に対応させてカラー表示を行う。
このような動作により、図２に示すステップＳ１は、所定の時間間隔で同一の音線方向に対して複数回の超音波送受信を行って複数の超音波データを得る超音波データ取得ステップとなっている。
図２に示すステップＳ２は、前記ステップＳ１で得られた複数の超音波データの中から所定の同一の深さに対応するデータを得るデータ抽出ステップとなっている。
図２に示すステップＳ３は、前記ステップＳ２で得られたで得られた所定の同一の深さに対応するデータを所定の特性を有するフィルタ手段に入力し、体内の運動部位からのエコーデータを抽出するフィルタリングステップとなっている。
図２に示すステップＳ４は、前記ステップＳ３で得られた前記体内の運動部位からのエコーデータを、第１の所定時間遅延させた場合の第１の自己相関の振幅値を得る第１の振幅値取得ステップと、
図２に示すステップＳ５は、前記ステップＳ３で得られた前記体内の運動部位からのエコーデータを、第２の所定時間遅延させた場合の第２の自己相関の振幅値を得る第２の振幅値取得ステップとなっている。
図２に示すステップＳ６は、前記第１及び第２の自己相関の振幅値を掛け合わせて得られた乗算値を所定の色相の輝度に対応させるカラー化ステップとなっている。
図２に示した超音波診断装置のデータ処理方法について、図６を用いて具体的な処理として説明する。
まず、図６に示すステップＳ１１において、超音波探触子１１と超音波送受信部１２は、複数の超音波データを取得する。
次に、図６に示すステップＳ１２において、直交変換器１４により複数の超音波データの直交変換を行う。
次に、図６に示すステップＳ１３において、超音波送受信部１２は、図６に示すステップＳ１１で得られた複数の超音波データの中から所定の同一の深さに対応するデータを抽出する。
このようなステップＳ１１，Ｓ１２，Ｓ１３の処理により、所定の同一深さのデータが図３Ａ〜図３Ｅ及び図４Ａ〜図４Ｅに示すように６個得られる。
次に、図６に示すステップＳ１４において、ＭＴＩフィルタ１６，１７は、直交変換器１４及び遅延器１５を介した所定の同一深さのデータに対してＭＴＩフィルタ演算を行いエコーデータを出力する。
次に、図６に示すステップＳ１５において、第１の自己相関演算部１８は、ＭＴＩフィルタ１６，１７のエコーデータに対して第１の血流情報データ演算を行う。
次に、図６に示すステップＳ１６において、第２の自己相関演算部１９は、ＭＴＩフィルタ１６，１７のエコーデータに対して第２の血流情報データ演算を行う。
次に、図６に示すステップＳ１７において、乗算器２０は、ステップＳ１５，Ｓ１６で生成した第１の血流情報データと第２の血流情報データを乗算し、乗算血流情報データを生成する。
次に、図６に示すステップＳ１８において、ＤＳＣ２２は、乗算器２０からの乗算血流情報データに対して閾値Ｂ１以上でかつ閾値Ｂ２以下の値を持つ乗算血流情報データを抽出する。この後、図６に示すステップＳ１９において、ＤＳＣ２２は、すべての計算対象データの演算が終了したかの判別を行い、演算が終了していない場合にはステップＳ１３の処理に戻り、演算が終了した場合にはステップＳ２０の処理に移行する。
図６に示すステップＳ２０において、ＤＳＣ２２は、ステップＳ１８で抽出された乗算血流情報データに色相と輝度の割り当てを行い、図６に示すステップＳ２１において、ＤＳＣ２２は、Ｂモード像演算部２１で演算されたＢモード画像データと前記乗算血流情報データが合成された超音波画像データを生成して、モニタ２３に表示する。この後、図６に示すステップＳ１１の処理に戻る。
以下、図７を用いてステップＳ１５の処理について説明する。
図７に示すステップＳ３１〜Ｓ３４は図６に示すステップＳ１５を構成している。
まず、図７に示すステップＳ３１において、第１の自己相関演算部１８は、図３Ｂに示す元データＺ（ｉ）を図３Ａに示すクロック信号ＣＬＫのＭクロック分だけ遅延させ図３Ｃに示す遅延データＺ_Ｍ（ｉ）を生成する。
次に、図７に示すステップＳ３２において、第１の自己相関演算部１８は、前記式（４）に示した演算を行い、ベクトル偏移データＲ_Ｍ（ｉ）を生成する。
次に、図７に示すステップＳ３３において、第１の自己相関演算部１８は、ベクトル偏移データＲ_Ｍ（ｉ）の積分を行い第１の自己相関値Ｒ_Ｍを生成する。
次に、図７に示すステップＳ３４において、第１の自己相関演算部１８は、第１の自己相関値Ｒ_Ｍの振幅値を演算すること第１の血流情報データで生成する。
以下、図８を用いてステップＳ１６の処理について説明する。
図８に示すステップＳ４１〜Ｓ４４は図６に示すステップＳ１６を構成している。
まず、図８に示すステップＳ４１において、第２の自己相関演算部１９は、図４Ｂに示す元データＺ（ｉ）を図４Ａに示すクロック信号ＣＬＫのＫクロック分だけ遅延させ図４Ｃに示す遅延データＺ_Ｋ（ｉ）を生成する。
次に、図８に示すステップＳ４２において、第２の自己相関演算部１９は、前記式（６）に示した演算を行い、ベクトル偏移データＲ_Ｋ（ｉ）を生成する。
次に、図８に示すステップＳ４３において、第２の自己相関演算部１９は、ベクトル偏移データＲ_Ｋ（ｉ）の積分を行い第２の自己相関値Ｒ_Ｋを生成する。
次に、図８に示すステップＳ４４において、第２の自己相関演算部１９は、第２の自己相関値Ｒ_Ｋの振幅値を演算すること第２の血流情報データで生成する。
（効果）
血流情報データが単純なパワー値ではホワイトノイズもある程度の値を持つのに対して、乗算器２０が出力する乗算血流情報データにおけるホワイトノイズの値は０に近づく。そのため、従来の血流のパワーを表示する方法と比べて第１実施形態における超音波診断装置のデータ処理方法では血流信号とホワイトノイズの比が大きくなる。すなわち、ホワイトノイズを抑制する効果に優れている。
これにより、第１実施形態の超音波診断装置は、血流可視化を行うとともにノイズの抑制効果を向上することができ、血流情報をより正確に検出できる。
（第２の実施形態）
図９〜図１１は本発明の第２実施形態の超音波診断装置を示したものである。
なお、図９〜図１１において、図１〜図８に示した上記第１実施形態の超音波診断装置と同様の構成要素には同じ符号を付するものとし、ここでの詳しい説明は省略する。
（構成）
図９に示すように、本発明の第２実施形態の超音波診断装置２では、上述した第１実施形態の超音波診断装置における乗算器２０の代わりに平均処理器３０を設けている。
この平均処理器３０は第１及び第２の自己相関演算部１８，１９がそれぞれ出力する第１の血流情報データと第２の血流情報データを振幅値を平均して得られた平均血流情報データを出力する。
ＤＳＣ２２は、平均処理器３０で演算された平均血流情報データ及びＢモード像演算部２１で演算されたＢモード画像データを基に超音波探触子１１のスキャン形状にあった形に座標変換及び補間処理し、平均血流情報データ及びＢモード画像データが合成された超音波画像データを生成する。
（作用）
図１０に示すように、第２実施形態の超音波診断装置２では、第１実施形態におけるステップＳ６の処理の代わりに、ステップＳ５６に示す前記第１及び第２の自己相関の振幅値を平均して得られた平均値を所定の色相の輝度に対応させるカラー化ステップを用いている。
また、図１１に示すように、第２実施形態の超音波診断装置２では、第１実施形態におけるステップＳ１７の処理の代わりに、ステップＳ６７に示す処理を用いている。
ステップＳ６７において、平均処理器３０は、ステップＳ１５，Ｓ１６で生成した第１の血流情報データと第２の血流情報データの平均を演算することで平均血流情報データを出力する。
上記第１実施形態の超音波診断装置について説明したように、第１の自己相関演算部１８が出力する第１の血流情報データ及び第２の自己相関演算部１９が出力する第２の血流情報データにおいて、ホワイトノイズ、血流からのエコー、生体組織からのエコーのそれぞれは図５に示すような大きさの関係にある。
従って、第１の血流情報データ及び第２の血流情報データを平均処理器３０で平均した平均血流情報データにおいても、血流からのエコー、生体組織からのエコーのそれぞれは図５とほぼ同じ大きさの関係にある。ただし、ホワイトノイズに関しては、平均をとることで突発的な大きな値を抑制する効果が得られるため、より低い値を安定してとるようになる。
従って、ステップＳ５６（図１１の場合ステップＳ６７）において特定の大きさの平均血流情報データを抽出して表示することで、血流からのエコーのみを表示することが可能となる。例えば、ホワイトノイズと血流エコーとの境界を定める閾値Ｂ１、及び生体組織からのエコーと血流からのエコーの境界を定める閾値Ｂ２を設け、閾値Ｂ１以上でかつ閾値Ｂ２以下の値を持つデータのみを画面表示することで、血流からのエコーのみを表示できる。
ＤＳＣ２２は、平均血流情報データを、振幅値を所定の色相の輝度に対応させて表示、もしくは振幅値を所定の色相に対応させてカラー表示を行う。
（効果）
血流情報データが単純なパワー値ではホワイトノイズもある程度の値を持つのに対して、平均処理器３０からの平均血流情報データにおけるホワイトノイズの値は０に近づく。これにより、第２実施形態の超音波診断装置は、第１実施形態の超音波診断装置と同様に、血流可視化を行うとともにノイズの抑制効果を向上することができ、血流情報をより正確に検出できる。
また、第２実施形態の超音波診断装置では、ホワイトノイズに関して突発的な大きな値を抑制する効果が得られるため、血流情報をより安定に検出できる。
（第３の実施形態）
図１２〜図１３は本発明の第３実施形態の超音波診断装置を示したものである。
なお、図１２〜図１３において、図１〜図８に示した上記第１実施形態および図９〜図１１に示した第２実施形態の超音波診断装置と同様の構成要素には同じ符号を付するものとし、ここでの詳しい説明は省略する。
（構成）
図１２に示すように、本発明の第３実施形態の超音波診断装置３は、上述した第２実施形態の超音波診断装置に対して、前記平均処理器３０の後段にさらに血流像演算部４０が配設されている。
この血流像演算部４０は、第１の自己相関演算部１８において演算された第１の自己相関値Ｒ_Ｍ、及び平均処理器３０で演算された平均血流情報データを入力し、所定の演算を経て血流像データを出力する。
この血流像データは、図示しない操作パネルから操作者が選択した動作モードに応じて血流速度データ、もしくは処理された平均血流情報データの何れかとなる。
また、本実施形態において、デジタルスキャンコンバータ２２は血流像演算部４０で演算された血流像データ及びＢモード像演算部２１で演算されたＢモード画像データを基に超音波探触子１１のスキャン形状にあった形に座標変換および補間処理を行い、血流像データ及びＢモード画像データが合成された超音波画像データを生成する。
（作用）
図１３に示すように、本第３実施形態の超音波診断装置３は、上記第２実施形態におけるステップＳ１８の代わりにステップＳ１０１に示すしきい値処理ステップを有し、さらに、ステップＳ１００に示す速度値の演算ステップを有する。
まず、ステップＳ１００について説明する。
血流の速度ｖは血流像演算部４０において以下の式により演算される。
ｖ＝｛ｃ／（４πｆ_０ＴＭ）｝ｔａｎ^−１（Ｒ_Ｍｙ／Ｒ_Ｍｘ）
ｃ：音速
ｆ_０：送受信した超音波パルスの中心周波数
Ｔ：同一の音線に送受信の繰り返しを行う際の時間周期
Ｍ：上記第１実施形態で述べた遅延を定める整数
Ｒ_Ｍｘ：自己相関値Ｒ_Ｍの実数成分
Ｒ_Ｍｙ：自己相関値Ｒ_Ｍの虚数成分
次にステップＳ１０１について説明する。
はじめに血流像演算部４０は、以下のしきい値処理を行って血流速度データを演算する。すなわち、平均血流情報データをしきい値Ｂ１、Ｂ２と比較し、平均血流情報データがしきい値Ｂ１以上でかつしきい値Ｂ２以下の値を持つ場合は、血流からのエコーであると判断できるため、血流速度データは速度ｖの値そのものとする。
平均血流情報データがしきい値Ｂ１より小さい場合、もしくはしきい値Ｂ２より大きい場合には、血流からのエコーではないと判断し、血流速度データはゼロに置き換える。一般に、血流速度データがゼロの場合はカラー表示を行わないため、血流からのエコーの速度値のみをカラー表示することができる。
また、血流像演算部４０は、処理された平均血流情報データを以下のように求める。
血流の速度ｖの絶対値を所定のしきい値ｖ_０と比較し、血流の速度ｖの絶対値がしきい値ｖ_０よりも大きい場合は、血流からのエコーであると判断し、処理された平均血流情報データは平均血流情報データそのものとする。血流の速度ｖの絶対値がしきい値ｖ_０よりも小さい場合は、生体組織からのエコーであると判断し、処理された平均血流情報データとしてゼロを出力する。
一般に、処理された平均血流情報データがゼロの場合はカラー表示を行わないため、血流からのエコーの平均血流情報データのみをカラー表示することができる。
血流像演算部４０は、図示しない操作パネルから操作者が選択した動作モードに応じて血流速度データもしくは処理された平均血流情報データの何れかを血流像データとして出力する。
なお、ステップＳ１００における速度値の演算ステップは第１の自己相関演算部１８において演算された第１の自己相関値Ｒ_Ｍを用いる場合について説明したが、第１の自己相関演算部１８の代わりに第２の自己相関演算部１９で演算された第２の自己相関値Ｒ_Ｋを用いてもよい。
また、上述した第１の実施形態および第２の実施形態における説明から、本第３の実施形態において平均処理器３０の代わりに乗算器２０を用いる構成でも良いことは明かである。
さらに、平均処理器３０での処理は必ずしも単純な平均処理に限定されず、所定の重みを掛けた平均もしくは加算であっても良い。
さらに、本発明におけるＭＴＩフィルタはデジタルＩＩＲフィルタに限定されず、デジタルＦＩＲで構成しても良い。
（効果）
血流像データをカラー表示する場合に、平均処理器３０から出力される平均血流情報データ及び血流像演算部４０で演算される速度値を用いたしきい値処理を行うことで、ホワイトノイズおよび生体組織からのエコーのカラー表示を抑制することができ、血流の速度情報をより正確に検出できる超音波診断装置を実現することができる。
なお、上述した各実施の形態等を部分的等で組み合わせる等して構成される実施の形態等も本発明に属する。
また、本発明は上述した実施形態に限定されるものではなく、本発明の主旨を逸脱しない範囲内において種々の変形や応用が可能であることは勿論である。

以上説明したように本発明によれば、血流可視化を行うとともにノイズの抑制効果を向上できるので、血流情報をより正確に検出できる超音波診断装置を実現できる。
関連出願へのクロスリファレンス
本出願は、２００３年５月２９日に日本国に出願された特願２００３−１５２９５８号を優先権主張の基礎として出願するものであり、上記の開示内容は、本願明細書、請求の範囲、図面に引用されたものとする。

Claims

超音波診断装置のデータ処理方法は以下を含む；
所定の時間間隔で同一の音線方向に対する複数回の超音波送受信を行って複数の超音波データを得る超音波送受信手段により得られた複数の超音波データの中から所定の同一の深さに対応するデータを得るデータ抽出ステップ、
前記データ抽出ステップで得られた所定の同一の深さに対応するデータを所定の特性を有するフィルタに入力し、体内の運動部位からのエコーデータを抽出するフィルタリングステップ、
前記フィルタリングステップで得られた前記体内の運動部位からのエコーデータを、第１の所定時間遅延させた場合の第１の自己相関の振幅値を得る第１の振幅値取得ステップ、
前記フィルタリングステップで得られた前記体内の運動部位からのエコーデータを、第２の所定時間遅延させた場合の第２の自己相関の振幅値を得る第２の振幅値取得ステップ、
前記第１、第２の自己相関の振幅値を掛け合わせて得られた乗算値を所定の色相の輝度に対応させるカラー化ステップ。
請求項１に記載の超音波診断装置のデータ処理方法において、
前記データ抽出ステップにおける超音波送受信は超音波探触子を用いて行う。
請求項１に記載の超音波診断装置のデータ処理方法において、
前記所定のフィルタは、ＭＴＩフィルタである。
超音波診断装置のデータ処理方法は以下を含む；
所定の時間間隔で同一の音線方向に対する複数回の超音波送受信を行って複数の超音波データを得る超音波送受信手段により得られた複数の超音波データの中から所定の同一の深さに対応するデータを得るデータ抽出ステップ、
前記データ抽出ステップで得られた所定の同一の深さに対応するデータを所定の特性を有するフィルタに入力し、体内の運動部位からのエコーデータを抽出するフィルタリングステップ、
前記フィルタリングステップで得られた前記体内の運動部位からのエコーデータを、第１の所定時間遅延させた場合の第１の自己相関の振幅値を得る第１の振幅値取得ステップ、
前記フィルタリングステップで得られた前記体内の運動部位からのエコーデータを、第２の所定時間遅延させた場合の第２の自己相関の振幅値を得る第２の振幅値取得ステップ、
前記第１、第２の自己相関の振幅値を平均して得られた平均値を所定の色相の輝度に対応させるカラー化ステップ。
請求項４に記載の超音波診断装置のデータ処理方法において、
前記データ抽出ステップにおける超音波送受信は超音波探触子を用いて行う。
請求項４に記載の超音波診断装置のデータ処理方法において、
前記所定のフィルタは、ＭＴＩフィルタである。
超音波診断装置は、以下を含む；
超音波振動子により所定の時間間隔で同一の音線方向に対する複数回の超音波送受信を行って複数の超音波データを得る超音波データ収得手段、
前記超音波データ取得手段により得られた複数の超音波データの中から所定の同一の深さに対応するデータを得るデータ抽出手段、
所定のフィルタ特性を有し、前記データ抽出手段により得られた所定の同一の深さに対応するデータを入力し、体内の運動部位からのエコーデータを抽出するフィルタ、
前記フィルタにより得られた前記体内の運動部位からのエコーデータを、第１の所定時間遅延させた場合の第１の自己相関の振幅値を得る第１の振幅値取得手段、
前記フィルタにより得られた前記体内の運動部位からのエコーデータを、第２の所定時間遅延させた場合の第２の自己相関の振幅値を得る第２の振幅値取得手段、
前記第１、第２の自己相関の振幅値を掛け合わせて得られた乗算値を所定の色相の輝度に対応させるカラー化手段。
請求項７に記載の超音波診断装置において、
前記フィルタは、ＭＴＩフィルタである。
超音波診断装置のデータ処理方法は以下を含む；
所定の時間間隔で同一の音線方向に対する複数回の超音波送受信を行って複数の超音波データを得る超音波送受信手段により得られた複数の超音波データの中から所定の同一の深さに対応するデータを得るデータ抽出ステップ、
前記データ抽出ステップで得られた所定の同一の深さに対応するデータを所定の特性を有するフィルタに入力し、体内の運動部位からのエコーデータを抽出するフィルタリングステップ、
前記フィルタリングステップで得られた前記体内の運動部位からのエコーデータを、第１の所定時間遅延させた場合の第１の自己相関の振幅値を得る第１の振幅値取得ステップ、
前記フィルタリングステップで得られた前記体内の運動部位からのエコーデータを、第２の所定時間遅延させた場合の第２の自己相関の振幅値を得る第２の振幅値取得ステップ、
前記第１、第２の自己相関の実数および虚数部から速度値を得るステップ、
前記第１、第２の自己相関の振幅値を掛け合わせて得られた平均値を用いて前記速度値をしきい値処理するしきい値処理ステップ、
前記しきい値処理ステップでしきい値処理された前記速度値を所定の色相の輝度に対応させるカラー化ステップ。