JP4545807B2 - 骨番号認識装置およびそのプログラム - Google Patents

Description

本発明は、医用画像より骨の解剖学的な位置を解析するための骨番号認識装置およびそのプログラムに関するものである。

近年、医療分野において、人体を胸から足の付け根までスキャンした多数の断層画像からなる１シリーズの三次元画像を取得して診断を行っているが、解剖学的な位置を把握するために椎骨や肋骨の位置の自動検出が欠かせない技術となってきている。

そこで、椎骨を分離する手法として、例えば、断層画像に対してモフォロージ演算を行うことにより、椎体終板と椎間を抽出し、椎骨を分離する手法が非特許文献１に提案されている。また、断層画像に対してモフォロージ演算を行うことにより、椎間部候補領域を検出し、検出された椎間部候補領域画素の投影値を用いて椎間領域を検出して、椎骨を分離する手法が非特許文献２に提案されている。

ところで、解剖学的な位置を正確に把握するためには椎骨の番号を知ることが重要である。しかし、椎骨はお互いに似通っているため、それぞれの椎骨だけから椎骨番号を特定するのは難しい。そこで、肋骨の第１２肋骨または第１肋骨が撮影されている画像であれば、肋骨番号を正確に知ることができるので、椎骨と第１２肋骨、または、椎骨と第1肋骨（あるいは、その両方）との接続から第１２胸椎か第１胸椎を特定してから他の椎骨の番号を特定している。
"マルチスライスCT画像を用いた骨粗しょう症診断支援アルゴリズムの構築",電子情報通信学会技術研究報告MI2005-90（2006-01）,p97-100 "Using mathematical morphology for the anatomical labeling of vertebrae from 3D CT-scan images", Computerized Medical Imaging and Graphics, Vol. 31, pp. 141-156, 2007

しかしながら、従来手法のように、第１２肋骨や第１肋骨と椎骨との接続関係から椎骨の番号を特定する手法では、第１２肋骨と第１肋骨が撮影されていない場合には、椎骨の番号を特定することは出来ない。また、肋骨自体の番号も特定することはできない。

そこで、本発明は、上記事情に鑑みてなされたものであり、第１肋骨と第１２肋骨が撮影されていない三次元画像から椎体の番号や肋骨の番号を特定することができる骨番号認識装置およびそのプログラムを提供することを目的とするものである。

本願発明の骨番号決定装置は、人体を撮影して得られた複数の断層画像からなる三次元画像を記憶する画像記憶手段と、
前記三次元画像より胸骨体を示す胸骨体領域を抽出する胸骨体領域抽出手段と、
前記三次元画像より個々の肋骨それぞれを示す肋骨領域を１以上抽出する肋骨領域抽出手段と、
前記胸骨体領域抽出手段により抽出された胸骨体領域のうちの胸骨体の下端部に対応する位置と各肋骨領域の位置から該肋骨領域の肋骨番号を決定する肋骨番号決定手段とを備えることを特徴とするものである。

また、本願発明のプログラムは、コンピュータを、
人体を撮影して得られた複数の断層画像からなる三次元画像より胸骨体を示す胸骨体領域を抽出する胸骨体領域抽出手段と、
前記三次元画像より前記人体の個々の肋骨それぞれを示す肋骨領域を１以上抽出する肋骨領域抽出手段と、
前記胸骨体領域抽出手段により抽出された胸骨体領域のうちの胸骨体の下端部に対応する位置と各肋骨領域の位置から該肋骨領域の肋骨番号を決定する肋骨番号決定手段として機能させることを特徴とするものである。

「肋骨領域」とは、複数の肋骨のうちの一つの肋骨が撮影された領域をいう。

また、前記骨番号決定装置が、前記三次元画像より前記人体の個々の椎骨それぞれを示す椎骨領域を１以上抽出する椎骨領域抽出手段と、
前記肋骨領域の位置と該肋骨領域の前記肋骨番号決定手段により決定された各肋骨番号、および、各椎骨領域の位置から該椎骨番号を決定する椎骨番号決定手段とをさらに備えるようにしてもよい。

あるいは、前記プログラムが、
前記三次元画像より前記人体の個々の椎骨それぞれを示す椎骨領域を１以上抽出する椎骨領域抽出手段と、
前記肋骨領域の位置と該肋骨領域の前記肋骨番号決定手段により決定された各肋骨番号、および、各椎骨領域の位置から該椎骨番号を決定する椎骨番号決定手段としてさらに機能するようにしてもよい。

