JP4495891B2

JP4495891B2 - 骨計測方法

Info

Publication number: JP4495891B2
Application number: JP2001501107A
Authority: JP
Inventors: 由明東; 善史原田; 典弘山田; 勉前田; 知裕太田
Original assignee: Teijin Ltd; Infocom Corp
Current assignee: Teijin Ltd; Infocom Corp
Priority date: 1999-06-03
Filing date: 2000-06-02
Publication date: 2010-07-07
Anticipated expiration: 2020-06-02
Also published as: WO2000074567A1; US6839457B1

Description

【０００１】
【発明の属する技術分野】
本発明は、被検骨あるいは関節断面撮影画像に基づいて骨形態を測定する骨計測方法に関し、更に詳しくは、被検骨あるいは関節断面の２値化画像に基づいて皮質骨と海綿骨とを分離する骨計測方法において、皮質骨と海綿骨とを自動的かつ高速かつ再現性よく分離することができ、非破壊的に３次元骨構造・骨強度などの解析の基礎となる皮質骨部と海綿骨部に分離した３次元画像が得られる骨計測方法に関する。更に、被検骨あるいは関節断面の２値化画像に基づいて皮質骨と海綿骨とを分離した３次元画像において、その分離断面部（以下、連結面とよぶ）に対し物体の境界面を定義する方法（改良マーチングキューブス法）を適用し、効率よく自動的かつ高速かつ再現性がよく面定義可能で、非破壊的に３次元骨構造・骨強度などの解析の指標となる境界面部位の３次元画像が得られる骨計測方法に関する。
【０００２】
【従来の技術】
骨の強度を評価する１つの手段として、骨構造について評価する事が指摘されている。生体の多くの骨組織は骨構造の異なる皮質骨と海綿骨に大別され、それぞれ役割が異なると考えられている。構造学的には骨に負荷がかかった場合、皮質骨部は大荷重に耐えられるように、海綿骨部は衝撃を吸収するように作用すると考えられている。
従って骨構造を評価するには皮質骨部と海綿骨部を別々かつ定量的に評価しなければならない。皮質骨と海綿骨は構造的に連続していることが多く、両者の境界を特定するのは難しく、皮質骨と海綿骨を厳密に再現よく分離できる骨計測方法が要望されている。
これに対して、特開平９−２９４７４０号公報には皮質骨部と海綿骨部を自動分離できる以下の骨計測方法が提案されている。この方法は、被検骨の断面撮影画像の２値化画像に基づいてその骨部の細線化画像を用いて分離するもので、具体的には骨部を細線化し、該細線化画像の外側骨部から皮質骨の平均骨幅を推定し、該推定値から細線化画像を膨張させて皮質骨部を求める方法である。
しかし、かかる方法では、皮質骨部あるいは海綿骨部が分断されたりする場合があり、この手法を用いて複数の被検骨を比較する場合には問題が生じる。
【０００３】
【特許文献１】
特開平９−２９４７４０号公報
【発明の開示】
【０００４】
【発明が解決しようとする課題】
本発明は、かかる問題に鑑み為されたもので、被検骨換言すれば被験者が異っても安定した皮質骨及び海綿骨の分離が可能で、さらに非破壊的に骨の３次元内部構造を構築することが可能で特に骨の３次元的評価がより正確に行える、高速かつ自動的に再現性のよく皮質骨と海綿骨が分離できる骨計測方法を第一の目的としている。さらには、被験者が異っても安定した皮質骨及び海綿骨の境界面生成が可能で、非破壊的に骨の境界を定義することができ、境界面特性を定量的に評価する特性（境界の個数、面積、周囲長、円形度）を、高速かつ自動的に再現性よく求めることができる骨計測方法を第二の目的としている。
【０００５】
【課題を解決するための手段】
上記第一の目的は以下の本発明で達成される。すなわち本発明は、被検骨あるいは関節断面の２値化画像に基づいて皮質骨と海綿骨とを分離する骨計測方法において、該２値化画像から皮質骨に囲まれた骨内部の単連結したテンプレート画像を求めるテンプレート画像抽出ステップを有し、該テンプレート画像と該２値化画像との積により皮質骨と海綿骨とを分離することを特徴とする骨計測方法である。
本発明は、上記テンプレート画像は骨部全体を抽出できるもので、従って再現性よく安定した分離ができる共に、その抽出ステップは、非常に簡単なアルゴリズムで構成できるので、高速処理でき、上記目的を達成できることを見出し、為されたものである。
本発明において、皮質骨と海綿骨に分離された分離画像を、得られた２値の分離画像と２値化前の濃淡画像との積により簡単に輝度値の分離画像とする事ができる。
また、被検骨を所定間隔で断層撮影した複数の連続した各画像から皮質骨部と海綿骨部とが分離された各部の分離画像を作成し、該分離画像を積み重ねることにより３次元の骨分離画像を生成する構成により、３次元画像が高速で得られる。そして、この分離画像に上記の輝度値換言すれば濃淡のある分離画像を用いれば、分離された皮質骨部と海綿骨部の立体画像が得られる。
このように、本発明により、非破壊的に高速かつ自動的かつ再現性のある分離された皮質骨と海綿骨の３次元画像を抽出することができるため、より正確な３次元骨関連情報の評価が可能である。
さらに、本発明は、皮質骨と海綿骨とを分離した画像をスライス方向（断面方向）へ積み上げた３次元画像上にて、その分離切断面部（以下、境界面とよぶ）に対し物体の境界面を定義する方法（改良マーチングキューブ法・特許登録番号０２７９８５６）を適用し、境界表面に三角形を連続する様に定義する画像処理手段を用いる。
通常、マーチングキューブス法は２値化された３次元物体画像に対し背景と物体との境界に多角形の面を定義するものであるが、これを皮質骨と海面骨の境界面の定義に適用するものである。この方法により、三角形で定義された境界面を定義することが可能である。
また、境界個数については三角形の頂点座標の連結状態を検索し、その繋がり性から境界面のラベル個数を数え上げる処理手段で求められる。
また、境界面積についてはラベル別三角形の面積を積算することを特徴とする骨計測方法である。
また、周囲長についてはラベル別三角形の最外側の辺の長さを積算することを特徴とする骨計測方法である。
また、円形度（複雑度）についてはラベル別三角形の面積と周囲長から境界面の複雑性を定量化し、その際、真円を１と定義することを特徴とする骨計測方法である。
このように、本発明により、非破壊的に高速かつ自動的かつ再現性のある分離された皮質骨と海綿骨の境界面画像を抽出することができるため、より正確な３次元骨関連情報の評価計測が可能である。
【０００６】
【発明を実施するための最良の形態】
以下、本発明の詳細を説明する。
以下、ラットの腰椎に適用した実施例に基いて、実施の形態共に順を追って説明する。また、ラットの大腿骨についても実施例を示し、長管骨へも適応可能であることを示す。
【０００７】
（１）皮質骨と海綿骨が分離できる骨計測方法
被検骨のラットの腰椎の断面は、前述の特開平９−２９４７４０号公報と同様に、本例では撮像手段に海綿骨の微細構造を計測するのに十分な焦点寸法及び解像度をもつ微小フォーカスＸ線断層撮影装置（以下、μＸ線ＣＴ）を用いて撮影した。なお、被検骨の撮影には、この他高解像度Ｘ線装置または核磁気共鳴画像診断装置（ＭＲＩ）またはフィルムスキャナー等の２次元情報を発生する装置等が適用可能である。
そして、撮像手段で得られた該被検骨の画像信号を以下の画像処理装置で処理して、骨計測を行う。本例の画像処理装置１０は、図１に示すように画像処理専用の画像処理ボード１２（本例ではシャープセミコンダクタ社製、ＧＰＢ−Ｋ）並びに画像記憶用に大容量のハードディスク１３を備え、画像表示用の外部モニター１４を搭載したパーソナルコンピュータ１１を用いた。そして、図２に示すフローチャートの処理手順を収納し、撮像手段から得られた被検骨の画像信号を処理して以下のようにその骨部を皮質骨部と海綿骨部に自動的に分離抽出する。
処理に先立ち、前述の特開平９−２９４７４０号公報と同様にして、まず、被検骨のラットの腰椎を、μＸ線ＣＴ装置によりその軸方向に所定間隔で連続した断層像を撮影して、その画像信号を取得する。なお、この画像は、骨梁を観察するのに十分な程度（１０〜２５μｍ）の空間解像度が得られるものがよい。本例において用いたμＸ線ＣＴ画像は、縦５１２画素×横５１２画素×高さ（軸方向）３０スライス（最大５１２スライスまで可能）であり、断面の１画素あたりの寸法は約２０μｍから２５μｍ、スライス間隔２０μｍ程度、各画素の輝度値ＣＴ（濃度値とも言う）は２^８階調で表現されている。
図３に、このμＸ線ＣＴ装置により得られた被検骨のラットの腰椎の断面の輝度値表示の濃淡画像（以下、「原画像」という）の代表例を示す。図３において、２０が骨部で、２１がその皮質骨部で、外側の輪郭線に見える外側部と内側の輪郭線に見える内側部とからなる。そして、皮質骨部２１に囲まれた骨内部の網目状の線状部分が海綿骨部２２である。図３の２３は骨髄腔部（以下、空間部）を示し、図３の２４は骨髄腔部（以下、空洞部）であり、骨部２０の外側は背景部２５である。図３の濃淡画像の原画像においては、皮質骨部２１と海綿骨部２２と空間部２３の境界が不鮮明であり、このままではこれらの分離は困難である。
この原画像をステップ１で画像処理装置１０に取り込む。そして、ステップ２の２値化ステップにおいて、この原画像を、以下のように上述の公報と同様にして判別分析法により２値化して、２値化画像に変換する。
判別分析法は、グループ内の濃度差が小さく、しかもグループ間の偏差が大きい場合に２値化の閾値を決定する際に使用することができ、本発明の対象である骨計測には好適である。
判別分析法では、画像の輝度値のヒストグラムにおいて輝度値の集合を閾値で二つのクラスに分割したと仮定して、下記する式２のクラス間分散および式３のクラス内分散を用いて式１の分散比（ｆ_０）が最大になるような閾値を求め、この閾値により２値化する。
【０００８】
【数１】
式１：ｆ_０＝δｂ^２／δｗ^２
但し、ｆ_０：分散比
δｂ^２：クラス間分散
δｗ^２：クラス内分散
【数２】
式２： δｂ^２＝ω_１ω_２（ｍ_１−ｍ_２）^２
但し、ω_１：クラス１の画素数
ω_２：クラス２の画素数
ｍ_１：クラス１の画素の輝度値の平均
ｍ_２：クラス２の画素の輝度値の平均
【数３】
式３： δｗ^２＝ω_１δ_１ ^２＋ω_２δ_２ ^２
但し、δ_１：クラス１の画素の輝度値の分散
δ_２：クラス２の画素の輝度値の分散
【０００９】
本例では、原画像の輝度値のヒストグラムを求め、このヒストグラフから閾値の候補値を選定し、それを変化させて最適値を求め、これを閾値とした。そして、この閾値より大きいものを１、小さいものを０として２値化して、２値化画像を得た。図４に図３の原画像を２値化した２値化画像（これを「２値化原画像」という。）を示す。図より、空間部２３と皮質骨部２１、海綿骨部２２との区分が明瞭となっていることが分かる。
そして、本例では２値画像処理として周知の８近傍に基づいてこの２値化原画像から以下のテンプレート画像抽出ステップにより海綿骨部２２の存在領域である外側と内側の皮質骨部２１で囲まれたその内側の領域をテンプレート画像として抽出する。
まず、この２値化原画像には通常孤立点等のノイズがあるので、ステップ３の平滑化ステップにおいて、以下のようなノイズ処理を行う。本例では、ノイズ処理として、該２値化原画像から８近傍の孤立点消去を行った。８近傍の孤立点消去とは、ある画素に善目した場合、その近傍８方向に関する画素を走査し全て輝度値が０の場合だけ該着目画素を０とする処理方法である。ノイズ処理後の画像を図５Ａに示す。
次に血管等により分断されている皮質骨２２等の該分断部分を仮想的に接続する接続処理を行う。この接続処理は、ステップ４の膨張ステップで画像の境界点を１層分太らせる膨張処理をＮ回施すことで行った。この回数は、該分断部分が接続される回数であればよく、本例ではその分断部分の間隔からＮ＝１０とした。処理後の画像を図５Ｂに示す。図より青部画像の中には皮質骨部２１で囲まれた内部領域の空間部２３に対応する穴部２６と内側の皮質骨部２１に囲まれた空洞部２４に対応する背景穴部部２７が残っていることが判る。
ところで、海綿骨部２２が存在する骨の内部領域画像すなわちテンプレート画像を求めるには、内部領域を全ての穴部２６を埋めて穴部の無い単連結した領域とする必要がある。従って、画像処理としては、内部領域の穴部２６のみを埋めて消去し、背景穴部２７はそのまま残す必要がある。この該背景穴部２７を残さない処理をすれば、大腿骨のような長管骨にも十分応用が可能である。
この穴埋め処理をステップ５の穴埋めステップで以下のように行った。前述の膨張ステップで得られた膨張画像を先ず反転して反転画像を得る。得られた反転画像を図５Ｃに示す。この反転画像において、各穴部２６及びは背景穴部２７に対応する各白部の内、残す背景穴部２７は最大面積となっているので、各白部の面積を求め、最大面積の白部を除く背景部除去処理を行って、背景穴部２７に対応する図の中心部の楕円状白部を消去した。なお、この背景部除去処理は、空洞部２４がある本例の被検骨の腰椎に必要なもので、かかる空洞部のない大腿骨等では必要なく、省略できる。この処理で得られた穴埋め用画像を図５Ｄに示す。図より、内部領域の全ての穴部２６に対応する白部を有することが判る。そして、この得られた図５Ｄの穴埋め用画像と先の膨張ステップで得られた図５Ｂの膨張画面の和（ＯＲ）画像を求める。この和処理で得られた画像を図５Ｅに示す。図より、穴部２６が無く、単連結した皮質骨部２１に囲まれた骨部２０の全領域が全面連続した膨張連結画像が得られたことが判る。
【００１０】
ところで、この膨張連結画像は前述の通りＮ回本例では１０回膨張処理を行ったものであり、大きさの面では２値化原画像を拡大したものとなっている。そこで、ステップ６の収縮ステップで、２値化原画像に相当する大きさまで収縮する。本例では、画像の境界点を取り除いて１層分小さくする収縮処理を（Ｎ＋１）回具体的には１１回行った。一般的に膨張させた画像を寸法面で原画像に戻すには膨張処理の回数と同数回の収縮処理を行えばよい。ところが、これらの処理は、周知の通り、単純な拡大縮小ではなく、平滑化、穴埋め等の所定の機能を有する。そこで本例では、テンプレート画像として皮質骨部２１で囲まれたその内部領域を抽出することが重要であり、皮質骨部２１の外側が残らないように膨張処理回数より１回多い回数の収縮処理を行った。これにより、皮質骨部２１の形状に影響されること無く、安定な内部領域の抽出が可能となった。この収縮ステップで得られた画像を図５Ｆに示す。この連結収縮画像は、以上の処理から輪郭部の皮質骨部２１を含むこれで囲まれた骨部２０の穴部の無い全面連続画像になっている。
これに対して、海綿骨部２２を抽出するテンプレート画像は皮質骨部２１を含まない骨部２０の内部の全面連続画像となる。そこで、皮質骨部２１を除くため、ステップ７の以下の皮質骨部除去ステップを施す。すなわち、ステップ６で得られた連結収縮画像からステップ３で得られた平滑処理された２値化原画像（図５Ａ）を差し引いてその差画像を求める。その画像を図５Ｇに示す。得られた差画像は、骨部２０の全面連結画像から皮質骨部２１と海綿骨部２２が除去された空間部の骨空間部画像となっている。
ところで、この骨空間部画像は、見方を変えると、海綿骨部２２の網状部を埋めれば、皮質骨部２１で囲まれた骨部２０から皮質骨部２１を除いたその内部領域すなわち海綿骨部２２の存在する領域の画像であるテンプレート画像となっていることが判る。そこで、このテンプレート画像を得るためステップ８の海綿骨部２２の間隙を除去する以下の間隙除去ステップを行う。すなわち、具体的には、図５Ｇで海綿骨部２２を差し引いて残され間隙を接続して全面連続して１つの領域とするため以下の接続する処理を施す。この接続処理は、先ずＭ回の前述の画像膨張処理を施して間隙を接続する。得られた画像を図５Ｈに示す。図のように間隙が除去されて全面連続した単連結の内部領域画像となっている。次いで、寸法を戻すために同じＭ回の前述画像収縮処理を施す。結果を図５Iに示す。図が海綿骨部２２を検出するための海綿骨部２２が存在する皮質骨部２１に囲まれたその内部領域の画像であるテンプレート画像である。なお、本例の場合には、接続する間隙幅からＭ＝１０とした。以上により海綿骨部２２を検出するための海綿骨用テンプレート画像が得られたことになる。
そして、ステップ１０の海綿骨部検出ステップで、この海綿骨用のテンプレート画像と前述の図５Ａに示した平滑化２値化原画像との積（ＡＮＤ）画像を求めることにより、２値化海綿骨画像を検出することができる。結果を図５Ｊに示す。次いで、ステップ１１の皮質骨部検出ステップで、前述の図５Ａに示した平滑化２値化原画像からこの２値化海綿骨画像を差し引いてその差画像を求めることにより、２値化皮質骨画像が検出できる。得られた画像を図５Ｋに示す。共に良好な画像が得られた。
以上の通り、本発明により、簡単な処理により、自動的に皮質骨部と海綿骨部に分離できることが判る。
また、上記の方法によって求めた該２値化皮質骨画像および２値化海綿骨画像から簡単な処理により入力原画像上の輝度値を保持し、他のバックグランド領域の輝度値を０にした半閾値化画像を求めることができる。具体的には、該２値化皮質骨部及び海綿骨部の各画像と図３に示す原画像とから輝度値の最小値画像を求めることによりそれぞれの半閾値化画像を得ることができる。この結果を図６に示す。図６Ａが海綿骨部、図６Ｂが皮質骨部の半閾値化画像である。
更に、図７に示すように、軸方向に所定のスライス間隔で撮影した各連続画像を以上のように処理して生成された図６の皮質骨部と海綿骨部の半閾値化画像データを各部毎にスライス方向に積み重ねていくことにより皮質骨部と海綿骨部の３次元画像を得ることができる。この結果を図８に示す。図８Ａが皮質骨部、図８Ｂが海綿骨部の３次元画像であり、非常に良好な立体画像が得られることが判る。
【００１１】
（２）皮質骨及び海綿骨の境界面生成および境界面特性の定量的評価方法
被検骨のラットの腰椎の断面は、特開平９−２９４７４０号公報と同様に、本例では撮像手段に油綿骨の微細構造を計測するのに十分な焦点寸法及び解像度をもつ微小フォーカスＸ線断層撮影装置（以下、μＸ線ＣＴ）を用いて撮影した。なお、被検骨の撮影には、この他高解像度Ｘ線装置または核磁気共鳴画像診断装置（ＭＲＩ）またはフィルムスキャナー等の２次元情報を発生する装置等が適用可能である。
そして、撮像手段で得られた該被検骨の画像信号を以下の画像処理装置で処理して、骨計測を行う。本例の画像処理装置１０は、図１に示すように画像処理専用の画像処理ボード１２（本例ではシャープセミコンダクタ社製、ＧＰＢ−Ｋ）並びに画像記憶用に大容量のハードディスク１３を備え、画像表示用の外部モニター１４を搭載したパーソナルコンピュータを用いた。そして、図９に示すフローチャートの処理手順を収納し、「発明の実施の形態」（１）記載の方法で骨部を皮質骨部と海綿骨部に自動的に分離抽出する。
まず、被検骨のラットの腰椎を、μＸ線ＣＴ装置によりその軸方向に所定間隔で連続した断層像を撮影して、画像処理装置１０へ取り込み、皮質骨と海面骨を分離した画像（以下、分離処理画像）を取得する。なお、この画像は、骨梁を観察するのに十分な程度（１０〜２５μｍ）の空間解像度が得られるものがよい。
本例において用いたμＸ線ＣＴ画像は、縦５１２画素×横５１２画素×高さ（軸方向）２００スライス（最大５１２スライスまで可能）であり、断面の１画素あたりの寸法は約２０μｍから２５μｍ、スライス間隔２０μｍ程度、各画素の輝度値ＣＴ（濃度値とも言う）は２^８階調で表現されている。この画像を、図９のステップ１の原画像入力により画像処理装置１０に読み込む。図１０は、ラットの腰椎の原画像である。
図９のステップ２の２値化ステップを実行する。その際、２値化の方法は、判別分析法を用いる。
皮質骨（図１１）と海面骨（図１２）の分離処理を図９のステップ３で行い、スライス方向（断面方向）に必要枚数分の画像を格子状に積み重ねる。代表例を図１３、図１４に示す。皮質骨は図１３、海面骨は図１４のような３次元画像になる。これらの画像（図１３と図１４）のＯＲ画像が全体骨画像である。図１５の様に２次元画像上の各ピクセルに対し厚みを持たせたものがボクセルである。
ここで、一般的な物体表面形状の定義はマーチングキューブス法を使用するのが一般的である。マーチングキューブス法とはボクセル集合の局所的な配置からボクセル集合を覆う面素を推定し、ボクセルの表面を多角形（本件では三角形）で近似して、物体の状面を構成する方法である。具体的には、図１６のように、８つの隣り合うボクセルを１つのキューブ上の頂点と考え、その頂点（ボクセル）にボクセル値（画素値）が含まれる（１）か、含まれない（０）か、で物体の表面形状を確認する。含まれる場合のボクセル値を１（□）、含まれない場合のボクセル値を０（なし）とすると、頂点の１と０の組み合わせは、２^８＝２５６種類できるが、回転や対象性を考慮すると、図１７の２３種類となる。この方法により三次元物体の境界表面に対し、パッチ面（本件の場合、三角形）を推定することができる。
これを、本件に適用する場合（図９−ステップ４）、を以下に示す。
前提条件として皮質骨（図１３）と海面骨（図１４）のボクセル値の取り扱いを以下のように定義する。
定義式４：皮質骨部→１
定義式５：海面骨部→２
定義式６：背景部→０
定義式４は図１３において皮質骨領域のボクセル値を１、定義式５は図１４において海面骨領域のボクセル値を２、定義式６は図１３および図１４の背景領域のボクセル値を０と定義する。
ここで、皮質骨と海面骨の間の境界面を推定することとは、既説明したキューブ上のボクセル値が次の必要条件をすべてみたしていることになる。
定義式７：１と２の両方のボクセル値で構成されている。
定義式８：背景値（０）は存在しない。
定義式７は、着目するキューブ上の頂点に、皮質骨あるいは海面骨の領域双方が存在することを示し、境界面が必ず存在することを意味している。（０）は、着目したキューブ上の頂点に、背景が存在しないことを意味している。
定義式７，８の条件を両方とも満たしている組み合わせを考え、回転や対象性を考慮すると、皮質骨と海面骨の境界キューブのパターンは２３−２＝２１種類に集約される。仮に皮質骨部を（□）、海面骨部（○）とすると、図１８の様なパターンとなる。この方法を適用して、境界面形状を作成し可視化すると図１９の様な画像を得ることが可能である。
また、ラットの大腿骨での適用例を示す。この場合、ラットの大腿骨のみでなく、一般的に長管骨へ適応も可能であることを示すのもである。ラットの大腿骨を所定の方法で分離し３次元表示したのもを図２０、図２１へ示す。また大腿骨の連結面形状を図２２へ示す。
こうして得られた境界面を用いて、図９のステップ５、ステップ６、ステップ７、ステップ８の順で、境界面のラベル数、面積、周囲長、円形度の骨計測を行う。
以下、定義された皮質骨と海面骨の境界面形状に対して、骨計測を行う方法と、手順を示す。図２３に、境界面として定義した三角形の、要素３０、節点３１、辺３２を示す。
まず、境界面のラベル数とは、空間上に独立して存在する面の個数のことである。具体的には、図２４に示す物体が存在しているとすると、左の物体と右の物体の２個の面が存在している。該説明した境界面生成法はあくまでも単位キューブに対して三角形状の要素を定義するものであり、広域としての面の数を画像処理装置１０で認識しているわけではない。骨面計測において、この面の個数を把握しておくことは、境界面の部位に対して定量化する際に基本の単位となり得る。
境界面の面の数（以下、ラベル数と呼ぶ）を計測するには、要素状に存在するある節点の座標に着目し、それと同じ座標を持つ節点が存在するかどうかを全節点に対し確認していく、該処理が終われば次の節点に着目し同処理をすべての節点に対し繰り返す。具体的には、図２５の流れによって境界面のラベル位置とラベル数が決定される。まず、図２５のステップ１にてラベル番号の設定処理を行う。ある要素に着目しラベル番号を付ける。その際、ラベル番号が定義されていなければラベル番号をＬ（Ｌ＝Ｌ＋１）とし、既に定義されているのであれば、その定義されているラベル番号をＬとする。次に図２５のステップ２にて着目した要素上に存在する節点の座標３組６点（Ｘ，Ｙ，Ｚ）を取得する。次に、図２５のステップ３にて着目要素以外の要素上の座標点３点（Ｘ，Ｙ，Ｚ）を取得し重複した節点が存在するかを全節点に対して確認する。図２５のステップ３にて重複点が存在する場合には、ラベル番号Ｌをその要素上に設定する。既に定義されている場合は、その祭儀されているラベル番号を持つ要素のラベルを、すべてＬに置き換える。この図２５のステップ４の処理を、すべての要素に対して同処理を施す。図２５のステップ４の処理が終了したら、図２５のステップ１へ戻り次の別の要素に対して、図２５のステップ１からステップ４の処理を施す。この方法で処理がすべて完了したら、図２５のステップ６にて、定義したラベル番号の昇順ソートを行う。結果的にソート後のラベル番号の最大値が、目的のラベル数（境界面の個数）となる。
この処理において要素数が非常に多い場合、ステップ３の処理が非常に時間がかかる処理となるケースがある。その場合には、全節点に対して重複した節点を確認するのではなく、着目した要素が存在するキューブと同じスライス上に存在する要素と、その上下の３スライス分のキューブ上に存在する要素を比較することによって解決される。
次に境界面の面積を計測するには、式９，式10，式11の計算を行うことによって生成された各々の三角形要素の面積を計算し、加算することによって求められる。
【００１２】
【数４】
式９：ｓ[i]＝（ａ[i]＋ｂ[i]＋ｃ[i]）／２

