JP6390516B2

JP6390516B2 - 超音波診断装置及び超音波診断装置の制御方法

Info

Publication number: JP6390516B2
Application number: JP2015107412A
Authority: JP
Inventors: 文平田路
Original assignee: Konica Minolta Inc
Current assignee: Konica Minolta Inc
Priority date: 2015-05-27
Filing date: 2015-05-27
Publication date: 2018-09-19
Anticipated expiration: 2035-05-27
Also published as: US20160345939A1; US10722215B2; JP2016220736A

Description

本発明は、超音波診断装置における画像化処理に関し、特に、せん断波を用いた組織の硬さ測定に関する。

近年、被検体内の組織の硬さを評価する機能を有する超音波診断装置が普及してきている。超音波診断装置を用いて、硬さを評価する方法は大きく分けて２つある。１つは、超音波探触子を用いて被検体内の組織を体表から圧迫した後に解放し、圧迫に対する被検体内の組織の歪みの大きさから、被検体内の組織の相対的な硬さを評価する方法である。この方法では、組織が周辺組織と比べて硬いか軟らかいかという、被検体内の相対的な硬さを評価できる。

もう１つは、被検体内の着目領域（ＲＯＩ；ＲｅｇｉｏｎｏｆＩｎｔｅｒｅｓｔ）にせん断波（ｓｈｅａｒｗａｖｅ）を生じさせ、着目領域内の組織の変位を時系列に取得することで、せん断波の伝播速度を計測する方法である。せん断波の伝播速度は組織の弾性率に応じて変化するため、この方法では組織の絶対的な硬さ（例えば、弾性率）を評価することができる。せん断波を生じさせる手法としては、例えば、ＡＲＦＩ（ＡｃｏｕｓｔｉｃＲａｄｉａｔｉｏｎＦｏｒｃｅＩｍｐｕｌｓｅ）と呼ばれる、プッシュパルスを焦点に集中させることで、超音波であるプッシュパルスの音圧により、焦点上にある被検体内の組織に変位を発生させる手法が用いられる。

特開２０１５−２３９１３号公報

しかしながら、せん断波を用いる方法では、せん断波による組織の変位が小さいと、せん断波の速度検出の精度が低下するため、組織の硬さの測定精度が低下する。例えばプッシュパルスを送信した先に硬い組織があったときなどは、せん断波の減衰や反射などが発生する場合があり、せん断波の強度が組織の変位を生じさせるのに不十分であったり届き難い領域が生じたりする場合がある。一方で、超音波探触子の特性上、および、組織に与える影響の観点から、プッシュパルスの音圧には上限がある。そのため、測定対象の構造により測定精度が低下してしまう場合があった。被検体内の組織の硬さの測定精度が低下し、または、測定ができなかった場合に、検査者は再測定のため、原因を推測し、それに基づき送信プロファイルの調整を試行錯誤する必要があった。

本開示は、上記課題に鑑みてなされたものであり、検査者にとって簡便に、組織の硬さ測定の精度を向上させることを目的とする。

本開示の一態様に係る超音波診断装置は、超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後超音波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の、前記被検体内の着目領域における伝播状態を検出する超音波診断装置であって、送信プロファイルに基づいてプッシュパルスを送信するプッシュパルス送信部と、前記プッシュパルスに続き超音波を被検体へ複数回送信し、送信された超音波に対応する被検体からの反射超音波を受信して複数の受信信号を時系列に取得し、前記着目領域内の組織の、前記プッシュパルスに起因するせん断波による変位をそれぞれ検出する変位検出部と、前記変位検出部による検出結果に基づいて、前記着目領域におけるせん断波の伝播の有無を評価する評価部と、せん断波が伝播していない領域が前記着目領域において占める位置に基づいて、前記着目領域において、前記せん断波が伝播していない領域より深い領域と浅い領域との双方に前記せん断波が伝播している場合、前記超音波探触子の素子配列方向にプッシュパルスが集中する位置を移動させることで、当該せん断波が伝播していない領域にせん断波が到達するように送信プロファイルを調整するプッシュパルス調整部とを備え、前記プッシュパルス調整部が送信プロファイルを調整した場合に、前記プッシュパルス送信部は、調整された後の送信プロファイルに基づいて新たなプッシュパルスを送信し、前記変位検出部は、前記新たなプッシュパルスに起因するせん断波による変位を検出することを特徴とする。

本開示によれば、上記構成により、せん断波の伝播状態に基づいてプッシュパルスの押圧位置やエネルギーを制御する送信プロファイルを最適化できる。そのため、せん断波の反射や屈折、減衰といった原因により被検体内の組織の硬さの測定精度が低下し、または、測定ができなかった場合に、送信プロファイルを最適化して再測定し、その結果を得ることができる。したがって、被検体内の組織の硬さの測定精度が低下し、または、測定ができなかった場合に、検査者が再測定のため、原因を推測し、それに基づき送信プロファイルの調整を試行錯誤する煩雑さを低減することができる。これにより、検査者にとって簡便な手段で高精度に組織の硬さ測定を行うことができる。

実施の形態に係る超音波診断装置１のブロック図である。実施の形態に係る超音波診断装置１全体の動作を示すフローチャートである。せん断波の発生と伝播の様子を示す模式図である。実施の形態に係るせん断波の伝播解析の動作を示す模式図である。実施の形態に係る弾性率測定の処理を示す模式図である。（ａ）は、実施の形態に係る伝播範囲図の例を示す模式図である。（ｂ）は、実施の形態に係るせん断波の速度分布図の例を示す模式図である。実施の形態に係るプッシュパルスの送信プロファイル調整の動作を示すフローチャートである。（ａ）は、せん断波の非伝播領域と送信プロファイル調整内容の一例を示す模式図である。（ｂ）は、せん断波の非伝播領域と送信プロファイル調整内容の一例を示す模式図である。（ｃ）は、せん断波の非伝播領域の一例を示す模式図である。（ｄ）は、図８（ｃ）の伝播範囲図に対する送信プロファイル調整内容の一例を示す模式図である。（ｅ）は、せん断波の非伝播領域の一例を示す模式図である。（ｆ）は、図８（ｅ）の伝播範囲図に対する送信プロファイル調整内容の一例を示す模式図である。（ｇ）は、せん断波の非伝播領域と送信プロファイル調整内容の一例を示す模式図である。（ｈ）は、せん断波の速度分布と送信プロファイル調整内容の一例を示す模式図である。（ａ）は、変形例１に係るせん断波の非伝播領域と送信プロファイル調整内容の一例を示す模式図である。（ｂ）は、図９（ａ）に対応するせん断波の速度分布と送信プロファイル調整内容の一例を示す模式図である。変形例３に係るせん断波の伝播解析の動作を示す模式図である。変形例４に係る超音波診断装置全体の動作を示すフローチャートである。変形例４に係る測定中の状態を示す表示例である。

≪発明を実施するための形態に至った経緯≫
発明者は、せん断波を用いて組織の硬さの測定を行う超音波診断装置において、測定精度を向上させるために各種の検討を行った。
上述したように、せん断波の伝播速度は組織の弾性率に応じて変化する。そのため、組織の硬さの測定精度を向上させるには、せん断波の伝播速度を精度よく算出する必要がある。そのためには、せん断波による組織の変位が大きく、せん断波が減衰せず着目領域全体に行きわたることが好ましい。一方で、被検体内の組織に必要以上に変位を生じさせることは、組織への影響の観点から好ましくない。したがって、被検体に負担をかけずに組織の硬さの測定精度を向上させるためには、組織に対して悪影響を与えない範囲内で変位を最大化することが好ましい。

また、着目領域に硬さの異なる複数の組織があるなど組織の硬さが不均一である場合、硬さが変化する箇所、例えば、組織の境界において、せん断波の反射や屈折が起きうる。このような場合、せん断波の進行方向が変化したり、進行波と反射波とが干渉したりすることで、せん断波の速度が算出できない、または、精度が低下することがある。
このような問題への対処は、組織の硬さの測定を行う前に行うことが難しい。せん断波の減衰、反射、屈折などは組織の硬さそのものに依存するからである。例えば、特許文献１に開示されているように、超音波画像（Ｂモード画像）に基づいて着目領域が適切であるか否かの判定を行うことはできるが、プッシュパルスの送信プロファイルの調整は、超音波画像（Ｂモード画像）や一度行った組織の硬さ測定の結果に基づいて、検査者が手探りで行う必要がある。

そこで、発明者は、上記課題に鑑み、プッシュパルスの送信プロファイルを最適化し組織の硬さ測定の精度を向上させる技術について検討を行い、実施の形態に係る超音波診断装置に相当するに至ったものである。
以下、実施の形態に係る超音波診断装置について図面を用いて詳細に説明する。
≪実施の形態≫
実施の形態に係る超音波診断装置１のブロック図を図１に示す。超音波診断装置１は、制御部１１、せん断波励起部１２、超音波信号取得部１３、変位検出部１４、伝播解析部１５、評価部１６、プッシュパルス調整部１７、断層画像記憶部１８、変位量記憶部１９、波面画像記憶部２０、弾性画像記憶部２１を備える。また、制御部１１には、超音波探触子２と表示部３とがそれぞれ接続可能に構成されている。図１は超音波診断装置１に超音波探触子２、表示部３が接続された状態を示している。

超音波探触子２は、例えば、一次元方向に配列された複数の振動子（不図示）を有する。各振動子は、例えば、ＰＺＴ（チタン酸ジルコン酸鉛）により構成される。超音波探触子２は、せん断波励起部１２で生成された電気信号（以下、「ＡＲＦＩ駆動信号」と呼ぶ）、または、超音波信号取得部１３で生成された電気信号（以下、「送信駆動信号」と呼ぶ）を制御部１１から受け取り、超音波に変換する。超音波探触子２は、超音波探触子２の振動子側外表面を被検体の皮膚表面など表面に接触させた状態で、ＡＲＦＩ駆動信号又は送信駆動信号から変換され、複数の振動子から発せられる複数の超音波からなる超音波ビームを被検体内の測定対象に向けて送信する。そして、超音波探触子２は、送信駆動信号に基づく送信超音波に対する、測定対象からの複数の反射超音波を受信し、複数の振動子によりこれら反射超音波をそれぞれ電気信号（以下、「素子受信信号」と呼ぶ）に変換し、制御部１１を介して素子受信信号を超音波信号取得部１３に供給する。なお、ここではせん断波励起部１２及び超音波信号取得部１３を別構成として説明しているが、超音波信号取得部１３の送信駆動信号を生成する構成と同一の構成により、ＡＲＦＩ駆動信号を生成してもよい。

せん断波励起部１２は、プッシュパルスを超音波探触子２に送出させるための電気信号であるＡＲＦＩ駆動信号を生成する。プッシュパルスとは、被検体内にせん断波を発生させるため、被検体内の組織に変位をもたらすためのパルス状の超音波である。具体的には、被検体内の着目領域内のある１点を焦点とする、波数が後述する送信超音波より多い超音波である。したがって、ＡＲＦＩ駆動信号は、超音波探触子２を構成する各振動素子から送出される超音波が焦点に届くように生成される、パルス状の電気信号である。せん断波励起部１２は、制御部１１から、プッシュパルスの連続送信回数、および各プッシュパルスの焦点位置、送信に用いる振動子、波数または送信時間長などを規定した送信プロファイルを受け取り、送信プロファイルに基づいて１以上のＡＲＦＩ駆動信号を生成する。

超音波信号取得部１３は、送信超音波を超音波探触子２に送出させるための電気信号である送信駆動信号を生成する。送信駆動信号は、超音波探触子２を構成する各振動素子から送出される送信超音波が送信フォーカス点に同時に届くように生成される、振動素子ごとにタイミングの異なるパルス状の電気信号である。また、超音波信号取得部１３は、反射超音波に基づく素子受信信号に整相加算を行って音響線信号を生成する。超音波信号取得部１３は、生成した音響線信号を、制御部１１を介して断層画像記憶部１８に出力する。

