JPH07284020A - 骨密度計測方法 - Google Patents
骨密度計測方法Info
- Publication number
- JPH07284020A JPH07284020A JP6074609A JP7460994A JPH07284020A JP H07284020 A JPH07284020 A JP H07284020A JP 6074609 A JP6074609 A JP 6074609A JP 7460994 A JP7460994 A JP 7460994A JP H07284020 A JPH07284020 A JP H07284020A
- Authority
- JP
- Japan
- Prior art keywords
- femur
- bone density
- region
- bone
- shaft
- 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.)
- Pending
Links
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】
【目的】大腿骨の解剖学的形状に基づいた計測再現性の
高い大腿骨中枢端部の骨密度計測を実現する。 【構成】高低２種類のエネルギで撮影した骨盤から大腿
骨を含む部分の２枚のＸ線画像を用いて、仮の差分画像
を作成し（１０）、計測対象領域と補正領域を設定し
（１１）、補正値を算出し（１２）、差分画像を作成し
（１３）、計測領域を抽出し（１４）、骨塩量および骨
密度を算出する（１５）。
高い大腿骨中枢端部の骨密度計測を実現する。 【構成】高低２種類のエネルギで撮影した骨盤から大腿
骨を含む部分の２枚のＸ線画像を用いて、仮の差分画像
を作成し（１０）、計測対象領域と補正領域を設定し
（１１）、補正値を算出し（１２）、差分画像を作成し
（１３）、計測領域を抽出し（１４）、骨塩量および骨
密度を算出する（１５）。
Description
【０００１】
【産業上の利用分野】本発明は、Ｘ線画像を用いて骨塩
量および骨密度を自動計測する方法、特にＤＥＸＡ（Du
al Energy X-ray Absorptiometry）法を用いた大腿骨中
枢端部の自動計測方法に関する。
量および骨密度を自動計測する方法、特にＤＥＸＡ（Du
al Energy X-ray Absorptiometry）法を用いた大腿骨中
枢端部の自動計測方法に関する。
【０００２】
【従来の技術】骨粗鬆症による骨折は、椎体や大腿骨中
枢端でおきるケースがほとんどである。ここでいう大腿
骨中枢端とは、大腿骨頸部，大転子，ウォード（Ward）
三角であり、大腿骨中枢端部とは、中枢端および、大転
子につながるシャフトの部分、および大腿骨頸部につな
がる骨盤を含めて総称する。これらの部位は、筋肉や脂
肪で構成される軟部組織で囲まれているため、骨密度を
計測する場合は、ＤＥＸＡ法またはＤＰＡ（Dual Photo
n Absorptiometry）法を用いる。
枢端でおきるケースがほとんどである。ここでいう大腿
骨中枢端とは、大腿骨頸部，大転子，ウォード（Ward）
三角であり、大腿骨中枢端部とは、中枢端および、大転
子につながるシャフトの部分、および大腿骨頸部につな
がる骨盤を含めて総称する。これらの部位は、筋肉や脂
肪で構成される軟部組織で囲まれているため、骨密度を
計測する場合は、ＤＥＸＡ法またはＤＰＡ（Dual Photo
n Absorptiometry）法を用いる。
【０００３】骨密度計測の主な目的は、骨粗鬆症である
かどうかを診断することと、骨粗鬆症の治療効果を観察
することである。特に治療効果の観察では、計測値の再
現性の高さが必要である。十分な再現性がない場合、数
回の計測にわたる計測値の変化が治療効果であるか、装
置の変動によるものであるかを判断できない。通常、年
に数％といわれる骨塩量の変化を観察できる高い再現性
を持つ計測方法が必要とされる。
かどうかを診断することと、骨粗鬆症の治療効果を観察
することである。特に治療効果の観察では、計測値の再
現性の高さが必要である。十分な再現性がない場合、数
回の計測にわたる計測値の変化が治療効果であるか、装
置の変動によるものであるかを判断できない。通常、年
に数％といわれる骨塩量の変化を観察できる高い再現性
を持つ計測方法が必要とされる。
【０００４】高い再現性を実現するには、（１）骨密度
を計測する領域の設定再現性，(２)脂肪による計測値の
ずれを補正する領域（以下、補正領域という）の設定再
現性、の二つが満たされなければならない。
を計測する領域の設定再現性，(２)脂肪による計測値の
ずれを補正する領域（以下、補正領域という）の設定再
現性、の二つが満たされなければならない。
【０００５】図２に示すように、補正領域２１は、通常
大腿骨中枢端部を囲む矩形の計測対象領域２０内の軟部
組織部分であり、補正値は補正領域２１の濃度値から算
出する。
大腿骨中枢端部を囲む矩形の計測対象領域２０内の軟部
組織部分であり、補正値は補正領域２１の濃度値から算
出する。
【０００６】従来は、操作者が計測対象領域２０を指示
し、それをもとに補正領域２１の設定，補正値の算出，
計測領域の抽出，各部の面積，骨塩量，骨密度を計測し
ている。計測領域は、図２に示すように大腿骨中枢端部
の、頸部２２，大転子部２３，ウォード三角部２４（以
下それぞれＮ部，Ｔ部，Ｗ部と表わす）の３ヵ所である
（ザ・ボーン（THE BONE）1991.9 vol5 No3 p73-79）。
し、それをもとに補正領域２１の設定，補正値の算出，
計測領域の抽出，各部の面積，骨塩量，骨密度を計測し
ている。計測領域は、図２に示すように大腿骨中枢端部
の、頸部２２，大転子部２３，ウォード三角部２４（以
下それぞれＮ部，Ｔ部，Ｗ部と表わす）の３ヵ所である
（ザ・ボーン（THE BONE）1991.9 vol5 No3 p73-79）。
【０００７】
【発明が解決しようとする課題】従来は、計測対象領域
２０を操作者が指示して決定するため、これに左右され
る補正領域２１の設定再現性が高くならなかった。この
ため、計測領域が再現性高く設定されても、計測値の再
現性につながらなかった。
２０を操作者が指示して決定するため、これに左右され
る補正領域２１の設定再現性が高くならなかった。この
ため、計測領域が再現性高く設定されても、計測値の再
現性につながらなかった。
【０００８】本発明の目的は、大腿骨中枢端部の骨密度
計測において、計測対象領域の設定と計測領域の抽出の
両方の処理を、大腿骨の解剖学的形状に基づいて、高い
再現性で行う骨密度の計測方法を提供することにある。
計測において、計測対象領域の設定と計測領域の抽出の
両方の処理を、大腿骨の解剖学的形状に基づいて、高い
再現性で行う骨密度の計測方法を提供することにある。
【０００９】
【課題を解決するための手段】本発明は上記目的を達成
するため、図５に示すように、大腿骨中枢端部の骨密度
計測時に通常撮影する、左右いずれかの骨盤から大腿骨
を含む部分を高低２種類のエネルギで撮影した２枚のＸ
線画像５０，５１から、予め定めた係数を用いて仮の差
分画像５２を作成する。次に、大腿骨中枢端部の解剖学
的形状と、その形状に伴う前記仮の差分画像５２の縦お
よび横方向の濃度値変化をもとに、大腿骨中枢端部を囲
む計測対象領域２０を設定し、前記計測対象領域２０内
の軟部組織部分の領域である補正領域２１の濃度値から
補正値Ｆを算出し、補正値Ｆと前記仮の差分画像５２よ
り骨密度を計測するための差分画像５４を作成する。更
に前記差分画像５４から、大腿骨のシャフトの主軸を算
出し、前記主軸を基準として大腿骨中枢端部の計測領域
すなわち前記Ｎ部，Ｔ部，Ｗ部を抽出し、各部の領域内
の差分画像５４の濃度値の和をもとに骨塩量を算出し、
平均濃度値をもとに骨密度を算出する。
するため、図５に示すように、大腿骨中枢端部の骨密度
計測時に通常撮影する、左右いずれかの骨盤から大腿骨
を含む部分を高低２種類のエネルギで撮影した２枚のＸ
線画像５０，５１から、予め定めた係数を用いて仮の差
分画像５２を作成する。次に、大腿骨中枢端部の解剖学
的形状と、その形状に伴う前記仮の差分画像５２の縦お
よび横方向の濃度値変化をもとに、大腿骨中枢端部を囲
む計測対象領域２０を設定し、前記計測対象領域２０内
の軟部組織部分の領域である補正領域２１の濃度値から
補正値Ｆを算出し、補正値Ｆと前記仮の差分画像５２よ
り骨密度を計測するための差分画像５４を作成する。更
に前記差分画像５４から、大腿骨のシャフトの主軸を算
出し、前記主軸を基準として大腿骨中枢端部の計測領域
すなわち前記Ｎ部，Ｔ部，Ｗ部を抽出し、各部の領域内
の差分画像５４の濃度値の和をもとに骨塩量を算出し、
平均濃度値をもとに骨密度を算出する。
【００１０】他の手段として、差分画像５４を作成せず
に、仮の差分画像上５２で、前記Ｎ部，Ｔ部，Ｗ部を抽
出し、各部の領域内の仮の差分画像５２の濃度値の和と
補正値Ｆより骨塩量を算出し、平均濃度値と補正値Ｆよ
り骨密度を算出することもできる。更に、骨盤と大腿骨
の連結部からＮ部、またはＴ部までの長さ，計測対象領
域内に含まれる大腿骨中枢端部全体の骨塩量，骨密度等
も算出することができる。
に、仮の差分画像上５２で、前記Ｎ部，Ｔ部，Ｗ部を抽
出し、各部の領域内の仮の差分画像５２の濃度値の和と
補正値Ｆより骨塩量を算出し、平均濃度値と補正値Ｆよ
り骨密度を算出することもできる。更に、骨盤と大腿骨
の連結部からＮ部、またはＴ部までの長さ，計測対象領
域内に含まれる大腿骨中枢端部全体の骨塩量，骨密度等
も算出することができる。
【００１１】
【作用】大腿骨周辺は、骨折などの著しい変化、または
撮影姿勢に大きな差がない限り、同一患者ならば撮影画
像における大腿骨部の解剖学的形状はほぼ同一と考えら
れる。上記手段は、この解剖学的形状をもとに計測対象
領域２０を設定するため、同一患者に対しては常に高い
再現性でこの領域を設定することができる。例えば、撮
影位置が縦方向または横方向にずれている場合でも、少
なくとも計測領域が撮影されているならば、解剖学的形
状をもとに画像のずれに影響されることなく計測対象領
域２０を設定することが可能である。
撮影姿勢に大きな差がない限り、同一患者ならば撮影画
像における大腿骨部の解剖学的形状はほぼ同一と考えら
れる。上記手段は、この解剖学的形状をもとに計測対象
領域２０を設定するため、同一患者に対しては常に高い
再現性でこの領域を設定することができる。例えば、撮
影位置が縦方向または横方向にずれている場合でも、少
なくとも計測領域が撮影されているならば、解剖学的形
状をもとに画像のずれに影響されることなく計測対象領
域２０を設定することが可能である。
【００１２】本発明は、図５のように、撮影したＬ画像
５０，Ｈ画像５１より作成する仮の差分画像５２と、骨
と軟部組織を区別するしきい濃度値ｔを用いて、仮の差
分画像５２の縦および横方向の濃度値変化をもとに、図
２に示すような計測対象領域２０を設定する。
５０，Ｈ画像５１より作成する仮の差分画像５２と、骨
と軟部組織を区別するしきい濃度値ｔを用いて、仮の差
分画像５２の縦および横方向の濃度値変化をもとに、図
２に示すような計測対象領域２０を設定する。
【００１３】例えば図２において、大腿骨と軟部組織の
縦方向の境界２０１は、しきい濃度値ｔ未満を０、以上
を１とすれば、骨部分は１、軟部組織や空気の部分の濃
度は０になるので、骨のない部分の画像の縦方向の累積
値は０であり、骨の部分を含むと累積値は０でなくな
る。したがって、縦方向の累積値が０である位置と０で
ない位置との境界を、境界２０１として設定できる。ま
た、骨盤と大腿骨の連結部分の横方向の境界２０２は、
仮の差分画像５２上では、連結部分の骨が薄いために少
し低濃度になっており、連結部分の上部の骨盤は、骨の
厚さが厚いため周囲に比べて高濃度になっている。した
がって、仮の差分画像５２における横方向の累積濃度値
または平均濃度値が最大になり、その少し下の位置で累
積濃度値が減少→増加になっている位置を、境界２０２
として設定できる。仮の差分画像５２の下端部の骨はシ
ャフトであるので、シャフトの左右の端点ｐ′２５，
ｑ′２６を抽出することができる。この端点を抽出する
ことにより、大腿骨中枢端部の位置や、大まかな大きさ
を決定でき、これに基づいて上記の境界２０１〜２０４
を設定する際の累積値を算出する範囲を決定できる。以
上の方法で設定した計測対象領域２０内で、しきい濃度
値ｔ未満の領域を補正領域２１と設定する。計測領域
は、図２に示すＮ部２２，Ｔ部２３，Ｗ部２４である。
縦方向の境界２０１は、しきい濃度値ｔ未満を０、以上
を１とすれば、骨部分は１、軟部組織や空気の部分の濃
度は０になるので、骨のない部分の画像の縦方向の累積
値は０であり、骨の部分を含むと累積値は０でなくな
る。したがって、縦方向の累積値が０である位置と０で
ない位置との境界を、境界２０１として設定できる。ま
た、骨盤と大腿骨の連結部分の横方向の境界２０２は、
仮の差分画像５２上では、連結部分の骨が薄いために少
し低濃度になっており、連結部分の上部の骨盤は、骨の
厚さが厚いため周囲に比べて高濃度になっている。した
がって、仮の差分画像５２における横方向の累積濃度値
または平均濃度値が最大になり、その少し下の位置で累
積濃度値が減少→増加になっている位置を、境界２０２
として設定できる。仮の差分画像５２の下端部の骨はシ
ャフトであるので、シャフトの左右の端点ｐ′２５，
ｑ′２６を抽出することができる。この端点を抽出する
ことにより、大腿骨中枢端部の位置や、大まかな大きさ
を決定でき、これに基づいて上記の境界２０１〜２０４
を設定する際の累積値を算出する範囲を決定できる。以
上の方法で設定した計測対象領域２０内で、しきい濃度
値ｔ未満の領域を補正領域２１と設定する。計測領域
は、図２に示すＮ部２２，Ｔ部２３，Ｗ部２４である。
【００１４】Ｎ部２２は、図３に示す大腿骨頸部の中心
軸ＣＤ３２に垂直、かつ最もくびれている領域である。
従って、中心軸ＣＤ３２の算出がＮ部２２の抽出に必要
である。しかし、大腿骨頸部は長さも短く、輪郭の形状
も滑らかとは限らないため、直接頸部の輪郭より中心軸
ＣＤ３２を算出するのは精度の面で問題がある。そこ
で、本発明では、図３に示す大腿骨のシャフトの中心軸
ＡＢ３１をもとに頸部の中心軸ＣＤ３２を算出する。
軸ＣＤ３２に垂直、かつ最もくびれている領域である。
従って、中心軸ＣＤ３２の算出がＮ部２２の抽出に必要
である。しかし、大腿骨頸部は長さも短く、輪郭の形状
も滑らかとは限らないため、直接頸部の輪郭より中心軸
ＣＤ３２を算出するのは精度の面で問題がある。そこ
で、本発明では、図３に示す大腿骨のシャフトの中心軸
ＡＢ３１をもとに頸部の中心軸ＣＤ３２を算出する。
【００１５】ラジオロジー（RADIOLOGY）87：p904−90
7，November 1966 によれば、中心軸ＡＢ３１とＣＤ３
２のなす角α３３は約１２７°である。中心軸ＡＢ３１
は、輪郭も十分長く、かつ輪郭の形状も滑らかであるた
め精度良く算出できる。このように、まずＡＢ３１を算
出し、ＡＢ３１と１２７°をなすＣＤ′を設定した後、
ＣＤ′に補正を加えてＣＤ３２を算出することで、ＣＤ
３２を精度良く設定することができる。
7，November 1966 によれば、中心軸ＡＢ３１とＣＤ３
２のなす角α３３は約１２７°である。中心軸ＡＢ３１
は、輪郭も十分長く、かつ輪郭の形状も滑らかであるた
め精度良く算出できる。このように、まずＡＢ３１を算
出し、ＡＢ３１と１２７°をなすＣＤ′を設定した後、
ＣＤ′に補正を加えてＣＤ３２を算出することで、ＣＤ
３２を精度良く設定することができる。
【００１６】Ｔ部２３は、Ｎ部２２と隣接し、シャフト
にかからない領域として抽出する。Ｔ部２３とシャフト
との境界は明確に設定することができないが、頸部の中
心軸ＣＤ３２より上部に境界を設定するのが適当である
ことから、頸部の中心軸CD32に対して予め定めた角β７
１をなす傾きで、Ｎ部２２の重心Ｇを通る直線を境界と
する。または、シャフトの中心軸ＡＢ３１に対して予め
定めた角β′７２をなす傾きで、Ｎ部２２の重心を通る
直線を境界とする。上述のとおりシャフトの中心軸ＡＢ
３１および頸部の中心軸ＣＤ３２は、精度よく算出され
るので、これに基づいて設定する境界も精度よく設定す
ることができる。
にかからない領域として抽出する。Ｔ部２３とシャフト
との境界は明確に設定することができないが、頸部の中
心軸ＣＤ３２より上部に境界を設定するのが適当である
ことから、頸部の中心軸CD32に対して予め定めた角β７
１をなす傾きで、Ｎ部２２の重心Ｇを通る直線を境界と
する。または、シャフトの中心軸ＡＢ３１に対して予め
定めた角β′７２をなす傾きで、Ｎ部２２の重心を通る
直線を境界とする。上述のとおりシャフトの中心軸ＡＢ
３１および頸部の中心軸ＣＤ３２は、精度よく算出され
るので、これに基づいて設定する境界も精度よく設定す
ることができる。
【００１７】Ｗ部２４は、皮質骨が多いＮ部２２とＴ部
２３の隣接部周辺にあり、海綿骨が多く、骨梁の狭間に
あるため比較的低濃度になっている。このことからＮ部
２２とＴ部２３の隣接部周辺で最も平均濃度が低い位置
をＷ部２４と設定する。Ｗ部２４の面積は予め定めてお
き、形状はＮ部の傾きと中心軸ＣＤ３２の傾きより設定
する。以上のように、精度良く設定できるシャフトの中
心軸ＡＢ３１や頸部の中心軸ＣＤ３２に基づいて各計測
領域を設定していくので領域の設定再現性も高くでき
る。
２３の隣接部周辺にあり、海綿骨が多く、骨梁の狭間に
あるため比較的低濃度になっている。このことからＮ部
２２とＴ部２３の隣接部周辺で最も平均濃度が低い位置
をＷ部２４と設定する。Ｗ部２４の面積は予め定めてお
き、形状はＮ部の傾きと中心軸ＣＤ３２の傾きより設定
する。以上のように、精度良く設定できるシャフトの中
心軸ＡＢ３１や頸部の中心軸ＣＤ３２に基づいて各計測
領域を設定していくので領域の設定再現性も高くでき
る。
【００１８】
【実施例】本実施例は、ＤＥＸＡ法を用いて、大腿骨中
枢端部の骨塩量および骨密度を計測する方法である。
枢端部の骨塩量および骨密度を計測する方法である。
【００１９】図４は本発明を実現するシステム構成の一
例である。以下、各部の機能概要を説明する。Ｘ線発生
器４１から高低２種類の電圧でファンビーム状のＸ線を
交互に発生する。Ｘ線は被写体４０を通過し、Ｘ線検出
器４２で検出され、それぞれＡＤ変換器４３で変換され
る。発生器４１と検出器４２は連動して撮影区間を移動
しながら撮影を行い、最終的に得られる２種類の画像が
画像処理装置４４に送られる。画像処理装置４４では、
以下に述べる骨塩量および骨密度計測処理等を行う。表
示装置４５は画像や操作メッセージ，計測結果等を表示
し、操作卓４６はキーボード，マウス，トラックボール
などを備えており指示，入力を行う。外部記憶装置４７
は画像データ等を記憶し、出力装置４８は表示装置４５
の表示内容のハードコピーを行う。
例である。以下、各部の機能概要を説明する。Ｘ線発生
器４１から高低２種類の電圧でファンビーム状のＸ線を
交互に発生する。Ｘ線は被写体４０を通過し、Ｘ線検出
器４２で検出され、それぞれＡＤ変換器４３で変換され
る。発生器４１と検出器４２は連動して撮影区間を移動
しながら撮影を行い、最終的に得られる２種類の画像が
画像処理装置４４に送られる。画像処理装置４４では、
以下に述べる骨塩量および骨密度計測処理等を行う。表
示装置４５は画像や操作メッセージ，計測結果等を表示
し、操作卓４６はキーボード，マウス，トラックボール
などを備えており指示，入力を行う。外部記憶装置４７
は画像データ等を記憶し、出力装置４８は表示装置４５
の表示内容のハードコピーを行う。
【００２０】以下、図１の処理手順のフローチャートに
従って大腿骨中枢端部の骨塩量および骨密度の計測方法
について説明する。
従って大腿骨中枢端部の骨塩量および骨密度の計測方法
について説明する。
【００２１】ステップ１０：仮の差分画像Ｓ′５２の作
成（図５） 高低２種類の電圧で撮影した画像（以下、Ｌ画像５０，
Ｈ画像５１）、と定数ｒ，Ｈ画像５１の濃度値関数とし
て定義される計測位置ごとの厚みに対する補正値Ｒ(Ｈ
(ｉ，ｊ)）を用いて、仮の差分画像Ｓ′５２をＳ′
(ｉ，ｊ)＝（Ｌ′(ｉ，ｊ)−ｒ・Ｈ′(ｉ，ｊ)）・Ｒ
（Ｈ′(ｉ，ｊ)）として作成する。ここで、撮影方法に
よってＬ画像５０とＨ画像５１は理論上０.５ 画素ずれ
ているため、少なくとも一方の画像をキュービック コ
ンボリューションCubic Convolution 法などで位置合わ
せした画像Ｌ′，Ｈ′を用いる。また、画像Ｌ′に、ビ
ームハードニング現象が起こっている場合には、その補
正も行う。ビームハードニングの補正は、補正値Ｒと同
様にＨ画像５１の関数として求めることができる。
成（図５） 高低２種類の電圧で撮影した画像（以下、Ｌ画像５０，
Ｈ画像５１）、と定数ｒ，Ｈ画像５１の濃度値関数とし
て定義される計測位置ごとの厚みに対する補正値Ｒ(Ｈ
(ｉ，ｊ)）を用いて、仮の差分画像Ｓ′５２をＳ′
(ｉ，ｊ)＝（Ｌ′(ｉ，ｊ)−ｒ・Ｈ′(ｉ，ｊ)）・Ｒ
（Ｈ′(ｉ，ｊ)）として作成する。ここで、撮影方法に
よってＬ画像５０とＨ画像５１は理論上０.５ 画素ずれ
ているため、少なくとも一方の画像をキュービック コ
ンボリューションCubic Convolution 法などで位置合わ
せした画像Ｌ′，Ｈ′を用いる。また、画像Ｌ′に、ビ
ームハードニング現象が起こっている場合には、その補
正も行う。ビームハードニングの補正は、補正値Ｒと同
様にＨ画像５１の関数として求めることができる。
【００２２】ＤＥＸＡ法で撮影する画像は積分像であ
り、骨密度計測のために作成する差分画像Ｓ５４は、軟
部組織の濃度を０とした画像で、骨の部分の画素濃度は
０より大きく、画素濃度がそのまま骨密度になるもので
ある。脂肪がない場合は、軟部組織は水とほとんど変わ
らないため、ここで作成したＳ′５２を差分画像Ｓ５４
として用いることができる。しかし、実際は脂肪による
影響が生じる。脂肪量が多いと、画素濃度は低下してし
まうため、各データごとに脂肪量に基づいた補正値Ｆを
求め、Ｓ′５２を補正する必要がある。
り、骨密度計測のために作成する差分画像Ｓ５４は、軟
部組織の濃度を０とした画像で、骨の部分の画素濃度は
０より大きく、画素濃度がそのまま骨密度になるもので
ある。脂肪がない場合は、軟部組織は水とほとんど変わ
らないため、ここで作成したＳ′５２を差分画像Ｓ５４
として用いることができる。しかし、実際は脂肪による
影響が生じる。脂肪量が多いと、画素濃度は低下してし
まうため、各データごとに脂肪量に基づいた補正値Ｆを
求め、Ｓ′５２を補正する必要がある。
【００２３】ステップ１１：計測対象領域２０と補正領
域２１の設定（図２および図５） 補正値Ｆは、計測領域の周辺の軟部組織の平均濃度を０
にする値として算出する。補正値Ｆを求めるための補正
領域２１を設定する。補正領域２１の設定は、計測対象
領域２０を設定する必要がある。まず、Ｓ′５２より軟
部組織と骨の部分を区別する２値画像５３を、例えば
Ｓ′の画素濃度が予め定めたしきい濃度値ｔ以上であれ
ば１、ｔ未満ならば０として作成する。しきい濃度値ｔ
は０を用いると、ほぼ軟部組織と骨の部分を区別でき
る。しかしこの処理だけでは軟部組織部分や骨にごま塩
状のノイズが発生するためノイズ除去を行って最終の２
値画像５３とする。ノイズ除去の方法は、平滑化や注目
点の周辺の濃度値から注目点が孤立点であるか判定して
除去する方法などがある。
域２１の設定（図２および図５） 補正値Ｆは、計測領域の周辺の軟部組織の平均濃度を０
にする値として算出する。補正値Ｆを求めるための補正
領域２１を設定する。補正領域２１の設定は、計測対象
領域２０を設定する必要がある。まず、Ｓ′５２より軟
部組織と骨の部分を区別する２値画像５３を、例えば
Ｓ′の画素濃度が予め定めたしきい濃度値ｔ以上であれ
ば１、ｔ未満ならば０として作成する。しきい濃度値ｔ
は０を用いると、ほぼ軟部組織と骨の部分を区別でき
る。しかしこの処理だけでは軟部組織部分や骨にごま塩
状のノイズが発生するためノイズ除去を行って最終の２
値画像５３とする。ノイズ除去の方法は、平滑化や注目
点の周辺の濃度値から注目点が孤立点であるか判定して
除去する方法などがある。
【００２４】次に、Ｓ′５２と２値画像５３を用いて計
測対象領域２０の各境界２０１〜２０４を設定する。境
界２０１は２値画像５３の縦方向の累積濃度値の変化か
ら定め、境界２０２はＳ′５２の横方向の累積濃度値の
変化から定めることができる。同様に境界２０３は、
Ｓ′５２の縦方向の累積濃度値から定める。これらの境
界２０１〜２０３は累積濃度値から定めるが、その累積
する範囲は撮影位置の下方の横のラインからシャフトの
左右端点ｐ′２５，ｑ′２６を求め、その端点の位置か
ら大腿骨の解剖学的形状に基づいて、累積範囲を設定す
る。また、２値画像５３より大腿骨の解剖学的形状に基
づいて、変曲点ｐ2〜ｐ7，ｑ2，ｑ3も求めることができ
るので、それらの位置も用いて累積範囲を設定すること
ができる。
測対象領域２０の各境界２０１〜２０４を設定する。境
界２０１は２値画像５３の縦方向の累積濃度値の変化か
ら定め、境界２０２はＳ′５２の横方向の累積濃度値の
変化から定めることができる。同様に境界２０３は、
Ｓ′５２の縦方向の累積濃度値から定める。これらの境
界２０１〜２０３は累積濃度値から定めるが、その累積
する範囲は撮影位置の下方の横のラインからシャフトの
左右端点ｐ′２５，ｑ′２６を求め、その端点の位置か
ら大腿骨の解剖学的形状に基づいて、累積範囲を設定す
る。また、２値画像５３より大腿骨の解剖学的形状に基
づいて、変曲点ｐ2〜ｐ7，ｑ2，ｑ3も求めることができ
るので、それらの位置も用いて累積範囲を設定すること
ができる。
【００２５】境界２０４は、簡単な方法としては撮影位
置の一番下方に設定したり、境界２０２から予め定めた
距離だけはなれた位置に設定できる。しかし、この方法
は計測対象領域２０が絶対的な大きさになってしまうた
め、撮影条件の違いによって計測対象領域２０の相対的
な大きさ，設定位置が異なってしまう。そこで、撮影条
件によらず相対的に計測対象領域２０を設定する方法と
して、例えばｐ2 を検出し、ｐ2 と境界２０２間の距離
を２倍した位置に境界２０４を設定する。さらに、計測
対象領域２０内の２値画像５３の濃度値が０である領域
を補正領域２１とする。
置の一番下方に設定したり、境界２０２から予め定めた
距離だけはなれた位置に設定できる。しかし、この方法
は計測対象領域２０が絶対的な大きさになってしまうた
め、撮影条件の違いによって計測対象領域２０の相対的
な大きさ，設定位置が異なってしまう。そこで、撮影条
件によらず相対的に計測対象領域２０を設定する方法と
して、例えばｐ2 を検出し、ｐ2 と境界２０２間の距離
を２倍した位置に境界２０４を設定する。さらに、計測
対象領域２０内の２値画像５３の濃度値が０である領域
を補正領域２１とする。
【００２６】以上の方法で、撮影条件によって被写体の
位置ずれが発生していたり、大きさが異なっていても相
対的に一定の計測対象領域２０および補正領域２１が設
定できる。他の方法として、シャフトの左右端点ｐ′２
５，ｑ′２６，変曲点ｐ2〜ｐ7，ｑ2，ｑ3 から、直接
境界２０１〜２０４を設定することもできる。例えば、
境界２０１は変曲点ｐ2，ｐ3の位置、境界２０２は変曲
点ｐ4，ｐ5の位置からｐ2とｐ4間の縦方向の長さ分上の
位置、境界２０３は変曲点ｑ3 の位置、を基準に設定
し、境界２０４は上記の方法で設定する。
位置ずれが発生していたり、大きさが異なっていても相
対的に一定の計測対象領域２０および補正領域２１が設
定できる。他の方法として、シャフトの左右端点ｐ′２
５，ｑ′２６，変曲点ｐ2〜ｐ7，ｑ2，ｑ3 から、直接
境界２０１〜２０４を設定することもできる。例えば、
境界２０１は変曲点ｐ2，ｐ3の位置、境界２０２は変曲
点ｐ4，ｐ5の位置からｐ2とｐ4間の縦方向の長さ分上の
位置、境界２０３は変曲点ｑ3 の位置、を基準に設定
し、境界２０４は上記の方法で設定する。
【００２７】ステップ１２：補正値Ｆ算出 Ｓ′５２において、補正領域２１の平均濃度値Ｆを求め
る。
る。
【００２８】ステップ１３：差分画像Ｓ５４作成 Ｓ′５２と補正値Ｆより脂肪による濃度値のずれを補正
し、骨のある位置の画素濃度がその位置の骨密度になる
ような差分画像Ｓ５４を作成する。
し、骨のある位置の画素濃度がその位置の骨密度になる
ような差分画像Ｓ５４を作成する。
【００２９】Ｓ(ｉ，ｊ)＝Ｓ′(ｉ，ｊ)−Ｆ ステップ１４：計測領域抽出 ２値画像５３を用いて計測対象領域２０の境界と大腿骨
の輪郭の接する点ｐ1,ｐ8，ｑ1，ｑ4 を求める。以下、
大腿骨中枢端部のＮ部２２，Ｔ部２３，Ｗ部２４の抽出
方法を順に説明する。
の輪郭の接する点ｐ1,ｐ8，ｑ1，ｑ4 を求める。以下、
大腿骨中枢端部のＮ部２２，Ｔ部２３，Ｗ部２４の抽出
方法を順に説明する。
【００３０】（１）Ｎ部２２の抽出（図６） （ａ）シャフトの中心軸ＡＢ３１の設定 大腿骨を囲む輪郭のｐ1とｐ2の中点ｐx、ｑ1とｑ2の中
点ｑxを求め、ｐ1とｐx間の輪郭座標と、ｑ1とｑx間の
輪郭座標を用いて、それぞれの近似直線を求め、その中
線を中心軸ＡＢ３１とする。他の方法として、領域ｐ1
ｐxｑxｑ1内の全座標から慣性の主軸を計算し、それを
中心軸ＡＢ３１とすることもできる。また、ｐx はｐ1
とｐ2の間の中点とは限らず、もちろんｐ1とｐ2間のど
こに設定してもかまわない。ｑxについても同様であ
る。しかし、ｐx がｐ2，ｑxがｑ2に近づくにつれてシ
ャフトが広がってくるので、シャフトが安定している位
置にｐxｑxを設定した方が良い。
点ｑxを求め、ｐ1とｐx間の輪郭座標と、ｑ1とｑx間の
輪郭座標を用いて、それぞれの近似直線を求め、その中
線を中心軸ＡＢ３１とする。他の方法として、領域ｐ1
ｐxｑxｑ1内の全座標から慣性の主軸を計算し、それを
中心軸ＡＢ３１とすることもできる。また、ｐx はｐ1
とｐ2の間の中点とは限らず、もちろんｐ1とｐ2間のど
こに設定してもかまわない。ｑxについても同様であ
る。しかし、ｐx がｐ2，ｑxがｑ2に近づくにつれてシ
ャフトが広がってくるので、シャフトが安定している位
置にｐxｑxを設定した方が良い。
【００３１】（ｂ）大腿骨頸部の中心軸ＣＤ３２の設定
とＮ部２２の抽出 中心軸ＡＢ３１と反時計回りに角α３３をなす直線の傾
きを、求める中心軸ＣＤ３２の初期傾きａ1 とする。角
α３３は、１２５°〜１３０°であり、ここでは１２７
°にする。次に、頸部の輪郭ｐ6〜ｐ7とｑ2〜ｑ3間で最
もくびれている対座標ＰN、ＱNを検出し、線分ＰNＱNに
垂直な傾きａ2 求め、ａ1とａ2の中間の傾きａを中心軸
ＣＤ３２の傾きとする。線分ＰNＱNの中点を中心とし、
傾きａに垂直な傾きｂで中心から予め定めた距離ｋ離れ
た線分Ｎ1Ｎ2、Ｎ3Ｎ4を設定し、この線分と輪郭で囲ま
れる領域をＮ部２２とする。更に、Ｎ部２２の設定位置
が不適切である場合は傾きａに垂直な方向に平行移動し
てＮ部の位置を補正する。Ｎ部２２の設定位置が不適切
とは、例えばＮ4 がｐ6よりｐ5側、またはＮ1 がｐ7よ
りｐ8側、またはＮ2がｑ3よりｑ4 側に設定されるよう
な場合である。補正したＮ部２２の重心Ｇを求め、傾き
ａでＧを通る直線を中心軸ＣＤ３２とする。
とＮ部２２の抽出 中心軸ＡＢ３１と反時計回りに角α３３をなす直線の傾
きを、求める中心軸ＣＤ３２の初期傾きａ1 とする。角
α３３は、１２５°〜１３０°であり、ここでは１２７
°にする。次に、頸部の輪郭ｐ6〜ｐ7とｑ2〜ｑ3間で最
もくびれている対座標ＰN、ＱNを検出し、線分ＰNＱNに
垂直な傾きａ2 求め、ａ1とａ2の中間の傾きａを中心軸
ＣＤ３２の傾きとする。線分ＰNＱNの中点を中心とし、
傾きａに垂直な傾きｂで中心から予め定めた距離ｋ離れ
た線分Ｎ1Ｎ2、Ｎ3Ｎ4を設定し、この線分と輪郭で囲ま
れる領域をＮ部２２とする。更に、Ｎ部２２の設定位置
が不適切である場合は傾きａに垂直な方向に平行移動し
てＮ部の位置を補正する。Ｎ部２２の設定位置が不適切
とは、例えばＮ4 がｐ6よりｐ5側、またはＮ1 がｐ7よ
りｐ8側、またはＮ2がｑ3よりｑ4 側に設定されるよう
な場合である。補正したＮ部２２の重心Ｇを求め、傾き
ａでＧを通る直線を中心軸ＣＤ３２とする。
【００３２】（２）Ｔ部２３の抽出（図７） 中心軸ＣＤ３２と線分Ｎ3Ｎ4の交点をＴ2とする。Ｔ2を
通り中心軸ＣＤ３２と時計回りに角β７１をなす直線の
大腿骨輪郭との交点をＴ1とし、Ｔ1Ｔ2Ｎ4とＴ1とＮ4間
の輪郭で囲まれる領域をＴ２３とする。角β７１は５°
〜２０°程度であり、ここでは１０°にする。別の方法
として、Ｔ1をＣとｐ2の中間輪郭位置にしたり、Ｔ2を
線分Ｔ1Ｇと線分Ｎ3Ｎ4の交点とすることもできる。更
に他の方法として、シャフトの中心軸ＡＢ３１と時計回
りに角β′７２をなす傾きで、Ｔ2を通る直線と大腿骨
輪郭との交点をＴ1とし、Ｔ部２３を抽出する。角β′7
2は６０°〜７０°程度である。
通り中心軸ＣＤ３２と時計回りに角β７１をなす直線の
大腿骨輪郭との交点をＴ1とし、Ｔ1Ｔ2Ｎ4とＴ1とＮ4間
の輪郭で囲まれる領域をＴ２３とする。角β７１は５°
〜２０°程度であり、ここでは１０°にする。別の方法
として、Ｔ1をＣとｐ2の中間輪郭位置にしたり、Ｔ2を
線分Ｔ1Ｇと線分Ｎ3Ｎ4の交点とすることもできる。更
に他の方法として、シャフトの中心軸ＡＢ３１と時計回
りに角β′７２をなす傾きで、Ｔ2を通る直線と大腿骨
輪郭との交点をＴ1とし、Ｔ部２３を抽出する。角β′7
2は６０°〜７０°程度である。
【００３３】（３）Ｗ部２４の抽出 Ｔ2を中心とし、半径がＴ2Ｇの長さの円形領域内におい
て、周辺の何点かで算出した差分画像５４の平均濃度が
最低になる位置Ｗ1を求め、Ｗ1を中心に予め定めた面積
Ｚになるように矩形領域Ｗ部２４を抽出する。Ｗ部２４
の各辺の傾きは中心軸ＣＤ３２および線分Ｎ3Ｎ4に平行
になるように設定する。ここでＷ1 を求める領域は円形
でなくとも構わないし、Ｗ1 は平均濃度で求めずに画素
濃度でもよい。さらに、Ｗ部２４は矩形でなく円形でも
よい。
て、周辺の何点かで算出した差分画像５４の平均濃度が
最低になる位置Ｗ1を求め、Ｗ1を中心に予め定めた面積
Ｚになるように矩形領域Ｗ部２４を抽出する。Ｗ部２４
の各辺の傾きは中心軸ＣＤ３２および線分Ｎ3Ｎ4に平行
になるように設定する。ここでＷ1 を求める領域は円形
でなくとも構わないし、Ｗ1 は平均濃度で求めずに画素
濃度でもよい。さらに、Ｗ部２４は矩形でなく円形でも
よい。
【００３４】ステップ１５：骨塩量，骨密度算出 上記で抽出したＮ部２２，Ｔ部２３，Ｗ部２４の差分画
像５４における各領域内の総座標数を面積とし、累積濃
度値を骨塩量、平均濃度値を骨密度とする。さらに、輪
郭ｐ1〜ｐ5〜Ｎ1と線分Ｎ1Ｎ2と輪郭ｑ1〜ｑ2〜Ｎ2と線
分ｐ1ｑ1で囲まれる領域を大腿骨中枢端部全体の領域と
して、面積，骨塩量，骨密度を算出する。または、輪郭
ｐ1〜ｐ8と線分ｐ8ｑ4と輪郭ｑ1〜ｑ4と線分ｐ1ｑ1で囲
まれる領域を大腿骨中枢端部全体の領域としてもよい。
また、他の測定値として線分CD32や、線分Ｔ2Ｄ の長さ
も算出し、骨密度と共に骨折の危険性を診断する指標と
して用いることもできる。
像５４における各領域内の総座標数を面積とし、累積濃
度値を骨塩量、平均濃度値を骨密度とする。さらに、輪
郭ｐ1〜ｐ5〜Ｎ1と線分Ｎ1Ｎ2と輪郭ｑ1〜ｑ2〜Ｎ2と線
分ｐ1ｑ1で囲まれる領域を大腿骨中枢端部全体の領域と
して、面積，骨塩量，骨密度を算出する。または、輪郭
ｐ1〜ｐ8と線分ｐ8ｑ4と輪郭ｑ1〜ｑ4と線分ｐ1ｑ1で囲
まれる領域を大腿骨中枢端部全体の領域としてもよい。
また、他の測定値として線分CD32や、線分Ｔ2Ｄ の長さ
も算出し、骨密度と共に骨折の危険性を診断する指標と
して用いることもできる。
【００３５】図１のフローチャートに従って、大腿骨中
枢端部の骨密度を測定する方法の他に、差分画像５４を
作成せずに、仮の差分画像５２と２値画像５３より計測
領域を抽出し、骨塩量を仮の差分画像５２の累積濃度値
と補正値Ｆより求めることができる。例えば、Ｎ部２２
の仮の差分画像５２における累積濃度値をＢＭＣＮ′と
すると、Ｎ部２２の骨塩量ＢＭＣＮは、ＢＭＣＮ＝ＢＭ
ＣＮ′−Ｆとなる。骨密度も同様に、Ｎ部２２の仮の差
分画像５２における平均濃度値をＢＭＤＮ′とすると、
Ｎ部２２の骨密度ＢＭＤＮは、ＢＭＤＮ＝ＢＭＤＮ′−
Ｆとなる。他の領域も同様にして求めることができる。
枢端部の骨密度を測定する方法の他に、差分画像５４を
作成せずに、仮の差分画像５２と２値画像５３より計測
領域を抽出し、骨塩量を仮の差分画像５２の累積濃度値
と補正値Ｆより求めることができる。例えば、Ｎ部２２
の仮の差分画像５２における累積濃度値をＢＭＣＮ′と
すると、Ｎ部２２の骨塩量ＢＭＣＮは、ＢＭＣＮ＝ＢＭ
ＣＮ′−Ｆとなる。骨密度も同様に、Ｎ部２２の仮の差
分画像５２における平均濃度値をＢＭＤＮ′とすると、
Ｎ部２２の骨密度ＢＭＤＮは、ＢＭＤＮ＝ＢＭＤＮ′−
Ｆとなる。他の領域も同様にして求めることができる。
【００３６】また、実施例では２値画像５３を作成する
が、仮の差分画像５２と、骨と軟部組織とを区別するし
きい濃度値ｔのみを用いて処理することもできるため、
２値画像５３を作成しなくてもよい。
が、仮の差分画像５２と、骨と軟部組織とを区別するし
きい濃度値ｔのみを用いて処理することもできるため、
２値画像５３を作成しなくてもよい。
【００３７】
【発明の効果】本発明は、大腿骨中枢端部の解剖学的形
状に基づいて、計測対象領域の設定と、再現性の高い補
正値の算出と計測領域の抽出が自動的にできるため、骨
密度の計測方法として有効である。
状に基づいて、計測対象領域の設定と、再現性の高い補
正値の算出と計測領域の抽出が自動的にできるため、骨
密度の計測方法として有効である。
【図１】本発明の一実施例の骨密度計測の処理フローチ
ャート。
ャート。
【図２】本発明の一実施例における画像処理の説明図。
【図３】大腿骨のシャフトの主軸と大腿骨頸部の主軸の
説明図。
説明図。
【図４】本発明を実現するシステム構成の一例を示す説
明図。
明図。
【図５】本発明の一実施例の差分画像作成の説明図。
【図６】Ｎ部２２の抽出方法の説明図。
【図７】Ｔ部２３の抽出方法の説明図。
２０…計測対象領域、２２…Ｎ部、２３…Ｔ部、２５，
２６…シャフトの左右端点ｐ′、３１…シャフトの主
軸、３２，３３…大腿骨頸部の主軸、５０…Ｌ画像、５
１…Ｈ画像、５４…差分画像。
２６…シャフトの左右端点ｐ′、３１…シャフトの主
軸、３２，３３…大腿骨頸部の主軸、５０…Ｌ画像、５
１…Ｈ画像、５４…差分画像。
───────────────────────────────────────────────────── フロントページの続き (72)発明者 熊本 三矢戒 東京都千代田区内神田一丁目１番14号 株 式会社日立メディコ内
Claims (9)
- 【請求項１】高低２種類のエネルギで撮影した大腿骨中
枢端部を含むＸ線画像より差分画像を作成し、前記差分
画像に基づいて大腿骨頸部の骨密度を計測する方法にお
いて、前記大腿骨中枢端部を囲む計測対象領域を自動設
定し、前記大腿骨頸部の領域を抽出する処理を含むこと
を特徴とする骨密度計測方法。 - 【請求項２】請求項１において、前記大腿骨頸部の領域
抽出処理は、大腿骨のシャフトの主軸を算出し、前記主
軸の傾きより予め定めた角αをなす傾きを基準に、前記
大腿骨頸部の主軸を算出する処理を含む骨密度計測方
法。 - 【請求項３】請求項２において、前記角αは、前記シャ
フトの主軸に対して反時計回りに１２５°〜１３０°で
ある骨密度計測方法。 - 【請求項４】高低２種類のエネルギで撮影した大腿骨中
枢端部を含むＸ線画像より差分画像を作成し、前記差分
画像に基づいて大転子部の骨密度を計測する方法におい
て、前記大腿骨中枢端部を囲む計測対象領域を自動設定
し、前記大転子部の領域を抽出する処理を含むことを特
徴とする骨密度計測方法。 - 【請求項５】請求項４において、前記大転子部の領域抽
出処理は、前記大腿骨頸部の主軸を算出し、前記主軸の
傾きより予め定めた角βをなす傾きを基準に、前記大転
子部の領域とシャフトの部分とを分離する境界線を算出
する処理を含む骨密度計測方法。 - 【請求項６】請求項５において、前記角βは、前記大腿
骨頸部の主軸に対して時計回りに５°〜２０°である骨
密度計測方法。 - 【請求項７】請求項４において、前記大転子部の領域抽
出処理は、前記大腿骨のシャフトの主軸を算出し、前記
主軸の傾きより予め定めた角β′をなす傾きを基準に、
前記大転子部の領域と前記シャフトの部分とを分離する
境界線を算出する処理を含む骨密度計測方法。 - 【請求項８】請求項６において、前記角β′は、前記シ
ャフトの主軸に対して時計回りに６０°〜７０°である
骨密度計測方法。 - 【請求項９】請求項１において、前記計測対象領域を自
動設定する処理は、前記差分画像の下端部より前記大腿
骨のシャフトの左右端点を抽出し、前記端点に基づいて
前記大腿骨中枢端部の位置，大きさを設定する処理を含
む骨密度計測方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6074609A JPH07284020A (ja) | 1994-04-13 | 1994-04-13 | 骨密度計測方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP6074609A JPH07284020A (ja) | 1994-04-13 | 1994-04-13 | 骨密度計測方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH07284020A true JPH07284020A (ja) | 1995-10-27 |
Family
ID=13552092
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP6074609A Pending JPH07284020A (ja) | 1994-04-13 | 1994-04-13 | 骨密度計測方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH07284020A (ja) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003061949A (ja) * | 2001-08-24 | 2003-03-04 | Ge Medical Systems Global Technology Co Llc | 領域抽出方法、領域抽出装置およびｘ線ｃｔ装置 |
| JP2007333471A (ja) * | 2006-06-13 | 2007-12-27 | Shimadzu Corp | 核医学診断装置 |
| US7502641B2 (en) | 2002-07-09 | 2009-03-10 | Aecc Enterprises Limited | Method for imaging the relative motion of skeletal segments |
| JP2013236962A (ja) * | 2007-08-15 | 2013-11-28 | Fujifilm Corp | 画像成分分離装置、方法、およびプログラム |
| JP2016198176A (ja) * | 2015-04-08 | 2016-12-01 | 株式会社日立製作所 | 骨塩量測定装置 |
| WO2016190327A1 (ja) * | 2015-05-28 | 2016-12-01 | 株式会社日立製作所 | 医療用ｘ線測定システム |
| JP2017131427A (ja) * | 2016-01-28 | 2017-08-03 | 株式会社日立製作所 | Ｘ線画像診断装置及び骨密度計測方法 |
| CN110840473A (zh) * | 2019-11-22 | 2020-02-28 | 江研伟 | 一种基于CT薄层扫描Hu值的骨密度测量系统 |
| US10959786B2 (en) | 2015-06-05 | 2021-03-30 | Wenzel Spine, Inc. | Methods for data processing for intra-operative navigation systems |
| JP2023082549A (ja) * | 2021-12-02 | 2023-06-14 | 富士フイルム株式会社 | 放射線画像処理装置、方法およびプログラム |
| WO2024122340A1 (ja) * | 2022-12-09 | 2024-06-13 | キヤノン株式会社 | 画像処理装置、画像処理装置の作動方法、およびプログラム |
-
1994
- 1994-04-13 JP JP6074609A patent/JPH07284020A/ja active Pending
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003061949A (ja) * | 2001-08-24 | 2003-03-04 | Ge Medical Systems Global Technology Co Llc | 領域抽出方法、領域抽出装置およびｘ線ｃｔ装置 |
| US7502641B2 (en) | 2002-07-09 | 2009-03-10 | Aecc Enterprises Limited | Method for imaging the relative motion of skeletal segments |
| JP2007333471A (ja) * | 2006-06-13 | 2007-12-27 | Shimadzu Corp | 核医学診断装置 |
| JP2013236962A (ja) * | 2007-08-15 | 2013-11-28 | Fujifilm Corp | 画像成分分離装置、方法、およびプログラム |
| JP2016198176A (ja) * | 2015-04-08 | 2016-12-01 | 株式会社日立製作所 | 骨塩量測定装置 |
| JP2016220850A (ja) * | 2015-05-28 | 2016-12-28 | 株式会社日立製作所 | 医療用ｘ線測定システム |
| WO2016190327A1 (ja) * | 2015-05-28 | 2016-12-01 | 株式会社日立製作所 | 医療用ｘ線測定システム |
| US10959786B2 (en) | 2015-06-05 | 2021-03-30 | Wenzel Spine, Inc. | Methods for data processing for intra-operative navigation systems |
| JP2017131427A (ja) * | 2016-01-28 | 2017-08-03 | 株式会社日立製作所 | Ｘ線画像診断装置及び骨密度計測方法 |
| CN110840473A (zh) * | 2019-11-22 | 2020-02-28 | 江研伟 | 一种基于CT薄层扫描Hu值的骨密度测量系统 |
| CN110840473B (zh) * | 2019-11-22 | 2023-05-26 | 江研伟 | 一种基于CT薄层扫描Hu值的骨密度测量系统 |
| JP2023082549A (ja) * | 2021-12-02 | 2023-06-14 | 富士フイルム株式会社 | 放射線画像処理装置、方法およびプログラム |
| US12533098B2 (en) | 2021-12-02 | 2026-01-27 | Fujifilm Corporation | Radiation image processing device, radiation image processing method, and radiation image processing program |
| WO2024122340A1 (ja) * | 2022-12-09 | 2024-06-13 | キヤノン株式会社 | 画像処理装置、画像処理装置の作動方法、およびプログラム |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Duryea et al. | Trainable rule‐based algorithm for the measurement of joint space width in digital radiographic images of the knee | |
| US7609871B2 (en) | Automatic region of interest locator for AP spinal images and for hip images in bone densitometry | |
| Stokes et al. | Measurement of axial rotation of vertebrae in scoliosis | |
| Sharp et al. | Computer‐based methods for measuring joint space and estimating erosion volume in the finger and wrist joints of patients with rheumatoid arthritis | |
| JP3134009B2 (ja) | 画像処理方法及び装置 | |
| US11957418B2 (en) | Systems and methods for pre-operative visualization of a joint | |
| US20090214099A1 (en) | Method of suppressing obscuring features in an image | |
| Safari et al. | A semi-automatic algorithm for estimating cobb angle | |
| EP1400203A2 (en) | Computer assisted bone densitometer | |
| JPH06261894A (ja) | 骨塩定量方法 | |
| US7280683B2 (en) | Method, code, and system for assaying joint deformity | |
| JPH07284020A (ja) | 骨密度計測方法 | |
| US20050119593A1 (en) | Method of viewing and controlling balance of the vertebral column | |
| WO2013174402A1 (en) | Locking screw length measurement | |
| JP2007130462A (ja) | 左右対称の器官系の三次元画像の評価方法、医用画像撮影システムおよびコンピュータプログラム | |
| US6560474B2 (en) | Method for evaluating structural strength of cancellous bone using X-ray image | |
| JP4717585B2 (ja) | 医用画像判定装置、医用画像判定方法およびそのプログラム | |
| US7046834B2 (en) | Method for measuring bone mineral density by using X-ray image | |
| JP7554008B2 (ja) | ２ｄのｘ線画像で見えない解剖学的構造の人工知能ベースの検出 | |
| JP2011056024A (ja) | 放射線撮影装置、放射線撮影方法、及びプログラム | |
| Kluck et al. | A 3D parameter can guide concave rod contour for the correction of hypokyphosis in adolescent idiopathic scoliosis | |
| Laine et al. | Image processing for femoral endosteal anatomy detection: description and testing of a computed tomography based program | |
| EP0755219A1 (en) | Method for periprosthetic bone mineral density measurement | |
| CN116309503A (zh) | 基于语义分割的医学影像图像质量评价系统 | |
| JP2024546917A (ja) | 自律的な自己較正手術ロボットのためのシステム及び方法 |