(57)【要約】【課題】被検体組織を精度よく判別可能な、被検体内部における超音波の減衰情報を提供する。【解決手段】振動子アレイから被検体の組織に、２つの周波数の超音波を送信し、被検体を透過した超音波の透過波を受信して受信信号を出力する。受信信号に基づいて、被検体の組織を透過する際の単位長さ当たりの超音波の減衰量である減衰率を２つの周波数についてそれぞれ求め、２つの周波数間の減衰率の変化率である２周波数間減衰率を求める。２周波数間減衰率を用いることにより、精度よく被検体組織を識別できる。【選択図】図１
2020099619.tif 000002
【特許請求の範囲】
【請求項１】
被検体の組織に、２つの周波数の超音波を送信し、前記被検体を透過した前記超音波の透過波を受信して受信信号を出力する振動子が複数配列された振動子アレイと、
前記受信信号を処理する信号処理部とを有し、
前記信号処理部は、前記受信信号に基づいて、前記被検体の組織を透過する際の単位長さ当たりの前記超音波の減衰量である減衰率を前記２つの周波数についてそれぞれ求め、前記２つの周波数間の減衰率の変化率である２周波数間減衰率を求める２周波数間減衰率算出部とを有することを特徴とする超音波送受信装置。
【請求項２】
請求項１に記載の超音波送受信装置であって、前記２周波数間減衰率を指標として、前記被検体の組織を判別する判別部をさらに有することを特徴とする超音波送受信装置。
【請求項３】
請求項１に記載の超音波送受信装置であって、前記２周波数間減衰率算出部は、前記組織の前記２周波数間減衰率を表示することを特徴とする超音波送受信装置。
【請求項４】
請求項１に記載の超音波送受信装置であって、前記２周波数間減衰率算出部は、求めた前記２つの周波数についての超音波の減衰率の差を求め、求めた減衰率の差を２つの周波数の差で除することにより、前記２周波数間減衰率を求めることを特徴とする超音波送受信装置。
【請求項５】
請求項１に記載の超音波送受信装置であって、前記２周波数間減衰率算出部は、前記２つの周波数にて前記減衰率を求め、減衰率と周波数とをそれぞれ軸とする２次元空間に、求めた前記減衰率と周波数の関係を示す複数の点をプロットし、前記点を所定の数式にフィッティングすることにより、前記２周波数減衰率を求めることを特徴とする超音波送受信装置。
【請求項６】
請求項５に記載の超音波送受信装置であって、前記２周波数間減衰率算出部は、前記所定の数式として、前記減衰率をAttn、前記周波数をfreqとし、ａ１、ａ２、ｂおよびｃは数値である場合、数式
Attn＝ａ１×freq＋ｂ（ただし、ａ１およびｂ≠０）、または、
Attn＝ａ２×（freq）ｃ（ただし、ａ２≠０、ｃ≠１）
を用い、前記数式におけるａ１またはａ２を求めて、前記２周波数間減衰率の値とすることを特徴とする超音波送受信装置。
【請求項７】
請求項１に記載の超音波送受信装置であって、前記信号処理部は、前記被検体の組織に関心領域を設定する関心領域設定部をさらに有し、
前記２周波数間減衰率算出部は、前記関心領域内の前記超音波の組織について、前記２周波数間減衰率を求めることを特徴とする超音波送受信装置。
【請求項８】
請求項７に記載の超音波送受信装置であって、前記振動子アレイは、前記被検体の組織に複数の方向から超音波を送信し、前記被検体を透過した前記透過波を前記複数の方向ごとに受信して透過信号を出力し、
前記信号処理部は、前記複数の方向ごとの前記透過信号を用いて、前記被検体における単位長さ当たりの減衰量である減衰率の値の分布を示す減衰画像を、前記２つの周波数についてそれぞれ生成する減衰画像生成部をさらに有し、
前記関心領域設定部は、前記減衰画像生成部が生成した前記２つの周波数についての前記減衰画像上の対応する位置に前記関心領域をそれぞれ設定し、
前記２周波数間減衰率算出部は、前記２つの周波数についての前記減衰画像の前記関心領域内の前記減衰率の値を用いて前記２周波数間減衰率を求めることを特徴とする超音波送受信装置。
【請求項９】
請求項８に記載の超音波送受信装置であって、前記振動子アレイは、前記透過波の他に、前記被検体によって反射された反射波を受信して反射信号を出力し、
前記信号処理部は、前記反射信号を用いて前記被検体における前記反射波の強度の分布を示す反射画像を生成する反射画像生成部、および、前記透過信号を用いて前記被検体を伝搬する前記超音波の音速の分布を求め、音速分布を示す音速画像を生成する音速画像生成部のうちの少なくとも一方をさらに有し、
前記関心領域設定部は、前記反射画像または前記音速画像を処理することにより前記被検体の組織の輪郭を抽出し、抽出した前記組織の輪郭内の領域を前記関心領域として前記減衰画像に設定することを特徴とする超音波送受信装置。
【請求項１０】
請求項８に記載の超音波送受信装置であって、前記関心領域設定部は、前記減衰画像を複数のセグメントに分割することにより、それぞれの前記セグメントを領域として設定し、
前記２周波数間減衰率算出部は、前記セグメントごとに前記２周波数間減衰率を算出し、前記２周波数間減衰率を前記セグメント内の画素の画素値とする２周波数間減衰率画像を生成することを特徴とする超音波送受信装置。
【請求項１１】
請求項２に記載の超音波送受信装置において、前記信号処理部は、前記透過波の受信信号を用いて前記被検体の組織を伝搬する前記超音波の音速を算出し、前記判別部は、前記２周波数間減衰率と前記音速とを指標として、前記組織を判別することを特徴とする超音波送受信装置。
【請求項１２】
請求項２に記載の超音波送受信装置おいて、前記判別部は、組織の種類が既知の複数の試料について予め求められた前記２周波数間減衰率と、前記試料についての１種以上の物性値との関係を示す点を、前記２周波数間減衰率と前記１種以上の物性値とをそれぞれ軸とする空間にプロットし、前記複数の試料ごとにプロットされた複数の前記点に基づき、前記空間上に前記組織の種類を隔てる線または面を設定することにより、前記組織の種類ごとの領域が生成された判別用空間を記憶しており、
前記判別部は、前記２周波数間減衰率算出部が前記被検体の組織について求めた前記２周波数間減衰率と、前記被検体の組織について求められた１種以上の物性値とを示す点を前記判別用空間にプロットし、前記プロットした点が含まれる前記領域の示す組織の種類を、前記被検体の組織の種類であると判定する組織種類識別部を備えることを特徴とする超音波送受信装置。
【請求項１３】
請求項１２に記載の超音波送受信装置において、前記判別部は、前記判定用空間に設定された前記線または面のうち、悪性の組織と正常の組織とを隔てる線または面から、前記被検体の組織を示す前記点までの距離を求め、前記距離を前記被検体の組織の悪性度として表示する悪性度算出部をさらに備えることを特徴とする超音波送受信装置。
【請求項１４】
被検体の組織に、複数の周波数成分を含む超音波を送信し、前記被検体を透過した前記超音波の透過波を受信して受信信号を出力する振動子が複数配列された振動子アレイと、
前記受信信号を処理する信号処理部とを有し、
前記信号処理部は、前記受信信号に基づいて、前記被検体の組織についての２種以上の物性値を求める物性値算出部と、前記２種以上の物性値を指標として用いて、前記被検体の組織の悪性度を判別する判別部とを含み、
前記判別部は、組織が悪性であるか正常であるかが既知の複数の試料について予め求められた２種以上の物性値の関係を示す点を、前記２種以上の物性値をそれぞれ軸とする空間にプロットし、前記複数の試料ごとにプロットされた前記点に基づき、前記空間上に悪性の組織と正常の組織とを隔てる線または面を設定することにより、悪性の組織と正常の組織の領域が設定された判別用空間を記憶しており、
前記判別部は、前記物性値算出部が前記被検体の組織について求めた前記２種以上の物性値の関係を示す点を前記判別用空間にプロットし、前記プロットした点が含まれる前記領域、および、前記線または面から前記点までの距離を用いて、前記被検体の組織が悪性かどうか、悪性である場合にはその悪性度を判定することを特徴とする超音波送受信装置。
【請求項１５】
請求項１４に記載の超音波送受信装置であって、前記２種以上の物性値のうちの一つは、予め定めた２つの周波数にて求めた、組織を超音波が透過する際の単位長さ当たりの減衰率の差を、前記２つの周波数の差で除して求めた２周波数間減衰率であることを特徴とする超音波送受信装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波送受信装置に関し、特に、被検体を透過した超音波を用いて被検体の組織を判別する装置に関する。
【背景技術】
【０００２】
医用画像診断装置として、超音波撮像装置が広く用いられている。この装置は、音響トランスデューサ（振動子）を配列したアレイから超音波を被検体に照射し、被検体内の各深さで反射された反射波信号を振動子により受信する。反射波信号は、受信ビームフォーミング処理されることにより、被検体内における反射波の強度分布が算出され、反射波の強度分布を示す反射画像（Ｂモード画像）が生成されて表示される。
【０００３】
また、媒質中に配置した被検体に対して、複数の方向から超音波を照射し、被検体を透過した透過波信号から被検体内部の物性値（音速や減衰率）の分布を求め、物性値の画像（断層画像）を生成する超音波ＣＴ（Computed Tomography）装置も知られている。超音波ＣＴ装置には、例えばリング状に振動子を並べたリングアレイが用いられ、リングアレイの開口に被検体が挿入される。リングアレイの一つの振動子から被検体に超音波が送信され、被検体を透過した透過波は他の振動子により受信される。超音波を送信する振動子をひとつずつずらしながら繰り返し超音波の送信を行い、得られた透過波の受信信号から被検体内の音速や減衰率を求めて再構成処理することにより、音速や減衰率の分布を示す画像を生成する。また、超音波ＣＴ装置では、被検体による反射波を受信することもできるため、反射画像を生成して表示することも可能である。
【０００４】
一方、超音波は、被検体を通過しながら減衰するため、特許文献１には、Ｂモード像とともに、減衰率の時間変化をグラフ等により表示することにより、ユーザが視覚的に減衰率の時間変化を認識できるようにした構成が開示されている。特許文献２には、超音波プローブの非線形的な特性により、反射信号（エコー信号）の減衰が被検体の深度に比例しない現象、すなわち深度方向に非線形に減衰する現象をキャンセルするために、ある周波数のエコー信号の深さに依存した減衰量を、別の周波数のエコー信号の深さに依存した減衰量によって割り、プローブ特性を除去したエコー信号データを得ることが開示されている。
【０００５】
また、懸濁液（スラリー）等の性状（粒子径分布や濃度等）をモニタリングするために超音波を用いる超音波減衰分光法と呼ばれる方法が広く知られている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１７−１１３１４５号公報
【特許文献２】特開２０１２−１７０５１２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
近年、乳癌等の早期癌のスクリーニングを非侵襲で行うために、超音波を用いて癌組織と正常組織とを識別することが望まれている。
【０００８】
しかしながら、超音波撮像装置で撮像されたＢモード像は、画質が操作者に依存し、定量的な画像を得ることが困難であるため、癌組織と正常組織をＢモード像に基づいて識別するのは難しい。
【０００９】
一方、従来の超音波ＣＴ装置では、物性値である音速の分布を示す音速画像が得られるが、組織の音速は被検体によってさまざまに異なるため、被検体によって癌組織と正常組織とを分ける音速閾値が異なり、閾値の設定が難しい。また、乳腺と腫瘍は、音速値の範囲が近く、一部は重なっているため、音速で識別するのは難しいという問題がある。また、超音波ＣＴ装置は、減衰率の分布を示す減衰画像も得られる。減衰率は、被検体の単位長さ当たりの減衰量である減衰率（単位ｄＢ／ｃｍ）を伝搬させた超音波の周波数（ＭＨｚ）により除することにより規格化した減衰率（単位ｄＢ／ｃｍ／ＭＨｚ）が用いられる。しかしながら、減衰画像には、散乱アーチファクトが多く含まれるため、癌組織と正常組織とを識別するのが難しい。
【００１０】
また、超音波減衰分光法等のように超音波を用いて物体の性状を測定する方法も、規格化した減衰率（単位ｄＢ／ｃｍ／ＭＨｚ）を測定して性状の判断を行っている。
【００１１】
本発明の目的は、被検体組織を精度よく判別可能な、被検体内部における超音波の減衰情報を提供することにある。
【課題を解決するための手段】
【００１２】
上記目的を達成するために、本発明は、被検体の組織に、２つの異なる周波数を少なくとも含む超音波を送信し、被検体を透過した超音波の透過波を受信して受信信号を出力する振動子が複数配列された振動子アレイと、受信信号を処理する信号処理部とを有する。信号処理部は、受信信号に基づいて、被検体の組織を透過する際の単位長さ当たりの超音波の減衰量である減衰率を２つの周波数についてそれぞれ求め、２つの周波数間の減衰率の変化率である２周波数間減衰率を求める２周波数間減衰率算出部とを有する。
【発明の効果】
【００１３】
本発明によれば、癌組織と正常組織とで値の差が大きい２周波数間減衰率を算出するため、被検体の組織を精度よく識別することができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施形態１の超音送受信装置の全体構成の概略を示すブロック図である。
【図２】正常組織と悪性組織の単位長さ当たりの減衰率と２周波数間減衰率とを示すグラフである。
【図３】正常組織と悪性組織の単位長さ当たりの減衰率と規格化減衰率とを示すグラフである。
【図４】実施形態１の超音波送受信装置の信号処理部の機能ブロック部である。
【図５】（ａ）異なる周波数ｆ１、ｆ２について生成された単位長さ当たりの減衰率を画素とする減衰画像の例を示す図、（ｂ）減衰画像に設定された関心領域の減衰率を示す説明図、（ｃ）反射画像の例を示す図、（ｄ）音速画像の例を示す図、である。
【図６】実施形態１において（ａ）周波数間減衰率と音速とをそれぞれ軸とする空間に、複数の試料の２周波数間減衰率と音速との関係を示す点をプロットした判別用空間を示す図であり、（ｂ）判別用空間から求めた悪性度の度合をバー上で示す表示例である。
【図７】実施形態２の超音波ＣＴ装置の構成を示すブロック図である。
【図８】実施形態２のリングアレイで受信される被検体１からの透過波と反射波とを示す説明図である。
【図９】実施形態２の超音波ＣＴ装置の動作例を示すフローチャートである。
【図１０】実施形態２において（ａ）振動子の出力した透過信号、（ｂ）透過信号をフーリエ変換した信号、をそれぞれ示すグラフである。
【図１１】実施形態２の超音波ＣＴ装置で生成された（ａ）周波数ごとの減衰画像、（ｂ）反射画像、（ｃ）音速画像、（ｄ）セグメントに分割した減衰画像、（ｅ）２周波数間減衰率画像、のそれぞれ一例を示す図である。
【図１２】乳癌と正常組織の、周波数で規格化した減衰率と音速とを示すグラフである。
【図１３】実施形態３の超音波ＣＴ装置の信号処理部の機能ブロック図である。
【図１４】実施形態３の判別用空間生成部が判別用空間を生成する処理を示すフローチャートである。
【図１５】実施形態３の組織種類識別部の被検体組織の識別処理を示すフローチャートである。
【図１６】実施形態３の超音波ＣＴ装置の処理の流れを示す図であり、（ａ）２周波数間減衰画像と音速画像が生成されることを示す図、（ｂ）２周波数間減衰画像と音速画像から組織識別を行うこと示す図、（ｃ）組織の悪性度合画像が生成されることを示す図、である。
【発明を実施するための形態】
【００１５】
以下、本発明の一実施形態について図面を用いて説明する。
【００１６】
＜実施形態１＞
まず、実施形態１の超音波送受信装置について図１、図２を用いて説明する。図１は、実施形態１の超音波送受信装置の構成を示している。図２は、生体組織（正常組織と悪性組織）の減衰率の測定値を周波数との関係で示したグラフの一例である。
【００１７】
本実施形態の超音波送受信装置は、図１に示すように振動子アレイ３０と信号処理部２０とを備えている。振動子アレイ３０は、複数の振動子３が配列された構造であり、振動子３は、被検体１の組織に、予め定めた２つの異なる周波数を少なくとも含む超音波を送信し、被検体１を透過した超音波の透過波を受信して受信信号を出力する。
【００１８】
信号処理部２０は、振動子３が出力する受信信号を処理し、被検体１の組織を判別するための構成として、２周波数間減衰率算出部２１を備えている。２周波数間減衰率算出部２１は、まず、振動子２の受信信号に基づいて、送信した超音波に含まれる予め定めた２つの周波数（第１周波数f1、第２周波数f2）について、被検体１の組織を透過する際の単位長さ当たりの超音波の減衰量である減衰率ACf1、ACf2をそれぞれ求める。さらに、２周波数間減衰率算出部２１は、求めた２つの周波数f1,f2間の減衰率ACf1、ACf2の減衰率の変化率a2fを求める。本実施形態では、２つの周波数間の減衰率の変化率を「２周波数間減衰率」と呼ぶ。
【００１９】
なお、第1周波数f1および第2周波数f2の両方が、振動子３の送信する超音波の帯域内に含まれている場合、超音波の送信回数を低減することができるため、組織の判別に係る時間を低減できるため好ましい。ただし、第1周波数f1を含む超音波を送信した後、第2周波数f2を含む超音波を送信する構成であっても構わない。
【００２０】
２周波数間減衰率算出部２１は、減衰率a2fを求めるために、例えば、２つの周波数f1,f2についての超音波の減衰率ACf1、ACf2の差（ACf2−ACf1）を求め、２つの周波数の差（f2−f1）で除することにより、２周波数間減衰率a2f＝（ACf2−ACf1）／（f2−f1）を求めることができる。
【００２１】
このように、本実施形態では、周波数で規格化した減衰率（例えば、ACf1／f1)ではなく、２周波数間減衰率a2fを算出する。２周波数間減衰率a2fは、周波数で規格化した減衰率よりも、悪性組織（癌組織）と正常組織との値の差が大きい。よって、２周波数間減衰率a2fの値をユーザに表示等することによりユーザは、被検体１の組織が悪性組織か正常組織化を、精度よく識別することができる。
【００２２】
具体的には、図２からわかるように、周波数f1と周波数f2の間の帯域において、正常組織（normal)の単位長さ当たりの減衰率(単位はdB／cm)の測定値は、周波数（ＭＨｚ）が変化することによる変化（傾き）が小さい。そのため、周波数f1と周波数f2との間の帯域における正常組織の２周波数間減衰率an2f（例えばan2f＝（ACnf2−ACnf1）／（f2−f1））は、値が小さい。これに対し、悪性組織（癌：cancer)の単位長さ当たりの減衰率は、周波数f1と周波数f2の間の帯域において、周波数が変化することによる変化（傾き）が正常組織より大きい。そのため、２周波数間減衰率ac2f（例えばac2f＝（ACcf2−ACcf1）／（f2−f1））の値は、正常組織より大きいという特徴を有する。
【００２３】
一方、図３に示したように、図２と同じ正常組織と悪性組織の単位長さ当たりの減衰率(dB／cm)の測定値について、周波数で規格化した減衰率を求めると、正常組織の周波数規格化減衰率ansf＝ACnf／fの方が、悪性組織の周波数規格化減衰率acsf＝ACcf／fよりも値が小さいという点で、２周波数間減衰率の傾向と同様である。しかしながら、図２と図３とを比較すると明らかなように、正常組織と悪性組織の２周波数間減衰率an2f、ac2fの値の差（ac2f−an2f）は、周波数規格化減衰率ansf、acsfの値の差（acsf−ansf）より大きく、正常組織と悪性組織の差を明確に判別する指標として適している。
【００２４】
このように、本実施形態では、悪性組織（癌組織）と正常組織とで値の差が大きい２周波数間減衰率を算出するため、被検体の組織を精度よく識別することができる。
【００２５】
なお、信号処理部２０に判別部２２をさらに配置し、判別部２２が２周波数間減衰率a2fを指標として、被検体１の組織を精度よく判別する構成としてもよい。また、２周波数間減衰率算出部２１が、算出した２周波数間減衰率の値を表示し、ユーザがその表示を見て組織を判別する構成とすることもできる。
【００２６】
なお、２周波数間減衰率算出部２１は、２周波数間減衰率を求めるために、２つの異なる周波数f1,f2を含む複数の周波数について、それぞれ単位長さ当たりの減衰率（dB／cm)を求め、図２に示すような減衰率と周波数とをそれぞれ軸とする２次元空間に、求めた複数の減衰率と周波数の関係を示す複数の点をプロットし、プロットした点を所定の数式にフィッティングすることにより、２周波数間減衰率a2fを求める構成としてもよい。
【００２７】
具体的には、例えば、単位長さ当たりの減衰率をAttn、周波数をfreqとした場合、プロットした点を下式（１）にフィッティングし、求めたａ１を２周波数間減衰率a2fとする。
【００２８】
Attn＝ａ１×freq＋ｂ  ・・・（１）
ただし、ａ１、ｂは数値であり、ａ１およびｂはいずれも０ではない（ａ１≠０およびｂ≠０）。
【００２９】
また、プロットした点を下式（２）にフィッティングし、求めたａ２を２周波数間減衰率a2fとしてもよい。
【００３０】
Attn＝ａ２×（freq）ｃ  ・・・（２）
ただし、ａ２、ｃは数値であり、ａ２≠０、ｃ≠１である。
【００３１】
また、図４に示すように、信号処理部２０は、被検体１の組織に関心領域を設定する関心領域設定部２３をさらに備えていてもよい。これにより、２周波数間減衰率算出部２１は、関心領域設定部２３に設定された関心領域内の超音波の組織について、２周波数間減衰率を求めることができる。
【００３２】
また、信号処理部２０は、図４に示すように、減衰画像生成部２４をさらに備えていてもよい。この場合、図１の振動子アレイ３０は、被検体１の組織に複数の方向から超音波を送信し、被検体１を透過した透過波を受信して受信信号を出力する構成とする。減衰画像生成部２４は、送信方向ごとに、振動子アレイ３０が透過波を受信した受信信号を用いて、２つの超音波周波数f1,f2のそれぞれについて被検体１における減衰率(dB／cm)の値の分布を示す減衰画像５１，５２を生成する（図５（ａ）参照）。
【００３３】
関心領域設定部２３は、減衰画像生成部２４が生成した２つの減衰画像５１，５２上の対応する位置に関心領域５３，５４をそれぞれ設定する（図５（ａ）参照）。２周波数間減衰率算出部２１は、減衰画像５１，５２の関心領域５３，５４内の減衰率ACf1、ACf2の値を用いて２周波数間減衰率a2fを求めることができる。例えば、減衰画像５１の関心領域５３内の画素値が示す減衰率の平均を求め、減衰率ACf1とする。同様に、減衰画像５２の関心領域５４内の画素値が示す減衰率の平均を求め、減衰率ACf2とする。
【００３４】
また、関心領域設定部２３は、減衰画像５１，５２を複数のセグメントに分割することにより、それぞれのセグメントを関心領域として設定してもよい。この場合、２周波数間減衰率算出部２１は、セグメントごとに２周波数間減衰率を算出する。例えば、減衰画像５１，５２の対応するセグメント内の全画素値（減衰率）の平均値等をそれぞれ求め、その値を用いてセグメントごとに２周波数間減衰率を算出する。2周波数間減衰率算出部２１は、２周波数間減衰率をセグメント内の画素の画素値とする２周波数間減衰率画像を生成することができる。ユーザは、２周波数間減衰率画像を見て、セグメントが悪性組織か正常組織かを判断することができる。
【００３５】
また、信号処理部２０は、図４に示すように、反射画像生成部２５を備えていてもよい。この場合、振動子アレイ３０は、透過波の他に、被検体１によって反射された反射波を受信して、その受信信号を出力する。反射画像生成部２５は、反射波の受信信号を用いて被検体における反射波の強度の分布を示す反射画像５５を生成する（図５（ｃ）。関心領域設定部２３は、反射画像５１を処理することにより、被検体１の組織の輪郭を抽出し、抽出した組織の輪郭内の領域を関心領域５３，５４として減衰画像５１，５２に設定することができる。
【００３６】
さらに、信号処理部２０は、図４に示すように音速画像生成部２６を備えていてもよい。音速画像生成部２６は、振動子アレイ３０から透過波の受信信号を受け取って、被検体１を伝搬する超音波の音速の分布を求め、音速分布を示す音速画像５６を生成する。関心領域設定部２３は、音速画像５６を用いて被検体の組織の輪郭を抽出し、抽出した組織の輪郭内を関心領域５３，５４として減衰画像５１，５２に設定することができる。
【００３７】
なお、関心領域設定部２３は、ユーザから受け付けた領域を関心領域５３，５４として設定する構成であってもよい。
【００３８】
判別部２２は、例えば、２周波数間減衰率a2fと、１種類以上の被検体１の物性値とを指標として、組織を判別することができる。物性値としては、音速、弾性率、粘性度等どのような値であってもよい。例えば、２周波数間減衰率a2fと音速とを指標として、被検体１の組織を判別することができる。
【００３９】
判別方法を具体的に説明する。まず、組織の種類（ここでは、癌、乳腺、脂肪）が既知の複数の試料を用意し、２周波数間減衰率a2fと、１種以上の物性値（例えば音速）とをそれぞれ求めておく。そして、図６（ａ）に示したように、２周波数間減衰率と１種以上の物性値（音速）とをそれぞれ軸とする２次元以上の空間に、複数の試料の２周波数間減衰率a2fと物性値（音速）との関係を示す点をプロットし、複数の点に基づき、空間上に組織の種類を隔てる線（３次元空間の場合は面）を設定する。これにより、線または面により分けられた領域が、組織の種類（癌、乳腺、脂肪）にそれぞれ対応する判別用空間（図６（ａ））を生成することができる。判別用空間は、判別部２２内のメモリ等に格納しておく。
【００４０】
判別部２２は、２周波数間減衰率算出部２１が被検体１の組織について求めた２周波数間減衰率a2fと、音速画像生成部２６等が求めた被検体１の組織についての音速等の１種以上の物性値とを示す点を、判別用空間（図６（ａ））にプロットする。判別部２２は、プロットした被検体１の組織の点が含まれる領域の示す組織の種類を、被検体１の組織の種類であると判定する。例えば、被検体１の組織の２周波数間減衰率と音速が示す点が、図６（ａ）の点６３である場合、判別部２２は、その組織は癌であると判定し、点６４である場合、乳腺であると判定する。
【００４１】
さらに、判別部２２は、判定用空間（図６（ａ））に設定された線（または面）のうち、悪性の組織と正常の組織（図６（ａ）では癌と乳腺）とを隔てる線６２（または面）から、被検体１の組織を示す点６３，５４までの距離を求め、距離を被検体１の組織の悪性度として示すことも可能である。例えば、点６３の組織は、癌の領域に位置し、線６２からの距離はＰ１である。よって、図６（ｂ）のように、悪性度の表示バー上で、悪性と正常組織の境界よりも悪性側に距離Ｐ１で表される悪性度であると表示することができる。一方、点６４ｂの組織は、乳腺の領域に位置し、線６２からの距離はＰ２であるので、悪性と正常組織の境界よりも正常組織側に距離Ｐ２で表される悪性度（マイナスの悪性度＝良性度）であると表示する。
【００４２】
＜実施形態２＞
本実施形態２として、実施形態１の超音波送受信装置を含む超音波ＣＴ装置について説明する。ここでは、超音波ＣＴ装置の一例として、乳房の組織を識別するのに適した構造の装置について説明する。
【００４３】
図７に示したように、本実施形態１の超音波ＣＴ装置は、すでに述べた振動子アレイ３０と信号処理部２０の他に、制御部６、入出力部９、記憶部８、被検体１を乗せるベッド２、水槽４、および予備タンク５を備えている。ベッド２には、被検体１を搭載する面に開口が設けられ、開口の下部には、胸部を挿入できる円柱状の水槽４が備えられている。その内部には、図１に示したようなリング状の振動子アレイ３０が、水槽４の軸方向に平行移動可能に備えられている。水槽４には、温水が満たされ、予備タンク５に接続されている。予備タンク５は、水槽４内の温水を浄化、加熱および脱気する。予備タンク５内や水槽４の下部には、温度計（不図示）が取り付けられ、制御部６に接続されている。
【００４４】
超音波ＣＴ装置の撮影条件は、ユーザにより、入出力部９のタッチパネルやキーボード等を通して設定される。設定された条件等は、記憶部８であるメモリやハードディスクドライブ等に保存される。
【００４５】
制御部６は、入出力部９から入力された条件や、記憶部８に保存された撮像条件をもとに、送信信号を生成し、振動子アレイ３０を構成する１つまたは複数の振動子３ａに出力する。これにより、図８に示したように送信信号を受け取った１つまたは複数の振動子３ａから超音波（入射波）が被検体１に向かって送信される。被検体１を透過した透過波は、超音波を送信した１つまたは複数の振動子３ａに対して、対向する所定の角度範囲にある複数の振動子３により受信される。被検体１により反射された反射波は、超音波を送信した振動子３ａを中心に所定の角度範囲にある複数の振動子３により受信される。透過波を受信した振動子３は、透過波の受信信号を、反射波を受信した振動子３は、反射波の受信信号をそれぞれ出力する。透過波および反射波の受信信号は制御部６に入力される。制御部３は、送信信号を出力する振動子３ａの位置をずらすことにより、所定の送信角度（ビュー）ごとに超音波の送信と受信を繰り返し行う。
【００４６】
また、制御部６は、送信信号と受信信号の送受の切り替えを行うスイッチング部と、振動子アレイ３０の上下動の制御を行うコントローラと、予備タンク５の水圧制御や温水の温度制御などを行うコントローラを備えている。
【００４７】
振動子アレイ３０の出力する透過波および反射波の受信信号は、記憶部８に記録されると共に、信号処理部２０で後述するように信号処理が実行される。信号処理部２０において生成した被検体１の画像は、入出力部９に備えられたモニタ等に表示される。なお、制御部６と信号処理部２０と記憶部８は、ベッド２の下部の空間に収めてもよい。
【００４８】
信号処理部２０は、図４に示した各部２１〜２６を含む。信号処理部２０は、プロセッサ（例えば、ＣＰＵ（Central Processing Unit）やＧＰＵ（Graphics Processing Unit））と、プログラムが予め格納されたメモリとを備えて構成され、プロセッサがプログラムを読み込んで実行することにより、信号処理部２０の各部２１〜２６の機能をソフトウエアにより実現する。なお、信号処理部２０は、機能の一部または全部をハードウエアによって実現することも可能である。例えば、ＡＳＩＣ（Application Specific Integrated Circuit）のようなカスタムＩＣや、ＦＰＧＡ（Field-Programmable Gate Array）のようなプログラマブルＩＣを用いて信号処理部２０を構成し、その動作を実現するように回路設計を行えばよい。
【００４９】
以下、本実施形態２の超音波ＣＴ装置の動作をさらに詳しく説明する。図９は、制御部６および信号処理部２０の動作を示すフローチャートである。
【００５０】
本実施形態の超音波ＣＴ装置は、電源が投入されたならば、制御部６は、予備タンク５内の水の温度を温度計から取り込んで、水の温度が所定の温度（体温程度）になるまで加熱し、脱気装置により脱気した後、ポンプを駆動させて、水槽４に移動させる。これにより、水槽４は、所定の温度に調整した脱気水で満たされる。被検体１が挿入されていない状態で、制御部６は、予め定めておいた条件で、超音波の送信および受信を行い、被検体１が挿入される前の受信データを取得しておく。
【００５１】
ユーザが、乳房の組織の識別開始を入出力部９を介して指示した場合、制御部６は、入出力部９のモニタに、被検体１にベッド２へうつ伏せになり、水槽４に片方の乳房を挿入するように促す表示を入出力部９のモニタに表示する（ステップ８０４）。被検体１の乳房が水槽４に挿入されたことを、被検体１による入出力部９の操作等により、制御部６が確認したならば、制御部６は、振動子アレイ３０に以下のように超音波の送受を行う（ステップ８０５、８０６）。ここでは、ピッチが０．５ｍｍで２０４８チャネルの振動子３がリング状に並べられた直径３２６ｍｍの振動子アレイ３０を用いる。なお、振動子３の水槽４の軸方向の厚みは１ｍｍとする。振動子アレイ３０から照射する超音波（入射波）の周波数は、第１周波数f1である１．５ＭＨｚ（水中での超音波の波長約１ｍｍ）と、第2周波数f2の３．０MHｚを含んでいる。
【００５２】
すなわち、ステップ８０５において、制御部６は、各ビューにおいて超音波を送信する１つまたは複数の振動子３ａの番号、振動子３ａが複数である場合には、振動子３ａごとの送信信号の遅延時間を読み込んで、内蔵するメモリから読み出して振動子３ａごとの送信信号と生成し、１つまたは複数の振動子３ａに出力する。これにより、１つまたは複数の振動子３ａから超音波（入射波）が送信される。ステップ８０６において、送信された超音波のうち被検体１を透過した透過波は、振動子３ａに対向する所定の角度範囲の振動子３により受信される。制御部６は、それらの振動子３から受信信号（透過信号）を受け取り、記憶部８に格納する。また、被検体１で反射された反射波は、振動子３ａを中心とする所定の角度範囲の振動子により受信される。制御部６は、それらの振動子から受信信号（反射信号）を受け取り、記憶部８に格納する（ステップ８０６）。制御部６は、超音波の送受信を振動子３ａを所定の角度ずつずらしながら、すべてのビューについて行う。
【００５３】
信号処理部２０の減衰画像生成部２４は、周波数f1,f2を含む所定の周波数帯域の複数の周波数（例えば、０．５MHzから４MHzまでの間の０．１MHzごとの周波数）についての減衰画像を求めるために、まず、記憶部８から各ビューにおける各振動子３の出力した透過信号（図１０（ａ））をそれぞれ読み出してフーリエ変換する。これにより、時間軸方向の信号強度の変化を表す透過信号が、図１０（ｂ）のように、周波数ごとの振幅を表す信号に変換される。これにより、被検体１が挿入された状態で受信した透過信号について、上記した各周波数の振幅が求められる（ステップ８０７）。
【００５４】
減衰画像生成部２４は、ステップ８０４の前に被検体１を挿入することなく受信しておいた各振動子３の受信信号についても同様にフーリエ変換を行い、各周波数の振幅を求める。
【００５５】
減衰画像生成部２４は、被検体１の挿入前後の周波数ごとの振幅の対数の差を、各ビュー、各受信チャネル毎にそれぞれ計算する。このデータの集まりはサイノグラムと呼ばれる。減衰画像生成部２４は、振幅の対数の差のサイノグラムを、Ｘ線ＣＴの分野で広く利用されているフィルタ補正逆投影法（Filtered Back Projection, FBP）等で処理することにより、断層画像を再構成する。これにより、被検体１の挿入前後の減衰率の差の分布画像が得られる。減衰画像生成部２４は、水の減衰率として予め定めておいた値（推定値）を用いることにより、上記減衰率の差の分布画像から、被検体１の減衰率（単位：dB/cm)の分布を示す画像（減衰画像）を生成する。減衰画像生成部２４は、上記処理を周波数ごとのサイノグラムにそれぞれ行って、図１１（ａ）のように周波数ごとの減衰画像を生成する（ステップ８０８）。減衰画像生成部２４は、生成した減衰画像を、記憶部８に格納する。
【００５６】
つぎに、反射画像生成部２５は、反射画像を生成する（図１１（ｂ）参照）。また、音速画像生成部２６は音速画像を生成する（図１１（ｃ）参照）。関心領域設定部２３は、反射画像または音速画像の一方を用いて、減衰率画像を組織ごとに分割する境界線を設定する（ステップ８０９）。以下、本実施形態では、組織ごとに分割する境界線を減衰画像に設定することをセグメンテーションと呼び、分割された領域をセグメントと呼ぶ。各セグメントが関心領域であり、２周波数間減衰率算出部２１は、各セグメントについて２周波数間減衰率を後述するように算出する。反射画像および音速画像は、ユーザが選択した一方のみを生成してもよいし、両方を生成して入出力部９の表示部に表示し、ユーザからセグメンテーションに使用する一方の画像の選択を受け付ける構成であってもよい。
【００５７】
反射画像生成部２５が、反射画像を生成する処理について説明する。反射画像生成部２５は、記憶部８から各ビューにおける各振動子３の出力した反射信号をそれぞれ読み出して、時間方向にヒルベルト変換を実施する。反射画像生成部２５は、振動子３ａから超音波（入射波）を送信してから被検体１からの反射波が振動子３に返って来る時間（タイミング）を、送信した振動子３ａから注目画素（被検体１内の点）までの距離と、注目画素から受信した振動子３までの距離の和を、超音波の音速（例えば水の音速）で割ることで求める。もしくは、このタイミングを予め求めておいて、記憶部８内に格納しておく。そして、反射画像生成部２５は、被検体１内の注目画素で反射した反射波が、反射波の受信を行う各振動子にそれぞれ到達するタイミングで各超音波素子が受信した信号を加算し、加算後の信号強度をその画素の値とする。この方法は、遅延加算法（Delay and Sum, DAS）と呼ばれている。これを視野内の全画素について行うことで、超音波エコー検査で広く利用されているＢモード画像が得られる。反射画像生成部２５は、ビューごとにＢモード画像を生成し、加算することで、被検体１の反射画像が得られる。反射画像生成部２５は、生成した反射画像を記憶部８に格納する。
【００５８】
音速画像生成部２６が、音速画像を生成する処理について説明する。まず、音速画像生成部２６は、記憶部８から各ビューにおける各振動子３の出力した透過信号をそれぞれ読み出し、時間方向にヒルベルト変換を実施し、受信波の最大振幅の受信タイミングを求める。音速画像生成部２６は、被検体１の挿入前に受信しておいた各振動子３の受信信号についても同様に最大振幅の受信タイミングを求める。音速画像生成部２６は、被検体１の挿入前後の受信タイミングの差を、各ビュー、各受信チャネル毎にそれぞれ計算し、サイノグラムを得る。音速画像生成部２６は、受信タイミングの差のサイノグラムをフィルタ補正逆投影法等で処理することにより、断層画像を再構成する。この断層画像は、被検体１の挿入前後の、超音波の「遅さ（Slowness）」の差の分布画像である。「遅さ」は、音速の逆数である。音速画像生成部２６は、水の音速値（推定値）を用いて、「遅さ（Slowness）」の差の分布画像から、被検体１の音速の分布画像（音速画像）を生成する。音速画像生成部２６は、生成した音速画像を、記憶部８に格納する。
【００５９】
関心領域設定部２３は、反射画像または音速画像を用いて、減衰画像のセグメンテーションを行う。反射画像および音速画像は、減衰画像よりも組織の境界が明確に表れるという特徴がある。関心領域設定部２３は、反射画像または音速画像から画像処理（２値化処理等）により組織の境界線を抽出し、抽出した境界線を、周波数ごとの減衰画像上の対応する位置に設定する。これにより、図１１（ｄ）のように、周波数ごとの減衰画像をそれぞれ複数のセグメントに分割する（ステップ８０９）。
【００６０】
２周波数間減衰率算出部２１は、減衰画像のセグメント内の減衰率を求める。例えば、セグメント内のすべての画素値（減衰率）の平均を求めて、これをセグメント内の減衰率とする。また、セグメント内の画素値の最大値や最小値や中央値等の所望の代表値を算出して、これをセグメント内の減衰率としてもよい。これを周波数ごとの減衰画像について行う。これにより、周波数ごとの減衰画像の対応するセグメントの減衰率を得ることができる。
【００６１】
２周波数間減衰率算出部２１は、一つのセグメントについて得られた周波数ごとの減衰率を示す点を、図２に示すように周波数と減衰率を２軸とする２次元空間にプロットし、第１の周波数f1と第２の周波数f2との間の周波数帯域の減衰率を示す点を示す点を所定の数式にフィッティングする。所定の数式は、上述した式（１）でもよいし、式（２）でもよい。式（１）と（２）を以下に再び示す。Attnは、単位長さ当たりの減衰率であり、freqは、周波数である。式（１）において、ａ１≠０およびｂ≠０であり、式（２）において、ａ２≠０、ｃ≠１である。
【００６２】
Attn＝ａ１×freq＋ｂ  ・・・（１）
Attn＝ａ２×（freq）ｃ  ・・・（２）
２周波数間減衰率算出部２１は、式（１）の場合は係数ａ１が、式（２）の場合は係数ａ２を求める。これにより、そのセグメントの組織について、第１の周波数f1と第２の周波数f2との間の周波数帯域における減衰率の変化率、すなわち２周波数間減衰率を求めることができる。
【００６３】
２周波数間減衰率算出部２１は、すべてのセグメントについて２周波数間減衰率を求める（ステップ８１０）。
【００６４】
２周波数間減衰率算出部２１は、２周波数間減衰率の値をセグメント内の全画素の画素値とする画像（２周波数間減衰率画像）を図１１（ｅ）のように生成し、入出力部９内の表示部に表示する（ステップ８１１）。
【００６５】
２周波数間減衰率は、実施形態１で述べたように、正常組織と悪性組織との値の差が、周波数で規格化された減衰率よりも大きい。よって、図１１（ｅ）の２周波数間減衰率画像に基づいて医師が診断を行うことにより、医師は、セグメント（組織）の画素値（２周波数間減衰率）を予め求めておいた閾値と比較する等により悪性組織か正常組織かを精度よく識別することができる。特に、乳癌と乳腺は、周波数で規格化された減衰率で比較すると、図１２に示したように、値の範囲が重なるため、識別することが困難であるが、本実施形態の２周波数間減衰率で比較すると値の差が大きいため、識別が可能になる。
【００６６】
また、本実施形態で得られる２周波数間減衰率画像は、セグメントの輪郭を反射画像または音速画像により生成しているため、医師は、悪性組織の輪郭を明確に把握することができる。
【００６７】
＜実施形態３＞
実施形態３の超音波ＣＴ装置について説明する。
【００６８】
実施形態３の超音波ＣＴ装置は、実施形態２の超音波ＣＴ装置と同様の構成であるが、実施形態２の装置の機能に加えて、判別部２２が、被検体１についての２種類以上の物性値を指標として、被検体１の組織の種類を識別する機能を備えている。指標となる２種類以上の物性値としては、どのような組み合わせであってもよいが、ここでは、２周波数間減衰率と音速とを用いる場合について以下説明する。
【００６９】
具体的には、図１３に示すように、判別部２２は、組織の種類を識別する組織種類識別部２２ａと、組織の悪性度を求めて表示する悪性度算出部２２ｂとを備えている。
【００７０】
また、判別部２２は、実施形態１において図６（ａ）を用いて説明した判別用空間を生成する判別用空間生成部２２ｃをさらに備えている。組織種類識別部２２ａは、被検体１の組織について求めた２周波数間減衰率と音速とを、判別用空間にプロットすることにより、組織の種類を識別する。悪性度算出部２２ｂは、判別用空間にプロットされた被検体１の組織を示す点と、悪性組織と正常組織とを隔てる線または面との距離を求め、悪性度を表示する。
【００７１】
図６（ａ）の判別用空間は、別の装置や方法によって予め求め、判別部２２内の記憶部に格納しておくことも可能であるが、本実施形態３では、判別用空間生成部２２ｃが、判別用空間を生成する。この処理について図１４を用いて説明する。
【００７２】
まず、操作者は、予め組織の種類が既知でその組織の領域が既知の複数のサンプルを用意する。例えば、乳癌と乳腺と脂肪の組織の領域が既知のサンプルを複数用意する。
【００７３】
判別用空間生成部２２ｃは、実施形態２のステップ８０４〜８０６と同様にサンプルを水槽４に配置するようにユーザに促し（ステップ１３１、１３２）、サンプルが配置されたならば、超音波を送信し、被検体１の透過波および被検体１による反射波の受信を行う（ステップ１３３、１３４）。
【００７４】
信号処理部２０の減衰画像生成部２４は、透過波の受信信号に対し、実施形態２のステップ８０７〜８０８の処理を行い、周波数ごとに減衰率画像を生成する（ステップ１３５）。
【００７５】
また、音速画像生成部２６は、実施形態２のステップ８０９で説明した処理により音速画像を生成する（ステップ１３６）。反射画像生成部２５は、実施形態２のステップ８０９で説明した処理により反射画像を生成する（ステップ１３７）。
【００７６】
関心領域設定部２３は、実施形態２のステップ８０９と同様の処理により、反射画像または音速画像を用いて、減衰画像をセグメントに分割し、サンプルの組織の種類が明らかな１以上の領域に対応するセグメントをそれぞれ関心領域（以下、ＲＯＩと称す）として設定する。２周波数間減衰率算出部２１は、実施形態２のステップ８０９と同様の処理により、減衰画像のＲＯＩとしたセグメント内の２周波数間減衰率を算出する（ステップ１３８）。
【００７７】
関心領域設定部２３は、音速画像において、ステップ１３８でＲＯＩに対応する位置にＲＯＩを設定する。判別用空間生成部２２ｃは、音速画像のＲＯＩ内の音速を求める。例えば、ＲＯＩ内の全画素が示す音速の平均値、最大値、最小値、中央値等のうち、予め定めた値をＲＯＩの音速として求める（ステップ１３９）。
【００７８】
判別用空間生成部２２ｃは、ステップ１３８、１３９で設定されたＲＯＩごとに、組織の種類をそのサンプルについて予めわかっている情報を読み出すことにより求める（ステップ１４０）。
【００７９】
判別用空間生成部２２ｃは、ステップ１３８、１３９でＲＯＩごとに求めた２周波数間減衰率と音速とを、ステップ１４０において読み出した組織の種類の値であるとして、２周波数間減衰率と音速とをそれぞれ軸とする２次元空間である判別用空間に図６（ａ）のようにプロットする（ステップ１４１）。
【００８０】
この処理を、すべてのサンプルについて繰り返し行う（ステップ１３１）。これにより、すべてのサンプルに含まれる種類が既知の組織の２周波数間減衰率と音速とを示す点が、図６（ａ）のようにプロットされる。
【００８１】
判別用空間生成部２２ｃは、図６（ａ）のように、プロットされたサンプル点を組織の種類によって隔てる境界線６１，６２を判別用空間に設定する。境界線６１、６２によって隔てられた領域は、被検体の組織の種類にそれぞれ対応している。プロットされた組織の種類を隔てる境界線６１，６２の設定方法としては、公知の種々の方法を用いることができるが、例えば線形判別分析法を用いることができる。
【００８２】
図６（ａ）の例では、指標が２つ（２周波数間減衰率と音速）であるため、判別領域が２次元平面であり、境界線６１，６２は線であるが、指標を３以上用いてもよく、その場合には、境界面を設定する。
【００８３】
つぎに、組織種類識別部２２ａと悪性度算出部２２ｂは、上述のようにして生成された判別用空間を用いて、被検体の組織を識別する。以下、その処理手順について図１５のフローを用いて説明する。
【００８４】
制御部６は、実施形態２のステップ８０４〜８０６と同様に、被検体１を水槽４に挿入するようにユーザに促し、被検体１が配置されたならば超音波を送信し、被検体１の透過波および被検体１による反射波の受信を行う（ステップ１５１〜１５３）。減衰画像生成部２４は、実施形態２のステップ８０７〜８０８の処理を行い、周波数ごとに減衰率画像を生成する（ステップ１５４）。音速画像生成部２６は、実施形態２のステップ８０９で説明した処理により音速画像を生成する（ステップ１５５）。反射画像生成部２５は、実施形態２のステップ８０９で説明した処理により反射画像を生成する（ステップ１５６）。
【００８５】
関心領域設定部２３は、実施形態２のステップ８０９と同様の処理により、反射画像または音速画像を用いて減衰画像をセグメントに分割し、それぞれを関心領域（ＲＯＩ）として設定する。２周波数間減衰率算出部２１は、実施形態２のステップ８０９と同様の処理により、減衰画像のＲＯＩとしたセグメント内の２周波数間減衰率を算出し、２周波数間減衰率画像を生成する（ステップ１５７）（図１６（ａ）参照）。
【００８６】
関心領域設定部２３は、音速画像にもステップ１５７の減衰画像に設定したＲＯＩと対応する位置にＲＯＩを設定する。組織種類識別部２２ａは、音速画像のＲＯＩ内の音速を求める（ステップ１５８）（図１６（ａ）参照）。
【００８７】
組織種類識別部２２ａは、ＲＯＩごとに求めた２周波数間減衰率と音速とを、判別用空間に図６（ａ）のようにプロットし、プロットした点が含まれる領域が示す組織の種類をＲＯＩの組織の種類であると識別する（ステップ１５９）（図１６（ｂ）参照）。例えば、点６３のＲＯＩの組織は、癌であり、点６４のＲＯＩの組織は、乳腺であると識別する。
【００８８】
悪性度算出部２２ｂは、悪性組織と正常組織との境界線６２と、ＲＯＩの点６３，６４との距離をそれぞれ算出し、距離の長さをそのＲＯＩの悪性度合とする（ステップ１６０）。
【００８９】
悪性度算出部２２ｂは、ＲＯＩ内の画素の画素値をステップ１６０で求めた悪性度合とする悪性度合画像を生成する（図１６（ｃ）参照）。
【００９０】
ユーザは、図１６（ｃ）の悪性度合画像と、図１６（ａ）の２周波数間減衰画像や音速画像を見ることにより、音速画像や反射画像に基づいて設定されたＲＯＩの位置とその組織の悪性度合を知ることができる。図１６（ａ）と図１６（ｃ）を比較すると明らかなように、ＲＯＩが設定された組織であっても、周囲の組織と悪性度合が同じ（正常である）組織は、悪性度合画像において同じ画素値となる。よって、ユーザは、悪性度合画像を見ることにより、悪性組織の位置と形状を明確に把握することができる。
【００９１】
また、図１６（ｂ）のようにステップ１５９で識別結果を参照することにより、ユーザは、ＲＯＩの組織の種類を認識することができる。
【００９２】
なお、実施形態３では、２周波数間減衰率と、他の物性値（音速）とを指標として組織の種類を識別したが、２周波数間減衰率を用いず、２以上の物性値を用いて判別空間を生成してもよい。２以上の物性値を用いて判別空間を生成することにより、組織の種類を識別し、悪性度合いを求めることが可能である。
【符号の説明】
【００９３】
１…被検体、２…ベッド、３…振動子、４…水槽、５…予備タンク、６…制御部、８…記憶部、９…入出力部、２０…信号処理部、２１…２周波数間減衰率算出部、２２…判別部、２２ａ…組織種類識別部、２２ｂ…悪性度算出部、２２ｃ…判別用空間生成部、２３…関心領域設定部、２４…減衰画像生成部、２５…反射画像生成部、２６…音速画像生成部、３０…振動子アレイ
【図１】
2020099619.tif 000003
【図２】
2020099619.tif 000004
【図３】
2020099619.tif 000005
【図４】
2020099619.tif 000006
【図５】
2020099619.tif 000007
【図６】
2020099619.tif 000008
【図７】
2020099619.tif 000009
【図８】
2020099619.tif 000010
【図９】
2020099619.tif 000011
【図１０】
2020099619.tif 000012
【図１１】
2020099619.tif 000013
【図１２】
2020099619.tif 000014
【図１３】
2020099619.tif 000015
【図１４】
2020099619.tif 000016
【図１５】
2020099619.tif 000017
【図１６】
2020099619.tif 000018