JP6457105B2 - 音響波診断装置およびその制御方法 - Google Patents

Description

この発明は，音響波たとえば超音波を利用した診断装置およびその制御方法に関する。

断層像や血流を観察するために超音波を利用した様々な診断装置が開発されている（特許文献１〜６）。近年では，生体組織の硬さまたは柔らかさに関する情報を，超音波を用いて計測し，これを医学的診断に利用するエラストグラフィも知られている。エラストグラフィでは，生体組織の硬さまたは柔らかさを表す歪み等の大きさに応じた色相（色）を有する弾性画像が診断に用いられる。

特開２０１４−３６７７８号公報 特開２０１３−５８７６号公報 特開２０１５−２３９１３号公報 特開２０１５−１６０６２号公報 特開２０１４−４２８２３号公報 特開２０１３−１２１４８３号公報

エラストグラフィでは，圧迫，拍動などに起因して生じる生体組織の移動量（変位）を微分することで生体組織の歪みを算出し，生体組織の硬さを推定する。ここで，生体組織が変位を起こせない程度に圧迫，拍動が小さい場合には生体組織の歪みを正しく算出することは難しい。この場合に求められる歪みは，生体組織の歪みではなく，信号のノイズのちらつきなどを歪みとして誤認識した結果である，すなわちノイズに起因するものであることがある。ノイズが歪みとして誤認識されてしまうと，生体組織の硬さに起因しないノイズがあたかも生体組織の硬さを表しているかのように弾性画像に現れてしまう。

上記特許文献のうち，たとえば特許文献１には，心臓壁の移動量に応じて色変換テーブルを変えることで，心拍による生体組織の圧迫と弛緩の度合にかかわらず，同じ弾性を有する部分を弾性画像において大きく異なることがない色で表示させるものが記載されている。しかしながら，特許文献１では，生体組織の圧迫と弛緩の度合が歪みを正しく求めることができない程度に小さい場合について考慮されていないため，生体組織の圧迫と弛緩の度合が小さいためにノイズが歪みとして誤認識されたときにそのノイズがあたかも正常に歪みとして計測できているものとして弾性画像に表示されてしまう。また，たとえば，特許文献２には，低Ｓ／Ｎ領域を信頼性の高い領域と区別して画像表示し，信頼性が高い領域を解析の対象とするものが開示されている。しかしながら，特許文献２では超音波を送信せずに受信処理のみを実行することで得られるＢモード画像から低Ｓ／Ｎ領域を決定しており，生体組織の移動量が小さい場合に歪みとして誤認識されるノイズの除去については何ら考慮されていない。

この発明は，生体組織の移動量が小さい場合にノイズが弾性画像において目立たないようにすることを目的とする。

この発明による音響波診断装置は，被検体の生体組織から反射される音響波エコーを表わす音響波エコー信号を用いて生成された音響波フレームデータを用いて移動量フレームデータを生成し，かつ，移動量フレームデータを用いて移動量の代表値を算出する移動量算出手段，フレーム間の移動量フレームデータを用いて歪みを算出する歪み算出手段，算出された移動量の代表値があらかじめ定められる閾値よりも大きい場合には，第１の色相変換ルックアップテーブルを用いて算出された歪みを表す弾性画像を生成し，算出された移動量の代表値があらかじめ定められる閾値よりも小さい場合には，第１の色相変換ルックアップテーブルよりも色相変化の程度が小さい第２の色相変換ルックアップテーブルを用いて算出された歪みを表す弾性画像を生成する弾性画像生成手段を備えている。

この発明は，音響波診断装置に適した制御方法も提供している。すなわち，この方法は，移動量算出手段が，被検体の生体組織から反射される音響波エコーを表わす音響波エコー信号を用いて生成された音響波フレームデータを用いて移動量フレームデータを生成し，かつ，移動量フレームデータを用いて移動量の代表値を算出し，歪み算出手段が，フレーム間の移動量フレームデータを用いて歪みを算出し，弾性画像生成手段が，算出された移動量の代表値があらかじめ定められる閾値よりも大きい場合には，第１の色相変換ルックアップテーブルを用いて算出された歪みを表す弾性画像を生成し，算出された移動量の代表値があらかじめ定められる閾値よりも小さい場合には，第１の色相変換ルックアップテーブルよりも色相変化の程度が小さい第２の色相変換ルックアップテーブルを用いて算出された歪みを表す弾性画像を生成する。