式10：Ａ[i]＝（ｓ[i]（ｓ[i]−ａ[i]）（ｓ[i]−ｂ[i]）（ｓ[i]−ｃ[i]））^１／２
_ｎ
式11：ＴＡ＝ΣＡ[i]
ⁱ ^＝ ¹
但し、
ｉ：要素番号（１〜最大要素番号）
ａ[i]：要素を構成する辺の長さ１
ｂ[i]：要素を構成する辺の長さ２
ｃ[i]：要素を構成する辺の長さ３
ｓ[i]：辺の長さの総和の１／２
Ａ[i]：各要素の面積
ＴＡ：総面積
【００１３】
式11にて、該境界面に定義したラベル番号ごとに積算を施すことにより、境界面ごとの面積を計測することも可能である。エレメントの総面積であるＴＡを求めることにより、境界面の総面積が計測可能である。
次に境界面の周囲長Ｓを計測するには、要素状に存在するある２つの節点で構成される辺座標に着目し、それと同じ座標を持つ辺が存在するかどうかを全辺に対し確認していく、該処理が終われば次の辺に着目し同処理をすべての辺に対し繰り返す。具体的には、具体的には、図２６の流れによって境界面の周囲長が求められる。まず、図２６のステップ１にてラベル番号Ｌが定義されている要素上の辺座標に着目する。辺座標とは、辺を構成する２点の節点座標（Ｘ，Ｙ，Ｚ）のことである。次に、図２６のステップ２で、ラベル番号Ｌが定義されている着目辺以外の要素の辺と重複がないかを確認する。この処理をステップ３にて、残りの要素に対しても同処理を施す。図２６のステップにて重複がない場合、その辺はラベル番号Ｌを構成する最外側の外周となりそれを、図２６のステップ４にてラベル番号別に積算することにより周囲長Ｓが求められる。最終的には、図２６のステップ５にて、Ｌが最大ラベル番号を超えるまで図２６のステップ１からステップ４の処理を施せば、ラベル番号ごとの周囲長Ｓ（Ｌ）を求めることができる。
次に境界面の円形度φを計測するには、求めた境界の面積Ａと周囲長Ｓから、境界面の境界面ラベル番号毎に次式12より表される。
【００１４】
【数５】
式12： φ［ｉ］：４πＡ［ｉ］／Ｓ［ｉ］^２
ここで、
ｉ：ラベル番号（１〜境界面ラベル数）
φ［ｉ］：境界面の円形度
π［ｉ］：円周率
Ａ［ｉ］：境界面の面積
Ｓ［ｉ］：境界面の周囲長

