JP5950599B2 - 被検体情報取得装置、被検体情報取得方法及びプログラム - Google Patents

Description

本発明は、被検体情報取得装置、被検体情報取得方法及びプログラムに関する。特に、被検体に弾性波を送信し、被検体内で反射した弾性波を受信して被検体内の情報を取得する技術に関する。

医療等の現場において被検体情報取得装置である超音波診断装置は広く用いられている。超音波診断装置は生体内の音響インピーダンス分布を反映した形態情報だけでなく、ドップラ技術を用いた血流速度情報などの対象物の移動情報も取得可能である。

血流速度情報を取得する際には、血液が存在する領域からの反射波を用いる。この血液が存在する領域からの反射波は、主に血液に含まれる赤血球から反射、散乱した弾性波（典型的には超音波）である。一般的な超音波診断装置で用いられる超音波の波長は赤血球のサイズよりも長いため、一つ一つの赤血球を区別することはできない。結果として、血液が存在する領域からの反射波は、複数の赤血球で形成される散乱体の集合体（散乱体群）の様相（個々の赤血球の形態や向き、相対位置など）を反映したものとなる。

赤血球で形成される散乱体群からの反射波を用いて血流速度情報を抽出するには、その反射波のドップラシフト周波数を求める技術がよく用いられる。しかしこのようなドップラシフト周波数を測定する手法は原理的には、血流速度の内、超音波ビームの送受信方向（走査線方向）への射影成分のみが測定可能である。つまり、本来の血流速度を求めるためには血流方向と超音波の送受信方向とで成す角度を考慮した補正をする必要がある。この補正は、血流方向と超音波の送受信方向とで成す角度の余弦によって、ドップラシフト周波数から推定された流速を除する処理であり、その角度が大きくなればなるほど誤差も大きくなる可能性が高まる。

そこで、特許文献１には、超音波の送受信方向に対して垂直な方向に進行する血流であっても流速が計測可能な例が開示されている。特許文献１は断層像内の２点について、超音波ビームの反射／散乱波形に関する２点間の相互相関値を計算し、２点間の相互相関値がピークとなる時間で２点間の距離を除することにより、血流速度を取得する。

特開昭６１−１８７８４３号公報

ところで、血液中の赤血球は血管内を流れると共にその分布を変化させてゆく。つまり、たとえ同じ赤血球が複数集まって構成される散乱体群であっても、赤血球の形態や向き、相対的な配置が変化していく。この変化により、異なる時間で取得された反射波の受信信号（反射波を受信することにより取得される受信信号）を用いて算出された相互相関値（相互相関係数）は、取得された時間の差が大きくなるとともに低下してゆく。

このような取得時間差の拡大と共に相互相関値が低下してゆく現象が生じる場合、その散乱体群の相対位置変化を含む散乱体群の時間変化を勘案して血流速度情報を推定しないと誤差を生じることがある。

本発明は、上記課題に鑑み、被検体内の対象物である散乱体群の時間変化をも考慮して血流速度等の移動情報を算出することで、対象物の移動情報の取得精度を向上させることを目的とする。

本発明の被検体情報取得装置は、被検体内の対象物の移動情報を取得する被検体情報取得装置であって、
被検体に弾性波が送信され、前記被検体内の各位置で反射された反射波を受信した複数の変換素子から出力される複数の受信信号を用いて、前記被検体内の各位置からの反射波に対応する信号を走査線信号として複数取得する走査線信号取得手段と、
複数の前記走査線信号を用いて、前記対象物の移動情報を取得する処理部を有し、
前記処理部は、異なる位置の走査線信号同士の相互相関値を複数用いて、時間差の軸と距離差の軸との２軸で表される平面上における相互相関値の等高線に基づいて前記対象物の移動情報を取得することを特徴とする。

また、本発明の被検体情報取得方法は、被検体内の対象物の移動情報を取得する被検体情報取得方法であって、
被検体に弾性波が送信され、前記被検体内の各位置で反射された反射波を受信した複数の変換素子から出力される複数の受信信号を用いて、前記被検体内の各位置からの反射波に対応する信号を走査線信号として複数取得するステップと、
複数の前記走査線信号を用いて、前記対象物の移動情報を取得するステップと、を有し前記移動情報を取得するステップでは、異なる位置の走査線信号同士の相互相関値を複数用いて、時間差の軸と距離差の軸との２軸で表される平面上における相互相関値の等高線に基づいて前記対象物の移動情報を取得することを特徴とする。

本発明によれば、被検体内の対象物である散乱体群の時間変化を考慮し、対象物の移動情報を高精度に取得することが可能となる。

本発明による被検体情報取得装置のシステム概要を示すブロック図である。 相互相関値の等高線を示す図である。 被検体及び散乱体群を模式的に示した図である。 本発明による処理の概念を説明するための図である。 本発明による処理ステップを示したフローチャートである。 被検体及び散乱体群を模式的に示した図である。 第１の実施形態の処理の概念を説明するための図である。 第１の実施形態の表示の一例を示した模式図である。 第２の実施形態の処理の概念を説明するための図である。 第２の実施形態の処理ステップを示したフローチャートである。 第３の実施形態のシステム概要を示すブロック図である。 速度分布を取得する処理の概念を説明するための図である。

本発明の実施の形態について、図面を用いて説明する。本発明において、弾性波とは、典型的には超音波であり、音波、超音波、音響波、と呼ばれる弾性波を含む。本発明の被検体情報取得装置とは、被検体に弾性波を送信し、被検体内部で反射した反射波（反射した弾性波）を受信して、被検体情報を画像データとして取得する装置を含む。取得される被検体情報とは、被検体内部の組織の音響インピーダンスの違いを反映した情報や、被検体内の対象物の移動情報である。被検体内の対象物の移動情報とは、赤血球で形成される散乱体群等の移動速度に相当する血流の速度や速度分布等の速度情報、被検体内の組織の変位や変位分布等の変位情報を含む。

（被検体情報取得装置の基本的な構成）
本発明が適用できる被検体情報取得装置の構成について図１を用いて説明する。図１は本発明による被検体情報取得装置のシステム概要を示す図である。本実施形態の被検体情報取得装置は、探触子００１、受信回路系００５、送信回路系００３、整相加算ブロック００６、相関算出ブロック００７、相関パターン処理ブロック００８、画像生成ブロック００９、画像処理ブロック０１０、システム制御部００４、画像表示手段０１１、を備える。探触子００１は弾性波を複数位置に送信して反射波を受信する送受信器であり、弾性波を電気信号に変換する変換素子００２を複数有する。なお、画像表示手段０１１は、本発明の被検体情報取得装置とは別に提供されていてもよい。