算出された移動量の代表値が閾値に一致する場合は，第１または第２の色相変換ルックアップテーブルのいずれかを，弾性画像の生成に用いることができる。

好ましい実施態様では，移動量算出手段は，音響波エコー信号から生成される取得時刻の異なる一対の音響波フレームデータの単位領域ごとの移動量の平均値を，移動量の代表値として算出するものである。

一実施態様では，音響波エコー信号から経時的に生成される音響波フレームデータのフレームレートに応じて閾値を変動させる閾値制御手段を備えている。

弾性画像において色相によって表される歪み（生体組織の硬さを表す指標）は，生体組織の移動量（変位）から算出される。移動量が小さい場合，移動量から算出される歪みはノイズの影響を受けやすくなる。この発明によると，移動量の代表値があらかじめ定められる閾値よりも小さい場合に，移動量の代表値が閾値よりも大きい場合に用いられる第１の色相変換ルックアップテーブルに比べて，色相変化の程度が小さい第２の色相変換ルックアップテーブルが用いられて弾性画像が生成される。したがって，移動量が小さいときに，ノイズが色相によって弾性画像中に顕著に現れてしまうことを防止するまたは低減することができる。

超音波診断装置の全体構成を示すブロック図である。 弾性画像の作成処理の流れを示すフローチャートである。 色相変換ルックアップテーブルを示す。 病変組織が存在する生体内部を模式的に示す。 圧迫前の周辺組織および病変組織の位置と，組織移動量が大きいときの圧迫後の周辺組織および病変組織の位置と，組織移動量が小さいときの圧迫後の周辺組織および病変組織の位置とを，それぞれ示す。 圧迫前後の周辺組織および病変組織の変位量を示すグラフである。 周辺組織および病変組織の歪みを示すグラフである。 組織移動量が大きいときおよび小さいときのそれぞれについて，従来技術における弾性画像と本願の超音波診断装置における弾性画像を，それぞれ示す。 色相変換ルックアップテーブルの他の例を示す。

以下に説明する実施例では，音響波として超音波を用いる例を説明する。超音波に限定されることはなく，被験対象，測定条件などに応じて適切な周波数が選択されれば，可聴周波数の音響波を用いるようにしてもよい。

図１は，超音波診断装置１の全体構成を示すブロック図を示している。超音波診断装置１は，超音波プローブ２，送受信ビームフォーマ３，エコーデータ処理装置４，画像制御装置５，表示装置６，操作装置７，制御装置８および記憶装置９を備えている。

制御装置８によって超音波診断装置１の動作が全体的に統括される。制御装置８に接続される記憶装置９に，超音波診断装置１を構成する以下に詳述する各種装置を制御する制御プログラム，後述する色相変換ルックアップテーブル，フレームデータなどが記憶される。操作者の指示，設定ないし調整される値などが，操作装置７から入力される。

超音波プローブ２が被検体（患者）の体表面に押し当てられる。超音波プローブ２は，被検体に向けて超音波ビームを送信し，かつ被検体内の生体組織から反射される超音波エコーを受信し，超音波エコーを表す超音波エコー信号を出力する。コンベックス型，セクタ型，リニア型など，任意の形状の超音波プローブ２を用いることができる。

送受信ビームフォーマ３は，超音波プローブ２を所定の走査条件で駆動して超音波ビームを走査する。セクタ走査，オフセット・セクタ走査，リニア走査など任意の走査方式を採用することができる。また送受信ビームフォーマ３は，超音波プローブ２からの超音波エコー信号に対し整相加算処理等の所定の信号処理を行ない，生体組織の一走査面（一断層）に対応する超音波フレームデータ（断層エコーデータ）を生成する。生成された超音波フレームデータは順次記憶装置９に記憶される。超音波フレームデータは所定のフレームレート（超音波フレームデータの単位時間あたりの生成数）によって経時的に生成される。超音波プローブ２によって送信される超音波ビームの送信タイミングを変更することで，超音波フレームデータのフレームレートを変更することができる。

超音波フレームデータは，次にエコーデータ処理装置４に入力する。エコーデータ処理装置４は，Ｂモードデータ作成部４Ａ，移動量計測部４Ｂおよび歪み演算部４Ｃを備える。