変位検出部１４は、変位検出の対象となる１つの断層画像に係る複数の音響線信号（以下、「断層画像信号」と呼ぶ）と、基準となる１つの断層画像に係る複数の音響線信号（以下、「基準断層画像信号」と呼ぶ）とを、制御部１１を介して断層画像記憶部１８から取得する。基準断層画像信号とは、断層画像信号からせん断波による変位を抽出するために用いるものであり、具体的には、プッシュパルス送出前に着目領域を撮像した断層画像信号である。そして、変位検出部１４は、断層画像信号と基準断層画像信号との差分から、断層画像信号の各画素の変位を検出し、変位を各画素の座標と関連付けて変位画像を生成する。変位検出部１４は、生成した変位画像を、制御部１１を介して変位量記憶部１９に出力する。

伝播解析部１５は、変位画像を、制御部１１を介して変位量記憶部１９から取得する。伝播解析部１５は、変位画像から、変位画像を取得した各時刻における、せん断波の波面の位置、進行方向および速度を検出し、波面画像を生成する。伝播解析部１５は、せん断波の波面の位置、進行方向および速度から、変位画像の各画素に対応する被検体組織の弾性率を算出し、弾性画像を生成する。伝播解析部１５は、生成した波面画像を波面画像記憶部２０に、弾性画像を弾性画像記憶部２１に、制御部１１を介してそれぞれ出力する。

評価部１６は、制御部１１を介して波面画像を波面画像記憶部２０から取得する。評価部１６は、せん断波の波面の位置、進行方向および速度を基に、例えば、着目領域におけるせん断波の伝播の有無、せん断波の伝播の一様性、せん断波の反射や屈折の有無を評価する。評価部１６は、評価結果を制御部１１に出力する。詳細は後述する。
プッシュパルス調整部１７は、制御部１１から評価部１６が出力した評価結果を取得する。プッシュパルス調整部１７は、評価結果を基に、プッシュパルスの送信プロファイルを最適化し、最適化した後の送信プロファイルを制御部１１に出力する。詳細は後述する。

制御部１１は、上述したような各構成要素の制御に加え、伝播解析部１５が生成した弾性画像を表示部３に出力する。
断層画像記憶部１８、変位量記憶部１９、波面画像記憶部２０、弾性画像記憶部２１は、それぞれ、断層画像、変位画像、波面画像、弾性画像データを記憶する。断層画像記憶部１８、変位量記憶部１９、波面画像記憶部２０、弾性画像記憶部２１のそれぞれは、例えば、ＲＡＭ、フラッシュメモリ、ハードディスク、光学ディスクなどの記憶媒体で実現される。なお、断層画像記憶部１８、変位量記憶部１９、波面画像記憶部２０、弾性画像記憶部２１のうち２以上を単一の記憶媒体で実現してもよい。また、断層画像記憶部１８、変位量記憶部１９、波面画像記憶部２０、弾性画像記憶部２１のうち１以上は、超音波診断装置１の外部に構成され、インターフェースを介して超音波診断装置１と接続されてもよいし、超音波診断装置１からネットワークを介してアクセス可能に構成された資源、例えば、ファイルサーバやＮＡＳ（ＮｅｔｗｏｒｋＡｔｔａｃｈｅｄＳｔｏｒａｇｅ）であってもよい。

制御部１１、せん断波励起部１２、超音波信号取得部１３、変位検出部１４、伝播解析部１５、評価部１６、プッシュパルス調整部１７のそれぞれは、例えば、ＦＰＧＡ（ＦｉｅｌｄＰｒｏｇｒａｍｍａｂｌｅＧａｔｅＡｒｒａｙ）、ＡＳＩＣ（ＡｐｐｌｉｃａｔｉｏｎＳｐｅｃｉｆｉｃＩｎｇｅｇｒａｔｅｄＣｉｒｃｕｉｔ）などのハードウェアにより実現される。なお、これらの一部または全部は、単一のＦＰＧＡ、または、ＡＳＩＣで実現されてもよい。また、これらは、それぞれ個別に、または、２以上を１まとめとして、メモリと、ＣＰＵ（ＣｅｎｔｒａｌＰｒｏｃｅｓｓｉｎｇＵｎｉｔ）、ＧＰＵ（ＧｒａｐｈｉｃＰｒｏｃｅｓｓｉｎｇＵｎｉｔ）などのプログラマブルデバイスとソフトウェアで実現されてもよい。

≪動作≫
実施の形態に係る超音波診断装置１の動作について説明する。図２は、超音波診断装置１全体の動作を示すフローチャートである。
まず、制御部１１が、被検体内に着目領域を設定する（ステップＳ１０）。着目領域を設定する方法は、例えば、表示部３に断層画像記憶部１８に記録されている最新の断層画像を表示し、タッチパネル、マウス、トラックボールなどの入力部（図示しない）を通して検査者に着目領域を指定させる。なお、着目領域の設定方法はこの場合に限られず、例えば、断層画像の全域を着目領域としてもよいし、あるいは、断層画像の中央部分を含む一定範囲を着目領域としてもよい。また、着目領域を設定する際に、断層画像を取得してもよい。

次に、着目領域に超音波を送受信し、取得した受信信号を記憶する（ステップＳ２０）。具体的には、次のような動作となる。まず、以下のように送信イベントを行う。最初に、超音波信号取得部１３が、パルス状の送信信号を生成する。次に、超音波信号取得部１３は、送信信号に対して、超音波探触子２の素子毎の遅延時間を設定する送信ビームフォーミングを行い、超音波探触子２の各素子に対応する複数の送信駆動信号を生成する。超音波探触子２の各振動子は対応する送信駆動信号を超音波に変換することで、超音波ビームが被検体内に送出される。次に、超音波探触子２の各振動子は、被検体内から反射された反射超音波を取得して素子受信信号に変換する。超音波信号取得部１３は、素子受信信号に整相加算を行い、音響線信号を生成する。制御部１１は、送信イベントごとに超音波信号取得部１３から音響線信号を取得し、１つの断層画像を構成する複数の音響線信号を、断層画像信号として、断層画像記憶部１８に記憶する。

次に、プッシュパルスを送信する（ステップＳ３０）。具体的には、せん断波励起部１２が、制御部１１に保持されている送信プロファイルを用いて、パルス状のＡＲＦＩ信号を生成する。次に、せん断波励起部１２は、ＡＲＦＩ信号に対して、超音波探触子２の素子毎の遅延時間を設定する送信ビームフォーミングを行い、超音波探触子２の各素子に対応する複数のＡＲＦＩ駆動信号を生成する。ここで、プッシュパルスを始めて送信する場合は、例えば、着目領域の中央を焦点位置とし、超音波探触子２の振動子列の中心を含む所定数の素子を用いて、１回だけ所定の時間長のプッシュパルスを送信するような初期送信プロファイルを用いる。超音波探触子２の各振動子は対応する送信駆動信号を超音波に変換することで、プッシュパルスが被検体内に送出される。なお、２以上のプッシュパルスを送信する送信プロファイルを用いる場合、それぞれのプッシュパルスを連続して送信する。

ここで、プッシュパルスによるせん断波の生成について、図３（ａ）〜（ｅ）の模式図を用いて説明する。図３（ａ）は、着目領域に対応した被検体内領域の、プッシュパルス印加前における組織を示した模式図である。図３（ａ）〜（ｅ）において、個々の“○”は、着目領域における被検体内の組織の一部を、破線の交点は、負荷がない場合の組織”○“の中心位置を、それぞれ示している。ここで、超音波探触子２を皮膚表面１００に密接させた状態で焦点１０１に対してプッシュパルスを印加すると、図３（ｂ）の模式図に示すように、焦点１０１に位置していた組織１３２が、プッシュパルスの進行方向に押されて移動する。また、組織１３２からプッシュパルスの進行方向側にある組織１３３は、組織１３２に押されてプッシュパルスの進行方向に移動する。次に、プッシュパルスの送信が終了すると、組織１３２、１３３が元の位置に復元しようとするので、図３（ｃ）の模式図に示すように、組織１３１〜１３３がプッシュパルスの進行方向に沿った振動を開始する。すると、図３（ｄ）の模式図に示すように、振動が組織１３１〜１３３に隣接する、組織１２１〜１２３および組織１４１〜１４３に伝播する。さらに、図３（ｅ）の模式図に示すように、振動がさらに組織１１１〜１１３および組織１５１〜１５３に伝播する。したがって、被検体内において、振動が振動の方向と直交する向きに伝播する。すなわち、せん断波がプッシュパルスの印加場所に発生し、被検体内を伝播する。

図２に戻って説明を続ける。次に、着目領域に超音波を複数回送受信し、取得した複数の超音波信号を記憶する（ステップＳ４０）。具体的には、プッシュパルスの送信終了の直後から、例えば、秒間１万回、ステップＳ２０と同じ動作を繰り返し行う。これにより、せん断波の発生直後から伝播が終わるまでの間、被検体の断層画像を繰り返し取得する。

次に、各画素の変位を検出する（ステップＳ５０）。具体的には、変位検出部１４が、ステップＳ２０で断層画像記憶部１８に記憶された断層画像信号を、基準断層信号として取得する。上述したように、基準断層画像信号は、プッシュパルスの送出前、すなわち、せん断波の発生前に取得された断層画像信号である。次に、変位検出部１４は、ステップＳ４０で断層画像記憶部１８に記憶された各断層画像信号に対し、基準断層画像信号との差分から、当該断層画像信号が取得された時刻における、各画素の変位を検出する。具体的には、例えば、断層画像信号を８ピクセル×８ピクセルなどの所定の大きさの領域に分割し、各領域と基準断層画像信号とをパターンマッチングすることで、断層画像信号の各画素の変位を検出する。パターンマッチングの方法としては、例えば、各領域と基準断層画像信号内の同サイズの基準領域との間で、対応する画素毎に輝度値の差分を算出してその絶対値の合計値を算出し、その合計値が最も小さくなる領域と基準領域との組み合わせについて、領域と基準領域とが同一の領域であるものとし、領域の基準点（例えば、左上の角）と基準領域の基準点との距離を変位として検出する。なお、領域のサイズは８ピクセル×８ピクセル以外であってもよいし、輝度値の差分の絶対値の合計値に替えて、例えば、輝度値の差分の２乗の合計値を用いてもよい。また、変位として、領域の基準点と基準領域の基準点とのｙ座標の差（深さの差）を算出してもよい。これにより、各断層画像信号の各画素に対応する被検体の組織が、プッシュパルスまたはせん断波によってどれだけ動いたかが変位として算出される。なお、変位の検出方法はパターンマッチングに限られず、例えば、断層画像信号と基準断層画像信号との相関処理など、２つの断層画像信号間の動き量を検出する任意の技術を用いてよい。変位検出部１４は、１の断層画像に係る各画素の変位を当該画素の座標と対応付けることで変位画像を生成し、生成した変位画像を変位量記憶部１９に出力する。

次に、せん断波の伝播解析を行う（ステップＳ６０）。具体的には、各変位画像からせん断波の波面を抽出して波面画像を生成する。この波面画像より、波面の位置、振幅、進行方向および速度を容易に検出することができる。波面画像の生成は、例えば、変位領域の抽出、細線化処理、空間フィルタリング、時間フィルタリングの手順によって行う。
図４を用いて具体的な処理を説明する。図４（ａ）は、変位画像の一例を示している。図３と同じく、図中の“○”は着目領域における被検体内の組織の一部を示しており、プッシュパルスを印加する前の位置は破線の交点である。また、ｘ軸は超音波探触子２における素子の並ぶ方向、ｙ軸は、被検体の深さ方向である。伝播解析部１５は、ｙ座標ごとに変位量δを座標ｘの関数として、動的閾値を用いることで変位量δが大きい領域を抽出する。また、ｘ座標ごとに変位量δを座標ｙの関数として、動的閾値を用いて、ある閾値を超える領域を変位量δが大きい領域として抽出する。動的閾値とは、対象領域内について信号解析又は画像解析を行って閾値を決定することである。閾値は一定値ではなく、対象領域の信号の幅や最大値などによって異なる値となる。図４（ａ）に、ｙ＝ｙ₁の直線２１０上における変位量をプロットしたグラフ２１１と、ｘ＝ｘ₁の直線２２０上における変位量をプロットしたグラフ２２１とを示す。これにより、例えば、変位量δが閾値より大きな変位領域２３０が抽出できる。

