JP2015131097A - 超音波診断装置、画像処理装置及び画像処理方法 - Google Patents

Abstract

【課題】生体組織の硬さの信頼性を超音波画像上に表すことができる超音波診断装置の提供。
【解決手段】超音波診断装置１０の送信部１１０は、音響放射力に基づいて生体組織に変位を発生させる変位発生用超音波を超音波プローブ１０１から送信させ、変位発生用超音波に基づいて発生する所定の走査領域内の生体組織の変位を観測する観測用超音波を超音波プローブ１０１から送信させる。受信部１２０は、超音波プローブ１０１で受信された反射波に基づいて、反射波データを生成する。画像生成部１４０は、反射波データに基づいて、走査領域内の複数位置それぞれにおける変位を複数の時相で算出し、算出された変位が略最大となる時相を複数位置それぞれで決定し、複数位置のうち決定された時相が略同一の位置同士を表す画像データを生成する。表示制御部１７３は、画像データに基づく画像を、走査領域を含む領域に対応する医用画像に重畳表示させる。
【選択図】図１

Description

本発明の実施形態は、超音波診断装置、画像処理装置及び画像処理方法に関する。

従来、生体組織の硬さを測定し、測定した硬さの分布を映像化するエラストグラフィー（Elastography）が知られている。エラストグラフィーは、例えば、肝硬変等、病変の進行度に応じて生体組織の硬さが変化する疾患の診断に利用されている。エラストグラフィーにおいて、生体組織を変位させて硬さを評価する方法は、以下の２つに大別される。

１つ目の方法は、超音波プローブ（probe）で体表から生体組織を圧迫・開放した際に観測される走査断面内の各点の歪みの大きさから、相対的な硬さを可視化する方法である。また、２つ目の方法は、体表から生体組織に音響放射力や機械的振動を与えてせん断波（shear wave）による変位を発生させ、走査断面内の各点における変位を経時的に観測することで、せん断波の伝播速度（propagation speed）を求めて、弾性率を求める方法である。前者の方法では、歪みの局所的な大きさは、超音波プローブを手動で動かす大きさに依存し、周囲と比較した時に着目領域が相対的に硬いか軟らかいかという評価が行われる。これに対して、後者の方法では、着目領域の絶対的な弾性率を求めることができる。

ところで、被検体の体動によって変位が発生した場合、或いはせん断波が組織内で反射・屈折した場合等には、生体組織の硬さが正確に測定できないことがある。このような場合、エラストグラフィーによって映像化された硬さ画像は、信頼性の低いものとなってしまう。

Ｔ． Ｄｅｆｆｉｅｕｘ， Ｊ．−Ｌ． Ｇｅｎｎｉｓｓｏｎ， Ｊ． Ｂｅｒｃｏｆｆ， ａｎｄ Ｍ． Ｔａｎｔｅｒ， 「Ｏｎ ｔｈｅ ｅｆｆｅｃｔｓ ｏｆ ｒｅｆｌｅｃｔｅｄ ｗａｖｅｓ ｉｎ ｔｒａｎｓｉｅｎｔ ｓｈｅａｒ ｗａｖｅ ｅｌａｓｔｏｇｒａｐｈｙ」， ＩＥＥＥ Ｔｒａｎｓ． ＵＦＦＣ， ｖｏｌ．５８， Ｉｓｓｕｅ １０， Ｏｃｔｏｂｅｒ ２０１１，ｐｐ．２０３２−２０３５

本発明が解決しようとする課題は、生体組織の硬さの信頼性を超音波画像上に表すことができる超音波診断装置、画像処理装置及び画像処理方法を提供することである。

実施形態に係る超音波診断装置は、送信部と、受信部と、画像生成部と、表示制御部とを備える。送信部は、音響放射力に基づいて生体組織に変位を発生させる変位発生用超音波を超音波プローブから送信させ、前記変位発生用超音波に基づいて発生する所定の走査領域内の生体組織の変位を観測する観測用超音波を前記超音波プローブから送信させる。受信部は、前記超音波プローブで受信された反射波に基づいて、反射波データを生成する。画像生成部は、前記反射波データに基づいて、前記走査領域内の複数位置それぞれにおける変位を複数の時相で算出し、算出された変位が略最大となる時相を前記複数位置それぞれで決定し、前記複数位置のうち決定された時相が略同一の位置同士を表す画像データを生成する。表示制御部は、前記画像データに基づく画像を、前記走査領域を含む領域に対応する医用画像に重畳表示させる。

図１は、第１の実施形態に係る超音波診断装置の構成例を示すブロック図である。 図２は、第１の実施形態に係る硬さ画像データの一例を示す図である。 図３は、第１の実施形態に係る線画像生成部の処理を説明するための図である。 図４は、第１の実施形態に係る線画像生成部の処理を説明するための図である。 図５は、第１の実施形態に係る表示制御部によって表示される表示画像の一例を示す図である。 図６は、第１の実施形態に係る超音波診断装置の処理手順の一例を説明するためのフローチャートである。 図７は、第２の実施形態に係る超音波診断装置の構成例を示すブロック図である。 図８は、第２の実施形態に係る算出部の処理を説明するための図である。 図９は、第２の実施形態に係る算出部の処理を説明するための図である。 図１０は、第２の実施形態に係る線画像生成部の処理を説明するための図である。 図１１は、第３の実施形態に係る線画像生成部の処理を説明するための図である。 図１２は、第３の実施形態に係る表示制御部によって表示される表示画像の一例を示す図である。 図１３は、第４の実施形態に係る線画像生成部の処理を説明するための図である。 図１４は、第４の実施形態に係る表示制御部によって表示される表示画像の一例を示す図である。 図１５は、変位画像データについて説明するための図である。 図１６は、明るさを変更する処理を説明するための図である。 図１７は、非表示とする処理を説明するための図である。 図１８は、各種パラメータの評価について説明するための図である。 図１９は、プッシュパルスの送信位置及び観測方向の表示について説明するための図である。 図２０は、並列表示について説明するための図である。 図２１は、並列表示される場合にモニタに表示される表示画像の一例を示す図である。

以下、図面を参照して、実施形態に係る超音波診断装置、画像処理装置及び画像処理方法を説明する。

（第１の実施形態）
まず、第１の実施形態に係る超音波診断装置の構成について説明する。図１は、第１の実施形態に係る超音波診断装置の構成例を示すブロック図である。図１に示すように、第１の実施形態に係る超音波診断装置１０は、装置本体１００と、超音波プローブ１０１と、入力装置１０２と、モニタ（monitor）１０３とを有する。

超音波プローブ１０１は、複数の振動子（例えば、圧電振動子）を有し、これら複数の振動子は、後述する装置本体１００が有する送信部１１０から供給される駆動信号に基づき超音波を発生する。また、超音波プローブ１０１が有する複数の振動子は、被検体Ｐからの反射波を受信して電気信号に変換する。また、超音波プローブ１０１は、振動子に設けられる整合層と、振動子から後方への超音波の伝播を防止するバッキング材等を有する。

超音波プローブ１０１から被検体Ｐに超音波が送信されると、送信された超音波は、被検体Ｐの体内組織における音響インピーダンス（impedance）の不連続面で次々と反射され、反射波信号として超音波プローブ１０１が有する複数の振動子にて受信される。受信される反射波信号の振幅は、超音波が反射される不連続面における音響インピーダンスの差に依存する。なお、送信された超音波パルスが、移動している血流や心臓壁等の表面で反射された場合の反射波信号は、ドプラ（Doppler）効果により、移動体の超音波送信方向に対する速度成分に依存して、周波数偏移を受ける。

なお、第１の実施形態は、図１に示す超音波プローブ１０１が、複数の圧電振動子が一列で配置された１次元超音波プローブである場合や、一列に配置された複数の圧電振動子が機械的に揺動される１次元超音波プローブである場合、複数の圧電振動子が格子状に２次元で配置された２次元超音波プローブである場合のいずれであっても適用可能である。

入力装置１０２は、マウス（mouse）、キーボード（keyboard）、ボタン（button）、パネルスイッチ（panel switch）、タッチコマンドスクリーン（touch command screen）、フットスイッチ（footswitch）、トラックボール（trackball）、ジョイスティック（joystick）等を有し、超音波診断装置１０の操作者からの各種設定要求を受け付け、装置本体１００に対して受け付けた各種設定要求を転送する。

モニタ１０３は、超音波診断装置１０の操作者が入力装置１０２を用いて各種設定要求を入力するためのＧＵＩ（Graphical User Interface）を表示したり、装置本体１００において生成された超音波画像データ等を表示したりする。

装置本体１００は、超音波プローブ１０１が受信した反射波信号に基づいて超音波画像データを生成する装置であり、図１に示すように、送信部１１０と、受信部１２０と、信号処理部１３０と、画像生成部１４０と、画像メモリ１５０と、内部記憶部１６０と、制御部１７０とを有する。

送信部１１０は、超音波送信における送信指向性を制御する。具体的には、送信部１１０は、レートパルサ（rate pulser）発生器、送信遅延部、送信パルサ等を有し、超音波プローブ１０１に駆動信号を供給する。レートパルサ発生器は、所定のレート周波数（ＰＲＦ：Pulse Repetition Frequency）で、送信超音波を形成するためのレートパルスを繰り返し発生する。レートパルスは、送信遅延部を通ることで異なる送信遅延時間を有した状態で送信パルサへ電圧を印加する。すなわち、送信遅延部は、超音波プローブ１０１から発生される超音波をビーム状に集束して送信指向性を決定するために必要な振動子ごとの送信遅延時間を、レートパルサ発生器が発生する各レートパルスに対し与える。送信パルサは、かかるレートパルスに基づくタイミング（timing）で、超音波プローブ１０１に駆動信号（駆動パルス）を印加する。送信方向或いは送信遅延時間は、後述する内部記憶部１６０に記憶されており、送信部１１０は、内部記憶部１６０を参照して、送信指向性を制御する。

駆動パルスは、送信パルサからケーブル（cable）を介して超音波プローブ１０１内の振動子まで伝達した後に、振動子において電気信号から機械的振動に変換される。この機械的振動は、生体内部で超音波として送信される。振動子ごとに異なる送信遅延時間を持った超音波は、収束されて、所定方向に伝搬していく。送信遅延部は、各レートパルスに対し与える送信遅延時間を変化させることで、振動子面からの送信方向を任意に調整する。送信部１１０は、超音波ビーム（beam）の送信に用いる振動子の数及び位置（送信開口）と、送信開口を構成する各振動子の位置に応じた送信遅延時間とを制御することで、送信指向性を与える。例えば、送信部１１０の送信遅延回路は、送信遅延時間をパルサ（pulser）回路が発生する各レートパルスに対し与えることで、超音波送信の深さ方向における集束点（送信フォーカス（focus））の位置を制御する。

なお、送信部１１０は、後述する制御部１７０の指示に基づいて、所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧等を瞬時に変更可能な機能を有している。特に、送信駆動電圧の変更は、瞬間にその値を切り替え可能なリニアアンプ（linear-amp）型の発信回路、または、複数の電源ユニットを電気的に切り替える機構によって実現される。

超音波プローブ１０１が送信した超音波の反射波は、超音波プローブ１０１内部の振動子まで到達した後、振動子において、機械的振動から電気的信号（反射波信号）に変換され、受信部１２０に入力される。

