JP5710566B2 - 超音波診断装置及びその制御プログラム - Google Patents

Description

本発明は、生体組織に対する超音波の送受信により得られたエコー信号に基づく画像を表示する超音波診断装置及びその制御プログラムに関する。

通常のＢモード画像と、生体組織の硬さ又は軟らかさを表す弾性画像とを合成して表示させる超音波診断装置がある。この種の超音波診断装置において、弾性画像は例えば次のようにして作成される。先ず、生体組織に対し、例えば超音波プローブによる圧迫とその弛緩を繰り返すなどして生体組織を変形させながら超音波の送受信を行なってエコーを取得する。そして、得られたエコーデータに基づいて、相関演算を行なって生体組織の弾性に関する物理量を算出し、この物理量を色情報に変換してカラーの弾性画像を作成する。ちなみに、生体組織の弾性に関する物理量としては、例えば生体組織の歪みなどを算出している。歪みの算出方法は、例えば特許文献１に開示されている。

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

ところで、生体組織にシスト（Ｃｙｓｔ、嚢胞）と呼ばれる液体部分が存在する場合がある。このシストの全体又はその一部分について、前記特許文献１の手法で算出される歪みの値は、生体組織が硬いことを示す値になる場合がある。ここで、弾性画像において、シストと同様に塊りとして表示される腫瘤についても歪みの値は硬いことを示す値になる。このため、弾性画像において腫瘤とシストとを区別することが困難な場合があった。

本願発明者は、シストにおいては、相関演算を行なって算出された生体組織の弾性に関する物理量として、符号がばらばらな値が得られ、その絶対値としては、シスト以外の部分よりも生体組織が軟らかいことを示す値が得られることがあり、かつ相関演算における相関係数が小さく、なおかつエコー信号の強度が小さいというパラメータ特性を有することに着目した。本願発明者は、これら物理量の絶対値、相関係数及びエコー信号の強度の三つのパラメータのうち、少なくとも二つのパラメータにおいてシストに見られるパラメータ特性を反映した画像を表示させる超音波診断装置を課題解決手段として提案する。具体的には、生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって前記生体組織における各部の弾性に関する物理量を算出する物理量算出部と、 この物理量算出部によって算出された物理量の絶対値、前記相関演算における相関係数及び前記生体組織から得られた超音波のエコー信号の強度の三つのパラメータのうち、少なくとも二つのパラメータを用いた演算式であって、前記生体組織におけるシストのパラメータ特性が強調されてシスト以外とは区別できる演算値が得られる演算式の演算値を演算する演算部と、この演算部の演算値に応じた表示形態を有する画像が表示される表示部と、を備えることを特徴とする超音波診断装置である。

上記観点の発明によれば、前記物理量算出部によって算出された物理量の絶対値、前記相関演算における相関係数及び前記生体組織から得られた超音波のエコー信号の強度の三つのパラメータのうち、少なくとも二つのパラメータ特性を用いた演算式であって、前記生体組織におけるシストのパラメータ特性が強調された演算値が得られる演算式の演算値を演算する演算部を備え、この演算部の演算値に応じた表示形態を有する画像が表示されるので、腫瘤とシストとを判別可能な画像を表示することができる。

本発明の第一実施形態の超音波診断装置の概略構成を示すブロック図である。 図１に示す第一実施形態の超音波診断装置における表示制御部の構成を示すブロック図である。 Ｂモード画像とカラー画像とが合成された合成超音波画像が表示された表示部を示す図である。 本発明の第二実施形態の超音波診断装置の概略構成を示すブロック図である。 図４に示す第二実施形態の超音波診断装置における表示制御部の構成を示すブロック図である。 Ｂモード画像と弾性画像とが合成された合成超音波画像が表示された表示部を示す図である。 Ｂモード画像とカラー画像とが合成された合成超音波画像と、Ｂモード画像と弾性画像とが合成された合成超音波画像とが表示された表示部を示す図である。

以下、本発明の実施形態について図面に基づいて詳細に説明する。
（第一実施形態）
先ず、第一実施形態について図１〜図３に基づいて説明する。図１に示す超音波診断装置１は、超音波プローブ２、送受信ビームフォーマ３、Ｂモードデータ作成部４、物理量データ作成部５、演算部６、表示制御部７、表示部８、操作部９、制御部１０及び記憶部１１を備える。

