JP6053860B2 - 医療用超音波イメージングシステムにおいて剪断波情報を求めるためのコンピュータプログラム - Google Patents

Description

本発明は、超音波イメージングに関し、とりわけ、超音波イメージングのための剪断波情報を求めることに関する。

剪断波速度検出は、たとえば肝臓疾患の評価等である種々の診断用途で使用される。剪断波速度によって組織の特性を知ることができ、腫瘍や他の領域の検出を行うことができる。しかしこの検出精度は、他の用途で快適に使用するには不十分である。

剪断波検出では、超音波プッシュパルスをスキャンラインに沿って送信し、このプッシュパルスが剪断波を発生させ、組織を変位させる。このようにして発生した組織の変位を検出する。剪断波速度を検出するためには、１つのスキャンラインにおいて複数のプッシュパルスを使用することができる。１つのスキャンラインにおいて複数のプッシュパルスを送信しなくても、複数のプッシュパルスを複数の隣接するスキャンラインに沿って順次送信する。同一位置または複数の隣接する位置において組織が多重に押されることにより、剪断波と組織共振との加算が引き起こされると、幾つかの用途では、アーティファクトが生じるかまたは十分な精度が得られないことがある。

本発明の課題は、上記従来技術の欠点を解消し、より高精度の剪断波推定を実現することである。

医療用超音波イメージングにおいて剪断波情報を求めるための方法の一実施形態を示すフローチャートである。 送信および変位検出のシーケンスを示す図である。 Ａ〜Ｃは、異なる送信シーケンスを示し、Ｄは、Ａ〜Ｃの異なる送信シーケンスが温度に影響したときの作用の一例を示す図である。 アキシャル方向、ラテラル方向および時間的な変位サンプリングの一例を示す図である。 到達範囲と遅延との関係の一例を示す図である。 医療用超音波イメージングにおいて剪断波情報を求めるためのシステムの一実施形態を示す。

まず最初に述べておきたいのは、以下に説明する好適な実施形態には、医療用超音波イメージングにおける剪断波情報を求めるための方法、命令およびシステムが含まれていることである。剪断波の発生源からの距離と各位置においてピーク変位に達するまでの時間とに基づいて各位置ごとに剪断波情報を求めるのではなく、発生源の位置と検出位置との複数の異なる組み合わせから得られる変位プロファイルの相関を求める。送信位置と検出位置との種々の空間的な組み合わせによって得られる変位を使用して、剪断波情報を検出する。複数のラテラル位置の遅延を同時に求めるために、１つの同じ関数において、この相関の組み合わせを使用する。空間的な多様性と、階層的な相関とにより、より高精度の剪断波推定を実現することができる。

第１の態様においては、不揮発性のコンピュータ読み出し可能記憶媒体が、医療用超音波イメージングにおける剪断波情報を求めるためにプログラミングされたプロセッサによって実行可能な命令を表すデータを記憶している。前記記憶媒体は、
第１のスキャンライン上で送信された第１の音波に応答して第１の剪断波によって引き起こされた、時間に依存する組織の第１の変位を検出し、ただし該第１の変位は、該第１のスキャンラインから少なくとも第２のスキャンラインおよび第３のスキャンラインの分だけ離隔している第４のスキャンライン上に位置し、
前記第３のスキャンライン上で送信された第２の音波に応答して第２の剪断波によって引き起こされた、時間に依存する組織の第２の変位を検出し、ただし、該第２の変位は前記第４のスキャンライン上に位置し、
前記第１のスキャンライン上で送信された第３の音波に応答して第３の剪断波によって引き起こされた、時間に依存する組織の第３の変位を検出し、ただし、前記第２の変位は前記第２のスキャンライン上に位置し、
複数の各ラテラル位置ごとに、
・前記第１の変位と前記第２の変位との最大相関に相当する遅延と、
・該第１の変位と前記第３の変位との最大相関に相当する遅延と、
・該第２の変位と該第３の変位との最大相関に相当する遅延と
を計算し、
前記各遅延に依存して、前記複数の各ラテラル位置ごとの組織中の剪断波速度を計算する
ための命令を含む。

第２の側面では、医療用超音波イメージングで剪断波情報を求めるための方法が提供される。複数の異なる位置において送信することによって組織内に発生した第１の剪断波は、１つの同じ位置で検出する。１つの同じ位置において複数送信を行うことによって組織内に発生した第２の剪断波が、複数の異なる位置で検出される。第１の剪断波および第２の剪断波の双方に依存して剪断波速度を計算する。画像は、前記剪断波速度に依存して組織を表現する。

第３の側面では、医療用超音波イメージングで剪断波情報を求めるためのシステムが提供される。受信ビームフォーマが、受信した音波信号に依存して空間的位置を表すデータを出力するように設けられている。前記受信ビームフォーマから出力されたデータに依存して組織変位を推定するためにプロセッサが設けられており、該プロセッサは、前記空間的位置のうち複数の異なる空間的位置から組織変位の少なくとも１対の相関係数を計算し、該相関係数に基づいて伝播時間を求め、該伝播時間に依存して画像を生成する。前記画像を表示するためにディスプレイが設けられている。

本発明は特許請求の範囲によって規定されたものであり、明細書の記載によって、特許請求の範囲に記載の発明が限定されることはない。以下、有利な実施形態を参酌して本発明の別の実施形態および利点を説明する。以下で説明される実施形態は、特許請求の範囲において独立して記載されているか、またはこれらの実施形態の組合せが記載されている。

送信位置と検出位置との種々の空間的な組み合わせによって得られる変位を使用して、剪断波情報を検出する。同一の空間位置において加振（励振）（プッシュパルス）を繰り返すことによって生じる剪断波波形の多義性（不明瞭）を回避するため、インタリーブシーケンスを使用する。このインタリーブにより、このシーケンスがビームの重なりを低減することもでき、温度上昇量を低減することもできる。このようなシーケンスの剪断波信号強度は、同一のスキャンライン上に複数送信するのと比較すると優れている。この剪断波信号強度により、より高信頼性かつロバストな検出を実現することができる。

剪断波のピーク到達時間を検出する代わりに、ビームシーケンスが行われた場所からの剪断波応答の相関を使用して、伝播時間を求める。伝播時間推定問題を、相関係数空間内における最適化問題と見なすことにより、ピーク到達時間の探索やヘルムホルツの式を解くのと比較して最適な解が得られる。

