WO2006121031A1

WO2006121031A1 - 超音波診断装置及び超音波画像表示方法

Info

Publication number: WO2006121031A1
Application number: PCT/JP2006/309290
Authority: WO
Inventors: Takeshi Matsumura; Toshihiko Kawano
Original assignee: Hitachi Medical Corporation
Priority date: 2005-05-09
Filing date: 2006-05-09
Publication date: 2006-11-16
Also published as: EP1880679A1; JPWO2006121031A1; US9060737B2; EP1880679B1; CN101160099A; EP1880679A4; JP5160227B2; CN100544678C; US20090216123A1

Abstract

　圧迫量の影響を排除した組織の硬さの空間分布を画像として表示することを目的とする。そのために、被検体の組織に圧力を加えて計測された超音波断層データに基づいて、被検体の組織の複数の計測点における組織の歪みに相関する物理量を求め、該物理量に基づいて組織の弾性画像を生成し、弾性画像に設定される基準領域の物理量を基準として、各計測点の物理量を指標化し、その指標値の分布を表す指標化弾性画像を生成する。

Description

明細書

超音波診断装置及び超音波画像表示方法

技術分野

[0001] 本発明は、超音波画像表示方法及び超音波診断装置に係り、具体的には、被検体に加えられる圧迫量に対する生体各部の糸且織の歪みに相関する物理量を指標化した画像をユーザに提供して、診断における組織の鑑別性を向上させることができる超音波画像表示方法及び超音波診断装置に関する。

背景技術

[0002] 超音波診断装置は、被検体に用手法又は機械的な方法により被検体に圧力を加えて生体各部の変位を求め、求めた変位に基づいて組織の硬さに関する情報を画像ィ匕することにより、ガン腫瘍等の病変組織を適切に鑑別可能にする装置である。

[0003] 例えば、組織の硬さ情報の一つである歪み画像は、変位分布を空間微分して生体各部の歪みを求め、歪みの大きさに応じて色相や輝度により階調を付けた画像である。しかし、歪みは、圧迫量に依存した定性的な物理量であり、圧迫力の加え方によつて歪みの大きさが変わってしまう。そのため、同一の糸且織であっても、圧迫量が大きいだけで、歪みが大きな値の色相や輝度として表示されることがあり、検者の経験や熟練度によって鑑別に誤認が生ずるおそれがある。

[0004] 一方、特許文献 1には、歪み画像に基づ!/、て圧迫量に相関しな、組織の硬さ特性を取得する試みが提案されている。すなわち、歪み画像内に病変組織とその周辺組織とにそれぞれ関心領域 ROI— l、 ROI— 2を設定し、各 ROI内の歪み ε 1、 ε 2の比 ε \/ ε 2を硬さ指標とすることが示唆されている。これにより、 2つの関心領域内の歪みの相対比は圧迫量に依存しないから、病変組織と周辺組織の歪みの違いを半ば定量的に識別できるという利点がある。

特許文献 1 :米国特許公開 WO 2006/013916号公報

発明の開示

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

[0005] ところで、ガン腫瘍等の病変部の超音波画像を観察して病変組織の良性や悪性を適切に鑑別する場合、硬さの違いによって病変組織と正常組織の境界を識別したり、病変組織の周辺組織の硬さ、あるいは病変組織と周辺組織との境界組織の硬さの変化が判れば、病変部位の組織を適切に鑑別できることが期待できる。

[0006] しかし、特許文献 1に開示された技術では、関心領域 ROIが設定された領域のみの圧迫量に相関しな、組織の硬さ特性を取得することが可能であって、関心領域 ROI が設定されな、領域にっ、ても圧迫量に相関しな、組織の硬さ特性を取得することは考慮されてない。また、関心領域 ROI同士の比では、歪みの 2次元分布情報を充分に抽出して、るわけではな、ので、さらに改善すべき余地が残されて、る。

[0007] 本発明は、圧迫量の影響を排除した組織の硬さの空間分布を画像として表示することを課題とする。

上記課題を解決するため、本発明の超音波診断装置は以下の様に構成される。即ち、超音波探触子と、被検体の組織に圧力を加えて断層部位の超音波断層データを計測する超音波断層データ計測手段と、超音波断層データから断層像を生成する断層像生成手段と、超音波断層データに基づいて断層部位の複数の計測点における組織の歪みに相関する物理量を求め、この物理量に基づいて断層部位における弾性画像を生成する弾性画像生成手段と、を備え、さらに、断層像又は弾性画像にぉ、て基準領域を選択する手段と、基準領域の物理量を基準として各計測点の物理量を指標値化し、指標値の分布を表す指標化弾性画像を生成する手段と、を備える。

