特許 5735718 超音波診断装置、及び弾性評価方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】5735718

(24)【登録日】2015年4月24日

(45)【発行日】2015年6月17日

(54)【発明の名称】超音波診断装置、及び弾性評価方法

(51)【国際特許分類】

   A61B 8/08 20060101AFI20150528BHJP

   A61B 8/14 20060101ALI20150528BHJP

【ＦＩ】

   A61B8/08

   A61B8/14

【請求項の数】13

【全頁数】19

(21)【出願番号】特願2014-547611(P2014-547611)

(86)(22)【出願日】2014年3月27日

(86)【国際出願番号】JP2014058806

(87)【国際公開番号】WO2014162966

(87)【国際公開日】20141009

【審査請求日】2014年10月1日

(31)【優先権主張番号】特願2013-77701(P2013-77701)

(32)【優先日】2013年4月3日

(33)【優先権主張国】JP

【早期審査対象出願】

(73)【特許権者】

【識別番号】390029791

【氏名又は名称】日立アロカメディカル株式会社

(74)【代理人】

【識別番号】110000350

【氏名又は名称】ポレール特許業務法人

(72)【発明者】

【氏名】吉川 秀樹

(72)【発明者】

【氏名】浅見 玲衣

(72)【発明者】

【氏名】田原 麻梨江

【審査官】 後藤 順也

(56)【参考文献】

【文献】 特開２０１２−２４９８５２（ＪＰ，Ａ）

【文献】 特表２０１３−５００７５２（ＪＰ，Ａ）

【文献】 特開２００７−１９５６６２（ＪＰ，Ａ）

【文献】 特開２０１４−０４２６３７（ＪＰ，Ａ）

【文献】 米国特許出願公開第２００８／０２４９４０８（ＵＳ，Ａ１）

【文献】 米国特許出願公開第２０１０／００１６７１８（ＵＳ，Ａ１）

【文献】 新編 画像解析ハンドブック，東京大学出版会，２００４年 ９月１０日，初版，第１２３４−１２３８頁

【文献】 吉川秀樹他，せん断波を利用した組織弾性評価における組織構造物の影響，超音波医学，２０１３年 ４月１５日，第４０巻、増刊号，第Ｓ４２５頁

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６１Ｂ ８／００−８／１５

医中誌ＷＥＢ

ＪＳＴＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）

(57)【特許請求の範囲】

【請求項1】

超音波を送受信する探触子と、
前記探触子から検査対象に超音波を送受信し、前記検査対象から得られる受信データを処理する処理部と、
表示部を備え、
前記処理部は、
前記探触子から第１の超音波を検査対象に送受信して得られた受信データにより、前記検査対象の組織構造を検出し、
検出した前記組織構造に基づいて、せん断波速度を検出する計測領域を決定し、
前記計測領域に第２の超音波を送信して、せん断波を発生させ、
前記計測領域に第３の超音波を送受信して得られた受信データにより、せん断波速度を算出し、
前記組織構造の位置を示す構造マップを作成し、前記構造マップを前記表示部に表示し、
前記計測領域の調整の際、前記構造マップに基づいて、前記組織構造を避ける位置に前記計測領域を限定、又は前記計測領域を移動させる、
ことを特徴とする超音波診断装置。

【請求項2】

請求項１に記載の超音波診断装置であって、
前記処理部は、
前記組織構造を検出する構造検出部と、
前記計測領域を調整する計測領域調整部と、
前記計測領域で計測した前記せん断波速度を用いて前記計測領域の弾性を評価する弾性評価部と、を備える、
ことを特徴とする超音波診断装置。

【請求項3】

請求項２に記載の超音波診断装置であって、
前記構造検出部は、前記構造マップを生成する、
ことを特徴とする超音波診断装置。

【請求項4】

請求項１に記載の超音波診断装置であって、
前記処理部は、
前記構造マップとして、前記第１の超音波を前記検査対象に送受信して得られた前記受信データに基づき生成した輝度モード（Ｂ）画像の輝度を反転した反転画像、前記組織構造に配色をしたカラー画像或いはドプラ画像、組織弾性/歪を示す画像、の少なくとも一つを利用する、
ことを特徴とする超音波診断装置。

【請求項5】

請求項１に記載の超音波診断装置であって、
前記処理部は、
前記構造マップを、前記組織構造の形状に合わせたマスクを有するフィルタ処理により作成する、
ことを特徴とする超音波診断装置。

【請求項6】

請求項１に記載の超音波診断装置であって、
前記処理部は、
前記構造マップを、空間的に異なる位置で取得する複数の前記受信データ、又は空間的に同じ位置で異なる時間で取得する複数の前記受信データを利用して、相関演算により算出する、
ことを特徴とする超音波診断装置。

【請求項7】

請求項１記載の超音波診断装置であって、
前記処理部は、
前記計測領域の正当性を、色の違い、色の濃淡、色の透明度で視覚的に表示する、
ことを特徴とする超音波診断装置。

【請求項8】

請求項２記載の超音波診断装置であって、
前記表示部に前記計測領域の設定を自動、又は手動で行うことを選択する領域設定ボタンを表示可能である、
ことを特徴とする超音波診断装置。

【請求項9】

請求項２に記載の超音波診断装置であって、
前記計測領域調整部は、前記計測領域が弾性評価に対して不適正な場合に、画像又は音としてアラート情報を出力する、
ことを特徴とする超音波診断装置。