前記超音波プローブ２は、生体組織に対して超音波を送信しそのエコーを受信する。この超音波プローブ２を生体組織の表面に当接させた状態で圧迫と弛緩を繰り返して、生体組織を変形させながら超音波の送受信を行なって取得されたエコーデータに基づいて、後述のようにＢモード画像及びカラー画像が作成される。

前記送受信ビームフォーマ３は、前記制御部１０からの制御信号に基づいて前記超音波プローブ２を所定の走査条件で駆動させて音線毎の超音波の走査を行なう。また、送受信ビームフォーマ３は、前記超音波プローブ２で受信したエコーについて、整相加算処理等の信号処理を行なう。前記送受信ビームフォーマ３で信号処理されたエコーデータは、前記Ｂモードデータ作成部４及び前記物理量データ作成部５に出力される。

前記Ｂモードデータ作成部４は、前記送受信ビームフォーマ３から出力されたエコーデータに対し、対数圧縮処理、包絡線検波処理等のＢモード処理を行い、Ｂモードデータを作成する。Ｂモードデータは、前記Ｂモードデータ作成部４から前記表示制御部７へ出力される。

前記物理量データ作成部５は、前記送受信ビームフォーマ３から出力されたエコーデータに基づいて、生体組織における各部の弾性に関する物理量のデータ（物理量データ）を作成する（物理量算出機能）。前記物理量データ作成部５は、例えば上記特許文献１に記載されているように、一の走査面における同一音線上の時間的に異なるエコーデータに相関ウィンドウを設定し、この相関ウィンドウ間で複素相互相関の虚数部を演算して前記弾性に関する物理量を算出し前記物理量データを作成する。より詳細には、前記物理量データ作成部は、前記超音波プローブ２による生体組織の変形の前後（圧迫前と圧迫後、弛緩前と弛緩後）におけるエコー信号の波形の圧縮率を求める。ここでの波形の圧縮率は、相関ウィンドウ間における波形の圧縮率である。この結果、前記弾性に関する物理量として、歪みが得られる。前記物理量データ作成部５は、音線方向において複数の相関ウィンドウを設定して各相関ウィンドウについて歪みを算出する。従って、歪みは、一音線上における複数点（相関ウィンドウに対応する部分）について算出される。

前記物理量データ作成部５は、本発明における物理量算出部の実施の形態の一例であり、また前記物理量算出機能は本発明における物理量算出機能の実施の形態の一例である。

ちなみに、算出される歪みは、生態組織の変位の方向に応じた符号を伴っている。前記超音波プローブによって圧迫を行なう時とその弛緩を行なう時とでは、符号は反対になるが、ある時間において、各部分（各点）において得られる演算結果は、同じ符号になるはずである。例えば、圧迫時においては各部において正の符号になり、弛緩時においては各部において負の符号になるはずである。しかし、シストにおいては、上述のようにエコー信号の信号強度が非常に小さい。これにより、Ｓ／Ｎ（ｓｉｇｎａｌ ｔｏ ｎｏｉｓｅ ｒａｔｉｏ）が悪いので、各部分において得られる演算結果の符号が同じにならない。

前記演算部６は、前記物理量データ作成部５で得られた歪み、前記相関演算における相関係数及び生体組織から得られた超音波のエコー信号の強度の三つのパラメータのうち、少なくとも二つのパラメータを用いた演算式の演算値を算出する（演算機能）。詳細は後述する。前記演算部６は、本発明における演算部の実施の形態の一例である。また、前記演算機能は、本発明における演算機能の実施の形態の一例である。

前記表示制御部７には、前記Ｂモードデータ作成部４からのＢモードデータ及び前記演算部６で得られた演算値のデータが入力されるようになっている。前記表示制御部７は、図２に示すように、Ｂモード画像データ作成部７１、カラー画像データ作成部７２、合成画像表示制御部７３を有している。

前記Ｂモード画像データ作成部７１は、前記Ｂモードデータ作成部４によって得られたＢモードデータを、エコーの信号強度に応じた輝度に対応する情報を有するＢモード画像データに変換する。また、前記カラー画像データ作成部７２は、前記演算部６によって得られた演算値のデータを、色に対応する情報に変換して、前記演算値に応じた色に対応する情報を有するカラー画像データを作成する。前記カラー画像データ作成部７２は、前記演算値のデータを階調化し、各階調に割り当てられた色に対応する情報からなるカラー画像データを作成する。階調数は、例えば２５６である。