１つの実施形態では、所与の種類の組織とイメージングシステム装置との間の少なくとも１つの距離値を求める。組織の熱作用を低減するため、加振パルスの空間インタリーブ方式を使用する。剪断波検出パルスの空間的位置の分布は分散しており、これらの分散された位置を使用してインターリーブ方式の送信によってデータが取得され、これにより、剪断波の遅延時間または伝播時間の相関による検出が可能になる。

別の実施形態では、剪断波変位の時間的な特性経過の少なくとも１対の相関係数を計算する。各所与の距離間における剪断波の伝播時間は、この相関係数から推定される。伝播時間の逆数を求め、その結果に距離を乗算して剪断波速度を求める。別の実施形態も考えられる。

図１は、医療用超音波イメージングにおいて剪断波を求めるための方法を示す。この方法は、図６のシステムまたは別のシステムによって実現される。付加的なステップ、別のステップを設けるか、またはステップを少なくすることも可能であり、たとえば、ステップ３５を実施せずに、その代わりに別の剪断波特性を求める。別の例として、ステップ３６を実施しない。別の実施形態では、ステップ３５からステップ３０に繋がる矢印によって示された繰り返しを行わない場合がある。またこれらのステップは、記載または図示された順序で実行されるが、別の順序で実行することも可能である。

ステップ３０において、剪断波を生成するために音波パルスを送信する。加振エネルギーが超音波トランスデューサから送信される。この加振は音響エネルギーであり、この音響エネルギーが集束されることにより、各送信ごとに１つまたは複数のビームが得られる。上記の加振は、フェーズドアレイおよび／または機械的な集束部を使用して集束される。この加振を１つの次元において非集束状態とし、例えば高さ方向において非集束状態とすることも可能である。加振エネルギーは、患者の組織内に送信される。この加振が１つの場所で集束されることにより、この加振によって生じた剪断波を検出することができ、たとえば、腫瘍である可能性のある場所を含む周辺の組織位置を検出することができる。

前記複数の送信は、複数の異なる場所において集束される。患者の領域上をスキャンして剪断波を発生させるためには、１つまたは複数の送信を同一位置に集束させる。別の送信は、別の場所において集束される。剪断波は集束領域において発生し、該集束領域からラテラル方向に伝播していく。いかなる送信であっても、剪断波は、集束領域に隣接する場所および／または該集束領域から離隔した場所で検出される。剪断波が組織中を伝播していくにつれて、該剪断波の振幅は低減していく。関心領域において、剪断波に対する組織応答を検出するためには、別の集束領域へ送信を行い、該別の集束領域の辺りで検出を行う。複数の異なる送信集束領域に由来する複数の検出領域は、相互に重なり合う。図２に一例として、１０個のスキャンラインｎ〜ｎ＋１０を示している。設けられるスキャンラインの数を増減することもできる。ここでは、各スキャンラインごとに２つ以上の送信を行って剪断波を発生させる。このことは、４０および４２に示されたプッシュパルスを表す矢印によって示されている。プッシュパルスを送信するためのスキャンラインは、検出用のスキャンラインと同じである。別の実施形態では上記実施形態と異なり、１つ以上のプッシュパルススキャンラインが、受信スキャンラインまたは検出スキャンラインと異なる。

剪断波を発生させるためには、高振幅または高パワーの加振が望ましい。たとえば、各加振の機械的指標は１．９に近く、かつ１．９を超えない。慎重に、かつ探触子のばらつきを考慮して機械的指標を選択するならば、１．７を機械的指標の上限とする場合もある。使用されるパワーを増大させたり（たとえば、１．９を超える機械的指標ＭＩ）、減少させたりすることもできる。同一のスキャンライン上または隣接するスキャンライン上において連続的に送信を行うと、組織の温度が時間とともに上昇してしまう。生物学的な作用には、約４１〜４５℃の高い組織温度、４５℃を上回る温度で生じるプロテイン変質、および５０℃を上回る温度で生じる組織壊死が含まる。組織のスティフネスは、４５℃以下の温度であっても影響を受けることがある。４５℃以上の温度では、粘度および／またはスティフネスが増大することがある。５０℃以上の温度では、組織のスティフネスが高くなり、および／または減衰が大きくなることがある。超音波の送信によって、生理学的な作用が発生することがある。その代わり、２℃を上回る温度上昇を阻止することにより、生理学的な作用が抑えられる。

図３Ａ〜３Ｃに、複数の異なる場所において剪断波を発生させるための３つの異なる送信シーケンスを示す。図３Ａは、複数の隣接するスキャンライン上で順次行われる送信を示す。これらのスキャンラインは、下方向を指す矢印として示されている。送信順序を番号１〜５によって示している。使用されるスキャンラインは増減することもできる。

図３Ｂおよび３Ｃに、時間とともに上昇する温度を下げるための送信シーケンスを示す。図示のように順次行われる送信は、すぐ隣のスキャンライン上で順次送信を行うのを回避するために、スキャンライン全体において空間インタリーブで送信することによって行われる。図３Ｂに、５つのスキャンライン上において加振が発信されるのを示す。この発信は、左から右への順に行われる。この観察面は２つに分かれており、右半分と左半分とに分かれている。最初の第１の送信は、図３Ｂの左半分のうち左側のスキャンライン上で行われ、次の第２の送信は、図３Ｂの右半分のうち左側のスキャンライン上で行われる。この送信パターンは、図３Ｂの右半分のうち右側において最後の送信が行われ、その次に、左半分のうち右側において最後の送信が行われるまで続けられる。この送信シーケンス全体において、時間的に連続する送信の相互間の間隔は、剪断波の検出のための観察面（たとえば関心領域）の半分になっている。

図３Ｃに、別の可能なシーケンスを示す。この実施例では、観察面は４等分されており、各４等分ごとに、インタリーブ周期中に順次スキャンが行われる。異なる分割形態の別のスキャン順序、スキャン順序の無作為化、時間的／空間的パターンが異なる別のスキャン順序を使用することもできる。

図３Ｄに、図３Ａ〜３Ｃの３つの異なるシーケンスで発生した温度上昇を示す。組織温度上昇は、ビームの音場が相互に著しく重なり合う近接した場で見られる。図３Ａのシーケンス１（破線）では、時間とともに大きな温度上昇が生じている。図３Ｂのシーケンス２（一点鎖線）および図３Ｃのシーケンス３（実線）では、温度上昇が小さくなっている。

