JP3745672B2

JP3745672B2 - 生体信号計測装置及び超音波診断装置

Info

Publication number: JP3745672B2
Application number: JP2001326783A
Authority: JP
Inventors: 孝岡田; 烈光原田
Original assignee: Aloka Co Ltd
Current assignee: Hitachi Ltd
Priority date: 2001-10-24
Filing date: 2001-10-24
Publication date: 2006-02-15
Anticipated expiration: 2021-10-24
Also published as: JP2003126090A

Description

【０００１】
【発明の属する技術分野】
本発明は、心拍に同期した生体信号を計測する装置に関し、特に血管壁の変位や血流速度を計測する超音波診断装置に関する。
【０００２】
【従来の技術】
血管の性状や心臓などの機能を診断するために超音波診断装置が利用される。超音波診断装置を利用して血管壁の変位を計測する場合、その変位信号には呼吸に起因した変動成分が含まれることがある。この呼吸性の変動成分を低減する方法として、血管壁の変位の時間変化を表す血管壁変位波形を心拍周期ごとのサンプル波形に分割し、複数のサンプル波形をアンサンブル平均して、平均波形を生成する方法がある。心拍周期は心電信号に一心拍ごとに現れる特徴的な波形を基準にして把握される。
【０００３】
特徴的な波形としては、一般に心電信号のＲ波が用いられる。従来は、計測者が血管壁変位信号において、心拍間の変動が小さい心拍区間を判断し、例えば、Ｒ波にカーソルを移動させ手動でアンサンブル平均に用いる複数のサンプル波形を取り出していた。
【０００４】
【発明が解決しようとする課題】
しかし、計測者の判断及び手動操作に基づいて、複数のサンプル波形を選択する作業は煩雑で時間を要するという問題があった。また、サンプル波形の選択が計測者の主観に委ねられ、計測結果のばらつきが大きくなるという問題があった。
【０００５】
本発明は上記問題点を解決するためになされたもので、自動的にサンプル波形を短時間で選択し、精度の高いアンサンブル平均波形が得られる生体信号装置及び超音波診断装置を提供することを目的とする。
【０００６】
【課題を解決するための手段】
本発明に係る生体信号計測装置は、心拍に同期した参照生体信号を計測する参照生体信号計測手段と、前記参照生体信号に現れる所定の特徴波形を検出する特徴波形検出手段と、生体から目的生体信号を計測する目的生体信号計測手段と、前記特徴波形に基づいて目的生体信号波形を心拍周期単位のサンプル波形に分割する波形分割手段と、前記サンプル波形相互の相関値を計算する相関値演算手段と、前記サンプル波形についての所定条件に基づいて、複数の前記サンプル波形から複数の代表サンプル波形を選択する代表選択手段と、前記複数の代表サンプル波形をアンサンブル平均して目的生体信号の平均波形を求めるアンサンブル平均手段と、を有し、前記相関値演算手段は、複数のサンプル波形のうちの互いに異なる二つのサンプル波形の全ての組み合わせについて、二つのサンプル波形相互の相関値を計算し、前記代表選択手段は、複数のサンプル波形のうちの各サンプル波形ごとに、そのサンプル波形が関与する複数の相関値のうち値の大きい所定個数の平均値を計算し、複数のサンプル波形についての複数の前記平均値のうちから最大値を検索し、その最大値を与える一群のサンプル波形を前記代表サンプル波形として選択する、ことを特徴とする。
【０００７】
本発明によれば、心拍に同期した波形を有する参照生体信号に基づいて、本来、目的とする生体信号である目的生体信号が心拍周期単位のサンプル波形に区分される。この区分は、特徴波形検出手段が、参照生体信号に心拍に同期して現れる特徴的な波形を検出し、その特徴波形を基準として自動的に行われる。サンプル波形をアンサンブル平均することにより、心拍周波数に比べて低周波の呼吸性変動成分などを低減・除去することができる。一方、アンサンブル平均の対象に他のサンプル波形と著しく異なる特異なサンプル波形が含まれると、平均波形にもその特異性が残存することとなり好ましくない。そこで、代表選択手段は、サンプル波形のうち、アンサンブル平均により平均波形を求める上で好適なものを、当該波形についての所定条件に基づいて選択する。そして、アンサンブル平均手段は、代表選択手段により選択された複数の代表サンプル波形をアンサンブル平均する。アンサンブル平均処理では、各代表サンプル波形の互いに対応する時相のデータの平均値が計算され、各時相での平均値から平均波形が形成される。
【０００８】
他の本発明に係る生体信号計測装置は、前記サンプル波形相互の相関値を計算する相関値演算手段を有し、前記代表選択手段は、前記相関値に基づいて、互いに類似した一群の前記サンプル波形を前記代表サンプル波形として選択する。
【０００９】
一過性のノイズ成分は、それを含むサンプル波形と含まないサンプル波形との相関を低くする。逆に、互いに相関を有する一群のサンプル波形は、そのようなノイズ成分が少ないことを期待することができる。そこで、本発明では、サンプル波形相互の相関値を計算し、その相関値に基づいて互いに類似した一群のサンプル波形を代表サンプル波形として選択する。
【００１０】
本発明の好適な態様は、前記参照信号が、心電信号である生体信号計測装置である。特に、前記特徴波形として、前記心電信号のＲ波を用いることができる。
【００１１】
本発明に係る超音波診断装置は、超音波パルスを送波し、エコー信号を受波する送受波手段と、前記エコー信号に基づいて血管壁の変位を計測する血管壁変位計測手段と、前記エコー信号に基づいて血流速度を計測する血流速度計測手段と、心拍に同期した参照生体信号を計測する参照生体信号計測手段と、前記参照生体信号に現れる所定の特徴波形を検出する特徴波形検出手段と、前記特徴波形に基づいて前記血管壁の変位及び前記血流速度それぞれの信号波形を心拍周期単位のサンプル波形に分割する波形分割手段と、前記サンプル波形についての所定条件に基づき、前記血管壁の変位及び前記血流速度それぞれについて、複数の前記サンプル波形から複数の代表サンプル波形を選択する代表選択手段と、前記血管壁の変位及び前記血流速度それぞれについて、前記複数の代表サンプル波形をアンサンブル平均して血管壁変位平均波形及び血流速度平均波形を求めるアンサンブル平均手段と、前記血管壁変位平均波形及び血流速度平均波形から評価値を演算する評価値演算手段とを有する。
【００１２】
本発明によれば、エコー信号に基づいて得られた血管壁の変位信号及び血流速度信号を目的生体信号として、それらの平均波形が生成される。そしてそれら平均波形に基づいて所定の評価値が演算される。その評価値としては、各種のものが想定されるが、その一つとして、循環動態指標であるウエーブインテンシティが挙げられる。
【００１３】
【発明の実施の形態】
次に、本発明の実施形態である超音波診断装置について図面を参照して説明する。本超音波診断装置は、血管の性状や心臓の機能などを評価するための評価値としてウエーブインテンシティを演算する機能を有している。
【００１４】
ウエーブインテンシティは、もともとは、心臓から末梢へ向かう前進脈波と末梢で反射して心臓へ向かう反射脈波のどちらの作用が優勢であるかを判別するための指標として提唱されたものである。具体的には、ウエーブインテンシティＩは、動脈中の局所部位の圧力をＰ、当該局所部位の血流速度をＵとして、それらのΔｔ間における時間変化分ΔＰ及びΔＵから、
Ｉ＝ΔＰ・ΔＵ ………（１）
として定義される。
【００１５】
つまり、ウエーブインテンシティは、一定の時間間隔Δｔにおける圧力Ｐの変化及び血流速度Ｕの変化の積として定義される。Δｔの取り方に依存しない、時間で正規化したウエーブインテンシティは、以下の（２）式で表現される。
【００１６】
ＷＩ＝（ｄＰ／ｄｔ）・（ｄＵ／ｄｔ） ………（２）
上記（２）式から分かるように、ウエーブインテンシティＷＩは、圧力Ｐの時間微分及び血流速度Ｕの時間微分の積として定義されている。
【００１７】
本超音波診断装置は超音波エコートラッキング法により、例えば、頸動脈の血管壁変位を自動的にトラッキングし、これにより血管径の時間変化が計測される。一方、超音波ドプラ法により血流速度の時間変化が計測される。
【００１８】
従来から、血管径の変化波形と血圧の変化波形との間には相似な関係があることが知られている。トラッキング法で測定した最大血管径及び最小血管径を、被検者の上腕に装着されたカフ型血圧計で測定した最大血圧及び最小血圧で較正すると、血管径変化を局所部位の血圧変化とみなすことができる。
【００１９】
そして、以上の血流速度の変化及び血圧の変化から、上記（２）式に従って、ウエーブインテンシティが算出される。
【００２０】
図１は、本超音波診断装置の概略のブロック図である。図１において、プローブ１０は、超音波パルスの送波及びエコーの受波を行う超音波探触子である。このプローブ１０はアレイ振動子を有しており、そのアレイ振動子の電子走査によって超音波ビームが電子的に走査される。その電子走査方式としては例えば電子リニア走査や電子セクタ走査などを挙げることができる。本実施形態に係る超音波診断装置は、上記のように、ウエーブインテンシティを計測する機能を有しており、そのウエーブインテンシティの計測にあたっては、血管１４の中心軸に超音波ビームの走査によって形成される走査面が一致するようにプローブ１０の生体表面１２への当接位置や当接姿勢が手動操作によって調整される。
【００２１】
送信回路１６は、プローブ１０に対して送信信号を供給する回路であり、その送信回路１６の動作は送受信制御部１８によって制御される。受信回路２０は、プローブ１０からの受信信号に対して増幅や整相加算などの処理を行う回路である。この受信回路２０も送受信制御部１８によって制御される。
【００２２】
送受信制御部１８は、送信ビームの形成及び受信ビームの形成を行うための送受信制御を実行している。
【００２３】
断層画像形成部２２は、断層画像すなわちＢモード画像を形成する回路である。形成された断層画像のイメージ情報は表示処理部２４へ出力される。変位演算部２６は、血管壁の位置、具体的にはプローブ１０から見て手前側にある前壁の位置とプローブ１０から見て奥側にある後壁の位置とを演算し、さらに、前壁の位置と後壁の位置とから血管径を演算する機能を有している。
【００２４】
具体的には、変位演算部２６は、後に示す計測ライン上においてユーザー設定されたトラッキングゲートにおいて、血管壁の位置をトラッキングする機能を有する。
【００２５】
血流速度演算部２８は、計測ライン上のエコー信号からカラードプラ法によって血流速度を求め、計測ライン上に設定されたサンプルゲート内の血流速度の平均値を演算する回路である。変位演算部２６で演算された血管径の変化信号１０２及び血流速度演算部２８で演算された血流速度を表す血流速度信号１０４は、表示処理部２４及びアンサンブル平均演算部３４へ出力される。
【００２６】
ちなみに、計測ライン設定器３０は後に示す計測ラインを設定するための手段であり、またトラッキングゲート設定器３２はトラッキングゲートをマニュアルで設定するための手段であり、それらは例えばキーボードやトラックボールなどのポインティングデバイスによって構成される。
【００２７】
心電計測部４２は、被検者の心電信号１０７を計測して、アンサンブル平均演算部３４へ出力する。
【００２８】
アンサンブル平均演算部３４は、血管径の変化信号１０２、血流速度信号１０４及び心電信号１０７をそれぞれ心拍周期単位のサンプル波形に区切り、サンプル波形をアンサンブル平均して、血管径の変化平均波形、血流速度平均波形及び心電平均波形を生成する。これについては後に図２を用いてより詳しく説明する。
【００２９】
評価値演算部３６は、血管径及び血流速度から、評価値としてのウエーブインテンシティを演算する回路である。後に示す図３にはその具体的な構成例が示されている。評価値演算部３６にて演算されたウエーブインテンシティの値は表示処理部２４へ出力される。
【００３０】
本実施形態においては、評価値演算部３６がウエーブインテンシティを計算するために、血圧計４０から出力される血圧信号１０６、具体的には最大血圧及び最小血圧を参照する。例えば、血圧計４０は被検者の上腕表面に巻き付けられたカフによって被検者の血圧を測定するものである。
【００３１】
表示処理部２４は、表示器３８に表示する表示画像を構成する回路である。表示処理部２４は画像合成機能などを有している。
【００３２】
なお、血流速度演算部２８は、従来の超音波ドプラ診断装置に搭載されている直交検波器や自己相関器などで構成するのが望ましい。
【００３３】
図２には、図１に示したアンサンブル平均演算部３４の具体的な構成例が示されている。Ｒ波検出器５０には心電信号１０７が入力され、Ｒ波検出器５０は心室収縮期の開始タイミングに心電信号に現れるＲ波を検出し、そのタイミングでトリガ信号を出力する。Ｒ波は心電信号においてひときわ高く鋭いピークを形成するので、ピーク検出の技法により検出することができる。
【００３４】
心拍分割器５２には、血管径の変化信号１０２、血流速度信号１０４、心電信号１０７及びＲ波検出器５０からのトリガ信号が入力される。心拍分割器５２は、トリガ信号のタイミングで血管径の変化信号１０２、血流速度信号１０４及び心電信号１０７をそれぞれ区切る。これにより、血管径の変化信号１０２、血流速度信号１０４及び心電信号１０７はそれぞれ、Ｒ波から次のＲ波までの心拍期間ごとのサンプル波形に分割される。心拍分割器５２は各心拍期間に指標を付し、この指標は各信号のサンプル波形のデータを識別するために用いられる。血管径の変化信号１０２、血流速度信号１０４及び心電信号１０７それぞれから生成されたサンプル波形１１２，１１４，１１６はアンサンブル平均演算器５４へ出力される。また血管径の変化信号１０２から生成されたサンプル波形１１２は心拍波形相関演算器５６へも出力される。
【００３５】
心拍波形相関演算器５６は、血管径の変化信号１０２のサンプル波形１１２を順次入力され、例えば、ある一定期間内に入力された複数のサンプル波形を１セットとして、そのセットに含まれる２つのサンプル波形の組み合わせごとに相関値を演算する。この相関演算は波形相互の類似の程度を調べる周知の方法を用いて行われる。
【００３６】
代表選択器５８は、心拍波形相関演算器５６で得られたサンプル波形の組み合わせごとの相関値に基づいて、互いに相関が強い所定数のサンプル波形を選択する。
【００３７】
アンサンブル平均演算器５４は、心拍分割器５２から順次入力される血管径の変化信号１０２、血流速度信号１０４及び心電信号１０７それぞれのサンプル波形１１２，１１４，１１６のうち、代表選択器５８で選択されたサンプル波形と同一心拍期間のものを代表サンプル波形として選択する。この選択は上述した指標に基づいて行われる。アンサンブル平均演算器５４は、血管径の変化信号１０２、血流速度信号１０４及び心電信号１０７についてそれぞれ別個に代表サンプル波形をアンサンブル平均する。アンサンブル平均演算器５４からは、アンサンブル平均によって得られた平均波形を有する血管径の変化平均信号、血流速度平均信号、及び心電平均信号が出力される。なお、これら各平均信号は評価値演算部３６にて利用され、また表示器３８に表示することができる。
【００３８】
図３には、図１に示した評価値演算部３６の具体的な構成例が示されている。アンサンブル平均演算部３４から出力される血流速度平均信号１２４は微分器７０に入力され、微分器７０によって血流速度の時間微分が演算される。その微分結果は乗算器７２へ出力される。
【００３９】
一方、スケーリング部７４には、アンサンブル平均演算部３４から出力される血管径の変化平均信号１２２及び血圧計４０から出力される血圧信号１０６が入力される。スケーリング部７４は、血管径の変化平均信号１２２の最大値、最小値をそれぞれ血圧信号１０６から得られる最大血圧、最小血圧に対応付けることにより、血管径の変化平均信号１２２の波形を血圧波形として較正する。すなわち単位変換を実行する。よって、スケーリング部７４から、換算された血圧信号が出力されることになる。微分器７６はその血圧信号に対する時間微分を実行し、その微分結果が乗算器７２へ出力される。
【００４０】
乗算器７２は、微分器７０及び微分器７６から出力される血流速度の微分結果及び血圧の微分結果を乗算することにより、評価値としてのウエーブインテンシティを求める回路である。すなわち、図３に示す構成は上記の（２）式を実行する回路である。
【００４１】
次に、図４に示すアンサンブル平均の処理フロー図を用いて、アンサンブル平均演算部３４での処理手順を説明する。まず、サンプル波形を採取するサンプル採取期間の長さと、アンサンブル平均される心拍数、すなわち代表選択器５８にて代表サンプル波形として選択されるサンプル波形の個数ｍを設定する（Ｓ２００）。アンサンブル平均演算部３４は、設定されたサンプル採取期間内に変位演算部２６、血流速度演算部２８及び心電計測部４２からそれぞれ入力される血管径の変化信号１０２、血流速度信号１０４及び心電信号１０７を保存する（Ｓ２０５）。これらの信号はデジタルデータとしてメモリに保存することができる。
【００４２】
Ｒ波検出器５０は、保存した心電信号１０７のデータを読み出して、Ｒ波のタイミングを検出する（Ｓ２１０）。心拍分割器５２は、血管径の変化信号１０２、血流速度信号１０４及び心電信号１０７の信号データ列を読み出す。信号データ列は、Ｒ波検出器５０により検出されたＲ波のタイミングで分割され、心拍周期単位の信号データ列が生成される。この一心拍分の信号データ列がそれぞれサンプル波形を定義する。心拍分割器５２は、各信号のサンプル波形に、その時系列順に応じた番号を指標として対応付ける。ここでは、サンプル採取期間の初めに位置するサンプル波形から最後のサンプル波形まで順番に、０番〜ｎ番の指標が付与されるとする（Ｓ２１５）。
【００４３】
心拍分割器５２で生成された各信号のサンプル波形のうち、血管径の変化信号１０２のサンプル波形が指標と共に心拍波形相関演算器５６へ渡される。心拍波形相関演算器５６では、血管径の変化信号１０２の２つのサンプル波形（ｉ番目及びｋ番目のサンプル波形）相互の相関値Ｓ_i,kがｉとｋとの全ての組み合わせについて計算される。この処理は、まずｉ，ｋをそれぞれ初期値０として開始される（Ｓ２２０）。そして、ｉを固定して、あるｉ，ｋについての相関値Ｓ_i,kが計算され、メモリ等の記憶手段に格納される（Ｓ２２５）。なお、ここで、ｉとｋとが等しい場合には、相関値の計算及び記録は行われない。
【００４４】
あるｉ，ｋについての処理Ｓ２２５が終わる度に、ｋが１ずつインクリメントされる（Ｓ２３０）。心拍波形相関演算器５６において、あるｉに対しｋ＝ｎについての相関値の計算及び記録が完了すると（Ｓ２３５）、代表選択器５８は当該ｉについてのｎ個の相関値Ｓ_i,k（ｋ＝０〜ｎ）のうち値の大きい上位ｍ個の平均値mean(Ｓ_i)を計算し、その平均値及び、平均されたサンプル波形の指標を保存する（Ｓ２４０）。一方、心拍波形相関演算器５６はｋ＝ｎまで処理が完了すると、ｉを１だけインクリメントして（Ｓ２４５）、処理Ｓ２２５〜Ｓ２３５を繰り返す。ｉ＝ｎまで処理を完了すると（Ｓ２５０）、代表選択器５８は、保存されているmean(Ｓ_i)の最大値を検索し、その最大値を与える一群のサンプル波形を代表サンプル波形と定め、その一群のサンプル波形の指標をメモリから読み出してアンサンブル平均演算器５４へ出力する（Ｓ２５５）。
【００４５】
アンサンブル平均演算器５４は、血管径の変化信号１０２、血流速度信号１０４及び心電信号１０７について心拍分割器５２で生成された一連のサンプル波形のデータ列のうち、代表選択器５８にて選択された指標に対応するサンプル波形を選んで、それらをアンサンブル平均して出力する（Ｓ２６０）。
【００４６】
なお、アンサンブル平均演算部３４によれば、心拍に同期した各種信号から、呼吸性の変動成分やランダムノイズが除去された信号波形を得ることができる。その信号波形は、本装置ではウエーブインテンシティの算出に利用したが、他の評価値の計測等においても用いることもでき、計測値の精度向上を図ることができる。
【００４７】
また、ここでは、心拍波形相関演算器５６が血管径の変化信号のサンプル波形について相関を計算する例を示したが、例えば血流速度信号のサンプル波形を用いて相関を計算するように構成してもよい。
【００４８】
【発明の効果】
本発明の生体信号計測装置及び超音波診断装置によれば、生体信号から、アンサンブル平均に好適な心拍期間が選択され、アンサンブル平均波形が生成される。これにより、心拍に同期した生体信号に対するアンサンブル平均処理が自動化され、省力化が図られると共に、心拍期間単位での相関値等の客観的基準に基づいて選択されたサンプル波形についてアンサンブル平均が行われるので、精度の高い平均波形が得られる。特に、ウエーブインテンシティを求める場合には、生体信号として血管径の変化信号から求めた血圧信号及び血流速度信号に対して微分処理が行われるが、アンサンブル平均波形を用いることでピークの位置の精度が向上し、ウエーブインテンシティの計測精度が向上する。
【図面の簡単な説明】
【図１】本発明に係る超音波診断装置の全体構成を示すブロック図である。
【図２】アンサンブル平均演算部の具体的な構成例を示すブロック図である。
【図３】評価値演算部の具体的な構成例を示すブロック図である。
【図４】アンサンブル平均の処理フロー図である。
【符号の説明】
１０プローブ、１６送信回路、１８送受信制御部、２０受信回路、２２断層画像形成部、２４表示処理部、２６変位演算部、２８血流速度演算部、３０計測ライン設定器、３２トラッキングゲート設定器、３４アンサンブル平均演算部、３６評価値演算部、３８表示器、４０血圧計、４２心電計測部、５０Ｒ波検出器、５２心拍分割器、５４アンサンブル平均演算器、５６心拍波形相関演算器、５８代表選択器、７０，７６微分器、７２乗算器、７４スケーリング部。