[0008] また、本発明の超音波画像表示方法は以下の様に構成される。即ち、被検体の組織に圧力を加えて超音波断層データを計測するステップと、超音波断層データに基づいて被検体の断層部位の複数の計測点における組織の歪みに相関する物理量を求めるステップと、物理量に基づいて断層部位の弾性画像を生成するステップと、弹性画像にお 1、て基準領域を選択するステップと、基準領域の前記物理量を基準として、各計測点の前記物理量を指標値化するステップと、指標値の分布を表す指標化弾性画像を生成するステップと、を有することを特徴とする。

上記本発明の超音波診断装置及び超音波画像表示方法によれば、圧迫量の影響を排除した組織の硬さの空間分布を画像として表示することができる。図面の簡単な説明

[図 1]本発明の超音波画像表示方法を適用可能な超音波診断装置の一実施形態のブロック構成図である。

[図 2]本発明の超音波画像表示方法の一実施形態の処理手順を示すフローチャートである。

[図 3]超音波の探触子の一例を示す外観図である。

[図 4]本発明の実施例の効果を説明するための従来例の歪み画像の問題点を説明する図である。

[図 5]本発明の一実施例の指標化弾性画像を示す図である。

[図 6]本発明の他の実施例の指標化弾性画像を示す図である。

[図 7]本発明の指標値に係る基準領域の自動設定法の一例を説明する図である。

[図 8]本発明の指標値に係る基準領域の自動設定法の他の例を説明する図である。

[図 9]本発明の他の実施例の指標化弾性画像を示す図である。

[図 10A]本発明の指標値に係る基準領域を、圧迫に伴う組織の変異に追従させて移動する方法を説明する図である。

[図 10B]本発明の指標値に係る基準領域 Rの構成する点の追従処理の具体的な実

0

施例を示す図である。

[図 10C]本発明の指標値に係る基準領域 Rの追従処理の具体的な実施例を示す図

0

である。

[図 11]本発明の指標化弾性画像の階調化の例を説明する図である。

[図 12]本発明の指標化弾性画像の階調化の他の例を説明する図である。

[図 13]本発明の指標化弾性画像の階調化のさらに他の例を説明する図である。

[図 14]本発明の指標化弾性画像を Bモード像に重ねて表示する例を説明する図である。

[図 15]本発明の指標化弾性画像に関心領域を設定し、その関心領域の硬さの指標値の平均値を数値表示する実施例を説明する図である。

[図 16]ラテラル方向に弾性画像を細分化して、細分化された画像毎に指標化を行うことによって全体として指標化弾性画像を生成する実施例を説明する図である。発明を実施するための最良の形態

[0010] 以下、本発明を実施形態に基づいて説明する。図 1に、本発明の超音波画像表示方法を実施するのに好適な一実施形態の超音波診断装置のブロック構成図を示す

。図 2に、本発明の特徴部に係る超音波画像表示方法の一実施形態のフローチヤ一卜を示す。

[0011] 図 1に示すように、被検体 1に当接して用いられる超音波の探触子 2は、図 3(A)に示すように、被検体 1との間で超音波を送信及び受信する複数の振動子が整列された超音波送受信面 21を有して形成されている。探触子 2は、送信回路 3から供給される超音波パルスにより駆動される。送受信制御回路 4は、探触子 2の複数の振動子を駆動する超音波パルスの送信タイミングを制御して、被検体 1内に設定される焦点に向けて超音波ビームを形成するようになっている。また、送受信制御回路 4は、探触子 2 の振動子の配列方向に電子的に超音波ビームを走査するようになって!/、る。

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

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

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

[0015] 歪み Z弾性率演算部 13は、入力される変位フレームデータに基づいて断層面における複数の計測点の歪みを求め、弾性フレームデータとして弾性データ処理部 14に出力するようになっている。また、歪み Z弾性率演算部 13は、圧力計側部 19から被検体に加えられた圧力計測データを取り込み、被検体各部の応力分布を求め、先に求めた歪みフレームデータと応力分布とから弾性率を求め、弾性フレームデータとして、弾性データ処理部 14に出力するようになっている。

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

[0017] カラースキャンコンバータ 15は、弾性データ処理部 14から出力される弾性フレームデータを変換してカラーの弾性画像を生成し、切替加算部 9を介して画像表示器 10 に表示させるようになつている。つまり、カラースキャンコンバータ 15は、予め設定された弾性 (歪み又は弾性率)の上限値及び下限値の範囲に基づいて、弾性画像に階調ィ匕 (例えば、 256階調)された赤、緑、青などの色相コードを付与する。例えば、弹性フレームデータの歪みが大きく計測された領域は赤色コードに変換し、逆に歪みが小さく計測された領域は青色コードに変換する。なお、カラースキャンコンバータ 15に代えて、白黒スキャンコンバータを用いることができる。この場合は、歪みが大きく計測された領域は輝度を明るぐ逆に歪みが小さく計測された領域は輝度を暗くするなどにより、歪みの分布を表すことができる。

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