本発明において、整相加算手段である整相加算ブロック００６は、被検体内の各位置の信号（つまり、各位置で反射した弾性波に対応する信号）を走査線信号として複数取得する走査線信号取得手段である。また、相関算出手段である相関算出ブロック００７と、相関処理手段である相関パターン処理ブロック００８と、で少なくとも本発明の処理部を構成する。

送信信号生成手段である送信回路系００３は、システム制御部００４からの制御信号に従って、注目位置や注目方向に応じた遅延時間や振幅を有する電圧波形を送信信号として生成する。この送信信号は複数の変換素子００２によって弾性波に変換され、探触子００１から被検体０００内へと送信される。被検体０００内部で反射された反射波（弾性波）は、複数の変換素子００２によって複数の受信信号に変換される。複数の受信信号は受信信号処理手段である受信回路系００５に入力される。

受信回路系００５では複数の受信信号を増幅し、複数のデジタル信号に変換する。受信回路系００５から出力された複数のデジタル信号は整相加算ブロック００６に入力される。整相加算ブロック００６では、弾性波を送信した方向や位置に応じて、複数のデジタル信号に対する遅延処理を行いさらに加算する、つまり整相加算処理を実行する。この整相加算処理は、被検体内の各位置で反射した反射波の音圧に対応する信号を走査線信号として生成する処理である。このように、複数の受信信号を用いて整相加算処理されることで得られる走査線信号は相関算出ブロック００７と画像生成ブロック００９とに入力される。

まず先に、画像生成ブロック００９での処理について説明する。画像生成ブロック００９では、入力された走査線信号に対して包絡線を算出し、包絡線信号として画像処理ブロック０１０へ出力する。この包絡線信号は夫々、被検体内の各位置の音響インピーダンスの違いを反映した信号である。また、画像生成ブロック００９では入力された走査線信号に対してバンドパスフィルタをかけるなどの各種処理を必要に応じて行ってもよい。画像処理ブロック０１０では入力された包絡線信号に対して強度調整や各種フィルタ処理を行い、組織の音響インピーダンスの違いを反映した分布情報として、インピーダンス輝度データ（いわゆるＢ−ｍｏｄｅ像）を画像表示系０１１に出力する。

次に、相関算出ブロック００７での処理について説明する。相関算出ブロック００７では、システム制御部００４からの指示に従い、入力された走査線信号同士の相互相関値（相互相関係数）を複数算出し、相関パターン処理ブロック００８へ出力する。相関パターン処理ブロック００８では入力された複数の相互相関値を用いて速度等の移動情報を算出し、画像処理ブロック０１０へ出力する。また、相関パターン処理ブロック００８は、速度の値そのものだけでなく、速度に関するパラメータや各位置の速度を示す速度分布を移動情報として取得しても良い。相関算出ブロック００７及び相関パターン処理ブロック００８の詳細な機能は、図３−５を用いて後述する。

画像処理ブロック０１０では、入力された速度等の移動情報を処理し、移動情報輝度データとして画像表示系０１１に出力する。画像表示系０１１は入力された移動情報輝度データならびにインピーダンス輝度データを、システム制御部００４からの指示に従い表示する。表示方法は、移動情報輝度データ及びインピーダンス輝度データを重畳して表示しても良いし、並べて表示してもよい。もちろん、夫々の輝度データのみを表示してもよい。表示モードは、システム制御部００４からの指示によって変更可能である。以上が、本実施形態の被検体情報取得装置の基本的な構成である。

（散乱体群の時間変化を加味した相互相関値の挙動）
次に、相関算出ブロック００７での具体的な処理を述べる前に、散乱体群からの反射波形の相互相関値の挙動に関して説明を行う。

散乱体群が時間変化（具体的には、散乱体の形態や向きの変化、相対的な位置変化）をする場合、ある注目位置に存在する散乱体群からの反射波形（走査線信号の波形に対応）と、そこからξだけ離れた注目位置に存在する散乱体群からの反射波形（走査線信号の波形に対応）とを考える。この場合、ある注目位置の走査線信号と、ある注目位置から距離ξだけ離れた注目位置の走査線信号と、の相互相関値は、ξが大きくなればなるほど、つまり２つの注目位置が離れるほど小さい値となる。また、散乱体群を構成する散乱体がその様相を変化（相対的な位置などを変化）させる場合、同じ散乱体で構成される散乱体群であっても異なる時間に反射された２つの反射波形同士の相互相関値は、時間の差τが大きくなればなるほど小さい値となる。つまり、注目位置が離れる、もしくは反射波の取得時間の差が大きくなるに従って相互相関値は低下してゆく。

図２（ａ）は、横軸を注目位置間の距離ξ、縦軸を反射波取得の時間差τにとって、相互相関値の等高線を引いた図を示す。散乱体群が移動していない場合、それぞれの軸に対称な楕円（円を含む）形状で近似できる等高線になる。軸の交点（原点）の位置、つまり注目位置間の距離０、時間差０の時に相互相関値は最大となる。

次に、散乱体群が速度ｖで移動している場合の相互相関値について考える。図２（ａ）の点Ｐ（注目位置間の距離差０、時間差τ０）に注目する。点Ｐは速度ｖで移動しているため、時間差τ０の間にｖ×τ０の距離を移動する。つまり、図２（ｂ）の点Ｐ’の位置（注目位置間の距離差ｖτ０、時間差τ０）に移動することになる。つまり、ＰとＰ’を結ぶ直線と同じ方向に散乱体群が移動すると、この点Ｐ’は点Ｐと同じ相互相関値を有する。結果として、図２（ａ）で示した相互相関値の等高線は、図２（ｂ）で示したように歪んだ楕円形状で近似できる。この時、図２（ｂ）内の一点鎖線で示した直線（原点とＰ’とを結ぶ直線）は傾き１／ｖを有する直線となる。つまり、この直線の傾きは速度の逆数であるため、実線で示した楕円形状に近似される形状を求めることで、一点鎖線で示した直線、ひいては散乱体群の速度を求めることができる。また、この一点鎖線で示した直線は、各時間差における相互相関値の最大値となる点が繋がれた線を示す。つまり、ある時間差における、相互相関値の最大値（つまりある時間差におけるピーク）に関する情報を求めることで、散乱体群の速度を求めることができる。ある時間差における相互相関値の「最大値に関する情報」とは、相互相関値が最大となる距離差の情報等、相互相関値のピークが特定できる情報であればよい。