Ｂモードデータ作成部４Ａは，超音波フレームデータに対し，対数圧縮処理，包絡線検波処理等を行い，Ｂモードデータを作成する。

移動量計測部４Ｂ（移動量算出手段）は，超音波フレームデータを用いて以下の移動量を算出する。すなわち，順次記憶装置９に記憶される取得時刻の異なる複数の超音波フレームデータのうちの所定時間間隔の一対のフレームデータに基づいて単位領域ごと（ピクセルごと）の生体組織の移動量を算出し，移動量フレームデータを作成する。移動量の大きさはピクセル（pix）を単位とすることができる。その後，移動量フレームデータにおける単位領域ごとの移動量の平均値を算出することで組織移動量の代表値を算出する。算出される組織移動量の代表値は，後述するように，色相変換ルックアップテーブルの設定（選択，切り換え）に用いられる。フレーム全体ではなくフレーム中心の単位領域を用いて組織移動量の代表値を算出してもよく，また平均値に代えて，最大値，最小値，分散などを組織移動量の代表値として用いてもよい。さらに，超音波画像の一部を関心領域（Region Of Interest，ＲＯＩ）として指定した場合には，ＲＯＩ中の単位領域を用いて組織移動量の代表値を算出してもよい。

歪み演算部４Ｃ（歪み算出手段）は，上述した移動量計測部４Ｂによって作成された移動量フレームデータに基づいて単位領域ごと（ピクセルごと）の歪みを求め，歪みフレームデータを作成する。歪みフレームデータは，移動量フレームデータを微分することによって作成（演算）される。

エコーデータ処理装置４のＢモードデータ作成部４Ａから出力されるＢモードデータおよび歪み演算部４Ｃから出力される歪みフレームデータは画像制御装置５に入力する。画像制御装置５は，Ｂモード画像データ作成部５Ａ，弾性画像データ作成部５Ｂおよび画像表示制御部５Ｃを備えている。

Ｂモード画像データ作成部５ＡはＢモードデータに対してスキャンコンバータによる走査変換を行い，表示装置６における表示に適する２次元の断層像データ（Ｂモード画像データ）を作成する。Ｂモード画像データでは信号強度が輝度によって表される。Ｂモード画像データはたとえばピクセルごとに２５６階調の輝度を表す情報を持つ。

弾性画像データ作成部５Ｂは，歪みフレームデータにおける単位領域ごとの歪みを，その大きさに応じた色相（色の違い）によって表すカラーの弾性画像データを作成する，すなわち歪みを色相化する処理を実行する。弾性画像データによって表される弾性画像では，歪みの大きさが色相の違いによって表される。歪みの大きさと色相との対応関係は，記憶装置９にあらかじめ記憶される色相変換ルックアップテーブルに基づく。色相変換ルックアップテーブルの詳細は後述する。

画像表示制御部５Ｃは，Ｂモード画像データと弾性画像データとを合成し，合成画像データを作成する。合成画像データが表示装置６に与えられることで，Ｂモード画像と弾性画像とが合成された合成超音波画像が表示装置６の表示画面に表示される。もちろん，Ｂモード画像と弾性画像と合成せずに，表示画面上に並べて表示してもよい。

図２は超音波診断装置１における弾性画像の作成処理を示すフローチャートである。

被検体に超音波プローブ２を軽く押し付ける，または心臓の拍動等を利用することで，生体組織を物理的に変位させつつ，断層（一走査面）に対応する超音波フレームデータを連続的に取得する（超音波の計測）（ステップ11）。超音波フレームデータは所定のフレームレートで取得されて記憶装置９に順次記録される。

記憶装置９に順次記録される超音波フレームデータのうちの一対の超音波フレームデータを用いて，エコーデータ処理装置４において超音波フレームデータの単位領域ごとの組織の移動量が計測されて移動量フレームデータが作成される。さらに上述した組織移動量の代表値が算出される（ステップ12）。

次に，算出された組織移動量の代表値の大きさに応じて，色相変換ルックアップテーブル（以下，色相変換ＬＵＴという）が設定される（ステップ13，14，15）。

図３は記憶装置９に記憶される色相変換ＬＵＴを，横軸を歪みσ，縦軸を色情報度数Ｎとするグラフによって示している。色情報度数Ｎはたとえば０〜２５５の間の値をとる。

