(57)【要約】  操作者に依存せず、安定した計測位置で弾性情報の指標値の演算を行うことができる超音波診断装置及び超音波計測方法を提供する。  超音波探触子12を介して被検体の診断部位の断層画像を構成する断層画像構成部24と、硬さを示す弾性情報を演算する弾性情報演算部34と、弾性情報演算部34において演算された弾性情報に基づいて弾性画像を構成する弾性画像構成部36と、断層画像と弾性画像を表示する画像表示部28とを備えた超音波診断装置において、弾性情報演算部34で演算された弾性情報から各計測点における時間ばらつきを解析する時間ばらつき解析部50と、時間ばらつき解析部で解析された時間ばらつきに基づいて、時間ばらつき画像を構成する時間ばらつき画像構成部52とを備え、画像表示部28は時間ばらつき画像を表示する。
2014061370.tif 000002
【特許請求の範囲】
【請求項１】
超音波探触子を介して被検体の診断部位の断層画像を構成する断層画像構成部と、硬さを示す弾性情報を演算する弾性情報演算部と、前記弾性情報演算部において演算された前記弾性情報に基づいて弾性画像を構成する弾性画像構成部と、前記断層画像と前記弾性画像を表示する画像表示部とを備えた超音波診断装置において、
前記弾性情報演算部で演算された前記弾性情報から各計測点における時間ばらつきを解析する時間ばらつき解析部と、前記時間ばらつき解析部で解析された時間ばらつきに基づいて、時間ばらつき画像を構成する時間ばらつき画像構成部とを備え、前記画像表示部は前記時間ばらつき画像を表示することを特徴とする超音波診断装置。
【請求項２】
前記時間ばらつき画像構成部は、前記時間ばらつきが閾値より大きい計測点については前記時間ばらつき画像を構成し、前記時間ばらつきが閾値より小さい計測点については前記時間ばらつき画像を構成しないことを特徴とする請求項１記載の超音波診断装置。
【請求項３】
前記時間ばらつき解析部は、複数の弾性情報のフレームを記憶する弾性情報フレーム記憶部と、複数の弾性情報のフレームから時間ばらつきを算出する時間ばらつき算出部を有していることを特徴とする請求項１記載の超音波診断装置。
【請求項４】
前記時間ばらつき算出部は、弾性情報の比、差、標準偏差、分散、変動係数のいずれか１つを用いて前記時間ばらつきを算出することを特徴とする請求項３記載の超音波診断装置。
【請求項５】
前記時間ばらつき算出部は、前記複数の弾性情報のフレームにおいて、前記断層画像と前記弾性画像の合成画像に設定された関心領域と各計測点における弾性情報の指標値を演算し、各計測点における弾性情報の指標値のばらつきを算出することを特徴とする請求項３記載の超音波診断装置。
【請求項６】
前記時間ばらつき算出部は、前記複数の弾性情報のフレームにおいて、各計測点における弾性情報の比の上限値と下限値の差を算出し、前記弾性情報の比の上限値と下限値の差が予め定められた閾値より大きい場合、時間ばらつきが大きいと算出することを特徴とする請求項３記載の超音波診断装置。
【請求項７】
前記時間ばらつき画像構成部は、前記弾性情報の比の上限値と下限値の差が予め定められた閾値より大きい場合、前記時間ばらつき画像を構成して前記画像表示部に表示させることを特徴とする請求項３記載の超音波診断装置。
【請求項８】
前記時間ばらつき算出部は、前記弾性情報の変位が予め定められた閾値より大きい場合、時間ばらつきが大きいと算出することを特徴とする請求項３記載の超音波診断装置。
【請求項９】
前記画像表示部は、前記断層画像と前記時間ばらつき画像の合成画像を表示することを特徴とする請求項１記載の超音波診断装置。
【請求項１０】
前記時間ばらつき画像が所定領域以上である場合、撮影し直すアラームを前記画像表示部に表示させる制御部を備えることを特徴とする請求項１記載の超音波診断装置。
【請求項１１】
前記時間ばらつき画像が表示された領域を認識し、前記時間ばらつき画像が表示された領域以外に関心領域を設定する制御部を備えることを特徴とする請求項１記載の超音波診断装置。
【請求項１２】
第１の関心領域と第２の関心領域を予め設定しておき、前記時間ばらつき画像構成部によって構成された前記時間ばらつき画像が前記第１の関心領域と前記第２の関心領域に重なる場合、前記弾性情報演算部によって前記第１の関心領域と前記第２の関心領域における弾性情報の指標値を演算しないことを特徴とする請求項１記載の超音波診断装置。
【請求項１３】
前記弾性情報演算部における前記弾性情報の指標値の演算を行う演算時間を設定する演算時間設定部を備えることを特徴とする請求項１記載の超音波診断装置。
【請求項１４】
断層画像と弾性画像を表示する画像表示方法において、
硬さを示す弾性情報から各計測点における時間ばらつきを解析し、
時間ばらつきに基づいて、時間ばらつき画像を構成し、
前記時間ばらつき画像を表示することを特徴とする画像表示方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波を利用して被検体内の断層画像を表示する超音波診断装置及び画像表示方法に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波探触子により被検体内部に超音波を送信し、被検体内部から生体組織の構造に応じた超音波の反射エコー信号を受信し、被検体内の断層画像を構成して表示する。【０００３】
また、被検体の断層部位における組織の歪みを求め、その歪みに基づいて断層部位における弾性画像を生成して表示装置に表示し、表示装置に表示された断層画像若しくは弾性画像に複数の関心領域を設定し、各関心領域における弾性情報を指標値化し、該指標値を表示装置に表示することが行われている(特許文献1参照)。【０００４】
特許文献1によれば、診断部位の生体組織の硬さを定量的に評価することができる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】国際公開第WO06/013916号
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、操作者が複数の関心領域を断層画像若しくは弾性画像に手動で設定していたため、算出された弾性情報の指標値は操作者の設定に依存してしまうことになる。また、操作者は安定した計測位置に関心領域を設定して、計測を行っているかどうか確認することができなかった。【０００７】
そこで、本発明は、弾性情報の時間ばらつきに基づく時間ばらつき画像を表示することができる超音波診断装置及び画像表示方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記の課題を解決するため、超音波探触子を介して被検体の診断部位の断層画像を構成する断層画像構成部と、硬さを示す弾性情報を演算する弾性情報演算部と、前記弾性情報演算部において演算された前記弾性情報に基づいて弾性画像を構成する弾性画像構成部と、前記断層画像と前記弾性画像を表示する画像表示部とを備えた超音波診断装置において、前記弾性情報演算部で演算された前記弾性情報から各計測点における時間ばらつきを解析する時間ばらつき解析部と、前記時間ばらつき解析部で解析された時間ばらつきに基づいて、時間ばらつき画像を構成する時間ばらつき画像構成部とを備え、前記画像表示部は前記時間ばらつき画像を表示する。【０００９】
また、断層画像と弾性画像を表示する画像表示方法において、硬さを示す弾性情報から各計測点における時間ばらつきを解析するステップと、時間ばらつきに基づいて、時間ばらつき画像を構成するステップと、前記時間ばらつき画像を表示するステップとを有する。
【発明の効果】
【００１０】
本発明によれば、弾性情報の時間ばらつきに基づく時間ばらつき画像を表示することができる。
【図面の簡単な説明】
【００１１】
【図１】本発明の超音波診断装置の構成を例示するブロック図。
【図２】本発明の時間ばらつき解析部の構成を例示するブロック図。
【図３】本発明の画像表示部の表示形態を例示する図。
【図４】本発明の画像表示部の表示形態を例示する図。
【図５】本発明の動作を示すフローチャート。
【図６】本発明の超音波診断装置の構成を例示するブロック図。
【図７】本発明の画像表示部の表示形態を例示する図。
【発明を実施するための形態】
【００１２】
以下、本発明の超音波診断装置について、図面を参照して説明する。
【００１３】
(第1の実施形態)
図1は、本発明の第1の実施形態の超音波診断装置を例示するブロック図である。
図1に示すように、本実施形態に係る超音波診断装置は、被検体10に当接させて用いる超音波探触子12と、超音波探触子12を介して被検体10に時間間隔をおいて超音波を繰り返し送信する送信部14と、被検体10から発生する時系列の反射エコー信号を受信する受信部20と、送信部14と受信部20を制御する超音波送受信制御部18と、受信された反射エコーを整相加算してRF信号フレームデータを時系列に生成する整相加算部22と、整相加算部22で生成されたRF信号フレームデータに基づいて断層画像を構成する断層画像構成部24と、断層画像と他の画像、数値情報などを合成する画像合成部26と、画像合成部26から出力される画像を表示する画像表示部28と、少なくとも2枚のRF信号フレームデータを選択するRF信号フレームデータ選択部30と、選択されたRF信号フレームデータを用いて被検体10の生体組織の変位を計測する変位計測部32と、変位計測部32で計測された変位から弾性情報を求める弾性情報演算部34と、弾性情報演算部34で演算した弾性情報から弾性画像を構成する弾性画像構成部36と、操作者が操作するための操作部40と、操作部40の操作に応じて各構成要素を制御する制御部42と、弾性情報演算部34で演算した弾性情報から時間ばらつきを解析する時間ばらつき解析部50と、時間ばらつき解析部50で解析された時間ばらつきに基づいて時間ばらつき画像を構成する時間ばらつき画像構成部52とを備えている。図1で示す破線は、超音波診断装置の本体を示すものである。
【００１４】
超音波探触子12は、複数の振動子を配設して形成されており、接触させた被検体10に振動子を介して超音波を送受信する。送信部14は、超音波探触子12を駆動して超音波を発生させるための送波パルスを生成するとともに、送信される超音波の収束点をある深さに設定し、超音波探触子12を介して被検体10に時間間隔をおいて超音波を繰り返し送信する。受信部20は、被検体10から発生する時系列の反射エコー信号を超音波探触子12を介して受信し、受信した反射エコー信号について所定のゲインで増幅してRF信号(受波信号)を生成する機能を有している。超音波送受信制御部18は、送信部14及び受信部20を制御し、超音波探触子12を介して被検体10に対して超音波を送受信させる。整相加算部22は、受信部20で受信された反射エコー信号を整相加算する。その際、整相加算部22は、受信部20で増幅されたRF信号を入力して位相制御し、一点又は複数の収束点に対し超音波ビームを形成して超音波断層データであるRF信号フレームデータを時系列に生成する。
【００１５】
断層画像構成部24は、被検体10の断層部位のデータ、具体的には整相加算部22からのRF信号フレームデータを入力してゲイン補正、ログ圧縮、検波、輪郭強調、フィルタ処理等の信号処理を行い、断層画像データ(例えば、被検体10の白黒の濃淡断層画像)を構成する。また、断層画像構成部24は、図示はしないが、断層画像データをデジタル信号に変換するA/D変換器と、変換された複数の断層画像データを時系列に記憶するフレームメモリと、制御コントローラを含んで構成されている。フレームメモリに記憶された被検体10内の断層画像データが1画像として取得され、取得された断層画像データがテレビ同期で読み出される。
【００１６】
RF信号フレームデータ選択部30は、整相加算部22から出力されるRF信号フレームデータを記憶し、記憶されたRF信号フレームデータ群から少なくとも2枚(一組)のRF信号フレームデータを選択する。例えば、RF信号フレームデータ選択部30は、整相加算部22から時系列、すなわち画像のフレームレートに基づいて生成されるRF信号フレームデータを順次記憶し、記憶されたRF信号フレームデータ(N)を第1のデータとして選択すると同時に、時間的に過去に記憶されたRF信号フレームデータ群(N-1、N-2、N-3・・・N-M)の中から1つのRF信号フレームデータ(X)を選択する。なお、N、M、Xは、RF信号フレームデータに付されたインデックス番号であり、自然数とする。
【００１７】
変位計測部32は、被検体10の生体組織の変位を計測する。具体的に説明すると、変位計測部32は、RF信号フレームデータ選択部30により選択された1組のデータ、すなわちRF信号フレームデータ(N)及びRF信号フレームデータ(X)から1次元あるいは2次元相関処理を行って、断層画像の各計測点に対応する生体組織における変位を示すベクトル、すなわち変位の方向と大きさに関する1次元又は2次元変位分布を求める。ここで、ベクトルの検出にはブロックマッチング法もしくは位相勾配法を用いる。ブロックマッチング法においては、画像を例えばN×N画素からなるブロックに分け、所定領域(例えば、後述するパラメータ取得領域)内のブロックに着目し、現フレーム中の着目しているブロックに最も近似しているブロックを前のフレームから探し、これを参照して予測符号化、すなわち差分により標本値を決定する処理を行う。これにより、断層画像の各計測点の変位を求め、ベクトルの検出を行う。位相勾配法においては、受信信号の波の位相情報からその波の移動量を算出して断層画像の各計測点の変位を求め、ベクトルの検出を行う。
【００１８】
図示しない圧力計測部は、超音波探触子12の超音波送受信面と被検体10との間に設けられた圧力センサ等により検出された圧力に基づいて、被検体10内部の計測点における応力を計測する。
【００１９】
弾性情報演算部34は、被検体10の断層部位の超音波断層データに基づいて、断層部位における組織の歪み又は弾性率を求める。本実施形態において、弾性情報演算部34は、RF信号フレームデータ選択部30により選択されたRF信号フレームデータを用いて変位計測部32で計測された生体組織の変位情報、例えば移動ベクトルに基づいて断層画像上の各計測点に対応する生体組織の歪みや弾性率を演算する。なお、生体組織の弾性率を演算する際、弾性情報演算部34は、圧力計測部から出力される圧力値も加味している。
【００２０】
このとき、歪みのデータは、生体組織の移動量、例えば変位を空間微分することによって算出される。また、弾性率のデータは、圧力の変化を歪みの変化で除することによって算出される。例えば、変位計測部32により計測された変位をL(X)、圧力計測部により計測された圧力をP(X)とすると、歪みΔS(X)は、L(X)を空間微分することによって算出することができるから、ΔS(X)＝ΔL(X)／ΔXという式を用いて求められる。
【００２１】
また、弾性率データのヤング率Ym(X)は、Ym＝ΔP(X)／ΔS(X)という式によって求められる。このヤング率Ymから断層画像の各計測点に相当する生体組織の弾性率が求められるので、2次元の弾性画像データを連続的に得ることができる。なお、ヤング率とは、物体に加えられた単純引張り応力と、引張りに平行に生じる歪みに対する比である。
【００２２】
弾性画像構成部36は、弾性情報演算部34で求めた弾性情報に基づいて断層部位における弾性画像を構成する。弾性画像構成部36は、フレームメモリと画像処理部とを含んで構成されており、弾性フレームデータをフレームメモリに記憶し、記憶されたフレームデータに対し画像処理を行うようになっている。
【００２３】
また、弾性画像構成部36は、弾性フレームデータに色相情報を付与する機能を有しており、弾性フレームデータに基づいて光の3原色である赤(R)、緑(G)、青(B)を付した画像データに変換する。例えば、弾性画像構成部36は、歪みが大きい弾性データを赤色コードに変換し、歪みが小さい弾性データを青色コードに変換する。
【００２４】
時間ばらつき解析部50は、弾性情報演算部34で演算した弾性情報から各計測点(各画素)における時間ばらつきを解析する。時間ばらつき解析部50は、複数の弾性情報のフレームから各計測点における時間ばらつきを解析し、時間ばらつきが閾値より大きい計測点を算出する。時間ばらつきが大きい計測点とは、複数の弾性情報のフレーム間における弾性情報の指標値が安定しない(ばらついてしまう)計測点であり、弾性情報の指標値による計測を行うことができない計測点である。弾性情報の指標値とは、弾性情報の比、弾性情報の差などで指標化された値である。
【００２５】
時間ばらつき画像構成部52は、時間ばらつき解析部50で解析された時間ばらつきに基づいて、時間ばらつき度合いを示す時間ばらつき画像を構成する。時間ばらつき画像構成部52は、時間ばらつきが閾値より大きい計測点については時間ばらつき画像を構成し、時間ばらつきが閾値より小さい計測点については時間ばらつき画像を構成しない。
【００２６】
画像合成部26は、フレームメモリと、画像処理部と、画像選択部を備えて構成されており、αブレンディングに代表されるような手法にて断層画像と弾性画像の合成画像、断層画像と時間ばらつき画像の合成画像を作成する。フレームメモリは、断層画像構成部24からの断層画像データと弾性画像構成部36からの弾性画像データと時間ばらつき画像構成部52からの時間ばらつき画像データを記憶する。
【００２７】
また、画像処理部は、フレームメモリに記憶された断層画像データと弾性画像データとを合成割合を変更して合成し、フレームメモリに記憶された断層画像データと時間ばらつき画像データとを合成割合を変更して合成する。合成画像の各画素の輝度情報及び色相情報は、それぞれ合成される画像の各情報を合成割合で加算したものとなる。さらに、画像選択部は、フレームメモリ内の断層画像データ、弾性画像データ、時間ばらつき画像データ、及び画像処理部の合成画像データのうちから表示する画像を選択し、画像表示部28に表示させる。
【００２８】
なお、画像合成部26は、操作部40を介して設定された画像表示条件等に基づいて制御部42によって制御されている。操作部40は、マウス、キーボード、トラックボール、タッチペン、ジョイスティック等の操作デバイスを備えており、画像表示条件等の設定を行うことができる。
【００２９】
画像表示部28は、画像合成部26の画像選択部により選択された断層画像及び弾性画像、時間ばらつき画像等の画像を表示する。
【００３０】
ここで、図2を用いて時間ばらつき解析部50を説明する。時間ばらつき解析部50は、複数の弾性情報のフレームを記憶する弾性情報フレーム記憶部60と、複数の弾性情報のフレームから時間ばらつきを算出する時間ばらつき算出部62を有している。
【００３１】
弾性情報演算部34は、被検体10の断層部位の超音波断層データに基づいて、断層部位の各計測点における組織の歪み又は弾性率の弾性情報をフレーム毎に求めている。弾性情報フレーム記憶部60は、弾性情報演算部34から出力される複数の弾性情報のフレーム、例えば10フレームの弾性情報を記憶する。弾性情報フレーム記憶部60に記憶するフレーム数は、操作部40を介して任意に設定することができる。
【００３２】
超音波探触子12の圧迫周期に合わせて弾性情報フレーム記憶部60に記憶するフレーム数を設定することもできる。超音波探触子12の圧迫周期に合わせるのであれば、弾性情報フレーム記憶部60は圧迫1周期分の弾性情報のフレームを記憶する。
【００３３】
また、超音波の放射圧を用いて組織を変形させ、変位から弾性情報を計測する手法がある。超音波の放射圧の照射時間に合わせて弾性情報フレーム記憶部60に記憶するフレーム数を設定することもできる。超音波の放射圧の照射時間に合わせるのであれば、弾性情報フレーム記憶部60は超音波の放射圧の照射中の弾性情報のフレームを記憶する。
【００３４】
図3を用いて、時間ばらつき算出部62の算出手法について説明する。図3の左側に示すように、画像表示部28は、断層画像と弾性画像の合成画像を表示する。なお、画像表示部28は断層画像のみ表示していてもよい。
【００３５】
断層画像と弾性画像の合成画像では、断層画像における形態情報である正常組織70と腫瘍組織72に加えて、弾性画像における硬さ情報である硬化組織74が、断層画像における腫瘍組織72に重ね合わせて表示されている。腫瘍組織72が黒色で表示され、硬化組織74が青色で表示されるとすると、腫瘍組織72と硬化組織74が重なる部位に関しては、深い青色で表示されることになる。操作者は、深い青色で表示される部分が悪性腫瘍である可能性があるとして、観察することができる。
【００３６】
操作者は、操作部40を用いて、腫瘍組織72を含む領域に関心領域80を設定する。断層画像と弾性画像の合成画像における各計測点(X1〜Xn,Y1〜Yn)には計測点82がそれぞれ設定される。nは自然数である。関心領域80は操作者が注目する領域であり、弾性情報の比の比較対象となる領域である。
【００３７】
操作部40で設定された関心領域80を、弾性情報フレーム記憶部60に記憶された複数の弾性情報のフレームに対応させることができる。すなわち、操作部40で設定された関心領域80の座標を、弾性情報フレーム記憶部60に記憶された複数の弾性情報のフレームの座標に対応させることができる。よって、複数の弾性情報のフレームの同じ位置に関心領域が設定されることになる。
【００３８】
時間ばらつき算出部62は、断層画像と弾性画像の合成画像に設定された関心領域80を弾性情報フレーム記憶部60に記憶された複数の弾性情報のフレームに対応させる。そして、時間ばらつき算出部62は、複数の弾性情報のフレームにおいて、断層画像と弾性画像の合成画像に設定された関心領域80と各計測点82における弾性情報の指標値をそれぞれ演算する。ここでは、弾性情報の指標値として弾性情報の比を例にあげて説明する。時間ばらつき算出部62は、弾性情報の比に基づいて各計測点において時間ばらつきを算出する。
【００３９】
画像表示部28は、図3の下側に示すように、各計測点において弾性情報の比の時間変化をグラフで表示することもできる。ここでは、弾性情報である歪みを例にあげて説明する。関心領域80の歪みはε1であり、計測点82の歪みはε2である。関心領域80の歪みε1と計測点82の歪みε2は、それぞれの領域における歪みの平均値である。時間ばらつき算出部62は、歪みの比ε2/ε1(εratio)を演算する。画像表示部28には、時間ばらつき算出部62で演算された歪みの比ε2/ε1が時間とともに表示される。なお、弾性情報の歪みを例にあげて説明したが、変位、弾性率、粘性などのパラメータであってもよい。
【００４０】
ここで、時間ばらつき算出部62は、複数の弾性情報のフレームにおいて、各計測点における弾性情報の指標値(弾性情報の比)のばらつきを算出する。時間ばらつき算出部62は、複数の弾性情報のフレームにおいて、各計測点における弾性情報の比の上限値と下限値の差D1を算出する。時間ばらつき算出部62は、弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より大きい場合、時間ばらつきが大きいと算出する。時間ばらつき算出部62は、弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より小さい場合、時間ばらつきが小さいと算出する。
【００４１】
具体的には、計測点(X1,Y1)の場合、弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より小さく、弾性情報の比のばらつきが小さい。時間ばらつき算出部62は、計測点(X1,Y1)を時間ばらつきが小さい計測点であると算出する。計測点(X2,Y2)の場合、弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より大きく、弾性情報の比のばらつきが大きい。時間ばらつき算出部62は、計測点(X2,Y2)を時間ばらつきが大きい計測点であると算出する。このように、時間ばらつき算出部62は、断層画像と弾性画像の合成画像の各計測点(X1〜Xn,Y1〜Yn)における弾性情報の比の上限値と下限値の差と予め定められた閾値Th1とを比較して、各計測点における弾性情報の比のばらつきを算出する。
【００４２】
時間ばらつき算出部62は、各計測点における弾性情報の比の上限値と下限値の差を用いて時間ばらつきを算出したが、各計測点において時間方向に弾性情報の標準偏差を取った下記式によって求めることもできる。
【００４３】
【数１】
2014061370.tif 000003
また、時間ばらつき算出部62は、上記の他に、分散、変動係数などの代用が可能である。時間ばらつき算出部62は、弾性情報の分散の大きさに基づいて時間ばらつきを算出したり、弾性情報の変動係数の大きさに基づいて時間ばらつきを算出する。
【００４４】
図3の右側に示すように、時間ばらつき画像構成部52は、時間ばらつき解析部50で解析された時間ばらつきが閾値より大きい計測点について、時間ばらつきが閾値より大きいことを示す時間ばらつき画像86を構成する。具体的には、時間ばらつき画像構成部52は、時間ばらつきが閾値より大きい計測点について識別される時間ばらつき表示バー84に基づいて、時間ばらつき画像86を構成する。
【００４５】
時間ばらつき表示バー84は、時間ばらつきが大きい計測点について時間ばらつき画像が表示されるように設定されている。例えば、時間ばらつき表示バー84は、弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より大きい場合、時間ばらつき画像が表示されるように設定され、弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より小さい場合、時間ばらつき画像が表示されないように設定されている。
【００４６】
すなわち、時間ばらつき画像構成部52は、弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より大きい場合、時間ばらつき画像を構成して画像表示部28に表示させる。弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より小さい場合、時間ばらつき画像を構成しない。
【００４７】
図3の右側に示すように、画像表示部28は、断層画像と時間ばらつき画像の合成画像を表示する。操作者は、時間ばらつき画像が表示された領域は時間ばらつきがある領域であると認識することができる。
【００４８】
図4を用いて、時間ばらつき算出部62の他の算出手法について説明する。図4の左側に示すように、画像表示部28は、断層画像と弾性画像の合成画像を表示する。
【００４９】
操作者は、操作部40を用いて、関心領域90を設定する。関心領域90を設定すると、断層画像と弾性画像の合成画像では、断層画像における形態情報である正常組織70と腫瘍組織72に加えて、弾性画像における硬さ情報である硬化組織74が、断層画像における腫瘍組織72に重ね合わせて表示されている。
【００５０】
断層画像と弾性画像の合成画像の関心領域90における各計測点(Xa〜Xm,Ya〜Ym)には計測点82がそれぞれ設定される。a、mは自然数であり、a＜mである。
【００５１】
時間ばらつき算出部62は、弾性情報フレーム記憶部60に記憶された複数の弾性情報のフレームにおいて、設定された関心領域90と計測点82における弾性情報を算出する。時間ばらつき算出部62は、弾性情報の変位に基づいて各計測点において時間ばらつきを算出する。
【００５２】
画像表示部28は、図4の下側に示すように、各計測点において弾性情報の変位をグラフで表示することもできる。ここでは、弾性情報である歪みを例にあげて説明する。計測点82の歪みはε3である。画像表示部28には、時間ばらつき算出部62で演算された歪みε3が時間とともに表示される。なお、弾性情報の歪みを例にあげて説明したが、変位、弾性率、粘性などのパラメータであってもよい。
【００５３】
ここで、時間ばらつき算出部62は、複数の弾性情報のフレームにおいて、各計測点における弾性情報のばらつきを算出する。時間ばらつき算出部62は、複数の弾性情報のフレームにおいて、各計測点における弾性情報の変位D2を算出する。時間ばらつき算出部62は、弾性情報の変位D2が予め定められた閾値Th2より大きい場合、時間ばらつきが大きいと算出する。時間ばらつき算出部62は、弾性情報の変位D2が予め定められた閾値Th2より小さい場合、時間ばらつきが小さいと算出する。
【００５４】
具体的には、計測点(X3,Y3)の場合、弾性情報の変位D2が予め定められた閾値Th2より小さく、弾性情報のばらつきが小さい。時間ばらつき算出部62は、計測点(X3,Y3)を時間ばらつきが小さい計測点であると算出する。計測点(X4,Y4)の場合、弾性情報の変位D2が予め定められた閾値Th2より大きく、弾性情報のばらつきが大きい。時間ばらつき算出部62は、計測点(X4,Y4)を時間ばらつきが大きい計測点であると算出する。このように、時間ばらつき算出部62は、弾性画像の表示領域である関心領域90の各計測点(Xa〜Xm,Ya〜Ym)における弾性情報の変位と予め定められた閾値Th2とを比較して、各計測点における弾性情報のばらつきを算出する。
【００５５】
図4の右側に示すように、時間ばらつき画像構成部52は、時間ばらつき解析部50で解析された時間ばらつきが大きい計測点について、時間ばらつきが大きいことを示す時間ばらつき画像96を構成する。時間ばらつき画像96を表示する表示領域92は、弾性画像の表示領域として表示された関心領域90と同じ領域である。
【００５６】
具体的には、時間ばらつき画像構成部52は、時間ばらつきが大きい計測点について識別される時間ばらつき表示バー94に基づいて、時間ばらつき画像96を構成する。時間ばらつき表示バー94は、時間ばらつきが大きい計測点について時間ばらつき画像が表示されるように設定されている。
【００５７】
図4の右側に示すように、画像表示部28は、断層画像と時間ばらつき画像の合成画像を表示する。操作者は、時間ばらつき画像が表示された領域は時間ばらつきがある領域であると認識することができる。
【００５８】
よって、操作者は、時間ばらつき画像が表示された領域は時間ばらつきがある領域であると認識することができるため、弾性情報フレーム記憶部60に記憶された複数の弾性情報のフレームにおいて、弾性情報演算部34は時間ばらつき画像が表示された領域以外で弾性情報の指標値の演算を行うことができる。よって、操作者は、時間ばらつき画像が表示された領域以外に弾性情報の指標値を演算するための関心領域を設定すれば、時間ばらつきの影響を受けることなく、弾性情報の指標値の演算を行うことができる。つまり、安定した計測位置で弾性情報の指標値の演算を行うことができる。
【００５９】
なお、時間ばらつき画像が所定領域以上(例えば、画面の半分以上)であり、弾性情報演算部34によって弾性情報の指標値を演算することが適切でない場合、制御部42は、撮影し直すアラームを画像表示部28に表示させる。撮影し直し、時間ばらつき画像が所定領域より小さい(例えば、画面の半分より小さい)場合、制御部42は弾性画像が安定したとみなし、弾性情報演算部34によって弾性情報の指標値を演算する。時間ばらつき画像構成部52は、撮影し直した際の時間ばらつき画像を構成し、画像表示部28に更新した時間ばらつき画像を表示させることもできる。
【００６０】
次に本実施形態の動作について図5を用いて説明する。
(S100)画像表示部28は、断層画像構成部24によって構成された断層画像と、弾性画像構成部36によって構成されたた歪み又は弾性率に基づく弾性画像の合成画像を表示する。
【００６１】
(S102)断層画像と弾性画像の合成画像上の所定の計測位置で弾性情報の指標値の演算を行うかどうかを選択する。弾性情報の指標値の演算を行なう場合、S104に進む。弾性情報の指標値の演算を行わない場合、終了する。
【００６２】
(S104)操作部40を用いて、関心領域80を設定する。弾性画像の表示領域として関心領域90を設定してもよい。
【００６３】
(S106)時間ばらつき解析部50は、複数の弾性情報のフレームにおいて、設定された関心領域と計測点における弾性情報の指標値若しくは変位を算出する。そして、弾性情報の指標値若しくは変位と予め定められた閾値とを比較して、各計測点における弾性情報のばらつきを算出する。
【００６４】
(S108)時間ばらつき画像構成部52は、時間ばらつき解析部50で解析された時間ばらつきが大きい計測点について、時間ばらつきが大きいことを示す時間ばらつき画像86、96を構成する。
【００６５】
(S110)操作者は、時間ばらつき画像が表示された領域以外に弾性情報の指標値を演算するための関心領域を設定し、弾性情報演算部34は弾性情報の指標値の演算を行なう。
【００６６】
なお、時間ばらつき画像が表示された領域以外に弾性情報の指標値を演算するための関心領域は、自動で設定することもできる。具体的には、制御部42において、時間ばらつき画像が表示された領域を認識し、時間ばらつき画像が表示された領域以外に関心領域を設定する。制御部42は、時間ばらつき画像が表示された領域以外の弾性画像における硬さ情報である硬化組織74に第1の関心領域を設定し、時間ばらつき画像が表示された領域以外であり硬化組織74以外に第2の関心領域を設定する。弾性情報演算部34は、弾性情報フレーム記憶部60に記憶された複数の弾性情報のフレームにおいて、制御部42によって設定された第1の関心領域と第2の関心領域における弾性情報の指標値(例えば、弾性情報の比)を算出する。弾性情報演算部34によって演算された弾性情報の指標値は、画像合成部26を介して、画像表示部28に表示される。
【００６７】
また、操作部40を介して第1の関心領域と第2の関心領域を予め設定しておき、時間ばらつき画像構成部52によって構成された時間ばらつき画像が第1の関心領域と第2の関心領域に重なる場合、弾性情報演算部34によって第1の関心領域と第2の関心領域における弾性情報の指標値を演算しない。制御部42は、撮影し直すアラームを画像表示部28に表示させる。時間ばらつき画像構成部52によって構成された時間ばらつき画像が第1の関心領域と第2の関心領域に重ならない場合、弾性情報演算部34によって第1の関心領域と第2の関心領域における弾性情報の指標値を演算する。弾性情報演算部34によって演算された弾性情報の指標値は、画像合成部26を介して、画像表示部28に表示される。
【００６８】
以上、本実施形態によれば、超音波探触子12を介して被検体の診断部位の断層画像を構成する断層画像構成部24と、硬さを示す弾性情報を演算する弾性情報演算部34と、弾性情報演算部34において演算された弾性情報に基づいて弾性画像を構成する弾性画像構成部36と、断層画像と弾性画像を表示する画像表示部28とを備えた超音波診断装置において、弾性情報演算部34で演算された弾性情報から各計測点における時間ばらつきを解析する時間ばらつき解析部50と、時間ばらつき解析部で解析された時間ばらつきに基づいて、時間ばらつき画像を構成する時間ばらつき画像構成部52とを備え、画像表示部28は時間ばらつき画像を表示する。
【００６９】
よって、弾性情報の時間ばらつきに基づく時間ばらつき画像を表示することができる。そして、時間ばらつき画像が表示された領域以外の安定した計測位置で弾性情報の指標値の演算を行うことができる。つまり、操作者に依存せず、安定した計測位置で弾性情報の指標値の演算を行うことができる。
【００７０】
(第2の実施形態)
第2の実施形態について、図6、7を用いて説明する。第1の実施形態と異なる点は、図6に示すように、弾性情報演算部34における弾性情報の指標値の演算を行う演算時間を設定する演算時間設定部100を備える点である。
【００７１】
図7に示すように、時間ばらつき解析部50の時間ばらつき算出部62は、第1の実施形態で示すように、複数の弾性情報のフレームにおいて、各計測点における弾性情報の指標値(弾性情報の比)のばらつきを算出する。
【００７２】
図3に示すように、計測点(X2,Y2)の場合、弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より大きく、弾性情報の比のばらつきが大きい。時間ばらつき算出部62は、計測点(X2,Y2)を時間ばらつきが大きい計測点であると算出したが、時間t1〜時間t2の時間帯は安定している。安定している時間t1〜時間t2の時間帯とは、弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より小さくなる時間帯である。
【００７３】
図7に示すように、時間t1〜時間t2の時間帯110は、弾性情報の比の上限値と下限値の差D1が予め定められた閾値Th1より小さくなる時間帯がある。つまり、時間t1〜時間t2の時間帯110を安定している時間帯とみなすことができる。そこで、演算時間設定部100は、安定している時間t1〜時間t2の時間帯を選択し、安定している時間t1〜時間t2の時間帯110を演算時間設定部100の演算時間として設定する。画像表示部28は、安定している時間t1〜時間t2の時間帯110を表示する。
【００７４】
弾性情報演算部34は、演算時間設定部100によって設定された時間t1〜時間t2の時間帯110おける弾性情報の指標値を演算する。時間t1〜時間t2の時間帯から複数フレームを選択し、選択されたフレーム112における弾性情報の指標値を演算することもできる。画像表示部28は、選択されたフレーム112を表示する。
【００７５】
弾性情報演算部34は、時間t1〜時間t2の時間帯で選択されるフレーム数を予め決めておけば、フレーム数に応じた弾性情報の指標値を演算することができる。例えば、フレーム数を10と設定しておけば、時間t1〜時間t2の時間帯で10フレームを選択し、弾性情報演算部34は、選択された10フレームから弾性情報の指標値を演算することができる。
【００７６】
本実施形態によれば、時間ばらつき画像が表示された領域以外の安定した計測位置であり、かつ安定した時間帯で弾性情報の指標値の演算を行うことができる。つまり、操作者に依存せず、安定した計測位置で弾性情報の指標値の演算を行うことができる。
【符号の説明】
【００７７】
10  被検体、12  超音波探触子、14  送信部、18  超音波送受信制御部、20  受信部、22  整相加算部、24  断層画像構成部、26  画像合成部、28  画像表示部、30  RF信号フレームデータ選択部、32  変位計測部、34  弾性情報演算部、36  弾性画像構成部、40  操作部、42  制御部、50  時間ばらつき解析部、52  時間ばらつき画像構成部、60  弾性情報フレーム記憶部、62  時間ばらつき算出部
【図１】
2014061370.tif 000004
【図２】
2014061370.tif 000005
【図３】
2014061370.tif 000006
【図４】
2014061370.tif 000007
【図５】
2014061370.tif 000008
【図６】
2014061370.tif 000009
【図７】
2014061370.tif 000010
【国際調査報告】
2014061370.tif 000011
2014061370.tif 000012
2014061370.tif 000013
2014061370.tif 000014