このように、散乱体群自身の時間変化まで考慮した場合、横軸を注目位置間の距離差、縦軸を反射波取得の時間差とし、その２軸で表される平面上での相互相関値の分布を取得することで、散乱体群の移動情報を取得することができる。

（基本的な処理フロー）
次に、図３から図５までを用いて、相関算出ブロック００７ならびに相関パターン処理ブロック００８での処理の詳細を説明する。

図３は、血管３０１内を流れる赤血球から構成される散乱体群３０２を模式的に示した図であり、図４は相関算出ブロック００７ならびに相関パターン処理ブロック００８における処理を説明するための図である。

相関算出ブロック００７の処理に入る前に、前段の整相加算ブロック００６において、複数の変換素子００２が弾性波を送受信して得られた複数の受信信号を用いて整相加算が行われ、走査線信号が算出される。本実施形態において整相加算ブロック００６は、図３における走査線３０３上と走査線３０４上の複数の走査線信号を算出する。本発明において走査線とは、探触子から送信される弾性波の進行方向に形成される仮想的な線を意味し、走査線信号は、走査線上のある位置における信号を意味する。つまり、１本の走査線は、弾性波の進行方向に走査線信号を並べて形成される仮想線である。

図５は相関算出ブロック００７及び相関パターン処理ブロック００８における処理の流れを示すフローチャートである。まず、ステップ５００において、整相加算手段００６で算出された、走査線３０３上と走査線３０４上の走査線信号が、相関算出ブロック００７に入力される。

ステップ５０１では、相関算出ブロック００７が、走査線３０３上の領域Ｑからの反射波の受信信号（反射波を受信した際に取得される受信信号）を用いた走査線信号Ｘ０、走査線３０４上の領域Ｑ’からの反射波の受信信号を用いた走査線信号Ｙ０を抽出する。つまり、領域Ｑの走査線信号Ｘ０と領域Ｑ’の走査線信号Ｙ０を抽出する。なお、領域Ｑと領域Ｑ’との間は距離Ｌ離れている。その後、複数回の弾性波送受信を繰り返し、領域Ｑの走査線信号Ｘ１、Ｘ２、Ｘ３、Ｘ４、・・・、ＸＮ、領域Ｑ’の走査線信号Ｙ０、Ｙ１、Ｙ２、Ｙ３、Ｙ４、・・・、ＹＮを抽出してゆく（Ｎは正の整数）。

ステップ５０２で、相関算出ブロック００７は、これら走査線信号同士の相互相関値Ｒ（ξ，τ）を、下記式を用いて算出し、出力する。ここでは例えば、走査線信号Ｘ０に対して、他の走査線信号（Ｘ１、Ｘ２、Ｘ３、・・・、ＸＮ、Ｙ０、Ｙ１、Ｙ２・・・、ＹＮ）との相互相関値を夫々算出する。

式中Ｓ（ｘ，ｔ）は、取得位置ｘ、取得時間ｔでの走査線信号である。また、ξ＝ｘ２−ｘ１、τ＝ｔ２−ｔ１である。なお、ｕは相互相関値の算出に用いる走査線信号の抽出位置の調整に用いており、波長に相当する程度の時間幅で最大値の探索を行うことで、より安定的に相互相関値を算出可能である。

ステップ５０３では、相関パターン処理ブロック００８に、相関算出ブロック００７から出力された相互相関値が複数入力される。相関パターン処理ブロック００８では、走査線信号Ｘ０と、その他の走査線信号と、の相互相関値を用いて領域Ｑと領域Ｑ’との間の移動情報を推定する。

以下、相関パターン処理ブロック００８の処理の一例を図４を用いて説明する。図４の横軸は２つの注目位置（領域Ｑ、領域Ｑ’）の距離差を示す。縦軸は走査線信号Ｘ０の取得時間からその他の走査線信号をそれぞれ取得した時間までの時間差（つまり、走査線信号Ｘ０が示す反射波の取得時間とその他の走査線信号が示す反射波の取得時間との差）を示す。点ＲＸ１、ＲＸ２、ＲＸ３、・・・、ＲＸＮ、ＲＹ０、ＲＹ１、ＲＹ２、・・・、ＲＹＮは、この２軸からなる平面上の座標軸に、それぞれ算出された相互相関値をプロットした点である。たとえば、点ＲＹ０は、横軸が領域Ｑと領域Ｑ’との距離（Ｌ）、縦軸が走査線信号Ｘ０を取得した時間と走査線信号Ｙ０を取得した時間との時間差、となる位置に、走査線信号Ｘ０と走査線信号Ｙ０との相互相関値をプロットしたものである。また、ＲＸ１、ＲＸ２、ＲＸ３、・・・、ＲＸＮは、同じ領域Ｑの走査線信号同士の相互相関値をそれぞれ取得時間差毎にプロットしたものである。

次にステップ５０４で、相関パターン処理ブロック００８は、点ＲＹ０から点ＲＹＮまでの間で相互相関値が最大となる点Ａ（Ｌ、τ’）を抽出する。この際、点ＲＹ０から点ＲＹＮまでのそれぞれの点の相互相関値そのものではなく、それぞれの相互相関値を補間してより正確に相互相関値が最大となる点を抽出しても構わない。

ステップ５０５では、点ＲＸ１から点ＲＸＮまでの間で、点Ａと同等の相互相関値を有する点Ｂを抽出する。ステップ５０４と同様に点ＲＸ１から点ＲＸＮまでのそれぞれの相互相関値を補間してより正確に求めても構わない。このようにして求められた、点Ａで接線（点ＲＹ０から点ＲＹ０で形成される）を有し点Ｂを通過する歪んだ楕円形状は、相互相関値の等高線を近似したものとなる。

さらにステップ５０６で、この歪んだ楕円形状上の、縦軸方向において最大値をとる点Ｃ（Ｌ０、τ０）を求める。この点Ｃと原点Ｏとを通る直線の傾きの逆数（Ｌ０／τ０）が、領域Ｑと領域Ｑ’との間を散乱体群が通過する速度となる。

ステップ５０７で、相関パターン処理ブロック００８はこのようにして求めた速度に関するパラメータを移動情報として画像処理ブロック０１０へ出力する。このように相関パターン処理ブロック００８は、走査線信号Ｘ０と、その他の走査線信号と、の相互相関値を複数用いて領域Ｑと領域Ｑ’との間の移動情報を取得することができる。