[0019] 次に、本発明の特徴部の実施形態である指標値演算部 16について説明する。指標値演算部 16は、装置制御インターフェイス部 17から出力される制御指令に従って処理を開始し、歪み Z弾性率演算部 13により演算された弾性フレームデータを取り込み、本発明に係る指標化弾性画像フレームデータを生成して、カラースキャンコンバータ 15に出力する。カラースキャンコンバータ 15は、前述したと同様に、指標化弹性画像フレームデータに基づいて階調化した指標化弾性画像を生成し、切替加算部 9を介して画像表示器 10に表示させるようになつている。指標値演算部 16の詳細については後述する。

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

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

[0022] 変位計測部 12で求められた変位フレームデータは、歪み Z弾性率演算部 13に入力され、各計測点の歪み又は弾性率を演算して弾性フレームデータとして弾性データ処理部 14に出力される。歪みの演算は、公知のように変位分布を空間微分することによって計算される。また、求めた歪みに基づいて各計測点の弾性率を演算する。弾性率を求める場合は、圧力計測部 19により計測された圧力の計測値を取り込み、これに基づいて各計測点における応力を演算する。圧力計測部 19は、例えば、図 3 の (B)に示すように、圧迫板 31の表面に参照変形体 33を設けて構成される。歪み Z弹性率演算部 13は、参照変形体 33の変形を基準にして被検体 1内部の計測点における応力を演算する。このようにして、歪み Z弾性率演算部 13は、各計測点における応力と、歪み Z弾性率演算部 13で求めた歪みフレームデータから断層面上の各計測点の弾性率 (例えば、ヤング率 Ym)を演算し、弾性データ処理部 14に出力する。弾性データ処理部 14は、歪み又は弾性率に基づいて弾性画像データを生成し、力ラースキャンコンバータ 15と切替加算部 9を介して画像表示器 10に弾性画像を表示する。

[0023] 次に、図 2を参照して、本実施形態の特徴部である指標値演算部 16に係る詳細構成を、実施例に分けて動作とともに説明する。

実施例 1

[0024] 図 2に、実施例 1の指標化弾性画像の生成手順をフローチャートを示す。本実施例は、図 4に示すように、弾性画像の一つである歪み画像 41に基準領域 Rを設定し、こ

0

れに基づいて図 5に示す指標化弾性画像 51を生成して表示する例である。

[0025] まず、図 2のステップ S1において、歪み Z弾性率演算部 13に指令を送り、歪み画像 41を画像表示器 10に表示させる。次いで、装置制御インターフェイス部 17を介してマウス等を用いて手動ある!、は自動で設定された基準領域 Rの

0 座標データを取り込む

(S2)。そして、歪み Z弾性率演算部 13力歪み画像フレームデータを取り込み、基準領域 Rに対応する領域の歪みの平均値を基準歪み ε として演算する (S3)。次いで、各歪み計測点 (i, j)における歪み ε と基準歪み ε との相対比である指標値 R を次

1, J i, J 式により求める (S)。求めた R 力もなる指標化弾性画像データは、カラースキャンコン

1，

バータ 15に出力される (S5)。これにより、画像表示器 10には、各計測点の歪みが基準歪み ε により正規化された指標化弾性画像が表示される。

[0026] R = ε / ε (1)

i, j i, j

なお、比に限らず、 ε と ε の違いが反映された指標であれば良い。例えば、

i, j

R =log( ε ) -1ο_§( ε )

でも良い。或いは、

R = ( ε - ε )/( ε + ε )

1,

又は

R = ( ε - ε )/ ε

i, j i, j

とすると、硬、領域の計測点では ε 〜0となり、軟らか、領域の計測点では ε 〜 ε

1,

となるので、

硬い計測点: R →1，軟らかい計測点: R →0

のように 0〜1の間の値を取るように規格ィ匕された指標とすることができる。

[0027] 本実施例の効果について、図 4及び図 5を参照して具体的に説明する。まず、生体糸且織を自由に圧迫してその歪みを計測する場合、歪みの大きさはそのときの圧迫量に依存する。例えば、図 4(A)に示すように、あるときの圧迫量では生体内の計測点 1と計測点 2の歪みが 0.8%、計測点 3の歪みが 0.1 %として計測されたとする。これに対し、別なときの圧迫量では、図 4(B)に示すように、計測点 1と計測点 2における歪みが 1.6 %となり、計測点 3における歪みが 0.2%として計測されることがある。このような歪み画像を、歪みの大きさ (％)に基づいて階調化すると、図 4(B)ように、計測点 1、 2の歪み力 S l .6%のときの圧迫量では、関心部位の腫瘍部 42に縁取りが生じるような画像パタンが得られる。し力し、図 4(A)のように、計測点 1、 2の歪みが 0.8%のときの圧迫量では、階調レベルの範囲が感度の悪い範囲に変わるので、関心部位の腫瘍部 42の縁取りの画像パタンを鮮明に表示することができなくなる。つまり、歪み画像に基づいた画像診断では、同一糸且織であっても、その時の圧迫量に応じて、色相や輝度の強度及び病変部位のパタンが変化してしまうから、病変部位の境界を適切に認識できなくなる。なお、図 4において、 (A), (B)の右側に示した図は、それぞれ計測点 1, 2を結んだ解析ライン上における歪み分布を表して、る。