前記合成画像表示制御部７３は、前記Ｂモード画像データ及び前記カラー画像データを合成し、前記表示部８に表示する合成超音波画像の画像データを作成する。また、前記合成画像表示制御部７３は、前記合成超音波画像の画像データを、図３に示すように、Ｂモード画像ＢＩとカラー画像ＣＩとが合成された合成超音波画像ＵＩとして前記表示部８に表示させる。前記カラー画像ＣＩは、前記Ｂモード画像ＢＩに設定された領域Ｒ内に表示される。前記カラー画像ＣＩは、半透明で（背景のＢモード画像が透けた状態で）表示される。

前記表示部８は、例えばＬＣＤ（Ｌｉｑｕｉｄ Ｃｒｙｓｔａｌ Ｄｉｓｐｌａｙ）やＣＲＴ（Ｃａｔｈｏｄｅ Ｒａｙ Ｔｕｂｅ）などで構成される。前記操作部９は、操作者が指示や情報を入力するためのキーボード及びポインティングデバイス（図示省略）などを含んで構成されている。

前記制御部１０は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）を有して構成され、前記記憶部１１に記憶された制御プログラムを読み出し、前記物理量算出機能、前記演算機能などをはじめとする前記超音波診断装置１の各部における機能を実行させる。

さて、本例の超音波診断装置１の作用について説明する。本例の超音波診断装置１は、生体組織において本例の抽出対象であるシストが抽出された画像を表示する。具体的に説明する。前記超音波プローブ２を被検体の体表面に当接させた状態で、生体組織に対する超音波の送受信を行なうことにより、エコー信号が取得される。超音波の送受信が行われる時には、前記超音波プローブ２によって生体組織に対する圧迫とその弛緩が行なわれる。

超音波の送受信は二種類であり、Ｂモード画像用の超音波の送受信と、前記物理量データ作成部５が歪みを算出するための超音波の送受信とが、それぞれの送受信パラメータ（ｐａｒａｍｅｔｅｒ）において行われる。

エコー信号が取得されると、前記Ｂモードデータ作成部４が前記Ｂモードデータを作成する。また、前記物理量データ作成部５が相関演算を行なって生体組織の歪みを算出し、前記物理量データを作成する。前記物理量データ作成部５は、一音線上の複数点について相関演算を行ない、各点について歪みを算出する。

前記物理量データが得られると、前記演算部６は、下記（式１）を用いて演算を行なう。
Ｆ（ｎ）＝Ｐ（ｎ）×Ｑ（ｎ） ・・・（式１）

上記（式１）において、ｎは歪みが算出される音線上の点（相関ウィンドウに対応する部分）を示す１〜Ｎまでの自然数である。また、Ｐ（ｎ）は歪みに関する関数であり、Ｑ（ｎ）は相関演算における相関係数に関する関数である。具体的には、Ｐ（ｎ）、Ｑ（ｎ）は下記（式２）及び（式３）である。
Ｐ（ｎ）＝ｋ１×Ａｂｓ（ｓｔｒａｉｎ（ｎ）） ・・・（式２）
Ｑ（ｎ）＝ｋ２×（１／ｘＣｏｒｒ（ｎ）） ・・・（式３）

これら（式２）及び（式３）において、ｋ１，ｋ２は任意に設定される重み係数である。これらｋ１，ｋ２はデフォルトで設定されていてもよいし、操作者が任意に設定できるようになっていてもよい。また、ｓｔｒａｉｎ（ｎ）は、前記音線上の点１〜Ｎにおける歪みの値であり、Ａｂｓ（ｓｔｒａｉｎ（ｎ））は、歪みの絶対値を返す関数である。さらに、ｘＣｏｒｒ（ｎ）は、前記音線上の点１〜Ｎにおける相関係数の値である。

Ｆ（ｎ）は、シストを抽出できる関数になっている。これについて具体的に説明する。例えば、圧迫時において、歪みとして、シストにおいては、「−０．１％」という値が得られ、腫瘤においては、「＋０．０１％」が得られ、それ以外の部分においては、「＋０．１％」、「＋０．３％」という値が得られたとする。また、弛緩時において、歪みとして逆の符号の値、すなわちシストにおいては、「＋０．１％」が得られ、腫瘤においては、「−０．０１％」が得られ、それ以外の部分においては、「−０．１％」、「−０．３％」という値が得られたとする。

