"体内の測定領域を加熱することが可能な加熱ビームを照射する加熱源と、測定領域に超音波信号を送波するとともに測定領域からの超音波エコー信号を受波する超音波送受機構と、加熱ビームを照射していないときの測定領域から受波した光非照射時超音波エコー信号と光照射後の測定領域から受波した光照射後超音波エコー信号とに基づいて測定領域内の各点における加熱ビーム照射前後の超音波速度変化を算出する超音波速度変化解析部と、算出された超音波速度変化の分布を輝度情報または色情報にして表示装置の画面上に断層画像として表示する超音波速度変化画像表示制御部と、前記断層画像の一部を関心領域として指定されることを待ち受ける関心領域指定部と、指定された関心領域内の輝度情報または色情報に基づいて、負の超音波速度変化を示す脂肪領域と正の超音波速度変化を示す正常領域とのヒストグラムを算出するヒストグラム算出部と、算出された脂肪領域および正常領域のヒストグラムから脂肪組織の割合の指標となる脂肪化指標を算出する脂肪化指標算出部とを備えたことを特徴とする脂肪組織画像表示装置。関心領域指定部は、関心領域として指定する領域の面積を一定にして指定する請求項1に記載の脂肪組織画像表示装置。脂肪化指標算出部は、脂肪化指標として、脂肪領域と正常領域とのヒストグラムの面積比を算出する請求項1または請求項2に記載の脂肪化像表示装置。超音波速度変化画像表示制御部は、画面上に表示する輝度情報または色情報は多段階にランク分けして表示し、脂肪化指標算出部は、脂肪化指標を算出する際に、ランク分けに基づいて重み付けを加味した脂肪化指標を算出する請求項1または請求項2に記載の脂肪組織画像表示装置。さらに、断層画像上で関心のある任意の直線位置を指定する関心ライン指定部と、指定された直線上の各点の超音波速度変化をグラフ表示する関心ライングラフ表示制御部とを備えた請求項1~請求項4のいずれかに記載の脂肪組織画像表示装置。特許文献2に開示された脂肪組織検出装置によれば、超音波速度変化画像上で、脂肪組織が含まれる領域を、脂肪が含まれない領域と明確に区別させて表示することができる。しかも、上述した脂肪組織検出装置自体は、従来の超音波診断装置の骨格構造を基本にしており、簡単な加熱源と、超音波速度変化の解析機能を追加するだけでよいので、従来の超音波診断装置に、後から追加することができるオプション機構として提供することもできる。しかしながら、生体での実際の診断においては、脂肪組織が含まれる領域の位置を特定するだけでは必ずしも十分ではなく、その領域内での脂肪組織の割合がどの程度であるかを客観的に評価診断できることが望まれている。例えば、肝臓を撮影した画像を用いて脂肪肝の有無を診断する際に、超音波速度変化画像を用いることで、肝臓内での脂肪組織が存在する領域を特定できるが、特定された領域に含まれる脂肪組織の割合脂肪肝細胞と正常肝細胞の割合までは測定できないため、脂肪肝の進行度合いを客観的に判断することは困難である。一般に、脂肪肝の場合、正常肝細胞と脂肪肝細胞とが混在する状態霜降り状態になっており、脂肪部分だけが偏在しているのではない。そのため、超音波速度変化画像では、脂肪が含まれる領域を画面上で、輝度差で表示すなわちグレースケールすると、白と黒との領域が混在した状態が表示されることになる。また、正常肝細胞を赤く、脂肪肝細胞を青く色分けして表示すると脂肪が存在する部分は青と赤とが混在した状態になる。このように、輝度表示や色分けによって脂肪肝の位置は画像上で特定できるが、進行した脂肪肝脂肪の割合が多いと、進行していない脂肪肝脂肪の割合が少ないとを画像で比較したとき、いずれも白青と黒赤とが混在しているように見えるので、脂肪肝の進行度合いを画像だけから客観的に把握することは困難であった。すなわち、超音波速度変化画像の目視による観察では、脂肪含有量の大小を定性的、主観的、感覚的に把握するため、診断する者の熟練度に左右されてしまうことになっていた。そのため、診断する者の熟練度に左右されず、定量的、客観的、合理的に安定して診断結果が得られることが望まれていた。例えば、被検体の脂肪含有量の定期的に計測し、その経時変化を確認する際に、数値的に把握できることが望まれていた。そこで、本発明は、超音波速度変化画像を用いて脂肪肝のような脂肪組織を画像表示する際に、脂肪化の進行度合いとなる指標を数値化して表示することができるようにした脂肪組織画像表示装置を提供することを目的とする。上記課題を解決するために本発明の脂肪組織画像表示装置は、体内の測定領域を加熱することが可能な加熱ビームを照射する加熱源と、測定領域に超音波信号を送波するとともに測定領域からの超音波エコー信号を受波する超音波送受機構と、加熱ビームを照射していないときの測定領域から受波した光非照射時超音波エコー信号と光照射後の測定領域から受波した光照射後超音波エコー信号とに基づいて測定領域内の各点における加熱ビーム照射前後の超音波速度変化を算出する超音波速度変化解析部と、算出された超音波速度変化の分布を輝度情報または色情報にして表示装置の画面上に断層画像として表示する超音波速度変化画像表示制御部と、前記断層画像の一部を関心領域として指定されることを待ち受ける関心領域指定部と、指定された関心領域内の輝度情報または色情報に基づいて負の超音波速度変化を示す脂肪領域と正の超音波速度変化を示す正常領域とのヒストグラムを算出するヒストグラム算出部と、算出された脂肪領域および正常領域のヒストグラムから脂肪組織の割合の指標となる脂肪化指標を算出する脂肪化指標算出部とを備えるようにしている。本発明では、超音波速度変化の分布を輝度情報または色情報にした超音波速度変化分布の断層画像を画面上に表示して、診断者によって関心領域ROIが指定されるのを待ち受ける。関心領域が指定されると、関心領域内の各点の超音波速度変化データから正の超音波速度変化データと負の超音波速度変化データとのヒストグラムを分けて算出する。そして正の超音波速度変化データは正常領域、負の超音波速度変化データは脂肪領域として扱い、これらを用いて脂肪組織の割合の指標となる脂肪化指標を算出する。これにより、算出された脂肪化指標を表示することで、関心領域における脂肪化の進行度合いを数値的に定量表示することができるので、脂肪領域の客観的な評価を行うことができる。これまで、画像上で脂肪組織を着色表示したりして脂肪含有量を把握しやすくしていたが、霜降り状態になっている脂肪組織では、脂肪の存在は確認できるものの、定量化が困難であった。本発明により、霜降り状態であっても脂肪含有量を数値的に出力できるので、客観的な診断を容易に実現できる。特に被検体を定期的に繰り返し計測することで、脂肪含有量の経時変化を数値変化として把握できるようになり、メタボリック症候群の診断に威力を発揮することができる。その他の課題を解決するための手段および効果上記発明において、関心領域指定部は、関心領域として指定する領域の面積を一定にして指定するようにしてもよい。これにより、測定する位置あるいは測定する時期が異なる複数の関心領域どうしを比較するときに、同じ面積で関心領域どうしを比較することができるので、脂肪化の進行度合いの比較を正確に行うことができる。上記発明において、脂肪化指標算出部は、脂肪化指標として、脂肪領域と正常領域とのヒストグラムの面積比を算出するようにしてもよい。すなわち、正の速度変化を示すヒストグラムの面積をS1、負の速度変化を示すヒストグラムの面積をS2として、脂肪化指標をF=S2/S1として定義すると、脂肪の割合が増大するとともに指標Fが大きくなる指標を得ることができ、簡単に進行度合いを数値表現することができるようになる。上記発明において、超音波速度変化画像表示制御部は、画面上に表示する輝度情報または色情報は多段階にランク分けして表示し、脂肪化指標算出部は、脂肪化指標を算出する際に、ランク分けに基づいて重み付けした脂肪化指標を算出するようにしてもよい。輝度情報や色情報の違いは各点の脂肪化の部分的な進行度合いを示すと考えられるので、輝度情報や色情報を脂肪化の進行度合いを表す係数として面積比の計算に加味させることで脂肪化の度合いをより正確に知ることができる。上記発明において、さらに、断層画像上で関心のある任意の直線位置を指定する関心ライン指定部と、指定された直線上の各点の超音波速度変化をグラフ表示する関心ライングラフ表示制御部とを備えるようにしてもよい。指定した直線に沿って超音波速度変化をグラフ表示することができるので、関心領域として指定する位置を定める前に中心位置を探索したりするときに参考にすることができる。以下、本発明の脂肪組織画像表示装置の一実施形態について図面に基づいて説明する。ここでは、加熱手段として赤外線レーザによる近赤外光を照射することにより、測定部位を加熱するものとする。超音波速度変化の測定原理最初に脂肪組織画像表示装置で表示する超音波速度変化画像の測定原理について説明する。この超音波速度変化のデータは以下の関係から求められる。図14は、光非照射時超音波エコー信号と光照射後超音波エコー信号とを示す模式図である。光非照射時の超音波速度をV、光照射後の超音波速度をV’とする。また、光非照射時にある境界間を超音波信号が伝播するときの時間をτとし、同じ境界間距離一定を光照射後に超音波信号の伝播時間をτ-Δτとする。すなわち、温度変化によりΔτだけパルス間隔が短くなるようにシフトしたとする。このとき、Vτ=V’τ-Δτ1の関係が成立し、したがって、2つのエコー信号におけるパルス間隔の時間変化から超音波速度変化は次式で算出することができる。V’/V=τ/τ-Δτ2実際の超音波速度変化データの計算方法についての具体的な説明は省略するが、測定した超音波エコー信号を細かく区分し、各区分に対し、相互相関法により照射前後の波形シフト量を算出することにより、各区分の超音波速度変化データを取得することで各点の超音波速度変化データを算出するようにしている。既述のように、温度変化に対する超音波速度変化を比較すると、水は2m/秒°Cであり、脂肪は-4m/秒°Cであるので、水分が多く含まれる筋肉や内臓腸、腎臓は温度が上がると超音波速度が速くなるのに対し、脂肪部分では超音波速度が減少することになり、超音波速度変化の極性が反転する。算出された超音波速度変化データから、光照射後に超音波速度が負の変化を示す領域を脂肪領域として検出し、正の変化を示す領域を正常領域として検出する。すなわち、式2で求められた超音波速度比が1より大きいか小さいかを判断することにより、1より小さければ、その領域は温度変化に対する超音波速度変化が負になる脂肪領域として検出する。脂肪組織画像表示装置の構成次に、脂肪組織画像表示装置の構成について説明する。図1は、本発明の一実施形態である脂肪組織画像表示装置の構成を示すブロック図である。脂肪組織画像表示装置1は、主に、アレイ型探触子2および赤外線レーザ光源3からなるプローブ4と、アレイ型探触子2を駆動する送受信部5と、制御系10と、DSC11デジタルスキャンコンバータと、表示装置液晶パネル12と、入力装置13キーボード13aおよびマウス13bとを備えている。制御系10は、CPU、メモリを有するコンピュータ装置からなり、メモリに蓄積されたソフトウェアプログラムをロードすることにより制御動作が作動し、脂肪組織画像表示装置1の各部を制御するとともに、得られた測定データの演算を行うようになる。これらの制御、演算を機能ブロックごとに分けて説明すると、照射制御部21、走査制御部22、Bモード信号処理部23、超音波速度変化解析部24、超音波速度変化画像表示制御部25、関心領域指定部26、ヒストグラム算出部27、脂肪化指標算出部28、関心ライン指定部29、関心ライングラフ表示制御部30からなる。プローブ4はアレイ型探触子2と赤外線レーザ光源3とが一体に固定してあり、被検体100の測定部位110に向けて超音波とレーザとを照射することができるようにしてある。アレイ型探触子2は、一方向に配列された複数の振動子を有しており、各振動子は、送受信部5からの駆動信号により励振されて超音波信号を発し、この超音波信号に対する被検体内からの超音波エコー信号を送受信部5に返す。赤外線レーザ光源3は、加熱源として700nm~1000nmの近赤外光を照射する。なお、照射光の波長を変更したり、光源の種類を変更したりすることにより、加熱する測定部位の表皮からの深さを調整することができる。照射制御部21は、予め設定したレーザ出力、設定した時間幅で、レーザ光源3から被検体にレーザ照射を行う制御を行う。レーザ出力および時間幅で測定部位の温度変化幅を調整することができる。走査制御部22は、送受信部5の振動子の駆動時期を制御して、アレイ型探触子2から複数例えば345本の走査線となる超音波信号を送受する制御を行う。照射制御部21と走査制御部22とは、超音波速度変化画像を得るとき、互いに連動して制御する。すなわち、まず、走査制御部22がレーザ照射を行う前の測定部位の超音波エコー信号を取得する制御を行い、続いて照射制御部22が測定部位にレーザを照射して加熱し、続いて再び走査制御部22がレーザ照射後の測定部位の超音波エコー信号を取得する制御を行う。被検体100からの受信信号は、送受信部5を介して取り込まれてメモリに蓄積され、Bモード信号処理部23、超音波速度解析部24に送られ、超音波によるBモード画像や超音波速度変化画像の形成に利用される。Bモード信号処理部23は、受信信号に対して周知のBモード断層画像形成処理を行って輝度情報に変換してビーム走査範囲の断層画像を形成し、DSC11に書き込む。図2は多数例えば345本の超音波ビームの受信信号超音波エコー信号から形成された鶏肉のBモード画像の一例を示す図であり、表示装置12に表示される。この画像からは隣接する組織の境界が明暗で示されるが、脂肪か正常かの差異の判断はつけにくい。超音波速度変化解析部24は、レーザ照射前と照射後の受信信号超音波エコー信号に基づいて超音波ビームを走査した範囲の超音波速度変化の分布の画像を形成する処理を行う。すなわち、式2の演算処理を行い、走査範囲の各点の超音波速度変化データを算出する。算出された超音波速度変化データは、超音波速度変化画像標示制御部25に送られる。超音波速度変化画像標示制御部25は、各点について速度変化データの変化幅の大きさによりランク分けし例えば正負の速度変化をそれぞれ5つの諧調にランク分けする、ランクごとに対応する色情報に変換する。このようにして色情報で正の速度変化を呈する部分と、負の速度変化を呈する部分とを分ける。そして、走査範囲の各点における超音波速度変化を色情報輝度情報として表した分布を、断層画像にしてDSC11に書き込む。DSC11は表示装置12に超音波速度変化の断層画像を表示する。例えば図3は、図2にBモード画像を示した鶏肉についての超音波速度変化画像を示した図であり、青が濃いほど超音波速度が負側に大きく変化していることを示し、赤が濃いほど超音波速度が正側に大きく変化していることを示している。関心領域指定部26は、表示装置12上に表示された断層画像の一部領域を、操作者が入力装置13を用いて関心領域ROIとして指定する入力がなされるのを待ち受け、入力を受け付けると、その指定された領域に対応する色情報輝度情報を、ヒストグラム算出部27に送る制御を行う。具体的には超音波速度変化の断層画像を表示した画面の端に、関心領域を指定するための入力を促す表示例えば関心領域を指定してくださいとの表示を行うを行い、操作者がマウス13b等を用いて画面上の一部領域をカーソルで指定するように促す。音声で指定を促してもよい。指定するカーソルは正方形、円、長方形のいずれかが好ましい。なお、このときカーソルで指定する関心領域の面積基準面積を、メモリに記憶させておき、指定が常に同じ基準カーソル26aになるようにしてもよい。指定される面積が一定になるようにすることにより、常に同じ面積での関心領域の指定が行われるようになり、異なる測定部位、あるいは、異なる測定時期の関心領域どうしを比較することが容易になる。他の測定データと比較する必要がないときは、関心領域として指定するカーソルの面積を、操作者が自由に設定できるようにしてもよい。ヒストグラム算出部27は、指定された関心領域ROI内の色情報または輝度情報に基づいて負の超音波速度変化を示す脂肪領域と正の超音波速度変化を示す正常領域とのヒストグラムを算出する演算を行う。具体的には、超音波速度変化が負に変化した画素と、正に変化した画素とを分離してヒストグラムを算出する。ここでは正の速度変化を5階調、負の速度変化を5階調に分離してヒストグラムを作成する。脂肪化指標算出部28は、正の速度変化をする5つの諧調の画素の総数と、負の速度変化をする5つの諧調の画素の総数との比を算出する。この場合は正常領域と脂肪領域との面積比28aを、脂肪化指標F=S2/S1として算出することになる。あるいは、正と負の各諧調のランクに応じて重み付けを行い、速度変化が大きいランクの重みを大きくするようにして算出した重み付け面積比28bを算出する。この場合は、正と負との各諧調の各ランクの重みをLi,Mii=1~5とすると脂肪化指標F=ΣLiS2/ΣMiS1を算出することになる。そして脂肪化指標Fの算出結果を、表示装置12に脂肪化指標値Fとして表示することもできる。また、関心ライン指定部29は、操作者が表示装置12上に表示された断層画像上で、特に関心のある任意の直線位置関心ラインを指定する入力がなされるのを待ち受け、入力を受け付けると、指定された直線上の各点について超音波速度変化の大きさを抽出し、関心ライングラフ表示制御部30に送る制御を行う。具体的には、上述した関心領域指定部26が超音波変化の断層画像を表示して関心領域を指定する入力を促す表示を行うときに、オプションとして、関心ラインを指定することもできる表示を行うことで、操作者がマウス13b等を用いて画面上の一部分を直線で指定するように促す。関心ライングラフ表示制御部30は、指定された関心ライン上の各点の超音波速度変化データをDSC11に書き込み、表示装置12にグラフ表示する制御を行う。画面表示例次に脂肪組織画像表示装置1によって表示される画面の表示例について説明する。図4は関心領域ROIを指定するときの画面を示す図である。超音波速度変化画像が表示された画面上で基準カーソル26aが表示されるので、所望の位置に移動して関心領域ROIを指定する。これにより、指定された領域内の各画素の色情報が抽出され、ヒストグラムが算出されることになる。生体における実際の脂肪組織は、正常領域と脂肪領域とが混在する霜降り状態になっているが、説明の便宜上、以下に示す種々の極端なダミーサンプルを測定したときに表示されるヒストグラムと算出される脂肪化指標Fについて説明する。1正常領域のみの場合図5は、基準カーソルで指定した関心領域ROI内の全体が、脂肪の含まれていない正常領域のみからなるダミーサンプルについて測定した場合のヒストグラムである。この場合は面積S1で正の速度変化のみが現れるヒストグラムになる。脂肪化指標F=S2/S1は、S2が0なので、0となる。2脂肪領域のみの場合図6は、基準カーソルで指定した関心領域ROI内の全体が、脂肪領域のみからなるダミーサンプルについて測定した場合のヒストグラムである。この場合は面積S2で負の速度変化のみが現れるヒストグラムになる。脂肪化指標F=S2/S1は、S1が0なので、無限大∞となる。3正常領域と脂肪領域とが半分ずつ細かく均等分散する場合図7は基準カーソルで指定した関心領域ROI内に、脂肪領域と正常領域とが細かく混在するダミーサンプルについて測定した場合のヒストグラムである。この場合は、速度変化0近傍で、面積S1で正の速度変化と、面積S2で負の速度変化とがほぼ半々のヒストグラムになる。脂肪化指標F=S2/S1は1近傍の値になる。4正常領域と脂肪領域とが半分ずつ細かく均等分散し、さらに脂肪が混在する場合図8は基準カーソルで指定した関心領域ROI内に、脂肪領域と正常領域とが細かく混在し、さらに脂肪が増えて点在するようになったダミーサンプルについて測定した場合のヒストグラムである。この場合は面積S1で正の速度変化と、面積S2ただしS1<S2で負の速度変化とが表れるヒストグラムになる。脂肪化指標F=S2/S1は1よりも少し大きな値になる。5正常領域と脂肪領域とが二分されている場合図9は、基準カーソルで指定した関心領域ROI内に、脂肪領域と正常領域とが2つに分かれて存在する場合のヒストグラムである。粗い霜降り状態が見られるときがこの状態に近いと考えられる。この場合は面積S1で正の速度変化と、面積S2で負の速度変化とが分離して表れるヒストグラムになる。脂肪化指標F=S2/S1は1近傍の値になる。6ほとんど正常領域とほとんど脂肪領域とが二分されている場合図10は、基準カーソルで指定した関心領域ROI内に、ほとんどが脂肪である領域とほとんどが正常である領域とが、2つに分かれて存在する場合のヒストグラムである。粗い霜降り状態から少し霜降りが細かくなった状態に近いと考えられる。この場合は面積S1で正の速度変化と、面積S2ただしS1<S2で負の速度変化とが分離して表れるヒストグラムになる。上述した5の場合よりもヒストグラムの形状がブロードになる。脂肪化指標F=S2/S1は1近傍の値になる。これらのダミーサンプルで類推されるように、正常領域だけの極端な場合から脂肪だけの極端な場合までを脂肪化指標Fで数値化して表現することができる。なお、上述した実施形態では脂肪化指標Fに面積比=S2/S1を用いたが、これに代えて階調ごとの重み付けを加味させた脂肪化指標F=ΣLiS2/ΣMiS1を採用してもよい。関心ラインの表示例次に、関心ラインによる超音波速度変化のグラフ表示について説明する。図11は、関心ラインを指定するときの画面を示す図である。超音波速度変化画像が表示された画面上で関心のあるラインの始点と終点とを指定すると、関心ラインLが画面上に表示される。そして、図12に示すように、指定された直線上の各点の超音波速度がグラフGとして表示される。このグラフ表示を参照して、超音波速度変化のライン状での大小を見やすくすることで、基準カーソルの位置を定めるときに、参考にすることができる。なお、図11では、関心ラインは1本のみ表示しているが、複数本表示させてもよいし、水平方向だけでなく垂直方向あるいは斜め方向に関心ラインを表示させてもよい。さらに、関心ラインの他の応用方法として、定期的に、同じ位置の関心ラインに沿った超音波速度変化のグラフ表示を取得しておき、比較するようにしてもよい。これにより、関心ラインに沿って脂肪含有量の経時的な変化を容易に把握することができるようになる。以上、本発明の一実施形態について説明したが、本発明は必ずしも上記の実施形態に限定されるものでなく、その目的を達成し、請求の範囲を逸脱しない範囲内で適宜修正、変更することが可能である。例えば、上記実施形態では測定部位を加熱するための加熱源として赤外線レーザ光源を用いたが、これに代えて、超音波送波器を用いた超音波加熱やマイクロ波などの電磁波による加熱を行うようにしてもよい。"