JP6157790B1 - 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム - Google Patents

Abstract

本発明にかかる超音波観測装置は、超音波信号に基づいて複数の特徴量を算出する特徴量算出部と、特徴量算出部が算出した複数の特徴量のうち、視覚情報と関連付けて超音波画像とともに表示する表示対象の特徴量を所定の表示仕様で配色した特徴量画像データを生成する特徴量画像データ生成部と、表示対象の特徴量以外の他の特徴量に基づいて、表示対象の特徴量の表示仕様を設定する表示仕様設定部と、を備えた。

Description

本発明は、超音波を用いて観測対象の組織を観測する超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラムに関する。

観測対象である生体組織または材料の特性を観測するために、超音波を適用することがある。具体的には、観測対象に超音波を送信し、その観測対象によって反射された超音波エコーに対して所定の信号処理を施すことにより、観測対象の特性に関する情報を取得する。

このうち、超音波を用いた被検体等の観測対象の組織性状を観測する技術として、受信した超音波信号の周波数スペクトルの特徴量を画像化する技術が知られている（例えば、特許文献１を参照）。この技術では、観測対象の組織性状を表す量として周波数スペクトルの特徴量を抽出した後、この特徴量に対応する視覚的な情報を付与した特徴量画像を生成して表示する。医師等のユーザは、表示された特徴量画像を見ることによって検体の組織性状を診断する。

例えば、特許文献１では、観測対象の組織の硬さを画像化した弾性画像が表示され、特徴量に応じて色を付与した色情報を表示する。この弾性画像は、一般にエラストグラフィと呼ばれ、設定された領域における観測対象の硬さに関する情報（特徴量）を取得し、超音波画像上に特徴量に応じた色情報を重畳する。具体的に、特許文献１では、画像内で選択された基準領域の歪みを基準として、計測対象の計測点における歪みと基準領域の歪みとの相対比に基づいて階調化した弾性画像を表示する。

特許第５１６０２２７号公報

ここで、特許文献１が開示する技術は、表示する特徴量について、基準領域の特徴量に対する相対的な色情報を付与することで硬さを表現するものである。このため、硬さを算出するために用いる特徴量によっては、色情報として表現できる範囲が狭くなって特徴量の差を明確に表現できない場合があった。

本発明は、上記に鑑みてなされたものであって、特徴量の差を明確に表現することができる超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラムを提供することを目的とする。

上述した課題を解決し、目的を達成するために、本発明に係る超音波観測装置は、観測対象へ超音波を送信し、該観測対象で反射された超音波を受信する超音波振動子を備えた超音波プローブが取得した超音波信号に基づいて超音波画像を生成する超音波観測装置であって、前記超音波信号に基づいて複数の特徴量を算出する特徴量算出部と、前記特徴量算出部が算出した前記複数の特徴量のうち、前記超音波画像とともに表示する表示対象の特徴量を所定の表示仕様で配色した特徴量画像データを生成する特徴量画像データ生成部と、前記表示対象の特徴量以外の他の特徴量に基づいて、前記表示対象の特徴量の表示仕様を設定する表示仕様設定部と、を備えたことを特徴とする。

本発明に係る超音波観測装置は、上記発明において、前記表示仕様設定部は、前記他の特徴量であって前記超音波信号に基づく画像内の対象領域の特徴量に基づいて、前記表示対象の特徴量の表示仕様を設定することを特徴とする。

本発明に係る超音波観測装置は、上記発明において、前記表示仕様設定部は、前記他の特徴量の代表値と閾値とを比較して、前記表示対象の特徴量の配色範囲、および／または色相の変化率を設定することを特徴とする。

本発明に係る超音波観測装置は、上記発明において、前記超音波信号に基づいて生成される信号の周波数を解析することによって複数の周波数スペクトルを算出する周波数解析部をさらに備え、前記特徴量算出部は、前記複数の特徴量の一つとして、前記周波数解析部が算出した周波数スペクトルに基づく周波数特徴量を算出することを特徴とする。

本発明に係る超音波観測装置は、上記発明において、前記表示対象の特徴量は、音速、硬さ、減衰量、前記周波数特徴量のうちのいずれかであり、前記他の特徴量は、音速、硬さ、減衰量、前記周波数特徴量のうちの前記表示対象の特徴量以外の少なくとも一つであることを特徴とする。

本発明に係る超音波観測装置は、上記発明において、前記特徴量算出部は、前記周波数解析部が算出した周波数スペクトルを直線近似することによって前記周波数特徴量を算出することを特徴とする。

本発明に係る超音波観測装置は、上記発明において、前記周波数特徴量は、前記周波数解析部が算出した周波数スペクトルを近似することによって得られた傾き、切片およびミッドバンドフィットのうち少なくとも一つであることを特徴とする。

本発明に係る超音波観測装置は、上記発明において、前記複数の特徴量のばらつきを算出するばらつき算出部と、前記ばらつき算出部が算出した各特徴量のばらつきをもとに、前記表示対象の特徴量を選定する特徴量選定部と、をさらに備え、前記表示仕様設定部は、前記特徴量選定部が選定した前記表示対象の特徴量以外の他の特徴量に基づいて、前記表示対象の特徴量の表示仕様を設定することを特徴とする。

本発明に係る超音波観測装置の作動方法は、観測対象へ超音波を送信し、該観測対象で反射された超音波を受信する超音波振動子を備えた超音波プローブが取得した超音波信号に基づいて超音波画像を生成する超音波観測装置の作動方法であって、特徴量算出部が、前記超音波信号に基づいて複数の特徴量を算出する特徴量算出ステップと、表示仕様設定部が、前記特徴量算出部が算出した前記複数の特徴量のうち、前記超音波画像とともに表示する表示対象の特徴量以外の他の特徴量に基づいて、前記表示対象の特徴量の表示仕様を設定する表示仕様設定ステップと、特徴量画像データ生成部が、前記表示対象の特徴量を所定の表示仕様で配色した特徴量画像データを生成する特徴量画像データ生成ステップと、を含むことを特徴とする。

本発明に係る超音波観測装置の作動プログラムは、観測対象へ超音波を送信し、該観測対象で反射された超音波を受信する超音波振動子を備えた超音波プローブが取得した超音波信号に基づいて超音波画像を生成する超音波観測装置の作動プログラムであって、特徴量算出部が、前記超音波信号に基づいて複数の特徴量を算出する特徴量算出手順と、表示仕様設定部が、前記特徴量算出部が算出した前記複数の特徴量のうち、前記超音波画像とともに表示する表示対象の特徴量以外の他の特徴量に基づいて、前記表示対象の特徴量の表示仕様を設定する表示仕様設定手順と、特徴量画像データ生成部が、前記表示対象の特徴量を所定の表示仕様で配色した特徴量画像データを生成する特徴量画像データ生成手順と、を前記超音波観測装置に実行させることを特徴とする。

本発明によれば、特徴量の差を明確に表現することができるという効果を奏する。