上述のような歪みの値が得られた場合、シストの歪みは他の部分と比較して硬いことを示す値になっているため、従来のように歪みの値に基づいて弾性画像を作成した場合、シストは硬いことを示す色で表示される。しかし、前記（式２）においては、歪みの絶対値が得られる。従って、Ｐ（ｎ）の値は腫瘤よりもシストの方が大きくなり、腫瘤よりも軟らかいことを示す値になる。また、シストにおいては、他の部分と比べて相関係数が小さい。従って、前記（式３）によって得られるＱ（ｎ）の値も大きくなる。以上のことから、シストにおいては、Ｆ（ｎ）の値が大きくなり、Ｆ（ｎ）は、シストに見られるパラメータ特性が強調されてシスト以外とは区別できる演算値が得られる関数になっている。

前記演算部６によって得られたＦ（ｎ）の演算値は、前記表示制御部７へ入力される。そして、前記カラー画像データ作成部７２は、Ｆ（ｎ）の演算値に応じた色に対応する情報を有するカラー画像データを作成する。

また、前記Ｂモード画像データ作成部７１は、前記Ｂモードデータ作成部４で作成された前記Ｂモードデータに基づいてＢモード画像データを作成する。

前記Ｂモード画像データ及び前記カラー画像データが作成されると、前記合成画像表示制御部７３は、これらＢモード画像データ及びカラー画像データに基づいて合成超音波画像の画像データを作成し、図３に示すように、前記表示部８に合成超音波画像ＵＩを表示させる。

前記合成超音波画像は、Ｂモード画像ＢＩとカラー画像ＣＩとが合成された画像である。前記カラー画像ＣＩは、Ｆ（ｎ）の演算値に応じた色を有する画像であり、青色、緑色、赤色などの色を有する。例えば、Ｆ（ｎ）の演算値が小さい部分は、前記カラー画像ＣＩにおいて青色で表示され、Ｆ（ｎ）の演算値が大きい部分は、前記カラー画像ＣＩにおいて赤色で表示される。この場合、シストにおいてはＦ（ｎ）の演算値が大きいので、前記カラー画像ＣＩにおいて赤色で表示される。従って、本例の超音波診断装置１によれば、シストが抽出された画像を表示することができる。

一方、腫瘤においては、歪みの絶対値が小さく、なおかつ相関係数は必ずしも小さくはならない。このため、腫瘤においては、シストよりもＦ（ｎ）の演算値が小さくなるので、前記カラー画像ＣＩにおいて、腫瘤はシストとは異なる色で表示される。したがって、前記カラー画像ＣＩにおいて、腫瘤とシストとを判別することができる

ちなみに、上述においては、シストにおいて、Ｆ（ｎ）の演算値が大きくなるようなＰ（ｎ）及びＱ（ｎ）になっているが、シストにおいて、Ｆ（ｎ）の演算値が小さくなるようなＰ（ｎ）及びＱ（ｎ）になっていてもよい。具体的には、Ｐ（ｎ）及びＱ（ｎ）は下記（式２′）及び（式３′）であってもよい。
Ｐ（ｎ）＝ｋ１×｛１／Ａｂｓ（ｓｔｒａｉｎ（ｎ））｝ ・・・（式２′）
Ｑ（ｎ）＝ｋ２×ｘＣｏｒｒ（ｎ） ・・・（式３′）

Ｐ（ｎ）及びＱ（ｎ）が、前記（式２′）及び（式３′）である場合、前記カラー画像ＣＩにおいて、シストは青色で表示される。一方、腫瘤においては、Ｆ（ｎ）の演算値はシストよりも大きくなるので、前記カラー画像ＣＩにおいて、腫瘤はシストとは異なる色で表示される。

次に、第一実施形態の変形例について説明する。先ず、第一変形例について説明する。本例では、前記演算部６は、下記（式１１）を用いて演算を行なう。
Ｆ（ｎ）＝Ｐ（ｎ）×Ｑ（ｎ）×Ｒ（ｎ） ・・・（式１１）

上記（式１１）において、Ｐ（ｎ）及びＱ（ｎ）は、上述の（式２）及び（式３）である。また、Ｒ（ｎ）は、下記（式１２）である。
Ｒ（ｎ）＝ｋ３×（１／Ｉｎｔｅｎｓｉｔｙ（ｎ）） ・・・（式１２）

