(57)【要約】  弾性ボリュームデータのクオリティを評価する手法を確立して、3次元弾性画像のクオリティを向上させるために、被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを算出する2次元弾性画像構成部116と、弾性フレームデータを複数収集して弾性ボリュームデータを生成する弾性ボリュームデータ生成部117と、弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成する3次元弾性画像構成部118と、弾性フレームデータ算出の基となった一対の断層フレームデータ間の自己相関値と、弾性フレームデータのクオリティを示すフレーム評価値に基づいて、弾性ボリュームデータのクオリティを示すボリューム評価値を算出するクオリティ算出部121を備えたことを特徴とする超音波診断装置。
2012029417.tif 000002
【特許請求の範囲】
【請求項１】
被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを複数収集して弾性ボリュームデータを生成する弾性ボリューム生成部と、前記弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成する3次元弾性画像構成部と、前記3次元弾性画像を表示する表示部とを備えた超音波診断装置において、
前記弾性フレームデータのクオリティを示すフレーム評価値に基づいて、前記弾性ボリュームデータのクオリティを示すボリューム評価値を算出するクオリティ算出部を備えたことを特徴とする超音波診断装置。
【請求項２】
請求項1に記載の超音波診断装置において、
前記クオリティ算出部は、前記弾性フレームデータ算出の基となった一対の断層フレームデータ間の自己相関値を前記弾性フレームデータのクオリティを示すフレーム評価値とし、前記弾性ボリュームデータを構成する全ての前記弾性フレームデータの前記フレーム評価値の加算値又は加算平均値に基づいて、前記弾性ボリュームデータのクオリティを示すボリューム評価値を求めることを特徴とする超音波診断装置。
【請求項３】
請求項2に記載の超音波診断装置において、
前記クオリティ算出部は、前記一対の断層フレームデータ間の自己相関値として、前記3次元スキャンにおける時間的に隣り合う一対の断層フレームデータの自己相関値を用いることを特徴とする超音波診断装置。
【請求項４】
請求項2に記載の超音波診断装置において、
前記クオリティ算出部は、前記一対の断層フレームデータ間の自己相関値として、前記弾性ボリューム生成部で繰り返し生成される複数の弾性ボリュームデータをそれぞれ構成する複数の前記弾性フレームデータのうち、スキャン面位置が同一の前記弾性フレームデータ算出の基となった前記一対の断層フレームデータの自己相関値を用いることを特徴とする超音波診断装置。
【請求項５】
請求項2に記載の超音波診断装置において、
前記クオリティ算出部は、前記一対の断層フレームデータ間の自己相関値として、現在の断層フレームデータと過去の複数の断層フレームデータと間で求めた前記自己相関値が最大の一対の断層フレームデータの自己相関値を用いることを特徴とする超音波診断装置。
【請求項６】
請求項5に記載の超音波診断装置において、
前記自己相関値が最大の一対の断層フレームデータに基づいて算出された前記弾性フレームデータが収集され前記弾性ボリュームデータが生成されることを特徴とする超音波診断装置。
【請求項７】
請求項1に記載の超音波診断装置において、
前記クオリティ算出部は、前記弾性フレームデータの弾性値の分布の平均、偏差、又はS/N比に基づいて、各弾性フレームデータのクオリティを示すフレーム評価値を算出し、該算出したフレーム評価値を前記弾性ボリュームデータ単位で加算平均することにより前記弾性ボリュームデータのクオリティを示すボリューム評価値を算出することを特徴とする超音波診断装置。
【請求項８】
請求項7に記載の超音波診断装置において、
前記弾性値は、変位、歪み、歪み比、粘性、弾性率のいずれか1つであることを特徴とする超音波診断装置。
【請求項９】
請求項1に記載の超音波診断装置において、
前記クオリティを示すボリューム評価値を前記3次元弾性画像に対応付けて前記表示部に表示することを特徴とする超音波診断装置。
【請求項１０】
請求項9に記載の超音波診断装置において、
前記クオリティを示すボリューム評価値を、該ボリューム評価値に対応して表示態様が異なるマーク、バーチャート又は円グラフで前記表示部に表示することを特徴とする超音波診断装置。
【請求項１１】
請求項1に記載の超音波診断装置において、
前記弾性ボリュームデータ生成部は、前記弾性ボリュームデータに基づいて構成した直交3断面弾性画像又は平行な複数の断面でスライスしたマルチ断面弾性画像を前記表示部に表示する断面弾性画像構成部を備え、
前記クオリティ算出部は、前記各断面弾性画像のクオリティを示すボリューム評価値を算出して前記表示部に表示する表示面クオリティ算出部を備えてなることを特徴とする超音波診断装置。
【請求項１２】
請求項1に記載の超音波診断装置において、
前記弾性ボリュームデータのクオリティを示す前記ボリューム評価値が表示基準値よりも高い場合、前記表示部は前記3次元弾性画像を表示することを特徴とする超音波診断装置。
【請求項１３】
請求項1に記載の超音波診断装置において、
前記弾性ボリュームデータのクオリティを示すボリューム評価値が記憶基準値よりも高い場合、前記3次元弾性画像を記憶する記憶部を備えることを特徴とする超音波診断装置。
【請求項１４】
被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを複数収集して弾性ボリュームデータを生成するステップと、弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成するステップと、前記3次元弾性画像を表示するステップと、前記弾性フレームデータのクオリティを示すフレーム評価値に基づいて、前記弾性ボリュームデータのクオリティを示すボリューム評価値を算出するステップを含む評価算出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波を利用して被検体の生体組織の硬さ又は軟らかさを示す3次元弾性画像を表示する超音波診断装置と評価算出方法に関する。
【背景技術】
【０００２】
超音波診断装置は、被検体に超音波を送信し、被検体内部の生体組織からの反射エコー信号を受信し、例えば、生体組織の硬さ又は軟らかさを示す3次元弾性画像などの超音波画像を生成してモニタに表示して診断に供するものである。また、被検体内部の形態を分かり易く表示して診断に寄与するために、3次元断層画像や3次元弾性画像の3次元超音波画像を生成して表示することが行われている。一方で、診断精度の向上のため、画像の分解能の向上、ノイズの低減、等々、超音波画像の画質を改善することが要請されている。【０００３】
例えば、特許文献1には、弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構築する際に、同じ変位もしくは圧力の弾性画像区間を合成して3次元弾性画像を構築することが提案されている。さらに、同文献には、相関係数が高い弾性フレームで弾性画像を合成して、3次元弾画像を構築することが提案されている。つまり、3次元弾性画像を構成する際に、連続して取得した複数の2次元弾性フレームデータから、圧迫量(変位)が同じ2次元弾性フレームデータにより弾性ボリュームデータを構成し、その弾性ボリュームデータをボリュームレンダリングすることにより、3次元弾性画像を構築するようにしている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開2008-259555号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、特許文献では、弾性ボリュームデータの中から所望の弾性フレームデータを選択して弾性ボリュームデータを合成しているが、取得した弾性ボリュームデータ全体についてのクオリティを評価することは考慮されていない。その結果、ボリュームレンダリングにより構築した3次元弾性画像にノイズ(例えば、スジ状のノイズ画像)が表れる可能性がある。【０００６】
本発明が解決しようとする課題は、弾性ボリュームデータのクオリティを評価する手法を確立して、3次元弾性画像のクオリティを向上させることにある。
【課題を解決するための手段】
【０００７】
上記の課題を解決するため、本発明は、被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを算出する弾性演算部と、前記弾性フレームデータを複数収集して弾性ボリュームデータを生成する弾性ボリューム生成部と、前記弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成する3次元弾性画像構成部と、前記3次元弾性画像を表示する表示部とを備えた超音波診断装置において、前記弾性フレームデータのクオリティを示すフレーム評価値基づいて、前記弾性ボリュームデータのクオリティを示すボリューム評価値を算出するクオリティ算出部を備えたことを特徴とする。【０００８】
以上のようにして求めたクオリティを示すボリューム評価値を、3次元弾性画像に対応付けて表示部に表示することが望ましい。これにより、表示された3次元弾性画像のクオリティがどの程度高いか否かを判断することができる。また、クオリティを示すボリューム評価値に対応して表示態様が異なるマーク、バーチャート又は円グラフを表示部に表示することにより、一目でクオリティを示すボリューム評価値を識別できるから、診断の精度向上に寄与することができる。【０００９】
さらに、本発明は、3次元弾性画像に限らず、弾性ボリュームデータに基づいて構成した直交3断面弾性画像又は平行な複数の断面でスライスしたマルチ断面弾性画像を前記表示部に表示する断面弾性画像構成部を備えて構成し、クオリティ算出部により各断面弾性画像のクオリティを示すボリューム評価値を算出して表示部に表示する表示面クオリティ算出部を備えて構成することができる。【００１０】
以上説明したように、本発明によれば、3次元弾性画像のクオリティを表示することができるから、検査者が3次元弾性画像の選択を簡便に行うことができる。また、クオリティの高い3次元弾性像を再構成して表示することもでき、高レベルな診断サポートが可能となる。
【発明の効果】
【００１１】
本発明によれば、弾性ボリュームデータのクオリティを評価する手法を確立できるから、3次元弾性画像のクオリティを向上させることができる。
【図面の簡単な説明】
【００１２】
【図１】本発明の実施例1の超音波診断装置の全体を示すブロック構成図
【図２】実施例1のクオリティ算出部のブロック構成図
【図３】実施例1の特徴部の処理手順を説明する図
【図４】本発明の実施例2の特徴部の処理手順を説明する図
【図５】本発明の実施例3の特徴部の処理手順を説明する図
【図６】本発明の実施例3の変形例の処理手順を説明する図
【図７】本発明の実施例4の断面弾性画像の表示例を説明する図
【図８】本発明の実施例4の断面弾性画像の表示を実現する主要部のブロック構成図
【図９】弾性ボリュームデータのクオリティを示すボリューム評価値の表示方法の一例を示す図
【図１０】弾性ボリュームデータのクオリティを示すボリューム評価値の表示方法の他の例を示す図
【発明を実施するための形態】
【００１３】
各図を用いて実施例について説明する。
【実施例１】
【００１４】
本発明を適用した一実施例の超音波診断装置100について、図1を用いて説明する。
図示のように、超音波診断装置100には、被検体101に当接させて用いる超音波探触子102と、超音波探触子102を介して被検体101に一定の時間間隔をおいて超音波を繰り返し送信させる送信部105と、被検体101の内部からの反射エコー信号を受信する受信部106と、送信部105と受信部106を制御する送受信制御部107と、受信部106で受信された反射エコーを整相加算する整相加算部108とが備えられている。【００１５】
超音波探触子102は、被検体101に振動子を介して超音波を送受信する機能を有している。超音波探触子102は、矩形又は扇形をなす複数の振動子を配列してなり、複数の振動子の配列方向(長軸方向)と直交する方向(短軸方向)に超音波探触子102を機械的に振り、超音波を3次元スキャンに可能に構成されている。なお、超音波の3次元スキャンは、超音波探触子102を機械的に短軸方向に振るものに限らず、例えば、複数の振動子を2次元に配列した超音波探触子102を用い、超音波を短軸方向に電子的にスキャンするものでもよい。【００１６】
送信部105は、超音波探触子102の振動子を駆動して超音波を発生させるための送波パルスを生成する。送信部105は、送信される超音波の収束点を任意の深さに設定する機能を有している。また、受信部106は、超音波探触子102で受信した反射エコー信号を所定のゲインで増幅してRF信号、すなわち受信信号を生成する。超音波送受信制御部107は、送信部105や受信部106を制御するためのものである。整相加算部108は、受信部106で増幅されたRF信号の位相を制御して、1点又は複数の収束点に対応する超音波ビームを形成することにより、RF信号フレームデータを生成する。【００１７】
整相加算部108で生成されたRF信号フレームデータはデータ記憶部109に記憶される。2次元断層画像構成部113は、データ記憶部109に記憶されたRF信号フレームデータに基づいて2次元断層画像を構成する。断層ボリュームデータ生成部114は、2次元断層画像構成部113で構成された2次元断層画像を取得位置に基づいて3次元座標変換を行なって断層ボリュームデータを生成する。3次元断層画像構成部115は、断層ボリュームデータの輝度と不透明度に基づいてボリュームレンダリングを行なって3次元断層画像を構成する。【００１８】
データ記憶部109に記憶された複数のRF信号フレームデータは、2次元弾性画像構成部116に適宜取り込まれ、2次元弾性画像が構成される。2次元弾性画像構成部116で構成された2次元弾性画像は、弾性ボリュームデータ生成部117に取り込まれ、2次元弾性画像の取得位置に基づいて3次元座標変換を行なって弾性ボリュームデータが生成される。3次元弾性画像構成部118は、弾性ボリュームデータの弾性値と不透明度に基づいてボリュームレンダリングを行って3次元弾性画像を構成する。合成処理部119は、2次元断層画像と2次元弾性画像を合成したり、3次元断層画像と3次元弾性画像を合成するようになっている。表示部120は、合成処理部119で合成された合成画像や2次元断層画像の超音波画像を表示するようになっている。また、超音波診断装置100には、上述した各構成要素を制御する制御部103と、制御部103に各種入力を行なう操作部104を備え、操作部104は、キーボードやトラックボール等を備えている。【００１９】
以下、図1の主要部の詳細構成について説明する。2次元断層画像構成部113は、制御部103における設定条件に基づいて、データ記憶部109から出力されるRF信号フレームデータを取り込んで、ゲイン補正、ログ圧縮、検波、輪郭強調、フィルタ処理等の信号処理を行なって2次元断層画像を構成する。このとき、超音波探触子102は、超音波の送受信と同時に送受信方向(θ、φ)を計測する。ここで、θは長軸方向のスキャン角度、φは短軸方向のスキャン(スイング)角度である。断層ボリュームデータ生成部114は、2次元断層画像の取得位置に相当する送受信方向(θ、φ)に基づいて、複数の2次元断層画像について3次元変換を行なって断層ボリュームデータを生成する。【００２０】
3次元断層画像構成部115は、断層ボリュームデータから3次元断層画像を構成する次式(1)〜(3)を用いて、ボリュームレンダリングを行なう。【００２１】
Cout(i)  =  Cout(i-1)＋(1−Aout(i-1))・A(i)・C(i)・S(i)   −(1)
Aout(i)  =  Aout(i-1)＋(1−Aout(i-1))・A(i)               −(2)
A(i)  =  Opacity[C(i)]                                    −(3)
ここで、C(i)は作成される2次元投影面上のある点から3次元断層画像を見た場合、視線上i番目に存在するボクセルの輝度値である。Cout(i)は、出力されるピクセル値である。例えば、視線上にＮボクセルの輝度値が並んだとき、i =0〜Ｎ-1までを積算した輝度値Cout(N-1)が最終的に出力されるピクセル値となる。Cout(i-1)はi-1番目までの積算値を示す。【００２２】
また、A(i)は視線上i番目に存在する輝度値の不透明度であり、式(3)に示すとおり、0〜1．0の値をとる断層不透明度テーブル(断層オパシティテーブル)である。断層不透明度テーブルは、輝度値から不透明度を参照することによって、出力する2次元投影面(3次元断層画像)上への寄与率を決定する。【００２３】
S(i)は、輝度C(i)とその周辺の画素値より求めた勾配より算定される陰影付けのための重み成分で、たとえば、光源とボクセルiを中心とした面の法線が一致する場合、もっとも強く反射するため1．0が与えられ、光源と法線が直交する場合には0．0が与えられるなどの強調効果を指し示す。【００２４】
Cout(i)、Aout(i)は、共に0を初期値としている。式(2)に示されるように、Aout(i)はボクセルを通過するたびに積算され1．0に収束される。よって、式(1)に示されるようにi-1番目までの不透明度の積算値Aout(i-1)およそ1．0となった場合、ｉ番目以降の輝度値C(i)は出力画像に反映されない。【００２５】
2次元弾性画像構成部116は、データ記憶部109に記憶された複数のRF信号フレームデータから変位を計測する。2次元弾性画像構成部116は被検体101に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを算出する弾性演算部を有している。そして、2次元弾性画像構成部116は、計測した変位に基づいて弾性値を演算し、2次元弾性画像を構成する。弾性値とは、歪み、弾性率、変位、粘性、歪み比等の弾性情報のいずれかである。弾性ボリュームデータ生成部117は、2次元弾性画像の取得位置に相当する送受信方向(θ、φ)に基づいて、複数の2次元弾性画像について3次元変換を行ない、弾性ボリュームデータを生成する。3次元弾性画像構成部118は、弾性値に基づいて弾性ボリュームデータを複数に分離し、分離された弾性ボリュームデータについてボリュームレンダリングを行ない、3次元弾性画像を構成する。【００２６】
以下、本発明の特徴部であるクオリティ算出部121について説明する。図2に、クオリティ算出部121のブロック構成図を示す。クオリティ算出部121は、基本的に、圧迫の情報もしくは3次元弾性画像の質であるクオリティ(Quality)を検出し、弾性ボリュームデータの質を演算して良し悪しを評価するものである。ここで、検出するクオリティとしては、2次元弾性画像を算出した際の一対の断層フレームデータの自己相関、変位、歪み、歪みのS/N比、等々を用いることができる。また、取得時間が異なる複数の弾性ボリュームデータ間の同一位置の弾性フレームデータの相関を用いることもできる。【００２７】
クオリティ算出部121は、フレーム相関処理部201、フレーム変位・歪み処理部203、ボリューム処理部205、クオリティ算出部207により構成される。フレーム相関処理部201は、2次元弾性画像構成部116で弾性フレームデータを算出した際に用いた一対の断層フレームデータ間の自己相関を保存しておき、フレーム単位の弾性フレームデータのクオリティを示すフレーム評価値を算出するようになっている。【００２８】
また、フレーム変位・歪み処理部203は、2次元弾性画像構成部116で算出された弾性値(歪み、弾性率、変位、粘性、歪み比)と圧力を弾性フレームデータに対応付けて保存し、フレーム単位の弾性フレームデータのクオリティを示すフレーム評価値を算出するようになっている。【００２９】
ボリューム処理部205は、フレーム相関処理部201、フレーム変位・歪み処理部203において、フレーム単位で算出されたフレーム評価値を加算平均等により、弾性ボリュームデータのクオリティを示すボリューム評価値を求める。ここで、ボリューム評価値は1つあればよいが、複数のボリューム評価値を求めることもできる。【００３０】
図3を参照して、本実施例1の詳細な処理手順について説明する。図3(a)に示すように、3次元スキャンは、長軸方向(θ)に電子スキャンを行って断層フレームデータを計測する操作を、短軸方向(φ)にスキャン断面をスライドさせながら複数(n)フレームの断層フレームデータFr.0―Fr.nを取得する。短軸方向(φ)のスキャンは、機械的なスキャン又は電子的なスキャンのいずれを適用してもよい。例えば、モータにより超音波探触子102を遥動させてスキャンする場合は、モータを低速又は高速に速度制御可能にする。【００３１】
3次元スキャンにおいて2次元弾性画像構成部116は、例えば、隣り合う一対の断層フレームデータFr.0、Fr.1間の自己相関演算を行って生体組織の変位分布を推定する。なお、図示していないが、図1の断層ボリュームデータ生成部114に記憶されている異なる時間に取得された同一部位の複数の断層ボリュームデータを参照して、現在(リアルタイム)取得された例えば断層フレームデータFr.0と、同一位置において過去に取得された断層フレームデータFr.0との間で自己相関演算を行って生体組織の変位分布を推定することもできる。このように、一対の断層フレームデータ間で求められた自己相関値が、その一対の断層フレームデータにより算出された弾性フレームデータのクオリティを示すフレーム評価値となる。なお、一対の断層フレーム間の自己相関値が高いということは、断層フレームデータを構成する複数本のRF信号の波形の一致度が高いことになり、これにより得られる弾性フレームデータのクオリティは高い傾向にある。【００３２】
そこで、フレーム相関処理部201は、2次元弾性画像構成部116から各対の断層フレーム間の自己相関値を取り込み、図3(a)、(b)に示すように、各弾性フレームデータについて画素単位で自己相関値のフレーム平均値Caveを求める。次いで、弾性ボリュームデータを構成する全部の弾性フレームデータに対応する自己相関値をフレーム評価値とし、その加算平均又は加算値を求めて、弾性ボリュームデータのクオリティを示すボリューム評価値を求める。例えば、自己相関値の弾性ボリュームの加算平均Cvol=0．95であればクオリティを示すボリューム評価値が高いとして、図3(c)のように「Quality High」が表示部120に表示される。また、その弾性ボリュームデータをボリュームレンダリングして得られる3次元弾性像画像は、崩れが少ないものとなる。一方、図3(d)のように自己相関値の加算平均Cvol=0．74であれば「Quality Low」が表示部120に表示される。この場合の3次元弾性像画像は、ノイズが多く、例えば本来なら円形の画像が変形して表示されるなどの先鋭さに欠ける画像となる。このようにして、弾性ボリュームデータごとにクオリティを示すボリューム評価値を求めることにより、弾性ボリュームデータのクオリティを評価することが可能である。なお、フレーム変位・歪み算出部203の詳細な実施例は、後述する。
【実施例２】
【００３３】
図4を参照して、本発明の実施例2の2次元弾性画像構成部116とクオリティ算出部121の詳細な処理手順について説明する。本実施例2が実施例1と異なる点は、2次元弾性画像構成部116とフレーム相関処理部201の処理手順にあり、その他の点は実施例1と同様であることから説明を省略する。実施例1においては、2次元弾性画像構成部116は隣り合う一対の断層フレームデータ間の自己相関演算を行って生体組織の変位分布を推定していた。これに対し、本実施例2では、2次元弾性画像構成部116は現在(リアルタイム)の断層フレームデータに対する過去数フレーム分の断層フレームデータとの間で自己相関演算を行って生体組織の変位分布を推定する。特に、過去数フレーム分との間の自己相関値が最も高い断層フレームデータを選択して、生体組織の変位分布を推定して弾性フレームデータを作成する。つまり、複数回の3次元スキャンで1フレーム分の自己相関値と変位分布データが算出される。その後の処理は、実施例1と同様であり、図4(c)、(d)に示す評価結果が表示される。【００３４】
本実施例2によれば、複数回の3次元スキャンで自己相関値が高い弾性フレームデータを選択できることから、最終的に質の高い弾性ボリュームデータを生成することができ、3次元弾性画像のクオリティを高めることができる。
【実施例３】
【００３５】
図5を参照して、本発明の実施例3のクオリティ算出部121の詳細な処理手順について説明する。本実施例3が実施例1，2と異なる点は、弾性フレームデータを算出した際の一対の断層フレームデータの自己相関値に代えて、弾性フレームデータの弾性値(変位、歪み、歪み比、粘性、弾性率)の平均、S/N比、偏差情報の少なくとも1つに基づいて、弾性ボリュームデータを構成する弾性フレームデータのクオリティの評価値を算出するようにしたことにある。【００３６】
図5(a)に示すように、短軸方向に2次元スキャンを繰り返して、フレーム間自己相関から検出された弾性(変位)フレームデータを複数(n)フレーム取得する。これにより、同図(b)に示すような2次元の弾性フレームデータの集合体として、3次元の弾性ボリュームデータが作成される。各弾性フレームデータの変位の例えば加算平均を求めて、弾性ボリュームデータのクオリティを示すボリューム評価値とする。また、各フレームの変位の加算平均の分散を用いることにより、画像のバラつきを検出できる。各フレームの加算平均の分散Saveは、式(4)により算出する。さらに、式(5)により、弾性ボリュームデータを構成する全ての弾性フレームデータの分散Saveの弾性ボリュームデータ平均Svolを求めて、弾性ボリュームデータのクオリティを示すボリューム評価値とする。
2012029417.tif 000003
ここで、Ｗは画像の横幅、Ｈは画像の高さ、d(i,j)は変位、d*は変位の平均である。
2012029417.tif 000004
本実施例によれば、バラつきが少ない方がノイズは少ないから、図5(c)に示すように、Svol＝0.001のように小さな値の方が評価が高いボリューム評価値「Quality High」と表示され、3次元弾性画像としても引き締まった画像となる。一方、ノイズの多い画像においては、バラつきが大きくなり、図5(d)に示すように、Svol＝0.02のように大きな値の方が評価が低いボリューム評価値「Quality Low」と表示され、画像も先鋭さを欠ける画像となる。【００３７】
図6は、本実施例の変形例であり、変位フレームデータに代えて、歪みフレームデータに基づいて弾性ボリュームデータのクオリティを評価するようにした例である。図6(a)に示すように、複数(n)フレームの歪みフレームデータにより弾性ボリュームデータが構成される。一般に、生体組織に圧迫が加えられていれば、負の歪みは存在しえない。したがって、図6(b)に示すように、一の歪みフレームデータ中に負の歪みが混在していると、歪みの平均が小さくなる。そこで、各フレームの歪みの平均をφ方向の各歪みフレームデータについて加算平均を求めて、弾性ボリュームデータのクオリティを示すボリューム評価値とする。これによれば、図6(c)に示すように評価が高い方が「Quality High」の表示になり、評価が低い場合は図6(d)に示すように「Quality Low」の表示になる。
【実施例４】
【００３８】
図7に、本発明の各実施例により求めた弾性ボリュームデータを用いて構成される各種の表示画像の例を示す。図7の例では、弾性ボリュームデータを用いて任意の断面における断面画像を表示する場合において、表示された各断面画像にそれぞれクオリティを示すボリューム評価値を表示するようにしている。具体的には、図7(a)に示すように、任意の直交3断面のアキシャル面A、シータ面T，コロナル面C、3次元弾性画像Renderのそれぞれに、クオリティを示すボリューム評価値を表示する。また、図7(b)に示すように、いわゆるマルチスライス断面画像の場合も、複数断面画像の表示にそれぞれに、クオリティを示すボリューム評価値を表示することができる。【００３９】
図7の画像表示例を実現するには、弾性ボリュームデータ生成部117とクオリティ算出部121を図8のように構成する。つまり、弾性ボリューム生成部117にて、2次元弾性像構成部116で算出された自己相関値は、弾性ボリューム生成部117において、自己相関値ボリューム保持部802又は変位・歪みボリューム保持部804に格納される。このデータは、超音波のRF配列データ曲座標(RTΦ)次元のデータであり、直交座標変換部806にて直交座標のデータに変換される。制御部103において、表示用に参照される座標面に対して表示面クオリティ算出部808で各断面のクオリティを示すボリューム評価値を算出する。
【実施例５】
【００４０】
図9及び図10を参照して、各実施例で評価した弾性ボリュームデータを示すボリューム評価値の表示方法の実施例について説明する。図9は、横軸が時間、縦軸がクオリティを示すボリューム評価値を示しており、連続的に取得した複数の弾性ボリュームデータV0−Vnのボリューム評価値が一目でわかる。また、最もボリューム評価値が高い弾性ボリュームデータと、最もボリューム評価値が低い弾性ボリュームデータに対応付けて、図3と同様のマークを表示する。【００４１】
図10は、弾性ボリュームデータ又は3次元弾性画像のクオリティを示すボリューム評価値を表すマーク、バーチャート、円グラフ、その他の変形例を示している。図10(a)は、円形のマークの表示色をボリューム評価値に対応付けて変えることにより、一目でボリューム評価値の高低を識別可能にした例である。同図(b)は、実施例1等で示したマークと同じであり、円形のマークの真円度合いに応じてボリューム評価値の高低を表した例である。同図(c)は、バーチャートのゲージ内の比率に応じてボリューム評価値の高低を表した例である。同図(d)は、円グラフの比率に応じてボリューム評価値の高低を表した例である。同図(e)は、複数の小円を並べてボリューム評価値に応じて、小円の表示態様の比率を異ならせた例である。
【実施例６】
【００４２】
本発明の各実施例により求めた弾性ボリュームデータを用いて構成される3次元弾性画像の表示形態と記憶形態を説明する。【００４３】
クオリティ算出部121によって算出された弾性ボリュームデータのクオリティを示すボリューム評価値に基づいて、表示部120は3次元弾性画像を表示する。ボリューム評価値が表示基準値よりも高い場合、制御部103は3次元弾性画像を表示するように指示を行い、表示部120はボリューム評価値が表示基準値よりも高い3次元弾性画像を表示する。ボリューム評価値が表示基準値よりも低い場合、制御部103は3次元弾性画像を表示しないように指示を行い、表示部120はボリューム評価値が表示基準値よりも低い3次元弾性画像を表示しない。【００４４】
表示基準値は、例えば、0.95とする。また、操作者は操作部104で表示基準値を設定することができる。【００４５】
本実施例によれば、弾性ボリュームデータのクオリティを示すボリューム評価値が表示基準値よりも高い場合、表示部120は3次元弾性画像を表示する。つまり、表示部120は弾性ボリュームデータのクオリティを示すボリューム評価値が表示基準値よりも高い3次元弾性画像のみを表示することができる。【００４６】
また、クオリティ算出部121によって算出された弾性ボリュームデータのクオリティを示すボリューム評価値に基づいて、記憶部(図示しない。)は3次元弾性画像を記憶する。ボリューム評価値が記憶基準値よりも高い場合、制御部103は3次元弾性画像を記憶部に記憶するように指示を行い、記憶部はボリューム評価値が記憶基準値よりも高い3次元弾性画像を記憶する。ボリューム評価値が記憶基準値よりも低い場合、制御部103は3次元弾性画像を記憶部に記憶しないように指示を行い、記憶部はボリューム評価値が記憶基準値よりも低い3次元弾性画像を記憶しない。【００４７】
本実施例によれば、弾性ボリュームデータのクオリティを示すボリューム評価値が記憶基準値よりも高い場合、3次元弾性画像を記憶する記憶部を備える。つまり、記憶部は弾性ボリュームデータのクオリティを示すボリューム評価値が記憶基準値よりも高い3次元弾性画像のみを記憶することができる。【００４８】
記憶基準値は、表示基準値と同様に設定することができ、操作者は操作部104で記憶基準値を設定することができる。表示基準値と記憶基準値を同一にすることもできる。【００４９】
以上、本発明によれば、被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを複数収集して弾性ボリュームデータを生成する弾性ボリューム生成部117と、弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成する3次元弾性画像構成部118と、3次元弾性画像を表示する表示部120とを備えた超音波診断装置において、弾性フレームデータのクオリティを示すフレーム評価値基づいて、弾性ボリュームデータのクオリティを示すボリューム評価値を算出するクオリティ算出部121を備えたことを特徴とする。また、被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを複数収集して弾性ボリュームデータを生成するステップと、弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成するステップと、前記3次元弾性画像を表示するステップと、前記弾性フレームデータのクオリティを示すフレーム評価値に基づいて、前記弾性ボリュームデータのクオリティを示すボリューム評価値を算出するステップを含む評価算出方法である。【００５０】
すなわち、弾性フレームデータのクオリティを示すフレーム評価値が高いことは、安定した計測状態で弾性フレームデータが算出されたことになる。これに鑑み、本発明は、弾性フレームデータのクオリティを評価するフレーム評価値を求め、さらに弾性ボリュームデータを構成する複数の弾性フレームデータを示すフレーム評価値に基づいて、弾性ボリュームデータのクオリティを示すボリューム評価値を求めて評価するようにしたのである。このように、弾性ボリュームデータのクオリティを評価する手法が確立されたことから、クオリティを示すボリューム評価値が高い弾性ボリュームデータを取得又は選択することにより、容易に3次元弾性画像のクオリティを向上させることができる。【００５１】
ここで、弾性フレームデータのクオリティあるいはクオリティ評価値とは、安定した適切な圧迫状態において弾性フレームデータが計測されたことを意味する。同様に、弾性ボリュームデータのクオリティあるいはクオリティ評価値とは、クオリティあるいはクオリティ評価値が高い弾性フレームデータの集合であることを意味する。したがって、クオリティあるいはクオリティ評価値が高いことは、結果としてノイズの少ない3次元弾性画像を生成することができる。【００５２】
ところで、弾性フレームデータのクオリティを示すフレーム評価値は、クオリティ評価方法により求めることができる。例えば、弾性フレームデータ算出の基となった一対の断層フレームデータ間の自己相関値が高いと、一対の断層フレームデータの一致度が高いことを意味し、その一対の断層フレームデータが安定した計測状態で計測されたことになる。そこで、弾性フレームデータ算出の基となった一対の断層フレームデータの自己相関値をフレーム評価値とし、これに基づいて弾性ボリュームデータのクオリティを示すボリューム評価値を求めて評価することができる。しかし、本発明に係るフレーム評価値は、後述するように、一対の断層フレームデータの自己相関値に限られるものではない。【００５３】
本発明において、クオリティ算出部121は、一対の断層フレームデータ間の自己相関値を弾性フレームデータのクオリティを示すフレーム評価値とし、弾性ボリュームデータを構成する全ての弾性フレームデータのフレーム評価値の加算値又は加算平均値に基づいて、弾性ボリュームデータのクオリティを示すボリューム評価値を求めるようにすることができる。この場合、一対の断層フレームデータ間の自己相関値として、3次元スキャンにおける時間的に隣り合う一対の断層フレームデータの自己相関値を用いることができる。また、これに代えて、一対の断層フレームデータ間の自己相関値として、弾性ボリューム生成部で繰り返し生成される複数の弾性ボリュームデータをそれぞれ構成する複数の弾性フレームデータのうち、スキャン面位置が同一の前記弾性フレームデータ算出の基となった前記一対の断層フレームデータの自己相関値を用いることができる。【００５４】
さらに、一対の断層フレームデータ間の自己相関値として、現在の断層フレームデータと過去の複数の断層フレームデータと間で求めた自己相関値が最大の一対の断層フレームデータの自己相関値を用いることができる。これによれば、クオリティの高い弾性フレームデータを選択して弾性ボリュームデータを構成できることから、一層、ボリューム評価値を高めることができる。つまり、自己相関値が最大の一対の断層フレームデータに基づいて算出された弾性フレームデータを収集して弾性ボリュームデータが生成される。なお、前述したように、本発明は、弾性フレームデータ算出の基となった一対の断層フレームデータの自己相関値に基づいてフレーム評価値を求めるものに限られない。これに代えて、弾性フレームデータの弾性値の分布の平均、偏差、又はS/N比に基づいて、各弾性フレームデータのクオリティを示すフレーム評価値を算出することができる。例えば、弾性値のS/N比が大きければ弾性フレームデータのクオリティが高いと評価することができる。また、弾性値としては、変位、歪み、歪み比、粘性、弾性率のいずれか1つを用いることができる。【００５５】
以上のようにして求めたクオリティを示すボリューム評価値を、3次元弾性画像に対応付けて表示部120に表示することが望ましい。これにより、表示された3次元弾性画像のクオリティがどの程度高いか否かを判断することができる。また、クオリティを示すボリューム評価値に対応して表示態様が異なるマーク、バーチャート又は円グラフを表示部に表示することにより、一目でクオリティを示すボリューム評価値を識別できるから、診断の精度向上に寄与することができる。【００５６】
さらに、本発明は、3次元弾性画像に限らず、弾性ボリュームデータに基づいて構成した直交3断面弾性画像又は平行な複数の断面でスライスしたマルチ断面弾性画像を前記表示部に表示する断面弾性画像構成部を備えて構成し、クオリティ算出部により各断面弾性画像のクオリティを示すボリューム評価値を算出して表示部に表示する表示面クオリティ算出部を備えて構成することができる。【００５７】
以上説明したように、本発明によれば、3次元弾性画像のクオリティを表示することができるから、検査者が3次元弾性画像の選択を簡便に行うことができる。また、クオリティの高い3次元弾性像を再構成して表示することもでき、高レベルな診断サポートが可能となる。
【符号の説明】
【００５８】
100  超音波診断装置、102  超音波探触子、103  制御部、104  操作部  、105  送信部、106  受信部、107  送受信制御部、108  整相加算部、109  データ記憶部、113  2次元断層画像構成部、114  断層ボリュームデータ生成部、115  3次元断層画像構成部、116  2次元弾性画像構成部、117  弾性ボリュームデータ生成部、118  3次元弾性画像構成部、119  合成処理部、120  表示部、121  クオリティ算出部
【図１】
2012029417.tif 000005
【図２】
2012029417.tif 000006
【図３】
2012029417.tif 000007
【図４】
2012029417.tif 000008
【図５】
2012029417.tif 000009
【図６】
2012029417.tif 000010
【図７】
2012029417.tif 000011
【図８】
2012029417.tif 000012
【図９】
2012029417.tif 000013
【図１０】
2012029417.tif 000014
【手続補正書】
【提出日】平成25年6月5日(2013.6.5)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを複数収集して弾性ボリュームデータを生成する弾性ボリューム生成部と、前記弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成する3次元弾性画像構成部と、前記3次元弾性画像を表示する表示部とを備えた超音波診断装置において、
前記弾性フレームデータのクオリティを示すフレーム評価値に基づいて、前記弾性ボリュームデータのクオリティを示すボリューム評価値を算出するクオリティ算出部を備えたことを特徴とする超音波診断装置。
【請求項２】
請求項1に記載の超音波診断装置において、
前記クオリティ算出部は、前記弾性フレームデータを算出する基となった一対の断層フレームデータ間の自己相関値を前記弾性フレームデータのクオリティを示すフレーム評価値とし、前記弾性ボリュームデータを構成する全ての前記弾性フレームデータの前記フレーム評価値の加算値又は加算平均値に基づいて、前記弾性ボリュームデータのクオリティを示すボリューム評価値を求めることを特徴とする超音波診断装置。
【請求項３】
請求項2に記載の超音波診断装置において、
前記クオリティ算出部は、前記一対の断層フレームデータ間の自己相関値として、前記3次元スキャンにおける時間的に隣り合う一対の断層フレームデータの自己相関値を用いることを特徴とする超音波診断装置。
【請求項４】
請求項2に記載の超音波診断装置において、
前記クオリティ算出部は、前記一対の断層フレームデータ間の自己相関値として、前記弾性ボリューム生成部で繰り返し生成される複数の弾性ボリュームデータをそれぞれ構成する複数の前記弾性フレームデータのうち、スキャン面位置が同一の前記弾性フレームデータを算出する基となった前記一対の断層フレームデータの自己相関値を用いることを特徴とする超音波診断装置。
【請求項５】
請求項2に記載の超音波診断装置において、
前記クオリティ算出部は、前記一対の断層フレームデータ間の自己相関値として、現在の断層フレームデータと過去の複数の断層フレームデータと間で求めた前記自己相関値が最大の一対の断層フレームデータの自己相関値を用いることを特徴とする超音波診断装置。
【請求項６】
請求項5に記載の超音波診断装置において、
前記自己相関値が最大の一対の断層フレームデータに基づいて算出された前記弾性フレームデータが収集され前記弾性ボリュームデータが生成されることを特徴とする超音波診断装置。
【請求項７】
請求項1に記載の超音波診断装置において、
前記クオリティ算出部は、前記弾性フレームデータの弾性値の分布の平均、偏差、又はS/N比に基づいて、各弾性フレームデータのクオリティを示すフレーム評価値を算出し、該算出したフレーム評価値を前記弾性ボリュームデータ単位で加算平均することにより前記弾性ボリュームデータのクオリティを示すボリューム評価値を算出することを特徴とする超音波診断装置。
【請求項８】
請求項7に記載の超音波診断装置において、
前記弾性値は、変位、歪み、歪み比、粘性、弾性率のいずれか1つであることを特徴とする超音波診断装置。
【請求項９】
請求項1に記載の超音波診断装置において、
前記弾性ボリュームデータのクオリティを示すボリューム評価値を前記3次元弾性画像に対応付けて前記表示部に表示することを特徴とする超音波診断装置。
【請求項１０】
請求項9に記載の超音波診断装置において、
前記弾性ボリュームデータのクオリティを示すボリューム評価値を、該ボリューム評価値に対応して表示態様が異なるマーク、バーチャート又は円グラフで前記表示部に表示することを特徴とする超音波診断装置。
【請求項１１】
請求項1に記載の超音波診断装置において、
前記弾性ボリュームデータ生成部は、前記弾性ボリュームデータに基づいて構成した直交3断面弾性画像又は平行な複数の断面でスライスしたマルチ断面弾性画像を前記表示部に表示する断面弾性画像構成部を備え、
前記クオリティ算出部は、前記各断面弾性画像のクオリティを示すボリューム評価値を算出して前記表示部に表示する表示面クオリティ算出部を備えてなることを特徴とする超音波診断装置。
【請求項１２】
請求項1に記載の超音波診断装置において、
前記弾性ボリュームデータのクオリティを示す前記ボリューム評価値が表示基準値よりも高い場合、前記表示部は前記3次元弾性画像を表示することを特徴とする超音波診断装置。
【請求項１３】
請求項1に記載の超音波診断装置において、
前記弾性ボリュームデータのクオリティを示すボリューム評価値が記憶基準値よりも高い場合、前記3次元弾性画像を記憶する記憶部を備えることを特徴とする超音波診断装置。
【請求項１４】
被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを複数収集して弾性ボリュームデータを生成するステップと、弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成するステップと、前記3次元弾性画像を表示するステップと、前記弾性フレームデータのクオリティを示すフレーム評価値に基づいて、前記弾性ボリュームデータのクオリティを示すボリューム評価値を算出するステップを含む評価算出方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波を利用して被検体の生体組織の硬さ又は軟らかさを示す3次元弾性画像を表示する超音波診断装置と評価算出方法に関する。
【背景技術】
【０００２】
超音波診断装置は、被検体に超音波を送信し、被検体内部の生体組織からの反射エコー信号を受信し、例えば、生体組織の硬さ又は軟らかさを示す3次元弾性画像などの超音波画像を生成してモニタに表示して診断に供するものである。また、被検体内部の形態を分かり易く表示して診断に寄与するために、3次元断層画像や3次元弾性画像の3次元超音波画像を生成して表示することが行われている。一方で、診断精度の向上のため、画像の分解能の向上、ノイズの低減、等々、超音波画像の画質を改善することが要請されている。
【０００３】
例えば、特許文献1には、弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構築する際に、同じ変位もしくは圧力の弾性画像区間を合成して3次元弾性画像を構築することが提案されている。さらに、同文献には、相関係数が高い弾性フレームで弾性画像を合成して、3次元弾性画像を構築することが提案されている。つまり、3次元弾性画像を構成する際に、連続して取得した複数の2次元弾性フレームデータから、圧迫量(変位)が同じ2次元弾性フレームデータにより弾性ボリュームデータを構成し、その弾性ボリュームデータをボリュームレンダリングすることにより、3次元弾性画像を構築するようにしている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開2008-259555号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、特許文献では、弾性ボリュームデータの中から所望の弾性フレームデータを選択して弾性ボリュームデータを合成しているが、取得した弾性ボリュームデータ全体についてのクオリティを評価することは考慮されていない。その結果、ボリュームレンダリングにより構築した3次元弾性画像にノイズ(例えば、スジ状のノイズ画像)が表れる可能性がある。
【０００６】
本発明が解決しようとする課題は、弾性ボリュームデータのクオリティを評価する手法を確立して、3次元弾性画像のクオリティを向上させることにある。
【課題を解決するための手段】
【０００７】
上記の課題を解決するため、本発明は、被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを算出する弾性演算部と、前記弾性フレームデータを複数収集して弾性ボリュームデータを生成する弾性ボリュームデータ生成部と、前記弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成する3次元弾性画像構成部と、前記3次元弾性画像を表示する表示部とを備えた超音波診断装置において、前記弾性フレームデータのクオリティを示すフレーム評価値に基づいて、前記弾性ボリュームデータのクオリティを示すボリューム評価値を算出するクオリティ算出部を備えたことを特徴とする。
【０００８】
以上のようにして求めたクオリティを示すボリューム評価値を、3次元弾性画像に対応付けて表示部に表示する。これにより、表示された3次元弾性画像のクオリティがどの程度高いか否かを判断することができる。また、クオリティを示すボリューム評価値に対応して表示態様が異なるマーク、バーチャート又は円グラフを表示部に表示することにより、一目でクオリティを示すボリューム評価値を識別できるから、診断の精度向上に寄与することができる。
【０００９】
さらに、本発明は、3次元弾性画像に限らず、弾性ボリュームデータに基づいて構成した直交3断面弾性画像又は平行な複数の断面でスライスしたマルチ断面弾性画像を前記表示部に表示する断面弾性画像構成部を備えて構成し、クオリティ算出部により各断面弾性画像のクオリティを示すボリューム評価値を算出して表示部に表示する表示面クオリティ算出部を備えて構成することができる。
【００１０】
以上説明したように、本発明によれば、3次元弾性画像のクオリティを表示することができるから、検査者が3次元弾性画像の選択を簡便に行うことができる。また、クオリティの高い3次元弾性画像を再構成して表示することもでき、高レベルな診断サポートが可能となる。
【発明の効果】
【００１１】
本発明によれば、弾性ボリュームデータのクオリティを評価する手法を確立できるから、3次元弾性画像のクオリティを向上させることができる。
【図面の簡単な説明】
【００１２】
【図１】本発明の実施例1の超音波診断装置の全体を示すブロック構成図
【図２】実施例1のクオリティ算出部のブロック構成図
【図３】実施例1の特徴部の処理手順を説明する図
【図４】本発明の実施例2の特徴部の処理手順を説明する図
【図５】本発明の実施例3の特徴部の処理手順を説明する図
【図６】本発明の実施例3の変形例の処理手順を説明する図
【図７】本発明の実施例4の断面弾性画像の表示例を説明する図
【図８】本発明の実施例4の断面弾性画像の表示を実現する主要部のブロック構成図
【図９】弾性ボリュームデータのクオリティを示すボリューム評価値の表示方法の一例を示す図
【図１０】弾性ボリュームデータのクオリティを示すボリューム評価値の表示方法の他の例を示す図
【発明を実施するための形態】
【００１３】
各図を用いて実施例について説明する。
【実施例１】
【００１４】
本発明を適用した一実施例の超音波診断装置100について、図1を用いて説明する。
図示のように、超音波診断装置100には、被検体101に当接させて用いる超音波探触子102と、超音波探触子102を介して被検体101に一定の時間間隔をおいて超音波を繰り返し送信させる送信部105と、被検体101の内部からの反射エコー信号を受信する受信部106と、送信部105と受信部106を制御する送受信制御部107と、受信部106で受信された反射エコーを整相加算する整相加算部108とが備えられている。
【００１５】
超音波探触子102は、被検体101に振動子を介して超音波を送受信する機能を有している。超音波探触子102は、矩形又は扇形をなす複数の振動子を配列してなり、複数の振動子の配列方向(長軸方向)と直交する方向(短軸方向)に超音波探触子102を機械的に振り、超音波を3次元スキャン可能に構成されている。なお、超音波の3次元スキャンは、超音波探触子102を機械的に短軸方向に振るものに限らず、例えば、複数の振動子を2次元に配列した超音波探触子102を用い、超音波を短軸方向に電子的にスキャンするものでもよい。
【００１６】
送信部105は、超音波探触子102の振動子を駆動して超音波を発生させるための送波パルスを生成する。送信部105は、送信される超音波の収束点を任意の深さに設定する機能を有している。また、受信部106は、超音波探触子102で受信した反射エコー信号を所定のゲインで増幅してRF信号、すなわち受信信号を生成する。超音波送受信制御部107は、送信部105や受信部106を制御するためのものである。整相加算部108は、受信部106で増幅されたRF信号の位相を制御して、1点又は複数の収束点に対応する超音波ビームを形成することにより、RF信号フレームデータを生成する。
【００１７】
整相加算部108で生成されたRF信号フレームデータはデータ記憶部109に記憶される。2次元断層画像構成部113は、データ記憶部109に記憶されたRF信号フレームデータに基づいて2次元断層画像を構成する。断層ボリュームデータ生成部114は、2次元断層画像構成部113で構成された2次元断層画像を取得位置に基づいて3次元座標変換を行なって断層ボリュームデータを生成する。3次元断層画像構成部115は、断層ボリュームデータの輝度と不透明度に基づいてボリュームレンダリングを行なって3次元断層画像を構成する。
【００１８】
データ記憶部109に記憶された複数のRF信号フレームデータは、2次元弾性画像構成部116に適宜取り込まれ、2次元弾性画像が構成される。2次元弾性画像構成部116で構成された2次元弾性画像は、弾性ボリュームデータ生成部117に取り込まれ、2次元弾性画像の取得位置に基づいて3次元座標変換を行なって弾性ボリュームデータが生成される。3次元弾性画像構成部118は、弾性ボリュームデータの弾性値と不透明度に基づいてボリュームレンダリングを行って3次元弾性画像を構成する。合成処理部119は、2次元断層画像と2次元弾性画像を合成したり、3次元断層画像と3次元弾性画像を合成するようになっている。表示部120は、合成処理部119で合成された合成画像や2次元断層画像の超音波画像を表示するようになっている。また、超音波診断装置100には、上述した各構成要素を制御する制御部103と、制御部103に各種入力を行なう操作部104を備え、操作部104は、キーボードやトラックボール等を備えている。
【００１９】
以下、図1の主要部の詳細構成について説明する。2次元断層画像構成部113は、制御部103における設定条件に基づいて、データ記憶部109から出力されるRF信号フレームデータを取り込んで、ゲイン補正、ログ圧縮、検波、輪郭強調、フィルタ処理等の信号処理を行なって2次元断層画像を構成する。このとき、超音波探触子102は、超音波の送受信と同時に送受信方向(θ、φ)を計測する。ここで、θは長軸方向のスキャン角度、φは短軸方向のスキャン(スイング)角度である。断層ボリュームデータ生成部114は、2次元断層画像の取得位置に相当する送受信方向(θ、φ)に基づいて、複数の2次元断層画像について3次元変換を行なって断層ボリュームデータを生成する。
【００２０】
3次元断層画像構成部115は、断層ボリュームデータから3次元断層画像を構成する次式(1)〜(3)を用いて、ボリュームレンダリングを行なう。
【００２１】
Cout(i)  =  Cout(i-1)＋(1−Aout(i-1))・A(i)・C(i)・S(i)   −(1)
Aout(i)  =  Aout(i-1)＋(1−Aout(i-1))・A(i)               −(2)
A(i)  =  Opacity[C(i)]                                    −(3)
ここで、C(i)は作成される2次元投影面上のある点から3次元断層画像を見た場合、視線上i番目に存在するボクセルの輝度値である。Cout(i)は、出力されるピクセル値である。例えば、視線上にNボクセルの輝度値が並んだとき、i =0〜Ｎ-1までを積算した輝度値Cout(N-1)が最終的に出力されるピクセル値となる。Cout(i-1)はi-1番目までの積算値を示す。
【００２２】
また、A(i)は視線上i番目に存在する輝度値の不透明度であり、式(3)に示すとおり、0〜1．0の値をとる断層不透明度テーブル(断層オパシティテーブル)である。断層不透明度テーブルは、輝度値から不透明度を参照することによって、出力する2次元投影面(3次元断層画像)上への寄与率を決定する。
【００２３】
S(i)は、輝度C(i)とその周辺の画素値より求めた勾配より算定される陰影付けのための重み成分で、たとえば、光源とボクセルiを中心とした面の法線が一致する場合、もっとも強く反射するため1．0が与えられ、光源と法線が直交する場合には0．0が与えられるなどの強調効果を指し示す。
【００２４】
Cout(i)、Aout(i)は、共に0を初期値としている。式(2)に示されるように、Aout(i)はボクセルを通過するたびに積算され1．0に収束される。よって、式(1)に示されるようにi-1番目までの不透明度の積算値Aout(i-1)がおよそ1.0となった場合、ｉ番目以降の輝度値C(i)は出力画像に反映されない。
【００２５】
2次元弾性画像構成部116は、データ記憶部109に記憶された複数のRF信号フレームデータから変位を計測する。2次元弾性画像構成部116は被検体101に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを算出する弾性演算部を有している。そして、2次元弾性画像構成部116は、計測した変位に基づいて弾性値を演算し、2次元弾性画像を構成する。弾性値とは、歪み、弾性率、変位、粘性、歪み比等の弾性情報のいずれかである。弾性ボリュームデータ生成部117は、2次元弾性画像の取得位置に相当する送受信方向(θ、φ)に基づいて、複数の2次元弾性画像について3次元変換を行ない、弾性ボリュームデータを生成する。3次元弾性画像構成部118は、弾性値に基づいて弾性ボリュームデータを複数に分離し、分離された弾性ボリュームデータについてボリュームレンダリングを行ない、3次元弾性画像を構成する。
【００２６】
以下、本発明の特徴部であるクオリティ算出部121について説明する。図2に、クオリティ算出部121のブロック構成図を示す。クオリティ算出部121は、基本的に、圧迫の情報もしくは3次元弾性画像の質であるクオリティ(Quality)を検出し、弾性ボリュームデータの質を演算して良し悪しを評価するものである。ここで、検出するクオリティとしては、2次元弾性画像を算出した際の一対の断層フレームデータの自己相関、変位、歪み、歪みのS/N比、等々を用いることができる。また、取得時間が異なる複数の弾性ボリュームデータ間の同一位置の弾性フレームデータの相関を用いることもできる。
【００２７】
クオリティ算出部121は、フレーム相関処理部201、フレーム変位・歪み処理部203、ボリューム処理部205、クオリティ算出部207により構成される。フレーム相関処理部201は、2次元弾性画像構成部116で弾性フレームデータを算出した際に用いた一対の断層フレームデータ間の自己相関を保存しておき、フレーム単位の弾性フレームデータのクオリティを示すフレーム評価値を算出するようになっている。
【００２８】
また、フレーム変位・歪み処理部203は、2次元弾性画像構成部116で算出された弾性値(歪み、弾性率、変位、粘性、歪み比)と圧力を弾性フレームデータに対応付けて保存し、フレーム単位の弾性フレームデータのクオリティを示すフレーム評価値を算出するようになっている。
【００２９】
ボリューム処理部205は、フレーム相関処理部201、フレーム変位・歪み処理部203において、フレーム単位で算出されたフレーム評価値を加算平均等により、弾性ボリュームデータのクオリティを示すボリューム評価値を求める。ここで、ボリューム評価値は1つあればよいが、複数のボリューム評価値を求めることもできる。
【００３０】
図3を参照して、本実施例1の詳細な処理手順について説明する。図3(a)に示すように、3次元スキャンは、長軸方向(θ)に電子スキャンを行って断層フレームデータを計測する操作を、短軸方向(φ)にスキャン断面をスライドさせながら複数(n)フレームの断層フレームデータFr.0―Fr.nを取得する。短軸方向(φ)のスキャンは、機械的なスキャン又は電子的なスキャンのいずれを適用してもよい。例えば、モータにより超音波探触子102を遥動させてスキャンする場合は、モータを低速又は高速に速度制御可能にする。
【００３１】
3次元スキャンにおいて2次元弾性画像構成部116は、例えば、隣り合う一対の断層フレームデータFr.0、Fr.1間の自己相関演算を行って生体組織の変位分布を推定する。なお、図示していないが、図1の断層ボリュームデータ生成部114に記憶されている異なる時間に取得された同一部位の複数の断層ボリュームデータを参照して、現在(リアルタイム)取得された例えば断層フレームデータFr.0と、同一位置において過去に取得された断層フレームデータFr.0との間で自己相関演算を行って生体組織の変位分布を推定することもできる。このように、一対の断層フレームデータ間で求められた自己相関値が、その一対の断層フレームデータにより算出された弾性フレームデータのクオリティを示すフレーム評価値となる。なお、一対の断層フレーム間の自己相関値が高いということは、断層フレームデータを構成する複数本のRF信号の波形の一致度が高いことになり、これにより得られる弾性フレームデータのクオリティは高い傾向にある。
【００３２】
そこで、フレーム相関処理部201は、2次元弾性画像構成部116から各対の断層フレーム間の自己相関値を取り込み、図3(a)、(b)に示すように、各弾性フレームデータについて画素単位で自己相関値のフレーム平均値Caveを求める。次いで、弾性ボリュームデータを構成する全部の弾性フレームデータに対応する自己相関値をフレーム評価値とし、その加算平均又は加算値を求めて、弾性ボリュームデータのクオリティを示すボリューム評価値を求める。例えば、自己相関値の弾性ボリュームの加算平均Cvol=0．95であればクオリティを示すボリューム評価値が高いとして、図3(c)のように「Quality High」が表示部120に表示される。また、その弾性ボリュームデータをボリュームレンダリングして得られる3次元弾性画像は、崩れが少ないものとなる。一方、図3(d)のように自己相関値の加算平均Cvol=0．74であれば「Quality Low」が表示部120に表示される。この場合の3次元弾性画像は、ノイズが多く、例えば本来なら円形の画像が変形して表示されるなどの先鋭さに欠ける画像となる。このようにして、弾性ボリュームデータごとにクオリティを示すボリューム評価値を求めることにより、弾性ボリュームデータのクオリティを評価することが可能である。なお、フレーム変位・歪み処理部203の詳細な実施例は、後述する。
【実施例２】
【００３３】
図4を参照して、本発明の実施例2の2次元弾性画像構成部116とクオリティ算出部121の詳細な処理手順について説明する。本実施例2が実施例1と異なる点は、2次元弾性画像構成部116とフレーム相関処理部201の処理手順にあり、その他の点は実施例1と同様であることから説明を省略する。実施例1においては、2次元弾性画像構成部116は隣り合う一対の断層フレームデータ間の自己相関演算を行って生体組織の変位分布を推定していた。これに対し、本実施例2では、2次元弾性画像構成部116は現在(リアルタイム)の断層フレームデータに対する過去数フレーム分の断層フレームデータとの間で自己相関演算を行って生体組織の変位分布を推定する。特に、過去数フレーム分との間の自己相関値が最も高い断層フレームデータを選択して、生体組織の変位分布を推定して弾性フレームデータを作成する。つまり、複数回の3次元スキャンで1フレーム分の自己相関値と変位分布データが算出される。その後の処理は、実施例1と同様であり、図4(c)、(d)に示す評価結果が表示される。
【００３４】
本実施例2によれば、複数回の3次元スキャンで自己相関値が高い弾性フレームデータを選択できることから、最終的に質の高い弾性ボリュームデータを生成することができ、3次元弾性画像のクオリティを高めることができる。
【実施例３】
【００３５】
図5を参照して、本発明の実施例3のクオリティ算出部121の詳細な処理手順について説明する。本実施例3が実施例1，2と異なる点は、弾性フレームデータを算出した際の一対の断層フレームデータの自己相関値に代えて、弾性フレームデータの弾性値(変位、歪み、歪み比、粘性、弾性率)の平均、S/N比、偏差情報の少なくとも1つに基づいて、弾性ボリュームデータを構成する弾性フレームデータのクオリティの評価値を算出するようにしたことにある。
【００３６】
図5(a)に示すように、短軸方向に2次元スキャンを繰り返して、フレーム間自己相関から検出された弾性(変位)フレームデータを複数(n)フレーム取得する。これにより、同図(b)に示すような2次元の弾性フレームデータの集合体として、3次元の弾性ボリュームデータが作成される。各弾性フレームデータの変位の例えば加算平均を求めて、弾性ボリュームデータのクオリティを示すボリューム評価値とする。また、各フレームの変位の加算平均の分散を用いることにより、画像のバラつきを検出できる。各フレームの加算平均の分散Saveは、式(4)により算出する。さらに、式(5)により、弾性ボリュームデータを構成する全ての弾性フレームデータの分散Saveの弾性ボリュームデータ平均Svolを求めて、弾性ボリュームデータのクオリティを示すボリューム評価値とする。
2012029417.tif 000015
ここで、Ｗは画像の横幅、Ｈは画像の高さ、d(i,j)は変位、d*は変位の平均である。
2012029417.tif 000016
本実施例によれば、バラつきが少ない方がノイズは少ないから、図5(c)に示すように、Svol＝0.001のように小さな値の方が評価が高いボリューム評価値「Quality High」と表示され、3次元弾性画像としても引き締まった画像となる。一方、ノイズの多い画像においては、バラつきが大きくなり、図5(d)に示すように、Svol＝0.02のように大きな値の方が評価が低いボリューム評価値「Quality Low」と表示され、画像も先鋭さを欠ける画像となる。
【００３７】
図6は、本実施例の変形例であり、変位フレームデータに代えて、歪みフレームデータに基づいて弾性ボリュームデータのクオリティを評価するようにした例である。図6(a)に示すように、複数(n)フレームの歪みフレームデータにより弾性ボリュームデータが構成される。一般に、生体組織に圧迫が加えられていれば、負の歪みは存在しえない。したがって、図6(b)に示すように、一の歪みフレームデータ中に負の歪みが混在していると、歪みの平均が小さくなる。そこで、各フレームの歪みの平均をφ方向の各歪みフレームデータについて加算平均を求めて、弾性ボリュームデータのクオリティを示すボリューム評価値とする。これによれば、図6(c)に示すように評価が高い方が「Quality High」の表示になり、評価が低い場合は図6(d)に示すように「Quality Low」の表示になる。
【実施例４】
【００３８】
図7に、本発明の各実施例により求めた弾性ボリュームデータを用いて構成される各種の表示画像の例を示す。図7の例では、弾性ボリュームデータを用いて任意の断面における断面画像を表示する場合において、表示された各断面画像にそれぞれクオリティを示すボリューム評価値を表示するようにしている。具体的には、図7(a)に示すように、任意の直交3断面のアキシャル面A、シータ面T，コロナル面C、3次元弾性画像Renderのそれぞれに、クオリティを示すボリューム評価値を表示する。また、図7(b)に示すように、いわゆるマルチスライス断面画像の場合も、複数断面画像の表示にそれぞれに、クオリティを示すボリューム評価値を表示することができる。
【００３９】
図7の画像表示例を実現するには、弾性ボリュームデータ生成部117とクオリティ算出部121を図8のように構成する。つまり、弾性ボリュームデータ生成部117にて、2次元弾性画像構成部116で算出された自己相関値は、弾性ボリュームデータ生成部117において、自己相関値ボリューム保持部802又は変位・歪みボリューム保持部804に格納される。このデータは、超音波のRF配列データ曲座標(RTΦ)次元のデータであり、直交座標変換部806にて直交座標のデータに変換される。制御部103において、表示用に参照される座標面に対して表示面クオリティ算出部808で各断面のクオリティを示すボリューム評価値を算出する。
【実施例５】
【００４０】
図9及び図10を参照して、各実施例で評価した弾性ボリュームデータを示すボリューム評価値の表示方法の実施例について説明する。図9は、横軸が時間、縦軸がクオリティを示すボリューム評価値を示しており、連続的に取得した複数の弾性ボリュームデータV0−Vnのボリューム評価値が一目でわかる。また、最もボリューム評価値が高い弾性ボリュームデータと、最もボリューム評価値が低い弾性ボリュームデータに対応付けて、図3と同様のマークを表示する。
【００４１】
図10は、弾性ボリュームデータ又は3次元弾性画像のクオリティを示すボリューム評価値を表すマーク、バーチャート、円グラフ、その他の変形例を示している。図10(a)は、円形のマークの表示色をボリューム評価値に対応付けて変えることにより、一目でボリューム評価値の高低を識別可能にした例である。同図(b)は、実施例1等で示したマークと同じであり、円形のマークの真円度合いに応じてボリューム評価値の高低を表した例である。同図(c)は、バーチャートのゲージ内の比率に応じてボリューム評価値の高低を表した例である。同図(d)は、円グラフの比率に応じてボリューム評価値の高低を表した例である。同図(e)は、複数の小円を並べてボリューム評価値に応じて、小円の表示態様の比率を異ならせた例である。
【実施例６】
【００４２】
本発明の各実施例により求めた弾性ボリュームデータを用いて構成される3次元弾性画像の表示形態と記憶形態を説明する。
【００４３】
クオリティ算出部121によって算出された弾性ボリュームデータのクオリティを示すボリューム評価値に基づいて、表示部120は3次元弾性画像を表示する。ボリューム評価値が表示基準値よりも高い場合、制御部103は3次元弾性画像を表示するように指示を行い、表示部120はボリューム評価値が表示基準値よりも高い3次元弾性画像を表示する。ボリューム評価値が表示基準値よりも低い場合、制御部103は3次元弾性画像を表示しないように指示を行い、表示部120はボリューム評価値が表示基準値よりも低い3次元弾性画像を表示しない。
【００４４】
表示基準値は、例えば、0.95とする。また、操作者は操作部104で表示基準値を設定することができる。
【００４５】
本実施例によれば、弾性ボリュームデータのクオリティを示すボリューム評価値が表示基準値よりも高い場合、表示部120は3次元弾性画像を表示する。つまり、表示部120は弾性ボリュームデータのクオリティを示すボリューム評価値が表示基準値よりも高い3次元弾性画像のみを表示することができる。
【００４６】
また、クオリティ算出部121によって算出された弾性ボリュームデータのクオリティを示すボリューム評価値に基づいて、記憶部(図示しない。)は3次元弾性画像を記憶する。ボリューム評価値が記憶基準値よりも高い場合、制御部103は3次元弾性画像を記憶部に記憶するように指示を行い、記憶部はボリューム評価値が記憶基準値よりも高い3次元弾性画像を記憶する。ボリューム評価値が記憶基準値よりも低い場合、制御部103は3次元弾性画像を記憶部に記憶しないように指示を行い、記憶部はボリューム評価値が記憶基準値よりも低い3次元弾性画像を記憶しない。
【００４７】
本実施例によれば、弾性ボリュームデータのクオリティを示すボリューム評価値が記憶基準値よりも高い場合、3次元弾性画像を記憶する記憶部を備える。つまり、記憶部は弾性ボリュームデータのクオリティを示すボリューム評価値が記憶基準値よりも高い3次元弾性画像のみを記憶することができる。
【００４８】
記憶基準値は、表示基準値と同様に設定することができ、操作者は操作部104で記憶基準値を設定することができる。表示基準値と記憶基準値を同一にすることもできる。
【００４９】
以上、本発明によれば、被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを複数収集して弾性ボリュームデータを生成する弾性ボリュームデータ生成部117と、弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成する3次元弾性画像構成部118と、3次元弾性画像を表示する表示部120とを備えた超音波診断装置において、弾性フレームデータのクオリティを示すフレーム評価値基づいて、弾性ボリュームデータのクオリティを示すボリューム評価値を算出するクオリティ算出部121を備えたことを特徴とする。また、被検体に超音波を3次元スキャンして計測される弾性分布を表す弾性フレームデータを複数収集して弾性ボリュームデータを生成するステップと、弾性ボリュームデータをボリュームレンダリングして3次元弾性画像を構成するステップと、前記3次元弾性画像を表示するステップと、前記弾性フレームデータのクオリティを示すフレーム評価値に基づいて、前記弾性ボリュームデータのクオリティを示すボリューム評価値を算出するステップを含む評価算出方法である。
【００５０】
すなわち、弾性フレームデータのクオリティを示すフレーム評価値が高いことは、安定した計測状態で弾性フレームデータが算出されたことになる。これに鑑み、本発明は、弾性フレームデータのクオリティを評価するフレーム評価値を求め、さらに弾性ボリュームデータを構成する複数の弾性フレームデータを示すフレーム評価値に基づいて、弾性ボリュームデータのクオリティを示すボリューム評価値を求めて評価するようにしたのである。このように、弾性ボリュームデータのクオリティを評価する手法が確立されたことから、クオリティを示すボリューム評価値が高い弾性ボリュームデータを取得又は選択することにより、容易に3次元弾性画像のクオリティを向上させることができる。
【００５１】
ここで、弾性フレームデータのクオリティあるいはクオリティ評価値とは、安定した適切な圧迫状態において弾性フレームデータが計測されたことを意味する。同様に、弾性ボリュームデータのクオリティあるいはクオリティ評価値とは、クオリティあるいはクオリティ評価値が高い弾性フレームデータの集合であることを意味する。したがって、クオリティあるいはクオリティ評価値が高いことは、結果としてノイズの少ない3次元弾性画像を生成することができる。
【００５２】
ところで、弾性フレームデータのクオリティを示すフレーム評価値は、クオリティ評価方法により求めることができる。例えば、弾性フレームデータ算出の基となった一対の断層フレームデータ間の自己相関値が高いと、一対の断層フレームデータの一致度が高いことを意味し、その一対の断層フレームデータが安定した計測状態で計測されたことになる。そこで、弾性フレームデータ算出の基となった一対の断層フレームデータの自己相関値をフレーム評価値とし、これに基づいて弾性ボリュームデータのクオリティを示すボリューム評価値を求めて評価することができる。しかし、本発明に係るフレーム評価値は、後述するように、一対の断層フレームデータの自己相関値に限られるものではない。
【００５３】
本発明において、クオリティ算出部121は、一対の断層フレームデータ間の自己相関値を弾性フレームデータのクオリティを示すフレーム評価値とし、弾性ボリュームデータを構成する全ての弾性フレームデータのフレーム評価値の加算値又は加算平均値に基づいて、弾性ボリュームデータのクオリティを示すボリューム評価値を求めるようにすることができる。この場合、一対の断層フレームデータ間の自己相関値として、3次元スキャンにおける時間的に隣り合う一対の断層フレームデータの自己相関値を用いることができる。また、これに代えて、一対の断層フレームデータ間の自己相関値として、弾性ボリュームデータ生成部で繰り返し生成される複数の弾性ボリュームデータをそれぞれ構成する複数の弾性フレームデータのうち、スキャン面位置が同一の前記弾性フレームデータ算出の基となった前記一対の断層フレームデータの自己相関値を用いることができる。
【００５４】
さらに、一対の断層フレームデータ間の自己相関値として、現在の断層フレームデータと過去の複数の断層フレームデータとの間で求めた自己相関値が最大の一対の断層フレームデータの自己相関値を用いることができる。これによれば、クオリティの高い弾性フレームデータを選択して弾性ボリュームデータを構成できることから、一層、ボリューム評価値を高めることができる。つまり、自己相関値が最大の一対の断層フレームデータに基づいて算出された弾性フレームデータを収集して弾性ボリュームデータが生成される。なお、前述したように、本発明は、弾性フレームデータ算出の基となった一対の断層フレームデータの自己相関値に基づいてフレーム評価値を求めるものに限られない。これに代えて、弾性フレームデータの弾性値の分布の平均、偏差、又はS/N比に基づいて、各弾性フレームデータのクオリティを示すフレーム評価値を算出することができる。例えば、弾性値のS/N比が大きければ弾性フレームデータのクオリティが高いと評価することができる。また、弾性値としては、変位、歪み、歪み比、粘性、弾性率のいずれか1つを用いることができる。
【００５５】
以上のようにして求めたクオリティを示すボリューム評価値を、3次元弾性画像に対応付けて表示部120に表示することが望ましい。これにより、表示された3次元弾性画像のクオリティがどの程度高いか否かを判断することができる。また、クオリティを示すボリューム評価値に対応して表示態様が異なるマーク、バーチャート又は円グラフを表示部に表示することにより、一目でクオリティを示すボリューム評価値を識別できるから、診断の精度向上に寄与することができる。
【００５６】
さらに、本発明は、3次元弾性画像に限らず、弾性ボリュームデータに基づいて構成した直交3断面弾性画像又は平行な複数の断面でスライスしたマルチ断面弾性画像を前記表示部に表示する断面弾性画像構成部を備えて構成し、クオリティ算出部により各断面弾性画像のクオリティを示すボリューム評価値を算出して表示部に表示する表示面クオリティ算出部を備えて構成することができる。
【００５７】
以上説明したように、本発明によれば、3次元弾性画像のクオリティを表示することができるから、検査者が3次元弾性画像の選択を簡便に行うことができる。また、クオリティの高い3次元弾性画像を再構成して表示することもでき、高レベルな診断サポートが可能となる。
【符号の説明】
【００５８】
100  超音波診断装置、102  超音波探触子、103  制御部、104  操作部  、105  送信部、106  受信部、107  送受信制御部、108  整相加算部、109  データ記憶部、113  2次元断層画像構成部、114  断層ボリュームデータ生成部、115  3次元断層画像構成部、116  2次元弾性画像構成部、117  弾性ボリュームデータ生成部、118  3次元弾性画像構成部、119  合成処理部、120  表示部、121  クオリティ算出部
【手続補正３】
【補正対象書類名】図面
【補正対象項目名】図１
【補正方法】変更
【補正の内容】
【図１】
2012029417.tif 000017
【手続補正４】
【補正対象書類名】図面
【補正対象項目名】図２
【補正方法】変更
【補正の内容】
【図２】
2012029417.tif 000018
【手続補正５】
【補正対象書類名】図面
【補正対象項目名】図８
【補正方法】変更
【補正の内容】
【図８】
2012029417.tif 000019
【国際調査報告】
2012029417.tif 000020
2012029417.tif 000021
2012029417.tif 000022
2012029417.tif 000023