図１は、本発明の実施の形態１に係る超音波観測装置を備えた超音波観測システムの構成を示すブロック図である。 図２は、本発明の実施の形態１に係る超音波観測装置の信号増幅部が行う増幅処理における受信深度と増幅率との関係を示す図である。 図３は、本発明の実施の形態１に係る超音波観測装置の増幅補正部が行う増幅補正処理における受信深度と増幅率との関係を示す図である。 図４は、超音波信号の１つの音線におけるデータ配列を模式的に示す図である。 図５は、本発明の実施の形態１に係る超音波観測装置の周波数解析部により算出された周波数スペクトルの例を示す図である。 図６は、本発明の実施の形態１に係る超音波観測装置の減衰補正部が補正した補正特徴量をパラメータとして有する直線を示す図である。 図７は、本発明の実施の形態１に係る超音波観測装置の表示仕様設定部が実行する処理を説明する図である。 図８は、本発明の実施の形態１に係る超音波観測装置が行う処理の概要を示すフローチャートである。 図９は、本発明の実施の形態１に係る超音波観測装置の周波数解析部が実行する処理の概要を示すフローチャートである。 図１０は、本発明の実施の形態１に係る超音波観測装置の表示装置における特徴量画像の表示例を模式的に示す図である。 図１１は、本発明の実施の形態１の変形例１に係る超音波観測装置の表示仕様設定部が実行する処理を説明する図である。 図１２は、本発明の実施の形態１の変形例２に係る超音波観測装置の表示仕様設定部が実行する処理を説明する図である。 図１３は、本発明の実施の形態１の変形例３に係る超音波観測装置の表示仕様設定部が実行する処理を説明する図である。 図１４は、本発明の実施の形態２に係る超音波観測装置を備えた超音波観測システムの構成を示すブロック図である。 図１５は、本発明の実施の形態３に係る超音波観測装置を備えた超音波観測システムの構成を示すブロック図である。 図１６は、本発明の実施の形態３に係る超音波観測装置の表示仕様設定部が実行する処理を説明する図である。

以下、添付図面を参照して、本発明を実施するための形態（以下、「実施の形態」という）を説明する。

（実施の形態１）
図１は、本発明の実施の形態１に係る超音波観測装置３を備えた超音波観測システム１の構成を示すブロック図である。同図に示す超音波観測システム１は、観測対象である被検体へ超音波を送信し、該被検体で反射された超音波を受信する超音波内視鏡２（超音波プローブ）と、超音波内視鏡２が取得した超音波信号に基づいて超音波画像を生成する超音波観測装置３と、超音波観測装置３が生成した超音波画像を表示する表示装置４と、を備える。

超音波内視鏡２は、その先端部に、超音波観測装置３から受信した電気的なパルス信号を超音波パルス（音響パルス）に変換して被検体へ照射するとともに、被検体で反射された超音波エコーを電圧変化で表現する電気的なエコー信号（超音波信号）に変換して出力する超音波振動子２１を有する。超音波振動子２１は、コンベックス振動子、リニア振動子およびラジアル振動子のいずれでも構わない。超音波内視鏡２は、超音波振動子２１をメカ的に走査させるものであってもよいし、超音波振動子２１として複数の素子をアレイ状に設け、送受信にかかわる素子を電子的に切り替えたり、各素子の送受信に遅延をかけたりすることで、電子的に走査させるものであってもよい。

超音波内視鏡２は、通常は撮像光学系および撮像素子を有しており、被検体の消化管（食道、胃、十二指腸、大腸）、または呼吸器（気管・気管支）へ挿入され、消化管や呼吸器、その周囲臓器（膵臓、胆嚢、胆管、胆道、リンパ節、縦隔臓器、血管等）を撮像することが可能である。また、超音波内視鏡２は、撮像時に被検体へ照射する照明光を導くライトガイドを有する。このライトガイドは、先端部が超音波内視鏡２の被検体への挿入部の先端まで達している一方、基端部が照明光を発生する光源装置に接続されている。なお、超音波内視鏡２に限らず、撮像光学系および撮像素子を有しない超音波プローブであってもよい。

超音波観測装置３は、超音波内視鏡２と電気的に接続され、所定の波形および送信タイミングに基づいて高電圧パルスからなる送信信号（パルス信号）を超音波振動子２１へ送信するとともに、超音波振動子２１から電気的な受信信号であるエコー信号を受信してデジタルの高周波（ＲＦ：Radio Frequency）信号のデータ（以下、ＲＦデータという）を生成、出力する送受信部３１と、送受信部３１から受信したＲＦデータをもとにデジタルのＢモード用受信データを生成する信号処理部３２と、送受信部３１から受信したＲＦデータに対して所定の演算を施す演算部３３と、各種画像データを生成する画像処理部３４と、キーボード、マウス、タッチパネル等のユーザインタフェースを用いて実現され、各種情報の入力を受け付ける入力部３５と、超音波観測システム１全体を制御する制御部３６と、超音波観測装置３の動作に必要な各種情報を記憶する記憶部３７と、を備える。

送受信部３１は、エコー信号を増幅する信号増幅部３１１を有する。信号増幅部３１１は、受信深度が大きいエコー信号ほど高い増幅率で増幅するＳＴＣ（Sensitivity Time Control）補正を行う。図２は、信号増幅部３１１が行う増幅処理における受信深度と増幅率との関係を示す図である。図２に示す受信深度ｚは、超音波の受信開始時点からの経過時間に基づいて算出される量である。図２に示すように、増幅率β（ｄＢ）は、受信深度ｚが閾値ｚ_thより小さい場合、受信深度ｚの増加に伴ってβ_０からβ_th（＞β₀）へ線型に増加する。また、増幅率β（ｄＢ）は、受信深度ｚが閾値ｚ_th以上である場合、一定値β_thをとる。閾値ｚ_thの値は、観測対象から受信する超音波信号がほとんど減衰してしまい、ノイズが支配的になるような値である。より一般に、増幅率βは、受信深度ｚが閾値ｚ_thより小さい場合、受信深度ｚの増加に伴って単調増加すればよい。なお、図２に示す関係は、予め記憶部３７に記憶されている。

送受信部３１は、信号増幅部３１１によって増幅されたエコー信号に対してフィルタリング等の処理を施した後、Ａ／Ｄ変換することによって時間ドメインのＲＦデータを生成し、信号処理部３２および演算部３３へ出力する。なお、超音波内視鏡２が複数の素子をアレイ状に設けた超音波振動子２１を電子的に走査させる構成を有する場合、送受信部３１は、複数の素子に対応したビーム合成用の多チャンネル回路を有する。

送受信部３１が送信するパルス信号の周波数帯域は、超音波振動子２１におけるパルス信号の超音波パルスへの電気音響変換の線型応答周波数帯域をほぼカバーする広帯域にするとよい。また、信号増幅部３１１におけるエコー信号の各種処理周波数帯域は、超音波振動子２１による超音波エコーのエコー信号への音響電気変換の線型応答周波数帯域をほぼカバーする広帯域にするとよい。これらにより、後述する周波数スペクトルの近似処理を実行する際、精度のよい近似を行うことが可能となる。

送受信部３１は、制御部３６が出力する各種制御信号を超音波内視鏡２に対して送信するとともに、超音波内視鏡２から識別用のＩＤを含む各種情報を受信して制御部３６へ送信する機能も有する。

信号処理部３２は、ＲＦデータに対してバンドパスフィルタ、包絡線検波、対数変換など公知の処理を施し、デジタルのＢモード用受信データを生成する。対数変換では、ＲＦデータを基準電圧Ｖ_cで除した量の常用対数をとってデシベル値で表現する。信号処理部３２は、生成したＢモード用受信データを、画像処理部３４へ出力する。信号処理部３２は、ＣＰＵ（Central Processing Unit）や各種演算回路等を用いて実現される。

演算部３３は、送受信部３１が生成したＲＦデータに対して受信深度ｚによらず増幅率β（ｄＢ）を一定とするよう増幅補正を行う増幅補正部３３１と、増幅補正を行ったＲＦデータに高速フーリエ変換（ＦＦＴ：Fast Fourier Transform）を施して周波数解析を行うことにより周波数スペクトルを算出する周波数解析部３３２と、周波数解析部３３２により算出された周波数スペクトルをもとに、該周波数スペクトルの特徴量を算出する特徴量算出部３３３と、表示装置４に表示させる表示対象の特徴量とは異なる特徴量に基づいて、表示対象の特徴量の表示仕様を設定する表示仕様設定部３３４と、を有する。演算部３３は、ＣＰＵや各種演算回路等を用いて実現される。