【請求項10】

超音波を送受信する探触子から第１の超音波を検査対象に送受信して得られた受信データにより、前記検査対象の組織構造を検出し、検出した前記組織構造に基づいて、せん断波速度を検出する計測領域を決定し、
前記計測領域に第２の超音波を送信して、せん断波を発生させ、
前記計測領域に第３の超音波を送受信して得られた受信データにより、せん断波速度を算出し、前記せん断波速度に基づいて前記検査対象の弾性評価を行い、
前記組織構造の位置を示す構造マップを作成し、前記構造マップを表示部に表示し、
前記計測領域の調整の際、前記構造マップに基づいて、前記組織構造を避ける位置に前記計測領域を限定、又は前記計測領域を移動させる、
ことを特徴とする弾性評価方法。

【請求項11】

請求項１０に記載の弾性評価方法であって、
前記構造マップとして、前記第１の超音波を前記検査対象に送受信して得られた前記受信データに基づき作成したＢ画像の輝度を反転した反転画像、前記組織構造に配色をしたカラー画像、或いはドプラ画像を利用する、
ことを特徴とする弾性評価方法。

【請求項12】

請求項１０に記載の弾性評価方法であって、
前記構造マップを、前記組織構造の形状に合わせたマスクを有するフィルタ処理により作成する、
ことを特徴とする弾性評価方法。

【請求項13】

請求項１０に記載の弾性評価方法であって、
前記構造マップを、空間的に異なる位置で取得する複数の前記受信データ、又は空間的に同じ位置で異なる時間で取得する複数の前記受信データを利用して、相関演算により算出する、
ことを特徴とする弾性評価方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は超音波診断装置に係る。特に、音響放射力を利用した生体の弾性を評価する技術に関する。

【背景技術】

【0002】

超音波やＭＲＩ（Magnetic Resonance Imaging）、Ｘ線ＣＴ（Computed Tomography）に代表される医療用の画像表示装置は、目視できない生体内の情報を数値または画像の形態で提示する装置として広く利用されている。中でも超音波を利用した画像表示装置は、他の装置と比較して高い時間分解能を備えており、拍動下の心臓を滲みなく画像化できる性能を持つ。

【0003】

生体内を伝搬する超音波は主に縦波と横波に区別され、製品に搭載されている多くの技術（組織形態を映像化する技術や血流速度を計測する技術）は、主に縦波（音速約1540m/s）の情報を利用している。

【0004】

近年、横波（以降、せん断波）を利用して組織の弾性率を評価する技術が注目されており、乳腺腫瘍や慢性肝疾患に対して臨床利用が進められている。この技術では、計測対象となる組織内部にせん断波を発生させ、その伝搬速度から弾性を評価する。せん断波を発生させる手法は、機械方式と放射圧方式に大別される。機械方式は、バイブレータ等を利用して体表面に1kHz程度の振動を与えてせん断波を発生させる方式で、振動源となる駆動装置が必要である。一方、放射圧方式は、超音波を組織内の局所に集中させる集束超音波を利用して生体内に音響放射圧を加え、それに伴う組織変位によってせん断波を発生させる。いずれの方式も、発生したせん断波の伝搬に伴う組織変位を超音波を利用して計測し、組織の硬さに関する情報を評価する技術である。

【0005】

これらに関連する先行技術文献としては、例えば、音響放射圧を利用した弾性評価の手法に関する特許文献１がある。

【先行技術文献】

【特許文献】

【0006】

【特許文献1】ＵＳ ８，１１８，７４４ Ｂ２

【発明の概要】

【発明が解決しようとする課題】

【0007】

特許文献１に記載の手法においては、集束超音波を利用して組織内に放射力を発生させ、組織内にせん断波を伝搬させる。そして、伝搬方向には超音波送受信を実施する複数の計測地点を設け、組織変位の時間変化を計測する。そして、変位の計測結果を利用して各計測地点におけるせん断波の到来時間を計測する。そして、各計測地点における到来時間を利用して、各計測地点間におけるせん断波の伝搬時間を算出し、速度を計測する。

【0008】

しかしながら、上述の特許文献１に記載の手法では、せん断波の伝搬経路に血管や線維組織などの組織構造が存在する場合、波面は回折や屈折の影響を受けて散乱し波面形状が乱れる。伝搬に伴う組織の変位量を利用してせん断波速度を推定する本手法においては、波面の乱れは弾性評価の誤差を増大させる大きな要因となる。

【0009】

本発明の目的は、組織構造による波面散乱の影響を軽減し、信頼性が高い弾性評価が可能な超音波診断装置、及び方法を提供することにある。

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

【0010】

上記目的を達成するため、本発明においては、超音波診断装置であって、
超音波を送受信する探触子と、探触子から検査対象に超音波を送受信し、検査対象から得られる受信データを処理する処理部を備え、処理部は、探触子から第１の超音波を検査対象に送受信して得られた受信データにより、検査対象の組織構造を検出し、検出した組織構造に基づいて、せん断波速度を検出する計測領域を決定し、計測領域に第２の超音波を送信して、せん断波を発生させ、
計測領域に第３の超音波を送受信して得られた受信データにより、せん断波速度を算出する超音波診断装置を提供する。

【0011】

