JPH0331061B2 - - Google Patents
Info
- Publication number
- JPH0331061B2 JPH0331061B2 JP59230077A JP23007784A JPH0331061B2 JP H0331061 B2 JPH0331061 B2 JP H0331061B2 JP 59230077 A JP59230077 A JP 59230077A JP 23007784 A JP23007784 A JP 23007784A JP H0331061 B2 JPH0331061 B2 JP H0331061B2
- Authority
- JP
- Japan
- Prior art keywords
- bone
- pattern
- ray
- image
- density
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/505—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of bone
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Description
【発明の詳細な説明】
〈技術分野〉
本発明は骨のＸ線写真像の評価方法に関する。
更に詳細には本発明は、四肢長管骨のＸ線写真陰
影度を測定して骨パターンを求め、得られる骨パ
ターンを修正し、修正された骨パターンに基い
て、各種の評価を行なう骨のＸ線写真像の評価方
法に関する。 〈従来技術〉 人間の骨の発育状態，老化度、あるいは骨粗鬆
症，骨軟化症等の骨病変の種類などを評価する方
法としてMD法が知られている（骨代謝，第13
巻，187〜195頁，1980；骨代謝，第14巻，91〜
104頁，1981等参照）。 MD法は、第２中手骨の骨幹部中央のＸ線像の
陰影度をマイクロデンシトメーターを用いて測定
し、コンピユーターを併用して陰影より、骨皮質
幅指数（MCL），骨幅(D)，最高骨濃度（GSnax），
最低骨濃度（GSnio），平均骨濃度（ΣGS／D.）
等の各指標を求め、この指標によつて、骨萎縮度
等を判定する方法である。 〈発明が解決しようとする問題点〉 MD法にあつては、MD法で得られる各指標、
例えば最高骨濃度（GSnax），骨皮質幅指数
（ACI）等の数値が、実際の骨の状態を反映しな
い場合があり、これらの指標によつては骨の評価
が困難なことも起り得る。たとえば最高骨濃度
（GSnax）は骨皮質の骨密度の指標とされている
が、骨密度が同じ骨の場合でも骨の直径や骨皮質
幅が異なれば違う数値を示し、また骨皮質が均質
でない場合には、その意味することが不明確とな
る。 フオトン・アブソープシオメトリーによる従来
の方法では、骨塩量（BMC）を求求めることし
か知られていない。 〈問題点を解決するための手段〉 本発明者らは、直接的に骨の状態を反映し得る
指標を求めることを目的として鋭意研究した結
果、第２中手骨のＸ線像の陰影度を、例えば従来
のMD法と同様にして、測定して骨パターンを求
め、該骨パターンの平滑化及び対称化を行なつて
骨パターンを修正し、該修正された骨パターンよ
り、各指標を求めることによつて得られる各指標
が、骨の状態をよく反映し得ることを見出した。 本発明者らはかかる知見に基いて更に鋭意研究
した結果、本発明を完成したものである。 本発明は第１に、 (a) 四肢長管骨をＸ線撮影して得られる骨のＸ線
写真像陰影度を測定して骨パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行なつて
骨パターンを修正し、 (c) 該修正された骨パターンより、骨皮質幅指数
（MCI），骨幅(D)，最高骨濃度（GSnax），最低
骨濃度（GSnio），骨塩量（ΣGS）及び平均骨濃
度（ΣGS／ＤあるいはΣGS／D2）より選ばれ
る少なくとも１種の指標を求めて、 (d) 該指標により骨の評価を行なうことを特徴と
する骨のＸ線写真像の評価方法。 第２に、 (a) 四肢長管骨をＸ線撮影して得られる骨のＸ線
写真像陰影度を測定して骨パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行なつて
骨パターンを修正し、 (c) 該修正された骨パターンより、骨の断面の外
形を楕円とし帯状の骨皮質をもちその内側に骨
密度が減少する部分を持つ骨モデルを設定し
て、骨密度分布を求め、 (d) 該骨密度分布により骨を評価することを特徴
とする骨のＸ線写真像の評価方法。 第３に、 (a) 四肢長管骨をＸ線撮影して得られる骨のＸ線
写真像陰影度を一定の間隔で測定して各部位の
骨パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行なつて
各部位の骨パターンを修正し、 (c) 該修正された骨パターンより、骨の断面の外
形を楕円とし帯状の骨皮質をもちその内側に骨
密度が減少する部分を持つ骨モデルを設定し
て、各部位の骨密度分布を求め、 (d) 各部位の骨密度分布を、骨密度分布ごとにそ
の値の大小により色分けして、 (e) 四肢長管骨のＸ線像を、骨密度分布の像に変
換して、 (f) 該像により骨の評価を行なうことを特徴とす
る骨のＸ線写真像の評価方法である。 以下第１の発明について説明する。 先ず、(a)、四肢長管骨をＸ線撮影して得られる
骨のＸ線写真陰影度を測定して骨パターンを求め
る。あるいは四肢長管骨からフオトン・アブソー
プシオメトリーにより骨パターンを求める。 四肢長管骨としては、骨のまわりの軟部組織の
厚さが薄く平均化している部分の骨が望ましく、
例えば第１〜第５の中手骨、上腕骨などが挙げら
れる。なかでも第２中手骨の骨幹部中央が好適で
ある。四肢長管骨をＸ線撮影してＸ線の像を得
る。Ｘ線撮影して得られる骨のＸ線写真陰影度を
測定して骨パターンを求める。骨パターンを求め
るには、例えばMD法によれば、Ｘ線像の陰影度
を、例えば１段の高さ１mmで20段（最低の高さ１
mm，最高の高さ20mm）のアルミ階段（alminium
step wedge）、あるいはアルミニウムロープのＸ
線像とともにデンシトメーターで測定する。骨の
近位端と速位端との中間点の陰影度を測定するの
が望ましい。得られる陰影度は、光学密度
（OD）を目盛つた用紙に拡大して記録し、骨の
各点における吸光度（OD）をアルミ階段に段数
に変換補正（GS値）して第１図に示すような骨
パターンを作図する（骨代謝，第13巻，187〜195
頁，1980；骨代謝，第14巻，91〜104頁，1981；
骨代謝，第14巻，319〜325頁，1981）。あるいは、
骨のＸ線像を、アルミ階段、アルミニウムスロー
プのＸ線像とともにテレビカメラにより読み取つ
て骨パターンを求めることもできる。 次いで(b)、得られる骨パターンの対称化及び平
滑化を行なつて骨パターンを修正する。 骨パターンの対称化を第２図に基いて説明す
る。先ず上述した方法で得られる骨パターンに中
心線（第２図の一点鎖線）を引く。中心線より左
右に等間隔（Ｘ）の位置におけるGS値（GSx1及
びGSx2）を求め、これらの平均値〔GSx＝
（GSx1＋GSx2）／２〕を、中心線より左右に等
間隔の位置のGS値（GSx）とする。Ｘを変えて
この操作をくり返すことにより、対称化された骨
パターン（第２図の点線）を求める。次いで骨パ
ターンの平滑化を行なう。平滑化の方法を第３図
に基いて説明する。先ず上述した方法で得られる
対称化された骨パターンを等分化する。次いで等
分化された骨パターンを、例えば５点ずつ区分す
る（Xi1，Xi2，Xi3，Xi4，Xi5）。次いで各位置
におけるGS値（GSx1，GSxi2，GSxi3，GSxi4，
GSxi5）を求めその平均値を算出して、その値
を、５区分の中心であるXi3におけるGS値とす
る。この操作を順次行ない得られるGS値を結ぶ
ことに平滑化された骨パターンが得られる。この
平滑化は移動平均法に準じた方法である。 本発明では、上記対称化，平滑化の操作の順は
いずれでもよいが、対称化を先ず行ない、次いで
平滑化するのが望ましい。 かくして得られる修正された骨パターンより、
各指標を求める。各指標について例えば第２中手
骨に本法を用いた場合に得られる第４図に基いて
説明する。 第４図において、Ｄは骨幅，ｄは骨髄幅，d1は
尺側骨皮質幅，d2は橈側骨皮質幅を示している。
骨幅(D)，尺側骨皮質幅（d1），橈側骨皮質幅（d2）
より骨皮質幅指数（MCI）が求まれる（MCI＝
（d1＋d2）／Ｄ）。GSnioは最低骨濃度であり、
GSmax1（尺側骨皮質におけるGS値の極大値）と
GSmax2（橈側骨皮質におけるGS値の極大値）と
の平均値（本発明の場合、GSmax1とGSmax2と
は、骨パターンの対称化により等しい）より最高
濃度（GSmax）〔GSmax＝1/2（GSmax1＋
GSmax2）〕が求まる。また各GS値を積分して得
られる面積が骨塩量（ΣGS）であり、これより
平均骨濃度（ΣGS／ＤまたはΣGS／D2）が求ま
る。 かくして得られる指標、骨皮質幅指数
（MCI），骨幅(D)，最高骨濃度（GSnax），最低骨
濃度（GSnio），骨塩量（ΣGS）及び平均骨濃度
（ΣGS／ＤまたはΣGS／D2）より選ばれる少なく
とも１種の指標によつて骨の評価を行なうことが
できる。これらの指標、例えば骨皮質幅指数
（MCI）は、従来のMD法で用いられるMCIに比
べて実際の骨の状態をより適確に表わした値であ
る。したがつてこれらの指標により、より正確な
骨の評価が可能である。骨の評価の例としては、
例えば骨萎縮度の判定（骨代謝，第13巻，187〜
195頁，1980），骨成長の判定（骨代謝，第14巻，
319〜325頁，1981），骨病変の類別判定（特開昭
59−8935，特開昭59−11840号公報），骨改善効果
の判定（特開昭59−49743号公報）などが挙げら
れる。本発明で得られる指標を用いてこれら判定
を行なうことが可能である。 次いで第２の発明、すなわち、 (a) 四肢長管骨をＸ線撮影して得られる骨のＸ線
写真陰影度を測定して骨パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行なつて
骨パターンを修正し、 (c) 該修正された骨パターンより、骨の断面の外
形を楕円とし帯状の骨皮質をもちその内側に骨
密度が減少する部分を持つ骨モデルを設定し
て、密度分布を求め、 (d) 該骨密度分布により骨を評価することを特徴
とする骨の評価方法について説明する。 第２の発明では先ず、(a)四肢長管骨をＸ線撮影
して得られる骨のＸ線写真陰影度を測定して骨パ
ターンを求め、(b)該骨パターンの平滑化及び対称
化を行なつて骨パターンを修正する。これら(a)，
(b)の操作は第１の発明において説明したと同様で
ある。 次いで(c)該修正された骨パターンより、骨の断
面の外形を楕円とし帯状の骨皮質をもちその内側
に骨密度が減少する部分を持つ骨モデルを設定し
て、骨密度分布を求める。骨密度分布の求め方を
第５図に基いて説明する。 前述した方法により得られる修正した骨パター
ンを第５図の如き楕円形の骨モデルにあてはめ
る。ここで用いる楕円形は骨の一断面の外形であ
る。楕円形は長径／短径が、１＜長径／短径≦
1.4の範囲にあるものが望ましく、中手骨を対象
とする場合には1.25あるいはその前後のものが、
現実の骨によく適合する点で望ましい。上腕骨あ
るいは前腕骨近位部を対象とする場合には、長
径／短径が１に近い値の楕円形が望ましい。次い
で骨パターンをｎ等分する。第５図では骨パター
ンの片側（尺側）を100等分した例が記載されて
いる。そして各層における骨密度をμ1′，μ2′μ3′
，
……μi′……μ′100とする。また第５図の如く、骨
モデルをｎ個の帯状の骨皮質を有するように設定
し、それぞれの帯における骨密度は一定とする。
なわち骨の断面の外周より等距離にある層の骨密
度は同一とする。ｎ等分した各位置における骨モ
デルの断面の長さは次のように表わすことができ
る。 ２・I1，1 ２・（I2，1＋I2，2） ２・（I3，1＋I3，2＋I3，3） 〓 〓 ２・（Ii，1＋Ii，2＋……＋Ii，ｊ＋……＋Ii，ｉ） 〓 〓 I1，1，I2，1……等の値は計算により求まる。こ
れらの値から以下の如き漸化式が求まる。この漸
化式は第５図に示すように100等分した場合の式
である。 GS1＝２・（i1，1・μ′1） GS2＝２・（I2，1・μ′2＋I2，2・μ′1） GSi＝２・（Ii，1・μi′＋…＋Ii，ｊ・μi−ｊ＋1 ＋…＋Ii，ｉ・μ′1） GS100＝２・（I100，1・μ′100＋…＋I100，ｋ・μ
′100 −ｋ−1＋I100，100・μ′1） この漸化式を解くことによつて各位置での骨密
度（μi′）が求まり、これをグラフ化することに
よつて骨密度分布が得られる。第５図の点線が骨
密度分布の例を示す。本発明では、上記漸化式を
解く際に、骨皮質の外側の部分の骨密度は一定
（μ1′）＝μ2′＝μ3′＝……）とみなして解くこと
がで
きる。ななわち骨皮質の外側の部分で骨密度が均
質と見られる部分の骨パターンを選び、その部分
の骨密度を一定として、上記漸化式を解くことも
できる。また骨パターンの両端部分は散乱線の影
響で誤差が大きいので、無視することができ、ま
た無視することによつて、より正しい値が求ま
る。本法を例えば橈骨速位部の如く、明らかに断
面が楕円として近似できない場合には、骨パター
ンの対象化を略し、橈側のみの骨パターンを用い
て、楕円として近似できる部分のみを骨モデルに
あてはめて計算することもできる。 ここで得られる骨密度（μi′）は、Ｘ線撮影し
て骨パターンを求めた場合には、Ｘ線像の陰影度
をアルミニウムステツプエツジあるいはアルミニ
ウムスロープの陰影度と比較した値より求まる
GS値に基いて得られるものであり、Ｘ線減弱係
数に相当するものである。そしてかかる骨密度す
なわちＸ線減弱係は、実際の骨のＸ線減弱係数
（μb）から、骨のまわりの軟部組織のＸ線減弱係
数（μs）を減じた値と、アルミニウムのＸ線減弱
係数（μal）との比として、すなわち μ′i＝（μb−μs）／μal として求まるものである。 骨の各部の正確な骨密度を求めるには、コンピ
ユータートモグラフイーのようなスキヤンニング
を要するが、本発明の如く骨パターンを修正し、
骨モデルにあてはめれば骨密度の分布を近似的に
求めることができる。 かくして得られる骨密度分布は、従来のMD法
によつて得られる骨密度分布に比べて、より正確
に実際の骨の状態を反映し得る指標であり、この
骨密度分布を用いることによつて、より客観的な
骨の評価が可能である。例えば骨疾患を有する患
者の骨密度分布と、健常人の骨密度分布を比較す
ることによつて骨病変の類別判定が可能であり、
また骨密度が不均一である腎性骨証の病態をより
適確に判定することができ、骨の成長度，老化度
等の判定にも応用できる。また得られる骨密度分
布より第６図に示した如き断面二次モーメント、
または断面二次極モーメント（Clin.Orthop.Rel.
Res.，149，268，1980）を求めることができる。
これらの指標は骨の強度の指標ともなり得るもの
であり、これらの指標を併用することによつてよ
り正確に骨の評価を行なうことができる。 次いで第３の発明、すなわち (a) 四肢長管骨をＸ線撮影して得られる骨のＸ線
写真陰影度を一定の間隔で測定して各部位の骨
パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行なつて
各部位の骨パターンを修正し、 (c) 該修正された骨パターンより、骨の断面の外
形を楕円として帯状の骨皮質をもちその内側に
骨密度が減少する部分を持つ骨モデルを設定し
て、各部位の骨密度分布を求め、 (d) 各部位の骨密度分布を、骨密度分布ごとにそ
の値の大小により色分けして、 (e) 四肢長管骨のＸ線像を、骨密度分布の色分け
された像に変換して、 (f) 該像により骨の評価を行なう ことを特徴とする骨の評価方法について説明す
る。 第３の発明では、第２の発明と全く同様にして
骨密度分布を求める。ただし第３の発明において
は、上記(a)で述べたように、四肢長管骨をＸ線撮
影して得られる骨のＸ線写真陰影度を一定の間隔
で測定して各部位の骨パターンを求め、この骨パ
ターンより各部位ごとに骨密度分布を求める。こ
れらの方法を７図に基いて説明する。 第７図のは第２中手骨のＸ線像を示す略図で
あり斜線はＸ線の陰影を示す。第３の発明では、
Ｘ線像を一定の間隔に区分し各部位ごとに、骨パ
ターンを求める。一定の間隔としては通常２mm以
下が採用される。またこの場合、特に末梢骨の骨
幹部を測定対象とするのが好ましい。各部位ごと
に得られる骨パターンより、第２の発明と全く同
様にして各部位ごとの骨密度分布を求める。第７
図のは骨密度分布の略図を示す。第３の発明で
は、得られた骨密度分布を、その骨密度（μ′）の
値の大小によつて色分けする。すなわち例えば骨
密度が高い部分（あるいは正常な部分）を例えば
赤系統の色、骨密度が低下していく部分（あるい
は軽度ないし中等度の萎縮を示す部分）、例えば
ピンク，黄色，緑と色別し、骨密度が極度に低下
した部分（あるいは高度の萎縮を示す部分）を、
例えば青系統の色で色別する。この操作を、各部
位で得られる骨密度分布について行なう。次いで
この色別化された骨密度分布により、四肢長管骨
のＸ線像を骨密度分布の色別された像に変換す
る。第７図のは変換された像の略図である。こ
れらの操作は、通常コンピユーターを用いてグラ
フイツクデイスプレイすることによつて行なわれ
る。 また上記(a)〜(c)の操作を、従来のMD法による
骨密度分布の求め方によつて行ない（骨代謝，第
14巻，91〜104頁，1981）、得られる骨密度分布を
本発明の(d)〜(e)の操作に付して、四肢の長管骨の
Ｘ線像等より得られる骨パターンを骨密度分布の
色別された像に変換し、この像に基いて骨を評価
することも可能である。 かくして得られる像は、骨の一断面ではなく、
骨全体の骨密度分布を表わしており、この像によ
つて骨をより正確に、り現実的に評価することが
できる。例えば本発明の像により、老人性骨粗鬆
症では骨皮質は菲薄化し、近部位から遠部位へほ
ぼ均質に骨密度の低下がみられる。慢性関節リウ
マチでは、関節近傍に高密度な骨萎縮がみられ
る。 〈発明の効果〉 本発明は、従来のMD法に改良を加えた骨の評
価方法を提供するものであり、本発明の方法によ
れば、骨の状態をより適確に客観的に評価するこ
とができ、本発明の意義は大きい。 〈実施例〉 実施例 １ 62才女性の骨粗鬆症患者の第２中手骨のＸ線写
真より求めた骨密度分布を第７図に示す。骨密度
を求める際に骨皮質の外側の部分の骨密度は一定
として求めた。第８図のａはデンシトメリーによ
り測定した値から得られた骨パターンである。ｂ
の「―」は骨幅を10mmへ換算して求めた骨パター
ンである。「―・―・―」は対称化及び平滑化し
て求めた修正された骨パターンである。 「……」は骨密度分布である。 実施例 ２ 20才から40才までの男子19名、女子18名の健常
者を対象として、実施例１と同様にして骨密度分
布を求めた。このとき得られた骨皮質外側の骨密
度（μ′）の結果を表にした。 【表】
更に詳細には本発明は、四肢長管骨のＸ線写真陰
影度を測定して骨パターンを求め、得られる骨パ
ターンを修正し、修正された骨パターンに基い
て、各種の評価を行なう骨のＸ線写真像の評価方
法に関する。 〈従来技術〉 人間の骨の発育状態，老化度、あるいは骨粗鬆
症，骨軟化症等の骨病変の種類などを評価する方
法としてMD法が知られている（骨代謝，第13
巻，187〜195頁，1980；骨代謝，第14巻，91〜
104頁，1981等参照）。 MD法は、第２中手骨の骨幹部中央のＸ線像の
陰影度をマイクロデンシトメーターを用いて測定
し、コンピユーターを併用して陰影より、骨皮質
幅指数（MCL），骨幅(D)，最高骨濃度（GSnax），
最低骨濃度（GSnio），平均骨濃度（ΣGS／D.）
等の各指標を求め、この指標によつて、骨萎縮度
等を判定する方法である。 〈発明が解決しようとする問題点〉 MD法にあつては、MD法で得られる各指標、
例えば最高骨濃度（GSnax），骨皮質幅指数
（ACI）等の数値が、実際の骨の状態を反映しな
い場合があり、これらの指標によつては骨の評価
が困難なことも起り得る。たとえば最高骨濃度
（GSnax）は骨皮質の骨密度の指標とされている
が、骨密度が同じ骨の場合でも骨の直径や骨皮質
幅が異なれば違う数値を示し、また骨皮質が均質
でない場合には、その意味することが不明確とな
る。 フオトン・アブソープシオメトリーによる従来
の方法では、骨塩量（BMC）を求求めることし
か知られていない。 〈問題点を解決するための手段〉 本発明者らは、直接的に骨の状態を反映し得る
指標を求めることを目的として鋭意研究した結
果、第２中手骨のＸ線像の陰影度を、例えば従来
のMD法と同様にして、測定して骨パターンを求
め、該骨パターンの平滑化及び対称化を行なつて
骨パターンを修正し、該修正された骨パターンよ
り、各指標を求めることによつて得られる各指標
が、骨の状態をよく反映し得ることを見出した。 本発明者らはかかる知見に基いて更に鋭意研究
した結果、本発明を完成したものである。 本発明は第１に、 (a) 四肢長管骨をＸ線撮影して得られる骨のＸ線
写真像陰影度を測定して骨パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行なつて
骨パターンを修正し、 (c) 該修正された骨パターンより、骨皮質幅指数
（MCI），骨幅(D)，最高骨濃度（GSnax），最低
骨濃度（GSnio），骨塩量（ΣGS）及び平均骨濃
度（ΣGS／ＤあるいはΣGS／D2）より選ばれ
る少なくとも１種の指標を求めて、 (d) 該指標により骨の評価を行なうことを特徴と
する骨のＸ線写真像の評価方法。 第２に、 (a) 四肢長管骨をＸ線撮影して得られる骨のＸ線
写真像陰影度を測定して骨パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行なつて
骨パターンを修正し、 (c) 該修正された骨パターンより、骨の断面の外
形を楕円とし帯状の骨皮質をもちその内側に骨
密度が減少する部分を持つ骨モデルを設定し
て、骨密度分布を求め、 (d) 該骨密度分布により骨を評価することを特徴
とする骨のＸ線写真像の評価方法。 第３に、 (a) 四肢長管骨をＸ線撮影して得られる骨のＸ線
写真像陰影度を一定の間隔で測定して各部位の
骨パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行なつて
各部位の骨パターンを修正し、 (c) 該修正された骨パターンより、骨の断面の外
形を楕円とし帯状の骨皮質をもちその内側に骨
密度が減少する部分を持つ骨モデルを設定し
て、各部位の骨密度分布を求め、 (d) 各部位の骨密度分布を、骨密度分布ごとにそ
の値の大小により色分けして、 (e) 四肢長管骨のＸ線像を、骨密度分布の像に変
換して、 (f) 該像により骨の評価を行なうことを特徴とす
る骨のＸ線写真像の評価方法である。 以下第１の発明について説明する。 先ず、(a)、四肢長管骨をＸ線撮影して得られる
骨のＸ線写真陰影度を測定して骨パターンを求め
る。あるいは四肢長管骨からフオトン・アブソー
プシオメトリーにより骨パターンを求める。 四肢長管骨としては、骨のまわりの軟部組織の
厚さが薄く平均化している部分の骨が望ましく、
例えば第１〜第５の中手骨、上腕骨などが挙げら
れる。なかでも第２中手骨の骨幹部中央が好適で
ある。四肢長管骨をＸ線撮影してＸ線の像を得
る。Ｘ線撮影して得られる骨のＸ線写真陰影度を
測定して骨パターンを求める。骨パターンを求め
るには、例えばMD法によれば、Ｘ線像の陰影度
を、例えば１段の高さ１mmで20段（最低の高さ１
mm，最高の高さ20mm）のアルミ階段（alminium
step wedge）、あるいはアルミニウムロープのＸ
線像とともにデンシトメーターで測定する。骨の
近位端と速位端との中間点の陰影度を測定するの
が望ましい。得られる陰影度は、光学密度
（OD）を目盛つた用紙に拡大して記録し、骨の
各点における吸光度（OD）をアルミ階段に段数
に変換補正（GS値）して第１図に示すような骨
パターンを作図する（骨代謝，第13巻，187〜195
頁，1980；骨代謝，第14巻，91〜104頁，1981；
骨代謝，第14巻，319〜325頁，1981）。あるいは、
骨のＸ線像を、アルミ階段、アルミニウムスロー
プのＸ線像とともにテレビカメラにより読み取つ
て骨パターンを求めることもできる。 次いで(b)、得られる骨パターンの対称化及び平
滑化を行なつて骨パターンを修正する。 骨パターンの対称化を第２図に基いて説明す
る。先ず上述した方法で得られる骨パターンに中
心線（第２図の一点鎖線）を引く。中心線より左
右に等間隔（Ｘ）の位置におけるGS値（GSx1及
びGSx2）を求め、これらの平均値〔GSx＝
（GSx1＋GSx2）／２〕を、中心線より左右に等
間隔の位置のGS値（GSx）とする。Ｘを変えて
この操作をくり返すことにより、対称化された骨
パターン（第２図の点線）を求める。次いで骨パ
ターンの平滑化を行なう。平滑化の方法を第３図
に基いて説明する。先ず上述した方法で得られる
対称化された骨パターンを等分化する。次いで等
分化された骨パターンを、例えば５点ずつ区分す
る（Xi1，Xi2，Xi3，Xi4，Xi5）。次いで各位置
におけるGS値（GSx1，GSxi2，GSxi3，GSxi4，
GSxi5）を求めその平均値を算出して、その値
を、５区分の中心であるXi3におけるGS値とす
る。この操作を順次行ない得られるGS値を結ぶ
ことに平滑化された骨パターンが得られる。この
平滑化は移動平均法に準じた方法である。 本発明では、上記対称化，平滑化の操作の順は
いずれでもよいが、対称化を先ず行ない、次いで
平滑化するのが望ましい。 かくして得られる修正された骨パターンより、
各指標を求める。各指標について例えば第２中手
骨に本法を用いた場合に得られる第４図に基いて
説明する。 第４図において、Ｄは骨幅，ｄは骨髄幅，d1は
尺側骨皮質幅，d2は橈側骨皮質幅を示している。
骨幅(D)，尺側骨皮質幅（d1），橈側骨皮質幅（d2）
より骨皮質幅指数（MCI）が求まれる（MCI＝
（d1＋d2）／Ｄ）。GSnioは最低骨濃度であり、
GSmax1（尺側骨皮質におけるGS値の極大値）と
GSmax2（橈側骨皮質におけるGS値の極大値）と
の平均値（本発明の場合、GSmax1とGSmax2と
は、骨パターンの対称化により等しい）より最高
濃度（GSmax）〔GSmax＝1/2（GSmax1＋
GSmax2）〕が求まる。また各GS値を積分して得
られる面積が骨塩量（ΣGS）であり、これより
平均骨濃度（ΣGS／ＤまたはΣGS／D2）が求ま
る。 かくして得られる指標、骨皮質幅指数
（MCI），骨幅(D)，最高骨濃度（GSnax），最低骨
濃度（GSnio），骨塩量（ΣGS）及び平均骨濃度
（ΣGS／ＤまたはΣGS／D2）より選ばれる少なく
とも１種の指標によつて骨の評価を行なうことが
できる。これらの指標、例えば骨皮質幅指数
（MCI）は、従来のMD法で用いられるMCIに比
べて実際の骨の状態をより適確に表わした値であ
る。したがつてこれらの指標により、より正確な
骨の評価が可能である。骨の評価の例としては、
例えば骨萎縮度の判定（骨代謝，第13巻，187〜
195頁，1980），骨成長の判定（骨代謝，第14巻，
319〜325頁，1981），骨病変の類別判定（特開昭
59−8935，特開昭59−11840号公報），骨改善効果
の判定（特開昭59−49743号公報）などが挙げら
れる。本発明で得られる指標を用いてこれら判定
を行なうことが可能である。 次いで第２の発明、すなわち、 (a) 四肢長管骨をＸ線撮影して得られる骨のＸ線
写真陰影度を測定して骨パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行なつて
骨パターンを修正し、 (c) 該修正された骨パターンより、骨の断面の外
形を楕円とし帯状の骨皮質をもちその内側に骨
密度が減少する部分を持つ骨モデルを設定し
て、密度分布を求め、 (d) 該骨密度分布により骨を評価することを特徴
とする骨の評価方法について説明する。 第２の発明では先ず、(a)四肢長管骨をＸ線撮影
して得られる骨のＸ線写真陰影度を測定して骨パ
ターンを求め、(b)該骨パターンの平滑化及び対称
化を行なつて骨パターンを修正する。これら(a)，
(b)の操作は第１の発明において説明したと同様で
ある。 次いで(c)該修正された骨パターンより、骨の断
面の外形を楕円とし帯状の骨皮質をもちその内側
に骨密度が減少する部分を持つ骨モデルを設定し
て、骨密度分布を求める。骨密度分布の求め方を
第５図に基いて説明する。 前述した方法により得られる修正した骨パター
ンを第５図の如き楕円形の骨モデルにあてはめ
る。ここで用いる楕円形は骨の一断面の外形であ
る。楕円形は長径／短径が、１＜長径／短径≦
1.4の範囲にあるものが望ましく、中手骨を対象
とする場合には1.25あるいはその前後のものが、
現実の骨によく適合する点で望ましい。上腕骨あ
るいは前腕骨近位部を対象とする場合には、長
径／短径が１に近い値の楕円形が望ましい。次い
で骨パターンをｎ等分する。第５図では骨パター
ンの片側（尺側）を100等分した例が記載されて
いる。そして各層における骨密度をμ1′，μ2′μ3′
，
……μi′……μ′100とする。また第５図の如く、骨
モデルをｎ個の帯状の骨皮質を有するように設定
し、それぞれの帯における骨密度は一定とする。
なわち骨の断面の外周より等距離にある層の骨密
度は同一とする。ｎ等分した各位置における骨モ
デルの断面の長さは次のように表わすことができ
る。 ２・I1，1 ２・（I2，1＋I2，2） ２・（I3，1＋I3，2＋I3，3） 〓 〓 ２・（Ii，1＋Ii，2＋……＋Ii，ｊ＋……＋Ii，ｉ） 〓 〓 I1，1，I2，1……等の値は計算により求まる。こ
れらの値から以下の如き漸化式が求まる。この漸
化式は第５図に示すように100等分した場合の式
である。 GS1＝２・（i1，1・μ′1） GS2＝２・（I2，1・μ′2＋I2，2・μ′1） GSi＝２・（Ii，1・μi′＋…＋Ii，ｊ・μi−ｊ＋1 ＋…＋Ii，ｉ・μ′1） GS100＝２・（I100，1・μ′100＋…＋I100，ｋ・μ
′100 −ｋ−1＋I100，100・μ′1） この漸化式を解くことによつて各位置での骨密
度（μi′）が求まり、これをグラフ化することに
よつて骨密度分布が得られる。第５図の点線が骨
密度分布の例を示す。本発明では、上記漸化式を
解く際に、骨皮質の外側の部分の骨密度は一定
（μ1′）＝μ2′＝μ3′＝……）とみなして解くこと
がで
きる。ななわち骨皮質の外側の部分で骨密度が均
質と見られる部分の骨パターンを選び、その部分
の骨密度を一定として、上記漸化式を解くことも
できる。また骨パターンの両端部分は散乱線の影
響で誤差が大きいので、無視することができ、ま
た無視することによつて、より正しい値が求ま
る。本法を例えば橈骨速位部の如く、明らかに断
面が楕円として近似できない場合には、骨パター
ンの対象化を略し、橈側のみの骨パターンを用い
て、楕円として近似できる部分のみを骨モデルに
あてはめて計算することもできる。 ここで得られる骨密度（μi′）は、Ｘ線撮影し
て骨パターンを求めた場合には、Ｘ線像の陰影度
をアルミニウムステツプエツジあるいはアルミニ
ウムスロープの陰影度と比較した値より求まる
GS値に基いて得られるものであり、Ｘ線減弱係
数に相当するものである。そしてかかる骨密度す
なわちＸ線減弱係は、実際の骨のＸ線減弱係数
（μb）から、骨のまわりの軟部組織のＸ線減弱係
数（μs）を減じた値と、アルミニウムのＸ線減弱
係数（μal）との比として、すなわち μ′i＝（μb−μs）／μal として求まるものである。 骨の各部の正確な骨密度を求めるには、コンピ
ユータートモグラフイーのようなスキヤンニング
を要するが、本発明の如く骨パターンを修正し、
骨モデルにあてはめれば骨密度の分布を近似的に
求めることができる。 かくして得られる骨密度分布は、従来のMD法
によつて得られる骨密度分布に比べて、より正確
に実際の骨の状態を反映し得る指標であり、この
骨密度分布を用いることによつて、より客観的な
骨の評価が可能である。例えば骨疾患を有する患
者の骨密度分布と、健常人の骨密度分布を比較す
ることによつて骨病変の類別判定が可能であり、
また骨密度が不均一である腎性骨証の病態をより
適確に判定することができ、骨の成長度，老化度
等の判定にも応用できる。また得られる骨密度分
布より第６図に示した如き断面二次モーメント、
または断面二次極モーメント（Clin.Orthop.Rel.
Res.，149，268，1980）を求めることができる。
これらの指標は骨の強度の指標ともなり得るもの
であり、これらの指標を併用することによつてよ
り正確に骨の評価を行なうことができる。 次いで第３の発明、すなわち (a) 四肢長管骨をＸ線撮影して得られる骨のＸ線
写真陰影度を一定の間隔で測定して各部位の骨
パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行なつて
各部位の骨パターンを修正し、 (c) 該修正された骨パターンより、骨の断面の外
形を楕円として帯状の骨皮質をもちその内側に
骨密度が減少する部分を持つ骨モデルを設定し
て、各部位の骨密度分布を求め、 (d) 各部位の骨密度分布を、骨密度分布ごとにそ
の値の大小により色分けして、 (e) 四肢長管骨のＸ線像を、骨密度分布の色分け
された像に変換して、 (f) 該像により骨の評価を行なう ことを特徴とする骨の評価方法について説明す
る。 第３の発明では、第２の発明と全く同様にして
骨密度分布を求める。ただし第３の発明において
は、上記(a)で述べたように、四肢長管骨をＸ線撮
影して得られる骨のＸ線写真陰影度を一定の間隔
で測定して各部位の骨パターンを求め、この骨パ
ターンより各部位ごとに骨密度分布を求める。こ
れらの方法を７図に基いて説明する。 第７図のは第２中手骨のＸ線像を示す略図で
あり斜線はＸ線の陰影を示す。第３の発明では、
Ｘ線像を一定の間隔に区分し各部位ごとに、骨パ
ターンを求める。一定の間隔としては通常２mm以
下が採用される。またこの場合、特に末梢骨の骨
幹部を測定対象とするのが好ましい。各部位ごと
に得られる骨パターンより、第２の発明と全く同
様にして各部位ごとの骨密度分布を求める。第７
図のは骨密度分布の略図を示す。第３の発明で
は、得られた骨密度分布を、その骨密度（μ′）の
値の大小によつて色分けする。すなわち例えば骨
密度が高い部分（あるいは正常な部分）を例えば
赤系統の色、骨密度が低下していく部分（あるい
は軽度ないし中等度の萎縮を示す部分）、例えば
ピンク，黄色，緑と色別し、骨密度が極度に低下
した部分（あるいは高度の萎縮を示す部分）を、
例えば青系統の色で色別する。この操作を、各部
位で得られる骨密度分布について行なう。次いで
この色別化された骨密度分布により、四肢長管骨
のＸ線像を骨密度分布の色別された像に変換す
る。第７図のは変換された像の略図である。こ
れらの操作は、通常コンピユーターを用いてグラ
フイツクデイスプレイすることによつて行なわれ
る。 また上記(a)〜(c)の操作を、従来のMD法による
骨密度分布の求め方によつて行ない（骨代謝，第
14巻，91〜104頁，1981）、得られる骨密度分布を
本発明の(d)〜(e)の操作に付して、四肢の長管骨の
Ｘ線像等より得られる骨パターンを骨密度分布の
色別された像に変換し、この像に基いて骨を評価
することも可能である。 かくして得られる像は、骨の一断面ではなく、
骨全体の骨密度分布を表わしており、この像によ
つて骨をより正確に、り現実的に評価することが
できる。例えば本発明の像により、老人性骨粗鬆
症では骨皮質は菲薄化し、近部位から遠部位へほ
ぼ均質に骨密度の低下がみられる。慢性関節リウ
マチでは、関節近傍に高密度な骨萎縮がみられ
る。 〈発明の効果〉 本発明は、従来のMD法に改良を加えた骨の評
価方法を提供するものであり、本発明の方法によ
れば、骨の状態をより適確に客観的に評価するこ
とができ、本発明の意義は大きい。 〈実施例〉 実施例 １ 62才女性の骨粗鬆症患者の第２中手骨のＸ線写
真より求めた骨密度分布を第７図に示す。骨密度
を求める際に骨皮質の外側の部分の骨密度は一定
として求めた。第８図のａはデンシトメリーによ
り測定した値から得られた骨パターンである。ｂ
の「―」は骨幅を10mmへ換算して求めた骨パター
ンである。「―・―・―」は対称化及び平滑化し
て求めた修正された骨パターンである。 「……」は骨密度分布である。 実施例 ２ 20才から40才までの男子19名、女子18名の健常
者を対象として、実施例１と同様にして骨密度分
布を求めた。このとき得られた骨皮質外側の骨密
度（μ′）の結果を表にした。 【表】
第１図は、Ｘ線像の陰影度からMD法により求
まる骨パターンを示す。第２図は本発明で用いら
れる対称化について示したものであり、第３図は
平滑化について示したものである。第４図は修正
された骨パターンより求まる各指標について示し
た。第５図は修正された骨パターンを骨モデルに
あてはめた図である。第６図は断面二次モーメン
ト及び断面二次極モーメントを示す図であり、第
７図は骨密度分布を色別けして、Ｘ線像を骨密度
分布の像に変換する様子を示す。第８図は本発明
方法の実例を示す。
まる骨パターンを示す。第２図は本発明で用いら
れる対称化について示したものであり、第３図は
平滑化について示したものである。第４図は修正
された骨パターンより求まる各指標について示し
た。第５図は修正された骨パターンを骨モデルに
あてはめた図である。第６図は断面二次モーメン
ト及び断面二次極モーメントを示す図であり、第
７図は骨密度分布を色別けして、Ｘ線像を骨密度
分布の像に変換する様子を示す。第８図は本発明
方法の実例を示す。
Claims (1)
- 【特許請求の範囲】 １ (a) 四肢長管骨をＸ線撮影して得られる骨の
Ｘ線写真像の陰影度を測定して骨パターンを求
め、 (b) 該骨パターンの平滑化及び対称化を行つて骨
パターンを修正し、 (c) 該修正された骨パターンより、骨皮質幅指数
（MCI），骨幅(D)，最高骨濃度（GSnax），最低
骨濃度（GSnio），骨塩量（ΣGS）及び平均骨濃
度（ΣGS／ＤあるいはΣGS／D2）より選ばれ
る少なくとも１種の指標を求めて、 (d) 該指標により骨の評価を行うことを特徴とす
る骨のＸ線写真像の評価方法。 ２ (a) 四肢長管骨をＸ線撮影して得られる骨の
Ｘ線写真像の陰影度を測定して骨パターンを求
め、 (b) 該骨パターンの平滑化及び対称化を行つて骨
パターンを修正し、 (c) 該修正された骨パターンより、骨の断面の外
形を楕円とし帯状の骨皮質をもちその内側に骨
密度が減少する部分を持つ骨モデルを設定して
骨密度分布を求め、 (d) 該骨密度分布により骨を評価することを特徴
とする骨のＸ線写真像の評価方法。 ３ (a) 四肢長管骨をＸ線撮影して得られる骨の
Ｘ線写真像の陰影度を一定の間隔で測定して各
部位の骨パターンを求め、 (b) 該骨パターンの平滑化及び対称化を行つて骨
パターンを修正し、 (c) 該修正された骨パターンより、骨の断面の外
形を楕円とし帯状の骨皮質をもちその内側に骨
密度が減少する部分を持つ骨モデルを設定し
て、各部位の骨密度分布を求め、 (d) 各部位の骨密度分布を、骨密度分布ごとにそ
の値の大小により色分けして、 (e) 四肢長管骨のＸ線写真像を、骨密度分布の色
分けされた像に変換して、 (f) 該像により骨の評価を行うことを特徴とする
骨のＸ線写真像の評価方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59230077A JPS61109557A (ja) | 1984-11-02 | 1984-11-02 | 骨のｘ線写真像の評価方法 |
| US06/793,387 US4721112A (en) | 1984-11-02 | 1985-10-31 | Bone evaluation method |
| EP85307943A EP0180482A3 (en) | 1984-11-02 | 1985-11-01 | Bone evaluation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59230077A JPS61109557A (ja) | 1984-11-02 | 1984-11-02 | 骨のｘ線写真像の評価方法 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1214163A Division JPH02126837A (ja) | 1989-08-22 | 1989-08-22 | 骨の評価装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS61109557A JPS61109557A (ja) | 1986-05-28 |
| JPH0331061B2 true JPH0331061B2 (ja) | 1991-05-02 |
Family
ID=16902184
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59230077A Granted JPS61109557A (ja) | 1984-11-02 | 1984-11-02 | 骨のｘ線写真像の評価方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US4721112A (ja) |
| EP (1) | EP0180482A3 (ja) |
| JP (1) | JPS61109557A (ja) |
Families Citing this family (57)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0228785B1 (en) * | 1985-11-11 | 1993-01-13 | Teijin Limited | Bone evaluation method |
| WO1990009761A1 (fr) * | 1989-02-23 | 1990-09-07 | Teijin Limited | Procede d'evaluation de la forme d'un os, appareil utilise et systeme d'evaluation d'os |
| US5090040A (en) * | 1989-03-10 | 1992-02-18 | Expert Image Systems, Inc. | Data acquisition system for radiographic imaging |
| US5138553A (en) * | 1989-03-10 | 1992-08-11 | Expert Image Systems, Inc. | Method and apparatus for bone mineral measurement using a measurement region defined at a selected distance from a styloid bone tip |
| JPH03230670A (ja) * | 1990-02-05 | 1991-10-14 | Teijin Ltd | 画像読み取り装置及び骨計測装置 |
| US5509042A (en) * | 1991-02-13 | 1996-04-16 | Lunar Corporation | Automated determination and analysis of bone morphology |
| US5533084A (en) * | 1991-02-13 | 1996-07-02 | Lunar Corporation | Bone densitometer with improved vertebral characterization |
| US5577089A (en) * | 1991-02-13 | 1996-11-19 | Lunar Corporation | Device and method for analysis of bone morphology |
| US5247934A (en) * | 1991-08-09 | 1993-09-28 | Trustees Of The University Of Pennsylvania | Method and apparatus for diagnosing osteoporosis with MR imaging |
| DE69329168T2 (de) * | 1992-05-20 | 2001-04-05 | Aloka Co. Ltd., Mitaka | Vorrichtung zur Bestimmung der Eigenschaften von Knochen |
| DE69320418T2 (de) * | 1992-06-04 | 1999-01-07 | Teijin Ltd., Osaka | Verfahren und Gerät zur Knochen-Messung |
| US5585603A (en) * | 1993-12-23 | 1996-12-17 | Design Systems, Inc. | Method and system for weighing objects using X-rays |
| AU3341295A (en) * | 1994-08-29 | 1996-03-22 | Eskofot A/S | A method of estimation |
| US5673303A (en) * | 1995-04-19 | 1997-09-30 | Wright State University | Method and apparatus for the evaluation of structural width and density by computed tomography |
| US5782763A (en) * | 1995-06-01 | 1998-07-21 | Cyberlogic Inc. | Electromagnetic bone-assessment apparatus and method |
| US6213934B1 (en) | 1995-06-01 | 2001-04-10 | Hyper3D Corp. | Electromagnetic bone-assessment and treatment: apparatus and method |
| US5696805A (en) * | 1996-09-17 | 1997-12-09 | Eastman Kodak Company | Apparatus and method for identifying specific bone regions in digital X-ray images |
| GB9702202D0 (en) * | 1997-02-04 | 1997-03-26 | Osteometer Meditech As | Diagnosis of arthritic conditions |
| AU8103198A (en) | 1997-07-04 | 1999-01-25 | Torsana Osteoporosis Diagnostics A/S | A method for estimating the bone quality or skeletal status of a vertebrate |
| US6252928B1 (en) | 1998-01-23 | 2001-06-26 | Guard Inc. | Method and device for estimating bone mineral content of the calcaneus |
| EP1078252A4 (en) | 1998-03-02 | 2007-05-30 | Image Analysis Inc | AUTOMATED BONE DENSITY MEASUREMENT BY X-RAY RADIATION |
| KR100338974B1 (ko) * | 1999-03-15 | 2002-05-31 | 최은백, 이찬경 | 턱뼈의 뼈밀도를 확인하기 위한 시뮬레이션 방법 및 이를 실현시키기 위한 프로그램이 기록된 기록매체 |
| DE69921087T2 (de) | 1999-04-20 | 2005-11-24 | Sectra Pronosco A/S | Vorrichtung zur Abschätzung der Knochenmineraldichte mittels Röntgenbildmessungen |
| AU4744500A (en) * | 1999-05-20 | 2000-12-12 | Torsana Osteoporosis Diagnostics A/S | Method and apparatus for selection and evaluation of substances in treatment of bone disorders |
| US6711282B1 (en) * | 1999-10-29 | 2004-03-23 | Compumed, Inc. | Method for automatically segmenting a target bone from a digital image |
| US6246745B1 (en) | 1999-10-29 | 2001-06-12 | Compumed, Inc. | Method and apparatus for determining bone mineral density |
| US6381301B1 (en) | 1999-12-01 | 2002-04-30 | Ronald E. Massie | Dental and orthopedic densitometry modeling system and method |
| US6944262B2 (en) | 1999-12-01 | 2005-09-13 | Massie Ronald E | Dental and orthopedic densitometry modeling system and method |
| US8073101B2 (en) | 1999-12-01 | 2011-12-06 | Massie Ronald E | Digital modality modeling for medical and dental applications |
| US8126112B2 (en) | 1999-12-01 | 2012-02-28 | Massie Ronald E | Osseo classification system and method |
| US6510197B1 (en) | 2000-01-11 | 2003-01-21 | Alara, Inc. | Method and apparatus for osteoporosis screening |
| US20020186818A1 (en) * | 2000-08-29 | 2002-12-12 | Osteonet, Inc. | System and method for building and manipulating a centralized measurement value database |
| US6904123B2 (en) * | 2000-08-29 | 2005-06-07 | Imaging Therapeutics, Inc. | Methods and devices for quantitative analysis of x-ray images |
| CA2419483A1 (en) * | 2000-08-29 | 2002-03-07 | Osteonet.Com, Inc. | Methods and devices for quantitative analysis of x-ray images |
| US7467892B2 (en) | 2000-08-29 | 2008-12-23 | Imaging Therapeutics, Inc. | Calibration devices and methods of use thereof |
| US7660453B2 (en) | 2000-10-11 | 2010-02-09 | Imaging Therapeutics, Inc. | Methods and devices for analysis of x-ray images |
| US8639009B2 (en) * | 2000-10-11 | 2014-01-28 | Imatx, Inc. | Methods and devices for evaluating and treating a bone condition based on x-ray image analysis |
| US8000766B2 (en) * | 2001-05-25 | 2011-08-16 | Imatx, Inc. | Methods to diagnose treat and prevent bone loss |
| US6990222B2 (en) * | 2001-11-21 | 2006-01-24 | Arnold Ben A | Calibration of tissue densities in computerized tomography |
| US8965075B2 (en) | 2002-09-16 | 2015-02-24 | Imatx, Inc. | System and method for predicting future fractures |
| US20040106868A1 (en) * | 2002-09-16 | 2004-06-03 | Siau-Way Liew | Novel imaging markers in musculoskeletal disease |
| US7840247B2 (en) * | 2002-09-16 | 2010-11-23 | Imatx, Inc. | Methods of predicting musculoskeletal disease |
| US8600124B2 (en) * | 2004-09-16 | 2013-12-03 | Imatx, Inc. | System and method of predicting future fractures |
| EP1605824A2 (en) | 2003-03-25 | 2005-12-21 | Imaging Therapeutics, Inc. | Methods for the compensation of imaging technique in the processing of radiographic images |
| EP1663002A4 (en) * | 2003-09-19 | 2007-11-28 | Imaging Therapeutics Inc | BONE STRUCTURE PROGNOSTIC METHOD AND SIMULATED BONE REMODELING |
| US8290564B2 (en) * | 2003-09-19 | 2012-10-16 | Imatx, Inc. | Method for bone structure prognosis and simulated bone remodeling |
| DE102004033989B4 (de) * | 2004-07-14 | 2015-08-13 | Siemens Aktiengesellschaft | Verfahren zur Messung der dreidimensionalen Dichteverteilung in Knochen |
| WO2006039358A2 (en) * | 2004-09-30 | 2006-04-13 | The Regents Of The University Of California | Method for assessment of the structure-function characteristics of structures in a human or animal body |
| US7959742B2 (en) * | 2007-07-11 | 2011-06-14 | Whirlpool Corporation | Outer support body for a drawer-type dishwasher |
| US7706501B2 (en) * | 2007-09-07 | 2010-04-27 | Carestream Health, Inc. | Method and apparatus for measuring long bone density of small-animals |
| US8597301B2 (en) * | 2007-10-19 | 2013-12-03 | David Mitchell | Cannula with lateral access and directional exit port |
| US20090105775A1 (en) * | 2007-10-19 | 2009-04-23 | David Mitchell | Cannula with lateral access and directional exit port |
| US8939917B2 (en) * | 2009-02-13 | 2015-01-27 | Imatx, Inc. | Methods and devices for quantitative analysis of bone and cartilage |
| US8280138B2 (en) * | 2009-09-10 | 2012-10-02 | General Electric Company | System and method for performing bone densitometer measurements |
| US9848818B1 (en) | 2013-08-09 | 2017-12-26 | O.N.Diagnostics, LLC | Clinical assessment of fragile bone strength |
| US11850061B2 (en) | 2013-08-09 | 2023-12-26 | O.N.Diagnostics, LLC | Clinical assessment of fragile bone strength |
| US10667777B2 (en) * | 2014-11-06 | 2020-06-02 | University Of Utah Research Foundation | Radiograph density detection device |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US2898466A (en) * | 1954-03-15 | 1959-08-04 | Halliburton Oil Well Cementing | Density determining device |
| US3180985A (en) * | 1962-05-14 | 1965-04-27 | Electronic Automation Systems | Standardization of radiation-absorption type density gages |
| US3483565A (en) * | 1968-05-21 | 1969-12-09 | Cedars Of Lebanon Mount Sinai | Color adapter for multidetector scanner |
| US3715588A (en) * | 1970-10-26 | 1973-02-06 | Norland Corp | Bone mineral analyzer |
| SE373279B (sv) * | 1973-05-11 | 1975-02-03 | Atomenergi Ab | Forfarande och anordning for undersokning av en del av en biologisk vevnad |
| US3996471A (en) * | 1975-03-11 | 1976-12-07 | Nasa | Method and system for in vivo measurement of bone tissue using a two level energy source |
| US4515165A (en) * | 1980-02-04 | 1985-05-07 | Energy Conversion Devices, Inc. | Apparatus and method for detecting tumors |
| JPS58206719A (ja) * | 1982-05-27 | 1983-12-02 | 帝人株式会社 | Ｘ線フィルムにおける像による骨年令評価方法 |
| JPS598935A (ja) * | 1982-07-08 | 1984-01-18 | 帝人株式会社 | Ｘ線フィルムにおける像の処理装置 |
-
1984
- 1984-11-02 JP JP59230077A patent/JPS61109557A/ja active Granted
-
1985
- 1985-10-31 US US06/793,387 patent/US4721112A/en not_active Expired - Fee Related
- 1985-11-01 EP EP85307943A patent/EP0180482A3/en not_active Ceased
Also Published As
| Publication number | Publication date |
|---|---|
| JPS61109557A (ja) | 1986-05-28 |
| US4721112A (en) | 1988-01-26 |
| EP0180482A2 (en) | 1986-05-07 |
| EP0180482A3 (en) | 1988-07-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH0331061B2 (ja) | ||
| Rüegsegger et al. | Quantification of bone mineralization using computed tomography | |
| Ortman et al. | Subtraction radiography and computer assisted densitometric analyses of standardized radiographs: a comparison study with 125I absorptiometry | |
| WO1997042602A1 (en) | A method of estimating skeletal status | |
| CN107485405B (zh) | 一种利用参考模块测量骨密度的装置 | |
| Budoff et al. | Trabecular bone mineral density measurement using thoracic and lumbar quantitative computed tomography | |
| KR100479207B1 (ko) | 엑스선 영상을 이용한 골밀도 측정 방법 | |
| JPH07102210B2 (ja) | 歯槽骨の骨萎縮度評価法 | |
| CN111481225A (zh) | 应用于x线成像设备骨密度测量的体模及x线成像设备 | |
| US7046834B2 (en) | Method for measuring bone mineral density by using X-ray image | |
| Jensen et al. | Assessment of bone mass in the radius by computed tomography | |
| JP2793502B2 (ja) | 骨塩量測定方法および装置 | |
| JPH053299B2 (ja) | ||
| Corcoran et al. | Calculation of cross-sectional geometry of bone from CT images with application in postmenopausal women | |
| JPS639461B2 (ja) | ||
| Bullen et al. | Effect of polyester and Plaster of Paris casts on determination of volumetric bone mineral density assessed by Peripheral Quantitative Computed Tomography (pQCT) | |
| JP2732979B2 (ja) | 骨計測装置 | |
| Trippi et al. | A computer-assisted method for the study of the trabecular bone of the distal radius on conventional radiographs | |
| JPH0340611B2 (ja) | ||
| JPS639462B2 (ja) | ||
| JPH0155016B2 (ja) | ||
| JPH0417052B2 (ja) | ||
| KR20000056709A (ko) | 요골 단순 엑스선 영상을 이용한 골다공증 지표 설정 방법 | |
| JPS639463B2 (ja) | ||
| JPH0484939A (ja) | 骨計測方法及び装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EXPY | Cancellation because of completion of term |