（式１２）において、ｋ３は任意に設定される重み係数である。このｋ３はデフォルトで設定されていてもよいし、操作者が任意に設定できるようになっていてもよい。また、Ｉｎｔｅｎｓｉｔｙ（ｎ）は、前記音線上の点１〜Ｎにおけるエコー信号の強度である。エコー信号の強度は、Ｂモード画像用の超音波の送受信によって得られたエコー信号の強度であってもよいし、前記物理量データ作成部５が歪みを算出するための超音波の送受信によって得られたエコー信号の強度であってもよい。

ここで、シストにおいては、エコー信号の強度は小さいので、Ｒ（ｎ）は大きい値になる。従って、シストにおいては、Ｆ（ｎ）の演算値は大きくなる。

また、前記（式１１）において、Ｐ（ｎ）及びＱ（ｎ）は、上述の（式２′）及び（式３′）であってもよい。この場合、Ｒ（ｎ）は下記（式１２′）である。
Ｒ（ｎ）＝ｋ３×Ｉｎｔｅｎｓｉｔｙ（ｎ） ・・・（式１２′）
この場合、シストにおいては、Ｆ（ｎ）の演算値は小さくなる。

以上説明した本例においても、シストにおいてはＦ（ｎ）の演算値が大きくなるか小さくなるので、シストが抽出された画像を表示することができる。

次に、第二変形例について説明する。本例では、前記演算部６は、下記（式２１）を用いて演算を行なう。
Ｆ（ｎ）＝Ｐ（ｎ）×Ｒ（ｎ） ・・・（式２１）

上記（式２１）において、Ｐ（ｎ）及びＲ（ｎ）は、上述の（式２）及び（式１２）である。この場合、シストにおいては、Ｆ（ｎ）の演算値は大きくなる。また、上記（式２１）において、Ｐ（ｎ）及びＲ（ｎ）は、上述の（式２′）及び（式１２′）であってもよい。この場合、シストにおいては、Ｆ（ｎ）の演算値は小さくなる。

以上説明した本例においても、シストにおいてはＦ（ｎ）の演算値が大きくなるか小さくなるので、シストが抽出された画像を表示することができる。

次に、第三変形例について説明する。本例では、前記演算部６は、下記（式３１）を用いて演算を行なう。
Ｆ（ｎ）＝Ｑ（ｎ）×Ｒ（ｎ） ・・・（式３１）

上記（式３１）において、Ｑ（ｎ）及びＲ（ｎ）は、上述の（式３）及び（式１２）である。この場合、シストにおいては、Ｆ（ｎ）の演算値は大きくなる。また、上記（式３１）において、Ｑ（ｎ）及びＲ（ｎ）は、上述の（式３′）及び（式１２′）であってもよい。この場合、シストにおいては、Ｆ（ｎ）の演算値は小さくなる。

以上説明した本例においても、シストにおいてはＦ（ｎ）の演算値が大きくなるか小さくなるので、シストが抽出された画像を表示することができる。

（第二実施形態）
次に、第二実施形態について説明する。なお、以下の説明では、第一実施形態と異なる事項について説明する。

図４に示す本例の超音波診断装置２０は、第一実施形態と同様に、超音波プローブ２、送受信ビームフォーマ３、Ｂモードデータ作成部４、物理量データ作成部５、演算部６、表示制御部７、表示部８、操作部９、制御部１０及び記憶部１１を備える。ただし、本例の超音波診断装置２０においては、物理量データ作成部５で作成された物理量データが、前記演算部６のみならず、前記表示制御部７にも入力される。また、前記表示制御部７は、図５に示すように、Ｂモード画像データ作成部７１、カラー画像データ作成部７２、合成画像表示制御部７３のほか、弾性画像データ作成部７４を有している。この弾性画像データ作成部７４は、前記物理量データを、色に対応する情報に変換して、歪みに応じた色に対応する情報を有する弾性画像データを作成する。前記弾性画像データ作成部７４は、前記演算値のデータを階調化し、各階調に割り当てられた色に対応する情報からなる弾性画像データを作成する。階調数は、例えば２５６である。弾性画像データに基づいて表示される弾性画像の色のバリエーションは、前記カラー画像ＣＩと同じでもよく、例えば歪みが小さい場合は青、歪みが大きい場合は赤で表示されるようになっていてもよい。