ここで、散乱体群の時間による変化を考慮せず、単純に領域Ｑと領域Ｑ’とで取得した相互相関値のみを用いて速度を推定する場合を考える。その場合、走査線信号Ｘ０と、領域Ｑ’で取得した走査線信号Ｙ０からＹＮまでと、の相互相関値が最大となった点、つまり点Ａ付近を、散乱体群が通過した時点であると想定することになる。その場合、算出される速度はＬ／τ’となり、真の値（Ｌ０／τ０）とは異なる。つまり、本発明は散乱体群の時間変化まで考慮することで、速度を算出する精度を向上させることが可能であることが分かる。

このように、複数の異なる位置で取得した相互相関値を少なくとも用いて移動情報を算出する場合、時間差の軸と距離差の軸との２軸で表される平面内での相互相関値の分布を取得することで、その算出精度を向上させることが可能となる。つまり、本発明においては、上記平面上での相互相関値の分布を取得することで、ある時間差における、相互相関値の最大値に関する情報を求めることができ、より精度の良い移動情報を求めることができる。

また、ここでは注目位置を領域Ｑと領域Ｑ’の２点としたが、その限りではない（実施形態１−３参照）。さらに複数の点に着目することで２次元平面内、３次元空間内の速度分布等を移動情報として取得することが可能である。なお、走査線方向に関してはドップラ法によって求めた速度を用い、走査線方向ではない方向（例えば走査線方向と垂直な方向）に関しては、本発明の方法を用いて速度を算出することで、２次元平面内、３次元空間内の速度分布を取得することも可能である。つまり、本発明の方法とドップラ法とを組み合わせて速度分布等の移動情報を取得しても良い。
さらに、ここでは領域Ｑと領域Ｑ’との走査線信号を別々のタイミングで取得する例を説明したが、両方の領域に対して同時に弾性波を送信し、両領域を同じタイミングで取得した場合であっても本発明は適用可能である。

以下、図面を用いて本発明による被検体情報取得装置の各実施形態を詳細に説明する。

＜第１の実施形態＞
本発明の第１の実施形態を、図１、図５、図６、図７を用いて説明する。本実施形態では、被検体情報を取得する範囲内において、弾性波のビーム（超音波パルス）を送受信する位置を順に移動させてゆく（つまり走査線信号を取得する位置を順に移動させていく）リニア電子走査等の制御を行った場合に、散乱体群の速度を求める例を述べる。本実施形態では、弾性波ビームを複数位置へ順番に送信する電子走査を２回以上行う。

本実施形態の被検体情報取得装置の装置構成は、図１を用いて説明した構成と同様のため各構成の説明は省略し、相関算出ブロック００７及び相関パターン処理ブロック００８での処理内容について前述とは異なる部分に絞って説明を行う。また、本実施形態の相関算出ブロック００７及び相関パターン処理ブロック００８での処理についても、図５を用いて説明したフローと基本的には同じであり、各ステップにおける前述とは異なる処理内容についてのみ詳細に説明を行う。図６は血管３０１内を流れる赤血球から構成される散乱体群３０２を模式的に示した図であり、図７は相関算出ブロック００７ならびに相関パターン処理ブロック００８における処理の概念を説明するための図である。

（第１の実施形態の処理フロー）
相関算出ブロック００７の処理に入る前に、前段の整相加算ブロック００６において、複数の変換素子００２が弾性波を送受信して得られた複数の受信信号を用いて整相加算が行われ、走査線６０１から走査線６０６上の走査線信号が複数算出される。ここでは、２回以上の電子走査による複数の走査線信号が取得される。そして、ステップ５００で、相関算出ブロック００７にそれらの走査線信号が入力される。

ステップ５０１では、入力された複数の走査線信号のうち、１回目の電子走査による走査線信号として以下の走査線信号が抽出される。なお、それぞれの領域は距離ｐだけ離れている。
・走査線６０１上の領域Ｑ１からの反射波の受信信号を用いた走査線信号Ｓ１＿０
・走査線６０２上の領域Ｑ２からの反射波の受信信号を用いた走査線信号Ｓ２＿０
・走査線６０３上の領域Ｑ３からの反射波の受信信号を用いた走査線信号Ｓ３＿０
・走査線６０４上の領域Ｑ４からの反射波の受信信号を用いた走査線信号Ｓ４＿０
・走査線６０５上の領域Ｑ５からの反射波の受信信号を用いた走査線信号Ｓ５＿０
・走査線６０６上の領域Ｑ６からの反射波の受信信号を用いた走査線信号Ｓ６＿０
その後、送受信を繰り返し、さらに２回目の電子走査による走査線信号として以下の走査線信号が抽出される。
・領域Ｑ１からの反射波の受信信号を用いた走査線信号Ｓ１＿１
・領域Ｑ２からの反射波の受信信号を用いた走査線信号Ｓ２＿１
・領域Ｑ３からの反射波の受信信号を用いた走査線信号Ｓ３＿１
・領域Ｑ４からの反射波の受信信号を用いた走査線信号Ｓ４＿１
・領域Ｑ５からの反射波の受信信号を用いた走査線信号Ｓ５＿１
・領域Ｑ６からの反射波の受信信号を用いた走査線信号Ｓ６＿１
ただし、ここでは、２回目の電子走査による走査線信号を示しているが、３回目の電子走査により得られる走査線信号でもよく、２回目以降であればよい。

ステップ５０２で、相関算出ブロック００７はこれらの走査線信号同士の相互相関値を算出し出力する。ここでは例えば、走査線信号Ｓ１＿０に関して、他の走査線信号（Ｓ１＿１、Ｓ２＿０、Ｓ２＿１、Ｓ３＿０，Ｓ３＿１、Ｓ４＿０，Ｓ４＿１、Ｓ５＿０，Ｓ５＿１、Ｓ６＿０，Ｓ６＿１）との相互相関値を夫々算出する。ここで、走査線信号Ｓ１＿０は、１回目の電子走査による走査線信号のうち所定の基準位置の走査線信号を示すものとする。つまり、本実施形態では、最初の送信位置である領域Ｑ１を基準位置としている。ただし、本発明においては、基準位置は、必ずしも最初の送信位置に限定されず、領域Ｑ１以外の位置（例えば領域Ｑ３等）を基準位置としてもよい。例えば領域Ｑ３を基準位置とした場合は、走査線信号Ｓ３＿０と、その他の走査線信号と、の相互相関値を夫々算出するとよい。

