(57)【特許請求の範囲】
【請求項１】
超音波探触子と、前記超音波探触子を介して被検体の断層部位の超音波断層データを計測する超音波断層データ計測手段と、計測された前記超音波断層データを処理して断層像を生成する信号処理手段と、計測された前記超音波断層データに基づいて前記断層部位における組織の弾性情報を求め、該弾性情報に基づいて前記断層部位における弾性像を生成する弾性像生成手段と、前記断層像及び／又は前記弾性像を表示する表示装置と、該表示装置に表示された前記弾性像に複数の関心領域を設定する入力手段と、設定された第1の関心領域と第2の関心領域における前記弾性情報の比を算出する指標化処理手段とを備えることを特徴とする超音波診断装置。
【請求項２】
前記指標化処理手段は、前記第１の関心領域及び前記第２の関心領域における前記弾性情報の平均値に基づいて、前記弾性情報の比を算出することを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
前記表示装置は、前記弾性情報の比を前記断層像及び／又は前記弾性像とともに表示することを特徴とする請求項１に記載の超音波診断装置。
【請求項４】
前記弾性像生成手段は、時系列的に計測される複数の前記超音波断層データに基づいて複数の前記弾性像を生成し、前記表示装置に表示された複数の前記弾性像の計測時相を指定する時相バーが入力手段によって操作されたとき、指定された計測時相に対応する前記弾性像を前記表示装置に表示させることを特徴とする請求項１に記載の超音波診断装置。
【請求項５】
前記指標化処理手段は、前記表示装置に前記弾性情報のスケールを表示させ、該スケールの前記第１の関心領域と前記第２の関心領域の前記弾性情報又は前記弾性情報の平均値に対応する位置に表示マークを表示させることを特徴とする請求項１又は２に記載の超音波診断装置。
【請求項６】
前記指標化処理手段は、前記各関心領域の前記弾性情報の平均値を求め、該平均値の最小値と最大値をカラー化の色相情報のダイナミックレンジの下限値と上限値に対応付けることを特徴とする請求項１に記載の超音波診断装置。
【請求項７】
前記指標化処理手段は、前記関心領域が３つ以上設定され、前記弾性情報の平均値の最小値と最大値との間に中間値があるとき、該中間値に対応する前記ダイナミックレンジのレベルを可変設定することを特徴とする請求項１に記載の超音波診断装置。
【請求項８】
前記表示装置は、前記弾性情報の比の時間変化を表示することを特徴とする請求項１に記載の超音波診断装置。
【請求項９】
前記指標化処理手段は、前記各関心領域における前記弾性情報を複数フレームにわたって順次加算した加算値の変化を前記表示装置に表示させる請求項１に記載の超音波診断装置。
【請求項１０】
前記第２の関心領域は、前記第１の関心領域を包囲する広い領域であることを特徴とする請求項１に記載の超音波診断装置。
【請求項１１】
前記第２の関心領域は、前記第１の関心領域から離して設定されていることを特徴とする請求項１に記載の超音波診断装置。
【請求項１２】
前記複数の関心領域は、病変部に設定された第２の関心領域と、該第２の関心領域を挟んで前記弾性像の横方向に離して設定された第１及び第３の関心領域とを含み、
前記指標化処理手段は、前記第２の関心領域の前記弾性情報と前記第１及び第３の関心領域の前記弾性情報の和とから前記弾性情報の比を算出することを特徴とする請求項１に記載の超音波診断装置。
【請求項１３】
前記第１の関心領域は、前記被検体の正常部位に設定されることを特徴とする請求項１に記載の超音波診断装置。
【請求項１４】
前記第１の関心領域は、前記被検体の脂肪層に設定されることを特徴とする請求項１に記載の超音波診断装置。
【請求項１５】
前記第１の関心領域は、前記被検体の体表から一定深度の領域に設定されることを特徴とする請求項１に記載の超音波診断装置。
【請求項１６】
前記弾性情報の比に基づいて前記関心領域の組織が正常か否かを判定する判定手段を備えたことを特徴とする請求項１に記載の超音波診断装置。
【請求項１７】
前記指標化処理手段は、２組の前記第１と前記第２の関心領域から算出される２つの前記弾性情報の比に基づいて、前記弾性情報の比の比を算出することを特徴とする請求項１に記載の超音波診断装置。
【請求項１８】
前記超音波探触子を支持すると共に該超音波探触子によって前記被検体に加える圧迫を加減する加圧制御手段と、該加圧制御手段を支持する固定手段とを備えたことを特徴とする請求項１に記載の超音波診断装置。
【請求項１９】
前記弾性情報は、前記歪み又は前記弾性率を含む物理量であることを特徴とする請求項１乃至１８のいずれか１項に記載の超音波診断装置。
【請求項２０】
超音波探触子と、前記超音波探触子を介して被検体の断層部位の超音波断層データを計測する超音波断層データ計測手段と、計測された前記超音波断層データを処理して断層像を生成する信号処理手段と、計測された前記超音波断層データに基づいて前記断層部位における組織の歪み又は弾性率を求め、該歪み又は弾性率に基づいて前記断層部位における弾性像を生成する弾性像生成手段と、前記断層像及び／又は前記弾性像を表示する表示装置と、該表示装置に表示された前記弾性像に複数の関心領域を設定する入力手段と、設定された第1の関心領域と第2の関心領域における前記歪み又は弾性率の比を算出する指標化処理手段とを備えることを特徴とする超音波診断装置。
【請求項２１】
被検体に圧迫を加えながら取得された超音波断層データに基づいて前記断層部位における組織の弾性情報を求めるステップと、該弾性情報に基づいて前記断層部位における弾性像を生成して表示するステップと、表示された前記弾性像に複数の関心領域を設定するステップと、設定された第1の関心領域と第2の関心領域における前記弾性情報の比を算出するステップとを含む弾性像表示方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断における診断部位の生体組織の硬さ（以下、弾性という。）を定量的に評価するのに好適な弾性像表示方法及び超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波探触子により被検体内部に超音波を送信し、被検体内部から生体組織の構造に応じた超音波の反射エコー信号を受信し、例えばＢモード像等の断層像を再構成して診断用に表示する。【０００３】
最近は、用手法又は機械的な方法により被検体に圧迫を加えて超音波画像データを計測し、計測時間が異なる２つの超音波画像のフレームデータに基づいて圧迫により生じた生体各部の変位を求め、その変位データに基づいて生体組織の弾性を表す弾性像を生成することが提案されている。生体組織の弾性に相関する物理量には、種々のものが知られているが、例えば、生体組織の歪み又は弾性率が用いられている。ここで、歪みは生体組織の移動量である変位を空間微分して得られる相対的な値であり、弾性率は生体組織の各部位に作用する応力変化を変位で除した定量的な値である。【０００４】
また、弾性像としては、例えば、特許文献１に示すように、歪み又は弾性率に応じて組織の各部位に赤、青、その他の色相情報を付与してカラー化した弾性像が知られている。【０００５】
ところで、生体組織の変位は圧迫力の大きさによって変わるから、同一部位の生体組織であっても、強く圧迫すると歪みが大きくなる。そのため、歪みを表すカラー弾性像は、そのカラー弾性像上の各部位の歪みの相対表示を表すに過ぎず、硬さを定量的に評価することができない。【０００６】
そこで、特許文献２には、Ｂモード像に２つの関心領域（ＲＯＩ）を設定し、その２つのＲＯＩの歪みを計測し、それぞれのＲＯＩに対応付けて歪みの数値を画像に表示することが提案されている。これによれば、２つのＲＯＩの歪みを比較することにより、関心部位の硬さを相対的に評価することができる。【０００７】
すなわち、特許文献２に記載の技術によれば、２つのＲＯＩを癌等の腫瘍と思われる部位と腫瘍が存在しない部位とに設定して、それら２つのＲＯＩの歪みを対比して、腫瘍が存在しない部位に対する腫瘍等の歪みを相対的に評価することができる。【０００８】
しかし、Ｂモード像は組織の歪みを表示するものではないので、歪みが異なる領域を認識し難いから、特許文献２に記載の技術によれば、歪みを対比しようとする２つのＲＯＩの大きさ及び位置を適切に設定するのが難しい。【０００９】
また、特許文献２では、例えば、正常組織の部位又は硬さに個体差が少ない生体組織の部位を基準として、腫瘍などの病変部と思われる部位の硬さを比などの指標により表すことについては考慮されていないから、硬さを定量的に評価をすることができない。【００１０】
特許文献１：ＪＰ２０００−６０８５３Ａ
特許文献２：ＵＳ２００２／０１７８８３３Ａ１
発明の開示
【発明の開示】
【００１１】
本発明は、診断部位の生体組織の硬さを定量的に評価することができる弾性像表示方法及び超音波診断装置を実現することを課題とする。
【００１２】
上記課題を解決する本発明の弾性像表示方法及びその方法を適用した装置は、超音波探触子と、前記超音波探触子を介して被検体の断層部位の超音波断層データを計測する超音波断層データ計測手段と、計測された前記超音波断層データを処理して断層像を生成する信号処理手段と、計測された前記超音波断層データに基づいて前記断層部位における組織の弾性情報を求め、該弾性情報に基づいて前記断層部位における弾性像を生成する弾性像生成手段と、前記断層像及び／又は前記弾性像を表示する表示装置と、該表示装置に表示された前記弾性像に複数の関心領域を設定する入力手段と、設定された第1の関心領域と第2の関心領域における前記弾性情報の比を算出する指標化処理手段とを備え、該弾性情報の比を前記表示装置に表示することを特徴とする。
【００１３】
本発明によれば、組織の弾性に相関する物理量、すなわち歪み又は弾性率などの弾性情報を対比しようとする複数の関心領域（ＲＯＩ）を弾性像上で設定するようにしているから、例えば、歪みが異なる領域を容易に認識でき、対比しようとする複数のＲＯＩの大きさ及び位置を適切に設定することができる。特に、対比に係る複数の関心領域における弾性情報を相互に指標値（例えば、相互の比）化し、その指標値を表示装置に表示するようにしているから、圧迫力に関係なく、関心部位の硬さの違いを定量的に評価をすることができ、精度の高い診断に結びつけることができる。例えば、正常組織の部位又は硬さに個体差が少ない生体組織の部位に設定したＲＯＩを基準として、腫瘍などの病変部と思われる部位に設定したＲＯＩの硬さを指標により表すことができる。
【００１４】
また、関心領域に関し、第１の関心領域と、該第１の関心領域を包囲する広い領域に設定された第２の関心領域とが設定される。この第２の関心領域は、第１の関心領域を含まない領域である。よって、悪性領域（第１の関心領域）と良性領域（第２の関心領域）とを区別して関心領域を設定できるため、良性領域に対する悪性領域の比較、及びその診断を正確に行うことができる。しかし、本発明はこれに限らず、第1の関心領域と第２の関心領域を離して設定することができる。いずれの場合も、指標値として、第１の関心領域と第２の関心領域の物理量の比を用いることができ、また、各関心領域の物理量の平均値の比を用いることができる。
【００１５】
さらに、表示装置に物理量のスケールを表示し、そのスケールの第１の関心領域と第２の関心領域の物理量が対応する位置に表示マークを表示することができる。これによれば、一目で、第１と第２の関心領域の組織の硬さの違いを定量的に認識することができる。
【００１６】
また、関心領域は、病変部に設定された第1の関心領域と、この第１の関心領域を挟んで弾性像の横方向に離して第２及び第３の関心領域を設定することができる。この場合、指標値として、第１の関心領域の物理量と、第２及び第３の関心領域の前記物理量の和との比を用いることができる。これによれば、各関心領域に対する加圧が一方に偏っても、加圧の弱い箇所と、加圧の強い箇所の平均を取ることにより、応力を均衡させることができる。
【００１７】
さらに、弾性像をカラー化して表示する場合は、各関心領域の物理量の平均値を求め、その平均値の最小値と最大値をカラー化の色相情報のダイナミックレンジの下限値と上限値に対応付けることが好ましい。これにより、色表示の解像度を向上できる。この場合、関心領域が３つ以上設定され、物理量の平均値の最小値と最大値との間に中間値があるときは、その中間値に対応するダイナミックレンジのレベルを可変設定することができる。これにより、或る関心領域間の階調を変化させ、多くの階調で表示させることにより、視覚的に観察しやすくすることができる。
【００１８】
また、本発明の弾性像表示方法は、同一断層部位における複数の関心領域の組織の硬さを対比することに限られるものではない。例えば、左右の手又は足、乳房などのように、被検体の対称性を有する左右の断層部位について別々に前記超音波断層データを計測し、前記左右の断層部位における弾性像を生成して表示装置に表示し、前記表示装置に表示された前記左右の前記弾性像にそれぞれ２つの関心領域を設定し、前記左右の２つの関心領域の前記物理量の比を求め、さらに前記左右の関心領域の前記物理量の比の比を求めて前記表示装置に表示することができる。これにより、左右の手又は足、乳房などの対応する部位の組織の硬さの違いを定量的に評価して、診断の精度を高めることができる。
【００１９】
さらに、本発明の弾性像表示方法は、被検体に圧迫を加えながら前記被検体の断層部位の超音波断層データを複数フレームにわたって計測し、該複数フレームの超音波断層データに基づいて前記断層部位における組織の弾性に相関する物理量の時間変化を求め、該物理量の時間変化に基づいて前記断層部位における弾性像の動画像を生成して表示装置に表示し、前記表示装置に表示された前記弾性像の静止画像上に少なくとも２つの関心領域を設定し、前記各関心領域における前記物理量の時間変化を前記表示装置に表示するようにすることができる。
【００２０】
これによれば、圧迫力の変化に応じて変化する歪みの各関心領域の定量的な対比を的確に行うことができる。この場合において、各関心領域における物理量を指標値化し、その指標値の時間変化を表示装置に表示するようにすることができる。また、各関心領域における指標値を複数フレームにわたって順次加算した加算値の変化を表示装置に表示するようにすることもできる。さらに、複数フレームにわたる各関心領域における指標値の時間平均を求めて、表示装置に表示するようにすることもできる。
【００２１】
上述した本発明の弾性像表示方法を実施する本発明の超音波診断装置は、被検体に超音波探触子によって圧迫を加えながら前記被検体の断層部位の超音波断層データを計測する超音波断層データ計測手段と、計測された前記超音波断層データを処理して断層像を生成する信号処理手段と、計測された前記超音波断層データに基づいて前記断層部位における組織の弾性に相関する物理量を求め、該物理量に基づいて前記断層部位における弾性像を生成する弾性情報演算手段と、前記断層像及び／又は前記弾性像を表示する表示装置と、該表示装置に表示された前記弾性像の静止画像に複数の関心領域を設定する入力手段と、設定された前記各関心領域における前記物理量を指標値化して前記表示装置に表示させる指標化処理手段とを含んで構成することができる。
【００２２】
また、超音波探触子を支持すると共に該超音波探触子によって被検体に加える圧迫を加減する加圧制御手段と、この加圧制御手段を支持する固定手段とを備えて構成することができる。これによれば、被検体の反発力によらず、同じ圧力、同じ速度で生体組織に圧迫をかけることができるので、再現性及び評価性に優れた歪み弾性像を得ることができ、本発明に係る弾性に相関する物理量の指標値を安定化することができる。
【図面の簡単な説明】
【００２３】
【図１】本発明の弾性像表示方法を適用してなる超音波診断装置の一実施形態の概略構成を示す図である。
【図２】本発明の弾性像表示方法の一実施形態の処理手順を示すフローチャートである。
【図３】弾性像の表示例を示す図である。
【図４】弾性像の他の表示例を示す図である。
【図５】図４の表示画像の一例におけるカラースケールの表示例を示す図である。
【図６】ＲＯＩの設定法の他の例を示す図である。
【図７】ＲＯＩの設定法のさらに他の例を示す図である。
【図８Ａ】各ＲＯＩの歪みの時間変化の表示例を示す図である。
【図８Ｂ】各ＲＯＩ間の歪み比の時間変化の表示例を示す図である。
【図８Ｃ】各ＲＯＩ間の歪み比の時間変化を加算して示す表示例図である。
【図９Ａ】各ＲＯＩの歪みor弾性率の時間変化の表示の一例を示す図である。
【図９Ｂ】各ＲＯＩの歪みor弾性率の時間変化の表示の他の一例を示す図である。
【図１０】弾性像と各ＲＯＩに適合した歪み平均値の比のグラフを同一画面上に表示する形態を示す図である。
【図１１】図１０の表示形態の変形例を示す図である。
【図１２Ａ】一般的なダイナミックレンジの設定例を示す図である。
【図１２Ｂ】本発明に係るダイナミックレンジの設定例の一実施形態を示す図である。
【図１２Ｃ】本発明に係るダイナミックレンジの設定例の他の実施形態を示す図である。
【図１３】本発明の弾性像表示方法のさらに他の実施形態の特徴を説明する図である。
【図１４】本発明の超音波診断装置の自動加圧装置の一実施形態を示す図である。
【図１５】図１４の探触子を固定する加圧システムの詳細構成を示す図である。
【発明を実施するための最良の形態】
【００２４】
以下、本発明の実施の形態を添付図面に基づいて説明する。
（実施形態１）
図１に本発明の一実施形態の超音波診断装置の概略構成を示し、図２に本発明の一実施形態の弾性像表示方法の処理手順にフローチャートを示す。図１に示すように、本実施形態の超音波診断装置は、被検体１に当接させて用いる超音波の探触子２と、超音波送受信部３と、整相加算回路４と、断層像構成部５と、歪み演算部６と、弾性画像構成部７と、グラフィック部８と、カラースケール発生部９と、画像の合成部１０と、画像表示器１１と、制御演算部１２と、キーボード１３とを備えて構成されている。【００２５】
探触子２は、複数の振動子を配設して形成されており、被検体１の所定の断層部位に電子的にビーム走査を行い、被検体１との間で超音波を送受信する機能を有している。【００２６】
超音波送受信部３は、探触子２を介して被検体１の所定の断層部位に、時間間隔をおいて超音波を繰り返し送信すると共に、その断層部位からの反射エコー信号を受信するように構成されている。つまり、送信時には、探触子２を駆動して超音波を発生させるための送波パルスを生成するとともに、送信される超音波の収束点をある深さに設定する機能を有している。また、受信時には、探触子２で受信した反射エコー信号を所定のゲインで増幅してＲＦ信号すなわち受波信号を生成するものである。【００２７】
整相加算回路４は、超音波送受信部３により増幅されたＲＦ信号を入力して位相を合わせて加算し、複数の収束点に対し収束させた超音波断層データであるＲＦ信号データを時系列に生成するようになっている。【００２８】
これらの探触子２、超音波送受信部３及び整相加算回路４によって、被検体１に探触子２によって圧迫を加えながら、被検体の断層部位の超音波断層データを計測する超音波断層データ計測手段が構成されている。【００２９】
断層像構成部５は、超音波断層データを処理して断層像を生成する信号処理手段を構成するものであり、整相加算回路４から出力されるＲＦ信号データに基づいて、被検体１のの断層部位の濃淡断層像、例えば白黒断層像を再構成するようになっている。つまり、断層像構成部５は、信号処理部と白黒スキャンコンバータを含んで構成されている。信号処理部は、整相加算回路４からのＲＦ信号データを入力してゲイン補正、ログ圧縮、検波、輪郭強調、フィルタ処理等の信号処理を行って、断層像データを得るものである。また、白黒スキャンコンバータは、信号処理部からの断層像データをデジタル信号に変換するＡ／Ｄ変換器と、変換された複数の断層像データを時系列に記憶するフレームメモリと、制御コントローラなどを含んで構成されている。また、断層像構成部５は、白黒スキャンコンバータやフレームメモリに格納された被検体１内の断層フレームデータを１画像として取得し、取得された断像フレームデータをテレビ同期で読み出すものである。【００３０】
歪み演算部６は、整相加算回路４から出力されるＲＦ信号データに基づいて、断層部位の被検体１の断層部位の生体組織の変位を計測し、後述するように、歪みを求めるようになっている。弾性画像構成部７は、弾性データ処理部及びカラースキャンコンバータを含んで構成され、歪み演算部６で求められた歪みに基づいてカラー弾性画像を再構成するようになっている。つまり、歪み演算部６と弾性画像構成部７によって、超音波断層データに基づいて断層部位における組織の弾性に相関する物理量の一つである歪みを求め、その物理量に基づいて断層部位における弾性像を生成する弾性情報演算手段が構成されている。【００３１】
グラフィック部８は、超音波信号以外の画像を描出するようになっている。カラースケール発生部９は、後述するように、弾性像の歪みのカラースケールを発生するようになっている。合成部１０は、断層像構成部５から出力される白黒断層像と、弾性画像構成部７から出力されるカラー弾性像と、グラフィック部８から出力される弾性データなどのデータ画像と、カラースケール発生部９から出力されるカラースケールを、所定のアドレス上に合成するようになっている。画像表示器１１は、合成部１０で合成された合成画像を表示するようになっている。なお、合成部１０は、フレームメモリと、画像処理部と、画像選択部とを備えて構成され、入力される白黒断層像と、カラー弾性像と、データ画像と、カラースケールとを、制御演算部１２から与えられる指令に応じて選択して合成するようになっている。【００３２】
制御演算部１２は、予め設定された処理手順にしたがって、また、キーボード１３から適宜入力される指令に従って、上述したした各構成要素を制御するようになっている。キーボード１３は、各種設定用のインタフェースとなるものである。特に、本実施形態のキーボード１３は、画像表示器１１に表示された弾性像の静止画像に複数の関心領域を設定する入力手段を構成している。そして、制御演算部１２は、キーボード１３を介して弾性像に設定された各関心領域における歪みを相互に指標値化し、グラフィック部８を介して画像表示器１１に表示させる指標化処理手段を含んで構成されている。【００３３】
次に、本実施形態の超音波診断装置の詳細構成及び動作を、図２のフローチャートに示した処理手順にしたがって説明する。超音波診断装置の基本的な動作は、周知であることから、本発明の弾性像表示方法に関係する機能及び動作を中心に説明する。
（ステップＳ１）
歪み演算部６は、ＲＦ信号選択部及び変位演算部を含んで構成されており、整相加算回路４の後段に分岐して設けられている。ＲＦ信号選択部は、フレームメモリと選択部とを含んで構成されている。このＲＦ信号選択部は、整相加算回路４からの複数のＲＦ信号データをフレームメモリに格納し、格納されたＲＦ信号フレームデータ群から選択部により１組のＲＦ信号フレームデータ、すなわち計測時間が異なる２つのＲＦ信号フレームデータを選び出すものである。【００３４】
例えば、ＲＦ信号選択部は、整相加算回路４から画像のフレームレートに基づいて時系列で生成されるＲＦ信号データをフレームメモリ内に順次確保する。そして、図示していない制御部からの指令に応じて、現在確保されたＲＦ信号フレームデータ（Ｎ）を第１のデータとして選択部で選択する。これと同時に、時間的に過去に確保されたＲＦ信号フレームデータ群（Ｎ−１、Ｎ−２、Ｎ−３・・・Ｎ−Ｍ）の中から１つのＲＦ信号フレームデータ（Ｘ）を選択する。ここで、符号Ｎ、Ｍ、Ｘは、ＲＦ信号フレームデータに付されたインデックス番号であり、自然数である。
（ステップＳ２）
歪み演算部６の変位演算部は、１組のＲＦ信号フレームデータから生体組織の変位などを求めるものである。例えば、変位演算部は、ＲＦ信号選択部により選択された１組のＲＦ信号フレームデータ（Ｎ）及びＲＦ信号フレームデータ（Ｘ）に１次元或いは２次元の相関処理を行って、断層像の各点に対応する生体組織の変位を表す移動ベクトル（変位の方向と大きさ）に関する１次元又は２次元変位分布を求める。ここで、移動ベクトルの検出にはブロックマッチング法を用いる。
（ステップＳ３）
ステップＳ２で求めた生体組織の移動量である変位を空間微分することによって  歪みを算出する。つまり、変位演算部より計測された変位をΔＬとすると、歪み（Ｓ）は、ΔＬを空間微分することによって算出することができるから、Ｓ＝ΔＬ／ΔＸという式を用いて求められる。【００３５】
ここで、ブロックマッチング法とは、画像を例えばＮ×Ｎ画素からなるブロックに分け、関心領域内のブロックに着目し、着目しているブロックに最も近似しているブロックを前のフレームから探し、これを参照して予測符号化すなわち差分により標本値を決定する処理であり、周知の方法である。【００３６】
また、本実施形態では、弾性情報として歪みを求める場合を例に説明するが、本発明はこれに限られるものではなく、弾性情報として弾性率を求めることができる。弾性率は、圧力の変化を移動量の変化で除することによって計算されることから、生体組織の各部位に作用する圧力（応力）を求める必要がある。すなわち、例えば、変位演算部より計測された変位をΔＬとし、図示していない圧力計測部により計測された圧力をΔＰとすると、弾性率の一つであるヤング率Ｙｍは、Ｙｍ＝（ΔＰ）／（ΔＬ／Ｌ）という式によって算出できる。このヤング率Ｙｍから断層像の各点に対応する生体組織の弾性率が求められるので、２次元の弾性像データを連続的に得ることができる。なお、ヤング率とは、物体に加えらえた単純引張り応力と、引張りに平行に生じるひずみに対する比である。
（ステップＳ４）
弾性画像構成部７の弾性データ処理部は、歪み演算部６から時系列に出力される弾性フレームデータをフレームメモリに確保し、確保されたフレームデータを図示していない制御部の指令に応じて画像処理を行う。弾性画像構成部７のカラースキャンコンバータは、弾性データ処理部から出力される弾性フレームデータに基づいて色相情報を付した画像データに変換する。つまり、カラースキャンコンバータは、弾性フレームデータに基づいて光の３原色である赤（Ｒ）、緑（Ｇ）、青（Ｂ）を付した画像データに変換する。例えば、歪みが大きい弾性データを赤色コードに変換すると同時に、歪みが小さい弾性データを青色コードに変換する。なお、赤（Ｒ）緑（Ｇ）青（Ｂ）の階調は２５６であり、２５５は、大輝度で表示することを、逆に０は全く表示されないことを意味する。このカラースキャンコンバータには、制御演算部１２を介してキーボード１３などの操作部に接続されており、この操作部によって弾性画像の色合いなどが制御される。また、探触子２に圧力計（図示しない）を付け、探触子２を被検体１に押しつける圧力を計測する圧力計測部（図示しない）が設けることができる。【００３７】
合成部１０は、断層像構成部５、弾性画像構成部７及びグラフィック部８から出力される画像データをフレームメモリに格納する。また、合成部１０の画像処理部は、フレームメモリに確保された断層像データと弾性画像データを制御部の指令に応じた設定割合で加算して合成する。したがって、合成画像の各画素の輝度情報及び色相情報は、白黒断層像とカラー弾性像の各情報を設定割合で加算したものとなる。さらに、合成部１０の画像選択部は、フレームメモリ内の断層像データと弾性画像データ及び画像処理部の合成画像データのうちから画像表示器１１に表示する画像を制御部の指令に応じて選択して表示する。
（ステップＳ５）
キーボード１３から指令を入力して、画像表示器１１に弾性像を表示させると共に、表示された弾性像を任意のタイミングでフリーズさせ、フリーズで静止された弾性像上に、例えば図３に示すように、弾性を対比したい２つの関心領域ＲＯＩ１、ＲＯＩ２を設定する。例えば、ＲＯＩ１は、ＲＯＩ２を包囲する広い領域に設定する。すなわち、弾性像中で生体組織の歪みが小さく（硬い）、腫瘍などの病変と思われる部位にＲＯＩ２を設定する。そして、ＲＯＩ２の歪みと対比するために、弾性像上で正常部位と思われる領域に指標の基準となるＲＯＩ１を、ＲＯＩ２を包囲する広い領域に設定する。ここで、ＲＯＩ２は、ＲＯＩ１内を含まない領域である。なお、制御演算部１２は、キーボード１３から入力設定されたＲＯＩ１、ＲＯＩ２を、図３のように、白い点線などにより表示する指令をグラフィック部８に出力する。また、ＲＯＩを設定する画像は、白黒のＢモード像にカラーの弾性像を重ねた合成画像とすることができる。なお、図３は白黒で示してあるが、実際には赤や青などの色付け表示されており、容易に腫瘍等の病変部の広がりや大きさを診断することができるようにしている。
（ステップＳ６）
制御演算部１２は、図３の弾性像に設定されたＲＯＩ１とＲＯＩ２内のそれぞれ歪みを算出し、さらにその歪み平均値ＲＯＩ１＿ＡＶＥとＲＯＩ２＿ＡＶＥを算出する指標化処理手段を含んで構成されている。また、制御演算部１２は、歪み演算に係る１組のＲＦ信号フレームデータ（Ｎ）、（Ｘ）のフレーム間の時間で除して、単位時間当たりの歪み平均値ＡＶＥを算出し、グラフィック部８に転送して、色と歪みの大きさの関係を示すカラースケールを画面表示する。このとき、カラースキャンコンバートの基準となる平均値の位置を、図４に示すように、カラースケール上又はその隣に矢印等で表示する。なお、カラースケールは、図５に示すように、歪み平均値ＡＶＥのＸ倍をカラー値の上限値とし、それ以下を線形的にカラースキャンコンバートしてカラースケールを構成する。【００３８】
次いで、制御演算部１２は、歪み平均値ＲＯＩ１＿ＡＶＥとＲＯＩ２＿ＡＶＥから歪み比ＣＯＭＰを算出する。すなわち、歪み比ＣＯＭＰは、組織の弾性に相関する物理量（硬さ）であるＲＯＩ１とＲＯＩ２における歪みの平均値の比であり、ＲＯＩ１＿ＡＶＥとＲＯＩ２＿ＡＶＥを相互に指標化した指標値である。
（ステップＳ７）
次いで、制御演算部１２は、演算結果をグラフィック部８に転送して、図３、図４に示すように、表示窓４０を画面の右下部に表示させ、表示窓４０内に歪み平均値及び歪み比を表示する。、つまり、表示窓４０内に、ＲＯＩ１＿ＡＶＥ＝０．５％、ＲＯＩ２＿ＡＶＥ＝０．１％、ＣＯＭＰ＝０．２が表示される。【００３９】
ここで、関心領域ＲＯＩ１、ＲＯＩ２はそれぞれ独立した関心領域であるが、それらが図３のように重なった場合は、ＲＯＩ２の部分をＲＯＩ１から除外して、歪み及びその平均値を算出するようにすることができる。
上述では、断層画像中で生体組織の歪みが小さく（硬い）、腫瘍などの病変と思われる部位にＲＯＩ2を設定したが、カラースケールを基準に取りＲＯＩ2を自動的に設定してもよい。カラースケールに矢印４１を設定し、設定された矢印４１よりも硬い側の色相を検出する。そして、その検出された色相に対応する領域をＲＯＩ2として設定する。このように腫瘍部は硬い領域である性質を利用し、ＲＯＩ2を効率的に設定することができる。【００４０】
以上説明したように、本実施形態によれば、ＲＯＩ１、ＲＯＩ２の歪み比（硬さ比）を算出して画面表示しているから、関心を持った部位の組織の硬さを比べることにより、病変部などの硬さを定量的に診断することができる。例えば、癌は、脂肪に比べて１０倍以上の歪み比を持つことがわかっているので、ＲＯＩ１を脂肪層に、ＲＯＩ２を癌等の病変部と思われる部位に設定し、その比率を出すことで、癌の的確な診断の一手法として役に立つ。また、比較的安定して脂肪層が存在する部位、例えば体表から２〜３ｃｍの深度の一定領域をＲＯＩ１として自動設定してもよい。【００４１】
また、図４に示したように、カラースケールの横に、もしくは重畳してＲＯＩ１とＲＯＩ２内の歪みの平均値に対応する位置に、矢印等の指示マークを表示することによって、数値ではなく視覚的にＲＯＩ１とＲＯＩ２の歪みの相対関係を認識することができる。つまり、ＲＯＩ２内の歪み平均値の大きさと色との関係を詳細に知ることができる。【００４２】
また、ＲＯＩ２を点として設定すれば、その点がどの位の歪みであるかが、画面上に表示されるから、また、カラースケールでその歪みの大きさと、色に対応する位置に矢印で表されるので、ユーザーがより細かい歪みの詳細を知ることができる。
（実施形態２）
上記の実施形態では、硬さ比を比べる対象のＲＯＩを２つ設定する場合を説明したが、本発明はこれに限らず、制御演算部１２及びキーボード１３において、対比対象のＲＯＩをいくつでも発生させることができる。また、その場合、複数のＲＯＩの歪みを相互に指標値化することができる。【００４３】
図６を参照して、ＲＯＩを３つ設定した場合を説明する。図示のように、弾性像又は弾性像とＢモード像の合成画像上で、ＲＯＩ１と、ＲＯＩ１から横方向に離れた位置に、上下にＲＯＩ２、ＲＯＩ３を設定する。そして、それらのＲＯＩにおける歪み又はその平均値を相互に指標値化し、その指標値を画像表示器１１に表示する。【００４４】
この場合の指標値の例を、下記に示す。【００４５】
（１）ＲＯＩ間の歪み平均値の比
0003991282.tif 000002
（２）ＲＯＩの歪みの加算値の比
0003991282.tif 000003
（３）複数ＲＯＩの組合せの比
0003991282.tif 000004
ここで、（３）の複数ＲＯＩの組合せの比で指標値化する一例を、図７を用いて説明する。図７は、圧迫を付与する探触子２と弾性像との関係を模式的に示したものである。つまり、探触子２の広がり範囲に対応させて、画像の横方向にＲＯＩ１、ＲＯＩ３、ＲＯＩ２を離して設定し、真中のＲＯＩ３を病変部と思われる部位に設定した例である。【００４６】
探触子２により用手的に圧迫を加える場合、図７に示すように、探触子２の長手方向に加わる圧迫力Ｆ１〜Ｆ３が均一にならない場合がある。圧迫力の大きさに違いがあると、生体組織の硬さが同じでも歪みが異なる値となってしまうから、例えば、ＲＯＩ１とＲＯＩ２の組織が同じであっても、ＲＯＩ１／ＲＯＩ３の比と、ＲＯＩ２／ＲＯＩ３の比は、異なる値になってしまう。そこで、上記の（３）のように、Ｒａｔｉｏ（５）＝（ＲＯＩ１＋ＲＯＩ２）／ＲＯＩ３を指標とすることにより、圧迫力Ｆ１、Ｆ３の大きさに違いを補償することができる。
（実施形態３）
上記の実施形態１、２では、対比対象の複数のＲＯＩの歪みを指標値化して、画像表示器１１に表示する例を説明したが、本発明はこれに限らず、図８Ａ又は図８Ｂに示すように、ＲＯＩ１、ＲＯＩ２、ＲＯＩ３の歪み又は歪み比の時間変化をグラフにより表示することができる。なお、本実施形態のＲＯＩの配置は、図６と同じである。【００４７】
つまり、図８Ａ又は図８Ｂのように、時系列で計測される複数フレームの超音波計測データに基づいて歪み又は歪み比を求め、それらの時間変化をグラフ化して表示することにより、各ＲＯＩにおける硬さの違いを誤りなく認識できる。
図８Bに関し、詳細に説明する。ここでは説明のため、ＲＯＩ１ とＲＯＩ3を正常部とし、ＲＯＩ2を腫瘍部とする。ＲＯＩ１ ／ＲＯＩ3は、正常部同士の比であるので、歪み平均値の比の値は、ほぼ１となる。また腫瘍部の歪みは、正常部の歪に対して小さいため、ＲＯＩ2 ／ＲＯＩ3の歪み平均値の比の値は、1より小さい値となる。このように、正常部であるＲＯＩ3を予め定めておき、ＲＯＩ3に対する各領域の比較を行うことにより、各領域が正常部であるか腫瘍部であるかを判別することができる。また、この判定したい領域の設定は、キーボード１３等で任意に行うことができ、ＲＯＩを設定すると、ＲＯＩ3を基準とする歪み平均値の比に基づいて正常部、腫瘍部の判定を制御演算部１２で行う。この判定結果は文字、或いはＲＯＩ内をカラー表示して画面上に表示する。カラー表示する場合、例えば、判定したい領域が正常部と判定(比の値がほぼ１)されたならば、その領域は緑色として表示し、また腫瘍部と判定（比の値が１より小さい値）されたならば、その領域はオレンジ色として表示する。すなわち、制御演算部１２は、各ＲＯＩ相互の歪み平均値の比に基づいて、各ＲＯＩの組織が正常か否かを判定して判定結果を表示する判定手段を備えて構成される。
また、理論的には、これらの歪み比は時系列的に一定の値となる。歪み比が時系列的に変化している場合、圧迫等がエラーを起こしていることになるため、圧迫を再度行うよう操作者に促すことができる。
また、図８Ｃに示すように、フレームごとに求まる歪み比を加算して、加算フレーム数に対応させてその加算値をグラフ表示すれば、ＲＯＩ１／ＲＯＩ３とＲＯＩ２／ＲＯＩ３の比の違いが明確になる。【００４８】
つまり、本実施形態は、被検体１に圧迫を加えながら断層部位の超音波断層データを複数フレームにわたって計測し、その複数フレームの超音波断層データに基づいて断層部位における組織の弾性に相関する物理量である歪みの時間変化を求めている。そして、その歪みの時間変化に基づいて断層部位における弾性像の動画像を生成して表示装置に表示し、表示された弾性像の静止画像上に少なくとも２つのＲＯＩを設定し、各ＲＯＩにおける歪みの時間変化を表示するようにしているのである。これによれば、圧迫力の変化に応じて変化する歪みの各ＲＯＩの定量的な対比を的確に行うことができる。【００４９】
図９Aは、歪みグラフの歪み or 弾性率の表示軸において、ＲＯＩ1の歪みがＲＯＩ2、ＲＯＩ3に対して掛け離れているので、ＲＯＩ1とＲＯＩ2の間の歪みor弾性率の表示を省略して、表示幅を自動調整する形態を示すものである。この形態より、正常部と腫瘍部の比較であっても１画面内に歪みor弾性率を表示させることができる。【００５０】
また、図９Bは、歪みグラフの歪み or 弾性率の表示軸において、すべてのＲＯＩの歪みが中央に集中しているので、上下の歪みの表示を省略して、表示幅を自動調整する形態を示すものである。この形態より、比較範囲を拡大して、１画面内に歪みor弾性率を表示させることができる。【００５１】
図１０は、画像表示器１１に、弾性像３０と各ＲＯＩに適合した歪み平均値の比のグラフ３１を同一画面上に表示する形態を示す図である。弾性像３０を画面の上部に表示させ、グラフ３１を下部に表示させる。よって、弾性像３０と対応させてＲＯＩに適合した歪みの平均値の比を相対的に認識することができ、時系列的に比の値を確認しながら弾性像を観察することができる。幅調整矢印３２は、弾性像３０の上下間隔を調整するものである。幅調整矢印の幅を縮めると、弾性像３０が縮小され、逆に幅調整矢印の幅を広げると、弾性像３０は拡大される。すなわち、弾性情報演算手段である歪み演算部６は、時系列的に計測される複数の超音波断層データに基づいて複数の弾性像を生成し、生成した弾性像を画像表示器１１に表示された断層像と同一画面上に重ねて表示するとともに、ＲＯＩの歪みの平均値の比の変化を同一画面上に表示するようになっている。【００５２】
図１１は、図１０の形態に歪み平均値の比のグラフ３１に時相バー３５を表示させ、時相バー３５に対応する弾性像３０を表示させる形態、また所定間隔のループ再生を実現させる形態を示す図である。キーボード１３等で時相バー３５を横方向にずらすことにより、制御演算部１２は任意に時相に対応する弾性像３０を表示させる。例えば、ＲＯＩの歪み平均値の比のグラフに基づいて圧迫が最適な時相をグラフの傾き、極大点或いは極小点から探し、圧迫が最適な時相に対応する弾性像３０を表示させることができる。つまり、弾性情報演算手段である歪み演算部６は、時系列的に計測される複数の超音波断層データに基づいて複数の弾性像を生成し、画像表示器１１に表示された複数の弾性像の計測時相を指定する時相バー３５が入力手段によって操作されたとき、指定された計測時相に対応する弾性像を画像表示器１１に表示するようになっている。【００５３】
また、所定間隔の開始フレームと終了フレームを指定して、時間軸の表示を拡大してループ再生を行う。矢印３６は、開始フレームと終了フレームを示す記号である。時間軸ですべてのフレームを表示しきれない場合、時相バー３５で時間軸をスクロールする。
（実施形態４）
図１２Ａ〜図１２Ｃを参照して、本発明の弾性像表示方法の他の実施形態の特徴を説明する。実施形態１〜３で求めた歪みに関係する各指標値を色付け表示する場合のダイナミックレンジを、各指標値に応じて可変設定することにより、色付け表示の解像度を向上させることができる。【００５４】
すなわち、従来のダイナミックレンジは、図１２Ａに示すように、歪みが取り得る理論的な値の範囲の全域にカラー値を線形的に割り付けていた。これに対し、本実施形態では、例えば、図３のように２つのＲＯＩ１、ＲＯＩ２を設定した場合、各ＲＯＩの歪みの平均値の小値と大値を、図１２Ｂに示すように、それぞれカラー化の色相情報のダイナミックレンジの下限値と上限値に対応付ける。これにより、色付け表示の解像度を向上できる。また、例えば、図６のように３つのＲＯＩ１、ＲＯＩ２、ＲＯＩ３を設定した場合は、図１２Ｃに示すように、各ＲＯＩの歪みの平均値の最小値と最大値をダイナミックレンジの下限値と上限値に対応付ける。そして、例えば歪みの平均値が中間値を示すＲＯＩ３については、カラー値を適宜のレベルに設定して、ダイナミックレンジを非線形とすることにより、注目部位の色付け表示の解像度を向上させることができる。
（実施形態５）
図１３を参照して、本発明の弾性像表示方法のさらに他の実施形態の特徴を説明する。上記の実施形態１〜４では、同一断層部位における複数のＲＯＩの組織の硬さを対比する場合を説明した。しかし、本発明は、これに限られるものではなく、例えば、左右の手又は足、乳房などのように、被検体の対称性を有する左右の断層部位について別々に超音波断層データを計測し、左右の断層部位における弾性像を生成して対比することにより、異常部位などを診断することに適用できる。【００５５】
すなわち、図１３に示すように、左右の弾性像に関心領域（ＲＯＩ）として、それぞれ２つの層１、層２を設定する。そして、左右の２つの層１、層２の歪みの比を次のように求める。【００５６】
Ｌ＿ｒａｔｉｏ＝左層１／左層２
Ｒ＿ｒａｔｉｏ＝右層１／右層２
そして、さらに左右のＲＯＩの歪みの比のさらに比を次のように求めて表示するようにすることができる。【００５７】
ｃｏｍｐ＝Ｌ＿ｒａｔｉｏ／Ｌ＿ｒａｔｉｏ
本実施形態によれば、左右の手又は足、乳房などの対応する部位の組織の硬さの違いを定量的に評価して、診断の精度を高めることができる。
（実施形態６）
図１４、１５を用いて、本発明の超音波診断装置の他の実施形態を説明する。前述したように、生体組織の歪みは圧迫条件によって変動する。そこで、本実施形態では、被検体１の反発力によらず、同じ圧力、同じ速度で生体組織に圧迫をかける自動加圧装置を設け、再現性及び評価性に優れた歪み弾性像を得るようにしたものである。【００５８】
図１４に、被検体１に加える探触子２による圧迫の加減を行うことができる加圧システム５０の一例を示す。加圧システム５０は、ベッド５２等に動かないように固定されたバー５１と、バー５１に支持された探触子保持部５３を備えて構成されている。また、加圧システム５０は、検査者が被検体１に触れることなく、一つのボタンを操作することによって探触子２による圧迫の加減を行うことができるようになっている。特に、加圧システム５０は、被検体１からの反発力を考慮したものであり、組織弾性像を取得したい探触子２を探触子保持部５３に固定し、制御命令に応じて探触子保持部５３がステッピングモータ等で一定速かつ一定圧で探触子２を上下運動するように繰り返し制御するようになっている。【００５９】
図１５は、探触子２を保持する探触子保持部５３の概念構成を示す図である。図から明かなように、探触子保持部５３はバー５１に固定されている。探触子保持部５３には、ステッピングモータ等の駆動手段５５が設けられている。探触子２は、駆動手段５５にネジなどの把持手段で固定されている。探触子２を上下動させる速度と圧力は、レンジ切り換えによって可変できるようになっている。つまり、探触子保持部５３は加圧制御手段を構成し、バー５１は加圧制御手段を支持する固定手段を構成している。【００６０】
このように構成される加圧システム５０を超音波診断装置に連動させて、超音波診断装置の方で切り換えるように構成することができる。また、一回指令を出せばストップボタン又は加圧が発生するまで上下し続けるようにプログラミングすることも可能である。【００６１】
従来、提案されている自動加圧システムは、被検体１からの反発力を考慮していないため、生体が硬ければそれだけの反力をうけ、得られる数値データは加工しにくいものであった。【００６２】
この点、本実施形態によれば、被検体１の生体からの反発力を受けても、同じ圧力及び同じ速度で被検体１に圧迫をかけることができるので、再現性や評価性に優れた組織弾性像やそのデータを取得することができる。その結果、各部位に関して、一定の速度で加圧したときに得られる歪みの値を安定化することができるから、例えば、癌の良性、悪性を鑑別する際の精度が向上する。【００６３】
以上、本発明の弾性像表示方法の特徴について、組織の歪みを例に説明したが、本発明はこれに限られるものではなく、組織の弾性に相関する物理量として弾性率を用いることができるのはいうまでもない。
【図１】
0003991282.tif 000005
【図２】
0003991282.tif 000006
【図３】
0003991282.tif 000007
【図４】
0003991282.tif 000008
【図５】
0003991282.tif 000009
【図６】
0003991282.tif 000010
【図７】
0003991282.tif 000011
【図８Ａ】
0003991282.tif 000012
【図８Ｂ】
0003991282.tif 000013
【図８Ｃ】
0003991282.tif 000014
【図９Ａ】
0003991282.tif 000015
【図９Ｂ】
0003991282.tif 000016
【図１０】
0003991282.tif 000017
【図１１】
0003991282.tif 000018
【図１２Ａ】
0003991282.tif 000019
【図１２Ｂ】
0003991282.tif 000020
【図１２Ｃ】
0003991282.tif 000021
【図１３】
0003991282.tif 000022
【図１４】
0003991282.tif 000023
【図１５】
0003991282.tif 000024