本例の作用について説明する。本例の超音波診断装置２０では、前記合成画像表示制御部７３が、前記図３に示すように、Ｂモード画像ＢＩとカラー画像ＣＩとが合成された合成超音波画像ＵＩと、図６に示すように、Ｂモード画像ＢＩと弾性画像ＥＩとが合成された合成超音波画像ＵＩ′とを、前記表示部８に切り替えて表示させる。前記合成画像表示制御部７３は、前記操作部９における操作者の入力に基づいて、前記合成超音波画像ＵＩと前記合成超音波画像ＵＩ′とを切り替えて表示させるようになっていてもよい。前記合成画像表示制御部７３は、本発明における画像表示制御部の実施の形態の一例である。

また、前記合成画像表示制御部７３は、図７に示すように、前記合成超音波画像ＵＩと前記合成超音波画像ＵＩ′とを並べて前記表示部８に表示させてもよい。

第二実施形態の超音波診断装置２０によれば、第一実施形態と同様の効果を得ることができるほか、操作者が、前記合成超音波画像ＵＩと前記合成超音波画像ＵＩ′とを切り替えて表示したり並べて表示したりすることにより、シストと腫瘤とをより確実に判別することができる。

以上、本発明を上記各実施形態によって説明したが、本発明はその主旨を変更しない範囲で種々変更実施可能なことはもちろんである。

１，２０ 超音波診断装置
５ 物理量データ作成部
６ 演算部
８ 表示部
７３ 合成画像表示制御部
ＣＩ カラー画像

Claims (8)

1. 生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって前記生体組織における各部の弾性に関する物理量を算出する物理量算出部と、
   該物理量算出部によって算出された物理量の絶対値、前記相関演算における相関係数及び前記生体組織から得られた超音波のエコー信号の強度の三つのパラメータのうち、少なくとも二つのパラメータを用いた演算式であって、前記生体組織におけるシストのパラメータ特性が強調されてシスト以外とは区別できる演算値が得られる演算式の演算値を演算する演算部と、
   該演算部の演算値に応じた表示形態を有する画像が表示される表示部と、
   を備えることを特徴とする超音波診断装置。
2. 前記演算式は、下記（式１）であることを特徴とする請求項１に記載の超音波診断装置。
   （式１）
   Ｆ（ｎ）＝Ｐ（ｎ）×Ｑ（ｎ）
   ただし、Ｐ（ｎ）＝ｋ１×Ａｂｓ（ｓｔｒａｉｎ（ｎ））
   Ｑ（ｎ）＝ｋ２×（１／ｘＣｏｒｒ（ｎ））
   または、Ｐ（ｎ）＝ｋ１×｛１／ａｂｓ（ｓｔｒａｉｎ（ｎ））｝
   Ｑ（ｎ）＝ｋ２×ｘＣｏｒｒ（ｎ）
   ｎ：前記物理量算出部によって歪みが算出される音線上の点を示す１〜Ｎまでの自然数
   ｋ１，ｋ２：重み係数
   ｓｔｒａｉｎ（ｎ）：前記物理量算出部によって算出された音線上の点１〜Ｎにおける歪みの値
   Ａｂｓ（ｓｔｒａｉｎ（ｎ））：前記物理量算出部によって算出された歪みの絶対値を返す関数
   ｘＣｏｒｒ（ｎ）：音線上の点１〜Ｎにおける相関係数の値
3. 前記演算式は、下記（式１１）であることを特徴とする請求項１に記載の超音波診断装置。
   （式１１）
   Ｆ（ｎ）＝Ｐ（ｎ）×Ｑ（ｎ）×Ｒ（ｎ）
   ただし、Ｐ（ｎ）＝ｋ１×Ａｂｓ（ｓｔｒａｉｎ（ｎ））
   Ｑ（ｎ）＝ｋ２×（１／ｘＣｏｒｒ（ｎ））
   Ｒ（ｎ）＝ｋ３×（１／Ｉｎｔｅｎｓｉｔｙ（ｎ））
   または、Ｐ（ｎ）＝ｋ１×｛１／ａｂｓ（ｓｔｒａｉｎ（ｎ））｝
   Ｑ（ｎ）＝ｋ２×ｘＣｏｒｒ（ｎ）
   Ｒ（ｎ）＝ｋ３×Ｉｎｔｅｎｓｉｔｙ（ｎ）
   ｎ：前記物理量算出部によって歪みが算出される音線上の点を示す１〜Ｎまでの自然数
   ｋ１，ｋ２，ｋ３：重み係数
   ｓｔｒａｉｎ（ｎ）：前記物理量算出部によって算出された音線上の点１〜Ｎにおける歪みの値
   Ａｂｓ（ｓｔｒａｉｎ（ｎ））：前記物理量算出部によって算出された歪みの絶対値を返す関数
   ｘＣｏｒｒ（ｎ）：超音波の音線上の点１〜Ｎにおける相関係数の値
   Ｉｎｔｅｎｓｉｔｙ（ｎ）：超音波の音線上の点１〜Ｎにおけるエコー信号の強度
