(57)【特許請求の範囲】
【請求項１】
被検体に対する放射線撮影により得られる放射線画像を解析する画像解析装置であって、
前記放射線画像に含まれる画素値の階調を複数の階調区間に区分けするときの各々の階調幅を設定する階調幅設定手段と、
画素値が互いに所定の階調区間に属している２つの画素のペアであって画素同士が所定の距離だけ離間しているものが前記放射線画像の一部に設けられた解析範囲に何回現れるかを階調区間の組み合わせごとに数えて同時生起行列を生成する同時生起行列生成手段と、
前記同時生起行列に基づいてテクスチャ解析指標を算出する指標算出手段とを備え、
前記階調幅設定手段は、前記放射線画像に含まれるノイズ成分の増加に応じて広くなるように前記階調幅を設定し、
前記階調幅設定手段で設定された階調幅によって、前記同時生起行列を生成する際の階調区間が決定されることを特徴とする画像解析装置。
【請求項２】
請求項１に記載の画像解析装置において、
前記階調幅設定手段が前記被検体の体厚に基づいて前記ノイズ成分を推定することを特徴とする画像解析装置。
【請求項３】
請求項２に記載の画像解析装置において、
術者に前記体厚を入力させる入力手段を備えることを特徴とする画像解析装置。
【請求項４】
請求項２に記載の画像解析装置において、
前記放射線画像に写り込む被検体像を構成する画素値のバラツキの大きさに基づいて前記体厚を推定することを特徴とする画像解析装置。
【請求項５】
請求項２ないし請求項４のいずれかに記載の画像解析装置において、
前記階調幅設定手段が参照する前記階調幅と前記ノイズ成分との関連性を記憶する記憶手段を備え、
前記記憶手段が記憶する関連性は、ある体厚の被検体を写し込んだ放射線画像をテクスチャ解析すると得られるテクスチャ解析指標が基準となる放射線画像についてのテクスチャ解析指標と同じとなるように定められていることを特徴とする画像解析装置。
【請求項６】
請求項１ないし請求項５のいずれかに記載の画像解析装置によって算出された前記テクスチャ解析指標に基づいて、前記被検体の骨梁の指標を算出することを特徴とする骨梁解析装置。
【請求項７】
請求項１ないし請求項６のいずれかに記載の画像解析装置において、
前記指標算出手段が算出する前記テクスチャ解析指標として、コリレーション、ディシミラレィティ、コントラスト、ホモジェネイティ、エントロピー、アングラーセカンドモーメント、バリアンス、インバースディファレンシャルモーメントのうちの１つまたは複数が選択されていることを特徴とする画像解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体が写り込んでいる画像をテクスチャ解析してテクスチャ解析指標を算出する画像解析装置に関する。
【背景技術】
【０００２】
医療機関には、被検体に放射線撮影をして得られた被検体像に対して画像解析を行って、被検体の状態を評価する評価値を算出する画像解析装置が備えられている。
【０００３】
画像解析を行う際に問題となるのが画像に写り込む被検体の体厚である。放射線には、体厚が厚いほど被検体を透過しにくくなる性質がある。したがって、被検体の体厚が厚いほど撮影が露光不足となり、得られる断層画像のＳ／Ｎ比が悪くなる。断層画像のＳ／Ｎ比が悪くなると、それだけ画像にノイズ成分が多く現れる。従って、画像解析を行うには、被検体の体厚についての補正を行う必要がある（例えば特許文献１参照）。
【０００４】
画像解析の手法の一つとしてテクスチャ解析がある。このテクスチャ解析は、画像を同時生起行列に変換して、同時生起行列を演算することにより評価値の一種であるテクスチャ解析指標を算出するというものである。同時生起行列は、画像を構成する画素値がどのようにバラついているかの統計を表しており、画像上の被検体像の状態によって変化する行列となっている（例えば、非特許文献１，非特許文献２参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００１−８６４０９号公報
【非特許文献１】Ｈａｒａｌｉｃｋ  ＲＭ．  Ｓｔａｔｉｓｔｉｃａｌ  ａｎｄ  ｓｔｒｕｃｔｕｒａｌ  ａｐｐｒｏａｃｈｅｓ  ｔｏ  ｔｅｘｔｕｒｅ．  Ｐｒｏｃ  ＩＥＥＥ  １９７９；６７（５）：７８６−８０４．
【非特許文献２】Ｈａｒａｌｉｃｋ  ＲＭ．  ｅｔ  ａｌ．  Ｔｅｘｔｕｒａｌ  Ｆｅａｔｕｒｅｓ  ｆｏｒ  Ｉｍａｇｅ  Ｃｌａｓｓｉｆｉｃａｔｉｏｎ．  ＩＥＥＥ    Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｓｙｓｔｅｍｓ  Ｍａｎ  ａｎｄ  Ｃｙｂｅｒｎｅｔｉｃｓ  １９７３；６：６１０−６２１．
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来構成に係る画像解析装置には次のような問題点がある。
すなわち、従来構成に係る画像解析装置によれば、テクスチャ解析を用いた画像解析において被検体の体厚の影響を十分に配慮していない。
【０００７】
従来構成の画像解析装置を説明するのに、ここでは、骨梁解析を例にとることにする。骨梁とは、骨の内部の海綿質を構成する細長状の構造である。従来構成の骨梁解析の対象となるのは、被検体に放射線撮影を行うことで得られた被検体の断層画像である。
【０００８】
テクスチャ解析を用いた骨梁解析では、骨梁についての情報のみを表したテクスチャ解析指標を算出するのが理想である。すなわち、骨梁解析のテクスチャ解析指標は、本来は骨梁の状態に応じて変動するのであり、他の条件によっては変動しないはずなのである。
【０００９】
しかし、実際の断層画像に骨梁解析を行ってみると、被検体の体厚が厚くなるのに応じて算出されるテクスチャ解析指標が変動してしまう。このような現象が起こる理由は、被検体の体厚が応じて断層画像のノイズ成分の量が変動することにある。ノイズ成分は、被検体の骨梁の状態を表すものではない。にもかかわらず、テクスチャ解析では、ノイズ成分が断層画像を構成する画素値のバラつきの一種であるとしてしまいテクスチャ解析指標の算出に加味する。テクスチャ解析は、断層画像におけるノイズ成分の含有量の変化を捉えてテクスチャ解析指標に反映させてしまう。
【００１０】
この体厚の影響をテクスチャ解析指標に及ぼさないようなテクスチャ解析の方法があれば、より正確なテクスチャ解析指標が算出できる。しかしながら、従来の画像解析装置によれば、体厚の違う被検体に対して同じテクスチャ解析を実行するようになっており、この様な配慮はなされていない。果たしてテクスチャ解析をどのように改変すれば体厚の影響を防げるのであろうか。なお、従来より知られている体厚の補正は、テクスチャ解析についてのものではない。
【００１１】
本発明は、この様な事情に鑑みてなされたものであって、その目的は、被検体の体厚に影響されずに適切なテクスチャ解析指標を算出するテクスチャ解析を実現できる画像解析装置を提供することにある。
【課題を解決するための手段】
【００１２】
本発明は上述の課題を解決するために次のような構成をとる。
すなわち、本発明に係る画像解析装置は、被検体に対する放射線撮影により得られる放射線画像を解析する画像解析装置であって、放射線画像に含まれる画素値の階調を複数の階調区間に区分けするときの各々の階調幅を設定する階調幅設定手段と、画素値が互いに所定の階調区間に属している２つの画素のペアであって画素同士が所定の距離だけ離間しているものが放射線画像の一部に設けられた解析範囲に何回現れるかを階調区間の組み合わせごとに数えて同時生起行列を生成する同時生起行列生成手段と、同時生起行列に基づいてテクスチャ解析指標を算出する指標算出手段とを備え、階調幅設定手段は、放射線画像に含まれるノイズ成分の増加に応じて広くなるように階調幅を設定することを特徴とするものである。
【００１３】
［作用・効果］本発明の構成によれば、被検体の体厚の違いが解析結果に影響しないテクスチャ解析の手法を具体的に示したものとなっている。すなわち、放射線画像に含まれるノイズ成分が増えるに応じて広くなるように階調幅を設定しておく。この階調幅は、同時生起行列生成時において同一視する階調の幅を表している。つまり、本発明によれば放射線画像に含まれるノイズ成分が増えると同時生起行列を構成する要素が少なくなるように工夫がされている。放射線画像のノイズ成分が増えると、それだけ放射線画像における画素値のバラつきが大きくなる。このような画像をテクスチャ解析するに当たり、従来構成のように同時生起行列の要素を多くしたままだと、ノイズ成分に由来する画素値のバラつきも同時生起行列に表現されるので、テクスチャ解析指標はノイズ成分を加味した値となってしまう。
本発明の構成のように、画像のノイズ成分が多くなるほど同時生起行列を構成する要素を少なくするようにすると、ノイズ成分に由来する画素値のバラつきの影響が解析途中の同時生起行列生成時に吸収されて、テクスチャ解析指標には現れない。したがって、本発明の画像解析装置によれば、より信頼性の高いテクスチャ解析指標を算出することができるのである。
【００１４】
また、上述の画像解析装置において、階調幅設定手段が被検体の体厚に基づいてノイズ成分を推定すればより望ましい。
【００１５】
［作用・効果］上述の構成は、本発明の画像解析装置をより具体的に示したものとなっている。放射線画像に含まれるノイズ成分の量は、被検体の体厚に相関して変動する。被検体の体厚によってノイズ成分の量を推定することができるのである。
【００１６】
また、上述の画像解析装置において、術者に体厚を入力させる入力手段を備えればより望ましい。
【００１７】
［作用・効果］上述の構成は、本発明の画像解析装置をより具体的に示したものとなっている。ノイズ成分の量を推定するのに用いる被検体の体厚を術者に入力させるような構成とすれば、より確実に体厚を取得することができるようになる。
【００１８】
また、上述の画像解析装置において、放射線画像に写り込んだ被検体像を構成する画素値のバラツキの大きさに基づいて、体厚を推定すればより望ましい。
【００１９】
［作用・効果］上述の構成は、本発明の画像解析装置をより具体的に示したものとなっている。被検体像を構成する画素の間で見られる画素値のバラツキの強度により体厚を推定すれば、より確実にノイズ成分の量を算出することができる。
【００２０】
また、上述の画像解析装置において、階調幅設定手段が参照する階調幅とノイズ成分との関連性を記憶する記憶手段を備え、記憶手段が記憶する関連性は、ある体厚の被検体を写し込んだ放射線画像をテクスチャ解析すると得られるテクスチャ解析指標が基準となる放射線画像についてのテクスチャ解析指標と同じとなるように定められていればより望ましい。
【００２１】
［作用・効果］上述の構成は、本発明の画像解析装置をより具体的に示したものとなっている。ある体厚の被検体を写し込んだ放射線画像をテクスチャ解析すると、得られるテクスチャ解析指標が基準となる放射線画像についてのテクスチャ解析指標と同じとなるように定められた階調幅とノイズ成分との関連性により階調幅設定手段が動作すれば、階調幅設定手段は、実測結果に基づいて動作することになるので、より信頼性の高いテクスチャ解析指標が取得できる。
【００２２】
また、上述の画像解析装置によって算出されたテクスチャ解析指標に基づいて、被検体の骨梁の指標を算出すればより望ましい。
【００２３】
［作用・効果］上述のように、本発明の構成は具体的には骨梁解析に適用することができる。
【００２４】
また、上述の骨梁解析装置において、指標算出手段が算出するテクスチャ解析指標として、コリレーション、ディシミラレィティ、コントラスト、ホモジェネイティ、エントロピー、アングラーセカンドモーメント、バリアンス、インバースディファレンシャルモーメントのうちの１つまたは複数が選択されていればより望ましい。
【００２５】
［作用・効果］上述の構成は、本発明の画像解析装置をより具体的に示したものとなっている。本発明の構成は多様なテクスチャ解析指標について適用することができる。
【発明の効果】
【００２６】
本発明の構成によれば、被検体の体厚の違いが解析結果に影響しないテクスチャ解析の手法を具体的に示したものとなっている。すなわち、被検体の体厚が増えるに応じて広くなるように階調幅を設定しておく。この階調幅は、同時生起行列生成時において同一視する階調の幅を表している。つまり、本発明によれば被検体の体厚が増えると同時生起行列を構成する要素が少なくなるように工夫がされている。このようにすると、ノイズ成分に由来する画素値のバラつきの影響が解析途中の同時生起行列生成時に吸収され、より信頼性の高いテクスチャ解析指標を算出することができる。
【図面の簡単な説明】
【００２７】
【図１】実施例１に係る画像解析装置の全体構成を説明する機能ブロック図である。
【図２】実施例１に係る背景処理を説明する模式図である。
【図３】実施例１に係る体厚値取得処理を説明する模式図である。
【図４】実施例１に係るテーブルを説明する模式図である。
【図５】実施例１に係る同時生起行列の生成処理を説明する模式図である。
【図６】実施例１に係る同時生起行列の生成処理を説明する模式図である。
【図７】実施例１に係る同時生起行列の生成処理を説明する模式図である。
【図８】実施例１に係る同時生起行列の生成処理を説明する模式図である。
【図９】実施例１に係る階調幅と体厚値の関係を示す模式図である。
【図１０】実施例１に係る構成の効果が生じる原理を説明する模式図である。
【図１１】実施例１に係る構成の効果の実証方法を説明する模式図である。
【図１２】実施例１に係る構成の効果の実証方法を説明する模式図である。
【図１３】実施例１に係る構成の効果の実証方法を説明する模式図である。
【図１４】実施例１に係る構成の効果の実証方法を説明する模式図である。
【図１５】実施例１に係る構成の効果を説明する模式図である。
【実施例１】
【００２８】
次に、本発明に係る画像解析装置の実施例について図面を参照しながら説明する。なお、本発明の画像解析処理装置の処理対象となる画像の例としてフィルタバックプロジェクション法により取得される被検体Ｍの断層画像を取り上げるものとする。断層画像は、本発明の放射線画像に相当する。実施例の画像解析装置は、入力された断層画像に対して画像処理を施して骨梁解析に係るテクスチャ解析の指標値を算出する。
【００２９】
図１は、本発明に係る画像解析装置の全体構成を示した機能ブロック図である。以下、図１に示された各部の動作の詳細について順を追って説明する。
【００３０】
＜背景処理部１１の動作＞
断層画像Ｄは、まず背景処理部１１に送出される。背景処理部１１は、断層画像Ｄにアンシャープネスマスキング処理などの背景処理を施し、背景処理がされた断層画像Ｄａを生成する。このような背景処理により異なる断層画像Ｄ同士で見られる撮影条件の相違による被検体像の写り込みの違いを消去することができる。背景処理部１１が行う背景処理は、アンシャープネスマスキング処理に限られず他の処理としてもよい。
【００３１】
図２は、アンシャープネスマスキング処理を簡単に説明している。アンシャープネスマスキング処理とは、断層画像Ｄから断層画像Ｄを平滑化して得られた平滑化画像を差し引く画像処理である。この様な画像処理により、断層画像Ｄは、よりくっきりした断層画像Ｄａに変換される。
【００３２】
＜バラツキ算出部１２の動作＞
断層画像Ｄａは、バラツキ算出部１２に送出される。バラツキ算出部１２は、断層画像Ｄａを解析して断層画像Ｄａに写り込む被検体Ｍの体厚を示す標準偏差ＳＤを算出する。図３は、バラツキ算出部１２が断層画像Ｄａの解析を行っている様子を示している。バラツキ算出部１２は、まず、術者が入力部１３を通じて設定した解析領域Ｒを認識する。解析領域Ｒは、画像における被検体像のうち骨が写り込んでいない部分（被検体像の軟組織が写り込んでいる部分）が解析領域Ｒの全域となるように選択される。したがって、解析領域Ｒには、被検体像の軟組織が写り込んでいる。入力部１３は、術者に入力値を入力させる構成であり、本発明の入力手段に相当する。
【００３３】
そして、バラツキ算出部１２は、認識した解析領域Ｒを構成する画素に対して統計的な解析を行って標準偏差ＳＤを算出しこれを被検体Ｍの体厚を示す体厚値とする。このときのバラツキ算出部１２は、断層画像Ｄａに写り込んだ被検体像を構成する画素値のバラつきの大きさ（標準偏差ＳＤ）に基づいて、体厚を推定しているともいえる。
【００３４】
体厚値（標準偏差ＳＤ）は、解析領域Ｒに属する画素における画素値のバラつきを示しており、被検体Ｍの軟組織像を構成する画素の間で見られる画素値のバラツキを示してもいる。標準偏差ＳＤが大きいほど、画素値のバラつきは大きいと評価できる。この標準偏差ＳＤにより、断層画像Ｄａに含まれるノイズ成分の量（ノイズ量）を知ることができる。断層画像Ｄａには、被検体像にノイズ成分が重畳している。このノイズ成分は具体的には砂嵐状のパターンとなる。画像に重畳するノイズ量が多くなると、画像には、砂嵐状のパターンがより強く表れる。つまり、ノイズ量が多いと、それだけ、画素の画素値のバラつきが大きくなるわけである。
【００３５】
つまり、体厚値（標準偏差ＳＤ）は、画像におけるノイズ量を示す指標ともなっている。すなわち、標準偏差ＳＤが大きいほどノイズ量は多いと評価ができるのである。一方、解析領域Ｒに写り込む被検体Ｍの軟組織像は、標準偏差ＳＤにさほど寄与しない。軟組織像は、互いに似通った画素値から構成されるからである。
【００３６】
解析領域Ｒを被検体Ｍの軟組織に設定したのには、次のような理由がある。すなわち、バラツキ算出部１２が算出する標準偏差ＳＤは、被検体Ｍの体厚によって変化するパラメータである。すなわち、標準偏差ＳＤが大きいと、それだけ断層画像Ｄは露光不足で撮影されたことになる。この露光不足は、被検体Ｍの体厚が厚いことによって引き起こされるから、標準偏差ＳＤは、間接的に被検体Ｍの体厚を表していることになるわけである。
【００３７】
ところで、被検体Ｍの骨部は、もともとＸ線を透過させにくいので多くのノイズ成分を含んでいる。したがって、被検体Ｍの骨部における画素値の標準偏差ＳＤは、被検体Ｍの体厚によらず大きな値をとる。したがって、体厚値を調べるには、被検体Ｍの軟組織についての標準偏差を算出したほうがよいのである。
【００３８】
＜階調幅設定部１４の動作＞
バラツキ算出部１２は、断層画像Ｄａにおける体厚値（標準偏差ＳＤ）を階調幅設定部１４に送出する。階調幅設定部１４は、送出された標準偏差ＳＤに対応する幅値Ｗａを図４に示すようなテーブルに基づいて取得し、取得した幅値Ｗａを階調幅Ｗに設定する。テーブルにおける各幅値Ｗａは、標準偏差ＳＤが増えると例えば線形的に増加するように規定されている。したがって、階調幅設定部１４は、断層画像Ｄａに含まれる画素値の階調を複数の階調区間の各々に区分けするときの各々の階調幅Ｗをノイズ成分の増加に応じて広くなるように設定することになる。この際、階調幅設定部１４は、体厚値（標準偏差ＳＤ）に基づいて放射線画像に含まれるノイズ成分の量の推定をしている。
【００３９】
階調幅設定部１４は、送出された標準偏差ＳＤがテーブルにない場合は、テーブルに補完処理を加えて新たな幅値Ｗａを取得するように構成される。また、階調幅設定部１４は、標準偏差ＳＤと幅値Ｗａとが関連した方程式に基づいて階調幅Ｗの設定を実行するようにしてもよい。上述のテーブルおよび方程式は、記憶部２０に記憶されている。階調幅設定部１４が設定した階調幅Ｗは、行列生成部１５に送出される。テーブルがどのように生成されるのかについては後述のものとする。階調幅設定部１４は、本発明の階調幅設定手段に相当し、行列生成部１５は、本発明の同時生起行列生成手段に相当する。
【００４０】
＜行列生成部１５の動作＞
テクスチャ解析を行う際に必要となる行列として同時生起行列（ＧＬＣＭ）がある。この行列は行列生成部１５により生成される。背景処理部１１が生成した断層画像Ｄａは、行列生成部１５に送出され、そこでＧＬＣＭに変換される。図５は、行列生成部１５が断層画像Ｄａに基づいてＧＬＣＭを生成する動作を説明している。図５の左側は、断層画像Ｄａを画素値の２次元配列として表している。説明の簡単のため、断層画像Ｄａを構成する各画素の画素値は、０から９までの１０通りの値をとるものとする。
【００４１】
図５の右側に示すように、断層画像Ｄａより生成されるＧＬＣＭの行数と列数は、いずれも画素の画素値がとりえる画素値の数と一致する。断層画像Ｄａを構成する各画素は、１０通りのうちのいずれかの画素値を有しているのであるから、断層画像Ｄａより生成されるＧＬＣＭは１０行１０列の２次元行列となる。行列生成部１５は、１０×１０行列となっているＧＬＣＭを構成する１００個の要素に数値を代入してＧＬＣＭを完成させる。各要素にどのような数値を入れるかは、断層画像Ｄａの画素値に基づいて判断される。
【００４２】
図５は、ＧＬＣＭの各行のうち０を意味する行、各列のうち１を意味する列に位置する要素ｐ（０，１）の数値を行列生成部１５が決めようとしているところを示している。行列生成部１５は、画素値０と画素値１とが隣り合って配列されている画素のペアが断層画像Ｄａに何組あるかを数えて、そのカウント数をＧＬＣＭの要素ｐ（０，１）とする。図５においては、画素値０と画素値１とが隣り合って配列されている画素のペアは２組あるので、要素ｐ（０，１）の値は、２となる。このＧＬＣＭにおける任意の要素ｐ（ａ，ｂ）は要素ｐ（ｂ，ａ）に等しいので、ＧＬＣＭの要素ｐ（１，０）の値も２となる。
【００４３】
行列生成部１５は、同様な動作をＧＬＣＭの全域に亘って行い、行列が有する要素の全てを断層画像Ｄａに基づいて決定する。こうして行列生成部１５は、断層画像Ｄａに基づいてＧＬＣＭを完成させる。
【００４４】
図６は、行列生成部１５が断層画像Ｄａに基づいてＧＬＣＭを生成する様子を示している。生成されるＧＬＣＭは、断層画像Ｄａの画素が取り得る画素値の数が多くなるほど大きくなる。ＧＬＣＭは、対称性を有する行列であり、図６の点線で示す対角線で２つ折りにすると、重なり合う要素同士の値が同じとなっているような行列である。
【００４５】
＜階調幅Ｗに基づいた動作の変更＞
続いて、本発明の最も特徴的な構成である階調幅設定部１４が設定した階調幅Ｗに基づいて行列生成部１５の動作が変更される様子について説明する。すなわち、行列生成部１５は、設定された階調幅Ｗに基づいて、生成するＧＬＣＭを変化させる。本発明に係るＧＬＣＭを構成する行列の要素は、階調幅Ｗに基づいて増減するように構成されている。
【００４６】
図７は、行列生成部１５が生成するＧＬＣＭの一種を示している。図５におけるＧＬＣＭは、０から９までの階調をすべて区別して生成されているが、図７のＧＬＣＭは、０と１を区別せずに生成されたものである。同様に、図７のＧＬＣＭは、２と３，４と５，６と７，８と９のそれぞれを区別してはいない。このＧＬＣＭは５行５列の２次元行列となる。行列生成部１５は、５×５行列となっているＧＬＣＭを構成する２５個の要素に数値を代入してＧＬＣＭを完成させる。
【００４７】
図７は、ＧＬＣＭの各行のうち０から１までの階調の区間（階調区間）を意味する行、各列のうち０から１までの階調区間を意味する列に位置する要素ｐ（０〜１，０〜１）の数値を行列生成部１５が決めようとしているところを示している。行列生成部１５は、一方が階調区間（０〜１）に属する画素値を有し、もう一方も階調区間（０〜１）に属する画素値を有する画素のペアが断層画像Ｄａに何組あるかを数えて、そのカウント数をＧＬＣＭの要素ｐ（０〜１，０〜１）とする。図７においては、画素値０と画素値１とが隣り合って配列されている画素のペアが２組、画素値１と画素値１とが隣り合って配列されている画素のペアが１組あるので、要素ｐ（０〜１，０〜１）の値は、３となる。
【００４８】
図８は、行列生成部１５が断層画像Ｄａに基づいてＧＬＣＭを生成する様子を示している。生成されるＧＬＣＭは、図６で説明したＧＬＣＭと比べて構成要素が少ない。
【００４９】
階調幅Ｗは、行列生成部１５がＧＬＣＭを生成するときに、同一視される階調区間の幅を示している。例えば、図５における行列生成部１５は、階調のすべてを区別していたので、このＧＬＣＭ生成動作における階調幅Ｗは１となる。また、図７における行列生成部１５は、２つの階調を同一視していたので、このＧＬＣＭ生成動作における階調幅Ｗは２となる。階調幅Ｗは、行列生成部１５がＧＬＣＭを生成する前に階調幅設定部１４が設定している。
【００５０】
行列生成部１５がＧＬＣＭを生成するときの階調幅Ｗは、断層画像Ｄａの体厚値（標準偏差ＳＤ）が増加すると大きくなるように設定されている。したがって、行列生成部１５が生成するＧＬＣＭは、図９に示すように断層画像Ｄａの標準偏差ＳＤに応じて変化する。断層画像Ｄａの標準偏差ＳＤが小さい場合は、図９の上段に示すように、階調幅Ｗが小さいので階調区間が多く設定され、生成されるＧＬＣＭは要素が多い行列となる。断層画像Ｄａの標準偏差ＳＤが大きくなってくると、階調幅Ｗは次第に広くなっていき、設定される階調区間の数が減少し、図９の中段、下段に示すように、生成されるＧＬＣＭは要素が次第に少ない行列となる。したがって、行列生成部１５は、断層画像Ｄａの標準偏差ＳＤの増加に応じてＧＬＣＭの規模が小さくなるように動作する。
【００５１】
行列生成部１５は、階調幅設定部１４から階調幅Ｗを受信して、断層画像Ｄａの階調を階調幅Ｗ単位で区切り、複数の階調区間を設定する。これにより、行列生成部１５は、画素値が互いに所定の階調区間に属している２つの画素のペアであって画素同士が所定の距離だけ離間しているものが断層画像Ｄの一部に設けられた解析範囲に何回現れるかを階調区間の組み合わせごとに数えてＧＬＣＭ（同時生起行列）を生成するわけである。解析範囲としては、骨の海綿骨の部分が選択される。
【００５２】
＜断層画像Ｄａの体厚値に応じてＧＬＣＭの規模を変化させる効果＞
図１０は、行列生成部１５が生成するＧＬＣＭの規模を可変にしたことによる効果について説明している。図１０の上段は、断層画像Ｄａの一部について階調幅Ｗを１としてＧＬＣＭを生成するとどうなるかを示している。階調幅Ｗは１だから、階調はすべて区別されてＧＬＣＭが生成される。断層画像Ｄａのうち図１０の上段に示す縦３×横３の領域は、ＧＬＣＭの生成時にどのような扱いを受けるかを考える。この領域には、縦方向に隣接する２つの画素のペアが６つあり、横方向に隣接する２つの画素のペアが６つある。そして、右斜め方向に隣接する２つの画素のペアが４つあり、左斜め方向に隣接する２つの画素のペアが４つある。これら２０組のペアは、すべて画素値が５００の画素の組み合わせだから、ＧＬＣＭを生成するときにすべて同じ扱いを受けることになる。すなわち、２０組のペアはすべてＧＬＣＭの要素ｐ（５００，５００）のカウント数に加味される。
【００５３】
図１０の中段は、図１０の上段に示す断層画像の一部にノイズ成分が重畳している状態を表している。ノイズ成分の重畳によりすべて５００だった画素値は、ノイズ成分によりバラつき、４９４となったり５０６となったりする。このようなノイズ成分が重畳した断層画像の一部について階調幅Ｗを１としてＧＬＣＭを生成するとどうなるか。領域上の２０組のペアにおける画素値の組み合わせは１つとして同じものがないから、２０組のペアはＧＬＣＭを生成するときにすべて別の扱いを受けることになる。たとえば、画素値が５０６，４９４の組み合わせとなっている画素のペアはＧＬＣＭの要素ｐ（５０６，４９４）のカウント数に加味され、残り１９組のペアは、この要素ｐ（５０６，４９４）のカウント数に加味されない。
【００５４】
以上のことから、ノイズ成分が断層画像に重畳することにより、生成されるＧＬＣＭの要素に影響が表れてしまうことがわかる。ノイズ成分により断層画像が乱れると、行列生成部１５のカウント動作に変化が現れ、ＧＬＣＭも乱されてしまうのである。
【００５５】
図１０の下段も中断と同様、図１０の上段に示す断層画像の一部にノイズ成分が重畳している状態を表している。このようなノイズ成分が重畳した断層画像の一部について階調幅Ｗを２０としてＧＬＣＭを生成するとどうなるか。領域上の２０組のペアにおける画素値の組み合わせは確かに１つとして同じものがない。しかし、階調幅Ｗは２０なので、行列生成部１５は、例えば４９０〜５１０までの画素値を同一視してＧＬＣＭを生成する。したがって、これら２０組のペアは、ＧＬＣＭを生成するときにすべて同じ扱いを受けることになる。すなわち、この例では、２０組のペアはすべてＧＬＣＭの要素ｐ（４９０〜５１０，４９０〜５１０）のカウント数に加味される。
【００５６】
以上のことから、階調幅Ｗを調整することで、生成されるＧＬＣＭの要素にノイズ成分の影響が表れてしまうことが防止できることが分かる。本来同じ扱いを受けるはずの図１０の下段に示す２０組のペアは、階調幅Ｗが調整されたことにより、別の扱いを受けることなくすべて同じ扱いを受けてＧＬＣＭの生成に加味されたからである。
【００５７】
断層画像Ｄａのノイズ成分が多いほどＧＬＣＭは大きく乱される。このようなＧＬＣＭの大きな乱れを防止するには、階調幅Ｗを広げる必要がある。断層画像Ｄａの体厚値（標準偏差ＳＤ）が大きくなると、画素値の乱れは、さらに大きくなる。すなわち、図１０の中段、下段の説明では、画素値のバラつきは、２０の階調幅Ｗに収まるものであったが、標準偏差ＳＤが大きくなると、画素値のこの幅に収まりきれなくなる。そこで、本発明によれば、標準偏差ＳＤの増加に伴ってＧＬＣＭ生成時の階調幅Ｗを広くするようにしている。このようにすれば、画素値が大きくばらついても、ＧＬＣＭ生成時において画素値の同一視される範囲が広がって、バラつきの影響が吸収される。結果として、ＧＬＣＭには断層画像Ｄａのノイズ成分の影響が現れない。
【００５８】
＜テクスチャ解析指標算出部１６の動作＞
ＧＬＣＭは、テクスチャ解析指標算出部１６に送出される。テクスチャ解析指標算出部１６は、ＧＬＣＭに種々の演算を実行することでテクスチャ解析指標を算出することが可能である。テクスチャ解析指標算出部１６が算出できるテクスチャ解析指標は、例えば次のようなものがある。式中のｐ（ｉ，ｊ）とは、ＧＬＣＭにおけるｉ行ｊ列目の要素の値、Σｉ，Σｊは、それぞれｉ行、ｊ列についての要素の合計、Ｎｇは、断層画像Ｄａの画素が取り得る画素値の数、μは平均値、μｘ，μｙは、それぞれ行方向、列方向の平均値、σｘ，σｙは、それぞれ行方向、列方向の標準偏差を表している。なお、これらテクスチャ解析指標ＡＳＭ（Ａｎｇｕｌａｒ  Ｓｅｃｏｎｄ  Ｍｏｍｅｎｔ：アングラーセカンドモーメント），ＣＮＴ（Ｃｏｎｔｒａｓｔ：コントラスト），ＣＯＲ（Ｃｏｒｒｅｌａｔｉｏｎ：コリレーション），ＶＡＲ（Ｖａｒｉａｎｃｅ：バリアンス），ＩＤＭ（Ｉｎｖｅｒｓｅ  Ｄｉｆｆｅｒｅｎｃｅ  Ｍｏｍｅｎｔ，インバースディファレンシャルモーメント），ＥＮＴ（Ｅｎｔｒｏｐｙ，エントロピー）の各々は、１９７３年にＨａｒｌｉｃｋらが上述の非特許文献２で提唱した１４種類のパラメータのうちの一部である。また、ＤＩＳは非類似度またはディシミラレィティと呼ばれるテクスチャ解析指標で、ＨＯＭは、一様性またはホモジェネイティと呼ばれるテクスチャ解析指標である。テクスチャ解析指標算出部１６は、本発明の指標算出手段に相当する。
【００５９】
【数１】
0006299504.tif 000002
【００６０】
テクスチャ解析指標算出部１６は、ＧＬＣＭに上述の種々の演算を行ってテクスチャ解析指標を算出する。テクスチャ解析指標算出部１６が算出するテクスチャ解析指標の種類と数は、適宜変更することができる。以上のように、テクスチャ解析指標算出部１６は、ＧＬＣＭ（同時生起行列）に基づいてテクスチャ解析を行いテクスチャ解析指標の算出を実行する。テクスチャ解析指標は、例えば、骨梁の状態変化を推測するのに用いられる。この場合、テクスチャ解析指標は、撮影時期の異なる同じ被検体Ｍの断層画像Ｄの各々について算出されることになる。
【００６１】
本発明を構成する各部１１，１２，１４，１５，１６は、ＣＰＵが各種のプログラムを実行することで実現される。なお、各部は、各部の機能を個別に実行する演算装置によって実現するようにしてもよい。
【００６２】
＜本発明の効果の実証＞
続いて、本発明の効果の実証を実際に行ったので、その結果について説明する。本発明の効果を理解するには、まずは断層画像Ｄの撮影方法を簡単に知っておく必要があるのでこれについて説明する。
【００６３】
図１１は、断層画像Ｄの生成に係るフィルタバックプロジェクション法の原理を説明している。例えば、天板２に平行な（鉛直方向に対して水平な）仮想平面（基準裁断面ＭＡ）について説明すると、図２に示すように、基準裁断面ＭＡに位置する点Ｐ，Ｑが、常にＦＰＤ４のＸ線検出面の不動点ｐ，ｑのそれぞれに投影されるように、Ｘ線管３によるコーン状のＸ線ビームＢの照射方向に合わせてＦＰＤ４をＸ線管３の反対方向に同期移動させながら一連のＸ線画像が生成される。
【００６４】
一連のＸ線画像には、被検体Ｍの投影像が位置を変えながら写り込んでいる。そして、この一連のＸ線画像を再構成すれば、基準裁断面ＭＡに位置する像（たとえば、不動点ｐ，ｑ）が集積され、Ｘ線断層画像としてイメージングされることになる。
【００６５】
一方、基準裁断面ＭＡに位置しない点Ｉは、ＦＰＤ４における投影位置を変化させながら一連の被検体Ｍ画像に点ｉとして写り込んでいる。この様な点ｉは、不動点ｐ，ｑとは異なり、Ｘ線投影画像を重ね合わせる段階で像を結ばずにボケる。このように、一連の投影画像の重ね合わせを行うことにより、被検体Ｍの基準裁断面ＭＡに位置する像のみが写り込んだＸ線断層画像が得られる。このように、投影画像を単純に重ね合わせると、基準裁断面ＭＡにおける断層画像が得られる。
【００６６】
さらに、フィルタバックプロジェクション法によれば、基準裁断面ＭＡに水平な任意の裁断面においても、同様な断層画像を得ることができる。撮影中、ＦＰＤ４において上記点ｉの投影位置は移動するが、投影前の点Ｉと基準裁断面ＭＡとの離間距離が大きくなるにしたがって、この移動速度は増加する。これを利用して、取得された一連の被検体Ｍ画像を所定のピッチで体軸方向Ａにずらしながら再構成を行うようにすれば、基準裁断面ＭＡに平行な裁断面における断層画像が得られるわけである。
【００６７】
つまり、断層画像Ｄは、撮影方向が異なるＸ線画像を重ね合わせた画像となっており、Ｘ線画像の各々が露光不足であると、断層画像Ｄも露光不足となる。このようなＸ線画像の各々が露光不足となる原因は、被検体Ｍの体厚が撮影時の想定よりも厚いことにある。撮影対象の被検体Ｍの体厚が厚いとそれだけＸ線が通過しにくくなるので被検体Ｍを通過してＦＰＤ４に検出されるＸ線の線量は少なくなる。これが、Ｘ線画像の各々が露光不足になる理由であり、断層画像Ｄが露光不足になる理由ともなっている。このような露光不足の断層画像には露光が適切な場合の断層画像と比べてより多くのノイズ成分を含んでしまっている。
【００６８】
＜テーブルの生成＞
本発明の効果を確認するには、図４で説明したテーブルを生成する必要がある。このテーブルは、断層画像に含まれる体厚値（標準偏差ＳＤ）と階調幅Ｗとが関連したものである。したがって、テーブルを生成するには被検体Ｍの体厚の条件が異なる複数の断層画像を用意する必要がある。この断層画像は、次のような方法で用意することができる。まず図１２の左側に示すように、天板２に被検体Ｍに見立てたファントムＰｈを載置し、図１１で説明した原理に基づいて断層画像Ｄ０を撮影する。
【００６９】
続いて、図１２の右側に示すように、ファントムＰｈの上側にアクリル板Ａｃを載置し、この状態で再び断層画像Ｄ１を撮影する。このアクリル板Ａｃは、透過しようとするＸ線の一部を透過させないので、被検体Ｍの体厚が厚くなったのと同じような効果を生み出すことができる。以降、図１３左側に示すように、ファントムＰｈの上側に置かれるアクリル板Ａｃの枚数を増やしながら断層画像Ｄ２，Ｄ３，Ｄ４を撮影していけば、体厚が異なる被検体Ｍをそれぞれ撮影したかのような断層画像を得ることができる。こうして得られた断層画像Ｄ０，Ｄ１，Ｄ２，Ｄ３，Ｄ４は、同じファントムＰｈを写し込んではいるものの、体厚値（標準偏差ＳＤ）が異なる。
【００７０】
アクリル板ＡｃをファントムＰｈの上に置かない状態で撮影した断層画像Ｄ０を用いて従来通りのテクスチャ解析を行うと、何らかのテクスチャ解析指標が算出される。一方、アクリル板ＡｃをファントムＰｈの上に置いた状態で撮影した断層画像Ｄ１を用いて従来通りのテクスチャ解析を行ってもテクスチャ解析指標は算出できる。ただし、断層画像Ｄ０で得られたテクスチャ解析指標と、断層画像Ｄ１で得られたテクスチャ解析指標とは一致しない。断層画像Ｄ０，Ｄ１により算出される体厚値（標準偏差ＳＤ）が異なり、この影響によりテクスチャ解析指標の値に食い違いが生じてしまうからである。従来通りのテクスチャ解析とは、一定の階調幅Ｗに基づいて同時生起行列を生成する解析をいう。
【００７１】
ここで、断層画像Ｄ１についてテクスチャ解析を施す時に同時生起行列を生成するときの階調幅Ｗを広くしてもう一度テクスチャ解析指標を生成してみる。すると、テクスチャ解析指標は、階調幅Ｗを広くしたことにより断層画像Ｄ０についてのテクスチャ解析指標の値に近づく。階調幅Ｗが広いとテクスチャ解析指標値が断層画像Ｄ１のノイズ成分の影響を受けなくなるからである。
【００７２】
階調幅Ｗを次第に広くしながら同じアクリル板Ａｃありの断層画像Ｄ１についてテクスチャ解析を繰り返していくと、テクスチャ解析指標がアクリル板Ａｃなしの断層画像Ｄ０についてのテクスチャ解析指標と同じになるときがくる。この時の階調幅Ｗは、当該断層画像Ｄ１のテクスチャ解析に適切な階調幅Ｗとなる。なお、階調幅Ｗをこれ以上広くすると、テクスチャ解析の手法が断層画像Ｄ０のものと不必要にかけ離れてしまう。このような事態は、テクスチャ解析指標の比較を行う上で問題を生じる可能性がある。テクスチャ解析指標の値同士が同じとは、断層画像Ｄ１に係るテクスチャ解析指標から断層画像Ｄ０に係るテクスチャ解析指標を除算して得られる数値が十分に１に近い状態をいう。
【００７３】
写り込むアクリル板Ａｃの枚数が異なる断層画像Ｄ２，Ｄ３，Ｄ４の各々について、上述した階調幅Ｗの調整を行えば、断層画像Ｄ２，Ｄ３，Ｄ４の各々について適切な階調幅Ｗが算出される。
【００７４】
ところで、断層画像Ｄ０，Ｄ１，Ｄ２，Ｄ３，Ｄ４について標準偏差ＳＤを個別に算出することができる。そこで、図１３の右側に示すように一方の軸に標準偏差、もう一方の軸に適切な階調幅Ｗをとった座標系に各断層画像Ｄ０，Ｄ１，Ｄ２，Ｄ３，Ｄ４をプロットすることができるはずである。これらプロットは、直線をなす近似式上に整列していると考えることができる。図４におけるテーブルが示す幅値Ｗａと標準偏差ＳＤの関連性は、実は、このような算出方法により得られたものである。
【００７５】
このように、テーブルが示す標準偏差ＳＤと階調幅Ｗとの関連性は、ある体厚の被検体を写し込み所定の標準偏差ＳＤを有する断層画像Ｄをテクスチャ解析すると、得られるテクスチャ解析指標が基準となる断層画像Ｄ０についてのテクスチャ解析指標と同じとなるように定められている。
【００７６】
今回の実証では、各断層画像のプロットに基づき、Ｗ＝０．０７７７ＳＤ−２．５９１５（Ｒ２＝０．９９８３）という近似式を得ることができた。テクスチャ解析指標のうち算出の対象となったのは、エントロピーである。なお、近似式を取得する際に算出された標準偏差ＳＤおよび適切な階調幅Ｗは、断層画像全体から算出された。
【００７７】
＜実際のテクスチャ解析への利用＞
このようにして得られた幅値Ｗａと標準偏差ＳＤの関連性に基づいて、実際にテクスチャ解析を行ったのでこの結果について説明する。テクスチャ解析は、図１４に示す大腿骨の異なる３つの部位について行われた。３つの部位とは、具体的には、主抗圧骨梁の分布域、主抗張骨梁の分布域、およびｗａｒｄ三角である。主抗圧骨梁とは圧力が強くかかる骨梁のことであり、主抗張骨梁とは張力が強くかかる骨梁のことである。ｗａｒｄ三角とは、大腿骨において骨梁が空いている特定の部位のことである。
【００７８】
図１５の左側は、本発明に依らない従来の方式で断層画像Ｄにテクスチャ解析を行った結果を示している。従来の方式とは、同時生起行列の生成時における階調の幅を一定に定める方式のことである。図１５の左側は、主抗圧骨梁、主抗張骨梁、ｗａｒｄ三角についてのテクスチャ解析指標を６枚の異なる断層画像Ｄの間で比較をしている。テクスチャ解析指標のうち算出の対象となったのは、エントロピーである。この算出の対象は、上述のテーブルの作成で用いられたものと同じテクスチャ解析指標とするのが望ましい。６枚の断層画像Ｄには、アクリル板Ａｃを設けない状態で撮影されたものが１枚、ファントムＰｈの上側に載置するアクリル板Ａｃの厚さが１ｃｍ，２ｃｍ，３ｃｍ，４ｃｍ，５ｃｍのそれぞれで撮影されたものが５枚ある。
【００７９】
撮影時のアクリル板Ａｃの条件が異なる６枚の断層画像Ｄに写り込んでいるファントムＰｈはすべて同じものであるから、テクスチャ解析指標は、どの断層画像Ｄにおいても同じ値となるのが理想である。この点、図１５の左側における主抗圧骨梁については、ほぼ理想通りとなっている。したがって、主抗圧骨梁について、従来の解析方法で得られるテクスチャ解析指標の信用性は高いとは言える。
【００８０】
しかし、主抗張骨梁、ｗａｒｄ三角については、断層画像に写り込むアクリル板Ａｃの枚数に依存してテクスチャ解析指標が変動してしまっている。したがって、主抗張骨梁、ｗａｒｄ三角について、従来の解析方法で得られるテクスチャ解析指標の信用性は高いとは言えない。被検体Ｍの体厚の違いによりテクスチャ解析指標が変動してしまっているからである。テクスチャ解析指標の変動は、６枚の断層画像Ｄの間で体厚値（標準偏差ＳＤ）が異なることで生じていると考えられる。
【００８１】
図１５の右側は、本発明の方式に従って断層画像Ｄにテクスチャ解析を行った結果を示している。本発明の方式とは、断層画像Ｄの体厚値（標準偏差ＳＤ）に応じて階調の幅を変える方式のことである。解析に使用した断層画像Ｄは、図１５の左側に係る解析と同じ６枚の画像であり、テクスチャ解析指標のうち算出の対象となったのは、エントロピーである。なお、図１５の左側における実証では、断層画像Ｄの被検体Ｍの軟組織像を解析して得られる標準偏差ＳＤとが線形的な関係となるように階調の幅を決定している。
【００８２】
図１５の右側によれば、本発明の方式の適用により主抗張骨梁、ｗａｒｄ三角について、撮影時のアクリルの条件に依存したテクスチャ解析指標の変動が抑制されていることがわかる。したがって、主抗張骨梁、ｗａｒｄ三角について、本発明に係る解析方法で得られるテクスチャ解析指標の信用性は高いとは言える。被検体Ｍの体厚の違いによりテクスチャ解析指標が変動しないからである。
【００８３】
一方、主抗圧骨梁については、上述の通り従来の解析方法で十分に信頼性の高いテクスチャ解析指標は得られている。この点、本発明の方式に従って主抗圧骨梁についてのテクスチャ解析を行ってみたところ、得られたテクスチャ解析指標が断層画像Ｄの間で大きく変動することはなかった。つまり、従来の解析方法に代えて本発明の解析方法を用いて主抗圧骨梁に対するテクスチャ解析を行ったとしても、テクスチャ解析指標の信頼性は低下しないというわけである。
【００８４】
これらのことからすると、従来の解析方法よりも本発明に係る解析方法を用いたほうが高い信頼性のテクスチャ解析指標が算出できるということが言える。
【００８５】
以上のように、本発明の構成によれば、被検体Ｍの体厚の違いが解析結果に影響しないテクスチャ解析の手法を具体的に示したものとなっている。すなわち、被検体Ｍの体厚が増えてノイズ成分が増えるに応じて広くなるように階調幅Ｗを設定しておく。この階調幅Ｗは、ＧＬＣＭ生成時において同一視する階調の幅を表している。つまり、本発明によれば被検体Ｍの体厚が増え、断層画像Ｄａに含まれるノイズ成分が増えるとＧＬＣＭを構成する要素が少なくなるように工夫がされている。断層画像Ｄのノイズ成分が増えると、それだけ断層画像Ｄにおける画素値のバラつきが大きくなる。このような画像をテクスチャ解析するに当たり、従来構成のようにＧＬＣＭの要素を多くしたままだと、ノイズ成分に由来する画素値のバラつきもＧＬＣＭに表現されるので、テクスチャ解析指標はノイズ成分を加味した値となってしまう。
【００８６】
本発明の構成のように、画像のノイズ成分が多くなるほどＧＬＣＭを構成する要素を少なくするようにすると、ノイズ成分に由来する画素値のバラつきの影響が解析途中のＧＬＣＭ生成時に吸収されて、テクスチャ解析指標には現れない。したがって、本発明の画像解析装置によれば、より信頼性の高いテクスチャ解析指標を算出することができるのである。
【００８７】
上述のようにバラツキ算出部１２は、被検体Ｍの軟組織像を構成する画素の間で見られる画素値のバラツキの強度により標準偏差ＳＤを算出すれば、より確実に標準偏差ＳＤを算出することができる。被検体Ｍの軟組織像ではない骨像は、元来Ｘ線を透過しにくい骨の像であることから、体厚に関係なく多くのノイズ成分を含んでしまう。
【００８８】
上述したように、ある体厚の被検体を写し込み所定の標準偏差ＳＤを有する断層画像Ｄをテクスチャ解析すると、得られるテクスチャ解析指標が基準となる断層画像Ｄについてのテクスチャ解析指標と同じとなるように定められた階調幅Ｗと標準偏差ＳＤとの関連性により階調幅設定部１４が動作すれば、階調幅設定部１４は、実測結果に基づいて動作することになるので、より信頼性の高いテクスチャ解析指標が取得できる。  
【００８９】
本発明の構成は、上述のものに限られず、下記のように変形実施が可能である。
【００９０】
（１）実施例１の構成では、標準偏差ＳＤを基準として画像の体厚値を認識していたが、本発明はこの構成に限られない。標準偏差ＳＤの代わりにコントラスト値を用いるようにしてもよい。その他、バラツキ算出部１２は、画素値のバラつきの指標として、標準偏差ＳＤ以外の指標を算出するように構成を変更することができる。
【００９１】
（２）実施例１の構成では、断層画像Ｄを解析することにより画像の体厚値（標準偏差ＳＤ）を認識していたが、本発明はこの構成に限られない。入力部１３を通じて術者より入力された情報または断層画像Ｄを構成する画像ファイルに埋め込まれた情報に基づいて実際の被検体Ｍの体厚または標準偏差ＳＤを認識するようにしてもよい。このときの参照される情報としては、例えば、被検体Ｍの体厚の計測結果を示すものであってもよいし、術者の目分量を示すものであってもよい。また、例えば、被検体ＭのＢＭＩ値であってもよいし、撮影時におけるホトタイマの駆動状況を示すものであってもよい。本変形例によれば、バラツキ算出部１２に代わり、体厚算出部が設けられている。体厚算出部は、上述の各種情報に基づいて被検体Ｍの体厚を算出する。本変形例によれば、体厚算出部は、術者の入力値に基づいて被検体Ｍの体厚を算出する。そして、階調幅設定部１４は、被検体Ｍの体厚と階調幅Ｗとが関連したテーブルを参照しながら動作することになる。
【００９２】
（３）本発明の画像解析装置は、骨梁解析専用というわけではない。肺の気管支像についての解析、肺の構造解析、腫瘍の形状解析などの他の解析についても適用が可能である。
【００９３】
（４）本発明の画像解析装置は、常に同じテーブルを使用していたが、本発明はこの構成に限られない。解析の目的や目的の部位によって固有のテーブルを使用するようにしてもよい。本変形例によれば、例えばｗａｒｄ三角の骨梁解析を行う時に、それ専用のテーブルが規定する体厚値（標準偏差ＳＤ）と幅値Ｗａとの関連性に基づいてテクスチャ解析指標の算出がなされることになる。
【００９４】
（５）上述の説明では断層画像全域を解析することで図４のテーブルを算出していたが、本発明はこの構成に限られない。断層画像の軟組織やファントムが写り込んでいない背景領域などについて画像解析を行うことでテーブルを算出するようにしてもよい。
【００９５】
（６）本発明の画像処理装置は、解析対象の画像に対して単一のテクスチャ解析指標を算出するように構成してもよいし、複数のテクスチャ解析指標を同時に算出するように構成してもよい。
【符号の説明】
【００９６】
ＳＤ   標準偏差（体厚値）
１３   入力部（入力手段）
１４   階調幅設定部（階調幅設定手段）
１５   行列生成部（同時生起行列生成手段）
１６   テクスチャ解析指標算出部（指標算出手段）
【図１】
0006299504.tif 000003
【図２】
0006299504.tif 000004
【図３】
0006299504.tif 000005
【図４】
0006299504.tif 000006
【図５】
0006299504.tif 000007
【図６】
0006299504.tif 000008
【図７】
0006299504.tif 000009
【図８】
0006299504.tif 000010
【図９】
0006299504.tif 000011
【図１０】
0006299504.tif 000012
【図１１】
0006299504.tif 000013
【図１２】
0006299504.tif 000014
【図１３】
0006299504.tif 000015
【図１４】
0006299504.tif 000016
【図１５】
0006299504.tif 000017