「椎骨」は、脊柱を形成する骨であり、個々の椎骨は隣接する椎骨に椎体部で椎間板を介して連結するものである。

「椎骨領域」とは、脊柱を形成する複数の椎骨のうちの一つの椎骨が撮影された領域をいう。

また、前記肋骨番号決定手段が、前記胸骨体の下端部に対応する位置より上にある胸骨体領域につながる肋骨領域の中で、最も下方にある肋骨領域の肋骨番号を第５番目として決定するものであってもよい。

さらに、前記椎骨番号決定手段が、前記椎骨領域の体軸方向に沿った中心線から所定の幅以内に存在する肋骨領域と前記椎骨領域とを体軸に投影したときの各椎骨領域との重なりが大きい肋骨領域の肋骨番号を該椎骨領域の椎骨番号として決定するものであってもよい。

本発明によれば、三次元画像より胸骨体が撮影された胸骨体領域と肋骨領域とを抽出して、胸骨体領域のうちの胸骨体の下端部に対応する位置と各肋骨領域の位置から肋骨番号を決定することにより、肋骨が１２本全て撮影されていない場合であっても、胸骨体の位置と肋骨の位置から肋骨番号を特定することが可能である。

また、胸骨体領域の下端部の位置から決められた肋骨番号とその肋骨領域の位置と椎骨領域の位置とから、椎骨番号を決定するようにしたことで、肋骨が１２本全て撮影されていない場合であっても、椎骨番号を特定することが可能である。

さらに、上記のように胸骨体領域のうちの胸骨体の下端部に対応する位置と各肋骨領域の位置から決められた肋骨番号を基準に、体軸方向の椎骨領域と肋骨領域の重なりから椎骨番号を特定することにより正確に椎骨番号を特定することができる。

以下、図面を参照しつつ、本発明の実施の形態について説明する。図１に、本発明の骨番号決定装置の概略構成図を示す。

骨番号決定装置１は、補助記憶装置に読み込まれたアプリケーションプログラムがコンピュータ（例えば、ワークステーション等）上で実行されることにより実現される。骨番号決定プログラムは、ＣＤ−ＲＯＭ等の記録媒体として配布され、記録媒体からコンピュータにインストールされても良いし、インターネット等のネットワーク経由で接続されたサーバからダウンロードされた後にインストールされてもよい。

本願発明の骨番号決定装置１は、人体を撮影して得られた複数の断層画像からなる三次元画像を記憶する画像記憶手段１０と、胸骨体領域を抽出する胸骨体領域抽出手段１１と、肋骨領域を抽出する肋骨領域抽出手段１２と、胸骨体領域の位置と各肋骨領域の位置から各肋骨領域の肋骨番号を決定する肋骨番号決定手段１３と、椎骨領域を抽出する椎骨領域抽出手段１４と、各肋骨領域の肋骨番号とその肋骨領域の位置と各椎骨領域の位置から各椎骨番号を決定する椎骨番号決定手段１５とを備える。

画像記憶手段１０は、ハードディスクなどの記憶装置であり、断層画像撮影装置で１回の検査で取得された１シリーズの複数の断層画像で構成される三次元画像を記憶する。断層画像撮影装置の具体例としては、ＣＴ（Computed Tomography）装置、ＭＲＩ（Magnetic Resonance Imaging）装置、ＰＥＴ（Positron Emission Tomography）装置などが挙げられる。以下、ＣＴ装置で撮影された三次元画像について具体的に説明する。また、以下単に断層画像という場合にはアキシャル像を指すものとして説明する。

椎骨領域抽出手段１４は、椎体位置決定手段２０を備え、椎体位置決定手段２０で決定された椎体の端の位置の情報を用いて、セグメンテーションされた椎体の端の位置から脊柱を各椎骨に分離した椎骨領域を抽出する。

椎体位置決定手段２０は、たとえばＣＴ装置により取得された断層画像に映し出された椎体の両端の位置を自動的に算出するものであって、図２に示すように、椎骨中心線算出手段２１、海綿体領域推定手段２２、中心検出手段２３、椎体位置決定手段２４を備えている。

椎骨中心線算出手段２１は、一本の椎骨の横断面を示す複数の断層画像から椎骨の中心線を求めるものである。