これは境界面が真円のときは１．０となり、複雑な形状になるほど小さい値（０に近い値）となる。計算の方法としては境界面ラベル毎に、境界面の面積Ａと境界面の周囲長Ｓを代入することにより円形度を計算し、その平均を平均円形度として計測する。その際、境界面の面積Ａと境界面の周囲長Ｓは該説明の境界面の面積Ａと周囲長Ｓを使用する。
表１にラットの腰椎（１０対）、表２にラットの大腿骨の連結形状に管するパラメータであるラベル数、総面積、総周囲長、平均円形度、総面積／ラベル数を計測した結果を示す。総面積／ラベル数は、総面積をラベル数で割った値である。すなわち、独立した単位境界面あたりの平均面積を示している。これらの計測値と骨強度指標との相関を調べれば骨強度の定量的な評価を行うことができる。
【００１５】
【表１】

【表２】

該骨強度指標は、図２７に示しているように、変位量と負荷荷重の関係を調べたときに、その最大荷重を用いて評価した。表３にラットの腰椎（１０対）の最大荷重を示す。
【表３】

この最大荷重と表１で示した境界特性パラメータとの相関を図２８から図３１に示す。結果から、周囲長および平均円形度との相関はあまり確認できなかったが、総面積と総面積／ラベル数との相関が、それぞれ約０．７６、約０．８５と高い相関を示すことがわかった。
【００１６】
【発明の効果】
本発明による手法では、簡単な処理で、安定した皮質骨部及び海綿骨部の分離が可能で、さらに非破壊的に骨の内部構造を構築することが可能であり、さらに高速かつ自動的かつ再現性のある３次元の皮質骨と海綿骨の骨関連情報を生成することができた。特に骨の３次元的解析に対する評価がより正確に行えることが期待できる。また、物体の表面形状を推定する処理の応用で、簡単に安定した皮質骨部と海綿骨部の間の境界面を推定可能で、さらに非破壊的に骨の内部構造を可視化することが可能であり、さらに再現性のある３次元の皮質骨と海綿骨の境界面特性を定量化することができた。特に骨の３次元的解析に対する評価がより正確に行えることが期待できる。これは、腰椎のみでなく大腿骨いわゆる長管骨への適用についても可能である。本発明は、μＸ線ＣＴ撮影装置をはじめＸ線フィルムまたはＭＲＩ等から得られる診断撮影写真を骨の３次元的評価の基礎となる画像情報全般に応用が期待できる。さらに、本発明は、骨関連疾患あるいは関節関連疾患動物モデルにおいて骨および関節を評価する評価方法や骨関連疾患あるいは関節関連疾患患者の骨および関節を評価する臨床評価方法として応用が可能であり、動物モデルを用いた薬剤や治療法の開発、および臨床試験、日常診療や正確な診断、薬剤や治療方法の効果判定への応用が期待できる。
【００１７】
【産業上の利用可能性】
本発明による手法は、コンピュータ用の種々のオペレーションシステム上で作動するソフトウエアとして利用することが可能であり、あらゆるコンピュータ言語で作成することができ、いかなる装置に対してもソフトウエアとして組み込むことができる。
【００１８】
【図面の簡単な説明】
図１は、本発明を実施する装置例の構成の説明図である。
図２は、皮質骨と海綿骨を分離する分離処理手段のフローチャートである。
図３は、ラット腰椎断面の原画像である。
図４は、ラット腰椎断面の２値化画像である。
図５は、皮質骨と海綿骨の分離方法の説明図である。
図６は、皮質骨と海綿骨の分離後の濃淡図である。
図７は、３次元の画像の生成を説明する説明図である。
図８は、生成された皮質骨と海綿骨の三次元画像を示す図である。
図９は、境界面特性を計算するための全体フローチャートである。
図１０は、ラット腰椎断面の原画像である。
図１１は、ラット腰椎断面の皮質骨の２値化画像である。
図１２は、ラット腰椎断面の海面骨の２値化画像である。
図１３は、ラット腰椎の皮質骨の３次元化可視化画像である。
図１４は、ラット腰椎の海面骨の３次元化可視化画像である。
図１５は、ボクセル定義を説明する図である。
図１６は、キューブ定義を説明する図である。
図１７は、マーチングキューブス法のパターンを示す図である。
図１８は、本件の境界面のパターンを示す図である。
図１９は、ラット腰椎の皮質骨と海綿骨の間の境界面を示す図である。
図２０は、ラット大腿骨の皮質骨の３次元化可視化画像である。
図２１は、ラット大腿骨の海面骨の３次元化可視化画像である。
図２２は、ラット大腿骨の皮質骨と海綿骨の間の境界面を示す図である。
図２３は、三角形要素の用語を説明する図である。
図２４は、ラベリング定義を示す図である。
図２５は、境界面の個数を求めるフローチャートである。
図２６は、周囲長計算のフローチャートである。
図２７は、骨強度を反映する指標を説明する図である。
図２８は、最大荷重−総面積の相関図である。
図２９は、最大荷重−総周囲長の相関図である。
図３０は、最大荷重−平均円形度の相関図である。
図３１は、最大荷重−総面積／ラベル数の相関図である。
【００１９】
【符号の説明】
図中の符号は、それぞれ以下のことを表わす。
１０画像処理装置
１１パーソナルコンピュータ
１２画像処理ボード
１３ハードディスク
２０骨部
２１皮質骨部
２２海綿骨部
２３骨髄腔部（空間部）
２４骨髄腔部（空洞部）
２５背景部
３０要素
３１節点
３２辺