受信部１２０は、超音波受信における受信指向性を制御する。具体的には、受信部１２０は、プリアンプ（pre-amp）、Ａ／Ｄ変換部、受信遅延部及び加算部等を有し、超音波プローブ１０１が受信した反射波信号に対して各種処理を行なって反射波データを生成する。プリアンプは、反射波信号をチャンネル（channel）ごとに増幅してゲイン（gain）補正処理を行なう。Ａ／Ｄ変換部は、ゲイン補正された反射波信号をＡ／Ｄ変換し、受信遅延部は、受信指向性を決定するのに必要な受信遅延時間をチャンネルごとに与える。加算部は、受信遅延時間が与えられた反射波信号（デジタル（digital）信号）を加算して、反射波データを生成する。加算器の加算処理により、反射波信号の受信指向性に応じた方向からの反射成分が強調される。受信方向或いは受信遅延時間は、後述する内部記憶部１６０に記憶されており、受信部１２０は、内部記憶部１６０を参照して、受信指向性を制御する。なお、第１の実施形態に係る受信部１２０は、並列同時受信を行なうことも可能である。

信号処理部１３０は、受信部１２０が反射波信号から生成した反射波データに対して各種の信号処理を行う。信号処理部１３０は、受信部１２０から受信した反射波データに対して、対数増幅、包絡線検波処理等を行って、サンプル（sample）点（観測点）ごとの信号強度が輝度の明るさで表現されるデータ（Ｂモードデータ（mode data））を生成する。

また、信号処理部１３０は、受信部１２０から受信した反射波データより、移動体のドプラ効果に基づく運動情報を、走査領域内の各サンプル点で抽出したデータ（ドプラデータ）を生成する。具体的には、信号処理部１３０は、移動体の運動情報として、平均速度（average velocity）、分散値、パワー値等を各サンプル点で抽出したドプラデータを生成する。ここで、移動体とは、例えば、血流や、心壁等の組織、造影剤である。

ここで、第１の実施形態に係る超音波診断装置１０は、生体組織の硬さを測定し、測定した硬さの分布を映像化するエラストグラフィーを実行可能な装置である。具体的には、第１の実施形態に係る超音波診断装置１０は、音響放射力を与えて生体組織に変位を発生させることで、エラストグラフィーを実行可能な装置である。

すなわち、第１の実施形態に係る送信部１１０は、音響放射力で生じるせん断波により変位を発生させるプッシュパルス（push pulse）（変位発生用バースト波）を超音波プローブ１０１から送信させる。そして、第１の実施形態に係る送信部１１０は、プッシュパルスにより発生する変位を観測する観測用パルスを、走査領域内の複数の走査線それぞれで超音波プローブ１０１から複数回送信させる。観測用パルスは、プッシュパルスにより発生したせん断波の伝播速度を、走査領域内の各サンプル点で観測するために送信される。通常、観測用パルスは、走査領域内の各走査線に対して、複数回（例えば、１００回）、送信される。受信部１２０は、走査領域内の各走査線で送信された観測用パルスの反射波信号から、反射波データを生成する。なお、プッシュパルスは、変位発生用超音波の一例である。また、観測用パルスは、観測用超音波の一例である。

言い換えると、送信部１１０は、音響放射力に基づいて生体組織に変位を発生させる変位発生用超音波を超音波プローブから送信させ、変位発生用超音波に基づいて発生する所定の走査領域内の生体組織の変位を観測する観測用超音波を超音波プローブから送信させる。また、受信部１２０は、超音波プローブで受信された反射波に基づいて、反射波データを生成する。

そして、信号処理部１３０は、走査領域内の各走査線で複数回送信された観測用パルスの反射波データを解析して走査領域の硬さの分布を示す硬さ分布情報を算出する。具体的には、信号処理部１３０は、プッシュパルスによって発生したせん断波の伝搬速度を各サンプル点で測定することで、走査領域の硬さ分布情報を生成する。

例えば、信号処理部１３０は、観測用パルスの反射波データを周波数解析する。これにより、信号処理部１３０は、各走査線の複数のサンプル点それぞれで、運動情報（組織ドプラデータ）を複数時相に渡って生成する。そして、信号処理部１３０は、各走査線の複数のサンプル点それぞれで得られた複数時相の組織ドプラデータの速度成分を時間積分する。これにより、信号処理部１３０は、各走査線の複数のサンプル点それぞれの変位を複数時相に渡って算出する。続いて、信号処理部１３０は、各サンプル点で変位が最大となる時間を求める。そして、信号処理部１３０は、各サンプル点で最大変位が得られた時間を、各サンプル点にせん断波が到達した到達時間として決定する。続いて、信号処理部１３０は、各サンプル点におけるせん断波の到達時間の空間的微分を行うことで、各サンプル点でのせん断波の伝搬速度を算出する。以下、「せん断波の伝播速度」を、「せん断速度（share wave speed）」と記載する。なお、上記の到達時間としては、各サンプル点で変位が最大となる時間ではなく、例えば、各サンプル点における変位の変化量が最大となる時間を用いても良い。

そして、信号処理部１３０は、せん断速度をカラーコード化し、対応するサンプル点にマッピング（mapping）することで、硬さ分布情報を生成する。硬い組織ではせん断速度が大きく、柔らかい組織ではせん断速度が小さくなる。すなわち、せん断速度の値は、組織の硬さ（弾性率）を示す値となる。上記の場合、観測用パルスは、組織ドプラ用の送信パルスである。なお、上記のせん断速度は、各サンプル点で変位が最大となる時間（到達時間）に基づくのではなく、例えば、信号処理部１３０が、隣接する走査線における組織の変位の相互相関により検出することで算出する場合であっても良い。

なお、信号処理部１３０は、せん断速度から、ヤング（Young）率又はせん断弾性率を算出し、算出したヤング率又はせん断弾性率により硬さ分布情報を生成しても良い。せん断速度、ヤング率及びせん断弾性率は、いずれも生体組織の硬さを表す物理量（指標値）として用いることができる。以下では、信号処理部１３０が、生体組織の硬さを表す物理量としてヤング率を用いる場合について説明する。

ここで、１回のプッシュパルスの送信で発生するせん断波は、伝播とともに、減衰する。広い領域に渡ってせん断速度を観測しようとした場合、ある特定の一つの走査線において送信されたプッシュパルスにより発生したせん断波は、伝播に伴い減衰し、プッシュパルスの送信位置から十分離れると、やがて観測不可能となる。

かかる場合、プッシュパルスを方位方向における複数の位置において送信する必要がある。具体的には、走査領域（或いは、関心領域）は、方位方向に沿って複数の領域に分割される。送信部１１０は、分割された各領域（以下、分割領域と表記）において観測用パルスを送受信する前に、それぞれ異なった走査線位置においてプッシュパルスを送信し、せん断波を発生させる。この際、典型的には、プッシュパルスの送信位置は、各分割領域の近傍に設定される。同時並列受信数が少数に限定されている場合、送信部１１０は、複数の分割領域それぞれにおいて、プッシュパルスを１回送信した後に、ある分割領域の各走査線で観測用パルスを複数回送信する処理を順次実行する。

画像生成部１４０は、信号処理部１３０が生成したデータから超音波画像データを生成する。画像生成部１４０は、信号処理部１３０が生成したＢモードデータから反射波の強度を輝度で表したＢモード画像データを生成する。また、画像生成部１４０は、信号処理部１３０が生成したドプラデータから移動体情報を表すドプラ画像データを生成する。ドプラ画像データは、速度画像データ、分散画像データ、パワー画像データ、又は、これらを組み合わせた画像データである。

また、画像生成部１４０は、信号処理部１３０によって生成された硬さ分布情報から、生体組織の硬さがカラー表示された硬さ画像データを生成する。例えば、画像生成部１４０は、走査領域内の各点（各サンプル点）におけるせん断速度に応じた画素値を各点に割り当てたせん断速度画像データ（share wave speed image data）を、硬さ画像データとして生成する。

図２は、第１の実施形態に係る硬さ画像データの一例を示す図である。図２に示すように、画像生成部１４０は、生体組織の硬さが色分けされた画像を生成する。なお、図２では、説明の都合上、エラストグラフィーにより映像化される対象として走査領域全体が指定された場合を例示したが、実施形態はこれに限定されるものではない。例えば、映像化される対象として走査領域の一部である関心領域（Region Of Interest：ＲＯＩ）が指定された場合には、ＲＯＩに対応する硬さ画像データが生成される。なお、第１の実施形態において、ＲＯＩは、エラストグラフィーによる映像化の対象となる対象領域に対応する。

ここで、画像生成部１４０は、一般的には、超音波走査の走査線信号列を、テレビ（TV）等に代表されるビデオフォーマット（video format）の走査線信号列に変換（スキャンコンバート（scan convert））し、表示用の超音波画像データを生成する。具体的には、画像生成部１４０は、超音波プローブ１０１による超音波の走査形態に応じて座標変換を行うことで、表示用の超音波画像データを生成する。また、画像生成部１４０は、スキャンコンバート以外に種々の画像処理として、例えば、スキャンコンバート後の複数の画像フレームを用いて、輝度の平均値画像を再生成する画像処理（平滑化処理）や、画像内で微分フィルタを用いる画像処理（エッジ強調処理）等を行う。また、画像生成部１４０は、超音波画像データに、付帯情報（種々のパラメータの文字情報、目盛り、ボディーマーク（body mark）等）を合成する。

すなわち、Ｂモードデータ、ドプラデータ及び硬さ分布情報は、スキャンコンバート処理前の超音波画像データであり、画像生成部１４０が生成するデータは、スキャンコンバート処理後の表示用の超音波画像データである。なお、画像生成部１４０は、信号処理部１３０が３次元のデータ（３次元Ｂモードデータ、３次元ドプラデータ及び３次元硬さ分布情報）を生成した場合、超音波プローブ１０１による超音波の走査形態に応じて座標変換を行うことで、ボリュームデータ（volume data）を生成する。そして、画像生成部１４０は、ボリュームデータに対して、各種レンダリング処理を行なって、表示用の２次元画像データを生成する。

画像メモリ１５０は、画像生成部１４０が生成した表示用の画像データを記憶するメモリである。また、画像メモリ１５０は、信号処理部１３０が生成したデータを記憶することも可能である。画像メモリ１５０が記憶するＢモードデータやドプラデータ、硬さ分布情報は、例えば、診断の後に操作者が呼び出すことが可能となっており、画像生成部１４０を経由して表示用の超音波画像データとなる。

内部記憶部１６０は、超音波送受信、画像処理及び表示処理を行うための制御プログラムや、診断情報（例えば、患者ＩＤ、医師の所見等）や、診断プロトコル（protocol）や各種ボディーマーク等の各種データを記憶する。また、内部記憶部１６０は、必要に応じて、画像メモリ１５０が記憶する画像データの保管等にも使用される。また、内部記憶部１６０が記憶するデータは、図示しないインタフェース部を介して、外部装置へ転送することができる。

また、内部記憶部１６０は、撮影済みのせん断速度画像データに関する情報を記憶する。例えば、内部記憶部１６０は、撮影済みのせん断速度画像データについて、各サンプル点におけるせん断波の到達時間を記憶する。

制御部１７０は、超音波診断装置１０の処理全体を制御する。具体的には、制御部１７０は、入力装置１０２を介して操作者から入力された各種設定要求や、内部記憶部１６０から読込んだ各種制御プログラム及び各種データに基づき、送信部１１０、受信部１２０、信号処理部１３０、画像生成部１４０の処理を制御する。また、制御部１７０は、画像メモリ１５０が記憶する表示用の超音波画像データをモニタ１０３にて表示するように制御する。