まず、三次元画像から脊髄領域を抽出する。脊髄領域は、図３に示すように、断層画像上に典型的なパターンで表れるため、断層画像から脊髄領域を安定的に検出することが可能である。具体的には、例えば、各断層画像中の各画素を中心として所定の大きさの領域を多数設定し、機械学習手法により生成された判別器を用いて各領域が脊髄領域であるか否かを判別しながら脊髄領域を検出する。

次に、検出された複数の脊髄領域内の中心点より脊髄中心線を生成する。中心点は、脊髄領域内の略真中に存在する画素を設定する。必ずしも脊髄領域内の厳密な中心である必要はなく、脊髄領域の周囲または所定の両端から略等距離にあるような点や、重心点でもよい。

さらに、断層画像を脊髄中心線に沿って体軸方向に切り直した縦断面画像を生成する。図４のように、各断層画像から得られた脊髄領域の中心点Ｐ５を基準として、Ｙ軸から逆時計周りにα度傾いた点Ｐ３,Ｐ４を結ぶ直線上の各画素の輝度値を抽出し、各断層画像から抽出された点Ｐ３,Ｐ４を結ぶ直線（右の図のｘ軸方向に対応）上の輝度値を断層画像のスライス番号の順（右の図のｙ軸方向に対応）に積み重ねて縦断面画像Ｖを生成する。心臓部は血管等が多く存在し、安定的なパターンが得られないため、点Ｐ３,Ｐ４を結ぶ直線は、中心点と脊椎領域を通るが心臓部は通らないように設定するのがよい。

この縦断面画像Ｖから椎骨の境界線を検出する。縦断面画像Ｖ上に脊髄中心線Ｌ３が一本の曲線として現れ、脊椎領域は脊髄中心線の左側に、二本のＣＴ値（画素値）の高い骨梁線とそれに囲まれたＣＴ値（画素値）の低い海綿骨領域が現れる。そこで、図４の縦断面画像Ｖに対し、Ｘ方向の画素値を見て脊椎領域の腹側で大きな正の差分値が表れるエッジを求め、その背側で大きな負な差分値が表れるエッジを求める。脊髄中心線Ｌ３を線形変換し、脊椎領域の背側のエッジ線Ｌ２にフィッティングさせてエッジ曲線Ｌ５を求める。同様に、腹側のエッジ曲線Ｌ１についても、同様にエッジ曲線Ｌ４を求める。

図５のように算出された曲線Ｌ４と曲線Ｌ５は、脊椎領域を含む椎骨の左右の境界線となる。そして算出された左右の境界線Ｌ４，Ｌ５から、椎骨の中心線と椎骨の幅を算出することができる。

海綿体領域推定手段２２は、椎骨中心線算出手段２１により求められた椎骨中心線の周辺の領域の画素を用いて海綿体領域を推定する。撮影された三次元画像には、椎骨の病変や変形が表れている部分が存在する可能性があるので、三次元画像全体を用いるよりも、椎骨中心線の周辺の領域だけを使うことでより安定的に海綿体を検出することができる。非特許文献１および非特許文献２は最初に椎間部を抽出してセグメンテーションをおこなっているが、画像パターンから比較的安定して検出することができる海綿体領域を先に抽出し、それを用いて椎間部を求めた方が確実に求めることができる。

図６は、椎骨および椎体等を概念的に示したものである。図６に示される円筒Ｓ１は椎体を表し、円筒Ｓ１のサーフェースは椎体の側方の皮質部分に相当する。円柱の天板と底板はそれぞれ椎体の上下の終板に相当する。突起Ｓ２は骨増殖部分を表している。円筒部分Ｔ１は、脊髄の中心付近の領域を表している。そこで、その円筒の部分Ｔ１のＣＴ値（画素値）を（横方向）加算した値を求めると、図６の右に示すようなグラフになる。断層画像上の椎間部はパーシャルボリューム効果で消失していることがあるが、図６に示すグラフのように、海綿体領域に対応するＣＴ値（画素値）の低い領域（例えば、Ｋ１〜３）が現れる。しかし、突起Ｓ２の骨増殖部分の影響で、椎間部の谷が顕著ではなくなる可能性がある。そこで、重度な圧迫骨折の人体を示す三次元画像のように、海綿体領域が潰されて海綿体領域のＣＴ値（画素値）の高くなる可能性がある場合には、ＣＴ値（画素値）の加算値から判定を行うのではなく、体軸方向のエッジを検出するようにすれば、圧迫骨折が存在する場合でも海綿体領域を検出することが可能である。

