(57)【特許請求の範囲】
【請求項１】
人又は動物の骨梁の撮影により得られ、骨部分が一定の大きさのボクセルから構成される所定の三次元の範囲の骨梁デジタルイメージモデルをデータの形で記憶装置に保存し骨折のリスクを評価するためにシミュレーションによって処理する骨折リスク評価のためのコンピュータの作動方法であって、
前記骨梁デジタルイメージモデルにストレス負荷をかけながら、骨梁表面のボクセルにおける応力が、周囲の平均応力よりも生成閾値レベル以上に大きければ当該ボクセルに接する新たな１以上のボクセルを生成し、周囲の平均応力よりも消失閾値レベル以下に小さければ当該ボクセルを含む１以上のボクセルを消失させる過程を繰り返して前記骨梁デジタルイメージモデルを再構築する前記コンピュータによる処理工程と、
前記再構築後の骨梁デジタルイメージモデルに少なくとも前記ストレス負荷と直交方向に測定負荷をかけて応力分布測定を行う前記コンピュータによる処理工程と、を有することを特徴とする骨折リスク評価のためのコンピュータの作動方法。
【請求項２】
請求項１に記載の骨折リスク評価のためのコンピュータの作動方法において、
患者の測定された骨密度に合わせて生成閾値レベルと消失閾値レベルを調整し、前記再構築により前記骨梁デジタルイメージモデルを作製し、当該骨密度に対応する応力分布を測定することを特徴とする骨折リスク評価のためのコンピュータの作動方法。
【請求項３】
請求項１又は２に記載の骨折リスク評価のためのコンピュータの作動方法において、
前記骨梁デジタルイメージモデルのボクセルにおけるＢＡｐ結晶の配向方向のヤング率を他の方向よりも大きくして応力を算出することを特徴とする骨折リスク評価のためのコンピュータの作動方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、骨粗鬆症の診断又は治療のために好適な骨折リスク評価のためのコンピュータの作動方法に関する。
【背景技術】
【０００２】
骨粗鬆症は、骨梁の骨密度の減少と微細構造の変化によって骨折リスクが高まった病態である。特に、高齢になってからの骨折は、寝たきりなどを誘発させ易くなる。骨粗鬆症の国内の患者は、１０００万人以上とも言われている。骨梁は、各微細部分が骨芽細胞により作られ（骨形成し）たり破骨細胞により壊され（骨吸収し）たりするという再構築（リモデリング）を長期的に繰り返しており、骨梁の微細構造は骨形成と骨吸収の量のバランスにより変化する。骨粗鬆症は、骨形成と骨吸収の量のバランスが崩れ、骨形成よりも骨吸収の量が相対的に多くなったものである。骨粗鬆症の診断又は治療のためには、現在、Ｘ線などによる骨密度の測定が行われているのであるが、より高度な診断又は治療のためには、骨粗鬆症の骨梁の力学的評価が必要である。
【０００３】
ところで、骨梁の力学的評価のために、非特許文献１に示すように、骨梁の微細構造を再現した骨梁デジタルイメージモデルに圧縮負荷をかけ、有限要素法解析（ＦＥＡ）によって骨梁表面の応力の不均一性を算出し、その不均一性に応じた再構築をシミュレートする力学的解析手法が提案されている。なお、非特許文献１に関連して、特許文献１には、人工的な骨代替物の構造を設計する手法が記載されている。
【０００４】
そして、これまでに、本願発明者らは、例えば非特許文献２に示すように、この力学的解析手法に多方面から改良を加えてきた。非特許文献２は、骨質の重要な因子であるＢＡｐ結晶の配向に起因する骨梁の材料異方性を考慮した力学的解析手法による骨折リスク評価の可能性を報告したものである。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００２−０６５７０８号公報
【非特許文献】
【０００６】
【非特許文献１】安達泰治、外２名、“デジタルイメージモデルを用いた海綿骨の力学的再構築シミュレーション”、日本機械学会論文集（Ａ編）、日本機械学会、２０００年、第６６巻、第６４８号、ｐ．２００−２０７
【非特許文献２】田原大輔、外３名、“骨質を考慮した椎体海綿骨の高分解能ミクロ応力解析”、日本臨床バイオメカニクス学会誌、日本臨床バイオメカニクス学会、２００８年、第２９巻、ｐ．７−１４
【発明の概要】
【発明が解決しようとする課題】
【０００７】
非特許文献１、２を含め、骨梁の力学的解析手法は、日常的に自重による圧縮負荷がかかる体軸方向に着目するのが基本である。しかし、骨粗鬆症の患者にあっては、転倒時などの非日常的な負荷がかかることによる骨折が現実的には多い。
【０００８】
本発明は係る事由に鑑みてなされたものであり、その目的は、より現実的な負荷を考慮した力学的解析手法による骨折リスク評価のためのコンピュータの作動方法を提供することにある。
【課題を解決するための手段】
【０００９】
上記目的を達成するために、請求項１に記載の骨折リスク評価のためのコンピュータの作動方法は、人又は動物の骨梁の撮影により得られ、骨部分が一定の大きさのボクセルから構成される所定の三次元の範囲の骨梁デジタルイメージモデルをデータの形で記憶装置に保存し骨折のリスクを評価するためにシミュレーションによって処理する骨折リスク評価のためのコンピュータの作動方法であって、前記骨梁デジタルイメージモデルにストレス負荷をかけながら、骨梁表面のボクセルにおける応力が、周囲の平均応力よりも生成閾値レベル以上に大きければ当該ボクセルに接する新たな１以上のボクセルを生成し、周囲の平均応力よりも消失閾値レベル以下に小さければ当該ボクセルを含む１以上のボクセルを消失させる過程を繰り返して前記骨梁デジタルイメージモデルを再構築する前記コンピュータによる処理工程と、前記再構築後の骨梁デジタルイメージモデルに少なくとも前記ストレス負荷と直交方向に測定負荷をかけて応力分布測定を行う前記コンピュータによる処理工程と、を有することを特徴とする。
【００１０】
請求項２に記載の骨折リスク評価のためのコンピュータの作動方法は、請求項１に記載の骨折リスク評価のためのコンピュータの作動方法において、患者の測定された骨密度に合わせて生成閾値レベルと消失閾値レベルを調整し、前記再構築により前記骨梁デジタルイメージモデルを作製し、当該骨密度に対応する応力分布を測定することを特徴とする。
【００１１】
請求項３に記載の骨折リスク評価のためのコンピュータの作動方法は、請求項１又は２に記載の骨折リスク評価のためのコンピュータの作動方法において、前記骨梁デジタルイメージモデルのボクセルにおけるＢＡｐ（Ｂｉｏｌｏｇｉｃａｌ  Ａｐａｔｉｔｅ  （生体アパタイト））結晶の配向方向のヤング率を他の方向よりも大きくして応力を算出することを特徴とする。
【発明の効果】
【００１２】
本発明によれば、骨梁デジタルイメージモデルを再構築し、ストレス負荷と直交方向に測定負荷をかけることにより、より現実的な負荷を考慮した力学的解析手法による骨折リスク評価のためのコンピュータの作動方法を提供することができる。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施形態に係る骨折リスク評価のためのコンピュータの作動方法に用いる骨梁デジタルイメージモデルをストレス負荷とともに例示する立体図である。
【図２】同上の骨梁デジタルイメージモデルの再構築を示すものであって、（ａ）がボクセルの生成・消失の模式図、（ｂ）が再構築の簡略フロー図である。
【図３】同上の骨梁デジタルイメージモデルの再構築を説明するものであって、（ａ）が骨梁表面のボクセルｘｃ近傍の模式図、（ｂ）が重み関数、（ｃ）が骨梁表面変化速度ｍの関数曲線である。
【図４】同上の骨梁デジタルイメージモデルを測定負荷とともに例示する立体図である。
【図５】同上の骨梁デジタルイメージモデルの再構築による骨密度の変化を示す特性グラフである。
【図６】同上の骨梁デジタルイメージモデルに測定負荷をかけたときの応力分布を示す特性グラフである。
【図７】同上の骨梁デジタルイメージモデルにストレス負荷をかけたときの主応力を模式的に例示する断面部分図である。
【発明を実施するための形態】
【００１４】
以下、本発明の実施形態を、図面を参照しながら説明する。本発明の実施形態に係る骨折リスク評価のためのコンピュータの作動方法は、記憶装置にデータの形で保存されてコンピュータによって処理される、骨梁の微細構造を再現した骨梁デジタルイメージモデル１を用いるものである（図１参照）。この骨梁デジタルイメージモデル１は、人の（場合によっては、動物の）骨梁の撮影により得られたデジタルイメージが初期形状時のものであり、骨部分が一定の大きさのボクセル（ｖｏｘｅｌ）の集合体から構成される所定の三次元の範囲のものである。そして、この骨折リスク評価のためのコンピュータの作動方法は、初期形状時の骨梁デジタルイメージモデル１にストレス負荷をかけながら骨梁デジタルイメージモデル１を再構築し、再構築後の骨梁デジタルイメージモデル１に測定負荷をかけて応力分布測定を行うシミュレーションによって骨折のリスクを評価するものである。
【００１５】
骨梁デジタルイメージモデル１を得るための骨梁の撮影は、高分解能（例えば、６４μｍ）のＸ線μＣＴシステムなどを用いることができる。撮影された画像は、画像処理されてから初期形状時の骨梁デジタルイメージモデル１となる。
【００１６】
骨梁デジタルイメージモデル１におけるボクセルの数は、骨密度に対応するものである。骨梁デジタルイメージモデル１は、例えば、４．７ｍｍ立方の範囲に切り取られたものである。また、ボクセルのサイズは、例えば、撮影された画像の分解能と同じ１辺６４μｍである。骨粗鬆症は、日本骨代謝学会の骨粗鬆症基準に基づくと、骨密度が健常例の７０％以下となっている。骨粗鬆症でない健常例の骨梁デジタルイメージモデル１Ａは、例えば、後述するように１８５，２５５個のボクセルから構成される。これに対し、骨粗鬆症例の骨梁デジタルイメージモデル１Ｂは、ボクセルの数が少なく、例えば、後述するように１２５，９４５個（骨梁デジタルイメージモデル１Ａの６８．０％）のボクセルから構成される。
【００１７】
骨梁デジタルイメージモデル１は、ヤング率（例えば、２０ＧＰａ）、ポアソン比（例えば、０．３０）を設定し、負荷がかけられると、有限要素法解析（ＦＥＡ）により骨梁表面のボクセルにおける応力が算出される。
【００１８】
ストレス負荷は、通常は、体軸方向の圧縮負荷であり、例えば、図１に示すように、骨梁デジタルイメージモデル１の下面を固定し、上面に１０ＭＰａをかけることによって、体軸方向（ｚ軸）の圧縮負荷とする。ここで、図１や後述する図４において、ｘ軸、ｙ軸、ｚ軸は、人の左右方向、人の前後方向、人の体軸方向に、それぞれ対応するものである。
【００１９】
骨梁デジタルイメージモデル１の再構築は、図２（ａ）に示すように、骨梁表面のボクセルにおける応力が、周囲の平均応力よりも生成閾値レベル以上に大きければ当該ボクセルに接する新たな１以上のボクセルを生成し、周囲の平均応力よりも消失閾値レベル以下に小さければ当該ボクセルを含む１以上のボクセルを消失させる。このボクセルの生成と消失は、実際の骨における骨形成と骨吸収に対応する。そして、骨梁デジタルイメージモデル１の再構築は、この過程を、図２（ｂ）に示すように繰り返す（例えば、３０〜５０ステップ）ことによって行う。この再構築により、骨梁デジタルイメージモデル１の骨梁の形態が変化する。なお、１ステップを、実験又は経験に基づいた実際の骨における骨形成と骨吸収の速度に対応させることも可能である。
【００２０】
より詳細には、以下に述べる再構築則を用いるのが好ましい。すなわち、図３（ａ）に示すように、骨梁デジタルイメージモデル１の骨梁表面のボクセルｘｃの応力値をσｃ、ボクセルｘｃ周囲の平均応力をσｄとし、応力不均一性Γを（１）式のように定義する。
【００２１】
【数１】
0005829921.tif 000002
【００２２】
ここで、σｄは、ボクセルｘｃから半径Ｒ（例えば、１０ボクセル）の範囲に含まれる骨梁表面上のボクセルｘｒにおける応力値σrの重み付き平均値により、（２）式のように表すことができる。Ｓは骨梁表面、ｗ（ｒ）はボクセルｘｃとｘｒの距離ｒに関する重み関数（図３（ｂ）参照）、をそれぞれ示している。
【００２３】
【数２】
0005829921.tif 000003
【００２４】
この応力不均一性Γは、ボクセルｘｃにおいて応力の一様化に向かう再構築駆動力となる。そして、骨梁表面変化速度ｍを、図３（ｃ）に示すような応力不均一性Γの連続関数として定義する。骨梁表面変化速度ｍは、応力不均一性Γが上記の生成閾値レベルに対応する閾値Γｕ以上のときに正の値となって、値の大きさに応じてボクセルｘｃに接する新たな１以上のボクセルが生成される。また、骨梁表面変化速度ｍは、応力不均一性Γが上記の消失閾値レベルに対応する閾値Γｌ以下のときに負の値となって、値の大きさに応じてボクセルｘｃを含む１以上のボクセルを消失させる。よって、閾値Γｕと閾値Γｌの間は、再構築の不感帯である。閾値Γｕと閾値Γｌは、骨梁表面の見かけの変化速度を調整する値である。
【００２５】
この再構築則によるボクセルの生成又は消失の過程を、所定のステップ数、繰り返す。
【００２６】
次に、再構築後の骨梁デジタルイメージモデル１の応力分布測定は、少なくともストレス負荷と直交方向に（ストレス負荷をかける方向に直交する方向に）測定負荷をかけて行う。例えば、図４に示すように、骨梁デジタルイメージモデル１の下面を固定し、上面に測定負荷を体軸方向（ｚ軸）に直交する方向にかけるのである。
【００２７】
以上説明した骨折リスク評価のためのコンピュータの作動方法の使用例として、再構築によって骨密度に応じた骨梁デジタルイメージモデル１を作製し、その骨密度に対応する応力分布を測定し、骨折のリスク評価を行う方法を説明する。
【００２８】
健常例の骨梁デジタルイメージモデル１Ａと、健常例に対し骨密度が約７０％の骨粗鬆症例の骨梁デジタルイメージモデル１Ｂは、健常な骨梁を６４μｍの分解能のＸ線μＣＴシステムを用いて撮影し、その画像処理によって４．７ｍｍ立方の範囲のデジタルイメージを初期形状時のものとした。健常例の骨梁デジタルイメージモデル１Ａは、閾値Γｕを４．０、閾値Γｌを−５．０として上記再構築則に従って作製した。骨粗鬆症例の骨梁デジタルイメージモデル１Ｂは、閾値Γｕを４．０、閾値Γｌを−０．５として上記再構築則に従って作製した。健常例の骨梁デジタルイメージモデル１Ａと骨粗鬆症例の骨梁デジタルイメージモデル１Ｂのおける閾値Γｕと閾値Γｌは、再構築後のものがほぼ目標の骨密度になるように調整したものである。また、ストレス負荷は、健常例の骨梁デジタルイメージモデル１Ａと骨粗鬆症例の骨梁デジタルイメージモデル１Ｂの上下に剛体板を付加し、上面に１０ＭＰａをかけることによって、体軸方向（ｚ軸）の圧縮負荷とした。
【００２９】
そうすると、図５に示すように、再構築の進行に応じて、健常例の骨梁デジタルイメージモデル１Ａはボクセルの数（すなわち、骨密度）が変化せず（１８５，２５５個）、骨粗鬆症例の骨梁デジタルイメージモデル１Ｂはボクセルの数が急激に減少している。骨粗鬆症例の骨梁デジタルイメージモデル１Ｂは、約５０ステップ程度で各ボクセルにおける応力がほぼ一様になって、骨密度が約７０％の状態でほぼ安定する。より詳細には、骨粗鬆症例の骨梁デジタルイメージモデル１Ｂは、骨密度が６８．０％で、ボクセルの数が１２５，９４５個のものが得られた。
【００３０】
図６に示すのは、１０ＭＰａの測定負荷をかけたときに、健常例の骨梁デジタルイメージモデル１Ａと骨粗鬆症例の骨梁デジタルイメージモデル１Ｂの各々のボクセルに発生する応力の大きさの存在割合（応力分布）を高応力側からの累積曲線で示したものである。図６において、曲線ａは骨粗鬆症例の骨梁デジタルイメージモデル１Ｂに測定負荷をストレス負荷と直交方向にかけたもの、曲線ｂは健常例の骨梁デジタルイメージモデル１Ａに測定負荷をストレス負荷と直交方向にかけたもの、曲線ｃは骨粗鬆症の骨梁デジタルイメージモデル１Ｂに測定負荷をストレス負荷と同じ方向にかけたもの、曲線ｄは健常例の骨梁デジタルイメージモデル１Ａに測定負荷をストレス負荷と同じ方向にかけたものである。
【００３１】
図６においては、高応力が発生するボクセルの数の割合は、曲線ａが最も大きく、７５ＭＰａ以上になるボクセルは約７０％も存在する。これに対し、曲線ｂ、ｃ、ｄではそれぞれ、応力が７５ＭＰａ以上になるボクセルは、約３５％、約１２％、約５％存在する。
【００３２】
このように、再構築によって骨密度に応じた骨梁デジタルイメージモデル１（健常例の骨梁デジタルイメージモデル１Ａ、骨粗鬆症例の骨梁デジタルイメージモデル１Ｂ）を得ることがでる。そして、骨粗鬆症例の骨梁デジタルイメージモデル１Ｂは、ストレス負荷と直交方向に測定負荷をかけると、高応力が発生するボクセルの数の割合が極めて多く、このことは、転倒時などの非日常的に負荷がかかる場合の骨折が多いことをよく表しており、しかも、骨折のリスクを理解し易い数値で示すことができている。このことは、健常例の骨梁デジタルイメージモデル１Ａと比較するとより明確になっている。こうして、骨密度に応じて、より現実的な負荷を考慮した力学的解析手法による骨折のリスク評価が可能になる。
【００３３】
更には、患者の測定された骨密度に合わせて、上述した骨粗鬆症の骨梁デジタルイメージモデル１Ｂのように、生成閾値レベル（閾値Γｕ）と消失閾値レベル（閾値Γｌ）を調整して再構築により骨梁デジタルイメージモデル１を作製し、その骨密度に対応する応力分布を測定し、骨折のリスク評価を行うことも可能である。
【００３４】
また、初期形状時の骨梁デジタルイメージモデル１として選択できるサンプルモデルが複数用意できる場合は、患者の性別や年齢などに応じて初期形状時の骨梁デジタルイメージモデル１を選ぶことも可能である。また、患者自身の骨梁の撮影ができて、生成閾値レベル（閾値Γｕ）、消失閾値レベル（閾値Γｌ）及びステップ数を実験又は経験に基づき、実際の骨における骨形成と骨吸収に対応して精度良く設定できる場合は、その患者の骨粗鬆症の進行を予測した骨折のリスク評価も可能である。
【００３５】
次に、骨梁デジタルイメージモデル１の各々のボクセルにおけるＢＡｐ（Ｂｉｏｌｏｇｉｃａｌ  Ａｐａｔｉｔｅ  （生体アパタイト））結晶の配向方向のヤング率を他の方向よりも大きくして応力を算出することを説明する。このような材料異方性を考慮した骨梁デジタルイメージモデル１を用いると、より現実の骨梁に近い力学的解析が可能になる。ＢＡｐ結晶配向は、図７に示すように、ＢＡｐ結晶配向を考慮していない等方性の骨梁デジタルイメージモデル１に体軸方向（ｚ軸方向）に一定のストレス負荷（圧縮負荷）をかけることによって各々のボクセルに応力を発生させ、主応力の方向（最大となる方向）（図７では、黒矢印で示す）に一致させると、適切かつ比較的容易に決定することができる。
【００３６】
健常例と骨粗鬆症例の骨梁デジタルイメージモデル１Ａ、１Ｂに対し、各方向に測定負荷をかけてマクロな剛性をＢＡｐ結晶配向の考慮有無で比較した結果を表１に示す。ここで、ＢＡｐ結晶の配向方向のヤング率は、その直交方向の２倍とした。マクロな剛性は、均質化法により得たヤング率としている。ＢＡｐ結晶配向を考慮したものは、体軸方向（ｚ軸方向）の剛性が高くなっているが、それに直交する方向（ｘ軸方向、ｙ軸方向）の剛性はさほど高くはなっていないことが分かる。よって、ＢＡｐ結晶配向を考慮した現実の骨梁により近い骨梁デジタルイメージモデル１は、転倒時などの非日常的に負荷がかかる場合の骨折が日常的な自重による負荷がかかる場合よりも多いことを更によく表している。
【００３７】
【表１】
0005829921.tif 000004
【００３８】
以上、本発明の実施形態に係る骨折リスク評価のためのコンピュータの作動方法について説明したが、本発明は、実施形態に記載したものに限られることなく、特許請求の範囲に記載した事項の範囲内でのさまざまな設計変更が可能である。例えば、再構築後の骨梁デジタルイメージモデルの応力分布測定は、前記ストレス負荷と直交方向に測定負荷をかけて行うこととあわせて、前記ストレス負荷と同じ方向に測定負荷をかけて行ってもよい。
【符号の説明】
【００３９】
１    骨梁デジタルイメージモデル
【図５】
0005829921.tif 000005
【図６】
0005829921.tif 000006
【図１】
0005829921.tif 000007
【図２】
0005829921.tif 000008
【図３】
0005829921.tif 000009
【図４】
0005829921.tif 000010
【図７】
0005829921.tif 000011