[0028] これに対し、本実施例によれば、図 5(A)又は図 5(B)に示すように、基準領域 Rを設

0 定し、その領域の歪みの平均値を基準歪み ε として、各計測点の歪みを正規化し

0 i, j

た歪みの指標値 R に基づ!/、て階調化して指標化弾性画像 51を生成して、ることか，

ら、圧迫量に依存しない組織固有の硬さの情報を取得することができる。その結果、圧迫量が異なることにより、同一の糸且織 1に設定された基準領域 R

0の基準歪み ε

0が

、 8%又は 16%のように異なっても、図 5(Α)、（Β)を対比して明らかなように、輝度分布などの階調分布が変わらない指標化弾性画像 51が得られる。したがって、関心部位の腫瘍部 42の縁取りの画像パタンを感度良く認識できる。特に、腫瘍部 42の周辺組織 2との境界組織 44の硬さの違いを識別することができる。また、 05(A), (Β)の右側に示した図は、計測点 1, 2を結んだ解析ライン上における指標値の分布を表している。

[0029] 例えば、乳癌検診における計測断面内においては、脂肪が最も柔らかい組織であり、この領域に基準領域 Rを設定すれば、それ以外の領域は脂肪より硬いため、歪

0

みの指標値 R である相対比は 1以上の値をもち、硬い組織ほどその相対比の値は大，

きく計測されることとなる。或いは、大胸筋に基準領域 R

0を設定しても良い。大胸筋は肋骨に支えられて固定されているために、感度良く基準歪みを検出することができる。また脂肪の厚さは個人差が大きく充分な広さの領域を基準領域として確保できないことがあるが、大胸筋は個人差が少なく必ず計測面に入れることができる。

また、腹部においては、例えば腎臓に基準領域を設定し、数式 (1)に基づいて肝臓の各計測点 G, j)の歪みとの比 R を求めて、肝臓の指標化弾性画像を生成しても良い。

また、前立腺部位を撮影する際には、前立腺皮膜部の脂肪組織を基準領域としても良い。

[0030] 本実施例においては、歪み画像に基づいて指標化弾性画像を生成する例を説明したが、本発明はこれに限らず、組織の硬さの特性に関係した情報である弾性率、粘性などの組織の特性データを用いて指標化弾性画像を生成することができ、本実施例と同様の効果を得ることができる。 [0031] また、本実施例では、基準領域 Rにおける歪みの平均値を基準歪み ε として指標

0 0

値 R を求めたが、本発明はこれに限らず、基準領域 Rにおける歪みの中央値や、ヒ i, j 0

ストグラム処理などで最も頻度が高く計上された歪みの値などを基準歪み ε として、

0 各計測点の歪み ε を正規ィ匕することができる。

，

また、基準領域 Rを、複数の計測点を含む一定範囲の領域とするのに代えて、基

0

準部位として 1点の計測点の歪みの値を基準歪み ε として設定してもよい。

0

[0032] 以上のように、本発明の指標化弾性画像は、設定された基準部位の歪みに相関する物理量を基準として、断層部位の複数の計測点における歪みに相関する物理量が正規化した指標値の画像として表される。その指標値は、同一の圧迫量を受けた各部位の相対値であるから、圧迫量が変化しても指標値は殆ど変わらない。したがって、本発明の指標化弾性画像は、圧迫量の影響を排除して、病変組織から周辺組織に至る硬さの変化、あるいは境界領域の硬さの違いを、基準部位の硬さを基準に半ば定量的に表示できる。その結果、病変組織の良性や悪性を客観的かつ普遍的に適切に鑑別することが容易になる。

実施例 2

[0033] 実施例 1においては、指標化の基準部位である基準領域 Rを生体組織内に設定す

0

る例を示した。しかし、本発明はこれに限らず、図 3(B)に示した弾性部材である参照変形体 33を有する探触子 2を用いて超音波断層データを取得すると、図 6に示すように、歪み画像では最も探触子 2に近い組織領域に、参照変形体 33の画像 53が表示される。この参照変形体 33の画像領域に基準領域 Rを設定して、実施例 1と同様に指