中心検出手段２３は、海綿体領域推定手段２２により推定された各椎体の海綿体領域内において、海綿体領域の三次元中心を検出する。海綿体領域の三次元中心は、海綿体領域を含む複数のサンプル画像の各画素の特徴量を予め機械学習することにより求めることができる。機械学習した判別器を用いて、椎骨中心線の周辺の領域の画素の特徴量から椎骨中心線の周辺の領域の中に海綿体領域が含まれる程度を示すスコアを算出し、このスコアに基づいて海綿体領域の三次元中心を検出する。

椎体位置決定手段２４は、海綿体領域検出手段２３により検出された各椎体の海綿体領域の三次元中心に基づいて、各椎体の体軸方向の両端の位置を算出することにより椎体位置を決定する。例えば、隣り合う二つの椎体の海綿体領域の三次元中心間の中間点をその二つの椎体の椎間部とする。その他の手法としては、隣り合う二つの椎体の海綿体領域の三次元中心間のＣＴ値（画素値）の和の極大点（又はエッジ強度）を検出することによって、椎体の皮質部の位置を検出し、椎体の端を決める手法も用いても良い。

肋骨領域抽出手段１２は、人体を表す三次元画像から人体の肋骨領域を1本ずつ分離して抽出する。

例えば、米国公開公報２００６−００６２４２５号Ａ１（以下、参考文献１という）の手法を用いることができる。参考文献１には、肋骨領域を三次元画像から抽出するための高速トレーシングベースの技術が開示されている。この技術は、肋骨を抽出して、各肋骨の一本一本を表す肋骨領域を得るのに用いることができる。

あるいは、肋骨領域の濃度ヒストグラムを用いて検出することもできる。図７に、椎骨及び肋骨部分の断層画像を示す。まず、椎骨の中心線と椎骨の横幅を求める。そのために、図７中の骨領域のうち、体表横方向の中心線Ａ（すなわち、頭と足を結ぶ体軸の方向に垂直な断面における体の前後方向の中央において体の前後方向に垂直な体の左右方向の線）より下側の領域を縦軸に投影し、濃度ヒストグラムを生成する。

図８に、このようにして生成した骨領域の濃度ヒストグラムを示す。図８に示すヒストグラムでは、中央に大きなピークがあり、左右に小さなピークが形成されている。図７に示す断層画像では、肋骨の部分が小さいピークを形成し、椎骨が一番骨の多く存在するところであり一番大きいピークが形成されると考えられる。そこで、図８のヒストグラムから、中央にある最大のピークを検出し、それに対応する部分として椎骨領域を検出する。

そして、図９に示すように、検出した椎骨領域Ｂの上部に対して円Ｓを当てはめて、椎骨の中心Ｃ及び幅ｄを求める。椎骨中心Ｃと椎骨幅ｄを用いて椎骨領域、左右肋骨領域及び胸骨領域を認識する。

図１０に示すように、上で求めた椎骨中心Ｃから上方に６０°の角度をなす２本の半直線を引くとともに、椎骨中心Ｃをその中心に含むようにして体表縦方向の中心線Ａに直交する２直線を引く。この２直線間の距離は、椎骨領域の幅であり、これは上記幅ｄの１.５倍となるようにする。このように、２本の半直線と体表縦方向の中心線Ａに直交する２直線により、図１０の断層画像の骨領域を４つの領域に分割し、それぞれ椎骨領域Ｂ、右肋骨領域Ｅ１、胸骨領域Ｅ２、左肋骨領域Ｅ３とする。

このようにして椎骨、右の肋骨、胸骨、左の肋骨のように骨をその部品単位で認識する。三次元画像の全ての断層画像に対して各部品の認識を行って、各断層画像から検出された肋骨の領域をスライス番号の順番に従って積み重ねた領域を肋骨領域として抽出する。

胸骨体領域抽出手段１１は、三次元画像より胸骨体領域を抽出する。胸骨体は、図１１に示すように、各断面画像上にほぼ同じ太さで内腔（海綿体）のある長方形の領域として表れ（図１１のＡの部分を参照）、パターン認識をすることにより検出することができる。胸骨体は胸部断面上の椎骨の位置とは反対の人体の前面に近い位置に表れるので、人体の前面に近い辺りから長方形の領域を各断層画像から検出する。各断層画像から検出した長方形の領域を各断層画像のスライス番号の順番に従って積み重ねた領域を、胸骨体領域として抽出する。

