(57)【要約】放射圧エラストグラフィのせん断波速度の計測において、計測時間が長引くことによる計測精度と再現性の劣化を低減し、高い診断能をもつ超音波画像を取得可能な技術を提供する。放射圧エラストグラフィにおいて、トラックパルスの照射によるエコー信号から、せん断波を検出しながら、計測領域の動き(変動)に関する情報を抽出し、計測結果の信頼度を示す信頼性情報としてユーザに提供する。また、抽出した情報から、変動の要因を特定し、ユーザに提示する。さらに、複数回の計測結果を加算平均する際、信頼性情報で重みづけを行う。
2015029651.tif 000002
【特許請求の範囲】
【請求項１】
被検体に集束したバースト超音波を送信することにより当該バースト超音波の焦点位置に生成されたせん断波を、複数のせん断波検出用パルスを繰り返し送信することにより得たエコー信号群を用いて検出するせん断波検出部と、
前記せん断波の伝搬速度であるせん断波速度を計算する速度計算部と、
前記せん断波の伝搬領域を含む計測領域の変動を評価し、評価結果を、前記せん断波速度の信頼度を示す信頼性情報として得る変動評価部と、
前記信頼性情報をユーザに提示する提示部と、を備えること
を特徴とする超音波診断装置。
【請求項２】
請求項１記載の超音波診断装置であって、
前記変動評価部は、前記計測領域内であって、前記せん断波の伝搬領域の深度とは異なる深度領域における前記変動に基づいて、前記信頼性情報を得ること
を特徴とする超音波診断装置。
【請求項３】
請求項１記載の超音波診断装置であって、
前記信頼性情報は、前記変動の大きさを示す指標であること
を特徴とする超音波診断装置。
【請求項４】
請求項１記載の超音波診断装置であって、
前記変動評価部は、さらに、前記変動の要因を特定し、
前記提示部は、前記特定された要因をさらにユーザに提示すること
を特徴とする超音波診断装置。
【請求項５】
請求項１記載の超音波診断装置であって、
前記バースト超音波の送信と前記複数のせん断波検出用パルスの繰り返しの送信と当該送信によるエコー信号の受信とからなる計測を、予め定めたパルスシーケンスに従って実行するシーケンス制御部と、
複数の前記せん断波速度の加算平均を計算する加算平均部と、をさらに備え、
前記シーケンス制御部は、前記計測を繰り返し、
前記せん断波検出部は、前記計測毎に前記せん断波を検出し、
前記速度計算部は、前記せん断波を検出する毎に前記せん断波速度を計算し、
前記加算平均部は、前記せん断波を検出する毎に計算される複数の前記せん断波速度の加算平均を計算し、
前記提示部は、前記加算平均結果を前記信頼性情報とともにユーザに提示すること
を特徴とする超音波診断装置。
【請求項６】
請求項５記載の超音波診断装置であって、
前記加算平均部は、前記加算平均時に前記信頼性情報を用いて重みづけを行うこと
を特徴とする超音波診断装置。
【請求項７】
請求項５記載の超音波診断装置であって、
前記提示部は、前記計測毎の前記信頼性情報および前記せん断波速度を散布図にし、さらに提示すること
を特徴とする超音波診断装置。
【請求項８】
請求項７記載の超音波診断装置であって、
ユーザから前記散布図上のプロット結果を介して指示を受け付ける受付部をさらに備え、
前記加算平均部は、前記指示に応じて、前記加算平均を再計算すること
を特徴とする超音波診断装置。
【請求項９】
請求項８記載の超音波診断装置であって、
前記受付部は、除外する前記せん断波速度の選択を受け付け、
前記加算平均部は、前記選択されたせん断波速度以外のせん断波速度を用いて、前記加算平均を再計算すること
を特徴とする超音波診断装置。
【請求項１０】
請求項８記載の超音波診断装置であって、
前記受付部は、前記せん断波速度を前記信頼性情報に応じて複数のグループに分割する指示を受け付け、
前記加算平均部は、グループごとに前記せん断波速度の加算平均を再計算すること
を特徴とする超音波診断装置。
【請求項１１】
請求項５記載の超音波診断装置であって、
ユーザから再計測の指示を受け付ける受付部をさらに備え、
前記シーケンス制御部は、前記指示に応じて、前記計測を実行すること
を特徴とする超音波診断装置。
【請求項１２】
請求項２記載の超音波診断装置であって、
前記変動評価部は、前記せん断波が伝搬する深度領域を、当該せん断波の発生位置と前記せん断波の振幅とにより特定すること
を特徴とする超音波診断装置。
【請求項１３】
請求項２記載の超音波診断装置であって、
前記変動評価部は、前記せん断波が伝搬する深度領域を、前記エコー信号群から得たデータに対して時間方向に相関演算を行って得た相関係数を用いて特定すること
を特徴とする超音波診断装置。
【請求項１４】
請求項４記載の超音波診断装置であって
前記変動評価部は、前記エコー信号群から得たデータに対し、時間方向に相関演算を行って得た相関係数の変化のパターンにより前記要因を特定すること
を特徴とする超音波診断装置。
【請求項１５】
被検体に集束したバースト超音波を送信することにより当該バースト超音波の焦点に生成されたせん断波を、複数のせん断波検出用パルスを繰り返し照射することにより得たエコー信号群を用いて検出し、
前記せん断波の伝搬速度であるせん断波速度を計算し、
前記せん断波の伝搬領域を含む計測領域の変動を評価し、評価結果を前記せん断波速度の信頼度を示す信頼性情報として得、
前記信頼性情報をユーザに提示すること
を特徴とする弾性評価方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波を用いて被験体内部の情報を非侵襲に取得する超音波イメージング技術に関し、特に、組織の硬さをイメージングするエラストグラフィ技術に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波を体外から照射し、体内で反射され返ってきた信号を、経過時間と信号強度とによって画像化する医用画像装置である。超音波は、音響インピーダンスの異なる界面でスネルの法則に従って反射される性質を有するため、生体内の組織によって微妙に異なる音響インピーダンスの違いを可視化することで、組織の構造を描写する。【０００３】
超音波診断装置を用い、組織の構造ではなく、組織の硬さを画像化する、エラストグラフィ技術がある。組織の硬さは病変と密接な関係にあり、診断に重要な情報を齎す。このエラストグラフィの手法のひとつに、せん断波を生じさせ、せん断波の伝播によって生じる変位からせん断速度を測定し、組織の硬さを得る放射圧エラストグラフィがある。組織のポアソン比を０．５として換算し、縦波速度は横波速度に比べて充分に大きいとすると、硬さの指標となるヤング率Ｅは、簡易的に以下の式(１)で表される。
【数１】
2015029651.tif 000003
ここで、ρは密度，Ｖｓは、せん断波速度である。この式(１)を用い、せん断波速度から硬さの絶対値を得る。【０００４】
せん断波は、集束超音波を一点に照射し、組織に放射圧を印加することで生じさせる。このとき、印加するパルスを、放射圧生成用パルス(プッシュパルス)と呼ぶ。プッシュパルスにより生成されたせん断波の変位は、せん断波検出用パルス(トラックパルス)で検出する。【０００５】
放射圧エラストグラフィでは、硬さの絶対値を算出するため、精度よくせん断波による変位を計測し、再現性の高いせん断速度を算出する必要がある。再現性を向上させるため、１回の測定で、計測領域内の複数個所で、複数回せん断波速度を計測し、得られた計測値の平均を測定値として提示する手法がある(例えば、特許文献１参照)。特許文献１に開示の手法では、測定値を、測定値からはずれた値の大きさによって評価し、評価結果とともに提示する。本手法によれば、加算平均時のはずれ値の影響が少なく抑えられるため測定精度も向上すると考えられる。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】米国特許出願公開２０１０／００１６７１８号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、上述のように、放射圧エラストグラフィでは、プッシュパルスおよびトラックパルスという２種類の超音波パルスを照射する。特許文献１に記載の手法では、これらを複数回繰り返し照射することとなるため、計測に時間を要する。そのため、被検体の呼気・心拍などに由来する体動や、ユーザの手ぶれなどによる撮像面のずれが生じ、計測の精度および再現性の劣化を引き起こす。【０００８】
具体的には、まず上述の動きによって、計測範囲にずれが生じ、測定したい部位とは異なる部位を測定することが考えられる。また、軽微なずれであっても、検出されたせん断波の本来の時相と測定した時相とがずれ、検出したせん断波速度が、本来の伝搬速度とずれることが考えられる。さらに、特に体動などによって面ずれが生じない場合においても、例えば、心拍による肝臓の“圧迫度合い(歪)”が、時相によって異なると考えられ、この圧迫度合いはせん断波速度に影響するため、計測精度劣化につながり得る。【０００９】
本発明は、上記事情に鑑みてなされたもので、放射圧エラストグラフィのせん断波速度の計測において、計測時間が長引くことによる計測精度と再現性との劣化を低減し、高い診断能をもつ超音波画像を取得可能な技術を提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明は、放射圧エラストグラフィにおいて、トラックパルスの照射によるエコー信号から、せん断波を検出しながら、計測領域の動き(変動)に関する情報を抽出し、計測結果の信頼度を示す信頼性情報としてユーザに提供する。また、抽出した情報から、変動の要因を特定し、ユーザに提示する。さらに、複数回の計測結果を加算平均する際、信頼性情報で重みづけを行う。
【発明の効果】
【００１１】
本発明により、放射圧エラストグラフィにおいて、計測時間が長引くことによる計測精度と再現性との劣化が低減され、高い診断能を持つ超音波画像を取得できる。
【図面の簡単な説明】
【００１２】
【図１】本発明の実施形態の超音波診断装置のブロック図である。
【図２】(ａ)は、本発明の実施形態のＢモード画像例を説明するための説明図であり、(ｂ)は、本発明の実施形態の計測領域（図２（ａ）の領域２２０）の拡大図である。
【図３】本発明の実施形態の計測領域内における、相関係数の深さ方向の変化を説明するための説明図である。
【図４】(ａ)は、本発明の実施形態の撮像領域のＢモード画像例を説明するための説明図である。(ｂ)〜(ｄ)は、非せん断波変動による相関係数の変化を説明するための説明図である。
【図５】本発明の実施形態の表示画面例を説明するための説明図である。
【図６】(ａ)〜(ｄ)は、本発明の実施形態の表示画面例を説明するための説明図である。
【図７】本発明の実施形態の撮像処理のフローチャートである。
【図８】本発明の実施形態の変形例の超音波診断装置のブロック図である。
【図９】(ａ)〜(ｃ)は、本発明の実施形態の変形例のユーザから受け付ける指示を説明するための説明図である。
【図１０】本発明の実施形態の変形例の表示後の処理のフローチャートである。
【発明を実施するための形態】
【００１３】
以下、本発明の実施形態の一例を、図面を用いて説明する。なお、各実施形態を説明するための全図において、特に断らない限り、同一機能を有するものは、同一名称および同一符号を付し、その繰り返しの説明は省略する。また、本明細書に記載した、せん断波速度とは、せん断波の伝搬速度を指す。本発明を適用する各実施形態では、せん断波速度に基づいて、例えば歪、ヤング率、粘性率、体積弾性率などの組織性状に関する情報に関する評価が可能となる。
【００１４】
まず、本実施形態の超音波診断装置１００を説明する。図１は、本実施形態の超音波診断装置１００のブロック図である。
【００１５】
本実施形態の超音波診断装置１００では、被検体の計測領域に放射圧を照射(送信)し、せん断波を発生させるための集束バースト超音波(以後、プッシュパルスと称す)と、プッシュパルスの送信により生じたせん断波の伝播を検出するパルス超音波(以後、トラックパルスと称す)とを送信する計測を行い、計測領域の組織の性質として、せん断波の伝搬速度を得る放射圧エラストグラフィ法を行う。また、信頼度および再現性を高めるため、上記計測を繰り返し、得られた結果を加算平均する。
【００１６】
このとき、本実施形態の超音波診断装置１００では、トラックパルスのエコー信号から、計測領域の動き(変動)に関する情報を抽出し、放射圧エラストグラフィで得た情報の信頼性に関する指針としてユーザに提供する。また、抽出した変動に関する情報から、変動の要因を特定し、ユーザに提示する。さらに、変動に関する情報を、上記加算平均時の重みづけにも用いる。
【００１７】
このため、本実施形態の超音波診断装置１００は、本図に示すように、送受ビームフォーマ１１０と、シーケンス制御部１２０と、送信条件設定部１３０と、画像生成部１４０と、弾性評価部１５０と、を備える。また、超音波診断装置１００には、探触子１６０、入力装置１７０、表示装置１８０が接続される。
【００１８】
＜送信ビームフォーマ＞
送受ビームフォーマ１１０は、シーケンス制御部１２０からの指示に従って、送信ビームを探触子１６０に送信するとともに、探触子１６０によって受信されたエコー信号を受信する。
【００１９】
具体的には、探触子１６０の各素子から送信する超音波パルスの電気信号を生成する。生成された電気信号は、送信ビームフォーマが備えるＤ／Ａ変換機によりアナログ信号に変換後、探触子１６０に送られ被検体に照射される。被検体中を伝搬する過程で音響インピーダンスが異なる界面から反射されてきた信号は、受信エコー信号として探触子１６０で受信され、送信時とは逆の過程を経てデジタル信号に変換され、整相加算などの加算処理が行われ、減衰補正等の処理がなされた後、複素のＲＦデータに変換される。
【００２０】
＜シーケンス制御部＞
シーケンス制御部１２０は、送信条件設定部１３０を介して設定された撮像条件に従って、超音波パルスを送信するタイミングおよびエコー信号を受信するタイミング、送信する超音波パルスの特性などをパルスシーケンスとして決定する。そして、決定したパルスシーケンスに従って、送受ビームフォーマ１１０を制御し、計測を実行する。本実施形態では、放射圧エラストグラフィを実行する。このため、本実施形態のシーケンス制御部１２０は、プッシュパルスの送信と、複数のトラックパルスの繰り返しの送信と、トラックパルスによるエコー信号の受信と、を行う計測を実行するようパルスシーケンスを生成する。
【００２１】
＜送信条件設定部＞
送信条件設定部１３０は、ユーザから受け付けたせん断波を発生させる位置(以後、計測領域と呼ぶ)に応じて、当該計測領域にプッシュパルスの送信条件、および、当該領域で発生するせん断波を検出するためのトラックパルスの送信条件を設定する。設定する送信条件は、集束位置、送信角度、バースト長、電圧、周波数、および、送信開口などの音圧パラメータを含む。
【００２２】
図２(ａ)および図２(ｂ)は、プッシュパルスおよびトラックパルス送信の概念を説明するための図である。図２(ａ)は、Ｂモード画像２１０の例であり、図２(ｂ)は、Ｂモード画像２１０内の計測領域２２０の拡大図である。矢印２３４は、深度方向である。
【００２３】
計測領域２２０内のプッシュパルスの焦点２２２で生成された、放射圧によるせん断波２２１は、組織中を伝播していく。トラックパルスは、このせん断波を検出するために送信される。このため、プッシュパルス一回に対し、トラックパルスは最短でもせん断波の伝播時間の間送信され続ける。
【００２４】
プッシュパルスの送信条件は、指定された計測領域２２０内の、所望の位置２２２にプッシュパルスを送信するよう設定され、トラックパルスの送信条件は、プッシュパルスにより発生するせん断波２２１を、そのエコー信号により計測可能なよう、設定される。また、トラックパルスについては、１回の計測における送信数、繰り返し回数、繰り返し毎の複数のトラックパルスの送信位置なども送信条件として設定される。
【００２５】
＜画像生成部＞
画像生成部１４０は、送受ビームフォーマ１１０で得た複素のＲＦデータを、シーケンス制御部１２０の制御に従って受信し、断層像を生成する。画像生成部１４０は、１つのエコー信号(ビーム)から得たＲＦデータの輝度値を、受信した時間に応じて深さ方向にプロットする。これを、探触子１６０の長軸方向に複数ビーム分並べることで２次元の情報を蓄積し、蓄積された情報から断層像を生成する。生成された断層像は、表示装置１８０に表示される。
【００２６】
例えば、エコー信号の強度を輝度として画像化するＢモード撮像では、探触子１６０長軸方向のビーム本数は撮像フレームレートに影響する。リアルタイム性を保持するため、通常、１枚のＢモード画像を取得するには、数十―数百ビームが用いられる。
【００２７】
＜探触子＞
探触子１６０は、上述のせん断波計測用のシーケンスが送受できる探触子１６０であれば良く、好適には、リニア、コンベックス、もしくは、セクタ形状の１Ｄアレイ探触子、または、３次元撮像用の１．５次元、もしくは、２次元アレイ探触子などが用いられる。
【００２８】
＜弾性評価部＞
弾性評価部１５０は、計測領域２２０の組織の硬さの情報を得る。本実施形態では、プッシュパルスの送信により生じたせん断波を検出し、その速度(せん断波速度)を得ることにより、組織の硬さを示す情報を得る。せん断波速度は、せん断波の伝搬によって生じる変位から算出する。さらに、本実施形態の弾性評価部１５０は、得られたせん断波速度の信頼度を示す情報(信頼性情報)を算出し、ユーザに提示する。
【００２９】
これを実現するため、本実施形態の弾性評価部１５０は、図１に示すように、相関演算部１５１と、せん断波検出部１５２と、速度計算部１５３と、変動評価部１５４と、加算平均部１５５と、提示部１５６と、を備える。
【００３０】
＜相関演算部＞
相関演算部１５１は、受信したエコー信号から得たＲＦデータに対し、時間方向に相関演算を行う。本実施形態では、ＲＦデータが複素ＲＦデータであるため、複素相互相関演算が実施される。複素相互相関演算は、時間的に隣接するＲＦデータ間で行ってもよいし、基準とするＲＦデータを決め、基準とするＲＦデータとの間で行ってもよい。
【００３１】
＜せん断波検出部＞
せん断波検出部１５２は、被検体１０１に集束したバースト超音波(プッシュパルス)を送信することによりプッシュパルスの焦点に生成されたせん断波を、複数のせん断波検出用パルス(トラックパルス)を繰り返し送信することにより得た受信エコー信号を用いて検出する。本実施形態では、相関演算部１５１における複素相互相関結果から、せん断波のピークを検出し、検出位置と検出時間とを得る。なお、本実施形態では、ピークの検出の前に、複素相互相関結果に対し、最適なフィルタ処理を行う。
【００３２】
＜速度計算部＞
速度計算部１５３は、せん断波の伝搬速度であるせん断波速度を計算する。本実施形態では、せん断波速度は、せん断波のピークが検出された時間、位置、および、せん断波生成位置により計算される。具体的には、プッシュパルスの焦点と、せん断波のピークを観測したトラックパルスの送信位置とから計算される。
【００３３】
＜変動評価部＞
変動評価部１５４は、せん断波の伝搬領域を含む計測領域２２０の変動を評価し、評価結果をせん断波速度の信頼度を示す信頼性情報として得る。評価する対象の変動は、放射圧エラストグラフィの計測精度や再現性に影響を与える変動である。本実施形態では、まず、計測領域２２０内で、この変動を検出する領域を特定し、当該領域の所定の位置（評価位置）の組織の変動を評価する。
【００３４】
まず、変動を検出する領域の特定手法について説明する。
【００３５】
放射圧エラストグラフィの計測精度や再現性に影響を与える変動は、２種類考えられる。一方は、ユーザの保持する探触子１６０が動くことに由来する撮像面全体の面ずれで、他方は、心拍、呼吸など、被検体の周期的体動に由来する計測部位のずれである。以後、前者を術者由来の変動、後者を体動由来の変動と呼ぶ。
【００３６】
ここで注意すべきことは、放射圧エラストグラフィにおける計測対象のせん断波そのものも、組織の微弱な動きを計測することにより得ていることである。従って、計測対象のせん断波による動きと、計測精度や再現性に影響を与える変動とを区別し、後者のみ検出する必要がある。以後、前者の動きをせん断波変動と呼び、後者の変動を、非せん断波変動と呼ぶ。本実施形態の変動評価部１５４は、せん断波変動の影響を受けない領域を特定し、当該領域内の所定の位置を評価位置とし、当該評価位置の非せん断波変動を評価する。
【００３７】
ここで、せん断波変動の影響を受けない位置について説明する。前述のように、図２(ａ)は、超音波画像(Ｂモード画像)２１０の概念図であり、図２(ｂ)は、図２(ａ)内の計測領域２２０の拡大図である。
【００３８】
図２(ａ)に示すように、超音波視野内では、組織が層構造で表示される。ここでは、層２１１、層２１２、層２１３の三層構造を例示する。また、図２(ｂ)に示すように、せん断波２２１は、放射圧が生成された部位(せん断波生成位置)２２２から横方向に伝搬する。すなわち、深さ方向(図中下向き)については、一定範囲(図中、領域ｂ２３２；以下、せん断波伝搬領域と呼ぶ。)内でのみ伝搬する。従って、せん断波伝搬領域２３２以外、例えば、図中、領域ａ２３１または領域ｃ２３３などが、せん断波の動きの影響を受けない位置といえる。これらの領域を、非伝搬領域２３１、２３３と呼ぶ。
【００３９】
本実施形態の変動評価部１５４は、計測領域２２０内であって、せん断波が伝搬する深度領域(せん断波伝搬領域２３２)の深度とは異なる深度領域(非伝搬領域２３１、２３３)内の所定の位置（評価位置）における変動に基づいて、信頼性情報を算出する。評価位置は、非伝搬領域２３１、２３３内で、できる限りせん断波生成位置に近い位置とする。
【００４０】
なお、変動評価部１５４は、せん断波検出部１５２の検出結果を用いて、せん断波伝搬領域２３２を決定し、非伝搬領域２３１、２３３を特定する。せん断波伝搬領域２３２は、当該せん断波の発生位置と前記せん断波の振幅とにより特定される。せん断波２２１の発生位置は、プッシュパルスにより放射圧が生成される位置である。この位置は、プッシュパルス生成に用いた素子数から求められる送信開口幅およびフォーカスの深度により決定される、焦点深度である。また、せん断波２２１の振幅は、せん断波検出部１５２が検出したせん断波のピークの位置と焦点深度との距離で特定される。
【００４１】
また、せん断波伝搬領域２３２は、せん断波検出部１５２の検出結果を用いず、独立して特定してもよい。例えば、相関演算部１５１による複素相互相関演算結果を用いる。すなわち、相関係数を用いて特定してもよい。
【００４２】
一般に、動き(変動)のある位置は、せん断波変動や非せん断波変動など、動きの種類にかかわらず相関係数が低下する。この中で、上述のように、せん断波２２１は、深さ方向に限定された位置、例えば、図２(ｂ)では、せん断波伝搬領域２３２にのみ生成される。このため、相関係数の深さ方向の変化をみると、図３に模式的に示すように、せん断波による相関係数の変化２４１において、低下は局所的である。一方、非せん断波変動、即ち探触子１６０の面ずれや体動による相関係数の変化２４２、２４３は、図２(ｂ)の領域によらず、一定である。なお、実際に得られる複素相互相関演算結果は、これらの相関係数の変化２４１、２４２、２４３の合成である。
【００４３】
これを利用し、相関係数が局所的に低下する領域を検出し、検出した領域をせん断波伝搬領域２３２とする。局所的に低下する領域は、例えば、微分演算などで検出する。せん断波伝搬領域２３２決定後、非伝搬領域２３１、２３３の特定手法、評価位置の決定手法は、上記と同様とする。
【００４４】
本実施形態の変動評価部１５４は、上述の手法で特定した位置（評価位置）の変動の激しさ(変動の大きさ)を指標化した変動インデックスｋｍを、信頼性情報として算出する。一般に、相互相関演算の相関係数は、動きが激しいほど低下する。これを利用し、本実施形態では、当該評価位置の、例えば、正規化した相互相関係数の平均値を変動インデックスｋｍとする。この場合、変動インデックスｋｍの値は、変動が激しいほど小さくなる。なお、ここで用いる相関演算は、せん断波検出時と共通した相関演算であってもよいし、異なる相関演算であってもよい。
【００４５】
なお、上述のように、本実施形態では、計測領域２２０内の複数個所で、複数回計測し、複数のせん断波速度を得る。また、各々の計測において、変動インデックスｋｍを算出する。このとき、変動評価部１５４は、複数回の計測において得られた各変動インデックスｋｍの、平均値に対するばらつき、例えば、標準偏差、をさらに算出し、信頼性情報としてもよい。また、各計測で得たせん断波速度の標準偏差を算出し、信頼性情報としてもよい。
【００４６】
次に、変動評価部１５４が、せん断波検出部１５２が検出した変動の要因を特定する手法について説明する。ここでは、要因として、体動由来であるか、探触子１６０保持の手技等の術者由来であるかを特定する。特定は、相関係数の変化のパターンを識別することにより行う。この特定手法を、図４(ａ)から図４(ｄ)を用いて説明する。図４(ａ)は、撮像領域のＢモード画像３１０である。図４(ｂ)から図４(ｄ)は、撮像領域内の、非せん断波変動によるトラックパルスの複素相関係数の変化のパターンを模式的に示したものである。
【００４７】
たとえば肝臓の撮像を行った場合、図４(ａ)に示すように、Ｂモード画像３１０上では、肝臓３１２の手前に皮膚・筋肉・脂肪などの表在組織３１１が、奥には横隔膜で隔てられた消化管など他組織３１３が観測される。矢印３１４は、深さ方向である。
【００４８】
術者由来の場合、その要因には、探触子１６０の横滑りと煽り方向のずれとの２種がある。探触子１６０の横滑りは、探触子１６０の設置位置自体がずれることにより発生する。また、煽り方向のずれは、探触子１６０の、設置位置はずれないが、角度がずれることにより発生する。
【００４９】
横滑り(探触子１６０の設置位置のずれ)の場合の相関係数の変化のパターン３４０を図４(ｂ)に示す。本図に示すように、横滑りの場合、計測領域３２０外も含む、撮像面全ての領域で同時に相関係数の低下が生ずる。
【００５０】
探触子１６０が煽り方向にずれた場合の相関係数の変化の３５０を図４(ｃ)に示す。本図に示すように、この場合、計測領域３２０外も含む撮像面全ての領域で相関係数の低下が生ずる。ただし、このとき、相関係数が低下するタイミングは同時ではなく、より深い部位、即ち探触子１６０表面から遠い位置ほど、早く低下が始まる。
【００５１】
対照的に、拍動や呼気などの周期的体動による変動は、組織によって異なる。具体的には、皮膚や脂肪層など表在組織３１１や、横隔膜を隔てて肝臓の奥にある他組織３１３は動かないのに対し、中間部に位置する肝臓３１２などの臓器が特徴的でかつ周期的な動きを示す。このため、この場合の相関係数の変化のパターン３６０は、図４(ｄ)に示すように、計測領域３２０を含む、中間部の臓器(肝臓３１２)の位置でのみ、周期的に変化するものとなる。
【００５２】
本実施形態の変動評価部１５４は、これらの相関係数の変化のパターンを検出し、変動が、術者由来か、あるいは、周期的体動由来かを判別する。なお、例えば、これらの変化の基準とする基準パターン、あるいは、基準パターンを特定可能な情報を、超音波診断装置１００が備える記憶装置に、予め保持しておき、変動評価部１５４は、検出したパターンと基準パターンとを比較することにより、変動の要因の判別を行ってもよい。
【００５３】
なお、判別結果は、提示部１５６によりユーザに提示される。また、このとき、要因として、術者由来と判別された場合、再計測を促すメッセージを表示するよう構成してもよい。また、上述の、変動インデックスｋｍの標準偏差が、予め定めた閾値以上となる場合、計測が不適切と判別し、同じく再計測を促すメッセージを表示するよう構成してもよい。
【００５４】
なお、変動評価部１５４が変動を評価する際に用いる信号は、前述のトラックパルスを送信して得た複素相互相関演算結果が望ましい。しかし、これに限定されない。例えば、Ｂモード撮像で得られたデータを用いてもよい。
【００５５】
＜加算平均部＞
加算平均部１５５は、計測領域２２０内の複数個所で、複数回計測することにより得た、複数のせん断波速度の平均値を算出する。
【００５６】
本実施形態では、このとき、変動インデックスｋｍを、重みとして用いてもよい。すなわち、各計測で得たせん断波速度を、その信頼度（ここでは、変動インデックスｋｍ）に応じて、平均値算出への寄与率を変える。これにより、得られる加算平均速度の信頼度が高まる。
【００５７】
重みづけ平均は、例えば、以下の式(２)により計算する。
【数２】
2015029651.tif 000004
ｎは、計測回数であり(ｎは２以上の整数)、Ｖｓiは、ｉ回目の計測により得たせん断波速度、ｋｍiは、ｉ回目の計測により得た変動インデックス、ｋｍmeanは、ｎ個の変動インデックスｋｍiの平均値、Ｖｓmeanは、重みづけ加算平均により得た加算平均速度である。
【００５８】
＜提示部＞
提示部１５６は、速度計算部１５３が計算した、計測毎のせん断波速度Ｖｓ、加算平均部１５５が算出した加算平均速度Ｖｓmean、および、信頼性情報などをユーザに提示する。本実施形態では、これらの計測結果、算出結果を用い、表示装置１８０に表示する表示情報を生成する。表示情報は、数値であってもよいし、定性的なグラフ、カラーマップ表示であってもよい。
【００５９】
図５に、提示部１５６が表示情報として作成する画面の例を示す。表示画面６００には、計測毎のせん断波速度Ｖｓと変動インデックスｋｍの逆数(１／ｋｍ)との散布図６１０と、参考情報６２０と、が表示される。
【００６０】
散布図６１０は、せん断波速度Ｖｓと変動インデックスｋｍの逆数(１／ｋｍ)とを、それぞれ軸とするグラフに、計測結果をプロットしたものである。
【００６１】
参考情報６２０には、例えば、加算平均部１５５が算出した加算平均速度Ｖmean、変動評価部１５４が算出した、せん断波速度の標準偏差ＳＤ、変動インデックスｋｍの標準偏差ｋｍ(ＳＤ)などが表示される。このとき、変動インデックスｋｍの表示形態は、ｋｍの標準偏差を平均値に対する割合で表示することが望ましいが、他の統計量および絶対値を用いて表示してもよい。
【００６２】
なお、表示画面６００は、再計測の指示を受け付ける受付ボタン６３０を表示するよう構成してもよい。
【００６３】
以下、具体的な表示例を説明する。
【００６４】
図６(ａ)は、撮像領域の変動が少ない場合の表示画面６１１例である。ほぼ動きがなく、良好な撮像の場合、変動インデックスｋｍの逆数は、相対的に小さく、変動インデックスｋｍの逆数が小さい範囲にプロット点がまとまって見られる。
【００６５】
図６(ｂ)に、何らかの動きによって計測のはずれ値がある場合に得られると予想される表示画面６１２の例を示す。プロット点の集団が二つ見られる。数が少ないほうの集団は、変動インデックスｋｍの逆数が大きいことから、動きが大きい場合に得られたせん断波速度であるとユーザに示唆することができる。
【００６６】
図６(ｃ)に、周期的な体動がある場合に得られると予想される表示画面６１３の例を示す。動きの大小によって、プロット点が２つのグループに分かれ、表示される。
【００６７】
図６(ｄ)に、変動インデックスｋｍの逆数の値が大きく、計測が不適切な場合に得られる表示画面６１４の例を示す。提示部１５６は、このような場合、再計測を促すメッセージも併せて表示するよう構成してもよい。
【００６８】
本実施形態の超音波診断装置１００は、ＣＰＵとメモリと記憶装置とを備え、記憶装置に予め保持したプログラムを、ＣＰＵがメモリにロードして実行することにより、シーケンス制御部１２０、画像生成部１４０、送信条件設定部１３０、弾性評価部１５０の機能を実現する。各機能の処理に用いる各種のデータ、処理中に生成される各種のデータは、記憶装置に格納される。なお、弾性評価部１５０の各機能の少なくとも１つは、超音波診断装置１００とデータの送受信が可能な、外部の情報処理装置などに備えられていてもよい。また、上記各部の全部または一部の機能は、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（ｆｉｅｌｄ−ｐｒｏｇｒａｍｍａｂｌｅ  ｇａｔｅ  ａｒｒａｙ）などのハードウェアによって実現してもよい。
【００６９】
＜撮像の流れ＞
次に、本実施形態の超音波診断装置１００による放射圧エラストグラフィ実行時の撮像処理の流れを、図７を用いて説明する。本処理は、ユーザからの指示を契機に開始する。ここでは、Ｎ回、プッシュパルスを送信するものとする。
【００７０】
まず、操作者がＢモード画像上で、せん断波の計測領域を指定する。操作者は、入力装置１７０を介して、この計測領域を指定する。送信条件設定部１３０は、指定された計測領域を受け付け(ステップＳ１００１)、プッシュパルスおよびトラックパルスの送信条件を設定する(ステップＳ１００２)。
【００７１】
プッシュパルスおよびトラックパルスの送信条件が設定されると、シーケンス制御部１２０は、放射圧エラストグラフィ計測を開始する。ここでは、まず、計測回数をカウントするカウンタｎを初期化(ｎ＝１)する(ステップＳ１００３)。そして、設定された条件に従って、プッシュパルスを送信する(ステップＳ１００４)。また、プッシュパルス送信直後から、トラックパルスの送信を開始する(ステップＳ１００５)。
【００７２】
シーケンス制御部１２０は、トラックパルスの送信により得られたエコー信号を複素ＲＦデータに変換し、相関演算部１５１は、当該データに対し、複素相互相関演算を行う(ステップＳ１００６)。複素相互相関演算結果は、せん断波検出部１５２と変動評価部１５４とに入力される。
【００７３】
せん断波検出部１５２は、複素相互相関演算結果から、せん断波のピーク位置とピーク検出時間とを算出し、せん断波を検出する(ステップＳ１００７)。そして、速度計算部１５３は、ピーク位置とピーク検出時間とから、せん断波速度を算出する(ステップＳ１００８)。なお、算出したせん断波速度は、計測回数ｎに対応づけて、記憶装置に保持する。
【００７４】
一方、変動評価部１５４では、複素相互相関演算結果から、信頼性情報を算出する(ステップＳ１００９)。なお、算出した信頼性情報は、計測回数ｎに対応づけて、記憶装置に保持する。
【００７５】
シーケンス制御部１２０は、Ｎ回計測を行ったか否かを判別し(ステップＳ１０１０)、行っていない場合は、カウンタｎを１インクリメントし(ステップＳ１０１１)、ステップＳ１００４へ戻り、処理を繰り返す。
【００７６】
一方、ステップＳ１０１２で、Ｎ回の計測を終えたと判別された場合、加算平均部１５５は、加算平均速度Ｖｓmeanを算出する(ステップＳ１０１２)。加算平均は、変動インデックスを用いて重みづけをしてもよい。このとき、併せて、標準偏差ＳＤを算出してもよい。また、変動評価部１５４が、併せて、信頼性情報の標準偏差値ｋｍ(ＳＤ)を算出してもよい。
【００７７】
提示部１５６は、算出結果を用いて表示画面を生成し、表示装置１８０に表示し(ステップＳ１０１３)、処理を終了する。
【００７８】
以上説明したように、本実施形態の超音波診断装置１００は、被検体に集束したバースト超音波を送信することにより当該バースト超音波の焦点位置に生成されたせん断波を、複数のせん断波検出用パルスを繰り返し送信することにより得たエコー信号群を用いて検出するせん断波検出部１５２と、前記せん断波の伝搬速度であるせん断波速度を計算する速度計算部１５３と、前記せん断波の伝搬領域を含む計測領域の変動を評価し、評価結果を前記せん断波速度の信頼度を示す信頼性情報として得る変動評価部１５４と、前記信頼性情報をユーザに提示する提示部１５６と、を備える。
このとき、前記変動評価部１５４は、前記エコー信号群を用いて前記変動の評価を行ってもよい。
また、前記変動評価部１５４は、さらに、前記変動の要因を特定し、前記提示部１５６は、前記特定された要因もさらにユーザに提示する。
【００７９】
このように本実施形態によれば、放射圧エラストグラフィに用いるトラックパルスから動きに関する情報を抽出し、せん断波を検出しながら計測対象の組織における動きや面ずれを検出する。検出した動きに関する情報をユーザに提供することで、ユーザに計測の信頼性に関する指針を提供できる。また、加算平均時に、動き情報に応じて、せん断波速度計測値の重みづけを行うことで、より信頼性の高い計測値を提供できる。また、動きのパターンによって手ぶれがあるか否かを検出し、手ぶれである場合はその旨ユーザに通知する。
【００８０】
これにより、ユーザは、提供された情報を見て、計測の信頼性を把握できる。そして、適宜、計測をやり直すことができる。このように、本実施形態によれば、より信頼性の高い計測手法を実現できる。このため、本実施形態によれば、放射圧エラストグラフィにおいて、計測精度と再現性の劣化を低減でき、高い診断能をもつ超音波画像(硬さの情報)をユーザに提供できる。
【００８１】
放射圧エラストグラフィによれば、例えば、乳癌などは、周辺組織に比べ優位に硬いため、硬い部分を描写することにより感度良く乳がんを検出することが可能である。また、肝硬変に至る肝炎などは肝臓の硬さが病気の進行度と強く関連しており、肝臓の硬さを計測することで、生検の数を最小限に抑えながら的確な診断と治療経過モニタリングをすることができる。本実施形態は、このような放射圧エラストグラフィの良さをそのまま生かし、新たな計測を追加することなく、上述の効果を得ることができる。
【００８２】
なお、上記実施形態では、提示部１５６が表示画面を生成し、ユーザにせん断波速度とその評価結果とを提示する構成としたが、これに限定されない。例えば、図８に示すように、上記実施形態の構成に加え、さらに、表示画面６００を介して、ユーザからの指示を受け付ける受付部１５７を備えてもよい。この場合、提示部１５６が生成する表示画面６００は、再計測の指示を受け付ける、受付ボタン６３０を備える。
【００８３】
例えば、図９(ａ)に示すように、表示画面６１２が表示された場合、ユーザは、当該表示画面６１２を介して、変動インデックスｋｍの逆数が小さいプロット点を選択する。選択は、例えば、本図に示すように、枠６３１で囲うなどにより行う。受付部１５７は、この選択を受け付け、対応するせん断波速度を特定し、加算平均部１５５に選択されたせん断波速度のみを用いて、加算平均速度を再計算させる。計算結果は、提示部１５６により表示される。
【００８４】
また、図９(ｂ)に示すように、表示画面６１３が表示された場合、ユーザは、例えば、任意に設定した変動インデックスｋｍの逆数で、プロット点群を分割する。分割の指示は、例えば、本図に示すように所定の変動インデックスｋｍの逆数６３２を表示画面６１３の散布図上で指定することにより行う。受付部１５７は、この分割する際に用いる変動インデックスｋｍの逆数の指定を受け付け、当該値で分割されたプロット点群を、それぞれグループ化し、加算平均部１５５にグループごとの加算平均速度を計算させる。計算結果は、提示部１５６により表示される。
【００８５】
また、図９(ｃ)に示すように、表示画面６１４が表示された場合、ユーザは、例えば、受付ボタン６３０を介して、再計測を指示する。受付部１５７は、この指示を受け付けると、シーケンス制御部１２０に再度計測を行うよう指示する。なお、このとき、探触子１６０を持ち直すよう指示を行ってもよい。
【００８６】
この場合の、表示後の処理の流れを、図１０を用いて説明する。
【００８７】
受付部１５７は、表示画面６００を介して、所定のプロット点群の選択を受け付けると(ステップＳ１１０１)、当該プロット点群に対応するせん断波速度データを除外する(ステップＳ１１０２)。そして、加算平均部１５５に、残りのせん断波速度データにより加算平均を再計算させる(ステップＳ１１０３)。このとき、変動評価部１５４が、残りのせん断波速度データの変動インデックスｋｍの標準偏差を算出するよう構成してもよい。そして、提示部１５６は、算出結果から表示画面を生成し、表示装置１８０に表示し(ステップＳ１１０４)、処理を終了する。
【００８８】
一方、受付部１５７は、表示画面６００を介して、プロット点群を分割する、変動インデックスｋｍの逆数の指定を受け付けると(ステップＳ１１０５)、プロット点を、当該変動インデックスｋｍの逆数以上とそれより小さいものとに分けて、それぞれグループ化し(ステップＳ１１０５)、加算平均部１５５に、グループごとに、そのせん断波速度データを用いて、加算平均を再計算させる(ステップＳ１１０６)。このとき、変動評価部１５４が、グループごとのせん断波速度データの変動インデックスｋｍの標準偏差を算出するよう構成してもよい。そして、提示部１５６は、算出結果から、各々のグループの表示画面を生成し、表示装置１８０に表示し(ステップＳ１１０７)、処理を終了する。
【００８９】
また、受付部１５７は、再計算の指示を行う指示ボタンの押下を受け付けた場合(ステップＳ１１０８)、シーケンス制御部１２０に再計測の指示を行う(ステップＳ１１０９)。いずれの指示もない場合、処理を終了する。
【００９０】
このように、受付部１５７を備え、信頼性情報に基づいて、ユーザからの指示を受け付けることにより、信頼性評価結果を計測にフィードバックでき、計測精度、再現性をさらに高めることができる。
【符号の説明】
【００９１】
１００：超音波診断装置、１１０：送受ビームフォーマ、１２０：シーケンス制御部、１３０：送信条件設定部、１４０：画像生成部、１５０：弾性評価部、１５１：相関演算部、１５２：せん断波検出部、１５３：速度計算部、１５４：変動評価部、１５５：加算平均部、１５６：提示部、１５７：受付部、１６０：探触子、１７０：入力装置、１８０：表示装置、２０１：計測領域、２１０：Ｂモード画像、２１１：組織層、２１２：組織層、２１３：組織層、２２０：計測領域、２２１：せん断波、２２２：せん断波発生位置、２３１：非伝搬領域、２３２：せん断波伝搬領域、２３３：非伝搬領域、２３４：矢印、２４１：相関係数の変化、２４２：相関係数の変化、２４３：相関係数の変化、３１０：Ｂモード画像、３１１：表在組織、３１２：肝臓、３１３：他組織、３１４：矢印、３２０：計測領域、３４０：相関関数の変化のパターン、３５０：相関関数の変化のパターン、３６０：相関関数の変化のパターン、６００：表示画面、６１０：散布図、６１１：表示画面、６１２：表示画面、６１３：表示画面、６１４：表示画面、６２０：参考情報、６３０：受付ボタン、６３１：枠、６３２：指定位置
【図１】
2015029651.tif 000005
【図２】
2015029651.tif 000006
【図３】
2015029651.tif 000007
【図４】
2015029651.tif 000008
【図５】
2015029651.tif 000009
【図６】
2015029651.tif 000010
【図７】
2015029651.tif 000011
【図８】
2015029651.tif 000012
【図９】
2015029651.tif 000013
【図１０】
2015029651.tif 000014
【手続補正書】
【提出日】平成28年12月2日(2016.12.2)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
被検体の所定の焦点位置に集束したバースト超音波を送信した後、当該バースト超音波により生成されて、前記被検体を伝搬するせん断波を検出するための複数のせん断波検出用パルスを前記被検体に対して繰り返し送信し、そのエコーを受信することにより、エコー信号群を得る送受ビームフォーマと、
前記エコー信号群から、前記被検体の前記せん断波による変動を検出するせん断波検出部と、
前記せん断波の速度を計算する速度計算部と、
前記エコー信号群から、前記被検体の前記せん断波以外の変動を、前記せん断波による変動とは区別して検出する変動評価部と、
前記変動評価部の検出結果をユーザに提示する提示部と、を備えること
を特徴とする超音波診断装置。
【請求項２】
請求項１記載の超音波診断装置であって、
前記変動評価部は、前記計測領域内であって、前記せん断波の伝搬領域の深度とは異なる深度領域における前記被検体の変動を検出することにより、前記せん断波以外の変動を検出すること
を特徴とする超音波診断装置。
【請求項３】
請求項１記載の超音波診断装置であって、
前記変動評価部は、前記せん断波以外の変動の大きさを示す指標を算出すること
を特徴とする超音波診断装置。
【請求項４】
被検体に集束したバースト超音波を送信することにより当該バースト超音波の焦点位置に生成されたせん断波を、複数のせん断波検出用パルスを繰り返し送信することにより得たエコー信号群を用いて検出するせん断波検出部と、
前記せん断波の伝搬速度であるせん断波速度を計算する速度計算部と、
前記せん断波の伝搬領域を含む計測領域の変動を評価し、評価結果を、前記せん断波速度の信頼度を示す信頼性情報として得る変動評価部と、
前記信頼性情報をユーザに提示する提示部と、を備え
前記変動評価部は、さらに、前記変動の要因を特定し、
前記提示部は、前記特定された要因をさらにユーザに提示すること
を特徴とする超音波診断装置。
【請求項５】
請求項１記載の超音波診断装置であって、
前記バースト超音波の送信と、前記複数のせん断波検出用パルスの繰り返しの送信と、前記エコー信号の受信とからなる計測を、予め定めたパルスシーケンスに従って実行するシーケンス制御部と、
複数の前記せん断波速度の加算平均を計算する加算平均部と、をさらに備え、
前記シーケンス制御部は、前記計測を繰り返し、
前記せん断波検出部は、前記計測毎に前記せん断波を検出し、
前記速度計算部は、前記せん断波を検出する毎に前記せん断波速度を計算し、
前記加算平均部は、前記せん断波を検出する毎に計算される複数の前記せん断波速度の加算平均を計算し、
前記提示部は、前記加算平均結果を前記せん断波以外の変動の検出結果とともにユーザに提示すること
を特徴とする超音波診断装置。
【請求項６】
請求項５記載の超音波診断装置であって、
前記変動評価部は、前記せん断波以外の変動の大きさを示す指標を算出し、
前記加算平均部は、前記加算平均時に前記指標を用いて重みづけを行うこと
を特徴とする超音波診断装置。
【請求項７】
請求項５記載の超音波診断装置であって、
前記変動評価部は、前記せん断波以外の変動の大きさを示す指標を算出し、
前記提示部は、前記計測毎の前記指標および前記せん断波速度を散布図にし、さらに提示すること
を特徴とする超音波診断装置。
【請求項８】
請求項７記載の超音波診断装置であって、
ユーザから前記散布図上のプロット結果を介して指示を受け付ける受付部をさらに備え、
前記加算平均部は、前記指示に応じて、前記加算平均を再計算すること
を特徴とする超音波診断装置。
【請求項９】
被検体に集束したバースト超音波を送信することにより当該バースト超音波の焦点位置に生成されたせん断波を、複数のせん断波検出用パルスを繰り返し送信することにより得たエコー信号群を用いて検出するせん断波検出部と、
前記せん断波の伝搬速度であるせん断波速度を計算する速度計算部と、
前記せん断波の伝搬領域を含む計測領域の変動を評価し、評価結果を、前記せん断波速度の信頼度を示す信頼性情報として得る変動評価部と、
前記信頼性情報をユーザに提示する提示部と、
前記バースト超音波の送信と前記複数のせん断波検出用パルスの繰り返しの送信と当該送信によるエコー信号の受信とからなる計測を、予め定めたパルスシーケンスに従って実行するシーケンス制御部と、
複数の前記せん断波速度の加算平均を計算する加算平均部と、
受付部と、を備え、
前記シーケンス制御部は、前記計測を繰り返し、
前記せん断波検出部は、前記計測毎に前記せん断波を検出し、
前記速度計算部は、前記せん断波を検出する毎に前記せん断波速度を計算し、
前記加算平均部は、前記せん断波を検出する毎に計算される複数の前記せん断波速度の加算平均を計算し、
前記提示部は、前記加算平均結果を前記信頼性情報とともにユーザに提示し、さらに、前記計測毎の前記信頼性情報および前記せん断波速度を散布図にして提示し、
前記受付部は、ユーザから前記散布図上のプロット結果を介して指示を受け付け、
前記加算平均部は、前記指示に応じて、前記加算平均を再計算し、
前記受付部は、除外する前記せん断波速度の選択をさらに受け付け、
前記加算平均部は、前記選択されたせん断波速度以外のせん断波速度を用いて、前記加算平均を再計算すること
を特徴とする超音波診断装置。
【請求項１０】
被検体に集束したバースト超音波を送信することにより当該バースト超音波の焦点位置に生成されたせん断波を、複数のせん断波検出用パルスを繰り返し送信することにより得たエコー信号群を用いて検出するせん断波検出部と、
前記せん断波の伝搬速度であるせん断波速度を計算する速度計算部と、
前記せん断波の伝搬領域を含む計測領域の変動を評価し、評価結果を、前記せん断波速度の信頼度を示す信頼性情報として得る変動評価部と、
前記信頼性情報をユーザに提示する提示部と、
前記バースト超音波の送信と前記複数のせん断波検出用パルスの繰り返しの送信と当該送信によるエコー信号の受信とからなる計測を、予め定めたパルスシーケンスに従って実行するシーケンス制御部と、
複数の前記せん断波速度の加算平均を計算する加算平均部と、
受付部と、を備え、
前記シーケンス制御部は、前記計測を繰り返し、
前記せん断波検出部は、前記計測毎に前記せん断波を検出し、
前記速度計算部は、前記せん断波を検出する毎に前記せん断波速度を計算し、
前記加算平均部は、前記せん断波を検出する毎に計算される複数の前記せん断波速度の加算平均を計算し、
前記提示部は、前記加算平均結果を前記信頼性情報とともにユーザに提示し、さらに、前記計測毎の前記信頼性情報および前記せん断波速度を散布図にして提示し、
前記受付部は、ユーザから前記散布図上のプロット結果を介して指示を受け付け、
前記加算平均部は、前記指示に応じて、前記加算平均を再計算し、
前記受付部は、前記せん断波速度を前記信頼性情報に応じて複数のグループに分割する指示をさらに受け付け、
前記加算平均部は、グループごとに前記せん断波速度の加算平均を再計算すること
を特徴とする超音波診断装置。
【請求項１１】
請求項５記載の超音波診断装置であって、
ユーザから再計測の指示を受け付ける受付部をさらに備え、
前記シーケンス制御部は、前記指示に応じて、前記計測を実行すること
を特徴とする超音波診断装置。
【請求項１２】
請求項２記載の超音波診断装置であって、
前記変動評価部は、前記せん断波が伝搬する深度領域を、当該せん断波の発生位置と前記せん断波の振幅とにより特定すること
を特徴とする超音波診断装置。
【請求項１３】
請求項２記載の超音波診断装置であって、
前記変動評価部は、前記せん断波が伝搬する深度領域を、前記エコー信号群から得たデータに対して時間方向に相関演算を行って得た相関係数を用いて特定すること
を特徴とする超音波診断装置。
【請求項１４】
被検体に集束したバースト超音波を送信することにより当該バースト超音波の焦点位置に生成されたせん断波を、複数のせん断波検出用パルスを繰り返し送信することにより得たエコー信号群を用いて検出するせん断波検出部と、
前記せん断波の伝搬速度であるせん断波速度を計算する速度計算部と、
前記せん断波の伝搬領域を含む計測領域の変動を評価し、評価結果を、前記せん断波速度の信頼度を示す信頼性情報として得る変動評価部と、
前記信頼性情報をユーザに提示する提示部と、を備え
前記変動評価部は、さらに、前記変動の要因を特定し、
前記提示部は、前記特定された要因をさらにユーザに提示し、
前記変動評価部は、前記エコー信号群から得たデータに対し、時間方向に相関演算を行って得た相関係数の変化のパターンにより前記要因を特定すること
を特徴とする超音波診断装置。
【請求項１５】
被検体の所定の焦点位置に集束したバースト超音波を送信し、生成されたせん断波を前記被検体に伝搬させ、複数のせん断波検出用パルスを前記被検体に対して繰り返し送信し、そのエコーを受信することにより、エコー信号群を得て、
前記エコー信号群から、前記被検体の前記せん断波による変動を検出して、前記せん断波の速度を計算するとともに、前記被検体の前記せん断波以外の変動を、前記せん断波による変動とは区別して検出し、
前記せん断波以外の変動の検出結果をユーザに提示すること
を特徴とする弾性評価方法。
【国際調査報告】
2015029651.tif 000015
2015029651.tif 000016