ステップ５０３で、相関算出ブロック００７から出力された相互相関値は、相関パターン処理ブロック００８に入力される。相関パターン処理ブロック００８は、走査線信号Ｓ１＿０（１回目の電子走査における基準位置の走査線信号）とその走査線信号以外の走査線信号（Ｓ１＿０以外の走査線信号）との相互相関値を用いて、領域Ｑ１から領域Ｑ５までの領域での移動情報を推定する。

以下、相関パターン処理ブロック００８の処理の一例を図７を用いて説明する。図７の横軸は領域Ｑ１と各領域Ｑ２、Ｑ３、Ｑ４、Ｑ５、Ｑ６との距離差を示す。縦軸は走査線信号Ｓ１＿０の取得時間とその他の走査線信号の取得時間との時間差（つまり、走査線信号Ｓ１＿０が示す反射波の取得時間とその他の走査線信号が示す反射波の取得時間との差）を示す。図７の点ＲＳ２＿０、ＲＳ３＿０、・・・、ＲＳ６＿０、ＲＳ１＿１、ＲＳ２＿１、・・・、ＲＳ６＿１は、この２軸からなる平面上の座標軸に、それぞれ算出された相互相関値をプロットした点である。

例えば、点ＲＳ２＿１は、横軸が領域Ｑ１と領域Ｑ２との距離ｐ、縦軸が走査線信号Ｓ１＿０と走査線信号Ｓ２＿１との取得時間差、となる位置に、走査線信号Ｓ１＿０と走査線信号Ｓ２＿１との相互相関値をプロットしたものである。

次に、ステップ５０４で、点ＲＳ１＿１から点ＲＳ６＿１までの間で相互相関値が最大となる点Ａ（ｐ’、τ’）を抽出する。この際、それぞれの点の相互相関値そのものではなく、それぞれの相互相関値を補間してより正確に相互相関値が最大となる位置を抽出しても構わない。

ステップ５０５では、点ＲＳ２＿０から点Ｒ６＿０までの間で、点Ａと同等の相互相関値を有する点Ｂを抽出する。ステップ５０４と同様に、それぞれの相互相関値を補間してより正確に点Ａと同等の相互相関値を有する点を求めても構わない。このようにして求められた、点Ａで接線（点ＲＳ１＿１から点Ｒ６＿１で形成される）を有し点Ｂを通過する歪んだ楕円形状は、相互相関値の等高線を近似したものとなる。

さらにステップ５０６で、この歪んだ楕円形状上の、縦軸方向において最大値をとる点Ｃ（ｐ０、τ０）を求める。この点Ｃと原点Ｏとを通る直線の傾きの逆数（ｐ０／τ０）が領域間を散乱体群が通過する速度となる。この直線は、上述したように、各時間差における相互相関値の最大値が繋がれた線を示している。

ここで、散乱体群の時間による変化を考慮せず、単純に取得した相互相関値のみを用いて速度を推定する場合を考える。その場合、走査線信号Ｓ１＿０と、走査線信号ＲＳ１＿１からＲＳ６＿１までとの相互相関値が最大となった点、つまり点Ａ付近で散乱体群が通過した時点であると想定することになる。その場合、算出される速度はｐ’／τ’となり、真の値（ｐ０／τ０）とは異なる。つまり、本実施形態は散乱体群の時間変化まで考慮することで、速度を推定する精度を向上させることが可能であることが分かる。具体的には、時間差の軸と距離差の軸との２軸で表される平面内での相互相関値の分布を取得することで、ある時間差における、相互相関値の最大値に関する情報を求めることができ、より精度の良い移動情報を求めることができる。

最後に、ステップ５０７で相関パターン処理ブロック００８は以上のステップにより求めた速度に関するパラメータを移動情報として画像処理ブロック０１０へ出力する。このように、相関パターン処理ブロック００８は、走査線信号Ｘ０と、その他の走査線信号と、の相互相関値を用いて散乱体群の移動情報を取得することができる。また、相関パターン処理ブロック００８は、速度に関するパラメータだけでなく各位置の速度を示す速度分布として移動情報を取得しても良い。

最終的に求められた速度に関するパラメータや速度分布等の移動情報は、画像処理部０１０で表示用に画像処理され、画像表示系０１１で表示される。図８は表示の一例を示した模式図である。図８では、求められた速度分布を用いて、速度の大きさを計測点からの線分の長さで示し、速度を有する方向（散乱体群が流れる方向）をその線分の方向で示している。図８では線分を用いて表示したが、矢印を用いても良く、また速度の大きさを輝度値、色、太さなどで表現することも可能である。さらにその速度を有する方向を色別で表現しても良い。また、さらに速度分布を用いて、流線で表現しても構わない。本実施形態のように注目位置を横方向に並べて設定するだけでなく、さらに複数の点に着目することで２次元平面内、３次元空間内の速度分布を算出することが可能となる。

また、本実施形態では弾性波を送受信する位置を順に移動させて（つまり、電子走査により走査線を移動させながら）走査線信号を取得しているため、通常のＢ−ｍｏｄｅ像を取得しながら、速度分布を求めることが可能である。なお、走査線を移動させる方向を反転させて取得したデータも用いることによって、さらに精度を向上させることが可能である。

またさらに、散乱体群が移動している方向に対して、走査線を移動させる方向が交差するよう（典型的には逆向き）に設定することで、速度分布の算出精度を向上させることが出来る。よって、被検体内の散乱体群の移動方向を検知した後に、その移動方向に対して逆方向に走査線を移動させる制御を行っても良い。

ここで、図７を用いて上記の方法によって速度分布の算出精度が向上する理由を説明する。散乱体群が移動している方向に対して走査線の移動方向を反対の向きに設定し、図７中の１点鎖線で示した直線上の相互相関値を算出した場合、その１点鎖線において相互相関値が最大となるのは点Ｗであり、点Ｗは歪んだ楕円の接点となる。また、図７中の点Ａは、散乱体群が移動している方向に対して走査線の移動方向を同じ向きに設定し、図７中の点線で示した直線上における相互相関値が最大となる点である（つまり楕円の接点である）。ここで、図７から分かるように、点Ａと点Ｗとの位置で歪んだ楕円の曲率を比較すると点Ｗの方が曲率が大きく、それぞれの点を含む接線上で相互相関値の変化が大きいのは点Ｗを含む接線である。つまり、直線上の相互相関値が最大となる点を求める場合、点Ｗを含む接線の方が、点Ａを含む接戦よりも高い精度で求めることができる可能性が高い。よって、散乱体群の移動方向に対して交差する向き（典型的には逆向き）に走査線を移動させることにより、結果的に速度分布の算出精度が向上することが期待できる。

