WO2007083745A1

WO2007083745A1 - 弾性画像表示方法及び弾性画像表示装置

Info

Publication number: WO2007083745A1
Application number: PCT/JP2007/050796
Authority: WO
Inventors: Takeshi Matsumura
Original assignee: Hitachi Medical Corporation
Priority date: 2006-01-20
Filing date: 2007-01-19
Publication date: 2007-07-26
Also published as: JPWO2007083745A1; EP1980210B1; US20100220901A1; JP4919972B2; EP1980210A1; CN101370431A; US8098921B2; CN101370431B; EP1980210A4

Abstract

　【課題】　弾性画像のグレイゾーンにおける組織の悪性と良性の鑑別の精度を向上する。　【解決手段】　異なる複数の弾性情報を縦軸と横軸とする２次元マップに従って、各画素の色相又は輝度をそれぞれの弾性情報に応じて割り付け、それらの弾性情報を合成して一つの弾性画像を生成する。これにより、複数の弾性情報を一つの弾性画像に合成して表示できるから、一つの弾性画像を観察しながらリアルタイムでスキャン操作することができ、一つの弾性情報では良性又は悪性のグレイゾーンであっても他の弾性情報によれば良性又は悪性が明確な場合があり、確度の高い組織の鑑別をすることができる。しかも、それらの複数の弾性情報を合成した弾性画像によれば、視線を移すことなく診断することができ、視認性が向上する。

Description

明細書

弾性画像表示方法及び弾性画像表示装置

技術分野

[0001] 本発明は、弾性画像表示方法及び弾性画像表示装置に係り、具体的には、生体組織の弾性情報を画像ィ匕した弾性画像を観察して病変部の良性又は悪性を診断するのに好適な弾性画像表示方法及び弾性画像表示装置に関する。

背景技術

[0002] 生体組織の腫瘍等の病変部は、他の正常な組織に比べて組織が硬、ことから、超音波診断装置又は磁気共鳴撮像装置などにより生体組織の弾性を計測した弾性情報に基づいて、関心部位を含む生体断面の弾性画像を生成して病変部の診断に供している。関連する先行技術としては、特許文献 1、 2がある。

[0003] 弾性画像としては、例えば、超音波診断装置の場合は、関心部位に圧迫力を加えたときの生体組織の変位を求めて画像ィ匕した変位分布画像、その変位分布に基づ V、て求めた歪み画像、生体組織の各部の歪みと応力に基づ!/、て求めた弾性率画像等が一般に知られている。これらの弾性画像は、リアルタイムで得ることができるから、スクリーニング検査に適している。このような弾性画像は、通常、弾性情報に応じてカラーマッピングにより色付けされる。

[0004] 特許文献 1 :特開平 5-317313号公報

特許文献 2：特開 2000-60853号公報

発明の開示

発明が解決しょうとする課題

[0005] しかし、例えば、癌細胞は、生成された部位又は個人差によって、癌細胞の弾性と正常な生体組織 (例えば、繊維組織)の弾性とが近似することがあり、弾性に応じて力ラーマッピングされた弾性画像を観察しても、悪性と良性の鑑別を的確にできな、場合がある。

[0006] 例えば、乳腺疾患の場合、良性腫瘍は柔らカゝく弾性率が 10〜300 [kPa]で、悪性腫瘍は硬く弾性率が 250〜750 [kPa]であることが多いとされている。したがって、良性と悪性の弾性率が重なる 250〜300 [kPa]の範囲は、良性と悪性のいずれとも決めかねるグレイゾーンとなる。このような場合は、弾性率以外の弾性情報を計測し、それらの所見を総合的に判断して、関心部位の組織の良性又は悪性の診断を行うことが考えられる。弾性率以外の弾性情報としては、例えば、圧迫力を加えた生体組織が元に戻るまでの時間を考慮した粘弾性率を表す弾性画像、ある!/ヽは弾性率が歪みに対して有する非線形性を指標とした弾性画像などがある。

[0007] しかし、一つの弾性画像で確証が得られなヽ場合に、他の弾性画像に切り替えて同一の関心部位を再度確認する必要があるから、操作が煩雑となって検査に時間がかかり、検査効率が低下するとともに、患者に負担力かかるという問題がある。

[0008] また、複数の弾性画像を同時に取得して同一の画面に並べて表示したとしても、リアルタイムで関心部位を探索するスキャン操作にぉ、ては、複数の弾性画像を同時に観察することは不可能である。したがって、結局は、一つの弾性画像を観察しながらスキャン操作し、その弾性画像による関心部位がグレイゾーンの弾性情報に属する場合は、他方の弾性画像に視線を移して確認する必要があり、見まちがえなどが起きる可能性がある。また、複数の弾性画像を限られた画面上に表示する必要があるから、それぞれの弾性画像が小さくなるので、視認性が悪いという問題がある。

[0009] 本発明は、弾性画像のグレイゾーンにおける組織の悪性と良性の鑑別の精度を向上でき、リアルタイム性を損なわずに視認性に優れた弾性画像表示方法及び弾性画像表示装置を実現することを課題とする。

課題を解決するための手段

[0010] 上記の課題を解決するため、本発明の弾性画像表示方法及び弾性画像表示装置の第 1の態様は、異なる複数の弾性情報を一つの弾性画像に合成して表示することを特徴とする。これによれば、複数の弾性情報を合成して一つの弾性画像として表示されるから、その弾性画像を観察しながらリアルタイムでスキャン操作することができる。し力も、一つの弾性情報によれば良性又は悪性のグレイゾーンに属しても、他の弹性情報によれば良性又は悪性が明確なときは、確度の高い組織の鑑別をすることができる。特に、それらの複数の弾性情報を合成した弾性画像によれば、視線を移すことなく診断することができ、視認性が向上する。 [0011] また、本発明によれば、リアルタイムによる一度のスキャン操作で、確度の高い組織の鑑別をすることができるから、検査効率を著しく向上させることができる。また、再検查などの精密検査を省略できる可能性があり、患者の負担を大幅に軽減することができる。さらに、精密検査に適用して検査効率を向上させることができる。

[0012] 本発明の第 1の態様において、前記一つの弾性画像は、前記弾性情報の 2つを縦軸と横

軸とする 2次元マップに従って、各画素の色相又は輝度を前記 2つの弾性情報に応じて割り付けて生成することができる。また、これに代えて、前記一つの弾性画像は、異なる複数の前記弾性情報を悪性を評価する評価値に変換し、各画素の色相又は輝度を各弾性情報の評価値に応じて割り付けて生成することができる。

[0013] また、本発明の弾性画像表示方法及び超音波診断装置の第 2の態様は、設定される同一の関心領域における異なる複数の前記弾性情報を、数値と線図と図形の少なくとも一つにより対比した一つの弾性画像とすることができる。