図３は、増幅補正部３３１が行う増幅補正処理における受信深度と増幅率との関係を示す図である。図３に示すように、増幅補正部３３１が行う増幅補正処理における増幅率βは、受信深度ｚがゼロのとき最大値β_th−β₀をとり、受信深度ｚがゼロから閾値ｚ_thに達するまで線型に減少し、受信深度ｚが閾値ｚ_th以上のときゼロである。このように定められる増幅率βによって増幅補正部３３１がデジタルＲＦ信号を増幅補正することにより、信号処理部３２におけるＳＴＣ補正の影響を相殺し、一定の増幅率β_thの信号を出力することができる。なお、増幅補正部３３１が行う受信深度ｚと増幅率βの関係は、信号処理部３２における受信深度と増幅率の関係に応じて異なることは勿論である。

このような増幅補正を行う理由を説明する。ＳＴＣ補正は、アナログ信号波形の振幅を全周波数帯域にわたって均一に、かつ、深度に対しては単調増加する増幅率で増幅させることで、アナログ信号波形の振幅から減衰の影響を排除する補正処理である。このため、エコー信号の振幅を輝度に変換して表示するＢモード画像を生成する場合、かつ、一様な組織を走査した場合には、ＳＴＣ補正を行うことによって深度によらず輝度値が一定になる。すなわち、Ｂモード画像の輝度値から減衰の影響を排除する効果を得ることができる。

一方、本実施の形態のように超音波の周波数スペクトルを算出して解析した結果を利用する場合、ＳＴＣ補正でも超音波の伝播に伴う減衰の影響を正確に排除できるわけではない。なぜなら、一般に減衰量は周波数によって異なるが（後述する式（１）を参照）、ＳＴＣ補正の増幅率は距離だけに応じて変化し、周波数依存性がないためである。

上述した問題、すなわち、超音波の周波数スペクトルを算出して解析した結果を利用する場合、ＳＴＣ補正でも超音波の伝播に伴う減衰の影響を正確に排除できるわけではない、という問題を解決するには、Ｂモード画像を生成する際にＳＴＣ補正を施した受信信号を出力する一方、周波数スペクトルに基づいた画像を生成する際に、Ｂモード画像を生成するための送信とは異なる新たな送信を行い、ＳＴＣ補正を施していない受信信号を出力することが考えられる。ところがこの場合には、受信信号に基づいて生成される画像データのフレームレートが低下してしまうという問題がある。

そこで、本実施の形態では、生成される画像データのフレームレートを維持しつつ、Ｂモード画像用にＳＴＣ補正を施した信号に対してＳＴＣ補正の影響を排除するために、増幅補正部３３１によって増幅率の補正を行う。

周波数解析部３３２は、増幅補正部３３１が増幅補正した各音線のＲＦデータ（ラインデータ）を所定の時間間隔でサンプリングし、サンプルデータを生成する。周波数解析部３３２は、サンプルデータ群にＦＦＴ処理を施すことにより、ＲＦデータ上の複数の箇所（データ位置）における周波数スペクトルを算出する。ここでいう「周波数スペクトル」とは、サンプルデータ群にＦＦＴ処理を施すことによって得られた「ある受信深度ｚにおける強度の周波数分布」を意味する。また、ここでいう「強度」とは、例えばエコー信号の電圧、エコー信号の電力、超音波エコーの音圧、超音波エコーの音響エネルギー等のパラメータ、これらパラメータの振幅や時間積分値やその組み合わせのいずれかを指す。

一般に、周波数スペクトルは、観測対象が生体組織である場合、超音波が走査された生体組織の性状によって異なる傾向を示す。これは、周波数スペクトルが、超音波を散乱する散乱体の大きさ、数密度、音響インピーダンス等と相関を有しているためである。ここでいう「生体組織の性状」とは、例えば悪性腫瘍（癌）、良性腫瘍、内分泌腫瘍、粘液性腫瘍、正常組織、嚢胞、脈管などのことである。

図４は、超音波信号の１つの音線におけるデータ配列を模式的に示す図である。同図に示す音線ＳＲ_kにおいて、白または黒の長方形は、１つのサンプル点におけるデータを意味している。また、音線ＳＲ_kにおいて、右側に位置するデータほど、超音波振動子２１から音線ＳＲ_kに沿って計った場合の深い箇所からのサンプルデータである（図４の矢印を参照）。音線ＳＲ_kは、送受信部３１が行うＡ／Ｄ変換におけるサンプリング周波数（例えば５０ＭＨｚ）に対応した時間間隔で離散化されている。図４では、番号ｋの音線ＳＲ_kの８番目のデータ位置を受信深度ｚの方向の初期値Ｚ^(k) ₀として設定した場合を示しているが、初期値の位置は任意に設定することができる。周波数解析部３３２による算出結果は複素数で得られ、記憶部３７に格納される。

図４に示すデータ群Ｆ_j（ｊ＝１、２、・・・、Ｋ）は、ＦＦＴ処理の対象となるサンプルデータ群である。一般に、ＦＦＴ処理を行うためには、サンプルデータ群が２のべき乗のデータ数を有している必要がある。この意味で、サンプルデータ群Ｆ_j（ｊ＝１、２、・・・、Ｋ−１）はデータ数が１６（＝２⁴）で正常なデータ群である一方、サンプルデータ群Ｆ_Kは、データ数が１２であるため異常なデータ群である。異常なデータ群に対してＦＦＴ処理を行う際には、不足分だけゼロデータを挿入することにより、正常なサンプルデータ群を生成する処理を行う。この点については、周波数解析部３３２の処理を説明する際に詳述する（図９を参照）。

図５は、周波数解析部３３２により算出された周波数スペクトルの例を示す図である。図５では、横軸が周波数ｆである。また、図５では、縦軸が、強度Ｉ₀を基準強度Ｉ_c（定数）で除した量の常用対数（デシベル表現）Ｉ＝１０ｌｏｇ₁₀（Ｉ₀／Ｉ_c）である。図５に示す直線Ｌ₁₀については後述する。なお、本実施の形態において、曲線および直線は、離散的な点の集合からなる。

図５に示す周波数スペクトルＣ₁において、以後の演算に使用する周波数帯域の下限周波数ｆ_Lおよび上限周波数ｆ_Hは、超音波振動子２１の周波数帯域、送受信部３１が送信するパルス信号の周波数帯域などをもとに決定されるパラメータである。以下、図５において、下限周波数ｆ_Lおよび上限周波数ｆ_Hによって定まる周波数帯域を「周波数帯域Ｆ」という。

特徴量算出部３３３は、設定されている関心領域（ＲＯＩ）内において、複数の周波数スペクトルの特徴量をそれぞれ算出する。特徴量算出部３３３は、周波数スペクトルを直線で近似することによって減衰補正処理を行う前の周波数スペクトルの特徴量（以下、補正前特徴量という）を算出する近似部３３３ａと、近似部３３３ａが算出した補正前特徴量に対して減衰補正を行うことによって特徴量を算出する減衰補正部３３３ｂと、を有する。

近似部３３３ａは、所定周波数帯域における周波数スペクトルの回帰分析を行って周波数スペクトルを一次式（回帰直線）で近似することにより、この近似した一次式を特徴付ける補正前特徴量を算出する。例えば、図５に示す周波数スペクトルＣ₁の場合、近似部３３３ａは、周波数帯域Ｆで回帰分析を行い周波数スペクトルＣ₁を一次式で近似することによって回帰直線Ｌ₁₀を得る。換言すると、近似部３３３ａは、回帰直線Ｌ₁₀の傾きａ₀、切片ｂ₀、および周波数帯域Ｆの中心周波数ｆ_M＝（ｆ_L＋ｆ_H）／２の回帰直線上の値であるミッドバンドフィット（Mid-band fit）ｃ₀＝ａ₀ｆ_M＋ｂ₀を補正前特徴量として算出する。

