(57)【特許請求の範囲】
【請求項１】
超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して、特定部位にある組織を物理的に押圧し、押圧により特定部位の組織が振動源となってせん断波が隣接組織へと伝播する過程を、プッシュパルスに続いて超音波の被検体への送信と、被検体からの反射超音波の受信とを繰り返すことによって検出する超音波診断装置であって、
前記プッシュパルスに続く超音波の被検体への送信と被検体からの反射超音波の受信により時系列に複数の受信信号を得て、前記複数の受信信号が取得された各時刻における前記被検体内の組織の変位をそれぞれ検出する変位検出部と、
前記複数の受信信号が取得された各時刻における組織の変位に基づいて、当該時刻における前記せん断波の波面を、変位量が閾値よりも大きい変位領域を抽出し、抽出した各変位領域に細線化処理を施すことで抽出する波面抽出部と、
前記複数の受信信号が取得された各時刻における波面を波面画像として生成する波面画像加工部と
を備えることを特徴とする超音波診断装置。
【請求項２】
超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して、特定部位にある組織を物理的に押圧し、押圧により特定部位の組織が振動源となってせん断波が隣接組織へと伝播する過程を、プッシュパルスに続いて超音波の被検体への送信と、被検体からの反射超音波の受信とを繰り返すことによって検出する超音波診断装置であって、
前記プッシュパルスに続く超音波の被検体への送信と被検体からの反射超音波の受信により時系列に複数の受信信号を得て、前記複数の受信信号が取得された各時刻における前記被検体内の組織の変位をそれぞれ検出する変位検出部と、
前記複数の受信信号が取得された各時刻における組織の変位に基づいて、当該時刻における前記せん断波の波面を、前記複数の受信信号が示す被検体内の各位置における変位が極値となる時刻を相関処理によって検出し、検出された被検体内の各位置における変位が極値となる時刻に基づいて前記受信信号のそれぞれにおける変位領域を抽出することで抽出する波面抽出部と、
前記複数の受信信号が取得された各時刻における波面を波面画像として生成する波面画像加工部と
を備えることを特徴とする超音波診断装置。
【請求項３】
前記波面画像加工部は、前記波面の時間変化を表示する画像をさらに生成する
ことを特徴とする請求項１または２記載の超音波診断装置。
【請求項４】
前記波面の時間変化を表示する画像は動画である
ことを特徴とする請求項３に記載の超音波診断装置。
【請求項５】
前記波面画像加工部は、前記波面画像を、前記被検体内の組織の各領域の弾性を示す弾性画像に重畳して出力する
ことを特徴とする請求項１乃至４のいずれか１項に記載の超音波診断装置。
【請求項６】
前記波面画像加工部は、前記波面画像を、前記被検体内の組織に対応した断層画像に重畳して出力する
ことを特徴とする請求項１乃至４のいずれか１項に記載の超音波診断装置。
【請求項７】
前記波面画像加工部は、前記波面抽出部が抽出した波面の表示態様を加工し、前記受信信号が取得された異なる時刻のそれぞれにおける波面について、波面の表示態様が前記時刻ごとに互いに異ならせるように波面の表示態様を加工する
ことを特徴とする請求項１乃至６のいずれか１項に記載の超音波診断装置。
【請求項８】
前記波面画像加工部は、時刻の異なる複数の波面に対し、線の太さ、色、輝度、および、種類のうち少なくとも１つを変化させることで表示態様を異ならせる
ことを特徴とする請求項７に記載の超音波診断装置。
【請求項９】
前記波面画像加工部は、波面のそれぞれを、伝播方向と伝播速度とをそれぞれ向きと長さで表現した矢印で表現し、時刻の異なる複数の波面に対し、矢印の太さ、色、および、輝度のうち少なくとも１つを変化させることで表示態様を異ならせる
ことを特徴とする請求項７に記載の超音波診断装置。
【請求項１０】
前記波面画像加工部は、前記波面抽出部が抽出した波面のうち、進行方向が同一であるせん断波に係る波面のみを表示する
ことを特徴とする請求項１乃至９のいずれか１項に記載の超音波診断装置。
【請求項１１】
前記波面画像加工部は、前記せん断波の進行方向ごとに、波面の表示態様を異ならせる
ことを特徴とする請求項１乃至９のいずれか１項に記載の超音波診断装置。
【請求項１２】
前記波面抽出部は、さらに、前記せん断波の進行方向が変化している界面を抽出し、
前記波面画像加工部は、さらに、前記界面を前記波面に重畳する
ことを特徴とする請求項１乃至１１のいずれか１項に記載の超音波診断装置。
【請求項１３】
前記プッシュパルスを第１のプッシュパルスとしたとき、さらに第２のプッシュパルスを前記超音波探触子に送信させるせん断波励起部をさらに備え、
前記波面画像加工部は、前記第１のプッシュパルスに起因するせん断波の波面と、前記第２のプッシュパルスに起因するせん断波の波面とで、表示態様を異ならせる
ことを特徴とする請求項１乃至１２のいずれか１項に記載の超音波診断装置。
【請求項１４】
前記波面抽出部が検出した時刻の異なる同一の波面の移動量を、当該時刻におけるせん断波の伝播速度とすることで、前記被検体内の組織のうち、前記波面が通過した各領域の弾性をそれぞれ推定する弾性推定部
をさらに備えることを特徴とする請求項１乃至１３のいずれか１項に記載の超音波診断装置。
【請求項１５】
超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して、特定部位にある組織を物理的に押圧し、押圧により特定部位の組織が振動源となってせん断波が隣接組織へと伝播する過程を、プッシュパルスに続いて超音波の被検体への送信と、被検体からの反射超音波の受信とを繰り返すことによって検出する超音波診断装置の制御方法であって、
前記プッシュパルスに続く超音波の被検体への送信と被検体からの反射超音波の受信により時系列に複数の受信信号を得て、前記複数の受信信号が取得された各時刻における前記被検体内の組織の変位をそれぞれ検出し、
前記複数の受信信号が取得された各時刻における組織の変位に基づいて、当該時刻における前記せん断波の波面を、変位量が閾値よりも大きい変位領域を抽出し、抽出した各変位領域に細線化処理を施すことで抽出し、
前記複数の受信信号が取得された各時刻における波面を波面画像として生成する
ことを特徴とする制御方法。
【請求項１６】
超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して、特定部位にある組織を物理的に押圧し、押圧により特定部位の組織が振動源となってせん断波が隣接組織へと伝播する過程を、プッシュパルスに続いて超音波の被検体への送信と、被検体からの反射超音波の受信とを繰り返すことによって検出する超音波診断装置の制御方法であって、
前記プッシュパルスに続く超音波の被検体への送信と被検体からの反射超音波の受信により時系列に複数の受信信号を得て、前記複数の受信信号が取得された各時刻における前記被検体内の組織の変位をそれぞれ検出し、
前記複数の受信信号が取得された各時刻における組織の変位に基づいて、当該時刻における前記せん断波の波面を、前記複数の受信信号が示す被検体内の各位置における変位が極値となる時刻を相関処理によって検出し、検出された被検体内の各位置における変位が極値となる時刻に基づいて前記受信信号のそれぞれにおける変位領域を抽出することで抽出し、
前記複数の受信信号が取得された各時刻における波面を波面画像として生成する
ことを特徴とする制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置における画像化処理に関し、特に、せん断波を用いた硬さ測定に関する。
【背景技術】
【０００２】
近年、被検体内の組織の硬さを評価する機能を有する超音波診断装置が普及してきている。超音波診断装置を用いて、硬さを評価する方法は大きく分けて２つある。１つは、超音波探触子を用いて被検体内の組織を体表から圧迫した後に解放し、圧迫に対する被検体内の組織の歪みの大きさから、被検体内の組織の相対的な硬さを評価する方法である。この方法では、組織が周辺組織と比べて硬いか軟らかいかという、被検体内の相対的な硬さを評価できる。
【０００３】
もう１つは、被検体内の着目領域（ＲＯＩ；Ｒｅｇｉｏｎ  ｏｆ  Ｉｎｔｅｒｅｓｔ）にせん断波（ｓｈｅａｒ  ｗａｖｅ）を生じさせ、着目領域内の組織の変位を時系列に取得することで、せん断波の伝播速度を計測する方法である。せん断波の伝播速度は組織の弾性率に応じて変化するため、この方法では組織の絶対的な硬さを評価することができる。せん断波を生じさせる手法としては、例えば、ＡＲＦＩ（Ａｃｏｕｓｔｉｃ  Ｒａｄｉａｔｉｏｎ  Ｆｏｒｃｅ  Ｉｍｐｕｌｓｅ）と呼ばれる、超音波の音圧により被検体内の組織に変位を発生させる手法が用いられる。
【０００４】
しかしながら、この方法において、せん断波の速度が正常に算出できない場合がある。例えば、着目領域内の組織の変位が不規則である場合である。このような場合、組織全体の動きに起因する変位が正確に抽出できないため、せん断波に起因する変位に組織全体の動きに起因する変位が混じって、せん断波の速度の算出精度が低下する。また、もう１つの例は、せん断波が組織の境界において屈折や反射を起こしている場合である。従来、組織内の各点において、変位が最初に最大となる時刻を、その点におけるせん断波の最初の波面の到達時刻であるとみなしてせん断波の速度を算出する。すなわち、せん断波が発生源から遠ざかるように進むことを前提としているため、それ以外の場合には、せん断波の速度の算出が行えない、または、算出結果の精度が悪い場合がある。
【０００５】
このような問題に対処するため、従来、組織内のある点において、せん断波の到達時刻がその周囲の点における到達時刻と大きく異なる場合やせん断波の速度が算出不能である場合、その点におけるせん断波の速度の算出を行わない、または、その点におけるせん断波の速度の信頼性が低いと検査者に提示することが提案されている。また、例えば、特許文献１の技術では、せん断波の到達時間に応じて色相を割り当てたせん断波到達画像でせん断波の到達時刻の分布を示すことで、せん断波の速度を視覚的に示すとともに、せん断波の速度についての信頼度が低い領域を明示している。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１４−２６０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、せん断波の速度の信頼性を検査者に示す従来の技術では、どのような理由でせん断波の速度が取得できなかったのかを具体的に知ることができない。特許文献１の技術では、せん断波が発生源から遠ざかるように進む前提でせん断波到達画像を生成し、まだら模様になるなど画像の乱れた領域を、せん断波の速度についての信頼度が低い領域としている。そのため、例えば、せん断波の速度の信頼度が低い理由が、せん断波の反射等である場合、その発生個所や、それによるせん断波の速度の測定結果への影響を具体的に知ることができない。したがって、従来の技術によれば、検査者は、どうすればせん断波の速度が取得できるのか分からないまま、手探りでせん断波の速度の再測定を行う必要がある。
【０００８】
本開示は、上記課題に鑑みてなされたものであり、せん断波が適切に伝播するために必要な情報を検査者に提供し、検査者の使い勝手を向上することを目的とする。
【課題を解決するための手段】
【０００９】
本開示の一態様に係る超音波診断装置は、超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して、特定部位にある組織を物理的に押圧し、押圧により特定部位の組織が振動源となってせん断波が隣接組織へと伝播する過程を、プッシュパルスに続いて超音波の被検体への送信と、被検体からの反射超音波の受信とを繰り返すことによって検出する超音波診断装置であって、前記プッシュパルスに続く超音波の被検体への送信と被検体からの反射超音波の受信により時系列に複数の受信信号を得て、前記複数の受信信号が取得された各時刻における前記被検体内の組織の変位をそれぞれ検出する変位検出部と、前記複数の受信信号が取得された各時刻における組織の変位に基づいて、当該時刻における前記せん断波の波面を、変位量が閾値よりも大きい変位領域を抽出し、抽出した各変位領域に細線化処理を施すことで抽出する波面抽出部と、前記複数の受信信号が取得された各時刻における波面を波面画像として生成する波面画像加工部とを備えることを特徴とする。
【発明の効果】
【００１０】
本開示によれば、上記構成により、せん断波が伝播している間の各時刻の波面を画像として検査者が確認することができるので、検査者は、せん断波の波面が反射や屈折をすることなく伝播したか否かを確認することができる。さらに、検査者は、せん断波の反射や屈折が発生した場合には、どの場所でそのような事態が起きたのかを、波面の形状や、時刻ごとの波面の位置およびそれが示す波面の伝播速度や方向に基づいて容易に推測できる。したがって、検査者は、せん断波の屈折や反射により被検体内の組織の硬さを取得できなかった場合に、せん断波の伝播の様子をもとに、超音波探触子の位置を変える、被検者の姿勢を変える、プッシュパルスを印加する位置を変える、などの適切な対処を行って再検査を行うことができる。
【図面の簡単な説明】
【００１１】
【図１】実施の形態１に係る超音波診断装置１のブロック図である。
【図２】実施の形態１に係る波面抽出部１５の機能ブロック図である。
【図３】実施の形態１に係る超音波診断装置１全体の動作を示すフローチャートである。
【図４】せん断波の発生と伝播の様子を示す模式図である。
【図５】実施の形態１に係る波面抽出の動作を示すフローチャートである。
【図６】実施の形態１に係る波面抽出の動作を示す模式図である。
【図７】実施の形態１に係る波面の表示例を示す図である。
【図８】実施の形態１に係る波面の伝播の様子を示す図である。
【図９】変形例１に係る波面抽出部２０の機能ブロック図である。
【図１０】変形例１に係る波面抽出の動作を示すフローチャートである。
【図１１】変形例１に係る波面抽出の動作を示す模式図である。
【図１２】変形例２に係る波面分類の動作を示すフローチャートである。
【図１３】変形例２に係る波面画像の一例である。
【図１４】変形例３に係る波面画像の一例である。
【図１５】実施の形態２に係る超音波診断装置４のブロック図である。
【図１６】実施の形態２に係る弾性推定の処理を示す模式図である。
【図１７】実施の形態２に係る超音波診断装置４の波面画像の表示形態の一例である。
【発明を実施するための形態】
【００１２】
以下、本発明の実施の形態について、図面を参照しながら説明する。
≪実施の形態１≫
実施の形態１に係る超音波診断装置１のブロック図を図１に示す。超音波診断装置１は、制御部１１、せん断波励起部１２、超音波信号取得部１３、変位検出部１４、波面抽出部１５、波面画像加工部１６、断層画像記憶部１７、変位量記憶部１８、波面画像記憶部１９を備える。また、制御部１１には、超音波探触子２と表示部３とがそれぞれ接続可能に構成されている。図１は超音波診断装置１に超音波探触子２、表示部３が接続された状態を示している。
【００１３】
超音波探触子２は、例えば、一次元方向に配列された複数の振動子（不図示）を有する。各振動子は、例えば、ＰＺＴ（チタン酸ジルコン酸鉛）により構成される。超音波探触子２は、せん断波励起部１２で生成された電気信号（以下、「ＡＲＦＩ駆動信号」と呼ぶ）、または、超音波信号取得部１３で生成された電気信号（以下、「送信駆動信号」と呼ぶ）を制御部１１から受け取り、超音波に変換する。超音波探触子２は、超音波探触子２の振動子側外表面を被検体の皮膚表面など表面に接触させた状態で、ＡＲＦＩ駆動信号又は送信駆動信号から変換され、複数の振動子から発せられる複数の超音波からなる超音波ビームを被検体内の測定対象に向けて送信する。そして、超音波探触子２は、送信駆動信号に基づく送信超音波に対する、測定対象からの複数の反射超音波を受信し、複数の振動子によりこれら反射超音波をそれぞれ電気信号（以下、「素子受信信号」と呼ぶ）に変換し、制御部１１を介して素子受信信号を超音波信号取得部１３に供給する。なお、ここではせん断波励起部１２及び超音波信号取得部１３を別構成として説明しているが、超音波信号取得部１３の送信駆動信号を生成する構成と同一の構成により、ＡＲＦＩ駆動信号を生成してもよい。
【００１４】
せん断波励起部１２は、プッシュパルスを超音波探触子２に送出させるための電気信号であるＡＲＦＩ駆動信号を生成する。プッシュパルスとは、被検体内にせん断波を発生させるため、被検体内の組織に変位をもたらすためのパルス状の超音波である。具体的には、被検体内の着目領域内のある１点を焦点とする、波数が後述する送信超音波より多い超音波である。したがって、ＡＲＦＩ駆動信号は、超音波探触子２を構成する各振動素子から送出される超音波が焦点に届くように生成される、パルス状の電気信号である。
【００１５】
超音波信号取得部１３は、送信超音波を超音波探触子２に送出させるための電気信号である送信駆動信号を生成する。送信駆動信号は、超音波探触子２を構成する各振動素子から送出される送信超音波が送信フォーカス点に同時に届くように生成される、振動素子ごとにタイミングの異なるパルス状の電気信号である。また、超音波信号取得部１３は、反射超音波に基づく素子受信信号に整相加算を行って音響線信号を生成する。超音波信号取得部１３は、生成した音響線信号を、制御部１１を介して断層画像記憶部１７に出力する。
【００１６】
変位検出部１４は、変位検出の対象となる１つの断層画像に係る複数の音響線信号（以下、「断層画像信号」と呼ぶ）と、基準となる１つの断層画像に係る複数の音響線信号（以下、「基準断層画像信号」と呼ぶ）とを、制御部１１を介して断層画像記憶部１７から取得する。基準断層画像信号とは、断層画像信号からせん断波による変位を抽出するために用いるものであり、具体的には、プッシュパルス送出前に着目領域を撮像した断層画像信号である。そして、変位検出部１４は、断層画像信号と基準断層画像信号との差分から、断層画像信号の各画素の変位を検出し、変位を各画素の座標と関連付けて変位画像を生成する。変位検出部１４は、生成した変位画像を、制御部１１を介して変位量記憶部１８に出力する。
【００１７】
波面抽出部１５は、変位画像を、制御部１１を介して変位量記憶部１８から取得する。波面抽出部１５は、変位画像から、せん断波の波面を抽出し、１の変位画像に係る波面を１の波面画像データとして生成する。図２に、波面抽出部１５の機能ブロック図を示す。波面抽出部１５は、変位領域抽出部３１、細線化部３２、空間フィルタリング部３３、時間フィルタリング部３４を備える。
【００１８】
変位領域抽出部３１は、変位画像から、変位の大きい変位領域を抽出する。細線化部３２は、変位領域に細線化処理を施して波面を抽出し、波面画像データを生成する。空間フィルタリング部３３は、波面の形状を用いて、波面画像において誤検出された波面を取り除く。時間フィルタリング部３４は、複数の波面画像から、波面の時間推移を検出し、誤検出された波面を取り除く。詳細は後述する。波面抽出部１５は、生成した波面画像データを、制御部を介して波面画像記憶部１９に出力する。
【００１９】
波面画像加工部１６は、波面画像データを、制御部を介して波面画像記憶部１９から取得する。波面画像加工部１６は、同一時刻の波面が同一の表示態様となるように各波面の表示態様を加工する。波面画像加工部１６は、加工後の波面画像データを出力する。
断層画像記憶部１７、変位量記憶部１８、波面画像記憶部１９は、それぞれ、断層画像、変位画像、波面画像データを記憶する。断層画像記憶部１７、変位量記憶部１８、波面画像記憶部１９のそれぞれは、例えば、ＲＡＭ、フラッシュメモリ、ハードディスクなどの記憶媒体で実現される。
【００２０】
制御部１１は、上述したような各構成要素の制御に加え、波面画像加工部１６が加工した複数の波面画像を動画として表示部３に出力することもできる。
表示部３は、制御部１１から出力される画像や動画を表示する。表示部３は、例えば、液晶ディスプレイ、有機ＥＬディスプレイなどの表示装置で実現される。
制御部１１、せん断波励起部１２、超音波信号取得部１３、変位検出部１４、波面抽出部１５、波面画像加工部１６のそれぞれは、例えば、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｇｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）などのハードウェアにより実現される。なお、制御部１１、せん断波励起部１２、超音波信号取得部１３、変位検出部１４、波面抽出部１５、波面画像加工部１６は、それぞれ個別に、または、２以上を１まとめとして、メモリと、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）などのプログラマブルデバイスとソフトウェアで実現されてもよい。
【００２１】
≪動作≫
実施の形態１に係る超音波診断装置１の動作について説明する。図３は、超音波診断装置１全体の動作を示すフローチャートである。
まず、制御部１１が、被検体内に着目領域を設定する（ステップＳ１０）。着目領域を設定する方法は、例えば、表示部３に断層画像記憶部１７に記録されている最新の断層画像を表示し、タッチパネル、マウス、トラックボールなどの入力部（図示しない）を通して検査者に着目領域を指定させる。なお、着目領域の設定方法はこの場合に限られず、例えば、断層画像の全域を着目領域としてもよいし、あるいは、断層画像の中央部分を含む一定範囲を着目領域としてもよい。
【００２２】
次に、着目領域に超音波を送受信し、取得した受信信号を記憶する（ステップＳ２０）。具体的には、次のような動作となる。まず、以下のように送信イベントを行う。最初に、超音波信号取得部１３が、パルス状の送信信号を生成する。次に、超音波信号取得部１３は、送信信号に対して、超音波探触子２の素子毎の遅延時間を設定する送信ビームフォーミングを行い、超音波探触子２の各素子に対応する複数の送信駆動信号を生成する。超音波探触子２の各振動子は対応する送信駆動信号を超音波に変換することで、超音波ビームが被検体内に送出される。次に、超音波探触子２の各振動子は、被検体内から反射された反射超音波を取得して素子受信信号に変換する。超音波信号取得部１３は、素子受信信号に整相加算を行い、音響線信号を生成する。制御部１１は、送信イベントごとに超音波信号取得部１３から音響線信号を取得し、１つの断層画像を構成する複数の音響線信号を、断層画像信号として、断層画像記憶部１７に記憶する。
【００２３】
次に、プッシュパルスを送信する（ステップＳ３０）。具体的には、せん断波励起部１２が、パルス状のＡＲＦＩ信号を生成する。次に、せん断波励起部１２は、ＡＲＦＩ信号に対して、超音波探触子２の素子毎の遅延時間を設定する送信ビームフォーミングを行い、超音波探触子２の各素子に対応する複数のＡＲＦＩ駆動信号を生成する。超音波探触子２の各振動子は対応する送信駆動信号を超音波に変換することで、プッシュパルスが被検体内に送出される。
【００２４】
ここで、プッシュパルスによるせん断波の生成について、図４の模式図を用いて説明する。図４（ａ）は、プッシュパルス印加前の着目領域の模式図である。図４（ａ）〜（ｅ）において、個々の“○”は、着目領域における被検体内の組織の一部を、破線の交点は、負荷がない場合の組織”○“の中心位置を、それぞれ示している。ここで、超音波探触子２を皮膚表面１００に密接させた状態で焦点１０１に対してプッシュパルスを印加すると、図４（ｂ）の模式図に示すように、焦点１０１に位置していた組織１３２が、プッシュパルスの進行方向に押されて移動する。また、組織１３２からプッシュパルスの進行方向側にある組織１３３は、組織１３２に押されてプッシュパルスの進行方向に移動する。次に、プッシュパルスの送信が終了すると、組織１３２、１３３が元の位置に復元しようとするので、図４（ｃ）の模式図に示すように、組織１３１〜１３３がプッシュパルスの進行方向に沿った振動を開始する。すると、図４（ｄ）の模式図に示すように、振動が組織１３１〜１３３に隣接する、組織１２１〜１２３および組織１４１〜１４３に伝播する。さらに、図４（ｅ）の模式図に示すように、振動がさらに組織１１１〜１１３および組織１５１〜１５３に伝播する。したがって、被検体内において、振動が振動の方向と直交する向きに伝播する。すなわち、せん断波がプッシュパルスの印加場所に発生し、被検体内を伝播する。
【００２５】
図３に戻って説明を続ける。次に、着目領域に超音波を複数回送受信し、取得した複数の超音波信号を記憶する（ステップＳ４０）。具体的には、プッシュパルスの送信終了の直後から、例えば、秒間１万回、ステップＳ２０と同じ動作を繰り返し行う。これにより、せん断波の発生直後から伝播が終わるまでの間、被検体の断層画像を繰り返し取得する。
【００２６】
次に、波面を抽出する（ステップＳ５０）。詳しくは、図５のフローチャートを用いて説明する。図５は、波面の抽出動作の詳細を示すフローチャートである。
まず、各画素の変位を検出する（ステップＳ５１）。具体的には、変位検出部１４が、ステップＳ２０で断層画像記憶部１７に記憶された断層画像信号を、基準断層信号として取得する。上述したように、基準断層画像信号は、プッシュパルスの送出前、すなわち、せん断波の発生前に取得された断層画像信号である。次に、変位検出部１４は、ステップＳ４０で断層画像記憶部１７に記憶された各断層画像信号に対し、基準断層画像信号との差分から、当該断層画像信号が取得された時刻における、各画素の変位を検出する。具体的には、例えば、断層画像信号を８ピクセル×８ピクセルなどの所定の大きさの領域に分割し、各領域と基準断層画像信号とをパターンマッチングすることで、断層画像信号の各画素の変位を検出する。これにより、各断層画像信号の各画素に対応する被検体の組織が、プッシュパルスまたはせん断波によってどれだけ動いたかが変位として算出される。なお、変位の検出方法はパターンマッチングに限られず、例えば、断層画像信号と基準断層画像信号との相関処理など、２つの断層画像信号間の動き量を検出する任意の技術を用いてよい。変位検出部１４は、１の断層画像に係る各画素の変位を当該画素の座標と対応付けることで変位画像を生成し、生成した変位画像を変位量記憶部１８に出力する。
【００２７】
次に、変位が大きい変位領域を抽出する（ステップＳ５２）。変位領域抽出部３１は、変位画像を変位量記憶部１８から取得し、変位の大きい変位領域を抽出する。以下、図６の模式図を用いて説明する。図６（ａ）は、変位画像の一例を示している。図４と同じく、図中の“○”は着目領域における被検体内の組織の一部を示しており、プッシュパルスを印加する前の位置は破線の交点である。また、ｘ軸は超音波探触子２における素子の並ぶ方向、ｙ軸は、被検体の深さ方向である。波面抽出部１５は、ｙ座標ごとに変位量δを座標ｘの関数として、動的閾値を用いることで変位量δが大きい領域を抽出する。また、ｘ座標ごとに変位量δを座標ｙの関数として、動的閾値を用いて、ある閾値を超える領域を変位量δが大きい領域として抽出する。動的閾値とは、対象領域内について信号解析又は画像解析を行って閾値を決定することである。閾値は一定値ではなく、対象領域の信号の幅や最大値などによって異なる値となる。図６（ａ）に、ｙ＝ｙ1の直線２１０上における変位量をプロットしたグラフ２１１と、ｘ＝ｘ1の直線２２０上における変位量をプロットしたグラフ２２１とを示す。これにより、例えば、変位量δが閾値より大きな変位領域２３０が抽出できる。
【００２８】
次に、細線化部３２は、変位領域に細線化処理をおこなって波面を抽出する（ステップＳ５３）。図６（ｂ）の模式図に示している変位領域２４０、２５０は、それぞれ、ステップＳ５２において変位領域として抽出された領域である。細線化部３２は、例えば、Ｈｉｌｄｉｔｃｈの細線化アルゴリズムを用いて、波面を抽出する。例えば、図６（ｂ）の模式図において、変位領域２４０から波面２４１が、変位領域２５０から波面２５１が、それぞれ抽出される。なお、細線化のアルゴリズムはＨｉｌｄｉｔｃｈに限らず、任意の細線化アルゴリズムを用いてよい。また、各変位領域に対して、変位量δが閾値以下の座標を変位領域から取り除く処理を、変位領域が幅１ピクセルの線になるまで、閾値を大きくしながら繰り返し行ってもよい。細線化部３２は、抽出した波面を波面画像データとして空間フィルタリング部３３に出力する。
【００２９】
次に、空間フィルタリング部３３は、細線化部３２が生成した波面画像データに対して空間フィルタリングを行い、長さが短い波面を除去する（ステップＳ５４）。例えば、ステップＳ５３で抽出した各波面の長さを検出し、全ての波面の長さの平均値の１／２よりも長さが短い波面を、ノイズとして削除する。具体的には、図６（ｃ）の波面画像に示すように、波面２６１〜２６４の長さの平均値を算出し、それよりも短い波面２６３、２６４を、ノイズとして消去する。これにより、誤検出された波面を消去できる。
【００３０】
波面抽出部１５は、ステップＳ５１〜Ｓ５４の動作を、全ての変位画像に対して行う（ステップＳ５５）。これにより、変位画像に対して１対１で波面画像データが生成される。
次に、時間フィルタリング部３４は、複数の波面画像データに対して時間フィルタリングを行い、伝播していない波面を除去する（ステップＳ５６）。具体的には、時間的に連続する２以上の波面画像において、波面位置の時間変化を検出し、速度が異常である波面をノイズとして除去する。波面抽出部１５は、例えば、時刻ｔ＝ｔ1の波面画像２７０、時刻ｔ＝ｔ1＋Δｔの波面画像２８０、時刻ｔ＝ｔ1＋２Δｔの波面画像２９０との間で、波面位置の時間変化を検出する。例えば、波面２７１に対して、波面画像２８０のうち、波面２７１と同じ位置を中心に、波面と垂直な向き（図６においてはｘ軸方向）にΔｔの間にせん断波が移動しうる領域２７６で、波面２７１との相関処理を行う。このとき、波面２７１のｘ軸の正方向（図の右側）と負方向（図の左側）の双方を含む範囲内で相関処理を行う。これは、透過波と反射波の両方を検出するためである。これにより、波面２７１の移動先が波面画像２８０内の波面２８１であると検出し、時間Δｔにおける波面２７１の移動距離を算出する。同様に、波面２７２、２７３のそれぞれについて、波面画像２８０において当該波面と同じ位置を中心に、波面と垂直な向きにΔｔの間にせん断波が移動しうる領域で相関処理を行う。これにより、波面２７２が波面２８３の位置に、波面２７３が波面２８２の位置に、それぞれ移動したことを検出する。波面画像２８０と波面画像２９０との間でも同様の処理を行い、波面２８１が波面２９１の位置に、波面２８２が波面２９２の位置に、波面２８３が波面２９３の位置に、それぞれ移動したことを検出する。ここで、波面２７３、波面２８２、波面２９２で示される１の波面については、他の波面と比べて移動距離が著しく小さい（伝播速度が著しく遅い）。このような波面は誤検知である可能性が高いので、ノイズとして消去する。これにより、図６（ｅ）の波面画像３００に示すように、波面３０１、３０２が検出できる。
【００３１】
これらの動作により、時刻ごとの波面画像データが生成できる。波面抽出部１５は、生成した複数の波面画像データを波面画像記憶部１９に出力する。このとき、生成した複数の波面の対応情報も波面画像記憶部１９へ出力してもよい。波面の対応情報とは、同一の波面が各波面画像のどの波面に対応するかを示した情報であり、例えば、波面２７２が波面２８３の位置に移動したことが検出された場合、波面２８３と波面２７２とが同一の波面であるという情報である。
【００３２】
図３に戻って説明を続ける。次に、波面画像加工部１６は、波面の表示態様を加工する（ステップＳ６０）。具体的には、波面画像加工部１６は、同一時刻の波面の表示態様が同一となるように、各波面画像データにおける波面の表示態様を変化させる。波面の表示態様の例を図７（ａ）〜（ｄ）に示す。なお、図７（ａ）〜（ｄ）では、時刻の異なる複数の波面画像を重畳合成して表示している。例えば、図７（ａ）の表示例に示すように、時刻ごとに、波面を示す線の太さを変えてもよい。または、例えば、図７（ｂ）の表示例に示すように、時刻ごとに、波面を示す線の輝度を変えてもよい。なお、時刻ごとに、波面を示す線の色を変えてもよいし、色と輝度の両方を変えるとしてもよい。このとき、例えば、ある時刻の波面画像について、当該時刻の波面を赤色で表示するとともに、他の時刻の波面（全ての時刻でもよいし、当該時刻より前の時刻全てでもよいし、当該時刻より遡って特定数のみでもよい）を黒色で当該波面画像内に重畳合成するとしてもよい。なお、赤色、黒色は一例であり、ほかの色でもよいことは言うまでもない。このようにすることで、波面の移動とともに、波面の移動の軌跡を確認しやすくなる。または、例えば、図７（ｃ）の表示例に示すように、時刻ごとに、実線、破線、一点破線など線の種類を変えてもよい。なお、これらの表示態様を組み合わせ、例えば、波面を示す線の太さと種類を変える、としてもよい。また、波面の表示態様としては、波面を線で示す場合に限らず、例えば、図７（ｄ）の表示例に示すように、１つの波面を同種の複数の矢印で示し、時刻ごとに、矢印の表示態様を変えてもよい。ここで、矢印の表示態様を変える手法としては、例えば、図７（ｄ）の表示例に示すように矢印の太さを変化させる、図７（ｅ）の表示例に示すように矢印の輝度や色を変化させる、矢印の長さを変化させる、などがある。波面を矢印で示す場合、波面の進行方向を矢印の向き、波面の伝播速度を矢印の長さとすれば、波面の進行方向と伝播速度がより明確となる。なお、波面の進行方向と伝播速度は、時間的に連続する２つの波面画像における波面の動きから取得できる。すなわち、時間的に連続する２つの波面画像における、波面の移動方向が波面の進行方向、移動距離が波面の伝播速度を示している。
【００３３】
最後に、制御部１１が波面を動画として表示する（ステップＳ７０）。具体的には、波面画像加工部１６が加工した複数の波面画像データを画面表示用の画像データとなるよう幾何変換し、幾何変換後の波面画像データを表示部３に連続表示する。
≪波面画像の例≫
以下に、硬い内包物を含む組織をせん断波が通過する場合の波面の例を示す。
【００３４】
図８に、せん断波の通過と反射が発生している波面画像の例を示す。図８（ａ）〜（ｈ）は、図面左側から右側に向かって伝播してきたせん断波の伝播の様子を、その順で示している。なお、図８（ａ）〜（ｈ）では、周辺組織より硬い内包物３２０を示しているが、波面画像には波面３１０または波面３１１、３１２しか描画されない。
図８（ａ）〜（ｃ）に示すように、組織の硬さがほぼ一定である場合には、せん断波の伝播速度がほぼ一定であるため、波面３１０は形状がほぼ変化せずに伝播する。一方、せん断波が硬い内包物３２０の表面に到達すると、図８（ｄ）に示すように、内包物３２０の内部に透過波が進行するとともに、内包物３２０の表面で反射波が発生する。図８（ｄ）に透過波の波面３１１と反射波の波面３１２とを示す。このとき、透過波については、せん断波の伝播速度は組織が硬いほど速いため、波面３１１は、内包物３２０の内部では波面が早く進行し、図８（ｅ）〜図８（ｇ）に示すように、波面の形状が変化する。図では内包物３２０は円形であり、波面が内包物３２０の内部を進行する距離が異なるため、波面は進行方向に凸な形状となっている。なお、透過波の波面３１１の進行方向への突出は、せん断波が再び内包物３２０の外に進行すればせん断波の伝播速度は一定となるため、図８（ｈ）に示すように、透過波の波面３１１の形状は変化しない。一方、反射波は透過波の波面３１１が内包物３２０の表面に接触する際に発生するため、反射波の波面３１２は、内包物３２０の表面のうち左半分の領域から、放射状に広がるように形成される。したがって、図８（ｅ）〜図８（ｇ）に示すように、内包物３２０から離れるほど、反射波の波面３１２の長さは長くなる。
【００３５】
≪まとめ≫
上記構成により、プッシュパルスに起因するせん断波の伝播の様子を確認することが可能となる。そのため、せん断波が反射等をすることなくプッシュパルスの印加場所から伝播したのか、内包物等によって反射が生じたのかが確認可能となり、検査者が組織の硬さの測定結果が信頼できるものか否かを判断できる。また、せん断波が反射して組織の硬さが取得できなかった場合には、組織の硬さの測定結果が信頼できないことが分かると同時に、せん断波がどのように伝播したのか知ることができる。そのため、せん断波の反射の原因となっている内包物の位置とおよその形状とを知ることができ、せん断波の反射の影響を減らすようにプッシュパルスの印加場所を変更することができる。したがって、検査者は、高い確率で結果が得られるような対処を行ったうえで再検査を行うことができる。
【００３６】
≪変形例１≫
実施の形態１では、変位画像のそれぞれについて変位の大きい変位領域を抽出して波面を検出してから、波面画像の時間変化に基づいてフィルタリングを行う場合について説明した。これに対し、本変形例は、組織内の各位置における変位の時間変化に基づいて、波面を検出することを特徴とする。
【００３７】
≪構成≫
図９は、本変形例に係る波面抽出部２０の機能ブロック図である。図９に示すように、波面抽出部２０は、変位ピーク検出部３６、時間フィルタリング部３４、空間フィルタリング部３３を備える。変位ピーク検出部３６は、変位画像における座標ごとの変位の時間変化に基づいて座標ごとに変位量がピークとなる時刻を算出し、これに基づいて各時刻における変位が大きい変位領域を算出し、変位領域それぞれを波面として抽出することで波面画像データを生成する。時間フィルタリング部３４は、波面画像データに対して時間フィルタリングを行う。空間フィルタリング部３３は、時間フィルタリング後の波面画像に対して空間フィルタリングを行う。
【００３８】
なお、本変形例に係る超音波診断装置の構成は、波面抽出部１５に替えて上述した波面抽出部２０を備える以外、超音波診断装置１と同様であるので、超音波診断装置１と共通の構成については説明を省略する。
≪動作≫
図１０は、本変形例に係る波面抽出の動作を示すフローチャートである。なお、本変形例１に係る超音波診断装置の動作は、以下に説明する波面抽出の動作が実施の形態１に係る超音波診断装置１と異なり、それ以外の動作は超音波診断装置１と同様であるので説明を省略する。
【００３９】
まず、各画像の変位を検出する（ステップＳ５１）。具体的な動作は、実施の形態１で説明したので省略する。変位検出部１４は、１の断層画像に係る各画素の変位を当該画素の座標と対応付けることで変位画像を生成し、生成した変位画像を変位量記憶部１８に出力する。
次に、各画素について時間領域毎に変位を相関処理し、変位のピークを検出する（ステップＳ１５１）。変位ピーク検出部３６は、変位画像を変位量記憶部１８から取得する。以下、図１１の模式図を用いて説明する。図１１（ａ）は、変位画像の一例を示している。図３、図５と同じく、図中の“○”は着目領域における被検体内の組織の一部を示しており、プッシュパルスを印加する前の位置は破線の交点である。また、ｘ軸は超音波探触子２における素子の並ぶ方向、ｙ軸は、被検体の深さ方向である。変位ピーク検出部３６は、座標ごとに変位量δが極値となる時刻を検出する。ここで、極値とは、極大値、および、極小値（変位量δが負の値で、絶対値が極大）のことである。本変形例においては、変位量δが極大値（変位量δが正の値）となる時刻のみを検出することとする。図１１（ａ）のグラフ４１０は、座標（ｘ1，ｙ1）の点４０１における変位量δを時間ｔの関数として示したものである。なお、ここでは、変位量δを１次元ベクトルとして扱うため、変位量δはｙ軸方向の変位のみを示す値であり、かつ、深い方向（ｙ座標が大きくなる方向）への変位は正の値、浅い方向（ｙ座標が小さくなる方向）への変位は負の値としている。波面抽出部は、例えば、断層画像信号が取得された時間を複数の区間４２１、４２２、４２３に分ける。そして、波面抽出部は、区間４２１における変位量δ、区間４２２における変位量δ、区間４２３における変位量δについて相関処理を行うことで、変位量δがピークとなる時刻ｔ1、ｔ2を検出する。ここで、断層画像信号が取得された時間を複数の区間に区切っているのは、着目点を反射波が通過するなど、波面が着目点を複数回通過することがあり得るからである。なお、本変形例では、変位量δが極大となる時刻のみを検出するとしたが、さらに、相関処理において相関値が負の値で絶対値が極大となる時刻、すなわち、変位量δが負の値で、その絶対値が極大となる時刻を検出してもよい。そのようにすることで、位相が反転している波面をさらに検出することができる。
【００４０】
次に、変位が大きい画素の集合を波面として抽出する（ステップＳ１５２）。具体的には、変位ピーク検出部３６は、各時刻について、変位量δがピークとなっている点の集合を変位領域として取得し、変位領域のそれぞれを波面として抽出する。すなわち、図１１（ｂ）の波面画像４３０に示すように、時刻ｔ1において、座標（ｘ1，ｙ1）の点を含む変位領域４３１や、変位領域４３２、４３３、４３４、４３５、４３６が、それぞれ波面として抽出される。なお、上述したように、変位量δが負の値で絶対値が極大値となっている画素の集合を基に、上記検出した波面とは位相が反転している位相反転波面をさらに抽出してもよい。変位ピーク検出部３６は、抽出した波面を時刻ごとの波面画像データとして時間フィルタリング部３４に出力する。
【００４１】
次に、時間フィルタリング部３４は、複数の波面画像データに対して時間フィルタリングを行い、伝播していない波面を除去する（ステップＳ１５３）。具体的には、実施の形態１に係るステップＳ５６と同様、時間的に連続する２以上の波面画像において、各波面位置の時間変化を検出し、速度が異常である波面をノイズとして除去する。波面抽出部は、例えば、時刻ｔ＝ｔ1の波面画像４４０と時刻ｔ＝ｔ1＋Δｔの波面画像４５０との間、時刻ｔ＝ｔ1＋Δｔの波面画像４５０と時刻ｔ＝ｔ1＋２Δｔの波面画像４６０との間のそれぞれにおいて、波面位置の時間変化の検出を行う。波面位置の時間変化の検出はステップＳ５６と同様、例えば、波面画像４４０内の各波面に対し、波面画像４５０において、当該波面と同じ位置を中心に、波面と垂直な向きにΔｔの間にせん断波が移動しうる領域で相関処理することで行う。例えば、波面画像４４０と波面画像４５０との間で、波面４４１が波面４５１の位置に、波面４４２が波面４５２の位置に、波面４４３が波面４５３の位置に、それぞれ移動したと検出する。ここで、波面４４４の移動先が検出されていないのは、以下の理由による。波面Ａと波面Ｂとを含む波面画像Ｘと、波面Ｃを含む波面画像Ｙとの相関処理において、波面画像Ｙで、波面Ａを中心とする領域と波面Ｂを中心とする領域とが重なっており、その重複領域に波面Ｃが存在する場合、以下のように移動先を判定する。すなわち、波面Ａと波面Ｃとの相関値と、波面Ｂと波面Ｃとの相関値とをそれぞれ算出し、（１）いずれか一方の相関値のみが閾値を超えていれば、波面Ａと波面Ｂのうち、波面Ｃとの相関値が閾値を超えている側を波面Ｃの移動元とする、（２）両方の相関値が閾値を超えていれば、波面Ａと波面Ｂのうち、波面Ｃとの相関値が高い方を波面Ｃの移動元とする、（３）いずれの相関値も閾値を下回っていれば、波面Ｃの移動元は存在しない、と判定する。波面４４２、４４３、４５２、４５３は略同じ長さであるのに対して波面４４４は他の波面より極端に短いため、波面４４４と波面４５３の相関値より波面４４３と波面４５３の相関値の方が高い。そのため、波面画像４５０において波面４４４と同じ位置を中心とした領域内に、波面４４４の移動先となる波面が存在しない。移動先が存在しない波面４４４は、ノイズとして消去する。同様に、波面画像４５０と波面画像４６０との間で、波面４５１が波面４６１の位置に、波面４５２が波面４６２の位置に、波面４５３が波面４６３の位置に、それぞれ移動したと検出する。時間フィルタリング部３４は、時間フィルタリング後の波面画像データを、空間フィルタリング部３３に出力する。なお、時間フィルタリング部３４は、実施の形態１で述べたように、さらに、波面の対応情報を空間フィルタリング部３３に出力する、としてもよい。
【００４２】
次に、空間フィルタリング部３３は、時間フィルタリング後の波面画像データに対して空間フィルタリングを行い、長さが短い波面を除去する（ステップＳ１５４）。具体的には、実施の形態１に係るステップＳ５４と同様、各波面画像において、波面の長さを検出し、所定の長さより短い波面、例えば、全ての波面の長さの平均値の１／２よりも長さが短い波面を、ノイズとして削除する。具体的には、図９（ｄ）の波面画像４７０に示すように、波面４７１〜４７３の長さの平均値を算出し、それよりも短い波面４７１を、ノイズとして消去する。これにより、誤検出された波面を消去できる。なお、空間フィルタリング部３３は、波面の対応情報を受け取っている場合、波面の対応情報から、消去した波面と同一である他の波面画像内の波面をまとめて消去してもよい。図９（ｅ）に示すように、空間フィルタリング後の波面画像４８０には、波面４７２、４７３のみが残る。
【００４３】
波面抽出部２０は、空間フィルタリング後の複数の波面画像データを波面画像記憶部１９に出力する。なお、波面抽出部２０は、波面の対応情報をさらに出力する、としてもよい。
≪まとめ≫
本変形例では、組織内の各位置における変位の時間変化に基づいて、波面を検出する場合について説明した。各座標について時系列に変位のピークを検出するため、時間変化のない変位を検出せず、基準断層画像信号内のノイズや、基準断層画像信号と断層画像信号との位置ずれなどによる、時間変化しないノイズを変位として誤検出することを防ぐことができる。また、各座標について時系列に変位のピークを複数回検出するため、１つの場所を波面が複数回通過する場合、特に、透過波と反射波が通過するような場合にも、波面を抽出することができる。また、各座標について時系列に変位のピークを検出する際、相関値が負の値で絶対値が極大となる時刻を検出することで、せん断波の位相が反転した場合に、その位相が反転した波面を検出することができ、固定端反射となるようなせん断波の反射に対しても、反射の波面を検出することができる。
【００４４】
≪変形例２≫
実施の形態１および変形例１では、透過波と反射波とで表示態様を変更しなかった。本変形例では、透過波と反射波とで表示態様を異ならせる、または、一方のみを表示することに特徴がある。
≪波面分類の動作≫
本変形例では、波面を透過波と反射波に分類する。ここで、透過波とは、プッシュパルスの印加によって発生し、プッシュパルスの印加場所から遠ざかるように伝播している波面のことである。一方、反射波は、被検体内において透過波が反射され、透過波とは異なる伝播をしている波面のことである。
【００４５】
図１２は、波面分類の動作を示すフローチャートである。まず、取得時間が最も早い波面画像中の各波面について、波面の位置と進行方向を推定する（ステップＳ２０１）。
次に、波面が、プッシュパルス印加位置の近傍から離れるように進行しているか否かを判定する（ステップＳ２０２）。具体的には、波面が、以下の条件の全てを満たしているか否かを判定する：（１）波面の位置が、プッシュパルスの印加位置の近傍である、（２）波面画像の取得された時刻が、プッシュパルスの印加時刻の直後、または、一定期間内である、（３）波面が、プッシュパルスの印加位置から遠ざかるように進行している。すべてを満たしている場合は、プッシュパルスによって発生したせん断波であると推定できるので、その波面は透過波であると判別し（ステップＳ２０３）、ステップＳ２０４に進む。一方、１つでも条件を満たしていない場合は、次の判定であるステップＳ２０５に進む。
【００４６】
ステップＳ２０５では、波面の位置を直前に透過波または反射波（以下、これらを「参照波面」と呼ぶ）が通過しているか否かを判断する。参照波面とは、具体的には、着目している波面画像が取得された時刻の直前、または、当該時刻の前であって時差が一定範囲内である時刻に取得された波面画像中の波面であって、分類の対象となっている波面の位置に存在する波面のことである。参照波面がある場合には、参照波面と同一の波面、または、参照波面に起因して発生した波面であることが推定できるので、次の判定であるステップＳ２０６に進む。一方、参照波面がない場合は、その波面の種類が特定できないので、当該波面の種類を未分類のままにしてステップＳ２０４に進む。
【００４７】
ステップＳ２０６では、波面と参照波面の進行方向がほぼ一致しているか逆方向であるかを判定する。波面と参照波面の進行方向がほぼ一致している場合、波面と参照波面が同一の波面であると推定できるので、波面を参照波面と同一の種類に分類する。すなわち、参照波面が透過波であれば当該波面は透過波であると判別し（ステップＳ２０７）、参照波面が反射波であれば当該波面は反射波であると判別する（ステップＳ２０８）。一方、波面と参照波面の進行方向が逆である場合、当該波面は、発生個所で参照波面が反射した波面であると推定できるので、当該波面は反射波であると判別する（ステップＳ２０８）。これにより、当該波面の種類が判別したので、ステップＳ２０４に進む。
【００４８】
ステップＳ２０４では、次の時刻に取得された波面画像データが存在するか否かを判別し、ある場合にはステップＳ２０１に戻って次の時刻に取得された波面画像データについて、当該波面画像中の波面の分類を行い、なければ処理を終了する。
このようにして波面を透過波と反射波に分類することで、透過波と反射波とのうち一方だけを表示したり、区別して表示したりすることが可能となる。例えば、図１３（ａ）に示す波面画像では、透過波５２０と反射波５３０とを区別なく両方表示しているが、図１３（ｂ）に示す波面画像のように、透過波５２０だけを表示してもよいし、図１３（ｃ）に示す波面画像のように、反射波５３０だけを表示してもよい。なお、図１３（ａ）〜（ｅ）において、反射波５３０の生成要因は内包物５４０であるが、内包物５４０は説明のために示しており、内包物５４０は波面画像上では直接確認できない。また、図１３（ｄ）に示すように、透過波５２０と反射波５３０とで輝度や色を変更するなど、表示態様を変更してもよい。なお、表示態様は輝度や色のみならず、線種や線の太さであってもよいし、透過波を線で、反射波を矢印で表示するとしてもよい。また、図１３（ｅ）に示すように、さらに、反射境界５５０を示してもよい。なお、反射境界５５０は、上述したステップＳ２０６において、波面と参照波面の進行方向が逆である場合に、当該波面の発生個所を記憶し、全ての反射波の発生個所を累積することで検出することができる。
【００４９】
≪まとめ≫
本変形例では、透過波と反射波とで表示態様を異ならせる、または、一方のみを表示する場合について説明した。このようにすることで、せん断波の反射が生じている場合に、反射波のみを表示して反射境界や反射波によるせん断波の速度が取得できなかった場所などを確認することができ、また、透過波と反射波の干渉により透過波の伝播状況が確認しづらくなることを防ぐことができる。したがって、反射波の発生位置および透過波の伝播状況をより的確に検査者に提示することが可能となる。
【００５０】
≪変形例３≫
実施の形態１および変形例１、２では、プッシュパルスを１度だけ１か所に送信する場合について説明した。これに対し本変形例では、プッシュパルスを２度、すなわち複数の箇所へ送信する場合について説明する。
プッシュパルスを２か所に送信する場合、ステップＳ２０においてプッシュパルスを１回の送信において２か所に送信してもよいし、１度目のプッシュパルスの送信の後、ステップＳ３０の超音波信号の取得を行い、その後に２度目のプッシュパルスを送信し、さらに、ステップＳ３０の超音波信号の取得を行うとしてもよい。
【００５１】
波面画像については、プッシュパルスが１か所では図１４（ａ）に示す波面画像となるのに対し、プッシュパルスが２か所では図１４（ｂ）に示す波面画像となる。ここで、図１４（ａ）および（ｂ）において、位置５０１に印加されたプッシュパルスに起因するせん断波の波面は５００である。また、図１４（ｂ）において、位置５０２に印加されたプッシュパルスに起因するせん断波の波面は５１０である。
【００５２】
なお、図１４（ｃ）の波面画像に示すように、起因するプッシュパルスごとに、波面の表示態様、例えば、輝度や色を変化させてもよい。ここで、波面をプッシュパルスと対応付ける方法について説明する。動作は、変形例２と同様であるが、さらに、いずれのプッシュパルスに起因する透過波または反射波であるかの判定を行う。具体的には、ステップＳ２０２において、第１のプッシュパルスの印加位置近傍で発生した透過波を第１のプッシュパルスに起因する透過波、第２のプッシュパルスの印加位置近傍で発生した透過波を第２のプッシュパルスに起因する透過波と判定する。さらに、ステップＳ２０７〜Ｓ２０８において、当該波面が参照波面と同じプッシュパルスに起因する透過波または反射波であると判定する。これにより、全ての波面について、いずれのプッシュパルスに起因する透過波または反射波であるかの判定ができる。なお、変形例２と同様、例えば、１つのプッシュパルスに起因する波面のみを表示するとしてもよい。
【００５３】
≪まとめ≫
本変形例では、プッシュパルスを複数回送信する場合について説明した。このようにすることで、プッシュパルスを同時に複数個所に印加した場合であっても、せん断波の波面を確認することができる。また、第１のプッシュパルスに起因する波面と第２のプッシュパルスに起因する波面との両方が通過する領域において、反射等が発生していないか確認することができ、また、透過波同士の干渉により伝播状況が確認しづらくなることを防ぐことができる。したがって、せん断波の伝播状況をより的確に検査者に提示することが可能となる。
【００５４】
≪実施の形態２≫
実施の形態１では、波面画像のみを提供してせん断波の伝播状態を示すだけの場合について説明した。これに対し、本実施の形態では、波面画像を提供するとともに、断層画像、組織の硬さを併せて提供する場合について説明する。
≪構成≫
図１５は、実施の形態２に係る超音波診断装置４のブロック図である。なお、図１と同じ構成については同じ符号を用い、説明を省略する。
【００５５】
弾性推定部２１は、波面画像記憶部１９に記憶されている、波面画像データを基に着目領域内の各座標における弾性率を推定し、座標ごとの弾性率を色で示した弾性画像として弾性画像記憶部２２に出力する。具体的には、時間的に連続する波面画像データから、波面の移動距離を算出する。波面の移動距離は波面が通過した場所におけるせん断波の速度を示しており、弾性率（例えば、ヤング率）はせん断波の速度から算出可能であるため、波面が通過した場所の弾性率が算出できる。
【００５６】
弾性画像記憶部２２は、弾性画像を記憶する。弾性画像記憶部２２は、例えば、ＲＡＭ、フラッシュメモリ、ハードディスクなどの記憶媒体で実現される。
≪動作≫
図１６に、せん断波の速度を算出する動作の模式図を示す。図１６（ａ）は、ある時刻ｔにおける波面画像と、時刻ｔ＋Δｔにおける波面画像を１つの波面画像５６０として合成したものである。ここで、時刻ｔにおける波面５６１と、時刻ｔ＋Δｔにおける波面５６２とが同一の波面であるとする。弾性推定部２１は、波面５６１と波面５６２との相関処理を行い、波面５６１上の座標（ｘ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）2＋（ｙt+Δt−ｙt）2｝／Δｔとなる。
【００５７】
弾性推定部２１は、全ての波面に対して上述の処理を行い、波面が通過した全座標についてせん断波の速度を取得し、せん断波の速度を基に、弾性率を算出する。弾性推定部２１は、弾性率に基づいて、色情報をマッピングした弾性画像を生成する。例えば、図１６（ｂ）のように、弾性率が一定値以上の座標は赤、弾性率が一定値未満の座標は緑というように色分けした弾性画像５７０を生成する。分類は二値化に限られず、所定の段階で分類及び色分けを行ってよい。図１６（ｂ）において、領域５８０は弾性率が一定値以上の領域であり、内包物５９０に対応する。なお、内包物５９０は弾性画像上には直接現れない。弾性推定部２１は、生成した弾性画像を弾性画像記憶部２２に出力する。
【００５８】
制御部１１は、例えば、断層画像記憶部１７に記憶されている断層画像信号に、弾性画像記憶部２２に保持されている弾性画像と、波面画像を重畳して表示する。なお、制御部は、断層画像信号、弾性画像、波面画像の全てに画面表示用の画像データとなるような幾何変換と、断層画像信号に対してはさらに、包絡線検波・対数圧縮による輝度変換を行う。図１７（ａ）は、断層画像信号を変換した断層画像上に、弾性画像と、波面画像とを重畳した例を示している。領域６０１は、弾性画像と、波面画像とが重畳されている着目領域である。領域６０３は、弾性画像上において、弾性率が一定値以上の領域である。また、領域６０１には波面６０２がさらに重畳されている。上述したように、弾性率が高い、すなわち硬い組織ほどせん断波の伝播速度が速いので、図１７（ａ）に示すように、波面６０２は領域６０３において進行方向に突出した形状となる。なお、制御部１１は、例えば、図１７（ｂ）に示すように、断層画像信号を変換した断層画像上に、波面画像を重畳してもよい。
【００５９】
≪まとめ≫
本実施の形態では、組織の弾性率の測定結果とせん断波の伝播の様子を重畳して表示する場合について説明した。このようにすることで、組織の弾性率の測定結果について、当該位置のせん断波の伝播の様子を確認することができるので、弾性率の測定結果における信頼性をより正確に評価できる。
【００６０】
≪実施の形態に係るその他の変形例≫
（１）各実施の形態および変形例では、変位検出部１４において、プッシュパルス送信前の断層画像信号を基準断層画像信号として、プッシュパルス送信後の断層画像信号との間でパターンマッチングや相関処理等による変位検出を行うとしたが、本発明は必ずしもこの場合に限られない。例えば、ステップＳ２０における断層画像信号の取得を行わず、基準断層画像信号として、着目領域を撮像した、予め断層画像記憶部１７に記憶されている断層画像信号、または、着目領域の代表的な断層画像信号を用いてもよい。また、プッシュパルス送信後の断層画像信号のうち、一定時間離れた異なる時間における２つの断層画像信号の片方を基準断層画像信号として用いてもよい。
【００６１】
（２）各実施の形態および変形例では、超音波信号取得部１３が音響線信号を生成し、断層画像記憶部１７が音響線信号を記憶し、変位検出部１４は１画像分の音響線信号である断層画像信号を用いて変位の検出を行う場合について説明したが、本発明は必ずしもこの場合に限られない。例えば、超音波信号取得部１３は断層画像信号に対して包絡線検波、対数圧縮による輝度変換や画面表示用の画像データとなるよう幾何変換を行うことで、断層画像信号を断層画像に変換し、断層画像記憶部１７は断層画像を記憶するとしてもよい。このようにすることで、制御部１１は波面画像の幾何変換を行う必要がなく、また、断層画像記憶部１７が記憶している断層画像に対し、対数圧縮による輝度変換や幾何変換をする必要がなく、そのまま表示部３に表示することが可能となる。
【００６２】
（３）各実施の形態、変形例２、３では、変位領域抽出部３１は、変位画像の全域に対して変位領域を抽出するとしたが、本発明は必ずしもこの場合に限られない。例えば、変位領域抽出部３１は、プッシュパルスの印加直後の変位画像に対して、プッシュパルスの印加位置の近傍のみ変位領域を抽出するとしてもよい。このようにすることで、演算量を低減させることができるとともに、明らかにせん断波が存在しない領域を抽出の対象としないことで、ノイズを波面として誤検出することを防ぐことができる。なお、変位領域の抽出範囲の限定はこの場合に限られず、例えば、あるｙ座標について全くせん断波が検出されない場合には、以降の時刻について当該ｙ座標を変位量位置の抽出範囲から除外してもよい。
【００６３】
また、同様に、時間フィルタリング部３４は、プッシュパルスの印加時刻を基に時間フィルタリング処理を行うとしてもよい。例えば、時間フィルタリング部３４は、プッシュパルスの印加直後の波面画像に対して、プッシュパルスの印加位置の近傍にない波面をノイズとして削除してもよい。このようにすることで、明らかにせん断波が存在しえない領域に存在するノイズを波面として誤検出することを防ぐことができる。
【００６４】
（４）各実施の形態および各変形例では、空間フィルタリング部３３と時間フィルタリング部３４とが個別にフィルタリング処理を行うとしたが、本発明は必ずしもこの場合に限られない。例えば、空間フィルタリング処理と時間フィルタリング処理は、同時に行ってもよい。
（５）各実施の形態および変形例２、３では、変位領域抽出部３１が、ｙ座標ごとに変位量δを座標ｘの関数として、および、ｘ座標ごとに変位量δを座標ｙの関数として、動的閾値を用いることで変位量δが大きい領域を抽出する場合について説明したが、本発明は必ずしもこの場合に限られない。例えば、変位領域抽出部３１は、ｙ座標ごとに変位量δを座標ｘの関数として変位量δが大きい領域を抽出する処理を全ｙ座標についてのみ行ってもよいし、逆に、ｘ座標ごとに変位量δを座標ｙの関数として変位量δが大きい領域を抽出する処理を全ｘ座標についてのみ行ってもよい。または、例えば、変位量δを座標ｘ、ｙの関数として、２次元的に変位量δが大きい領域を抽出するとしてもよい。
【００６５】
（６）各実施の形態および変形例２、３では、細線化部３２が変位領域の細線化処理を行う場合について説明したが、本発明は必ずしもこの場合に限られない。例えば、変位領域抽出部３１が変位領域における変位の大きい画素を線で繋ぎ、これを波面としてもよい。このようにすることで、細線化の処理が不要となる。
（７）変形例１では、変位ピーク検出部３６は、相関処理により、変位画像から変位量がピークとなる時刻を検出する場合について説明したが、本発明は必ずしもこの場合に限られない。例えば、相関処理による相関値に対してしきい値を設け、しきい値を上回る相関値が検出された時間帯を検出するとしてもよい。このようにすることで、変位量がピークとなる時刻をピンポイントで検出するだけではなく、その周辺時刻を検出して波面を正確に検出することが可能となる。また、逆に、相関処理による相関値に対して下限しきい値を設け、下限しきい値を下回る相関値が検出された場合、ピークとして検出しないとしてもよい。このようにすることで、波面が通過していない時刻における揺動などを波面と誤検出することを防ぐことができる。
【００６６】
（８）各実施の形態および各変形例では、波面画像加工部１６が波面の表示態様を加工する場合について説明したが、本発明は必ずしもこの場合に限られない。例えば、波面画像加工部１６は、波面画像記憶部１９に記憶されている波面画像データの表示形態を変更しないとしてもよい。
（９）各実施の形態および変形例では、ステップＳ２０において、超音波信号取得部１３が、送信超音波が送信フォーカス点に同時に届くように送信ビームフォーミングを行い、送信イベントごとに素子受信信号に整相加算を行い、音響線信号を生成する場合について説明したが、本発明は必ずしもこの場合に限られない。例えば、超音波信号取得部１３は、超音波探触子２を構成する各振動素子から同時に超音波を送出することで、平面波状の送信超音波を送出してもよい。この場合、受信フォーカス点を送信フォーカス点の近傍や同一線上に限定する必要がないため、１つの断層画像を構成する音響線信号を、１送信イベントごとに１つではなく１画像分まとめて取得することができる。そのため、断層画像信号のフレームレートを向上させることができる。
【００６７】
（１０）各実施の形態および変形例では、超音波診断装置は表示部３と接続される構成であるとしたが、本発明は必ずしもこの場合に限られない。例えば、超音波診断装置１は表示部３を内蔵しているとしてもよいし、または、超音波診断装置１は表示部３と接続されず、波面画像加工部１６が生成した波面画像データを記憶媒体に保存、あるいは、ネットワークを通じて他の装置に出力するとしてもよい。
【００６８】
また、同様に、超音波診断装置は超音波探触子２を内蔵しているとしてもよいし、あるいは、超音波探触子２が超音波信号取得部１３を備え、超音波信号取得部１３を有さない超音波診断装置が超音波探触子２から音響線信号を取得するとしてもよい。または、超音波診断装置は、音響線データあるいは断層画像を、記録媒体から読み取り、もしくは、ネットワークを通じて他の装置から取得するとしてもよい。同様に、超音波診断装置は波面画像を重畳するための断層画像あるいは弾性画像を、記録媒体から取得もしくはネットワークを通じて他の装置から取得するとしてもよい。
【００６９】
（１１）各実施の形態および各変形例では、制御部１１が波面を動画として表示するとしたが、本発明は必ずしもこの場合に限られない。例えば、制御部１１は、複数の波面画像データのいくつかを画面表示用の画像データとなるよう幾何変換し、幾何変換後の波面画像データを表示部３に表示してもよい。または、例えば、制御部１１は、複数の波面画像データの１つを画面表示用の画像データとなるよう幾何変換し、幾何変換後の波面画像データを表示部３に表示してもよい。あるいは、例えば、制御部１１は、複数の波面画像データを１つの画像データとして重畳合成した後、画面表示用の画像データとなるよう幾何変換し、幾何変換後の波面画像データを表示部３に表示してもよい。
【００７０】
（１２）各実施の形態および各変形例に係る超音波診断装置は、その構成要素の全部又は一部を、１チップ又は複数チップの集積回路で実現してもよいし、コンピュータのプログラムで実現してもよいし、その他どのような形態で実施してもよい。例えば、波面抽出部を１チップで実現してもよいし、変位領域抽出部のみを１チップで実現し、細線化部等を別のチップで実現してもよい。
【００７１】
集積回路で実現する場合、典型的には、ＬＳＩ（Ｌａｒｇｅ  Ｓｃａｌｅ  Ｉｎｔｅｇｒａｔｉｏｎ）として実現される。ここでは、ＬＳＩとしたが、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路、又は汎用プロセッサで実現してもよい。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサを利用してもよい。
【００７２】
さらには、半導体技術の進歩、又は派生する別技術により、ＬＳＩに置き換わる集積回路化の技術が登場すれば、当然その技術を用いて機能ブロックの集積化を行ってもよい。
また、各実施の形態および各変形例に係る超音波診断装置は、記憶媒体に書き込まれたプログラムと、プログラムを読み込んで実行するコンピュータとで実現されてもよい。記憶媒体は、メモリカード、ＣＤ−ＲＯＭなどいかなる記録媒体であってもよい。また、本発明に係る超音波診断装置は、ネットワークを経由してダウンロードされるプログラムと、プログラムをネットワークからダウンロードして実行するコンピュータとで実現されてもよい。
【００７３】
（１３）以上で説明した実施の形態は、いずれも本発明の好ましい一具体例を示すものである。実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、工程、工程の順序などは一例であり、本発明を限定する主旨ではない。また、実施の形態における構成要素のうち、本発明の最上位概念を示す独立請求項に記載されていない工程については、より好ましい形態を構成する任意の構成要素として説明される。
【００７４】
また、発明の理解の容易のため、上記各実施の形態で挙げた各図の構成要素の縮尺は実際のものと異なる場合がある。また本発明は上記各実施の形態の記載によって限定されるものではなく、本発明の要旨を逸脱しない範囲において適宜変更可能である。
さらに、超音波診断装置においては基板上に回路部品、リード線等の部材も存在するが、電気的配線、電気回路について当該技術分野における通常の知識に基づいて様々な態様を実施可能であり、本発明の説明として直接的には無関係のため、説明を省略している。尚、上記示した各図は模式図であり、必ずしも厳密に図示したものではない。
【００７５】
≪補足≫
（１）実施の形態に係る超音波診断装置は、超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して、特定部位にある組織を物理的に押圧し、押圧により特定部位の組織が振動源となってせん断波が隣接組織へと伝播する過程を、プッシュパルスに続いて超音波の被検体への送信と、被検体からの反射超音波の受信とを繰り返すことによって検出する超音波診断装置であって、前記プッシュパルスに続く超音波の被検体への送信と被検体からの反射超音波の受信により時系列に複数の受信信号を得て、前記複数の受信信号が取得された各時刻における前記被検体内の組織の変位をそれぞれ検出する変位検出部と、前記複数の受信信号が取得された各時刻における組織の変位に基づいて、当該時刻における前記せん断波の波面を抽出する波面抽出部と、前記複数の受信信号が取得された各時刻における波面を波面画像として生成する波面画像加工部とを備えることを特徴とする。
【００７６】
また、実施の形態に係る超音波診断装置の制御方法は、超音波探触子を用い、被検体内の特定部位に超音波を集中させるプッシュパルスを送信して、特定部位にある組織を物理的に押圧し、押圧により特定部位の組織が振動源となってせん断波が隣接組織へと伝播する過程を、プッシュパルスに続いて超音波の被検体への送信と、被検体からの反射超音波の受信とを繰り返すことによって検出する超音波診断装置の制御方法であって、前記プッシュパルスに続く超音波の被検体への送信と被検体からの反射超音波の受信により時系列に複数の受信信号を得て、前記複数の受信信号が取得された各時刻における前記被検体内の組織の変位をそれぞれ検出し、前記複数の受信信号が取得された各時刻における組織の変位に基づいて、当該時刻における前記せん断波の波面を抽出し、前記複数の受信信号が取得された各時刻における波面を波面画像として生成することを備えることを特徴とする
本開示によれば、せん断波が伝播している間の各時刻の波面を画像として検査者が確認することができるので、検査者は、せん断波の波面が反射や屈折をすることなく伝播したか否かを確認することができる。さらに、検査者は、せん断波の反射や屈折が発生した場合には、どの場所でそのような事態が起きたのかを、波面の形状や、時刻ごとの波面の位置およびそれが示す波面の伝播速度や方向に基づいて容易に推測できる。したがって、検査者は、せん断波の屈折や反射により被検体内の組織の硬さを取得できなかった場合に、せん断波の伝播の様子をもとに、超音波探触子の位置を変える、被検者の姿勢を変える、プッシュパルスを印加する位置を変える、などの適切な対処を行って再検査を行うことができる。
【００７７】
（２）また、実施の形態に係る上記（１）の超音波診断装置は、前記波面画像加工部は、前記波面の時間変化を表示する画像をさらに生成する、としてもよい。
このようにすることで、時系列に波面の位置と伝播方向を確認することができ、検査者はより的確にせん断波の伝播の様子を確認できる。
（３）また、実施の形態に係る上記（２）の超音波診断装置は、前記波面の時間変化を表示する画像は動画である、としてもよい。
【００７８】
このようにすることで、検査者はより直感的にせん断波の伝播の様子を確認できる。
（４）また、実施の形態に係る上記（１）〜（２）の超音波診断装置は、前記波面画像加工部は、前記波面画像を、前記被検体内の組織の各領域の弾性を示す弾性画像に重畳して出力する、としてもよい。
このようにすることで、弾性画像と波面の伝播の様子を同時に確認することができ、検査者は、弾性画像の各領域の信頼性を具体的に知ることができる。
【００７９】
（５）また、実施の形態に係る上記（１）〜（２）の超音波診断装置は、前記波面画像加工部は、前記波面画像を、前記被検体内の組織に対応した断層画像に重畳して出力する、としてもよい。
このようにすることで、断層画像と波面の伝播の様子を同時に確認することができ、検査者は、断層画像とせん断波の伝播の様子との双方を加味して被検体の組織の状態をより正確に推測することができる。
【００８０】
（６）また、実施の形態に係る上記（１）〜（５）の超音波診断装置は、前記波面抽出部は、前記受信信号のそれぞれにおいて、変位量が閾値よりも大きい変位領域を抽出し、抽出した各変位領域に細線化処理を施すことで、当該受信信号が取得された時刻の波面を抽出する、としてもよい。
このようにすることで、変位の大きさに基づいて、波面を抽出することができる。
【００８１】
（７）また、実施の形態に係る上記（１）〜（５）の超音波診断装置は、前記波面抽出部は、前記複数の受信信号が示す被検体内の各位置における変位が極値となる時刻を相関処理によって検出し、検出された被検体内の各位置における変位が極値となる時刻に基づいて前記受信信号のそれぞれにおける変位領域を抽出することで、当該受信信号が取得された時刻の波面を抽出する、としてもよい。
【００８２】
このようにすることで、変位の時系列変化に基づいて、波面の通過回数や波面の進行向きにかかわらず各波面を検出することができる。
（８）また、実施の形態に係る上記（１）〜（７）の超音波診断装置は、前記波面画像加工部は、前記波面抽出部が抽出した波面の表示態様を加工し、１の受信信号が取得された異なる時刻における波面について、波面の表示態様が前記時刻ごとに互いに異ならせるように波面の表示態様を加工する、としてもよい。
【００８３】
このようにすることで、波面の区別がしやすくなり、より的確にせん断波の伝播の様子を確認することができる。
（９）また、実施の形態に係る上記（８）の超音波診断装置は、前記波面画像加工部は、時刻の異なる複数の波面に対し、線の太さ、色、輝度、および、種類のうち少なくとも１つを変化させることで表示態様を異ならせる、としてもよい。
【００８４】
（１０）また、実施の形態に係る上記（８）の超音波診断装置は、前記波面画像加工部は、波面のそれぞれを、伝播方向と伝播速度とをそれぞれ向きと長さで表現した矢印で表現し、時刻の異なる複数の波面に対し、矢印の太さ、色、および、輝度のうち少なくとも１つを変化させることで表示態様を異ならせる、としてもよい。
これら（９）〜（１０）のようにすることで、時刻の異なる波面をより明確に区別することが可能となり、より的確にせん断波の伝播の様子を確認することができる。
【００８５】
（１１）また、実施の形態に係る上記（１）〜（１０）の超音波診断装置は、前記波面画像加工部は、前記波面抽出部が抽出した波面のうち、進行方向が同一であるせん断波に係る波面のみを表示する、としてもよい。
（１２）また、実施の形態に係る上記（１）〜（１０）の超音波診断装置は、前記波面画像加工部は、前記せん断波の進行方向ごとに、波面の表示態様を異ならせる、としてもよい。
【００８６】
これら（１１）〜（１２）のようにすることで、せん断波の反射が発生しているか否かを容易に確認することができ、また、透過波と反射波が干渉していても、一方のみを確認することができる。
（１３）また、実施の形態に係る上記（１）〜（１２）の超音波診断装置は、前記波面抽出部は、さらに、前記せん断波の進行方向が変化している界面を抽出し、前記波面画像加工部は、さらに、前記界面を前記波面に重畳する、としてもよい。
【００８７】
このようにすることで、せん断波の反射の原因がどこに存在するのかを容易に確認することが可能となる。
（１４）また、実施の形態に係る上記（１）〜（１３）の超音波診断装置は、前記せん断波励起部は、前記プッシュパルスを第１のプッシュパルスとしたとき、さらに第２のプッシュパルスを前記超音波探触子に送信させ、前記波面画像加工部は、前記第１のプッシュパルスに起因するせん断波の波面と、前記第２のプッシュパルスに起因するせん断波の波面とで、表示態様を異ならせる、としてもよい。
【００８８】
このようにすることで、プッシュパルスの送信が２度以上あった場合に、せん断波同士の干渉があっても、各せん断波の伝播の様子を容易に確認できる。
（１５）また、実施の形態に係る上記（１）〜（１４）の超音波診断装置は、前記波面抽出部が検出した時刻の異なる同一の波面の移動量を、当該時刻におけるせん断波の伝播速度とすることで、前記被検体内の組織のうち、前記波面が通過した各領域の弾性をそれぞれ推定する弾性推定部をさらに備える、としてもよい。
【００８９】
このようにすることで、波面の検出と組織の弾性の測定とをまとめて実施することができ、せん断波の伝播の様子から、組織の弾性に対する信頼性を容易に確認できるとともに組織の弾性に対する信頼性が低い原因を確認することができる。
【産業上の利用可能性】
【００９０】
本開示に係る超音波診断装置は、超音波を用いた組織の硬さの測定において、せん断波の伝播の様子の確認に有用である。そのため、せん断波の伝播に関する情報を組織の硬さの再測定にフィードバックすることが可能となり、医療診断機器等において高い利用可能性を持つ。
【符号の説明】
【００９１】
１、４  超音波診断装置
２  超音波探触子
３  表示部
１１  制御部
１２  せん断波励起部
１３  超音波信号取得部
１４  変位検出部
１５  波面抽出部
１６  波面画像加工部
１７  断層画像記憶部
１８  変位量記憶部
１９  波面画像記憶部
２０  波面抽出部
２１  弾性推定部
２２  弾性画像記憶部
３１  変位領域抽出部
３２  細線化部
３３  空間フィルタリング部
３４  時間フィルタリング部
３６  変位ピーク検出部
【図１】
0006384340.tif 000002
【図２】
0006384340.tif 000003
【図３】
0006384340.tif 000004
【図４】
0006384340.tif 000005
【図５】
0006384340.tif 000006
【図６】
0006384340.tif 000007
【図７】
0006384340.tif 000008
【図８】
0006384340.tif 000009
【図９】
0006384340.tif 000010
【図１０】
0006384340.tif 000011
【図１１】
0006384340.tif 000012
【図１２】
0006384340.tif 000013
【図１３】
0006384340.tif 000014
【図１４】
0006384340.tif 000015
【図１５】
0006384340.tif 000016
【図１６】
0006384340.tif 000017
【図１７】
0006384340.tif 000018