[0014] この場合にお、て、前記一つの弾性画像は、前記関心領域における異なる複数の前記弾性情報を悪性の評価値に変換し、複数の前記弾性情報の評価値をバーチヤートにより対比させて表示する画像とすることができる。また、これに代えて、前記一つの弾性画像は、前記関心領域における異なる複数の弾性情報の時間変化を線図により対比させて表示する画像とすることができる。

[0015] ここで、本発明の異なる複数の弾性情報には、弾性率、粘弾性率、歪み量に対する弾性率の非線形性に係る非線形パラメタ、前記変位の演算結果に含まれる局所的な変位分布バラツキ (分散値)又は相関係数等の演算結果の信頼度に係る S/N判定情報、各計測点に加わる圧迫方向、変位、応力、歪み、ポアソン比、 2つの関心領域の歪みの比率などを適用できる。

発明の効果

[0016] 本発明によれば、弾性画像のグレイゾーンにおける組織の悪性と良性の鑑別の精度を向上でき、リアルタイム性を損なわずに視認性に優れた弾性画像表示方法及び弾性画像表示装置を実現することができる。

図面の簡単な説明 [0017] [図 1]本発明の弾性画像表示方法を実施する実施形態 1の超音波診断装置の構成図である。

[図 2]本発明の実施形態 1のカラー弾性像を生成する 2次元カラーマップの構成を説明する図である。

[図 3]乳腺組織の各部位における弾性率との関係の一例を示す図である。

[図 4]実施形態 1の実施例 1の 2次元カラーマップとカラー弾性画像を示す図である。

[図 5]歪み量に対する弾性率の非線形性を説明する図である。

[図 6]実施形態 1の実施例 2の 2次元カラーマップとカラー弾性画像を示す図である。

[図 7]変位バラツキを説明する図である。

[図 8]実施形態 1の実施例 3の 2次元カラーマップとカラー弾性画像を示す図である。

[図 9]実施形態 2の一実施例であり、複数の弾性情報を共通の指標値に変換して合成したカラー弾性画像を生成する説明図である。

[図 10]実施形態 3の一実施例を示し、複数の弾性情報を異なる 1次元カラーマップにより色付け、隣接する画素領域の弾性情報を交互に変えてモザイク状のカラー弾性画像を生成する説明図である。

[図 11]実施形態 3の他の実施例の複数の弾性情報を異なる 1次元カラーマップにより色付け、隣接する画素領域の弾性情報を交互に変えてモザイク状のカラー弾性画像を示す図である。

[図 12]設定された関心領域 ROIにおける複数の弾性情報を、数値、バーチャート長による表示、時間変化のグラフとして表示するようにした実施形態 4の説明図である。発明を実施するための最良の形態

[0018] 以下、本発明の弾性画像表示装置及び弾性画像表示方法を実施形態に基づ!/ヽて説明する。

(実施形態 1)

図 1に、本発明の弾性画像表示装置の一例として、例えば、超音波診断装置の一実施形態のブロック構成図を示す。図に示すように、被検体 1に当接して用いられる超音波の探触子 2は、被検体 1との間で超音波を送信及び受信する複数の振動子が整列された超音波送受信面を有して形成されている。探触子 2は、送信回路 3から供給される超音波パルスにより駆動される。送受信制御回路 4は、探触子 2の複数の振動子を駆動する超音波パルスの送信タイミングを制御して、被検体 1内に設定される焦点に向けて超音波ビームを形成するようになっている。また、送受信制御回路 4は、探触子 2の振動子の配列方向に電子的に超音波ビームを走査するようになって!/、る。

[0019] 一方、探触子 2は、被検体 1内から発生する反射エコー信号を受信して受信回路 5 に出力する。受信回路 5は、送受信制御回路 4から入力されるタイミング信号に従って、反射エコー信号を取り込んで増幅などの受信処理を行う。受信回路 5により受信処理された反射エコー信号は、整相加算回路 6において複数の振動子により受信された反射エコー信号の位相を合わせて加算することにより増幅される。整相加算回路 6 において整相加算された反射エコー信号の RF信号は、信号処理部 7に入力され、ゲイン補正、ログ圧縮、検波、輪郭強調、フィルタ処理等の信号処理がなされる。なお、整相加算回路 6において生成される RF信号は、複素復調した I、 Q信号であっても良い。

[0020] 信号処理部 7により処理された RF信号は白黒スキャンコンバータ 8に導かれ、ここにおいてディジタル信号に変換されるとともに、超音波ビームの走査面に対応した 2次元の断層像データに変換される。これらの信号処理部 7と白黒スキャンコンバータ 8によって断層像 (Bモード像)の画像再構成手段が構成される。白黒スキャンコンバータ 8 カゝら出力される断層像データは、切替加算部 9を介して画像表示器 10に供給されて B モード像が表示されるようになって!/、る。

[0021] 一方、整相加算回路 6から出力される RF信号は、 RF信号フレームデータ取得部 11 に導かれる。 RF信号フレームデータ取得部 11は、超音波ビームの走査面 (断層面)に対応する RF信号群を、フレームデータとして複数フレーム分を取得してメモリなどに格納する。変位計測部 12は、 RF信号フレームデータ取得部 11に格納されている取得時刻が異なる複数対のフレームデータを順次取り込み、取り込んだ一対のフレームデータに基づ、て断層面における複数の計測点の変位ベクトルを求め、変位フレームデータとして弾性情報演算部 13に出力するようになっている。