上述の実施形態の代わりに、または付加的に、剪断波推定において付加的な相関を求めて剪断波速度を求めることにより、温度を制御することができる。この付加的な相関を求めて、複数のラテラル位置における遅延を同時に求めることができるようにするためには、１つの同じ場所に対して限られた数のプッシュパルスを有する取得シーケンスを段階的に行い、たとえば２つ、３つまたは４つのプッシュパルスを有する取得シーケンスを段階的に行う。図２には、符号４０および４２によって、剪断波を発生させるためのプッシュパルスまたは送信の２つの段階を示す。短時間のパルス繰り返し周期（ＰＲＩ）で１つの同じ位置において複数のプッシュ時間（たとえば５個または１０個以上）を使用する代わりに、２つ以上の取得段階が使用される。各段階４０，４２がスキャンライン全体を任意のパターンでスキャンし、たとえば図３Ａ〜３Ｃに示されたパターンでスキャンする。２つの段階を使用するだけで、１つの空間的場所における加振ＰＲＩは非常に小さくなる（１００ｍｓのオーダ）。

別の送信も行われるが、この別の送信は、剪断波を生じさせるのではなく剪断波を検出するためのものである。検出のための送信のパワーをより低くし、かつ／または検出用の送信のパルス時間を短時間にすることができ、検出用の送信では、剪断波を発生させるためのプッシュパルスと同じスキャンラインを使用するか、または異なるスキャンラインを使用することができる。複数のスキャンライン上で受信サンプルを形成するために、検出用の送信ビームの少なくとも一方向における形状、たとえばラテラル方向の形状を、より広幅にすることができる。

ステップ３２において、上記の組織の変位が検出される。剪断波によって組織が動く。図２に、剪断波をモニタリングするために使用される変位を示す。これは、符号４４の矢印によって示されている。これらの矢印は１方向であるが、剪断波は複数の方向に伝播する。１方向、２方向またはそれ以上の数の方向において、剪断波をモニタリングすることができる。

剪断波によって引き起こされた組織の変位は、ある程度の時間にわたって検出される。剪断波が発生した場合、この剪断波の応力が伝播するまで、隣接する組織は動かない。剪断波が或る場所を通ると、組織はある程度の量または距離だけ変位する。この変位量または変位距離は、ピーク変位量に達するまで増大し、組織が静止状態に戻っていくにつれて変位量または変位距離は減少していく。このピーク変位量は、集束点における剪断波の発生時点から該剪断波のピークが別の場所を通るまでの遅延時間または伝播時間に相当する。ある程度の時間にわたる変位経過は、ある程度の時間にわたって組織の動きの大きさをトラッキングすることによって求められる。

組織の変位は、超音波スキャンによって検出される。この変位を検出するためには、変位中の組織に超音波エネルギーが送信され、このエネルギーの反射が受信される。図２にこの検出が４６で示されている。検出用の送信および受信は複数回行われ、変位に起因する変化が求められる。送信および受信には任意のシーケンスを使用することができる。変位の検出は別のスキャンによってインタリーブすることができ、たとえば、異なる変位領域を別個にスキャンすることができる。

関心領域をモニタリングしてこの剪断波を検出する。この関心領域のサイズは、例えばラテラル方向に6mm，アキシャル方向に10mmなど任意である。この検出領域は超音波によってモニタリングされ、例えば、Ｂモードスキャンを実行して、上記の剪断波によって引き起こされた組織の変位を検出する。ドップラモード、カラーフローモード、または別の超音波モードを使用して上記の剪断波をモニタリングすることができる。

このモニタリングは、任意の数のスキャンラインで実行され、例えば、各送信ごとに４つの受信ビームが応答して形成される。加振エネルギーを送信して剪断波を発生させた後、１つのスキャンライン上でＢモード送信が繰り返し行われ、また、４つの隣接するスキャンライン上で受信が行われる。別の実施形態では、各送信ごとに１つの受信ビームだけが応答して形成されるか、または別の個数の受信ビームが形成される。任意の回数の繰り返し、例えば、約１０〜１２０回の繰り返しを使用することができる。上記超音波データの中には、例えば、上記の繰り返しの開始時または終了時において生じる超音波データの中には、上記剪断波に応答して生じたものではない超音波データがある。

上記剪断波はスキャンラインを通って伝播するため、Ｂモード強度は変化することがある。上記モニタリングされるスキャンラインでは、剪断波によって発生した組織の動きの時間特性を表すデータ系列が得られ、例えば、複数の空間位置から得られたデータ（例えば、複数のスキャンライン上で得られたデータ）は時間に依存して、データまたは画像の参照フレームに相関付けされる。任意の弾性検出を使用することができる。複数の奥行位置または空間位置にわたって、各奥行位置または空間位置毎に相関が行われる（例えば、プロファイル計算が行われる対象であるポイントである奥行の中心を有する６４個の奥行のカーネル）。空間における２次元または３次元の変位を使用することができ、また、スキャンラインまたはビームと異なる方向の１次元の変位も使用することができる。

特定の時点における最大の相関または十分な相関を有する空間的なずれが、変位の量を示す。１つの位置において複数の異なる時点の変位が求められ、各位置ごとの時間特性により、剪断波の検出結果が分かる。この時間特性にノイズが無いか否か、または、この時間特性の変化は単独のものであるか否かが検査される。時間的なローパスフィルタリングの有無にかかわらず、時間特性におけるピークによって、剪断波の波面が通り過ぎたことが分かる。最大変位が選択されるが、平均または別の変位統計結果を使用することができる。図２に、４８において検出される変位の時間的応答を示す。

より大きな領域をモニタリングするためには、１つの同じプッシュパルスに応じて形成される受信ビームを追加する。その代わりに、さらに別の剪断波を発生させ、送信ビームが送信される場所からこの剪断波の発生地点までの距離と、受信ビームが送信される場所からこの剪断波発生地点までの距離とが異なるようにする。上記の実施形態を使用すると、スキャンラインあたり２つ以上のプッシュパルスを使用するだけで、複数のラテラル位置の変位が得られる。受信ビーム位置毎に、超音波データによって表される動き情報の時間特性が得られる。時間特性を形成する間、より高い時間分解能を実現するためには、１つの同じ剪断波をモニタリングするために複数の異なるスキャンライン上で送信を行うのを避け、スキャン位置をインタリーブまたはシフトすることができる。

上記の説明は、奥行が１つである場合についての説明である。上記の関心領域のアキシャル方向の拡がりをすべてカバーする１つのゲートが得られるように、上記のサンプリングを配置することができる。別の１実施形態では、各受信ビーム毎に複数の奥行においてサンプルが取得される。アキシャル方向の奥行でも、ラテラル位置でも、各位置ごとに別個の時間特性が得られる。任意の数の奥行を使用することができ、たとえば５ｍｍで約２００個の奥行を使用し、１０ｍｍで約４００個の奥行を使用することができる。

