JP5599831B2 - 骨塩定量分析における分析対象部位の特定方法および、この方法を実施するための画像処理装置並びに記録媒体 - Google Patents

Description

本発明は骨塩定量分析方法、特に分析対象の骨部を撮影した放射線画像を利用して骨塩量定量分析を行う方法において、分析対象部位を特定する方法に関するものである。

また本発明は、そのような方法を実施するための画像処理装置並びに、そのような方法をコンピュータに実行させるプログラムを記録したコンピュータ読取り可能な記録媒体に関するものである。

従来、骨粗鬆症の診断等のために、分析対象の骨部の放射線画像を利用して、その骨部の骨塩定量を求める分析方法が知られている。そのような骨塩定量分析方法の中で、比較的簡便に実施可能な方法の一つとして、ＭＤ（Microdensitometry）法と呼ばれる方法が知られている。このＭＤ法は基本的に、放射線管球から発生させた放射線を、分析対象の骨部と、放射線透過特性が互いに異なる複数の部分を有する標準物質とに同時に照射し、骨部および標準物質を透過した放射線をＸ線フィルム等の放射線検出体で検出して該骨部および標準物質を示す放射線画像を得、この放射線画像において、分析対象の骨部と同じ濃度を示す標準物質の部分の放射線透過特性に基づいて該骨部の骨塩定量を求めるものである。

なお上記の標準物質としては一般に、厚さが連続的に変化するアルミスロープが用いられ、その場合は上記放射線透過特性に対応するアルミスロープの厚さを、骨塩定量を示す指標として定義することが多い。

また上述のＭＤ法の中でも、放射線検出体として特に、放射線画像を示すデジタル画像信号を得ることができるものを用い、そのデジタル画像信号を処理して骨塩定量を求めるようにしたＤＩＰ（Digital Image Processing）法が広く知られている（例えば特許文献１、２および３参照）。このＤＩＰ法による骨塩定量分析は、操作が簡単で短時間に実行できることから近時広く普及しつつある。

上記ＤＩＰ法のように放射線画像を用いる骨塩定量分析においては、分析結果に一般性を持たせて、他の分析結果と比較したり、骨粗鬆症等の診断のために蓄積されている診断データと比較することができるように、一定の骨の一定の部位を分析対象として定めることが多い。一般にこの部位は、左手の第２中手骨つまり人差し指の中手骨の所定範囲（例えば該第２中手骨の長さ方向中央位置に有る、骨全長の１／１０の長さに亘る範囲）とされる場合が多い。また、左手に限らず、右手の第２中手骨の所定範囲とすることも考えられている。

特開２００６−３３４０４６号公報 再表２００８−０４４４３９号 特開２０１０−２００８２４号公報

上述したように、左手あるいは右手の第２中手骨の所定範囲を分析対象部位として骨塩定量分析を行うには、まず、放射線画像においてその分析対象部位を特定する必要がある。そのために従来は、分析装置使用者が例えば表示手段に表示されている放射線画像上で分析対象部位の端部等をマウスの操作等によって指定し、その指定入力に基づいて分析対象部位を特定することが多くなされている。

しかし、このように装置使用者のマニュアル操作によって分析対象部位を特定する場合は、その操作が、装置使用者が気付かないまま多少不正確になされたり、ばらつきを持って行われることもある。そうであると、本来分析対象部位とすべき部位とは多少異なる部位について骨塩定量分析がなされてしまうこともあり、そのときに得られた分析結果を、蓄積されている診断データと比較して診断するような場合は、誤った診断が下される可能性もある。

また、分析対象部位を特定するマニュアル操作が不安定であると、同じ放射線画像を使用していても骨塩定量分析結果にばらつきが生じるので、分析の再現性、信頼性が低いものになってしまう。

本発明は上記の事情に鑑みてなされたものであり、骨塩定量分析に用いられる放射線画像から、分析対象となる骨（左手あるいは右手の第２中手骨等）の中の所定領域を、安定して再現性良く特定できる方法を提供することを目的とする。

また本発明は、そのような方法を実施するための画像処理装置並びに記録媒体を提供することを目的とするものである。

本発明による骨塩定量分析における分析対象部位の特定方法は、前述したＤＩＰ法のように、骨部を撮影した放射線画像から該骨部の骨塩定量を求める骨塩定量分析方法において、放射線画像から、画像処理によって分析対象となる骨の中の所定領域を分析対象部位として特定する方法であって、
手が撮影された放射線画像において分析対象となる骨の部分を特定し、
特定された分析対象となる骨の中の特徴点を特定し、
特定されたこの特徴点に基づいて前記所定領域を特定することを特徴とするものである。

なお本発明による分析対象部位の特定方法においては、分析対象となる骨の部分を特定するために、放射線画像において一つの手の第１中手骨、第２中手骨および分析対象となる骨が存在する領域をテンプレートにし、該テンプレートと上記手が撮影された放射線画像とに基づいて分析対象となる骨の部分を特定することが望ましい。

ここで、上記の「第１中手骨、第２中手骨および分析対象となる骨が存在する領域」とは、第１中手骨と第２中手骨の一方が分析対象となる骨である場合は、その骨と他方の骨が存在する領域となる。つまり、例えば第２中手骨が分析対象となる骨である場合は、該第２中手骨と第１中手骨が存在する領域が「第１中手骨、第２中手骨および分析対象となる骨が存在する領域」となる。それに対して、分析対象となる骨が第１中手骨と第２中手骨のいずれでもない場合は、当然、それら３つの骨を含む領域が「第１中手骨、第２中手骨および分析対象となる骨が存在する領域」となる。例えば、分析対象となる骨が第３基節骨（中指の指骨）であれば、第１中手骨、第２中手骨および第３基節骨が存在する領域が「第１中手骨、第２中手骨および分析対象となる骨が存在する領域」となる。これは、以下においても同様である。

また、本発明による分析対象部位の特定方法においては、分析対象となる骨の中の特徴点を特定するために、放射線画像において一つの手の第１中手骨、第２中手骨および分析対象となる骨が存在する領域をテンプレートにし、該テンプレートと上記手が撮影された放射線画像とに基づいて前記特徴点を特定することが望ましい。

そして、そのように一つの手の第１中手骨、第２中手骨および分析対象となる骨が存在する領域をテンプレートとし、上記手が撮影された放射線画像のテンプレートと対応する部分を探索する際には、第１および第２の中手骨どうしが所定角度（例えば２０°）以上の角度をなしていることに基づいて、分析対象となる骨を認識することが望ましい。