次に、伝播解析部１５は、変位領域に細線化処理をおこなって波面を抽出する。図４（ｂ）の模式図に示している変位領域２４０、２５０は、それぞれ、ステップＳ５２において変位領域として抽出された領域である。伝播解析部１５は、例えば、Ｈｉｌｄｉｔｃｈの細線化アルゴリズムを用いて、波面を抽出する。例えば、図４（ｂ）の模式図において、変位領域２４０から波面２４１が、変位領域２５０から波面２５１が、それぞれ抽出される。なお、細線化のアルゴリズムはＨｉｌｄｉｔｃｈに限らず、任意の細線化アルゴリズムを用いてよい。また、各変位領域に対して、変位量δが閾値以下の座標を変位領域から取り除く処理を、変位領域が幅１ピクセルの線になるまで、閾値を大きくしながら繰り返し行ってもよい。

次に、伝播解析部１５は、細線化処理後の波面画像データに対して空間フィルタリングを行い、長さが短い波面を除去する。例えば、ステップＳ５３で抽出した各波面の長さを検出し、全ての波面の長さの平均値の１／２よりも長さが短い波面を、ノイズとして削除する。具体的には、図４（ｃ）の波面画像に示すように、波面２６１〜２６４の長さの平均値を算出し、それよりも短い波面２６３、２６４を、ノイズとして消去する。これにより、誤検出された波面を消去できる。

伝播解析部１５は、変位領域の抽出、細線化処理、空間フィルタリングの動作を、全ての変位画像に対して行う。これにより、変位画像に対して１対１で波面画像データが生成される。
最後に、伝播解析部１５は、複数の波面画像データに対して時間フィルタリングを行い、伝播していない波面を除去する。具体的には、時間的に連続する２以上の波面画像において、波面位置の時間変化を検出し、速度が異常である波面をノイズとして除去する。伝播解析部１５は、例えば、時刻ｔ＝ｔ₁の波面画像２７０、時刻ｔ＝ｔ₁＋Δｔの波面画像２８０、時刻ｔ＝ｔ₁＋２Δｔの波面画像２９０との間で、波面位置の時間変化を検出する。例えば、波面２７１に対して、波面画像２８０のうち、波面２７１と同じ位置を中心に、波面と垂直な向き（図４においてはｘ軸方向）にΔｔの間にせん断波が移動しうる領域２７６で、波面２７１との相関処理を行う。このとき、波面２７１のｘ軸の正方向（図の右側）と負方向（図の左側）の双方を含む範囲内で相関処理を行う。これは、透過波と反射波の両方を検出するためである。これにより、波面２７１の移動先が波面画像２８０内の波面２８１であると検出し、時間Δｔにおける波面２７１の移動距離を算出する。同様に、波面２７２、２７３のそれぞれについて、波面画像２８０において当該波面と同じ位置を中心に、波面と垂直な向きにΔｔの間にせん断波が移動しうる領域で相関処理を行う。これにより、波面２７２が波面２８３の位置に、波面２７３が波面２８２の位置に、それぞれ移動したことを検出する。波面画像２８０と波面画像２９０との間でも同様の処理を行い、波面２８１が波面２９１の位置に、波面２８２が波面２９２の位置に、波面２８３が波面２９３の位置に、それぞれ移動したことを検出する。ここで、波面２７３、波面２８２、波面２９２で示される１の波面については、他の波面と比べて移動距離が著しく小さい（伝播速度が著しく遅い）。このような波面は誤検知である可能性が高いので、ノイズとして消去する。これにより、図４（ｅ）の波面画像３００に示すように、波面３０１、３０２が検出できる。

伝播解析部１５は、生成した時刻ごとの波面画像データと、波面の対応情報とを波面画像記憶部２０に出力する。ここで、波面の対応情報とは、同一の波面が各波面画像のどの波面に対応するかを示した情報であり、例えば、波面２７２が波面２８３の位置に移動したことが検出された場合、波面２８３と波面２７２とが同一の波面であるという情報である。

次に、弾性画像を生成して表示する（ステップＳ７０）。伝播解析部１５は、弾性画像を生成する。具体的には、時刻ごとの波面画像データと、波面の対応情報とから、各時刻における波面の位置、速度を検出する。さらに、変位画像と断層画像との関係から、断層画像の各画素について複数の変位画像におけるせん断波の最大速度から弾性率を算出し、断層画像の各画素と弾性率とを対応付けて弾性画像を生成する。図５を用いて弾性画像の生成について説明する。図５（ａ）は、ある時刻ｔにおける波面画像と、時刻ｔ＋Δｔにおける波面画像を１つの波面画像３１０として合成したものである。ここで、時刻ｔにおける波面３１１と、時刻ｔ＋Δｔにおける波面３１２とが同一の波面であるとする対応情報が存在するものとする。伝播解析部１５は、対応情報から、波面３１１上の座標（ｘ_t、ｙ_t）に対応する波面３１２上の座標（ｘ_t+Δt、ｙ_t+Δt）を検出する。これにより、時刻ｔに座標（ｘ_t、ｙ_t）を通過したせん断波が、時刻ｔ＋Δｔに座標（ｘ_t+Δt、ｙ_t+Δt）に到達していると推定できる。したがって、座標（ｘ_t、ｙ_t）を通過したせん断波の速度ｖ（ｘ_t、ｙ_t）は、座標（ｘ_t、ｙ_t）と座標（ｘ_t+Δt、ｙ_t+Δt）との間の距離ｍを所要時間Δｔで割った値と推定できる。すなわち、ｖ（ｘ_t、ｙ_t）＝ｍ／Δｔ＝√｛（ｘ_t+Δt−ｘ_t）²＋（ｙ_t+Δt−ｙ_t）²｝／Δｔとなる。伝播解析部１５は、全ての波面に対して上述の処理を行い、波面が通過した全座標についてせん断波の速度を取得し、せん断波の速度を基に、各座標における弾性率を算出する。伝播解析部１５は、弾性率に基づいて、色情報をマッピングした弾性画像を生成する。例えば、図５（ｂ）のように、弾性率が一定値以上の座標は赤、弾性率が一定値未満の座標は緑、弾性率が取得できなかった座標は黒、というように色分けした弾性画像３２０を生成する。分類は二値化に限られず、所定の段階で分類及び色分けを行ってよい。図５（ｂ）において、領域３２２は弾性率が一定値以上の領域であり、内包物３２１に対応する。なお、図５（ｂ）では説明のために内包物３２１を明示しているが、内包物３２１は実際の弾性画像上には直接現れない。伝播解析部１５は、生成した弾性画像を制御部１１に出力し、制御部１１は弾性画像を弾性画像記憶部２１に出力する。さらに、制御部１１は、弾性画像に対して画面表示用の画像データとなるよう幾何変換を行い、幾何変換後の弾性画像を表示部３に出力する。

次に、せん断波の伝播を評価する（ステップＳ８０）。評価部１６は、制御部１１を介して波面画像記憶部２０から波面画像を取得し、波面の位置、振幅、進行方向および速度から、せん断波がどのように伝播しているのかを評価する。評価部１６は、例えば、せん断波が着目領域の全域に伝播しているか否か、せん断波の伝播速度や向きに一様性のない部分があるか否かを評価する。評価の方法としては、例えば、せん断波の伝播の有無を示した伝播範囲図、せん断波の速度分布を示した速度分布図などを生成することで行う。以下、図６に示すこれらの具体例を用いて説明する。

図６（ａ）は、せん断波の伝播の有無を示した伝播範囲図であり、所定の振幅以上のせん断波が通過した領域を、せん断波伝播領域として示している。伝播範囲図３３０においては、着目領域３３１（ｘ_L≦ｘ≦ｘ_R、ｙ_S≦ｙ≦ｙ_D）の全域が、せん断波伝播領域３３２となっている。一方、伝播範囲図３４０においては、着目領域３４１（ｘ_L≦ｘ≦ｘ_R、ｙ_S≦ｙ≦ｙ_D）のうち、振幅が十分なせん断波が通過しているのはせん断波伝播領域３４２（ｙ_S≦ｙ≦ｙ_p）だけで、非伝播領域３４３（ｙ_p≦ｙ≦ｙ_D）は、せん断波が通過していないか、通過したせん断波の振幅が所定の振幅未満である領域である。また、図６（ｂ）は、せん断波の速度分布を示した速度分布図であり、せん断波の伝播速度および向きを色で示したものである。速度分布図３５０では、せん断波の伝播向きがｘ軸の負の方向（左向き）で速度が所定範囲の領域３５２を緑、せん断波の伝播向きがｘ軸の正の方向（右向き）の領域のうち、速度が所定範囲の領域３５３を青、速度が所定範囲より速い領域３５４を赤で示している。一方、速度分布図３６０では、せん断波の伝播向きが左向きで速度が所定範囲の領域３６２を緑、せん断波の伝播向きが右向きの領域のうち、速度が所定範囲の領域３６３を青、速度が所定範囲より速い領域３６４を赤で示しており、さらに、せん断波の速度が取得できなかった領域３６５を黒で示している。なお、分類はこの場合に限られず、例えば、伝播速度をスカラー的に評価して向きに関係なく大きさ（絶対値）のみに依存して色分けを行ってもよいし、速度に対する段階分けをより細かく行ってもよい。

次に、プッシュパルスの送信プロファイルを調整して再測定を行う必要があるか否かを判断する（ステップＳ９０）。具体的には、制御部１１が、評価部１６の評価結果を基に、再測定の要否を判断する。例えば、伝播範囲図中に非伝播領域が存在する、速度分布図中にせん断波の速度が取得できなかった領域が存在する、のいずれかに該当する場合、制御部１１は、プッシュパルスの送信プロファイルを調整して再測定を行う必要があると判断する。

制御部１１は、プッシュパルスの送信プロファイルを調整して再測定を行う必要がないと判断した場合、そのまま処理を終了し、弾性画像を表示したまま、ユーザからの入力を受け付ける。
一方、制御部１１がプッシュパルスの送信プロファイルを調整して再測定を行う必要があると判断した場合、送信プロファイルの調整を行う（ステップＳ１００）。プッシュパルス調整部１７は、ステップＳ３０で用いたプッシュパルスの送信プロファイルと、ステップＳ９０で評価部１６が評価した内容とを、制御部１１から取得する。そして、評価部１６が評価した内容に基づき、プッシュパルスの送信プロファイルを最適化する。以下、フローチャート及び具体例を用いて送信プロファイルの最適化について説明する。図７は、プッシュパルスの送信プロファイル動作を示すフローチャートである。

まず、プッシュパルス調整部１７は、着目領域内にせん断波が伝播していない非伝播領域が存在するか否かを判定する（ステップＳ１１０）。着目領域内に非伝播領域が存在する場合には、次に、非伝播領域がどのような形状であるかを判断する。具体的には、まず、非伝播領域が素子列方向に伸びている（ｘ座標にほぼ関わりなくｙ座標が特定の範囲の領域）か、深さ方向に伸びている（ｙ座標にほぼ関わりなくｘ座標が特定の範囲の領域）かを判定する（ステップＳ１１１）。ここで、非伝播領域がせん断波伝播領域に囲まれている場合（ｘ座標が特定の範囲、かつ、ｙ座標が特定の範囲の領域）は、深さ方向に伸びているとして判定を行う。理由は後述する。