また、上記の目的を達成するため、本発明においては、弾性評価方法であって、超音波を送受信する探触子から第１の超音波を検査対象に送受信して得られた受信データにより、検査対象の組織構造を検出し、検出した組織構造に基づいて、せん断波速度を検出する計測領域を決定し、計測領域に第２の超音波を送信して、せん断波を発生させ、計測領域に第３の超音波を送受信して得られた受信データにより、せん断波速度を算出する弾性評価方法を提供する。

【発明の効果】

【0012】

本発明により、組織構造による波面乱れを抑制し、信頼性が高い組織弾性評価が実現される。

【図面の簡単な説明】

【0013】

【図1】実施例１の超音波診断装置の一構成例を示すブロック図である。

【図2】実施例１の超音波診断装置の他の構成例を示すブロック図である。

【図3】実施例１の処理工程の全体を示すフローチャート図である。

【図4】実施例１の弾性評価部の機能を示す図である。

【図5】実施例１の組織構造による波面乱れを説明する図である。

【図6】実施例１の処理方式に係る数式を示す図である。

【図7】実施例１の構造マップの一例を示す図である。

【図8】実施例１の構造マップの他の例を示す図である。

【図9】実施例１の構造マップの他の例を示す図である。

【図10】実施例１の構造抽出フィルタの一例を示す図である。

【図11】実施例１の組織構造の検出の一例を示す図である。

【図12】実施例１のマスクの一例を示す図である。

【図13】実施例１の構造マップの他の例を示す図である。

【図14】実施例１のマスクの他の例を示す図である。

【図15】実施例１の構造マップの他の例を示す図である。

【図16】実施例１の相関演算の一例を説明するための図である。

【図17】実施例１の計測領域設定のＧＵＩの一例を示す図である。

【図18】実施例１の計測領域の設定を説明する図である。

【図19】実施例１の計測領域の調整を説明する図である。

【図20】実施例１の計測領域の適正程度を示すアラート機能を説明する図である。

【発明を実施するための形態】

【0014】

以下、本発明の実施形態を図面に従い説明する。なお、本明細書において、組織の弾性情報とは、歪、せん断波速度、縦波速度、ヤング率、剛性率、体積弾性率、ポアソン比、粘性率など組織の物性値全般を指すこととする。また、本発明に係る超音波診断装置でせん断波速度を測定する検査対象である、血管や、局所的な線維組織や、のう胞等の構造物を、組織構造、又は構造物と呼ぶこととする。

【実施例1】

【0015】

実施例１は、超音波を送受信する探触子と、探触子から検査対象に超音波を送受信し、検査対象から得られる受信データを処理する処理部を備え、処理部は、探触子から第１の超音波を検査対象に送受信して得られた受信データにより、検査対象の組織構造を検出し、検出した組織構造に基づいて、せん断波速度を検出する計測領域を決定し、計測領域に第２の超音波を送信して、せん断波を発生させ、計測領域に第３の超音波を送受信して得られた受信データにより、せん断波速度を算出する超音波診断装置の実施例である。

【0016】

また、実施例１は、弾性評価方法であって、超音波を送受信する探触子から第１の超音波を検査対象に送受信して得られた受信データにより、検査対象の組織構造を検出し、検出した組織構造に基づいて、せん断波速度を検出する計測領域を決定し、計測領域に第２の超音波を送信して、せん断波を発生させ、計測領域に第３の超音波を送受信して得られた受信データにより、せん断波速度を算出する弾性評価方法の実施例である。ここで、検出した組織構造に基づいてとは、例えば、検出した組織構造を除く領域で、などである。

【0017】

第１の実施例であるせん断波速度の計測法、及び超音波診断装置の構成例について、図１のブロック図を用いて説明する。本実施例の超音波診断装置の構成において、上記の処理部は、後で詳述するように、制御部１５と信号処理部２４の総称である。

【0018】

まず本実施例で利用するＲＦデータおよび画像データの生成に関わる構成について説明する。図１に記載の検査対象の体表面に設置させた、超音波を送受信する探触子１１に対して、超音波信号を生成する送信ビームフォーマ（送信ＢＦ）１３から送信パルス用の電気信号が、図示を省略したデジタルアナログ（Ｄ／Ａ）変換器を経て探触子１１に送られる。探触子１１に入力された電気信号は内部に設置されたセラミック素子にて、電気信号から音響信号に変換され、被検体内に送信される。送信は複数のセラミック素子で行ない、被検体内の所定の深度で集束するように、各素子には所定の時間遅延が掛けられる。

【0019】

検査対象の内部を伝播する過程で反射した音響信号は再び探触子１１にて受信され、送信時とは逆に電気信号に変換され、図示を省略したアナログデジタル（Ａ／Ｄ）変換器を経て、受信した超音波信号から複素のRFデータを生成する受信ビームフォーマ（受信ＢＦ）１４に受信データとして送られる。送受信の切替は、処理部である制御部１５の制御に基づき、送受切替ＳＷ１２で行なわれる。受信ＢＦ１４では、複数の素子で受信した信号に対して、送信時に掛けた時間遅延を考慮した加算処理（整相加算）が行なわれ、減衰補正等の処理がなされた後、複素のＲＦデータとして、処理部である信号処理部２４の血流の速度や方向を示すドプラ画像を生成するドプラ画像生成部１６、またはＲＦデータから組織の形態情報を示す輝度モード（Ｂ）画像を生成するＢ画像生成部１７に送られる。