上記の関心領域において、異なる位置を表す超音波データが得られる。この超音波データは、スキャンによってリアルタイムで得られるか、またはメモリから得られる。各位置ごとに、動き情報は複数の異なる時点における応答を表す。別のスキャン手法、別のモニタリング手法、または別の技術を使用して、超音波データを得て、剪断波の大きさを推定することができる。

時系列でラテラル位置に配列した複数の変位を使用することにより、スキャンラインごとに発信される２つの少ないプッシュパルスで、剪断波速度を推定することができる。たとえば、複数の異なる位置からの送信によって、組織中の剪断波が１つの同じ位置で引き起こされて検出され、１つの同じ位置で複数の送信が行われることにより、複数の異なる位置において組織中の剪断波が引き起こされて検出され、別の複数の位置において送信が行われることにより、組織中の剪断波が複数の異なる位置において引き起こされて検出される。このように空間的および時間的に多様であることにより、多義性を小さくして剪断波を求めることができる。

図２に、このような種々の位置を示す。このパルスパターンでは、剪断波を発生させるために各指数の位置ごと（たとえば各スキャンラインごと）に加振力が２回印加され、検出パルス列もまた、別の位置で剪断波を検出するために２回印加される。しかし、各検出パルス列は複数の位置に対応するので、解を検出する剪断波の３つの可能なセットを有するシステムが得られる。そのうち１つのセットは、同じ１つの位置で検出されるが２つの異なる位置で加振されるシーケンスであり、次の第２のセットは、複数の異なる位置で検出されるが１つの同じ位置で加振されるシーケンスであり、最後の第３のセットでは、検出も加振も複数の異なる位置で行われる。図２にこの基本的構成が、４６における指数オフセットパラメータｐおよびｑと、２つの管形を結ぶ接続部分とによって示されている。時間的および空間的な多様性を加減することができ、たとえば、位置および／または時点を追加し、これらの組み合わせを追加することができる。

図２中、４８において検出された変位応答に描かれた楕円によって、別の成分を含む領域を示している。この変位の時間的応答は、この別の成分を含む領域の内部の位置と外部の位置とで異なる。媒質が異種を含むことにより、上記の解法のどのセットも完璧でなくなる。

図４に、変位の時間的応答の空間的位置および時間的位置と、この変位に相当する加振の一例を示す。Ａ，ＢおよびＣは剪断波であり、これらの剪断波は、それぞれに対応する加振位置（破線）で発生したものである。Ａ′，Ｂ′およびＣ′は、検出位置で検出された剪断波である。これらの検出位置は、相互に等間隔で離隔されたスキャンライン上にあるが、これらの検出位置を別の分布にすることもできる。上記の例では、剪断波Ａ，ＢおよびＣを発生させるための３つのプッシュパルス、および、各剪断波に対応して検出され時間に依存する組織変位（Ａ′，Ｂ′およびＣ′）である。図４に、１つの奥行におけるプロセスを示す。横軸はラテラル方向の間隔を示しており、ここでは、剪断波ＡおよびＣを発生させるためのプッシュパルスが第１のスキャンライン上で送信され、剪断波Ｂを発生させるためのプッシュパルスが第３のスキャンライン上で送信される。別の横軸が時間を表しており、複数の異なる時点で第１のスキャンライン上で送信され剪断波ＡおよびＣを発生させるためのプッシュパルスを示している。たとえば、剪断波ＡおよびＣを発生させるためのプッシュパルスは、図２中の１つの同じスキャンラインｎ＋３上のプッシュパルス４０，４２に相当する。剪断波Ｂを発生させるためのプッシュパルスが同時点で送信されるのが示されているが、このことは比較的そうなっているということであり、実際の送信は、第３のスキャンライン上（たとえば、図２中のスキャンラインｎ＋６）で順次行われる。同じ第４のスキャンライン上で検出される変位Ａ′およびＢ′の関係をより良好に示すため、プッシュパルスＡおよびＢに応答して行われる変位検出のタイミングが同時点で行われるように示している。

図２および４の例では、スキャンラインｎ＋３上の音波送信に応答して剪断波Ａが発生することによって生じた変位Ａ′が検出される。変位Ａ′はスキャンラインｎ＋６上に位置し、図２ではこの変位Ａ′を、第２のプッシュパルスからの矢印によって示す。この検出は、スキャンラインを２つ介在させて行われる。他に介在するスキャンラインは無い。他の実施形態では、上記実施形態の代わりに、介在するスキャンラインを１つとするか、または２つより多くする。スキャンラインｎ＋５上の音波送信に応答して剪断波Ｂが発生することによって引き起こされた変位Ｂ′が検出される。変位Ｂ′はスキャンラインｎ＋６上で発生する。このスキャンラインｎ＋６は、剪断波Ｂを発生させるためのプッシュパルスで使用されたスキャンラインｎ＋５の直ぐ隣である。介在するスキャンラインを追加することができる。剪断波Ａを発生させるためのプッシュパルスと同じスキャンラインｎ＋３上で送信された音波に応答して剪断波Ｃが発生することによって引き起こされた変位Ｃ′が検出される。この変位Ｃ′は、直ぐ隣のスキャンラインｎ＋４上で検出されるが、介在するスキャンラインを１つとするか、またはより多くすることができる。図２中、ＣおよびＣ′を表すために、スキャンラインｎ＋３上のプッシュパルス４０からスキャンラインｎ＋４上の検出パルス４６の方向を指す矢印が示されている。

別の空間的関係および時間的関係を設けることもできる。図２および４は一例であり、スキャンラインｎ＋３および／またはｎ＋５上においてさらに別のプッシュパルスを使用することもできる。１つの同じスキャンラインまたは異なるスキャンライン上で、プッシュパルス送信を追加するかまたは追加せずに、さらなる変位を検出することもできる（たとえば、同じプッシュパルスＣに応答してｎ＋４上でＣ′の他に発生する別のＣ″も検出する）。図２および４には、送信の右側から離れたスキャンライン上で検出が行われるのが示されているが、どのような距離に配置してもよい（たとえば、２次元における左側から離れたスキャンライン、または、３次元における前面／背面から離れたスキャンライン）。

図４では、ＡＡ′およびＢＢ′が、範囲＿１における伝播時間パラメータを検出するための対を成し、ＡＡ′およびＣＣ′がそれぞれ、範囲＿２における伝播時間パラメータを検出するための別の一対を成す。ＢＢ′およびＣＣ′は、解を限定するための別の一対を成す。前記範囲パラメータは、スキャンライン相互間の既知の距離を表す。範囲＿１および範囲＿２を伝播したときの伝播時間を、前記空間的／時間的な対から遅延１および遅延２として求めることができる。剪断波の速度または伝播時間を求めるために、この遅延情報を使用することができる。