また、本発明による分析対象部位の特定方法においては、分析対象となる骨の両端を特定し、これらの両端に基づいて前記特徴点を特定することが望ましい。そして、そのようにして特徴点を特定する場合は、分析対象となる骨の両端を、該分析対象となる骨の関節面を認識することによって特定することが望ましい。

また、本発明による分析対象部位の特定方法においては、手が撮影された放射線画像において、手と共に撮影された骨塩定量分析用の標準物質を認識し、認識されたこの標準物質の位置および／または形状に基づいて、撮影された手が左右どちらの手であるか特定することが望ましい。

他方、本発明による分析対象部位の特定方法においては、前記特徴点として、分析対象となる骨の長さ方向の中央点を特定することが望ましい。

また、本発明による分析対象部位の特定方法においては、前記所定領域として、分析対象となる骨の長さ方向の中央部分に有る、該分析対象となる骨の長さの１／１０に相当する領域を特定することが望ましい。

さらに本発明による分析対象部位の特定方法においては、上に挙げたような特徴点や所定領域を、分析対象となる骨の放射線画像と共に表示手段に表示させることが好ましい。

本発明による骨塩定量分析における分析対象部位の特定方法は、上述した通り画像処理により、手が撮影された放射線画像において分析対象となる骨の部分を特定し、特定された分析対象となる骨の中の特徴点を特定し、特定されたこの特徴点に基づいて所定領域を特定するようにしたので、この方法によれば左手あるいは右手の第２中手骨等の分析対象となる骨の中の所定領域を、安定して再現性良く特定することが可能になる。

なお、人間の一つの手の５本の中手骨のうち隣り合う２本の中手骨どうしがなす角度に着目すると、骨塩定量分析のために普通に手を放射線検出体の上に置いた状態では、第１中手骨と第２中手骨とがなす角度が、その他の中手骨の組合せにおける角度と比べて顕著に大きくなっている。そこで、本発明の方法において特に、一つの手の第１中手骨、第２中手骨および分析対象となる骨が存在する領域をテンプレートとし、上記手が撮影された放射線画像のテンプレートと対応する部分を探索する際に、第１および第２の中手骨どうしが所定角度（例えば２０°）以上の角度をなしていることに基づいて、分析対象となる骨を認識すれば、この認識が極めて高い確率で正しく行われるようになる。

他方、本発明による画像処理装置は、
骨部を撮影した放射線画像から該骨部の骨塩定量を求める骨塩定量分析方法のために、該放射線画像から、分析対象となる骨の中の所定領域を分析対象部位として特定する画像処理装置であって、
手が撮影された放射線画像において分析対象となる骨の部分を特定する手段と、
特定された分析対象となる骨の中の特徴点を特定する手段と、
特定されたこの特徴点に基づいて前記所定領域を特定する手段とを備えたことを特徴とするものである。

また、本発明による記録媒体は、
骨部を撮影した放射線画像から該骨部の骨塩定量を求める骨塩定量分析方法のために、該放射線画像から、分析対象となる骨の中の所定領域を分析対象部位として特定する画像処理方法をコンピュータに実行するためのプログラムを記録したコンピュータ読取り可能な記録媒体であって、
前記プログラムが、手が撮影された放射線画像において分析対象となる骨の部分を特定する手順と、
特定された分析対象となる骨の中の特徴点を特定する手順と、
特定されたこの特徴点に基づいて前記所定領域を特定する手順とを有することを特徴とするものである。

本発明の一実施形態の方法を実施する骨塩定量分析システムを示す概略構成図 骨塩定量分析用に撮影された放射線画像の例を示す概略図 骨塩定量分析のために画像信号が抽出される所定領域を説明する図 上記所定領域における放射線画像の濃度プロファイル例を示す図 ある撮影装置で撮影された放射線画像におけるアルミスロープの厚さ（座標）と、そこを透過した放射線による像濃度との関係を、放射線管球の管電圧の値毎に示すグラフ 放射線管球の管電圧の違いによる画像濃度の違いを補正する処理を説明する図 管電圧とΣＧＳ／Ｄ値との関係を撮影装置の種類毎に示すグラフ 別の撮影装置で撮影された放射線画像におけるアルミスロープの厚さ（座標）と、そこを透過した放射線による像濃度との関係を、放射線管球の管電圧の値毎に示すグラフ さらに別の撮影装置で撮影された放射線画像におけるアルミスロープの厚さ（座標）と、そこを透過した放射線による像濃度との関係を、放射線管球の管電圧の値毎に示すグラフ 本発明の一実施形態による分析対象部位の特定方法の流れを示すフローチャート

以下、本発明の実施形態について図面を参照して詳細に説明する。図１は、本発明の一実施形態の方法により分析対象部位を特定する骨塩定量分析システムの概略構成を示すものである。本システムは前述したＤＩＰ法により骨塩定量分析を行うものであり、図示のように、骨塩定量分析対象の骨部を含む被写体の放射線画像を撮影する第１撮影装置１０と、この第１撮影装置１０で放射線画像が撮影された記録媒体から放射線画像を読み取って、その放射線画像を示すデジタル画像信号Ｐｃｒを出力する読取装置２０と、骨塩定量分析対象の骨部を含む被写体の放射線画像を撮影し、直接、その放射線画像を示すデジタル画像信号Ｐｄｒを出力する第２撮影装置３０と、上記デジタル画像信号ＰｃｒあるいはＰｄｒに基づいて分析対象の骨部の骨塩量を求める信号処理装置４０と、上記信号処理手段４０に各種指示を与えるための入力部５０と、骨塩定量分析結果を表示する表示部６０とから構成されている。

第１撮影装置１０は一例として、特開平８−２６６５２９号公報、特開平９−２４０３９号公報等に示される放射線検出体としての蓄積性蛍光体シートに被写体の放射線画像情報を蓄積記録するものであり、ここでは特に、蓄積性蛍光体シートを収容したカセッテ１１が用いられる。すなわち本装置１０は、カセッテ１１が略水平状態に載置される撮影台１４と、このカセッテ１１に向けて上方から放射線Ｒ（一例としてＸ線）を照射する放射線管球１２と、この放射線管球１２の駆動を制御する撮影制御部１３とを有している。

この第１撮影装置１０において、カセッテ１１の上に被写体Ｈが載置された状態下で放射線管球１２が駆動され、そこから放射線Ｒがカセッテ１１に向けて照射されると、被写体Ｈを透過した放射線Ｒのエネルギーがカセッテ内の蓄積性蛍光体シートに蓄積される。つまり該蓄積性蛍光体シートには、被写体Ｈの透過放射線画像情報が記録される。