Claims

被検骨断面の２値化画像に基づいて皮質骨と海綿骨とを分離した画像上にて、その分離切断面部（以下、境界面とよぶ）に対し物体の境界面を定義する方法を適用し、皮質骨と海面骨の境界面を定義する画像処理方法であり、３次元の骨分離画像を、被検骨を所定間隔で断層撮影した複数の連続した各画像から皮質骨部と海綿骨部とが分離された各部の分離画像を作成し、該分離画像を積み重ねることにより生成することを特徴とする、画像処理方法にて生成される境界面に対し、境界面成分のラベリングを行い、３次元空間上に独立している境界面の個数に関する情報を得るための骨計測方法。
被検骨断面の２値化画像に基づいて皮質骨と海綿骨とを分離した画像上にて、その分離切断面部（以下、境界面とよぶ）に対し物体の境界面を定義する方法を適用し、皮質骨と海面骨の境界面を定義する画像処理方法であり、３次元の骨分離画像を、被検骨を所定間隔で断層撮影した複数の連続した各画像から皮質骨部と海綿骨部とが分離された各部の分離画像を作成し、該分離画像を積み重ねることにより生成することを特徴とする、画像処理方法にて生成される境界面に対し、３次元空間上に独立している境界面の面積に関する情報を得るための計測方法。
被検骨断面の２値化画像に基づいて皮質骨と海綿骨とを分離した画像上にて、その分離切断面部（以下、境界面とよぶ）に対し物体の境界面を定義する方法を適用し、皮質骨と海面骨の境界面を定義する画像処理方法であり、３次元の骨分離画像を、被検骨を所定間隔で断層撮影した複数の連続した各画像から皮質骨部と海綿骨部とが分離された各部の分離画像を作成し、該分離画像を積み重ねることにより生成することを特徴とする、画像処理方法にて生成される境界面に対し、３次元空間上に独立している境界面の周囲長に関する情報を得るための計測方法。
被検骨断面の２値化画像に基づいて皮質骨と海綿骨とを分離した画像上にて、その分離切断面部（以下、境界面とよぶ）に対し物体の境界面を定義する方法を適用し、皮質骨と海面骨の境界面を定義する画像処理方法であり、３次元の骨分離画像を、被検骨を所定間隔で断層撮影した複数の連続した各画像から皮質骨部と海綿骨部とが分離された各部の分離画像を作成し、該分離画像を積み重ねることにより生成することを特徴とする、画像処理方法にて生成される境界面に対し、３次元空間上に独立している境界面の円形度に関する情報を得るための計測方法。