非伝播領域が素子列方向に伸びている場合、深さ方向における非伝播領域の位置を判定する（ステップＳ１１２）。非伝播領域が浅い側の端、または、深い側の端である場合（着目領域における一定の深さより浅い領域全て、または、一定の深さより深い領域全てが非伝播領域である場合）、プッシュパルスの焦点位置の深さ（ｙ座標）を、非伝播領域に近づけるように調整する(ステップＳ１１３)。より詳しくは、深さが一定の深さ以上の領域が非伝播領域である場合は、プッシュパルスの焦点位置をより深い位置にずらし、深さが一定の深さ以下の領域が非伝播領域である場合は、プッシュパルスの焦点位置をより浅い位置にずらす。これは、図８（ａ）の伝播範囲図４００に示すように、プッシュパルスの焦点位置４０２が着目領域４０１に対して浅すぎるため、着目領域４０１内に非伝播領域４０４が発生しているためである。または、図８（ｂ）の伝播範囲図４１０に示すように、プッシュパルスの焦点位置４１２が着目領域４１１に対して深すぎるため、着目領域４１１内に非伝播領域４１４が発生しているためである。なお、伝播範囲図４００および伝播範囲図４１０では、説明のために波面を表示している。図８（ａ）の伝播範囲図４００に対しては、せん断波伝播領域４０３と非伝播領域４０４との境界線４０５におけるｙ座標の平均値ｙ_pを用いて、例えば、新たなプッシュパルスの焦点位置４０６を、もとのプッシュパルスの焦点位置４０２から、ｙ軸の正の方向（深さが増す方向）に（ｙ_D−ｙ_p）／２だけ移動させる。同様に、図８（ｂ）の伝播範囲図４１０に対しては、せん断波伝播領域４１３と非伝播領域４１４との境界線４１５におけるｙ座標の平均値ｙ_pを用いて、例えば、新たなプッシュパルスの焦点位置４１６を、もとのプッシュパルスの焦点位置４１２から、ｙ軸の負の方向（深さが減る方向）に（ｙ_p−ｙ_S）／２だけ移動させる。すなわち、非伝播領域の深さ方向（ｙ方向）の平均幅の１／２だけ、プッシュパルスの焦点位置を、非伝播領域が存在する側に移動させる。送信プロファイルとしては、「深さ」のパラメータを変更することとなる。なお、せん断波伝播領域と非伝播領域との境界のｙ座標ｙ_pとして、非伝播領域がせん断波伝播領域より深い場合はせん断波伝播領域と非伝播領域との境界線のｙ座標の最大値、非伝播領域がせん断波伝播領域より浅い場合はせん断波伝播領域と非伝播領域との境界線のｙ座標の最小値を用いてもよい。また、プッシュパルスの焦点位置のｙ座標の移動距離は、例えば、（ｙ_D−ｙ_p）／４でもよいし、（ｙ_D−ｙ_p）でもよいが、（ｙ_D−ｙ_p）／２が、再度のプッシュパルスの焦点位置の深さ調整を行う必要がなくなる可能性が高いため、より好ましい。

一方、非伝播領域が素子列方向に延びており、かつ、非伝播領域がせん断波伝播領域に挟まれている場合（非伝播領域より浅い側と深い側の双方にせん断波伝播領域が存在する場合）、プッシュパルスの焦点位置を素子の並び方向（ｘ方向）にずらす。これは、図８(ｃ)の伝播範囲図４２０に示すように、プッシュパルスの焦点位置４２２が硬い組織４２５内に存在するため、プッシュパルスによる組織の移動が十分に発生せず、着目領域４２１においてせん断波伝播領域４２３に挟まれるように非伝播領域４２４が生じていることが想定されるためである。この場合、プッシュパルスの焦点位置を素子の並び方向（ｘ方向）にずらすことで、新たな焦点を硬い組織の外部に設定する。例えば、伝播範囲図４３０に示すように、着目領域４２１の中心を通り深さ方向（ｙ軸方向）と平行な直線４３１を仮定し、直線４３１を基準にプッシュパルスの焦点位置４２２と線対称となる位置を新たなプッシュパルスの焦点位置４３２とする。具体的には、着目領域４２１のｘ座標の範囲がｘ_L≦ｘ≦ｘ_Rであり、プッシュパルスの焦点位置４２２のｘ座標がｘ_fであるときに、調整後のプッシュパルスの焦点位置４３２のｘ座標をｘ_nf＝ｘ_L＋ｘ_R−ｘ_fとする。送信プロファイルとしては、「深さ」並びに「プッシュパルスの伝播方向と超音波探触子とがなす角度」を変更することで、焦点位置のｘ座標が変更させる。これにより、プッシュパルスの焦点位置が硬い組織４２５の外部に設定され、せん断波伝播領域４３３が着目領域４２１全域となることが期待できる。なお、新たなプッシュパルスの焦点位置のｘ座標は上述のｘ_nfに限られず、例えば、ｘ_fから所定の量（具体的には、着目領域４２１のｘ方向の幅の１／４など）だけ着目領域４２１の中心方向にずらした位置、としてもよい。

一方で、非伝播領域が深さ方向に伸びている場合（非伝播領域がせん断波伝播領域に囲まれている場合を含む）、素子列方向における非伝播領域の位置を判定する（ステップＳ１１５）。なお、非伝播領域がせん断波伝播領域に囲まれている場合、非伝播領域がせん断波伝播領域に挟まれているとして判定する。
非伝播領域が深さ方向に伸びている場合において、非伝播領域が素子列方向の一方の端である場合（着目領域における一定のｘ座標よりｘが増加する側全て、または、一定のｘ座標よりｘが減少する側全てが非伝播領域である場合）、プッシュパルスを強めるように調整する（ステップＳ１１６）。これは、図８（ｅ）の伝播範囲図４４０に示すように、プッシュパルスの焦点位置４４２から伝播したせん断波がせん断波伝播領域４４３の端で減衰してしまい、せん断波が着目領域４４１全域に行き渡らず、非伝播領域４４４が生じているためである。具体的には、以下のような調整のいずれかを行う。
（１）プッシュパルスの送出に用いる振動子数、プッシュパルスの送信時間、振動子に印可するＡＲＦＩ駆動信号の電圧値の少なくとも１つを増加させ、プッシュパルスによって生じるせん断波の振幅を増加させる。具体的には、例えば、プッシュパルスの送出に用いる振動子数を４つ増加させる、プッシュパルスの送信時間を２０μｓだけ増加させる、ＡＲＦＩ駆動信号の電圧値を１０％上昇させる、などである。なお、これらの２つ以上を同時に行ってもよい。
（２）連続して送出するプッシュパルスの数を増加させ、非伝播領域内または近傍に新たなプッシュパルスの焦点位置を設けることで、非伝播領域にせん断波を伝播させる。より具体的には、例えば、図８（ｆ）の伝播範囲図４５０に示すように、着目領域４４１の中心を通り深さ方向（ｙ軸方向）と平行な直線４５１を仮定し、直線４５１を基準にプッシュパルスの焦点位置４４２と線対称となる位置を２つ目のプッシュパルスの焦点位置４５２とする。

なお、調整内容としては、プッシュパルスの連続送出回数を増加させることが最も好ましく、プッシュパルスの送出に用いる素子数の増加、プッシュパルスの送信時間の延長、ＡＲＦＩ駆動信号の電圧値の上昇では、前者が後者より好ましい。これは、以下の理由による。被検体に与える影響を考えると、組織内の変位は小さいことが好ましく、特に、組織内の変位の最大値を小さくすることが好ましい。一方で、硬さの測定精度の向上のためには変位の平均値や最小値が大きいことが好ましい。したがって、この２つを両立するためには、変位の最大値と変位の平均値との差が小さい方が好ましく、強いプッシュパルスを少ない回数送信するよりは、強すぎないプッシュパルスを数多く送信する方がより好ましいからである。また、強いプッシュパルスを送信する場合は、超音波の音圧を上昇させるより超音波の送信時間を長くする方が、さらに、超音波の送信時間を長くするより超音波の通過面積と集束度を上昇させる方が、焦点以外のプッシュパルスの通過範囲にもたらす影響を小さくできるためより好ましい。

一方で、非伝播領域が深さ方向に伸びている場合において、非伝播領域がせん断波伝播領域に挟まれている、または囲まれている場合、焦点位置の異なる２つ目のプッシュパルスを送出するように調整する（ステップＳ１１７）。このとき、プッシュパルス調整部１７は、非伝播領域を素子列方向に挟むように、せん断波伝播領域内に２つのプッシュパルスの焦点位置を設ける。この場合、非伝播領域の内部にせん断波の反射源（例えば、硬い組織）がある可能性が高い。そのため、（ｉ）図８（ｇ）の伝播分布図４６０に示すように、プッシュパルスの焦点４６２から伝播するせん断波の進行波と、着目領域４６１内に存在する反射源４６５から拡がるせん断波の反射波との干渉により、非伝播領域４６４内でせん断波が正常に検出できない、（ｉｉ）反射源が周囲の組織に比べて非常に硬い場合などにおいて、反射源がせん断波の多くを反射してしまい、反射源の内部にせん断波がほとんど伝播しない（この場合、反射源が存在する領域が非伝播領域となる）、のいずれかが考えられる。いずれの場合も、非伝播領域にせん断波を伝播させるためには、非伝播領域の近傍の複数個所をプッシュパルスの焦点位置とすることが好ましく、非伝播領域を素子列方向に挟むように、せん断波伝播領域内に２つのプッシュパルスの焦点位置を設ければよい。例えば、図８（ｇ）の伝播分布図４６０に対しては、せん断波伝播領域４６３内とせん断波伝播領域４６６内とにそれぞれ１つずつプッシュパルスの焦点位置を設ける。例えば、新たなプッシュパルスの焦点位置４６７を断波伝播領域４６６内に追加する。このとき、せん断波伝播領域４６３内のプッシュパルスの焦点位置は焦点位置４６２をそのまま用いてもよい。

また、着目領域内にせん断波が伝播していない領域が存在しない場合、せん断波の伝播速度が局所的に速い領域の近傍、および／または、せん断波の速度が一様である領域の内部に、プッシュパルスの焦点位置を設定する（ステップＳ１１８）。このとき、プッシュパルスの送信回数を１増加させ、ステップＳ３０におけるプッシュパルスの焦点位置に加えて、新たに設定した焦点位置にもプッシュパルスを送信するようにする。例えば、図８（ｈ）の速度分布図４７０では、着目領域４７１に対し、せん断波がプッシュパルスの焦点４７６からｘ方向（左右）に伝播し、領域４７２と領域４７３ではせん断波の伝播速度が一様であるものの、領域４７４ではせん断波の伝播速度が局所的に速く、領域４７５ではせん断波の伝播速度が取得できていない。このとき、既存のプッシュパルスの焦点位置４７６に加え、せん断波の伝播速度が速い領域４７４の近傍であって、せん断波の伝播速度が一様な領域４６３の内部にプッシュパルスの焦点４７７を追加する。このとき、新たなプッシュパルスの焦点４７７は、次の２条件を満たしていることが好ましい。１つは、せん断波の伝播速度が局所的に変化している領域の近傍にプッシュパルスの焦点が存在していることである。これは、せん断波の伝播速度が局所的に変化している領域に対応する、硬さが周囲と異なる組織に対し、せん断波の振動源であるプッシュパルスの焦点を近づけることで、当該組織に振幅の大きいせん断波を伝播させ、当該組織の硬さを精度よく取得するためである。もう１つは、せん断波の伝播速度が一様である領域の内部にプッシュパルスの焦点が存在していることである。これは、せん断波の伝播速度が局所的に変化している領域以外を原因とするせん断波の屈折や反射をできるだけ防ぎ、せん断波の伝播速度の信頼性を高めるためである。なお、上記２条件は、いずれか一方のみを満たす、としてもよい。また、プッシュパルスの焦点位置を追加した結果、２以上のプッシュパルスの焦点位置が互いに近接することとなった場合は、新たに設定したプッシュパルスの焦点位置の近傍に存在する以前のプッシュパルスの焦点位置にはプッシュパルスを送信しない、としてもよい。

図２に戻って説明を続ける。制御部１１は、調整後のプッシュパルスの送信プロファイルをプッシュパルス調整部１７から受け取ると、調整後のプッシュパルスの送信プロファイルを用いてプッシュパルスの送信と弾性画像の生成を行うため、ステップＳ３０からステップＳ９０を再度実行する。
＜まとめ＞
上記構成により、せん断波の伝播状況に基づいてプッシュパルスの送信プロファイルを最適化することが可能となる。そのため、検査者が弾性画像に基づいて試行錯誤することなく、着目領域における被検体の組織の硬さを高精度に測定することができる。また、上記構成により、プッシュパルスの送信プロファイルを最適化した後そのまま組織の硬さの再測定が行われる。そのため、検査者は、弾性画像が確定するまで超音波探触子を被検体に当てておくだけでよい。

≪変形例１≫
実施の形態では、せん断波が着目領域全体に伝播していない場合について、非伝播領域の位置及び形状に基づいてプッシュパルスの送信プロファイルを調整する場合について説明した。
しかしながら、せん断波が着目領域全体に伝播している場合においても、例えば、せん断波の伝播速度の分布も考慮した上で、プッシュパルスの送信プロファイルを調整してもよい。

