"検体に対して超音波を送信するとともに前記検体によって反射された超音波を受信する超音波観測装置であって、受信した超音波の周波数を解析することによって周波数スペクトルを算出する周波数解析部と、前記周波数解析部が算出した周波数スペクトルを近似することによって前記周波数スペクトルの複数の特徴量を抽出する特徴量抽出部と、複数の既知検体に対してそれぞれ抽出される複数の特徴量の少なくとも一部を座標成分とする第1座標系を有する特徴量空間上において、前記複数の既知検体を各既知検体の情報に基づいて分類した複数のグループをそれぞれ代表する複数の代表点を所定の軸へ射影したときに、前記第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が大きくなる新座標軸を一つの座標軸とする第2座標系における前記検体の特徴量の座標値を算出し、算出した座標値に応じた表示パラメータを付与する表示パラメータ付与部と、を備えたことを特徴とする超音波観測装置。前記新座標軸は、前記第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が最大の値になる座標軸であることを特徴とする請求項1に記載の超音波観測装置。前記表示パラメータ付与部が前記検体の特徴量に付与した表示パラメータによって定まる画素値を有する特徴量画像データを生成する画像処理部と、前記画像処理部が生成した特徴量画像データに対応する画像を表示可能な表示部と、を備えたことを特徴とする請求項1または2に記載の超音波観測装置。前記表示パラメータは、色空間を構成する変数であることを特徴とする請求項1~3のいずれか1項に記載の超音波観測装置。前記色空間を構成する変数は、原色系の特定成分、補色系の特定成分、色相、彩度、明度のいずれかであることを特徴とする請求項4に記載の超音波観測装置。前記画像処理部は、前記新座標軸と直交する第2新座標軸上における前記検体の特徴量の座標値に対応付けられ、前記表示パラメータとは独立に画像の表示態様を決定する第2表示パラメータをさらに用いて特徴量画像データを生成することを特徴とする請求項1~5のいずれか1項に記載の超音波観測装置。前記第2表示パラメータは、色空間を構成する変数であることを特徴とする請求項6に記載の超音波観測装置。前記色空間を構成する変数は、原色系の特定成分、補色系の特定成分、色相、彩度、明度のいずれかであることを特徴とする請求項7に記載の超音波観測装置。同一のグループに属する特徴量をそれぞれ有する複数の検体は、組織性状が互いに同じであることを特徴とする請求項1~8のいずれか1項に記載の超音波観測装置。前記特徴量抽出部は、前記周波数解析部が算出した周波数スペクトルに対して前記近似処理を行うことにより、前記減衰補正処理を行う前の補正前特徴量を抽出する近似部と、前記近似部が抽出した補正前特徴量に対して前記減衰補正処理を行うことにより、前記周波数スペクトルの特徴量を抽出する減衰補正部と、を有することを特徴とする請求項1~9のいずれか1項に記載の超音波観測装置。前記特徴量抽出部は、前記周波数スペクトルに対して前記減衰補正処理を行う減衰補正部と、前記減衰補正部が補正した周波数スペクトルに対して前記近似処理を行うことにより、前記周波数スペクトルの特徴量を抽出する近似部と、を有することを特徴とする請求項1~9のいずれか1項に記載の超音波観測装置。前記特徴量抽出部は、回帰分析によって前記周波数スペクトルを多項式で近似することを特徴とする請求項10または11に記載の超音波観測装置。前記近似部は、前記周波数スペクトルを一次式で近似し、前記一次式の傾き、前記一次式の切片、および前記傾きと前記切片と前記周波数スペクトルの周波数帯域に含まれる特定の周波数とを用いて定まる強度、のうち少なくとも2つを含む複数の特徴量を抽出することを特徴とする請求項12に記載の超音波観測装置。前記減衰補正部は、超音波の受信深度が大きいほど大きな補正を行うことを特徴とする請求項10~13のいずれか1項に記載の超音波観測装置。検体に対して超音波を送信するとともに前記検体によって反射された超音波を受信する超音波観測装置が行う超音波観測装置の作動方法であって、受信した超音波の周波数を解析することによって周波数スペクトルを周波数解析部により算出する周波数解析ステップと、前記周波数解析ステップで算出した周波数スペクトルを近似することによって前記周波数スペクトルの特徴量を特徴量抽出部により複数抽出する特徴量抽出ステップと、複数の既知検体に対してそれぞれ抽出される複数の特徴量の少なくとも一部を座標成分とする第1座標系を有する特徴量空間上において、前記複数の既知検体を各既知検体の情報に基づいて分類した複数のグループをそれぞれ代表する複数の代表点を所定の軸へ射影したときに、前記第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が大きくなる新座標軸を一つの座標軸とする第2座標系における前記検体の特徴量の座標値を算出し、算出した座標値に応じた表示パラメータを表示パラメータ付与部により付与する表示パラメータ付与ステップと、を有することを特徴とする超音波観測装置の作動方法。検体に対して超音波を送信するとともに前記検体によって反射された超音波を受信する超音波観測装置に、受信した超音波の周波数を解析することによって周波数スペクトルを周波数解析部により算出する周波数解析ステップと、前記周波数解析ステップで算出した周波数スペクトルを近似することによって前記周波数スペクトルの特徴量を特徴量抽出部により複数抽出する特徴量抽出ステップと、複数の既知検体に対してそれぞれ抽出される複数の特徴量の少なくとも一部を座標成分とする第1座標系を有する特徴量空間上において、前記複数の既知検体を各既知検体の情報に基づいて分類した複数のグループをそれぞれ代表する複数の代表点を所定の軸へ射影したときに、前記第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が大きくなる新座標軸を一つの座標軸とする第2座標系における前記検体の特徴量の座標値を算出し、算出した座標値に応じた表示パラメータを表示パラメータ付与部により付与する表示パラメータ付与ステップと、を実行させることを特徴とする超音波観測装置の作動プログラム。しかしながら、上述した超音波エラストグラフィでは、血管やリンパ管などの脈管の下部には押し付ける圧力が伝わりにくいという問題があった。そのため、脈管の近傍に腫瘍が形成されている場合、腫瘍の境界が不明りょうであり、脈管内への腫瘍の浸潤の鑑別も難しかった。このように、超音波エラストグラフィでは、組織性状の鑑別等の検体の観測を精度よく行うことができない場合があった。また、超音波エラストグラフィでは、検査者が検査箇所を圧迫する際の圧力や圧迫速度に個人差が生じやすいため、測定結果の信頼性が低いという問題もあった。本発明は、上記に鑑みてなされたものであって、検体を精度よく観測することができるとともに、観測結果の信頼性を向上させることができる超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラムを提供することを目的とする。上述した課題を解決し、目的を達成するために、本発明に係る超音波観測装置は、検体に対して超音波を送信するとともに前記検体によって反射された超音波を受信する超音波観測装置であって、受信した超音波の周波数を解析することによって周波数スペクトルを算出する周波数解析部と、前記周波数解析部が算出した周波数スペクトルを近似することによって前記周波数スペクトルの複数の特徴量を抽出する特徴量抽出部と、複数の既知検体に対してそれぞれ抽出される複数の特徴量の少なくとも一部を座標成分とする第1座標系を有する特徴量空間上において、前記複数の既知検体を各既知検体の情報に基づいて分類した複数のグループをそれぞれ代表する複数の代表点を所定の軸へ射影したときに、前記第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が大きくなる新座標軸を一つの座標軸とする第2座標系における前記検体の特徴量の座標値を算出し、算出した座標値に応じた表示パラメータを付与する表示パラメータ付与部と、を備えたことを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記新座標軸は、前記第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が最大の値になる座標軸であることを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記表示パラメータ付与部が前記検体の特徴量に付与した表示パラメータによって定まる画素値を有する特徴量画像データを生成する画像処理部と、前記画像処理部が生成した特徴量画像データに対応する画像を表示可能な表示部と、を備えたことを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記表示パラメータは、色空間を構成する変数であることを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記色空間を構成する変数は、原色系の特定成分、補色系の特定成分、色相、彩度、明度のいずれかであることを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記画像処理部は、前記新座標軸と直交する第2新座標軸上における前記検体の特徴量の座標値に対応付けられ、前記表示パラメータとは独立に画像の表示態様を決定する第2表示パラメータをさらに用いて特徴量画像データを生成することを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記第2表示パラメータは、色空間を構成する変数であることを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記色空間を構成する変数は、原色系の特定成分、補色系の特定成分、色相、彩度、明度のいずれかであることを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、同一のグループに属する特徴量をそれぞれ有する複数の検体は、組織性状が互いに同じであることを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記特徴量抽出部は、前記周波数解析部が算出した周波数スペクトルに対して前記近似処理を行うことにより、前記減衰補正処理を行う前の補正前特徴量を抽出する近似部と、前記近似部が抽出した補正前特徴量に対して前記減衰補正処理を行うことにより、前記周波数スペクトルの特徴量を抽出する減衰補正部と、を有することを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記特徴量抽出部は、前記周波数スペクトルに対して前記減衰補正処理を行う減衰補正部と、前記減衰補正部が補正した周波数スペクトルに対して前記近似処理を行うことにより、前記周波数スペクトルの特徴量を抽出する近似部と、を有することを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記特徴量抽出部は、回帰分析によって前記周波数スペクトルを多項式で近似することを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記近似部は、前記周波数スペクトルを一次式で近似し、前記一次式の傾き、前記一次式の切片、および前記傾きと前記切片と前記周波数スペクトルの周波数帯域に含まれる特定の周波数とを用いて定まる強度、のうち少なくとも2つを含む複数の特徴量を抽出することを特徴とする。また、本発明に係る超音波観測装置は、上記発明において、前記減衰補正部は、超音波の受信深度が大きいほど大きな補正を行うことを特徴とする。また、本発明に係る超音波観測装置の作動方法は、検体に対して超音波を送信するとともに前記検体によって反射された超音波を受信する超音波観測装置が行う超音波観測装置の作動方法であって、受信した超音波の周波数を解析することによって周波数スペクトルを周波数解析部により算出する周波数解析ステップと、前記周波数解析ステップで算出した周波数スペクトルを近似することによって前記周波数スペクトルの特徴量を特徴量抽出部により複数抽出する特徴量抽出ステップと、複数の既知検体に対してそれぞれ抽出される複数の特徴量の少なくとも一部を座標成分とする第1座標系を有する特徴量空間上において、前記複数の既知検体を各既知検体の情報に基づいて分類した複数のグループをそれぞれ代表する複数の代表点を所定の軸へ射影したときに、前記第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が大きくなる新座標軸を一つの座標軸とする第2座標系における前記検体の特徴量の座標値を算出し、算出した座標値に応じた表示パラメータを表示パラメータ付与部により付与する表示パラメータ付与ステップと、を有することを特徴とする。また、本発明に係る超音波観測装置の作動プログラムは、検体に対して超音波を送信するとともに前記検体によって反射された超音波を受信する超音波観測装置に、受信した超音波の周波数を解析することによって周波数スペクトルを周波数解析部により算出する周波数解析ステップと、前記周波数解析ステップで算出した周波数スペクトルを近似することによって前記周波数スペクトルの特徴量を特徴量抽出部により複数抽出する特徴量抽出ステップと、複数の既知検体に対してそれぞれ抽出される複数の特徴量の少なくとも一部を座標成分とする第1座標系を有する特徴量空間上において、前記複数の既知検体を各既知検体の情報に基づいて分類した複数のグループをそれぞれ代表する複数の代表点を所定の軸へ射影したときに、前記第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が大きくなる新座標軸を一つの座標軸とする第2座標系における前記検体の特徴量の座標値を算出し、算出した座標値に応じた表示パラメータを表示パラメータ付与部により付与する表示パラメータ付与ステップと、を実行させることを特徴とする。本発明によれば、複数の既知検体に対してそれぞれ抽出される複数の特徴量の少なくとも一部を座標成分とする第1座標系を有する特徴量空間上において、複数の既知検体を各既知検体の情報に基づいて分類した複数のグループをそれぞれ代表する複数の代表点を所定の軸へ射影したときに、第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が大きくなる新座標軸を一つの座標軸とする第2座標系における検体の特徴量の座標値を算出し、算出した座標値に応じた表示パラメータを付与するため、グループの違いを明確に峻別することができる。したがって、検体を精度よく観測することができるとともに、観測結果の信頼性を向上させることができる。以下、添付図面を参照して、本発明を実施するための形態以下、実施の形態というを説明する。実施の形態1図1は、本発明の実施の形態1に係る超音波観測装置の構成を示すブロック図である。同図に示す超音波観測装置1は、超音波を用いて検体を観測する装置である。超音波観測装置1は、外部へ超音波パルスを出力するとともに、外部で反射された超音波エコーを受信する超音波探触子2と、超音波探触子2との間で電気信号の送受信を行う送受信部3と、超音波エコーを変換した電気的なエコー信号に対して所定の演算を施す演算部4と、超音波エコーを変換した電気的なエコー信号に対応する画像データの生成を行う画像処理部5と、キーボード、マウス、タッチパネル等のインタフェースを用いて実現され、各種情報の入力を受け付ける入力部6と、液晶または有機EL等からなる表示パネルを用いて実現され、画像処理部5が生成した画像を含む各種情報を表示可能な表示部7と、複数の既知検体に関する情報を記憶するとともに複数の既知検体を複数のグループに分けて記憶する記憶部8と、超音波観測装置1の動作制御を行う制御部9と、を備える。超音波探触子2は、送受信部3から受信した電気的なパルス信号を超音波パルス音響パルス信号に変換するとともに、外部の検体で反射された超音波エコーを電気的なエコー信号に変換する信号変換部21を有する。超音波探触子2は、超音波振動子をメカ的に走査させるものであってもよいし、複数の超音波振動子を電子的に走査させるものであってもよい。送受信部3は、超音波探触子2と電気的に接続され、パルス信号を超音波探触子2へ送信するとともに、超音波探触子2から受信信号であるエコー信号を受信する。具体的には、送受信部3は、予め設定された波形および送信タイミングに基づいてパルス信号を生成し、この生成したパルス信号を超音波探触子2へ送信する。送受信部3は、超音波探触子2と電気的に接続され、パルス信号を超音波探触子2へ送信するとともに、超音波探触子2からエコー信号を受信する。具体的には、送受信部3は、予め設定された波形および送信タイミングに基づいてパルス信号を生成し、この生成したパルス信号を超音波探触子2へ送信する。また、送受信部3は、受信したエコー信号に増幅、フイルタリング等の処理を施した後、A/D変換することによってデジタルRF信号を生成して出力する。なお、超音波探触子2が複数の超音波振動子を電子的に走査させるものである場合、送受信部3は、複数の超音波振動子に対応したビーム合成用の多チャンネル回路を有する。演算部4は、送受信部3が出力したデジタルRF信号に高速フーリエ変換FFTを施すことによってエコー信号の周波数解析を行う周波数解析部41と、周波数解析部41が算出した周波数スペクトルパワースペクトルに対し、超音波が伝播する際に該超音波の受信深度および周波数に応じて発生する減衰の寄与を削減する減衰補正処理と近似処理とを行うことにより、検体の特徴量を抽出する特徴量抽出部42と、検体の特徴量に対して画像表示する際の表示パラメータを付与する表示パラメータ付与部43と、を有する。周波数解析部41は、各音線ラインデータに対し、所定のデータ量からなるFFTデータ群を高速フーリエ変換することによって周波数スペクトルを算出する。周波数スペクトルは、検体の組織性状によって異なる傾向を示す。これは、周波数スペクトルが、超音波を散乱する散乱体としての検体の大きさ、密度、音響インピーダンス等と相関を有しているためである。特徴量抽出部42は、周波数解析部41が算出した周波数スペクトルに対し、近似処理を行うことによって減衰補正処理を行う前の補正前特徴量を算出する近似部421と、近似部421が近似した補正前特徴量に対して減衰補正処理を行うことによって特徴量を抽出する減衰補正部422と、を有する。近似部421は、回帰分析によって周波数スペクトルを一次式で近似し、この近似した一次式を特徴付ける特徴量を抽出する。具体的には、近似部421は、回帰分析によって一次式の傾きa0および切片b0を算出するとともに、周波数スペクトルにおける周波数帯域内の特定周波数における強度を補正前特徴量として算出する。本実施の形態1において、近似部421は、中心周波数fMID=fLOWfHIGH/2における強度Mid-bandfitc0=a0fMIDb0を算出するものとするが、これはあくまでも一例に過ぎない。ここでいう強度とは、電圧、電力、音圧、音響エネルギー等のパラメータのいずれかを指す。三つの特徴量のうち、傾きa0は、超音波の散乱体の大きさと相関を有し、一般に散乱体が大きいほど傾きが小さな値を有すると考えられる。また、切片b0は、散乱体の大きさ、音響インピーダンスの差、散乱体の密度濃度等と相関を有している。具体的には、切片b0は、散乱体が大きいほど大きな値を有し、音響インピーダンスが大きいほど大きな値を有し、散乱体の密度濃度が大きいほど大きな値を有すると考えられる。中心周波数fMIDにおける強度以下、単に強度というc0は、傾きa0と切片b0から導出される間接的なパラメータであり、有効な周波数帯域内の中心におけるスペクトル強度を与える。このため、強度c0は、散乱体の大きさ、音響インピーダンスの差、散乱体の密度に加えて、Bモード画像の輝度とある程度の相関を有していると考えられる。なお、特徴量抽出部42が算出する近似多項式は一次式に限定されるわけではなく、二次以上の近似多項式を用いることも可能である。減衰補正部422が行う補正について説明する。超音波の減衰量Aは、A=2αzf1と表すことができる。ここで、αは減衰率であり、zは超音波の受信深度であり、fは周波数である。式1からも明らかなように、減衰量Aは、周波数fに比例している。減衰率αの具体的な値は、生体の場合、0~1.0dB/cm/MHz、より好ましくは0.3~0.7dB/cm/MHzであり、観察対象の臓器の種類に応じて定まる。例えば、観察対象の臓器が膵臓である場合、α=0.6dB/cm/MHzと定められる。なお、本実施の形態1において、減衰率αの値を入力部6からの入力によって変更することが可能な構成とすることも可能である。減衰補正部422は、近似部421が抽出した補正前特徴量傾きa0,切片b0,強度c0を、次のように補正する。a=a02αz2b=b03c=c02αzfMID=afMIDb4式2、4からも明らかなように、減衰補正部422は、超音波の受信深度zが大きいほど、補正量が大きい補正を行う。また、式3によれば、切片に関する補正は恒等変換である。これは、切片が周波数0Hzに対応する周波数成分であって減衰を受けないためである。表示パラメータ付与部43は、特徴量抽出部42が抽出し、減衰補正部422が補正を加えた特徴量を座標成分とする特徴量空間において、検体の組織性状を反映させるために所定の条件を満たすように設定される第2座標系における特徴量の座標値を算出し、算出した座標値に応じて画像の表示態様を決定する表示パラメータを付与する。ここでいう組織性状とは、例えば癌、内分泌腫瘍、粘液性腫瘍、正常組織、脈管などのいずれかであり、検体が膵臓である場合には、組織性状として慢性膵炎、自己免疫性膵炎なども含まれる。第2座標系を含む特徴量空間の情報は、後述する記憶部8が記憶している。画像処理部5は、エコー信号からBモード画像データを生成するBモード画像データ生成部51と、Bモード画像データ生成部51が生成したBモード画像データおよび表示パラメータ付与部43が検体の特徴量に付与した表示パラメータによって定まる画素値を有する特徴量画像データを生成する特徴量画像データ生成部52と、を有する。Bモード画像データ生成部51は、デジタル信号に対してバンドパスフィルタ、対数変換、ゲイン処理、コントラスト処理等の公知の技術を用いた信号処理を行うとともに、表示部7における画像の表示レンジに応じて定まるデータステップ幅に応じたデータの間引き等を行うことによってBモード画像データを生成する。特徴量画像データ生成部52は、Bモード画像データ生成部51が生成したBモード画像データ、特徴量抽出部42が算出した後で減衰補正部422が補正した特徴量、および表示パラメータ付与部43が特徴量に対して付与した表示パラメータを用いることにより、特徴量画像データを生成する。記憶部8は、既知検体の特徴量を含む既知検体情報を記憶する既知検体情報記憶部81と、周波数解析部41が行う周波数解析処理の際に使用する窓関数を記憶する窓関数記憶部82と、減衰補正部422が処理を行う際に参照する補正情報を記憶する補正情報記憶部83と、既知検体情報記憶部81が記憶する既知検体の特徴量をもとに設定される特徴量空間に関する情報を記憶する特徴量空間情報記憶部84と、表示パラメータ付与部43が算出する新座標軸の座標値と表示パラメータとの関係を含む表示パラメータ情報を記憶する表示パラメータ情報記憶部85と、を有する。既知検体情報記憶部81は、既知検体に対して抽出された周波数スペクトルの特徴量を既知検体の組織性状と関連付けて記憶している。ここで、既知検体の特徴量は、本実施の形態1と同様の処理によって抽出されたものである。ただし、既知検体の特徴量抽出処理を超音波観測装置1で行う必要はない。また、既知検体情報記憶部81は、既知検体に関連した周波数スペクトルの特徴量に対し、既知検体の組織性状を含む情報に基づいて分類されたグループごとに算出された平均および標準偏差を、既知検体の特徴量の全データとともに記憶している。本実施の形態1では、超音波受信信号の周波数スペクトルの特徴量の平均および標準偏差が、検体における核の腫大や異形などの細胞レベルの変化、間質における線維の増生や実質組織の線維への置換などの組織的な変化を反映しており、組織性状に応じて特有の値を示すことに鑑みて、既知検体の周波数スペクトルの特徴量の平均および標準偏差を用いて組織性状の分類を行っている。窓関数記憶部82は、Hamming,Hanning,Blackmanなどの窓関数のうち少なくともいずれか一つの窓関数を記憶している。補正情報記憶部83は、式2~4の変換に関する情報を記憶している。特徴量空間情報記憶部84は、既知検体情報記憶部81が記憶する既知検体情報に基づいて設定される特徴量空間に関する情報として、複数の既知検体を特徴量に基づいて分類することによって得られる複数のグループをそれぞれ代表する複数の代表点を所定の軸へ射影したときに、該座標軸に沿って隣接する代表点間の距離の和と比較して、自身が延びる方向に沿って隣接する代表点間の距離の和が大きい新座標軸を一つの座標軸とする第2座標系に関する情報を記憶する。図2は、特徴量空間情報記憶部84が記憶する特徴量空間情報の概要を示す図である。図2に示す特徴量空間は、横軸が切片b、縦軸が強度cである式3、4を参照。また、図2に示す領域Gμ,Gν,Gρは、既知検体情報記憶部81が記憶する既知検体の組織性状が、それぞれμ,ν,ρであるグループを示している。図2に示す場合、三つのグループGμ,Gν,Gρは、特徴量空間上において、互いに他のグループと交わりを有しない領域に存在している。このように、本実施の形態1では、周波数解析によって得られた周波数スペクトルの特徴量を指標としてグループ分けを行うため、互いに異なるグループを峻別することができる。特に、本実施の形態1では、超音波のエコー信号に対して減衰補正を行っているため、減衰補正を行わない場合と比較して、特徴量空間における各グループの領域もより明確に分離した状態で得ることができる。なお、特徴量空間におけるb軸成分とc軸成分のスケールが大きく異なる場合には、各距離の寄与を略均等にするための重み付けを適宜行うことが望ましい。図2では、第1座標系b,cとは異なる第2座標系h,vも記載されている。この第2座標系の座標軸hは、グループGμ,Gν,Gρにそれぞれ含まれるFFTデータ群の周波数スペクトルの切片bおよび強度cの各平均を特徴量空間の座標として有する点μ0,ν0,ρ0以下、これらの点を代表点というを射影したときにその射影方向に沿って隣接する代表点同士の距離の和が最大となる軸新座標軸である。すなわち、h軸は、代表点μ0,ν0,ρ0のh軸方向の成分をμ0h,ν0h,ρ0hとしたとき、|μ0-ν0||ν0-ρ0|が最大となるような方向に沿った軸として定められる。なお、h軸と直交するv軸第2新座標軸は、設定しなくてもよい。表示パラメータ情報記憶部85は、特徴量空間情報記憶部84が記憶する特徴量空間において、上述した新座標軸の座標値と画像の表示態様を決定する表示パラメータとを対応付けて記憶する。本実施の形態1においては、例えばh軸における座標値を、光の3属性の一つである色相と対応付けて記憶している。また、本実施の形態1においては、v軸の座標値に対して光の3属性で色相と独立に定められる明度を対応付けて記憶している。なお、表示パラメータは、上述したものに限られるわけではない。例えば、光の3属性の残りの一つである彩度を表示パラメータとしてもよいし、より一般に色空間を構成する変数例えば、RGB表色系や補色系の変数を表示パラメータとしてもよい。また、模様を表示パラメータとしてもよい。模様を表示パラメータとする場合には、座標値の帯域ごとに模様が変化するように設定しておけばよい。記憶部8は、本実施の形態1に係る超音観測装置の作動プログラムや所定のOSを起動するプログラム等が予め記憶されたROM、および各処理の演算パラメータやデータ等を記憶するRAM等を用いて実現される。以上の機能構成を有する超音波観測装置1の超音波探触子2以外の構成要素は、演算および制御機能を有するCPUを備えたコンピュータを用いて実現される。超音波観測装置1が備えるCPUは、記憶部8が記憶、格納する情報および上述した超音波観測装置の作動プログラムを含む各種プログラムを記憶部8から読み出すことにより、本実施の形態1に係る超音波観測装置の作動方法に関連した演算処理を実行する。なお、本実施の形態1に係る超音波観測装置の作動プログラムは、ハードディスク、フラッシュメモリ、CD-ROM、DVD-ROM、フレキシブルディスク等のコンピュータ読み取り可能な記録媒体に記録して広く流通させることも可能である。図3は、以上の構成を有する超音波観測装置1の処理の概要を示すフローチャートである。図3において、超音波観測装置1は、まず超音波探触子2によって新規の検体の測定を行うステップS1。この後、Bモード画像データ生成部51は、送受信部3から出力されたBモード画像用エコー信号を用いてBモード画像データを生成するステップS2。続いて、制御部9は、Bモード画像データ生成部51が生成したBモード画像データに対応するBモード画像を表示部7に表示させる制御を行うステップS3。図4は、表示部7におけるBモード画像の表示例を示す図である。同図に示すBモード画像100は、色空間としてRGB表色系を採用した場合の変数であるR赤、G緑、B青の値を一致させたグレースケール画像である。その後、周波数解析部41は、FFT演算による周波数解析を行うことによって周波数スペクトルを算出するステップS4。ここで、周波数解析部41が行う処理ステップS4について、図5に示すフローチャートを参照して詳細に説明する。まず、周波数解析部41は、最初に解析対象とする音線の音線番号Lを初期値L0とするステップS21。初期値L0は、例えば送受信部3が最初に受信する音線に対して付与してもよいし、入力部6によって設定される関心領域の左右の一方の境界位置に対応する音線に対して付与してもよい。続いて、周波数解析部41は、一つの音線上に設定した複数のデータ位置全ての周波数スペクトルを算出する。まず、周波数解析部41は、FFT演算用に取得する一連のデータ群FFTデータ群を代表するデータ位置Z受信深度に相当の初期値Z0を設定するステップS22。図6は、一つの音線のデータ配列を模式的に示す図である。同図に示す音線LDにおいて、白または黒の長方形は、一つのデータを意味している。音線LDは、送受信部3が行うA/D変換におけるサンプリング周波数例えば50MHzに対応した時間間隔で離散化されている。図6では、音線LDの1番目のデータをデータ位置Zの初期値Z0として設定した場合を示している。なお、図6はあくまでも一例に過ぎず、初期値Z0の位置は任意に設定することができる。例えば、関心領域の上端位置に対応するデータ位置Zを初期値Z0として設定してもよい。その後、周波数解析部41は、データ位置ZのFFTデータ群を取得しステップS23、取得したFFTデータ群に対し、窓関数記憶部82が記憶する窓関数を作用させるステップS24。このようにFFTデータ群に対して窓関数を作用させることにより、FFTデータ群が境界で不連続になることを回避し、アーチファクトが発生するのを防止することができる。続いて、周波数解析部41は、データ位置ZのFFTデータ群が正常なデータ群であるか否かを判定するステップS25。ここで、FFTデータ群は、2のべき乗のデータ数を有している必要がある。以下、FFTデータ群のデータ数を2nnは正の整数とする。FFTデータ群が正常であるとは、データ位置ZがFFTデータ群で前から2n-1番目の位置であること意味する。換言すると、FFTデータ群が正常であるとは、データ位置Zの前方に2n-1-1=Nとする個のデータがあり、データ位置Zの後方に2n-1=Mとする個のデータがあることを意味する。図6に示す場合、FFTデータ群F2、F3、FK-1は正常である一方、FFTデータ群F1、FKは異常である。ただし、図6ではn=4N=7,M=8としている。ステップS25における判定の結果、データ位置ZのFFTデータ群が正常である場合ステップS25:Yes、周波数解析部41は、後述するステップS27へ移行する。ステップS25における判定の結果、データ位置ZのFFTデータ群が正常でない場合ステップS25:No、周波数解析部41は、不足分だけゼロデータを挿入することによって正常なFFTデータ群を生成するステップS26。ステップS25において正常でないと判定されたFFTデータ群は、ゼロデータを追加する前に窓関数が作用されている。このため、FFTデータ群にゼロデータを挿入してもデータの不連続は生じない。ステップS26の後、周波数解析部41は、後述するステップS27へ移行する。ステップS27において、周波数解析部41は、FFTデータ群を用いてFFT演算を行うことにより、周波数スペクトルを得るステップS27。図7および図8は、周波数解析部41が算出した周波数スペクトルの例を示す図である。図7および図8では、横軸fが周波数であり、縦軸Iが強度である。図7および図8にそれぞれ示す周波数スペクトル曲線C1およびC2において、周波数スペクトルの下限周波数fLOWおよび上限周波数fHIGHは、超音波探触子2の周波数帯域、送受信部3が送信するパルス信号の周波数帯域などをもとに決定されるパラメータであり、例えばfLOW=3MHz、fHIGH=10MHzである。なお、図7に示す直線L1および図8に示す直線L2については、後述する特徴量抽出処理で説明する。ところで、本実施の形態1において、曲線および直線は、離散的な点の集合からなる。この点については、後述する実施の形態においても同様である。続いて、周波数解析部41は、データ位置Zに所定のデータステップ幅Dを加算して次の解析対象のFFTデータ群のデータ位置Zを算出するステップS28。ここでのデータステップ幅Dは、Bモード画像データ生成部51がBモード画像データを生成する際に利用するデータステップ幅と一致させることが望ましいが、周波数解析部41における演算量を削減したい場合には、Bモード画像データ生成部51が利用するデータステップ幅より大きい値を設定してもよい。図6では、D=15の場合を示している。その後、周波数解析部41は、データ位置Zが最終データ位置Zmaxより大きいか否かを判定するステップS29。ここで、最終データ位置Zmaxは、音線LDのデータ長としてもよいし、関心領域の下端に対応するデータ位置としてもよい。判定の結果、データ位置Zが最終データ位置Zmaxより大きい場合ステップS29:Yes、周波数解析部41は、音線番号Lを1だけ増加するステップS30。一方、データ位置Zが最終データ位置Zmax以下である場合ステップS29:No、周波数解析部41はステップS23へ戻る。このようにして、周波数解析部41は、一つの音線LDに対して、[{Zmax-Z0/D}1]=K個のFFTデータ群に対するFFT演算を行う。ここで、[X]は、Xを超えない最大の整数を表す。ステップS30で増加した後の音線番号Lが最終音線番号Lmaxより大きい場合ステップS31:Yes、周波数解析部41は図3に示すメインルーチンへ戻る。一方、ステップS30で増加した後の音線番号Lが最終音線番号Lmax以下である場合ステップS31:No、周波数解析部41はステップS22へ戻る。このようにして、周波数解析部41は、Lmax-L01本の音線の各々についてK回のFFT演算を行う。なお、最終音線番号Lmaxは、例えば送受信部3が受信する最終の音線に付与してもよいし、関心領域の左右のいずれか一方の境界に対応する音線に付与してもよい。以下、周波数解析部41が全ての音線に対して行うFFT演算の総数Lmax-L01×KをPとおく。以上説明したステップS4の周波数解析処理に続いて、近似部421は、近似処理として周波数解析部41が算出したP個の周波数スペクトルを回帰分析することにより、補正前特徴量を抽出するステップS5。具体的には、近似部421は、周波数帯域fLOW<f<fHIGHの周波数スペクトルを近似する一次式を回帰分析によって算出することにより、この一次式を特徴づける傾きa0,切片b0,強度c0を補正前特徴量として抽出する。図7に示す直線L1および図8に示す直線L2は、このステップS5において、周波数スペクトル曲線C1およびC2に対して回帰分析をそれぞれ行うことによって得られる回帰直線である。この後、減衰補正部422は、近似部421が抽出した補正前特徴量に対して減衰補正処理を行うステップS6。例えば、データのサンプリング周波数が50MHzである場合、データのサンプリングの時間間隔は20nsecである。ここで、音速を1530m/secとすると、データのサンプリング距離間隔は、1530m/sec×20nsec/2=0.0153mmとなる。処理対象のFFTデータ群のデータ位置までの音線LDの1番目のデータからのデータステップ数がkであるとすると、そのデータ位置Zは0.0153kmmとなる。減衰補正部422は、このようにして求まるデータ位置Zの値を上述した式2~4の受信深度zへ代入することにより、周波数スペクトルの特徴量である傾きa,切片b,強度cを算出する。図9は、図7に示す直線L1に関連する特徴量に対して減衰補正を行った後の特徴量から定まる直線を示す図である。図9に示す直線L1’を表す式は、I=afb=a02αZfb05である。この式5からも明らかなように、直線L1’は、直線L1と比較して、傾きが大きく、かつ切片の値が同じである。この後、表示パラメータ付与部43は、特徴量抽出部42によって抽出され、減衰補正部422によって補正された特徴量、既知検体情報記憶部81が記憶する既知検体情報および特徴量空間情報記憶部84が記憶する特徴量空間情報に基づいて、検体の特徴量の第2座標系における座標値を算出するステップS7。図10は、この場合の座標値算出処理の概要を示す図である。具体的には、図10は、図2に示す特徴量空間上で、観測対象検体に対して抽出された特徴量を示す点以下、検体点というSpの第2座標系における座標値hS,vSを算出する状況を示している。続いて、表示パラメータ付与部43は、ステップS7で算出した第2座標系の座標値に応じた表示パラメータを付与するステップS8。この際、表示パラメータ付与部43は、表示パラメータ情報記憶部85がそれぞれ記憶する情報を参照して表示パラメータを付与する。以上説明したステップS8の後、特徴量画像データ生成部52は、Bモード画像データ生成部51が生成したBモード画像データ、および表示パラメータ付与部43が画素ごとに付与した表示パラメータを用いることにより、特徴量画像データを生成するステップS9。続いて、表示部7は、特徴量画像データ生成部52が生成した特徴量画像を表示するステップS10。図11は、表示部7が表示する特徴量画像の表示例を示す図である。図12は、図11に示す画像を白黒で模式的に示す図である。これらの図に示す特徴量画像200は、Bモード画像100と比較して、画像がカラー化しており、グループに応じた色の違いが鮮明になっている。特徴量画像200は、大別して緑色系領域200gおよび赤色系領域200rからなり、2つの領域の境界部は黄色系の色で表示されている図12では図示せず。図11に示すように、各領域は単一の色によって構成されているわけではない。例えば、緑色系領域200gは、緑色に近い色からなる画素が集まっている領域である。同様に、赤色系領域200rは、赤色に近い色からなる画素が集まっている領域である。このような特徴量画像200を見た観測者は、グループの違いすなわち組織性状の違いを明確に認識することができる。図13は、本実施の形態1における減衰補正の効果を説明する図である。図13に示す画像300は、減衰補正を行わなかった場合の特徴量画像である。この場合の特徴量画像は、Bモード画像データ生成部51が生成したBモード画像に対して、切片bをR赤,G緑,B青に対して均等に割り当てたグレースケール画像である。特徴量画像300は、受信深度が大きい領域図の下方領域で減衰の影響によって信号強度が下がり、画像が暗くなっている。これに対し、同じBモード画像を用いて減衰補正を行った特徴量画像400では、画面全体にわたって均一な明るさの画像が得られていることがわかる。以上説明した本発明の実施の形態1によれば、複数の既知検体に対してそれぞれ抽出される複数の特徴量の少なくとも一部を座標成分とする第1座標系を有する特徴量空間上において、複数の既知検体を各既知検体の情報に基づいて分類した複数のグループをそれぞれ代表する複数の代表点を所定の軸へ射影したときに、第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が大きくなる新座標軸を一つの座標軸とする第2座標系における検体の特徴量の座標値を算出し、算出した座標値に応じた表示パラメータを付与するため、グループの違いを明確に峻別することができる。したがって、検体を精度よく観測することができるとともに、観測結果の信頼性を向上させることができる。また、本実施の形態1によれば、第2座標系における検体の特徴量の座標値を算出し、座標値に応じた表示パラメータを付与し、付与した表示パラメータによって定まる画素値を有する特徴量画像データを生成して表示するため、利用者は画像中におけるグループの違いを明確に認識することができる。また、本実施の形態1によれば、周波数スペクトルから抽出した補正前特徴量に対して減衰補正を施しているため、超音波の伝播に伴う減衰の影響を取り除くことができ、一段と高精度な観測を行うことができる。なお、本実施の形態1において、図3のステップS6で減衰補正部422が特徴量の減衰補正を行った後、ステップS7へ移行する前に、特徴量に基づいて検体の組織性状を判定するようにしてもよい。ただし、この場合には、既知検体情報記憶部81が、既知検体の特徴量を組織性状と関連付けて記憶しておく必要がある。ここで、超音波観測装置1が組織性状判定を行う場合の具体的な処理を説明する。超音波観測装置1の演算部4は、特徴量空間における検体点Spと各グループの代表点μ0,ν0,ρ0との距離を算出し、その距離が最小となるグループに検体点Spが属していると判定する。なお、座標値と代表点との距離が極端に離れている場合には、たとえ最小値が求まったとしても組織性状の判定結果の信頼度は低い。そこで、検体点Spと代表点との距離が所定の閾値より大きい場合、超音波観測装置1がエラー信号を出力するようにしてもよい。また、検体点Spと代表点との距離の最小値が二つ以上ある場合、超音波観測装置1は、最小値に対応するすべての組織性状を候補として選択してもよいし、所定の規則にしたがっていずれか一つの組織性状を選択してもよい。後者の場合、例えば癌などの悪性の高い組織の優先順位を高く設定する方法を挙げることができる。また、検体点Spと代表点との距離の最小値が二つ以上ある場合、超音波観測装置1はエラー信号を出力してもよい。図14は、超音波観測装置1が特徴量に基づく組織性状の判定を行う場合に、表示部7が表示する判定結果表示画像の表示例を示す図である。図15は、図14に示す画像を白黒で模式的に示す図である。これらの図に示す判定結果表示画像500は、組織性状の判定結果を含む各種関連情報を表示する情報表示部501と、特徴量画像を表示する画像表示部502とを有する。図14および図15では、画像表示部502が、図11および図12に示す特徴量画像200を表示している場合を示している。情報表示部501には、例えば検体の識別情報ID番号、名前、性別等、組織性状判定結果、組織性状判定を行う際の特徴量に関する情報、ゲインやコントラスト等の超音波画質情報が表示される。ここで、特徴量に関する情報として、関心領域の内部に位置するQ組のFFTデータ群の周波数スペクトルの特徴量の平均、標準偏差を利用した表示を行うことが可能である。具体的には、情報表示部501では、例えば傾き=1.5±0.3dB/MHz、切片=-60±2dB、強度=-50±1.5dB、と表示することができる。以上の構成を有する判定結果表示画像500を表示部7が表示することにより、操作者はより正確に関心領域の組織性状を把握することが可能となる。なお、判定結果表示画像は上述した構成に限られるわけではない。たとえば、判定結果表示画像として、組織性状強調画像とBモード画像とを並べて表示するようにしてもよい。これにより、両画像の違いを一つの画面上で認識することができる。実施の形態2本発明の実施の形態2は、特徴量抽出部が行う特徴量抽出処理が、実施の形態1と異なる。本実施の形態2に係る超音波観測装置の構成は、実施の形態1で説明した超音波観測装置1の構成と同様である。そこで、以下の説明において、超音波観測装置1の構成要素と対応する構成要素には同一の符号を付すものとする。本実施の形態2における特徴量抽出処理においては、まず減衰補正部422が、周波数解析部41によって算出された周波数スペクトルに対して減衰補正処理を行う。その後、近似部421は、減衰補正部422によって減衰補正された周波数スペクトルに対して近似処理を行うことにより、周波数スペクトルの特徴量を抽出する。図16は、本実施の形態2に係る超音波観測装置の処理の概要を示すフローチャートである。図16において、ステップS41~S44の処理は、図3のステップS1~S4の処理に順次対応している。ステップS45において、減衰補正部422は、周波数解析部41がFFT演算によって算出した周波数スペクトルに対して減衰補正を行うステップS45。図17は、このステップS45の処理の概要を模式的に示す図である。図17に示すように、減衰補正部422は、周波数スペクトル曲線C3に対し、上述した式1の減衰量Aを強度Iに加える補正を全ての周波数fに対して行うことにより、新たな周波数スペクトル曲線C3’を得る。これにより、超音波の伝播に伴う減衰の寄与を削減した周波数スペクトルを得ることができる。この後、近似部421は、減衰補正部422によって減衰補正された全ての周波数スペクトルを回帰分析することによって周波数スペクトルの特徴量を抽出するステップS46。具体的には、近似部421は、回帰分析によって一次式の傾きa、切片bおよび中心周波数fMIDにおける強度cを算出する。図17に示す直線L3は、このステップS46で周波数スペクトル曲線C3に対して特徴量抽出処理を行うことによって得られる回帰直線切片b3である。ステップS47~S50の処理は、図3のステップS7~S10の処理に順次対応している。以上説明した本発明の実施の形態2によれば、複数の既知検体に対してそれぞれ抽出される複数の特徴量の少なくとも一部を座標成分とする第1座標系を有する特徴量空間上において、複数の既知検体を各既知検体の情報に基づいて分類した複数のグループをそれぞれ代表する複数の代表点を所定の軸へ射影したときに、第1座標系の所定の座標軸に沿って隣接する代表点間の距離の和が大きくなる新座標軸を一つの座標軸とする第2座標系における検体の特徴量の座標値を算出し、算出した座標値に応じた表示パラメータを付与するため、グループの違いを明確に峻別することができる。したがって、検体を精度よく観測することができるとともに、観測結果の信頼性を向上させることができる。また、本実施の形態2によれば、第2座標系における検体の特徴量の座標値を算出し、座標値に応じた表示パラメータを付与し、付与した表示パラメータによって定まる画素値を有する特徴量画像データを生成して表示するため、利用者は画像中におけるグループの違いを明確に認識することができる。また、本実施の形態2によれば、周波数スペクトルに対して減衰補正を施した後、特徴量を抽出しているため、超音波の伝播に伴う減衰の影響を取り除くことができ、一段と高精度な観測を行うことができる。ここまで、本発明を実施するための形態を説明してきたが、本発明は上述した実施の形態1、2によってのみ限定されるべきものではない。例えば、特徴量空間上で代表点の分布に偏りがある場合には、γ補正を施してもよい。図18は、本発明の別な実施の形態に係る超音波観測装置が行うγ補正処理の概要を模式的に示す図である。図18において、点Sp1は検体点であり、点μ1、ν1、ρ1は、複数の既知検体を特徴量空間上でグループ分けした時の各グループの代表点である。図18の上側の図に示すように、各グループの代表点μ1、ν1、ρ1のうち、代表点μ1と代表点ν1とのh軸方向の距離は、代表点ν1と代表点ρ1とのh軸方向の距離よりも小さい。そこで、図18に示す場合、超音波観測装置は、特徴量空間上の点に対し、特徴量hが小さい領域の分布をばらつかせるようなγ補正プラスのγ補正を施す。これにより、図18の下側の図に示すように、代表点μ’1と代表点ν’1とのh軸方向の距離が大きくなる一方、代表点ν’1と代表点ρ’1とのh軸方向の距離が小さくなり、3つの代表点がh軸方向に沿って略均一にばらついた分布が得られる。このようにして、超音波観測装置が適切なγ補正処理を行うことにより、特徴量空間における複数のグループの分布に偏りがあるような場合であっても、グループの違いを明確に峻別することが可能となる。なお、ここでは特徴量にガンマ補正を施す場合を説明したが、予めh軸にγ補正を施すようにしてもよい。また、h軸方向の成分のみならず、v軸方向の成分にγ補正を施してもよい。"