肋骨番号決定手段１３は、肋骨領域が胸骨体領域とどのように接続しているかに応じて肋骨番号を決定する。図１２に示すように、解剖学的には第２肋骨〜第１０肋骨は胸骨体に繋がっている。また、第５肋骨までは、肋骨のほぼ全体が胸骨体の下端Ｌより上にある。したがって、胸骨体領域に接続している肋骨領域の中、ほぼ全体が胸骨体の下端Ｌより上にあり、かつ、最も下にある肋骨領域の肋骨番号を第５番の肋骨として特定することが出来る。この第５番の肋骨領域を基準に各肋骨領域の肋骨番号を肋骨領域の並びの順番にしたがって決定する。

椎骨番号決定手段１５は、肋骨領域の肋骨番号とその肋骨領域が椎骨領域とどのように接続しているかから、各椎骨領域の椎骨番号を決定する。

肋骨領域がどの椎骨領域と接続しているかを求める一つの手法としては、図１３に示すように、まず椎骨領域の中心線からある幅Ｗを決め（例えば、Ｗとして椎骨の幅を用いることができる。）、その範囲内の肋骨番号ｋの肋骨領域をＺ軸上に投影し、Ｚ方向に肋骨領域が表れる範囲を求める。同様に、椎骨領域をＺ軸上に投影してＺ軸における各椎骨領域の範囲を求める。各椎骨領域の範囲と重なる範囲が最も大きい肋骨番号ｋの肋骨領域が、その椎骨領域に接続する椎骨であると判定して、その椎骨領域の椎骨番号をｋとする。

次に、図１４のフローチャートに従って、骨番号決定装置の処理の流れについて説明する。

まず、肋骨領域抽出手段１２で、画像記憶手段１０に記憶されている三次元画像に撮影されている全ての肋骨領域を抽出する（ＳＴ１００）。次に、椎骨領域抽出手段１４で、三次元画像から脊髄領域を抽出し、椎体の両端の位置を自動的に算出して各椎骨領域に分離する（ＳＴ１０１）。

抽出された肋骨領域の数をカウントする（ＳＴ１０２）。肋骨領域の数が１２本あれば、第１２番目の肋骨領域を特定して（ＳＴ１０３）、肋骨番号が１２番の肋骨領域を基準に上の肋骨領域へ行くに従って肋骨番号を１１番、１０番、・・・、１番というように順に割り当てる（ＳＴ１０４）。

肋骨領域の数が１２本以下の場合には、３次元画像に腰骨が撮影されているか否かを判定して、腰骨が撮影されていれば（ＳＴ１０５）、肋骨の一番下の肋骨は撮影されているので、肋骨領域の一番下の肋骨番号を１２番と特定して（ＳＴ１０３）、肋骨番号が１２番の肋骨領域を基準に上の肋骨領域へ行くに従って肋骨番号を１１番、１０番、・・・というように順に割り当てる（ＳＴ１０４）。肋骨領域の数が１２本以下で腰骨も撮影されていないときには、頚骨が撮影されているか否かの判定を行う。頚骨が撮影されているときには肋骨の一番上の肋骨は撮影されているので（ＳＴ１０６）、肋骨領域の一番上の肋骨番号を１番と特定して（ＳＴ１０７）、肋骨番号が１番の肋骨領域を基準に下の肋骨領域へ行くに従って肋骨番号を２番、３番、・・・というように順に割り当てる（ＳＴ１０８）。

さらに、椎骨番号決定手段１５で肋骨領域の肋骨番号とその肋骨領域が椎骨領域とどのように接続しているかから、各椎骨領域の椎骨番号を決定する（ＳＴ１１１）。

肋骨領域の数が１２本以下で腰骨も頚骨も撮影されていないときには、胸骨体領域抽出手段１１で、三次元画像から胸骨体領域を抽出する（ＳＴ１０９）。肋骨番号決定手段１３で肋骨領域が胸骨体領域とどのように接続しているかに応じて肋骨番号を決定する（ＳＴ１１０）。

次に、椎骨番号決定手段１５で肋骨領域の肋骨番号とその肋骨領域が椎骨領域とどのように接続しているかから、各椎骨領域の椎骨番号を決定する（ＳＴ１１１）。

上述では、肋骨領域が１２本全部確認できるときや、第１肋骨や第１２肋骨が確認できる時には胸体骨の認識を行わずに肋骨番号を特定する場合について説明したが、胸体骨の認識を必ず行って胸体骨と第５肋骨の位置を特定して、他の肋骨の肋骨番号を特定するようにしてもよい。