＜第２の実施形態＞
本発明の第２の実施形態を、図１、図９、図１０を用いて説明する。本実施形態ではランダムな複数の走査線信号の組合せで相互相関値を算出し、さらに精度良く安定に散乱体群の速度を求める例について説明する。

本実施形態の被検体情報取得装置の装置構成は、図１を用いて説明した構成と同様のため各構成の説明は省略し、相関算出ブロック００７及び相関パターン処理ブロック００８での処理内容について前述までとは異なる部分に絞って説明を行う。図９は、相関算出ブロック００７ならびに相関パターン処理ブロック００８における処理を説明するための図である。図１０は相関算出ブロック００７ならびに相関パターン処理ブロック００８における処理フローを示すフローチャートである。

ステップ１０００において、相関算出ブロック００７に整相加算ブロック００６から走査線信号が入力される。ここまでは基本的な処理フロー及び第１の実施形態の処理フローと同様である。

ステップ１００１で、相関算出ブロック００７は、入力された複数の走査線信号から注目位置の走査線信号を抽出する。図９（ａ）は相関算出ブロック００７に入力された走査線信号の中で、抽出する注目位置の走査線信号の取得タイミングと取得位置とを示している。図９（ａ）の横軸は取得位置、縦軸は取得タイミング（取得時刻）である。例えば、点ｆ３で表される走査線信号は取得位置が３Δξで、取得タイミングが２Δτである。

ステップ１００２で、相関算出ブロック００７は、このようにして取得されたＭ種類の走査線信号から２つの走査線信号を抽出し相互相関値を算出する。Ｍは正の整数を示す。Ｍ種類から選択するため、Ｍ×（Ｍ−１）／２の数の相互相関値を算出することが可能である。

ステップ１００３で、相関算出ブロック００７から出力された相互相関値は、相関パターン処理ブロック００８に入力される。

ステップ１００４で、相関パターン処理ブロック００８は、入力された複数の相互相関値を用いて相互相関値の分布を推定する。図９（ｂ）は、図９（ａ）の点ｆ１から点ｆ７の点で表せる走査線信号の組み合わせから算出できる相互相関値がプロットされている様子を示している。以下、いくつか例をあげて説明する。

点ｆ１で表せる走査線信号と点ｆ２で表せる走査線信号との相互相関値は、図９（ｂ）の横軸Δξ、縦軸Δτの位置にプロットされる。点ｆ１で表せる走査線信号と点ｆ３で表せる走査線信号との相互相関値は、図９（ｂ）の横軸３Δξ、縦軸２Δτの位置にプロットされる。また、点ｆ４で表せる走査線信号と点ｆ５で表せる走査線信号との相互相関値は、横軸Δξ、縦軸Δτの位置にプロットされる。点ｆ１で表せる走査線信号と点ｆ３で表せる走査線信号との相互相関値と、点ｆ４で表せる走査線信号と点ｆ６で表せる走査線信号との相互相関値とは図９（ｂ）において、同じ位置にプロットされる。このように、同じ位置にプロットされた相互相関値を平均していくことで、より安定的に相互相関値を算出することができ、速度推定の精度をさらに向上させることが出来る。このように、多数の点における相互相関値を入力として、以下の式を用いて相互相関値の分布を推定する。

ρｉは、距離差ξｉ、時間差τｉの位置での相互相関値を示し、ρｉ’（ξ０，τ０，ｖ）は、距離差ξ０、時間差τ０、速度ｖにおける相互相関値の理論値を示す。Ｎは相互相関値の分布推定に用いる相互相関値の組合せ数（プロットされた位置の種類数）である。例えば図９（ｂ）の場合はＮ＝１３となる。また、αｉは重みであり、例えば

で求められるｄｉを用いて、αｉ＝ｄｉや、αｉ＝ρｉ×ｄｉなどのように設定可能である。時間差の軸と距離差の軸との２軸で表される平面上での原点から遠いほど楕円形状の推定精度が高い、また、相互相関値が高いほど信頼度が高い、と考えられるため、このような重みαｉを用いることでさらに楕円形状の推定精度の向上が見込める。ここでの理論値とは、時間差の軸と距離差の軸との２軸で表される平面上での相互相関値の等高線が、楕円形状で近似できるとして、算出することが可能である。

またさらに、αｉは同じ位置にプロットされた相互相関値の個数もしくはその個数を反映した係数（個数が増加するにつれて増加する係数）に設定しても良い。このように設定した場合、相互相関値を平均することによって精度が向上した相互相関値に、夫々重みをつけて相互相関値の分布推定を行うことができるため、さらに楕円形状の推定精度が向上する。

このように、本実施形態においても、時間差の軸と距離差の軸との２軸で表される平面内での相互相関値の分布を取得することで、ある時間差における、相互相関値の最大値に関する情報を求めることができ、より精度の良い移動情報を求めることができる。ステップ１００５において、相関パターン処理ブロック００８は、以上のステップで求めた速度に関するパラメータを移動情報として画像処理ブロック０１０へ出力する。画像処理ブロック０１０以降の処理は第１の実施の形態と同様であるため、説明を省略する。

なお、図１２（ａ）で示すように、時間差の軸と距離差の軸との２軸で表される平面内において、例えば原点に対して複数の方向に位置する位置（例えば図中の黒点に相当）での相互相関値を算出することで、散乱体群の進行方向によらず安定的に速度を算出することが出来る。また、図１２（ｂ）で示すように、時間差の軸と、ｘ方向の距離差ξｘの軸、ｘ方向に垂直なｙ方向の距離差ξｙの軸、の３軸で表される空間内において、例えば点線で示した半球面上の複数の点（図示しない）のような、原点に対して複数の方向に位置する位置での相互相関値を算出することで、散乱体群の移動情報を２次元平面内で取得することが可能である。

本実施の形態では、複数の走査線信号を組み合わせることで、より多くの相互相関値を算出することができるため、より正確に散乱体群の時間変化を考慮でき、安定して精度良く速度を求めることができる。さらに、ＳＮ比が低い信号を用いた相互相関値や、値が小さい相互相関値を、相関パターン処理ブロック００８における相関値分布の推定に使用しないことで、さらに安定性を高めることができる。