0

標化弾性画像 52を生成して表示することができる。

[0034] 特に、参照変形体 33の弾性率を予め計測しておくことにより、各計測点における組織のおおまかな弾性率を推定することができる。つまり、参照変形体 33の弾性率を Er e kPa)とすると、次式 (2)により、各計測点の大まかな弾性率 E を求めることができる，

。これを画像化することにより、より精細な硬さの画像情報を取得することができる。

E =ErefX R (kPa) (2)

i, j i, j

[0035] 上記説明にお!/、ては、体外の参照変形体を利用する例を示した力本発明はこれに限らず、被検体内に参照変形体を配置してもよい。例えば、前立腺部位を撮影する場合、足の付け根力も挿入するカテーテルを参照変形体とする。このカテーテルは

、中空ではなくゴム等で満たされて構成されている。カテーテルの外径は、 l〜2mm である。

前立腺部位を撮影する際、前立腺用探触子を用いて撮影を行う。前立腺の近傍にカテーテルを挿入した血管があり、得られる画像にはカテーテルの断面が表示される

。このカテーテル部位を参照変形体とする。

[0036] また、中空のカテーテルの先端に伸縮するノレーンを設け、体内に挿入されたバルーンを参照変形体としてもよい。このバルーンは中空であり、中に液体ゴムゃゼリ一等の液体を挿入することにより、伸縮可能となっている。膨張するときの外径は、 0. 5〜5mmであり、液体の挿入量の調整により、任意に外径を設定することができる。実施例 3

[0037] 実施例 1においては、基準領域 Rを装置制御インターフェイス 17力入力設定する

0

例を説明したが、次に説明するように、基準領域 R

0を自動設定することができる。例えば、歪み画像などにおいて最も軟らかい領域を検出し、その領域に基準領域 R

0を設定することができる。

[0038] また、例えば、図 7(A)に示す図は、乳腺部位の Bモード断層像 54であり、探触子 2に最も近、位置に脂肪組織が表れ、その脂肪組織が計測部位における最も軟らか!/ヽ組織でもある。ここで、 Bモード画像 54でも脂肪 61と乳腺 62の境界を検出することができる力、体表近傍力乳腺との境界までの領域を抽出し、基準領域 R

0として自動設定することが可能である。

この場合は、図 7(B)に示すように、脂肪 61の組織 1の全体を基準領域 Rとする指標

0 化弾性画像 55が得られる。

[0039] また、例えば、図 8(A)に示すように、マウスなどの入力デバイスを操作して表示画面上のカーソル 63を所望の基準領域に移動してクリックして指定することにより、図 8(B) に示すように、歪み値などの弾性情報や、 Bモード断層像の輝度情報やパタン情報に基づいてその領域を含む同一組織領域を同定して、自動で基準領域 Rを設定す

0 ることがでさる。

このように基準領域 Rを自動設定することにより、基準領域 Rの設定にあたっての任意性を排除できるから、客観的な硬さの指標化弾性画像による診断が可能となる。

[0040] なお、上記の説明では、乳腺を例に示したが、本発明はこれに限らず、図 9に示す指標化弾性画像 56のように、頸動脈部位 57におけるプラーク 58の弾性評価などにおいても同様に適用できる。この場合は、例えば、筋 59の領域を基準領域 Rに設定し、

0

血管壁 60やプラーク 58の歪みの指標値を画像ィ匕して診断することができる。

実施例 4

[0041] 上記の各実施例において設定した基準領域 Rは、圧迫過程において組織が移動

0

すると、設定した同一組織の領域力基準領域 R

0が外れる場合がある。つまり、図 10

A(a)〜(c)に示すように、基準領域 Rを常に一定の座標領域として固定して設定して

0

いる場合、圧迫過程において、別の組織が基準領域 R

0の領域に入ってくることがあるこの場合は、同図 (d)〜(i)に示すように、常に同一の組織領域を追従して、基準領域 R

0を変更設定することが好ましい。このようにした場合、圧迫過程において別の組織が基準領域 Rに入ってくるようなことがなぐ常に同一組織内において計測された

0

歪みの値を基準歪み _ε とすることができ、外乱を排除して指標値を精度良く求める

0

ことができる力、診断の信頼性が向上する。

[0042] 上記基準領域 Rの追従処理の具体的な実施例を以下に示す。

0