なお、装置本体１００に内蔵される送信部１１０、受信部１２０、及び制御部１７０等は、プロセッサ（ＣＰＵ（Central Processing Unit）、ＭＰＵ（Micro-Processing Unit）、集積回路等）のハードウェア（hardware）により構成されることもあるが、ソフトウェア（software）的にモジュール（module）化されたプログラムにより構成される場合もある。

ところで、被検体の体動によって変位が発生した場合、或いはせん断波が組織内で反射・屈折した場合等には、生体組織の硬さが正確に測定できないことがある。このような場合、エラストグラフィーによって映像化された硬さ画像は、信頼性の低いものとなってしまう。

そこで、第１の実施形態に係る超音波診断装置１０は、硬さ画像上にその信頼性を表すことができる。この機能を実現するための超音波診断装置１０の構成について、以下に説明する。

第１の実施形態に係る超音波診断装置１０において、制御部１７０は、送信制御部１７１と、線画像生成部１７２と、表示制御部１７３とを備える。

送信制御部１７１は、送信部１１０によるプッシュパルスの送信や、観測用パルスの送受信を制御する。例えば、送信制御部１７１は、操作者からＲＯＩを決定する指示を受け付ける。そして、送信制御部１７１は、受け付けた指示に基づいて、ＲＯＩに対応する硬さ画像データを生成するためのプッシュパルスの送信位置、送信位置の数、ＲＯＩの位置、ＲＯＩの範囲、及びＲＯＩの数等を設定する。送信制御部１７１の制御により、送信部１１０は、プッシュパルスを超音波プローブ１０１から送信させる。また、送信制御部１７１の制御により、送信部１１０は、送信させたプッシュパルスにより発生する変位を観測する観測用パルスを、走査領域内の複数の走査線それぞれで超音波プローブ１０１から複数回送信させる。

線画像生成部１７２は、走査領域内の各点にせん断波が到達した到達時間が略同一の線を表す線画像データを生成する。この線画像データは、例えば、到達時間が略同一の位置を線で結んだ等高線を表示するための情報である。例えば、線画像生成部１７２は、走査領域内の各点（各サンプル点）のうち、所定の到達時間を有する点を複数抽出する。そして、線画像生成部１７２は、抽出した複数の点を結ぶことで、線画像データを生成する。なお、所定の到達時間とは、超音波診断装置１０の操作者若しくは設計者によって予め指定された値であり、通常、複数の値が指定される。なお、線画像生成部１７２は、到達時間画像生成部、若しくは画像生成部の一例である。

言い換えると、線画像生成部１７２は、反射波データに基づいて、走査領域内の複数位置それぞれにおける変位を複数の時相で算出し、算出された変位が略最大となる時相を複数位置それぞれで決定し、複数位置のうち決定された時相が略同一の位置同士を表す画像データを生成する。具体的には、線画像生成部１７２は、離散的な複数時相に対応する前記画像データを生成する。

図３及び図４は、第１の実施形態に係る線画像生成部１７２の処理を説明するための図である。図３には、走査領域内の各点に対して、各点における到達時間に応じた画素値を割り当てた到達時間画像データの一例を示す。図４には、線画像生成部１７２によって生成された線画像データ２０ａ，２０ｂ，２０ｃ，２０ｄが、図３の到達時間画像データ上に表示された場合を例示する。なお、線画像データ２０ａは、到達時間が１０［ｍｓｅｃ］の位置を示し、線画像データ２０ｂは、到達時間が２０［ｍｓｅｃ］の位置を示し、線画像データ２０ｃは、到達時間が３０［ｍｓｅｃ］の位置を示し、線画像データ２０ｄは、到達時間が４０［ｍｓｅｃ］の位置を示す。また、線画像データ２０ａ，２０ｂ，２０ｃ，２０ｄを区別無く総称する場合に、「線画像データ２０」と表記する。この線画像データ２０は、到達時間位置画像データの一例である。

ここで、まず、図３の到達時間画像データについて説明する。到達時間画像データは、例えば、画像生成部１４０によって生成される。具体的には、画像生成部１４０は、信号処理部１３０によって決定された各点の到達時間に応じた画素値を、走査領域内の各点に対して割り当てることで、到達時間画像データを生成する。

そして、線画像生成部１７２は、画像生成部１４０によって生成された到達時間画像データを取得する。そして、線画像生成部１７２は、取得した到達時間画像データから、到達時間が１０［ｍｓｅｃ］である点を複数抽出する。そして、線画像生成部１７２は、抽出した複数の点を結ぶことで、線画像データ２０ａを生成する。また、線画像生成部１７２は、線画像データ２０ａと同様に、２０［ｍｓｅｃ］の点を結ぶことで線画像データ２０ｂを生成し、３０［ｍｓｅｃ］の点を結ぶことで線画像データ２０ｃを生成し、４０［ｍｓｅｃ］の点を結ぶことで線画像データ２０ｄを生成する。

また、線画像生成部１７２は、例えば、それぞれの線画像データ２０に対して、到達時間に応じた画素値を割り当てる。図４に示す例では、線画像生成部１７２は、線画像データ２０ａ，２０ｂ，２０ｃ，２０ｄそれぞれに対して、その到達時間に応じた画素値を割り当てる。一例としては、線画像生成部１７２は、図３の到達時間画像データにおいて、各点の到達時間に応じて割り当てられた色を、線画像データ２０ａ，２０ｂ，２０ｃ，２０ｄに割り当てる。具体的には、線画像生成部１７２は、線画像データ２０ａに青色を割り当て、線画像データ２０ｂに緑色を割り当て、線画像データ２０ｃに黄色を割り当て、線画像データ２０ｄに赤色を割り当てる。

このように、線画像生成部１７２は、到達時間が１０，２０，３０，４０［ｍｓｅｃ］である位置をそれぞれ示す線画像データ２０ａ，２０ｂ，２０ｃ，２０ｄをそれぞれ生成する。

なお、図４はあくまで一例に過ぎない。例えば、図４の例では、到達時間が１０，２０，３０，４０［ｍｓｅｃ］である位置に線画像データ２０を生成する場合を説明したが、実施形態はこれに限定されるものではない。例えば、線画像生成部１７２は、任意の到達時間の位置に線画像データ２０を生成しても良いし、任意の本数の線画像データ２０を生成しても良い。

また、図４の例では、それぞれの線画像データ２０に対して、到達時間に応じた画素値を割り当てる場合を説明したが、実施形態はこれに限定されるものではない。例えば、線画像生成部１７２は、到達時間がそれぞれ異なる線画像データ２０ごとに、それぞれ異なる線種（実線、破線、点線、一点鎖線等）を割り当てても良い。

また、ここでは説明の都合上、到達時間画像データ（図３）が生成された後に線画像データ２０が生成される場合を説明したが、実施形態はこれに限定されるものではない。すなわち、線画像生成部１７２は、必ずしも到達時間画像データが生成されていなくとも、線画像データ２０を生成することができる。例えば、線画像生成部１７２は、信号処理部１３０で決定された各点の到達時間のうち、所定の到達時間を有する点を複数抽出することで、線画像データ２０を生成することができる。

また、ここでは説明を省略したが、線画像生成部１７２は、線画像データ２０を生成する場合に、例えば、線画像データ２０を滑らかにするための平滑化処理を行っても良いし、線画像データ２０の太さを均一にする処理を行っても良い。また、線画像生成部１７２が所定の到達時間を有する点を抽出する場合、必ずしも所定の到達時間と厳密に一致する点のみを抽出するのではなく、例えば、抽出された複数の点が線形を成す程度に誤差を許容して抽出しても良い。すなわち、上記の「到達時間が略同一の線」とは、抽出された複数の点が線形を成す程度の誤差の範囲内で同一の到達時間を有する複数の点に基づく線であり、平滑化処理や太さを均一にする処理が行われた線である。

表示制御部１７３は、線画像データを超音波画像データ上に重畳表示させる。例えば、表示制御部１７３は、線画像生成部１７２によって生成された線画像データ２０を、硬さ画像データ上に重畳表示させる。

言い換えると、表示制御部１７３は、線画像生成部１７２によって生成された画像データに基づく画像を、走査領域を含む領域に対応する医用画像に重畳表示させる。

図５は、第１の実施形態に係る表示制御部１７３によって表示される表示画像の一例を示す図である。図５に示すように、表示制御部１７３は、図２に例示の硬さ画像データ上に、図４で生成した線画像データ２０ａ，２０ｂ，２０ｃ，２０ｄをそれぞれ表示させる。

なお、図５はあくまで一例に過ぎない。例えば、図５では、硬さ画像データを背景画像として線画像データ２０を重畳表示する場合を説明したが、背景画像はこれに限定されない。例えば、背景画像としては、Ｂモード画像データや、図３の到達時間画像データ、更には後述する変位画像データ及び分散値画像データを用いても良い。

図６は、第１の実施形態に係る超音波診断装置１０の処理手順の一例を説明するためのフローチャートである。なお、図６に示す例では、超音波診断装置１０は、最初はＢモードであり、Ｂモード画像データを表示している。

図６に示すように、第１の実施形態に係る超音波診断装置１０の送信制御部１７１は、硬さ画像を生成するための硬さ画像生成モードを開始させる開始指示を操作者から受け付けたか否かを判定する（ステップＳ１０１）。この硬さ画像生成モードとは、例えば、硬さ画像を生成するＲＯＩを設定したり、ＲＯＩを設定した後にプッシュパルスを送信して硬さ画像を生成したりするための状態である。ここで、開始指示を受け付けない場合（ステップＳ１０１否定）、送信制御部１７１は、開始指示を受け付けるまで待機する。

一方、開始指示を受け付けた場合（ステップＳ１０１肯定）、送信制御部１７１の制御により、モニタ１０３は、ＲＯＩ設定用ＧＵＩを表示する（ステップＳ１０２）。このＲＯＩ設定用ＧＵＩとしては、例えば、エラストグラフィーによる映像化の対象となる対象領域を指定するＲＯＩが、Ｂモード画像データ上に表示されている。このＲＯＩの位置及び大きさはプリセットされている。そして、送信制御部１７１は、操作者からＲＯＩの位置及び大きさを変更する指示を受け付けて、受け付けた指示に応じてＲＯＩの位置及び大きさを変更する（ステップＳ１０３）。

そして、送信制御部１７１は、ＲＯＩが決定されたか否かを判定する（ステップＳ１０４）。例えば、送信制御部１７１は、操作者からＲＯＩの位置及び大きさを決定する旨の指示を受け付けたか否かに応じて、ＲＯＩが決定されたか否かを判定する。ここで、ＲＯＩが決定されない場合（ステップＳ１０４否定）、送信制御部１７１は、ＲＯＩが決定されるまで、ステップＳ１０３の処理を実行する。

一方、ＲＯＩが決定された場合（ステップＳ１０４肯定）、送信制御部１７１は、決定されたＲＯＩに基づいて、プッシュパルスの送信位置を決定する（ステップＳ１０５）。例えば、送信制御部１７１は、ＲＯＩの位置及び大きさに基づいて、プッシュパルスの送信位置を１箇所又は複数箇所決定する。これは、プッシュパルスにより発生するせん断波が伝播の過程で減衰することが知られており、この減衰によってＲＯＩの全域を走査できない事態を防ぐためである。

そして、送信制御部１７１は、操作者から硬さ画像データの撮影開始要求を受け付けたか否かを判定する（ステップＳ１０６）。ここで、撮影開始要求を受け付けない場合（ステップＳ１０６否定）、送信制御部１７１は、撮影開始要求を受け付けるまで待機する。