以上、詳細に説明したように、肋骨が１２本全て撮影されていない場合であっても、胸体骨から肋骨番号を決定することができ、さらに、肋骨がどの椎骨に接続しているかに応じて椎骨番号を決定することができる。

骨番号決定装置の概略構成図 椎体位置決定手段の概略構成図 断層画像に撮影された椎骨の一例 脊髄中心線の算出方法を説明するための図 椎骨の左右の境界線の算出方法を説明するための図 椎体位置の決定方法を説明するための図 椎骨及び肋骨部分の断層画像の一例 骨領域の濃度ヒストグラムを表す図 椎骨の中心及び幅を表す図 断面画像を椎骨領域、右肋骨領域、胸骨領域、左肋骨領域に分けた図 断層画像に撮影された胸骨体の一例 胸体骨と肋骨との接続関係を説明するための図 椎骨と肋骨との接続関係を説明するための図 骨番号決定装置の処理の流れを示すフローチャート

符号の説明

１ 骨番号決定装置
１０ 画像記憶手段
１１ 胸骨体領域抽出手段
１２ 肋骨領域抽出手段
１３ 肋骨番号決定手段
１４ 椎骨領域抽出手段
１５ 椎骨番号決定手段
２０ 椎体位置決定手段
２１ 椎骨中心線算出手段
２２ 海綿体領域推定手段
２３ 中心検出手段
２４ 椎体位置決定手段

Claims (7)

1. 人体を撮影して得られた複数の断層画像からなる三次元画像を記憶する画像記憶手段と、
   前記三次元画像より胸骨体を示す胸骨体領域を抽出する胸骨体領域抽出手段と、
   前記三次元画像より個々の肋骨それぞれを示す肋骨領域を１以上抽出する肋骨領域抽出手段と、
   前記胸骨体領域抽出手段により抽出された胸骨体領域のうちの胸骨体の下端部に対応する位置と各肋骨領域の位置から該肋骨領域の肋骨番号を決定する肋骨番号決定手段とを備えたことを特徴とする骨番号決定装置。
2. 前記三次元画像より前記人体の個々の椎骨それぞれを示す椎骨領域を１以上抽出する椎骨領域抽出手段と、
   前記肋骨領域の位置と該肋骨領域の前記肋骨番号決定手段により決定された各肋骨番号、および、各椎骨領域の位置から該椎骨番号を決定する椎骨番号決定手段とをさらに備えたことを特徴とする請求項１記載の骨番号決定装置。
3. 前記肋骨番号決定手段が、前記胸骨体の下端部に対応する位置より上にある胸骨体領域につながる肋骨領域の中で、最も下方にある肋骨領域の肋骨番号を第５番目として決定するものであることを特徴とする請求項１記載の骨番号決定装置。
4. 前記肋骨番号決定手段が、前記胸骨体の下端部に対応する位置より上にある胸骨体領域につながる肋骨領域の中で、最も下方にある肋骨領域の肋骨番号を第５番目として決定するものであることを特徴とする請求項２記載の骨番号決定装置。
5. 前記椎骨番号決定手段が、前記椎骨領域の体軸方向に沿った中心線から所定の幅以内に存在する肋骨領域と前記椎骨領域とを体軸に投影したときの各椎骨領域との重なりが大きい肋骨領域の肋骨番号を該椎骨領域の椎骨番号として決定するものであることを特徴とする請求項２または４記載の骨番号決定装置。
6. コンピュータを、
   人体を撮影して得られた複数の断層画像からなる三次元画像より胸骨体を示す胸骨体領域を抽出する胸骨体領域抽出手段と、
   前記三次元画像より個々の肋骨それぞれを示す肋骨領域を１以上抽出する肋骨領域抽出手段と、
   前記胸骨体領域抽出手段により抽出された胸骨体領域のうちの胸骨体の下端部に対応する位置と各肋骨領域の位置から該肋骨領域の肋骨番号を決定する肋骨番号決定手段として機能させるプログラム。
7. 前記三次元画像より前記人体の個々の椎骨それぞれを示す椎骨領域を１以上抽出する椎骨領域抽出手段と、
   前記肋骨領域の位置と該肋骨領域の前記肋骨番号決定手段により決定された各肋骨番号、および、各椎骨領域の位置から該椎骨番号を決定する椎骨番号決定手段としてさらに機能させる請求項６記載のプログラム。