色相変換ＬＵＴは，歪みを，その大きさに応じた色相に変換するために用いられる。歪みが大きい生体組織（柔らかい生体組織）には一般に赤色が対応付けられる。逆に歪みが小さい生体組織（硬い生体組織）には一般に青色が対応づけられる。中間の大きさの歪みには緑色が対応づけられる。これにより生体組織の硬さが色相によって分かりやすく可視化されることになる。

記憶装置９には，色相変換ＬＵＴ_largeと，色相変換ＬＵＴ_smallの２つの色相変換ＬＵＴが記憶されている。図３を参照して，色相変換ＬＵＴ_largeのグラフ（実線）の勾配（歪みの変化量に対する色相変化の程度の割合，グラフの傾き）は大きく，他方，色相変換ＬＵＴ_smallのグラフ（一点鎖線）の勾配は非常に小さい。勾配の大きい色相変換ＬＵＴ_largeは歪みの変化量に対する色相変化の程度が大きく，色表現範囲が広いことを意味する。他方，勾配の小さい色相変換ＬＵＴ_smallは歪みの変化量に対する色相変化の程度が小さく，色表現範囲が狭いことを意味する。色相変換ＬＵＴ_largeのグラフと色相変換ＬＵＴ_smallのグラフはＮ／２（緑色）において交差している。

図２に戻って，エコーデータ処理装置４において算出される組織移動量の代表値の大きさに応じて，色相変換ＬＵＴ_largeと色相変換ＬＵＴ_smallのいずれかが選択される。制御装置８において組織移動量の代表値と所定の閾値とが比較される（ステップ13）。代表値が所定の閾値以上の場合には色相変換ＬＵＴ_largeが設定される（ステップ13でYES ，ステップ14）。すなわち，色相変換ＬＵＴ_largeが記憶装置９から読み出されて画像制御装置５に与えられる。他方，代表値が所定の閾値未満の場合には，色相変換ＬＵＴ_smallが設定される（ステップ13でNO，ステップ15）。もちろん，代表値が所定の閾値を超える場合に色相変換ＬＵＴ_largeを設定し，かつ代表値が所定の閾値以下の場合に色相変換ＬＵＴ_smallを設定してもよい。

組織移動量の代表値が閾値以上である場合，その組織移動量は生体組織の硬さに起因していると考えられ，信頼性は高い。この場合には，従来通りの勾配を持つ色相変換ＬＵＴ_largeを用いることで，従来と変わらない弾性画像（エラストグラフィ画像）を表現することができる。他方，組織移動量の代表値が閾値よりも小さい場合，その組織移動量は生体組織の硬さに起因するのではなく，ノイズによって支配されている可能性がある。この場合に勾配の小さい色相変換ＬＵＴ_smallを用いることで，ノイズが目立たない弾性画像とすることができる（詳細は後述する）。

エコーデータ処理装置４において移動量フレームデータを微分することによって単位領域ごとの歪みが算出され，歪みフレームデータが作成される（ステップ16）。次に，上述した色相変換ＬＵＴ_largeまたは色相変換ＬＵＴ_smallを用いて歪みフレームデータを色相変換して，色相化（カラー化）することで弾性画像データが作成される（ステップ17）。歪みをその大きさに応じた色相によって表わす弾性画像が，画像制御装置５の制御の下，表示装置６に表示される（ステップ18）。

図４から図８を参照して，組織移動量に応じて色相変換ＬＵＴ_largeまたは色相変換ＬＵＴ_smallのいずれかを選択する効果を説明する。図４は柔らかい周辺組織（正常な生体組織）に囲まれるようにして硬い病変組織が存在する生体内部を模式的に示している。図５は，図４のＶ−Ｖ線に沿う，すなわち病変組織を通る断面位置における周辺組織および病変組織を模式的に示すもので，符号21は圧迫前（弛緩時）の周辺組織および病変組織を，符号22および23はいずれも圧迫後（圧縮時）の周辺組織および病変組織を示している。ここで符号22は組織移動量が大きいときを，符号23は組織移動量が小さいときを，それぞれ模式的に示している。分かりやすくするために，図５の符号21〜23の周辺組織および病変組織の模式図において，周辺組織の左端は固定されている（変位がない）ものとして示されている。