一方、撮影開始要求を受け付けた場合（ステップＳ１０６肯定）、送信制御部１７１は、決定した送信位置において、超音波プローブ１０１からプッシュパルスを送信させる（ステップＳ１０７）。そして、送信部１１０及び受信部１２０の制御により、超音波プローブ１０１は、ＲＯＩ内で観測用パルスの送受信を実行する（ステップＳ１０８）。例えば、観測用パルスは、ＲＯＩ内のある走査線に対して複数回（１００回程度）送受信される。これにより、各点における変位の時間変化を算出する。１つのパルスに対して多数の受信を行うことができるシステムを有している場合、１度のプッシュパルスの送信でＲＯＩ内の全領域の変位の時間変化を知ることができるが、同時受信本数が限られている場合は、複数回の観測用パルスの送受信を、ラスタ位置を変えて複数回行われる。その際は、ラスタ位置を変えて観測用パルスを送信するたびに、プッシュパルスを送信する。

続いて、信号処理部１３０は、ＲＯＩの各点の変位を複数時相に渡って算出する（ステップＳ１０９）。そして、信号処理部１３０は、ＲＯＩの各点におけるせん断波の到達時間を決定する（ステップＳ１１０）。例えば、信号処理部１３０は、各点で最大変位が得られた時間を、各点におけるせん断波の到達時間として決定する。そして、信号処理部１３０は、各点におけるせん断波の到達時間の空間的微分を行うことで、各点でのせん断波の伝搬速度（せん断速度）を算出する（ステップＳ１１１）。そして、信号処理部１３０は、ＲＯＩに対応するせん断速度をカラーコード化し、対応する各点にマッピングすることで、硬さ分布情報を生成する（ステップＳ１１２）。

続いて、画像生成部１４０は、信号処理部１３０によって生成された硬さ分布情報から、生体組織の硬さがカラー表示された硬さ画像データを生成する（ステップＳ１１３）。例えば、画像生成部１４０は、走査領域内の各点（各サンプル点）におけるせん断速度に応じた画素値を各点に割り当てたせん断速度画像データを、硬さ画像データとして生成する。そして、表示制御部１７３は、画像生成部１４０によって生成された硬さ画像データをモニタ１０３に表示させる（ステップＳ１１４）。

そして、線画像生成部１７２は、信頼性を表示する旨の信頼性表示要求を受け付けたか否かを判定する（ステップＳ１１５）。ここで、信頼性表示要求受け付けた場合（ステップＳ１１５肯定）、線画像生成部１７２は、ＲＯＩ内の各点のうち所定の到達時間を有する複数の点を用いて、線画像データ２０を生成する（ステップＳ１１６）。例えば、線画像生成部１７２は、ＲＯＩ内の各点のうち、所定の到達時間を有する点を複数抽出し、抽出した複数の点を結ぶことで、線画像データ２０を生成する。そして、表示制御部１７３は、線画像生成部１７２によって生成された線画像データ２０を、硬さ画像データ上に重畳表示させる（ステップＳ１１７）。なお、信頼性表示要求を受け付けない場合（ステップＳ１１５否定）、ステップＳ１１８に移行する。

そして、送信制御部１７１は、硬さ画像生成モードを終了させる終了指示を操作者から受け付けたか否かを判定する（ステップＳ１１８）。ここで、終了指示を受け付けない場合（ステップＳ１１８否定）、送信制御部１７１は、終了指示を受け付けるまで、硬さ画像データを表示させたまま、信頼性表示要求の受け付けを待機する。

一方、送信制御部１７１は、終了指示を受け付けた場合（ステップＳ１１８肯定）、硬さ画像生成モードを終了させる。例えば、送信制御部１７１は、硬さ画像データを非表示とし、Ｂモードに戻る。

なお、上記の処理手順はあくまで一例に過ぎず、実施形態は図６の処理手順に限定されるものではない。例えば、ＲＯＩの位置及び大きさを変更する処理（ステップＳ１０３）は、必ずしも実行されなくても良い。例えば、プリセットのままのＲＯＩが利用される場合には、ステップＳ１０３の処理は実行されなくても良い。また、例えば、ステップＳ１０３の処理は、ＲＯＩの位置又は大きさのいずれか一方のみが変更される場合であっても良い。

また、上記の処理手順では、信頼性表示要求を受け付けた後に（ステップＳ１１５肯定）、線画像データ２０を生成する処理（ステップＳ１１６）が実行される場合を説明したが、これに限定されるものではない。例えば、線画像データ２０を生成する処理は、バックグラウンドの処理として実行しておき、信頼性表示要求を受け付けた場合に、線画像データ２０が表示されても良い。すなわち、線画像データ２０を生成する処理は、到達時間を決定する処理（ステップＳ１１０）の実行後であれば、任意のタイミングで実行されて良い。

また、上記の処理手順では、信頼性表示要求を受け付けた後に（ステップＳ１１５肯定）、線画像データ２０を表示する処理（ステップＳ１１７）が実行される場合を説明したが、これに限定されるものではない。例えば、信頼性表示要求を受け付けなくとも、自動的に線画像データ２０を表示しても良い。すなわち、ステップＳ１１５の処理は実行されなくても良い。

上述してきたように、第１の実施形態に係る超音波診断装置１０は、走査領域内の各点にせん断波が到達した到達時間が略同一の線を表す線画像データを生成する。そして、超音波診断装置１０は、線画像データを超音波画像データ上に重畳表示させる。

例えば、超音波診断装置１０は、ＲＯＩ内の各点のうち、所定の到達時間を有する点を複数抽出し、抽出した複数の点を結ぶことで、線画像データ２０を生成する。これによれば、超音波診断装置１０は、あたかも地図上に等高線を引くかの如く、到達時間画像データ上に到達時間の目盛り線となる線画像データ２０を表示することができる。

ここで、到達時間を用いて線画像データ２０を表示するのは、到達時間が生体組織の硬さの信頼性の指標となるからである。例えば、被検体の体動による変位の発生、或いは、せん断波の反射・屈折等が無く、生体組織の硬さが正確に測定できる環境下であれば、せん断波は、プッシュパルスの送信位置からほぼ均等に伝播するはずである。この結果、到達時間の線画像データ２０は、プッシュパルスの送信方向とほぼ並行となり、生体組織の硬さに応じて曲線を描くこととなる。これに対して、生体組織の硬さが正確に測定できない環境下であれば、せん断波の伝播が極端に早く観測されたり、極端に遅く観測されたりしてしまう。この結果、到達時間の線画像データ２０は、大きく湾曲してしまう。

したがって、操作者は、超音波画像データに重畳表示された線画像データ２０の曲がり方の程度を確認することで、線画像データ２０の近傍の領域における生体組織の硬さの信頼性を判断することができる。

また、例えば、超音波診断装置１０は、それぞれの線画像データ２０に対して、到達時間に応じた画素値を割り当てる。このため、操作者は、それぞれの線画像データ２０の画素値（色）を確認することで、その線画像データ２０付近の到達時間の値を知ることができる。

（第２の実施形態）
第１の実施形態では、線画像データ２０に対して、その到達時間に応じた画素値を割り当てる場合を説明したが、実施形態はこれに限定されるものではない。例えば、超音波診断装置１０は、線画像データ２０に対して、他のパラメータに応じた画素値を割り当てても良い。そこで、第２の実施形態では、超音波診断装置１０が、線画像データ２０に対して、他のパラメータに応じた画素値を割り当てる場合を説明する。

図７は、第２の実施形態に係る超音波診断装置１０の構成例を示すブロック図である。第２の実施形態に係る超音波診断装置１０は、図１に例示した超音波診断装置１０と同様の構成を備え、算出部１７４を更に備える点と、線画像生成部１７２の処理の一部とが相違する。そこで、第２の実施形態では、第１の実施形態と相違する点を中心に説明することとし、第１の実施形態において説明した構成と同様の機能を有する点については、図１と同一の符号を付し、説明を省略する。

第２の実施形態に係る算出部１７４は、線画像データ２０に含まれる点ごとに、その点を含む所定領域内の各点の到達時間の分散値を算出する。

図８及び図９は、第２の実施形態に係る算出部１７４の処理を説明するための図である。図８には、図４と同様の到達時間画像データ及び線画像データ２０の上に、分散値算出用領域２１が配置された場合を例示する。ここで、分散値算出用領域２１は、分散値を算出するために、線画像データ２０に含まれる各点を中心とし、所定の大きさを有する領域である。この大きさはプリセットされているが、プリセットの値は任意に変更可能である。また、図９には、図８の分散値算出用領域２１付近の拡大図を例示する。図９において、黒の矩形は、線画像データ２０ｃに含まれる各点（画素）を示す。また、網掛けの矩形は、分散値算出用領域２１の中心点（中心画素）を示す。

図８に示すように、算出部１７４は、例えば、線画像データ２０ｃに含まれる点を中心として、分散値算出用領域２１を配置する。そして、図９に示すように、算出部１７４は、配置された分散値算出用領域２１に含まれる各点の到達時間の分散値を、中心点の分散値として算出する。そして、算出部１７４は、線画像データ２０ｃに含まれる他の点についても同様に、各点を中心として分散値算出用領域２１を配置し、その中に含まれる各点の分散値を算出する。これにより、算出部１７４は、線画像データ２０ｃに含まれる各点の分散値を算出する。

また、算出部１７４は、線画像データ２０ａ，２０ｂ，２０ｄについても同様に、各線画像データ２０ａ，２０ｂ，２０ｄに含まれる各点の分散値をそれぞれ算出する。これにより、線画像データ２０の各点は、その周辺の点の到達時間のバラツキに対応するパラメータとして、到達時間の分散値を有することとなる。

なお、図８及び図９はあくまで一例に過ぎない。例えば、ここでは説明の都合上、到達時間画像データ（図３）上で各点の分散値が算出される場合を説明したが、実施形態はこれに限定されるものではない。すなわち、算出部１７４は、必ずしも到達時間画像データ上で処理を行わなくとも、線画像データ２０に含まれる各点の分散値を算出することができる。

第２の実施形態に係る線画像生成部１７２は、第１の実施形態において説明した機能と同様の機能を有し、更に、線画像データ２０に含まれる各点に、その点の分散値に応じた画素値を割り当てる。例えば、線画像生成部１７２は、算出部１７４によって算出された各点の分散値に応じた画素値を、線画像データ２０に含まれる各点に割り当てる。

図１０は、第２の実施形態に係る線画像生成部１７２の処理を説明するための図である。図１０には、図９の線画像データ２０ｃに含まれる各点に対して、各点の分散値に応じた画素値を割り当てたもの一例を示す。

図１０に示すように、線画像生成部１７２は、図９の線画像データ２０ｃに含まれる各点に対して、算出部１７４によって算出された各点の分散値に応じた画素値を割り当てる。一例としては、線画像生成部１７２は、各点のうち分散値が低いものから順に、青色、青緑色、黄緑色、黄色、橙色、赤色を割り当てる。

なお、図１０はあくまで一例に過ぎない。例えば、図１０の例では、線画像データ２０に含まれる各点に対して、各点の分散値に応じた画素値を割り当てる場合を説明したが、実施形態はこれに限定されるものではない。例えば、線画像生成部１７２は、線画像データ２０に含まれる各点に対して、各点の分散値ごとに異なるハッチングパターンを割り当てても良い。