また、速度を求めたい領域付近の走査線（例えば３本、もしくは５本など）上の走査線信号を多く用いることで、さらに精度の良い相互相関値ならびに速度の算出が可能となる。理論的には、横軸を注目位置間の距離差とし、縦軸を走査線信号の取得時間差とした座標空間上の相互相関値分布において、相互相関値を速度で偏微分した時に、その偏微分値が他の領域（座標点群）より大きくなる領域（座標点群）の相互相関値を用いることにより精度の良い速度算出が可能となる。そして、注目位置間の距離差が小さい領域に、大きな偏微分値をもった領域が現れることが多い。よって、観察する領域内において、近接する位置で得られた走査線信号同士を用いることが好ましい。特に、隣接する３本の走査線上の３点の走査線信号を用いることが好ましい。また、速度を求めたい領域付近の走査線をより多く用いるには、先述した重みを、その走査線信号の組合せによって求められた相互相関値に対して大きくしたり、単位時間当たりの送受信回数を他の領域と比較して多くしたりするなどの手法を用いることができる。

＜第３の実施形態＞
次に第３の実施形態について説明する。本実施形態は、血液の存在する領域のみ相互相関値の算出を行うことを特徴とする。図１１は本実施の形態における被検体情報取得装置のシステム概略図を示す。本実施形態の装置構成は、図１の装置構成に加えて血流位置抽出ブロック０１２が加わっているが、それ以外の構成は基本的には図１の装置構成と同じであるため、前述までの実施形態とは異なる部分に絞って説明を行う。

本実施形態では、整相加算ブロック００６から走査線信号が画像生成ブロック００９、相関算出ブロック００７、血管位置抽出ブロック０１２へそれぞれ出力される。画像生成ブロック００９へ走査線信号が入力された後の処理は前述までの実施形態と同様であるためここでは説明を省略する。

血流位置抽出ブロック０１２では、整相加算ブロック００６から入力された複数の走査線信号を用いて、パワードップラ処理を行う。このパワードップラ処理は、被検体内の血液（つまり、複数の赤血球からなる散乱体等の対象物）が存在する対象範囲を特定するための処理である。つまり血流位置抽出ブロック０１２は、本発明における対象範囲特定手段である。

ただし、本発明において、対象物が存在する対象範囲を特定する方法は、パワードップラ処理に限定されず、カラードップラー処理でもよい。また、通常のインピーダンス輝度データ（通常のＢ−ｍｏｄｅ画像等、被検体内の音響インピーダンスの違いを反映した分布の輝度データ）をもとに、血管等の対象範囲を抽出する方法でもよい。この場合は、例えば画像生成ブロック００９と同様に、複数の走査線信号をもとに複数の包絡線信号を算出し、この包絡線信号を用いてインピーダンス輝度データを取得する。そして、この分布から、輝度がある一定以上の位置を抽出する等の処理を行うことで血管壁を抽出し、血管の範囲を対象範囲として取得することができる。さらに、インピーダンス輝度データの画像を画像表示系０１１に表示し、ユーザーにより指定された範囲を対象範囲として取得しても良い。

本実施形態では、対象範囲特定手段である血流位置抽出ブロック０１２は、パワードップラ処理によって抽出された血流が存在する対象範囲の情報をシステム制御部００４へ出力する。

相関算出ブロック００７ではシステム制御部００４から入力された血流が存在する対象範囲の情報に従って、整相加算ブロック００６から入力された走査線信号の内、血流が存在する対象範囲内の走査線信号のみを用いて相互相関値の算出を行う。その後の処理は基本的な処理フロー及び第１の実施形態、第２の実施形態と同様であるため省略する。

本実施の形態によれば、血流が存在する対象範囲のみ相互相関値の算出を行うため、計算の負荷が少なく、より小型の装置を提供することができる。

＜第４の実施形態＞
本実施形態は、同じ時刻に複数の位置の走査線信号を取得し、それらを用いることを特徴とする。本実施形態の被検体情報取得装置の構成は図１や図１１を用いて説明した形態と基本的には同様のため、各構成の異なる機能部分に絞って説明を行う。

本実施形態では、同じ時刻に複数の方向もしくは位置に弾性波を送信し、その反射波を受信する処理を複数回行う。本発明において、「同じ時刻」に弾性波を送信するとは、厳密に同じ時刻の場合だけでなく、同じ時刻に弾性波を送信したとして受信信号を扱える程度の誤差を含んでいる場合も含む。

整相加算ブロック００６は受信信号を用いて整相加算を行って走査線信号を複数算出する。相関算出ブロック００７は、異なる時刻に送信された弾性波に対する（つまり、異なる送信時刻に送信した弾性波に基づく反射波に基づいており）、複数の異なる位置の走査線信号の相互相関値を複数算出する。

相関パターン処理ブロック００８は、複数の相互相関値の中にピークが存在するかどうかを判定し、ピークが存在する場合は、その相互相関値がピークとなる位置を推定する。そして、そのピークとなる時間差と距離差とを基に、移動情報を算出する。また、複数の相互相関値の中に、ピークが存在しない場合は、前述した実施形態のいずれかの方法を用いて移動情報を算出する。なお、相互相関値の中にピークがある場合であっても、先述した実施形態のいずれかの手法を用いて移動情報を算出することも可能である。なお、本実施形態においては、ある領域に弾性波を送信し、その送信した領域の中で複数の方向もしくは位置に関して整相加算を行い、複数の走査線信号を算出することによっても適応可能である。

本実施形態の場合、前述の実施形態のように時間差の軸と距離差の軸との２軸で表される平面内での相互相関値の分布を取得する場合と、相互相関値のピークを用いる場合と、を切り替え、移動情報を求めることができる。

＜第５の実施形態＞
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した各実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。

０００ 被検体
００１ 探触子
００２ 変換素子
００３ 送信回路系
００４ システム制御部
００５ 受信回路系
００６ 整相加算ブロック
００７ 相関算出ブロック
００８ 相関パターン処理ブロック
００９ 画像生成ブロック
０１０ 画像処理ブロック
０１１ 画像表示系
０１２ 血流位置抽出ブロック
３０１ 血管
３０２ 散乱体群
３０３、３０４ 走査線
６０１〜６０６ 走査線

Claims (18)