＜送信プロファイルの調整＞
図９に具体例を示す。図９（ａ）は、図８（ｅ）と同様に、伝播範囲図４８０において、非伝播領域４８３が着目領域４８１の素子列方向の一方の端で深さ方向に伸びている場合を示している。なお、プッシュパルスの焦点位置は４８４であり、領域４８２はせん断波の伝播領域である。一方、図９（ｂ）は、図９（ａ）の伝播範囲図４８０に対応する、せん断波の速度分布図４９０である。着目領域４８１において、非伝播領域４８３に対応する領域４９３ではせん断波の速度が取得できていない。一方、せん断波の伝播領域４８３に対応する領域には、せん断波の速度が一様である領域４９１と、せん断波の速度が局所的に高い領域４９２とが存在する。このような場合、新たに追加するプッシュパルスの焦点位置としては、着目領域４８１の中心を通り深さ方向（ｙ軸方向）と平行な直線を基準にプッシュパルスの焦点位置４８４と線対称となる位置４８５よりも、ａ）せん断波の速度が一様である領域４９１の内部であって、ｂ）せん断波の速度が局所的に高い領域４９２の近傍である、の２条件を同時に満たす位置４９４の方が好ましい。なぜならば、せん断波の速度が局所的に高い領域４９２の存在から、領域４９２を含む硬い組織４９５の存在が推測できる一方で、非伝播領域４８３（せん断波の速度が取得できていない領域４９３）における、組織４９５が占める領域が明確ではない。そのため、非伝播領域４８３内に新たな焦点位置４８５を設けた場合、焦点位置４８５が組織４９５の内部に存在する危険性がある。したがって、このような場合は、非伝播領域４８３に近接したせん断波の伝播領域４８２内に新たなプッシュパルスの焦点位置４９４を設けることが好ましく、さらに、焦点位置４９４は、せん断波の速度が一様である領域４９１の内部であることが好ましく、また、せん断波の速度が局所的に高い領域４９２の近傍であることが好ましい。

なお、この場合に限らず、プッシュパルスの焦点位置を変更、または、追加する場合において、せん断波の伝播速度の分布を考慮し、せん断波の伝播速度が局所的に変化している領域の近傍に、および／または、せん断波の伝播速度が一様である領域の内部に、プッシュパルスの焦点が存在するように制御してよい。
＜まとめ＞
このようにすることで、プッシュパルスの焦点位置を変更、または、追加する場合において、せん断波が着目領域全体に伝播させるという観点のみならず、硬さの異なる領域における測定精度の向上、および／または、せん断波の伝播速度の信頼性向上を図ることで、再測定後にさらに再測定の必要が生じる可能性を減少させ、組織の硬さの再測定の回数を減らし、測定時間の短縮および被検体への負担を軽減することができる。

≪変形例２≫
実施の形態および変形例１では、せん断波が着目領域全域に伝播したか否か、および、せん断波の伝播速度の分布を、せん断波の波面を用いて判断した。
本変形例では、せん断波の伝播有無、および／または、速度分布をより簡易な方法で判断することに特徴がある。

＜せん断波の伝播有無の解析＞
本変形例では、伝播範囲図の作成にあたって、着目領域内の各位置をせん断波が通過したか否かを波面画像から生成せず、代わりに、変位画像から生成する。すなわち、いずれかの時刻において変位が検出された場所は、波面が通過したものとみなす。より具体的には、着目領域内の各位置に対して、複数の変位画像における同一位置の変位の絶対値（大きさ）の最大値を検出し、最大値が所定の閾値以上であれば当該位置をせん断波が通過したものとし、最大値が所定の閾値未満であれば当該位置をせん断波が通過しなかったものとする。

したがって、いずれの時刻においても変位が検出できなかった領域があれば、これを非伝播領域であるものとし、この非伝播領域にせん断波が伝播するようにプッシュパルスの送信プロファイルを調整する。
＜せん断波の伝播速度の解析＞
本変形例では、せん断波の速度の解析にあたって、着目領域内の各位置をせん断波が通過した速度を示す速度分布図を作成せず、代わりに、断層画像をそのまま用いる。すなわち、断層画像において、輝度が局所的に低い領域を、せん断波の速度が局所的に異なる領域とみなす。これは、腫瘍等の硬さの異なる領域は周囲に比べて輝度が低くなることが多く、当該領域が硬さの異なる組織（腫瘍等）である可能性が高いからである。

したがって、プッシュパルスの焦点位置を変更、または、追加する場合には、断層画像において輝度が局所的に低い領域ではない箇所に設定する、および／または、断層画像において輝度が局所的に低い領域の近傍に設定する。
＜まとめ＞
このようにすることで、波面を検出せずに評価を行うことができる。そのため、波面画像記憶部が必要ない。さらに、弾性画像を生成する際、波面画像に基づく方法ではなく、変位画像に直接基づく方法を用いた場合、波面の検出を行う必要がなく、伝播解析部がなくても本発明の実施をすることができる。なお、変位画像に直接基づいて弾性画像を生成する方法としては、例えば、変位画像内の各画素において変位が最初に検出された時刻をせん断波の通過時刻とみなす従来手法を用いることができる。

また、断層画像の輝度を用いて腫瘍領域の有無を予測することで、波面速度や弾性が不明でもプッシュパルスの送信プロファイルを調整することができる。そのため、波面速度の算出や弾性測定を行わずにプッシュパルスの送信プロファイルを調整したり、断層画像に基づいて最初のプッシュパルスの送信プロファイルを調整したりすることが可能となる。

または、本変形例による簡易評価に直接基づいて送信プロファイルの調整を行うのではなく、簡易評価により再測定が必要であると判定された場合に、波面を検出して波面に基づいた評価を行ったうえでプッシュパルスの送信プロファイルを変更するものとしてもよい。
≪変形例３≫
実施の形態、および変形例１、２では、ステップＳ６０におけるせん断波の伝播解析を、変位領域の抽出、細線化処理、空間フィルタリング、時間フィルタリングの手順によって行った。

本変形例では、せん断波の伝播解析を、各場所の変位のピーク時刻検出、時間フィルタリング、空間フィルタリングの手順で行うことに特徴がある。
ステップＳ６０におけるせん断波の伝播解析以外は実施の形態と同様であるため、以下、差分である、せん断波の伝播解析についてのみ説明する。
＜せん断波の伝播解析＞
最初に、各画素について時間領域毎に変位を相関処理し、変位のピークを検出する。まず、伝播解析部は変位画像を変位量記憶部１９から取得する。以下、図１０の模式図を用いて説明する。図１０（ａ）は、変位画像の一例を示している。図３、図４と同じく、図中の“○”は着目領域における被検体内の組織の一部を示しており、プッシュパルスを印加する前の位置は破線の交点である。また、ｘ軸は超音波探触子２における素子の並ぶ方向、ｙ軸は、被検体の深さ方向である。伝播解析部は、座標ごとに変位量δが極値となる時刻を検出する。ここで、極値とは、極大値、および、極小値（変位量δが負の値で、絶対値が極大）のことである。本変形例においては、変位量δが極大値（変位量δが正の値）となる時刻のみを検出することとする。図１０（ａ）のグラフ５１０は、座標（ｘ₁，ｙ₁）の点５０１における変位量δを時間ｔの関数として示したものである。なお、ここでは、変位量δを１次元ベクトルとして扱うため、変位量δはｙ軸方向の変位のみを示す値であり、かつ、深い方向（ｙ座標が大きくなる方向）への変位は正の値、浅い方向（ｙ座標が小さくなる方向）への変位は負の値としている。伝播解析部は、例えば、断層画像信号が取得された時間を複数の区間５２１、５２２、５２３に分ける。そして、伝播解析部は、区間５２１における変位量δ、区間５２２における変位量δ、区間５２３における変位量δについて相関処理を行うことで、変位量δがピークとなる時刻ｔ₁、ｔ₂を検出する。ここで、断層画像信号が取得された時間を複数の区間に区切っているのは、着目点を反射波が通過するなど、波面が着目点を複数回通過することがあり得るからである。なお、本変形例では、変位量δが極大となる時刻のみを検出するとしたが、さらに、相関処理において相関値が負の値で絶対値が極大となる時刻、すなわち、変位量δが負の値で、その絶対値が極大となる時刻を検出してもよい。そのようにすることで、位相が反転している波面をさらに検出することができる。

次に、伝播解析部は、変位が大きい画素の集合を波面として抽出する。具体的には、伝播解析部は、各時刻について、変位量δがピークとなっている点の集合を変位領域として取得し、変位領域のそれぞれを波面として抽出する。すなわち、図１０（ｂ）の波面画像５３０に示すように、時刻ｔ₁において、座標（ｘ₁，ｙ₁）の点を含む変位領域５３１や、変位領域５３２、５３３、５３４、５３５、５３６が、それぞれ波面として抽出される。なお、上述したように、変位量δが負の値で絶対値が極大値となっている画素の集合を基に、上記検出した波面とは位相が反転している位相反転波面をさらに抽出してもよい。伝播解析部は、抽出した波面を時刻ごとの波面画像データとして保持する。

次に、伝播解析部は、複数の波面画像データに対して時間フィルタリングを行い、伝播していない波面を除去する。具体的には、ステップＳ６０の時間フィルタリングと同様、時間的に連続する２以上の波面画像において、各波面位置の時間変化を検出し、速度が異常である波面をノイズとして除去する。伝播解析部は、例えば、時刻ｔ＝ｔ₁の波面画像５４０と時刻ｔ＝ｔ₁＋Δｔの波面画像５５０との間、時刻ｔ＝ｔ₁＋Δｔの波面画像５５０と時刻ｔ＝ｔ₁＋２Δｔの波面画像５６０との間のそれぞれにおいて、波面位置の時間変化の検出を行う。波面位置の時間変化の検出はステップＳ６０の時間フィルタリングと同様、例えば、波面画像５４０内の各波面に対し、波面画像５５０において、当該波面と同じ位置を中心に、波面と垂直な向きにΔｔの間にせん断波が移動しうる領域で相関処理することで行う。例えば、波面画像５４０と波面画像５５０との間で、波面５４１が波面５５１の位置に、波面５４２が波面５５２の位置に、波面５４３が波面５５３の位置に、それぞれ移動したと検出する。ここで、波面５４４の移動先が検出されていないのは、以下の理由による。波面Ａと波面Ｂとを含む波面画像Ｘと、波面Ｃを含む波面画像Ｙとの相関処理において、波面画像Ｙで、波面Ａを中心とする領域と波面Ｂを中心とする領域とが重なっており、その重複領域に波面Ｃが存在する場合、以下のように移動先を判定する。すなわち、波面Ａと波面Ｃとの相関値と、波面Ｂと波面Ｃとの相関値とをそれぞれ算出し、（１）いずれか一方の相関値のみが閾値を超えていれば、波面Ａと波面Ｂのうち、波面Ｃとの相関値が閾値を超えている側を波面Ｃの移動元とする、（２）両方の相関値が閾値を超えていれば、波面Ａと波面Ｂのうち、波面Ｃとの相関値が高い方を波面Ｃの移動元とする、（３）いずれの相関値も閾値を下回っていれば、波面Ｃの移動元は存在しない、と判定する。波面５４２、５４３、５５２、５５３は略同じ長さであるのに対して波面５４４は他の波面より極端に短いため、波面５４４と波面５５３の相関値より波面５４３と波面５５３の相関値の方が高い。そのため、波面画像５５０において波面５４４と同じ位置を中心とした領域内に、波面５４４の移動先となる波面が存在しない。移動先が存在しない波面５４４は、ノイズとして消去する。同様に、波面画像５５０と波面画像５６０との間で、波面５５１が波面５６１の位置に、波面５５２が波面５６２の位置に、波面５５３が波面５６３の位置に、それぞれ移動したと検出する。伝播解析部は、時間フィルタリング後の波面画像データを保持する。なお、伝播解析部は、実施の形態で述べたように、さらに、波面の対応情報を保持してもよい。