また、線画像生成部１７２は、線画像データ２０ａ，２０ｂ，２０ｄについても同様に、各線画像データ２０ａ，２０ｂ，２０ｄに含まれる各点の分散値に応じた画素値を割り当てる。これにより、モニタ１０３に表示される線画像データ２０の各点は、その周辺の点の到達時間のバラツキに応じた色を有することとなる。

このように、第２の実施形態に係る超音波診断装置１０は、線画像データ２０に含まれる点ごとに、その点を含む所定領域内の各点の到達時間の分散値を算出する。そして、超音波診断装置１０は、線画像データ２０に含まれる各点に、その点の分散値に応じた画素値を割り当てる。

ここで、分散値を用いて線画像データ２０に含まれる各点に画素値を割り当てる（色分けする）のは、分散値が生体組織の硬さの信頼性の指標となるからである。例えば、第２の実施形態において、モニタ１０３に表示される線画像データ２０の各点は、その周辺の点の到達時間のバラツキに応じた色を有する。ここで、生体組織の硬さが正確に測定できる環境下であれば、せん断波がほぼ均等に伝播する結果、分散値は高くはならない。これに対して、せん断波の伝搬が生体組織内で乱れた場合、分散値は高くなる。この分散値が高くなった領域での硬さの信頼性は低いと判断できる。

したがって、操作者は、線画像データ２０に含まれる各点に割り当てられた画素値（色）を確認することで、線画像データ２０の近傍の領域における生体組織の硬さの信頼性を判断することができる。

なお、第２の実施形態では、線画像データ２０に含まれる各点に、その点の分散値に応じた画素値を割り当てる場合を説明したが、実施形態はこれに限定されるものではない。例えば、超音波診断装置１０は、線画像データ２０の各点に対して、各点の生体組織の硬さの指標値（せん断速度、ヤング率、せん断弾性率）、或いは、各点における変位の大きさ等のパラメータに応じた画素値を割り当てても良い。これにより、操作者は、線画像データ２０の各点の色を確認することで、各種のパラメータによって示唆される情報を容易に認知することができる。

（第３の実施形態）
上述した第１及び第２の実施形態では、到達時間を用いて線画像データ２０を表示する場合を例示したが、実施形態はこれに限定されるものではない。例えば、超音波診断装置１０は、他のパラメータを用いて線画像データ２０のような目盛り線となる画像を表示しても良い。そこで、第３の実施形態では、超音波診断装置１０が、他のパラメータを用いて線画像データ２０のような目盛り線となる画像を表示する場合を説明する。

第３の実施形態に係る超音波診断装置１０は、図１に例示した超音波診断装置１０と同様の構成を備え、線画像生成部１７２及び表示制御部１７３の処理の一部が相違する。そこで、第３の実施形態では、第１の実施形態と相違する点を中心に説明することとし、第１の実施形態において説明した構成と同様の機能を有する点については、図１と同一の符号を付し、説明を省略する。

第３の実施形態に係る線画像生成部１７２は、第１の実施形態において説明した機能と同様の機能を有し、更に、走査領域内の各点のせん断波に基づく生体組織の硬さの指標値が略同一の線を表す硬さ線画像データを生成する。例えば、線画像生成部１７２は、走査領域内の各点のうち、所定の大きさの指標値を有する点を複数抽出する。そして、線画像生成部１７２は、抽出した複数の点を結ぶことで、硬さ線画像データを生成する。なお、所定の大きさの指標値とは、超音波診断装置１０の操作者若しくは設計者によって予め指定された値であり、通常、複数の値が指定される。

図１１は、第３の実施形態に係る線画像生成部１７２の処理を説明するための図である。図１１には、線画像生成部１７２によって生成された線画像データ２２ａ，２２ｂが、図２の硬さ画像データ上に表示された場合を例示する。なお、線画像データ２２ａは、ヤング率が３０［ｋＰａ］の位置を示し、線画像データ２２ｂは、ヤング率が６０［ｋＰａ］の位置を示す。また、線画像データ２２ａ，２２ｂを区別無く総称する場合に、「線画像データ２２」と表記する。この線画像データ２２は、硬さ位置画像データの一例である。

図１１に示すように、線画像生成部１７２は、画像生成部１４０によって生成された硬さ画像データを取得する。そして、線画像生成部１７２は、取得した硬さ画像データから、ヤング率が３０［ｋＰａ］である点を複数抽出する。そして、線画像生成部１７２は、抽出した複数の点を結ぶことで、線画像データ２２ａを生成する。また、線画像生成部１７２は、線画像データ２２ａと同様に、６０［ｋＰａ］の点を結ぶことで線画像データ２２ｂを生成する。

また、線画像生成部１７２は、例えば、それぞれの線画像データ２２に対して、硬さの指標値（ここではヤング率）に応じた画素値を割り当てる。図１１に示す例では、線画像生成部１７２は、線画像データ２２ａ，２２ｂそれぞれに対して、そのヤング率に応じた画素値を割り当てる。具体的には、線画像生成部１７２は、線画像データ２２ａに対して黄緑色を割り当て、線画像データ２２ｂに対して橙色を割り当てる。

このように、線画像生成部１７２は、ヤング率が３０，６０［ｋＰａ］である位置をそれぞれ示す線画像データ２２ａ，２２ｂをそれぞれ生成する。

なお、図１１はあくまで一例に過ぎない。例えば、図１１の例では、ヤング率が３０，６０［ｋＰａ］である位置に線画像データ２２を生成する場合を説明したが、実施形態はこれに限定されるものではない。例えば、線画像生成部１７２は、任意のヤング率の位置に線画像データ２２を生成しても良いし、任意の本数の線画像データ２２を生成しても良い。

また、図１１の例では、それぞれの線画像データ２２に対して、ヤング率に応じた画素値を割り当てる場合を説明したが、実施形態はこれに限定されるものではない。例えば、線画像生成部１７２は、ヤング率がそれぞれ異なる線画像データ２２ごとに、それぞれ異なる線種（実線、破線、点線、一点鎖線等）を割り当てても良い。

また、ここでは説明の都合上、硬さ画像データ（図２）が生成された後に線画像データ２２が生成される場合を説明したが、実施形態はこれに限定されるものではない。すなわち、線画像生成部１７２は、必ずしも硬さ画像データが生成されていなくとも、線画像データ２２を生成することができる。例えば、線画像生成部１７２は、信号処理部１３０で決定された各点のヤング率のうち、所定のヤング率を有する点を複数抽出することで、線画像データ２２を生成することができる。

また、ここでは説明を省略したが、線画像生成部１７２は、線画像データ２２を生成する場合に、例えば、線画像データ２２を滑らかにするための平滑化処理を行っても良いし、線画像データ２２の太さを均一にする処理を行っても良い。また、線画像生成部１７２が所定のヤング率を有する点を抽出する場合、必ずしも所定のヤング率と厳密に一致する点のみを抽出するのではなく、例えば、抽出された複数の点が線形を成す程度に誤差を許容して抽出しても良い。すなわち、上記の「硬さの指標値が略同一の線」とは、抽出された複数の点が線形を成す程度の誤差の範囲内で同一の硬さの指標値を有する複数の点に基づく線であり、平滑化処理や太さを均一にする処理が行われた線である。

第３の実施形態に係る表示制御部１７３は、第１の実施形態において説明した機能と同様の機能を有し、更に、線画像データ２２を超音波画像データ上に重畳表示させる。例えば、表示制御部１７３は、線画像生成部１７２によって生成された線画像データ２２を、到達時間画像データ上に重畳表示させる。

図１２は、第３の実施形態に係る表示制御部１７３によって表示される表示画像の一例を示す図である。図１２に示すように、表示制御部１７３は、図３に例示の到達時間画像データ上に、図１１で生成した線画像データ２２ａ，２２ｂをそれぞれ表示させる。

なお、図１２はあくまで一例に過ぎない。例えば、図１２では、到達時間画像データを背景画像として線画像データ２２を重畳表示する場合を説明したが、背景画像はこれに限定されない。例えば、背景画像としては、Ｂモード画像データや、図２の硬さ画像データ、更には後述する変位画像データ及び分散値画像データを用いても良い。

このように、第３の実施形態に係る超音波診断装置１０は、硬さの指標値を用いて、線画像データ２２を生成する。そして、超音波診断装置１０は、生成した線画像データ２２を超音波画像データ上に重畳表示させる。これによれば、操作者は、例えば、硬さの指標値と背景画像のパラメータとの関連性を視覚的に把握することができる。

なお、第３の実施形態では、到達時間以外のパラメータとして、硬さの指標値を用いて目盛り線（等高線）を生成する場合を説明したが、実施形態はこれに限定されるものではない。例えば、線画像生成部１７２は、各点の変位の大きさを用いて目盛り線を生成しても良い。具体的には、線画像生成部１７２は、走査領域内の各点の変位の大きさが略同一の線を表す画像データ（変位位置画像データ）を生成しても良い。そして、表示制御部１７３は、生成した変位位置画像データを超音波画像データ上に表示させる。これにより、超音波診断装置１０は、背景画像と異なるパラメータを用いて線画像データ２０のような目盛り線を表示するので、操作者は、これらの異なるパラメータの関連性を視覚的に把握することができる。

また、第３の実施形態に係る超音波診断装置１０は、到達時間の線画像データ２０、線画像データ２２、及び変位位置画像データのうち、任意の画像データを選択し、選択した全ての画像データを超音波画像データ上に重畳表示させても良い。すなわち、超音波診断装置１０は、例えば、到達時間の線画像データ２０及び線画像データ２２を、一つの超音波画像データ上に重畳表示させても良い。

また、第３の実施形態に係る超音波診断装置１０が、他のパラメータを用いて線画像データ２０のような目盛り線となる画像を表示する処理は、到達時間の線画像データ２０を生成したり表示したりする処理とは別に独立して実行可能である。

（第４の実施形態）
また、例えば、超音波診断装置１０は、所定領域内のせん断速度で走査領域全体にせん断波が伝播する場合の到達時間に基づいて、線画像データ２０のような目盛り線となる画像を表示しても良い。そこで、第４の実施形態では、超音波診断装置１０が、所定領域内のせん断速度で走査領域全体にせん断波が伝播する場合の到達時間に基づいて、線画像データ２０のような目盛り線となる画像を表示する場合を説明する。

第４の実施形態に係る超音波診断装置１０は、図１に例示した超音波診断装置１０と同様の構成を備え、線画像生成部１７２及び表示制御部１７３の処理の一部が相違する。そこで、第４の実施形態では、第１の実施形態と相違する点を中心に説明することとし、第１の実施形態において説明した構成と同様の機能を有する点については、図１と同一の符号を付し、説明を省略する。

第４の実施形態に係る線画像生成部１７２は、第１の実施形態において説明した機能と同様の機能を有する。更に、線画像生成部１７２は、所定領域におけるせん断波の伝播速度で、せん断波が走査領域内の各点に到達する到達時間を算出し、算出した到達時間が略同一の線を表す所定領域位置画像データを生成する。