【0020】

信号処理部２４が受信した受信データであるＲＦデータは、最終的に表示部２０に表示される画像データのうち、超音波の送受信方向に沿った特定の１ラインの要素データとなる。検査対象に対する超音波の送受信を、探触子１１を構成するセラミック素子の配列方向に順次切り替えて実施することにより、画像データの構成要素となる全ての受信データとしてＲＦデータが取得される。

【0021】

取得した受信データであるＲＦデータは信号処理部２４のＢ画像生成部１７において、ゲイン制御、対数圧縮、包絡線検波など、普及している超音波診断装置で一般的に用いられている画像生成処理が実施され、検査対象の内部の形態情報を示すＢ画像が生成される。

【0022】

一方、信号処理部２４のドプラ画像生成部１６では相関演算により血流情報（速度や方向）が算出され、ドプラ画像が生成される。なお、ドプラ画像の生成に関しては超音波の送受シーケンスがＢ画像の方式と異なるが、技術内容は一般的に既知の内容であるため詳細説明は省略する。Ｂ画像およびドプラ画像は、シネメモリ（cine memory）１８に記憶される。生成されたＢ画像やドプラ画像は、スキャンコンバータ１９にて探触子の種類に応じた座標変換と画素補間が為され、これらの画像、評価した画像や数値を表示する表示部２０に表示される。

【0023】

図１に示すように、信号処理部２４は、更に組織構造を検出する構造検出部２１、計測領域を調整する計測領域調整部２２、計測領域で計測したせん断波速度を用いて計測領域の弾性を評価する弾性評価部２３を備える。

【0024】

信号処理部２４の構造検出部２１では受信ＢＦ１４からの受信データを利用して、せん断波の散乱に影響する組織構造を検出し、組織構造を強調した構造マップを生成する。すなわち、信号処理部２４の構造検出部２１は、組織構造の位置と大きさを示す構造マップを作成する。作成された構造マップは表示部２０に表示され、術者がせん断波を計測する領域を指定する際のサポート情報となる。

【0025】

なお、図１に示した超音波診断装置の構成にあっては、構造検出部２１はＦＲデータを利用して構造マップを生成しているが、後で順次説明するように、図２の超音波診断装置の変形構成例に示す通りシネメモリ１８に保存される画像データを用いることもできる。

【0026】

信号処理部２４の計測領域調整部２２は、生成された構造マップの情報を利用して、組織構造の影響が少ない領域を自動的に判定する。計測領域が組織構造を含む場合、計測領域から組織構造の影響が少ない領域を判別し、組織変位を計測する範囲を限定する。

【0027】

信号処理部２４の弾性評価部２３は、計測領域調整部２２で判定した計測領域にてバースト波送信とトラックパルスの送受信を行ない、組織弾性に関係する情報を評価する。

【0028】

図１、２に示した本実施例の超音波診断装置の構成において、装置本体内のデータフローや処理全般を制御する制御部１５と信号処理部２４は、中央処理部（Central Processing Unit：ＣＰＵ）とメモリを備える通常のコンピュータ構成によって実現可能である。すなわち、シネメモリ１８と、必要に応じてハードウェアで実現されるスキャンコンバータ１９を除き、ＣＰＵのプログラム処理で実現することができる。そこで、本明細書においては、信号処理部２４に制御部１５を併せて処理部と呼ぶこととする。

【0029】

次に、図４の機能説明図を利用して、本実施例の構成における信号処理部２４の弾性評価部２３の詳細について説明する。上述のように、本機能処理は、ＣＰＵにおけるプログラム処理によって実現できる。

【0030】

同図に示すように、弾性評価部２３は、前記計測領域制御にて確定した計測領域内にて、放射圧を発生させるための超音波バースト波の送波条件である、集束位置、送信角度、バースト長、電圧、周波数、駆動素子数など送波に必要な音響パラメータ等を確定する第２超音波制御機能２５と、座標情報に基づき、組織の変位を計測するための超音波パルス波（トラックパルス）の送波条件である、集束位置、送信角度、波数、電圧、周波数、送受信回数、駆動素子数など送波に必要な音響パラメータ等を確定する第３超音波制御機能２６と、超音波送受信部から出力されるＲＦデータを利用して組織の変位を計測する変位計測機能２７と、変位計測機能の結果を利用してせん断波速度を計測する速度計測機能２８と、速度計測機能の結果を利用して組織の弾性情報（ここでは、歪、せん断波速度、縦波速度、ヤング率、剛性率、体積弾性率、ポアソン比、粘性率など物質の変形や流動に関する物性値全般を指す）を評価する弾性評価機能２９を備える。なお、せん断波の到来時間は、変位計測機能２７において計測した変位の時間変化から、最大値、最小値、最大値と最小値の中間値などを利用して算出することができる。

【0031】

同図において、まず第２超音波制御機能２５により、指定された計測領域の位置座標に基づきプッシュパルスの送波条件が確定される。生体への影響がなく、効果的にせん断波を発生させる送波条件は、概ねFナンバ=1~2（口径の幅を焦点深度で割った値）の集束条件が適当であり、強度およびバースト長として強度は0.1~1kW/cm²、バースト長は100〜1000μsの範囲が適している。

【0032】