次に、伝播解析部は、時間フィルタリング後の波面画像データに対して空間フィルタリングを行い、長さが短い波面を除去する。具体的には、ステップＳ６０における空間フィルタリングと同様、各波面画像において、波面の長さを検出し、所定の長さより短い波面、例えば、全ての波面の長さの平均値の１／２よりも長さが短い波面を、ノイズとして削除する。具体的には、図１０（ｄ）の波面画像５７０に示すように、波面５７１〜５７３の長さの平均値を算出し、それよりも短い波面５７１を、ノイズとして消去する。これにより、誤検出された波面を消去できる。なお、伝播解析部は、波面の対応情報から、消去した波面と同一である他の波面画像内の波面をまとめて消去してもよい。図１０（ｅ）に示すように、空間フィルタリング後の波面画像４８０には、波面４７２、４７３のみが残る。

伝播解析部は、空間フィルタリング後の複数の波面画像データと、波面の対応情報とを制御部を介して、波面画像記憶部２０に出力する。
＜まとめ＞
本変形例では、組織内の各位置における変位の時間変化に基づいて、伝播解析を行う場合について説明した。各座標について時系列に変位のピークを検出するため、時間変化のない変位を検出せず、基準断層画像信号内のノイズや、基準断層画像信号と断層画像信号との位置ずれなどによる、時間変化しないノイズを変位として誤検出することを防ぐことができる。また、各座標について時系列に変位のピークを複数回検出するため、１つの場所を波面が複数回通過する場合、特に、透過波と反射波が通過するような場合にも、波面を抽出することができる。また、各座標について時系列に変位のピークを検出する際、相関値が負の値で絶対値が極大となる時刻を検出することで、せん断波の位相が反転した場合に、その位相が反転した波面を検出することができ、固定端反射となるようなせん断波の反射に対しても、反射の波面を検出することができる。

≪変形例４≫
実施の形態、および、変形例１〜３では、ステップＳ７０において、生成した弾性画像を表示部３に出力するものとした。
本変形例では、プッシュパルスの送信プロファイルが確定するまで、ユーザに通知画面を出力することで、超音波診断装置による被検体の硬さ測定が完了しているか、未完了であるかを通知することに特徴がある。

＜動作＞
本変形例に係るフローチャートを図１１に示す。なお、図１１において、図２と同様の動作については同じステップ番号を用い、説明を省略する。
ステップＳ２１０において、図１２（ａ）に示すように、測定中を示す表示６０１を表示部３の画面６００に出力する。

ステップＳ２７０では、ステップＳ７０と同様に弾性画像を生成するが、表示部３には出力しない。
さらに、ステップＳ１００で再測定を行わないとした場合において、制御部は、測定中を示す表示に替えて、最近のステップＳ７０で生成した弾性画像を表示部３に出力する（ステップＳ２２０）。

＜まとめ＞
このようにすることで、弾性画像を生成したプッシュパルスが再調整不要となるまで、表示部３に弾性画像が表示されない。したがって、ユーザは、弾性画像が出力されるまでは、超音波探触子を被検体に押し当てたまま待っているだけでよく、弾性画像の表示によって測定終了を知ることができる。

なお、ステップＳ２１０において表示部３に何も出力せず、ステップＳ２７０においても表示部３に何も出力しない、としてもよい。この場合、弾性画像の表示があることでユーザは測定終了を知ることができるので、ユーザは、弾性画像が出力されるまでは、超音波探触子を被検体に押し当てたまま待っているだけでよい。
また、ステップＳ２１０において、図１２（ｂ）に示すようなアイコン６１１を表示部３に出力する画面６１０に重畳し、ステップＳ２７０では、ステップＳ７０と同様に生成した弾性画像６１２を出力し、ステップＳ２２０において、アイコン６１１の重畳表示を終了するとしてもよい。このようにすれば、アイコンが画面から消えることによりユーザは測定終了を知ることができ、さらに、弾性画像の取得状況を確認することができる。

また、ステップＳ２７０を、ステップＳ６０とＳ８０の間ではなく、ステップＳ１００で再測定を行わないとした場合において、ステップＳ２２０の前に行ってもよい。このようにすることで、表示しない弾性画像を生成する必要がない。この場合、弾性画像記憶部がなくてもよい。
≪実施の形態に係るその他の変形例≫
（１）実施の形態および各変形例では、プッシュパルスの焦点位置を超音波診断装置１が送信プロファイルの初期値に基づいて決定し、せん断波の伝播状況の評価に応じて変更または追加するとした。しかしながら、プッシュパルスの焦点位置をユーザが調整してもよい。例えば、プッシュパルスの焦点位置を表示部３に表示し、ユーザが超音波探触子の位置を変更する、または、超音波診断装置１がプッシュパルスの焦点位置をユーザから受け取る、としてもよい。このとき、超音波診断装置１が受け取るプッシュパルスの焦点位置は、プッシュパルスの焦点位置そのものでもよいし、例えば、焦点が存在しうる範囲であってもよい。このようにすることで、ユーザが断層画像や弾性画像等に基づいて焦点位置を選択できる場合、ユーザの所望する条件で断層画像を生成することができる。

（２）実施の形態および各変形例では、２以上のプッシュパルスを送信する送信プロファイルに対して、２以上のプッシュパルスを連続的に送信する場合について説明した。しかしながら、プッシュパルスの送信と、超音波の複数回送受信との組み合わせをプッシュパルスの数だけ繰り返すとしてもよい。このようにすることで、最初のプッシュパルスを振動源とするせん断波の伝播解析において後のプッシュパルスを振動源とするせん断波の影響を防ぐことができ、より高精度な硬さの測定が可能となる。

（３）実施の形態および各変形例では、ステップＳ９０でプッシュパルスの送信プロファイルを調整して再測定を行うとした場合に、ステップＳ１００で調整した送信プロファイルを用いてステップＳ３０〜Ｓ８０を行い、再度、ステップＳ９０でプッシュパルスの送信プロファイルを調整して再測定を行うか否を判断する場合について説明した。しかしながら、例えば、２度目のステップＳ３０〜Ｓ８０の後については、ステップＳ９０を行わず、そのまま測定を終了するとしてもよい。このようにすることで、例えば、超音波探触子２と着目領域との位置関係上、送信プロファイルをどう調整しても伝播解析で好ましくない結果が出るなどの場合、プッシュパルスの送信プロファイルの調整と再測定とを過度に繰り返すことを防ぐことができる。なお、プッシュパルスの再送信と再測定との上限回数は１回に限らず、被検体への影響や測定時間等を考慮した回数であればよく、例えば、３回であってもよい。この場合においては、例えば、制御部１１は、評価部１６による評価が最も好ましかった測定の断層画像を表示するとしてもよい。

または、評価部１６は、評価結果を保持しておき、２度目以降のステップＳ９０において、評価の結果が以前のステップＳ９０における評価の結果より好ましくない場合、評価の結果が最も好ましい送信プロファイルに係る弾性画像を表示して測定を終了するとしてもよいし、あるいは、最後の評価結果を考慮して同一結果とならないようにしながら、保持している評価結果のうち結果が最も好ましいものに基づいて再度送信プロファイルを調整するとしてもよい。このとき、最後の評価結果が所定の基準に達していれば再度送信プロファイルの調整を行い、最後の評価結果が所定の基準に達していなければ評価の結果が最も好ましい送信プロファイルに係る弾性画像を表示して測定を終了するとしてもよい。

（４）実施の形態および各変形例では、波面画像を基に生成した伝播範囲図、または、変位画像を基に簡易作成した伝播範囲図に基づいて、せん断波の伝播していない領域に着目してプッシュパルスの送信プロファイルを調整する場合について説明した。しかしながら、せん断波の伝播していない領域を検出せず、例えば、せん断波の伝播の一様性のみに基づいてプッシュパルスの送信プロファイルを調整するとしてもよい。例えば、ａ）せん断波の伝播が一様な領域内にプッシュパルスの焦点がある、ｂ）せん断波の伝播が局所的に異なる領域の近傍にプッシュパルスの焦点がある、のいずれか一方を、または、両方を同時に、満たすように送信プロファイルを調整する。さらに、せん断波の一様性を評価する際に、実施の形態で説明したように、せん断波の速度と向きを用いて評価を行ってもよいし、変形例２で説明したように断層画像を用いた簡易評価を行ってもよい。あるいは、せん断波の伝播の一様性を評価する際に、弾性画像を用いてもよい。なぜならば、弾性値からせん断波の速度の絶対値（大きさ）を評価することができるからである。したがって、弾性値の一様な領域はせん断波速度の一様な領域、弾性値の局所的に異なる領域はせん断波速度の局所的に異なる領域、と評価することができる。

（５）変形例３では、伝播解析部は、相関処理により、変位画像から変位量がピークとなる時刻を検出する場合について説明したが、本発明は必ずしもこの場合に限られない。例えば、相関処理による相関値に対してしきい値を設け、しきい値を上回る相関値が検出された時間帯を検出するとしてもよい。このようにすることで、変位量がピークとなる時刻をピンポイントで検出するだけではなく、その周辺時刻を検出して波面を正確に検出することが可能となる。また、逆に、相関処理による相関値に対して下限しきい値を設け、下限しきい値を下回る相関値が検出された場合、ピークとして検出しないとしてもよい。このようにすることで、波面が通過していない時刻における揺動などを波面と誤検出することを防ぐことができる。

（６）実施の形態および各変形例では、ステップＳ２０およびＳ４０において、超音波信号取得部１３が、送信超音波が送信フォーカス点に同時に届くように送信ビームフォーミングを行い、送信イベントごとに素子受信信号に整相加算を行い、音響線信号を生成する場合について説明したが、本発明は必ずしもこの場合に限られない。例えば、超音波信号取得部１３は、超音波探触子２を構成する各振動素子から同時に超音波を送出することで、平面波状の送信超音波を送出してもよい。この場合、受信フォーカス点を送信フォーカス点の近傍や同一線上に限定する必要がないため、１つの断層画像を構成する音響線信号を、１送信イベントごとに１つではなく１画像分まとめて取得することができる。そのため、断層画像信号のフレームレートを向上させることができ、せん断波の伝播解析や速度測定に対する時間解像度を向上させることができる。なお、送信ビームフォーミングおよび受信ビームフォーミングは例示した方法に限らず、せん断波の伝播状態を取得できるフレームレートで断層画像が得られるような任意のビームフォーミングを行ってよい。

（７）実施の形態および各変形例では、変位検出部１４において、プッシュパルス送信前の断層画像信号を基準断層画像信号として、プッシュパルス送信後の断層画像信号との間でパターンマッチングや相関処理等による変位検出を行うとしたが、本発明は必ずしもこの場合に限られない。例えば、ステップＳ２０における断層画像信号の取得を行わず、基準断層画像信号として、着目領域を撮像した、断層画像記憶部１８に予め記憶されている断層画像信号、または、着目領域の代表的な断層画像信号を用いてもよい。また、プッシュパルス送信後の断層画像信号のうち、一定時間離れた異なる時間における２つの断層画像信号の片方を基準断層画像信号として用いてもよい。

（８）実施の形態および各変形例では、超音波信号取得部１３が音響線信号を生成し、断層画像記憶部１８が音響線信号を記憶し、変位検出部１４は１画像分の音響線信号である断層画像信号を用いて変位の検出を行う場合について説明したが、本発明は必ずしもこの場合に限られない。例えば、超音波信号取得部１３は断層画像信号に対して包絡線検波、対数圧縮による輝度変換や画面表示用の画像データとなるよう幾何変換を行うことで、断層画像信号を断層画像データに変換し、断層画像記憶部１８は断層画像データを記憶するとしてもよい。このようにすることで、制御部１１は弾性画像の幾何変換を行わずにそのまま表示部３に表示することが可能となる。また、断層画像記憶部１８が記憶している断層画像データについても、対数圧縮による輝度変換や幾何変換をする必要がなく、そのまま表示部３に表示することができる。