図６は，図５の符号21と符号22，符号21と符号23のそれぞれに示す圧迫前後の周辺組織および病変組織の変位量をグラフによって示している。図７は図６に示す変位を微分したものである。図６および図７において，組織移動量が大きい（代表値が閾値以上）ときのグラフ（図５の符号21と符号22の関係）が太実線によって示され，組織移動量が小さい（代表値が閾値未満）ときのグラフ（図５の符号21と符号23の関係）が細実線によって示されている。

図８は，組織移動量が大きいときおよび組織移動量が小さいときのそれぞれについて，色相変換ＬＵＴ_largeのみが用いられる従来の弾性画像と，組織移動量が大きいときに色相変換ＬＵＴ_largeが，組織移動量が小さいときに色相変換ＬＵＴ_smallがそれぞれ設定される，本願の超音波診断装置１において表示される弾性画像とを並べて示している。図８に示す弾性画像において，青色で表示される組織範囲を「Ｂ」により，赤色で表示される組織範囲を「Ｒ」により，緑色で表示される組織範囲を「Ｇ」により，それぞれ模式的に示す。

まず，組織移動量が大きい場合，すなわち組織移動量の代表値が閾値以上である場合を説明する。図５の符号21と符号22の関係，および図６の太実線のグラフを参照して，組織移動量が大きい場合には，柔らかい周辺組織の変位量と硬い病変組織の変位量に大きな相違が生じる（グラフの傾きが全く異なる）。図７の太実線を参照して，組織移動量が大きい場合，変位を微分（深さ方向に微分）することで算出される歪みは，柔らかい周辺組織について大きな値が算出され，硬い病変組織について小さい値が算出される。図８を参照して，組織移動量が大きい場合には，従来例においても，本願においても，色相変換ＬＵＴ_largeが用いられる。弾性画像では，病変組織については硬い（歪みが小さい）ことを表す青色Ｂが示される。周辺組織については病変組織よりも柔らかいことを表す色（ここでは緑色Ｇ）によって表現される。硬い生体組織（病変組織）の存在を視覚的に認識することができる。

次に組織移動量が小さい，すなわち組織移動量の代表値が閾値未満である場合を説明する。超音波プローブによる圧迫によって生体組織を変位させる場合には，超音波プローブによる圧迫が非常に小さい場合や，被検体の計測箇所が深い（超音波プローブからの距離が遠い）場合に組織移動量は小さくなる。拍動によって生体組織を変位させる場合であれば，被検体の計測箇所が心臓や大動脈から離れており拍動が伝わりにくい場合に組織移動量が小さくなる。

図５の符号21と符号23の関係，および図６の細線を参照して，組織移動量が小さい場合，周辺組織の変位量と病変組織の変位量との間に変位の差が認められなくなり，さらにノイズ信号が現れることがある。これを微分すると，図７の細線のグラフを参照して，病変組織の硬さと周辺組織の硬さの違いに起因しない値（ノイズ）が大きな歪みとして算出されてしまう。

図８を参照して，組織移動量の大きさにかかわらずに同じ色相変換ＬＵＴ_largeが用いられる従来技術では，組織移動量が小さいとき，ノイズに起因する，生体組織の硬さを表さない色相（ここでは赤色（Ｒ）：柔らかいことを示す色）が弾性画像に現れてしまう。他方，組織移動量の代表値が閾値よりも小さい場合に色相変換ＬＵＴ_smallを用いる本願においては，色相変換ＬＵＴ_smallの色表現範囲が狭いので（図３），生体組織の硬さに起因しないノイズに基づく色相が弾性画像中に顕著に表示されることはない。このように，組織移動量が小さいときに，ノイズが目立つ弾性画像が表示されないようにすることができる。

上述した閾値はたとえば0.2〜0.5ピクセルの間で設定するのが望ましい。この値は，ノイズのちらつきを変位量として計算したときの変動範囲である。ノイズのちらつきは，使用するシステム，プローブなどの性能によって変わるので，使用する条件ごとに，ノイズのちらつきがどの程度の変位量に相当するのかを調べて，適切な閾値を設定するのが望ましい。

図９は記憶装置９に記憶される色相変換ＬＵＴの他の例を示している。色相変換ＬＵＴ_largeおよび色相変換ＬＵＴ_smallに加えて，中間の色表現範囲を持つ色相変換ＬＵＴ_middleが記憶されている。