３つの補正前特徴量のうち、傾きａ₀は、超音波の散乱体の大きさと相関を有し、一般に散乱体が大きいほど傾きが小さな値を有すると考えられる。また、切片ｂ₀は、散乱体の大きさ、音響インピーダンスの差、散乱体の数密度（濃度）等と相関を有している。具体的には、切片ｂ₀は、散乱体が大きいほど大きな値を有し、音響インピーダンスの差が大きいほど大きな値を有し、散乱体の数密度が大きいほど大きな値を有すると考えられる。ミッドバンドフィットｃ₀は、傾きａ₀と切片ｂ₀から導出される間接的なパラメータであり、有効な周波数帯域内の中心におけるスペクトルの強度を与える。このため、ミッドバンドフィットｃ₀は、散乱体の大きさ、音響インピーダンスの差、散乱体の数密度に加えて、Ｂモード画像の輝度とある程度の相関を有していると考えられる。なお、特徴量算出部３３３は、回帰分析によって二次以上の多項式で周波数スペクトルを近似するようにしてもよい。

減衰補正部３３３ｂが行う補正について説明する。一般に、超音波の減衰量Ａ（ｆ，ｚ）は、超音波が受信深度０と受信深度ｚとの間を往復する間に生じる減衰であり、往復する前後の強度変化（デシベル表現での差）として定義される。減衰量Ａ（ｆ，ｚ）は、一様な組織内では周波数に比例することが経験的に知られており、以下の式（１）で表現される。
Ａ（ｆ，ｚ）＝２αｚｆ ・・・（１）
ここで、比例定数αは減衰率と呼ばれる量である。また、ｚは超音波の受信深度であり、ｆは周波数である。減衰率αの具体的な値は、観測対象が生体である場合、生体の部位に応じて定まる。減衰率αの単位は、例えばｄＢ／ｃｍ／ＭＨｚである。なお、本実施の形態において、減衰率αの値を入力部３５からの入力によって変更できる構成とすることも可能である。

減衰補正部３３３ｂは、近似部３３３ａが抽出した補正前特徴量（傾きａ₀、切片ｂ₀、ミッドバンドフィットｃ₀）に対し、以下に示す式（２）〜（４）にしたがって減衰補正を行うことにより、特徴量ａ、ｂ、ｃを算出する。
ａ＝ａ₀＋２αｚ ・・・（２）
ｂ＝ｂ₀ ・・・（３）
ｃ＝ｃ₀＋Ａ（ｆ_M，ｚ）＝ｃ₀＋２αｚｆ_M（＝ａｆ_M＋ｂ） ・・・（４）
式（２）、（４）からも明らかなように、減衰補正部３３３ｂは、超音波の受信深度ｚが大きいほど、補正量が大きい補正を行う。また、式（３）によれば、切片に関する補正は恒等変換である。これは、切片が周波数０（Ｈｚ）に対応する周波数成分であって減衰の影響を受けないためである。

図６は、減衰補正部３３３ｂが算出した特徴量ａ、ｂ、ｃをパラメータとして有する直線を示す図である。直線Ｌ₁の式は、
Ｉ＝ａｆ＋ｂ＝（ａ₀＋２αｚ）ｆ＋ｂ₀ ・・・（５）
で表される。この式（５）からも明らかなように、直線Ｌ₁は、減衰補正前の直線Ｌ₁₀と比較して、傾きが大きく（ａ＞ａ₀）、かつ切片が同じ（ｂ＝ｂ₀）である。

表示仕様設定部３３４は、表示装置４に表示させる表示対象の特徴量とは異なる特徴量に基づいて、表示対象の特徴量の表示仕様を設定する。具体的に、表示仕様設定部３３４は、本実施の形態１では、傾きを示す特徴量ａの表示仕様である色相スケールを、ミッドバンドフィットを示す特徴量ｃに基づいて設定する。本明細書でいう色相スケールとは、光の波長が異なる色（色相）を連続的（多段階的を含む）に並べてなり、色相と特徴量の値とを対応付けたものである。図７では、左にいくほど赤色に近づき、右にいくほど青色に近づく。具体的には、左から、可視光線の波長が長い順に、赤、橙、黄、緑、青となっている。例えば、最も短い波長が３８０ｎｍであり、最も長い波長が７５０ｎｍである。図７は、本発明の実施の形態１に係る超音波観測装置３の表示仕様設定部３３４が実行する処理を説明する図である。図７では、横軸が特徴量ａである。また、図７では、縦軸が特徴量ｃである。図７は、特徴量ａおよび特徴量ｃの分布を示すグラフとなっている。

ここで、図７に示すように、生体組織の性状（以下、組織性状という）の種類によって、特徴量ａおよび特徴量ｃの分布が異なる。例えば、ある組織性状は領域Ｑ₁に分布され、それとは異なる組織性状は領域Ｑ₂に分布される。この際、同じ色相スケールで特徴量ａを表現すると、一方の組織性状は特徴量ａの値に応じて明確に表現されるものの、他方の組織性状は特徴量ａの値の差が明確に表現されなかったり、色相スケールから外れた領域に特徴量ａの値が分布されて特徴量ａが表示されなかったりしてしまう。

本実施の形態１では、表示仕様設定部３３４が、特徴量ｃに基づいて、表示対象の特徴量ａの表示仕様を設定する。具体的に、表示仕様設定部３３４は、特徴量算出部３３３が算出した複数の特徴量ｃの平均値、例えば、設定されている関心領域における特徴量の平均値を求め、該平均値と閾値ｃ_thとを比較して、比較結果に基づいて色相スケールを設定する。例えば、表示仕様設定部３３４は、特徴量ｃの平均値が閾値ｃ_thより小さい場合、色相スケールを色相スケールＣＢ₁に設定する。これに対し、表示仕様設定部３３４は、特徴量ｃの平均値が閾値ｃ_th以上である場合、色相スケールＣＢ₁とは配色範囲が異なる色相スケールＣＢ₂に設定する。このようにして、特徴量ａの色相スケールを特徴量ｃに基づいて設定することによって、特徴量ａおよび特徴量ｃの分布に応じた組織性状に対応する色相スケールを設定することができる。

画像処理部３４は、エコー信号の振幅を輝度に変換して表示する超音波画像であるＢモード画像データを生成するＢモード画像データ生成部３４１と、減衰補正部３３３ｂが算出した特徴量を視覚情報と関連付けてＢモード画像とともに表示する特徴量画像データを生成する特徴量画像データ生成部３４２と、を有する。

Ｂモード画像データ生成部３４１は、信号処理部３２から受信したＢモード用受信データに対してゲイン処理、コントラスト処理等の公知の技術を用いた信号処理を行うとともに、表示装置４における画像の表示レンジに応じて定まるデータステップ幅に応じたデータの間引き等を行うことによってＢモード画像データを生成する。Ｂモード画像は、色空間としてＲＧＢ表色系を採用した場合の変数であるＲ（赤）、Ｇ（緑）、Ｂ（青）の値を一致させたグレースケール画像である。

Ｂモード画像データ生成部３４１は、信号処理部３２からのＢモード用受信データに対して走査範囲を空間的に正しく表現できるよう並べ直す座標変換を施した後、Ｂモード用受信データ間の補間処理を施すことによってＢモード用受信データ間の空隙を埋め、Ｂモード画像データを生成する。Ｂモード画像データ生成部３４１は、生成したＢモード画像データを特徴量画像データ生成部３４２へ出力する。

特徴量画像データ生成部３４２は、特徴量算出部３３３が算出した特徴量に関連する視覚情報をＢモード画像データにおける画像の各画素に対して重畳することによって特徴量画像データを生成する。特徴量画像データ生成部３４２は、例えば図４に示す１つのサンプルデータ群Ｆ_j（ｊ＝１、２、・・・、Ｋ）のデータ量に対応する画素領域に対し、そのサンプルデータ群Ｆ_jから算出される周波数スペクトルの特徴量に対応する視覚情報を割り当てる。特徴量画像データ生成部３４２は、例えば上述した傾き、切片、ミッドバンドフィットのいずれか一つに視覚情報としての色相を対応付けることによって特徴量画像データを生成する。具体的に、特徴量画像データ生成部３４２は、特徴量ａに視覚情報としての色相を対応付ける場合、表示仕様設定部３３４が設定した色相スケールを参照して視覚情報を割り当てる。特徴量に関連する視覚情報としては、色相のほか、例えば彩度、明度、輝度値、Ｒ（赤）、Ｇ（緑）、Ｂ（青）などの所定の表色系を構成する色空間の変数を挙げることができる。