ここで口径の幅は、実際には駆動するセラミック素子の範囲であり、素子間隔の離散的な値を取る。そして理想的な焦点領域を形成するため、各素子への印加電圧には口径重み（アポダイゼーション）が掛けられる。口径の中心から隅にかけて重みを減らすことにより、回折の影響による焦点領域の乱れを抑制する。但し、口径重みは強度を低下させる短所もあるため、評価位置が深部で減衰の影響が大きい場合には、領域形成よりも強度を優先させ、口径重みを軽くする場合もある。また、送信周波数は探触子１１の感度帯域の中心周波数近傍とするのが効果的である。プッシュパルスの送波条件は直ちに制御部１５を介して送信ＢＦ１３に送られ、探触子１１から生体内に照射される。

【0033】

続いて、第３超音波制御機能２６により、トラックパルスの送波条件が確定される。周波数、波数、Ｆナンバなどの音響パラメータは画像データを生成する時の条件とほぼ同じとなる。検査対象が腹部の場合、周波数は1~5MHz、波数は1~3波、Ｆナンバは1~2の条件が利用される。

【0034】

トラックパルスの送信により取得される生体からの反射信号は、探触子１１を介して受信ＢＦ１４に送られ、複素のＲＦデータが生成される。ＲＦデータは変位計測機能２７に入力され、せん断波の伝搬に伴う組織変位が計測される。変位計測機能２７は、パルス繰り返し周期（Pulse Repetition Time：ＰＲＴ）の時間間隔で取得したＲＦデータ間による複素相関演算により実施される。この場合、本実施例においては、単位時間における変位として粒子速度を算出する。プッシュパルスの送波前のＲＦデータを基準に、変位の絶対値を算出する方式もあるが、粒子速度は、探触子の振れや生体組織の自然な動きに伴う低周波成分を除去し、せん断波を高感度に計測する効果がある。

【0035】

変位計測機能２７における演算は取得した全てのＲＦ信号に対して行なわれ、算出された変位情報に基づいて、速度計測機能２８にてせん断波の速度が計測される。

【0036】

最後に弾性評価機能２９にて、計測したせん断波の速度に基づき、組織の弾性情報（ここでは、歪、せん断波速度、縦波速度、ヤング率、剛性率、体積弾性率、ポアソン比、粘性率など組織の物性値全般を指す）が評価される。

【0037】

次に、図３のフローチャートに基づいて、本実施例の装置の信号処理部２４の構造検出部２１と計測領域調整部２２において、構造マップの生成から計測領域の自動制御に至る詳細内容を説明する。

【0038】

同図において、工程３１は第１超音波（パルス信号）の送受信によりＢ像を取得し、術者は検査対象の全体像の観察と弾性評価計測に適した撮像面、例えば、1cm程度の大血管など計測に明らかに影響する組織構造が含まれない面の探索を行なう。

【0039】

工程３２ではせん断波の波面形状を乱す原因となる組織構造を検出する。
図５は組織構造を通過する前後の波面を模式的に表わしたものである。同図の左側に示すように、検査対象５１の内部には組織構造５２が存在し、これに向けてせん断波５３が進行している。同図の右側に示すように、組織構造５２を通過した後のせん断波５４は回折や屈折の作用により波面形状が変化している。

【0040】

せん断波の速度を算出する際には、空間上に設置する計測点に波面が到達する時間を利用する。すなわち、波面形状が乱れている場合には、波面の場所によって到達する時間が変わることになり、速度計測の誤差に強く影響する。組織構造検出工程３２で検出の対象となる組織構造は、せん断波の特性上、波が透過しない血管が主な対象となるが、他にも局所的な線維組織やのう胞が対象に含まれる。ここでは簡便のため組織構造として血管を例に説明を続ける。

【0041】

以下に、組織構造の検出の第１〜第４の方法と、その結果得られる構造マップについて説明する。第１の方法は、人間の視覚効果を利用する方法であり、第１の超音波を検査対象に送受信して得られた受信データに基づき、Ｂ画像生成部１７で生成されたＢ画像の輝度（グレースケールでは白黒の濃淡）を反転した反転画像を構造マップとして利用する。画像データが８ビット(bit)の場合、図６の数式１の処理で反転画像が構成される。数式１において、Ｒは画像データを表わす。血管はＢ像上では黒に配色されるのが一般的であるが、反転画像上では白に配色される。白は黒に比べて進行色もしくは膨張色であるため、反転画像上では細かい血管が強調または拡大されて見えるため、血管の視認が容易になる。

【0042】

図７に、第１の方法で用いるＢ画像と反転画像の一例を示す。同図の左側に示すＢ画像７１には白の検査対象７２と黒の血管７３が示されている。一方、同図の右側に示す反転画像７４では配色が反転し、検査対象７５と血管７６は各々黒と白に配色される。図面では判断が難しいが、反転画像７４上では、Ｂ画像７１に比べて血管が強調されて見える。この効果は、超音波診断装置が通常使用される暗室では更に強調される。

【0043】

また、この第１の方法においては、反転画像７４に代え、組織構造に配色をしたカラ−画像を構造マップとして用いることもできる。図８は同図左側のＢ画像８１を、同図右側のカラー画像８４に変更したものであり、Ｂ画像８１には白の検査対象８２と黒の血管８３が示されているが、例えば血管８６を膨張色である赤、検査対象８５を縮小色である青で配色することにより、血管の視認性を向上できる。