[0022] 本実施例の弾性情報演算部 13は、変位フレームデータに基づいて各計測点の生体糸且織の歪みを求めて歪みフレームデータを生成する歪み演算部と、歪みフレームデータに基づいて各計測点の生体組織の弾性率を求めて弾性率フレームデータを生成する弾性率演算部と、歪みや弾性率とは異なる他の弾性情報 (以下、関連弾性情報という。）を求めて関連弾性情報のフレームデータを生成する関連弾性情報演算部を有して構成されている。弾性情報演算部 13にて求められた歪みと弾性率と関連弾性情報のフレームデータは弾性情報処理部 14に出力されるようになっている。

[0023] 弾性情報処理部 14は、弾性情報演算部 13から入力される各弾性情報のフレームデータに対して、座標平面内におけるスムージング処理、コントラスト最適化処理、フレーム間における時間軸方向のスムージング処理などの様々な画像処理を施して、カラースキャンコンバータ 15に送出するようになって!/、る。

[0024] カラースキャンコンバータ 15は、本実施形態の特徴に係る異なる複数の弾性情報を合成した一つの弾性画像を生成する合成弾性画像生成部の機能を有する部分である。詳細は後述するが、弾性情報処理部 14により処理された弾性率と関連弾性情報のフレームデータを取り込み、設定された弾性率と関連弾性情報の 2次元カラーマツプに従って、フレームデータの画素ごとに色調コードを付与して弾性率と関連弹性情報を合成したカラー弾性像を生成するようになってヽる。

[0025] カラースキャンコンバータ 15により生成されたカラー弾性画像は、切替加算部 9を介して画像表示器 10に表示されるようになっている。また、切替加算部 9は、白黒スキヤンコンバータ 8から出力される白黒の断層像と、カラースキャンコンバータ 15から出力されるカラー弾性画像とを入力し、両画像を切り替えていずれか一方を表示させる機能と、両画像の一方を半透明にして加算合成して画像表示器 10に重ねて表示させる機能と、両画像を並べて表示させる機能を有して形成されている。なお、図示していないが、切替加算部 9から出力される画像データを格納するシネメモリ部を設け、装置制御インターフェイス部 18からの指令に従って、過去の画像データを呼び出して画像表示器 10に表示させることができる。

[0026] このように構成される本実施形態の基本的な動作につ!、て説明する。まず、探触子 2により被検体 1に加える圧力を変化させながら、被検体 1に超音波ビームを走査するとともに、走査面からの反射エコー信号を連続的に受信する。そして、整相加算回路 6から出力される RF信号に基づいて、信号処理部 7及び白黒スキャンコンバータ 8により断層像が再構成され、切替加算器 9を介して画像表示器 10に表示される。

[0027] 一方、 RF信号フレームデータ取得部 11は、被検体 1に加えられる圧迫力が変化する過程で、 RF信号を取り込んでフレームレートに同期させてフレームデータを繰り返し取得し、内蔵されたフレームメモリ内に時系列順に保存する。そして、取得時刻が異なる一対のフレームデータを単位として、連続的に複数対のフレームデータを選択して変位計測部 12に出力する。変位計測部 12は、選択された一対のフレームデータを 1次元もしくは 2次元相関処理し、走査面における各計測点の変位を計測して変位フレームデータを生成する。この変位ベクトルの検出法としては、例えば特開平 5-31 7313号公報等に記載されているブロックマッチング法又はグラジェント法が知られている。ブロックマッチング法は、画像を例えば N X N画素からなるブロックに分け、現フレーム中の着目しているブロックに最も近似しているブロックを前フレーム力も探索し、これに基づいて計測点の変位を求める。また、一対の RF信号フレームデータの同一領域における自己相関を計算して変位を算出することができる。

[0028] 変位計測部 12で求められた変位フレームデータは、弾性情報演算部 13に入力される。弾性情報演算部 13は、各計測点の歪み、弾性率、関連弾性情報といった予め設定された複数の弾性情報を演算して、必要な弾性情報フレームデータを弾性情報処理部 14に出力する。歪みの演算は、公知のように変位を空間微分することによって計算される。また、求めた歪みに基づいて各計測点の弾性率を演算する。弾性率を求める場合は、圧力計測部 17により計測された圧力の計測値を取り込み、これに基づいて各計測点における応力を演算する。圧力計側部 17は、探触子 2の超音波送受信面と被検体 1との間に設けられた圧力センサ 16により検出された圧力に基づいて、被検体 1内部の計測点における応力を演算する。

[0029] 弾性情報演算部 13の弾性率演算部は、各計測点における応力と、弾性情報演算部 13の歪み演算部で求めた歪みフレームデータから走査面上の各計測点の弾性率 E (例えば、ヤング率)を演算し、弾性情報処理部 14に出力する。ここで、図 3に、乳腺組織の各部位の弾性率を計測した一例を示す (出典は、 T.A.Krouskop et.al.Ultraso nic Imaging, 1998)。図において、横軸は乳腺組織の各部位を示し、縦軸は弾性率 [k Pa]を示している。図示のように、脂肪及び良性の乳腺組織は弾性率が十分小さく、浸潤性乳管癌の弾性率は概ね 380〜730となっている。一方、良性の繊維組織は概ね 130〜300であるのに対し、悪性である一部の乳管癌は概ね 220〜360程度であり、両者が重なる 220〜300の範囲は、良性と悪性が混在するグレイゾーンとなっている。したがって、弾性率のみを指標とすると、グレイゾーンの弾性率の部位については、良性と悪性を鑑別できない。

[0030] そこで、本実施形態では、弾性情報演算部 13の関連弾性情報演算部において、後述する各実施例で説明するように、歪みと弾性率とは異なる関連弾性情報を求めて弾性情報処理部 14に出力するようにしている。関連弾性情報としては、後述するように種々の弾性に相関する情報を適用できる。弾性情報処理部 14は、入力される弾性率及び関連弾性情報をスムージング処理などの処理を行って、カラースキャンコンパータ 15に出力する。

[0031] カラースキャンコンバータ 15は、異なる複数の弾性情報を合成して一つの弾性画像を生成する。カラースキャンコンバータ 15には、図 2に示すような 2次元カラーマップが設定されている。すなわち、縦軸に弾性率 E[kPa]、横軸に関連弾性情報を割り当て、縦軸は、下から上に向って弾性率 Eが大きくなるように、つまり下から上に向うにつれて生体組織の良性力も悪性を評価できるようになつている。また、横軸は、左から右に向って関連弾性情報が大きくなるように、つまり左力右に向うにつれ生体組織の良性力も悪性を評価できるようになつている。また、縦軸及び横軸に沿って例えば 2 56階調化による色調のグラデーションが設定されている。また、色調のグラデーションは、縦軸の左右軸と、横軸の上下軸のそれぞれに異なる色調のグラデーションが設定されている。図示例では、弾性率の左軸は「黄〜緑」のグラデーション、右軸は「赤〜青」のグラデーション、関連弾性情報の下軸は「黄〜赤」のグラデーション、上軸は「緑〜青」のグラデーションが設定されている。つまり、縦軸方向、横軸方向に限らず、斜め方向を含む任意の方向にグラデーションがかかり、 2次元マップの 2次元平面内に割り当てられた色相はすべての点において固有の色相を持っている。

[0032] したがって、図 2の 2次元カラーマップによれば、図中に矢印で示したように、右上隅が弾性率及び関連弾性情報ともに悪性を示し、左下隅が弾性率及び関連弾性情報ともに良性を示すことになる。また、良性又は悪性を判断しにくい弾性率 Eのグレイゾーンは領域 R21, R22, R23、関連弾性情報のグレイゾーンは R12, R22, R32となっている。しかし、弾性率と関連弾性情報の両方の評価を考慮すれば、領域 R21、 R12は良性の可能性が高ぐ領域 R23、 R32は悪性の可能性が高いことを鑑別でき、グレイゾーンを領域 R22に限定できる。

[0033] ここで、カラースキャンコンバータ 15において図 2の 2次元カラーマップを用いて弾性率と関連弾性情報の 2つの弾性情報を合成したカラー弾性画像を生成する手順を説明する。カラースキャンコンバータ 15は、弾性情報処理部 14によりスムージングなどの処理がされた弾性率と関連弾性情報のフレームデータを取り込み、 2次元カラーマツプ上で画素単位ごとにフレームデータの弾性率と関連弾性情報が交叉する位置の色調コードを読み取り、各計測点に対応する各画素の色調として設定してカラー弾性画像を生成する。

[0034] 上述したように、本実施形態によれば、 1つの弾性情報では鑑別できな力つたグレイゾーンを、関連弾性情報を考慮することによって大幅に低減できる。その結果、弾性率だけでは鑑別できな力つたグレイゾーンの良悪性の鑑別の確度を向上させることができる。

[0035] なお、カラースキャンコンバータ 15に代えて、白黒スキャンコンバータを用いることができる。この場合は、弾性率と関連弾性情報が大きい領域は輝度を明るぐ逆に小さい領域は輝度を暗くするなどにより、良性又は悪性を鑑別できるようにすることができる。

[0036] 以下に、本実施形態の特徴部である弾性情報演算部 13とカラースキャンコンパ一タ 15の詳細構成と動作を、関連弾性情報の具体的な実施例に基づいて説明する。実施例 1

[0037] 本発明の関連性情報として、歪み量に対する弾性率の非線形性を表す非線形パラメタを適用した実施例 1について説明する。本実施例の弾性情報演算部 13を構成する関連弾性情報演算部は、図 4(A)に示した 2次元カラーマップのように、関心部位の生体組織の弾性特性の非線形性を表す非線形パラメタ αを求めるようになつている [0038] すなわち、一般に、図 5に示すように、生体組織は圧迫の圧力もしくは歪み量が増すほど硬くなる特性があり、脂肪組織は大きな歪み量の範囲までほぼ一定の弾性率が計測される線形性を示すのに対し、繊維組織や浸潤性腫瘍は歪み量の増加に伴い弾性率も顕著に増加し、歪み硬化の現象が認められている。また、繊維組織よりも浸潤性の悪性腫瘍の方が歪み硬化の程度が大きぐ後述する非線形パラメタも大きくなる。ここで、歪み量とは、圧迫開始力計測時点までの歪みの積分値であり、歪みフレームデータの歪みとは区別して用いる。

[0039] 図 5に示すように、一般に生体組織の歪み量∑ ε [%] 弾性率 E [kPa]の関係は、非線形の弾性応答を示す。そこで、歪みフレームデータに基づいて圧迫開始力も計測時点までの歪み量を求め、上述した弾性率フレームデータに基づいて歪み量弾性率の関係データを取得する。そして、歪み量弾性率の非線形曲線を関数で近似し、最小二乗法などにより最もよく近似する曲線に基づいて非線形パラメタ OCを求めて、その組織の非線形性を表す情報として評価することが提案されて、る (例えば、電子情報通信学会論文誌 2001/12 Vol.J84-A No.12、「超音波による組織の非線形弾性特性の画像化」新田他)。

[0040] しかし、従来は、弾性率 Eとしてヤング率を用い、歪み量∑ εに対する弾性率 Εの関係の曲線の傾きを 1次のパラメタで仮定して、ることから非線形性が顕著な実際の組織の場合には、定量情報とは成り得ないため、確定診断に適用できない。

[0041] そこで、本実施例では、組織の非線形応答を忠実に近似する曲線として、次式 (1) に示す高次の関数を適用する。

Ε = Ε0 + α Χ ( Σ ε ) ^ν ( γは、 2以上の自然数） (1)

[0042] 式 (1)の非線形パラメタ aは、曲線の傾きの変化に相当する値であるから、組織の非線形性が大きいほど大きな値となり、線形に近づくほど限りなくゼロに近くなる。なお、式 (1)をさらに一般ィ匕して次式 (2)により指数関数型で近似することができる。また、非線形パラメタ OCが非線形性を表すものと設定することが可能である。

E = (E0 - l) + exp( a X∑ ε ) (2)

[0043] 式 (1)又は (2)で示される近似関数によれば、圧迫の程度にかかわらず、任意の計測区間において非線形パラメタ αを評価しても、高精度で確定した非線形パラメタ αを得ることができ、定量的な評価が実現できる。つまり、現時刻 t(0)において計測された歪み量—弾性率の関係 (∑ ε (0)、 Υ(0))と、過去の時刻 t(— 1)、 t(— 2)、 · · ·において得られた歪み量弾性率の関係 (∑ ε (一 1)、 Υ (— 1》、 (∑ ε (一 2)、 Υ (— 2》、…を歪み量—弾性率の関係の蓄積データとして適用し、現時刻 t(0)にお、てリアルタイムに非線形パラメタ αを評価しても、その非線形パラメタの数値情報、画像情報は圧迫程度に依存しない客観的な情報となる。

[0044] したがって、非線形パラメタ aをリアルタイムの弾性画像として表示することが可能であり、弾性率の弾性画像に良性と悪性が不確かなグレイゾーンであっても、非線形ノラメタ OCの弾性画像を合成すれば、良悪性を確度高く鑑別することができる。

[0045] 図 4に、本実施例により生成されたカラー弾性画像の一例を示す。同図 (A)は 2次元カラーマップの設定例であり、同図 (B)はその 2次元カラーマップを用いて作成した力ラー弾性画像である。図 4(A)に示すように、本実施例の 2次元カラーマップは、縦軸に弾性率 E [kPa]力横軸に非線形パラメタ (Xが割り当てられている。また、縦軸と横軸の色調のグラデーションは図 2と同様に設定される。

[0046] このような 2次元カラーマップを用い、カラースキャンコンバータ 15において弾性率と非線形パラメタの両方の弾性情報を合成したカラー弾性画像の一例は、図 4(B)のようになる。同図において、組織 3は 2次元カラーマップ上の P3の位置の色相が与えられ、弾性率 Eが高ぐ非線形パラメタ αも高いので、悪性を強く疑うことができる。一方、組織 4は、 Ρ4の位置の色相が与えられ、弾性率 Εはグレイゾーンの値である力非線形パラメタ αは十分大きいことから、良性よりは悪性の疑いが促される。その他の組織 1、 2、 5については、良性と判断される。

[0047] 以上説明したように、本実施例によれば、 1つの弾性情報では鑑別できな力つたグレイゾーンを大幅に低減できるから、良性又は悪性の鑑別の確度を向上させることができる。特に、弾性率の画像と、非線形パラメタの画像の 2つの画像を見比べる必要が無ぐ 1つの弾性画像で良性又は悪性の鑑別の確度を向上させることができるから、視認性及び検査効率を向上できる。

実施例 2

[0048] 図 6(A)に、弾性率 Εの他に、関連弾性情報として演算結果の信頼度に係る S/N判定情報を適用した実施例 2の 2次元カラーマップを示し、図 6(B)に、弾性率と S/N判定情報の 2つの弾性情報に基づいて合成したカラー弾性画像の一例を示す。

[0049] 本実施例の弾性情報演算部 13の関連弾性情報演算部は、 S/N判定情報演算部を備えている。 S/N判定情報演算部は、変位の演算結果に含まれる変位分布の局所的なバラツキ (分散)、又は変位計測部 12の変位演算で得られる相関係数により表される S/N判定情報を求める。

[0050] 変位計測部 12で求められた変位フレームデータは、図 7(A)、（B)に示すように、各計測点における 2次元の変位ベクトル (Xi,j、 Yi,j)で表される。つまり、画像縦軸方向 (y方向)成分と画像横軸方向 (X方向)成分の要素データ群を、それぞれ、

Xi,j(t) (ί=1,2,3，· · ·、 N、】=1,2,3, · · ·Μ)

Yi,j(t) (ί=1,2,3，· · ·、 Ν、』'=1,2,3, · · ·Μ)

として表す。添え字の iは横軸方向の座標を、 jは縦軸方向の座標を示している。 S/N 判定情報演算部は、図 7の変位フレームデータに基づいて、例えば、本実施例の特徴の S/N判定情報の一つである変位バラツキのフレームデータを生成するようになつている。

[0051] つまり、被検体 1に一定方向の圧迫力を加えたときの生体各部の組織の動的振舞いは、（A)非常に硬い充実性組織 (例えば、癌腫瘍等)、（B)柔らかい充実性組織 (例えば、繊維腺腫)、あるいは (C)流動的な嚢胞組織 (例えば、嚢胞)等に応じて顕著な相違がある。例えば、（A)の癌腫瘍等を形成する組織は、細胞の密度及び組織の結合度が高いことから纏まって同じ動きをする。これに対し、癌腫瘍等の周辺の組織は必ずしも組織結合度が高くないので、癌腫瘍等に比べて動きの大きさや方向にバラッキを有する。

[0052] このような組織の動的な振舞、の違、を変位のバラツキ度合として評価することによって、癌腫瘍等の存在を鑑別したり、癌腫瘍等の大きさを鑑別することができる。また、（B)の繊維腺腫のように比較的柔らかい組織の場合は、圧迫力を受けた位置に応じて変位の大きさや方向が分散して、圧迫方向に対して傾斜を持つバラツキが現れ、圧迫の中心から外れるほどその傾きが大きくなることにより識別できる。さらに、（C) の嚢胞のような非細胞組織の場合は、胞内には組織の結合がないので、各部位の動きが完全に自由であるから、変位の大きさや方向のバラツキが顕著に現れた分布になる。そこで、変位バラツキ演算部において、各計測点の変位の大きさ及び方向を示す変位バラツキを演算して、その変位バラツキフレームデータを関連弾性情報として生成する。

[0053] また、変位は圧迫力の大きさに相関するから、変位の大きさ及び方向を示す変位バラツキは、組織が同じであっても、異なる組織として誤認する場合が生ずる。そこで、変位の大きさ及び方向を定量的な値にするため、各計測点を取巻く局所領域に含まれる複数点の変位の平均値で正規ィ匕して、変位バラツキ度合を求めることが好ましい。これにより、関心部位とその周囲の組織との結合関係を客観的、定量的に評価することが可能となり、例えば、組織の充実性ゃ嚢胞性などの組織性状の識別が容易になるなど、組織鑑別を行うための新たな弾性情報とすることができる。また、変位ノラツキ度合は、局所的な変位バラツキを規格ィ匕したものであるから、組織性状を圧迫の方法に依存しな、定量的な評価をすることができ、 S/N判定情報としての変位バラツキの逆数 (変位バラツキが大き、ほど信頼性が低、ため)が、ある設定閾値より小さ、計測点にぉ、ては、圧迫前後にお、て計測点断面がずれてしまった場合であるなど、その計測点に何らかの理由により変位演算にノイズが含まれていると判断し、変位データを利用した弾性率やその他の弾性データの信頼性が劣るものとして扱うことができる。

[0054] なお、 S/N判定情報演算部は、変位バラツキに代えて、 S/N判定情報の他の一例である相関係数を演算することができる。この相関係数は、変位分布を相関法により演算する場合に得られる係数である。相関法による変位分布演算は、圧迫前における組織部位の RF信号波形によりテンプレートを設定し、そのテンプレートに最も近似した波形を有する領域を圧迫後の RF信号を探索する。そして、設定された探索範囲において、最も大きな相関係数が演算された点が、圧迫後におけるテンプレートの行き先として認識され、変位を求める方法である。したがって、相関係数の絶対値は、テンプレートのマッチングの度合いを表す値であり、変位演算の信頼度が高いほど「1」に近い値となり、信頼度が低いほど「0」に近い値となる。そこで、相関係数が設定閾値 (例えば、 0.8)を越えない計測点は、その計測点に何らかの理由により変位演算にノイズが含まれてヽると判断し、変位データを利用した弾性率やその他の弾性データの信頼性が劣るものとして扱うことができる。

[0055] 図 6(A)、（B)は、このようにして求められた変位バラツキ (の逆数)又は相関係数に基づ!、た S/N判定情報を信頼性の尺度として、弾性率と S/N判定情報を 2軸とする 2次元カラーマップを設定してカラー弾性画像を生成した例を示す。同図 (A)に示すように、本実施例の 2次元カラーマップは、縦軸を弾性率の大きさ、横軸を S/N判定情報力もみた弾性率の信頼性としている。弾性率は、図 2と同様に設定され、 S/N判定情報は左力も右に向うにつれて信頼性低力も信頼性高に設定されている。また、領域 R 25、 R26は弾性率のグレイゾーンであり、 R15、 R25、 R35は S/N判定情報のグレイゾーンである。しかし、領域 R26は S/N判定情報力見て信頼性が高いので、弾性率のグレイゾーン力も外して評価することが可能になる。

[0056] 例えば、図 6(B)に示すように、組織 3は 2次元カラーマップ上の P3の位置の色相が与えられ、弾性率 Eが高ぐ S/N判定情報の信頼性も高いので、悪性を強く疑うことができる。一方、組織 4は、組織 3と同程度に弾性率が高いが、 S/N判定情報の信頼性が低いので、 P4の位置の色相が与えられ、弾性率 Eはグレイゾーンではないが、信頼性が疑われるので、再検査が促される。

[0057] 本実施例によれば、実施例 1と同様に、 1つの弾性情報では鑑別できな力つたダレィゾーンを大幅に低減できるから、良性又は悪性の鑑別の確度を向上させることができる。特に、弾性率の画像と、 S/N判定情報の画像の 2つの画像を見比べる必要が無ぐ 1つの弾性画像で良性又は悪性の鑑別の確度を向上させることができるから、視認性及び検査効率を向上できる。

実施例 3

[0058] 図 8(A)に、弾性率 Eの他に、関連弾性情報として圧迫方向に対する角度を適用した実施例 3の 2次元カラーマップを示し、図 8(B)に、圧迫方向に対する角度の説明図を示す。図 8(B)に示すように、圧迫を加える探触子の真下の方向 (圧迫方向)は適切に圧迫が加わる。しかし、圧迫方向から角度 Θ外れた位置には適正な圧迫が加わらな、から、その位置の変位計測データに基づく弾性情報は適正でな、ことが予想される。特に、前立腺癌の計測の場合は、曲率の大きな円柱状の探触子を用いるから、全視野方向に対して同等の圧迫を加えることは困難である。したがって、同一の計測断面において、圧迫が適切な領域と、不適切な領域が存在する。

[0059] そこで、本実施例の 2次元カラーマップは、縦軸を弾性率の大きさ、横軸を圧迫方向に対する角度 0 (— π Ζ2〜+ π Ζ2)としている。このように設定することにより、注目している関心部位における弾性率の値が、どの程度適正な圧迫条件下で計測された結果であるかを一目瞭然に評価することができる。したがって、不適正な Θの範囲で計測された弾性率は確度が低ヽので、疑わしヽ場合は再検査を促すことができる。具体的には、 Ρ3の色相が与えられた領域 R32の組織は、適正な圧迫条件の下で計測された弾性率であることがわかる。しかし、 Ρ4の色相が与えられた領域 R33の組織は、弾性率が高くて悪性が疑われるが、不適正な圧迫条件下で計測された値であることがわ力る。したがって、本実施例のカラー画像を観察し、 R33の領域の色相が付与された関心部位については、リアルタイム中において、再度その関心部位に向つて圧迫することにより、適正な圧迫下で良悪性を鑑別することが可能になる。

したがって、本実施例によれば、実施例 2と同一の効果を奏することができる。実施例 4

[0060] 上記の実施例 1〜3は、それぞれ、弾性率非線形パラメタ (X、弾性率 S/N判定情報、弾性率圧迫方向の 2次元カラーマップを用いて 2つの弾性情報を合成した力ラー弾性画像を生成する例について説明したが、本発明はこれに限らない。例えば、弾性率と組み合わせて、粘弾性率、変位、応力、歪み、ポアソン比、 2つの関心領域の歪みの比率などの弾性情報を演算して、 2次元カラーマップを設定することができ、実施例 1〜3と同一の効果を奏することができる。

また、弾性率との組み合わせに限ることなぐ上記の弾性情報を任意に組み合わせて用いることができ、実施例 1〜3と同一の効果を奏することができる。

[0061] (実施形態 2)

上記の実施形態 1では、 2次元カラーマップを用いて 2つの弾性情報を合成したカラ一弾性画像を生成する本発明を説明したが、本発明はこれに限らない。すなわち、図 1の弾性情報演算部 13にて、複数の異なる弾性情報を演算してフレームデータを生成し、カラースキャンコンバータ 15にて、それらの弾性情報のフレームデータを重ねて一つのカラー弾性画像を生成するようにすることができる。

[0062] 本実施形態によれば、実施形態 1と同様に、 1つの弾性情報では鑑別できなかったグレイゾーンを大幅に低減できるから、良性又は悪性の鑑別の確度を向上させることができる。特に、弾性率の画像と、他の関連弾性情報の画像の 2つの画像を見比べる必要が無ぐ 1つの弾性画像で良性又は悪性の鑑別の確度を向上させることができるから、視認性及び検査効率を向上できる。

[0063] 実施形態 2を適用して生成するカラー弾性画像の一例を図 9を用いて説明する。本実施例は、弾性情報演算部 13にて、弾性率と、粘弾性率と、非線形パラメタ等の複数の弾性情報のフレームデータを演算するとともに、各弾性情報のフレームデータの計測点ごとに、組織の悪性の確率を (悪性度)を求める。ここで、粘弾性率とは、圧迫力を加えた生体組織が元に戻るまでの時間を考慮した弾性率である。これにより弹性率による悪性確率画像と、粘弾性率による悪性確率画像と、非線形パラメタによる悪性確率画像の一例を、それぞれ図 9(A)〜(C)に示す。そして、弾性情報演算部 13 又はカラースキャンコンバータ 15にて、複数の悪性確率画像 x(ここで、 x= l、 2、 3、 · · の計測点 G、 j)ごとに悪性確率 Pij(x)を求める。ここで、 1= 1, 2, 3, · · · , N、 j= l, 2, 3, · · · , Mである。

[0064] さらに、計測点 (i、 j)ごとに、複数の悪性確率画像の悪性確率 Pij(x)の積を総合的な診断による悪性確率 Pijとして、下式により求める。

Pij = Pij(l) X Pij(2) X Pij(3) X · · ·

求めた値 Pij(x)を画素値として図 9の左側に示すカラーバーに設定された色相を付与して、図 9(D)に示すカラー弾性画像を生成する。

[0065] 本実施例によれば、異なる複数の弾性情報を悪性を評価する評価値に変換し、各画素の色相 (又は輝度)を各弾性情報の評価値に応じて割り付けて弾性画像を合成して、ることから、 1つの弾性情報では鑑別できな力つたグレイゾーンを大幅に低減できる。その結果、良性又は悪性の鑑別の確度を向上させることができる。特に、複数の弾性画像を見比べる必要が無ぐ 1つのカラー弾性画像で良性又は悪性の鑑別の確度を向上させることができるから、視認性及び検査効率を向上できる。

[0066] なお、本実施例では、計測点 (i, j)ごとの悪性確率 Pij(x)の積を求めて画素値とする例を示したが、これに代えて、悪性確率 Pijを多次元ベクトルの長さとして、下式により求める悪性確率 Pijを指標とすることができる。

Pij = ^(Pij(l) X Pij(l) + Pij(2) X Pij(2)

+ Pij(3) X Pij(3) X - - -)

さらに、これに代えて、下式に示すように、悪性確率 Pijを悪性確率 Pij(x)の平均値とすることができる。

Pij = (Pij(l) + Pij(2) + Pij(3) + - · Z (弾性画像の個数）

[0067] (実施形態 3)

実施形態 1では、 2次元カラーマップを用いて 2つの弾性情報を合成したカラー弾性画像を生成する本発明を説明したが、本実施形態では、弾性情報ごとに 1次元カラ一マップを設定し、それら 2つのカラーマップを例えば暖色系と寒色系など、互いに混合しない、つまり相容れない色相のグラデーションを設定する。そして、画素単位又は隣接する複数の画素グループ単位で交互に 2つのカラーマップで色付けして、モザイク状のカラー弾性画像を合成することを特徴とする。

[0068] 図 10に、本実施形態により生成されるモザイク状のカラー弾性画像の一例を示す。

図示のように、弾性率 E[kPa]と粘弾性率 7? [kPa's]ごとに 1次元カラーマップが設定されている。それらのカラーマップは、いずれも下から上に向って良性力も悪性に変化させて対応付けられ、図示では明らかでないが、例えば暖色系と寒色系の互いに混合しない、つまり相容れない色相のグラデーションが設定されている。また、それらのカラーマップの中間値の近傍は、良悪性の判別が困難なグレイゾーンとなって!/、る

[0069] このようなカラーマップを用い、例えば、隣接する複数の画素グループ単位で交互に 2つのカラーマップで色付けして、モザイク状のカラー弾性画像を合成すると、同図右側のカラー弾性画像が生成される。図から明らかなように、組織 3は、弾性率及び粘弾性率は強く悪性が疑われることから、弾性率と粘弾性率の悪性を表す異なる色相がモザイク状になって表示される。一方、組織 4は、グレイゾーンの領域であり、診断に迷うが、粘弾性率では強く悪性を疑えることが認識できる。

[0070] 図 11に、本実施形態により生成されるモザイク状のカラー弾性画像の他の一例を示す。図示のように、弾性率 E[kPa]と粘弾性率 7? [kPa's]ごとに 1次元カラーマップが設定されている点は、図 10と同じである。しかし、それらのカラーマップは、例えば、良性方向の色は同一 (例えば、黄色)に設定し、悪性方向は識別し易くするように、相反する色相を割り当てて設定されている。図示では明らかでないが、それらの色相にはグラデーションが設定されている。また、それらのカラーマップの中間値の近傍は、良悪性の判別が困難なグレイゾーンとなっている。

[0071] このようなカラーマップを用い、例えば、隣接する複数の画素グループ単位で交互に 2つのカラーマップで色付けして、モザイク状のカラー弾性画像を合成すると、同図右側のカラー弾性画像が生成される。図から明らかなように、組織 4の領域のモザイク模様は色が黄色であまり色相に差がないから、弾性率及び粘弾性率共に、良性であることを直ちに認識できる。また、組織 3は、隣接する画素グループの色が大きく異なつたモザイク模様となることから、総合的に判断して強く悪性を疑うべきであることが認識できる。

[0072] (実施形態 4)

上述の実施形態 1では、 2次元カラーマップにより 2つの弾性情報を合成して一つのカラー弾性画像を生成する態様、実施形態 2では複数の異なる弾性情報のフレームデータを重ねて一つのカラー弾性画像を生成する態様、実施形態 3では弾性情報ごとに 1次元カラーマップを設定し、画素単位又は隣接する複数の画素グループ単位で交互に色付けしてモザイク状のカラー弾性画像を合成する態様について説明した

[0073] しかし、本発明の弾性画像表示方法は、それらに限られるものではなぐ図 12に示すカラー弾性画像を合成する実施形態 4をも包含している。すなわち、同図に示すように、本実施形態は、 1つの弾性画像又は Bモード断層像において、関心領域 ROIを設定すると、その ROIにおける複数の弾性情報の例えば平均値を、数値、バーチヤート長による表示、時間変化のグラフとして表示するようにした例である。図示例では、組織 3に設定された ROI内の弾性率と粘弾性率の平均値を演算し、それらの値をリアルタイムで、数値表示、バーチャート長表示、時間変化のグラフ表示するようにしている。また、弾性画像などの表示画像には、弾性率などの弾性情報に応じた色調の色を付して表示するようにしてヽる。

[0074] すなわち、本実施形態によれば、関心領域における異なる複数の弾性情報を悪性の評価値に変換し、複数の弾性情報の評価値をバーチャートにより対比させて表示して、るから、一目で複数の弾性情報を対比して良悪性の鑑別を確度高く行うことができる。同一の関心領域における異なる複数の弾性情報を、数値と線図と図形の少なくとも一つにより対比した場合も、同様である。

[0075] 以上、本発明を実施形態 1〜4により説明したが、本発明の弾性画像表示方法は、それらに限られるものではなぐ弾性率等に組み合わせる他の弾性情報には、 Bモード像、磁気共鳴 (MR)画像、 CT画像、血流カラードプラ像、ティシュドプラ像などを適用できる。 Bモード像や MR画像は弾性の高い部位が例えば暗く表示されるから、弹性情報として利用できる。また、血流カラードプラ像は、悪性の腫瘍 (癌細胞など)は血管の密度が高いことから、血流情報を弾性情報として利用できる。ティシュドプラ像は、例えば、心筋などの弾性を診断する際に弾性情報として利用できる。

[0076] また、本発明は、超音波診断装置上で実現するのに限られるものではなぐ CT装置や MR装置上で実現することができる。さらに、それらの診断装置力も分離して、パ一ソナルコンピュータ上で実現するようにしてもょ、。

[0077] また、 2次元カラーマップを例として説明した力弾性率、粘弾性率、歪み量に対する弾性率の非線形性に係る非線形パラメタ、変位の演算結果の信頼度に係る S/N判定情報、各計測点に加わる圧迫方向、変位、応力、歪み、ポアソン比等から 3つ (又は 3つ以上)を選択し、 3次元 (又は高次元)カラーマップを作成して、本発明を実施することができる。選択する情報が多ければ多いほど、グレーゾーンは減少する。

符号の説明

[0078] 1 被検体

2 探触子

3 送信回路

4 送受信制御回路

5 受信回路

6 整相加算回路信号処理部

白黒スキャンコンバータ切替加算部

画像表示器

RF信号フレームデータ取得部変位計測部

弾性情報演算部

カラースキャンコンバータ圧力センサ

圧力計測部

Claims

請求の範囲

[1] 被検体から得られた受信信号のフレームデータを取得し、取得時刻が異なる一対の前記フレームデータに基づ、て複数の計測点における変位を求め、前記各計測点の変位に基づいて前記各計測点における異なる複数の弾性情報を求め、求められた複数の弾性情報に基づいて一つの弾性画像を生成して表示画面に表示する弾性画像表示方法。

[2] 請求項 1に記載の弾性画像表示方法にお!、て、

前記複数の弾性情報は、弾性率、粘弾性率、歪み量に対する弾性率の非線形性に係る非線形パラメタ、前記変位の演算結果に含まれる局所的なバラツキ又は相関係数等の演算結果の信頼度に係る SZN判定情報、前記各計測点に加わる圧迫方向、変位、応力、歪み、ポアソン比力選択される少なくとも 2つであることを特徴とする弾性画像表示方法。

[3] 請求項 1に記載の弾性画像表示方法にお!、て、

前記一つの弾性画像は、前記弾性情報の 2つを縦軸と横軸とする 2次元マップに従って、各画素の色相又は輝度を前記 2つの弾性情報に応じて割り付けて生成されることを特徴とする弾性画像表示方法。

[4] 請求項 3に記載の弾性画像表示方法において、

前記 2次元マップは、 2次元カラーマップであり、前記縦軸の左右軸と前記横軸の上下軸は、それぞれ異なる色調のグラデーションが設定されて、ることを特徴とする弾性画像表示方法。

[5] 請求項 4に記載の弾性画像表示方法におヽて、

前記縦軸の弾性情報は弾性率であり、前記横軸の弾性情報は前記歪み量に対する弾性率の非線形性に係る非線形パラメタであることを特徴とする弾性画像表示方法。

[6] 請求項 4に記載の弾性画像表示方法におヽて、

前記縦軸の弾性情報は弾性率であり、前記横軸の弾性情報は粘弾性率であることを特徴とする弾性画像表示方法。

[7] 請求項 4に記載の弾性画像表示方法におヽて、前記縦軸の弾性情報は弾性率であり、前記横軸の弾性情報は前記各計測点にカロわる圧迫方向であることを特徴とする弾性画像表示方法。

[8] 請求項 1に記載の弾性画像表示方法にお!、て、

前記一つの弾性画像は、前記各計測点を隣接する複数の計測点を有する複数の設定領域に分けて、該各設定領域における 2つの前記弾性情報の代表値を求め、 2 つの前記弾性情報に対してそれぞれ異なる 2つの色相を割り付けてなる 2つの 1次元カラーマップに従って、前記各設定領域の各画素に色付けして生成されることを特徴とする弾性画像表示方法。

[9] 請求項 8に記載の弾性画像表示方法にぉ、て、

前記 2つの前記弾性情報は、弾性率と粘弾性率であることを特徴とする弾性画像表示方法。

[10] 請求項 1に記載の弾性画像表示方法にぉ、て、

前記一つの弾性画像は、異なる複数の前記弾性情報を悪性を評価する評価値に変換し、各画素の色相又は輝度を各弾性情報の前記評価値に応じて割り付けて生成されることを特徴とする弾性画像表示方法。

[11] 請求項 1に記載の弾性画像表示方法にぉ、て、

前記一つの弾性画像は、設定される同一の関心領域における異なる複数の前記弾性情報を、数値と線図と図形の少なくとも一つにより対比して表した弾性画像であることを特徴とする弾性画像表示方法。

[12] 請求項 11に記載の弾性画像表示方法にぉ、て、

前記一つの弾性画像は、前記関心領域における異なる複数の前記弾性情報を悪性の評価値に変換し、複数の前記弾性情報の評価値をバーチャートにより対比させて表示する画像であることを特徴とする弾性画像表示方法。

[13] 請求項 11に記載の弾性画像表示方法にお!、て、

前記一つの弾性画像は、前記関心領域における異なる複数の弾性情報の時間変化を線図により対比させて表示する画像であることを特徴とする弾性画像表示方法。

[14] 被検体から得られた受信信号からなるフレームデータを前記被検体に加わる圧迫が変化する過程で取得するフレームデータ取得手段と、取得時刻が異なる一対の前記フレームデータに基づいて前記走査面の複数の計測点における変位を求める変位計測手段と、前記各計測点の変位に基づいて前記各計測点における異なる複数の弾性情報を求める弾性情報演算手段と、求められた複数の弾性情報に基づいて一つの弾性画像を生成する弾性画像生成手段と、該生成された弾性画像を表示する表示手段とを有してなる弾性画像表示装置。

[15] 請求項 14に記載の弾性画像表示装置において、

前記複数の弾性情報は、弾性率、粘弾性率、歪み量に対する弾性率の非線形性に係る非線形パラメタ、前記変位の演算結果に含まれる局所的なバラツキ又は相関係数等の演算結果の信頼度に係る SZN判定情報、前記各計測点に加わる圧迫方向、変位、応力、歪み、ポアソン比力選択される少なくとも 2つであることを特徴とする弾性画像表示装置。

[16] 請求項 14に記載の弾性画像表示装置において、

前記弾性画像生成手段は、前記弾性情報の 2つを縦軸と横軸とする 2次元マップに従って、各画素の色相又は輝度を前記 2つの弾性情報に応じて割り付けて前記一つの弾性画像を生成することを特徴とする弾性画像表示装置。

[17] 請求項 16に記載の弾性画像表示装置において、

前記 2次元マップは、 2次元カラーマップであり、前記縦軸の左右軸と前記横軸の上下軸は、それぞれ異なる色調のグラデーションが設定されて、ることを特徴とする弾性画像表示装置。

[18] 請求項 16に記載の弾性画像表示装置において、

前記縦軸の弾性情報は弾性率であり、前記横軸の弾性情報は前記歪み量に対する弾性率の非線形性に係る非線形パラメタであることを特徴とする弾性画像表示装置。

[19] 請求項 16に記載の弾性画像表示装置において、

前記縦軸の弾性情報は弾性率であり、前記横軸の弾性情報は粘弾性率であることを特徴とする弾性画像表示装置。

[20] 請求項 16に記載の弾性画像表示装置において、

前記縦軸の弾性情報は弾性率であり、前記横軸の弾性情報は前記各計測点に加わる圧迫方向であることを特徴とする弾性画像表示装置。

[21] 請求項 14に記載の弾性画像表示装置において、

前記弾性画像生成手段は、前記各計測点を隣接する複数の計測点を有する複数の設定領域に分けて、該各設定領域における 2つの前記弾性情報の代表値を求め、 2つの前記弾性情報に対してそれぞれ異なる 2つの色相を割り付けてなる 2つの 1 次元カラーマップに従って、前記各設定領域の各画素に色付けして前記一つの弹性画像を生成することを特徴とする弾性画像表示装置。

[22] 請求項 21に記載の弾性画像表示装置にぉ、て、

前記 2つの前記弾性情報は、弾性率と粘弾性率であることを特徴とする弾性画像表示装置。

[23] 請求項 14に記載の弾性画像表示装置において、

前記弾性画像生成手段は、異なる複数の前記弾性情報を悪性を評価する評価値に変換し、各画素の色相又は輝度を各弾性情報の前記評価値に応じて割り付けて前記一つの弾性画像を生成することを特徴とする弾性画像表示装置。

[24] 請求項 14に記載の弾性画像表示装置において、

前記弾性画像生成手段は、設定される同一の関心領域における異なる複数の前記弾性情報を、数値と線図と図形の少なくとも一つにより対比して表した前記一つの弾性画像を生成することを特徴とする弾性画像表示装置。

[25] 請求項 24に記載の弾性画像表示装置において、

前記一つの弾性画像は、前記関心領域における異なる複数の前記弾性情報を悪性の評価値に変換し、複数の前記弾性情報の評価値をバーチャートにより対比させて表示する画像であることを特徴とする弾性画像表示装置。

[26] 請求項 24に記載の弾性画像表示装置において、

前記一つの弾性画像は、前記関心領域における異なる複数の弾性情報の時間変化を線図により対比させて表示する画像であることを特徴とする弾性画像表示装置。