制御部３６は、演算および制御機能を有するＣＰＵ（Central Processing Unit）や各種演算回路等を用いて実現される。制御部３６は、記憶部３７が記憶、格納する情報を記憶部３７から読み出し、超音波観測装置３の作動方法に関連した各種演算処理を実行することによって超音波観測装置３を統括して制御する。なお、制御部３６を信号処理部３２および演算部３３と共通のＣＰＵ等を用いて構成することも可能である。

記憶部３７は、減衰補正部３３３ｂが周波数スペクトルごとに算出した複数の特徴量や、画像処理部３４が生成した画像データを記憶する。また、記憶部３７は、図７に示したような、複数の特徴量のうちの二つの特徴量の分布と、該分布における色相スケールの割り当て設定を行うための閾値と、複数の色相スケールとを特徴量の組み合わせごとに記憶する表示仕様情報記憶部３７１を有する。

記憶部３７は、上記以外にも、例えば増幅処理に必要な情報（図２に示す増幅率と受信深度との関係）、増幅補正処理に必要な情報（図３に示す増幅率と受信深度との関係）、減衰補正処理に必要な情報（式（１）参照）、周波数解析処理に必要な窓関数（Hamming、Hanning、Blackman等）の情報等を記憶する。

また、記憶部３７は、超音波観測装置３の作動方法を実行するための作動プログラムを含む各種プログラムを記憶する。作動プログラムは、ハードディスク、フラッシュメモリ、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、フレキシブルディスク等のコンピュータ読み取り可能な記録媒体に記録して広く流通させることも可能である。なお、上述した各種プログラムは、通信ネットワークを介してダウンロードすることによって取得することも可能である。ここでいう通信ネットワークは、例えば既存の公衆回線網、ＬＡＮ（Local Area Network）、ＷＡＮ（Wide Area Network）などによって実現されるものであり、有線、無線を問わない。

以上の構成を有する記憶部３７は、各種プログラム等が予めインストールされたＲＯＭ（Read Only Memory）、および各処理の演算パラメータやデータ等を記憶するＲＡＭ（Random Access Memory）等を用いて実現される。

図８は、以上の構成を有する超音波観測装置３が行う処理の概要を示すフローチャートである。まず、超音波観測装置３は、超音波内視鏡２から超音波振動子２１による観測対象の測定結果としてのエコー信号を受信する（ステップＳ１）。

超音波振動子２１からエコー信号を受信した信号増幅部３１１は、そのエコー信号の増幅を行う（ステップＳ２）。ここで、信号増幅部３１１は、例えば図２に示す増幅率と受信深度との関係に基づいてエコー信号の増幅（ＳＴＣ補正）を行う。

続いて、Ｂモード画像データ生成部３４１は、信号増幅部３１１が増幅したエコー信号を用いてＢモード画像データを生成して、表示装置４へ出力する（ステップＳ３）。Ｂモード画像データを受信した表示装置４は、そのＢモード画像データに対応するＢモード画像を表示する（ステップＳ４）。

増幅補正部３３１は、送受信部３１から出力された信号に対して受信深度によらず増幅率が一定となる増幅補正を行う（ステップＳ５）。ここで、増幅補正部３３１は、例えば図３に示す増幅率と受信深度との関係が成立するように増幅補正を行う。

この後、周波数解析部３３２は、ＦＦＴ処理による周波数解析を行うことによって全てのサンプルデータ群に対する周波数スペクトルを算出する（ステップＳ６：周波数解析ステップ）。図９は、ステップＳ６において周波数解析部３３２が実行する処理の概要を示すフローチャートである。以下、図９に示すフローチャートを参照して、周波数解析処理を詳細に説明する。

まず、周波数解析部３３２は、解析対象の音線を識別するカウンタｋをｋ₀とする（ステップＳ２１）。

続いて、周波数解析部３３２は、ＦＦＴ処理用に取得する一連のデータ群（サンプルデータ群）を代表するデータ位置（受信深度に相当）Ｚ^(k)の初期値Ｚ^(k) ₀を設定する（ステップＳ２２）。例えば、図４では、上述したように、音線ＳＲ_kの８番目のデータ位置を初期値Ｚ^(k) ₀として設定した場合を示している。

その後、周波数解析部３３２は、サンプルデータ群を取得し（ステップＳ２３）、取得したサンプルデータ群に対し、記憶部３７が記憶する窓関数を作用させる（ステップＳ２４）。このようにサンプルデータ群に対して窓関数を作用させることにより、サンプルデータ群が境界で不連続になることを回避し、アーチファクトが発生するのを防止することができる。

続いて、周波数解析部３３２は、データ位置Ｚ^(k)のサンプルデータ群が正常なデータ群であるか否かを判定する（ステップＳ２５）。図４を参照した際に説明したように、サンプルデータ群は、２のべき乗のデータ数を有している必要がある。以下、正常なサンプルデータ群のデータ数を２ⁿ（ｎは正の整数）とする。本実施の形態では、データ位置Ｚ^(k)が、できるだけＺ^(k)が属するサンプルデータ群の中心になるよう設定される。具体的には、サンプルデータ群のデータ数は２ⁿであるので、Ｚ^(k)はそのサンプルデータ群の中心に近い２ⁿ／２（＝２^n-1）番目の位置に設定される。この場合、サンプルデータ群が正常であるとは、データ位置Ｚ^(k)の前方に２^n-1−１（＝Ｎとする）個のデータがあり、データ位置Ｚ^(k)の後方に２^n-1（＝Ｍとする）個のデータがあることを意味する。図４に示す場合、サンプルデータ群Ｆ₁、Ｆ₂、Ｆ₃、・・・、Ｆ_K-1はともに正常である。なお、図４ではｎ＝４（Ｎ＝７，Ｍ＝８）の場合を例示している。

ステップＳ２５における判定の結果、データ位置Ｚ^(k)のサンプルデータ群が正常である場合（ステップＳ２５：Ｙｅｓ）、周波数解析部３３２は、後述するステップＳ２７へ移行する。

ステップＳ２５における判定の結果、データ位置Ｚ^(k)のサンプルデータ群が正常でない場合（ステップＳ２５：Ｎｏ）、周波数解析部３３２は、不足分だけゼロデータを挿入することによって正常なサンプルデータ群を生成する（ステップＳ２６）。ステップＳ２５において正常でないと判定されたサンプルデータ群（例えば図４のサンプルデータ群Ｆ_K）は、ゼロデータを追加する前に窓関数が作用されている。このため、サンプルデータ群にゼロデータを挿入してもデータの不連続は生じない。ステップＳ２６の後、周波数解析部３３２は、後述するステップＳ２７へ移行する。

ステップＳ２７において、周波数解析部３３２は、サンプルデータ群を用いてＦＦＴ処理を行うことにより、振幅の周波数分布である周波数スペクトルを得る（ステップＳ２７）。

続いて、周波数解析部３３２は、データ位置Ｚ^(k)をステップ幅Ｄで変化させる（ステップＳ２８）。ステップ幅Ｄは、記憶部３７が予め記憶しているものとする。図４では、Ｄ＝１５の場合を例示している。ステップ幅Ｄは、Ｂモード画像データ生成部３４１がＢモード画像データを生成する際に利用するデータステップ幅と一致させることが望ましいが、周波数解析部３３２における演算量を削減したい場合には、ステップ幅Ｄとしてデータステップ幅より大きい値を設定してもよい。