1. 被検体内の対象物の移動情報を取得する被検体情報取得装置であって、
   被検体に弾性波が送信され、前記被検体内の各位置で反射された反射波を受信した複数の変換素子から出力される複数の受信信号を用いて、前記被検体内の各位置からの反射波に対応する信号を走査線信号として複数取得する走査線信号取得手段と、
   複数の前記走査線信号を用いて、前記対象物の移動情報を取得する処理部を有し、
   前記処理部は、異なる位置の走査線信号同士の相互相関値を複数用いて、時間差の軸と距離差の軸との２軸で表される平面上における相互相関値の等高線に基づいて前記対象物の移動情報を取得することを特徴とする被検体情報取得装置。
2. 前記処理部は、前記相互相関値の等高線の形状を楕円に近似した形状として求め、前記時間差の軸方向のある時間差において前記形状上の相関値が最大値となる点と、原点と、を結ぶ線の傾きの逆数を前記移動情報として取得することを特徴とする請求項１に記載の被検体情報取得装置。
3. 前記処理部は、同じ位置の走査線信号同士の相互相関値をさらに用いて、前記移動情報を取得することを特徴とする請求項１または２に記載の被検体情報取得装置。
4. 前記走査線信号取得手段は、前記弾性波のビームを複数位置に順番に送信する電子走査を２回以上行うことで得られる複数の走査線信号を取得し、
   前記処理部は、
   ２回目以降の電子走査における走査線信号と、１回目の電子走査における基準位置の走査線信号と、の相互相関値と、
   前記１回目の電子走査における基準位置の走査線信号と、該走査線信号以外の走査線信号と、の相互相関値と、
   を用いて前記相互相関値の分布を求めることを特徴とする請求項１または２に記載の被検体情報取得装置。
5. 前記処理部は、
   ランダムな組み合わせの走査線信号同士の相互相関値を用いて、前記相互相関値の分布を求めることを特徴とする請求項１または２に記載の被検体情報取得装置。
6. 前記対象物の移動情報を取得する対象範囲を特定する対象範囲特定手段をさらに有し、前記対象範囲特定手段より特定された対象範囲内の走査線信号のみを用いて前記対象範囲内の前記移動情報を取得することを特徴とする請求項１乃至５のいずれか１項に記載の被検体情報取得装置。
7. 前記処理部は、前記移動情報として前記対象物の速度分布を取得することを特徴とする請求項１乃至６のいずれか１項に記載の被検体情報取得装置。
8. 被検体内の対象物の移動情報を取得する被検体情報取得装置であって、
   被検体に弾性波が送信され、前記被検体内の各位置で反射された反射波を受信した複数の変換素子から出力される複数の受信信号を用いて、前記被検体内の各位置からの反射波に対応する信号を走査線信号として複数取得する走査線信号取得手段と、
   複数の前記走査線信号を用いて、前記対象物の移動情報を取得する処理部と、を有し
   前記処理部は、異なる位置の走査線信号同士の複数の相互相関値のうち、ランダムな組み合わせの走査線信号同士の相互相関値を用いて、前記対象物の時間変化による相互相関値の低下が考慮された前記対象物の移動情報を取得することを特徴とする被検体情報取得装置。
9. 前記処理部は、異なる位置の走査線信号同士の相互相関値を複数用いて、走査線信号同士の取得時間差のうちある時間差における相互相関値の最大値に関する情報を基に、前記対象物の移動情報を取得することにより、前記対象物の移動情報を取得することを特徴とする請求項８に記載の被検体情報取得装置。
10. 前記処理部は、同じ位置の走査線信号同士の相互相関値をさらに用いて、前記移動情報を取得することを特徴とする請求項８又は９に記載の被検体情報取得装置。
11. 前記走査線信号取得手段は、前記弾性波のビームを複数位置に順番に送信する電子走査を２回以上行うことで得られる複数の走査線信号を取得し、
    前記処理部は、
    ２回目以降の電子走査における走査線信号と、１回目の電子走査における基準位置の走査線信号と、の相互相関値と、
    前記１回目の電子走査における基準位置の走査線信号と、該走査線信号以外の走査線信号と、の相互相関値と、
    を用いて前記移動情報を取得することを特徴とする請求項８又は９に記載の被検体情報取得装置。
12. 前記処理部は、複数の位置に同時に弾性波を送信する処理が複数回行われることで得られる、異なる時刻に送信された弾性波に基づいた、且つ、異なる位置の走査線信号同士の相互相関値を用いて前記移動情報を取得することを特徴とする請求項８に記載の被検体情報取得装置。
13. 前記対象物の移動情報を取得する対象範囲を特定する対象範囲特定手段をさらに有し、前記対象範囲特定手段より特定された対象範囲内の走査線信号のみを用いて前記対象範囲内の前記移動情報を取得することを特徴とする請求項８乃至１２のいずれか１項に記載の被検体情報取得装置。
14. 前記処理部は、前記移動情報として前記対象物の速度分布を取得することを特徴とする請求項８乃至１３のいずれか１項に記載の被検体情報取得装置。
15. 前記被検体に弾性波を送信し、前記被検体内の各位置で反射された反射波を受信する変換素子を複数有する送受信器をさらに備えることを特徴とする請求項１乃至１４のいずれか１項に記載の被検体情報取得装置。
16. 被検体内の対象物の移動情報を取得する被検体情報取得方法であって、
    被検体に弾性波が送信され、前記被検体内の各位置で反射された反射波を受信した複数の変換素子から出力される複数の受信信号を用いて、前記被検体内の各位置からの反射波に対応する信号を走査線信号として複数取得するステップと、
    複数の前記走査線信号を用いて、前記対象物の移動情報を取得するステップと、を有し前記移動情報を取得するステップでは、異なる位置の走査線信号同士の相互相関値を複数用いて、時間差の軸と距離差の軸との２軸で表される平面上における相互相関値の等高線に基づいて前記対象物の移動情報を取得することを特徴とする被検体情報取得方法。
17. 被検体内の対象物の移動情報を取得する被検体情報取得方法であって、
    被検体に弾性波が送信され、前記被検体内の各位置で反射された反射波を受信した複数の変換素子から出力される複数の受信信号を用いて、前記被検体内の各位置からの反射波に対応する信号を走査線信号として複数取得するステップと、
    複数の前記走査線信号を用いて、前記対象物の移動情報を取得するステップと、を有し前記移動情報を取得するステップでは、異なる位置の複数の走査線信号同士の相互相関値のうち、ランダムな組み合わせの走査線信号同士の相互相関値を用いて、前記対象物の時間変化による相互相関値の低下が考慮された前記対象物の移動情報を取得することを特徴とする被検体情報取得方法。
18. 請求項１６又は１７に記載の被検体情報取得方法の各ステップをコンピュータに実行させるためのプログラム。

Images