Claims

心拍に同期した参照生体信号を計測する参照生体信号計測手段と、
前記参照生体信号に現れる所定の特徴波形を検出する特徴波形検出手段と、
生体から目的生体信号を計測する目的生体信号計測手段と、
前記特徴波形に基づいて目的生体信号波形を心拍周期単位のサンプル波形に分割する波形分割手段と、
前記サンプル波形相互の相関値を計算する相関値演算手段と、
前記サンプル波形についての所定条件に基づいて、複数の前記サンプル波形から複数の代表サンプル波形を選択する代表選択手段と、
前記複数の代表サンプル波形をアンサンブル平均して目的生体信号の平均波形を求めるアンサンブル平均手段と、
を有し、
前記相関値演算手段は、複数のサンプル波形のうちの互いに異なる二つのサンプル波形の全ての組み合わせについて、二つのサンプル波形相互の相関値を計算し、
前記代表選択手段は、複数のサンプル波形のうちの各サンプル波形ごとに、そのサンプル波形が関与する複数の相関値のうち値の大きい所定個数の平均値を計算し、複数のサンプル波形についての複数の前記平均値のうちから最大値を検索し、その最大値を与える一群のサンプル波形を前記代表サンプル波形として選択する、
ことを特徴とする生体信号計測装置。
請求項１記載の生体信号計測装置において、
前記相関値演算手段は、０番からｎ番までの複数のサンプル波形のうちの互いに異なる二つのサンプル波形（ｉ番目及びｋ番目のサンプル波形）相互の相関値Ｓ _i,k を、ｉとｋとの全ての組み合わせについて計算し、
前記代表選択手段は、番号ｉについてのｎ個の相関値Ｓ _i,k のうち値の大きい上位ｍ個の平均値 mean( Ｓ _i ) を計算し、番号ｉが０番からｎ番までの複数の前記平均値 mean( Ｓ _i ) のうちから最大値を検索し、その最大値を与える一群のサンプル波形を前記代表サンプル波形として選択する、
ことを特徴とする生体信号計測装置。
請求項１又は請求項２に記載の生体信号計測装置において、
前記参照生体信号は、心電信号であることを特徴とする生体信号計測装置。
請求項３記載の生体信号計測装置において、
前記特徴波形は、前記心電信号のＲ波であることを特徴とする生体信号計測装置。
超音波パルスを送波し、エコー信号を受波する送受波手段と、
前記エコー信号に基づいて血管壁の変位を計測する血管壁変位計測手段と、
前記エコー信号に基づいて血流速度を計測する血流速度計測手段と、
心拍に同期した参照生体信号を計測する参照生体信号計測手段と、
前記参照生体信号に現れる所定の特徴波形を検出する特徴波形検出手段と、
前記特徴波形に基づいて前記血管壁の変位及び前記血流速度それぞれの信号波形を心拍周期単位のサンプル波形に分割する波形分割手段と、
前記サンプル波形相互の相関値を計算する相関値演算手段と、
前記サンプル波形についての所定条件に基づき、前記血管壁の変位及び前記血流速度それぞれについて、複数の前記サンプル波形から複数の代表サンプル波形を選択する代表選択手段と、
前記血管壁の変位及び前記血流速度それぞれについて、前記複数の代表サンプル波形をアンサンブル平均して血管壁変位平均波形及び血流速度平均波形を求めるアンサンブル平均手段と、
前記血管壁変位平均波形及び血流速度平均波形から評価値を演算する評価値演算手段と、
を有し、
前記相関値演算手段は、複数のサンプル波形のうちの互いに異なる二つのサンプル波形の全ての組み合わせについて、二つのサンプル波形相互の相関値を計算し、
前記代表選択手段は、複数のサンプル波形のうちの各サンプル波形ごとに、そのサンプル波形が関与する複数の相関値のうち値の大きい所定個数の平均値を計算し、複数のサンプル波形についての複数の前記平均値のうちから最大値を検索し、その最大値を与える一群のサンプル波形を前記代表サンプル波形として選択する、
ことを特徴とする超音波診断装置。
請求項５記載の超音波診断装置において、
前記アンサンブル平均手段は、前記血管壁の変位及び前記血流速度それぞれについて、前記複数の代表サンプル波形のうちの各代表サンプル波形の互いに対応する時相のデータの平均値を計算し、各時相での平均値から平均波形を求める、
ことを特徴とする超音波診断装置。