その後、周波数解析部３３２は、データ位置Ｚ^(k)が音線ＳＲ_kにおける最大値Ｚ^(k) _maxより大きいか否かを判定する（ステップＳ２９）。データ位置Ｚ^(k)が最大値Ｚ^(k) _maxより大きい場合（ステップＳ２９：Ｙｅｓ）、周波数解析部３３２はカウンタｋを１増加させる（ステップＳ３０）。これは、処理をとなりの音線へ移すことを意味する。一方、データ位置Ｚ^(k)が最大値Ｚ^(k) _max以下である場合（ステップＳ２９：Ｎｏ）、周波数解析部３３２はステップＳ２３へ戻る。このようにして、周波数解析部３３２は、音線ＳＲ_kに対して、［（Ｚ^(k) _max−Ｚ^(k) ₀＋１）／Ｄ＋１］個のサンプルデータ群に対するＦＦＴ処理を行う。ここで、［Ｘ］は、Ｘを超えない最大の整数を表す。

ステップＳ３０の後、周波数解析部３３２は、カウンタｋが最大値ｋ_maxより大きいか否かを判定する（ステップＳ３１）。カウンタｋが最大値ｋ_maxより大きい場合（ステップＳ３１：Ｙｅｓ）、周波数解析部３３２は一連の周波数解析処理を終了する。一方、カウンタｋが最大値ｋ_max以下である場合（ステップＳ３１：Ｎｏ）、周波数解析部３３２はステップＳ２２に戻る。この最大値ｋ_maxは、医師等のユーザが入力部３５を通じて任意に指示入力した値、もしくは、記憶部３７にあらかじめ設定された値とする。

このようにして、周波数解析部３３２は、解析対象領域内の（ｋ_max−ｋ₀＋１）本の音線の各々について複数回のＦＦＴ処理を行う。ＦＦＴ処理の結果は、受信深度および受信方向とともに記憶部３７に格納される。

なお、以上の説明では、周波数解析部３３２が超音波信号を受信したすべての領域に対して周波数解析処理を行うものとしたが、設定された関心領域内においてのみ周波数解析処理を行うようにすることも可能である。

以上説明したステップＳ６の周波数解析処理に続いて、特徴量算出部３３３は、複数の周波数スペクトルの補正前特徴量をそれぞれ算出し、各周波数スペクトルの補正前特徴量に対して超音波の減衰の影響を排除する減衰補正を行うことによって各周波数スペクトルの補正特徴量を算出する（ステップＳ７〜Ｓ８：特徴量算出ステップ）。

ステップＳ７において、近似部３３３ａは、周波数解析部３３２が生成した複数の周波数スペクトルをそれぞれ回帰分析することにより、各周波数スペクトルに対応する補正前特徴量を算出する（ステップＳ７）。具体的には、近似部３３３ａは、各周波数スペクトルを回帰分析することによって一次式で近似し、補正前特徴量として傾きａ₀、切片ｂ₀、ミッドバンドフィットｃ₀を算出する。例えば、図５に示す直線Ｌ₁₀は、近似部３３３ａが周波数帯域Ｆの周波数スペクトルＣ₁に対し回帰分析によって近似した回帰直線である。

続いて、減衰補正部３３３ｂは、近似部３３３ａが各周波数スペクトルに対して近似した補正前特徴量に対し、減衰率αを用いて減衰補正を行うことにより、補正特徴量を算出し、算出した補正特徴量を記憶部３７に格納する（ステップＳ８）。図６に示す直線Ｌ₁は、減衰補正部３３３ｂが減衰補正処理を行うことによって得られる直線の例である。

ステップＳ８において、減衰補正部３３３ｂは、上述した式（２）、（４）における受信深度ｚに、超音波信号の音線のデータ配列を用いて得られるデータ位置Ｚ＝（ｆ_sp／２ｖ_s）Ｄｎを代入することによって算出する。ここで、ｆ_spはデータのサンプリング周波数、ｖ_sは音速、Ｄはデータステップ幅、ｎは処理対象のサンプルデータ群のデータ位置までの音線の１番目のデータからのデータステップ数である。例えば、データのサンプリング周波数ｆ_spを５０ＭＨｚとし、音速ｖ_sを１５３０ｍ／ｓｅｃとし、図４に示すデータ配列を採用してステップ幅Ｄを１５とすると、ｚ＝０．２２９５ｎ（ｍｍ）となる。

その後、Ｂモード画像データ生成部３４１が生成したＢモード画像データにおける各画素に対し、ステップＳ８で算出された特徴量のうち、表示対象の特徴量とは異なる特徴量に基づいて表示対象の特徴量の表示仕様（色相スケール）を設定する（ステップＳ９：表示仕様設定ステップ）。例えば、傾きを示す特徴量ａの表示仕様である色相スケールを、ミッドバンドフィットを示す特徴量ｃに基づいて設定する。

特徴量画像データ生成部３４２は、Ｂモード画像データ生成部３４１が生成したＢモード画像データにおける各画素に対し、ステップＳ８で算出された特徴量に関連づけた視覚情報であって、ステップＳ９で設定された色相スケールを用いて、視覚情報（例えば色相）を重畳することによって特徴量画像データを生成する（ステップＳ１０：特徴量画像データ生成ステップ）。

この後、表示装置４は、制御部３６の制御のもと、特徴量画像データ生成部３４２が生成した特徴量画像データに対応する特徴量画像を表示する（ステップＳ１１）。図１０は、表示装置４における特徴量画像の表示例を模式的に示す図である。同図に示す特徴量画像２０１は、Ｂモード画像に特徴量に関する視覚情報が重畳された画像を表示する重畳画像表示部２０２と、観測対象の識別情報などを表示する情報表示部２０３とを有する。なお、情報表示部２０３に、特徴量の情報、近似式の情報、ゲインやコントラスト等の画像情報等をさらに表示するようにしてもよい。また、特徴量画像に対応するＢモード画像を特徴量画像と並べて表示してもよい。

以上説明してきた一連の処理（ステップＳ１〜Ｓ１１）において、ステップＳ２の処理とステップＳ４〜Ｓ１１の処理とを並行して行うようにしてもよい。

以上説明した本発明の実施の形態１によれば、表示仕様設定部３３４が、表示装置４に表示させる表示対象の特徴量ａとは異なる特徴量ｃに基づいて、表示対象の特徴量ａの表示仕様を設定するようにしたので、二つの特徴量の分布から最適な色相スケールを選択して、特徴量の差を明確に表現することができる。

なお、上述した実施の形態１では、表示仕様設定部３３４が、表示装置４に表示させる表示対象の特徴量とは異なる特徴量の平均値に基づいて、表示対象の特徴量の表示仕様を設定するものとして説明したが、これに限らず、最頻値や中央値に基づいて色相スケールを設定するものであってもよい。

また、上述した実施の形態１において、例えば、特徴量ａおよび特徴量ｃの分布において組織性状間で領域が重複するなど、一つの閾値のみで判断することが難しい場合は、領域ごとに設定する色相スケールの選択範囲を予め設定しておき、表示仕様設定部３３４が、表示装置４に表示させる表示対象の特徴量とは異なる特徴量の平均値がどの選択範囲に含まれているかを判断して、表示対象の特徴量の表示仕様を設定してもよい。

（実施の形態１の変形例１）
図１１は、本発明の実施の形態１の変形例１に係る超音波観測装置３の表示仕様設定部３３４が実行する処理を説明する図である。上述した実施の形態１では、各色相スケールの一端側が揃っている、例えば、特徴量ａの最小値をゼロ［ｄＢ／ＭＨｚ］に揃えている（図７参照）ものとして説明したが、各色相スケールの両端の値が異なっていてもよい。例えば、表示仕様設定部３３４が、図１１の（ａ）に示すような特徴量ａの値の範囲（配色範囲）が０〜４［ｄＢ／ＭＨｚ］の色相スケールＣＢ₃と、図１１の（ｂ）に示すような特徴量ａの値の範囲が１〜３［ｄＢ／ＭＨｚ］の色相スケールＣＢ₄とを、特徴量ｃの値に応じて設定するようにしてもよい。この際、表示仕様設定部３３４は、特徴量ｃの代表値（例えば平均値や最頻値）と、閾値とを比較し、比較結果に基づいて配色範囲を設定する。