読取装置２０は、蓄積性蛍光体シートから上記被写体Ｈの放射線画像情報を読み取るものである。この種の読取装置については、例えば特開平５−２９７４８９号公報等に詳しい記載がなされているが、以下、基本的なことを簡単に説明する。この読取装置２０においては、カセッテ１１から取り出された蓄積性蛍光体シートがレーザ光等の読取光で２次元走査され、その読取光の照射を受けた蓄積性蛍光体シートの部分から発せられた輝尽発光光が光電変換手段により読み取られて、該シートに記録されていた放射線画像情報を示す画像信号が得られる。この画像信号は、後の信号処理のためにＡ／Ｄ変換処理されて、上記のデジタル画像信号Ｐｃｒとされる。

第２撮影装置３０は、上記第１撮影装置１０における放射線管球１２、撮影制御部１３および撮影台１４とそれぞれ同様の放射線管球３２、撮影制御部３３および撮影台３４を有するものであるが、前述のカセッテ１１に代えて放射線検出器３１に放射線Ｒが照射されるように構成されている点において第１撮影装置１０と基本的に異なっている。上記放射線検出器３１は、マトリクス状に配置された画素毎に照射放射線のエネルギーレベルに対応した放射線検出信号を出力するものであり、この検出信号はＡ／Ｄ変換処理され、被写体の透過放射線画像を示すデジタル画像信号Ｐｄｒとして出力される。

なお、以上のような放射線検出器３１としては、例えば特開平７−７２２５３号公報に記載がなされているように、放射線の照射を受けて可視光を発するシンチレータと、その可視光を検出する固体光検出素子とが積層されてなるものや、あるいは例えば特開２０１０−２０６０６７号公報に記載がなされているように、放射線の照射を受けてそのエネルギーに対応した電気信号を出力する放射線光導電層を有してなるものを適用することができる。

信号処理装置４０は、上述のデジタル画像信号ＰｃｒやＰｄｒが入力される前処理部４１と、その後段に順次接続された部位抽出部４２、濃度分析部４３、管電圧補正部４４、撮影装置特性補正部４５、骨塩定量分析部４６および表示制御部４７とを有している。さらにこの信号処理装置４０は、上記管電圧補正部４４および撮影装置特性補正部４５に接続された記憶部４８を有している。

入力部５０は、例えばキーボード５１やマウス５２等の入力手段を備えてなるものであり、それらの入力手段により、信号処理装置４０が行う処理の指示を与える。

表示部６０は、例えば液晶表示装置やＣＲＴ表示装置等の表示手段６１からなるものであり、後述のようにして入力される情報に基づいて、骨塩定量分析の結果や、撮影された被写体の放射線画像を必要に応じて表示する。

以上述べた信号処理装置４０、入力部５０および表示部６０は、例えば一般的なパーソナルコンピュータ等のコンピュータシステムから構成することができる。

次に図２を参照して、骨塩定量分析のための放射線画像の撮影について説明する。ここではまず、第１撮影装置１０における撮影について説明する。この撮影に際しては、蓄積性蛍光体シートを収容したカセッテ１１が第１撮影装置１０の撮影台１４の上に載置され、その上に被検者の左手および右手が置かれ、またそれら両手の間に基準物質としてのアルミスロープが置かれる。このアルミスロープは厚さが連続的に変化するアルミニウム製の板状部材である。なおこの種のアルミスロープに代えて、厚さが段階的に変化するアルミニウム製の板状部材が用いられてもよい。

この状態で撮影制御部１３が操作されることにより放射線管球１２が駆動され、そこから発せられた放射線Ｒが上記左手および右手、並びにアルミスロープを透過してカセッテ１１内の蓄積性蛍光体シートに照射される。なおＤＩＰ法においては通常、放射線管球１２は管電圧を５０ｋＶとして撮影を行うようにしており、本実施形態においても撮影制御部１３により、管電圧は５０ｋＶに設定される。ただし、実効管電圧は経時により低下する傾向があるので、上述のように設定しても実効管電圧が５０ｋＶにならないこともある。本実施形態では、そのために分析誤差が生じることを防止するようにしているが、その点に関しては後述する。

撮影が終了すると、カセッテ１１が第１撮影装置１０から取り出されて、読取装置２０にセットされる。読取装置２０では前述の通りにして、カセッテ１１内の蓄積性蛍光体シートからそこに蓄積記録された放射線画像情報が読み取られ、その放射線画像情報を示すデジタル画像信号Ｐｃｒが得られる。このデジタル画像信号Ｐｃｒが担持する放射線画像は、該信号Ｐｃｒを前記表示部６０に入力する等により再生表示可能であり、もし表示した場合その放射線画像は図２に示すようなものとなる。すなわちこの放射線画像ＰＲには、被検者の左手ＬＨ、右手ＲＨおよびアルミスロープＡＳが記録されている。なおアルミスロープＡＳは、前述のカセッテ１１上において、左手ＬＨおよび右手ＲＨの指先方向（図２中で上方）に行くに従って次第に薄くなる状態にしてセットされる。

上記デジタル画像信号Ｐｃｒは、それを得たカセッテ１１を示す識別情報と共に信号処理装置４０の前処理部４１に入力される。また第１撮影装置１０において前述の撮影がなされたとき、撮影制御部１３からは第１撮影装置１０を示す識別情報、カセッテ１１を示す識別情報、撮影順番を示す情報等の撮影情報Ｓｃｒが、信号処理装置４０の撮影装置特性補正部４５に入力される。

次に、信号処理装置４０における処理について説明する。この信号処理装置４０に入力されたデジタル画像信号Ｐｃｒは、まず前処理部４１において、放射線の照射ムラや、読取装置２０の読取特性のムラ等に起因する信号値の変動を補正する処理や、必要に応じて適宜なされるその他の処理を受け、次に部位抽出部４２に入力される。

部位抽出部４２は、デジタル画像信号Ｐｃｒが示す画像の中から、画像処理により自動的に、あるいは前記入力部５０による指示に基づいて、骨塩定量分析を行う部位を抽出する。ＤＩＰ法においては通常、図２に示す左手第２中手骨Ｂ２Ｌについて骨塩定量分析を行うようにしているので、本実施形態においても左手第２中手骨Ｂ２Ｌが抽出される。そしてより詳細には、この左手第２中手骨Ｂ２Ｌの全長の中央部分に有るｈ／１０の領域（図３にハッチングを付して示す領域）が抽出される。