たとえば，算出される組織移動量の代表値が0.5ピクセル未満であれば，最もグラフの傾きの小さい色相変換ＬＵＴ_small（一点鎖線）が，0.5ピクセル以上１ピクセル未満であれば中間のグラフの傾きを持つ色相変換ＬＵＴ_middle（二点鎖線）が，１ピクセル以上であれば色相変換ＬＵＴ_large（実線）が，それぞれ設定される。

４つ以上の色表現範囲の異なる色相変換ＬＵＴを記憶装置９に記憶しておき，算出される組織移動量の代表値の大きさに応じて設定する色相変換ＬＵＴを細かく切り換えるようにしてもよい。もちろん，色相変換ＬＵＴを組織移動量の代表値の関数として規定して色相変換ＬＵＴをさらに細かく設定してもよい。また，色相はグレースケールによる色表示としてもよく，この場合には白から黒までの階調変化によって弾性画像が表示される。

色相変換ＬＵＴを切り換えるための上述の閾値は，超音波診断装置１のフレームレートに応じて変更できるようにしてもよい。たとえばフレームレートが２倍になると，一対のフレームによって算出される組織移動量は半分になる。あるフレームレートにおける閾値をＡとすると，フレームレートが２倍にされたときに閾値をＡ／２にすれば，フレームレートの変更が色相変換ＬＵＴの切り換えに影響を及ぼさないようにすることができる。たとえば，フレームレートが20fr／sec以上の場合に閾値を0.1〜0.5ピクセルの間に設定し，フレームレートが20fr／sec未満の場合に閾値を0.2〜1.0ピクセルの間に設定してもよい。フレームレートの変更は操作装置７から入力され，フレームレートに応じた閾値の変更は制御装置８によって実行することができる。

１ 超音波診断装置（音響波診断装置）
２ 超音波プローブ
３ 送受信ビームフォーマ
４ エコーデータ処理装置（移動量算出手段，歪み算出手段）
５ 画像制御装置（弾性画像生成手段）
６ 表示装置
７ 操作装置
８ 制御装置（閾値制御手段）
９ 記憶装置

Claims (4)

1. 被検体の生体組織から反射される音響波エコーを表わす音響波エコー信号を用いて生成された音響波フレームデータを用いて移動量フレームデータを生成し，かつ，該移動量フレームデータを用いて移動量の代表値を算出する移動量算出手段，
   フレーム間の移動量フレームデータを用いて歪みを算出する歪み算出手段，および
   上記算出された移動量の代表値があらかじめ定められる閾値よりも大きい場合には，第１の色相変換ルックアップテーブルを用いて上記算出された歪みを表す弾性画像を生成し，上記算出された移動量の代表値が上記あらかじめ定められる閾値よりも小さい場合には，上記第１の色相変換ルックアップテーブルよりも色相変化の程度が小さい第２の色相変換ルックアップテーブルを用いて上記算出された歪みを表す弾性画像を生成する弾性画像生成手段，
   を備える音響波診断装置。
2. 上記移動量算出手段は，上記音響波エコー信号から生成される取得時刻の異なる一対の音響波フレームデータの単位領域ごとの移動量の平均値を，移動量の代表値として算出するものである，
   請求項１に記載の音響波診断装置。
3. 上記音響波エコー信号から経時的に生成される音響波フレームデータのフレームレートに応じて上記閾値を変動させる閾値制御手段を備えている，
   請求項１または２に記載の音響波診断装置。
4. 移動量算出手段が，被検体の生体組織から反射される音響波エコーを表わす音響波エコー信号を用いて生成された音響波フレームデータを用いて移動量フレームデータを生成し，かつ，移動量フレームデータを用いて移動量の代表値を算出し，
   歪み算出手段が，フレーム間の移動量フレームデータを用いて歪みを算出し，
   弾性画像生成手段が，算出された移動量の代表値があらかじめ定められる閾値よりも大きい場合には，第１の色相変換ルックアップテーブルを用いて算出された歪みを表す弾性画像を生成し，算出された移動量の代表値が上記あらかじめ定められる閾値よりも小さい場合には，第１の色相変換ルックアップテーブルよりも色相変化の程度が小さい第２の色相変換ルックアップテーブルを用いて算出された歪みを表す弾性画像を生成する，
   音響波診断装置の制御方法。

Images