【0044】

以上説明した組織構造検出の第１の方法は、シネメモリ１８に記憶されたＢ画像生成部１７で生成したＢ画像であるスキャンコンバージョン前、或いは後の信号を利用するため、装置構成は先に説明した図２に対応する構成を利用する。

【0045】

組織構造検出の第２の方法は、ドプラ画像生成部１６で生成したドプラ画像を構造マップとして利用する。血管を検出する上では最も直接的なアプローチであるが、低速血流に対する感度が低いことや、画像の取得に時間を要する問題がある。時間の問題は重要であり、構造マップの表示に時間が掛かると、計測領域を設定する段階で既に別の撮像面に移動している可能性があり、血管を避けて計測領域を設定する目的を達成できない。撮像技術の進展により前記問題が解決する場合には、構造マップとして有用性が高くなる。すなわち、図９の左側のＢ画像９１を、同図右側のドプラ画像９４に変更したものであり、検査対象９２に対する血管９３に比較し、検査対象９５に対する血管９６の視認性を向上できる。

【0046】

この第２の方法はドプラ画像生成部１６から出力され、シネメモリ１８に記憶されたドプラ画像を利用するため、装置構成は第１の方法同様、図２に対応する構成となる。
この他、超音波診断装置には組織構造を示す表示形態が多数備えられており、この表示を構造マップとして利用することも本実施形態に含む。例えば組織に対して外部から圧力を加え、その応答である歪を計測してカラーマップ表示する機能（エラストグラフィ）も広く知られている。このように、組織弾性に関連する表示画像や、非線形信号のみを抽出した高調波画像も組織としての均質性を評価する上で有用な情報であり、計測領域の自動選定に有用な構造マップとなる。

【0047】

組織構造検出の第３の方法はフィルタ処理を利用する。この方法は、構造マップを、組織構造の形状に合わせたマスクを有するフィルタ処理により作成するものである。
図１０は、このフィルタ処理を説明するための図である。図１０上段の（ａ）はＢ画像１０１、図１０中段の（ｂ）は上段の（ａ）のＡで示すラインの輝度プロファイル１０２、図１０下段の（ｃ）はフィルタ処理で利用するマスク１０３を示す。

【0048】

このマスク１０３は、中心の１画素とそれから等距離（図１０ではＤ）だけ離れた画素を１とし、その他の画素値を０とするリング状の形状を持つ。本フィルタの特性としては、図１０中段の（ｂ）に示すように、マスク１０３の中心画素で輝度プロファイル１０２のピーク値を検知し、その周囲にリング状に配置した画素（リング画素）で輝度プロファイル１０２の形状を検知する。輝度プロファイル１０２が円形でない場合や、Ｄより広い半値幅を持つ場合、中心画素とリング画素で検知する輝度プロファイル１０２の輝度値の差異が小さくなる。つまり円形の形状で、かつＤより狭い半値幅を持つ組織構造が、中心画素とリング画素で検知する画素の輝度差が大きくなる。この輝度差を利用することで、所望の形状を有する組織構造である組織構造物のみを検知した構造マップが構成される。

【0049】

図１１を利用して、第３の方法の具体的な処理フローを説明する。同図において、まず弾性評価に利用するＢ画像１１１上に関心領域(Region of Interest：ＲＯＩ)１１２を設置し、抽出画像１１３を得る。図６の数式１に従って輝度を反転させて反転画像１１４を得、これを行列Ｒと表わす。続いて、マスクの中心画素に相当する第１マスク１１５（Ｍ１）とマスクのリング画素に相当する第２マスク１１６（Ｍ２）を用意する。両マスクサイズは、反転画像１１４と同じにする。

【0050】

図６の数式２および数式３に従って、反転画像１１４と第１マスク１１５、反転画像１１４と第２マスク１１６の成分同士の積算を実施し、その最大値（Ｉ_ｃおよびＩ_ｎ）を算出する。次に図６の数式（４）によりＩ_ｃおよびＩ_ｎの差分値を算出する。この差分値を設置したＲＯＩ１１２の代表値とする。

【0051】

以上の処理をＢ画像１１１の全ての画素で実施し、設置するＲＯＩの代表値を２次元プロットすることで第３の方法の構造マップが形成される。構造マップは、半径Ｄの円形形状の血管に強い輝度を持ち、Ｂ画像を直接利用する場合に比べて簡易に血管位置を判定できる。

【0052】

上述の説明では血管検出に特化した内容であるが、輝度と形態の情報を利用して特定の組織構造である組織構造物を抽出するという意味では、腫瘍や線維化、石灰化の検出にも適用可能である。なお、腫瘍、線維化、石灰化は、血管とは逆で画像上で高輝度となるため、図１１に示す処理フローにおいて、輝度を反転させる反転処理は不要である。

【0053】

図１２左側に、マスクの中心画素を広げた形態１２１を示す。この場合、輝度プロファイルのピーク値の位置に対して冗長性を持ち、構造マップで描出される血管をより強調できる。また、図１２右側はマスクを簡略化した形態１２２であり、この場合は処理時間の短縮が実現する。

【0054】

図１３は第３の方法のフィルタ処理を実施した結果の一例を示す。Ｂ画像１３１には検査対象１３２とその構造物１３３、１３４が模式的に表わされている。構造物１３３に対応するマスクを設置した場合、これと異なる形状や大きさを有する構造、例えば構造物１３４はフィルタ処理によって検知されず、構造マップ１３５が取得される。