この領域の特定は部位抽出部４２において、本発明による分析対象部位の特定方法を適用してなされるが、その点に関しては後に詳しく説明する。

次に濃度分析部４３は、上記抽出された領域の平均的な濃度を求める。より具体的にこの濃度分析部４３は、上記領域において左手第２中手骨Ｂ２Ｌを横切る方向の濃度プロファイルを求める。この濃度プロファイルは、濃度に代えて輝度を用いて示すと、図４に曲線Ｑで示すようなものとなる。なお、同図に示すＤが骨幅になる。このような濃度プロファイルは、まず上記領域において骨の長さ方向に亘って分布する例えば１０数箇所程度について求められ、次にそれらの平均的なプロファイルを演算することによって求められる。

従来は、この平均的な濃度プロファイルにおける濃度がそのままアルミスロープの厚さ（アルミ厚）に換算され、つまり放射線画像においてプロファイルの各点濃度と同じ濃度となるアルミスロープ部分の厚さが求められ、そのアルミ厚換算値の積分値（図４における斜線部）ΣＧＳを骨幅Ｄで除した値ΣＧＳ／Ｄ［単位：ｍｍＡＬ（アルミニウム）］を、骨塩量を示すＤＩＰ値としていた。このＤＩＰ値については、例えば日本骨代謝学会から性別および年齢層毎の基準値が公表されており、その基準値の１００〜８０％の範囲に有れば骨塩量は正常範囲にある、といった診断が下されるようになっている。

ただし上記のＤＩＰ値＝ΣＧＳ／Ｄは、第１撮影装置１０のように蓄積性蛍光体シートを撮影記録媒体とする撮影装置に代えて、蓄積性蛍光体シートとは放射線吸収特性が異なる放射線検出器３１を用いる第２撮影装置３０が撮影に使用されたり、あるいは実効管電圧が５０ｋＶ以外になっていた場合は、同じ骨部を撮影しても、以上述べたようにして求められた場合とは異なる値を示すことがある。上に述べた基準値は、蓄積性蛍光体シートを用い、管電圧を５０ｋＶに設定して放射線画像を撮影したときのＤＩＰ値に対して定められたものであるので、この基準値を適用して骨塩量に関する診断を下すのであれば、上述のように異なる値を示すＤＩＰ値を、蓄積性蛍光体シートを用いて管電圧を５０ｋＶに設定したときのＤＩＰ値に相当するように補正する必要が有る。

以下、その補正について説明する。本発明者は、第１撮影装置１０のように蓄積性蛍光体シートに放射線画像を撮影記録する撮影装置においてアルミスロープを撮影したとき、その厚さと、撮影された放射線画像における濃度との関係が、放射線管球の管電圧に応じてどのように変化するか調べた。図５はその結果を示すものであり、これは骨塩定量分析するための放射線画像撮影に先行して、予め求められている。

なおこの図５において、横軸がアルミスロープの厚さと一義的に対応するアルミスロープの長手方向位置の座標であり、縦軸が放射線画像における濃度（相対値）である。そしてここに示す１０本の特性曲線は、上から順に管電圧が４７ｋＶ、４８ｋＶ、４９ｋＶ、５０ｋＶ、５１ｋＶ、５２ｋＶ、５３ｋＶ、５４ｋＶ、５５ｋＶ、５６ｋＶのときのものである。ここに挙げた管電圧の値は、単に撮影制御装置において設定しただけでなく、管電圧計で測定して設定通りの値になっていることを確認したものである。

同図に示される通り上記特性曲線の傾き（濃度勾配）は、管電圧の値毎に明確に相違している。このことを利用して本実施形態では、図１の管電圧補正部４４がまず、アルミスロープを示すデジタル画像信号Ｐｃｒから所定座標間の特性曲線の傾きに基づいて、第１撮影装置１０における撮影時の実効管電圧が何ｋＶであったかを求める。そのために記憶部４８には、上記傾きと管電圧との対応関係が記憶されており、管電圧補正部４４はデジタル画像信号Ｐｃｒから求めた傾きに対応する管電圧を読み出す。そして次に管電圧補正部４４は、その読み出した管電圧と管電圧５０ｋＶとに基づいて、デジタル画像信号Ｐｃｒを補正する。なお上記傾きを知るためには、アルミスロープの少なくとも２点に関して前記座標と濃度との関係を求めればよい。

上記補正は、例えば実効管電圧が４８ｋＶであった場合を例に取ると、図６に示すようになされる。すなわち、アルミ厚（もしくは骨部の放射線吸収特性）が図示のようにＥ１であった場合、管電圧が所定値の５０ｋＶであったなら本来画像濃度はＤｅ１’であるところ、実効管電圧が４８ｋＶになっているため画像濃度がＤｅ１になっているので、画像濃度Ｄｅ１を示しているデジタル画像信号Ｐｃｒを、画像濃度Ｄｅ１’を示す値に補正する。同様にして、例えば画像濃度Ｄｅ２を示しているデジタル画像信号Ｐｃｒを、画像濃度Ｄｅ２’を示す値に補正する。この補正前と補正後のデジタル画像信号Ｐｃｒの対応関係は、実効管電圧の値毎にＬＵＴ（ルックアップテーブル）の形で記憶部４８に記憶されており、管電圧補正部４４は求めた実効管電圧に関するＬＵＴを参照して、補正後の、つまり補正値とすべきデジタル画像信号Ｐｃｒの値を求める。この補正は、前述した平均的な濃度プロファイルを示すデジタル画像信号Ｐｃｒ、およびアルミスロープの部分を示すデジタル画像信号Ｐｃｒに関して全て行われ、その補正後のデジタル画像信号Ｐｃｒ’が撮影装置特性補正部４５に入力される。

前述した通り撮影装置特性補正部４５には、第１撮影装置１０の撮影制御部１３から、第１撮影装置１０を示す識別情報、カセッテ１１を示す識別情報、撮影順番を示す情報等の撮影情報Ｓｃｒが入力されている。撮影装置特性補正部４５はこの撮影情報Ｓｃｒに基づいて、あるいはデジタル画像信号Ｐｃｒに付帯された撮影情報に基づいて、入力されたデジタル画像信号Ｐｃｒ’の元のデジタル画像信号Ｐｃｒが第１撮影装置１０で生成されたものであると判別した場合は、入力されたデジタル画像信号Ｐｃｒ’を素通りさせて骨塩定量分析部４６に入力させる。