ステップ３４において、空間的な相違を利用して遅延を計算する。変位を相関づけて、遅延を計算する。発生する可能性のある様々な遅延の範囲を考慮することにより、変位のプロファイルを使用して、同じ計算で複数の異なる位置に関連する遅延を求めることができる。１つの位置については、空間的に異なる複数の変位の最大相関に対応する遅延が計算される。１つの相関を行うというよりも、この最大相関は複数の相関に依存し、たとえば複数の異なる変位対の相関に依存する。剪断波速度または別の剪断波情報は、相関関係を利用して、検出された剪断波との関係に基づいて計算される。たとえばスペックルトラッキングを行うために、相関関係を利用して変位を求める場合、さらなる相関または相関関係の階層化を使用して、変位から剪断波情報を求める。

通常の波伝播時間は、以下の関数によって求められる：

上記式では、τは遅延であり、ｃは相関性であり、ｔは時点であり、ｄは変位であり、ｍ，ｎ，ｕおよびｖは、２次元における空間的座標である。たとえば、これら２つの変位関数は変位曲線または変位プロファイルＡ′およびＢ′である。この最大相関は、十分な大きさの信号雑音比を実現し、信号波形の歪みを小さくするのに十分であるが、別の状況では十分でない場合もある。

図２において示された異質混合状態の組織領域のような、移行区域で生じる問題を解消するために、空間的な相違が使用される。複数の相関対を使用して遅延を求める。空間的相違および時間的相違を使用して、複数の各ラテラル位置ごとに、剪断波の空間的分布と、剪断波の持続時間とを求めることができる。最大相関を見つけるために、相互相関の群の挙動を使用する。各ラテラル位置ｙごとに、可能な遅延範囲を使用して、最大相関に関連づけられる遅延を求める。同じ奥行におけるラテラル方向の拡がり全体の相関関数を、２次元行列ｃ（ｙ，τ）として得ることができる。複数の各ラテラル位置における遅延が同時に求められる。複数のスキャンラインにおける同時点の剪断波速度に相当する遅延を、同じ関数の一部として求めることができる。剪断波変位時間特性の少なくとも１つの対で、複数の異なる位置の変位の相関が計算される。対応する位置の遅延の行列を出力するために１つの関数が使用される。

変位の時間的応答はｄ（ｉ，ｊ，ｔ）で表され、指数ｉは加振位置指数（たとえばスキャンライン）を表し、ｊはトラッキング位置指数（たとえばスキャンライン）を表す。変位曲線の２対間のパルス遅延ｄ（ｉ，ｊ，ｔ）を検出するために相関を使用する。すなわち、｛ｄ（ｉ，ｉ＋ｐ，ｔ），ｄ（ｉ，ｉ＋ｑ），ｔ）｝、｛ｄ（ｉ，ｉ＋ｐ，ｔ），ｄ（ｉ＋ｐ−ｑ，ｉ＋ｐ，ｔ）｝である。異質混在区域を特定するためには、別の対｛ｄ（ｉ，ｉ＋ｑ），ｔ），ｄ（ｉ＋ｐ−ｑ，ｉ＋ｐ，ｔ）｝が使用される。図４の例では、一対はＡ′Ｂ′であり、別の一対はＡ′Ｃ′であり、異質混在区域を特定するために使用される対はＢ′Ｃ′である。

たとえば、行われるプッシュパルス送信および／または検出を追加する場合には、別の対を使用することができる。空間的に相違する変位対であるか、または時間的に相違する変位対であるか、または空間的および時間的に相違する変位対であれば、どのような変位対でも使用することができる。この計算は、スキャンライン相互間の距離との関係に基づいて行われ、この距離はたとえば、１つの同じスキャンライン上のプッシュパルスに対して複数の異なるスキャンラインにおいて生じた変位を検出することによって得られる。この到達範囲を決定するために、スキャンライン間の範囲が使用される。この到達範囲と遅延とを使用して、剪断波の速度が求められる。

一例の関数は、３つの異なる変位対の相関度の和を計算する。１つの同じスキャンライン上の複数の送信によって発生して検出されかつ／または複数の異なる検出スキャンラインで検出された剪断波の相互間の相関が求められる。同じ奥行においてτが連続的に変化し続けると仮定すると、同じ奥行におけるラテラル方向の拡がり全体で得られるτ値は境界線または境界曲線Гτと称され、以下のように定義された費用関数を解くことによって得られる：

上記式では、

この曲線は遅延τおよびラテラル位置ｙの行列として表され、τ（ｙ）は、ｙにおけるτの１階微分であるか、または２階微分であるか、または３階微分であるか、または別の微分である。α、β、γは重み付けパラメータである。これらの重み付けパラメータは、実数の１．０とするか、または別の数にするか、またはτに依存することができる。任意の範囲の遅延値および任意の数の遅延値を使用することができ、たとえば、考察されている０〜１００の単位範囲における可能なτ値は３２個である。使用可能な剪断波伝播から得られたすべての情報を組み合わせることにより、各奥行ごとに最大相関係数の曲線が求められる。上記実施形態の代わりの別の実施形態では、使用される対の項または相関を追加することができる。異なる多様性の組み合わせを使用することができる。

数式（２）を使用するか、または、複数の相違対の項を含む別の相関関数を使用することにより、弾性率が異なる２つの位置における加振に起因するばらつきを低減することができ、相関項ｃ_３を導入して（たとえばＢ′Ｃ′の相関）媒質の異質混合状態に適応的に対応することができ、隣接の制約を適用することによって誤差を低減することができる。

図４に示されているように、遅延τの推定は、相関関数を計算するのに使用されるデータ対では異なる。対データ｛ｄ（ｉ，ｉ＋ｐ，ｔ），ｄ（ｉ，ｉ＋ｑ），ｔ）｝（Ａ′Ｂ′相関）が範囲＿１における平均遅延時間を推定し、対｛ｄ（ｉ，ｉ＋ｐ，ｔ），ｄ（ｉ＋ｐ−ｑ，ｉ＋ｐ，ｔ）｝（Ａ′Ｃ′相関）が範囲＿２における平均遅延時間を推定する。これら２つの推定結果には、共有される範囲がある。複数の隣接するラインの各群の間の高精度の推定結果（τ）は、両端の値が既知である場合、ラテラル方向の拡がり全体においてΓτをアンラッピングすることによって得られる。このラテラル方向全体の拡がりは関心領域全体（たとえば１０個のスキャンライン）に及ぶものとすることができ、また、関心領域のサブセットに及ぶものとすることもできる。前記既知の値は、両端に対して僅かに異なる、加振パルスおよび検出パルスのシーケンスによって得られる。撮像視野の両端では、検出Ｃ′はＣと同じ位置にあり、加振ＢはトラッキングパルスＢ′と同じ位置にある。このことにより、推定伝播時間はＡからＡ′までの平均値となる。各ビーム群間の間隔のτ値を解くために、以下の数式が使用される：