【0055】

なお、マスクの形状は任意に変更可能であり、例えば図１４左側のマスク１４１、または図１４右側に示すマスク１４２の形状によれば、図１５左側のＢ画像１５１から、同図右側の構造マップ１５２のようにライン形状の構造を検出できる。以上説明した第３の方法は、受信ＢＦ１４の直後のＲＦデータ、またはスキャンコンバータ１９の前後のビデオデータのいずれのデータ形式にも対応する。そのため、第３の方法の装置構成は、図１および図２の両構成に対応する。

【0056】

組織構造検出の第４の方法は、相関演算を利用する。この方法は、構造マップを、空間的に異なる位置で取得する複数の前記受信データ、又は空間的に同じ位置で異なる時間で取得する複数の受信データを利用して、相関演算により算出するものである。
図１６上左側に示すように、検査対象の特定の位置にて超音波信号を送受信し、続いて位置を僅かに移動させて再度超音波信号を送受信する。そして、取得した各超音波信号に対して、図６の数式５の相関演算を実施する。血管内部など輝度が低い構造物では空間上での輝度や形状の相関性が低くなるため、演算結果の相関値は低くなる。

【0057】

別の方法として、図１６上右側に示すように時間を変え、超音波の送受信を同一位置にて二回実施し、相関演算を行なう方式でも同様の結果が得られる。図６の数式５に記載の方法が空間上での相関性を利用しているのに対し、図６の数式６に記載の方法は時間軸上での相関性を利用していることとなる。

【0058】

この第４の方法において、相関値が低い位置を強調した相関値マップをＢ画像に重畳表示することにより、血管位置を示す構造マップが形成される。第４の方法は、第３の方法と同様にＲＦデータおよびビデオデータの双方のデータ形式に対応する。したがって、装置構成は図１および図２に対応する。この第４の方法は、信号輝度に対する感度はＲＦデータの方が高く、構造マップの構成に適している。

【0059】

さて、図３のフローチャートの工程３４では工程３３で作成した構造マップを利用して計測領域の位置を決める。図１７は、図１の超音波診断装置の表示部２０などに表示される、Ｂ画像と構造マップの表示形態の一例を示す。表示画面１７１に、Ｂ画像１７２と構造マップ１７３が並列表示されている。並列表示することにより、術者は検査対象の全体像をＢ画像１７２で把握しながら、計測領域の位置は関係する構造物を抽出した構造マップ１７３を利用して設定できる。また、Ｂ画像１７２上に半透明化した構造マップ１７３を重畳表示する形態の場合は、視線移動が少ないため有用な形態の一例である。なお、同図の表示画面１７１には、術者がマニュアルで領域設定を行うか、或いは自動設定を行うかの選択を行う領域設定ボタン１７４を設ける。
なお、設定した計測領域は、ポインタ等の一般的に知られる外部入力機能により、術者が自由に微調整できる。微調整した結果、計測領域に本来は避けるべき構造物が含まれた場合、計測領域の色や明滅や音といった感覚的に判る形態で術者にアラートとして認識させる機能を有する。構造マップして利用する情報は術者が自由に選択でき、例えばドプラ画像などの血管抽出画像と弾性画像が選択された場合、均質領域として自動認識される領域は、血管構造がなく、歪分布が均質な領域が選択される。これに伴い、微調整によるアラートは例えば多段階になる。例えば図２０には、血管画像２０１と歪画像２０３を重畳した構造マップ２０５を図示している。これに計測領域を設定する場合、領域２０６は最も高い適正を示す太枠、領域２０７は中間の適正を示す中太の枠線、領域２０８は不適を示す細枠として、術者に弾性評価に対する適正度合を認識させる画像を画面に表示する。このような画面表示や音といったアラート機能は、例えば計測領域調整部２２により実現され、構造マップを表示しない、もしくは重畳しない並列表示するときに特に効果的である。構造マップをＢ像に重畳させた場合、時には情報過多となり、Ｂ像による全体認識が難しくなる場合が考えられる。この場合には、構造マップを表示しないか、並列表示させる必要があり、領域設定の適正度合は、アラート機能のみで効果的に実行される。

【0060】

図１８に本実施例の装置により、計測領域を設定する際の形態の一例を示す。同図に示すように、バースト波送信およびトラックパルスを送受信する範囲を示す計測領域１８１を、Ｂ画像または構造マップ上で設置する。弾性評価を実施する位置は、術者が手動で行なうか、または、構造マップ上で表示されている構造物を避ける形で自動的に設置する。この計測領域１８１が構造物を避けた形で設定されたもの、言い換えるなら計測領域の正当性を示すため、色の違い、色の濃淡、色の透明度で視覚的に表示することができる。

【0061】

上述の通り、図１７の表示画面１７１上に領域設定ボタン１７４が備えられており、手動か自動かを任意に設定できる仕組みになっている。弾性評価を実施する位置が確定すると同時に、せん断波を計測する計測する計測領域１７５が確定する。

【0062】