骨塩定量分析部４６は入力されたデジタル画像信号Ｐｃｒ’から前述のＤＩＰ値＝ΣＧＳ／Ｄを求める。すなわち骨塩定量分析部４６は、デジタル画像信号Ｐｃｒ’が示す平均的な濃度プロファイル（図４のＱ）における濃度をアルミスロープの厚さ（アルミ厚）に換算し、そのアルミ厚換算値の積分値ΣＧＳを骨幅Ｄで除した値ΣＧＳ／ＤをＤＩＰ値とする。骨塩定量分析部４６は、こうして求めたＤＩＰ値＝ΣＧＳ／Ｄを示す情報を表示制御部４７に入力する。表示制御部４７はこのＤＩＰ値を表示部６０の表示手段６１において表示させる。

以上のようにして表示手段６１に表示されるＤＩＰ値は、管電圧補正部４４による補正がなされた後のデジタル画像信号Ｐｃｒ’に基づくものであるから、管電圧を５０ｋＶに設定して得られたときのＤＩＰ値と等しくなる。よって、前述の基準値を利用してなされる骨塩量に関する診断も、信頼性が高いものとなり得る。なお表示部６０の表示手段６１においては、ＤＩＰ値の表示だけでなく、上記基準値との比較に基づく診断結果の表示、例えば基準値に対する比率の表示や、「骨粗鬆症の心配はありません」等の表示を併せて行うようにしてもよい。

ここで、上記デジタル画像信号Ｐｃｒをデジタル画像信号Ｐｃｒ’に補正する代わりに、まずデジタル画像信号Ｐｃｒに基づいてＤＩＰ値＝ΣＧＳ／Ｄを求め、その求められたＤＩＰ値を、デジタル画像信号Ｐｃｒ’から求められるＤＩＰ値に相当するように補正してもよい。

なお、先に述べた基準値は、蓄積性蛍光体シートを用い、管電圧を５０ｋＶに設定して撮影したときのＤＩＰ値に対して定められたものであるので、蓄積性蛍光体シートを用いる第１撮影装置１０で放射線画像を撮影した場合は、装置間の特性の違いによるＤＩＰ値の相違については考慮する必要がない。そこでこの場合は前述した通り、デジタル画像信号Ｐｃｒ’を撮影装置特性補正部４５は素通りさせているものである。また、実効管電圧が５０ｋＶであると認められた場合も、管電圧の違いによるＤＩＰ値の相違を補正する処理は不要であるので、管電圧補正部４４における上述の補正はなされない。

次に、図１の第２撮影装置３０を用いて放射線画像を撮影した場合に、信号処理装置４０でなされる処理等について説明する。まず、この撮影がなされた場合、第２撮影装置３０が出力するデジタル画像信号Ｐｄｒは前処理部４１に入力され、そこで前述と同様の処理を受ける。またこの撮影がなされたとき、撮影制御部３３からは、前述の撮影情報Ｓｃｒと同様の撮影情報Ｓｄｒが撮影装置特性補正部４５に入力される。

そしてこの場合も、管電圧補正部４４において、実効管電圧が５０ｋＶになっていなかったことが検出されると、第１撮影装置１０により放射線画像を撮影した場合と同様の補正処理がなされる。この管電圧補正部４４による補正処理は前述と同じものであるので、ここでは詳しい説明を省略する。なお図１においては、第２撮影装置３０から出力されたデジタル画像信号Ｐｄｒが上記補正処理を受けた場合、その処理済みのデジタル画像信号をＰｄｒ’と表している。

ただし、図５に示したような座標と濃度との関係は、撮影装置毎に固有のものとなるので、第２撮影装置３０によって撮影を行った場合のこの関係を図８に示す。ここで第２撮影装置３０の放射線検出器３１は、ＧｏＳ（ガドリニウムオキサイドサルファ）からなるシンチレータおよび固体光検出素子が積層されてなるものである。なお、図１には示していないが、上述のような放射線検出器として、ＣｓＩ（ヨウ化セシウム）からなるシンチレータおよび固体光検出素子が積層されてなるものも適用可能であり、その種の放射線検出器が適用された撮影装置を以下、第３撮影装置と称することとする。図９には、そのような第３撮影装置で撮影を行った場合の、上記座標と濃度との関係を示してある。

一方、本発明者は、ある共通の骨部（これは一定の厚さを持つアルミニウム板材に代えられてもよい）の放射線画像を第１撮影装置１０および第２撮影装置３０により管電圧を種々に変えて撮影し、その撮影された骨部について前述のＤＩＰ値＝ΣＧＳ／Ｄを求めた。図７はその結果を示すものである。この場合、ＤＩＰ値は前述した通り撮影時の管電圧の値に応じて変化するが、図示されるようにその変化特性は、撮影装置毎に異なっていることが分かった。なお同図において、第１撮影装置１０により撮影を行った場合の特性がＡ、第２撮影装置３０で撮影を行った場合の特性がＢ、そして前記第３撮影装置で撮影を行ったときの特性がＣである。

本実施形態においては、撮影時の管電圧を５０ｋＶに設定した場合のΣＧＳ／Ｄ値を求めるようにしているので、図７において、この管電圧が５０ｋＶのときの各ΣＧＳ／Ｄ値に着目すると、共通の骨部を撮影しているのにも拘わらず、撮影装置毎にΣＧＳ／Ｄ値が異なっている。これは、上記蓄積性蛍光体シートや、放射線検出器３１等の２種の放射線検出器における放射線吸収特性が互いに異なることに起因している。

ここで、図７の特性が得られた骨部とは異なる骨部を撮影した場合も、管電圧が５０ｋＶのときの３つのΣＧＳ／Ｄ値の間の比率は、図７の特性における比率と略同じになると考えられる。また、ΣＧＳ／Ｄ値は画像濃度と対応している。以上のことに鑑みて図１の撮影装置特性補正部４５は、入力されたデジタル画像信号Ｐｄｒ’の元になるデジタル画像信号Ｐｄｒが、第２撮影装置３０での撮影により得られたものであると判別した場合は、デジタル画像信号Ｐｄｒ’ のうちアルミスロープＡＳの部分に関する信号を、それが示す濃度Ｄｄ’が濃度Ｄｄ”＝ｋＤｄ’となるように変換する。なおｋは、図７において管電圧が５０ｋＶのときの特性ＢにおけるΣＧＳ／Ｄ値に対する、特性ＡにおけるΣＧＳ／Ｄ値の比率である。なお上記の判別は、撮影制御部３３から撮影装置特性補正部４５に入力される撮影情報Ｓｄｒや、あるいはデジタル画像信号Ｐｄｒに付帯される撮影情報等に基づいて行うことができる。