図１３は、第４の実施形態に係る線画像生成部１７２の処理を説明するための図である。図１３には、図２の硬さ画像データ上に、操作者により指定されたＲＯＩ２３と、線画像データ２４ａ〜２４ｆとが表示された場合を例示する。なお、図１３において、線画像データ２４ａは、到達時間が５［ｍｓｅｃ］の位置を示し、線画像データ２４ｂは、到達時間が１０［ｍｓｅｃ］の位置を示し、線画像データ２４ｃは、到達時間が１５［ｍｓｅｃ］の位置を示し、線画像データ２４ｄは、到達時間が２０［ｍｓｅｃ］の位置を示し、線画像データ２４ｅは、到達時間が２５［ｍｓｅｃ］の位置を示し、線画像データ２４ｆは、到達時間が３０［ｍｓｅｃ］の位置を示す。また、線画像データ２４ａ〜２４ｆを区別無く総称する場合に、「線画像データ２４」と表記する。この線画像データ２４は、所定領域位置画像データの一例である。

図１３に示すように、線画像生成部１７２は、例えば、操作者からＲＯＩ２３を指定する操作を受け付ける。そして、線画像生成部１７２は、受け付けたＲＯＩ２３に含まれる各点におけるせん断波速度の平均値を算出する。そして、線画像生成部１７２は、算出したせん断速度の平均値で、せん断波が走査領域内の各点に到達する到達時間を算出する。言い換えると、線画像生成部１７２は、走査領域内の全ての点が均一な硬さ（ＲＯＩ２３内の平均的な硬さ）であり、その走査領域全体をせん断波が一定の速度（算出したせん断速度の平均値）で伝播するものとして、走査領域内の各点の到達時間を算出する。例えば、線画像生成部１７２は、算出したせん断速度の平均値と、プッシュパルスの送信位置から各点までの距離とを用いて、各点の到達時間を算出する。そして、線画像生成部１７２は、算出した到達時間が略同一の線を表す線画像データ２４ａ〜２４ｆをそれぞれ生成する。ここで、線画像データ２４ａ〜２４ｆは、走査領域内の全ての点がＲＯＩ２３内の平均的な硬さであるものとして生成されるので、プッシュパルスの送信位置に対して平行な線となる。なお、図１３では、走査領域の左端にプッシュパルスの送信位置があるものとする。また、線画像生成部１７２が、算出した到達時間が略同一の位置を表す線画像データ２４ａ〜２４ｆをそれぞれ生成する処理は、第１の実施形態で説明した内容と同様であるので、詳細な説明は省略する。

なお、図１３はあくまで一例に過ぎない。例えば、線画像生成部１７２は、任意の到達時間の位置に線画像データ２４を生成しても良いし、任意の本数の線画像データ２４を生成しても良い。

また、例えば、線画像生成部１７２は、それぞれの線画像データ２４に対して、算出した到達時間に応じた画素値を割り当てても良い。また、例えば、線画像生成部１７２は、算出した到達時間がそれぞれ異なる線画像データ２４ごとに、それぞれ異なる線種（実線、破線、点線、一点鎖線等）を割り当てても良い。

また、例えば、ＲＯＩ２３の位置及び大きさがプリセットされている場合には、必ずしも操作者によるＲＯＩ２３の指定を要しない。例えば、プリセットのままのＲＯＩ２３が利用される場合には、操作者によるＲＯＩ２３の指定は受け付けなくても良い。

また、図１３では、硬さ画像データが背景画像として用いられる場合を説明したが、実施形態はこれに限定されるものではない。例えば、線画像生成部１７２は、背景画像の有無、及び、背景画像となる超音波画像データの種類に関係なく、上記の処理を実行することができる。

また、図１３では、ＲＯＩ２３に含まれる各点のせん断速度の平均値を用いる場合を説明したが、これに限定されるものではなく、例えば、ＲＯＩ２３の中心点のせん断速度を用いる場合であっても良い。

第４の実施形態に係る表示制御部１７３は、第１の実施形態において説明した機能と同様の機能を有し、更に、線画像データ２４を超音波画像データ上に表示させる。例えば、表示制御部１７３は、線画像生成部１７２によって生成された線画像データ２４を、到達時間画像データ上に重畳表示させる。

図１４は、第４の実施形態に係る表示制御部１７３によって表示される表示画像の一例を示す図である。図１４に示すように、表示制御部１７３は、図３に例示の到達時間画像データ上に、図１３で生成した線画像データ２４をそれぞれ表示させる。

なお、図１４はあくまで一例に過ぎない。例えば、図１４では、到達時間画像データを背景画像として線画像データ２４を重畳表示する場合を説明したが、背景画像はこれに限定されない。例えば、背景画像としては、Ｂモード画像データや、図２の硬さ画像データ、更には後述する変位画像データ及び分散値画像データを用いても良い。

このように、第４の実施形態に係る超音波診断装置１０は、ＲＯＩ２３内のせん断速度で走査領域全体にせん断波が伝播する場合の到達時間に基づいて、線画像データ２０のような目盛り線となる画像を表示する。これによれば、操作者は、ＲＯＩ２３内のせん断速度で走査領域全体にせん断波が伝播する場合の到達時間を、他のパラメータと比較検討することが容易となる。

また、第４の実施形態に係る超音波診断装置１０は、ここで表示される線画像データ２４を、到達時間の線画像データ２０、線画像データ２２、及び変位位置画像データのうち任意の画像データとともに超音波画像データ上に重畳表示されても良い。これにより、操作者は、複数のパラメータを同一の超音波画像データ上で比較検討することができる。

また、第４の実施形態に係る超音波診断装置１０が、所定領域内のせん断速度で走査領域全体にせん断波が伝播する場合の到達時間に基づいて、線画像データ２０のような目盛り線となる画像を表示する処理は、到達時間の線画像データ２０を生成したり表示したりする処理とは別に独立して実行可能である。

（その他の実施形態）
さて、これまで第１〜第４の実施形態について説明したが、これ以外にも、種々の異なる形態にて実施されて良い。なお、以下に説明する各処理は、到達時間の線画像データ２０を生成したり表示したりする処理とは別に独立して実行可能である。

（変位画像データの生成）
例えば、超音波診断装置１０は、走査領域内の各点に対して、各点における変位の大きさに応じた画素値を割り当てた変位画像データを生成しても良い。

例えば、画像生成部１４０は、走査領域内の各点に対して、各点における変位の大きさに応じた画素値を割り当てることで、変位画像データを生成する。具体的には、画像生成部１４０は、信号処理部１３０によって複数時相に渡って算出された走査領域の各点における変位を取得する。そして、画像生成部１４０は、各点について、複数時相に渡って算出された変位のうちで最大となる変位量（最大変位量）を特定する。そして、画像生成部１４０は、特定した各点の最大変位量に応じた画素値を各点に割り当てることで、変位画像データを生成する。

図１５は、変位画像データについて説明するための図である。図１５に示すように、画像生成部１４０は、走査領域内の各点が最大変位量によって色分けされた変位画像データを生成する。

このように、画像生成部１４０は、走査領域内の各点に対して、各点における変位の大きさに応じた画素値を割り当てることで、変位画像データを生成する。なお、ここで、変位画像データを生成するのは、変位の大きさが生体組織の硬さの信頼性の指標となるからである。例えば、走査領域内の各点において、観測用パルスで観測するのに十分な大きさの変位が発生していれば、観測された生体組織の硬さの信頼性は高くなる。これに対して、走査領域内の各点において、十分な大きさの変位が発生していなければ、観測用パルスで観測するのが困難であるので、観測された生体組織の硬さの信頼性は低くなる。したがって、操作者は、変位画像データを確認することで、観測された生体組織の硬さの信頼性を容易に把握することができる。

図１５に示す例では、走査領域内の上方の点ほど変位が大きいので、信頼性が高いことがわかる。一方、走査領域内の下方の点ほど変位が小さいので、信頼性が低いことがわかる。

なお、図１５の例では、各点における変位の大きさとして最大変位量を用いる場合を説明したが、実施形態はこれに限定されるものではない。例えば、時間−変位曲線における積分値（一定時間内の変位の合計値）に応じた画素値を各点に割り当てることで、変位画像データを生成しても良い。また、例えば、時間−変位曲線における微分値（傾きの大きさ）に応じた画素値を各点に割り当てることで、変位画像データを生成しても良い。すなわち、せん断波の伝搬速度が最大となった時刻に応じた画素値を各点に割り当てることで、変位画像データを生成しても良い。

（分散値画像データの生成）
例えば、超音波診断装置１０は、走査領域内の各点に対して、各点における分散値に応じた画素値を割り当てた分散値画像データを生成しても良い。

例えば、算出部１７４は、走査領域内の各点について、その点を含む所定領域内の各点の到達時間の分散値を算出する。一例としては、算出部１７４は、図８の分散値算出用領域２１を用いて、走査領域内の各点の分散値を算出する。なお、分散値算出用領域２１を用いて分散値を算出する処理は、第２の実施形態で説明した内容と同様であるので、詳細な説明は省略する。

そして、画像生成部１４０は、走査領域内の各点に対して、各点における分散値に応じた画素値を割り当てることで、分散値画像データを生成する。具体的には、画像生成部１４０は、算出部１７４によって算出された走査領域内の各点の分散値を取得する。そして、画像生成部１４０は、取得した分散値に応じた画素値を各点に割り当てることで、分散値画像データを生成する。

このように、画像生成部１４０は、走査領域内の各点に対して、各点における分散値に応じた画素値を割り当てることで、分散値画像データを生成する。これにより、操作者は、分散値画像データを確認することで、観測された生体組織の硬さの信頼性を容易に把握することができる。

（各種パラメータによる画像データの変更）
また、例えば、超音波診断装置１０は、超音波画像データに含まれる各点の画素値を、各点における硬さの指標値、せん断波の到達時間、変位の大きさ、及び、分散値の各パラメータのうちいずれか一つのパラメータを用いて変更しても良い。

画像生成部１４０は、超音波画像データに含まれる各点の画素値を、各点における硬さの指標値、せん断波の到達時間、変位の大きさ、及び、分散値の各パラメータのうちいずれか一つのパラメータを用いて変更する。

例えば、画像生成部１４０は、超音波画像データの明るさを変更する処理を行う。図１６は、明るさを変更する処理を説明するための図である。図１６では、変位の大きさを用いて硬さ画像データの明るさを変更する場合の処理を説明する。図１６に示す例では、硬さ画像データの領域２５ａが暗くなり、領域２５ｂが黒く塗りつぶされている。なお、硬さ画像データの上方の領域（領域２５ａよりも上の領域）は、元の硬さ画像データの明るさのまま表示されている。

例えば、画像生成部１４０は、走査領域内の各点における変位の大きさに基づいて、信頼性を３段階で評価する。具体例を挙げると、画像生成部１４０は、最大変位量が第１閾値以上である点の信頼性を「高」と判定する。また、画像生成部１４０は、最大変位量が第１閾値未満、かつ、第２閾値以上である点の信頼性を「中」と判定する。ここで、第２閾値は第１閾値より小さい値である。また、画像生成部１４０は、最大変位量が第２閾値未満である点の信頼性を「低」と判定する。図１６に示す例では、画像生成部１４０は、硬さ画像データの上方の領域を「高」と判定し、領域２５ａを「中」と判定し、領域２５ｂを「低」と判定する。

そして、画像生成部１４０は、信頼性が「高」と判定された領域については、元の画像の明るさのままとし、変更しない。また、画像生成部１４０は、信頼性が「中」と判定された領域２５ａについては、元の画像を暗くする。例えば、画像生成部１４０は、領域２５ａに含まれる点のＲＧＢ値を所定数減算することで、暗くする。また、画像生成部１４０は、信頼性が「低」と判定された領域２５ｂについては、元の画像を黒く塗りつぶす。例えば、画像生成部１４０は、領域２５ｂに含まれる点のＲＧＢ値を０とすることで、黒く塗りつぶす。このように、画像生成部１４０は、図２の硬さ画像データの明るさを、各点の変位の大きさに基づいて変更する。