（９）実施の形態、変形例１、２、４では、伝播解析部１５は、変位画像の全域に対して変位領域を抽出するとしたが、本発明は必ずしもこの場合に限られない。例えば、伝播解析部１５は、プッシュパルスの印加直後の変位画像に対して、プッシュパルスの印加位置の近傍のみ変位領域を抽出するとしてもよい。このようにすることで、演算量を低減させることができるとともに、明らかにせん断波が存在しない領域を抽出の対象としないことで、ノイズを波面として誤検出することを防ぐことができる。なお、変位領域の抽出範囲の限定はこの場合に限られず、例えば、あるｙ座標について全くせん断波が検出されない場合には、以降の時刻について当該ｙ座標を変位量位置の抽出範囲から除外してもよい。

また、同様に、伝播解析部１５は、時間フィルタリングにおいて、プッシュパルスの印加時刻を基に時間フィルタリング処理を行うとしてもよい。例えば、プッシュパルスの印加直後の波面画像に対して、プッシュパルスの印加位置の近傍にない波面をノイズとして削除してもよい。このようにすることで、明らかにせん断波が存在しえない領域に存在するノイズを波面として誤検出することを防ぐことができる。

また、伝播解析部１５は、空間フィルタリング処理と時間フィルタリング処理とを特定の順に行うのではなく、同時に行ってもよい。
（１０）変形例３では、伝播解析部１５は、相関処理により、変位画像から変位量がピークとなる時刻を検出する場合について説明したが、本発明は必ずしもこの場合に限られない。例えば、相関処理による相関値に対してしきい値を設け、しきい値を上回る相関値が検出された時間帯を検出するとしてもよい。このようにすることで、変位量がピークとなる時刻をピンポイントで検出するだけではなく、その周辺時刻を検出して波面を正確に検出することが可能となる。また、逆に、相関処理による相関値に対して下限しきい値を設け、下限しきい値を下回る相関値が検出された場合、ピークとして検出しないとしてもよい。このようにすることで、波面が通過していない時刻における揺動などを波面と誤検出することを防ぐことができる。

（１１）実施の形態および各変形例では、超音波診断装置は表示部３と接続される構成であるとしたが、本発明は必ずしもこの場合に限られない。例えば、超音波診断装置１は表示部３を内蔵しているとしてもよいし、または、超音波診断装置１は表示部３と接続されず、伝播解析部１５が生成し弾性画像記憶部２１に記憶されている弾性画像データを他の記憶媒体に保存、あるいは、ネットワークを通じて他の装置に出力するとしてもよい。

また、同様に、超音波診断装置は超音波探触子２を内蔵しているとしてもよいし、あるいは、超音波探触子２が超音波信号取得部１３を備え、超音波信号取得部１３を有さない超音波診断装置が超音波探触子２から音響線信号を取得するとしてもよい。
（１２）実施の形態および各変形例に係る超音波診断装置は、その構成要素の全部又は一部を、１チップ又は複数チップの集積回路で実現してもよいし、コンピュータのプログラムで実現してもよいし、その他どのような形態で実施してもよい。例えば、伝播解析部と評価部とを１チップで実現してもよいし、超音波信号取得部のみを１チップで実現し、変位検出部等を別のチップで実現してもよい。

集積回路で実現する場合、典型的には、ＬＳＩ（ＬａｒｇｅＳｃａｌｅＩｎｔｅｇｒａｔｉｏｎ）として実現される。ここでは、ＬＳＩとしたが、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路、又は汎用プロセッサで実現してもよい。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（ＦｉｅｌｄＰｒｏｇｒａｍｍａｂｌｅＧａｔｅＡｒｒａｙ）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサを利用してもよい。

さらには、半導体技術の進歩、又は派生する別技術により、ＬＳＩに置き換わる集積回路化の技術が登場すれば、当然その技術を用いて機能ブロックの集積化を行ってもよい。
また、各実施の形態および各変形例に係る超音波診断装置は、記憶媒体に書き込まれたプログラムと、プログラムを読み込んで実行するコンピュータとで実現されてもよい。記憶媒体は、メモリカード、ＣＤ−ＲＯＭなどいかなる記録媒体であってもよい。また、本発明に係る超音波診断装置は、ネットワークを経由してダウンロードされるプログラムと、プログラムをネットワークからダウンロードして実行するコンピュータとで実現されてもよい。

（１３）以上で説明した実施の形態は、いずれも本発明の好ましい一具体例を示すものである。実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、工程、工程の順序などは一例であり、本発明を限定する主旨ではない。また、実施の形態における構成要素のうち、本発明の最上位概念を示す独立請求項に記載されていない工程については、より好ましい形態を構成する任意の構成要素として説明される。

また、発明の理解の容易のため、上記各実施の形態で挙げた各図の構成要素の縮尺は実際のものと異なる場合がある。また本発明は上記各実施の形態の記載によって限定されるものではなく、本発明の要旨を逸脱しない範囲において適宜変更可能である。
さらに、超音波診断装置においては基板上に回路部品、リード線等の部材も存在するが、電気的配線、電気回路について当該技術分野における通常の知識に基づいて様々な態様を実施可能であり、本発明の説明として直接的には無関係のため、説明を省略している。尚、上記示した各図は模式図であり、必ずしも厳密に図示したものではない。

≪補足≫
（１）実施の形態に係る超音波診断装置は、超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後超音波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の、前記被検体内の着目領域における伝播状態を検出する超音波診断装置であって、送信プロファイルに基づいてプッシュパルスを送信するプッシュパルス送信部と、前記プッシュパルスに続き超音波を被検体へ複数回送信し、送信された超音波に対応する被検体からの反射超音波を受信して複数の受信信号を時系列に取得し、前記複数の受信信号が取得された時刻における前記着目領域内の組織の、前記プッシュパルスに起因するせん断波による変位をそれぞれ検出する変位検出部と、前記変位検出部による検出結果に基づいて、前記着目領域におけるせん断波の伝播状態を評価する評価部と、前記評価部が評価した結果に基づいて、前記送信プロファイルを調整するプッシュパルス調整部とを備え、前記プッシュパルス調整部が送信プロファイルを調整した場合に、前記プッシュパルス送信部は、調整された後の送信プロファイルに基づいて新たなプッシュパルスを送信し、前記変位検出部は、前記新たなプッシュパルスに起因するせん断波による変位を検出することを特徴とする。

また、実施の形態に係る超音波診断装置の制御方法は、超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後超音波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の、前記被検体内の着目領域における伝播状態を検出する超音波診断装置の制御方法であって、送信プロファイルに基づいてプッシュパルスを送信し、前記プッシュパルスに続き超音波を被検体へ複数回送信し、送信された超音波に対応する被検体からの反射超音波を受信して複数の受信信号を時系列に取得し、前記複数の受信信号が取得された時刻における前記着目領域内の組織の、前記プッシュパルスに起因するせん断波による変位をそれぞれ検出し、前記せん断波による変位の検出結果に基づいて、前記着目領域におけるせん断波の伝播状態を評価し、前記せん断波の伝播状態の評価結果に基づいて、前記送信プロファイルを調整し、前記送信プロファイルを調整した場合、さらに、調整された後の送信プロファイルに基づいて新たなプッシュパルスを送信し、前記新たなプッシュパルスに起因するせん断波による変位を検出することを特徴とする。

（２）また、上記（１）の超音波診断装置は、前記変位検出部が検出した前記複数の受信信号が取得された各時刻に対応する前記着目領域内の組織の変位を用いて、前記着目領域におけるせん断波の伝播状態を解析する伝播解析部をさらに備え、前記評価部は、前記伝播解析部が前記変位検出部による検出結果を解析した結果に基づいて、前記着目領域におけるせん断波の伝播状態を評価する、としてもよい。

このようにすることで、せん断波の向きや速度、伝播範囲等を解析した結果に基づいて、せん断波の伝播状態の評価やプッシュパルスの送信プロファイルの最適化を行うことができる。
（３）また、上記（２）の超音波診断装置は、前記伝播解析部は、前記変位検出部による検出結果に基づいて、前記複数の受信信号が取得された各時刻に対応する、前記着目領域内の各位置におけるせん断波の伝播速度を算出する、としてもよい。

このようにすることで、各時刻におけるせん断波の速度に基づいて、せん断波の伝播状態の評価やプッシュパルスの送信プロファイルの最適化を行うことができる。
（４）また、上記（１）〜（３）の超音波診断装置は、前記評価部は、前記変位検出部が検出した前記複数の受信信号が取得された各時刻に対応する前記着目領域内の組織の変位を用いて、前記着目領域内におけるせん断波の伝播状態を評価する、としてもよい。

このようにすることで、せん断波の向きや速度、伝播範囲等を組織の変位に直接基づいて評価を行うことができる。
（５）また、上記（４）の超音波診断装置は、前記評価部は、前記複数の受信信号が取得された時刻における前記着目領域内の各位置について、変位の有無に基づいてせん断波が伝播しているか否かを評価する、としてもよい。

このようにすることで、せん断波の伝播解析を行わなくても、せん断波の伝播範囲を評価することができる。
（６）また、上記（１）〜（５）の超音波診断装置は、前記送信プロファイルは、プッシュパルスが集中する位置、プッシュパルスの時間長、プッシュパルスの送信に用いる前記超音波探触子の素子数、又はプッシュパルスの送信回数の少なくとも１つを制御パラメータとして含み、前記プッシュパルス調整部は、前記制御パラメータのうち少なくとも１つを変更する、としてもよい。

（７）また、上記（６）の超音波診断装置は、前記送信プロファイルは、プッシュパルスが集中する位置を示す制御パラメータとして、プッシュパルスが集中する深さと、プッシュパルスの伝播する向きと前記超音波探触子との間の角度とを含む、としてもよい。
これらの構成により、プッシュパルスの焦点位置およびプッシュパルスのエネルギー、および、プッシュパルスの焦点の数を調整することができる。

（８）また、上記（６）または（７）の超音波診断装置は、前記プッシュパルス調整部は、プッシュパルスの送信回数を１から２に変更する場合、被検体内の異なる２つの位置のうち一方を第１のプッシュパルスが集中する位置として、他方を第２のプッシュパルスが集中する位置として選択する、としてもよい。
このようにすることで、プッシュパルス１つあたりのエネルギーを上昇させることなく、被検体組織への影響を抑制しながらせん断波の伝播範囲を拡大することができる。

（９）また、上記（１）〜（７）の超音波診断装置は、前記評価部は、前記着目領域においてせん断波の伝播の有無を評価し、前記プッシュパルス調整部は、せん断波が伝播していない領域が前記着目領域において占める位置に基づいて、当該せん断波が伝播していない領域にせん断波が到達するように送信プロファイルを調整する、としてもよい。
このようにすることで、せん断波が着目領域全体に行き渡るようにプッシュパルスを送信することができ、弾性値が取得できない領域の発生を抑制することができる。

（１０）また、上記（９）の超音波診断装置は、前記プッシュパルス調整部は、前記着目領域において、前記せん断波が伝播していない領域より浅い領域と深い領域とのいずれか一方のみに前記せん断波が伝播している場合、前記せん断波が伝播していない領域にプッシュパルスが集中する位置を近づける、としてもよい。
このようにすることで、例えば、プッシュパルスの焦点の深さが適切でない場合に、せん断波が着目領域全体に行き渡るようにプッシュパルスを送信することができる。

（１１）また、上記（９）の超音波診断装置は、前記プッシュパルス調整部は、前記着目領域において、前記せん断波が伝播していない領域より深い領域と浅い領域との双方に前記せん断波が伝播している場合、前記超音波探触子の素子配列方向にプッシュパルスが集中する位置を移動させる、としてもよい。
このようにすることで、例えば、せん断波の屈折等により非伝播領域が生じている場合に、せん断波が着目領域全体に行き渡るようにプッシュパルスを送信することができる。

（１２）また、上記（９）の超音波診断装置は、前記プッシュパルス調整部は、前記超音波探触子の素子配列方向において、前記せん断波が伝播していない領域が前記着目領域の一端を含み、かつ、前記せん断波が伝播していない領域より前記着目領域の他端側に前記せん断波が伝播している場合、プッシュパルスのエネルギーを上昇させる、前記せん断波が伝播していない領域の内部または近傍に新たなプッシュパルスが集中する位置を追加する、のいずれかを行う、としてもよい。