（実施の形態１の変形例２）
図１２は、本発明の実施の形態１の変形例２に係る超音波観測装置３の表示仕様設定部３３４が実行する処理を説明する図である。上述した実施の形態１では、特徴量ａの値の配色範囲が異なる色相スケールであって、色相の変化率が同じ割合の色相スケールを用いる（図７参照）ものとして説明したが、色相スケールにおける色相の変化率が異なっていてもよい。例えば、表示仕様設定部３３４が、図１２の（ａ）および図１２の（ｂ）に示すように、特徴量ａの値の範囲が０〜４［ｄＢ／ＭＨｚ］で同じ色相スケールＣＢ₅および色相スケールＣＢ₆であって、色相の階調の変化率が異なる色相スケールＣＢ₅および色相スケールＣＢ₆を、特徴量ｃの値に応じて設定するようにしてもよい。この際、表示仕様設定部３３４は、特徴量ｃの代表値（例えば平均値や最頻値）と、閾値とを比較し、比較結果に基づいて色相の変化率を設定する。

（実施の形態１の変形例３）
図１３は、本発明の実施の形態１の変形例３に係る超音波観測装置３の表示仕様設定部３３４が実行する処理を説明する図である。上述した実施の形態１では、表示仕様設定部３３４が、関心領域内の特徴量ｃの値に応じて色相スケールを設定するものとして説明したが、関心領域のうち、対象部位の特徴量ｃの値に応じて色相スケールを設定するようにしてもよい。例えば、表示仕様設定部３３４が、図１３に示すような関心領域Ｒのうち腫瘍などの対象部位Ｓにかかる特徴量ｃの値に応じて色相スケールを設定する。本変形例３では、例えば、入力部３５を介してユーザによって選択された対象部位Ｓ内の位置における特徴量ｃの値に応じて、色相スケールが設定される。これにより、視覚情報に対応する特徴量に基づいて色相スケールを設定することになり、表示対象の特徴量について、一層高精度な表示を行うことができる。

なお、関心領域内の特徴量の値、入力部３５を介してユーザによって選択された対象部位内の位置における特徴量の値のほか、Ｂモード画像全体の特徴量の値を用いてもよい。

（実施の形態２）
図１４は、本発明の実施の形態２に係る超音波観測装置３ａを備えた超音波観測システム１ａの構成を示すブロック図である。同図に示す超音波観測システム１ａは、上述した実施の形態１に係る超音波観測システム１の超音波観測装置３に代えて、超音波観測装置３ａを備える。上述した実施の形態１では、表示仕様設定部３３４が、予め設定された特徴量の色相スケールを設定するものとして説明したが、本実施の形態２では、特徴量のばらつきに応じて色相スケールを設定する特徴量を選定する。

超音波観測装置３ａの演算部３３ａは、上述した演算部３３の構成に加えてばらつき算出部３３５および特徴量選定部３３６を有する。ばらつき算出部３３５は、特徴量算出部３３３が算出した複数の補正特徴量のばらつきを算出する。具体的には、ばらつき算出部３３５は、色相スケールの設定に用いる特徴量が特徴量ａおよび特徴量ｃである場合、特徴量ａおよび特徴量ｃについて、それぞればらつきを求める。ばらつきとしては、分散や、特徴量の最大値と最小値との差などが挙げられる。

特徴量選定部３３６は、ばらつきを求めた二つの特徴量について、ばらつきを比較して表示対象とする特徴量を選定する。具体的には、特徴量選定部３３６は、色相スケールの設定に用いる特徴量が特徴量ａおよび特徴量ｃである場合、特徴量ａおよび特徴量ｃのばらつきを比較して、ばらつきの大きい方の特徴量を表示対象の特徴量に選定する。特徴量選定部３３６は、例えば、特徴量ｃのばらつきが特徴量ａのばらつきよりも大きい場合は、特徴量ｃを表示対象の特徴量とし、特徴量ａに基づいて色相スケールを設定する。なお、ばらつきの大きい方を表示対象として選定するものとして説明したが、ばらつきの小さい方であってもよく、明確に表示できる特徴量を選定することが好ましい。

表示仕様設定部３３４は、特徴量選定部３３６が選定した特徴量の表示スケールについて、ばらつきが大きい特徴量に基づいて設定する。

以上説明した本発明の実施の形態２によれば、ばらつき算出部３３５が算出したばらつきをもとに、特徴量選定部３３６が表示対象の特徴量を選定し、表示仕様設定部３３４が、選定された特徴量（表示対象の特徴量）とは異なる特徴量に基づいて、表示対象の特徴量の表示仕様を設定するようにしたので、差をより明確に表現できる特徴量を選定して表示装置４に表示させることができる。

（実施の形態３）
図１５は、本発明の実施の形態３に係る超音波観測装置３ｂを備えた超音波観測システム１ｂの構成を示すブロック図である。同図に示す超音波観測システム１ｂは、上述した実施の形態１に係る超音波観測システム１の超音波観測装置３に代えて、超音波観測装置３ｂを備える。上述した実施の形態１では、表示仕様設定部３３４が、周波数スペクトルをもとに算出された特徴量ａ、ｃの値に応じて色相スケールを設定するものとして説明したが、本実施の形態３では、硬さの値に応じて音速を示す色相スケールを設定する。

超音波観測装置３ｂにおける演算部３３ｂは、上述した演算部３３の増幅補正部３３１、周波数解析部３３２および特徴量算出部３３３に代えて、特徴量算出部３３７を有する。特徴量算出部３３７は、送受信部３１が生成したＲＦデータをもとに、各位置における音速および硬さを算出する。特徴量算出部３３７は、ＲＦデータをもとに、超音波振動子２１から送信された超音波が観測対象で反射して戻ってくるまでの時間を算出することで、音速を求める。また、特徴量算出部３３７は、各位置における複数のＲＦデータをもとに、時間当たりの変化量を算出することで、硬さを求める。

本実施の形態３に係る表示仕様設定部３３４は、特徴量としての音速の表示仕様である色相スケールを、これとは異なる特徴量である硬さに基づいて設定する。図１６は、本発明の実施の形態３に係る超音波観測装置３ｂの表示仕様設定部３３４が実行する処理を説明する図である。図１６では、横軸が音速（ｃｍ／ｓ）である。また、図１６では、縦軸が硬さ（ｋＰａ）である。図１６は、音速および硬さの分布を示すグラフとなっている。

ここで、図１６に示すように、生体組織の性状（以下、組織性状という）の種類によって、音速および硬さの分布が異なる。例えば、ある組織性状は領域Ｑ₃に分布され、それとは異なる組織性状は領域Ｑ₄に分布される。この際、同じ色相スケールで音速を表現すると、一方の組織性状は音速の値に応じて明確に表現されるものの、他方の組織性状は音速の値の差が明確に表現されなかったり、色相スケールから外れた領域に音速の値が分布されて音速が表示されなかったりしてしまう。

本実施の形態３では、表示仕様設定部３３４が、硬さに基づいて、表示対象の音速の表示仕様を設定する。具体的に、表示仕様設定部３３４は、特徴量算出部３３７が算出した複数の硬さの平均値を求め、該平均値と閾値Ｈ_thとを比較して、比較結果に基づいて色相スケールを設定する。例えば、表示仕様設定部３３４は、硬さの平均値が閾値Ｈ_thより小さい場合、表示仕様情報記憶部３７１を参照して、領域Ｑ₃における音速の範囲に応じた色相スケールＣＢ₇の設定を行う。これに対し、表示仕様設定部３３４は、硬さの平均値が閾値Ｈ_th以上である場合、表示仕様情報記憶部３７１を参照して、領域Ｑ₄における音速の範囲に応じた色相スケールＣＢ₈の設定を行う。このようにして、音速の色相スケールを硬さに基づいて設定することによって、音速および硬さの分布に応じた組織性状に対応する色相スケールを設定することができる。