また、例えば、画像生成部１４０は、超音波画像データの一部を非表示とする処理を行う。図１７は、非表示とする処理を説明するための図である。図１７では、変位の大きさを用いて硬さ画像データの一部を非表示とする場合の処理を説明する。図１７に示す例では、硬さ画像データの領域２６が白くなり、何も表示されていない。

例えば、画像生成部１４０は、走査領域内の各点における変位の大きさに基づいて、信頼性を３段階で評価する。この評価は、例えば、図１６で説明した方法と同様で良い。例えば、画像生成部１４０は、硬さ画像データの上方の領域を「高」と判定し、領域２５ａを「中」と判定し、領域２５ｂを「低」と判定する。

そして、画像生成部１４０は、信頼性が「高」と判定された領域のみ、硬さ画像データを生成し、信頼性が「中」及び「低」と判定された領域２６については生成しない。このように、画像生成部１４０は、領域２６に対応する硬さ画像データを生成しないことで、領域２６を非表示とする。

このように、画像生成部１４０は、例えば、硬さ画像データに含まれる各点の画素値を、各点における変位の大きさを用いて変更する。なお、上記の例では、信頼性を３段階で評価する場合を説明したが、これに限定されるものではない。例えば、画像生成部１４０は、一つの閾値を用いることで、信頼性を２段階で評価しても良いし、複数の閾値を用いることで、信頼性を多段階評価しても良い。

また、ここでは、変位の大きさを用いて硬さ画像データを変更する場合を説明したが、これに限定されるものではない。すなわち、上述した各種の超音波画像データに含まれる各点の画素値を、各点における硬さの指標値、せん断波の到達時間、変位の大きさ、及び、分散値の各パラメータのうちいずれか一つのパラメータを用いて変更することができる。

これにより、超音波診断装置１０は、生体組織の硬さやその信頼性など、各種のパラメータによって表現される情報を、他の超音波画像データ上に表現することができる。

（各種パラメータの評価）
また、例えば、超音波診断装置１０は、関心領域内の各点における各種のパラメータを評価し、評価結果を表示させても良い。

例えば、表示制御部１７３は、関心領域内の各点におけるせん断波に基づく生体組織の硬さの指標値、到達時間、変位の大きさ、及び到達時間の分散値のうち少なくとも一つのパラメータを評価し、評価結果を表示させる。

図１８は、各種パラメータの評価について説明するための図である。図１８には、図２の硬さ画像データと、その画像上に指定された信頼性評価用のＲＯＩ２７とを例示する。

例えば、表示制御部１７３は、操作者によりＲＯＩ２７の指定を受け付ける。そして、表示制御部１７３は、受け付けたＲＯＩ２７内の各点における各種パラメータを評価する。

一例としては、表示制御部１７３は、ＲＯＩ２７内の各点における分散値と、閾値とを比較することで、信頼性を評価する。具体的には、表示制御部１７３は、ＲＯＩ２７内の各点における分散値の平均値を算出する。そして、表示制御部１７３は、算出した平均値と閾値とを比較して、平均値が閾値未満であれば、信頼度「高」と評価し、平均値が閾値以上であれば、信頼度「低」と評価する。例えば、表示制御部１７３は、図１８に示すように、評価結果「信頼度：高い」をモニタ１０３に表示させる。

なお、図１８の例はあくまで一例である。例えば、表示制御部１７３は、ＲＯＩ２７内の各点における分散値の平均値ではなく、ＲＯＩ２７内の分散値の最大値を閾値と比較することで、評価しても良い。また、モニタ１０３には評価結果のみならず、パラメータそのものが表示されても良い。例えば、表示制御部１７３は、ＲＯＩ２７の各点の到達時間の平均値を算出する。そして、表示制御部１７３は、算出した平均値「８．０ｍｓｅｃ」をモニタ１０３上に表示させても良い。

（プッシュパルスの送信位置及び観測方向の表示）
また、例えば、超音波診断装置１０は、プッシュパルスの送信位置及び観測方向を表示しても良い。

図１９は、プッシュパルスの送信位置及び観測方向の表示について説明するための図である。図１９には、図２の硬さ画像データ上に、フラグ２８ａ，２８ｂ，２８ｃと、ライン２９ａ，２９ｂ，２９ｃとを例示する。なお、図１９に示す例では、プッシュパルスがライン２９ａ，２９ｂ，２９ｃ上で送信され、それぞれの送信位置から右方向に伝播するせん断波が観測された場合を説明する。

図１９に示すように、画像生成部１４０は、送信制御部１７１からプッシュパルスの送信位置と、その観測方向とをそれぞれ示す情報を取得する。そして、画像生成部１４０は、左端のプッシュパルスの送信位置に対応する画像データとして、ライン２９ａを生成する。このライン２９ａは、左端のプッシュパルスの送信位置に配置される。また、画像生成部１４０は、左端のプッシュパルスの観測方向を示すために、フラグ２８ａを生成する。この場合、観測方向はプッシュパルスの送信位置の右側であるので、画像生成部１４０は、フラグ２８ａをライン２９ａの右側に配置する。このとき、画像生成部１４０は、フラグ２８ａの短辺をライン２９ａに重ねる。これにより、ライン２９ａから右方向に向かうフラグ２８ａによって、ライン２９ａの送信位置から右方向に伝播するせん断波が観測されたことを明示することができる。

（並列表示）
また、例えば、超音波診断装置１０は、上述した各種の超音波画像データと、それに重畳表示される目盛り線（線画像データ２０等）とを並列表示しても良い。

例えば、表示制御部１７３は、硬さ画像データ、到達時間画像データ、変位画像データ、及び分散値画像データのうち、少なくとも二つの画像データを並列表示させる。そして、表示制御部１７３は、並列表示させた画像データのうち、少なくとも一つの画像データに、到達時間の線画像データ２０、線画像データ２２、変位位置画像データ、及び線画像データ２４のうち少なくとも一つを重畳表示させる。

図２０は、並列表示について説明するための図である。図２０の左側には、硬さ画像データ上に到達時間の線画像データ２０を重畳表示させた画像を示し、右側には、到達時間画像データ上に線画像データ２４を重畳表示させた画像を示す。なお、ここに例示した画像データは一例に過ぎず、表示制御部１７３は、各種の超音波画像データや、各種パラメータにより変更した画像データ（図１６，１７等）を並列表示して良い。また、表示制御部１７３は、上述した到達時間の線画像データ２０、線画像データ２２、変位位置画像データ、及び線画像データ２４のうち任意の画像データを重畳表示して良い。これにより、各種画像データの並列表示が可能となる。

（線画像データ２０の動画表示）
また、例えば、超音波診断装置１０は、線画像データ２０を動画として表示しても良い。

例えば、画像生成部１４０は、それぞれ異なる複数の到達時間ごとに、到達時間の線画像データ２０が複数生成する。そして、表示制御部１７３は、到達時間が早い順に、複数の線画像データ２０をそれぞれ表示させる。これにより、超音波診断装置１０は、線画像データ２０を動画として表示する。

なお、超音波診断装置１０は、上述した線画像データ２２、変位位置画像データ、及び、線画像データ２４についても同様に、動画として表示することができる。

（信頼性の指標となる独自パラメータの算出）
また、上記の実施形態では、硬さの信頼性の指標として、到達時間、変位の大きさ（最大変位量）、及び到達時間の分散値のうち少なくとも一つを画像化して表示する場合を説明したが、実施形態はこれに限定されるものではない。例えば、上記の硬さの信頼性の指標となるパラメータを複数組み合わせることで、信頼性の指標となる独自のパラメータを算出しても良い。

例えば、算出部１７４は、変位の大きさと分散値とに基づく値を、硬さの信頼性の指標として算出する。ここで、変位の大きさは、例えば、最大変位量であり、この値が大きいほど信頼性が高くなる。また、分散値は、この値が大きいほど、信頼性が低くなる。したがって、例えば、算出部１７４は、最大変位量と、分散値の逆数とをスコア化することで、信頼性の指標となる独自のパラメータを算出する。

そして、線画像生成部１７２及び画像生成部１４０は、到達時間、変位の大きさ、又は到達時間の分散値に代えて、算出された独自のパラメータを用いて、上記の実施形態にて説明した処理を実行する。例えば、線画像生成部１７２は、線画像データ２０に含まれる各位置に、独自のパラメータに応じた画素値を割り当てる。また、例えば、画像生成部１４０は、走査領域内の各位置に対して、独自のパラメータに応じた画素値を割り当てた画像データを背景画像として生成する。また、例えば、画像生成部１４０は、背景画像に含まれる各位置の画素値を、走査領域内の各位置における独自のパラメータに応じて変更する。

なお、これに限らず、算出部１７４は、到達時間、変位の大きさ（最大変位量）、及び到達時間の分散値のうち、任意のパラメータを任意数組み合わせることで、信頼性の指標となる独自のパラメータを算出しても良い。

（並列表示における信頼性表示の連動）
また、上述した並列表示においては、更に、並列表示の対象となる画像同士で信頼性表示を連動させても良い。すなわち、表示制御部１７３は、並列表示において、硬さ画像データに含まれる各位置における信頼性の指標となるパラメータ（例えば、到達時間、変位の大きさ、分散値、独自のパラメータ等）が所定の条件に満たない場合に、その位置の画像を非表示にするとともに、並列表示の対象となる線画像データ２０のうち、その位置に対応する画像を非表示にする。

図２１は、並列表示される場合にモニタに表示される表示画像３０の一例を示す図である。表示画像３０は、重畳画像３１と、重畳画像３２とを含む。重畳画像３１及び重畳画像３２は、ある走査領域の反射波データに基づいて生成されたＢモード画像３３上に、硬さ画像３４又は線画像群３５が重畳された画像である。つまり、硬さ画像３４は、Ｂモード画像３３において、対応する領域の硬さを表す画像である。また、線画像群３５は、線画像３５ａ，３５ｂ，３５ｃ・・・等、複数の線画像を含み、Ｂモード画像３３において、対応する領域の到達時間に基づいて生成された画像である。

ここで、硬さ画像３４のうち、信頼性の低い領域３６は、非表示となる。図２１の例では、領域３６は、黒い線で囲まれた白色の領域で示される。例えば、表示制御部１７３は、硬さ画像３４に含まれる点のうち、分散値が所定の閾値以上となる点の画素値を変更することで、領域３６を非表示とする。なお、ここでは、領域３６を黒い線で囲まれた白色の領域で示したが、これに限定されるものではない。例えば、領域３６は、黒色の領域で示されても良いし、無色の領域で示されても良い。なお、無色の領域で示される場合には、背景であるＢモード画像３３がそのまま表示されることとなる。

また、図２１において、線画像群３５のうち、領域３６に対応する領域３７は、他の領域とは異なる色で表示される。図２１の例では、領域３７は、黒い線で囲まれた白色の領域で示される。例えば、表示制御部１７３は、線画像群３５に含まれる点のうち、分散値が所定の閾値以上となる点の画素値を変更することで、領域３７を他の領域とは異なる色で表示する。

このように、表示制御部１７３は、並列表示を行う場合に、並列表示の対象となる複数の画像それぞれを、同一の信頼性の指標値の閾値を用いて変更する。これにより、表示制御部１７３は、並列表示の対象となる複数の画像それぞれを連動させて非表示とすることが可能となる。