両端におけるシーケンスにより、τ_０およびτ_Ｎ−１の推定結果が得られる。

図５に、数式（４）の中間部分を示す。両端の時点τ_０およびτ_Ｎ−１において値が既知である場合の遅延τ値をアンラップするための方法は、再帰的解法方程式と称される。与えられたτ_０の低い方から高い方へ順に、以下の式が適用される：

ａ（ｎ）もｂ（ｎ）もすべて計算することができるので、数式（４）は２つの未知の変数にまで縮小される。この方程式系と、τ_Ｎ−１が既知であるということとから、指数の高い方から低い方へ順に数式（５）を解いて、複数の隣接するビームの各群ごとにτを得ることができる。平行なビームを受け取る場合、同じ方法を使用して個々のτに分解する。たとえば、指数０〜４が付された４つの平行なビームをトラッキングパルスに使用する場合、数式（４）ないし（６）によって、この指数を有する各ビームを求めることができる。指数が付された各ビーム間の距離が異なることを考慮して、数式（４）において重み付けされる係数１／２を変化させている。

ステップ３４において、遅延およびラテラル位置に依存して剪断波速度を計算する。組織の複数の異なる空間的位置の剪断波速度が求められる。複数の各ラテラル位置ごとの組織中の剪断波速度は、各遅延に依存して計算される。

剪断波速度は、伝播時間または遅延（τ）と、スキャンライン間の既知の距離（たとえば範囲）とから求められる。位置までの時間および距離が速度を決定する。この距離はスキャンライン間隔（すなわち、剪断波を発生させるための送信ビーム位置と、該剪断波を検出するための受信ビーム位置との間隔）から分かる。

剪断波速度の代わりに、または剪断波速度の他に付加的に、弾性率または剪断波に関する任意の値を推定することができる。組織弾性率値は、その位置における硬さまたはスティフネスを表す。例えば、組織の剪断弾性率が推定される。その代わりに、他の実施形態では、ヤング弾性率が推定される。別の実施形態では、剪断波に関する別の定量的または定性的な値が推定される。

上記の剪断弾性率は、ｇ＝ρｖ_s ²によって得られ、ここでρは密度、ｖ_sは剪断波速度の推定値である。１実施形態では、上記の組織弾性率または剪断波情報、例えば剪断弾性率は、ひずみまたは変位と、上記の弾性率または剪断波情報とに依存して求められる。例えば、各サンプル位置の剪断弾性率が、拡散方程式を繰り返して解くことによって求められる。0.5のポワソン比を仮定するかまたは既知のポワソン比を使用することにより、特定の剪断弾性率の場合に複数の異なる時点で発生した歪み場に依存して、または、特定の剪断弾性率の場合に複数の異なる位置において異なる応力がかかった場合の歪み場に依存して、複数の異なる位置における剪断弾性率が繰り返して計算される。

別のスキャンラインおよび／または別の奥行において上記のステップを繰り返す。たとえば、１つまたは複数の各奥行において、ステップ３０，３２，３４および３５を再び行う。これらのステップは、ラテラル方向に離隔された複数の別の位置に対して繰り返すことができる。図４は一例として、シーケンス全体のうち一部を示す。数式（３）を解く際に、複数の異なるラテラル位置ごとに異なる複数の相関が準備される。その代わりに、別の各ラテラル部分を別個に扱い、たとえば４つのスキャンラインの一意の群で扱う。その代わりに、関心領域の別の解法または別の部分に、プッシュパルスと、検出された変位対または変位情報を使用することができる。たとえば１つの解法では、ＢおよびＢ′をＣおよびＣ′の対として使用することができる。

ステップ３６において、画像が形成される。この画像は、前記剪断波速度または別の剪断波情報に依存して組織を表現する。剪断波情報を求める対象である領域は、患者の比較的大きな面積または体積を表す画像中のサブセットまたは関心領域とすることができる。たとえば剪断波速度は、一領域内の画素の色を変調して濃淡Ｂモード画像を生成する。この画像は、複数の異なる位置における剪断波または弾性率（例えば剪断弾性率）などの変位情報を表現することができる。ディスプレイグリッドは、スキャングリッドおよび／または変位が計算されるグリッドとは異なるものとすることも可能である。色、明るさ、輝度、色相または別の特性を、剪断波情報に依存して変調する。

図６に、医療用超音波イメージングにおいて剪断波情報を求めるためのシステム１０の一実施形態を示す。システム１０により、図１，図２，図４の方法または別の方法が実現される。システム１０には、送信ビームフォーマ１２と、トランスデューサ１４と、受信ビームフォーマ１６と、画像処理プロセッサ１８と、ディスプレイ２０と、記憶装置２２とが含まれている。別の構成要素を付け加えてもよいし、上記の構成要素と異なる構成要素を設けることができ、また設けられる構成要素の数を低減してもよい。たとえば、剪断波情報を得る対象である関心領域を手動で指定するか、または該関心領域の指定を支援するために、ユーザ入力部が設けられる。別の例として、組織を治療するために付加的にＨＩＦＵトランスデューサ２４を設ける。システム１０は、医療用超音波診断画像システムである。別の実施形態ではその代わりに、システム１０は、パーソナルコンピュータ、ワークステーション、ＰＡＣＳステーションであるか、または、リアルタイム画像構築または取得後画像構築を行うために同じ場所に配置されるかまたはネットワーク上で分散されて配置される別の構成の装置である。

