(57)【要約】【課題】二次元画像を基にした骨梁解析をより正確に実行できる骨梁解析装置を提供する。【解決手段】断層画像Ｄａに画像解析を施すことで得られるテクスチャ解析指標ｔ１は、骨梁の状態により変動するパラメータではあるものの、骨梁の三次元構造の状態を盛り込んでいないから骨梁の状態を正確に表したものとはならない。そこで、本発明では、断層画像Ｄａで得られたテクスチャ解析指標に基づいて、断層画像Ｄａの被写体が仮に三次元解析されたとしたときに得られる評価値を推定するようにしている。断層画像Ｄａから容易に取得できるテクスチャ解析指標により、従来では行えなかった正確な骨梁の状態を把握できる骨梁解析装置が提供できる。【選択図】図１４
2015107142.tif 000002
【特許請求の範囲】
【請求項１】
生体の骨を写し込んだ二次元画像の一部に設けられた解析範囲を構成する各画素のうち所定の画素値の組み合わせを有する２つの画素のペアで画素同士が所定の距離だけ離間しているものが前記解析範囲において何回現れるかを各画素値の組み合わせごとに数えて同時生起行列を生成する同時生起行列生成手段と、
前記同時生起行列に基づいてテクスチャ解析を行いテクスチャ解析指標を算出する指標算出手段と、
前記テクスチャ解析指標と生体での実測が難しい評価値との関連性に基づいて前記評価値の推定値を算出する評価値推定手段とを備えることを特徴とする骨梁解析装置。
【請求項２】
請求項１に記載の骨梁解析装置において、
前記評価値推定手段が前記推定値の算出に用いる前記テクスチャ解析指標と評価値との関連性は、標本骨における標本解析範囲を解析して得られる前記評価値と前記テクスチャ解析指標とを前記標本骨における前記標本解析範囲を定める位置を変えながら次々と算出することで得られたものであることを特徴とする骨梁解析装置。
【請求項３】
請求項２に記載の骨梁解析装置において、
前記評価値推定手段が前記推定値の算出に用いる前記テクスチャ解析指標と評価値との関連性は、前記標本骨から得られた結果を回帰分析することによって算出される回帰式である推定式として表されており、回帰分析をするときの目的変数は、前記標本解析範囲について三次元解析して得られる前記評価値であり、独立変数は、前記標本解析範囲に係る二次元画像を解析して得られる前記テクスチャ解析指標であることを特徴とする骨梁解析装置。
【請求項４】
請求項１ないし請求項３のいずれかに記載の骨梁解析装置において、
前記評価値推定手段が推定する評価値は、解析範囲内の骨成分の体積を示すＢＶ／ＴＶ値、解析範囲内の骨のつまり具合を示すＴＢＰｆ値、三次元解析により得られる骨梁総延長を表すＴＳＬ値、骨梁の幅を表すＴｂＴｈ値、骨の強度を表す骨強度値のいずれかであることを特徴とする骨梁解析装置。
【請求項５】
請求項１ないし請求項４のいずれかに記載の骨梁解析装置において、
前記指標算出手段が算出するテクスチャ解析指標として、コリレーション、ディシミラレィティ、コントラスト、ホモジェネイティ、エントロピー、アングラーセカンドモーメント、バリアンス、インバースディファレンシャルモーメントのうちの１つまたは複数が選択されていることを特徴とする骨梁解析装置。
【請求項６】
請求項１ないし請求項５のいずれかに記載の骨梁解析装置において、
放射線を照射する放射線源と、
前記放射線源を被検体に対し移動させる放射線源移動手段と、
前記放射線源移動手段を制御する放射線源移動制御手段と、
被検体を透過した放射線を検出する検出手段と、
前記検出手段を被検体に対し移動させる検出器移動手段と、
前記検出器移動手段を制御する検出器移動制御手段と、
前記検出手段の出力を基に画像を生成する画像生成手段と、
前記放射線源および前記検出手段を被検体に対して移動させながら連写された画像を基に断層画像を生成する断層画像生成手段とを備え、
前記同時生起行列生成手段が解析する二次元画像は、断層画像であることを特徴とする骨梁解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体の骨梁を解析する骨梁解析装置に関し、特に被検体の放射線透過画像を取得して放射線透過画像の解析を行うことで骨梁解析を実行する骨梁解析装置に関する。
【背景技術】
【０００２】
骨梁とは、骨の内部の海綿質を構成する細長状の構造である。この骨梁が骨の内部で充実しているかを知ることで被検体Ｍの健康の診断や疾病の診断ができる。また、被検体Ｍの骨梁を解析すれば被検体Ｍの骨強度も知ることができる（例えば、非特許文献１，非特許文献２参照）。
【０００３】
従来の骨梁解析手法について説明する。従来の骨梁解析装置５０は、図２０に示すように、被検体Ｍを載置する天板５２と、天板５２の上側に設けられている放射線源５３と、天板５２の下側に設けられている検出器５４とを備えている。骨梁解析を行うには、図２０の様な装置を用いて、断層撮影を行い、取得された断層画像に骨梁解析が施される。
【０００４】
従来構成における骨梁解析の具体的手法について説明する。従来構成によれば、骨梁解析に先立って、放射線源５３と検出器５４とを被検体Ｍに対して移動させながら撮影を行うことにより、被検体Ｍの断層画像が取得が行われる。そして、取得された断層画像に写り込む骨梁の解析が実行される（例えば特許文献１参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１３−０２７６０８号公報
【非特許文献】
【０００６】
【非特許文献１】Ｈａｒａｌｉｃｋ  ＲＭ．  Ｓｔａｔｉｓｔｉｃａｌ  ａｎｄ  ｓｔｒｕｃｔｕｒａｌ  ａｐｐｒｏａｃｈｅｓ  ｔｏ  ｔｅｘｔｕｒｅ．  Ｐｒｏｃ  ＩＥＥＥ  １９７９；６７（５）：７８６−８０４．
【非特許文献２】Ｈａｒａｌｉｃｋ  ＲＭ．  ｅｔ  ａｌ．  Ｔｅｘｔｕｒａｌ  Ｆｅａｔｕｒｅｓ  ｆｏｒ  Ｉｍａｇｅ  Ｃｌａｓｓｉｆｉｃａｔｉｏｎ．  ＩＥＥＥ    Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｓｙｓｔｅｍｓ  Ｍａｎ  ａｎｄ  Ｃｙｂｅｒｎｅｔｉｃｓ  １９７３；６：６１０−６２１．
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、上述のような従来構成によれば、次のような問題点がある。
すなわち、従来装置では、骨梁解析の方法が最適化がされていない。
【０００８】
従来方法における断層画像に基づいた骨梁解析によれば、被検体の骨の状態を必ずしも正確に表した指標が得られない。被検体の骨を形作る骨梁は、そもそも３次元的なネットワーク構造を有している。したがって本来ならば、被検体の骨梁の状態は、３次元構造を考慮して行わなければ正確に知ることができない。したがって、被検体の骨梁の状態を知ろうとすれば、骨梁の３次元構造を取得するような撮影を被検体に対して行わなければならない。
【０００９】
しかし、実際には、被検体の骨梁の三次元像を正確に取得することは難しい。骨梁は三次元像を撮影するにはあまりにも微細だからである。そこで、従来構成においては、骨梁の三次元像の代わりに骨の断層画像などの２次元画像を骨梁解析に用いるような構成としている。骨の断層画像のほうが骨梁の三次元像よりも撮影しやすいのである。
【００１０】
この様な事情があるので、断層画像を基にして実行される骨梁解析の結果は、被検体の骨の状態を必ずしも正確に表したものとならない。この骨梁解析は、骨梁の三次元構造を考慮せずになされたものであり、実際の骨梁が有する様々な情報のうち断層画像を解析することで分かる情報しか用いずに結果を導くものだからである。
【００１１】
したがって、従来の断層画像を基にして実行される骨梁解析においては、解析方法に改良の余地があるということになる。被検体の骨の状態をより正確に知るには、従来行われている方法よりも多角的な観点を盛り込んで骨梁解析を行う必要がある。しかし、それを具体的にどのように行えばよいのかは、十分に検討されてこなかった。
【００１２】
本発明は、この様な事情に鑑みてなされたものであって、断層画像を基にした骨梁解析をより正確に実行できる骨梁解析装置を提供することにある。
【課題を解決するための手段】
【００１３】
本発明は上述の課題を解決するために次のような構成をとる。
すなわち、本発明に係る骨梁解析装置は、生体の骨を写し込んだ二次元画像の一部に設けられた解析範囲を構成する各画素のうち所定の画素値の組み合わせを有する２つの画素のペアで画素同士が所定の距離だけ離間しているものが解析範囲において何回現れるかを各画素値の組み合わせごとに数えて同時生起行列を生成する同時生起行列生成手段と、同時生起行列に基づいてテクスチャ解析を行いテクスチャ解析指標を算出する指標算出手段と、テクスチャ解析指標と生体での実測が難しい評価値との関連性に基づいて評価値の推定値を算出する評価値推定手段とを備えることを特徴とするものである。
【００１４】
［作用・効果］本発明の骨梁解析装置によれば、二次元画像を基にした骨梁解析でありながらも、骨梁の状態をより正確に把握できる。すなわち、本発明に係る骨梁解析装置は、二次元画像に画像解析を施すことで骨梁の大まかな状態を表す数値を取得する。このときの数値とは、本発明においてはテクスチャ解析指標である。この数値は、骨梁の状態により変動するパラメータではあるものの、骨梁の三次元構造の状態を盛り込んでいないから骨梁の状態を正確に表したものとはならない。そこで、本発明では、この二次元画像で得られたテクスチャ解析指標に基づいて、二次元画像の被写体が仮に三次元解析されたとしたときに得られる評価値を推定するようにしている。二次元画像から容易に取得できるテクスチャ解析指標により、従来では行えなかった正確な骨梁の状態を把握できる骨梁解析装置が提供できる。
【００１５】
また、上述の骨梁解析装置において、評価値推定手段が推定値の算出に用いるテクスチャ解析指標と評価値との関連性は、標本骨における標本解析範囲を解析して得られる評価値とテクスチャ解析指標とを標本骨における標本解析範囲を定める位置を変えながら次々と算出することで得られたものであればより望ましい。
【００１６】
［作用・効果］上述の構成は本発明の骨梁解析装置の具体的構成を表したものとなっている。評価値推定手段が推定値の算出に用いるテクスチャ解析指標と評価値との関連性は、標本骨における標本解析範囲を解析して得られる評価値とテクスチャ解析指標とを標本骨における標本解析範囲を定める位置を変えながら次々と算出することで得られれば、より確実な推定値を求めることができる。
【００１７】
また、上述の骨梁解析装置において、評価値推定手段が推定値の算出に用いるテクスチャ解析指標と評価値との関連性は、標本骨から得られた結果を回帰分析することによって算出される回帰式である推定式として表されており、回帰分析をするときの目的変数は、標本解析範囲について三次元解析して得られる評価値であり、独立変数は、標本解析範囲に係る二次元画像を解析して得られるテクスチャ解析指標であればより望ましい。
【００１８】
［作用・効果］上述の構成は本発明の骨梁解析装置の具体的構成を表したものとなっている。評価値推定手段が推定値の算出に用いるテクスチャ解析指標と評価値との関連性を示す推定式は標本骨について実際に三次元解析および二次元画像解析の両方をして、目的変数を三次元解析に係る評価値とし、独立変数を二次元画像解析に係るテクスチャ解析指標としたときの回帰分析によって得られたものである。この様にすることで、より確実に評価値の推定を実行できる。
【００１９】
また、上述の骨梁解析装置において、評価値推定手段が推定する評価値は、解析範囲内の骨成分の体積を示すＢＶ／ＴＶ値、解析範囲内の骨のつまり具合を示すＴＢＰｆ値、三次元解析により得られる骨梁総延長を表すＴＳＬ値、骨梁の幅を表すＴｂＴｈ値、骨の強度を表す骨強度値のいずれかであればより望ましい。
【００２０】
［作用・効果］上述の構成は本発明の骨梁解析装置の具体的構成を表したものとなっている。評価値推定手段が推定する評価値は、ＢＶ／ＴＶ値、ＴＢＰｆ値、ＴＳＬ値、ＴｂＴｈ値、骨強度値のいずれかであれば、骨梁の評価をするときに必要な評価値が得られる骨梁解析装置が提供できる。
【００２１】
また、上述の骨梁解析装置において、指標算出手段が算出するテクスチャ解析指標として、コリレーション、ディシミラレィティ、コントラスト、ホモジェネイティ、エントロピー、アングラーセカンドモーメント、バリアンス、インバースディファレンシャルモーメントのうちの１つまたは複数が選択されていればより望ましい。
【００２２】
［作用・効果］上述の構成は本発明の骨梁解析装置の具体的構成を表したものとなっている。上述のテクスチャ指標値は、既知のものであり、比較的容易に算出できる。従って、上述の構成によれば、本発明の骨梁解析装置をより確実に実現できる。
【００２３】
また、上述の骨梁解析装置において、放射線を照射する放射線源と、放射線源を被検体に対し移動させる放射線源移動手段と、放射線源移動手段を制御する放射線源移動制御手段と、被検体を透過した放射線を検出する検出手段と、検出手段を被検体に対し移動させる検出器移動手段と、検出器移動手段を制御する検出器移動制御手段と、検出手段の出力を基に画像を生成する画像生成手段と、放射線源および検出手段を被検体に対して移動させながら連写された画像を基に断層画像を生成する断層画像生成手段とを備え、同時生起行列生成手段が解析する二次元画像は、断層画像であればより望ましい。
【００２４】
［作用・効果］上述の構成は本発明の骨梁解析装置の具体的構成を表したものとなっている。本発明は、上述のようなデジタルトモシンセシス装置にも適用できる。
【発明の効果】
【００２５】
本発明の骨梁解析装置によれば、二次元画像を基にした骨梁解析でありながらも、骨梁の状態をより正確に把握できる。二次元画像に画像解析を施すことで得られるテクスチャ解析指標は、骨梁の状態により変動するパラメータではあるものの、骨梁の三次元構造の状態を盛り込んでいないから骨梁の状態を正確に表したものとはならない。そこで、本発明では、二次元画像で得られたテクスチャ解析指標に基づいて、二次元画像の被写体が仮に三次元解析されたとしたときに得られる評価値を推定するようにしている。二次元画像から容易に取得できるテクスチャ解析指標により、従来では行えなかった正確な骨梁の状態を把握できる骨梁解析装置が提供できる。
【図面の簡単な説明】
【００２６】
【図１】実施例１に係る骨梁解析装置の全体構成を説明する機能ブロック図である。
【図２】実施例１に係る断層画像の撮影原理を説明する模式図である。
【図３】実施例１に係る骨梁解析部の詳細を説明する機能ブロック図である。
【図４】実施例１に係るアンシャープネスマスキング処理を説明する模式図である。
【図５】実施例１に係る同時生起行列の生成方法を説明する模式図である。
【図６】実施例１に係る同時生起行列が生成される様子を説明する模式図である。
【図７】実施例１に係るＴＢＰｆ値を説明する模式図である。
【図８】実施例１に係る推定式について説明する模式図である。
【図９】実施例１に係る推定式について説明する模式図である。
【図１０】実施例１に係る推定式について説明する模式図である。
【図１１】実施例１に係るシミュレーションについて説明する模式図である。
【図１２】実施例１に係るシミュレーションについて説明する模式図である。
【図１３】実施例１に係るシミュレーションについて説明する模式図である。
【図１４】実施例１に係る骨梁解析部の動作について説明する模式図である
【図１５】実施例２に係る断層画像の撮影原理を説明する模式図である。
【図１６】実施例２に係る断層画像の撮影原理を説明する模式図である。
【図１７】実施例２に係る断層画像の撮影原理を説明する模式図である。
【図１８】実施例２に係る断層画像の撮影原理を説明する模式図である。
【図１９】本発明の１変形例の構成を説明する模式図である。
【図２０】従来構成の骨梁解析装置の構成を説明する模式図である。
【実施例１】
【００２７】
次に、本発明に係る骨梁解析装置の実施例について図面を参照しながら説明する。なお、実施例におけるＸ線は、本発明の構成の放射線に相当する。なお、ＦＰＤは、フラットパネル型Ｘ線検出器（フラット・パネル・ディテクタ）の略である。そして、断層画像Ｄは、本発明の二次元画像に相当する。
【００２８】
図１は、実施例１に係る骨梁解析装置の構成を説明する機能ブロック図である。図１に示すように、実施例１に係る骨梁解析装置１は、Ｘ線断層撮影の対象である被検体Ｍを載置する天板２と、天板２の上部（天板２の１面側）に設けられた被検体Ｍに対してコーン状のＸ線ビームを照射するＸ線管３と、天板２の下部（天板の他面側）に設けられ、被検体Ｍを透過したＸ線を検出するＦＰＤ４と、コーン状のＸ線ビームの中心軸とＦＰＤ４の中心点とが常に一致する状態でＸ線管３とＦＰＤ４との各々を被検体Ｍの関心部位を挟んで互いに反対方向に同期移動させる同期移動機構７と、これを制御する同期移動制御部８と、ＦＰＤ４のＸ線を検出するＸ線検出面を覆うように設けられた散乱Ｘ線を吸収するＸ線グリッド５とを備えている。この様に、天板２は、Ｘ線管３とＦＰＤ４とに挟まれる位置に配置されている。Ｘ線管３は、本発明の放射線源に相当し、ＦＰＤ４は、本発明の検出手段に相当する。
【００２９】
同期移動機構７は、Ｘ線管３を被検体Ｍに対して体軸方向Ａに移動させるＸ線管移動機構７ａと、ＦＰＤ４を被検体Ｍに対して体軸方向Ａに移動させるＦＰＤ移動機構７ｂとを備えている。また、同期移動制御部８は、Ｘ線管移動機構７ａを制御するＸ線管移動制御部８ａとＦＰＤ移動機構７ｂを制御するＦＰＤ移動制御部８ｂとを備えている。Ｘ線管移動機構７ａは、本発明の放射線源移動手段に相当し、ＦＰＤ移動機構７ｂは、本発明の検出器移動手段に相当する。また、Ｘ線管移動制御部８ａは、本発明の放射線源移動制御手段に相当し、ＦＰＤ移動制御部８ｂは、本発明の検出器移動制御手段に相当する。
【００３０】
Ｘ線管３は、Ｘ線管制御部６の制御にしたがってコーン状でパルス状のＸ線ビームを被検体Ｍに対して繰り返し照射する構成となっている。このＸ線管３には、Ｘ線ビームを角錐となっているコーン状にコリメートするコリメータが付属している。そして、このＸ線管３と、ＦＰＤ４はＸ線透過画像を撮像する撮像系３，４を生成している。
【００３１】
同期移動機構７は、Ｘ線管３とＦＰＤ４とを同期させて移動させる構成となっている。この同期移動機構７は、同期移動制御部８の制御にしたがって被検体Ｍの体軸方向Ａに平行な直線軌道（天板２の長手方向）に沿ってＸ線管３を直進移動させる。このＸ線管３とＦＰＤ４との移動方向は、天板２の長手方向に一致している。しかも、検査中、Ｘ線管３の照射するコーン状のＸ線ビームは、常に被検体Ｍの関心部位に向かって照射されるようになっており、このＸ線照射角度は、Ｘ線管３の角度を変更することによって、たとえば初期角度−２０°から最終角度２０°まで変更される。この様なＸ線照射角度の変更は、Ｘ線管傾斜機構９が行う。Ｘ線管傾斜制御部１０は、Ｘ線管傾斜機構９を制御する目的で設けられている。
【００３２】
そして、さらに実施例１に係る骨梁解析装置１は、各制御部６，８，１０を統括的に制御する主制御部２５と、断層画像Ｄを表示する表示部２７とを備えている。この主制御部２５は、ＣＰＵによって構成され、各種のプログラムを実行することにより各制御部６，８，１０および後述の各部１１，１２，１４，１５，１６，１７を実現している。記憶部２３は、各部の制御様式や、後述の評価値推定部１７が参照する推定式などの骨梁解析に関するデータの一切を記憶する。操作卓２６は、術者の骨梁解析装置１に対する各操作を入力させるものである。評価値推定部１７は、本発明の評価値推定手段に相当する。
【００３３】
また、同期移動機構７は、上述のＸ線管３の直進移動に同期して、天板２の下部に設けられたＦＰＤ４を被検体Ｍの体軸方向Ａ（天板２の長手方向）に直進移動させる。そして、その移動方向は、Ｘ線管３の移動方向と反対方向となっている。つまり、Ｘ線管３が移動することによってＸ線管３の焦点の位置と照射方向が変化するコーン状のＸ線ビームは、常にＦＰＤ４のＸ線検出面の全面で受光される構成となっている。このように、一度の検査において、ＦＰＤ４は、Ｘ線管３と互いに反対方向に同期して移動しながら、たとえば７４枚の透視画像Ｐ０を取得するようになっている。具体的には、撮像系３，４は、実線の位置を初期位置として、破線で示した位置を介して、図１に示した一点鎖線で示す位置まで対向移動する。すなわち、Ｘ線管３とＦＰＤ４の位置を変化させながら複数のＸ線透過画像が撮影されることになる。ところで、コーン状のＸ線ビームは常にＦＰＤ４のＸ線検出面の全面で受光されるので、撮影中コーン状のＸ線ビームの中心軸は、常にＦＰＤ４の中心点と一致している。また、撮影中、ＦＰＤ４の中心は、直進移動するが、この移動はＸ線管３の移動の反対方向となっている。つまり、体軸方向ＡにＸ線管３とＦＰＤ４とを同期的、かつ互いに反対方向に移動させる構成となっている。図１における符号Ｓは被検体Ｍの体側方向を表している。
【００３４】
すなわち、同期移動機構７は、Ｘ線管３を天板２の長手方向における一端側に向けて移動させるのに同期してＦＰＤ４を天板２の長手方向における他端側に向けて移動させるような動作をする。
【００３５】
また、ＦＰＤ４の後段には、そこから出力される検出信号を基に透視画像Ｐ０を生成する画像生成部１１が備えられており（図１参照），この画像生成部１１の更に後段には、透視画像Ｐ０を合成して断層画像Ｄを生成する断層画像生成部１２とを備えている。画像生成部１１は、本発明の画像生成手段に相当し、断層画像生成部１２は、本発明の断層画像生成手段に相当する。
【００３６】
続いて、実施例１に係る骨梁解析装置１の断層画像の取得原理について説明する。図２は、実施例１に係るＸ線撮影装置の断層画像の取得方法を説明する図である。例えば、天板２に平行な（鉛直方向に対して水平な）仮想平面（基準裁断面ＭＡ）について説明すると、図２に示すように、基準裁断面ＭＡに位置する点Ｐ，Ｑが、常にＦＰＤ４のＸ線検出面の不動点ｐ，ｑのそれぞれに投影されるように、Ｘ線管３によるコーン状のＸ線ビームＢの照射方向に合わせてＦＰＤ４をＸ線管３の反対方向に同期移動させながら一連の透視画像Ｐ０が画像生成部１１にて生成される。一連の透視画像Ｐ０には、被検体Ｍの投影像が位置を変えながら写り込んでいる。そして、この一連の透視画像Ｐ０を断層画像生成部１２にて再構成すれば、基準裁断面ＭＡに位置する像（たとえば、不動点ｐ，ｑ）が集積され、Ｘ線断層画像としてイメージングされることになる。一方、基準裁断面ＭＡに位置しない点Ｉは、ＦＰＤ４における投影位置を変化させながら一連の被検体画像に点ｉとして写り込んでいる。この様な点ｉは、不動点ｐ，ｑとは異なり、断層画像生成部１２でＸ線透過画像を重ね合わせる段階で像を結ばずにボケる。このように、一連の透視画像Ｐ０の重ね合わせを行うことにより、被検体Ｍの基準裁断面ＭＡに位置する像のみが写り込んだＸ線断層画像が得られる。このように、透視画像Ｐ０を単純に重ね合わせると、基準裁断面ＭＡにおける断層画像Ｄが得られる。
【００３７】
さらに、断層画像生成部１２の設定を変更することにより、基準裁断面ＭＡに水平な任意の裁断面においても、同様な断層画像を得ることができる。撮影中、ＦＰＤ４において上記点ｉの投影位置は移動するが、投影前の点Ｉと基準裁断面ＭＡとの離間距離が大きくなるにしたがって、この移動速度は増加する。これを利用して、取得された一連の被検体画像を所定のピッチで体軸方向Ａにずらしながら再構成を行うようにすれば、基準裁断面ＭＡに平行な裁断面における断層画像Ｄが得られる。このような一連の被検体画像の再構成は、断層画像生成部１２が行う。このように、断層画像生成部１２は、Ｘ線管３およびＦＰＤ４を被検体Ｍに対して移動させながら連写された画像を基に被検体Ｍを載置する天板に平行な断面に係る断層画像Ｄを生成する。
【００３８】
＜骨梁解析部の構成＞
生成された断層画像Ｄは、骨梁解析部１４，１５，１６，１７に送られる。この骨梁解析部１４，１５，１６，１７は、図３に示す背景処理部１４，行列生成部１５，テクスチャ解析指標算出部１６および評価値推定部１７をまとめ機能ブロックの一つとして表現したものとなっている。骨梁解析部１４，１５，１６，１７は、断層画像Ｄに種々の画像処理を施して骨梁解析を行う。行列生成部１５は、本発明の同時生起行列生成手段に相当し、テクスチャ解析指標算出部１６は、本発明の指標算出手段に相当する。
【００３９】
本発明の最も特徴的な構成は、テクスチャ解析指標ｔ１，ｔ２，ｔ３……とを用いて骨梁の状態を表す評価値を推定する構成にある。すなわち、骨梁解析部１４，１５，１６，１７は、断層画像Ｄよりテクスチャ解析指標ｔ１，ｔ２，ｔ３……とを算出して、これらを用いて種々の推定値Ｐを算出する。図３において推定される対象はＢＶ／ＴＶ値となっている。ＢＶ／ＴＶ値とは骨梁が所定の空間を占める割合を表す骨梁を評価するときの評価値の一種である。術者はこのＢＶ／ＴＶ値の推定値Ｐを被検体Ｍの骨梁がどのような状態になっているかの判断材料とすることができる。
【００４０】
テクスチャ解析指標ｔ１，ｔ２，ｔ３……とは、骨梁の状態を数値化したパラメータである。したがって、断層画像Ｄから得られたこれらパラメータで骨梁の状態の評価できるものと予想される。実際はこのような予想に反し、これらパラメータでは骨梁の状態を正確に把握することができない。これらのパラメータは、断層画像Ｄから取得されるものだからである。骨梁の状態を評価するには骨梁が３次元空間内にどのように張り巡らされているかを考慮する必要がある。テクスチャ解析指標ｔ１，ｔ２，ｔ３……とは、断層画像Ｄから算出されるものであるので、骨梁の三次元構造に対して考慮がされていない。したがって、これらパラメータは、骨梁の状態を正確に知るには、力不足なのである。
【００４１】
一方ＢＶ／ＴＶ値は、骨の海綿質の３Ｄデータを基に算出される骨梁の評価値であるが、生体での実測は困難である。このＢＶ／ＴＶ値は、本来ならば高分解能のＣＴ装置により算出されるものであるので、本発明のような断層撮影装置では撮影することは原理上できない。したがって、本発明に係る装置は、実際にＢＶ／ＴＶ値を算出する代わりに断層画像Ｄ由来のテクスチャ解析指標ｔ１，ｔ２，ｔ３……からＢＶ／ＴＶ値の推定を行うようにしている。この推定値Ｐの算出には、テクスチャ解析指標ｔ１，ｔ２，ｔ３……のみならず実際ＣＴ撮影でられるＢＶ／ＴＶ値も用いるようにしているので、本発明に係る装置は、生体での実測が難しいＢＶ／ＴＶ値の推定を正確に行うことができるようになっている。
【００４２】
続いて、骨梁解析部１４，１５，１６，１７を構成する各部の詳細について説明する。
【００４３】
＜背景処理部１４＞
断層画像Ｄは、まず背景処理部１４に送出される。背景処理部１４は、断層画像Ｄにアンシャープネスマスキング処理などの背景処理を施し、背景処理がされた断層画像Ｄａを生成する。このような背景処理により異なる断層画像Ｄ同士で見られる撮影条件の相違による被検体像の写り込みの違いを消去することができる。背景処理部１４が行う背景処理は、アンシャープネスマスキング処理に限られず他の処理としてもよい。
【００４４】
図４は、アンシャープネスマスキング処理を簡単に説明している。アンシャープネスマスキング処理とは、断層画像Ｄから断層画像Ｄを平滑化して得られた平滑化画像を差し引く画像処理である。この様な画像処理により、断層画像Ｄは、よりくっきりした断層画像Ｄａに変換される。
【００４５】
＜行列生成部１５＞
テクスチャ解析を行う際に必要となる行列として同時生起行列（ＧＬＣＭ）がある。この行列は行列生成部１５により生成される。背景処理部１４が生成した断層画像Ｄａは、行列生成部１５に送出され、そこでＧＬＣＭに変換される。図５は、行列生成部１５が断層画像Ｄａに基づいてＧＬＣＭを生成する動作を説明している。図５の左側は、断層画像Ｄａを画素値の２次元配列として表している。説明の簡単のため、断層画像Ｄａを構成する各画素の画素値は、０から９までの１０通りの値をとるものとする。
【００４６】
図５の右側に示すように、断層画像Ｄａより生成されるＧＬＣＭの行数と列数は、いずれも画素の画素値がとりえる画素値の数と一致する。断層画像Ｄａを構成する各画素は、１０通りのうちのいずれかの画素値を有しているのであるから、断層画像Ｄａより生成されるＧＬＣＭは１０行１０列の２次元行列となる。行列生成部１５は、１０×１０行列となっているＧＬＣＭを構成する１００個の要素に数値を代入してＧＬＣＭを完成させる。各要素にどのような数値を入れるかは、断層画像Ｄａの画素値に基づいて判断される。
【００４７】
図５は、ＧＬＣＭの各行のうち０を意味する行、各列のうち１を意味する行に位置する要素ｐ（０，１）の数値を行列生成部１５が決めようとしているところを示している。行列生成部１５は、画素値０と画素値１とが隣り合って配列されている画素のペアが断層画像Ｄａに何組あるかを数えて、そのカウント数をＧＬＣＭの要素ｐ（０，１）とする。図５においては、画素値０と画素値１とが隣り合って配列されている画素のペアは２組あるので、要素ｐ（０，１）の値は、２となる。このＧＬＣＭにおける任意の要素ｐ（ａ，ｂ）は要素ｐ（ｂ，ａ）に等しいので、ＧＬＣＭの要素ｐ（１，０）の値も２となる。
【００４８】
行列生成部１５は、同様な動作をＧＬＣＭの全域に亘って行い、行列が有する要素の全てを断層画像Ｄａに基づいて決定する。こうして行列生成部１５は、断層画像Ｄａに基づいてＧＬＣＭを完成させる。
【００４９】
図６は、行列生成部１５が断層画像Ｄａに基づいてＧＬＣＭを生成する様子を示している。生成されるＧＬＣＭは、断層画像Ｄａの画素が取り得る画素値の数が多くなるほど大きくなる。ＧＬＣＭは、対称性を有する行列であり、図６の点線で示す対角線で２つ折りにすると、重なり合う要素同士の値が同じとなっているような行列である。
【００５０】
このように、行列生成部１５は、断層画像Ｄａの一部に設けられた解析範囲を構成する各画素のうち所定の画素値の組み合わせを有する２つの画素のペアで画素同士が所定の距離だけ離間しているものが解析範囲において何回現れるかを各画素値の組み合わせごとに数えてＧＬＣＭ（同時生起行列）を生成する。解析範囲としては、骨の海綿骨の部分が選択される。行列生成部１５の解析範囲は骨梁総延長算出部１４の解析範囲Ｒに一致してい
る。
【００５１】
＜テクスチャ解析指標算出部１６＞
ＧＬＣＭは、テクスチャ解析指標算出部１６に送出される。テクスチャ解析指標算出部１６は、ＧＬＣＭに種々の演算を実行することでテクスチャ解析指標ｔ１，ｔ２，ｔ３……を算出することが可能である。テクスチャ解析指標算出部１６が算出できるテクスチャ解析指標は、例えば次のようなものがある。式中のｐ（ｉ，ｊ）とは、ＧＬＣＭにおけるｉ行ｊ列目の要素の値、Σｉ，Σｊは、それぞれｉ行、ｊ列についての要素の合計、Ｎｇは、断層画像Ｄａの画素が取り得る画素値の数、μは平均値、μｘ，μｙは、それぞれ行方向、列方向の平均値、σｘ，σｙは、それぞれ行方向、列方向の標準偏差を表している。なお、これらテクスチャ解析指標ＡＳＭ（Ａｎｇｕｌａｒ  Ｓｅｃｏｎｄ  Ｍｏｍｅｎｔ：アングラーセカンドモーメント），ＣＮＴ（Ｃｏｎｔｒａｓｔ：コントラスト），ＣＯＲ（Ｃｏｒｒｅｌａｔｉｏｎ：コリレーション），ＶＡＲ（Ｖａｒｉａｎｃｅ：バリアンス），ＩＤＭ（Ｉｎｖｅｒｓｅ  Ｄｉｆｆｅｒｅｎｃｅ  Ｍｏｍｅｎｔ，インバースディファレンシャルモーメント），ＥＮＴ（Ｅｎｔｒｏｐｙ，エントロピー）の各々は、１９７３年にＨａｒｌｉｃｋらが上述の非特許文献２で提唱した１４種類のパラメータのうちの一部である。また、ＤＩＳは非類似度またはディシミラレィティと呼ばれるテクスチャ解析指標で、ＨＯＭは、一様性またはホモジェネイティと呼ばれるテクスチャ解析指標である。
【００５２】
【数１】
2015107142.tif 000003
【００５３】
テクスチャ解析指標算出部１６は、ＧＬＣＭに上述の種々の演算を行ってテクスチャ解析指標ｔ１，ｔ２，ｔ３……を算出する。テクスチャ解析指標算出部１６が算出するテクスチャ指標の種類と数は、適宜変更することができる。テクスチャ解析指標の数は３つ以下でもよい。以上のように、テクスチャ解析指標算出部１６は、ＧＬＣＭ（同時生起行列）に基づいてテクスチャ解析を行いテクスチャ解析指標ｔ１，ｔ２，ｔ３……を算出する。
【００５４】
＜評価値推定部１７＞
テクスチャ解析指標ｔ１，ｔ２，ｔ３……は、評価値推定部１７に送出される。評価値推定部１７は、テクスチャ解析指標ｔ１，ｔ２，ｔ３……を入力すると生体での実測が難しい評価値の推定値Ｐが出力される推定式に基づいて推定値Ｐを算出する。評価値推定部１７が推定する評価値としては、ＢＶ／ＴＶ値、ＴＢＰｆ値、ＴＳＬ値などがある。これらは、本来は骨を３次元解析しなければ得られない骨梁の評価値である。評価値推定部１７は、断層画像Ｄａを解析して得られた各テクスチャ解析指標ｔ１，ｔ２，ｔ３……を用いて、この断層画像Ｄａに写り込んでいる骨を３次元解析したとすれば得られるであろう評価値を推定するのである。３次元解析は、正確ではあるが生体に実施するのは相当難しい。本発明は、この様な事情に鑑みて断層画像Ｄａから３次元解析の結果を推定するようにしている。
【００５５】
評価値の各々の意味について説明する。ＴｂＴｈ値とは骨の海綿組織をＣＴ撮影して３次元解析をすることで得られるものである。ＴｂＴｈ値の意味は、骨梁の幅である。この値が高いほど、海綿組織がより骨の成分で満たされていることを示し、強い骨であることを表す。また、ＢＶ／ＴＶ値は、関心領域中における骨梁の部分と中空の部分との割合である。この値が高いほど、海綿組織がより骨の成分で満たされていることを示し、強い骨であることを表す。このＢＶ／ＴＶ値もＣＴ撮影による３次元解析で得られる値である。
【００５６】
ＴＢＰｆ値は、体積／表面積で算出される値で、ＢＶ／ＴＶ値と同じく海綿組織の詰まり具合を数値化するものである。このＴＢＰｆ値もＣＴ撮影による３次元解析で得られる値である。図７は、ＴＢＰｆ値の意味を説明している。図７の左側は、ＣＴ撮影で得られた３次元画像の一部（関心領域）であるとする。この立方体の形状をしている関心領域は周縁部が骨成分で構成されるものの、中心には球形の中空が包含されているものとする。この関心領域にグロー処理を施して、球形の中空を小さくしたとする。グロー処理とは、中空の周辺を骨成分に置換することにより、骨成分と中空との界面に骨成分の層を積層するような画像処理である。グロー処理後の関心領域は、図７の右側に示されている。
【００５７】
この様な画像処理を施すと、骨成分と中空との海面の表面積と関心領域に占める骨成分の体積とが変化するはずである。それぞれの変化量をΔＳ，ΔＶとすると、ＴＢＰｆ値は、ΔＳ／ΔＶで表される値である。骨梁解析においては、一概には言えないが、値が低いほど骨が強いものであることを表す場合が多い。
【００５８】
ＴＳＬ値は、海綿骨をＣＴ撮影して得られた３次元画像を基にした骨梁総延長である。なお、ＴＳＬ値は大きいほど骨が強いものであることを示す。
【００５９】
評価値推定部１７の動作についてより詳細に説明する。評価値推定部１７には、記憶部２３に記憶されている推定式を読み出して評価値を推定する。このときの数式は次のようなものである。
Ｐ＝ｋ１・ｔ１＋Ｃ  ……（１）
ここで、Ｐは推定値、ｔ１はテクスチャ解析指標、ｋ１は係数、Ｃは定数を表している。このように、本発明によれば、テクスチャ解析指標ｔ１，ｔ２，ｔ３……と各評価値との関連性は（１）式のような推定式によって表現されている。
【００６０】
上述の１式は、推定値Ｐを１つのテクスチャ解析指標から求めていたが、本発明はこれに限らず、推定値Ｐの算出に用いるテクスチャ解析指標の個数を適宜増減することができる。以降、１つのテクスチャ解析指標から推定値Ｐを求めるものとして説明する。
【００６１】
＜推定式の決定＞
評価値推定部１７が動作に用いる推定式をどのように決定するのかについて説明する。推定式を完成させるには、数あるテクスチャ解析指標のうちどれを用いるのかと、各係数と定数の決定とを行わなければならない。このような推定式は、被検体Ｍの骨梁解析に先立って標本骨を用いて取得される。標本骨とは、生体より摘出した海綿骨のサンプルである。標本骨は、生体から取り出されたものなので、撮影視野が狭いＣＴ装置や、Ｘ線量が非常に強いＣＴ装置などにセットして骨梁がはっきり写り込んだ三次元画像が取得できる。
【００６２】
標本骨を用いて推定式を決定するには、同じ標本骨にＣＴ撮影と断層画像撮影との両方を行って、それぞれの解析結果を得る必要がある。以下の表は、標本骨についてμＣＴ撮影と本発明に係る装置を用いて断層撮影を行った時の解析結果を示す。μＣＴ撮影では、ＢＶ／ＴＶ値、ＴＢＰｆ値、ＴｂＴｈ値、ＴＳＬ値が算出される。断層撮影では、各テクスチャ解析指標、ＡＳＭ値、ＨＯＭ値、ＥＮＴ値、ＣＯＲ値、ＤＩＳ値、ＣＮＴ値、ＶＡＲ値の各々が算出される。表１では、μＣＴ撮影で得られる値と断層撮影で得られる値との相関性（Ｒ２値）が表されている。値が１に近づくほど相関が大きい。なお、表中のＴＢＰｆの単位はｍｍ−１である。
【００６３】
【表１】
2015107142.tif 000004
【００６４】
このような値同士の相関性をどのように得ているのかについて説明する。例えば、ＢＶ／ＴＶ値とＡＳＭ値（断層画像解析で得られた骨梁総延長）の相関性を得るには、次のような操作が行われる。まず、標本骨の関心領域を変えながらμＣＴを行う。こうして、例えば標本骨の５カ所について５通りのＢＶ／ＴＶ値が得られる。続いて、μＣＴで設定した関心領域の各々について断層撮影を行う。こうして断層画像上においてμＣＴの関心領域に対応する部分のＡＳＭ値が５通り得られる。これで、５通りのＢＶ／ＴＶ値のそれぞれには、対応するＡＳＭ値が用意されたことになる。このＢＶ／ＴＶ値とＡＳＭ値とのペアのそれぞれを、２次元座標上でプロットしてみる。図８は、その結果である。図８では、縦軸が３次元解析由来のＢＶ／ＴＶ値、横軸が断層像由来のＡＳＭ値となっている。
【００６５】
図８におけるプロットを基に、回帰分析行うことができる。図９は、回帰分析により近似式（一次近似式）を求めた結果である。このとき相関性を示す相関係数の２乗値、Ｒ２値が算出できる。図９におけるＲ２値は、０．７５である。
【００６６】
この様な要領で他の値同士についてＲ２値を算出することができる。表１は、算出されたＲ２値をテーブル化したものである。この表１を見ると気が付くことがある。断層撮影で得られた各値は、μＣＴで得られるＴｂＴｈ値との間で有意な相関が見られないということである。例えば、ＤＩＳ値とＴｂＴｈ値についてのＲ２値が０．９９であったとすると、図１０に示すように、ＴｂＴｈ値とＤＩＳ値とは非常に高い相関を有している。このことは、ＤＩＳ値さえ求めれば、ＴｂＴｈ値を求めなくても、ＴｂＴｈ値を高い信頼性のうちに推定できるということを意味している。例えば、ある被検体Ｍに断層撮影を行い、この結果、図１０のａというＤＩＳ値が得られたとすると、このときのＴｂＴｈ値は、実際に測定してみなくても近似曲線からｂという値になっていることが容易に予想できる。
【００６７】
ところが、ＤＩＳ値とＴｂＴｈ値についての実際のＲ２値は、０．００であり、全く相関がない。従って、いくらＤＩＳ値を正確に多数取得しても、ＴｂＴｈ値を予想することはできないわけである。この例のように、表１が示す結果は、断層撮影から取得される各値は、単独ではＴｂＴｈ値等の推定に用いることができないことを表している。
【００６８】
この様な原理に基づいて、各評価値を推定するのには、どのような回帰式が有効であるかを決定することができる。すなわち、各評価値を推定するのには、各評価値についてＲ２値が最も高いテクスチャ解析指標を用いれば、信頼性がより高くなる。例えば、ＢＶ／ＴＶ値を推定するには、ＣＮＴ値に係る回帰分析の結果得られた回帰式がふさわしく、ＴＢＰｆ値を推定するには、ＶＡＲ値に係る回帰分析の結果得られた回帰式がふさわしいことがわかる。なお、Ｒ２値が０．５以上であれば、回帰分析で得られた回帰式は、実測が難しい各評価値を予想する上で有効であるものとされる。
【００６９】
各評価値と各テクスチャ解析指標との相関性は、表１に係る方法とは異なる方法によっても知ることもできる。すなわち、シミュレーションにより両者の相関を知ることができる。図１１，１２，１３は、行ったシミュレーションの例を示している。図１１は、空間上の骨梁を次第に太くなっていくような一連の画像を生成して、ＴｂＴｈ値とテクスチャ解析指標との相関を知ろうとしている。各画像についてのＴｂＴｈ値の算出は容易に行え、各画像についてのテクスチャ解析指標の算出も容易に行える。この様にして得られたデータ群についてＴｂＴｈ値を目的変数とし、テクスチャ解析指標の各々を独立変数とした単回帰分析を行い、得られたＲ２値を比較することで、ＴｂＴｈ値と各テクスチャ解析指標との間に相関が見られるかを知ることができる。
【００７０】
図１２は、空間上の骨の中空が次第に小さくなっていくような一連の画像を生成して、ＴＢＰｆ値とテクスチャ解析指標との相関を知ろうとしている。各画像についてのＴＢＰｆ値の算出は容易に行え、各画像についてのテクスチャ解析指標の算出も容易に行える。この様にして得られたデータ群についてＴＢＰｆ値を目的変数とし、テクスチャ解析指標の各々を独立変数とした単回帰分析を行い、得られたＲ２値を比較することで、ＴＢＰｆ値と各テクスチャ解析指標との間に相関が見られるかを知ることができる。
【００７１】
図１３は、空間上の骨を次第に長くしていくような一連の画像を生成して、ＴＳＬ値とテクスチャ解析指標との相関を知ろうとしている。各画像についてのＴＳＬ値の算出は容易に行え、各画像についてのテクスチャ解析指標の算出も容易に行える。この様にして得られたデータ群についてＴＳＬ値を目的変数とし、テクスチャ解析指標の各々を独立変数とした単回帰分析を行い、得られたＲ２値を比較することで、ＴＳＬ値と各テクスチャ解析指標との間に相関が見られるかを知ることができる。
【００７２】
次に、示す表は、この様なシミュレーションの結果を表している。表中の○は、互いに相関があると判定されたパラメータの組み合わせを表し、×は、互いに相関がないと判定されたパラメータの組み合わせを表し、△は、互いの相関が弱いと判定されたパラメータの組み合わせを表している。
【００７３】
【表２】
2015107142.tif 000005
【００７４】
本発明によれば、シミュレーションにより複数の回帰式を取得し、表２の結果に従って、複数の回帰式の中から評価値の推定に用いる推定式を決定できる。
【００７５】
図１４は、本発明の概要をまとめたものとなっている。本発明に係る骨梁解析の下準備として、まず標本骨に対してμＣＴ撮影および断層画像撮影が行われ、得られたμＣＴ像および断層画像の画像解析が行われる。画像解析結果のうち、断層画像由来のパラメータを独立変数とし、μＣＴ像由来のパラメータを目的変数として回帰分析を行い、推定式を算出する。この推定式は、独立変数を変えて算出された数ある推定式のうちから最も信頼性が高い（Ｒ２値が高い）ものとなっている。
【００７６】
下準備により用意された推定式は、記憶部２３に記憶される。被検体Ｍの骨梁解析を行うときは、まず、断層画像Ｄが撮影され、画像処理により断層画像Ｄａが生成される。この断層画像Ｄａは骨梁解析およびテクスチャ解析がなされ、解析目的の骨の部位についてびテクスチャ解析指標ｔ１が算出される。評価値推定部１７はＢＶ／ＴＶ値の推定値Ｐは記憶部２３に記憶されている推定式に断層画像Ｄａより実際に算出されたテクスチャ解析指標ｔ１を代入することで求められる。
【００７７】
以上のように、本発明の骨梁解析装置１によれば、断層画像Ｄａを基にした骨梁解析でありながらも、骨梁の状態をより正確に把握できる。すなわち、本発明に係る骨梁解析装置１は、断層画像Ｄａに画像解析を施すことで骨梁の大まかな状態を表す数値を取得する。このときの数値とは、本発明においてテクスチャ解析指標ｔ１，ｔ２，ｔ３……である。これら数値は、骨梁の状態により変動するパラメータではあるものの、骨梁の三次元構造の状態を盛り込んでいないから骨梁の状態を正確に表したものとはならない。そこで、本発明では、この断層画像Ｄａで得られた数値に基づいて、断層画像Ｄａの被写体が仮に三次元解析されたとしたときに得られる評価値を推定するようにしている。この推定式は標本骨について実際に三次元解析および断層画像解析の両方をして、目的変数を三次元解析に係る評価値とし、独立変数を断層画像解析に係るテクスチャ解析指標ｔ１，ｔ２，ｔ３……としたときの回帰分析によって得られたものである。この様な構成とすることにより、断層画像解析をするだけで三次元解析に裏付けられた推定値Ｐが算出可能となり、より骨梁の状態を正確に把握できる骨梁解析装置１が提供できる。断層画像Ｄａから容易に取得できるテクスチャ解析指標ｔ１，ｔ２，ｔ３……により、従来では行えなかった正確な骨梁の状態を把握できる骨梁解析装置が提供できる。
【実施例２】
【００７８】
続いて、実施例２に係る骨梁解析装置について説明する。実施例２の構成は、図１５に示すように、Ｘ線管３とＦＰＤ４とが互いの位置関係を保った状態で被検体Ｍの体軸方向Ａに移動されながら断層画像を撮影することができる構成である。すなわち、同期移動機構７は、Ｘ線管３を天板２の長手方向における一端側に向けて移動させるのに同期してＦＰＤ４を天板２の長手方向における一端側に向けて移動させるような動作をする。
【００７９】
実施例２に係るＸ線撮影装置の構成は図１における機能ブロック図と同様である。図１に関して実施例２の構成が実施例１と異なる点は、ＦＰＤ４がＸ線管３に追従して移動すること（図１５参照），Ｘ線管３が傾斜しないことである。したがって、実施例２においては図１におけるＸ線管傾斜機構９，Ｘ線管傾斜制御部１０は必ずしも必要とされない。
【００８０】
実施例２に係る断層画像の撮影の原理について説明する。まず、図１５に示すように撮像系３，４が相対位置を保った状態で被検体Ｍに対して移動しながら間歇的にＸ線を照射する。つまり一度の照射が終了する毎にＸ線管３は被検体Ｍの体軸方向Ａに移動し、再びＸ線の照射を行う。こうして複数枚の透過画像が取得され、透過画像の加工画像（後述の長尺透過画像）がフィルタバックプロジェクション法により断層画像に再構成される。完成した断層画像は、被検体Ｍをある裁断面で裁断したときの断層像が写りこんだ画像となっている。
【００８１】
断層画像を生成するには、異なる方向から被検体Ｍを透視したときの画像が必要となる。実施例２に係る骨梁解析装置は、得られた透過画像を分割してつなぎ合わせてこの画像を生成するようにしている。この動作について説明する。図１６は、Ｘ線管３のＸ線を照射する焦点がｄ１の位置にあるときのＦＰＤ４の位置を表している。この撮影において、被検体Ｍの体軸方向ＡにおけるＦＰＤ４の１／５の幅だけＸ線管３およびＦＰＤ４が天板２に対してこの方向に移動する度に透過画像の撮影が行われるものとする。
【００８２】
Ｘ線はＸ線管３から放射状に広がってＦＰＤ４に到達するので、生成された透過画像を被検体Ｍの体軸方向Ａに５分割すると、ＦＰＤ４に対するＸ線の入射角度は、矢印に示すように、その分割区の間で互いに異なっている。そのうちのあるの１つの方向ｋに注目する。この方向ｋに進んできたＸ線は、被検体Ｍの斜線の部分を通過してＦＰＤ４に写り込んでいるので、方向ｋのＸ線が入射したＦＰＤ４の分割区には、被検体Ｍの斜線部が写り込んでいる。透過画像において、この分割区に相当する部分を断片Ｒ１とする。
【００８３】
図１７は、Ｘ線管３のＸ線を照射する焦点がｄ１からＦＰＤ４の１／５の幅だけ移動したｄ２の位置にあるときのＦＰＤ４の位置を表している。Ｘ線管３とＦＰＤ４の位置関係は変化しないので、このときの撮影においてもＦＰＤ４には、方向ｋに進んできたＸ線が写り込んでいる分割区があるはずであり、方向ｋのＸ線が入射したＦＰＤ４の分割区には、被検体Ｍの斜線部が写り込んでいる。透過画像において、この分割区に相当する部分を断片Ｒ２とする。
【００８４】
断片Ｒ１と断片Ｒ２とを比較すると、撮像系３，４に対する被検体Ｍの位置が異なるので、両断片Ｒ１，Ｒ２に写り込んでいる被検体Ｍの部分は互いに異なっている。Ｘ線管３をＦＰＤ４の１／５の幅だけずらすことにより、焦点ｄ１〜ｄ９において９回の撮影を行ったとして、そのときの方向ｋのＸ線が入射したＦＰＤ４の分割区における透過画像の各断片Ｒ１〜Ｒ９には、それぞれ異なる被検体Ｍの位置が写り込んでいる。そこで、図１８に示すように透過画像の各断片Ｒ１〜Ｒ９をこの順に被検体Ｍの体軸方向Ａにつなぎ合わせれば、ある方向ｋで被検体Ｍの全身にＸ線を照射したときに撮影される画像を得ることができる。この画像を長尺透過画像と呼ぶことにする。
【００８５】
実施例２に係る骨梁解析装置は、断層画像生成部１２において方向ｋ以外の方向についても長尺透過画像を生成する。そして、断層画像生成部１２は、被検体Ｍを投影した方向が異なる複数の長尺透過画像を基に被検体Ｍを所定の裁断位置で裁断したときの断層画像Ｄを生成するのである。
【００８６】
実施例２における断層画像Ｄについて行われる解析は、実施例１と同様であり、最終的に骨梁の評価値の推定値Ｐが算出される。
【００８７】
以上のように、実施例２の構成によれば、スロット撮影を仮想的に行うことにより取得された長尺画像を撮影しこれらから断層画像Ｄを撮影する構成となっている。この様な撮影を行うようにすれば、広範囲に亘って撮影された断層画像を取得できる放射線撮影装置を提供できる。
【００８８】
本発明は、上述の構成に限られず、下記のように変形実施することが可能である。
【００８９】
（１）上述の実施例では、表１に示すような８種類のテクスチャ解析指標が用いられていたが、本発明はこの構成に限られない。本発明に同時生起行列より導き出せる他のテクスチャ解析指標を用いることもできる。また、本発明においては、非特許文献２で提唱されている以外のテクスチャ解析指標を用いることもできる。
【００９０】
（２）上述の実施例における評価値推定部１７が推定値の算出に用いるテクスチャ解析指標と評価値との関連性は、推定式により表されていたが、本発明はこの構成に限られない。本発明をテクスチャ解析指標と評価値の推定値とが関連したテーブルとして関連性を表す構成としてもよい。
【００９１】
（３）上述の実施例では、骨の部位について何ら言及がなかったが、骨の部位ごとに異なる推定式を用いて評価値の推定を行うようにしてもよい。本変形例では、例えば、大腿骨の頚部と骨幹部のそれぞれについて回帰分析を行い、各部位に対応した推定式を算出するようにしてもよい。この頚部の推定式は、頚部の解析結果で得られた複数の推定式のうち最もＲ２値が高かったものであり、骨幹部の推定式は、骨幹部の解析結果で得られた複数の推定式のうち最もＲ２値が高かったものとなっている。
【００９２】
（４）上述の実施例の重回帰分析は、１次近似法によって行っていたが、本発明はこの構成に限られない。重回帰分析を２次近似法によって行うようにしてもよい。また、重回帰分析をより高次の近似法により行うようにしてもよい。
【００９３】
（５）上述の実施例では、独立変数と目的変数とが１対１となっている単回帰分析により推定式を算出していたが、本発明はこれに代えて、複数の独立変数を用いた重回帰分析により推定式を算出するようにしてもよい。この場合、独立変数としては、各テクスチャ解析指標のうちから選んで使用することができる。この場合においても、様々な独立変数の組み合わせで回帰式を算出してみてＲ２値が最も高いものを探索することで推定式を決定することができる。
【００９４】
（６）上述の実施例では、骨梁の評価値は、ＢＶ／ＴＶ値、ＴＢＰｆ値、ＴＳＬ値、ＴｂＴｈ値であるとしていたが、本発明はこの構成に限られない。本発明は、骨梁数や異方性などの他の骨梁の評価値についても適用できる。
【００９５】
（７）上述の実施例の行列生成部１５は、断層画像Ｄａにおいて互いに隣り合った画素のペアの個数を数えるように動作していたが、本発明はこの構成に限られない。すなわち、図１９に示すように、所定の距離だけ離間した画素のペアの個数を数えて同時生起行列を生成するようにしてもよい。図１９の例では、両方の画素値が４になっている１画素の幅だけ離間した画素のペアを行列生成部１５がカウントしている様子を表している。
【００９６】
（８）上述した実施例では、被検体Ｍに既に骨梁解析が行われていることが前提となっていたが本発明はこの構成に限られない。すなわち、本発明の骨梁解析装置は、以前に骨梁解析が行われていない被検体Ｍに対しても推定値Ｐを算出するような動作をしてもよい。このときの断層画像Ｄは、過去または、別の装置により撮影されたものである。この様な構成とすることで、術者が所望する裁断面における断層像を用いて骨梁解析ができるようになり、解析の自由度が向上した装置が提供できる。
【００９７】
（９）上述した実施例は、医用の装置であったが、本発明は、工業用や、原子力用の装置に適用することもできる。
【００９８】
（１０）上述した実施例のいうＸ線は、本発明における放射線の一例である。したがって、本発明は、Ｘ線以外の放射線にも適応できる。
【符号の説明】
【００９９】
３     Ｘ線管（放射線源）
４     ＦＰＤ（検出手段）
７ａ   Ｘ線管移動機構（放射線源移動手段）
７ｂ   ＦＰＤ移動機構（検出器移動手段）
８ａ   Ｘ線管移動制御部（放射線源移動制御手段）
８ｂ   ＦＰＤ移動制御部（検出器移動制御手段）
１１   画像生成部（画像生成手段）
１２   断層画像生成部（断層画像生成手段）
１５   行列生成部（同時生起行列生成手段）
１６   テクスチャ解析指標算出部（指標算出手段）
１７   評価値推定部（評価値推定手段）
【図１】
2015107142.tif 000006
【図２】
2015107142.tif 000007
【図３】
2015107142.tif 000008
【図４】
2015107142.tif 000009
【図５】
2015107142.tif 000010
【図６】
2015107142.tif 000011
【図７】
2015107142.tif 000012
【図８】
2015107142.tif 000013
【図９】
2015107142.tif 000014
【図１０】
2015107142.tif 000015
【図１１】
2015107142.tif 000016
【図１２】
2015107142.tif 000017
【図１３】
2015107142.tif 000018
【図１４】
2015107142.tif 000019
【図１５】
2015107142.tif 000020
【図１６】
2015107142.tif 000021
【図１７】
2015107142.tif 000022
【図１８】
2015107142.tif 000023
【図１９】
2015107142.tif 000024
【図２０】
2015107142.tif 000025