上述した変換処理がなされた後のデジタル画像信号Ｐｄｒ”は、骨塩定量分析部４６に入力される。骨塩定量分析部４６ではこのデジタル画像信号Ｐｄｒ”に基づいて、前述と同様にしてＤＩＰ値＝ΣＧＳ／Ｄが求められ、表示制御部４７はこのＤＩＰ値を表示部６０の表示手段６１において表示させる。なお上記の変換処理は、処理の都度演算を行って実行してもよいし、あるいは変換処理の前後の信号値の組合せをＬＵＴの形で記憶手段に記憶しておき、そのＬＵＴを参照して実行してもよい。

また、デジタル画像信号Ｐｄｒ’のうちアルミスロープＡＳの部分に関する信号を上記のように変換する代わりに、同じ効果が得られるように、反対に、デジタル画像信号Ｐｄｒ’のうち左手第２中手骨Ｂ２Ｌに関する信号だけを変換処理するようにしてもよい。またこのような変換処理は、アルミスロープＡＳまたは左手第２中手骨Ｂ２Ｌの濃度プロファイルを抽出してから、その濃度プロファイルを示すデジタル画像信号Ｐｄｒ’に対して施してもよいし、あるいは、抽出前の画像内該当エリアのデジタル画像信号Ｐｄｒ’に対して施してもよい。

こうして表示部６０に表示されるＤＩＰ値は、上記の変換処理がなされていることにより、撮影対象が同じ骨部であれば、第１撮影装置１０により放射線画像の撮影がなされた場合の値と同じ値となる。そこで、前述の基準値を利用してなされる骨塩量に関する診断も、信頼性が高いものとなり得る。

ここで本実施形態では、デジタル画像信号Ｐｄｒ’をデジタル画像信号Ｐｄｒ”に変換し、その変換後のデジタル画像信号Ｐｄｒ”からＤＩＰ値＝ΣＧＳ／Ｄを求めるようにしているが、デジタル画像信号Ｐｄｒ’からＤＩＰ値＝ΣＧＳ／Ｄを求め、その求められたＤＩＰ値＝ΣＧＳ／Ｄを図７の関係に基づいて、第１撮影装置１０により撮影した場合のＤＩＰ値に相当するように変換処理しても構わない。

また、管電圧補正部４４による補正と撮影装置特性補正部４５による補正を双方とも行う場合のために、それら双方の補正を同時に行うための変換値を規定したＬＵＴを作成して記憶手段に記憶させておき、そのＬＵＴを用いて上記双方の補正を一度に行うようにしてもよい。

なお、前述した第３撮影装置により放射線画像の撮影がなされた場合も、そこから出力されて前処理部４１に入力され、必要に応じて管電圧補正処理を受けた後のデジタル画像信号に対して、上記デジタル画像信号Ｐｄｒ’をデジタル画像信号Ｐｄｒ”に変換した処理と同様の変換処理を施せばよい。ただしその場合は、前記ｋの値として、図７において管電圧が５０ｋＶのときの特性ＣにおけるΣＧＳ／Ｄ値に対する、特性ＡにおけるΣＧＳ／Ｄ値の比率を適用する。

なお、図１に示すデジタル画像信号ＰｃｒやＰｄｒから、図２に示すアルミスロープＡＳの長さ方向に亘る（すなわちその厚さに対応した）濃度プロファイルを示すデジタル画像信号を抽出するには、単純に、矩形の放射線画像ＰＲの左右側縁と平行な方向（図２の上下方向）に並ぶ画素についての信号を抽出するのが簡単で好ましい。ただし、そのようにする場合、アルミスロープＡＳが斜めに、つまり放射線画像ＰＲの左右側縁に対して平行ではない状態に撮影されていると、アルミスロープＡＳの長さ方向に亘る濃度プロファイルを正確に示していない画像信号が抽出されることになる。このような不具合を防止するには、放射線画像ＰＲの左右側縁に対するアルミスロープＡＳの傾き角度を検出し、その角度だけ放射線画像を回転させた画像エリアを設定し、そのエリアの左右側縁と平行な方向に並ぶ画素についての信号を抽出すればよい。

次に、図３にハッチングを付して示した領域、つまり分析対象部位を特定する方法について詳しく説明する。図１０は、部位抽出部４２においてなされるこの方法の処理の流れを示すものである。以下、この図１０を参照して説明する。

部位抽出部４２は、本発明の画像処理装置において、分析対象となる骨の一例である第２中手骨の部分を特定する手段と、特定された第２中手骨の中の特徴点を特定する手段と、特定されたこの特徴点に基づいて所定領域を特定する手段とを構成するものであり、まず、デジタル画像信号ＰＣｒあるいはデジタル画像信号Ｐｄｒが担持している放射線画像から、左手を特定する（ステップＳＴ１）。この左手の特定は、図２の放射線画像ＰＲに示されているようにアルミスロープＡＳはそのギザギザ部分ＣＣが左手ＬＨの方を向くように配置するという撮影時のルールに則って、画像処理により、アルミスロープＡＳのギザギザ部分ＣＣに対向している方の手を左手として認識することによってなされる。

ここで、アルミスロープＡＳの存在位置および向きを認識する方法について、一例を挙げて説明する。例えば、アルミスロープＡＳの形状パターンを記憶手段に記憶しておき、公知のパターン認識処理を適用すれば、アルミスロープＡＳの存在の有無、および存在が確認されたときの位置を認識することができる。また、アルミスロープＡＳの向き、つまりそのギザギザ部分ＣＣがどの方向を向いているかということは、ギザギザ部分ＣＣのパターンを含めた位置合わせの手法を適用し、あるいは、アルミスロープＡＳの２つの長辺を探索して、その辺縁の合計長さがより長い方をギザギザ部分ＣＣが有る方の長辺とする、という手法を適用して特定することができる。

部位抽出部４２は次に、認識された左手に有る複数の骨の中から、第２中手骨Ｂ２Ｌ（図２参照）の部分を特定する（ステップＳＴ２）。この特定は、Ａｆｉｎ変換のような線形的な位置合わせでテンプレート探索する手法、あるいはモーフィングのような非線形的な位置合わせ、特に特開２０１１−２５５０６０号公報に記載されているような多重解像度位置合わせによるモデルフィッテイングでテンプレート探索する手法等を適用して行うことができる。