送信ビームフォーマ１２は、超音波送信器、記憶装置、パルス発生器、アナログ回路、デジタル回路またはこれらの組み合わせである。送信ビームフォーマ１２は複数のチャネルに対して、異なる振幅または相対的な振幅、遅延および／または位相調整で波形を形成する。生成された音波に応答してトランスデューサ１４から音波が送信されると、１つまたは複数のビームが形成される。複数の送信ビームが順次形成されて２次元または３次元領域がスキャンされる。セクタ、Ｖｅｃｔｏｒ（登録商標）、線形スキャンフォーマット、または別のスキャンフォーマットを使用することができる。１つの領域が複数回スキャンされる。フローイメージングまたはドップラーイメージングにおいて、また歪みイメージングにおいて、複数のスキャンを順次行う。ドップラーイメージングでは、上述の順次行われる複数のスキャンには、隣のスキャンラインをスキャンする前に同じスキャンラインに沿って生成される複数のビームを含むことができる。剪断波イメージングには任意のスキャンシーケンスを使用することができ、たとえば図３Ａ〜３Ｃに示されたようなスキャンシーケンスを使用することができる。別のスキャンラインを走査する前に形成される、剪断波発生用の送信ビームを、１つまたは２つまたはそれ以上とすることができ、また、送信ビームをインタリーブすることもできる。たとえば取得の順序は、相関の計算と別個である。記憶装置が、相関を行うために変位情報を得るためのデータを記憶する。別の実施形態ではその代わりに、より高速のスキャンを実現するために送信ビームフォーマ１２が平面波または発散波を生成する。

トランスデューサ１４は、圧電ダイアフラム素子または容量性ダイアフラム素子の１次元、１．２５次元、１．５次元、１．７５次元または２次元のアレイである。トランスデューサ１４には、音響エネルギーを電気エネルギーに変換し電気エネルギーを音響エネルギーに変換するための複数の素子が含まれている。たとえばトランスデューサ１４は、ＰＺＴ素子の２次元アレイである（たとえば３，０００個のＰＺＴ素子）。受信信号は、トランスデューサ１４の素子に衝突する超音波エネルギー（エコー）に応答して生成される。これらの素子は、送信ビームフォーマ１２および受信ビームフォーマ１６のチャネルに接続される。

受信ビームフォーマ１６には、増幅器、遅延要素および／または位相回転器および１つまたは複数の加算器を有する複数のチャネルが含まれている。各チャネルは、１つまたは複数のトランスデューサ素子に接続される。受信ビームフォーマ１６は、相対遅延量、位相および／またはアポダイゼーション量に基づき、検出のための各送信ごとに応じて１つまたは複数の受信ビームを形成する。受信ビームフォーマ１６は、上記の受信音響信号を使用して空間位置を表すデータを出力する。複数の異なる素子からの信号の相対的な遅延および／または位相調整および加算により、ビーム成形が行われる。別の実施形態ではその代わりに、受信ビームフォーマ１６は、フーリエ変換または別の変換を用いてサンプルを形成するプロセッサである。

受信ビームフォーマ１６は、フィルタを含むことができ、送信周波数帯域に対する２次高調波または別の周波数帯域における情報を分離するフィルタなどを含むことができる。このような周波数の情報の方が、所望の組織情報、造影剤情報および／または流れ情報を含んでいる確率が高い。別の実施形態では受信ビームフォーマ１６にメモリまたはバッファと、フィルタまたは加算器とが含まれる。所望の周波数帯域の情報を分離するためには、たとえば２次高調波、３次基本波または別の帯域を分離するためには、２つ以上の受信ビームを組み合わせる。

受信ビームフォーマ１６は、複数の空間位置を表すビームを加算したデータを出力する。１つの位置、線に沿った複数の位置、面積における複数の位置または体積における複数の位置に対応するデータが出力される。動的な集束を行うことも可能である。これらのデータの目的は異なっていてもよい。例えば、Ｂモードデータや組織データを得るためには、剪断波検出を行う場合と異なるスキャンが行われる。択一的には上記のＢモードデータは、剪断波によって引き起こされた変位を求めるためにも使用される。

プロセッサ１８は、Ｂモード検出器、ドップラ検出器、パルス波ドップラ検出器、相関プロセッサ、フーリエ変換プロセッサ、特定用途向け集積回路、汎用プロセッサ、制御プロセッサ、画像処理プロセッサ、フィールドプログラマブルゲートアレイ、デジタル信号処理プロセッサ、アナログ回路、デジタル回路またはこれらの組み合わせ、あるいはビーム成形された超音波サンプルからディスプレイ用の情報を検出および処理して表示する現在公知の別の装置または将来開発される装置である。

１実施形態では、プロセッサ１８には１つまたは複数の検出器と、別体のプロセッサとが含まれている。この別体のプロセッサは、制御プロセッサ、汎用プロセッサ、デジタル信号処理プロセッサ、特定用途向け集積回路、フィールドプログラマブルゲートアレイ、ネットワーク、サーバ、プロセッサ群、データパス、それらの組み合わせ、または変位を検出して組織特性を計算するための現在公知の別の装置または将来開発される装置である。例えば、上記の別体のプロセッサにより、図１に示したステップのうちの１つまたは複数の任意の組み合わせが実行される。プロセッサ１８は、上記ステップを実施するためのソフトウェアおよび／またはハードウェアによって構築される。

１つの実施形態ではプロセッサ１８は、受信ビームフォーマ１６からの出力データに依存して組織変位を推定する。この組織変位は、時間に依存する変位の大きさのプロファイルとして推定されるか、または、時間に依存する変位の大きさの曲線を表すデータとして推定される。前記変位プロファイルは、参照データと、剪断波が通過する間に得られたデータとの相関によって得られるか、または両データの類似度を別の手法で求めることによって得られる。プロセッサ１８は、複数の位置から得られた変位情報を使用し、前記空間的位置のうち複数の異なる空間的位置の組織変位の少なくとも一対の相関係数が求められる。伝播時間を求めるためには、複数の異なる対から求められた相関係数が使用される。複数の異なる空間的位置の伝播時間、たとえば、同じ入力によって方程式から得られた出力である、複数の異なる空間的位置の伝播時間が、同時に求められる。複数の異なる空間的位置の伝播時間は、前記空間的位置のうち複数の異なる空間的位置から得られたデータを使用して関数から求められた出力セットの形で得られる。前記出力セットの各値は、複数の空間的位置の組織変位に依存する。この複数の空間的位置には、伝播時間の計算が行われる対象である距離に含まれない空間的位置も含まれる。前記空間的位置のうち複数の異なる空間的位置の変位が、所与の伝播時間を求めるのに使用することができる。たとえば、複数の異なる位置で送信が行われることによって１つの同じ位置で引き起こされた組織変位の一対の相関係数と、１つの同じ位置で送信が行われることによって複数の異なる位置において引き起こされた組織変位の一対の相関係数と、別の複数の異なる位置において送信が行われることによって複数の異なる位置で引き起こされた組織変位の一対の相関係数の、合計３対の相関係数を計算する。

プロセッサ１８は、前記伝播時間から剪断波速度を計算する。別の剪断波情報を求めることもできる。特定のスキャンラインの遅延が求められ、該遅延と、該当の距離または範囲とから、剪断波速度が計算される。