以上説明した本発明の実施の形態３によれば、表示仕様設定部３３４が、表示装置４に表示させる表示対象の特徴量（音速）とは異なる特徴量（硬さ）に基づいて、表示対象の特徴量（硬さ）の表示仕様を設定するようにしたので、二つの特徴量の分布から最適な色相スケールを選択して、特徴量の差を明確に表現することができる。

ここまで、本発明を実施するための形態を説明してきたが、本発明は上述した実施の形態によってのみ限定されるべきものではない。上述した実施の形態１〜３では、二つの特徴量のうちの一方の特徴量に基づいて他方の特徴量の色相スケールを設定するものとして説明したが、これに限らず、三つの特徴量を用いて、例えば各特徴量を直交する三軸の座標系で表現し、表示対象の特徴量の色相スケールを、他の二つの特徴量に基づいて設定するものであってもよい。この場合、閾値は、他の二つの特徴量により空間的に設定されるものとなる。また、周波数特徴量と、音速や硬さとを組としてもよいし、減衰量を特徴量として用いてもよい。本実施の形態においては、周波数特徴量である特徴量ａ、特徴量ｂおよび特徴量ｃ、ならびに音速、硬さ、減衰量のうちの少なくとも二つの特徴量を用いて色相スケールが設定される。

また、上述した実施の形態１〜３では、組織性状の分布に応じて二つの色相スケールを設定するものとして説明したが、これに限らず、組織性状に応じた分布（領域）が三つ以上あれば、設定されうる色相スケールもその数に応じて表示仕様情報記憶部３７１に記憶される。また、設定対象の複数の色相スケールにおいて色相の変化率が同じであり、複数の組織性状の各領域の中心値（平均値や中央値、最頻値）が直線上にあれば、閾値に応じて、配色範囲をスライドして設定するようにしてもよい。

また、上述した実施の形態１〜３では、ある閾値に応じて色相スケールを設定するものとして説明したが、組織性状の分布に応じた特徴量の範囲を設定し、表示対象とは異なる特徴量の値がどの範囲に含まれているかを判断して色相スケールを設定してもよい。

このように、本発明は、請求の範囲に記載した技術的思想を逸脱しない範囲内において、様々な実施の形態を含みうるものである。

以上のように、本発明にかかる超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラムは、特徴量の差を明確に表現するのに有用である。

１，１ａ，１ｂ 超音波観測システム
２ 超音波内視鏡
３，３ａ，３ｂ 超音波観測装置
４ 表示装置
２１ 超音波振動子
３１ 送受信部
３２ 信号処理部
３３，３３ａ，３３ｂ 演算部
３４ 画像処理部
３５ 入力部
３６ 制御部
３７ 記憶部
２０１ 特徴量画像
２０２ 重畳画像表示部
２０３ 情報表示部
３３１ 増幅補正部
３３２ 周波数解析部
３３３，３３７ 特徴量算出部
３３３ａ 近似部
３３３ｂ 減衰補正部
３３４ 表示仕様設定部
３３５ ばらつき算出部
３３６ 特徴量選定部
３４１ Ｂモード画像データ生成部
３４２ 特徴量画像データ生成部
３７１ 表示仕様情報記憶部
Ｃ₁ 周波数スペクトル

Claims (9)

1. 観測対象へ超音波を送信し、該観測対象で反射された超音波を受信する超音波振動子を備えた超音波プローブが取得した超音波信号に基づいて超音波画像を生成する超音波観測装置であって、
   前記超音波信号に基づいて複数の特徴量を算出する特徴量算出部と、
   前記特徴量算出部が算出した前記複数の特徴量のうち、前記超音波画像とともに表示する表示対象の特徴量を所定の表示仕様で配色した特徴量画像データを生成する特徴量画像データ生成部と、
   前記表示対象の特徴量以外の他の特徴量の代表値と閾値とを比較して、前記表示対象の特徴量の表示仕様である前記表示対象の特徴量の配色範囲、および／または色相の変化率を設定する表示仕様設定部と、
   を備えたことを特徴とする超音波観測装置。
2. 前記表示仕様設定部は、前記他の特徴量であって前記超音波信号に基づく画像内の対象領域の特徴量に基づいて、前記表示対象の特徴量の表示仕様を設定する
   ことを特徴とする請求項１に記載の超音波観測装置。
3. 前記超音波信号に基づいて生成される信号の周波数を解析することによって複数の周波数スペクトルを算出する周波数解析部をさらに備え、
   前記特徴量算出部は、前記複数の特徴量の一つとして、前記周波数解析部が算出した周波数スペクトルに基づく周波数特徴量を算出する
   ことを特徴とする請求項１に記載の超音波観測装置。
4. 前記表示対象の特徴量は、音速、硬さ、減衰量、前記周波数特徴量のうちのいずれかであり、
   前記他の特徴量は、音速、硬さ、減衰量、前記周波数特徴量のうちの前記表示対象の特徴量以外の少なくとも一つである
   ことを特徴とする請求項３に記載の超音波観測装置。
5. 前記特徴量算出部は、前記周波数解析部が算出した周波数スペクトルを直線近似することによって前記周波数特徴量を算出する
   ことを特徴とする請求項３に記載の超音波観測装置。
6. 前記周波数特徴量は、前記周波数解析部が算出した周波数スペクトルを近似することによって得られた傾き、切片およびミッドバンドフィットのうち少なくとも一つである
   ことを特徴とする請求項５に記載の超音波観測装置。
7. 前記複数の特徴量のばらつきを算出するばらつき算出部と、
   前記ばらつき算出部が算出した各特徴量のばらつきをもとに、前記表示対象の特徴量を選定する特徴量選定部と、
   をさらに備え、
   前記表示仕様設定部は、前記特徴量選定部が選定した前記表示対象の特徴量以外の他の特徴量に基づいて、前記表示対象の特徴量の表示仕様を設定する
   ことを特徴とする請求項１に記載の超音波観測装置。
8. 観測対象へ超音波を送信し、該観測対象で反射された超音波を受信する超音波振動子を備えた超音波プローブが取得した超音波信号に基づいて超音波画像を生成する超音波観測装置の作動方法であって、
   特徴量算出部が、前記超音波信号に基づいて複数の特徴量を算出する特徴量算出ステップと、
   表示仕様設定部が、前記特徴量算出部が算出した前記複数の特徴量のうち、前記超音波画像とともに表示する表示対象の特徴量以外の他の特徴量の代表値と閾値とを比較して、前記表示対象の特徴量の表示仕様である前記表示対象の特徴量の配色範囲、および／または色相の変化率を設定する表示仕様設定ステップと、
   特徴量画像データ生成部が、前記表示対象の特徴量を所定の表示仕様で配色した特徴量画像データを生成する特徴量画像データ生成ステップと、
   を含むことを特徴とする超音波観測装置の作動方法。
9. 観測対象へ超音波を送信し、該観測対象で反射された超音波を受信する超音波振動子を備えた超音波プローブが取得した超音波信号に基づいて超音波画像を生成する超音波観測装置の作動プログラムであって、
   特徴量算出部が、前記超音波信号に基づいて複数の特徴量を算出する特徴量算出手順と、
   表示仕様設定部が、前記特徴量算出部が算出した前記複数の特徴量のうち、前記超音波画像とともに表示する表示対象の特徴量以外の他の特徴量の代表値と閾値とを比較して、前記表示対象の特徴量の表示仕様である前記表示対象の特徴量の配色範囲、および／または色相の変化率を設定する表示仕様設定手順と、
   特徴量画像データ生成部が、前記表示対象の特徴量を所定の表示仕様で配色した特徴量画像データを生成する特徴量画像データ生成手順と、
   を前記超音波観測装置に実行させることを特徴とする超音波観測装置の作動プログラム。