このようにすることで、例えば、せん断波が伝播途上に減衰している場合に、せん断波が着目領域全体に行き渡るようにプッシュパルスを送信することができる。
（１３）また、上記（９）の超音波診断装置は、前記プッシュパルス調整部は、前記超音波探触子の素子配列方向において、せん断波が伝播している領域に前記せん断波が伝播していない領域が挟まれている場合、前記せん断波が伝播していない領域の近傍に新たなプッシュパルスが集中する位置を追加する、としてもよい。

このようにすることで、例えば、反射源により反射したせん断波が伝播しているせん断波と干渉を起こしている場合や、反射源の内部にせん断波が伝播していない場合、せん断波が着目領域全体に行き渡るようにプッシュパルスを送信することができる。
（１４）また、上記（１）〜（１３）の超音波診断装置は、前記評価部は、前記着目領域におけるせん断波の伝播の一様性を評価し、前記プッシュパルス調整部は、せん断波が一様に伝播している領域内にプッシュパルスが集中する位置が存在するように、送信プロファイルを調整する、としてもよい。

このようにすることで、せん断波の屈折や反射が発生することを抑制し、せん断波速度の測定精度を向上させることができる。
（１５）また、上記（１）〜（１４）の超音波診断装置は、前記評価部は、前記着目領域におけるせん断波の伝播の一様性を評価し、前記プッシュパルス調整部は、せん断波の伝播状態が局所的に異なる領域の近傍にプッシュパルスが集中する位置が存在するように、送信プロファイルを調整する、としてもよい。

このようにすることで、せん断波の伝播状態が局所的に異なる領域に対応する、弾性値が周囲と異なる組織（例えば、腫瘍など）について、弾性測定の精度を向上させることができる。
（１６）また、上記（１）〜（１５）の超音波診断装置は、前記変位検出部による検出結果に基づいて、被検体内の各組織の弾性を計測する弾性計測部をさらに備える、としてもよい。

このようにすることで、プッシュパルス送信に基づいて、被検体の各組織の弾性を計測することができる。

本開示に係る超音波診断装置は、超音波を用いた組織の硬さの測定に有用である。そのため、組織の硬さの測定精度を向上させることが可能となり、医療診断機器等において高い利用可能性を持つ。

１超音波診断装置
２超音波探触子
３表示部
１１制御部
１２せん断波励起部
１３超音波信号取得部
１４変位検出部
１５伝播解析部
１６評価部
１７プッシュパルス調整部
１８断層画像記憶部
１９変位量記憶部
２０波面画像記憶部
２１弾性画像記憶部

Claims

超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後超音波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の、前記被検体内の着目領域における伝播状態を検出する超音波診断装置であって、
送信プロファイルに基づいてプッシュパルスを送信するプッシュパルス送信部と、
前記プッシュパルスに続き超音波を被検体へ複数回送信し、送信された超音波に対応する被検体からの反射超音波を受信して複数の受信信号を時系列に取得し、前記着目領域内の組織の、前記プッシュパルスに起因するせん断波による変位をそれぞれ検出する変位検出部と、
前記変位検出部による検出結果に基づいて、前記着目領域におけるせん断波の伝播の有無を評価する評価部と、
せん断波が伝播していない領域が前記着目領域において占める位置に基づいて、前記着目領域において、前記せん断波が伝播していない領域より深い領域と浅い領域との双方に前記せん断波が伝播している場合、前記超音波探触子の素子配列方向にプッシュパルスが集中する位置を移動させることで、当該せん断波が伝播していない領域にせん断波が到達するように送信プロファイルを調整するプッシュパルス調整部とを備え、
前記プッシュパルス調整部が送信プロファイルを調整した場合に、
前記プッシュパルス送信部は、調整された後の送信プロファイルに基づいて新たなプッシュパルスを送信し、
前記変位検出部は、前記新たなプッシュパルスに起因するせん断波による変位を検出する
ことを特徴とする超音波診断装置。
超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後超音波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の、前記被検体内の着目領域における伝播状態を検出する超音波診断装置であって、
送信プロファイルに基づいてプッシュパルスを送信するプッシュパルス送信部と、
前記プッシュパルスに続き超音波を被検体へ複数回送信し、送信された超音波に対応する被検体からの反射超音波を受信して複数の受信信号を時系列に取得し、前記着目領域内の組織の、前記プッシュパルスに起因するせん断波による変位をそれぞれ検出する変位検出部と、
前記変位検出部による検出結果に基づいて、前記着目領域におけるせん断波の伝播の有無を評価する評価部と、
せん断波が伝播していない領域が前記着目領域において占める位置に基づいて、前記超音波探触子の素子配列方向において、前記せん断波が伝播していない領域が前記着目領域の一端を含み、かつ、前記せん断波が伝播していない領域より前記着目領域の他端側に前記せん断波が伝播している場合、プッシュパルスのエネルギーを上昇させることで、当該せん断波が伝播していない領域にせん断波が到達するように送信プロファイルを調整するプッシュパルス調整部とを備え、
前記プッシュパルス調整部が送信プロファイルを調整した場合に、
前記プッシュパルス送信部は、調整された後の送信プロファイルに基づいて新たなプッシュパルスを送信し、
前記変位検出部は、前記新たなプッシュパルスに起因するせん断波による変位を検出する
ことを特徴とする超音波診断装置。
超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後超音波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の、前記被検体内の着目領域における伝播状態を検出する超音波診断装置であって、
送信プロファイルに基づいてプッシュパルスを送信するプッシュパルス送信部と、
前記プッシュパルスに続き超音波を被検体へ複数回送信し、送信された超音波に対応する被検体からの反射超音波を受信して複数の受信信号を時系列に取得し、前記着目領域内の組織の、前記プッシュパルスに起因するせん断波による変位をそれぞれ検出する変位検出部と、
前記変位検出部による検出結果に基づいて、前記着目領域におけるせん断波の伝播の一様性を評価する評価部と、
せん断波が一様に伝播している領域内にプッシュパルスが集中する位置が存在するように、送信プロファイルを調整するプッシュパルス調整部とを備え、
前記プッシュパルス調整部が送信プロファイルを調整した場合に、
前記プッシュパルス送信部は、調整された後の送信プロファイルに基づいて新たなプッシュパルスを送信し、
前記変位検出部は、前記新たなプッシュパルスに起因するせん断波による変位を検出する
ことを特徴とする超音波診断装置。
前記変位検出部が検出した前記複数の受信信号が取得された各時刻に対応する前記着目領域内の組織の変位を用いて、前記着目領域におけるせん断波の伝播状態を解析する伝播解析部をさらに備え、
前記評価部は、前記伝播解析部が前記変位検出部による検出結果を解析した結果に基づいて、前記着目領域におけるせん断波の伝播状態を評価する
ことを特徴とする請求項１から３のいずれか１項に記載の超音波診断装置。
前記伝播解析部は、前記変位検出部による検出結果に基づいて、前記複数の受信信号が取得された各時刻に対応する、前記着目領域内の各位置におけるせん断波の伝播速度を算出する
ことを特徴とする請求項４に記載の超音波診断装置。
前記評価部は、前記変位検出部が検出した前記複数の受信信号が取得された各時刻に対応する前記着目領域内の組織の変位を用いて、前記着目領域内におけるせん断波の伝播状態を評価する
ことを特徴とする請求項１から５のいずれか１項に記載の超音波診断装置。
前記評価部は、前記複数の受信信号が取得された時刻における前記着目領域内の各位置について、変位の有無に基づいてせん断波が伝播しているか否かを評価する
ことを特徴とする請求項６に記載の超音波診断装置。
前記送信プロファイルは、プッシュパルスが集中する位置、プッシュパルスの時間長、プッシュパルスの送信に用いる前記超音波探触子の素子数、又はプッシュパルスの送信回数の少なくとも１つを制御パラメータとして含み、
前記プッシュパルス調整部は、前記制御パラメータのうち少なくとも１つを変更する
ことを特徴とする請求項１から７のいずれか１項に記載の超音波診断装置。
前記送信プロファイルは、プッシュパルスが集中する位置を示す制御パラメータとして、プッシュパルスが集中する深さと、プッシュパルスの伝播する向きと前記超音波探触子との間の角度とを含む
ことを特徴とする請求項８に記載の超音波診断装置。
前記プッシュパルス調整部は、プッシュパルスの送信回数を１から２に変更する場合、被検体内の異なる２つの位置のうち一方を第１のプッシュパルスが集中する位置として、他方を第２のプッシュパルスが集中する位置として選択する
ことを特徴とする請求項８または９に記載の超音波診断装置。
前記変位検出部による検出結果に基づいて、被検体内の各組織の弾性を計測する弾性計測部をさらに備える
ことを特徴とする請求項１から１０のいずれか１項に記載の超音波診断装置。
超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後超音波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の、前記被検体内の着目領域における伝播状態を検出する超音波診断装置の制御方法であって、
送信プロファイルに基づいてプッシュパルスを送信し、
前記プッシュパルスに続き超音波を被検体へ複数回送信し、送信された超音波に対応する被検体からの反射超音波を受信して複数の受信信号を時系列に取得し、前記着目領域内の組織の、前記プッシュパルスに起因するせん断波による変位の有無をそれぞれ検出し、
前記せん断波による変位の検出結果に基づいて、せん断波が伝播していない領域が前記着目領域において占める位置に基づいて、前記着目領域において、前記せん断波が伝播していない領域より深い領域と浅い領域との双方に前記せん断波が伝播している場合、前記超音波探触子の素子配列方向にプッシュパルスが集中する位置を移動させることで、当該せん断波が伝播していない領域にせん断波が到達するように前記着目領域におけるせん断波の伝播状態を評価し、
前記せん断波の伝播状態の評価結果に基づいて、前記送信プロファイルを調整し、
前記送信プロファイルを調整した場合、さらに、
調整された後の送信プロファイルに基づいて新たなプッシュパルスを送信し、
前記新たなプッシュパルスに起因するせん断波による変位を検出する
ことを特徴とする制御方法。
超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後超音波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の、前記被検体内の着目領域における伝播状態を検出する超音波診断装置の制御方法であって、
送信プロファイルに基づいてプッシュパルスを送信し、
前記プッシュパルスに続き超音波を被検体へ複数回送信し、送信された超音波に対応する被検体からの反射超音波を受信して複数の受信信号を時系列に取得し、前記着目領域内の組織の、前記プッシュパルスに起因するせん断波による変位の有無をそれぞれ検出し、
前記せん断波による変位の検出結果に基づいて、せん断波が伝播していない領域が前記着目領域において占める位置に基づいて、前記超音波探触子の素子配列方向において、前記せん断波が伝播していない領域が前記着目領域の一端を含み、かつ、前記せん断波が伝播していない領域より前記着目領域の他端側に前記せん断波が伝播している場合、プッシュパルスのエネルギーを上昇させることで、当該せん断波が伝播していない領域にせん断波が到達するように前記着目領域におけるせん断波の伝播状態を評価し、
前記せん断波の伝播状態の評価結果に基づいて、前記送信プロファイルを調整し、
前記送信プロファイルを調整した場合、さらに、
調整された後の送信プロファイルに基づいて新たなプッシュパルスを送信し、
前記新たなプッシュパルスに起因するせん断波による変位を検出する
ことを特徴とする制御方法。
超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して特定部位にある組織を物理的に押圧した後超音波の送受信を被検体内に対して繰り返し行うことで、押圧された特定部位の組織を振動源とするせん断波の、前記被検体内の着目領域における伝播状態を検出する超音波診断装置の制御方法であって、
送信プロファイルに基づいてプッシュパルスを送信し、
前記プッシュパルスに続き超音波を被検体へ複数回送信し、送信された超音波に対応する被検体からの反射超音波を受信して複数の受信信号を時系列に取得し、前記着目領域内の組織の、前記プッシュパルスに起因するせん断波の伝播の一様性をそれぞれ検出し、
せん断波が一様に伝播している領域内にプッシュパルスが集中する位置が存在するように、前記着目領域におけるせん断波の伝播状態を評価し、
前記せん断波の伝播状態の評価結果に基づいて、前記送信プロファイルを調整し、
前記送信プロファイルを調整した場合、さらに、
調整された後の送信プロファイルに基づいて新たなプッシュパルスを送信し、
前記新たなプッシュパルスに起因するせん断波による変位を検出する
ことを特徴とする制御方法。