(57)【要約】      （修正有）【課題】計算量が少なく、かつ高速の速度まで算出が可能なドップラ演算法を提供する。【解決手段】受信信号の互いに隣接する２つのゼロクロス点に挟まれた半波形の代表振幅値および代表時刻を各半波形ごとに算出することにより、半波形ごとの代表振幅値および代表時刻からなる圧縮データを生成し、生成された圧縮データを構成する代表振幅値を対応する代表時刻に並べたときの、波形どうしが近似している領域のずれ量に基づいて、被検体内の被検体構成要素の動きを算出する。【選択図】図３
2015084977.tif 000002
【特許請求の範囲】
【請求項１】
被検体内への超音波ビームの送波と該被検体内で反射して戻ってきた超音波の受波とを担うプローブと、
複数のパルス信号を前記プローブに送信して該プローブに超音波ビームを送波させる送信部と、
前記プローブでの超音波の受波を捉えて超音波ビームを表わすＲＦ信号としての受信信号を生成する受信部と、
前記受信部で生成された受信信号の互いに隣接する２つのゼロクロス点に挟まれた半波形の代表振幅値および代表時刻を各半波形ごとに算出することにより、半波形ごとの代表振幅値
【数１】
2015084977.tif 000091
および代表時刻
【数２】
2015084977.tif 000092
但し、ｍは、ｍ番目の受信信号であることを表わし、ｉは、そのｍ番目の受信信号中  のｉ番目の半波形であることを表わす。
からなる圧縮データを生成する圧縮部と、該圧縮部で生成された圧縮データを構成する代表振幅値
【数３】
2015084977.tif 000093
を対応する代表時刻
【数４】
2015084977.tif 000094
に並べたときの、波形どうしが近似している領域のずれ量に基づいて、被検体構成要素の動きを算出する算出部とを有する演算部と、
前記演算部で得られた前記被検体構成要素の動きを表示する表示部とを備えたことを特徴とする超音波診断装置。
【請求項２】
前記受信部が、所定のサンプリングレートでのサンプリングにより得られたデータ列としてのＲＦ信号からなる受信信号を生成するものであって、
前記圧縮部が、前記受信信号を、
【数５】
2015084977.tif 000095
但し、ｍは、ｍ番目の受信信号であることを表わし、ｎは、そのｍ番目の受信信号中
のｎ番目のサンプリングデータであることを表わす。
としたとき、半波形ごとの代表振幅値
【数６】
2015084977.tif 000096
および半波形ごとの代表時刻
【数７】
2015084977.tif 000097
但し、ｍは、ｍ番目の受信信号であることを表わし、ｉは、そのｍ番目の受信信号中のｉ番目の半波形であることを表わす。
を、式
【数８】
2015084977.tif 000098
【数９】
2015084977.tif 000099
但し、
【数１０】
2015084977.tif 000100
は、半波形の、それぞれ開始点、終了点を表わし、
【数１１】
2015084977.tif 000101
は、受信信号のサンプリング周期を表わす。
に従って算出するものであることを特徴とする請求項１記載の超音波診断装置。
【請求項３】
前記算出部が、
【数１２】
2015084977.tif 000102
但し、ｉは観測点、
Ｌは奇数であって、
【数１３】
2015084977.tif 000103
は、波形のセグメントの範囲、
【数１４】
2015084977.tif 000104
は、サーチ範囲、
【数１５】
2015084977.tif 000105
は、サーチ範囲内での最小誤差を表わす。
に従って、最小誤差
【数１６】
2015084977.tif 000106
を満たす
【数１７】
2015084977.tif 000107
を求め、
【数１８】
2015084977.tif 000108
に従って時間シフト
【数１９】
2015084977.tif 000109
を求め、
【数２０】
2015084977.tif 000110
但し、
【数２１】
2015084977.tif 000111
は音速、
【数２２】
2015084977.tif 000112
は送受信の繰返し周期を表わす。
に従って、前記被検体構成要素の動きの速度
【数２３】
2015084977.tif 000113
を算出するものであることを特徴とする請求項１又は２記載の超音波診断装置。
【請求項４】
前記演算部が、前記被検体構成要素としての血液速度を求めるものであることを特徴とする請求項１から３のうちのいずれか１項記載の超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体内に超音波ビームを送波し反射して戻ってきた超音波を受波して受信信号を得、その受信信号に基づいて、超音波のドプラ遷移に基づく、被検体内の被検体構成要素の動きを算出して表示する超音波診断装置に関する。
【背景技術】
【０００２】
被検体、とりわけ人体の内部の画像を写し出す装置の１つとして、超音波を送受信して得た受信信号に基づいて画像を表示する超音波診断装置が知られている。その超音波診断装置には、通常、受信信号に基づいて被検体内の観察領域内の血流分布をカラー表示する機能が備えられている。
【０００３】
この超音波診断装置では、電圧印加を受けて振動して超音波を送波し、また超音波による振動を受けて電圧信号を発生する超音波振動子が多数個配列された超音波探触子を備えたプローブが使われる。このプローブの超音波探触子を被検体の体表に宛てがい、超音波探触子を構成している多数個の超音波振動子のそれぞれに、所定の遅延パターンに従ってそれぞれ遅延された、中心周波数ｆ０の複数のパルスからなるバースト波信号を印加する。すると、その超音波探触子から被検体内に、中心周波数ｆ０の、所定の深さ位置に焦点を結ぶ超音波ビームが所定の方向に送波される。そしてその反射超音波を、超音波探触子を構成する複数の超音波振動子のそれぞれでピックアップして複数の信号を得、それら複数の信号を所定の遅延パターンに従ってそれぞれ遅延させて互いに加算する。これにより、被検体内に延びる超音波ビームを表わす、ＲＦ信号としての受信信号が得られる。この超音波送受信が複数回繰り返され、その間の超音波のドプラ遷移による位相の変化Φ（ｔ）が求められて、その位相の変化Φ（ｔ）と中心周波数ｆ０とから、
【０００４】
【数１】
2015084977.tif 000003
【０００５】
但し、Ｔは送受信の繰返し周期
Ｃは音速
Ｖｄ（ｔ）は超音波ビーム方向のドプラ速度（血流速度）
ｔは時刻
である。
が算出される。そして、観察領域内の各点の速度Ｖｄを、通常は体表に近づく向きの血流が赤、遠ざかる向きの血流が青で、かつＶｄの大きさを色の輝度で表現する。
【０００６】
この演算法は、複素自己相関法と呼ばれる。
【０００７】
ここで、（１）式において、位相差Φ（Ａ）は、
【０００８】
【数２】
2015084977.tif 000004
【０００９】
の制限を受ける。この制限により、（１）式で算出される速度Ｖｄ（ｔ）に関し、算出可能な最高速度が制限される。
【００１０】
この複素自己相関法に代えて、時間領域の相互相関法（ＣＣＭ）を採用することが考えられる。このＣＣＭを採用すると、上記のような制限を受けず高速な血流速度も算出可能である。
【００１１】
このＣＣＭでは、２つの受信信号を
【００１２】
【数３】
2015084977.tif 000005
【００１３】
但し、ｍは、ｍ番目の受信信号であることを表わし、ｎは、そのｍ番目の受信信号中  のｎ番目のサンプリングデータであることを表わす。
としたとき、自己相関演算
【００１４】
【数４】
2015084977.tif 000006
【００１５】
但し、
【００１６】
【数５】
2015084977.tif 000007
【００１７】
は、波形のセグメント、
【００１８】
【数６】
2015084977.tif 000008
【００１９】
は、サーチ範囲、
【００２０】
【数７】
2015084977.tif 000009
【００２１】
は、サーチ範囲内での最大の相関値を表わす。
に従って、最大相関値
【００２２】
【数８】
2015084977.tif 000010
【００２３】
を満たす
【００２４】
【数９】
2015084977.tif 000011
【００２５】
を求め、
【００２６】
【数１０】
2015084977.tif 000012
【００２７】
但し、
【００２８】
【数１１】
2015084977.tif 000013
【００２９】
は、サンプリング周期を表わす。
に従って時間シフト
【００３０】
【数１２】
2015084977.tif 000014
【００３１】
を求め、
【００３２】
【数１３】
2015084977.tif 000015
【００３３】
但し、
【００３４】
【数１４】
2015084977.tif 000016
【００３５】
は音速、
【００３６】
【数１５】
2015084977.tif 000017
【００３７】
は、送受信の繰り返し周期を表わす。
に従って、速度
【００３８】
【数１６】
2015084977.tif 000018
【００３９】
を算出する。
【００４０】
この演算法を採用すれば、前述の複素自己相関法を採用したときのような制限を受けずに高速な速度も算出可能である。
【００４１】
ただし、この時間領域の相互相関法（ＣＣＭ）も従来から知られている演算法であるが、速度の検出が信号振幅の強さに依存し、弱い信号に対する正しい速度検出率が低い。また、計算量も膨大であり演算に時間がかかり過ぎるため、実際の装置では、採用され難い。
【００４２】
信号振幅の差を用いて、血流速度の検出精度が高い方法として、上記の（３）式に代えて、
【００４３】
【数１７】
2015084977.tif 000019
【００４４】
但し、
【００４５】
【数１８】
2015084977.tif 000020
【００４６】
は、サーチ範囲内での最小誤差を表わす。
に従って、最小誤差
【００４７】
【数１９】
2015084977.tif 000021
【００４８】
を満たす
【００４９】
【数２０】
2015084977.tif 000022
【００５０】
を求め、その後、上記の（４），（５）式に従って速度
【００５１】
【数２１】
2015084977.tif 000023
【００５２】
を算出することが知られている（特許文献１参照）。
【００５３】
（６）式を採用すると、（３）式を採用する場合と比べ、血流速度の検出精度が上げられる。
【００５４】
しかしながら、この（６）式を採用したとしても、計算量はまだまだ膨大であり、複素自己相関法を採用したときのようなリアルタイム性は望めない。
【先行技術文献】
【特許文献】
【００５５】
【特許文献１】特開２００１−２８６４７２号公報
【発明の概要】
【発明が解決しようとする課題】
【００５６】
本発明は、上記事情に鑑み、計算量が少なく、かつ速い速度まで算出が可能な演算法を採用した超音波診断装置を提供することを目的とする。
【課題を解決するための手段】
【００５７】
上記目的を達成する本発明の超音波診断装置は、
被検体内への超音波ビームの送波と該被検体内で反射して戻ってきた超音波の受波とを担うプローブと、
複数のパルス信号をプローブに送信してプローブに超音波ビームを送波させる送信部と、
プローブでの超音波の受波を捉えて超音波ビームを表わすＲＦ信号としての受信信号を生成する受信部と、
受信部で生成された受信信号の互いに隣接する２つのゼロクロス点に挟まれた半波形の代表振幅値および代表時刻を各半波形ごとに算出することにより、半波形ごとの代表振幅値
【００５８】
【数２２】
2015084977.tif 000024
【００５９】
および代表時刻
【００６０】
【数２３】
2015084977.tif 000025
【００６１】
但し、ｍは、ｍ番目の受信信号であることを表わし、ｉは、そのｍ番目の受信信号中  のｉ番目の半波形であることを表わす。
からなる圧縮データを生成する圧縮部と、その圧縮部で生成された圧縮データを構成する代表振幅値
【００６２】
【数２４】
2015084977.tif 000026
【００６３】
を対応する代表時刻
【００６４】
【数２５】
2015084977.tif 000027
【００６５】
に並べたときの、波形どうしが近似している領域のずれ量に基づいて、被検体構成要素の動きを算出する算出部とを有する演算部と、
演算部で得られた被検体構成要素の動きを表示する表示部とを備えたことを特徴とする。
【００６６】
本発明の超音波診断装置では、半波形ごとの代表振幅値
【００６７】
【数２６】
2015084977.tif 000028
【００６８】
と代表時刻
【００６９】
【数２７】
2015084977.tif 000029
【００７０】
が算出され、これによりデータ量が大幅に圧縮される。そのデータ量が大幅に圧縮されたデータを使って被検体構成要素の動きを算出することにより、時間領域で演算を行なう従来の演算法と比べ、遥かに少ない計算量で済む。また、時間領域で演算を行なうことにより、複素自己相関法のような算出可能速度の制限を受けずに済み、速い速度まで算出可能である。
【００７１】
ここで、本発明の超音波診断装置において、上記受信部が、所定のサンプリングレートでのサンプリングにより得られたデータ列としてのＲＦ信号からなる受信信号を生成するものであって、
上記圧縮部が、受信信号を、
【００７２】
【数２８】
2015084977.tif 000030
【００７３】
但し、ｍは、ｍ番目の受信信号であることを表わし、ｎは、そのｍ番目の受信信号中
のｎ番目のサンプリングデータであることを表わす。
としたとき、半波形ごとの代表振幅値
【００７４】
【数２９】
2015084977.tif 000031
【００７５】
および半波形ごとの代表時刻
【００７６】
【数３０】
2015084977.tif 000032
【００７７】
但し、ｍは、ｍ番目の受信信号であることを表わし、ｉは、そのｍ番目の受信信号中のｉ番目の半波形であることを表わす。
を、式
【００７８】
【数３１】
2015084977.tif 000033
【００７９】
【数３２】
2015084977.tif 000034
【００８０】
但し、
【００８１】
【数３３】
2015084977.tif 000035
【００８２】
は、半波形の、それぞれ開始点、終了点を表わし、
【００８３】
【数３４】
2015084977.tif 000036
【００８４】
は、受信信号のサンプリング周期を表わす。
に従って算出するものであることが好ましい。
【００８５】
（７）式に基づく代表振幅値および（８）式に基づく代表時刻を採用すると、後述するように、このような代表振幅値および代表時刻を算出することなく元の受信信号をそのまま使って速度を算出した場合とほぼ同程度の精度で速度が算出される。
【００８６】
また、本発明の超音波診断装置において、上記算出部が、
【００８７】
【数３５】
2015084977.tif 000037
【００８８】
但し、ｉは観測点、
Ｌは奇数であって、
【００８９】
【数３６】
2015084977.tif 000038
【００９０】
は、波形のセグメントの範囲、
【００９１】
【数３７】
2015084977.tif 000039
【００９２】
は、サーチ範囲、
【００９３】
【数３８】
2015084977.tif 000040
【００９４】
は、サーチ範囲内での最小誤差を表わす。
に従って、最小誤差
【００９５】
【数３９】
2015084977.tif 000041
【００９６】
を満たす
【００９７】
【数４０】
2015084977.tif 000042
【００９８】
を求め、
【００９９】
【数４１】
2015084977.tif 000043
【０１００】
に従って時間シフト
【０１０１】
【数４２】
2015084977.tif 000044
【０１０２】
を求め、
【０１０３】
【数４３】
2015084977.tif 000045
【０１０４】
但し、
【０１０５】
【数４４】
2015084977.tif 000046
【０１０６】
は音速、
【０１０７】
【数４５】
2015084977.tif 000047
【０１０８】
は送受信の繰返し周期を表わす。
に従って、被検体構成要素の動きの速度
【０１０９】
【数４６】
2015084977.tif 000048
【０１１０】
を算出するものであることが好ましい。
【０１１１】
（９）式は、代表振幅値および代表時刻を採用している点を除き、演算法自体は前述の（６）式と同じである。（６）式の最小自乗法を採用すると、（３）式の相互相関法を採用した場合と比べ、信号振幅の強さに依存せず、血流速度等、被検体構成要素の動きの速度の検出精度が高くなる。そこで、代表振幅値および代表時刻をこの（６）式の演算法にあてはめた（９）式を採用すると、代表振幅値および代表時刻を（３）式にあてはめるよりも計算量が少なくて済み、さらなる高速演算が可能となる。
【０１１２】
ここで、本発明の超音波診断装置において、上記演算部が、被検体構成要素としての血流速度を求めるものであることが好ましい。
【０１１３】
超音波診断装置では、血流速度が１つの大きな診断材料であり、上記の演算法を血流速度の算出に採用することが好ましい。
【発明の効果】
【０１１４】
以上の本発明によれば、計算量が少なく、かつ速い速度まで算出が可能な演算法を採用した超音波診断装置が実現する。
【図面の簡単な説明】
【０１１５】
【図１】本実施形態の超音波診断装置の構成を表わすブロック図である。
【図２】複数回（ここではＭ回）の超音波ビームを送波するときに送信部で生成されるパルス信号送受信間隔の模式図である。
【図３】受信信号（Ａ）とその受信信号から算出された代表振幅値および代表時刻（Ｂ）を示した図である。
【図４】受信信号からの代表振幅値および代表時刻の算出方法の説明図である。
【図５】（１０）式に従う演算の説明図である。
【図６】シミレーション結果を示した図である。
【発明を実施するための形態】
【０１１６】
以下、本発明の実施の形態を説明する。
【０１１７】
図１は、本実施形態の超音波診断装置の構成を表わすブロック図である。
【０１１８】
この超音波診断装置１には、プローブ１０と、送信部１１と、受信部１２と、制御部１３が設けられている。
【０１１９】
プローブ１０には、多数の超音波振動子が配列された超音波探触子（不図示）が備えられている。その超音波探触子が、被検体としての人体の体表に宛がわれる。
【０１２０】
送信部１１はパルス信号を生成し、プローブ１０の超音波探触子を構成する多数の超音波振動子それぞれに向けてパルス信号を送信する。このパルス信号は、１つのパルス信号につき、中心周波数ｆ０であって長さが複数波長のバースト波信号である。この生成されたバースト波信号は、人体内の超音波ビームの延びる向き、およびその超音波ビームの焦点の深さに応じて定められる遅延パターンに従ってそれぞれ遅延されて、各超音波振動子に印加される。すると、それら多数の超音波振動子それぞれから超音波が送波され、それらの超音波の干渉作用により、所望の向きに延び、かつ所望の深さ位置に焦点を持つ超音波ビームが人体内に送り込まれる。人体内に送り込まれた超音波ビームは、その超音波ビームが人体内の浅い位置から深い位置へと進むに従って人体内の各深さ位置で反射されて超音波探触子に戻り、多数の超音波振動子のそれぞれで受波される。したがってこの受波により得られる信号は、時間軸が超音波探触子からの人体内の深さに対応している。
【０１２１】
受信部１２では、それら多数の超音波振動子での受波により、各超音波振動子で得られた信号を、これも人体内において超音波ビームの延びる向きや焦点の深さ等に応じた遅延パターンに従ってそれぞれ遅延させて互いに加算する。こうすることにより、人体内を所望の向きに延びる超音波ビームを表わす、ＲＦ（Ｒａｄｉｏ  Ｆｒｅｑｕｅｎｃｙ）信号としての受信信号が生成される。この受信部１２で生成された受信信号は、Ｂモード処理部１４とＭＴＩフィルタ１５に入力される。
【０１２２】
制御部１３は、送信部１１からのバースト波信号の送信タイミングや、受信部１２での受信タイミングを制御する。また制御部１３は、この他にも、この超音波診断装置１の動作全体の制御を担っている。
【０１２３】
表示部１９にＢモード像、すなわち人体内の超音波反射率分布に基づく画像を表示するときは、送信部１１では、プローブ１０から人体内に向けて、人体内の観察領域内において順次に向きの異なる超音波ビームが送波されるように、順次異なる遅延パターンに基づいて遅延させたバースト波信号がプローブ１０に送信される。受信部１２でも同様に、その超音波ビームが送波された向きの超音波ビームが生成されるように遅延加算される。
【０１２４】
受信部１２でこのようにして得られた受信信号は、Ｂモード処理部１４に入力されてＢモード用の画像処理が施され、座標変換部１７に入力される。
【０１２５】
受信部１２で得られる受信信号は、各超音波ビームに沿う深さ方向を時間軸とする信号である。これに対して、表示部１９では、ラスタスキャンの向きに並ぶピクセルデータからなる画像信号に基づいて画像が表示される。そこで、座標変換部１７では、Ｂモード処理部１４から入力されてきた信号が表示部１９での表示に適したピクセルデータの配列からなる信号に変換される。
【０１２６】
この座標変換部１７から出力された信号は画像合成部１８に入力される。この画像合成部１８では、Ｂモード像と人体内の血流を表わすカラードプラ画像を合成する。人体内の血流を表わすカラードプラ画像は、通常は、Ｂモード像に重ねて表示される。カラードプラ画像を表示することなく、Ｂモード像のみ表示することもある。
【０１２７】
この画像合成部１８で合成された画像は表示部１９に入力され、図示しない表示画面上に画像が表示され、診断に供される。
【０１２８】
また、人体内の血流を表わすカラードプラ画像を表示させるときは、上記のＢモード像の生成に加えて以下の処理が行われる。
【０１２９】
プローブ１０からは人体内の同一の向きに延びる超音波ビームが複数回（例えば８回）送波され、かつプローブ１０で受波されるように、送信部１１によるバースト波信号の遅延、受信部１２での遅延加算が行われる。
【０１３０】
この人体内の同一の向きへの送受信が、観察領域内に延びる超音波ビームの各向きそれぞれについて行われる。
【０１３１】
このようにして得られた受信信号はＭＴＩフィルタ１５に入力される。ＭＴＩフィルタ１５はハイパスフィルタの一種であり、ここでは血流のドプラ遷移に起因する成分が通過され、臓器等の動きに起因する成分がカットされる。
【０１３２】
ＭＴＩフィルタ１５を通過した受信信号は、カラードプラ処理部１６に入力される。このカラードプラ処理部１６では入力されてきた受信信号に基づいて、後述する演算により血流速度等が算出される。この血流速度等は、超音波ビーム上の各点ごとに行われる。超音波ビームは観察領域内で順次向きを変えるため、観察領域内の２次元的な各点について血流速度等が算出される。この算出された血流速度等を表わすデータは座標変換部１７に入力されて座標変換を受け、画像合成部１８においてＢモード像に重ねられる。そして表示部１９でＢモード像に重ねられた血流分布が表示される。この血流分布は、通常は、プローブ１０の、体表に宛てがわれている超音波探触子に向かう向きの血流が赤、遠ざかる向きの血流が青で表示される。さらに血流の速度が、赤あるいは青の色の輝度で表示される。このようにして、観察領域内の血流分布が表示されて診断に供される。
【０１３３】
以下、このカラードプラ処理部１６での処理を詳述する。
【０１３４】
図２は、複数回（ここではＭ回）の超音波ビームを送波するときに送信部で生成されるパルス信号送受信間隔の模式図である。
【０１３５】
ここでは、バースト波信号１，２，…，Ｍ−１，ＭからなるＭ個のパルス信号が示されている。この図では、パルス信号が単純化されて示されているが、パルス信号１つずつが、中心周波数ｆ０であって長さが複数波長のバースト波信号である。ここでは、パルス信号の繰返し周期をＴ、繰返し周波数をＰＲＦと称する。Ｔ＝１／ＰＲＦの関係にある。
【０１３６】
送信部１１でこの図２に模式的に示すバースト波信号がプローブ１０に同一の遅延パターンで送信されたときの受信部１２でのＲＦ（Ｒａｄｉｏ  Ｆｒｅｑｕｅｎｃｙ）としての受信信号が、所定のサンプリングレート（例えばサンプリング周波数２４ＭＨｚ）でサンプリングされ、ＭＴＩフィルタ１５でクラッタ成分除去用のフィルタリング処理がなされて、カラードプラ処理部１６に入力される。ここでは、カラードプラ処理部１６に入力されてきた受信信号を、
【０１３７】
【数４７】
2015084977.tif 000049
【０１３８】
但し、ｍは、ｍ番目の受信信号であることを表わし、ｎは、そのｍ番目の受信信号中
のｎ番目のサンプリングデータであることを表わす。
とする。
【０１３９】
図３は、受信信号（Ａ）とその受信信号から算出された代表振幅値および代表時刻（Ｂ）を示した図である。
【０１４０】
図１に示すカラードプラ処理部１６では、入力されてきた受信信号
【０１４１】
【数４８】
2015084977.tif 000050
【０１４２】
に基づいて、その受信信号の互いに隣接する２つのゼロクロス点に挟まれた半波形の代表振幅値および代表時刻を各半波形ごとに算出することにより、半波形ごとの代表振幅値
【０１４３】
【数４９】
2015084977.tif 000051
【０１４４】
および代表時刻
【０１４５】
【数５０】
2015084977.tif 000052
【０１４６】
但し、ｍは、ｍ番目の受信信号であることを表わし、ｉは、そのｍ番目の受信信号中  のｉ番目の半波形であることを表わす。
が算出される。この算出にあたっては、本実施形態では、式
【０１４７】
【数５１】
2015084977.tif 000053
【０１４８】
【数５２】
2015084977.tif 000054
【０１４９】
但し、
【０１５０】
【数５３】
2015084977.tif 000055
【０１５１】
は、半波形の、それぞれ開始点、終了点を表わし、
【０１５２】
【数５４】
2015084977.tif 000056
【０１５３】
は、受信信号のサンプリング周期を表わす。
が採用されており、これら（７），（８）式に従って算出される。
【０１５４】
図４は、受信信号からの代表振幅値および代表時刻の算出方法の説明図である。
【０１５５】
ここでは、図４（Ａ）に実線で示した半波形についての代表振幅値および代表時刻を例に挙げて説明する。
【０１５６】
ここで、半波形は、２つのゼロクロス点、すなわち、図４（Ａ）の例では、ｎ＝ｎｓの開始点からｎ＝ｎｅの終了点に挟まれた部分の波形をいう。図４（Ａ）にサンプリング周期
【０１５７】
【数５５】
2015084977.tif 000057
【０１５８】
間隔の縦線の長さが、サンプリングデータ
【０１５９】
【数５６】
2015084977.tif 000058
【０１６０】
である。
【０１６１】
ここでは、（８）式に従って、そのサンプリングデータ
【０１６２】
【数５７】
2015084977.tif 000059
【０１６３】
を開始点ｎ＝ｎｓから終了点ｎ＝ｎｅまで加算し、その加算値をこの半波形ｉの代表振幅値
【０１６４】
【数５８】
2015084977.tif 000060
【０１６５】
とする。
【０１６６】
また、（８）式に従う、サンプリングデータ
【０１６７】
【数５９】
2015084977.tif 000061
【０１６８】
によって重み付けられた時刻をこの半波形ｉを代表する代表時刻
【０１６９】
【数６０】
2015084977.tif 000062
【０１７０】
とする。
【０１７１】
すなわち、この演算法により、この半波形ｉのデータは、図４（Ｂ）に示すように、代表時刻
【０１７２】
【数６１】
2015084977.tif 000063
【０１７３】
における代表振幅値
【０１７４】
【数６２】
2015084977.tif 000064
【０１７５】
というデータに圧縮される。
【０１７６】
図３（Ｂ）は、図３（Ａ）に示す受信信号の各半波形について上記の演算を行い、各代表振幅値を対応する各代表時刻に配置したデータ列をあらわしている。
【０１７７】
図１のカラードプラ処理部１６では、受信信号（図３（Ａ））から代表振幅値および代表時刻からなる圧縮データ（図３（Ｂ））を算出した後、次に、その算出された圧縮データを構成する代表振幅値
【０１７８】
【数６３】
2015084977.tif 000065
【０１７９】
を対応する代表時刻
【０１８０】
【数６４】
2015084977.tif 000066
【０１８１】
に並べたときの、波形どうしが近似している領域のずれ量に基づいて、人体内の血流速度が算出される。
【０１８２】
この血流速度の算出にあたり、本実施形態のカラードプラ処理部１６では、先ず、ｍ番目とｍ＋１番目の受信信号の代表振幅値
【０１８３】
【数６５】
2015084977.tif 000067
【０１８４】
を使い、
【０１８５】
【数６６】
2015084977.tif 000068
【０１８６】
但し、ｉは観測点、
Ｌは奇数であって、
【０１８７】
【数６７】
2015084977.tif 000069
【０１８８】
は、波形のセグメントの範囲、
【０１８９】
【数６８】
2015084977.tif 000070
【０１９０】
は、サーチ範囲、
【０１９１】
【数６９】
2015084977.tif 000071
【０１９２】
は、サーチ範囲内での最小誤差を表わす。
に従って、最小誤差
【０１９３】
【数７０】
2015084977.tif 000072
【０１９４】
を満たす
【０１９５】
【数７１】
2015084977.tif 000073
【０１９６】
を求める。さらに、
【０１９７】
【数７２】
2015084977.tif 000074
【０１９８】
に従って時間シフト
【０１９９】
【数７３】
2015084977.tif 000075
【０２００】
を求める。そして、
【０２０１】
【数７４】
2015084977.tif 000076
【０２０２】
但し、
【０２０３】
【数７５】
2015084977.tif 000077
【０２０４】
は音速、
【０２０５】
【数７６】
2015084977.tif 000078
【０２０６】
は送受信の繰り返し周期を表わす。
に従って、血流の速度
【０２０７】
【数７７】
2015084977.tif 000079
【０２０８】
を算出する。
【０２０９】
図５は、（１０）式に従う演算の説明図である。
【０２１０】
図５（Ａ），（Ｂ）は、それぞれ、ｍ番目、ｍ＋１番目の受信信号の代表振幅値をそれに対応する代表時刻に並べた波形を示した図である。
【０２１１】
ここでは、図５（Ａ）に示すデータ列の、枠ａ内の波形と近似した波形を図５（Ｂ）に示すデータ列の中から探索する。この枠ａ中の波形（データ列）は、近似した波形を探索するためのデータ列のセグメントであり、ここでは、
【０２１２】
【数７８】
2015084977.tif 000080
【０２１３】
の幅の枠ａが設定されている。この枠ａの幅は、この中に代表振幅値が３又は５程度含まれる幅である。
【０２１４】
図５（Ｂ）に示すデータ列においても、枠ａと同じ幅の枠ｂが設定される。
【０２１５】
そして、（１０）式に従い、枠ａ内のセグメントと枠ｂ内のセグメントの一致の程度が算出される。
【０２１６】
誤差
【０２１７】
【数７９】
2015084977.tif 000081
【０２１８】
が小さいほど一致していることを意味している。上記の枠ｂを、サーチ範囲
【０２１９】
【数８０】
2015084977.tif 000082
【０２２０】
の中で順次移動させて移動した各位置において、誤差
【０２２１】
【数８１】
2015084977.tif 000083
【０２２２】
が算出され、それらの中の最小誤差
【０２２３】
【数８２】
2015084977.tif 000084
【０２２４】
が得られる位置
【０２２５】
【数８３】
2015084977.tif 000085
【０２２６】
が検出される。このようにして、
【０２２７】
【数８４】
2015084977.tif 000086
【０２２８】
が検出されると、次に上記の（１１）式に従って、時間シフト
【０２２９】
【数８５】
2015084977.tif 000087
【０２３０】
が算出され、（１２）式に従って血流速度
【０２３１】
【数８６】
2015084977.tif 000088
【０２３２】
が算出される。
【０２３３】
図６は、シミレーション結果を示した図である。
【０２３４】
図６（Ａ）は、カラードプラ処理部１６に入力されてきた受信信号
【０２３５】
【数８７】
2015084977.tif 000089
【０２３６】
をそのまま使い、（６）式に示す最小自乗法で計算したときの血流速度を示している。
【０２３７】
また図６（Ｂ）は、カラードプラ処理部１６に入力されてきた受信信号
【０２３８】
【数８８】
2015084977.tif 000090
【０２３９】
から（７）式，（８）式に従って半波形ごとの代表振幅値と代表時刻を算出し、（９）式に示す最小自乗法で計算したときの血流速度を示している。
【０２４０】
また、図６（Ａ），（Ｂ）とも、左側の図（ａ）は一定の流速についてのシミュレーション結果、右側の図（ｂ）は、血管壁付近の流れ速度が遅く、血管中心の速度が速い層流の場合のシミュレーション結果を示している。
【０２４１】
ここでは、送信周波数ｆ０＝２ＭＨｚ、サンプリング周波数Ｆｓ＝２４ＭＨｚを採用している。
【０２４２】
図６（Ａ）と図６（Ｂ）を比較すると、血流速度の分布と算出精度（縦軸）は、ほぼ同じである。更に、圧縮された代表時刻の間の血流速度を検出する必要がある場合、その代表時刻の隣接の二つの速度値の線形補間により高速に求めることが可能である。
【０２４３】
一方、（６）式、（１０）式に従って算出される差の自乗和を求める計算量は、１回の計算（図５に示した枠ａと枠ｂとの間での計算）につき、ほぼ１／６となり、また、サーチ範囲内で枠ｂを動かしながら行なう計算の回数も１／６となる。
【０２４４】
尚、本実施形態では、半波形ごとの代表振幅値および代表時刻を、（７）式、（８）式に従って算出する例を示したが、代表振幅値および代表時刻は、（７）式、（８）式に限定されるものではない。例えば半波形ごとのサンプリングデータのピークの値を代表振幅値とし、そのピークのサンプリングデータのサンプリング時刻を代表時刻としてもよい。
【０２４５】
また、本実施形態では、（１０）式に従う最小自乗和の演算に従って時間シフト（（１１）式参照）を求めているが、本発明は、代表振幅値および代表軸を求めることによってデータ量を圧縮することを旨とするものであり、時間シフトを求めるにあたり、必ずしも（１０）式に従う最小自乗和の演算を採用する必要はなく、時間シフトを求めるための様々な演算法を採用することができる。
【０２４６】
さらに、本実施形態では、血流速度を算出することを念頭に置いて説明したが、例えば臓器の動き速度等、血流以外の被検体構成要素の動きの算出のために、上記の演算を採用してもよい。
【符号の説明】
【０２４７】
１    超音波診断装置
１０    プローブ
１１    送信部
１２    受信部
１３    制御部
１４    Ｂモード処理部
１５    ＭＴＩフィルタ
１６    カラードプラ処理部
１７    座標変換部
１８    画像合成部
１９    表示部
【図１】
2015084977.tif 000114
【図２】
2015084977.tif 000115
【図３】
2015084977.tif 000116
【図４】
2015084977.tif 000117
【図５】
2015084977.tif 000118
【図６】
2015084977.tif 000119