なお、図２１では、領域３６を黒い線で囲まれた白色の領域で示したが、これに限定されるものではない。例えば、領域３６は、黒色の領域で示されても良いし、無色の領域で示されても良い。なお、無色の領域で示される場合には、背景であるＢモード画像３３がそのまま表示されることとなる。また、図２１では、信頼性の指標値として、分散値が用いられる場合を説明したが、実施形態はこれに限定されるものではない。例えば、信頼性の指標となるパラメータとしては、到達時間、変位の大きさ、分散値、及び上記の独自のパラメータのうち、任意のパラメータが用いられれば良い。

なお、上記の実施形態の説明で図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示のように構成されていなくても良い。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行なわれる各処理機能は、その全部又は任意の一部が、ＣＰＵおよび当該ＣＰＵにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現され得る。

また、上記の実施形態で説明した超音波診断装置１０における各処理は、予め用意された画像処理プログラムを実行することによって実現することができる。この画像処理プログラムは、インターネット等のネットワークを介して配布することができる。また、この画像処理プログラムは、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤ等のコンピュータで読み取り可能な非一時的な記録媒体に記録され、コンピュータによって非一時的な記録媒体から読み出されることによって実行することもできる。

以上説明した少なくともひとつの実施形態によれば、生体組織の硬さの信頼性を超音波画像上に表すことができる。

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

１０ 超音波診断装置
１１０ 送信部
１４０ 画像生成部
１７０ 制御部
１７２ 線画像生成部
１７３ 表示制御部

Claims (22)

1. 音響放射力に基づいて生体組織に変位を発生させる変位発生用超音波を超音波プローブから送信させ、前記変位発生用超音波に基づいて発生する所定の走査領域内の生体組織の変位を観測する観測用超音波を前記超音波プローブから送信させる送信部と、
   前記超音波プローブで受信された反射波に基づいて、反射波データを生成する受信部と、
   前記反射波データに基づいて、前記走査領域内の複数位置それぞれにおける変位を複数の時相で算出し、算出された変位が略最大となる時相を前記複数位置それぞれで決定し、前記複数位置のうち決定された時相が略同一の位置同士を表す画像データを生成する画像生成部と、
   前記画像データに基づく画像を、前記走査領域を含む領域に対応する医用画像に重畳表示させる表示制御部と、
   を備える、超音波診断装置。
2. 前記画像生成部は、離散的な複数時相に対応する前記画像データを生成する、請求項１に記載の超音波診断装置。
3. 前記画像データに含まれる位置ごとに、当該位置を含む所定領域内の各位置の前記決定された時相の分散値、及び、当該分散値と前記変位の大きさとに基づく値のうち、少なくとも一方を算出する算出部を更に備え、
   前記画像生成部は、更に、前記画像データに含まれる各位置に、当該位置の前記分散値又は前記値に応じた画素値を割り当てる、請求項１又は２に記載の超音波診断装置。
4. 前記画像生成部は、更に、前記画像データに含まれる各位置に、当該位置のせん断波に基づく生体組織の硬さの指標値、又は、当該位置における前記変位の大きさのうちのいずれかに応じた画素値を割り当てる、請求項１に記載の超音波診断装置。
5. 前記画像生成部は、それぞれ異なる複数の前記決定された時相ごとに前記画像データを複数生成し、当該画像データごとに当該決定された時相に応じた画素値を割り当てる、請求項１に記載の超音波診断装置。
6. 前記表示制御部は、それぞれ異なる複数の前記決定された時相ごとに前記画像データが複数生成された場合に、当該決定された時相が早い順に、当該画像データをそれぞれ表示させる、請求項５に記載の超音波診断装置。
7. 前記画像生成部は、更に、前記走査領域内の各位置のせん断波に基づく生体組織の硬さの指標値が略同一の位置同士を表す第２画像データを生成し、
   前記表示制御部は、更に、前記第２画像データに基づく画像を前記医用画像上に表示させる、請求項１に記載の超音波診断装置。
8. 前記画像生成部は、更に、前記走査領域内の各位置の前記変位の大きさが略同一の位置同士を表す第３画像データを生成し、
   前記表示制御部は、更に、前記第３画像データに基づく画像を前記医用画像上に表示させる、請求項１に記載の超音波診断装置。
9. 前記画像生成部は、更に、所定領域におけるせん断波の伝播速度で、前記せん断波が前記走査領域内の各位置に到達する到達時間を算出し、算出した到達時間が略同一の位置同士を表す第４画像データを生成し、
   前記表示制御部は、更に、前記第４画像データに基づく画像を前記医用画像上に表示させる、請求項１に記載の超音波診断装置。
10. 前記走査領域内の各位置に対してＢモードにおける信号強度に応じた画素値を割り当てた第５画像データ、前記走査領域内の各位置に対してせん断波に基づく生体組織の硬さの指標値に応じた画素値を割り当てた第６画像データ、前記走査領域内の各位置に対して当該位置における前記決定された時相に応じた画素値を割り当てた第７画像データ、前記走査領域内の各位置に対して当該位置における前記変位の大きさに応じた画素値を割り当てた第８画像データ、及び、前記走査領域内の各位置に対して前記決定された時相の分散値に応じた画素値を割り当てた第９画像データ、前記走査領域内の各位置に対して前記決定された時相の分散値と前記変位の大きさとに基づく値に応じた画素値を割り当てた第１０画像データのうち、少なくとも一つの画像データに基づく画像を、前記医用画像として生成する第２画像生成部を更に備える、請求項１に記載の超音波診断装置。
11. 前記第２画像生成部は、前記医用画像に含まれる各位置の画素値を、前記走査領域内の各位置における前記指標値、前記決定された時相、前記変位の大きさ、前記分散値、及び、前記値の各パラメータのうちいずれか一つのパラメータを用いて変更する、請求項１０に記載の超音波診断装置。
12. 前記表示制御部は、前記第５画像データ、前記第６画像データ、前記第７画像データ、前記第８画像データ、前記第９画像データ、及び前記第１０画像データのうち、少なくとも二つの画像データに基づく画像を並列表示させるとともに、並列表示させた画像データのうち、少なくとも一つの画像データに基づく画像上に、前記画像生成部により生成された画像データに基づく画像を重畳表示させる、請求項１０に記載の超音波診断装置。
13. 前記表示制御部は、前記並列表示において、前記第６画像データに含まれる各位置の前記決定された時相、前記変位の大きさ、前記分散値、及び、前記値の各パラメータのうちいずれか一つのパラメータが所定の条件に満たない場合に、当該位置の画像を非表示にするとともに、当該並列表示の対象となる前記画像生成部により生成された画像データのうち、当該位置に対応する画像を非表示にする、請求項１２に記載の超音波診断装置。
14. 前記表示制御部は、関心領域内の各位置におけるせん断波に基づく生体組織の硬さの指標値、前記決定された時相、前記変位の大きさ、前記決定された時相の分散値、及び当該分散値と前記変位の大きさとに基づく値のうち少なくとも一つのパラメータを評価し、評価結果を表示させる、請求項１に記載の超音波診断装置。
15. 反射波データに基づいて、走査領域内の複数位置それぞれにおける変位を複数の時相で算出し、算出された変位が略最大となる時相を前記複数位置それぞれで決定し、前記複数位置のうち決定された時相が略同一の位置同士を表す画像データを生成する画像生成部と、
    前記画像データに基づく画像を、前記走査領域を含む領域に対応する医用画像に重畳表示させる表示制御部と
    を備える、画像処理装置。
16. 反射波データに基づいて、走査領域内の複数位置それぞれにおける変位を複数の時相で算出し、算出された変位が略最大となる時相を前記複数位置それぞれで決定し、前記複数位置のうち決定された時相が略同一の位置同士を表す画像データを生成し、
    前記画像データに基づく画像を、前記走査領域を含む領域に対応する医用画像に重畳表示させる
    ことを含む、画像処理方法。
17. 音響放射力に基づいて生体組織に変位を発生させる変位発生用超音波を超音波プローブから送信させ、前記変位発生用超音波に基づいて発生する所定の走査領域内の生体組織の変位を観測する観測用超音波を前記超音波プローブから送信させる送信部と、
    前記超音波プローブで受信された反射波に基づいて、反射波データを生成する受信部と、
    前記走査領域内の各位置のせん断波に基づく生体組織の硬さの指標値が略同一の位置同士を表す画像データを生成する画像生成部と、
    前記画像データに基づく画像を、前記走査領域を含む領域に対応する医用画像に重畳表示させる表示制御部と
    を備える、超音波診断装置。
18. 音響放射力に基づいて生体組織に変位を発生させる変位発生用超音波を超音波プローブから送信させ、前記変位発生用超音波に基づいて発生する所定の走査領域内の生体組織の変位を観測する観測用超音波を前記超音波プローブから送信させる送信部と、
    前記超音波プローブで受信された反射波に基づいて、反射波データを生成する受信部と、
    前記走査領域内の各位置の前記変位の大きさが略同一の位置同士を表す画像データを生成する画像生成部と、
    前記画像データに基づく画像を、前記走査領域を含む領域に対応する医用画像に重畳表示させる表示制御部と
    を備える、超音波診断装置。
19. 音響放射力に基づいて生体組織に変位を発生させる変位発生用超音波を超音波プローブから送信させ、前記変位発生用超音波に基づいて発生する所定の走査領域内の生体組織の変位を観測する観測用超音波を前記超音波プローブから送信させる送信部と、
    前記超音波プローブで受信された反射波に基づいて、反射波データを生成する受信部と、
    所定領域におけるせん断波の伝播速度で、前記せん断波が前記走査領域内の各位置に到達する到達時間を算出し、算出した到達時間が略同一の位置同士を表す画像データを生成する画像生成部と、
    前記画像データに基づく画像を、前記走査領域を含む領域に対応する医用画像に重畳表示させる表示制御部と
    を備える、超音波診断装置。
20. 音響放射力に基づいて生体組織に変位を発生させる変位発生用超音波を超音波プローブから送信させ、前記変位発生用超音波に基づいて発生する所定の走査領域内の生体組織の変位を観測する観測用超音波を前記超音波プローブから送信させる送信部と、
    前記走査領域内の各位置に対して当該位置における前記変位の大きさに応じた画素値を割り当てた画像データを生成する画像生成部と
    を備えたことを特徴とする超音波診断装置。
21. 音響放射力に基づいて生体組織に変位を発生させる変位発生用超音波を超音波プローブから送信させ、前記変位発生用超音波に基づいて発生する所定の走査領域内の生体組織の変位を観測する観測用超音波を前記超音波プローブから送信させる送信部と、
    前記走査領域内の各位置に対して到達時間の分散値に応じた画素値を割り当てた画像データを生成する画像生成部と
    を備えたことを特徴とする超音波診断装置。
22. 音響放射力に基づいて生体組織に変位を発生させる変位発生用超音波を超音波プローブから送信させ、前記変位発生用超音波に基づいて発生する所定の走査領域内の生体組織の変位を観測する観測用超音波を前記超音波プローブから送信させる送信部と、
    前記超音波プローブで受信された反射波に基づいて、反射波データを生成する受信部と、
    前記反射波データに基づいて、超音波画像データを生成し、当該超音波画像データに含まれる各位置の画素値を、前記走査領域内の各位置におけるせん断波に基づく生体組織の硬さの指標値、到達時間、前記変位の大きさ、及び、前記到達時間の分散値の各パラメータのうちいずれか一つのパラメータを用いて変更する画像生成部と
    を備えたことを特徴とする超音波診断装置。