(57)【要約】  第１の超音波パルスを生体組織に印加して、前記生体組織内にせん断波を生成し、超音波パルスを前記生体組織内に送信し、１以上の超音波信号を前記生体組織から受信し、前記生体組織内のせん断波を、前記受信された１以上の超音波信号に基づいて検出する。前記検出されたせん断波に関連する少なくとも１つの伝搬特性が決定され、前記決定された少なくとも１つの伝搬特性が表示される。
2015512273.tif 000002
【特許請求の範囲】
【請求項１】
第１の超音波パルスを生体組織に印加して、前記生体組織内にせん断波を生成するステップと、
第２の超音波パルスを前記生体組織内に送信するステップと、
前記第２超音波パルスに応答して発生された１以上の超音波信号を前記生体組織から受信するステップと、
前記生体組織内のせん断波を、前記受信された１以上の超音波信号に基づいて検出するステップと、
前記生体組織における複数の位置の各々におけるせん断波のピーク変位までの時間を決定するステップと、
前記複数の位置の各々におけるせん断波のピーク変位までの前記決定された時間に基づいて、せん断波伝搬速度の２乗を決定するステップと、
を含む方法。
【請求項２】
前記せん断波伝搬速度の２乗を決定するステップが、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出するステップと、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗を算出するステップと、
前記複数の位置の各々に関して、前記距離の２乗を前記時間の２乗で割算するステップと、
を含む請求項１に記載の方法。
【請求項３】
前記せん断波のピーク変位までの時間を決定するステップが、
前記受信された超音波信号間の相互相関、差の絶対値の和（ＳＡＤ）、差の２乗の和（ＳＳＤ）、差の絶対値の３乗の和（ＳＣＤ）、又は、差の絶対値の累乗の和（ＳＰＤ）を算出するステップ
を含む、請求項１に記載の方法。
【請求項４】
前記せん断波伝搬速度の２乗を決定するステップが、
前記複数の位置の各々における前記せん断波のピーク変位までの時間の２乗を算出するステップと、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗を算出するステップと、
前記複数の位置の各々の、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線を算出するステップと、
を含む請求項１に記載の方法。
【請求項５】
前記せん断波伝搬速度の２乗を決定するステップが、
前記複数の位置の各々における前記せん断波のピーク変位までの時間の２乗を算出するステップと、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗を算出するステップと、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線及び相関係数を算出するステップと、
を含む請求項１に記載の方法。
【請求項６】
前記せん断波伝搬速度の２乗を決定するステップが、
前記複数の位置のうちの、前記ピーク変位までの時間と距離の２乗との相関係数が予め設定された閾値よりも大きい複数の位置を決定するステップと、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との第２の回帰直線を算出するステップと、
前記第２の回帰直線の勾配を決定するステップと、
を含む請求項５に記載の方法。
【請求項７】
前記せん断波伝搬速度の２乗を決定するステップが、
前記複数の位置のうちの２つの位置におけるせん断波のピーク変位までの時間の差の２乗を算出するステップと、
前記複数の位置のうちの前記２つの位置間の距離の２乗を算出するステップと、
前記距離の２乗を前記時間の差の２乗で割算するステップと、
を含む請求項１に記載の方法。
【請求項８】
前記せん断波伝搬速度の２乗を決定するステップが、
前記せん断波の起源におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算するステップと、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出するステップと、
前記複数の位置の各々の、前記せん断波の前記起源からの距離の２乗を算出するステップと、
前記複数の位置の各々に関し、前記距離の２乗を前記時間の２乗で割算するステップと、
を含む、請求項１に記載の方法。
【請求項９】
前記せん断波伝搬速度の２乗を決定するステップが、
前記せん断波の起源におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算するステップと、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出するステップと、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗を算出するステップと、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線を算出するステップと、
を含む請求項１に記載の方法。
【請求項１０】
前記せん断波伝搬速度の２乗を決定するステップが、
前記せん断波の起源におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算するステップと、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出するステップと、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗を算出するステップと、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線及び相関係数を算出するステップと、
を含む請求項１に記載の方法。
【請求項１１】
前記せん断波伝搬速度の２乗を決定するステップが、
前記複数の位置の１つを選択するステップと、
前記選択された位置におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算するステップと、
複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出するステップと、
前記複数の位置の各々の、前記選択された位置からの距離の２乗を算出するステップと、
前記複数の位置の各々に関し、前記距離の２乗を前記時間の２乗で割算するステップと、
を含む請求項１に記載の方法。
【請求項１２】
前記せん断波伝搬速度の２乗を決定するステップが、
前記複数の位置の１つを選択するステップと、
前記選択された位置におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算するステップと、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出するステップと、
前記複数の位置の各々の、前記選択された位置からの距離の２乗を算出するステップと、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記選択された位置からの距離の２乗との回帰直線を算出するステップと、
を含む請求項１に記載の方法。
【請求項１３】
前記せん断波伝搬速度の２乗を決定するステップが、
前記複数の位置の１つを選択するステップと、
前記選択された位置におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算するステップと、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出するステップと、
前記複数の位置の各々の、前記選択された位置からの距離の２乗を算出するステップと、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記選択された位置からの距離の２乗との回帰直線及び相関係数を算出するステップと、
を含む請求項１に記載の方法。
【請求項１４】
前記せん断波伝搬速度の２乗を決定するステップが、
前記複数の位置のうちの、前記ピーク変位までの時間と前記距離の２乗との相関係数が予め設定された閾値よりも大きい複数の位置を決定するステップと、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記選択された位置からの距離の２乗との第２の回帰直線を算出するステップと、
前記第２の回帰直線の勾配を決定するステップと、
を含む請求項１３に記載の方法。
【請求項１５】
前記せん断波伝搬速度の２乗を決定するステップが、
前記複数の位置のうちの、前記ピーク変位までの時間と前記距離の２乗との相関係数が予め設定された閾値よりも大きい複数の位置を決定するステップと、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記せん断波の起源からの距離の２乗との第２の回帰直線を算出するステップと、
前記第２の回帰直線の勾配を決定するステップと、
を含む請求項１０に記載の方法。
【請求項１６】
コンピュータ実行可能なプログラムコードを格納している一時的でない媒体であって、前記プログラムコードが、演算装置により、
第１の超音波パルスを生体組織に印加して、前記生体組織内にせん断波を生成し、
第２の超音波パルスを前記生体組織内に送信し、
前記第２超音波パルスに応答して発生された１以上の超音波信号を前記生体組織から受信し、
前記生体組織内のせん断波を、前記受信された１以上の超音波信号に基づいて検出し、
前記生体組織における複数の位置の各々におけるせん断波のピーク変位までの時間を決定し、
前記複数の位置の各々におけるせん断波のピーク変位までの前記測定された時間に基づいてせん断波伝搬速度の２乗を決定する、
ことが実行可能である、媒体。
【請求項１７】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出し、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗を算出し、且つ、
前記複数の位置の各々に関して、前記距離の２乗を前記時間の２乗で割算する、
ように演算装置により実行可能なプログラムコードを含む請求項１６に記載の媒体。
【請求項１８】
前記せん断波のピーク変位までの時間を決定するために演算装置により実行可能な前記プログラムコードが、
前記受信された超音波信号間の相互相関、差の絶対値の和（ＳＡＤ）、差の２乗の和（ＳＳＤ）、差の絶対値の３乗の和（ＳＣＤ）、又は、差の絶対値の累乗の和（ＳＰＤ）を算出する、
ように演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項１９】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記複数の位置の各々における前記せん断波のピーク変位までの時間の２乗を算出し、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗を算出し、且つ、
前記複数の位置の各々の、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線を算出する、
ように演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項２０】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記複数の位置の各々における前記せん断波のピーク変位までの時間の２乗を算出し、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗を算出し、且つ、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線及び相関係数を算出する、
ように演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項２１】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記複数の位置のうちの、前記ピーク変位までの時間と距離の２乗との相関係数が予め設定された閾値よりも大きい複数の位置を決定し、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との第２の回帰直線を算出し、且つ、
前記第２の回帰直線の勾配を決定する、
ように演算装置により実行可能なプログラムコードを含む、請求項２０に記載の媒体。
【請求項２２】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記複数の位置のうちの２つの位置におけるせん断波のピーク変位までの時間の差の２乗を算出し、
前記複数の位置のうちの前記２つの位置間の距離の２乗を算出し、且つ、
前記距離の２乗を前記時間の差の２乗で割算する、
ように演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項２３】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記せん断波の起源におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算し、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出し、
前記複数の位置の各々の、前記せん断波の前記起源からの距離の２乗を算出し、且つ、
前記複数の位置の各々に関し、前記距離の２乗を前記時間の２乗で割算する、
ように演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項２４】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記せん断波の起源におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算し、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出し、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗を算出し、且つ、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線を算出する、
ように演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項２５】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記せん断波の起源におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算し、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出し、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗を算出し、且つ、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線及び相関係数を算出する、
ように演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項２６】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記複数の位置の１つを選択し、
前記選択された位置におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算し、
複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出し、
前記複数の位置の各々の、前記選択された位置からの距離の２乗を算出し、且つ、
前記複数の位置の各々に関し、前記距離の２乗を前記時間の２乗で割算する、
ように演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項２７】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記複数の位置の１つを選択し、
前記選択された位置におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算し、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出し、
前記複数の位置の各々の、前記選択された位置からの距離の２乗を算出し、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記選択された位置からの距離の２乗との回帰直線を算出する、
ように演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項２８】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記複数の位置の１つを選択し、
前記選択された位置におけるピーク変位までの時間を、前記複数の位置の各々におけるせん断波のピーク変位までの時間から減算し、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗を算出し
前記複数の位置の各々の、前記選択された位置からの距離の２乗を算出し、且つ、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記選択された位置からの距離の２乗との回帰直線及び相関係数を算出する、
ように演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項２９】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記複数の位置のうちの、前記ピーク変位までの時間と前記距離の２乗との相関係数が予め設定された閾値よりも大きい複数の位置を決定し、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記選択された位置からの距離の２乗との第２の回帰直線を算出し、且つ、
前記第２の回帰直線の勾配を決定する、
ように演算装置により実行可能なプログラムコードを含む、請求項２８に記載の媒体。
【請求項３０】
せん断波伝搬速度の２乗を決定するために演算装置により実行可能な前記プログラムコードが、
前記複数の位置のうちの、前記ピーク変位までの時間と前記距離の２乗との相関係数が予め設定された閾値よりも大きい複数の位置を決定し、
前記複数の位置の各々に関し、前記ピーク変位までの時間の２乗と、前記せん断波の起源からの距離の２乗との第２の回帰直線を算出し、且つ、
前記第２の回帰直線の勾配を決定する、
ように演算装置により実行可能なプログラムコードを含む、請求項２５に記載の媒体。
【請求項３１】
システムであって、
メモリと、
前記メモリに格納されたコンピュータ実行可能なプログラムコードを実行するためのプロセッサとを備え、前記プログラムコードが、前記システムに、
第１の超音波パルスを生体組織に印加して、前記生体組織内にせん断波を生成させ、
第２の超音波パルスを前記生体組織内に送信させ、
前記第２の超音波パルスに応答して発生された１以上の超音波信号を前記生体組織から受信させ、
前記生体組織内の前記せん断波を、前記受信した１以上の超音波信号に基づいて検出させ、
前記生体組織における複数の位置の各々における前記せん断波のピーク変位までの時間を決定させ、
前記複数の位置の各々における前記せん断波の前記決定されたピーク変位までの時間に基づいて、せん断波伝搬速度の２乗を決定させる、
システム。
【請求項３２】
前記せん断波伝搬速度の２乗の決定が、
前記複数の位置の各々における前記せん断波のピーク変位までの時間の２乗の算出と、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗の算出と、
前記複数の位置の各々に関する、前記距離の２乗の、前記時間の２乗による割算と、
を含む、請求項３１に記載のシステム。
【請求項３３】
前記せん断波のピーク変位までの時間の決定が、
前記受信された超音波信号間の相互相関、差の絶対値の和（ＳＡＤ）、差の２乗の和（ＳＳＤ）、差の絶対値の３乗の和（ＳＣＤ）、又は、差の絶対値の累乗の和（ＳＰＤ）の算出、
を含む、請求項３１に記載のシステム。
【請求項３４】
前記せん断波伝搬速度の２乗の決定が、
前記複数の位置の各々における前記せん断波のピーク変位までの時間の２乗の算出と、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗の算出と、
前記複数の位置の各々の、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線の算出と、
を含む、請求項３１に記載のシステム。
【請求項３５】
前記せん断波伝搬速度の２乗の決定が、
前記複数の位置の各々における前記せん断波のピーク変位までの時間の２乗の算出と、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗の算出と、
前記複数の位置の各々に関する、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線及び相関係数の算出と、
を含む、請求項３１に記載のシステム。
【請求項３６】
前記せん断波伝搬速度の２乗の決定が、
前記複数の位置のうちの、前記ピーク変位までの時間と距離の２乗との相関係数が予め設定された閾値よりも大きい複数の位置の決定と、
前記複数の位置の各々に関する、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との第２の回帰直線の算出と、
前記第２の回帰直線の勾配の決定と、
を含む、請求項３５に記載のシステム。
【請求項３７】
前記せん断波伝搬速度の２乗の決定が、
前記複数の位置のうちの２つの位置におけるせん断波のピーク変位までの時間の差の２乗の算出と、
前記複数の位置のうちの前記２つの位置間の距離の２乗の算出と、
前記距離の２乗の、前記時間の差の２乗による割算と、
を含む、請求項３１に記載のシステム。
【請求項３８】
前記せん断波伝搬速度の２乗の決定が、
前記せん断波の起源におけるピーク変位までの時間の、前記複数の位置の各々におけるせん断波のピーク変位までの時間からの減算と、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗の算出と、
前記複数の位置の各々の、前記せん断波の前記起源からの距離の２乗の算出と、
前記複数の位置の各々に関する、前記距離の２乗の、前記時間の２乗による割算と、
を含む、請求項３１に記載のシステム。
【請求項３９】
前記せん断波伝搬速度の２乗の決定が、
前記せん断波の起源におけるピーク変位までの時間の、前記複数の位置の各々におけるせん断波のピーク変位までの時間からの減算と、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗の算出と、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗の算出と、
前記複数の位置の各々に関する、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線の算出と、
を含む、請求項３１に記載のシステム。
【請求項４０】
前記せん断波伝搬速度の２乗の決定が、
前記せん断波の起源におけるピーク変位までの時間の、前記複数の位置の各々におけるせん断波のピーク変位までの時間からの減算と、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗の算出と、
前記複数の位置の各々の、前記せん断波の起源からの距離の２乗の算出と、
前記複数の位置の各々に関する、前記ピーク変位までの時間の２乗と、前記せん断波の前記起源からの距離の２乗との回帰直線及び相関係数の算出と、
を含む、請求項３１に記載のシステム。
【請求項４１】
前記せん断波伝搬速度の２乗の決定が、
前記複数の位置の１つの選択と、
前記選択された位置におけるピーク変位までの時間の、前記複数の位置の各々におけるせん断波のピーク変位までの時間からの減算と、
複数の位置の各々におけるせん断波のピーク変位までの時間の２乗の算出と、
前記複数の位置の各々の、前記選択された位置からの距離の２乗の算出と、
前記複数の位置の各々に関する、前記距離の２乗の、前記時間の２乗による割算と、
を含む、請求項３１に記載のシステム。
【請求項４２】
前記せん断波伝搬速度の２乗の決定が、
前記複数の位置の１つの選択と、
前記選択された位置におけるピーク変位までの時間の、前記複数の位置の各々におけるせん断波のピーク変位までの時間からの減算と、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗の算出と、
前記複数の位置の各々の、前記選択された位置からの距離の２乗の算出と、
前記複数の位置の各々に関する、前記ピーク変位までの時間の２乗と、前記選択された位置からの距離の２乗との回帰直線の算出と、
を含む、請求項３１に記載のシステム。
【請求項４３】
前記せん断波伝搬速度の２乗の決定が、
前記複数の位置の１つの選択と、
前記選択された位置におけるピーク変位までの時間の、前記複数の位置の各々におけるせん断波のピーク変位までの時間からの減算と、
前記複数の位置の各々におけるせん断波のピーク変位までの時間の２乗の算出と、
前記複数の位置の各々の、前記選択された位置からの距離の２乗の算出と、
前記複数の位置の各々に関する、前記ピーク変位までの時間の２乗と、前記選択された位置からの距離の２乗との回帰直線及び相関係数の算出と、
を含む、請求項３１に記載のシステム。
【請求項４４】
前記せん断波伝搬速度の２乗の決定が、
前記複数の位置のうちの、前記ピーク変位までの時間と前記距離の２乗との相関係数が予め設定された閾値よりも大きい複数の位置の決定と、
前記複数の位置の各々に関する、前記ピーク変位までの時間の２乗と、前記選択された位置からの距離の２乗との第２の回帰直線の算出と、
前記第２の回帰直線の勾配の決定と、
を含む、請求項４３に記載のシステム。
【請求項４５】
前記せん断波伝搬速度の２乗の決定が、
前記複数の位置のうちの、前記ピーク変位までの時間と前記距離の２乗との相関係数が予め設定された閾値よりも大きい複数の位置の決定と、
前記複数の位置の各々に関する、前記ピーク変位までの時間の２乗と、前記せん断波の起源からの距離の２乗との第２の回帰直線の算出と、
前記第２の回帰直線の勾配の決定と、
を含む、請求項４０に記載のシステム。
【請求項４６】
前記せん断波のピーク変位までの時間を決定するステップが、
カラードプラシフト周波数、カラードプラ位相、又は、カラードプラ速度を、前記受信された超音波信号のＩ−Ｑ信号に基づいて算出するステップを含む、請求項１に記載の方法。
【請求項４７】
前記せん断波のピーク変位までの時間を決定するために演算装置により実行可能な前記プログラムコードが、
カラードプラシフト周波数、カラードプラ位相、又は、カラードプラ速度を、前記受信された超音波信号のＩ−Ｑ信号に基づいて算出するために演算装置により実行可能なプログラムコードを含む、請求項１６に記載の媒体。
【請求項４８】
前記せん断波のピーク変位までの時間の決定が、
カラードプラシフト周波数、カラードプラ位相、又は、カラードプラ速度の、前記受信された超音波信号のＩ−Ｑ信号に基づいた算出を含む、請求項３１に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本明細書に記載するシステム及び方法は、概して、超音波イメージング（撮像）の分野に関する。より詳細には、以下に記載する実施形態は、組織内でのせん断波速度の２乗を測定するための方法及びシステムに関する。
【背景技術】
【０００２】
病理的状態においては、生理学的状態下でよりも硬い軟部組織が生じ得る。従って、医師は触診により身体内の硬組織の位置を特定し、それにより病理的状態を識別する。例えば、乳癌は、健全な乳房組織よりも一般的に硬いことが知られており、触診法により、硬いしこりとして検知され得る。
【０００３】
組織内でのせん断波伝搬速度の２乗は、以下の式に示すように、組織の硬さ（ヤング率）に関連している。
【数１】
2015512273.tif 000003
式中、ｃはせん断波の伝搬速度であり、Ｅはヤング率であり、ρは組織密度である。従って、組織を通過するせん断波の伝搬速度の２乗を測定することにより、組織内の癌又はその他の病理的状態を検出し得る。
【０００４】
せん断波は、強い超音波パルスを組織に印加することにより組織内に生成され得る。超音波パルスは、高い振幅及び長い持続時間（例えば、約１００マイクロ秒）を呈し得る。超音波パルスは音響放射力を発生し、この力が組織を押圧し、それにより、組織の層を超音波パルスの方向に沿ってスライドさせる。このような組織のスライド（せん断）移動を、せん断波とみなすことができ、せん断波は、低周波数（例えば、１０Ｈｚ〜５００Ｈｚ）であり、超音波パルスの方向に対して垂直な方向に伝搬し得る。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
組織の運動が一般に軸方向（すなわち、超音波パルスの方向）であるため、せん断波は慣用の超音波ドプラ技術を用いて検出され得る。この点で、超音波ドプラ技術は、軸方向の速度を検出するのに最適である。或いは、せん断波を、音響放射力により生じる組織変位を測定することにより検出することも可能である。
【課題を解決するための手段】
【０００６】
せん断波伝搬速度の２乗を正確に測定するために、せん断波を複数の位置にて検出する必要がある。せん断波伝搬速度の２乗の正確な測定値を得るために、新しいアルゴリズムが必要であろう。
【図面の簡単な説明】
【０００７】
【図１】音響放射力による、せん断波生成の図である。
【図２】幾つかの実施形態の超音波イメージングシステムの図である。
【図３】慣用の超音波イメージングシステムの図である。
【図４】複数の超音波送信／受信ビームを示す図である。
【図５】超音波送信ビーム及び複数の超音波受信ビームを示す図である。
【図６】せん断波伝搬速度の２乗のカラーコーディング（色分け）の図である。
【図７】せん断波伝搬速度の２乗のカラーコーディングの図である。
【図８】音響放射力によるせん断波の生成及びせん断波の伝搬を示す図である。
【図９】せん断波のスライド移動を示す図である。
【図１０】せん断波の伝搬を示す図である。
【図１１】せん断波の伝搬を示す図である。
【図１２】組織内での、せん断波伝搬速度の２乗のカラーコーディング画像の例である。
【図１３】音響放射力により生じた組織変位を示す図である。
【図１４】ＲＧＢ表示で構成されたカラーコーディングバーによる、せん断波速度の２乗ｃ２のスケール図である。
【図１５】超音波振動子に対する超音波座標系を示す図である。
【図１６】複数の位置における時間の経過による変位変化を示す図である。
【図１７】均一の領域における、せん断波の起源からの距離の２乗に対するピーク変位までの時間の２乗をプロットした図である。
【図１８】不均一の領域における、せん断波の起源からの距離の２乗に対するピーク変位までの時間の２乗をプロットした図である。
【図１９】せん断波の起源及び測定ポイントを示す図である。
【図２０】幾つかの実施形態によるプロセスのフローチャートである。
【図２１】幾つかの実施形態によるプロセスのフローチャートである。
【図２２】幾つかの実施形態によるプロセスのフローチャートである。
【発明を実施するための形態】
【０００８】
本発明の実施形態を、添付図面を参照しつつ以下に記載する。図面全体を通じて、類似の番号は類似の要素を示す。実施形態を詳細に説明する前に、これらの実施形態の適用が、以下の記載にて説明され又は図面に示されている例の詳細に限定されないことを理解されたい。その他の実施形態も、様々な用途及び様々な方法で実施又は実行し得る。また、本明細書中で使用される用語及び術語が、説明のためのものであり、限定的なものとみなされるべきではないことが理解されよう。本明細書中の「含む」(”including”)、「備える」(”comprising”)、又は「有する」(”having”)、及びこれらの変化形を用いることは、その後に記載される項目及びその等価物、並びに追加項目を包含することを意味する。用語「取り付けられる」（”mounted”）、「接続される」（”connected”）、及び「結合される」（”coupled”）は広義に使用され、直接的な及び間接的な取り付け、接続、及び結合を含む。さらに、「接続される」及び「結合される」は、物理的な又は機械的な接続又は結合に限定されない。
【０００９】
図１に示されているように、強い超音波パルス１２０により音響放射力が生成される。超音波パルス１２０は、高い振幅、及び、長い持続時間（例えば、約１００マイクロ秒）を有する。超音波パルス１２０は超音波振動子アレイ１１０から送信される。超音波パルス１２０は生体組織１６０内の焦点１３０に集束し、焦点１３０にて、組織１６０を押圧する音響放射力となる。超音波パルス１２０は複数回送信されてもよく、また、送信された複数の超音波パルスの各々が、異なる焦点に集束し得る。
【００１０】
組織１６０は、主に超音波パルス１２０の軸方向に押されて、せん断波１４０，１５０を生成し、これらのせん断波は、横方向に伝搬するか、又は、軸方向（すなわち鉛直方向）以外の方向に伝搬し得る。せん断波１４０，１５０の伝搬速度の２乗は、組織１６０の硬さ（すなわち、ヤング率）に依存する。式１に示したように、組織の硬さが大きくなると、せん断波の伝搬速度の２乗が大きくなる。癌などの病理的状態が組織の硬さを増大させることがあり、従って、このような状態は、伝搬速度の２乗を決定することにより診断され得る。例えば、せん断波の伝搬速度の２乗は、組織の状態に応じて１ｍ２／ｓ２から１００ｍ２／ｓ２まで変化し得る。
【００１１】
せん断波が組織の移動（又は運動）により特徴付けられるため、せん断波は、超音波ドプラ技術により検出され得る（例えば、米国特許第４５７３４７７号、米国特許第４６２２９７７号、米国特許第４６４１６６８号、米国特許第４６５１７４２号、米国特許第４６５１７４５号、米国特許第４７５９３７５号、米国特許第４７６６９０５号、米国特許第４７６８５１５号、米国特許第４７７１７８９号、米国特許第４７８０８３７号、米国特許第４７９９４９０号、及び、米国特許第４９６１４２７号を参照されたい）。この組織の移動（又は運動）を検出するために、超音波パルスを組織に複数回送信し、超音波は組織内の散乱体により散乱され、そして、超音波振動子により受信超音波信号として受信される。超音波アレイ振動子から受信された超音波信号に、集束及びステアリングのために遅延及び／又は位相回転を加えた後、これらの信号を、フィルタリング、増幅、デジタル化、円滑化（アポダイズ）(apotize)、及び、ビーム形成（すなわち、合算）する。これらの処理ステップの順序を入れ替えてもよい。ビーム形成された受信ＲＦ超音波信号に直交復調（直交検波）を行い、これにより、複素ドプラＩ−Ｑ信号を得る。カラードプラ技術において、超音波は、パルス繰返し周波数（ＰＲＦ）で送信され、その速度が、受信超音波信号における周波数のシフト（ドプラシフト周波数）として検出される。受信された超音波は、送信された超音波周波数と同一周波数の同相（０度）及び直角（９０度）の参照信号と混合されて、複素Ｉ−Ｑドプラ信号になる。
【００１２】
一般に、ドプラシフト周波数を導出するために複素Ｉ−Ｑ信号を用いる。これは、ドプラシフト周波数と血流速度とが以下の関係を有するためである。
【数２】
2015512273.tif 000004
式中、Δｆはドプラシフト周波数であり、ｆｔは送信周波数であり、ｖは血流速度であり、θは、超音波ビーム方向と速度ベクトルとが成す角度であり、ｃＳは音速である。従って、ドプラシフト周波数は、速度方向と超音波ビーム方向との間の角度に依存し、超音波カラードプラシステムにより得られる測定値である。
【００１３】
カラードプラの場合、サンプリング信号の個数が数個に限られるであろう。従って、通常、自己相関技術を用いてＩ−Ｑ信号間の位相差を特定し、そしてそれによりドプラシフト周波数及び速度を、以下のように決定する。カラードプラのＩ−Ｑ信号ｚ（ｍ）＝ｘ（ｍ）＋ｊｙ（ｍ）を用いて「自己相関」(“auto-correlation”)Ｒを、以下の式に示すように算出する。式中、ｚ（ｍ）は複素Ｉ−Ｑドプラ信号であり、ｘ（ｍ）は同相（実数）信号であり、ｙ（ｍ）は直角位相（虚数）信号であり、ｍは信号番号を示し、ｊは虚数単位であり、＊は複素共役を示す。
【数３】
2015512273.tif 000005
【００１４】
Ｒの実数（Ｒｅａｌ（Ｒ））部及び虚数（Ｉｍａｇ（Ｒ））部を用いて、以下の式に示すように、位相φを求める。
【数４】
2015512273.tif 000006
【００１５】
ｔａｎ−１が、通常、−０．５π〜０．５πに限られるため、−π〜πの範囲のφを導出するために、複素座標における複素数値Ｒの位置を用いてもよい。次いで、位相（すなわち、カラードプラ位相）φを、以下の式に示すように、ドプラシフト周波数（すなわち、カラードプラシフト周波数）に関連付ける。
【数５】
2015512273.tif 000007
【００１６】
このようにして、受信した複素ベースバンド超音波信号間の自己相関Ｒを求め、それにより、組織の速度又は移動を検出する。
【００１７】
組織の移動は、組織領域のフィールド内の横方向の複数のポイントにて複数の超音波ビーム（例えば、図５の５４０，５４５，５５０）により検出され、これにより移動が監視される。この移動は、これらの複数の横方向ポイント（又は複数の超音波ビーム）におけるせん断波の動きを反映している。従って、せん断波の横方向の伝搬速度を、検出された組織移動から決定し得る。
【００１８】
別の方法として、図１３に示されているように、せん断波を、強力な超音波パルスにより生じた音響放射力により引き起こされた組織変位を測定することにより検出することも可能である。組織１３１０は、音響放射が印加される前は位置１３２０にあり、音響放射力が印加された後、位置１３３０に移動する。強力な超音波パルスにより生じた組織変位を測定するために、超音波パルスが超音波振動子１３０５から組織に送信され、次いで、超音波パルスは組織内の散乱体から散乱され、振動子１３０５に戻され、そして、受信超音波信号として振動子１３０５により受信される。超音波パルスを所定の深さにて集束させる。これは、得られた受信超音波信号の信号対ノイズ比を、非集束超音波パルスと比較して高めるためである。組織から受信した超音波信号の相関を用いて、音響放射力による組織１３１０の（位置１３２０から位置１３３０への）変位１３４０を求めることができ、そしてその後、組織１３１０を追跡（トラッキング）し得る。これにより、超音波パルスは、音響放射力によりせん断波が生成された後に、せん断波を追跡し得る。
【００１９】
第１の超音波パルスにより生じ、且つ、音響放射力が印加される前に組織１３１０から受信された超音波信号と、音響放射力が印加された後に第２の超音波パルスから得られた受信超音波信号とを相互相関させ、それにより、受信超音波信号間の最適マッチングを見つける。この最適マッチングは、組織と、音響放射力による組織変位とを追跡するために最大相関値を求めることにより見出され得る。従って、組織変位が観察又は測定されたとき、せん断波が検出される。変位と組織速度とは、変位が組織速度ｖＳの時間積分値∫ｖＳｄｔであることにより関連し得る。従って、組織変位は、カラードプラ速度の時間積分を算出することにより求められ得る。受信超音波信号は、ＲＦ（無線周波数）信号、ＩＦ（中間周波数）信号、又は、復調後のベースバンド信号であり得る。或いは、組織の歪みを求めるために変位をさらに微分してもよく、そしてこの歪みを用いることにより、せん断波の伝搬速度の２乗を検出し得る。
【００２０】
上記段落における信号の相互相関ＣＣ（ｔ，τ）を、以下のように数式で示し得る。
【数６】
2015512273.tif 000008
式中、ＣＣ（ｔ，τ）は相互相関であり、Ｓ１（ｔ’）は、第１の超音波送信から受信した信号であり、Ｓ２（ｔ’−τ）は、第２の超音波送信からの受信超音波信号であり、Ｗはウィンドウ（窓）長さ、ｔは時間、ｔ’は時間、τは時間変位である。最大相互相関（又は最適マッチング）をもたらす時間変位値τにより組織変位が決定される。空間分解能を高めるために、相互相関前に、補間関数（例えば、３次スプライン）を用いて信号の補間を実行してもよい。
【００２１】
相互相関は、以下のように、差の絶対値の和（ＳＡＤ）、差の２乗の和（ＳＳＤ）、差の絶対値の３乗の和（ＳＣＤ）、又は、差の絶対値の累乗の和（ＳＰＤ）で置き換えられ得る。
【数７】
2015512273.tif 000009
【数８】
2015512273.tif 000010
【数９】
2015512273.tif 000011
【数１０】
2015512273.tif 000012
Ｓ１は、第１の超音波送信からの、変位前の受信超音波信号であり、Ｓ２は、第２の超音波送信からの、変位後の受信超音波信号である。Ｎは、信号ウィンドウ内の信号数である。ｋは、信号数によるウィンドウ変位であり、τと等価である。ｌはウィンドウの位置であり、ｐは実数である。ＳＡＤ，ＳＳＤ，ＳＣＤ，ＳＰＤに関し、組織変位は、ＳＡＤ，ＳＳＤ，ＳＣＤ，ＳＰＤの各々を最小にする（最適マッチング）ｋの値に基づいて決定される。
【００２２】
図８及び図９は、せん断波の生成及び検出を詳細に示す。強力な超音波パルス８２０を超音波振動子８１０，９１０から組織８６０，９６０に１回又は複数回印加し、それにより、せん断波（超音波パルスにより生じる音響放射力により引き起こされる）の振幅を増大させる。せん断波は、組織内で非常に急速に減衰するため、振幅が大きいほど伝搬距離が大きくなる。１又は複数の超音波パルスを、１つの焦点、又は、異なる複数の焦点に集束させ得る。超音波パルスは、組織の層を押圧する音響放射力を生成し、これにより、図９に示されているように、主に軸（鉛直）方向における組織移動８３０，９１０が生じる。組織層の移動９１０は、隣接する組織層の主に軸方向の移動９２０，９２５を生じさせる。そして、組織層の移動９２０，９２５が、隣の組織層の移動９３０，９３５を生じ、そして、この移動が、隣接する組織層の移動９４０，９４５を生じさせる。このような連続した組織移動が、図８に示されているような、せん断波８４０，８５０の横（水平）方向の伝搬として表れる。音響放射力により引き起こされる組織の移動（又は運動）は、主に軸方向であるため、この運動を、軸方向運動に感受性のカラードプラ技術により検出し得る。
【００２３】
例えば、カラードプラ技術は、先に述べたように、また当分野にて知られているように、数個の超音波パルスを送受信し、受信された超音波信号間の位相差を特定し、そして、組織又は血流の速度を、自己相関技術を用いて算出する。速度だけでなく、カラードプラ信号の分散及びパワーも算出し得る。これらのパラメータの１つを用いて、移動している組織及び血液の従来の表示と同様に、せん断波を表示し得る。これを、図１０，図１１に示す。せん断波１０４０（１１４０），１０５０（１１５０）が、所定の時間を表すカラードプラフレーム内で決定されると、次の瞬間に又は次のフレームにおいて、せん断波１０６０（１１６０），１０７０（１１７０）が決定されると考えられよう。せん断波を追跡し、そしてせん断波伝搬の動画を作成するために、せん断波のより多くの画像フレームを取得してもよい。別の実施形態においては、音響放射力による組織変位も検出し得る。
【００２４】
図１０及び図１１は、２つの時点におけるせん断波の伝搬を示す。矢印１０８０，１０９０で示されているような局所せん断波伝搬速度の２乗を、２つの時点におけるせん断波の２つの画像を相関させることにより導出し得る。以下に記載するように局所せん断波伝搬速度の２乗を２次元画像で提示するために、より多数のせん断波の画像フレームを用いて、より多くの画像領域におけるせん断波の伝搬を追跡することも可能である。
【００２５】
第１のフレーム信号Ｓ１と第２のフレーム信号Ｓ２との相関係数（ＣＣＶ）を、以下に示すように、スペックルトラッキングとして求め得る。
【数１１】
2015512273.tif 000013
2015512273.tif 000014
座標系（ｘ，ｙ，ｚ）は、図１５において、超音波振動子１５１０に対して示されている。縦軸ｙは図１５の紙面に対して垂直であるが、説明のためにわずかに異なる状態で示されている。
【００２６】
最大相関係数をもたらす変位Ｘ，Ｚを決定し、次いで、変位の２乗を決定して、第１フレーム信号と第２フレーム信号との間の時間の２乗で割算することにより、せん断波伝搬速度の２乗が得られる。
【００２７】
１次元（１Ｄ）の場合と同様に、相関係数を、以下に示すように、差の絶対値の和（ＳＡＤ）、差の２乗の和（ＳＳＤ）、差の絶対値の３乗の和（ＳＣＤ）、及び、差の絶対値の累乗の和（ＳＰＤ）で置き換え得る。
【数１２】
2015512273.tif 000015
【数１３】
2015512273.tif 000016
【数１４】
2015512273.tif 000017
【数１５】
2015512273.tif 000018
ｐは実数であり、ｍ及びｎは整数である。２次元（２Ｄ）スペックルトラッキングを１Ｄスペックルトラッキングに近似化して、せん断波伝搬速度の２乗を得ることもできる。数式は、変位測定で用いたものに類似するであろう。
【００２８】
或いは、せん断波の伝搬速度の２乗を導出するために、以下に示すようなせん断波の式（１６）を用いてもよい。
【数１６】
2015512273.tif 000019
式中、ｉ＝ｘ，ｙ，ｚであり、ρは組織密度であり、μはせん断弾性率であり、ｕｉは変位ベクトルであり、図１５に示すように、ｘは横座標、ｙは縦座標、ｚは軸座標である。非圧縮性材料に関し、ヤング率Ｅ及びせん断弾性率μは以下の関係にある。
【数１７】
2015512273.tif 000020
従って、せん断波の伝搬速度の２乗は、以下の式に示すように、密度に対するせん断弾性率の比として求められる。
【数１８】
2015512273.tif 000021
式１６における変位成分のうちの１つｕｚを、先に述べたように相互相関により決定し得る。式１６のｚ成分と式１８とを組み合わせることにより、せん断波伝搬速度の２乗が、以下のように得られる。
【数１９】
2015512273.tif 000022
従って、せん断波伝搬速度の２乗が、変位の２次時間導関数と変位の２次空間導関数との比として得られる。縦方向における変位の空間導関数∂2uz／∂y2が、その他の空間導関数と比較して無視できる値と考えられるため、その他の測定値からせん断波伝搬速度の２乗を得ることができよう。
【００２９】
せん断波を、頻繁に、すなわち速い速度又はフレームレートで監視及び追跡することが望ましい。フレームレートを高速化するために、図５に示したように、幅広の集束超音波パルス５２０を送信でき、そして、複数の超音波信号５４０，５４５，５５０を同時に受信し得る。受信した超音波ビームを、先に述べたように用いてせん断波を検出し、そして、そのせん断波から、せん断波伝搬特性（すなわち、速度の２乗）を導き出す。集束された送信超音波ビーム５２０は、得られた受信超音波ビームの良好な信号対ノイズ比を、せん断波の検出中に維持するのに特に適し得る。
【００３０】
幾つかの実施形態においてフレームレートを上げるために、図４に示したように、複数の超音波ビーム（パルス）を組織フィールドに同時に印加及び送信すると、送信された超音波パルスごとに複数の超音波ビーム（パルス）が受信される。図４において、超音波パルス４２０，４３０が、超音波振動子アレイ４１０から生体組織４８０に向けて同時に送信されている。送信された超音波パルス４２０，４３０の各々に対して複数の超音波受信信号４４０，４４５，４６５，４６０，４６５，４７０が同時に受信される。複数の超音波パルスは、同時に、又はほぼ同時期に送信され得る。複数の超音波パルスは同時に送信されてよい。或いは、第２の超音波パルスを、第１の超音波パルスが送信された後、この第１の超音波パルスが超音波場の最深部から超音波振動子に戻るまでの間に前に送信してもよい。この送信方法によりフレームレートが向上する。
【００３１】
図４は、同時に送信された２つの超音波パルスの例を示しているが、２つより多くの送信超音波パルスも用いられ得る。幾つかの実施形態において、複数の同時の超音波信号の、より良好な分離のために、コード化された超音波波形を送信し得る。例えば、チャープコード、バーカーコード、ゴーレイコード、又はアダマールコードを用いて、超音波パルスをより良好に分離し得る。この場合もまた、受信された信号を、先に述べた方法を用いて分析し、それにより、複数のポイントにおける組織移動を特定し、そしてこの結果から、せん断波伝搬特性を導出する。
【００３２】
せん断波の画像を、イメージングフィールド（撮像領域）内の複数のポイントにて検出された運動（又は速度）に基づいて生成できる。その後の超音波の送信／受信シーケンスが、複数の時点におけるせん断波の複数の画像を生成し得る。次いで、先に論じたように、せん断波の画像間の相関を計算して、せん断波伝搬速度の２乗を求める。或いは、音響放射力により生じた組織変位を特定して、せん断波伝搬速度の２乗を、変位の２次時間導関数と変位の２次空間導関数との比として算出する。
【００３３】
幾つかの実施形態において、検出したせん断波の伝搬速度の２乗（ｃ２）を表示し得る。有利なことに、式１に示したように、せん断波伝搬速度の２乗（ｃ２）をヤング率又はせん断弾性率に、伝搬速度（ｃ）よりも密接に関連付けることが可能である。従って、伝搬速度の２乗（ｃ２）は、実際の硬さに関する有効な代用値(proxy)を提供し得る。幾つかの実施形態において、伝搬速度の２乗（ｃ２）を３倍して、その値を表示し得る。組織密度が１ｇ／ｃｍ３に近い場合、この数値（すなわち、３ｃ２）が実際のヤング率に近い可能性がある。幾つかの実施形態において、任意の実数（ｂ）と伝搬速度の２乗（ｃ２）との積（ｂｃ２）を表示し得る。組織の密度は未知であり、推測に頼らざるを得ないため、実際の硬さの決定は難しく、誤りが生じやすい。
【００３４】
せん断波の伝搬特性（すなわち、速度の２乗ｃ２）をユーザに提示するために、カラーコーディング（色分け）技術、グレースケール技術、又は、グラフィックコーディング技術を用い得る。幾つかの実施形態において、組織内のせん断波の伝搬速度の２乗（ｃ２）が２次元カラー画像で表示される。幾つかの実施形態において、グラフィックコーディング画像及び／又は２次元画像も、伝搬速度の２乗ｃ２を表示するために用いられ得る。
【００３５】
せん断波伝搬速度の２乗ｃ２のうち低い値を、赤色を用いてコード化でき、高いｃ２値は青色を用いてコード化し得る。例えば、図６に示されている色凡例は、赤色の組織領域が、低いｃ２値（例えば、１ｍ２／ｓ２）に関連付けられたせん断波を含み、青色の組織領域が、高いｃ２値（例えば、１００ｍ２／ｓ２）に関連付けられたせん断波を含むことを示す。実施形態は、カラーベースのコード化に限定されない。組織内でのせん断波伝搬特性の画像を、グレースケール、又は図形パターン（例えば、鉛直線、水平線、斜交平行線、密度が異なるドットパターンなど）と色との任意の組合せを使用してコード化し得る。
【００３６】
伝搬速度の２乗（ｃ２）を決定した後、図６に示されているように、ｃ２を色波長に対して線形にコード化し得る。例えば、組織領域内のｃ２が５０ｍ２／ｓ２であると判断された場合、その組織領域を、黄色６３０を用いて表示し得る。
【００３７】
別の方法として、せん断波の伝搬速度の２乗（ｃ２）のカラーコーディングを、図７に示されているようにも定め得る。せん断波の伝搬速度の２乗のうち低い値に関連付けられた組織領域を青色７１０として表示し、速度の２乗のうち高い値に関連付けられた領域を赤色７２０として表示し得る。異なるカラーコーディング方法を用いてせん断波伝搬速度の２乗（ｃ２）を表してもよい。例えば、カラーコーディングが、色相、輝度、及び、その他の色特性に基づくものであってよい。カラーコード化されたスケールは、せん断波伝搬速度の２乗の、図６及び図７に示したものとは異なる最大値及び最小値を表し得る。これに関し、図６及び図７における速度の２乗の最大値１００ｍ２／ｓ２、及び、速度の２乗の最小値１ｍ２／ｓ２は例示のためのものに過ぎず、特許請求の範囲を限定するものではない。その他の値も、コード化スケールの最大値又は最小値を表し得る。
【００３８】
図１４に示されているように、赤色、緑色及び青色（ＲＧＢ）の値に基づいたカラーコーディングを用いて、せん断波の伝搬速度ｃ又は速度の２乗（ｃ２）を表してもよい。この例（図１４）において、組織内でのせん断波伝搬速度の２乗（ｃ２）は、ＲＧＢの値１４２０，１４３０，１４４０に基づいたカラーコーディングバー１４１０に従って示されている。この例において、せん断波伝搬速度の２乗（ｃ２）の値はカラーコーディングバー１４１０にて２５６色で表されるため、２５６個であり得る。速度の２乗（ｃ２）の最小値（０）１４１２は、Ｒ（０）１４２２と、Ｇ（０）１４３２と、Ｂ（０）１４４２との組み合わせから構成された色により表される。速度の２乗（ｃ２）の中間値（１２７）１４１５は、Ｒ（１２７）１４２５と、Ｇ（１２７）１４３５と、Ｂ（１２７）１４４５との組み合わせから構成された色で表される。速度の２乗（ｃ２）の最大値（２５５）１４１８は、Ｒ（２５５）１４２８と、Ｇ（２５５）１４３８と、Ｂ（２５５）１４４８との組み合わせから構成された色で表される。この例において、Ｒ（２５５）は、赤色指標２５５に相当する赤色を示しているに過ぎず、赤色値２５５が最高輝度の赤色であることを必ずしも示すわけではない。同様に、Ｇ（２５５）は、緑色指標２５５に対応する緑色を示し、Ｂ（２５５）は、青色指標２５５に相当する青色を示す。
【００３９】
別の方法として、赤色、緑色、青色及び黄色を用いてカラーコーディングバーを定義してもよい。或いは、色相に基づいたカラーコーディングバーを用いてもよい。
【００４０】
図１２は、ヒトの軟組織（例えば、乳房）内のせん断波の伝搬速度の２乗ｃ２を表示するカラーコード化画像１２６０の例を示す。カラーコーディングスケール１２５０が図示されており、このスケールにおいて、カラーコード１２１０（すなわち、赤色を表すコードであるが、この白黒文書では白として表示されている）が、せん断波伝搬速度の２乗の小さい値を表し、カラーコード１２２０（すなわち青色を表すコードであるが、この白黒文書では斜線で表示されている）は、せん断波伝搬速度の２乗の、より大きい値を表している。
【００４１】
コーディングスケール１２５０に基づくと、カラーコード化された画像１２６０が、伝搬速度の２乗ｃ２の大きい値の領域１２８０を含むことが分かる。せん断波伝搬速度の２乗ｃ２はヤング率に比例するため、領域１２８０に対応する組織領域は硬くなっている可能性がある。腫瘍は一般に硬いため、画像１２６０が病理的状態を示すと推測し得る。
【００４２】
カラーコーディング方法は、伝搬速度の２乗の大きい値を有するせん断波を含む領域と、伝搬速度の２乗の小さい値を有するせん断波を含むその他の領域との効率的な区別を提供する。従って、カラーコーディング方法は、軟組織領域内の硬い組織領域の効率的な識別を可能にする。せん断波の伝搬速度の２乗を表示する画像を、通常の超音波画像（例えば、Ｂモード画像）と組み合わせても、或いは、Ｂモード画像及びカラードプラ画像、並びに／又はスペクトルドプラ画像の組合せと組み合わせてもよい（たとえば、重畳により）。或いは、せん断波の伝搬速度の２乗を数値的に表示してもよい。幾つかの実施形態において、せん断波の伝搬速度の２乗を、グレースケールで表示しても、又は、その他のグラフィックコーディング方法（例えば、色ではないパターンを用いた方法）に基づいて表示してもよい。例えば、グレースケールコーディング方法を用いて、せん断伝搬速度の２乗の値のうち小さいものを黒色又は濃いグレーで表示し、せん断伝搬速度の２乗の値のうち大きいものを薄いグレー又は白色で表示し得る。
【００４３】
幾つかの実施形態において、第１の超音波パルスを生体組織に印加して、せん断波を生成する。第２の超音波パルスをこの生体組織に送信する。生体組織にて発生された１以上の超音波信号を生体組織から受信する。生体組織内のせん断波を、受信した１以上の超音波信号に基づいて検出する。生体組織内の複数の位置の各々におけるせん断波のピーク変位までの時間を決定する。こうして決定された、複数の位置の各々におけるせん断波のピーク変位までの時間に基づいて、せん断波伝搬速度の２乗を決定する。
【００４４】
幾つかの実施形態において、第１の超音波パルスを生体組織に印加して、せん断波を生成する。第２の超音波パルスを生体組織に送信する。第２の超音波に応答して生体組織にて発生された第１の超音波信号を生体組織から受信する。第３の超音波パルスを生体組織に送信する。第３の超音波パルスに応答して生体組織内にて発生された第２の超音波信号を生体組織から受信する。生体組織におけるせん断波を、受信した第１及び第２の超音波信号に基づいて検出する。生体組織内の複数の位置の各々におけるせん断波のピーク変位までの時間を決定する。こうして決定された、複数の位置の各々におけるせん断波のピーク変位までの時間に基づいて、せん断波伝搬速度の２乗を決定する。
【００４５】
幾つかの実施形態において、せん断波伝搬速度の２乗は、２つのポイントのみにてせん断波を検出することにより得られる。図１９は、せん断波の起源１９１０を示す。幾つかの実施形態によれば、２つのポイント１９３０と１９４０との間のせん断波の飛行時間（ｔ）が決定される。例えば、ポイント１９３０におけるピーク組織変位とポイント１９４０におけるピーク組織変位との間で経過した時間量を測定することにより、飛行時間（ｔ）を特定し得る。また、ポイント１９３０とポイント１９４０との間の距離（ｄ）を決定し得る。次いで、以下のように距離の２乗ｄ２を時間の２乗ｔ２で割算することにより、せん断波伝搬速度の２乗（ｃ２）を決定し得る。
【数２０】
2015512273.tif 000023
【００４６】
幾つかの実施形態において、せん断波伝搬速度の２乗を正確に求めることが可能であり、これは、せん断波を複数のポイントにて検出して、これらのポイントのうちの２つのポイント間の距離の飛行時間を決定することによる。図１９は、上述のように、せん断波の起源１９１０、並びに、複数のポイント１９１０，１９３０，１９４０，１９５０、及び１９６０を示す。せん断波１９２０は、図１９に示されているせん断波の起源１９１０における強い送信超音波パルスにより生成され得る。
【００４７】
第１の受信超音波信号と第２の受信超音波信号との相互相関ＳＡＤ、ＳＳＤ、ＳＣＤ又はＳＰＤを、先に論じたように、組織変位を追跡するために用い得る。第１の受信超音波信号は、せん断波が生成される前の基準信号として取得でき、第２の受信超音波信号は、せん断波が生成された後に取得され得る。組織変位は、式（６）における相互相関ＣＣ（ｔ，τ）を最大にする時間変位に対応する距離である。先に述べたように、式（７），（８），（９）及び（１０）におけるＳＡＤ、ＳＳＤ、ＳＣＤ又はＳＰＤに関して、ＳＡＤ、ＳＳＤ、ＳＣＤ又はＳＰＤを最小値にするウィンドウ変位に対応する距離が組織変位である。幾つかの実施形態において、基準信号を、せん断波が生成された後に取得してもよい。また、せん断波を、せん断波の起源１９１０にて検出及び監視することも可能である。せん断波の起源から各測定ポイントまでの距離を、図１９に示されているようにｄ０，ｄ１，ｄ２，ｄ３及びｄ４で示す。例えば、ｄ０はゼロである。
【００４８】
幾つかの実施形態において、カラードプラ速度、カラードプラシフト周波数、又はカラードプラ位相（式２、式３、式４及び式５に示したような）を、先に論じたように受信超音波信号のＩ−Ｑ信号に基づいて算出して組織変位を追跡し、そしてそれにより、ピーク変位までの時間を特定し得る。先に論じたように、カラードプラ速度（又は、カラードプラシフト周波数、カラードプラ位相）の時間積分により、組織変位が求められる。この場合、式２における角度θは、ゼロ（０）度とみなされ得る。
【００４９】
図１６は、せん断波による、時間経過（Ｘ軸）における組織変位の変化（Ｙ軸）を、組織内の複数のポイントにて観察したグラフを示す（各曲線に、それぞれ符号が付されている）。せん断波の起源から複数のポイントまでの、ピーク変位までの時間が図１６に示されている。図１７に示されているように、せん断波の起源から複数のポイントまでの距離（すなわちｄ０，ｄ１，ｄ２，ｄ３，ｄ４）を２乗し（ｄ２、Ｘ軸）、そして、ピーク変位までの時間の２乗（ｔ２、Ｙ軸）に対してプロットし、それにより、せん断波伝搬速度の２乗を測定する。これに関し、式１に示したように、せん断波伝搬速度の２乗は、組織の硬さ又はヤング率に、組織密度を介して直接関連付けられる。図１７に示されているように、測定精度を高めるために、距離の２乗と時間の２乗との回帰直線を得る。
【００５０】
回帰直線は、誤差の最小化により変数Ｘと変数Ｙとを密接に関連付ける線である。所与の画像ポイント（ピクセル）における変数Ｘ（Ｘｉは、実際のせん断データサンプル）と変数Ｙ（Ｙｉは、実際の引張データサンプル）とが回帰直線により関連付けられると考えられる。例えば、
【数２１】
2015512273.tif 000024
ａ及びｂは、誤差の合計（例えば誤差の２乗）を最小化するために、以下のように求められる。
【数２２】
2015512273.tif 000025
誤差は、回帰直線を用いた所与のＸｉに対して回帰直線から推定されたＹ値Ｙｅｓｔと実際のデータサンプルＹｉとの差である。これを以下の式に示す。
【数２３】
2015512273.tif 000026
【数２４】
2015512273.tif 000027
【００５１】
せん断波を２つのポイントのみにて検出し、そして、既知の距離ｄだけ離れているこれらの２つのポイント間の時間（ｔ）を測定することにより、せん断波伝搬速度の２乗を決定することができる。しかし、組織が均一で、且つ、組織全体を通じて同一のせん断特性（すなわち、せん断波伝搬速度の２乗）を有する場合、２つより多数のポイントを用いれば、測定精度が増大するであろう。従って、距離の２乗ｄ２と時間の２乗ｔ２との回帰直線を得て、これを用いて、せん断波伝搬速度の２乗ｃ２を決定し得る。回帰直線の勾配（又は勾配の逆数）が、せん断波伝搬速度の２乗を示し得る。
【００５２】
しかし、せん断波の起源１９１０においてさえも、せん断波生成の開始からピーク変位までの時間が正確に０にはならない場合があり、これにより、オフセット（ずれ）時間が生じる。従って、回帰直線は、図１７に示すように、グラフの起源（０，０）と交差しない場合があり、これにより、グラフの起源（０，０）からのずれを生じる。従って、回帰直線は、以下の式（ｔ２をＹ軸とし、ｄ２をＸ軸とする）に従うことになろう。
【数２５】
2015512273.tif 000028
又は  
【数２６】
2015512273.tif 000029
こうして、回帰直線勾配の逆数１／ａは、せん断波伝搬速度の２乗を示し得る。
【００５３】
幾つかの実施形態において、複数の位置の各々におけるピーク変位までの時間の２乗を算出する前に、せん断波起源におけるピーク変位まで時間を、複数の位置の各々におけるピーク変位までの時間から減算し得る。このプロセスは、上記の段落で論じたような、起源における回帰直線のオフセット（ｂ）を排除し得る。このプロセスが測定精度を向上させ得る。
【００５４】
組織は不均一であり得る。また、組織は、測定ポイントごとに、様々な硬さ、及び、せん断波伝搬速度の２乗に関する様々な値を呈し得る。図１８は、そのような不均一な組織の場合を示す。時間の２乗と距離の２乗との相関係数が計算される。
【００５５】
データサンプルＸｉとデータサンプルＹｉとの相関係数ｒが、以下のように得られる
【数２７】
2015512273.tif 000030
式中、
2015512273.tif 000031
及び
2015512273.tif 000032
は、それぞれ、Ｘｉ及びＹｉの平均値である。
【００５６】
相関係数が、予め設定された閾値よりも低い場合には、相関係数の算出のために用いられるデータ間では回帰直線が得られない。回帰直線は、適切に高い相関係数に関連したデータに対してのみ得られる。図１８において、右方に示されている１つのデータ値１８１０は、左方に示されている４つのデータとは異なる硬さ、及び、異なるせん断波伝搬速度の２乗を示し、その結果、図１７に示した相関係数よりも低い相関係数となっている。従って、回帰直線及び相関係数を得るために左方の４つのデータ値のみを用い、右方に示された１つのデータ値１８１０は除外し得る。これにより、せん断波伝搬速度の２乗の改善された測定を、幾つかの実施形態において達成し得る。
【００５７】
図２０は、幾つかの実施形態によるプロセスのフローチャートである。２０１０にてせん断波を生成する。２０２０にて、せん断波を複数の位置にて検出する。２０３０にて、これらの位置の各々に関してピーク変位までの時間を決定する。この測定されたピーク変位までの時間に基づいて、２０４０にて、せん断波伝搬速度の２乗を決定する。
【００５８】
図２１は、幾つかの実施形態によるプロセスのフローチャートである。２１１０にて、せん断波を生成する。２１２０にて、せん断波を複数の位置にて検出する。２１３０にて、位置の各々に関してピーク変位までの時間を決定する。２１４０にて、複数の位置の各々に関し、ピーク変位までの時間の２乗ｔ２を決定する。２１４０にて、測定位置の各々までの距離の２乗ｄ２を決定する。２１５０にて、複数の位置の各々に関してｄ２をｔ２で割算し、それにより、複数の位置の各々に関するせん断波伝搬速度の２乗ｃ２を得る。
【００５９】
幾つかの実施形態において、図２２のフローチャートの２２１０にてせん断波を生成する。２２２０にてせん断波を複数の位置にて検出する。位置の各々に関し、ピーク変位までの時間が測定される。複数の位置の各々に関してピーク変位までの時間の２乗ｔ２を決定する。２２５０にて、せん断波の起源から測定位置までの距離の２乗ｄ２を決定する。２２５０にて、決定された値ｔ２と値ｄ２との回帰直線及び相関係数を得る。相関係数が閾値よりも大きければ（２２６０）、２２８０にて回帰直線の勾配を計算することより、せん断波伝搬速度の２乗ｃ２が決定される。相関係数が閾値よりも大きくなければ（２２６０）、２２７０にて、データセットにおける、せん断波の起源から最も遠い位置にて得られたピーク変位までの時間が破棄される。フローは、２２５０に戻り、新しいデータセットに基づいた新しい回帰直線及び新しい相関係数を得て、上述のように継続する。それにより、起源から２番目に遠い位置にて得られたピーク変位までの時間が破棄される場合もあり、その場合は、また２２５０に戻る。
【００６０】
幾つかの実施形態において、せん断波起源付近の測定ポイント（例えば測定ポイント１９３０）を基準位置として用いることができ、そして、せん断波伝搬速度の２乗を得るための新しい演算を以下のように行い得る。
【００６１】
最初に、測定位置の１つを選択し得る。次に、選択した位置から複数の位置の各々までの距離を測定する。選択した位置から複数の位置の各々までの距離の２乗を計算する。このプロセスを、以下の段落において、より詳細に説明する。せん断波起源１９１０から測定ポイント１９３０までの距離ｄ１を距離から差し引いて、選択された位置からの新しい距離０を作成する。測定ポイント１９４０，１９５０及び１９６０の各々までの距離からも距離ｄ１を差し引く。同様に、測定ポイント１９３０（又は選択された位置）におけるピーク変位までの時間も、先に述べたように、測定ポイント１９３０，１９４０，１９５０，１９６０の各々におけるピーク変位までの時間から減算する。次いで、測定ポイント１９３０，１９４０，１９５０，１９６０の各々におけるピーク変位までの時間の２乗を計算する。また、選択された位置１９３０から測定ポイントの各々までの距離の２乗も計算する。幾つかの実施形態において、２１５０のプロセスと同様に、いずれの測定ポイントまでの距離の２乗も、同一測定ポイントにおけるピーク変位までの時間の２乗で割算し、これにより、せん断波伝搬速度の２乗を得る。
【００６２】
選択された位置にせん断波が到着する前に取得された測定データ（例えば、せん断波起源１９１０のデータ）は、現在説明しているプロセスには必要ないため、データセットから破棄し得る。幾つかの実施形態において、先に述べたように、回帰直線を、新しいデータセットに関して決定する。せん断波伝搬速度の２乗を求めるために、回帰直線の勾配の逆数を得る。幾つかの実施形態において、先に論じたように、回帰直線及び相関係数を、新しいデータセットに関して決定する。相関係数が、予め設定された閾値よりも高ければ、回帰直線勾配の逆数を求め、それにより、せん断波伝搬速度の２乗を得る。相関係数が予め設定された閾値よりも低い場合は、選択された位置から最も遠いデータを破棄する。次いで、新しい回帰直線及び新しい相関係数を決定する。このプロセスは、閾値よりも大きい相関係数が決定されるまで繰り返されることができ、その後、図２２のプロセスに示されているように、せん断波伝搬速度の２乗が決定される。
【００６３】
上記の説明において、せん断波の、１つの位置から別の位置への飛行時間を測定するためにピーク変位を用いた。飛行時間を測定するために、その他の変位特性（例えば、図１６に示されているような、変位カーブの先行の(leading)右上がり斜線）を用いてもよい。せん断波の、１つの位置から別の位置への飛行時間を測定するために、変位の最大時間導関数を用いてもよい。
【００６４】
上記の説明において、せん断波のピーク変位までの時間は、せん断波発生の開始から測定される。この時間は、せん断波を生成する超音波の送信時点から正確に測定され得る。幾つかの実施形態において、この時間は、超音波パルスがパルスの焦点又はせん断波の起源に到達した時点から測定され得る。幾つかの実施形態において、せん断波のピーク変位までの時間が任意の時点から測定され得る。
【００６５】
図３は、Ｂモードイメージング（撮像）、ドプラスペクトル、及び、カラードプライメージングを用いた慣用の超音波診断撮像システムの図である。このシステムは、その他のイメージングモード、例えば、弾性イメージング、３Ｄイメージング、リアルタイム３Ｄイメージング、組織ドプライメージング、組織ハーモニックイメージング、コントラストイメージングなどを含み得る。超音波信号が、送信機／送信ビームフォーマ３１０により駆動される超音波プローブ３３０から、送受信スイッチ３２０を介して送信される。プローブ３２０は、超音波振動子素子のアレイから成っていてよく、これらの素子は、送信機／送信ビームフォーマ３１０により、送信超音波ビームが集束及びステアリングされるように、異なる遅延時間を有して別々に駆動される。受信ビームフォーマ３４０が、プローブ３３０からスイッチ３２０を介して受信超音波信号を受信し、この信号３２５を処理する。受信ビームフォーマ３４０は、遅延及び／又は位相を信号に加え、これにより得られた信号が、受信された超音波ビームを集束及びステアリングするために合計される。受信ビームフォーマ３４０は、アポダイゼーション、増幅、及びフィルタリングを適用し得る。
【００６６】
処理された信号３４５は、ドプラスペクトルプロセッサ３５０、カラードプラプロセッサ３６０、及び、Ｂモード画像プロセッサ３７０に結合される。ドプラスペクトルプロセッサ３５０は、ドプラ信号プロセッサ及びスペクトル分析器を含み、ドプラ流速信号を処理して、ドプラスペクトル３５５を算出及び出力する。カラードプラプロセッサ３６０は、受信信号３４５を処理して、速度信号、パワー信号、及び分散信号３６５を算出及び出力する。Ｂモード画像プロセッサ３７０は、受信信号３４５を処理して、Ｂモード画像３７５、又は、振幅検出による信号の振幅を算出及び出力する。
【００６７】
ドプラスペクトル信号３５５、カラードプラプロセッサ信号（速度、出力、及び分散）３６５、及び、Ｂモードプロセッサ信号３７５は、これらの信号を走査変換信号に変換するスキャンコンバータ３８０に結合される。スキャンコンバータ３８０の出力は、超音波画像を表示するためのディスプレイモニタ３９０に結合される。
【００６８】
図２は、幾つかの実施形態によるせん断波プロセッサ２９５を含む超音波イメージングシステムの要素の図である。図２における超音波システムは、強い超音波パルスを生体組織に送信して音響放射力を生成し、この音響放射力が生体組織を押圧する。生体組織が押圧された後、せん断波が生成され、組織内で伝搬する。次いで、超音波システムは超音波パルスを送信及び受信し、それにより、生体組織内でせん断波が伝搬するときにこれらのせん断波を追跡する。複数の受信超音波ビームを、受信ビームフォーマ２４０により同時に形成し得る。同様に、複数の送信超音波ビームを、送信機／送信ビームフォーマ２１０により同時に形成し得る。先に述べたように、受信ビームフォーマ２４０から受信した超音波信号を処理して、組織変位、ドプラ速度、相関、及び、せん断波伝搬速度の２乗を求める。せん断波プロセッサ２９５が、先に述べたせん断波及び歪み処理方法を実行し得る。せん断波プロセッサ２９５は、受信ビームフォーマ２４０からの出力２４５を受信する。出力２９７は、せん断波速度の２乗のデータ、又は、その他のせん断波特性を含む。例えば、せん断波プロセッサ２９５が、せん断波伝搬速度の２乗をスキャンコンバータ２８０に出力すると、せん断波伝搬速度の２乗の表示が、Ｂモード画像、カラードプラ画像、又はスペクトルドプラ画像と共に、ディスプレイモニタに出力される。
【００６９】
せん断波プロセッサ２９５は、汎用の中央処理装置（ＣＰＵ）、デジタル信号プロセッサ（ＤＳＰ）、フィールドプログラマブルアレイ（ＦＰＧＡ）、画像処理ユニット（ＧＰＵ）、及び／又は別の（discreet）電子デバイスを備え得る。
【００７０】
図２は幾つかの実施形態による論理アーキテクチャを示す。実際の実装は、異なるように構成された、より多くの、又は異なる要素を含み得る。他のトポロジをその他の実施形態と組み合わせて用いてもよい。また、図２のシステムの各要素を、互いに通信する任意の個数の演算デバイスにより、任意の個数のその他の公衆及び／又は私設ネットワークを介して実装してもよい。このような演算デバイスの２つ以上を互いに離して配置することも、任意の既知の方式のネットワーク及び／又は専用接続を介して互いに通信させることも可能である。このシステムは、本明細書に記載した機能、及び、その他の任意の機能を提供するのに適したハードウェア及び／又はソフトウェア要素を含み得る。例えば、図２のシステムの実装にて用いられるいずれの演算デバイスも、本明細書に記載したように演算装置を動作させるプログラムコードを実行するプロセッサを含み得る。
【００７１】
本明細書中で論じた全てのシステム及びプロセッサは、１以上の一時的でないコンピュータ可読媒体に格納されたプログラムコードで実現され得る。このような媒体は、例えば、フロッピィディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、ブルーレイディスク、フラッシュドライブ、磁気テープ、及びソリッドステートランダムアクセスメモリ（ＲＡＭ）又はリードオンリメモリ（ＲＯＭ）記憶装置を含み得る。プログラムコードは、１以上の演算デバイスにより実行され、これにより、本文中に記載したプロセスのいずれもが実行され得る。これに関し、１以上の演算デバイスは、１以上のメモリデバイス（例えば、ＲＡＭ記憶装置又はＲＯＭ記憶装置）も含むことができ、これにより、プログラムコードの実行前に、これらのプログラムコードを記憶し得る。従って、実施形態は、ハードウェアとソフトウェアとの何ら特定の組み合わせにも限定されない。
【００７２】
１以上の実施形態について説明してきた。しかし、当業者には、様々な変更が明らかであろう。
【符号の説明】
【００７３】
１１０  超音波振動子アレイ
１２０  超音波パルス
１３０  焦点
１４０  せん断波
１５０  せん断波
１６０  生体組織
２４５  出力
３４５  受信信号
４１０  超音波振動子アレイ
４２０  超音波パルス
４３０  超音波パルス
４８０  生体組織
４４０  超音波受信信号
４４５  超音波受信信号
４６０  超音波受信信号
４６５  超音波受信信号
４７０  超音波受信信号
５２０  送信超音波ビーム
５４０  超音波信号
５４５  超音波信号
５５０  超音波信号
８１０  超音波振動子
８２０  超音波パルス
８３０  組織移動
８４０  せん断波
８５０  せん断波
８６０  生体組織
９１０  超音波振動子
９２０  組織の軸方向移動
９６０  生体組織
１０４０  せん断波
１０５０  せん断波
１０６０  せん断波
１０７０  せん断波
１１４０  せん断波
１１５０  せん断波
１１６０  せん断波
１０８０  局所せん断波伝搬
１１７０  せん断波
１４１０  カラーコーディングバー
１９１０  せん断波の起源
ｄ１  せん断波の起源から測定ポイントまでの距離
ｄ２  せん断波の起源から測定ポイントまでの距離
ｄ３  せん断波の起源から測定ポイントまでの距離
ｄ４  せん断波の起源から測定ポイントまでの距離
【図１】
2015512273.tif 000033
【図２】
2015512273.tif 000034
【図３】
2015512273.tif 000035
【図４】
2015512273.tif 000036
【図５】
2015512273.tif 000037
【図６】
2015512273.tif 000038
【図７】
2015512273.tif 000039
【図８】
2015512273.tif 000040
【図９】
2015512273.tif 000041
【図１０】
2015512273.tif 000042
【図１１】
2015512273.tif 000043
【図１２】
2015512273.tif 000044
【図１３】
2015512273.tif 000045
【図１４】
2015512273.tif 000046
【図１５】
2015512273.tif 000047
【図１６】
2015512273.tif 000048
【図１７】
2015512273.tif 000049
【図１８】
2015512273.tif 000050
【図１９】
2015512273.tif 000051
【図２０】
2015512273.tif 000052
【図２１】
2015512273.tif 000053
【図２２】
2015512273.tif 000054
【国際調査報告】
2015512273.tif 000055