上記の多重解像度位置合わせは従来公知の手法であるが、概要を説明すると下に述べる（１）〜（５）の手順によって行われる。

（１）まず、テンプレート画像（本実施形態での第１および第２中手骨が含まれる画像）と参照画像（本実施形態での分析対象となる左手の画像）をそれぞれ複数の解像度に分解(ラプラシアンピラミッドやWavelet分解など)する。

（２）それぞれの解像度の画像に対して、画素探索しながら最大となる相互相関係数となる場所を選択することでシフトベクトルを求める。

（３）シフトベクトルを求める時に、最も低解像度の(粗い)画像で相互相関などを用いてシフトベクトルを算出し、次に低い解像度画像は、最も低解像度の画像で用いたシフトベクトルの位置を初期値にして、同様に相互相関などにより、次に低い解像度の画像に対するシフトベクトルを算出する。

（４）上記（１）〜（３）の処理を所定の解像度まで繰り返すことで、対象画像の全ての画素に対するシフトベクトルを算出する。

（５）さらに、各解像度で算出したシフトベクトルに基づき、変形対象となる各解像度画像(テンプレート画像か参照画像のどちらか一方）を変形して、再構成することで変形画像を得る。本実施形態においては、テンプレート画像を変形した方が、骨塩定量分析の精度を確保する上でより都合が良い。なお別の方法として、（４）で最終的に求めたシフトベクトルで画像を変形させる方法もある。そのような方法としては、例えば下記参考文献に示されるように、最終的なシフトベクトルをまとめて変形するものなどが挙げられる。

《参考文献》
G. J. Gang, C. A. Varon, H. Kashani, S. Richard, N. S. Paul, R. Van Metter, J. Yorkston, J. H. Siewerdsen, "Multiscale deformable registration for dual-energy x-ray imaging," Medical Physics 36:351-363 (2009).
以上のように、第１および第２中手骨画像が含まれたテンプレート画像を、参照画像に対して探索し、テンプレートを変形して位置合わせすることにより、参照画像における第１および第２中手骨の部分を特定できる。

部位抽出部４２は次に、この第２中手骨Ｂ２Ｌにおける特徴点を特定する（ステップＳＴ３）。本実施形態においてこの特徴点は、例えば図３に示す第２中手骨Ｂ２Ｌの中心点、両端点、重心点等である。このような特徴点の特定は、例えば左手の第２中手骨と第１中手骨が含まれる領域をテンプレートとし、そのテンプレートにランドマークを付け、変形により所定の特徴点となるランドマークを特定する手法や、テンプレートにモデル（骨の輪郭）を付け、輪郭を変形して所定の特徴点を特定する手法等を適用して行うことができる。

ここで、上記ランドマークを用いる手法について概略を説明する。例えば特徴点を上記中心点とする場合は、第２中手骨のテンプレートにランドマークとして中心点を定義し、テンプレートと参照画像との間で求められたシフトベクトルに基づいて中心点(１つの特定画素位置)を移動させることで特定できる。同様に、ランドマークとして両端点(別の特定画素位置)を定義すると、算出されたシフトベクトルに基づいて２つの画素位置が移動し、参照画像での第２中手骨の両端点となる。またランドマークとして、テンプレートの第２中手骨の輪郭線を定義すると、シフトベクトルにより輪郭線が移動変形し、参照画像の第２中手骨の輪郭線になる。ここから特徴点として重心点を求める場合は、この輪郭線(閉曲線)の重心を求める。

なお、上述のように左手の第２中手骨と第１中手骨が含まれる領域をテンプレートとするに当たっては、２つの中手骨どうしが所定角度（例えば２０°）以上の角度をなしていることに基づいて、第１中手骨および第２中手骨を認識するのが望ましい。そうすれば、先に詳述した理由により、この認識が極めて高い確率で正しく行われるようになる。

本発明の方法では、上述のようにして求めた特徴点に基づいて図１０のステップＳＴ５により、分析対象部位とする所定領域（図３にハッチングを付して示す領域）を特定してもよいが、本実施形態では分析対象部位をより高精度で特定するために、適宜、図１０のステップＳＴ４を加えて第２の特徴点を特定し、その第２の特徴点に基づいて所定領域を特定することも選択可能としている。以下、この処理について説明する。

ステップＳＴ４において部位抽出部４２は、ステップＳＴ２で特定された第２中手骨Ｂ２Ｌの部分において、まずその両端を特定する。この両端の特定は、例えば関節面を特定することによってなされる。より具体的には、第２中手骨Ｂ２Ｌの形状を考慮して、指先側を凸形状とし手首側を凹形状とした関節面のモデルを適用し、そのモデル形状を拘束（維持）しながら特定対象画像の関節に合うように変形させて行くモデルフィッテイングの手法により関節面が、つまりは第２中手骨Ｂ２Ｌの両端が特定される。なお、このモデルフィッテイングの手法として具体的には、先に述べた多重解像度位置合わせによるモデルフィッテイング等を適用することができる。

そして部位抽出部４２は、上記両端の中間点である第２中手骨Ｂ２Ｌの中央点を、第２の特徴点として求める。この中央点は、図３に示す骨の長さ方向の中央線ＣＬを求めるためのものである。

部位抽出部４２は次に、特定された中央点に基づいて、図３にハッチングを付して示す領域を特定する（ステップＳＴ５）。この特定は、例えば上記中央点から骨の長さ方向の中央線ＣＬを求め、この中央線ＣＬから骨の両端側にそれぞれｈ／２０の領域を設定する、等の手法によってなされる。

なお、ステップＳＴ３で求めた特徴点から上記ｈ／２０の領域を特定するには、例えば第２中手骨Ｂ２Ｌの両端点等から計算により中央点を求め、あるいは前記中心点をそのまま中央点として扱い、以後は上記と同様にすればよい。

また、以上説明したようにして第２中手骨の部分を特定する手順と、特定された第２中手骨の中の特徴点を特定する手順と、特定されたこの特徴点に基づいて分析対象部位となる所定領域を特定する手順とを有するプログラムをコンピュータ読取り可能な記録媒体に記録しておき、その記録媒体を用いて各手順をコンピュータに実行させることも可能である。

また、上述のようにして特定した特徴点は、第２中手骨Ｂ２Ｌを含む左手の放射線画像と共に、さらには四角形等で規定したＲＯＩ（関心領域）と共に、図１の表示手段６１において表示するようにしてもよい。それにより、特徴点の特定が正しくなされているかどうかを確認することが可能になる。また、そうして表示された特徴点を、図１の入力部５０からの指示に基づいて移動させたり、打ち直したり、骨部の輪郭やＲＯＩを変形させたりするようにしてもよい。

また、以上説明した実施形態では、左手の第２中手骨の中の所定領域を特定するようにしているが、この実施形態におけるのと同様の手法によって右手の第２中手骨の中の所定領域を特定することも可能である。

さらに本発明の分析対象部位の特定方法は、分析対象となる骨が第２中手骨以外である場合にも適用可能である。すなわち、例えばテンプレート上の分析対象となる骨の部分に、中心点や両端点などのランドマークを定義しておけば、前述したような第１中手骨と第２中手骨との位置関係から、それらおよび他の指を正確に認識してランドマークを変形移動する等により、第２中手骨以外の分析対象となる骨を正しく特定可能となる。

１０ 第１撮影装置
１１ カセッテ
１２、３２ 放射線管球
１３、３３ 撮影制御部
１４、３４ 撮影台
２０ 読取装置
３０ 第２撮影装置
４０ 信号処理装置
４１ 前処理部
４２ 部位抽出部
４３ 濃度分析部
４４ 管電圧補正部
４５ 撮影装置特性補正部
４６ 骨塩定量分析部
４７ 表示制御部
４８ 記憶部
５０ 入力部
６０ 表示部

Claims (12)

1. 骨部を撮影した放射線画像から該骨部の骨塩定量を求める骨塩定量分析方法において、該放射線画像から、画像処理によって分析対象となる骨の中の所定領域を分析対象部位として特定する方法であって、
   手が撮影された放射線画像において分析対象となる骨の部分を特定し、
   特定された分析対象となる骨の中の特徴点を特定し、
   特定されたこの特徴点に基づいて前記所定領域を特定するようにした、骨塩定量分析における分析対象部位の特定方法において、
   放射線画像において一つの手の第１中手骨、第２中手骨および分析対象となる骨が存在する領域をテンプレートとし、
   前記手が撮影された放射線画像の前記テンプレートと対応する部分を探索し、
   この探索により、第１および第２の中手骨どうしが所定角度以上の角度をなしていることに基づいて、分析対象となる骨を認識する、
   ことを特徴とする、骨塩定量分析における分析対象部位の特定方法。
2. 前記分析対象となる骨の部分を特定するために、前記テンプレートと前記手が撮影された放射線画像とに基づいて分析対象となる骨の部分を特定することを特徴とする請求項１記載の骨塩定量分析における分析対象部位の特定方法。
3. 前記分析対象となる骨の中の特徴点を特定するために、前記テンプレートと前記手が撮影された放射線画像とに基づいて前記特徴点を特定することを特徴とする請求項１または２記載の骨塩定量分析における分析対象部位の特定方法。
4. 前記分析対象となる骨の両端を特定し、これらの両端に基づいて前記特徴点を特定することを特徴とする請求項１から３いずれか１項記載の骨塩定量分析における分析対象部位の特定方法。
5. 前記分析対象となる骨の両端を、該分析対象となる骨の関節面を認識することによって特定することを特徴とする請求項４記載の骨塩定量分析における分析対象部位の特定方法。
6. 手が撮影された放射線画像において、手と共に撮影された骨塩定量分析用の標準物質を認識し、認識されたこの標準物質の位置および／または形状に基づいて、撮影された手が左右どちらの手であるか特定することを特徴とする請求項１から５いずれか１項記載の骨塩定量分析における分析対象部位の特定方法。
7. 前記特徴点として、前記分析対象となる骨の長さ方向の中央点を特定することを特徴とする請求項１から６いずれか１項記載の骨塩定量分析における分析対象部位の特定方法。
8. 前記所定領域として、前記分析対象となる骨の長さ方向の中央部分に有る、該分析対象となる骨の長さの１／１０に相当する領域を特定することを特徴とする請求項１から７いずれか１項記載の骨塩定量分析における分析対象部位の特定方法。
9. 前記特徴点を、前記分析対象となる骨の放射線画像と共に表示手段に表示させることを特徴とする請求項１から８いずれか１項記載の骨塩定量分析における分析対象部位の特定方法。
10. 前記所定領域を、前記分析対象となる骨の放射線画像と共に表示手段に表示させることを特徴とする請求項１から９いずれか１項記載の骨塩定量分析における分析対象部位の特定方法。
11. 骨部を撮影した放射線画像から該骨部の骨塩定量を求める骨塩定量分析方法のために、該放射線画像から、分析対象となる骨の中の所定領域を分析対象部位として特定する画像処理装置であって、
    手が撮影された放射線画像において分析対象となる骨の部分を特定する手段と、
    特定された分析対象となる骨の中の特徴点を特定する手段と、
    特定されたこの特徴点に基づいて前記所定領域を特定する手段とを備えてなる画像処理装置において、
    放射線画像において一つの手の第１中手骨、第２中手骨および分析対象となる骨が存在する領域を含むテンプレートを用いて、前記手が撮影された放射線画像の該テンプレートと対応する部分を探索する手段と、
    前記探索により、第１および第２の中手骨どうしが所定角度以上の角度をなしていることに基づいて、分析対象となる骨を認識する手段とを有することを特徴とする画像処理装置。
12. 骨部を撮影した放射線画像から該骨部の骨塩定量を求める骨塩定量分析方法のために、該放射線画像から、分析対象となる骨の中の所定領域を分析対象部位として特定する画像処理方法をコンピュータに実行するためのプログラムを記録したコンピュータ読取り可能な記録媒体であって、
    前記プログラムが、
    手が撮影された放射線画像において分析対象となる骨の部分を特定する手順と、
    特定された分析対象となる骨の中の特徴点を特定する手順と、
    特定されたこの特徴点に基づいて前記所定領域を特定する手順とを有し、さらに、
    放射線画像において一つの手の第１中手骨、第２中手骨および分析対象となる骨が存在する領域を含むテンプレートを用いて、前記手が撮影された放射線画像の該テンプレートと対応する部分を探索する手順と、
    前記探索により、第１および第２の中手骨どうしが所定角度以上の角度をなしていることに基づいて、分析対象となる骨を認識する手順とを有することを特徴とするコンピュータ読取り可能な記録媒体。

Images