4. 前記演算式は、下記（式２１）であることを特徴とする請求項１に記載の超音波診断装置。
   （式２１）
   Ｆ（ｎ）＝Ｐ（ｎ）×Ｒ（ｎ）
   ただし、Ｐ（ｎ）＝ｋ１×Ａｂｓ（ｓｔｒａｉｎ（ｎ））
   Ｒ（ｎ）＝ｋ３×（１／Ｉｎｔｅｎｓｉｔｙ（ｎ））
   または、Ｐ（ｎ）＝ｋ１×｛１／ａｂｓ（ｓｔｒａｉｎ（ｎ））｝
   Ｒ（ｎ）＝ｋ３×Ｉｎｔｅｎｓｉｔｙ（ｎ）
   ｎ：前記物理量算出部によって歪みが算出される音線上の点を示す１〜Ｎまでの自然数
   ｋ１，ｋ３：重み係数
   ｓｔｒａｉｎ（ｎ）：前記物理量算出部によって算出された音線上の点１〜Ｎにおける歪みの値
   Ａｂｓ（ｓｔｒａｉｎ（ｎ））：前記物理量算出部によって算出された歪みの絶対値を返す関数
   Ｉｎｔｅｎｓｉｔｙ（ｎ）：超音波の音線上の点１〜Ｎにおけるエコー信号の強度
5. 前記演算式は、下記（式３１）であることを特徴とする請求項１に記載の超音波診断装置。
   （式３１）
   Ｆ（ｎ）＝Ｑ（ｎ）×Ｒ（ｎ）
   ただし、Ｑ（ｎ）＝ｋ２×（１／ｘＣｏｒｒ（ｎ））
   Ｒ（ｎ）＝ｋ３×（１／Ｉｎｔｅｎｓｉｔｙ（ｎ））
   または、Ｑ（ｎ）＝ｋ２×ｘＣｏｒｒ（ｎ）
   Ｒ（ｎ）＝ｋ３×Ｉｎｔｅｎｓｉｔｙ（ｎ）
   ｎ：前記物理量算出部によって歪みが算出される音線上の点を示す１〜Ｎまでの自然数
   ｋ２，ｋ３：重み係数
   ｘＣｏｒｒ（ｎ）：超音波の音線上の点１〜Ｎにおける相関係数の値
   Ｉｎｔｅｎｓｉｔｙ（ｎ）：超音波の音線上の点１〜Ｎにおけるエコー信号の強度
6. 前記演算部の演算値に応じた表示形態を有する画像と、前記物理量に応じた表示形態を有する弾性画像とを、前記表示部に切り替えて表示させる画像表示制御部を備えることを特徴とする請求項１〜５のいずれか一項に記載の超音波診断装置。
7. 前記物理量算出部は、前記相関ウィンドウ間で複素相互相関の虚数部を演算して前記生体組織の歪みを算出することを特徴とする請求項１〜６のいずれか一項に記載の超音波診断装置。
8. コンピュータに、
   生体組織に対する超音波の送受信により得られた同一音線上の時間的に異なる二つのエコー信号に相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって前記生体組織における各部の弾性に関する物理量を算出する物理量算出機能と、
   該物理量算出機能によって算出された物理量の絶対値、前記相関演算における相関係数及び前記生体組織から得られた超音波のエコー信号の強度の三つのパラメータのうち、少なくとも二つのパラメータを用いた演算式であって、前記生体組織におけるシストのパラメータ特性が強調されてシスト以外とは区別できる演算値が得られる演算式の演算値を演算する演算機能と、
   該演算機能によって得られた演算値に応じた表示形態を有する画像を表示させる表示画像制御機能と、
   を実行させることを特徴とする超音波診断装置の制御プログラム。

Images