プロセッサ１８は、上記の組織特性から、ディスプレイ２０にマッピングされた画像または表示値を生成して出力する。例えば各位置ごとに、剪断波速度、剪断弾性率または別の値が求められる。上記の値の大きさにより、色、色相、輝度および／または別の表示特性が変調される。プロセッサ１８は画素値（たとえばＲＧＢ）を決定するか、または画素値に変換されるスカラ値を決定する。画像はスカラ値または画素値として生成される。画像をビデオプロセッサ、ルックアップテーブル、カラーマップへ出力するか、またはディスプレイ２０へ直接出力することができる。

プロセッサ１８は、記憶装置２２または別の記憶装置に記憶されている命令にしたがって動作する。プロセッサ１８は、医療用超音波イメージングにおいて剪断波情報を求めるためにプログラミングされている。記憶装置２２は、不揮発性のコンピュータ読み出し可能記憶媒体である。本明細書において説明した処理、方法および/または技術を実現するための命令は、コンピュータ読み出し可能記憶媒体または記憶装置、例えば、キャッシュ、バッファ、RAM、リムーバル媒体、ハードディスクドライブまたは別のコンピュータ読み出し可能記憶媒体に記憶される。コンピュータで読み出し可能な記憶媒体には、種々のタイプの揮発性および非揮発性記憶媒体が含まれる。図面または本明細書において説明した機能、動作またはタスクは、コンピュータ読み出し可能記憶媒体に記憶されている１つまたは複数の命令セットに応答して実行される。この機能、ステップまたはタスクは、命令セット、記憶媒体、プロセッサまたは処理ストラテジなどの特定のタイプに依存することがなく、ソフトウェア、ハードウェア、集積回路、ファームウェア、マイクロコードその他などを単独または組み合わせて用いて実行される。同様に処理ストラテジーは、マルチプロセシング、マルチタスク処理、並列処理などを含んでいてもよい。１つの実施形態において、命令はローカルシステムまたはリモートシステムによって読み出せるように、リムーバブル記憶媒体装置に記憶されている。別の実施形態では、命令はコンピュータネットワークまたは電話回線を介して転送するために、遠隔地に記憶されている。さらに別の実施形態では、命令は所定のコンピュータ、ＣＰＵ、ＧＰＵまたはシステム内に記憶されている。

ディスプレイ２０は、２次元画像または３次元表現を表示するＣＲＴ，ＬＣＤ，プロジェクタ、プラズマディスプレイまたは別のディスプレイである。ディスプレイ２０は、剪断波速度等の剪断波情報を表現する１つまたは複数の画像を表示する。

これまで本発明を複数の実施形態に基づいて説明してきたが、本発明の趣旨から逸脱することなく種々の変更を加えることが可能であるものと理解されたい。したがって上記における詳細な説明は制限を意図しているもののではなく説明を意図しているものであり、本発明の精神および範囲を定めるのは、全ての同等のものを含む添付の特許請求の範囲であると解される。

Ｓｅｑ１，Ｓｅｑ２，Ｓｅｑ３ 加振信号の送信シーケンス
１０ 本発明のシステム
１２ 送信ビームフォーマ
１４ トランスデューサ
１６ 受信ビームフォーマ
１８ 画像処理プロセッサ
２０ ディスプレイ

Claims (9)

1. 医療用超音波イメージングシステムに備えられたプロセッサ（１８）によって剪断波情報を求めて画像を生成するための、該プロセッサ（１８）によって実行可能な命令を含むコンピュータプログラムであって、
   前記コンピュータプログラムは前記プロセッサ（１８）に対し、
   ・第１のスキャンライン上で送信された第１の音波に応答して第１の剪断波によって引き起こされた、時間に依存する組織の第１の変位を検出し（３２）、ただし該第１の変位は、該第１のスキャンラインから少なくとも第２のスキャンラインおよび第３のスキャンラインの分だけ離隔している第４のスキャンライン上に位置し、
   ・前記第３のスキャンライン上で送信された第２の音波に応答して第２の剪断波によって引き起こされた、時間に依存する組織の第２の変位を検出し（３２）、ただし、該第２の変位は前記第４のスキャンライン上に位置し、
   ・前記第１のスキャンライン上で送信された第３の音波に応答して第３の剪断波によって引き起こされた、時間に依存する組織の第３の変位を検出し（３２）、ただし、前記第３の変位は前記第２のスキャンライン上に位置し、
   ・複数の各ラテラル位置ごとに、
   前記第１の変位と前記第２の変位との最大相関に相当する遅延と、
   該第１の変位と前記第３の変位との最大相関に相当する遅延と、
   該第２の変位と該第３の変位との最大相関に相当する遅延と
   を計算し（３４）、
   ・前記各遅延に依存して、前記複数の各ラテラル位置ごとの組織中の剪断波速度を計算する（３５）
   ことを実行させることを特徴とする、コンピュータプログラム。
2. 時間に依存する前記第１の変位、前記第２の変位、および前記第３の変位を求める際に（３２）、所定の時間にわたって組織の動きの大きさをトラッキングする、請求項１記載のコンピュータプログラム。
3. 前記第１のスキャンライン、前記第２のスキャンライン、前記第３のスキャンライン、および前記第４のスキャンラインは、相互に隣接するスキャンラインである、請求項１記載のコンピュータプログラム。
4. 前記最大相関に相当する遅延を計算する際に（３４）、前記複数のラテラル位置と複数の可能性のある遅延とに依存して該遅延を計算する、請求項１記載のコンピュータプログラム。
5. 前記最大相関に相当する遅延を計算する際に（３４）、前記第１の変位と前記第２の変位との相関と、該第１の変位と前記第３の変位との相関と、該第２の変位と該第３の変位との相関との和を計算する（３４）、請求項１記載のコンピュータプログラム。
6. 前記剪断波速度を計算する際に（３５）、前記遅延および前記ラテラル位置に依存して該剪断波速度を計算する（３５）、請求項１記載のコンピュータプログラム。
7. 前記遅延を計算する際に（３４）、前記スキャンラインのうち複数のスキャンラインを同時に求める、請求項１記載のコンピュータプログラム。
8. 別のスキャンラインおよび別の奥行において、前記変位の検出（３２）と前記遅延の計算（３４）とを繰り返す、請求項１記載のコンピュータプログラム。
9. 直ぐ隣のスキャンライン上で連続して送信を行うのを回避するように、前記第１の音波と前記第２の音波と前記第３の音波とを送信する（３０）、請求項１記載のコンピュータプログラム。