[図 1]における変位計測部 12において、変位分布 (変位フレームデータ）を求める処理を具体的に [図 10Β]を用、て示す。

例えば、 Βモード画像上の横方向に y座標、深度方向に X座標を設定した座標系を設け、ライン ylに沿った方向にある計測点

(xl, yl)、 (x2, yl)

の 2点に注目する。図示したように、時刻 t-l(過去)から時刻 t (現在)の時間変化の間に、対象組織が圧迫されたものとする。

時刻 t- 1において、座標 (xl(t- 1), yl(t- 1》にあった組織は、時刻 tにおいて座標 (xl(t ), yl(t》に移動し、同様に、時刻 t- 1において、座標 (x2(t- 1), yl(t- 1》にあった組織は、時刻 t〖こおいて座標 (x2(t), yl(t》に移動する。

このとき、変位計測部 12において、全計測点 (X, y)における変位 d(x, y)を演算しており、例えば、上記時刻 t-1から tにおけるライン ylに沿った変位の分布 d( x(t-l) )は、図に示すように求められ、現時刻 tにおける上記移動後の X座標、 xl(t)、 x2(t)は、それぞれ、

xl(t)=xl(t- l) + d(xl(t- 1》

x2(t)=x2(t-l) + d(x2(t-l))

として求めることができる。

同様にして、 y座標方向の移動後の座標 yl(t), y2(t)も求められる。

[0043] 以上の方法を基づけば、例えば [図 10C]に示すように、時刻 t-1において

(xl(t-l), yl(t- 1》、 (x2(t-l), yl(t- 1》、 (xl(t- 1), y2(t- 1》、 (x2(t- 1), y2(t- 1)) の 4点力もなる長方形の基準領域 Rを設定したとき、上記 4点における組織が時刻 tに

0

おいてどこに移動したかを求めることができ、これにより、 Rの内部における組織は同

0

一の組織領域が指定されて、ることになる。

上記の処理をリアルタイムに行うことができ、 [図 10A]に示すように、基準領域 R

0として、同一組織領域をトラッキングして追従するようにすることができる。

上記の説明では、基準領域 Rを指定する角に配置された 4点の座標における変位

0

の情報を用いて追従する処理を示したが、その方法に限らず、基準領域 Rの

0 境界を示す線上に沿って配置された座標における変位の情報を用いて、より多くの情報に基づくことにより高精度に追従処理を行うようにしてもよい。

また、基準領域 Rの内部の座標における変位の情報を用いて、更に高精度に追従

0

処理を行うようにしてもよい。

実施例 5

[0044] 上記各実施例においては、各計測点の指標値に基づいて指標化弾性画像を階調化する際の指標値と階調レベルとの対応については特に説明しな力つたが、指標値と階調レベルとの関係を規定する階調化関数は、比例関数、対数関数、及び組織の良性又は悪性を判断する閾値近傍の分解能を高くした設定関数等、種々のものを適用することができる。

具体的には、図 11(A)のバー表示 71に示すように、階調化関数を比例関数とし、指標値と色相や輝度で表現された例えば 8ビットの 256段階の階調マップとをリニアの関係に割り当てることができる。

[0045] また、図 11(B)のバー表示 72に示すように、階調化関数を対数関数で規定することができる。特に、色相、輝度が敏感に変化する範囲 73に、閾値 Rthが設定される特定の関数を定義することも可能である。例えば、良性と悪性を効果的に識別する指標値の閾値 Rthが得られて、る場合、この閾値 Rthの近傍にぉ、て色相や輝度が敏感に変化する階調化関数を定義することにより、高感度に良悪性を画像識別できるようにすることができる。

[0046] また、図 11(C)のバー表示 80に示すように、弾性率の指標としてある特定部位を表示させてもよい。例えば、正常部位である脂肪は 20kPaであり、それを指標として表示させる。他の部位も弾性率に換算して表示させる。したがって、正常部位に対して高感度に良悪性を画像で認識できるようにすることができる。

[0047] さらに、図 12(A)のバー表示 74の範囲 75に示すように、閾値 Rthより大きい場合と小さい場合とで色相や輝度が大きく変化する特別な階調マップの階調化関数を準備し、それによりマッピングするようにすることができる。また、同図 (B)に示すように、バー表示 76の指標値力 S小さい範囲 77に、階調マップの色相や輝度が段階的に不連続的に変化するようにすることができる。

[0048] また、図 13に示すように、バー表示 78の例えば閾値 Rthの近傍を、ある一定の色づけ範囲 79を設定し、その色づけ範囲 79をポインティングデバイス等のより検者が上下させて閾値 Rthを変えることにより、関心領域の指標値の大きさがどの範囲に入っているのかが簡単に解析できるようにすることができる。また、同様に、検者が閾値 Rth をリアルタイムに上下して、その閾値 Rthを超える領域にだけ色をつけるなどの処理を行うようにすることができる。

その結果、指標化弾性画像について関心のある病変部位に、特に関心領域を設定し、その関心領域の指標値を求めて表示させることなぐ階調マップの閾値 Rthの位置を自由に移動させることにより、関心領域の組織の硬さの指標値を正確かつ容易に把握することができる。

実施例 6

[0049] 上記の各実施例で形成される指標化弾性画像は、 Bモード画像、歪み画像などの弾性画像に並べて、あるいは独立に表示することができる。さらに、例えば、図 14に示すように、白黒 Bモード画像に指標化弾性画像をカラーの半透明画像として重畳して表示するようにしてもょヽ。

また、例えば、 Bモード画像と指標化弾性画像を 2画面で独立に表示する場合は、 B モード画像をガイドとして、指標化弾性画像の基準領域 Rのサイズや位置などを設

0

定可能にすることができる。また、 Bモード画像に限らず、指標化弾性画像とペアになつて表示される画像をガイドとして、基準領域 R

0を設定するようにしてもよい。

実施例 7

[0050] さらに、本発明は、表示された指標化弾性画像において、関心領域 ROIを設定して、基準領域 R

0との組織の硬さの解析を行わせることができる。例えば、図 15に示すように、指標化弾性画像上に関心領域 ROIを設定することにより、 ROI内の歪みの平均値などを求めて、例えば、その関心領域 ROIの近傍に対応付けて表示する。また、このような関心領域 ROIは、複数設定できる。

実施例 8

[0051] 上記の各実施例は、 Bモード像あるいは歪み画像とリアルタイムで指標化弾性画像を生成して表示することを説明したが、本発明はこれに限らず、フリーズ後にオフラインで諸設定を行い、シネメモリに確保された弾性フレームデータを利用して、指標化弾性画像を生成して表示するようにしてもょヽ。

実施例 9

[0052] 上記の各実施例は、探触子 2により被検体 1の生体組織が均一に圧迫されることを想定しているが、被検体 1の体表面に対して探触子 2が傾いた状態で圧迫される等の圧迫操作が不適切である場合には、生体組織が不均一に圧迫されることになる。特に送受信方向 (つまり音線方向)に垂直な方向 (以下、ラテラル方向という)に圧迫が不均一となってしまう。その結果、基準領域の歪みがラテラル方向に不均一となり指標値が不正確となってしまう可能性がある。そこで本実施例では、基準領域の歪みが不均一となる場合でも、その影響を排除して精度良く指標値を求める。そのために、画像領域全体をラテラル方向に細分化して、細分化された各領域内では均一に圧迫されているとみなすことができるようにする。そして、細分化された領域毎に指標値を算出する。

[0053] 図 16に本実施例を具体的に示す。図 16は、画像表示器 10に表示された指標化弹性画像であって、基準領域とそれ以外の領域を含む指標化弾性画像を表して!/ヽる例である。この指標化弾性画像をラテラル方向に短冊状に例えば M個に細分ィ匕して、音線方向に平行な細長い複数の細分領域 R , R , · ··, R

1 2 k , を得る。それぞれの細分領域内における基準領域を細分基準領域 R , R , · ··, R ,〜R

01 02 Ok 0として得る。各細分基準領域のラテラル方向の幅は、例えば数 mm程度とすることができる。そして、各計測点の歪みの指標値を、その計測点の細分領域における細分基準領域の基準歪み ε , ε , · ··, ε ,… ε を用いて求める。つまり、例えば、図 16に示すよう

01 02 Ok 0

に、計測点 G, j)にお、て計測された歪み ε は、その細分領域内にある細分基準領域 R を介して与えられたと考えられる。従い、計測点 (i, j)における指標値 R を次式（

Ok i, j

3)により求める。

R = ε / ε (3)

O k

[0054] これにより、圧迫操作が不適切で基準領域の歪みがラテラル方向に不均一となっても、細分ィ匕された各領域内では均一に圧迫されているとみなすことができるので、精度良く指標値を細分化された領域毎に求めることができる様になる。その結果、指標化弾性画像全体として精度良く指標値を求める事が出来るようになる。

[0055] 上記説明においては、細分領域をラテラル方向に比較的広くとるように説明したが、細分領域を各ライン毎に独立して設定し、計測点 G, j)における指標値 R を次式 (4) ， j により求めるようになって!/、てもよ!、。

R = ε / ε (4)

0

Claims

請求の範囲

[1] 超音波探触子と、被検体の組織に圧力を加えて該組織の超音波断層データを計測する超音波断層データ計測手段と、前記超音波断層データから断層像を生成する断層像生成手段と、前記超音波断層データに基づいて複数の計測点における組織の歪みに相関する物理量を求め、該物理量に基づいて前記組織の弾性画像を生成する弾性画像生成手段と、を備えた超音波診断装置において、

前記断層像又は前記弾性画像において基準領域を選択する手段と、

前記基準領域の前記物理量を基準として前記各計測点の前記物理量を指標値ィ匕し、前記指標値の分布を表す指標化弾性画像を生成する手段と、

を備えたことを特徴とする超音波診断装置。

[2] 請求項 1記載の超音波診断装置において、

前記指標化弾性画像生成手段は、前記指標値として、前記基準領域の物理量と前記計測点の物理量との比を求めることを特徴とする超音波診断装置。

[3] 請求項 1記載の超音波診断装置において、

前記指標化弾性画像生成手段は、前記指標値として、前記基準領域の物理量と前記計測点の物理量とを用いて、固さを表す規格化された値を求めることを特徴とする超音波画像表示方法。

[4] 請求項 1記載の超音波診断装置において、

前記指標化弾性画像生成手段は、前記指標値として、前記基準領域の物理量の対数と前記計測点の物理量の対数との差を求めることを特徴とする超音波診断装置

[5] 請求項 1記載の超音波診断装置において、

前記基準領域選択手段は、前記断層像又は前記弾性画像上で選択された基準点に基づいて前記基準領域を選択することを特徴とする超音波診断装置。

[6] 請求項 5記載の超音波診断装置において、

前記基準領域選択手段は、前記基準領域として、前記弾性画像上の各点の前記物理量に基づヽて、前記基準点が設定された組織と同一組織を有する領域を選択することを特徴とする超音波診断装置。

[7] 請求項 5又は 6に記載の超音波診断装置において、

前記基準領域選択手段は、前記被検体に加えられる圧力に応じて変位する前記組織の画像上の移動先に合わせて前記基準領域を移動させることを特徴とする超音波診断装置。

[8] 請求項 1記載の超音波診断装置において、

前記指標値と階調との関係を規定する階調化関数に基づ!、て、前記指標化弾性画像を階調化する階調化手段と、

前記階調化された指標化弾性画像を表示する表示手段と、

を備えたことを特徴とする超音波診断装置。

[9] 請求項 8記載の超音波診断装置において、

前記階調化関数の前記指標値と前記階調との関係をバー表示することを特徴とする超音波診断装置。

[10] 請求項 9記載の超音波診断装置において、

前記表示手段は、白黒の前記断層像と前記階調化された指標化弾性画像とを重ね合わせて表示することを特徴とする超音波診断装置。

[11] 請求項 1記載の超音波診断装置において、

前記表示手段は、前記指標化弾性画像に設定された関心領域の前記物理量の平均値と、前記基準領域の前記物理量との比を数値により前記関心領域に対応付けて表示することを特徴とする超音波診断装置。

[12] 請求項 1に記載の超音波診断装置において、

前記指標化弾性画像生成手段は、前記超音波探触子の送受信方向に垂直な方向に前記弾性画像を複数の領域に分割し、各領域毎に前記指標値化を行うことを特徴とする超音波診断装置。

[13] 請求項 1に記載の超音波診断装置において、

前記弾性画像生成手段は、前記組織の歪みに相関する物理量として、前記組織の変位、歪み、弾性率及び粘性のいずれか一つを求めることを特徴とする超音波診断装置。

[14] 被検体の組織に圧力を加えて超音波断層データを計測するステップと、前記超音波断層データに基づいて複数の計測点における組織の歪みに相関する物理量を求めるステップと、

前記物理量に基づいて前記組織の弾性画像を生成するステップと、

前記弾性画像において基準領域を選択するステップと、

前記基準領域の前記物理量を基準として、前記各計測点の前記物理量を指標値化するステップと、

前記指標値の分布を表す指標化弾性画像を生成するステップと、

を有することを特徴とする超音波画像表示方法。

[15] 請求項 14記載の超音波画像表示方法にぉ、て、

前記基準領域設定ステップでは、前記基準領域として脂肪組織が選択されることを特徴とする超音波画像表示方法。

[16] 請求項 14記載の超音波画像表示方法にお、て、

前記基準領域設定ステップでは、前記基準領域として筋組織が選択されることを特徴とする超音波画像表示方法。

[17] 請求項 14記載の超音波画像表示方法にぉ、て、

前記基準領域設定ステップでは、前記被検体の外に配置された参照変形体に対応する画像上の領域内に前記基準領域が選択されることを特徴とする超音波画像表示方法。

[18] 請求項 17記載の超音波画像表示方法において、

前記参照変形体は、前記探触子と前記被検体との間に挿入されて用いられる弾性部材であることを特徴とする超音波画像表示方法。

[19] 請求項 14記載の超音波画像表示方法にお、て、

前記基準領域設定ステップでは、前記被検体の内に挿入された参照変形体に対応する画像上の領域内に前記基準領域が選択されることを特徴とする超音波画像表示方法。

[20] 請求項 19記載の超音波画像表示方法にお、て、

前記参照変形体は、前記被検体の内に挿入されたカテーテル又はバルーンであることを特徴とする超音波画像表示方法。