図１８の（ａ）に示すように、この計測領域１７５内に構造物を含まない場合と、図１８の（ｂ）に示すように、構造物を含む場合とでそれぞれ表示する計測領域１８１、１８２の配色や濃度、透明度などを変化させることにより、計測領域の正当性を視覚的に認識させ、より最適な位置決めが可能となる。更に、構造物がせん断波の波面に与える影響は、伝搬するせん断波の波長に強く依存する。そのため、計測領域１８２に構造物が含まれる場合、その構造物の大きさに応じて配色等を変化させる方法も、位置決めを行なう上で重要な情報となる。

【0063】

図３のフローチャートに戻り、工程３５では、計測領域の調整を行なう。この計測領域の調整は、信号処理部２４において、得られた構造マップに基づいて、組織構造を避ける位置に計測領域を限定し、または、計測領域を移動させることによって行われる。すなわち、工程３４での位置決めの際に、計測領域が構造物を含む場合の処理工程である。第１の調整方法は、計測領域を限定する方法である。図１９左側には破線で示す調整前の計測領域１９１と、調整後の計測領域１９２、構造物１９３を示す。調整後の計測領域１９２は、構造物１９３を含む伝搬経路を除く範囲に限定して設置されている。その結果、組織変位を計測する範囲は狭くなり、その分のデータ数も低減するが、予め確定されている計測位置を変更せずに計測が実施される。

【0064】

第２の調整方法は、計測領域を移動させる方法である。図１９右側には、破線で示す調整前の計測領域１９１、構造物１９３、調整前の焦域１９４、および、調整後の計測領域１９６、調整後の焦域１９５が表わされている。ここで焦域とは音響放射圧が発生する領域をさす。この第２の調整方法では、計測領域のサイズを変えずに、構造物を含まない位置に移動させる。具体的には、計測領域の中心位置と構造物の中心位置とのベクトル、すなわち距離と方向の情報を算出し、このベクトルと逆方向に移動することで、移動量を最小化できる。第２の調整方法では、焦域も移動させるため、ベクトル情報は装置本体の制御部１５に送られ、探触子１１の送信条件である駆動チャンネル数やフォーカス深度の制御が実施される。

【0065】

続いて、図３の工程３６では音響放射圧およびせん断波を発生させるための第２超音波（バースト波）の送信が行なわれ、工程３７ではせん断波の伝搬に伴う組織の変位量を計測するための第３超音波送受信が実施される。工程３８ではせん断波の速度情報に基づき、検査対象の弾性評価が実施され、表示部２０に数値が提示される（工程３９）。

【0066】

以上、説明してきた種々の実施の形態により、組織構造による波面乱れを抑制し、信頼性が高い組織弾性評価が実現される。その結果、超音波診断装置の検査時間の短縮、術者や患者負担の軽減、診断確度の向上が期待できる。
なお、本発明は上述した実施の形態に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施の形態は本発明のより良い理解のために詳細に説明したのであり、必ずしも説明の全ての構成を備えるものに限定されものではない。また、ある実施の形態の構成の一部を他の実施の形態の構成に置き換えることが可能であり、また、ある実施の形態の構成に他の実施の形態の構成を加えることが可能である。例えは、図１の構成と図２の構成を合わせ持つ信号処理部を利用することができる。また、各実施の形態の構成の一部について、他の構成の追加・削除・置換をすることが可能である。

【0067】

更に、上述した各構成、機能、処理部等は、それらの一部又は全部を実現するプログラムを作成する例を中心に説明したが、例えばスキャンコンバータなど、それらの一部又は全部を例えば集積回路で設計する等によりハードウェアで実現しても良いことは言うまでもない。

【符号の説明】

【0068】

１１ 探触子
１２ 送受切替ＳＷ
１３ 送信ＢＦ
１４ 受信ＢＦ
１５ 制御部
１６ ドプラ画像生成部
１７ Ｂ画像生成部
１８ シネメモリ
１９ スキャンコンバータ
２０ 表示部
２１ 構造検出部
２２ 計測領域調整部
２３ 弾性評価部
２４ 信号処理部
５１、７２、８２、９２、１３２、１６１ 検査対象
５２、７３、８３、９３、１６２、１９３ 構造物
５３ 通過前のせん断波の波面
５４ 通過後のせん断波の波面
７１、８１、９１、１１１、１３１、１５１、１７２ Ｂ画像
７４、１１４ 反転画像
７５ 反転画像の検査対象
７６ 反転画像の構造物
８４ カラー画像
８５ カラー画像の検査対象
８６ カラー画像の構造物
９４ ドプラ画像
９５ ドプラ画像の検査対象
９６ ドプラ画像の構造物
１１２ ＲＯＩ
１１３ 抽出画像
１１５ 第１マスク
１１６ 第２マスク
１３３ 検出構造物
１３４ 非検出対象物
１３５、１５２ 構造マップ
１７１ 表示画面
１７３ 構造マップ
１７４ 領域設定ボタン
１７５ 計測領域
１８１ 構造物を含まない計測領域
１８２ 構造物を含む計測領域
１９１ 調整前の計測領域
１９２ 調整後の計測領域
１９４ 調整前の焦域
１９５ 調整後の焦域
１９６ 調整後の計測領域
２０１ 血管画像
２０２ 血管
２０３ 歪画像
２０４ 歪均質領域
２０５ 血管と歪の複合情報を含む構造マップ
２０６ 最適な計測領域
２０７ 中間適正の計測領域
２０８ 不適正な計測領域

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】

【図19】

【図20】